Convert a bounding box to polygon
bbox2poly(bbox)
bbox: object of class bbox from sfAn sfc object from the sf package
## Not run: library(sf) wk <- wikilake::lake_wiki("Gull Lake (Michigan)") pnt <- st_as_sf(wk, coords = c("Lon", "Lat"), crs = 4326) pnt <- st_transform(pnt, st_crs(vpu_shp)) qry <- nhd_plus_query(wk$Lon, wk$Lat, dsn = c("NHDWaterbody"), buffer_dist = 0.05) wbd <- qry$sp$NHDWaterbody[which.max(st_area(qry$sp$NHDWaterbody)), ] bbox2poly(st_bbox(wbd)) ## End(Not run)
Related functions from the same R package