The goal of this vignette is to demonstrate a simple and lightweight approach to building maps with NHDPlus data.
plot_nhdplus is a work in progress. Not all
inputs in the function have been implemented as of 11/18/2019 and
additional functionality will be added later. Please leave feature
requests and issues you find in an issue
plot_nhdplus is a function that makes getting a simple
plot of NHDPlus data as easy as possible. It works with other functions
nhdplusTools for identifying and retrieving watershed
outlet locations. See the
plot_nhdplus documentation for
If we pass
plot_nhdplus a single NWIS site id,
nhdplusTools uses web services to get data and we get a
plot like this:
If we want to add other watersheds, we can use any outlet available
from the Network Linked Data Index. See “nldi” functions elsewhere in
If we don’t know a site id, we can just pass in one or more latitude / longitude locations.
plot_nhdplus also allows modification of streamorder (if
you have data available locally) and styles. This plot request shows how
to get a subset of data for a plot and the range of options. See
documentation for more details.
source(system.file("extdata/sample_data.R", package = "nhdplusTools")) plot_nhdplus(list(list("comid", "13293970"), list("nwissite", "USGS-05428500"), list("huc12pp", "070900020603"), list("huc12pp", "070900020602")), streamorder = 2, nhdplus_data = sample_data, plot_config = list(basin = list(lwd = 2), outlets = list(huc12pp = list(cex = 1.5), comid = list(col = "green"))))
We can also plot NHDPlus data without an outlet at all.