diff --git a/docs/reference/add_plus_network_attributes.html b/docs/reference/add_plus_network_attributes.html index 976bdd78..61eff878 100644 --- a/docs/reference/add_plus_network_attributes.html +++ b/docs/reference/add_plus_network_attributes.html @@ -1,266 +1,261 @@ - - - -
- - - - -R/run_plus_attributes.R
- add_plus_network_attributes.Rd
Given a river network with required base attributes, adds the -NHDPlus network attributes: hydrosequence, levelpath, terminalpath, pathlength, -down levelpath, down hydroseq, total drainage area, and terminalflag. -The function implements two parallelization schemes for small and large basins -respectively. If a number of cores is specified, parallel execution will be used.
-add_plus_network_attributes(
- net,
- override = 5,
- cores = NULL,
- split_temp = NULL,
- status = TRUE
-)
net | -data.frame containing comid, tocomid, nameID, lengthkm, and areasqkm. -Additional attributes will be passed through unchanged. -tocomid == 0 is the convention used for outlets. -If a "weight" column is provided, it will be used in get_levelpaths -otherwise, arbolate sum is calculated for the network and used as the weight. |
-
---|---|
override | -numeric factor to be passed to get_levelpaths |
-
cores | -integer number of processes to spawn if run in parallel. |
-
split_temp | -character path to optional temporary copy of the network -split into independent sub-networks. If it exists, it will be read from disk -rather than recreated. |
-
status | -logical should progress be printed? |
-
-source(system.file("extdata", "walker_data.R", package = "nhdplusTools"))
-
-test_flowline <- prepare_nhdplus(walker_flowline, 0, 0, FALSE)
-#> Warning: removing geometry
-#> Warning: Removed 0 flowlines that don't apply.
-#> Includes: Coastlines, non-dendritic paths,
-#> and networks with drainage area less than 0 sqkm, and drainage basins smaller than FALSE
-
-test_flowline <- data.frame(
- comid = test_flowline$COMID,
- tocomid = test_flowline$toCOMID,
- nameID = walker_flowline$GNIS_ID,
- lengthkm = test_flowline$LENGTHKM,
- areasqkm = walker_flowline$AreaSqKM)
-
-add_plus_network_attributes(test_flowline)
-#> # A tibble: 62 x 14
-#> comid tocomid nameID lengthkm areasqkm weight terminalpa hydroseq levelpathi
-#> <dbl> <dbl> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
-#> 1 5.33e6 0 255208 1.20 0.868 137. 1 1 1
-#> 2 5.33e6 5329303 255208 0.582 0.260 133. 1 2 1
-#> 3 5.33e6 5329293 255208 4.02 3.76 127. 1 3 1
-#> 4 5.33e6 5329305 255208 4.44 4.51 120. 1 4 1
-#> 5 5.33e6 5329317 255208 2.68 3.64 112. 1 5 1
-#> 6 5.33e6 5329315 255208 4.79 5.96 75.3 1 6 1
-#> 7 5.33e6 5329339 255208 1.03 1.27 67.4 1 7 1
-#> 8 5.33e6 5329343 255208 3.60 6.94 61.9 1 8 1
-#> 9 5.33e6 5329357 255208 3.14 5.13 51.1 1 9 1
-#> 10 5.33e6 5329365 254577 1.46 1.69 37.8 1 10 1
-#> # ... with 52 more rows, and 5 more variables: pathlength <dbl>,
-#> # dnlevelpat <dbl>, dnhydroseq <dbl>, totdasqkm <dbl>, terminalfl <dbl>
-
R/run_plus_attributes.R
+ add_plus_network_attributes.Rd
Given a river network with required base attributes, adds the +NHDPlus network attributes: hydrosequence, levelpath, terminalpath, pathlength, +down levelpath, down hydroseq, total drainage area, and terminalflag. +The function implements two parallelization schemes for small and large basins +respectively. If a number of cores is specified, parallel execution will be used.
+add_plus_network_attributes(
+ net,
+ override = 5,
+ cores = NULL,
+ split_temp = NULL,
+ status = TRUE
+)
net | +data.frame containing comid, tocomid, nameID, lengthkm, and areasqkm. +Additional attributes will be passed through unchanged. +tocomid == 0 is the convention used for outlets. +If a "weight" column is provided, it will be used in get_levelpaths +otherwise, arbolate sum is calculated for the network and used as the weight. |
+
---|---|
override | +numeric factor to be passed to get_levelpaths |
+
cores | +integer number of processes to spawn if run in parallel. |
+
split_temp | +character path to optional temporary copy of the network +split into independent sub-networks. If it exists, it will be read from disk +rather than recreated. |
+
status | +logical should progress be printed? |
+
+source(system.file("extdata", "walker_data.R", package = "nhdplusTools"))
+
+test_flowline <- prepare_nhdplus(walker_flowline, 0, 0, FALSE)
+#> Warning: removing geometry
+#> Warning: Removed 0 flowlines that don't apply.
+#> Includes: Coastlines, non-dendritic paths,
+#> and networks with drainage area less than 0 sqkm, and drainage basins smaller than FALSE
+
+test_flowline <- data.frame(
+ comid = test_flowline$COMID,
+ tocomid = test_flowline$toCOMID,
+ nameID = walker_flowline$GNIS_ID,
+ lengthkm = test_flowline$LENGTHKM,
+ areasqkm = walker_flowline$AreaSqKM)
+
+add_plus_network_attributes(test_flowline)
+#>
| | 0 % ~calculating
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s
+#>
| | 0 % ~calculating
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s
+#> # A tibble: 62 x 14
+#> comid tocomid nameID lengthkm areasqkm weight terminalpa hydroseq levelpathi
+#> <dbl> <dbl> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
+#> 1 5.33e6 0 255208 1.20 0.868 137. 1 1 1
+#> 2 5.33e6 5329303 255208 0.582 0.260 133. 1 2 1
+#> 3 5.33e6 5329293 255208 4.02 3.76 127. 1 3 1
+#> 4 5.33e6 5329305 255208 4.44 4.51 120. 1 4 1
+#> 5 5.33e6 5329317 255208 2.68 3.64 112. 1 5 1
+#> 6 5.33e6 5329315 255208 4.79 5.96 75.3 1 6 1
+#> 7 5.33e6 5329339 255208 1.03 1.27 67.4 1 7 1
+#> 8 5.33e6 5329343 255208 3.60 6.94 61.9 1 8 1
+#> 9 5.33e6 5329357 255208 3.14 5.13 51.1 1 9 1
+#> 10 5.33e6 5329365 254577 1.46 1.69 37.8 1 10 1
+#> # ... with 52 more rows, and 5 more variables: pathlength <dbl>,
+#> # dnlevelpat <dbl>, dnhydroseq <dbl>, totdasqkm <dbl>, terminalfl <dbl>
+
library(dplyr)
#>
-#> Attaching package: 'dplyr'
-#> The following objects are masked from 'package:stats':
+#> Attaching package: ‘dplyr’
+#> The following object is masked from ‘package:testthat’:
+#>
+#> matches
+#> The following objects are masked from ‘package:stats’:
#>
#> filter, lag
-#> The following objects are masked from 'package:base':
+#> The following objects are masked from ‘package:base’:
#>
#> intersect, setdiff, setequal, union
source(system.file("extdata", "walker_data.R", package = "nhdplusTools"))
diff --git a/docs/reference/calculate_total_drainage_area.html b/docs/reference/calculate_total_drainage_area.html
index 20e65c4a..c4dfd013 100644
--- a/docs/reference/calculate_total_drainage_area.html
+++ b/docs/reference/calculate_total_drainage_area.html
@@ -8,13 +8,6 @@
Total Drainage Area — calculate_total_drainage_area • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/disambiguate_flowline_indexes.html b/docs/reference/disambiguate_flowline_indexes.html
index d9d58af5..6f2f2b99 100644
--- a/docs/reference/disambiguate_flowline_indexes.html
+++ b/docs/reference/disambiguate_flowline_indexes.html
@@ -8,13 +8,6 @@
Disambiguate Flowline Indexes — disambiguate_flowline_indexes • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/discover_nhdplus_id.html b/docs/reference/discover_nhdplus_id.html
index 07cd1a8a..15fb4b13 100644
--- a/docs/reference/discover_nhdplus_id.html
+++ b/docs/reference/discover_nhdplus_id.html
@@ -8,13 +8,6 @@
Discover NHDPlus ID — discover_nhdplus_id • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/discover_nldi_characteristics.html b/docs/reference/discover_nldi_characteristics.html
index 1f5673be..bdb78fa8 100644
--- a/docs/reference/discover_nldi_characteristics.html
+++ b/docs/reference/discover_nldi_characteristics.html
@@ -8,13 +8,6 @@
Discover Characteristics Metadata — discover_nldi_characteristics • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/download_nhdplushr.html b/docs/reference/download_nhdplushr.html
index cdc86d05..a1dcf8fa 100644
--- a/docs/reference/download_nhdplushr.html
+++ b/docs/reference/download_nhdplushr.html
@@ -8,13 +8,6 @@
Download NHDPlus HiRes — download_nhdplushr • nhdplusTools
-
-
-
-
-
-
-
@@ -181,6 +174,7 @@ Value
Examples
# \donttest{
hu <- nhdplusTools::get_huc8(sf::st_sfc(sf::st_point(c(-73, 42)), crs = 4326))
+#> Found invalid geometry, attempting to fix.
#> although coordinates are longitude/latitude, st_intersects assumes that they are planar
(hu <- substr(hu$huc8, 1, 2))
diff --git a/docs/reference/download_nhdplusv2.html b/docs/reference/download_nhdplusv2.html
index f711b812..837c6555 100644
--- a/docs/reference/download_nhdplusv2.html
+++ b/docs/reference/download_nhdplusv2.html
@@ -8,13 +8,6 @@
Download seamless National Hydrography Dataset Version 2 (NHDPlusV2) — download_nhdplusv2 • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/download_rf1.html b/docs/reference/download_rf1.html
index 18997a07..8e891a3a 100644
--- a/docs/reference/download_rf1.html
+++ b/docs/reference/download_rf1.html
@@ -8,13 +8,6 @@
Download the seamless Reach File (RF1) Database — download_rf1 • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/download_vaa.html b/docs/reference/download_vaa.html
index 326ee2af..b05f1f67 100644
--- a/docs/reference/download_vaa.html
+++ b/docs/reference/download_vaa.html
@@ -8,13 +8,6 @@
Download nhdplusVAA data from HydroShare — download_vaa • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/download_wbd.html b/docs/reference/download_wbd.html
index 7b48843e..23c169d3 100644
--- a/docs/reference/download_wbd.html
+++ b/docs/reference/download_wbd.html
@@ -8,13 +8,6 @@
Download the seamless Watershed Boundary Dataset (WBD) — download_wbd • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_DD.html b/docs/reference/get_DD.html
index a526ec5c..b9aaca61 100644
--- a/docs/reference/get_DD.html
+++ b/docs/reference/get_DD.html
@@ -8,13 +8,6 @@
Navigate Downstream with Diversions — get_DD • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_DM.html b/docs/reference/get_DM.html
index 9d0d9cc7..7cf6b4cf 100644
--- a/docs/reference/get_DM.html
+++ b/docs/reference/get_DM.html
@@ -8,13 +8,6 @@
Navigate Downstream Mainstem — get_DM • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_UM.html b/docs/reference/get_UM.html
index 8d0808e5..c954c76f 100644
--- a/docs/reference/get_UM.html
+++ b/docs/reference/get_UM.html
@@ -8,13 +8,6 @@
Navigate Upstream Mainstem — get_UM • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_UT.html b/docs/reference/get_UT.html
index 8b747d0a..d52e6c28 100644
--- a/docs/reference/get_UT.html
+++ b/docs/reference/get_UT.html
@@ -8,13 +8,6 @@
Navigate Upstream with Tributaries — get_UT • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_flowline_index.html b/docs/reference/get_flowline_index.html
index c9a0d6c1..c731cdf1 100644
--- a/docs/reference/get_flowline_index.html
+++ b/docs/reference/get_flowline_index.html
@@ -8,13 +8,6 @@
Get Flowline Index — get_flowline_index • nhdplusTools
-
-
-
-
-
-
-
@@ -255,16 +248,16 @@ Examp
#> # A tibble: 30 x 5
#> id COMID REACHCODE REACH_meas offset
#> <int> <int> <chr> <dbl> <dbl>
-#> 1 1 11688298 02060003000579 0 0.0000603
-#> 2 2 11688808 02060003000519 53.6 0.000564
-#> 3 3 11688980 02060003000253 75.4 0.000310
-#> 4 1 11689926 02060003001467 100 0.0000603
-#> 5 3 11688950 02060003000254 18.5 0.00113
-#> 6 1 11689928 02060003001468 0 0.00203
-#> 7 1 11689978 02060003001472 100 0.00203
-#> 8 1 11690530 02060003000585 0 0.00285
-#> 9 1 11690490 02060003000580 0 0.00285
-#> 10 3 11688948 02060003000516 0 0.00321
+#> 1 1 11688298 02060003000579 1.42e-14 0.0000603
+#> 2 2 11688808 02060003000519 5.36e+ 1 0.000564
+#> 3 3 11688980 02060003000253 7.54e+ 1 0.000310
+#> 4 1 11689926 02060003001467 1 e+ 2 0.0000603
+#> 5 3 11688950 02060003000254 1.85e+ 1 0.00113
+#> 6 1 11689928 02060003001468 1.42e-14 0.00203
+#> 7 1 11689978 02060003001472 1 e+ 2 0.00203
+#> 8 1 11690530 02060003000585 2.64e-15 0.00285
+#> 9 1 11690490 02060003000580 0 0.00285
+#> 10 3 11688948 02060003000516 0 0.00321
#> # ... with 20 more rows
# }
diff --git a/docs/reference/get_gagesII.html b/docs/reference/get_gagesII.html
index 30ab9cee..62035638 100644
--- a/docs/reference/get_gagesII.html
+++ b/docs/reference/get_gagesII.html
@@ -8,13 +8,6 @@
Find GAGESII Features — get_gagesII • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_hr_data.html b/docs/reference/get_hr_data.html
index 7d52ad62..e84718f1 100644
--- a/docs/reference/get_hr_data.html
+++ b/docs/reference/get_hr_data.html
@@ -8,13 +8,6 @@
Get NHDPlus HiRes Data — get_hr_data • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_huc12.html b/docs/reference/get_huc12.html
index b1b1e9e3..943377b9 100644
--- a/docs/reference/get_huc12.html
+++ b/docs/reference/get_huc12.html
@@ -8,13 +8,6 @@
Find WBD HUC 12 unit subsets — get_huc12 • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_huc8.html b/docs/reference/get_huc8.html
index 0a5a7dd4..d6e78f05 100644
--- a/docs/reference/get_huc8.html
+++ b/docs/reference/get_huc8.html
@@ -8,13 +8,6 @@
Find WBD HUC 08 unit subsets — get_huc8 • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_hydro_location.html b/docs/reference/get_hydro_location.html
index 3f2d4132..051a95d2 100644
--- a/docs/reference/get_hydro_location.html
+++ b/docs/reference/get_hydro_location.html
@@ -8,13 +8,6 @@
get hydro location — get_hydro_location • nhdplusTools
-
-
-
-
-
-
-
@@ -161,14 +154,12 @@ Arg
indexes
- data.frame as output from get_flowline_index with more than
-one hydrologic location per indexed point.
+ data.frame as output from get_flowline_index.
flowpath
data.frame with two columns. The first should join to the COMID
-field of the indexes and the second should be the numeric or ascii metric such as drainage
-area or GNIS Name. Names of this data.frame are not used.
+field of the indexes and the second should be linear geometry.
diff --git a/docs/reference/get_levelpaths.html b/docs/reference/get_levelpaths.html
index 3b1f9173..758d46a8 100644
--- a/docs/reference/get_levelpaths.html
+++ b/docs/reference/get_levelpaths.html
@@ -8,13 +8,6 @@
Get Level Paths — get_levelpaths • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_nhdarea.html b/docs/reference/get_nhdarea.html
index 49d5864c..8204c22c 100644
--- a/docs/reference/get_nhdarea.html
+++ b/docs/reference/get_nhdarea.html
@@ -8,13 +8,6 @@
Find NHD Areas — get_nhdarea • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_nhdplus.html b/docs/reference/get_nhdplus.html
index 5b963ffc..51f7a766 100644
--- a/docs/reference/get_nhdplus.html
+++ b/docs/reference/get_nhdplus.html
@@ -8,13 +8,6 @@
Get National Hydrography Dataset V2 Subsets (Multirealization) — get_nhdplus • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/get_nhdplushr.html b/docs/reference/get_nhdplushr.html
index f97f6394..a8ebf087 100644
--- a/docs/reference/get_nhdplushr.html
+++ b/docs/reference/get_nhdplushr.html
@@ -8,13 +8,6 @@
Get NHDPlus HiRes — get_nhdplushr • nhdplusTools
-
-
-
-
-
-
-
@@ -224,7 +217,7 @@ Details
Set layers to "NULL" to get all layers.
Examples
- # \donttest{
+ if (FALSE) {
# Note this will download a lot of data to a temp directory.
# Change 'temp_dir' to your directory of choice.
temp_dir <- file.path(nhdplusTools_data_dir(), "temp_hr_cache")
@@ -232,537 +225,20 @@ Examp
download_dir <- download_nhdplushr(temp_dir, c("0302", "0303"))
get_nhdplushr(download_dir, file.path(download_dir, "nhdplus_0302-03.gpkg"))
-#> $NHDFlowline
-#> Simple feature collection with 211528 features and 57 fields (with 2150 geometries empty)
-#> Geometry type: MULTILINESTRING
-#> Dimension: XY
-#> Bounding box: xmin: -80.06746 ymin: 33.85333 xmax: -75.48454 ymax: 36.44617
-#> Geodetic CRS: NAD83
-#> # A tibble: 211,528 x 58
-#> Permanent_Identif~ FDate Resolution GNIS_ID GNIS_Name LENGTHKM
-#> * <chr> <dttm> <int> <chr> <chr> <dbl>
-#> 1 31395075 2012-03-10 19:59:14 2 NA NA 0.131
-#> 2 31395981 2012-03-10 21:05:34 2 NA NA 0.031
-#> 3 31395993 2012-03-10 21:05:34 2 NA NA 0.062
-#> 4 31393833 2012-03-10 21:01:18 2 NA NA 0.223
-#> 5 31394775 2012-03-10 21:13:48 2 NA NA 0.017
-#> 6 31374279 2012-03-10 19:59:53 2 009821~ Buffalo C~ 0.029
-#> 7 31374951 2012-03-10 21:00:49 2 NA NA 0.074
-#> 8 31374695 2012-03-10 21:06:58 2 NA NA 0.684
-#> 9 31374725 2012-03-10 21:06:58 2 NA NA 0.099
-#> 10 31373805 2012-03-10 19:59:53 2 009821~ Buffalo C~ 0.092
-#> # ... with 211,518 more rows, and 52 more variables: REACHCODE <chr>,
-#> # FlowDir <int>, WBArea_Permanent_Identifier <chr>, FTYPE <int>, FCODE <int>,
-#> # MainPath <int>, InNetwork <int>, VisibilityFilter <int>,
-#> # Shape_Length <dbl>, COMID <dbl>, VPUID <chr>, Enabled <int>,
-#> # Shape <MULTILINESTRING [°]>, StreamLeve <int>, StreamOrde <int>,
-#> # StreamCalc <int>, FromNode <dbl>, ToNode <dbl>, Hydroseq <dbl>,
-#> # LevelPathI <dbl>, Pathlength <dbl>, TerminalPa <dbl>, ArbolateSu <dbl>,
-#> # Divergence <int>, StartFlag <int>, TerminalFl <dbl>, UpLevelPat <dbl>,
-#> # UpHydroSeq <dbl>, DnLevel <dbl>, DnLevelPat <dbl>, DnHydroseq <dbl>,
-#> # DnMinorHyd <dbl>, DnDrainCou <int>, FromMeas <dbl>, ToMeas <dbl>,
-#> # RtnDiv <int>, Thinner <int>, VPUIn <int>, VPUOut <int>, AreaSqKM <dbl>,
-#> # TotDASqKM <dbl>, DivDASqKm <dbl>, MaxElevRaw <dbl>, MinElevRaw <dbl>,
-#> # MaxElevSmo <dbl>, MinElevSmo <dbl>, Slope <dbl>, SlopeLenKm <dbl>,
-#> # ElevFixed <int>, HWType <int>, HWNodeSqKm <dbl>, StatusFlag <chr>
-#>
-#> $NHDPlusCatchment
-#> Simple feature collection with 185946 features and 7 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -80.07324 ymin: 33.84224 xmax: -75.48088 ymax: 36.45034
-#> Geodetic CRS: NAD83
-#> # A tibble: 185,946 x 8
-#> FEATUREID SourceFC GridCode AreaSqKM VPUID SHAPE_Length SHAPE_Area
-#> * <dbl> <chr> <int> <dbl> <chr> <dbl> <dbl>
-#> 1 1.50e13 NHDPlusBurnLineEve~ 9810 0.574 0302 0.0458 0.0000576
-#> 2 1.50e13 NHDPlusBurnLineEve~ 11756 0.194 0302 0.0263 0.0000194
-#> 3 1.50e13 NHDPlusBurnLineEve~ 69563 0.201 0302 0.0253 0.0000202
-#> 4 1.50e13 NHDPlusBurnLineEve~ 100062 0.250 0302 0.0340 0.0000251
-#> 5 1.50e13 NHDPlusBurnLineEve~ 90919 0.176 0302 0.0215 0.0000177
-#> 6 1.50e13 NHDPlusBurnLineEve~ 13162 0.0281 0302 0.00827 0.00000282
-#> 7 1.50e13 NHDPlusBurnLineEve~ 83340 0.258 0302 0.0314 0.0000259
-#> 8 1.50e13 NHDPlusBurnLineEve~ 12802 0.386 0302 0.0360 0.0000388
-#> 9 1.50e13 NHDPlusBurnLineEve~ 9775 0.791 0302 0.0526 0.0000794
-#> 10 1.50e13 NHDPlusBurnLineEve~ 12043 0.357 0302 0.0415 0.0000359
-#> # ... with 185,936 more rows, and 1 more variable: SHAPE <MULTIPOLYGON [°]>
-#>
get_nhdplushr(download_dir,
file.path(download_dir, "nhdplus_0302-03.gpkg"),
layers = NULL, overwrite = TRUE)
-#> Warning: GDAL Message 1: organizePolygons() received a polygon with more than 100 parts. The processing may be really slow. You can skip the processing by setting METHOD=SKIP, or only make it analyze counter-clock wise parts by setting METHOD=ONLY_CCW if you can assume that the outline of holes is counter-clock wise defined
-#> $NHDPlusBurnLineEvent
-#> Simple feature collection with 186937 features and 12 fields
-#> Geometry type: MULTILINESTRING
-#> Dimension: XY
-#> Bounding box: xmin: -80.06698 ymin: 33.84234 xmax: -75.48283 ymax: 36.44529
-#> Geodetic CRS: NAD83
-#> # A tibble: 186,937 x 13
-#> FEATUREID GridCode REACHCODE FromMeas ToMeas BurnLenKm RPUID StatusFlag
-#> * <dbl> <int> <chr> <dbl> <dbl> <dbl> <chr> <chr>
-#> 1 1.50e13 960 03020101000625 79.4 84.3 0.334 0302 A
-#> 2 1.50e13 961 03020101000933 0 2.51 0.063 0302 A
-#> 3 1.50e13 962 03020101000170 62.7 78.3 0.385 0302 A
-#> 4 1.50e13 963 03020101000032 66.0 80.9 0.47 0302 A
-#> 5 1.50e13 1188 03020101000815 0 18.1 0.676 0302 A
-#> 6 1.50e13 1189 03020101001059 36.7 54.2 0.183 0302 A
-#> 7 1.50e13 1190 03020101000366 76.8 100 0.464 0302 A
-#> 8 1.50e13 1191 03020101000717 51.4 54.6 0.051 0302 A
-#> 9 1.50e13 1192 03020101000601 4.96 11.9 0.178 0302 A
-#> 10 1.50e13 1193 03020101000535 0 19.6 0.191 0302 A
-#> # ... with 186,927 more rows, and 5 more variables: Catchment <int>,
-#> # Burn <int>, VPUID <chr>, SHAPE_Length <dbl>, SHAPE <MULTILINESTRING [°]>
-#>
-#> $NHDPlusBurnWaterbody
-#> Simple feature collection with 51147 features and 8 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -80.06802 ymin: 33.8532 xmax: -75.48533 ymax: 36.4423
-#> Geodetic CRS: NAD83
-#> # A tibble: 51,147 x 9
-#> COMID SourceFC OnOffNet PurpCode Burn VPUID SHAPE_Length SHAPE_Area
-#> * <dbl> <chr> <int> <chr> <int> <chr> <dbl> <dbl>
-#> 1 1.50e13 NHDWaterbody 1 WB 1 0302 0.00320 0.000000613
-#> 2 1.50e13 NHDWaterbody 1 WB 1 0302 0.00151 0.000000151
-#> 3 1.50e13 NHDWaterbody 0 WB 1 0302 0.00213 0.000000261
-#> 4 1.50e13 NHDWaterbody 0 WB 1 0302 0.000734 0.0000000265
-#> 5 1.50e13 NHDWaterbody 1 WB 1 0302 0.00170 0.000000149
-#> 6 1.50e13 NHDWaterbody 1 WB 1 0302 0.00514 0.00000127
-#> 7 1.50e13 NHDWaterbody 1 WB 1 0302 0.00246 0.000000282
-#> 8 1.50e13 NHDWaterbody 0 WB 1 0302 0.000671 0.0000000318
-#> 9 1.50e13 NHDWaterbody 0 WB 1 0302 0.00180 0.000000193
-#> 10 1.50e13 NHDWaterbody 0 WB 1 0302 0.00358 0.000000948
-#> # ... with 51,137 more rows, and 1 more variable: SHAPE <MULTIPOLYGON [°]>
-#>
-#> $NHDPlusCatchment
-#> Simple feature collection with 185946 features and 7 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -80.07324 ymin: 33.84224 xmax: -75.48088 ymax: 36.45034
-#> Geodetic CRS: NAD83
-#> # A tibble: 185,946 x 8
-#> FEATUREID SourceFC GridCode AreaSqKM VPUID SHAPE_Length SHAPE_Area
-#> * <dbl> <chr> <int> <dbl> <chr> <dbl> <dbl>
-#> 1 1.50e13 NHDPlusBurnLineEve~ 9810 0.574 0302 0.0458 0.0000576
-#> 2 1.50e13 NHDPlusBurnLineEve~ 11756 0.194 0302 0.0263 0.0000194
-#> 3 1.50e13 NHDPlusBurnLineEve~ 69563 0.201 0302 0.0253 0.0000202
-#> 4 1.50e13 NHDPlusBurnLineEve~ 100062 0.250 0302 0.0340 0.0000251
-#> 5 1.50e13 NHDPlusBurnLineEve~ 90919 0.176 0302 0.0215 0.0000177
-#> 6 1.50e13 NHDPlusBurnLineEve~ 13162 0.0281 0302 0.00827 0.00000282
-#> 7 1.50e13 NHDPlusBurnLineEve~ 83340 0.258 0302 0.0314 0.0000259
-#> 8 1.50e13 NHDPlusBurnLineEve~ 12802 0.386 0302 0.0360 0.0000388
-#> 9 1.50e13 NHDPlusBurnLineEve~ 9775 0.791 0302 0.0526 0.0000794
-#> 10 1.50e13 NHDPlusBurnLineEve~ 12043 0.357 0302 0.0415 0.0000359
-#> # ... with 185,936 more rows, and 1 more variable: SHAPE <MULTIPOLYGON [°]>
-#>
-#> $NHDPlusLandSea
-#> Simple feature collection with 367 features and 5 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -79.28467 ymin: 33.72576 xmax: -75.42428 ymax: 36.49496
-#> Geodetic CRS: NAD83
-#> # A tibble: 367 x 6
-#> COMID Land VPUID SHAPE_Length SHAPE_Area SHAPE
-#> * <dbl> <int> <chr> <dbl> <dbl> <MULTIPOLYGON [°]>
-#> 1 1.50e13 1 0302 0.381 1.63e-3 (((-77.81398 34.18475, -77.8138~
-#> 2 1.50e13 1 0302 0.0621 1.50e-4 (((-77.82168 34.19019, -77.8216~
-#> 3 1.50e13 1 0302 0.00464 1.31e-6 (((-77.80739 34.19895, -77.8073~
-#> 4 1.50e13 1 0302 0.00726 2.67e-6 (((-77.80622 34.2013, -77.80614~
-#> 5 1.50e13 1 0302 0.00135 9.27e-8 (((-77.74916 34.30372, -77.7491~
-#> 6 1.50e13 1 0302 0.0219 1.25e-5 (((-77.72203 34.31143, -77.7219~
-#> 7 1.50e13 1 0302 0.229 1.52e-4 (((-77.72629 34.31415, -77.7261~
-#> 8 1.50e13 1 0302 0.0608 7.12e-5 (((-77.71728 34.32186, -77.7171~
-#> 9 1.50e13 1 0302 0.0360 9.85e-6 (((-77.34969 34.54426, -77.3493~
-#> 10 1.50e13 1 0302 0.00968 2.29e-6 (((-77.33952 34.54479, -77.3395~
-#> # ... with 357 more rows
-#>
-#> $NHDPlusSink
-#> Simple feature collection with 369 features and 10 fields
-#> Geometry type: POINT
-#> Dimension: XY
-#> Bounding box: xmin: -79.81687 ymin: 33.93375 xmax: -75.77991 ymax: 36.25015
-#> Geodetic CRS: NAD83
-#> # A tibble: 369 x 11
-#> COMID GridCode PurpCode FEATUREID SourceFC RPUID StatusFlag Catchment Burn
-#> * <dbl> <int> <chr> <dbl> <chr> <chr> <chr> <int> <int>
-#> 1 1.50e13 78587 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 2 1.50e13 73682 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 3 1.50e13 62479 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 4 1.50e13 50401 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 5 1.50e13 87344 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 6 1.50e13 43850 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 7 1.50e13 18095 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 8 1.50e13 90089 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 9 1.50e13 45221 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> 10 1.50e13 4438 SE 1.50e13 NHDFlow~ 0302 N 1 1
-#> # ... with 359 more rows, and 2 more variables: VPUID <chr>, SHAPE <POINT [°]>
-#>
-#> $NHDPlusWall
-#> Simple feature collection with 2320 features and 5 fields
-#> Geometry type: MULTILINESTRING
-#> Dimension: XY
-#> Bounding box: xmin: -80.07334 ymin: 33.76852 xmax: -75.47986 ymax: 36.45033
-#> Geodetic CRS: NAD83
-#> # A tibble: 2,320 x 6
-#> COMID WallSource VPUID SHAPE_Length Burn SHAPE
-#> * <dbl> <chr> <chr> <dbl> <int> <MULTILINESTRING [°]>
-#> 1 1.50e13 WBDLine 0302 0.0348 1 ((-77.24908 35.3333, -77.24842 3~
-#> 2 1.50e13 WBDLine 0302 0.0607 1 ((-76.61913 35.23761, -76.61766 ~
-#> 3 1.50e13 WBDLine 0302 0.0497 1 ((-77.85314 35.47356, -77.85306 ~
-#> 4 1.50e13 WBDLine 0302 0.128 1 ((-77.54543 35.50729, -77.54508 ~
-#> 5 1.50e13 WBDLine 0302 0.108 1 ((-77.61271 35.4832, -77.61223 3~
-#> 6 1.50e13 WBDLine 0302 0.0921 1 ((-78.84642 35.75455, -78.84617 ~
-#> 7 1.50e13 WBDLine 0302 0.0368 1 ((-76.83527 35.06029, -76.83515 ~
-#> 8 1.50e13 WBDLine 0302 0.0648 0 ((-76.51809 35.13499, -76.51724 ~
-#> 9 1.50e13 WBDLine 0302 0.0214 1 ((-77.00828 35.14243, -77.00822 ~
-#> 10 1.50e13 WBDLine 0302 0.111 1 ((-77.44302 35.11485, -77.44265 ~
-#> # ... with 2,310 more rows
-#>
-#> $NHDArea
-#> Simple feature collection with 1129 features and 14 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XYZ
-#> Bounding box: xmin: -79.9673 ymin: 33.75017 xmax: -75.37459 ymax: 37.00015
-#> z_range: zmin: 0 zmax: 0
-#> Geodetic CRS: NAD83
-#> # A tibble: 1,129 x 15
-#> Permanent_Identifi~ FDate Resolution GNIS_ID GNIS_Name AreaSqKM
-#> * <chr> <dttm> <int> <chr> <chr> <dbl>
-#> 1 139013696 2005-08-19 23:44:15 2 NA NA 0.042
-#> 2 139010028 2005-08-19 23:39:26 2 NA NA 0.011
-#> 3 139024288 2005-08-19 23:58:58 2 NA NA 0.009
-#> 4 139024269 2005-08-19 23:58:53 2 NA NA 0.001
-#> 5 136049893 2005-06-07 04:39:34 2 NA NA 0.013
-#> 6 132126143 2005-02-10 09:47:21 2 NA NA 0.041
-#> 7 141615988 2005-10-25 12:09:29 2 NA NA 0.008
-#> 8 157525776 2007-06-22 03:31:36 2 NA NA 0.003
-#> 9 157526605 2007-06-22 03:32:24 2 NA NA 0.367
-#> 10 129758127 2004-12-02 10:38:02 2 NA NA 2.96
-#> # ... with 1,119 more rows, and 9 more variables: Elevation <dbl>, FTYPE <int>,
-#> # FCODE <int>, VisibilityFilter <int>, Shape_Length <dbl>, Shape_Area <dbl>,
-#> # COMID <dbl>, VPUID <chr>, Shape <MULTIPOLYGON [°]>
-#>
-#> $NHDFlowline
-#> Simple feature collection with 211528 features and 57 fields (with 2150 geometries empty)
-#> Geometry type: MULTILINESTRING
-#> Dimension: XY
-#> Bounding box: xmin: -80.06746 ymin: 33.85333 xmax: -75.48454 ymax: 36.44617
-#> Geodetic CRS: NAD83
-#> # A tibble: 211,528 x 58
-#> Permanent_Identif~ FDate Resolution GNIS_ID GNIS_Name LENGTHKM
-#> * <chr> <dttm> <int> <chr> <chr> <dbl>
-#> 1 31395075 2012-03-10 19:59:14 2 NA NA 0.131
-#> 2 31395981 2012-03-10 21:05:34 2 NA NA 0.031
-#> 3 31395993 2012-03-10 21:05:34 2 NA NA 0.062
-#> 4 31393833 2012-03-10 21:01:18 2 NA NA 0.223
-#> 5 31394775 2012-03-10 21:13:48 2 NA NA 0.017
-#> 6 31374279 2012-03-10 19:59:53 2 009821~ Buffalo C~ 0.029
-#> 7 31374951 2012-03-10 21:00:49 2 NA NA 0.074
-#> 8 31374695 2012-03-10 21:06:58 2 NA NA 0.684
-#> 9 31374725 2012-03-10 21:06:58 2 NA NA 0.099
-#> 10 31373805 2012-03-10 19:59:53 2 009821~ Buffalo C~ 0.092
-#> # ... with 211,518 more rows, and 52 more variables: REACHCODE <chr>,
-#> # FlowDir <int>, WBArea_Permanent_Identifier <chr>, FTYPE <int>, FCODE <int>,
-#> # MainPath <int>, InNetwork <int>, VisibilityFilter <int>,
-#> # Shape_Length <dbl>, COMID <dbl>, VPUID <chr>, Enabled <int>,
-#> # Shape <MULTILINESTRING [°]>, StreamLeve <int>, StreamOrde <int>,
-#> # StreamCalc <int>, FromNode <dbl>, ToNode <dbl>, Hydroseq <dbl>,
-#> # LevelPathI <dbl>, Pathlength <dbl>, TerminalPa <dbl>, ArbolateSu <dbl>,
-#> # Divergence <int>, StartFlag <int>, TerminalFl <dbl>, UpLevelPat <dbl>,
-#> # UpHydroSeq <dbl>, DnLevel <dbl>, DnLevelPat <dbl>, DnHydroseq <dbl>,
-#> # DnMinorHyd <dbl>, DnDrainCou <int>, FromMeas <dbl>, ToMeas <dbl>,
-#> # RtnDiv <int>, Thinner <int>, VPUIn <int>, VPUOut <int>, AreaSqKM <dbl>,
-#> # TotDASqKM <dbl>, DivDASqKm <dbl>, MaxElevRaw <dbl>, MinElevRaw <dbl>,
-#> # MaxElevSmo <dbl>, MinElevSmo <dbl>, Slope <dbl>, SlopeLenKm <dbl>,
-#> # ElevFixed <int>, HWType <int>, HWNodeSqKm <dbl>, StatusFlag <chr>
-#>
-#> $NHDLine
-#> Simple feature collection with 3455 features and 12 fields
-#> Geometry type: MULTILINESTRING
-#> Dimension: XYZ
-#> Bounding box: xmin: -80.06293 ymin: 33.89045 xmax: -75.50004 ymax: 36.43778
-#> z_range: zmin: 0 zmax: 0
-#> Geodetic CRS: NAD83
-#> # A tibble: 3,455 x 13
-#> Permanent_Identifi~ FDate Resolution GNIS_ID GNIS_Name LENGTHKM
-#> * <chr> <dttm> <int> <chr> <chr> <dbl>
-#> 1 139025426 2005-08-20 00:00:37 2 NA NA 0.078
-#> 2 139012370 2005-08-19 23:42:34 2 NA NA 6.82
-#> 3 139017111 2005-08-19 23:48:47 2 NA NA 0.148
-#> 4 139016936 2005-08-19 23:48:34 2 NA NA 0.848
-#> 5 139018701 2005-08-19 23:50:51 2 NA NA 1.27
-#> 6 139019008 2005-08-19 23:51:25 2 NA NA 0.969
-#> 7 139023859 2005-08-19 23:58:26 2 NA NA 0.104
-#> 8 139023367 2005-08-19 23:57:37 2 NA NA 0.495
-#> 9 139023626 2005-08-19 23:58:04 2 NA NA 1.51
-#> 10 139023588 2005-08-19 23:58:02 2 NA NA 1.92
-#> # ... with 3,445 more rows, and 7 more variables: FTYPE <int>, FCODE <int>,
-#> # VisibilityFilter <int>, Shape_Length <dbl>, COMID <dbl>, VPUID <chr>,
-#> # Shape <MULTILINESTRING [°]>
-#>
-#> $NHDPoint
-#> Simple feature collection with 205 features and 10 fields
-#> Geometry type: POINT
-#> Dimension: XYZ
-#> Bounding box: xmin: -79.95278 ymin: 33.89692 xmax: -76.74611 ymax: 36.19512
-#> z_range: zmin: 0 zmax: 0
-#> Geodetic CRS: NAD83
-#> # A tibble: 205 x 11
-#> Permanent_Identif~ FDate Resolution GNIS_ID GNIS_Name REACHCODE
-#> * <chr> <dttm> <int> <chr> <chr> <chr>
-#> 1 {561E554C-9ACA-40~ 2015-03-06 04:56:16 2 NA NA NA
-#> 2 {E47148F4-D87A-46~ 2015-03-06 04:56:38 2 NA NA NA
-#> 3 {E4AFC45B-59F8-4E~ 2015-03-06 04:54:47 2 NA NA NA
-#> 4 {3383D3AE-8A45-44~ 2015-03-06 04:53:19 2 NA NA NA
-#> 5 {4E943A20-7DFB-48~ 2015-03-06 04:55:30 2 NA NA NA
-#> 6 {0D8CF1D6-1633-4D~ 2015-03-06 04:55:30 2 NA NA NA
-#> 7 {5E0C8165-0C5D-43~ 2015-03-06 04:53:19 2 NA NA NA
-#> 8 {D212825B-BD81-4D~ 2015-03-06 04:56:38 2 NA NA NA
-#> 9 {75E55D94-E988-40~ 2015-03-06 04:56:16 2 NA NA NA
-#> 10 {A8BD10FB-F9C7-4C~ 2015-03-06 04:56:16 2 NA NA NA
-#> # ... with 195 more rows, and 5 more variables: FTYPE <int>, FCODE <int>,
-#> # COMID <dbl>, VPUID <chr>, Shape <POINT [°]>
-#>
-#> $NHDWaterbody
-#> Simple feature collection with 63998 features and 15 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XYZ
-#> Bounding box: xmin: -80.06802 ymin: 33.85928 xmax: -75.48223 ymax: 36.4423
-#> z_range: zmin: 0 zmax: 0
-#> Geodetic CRS: NAD83
-#> # A tibble: 63,998 x 16
-#> Permanent_Identif~ FDate Resolution GNIS_ID GNIS_Name AreaSqKM
-#> * <chr> <dttm> <int> <chr> <chr> <dbl>
-#> 1 31414233 2002-01-09 18:00:00 2 NA NA 0.001
-#> 2 31413699 2002-01-09 18:00:00 2 NA NA 0.006
-#> 3 31425323 2003-09-18 19:00:00 2 009934~ Robertson~ 0.618
-#> 4 31412701 2002-01-09 18:00:00 2 NA NA 0.007
-#> 5 31411617 2002-01-09 18:00:00 2 NA NA 0.007
-#> 6 129759890 2004-12-02 23:01:04 2 NA NA 0.005
-#> 7 129759009 2004-12-02 22:46:36 2 NA NA 0.002
-#> 8 129758964 2004-12-02 22:45:47 2 NA NA 0.003
-#> 9 129758091 2004-12-02 22:29:04 2 NA NA 0.003
-#> 10 129756601 2004-12-02 22:02:43 2 NA NA 0.004
-#> # ... with 63,988 more rows, and 10 more variables: Elevation <dbl>,
-#> # REACHCODE <chr>, FTYPE <int>, FCODE <int>, VisibilityFilter <int>,
-#> # Shape_Length <dbl>, Shape_Area <dbl>, COMID <dbl>, VPUID <chr>,
-#> # Shape <MULTIPOLYGON [°]>
-#>
-#> $WBDHU10
-#> Simple feature collection with 127 features and 16 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -80.07334 ymin: 33.76852 xmax: -75.47986 ymax: 36.45033
-#> Geodetic CRS: NAD83
-#> # A tibble: 127 x 17
-#> TNMID MetaSourceID SourceDataDesc SourceOriginator SourceFeatureID
-#> * <chr> <chr> <chr> <chr> <chr>
-#> 1 {F9D1A928-B~ NA NA NA NA
-#> 2 {462ADC6F-D~ NA NA NA NA
-#> 3 {554A9BCF-9~ {6CBC8E85-18DD-~ NA NA NA
-#> 4 {BD302854-D~ NA NA NA NA
-#> 5 {9D5CB3CC-4~ NA NA NA NA
-#> 6 {7577A691-E~ NA NA NA NA
-#> 7 {19F9F3DE-6~ NA NA NA NA
-#> 8 {5B2ED16F-B~ NA NA NA NA
-#> 9 {12C74C5E-2~ NA NA NA NA
-#> 10 {A3E6670B-8~ NA NA NA NA
-#> # ... with 117 more rows, and 12 more variables: LoadDate <dttm>,
-#> # AreaSqKM <dbl>, AreaAcres <dbl>, GNIS_ID <int>, Name <chr>, States <chr>,
-#> # HUC10 <chr>, HUType <chr>, HUMod <chr>, Shape_Length <dbl>,
-#> # Shape_Area <dbl>, Shape <MULTIPOLYGON [°]>
-#>
-#> $WBDHU12
-#> Simple feature collection with 728 features and 21 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -80.07334 ymin: 33.76852 xmax: -75.47986 ymax: 36.45033
-#> Geodetic CRS: NAD83
-#> # A tibble: 728 x 22
-#> TNMID MetaSourceID SourceDataDesc SourceOriginator SourceFeatureID
-#> * <chr> <chr> <chr> <chr> <chr>
-#> 1 {C248C033-2~ NA NA NA NA
-#> 2 {4C5C83E6-0~ NA NA NA NA
-#> 3 {31D10A49-5~ NA NA NA NA
-#> 4 {819F9C49-9~ NA NA NA NA
-#> 5 {E7B23A13-0~ NA NA NA NA
-#> 6 {DCE3E791-1~ NA NA NA NA
-#> 7 {A7BD3900-B~ NA NA NA NA
-#> 8 {8D6652B7-D~ {F85C7E2B-B7BA-~ NA NA NA
-#> 9 {108614E0-F~ NA NA NA NA
-#> 10 {040938EA-E~ NA NA NA NA
-#> # ... with 718 more rows, and 17 more variables: LoadDate <dttm>,
-#> # NonContributingAreaAcres <dbl>, NonContributingAreaSqKm <dbl>,
-#> # AreaSqKM <dbl>, AreaAcres <dbl>, GNIS_ID <int>, Name <chr>, States <chr>,
-#> # HUC12 <chr>, HUType <chr>, HUMod <chr>, TOHUC <chr>, Shape_Length <dbl>,
-#> # Shape_Area <dbl>, COMID <dbl>, VPUID <chr>, Shape <MULTIPOLYGON [°]>
-#>
-#> $WBDHU2
-#> Simple feature collection with 2 features and 14 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -90.6235 ymin: 24.39533 xmax: -75.3981 ymax: 37.52103
-#> Geodetic CRS: NAD83
-#> # A tibble: 2 x 15
-#> TNMID MetaSourceID SourceDataDesc SourceOriginator SourceFeatureID
-#> * <chr> <chr> <chr> <chr> <chr>
-#> 1 {9C77B5BA-946A-4~ NA NA NA NA
-#> 2 {9C77B5BA-946A-4~ NA NA NA NA
-#> # ... with 10 more variables: LoadDate <dttm>, AreaSqKM <dbl>, AreaAcres <dbl>,
-#> # GNIS_ID <int>, Name <chr>, States <chr>, HUC2 <chr>, Shape_Length <dbl>,
-#> # Shape_Area <dbl>, Shape <MULTIPOLYGON [°]>
-#>
-#> $WBDHU4
-#> Simple feature collection with 2 features and 14 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -80.07334 ymin: 33.76852 xmax: -75.47986 ymax: 36.45033
-#> Geodetic CRS: NAD83
-#> # A tibble: 2 x 15
-#> TNMID MetaSourceID SourceDataDesc SourceOriginator SourceFeatureID
-#> * <chr> <chr> <chr> <chr> <chr>
-#> 1 {E3A2BA98-7EAC-4~ NA NA NA NA
-#> 2 {3F7BF113-6CB0-4~ NA NA NA NA
-#> # ... with 10 more variables: LoadDate <dttm>, AreaSqKM <dbl>, AreaAcres <dbl>,
-#> # GNIS_ID <int>, Name <chr>, States <chr>, HUC4 <chr>, Shape_Length <dbl>,
-#> # Shape_Area <dbl>, Shape <MULTIPOLYGON [°]>
-#>
-#> $WBDHU6
-#> Simple feature collection with 4 features and 14 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -80.07334 ymin: 33.76852 xmax: -75.47986 ymax: 36.45033
-#> Geodetic CRS: NAD83
-#> # A tibble: 4 x 15
-#> TNMID MetaSourceID SourceDataDesc SourceOriginator SourceFeatureID
-#> * <chr> <chr> <chr> <chr> <chr>
-#> 1 {CE8CC23B-6615-4~ NA NA NA NA
-#> 2 {B327329B-9A39-4~ NA NA NA NA
-#> 3 {98935E04-7B9A-4~ NA NA NA NA
-#> 4 {417E4E48-D599-4~ NA NA NA NA
-#> # ... with 10 more variables: LoadDate <dttm>, AreaSqKM <dbl>, AreaAcres <dbl>,
-#> # GNIS_ID <int>, Name <chr>, States <chr>, HUC6 <chr>, Shape_Length <dbl>,
-#> # Shape_Area <dbl>, Shape <MULTIPOLYGON [°]>
-#>
-#> $WBDHU8
-#> Simple feature collection with 17 features and 14 fields
-#> Geometry type: MULTIPOLYGON
-#> Dimension: XY
-#> Bounding box: xmin: -80.07334 ymin: 33.76852 xmax: -75.47986 ymax: 36.45033
-#> Geodetic CRS: NAD83
-#> # A tibble: 17 x 15
-#> TNMID MetaSourceID SourceDataDesc SourceOriginator SourceFeatureID
-#> * <chr> <chr> <chr> <chr> <chr>
-#> 1 {DA444DC2-2005-~ NA NA NA NA
-#> 2 {30075BED-99D6-~ NA NA NA NA
-#> 3 {BE9AF0AF-35BA-~ NA NA NA NA
-#> 4 {552AF5D8-63AD-~ NA NA NA NA
-#> 5 {9B302DAB-C9D6-~ NA NA NA NA
-#> 6 {90893D3B-D3DC-~ NA NA NA NA
-#> 7 {626DC6A5-FE7C-~ NA NA NA NA
-#> 8 {F520399F-E70A-~ NA NA NA NA
-#> 9 {FBD90088-665C-~ NA NA NA NA
-#> 10 {284C685F-1D0F-~ NA NA NA NA
-#> 11 {DDBE7867-F9AA-~ NA NA NA NA
-#> 12 {5C726AB4-E34B-~ NA NA NA NA
-#> 13 {7B35EF71-C30B-~ NA NA NA NA
-#> 14 {E67F756E-4BC3-~ NA NA NA NA
-#> 15 {004E87ED-EA4E-~ NA NA NA NA
-#> 16 {EE4CF6AB-BDFA-~ NA NA NA NA
-#> 17 {9FE761E6-5B59-~ NA NA NA NA
-#> # ... with 10 more variables: LoadDate <dttm>, AreaSqKM <dbl>, AreaAcres <dbl>,
-#> # GNIS_ID <int>, Name <chr>, States <chr>, HUC8 <chr>, Shape_Length <dbl>,
-#> # Shape_Area <dbl>, Shape <MULTIPOLYGON [°]>
-#>
-#> $WBDLine
-#> Simple feature collection with 2320 features and 9 fields
-#> Geometry type: MULTILINESTRING
-#> Dimension: XY
-#> Bounding box: xmin: -80.07334 ymin: 33.76852 xmax: -75.47986 ymax: 36.45033
-#> Geodetic CRS: NAD83
-#> # A tibble: 2,320 x 10
-#> TNMID HUDigit HUMod LineSource MetaSourceID LoadDate Shape_Length
-#> * <chr> <int> <chr> <chr> <chr> <dttm> <dbl>
-#> 1 {7961~ 10 NA TOPO24,DR~ NA 2012-10-15 07:56:53 0.0983
-#> 2 {3E17~ 10 NA TOPO24,DR~ NA 2012-10-15 07:56:53 0.255
-#> 3 {402D~ 10 NA TOPO24,DR~ NA 2012-10-15 07:56:42 0.228
-#> 4 {AB0B~ 10 NA TOPO24,DR~ NA 2012-10-15 07:56:42 0.109
-#> 5 {1119~ 12 NA TOPO24,DR~ NA 2012-10-15 07:56:42 0.0689
-#> 6 {46A9~ 12 NA TOPO24,DR~ NA 2012-10-15 07:56:42 0.0272
-#> 7 {990D~ 6 NA TOPO24,DR~ NA 2012-10-15 07:56:42 0.0286
-#> 8 {C4E8~ 12 NA TOPO24,DR~ NA 2012-10-15 07:56:42 0.0623
-#> 9 {F91F~ 10 NA TOPO24,DR~ NA 2012-10-15 07:56:53 0.0268
-#> 10 {5B09~ 10 NA TOPO24,DR~ NA 2012-10-15 07:56:53 0.0309
-#> # ... with 2,310 more rows, and 3 more variables: COMID <dbl>, VPUID <chr>,
-#> # Shape <MULTILINESTRING [°]>
-#>
-#> $HYDRO_NET_Junctions
-#> Simple feature collection with 205950 features and 1 field
-#> Geometry type: POINT
-#> Dimension: XY
-#> Bounding box: xmin: -80.06746 ymin: 33.8536 xmax: -75.48296 ymax: 36.44617
-#> Geodetic CRS: NAD83
-#> # A tibble: 205,950 x 2
-#> Enabled Shape
-#> * <int> <POINT [°]>
-#> 1 1 (-79.22356 36.06807)
-#> 2 1 (-79.22117 36.07732)
-#> 3 1 (-79.21987 36.07227)
-#> 4 1 (-79.21977 36.04575)
-#> 5 1 (-79.21945 36.04558)
-#> 6 1 (-79.21937 36.08515)
-#> 7 1 (-79.21504 36.06958)
-#> 8 1 (-79.21484 36.05786)
-#> 9 1 (-79.21476 36.06924)
-#> 10 1 (-79.21457 36.05783)
-#> # ... with 205,940 more rows
-#>
get_nhdplushr(download_dir,
file.path(download_dir, "nhdplus_0302-03.gpkg"),
layers = "NHDFlowline", overwrite = TRUE,
min_size_sqkm = 10, simp = 10, proj = "+init=epsg:5070")
-#> Warning: GDAL Message 1: +init=epsg:XXXX syntax is deprecated. It might return a CRS with a non-EPSG compliant axis order.
-#> $NHDFlowline
-#> Simple feature collection with 45242 features and 57 fields
-#> Geometry type: LINESTRING
-#> Dimension: XY
-#> Bounding box: xmin: 1414457 ymin: 1355556 xmax: 1795513 ymax: 1640020
-#> Projected CRS: NAD83 / Conus Albers
-#> # A tibble: 45,242 x 58
-#> Permanent_Identif~ FDate Resolution GNIS_ID GNIS_Name LENGTHKM
-#> * <chr> <dttm> <int> <chr> <chr> <dbl>
-#> 1 31374279 2012-03-10 19:59:53 2 009821~ Buffalo C~ 0.029
-#> 2 31373805 2012-03-10 19:59:53 2 009821~ Buffalo C~ 0.092
-#> 3 31368135 2012-03-10 21:01:26 2 009920~ Perry Cre~ 0.513
-#> 4 129726996 2012-03-10 16:04:41 2 010239~ Crooked C~ 0.74
-#> 5 129744832 2012-03-10 16:03:34 2 010241~ Lynch Cre~ 0.059
-#> 6 129725633 2012-03-10 16:00:22 2 NA NA 1.44
-#> 7 129723709 2012-03-10 16:03:44 2 009814~ Billys Cr~ 0.889
-#> 8 31402583 2012-03-10 19:58:45 2 NA NA 0.371
-#> 9 129721510 2012-03-10 16:04:01 2 010238~ Buffalo C~ 0.39
-#> 10 129720736 2012-03-10 16:01:17 2 010238~ Buffalo C~ 0.116
-#> # ... with 45,232 more rows, and 52 more variables: REACHCODE <chr>,
-#> # FlowDir <int>, WBArea_Permanent_Identifier <chr>, FTYPE <int>, FCODE <int>,
-#> # MainPath <int>, InNetwork <int>, VisibilityFilter <int>,
-#> # Shape_Length <dbl>, COMID <dbl>, VPUID <chr>, Enabled <int>,
-#> # Shape <LINESTRING [m]>, StreamLeve <int>, StreamOrde <int>,
-#> # StreamCalc <int>, FromNode <dbl>, ToNode <dbl>, Hydroseq <dbl>,
-#> # LevelPathI <dbl>, Pathlength <dbl>, TerminalPa <dbl>, ArbolateSu <dbl>,
-#> # Divergence <int>, StartFlag <int>, TerminalFl <dbl>, UpLevelPat <dbl>,
-#> # UpHydroSeq <dbl>, DnLevel <dbl>, DnLevelPat <dbl>, DnHydroseq <dbl>,
-#> # DnMinorHyd <dbl>, DnDrainCou <int>, FromMeas <dbl>, ToMeas <dbl>,
-#> # RtnDiv <int>, Thinner <int>, VPUIn <int>, VPUOut <int>, AreaSqKM <dbl>,
-#> # TotDASqKM <dbl>, DivDASqKm <dbl>, MaxElevRaw <dbl>, MinElevRaw <dbl>,
-#> # MaxElevSmo <dbl>, MinElevSmo <dbl>, Slope <dbl>, SlopeLenKm <dbl>,
-#> # ElevFixed <int>, HWType <int>, HWNodeSqKm <dbl>, StatusFlag <chr>
-#>
# Cleanup
unlink(temp_dir, recursive = TRUE)
-# }
+}
diff --git a/docs/reference/nhdplus_path.html b/docs/reference/nhdplus_path.html
index 2aac1f9b..d2ff09c3 100644
--- a/docs/reference/nhdplus_path.html
+++ b/docs/reference/nhdplus_path.html
@@ -8,13 +8,6 @@
NHDPlus Data Path — nhdplus_path • nhdplusTools
-
-
-
-
-
-
-
diff --git a/docs/reference/plot_nhdplus.html b/docs/reference/plot_nhdplus.html
index 6151978a..3d6ba6b7 100644
--- a/docs/reference/plot_nhdplus.html
+++ b/docs/reference/plot_nhdplus.html
@@ -8,13 +8,6 @@
Plot NHDPlus — plot_nhdplus • nhdplusTools
-
-
-
-
-
-
-
@@ -257,12 +250,6 @@ Examp
options("rgdal_show_exportToProj4_warnings"="none")
plot_nhdplus("05428500")
-#> Request failed [503]. Retrying in 1 seconds...
-#> Request failed [503]. Retrying in 3 seconds...
-#> Request failed [503]. Retrying in 2.4 seconds...
-#> Request failed [503]. Retrying in 6.1 seconds...
-#> Request failed [503]. Retrying in 5.8 seconds...
-#> Request failed [503]. Retrying in 40.8 seconds...
#> although coordinates are longitude/latitude, st_intersects assumes that they are planar
#> Zoom: 11
#> Map tiles by Carto, under CC BY 3.0. Data by OpenStreetMap, under ODbL.
@@ -377,8 +364,9 @@ Examp
#> Audotdetect projection: assuming Google Mercator (epsg 3857)
-# Cleanup downloaded open street map cache if desired.
-unlink(nhdplusTools:::osm_cache_dir(), recursive = TRUE)
+# Cleanup downloaded open street map cache and temp data dir if desired.
+# This is included for CRAN checks primarily.
+unlink(nhdplusTools::nhdplusTools_data_dir(), recursive = TRUE)
# }