Utils.h

Go to the documentation of this file.

 
/*!
\file terralib/mnt/core/Utils.h
 
\brief Utility functions for MNT support.
*/
 
#ifndef __TERRALIB_MNT_INTERNAL_UTILS_H
#define __TERRALIB_MNT_INTERNAL_UTILS_H
 
#include "Config.h"
#include "Enums.h"
 
#include "../../common/UnitsOfMeasureManager.h"
 
#include "../../dataaccess/dataset/DataSet.h"
#include "../../dataaccess/dataset/DataSetType.h"
#include "../../dataaccess/datasource/DataSource.h"
 
#include "../../geometry/Coord2D.h"
#include "../../geometry/LineString.h"
#include "../../geometry/MultiLineString.h"
#include "../../geometry/MultiPoint.h"
#include "../../geometry/Point.h"
 
 
#ifndef MIN
#define MIN(x, y) ( ((x) < (y)) ? (x) : (y) )
#endif
 
#ifndef MAX
#define MAX(x, y) ( ((x) > (y)) ? (x) : (y) )
#endif
 
 
namespace te
{
  namespace mnt
  {
    TEMNTEXPORT size_t ReadPoints(  std::string &inDsetName, te::da::DataSourcePtr &inDsrc, std::string &atrZ, double tol, te::gm::MultiPoint &mpt, std::string &geostype, 
                        te::gm::Envelope &env);
 
  TEMNTEXPORT size_t ReadSamples(std::string &inDsetName, te::da::DataSourcePtr &inDsrc, std::string &atrZ, double tol, double max, Simplify alg, te::gm::MultiPoint &mpt,
                        te::gm::MultiLineString &isolines, std::string &geostype, te::gm::Envelope &env, int srid = 0);
 
  TEMNTEXPORT double Distance(const te::gm::Coord2D &pt1, const te::gm::Coord2D &pt2);
 
  TEMNTEXPORT double Distance(const double pt1x, const double pt1y, const double pt2x, const double pt2y);
 
  TEMNTEXPORT bool Equal(te::gm::Point &p1, te::gm::Point &p2, double &tol);
 
    //DouglasPeucker from TerraAmazon 4
  TEMNTEXPORT te::gm::LineString* DouglasPeuckerTA(te::gm::LineString *lineIn, double simpFactor, double Zvalue);
 
    //DouglasPeucker from libgeos
  TEMNTEXPORT te::gm::LineString* GEOS_DouglasPeucker(te::gm::LineString *ls, double snap, double Zvalue);
 
    TEMNTEXPORT te::gm::LineString* pointListSimplify(te::gm::LineString *ls, double snap, double maxdist, double Zvalue);
 
  TEMNTEXPORT double SegmentDistance(double fx, double fy, double lx, double ly, double ptx, double pty, double *pix, double *piy);
      
  TEMNTEXPORT double pointToSegmentDistance(te::gm::Point &fseg, te::gm::Point &lseg, te::gm::Point &pt, te::gm::Point *pti);
 
  TEMNTEXPORT double coordToSegmentDistance(te::gm::Coord2D &fseg, te::gm::Coord2D &lseg, te::gm::Coord2D &pt, te::gm::Coord2D *pti);
 
    // Find center point of triangle using its vertices.
  TEMNTEXPORT short findCenter(te::gm::Point* vert, double* pcx, double* pcy);
 
  TEMNTEXPORT short segIntersect(te::gm::Point &pfr, te::gm::Point &pto, te::gm::Point &lfr, te::gm::Point &lto);
  TEMNTEXPORT bool segInterPoint(te::gm::Point &pfr, te::gm::Point &pto, te::gm::Point &lfr, te::gm::Point &lto, te::gm::Point *pt);
 
    // Test Vertex values for valid and special cases, returns TRUE if normal, FALSE otherwise.
  TEMNTEXPORT bool testVertexValues(double, te::gm::Point*);
 
    // Define intersections of contour with edges of triangle
  TEMNTEXPORT bool defineInterTriangle(std::vector<te::gm::Point>&, te::gm::Point*, double);
 
    // Define intersection point of contour with edge between two points
  TEMNTEXPORT bool defineInterEdge(std::vector<te::gm::Point>&, te::gm::Point&, te::gm::Point&, double);
 
    //Extract lines from pline and store connected ones in clinlist.
  TEMNTEXPORT bool extractLines(std::vector<te::gm::Point> &pline, std::vector< std::unique_ptr< te::gm::LineString > > &clinlist, double scale);
 
  TEMNTEXPORT bool initLineVector(std::vector<te::gm::Point> &pline, std::vector<te::gm::Point> &vect);
  TEMNTEXPORT bool assembLine(std::vector< std::unique_ptr< te::gm::LineString > > &linlout, std::vector<te::gm::Point> &vect);
 
  TEMNTEXPORT bool equalFptSpt(te::gm::Point & fpt, te::gm::Point &spt, double scale);
 
  TEMNTEXPORT te::gm::Point Max(te::gm::Point& p1, te::gm::Point& p2);
  TEMNTEXPORT te::gm::Point Min(te::gm::Point& p1, te::gm::Point& p2);
 
    // Calculate Normal vector of Triangle
  TEMNTEXPORT bool triangleNormalVector(te::gm::Point*, double*);
 
  TEMNTEXPORT bool normalizeVector(double*);
 
  TEMNTEXPORT bool onSegment(te::gm::Point& pt, te::gm::Point& fseg, te::gm::Point& lseg, double tol);
 
  TEMNTEXPORT int onSameSide(te::gm::Coord2D pt1, te::gm::Coord2D pt2, te::gm::Coord2D fseg, te::gm::Coord2D lseg);
 
  TEMNTEXPORT bool point3dListFilter(std::vector<te::gm::Point> &p3dl, std::vector<bool> &fixed, double tol);
 
  TEMNTEXPORT bool Least_square_fitting(double *vectx, double *vecty, short np, short deg, double *coef);
 
  TEMNTEXPORT bool Gauss_elimination(short m, short n, double mat[6][6]);
 
  TEMNTEXPORT bool SaveIso(std::string& outDsetName, te::da::DataSourcePtr &outDsrc, std::vector< std::unique_ptr< te::gm::LineString > > &isolist, std::vector<double> &guidevalues, int srid);
 
    TEMNTEXPORT bool convertAngleToPlanar(double& val, te::common::UnitOfMeasurePtr planar, int type);
 
    TEMNTEXPORT bool convertPlanarToAngle(double& val, te::common::UnitOfMeasurePtr unit);
 
    TEMNTEXPORT void getMinMax(te::rst::Raster* inputRst, double &vmin, double &vmax);
 
  TEMNTEXPORT double PerpendicularDistance(te::gm::Coord2D& first, te::gm::Coord2D& last, te::gm::Coord2D& pin, te::gm::Coord2D& pinter);
 
    TEMNTEXPORT te::mnt::mntType getMNTType(const te::da::DataSetType* dt);
  }
}
 
#endif