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

Joaquim Manuel Freire Luís jluis at ualg.pt
Sun Jun 25 10:38:42 PDT 2023


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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20230625/14c51650/attachment.htm>


More information about the gdal-dev mailing list