<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<title>GDAL: thinplatespline.h Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<link href="doxygen.css" rel="stylesheet" type="text/css"/>
</head>
<body>
<!-- Generated by Doxygen 1.6.2-20100208 -->
<div class="navigation" id="top">
<div class="tabs">
<ul>
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="pages.html"><span>Related Pages</span></a></li>
<li><a href="annotated.html"><span>Classes</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
</ul>
</div>
<div class="tabs">
<ul>
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>File Members</span></a></li>
</ul>
</div>
<h1>thinplatespline.h</h1><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/******************************************************************************</span>
<a name="l00002"></a>00002 <span class="comment"> * $Id: thinplatespline.h 14122 2008-03-30 10:33:39Z rouault $</span>
<a name="l00003"></a>00003 <span class="comment"> *</span>
<a name="l00004"></a>00004 <span class="comment"> * Project: GDAL Warp API</span>
<a name="l00005"></a>00005 <span class="comment"> * Purpose: Declarations for 2D Thin Plate Spline transformer. </span>
<a name="l00006"></a>00006 <span class="comment"> * Author: VIZRT Development Team.</span>
<a name="l00007"></a>00007 <span class="comment"> *</span>
<a name="l00008"></a>00008 <span class="comment"> * This code was provided by Gilad Ronnen (gro at visrt dot com) with </span>
<a name="l00009"></a>00009 <span class="comment"> * permission to reuse under the following license.</span>
<a name="l00010"></a>00010 <span class="comment"> * </span>
<a name="l00011"></a>00011 <span class="comment"> ******************************************************************************</span>
<a name="l00012"></a>00012 <span class="comment"> * Copyright (c) 2004, VIZRT Inc.</span>
<a name="l00013"></a>00013 <span class="comment"> *</span>
<a name="l00014"></a>00014 <span class="comment"> * Permission is hereby granted, free of charge, to any person obtaining a</span>
<a name="l00015"></a>00015 <span class="comment"> * copy of this software and associated documentation files (the "Software"),</span>
<a name="l00016"></a>00016 <span class="comment"> * to deal in the Software without restriction, including without limitation</span>
<a name="l00017"></a>00017 <span class="comment"> * the rights to use, copy, modify, merge, publish, distribute, sublicense,</span>
<a name="l00018"></a>00018 <span class="comment"> * and/or sell copies of the Software, and to permit persons to whom the</span>
<a name="l00019"></a>00019 <span class="comment"> * Software is furnished to do so, subject to the following conditions:</span>
<a name="l00020"></a>00020 <span class="comment"> *</span>
<a name="l00021"></a>00021 <span class="comment"> * The above copyright notice and this permission notice shall be included</span>
<a name="l00022"></a>00022 <span class="comment"> * in all copies or substantial portions of the Software.</span>
<a name="l00023"></a>00023 <span class="comment"> *</span>
<a name="l00024"></a>00024 <span class="comment"> * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS</span>
<a name="l00025"></a>00025 <span class="comment"> * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,</span>
<a name="l00026"></a>00026 <span class="comment"> * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL</span>
<a name="l00027"></a>00027 <span class="comment"> * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER</span>
<a name="l00028"></a>00028 <span class="comment"> * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING</span>
<a name="l00029"></a>00029 <span class="comment"> * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER</span>
<a name="l00030"></a>00030 <span class="comment"> * DEALINGS IN THE SOFTWARE.</span>
<a name="l00031"></a>00031 <span class="comment"> ****************************************************************************/</span>
<a name="l00032"></a>00032
<a name="l00033"></a>00033 <span class="preprocessor">#include "<a class="code" href="gdal__alg_8h.html">gdal_alg.h</a>"</span>
<a name="l00034"></a>00034 <span class="preprocessor">#include "<a class="code" href="cpl__conv_8h.html">cpl_conv.h</a>"</span>
<a name="l00035"></a>00035
<a name="l00036"></a>00036 <span class="keyword">typedef</span> <span class="keyword">enum</span>
<a name="l00037"></a>00037 {
<a name="l00038"></a>00038 VIZ_GEOREF_SPLINE_ZERO_POINTS,
<a name="l00039"></a>00039 VIZ_GEOREF_SPLINE_ONE_POINT,
<a name="l00040"></a>00040 VIZ_GEOREF_SPLINE_TWO_POINTS,
<a name="l00041"></a>00041 VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL,
<a name="l00042"></a>00042 VIZ_GEOREF_SPLINE_FULL,
<a name="l00043"></a>00043
<a name="l00044"></a>00044 VIZ_GEOREF_SPLINE_POINT_WAS_ADDED,
<a name="l00045"></a>00045 VIZ_GEOREF_SPLINE_POINT_WAS_DELETED
<a name="l00046"></a>00046
<a name="l00047"></a>00047 } vizGeorefInterType;
<a name="l00048"></a>00048
<a name="l00049"></a>00049 <span class="comment">//#define VIZ_GEOREF_SPLINE_MAX_POINTS 40</span>
<a name="l00050"></a>00050 <span class="preprocessor">#define VIZGEOREF_MAX_VARS 2</span>
<a name="l00051"></a>00051 <span class="preprocessor"></span>
<a name="l00052"></a><a class="code" href="classVizGeorefSpline2D.html">00052</a> <span class="keyword">class </span><a class="code" href="classVizGeorefSpline2D.html">VizGeorefSpline2D</a>
<a name="l00053"></a>00053 {
<a name="l00054"></a>00054 <span class="keyword">public</span>:
<a name="l00055"></a>00055
<a name="l00056"></a>00056 <a class="code" href="classVizGeorefSpline2D.html">VizGeorefSpline2D</a>(<span class="keywordtype">int</span> nof_vars = 1){
<a name="l00057"></a>00057 x = y = u = NULL;
<a name="l00058"></a>00058 unused = index = NULL;
<a name="l00059"></a>00059 <span class="keywordflow">for</span>( <span class="keywordtype">int</span> i = 0; i < nof_vars; i++ )
<a name="l00060"></a>00060 {
<a name="l00061"></a>00061 rhs[i] = NULL;
<a name="l00062"></a>00062 coef[i] = NULL;
<a name="l00063"></a>00063 }
<a name="l00064"></a>00064
<a name="l00065"></a>00065 _tx = _ty = 0.0;
<a name="l00066"></a>00066 _ta = 10.0;
<a name="l00067"></a>00067 _nof_points = 0;
<a name="l00068"></a>00068 _nof_vars = nof_vars;
<a name="l00069"></a>00069 _max_nof_points = 0;
<a name="l00070"></a>00070 _AA = NULL;
<a name="l00071"></a>00071 _Ainv = NULL;
<a name="l00072"></a>00072 grow_points();
<a name="l00073"></a>00073 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> v = 0; v < _nof_vars; v++ )
<a name="l00074"></a>00074 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i = 0; i < 3; i++ )
<a name="l00075"></a>00075 rhs[i][v] = 0.0;
<a name="l00076"></a>00076 type = VIZ_GEOREF_SPLINE_ZERO_POINTS;
<a name="l00077"></a>00077 }
<a name="l00078"></a>00078
<a name="l00079"></a>00079 ~<a class="code" href="classVizGeorefSpline2D.html">VizGeorefSpline2D</a>(){
<a name="l00080"></a>00080 <span class="keywordflow">if</span> ( _AA )
<a name="l00081"></a>00081 CPLFree(_AA);
<a name="l00082"></a>00082 <span class="keywordflow">if</span> ( _Ainv )
<a name="l00083"></a>00083 CPLFree(_Ainv);
<a name="l00084"></a>00084
<a name="l00085"></a>00085 CPLFree( x );
<a name="l00086"></a>00086 CPLFree( y );
<a name="l00087"></a>00087 CPLFree( u );
<a name="l00088"></a>00088 CPLFree( unused );
<a name="l00089"></a>00089 CPLFree( index );
<a name="l00090"></a>00090 <span class="keywordflow">for</span>( <span class="keywordtype">int</span> i = 0; i < _nof_vars; i++ )
<a name="l00091"></a>00091 {
<a name="l00092"></a>00092 CPLFree( rhs[i] );
<a name="l00093"></a>00093 CPLFree( coef[i] );
<a name="l00094"></a>00094 }
<a name="l00095"></a>00095 }
<a name="l00096"></a>00096
<a name="l00097"></a>00097 <span class="keywordtype">int</span> get_nof_points(){
<a name="l00098"></a>00098 <span class="keywordflow">return</span> _nof_points;
<a name="l00099"></a>00099 }
<a name="l00100"></a>00100
<a name="l00101"></a>00101 <span class="keywordtype">void</span> set_toler( <span class="keywordtype">double</span> tx, <span class="keywordtype">double</span> ty ){
<a name="l00102"></a>00102 _tx = tx;
<a name="l00103"></a>00103 _ty = ty;
<a name="l00104"></a>00104 }
<a name="l00105"></a>00105
<a name="l00106"></a>00106 <span class="keywordtype">void</span> get_toler( <span class="keywordtype">double</span>& tx, <span class="keywordtype">double</span>& ty) {
<a name="l00107"></a>00107 tx = _tx;
<a name="l00108"></a>00108 ty = _ty;
<a name="l00109"></a>00109 }
<a name="l00110"></a>00110
<a name="l00111"></a>00111 vizGeorefInterType get_interpolation_type ( ){
<a name="l00112"></a>00112 <span class="keywordflow">return</span> type;
<a name="l00113"></a>00113 }
<a name="l00114"></a>00114
<a name="l00115"></a>00115 <span class="keywordtype">void</span> dump_data_points()
<a name="l00116"></a>00116 {
<a name="l00117"></a>00117 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i = 0; i < _nof_points; i++ )
<a name="l00118"></a>00118 {
<a name="l00119"></a>00119 fprintf(stderr, <span class="stringliteral">"X = %f Y = %f Vars = "</span>, x[i], y[i]);
<a name="l00120"></a>00120 <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> v = 0; v < _nof_vars; v++ )
<a name="l00121"></a>00121 fprintf(stderr, <span class="stringliteral">"%f "</span>, rhs[i+3][v]);
<a name="l00122"></a>00122 fprintf(stderr, <span class="stringliteral">"\n"</span>);
<a name="l00123"></a>00123 }
<a name="l00124"></a>00124 }
<a name="l00125"></a>00125 <span class="keywordtype">int</span> delete_list()
<a name="l00126"></a>00126 {
<a name="l00127"></a>00127 _nof_points = 0;
<a name="l00128"></a>00128 type = VIZ_GEOREF_SPLINE_ZERO_POINTS;
<a name="l00129"></a>00129 <span class="keywordflow">if</span> ( _AA )
<a name="l00130"></a>00130 {
<a name="l00131"></a>00131 CPLFree(_AA);
<a name="l00132"></a>00132 _AA = NULL;
<a name="l00133"></a>00133 }
<a name="l00134"></a>00134 <span class="keywordflow">if</span> ( _Ainv )
<a name="l00135"></a>00135 {
<a name="l00136"></a>00136 CPLFree(_Ainv);
<a name="l00137"></a>00137 _Ainv = NULL;
<a name="l00138"></a>00138 }
<a name="l00139"></a>00139 <span class="keywordflow">return</span> _nof_points;
<a name="l00140"></a>00140 }
<a name="l00141"></a>00141
<a name="l00142"></a>00142 <span class="keywordtype">void</span> grow_points();
<a name="l00143"></a>00143 <span class="keywordtype">int</span> add_point( <span class="keyword">const</span> <span class="keywordtype">double</span> Px, <span class="keyword">const</span> <span class="keywordtype">double</span> Py, <span class="keyword">const</span> <span class="keywordtype">double</span> *Pvars );
<a name="l00144"></a>00144 <span class="keywordtype">int</span> delete_point(<span class="keyword">const</span> <span class="keywordtype">double</span> Px, <span class="keyword">const</span> <span class="keywordtype">double</span> Py );
<a name="l00145"></a>00145 <span class="keywordtype">int</span> get_point( <span class="keyword">const</span> <span class="keywordtype">double</span> Px, <span class="keyword">const</span> <span class="keywordtype">double</span> Py, <span class="keywordtype">double</span> *Pvars );
<a name="l00146"></a>00146 <span class="keywordtype">bool</span> get_xy(<span class="keywordtype">int</span> index, <span class="keywordtype">double</span>& x, <span class="keywordtype">double</span>& y);
<a name="l00147"></a>00147 <span class="keywordtype">bool</span> change_point(<span class="keywordtype">int</span> index, <span class="keywordtype">double</span> x, <span class="keywordtype">double</span> y, <span class="keywordtype">double</span>* Pvars);
<a name="l00148"></a>00148 <span class="keywordtype">void</span> reset(<span class="keywordtype">void</span>) { _nof_points = 0; }
<a name="l00149"></a>00149 <span class="keywordtype">int</span> solve(<span class="keywordtype">void</span>);
<a name="l00150"></a>00150
<a name="l00151"></a>00151 <span class="keyword">private</span>:
<a name="l00152"></a>00152 <span class="keywordtype">double</span> base_func( <span class="keyword">const</span> <span class="keywordtype">double</span> x1, <span class="keyword">const</span> <span class="keywordtype">double</span> y1,
<a name="l00153"></a>00153 <span class="keyword">const</span> <span class="keywordtype">double</span> x2, <span class="keyword">const</span> <span class="keywordtype">double</span> y2 );
<a name="l00154"></a>00154
<a name="l00155"></a>00155 vizGeorefInterType type;
<a name="l00156"></a>00156
<a name="l00157"></a>00157 <span class="keywordtype">int</span> _nof_vars;
<a name="l00158"></a>00158 <span class="keywordtype">int</span> _nof_points;
<a name="l00159"></a>00159 <span class="keywordtype">int</span> _max_nof_points;
<a name="l00160"></a>00160 <span class="keywordtype">int</span> _nof_eqs;
<a name="l00161"></a>00161
<a name="l00162"></a>00162 <span class="keywordtype">double</span> _tx, _ty;
<a name="l00163"></a>00163 <span class="keywordtype">double</span> _ta;
<a name="l00164"></a>00164 <span class="keywordtype">double</span> _dx, _dy;
<a name="l00165"></a>00165
<a name="l00166"></a>00166 <span class="keywordtype">double</span> *x; <span class="comment">// [VIZ_GEOREF_SPLINE_MAX_POINTS+3];</span>
<a name="l00167"></a>00167 <span class="keywordtype">double</span> *y; <span class="comment">// [VIZ_GEOREF_SPLINE_MAX_POINTS+3];</span>
<a name="l00168"></a>00168
<a name="l00169"></a>00169 <span class="comment">// double rhs[VIZ_GEOREF_SPLINE_MAX_POINTS+3][VIZGEOREF_MAX_VARS];</span>
<a name="l00170"></a>00170 <span class="comment">// double coef[VIZ_GEOREF_SPLINE_MAX_POINTS+3][VIZGEOREF_MAX_VARS];</span>
<a name="l00171"></a>00171 <span class="keywordtype">double</span> *rhs[VIZGEOREF_MAX_VARS];
<a name="l00172"></a>00172 <span class="keywordtype">double</span> *coef[VIZGEOREF_MAX_VARS];
<a name="l00173"></a>00173
<a name="l00174"></a>00174 <span class="keywordtype">double</span> *u; <span class="comment">// [VIZ_GEOREF_SPLINE_MAX_POINTS];</span>
<a name="l00175"></a>00175 <span class="keywordtype">int</span> *unused; <span class="comment">// [VIZ_GEOREF_SPLINE_MAX_POINTS];</span>
<a name="l00176"></a>00176 <span class="keywordtype">int</span> *index; <span class="comment">// [VIZ_GEOREF_SPLINE_MAX_POINTS];</span>
<a name="l00177"></a>00177
<a name="l00178"></a>00178 <span class="keywordtype">double</span> *_AA, *_Ainv;
<a name="l00179"></a>00179 };
<a name="l00180"></a>00180
<a name="l00181"></a>00181
</pre></div></div>
<hr>
Generated for GDAL by
<a href="http://www.doxygen.org/index.html"><img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.6.2-20100208.
</body>
</html>
