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