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 9 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 × 10
#> sourceName identifier comid name reachcode measure mainstem X Y
#> <chr> <chr> <int> <chr> <chr> <dbl> <chr> <dbl> <dbl>
#> 1 NWIS Surface W… USGS-0542… 1.33e7 YAHA… 07090002… 42.0 https:/… -89.4 43.1
#> # ℹ 1 more variable: geometry <POINT [°]>
# }