[geos-devel] Order of results in GEOSDelaunayTriangulation

[梁天辰]

a possible method to determine a mapping from triangles to nodes: 
1.edit geos-3.6.1\src\triangulate\quadedge\QuadEdgeSubdivision.cpp line 588 ,add:	
		cellPoly->setUserData(reinterpret_cast<void*>(geomFact.createPoint(c)));
		return cellPoly;
2.edit geos-3.6.1\src\triangulate\VoronoiDiagramBuilder.cpp line 136 ,moved one line:
		else if(clipEnv.intersects(g->getEnvelopeInternal()))
		{
			result.reset( clipPoly->intersection(g) );
		}
		
		if(result.get() && !result->isEmpty() )
		{result->setUserData(((Geometry*)g)->getUserData()); // moved
			clipped->push_back(result.release());
		}

hope this will help.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/geos-devel/attachments/20170908/161b002a/attachment.html>