ZeroBufferFit.h

Go to the documentation of this file.

#ifndef ZeroBufferFit_h
#define ZeroBufferFit_h

#include <string>
#include <vector>

#include "HiCalTypes.h"
#include "HiCalUtil.h"
#include "HiCalConf.h"
#include "NonLinearLSQ.h"
#include "Module.h"
#include "IException.h"

namespace Isis {

  class ZeroBufferFit : public NonLinearLSQ, public Module {

    public:
    ZeroBufferFit(const HiCalConf &conf);
    virtual ~ZeroBufferFit() { }

    void setBin(int bin) { _timet.setBin(bin); }

    void setLineTime(double ltime) { _timet.setLineTime(ltime); }

     inline int size() const { return (_data.dim()); }

    int nSize()  const { return (_b2.dim()); }

    int nParms() const { return (4); }

    void setabsErr(double absError) { _absErr = absError; }
    void setrelErr(double relError) { _relErr = relError; }
    double absErr() const { return (_absErr); }
    double relErr() const { return (_relErr); }

    HiVector Solve(const HiVector &d);
    NLVector guess();
    int checkIteration(const int Iter, const NLVector &fitcoefs, 
                       const NLVector &uncerts, double cplxconj,
                       int Istatus);

    NLVector f_x(const NLVector &a);
    NLMatrix df_x(const NLVector &a);

    double Chisq() const { return (_chisq); }
    int    DoF()   const { return (nSize() - nParms()); }

    HiVector Yfit() const; 
    HiVector Normalize(const HiVector &v);

  private:
     HiLineTimeEqn _timet;    //  This is the X data set
     HiVector      _data;     //  Typically will be the HiRISE buffer data
     HiVector      _b2;       //  Data buffer used in fitting
     double        _absErr;   //  Absolute error convergence test
     double        _relErr;   //  Relative error convergence test
     double        _maxLog;   //  Maximum log value to constrain
     int           _badLines; //  Exclude lines at end of buffer
     int           _sWidth;   //  Width of guestimate filter
     int           _sIters;   //  Filter interations
     bool          _skipFit;  //  Skip fitting and pass input through
     bool          _useLinFit; // Use linear fit on failure of LM, else Zf
     int           _minLines; //  Minimum number of lines to fit (Default: 100)
     HiVector      _cc;       //  Parameter of 2-D fit
     HiVector      _guess;    //  Initial guestimate of solutions
     HiVector      _coefs;    //  Coefficients of solution
     HiVector      _uncert;   //  Uncertanties
     double        _chisq;    //  ChiSq of NonLinear equation

     HiVector poly_fit(const HiVector &d, const double line0 = 0.0) const;
     virtual void printOn(std::ostream &o) const;
     int goodLines(const HiVector &d) const { return (d.dim() - _badLines); }
     bool gotGoodLines(const HiVector &d) const { 
       return (goodLines(d) >= _minLines);
     }
  };               

}     // namespace Isis
#endif