hypertidy / anglr

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

why doesnt copy_down auto-reproject? #92

Closed mdsumner closed 4 years ago

mdsumner commented 4 years ago

WHY

elev_raster <- quadmesh::etopo
library(sf)
poly <- subset(silicate::inlandwaters, Province %in% c("Victoria", "South Australia"))
## should get a warning here as it detects LCC on poly, longlat on elev_raster
copy_down(anglr::DEL(poly, max_area = 1e12), elev_raster)

89

mdsumner commented 4 years ago

it's because of #91 - when you subset with sf it rebuilds and moves your projstring to $input, and adds WKT2

so #91 becomes absolutely essential, until we can get PROJ sorted

mdsumner commented 4 years ago

FIXED. LOTS OF TRAPS STAY AWAY. Actually things will settle down in a year or two