[geos-commits] r3951 - in trunk: include/geos/triangulate src/triangulate tests/unit tests/unit/triangulate
svn_geos at osgeo.org
svn_geos at osgeo.org
Tue Sep 10 08:08:01 PDT 2013
Author: strk
Date: 2013-09-10 08:08:01 -0700 (Tue, 10 Sep 2013)
New Revision: 3951
Added:
trunk/include/geos/triangulate/VoronoiDiagramBuilder.h
trunk/src/triangulate/VoronoiDiagramBuilder.cpp
trunk/tests/unit/triangulate/VoronoiTest.cpp
Modified:
trunk/include/geos/triangulate/Makefile.am
trunk/src/triangulate/Makefile.am
trunk/tests/unit/Makefile.am
Log:
VoronoiDigramBuilder class + test
Contributed by Vishal Tiwari
See https://github.com/libgeos/libgeos/pull/25
Modified: trunk/include/geos/triangulate/Makefile.am
===================================================================
--- trunk/include/geos/triangulate/Makefile.am 2013-09-10 07:18:48 UTC (rev 3950)
+++ trunk/include/geos/triangulate/Makefile.am 2013-09-10 15:08:01 UTC (rev 3951)
@@ -10,4 +10,5 @@
geos_HEADERS = \
IncrementalDelaunayTriangulator.h \
- DelaunayTriangulationBuilder.h
+ DelaunayTriangulationBuilder.h \
+ VoronoiDiagramBuilder.h
Added: trunk/include/geos/triangulate/VoronoiDiagramBuilder.h
===================================================================
--- trunk/include/geos/triangulate/VoronoiDiagramBuilder.h (rev 0)
+++ trunk/include/geos/triangulate/VoronoiDiagramBuilder.h 2013-09-10 15:08:01 UTC (rev 3951)
@@ -0,0 +1,124 @@
+/**********************************************************************
+ *
+ * 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
+ *
+ **********************************************************************/
+
+#ifndef GEOS_TRIANGULATE_VORONOIDIAGRAMBUILDER_H
+#define GEOS_TRIANGULATE_VORONOIDIAGRAMBUILDER_H
+
+#include <geos/triangulate/quadedge/QuadEdgeSubdivision.h>
+#include <memory>
+#include <iostream>
+
+namespace geos {
+ namespace geom{
+ class Geometry;
+ class CoordinateSequence;
+ class GeometryCollection;
+ class GeometryFactory;
+ class Envelope;
+}
+namespace triangulate { //geos.triangulate
+
+/**
+ * A utility class which creates Voronoi Diagrams
+ * from collections of points.
+ * The diagram is returned as a {@link GeometryCollection} of {@link Polygon}s,
+ * clipped to the larger of a supplied envelope or to an envelope determined
+ * by the input sites.
+ *
+ * @author Martin Davis
+ *
+ */
+class GEOS_DLL VoronoiDiagramBuilder{
+private:
+ geom::CoordinateSequence* siteCoords;
+ double tolerance;
+ quadedge::QuadEdgeSubdivision* subdiv;
+ geom::Envelope* clipEnv;
+ geom::Envelope* diagramEnv;
+
+public:
+ /**
+ * Creates a new Voronoi diagram builder.
+ *
+ */
+ VoronoiDiagramBuilder();
+
+ ~VoronoiDiagramBuilder();
+
+ /**
+ * Sets the sites (point or vertices) which will be diagrammed.
+ * All vertices of the given geometry will be used as sites.
+ *
+ * @param geom the geometry from which the sites will be extracted.
+ */
+ void setSites(const geom::Geometry& geom);
+
+ /**
+ * Sets the sites (point or vertices) which will be diagrammed
+ * from a collection of {@link Coordinate}s.
+ *
+ * @param coords a collection of Coordinates.
+ */
+ void setSites(const geom::CoordinateSequence& coords);
+
+ /**
+ * Sets the envelope to clip the diagram to.
+ * The diagram will be clipped to the larger
+ * of this envelope or an envelope surrounding the sites.
+ *
+ * @param clipEnv the clip envelope.
+ */
+ void setClipEnvelope(const geom::Envelope& clipEnv);
+
+ /**
+ * Sets the snapping tolerance which will be used
+ * to improved the robustness of the triangulation computation.
+ * A tolerance of 0.0 specifies that no snapping will take place.
+ *
+ * @param tolerance the tolerance distance to use
+ */
+ void setTolerance(const double tolerance);
+
+ /**
+ * Gets the {@link QuadEdgeSubdivision} which models the computed diagram.
+ *
+ * @return the subdivision containing the triangulation
+ */
+ quadedge::QuadEdgeSubdivision* getSubdivision();
+
+ /**
+ * Gets the faces of the computed diagram as a {@link GeometryCollection}
+ * of {@link Polygon}s, clipped as specified.
+ *
+ * @param geomFact the geometry factory to use to create the output
+ * @return the faces of the diagram
+ */
+ std::auto_ptr<geom::GeometryCollection> getDiagram(const geom::GeometryFactory& geomFact);
+
+private:
+ void create();
+
+ static std::auto_ptr<geom::GeometryCollection>
+ clipGeometryCollection(const geom::GeometryCollection& geom, const geom::Envelope& clipEnv);
+
+};
+
+} //namespace geos.triangulate
+} //namespace geos
+
+#endif //GEOS_TRIANGULATE_VORONOIDIAGRAMBUILDER_H
Modified: trunk/src/triangulate/Makefile.am
===================================================================
--- trunk/src/triangulate/Makefile.am 2013-09-10 07:18:48 UTC (rev 3950)
+++ trunk/src/triangulate/Makefile.am 2013-09-10 15:08:01 UTC (rev 3951)
@@ -10,7 +10,8 @@
libtriangulate_la_SOURCES = \
IncrementalDelaunayTriangulator.cpp \
- DelaunayTriangulationBuilder.cpp
+ DelaunayTriangulationBuilder.cpp \
+ VoronoiDiagramBuilder.cpp
libtriangulate_la_LIBADD = \
quadedge/libquadedge.la
Added: trunk/src/triangulate/VoronoiDiagramBuilder.cpp
===================================================================
--- trunk/src/triangulate/VoronoiDiagramBuilder.cpp (rev 0)
+++ trunk/src/triangulate/VoronoiDiagramBuilder.cpp 2013-09-10 15:08:01 UTC (rev 3951)
@@ -0,0 +1,151 @@
+/**********************************************************************
+ *
+ * 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 rev. r524
+ *
+ **********************************************************************/
+
+#include <geos/triangulate/VoronoiDiagramBuilder.h>
+
+#include <algorithm>
+#include <math.h>
+#include <vector>
+#include <iostream>
+
+#include <geos/geom/GeometryFactory.h>
+#include <geos/geom/Coordinate.h>
+#include <geos/geom/Envelope.h>
+#include <geos/geom/CoordinateSequence.h>
+#include <geos/triangulate/IncrementalDelaunayTriangulator.h>
+#include <geos/triangulate/DelaunayTriangulationBuilder.h>
+#include <geos/triangulate/quadedge/QuadEdgeSubdivision.h>
+// WKT writer
+#include <geos/io/WKTWriter.h>
+
+namespace geos {
+namespace triangulate { //geos.triangulate
+
+using namespace geos::geom;
+
+
+VoronoiDiagramBuilder::VoronoiDiagramBuilder() :
+ siteCoords(NULL), tolerance(0.0), subdiv(NULL) , clipEnv(NULL), diagramEnv(NULL)
+{
+}
+
+VoronoiDiagramBuilder::~VoronoiDiagramBuilder()
+{
+ if(siteCoords)
+ delete siteCoords;
+ if(subdiv)
+ delete subdiv;
+ if(clipEnv)
+ delete clipEnv;
+ if(diagramEnv)
+ delete diagramEnv;
+}
+
+void
+VoronoiDiagramBuilder::setSites(const geom::Geometry& geom)
+{
+ siteCoords = DelaunayTriangulationBuilder::extractUniqueCoordinates(geom);
+}
+
+void
+VoronoiDiagramBuilder::setSites(const geom::CoordinateSequence& coords)
+{
+ siteCoords = coords.clone();
+ DelaunayTriangulationBuilder::unique(*siteCoords);
+}
+
+void
+VoronoiDiagramBuilder::setClipEnvelope(const geom::Envelope& clipEnv)
+{
+ *(this->clipEnv) = clipEnv;
+}
+
+void
+VoronoiDiagramBuilder::setTolerance(const double tolerance)
+{
+ this->tolerance = tolerance;
+}
+
+void
+VoronoiDiagramBuilder::create()
+{
+ if(subdiv!=NULL)
+ return;
+ geom::Envelope siteEnv = DelaunayTriangulationBuilder::envelope(*siteCoords);
+ diagramEnv = new Envelope();
+ *diagramEnv = siteEnv;
+ //adding buffer around the final envelope
+ double expandBy = fmax(diagramEnv->getWidth() , diagramEnv->getHeight());
+ diagramEnv->expandBy(expandBy);
+ if(clipEnv!=NULL)
+ diagramEnv->expandToInclude(clipEnv);
+
+ IncrementalDelaunayTriangulator::VertexList* vertices = DelaunayTriangulationBuilder::toVertices(*siteCoords);
+
+ subdiv = new quadedge::QuadEdgeSubdivision(*diagramEnv,tolerance);
+ IncrementalDelaunayTriangulator triangulator(subdiv);
+ triangulator.insertSites(*vertices);
+ delete vertices;
+}
+
+quadedge::QuadEdgeSubdivision*
+VoronoiDiagramBuilder::getSubdivision()
+{
+ create();
+ return subdiv;
+}
+
+std::auto_ptr<geom::GeometryCollection>
+VoronoiDiagramBuilder::getDiagram(const geom::GeometryFactory& geomFact)
+{
+ create();
+ std::auto_ptr<geom::GeometryCollection> polys = subdiv->getVoronoiDiagram(geomFact);
+ return clipGeometryCollection(*polys,*diagramEnv);
+}
+
+std::auto_ptr<geom::GeometryCollection>
+VoronoiDiagramBuilder::clipGeometryCollection(const geom::GeometryCollection& geom, const geom::Envelope& clipEnv)
+{
+ geom::Geometry* clipPoly = geom.getFactory()->toGeometry(&clipEnv);
+ std::auto_ptr< std::vector<Geometry*> >clipped(new std::vector<Geometry*>);
+ for(std::size_t i=0 ; i < geom.getNumGeometries() ; i++)
+ {
+ Geometry* g = (Geometry*)geom.getGeometryN(i);
+ Geometry* result=NULL;
+ // don't clip unless necessary
+ if(clipEnv.contains(g->getEnvelopeInternal()))
+ {
+ result = g->clone();
+ }
+ else if(clipEnv.intersects(g->getEnvelopeInternal()))
+ {
+ result = clipPoly->intersection(g);
+ result->setUserData(g->getUserData());
+ }
+
+ if(result!=NULL && !result->isEmpty() )
+ {
+ clipped->push_back(result);
+ }
+ }
+ delete clipPoly;
+ return std::auto_ptr<GeometryCollection>(geom.getFactory()->createGeometryCollection(clipped.release()));
+}
+
+} //namespace geos.triangulate
+} //namespace geos
Modified: trunk/tests/unit/Makefile.am
===================================================================
--- trunk/tests/unit/Makefile.am 2013-09-10 07:18:48 UTC (rev 3950)
+++ trunk/tests/unit/Makefile.am 2013-09-10 15:08:01 UTC (rev 3951)
@@ -107,6 +107,7 @@
triangulate/quadedge/QuadEdgeSubdivisionTest.cpp \
triangulate/quadedge/VertexTest.cpp \
triangulate/DelaunayTest.cpp \
+ triangulate/VoronoiTest.cpp \
util/UniqueCoordinateArrayFilterTest.cpp \
capi/GEOSCoordSeqTest.cpp \
capi/GEOSDelaunayTriangulationTest.cpp \
Added: trunk/tests/unit/triangulate/VoronoiTest.cpp
===================================================================
--- trunk/tests/unit/triangulate/VoronoiTest.cpp (rev 0)
+++ trunk/tests/unit/triangulate/VoronoiTest.cpp 2013-09-10 15:08:01 UTC (rev 3951)
@@ -0,0 +1,187 @@
+
+//
+// Test Suite for geos::triangulate::Voronoi
+//
+// tut
+#include <tut.hpp>
+// geos
+#include <geos/triangulate/quadedge/QuadEdge.h>
+#include <geos/triangulate/quadedge/QuadEdgeSubdivision.h>
+#include <geos/triangulate/IncrementalDelaunayTriangulator.h>
+#include <geos/triangulate/VoronoiDiagramBuilder.h>
+
+#include <geos/io/WKTWriter.h>
+#include <geos/io/WKTReader.h>
+#include <geos/geom/GeometryCollection.h>
+#include <geos/geom/GeometryFactory.h>
+
+#include <geos/geom/CoordinateArraySequence.h>
+//#include <stdio.h>
+#include <iostream>
+using namespace std;
+using namespace geos::triangulate;
+using namespace geos::triangulate::quadedge;
+using namespace geos::geom;
+using namespace geos::io;
+
+namespace tut
+{
+ //
+ // Test Group
+ //
+
+ // dummy data, not used
+ struct test_voronoidiag_data
+ {
+ test_voronoidiag_data()
+ {
+ }
+ };
+
+ typedef test_group<test_voronoidiag_data> group;
+ typedef group::object object;
+
+ group test_voronoidiag_group("geos::triangulate::Voronoi");
+
+ //helper function for funning triangulation
+ void runVoronoi(const char *sitesWkt, const char *expectedWkt , const double tolerance)
+ {
+ WKTReader reader;
+ WKTWriter writer;
+ geos::triangulate::VoronoiDiagramBuilder builder;
+ Geometry* sites = reader.read(sitesWkt);
+ Geometry* expected = reader.read(expectedWkt);
+ std::auto_ptr<GeometryCollection> results;
+ GeometryFactory geomFact;
+ builder.setSites(*sites);
+
+ //set Tolerance:
+ builder.setTolerance(tolerance);
+ results = builder.getDiagram(geomFact);
+
+ results->normalize();
+ expected->normalize();
+
+ ensure(results->equalsExact(expected, 1e-7));
+ ensure(results->getCoordinateDimension() == expected->getCoordinateDimension());
+
+ delete sites;
+ delete expected;
+ }
+
+ // Test Cases
+ // Basic Tests.
+ template<>
+ template<>
+ void object::test<1>()
+ {
+ using geos::geom::CoordinateSequence;
+ Coordinate a(180,300);
+ Coordinate b(300,290);
+ Coordinate c(230,330);
+ Coordinate d(244,284);
+
+ geos::triangulate::VoronoiDiagramBuilder builder;
+ std::auto_ptr< std::vector<Coordinate> > v(new std::vector<Coordinate>());
+ v->push_back(a);
+ v->push_back(b);
+ v->push_back(c);
+ v->push_back(d);
+
+ geos::geom::CoordinateArraySequence *seq = new CoordinateArraySequence(v.release());
+ builder.setSites(*seq);
+
+ //getting the subdiv()
+ QuadEdgeSubdivision* subdiv = builder.getSubdivision();
+
+ ensure_equals(subdiv->getTolerance() , 0);
+ ensure_equals(subdiv->getEnvelope().toString(),"Env[-3540:4020,-3436:4050]");
+ delete seq;
+
+ }
+ // 1 - Case with a single point
+ template<>
+ template<>
+ void object::test<2>()
+ {
+ const char *wkt = "MULTIPOINT ((150 200))";
+ const char *expected = "GEOMETRYCOLLECTION EMPTY";
+ runVoronoi(wkt,expected,0);
+ }
+
+ template<>
+ template<>
+ void object::test<3>()
+ {
+ const char *wkt = "MULTIPOINT ((150 200), (180 270), (275 163))";
+
+ const char *expected = "GEOMETRYCOLLECTION (POLYGON ((25 38, 25 295, 221.20588235294116 210.91176470588235, 170.024 38, 25 38)), POLYGON ((400 369.6542056074766, 400 38, 170.024 38, 221.20588235294116 210.91176470588235, 400 369.6542056074766)), POLYGON ((25 295, 25 395, 400 395, 400 369.6542056074766, 221.20588235294116 210.91176470588235, 25 295)))";
+
+ runVoronoi(wkt,expected,0);
+ }
+
+ template<>
+ template<>
+ void object::test<4>()
+ {
+ const char *wkt = "MULTIPOINT ((280 300), (420 330), (380 230), (320 160))";
+
+ const char *expected = "GEOMETRYCOLLECTION (POLYGON ((110 175.71428571428572, 110 500, 310.35714285714283 500, 353.515625 298.59375, 306.875 231.96428571428572, 110 175.71428571428572)), POLYGON ((590 204, 590 -10, 589.1666666666666 -10, 306.875 231.96428571428572, 353.515625 298.59375, 590 204)), POLYGON ((110 -10, 110 175.71428571428572, 306.875 231.96428571428572, 589.1666666666666 -10, 110 -10)), POLYGON ((310.35714285714283 500, 590 500, 590 204, 353.515625 298.59375, 310.35714285714283 500)))";
+ runVoronoi(wkt,expected,0);
+ }
+
+ template<>
+ template<>
+ void object::test<5>()
+ {
+ const char *wkt = "MULTIPOINT ((320 170), (366 246), (530 230), (530 300), (455 277), (490 160))";
+
+ const char *expected = "GEOMETRYCOLLECTION (POLYGON ((110 -50, 110 349.02631578947364, 405.31091180866963 170.28550074738416, 392.35294117647055 -50, 110 -50)), POLYGON ((740 63.57142857142859, 740 -50, 392.35294117647055 -50, 405.31091180866963 170.28550074738416, 429.9147677857019 205.76082797008175, 470.12061711079946 217.7882187938289, 740 63.57142857142859)), POLYGON ((110 349.02631578947364, 110 510, 323.9438202247191 510, 429.9147677857019 205.76082797008175, 405.31091180866963 170.28550074738416, 110 349.02631578947364)), POLYGON ((323.9438202247191 510, 424.57333333333327 510, 499.70666666666665 265, 470.12061711079946 217.7882187938289, 429.9147677857019 205.76082797008175, 323.9438202247191 510)),POLYGON ((740 265, 740 63.57142857142859, 470.12061711079946 217.7882187938289, 499.70666666666665 265, 740 265)), POLYGON ((424.57333333333327 510, 740 510, 740 265, 499.70666666666665 265, 424.57333333333327 510)))";
+ runVoronoi(wkt,expected,0);
+ }
+ //6. A little larger number of points
+ template<>
+ template<>
+ void object::test<6>()
+ {
+ const char *wkt = "MULTIPOINT ((280 200), (406 285), (580 280), (550 190), (370 190), (360 90), (480 110), (440 160), (450 180), (480 180), (460 160), (360 210), (360 220), (370 210), (375 227))";
+
+ const char *expected = "GEOMETRYCOLLECTION (POLYGON ((-20 -102.27272727272727, -20 585, 111.9484126984127 585, 293.54906542056074 315.803738317757, 318.75 215, 323.2352941176471 179.11764705882354, 319.3956043956044 144.56043956043956, -20 -102.27272727272727)), POLYGON ((365 200, 365 215, 369.40909090909093 219.4090909090909, 414.2119205298013 206.2317880794702, 411.875 200, 365 200)), POLYGON ((365 215, 365 200, 323.2352941176471 179.11764705882354, 318.75 215, 365 215)), POLYGON ((-20 -210, -20 -102.27272727272727, 319.3956043956044 144.56043956043956, 388.972602739726 137.60273972602738, 419.55882352941177 102.64705882352942, 471.66666666666674 -210, -20 -210)), POLYGON ((319.3956043956044 144.56043956043956, 323.2352941176471 179.11764705882354, 365 200, 411.875 200, 410.29411764705884 187.35294117647058, 388.972602739726 137.60273972602738, 319.3956043956044 144.56043956043956)), POLYGON ((410.29411764705884 187.35294117647058, 411.875 200, 414.2119205298013 206.2317
880794702, 431.62536593766146 234.01920096435336, 465 248.0047619047619, 465 175, 450 167.5, 410.29411764705884 187.35294117647058)), POLYGON ((365 215, 318.75 215, 293.54906542056074 315.803738317757, 339.6500765696784 283.1784073506891, 369.40909090909093 219.4090909090909, 365 215)), POLYGON ((111.9484126984127 585, 501.69252873563215 585, 492.5670391061452 267.4329608938547, 465 248.0047619047619, 431.62536593766146 234.01920096435336, 339.6500765696784 283.1784073506891, 293.54906542056074 315.803738317757, 111.9484126984127 585)), POLYGON ((369.40909090909093 219.4090909090909, 339.6500765696784 283.1784073506891, 431.62536593766146 234.01920096435336, 414.2119205298013 206.2317880794702, 369.40909090909093 219.4090909090909)), POLYGON ((388.972602739726 137.60273972602738, 410.29411764705884 187.35294117647058, 450 167.5, 450 127, 419.55882352941177 102.64705882352942, 388.972602739726 137.60273972602738)), POLYGON ((465 175, 465 248.0047619047619, 492.5670391061452
267.4329608938547, 505 255, 520.7142857142857 145, 495 145, 465 175)),POLYGON ((880 -169.375, 880 -210, 471.66666666666674 -210, 419.55882352941177 102.64705882352942, 450 127, 495 145, 520.7142857142857 145, 880 -169.375)), POLYGON ((450 167.5, 465 175, 495 145, 450 127, 450 167.5)), POLYGON ((501.69252873563215 585, 880 585, 880 130.00000000000006, 505 255, 492.5670391061452 267.4329608938547, 501.69252873563215 585)), POLYGON ((880 130.00000000000006, 880 -169.375, 520.7142857142857 145, 505 255, 880 130.00000000000006)))";
+
+ runVoronoi(wkt,expected,0);
+ }
+ //Test with case of Tolerance:
+ template<>
+ template<>
+ void object::test<7>()
+ {
+ const char *wkt = "MULTIPOINT ((100 200), (105 202), (110 200), (140 230), (210 240), (220 190), (170 170), (170 260), (213 245), (220 190))";
+
+ const char *expected = "GEOMETRYCOLLECTION (POLYGON ((-20 50, -20 380, -3.75 380, 105 235, 105 115, 77.14285714285714 50, -20 50)), POLYGON ((247 50, 77.14285714285714 50, 105 115, 145 195, 178.33333333333334 211.66666666666666, 183.51851851851853 208.7037037037037, 247 50)), POLYGON ((-3.75 380, 20.000000000000007 380, 176.66666666666666 223.33333333333334, 178.33333333333334 211.66666666666666, 145 195, 105 235, -3.75 380)), POLYGON ((105 115, 105 235, 145 195, 105 115)), POLYGON ((20.000000000000007 380, 255 380, 176.66666666666666 223.33333333333334, 20.000000000000007 380)), POLYGON ((255 380, 340 380, 340 240, 183.51851851851853 208.7037037037037, 178.33333333333334 211.66666666666666, 176.66666666666666 223.33333333333334, 255 380)), POLYGON ((340 240, 340 50, 247 50, 183.51851851851853 208.7037037037037, 340 240)))";
+
+ runVoronoi(wkt,expected,6);
+ }
+ template<>
+ template<>
+ void object::test<8>()
+ {
+ const char *wkt = "MULTIPOINT ((170 270), (177 275), (190 230), (230 250), (210 290), (240 280), (240 250))";
+
+ const char *expected = "GEOMETRYCOLLECTION (POLYGON ((100 210, 100 360, 150 360, 200 260, 100 210)), POLYGON ((150 360, 250 360, 220 270, 200 260, 150 360)), POLYGON ((100 160, 100 210, 200 260, 235 190, 247 160, 100 160)), POLYGON ((220 270, 235 265, 235 190, 200 260, 220 270)), POLYGON ((250 360, 310 360, 310 265, 235 265, 220 270, 250 360)), POLYGON ((310 265, 310 160, 247 160, 235 190, 235 265, 310 265)))";
+
+ runVoronoi(wkt,expected,10);
+ }
+ //Taking tolerance very very high
+ template<>
+ template<>
+ void object::test<9>()
+ {
+ const char *wkt = "MULTIPOINT ((155 271), (150 360), (260 360), (271 265), (280 260), (270 370), (154 354), (150 260))";
+
+ const char *expected = "GEOMETRYCOLLECTION (POLYGON ((20 130, 20 310, 205 310, 215 299, 215 130, 20 130)), POLYGON ((205 500, 410 500, 410 338, 215 299, 205 310, 205 500)), POLYGON ((20 310, 20 500, 205 500, 205 310, 20 310)), POLYGON ((410 338, 410 130, 215 130, 215 299, 410 338)))";
+
+ runVoronoi(wkt,expected,100);
+ }
+
+} // namespace tut
+
More information about the geos-commits
mailing list