2023-10-16 12:17:24 -05:00
|
|
|
# load libraries
|
|
|
|
library(sf)
|
|
|
|
|
2023-10-16 13:39:07 -05:00
|
|
|
# load and crop hydrologic data
|
2023-10-16 12:17:24 -05:00
|
|
|
data <- list(NULL)
|
|
|
|
sf_use_s2(FALSE)
|
|
|
|
for (layer in layers){
|
2023-10-16 12:53:23 -05:00
|
|
|
data[[layer]] <- st_crop(st_read(paste0(nhdplusTools_data_dir(),"/data.gpkg"), layer = layer), y = extent_bbox)
|
2023-10-16 12:17:24 -05:00
|
|
|
}
|
2023-10-16 13:39:07 -05:00
|
|
|
|
|
|
|
# load political boundaries
|
|
|
|
political <- list(NULL)
|
2023-10-16 13:59:59 -05:00
|
|
|
political_boundaries <- list.files(path = paste0(data_dir, "/political_boundaries"))
|
2023-10-16 13:39:07 -05:00
|
|
|
for (boundary in political_boundaries) {
|
2023-10-16 13:59:59 -05:00
|
|
|
political[[boundary]] <- st_read(paste0(data_dir, "/political_boundaries/", boundary))
|
2023-10-16 13:39:07 -05:00
|
|
|
}
|