source(system.file("extdata", "walker_data.R", package = "nhdplusTools"))
hydro_location <- list(comid = 5329339,
reachcode = "18050005000078",
reach_meas = 30)
(pl <- get_partial_length(hydro_location, walker_flowline))
hydro_location <- sf::st_sf(hydro_location,
geom = nhdplusTools::get_hydro_location(data.frame(hydro_location),
walker_flowline))
net <- navigate_network(hydro_location,
mode = "DM", network = walker_flowline,
distance_km = 4, trim_start = TRUE)
plot(sf::st_geometry(walker_flowline[walker_flowline$COMID == hydro_location$comid,]))
plot(sf::st_geometry(hydro_location), add = TRUE)
plot(sf::st_geometry(net), add = TRUE, col = "blue", lwd = 2)
sf::st_length(net)
pl$dn
Run the code above in your browser using DataLab