This is more of an example than anything
Arguments
- x
The result of get_nhd.
- template
An
Simple FeatureorSpatRasterobject to serve as a template for cropping.
Examples
# \donttest{
# Get the NHD (USA ONLY)
NHD <- get_nhd(
template = FedData::meve,
label = "meve"
)
NHD
#> $Point
#> Simple feature collection with 2 features and 10 fields
#> Geometry type: POINT
#> Dimension: XY
#> Bounding box: xmin: -12083250 ymin: 4469345 xmax: -12062340 ymax: 4474361
#> Projected CRS: WGS 84 / Pseudo-Mercator
#> # A tibble: 2 × 11
#> OBJECTID FDATE RESOLUTION GNIS_ID GNIS_NAME REACHCODE FTYPE
#> <dbl> <dttm> <dbl> <chr> <chr> <chr> <dbl>
#> 1 16724 2003-07-02 00:00:00 2 00179195 Rock Springs NA 458
#> 2 481568 2015-03-06 10:56:39 2 NA NA NA 367
#> # ℹ 4 more variables: FCODE <dbl>, PERMANENT_IDENTIFIER <chr>, GLOBALID <chr>,
#> # geom <POINT [m]>
#>
#> $Flowline
#> Simple feature collection with 424 features and 17 fields
#> Geometry type: MULTILINESTRING
#> Dimension: XY
#> Bounding box: xmin: -12084260 ymin: 4460933 xmax: -12060340 ymax: 4487920
#> Projected CRS: WGS 84 / Pseudo-Mercator
#> # A tibble: 424 × 18
#> OBJECTID permanent_identifier fdate resolution gnis_id
#> <dbl> <chr> <dttm> <dbl> <chr>
#> 1 25888 89636109 2012-02-18 12:02:41 2 NA
#> 2 25988 89636095 2012-02-18 12:00:10 2 00902295
#> 3 143005 89639153 2015-09-08 00:00:00 2 00902295
#> 4 471453 89639251 2012-02-18 12:01:07 2 00902295
#> 5 472018 89636103 2012-02-18 12:00:10 2 00902295
#> 6 472099 89638687 2012-02-18 12:01:27 2 NA
#> 7 586689 89635471 2012-02-18 12:00:33 2 NA
#> 8 586796 89638705 2012-02-18 12:00:14 2 NA
#> 9 809460 89634465 2015-07-01 00:00:00 2 NA
#> 10 887363 82523370 2020-10-09 00:00:00 2 NA
#> # ℹ 414 more rows
#> # ℹ 13 more variables: gnis_name <chr>, lengthkm <dbl>, reachcode <chr>,
#> # flowdir <dbl>, wbarea_permanent_identifier <chr>, ftype <dbl>, fcode <dbl>,
#> # innetwork <dbl>, mainpath <dbl>, visibilityfilter <dbl>,
#> # Shape_Length <dbl>, globalid <chr>, geom <MULTILINESTRING [m]>
#>
#> $Waterbody
#> Simple feature collection with 10 features and 15 fields
#> Geometry type: POLYGON
#> Dimension: XY
#> Bounding box: xmin: -12081810 ymin: 4465756 xmax: -12067270 ymax: 4484894
#> Projected CRS: WGS 84 / Pseudo-Mercator
#> # A tibble: 10 × 16
#> OBJECTID PERMANENT_IDENTIFIER FDATE RESOLUTION GNIS_ID
#> <dbl> <chr> <dttm> <dbl> <chr>
#> 1 1326666 89638655 2003-07-02 00:00:00 2 NA
#> 2 2436074 89634565 2003-07-02 00:00:00 2 NA
#> 3 2532090 89638207 2003-07-02 00:00:00 2 NA
#> 4 2893566 89634585 2003-07-02 00:00:00 2 NA
#> 5 3640339 89634605 2003-07-02 00:00:00 2 NA
#> 6 4339039 89634591 2003-07-02 00:00:00 2 NA
#> 7 4459618 89634587 2003-07-02 00:00:00 2 NA
#> 8 4821810 89638209 2003-07-02 00:00:00 2 NA
#> 9 5304625 89634583 2003-07-02 00:00:00 2 NA
#> 10 5906925 89634589 2003-07-02 00:00:00 2 NA
#> # ℹ 11 more variables: GNIS_NAME <chr>, AREASQKM <dbl>, ELEVATION <dbl>,
#> # REACHCODE <chr>, FTYPE <dbl>, FCODE <dbl>, VISIBILITYFILTER <dbl>,
#> # SHAPE_Length <dbl>, SHAPE_Area <dbl>, GLOBALID <chr>, geom <POLYGON [m]>
#>
NHD %>%
plot_nhd(template = FedData::meve)
# }