107 lines
3.4 KiB
C++
107 lines
3.4 KiB
C++
/**********************************************************************
|
|
*
|
|
* 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/IncrementalDelaunayTriangulator.java r524
|
|
*
|
|
**********************************************************************/
|
|
|
|
#pragma once
|
|
|
|
#include <list>
|
|
|
|
#include <geos/triangulate/quadedge/Vertex.h>
|
|
|
|
|
|
namespace geos {
|
|
namespace triangulate { //geos.triangulate
|
|
|
|
namespace quadedge {
|
|
class QuadEdge;
|
|
class QuadEdgeSubdivision;
|
|
}
|
|
|
|
/** \brief
|
|
* Computes a Delauanay Triangulation of a set of {@link quadedge::Vertex}es,
|
|
* using an incrementatal insertion algorithm.
|
|
*
|
|
* @author JTS: Martin Davis
|
|
* @author Benjamin Campbell
|
|
*/
|
|
class GEOS_DLL IncrementalDelaunayTriangulator {
|
|
private:
|
|
quadedge::QuadEdgeSubdivision* subdiv;
|
|
bool isUsingTolerance;
|
|
bool m_isForceConvex;
|
|
|
|
public:
|
|
/**
|
|
* Creates a new triangulator using the given {@link quadedge::QuadEdgeSubdivision}.
|
|
* The triangulator uses the tolerance of the supplied subdivision.
|
|
*
|
|
* @param subdiv
|
|
* a subdivision in which to build the TIN
|
|
*/
|
|
IncrementalDelaunayTriangulator(quadedge::QuadEdgeSubdivision* subdiv);
|
|
|
|
typedef std::vector<quadedge::Vertex> VertexList;
|
|
|
|
/**
|
|
* Sets whether the triangulation is forced to have a convex boundary. Because
|
|
* of the use of a finite-size frame, this condition requires special logic to
|
|
* enforce. The default is true, since this is a requirement for some uses of
|
|
* Delaunay Triangulations (such as Concave Hull generation). However, forcing
|
|
* the triangulation boundary to be convex may cause the overall frame
|
|
* triangulation to be non-Delaunay. This can cause a problem for Voronoi
|
|
* generation, so the logic can be disabled via this method.
|
|
*
|
|
* @param isForceConvex true if the triangulation boundary is forced to be convex
|
|
*/
|
|
void forceConvex(bool isForceConvex);
|
|
|
|
/**
|
|
* Inserts all sites in a collection. The inserted vertices <b>MUST</b> be
|
|
* unique up to the provided tolerance value. (i.e. no two vertices should be
|
|
* closer than the provided tolerance value). They do not have to be rounded
|
|
* to the tolerance grid, however.
|
|
*
|
|
* @param vertices a Collection of Vertex
|
|
*
|
|
* @throws LocateFailureException if the location algorithm
|
|
* fails to converge in a reasonable number of iterations
|
|
*/
|
|
void insertSites(const VertexList& vertices);
|
|
|
|
/**
|
|
* Inserts a new point into a subdivision representing a Delaunay
|
|
* triangulation, and fixes the affected edges so that the result
|
|
* is still a Delaunay triangulation.
|
|
* <p>
|
|
*
|
|
* @return a quadedge containing the inserted vertex
|
|
*/
|
|
quadedge::QuadEdge& insertSite(const quadedge::Vertex& v);
|
|
|
|
private:
|
|
|
|
bool isConcaveBoundary(const quadedge::QuadEdge& e) const;
|
|
|
|
bool isConcaveAtOrigin(const quadedge::QuadEdge& e) const;
|
|
|
|
bool isBetweenFrameAndInserted(const quadedge::QuadEdge& e, const quadedge::Vertex& vInsert) const;
|
|
};
|
|
|
|
} //namespace geos.triangulate
|
|
} //namespace goes
|
|
|