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

@@ -17,13 +17,19 @@ 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$lightblue,
+       border = NA)
+  plot(sf::st_geometry(political$geometry), 
        col = colors$lightgrey,
-     border = "black")
-plot(political$Wisconsin_State_Boundary_24K.geojson, 
+       border = NA,
+       add = TRUE)
+  plot(sf::st_geometry(political[political$NAME_En %in% state, ]), 
        col = colors$darkblue,
        border = NA,
        add = TRUE)
@@ -40,6 +46,11 @@ plot(sf::st_geometry(data$NHDFlowline),
        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

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