GeoGridParameters.h

#pragma once
#include <isce3/core/forward.h>
#include <isce3/io/forward.h>
 
#include <limits>
 
#include <pyre/journal.h>
 
#include <isce3/core/Common.h>
#include <isce3/core/Constants.h>
#include <isce3/geometry/detail/Rdr2Geo.h>
 
#include "RadarGridParameters.h"
 
namespace isce3 { namespace product {
 
std::string to_string(const GeoGridParameters& geogrid);
 
std::ostream& operator<<(std::ostream& out, const GeoGridParameters& geogrid);
 
class GeoGridParameters {
public:
    GeoGridParameters() = default;

    GeoGridParameters(double geoGridStartX, double geoGridStartY,
                             double geoGridSpacingX, double geoGridSpacingY,
                             int width, int length, int epsgcode);

    void print() const;

    void startX(double x0) { _startX = x0; }

    void startY(double y0) { _startY = y0; }

    void spacingX(double dx) { _spacingX = dx; }

    void spacingY(double dy) { _spacingY = dy; }

    void length(int l) { _length = l; };

    void width(int w) { _width = w; };

    void epsg(int e) { _epsg = e; };

    CUDA_HOSTDEV constexpr double startX() const { return _startX; };

    CUDA_HOSTDEV constexpr double startY() const { return _startY; };

    CUDA_HOSTDEV constexpr double endX() const {
        return _startX + _spacingX * _width;
    };

    CUDA_HOSTDEV constexpr double endY() const {
        return _startY + _spacingY * _length;
    };

    CUDA_HOSTDEV constexpr double spacingX() const { return _spacingX; };

    CUDA_HOSTDEV constexpr double spacingY() const { return _spacingY; };

    CUDA_HOSTDEV constexpr int width() const { return _width; };

    CUDA_HOSTDEV constexpr int length() const { return _length; };

    CUDA_HOSTDEV constexpr int epsg() const { return _epsg; };
 
 
protected:
    double _startX = 0.0;

    double _startY = 0.0;

    double _spacingX = 0.0;

    double _spacingY = 0.0;

    int _width = 0;

    int _length = 0;

    int _epsg = 4326;
};

GeoGridParameters
bbox2GeoGrid(const isce3::product::RadarGridParameters& radar_grid,
             const isce3::core::Orbit& orbit,
             const isce3::core::LUT2d<double>& doppler, double spacing_x,
             double spacing_y, int epsg,
             double min_height = isce3::core::GLOBAL_MIN_HEIGHT,
             double max_height = isce3::core::GLOBAL_MAX_HEIGHT,
             const double margin = 0.0, const int points_per_edge = 11,
             const double threshold = isce3::geometry::detail::DEFAULT_TOL_HEIGHT,
             const double height_threshold = 100);

GeoGridParameters
bbox2GeoGridScaled(const isce3::product::RadarGridParameters& radar_grid,
                   const isce3::core::Orbit& orbit,
                   const isce3::core::LUT2d<double>& doppler,
                   const isce3::io::Raster& dem_raster,
                   double spacing_scale = 1.0,
                   double min_height = isce3::core::GLOBAL_MIN_HEIGHT,
                   double max_height = isce3::core::GLOBAL_MAX_HEIGHT,
                   const double margin = 0.0, const int points_per_edge = 11,
                   const double threshold = isce3::geometry::detail::DEFAULT_TOL_HEIGHT,
                   const double height_threshold = 100);
 
 
}} // namespace isce3::product