hypertidy / anglr

Mesh creation and topology for spatial data (and not just geographic)
https://hypertidy.github.io/anglr/
83 stars 10 forks source link

DEL0 is broken for polygons #138

Closed mdsumner closed 4 years ago

mdsumner commented 4 years ago

It's not working at all

mesh_plot(DEL0(silicate::minimal_mesh, max_area = .0005))

image

DEL is fine

mdsumner commented 4 years ago

I had completely left out the segments so we got a rough triangle cull from a raw-point Delaunay! Should make than an option actually just for kicks

mdsumner commented 4 years ago

Another manifestation, in case it comes up - makes multipolygons look more busted (top right)

library(gibble)
par(mfrow = c(2, 2)); plot(DEL0(mpoly[1, ], max_area =  0.06)); plot(DEL0(mpoly[1, ]), col = "grey") > plot(DEL(mpoly[1, ], max_area = 0.06)); plot(DEL(mpoly[1, ], max_area = 0.06), col = "grey")

image