QtiPlot
0.9.8.2
|
#include <SmoothFilter.h>
Public Types | |
enum | SmoothMethod { SavitzkyGolay = 1, FFT = 2, Average = 3, Lowess = 4 } |
Public Member Functions | |
int | method () |
void | setLowessParameter (double f, int iterations) |
void | setMethod (int m) |
void | setPolynomOrder (int order) |
Sets the polynomial order in the Savitky-Golay algorithm. More... | |
void | setSmoothPoints (int points, int left_points=0) |
SmoothFilter (ApplicationWindow *parent, QwtPlotCurve *c, int m=3) | |
SmoothFilter (ApplicationWindow *parent, QwtPlotCurve *c, double start, double end, int m=3) | |
SmoothFilter (ApplicationWindow *parent, Graph *g, const QString &curveTitle, int m=3) | |
SmoothFilter (ApplicationWindow *parent, Graph *g, const QString &curveTitle, double start, double end, int m=3) | |
SmoothFilter (ApplicationWindow *parent, Table *t, const QString &xCol, const QString &yCol, int start=0, int end=-1, int m=3) | |
![]() | |
int | dataSize () |
Returns the size of the input data set. More... | |
virtual void | enableGraphicsDisplay (bool on=true, Graph *g=0) |
bool | error () |
Filter (ApplicationWindow *parent, Table *t=0, const QString &name=QString()) | |
Filter (ApplicationWindow *parent, Graph *g=0, const QString &name=QString()) | |
Filter (ApplicationWindow *parent, Matrix *m, const QString &name=QString()) | |
Filter (ApplicationWindow *parent, QwtPlotCurve *c) | |
virtual QString | legendInfo () |
Output string added to the plot as a new legend. More... | |
Graph * | outputGraph () |
Returns a pointer to the graph where the result curve should be displayed. More... | |
Table * | resultTable () |
Returns a pointer to the table created to display the results. More... | |
virtual bool | run () |
Actually does the job. Should be reimplemented in derived classes. More... | |
void | setColor (int colorId) |
Obsolete: sets the color of the output fit curve. More... | |
void | setColor (const QColor &color) |
Sets the color of the output fit curve. More... | |
void | setColor (const QString &colorName) |
Sets the color of the output fit curve. Provided for convenience. To be used in scripts only! More... | |
virtual void | setDataCurve (QwtPlotCurve *curve, double start, double end) |
bool | setDataFromCurve (QwtPlotCurve *c) |
bool | setDataFromCurve (QwtPlotCurve *c, double from, double to) |
bool | setDataFromCurve (const QString &curveTitle, Graph *g=0) |
bool | setDataFromCurve (const QString &curveTitle, double from, double to, Graph *g=0) |
virtual bool | setDataFromTable (Table *, const QString &, const QString &, int=1, int=-1, bool=false) |
void | setError (bool on=true) |
void | setInterval (double from, double to) |
Changes the data range if the source curve was already assigned. Provided for convenience. More... | |
void | setMaximumIterations (int iter) |
Sets the maximum number of iterations to be performed during an iterative session. More... | |
void | setOutputPoints (int points) |
Sets the number of points in the output curve. More... | |
void | setOutputPrecision (int digits) |
Sets the precision used for the output. More... | |
void | setSortData (bool on=true) |
void | setTolerance (double eps) |
Sets the tolerance used by the GSL routines. More... | |
void | setUpdateOutputGraph (bool update=true) |
virtual void | showLegend () |
Adds a new legend to the plot. Calls virtual legendInfo() More... | |
double * | x () |
Returns the x values of the input data set. More... | |
double * | y () |
Returns the y values of the input data set. More... | |
~Filter () | |
Private Member Functions | |
void | calculateOutputData (double *x, double *y) |
Calculates the data for the output curve and store it in the X an Y vectors. More... | |
void | init (int m) |
void | smoothAverage (double *x, double *y) |
void | smoothFFT (double *x, double *y) |
void | smoothLowess (double *x, double *y) |
void | smoothSavGol (double *x, double *y) |
Savitzky-Golay smoothing of (uniformly distributed) data. More... | |
Static Private Member Functions | |
static int | savitzkyGolayCoefficients (int points, int polynom_order, gsl_matrix *h) |
Compute Savitzky-Golay coefficients and store them into #h. More... | |
Private Attributes | |
double | d_f |
Parameter f for the Lowess algorithm. More... | |
int | d_iterations |
Number of iterations for the Lowess algorithm. More... | |
SmoothMethod | d_method |
The smooth method. More... | |
int | d_polynom_order |
Polynomial order in the Savitky-Golay algorithm (see Numerical Receipes in C for details). More... | |
int | d_sav_gol_points |
The number of left adjacents points used by the Savitky-Golay algorithm. More... | |
int | d_smooth_points |
The number of adjacents points used to smooth the data set. More... | |
Additional Inherited Members | |
![]() | |
QwtPlotCurve * | addResultCurve (double *x, double *y) |
Adds the result curve to the target output plot window. Creates a hidden table and frees the input data from memory. More... | |
MultiLayer * | createOutputGraph () |
virtual int | curveData (QwtPlotCurve *c, double start, double end, double **x, double **y) |
int | curveIndex (const QString &curveTitle, Graph *g) |
Performs checks and returns the index of the source data curve if OK, -1 otherwise. More... | |
int | curveRange (QwtPlotCurve *c, double start, double end, int *iStart, int *iEnd) |
virtual void | freeMemory () |
Frees the memory allocated for the X and Y data sets. More... | |
void | init () |
virtual QString | logInfo () |
Output string added to the log pannel of the application. More... | |
void | memoryErrorMessage () |
virtual void | output () |
Performs the data analysis and takes care of the output. More... | |
virtual int | sortedCurveData (QwtPlotCurve *c, double start, double end, double **x, double **y) |
Same as curveData, but sorts the points by their x value. More... | |
![]() | |
QwtPlotCurve * | d_curve |
The curve to be analysed. More... | |
QColor | d_curveColor |
Color index of the result curve. More... | |
QString | d_explanation |
String explaining the operation in the comment of the result table and in the project explorer. More... | |
double | d_from |
Data interval. More... | |
Graph * | d_graph |
The source graph with the curve to be analyzed. More... | |
bool | d_graphics_display |
Specifies if the filter should display a result curve. More... | |
bool | d_init_err |
Error flag telling if something went wrong during the initialization phase. More... | |
Matrix * | d_matrix |
int | d_max_iterations |
Maximum number of iterations per fit. More... | |
int | d_min_points |
Minimum number of data points necessary to perform the operation. More... | |
int | d_n |
Size of the data arrays. More... | |
Graph * | d_output_graph |
The graph where the result curve should be displayed. More... | |
int | d_points |
Number of result points to de calculated and displayed in the output curve. More... | |
int | d_prec |
Precision (number of significant digits) used for the results output. More... | |
Table * | d_result_table |
The table displaying the results of the filtering operation (not alvays valid!) More... | |
bool | d_sort_data |
Specifies if the filter needs sorted data as input. More... | |
Table * | d_table |
A table source of data. More... | |
double | d_to |
double | d_tolerance |
GSL Tolerance, if ever needed... More... | |
bool | d_update_output_graph |
Specifies if the filter should update the output graph. More... | |
double * | d_x |
x data set to be analysed More... | |
double * | d_y |
y data set to be analysed More... | |
QString | d_y_col_name |
SmoothFilter::SmoothFilter | ( | ApplicationWindow * | parent, |
QwtPlotCurve * | c, | ||
int | m = 3 |
||
) |
References Filter::init(), and Filter::setDataFromCurve().
SmoothFilter::SmoothFilter | ( | ApplicationWindow * | parent, |
QwtPlotCurve * | c, | ||
double | start, | ||
double | end, | ||
int | m = 3 |
||
) |
References Filter::init(), and Filter::setDataFromCurve().
SmoothFilter::SmoothFilter | ( | ApplicationWindow * | parent, |
Graph * | g, | ||
const QString & | curveTitle, | ||
int | m = 3 |
||
) |
References Filter::init(), and Filter::setDataFromCurve().
SmoothFilter::SmoothFilter | ( | ApplicationWindow * | parent, |
Graph * | g, | ||
const QString & | curveTitle, | ||
double | start, | ||
double | end, | ||
int | m = 3 |
||
) |
References Filter::init(), and Filter::setDataFromCurve().
SmoothFilter::SmoothFilter | ( | ApplicationWindow * | parent, |
Table * | t, | ||
const QString & | xCol, | ||
const QString & | yCol, | ||
int | start = 0 , |
||
int | end = -1 , |
||
int | m = 3 |
||
) |
References Filter::init(), and Filter::setDataFromTable().
|
privatevirtual |
Calculates the data for the output curve and store it in the X an Y vectors.
Reimplemented from Filter.
References Filter::d_explanation, d_f, d_iterations, d_method, Filter::d_points, d_smooth_points, Filter::d_x, Filter::d_y, smoothAverage(), smoothFFT(), smoothLowess(), and smoothSavGol().
|
private |
References Filter::d_n, Filter::d_points, d_polynom_order, d_sav_gol_points, d_smooth_points, and setMethod().
|
inline |
References d_method.
|
staticprivate |
Compute Savitzky-Golay coefficients and store them into #h.
This function follows GSL conventions in that it writes its result into a matrix allocated by the caller and returns a non-zero result on error.
The coefficient matrix is defined as the matrix H mapping a set of input values to the values of the polynomial of order #polynom_order which minimizes squared deviations from the input values. It is computed using the formula $H=V(V^TV)^(-1)V^T$, where $V$ is the Vandermonde matrix of the point indices.
For a short description of the mathematical background, see http://www.statistics4u.info/fundstat_eng/cc_filter_savgol_math.html
References Filter::error().
Referenced by smoothSavGol().
void SmoothFilter::setLowessParameter | ( | double | f, |
int | iterations | ||
) |
References d_f, Filter::d_init_err, d_iterations, d_method, and Lowess.
Referenced by SmoothCurveDialog::smooth().
void SmoothFilter::setMethod | ( | int | m | ) |
References Filter::d_init_err, and d_method.
Referenced by init().
void SmoothFilter::setPolynomOrder | ( | int | order | ) |
Sets the polynomial order in the Savitky-Golay algorithm.
References Filter::d_init_err, d_method, d_polynom_order, d_sav_gol_points, d_smooth_points, and SavitzkyGolay.
Referenced by SmoothCurveDialog::smooth().
void SmoothFilter::setSmoothPoints | ( | int | points, |
int | left_points = 0 |
||
) |
References Filter::d_init_err, d_polynom_order, d_sav_gol_points, and d_smooth_points.
Referenced by SmoothCurveDialog::smooth().
|
private |
References Filter::d_n, and d_smooth_points.
Referenced by calculateOutputData().
|
private |
References Filter::d_n, d_smooth_points, and Filter::d_x.
Referenced by calculateOutputData().
|
private |
References clowess(), d_f, d_iterations, and Filter::d_n.
Referenced by calculateOutputData().
|
private |
Savitzky-Golay smoothing of (uniformly distributed) data.
When the data is not uniformly distributed, Savitzky-Golay looses its interesting conservation properties. On the other hand, a central point of the algorithm is that for uniform data, the operation can be implemented as a convolution.
There are at least three possible approaches to handling edges of the data vector (cutting them off, zero padding and using the left-/rightmost smoothing polynomial for computing smoothed values near the edges). Zero-padding is a particularly bad choice for signals with a distinctly non-zero baseline and cutting off edges makes further computations on the original and smoothed signals more difficult; therefore, deviating from the user-specified number of left/right adjacent points (by smoothing over a fixed window at the edges) would be the least annoying alternative;
References Filter::d_n, d_polynom_order, d_sav_gol_points, d_smooth_points, Filter::error(), Filter::memoryErrorMessage(), and savitzkyGolayCoefficients().
Referenced by calculateOutputData().
|
private |
Parameter f for the Lowess algorithm.
Referenced by calculateOutputData(), setLowessParameter(), and smoothLowess().
|
private |
Number of iterations for the Lowess algorithm.
Referenced by calculateOutputData(), setLowessParameter(), and smoothLowess().
|
private |
The smooth method.
Referenced by calculateOutputData(), method(), setLowessParameter(), setMethod(), and setPolynomOrder().
|
private |
Polynomial order in the Savitky-Golay algorithm (see Numerical Receipes in C for details).
Referenced by init(), setPolynomOrder(), setSmoothPoints(), and smoothSavGol().
|
private |
The number of left adjacents points used by the Savitky-Golay algorithm.
Referenced by init(), setPolynomOrder(), setSmoothPoints(), and smoothSavGol().
|
private |
The number of adjacents points used to smooth the data set.
Referenced by calculateOutputData(), init(), setPolynomOrder(), setSmoothPoints(), smoothAverage(), smoothFFT(), and smoothSavGol().