QtiPlot  0.9.8.2
SmoothFilter.h
Go to the documentation of this file.
1 /***************************************************************************
2  File : SmoothFilter.h
3  Project : QtiPlot
4  --------------------------------------------------------------------
5  Copyright : (C) 2007 by Ion Vasilief
6  Email (use @ for *) : ion_vasilief*yahoo.fr
7  Description : Numerical smoothing of data sets
8 
9  ***************************************************************************/
10 
11 /***************************************************************************
12  * *
13  * This program is free software; you can redistribute it and/or modify *
14  * it under the terms of the GNU General Public License as published by *
15  * the Free Software Foundation; either version 2 of the License, or *
16  * (at your option) any later version. *
17  * *
18  * This program is distributed in the hope that it will be useful, *
19  * but WITHOUT ANY WARRANTY; without even the implied warranty of *
20  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
21  * GNU General Public License for more details. *
22  * *
23  * You should have received a copy of the GNU General Public License *
24  * along with this program; if not, write to the Free Software *
25  * Foundation, Inc., 51 Franklin Street, Fifth Floor, *
26  * Boston, MA 02110-1301 USA *
27  * *
28  ***************************************************************************/
29 #ifndef SMOOTHFILTER_H
30 #define SMOOTHFILTER_H
31 
32 #include "Filter.h"
33 #include <gsl/gsl_matrix_double.h>
34 
35 class SmoothFilter : public Filter
36 {
37 Q_OBJECT
38 
39 public:
40  SmoothFilter(ApplicationWindow *parent, QwtPlotCurve *c, int m = 3);
41  SmoothFilter(ApplicationWindow *parent, QwtPlotCurve *c, double start, double end, int m = 3);
42  SmoothFilter(ApplicationWindow *parent, Graph *g, const QString& curveTitle, int m = 3);
43  SmoothFilter(ApplicationWindow *parent, Graph *g, const QString& curveTitle, double start, double end, int m = 3);
44  SmoothFilter(ApplicationWindow *parent, Table *t, const QString& xCol, const QString& yCol, int start = 0, int end = -1, int m = 3);
45 
46  enum SmoothMethod{SavitzkyGolay = 1, FFT = 2, Average = 3, Lowess = 4};
47 
48  int method(){return (int)d_method;};
49  void setMethod(int m);
50 
51  void setSmoothPoints(int points, int left_points = 0);
53  void setPolynomOrder(int order);
54  void setLowessParameter(double f, int iterations);
55 
56 private:
57  void init(int m);
58  void calculateOutputData(double *x, double *y);
59  void smoothFFT(double *x, double *y);
60  void smoothAverage(double *x, double *y);
61  void smoothSavGol(double *x, double *y);
62  void smoothLowess(double *x, double *y);
63  static int savitzkyGolayCoefficients(int points, int polynom_order, gsl_matrix *h);
64 
67 
70 
73 
76 
78  double d_f;
79 
82 };
83 
84 #endif
Definition: SmoothFilter.h:46
SmoothMethod d_method
The smooth method.
Definition: SmoothFilter.h:66
void smoothFFT(double *x, double *y)
Definition: SmoothFilter.cpp:127
int method()
Definition: SmoothFilter.h:48
A 2D-plotting widget.
Definition: Graph.h:123
static int savitzkyGolayCoefficients(int points, int polynom_order, gsl_matrix *h)
Compute Savitzky-Golay coefficients and store them into #h.
Definition: SmoothFilter.cpp:200
double d_f
Parameter f for the Lowess algorithm.
Definition: SmoothFilter.h:78
void setSmoothPoints(int points, int left_points=0)
Definition: SmoothFilter.cpp:333
void setPolynomOrder(int order)
Sets the polynomial order in the Savitky-Golay algorithm.
Definition: SmoothFilter.cpp:351
Abstract base class for data analysis operations.
Definition: Filter.h:43
SmoothMethod
Definition: SmoothFilter.h:46
void calculateOutputData(double *x, double *y)
Calculates the data for the output curve and store it in the X an Y vectors.
Definition: SmoothFilter.cpp:99
Definition: SmoothFilter.h:46
void smoothAverage(double *x, double *y)
Definition: SmoothFilter.cpp:149
void init()
Definition: Filter.cpp:79
void setLowessParameter(double f, int iterations)
Definition: SmoothFilter.cpp:368
int d_polynom_order
Polynomial order in the Savitky-Golay algorithm (see Numerical Receipes in C for details).
Definition: SmoothFilter.h:75
QtiPlot's main window.
Definition: ApplicationWindow.h:133
Definition: FFT.h:34
void setMethod(int m)
Definition: SmoothFilter.cpp:88
Definition: SmoothFilter.h:35
int d_sav_gol_points
The number of left adjacents points used by the Savitky-Golay algorithm.
Definition: SmoothFilter.h:72
int d_smooth_points
The number of adjacents points used to smooth the data set.
Definition: SmoothFilter.h:69
Definition: SmoothFilter.h:46
SmoothFilter(ApplicationWindow *parent, QwtPlotCurve *c, int m=3)
Definition: SmoothFilter.cpp:41
double * x()
Returns the x values of the input data set.
Definition: Filter.h:98
double * y()
Returns the y values of the input data set.
Definition: Filter.h:100
void smoothSavGol(double *x, double *y)
Savitzky-Golay smoothing of (uniformly distributed) data.
Definition: SmoothFilter.cpp:261
MDI window providing a spreadsheet table with column logic.
Definition: Table.h:57
int d_iterations
Number of iterations for the Lowess algorithm.
Definition: SmoothFilter.h:81
void smoothLowess(double *x, double *y)
Definition: SmoothFilter.cpp:393