| autoclean | R Documentation | 
Automatically clean a rivernetwork object
autoclean(rivernetwork, mouthseg, mouthvert, crs, dissolve = FALSE)
| rivernetwork | rivernetwork object from the riverdist package | 
| mouthseg | integer | 
| mouthvert | integer | 
| crs | crs string or epsg code | 
| dissolve | logical combine runs of segments with no other connections? | 
## Not run: 
library(nhdR)
library(sf)
# Generate test dataset
bbox <- data.frame(xmin = -73.33838, ymin = 41.35841,
xmax = -73.14540, ymax = 41.48593)
b0 <- sf::st_sfc(sf::st_polygon(list(rbind(
  c(bbox$xmin, bbox$ymin),
  c(bbox$xmax, bbox$ymin),
  c(bbox$xmax, bbox$ymax),
  c(bbox$xmin, bbox$ymax),
  c(bbox$xmin, bbox$ymin)))))
sf::st_crs(b0) <- 4326
b0 <- st_transform(b0, sf::st_crs(nhdR::vpu_shp))
nhd <- nhd_plus_query(poly = b0, dsn = c("NHDFlowLine"))$sp$NHDFlowLine
outlet_reach <- terminal_reaches(network = nhd,
  approve_all_dl = TRUE)
outlet_point <- st_cast(st_line_sample(outlet_reach, sample = 1), "POINT")
mouthseg  <- which(nhd$COMID == outlet_reach$comid)
res <- autoclean(nhd_rv, mouthseg, mouthvert, crs = st_crs(nhd))
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.