diff --git a/02_process_data.R b/02_process_data.R index ce84362..3418dab 100644 --- a/02_process_data.R +++ b/02_process_data.R @@ -36,8 +36,9 @@ for (layer in layers){ } # 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 <- st_crop( + st_transform( + st_read( + paste0(data_dir, "/political_boundaries/bound_p/boundaries_p_2021_v3.shp")), + crs=crs), + y = extent_bbox) \ No newline at end of file diff --git a/03_make_figures.R b/03_make_figures.R index 097bb63..f6e03b5 100644 --- a/03_make_figures.R +++ b/03_make_figures.R @@ -7,8 +7,8 @@ ifelse(!dir.exists(file.path(getwd(), figure_dir)), dir.create(file.path(getwd() # set colors colors <- list(darkblue = "#062e57", - lightblue = "#b1dcf3", - lightgrey = "#6D6D6D") + lightblue = "#b1dcf3", + lightgrey = "#6D6D6D") # set dimensions of figures width <- 8 @@ -17,29 +17,40 @@ units <- "in" res <- 600 margin <- 0 +# set state to highlight +state <- c("Wisconsin") # plot map plot_map <- function () { -plot(sf::st_geometry(extent_poly), - col = colors$lightgrey, - border = "black") -plot(political$Wisconsin_State_Boundary_24K.geojson, - col = colors$darkblue, - border = NA, - add = TRUE) -plot(sf::st_geometry(data$NHDWaterbody), - col = colors$lightblue, - border = NA, - add = TRUE) -plot(sf::st_geometry(data$NHDArea), - col = colors$lightblue, - border = NA, - add = TRUE) -plot(sf::st_geometry(data$NHDFlowline), - col = colors$lightblue, - lwd = data$NHDFlowline$TotDASqKM^0.3204*0.0446, - border = NA, - add = TRUE) + plot(sf::st_geometry(extent_poly), + col = colors$lightblue, + border = NA) + plot(sf::st_geometry(political$geometry), + col = colors$lightgrey, + border = NA, + add = TRUE) + plot(sf::st_geometry(political[political$NAME_En %in% state, ]), + col = colors$darkblue, + border = NA, + add = TRUE) + plot(sf::st_geometry(data$NHDWaterbody), + col = colors$lightblue, + border = NA, + add = TRUE) + plot(sf::st_geometry(data$NHDArea), + col = colors$lightblue, + border = NA, + add = TRUE) + plot(sf::st_geometry(data$NHDFlowline), + col = colors$lightblue, + lwd = data$NHDFlowline$TotDASqKM^0.3204*0.0446, + border = NA, + add = TRUE) + plot(sf::st_geometry(extent_poly), + col = NA, + border = "black", + lwd = 1, + add = TRUE) } # save figure @@ -66,10 +77,10 @@ plot(sf::st_geometry(data$NHDFlowline), # dev.off() png(filename = paste0(figure_dir, "/map.png"), - width = width, - height = height, - units = units, - res = res) + width = width, + height = height, + units = units, + res = res) par(mai=c(margin, margin, margin, margin)) par(mar=c(1,1,1,1)) plot_map() diff --git a/extent_olympic.csv b/extent_olympic.csv new file mode 100644 index 0000000..e5fc94e --- /dev/null +++ b/extent_olympic.csv @@ -0,0 +1,3 @@ +latitude,longitude +48.50775,-124.91972 +46.88398,-122.07015 diff --git a/extent_wa.csv b/extent_wa.csv new file mode 100644 index 0000000..9ab3c21 --- /dev/null +++ b/extent_wa.csv @@ -0,0 +1,3 @@ +latitude,longitude +49.17657,-125.43151 +45.23644,-115.81455 diff --git a/extent_wi.csv b/extent_wi.csv new file mode 100644 index 0000000..a630d20 --- /dev/null +++ b/extent_wi.csv @@ -0,0 +1,3 @@ +latitude,longitude +47.73208,-93.55154 +41.49273,-86.64896