geosGeometry.h

Go to the documentation of this file.

/*******************************************************************************
** Copyright (c) 2020 MAK Technologies, Inc.
** All rights reserved.
*******************************************************************************/


#pragma once

#include <features/featureGeometry.h>

#include <geos.h>
#include <geos/geom/prep/PreparedGeometryFactory.h>

#include <boost/optional.hpp>

#include <vrfutil/objectCounter.h>

namespace MAKVRinTerra
{

class DT_DLL_features DtGeosGeometry
   : public DtFeatureGeometry::Implementation
#ifdef DtObjectDebugger_Enable
   , public DtObjectDebugger<DtGeosGeometry>
#endif
{
public:
   typedef boost::shared_ptr<DtGeosGeometry>       Ptr;
   typedef boost::shared_ptr<const DtGeosGeometry> CPtr;

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

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

   DtGeosGeometry() {};
   DtGeosGeometry(const DtGeosGeometry& orig);

   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 DtGeosGeometry::Ptr Convert(DtFeatureGeometry);

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

   static DtFeatureGeometry MakeSet(const GeometryVector&, 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 MakeBox(const DtProj::Bounds&);

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

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

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

   static std::vector<DtFeatureGeometry::Path> MergeLines(DtFeatureGeometry);

   std::string description() const;
   std::string toWKT() const;

   std::auto_ptr<PointVector> points() const;

   unsigned numberOfPoints() const;

   bool isValid() const;
   bool isEmpty() const;

   unsigned dimension() const;

   bool equals(DtFeatureGeometry) const;
   bool intersects(DtFeatureGeometry) const;
   bool overlaps(DtFeatureGeometry) const;
   bool crosses(DtFeatureGeometry) const;
   bool covers(DtFeatureGeometry) const;
   bool within(DtFeatureGeometry) const;
   bool touches(DtFeatureGeometry) const;
   bool contains(DtFeatureGeometry) const;

   std::auto_ptr<DtFeatureGeometry::RelationshipMatrix> relate(DtFeatureGeometry) const;

   double length() const;

   double area() const;

   DtFeatureGeometry::Point centroid() const;

   double distanceTo(DtFeatureGeometry, const DtProj&) const;

   DtFeatureGeometry unionWith(DtFeatureGeometry) const;

   DtFeatureGeometry intersectionWith(DtFeatureGeometry) const;

   DtFeatureGeometry differenceWith(DtFeatureGeometry) const;

   DtFeatureGeometry buffer(double width, unsigned segments) const;

   DtFeatureGeometry bufferWithFlatCap(double width, unsigned segments) const;

   DtFeatureGeometry convexHull() const;

   DtFeatureGeometry boundary() const;

   DtFeatureGeometry simplify(double tolerance) const;

   DtFeatureGeometry::Point closestPointTo(
      DtFeatureGeometry, const DtProj&, double* distanceTo = 0) const;

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

   virtual const geos::geom::prep::PreparedGeometry* geosPreparedGeometry() const = 0;

   boost::shared_ptr<DtGeosGeometry> shared_from_this();
   boost::shared_ptr<const DtGeosGeometry> shared_from_this() const;

   void markValid() { myValidFlag = true; };

protected:
   boost::optional<bool> myValidFlag;
};

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
#ifdef DtObjectDebugger_Enable
   , public DtObjectDebugger< DtGeosGeometryImpl<T> >
#endif
{
public:

   typedef boost::shared_ptr<DtGeosGeometryImpl>   Ptr;

   DtGeosGeometryImpl() {};
   DtGeosGeometryImpl(const DtGeosGeometryImpl& orig) : DtGeosGeometry(orig) {};

   static DtGeosGeometry::Ptr Convert(const DtFeatureGeometry::Handle<T>&);

   Type type() const;

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

   const typename GeosType<T>::type& geosGeometry() const
   {
      DtTHROW_NEW(DtFeatureGeometry::InvalidGeometry, "Invalid geosGeometry call");
   }
};

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

}

DT_DEFINE_OBJECT_DEBUGGER_NAME(MAKVRinTerra::DtGeosGeometry);
DT_DEFINE_OBJECT_DEBUGGER_NAME(MAKVRinTerra::DtGeosUnknown);
DT_DEFINE_OBJECT_DEBUGGER_NAME(MAKVRinTerra::DtGeosPoint);
DT_DEFINE_OBJECT_DEBUGGER_NAME(MAKVRinTerra::DtGeosPath);
DT_DEFINE_OBJECT_DEBUGGER_NAME(MAKVRinTerra::DtGeosArea);