geometry.h

Go to the documentation of this file.

// -*- C++ -*-
// -*- coding: utf-8 -*-
//
// Author: Bryan Riel , Hirad Ghaemi
// Copyright 2017-2018 , 2020-2021
//
 
#pragma once
 
#include "forward.h"
#include <isce3/core/forward.h>
#include <isce3/product/forward.h>
 
#include <optional>
#include <tuple>
 
#include <Eigen/Dense>
 
#include <isce3/core/Constants.h>
#include <isce3/core/Ellipsoid.h>
#include <isce3/geometry/DEMInterpolator.h>
 
// Declaration
namespace isce3 {
namespace geometry {

int rdr2geo(double aztime, double slantRange, double doppler,
        const isce3::core::Orbit& orbit,
        const isce3::core::Ellipsoid& ellipsoid,
        const DEMInterpolator& demInterp, isce3::core::Vec3& targetLLH,
        double wvl, isce3::core::LookSide side, double threshold, int maxIter,
        int extraIter);

int rdr2geo(const isce3::core::Pixel& pixel, const isce3::core::Basis& TCNbasis,
        const isce3::core::Vec3& pos, const isce3::core::Vec3& vel,
        const isce3::core::Ellipsoid& ellipsoid,
        const DEMInterpolator& demInterp, isce3::core::Vec3& targetLLH,
        isce3::core::LookSide side, double threshold, int maxIter,
        int extraIter);

int rdr2geo(const isce3::core::Vec3& radarXYZ, const isce3::core::Vec3& axis,
        double angle, double range, const DEMInterpolator& dem,
        isce3::core::Vec3& targetXYZ, isce3::core::LookSide side,
        double threshold, int maxIter, int extraIter);

int geo2rdr(const isce3::core::Vec3& inputLLH,
        const isce3::core::Ellipsoid& ellipsoid,
        const isce3::core::Orbit& orbit, const isce3::core::Poly2d& doppler,
        double& aztime, double& slantRange, double wavelength,
        double startingRange, double rangePixelSpacing, size_t rwidth,
        isce3::core::LookSide side, double threshold, int maxIter,
        double deltaRange);

int geo2rdr(const isce3::core::Vec3& inputLLH,
        const isce3::core::Ellipsoid& ellipsoid,
        const isce3::core::Orbit& orbit,
        const isce3::core::LUT2d<double>& doppler, double& aztime,
        double& slantRange, double wavelength, isce3::core::LookSide side,
        double threshold, int maxIter, double deltaRange);

void computeDEMBounds(const isce3::core::Orbit& orbit,
        const isce3::core::Ellipsoid& ellipsoid,
        const isce3::core::LUT2d<double>& doppler,
        const isce3::product::RadarGridParameters& radarGrid, size_t xoff,
        size_t yoff, size_t xsize, size_t ysize, double margin, double& min_lon,
        double& min_lat, double& max_lon, double& max_lat);
 
template<class T>
double _compute_doppler_aztime_diff(isce3::core::Vec3 dr,
        isce3::core::Vec3 satvel, T& doppler, double wavelength, double aztime,
        double slantRange, double deltaRange);

isce3::core::Vec3 nedVector(
        double lon, double lat, const isce3::core::Vec3& vel);

isce3::core::Vec3 nwuVector(
        double lon, double lat, const isce3::core::Vec3& vel);

isce3::core::Vec3 enuVector(
        double lon, double lat, const isce3::core::Vec3& vel);

double heading(double lon, double lat, const isce3::core::Vec3& vel);

double slantRangeFromLookVec(const isce3::core::Vec3& pos,
        const isce3::core::Vec3& lkvec,
        const isce3::core::Ellipsoid& ellips = {});

std::pair<int, double> srPosFromLookVecDem(double& sr,
        isce3::core::Vec3& tg_pos, isce3::core::Vec3& llh,
        const isce3::core::Vec3& sc_pos, const isce3::core::Vec3& lkvec,
        const DEMInterpolator& dem_interp = {}, double hgt_err = 0.5,
        int num_iter = 10, const isce3::core::Ellipsoid& ellips = {},
        std::optional<double> initial_height = {});

std::tuple<double, double> lookIncAngFromSlantRange(double slant_range,
        const isce3::core::Orbit& orbit, std::optional<double> az_time = {},
        const DEMInterpolator& dem_interp = {},
        const isce3::core::Ellipsoid& ellips = {});

std::tuple<Eigen::ArrayXd, Eigen::ArrayXd> lookIncAngFromSlantRange(
        const Eigen::Ref<const Eigen::ArrayXd>& slant_range,
        const isce3::core::Orbit& orbit, std::optional<double> az_time = {},
        const DEMInterpolator& dem_interp = {},
        const isce3::core::Ellipsoid& ellips = {});

double compute_mean_dem(const DEMInterpolator& dem);
 
} // namespace geometry
} // namespace isce3