DRAFT: Add support for PFA RTC

cxx/isce3/Headers.cmake

@@ -158,7 +158,9 @@ product/Metadata.h
 product/ProcessingInformation.h
 product/RadarGridProduct.h
 product/GeoGridProduct.h
+product/RngAzmGridParameters.h
 product/RadarGridParameters.h
+product/PolarGridParameters.h
 product/Serialization.h
 product/SubSwaths.h
 product/Swath.h

cxx/isce3/Sources.cmake

@@ -91,6 +91,8 @@ math/polyfunc.cpp
 math/RootFind1dNewton.cpp
 math/RootFind1dSecant.cpp
 polsar/symmetrize.cpp
+product/RngAzmGridParameters.cpp
+product/PolarGridParameters.cpp
 product/RadarGridParameters.cpp
 product/GeoGridParameters.cpp
 product/RadarGridProduct.cpp

cxx/isce3/geocode/GeocodeCov.h

@@ -20,6 +20,7 @@
 // isce3::product
 #include <isce3/product/RadarGridProduct.h>
 #include <isce3/product/RadarGridParameters.h>
+#include <isce3/product/PolarGridParameters.h>
 
 // isce3::geometry
 #include <isce3/geometry/RTC.h>
@@ -32,8 +33,17 @@ enum geocodeOutputMode {
     AREA_PROJECTION = 1,
 };
 
+/** Wrapper function to simplify calling geo2rdr with different grid types */
+static int _geo2rdrGrid(const Vec3&, const Ellipsoid&, const Orbit&, const LUT2d<double>&,
+                        double&, double&, const isce3::product::RadarGridParameters&,
+                        double, int, double, bool);
 
-template<class T>
+/** Wrapper function to simplify calling geo2rdr with different grid types */
+static int _geo2rdrGrid(const Vec3&, const Ellipsoid&, const Orbit&, const LUT2d<double>&,
+                        double&, double&, const isce3::product::PolarGridParameters&,
+                        double, int, double, bool);
+
+template<class T, class T_grid = isce3::product::RadarGridParameters>
 class Geocode {
 public:
     /** Geocode data from slant-range to map coordinates
@@ -145,7 +155,7 @@ class Geocode {
      * @param[in]  max_block_size      Maximum block size (per thread)
      * @param[in]  dem_interp_method   DEM interpolation method
      */
-    void geocode(const isce3::product::RadarGridParameters& radar_grid,
+    void geocode(const T_grid& radar_grid,
             isce3::io::Raster& input_raster, isce3::io::Raster& output_raster,
             isce3::io::Raster& dem_raster,
             geocodeOutputMode output_mode = geocodeOutputMode::INTERP,
@@ -286,7 +296,7 @@ class Geocode {
      * @param[in]  dem_interp_method   DEM interpolation method
      */
     template<class T_out>
-    void geocodeInterp(const isce3::product::RadarGridParameters& radar_grid,
+    void geocodeInterp(const T_grid& radar_grid,
             isce3::io::Raster& input_raster, isce3::io::Raster& output_raster,
             isce3::io::Raster& dem_raster, bool flag_apply_rtc = false,
             bool flag_az_baseband_doppler = false, bool flatten = false,
@@ -429,7 +439,7 @@ class Geocode {
      */
     template<class T_out>
     void geocodeAreaProj(
-            const isce3::product::RadarGridParameters& radar_grid,
+            const T_grid& radar_grid,
             isce3::io::Raster& input_raster, isce3::io::Raster& output_raster,
             isce3::io::Raster& dem_raster,
             double geogrid_upsampling = 1,
@@ -494,8 +504,7 @@ class Geocode {
      * @param[in]  radar_grid          Radar grid
      * @param[in]  dem_raster          Input DEM raster
      */
-    void updateGeoGrid(const isce3::product::RadarGridParameters& radar_grid,
-                       isce3::io::Raster& dem_raster);
+    void updateGeoGrid(const T_grid& radar_grid, isce3::io::Raster& dem_raster);
 
     // Get/set data interpolator
     isce3::core::dataInterpMethod dataInterpolator() const 
@@ -555,7 +564,7 @@ class Geocode {
     the Geocode object geogrid attributes.
     */
     void _getRadarGridBoundaries(
-            const isce3::product::RadarGridParameters& radar_grid,
+            const T_grid& radar_grid,
             isce3::io::Raster& input_raster, isce3::io::Raster& dem_raster,
             isce3::core::ProjectionBase* proj, double geogrid_upsampling,
             bool flag_upsample_radar_grid,
@@ -572,7 +581,7 @@ class Geocode {
             const int k_end, double geogrid_upsampling, double* a11,
             double* r11, double* y_min, double* x_min, double* y_max,
             double* x_max,
-            const isce3::product::RadarGridParameters& radar_grid,
+            const T_grid& radar_grid,
             isce3::core::ProjectionBase* proj,
             isce3::geometry::DEMInterpolator& dem_interp_block,
             const std::function<Vec3(double, double,
@@ -591,7 +600,7 @@ class Geocode {
     */
     bool _checkLoadEntireRslcCorners(const double y0, const double x0,
             const double yf, const double xf,
-            const isce3::product::RadarGridParameters& radar_grid,
+            const T_grid& radar_grid,
             isce3::core::ProjectionBase* proj,
             const std::function<Vec3(double, double,
                     const isce3::geometry::DEMInterpolator&,
@@ -605,7 +614,7 @@ class Geocode {
     void _getRadarPositionBorder(double geogrid_upsampling, const double dem_y1,
             const double dem_x1, const double dem_yf, const double dem_xf,
             double* a_min, double* r_min, double* a_max, double* r_max,
-            const isce3::product::RadarGridParameters& radar_grid,
+            const T_grid& radar_grid,
             isce3::core::ProjectionBase* proj,
             const std::function<Vec3(double, double,
                     const isce3::geometry::DEMInterpolator&,
@@ -615,7 +624,7 @@ class Geocode {
             const isce3::core::LUT2d<double>& slant_range_correction = {});
 
     template<class T2, class T_out>
-    void _runBlock(const isce3::product::RadarGridParameters& radar_grid,
+    void _runBlock(const T_grid& radar_grid,
             bool is_radar_grid_single_block,
             std::vector<std::unique_ptr<isce3::core::Matrix<T2>>>& rdrData,
             int block_size_y, int block_size_with_upsampling_y, int block_y,
@@ -652,7 +661,7 @@ class Geocode {
 
     std::string _get_nbytes_str(long nbytes);
 
-    int _geo2rdr(const isce3::product::RadarGridParameters& radar_grid,
+    int _geo2rdr(const T_grid& radar_grid,
             double x, double y, double& azimuthTime, double& slantRange,
             isce3::geometry::DEMInterpolator& demInterp,
             isce3::core::ProjectionBase* proj, float& dem_value);
@@ -720,7 +729,7 @@ class Geocode {
             const int radarBlockLength, const int azimuthFirstLine,
             const int rangeFirstPixel,
             const isce3::core::Interpolator<T_out>* interp,
-            const isce3::product::RadarGridParameters& radarGrid,
+            const T_grid& radarGrid,
             const bool flag_az_baseband_doppler, const bool flatten,
             isce3::io::Raster* phase_screen_raster,
             isce3::core::Matrix<float>& phase_screen_array,