SavitzkyGolayFilter.h

Go to the documentation of this file.

#ifndef GRT_SAVITZKY_GOLAY_FILTER_HEADER
#define GRT_SAVITZKY_GOLAY_FILTER_HEADER

#include "../CoreModules/PreProcessing.h"
#include "../Util/LUDecomposition.h"

GRT_BEGIN_NAMESPACE

class GRT_API SavitzkyGolayFilter : public PreProcessing{
public:
 SavitzkyGolayFilter(const UINT numLeftHandPoints=10,const UINT numRightHandPoints=10,const UINT derivativeOrder=0,const UINT smoothingPolynomialOrder=2,const UINT numDimensions = 1);
 
 SavitzkyGolayFilter(const SavitzkyGolayFilter &rhs);
 
 virtual ~SavitzkyGolayFilter();
 
 SavitzkyGolayFilter& operator=(const SavitzkyGolayFilter &rhs);
 
 virtual bool deepCopyFrom(const PreProcessing *preProcessing);
 
 virtual bool process(const VectorFloat &inputVector);
 
 virtual bool reset();
 
 virtual bool save(std::fstream &file) const;
 
 virtual bool load(std::fstream &file);
 
 bool init(const UINT numLeftHandPoints,const UINT numRightHandPoints,const UINT derivativeOrder,const UINT smoothingPolynomialOrder,const UINT numDimensions);
 
 Float filter(const Float x);
 
 VectorFloat filter(const VectorFloat &x);
 
 VectorFloat getFilteredData() const;

 static std::string getId();
 
 //Tell the compiler we are using the following functions from the MLBase class to stop hidden virtual function warnings
 using MLBase::save;
 using MLBase::load;
 
protected:
 inline int min_(int a,int b) {return b < a ? (b) : (a);}
 inline Float min_(Float a,Float b) {return b < a ? (b) : (a);}
 bool calCoeff();
 
 UINT numPoints; //The physical length of the output array
 UINT numLeftHandPoints; //Num of leftward (past) points to use
 UINT numRightHandPoints; //Num of rightward (future) points to use
 UINT derivativeOrder; //Order of the derivative desired
 UINT smoothingPolynomialOrder; //Order of smoothing polynomial
 CircularBuffer< VectorFloat > data; //A buffer to hold the input data
 VectorFloat yy; //The filtered values
 VectorFloat coeff; //Buffer for the filter coefficients
 
private:
 static const std::string id; 
 static RegisterPreProcessingModule< SavitzkyGolayFilter > registerModule;
 
};

GRT_END_NAMESPACE

#endif //GRT_SAVITZKY_GOLAY_FILTER_HEADER