Given a list of outlets, get their basin boundaries and network and return a leaflet map in EPSG:4326.
map_nhdplus(
outlets = NULL,
bbox = NULL,
streamorder = NULL,
nhdplus_data = NULL,
gpkg = NULL,
flowline_only = NULL,
plot_config = NULL,
overwrite = TRUE,
cache_data = NULL,
return_map = FALSE
)
list of nldi outlets. Other inputs are coerced into nldi outlets, see details.
object of class bbox with a defined crs. See examples.
integer only streams of order greater than or equal will be returned
geopackage containing source nhdplus data (omit to download)
path and file with .gpkg ending. If omitted, no file is written.
boolean only subset and plot flowlines only, default=FALSE
list containing plot configuration, see details.
passed on the subset_nhdplus.
character path to rds file where all plot data can be cached. If file doesn't exist, it will be created. If set to FALSE, all caching will be turned off -- this includes basemap tiles.
if FALSE (default), a data.frame of plot data is returned invisibly in NAD83 Lat/Lon, if TRUE the leaflet object is returned
data.frame or leaflet map (see return_map)
map_nhdplus supports several input specifications. An unexported function "as_outlet" is used to convert the outlet formats as described below.
if outlets is omitted, the bbox input is required and all nhdplus data in the bounding box is plotted.
If outlets is a list of integers, it is assumed to be NHDPlus IDs (comids) and all upstream tributaries are plotted.
if outlets is an integer vector, it is assumed to be all NHDPlus IDs (comids) that should be plotted. Allows custom filtering.
If outlets is a character vector, it is assumed to be NWIS site ids.
if outlets is a list containing only characters, it is assumed to be a list of nldi features and all upstream tributaries are plotted.
if outlets is a data.frame with point geometry, a point in polygon match is performed and upstream with tributaries from the identified catchments is plotted.
See plot_nhdplus for details on plot configuration.
# \donttest{
map_nhdplus("05428500")
#> Spherical geometry (s2) switched off
#> although coordinates are longitude/latitude, st_intersects assumes that they
#> are planar
#> Spherical geometry (s2) switched on
map_nhdplus("05428500", streamorder = 2)
#> Spherical geometry (s2) switched off
#> although coordinates are longitude/latitude, st_intersects assumes that they
#> are planar
#> Spherical geometry (s2) switched on
map_nhdplus(list(13293970, 13293750))
#> Spherical geometry (s2) switched off
#> although coordinates are longitude/latitude, st_intersects assumes that they
#> are planar
#> Spherical geometry (s2) switched on
source(system.file("extdata/sample_data.R", package = "nhdplusTools"))
map_nhdplus(list(13293970, 13293750), streamorder = 3, nhdplus_data = sample_data)
#return leaflet object
map_nhdplus("05428500", return_map = TRUE)
# }