2023-09-27_GES668_in-class-script.R

library(sf)
library(dplyr)

# Can you convert `storms_sf` to make the feature geometry represent
# storm tracks instead of storm observations?
  
storms_sf <- storms |> 
  st_as_sf(
    coords = c("long", "lat"),
    crs = 4326
  )

storms_sf <- storms_sf |> 
  st_transform(crs = 4269)

storm_tracks <- storms_sf |>
  arrange(year, month, day, hour) |> 
  group_by(year, name) |> 
  summarise(
    avg_wind = mean(wind),
    geometry = st_cast(st_union(geometry), to = "LINESTRING")
  ) |> 
  ungroup()

md_storm_tracks <- st_filter(storm_tracks, maryland)

storm_tracks <- smoothr::smooth(storm_tracks)

md_bbox <- st_bbox(maryland)

md_storm_map <- md_storm_tracks |> 
  ggplot() +
  geom_sf(aes(color = avg_wind)) +
  geom_sf(data = lower48, fill = NA, color = "black") +
  coord_sf(
    xlim = c(md_bbox$xmin, md_bbox$xmax),
    ylim = c(md_bbox$ymin, md_bbox$ymax)
  ) +
  theme_void()