source(system.file("extdata/new_hope_data.R", package = "nhdplusTools"))
x <- dplyr::select(get_tocomid(
dplyr::select(new_hope_flowline, COMID, FromNode, ToNode, Divergence, FTYPE,
AreaSqKM, LENGTHKM, GNIS_ID)
), -tonode, -fromnode)
head(y <- make_node_topology(x))
# just the divergences which have unique fromids in x but don't in new hope.
div <- get_tocomid(dplyr::select(new_hope_flowline, COMID, FromNode, ToNode),
return_dendritic = FALSE,
remove_coastal = FALSE)
div <- div[div$tocomid %in%
new_hope_flowline$COMID[new_hope_flowline$Divergence == 2],]
y <- make_node_topology(x, div)
Run the code above in your browser using DataLab