source(system.file("extdata", "walker_data.R", package = "nhdplusTools"))
test_flowline <- prepare_nhdplus(walker_flowline, 0, 0, FALSE)
test_flowline <- data.frame(
ID = test_flowline$COMID,
toID = test_flowline$toCOMID)
(order <- get_streamorder(test_flowline))
walker_flowline$order <- order
plot(sf::st_geometry(walker_flowline), lwd = walker_flowline$order, col = "blue")
Run the code above in your browser using DataLab