Subsets NHD waterbody features by location (POINT), area (POLYGON), or set of IDs.

get_waterbodies(AOI = NULL, id = NULL, t_srs = NULL, buffer = 0.5)



sf (MULTI)POINT or (MULTI)POLYGON. An 'area of interest' can be provided as either a location (sf POINT) or area (sf POLYGON) in any Spatial Reference System.


NHD Waterbody COMID(s)


character (PROJ string or EPSG code) or numeric (EPSG code). A user specified - target -Spatial Reference System (SRS/CRS) for returned objects. Will default to the CRS of the input AOI if provided, and to 4326 for ID requests.


numeric. The amount (in meters) to buffer a POINT AOI by for an extended search. Default = 0.5


a simple features (sf) object


The returned object(s) will have the same Spatial Reference System (SRS) as the input AOI. If a individual or set of IDs are used to query, then the default geoserver CRS of EPSG:4326 is preserved. In all cases, a user-defined SRS can be passed to t_srs which will override all previous SRS's (either input or default). All buffer and distance operations are handled internally using in EPSG:5070 Albers Equal Area projection