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

package/dist/prebuilt/linuxmusl-x64-mt-release/include/geos/operation/buffer/BufferOp.h

@@ -0,0 +1,307 @@
+/**********************************************************************
+ *
+ * GEOS - Geometry Engine Open Source
+ * http://geos.osgeo.org
+ *
+ * Copyright (C) 2009-2011 Sandro Santilli <strk@kbt.io>
+ * Copyright (C) 2005-2007 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/buffer/BufferOp.java r378 (JTS-1.12)
+ *
+ **********************************************************************/
+
+#pragma once
+
+#include <memory> // for unique_ptr
+#include <vector> // for vector
+
+#include <geos/export.h>
+#include <geos/operation/buffer/BufferParameters.h> // for enum values
+#include <geos/geom/Geometry.h> // to materialize Geometry
+
+#include <geos/util/TopologyException.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;
+}
+}
+
+namespace geos {
+namespace operation { // geos.operation
+namespace buffer { // geos.operation.buffer
+
+/**
+ * \class BufferOp
+ *
+ * \brief
+ * Computes the buffer of a geometry, for both positive and negative
+ * buffer distances.
+ *
+ * In GIS, the positive (or negative) buffer of a geometry is defined as
+ * the Minkowski sum (or difference) of the geometry
+ * with a circle with radius equal to the absolute value of the buffer
+ * distance.
+ * In the CAD/CAM world buffers are known as *offset curves*.
+ * In morphological analysis the operation of positive and negative buffering
+ * is referred to as *erosion* and *dilation*.
+ *
+ * The buffer operation always returns a polygonal result.
+ * The negative or zero-distance buffer of lines and points is always
+ * an empty Polygon.
+ *
+ * Since true buffer curves may contain circular arcs,
+ * computed buffer polygons can only be approximations to the true geometry.
+ * The user can control the accuracy of the curve approximation by specifying
+ * the number of linear segments with which to approximate a curve.
+ *
+ * The end cap style of a linear buffer may be specified.
+ * The following end cap styles are supported:
+ * - CAP_ROUND - the usual round end caps
+ * - CAP_BUTT - end caps are truncated flat at the line ends
+ * - CAP_SQUARE - end caps are squared off at the buffer distance
+ *   beyond the line ends
+ *
+ */
+class GEOS_DLL BufferOp {
+
+
+private:
+
+    /** \brief
+     *  A number of digits of precision which leaves some computational "headroom"
+     *  for floating point operations.
+     *
+     *  This value should be less than the decimal precision of double-precision values (16).
+     */
+    static constexpr int MAX_PRECISION_DIGITS = 12;
+
+    const geom::Geometry* argGeom;
+
+    util::TopologyException saveException;
+
+    double distance;
+
+    BufferParameters bufParams;
+
+    std::unique_ptr<geom::Geometry> resultGeometry;
+
+    bool isInvertOrientation = false;
+
+    /**
+     * Compute a reasonable scale factor to limit the precision of
+     * a given combination of Geometry and buffer distance.
+     * The scale factor is based on a heuristic.
+     *
+     * @param g the Geometry being buffered
+     *
+     * @param distance the buffer distance
+     *
+     * @param maxPrecisionDigits the mzx # of digits that should be
+     *        allowed by the precision determined by the
+     *        computed scale factor
+     *
+     * @return a scale factor that allows a reasonable amount of
+     *         precision for the buffer computation
+     */
+    static double precisionScaleFactor(const geom::Geometry* g,
+                                       double distance, int maxPrecisionDigits);
+
+
+    void computeGeometry();
+
+    void bufferOriginalPrecision();
+
+    void bufferReducedPrecision(int precisionDigits);
+
+    void bufferReducedPrecision();
+
+    void bufferFixedPrecision(const geom::PrecisionModel& fixedPM);
+
+    static void extractPolygons(
+        geom::Geometry* poly0,
+        std::vector<std::unique_ptr<geom::Geometry>>& polys);
+
+public:
+
+    enum {
+        /// Specifies a round line buffer end cap style.
+        /// @deprecated use BufferParameters
+        CAP_ROUND = BufferParameters::CAP_ROUND,
+
+        /// Specifies a butt (or flat) line buffer end cap style.
+        /// @deprecated use BufferParameters
+        CAP_BUTT = BufferParameters::CAP_FLAT,
+
+        /// Specifies a square line buffer end cap style.
+        /// @deprecated use BufferParameters
+        CAP_SQUARE = BufferParameters::CAP_SQUARE
+    };
+
+    /** \brief
+     * Computes the buffer for a geometry for a given buffer distance
+     * and accuracy of approximation.
+     *
+     * @param g the geometry to buffer
+     * @param distance the buffer distance
+     * @param quadrantSegments the number of segments used to
+     *        approximate a quarter circle
+     * @param endCapStyle line buffer end cap style (default: BufferParameters::CAP_ROUND)
+     * @return the buffer of the input geometry
+     *
+     */
+    static std::unique_ptr<geom::Geometry> bufferOp(
+        const geom::Geometry* g,
+        double distance,
+        int quadrantSegments = BufferParameters::DEFAULT_QUADRANT_SEGMENTS,
+        int endCapStyle = BufferParameters::CAP_ROUND);
+
+    static std::unique_ptr<geom::Geometry> bufferOp(
+        const geom::Geometry* g,
+        double distance,
+        BufferParameters& bufParms);
+
+    /** \brief
+     * Initializes a buffer computation for the given geometry.
+     *
+     * @param g the geometry to buffer
+     */
+    BufferOp(const geom::Geometry* g)
+        :
+        argGeom(g),
+        bufParams(),
+        resultGeometry(nullptr),
+        isInvertOrientation(false)
+    {
+    }
+
+    /** \brief
+     * Initializes a buffer computation for the given geometry
+     * with the given set of parameters.
+     *
+     * @param g the geometry to buffer
+     * @param params the buffer parameters to use. This class will
+     *               copy it to private memory.
+     */
+    BufferOp(const geom::Geometry* g, const BufferParameters& params)
+        :
+        argGeom(g),
+        bufParams(params),
+        resultGeometry(nullptr),
+        isInvertOrientation(false)
+    {
+    }
+
+    /** \brief
+     * Specifies the end cap style of the generated buffer.
+     *
+     * The styles supported are CAP_ROUND, CAP_BUTT, and CAP_SQUARE.
+     * The default is CAP_ROUND.
+     *
+     * @param nEndCapStyle the end cap style to specify
+     */
+    inline void setEndCapStyle(int nEndCapStyle);
+
+    /** \brief
+     * Sets the number of segments used to approximate a angle fillet
+     *
+     * @param nQuadrantSegments the number of segments in a fillet for a quadrant
+     */
+    inline void setQuadrantSegments(int nQuadrantSegments);
+
+    /** \brief
+     * Sets whether the computed buffer should be single-sided.
+     *
+     * A single-sided buffer is constructed on only one side
+     * of each input line.
+     *
+     * The side used is determined by the sign of the buffer distance:
+     * - a positive distance indicates the left-hand side
+     * - a negative distance indicates the right-hand side
+     *
+     * The single-sided buffer of point geometries is
+     * the same as the regular buffer.
+     *
+     * The End Cap Style for single-sided buffers is
+     * always ignored,
+     * and forced to the equivalent of `CAP_FLAT`.
+     *
+     * @param isSingleSided `true` if a single-sided buffer
+     *                      should be constructed
+     */
+    inline void setSingleSided(bool isSingleSided);
+
+    /** \brief
+     * Returns the buffer computed for a geometry for a given buffer
+     * distance.
+     *
+     * @param nDistance the buffer distance
+     * @return the buffer of the input geometry
+     */
+    std::unique_ptr<geom::Geometry> getResultGeometry(double nDistance);
+
+    /**
+    * Buffers a geometry with distance zero.
+    * The result can be computed using the maximum-signed-area orientation,
+    * or by combining both orientations.
+    *
+    * This can be used to fix an invalid polygonal geometry to be valid
+    * (i.e. with no self-intersections).
+    * For some uses (e.g. fixing the result of a simplification)
+    * a better result is produced by using only the max-area orientation.
+    * Other uses (e.g. fixing geometry) require both orientations to be used.
+    *
+    * This function is for INTERNAL use only.
+    *
+    * @param geom the polygonal geometry to buffer by zero
+    * @param isBothOrientations true if both orientations of input rings should be used
+    * @return the buffered polygonal geometry
+    */
+    static std::unique_ptr<geom::Geometry> bufferByZero(
+        const geom::Geometry* geom,
+        bool isBothOrientations);
+
+};
+
+// BufferOp inlines
+void
+BufferOp::setQuadrantSegments(int q)
+{
+    bufParams.setQuadrantSegments(q);
+}
+
+void
+BufferOp::setEndCapStyle(int s)
+{
+    bufParams.setEndCapStyle((BufferParameters::EndCapStyle)s);
+}
+
+void
+BufferOp::setSingleSided(bool isSingleSided)
+{
+    bufParams.setSingleSided(isSingleSided);
+}
+
+} // namespace geos::operation::buffer
+} // namespace geos::operation
+} // namespace geos
+
+#ifdef _MSC_VER
+#pragma warning(pop)
+#endif
+

package/dist/prebuilt/linuxmusl-x64-mt-release/include/geos/operation/buffer/BufferParameters.h

@@ -0,0 +1,307 @@
+/**********************************************************************
+ *
+ * GEOS - Geometry Engine Open Source
+ * http://geos.osgeo.org
+ *
+ * Copyright (C) 2009  Sandro Santilli <strk@kbt.io>
+ *
+ * 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/BufferParameters.java r378 (JTS-1.12)
+ *
+ **********************************************************************/
+
+#pragma once
+
+#include <geos/export.h>
+
+//#include <vector>
+
+//#include <geos/algorithm/LineIntersector.h> // for composition
+//#include <geos/geom/Coordinate.h> // for composition
+//#include <geos/geom/LineSegment.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 CoordinateSequence;
+class PrecisionModel;
+}
+namespace operation {
+namespace buffer {
+class OffsetCurveVertexList;
+}
+}
+}
+
+namespace geos {
+namespace operation { // geos.operation
+namespace buffer { // geos.operation.buffer
+
+/** \brief
+ * Contains the parameters which
+ * describe how a buffer should be constructed.
+ *
+ */
+class GEOS_DLL BufferParameters {
+
+public:
+
+    /// End cap styles
+    enum EndCapStyle {
+
+        /// Specifies a round line buffer end cap style.
+        CAP_ROUND = 1,
+
+        /// Specifies a flat line buffer end cap style.
+        CAP_FLAT = 2,
+
+        /// Specifies a square line buffer end cap style.
+        CAP_SQUARE = 3
+    };
+
+    /// Join styles
+    enum JoinStyle {
+
+        /// Specifies a round join style.
+        JOIN_ROUND = 1,
+
+        /// Specifies a mitre join style.
+        JOIN_MITRE = 2,
+
+        /// Specifies a bevel join style.
+        JOIN_BEVEL = 3
+    };
+
+    /// \brief
+    /// The default number of facets into which to divide a fillet
+    /// of 90 degrees.
+    ///
+    /// A value of 8 gives less than 2% max error in the buffer distance.
+    /// For a max error of < 1%, use QS = 12.
+    /// For a max error of < 0.1%, use QS = 18.
+    ///
+    static const int DEFAULT_QUADRANT_SEGMENTS = 8;
+
+    /// The default mitre limit
+    ///
+    /// Allows fairly pointy mitres.
+    ///
+    static const double DEFAULT_MITRE_LIMIT; // 5.0 (in .cpp file)
+
+    /// Creates a default set of parameters
+    BufferParameters();
+
+    /// Creates a set of parameters with the given quadrantSegments value.
+    ///
+    /// @param quadrantSegments the number of quadrant segments to use
+    ///
+    BufferParameters(int quadrantSegments);
+
+    /// \brief
+    /// Creates a set of parameters with the
+    /// given quadrantSegments and endCapStyle values.
+    ///
+    /// @param quadrantSegments the number of quadrant segments to use
+    /// @param endCapStyle the end cap style to use
+    ///
+    BufferParameters(int quadrantSegments, EndCapStyle endCapStyle);
+
+    /// \brief
+    /// Creates a set of parameters with the
+    /// given parameter values.
+    ///
+    /// @param quadrantSegments the number of quadrant segments to use
+    /// @param endCapStyle the end cap style to use
+    /// @param joinStyle the join style to use
+    /// @param mitreLimit the mitre limit to use
+    ///
+    BufferParameters(int quadrantSegments, EndCapStyle endCapStyle,
+                     JoinStyle joinStyle, double mitreLimit);
+
+    /// Gets the number of quadrant segments which will be used
+    ///
+    /// @return the number of quadrant segments
+    ///
+    int
+    getQuadrantSegments() const
+    {
+        return quadrantSegments;
+    }
+
+    /// \brief
+    /// Sets the number of line segments used to approximate
+    /// an angle fillet in round joins.
+    ///
+    /// <tt>quadSegs</tt> determines the maximum
+    /// error in the approximation to the true buffer curve.
+    /// The default value of 8 gives less than 2% max error in the
+    /// buffer distance.
+    /// For a max error of < 1%, use QS = 12.
+    /// For a max error of < 0.1%, use QS = 18.
+    /// The error is always less than the buffer distance
+    /// (in other words, the computed buffer curve is always inside
+    ///  the true curve).
+    ///
+    /// @param quadSegs the number of segments in a fillet for a quadrant
+    void setQuadrantSegments(int quadSegs);
+
+    /// \brief
+    /// Computes the maximum distance error due to a given level
+    /// of approximation to a true arc.
+    ///
+    /// @param quadSegs the number of segments used to approximate
+    ///                 a quarter-circle
+    /// @return the error of approximation
+    ///
+    static double bufferDistanceError(int quadSegs);
+
+    /// Gets the end cap style.
+    ///
+    /// @return the end cap style
+    ///
+    EndCapStyle
+    getEndCapStyle() const
+    {
+        return endCapStyle;
+    }
+
+    /// Specifies the end cap style of the generated buffer.
+    ///
+    /// The styles supported are CAP_ROUND, CAP_BUTT,
+    /// and CAP_SQUARE.
+    ///
+    /// The default is CAP_ROUND.
+    ///
+    /// @param style the end cap style to specify
+    ///
+    void
+    setEndCapStyle(EndCapStyle style)
+    {
+        endCapStyle = style;
+    }
+
+    /// Gets the join style.
+    ///
+    /// @return the join style
+    ///
+    JoinStyle
+    getJoinStyle() const
+    {
+        return joinStyle;
+    }
+
+    /// \brief
+    /// Sets the join style for outside (reflex) corners between
+    /// line segments.
+    ///
+    /// Allowable values are JOIN_ROUND (which is the default),
+    /// JOIN_MITRE and JOIN_BEVEL.
+    ///
+    /// @param style the code for the join style
+    ///
+    void
+    setJoinStyle(JoinStyle style)
+    {
+        joinStyle = style;
+    }
+
+    /// Gets the mitre ratio limit.
+    ///
+    /// @return the limit value
+    ///
+    double
+    getMitreLimit() const
+    {
+        return mitreLimit;
+    }
+
+    /// Sets the limit on the mitre ratio used for very sharp corners.
+    ///
+    /// The mitre ratio is the ratio of the distance from the corner
+    /// to the end of the mitred offset corner.
+    /// When two line segments meet at a sharp angle,
+    /// a miter join will extend far beyond the original geometry.
+    /// (and in the extreme case will be infinitely far.)
+    /// To prevent unreasonable geometry, the mitre limit
+    /// allows controlling the maximum length of the join corner.
+    /// Corners with a ratio which exceed the limit will be beveled.
+    ///
+    /// @param limit the mitre ratio limit
+    ///
+    void
+    setMitreLimit(double limit)
+    {
+        mitreLimit = limit;
+    }
+
+    /**
+     * Sets whether the computed buffer should be single-sided.
+     * A single-sided buffer is constructed on only one side of each input line.
+     *
+     * The side used is determined by the sign of the buffer distance:
+     * - a positive distance indicates the left-hand side
+     * - a negative distance indicates the right-hand side
+     *
+     * The single-sided buffer of point geometries is
+     * the same as the regular buffer.
+     *
+     * The End Cap Style for single-sided buffers is
+     * always ignored,
+     * and forced to the equivalent of <tt>CAP_FLAT</tt>.
+     *
+     * @param p_isSingleSided true if a single-sided buffer should be constructed
+     */
+    void
+    setSingleSided(bool p_isSingleSided)
+    {
+        _isSingleSided = p_isSingleSided;
+    }
+
+    /**
+     * Tests whether the buffer is to be generated on a single side only.
+     *
+     * @return true if the generated buffer is to be single-sided
+     */
+    bool
+    isSingleSided() const
+    {
+        return _isSingleSided;
+    }
+
+
+private:
+
+    /// Defaults to DEFAULT_QUADRANT_SEGMENTS;
+    int quadrantSegments;
+
+    /// Defaults to CAP_ROUND;
+    EndCapStyle endCapStyle;
+
+    /// Defaults to JOIN_ROUND;
+    JoinStyle joinStyle;
+
+    /// Defaults to DEFAULT_MITRE_LIMIT;
+    double mitreLimit;
+
+    bool _isSingleSided;
+};
+
+} // namespace geos::operation::buffer
+} // namespace geos::operation
+} // namespace geos
+
+#ifdef _MSC_VER
+#pragma warning(pop)
+#endif
+