Given a list of outlets, get their basin boundaries and network and return a plot in EPSG:3857 Web Mercator Projection.

plot_nhdplus(
  outlets = NULL,
  bbox = NULL,
  streamorder = NULL,
  nhdplus_data = NULL,
  gpkg = NULL,
  plot_config = NULL,
  basemap = "Esri.NatGeoWorldMap",
  zoom = NULL,
  add = FALSE,
  actually_plot = TRUE,
  overwrite = TRUE,
  flowline_only = NULL,
  cache_data = NULL
)

Arguments

outlets

list of nldi outlets. Other inputs are coerced into nldi outlets, see details.

bbox

object of class bbox with a defined crs. See examples.

streamorder

integer only streams of order greater than or equal will be returned

nhdplus_data

geopackage containing source nhdplus data (omit to download)

gpkg

path and file with .gpkg ending. If omitted, no file is written.

plot_config

list containing plot configuration, see details.

basemap

character indicating which basemap type to use. Chose from: get_tiles.

zoom

integer passed on to get_tiles. This value will override the default set by the package.

add

boolean should this plot be added to an already built map.

actually_plot

boolean actually draw the plot? Use to get data subset only.

overwrite

passed on the subset_nhdplus.

flowline_only

boolean only subset and plot flowlines only, default=FALSE

cache_data

character path to rds file where all plot data can be cached. If file doesn't exist, it will be created. If set to FALSE, all caching will be turned off – this includes basemap tiles.

Value

data.frame plot data is returned invisibly in NAD83 Lat/Lon.

Details

plot_nhdplus supports several input specifications. An unexported function "as_outlet" is used to convert the outlet formats as described below.

  1. if outlets is omitted, the bbox input is required and all nhdplus data in the bounding box is plotted.

  2. If outlets is a list of integers, it is assumed to be NHDPlus IDs (comids) and all upstream tributaries are plotted.

  3. if outlets is an integer vector, it is assumed to be all NHDPlus IDs (comids) that should be plotted. Allows custom filtering.

  4. If outlets is a character vector, it is assumed to be NWIS site ids.

  5. if outlets is a list containing only characters, it is assumed to be a list of nldi features and all upstream tributaries are plotted.

  6. if outlets is a data.frame with point geometry, a point in polygon match is performed and upstream with tributaries from the identified catchments is plotted.

The plot_config parameter is a list with names "basin", "flowline", "outlets", "network_wtbd", and "off_network_wtbd". The following shows the defaults that can be altered.

  1. basin

    
        list(lwd = 1, col = NA, border = "black")
  2. flowline

    
        list(lwd = 1, col = "blue")
  3. outlets

    
       list(default = list(col = "black", border = NA, pch = 19, cex = 1),
            nwissite = list(col = "grey40", border = NA, pch = 17, cex = 1),
            huc12pp = list(col = "white", border = "black", pch = 22, cex = 1),
            wqp = list(col = "red", border = NA, pch = 20, cex = 1))
  4. network_wtbd list(lwd = 1, col = "lightblue", border = "black")

  5. off_network_wtbd list(lwd = 1, col = "darkblue", border = "black")

If adding additional layers to the plot, data must be projected to EPSG:3857 with `sf::st_transform(x, 3857)` prior to adding to the plot.

Examples

# \donttest{
options("rgdal_show_exportToProj4_warnings"="none")
# Beware plot_nhdplus caches data to the default location.
# If you do not want data in "user space" change the default.
old_dir <- nhdplusTools::nhdplusTools_data_dir()
nhdplusTools_data_dir(tempdir())

plot_nhdplus("05428500")
#> 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


plot_nhdplus("05428500", streamorder = 2)
#> 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>1</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


plot_nhdplus(list(13293970, 13293750))
#> https://labs-beta.waterdata.usgs.gov/api/nldi/linked-data/comid/13293970/basin?simplify=true&splitCatchment=false
#> https://labs-beta.waterdata.usgs.gov/api/nldi/linked-data/comid/13293750/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>3</ogc:Literal>
#>         </ogc:PropertyIsGreaterThan>
#>         <ogc:BBOX>
#>           <ogc:PropertyName>the_geom</ogc:PropertyName>
#>           <gml:Envelope srsName="urn:x-ogc:def:crs:EPSG:4326">
#>             <gml:lowerCorner>42.950403765 -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: 10


source(system.file("extdata/sample_data.R", package = "nhdplusTools"))

plot_nhdplus(list(13293970, 13293750), streamorder = 3, nhdplus_data = sample_data)
#> Zoom set to: 10


plot_nhdplus(list(list("comid", "13293970"),
                  list("nwissite", "USGS-05428500"),
                  list("huc12pp", "070900020603"),
                  list("huc12pp", "070900020602")),
             streamorder = 2,
             nhdplus_data = sample_data)
#> Zoom set to: 10



plot_nhdplus(sf::st_as_sf(data.frame(x = -89.36083,
                                     y = 43.08944),
                          coords = c("x", "y"), crs = 4326),
             streamorder = 2,
             nhdplus_data = sample_data)
#> Zoom set to: 11


plot_nhdplus(list(list("comid", "13293970"),
                  list("nwissite", "USGS-05428500"),
                  list("huc12pp", "070900020603"),
                  list("huc12pp", "070900020602")),
             streamorder = 2,
             nhdplus_data = sample_data,
             plot_config = list(basin = list(lwd = 2),
                                outlets = list(huc12pp = list(cex = 1.5),
                                               comid = list(col = "green"))))
#> Zoom set to: 10


bbox <- sf::st_bbox(c(xmin = -89.43, ymin = 43, xmax = -89.28, ymax = 43.1),
                    crs = "+proj=longlat +datum=WGS84 +no_defs")

fline <- sf::read_sf(sample_data, "NHDFlowline_Network")
comids <- nhdplusTools::get_UT(fline, 13293970)

plot_nhdplus(comids)
#> 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>13293970</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>13293504</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>13294128</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>13293454</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>13293424</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>13293398</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>13293388</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>13293380</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>13294288</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>13293506</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>13294318</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294320</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294328</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294336</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294340</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294342</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294358</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294360</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294390</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13302592</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>13294298</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>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>13294346</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294002</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13297098</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13296564</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293996</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293998</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293990</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13296566</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293876</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293858</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294344</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294224</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294376</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293826</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294374</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>13293694</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>13293676</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>13293612</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>13294166</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>13293540</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>13293520</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>13294132</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>13293550</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>13293508</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>13293526</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>13294138</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>13293488</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>13293474</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>13293440</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>13293458</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>13294274</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>13293416</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>13293382</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>13293376</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>13293394</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>13293404</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>13294366</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>13293432</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>13293456</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>13294158</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>13293634</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>13293744</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293728</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293706</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293848</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293844</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293846</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293910</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294242</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293940</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293938</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294248</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294250</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293986</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13297324</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13296562</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293898</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293892</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293934</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>937090090</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>937090091</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>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>13293434</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>13294154</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>13294316</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294322</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294324</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294330</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294332</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294334</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294348</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294354</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294356</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294364</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13297162</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13297164</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>13294300</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>13294314</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>13294276</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>13294278</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>13293494</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>937090064</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13297104</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13297106</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293992</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13296582</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293994</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294228</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294396</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294378</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>13294184</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>13293586</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>13293624</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>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>13294176</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>13293578</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>13293548</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>13293476</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>13293442</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>13293472</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>13293568</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>13293552</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>13293522</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>13293532</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>13294148</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>13293516</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>13293498</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>13294122</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>13294112</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>13293486</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>13293730</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294194</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293592</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293594</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293852</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294238</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293946</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293972</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294254</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293958</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293942</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293896</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13293930</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>13293428</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>13293412</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>13293490</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>13294152</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294326</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294338</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13294350</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>         <ogc:PropertyIsEqualTo>
#>           <ogc:PropertyName>comid</ogc:PropertyName>
#>           <ogc:Literal>13297166</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>13294270</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>13293688</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>13294174</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>13293562</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>13293448</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>13294120</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>13293602</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>13293418</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>13293590</ogc:Literal>
#>         </ogc:PropertyIsEqualTo>
#>       </ogc:Or>
#>     </ogc:Filter>
#>   </wfs:Query>
#> </wfs:GetFeature>
#> Spherical geometry (s2) switched on
#> Zoom set to: 11


#' # With Local Data
plot_nhdplus(bbox = bbox, nhdplus_data = sample_data)
#> Zoom set to: 11


# With downloaded data
plot_nhdplus(bbox = bbox, streamorder = 3)
#> 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 -89.43</gml:lowerCorner>
#>             <gml:upperCorner>43.1 -89.28</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: 12


# Can also plot on top of the previous!
plot_nhdplus(bbox = bbox, nhdplus_data = sample_data,
             plot_config = list(flowline = list(lwd = 0.5)))
#> Zoom set to: 11
plot_nhdplus(comids, nhdplus_data = sample_data, streamorder = 3, add = TRUE,
             plot_config = list(flowline = list(col = "darkblue")))


nhdplusTools::nhdplusTools_data_dir(old_dir)
# }