Merge branch 'bvarick:main' into sgy

Makefile

@@ -1,5 +1,5 @@
 route_analysis: R/route_analysis.Rmd
-	R -e 'library("rmarkdown"); old_path <- Sys.getenv("PATH"); Sys.setenv(PATH = paste(old_path, "/usr/local/bin", sep = ":")); rmarkdown::render(knit_root_dir = "~/route_analysis/", output_dir = "~/route_analysis/html", input = "./R/route_analysis.Rmd", output_file = "./html/route_analysis.html")'
+	R -e 'library("rmarkdown"); old_path <- Sys.getenv("PATH"); Sys.setenv(PATH = paste(old_path, "/usr/local/bin", sep = ":")); rmarkdown::render(knit_root_dir = "./", output_dir = "./html", input = "./R/route_analysis.Rmd", output_file = "./html/route_analysis.html")'
 
 clean: clean-data clean-figure clean-script
 

R/route_analysis.Rmd

@@ -25,7 +25,7 @@ set.seed(1)
 
 # Main R script
 
-```{r Rscript, eval = FALSE, echo = TRUE, results = "show", warning = FALSE, error = FALSE, message = FALSE}
+```{r Rscript, eval = TRUE, echo = TRUE, results = "show", warning = FALSE, error = TRUE, message = FALSE}
 ## school focus
 school_focus <- data.frame(name = c("East High School"), NCES_CODE = c("550852000925"))
 
@@ -33,10 +33,10 @@ school_focus <- data.frame(name = c("East High School"), NCES_CODE = c("55085200
 walk_boundary_m <- 1.5 * 1609
 
 ## load school locations
-WI_schools <- st_read(dsn = "data/Schools/WI_schools.gpkg")
+WI_schools <- st_read(dsn = "../data/Schools/WI_schools.gpkg")
 
 ## load addresses
-addresses <- read_csv(file="data/addresses/Addresses_Students_EastHS_2024_GeocodeResults.csv") %>%
+addresses <- read_csv(file="../data/addresses/Addresses_Students_EastHS_2024_GeocodeResults.csv") %>%
   filter(lat > 0) %>%
   st_as_sf(coords=c("lon","lat"), crs=4326) # remember x=lon and y=lat
 
@@ -44,7 +44,7 @@ addresses <- read_csv(file="data/addresses/Addresses_Students_EastHS_2024_Geocod
 options(osrm.server = "http://127.0.0.1:5000/")
 options(osrm.profile = "walk")
 
-register_stadiamaps(key = substr(read_file(file = "api_keys/stadia_api_key"), 1, 36))
+register_stadiamaps(key = substr(read_file(file = "../api_keys/stadia_api_key"), 1, 36))
 
 ## subset addresses within 1.5 miles
 walk_boundary_poly <- fill_holes(st_make_valid(osrmIsodistance(
@@ -56,7 +56,7 @@ walk_boundary_poly <- fill_holes(st_make_valid(osrmIsodistance(
 addresses_near <- st_intersection(addresses, walk_boundary_poly)
 
 ## load bike tls
-bike_lts <- st_read("data/bike_lts/bike_lts_DANE.geojson")
+bike_lts <- st_read("../data/bike_lts/bike_lts_DANE.geojson")
 bike_lts[["lts"]] <- as.factor(bike_lts$LTS_F)
 
 bike_lts_scale <- data.frame(code = c(1, 2, 3, 4, 9),
@@ -77,7 +77,7 @@ for(i in addresses_near$number) {
 routes <- bind_rows(routes)
 
 ## combine routes
-bike_lts_buffer <- st_buffer(st_intersection(bike_lts, walk_boundary_poly), 10)
+bike_lts_buffer <- st_buffer(st_intersection(bike_lts, walk_boundary_poly), 20)
 
 bike_lts_buffer["student_use"] <- unlist(lapply(st_intersects(bike_lts_buffer, routes), length))
 
@@ -85,8 +85,8 @@ bike_lts <- st_join(bike_lts, bike_lts_buffer %>% select(OBJECTID, student_use))
 
 ## make maps
 # load logo
-logo <- image_read(path = "other/BFW_Logo_180_x_200_transparent_background.png")
-school_symbol <- image_read_svg(path = "other/school_FILL0_wght400_GRAD0_opsz24.svg")
+logo <- image_read(path = "../other/BFW_Logo_180_x_200_transparent_background.png")
+school_symbol <- image_read_svg(path = "../other/school_FILL0_wght400_GRAD0_opsz24.svg")
 
 
 bbox <- st_bbox(st_transform(st_buffer(addresses_near, dist = 500), crs = 4326))
@@ -106,7 +106,7 @@ ggmap(basemap) +
     x = NULL,
     y = NULL,
     color = NULL,
-    linewidth = "How many students can use road") +
+    linewidth = "Potential student walkers") +
   theme(axis.text=element_blank(),
         axis.ticks=element_blank(),
         plot.caption = element_text(color = "grey")) +
@@ -136,7 +136,7 @@ ggmap(basemap) +
                 label.size = 0.04,
                 size = 2)
 
-ggsave(file = paste0("figures/",
+ggsave(file = paste0("../figures/",
                      school_focus %>% pull(name),
                      " Routes.pdf"),
        title = paste0(school_focus %>% pull(name), " Walking Routes"),
@@ -150,11 +150,11 @@ ggsave(file = paste0("figures/",
 ggmap(basemap) +
   labs(title = paste0("Walking routes for students at ",
                       school_focus %>% pull(name)),
-       subtitle = "only showing routes within the 1.5 walk boundary",
+       subtitle = "only showing routes within the walk boundary",
        x = NULL,
        y = NULL,
        color = NULL,
-       linewidth = "How many students can use road") +
+       linewidth = "Potential student walkers") +
   theme(axis.text=element_blank(),
         axis.ticks=element_blank(),
         plot.caption = element_text(color = "grey")) +
@@ -183,7 +183,7 @@ ggmap(basemap) +
                 nudge_y = 0.0015,
                 label.size = 0.04,
                 size = 2)
-ggsave(file = paste0("figures/",
+ggsave(file = paste0("../figures/",
                      school_focus %>% pull(name),
                      " Routes - Traffic Stress.pdf"),
        title = paste0(school_focus %>% pull(name), " Walking Routes - Traffic Stress"),
@@ -227,7 +227,7 @@ ggmap(basemap) +
                 nudge_y = 0.0015,
                 label.size = 0.04,
                 size = 2)
-ggsave(file = paste0("figures/",
+ggsave(file = paste0("../figures/",
                      school_focus %>% pull(name),
                      " Addresses.pdf"),
        title = paste0(school_focus %>% pull(name), " Addresses"),