Automatically clean a rivernetwork object

autoclean(rivernetwork, mouthseg, mouthvert, crs, dissolve = FALSE)

Arguments

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?

Examples

if (FALSE) {
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))
}