geocodeSlc.h

#pragma once
#include <cstddef>
#include <vector>
#include <optional>
#include <isce3/core/EMatrix.h>
#include <isce3/core/forward.h>
#include <isce3/core/LUT2d.h>
#include <isce3/core/Poly2d.h>
#include <isce3/io/forward.h>
#include <isce3/product/forward.h>
 
namespace isce3 { namespace geocode {
 
typedef Eigen::Ref<isce3::core::EArray2D<std::complex<float>>> EArray2dc64;
typedef Eigen::Ref<isce3::core::EArray2D<double>> EArray2df64;
typedef Eigen::Ref<isce3::core::EArray2D<unsigned char>> EArray2duc8;

template<typename AzRgFunc = isce3::core::Poly2d>
void geocodeSlc(isce3::io::Raster& outputRaster, isce3::io::Raster& inputRaster,
        isce3::io::Raster& demRaster,
        const isce3::product::RadarGridParameters& radarGrid,
        const isce3::product::GeoGridParameters& geoGrid,
        const isce3::core::Orbit& orbit,
        const isce3::core::LUT2d<double>& nativeDoppler,
        const isce3::core::LUT2d<double>& imageGridDoppler,
        const isce3::core::Ellipsoid& ellipsoid,
        const double& thresholdGeo2rdr, const int& numiterGeo2rdr,
        const size_t& linesPerBlock,
        const bool flatten = true,
        const bool reramp = true,
        const AzRgFunc& azCarrier = AzRgFunc(),
        const AzRgFunc& rgCarrier = AzRgFunc(),
        const isce3::core::LUT2d<double>& azTimeCorrection = {},
        const isce3::core::LUT2d<double>& sRangeCorrection = {},
        const bool flattenWithCorrectedSRng = false,
        const std::complex<float> invalidValue =
            std::complex<float>(std::numeric_limits<float>::quiet_NaN(),
                                std::numeric_limits<float>::quiet_NaN()),
        isce3::io::Raster* carrierPhaseRaster = nullptr,
        isce3::io::Raster* flattenPhaseRaster = nullptr);

template<typename AzRgFunc = isce3::core::Poly2d>
void geocodeSlc(isce3::io::Raster& outputRaster, isce3::io::Raster& inputRaster,
        isce3::io::Raster& demRaster,
        const isce3::product::RadarGridParameters& radarGrid,
        const isce3::product::RadarGridParameters& slicedRadarGrid,
        const isce3::product::GeoGridParameters& geoGrid,
        const isce3::core::Orbit& orbit,
        const isce3::core::LUT2d<double>& nativeDoppler,
        const isce3::core::LUT2d<double>& imageGridDoppler,
        const isce3::core::Ellipsoid& ellipsoid,
        const double& thresholdGeo2rdr, const int& numiterGeo2rdr,
        const size_t& linesPerBlock,
        const bool flatten = true,
        const bool reramp = true,
        const AzRgFunc& azCarrier = AzRgFunc(),
        const AzRgFunc& rgCarrier = AzRgFunc(),
        const isce3::core::LUT2d<double>& azTimeCorrection = {},
        const isce3::core::LUT2d<double>& sRangeCorrection = {},
        const bool flattenWithCorrectedSRng = false,
        const std::complex<float> invalidValue =
            std::complex<float>(std::numeric_limits<float>::quiet_NaN(),
                                std::numeric_limits<float>::quiet_NaN()),
        isce3::io::Raster* carrierPhaseRaster = nullptr,
        isce3::io::Raster* flattenPhaseRaster = nullptr);
 

template<typename AzRgFunc = isce3::core::Poly2d>
void geocodeSlc(
        std::vector<EArray2dc64>& geoDataBlocks,
        EArray2df64 carrierPhaseBlock,
        EArray2df64 flattenPhaseBlock,
        const std::vector<EArray2dc64>& rdrDataBlocks,
        isce3::io::Raster& demRaster,
        const isce3::product::RadarGridParameters& radarGrid,
        const isce3::product::RadarGridParameters& slicedRadarGrid,
        const isce3::product::GeoGridParameters& geoGrid,
        const isce3::core::Orbit& orbit,
        const isce3::core::LUT2d<double>& nativeDoppler,
        const isce3::core::LUT2d<double>& imageGridDoppler,
        const isce3::core::Ellipsoid& ellipsoid,
        const double& thresholdGeo2rdr, const int& numiterGeo2rdr,
        std::optional<isce3::geocode::EArray2duc8> maskBlock = std::nullopt,
        const size_t& azimuthFirstLine = 0,
        const size_t& rangeFirstPixel = 0,
        const bool flatten = true,
        const bool reramp = true,
        const AzRgFunc& azCarrier = AzRgFunc(),
        const AzRgFunc& rgCarrier = AzRgFunc(),
        const isce3::core::LUT2d<double>& azTimeCorrection = {},
        const isce3::core::LUT2d<double>& sRangeCorrection = {},
        const bool flattenWithCorrectedSRng = false,
        const std::complex<float> invalidValue =
            std::complex<float>(std::numeric_limits<float>::quiet_NaN(),
                                std::numeric_limits<float>::quiet_NaN()),
        const isce3::product::SubSwaths* subswaths = nullptr);
 
}} // namespace isce3::geocode