Rdr2Geo.h

#pragma once
 
#include <isce3/core/forward.h>
 
#include <cmath>
 
#include <isce3/core/Common.h>
#include <isce3/core/LookSide.h>
#include <isce3/error/ErrorCode.h>
 
namespace isce3 { namespace geometry { namespace detail {

struct Rdr2GeoParams {
    double threshold = 1e-8;

    int maxiter = 25;

    int extraiter = 15;
};

template<class Orbit, class DEMInterpolator>
CUDA_HOSTDEV isce3::error::ErrorCode
rdr2geo(isce3::core::Vec3* llh, double t, double r, double fd,
        const Orbit& orbit, const DEMInterpolator& dem,
        const isce3::core::Ellipsoid& ellipsoid, double wvl,
        isce3::core::LookSide side, double h0 = 0.,
        const Rdr2GeoParams& params = {});

template<class DEMInterpolator>
CUDA_HOSTDEV isce3::error::ErrorCode
rdr2geo(isce3::core::Vec3* llh, const isce3::core::Pixel& pixel,
        const isce3::core::Basis& tcnbasis, const isce3::core::Vec3& pos,
        const isce3::core::Vec3& vel, const DEMInterpolator& dem,
        const isce3::core::Ellipsoid& ellipsoid, isce3::core::LookSide side,
        double h0 = 0., const Rdr2GeoParams& params = {});
 

inline constexpr double DEFAULT_TOL_HEIGHT = 1e-5;

struct Rdr2GeoBracketParams {
    double tol_height = DEFAULT_TOL_HEIGHT;

    double look_min = 0.0;

    double look_max = M_PI / 2;
};

template<class Orbit, class DEMInterpolator>
CUDA_HOSTDEV isce3::error::ErrorCode
rdr2geo_bracket(isce3::core::Vec3* xyz,
        double aztime, double slantRange, double doppler, const Orbit& orbit,
        const DEMInterpolator& dem, const isce3::core::Ellipsoid& ellipsoid,
        double wavelength, isce3::core::LookSide side,
        const Rdr2GeoBracketParams& params = {});
 
}}} // namespace isce3::geometry::detail
 
#include "Rdr2Geo.icc"