SplineGrassMitasova.h

Go to the documentation of this file.

 
/*!
\file terralib/mnt/core/SplineGrassMitasova.h
 
\brief This file contains a class to calculate retangular grid from Samples using Spline Grass Mitasova Interpolation.
Adapted from SPRING
 
*/
 
#ifndef __TERRALIB_MNT_INTERNAL_SPLINEGRASSMITASOVA_H
#define __TERRALIB_MNT_INTERNAL_SPLINEGRASSMITASOVA_H
 
#include "CalculateGrid.h"
 
namespace te
{
  namespace mnt
  {
 
    class TEMNTEXPORT SplineInterpolationGrassMitasova : public CalculateGrid
    {
    public:
      //! Contructors
      SplineInterpolationGrassMitasova(int minPoints = mitasovaMINPOINTS, double tension = mitasovaTENSION, double smooth = mitasovaSMOOTH);
 
      //! Destructor
      virtual ~SplineInterpolationGrassMitasova();
 
      void setControlPoints(double xMin, double xMax, double yMin, double yMax, int nMax, double xOrigem, double yOrigem, double dnorm);
 
      bool calculateGrid();
 
    private:
      int m_minimoPontos;
      double m_tension;
      double m_smooth;
 
      int G_ludcmp(std::vector< std::vector<double> > &, int, std::vector<int> &, double &);
      void G_lubksb(std::vector< std::vector<double> > &a, int n, std::vector<int> &indx, std::vector<double> &b);
      double amax1(double arg1, double arg2);
      double amin1(double arg1, double arg2);
      double IL_crst(double r, double fi);
      int IL_matrix_create(std::vector< std::vector<double> > &matrix, std::vector<int> &indx, double fi, double rsm, int KMAX2);
      bool IL_grid_calc(double xMin, double xMax, double yMin, double yMax, double zMin, double zMax, int KMAX2, std::vector<double>& b, double fi, double dnorm);
      int IL_check_at_points_2d(std::vector<double> &b, double &ertot, double zmin, double dnorm, double fi);
 
      std::vector<te::gm::Point> m_pointsRegion;  //!< Control Points of grid region
 
    }; //class SplineInterpolationGrassMitasova
  } //mnt
} //te
 
#endif