TPVWSimplifier.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2021 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 <queue>
 
#include <geos/index/strtree/TemplateSTRtree.h>
#include <geos/index/VertexSequencePackedRtree.h>
#include <geos/simplify/LinkedLine.h>
#include <geos/coverage/Corner.h>
#include <geos/export.h>
#include <geos/util/Progress.h>
 
 
namespace geos {
namespace geom {
class Coordinate;
class CoordinateSequence;
class Envelope;
class Geometry;
class GeometryFactory;
class LineString;
class MultiLineString;
}
}
 
namespace geos {
namespace coverage { // geos::coverage
 

class GEOS_DLL TPVWSimplifier
{
    using Corner = geos::coverage::Corner;
    using Coordinate = geos::geom::Coordinate;
    using CoordinateSequence = geos::geom::CoordinateSequence;
    using Envelope = geos::geom::Envelope;
    using Geometry = geos::geom::Geometry;
    using GeometryFactory = geos::geom::GeometryFactory;
    using LineString = geos::geom::LineString;
    using MultiLineString = geos::geom::MultiLineString;
    using VertexSequencePackedRtree = geos::index::VertexSequencePackedRtree;
    template<typename ItemType>
    using TemplateSTRtree = geos::index::strtree::TemplateSTRtree<ItemType>;
    using LinkedLine = geos::simplify::LinkedLine;
 
public:
 
    // Prototype
    class EdgeIndex;
 
    /* private static */
    class Edge {
 
    public:
 
        // Members
        double areaTolerance;
        bool isFreeRing;
        const Envelope* envelope;
        std::size_t nbPts;
        LinkedLine linkedLine;
        VertexSequencePackedRtree vertexIndex;
        std::size_t minEdgeSize;

        Edge(const LineString* p_inputLine, bool p_isFreeRing, double p_areaTolerance);
 
        const Coordinate& getCoordinate(std::size_t index) const;
 
        const Envelope* getEnvelopeInternal() const;
 
        std::size_t size() const;
 
        std::unique_ptr<CoordinateSequence> simplify(EdgeIndex& edgeIndex);
 
        void createQueue(Corner::PriorityQueue& pq);
 
        void addCorner(std::size_t i, Corner::PriorityQueue& cornerQueue);
 
        bool isRemovable(Corner& corner, EdgeIndex& edgeIndex) const;

        bool hasIntersectingVertex(const Corner& corner,
            const Envelope& cornerEnv,
            const Edge& edge) const;
 
        std::vector<std::size_t> query(const Envelope& cornerEnv) const;

        void removeCorner(
            Corner& corner,
            Corner::PriorityQueue& cornerQueue);
 
    }; // Edge
 
    class EdgeIndex
    {
        public:
 
            TemplateSTRtree<const Edge*> index;
 
            void add(std::vector<Edge>& edges);
 
            std::vector<const Edge*> query(const Envelope& queryEnv);
 
    }; // EdgeIndex
 

    static std::unique_ptr<MultiLineString> simplify(
        const MultiLineString* lines,
        double distanceTolerance,
        geos::util::ProgressFunction* progressFunction);

    static std::unique_ptr<MultiLineString> simplify(
        const MultiLineString* lines,
        std::vector<bool>& freeRings,
        const MultiLineString* constraintLines,
        double distanceTolerance,
        geos::util::ProgressFunction* progressFunction);
 
    // Constructor
    TPVWSimplifier(const MultiLineString* lines,
        double distanceTolerance);
 
 
private:
 
    // Members
    const MultiLineString* inputLines;
    std::vector<bool> isFreeRing;
    double areaTolerance;
    const GeometryFactory* geomFactory;
    const MultiLineString* constraintLines;
 
 
    // Methods
    void setFreeRingIndices(std::vector<bool>& freeRing);
 
    void setConstraints(const MultiLineString* constraints);
 
    std::unique_ptr<MultiLineString> simplify(geos::util::ProgressFunction* progressFunction);
 
    std::vector<Edge> createEdges(
        const MultiLineString* lines,
        std::vector<bool>& freeRing,
        geos::util::ProgressFunction* progressFunction);
 
 
}; // TPVWSimplifier
 
 
} // geos::coverage
} // geos