geosGeometry.h

Go to the documentation of this file.

#pragma once

#include <features/featureGeometry.h>

#include <geos.h>

#include <tbb/recursive_mutex.h>

namespace MAKVRinTerra
{

class DT_DLL_features DtGeosGeometry : public DtFeatureGeometry::Implementation
{
public:
   typedef boost::shared_ptr<DtGeosGeometry>       Ptr;
   typedef boost::shared_ptr<const DtGeosGeometry> CPtr;

   typedef DtFeatureGeometry::Implementation Implementation;
   typedef DtFeatureGeometry::Type           Type;

   typedef tbb::recursive_mutex Mutex;

   template <Type T>
   struct GeosType
   {
      typedef geos::geom::Geometry type;
   };

   typedef std::vector<geos::geom::Coordinate> CoordinateVector;

   static DtFeatureGeometry MakePoint(const geos::geom::Coordinate&, const DtProj&);
   static DtFeatureGeometry MakeLineString(std::auto_ptr<CoordinateVector>, const DtProj&);
   static DtFeatureGeometry MakeLinearRing(std::auto_ptr<CoordinateVector>, const DtProj&);
   static DtFeatureGeometry MakePolygon(
      std::auto_ptr<CoordinateVector>, std::auto_ptr< std::vector<CoordinateVector> >,
      const DtProj&);
   static DtFeatureGeometry MakeMultiPoint(std::auto_ptr<CoordinateVector>, const DtProj&);

   static DtFeatureGeometry Convert(DtFeatureGeometry);

   static DtFeatureGeometry Make(geos::geom::Geometry*, const DtProj&);

   static DtFeatureGeometry::Area MakeBox(
      DtFeatureGeometry::Point, double length, double width, double rotation);
   static DtFeatureGeometry::Area MakeBox(
      const DtProj&, double xmin, double ymin, double xmax, double ymax);

   static DtFeatureGeometry::Area MakeArea(
      const DtProj&, std::vector<DtPoint> points);
   static DtFeatureGeometry::Area MakeArea(
      const DtProj&, std::vector<DtVector> points);

   static DtFeatureGeometry::Point MakePoint(const DtPoint &, const DtProj &);

   static DtFeatureGeometry::Path MakePath(std::auto_ptr<CoordinateVector>, const DtProj &);

   unsigned numberOfPoints() const;

   DtPoint centroid() const;

   bool intersects(DtFeatureGeometry) const;

   bool overlaps(DtFeatureGeometry) const;

   bool covers(DtFeatureGeometry) const;

   bool within(DtFeatureGeometry) const;

   double distanceTo(DtFeatureGeometry, const DtProj&) const;

   DtFeatureGeometry unionWith(DtFeatureGeometry) const;

   DtFeatureGeometry intersectionWith(DtFeatureGeometry) const;

   DtFeatureGeometry::Point closestPointTo(DtFeatureGeometry, double* distanceTo = NULL) const;



   virtual const geos::geom::Geometry& geosGeometry() const = 0;
};

template <>
struct DtGeosGeometry::GeosType<DtFeatureGeometry::PointType>
{
   typedef geos::geom::Point type;
};

template <>
struct DtGeosGeometry::GeosType<DtFeatureGeometry::PathType>
{
   typedef geos::geom::LineString type;
};

template <>
struct DtGeosGeometry::GeosType<DtFeatureGeometry::AreaType>
{
   typedef geos::geom::Polygon type;
};

template <>
struct DtGeosGeometry::GeosType<DtFeatureGeometry::UnknownType>
{
   typedef geos::geom::GeometryCollection type;
};

template <DtFeatureGeometry::Type T>
class DT_DLL_features DtGeosGeometryImpl : public DtGeosGeometry
{
public:
   typedef boost::shared_ptr<DtGeosGeometryImpl>   Ptr;

   static DtFeatureGeometry Convert(DtFeatureGeometry::Handle<T>);

   Type type() const;

   std::auto_ptr<PointVector> points() const;
   unsigned numberOfSubGeometries() const;
   std::auto_ptr<GeometryVector> subGeometries() const;

//protected:
   const typename GeosType<T>::type& geosGeometry() const = 0;

   virtual boost::shared_ptr<const typename GeosType<T>::type> geosPtr() const = 0;
};

typedef DtGeosGeometryImpl<DtFeatureGeometry::UnknownType>  DtGeosUnknown;
typedef DtGeosGeometryImpl<DtFeatureGeometry::PointType>    DtGeosPoint;
typedef DtGeosGeometryImpl<DtFeatureGeometry::PathType>     DtGeosPath;
typedef DtGeosGeometryImpl<DtFeatureGeometry::AreaType>     DtGeosArea;

}