VoronoiDiagramBuilder.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2012 Excensus LLC.
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU Lesser General Licence as published
 * by the Free Software Foundation.
 * See the COPYING file for more information.
 *
 **********************************************************************
 *
 * Last port: triangulate/VoronoiDiagramBuilder.java r524
 *
 **********************************************************************/
 
#pragma once
 
#include <geos/triangulate/quadedge/QuadEdgeSubdivision.h>
#include <geos/geom/Envelope.h> // for composition
#include <memory>
#include <iostream>
#include <unordered_map>
 
namespace geos {
namespace geom {
class Geometry;
class CoordinateSequence;
class GeometryCollection;
class GeometryFactory;
}
namespace triangulate { //geos.triangulate
 
class GEOS_DLL VoronoiDiagramBuilder {
public:
    VoronoiDiagramBuilder();
 
    ~VoronoiDiagramBuilder() = default;
 
    void setSites(const geom::Geometry& geom);
 
    void setSites(const geom::CoordinateSequence& coords);
 
    void setOrdered(bool isOrdered);
 
    void setClipEnvelope(const geom::Envelope* clipEnv);
 
    void setTolerance(double tolerance);
 
    std::unique_ptr<quadedge::QuadEdgeSubdivision> getSubdivision();
 
    std::unique_ptr<geom::GeometryCollection> getDiagram(const geom::GeometryFactory& geomFact);
 
    std::unique_ptr<geom::MultiLineString> getDiagramEdges(const geom::GeometryFactory& geomFact);
 
    void reorderCellsToInput(std::vector<std::unique_ptr<geom::Geometry>> & polys) const;
 
private:
    using CoordinateCellMap = std::unordered_map<geom::CoordinateXY, std::unique_ptr<geom::Geometry>, geom::Coordinate::HashCode>;
 
    std::unique_ptr<geom::CoordinateSequence> siteCoords;
    double tolerance;
    std::unique_ptr<quadedge::QuadEdgeSubdivision> subdiv;
    const geom::Envelope* clipEnv; // externally owned
    const geom::Geometry* inputGeom;
    const geom::CoordinateSequence* inputSeq;
    geom::Envelope diagramEnv;
    bool isOrdered;
 
    void create();
 
    std::size_t getNumInputPoints() const;
 
    static std::unique_ptr<geom::GeometryCollection>
    clipGeometryCollection(std::vector<std::unique_ptr<geom::Geometry>> & geoms, const geom::Envelope& clipEnv);
 
 
    static void addCellsForCoordinates(CoordinateCellMap& cellMap, const geom::Geometry& g, std::vector<std::unique_ptr<geom::Geometry>> & polys);
    static void addCellsForCoordinates(CoordinateCellMap& cellMap, const geom::CoordinateSequence& g, std::vector<std::unique_ptr<geom::Geometry>> & polys);
};
 
} //namespace geos.triangulate
} //namespace geos