Multipurpose function to find a COMID of interest.
discover_nhdplus_id(point = NULL, nldi_feature = NULL, raindrop = FALSE)sfc POINT including crs as created by:
sf::st_sfc(sf::st_point(.. ,..), crs)
list with names `featureSource` and `featureID` where `featureSource` is derived from the "source" column of the response of get_nldi_sources and the `featureSource` is a known identifier from the specified `featureSource`.
logical if TRUE will call a raindrop trace web service and
return will be the same as get_raindrop_trace with direction "none".
integer COMID or list containing COMID and raindrop trace.
# \donttest{
point <- sf::st_sfc(sf::st_point(c(-76.874, 39.482)), crs = 4326)
discover_nhdplus_id(point)
#> [1] 11689978
discover_nhdplus_id(point, raindrop = TRUE)
#> Simple feature collection with 2 features and 7 fields
#> Geometry type: LINESTRING
#> Dimension: XY
#> Bounding box: xmin: -76.88292 ymin: 39.47138 xmax: -76.86733 ymax: 39.49148
#> Geodetic CRS: WGS 84
#> # A tibble: 2 × 8
#> id gnis_name comid reachcode raindrop_pathDist measure intersection_point
#> <chr> <chr> <int> <chr> <dbl> <dbl> <list>
#> 1 nhdF… North Br… 1.17e7 02060003… 62.6 53.9 <dbl [2]>
#> 2 rain… NA NA NA NA NA <dbl [0]>
#> # ℹ 1 more variable: geometry <LINESTRING [°]>
nldi_nwis <- list(featureSource = "nwissite", featureID = "USGS-08279500")
discover_nhdplus_id(nldi_feature = nldi_nwis)
#> [1] 17864756
# }