PolygonTools.h

#ifndef Polygontools_h
#define Polygontools_h
 
/* SPDX-License-Identifier: CC0-1.0 */
#include <geos/geom/GeometryFactory.h>
#include <geos/geom/MultiPolygon.h>
#include <geos/geom/CoordinateSequence.h>
 
#include <QString>
 
namespace Isis {
 
  class UniversalGroundMap;
  class TProjection;
  
  
  static geos::geom::GeometryFactory::Ptr globalFactory = geos::geom::GeometryFactory::create();
 
  class PolygonTools {
 
    public:
      PolygonTools();
      ~PolygonTools();
 
      static geos::geom::MultiPolygon *LatLonToXY(
        const geos::geom::MultiPolygon &lonLatPoly, TProjection *proj);
 
      static geos::geom::MultiPolygon *XYToLatLon(
        const geos::geom::MultiPolygon &xYPoly, TProjection *proj);
 
      static geos::geom::MultiPolygon *LatLonToSampleLine(
        const geos::geom::MultiPolygon &lonLatPoly, UniversalGroundMap *ugm);
 
      // Return a deep copy of a multpolygon
      static geos::geom::MultiPolygon *CopyMultiPolygon(const geos::geom::MultiPolygon *mpolygon);
      static geos::geom::MultiPolygon *CopyMultiPolygon(const geos::geom::MultiPolygon &mpolygon);
 
      static geos::geom::MultiPolygon *Despike(const geos::geom::Geometry *geom);
      static geos::geom::MultiPolygon *Despike(const geos::geom::MultiPolygon *multiPoly);
      static geos::geom::LinearRing *Despike(const geos::geom::LineString *linearRing);
 
      //  Return polygon in -180/180 coordinated system and merge split polys
      static geos::geom::MultiPolygon *To180(geos::geom::MultiPolygon *poly360);
 
      //Return a polygon in GML format
      static QString ToGML(const geos::geom::MultiPolygon *mpolygon, 
                           QString idString = QString("0"), 
                           QString schema = QString(""));
 
      //Return the GML schema for a polygon
      static QString GMLSchema();
 
      //Return the thickness of a polygon
      static double Thickness(const geos::geom::MultiPolygon *mpolygon);
 
      static geos::geom::Geometry *Intersect(const geos::geom::Geometry *geom1,
                                             const geos::geom::Geometry *geom2);
      static geos::geom::Geometry *Difference(const geos::geom::Geometry *geom1,
                                              const geos::geom::Geometry *geom2);
 
      static geos::geom::MultiPolygon *MakeMultiPolygon(const geos::geom::Geometry *geom);
 
      static QString GetGeometryName(const geos::geom::Geometry *geom);
 
      static bool Equal(const geos::geom::MultiPolygon *poly1,
                        const geos::geom::MultiPolygon *poly2);
      static bool Equal(const geos::geom::Polygon *poly1, const geos::geom::Polygon *poly2);
      static bool Equal(const geos::geom::LineString *lineString1,
                        const geos::geom::LineString *lineString2);
      static bool Equal(const geos::geom::Coordinate &coord1, const geos::geom::Coordinate &coord2);
      static bool Equal(const double d1, const double d2);
 
      static geos::geom::MultiPolygon *FixSeam(const geos::geom::MultiPolygon *poly);
      static geos::geom::MultiPolygon *FixSeam(const geos::geom::Polygon *polyA,
                                               const geos::geom::Polygon *polyB);
 
      static geos::geom::Geometry     *ReducePrecision(const geos::geom::Geometry *geom,
                                                       unsigned int precision);
      static geos::geom::MultiPolygon *ReducePrecision(const geos::geom::MultiPolygon *poly,
                                                       unsigned int precision);
      static geos::geom::Polygon      *ReducePrecision(const geos::geom::Polygon *poly,
                                                       unsigned int precision);
      static geos::geom::LinearRing   *ReducePrecision(const geos::geom::LinearRing *ring,
                                                       unsigned int precision);
      static geos::geom::Coordinate   *ReducePrecision(const geos::geom::Coordinate *coord,
                                                       unsigned int precision);
      static double ReducePrecision(double num, unsigned int precision);
 
      static geos::geom::MultiPolygon *FixPolePolygon(const geos::geom::MultiPolygon *polePolygon,
                                                      UniversalGroundMap *ugm);
      static geos::geom::MultiPolygon *SplitPolygonOn360(const geos::geom::Polygon *inPoly);
 
 
    private:
      static bool IsSpiked(geos::geom::Coordinate first,
                           geos::geom::Coordinate middle, geos::geom::Coordinate last);
      static bool TestSpiked(geos::geom::Coordinate first, geos::geom::Coordinate middle,
                             geos::geom::Coordinate last);
 
      static geos::geom::Geometry     *FixGeometry(const geos::geom::Geometry *geom);
      static geos::geom::MultiPolygon *FixGeometry(const geos::geom::MultiPolygon *poly);
      static geos::geom::Polygon      *FixGeometry(const geos::geom::Polygon *poly);
      static geos::geom::LinearRing   *FixGeometry(const geos::geom::LinearRing *ring);
 
      static geos::geom::Geometry *Operate(const geos::geom::Geometry *geom1,
                                           const geos::geom::Geometry *geom2, unsigned int opcode);
 
      static int DecimalPlace(double);
 
      geos::geom::MultiPolygon *p_polygons;
 
  };
};
 
#endif