<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-ligatures:standardcontextual;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span lang="PT">Hi,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT"><o:p> </o:p></span></p>
<p class="MsoNormal">This (GMT Julia wrapper) works fine:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">julia> D1 = mat2ds([0 0;10 0;10 10;11 10], geom=wkbMultiPoint);<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">julia> delaunay(D1,0.0,false)<o:p></o:p></p>
<p class="MsoNormal">Vector{GMTdataset} with 2 segments<o:p></o:p></p>
<p class="MsoNormal">Show first segment. To see other segments just type its element number. E.g. D[7]<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">BoundingBox: [0.0, 10.0, 0.0, 10.0]<o:p></o:p></p>
<p class="MsoNormal">4×2 GMTdataset{Float64, 2}<o:p></o:p></p>
<p class="MsoNormal">Row │ col.1 col.2<o:p></o:p></p>
<p class="MsoNormal"> │ Float64 Float64<o:p></o:p></p>
<p class="MsoNormal"><span lang="PT">─────</span><span lang="PT" style="font-family:"Arial",sans-serif">┼</span><span lang="PT">──────────────────</span><span lang="PT"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT"> 1 │ 0.0 0.0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT"> 2 │ 10.0 0.0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT"> 3 │ 10.0 10.0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT"> 4 │ 0.0 0.0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT"><o:p> </o:p></span></p>
<p class="MsoNormal">But when I try with 3D points<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">julia> D1 = mat2ds([0 0 0;10 0 1;10 10 2;11 10 3], geom=wkbMultiPoint);<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">julia> delaunay(D1,0.0,false)<o:p></o:p></p>
<p class="MsoNormal">ERROR 6: Incompatible geometry for operation<o:p></o:p></p>
<p class="MsoNormal">ERROR 6: Incompatible geometry for operation<o:p></o:p></p>
<p class="MsoNormal">ERROR 6: Incompatible geometry for operation<o:p></o:p></p>
<p class="MsoNormal">ERROR 6: Incompatible geometry for operation<o:p></o:p></p>
<p class="MsoNormal">ERROR 6: Incompatible geometry for operation<o:p></o:p></p>
<p class="MsoNormal">ERROR 6: Incompatible geometry for operation<o:p></o:p></p>
<p class="MsoNormal">BoundingBox: [0.0, 0.0, 0.0, 0.0, 0.0, 0.0]<o:p></o:p></p>
<p class="MsoNormal">2×3 GMTdataset{Float64, 2}<o:p></o:p></p>
<p class="MsoNormal">Row │ col.1 col.2 col.3<o:p></o:p></p>
<p class="MsoNormal"> │ Float64 Float64 Float64<o:p></o:p></p>
<p class="MsoNormal">─────<span style="font-family:"Arial",sans-serif">┼</span>───────────────────────────<o:p></o:p></p>
<p class="MsoNormal"> 1 │ 0.0 0.0 0.0<o:p></o:p></p>
<p class="MsoNormal"> 2 │ 0.0 0.0 0.0<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I’ve tried with wkbMultiPointZ, wkbMultiPoint25D and same thing.
<o:p></o:p></p>
<p class="MsoNormal">Looked at the GDAL tests and there it also uses only 2D (a 3D test is commented).<o:p></o:p></p>
<p class="MsoNormal">I am doing something wrong or delaunay in GDAL is only possible with with 2D points?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Joaquim<o:p></o:p></p>
</div>
</body>
</html>