source(system.file("extdata/sample_data.R", package = "nhdplusTools"))
fline <- sf::read_sf(sample_data, "NHDFlowline_Network")
start <- get_node(fline, "start")
end <- get_node(fline, "end")
plot(sf::st_zm(fline$geom),
lwd = fline$StreamOrde, col = "blue")
plot(sf::st_geometry(start), add = TRUE)
plot(sf::st_zm(fline$geom),
lwd = fline$StreamOrde, col = "blue")
plot(sf::st_geometry(end), add = TRUE)
Run the code above in your browser using DataLab