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

added layers variable to 2nd script

Browse Source
This commit is contained in:
Ben Varick 2023-10-16 16:46:49 -05:00
parent 3b77f65b7b
commit 27ada9f3bb
Signed by: ben
SSH Key Fingerprint: SHA256:jWnpFDAcacYM5aPFpYRqlsamlDyKNpSj3jj+k4ojtUo
1 changed files with 15 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
02_process_data.R
View File

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