library(wikilake)
library(sf)
library(nhdR)
library(dplyr)
library(magrittr)
# ---- pull_geom ----
# wb <- nhd_plus_load(vpu = 7, dsn = "NHDWaterbody")
# fl <- nhd_plus_load(vpu = 7, dsn = "NHDFlowLine")
sunapee_meta <- wikilake::lake_wiki("Lake Sunapee")
sunapee <- st_sfc(st_point(c(sunapee_meta$Lon, sunapee_meta$Lat)))
sunapee <- nhd_plus_query(lon = sunapee_meta$Lon, lat = sunapee_meta$Lat,
dsn = c("NHDWaterbody", "NHDFlowLine"), buffer_dist = 0.2)
vpu <- find_vpu(sunapee$pnt)
# ---- pull_sunapee_network ----
sunapee_network <- extract_network(lon = sunapee_meta$Lon, lat = sunapee_meta$Lat,
maxsteps = Inf)
usethis::use_data(sunapee, overwrite = TRUE)
usethis::use_data(sunapee_network, overwrite = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.