From f755308d145e69b0405fa791e96a5ed60a703123 Mon Sep 17 00:00:00 2001 From: Ben Varick Date: Mon, 16 Oct 2023 16:59:19 -0500 Subject: [PATCH] changed spacing/formatting of scripts --- 01_download_data.R | 1 - 02_process_data.R | 26 +++++++++++++------------- 2 files changed, 13 insertions(+), 14 deletions(-) diff --git a/01_download_data.R b/01_download_data.R index 72843a3..bdd4de3 100644 --- a/01_download_data.R +++ b/01_download_data.R @@ -30,7 +30,6 @@ extent <- list(longitude_max = max(extent$longitude), latitude_max = max(extent$latitude), latitude_min = min(extent$latitude)) - # get watershed areas extent_huc <- get_huc(AOI = extent_poly, buffer = 0, type = "huc04") diff --git a/02_process_data.R b/02_process_data.R index 0a06587..ce84362 100644 --- a/02_process_data.R +++ b/02_process_data.R @@ -10,34 +10,34 @@ extent <- read.csv(file = "extent.csv") crs <- 4269 extent_poly <- st_polygon( -x = list( - cbind( - extent$longitude[c(1,2,2,1,1)], - extent$latitude[c(1,1,2,2,1)]) -) + 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) extent_bbox <- st_bbox(st_transform(x = extent_poly, crs)) extent <- list(longitude_max = max(extent$longitude), - longitude_min = min(extent$longitude), - latitude_max = max(extent$latitude), - latitude_min = min(extent$latitude)) + longitude_min = min(extent$longitude), + latitude_max = max(extent$latitude), + latitude_min = min(extent$latitude)) # load and crop hydrologic data data <- list(NULL) sf_use_s2(FALSE) layers <- c("NHDArea", - "NHDFlowline", - "NHDWaterbody", - "NHDPlusLandSea") + "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)) }