Stats.h

#include <isce3/core/TypeTraits.h>
#include <isce3/core/Constants.h>
#include <isce3/core/blockProcessing.h>
#include <isce3/io/Raster.h>
 
namespace isce3 { namespace math {

template<class T>
struct Stats {
    using type = T;
    using T_real = typename isce3::real<T>::type;
    T min = std::numeric_limits<T_real>::quiet_NaN();
    T max = std::numeric_limits<T_real>::quiet_NaN();
    T mean = 0;
 
    long long n_valid = 0;
 
    T_real sample_stddev() const;

    void update(const Stats<T>& other);

    void update(const T& value);

    void update(const T* values, size_t size, size_t stride = 1);

    Stats(const T* values, size_t size, size_t stride = 1);
 
    Stats() = default;
 
private:
    double real_valued_mean = 0;
    double square_diff_sum = 0;
};

template<class T>
struct StatsRealImag {
    using type = std::complex<T>;
 
    Stats<T> real;
    Stats<T> imag;
    long long n_valid = 0;

    void update(const StatsRealImag<T>& other);

    void update(const std::complex<T>& value);

    void update(const std::complex<T>* values, size_t size, size_t stride = 1);

    StatsRealImag(const std::complex<T>* values, size_t size, size_t stride = 1);
 
    StatsRealImag() = default;
};

template<class T>
std::vector<isce3::math::Stats<T>> computeRasterStats(
    isce3::io::Raster& input_raster,
    isce3::core::MemoryModeBlocksY memory_mode = 
        isce3::core::MemoryModeBlocksY::AutoBlocksY);

template<class T>
std::vector<isce3::math::StatsRealImag<T>> computeRasterStatsRealImag(
    isce3::io::Raster& input_raster, 
    isce3::core::MemoryModeBlocksY memory_mode = 
        isce3::core::MemoryModeBlocksY::AutoBlocksY);
 
}}