LeastSquareDecon.h

#ifndef __SIMPLE_LEAST_SQUARE_DECON_H__
#define __SIMPLE_LEAST_SQUARE_DECON_H__
#include "mspass/algorithms/deconvolution/FFTDeconOperator.h"
#include "mspass/algorithms/deconvolution/ScalarDecon.h"
#include "mspass/algorithms/deconvolution/ShapingWavelet.h"
#include "mspass/seismic/CoreTimeSeries.h"
#include "mspass/utility/Metadata.h"
#include <boost/archive/text_iarchive.hpp>
#include <boost/archive/text_oarchive.hpp>
#include <boost/serialization/base_object.hpp>
#include <vector>
namespace mspass::algorithms::deconvolution {
class LeastSquareDecon : public FFTDeconOperator, public ScalarDecon {
public:
  LeastSquareDecon() : FFTDeconOperator(), ScalarDecon() { damp = 1.0; };
  LeastSquareDecon(const LeastSquareDecon &parent);
  LeastSquareDecon(const mspass::utility::Metadata &md);
  LeastSquareDecon(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 damp;
  friend boost::serialization::access;
  template <class Archive>
  void serialize(Archive &ar, const unsigned int version) {
    ar &boost::serialization::base_object<FFTDeconOperator>(*this);
    ar &boost::serialization::base_object<ScalarDecon>(*this);
    ar & damp;
  }
};
} // namespace mspass::algorithms::deconvolution
#endif