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/Product.h>
#include <isce3/product/RadarGridParameters.h>

// Declaration
namespace isce3 {
    namespace geometry {
        class Geo2rdr;
    }
}

class isce3::geometry::Geo2rdr {
public:

    Geo2rdr(const isce3::product::Product &,
            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 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 = -1.0e6;

    // 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; }

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;
    size_t _linesPerBlock = 1000;
};

// Get inline implementations for Geo2rdr
#define ISCE_GEOMETRY_GEO2RDR_ICC
#include "Geo2rdr.icc"
#undef ISCE_GEOMETRY_GEO2RDR_ICC