loadDem.h

#pragma once
#include <isce3/geometry/DEMInterpolator.h>
#include <isce3/io/Raster.h>
#include <isce3/product/GeoGridParameters.h>
#include <isce3/core/forward.h>
 
namespace isce3 { namespace geometry {
 

DEMInterpolator DEMRasterToInterpolator(
        isce3::io::Raster& demRaster,
        const isce3::product::GeoGridParameters& geoGrid,
        const int demMarginInPixels = 50,
        const isce3::core::dataInterpMethod demInterpMethod =
            isce3::core::BIQUINTIC_METHOD);
 

isce3::geometry::DEMInterpolator DEMRasterToInterpolator(
        isce3::io::Raster& demRaster,
        const isce3::product::GeoGridParameters& geoGrid, const int lineStart,
        const int blockLength, const int blockWidth,
        const int demMarginInPixels = 50,
        const isce3::core::dataInterpMethod demInterpMethod =
            isce3::core::BIQUINTIC_METHOD);
 

isce3::error::ErrorCode loadDemFromProj(
    isce3::io::Raster& dem_raster,
    const double minX, const double maxX, const double minY,
    const double maxY, isce3::geometry::DEMInterpolator* dem_interp,
    isce3::core::ProjectionBase* proj = nullptr,
    const int dem_margin_x_in_pixels = 100,
    const int dem_margin_y_in_pixels = 100,
    const int dem_raster_band = 1);
 
/*
 Interpolate DEM at position (x, y) considering that input_proj and
 dem_interp have same coordinate systems. The function is written to
 have the same interface of getDemCoordsDiffEpsg()
 * @param[in]  x           X-coordinate in input coordinates
 * @param[in]  y           Y-coordinate in input coordinates
 * @param[in]  dem_interp  DEM interpolation object
 * @param[in]  input_proj  Input projection object
 * @returns                3-elements vector containing the x and
 * y coordinates over DEM projection coordinates, and interpolated
 * DEM value at that position: {x_dem, y_dem, z_dem}
 */
isce3::core::Vec3 getDemCoordsSameEpsg(double x, double y,
        const DEMInterpolator& dem_interp,
        isce3::core::ProjectionBase* input_proj);
 
/*
 Convert x and y coordinates to from input_proj coordinates to
 DEM (dem_interp) coordinates and interpolate DEM at that position.
 3-elements vector containing the
 * @param[in]  x           X-coordinate in input coordinates
 * @param[in]  y           Y-coordinate in input coordinates
 * @param[in]  dem_interp  DEM interpolation object
 * @param[in]  input_proj  Input projection object
 * @returns                3-elements vector containing the x and
 * y coordinates over DEM projection coordinates, and interpolated
 * DEM value at that position: {x_dem, y_dem, z_dem}
 */
isce3::core::Vec3 getDemCoordsDiffEpsg(double x, double y,
        const DEMInterpolator& dem_interp,
        isce3::core::ProjectionBase* input_proj);
 
 
}} // namespace isce3::geometry