if (FALSE) {
coords <- data.frame(lat = 41.42217, lon = -73.24189)
network <- nhd_plus_query(lon = coords$lon, lat = coords$lat,
dsn = "NHDFlowline", buffer_dist = units::as_units(5, "km"))$sp$NHDFlowline
t_reaches <- tip_reaches(network = network)
plot(network$geometry)
plot(t_reaches$geometry, col = "red", add = TRUE)
}
Run the code above in your browser using DataLab