gpuDEMInterpolator.h

//
// Author: Bryan Riel
// Copyright 2018
//

#pragma once

#include "forward.h"
#include <isce3/geometry/forward.h>
#include <isce3/cuda/core/forward.h>

#include <isce3/cuda/core/gpuProjections.h>

// DEMInterpolator declaration
class isce3::cuda::geometry::gpuDEMInterpolator {

    typedef isce3::core::Vec3 Vec3;

    public:
        CUDA_HOSTDEV inline gpuDEMInterpolator() :
            _haveRaster(false), _refHeight(0.0), _interpMethod(isce3::core::BILINEAR_METHOD),
            _xstart(0.0), _ystart(0.0), _deltax(1.0), _deltay(1.0), _owner(false) {}

        CUDA_HOSTDEV inline gpuDEMInterpolator(float height) :
            _haveRaster(false), _refHeight(height), _interpMethod(isce3::core::BILINEAR_METHOD),
            _xstart(0.0), _ystart(0.0), _deltax(1.0), _deltay(1.0), _owner(false) {}

        ~gpuDEMInterpolator();

        CUDA_HOST gpuDEMInterpolator(const isce3::geometry::DEMInterpolator&);

        CUDA_HOSTDEV gpuDEMInterpolator(gpuDEMInterpolator &);

        CUDA_DEV float interpolateLonLat(double lon, double lat) const;

        CUDA_DEV float interpolateXY(double x, double y) const;

        CUDA_DEV inline double midX() const { return _xstart + 0.5*_width*_deltax; }
        CUDA_DEV inline double midY() const { return _ystart + 0.5*_length*_deltay; }

        CUDA_DEV Vec3 midLonLat() const;

        CUDA_HOSTDEV inline double xStart() const { return _xstart; }

        CUDA_HOSTDEV inline double yStart() const { return _ystart; }

        CUDA_HOSTDEV inline double deltaX() const { return _deltax; }

        CUDA_HOSTDEV inline double deltaY() const { return _deltay; }

        CUDA_HOSTDEV inline bool haveRaster() const { return _haveRaster; }

        CUDA_HOSTDEV inline float refHeight() const { return _refHeight; }

        CUDA_HOSTDEV inline size_t length() const { return _length; }

        CUDA_HOSTDEV inline size_t width() const { return _width; }

        CUDA_HOSTDEV inline int epsgCode() const { return _epsgcode; }

        CUDA_HOSTDEV inline isce3::core::dataInterpMethod interpMethod() const {
            return _interpMethod;
        }

        CUDA_HOSTDEV inline isce3::cuda::core::ProjectionBase ** proj() const {
            return _proj;
        }

        CUDA_HOSTDEV inline isce3::cuda::core::gpuInterpolator<float> ** interp() const {
            return _interp;
        }

        CUDA_HOST void initProjInterp();

        CUDA_HOST void finalizeProjInterp();

        // Make DEM pointer data public for now
        float * _dem;

    private:
        // Flag indicating whether we have access to a DEM raster
        bool _haveRaster;
        // Constant value if no raster is provided
        float _refHeight;
        // Pointers to ProjectionBase
        int _epsgcode;
        isce3::cuda::core::ProjectionBase ** _proj;
        // Pointer to an Interpolator
        isce3::core::dataInterpMethod _interpMethod;
        isce3::cuda::core::gpuInterpolator<float> ** _interp;
        //isce3::core::Matrix<float> _dem;
        // DEM dimensions
        size_t _length, _width;
        // Starting x/y for DEM subset and spacing
        double _xstart, _ystart, _deltax, _deltay;
        // Boolean for owning memory
        bool _owner;
};