pretty_rivers/02_process_data.R
Ben Varick 67c794c419
added political boundaries
adjusted data dir
added steps and example to readme
2023-10-18 14:11:35 -05:00

17 lines
487 B
R

# load libraries
library(sf)
# load and crop hydrologic data
data <- list(NULL)
sf_use_s2(FALSE)
for (layer in layers){
data[[layer]] <- st_crop(st_read(paste0(nhdplusTools_data_dir(),"/data.gpkg"), layer = layer), y = extent_bbox)
}
# load political boundaries
political <- list(NULL)
political_boundaries <- list.files(path = "data/political_boundaries")
for (boundary in political_boundaries) {
political[[boundary]] <- st_read(paste0("data/political_boundaries/", boundary))
}