PolygonTopologyAnalyzer.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2021 Paul Ramsey <pramsey@cleverelephant.ca>
 * Copyright (C) 2021 Martin Davis
 *
 * 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/operation/valid/PolygonIntersectionAnalyzer.h>
#include <geos/operation/valid/PolygonRing.h>
#include <geos/noding/BasicSegmentString.h>
 
#include <memory>
 
// Forward declarations
namespace geos {
namespace geom {
class Geometry;
class Coordinate;
}
}
 
namespace geos {      // geos.
namespace operation { // geos.operation
namespace valid {     // geos.operation.valid
 
using geos::geom::Coordinate;
using geos::geom::CoordinateSequence;
using geos::geom::Geometry;
using geos::geom::LinearRing;
 
class GEOS_DLL PolygonTopologyAnalyzer {
 
private:
 
    // const Geometry* inputGeom;
    bool isInvertedRingValid = false;
    PolygonIntersectionAnalyzer segInt;
    std::vector<PolygonRing*> polyRings;
    geom::Coordinate disconnectionPt;
 
 
    // holding area for PolygonRings and SegmentStrings so we
    // can pass around pointers with abandon
    std::deque<PolygonRing> polyRingStore;
    std::deque<noding::BasicSegmentString> segStringStore;
    // when building SegmentStrings we sometimes want
    // to use deduped CoordinateSequences so we will
    // keep the deduped ones here so they get cleaned
    // up when processing is complete
    std::vector<std::unique_ptr<CoordinateSequence>> coordSeqStore;
 
    PolygonRing* createPolygonRing(const LinearRing* p_ring);
    PolygonRing* createPolygonRing(const LinearRing* p_ring, int p_index, PolygonRing* p_shell);
 
    static std::size_t intersectingSegIndex(const CoordinateSequence* ringPts, const Coordinate* pt);
 
    static std::size_t ringIndexPrev(const CoordinateSequence* ringPts, std::size_t index);
 
    std::vector<SegmentString*> createSegmentStrings(const Geometry* geom, bool isInvertedRingValid);
 
    std::vector<PolygonRing*> getPolygonRings(const std::vector<SegmentString*>& segStrings);
 
    SegmentString* createSegString(const LinearRing* ring, const PolygonRing* polyRing);
 
    // Declare type as noncopyable
    PolygonTopologyAnalyzer(const PolygonTopologyAnalyzer& other) = delete;
    PolygonTopologyAnalyzer& operator=(const PolygonTopologyAnalyzer& rhs) = delete;
 
public:
 
    /* public */
    PolygonTopologyAnalyzer(const Geometry* geom, bool p_isInvertedRingValid);
 
    static Coordinate findSelfIntersection(const LinearRing* ring);
 
    static bool
    isSegmentInRing(const Coordinate* p0, const Coordinate* p1,
        const LinearRing* ring);
 
    static bool isIncidentSegmentInRing(const Coordinate* p0, const Coordinate* p1,
        const CoordinateSequence* ringPts);
 
 
    bool hasInvalidIntersection() {
        return segInt.isInvalid();
    }
 
    int getInvalidCode() {
        return segInt.getInvalidCode();
    }
 
    const Coordinate& getInvalidLocation() {
        return segInt.getInvalidLocation();
    }
 
    bool isInteriorDisconnected();
 
    const Coordinate& getDisconnectionLocation() const
    {
        return disconnectionPt;
    };
 
 
    void checkInteriorDisconnectedByHoleCycle();
 
    void checkInteriorDisconnectedBySelfTouch();
 
};
 
 
} // namespace geos.operation.valid
} // namespace geos.operation
} // namespace geos