[gdal-dev] Can we do Delaunay with 3D points?

[Joaquim Manuel Freire Luís via gdal-dev]

Hi,

This (GMT Julia wrapper) works fine:

julia> D1 = mat2ds([0 0;10 0;10 10;11 10], geom=wkbMultiPoint);

julia> delaunay(D1,0.0,false)
Vector{GMTdataset} with 2 segments
Show first segment. To see other segments just type its element number. E.g. 
D[7]

BoundingBox: [0.0, 10.0, 0.0, 10.0]
4×2 GMTdataset{Float64, 2}
Row │   col.1    col.2
     │ Float64  Float64
─────┼──────────────────
   1 │     0.0      0.0
   2 │    10.0      0.0
   3 │    10.0     10.0
   4 │     0.0      0.0

But when I try with 3D points

julia> D1 = mat2ds([0 0 0;10 0 1;10 10 2;11 10 3], geom=wkbMultiPoint);

julia> delaunay(D1,0.0,false)
ERROR 6: Incompatible geometry for operation
ERROR 6: Incompatible geometry for operation
ERROR 6: Incompatible geometry for operation
ERROR 6: Incompatible geometry for operation
ERROR 6: Incompatible geometry for operation
ERROR 6: Incompatible geometry for operation
BoundingBox: [0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
2×3 GMTdataset{Float64, 2}
Row │   col.1    col.2    col.3
     │ Float64  Float64  Float64
─────┼───────────────────────────
   1 │     0.0      0.0      0.0
   2 │     0.0      0.0      0.0


I’ve tried with wkbMultiPointZ, wkbMultiPoint25D and same thing.
Looked at the GDAL tests and there it also uses only 2D (a 3D test is 
commented).
I am doing something wrong or delaunay in GDAL is only possible with with 2D 
points?

Joaquim

_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev