1
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
2
<html><head><meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1">
3
<title>Qwt User's Guide: qwt_curve_fitter.cpp Source File</title>
4
<link href="doxygen.css" rel="stylesheet" type="text/css">
5
<link href="tabs.css" rel="stylesheet" type="text/css">
7
<!-- Generated by Doxygen 1.4.6 -->
10
<li><a href="index.html"><span>Main Page</span></a></li>
11
<li><a href="classes.html"><span>Classes</span></a></li>
12
<li id="current"><a href="files.html"><span>Files</span></a></li>
13
<li><a href="pages.html"><span>Related Pages</span></a></li>
15
<h1>qwt_curve_fitter.cpp</h1><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/* -*- mode: C++ ; c-file-style: "stroustrup" -*- *****************************</span>
16
<a name="l00002"></a>00002 <span class="comment"> * Qwt Widget Library</span>
17
<a name="l00003"></a>00003 <span class="comment"> * Copyright (C) 1997 Josef Wilgen</span>
18
<a name="l00004"></a>00004 <span class="comment"> * Copyright (C) 2002 Uwe Rathmann</span>
19
<a name="l00005"></a>00005 <span class="comment"> * </span>
20
<a name="l00006"></a>00006 <span class="comment"> * This library is free software; you can redistribute it and/or</span>
21
<a name="l00007"></a>00007 <span class="comment"> * modify it under the terms of the Qwt License, Version 1.0</span>
22
<a name="l00008"></a>00008 <span class="comment"> *****************************************************************************/</span>
23
<a name="l00009"></a>00009
24
<a name="l00010"></a>00010 <span class="preprocessor">#include "qwt_math.h"</span>
25
<a name="l00011"></a>00011 <span class="preprocessor">#include "qwt_spline.h"</span>
26
<a name="l00012"></a>00012 <span class="preprocessor">#include "qwt_curve_fitter.h"</span>
27
<a name="l00013"></a>00013
28
<a name="l00015"></a><a class="code" href="class_qwt_curve_fitter.html#d804660017ae0f3ede604470724b3df3">00015</a> <a class="code" href="class_qwt_curve_fitter.html#d804660017ae0f3ede604470724b3df3">QwtCurveFitter::QwtCurveFitter</a>()
29
<a name="l00016"></a>00016 {
30
<a name="l00017"></a>00017 }
31
<a name="l00018"></a>00018
32
<a name="l00020"></a><a class="code" href="class_qwt_curve_fitter.html#62919f309f0441d94adea94573d0f778">00020</a> <a class="code" href="class_qwt_curve_fitter.html#62919f309f0441d94adea94573d0f778">QwtCurveFitter::~QwtCurveFitter</a>()
33
<a name="l00021"></a>00021 {
34
<a name="l00022"></a>00022 }
35
<a name="l00023"></a>00023
36
<a name="l00024"></a>00024 <span class="keyword">class </span>QwtSplineCurveFitter::PrivateData
37
<a name="l00025"></a>00025 {
38
<a name="l00026"></a>00026 <span class="keyword">public</span>:
39
<a name="l00027"></a>00027 PrivateData():
40
<a name="l00028"></a>00028 <a class="code" href="class_qwt_spline_curve_fitter.html#473468f2151839f3290975f6b18f1c4f">fitMode</a>(<a class="code" href="class_qwt_spline_curve_fitter.html">QwtSplineCurveFitter</a>::Auto),
41
<a name="l00029"></a>00029 <a class="code" href="class_qwt_spline_curve_fitter.html#09e5cb291c90db0aa8e6e51377bc9308">splineSize</a>(250)
42
<a name="l00030"></a>00030 {
43
<a name="l00031"></a>00031 }
44
<a name="l00032"></a>00032
45
<a name="l00033"></a>00033 <a class="code" href="class_qwt_spline.html">QwtSpline</a> <a class="code" href="class_qwt_spline_curve_fitter.html#7e5211e9e84df705d8bbbc7a57fbfa6d">spline</a>;
46
<a name="l00034"></a>00034 <a class="code" href="class_qwt_spline_curve_fitter.html#8c5e6858f885b5691c30092a950879a8">QwtSplineCurveFitter::FitMode</a> <a class="code" href="class_qwt_spline_curve_fitter.html#473468f2151839f3290975f6b18f1c4f">fitMode</a>;
47
<a name="l00035"></a>00035 <span class="keywordtype">int</span> <a class="code" href="class_qwt_spline_curve_fitter.html#09e5cb291c90db0aa8e6e51377bc9308">splineSize</a>;
48
<a name="l00036"></a>00036 };
49
<a name="l00037"></a>00037
50
<a name="l00038"></a><a class="code" href="class_qwt_spline_curve_fitter.html#98ae80240b254df85dcc44e1f3e4e830">00038</a> <a class="code" href="class_qwt_spline_curve_fitter.html#98ae80240b254df85dcc44e1f3e4e830">QwtSplineCurveFitter::QwtSplineCurveFitter</a>()
51
<a name="l00039"></a>00039 {
52
<a name="l00040"></a>00040 d_data = <span class="keyword">new</span> PrivateData;
53
<a name="l00041"></a>00041 }
54
<a name="l00042"></a>00042
55
<a name="l00043"></a><a class="code" href="class_qwt_spline_curve_fitter.html#933ecbf859698978104d0dd4ec2a2f6a">00043</a> <a class="code" href="class_qwt_spline_curve_fitter.html#933ecbf859698978104d0dd4ec2a2f6a">QwtSplineCurveFitter::~QwtSplineCurveFitter</a>()
56
<a name="l00044"></a>00044 {
57
<a name="l00045"></a>00045 <span class="keyword">delete</span> d_data;
58
<a name="l00046"></a>00046 }
59
<a name="l00047"></a>00047
60
<a name="l00048"></a><a class="code" href="class_qwt_spline_curve_fitter.html#8381be57ee16b5a2bdacafbd5d71908b">00048</a> <span class="keywordtype">void</span> <a class="code" href="class_qwt_spline_curve_fitter.html#8381be57ee16b5a2bdacafbd5d71908b">QwtSplineCurveFitter::setFitMode</a>(<a class="code" href="class_qwt_spline_curve_fitter.html#8c5e6858f885b5691c30092a950879a8">FitMode</a> mode)
61
<a name="l00049"></a>00049 {
62
<a name="l00050"></a>00050 d_data->fitMode = mode;
63
<a name="l00051"></a>00051 }
64
<a name="l00052"></a>00052
65
<a name="l00053"></a><a class="code" href="class_qwt_spline_curve_fitter.html#473468f2151839f3290975f6b18f1c4f">00053</a> <a class="code" href="class_qwt_spline_curve_fitter.html#8c5e6858f885b5691c30092a950879a8">QwtSplineCurveFitter::FitMode</a> <a class="code" href="class_qwt_spline_curve_fitter.html#473468f2151839f3290975f6b18f1c4f">QwtSplineCurveFitter::fitMode</a>()<span class="keyword"> const</span>
66
<a name="l00054"></a>00054 <span class="keyword"></span>{
67
<a name="l00055"></a>00055 <span class="keywordflow">return</span> d_data->fitMode;
68
<a name="l00056"></a>00056 }
69
<a name="l00057"></a>00057
70
<a name="l00058"></a><a class="code" href="class_qwt_spline_curve_fitter.html#7f819ad010b19d58179655e4ceb1c6f1">00058</a> <span class="keywordtype">void</span> <a class="code" href="class_qwt_spline_curve_fitter.html#7f819ad010b19d58179655e4ceb1c6f1">QwtSplineCurveFitter::setSpline</a>(<span class="keyword">const</span> <a class="code" href="class_qwt_spline.html">QwtSpline</a> &spline)
71
<a name="l00059"></a>00059 {
72
<a name="l00060"></a>00060 d_data->spline = <a class="code" href="class_qwt_spline_curve_fitter.html#7e5211e9e84df705d8bbbc7a57fbfa6d">spline</a>;
73
<a name="l00061"></a>00061 d_data->spline.reset();
74
<a name="l00062"></a>00062 }
75
<a name="l00063"></a>00063
76
<a name="l00064"></a><a class="code" href="class_qwt_spline_curve_fitter.html#7e5211e9e84df705d8bbbc7a57fbfa6d">00064</a> <span class="keyword">const</span> <a class="code" href="class_qwt_spline.html">QwtSpline</a> &<a class="code" href="class_qwt_spline_curve_fitter.html#7e5211e9e84df705d8bbbc7a57fbfa6d">QwtSplineCurveFitter::spline</a>()<span class="keyword"> const</span>
77
<a name="l00065"></a>00065 <span class="keyword"></span>{
78
<a name="l00066"></a>00066 <span class="keywordflow">return</span> d_data->spline;
79
<a name="l00067"></a>00067 }
80
<a name="l00068"></a>00068
81
<a name="l00069"></a><a class="code" href="class_qwt_spline_curve_fitter.html#c501260a25953e1ded6bbc84c3250fa8">00069</a> <a class="code" href="class_qwt_spline.html">QwtSpline</a> &<a class="code" href="class_qwt_spline_curve_fitter.html#7e5211e9e84df705d8bbbc7a57fbfa6d">QwtSplineCurveFitter::spline</a>()
82
<a name="l00070"></a>00070 {
83
<a name="l00071"></a>00071 <span class="keywordflow">return</span> d_data->spline;
84
<a name="l00072"></a>00072 }
85
<a name="l00073"></a>00073
86
<a name="l00074"></a><a class="code" href="class_qwt_spline_curve_fitter.html#8ecea80fb540b92cd22b6b0703636460">00074</a> <span class="keywordtype">void</span> <a class="code" href="class_qwt_spline_curve_fitter.html#8ecea80fb540b92cd22b6b0703636460">QwtSplineCurveFitter::setSplineSize</a>(<span class="keywordtype">int</span> splineSize)
87
<a name="l00075"></a>00075 {
88
<a name="l00076"></a>00076 d_data->splineSize = qwtMax(splineSize, 10);
89
<a name="l00077"></a>00077 }
90
<a name="l00078"></a>00078
91
<a name="l00079"></a><a class="code" href="class_qwt_spline_curve_fitter.html#09e5cb291c90db0aa8e6e51377bc9308">00079</a> <span class="keywordtype">int</span> <a class="code" href="class_qwt_spline_curve_fitter.html#09e5cb291c90db0aa8e6e51377bc9308">QwtSplineCurveFitter::splineSize</a>()<span class="keyword"> const</span>
92
<a name="l00080"></a>00080 <span class="keyword"></span>{
93
<a name="l00081"></a>00081 <span class="keywordflow">return</span> d_data->splineSize;
94
<a name="l00082"></a>00082 }
95
<a name="l00083"></a>00083
96
<a name="l00084"></a>00084 <span class="preprocessor">#if QT_VERSION < 0x040000</span>
97
<a name="l00085"></a>00085 <span class="preprocessor"></span>QwtArray<QwtDoublePoint> <a class="code" href="class_qwt_spline_curve_fitter.html#14d64180d7777d0cd9c2adf81e120140">QwtSplineCurveFitter::fitCurve</a>(
98
<a name="l00086"></a>00086 <span class="keyword">const</span> QwtArray<QwtDoublePoint> & points) <span class="keyword">const</span>
99
<a name="l00087"></a>00087 <span class="preprocessor">#else</span>
100
<a name="l00088"></a><a class="code" href="class_qwt_spline_curve_fitter.html#14d64180d7777d0cd9c2adf81e120140">00088</a> <span class="preprocessor"></span>QPolygonF <a class="code" href="class_qwt_spline_curve_fitter.html#14d64180d7777d0cd9c2adf81e120140">QwtSplineCurveFitter::fitCurve</a>(
101
<a name="l00089"></a>00089 <span class="keyword">const</span> QPolygonF &points) <span class="keyword">const</span>
102
<a name="l00090"></a>00090 <span class="preprocessor">#endif</span>
103
<a name="l00091"></a>00091 <span class="preprocessor"></span>{
104
<a name="l00092"></a>00092 <span class="keyword">const</span> <span class="keywordtype">int</span> size = (int)points.size();
105
<a name="l00093"></a>00093 <span class="keywordflow">if</span> ( size <= 2 )
106
<a name="l00094"></a>00094 <span class="keywordflow">return</span> points;
107
<a name="l00095"></a>00095
108
<a name="l00096"></a>00096 <a class="code" href="class_qwt_spline_curve_fitter.html#8c5e6858f885b5691c30092a950879a8">FitMode</a> fitMode = d_data->fitMode;
109
<a name="l00097"></a>00097 <span class="keywordflow">if</span> ( fitMode == Auto )
110
<a name="l00098"></a>00098 {
111
<a name="l00099"></a>00099 fitMode = Spline;
112
<a name="l00100"></a>00100
113
<a name="l00101"></a>00101 <span class="keyword">const</span> QwtDoublePoint *p = points.data();
114
<a name="l00102"></a>00102 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i = 1; i < size; i++ )
115
<a name="l00103"></a>00103 {
116
<a name="l00104"></a>00104 <span class="keywordflow">if</span> ( p[i].x() <= p[i-1].x() )
117
<a name="l00105"></a>00105 {
118
<a name="l00106"></a>00106 fitMode = ParametricSpline;
119
<a name="l00107"></a>00107 <span class="keywordflow">break</span>;
120
<a name="l00108"></a>00108 }
121
<a name="l00109"></a>00109 };
122
<a name="l00110"></a>00110 }
123
<a name="l00111"></a>00111
124
<a name="l00112"></a>00112 <span class="keywordflow">if</span> ( fitMode == ParametricSpline )
125
<a name="l00113"></a>00113 <span class="keywordflow">return</span> fitParametric(points);
126
<a name="l00114"></a>00114 <span class="keywordflow">else</span>
127
<a name="l00115"></a>00115 <span class="keywordflow">return</span> fitSpline(points);
128
<a name="l00116"></a>00116 }
129
<a name="l00117"></a>00117
130
<a name="l00118"></a>00118 <span class="preprocessor">#if QT_VERSION < 0x040000</span>
131
<a name="l00119"></a>00119 <span class="preprocessor"></span>QwtArray<QwtDoublePoint> QwtSplineCurveFitter::fitSpline(
132
<a name="l00120"></a>00120 <span class="keyword">const</span> QwtArray<QwtDoublePoint> &points) <span class="keyword">const</span>
133
<a name="l00121"></a>00121 <span class="preprocessor">#else</span>
134
<a name="l00122"></a>00122 <span class="preprocessor"></span>QPolygonF QwtSplineCurveFitter::fitSpline(
135
<a name="l00123"></a>00123 <span class="keyword">const</span> QPolygonF &points) <span class="keyword">const</span>
136
<a name="l00124"></a>00124 <span class="preprocessor">#endif</span>
137
<a name="l00125"></a>00125 <span class="preprocessor"></span>{
138
<a name="l00126"></a>00126 d_data->spline.setPoints(points);
139
<a name="l00127"></a>00127 <span class="keywordflow">if</span> ( !d_data->spline.isValid() )
140
<a name="l00128"></a>00128 <span class="keywordflow">return</span> points;
141
<a name="l00129"></a>00129
142
<a name="l00130"></a>00130 <span class="preprocessor">#if QT_VERSION < 0x040000</span>
143
<a name="l00131"></a>00131 <span class="preprocessor"></span> QwtArray<QwtDoublePoint> fittedPoints(d_data->splineSize);
144
<a name="l00132"></a>00132 <span class="preprocessor">#else</span>
145
<a name="l00133"></a>00133 <span class="preprocessor"></span> QPolygonF fittedPoints(d_data->splineSize);
146
<a name="l00134"></a>00134 <span class="preprocessor">#endif</span>
147
<a name="l00135"></a>00135 <span class="preprocessor"></span>
148
<a name="l00136"></a>00136 <span class="keyword">const</span> <span class="keywordtype">double</span> x1 = points[0].x();
149
<a name="l00137"></a>00137 <span class="keyword">const</span> <span class="keywordtype">double</span> x2 = points[int(points.size() - 1)].x();
150
<a name="l00138"></a>00138 <span class="keyword">const</span> <span class="keywordtype">double</span> dx = x2 - x1;
151
<a name="l00139"></a>00139 <span class="keyword">const</span> <span class="keywordtype">double</span> delta = dx / (d_data->splineSize - 1);
152
<a name="l00140"></a>00140
153
<a name="l00141"></a>00141 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i < d_data->splineSize; i++)
154
<a name="l00142"></a>00142 {
155
<a name="l00143"></a>00143 QwtDoublePoint &p = fittedPoints[i];
156
<a name="l00144"></a>00144
157
<a name="l00145"></a>00145 <span class="keyword">const</span> <span class="keywordtype">double</span> v = x1 + i * delta;
158
<a name="l00146"></a>00146 <span class="keyword">const</span> <span class="keywordtype">double</span> sv = d_data->spline.value(v);
159
<a name="l00147"></a>00147
160
<a name="l00148"></a>00148 p.setX(qRound(v));
161
<a name="l00149"></a>00149 p.setY(qRound(sv));
162
<a name="l00150"></a>00150 }
163
<a name="l00151"></a>00151 d_data->spline.reset();
164
<a name="l00152"></a>00152
165
<a name="l00153"></a>00153 <span class="keywordflow">return</span> fittedPoints;
166
<a name="l00154"></a>00154 }
167
<a name="l00155"></a>00155
168
<a name="l00156"></a>00156 <span class="preprocessor">#if QT_VERSION < 0x040000</span>
169
<a name="l00157"></a>00157 <span class="preprocessor"></span>QwtArray<QwtDoublePoint> QwtSplineCurveFitter::fitParametric(
170
<a name="l00158"></a>00158 <span class="keyword">const</span> QwtArray<QwtDoublePoint> &points) <span class="keyword">const</span>
171
<a name="l00159"></a>00159 <span class="preprocessor">#else</span>
172
<a name="l00160"></a>00160 <span class="preprocessor"></span>QPolygonF QwtSplineCurveFitter::fitParametric(
173
<a name="l00161"></a>00161 <span class="keyword">const</span> QPolygonF &points) <span class="keyword">const</span>
174
<a name="l00162"></a>00162 <span class="preprocessor">#endif</span>
175
<a name="l00163"></a>00163 <span class="preprocessor"></span>{
176
<a name="l00164"></a>00164 <span class="keywordtype">int</span> i;
177
<a name="l00165"></a>00165 <span class="keyword">const</span> <span class="keywordtype">int</span> size = points.size();
178
<a name="l00166"></a>00166
179
<a name="l00167"></a>00167 <span class="preprocessor">#if QT_VERSION < 0x040000</span>
180
<a name="l00168"></a>00168 <span class="preprocessor"></span> QwtArray<QwtDoublePoint> fittedPoints(d_data->splineSize);
181
<a name="l00169"></a>00169 QwtArray<QwtDoublePoint> splinePointsX(size);
182
<a name="l00170"></a>00170 QwtArray<QwtDoublePoint> splinePointsY(size);
183
<a name="l00171"></a>00171 <span class="preprocessor">#else</span>
184
<a name="l00172"></a>00172 <span class="preprocessor"></span> QPolygonF fittedPoints(d_data->splineSize);
185
<a name="l00173"></a>00173 QPolygonF splinePointsX(size);
186
<a name="l00174"></a>00174 QPolygonF splinePointsY(size);
187
<a name="l00175"></a>00175 <span class="preprocessor">#endif</span>
188
<a name="l00176"></a>00176 <span class="preprocessor"></span>
189
<a name="l00177"></a>00177 <span class="keyword">const</span> QwtDoublePoint *p = points.data();
190
<a name="l00178"></a>00178 QwtDoublePoint *spX = splinePointsX.data();
191
<a name="l00179"></a>00179 QwtDoublePoint *spY = splinePointsY.data();
192
<a name="l00180"></a>00180
193
<a name="l00181"></a>00181 <span class="keywordtype">double</span> param = 0.0;
194
<a name="l00182"></a>00182 <span class="keywordflow">for</span> (i = 0; i < size; i++)
195
<a name="l00183"></a>00183 {
196
<a name="l00184"></a>00184 <span class="keyword">const</span> <span class="keywordtype">double</span> x = p[i].x();
197
<a name="l00185"></a>00185 <span class="keyword">const</span> <span class="keywordtype">double</span> y = p[i].y();
198
<a name="l00186"></a>00186 <span class="keywordflow">if</span> ( i > 0 )
199
<a name="l00187"></a>00187 {
200
<a name="l00188"></a>00188 <span class="keyword">const</span> <span class="keywordtype">double</span> delta = sqrt( qwtSqr(x - spX[i-1].y())
201
<a name="l00189"></a>00189 + qwtSqr( y - spY[i-1].y() ) );
202
<a name="l00190"></a>00190 param += qwtMax(delta, 1.0);
203
<a name="l00191"></a>00191 }
204
<a name="l00192"></a>00192 spX[i].setX(param);
205
<a name="l00193"></a>00193 spX[i].setY(x);
206
<a name="l00194"></a>00194 spY[i].setX(param);
207
<a name="l00195"></a>00195 spY[i].setY(y);
208
<a name="l00196"></a>00196 }
209
<a name="l00197"></a>00197
210
<a name="l00198"></a>00198 d_data->spline.setPoints(splinePointsX);
211
<a name="l00199"></a>00199 <span class="keywordflow">if</span> ( !d_data->spline.isValid() )
212
<a name="l00200"></a>00200 <span class="keywordflow">return</span> points;
213
<a name="l00201"></a>00201
214
<a name="l00202"></a>00202 <span class="keyword">const</span> <span class="keywordtype">double</span> deltaX =
215
<a name="l00203"></a>00203 splinePointsX[size - 1].x() / (d_data->splineSize-1);
216
<a name="l00204"></a>00204 <span class="keywordflow">for</span> (i = 0; i < d_data->splineSize; i++)
217
<a name="l00205"></a>00205 {
218
<a name="l00206"></a>00206 <span class="keyword">const</span> <span class="keywordtype">double</span> dtmp = i * deltaX;
219
<a name="l00207"></a>00207 fittedPoints[i].setX(qRound(d_data->spline.value(dtmp)));
220
<a name="l00208"></a>00208 }
221
<a name="l00209"></a>00209
222
<a name="l00210"></a>00210 d_data->spline.setPoints(splinePointsY);
223
<a name="l00211"></a>00211 <span class="keywordflow">if</span> ( !d_data->spline.isValid() )
224
<a name="l00212"></a>00212 <span class="keywordflow">return</span> points;
225
<a name="l00213"></a>00213
226
<a name="l00214"></a>00214 <span class="keyword">const</span> <span class="keywordtype">double</span> deltaY =
227
<a name="l00215"></a>00215 splinePointsY[size - 1].x() / (d_data->splineSize-1);
228
<a name="l00216"></a>00216 <span class="keywordflow">for</span> (i = 0; i < d_data->splineSize; i++)
229
<a name="l00217"></a>00217 {
230
<a name="l00218"></a>00218 <span class="keyword">const</span> <span class="keywordtype">double</span> dtmp = i * deltaY;
231
<a name="l00219"></a>00219 fittedPoints[i].setY(qRound(d_data->spline.value(dtmp)));
232
<a name="l00220"></a>00220 }
233
<a name="l00221"></a>00221
234
<a name="l00222"></a>00222 <span class="keywordflow">return</span> fittedPoints;
235
<a name="l00223"></a>00223 }
236
</pre></div><hr size="1"><address style="align: right;"><small>Generated on Mon Feb 26 21:22:36 2007 for Qwt User's Guide by
237
<a href="http://www.doxygen.org/index.html">
238
<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.4.6 </small></address>