From 51806c81d477cf9394900c750e69e9befa409a85 Mon Sep 17 00:00:00 2001
From: Ben Varick <ben@dendroalsia.net>
Date: Thu, 21 Nov 2024 17:33:44 -0600
Subject: [PATCH] changed from segments to lat lon buckets

---
 madison-metro.R | 34 ++++++++++++++++++++++++++++++----
 1 file changed, 30 insertions(+), 4 deletions(-)

diff --git a/madison-metro.R b/madison-metro.R
index ea1ad3d..c33f3db 100644
--- a/madison-metro.R
+++ b/madison-metro.R
@@ -146,11 +146,17 @@ for (route in unique(routes_categorized$name)){
     theme(axis.text=element_blank(),
           axis.ticks=element_blank(),
           plot.caption = element_text(color = "grey")) +
-    geom_sf(data = metro_segments %>% filter(name %in% route),
+    geom_tile(data = metro_summary %>% 
+                filter(name %in% route) %>% 
+                group_by(lon_bucket, lat_bucket) %>% 
+                summarise(spd_calc = weighted.mean(spd_calc, trip_count)),
             inherit.aes = FALSE,
-            aes(color = spd_calc),
-            linewidth = 1) +
-    scale_color_distiller(palette = "RdYlGn", direction = "reverse", limits = c(0,70), name = "Average speed or segment\n(calculated with locations, not reported speed)")
+            aes(x = lon_bucket,
+                y = lat_bucket,
+                fill = spd_calc,
+                height = lat_round,
+                width = lon_round)) +
+    scale_fill_distiller(palette = "RdYlGn", direction = "reverse", limits = c(0,70), name = "Average speed or segment\n(calculated with locations, not reported speed)")
   ggsave(file = paste0("figures/",
                        route,
                        "_map.pdf"),
@@ -175,3 +181,23 @@ for (route in unique(routes_categorized$name)){
          units = "in",
          create.dir = TRUE)
 }
+
+
+# ggmap(basemap) +
+#   labs(title = paste0("Metro Route Speed - ", route),
+#        subtitle = paste0("averaged between ",
+#                          sum(route_counts %>% filter(pid %in% route_focus) %>% pull(route_count)),
+#                          " bus trips - ",
+#                          min(date(metro_data$time)),
+#                          " to ",
+#                          max(date(metro_data$time))),
+#        x = NULL,
+#        y = NULL) +
+#   theme(axis.text=element_blank(),
+#         axis.ticks=element_blank(),
+#         plot.caption = element_text(color = "grey")) +
+#   geom_sf(data = metro_segments %>% filter(name %in% route),
+#           inherit.aes = FALSE,
+#           aes(color = spd_calc),
+#           linewidth = 1) +
+#   scale_color_distiller(palette = "RdYlGn", direction = "reverse", limits = c(0,70), name = "Average speed or segment\n(calculated with locations, not reported speed)")