if (FALSE) {
library(sf)
# vpu centers
pnt <- st_cast(st_point_on_surface(nhdR::vpu_shp), "POINT")
find_vpu(pnt[1, ])
find_vpu(pnt)
find_vpu(nhdR::gull$sp$NHDWaterbody[1, ])
find_vpu(nhdR::gull$sp$NHDWaterbody)
}
Run the code above in your browser using DataLab