changed political boundary data

added a way to select state in 03_make_figures.R

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)

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()

extent_olympic.csv

@@ -0,0 +1,3 @@
+latitude,longitude
+48.50775,-124.91972
+46.88398,-122.07015

extent_wa.csv

@@ -0,0 +1,3 @@
+latitude,longitude
+49.17657,-125.43151
+45.23644,-115.81455

extent_wi.csv

@@ -0,0 +1,3 @@
+latitude,longitude
+47.73208,-93.55154
+41.49273,-86.64896