PolygonTools.h

Go to the documentation of this file.

#ifndef Polygontools_h
#define Polygontools_h

#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 globalFactory;

  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