📄 qwt__curve__fitter_8cpp-source.html
字号:
<a name="l00110"></a>00110 }<a name="l00111"></a>00111 <a name="l00112"></a>00112 <span class="keywordflow">if</span> ( fitMode == ParametricSpline )<a name="l00113"></a>00113 <span class="keywordflow">return</span> fitParametric(points);<a name="l00114"></a>00114 <span class="keywordflow">else</span><a name="l00115"></a>00115 <span class="keywordflow">return</span> fitSpline(points);<a name="l00116"></a>00116 }<a name="l00117"></a>00117 <a name="l00118"></a>00118 <span class="preprocessor">#if QT_VERSION < 0x040000</span><a name="l00119"></a>00119 <span class="preprocessor"></span>QwtArray<QwtDoublePoint> QwtSplineCurveFitter::fitSpline(<a name="l00120"></a>00120 <span class="keyword">const</span> QwtArray<QwtDoublePoint> &points) <span class="keyword">const</span><a name="l00121"></a>00121 <span class="preprocessor">#else</span><a name="l00122"></a>00122 <span class="preprocessor"></span>QPolygonF QwtSplineCurveFitter::fitSpline(<a name="l00123"></a>00123 <span class="keyword">const</span> QPolygonF &points) <span class="keyword">const</span><a name="l00124"></a>00124 <span class="preprocessor">#endif</span><a name="l00125"></a>00125 <span class="preprocessor"></span>{<a name="l00126"></a>00126 d_data->spline.setPoints(points);<a name="l00127"></a>00127 <span class="keywordflow">if</span> ( !d_data->spline.isValid() )<a name="l00128"></a>00128 <span class="keywordflow">return</span> points;<a name="l00129"></a>00129 <a name="l00130"></a>00130 <span class="preprocessor">#if QT_VERSION < 0x040000</span><a name="l00131"></a>00131 <span class="preprocessor"></span> QwtArray<QwtDoublePoint> fittedPoints(d_data->splineSize);<a name="l00132"></a>00132 <span class="preprocessor">#else</span><a name="l00133"></a>00133 <span class="preprocessor"></span> QPolygonF fittedPoints(d_data->splineSize);<a name="l00134"></a>00134 <span class="preprocessor">#endif</span><a name="l00135"></a>00135 <span class="preprocessor"></span><a name="l00136"></a>00136 <span class="keyword">const</span> <span class="keywordtype">double</span> x1 = points[0].x();<a name="l00137"></a>00137 <span class="keyword">const</span> <span class="keywordtype">double</span> x2 = points[int(points.size() - 1)].x();<a name="l00138"></a>00138 <span class="keyword">const</span> <span class="keywordtype">double</span> dx = x2 - x1;<a name="l00139"></a>00139 <span class="keyword">const</span> <span class="keywordtype">double</span> delta = dx / (d_data->splineSize - 1);<a name="l00140"></a>00140 <a name="l00141"></a>00141 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i < d_data->splineSize; i++)<a name="l00142"></a>00142 {<a name="l00143"></a>00143 QwtDoublePoint &p = fittedPoints[i];<a name="l00144"></a>00144 <a name="l00145"></a>00145 <span class="keyword">const</span> <span class="keywordtype">double</span> v = x1 + i * delta;<a name="l00146"></a>00146 <span class="keyword">const</span> <span class="keywordtype">double</span> sv = d_data->spline.value(v);<a name="l00147"></a>00147 <a name="l00148"></a>00148 p.setX(qRound(v));<a name="l00149"></a>00149 p.setY(qRound(sv));<a name="l00150"></a>00150 }<a name="l00151"></a>00151 d_data->spline.reset();<a name="l00152"></a>00152 <a name="l00153"></a>00153 <span class="keywordflow">return</span> fittedPoints;<a name="l00154"></a>00154 }<a name="l00155"></a>00155 <a name="l00156"></a>00156 <span class="preprocessor">#if QT_VERSION < 0x040000</span><a name="l00157"></a>00157 <span class="preprocessor"></span>QwtArray<QwtDoublePoint> QwtSplineCurveFitter::fitParametric(<a name="l00158"></a>00158 <span class="keyword">const</span> QwtArray<QwtDoublePoint> &points) <span class="keyword">const</span><a name="l00159"></a>00159 <span class="preprocessor">#else</span><a name="l00160"></a>00160 <span class="preprocessor"></span>QPolygonF QwtSplineCurveFitter::fitParametric(<a name="l00161"></a>00161 <span class="keyword">const</span> QPolygonF &points) <span class="keyword">const</span><a name="l00162"></a>00162 <span class="preprocessor">#endif</span><a name="l00163"></a>00163 <span class="preprocessor"></span>{<a name="l00164"></a>00164 <span class="keywordtype">int</span> i;<a name="l00165"></a>00165 <span class="keyword">const</span> <span class="keywordtype">int</span> size = points.size();<a name="l00166"></a>00166 <a name="l00167"></a>00167 <span class="preprocessor">#if QT_VERSION < 0x040000</span><a name="l00168"></a>00168 <span class="preprocessor"></span> QwtArray<QwtDoublePoint> fittedPoints(d_data->splineSize);<a name="l00169"></a>00169 QwtArray<QwtDoublePoint> splinePointsX(size);<a name="l00170"></a>00170 QwtArray<QwtDoublePoint> splinePointsY(size);<a name="l00171"></a>00171 <span class="preprocessor">#else</span><a name="l00172"></a>00172 <span class="preprocessor"></span> QPolygonF fittedPoints(d_data->splineSize);<a name="l00173"></a>00173 QPolygonF splinePointsX(size);<a name="l00174"></a>00174 QPolygonF splinePointsY(size);<a name="l00175"></a>00175 <span class="preprocessor">#endif</span><a name="l00176"></a>00176 <span class="preprocessor"></span><a name="l00177"></a>00177 <span class="keyword">const</span> QwtDoublePoint *p = points.data();<a name="l00178"></a>00178 QwtDoublePoint *spX = splinePointsX.data();<a name="l00179"></a>00179 QwtDoublePoint *spY = splinePointsY.data();<a name="l00180"></a>00180 <a name="l00181"></a>00181 <span class="keywordtype">double</span> param = 0.0;<a name="l00182"></a>00182 <span class="keywordflow">for</span> (i = 0; i < size; i++)<a name="l00183"></a>00183 {<a name="l00184"></a>00184 <span class="keyword">const</span> <span class="keywordtype">double</span> x = p[i].x();<a name="l00185"></a>00185 <span class="keyword">const</span> <span class="keywordtype">double</span> y = p[i].y();<a name="l00186"></a>00186 <span class="keywordflow">if</span> ( i > 0 )<a name="l00187"></a>00187 {<a name="l00188"></a>00188 <span class="keyword">const</span> <span class="keywordtype">double</span> delta = sqrt( qwtSqr(x - spX[i-1].y())<a name="l00189"></a>00189 + qwtSqr( y - spY[i-1].y() ) );<a name="l00190"></a>00190 param += qwtMax(delta, 1.0);<a name="l00191"></a>00191 }<a name="l00192"></a>00192 spX[i].setX(param);<a name="l00193"></a>00193 spX[i].setY(x);<a name="l00194"></a>00194 spY[i].setX(param);<a name="l00195"></a>00195 spY[i].setY(y);<a name="l00196"></a>00196 }<a name="l00197"></a>00197 <a name="l00198"></a>00198 d_data->spline.setPoints(splinePointsX);<a name="l00199"></a>00199 <span class="keywordflow">if</span> ( !d_data->spline.isValid() )<a name="l00200"></a>00200 <span class="keywordflow">return</span> points;<a name="l00201"></a>00201 <a name="l00202"></a>00202 <span class="keyword">const</span> <span class="keywordtype">double</span> deltaX =<a name="l00203"></a>00203 splinePointsX[size - 1].x() / (d_data->splineSize-1);<a name="l00204"></a>00204 <span class="keywordflow">for</span> (i = 0; i < d_data->splineSize; i++)<a name="l00205"></a>00205 {<a name="l00206"></a>00206 <span class="keyword">const</span> <span class="keywordtype">double</span> dtmp = i * deltaX;<a name="l00207"></a>00207 fittedPoints[i].setX(qRound(d_data->spline.value(dtmp)));<a name="l00208"></a>00208 }<a name="l00209"></a>00209 <a name="l00210"></a>00210 d_data->spline.setPoints(splinePointsY);<a name="l00211"></a>00211 <span class="keywordflow">if</span> ( !d_data->spline.isValid() )<a name="l00212"></a>00212 <span class="keywordflow">return</span> points;<a name="l00213"></a>00213 <a name="l00214"></a>00214 <span class="keyword">const</span> <span class="keywordtype">double</span> deltaY =<a name="l00215"></a>00215 splinePointsY[size - 1].x() / (d_data->splineSize-1);<a name="l00216"></a>00216 <span class="keywordflow">for</span> (i = 0; i < d_data->splineSize; i++)<a name="l00217"></a>00217 {<a name="l00218"></a>00218 <span class="keyword">const</span> <span class="keywordtype">double</span> dtmp = i * deltaY;<a name="l00219"></a>00219 fittedPoints[i].setY(qRound(d_data->spline.value(dtmp)));<a name="l00220"></a>00220 }<a name="l00221"></a>00221 <a name="l00222"></a>00222 <span class="keywordflow">return</span> fittedPoints;<a name="l00223"></a>00223 }</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 <a href="http://www.doxygen.org/index.html"><img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.4.6 </small></address></body></html>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -