IsValidOp.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/PolygonTopologyAnalyzer.h>
#include <geos/operation/valid/TopologyValidationError.h>
#include <geos/util.h>
 
 
// Forward declarations
namespace geos {
namespace geom {
class CoordinateXY;
class Geometry;
class Point;
class MultiPoint;
class LineString;
class LinearRing;
class Polygon;
class MultiPolygon;
class GeometryCollection;
}
namespace algorithm {
namespace locate {
class IndexedPointInAreaLocator;
}
}
}
 
 
namespace geos {      // geos.
namespace operation { // geos.operation
namespace valid {     // geos.operation.valid
 
class GEOS_DLL IsValidOp {
    using CoordinateXY = geos::geom::CoordinateXY;
 
private:
 
    static constexpr int MIN_SIZE_LINESTRING = 2;
    static constexpr int MIN_SIZE_RING = 4;
 
    const geom::Geometry* inputGeometry;
    bool isInvertedRingValid = false;
    std::unique_ptr<TopologyValidationError> validErr;
 
    bool hasInvalidError()
    {
        return validErr != nullptr;
    }
 
    void logInvalid(int code, const geom::CoordinateXY& pt);
 
    bool isValidGeometry(const geom::Geometry* g);
 
    bool isValid(const geom::Point* g);
 
    bool isValid(const geom::MultiPoint* g);
 
    bool isValid(const geom::LineString* g);
 
    bool isValid(const geom::LinearRing* g);
 
    bool isValid(const geom::Polygon* g);
 
    bool isValid(const geom::MultiPolygon* g);
 
    bool isValid(const geom::GeometryCollection* gc);
 
    void checkCoordinatesValid(const geom::CoordinateSequence* coords);
    void checkCoordinatesValid(const geom::Polygon* poly);
    void checkRingClosed(const geom::LinearRing* ring);
    void checkRingsClosed(const geom::Polygon* poly);
    void checkRingsPointSize(const geom::Polygon* poly);
    void checkRingPointSize(const geom::LinearRing* ring);
 
    void checkTooFewPoints(const geom::LineString* line, std::size_t minSize);
 
    bool isNonRepeatedSizeAtLeast(const geom::LineString* line, std::size_t minSize);
 
    void checkAreaIntersections(PolygonTopologyAnalyzer& areaAnalyzer);
 
    void checkRingSimple(const geom::LinearRing* ring);
 
 
    void checkHolesInShell(const geom::Polygon* poly);
 
    const CoordinateXY* findHoleOutsideShellPoint(
        const geom::LinearRing* hole,
        const geom::LinearRing* shell);
 
    void checkHolesNotNested(const geom::Polygon* poly);
 
    void checkShellsNotNested(const geom::MultiPolygon* mp);
 
    void checkInteriorConnected(PolygonTopologyAnalyzer& areaAnalyzer);
 
 
public:
 
    IsValidOp(const geom::Geometry* p_inputGeometry)
        : inputGeometry(p_inputGeometry)
        , validErr(nullptr)
        {};
 
    void setSelfTouchingRingFormingHoleValid(bool p_isValid)
    {
        isInvertedRingValid = p_isValid;
    };
 
    static bool isValid(const geom::Geometry* geom)
    {
        IsValidOp ivo(geom);
        return ivo.isValid();
    };
 
    static bool isValid(const geom::CoordinateXY& coord)
    {
        return isValid(&coord);
    }
 
    bool isValid();
 
    static bool isValid(const geom::CoordinateXY* coord);
 
    const TopologyValidationError* getValidationError();
 
 
};
 
 
} // namespace geos.operation.valid
} // namespace geos.operation
} // namespace geos