thinplatespline.h

00001 /******************************************************************************
00002  * $Id: thinplatespline.h 10646 2007-01-18 02:38:10Z warmerdam $
00003  *
00004  * Project:  GDAL Warp API
00005  * Purpose:  Declarations for 2D Thin Plate Spline transformer. 
00006  * Author:   VIZRT Development Team.
00007  *
00008  * This code was provided by Gilad Ronnen (gro at visrt dot com) with 
00009  * permission to reuse under the following license.
00010  * 
00011  ******************************************************************************
00012  * Copyright (c) 2004, VIZRT Inc.
00013  *
00014  * Permission is hereby granted, free of charge, to any person obtaining a
00015  * copy of this software and associated documentation files (the "Software"),
00016  * to deal in the Software without restriction, including without limitation
00017  * the rights to use, copy, modify, merge, publish, distribute, sublicense,
00018  * and/or sell copies of the Software, and to permit persons to whom the
00019  * Software is furnished to do so, subject to the following conditions:
00020  *
00021  * The above copyright notice and this permission notice shall be included
00022  * in all copies or substantial portions of the Software.
00023  *
00024  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
00025  * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
00026  * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
00027  * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
00028  * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
00029  * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
00030  * DEALINGS IN THE SOFTWARE.
00031  ****************************************************************************/
00032 
00033 #include "gdal_alg.h"
00034 #include "cpl_conv.h"
00035 
00036 typedef enum
00037 {
00038         VIZ_GEOREF_SPLINE_ZERO_POINTS,
00039         VIZ_GEOREF_SPLINE_ONE_POINT,
00040         VIZ_GEOREF_SPLINE_TWO_POINTS,
00041         VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL,
00042         VIZ_GEOREF_SPLINE_FULL,
00043         
00044         VIZ_GEOREF_SPLINE_POINT_WAS_ADDED,
00045         VIZ_GEOREF_SPLINE_POINT_WAS_DELETED
00046 
00047 } vizGeorefInterType;
00048 
00049 //#define VIZ_GEOREF_SPLINE_MAX_POINTS 40
00050 #define VIZGEOREF_MAX_VARS 2
00051 
00052 class VizGeorefSpline2D
00053 {
00054   public:
00055 
00056     VizGeorefSpline2D(int nof_vars = 1){
00057         x = y = u = NULL;
00058         unused = index = NULL;
00059         for( int i = 0; i < nof_vars; i++ )
00060         {
00061             rhs[i] = NULL;
00062             coef[i] = NULL;
00063         }
00064           
00065         _tx = _ty = 0.0;                
00066         _ta = 10.0;
00067         _nof_points = 0;
00068         _nof_vars = nof_vars;
00069         _max_nof_points = 0;
00070         _AA = NULL;
00071         _Ainv = NULL;
00072         grow_points();
00073         for ( int v = 0; v < _nof_vars; v++ )
00074             for ( int i = 0; i < 3; i++ )
00075                 rhs[i][v] = 0.0;
00076         type = VIZ_GEOREF_SPLINE_ZERO_POINTS;
00077     }
00078 
00079     ~VizGeorefSpline2D(){
00080         if ( _AA )
00081             delete _AA;
00082         if ( _Ainv )
00083             delete _Ainv;
00084 
00085         CPLFree( x );
00086         CPLFree( y );
00087         CPLFree( u );
00088         CPLFree( unused );
00089         CPLFree( index );
00090         for( int i = 0; i < _nof_vars; i++ )
00091         {
00092             CPLFree( rhs[i] );
00093             CPLFree( coef[i] );
00094         }
00095     }
00096 
00097     int get_nof_points(){
00098         return _nof_points;
00099     }
00100 
00101     void set_toler( double tx, double ty ){
00102         _tx = tx;
00103         _ty = ty;
00104     }
00105 
00106     void get_toler( double& tx, double& ty) {
00107         tx = _tx;
00108         ty = _ty;
00109     }
00110 
00111     vizGeorefInterType get_interpolation_type ( ){
00112         return type;
00113     }
00114 
00115     void dump_data_points()
00116         {
00117             for ( int i = 0; i < _nof_points; i++ )
00118             {
00119                 fprintf(stderr, "X = %f Y = %f Vars = ", x[i], y[i]);
00120                 for ( int v = 0; v < _nof_vars; v++ )
00121                     fprintf(stderr, "%f ", rhs[i+3][v]);
00122                 fprintf(stderr, "\n");
00123             }
00124         }
00125     int delete_list()
00126         {
00127             _nof_points = 0;
00128             type = VIZ_GEOREF_SPLINE_ZERO_POINTS;
00129             if ( _AA )
00130             {
00131                 delete _AA;
00132                 _AA = NULL;
00133             }
00134             if ( _Ainv )
00135             {
00136                 delete _Ainv;
00137                 _Ainv = NULL;
00138             }
00139             return _nof_points;
00140         }
00141 
00142     void grow_points();
00143     int add_point( const double Px, const double Py, const double *Pvars );
00144     int delete_point(const double Px, const double Py );
00145     int get_point( const double Px, const double Py, double *Pvars );
00146     bool get_xy(int index, double& x, double& y);
00147     bool change_point(int index, double x, double y, double* Pvars);
00148     void reset(void) { _nof_points = 0; }
00149     int solve(void);
00150 
00151   private:      
00152     double base_func( const double x1, const double y1,
00153                       const double x2, const double y2 );
00154 
00155     vizGeorefInterType type;
00156 
00157     int _nof_vars;
00158     int _nof_points;
00159     int _max_nof_points;
00160     int _nof_eqs;
00161 
00162     double _tx, _ty;
00163     double _ta;
00164     double _dx, _dy;
00165 
00166     double *x; // [VIZ_GEOREF_SPLINE_MAX_POINTS+3];
00167     double *y; // [VIZ_GEOREF_SPLINE_MAX_POINTS+3];
00168 
00169 //    double rhs[VIZ_GEOREF_SPLINE_MAX_POINTS+3][VIZGEOREF_MAX_VARS];
00170 //    double coef[VIZ_GEOREF_SPLINE_MAX_POINTS+3][VIZGEOREF_MAX_VARS];
00171     double *rhs[VIZGEOREF_MAX_VARS];
00172     double *coef[VIZGEOREF_MAX_VARS];
00173 
00174     double *u; // [VIZ_GEOREF_SPLINE_MAX_POINTS];
00175     int *unused; // [VIZ_GEOREF_SPLINE_MAX_POINTS];
00176     int *index; // [VIZ_GEOREF_SPLINE_MAX_POINTS];
00177         
00178     double *_AA, *_Ainv;
00179 };
00180 
00181 

Generated for GDAL by doxygen 1.4.7.