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/Point.h>
#include <geos/geom/LineString.h>
#include <geos/geom/Polygon.h>
#include <geos/geom/GeometryCollection.h>
#include <geos/geom/prep/PreparedGeometry.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 std::shared_ptr<DtGeosGeometry>       Ptr;
   typedef std::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::unique_ptr<CoordinateVector>, const DtProj&);
   static DtFeatureGeometry MakeLinearRing(std::unique_ptr<CoordinateVector>, const DtProj&);
   static DtFeatureGeometry MakePolygon(
      std::unique_ptr<CoordinateVector>, std::unique_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::unique_ptr<CoordinateVector>, const DtProj &);

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

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

   virtual std::unique_ptr<PointVector> points() const override;

   virtual unsigned numberOfPoints() const override;

   virtual bool isValid() const override;
   virtual bool isEmpty() const override;

   virtual unsigned dimension() const override;

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

   virtual std::unique_ptr<DtFeatureGeometry::RelationshipMatrix> relate(DtFeatureGeometry) const override;

   virtual double length() const override;

   virtual double area() const override;

   virtual DtFeatureGeometry::Point centroid() const override;

   virtual double distanceTo(DtFeatureGeometry, const DtProj&) const override;

   virtual DtFeatureGeometry unionWith(DtFeatureGeometry) const override;

   virtual DtFeatureGeometry intersectionWith(DtFeatureGeometry) const override;

   virtual DtFeatureGeometry differenceWith(DtFeatureGeometry) const override;

   virtual DtFeatureGeometry buffer(double width, unsigned segments) const override;

   virtual DtFeatureGeometry bufferWithFlatCap(double width, unsigned segments) const override;

   virtual DtFeatureGeometry convexHull() const override;

   virtual DtFeatureGeometry boundary() const override;

   virtual DtFeatureGeometry simplify(double tolerance) const override;

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

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

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

   std::shared_ptr<DtGeosGeometry> shared_from_this();
   std::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 std::shared_ptr<DtGeosGeometryImpl>   Ptr;

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

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

   virtual Type type() const override;

   virtual unsigned numberOfSubGeometries() const override;
   virtual std::unique_ptr<GeometryVector> subGeometries() const override;

   virtual const typename GeosType<T>::type& geosGeometry() const override
   {
      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);