x <- sf::read_sf(system.file("extdata/new_hope.gpkg", package = "hydroloom"))
y <- dplyr::select(add_toids(x), -ToNode, -FromNode)
y <- make_node_topology(y)
# just the divergences which have unique fromids in x but don't in new hope.
div <- add_toids(dplyr::select(x, COMID, FromNode, ToNode),
return_dendritic = FALSE)
div <- div[div$toid %in%
x$COMID[x$Divergence == 2],]
y <- dplyr::select(add_toids(x), -ToNode, -FromNode)
y <- make_node_topology(y, add_div = div)
Run the code above in your browser using DataLab