@crossbind/port-geos-ios 2.0.0-beta.50

package/geos.xcframework/ios-arm64/Headers/geos/noding/snapround/SnapRoundingNoder.h

@@ -0,0 +1,167 @@
+/**********************************************************************
+ *
+ * GEOS - Geometry Engine Open Source
+ * http://geos.osgeo.org
+ *
+ * Copyright (C) 2006      Refractions Research Inc.
+ *
+ * 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.
+ *
+ **********************************************************************
+ *
+ * Last port: noding/snapround/SnapRoundingNoder.java r320 (JTS-1.12)
+ *
+ **********************************************************************/
+
+#pragma once
+
+#include <geos/export.h>
+
+#include <geos/geom/Coordinate.h>
+#include <geos/geom/Envelope.h>
+#include <geos/noding/snapround/HotPixelIndex.h>
+#include <geos/noding/Noder.h>
+
+
+// Forward declarations
+namespace geos {
+namespace geom {
+class Envelope;
+class PrecisionModel;
+}
+namespace noding {
+class NodedSegmentString;
+}
+}
+
+namespace geos {
+namespace noding { // geos::noding
+namespace snapround { // geos::noding::snapround
+
+/**
+ * Uses Snap Rounding to compute a rounded,
+ * fully noded arrangement from a set of {@link noding::SegmentString}s,
+ * in a performant way, and avoiding unnecessary noding.
+ *
+ * Implements the Snap Rounding technique described in
+ * the papers by Hobby, Guibas &amp; Marimont, and Goodrich et al.
+ * Snap Rounding enforces that all output vertices lie on a uniform grid,
+ * which is determined by the provided {@link geom::PrecisionModel}.
+ *
+ * Input vertices do not have to be rounded to the grid beforehand;
+ * this is done during the snap-rounding process.
+ * In fact, rounding cannot be done a priori,
+ * since rounding vertices by themselves can distort the rounded topology
+ * of the arrangement (i.e. by moving segments away from hot pixels
+ * that would otherwise intersect them, or by moving vertices
+ * across segments).
+ *
+ * To minimize the number of introduced nodes,
+ * the Snap-Rounding Noder avoids creating nodes
+ * at edge vertices if there is no intersection or snap at that location.
+ * However, if two different input edges contain identical segments,
+ * each of the segment vertices will be noded.
+ * This still provides fully-noded output.
+ * This is the same behaviour provided by other noders,
+ * such as {@link noding::MCIndexNoder} and {@link noding::snap::SnappingNoder}.
+ */
+class GEOS_DLL SnapRoundingNoder : public Noder {
+
+private:
+    /**
+    * The division factor used to determine
+    * nearness distance tolerance for interior intersection detection.
+    */
+    static constexpr int INTERSECTION_NEARNESS_FACTOR = 100;
+
+    // Members
+    const geom::PrecisionModel* pm;
+    noding::snapround::HotPixelIndex pixelIndex;
+    std::vector<SegmentString*> snappedResult;
+
+    // Methods
+    void snapRound(std::vector<SegmentString*>& inputSegStrings, std::vector<SegmentString*>& resultNodedSegments);
+
+    /**
+    * Creates HotPixels for each vertex in the input segStrings.
+    * The HotPixels are not marked as nodes, since they will
+    * only be nodes in the final line arrangement
+    * if they interact with other segments (or they are already
+    * created as intersection nodes).
+    */
+    void addVertexPixels(std::vector<SegmentString*>& segStrings);
+
+    /**
+    * Detects interior intersections in the collection of {@link SegmentString}s,
+    * and adds nodes for them to the segment strings.
+    * Also creates HotPixel nodes for the intersection points.
+    */
+    void addIntersectionPixels(std::vector<SegmentString*>& segStrings);
+
+    /**
+    * Gets a list of the rounded coordinates.
+    * Duplicate (collapsed) coordinates are removed.
+    *
+    * @param pts the coordinates to round
+    * @return array of rounded coordinates
+    */
+    std::unique_ptr<geom::CoordinateSequence> round(const geom::CoordinateSequence& pts) const;
+
+    /**
+    * Computes new segment strings which are rounded and contain
+    * intersections added as a result of snapping segments to snap points (hot pixels).
+    *
+    * @param segStrings segments to snap
+    * @return the snapped segment strings
+    */
+    void computeSnaps(const std::vector<SegmentString*>& segStrings, std::vector<SegmentString*>& snapped);
+    NodedSegmentString* computeSegmentSnaps(NodedSegmentString* ss);
+
+    /**
+    * Snaps a segment in a segmentString to HotPixels that it intersects.
+    *
+    * @param p0 the segment start coordinate
+    * @param p1 the segment end coordinate
+    * @param ss the segment string to add intersections to
+    * @param segIndex the index of the segment
+    */
+    void snapSegment(const geom::CoordinateXY& p0, const geom::CoordinateXY& p1, NodedSegmentString* ss, std::size_t segIndex);
+
+    /**
+    * Add nodes for any vertices in hot pixels that were
+    * added as nodes during segment noding.
+    */
+    void addVertexNodeSnaps(NodedSegmentString* ss);
+
+    void snapVertexNode(const geom::CoordinateXY& p0, NodedSegmentString* ss, std::size_t segIndex);
+
+public:
+
+    SnapRoundingNoder(const geom::PrecisionModel* p_pm)
+        : pm(p_pm)
+        , pixelIndex(p_pm)
+        {}
+
+    /**
+    * @return a Collection of NodedSegmentStrings representing the substrings
+    */
+    std::vector<SegmentString*>* getNodedSubstrings() const override;
+
+    /**
+    * Computes the nodes in the snap-rounding line arrangement.
+    * The nodes are added to the {@link NodedSegmentString}s provided as the input.
+    */
+    void computeNodes(std::vector<SegmentString*>* inputSegStrings) override; //override
+
+};
+
+
+} // namespace geos::noding::snapround
+} // namespace geos::noding
+} // namespace geos
+
+
+

package/geos.xcframework/ios-arm64/Headers/geos/operation/BoundaryOp.h

@@ -0,0 +1,120 @@
+/**********************************************************************
+ *
+ * GEOS - Geometry Engine Open Source
+ * http://geos.osgeo.org
+ *
+ * Copyright (C) 2022 ISciences LLC
+ *
+ * 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.
+ *
+ **********************************************************************
+ *
+ * Last port: operation/BoundaryOp.java fd5aebb
+ *
+ **********************************************************************/
+
+#pragma once
+
+#include <geos/algorithm/BoundaryNodeRule.h>
+#include <geos/geom/Geometry.h>
+#include <map>
+
+namespace geos {
+namespace geom {
+class LineString;
+class MultiLineString;
+}
+}
+
+namespace geos {
+namespace operation {
+
+/**
+ * Computes the boundary of a Geometry.
+ * Allows specifying the BoundaryNodeRule to be used.
+ * This operation will always return a Geometry of the appropriate
+ * dimension for the boundary (even if the input geometry is empty).
+ * The boundary of zero-dimensional geometries (Points) is
+ * always the empty GeometryCollection.
+ *
+ * @author Martin Davis
+ * @version 1.7
+ */
+class GEOS_DLL BoundaryOp {
+
+public:
+  /**
+   * Creates a new instance for the given geometry.
+   *
+   * @param geom the input geometry
+   */
+  BoundaryOp(const geom::Geometry& geom);
+
+  /**
+   * Creates a new instance for the given geometry.
+   *
+   * @param geom the input geometry
+   * @param bnRule the Boundary Node Rule to use
+   */
+  BoundaryOp(const geom::Geometry& geom, const algorithm::BoundaryNodeRule& bnRule);
+
+ /**
+   * Computes a geometry representing the boundary of a geometry.
+   *
+   * @param g the input geometry
+   * @return the computed boundary
+   */
+  static std::unique_ptr<geom::Geometry> getBoundary(const geom::Geometry& g);
+
+  /**
+   * Computes a geometry representing the boundary of a geometry,
+   * using an explicit BoundaryNodeRule.
+   *
+   * @param g the input geometry
+   * @param bnRule the Boundary Node Rule to use
+   * @return the computed boundary
+   */
+  static std::unique_ptr<geom::Geometry> getBoundary(const geom::Geometry& g, const algorithm::BoundaryNodeRule& bnRule);
+
+  /**
+   * Tests if a geometry has a boundary (it is non-empty).
+   * The semantics are:
+   * <ul>
+   * <li>Empty geometries do not have boundaries.
+   * <li>Points do not have boundaries.
+   * <li>For linear geometries the existence of the boundary
+   * is determined by the BoundaryNodeRule.
+   * <li>Non-empty polygons always have a boundary.
+   * </ul>
+   *
+   * @param geom the geometry providing the boundary
+   * @param boundaryNodeRule  the Boundary Node Rule to use
+   * @return true if the boundary exists
+   */
+  static bool hasBoundary(const geom::Geometry& geom, const algorithm::BoundaryNodeRule& boundaryNodeRule);
+
+  /**
+   * Gets the computed boundary.
+   *
+   * @return the boundary geometry
+   */
+  std::unique_ptr<geom::Geometry> getBoundary();
+
+private:
+  const geom::Geometry& m_geom;
+  const geom::GeometryFactory& m_geomFact;
+  const algorithm::BoundaryNodeRule& m_bnRule;
+
+  std::unique_ptr<geom::Geometry> boundaryMultiLineString(const geom::MultiLineString& mLine);
+
+  std::unique_ptr<geom::CoordinateSequence> computeBoundaryCoordinates(const geom::MultiLineString& mLine);
+
+  std::unique_ptr<geom::Geometry> boundaryLineString(const geom::LineString& line);
+};
+
+}
+}
+

package/geos.xcframework/ios-arm64/Headers/geos/operation/GeometryGraphOperation.h

@@ -0,0 +1,92 @@
+/**********************************************************************
+ *
+ * GEOS - Geometry Engine Open Source
+ * http://geos.osgeo.org
+ *
+ * Copyright (C) 2005-2006 Refractions Research Inc.
+ * Copyright (C) 2001-2002 Vivid Solutions Inc.
+ *
+ * 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.
+ *
+ **********************************************************************
+ *
+ * Last port: operation/GeometryGraphOperation.java rev. 1.18 (JTS-1.10)
+ *
+ **********************************************************************/
+
+#pragma once
+
+#include <geos/export.h>
+#include <geos/algorithm/LineIntersector.h> // for composition
+
+#include <vector>
+
+#ifdef _MSC_VER
+#pragma warning(push)
+#pragma warning(disable: 4251) // warning C4251: needs to have dll-interface to be used by clients of class
+#endif
+
+// Forward declarations
+namespace geos {
+namespace algorithm {
+class BoundaryNodeRule;
+}
+namespace geom {
+class Geometry;
+class PrecisionModel;
+}
+namespace geomgraph {
+class GeometryGraph;
+}
+}
+
+
+namespace geos {
+namespace operation { // geos.operation
+
+/// The base class for operations that require GeometryGraph
+class GEOS_DLL GeometryGraphOperation {
+
+public:
+
+    GeometryGraphOperation(const geom::Geometry* g0,
+                           const geom::Geometry* g1);
+
+    GeometryGraphOperation(const geom::Geometry* g0,
+                           const geom::Geometry* g1,
+                           const algorithm::BoundaryNodeRule& boundaryNodeRule);
+
+    GeometryGraphOperation(const geom::Geometry* g0);
+
+    virtual ~GeometryGraphOperation();
+
+    const geom::Geometry* getArgGeometry(unsigned int i) const;
+
+protected:
+
+    algorithm::LineIntersector li;
+
+    const geom::PrecisionModel* resultPrecisionModel;
+
+    /** \brief
+     * The operation args into an array so they can be accessed by index
+     */
+    std::vector<std::unique_ptr<geomgraph::GeometryGraph>> arg;
+
+    void setComputationPrecision(const geom::PrecisionModel* pm);
+
+    // Declare type as noncopyable
+    GeometryGraphOperation(const GeometryGraphOperation& other) = delete;
+    GeometryGraphOperation& operator=(const GeometryGraphOperation& rhs) = delete;
+};
+
+} // namespace geos.operation
+} // namespace geos
+
+#ifdef _MSC_VER
+#pragma warning(pop)
+#endif
+

package/geos.xcframework/ios-arm64/Headers/geos/operation/buffer/BufferBuilder.h

@@ -0,0 +1,269 @@
+/**********************************************************************
+ *
+ * GEOS - Geometry Engine Open Source
+ * http://geos.osgeo.org
+ *
+ * Copyright (C) 2009-2011 Sandro Santilli <strk@kbt.io>
+ * Copyright (C) 2008-2010 Safe Software Inc.
+ * Copyright (C) 2006-2007 Refractions Research Inc.
+ *
+ * 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.
+ *
+ **********************************************************************
+ *
+ * Last port: operation/buffer/BufferBuilder.java 149b38907 (JTS-1.12)
+ *
+ **********************************************************************/
+
+#pragma once
+
+#include <geos/export.h>
+
+#include <vector>
+
+#include <geos/operation/buffer/BufferOp.h> // for inlines (BufferOp enums)
+#include <geos/operation/buffer/OffsetCurveBuilder.h> // for inline (OffsetCurveBuilder enums)
+#include <geos/geomgraph/EdgeList.h> // for composition
+
+#ifdef _MSC_VER
+#pragma warning(push)
+#pragma warning(disable: 4251) // warning C4251: needs to have dll-interface to be used by clients of class
+#endif
+
+// Forward declarations
+namespace geos {
+namespace geom {
+class PrecisionModel;
+class Geometry;
+class GeometryFactory;
+}
+namespace algorithm {
+class LineIntersector;
+}
+namespace noding {
+class Noder;
+class SegmentString;
+class IntersectionAdder;
+}
+namespace geomgraph {
+class Edge;
+class Label;
+class PlanarGraph;
+}
+namespace operation {
+namespace buffer {
+class BufferSubgraph;
+class PolygonBuilder;
+}
+
+}
+}
+
+namespace geos {
+namespace operation { // geos.operation
+namespace buffer { // geos.operation.buffer
+
+/**
+ *
+ * \brief
+ * Builds the buffer geometry for a given input geometry and precision model.
+ *
+ * Allows setting the level of approximation for circular arcs,
+ * and the precision model in which to carry out the computation.
+ *
+ * When computing buffers in floating point double-precision
+ * it can happen that the process of iterated noding can fail to converge
+ * (terminate).
+ *
+ * In this case a TopologyException will be thrown.
+ * Retrying the computation in a fixed precision
+ * can produce more robust results.
+ *
+ */
+class GEOS_DLL BufferBuilder {
+
+public:
+    /**
+     * Creates a new BufferBuilder
+     *
+     * @param nBufParams buffer parameters, this object will
+     *                   keep a reference to the passed parameters
+     *                   so caller must make sure the object is
+     *                   kept alive for the whole lifetime of
+     *                   the buffer builder.
+     */
+    BufferBuilder(const BufferParameters& nBufParams)
+        :
+        bufParams(nBufParams),
+        workingPrecisionModel(nullptr),
+        li(nullptr),
+        intersectionAdder(nullptr),
+        workingNoder(nullptr),
+        geomFact(nullptr),
+        edgeList(),
+        isInvertOrientation(false)
+    {}
+
+    ~BufferBuilder();
+
+
+    /**
+     * Sets the precision model to use during the curve computation
+     * and noding,
+     * if it is different to the precision model of the Geometry.
+     * If the precision model is less than the precision of the
+     * Geometry precision model,
+     * the Geometry must have previously been rounded to that precision.
+     *
+     * @param pm the precision model to use
+     */
+    void
+    setWorkingPrecisionModel(const geom::PrecisionModel* pm)
+    {
+        workingPrecisionModel = pm;
+    }
+
+    /**
+     * Sets the {@link noding::Noder} to use during noding.
+     * This allows choosing fast but non-robust noding, or slower
+     * but robust noding.
+     *
+     * @param newNoder the noder to use
+     */
+    void
+    setNoder(noding::Noder* newNoder)
+    {
+        workingNoder = newNoder;
+    }
+
+    /**
+    * Sets whether the offset curve is generated
+    * using the inverted orientation of input rings.
+    * This allows generating a buffer(0) polygon from the smaller lobes
+    * of self-crossing rings.
+    *
+    * @param p_isInvertOrientation true if input ring orientation should be inverted
+    */
+    void
+    setInvertOrientation(bool p_isInvertOrientation)
+    {
+        isInvertOrientation = p_isInvertOrientation;
+    }
+
+
+    std::unique_ptr<geom::Geometry> buffer(const geom::Geometry* g, double distance);
+
+    /**
+     * Generates offset curve for linear geometry.
+     *
+     * @param g non-areal geometry object
+     * @param distance width of offset
+     * @param leftSide controls on which side of the input geometry
+     *        offset curve is generated.
+     *
+     * @note For left-side offset curve, the offset will be at the left side
+     *       of the input line and retain the same direction.
+     *       For right-side offset curve, it'll be at the right side
+     *       and in the opposite direction.
+     *
+     * @note BufferParameters::setSingleSided parameter, which is specific to
+     *       areal geometries only, is ignored by this routine.
+     *
+     * @note Not in JTS: this is a GEOS extension
+     */
+    std::unique_ptr<geom::Geometry> bufferLineSingleSided(
+        const geom::Geometry* g,
+        double distance, bool leftSide);
+
+private:
+
+    /**
+     * Compute the change in depth as an edge is crossed from R to L
+     */
+    static int depthDelta(const geomgraph::Label& label);
+
+    const BufferParameters& bufParams;
+
+    const geom::PrecisionModel* workingPrecisionModel;
+
+    algorithm::LineIntersector* li;
+
+    noding::IntersectionAdder* intersectionAdder;
+
+    noding::Noder* workingNoder;
+
+    const geom::GeometryFactory* geomFact;
+
+    geomgraph::EdgeList edgeList;
+
+    std::vector<geomgraph::Label*> newLabels;
+
+    bool isInvertOrientation;
+
+    void computeNodedEdges(std::vector<noding::SegmentString*>& bufSegStr,
+                           const geom::PrecisionModel* precisionModel);
+    // throw(GEOSException);
+
+    /**
+     * Inserted edges are checked to see if an identical edge already
+     * exists.
+     * If so, the edge is not inserted, but its label is merged
+     * with the existing edge.
+     *
+     * The function takes responsibility of releasing the Edge parameter
+     * memory when appropriate.
+     */
+    void insertUniqueEdge(geomgraph::Edge* e);
+
+    void createSubgraphs(geomgraph::PlanarGraph* graph,
+                         std::vector<BufferSubgraph*>& list);
+
+    /**
+     * Completes the building of the input subgraphs by
+     * depth-labelling them,
+     * and adds them to the PolygonBuilder.
+     * The subgraph list must be sorted in rightmost-coordinate order.
+     *
+     * @param subgraphList the subgraphs to build
+     * @param polyBuilder the PolygonBuilder which will build
+     *        the final polygons
+     */
+    void buildSubgraphs(const std::vector<BufferSubgraph*>& subgraphList,
+                        PolygonBuilder& polyBuilder);
+
+    /// \brief
+    /// Return the externally-set noding::Noder OR a newly created
+    /// one using the given precisionModel.
+    ///
+    /// NOTE: if an externally-set noding::Noder is available no
+    /// check is performed to ensure it will use the
+    /// given PrecisionModel
+    ///
+    noding::Noder* getNoder(const geom::PrecisionModel* precisionModel);
+
+
+    /**
+     * Gets the standard result for an empty buffer.
+     * Since buffer always returns a polygonal result,
+     * this is chosen to be an empty polygon.
+     *
+     * @return the empty result geometry, transferring ownership to caller.
+     */
+    std::unique_ptr<geom::Geometry> createEmptyResultGeometry() const;
+
+    // Declare type as noncopyable
+    BufferBuilder(const BufferBuilder& other) = delete;
+    BufferBuilder& operator=(const BufferBuilder& rhs) = delete;
+};
+
+} // namespace geos::operation::buffer
+} // namespace geos::operation
+} // namespace geos
+
+#ifdef _MSC_VER
+#pragma warning(pop)
+#endif
+