ossimThinPlateSpline Class Reference

#include <ossimThinPlateSpline.h>

Public Types
enum	vizGeorefInterType {   VIZ_GEOREF_SPLINE_ZERO_POINTS, VIZ_GEOREF_SPLINE_ONE_POINT, VIZ_GEOREF_SPLINE_TWO_POINTS, VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL,   VIZ_GEOREF_SPLINE_FULL, VIZ_GEOREF_SPLINE_POINT_WAS_ADDED, VIZ_GEOREF_SPLINE_POINT_WAS_DELETED }

Public Member Functions
	ossimThinPlateSpline (int nof_vars=1)
	~ossimThinPlateSpline ()
int	getNumberOfPoints () const
void	setToler (double tx, double ty)
void	getToler (double &tx, double &ty) const
vizGeorefInterType	getInterpolationType () const
int	deleteList ()
void	growPoints ()
int	addPoint (const double Px, const double Py, const double *Pvars)
int	deletePoint (const double Px, const double Py)
int	getPoint (const double Px, const double Py, double *Pvars) const
bool	getXy (int index, double &x, double &y) const
bool	changePoint (int index, double x, double y, double *Pvars)
void	reset (void)
int	solve (void)

Private Member Functions
double	baseFunc (const double x1, const double y1, const double x2, const double y2) const

Private Attributes
vizGeorefInterType	type
int	_nof_vars
int	_nof_points
int	_max_nof_points
int	_nof_eqs
double	_tx
double	_ty
double	_ta
double	_dx
double	_dy
std::vector< double >	x
std::vector< double >	y
std::vector< std::vector< double > >	rhs
std::vector< std::vector< double > >	coef
std::vector< double >	u
std::vector< int >	unused
std::vector< int >	index
std::vector< double >	_AA
std::vector< double >	_Ainv

Detailed Description

Definition at line 37 of file ossimThinPlateSpline.h.

Member Enumeration Documentation

vizGeorefInterType

enum ossimThinPlateSpline::vizGeorefInterType

Enumerator
VIZ_GEOREF_SPLINE_ZERO_POINTS
VIZ_GEOREF_SPLINE_ONE_POINT
VIZ_GEOREF_SPLINE_TWO_POINTS
VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL
VIZ_GEOREF_SPLINE_FULL
VIZ_GEOREF_SPLINE_POINT_WAS_ADDED
VIZ_GEOREF_SPLINE_POINT_WAS_DELETED

Definition at line 40 of file ossimThinPlateSpline.h.

   {
      VIZ_GEOREF_SPLINE_ZERO_POINTS,
      VIZ_GEOREF_SPLINE_ONE_POINT,
      VIZ_GEOREF_SPLINE_TWO_POINTS,
      VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL,
      VIZ_GEOREF_SPLINE_FULL,
      
      VIZ_GEOREF_SPLINE_POINT_WAS_ADDED,
      VIZ_GEOREF_SPLINE_POINT_WAS_DELETED
      
   } vizGeorefInterType;

Constructor & Destructor Documentation

ossimThinPlateSpline()

ossimThinPlateSpline::ossimThinPlateSpline ( int  nof_vars = 1 )

inline

Definition at line 53 of file ossimThinPlateSpline.h.

                                         {
      _tx = _ty = 0.0;      
      _ta = 10.0;
      _nof_points = 0;
      _nof_vars = nof_vars;
      _max_nof_points = 0;
      rhs.resize(_nof_vars);
      coef.resize(_nof_vars);
      growPoints();
      for ( int v = 0; v < _nof_vars; v++ )
         for ( int i = 0; i < 3; i++ )
            //            rhs[i][v] = 0.0;
            rhs[v][i] = 0.0;
      type = VIZ_GEOREF_SPLINE_ZERO_POINTS;
   }

~ossimThinPlateSpline()

ossimThinPlateSpline::~ossimThinPlateSpline ( )

inline

Definition at line 69 of file ossimThinPlateSpline.h.

References x, and y.

                          {
      _AA.clear();
      _Ainv.clear();
      x.clear();
      y.clear();
      u.clear();
      unused.clear();
      index.clear();
      rhs.clear();
      coef.clear();
   }

Member Function Documentation

addPoint()

int ossimThinPlateSpline::addPoint	(	const double	Px,
		const double	Py,
		const double *	Pvars
	)

Definition at line 77 of file ossimThinPlateSpline.cpp.

References _max_nof_points, _nof_points, _nof_vars, growPoints(), rhs, type, VIZ_GEOREF_SPLINE_POINT_WAS_ADDED, x, and y.

Referenced by ossimHistogram::fillInteriorEmptyBins().

{
   type = VIZ_GEOREF_SPLINE_POINT_WAS_ADDED;
   int i;
   
   if( _nof_points == _max_nof_points )
      growPoints();
   
   i = _nof_points;
   //A new point is added
   x[i] = Px;
   y[i] = Py;
   for ( int j = 0; j < _nof_vars; j++ )
      rhs[j][i+3] = Pvars[j];
   _nof_points++;
   return 1;
}

baseFunc()

double ossimThinPlateSpline::baseFunc ( const double  x1, const double  y1, const double  x2, const double  y2  ) const

private

Definition at line 391 of file ossimThinPlateSpline.cpp.

Referenced by getPoint(), and solve().

{
    if ( ( x1 == x2 ) && (y1 == y2 ) )
        return 0.0;
    
    double dist  = ( x2 - x1 ) * ( x2 - x1 ) + ( y2 - y1 ) * ( y2 - y1 );
    
    return dist * log( dist );  
}

changePoint()

bool ossimThinPlateSpline::changePoint	(	int	index,
		double	x,
		double	y,
		double *	Pvars
	)

Definition at line 95 of file ossimThinPlateSpline.cpp.

References _nof_points, _nof_vars, index, rhs, x, and y.

{
   if ( index < _nof_points )
   {
      int i = index;
      x[i] = Px;
      y[i] = Py;
      for ( int j = 0; j < _nof_vars; j++ )
         rhs[j][i+3] = Pvars[j];
   }
   
   return( true );
}

deleteList()

int ossimThinPlateSpline::deleteList ( )

inline

Definition at line 99 of file ossimThinPlateSpline.h.

    {
      _nof_points = 0;
      type = VIZ_GEOREF_SPLINE_ZERO_POINTS;
      _AA.clear();
      _Ainv.clear();
      return _nof_points;
    }

deletePoint()

int ossimThinPlateSpline::deletePoint	(	const double	Px,
		const double	Py
	)

Definition at line 128 of file ossimThinPlateSpline.cpp.

References _nof_points, _nof_vars, _tx, _ty, ossim::abs(), rhs, type, VIZ_GEOREF_SPLINE_POINT_WAS_DELETED, x, and y.

{
   for ( int i = 0; i < _nof_points; i++ )
   {
      if ( ( ossim::abs(Px - x[i]) <= _tx ) && ( ossim::abs(Py - y[i]) <= _ty ) )
      {
         for ( int j = i; j < _nof_points - 1; j++ )
         {
            x[j] = x[j+1];
            y[j] = y[j+1];
            for ( int k = 0; k < _nof_vars; k++ )
               rhs[k][j+3] = rhs[k][j+3+1];
         }
         _nof_points--;
         type = VIZ_GEOREF_SPLINE_POINT_WAS_DELETED;
         return(1);
      }
   }
   return(0);
}

getInterpolationType()

vizGeorefInterType ossimThinPlateSpline::getInterpolationType ( ) const

inline

Definition at line 95 of file ossimThinPlateSpline.h.

                                                   {
      return type;
   }

getNumberOfPoints()

int ossimThinPlateSpline::getNumberOfPoints ( ) const

inline

Definition at line 81 of file ossimThinPlateSpline.h.

   {
      return _nof_points;
   }

getPoint()

int ossimThinPlateSpline::getPoint	(	const double	Px,
		const double	Py,
		double *	Pvars
	)		const

Definition at line 309 of file ossimThinPlateSpline.cpp.

References _dx, _dy, _nof_points, _nof_vars, baseFunc(), coef, index, rhs, type, u, VIZ_GEOREF_SPLINE_FULL, VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL, VIZ_GEOREF_SPLINE_ONE_POINT, VIZ_GEOREF_SPLINE_POINT_WAS_ADDED, VIZ_GEOREF_SPLINE_POINT_WAS_DELETED, VIZ_GEOREF_SPLINE_TWO_POINTS, VIZ_GEOREF_SPLINE_ZERO_POINTS, x, and y.

Referenced by ossimHistogram::fillInteriorEmptyBins().

{
    int v, r;
    double tmp, Pu;
    double fact;
    int leftP=0, rightP=0, found = 0;
    
    switch ( type )
    {
      case VIZ_GEOREF_SPLINE_ZERO_POINTS :
         for ( v = 0; v < _nof_vars; v++ )
            vars[v] = 0.0;
         break;
      case VIZ_GEOREF_SPLINE_ONE_POINT :
         for ( v = 0; v < _nof_vars; v++ )
            vars[v] = rhs[v][3];
         break;
      case VIZ_GEOREF_SPLINE_TWO_POINTS :
         fact = _dx * ( Px - x[0] ) + _dy * ( Py - y[0] );
         for ( v = 0; v < _nof_vars; v++ )
            vars[v] = ( 1 - fact ) * rhs[v][3] + fact * rhs[v][4];
         break;
      case VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL :
         Pu = _dx * ( Px - x[0] ) + _dy * ( Py - y[0] );
         if ( Pu <= u[index[0]] )
         {
            leftP = index[0];
            rightP = index[1];
         }
         else if ( Pu >= u[index[_nof_points-1]] )
         {
            leftP = index[_nof_points-2];
            rightP = index[_nof_points-1];
         }
         else
         {
            for ( r = 1; !found && r < _nof_points; r++ )
            {
               leftP = index[r-1];
               rightP = index[r];                   
               if ( Pu >= u[leftP] && Pu <= u[rightP] )
                  found = 1;
            }
         }
         
         fact = ( Pu - u[leftP] ) / ( u[rightP] - u[leftP] );
         for ( v = 0; v < _nof_vars; v++ )
            vars[v] = ( 1.0 - fact ) * rhs[v][leftP+3] +
            fact * rhs[v][rightP+3];
         break;
      case VIZ_GEOREF_SPLINE_FULL :
         for ( v = 0; v < _nof_vars; v++ )
            vars[v] = coef[v][0] + coef[v][1] * Px + coef[v][2] * Py;
         
         for ( r = 0; r < _nof_points; r++ )
         {
            tmp = baseFunc( Px, Py, x[r], y[r] );
            for ( v= 0; v < _nof_vars; v++ )
               vars[v] += coef[v][r+3] * tmp;
         }
         break;
      case VIZ_GEOREF_SPLINE_POINT_WAS_ADDED :
         fprintf(stderr, " A point was added after the last solve\n");
         fprintf(stderr, " NO interpolation - return values are zero\n");
         for ( v = 0; v < _nof_vars; v++ )
            vars[v] = 0.0;
         return(0);
         break;
      case VIZ_GEOREF_SPLINE_POINT_WAS_DELETED :
         fprintf(stderr, " A point was deleted after the last solve\n");
         fprintf(stderr, " NO interpolation - return values are zero\n");
         for ( v = 0; v < _nof_vars; v++ )
            vars[v] = 0.0;
         return(0);
         break;
      default :
         return(0);
         break;
    }
    return(1);
}

getToler()

void ossimThinPlateSpline::getToler ( double &  tx, double &  ty  ) const

inline

Definition at line 90 of file ossimThinPlateSpline.h.

                                               {
      tx = _tx;
      ty = _ty;
   }

getXy()

bool ossimThinPlateSpline::getXy	(	int	index,
		double &	x,
		double &	y
	)		const

Definition at line 109 of file ossimThinPlateSpline.cpp.

References _nof_points, index, x, and y.

{
   bool ok;
   
   if ( index < _nof_points )
   {
      ok = true;
      outX = x[index];
      outY = y[index];
   }
   else
   {
      ok = false;
      outX = outY = 0.0f;
   }
   
   return(ok);
}

growPoints()

void ossimThinPlateSpline::growPoints

(

)

Definition at line 57 of file ossimThinPlateSpline.cpp.

References _max_nof_points, _nof_vars, coef, index, rhs, u, unused, x, and y.

Referenced by addPoint().

{
   int new_max = _max_nof_points*2 + 2 + 3;
   int i;
   
   x.resize( new_max );
   y.resize( new_max );
   u.resize( new_max );
   unused.resize( new_max );
   index.resize( new_max );
   for( i = 0; i < _nof_vars; i++ )
   {
      rhs[i].resize( new_max );
      coef[i].resize( new_max );
   }
   
   _max_nof_points = new_max - 3;
}

reset()

void ossimThinPlateSpline::reset ( void  )

inline

Definition at line 114 of file ossimThinPlateSpline.h.

{ _nof_points = 0; }

setToler()

void ossimThinPlateSpline::setToler ( double  tx, double  ty  )

inline

Definition at line 85 of file ossimThinPlateSpline.h.

                                        {
      _tx = tx;
      _ty = ty;
   }

solve()

int ossimThinPlateSpline::solve

(

void

)

Definition at line 149 of file ossimThinPlateSpline.cpp.

References _AA, _Ainv, _dx, _dy, _nof_eqs, _nof_points, A, baseFunc(), index, ossim::max(), ossim::min(), status, type, u, unused, VIZ_GEOREF_SPLINE_FULL, VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL, VIZ_GEOREF_SPLINE_ONE_POINT, VIZ_GEOREF_SPLINE_TWO_POINTS, VIZ_GEOREF_SPLINE_ZERO_POINTS, x, and y.

Referenced by ossimHistogram::fillInteriorEmptyBins().

{
   int r, c, v;
   int p;
    
   //   No points at all
   if ( _nof_points < 1 )
   {
      type = VIZ_GEOREF_SPLINE_ZERO_POINTS;
      return(0);
   }
    
   // Only one point
   if ( _nof_points == 1 )
   {
      type = VIZ_GEOREF_SPLINE_ONE_POINT;
      return(1);
   }
   // Just 2 points - it is necessarily 1D case
   if ( _nof_points == 2 )
   {
      _dx = x[1] - x[0];
      _dy = y[1] - y[0];     
      double fact = 1.0 / ( _dx * _dx + _dy * _dy );
      _dx *= fact;
      _dy *= fact;
        
      type = VIZ_GEOREF_SPLINE_TWO_POINTS;
      return(2);
   }
    
   // More than 2 points - first we have to check if it is 1D or 2D case
   
   double xmax = x[0], xmin = x[0], ymax = y[0], ymin = y[0];
   double delx, dely;
   double xx, yy;
   double sumx = 0.0f, sumy= 0.0f, sumx2 = 0.0f, sumy2 = 0.0f, sumxy = 0.0f;
   double SSxx, SSyy, SSxy;
    
   for ( p = 0; p < _nof_points; p++ )
   {
      xx = x[p];
      yy = y[p];
        
      xmax = ossim::max( xmax, xx );
      xmin = ossim::min( xmin, xx );
      ymax = ossim::max( ymax, yy );
      ymin = ossim::min( ymin, yy );
        
      sumx  += xx;
      sumx2 += xx * xx;
      sumy  += yy;
      sumy2 += yy * yy;
      sumxy += xx * yy;
   }
   delx = xmax - xmin;
   dely = ymax - ymin;
    
   SSxx = sumx2 - sumx * sumx / _nof_points;
   SSyy = sumy2 - sumy * sumy / _nof_points;
   SSxy = sumxy - sumx * sumy / _nof_points;
    
   if ( delx < 0.001 * dely || dely < 0.001 * delx || 
       fabs ( SSxy * SSxy / ( SSxx * SSyy ) ) > 0.99 )
   {
      int p1;
        
      type = VIZ_GEOREF_SPLINE_ONE_DIMENSIONAL;
        
      _dx = _nof_points * sumx2 - sumx * sumx;
      _dy = _nof_points * sumy2 - sumy * sumy;
      double fact = 1.0 / sqrt( _dx * _dx + _dy * _dy );
      _dx *= fact;
      _dy *= fact;
        
      for ( p = 0; p < _nof_points; p++ )
      {
         double dxp = x[p] - x[0];
         double dyp = y[p] - y[0];
         u[p] = _dx * dxp + _dy * dyp;
         unused[p] = 1;
      }
        
      for ( p = 0; p < _nof_points; p++ )
      {
         int min_index = -1;
         double min_u = 0;
         for ( p1 = 0; p1 < _nof_points; p1++ )
         {
            if ( unused[p1] )
            {
               if ( min_index < 0 || u[p1] < min_u )
               {
                  min_index = p1;
                  min_u = u[p1];
               }
            }
         }
         index[p] = min_index;
         unused[min_index] = 0;
      }
        
      return(3);
   }
    
   type = VIZ_GEOREF_SPLINE_FULL;
   
   _nof_eqs = _nof_points + 3;
    
   _AA.resize( _nof_eqs * _nof_eqs );
   _Ainv.resize(_nof_eqs * _nof_eqs );
    
   // Calc the values of the matrix A
   for ( r = 0; r < 3; r++ )
      for ( c = 0; c < 3; c++ )
         A(r,c) = 0.0;
   
   for ( c = 0; c < _nof_points; c++ )
   {
      A(0,c+3) = 1.0;
      A(1,c+3) = x[c];
      A(2,c+3) = y[c];
      
      A(c+3,0) = 1.0;
      A(c+3,1) = x[c];
      A(c+3,2) = y[c];
   }
   
   for ( r = 0; r < _nof_points; r++ )
      for ( c = r; c < _nof_points; c++ )
      {
         A(r+3,c+3) = baseFunc( x[r], y[r], x[c], y[c] );
         if ( r != c )
            A(c+3,r+3 ) = A(r+3,c+3);
      }
   
   
   // Invert the matrix
   int status = matrixInvert( _nof_eqs, &_AA.front(), &_Ainv.front() );
   
   if ( !status )
   {
      ossimNotify(ossimNotifyLevel_WARN)
         << " There is a problem to invert the interpolation matrix\n"
         << std::endl;
      return 0;
   }
   
   // calc the coefs
   for ( v = 0; v < _nof_vars; v++ )
      for ( r = 0; r < _nof_eqs; r++ )
      {
         coef[v][r] = 0.0;
         for ( c = 0; c < _nof_eqs; c++ )
            coef[v][r] += Ainv(r,c) * rhs[v][c];
      }
   
   return(4);
}

Member Data Documentation

_AA

std::vector<double> ossimThinPlateSpline::_AA

private

Definition at line 142 of file ossimThinPlateSpline.h.

Referenced by solve().

_Ainv

std::vector<double> ossimThinPlateSpline::_Ainv

private

Definition at line 142 of file ossimThinPlateSpline.h.

Referenced by solve().

_dx

double ossimThinPlateSpline::_dx

private

Definition at line 130 of file ossimThinPlateSpline.h.

Referenced by getPoint(), and solve().

_dy

double ossimThinPlateSpline::_dy

private

Definition at line 130 of file ossimThinPlateSpline.h.

Referenced by getPoint(), and solve().

_max_nof_points

int ossimThinPlateSpline::_max_nof_points

private

Definition at line 125 of file ossimThinPlateSpline.h.

Referenced by addPoint(), and growPoints().

_nof_eqs

int ossimThinPlateSpline::_nof_eqs

private

Definition at line 126 of file ossimThinPlateSpline.h.

Referenced by solve().

_nof_points

int ossimThinPlateSpline::_nof_points

private

Definition at line 124 of file ossimThinPlateSpline.h.

Referenced by addPoint(), changePoint(), deletePoint(), getPoint(), getXy(), and solve().

_nof_vars

int ossimThinPlateSpline::_nof_vars

private

Definition at line 123 of file ossimThinPlateSpline.h.

Referenced by addPoint(), changePoint(), deletePoint(), getPoint(), and growPoints().

_ta

double ossimThinPlateSpline::_ta

private

Definition at line 129 of file ossimThinPlateSpline.h.

_tx

double ossimThinPlateSpline::_tx

private

Definition at line 128 of file ossimThinPlateSpline.h.

Referenced by deletePoint().

_ty

double ossimThinPlateSpline::_ty

private

Definition at line 128 of file ossimThinPlateSpline.h.

Referenced by deletePoint().

coef

std::vector<std::vector<double> > ossimThinPlateSpline::coef

private

Definition at line 136 of file ossimThinPlateSpline.h.

Referenced by getPoint(), and growPoints().

index

std::vector<int> ossimThinPlateSpline::index

private

Definition at line 140 of file ossimThinPlateSpline.h.

Referenced by changePoint(), getPoint(), getXy(), growPoints(), and solve().

rhs

std::vector<std::vector<double> > ossimThinPlateSpline::rhs

private

Definition at line 135 of file ossimThinPlateSpline.h.

Referenced by addPoint(), changePoint(), deletePoint(), getPoint(), and growPoints().

type

vizGeorefInterType ossimThinPlateSpline::type

private

Definition at line 121 of file ossimThinPlateSpline.h.

Referenced by addPoint(), deletePoint(), getPoint(), and solve().

u

std::vector<double> ossimThinPlateSpline::u

private

Definition at line 138 of file ossimThinPlateSpline.h.

Referenced by getPoint(), growPoints(), and solve().

unused

std::vector<int> ossimThinPlateSpline::unused

private

Definition at line 139 of file ossimThinPlateSpline.h.

Referenced by growPoints(), and solve().

x

std::vector<double> ossimThinPlateSpline::x

private

Definition at line 132 of file ossimThinPlateSpline.h.

Referenced by addPoint(), changePoint(), deletePoint(), getPoint(), getXy(), growPoints(), and solve().

y

std::vector<double> ossimThinPlateSpline::y

private

Definition at line 133 of file ossimThinPlateSpline.h.

Referenced by addPoint(), changePoint(), deletePoint(), getPoint(), getXy(), growPoints(), and solve().

---

The documentation for this class was generated from the following files:

• ossimThinPlateSpline.h
• ossimThinPlateSpline.cpp