Crossmul.h

#pragma once
 
#include "forward.h"
 
#include <complex>
#include <isce3/core/LUT1d.h>
#include <isce3/io/forward.h>

class isce3::signal::Crossmul {
    public:
        // Constructor from product
        Crossmul() {};
 
        ~Crossmul() {};

        void crossmul(isce3::io::Raster& refSlcRaster,
                    isce3::io::Raster& secSlcRaster,
                    isce3::io::Raster& ifgRaster,
                    isce3::io::Raster& coherence,
                    isce3::io::Raster* rngOffsetRaster = nullptr) const;

        inline void doppler(isce3::core::LUT1d<double>,
                            isce3::core::LUT1d<double>);

        inline void refDoppler(isce3::core::LUT1d<double> refDopp) { _refDoppler = refDopp; }

        inline const isce3::core::LUT1d<double> & refDoppler() const { return _refDoppler; }

        inline void secDoppler(isce3::core::LUT1d<double> secDopp) { _secDoppler = secDopp; }

        inline const isce3::core::LUT1d<double> & secDoppler() const { return _secDoppler; }

        inline void startingRangeShift(double rng_shift) { _offsetStartingRangeShift = rng_shift; }

        inline double startingRangeShift() const { return _offsetStartingRangeShift; }

        inline void rangePixelSpacing(double rgPxlSpacing) { _rangePixelSpacing = rgPxlSpacing; }

        inline double rangePixelSpacing() const { return _rangePixelSpacing; }

        inline void wavelength(double wvl) { _wavelength = wvl; }

        inline double wavelength() const { return _wavelength; }

        inline void rangeLooks(int);

        inline int rangeLooks() const { return _rangeLooks; }

        inline void azimuthLooks(int);

        inline int azimuthLooks() const { return _azimuthLooks; }

        inline void oversampleFactor(size_t oversamp) { _oversampleFactor = oversamp; }

        inline size_t oversampleFactor() const { return _oversampleFactor; }

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

        inline size_t linesPerBlock() const { return _linesPerBlock; }

        inline bool multiLookEnabled() const { return _multiLookEnabled; }

        inline void rangeFrequencyShift(std::valarray<std::complex<float>> &refAvgSpectrum,
                std::valarray<std::complex<float>> &secAvgSpectrum,
                std::valarray<double> &rangeFrequencies,
                size_t linesPerBlockData,
                size_t fft_size,
                double &frequencyShift);

        inline void getPeakIndex(std::valarray<float> data,
                                size_t &peakIndex);
 
    private:
        //Doppler LUT for the refernce SLC
        isce3::core::LUT1d<double> _refDoppler;
 
        //Doppler LUT for the secondary SLC
        isce3::core::LUT1d<double> _secDoppler;
 
        // starting range shifts between the secondary and reference RSLC in meters
        double _offsetStartingRangeShift = 0.0;
 
        // range pixel spacing
        double _rangePixelSpacing;
 
        // radar wavelength
        double _wavelength;
 
        // number of range looks
        int _rangeLooks = 1;
 
        // number of azimuth looks
        int _azimuthLooks = 1;
 
        bool _multiLookEnabled = false;
 
        // number of lines per block
        size_t _linesPerBlock = 1024;
 
        // upsampling factor
        size_t _oversampleFactor = 1;
 
 
};
 
// Get inline implementations for Crossmul
#define ISCE_SIGNAL_CROSSMUL_ICC
#include "Crossmul.icc"
#undef ISCE_SIGNAL_CROSSMUL_ICC