<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<o:p></o:p>Joaquim,<br>
<blockquote type="cite"
cite="mid:GVXPR04MB9926A8A43EC6FF536068F952A621A@GVXPR04MB9926.eurprd04.prod.outlook.com">
<div class="WordSection1">
<p class="MsoNormal">Looked at the GDAL tests and there it also
uses only 2D (a 3D test is commented).</p>
</div>
</blockquote>
I assume you're speaking about the PolyhedralSurface / TIN tests in
ogr_geom.py. The reason is that the code (and in particularly GEOS)
is not ready for those geometry types.<br>
<blockquote type="cite"
cite="mid:GVXPR04MB9926A8A43EC6FF536068F952A621A@GVXPR04MB9926.eurprd04.prod.outlook.com">
<div class="WordSection1">
<p class="MsoNormal"><o:p></o:p></p>
<p class="MsoNormal">I am doing something wrong or delaunay in
GDAL is only possible with with 2D points?</p>
</div>
</blockquote>
<p>Quick testing shows that it works with a multipoint z with a
recent GEOS master version (but might work with older ones
perhaps)<br>
</p>
<p>>>> g = ogr.CreateGeometryFromWkt("MULTIPOINT Z (0 0
10,0 1 10,1 1 12,1 0 10)")</p>
<p>>>> triangulation = g.DelaunayTriangulation()</p>
<p>>>> triangulation.ExportToWkt()</p>
<p>'GEOMETRYCOLLECTION (POLYGON ((0 1 10,0 0 10,1 0 10,0 1
10)),POLYGON ((0 1 10,1 0 10,1 1 12,0 1 10)))'<br>
<br>
</p>
<p>Even</p>
<p><br>
</p>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>