WaterLevelDecon.h

#ifndef __SIMPLE_WATER_LEVEL_DECON_H__
#define __SIMPLE_WATER_LEVEL_DECON_H__
#include <vector>
#include "mspass/utility/Metadata.h"
#include "mspass/algorithms/deconvolution/ScalarDecon.h"
#include "mspass/algorithms/deconvolution/FFTDeconOperator.h"
#include "mspass/algorithms/deconvolution/ShapingWavelet.h"
#include "mspass/seismic/CoreTimeSeries.h"
namespace mspass::algorithms::deconvolution{
class WaterLevelDecon : public FFTDeconOperator, public ScalarDecon
{
public:
    WaterLevelDecon(const WaterLevelDecon &parent);
    WaterLevelDecon(const mspass::utility::Metadata &md);
    WaterLevelDecon(const mspass::utility::Metadata &md,const std::vector<double> &wavelet,const std::vector<double> &data);
    void changeparameter(const mspass::utility::Metadata &md);
    void process();
    mspass::seismic::CoreTimeSeries actual_output();
    mspass::seismic::CoreTimeSeries inverse_wavelet(const double t0parent=0.0);
    mspass::seismic::CoreTimeSeries inverse_wavelet() ;
    mspass::utility::Metadata QCMetrics();
private:
    int read_metadata(const mspass::utility::Metadata &md);
    int apply();
    double wlv;
    /* QC metrics.   */
    /* This is the fraction of frequencies below the water level */
    double regularization_fraction;
};
}
#endif