@crossbind/port-geos-linuxmusl 0.0.0-stage → 2.0.0-beta.63

package/dist/prebuilt/linuxmusl-x64-mt-release/include/geos/noding/snapround/HotPixelIndex.h

@@ -0,0 +1,116 @@
+/**********************************************************************
+ *
+ * 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/export.h>
+#include <geos/geom/Envelope.h> // for unique_ptr
+#include <geos/geom/Coordinate.h> // for composition
+#include <geos/algorithm/LineIntersector.h>
+#include <geos/noding/snapround/HotPixel.h>
+#include <geos/geom/PrecisionModel.h>
+#include <geos/util/IllegalArgumentException.h>
+#include <geos/io/WKTWriter.h>
+#include <geos/index/kdtree/KdTree.h>
+#include <geos/index/kdtree/KdNodeVisitor.h>
+#include <geos/util.h>
+
+#include <array>
+#include <map>
+#include <memory>
+
+
+#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 LineIntersector;
+}
+namespace index {
+class ItemVisitor;
+}
+namespace noding {
+namespace snapround {
+class HotPixel;
+}
+}
+}
+
+
+namespace geos {
+namespace noding { // geos::noding
+namespace snapround { // geos::noding::snapround
+
+
+class GEOS_DLL HotPixelIndex {
+
+private:
+
+    /* members */
+    const geom::PrecisionModel* pm;
+    double scaleFactor;
+    std::unique_ptr<geos::index::kdtree::KdTree> index;
+    std::deque<HotPixel> hotPixelQue;
+
+    /* methods */
+    template<typename CoordType>
+    geom::CoordinateXYZM round(const CoordType& pt) {
+        geom::CoordinateXYZM p2(pt);
+        pm->makePrecise(p2);
+        return p2;
+    }
+
+    HotPixel* find(const geom::Coordinate& pixelPt);
+
+public:
+
+    HotPixelIndex(const geom::PrecisionModel* p_pm);
+    HotPixel* addRounded(const geom::CoordinateXYZM& pt);
+
+    template<typename CoordType>
+    HotPixel* add(const CoordType& p) {
+        static_assert(std::is_base_of<geom::CoordinateXY, CoordType>(), "Only valid for Coordinate types");
+
+        auto pRound = round(p);
+        return addRounded(pRound);
+    }
+
+    void add(const geom::CoordinateSequence* pts);
+    void add(const std::vector<geom::Coordinate>& pts);
+    void addNodes(const geom::CoordinateSequence* pts);
+    void addNodes(const std::vector<geom::Coordinate>& pts);
+
+    /**
+    * Visits all the hot pixels which may intersect a segment (p0-p1).
+    * The visitor must determine whether each hot pixel actually intersects
+    * the segment.
+    */
+    void query(const geom::CoordinateXY& p0, const geom::CoordinateXY& p1,
+               index::kdtree::KdNodeVisitor& visitor);
+
+};
+
+} // namespace geos::noding::snapround
+} // namespace geos::noding
+} // namespace geos
+
+#ifdef _MSC_VER
+#pragma warning(pop)
+#endif
+
+

package/dist/prebuilt/linuxmusl-x64-mt-release/include/geos/noding/snapround/MCIndexPointSnapper.h

@@ -0,0 +1,98 @@
+/**********************************************************************
+ *
+ * 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/MCIndexPointSnapper.java r486 (JTS-1.12+)
+ *
+ **********************************************************************/
+
+#pragma once
+
+#include <cstddef>
+#include <geos/export.h>
+
+// Forward declarations
+namespace geos {
+namespace geom {
+class Envelope;
+}
+namespace index {
+class SpatialIndex;
+}
+namespace noding {
+class SegmentString;
+namespace snapround {
+class HotPixel;
+}
+}
+}
+
+namespace geos {
+namespace noding { // geos::noding
+namespace snapround { // geos::noding::snapround
+
+/** \brief
+ * "Snaps" all [SegmentStrings](@ref SegmentString) in a [SpatialIndex](@ref index::SpatialIndex) containing
+ * [MonotoneChains](@ref index::chain::MonotoneChain) to a given HotPixel.
+ *
+ */
+class GEOS_DLL MCIndexPointSnapper {
+
+public:
+
+
+    MCIndexPointSnapper(index::SpatialIndex& nIndex)
+        :
+        index(nIndex)
+    {}
+
+    /** \brief
+     * Snaps (nodes) all interacting segments to this hot pixel.
+     *
+     * The hot pixel may represent a vertex of an edge, in which case this
+     * routine uses the optimization of not noding the vertex itself.
+     *
+     * @param hotPixel the hot pixel to snap to
+     * @param parentEdge the edge containing the vertex,
+     *                   if applicable, or `null`
+     * @param vertexIndex the index of the vertex, if applicable, or -1
+     * @return `true` if a node was added for this pixel
+     */
+    bool snap(HotPixel& hotPixel, SegmentString* parentEdge,
+              std::size_t vertexIndex);
+
+    bool
+    snap(HotPixel& hotPixel)
+    {
+        return snap(hotPixel, nullptr, 0);
+    }
+
+    geom::Envelope getSafeEnvelope(const HotPixel& hp) const;
+
+
+private:
+
+    static constexpr double SAFE_ENV_EXPANSION_FACTOR = 0.75;
+
+    index::SpatialIndex& index;
+
+    // Declare type as noncopyable
+    MCIndexPointSnapper(const MCIndexPointSnapper& other) = delete;
+    MCIndexPointSnapper& operator=(const MCIndexPointSnapper& rhs) = delete;
+};
+
+
+} // namespace geos::noding::snapround
+} // namespace geos::noding
+} // namespace geos
+

package/dist/prebuilt/linuxmusl-x64-mt-release/include/geos/noding/snapround/MCIndexSnapRounder.h

@@ -0,0 +1,157 @@
+/**********************************************************************
+ *
+ * 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/MCIndexSnapRounder.java r486 (JTS-1.12+)
+ *
+ **********************************************************************/
+
+#pragma once
+
+#include <geos/export.h>
+
+#include <geos/noding/Noder.h> // for inheritance
+#include <geos/noding/NodedSegmentString.h> // for inlines
+#include <geos/noding/snapround/MCIndexPointSnapper.h> // for inlines
+#include <geos/algorithm/LineIntersector.h> // for composition
+#include <geos/geom/Coordinate.h> // for use in vector
+#include <geos/geom/PrecisionModel.h> // for inlines
+
+#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 LineIntersector;
+}
+namespace noding {
+class SegmentString;
+class MCIndexNoder;
+}
+}
+
+namespace geos {
+namespace noding { // geos::noding
+namespace snapround { // geos::noding::snapround
+
+
+/** \brief
+ * Uses Snap Rounding to compute a rounded,
+ * fully noded arrangement from a set of SegmentString
+ *
+ * Implements the Snap Rounding technique described in Hobby, Guibas & Marimont,
+ * and Goodrich et al.
+ *
+ * Snap Rounding assumes that all vertices lie on a uniform grid
+ * (hence the precision model of the input must be fixed precision,
+ * and all the input vertices must be rounded to that precision).
+ *
+ * This implementation uses a monotone chains and a spatial index to
+ * speed up the intersection tests.
+ *
+ * This implementation appears to be fully robust using an integer
+ * precision model.
+ *
+ * It will function with non-integer precision models, but the
+ * results are not 100% guaranteed to be correctly noded.
+ */
+class GEOS_DLL MCIndexSnapRounder: public Noder { // implements Noder
+
+public:
+
+    MCIndexSnapRounder(const geom::PrecisionModel& nPm)
+        :
+        pm(nPm),
+        scaleFactor(nPm.getScale()),
+        pointSnapper(nullptr)
+    {
+        li.setPrecisionModel(&pm);
+    }
+
+    std::vector<std::unique_ptr<SegmentString>>
+    getNodedSubstrings() override
+    {
+        return NodedSegmentString::getNodedSubstrings(nodedSegStrings);
+    }
+
+    void computeNodes(const std::vector<SegmentString*>& segStrings) override;
+
+    /**
+     * Computes nodes introduced as a result of
+     * snapping segments to vertices of other segments
+     *
+     * @param edges the list of segment strings to snap together
+     *        NOTE: they *must* be instances of NodedSegmentString, or
+     * 	            an assertion will fail.
+     */
+    void computeVertexSnaps(const std::vector<SegmentString*>& edges);
+
+private:
+
+    /// externally owned
+    const geom::PrecisionModel& pm;
+
+    algorithm::LineIntersector li;
+
+    double scaleFactor;
+
+    std::vector<SegmentString*> nodedSegStrings;
+
+    std::unique_ptr<MCIndexPointSnapper> pointSnapper;
+
+    void snapRound(MCIndexNoder& noder, const std::vector<SegmentString*>& segStrings);
+
+
+    /**
+     * Computes all interior intersections in the collection of SegmentStrings,
+     * and push their Coordinate to the provided vector.
+     *
+     * Does NOT node the segStrings.
+     *
+     */
+    void findInteriorIntersections(MCIndexNoder& noder,
+                                   const std::vector<SegmentString*>& segStrings,
+                                   std::vector<geom::Coordinate>& intersections);
+
+    /**
+     * Computes nodes introduced as a result of snapping
+     * segments to snap points (hot pixels)
+     */
+    void computeIntersectionSnaps(std::vector<geom::Coordinate>& snapPts);
+
+    /**
+     * Performs a brute-force comparison of every segment in each {@link SegmentString}.
+     * This has n^2 performance.
+     */
+    void computeVertexSnaps(NodedSegmentString* e);
+
+    void checkCorrectness(std::vector<SegmentString*>& inputSegmentStrings);
+
+    // Declare type as noncopyable
+    MCIndexSnapRounder(const MCIndexSnapRounder& other) = delete;
+    MCIndexSnapRounder& operator=(const MCIndexSnapRounder& rhs) = delete;
+};
+
+} // namespace geos::noding::snapround
+} // namespace geos::noding
+} // namespace geos
+
+#ifdef _MSC_VER
+#pragma warning(pop)
+#endif
+

package/dist/prebuilt/linuxmusl-x64-mt-release/include/geos/noding/snapround/SnapRoundingIntersectionAdder.h

@@ -0,0 +1,126 @@
+/**********************************************************************
+ *
+ * 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/export.h>
+
+#include <vector>
+#include <memory>
+
+#include <geos/noding/Noder.h> // for inheritance
+#include <geos/algorithm/LineIntersector.h> // for composition
+#include <geos/geom/Coordinate.h> // for use in vector
+#include <geos/geom/CoordinateSequence.h>
+#include <geos/geom/PrecisionModel.h> // for inlines (should drop)
+#include <geos/noding/SegmentIntersector.h>
+
+
+// Forward declarations
+namespace geos {
+namespace geom {
+class PrecisionModel;
+}
+namespace noding {
+class SegmentString;
+class NodedSegmentString;
+namespace snapround {
+class HotPixel;
+}
+}
+}
+
+namespace geos {
+namespace noding { // geos::noding
+namespace snapround { // geos::noding::snapround
+
+/**
+ * Finds intersections between line segments which will be snap-rounded,
+ * and adds them as nodes to the segments.
+ *
+ * Intersections are detected and computed using full precision.
+ * Snapping takes place in a subsequent phase.
+ *
+ * The intersection points are recorded, so that HotPixels can be created for them.
+ *
+ * To avoid robustness issues with vertices which lie very close to line segments
+ * a heuristic is used:
+ * nodes are created if a vertex lies within a tolerance distance
+ * of the interior of a segment.
+ * The tolerance distance is chosen to be significantly below the snap-rounding grid size.
+ * This has empirically proven to eliminate noding failures.
+ */
+class GEOS_DLL SnapRoundingIntersectionAdder: public SegmentIntersector { // implements SegmentIntersector
+
+private:
+
+    algorithm::LineIntersector li;
+    geom::CoordinateSequence intersections;
+    // const geom::PrecisionModel* pm;
+    double nearnessTol;
+
+    /**
+    * If an endpoint of one segment is near
+    * the interior of the other segment, add it as an intersection.
+    * EXCEPT if the endpoint is also close to a segment endpoint
+    * (since this can introduce "zigs" in the linework).
+    *
+    * This resolves situations where
+    * a segment A endpoint is extremely close to another segment B,
+    * but is not quite crossing.  Due to robustness issues
+    * in orientation detection, this can
+    * result in the snapped segment A crossing segment B
+    * without a node being introduced.
+    */
+    void processNearVertex(const geom::CoordinateSequence& seq0,
+                           std::size_t ptIndex,
+                           const geom::CoordinateSequence& seq1,
+                           std::size_t segIndex,
+                           SegmentString* edge);
+
+    bool isNearSegmentInterior(const geom::CoordinateXY& p, const geom::CoordinateXY& p0, const geom::CoordinateXY& p1) const;
+
+public:
+
+    SnapRoundingIntersectionAdder(double p_nearnessTol)
+        : SegmentIntersector()
+        , intersections(geom::CoordinateSequence::XYZM(0))
+        , nearnessTol(p_nearnessTol)
+    {}
+
+    geom::CoordinateSequence getIntersections() { return std::move(intersections); };
+
+    /**
+    * This method is called by clients
+    * of the {@link SegmentIntersector} class to process
+    * intersections for two segments of the {@link SegmentString}s being intersected.
+    * Note that some clients (such as MonotoneChains) may optimize away
+    * this call for segment pairs which they have determined do not intersect
+    * (e.g. by an disjoint envelope test).
+    */
+    void processIntersections(SegmentString* e0, std::size_t segIndex0, SegmentString* e1, std::size_t segIndex1) override;
+
+    /**
+    * Always process all intersections
+    *
+    */
+    bool isDone() const override { return false; }
+
+
+};
+
+} // namespace geos::noding::snapround
+} // namespace geos::noding
+} // namespace geos
+

package/dist/prebuilt/linuxmusl-x64-mt-release/include/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(const 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(const 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(const 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<std::unique_ptr<SegmentString>> getNodedSubstrings() 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(const std::vector<SegmentString*>& inputSegStrings) override; //override
+
+};
+
+
+} // namespace geos::noding::snapround
+} // namespace geos::noding
+} // namespace geos
+
+
+