OverlayMixedPoints.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2020 Paul Ramsey <pramsey@cleverelephant.ca>
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU Lesser General Public Licence as published
 * by the Free Software Foundation.
 * See the COPYING file for more information.
 *
 **********************************************************************/
 
#pragma once
 
 
#include <geos/geom/Point.h>
#include <geos/geom/Polygon.h>
#include <geos/geom/LineString.h>
 
#include <geos/export.h>
 
#include <set>
#include <memory>
#include <vector>
 
 
// Forward declarations
namespace geos {
namespace geom {
class GeometryFactory;
class PrecisionModel;
class Geometry;
class Coordinate;
class CoordinateSequence;
}
namespace algorithm {
namespace locate {
class PointOnGeometryLocator;
}
}
namespace noding {
class Noder;
}
}
 
namespace geos {      // geos.
namespace operation { // geos.operation
namespace overlayng { // geos.operation.overlayng
 
class GEOS_DLL OverlayMixedPoints {
    using Curve = geos::geom::Curve;
    using GeometryFactory = geos::geom::GeometryFactory;
    using PrecisionModel = geos::geom::PrecisionModel;
    using Geometry = geos::geom::Geometry;
    using Coordinate = geos::geom::Coordinate;
    using CoordinateXY = geos::geom::CoordinateXY;
    using CoordinateSequence = geos::geom::CoordinateSequence;
    using Point = geos::geom::Point;
    using Polygon = geos::geom::Polygon;
    using LineString = geos::geom::LineString;
    using Noder = geos::noding::Noder;
    using Surface = geos::geom::Surface;
    using PointOnGeometryLocator = algorithm::locate::PointOnGeometryLocator;
 
private:
 
    // Members
    int opCode;
    const PrecisionModel* pm;
    const Geometry* geomPoint;
    const Geometry* geomNonPointInput;
    const GeometryFactory* geometryFactory;
    bool isPointRHS;
    Noder* noder;
 
    std::unique_ptr<Geometry> geomNonPoint;
    int geomNonPointDim;
    std::unique_ptr<PointOnGeometryLocator> locator;
    int resultDim;
 
    // Methods
    std::unique_ptr<PointOnGeometryLocator> createLocator(const Geometry* geomNonPoint);
 
    std::unique_ptr<Geometry> prepareNonPoint(const Geometry* geomInput);
 
    std::unique_ptr<Geometry> computeIntersection(const CoordinateSequence* coords) const;
 
    std::unique_ptr<Geometry> computeUnion(const CoordinateSequence* coords);
 
    std::unique_ptr<Geometry> computeDifference(const CoordinateSequence* coords);
 
    std::unique_ptr<Geometry> createPointResult(std::vector<std::unique_ptr<Point>>& points) const;
 
    std::vector<std::unique_ptr<Point>> findPoints(bool isCovered, const CoordinateSequence* coords) const;
 
    std::vector<std::unique_ptr<Point>> createPoints(const CoordinateSequence& coords) const;
 
    bool hasLocation(bool isCovered, const CoordinateXY& coord) const;
 
    std::unique_ptr<Geometry> copyNonPoint() const;
 
    static std::unique_ptr<CoordinateSequence> extractCoordinates(const Geometry* points, const PrecisionModel* pm);
 
    static std::vector<std::unique_ptr<Surface>> extractPolygons(const Geometry* geom);
 
    static std::vector<std::unique_ptr<Curve>> extractLines(const Geometry* geom);
 
 
 
public:
 
    OverlayMixedPoints(int p_opCode, const Geometry* geom0, const Geometry* geom1, const PrecisionModel* p_pm, Noder* p_noder);
 
    static std::unique_ptr<Geometry> overlay(int opCode, const Geometry* geom0, const Geometry* geom1, const PrecisionModel* pm, Noder* p_noder);
 
    std::unique_ptr<Geometry> getResult();
 
 
 
 
};
 
 
} // namespace geos.operation.overlayng
} // namespace geos.operation
} // namespace geos