added layers variable to 2nd script

02_process_data.R

@@ -1,4 +1,5 @@
 # load libraries
+library(nhdplusTools)
 library(sf)
 
 # reload extent data so you don't need to run the download data script again
@@ -9,11 +10,11 @@ extent <- read.csv(file = "extent.csv")
 crs <- 4269
 
 extent_poly <- st_polygon(
-  x = list(
+x = list(
   cbind(
     extent$longitude[c(1,2,2,1,1)], 
     extent$latitude[c(1,1,2,2,1)])
-  )
+)
 )
 
 extent_poly <- st_sfc(extent_poly, crs=4326)
@@ -26,13 +27,17 @@ extent <- list(longitude_max = max(extent$longitude),
 # load and crop hydrologic data
 data <- list(NULL)
 sf_use_s2(FALSE)
+layers <- c("NHDArea",
+          "NHDFlowline",
+          "NHDWaterbody",
+          "NHDPlusLandSea")
 for (layer in layers){
-  data[[layer]] <- st_crop(st_read(paste0(nhdplusTools_data_dir(),"/data.gpkg"), layer = layer), y = extent_bbox)
+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 = paste0(data_dir, "/political_boundaries"))
 for (boundary in political_boundaries) {
-  political[[boundary]] <- st_read(paste0(data_dir, "/political_boundaries/", boundary))
+political[[boundary]] <- st_read(paste0(data_dir, "/political_boundaries/", boundary))
 }