[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