Navigate upstream fails for comid 4651042

[hdugan]

Other NHD functions work for this comid. This fails, and so this site is not mapped to salt use.

comids_ut <- navigate_nldi(
  list(featureSource = "comid",
       featureID = 4651042),
  mode = "upstreamTributaries",
  distance_km = 1000)

[lindsayplatt]

Hmmm that command did not fail for me. I am using nhdplusTools v1.1.0.

> comids_ut 
$origin
Simple feature collection with 1 feature and 3 fields
Geometry type: LINESTRING
Dimension:     XY
Bounding box:  xmin: -75.78338 ymin: 39.82052 xmax: -75.76803 ymax: 39.82447
Geodetic CRS:  WGS 84
# A tibble: 1 × 4
  sourceName    identifier comid                                                                      geometry
  <chr>         <chr>      <chr>                                                              <LINESTRING [°]>
1 NHDPlus comid 4651042    4651042 (-75.76803 39.82382, -75.7682 39.82394, -75.76927 39.82405, -75.77001 39.8…

$UT_flowlines
Simple feature collection with 4 features and 1 field
Geometry type: LINESTRING
Dimension:     XY
Bounding box:  xmin: -75.78338 ymin: 39.82052 xmax: -75.76114 ymax: 39.82447
Geodetic CRS:  WGS 84
# A tibble: 4 × 2
  nhdplus_comid                                                                                  geometry
  <chr>                                                                                  <LINESTRING [°]>
1 4651042       (-75.76803 39.82382, -75.7682 39.82394, -75.76927 39.82405, -75.77001 39.82433, -75.77...
2 4652086                (-75.76598 39.82362, -75.76606 39.82352, -75.76784 39.82379, -75.76803 39.82382)
3 4651036                 (-75.76479 39.8235, -75.76491 39.82373, -75.76536 39.82373, -75.76598 39.82362)
4 4652354                                                         (-75.76114 39.82251, -75.76479 39.8235)