Geo2rdr.h

//-*- C++ -*-
//-*- coding: utf-8 -*-
//
// Author: Bryan Riel, Joshua Cohen
// Copyright 2017-2018
 
#pragma once
 
// pyre
#include <pyre/journal.h>
 
// isce3::core
#include <isce3/core/Metadata.h>
#include <isce3/core/Orbit.h>
#include <isce3/core/Poly2d.h>
#include <isce3/core/LUT1d.h>
#include <isce3/core/Ellipsoid.h>
#include <isce3/core/Peg.h>
#include <isce3/core/Projections.h>
 
// isce3::io
#include <isce3/io/Raster.h>
 
// isce3::product
#include <isce3/product/RadarGridProduct.h>
#include <isce3/product/RadarGridParameters.h>
 
#include <limits>
 
// Declaration
namespace isce3 {
    namespace geometry {
        class Geo2rdr;
    }
}

class isce3::geometry::Geo2rdr {
public:

    Geo2rdr(const isce3::product::RadarGridProduct &,
            char frequency = 'A',
            bool nativeDoppler = false);

    Geo2rdr(const isce3::core::Ellipsoid &,
            const isce3::core::Orbit &,
            const isce3::core::LUT2d<double> &,
            const isce3::core::Metadata &);

    Geo2rdr(const isce3::product::RadarGridParameters & radarGrid,
            const isce3::core::Orbit & orbit,
            const isce3::core::Ellipsoid & ellipsoid,
            const isce3::core::LUT2d<double> & doppler = {});

    void threshold(double t) { _threshold = t; }

    void numiter(int n) { _numiter = n; }

    void linesPerBlock(size_t linesPerBlock) { _linesPerBlock = linesPerBlock; }

    void geo2rdr(isce3::io::Raster & topoRaster,
                 isce3::io::Raster & rgoffRaster,
                 isce3::io::Raster & azoffRaster,
                 double azshift=0.0, double rgshift=0.0);

    void geo2rdr(isce3::io::Raster & topoRaster,
                 const std::string & outdir,
                 double azshift=0.0, double rgshift=0.0);

    const double NULL_VALUE =  -1000000.0;
 
    // Getters for isce objects

    const isce3::core::Orbit & orbit() const { return _orbit; }

    const isce3::core::Ellipsoid & ellipsoid() const { return _ellipsoid; }

    const isce3::core::LUT2d<double> & doppler() const { return _doppler; }

    const isce3::product::RadarGridParameters & radarGridParameters() const { return _radarGrid; }
 
    // Get geo2rdr processing options

    double threshold() const { return _threshold; }

    int numiter() const { return _numiter; }

    size_t linesPerBlock() const { return _linesPerBlock; }
 
private:

    void _printExtents(pyre::journal::info_t &,
                       double, double, double,
                       double, double, double,
                       size_t, size_t);

    void _checkOrbitInterpolation(double);
 
    // isce3::core objects
    isce3::core::Ellipsoid _ellipsoid;
    isce3::core::Orbit _orbit;
    isce3::core::LUT2d<double> _doppler;
 
    // RadarGridParameters
    isce3::product::RadarGridParameters _radarGrid;
 
    // Projection related data
    isce3::core::ProjectionBase * _projTopo;
 
    // Processing parameters
    int _numiter;
    double _threshold = 1e-8;
    size_t _linesPerBlock = 1000;
};
 
// Get inline implementations for Geo2rdr
#define ISCE_GEOMETRY_GEO2RDR_ICC
#include "Geo2rdr.icc"
#undef ISCE_GEOMETRY_GEO2RDR_ICC