ben/pretty_rivers
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

changed spacing/formatting of scripts

Browse Source
This commit is contained in:
Ben Varick 2023-10-16 16:59:19 -05:00
parent 4b33f2d71a
commit 4e4576d93d
2 changed files with 13 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
01_download_data.R
View File

@ -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")

8
02_process_data.R
View File

@ -10,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)
@ -32,12 +32,12 @@ layers <- c("NHDArea",
"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))
}