Get a single feature from the NLDI
get_nldi_feature(nldi_feature)list with names `featureSource` and `featureID` where `featureSource` is derived from the "source" column of the response of get_nldi_sources and the `featureID` is a known identifier from the specified `featureSource`.
sf data.frame with one feature
# \donttest{
get_nldi_feature(list("featureSource" = "nwissite", featureID = "USGS-05428500"))
#> Simple feature collection with 1 feature and 8 fields
#> Geometry type: POINT
#> Dimension: XY
#> Bounding box: xmin: -89.36083 ymin: 43.08944 xmax: -89.36083 ymax: 43.08944
#> Geodetic CRS: WGS 84
#> # A tibble: 1 × 9
#> sourceName identifier comid measure reachcode name X Y
#> <chr> <chr> <chr> <dbl> <chr> <chr> <dbl> <dbl>
#> 1 NWIS Surface Water Sites USGS-05428… 1329… 42.9 07090002… YAHA… -89.4 43.1
#> # ℹ 1 more variable: geometry <POINT [°]>
# }