# U.S. Data and Package Basics

#### dblodgett@usgs.gov

## TL;DR

Pick an outlet location and download some data.

    # Uncomment to install!
    # install.packages("nhdplusTools")

    library(nhdplusTools)
    library(sf)

    start_point <- st_sfc(st_point(c(-89.362239, 43.090266)), crs = 4269)
    start_comid <- discover_nhdplus_id(start_point)

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

    subset_file <- tempfile(fileext = ".gpkg")
    subset <- subset_nhdplus(comids = as.integer(flowline$UT$nhdplus_comid),
                             output_file = subset_file,
                             nhdplus_data = "download",
                             flowline_only = FALSE,
                             return_data = TRUE, overwrite = TRUE)
    #> All intersections performed in latitude/longitude.
    #> Reading NHDFlowline_Network
    #> Spherical geometry (s2) switched off
    #> https://api.water.usgs.gov/geoserver/wmadata/ows
    #> <?xml version="1.0" encoding="UTF-8"?>
    #> <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" service="WFS" version="1.1.0" outputFormat="application/json" xsi:schemaLocation="http://www.opengis.net/wfs                         http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    #>   <wfs:Query xmlns:feature="https://api.water.usgs.gov/wmadata" typeName="feature:nhdflowline_network" srsName="EPSG:4269">
    #>     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
    #>       <ogc:Or>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293376</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293378</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293380</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293382</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293384</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293386</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293388</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293390</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293392</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293394</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293396</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293398</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293400</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293404</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293406</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293410</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293412</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293416</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293418</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293420</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293422</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293424</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293426</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293428</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293430</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293432</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293434</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293436</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293438</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293440</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293442</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293444</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293446</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293448</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293450</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293452</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293454</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293456</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293458</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293460</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293462</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293464</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293468</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293472</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293474</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293476</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293478</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293480</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293484</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293486</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293488</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293490</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293492</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293494</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293496</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293498</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293502</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293504</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293506</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293508</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293512</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293514</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293516</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293518</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293520</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293522</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293524</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293526</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293530</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293532</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293536</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293540</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293542</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293548</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293550</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293552</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293554</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293556</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293558</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293562</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293564</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293568</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293570</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293572</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293574</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293576</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293578</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293584</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293586</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293588</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293590</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293600</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293612</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293614</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293620</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293624</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293634</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293646</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293672</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293674</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293676</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293678</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293688</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293690</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293692</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293694</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293696</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293750</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294110</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294112</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294114</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294116</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294118</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294120</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294122</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294124</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294128</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294130</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294132</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294134</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294138</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294140</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294144</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294148</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294150</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294152</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294154</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294158</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294166</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294168</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294174</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294176</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294178</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294184</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294264</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294268</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294270</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294272</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294274</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294276</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294278</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294280</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294282</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294284</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294286</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294288</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294290</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294292</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294294</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294296</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294298</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294300</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294304</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294306</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294308</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294310</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294312</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294314</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294366</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294368</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294382</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294384</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294386</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294394</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13302588</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13302590</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>937090090</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>937090091</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>       </ogc:Or>
    #>     </ogc:Filter>
    #>   </wfs:Query>
    #> </wfs:GetFeature>
    #> Spherical geometry (s2) switched on
    #> Writing NHDFlowline_Network
    #> Reading CatchmentSP
    #> Spherical geometry (s2) switched off
    #> https://api.water.usgs.gov/geoserver/wmadata/ows
    #> <?xml version="1.0" encoding="UTF-8"?>
    #> <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" service="WFS" version="1.1.0" outputFormat="application/json" xsi:schemaLocation="http://www.opengis.net/wfs                         http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    #>   <wfs:Query xmlns:feature="https://api.water.usgs.gov/wmadata" typeName="feature:catchmentsp" srsName="EPSG:4269">
    #>     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
    #>       <ogc:Or>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293376</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293378</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293380</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293382</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293384</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293386</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293388</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293390</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293392</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293394</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293396</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293398</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293400</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293404</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293406</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293410</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293412</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293416</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293418</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293420</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293422</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293424</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293426</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293428</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293430</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293432</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293434</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293436</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293438</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293440</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293442</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293444</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293446</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293448</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293450</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293452</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293454</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293456</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293458</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293460</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293462</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293464</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293468</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293472</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293474</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293476</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293478</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293480</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293484</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293486</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293488</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293490</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293492</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293494</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293496</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293498</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293502</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293504</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293506</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293508</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293512</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293514</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293516</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293518</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293520</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293522</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293524</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293526</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293530</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293532</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293536</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293540</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293542</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293548</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293550</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293552</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293554</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293556</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293558</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293562</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293564</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293568</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293570</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293572</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293574</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293576</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293578</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293584</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293586</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293588</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293590</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293600</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293612</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293614</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293620</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293624</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293634</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293646</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293672</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293674</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293676</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293678</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293688</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293690</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293692</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293694</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293696</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13293750</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294110</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294112</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294114</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294116</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294118</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294120</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294122</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294124</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294128</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294130</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294132</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294134</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294138</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294140</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294144</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294148</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294150</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294152</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294154</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294158</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294166</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294168</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294174</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294176</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294178</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294184</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294264</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294268</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294270</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294272</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294274</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294276</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294278</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294280</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294282</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294284</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294286</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294288</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294290</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294292</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294294</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294296</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294298</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294300</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294304</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294306</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294308</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294310</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294312</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294314</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294366</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294368</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294382</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294384</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294386</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13294394</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13302588</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>13302590</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>937090090</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>featureid</ogc:PropertyName>
    #>           <ogc:Literal>937090091</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>       </ogc:Or>
    #>     </ogc:Filter>
    #>   </wfs:Query>
    #> </wfs:GetFeature>
    #> Spherical geometry (s2) switched on
    #> Writing CatchmentSP
    #> Spherical geometry (s2) switched off
    #> https://api.water.usgs.gov/geoserver/wmadata/ows
    #> <?xml version="1.0" encoding="UTF-8"?>
    #> <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" service="WFS" version="1.1.0" outputFormat="application/json" xsi:schemaLocation="http://www.opengis.net/wfs                         http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    #>   <wfs:Query xmlns:feature="https://api.water.usgs.gov/wmadata" typeName="feature:nhdarea" srsName="EPSG:4269">
    #>     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
    #>       <ogc:And>
    #>         <ogc:BBOX>
    #>           <ogc:PropertyName>the_geom</ogc:PropertyName>
    #>           <gml:Envelope srsName="urn:x-ogc:def:crs:EPSG:4326">
    #>             <gml:lowerCorner>43.035066 -89.604794</gml:lowerCorner>
    #>             <gml:upperCorner>43.366068 -89.203781</gml:upperCorner>
    #>           </gml:Envelope>
    #>         </ogc:BBOX>
    #>       </ogc:And>
    #>     </ogc:Filter>
    #>   </wfs:Query>
    #> </wfs:GetFeature>
    #> although coordinates are longitude/latitude, st_intersects assumes that they
    #> are planar
    #> Spherical geometry (s2) switched on
    #> Spherical geometry (s2) switched off
    #> https://api.water.usgs.gov/geoserver/wmadata/ows
    #> <?xml version="1.0" encoding="UTF-8"?>
    #> <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" service="WFS" version="1.1.0" outputFormat="application/json" xsi:schemaLocation="http://www.opengis.net/wfs                         http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    #>   <wfs:Query xmlns:feature="https://api.water.usgs.gov/wmadata" typeName="feature:nhdwaterbody" srsName="EPSG:4269">
    #>     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
    #>       <ogc:And>
    #>         <ogc:BBOX>
    #>           <ogc:PropertyName>the_geom</ogc:PropertyName>
    #>           <gml:Envelope srsName="urn:x-ogc:def:crs:EPSG:4326">
    #>             <gml:lowerCorner>43.035066 -89.604794</gml:lowerCorner>
    #>             <gml:upperCorner>43.366068 -89.203781</gml:upperCorner>
    #>           </gml:Envelope>
    #>         </ogc:BBOX>
    #>       </ogc:And>
    #>     </ogc:Filter>
    #>   </wfs:Query>
    #> </wfs:GetFeature>
    #> although coordinates are longitude/latitude, st_intersects assumes that they
    #> are planar
    #> Spherical geometry (s2) switched on
    #> Spherical geometry (s2) switched off
    #> https://api.water.usgs.gov/geoserver/wmadata/ows
    #> <?xml version="1.0" encoding="UTF-8"?>
    #> <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" service="WFS" version="1.1.0" outputFormat="application/json" xsi:schemaLocation="http://www.opengis.net/wfs                         http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    #>   <wfs:Query xmlns:feature="https://api.water.usgs.gov/wmadata" typeName="feature:nhdflowline_nonnetwork" srsName="EPSG:4269">
    #>     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
    #>       <ogc:And>
    #>         <ogc:BBOX>
    #>           <ogc:PropertyName>the_geom</ogc:PropertyName>
    #>           <gml:Envelope srsName="urn:x-ogc:def:crs:EPSG:4326">
    #>             <gml:lowerCorner>43.035066 -89.604794</gml:lowerCorner>
    #>             <gml:upperCorner>43.366068 -89.203781</gml:upperCorner>
    #>           </gml:Envelope>
    #>         </ogc:BBOX>
    #>       </ogc:And>
    #>     </ogc:Filter>
    #>   </wfs:Query>
    #> </wfs:GetFeature>
    #> although coordinates are longitude/latitude, st_intersects assumes that they
    #> are planar
    #> Spherical geometry (s2) switched on

    flowline <- subset$NHDFlowline_Network
    catchment <- subset$CatchmentSP
    waterbody <- subset$NHDWaterbody

    ## Or using a file:

    flowline <- sf::read_sf(subset_file, "NHDFlowline_Network")
    catchment <- sf::read_sf(subset_file, "CatchmentSP")
    waterbody <- sf::read_sf(subset_file, "NHDWaterbody")

    plot(sf::st_geometry(flowline), col = "blue")
    plot(start_point, cex = 1.5, lwd = 2, col = "red", add = TRUE)
    plot(sf::st_geometry(catchment), add = TRUE)
    plot(sf::st_geometry(waterbody), col = rgb(0, 0, 1, alpha = 0.5), add = TRUE)

Or fetch NWIS an site as the starting point and generate a plot. Data is returned and/or stored in a local file for later use.

    # ?plot_nhdplus for more
    plot_data <- plot_nhdplus(
      outlets = list(featureSource = "nwissite", featureID = "USGS-05428500"),
      gpkg = subset_file, overwrite = TRUE)
    #> https://labs-beta.waterdata.usgs.gov/api/nldi/linked-data/nwissite/USGS-05428500/basin?simplify=true&splitCatchment=false
    #> Spherical geometry (s2) switched off
    #> https://api.water.usgs.gov/geoserver/wmadata/ows
    #> <?xml version="1.0" encoding="UTF-8"?>
    #> <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" service="WFS" version="1.1.0" outputFormat="application/json" xsi:schemaLocation="http://www.opengis.net/wfs                         http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    #>   <wfs:Query xmlns:feature="https://api.water.usgs.gov/wmadata" typeName="feature:nhdflowline_network" srsName="EPSG:4269">
    #>     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
    #>       <ogc:And>
    #>         <ogc:PropertyIsGreaterThan>
    #>           <ogc:PropertyName>streamorde</ogc:PropertyName>
    #>           <ogc:Literal>2</ogc:Literal>
    #>         </ogc:PropertyIsGreaterThan>
    #>         <ogc:BBOX>
    #>           <ogc:PropertyName>the_geom</ogc:PropertyName>
    #>           <gml:Envelope srsName="urn:x-ogc:def:crs:EPSG:4326">
    #>             <gml:lowerCorner>43.035066012 -89.604654655</gml:lowerCorner>
    #>             <gml:upperCorner>43.366068182 -89.203780859</gml:upperCorner>
    #>           </gml:Envelope>
    #>         </ogc:BBOX>
    #>       </ogc:And>
    #>     </ogc:Filter>
    #>   </wfs:Query>
    #> </wfs:GetFeature>
    #> although coordinates are longitude/latitude, st_intersects assumes that they
    #> are planar
    #> Spherical geometry (s2) switched on
    #> Zoom set to: 11

This vignette covers a range of utilities the nhdplusTools package offers for working with data in a U.S. context.

The first thing you are going to need to do is go get some data to work with. `nhdplusTools` provides the ability to download small subsets of the NHDPlus directly from web services. For large subsets, greater than a few thousand square kilometers, you can use `download_nhdplusv2()`.

If you are working with the whole National Seamless database, `nhdplusTools` has some convenience functions you should be aware of. Once you have it downloaded and extracted, you can tell the nhdplusTools package where it is with the `nhdplus_path()` function.

    nhdplus_path(file.path(work_dir, "natseamless.gpkg"))

    basename(nhdplus_path())
    #> [1] "natseamless.gpkg"

If you don’t need or want all the geometry for the flowlines, consider using `get_vaa()`. It allows you to retrieve specific NHDPlus attributes with very little overhead. It also supports access to an updated set of network attributes that incorporate numerous network updates from national hydrologic modeling projects.

    vaa <- get_vaa()
    names(vaa)
    #>  [1] "comid"      "streamleve" "streamorde" "streamcalc" "fromnode"
    #>  [6] "tonode"     "hydroseq"   "levelpathi" "pathlength" "terminalpa"
    #> [11] "arbolatesu" "divergence" "startflag"  "terminalfl" "dnlevel"
    #> [16] "thinnercod" "uplevelpat" "uphydroseq" "dnlevelpat" "dnminorhyd"
    #> [21] "dndraincou" "dnhydroseq" "frommeas"   "tomeas"     "reachcode"
    #> [26] "lengthkm"   "fcode"      "vpuin"      "vpuout"     "areasqkm"
    #> [31] "totdasqkm"  "divdasqkm"  "totma"      "wbareatype" "pathtimema"
    #> [36] "slope"      "slopelenkm" "ftype"      "gnis_name"  "gnis_id"
    #> [41] "wbareacomi" "hwnodesqkm" "rpuid"      "vpuid"      "roughness"
    nrow(vaa)
    #> [1] 2691339

### NHDPlus HiRes

NHDPlus HiRes is an in-development dataset that introduces much more dense flowlines and catchments. In the long run, `nhdplusTools` will have complete support for NHDPlus HiRes. So far, `nhdplusTools` will help download and interface NHDPlus HiRes data with existing `nhdplusTools` functionality. It’s important to note that `nhdplusTools` was primarily implemented using NHDPlusV2 and any use of HiRes should be subject to scrutiny.

For the demo below, a small sample of HiRes data that has been loaded into `nhdplusTools` is used. The first line shows how you can download additional data (just change `download_files` to `TRUE`).

    download_nhdplushr(nhd_dir = "download_dir",
                       hu_list = c("0101"), # can mix hu02 and hu04 codes.
                       download_files = FALSE) # TRUE will download files.
    #> [1] "https://prd-tnm.s3.amazonaws.com/StagedProducts/Hydrography/NHDPlusHR/VPU/Current/GDB/NHDPLUS_H_0101_HU4_20220901_GDB.zip"

    out_gpkg <- file.path(work_dir, "nhd_hr.gpkg")
    hr_data <- get_nhdplushr(work_dir,
                             out_gpkg = out_gpkg)
    (layers <- st_layers(out_gpkg))
    #> Driver: GPKG
    #> Available layers:
    #>         layer_name geometry_type features fields             crs_name
    #> 1      NHDFlowline   Line String     2691     57 GRS 1980(IUGG, 1980)
    #> 2 NHDPlusCatchment Multi Polygon     2603      7 GRS 1980(IUGG, 1980)
    #> 3     NHDWaterbody       Polygon     1044     15 GRS 1980(IUGG, 1980)
    #> 4          NHDArea       Polygon       10     14 GRS 1980(IUGG, 1980)
    #> 5          NHDLine   Line String      142     12 GRS 1980(IUGG, 1980)
    #> 6      NHDPlusSink         Point        1     10 GRS 1980(IUGG, 1980)
    #> 7         NHDPoint      3D Point        7     10 GRS 1980(IUGG, 1980)
    names(hr_data)
    #> [1] "NHDFlowline"      "NHDPlusCatchment"
    unlink(out_gpkg)

    hr_data <- get_nhdplushr(work_dir,
                             out_gpkg = out_gpkg,
                             layers = NULL)
    (layers <- st_layers(out_gpkg))
    #> Driver: GPKG
    #> Available layers:
    #>         layer_name geometry_type features fields             crs_name
    #> 1      NHDFlowline   Line String     2691     57 GRS 1980(IUGG, 1980)
    #> 2 NHDPlusCatchment Multi Polygon     2603      7 GRS 1980(IUGG, 1980)
    #> 3     NHDWaterbody       Polygon     1044     15 GRS 1980(IUGG, 1980)
    #> 4          NHDArea       Polygon       10     14 GRS 1980(IUGG, 1980)
    #> 5          NHDLine   Line String      142     12 GRS 1980(IUGG, 1980)
    #> 6      NHDPlusSink         Point        1     10 GRS 1980(IUGG, 1980)
    #> 7         NHDPoint      3D Point        7     10 GRS 1980(IUGG, 1980)
    names(hr_data)
    #> [1] "NHDFlowline"      "NHDPlusCatchment" "NHDWaterbody"     "NHDArea"
    #> [5] "NHDLine"          "NHDPlusSink"      "NHDPoint"

## Discovery and Subsetting

One of the primary workflows `nhdplusTools` is designed to accomplish can be described in three steps:

1. what NHDPlus catchment is at the outlet of a watershed,
2. figure out what catchments are up or downstream of that catchment, and
3. create a stand alone subset for that collection of catchments.

Say we want to get a subset of the NHDPlus upstream of a given location. We can start with `discover_nhdplus_id()` First, let’s look at a given point location. Then see where it is relative to our flowlines.

    lon <- -89.36
    lat <- 43.09

    start_point <- sf::st_sfc(sf::st_point(c(lon, lat)),
                              crs = 4269)

    plot(sf::st_geometry(flowline))
    plot(start_point, cex = 1.5, lwd = 2, col = "red", add = TRUE)

OK, so we have a point location near a river and we want to figure out what catchment is at its outlet. We can use the `discover_nhdplus_id()` function which calls out to a web service and returns an NHDPlus catchment identifier, typically called a COMID.

If we set `raindrop = TRUE`, we can also get a elevation derived downslope trace to the nearest flowline and some additional info. See `get_raindrop_trace()` for more on this functionality.

    start_comid <- discover_nhdplus_id(start_point, raindrop = TRUE)
    #> https://labs-beta.waterdata.usgs.gov/api/nldi/pygeoapi/processes/nldi-flowtrace/execution
    #> {
    #>   "inputs": [
    #>     {
    #>       "id": "lat",
    #>       "type": "text/plain",
    #>       "value": "43.09"
    #>     },
    #>     {
    #>       "id": "lon",
    #>       "type": "text/plain",
    #>       "value": "-89.36"
    #>     },
    #>     {
    #>       "id": "raindroptrace",
    #>       "type": "text/plain",
    #>       "value": "true"
    #>     },
    #>     {
    #>       "id": "direction",
    #>       "type": "text/plain",
    #>       "value": "none"
    #>     }
    #>   ]
    #> }
    start_comid
    #> Simple feature collection with 2 features and 7 fields
    #> Geometry type: LINESTRING
    #> Dimension:     XY
    #> Bounding box:  xmin: -89.37037 ymin: 43.08522 xmax: -89.35393 ymax: 43.09491
    #> 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… Yahara R…  1.33e7 07090002…              91.8    42.4 <dbl [2]>
    #> 2 rain… <NA>      NA      <NA>                   NA      NA   <dbl [0]>
    #> # ℹ 1 more variable: geometry <LINESTRING [°]>

    plot(sf::st_geometry(start_comid))
    plot(sf::st_geometry(flowline), add = TRUE, col = "blue", lwd = 2)
    plot(start_point, cex = 1.5, lwd = 2, col = "red", add = TRUE)

If you have the whole National Seamless database and want to work at regional to national scales, skip down the the Local Data Subsetting section.

### Web Service Data Subsetting

`nhdplusTools` supports discovery and data subsetting using web services made available through the Network Linked Data Index (NLDI) and the Water Mission Area Geoserver. The code below shows how to use the NLDI functions to build a dataset upstream of our `start_comid` that we found above.

The NLDI can be queried with any set of watershed outlet locations that it has in its index. We call these “featureSources”. We can query the NLDI for an identifier of a given feature from any of its “featureSources” and find out what our navigation options are as shown below.

    dataRetrieval::get_nldi_sources()$source
    #>  [1] "comid"             "HILARRI"           "GRAND"
    #>  [4] "wade_timeseries"   "wade_rights"       "ref_dams"
    #>  [7] "huc12pp_102020"    "npdes"             "census2020-nhdpv2"
    #> [10] "epa_nrsa"          "geoconnex-demo"    "nmwdi-st"
    #> [13] "nwisgw"            "vigil"             "gfv11_pois"
    #> [16] "ca_gages"          "ref_gage"          "wade"
    #> [19] "nwissite"          "huc12pp"           "WQP"

    nldi_feature <- list(featureSource = "comid",
                         featureID = as.integer(start_comid$comid)[1])

    get_nldi_feature(nldi_feature)
    #> Simple feature collection with 1 feature and 4 fields
    #> Geometry type: LINESTRING
    #> Dimension:     XYZ
    #> Bounding box:  xmin: -89.37037 ymin: 43.08521 xmax: -89.35393 ymax: 43.09491
    #> z_range:       zmin: 0 zmax: 100
    #> Geodetic CRS:  WGS 84
    #> # A tibble: 1 × 5
    #>   sourceName    identifier    comid reachcode                           geometry
    #>   <chr>         <chr>         <int> <chr>                       <LINESTRING [°]>
    #> 1 NHDPlus comid 13293750   13293750 07090002007373 Z (-89.37037 43.09491 100, -…

We can use `get_nldi_feature()` as a way to make sure the featureID is available for the chosen “featureSource”. Now that we know the NLDI has our comid, we can use the “upstreamTributaries” navigation option to get all the flowlines upstream or all the features from any of the “featureSources” as shown below.

    flowline_nldi <- navigate_nldi(nldi_feature,
                                   mode = "upstreamTributaries",
                                   distance_km = 1000)

    plot(sf::st_geometry(flowline), lwd = 3, col = "black")
    plot(sf::st_geometry(flowline_nldi$origin), lwd = 3, col = "red", add = TRUE)
    plot(sf::st_geometry(flowline_nldi$UT), lwd = 1, col = "red", add = TRUE)

The NLDI only provides geometry and a comid for each of the flowlines. The `subset_nhdplus()` function has a “download” option that allows us to download four layers and all attributes as shown below. There is also a `navigate_network()` function that will replace `navigate_nldi()` and `subset_nhdplus()` for many use cases.

    output_file_download <- file.path(work_dir, "subset_download.gpkg")

    output_file_download <-subset_nhdplus(comids = as.integer(flowline_nldi$UT$nhdplus_comid),
                                          output_file = output_file_download,
                                          nhdplus_data = "download", return_data = FALSE,
                                          overwrite = TRUE)
    #> All intersections performed in latitude/longitude.
    #> Reading NHDFlowline_Network
    #> Spherical geometry (s2) switched off
    #> https://api.water.usgs.gov/geoserver/wmadata/ows
    #> <?xml version="1.0" encoding="UTF-8"?>
    #> <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" service="WFS" version="1.1.0" outputFormat="application/json" xsi:schemaLocation="http://www.opengis.net/wfs                         http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    #>   <wfs:Query xmlns:feature="https://api.water.usgs.gov/wmadata" typeName="feature:nhdflowline_network" srsName="EPSG:4269">
    #>     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
    #>       <ogc:Or>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293376</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293378</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293380</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293382</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293384</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293386</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293388</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293390</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293392</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293394</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293396</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293398</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293400</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293404</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293406</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293410</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293412</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293416</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293418</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293420</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293422</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293424</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293426</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293428</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293430</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293432</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293434</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293436</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293438</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293440</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293442</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293444</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293446</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293448</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293450</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293452</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293454</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293456</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293458</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293460</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293462</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293464</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293468</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293472</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293474</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293476</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293478</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293480</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293484</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293486</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293488</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293490</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293492</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293494</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293496</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293498</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293502</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293504</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293506</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293508</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293512</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293514</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293516</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293518</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293520</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293522</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293524</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293526</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293530</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293532</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293536</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293540</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293542</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293548</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293550</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293552</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293554</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293556</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293558</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293562</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293564</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293568</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293570</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293572</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293574</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293576</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293578</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293584</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293586</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293588</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293590</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293600</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293612</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293614</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293620</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293624</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293634</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293646</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293672</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293674</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293676</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293678</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293688</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293690</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293692</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293694</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293696</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293750</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294110</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294112</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294114</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294116</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294118</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294120</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294122</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294124</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294128</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294130</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294132</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294134</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294138</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294140</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294144</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294148</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294150</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294152</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294154</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294158</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294166</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294168</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294174</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294176</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294178</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294184</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294264</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294268</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294270</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294272</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294274</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294276</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294278</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294280</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294282</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294284</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294286</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294288</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294290</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294292</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294294</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294296</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294298</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294300</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294304</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294306</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294308</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294310</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294312</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294314</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294366</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294368</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294382</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294384</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294386</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294394</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13302588</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13302590</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>937090090</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>937090091</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>       </ogc:Or>
    #>     </ogc:Filter>
    #>   </wfs:Query>
    #> </wfs:GetFeature>
    #> Spherical geometry (s2) switched on
    #> Writing NHDFlowline_Network

    sf::st_layers(output_file_download)
    #> Driver: GPKG
    #> Available layers:
    #>            layer_name geometry_type features fields crs_name
    #> 1 NHDFlowline_Network   Line String      168    137    NAD83

    flowline_download <- sf::read_sf(file.path(work_dir, "subset_download.gpkg"),
                                     "NHDFlowline_Network")

    plot(sf::st_geometry(dplyr::filter(flowline_download,
                                       streamorde > 2)),
         lwd = 7, col = "darkgrey")
    plot(sf::st_geometry(flowline_nldi$UT),
         lwd = 3, col = "red", add = TRUE)

This plot illustrates the kind of thing that’s possible (filtering to specific stream orders) using the attributes that are downloaded.

Before moving on, one more demonstration of what can be done using the NLDI. Say we knew the USGS gage ID that we want NHDPlus data upstream of. We can use the NLDI to navigate from the gage the same as we did for our comid. We can also get back all the nwis sites the NLDI knows about upstream of the one we chose!

    nldi_feature <- list(featureSource = "nwissite",
                         featureID = "USGS-05428500")

    flowline_nldi <- navigate_nldi(nldi_feature,
                                   mode = "upstreamTributaries",
                                   distance_km = 1000)

    output_file_nwis <- file.path(work_dir, "subset_download_nwis.gpkg")

    output_file_nwis <-subset_nhdplus(comids = as.integer(flowline_nldi$UT$nhdplus_comid),
                                      output_file = output_file_nwis,
                                      nhdplus_data = "download",
                                      return_data = FALSE, overwrite = TRUE)
    #> All intersections performed in latitude/longitude.
    #> Reading NHDFlowline_Network
    #> Spherical geometry (s2) switched off
    #> https://api.water.usgs.gov/geoserver/wmadata/ows
    #> <?xml version="1.0" encoding="UTF-8"?>
    #> <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" service="WFS" version="1.1.0" outputFormat="application/json" xsi:schemaLocation="http://www.opengis.net/wfs                         http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    #>   <wfs:Query xmlns:feature="https://api.water.usgs.gov/wmadata" typeName="feature:nhdflowline_network" srsName="EPSG:4269">
    #>     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">
    #>       <ogc:Or>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293376</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293378</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293380</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293382</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293384</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293386</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293388</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293390</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293392</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293394</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293396</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293398</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293400</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293404</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293406</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293410</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293412</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293416</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293418</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293420</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293422</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293424</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293426</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293428</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293430</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293432</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293434</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293436</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293438</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293440</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293442</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293444</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293446</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293448</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293450</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293452</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293454</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293456</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293458</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293460</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293462</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293464</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293468</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293472</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293474</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293476</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293478</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293480</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293484</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293486</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293488</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293490</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293492</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293494</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293496</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293498</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293502</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293504</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293506</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293508</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293512</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293514</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293516</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293518</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293520</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293522</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293524</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293526</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293530</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293532</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293536</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293540</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293542</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293548</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293550</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293552</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293554</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293556</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293558</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293562</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293564</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293568</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293570</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293572</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293574</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293576</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293578</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293584</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293586</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293588</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293590</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293600</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293612</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293614</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293620</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293624</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293634</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293646</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293672</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293674</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293676</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293678</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293688</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293690</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293692</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293694</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293696</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13293750</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294110</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294112</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294114</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294116</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294118</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294120</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294122</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294124</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294128</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294130</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294132</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294134</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294138</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294140</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294144</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294148</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294150</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294152</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294154</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294158</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294166</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294168</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294174</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294176</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294178</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294184</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294264</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294268</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294270</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294272</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294274</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294276</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294278</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294280</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294282</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294284</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294286</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294288</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294290</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294292</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294294</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294296</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294298</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294300</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294304</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294306</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294308</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294310</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294312</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294314</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294366</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294368</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294382</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294384</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294386</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13294394</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13302588</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>13302590</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>937090090</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>         <ogc:PropertyIsEqualTo>
    #>           <ogc:PropertyName>comid</ogc:PropertyName>
    #>           <ogc:Literal>937090091</ogc:Literal>
    #>         </ogc:PropertyIsEqualTo>
    #>       </ogc:Or>
    #>     </ogc:Filter>
    #>   </wfs:Query>
    #> </wfs:GetFeature>
    #> Spherical geometry (s2) switched on
    #> Writing NHDFlowline_Network

    sf::st_layers(output_file_download)
    #> Driver: GPKG
    #> Available layers:
    #>            layer_name geometry_type features fields crs_name
    #> 1 NHDFlowline_Network   Line String      168    137    NAD83

    flowline_nwis <- sf::read_sf(output_file_nwis,
                                     "NHDFlowline_Network")

    upstream_nwis <- navigate_nldi(nldi_feature,
                                   mode = "upstreamTributaries",
                                   data_source = "nwissite",
                                   distance_km = 1000)

    plot(sf::st_geometry(flowline_nwis),
         lwd = 3, col = "blue")
    plot(sf::st_geometry(upstream_nwis$UT_nwissite),
         cex = 1, lwd = 2, col = "red", add = TRUE)

### Local Data Subsetting

While web service data access is very convenient, some use cases make working with web services impossible or cumbersome such that working with local data is preferable. nhdplusTools supports such workflows with hybrid, web-service and local, workflows.

With the starting COMID we found with `discover_nhdplus_id()` above, we can use one of the network navigation functions, `get_UM()`, `get_UT()`, `get_DM()`, or `get_DD()` to retrieve a collection of comids along the upstream mainstem, upstream with tributaries, downstream mainstem, or downstream with diversions network paths. Here we’ll use upstream with tributaries.

    UT_comids <- get_UT(flowline, start_comid$comid[1])
    UT_comids
    #>   [1]  13293380  13293384  13293388  13293392  13293398  13293424  13293430
    #>   [8]  13293454  13293504  13293506  13293576  13293750  13294110  13294128
    #>  [15]  13294134  13294280  13294284  13294288  13294290  13294298  13294304
    #>  [22]  13294310  13294312  13294394  13293376  13293382  13293386  13293390
    #>  [29]  13293394  13293396  13293400  13293404  13293406  13293416  13293422
    #>  [36]  13293426  13293432  13293434  13293440  13293444  13293446  13293452
    #>  [43]  13293456  13293458  13293464  13293474  13293480  13293484  13293488
    #>  [50]  13293492  13293494  13293496  13293508  13293520  13293524  13293526
    #>  [57]  13293530  13293536  13293540  13293542  13293550  13293554  13293574
    #>  [64]  13293584  13293588  13293612  13293620  13293634  13293676  13293694
    #>  [71]  13293696  13294118  13294124  13294130  13294132  13294138  13294154
    #>  [78]  13294158  13294166  13294264  13294268  13294272  13294274  13294276
    #>  [85]  13294278  13294282  13294286  13294292  13294294  13294300  13294308
    #>  [92]  13294314  13294366  13294368  13294382  13294384  13294386 937090090
    #>  [99] 937090091  13293378  13293412  13293420  13293428  13293436  13293438
    #> [106]  13293442  13293450  13293460  13293462  13293468  13293472  13293476
    #> [113]  13293478  13293486  13293490  13293498  13293502  13293512  13293514
    #> [120]  13293516  13293518  13293522  13293532  13293548  13293552  13293556
    #> [127]  13293558  13293564  13293568  13293572  13293578  13293586  13293614
    #> [134]  13293624  13293672  13293674  13293678  13293690  13293692  13294112
    #> [141]  13294122  13294140  13294144  13294148  13294150  13294152  13294168
    #> [148]  13294176  13294184  13294270  13294296  13294306  13302588  13302590
    #> [155]  13293410  13293418  13293448  13293562  13293570  13293600  13293646
    #> [162]  13293688  13294114  13294116  13294120  13294174  13294178  13293590

If you are familiar with the NHDPlus, you will recognize that now that we have this list of COMIDs, we could go off and do all sorts of things with the various flowline attributes. For now, let’s just use the COMID list to filter our `fline` `sf` `data.frame` and plot it with our other layers.

    plot(sf::st_geometry(flowline))
    plot(start_point, cex = 1.5, lwd = 2, col = "red", add = TRUE)
    plot(sf::st_geometry(dplyr::filter(flowline, comid %in% UT_comids)),
         add=TRUE, col = "red", lwd = 2)

Say you want to save the network subset for later use in R or in some other GIS. The `subset_nhdplus()` function is your friend. If you have the whole national seamless database downloaded, you can pull out large subsets of it like shown below (this queries for data from the local geodatabase without loading the whole thing into memory). If you don’t have the whole national seamless, look at the second example in this section.

    output_file <- file.path(work_dir, "subset.gpkg")

    output_file <-subset_nhdplus(comids = UT_comids,
                                 output_file = output_file,
                                 nhdplus_data = nhdplus_path(),
                                 return_data = FALSE, overwrite = TRUE)
    #> All intersections performed in latitude/longitude.
    #> Reading NHDFlowline_Network
    #> 168 comids of 168
    #> Writing NHDFlowline_Network
    #> Reading CatchmentSP
    #> 168 comids of 168
    #> Writing CatchmentSP
    #> Reading NHDArea
    #> Writing NHDArea
    #> Reading NHDWaterbody
    #> Writing NHDWaterbody
    #> Reading NHDFlowline_NonNetwork
    #> Writing NHDFlowline_NonNetwork
    #> Reading Gage
    #> Writing Gage
    #> Reading Sink
    #> No features to write in Sink

    sf::st_layers(output_file)
    #> Driver: GPKG
    #> Available layers:
    #>               layer_name geometry_type features fields             crs_name
    #> 1    NHDFlowline_Network   Line String      168    136 GRS 1980(IUGG, 1980)
    #> 2            CatchmentSP Multi Polygon      167      6 GRS 1980(IUGG, 1980)
    #> 3                NHDArea       Polygon        1     14 GRS 1980(IUGG, 1980)
    #> 4           NHDWaterbody       Polygon       90     21 GRS 1980(IUGG, 1980)
    #> 5 NHDFlowline_NonNetwork   Line String       45     12 GRS 1980(IUGG, 1980)
    #> 6                   Gage         Point       33     19 GRS 1980(IUGG, 1980)

Now we have an output geopackage that can be used later. It contains the network subset of catchments and flowlines as well as a spatial subset of other layers as shown in the status output above. To complete the demonstration, here are a couple more layers plotted up.

    catchment <- sf::read_sf(output_file, "CatchmentSP")
    waterbody <- sf::read_sf(output_file, "NHDWaterbody")

    plot(sf::st_geometry(flowline))
    plot(start_point, cex = 1.5, lwd = 2, col = "red", add = TRUE)
    plot(sf::st_geometry(dplyr::filter(flowline, comid %in% UT_comids)),
         add=TRUE, col = "red", lwd = 2)
    plot(sf::st_geometry(catchment), add = TRUE)
    plot(sf::st_geometry(waterbody), col = rgb(0, 0, 1, alpha = 0.5), add = TRUE)

## Indexing

nhdplusTools supports a number of indexing use cases. See the function index for specifics.

Using the data above, we can use the `get_flowline_index()` function to get the comid, reachcode, and measure of our starting point like this.

    get_flowline_index(flowline, start_point)
    #>   id    COMID      REACHCODE REACH_meas       offset
    #> 1  1 13293750 07090002007373       41.8 0.0009619865

`get_flowline_index()` will work with a list of points too. For demonstration purposes, we can use the gages in our subset from above.

    gage <- sf::read_sf(output_file, "Gage")

    get_flowline_index(flowline, sf::st_geometry(gage), precision = 10)
    #>    id    COMID      REACHCODE REACH_meas         offset
    #> 1   2 13294276 07090002008387    14.8344 0.000012415352
    #> 2   3 13294264 07090002007650    56.4452 0.000020761363
    #> 3   4 13293750 07090002007373    42.5086 0.000025174363
    #> 4   5 13294312 07090002008383     1.2149 0.000002608798
    #> 5   6 13294264 07090002007650    41.0537 0.000032819820
    #> 6   7 13294264 07090002007650     2.0882 0.000002499335
    #> 7   9 13294300 07090002008379    85.4447 0.000016089689
    #> 8  10 13293690 07090002007648     0.7648 0.000014026465
    #> 9  11 13294264 07090002007650    71.4057 0.008205090544
    #> 10 13 13294176 07090002007664     6.1175 0.000032343949
    #> 11 15 13294290 07090002008374    88.0159 0.000028377490
    #> 12 16 13294138 07090002007709    73.3730 0.000011168868
    #> 13 18 13293486 07090002007724     8.8510 0.000012021333
    #> 14 19 13293512 07090002007723     5.9875 0.000010197748
    #> 15 20 13294176 07090002007664    16.9945 0.000012367772
    #> 16 22 13293474 07090002007713    72.0840 0.000005833049
    #> 17 23 13293520 07090002007676    18.1948 0.000030297908
    #> 18 26 13293512 07090002007723    68.1727 0.000026036784
    #> 19 27 13294264 07090002007650    97.1682 0.000014164679
    #> 20 29 13293456 07090002007738    90.8329 0.000036145292
    #> 21 31 13293512 07090002007723    67.6234 0.000022194820
    #> 22 33 13293576 07090002008236     2.7112 0.000009617344

For more info about `get_flowline_index()` and other indexing functions, see the article `vignette("indexing")` about it or the reference page that describes it.