rogers_centre_home_runs.R

###########################################################################################-
###########################################################################################-
##
## Rogers Centre Home Runs ----
##
###########################################################################################-
###########################################################################################-

#=========================================================================================#
# Setting up ----
#=========================================================================================#

#-----------------------------------------------------------------------------------------#
# Loading libraries
#-----------------------------------------------------------------------------------------#

library(tidyverse)
library(lubridate)
library(magrittr)
library(dbplyr)
library(DBI)
library(RSQLite)
library(glue)

#-----------------------------------------------------------------------------------------#
# Loading Nationals Park trace
#-----------------------------------------------------------------------------------------#

rogers.centre.df <- read_csv("data/rogers.centre.trace.df.csv")

#-----------------------------------------------------------------------------------------#
# Connecting to database
#-----------------------------------------------------------------------------------------#

statcast_db <- dbConnect(SQLite(), "data/statcast_db_rebuilt.sqlite3")

#-----------------------------------------------------------------------------------------#
# Pulling data from database ----
#-----------------------------------------------------------------------------------------#

rogers_centre_data <-
    tbl(statcast_db, "statcast_data") %>%
    filter(ballpark == "Rogers Centre" & events == "home_run") %>%
    filter(!is.na(events)) %>% 
    select(game_date, 
           game_year,
           bb_type,
           events,
           player_name,
           batter,
           hit_distance_sc,
           launch_speed,
           horiz_angle,
           horiz_angle_deg
    ) %>% 
    collect() %>% 
    mutate(
        horiz_angle_3 = case_when(
            horiz_angle_deg >= 0  & horiz_angle_deg <  30 ~ "Right",
            horiz_angle_deg >= 30 & horiz_angle_deg <= 60 ~ "Center",
            horiz_angle_deg >  60 & horiz_angle_deg <= 90 ~ "Left"
        )
    )

most_recent_day <-
    rogers_centre_data %>%
    select(game_date) %>%
    arrange(desc(game_date)) %>%
    head(1) %>%
    pull(game_date) %>% 
    as_date()


#=========================================================================================#
# Plotting ----
#=========================================================================================#

diamond.plot.field.RC <-
    rogers_centre_data %>% 
    ggplot(
        aes(
            x = 0 + hit_distance_sc * cos(horiz_angle), 
            y = 0 + hit_distance_sc * sin(horiz_angle))) +
    
    coord_fixed(xlim = c(0, 500), ylim = c(0, 500), expand = TRUE) +
    
    labs(
        title = "Home Runs @ The Rogers Centre",
        subtitle = glue(
            "{nrow(rogers_centre_data)} home runs",
            " ({rogers_centre_data %>% drop_na(hit_distance_sc, horiz_angle) %>% nrow()} plotted)",
            " as of {format(most_recent_day, '%m/%d/%Y')}")
    ) +
    
    
    ## Dividing the fields ##
    
    ### Left ###
    
    annotate(
        x = c(
            0,
            (0 + 550 * cos((60:90 * pi) / 180)),
            0),
        y = c(
            0,
            (0 + 550 * sin((60:90 * pi) / 180)),
            0),
        geom = "polygon",
        fill = "white",
        alpha = 1) +
    
    ### Center ###
    
    annotate(
        x = c(
            0,
            (0 + 550 * cos((30:60 * pi) / 180)),
            0), 
        y = c(
            0,
            (0 + 550 * sin((30:60 * pi) / 180)),
            0), 
        geom = "polygon",
        fill = "ivory2",
        alpha = 1) +
    
    
    ### Right ###
    
    annotate(
        x = c(
            0,
            (0 + 550 * cos((0:30 * pi) / 180)),
            0),
        y = c(
            0,
            (0 + 550 * sin((0:30 * pi) / 180)),
            0),
        geom = "polygon",
        fill = "white",
        alpha = 1) +
    
    
    ## Grass ##
    
    annotate(
        x = rogers.centre.df$x,
        y = rogers.centre.df$y,
        geom = "polygon",
        fill = "#9ACB9A",
        alpha = 1) +
    
    
    ## Infield dirt ##
    
    annotate(
        x = c(
            (59 * cos((45 * pi) / 180)) + 95 * cos((-26:116 * pi) / 180),
            0,
            90,
            90),
        y = c(
            (59 * cos((45 * pi) / 180)) + 95 * sin((-26:116 * pi) / 180),
            90,
            90,
            0),
        geom = "polygon",
        fill = "#DBAF84",
        alpha = 1) +
    
    ## Infield grass ##
    
    annotate(
        x = c(90, 90, 0, 0), 
        y = c(0, 90, 90, 0),
        geom = "polygon", 
        fill = "#9ACB9A",
        alpha = 1) + 
    
    ## Pitcher's mound dirt ##
    
    annotate(
        x = (59 * cos((45 * pi) / 180)) + (10 * cos((seq(0, 360, .1) * pi) / 180)),
        y = (59 * sin((45 * pi) / 180)) + (10 * sin((seq(0, 360, .1) * pi) / 180)),
        geom = "polygon",
        fill = "#DBAF84",
        alpha = 1) +
    
    ## Pitcher's mound cut ##
    
    annotate(
        x = (59 * cos((45 * pi) / 180)) + (10 * cos((seq(0, 360, .1) * pi) / 180)),
        y = (59 * sin((45 * pi) / 180)) + (10 * sin((seq(0, 360, .1) * pi) / 180)),
        geom = "path",
        size = .25,
        color = "gray20"
    ) +
    
    ## Pitcher's mound rubber ##
    
    annotate(
        geom = "text",
        label = "-",
        x = (58 * cos((45 * pi) / 180)),
        y = (58 * sin((45 * pi) / 180)), 
        color = "white",
        size = 8,
        angle = 135
    ) +
    
    ## Wall ##
    
    annotate(
        geom = "path",
        x = rogers.centre.df$x,
        y = rogers.centre.df$y,
        color = "black",
        size = 1
    ) +
    
    ## Grass cut ##
    
    annotate(
        geom = "path",
        x = (59 * cos((45 * pi) / 180)) + 95 * cos((-26:116 * pi) / 180),
        y = (59 * sin((45 * pi) / 180)) + 95 * sin((-26:116 * pi) / 180),
        color = "gray20",
        size = .5, 
        lineend = "round"
    ) +
    
    ## Basepaths ##
    
    annotate(
        geom = "path",
        x = c(
            home1  = 0,
            first  = 90,
            second = 90,
            third  = 0,
            home2  = 0),
        y = c(
            home1  = 0,
            first  = 0,
            second = 90,
            third  = 90,
            home2  = 0), 
        color = "gray20",
        size = .5
    ) +
    
    ## Extended foul lines ##
    
    annotate(
        geom = "segment", 
        x = c(0, 0), 
        y = c(0, 0),
        xend = c(550, 0),
        yend = c(0, 550),
        color = "black",
        size = .5
    ) +
    
    ## Foul line up to foul poles ##
    
    annotate(
        geom = "segment",
        x = 0,
        y = 0,
        xend = c(rogers.centre.df$y[2], 0),
        yend = c(0, rogers.centre.df$x[length(rogers.centre.df$x) - 1]),
        color = "black",
        size = .75
    ) +
    
    ## Bases ##
    
    annotate(
        geom = "point",
        x = c(
            first  = 90,
            second = 90,
            third  = 0),
        y = c(
            first  = 0,
            second = 90,
            third  = 90), 
        shape = 22,
        color = "gray20",
        fill = "white",
        size = 4
    ) +
    
    annotate(
        geom = "point",
        x = c(home = -2),
        y = c(home = -2), 
        shape = 23,
        color = "gray20",
        fill = "white",
        size = 4
    ) +
    
    ## Range rings ##
    
    annotate(
        geom = "path",
        x = 0 + 300 * cos((0:90 * pi) / 180),
        y = 0 + 300 * sin((0:90 * pi) / 180),
        color = "gray50",
        size = .25,
        linetype = "dashed"
    ) +
    annotate(
        geom = "path",
        x = 0 + 350 * cos((0:90 * pi) / 180),
        y = 0 + 350 * sin((0:90 * pi) / 180),
        color = "gray50",
        size = .25,
        linetype = "dashed"
    ) +
    annotate(
        geom = "path",
        x = 0 + 400 * cos((0:90 * pi) / 180),
        y = 0 + 400 * sin((0:90 * pi) / 180),
        color = "gray50",
        size = .25,
        linetype = "dashed"
    ) +
    annotate(
        geom = "path",
        x = 0 + 450 * cos((0:90 * pi) / 180),
        y = 0 + 450 * sin((0:90 * pi) / 180),
        color = "gray50",
        size = .25,
        linetype = "dashed"
    ) +
    annotate(
        geom = "path",
        x = 0 + 500 * cos((0:90 * pi) / 180),
        y = 0 + 500 * sin((0:90 * pi) / 180),
        color = "gray50",
        size = .25,
        linetype = "dashed"
    ) +
    annotate(
        geom = "path",
        x = 0 + 550 * cos((0:90 * pi) / 180),
        y = 0 + 550 * sin((0:90 * pi) / 180),
        color = "gray50",
        size = .25,
        linetype = "dashed"
    ) +
    
    ## Text ##
    
    annotate(
        geom = "text",
        label   = 300,
        x = -10 + c(300, 300) * cos((c(6, 84) * pi) / 180),
        y = -10 + c(300, 300) * sin((c(6, 84) * pi) / 180)
    ) +
    annotate(
        geom = "text",
        label   = 350,
        x = -10 + c(350, 350) * cos((c(6, 84) * pi) / 180),
        y = -10 + c(350, 350) * sin((c(6, 84) * pi) / 180)
    ) +
    annotate(
        geom = "text",
        label   = 400,
        x = -10 + c(400, 400) * cos((c(6, 84) * pi) / 180),
        y = -10 + c(400, 400) * sin((c(6, 84) * pi) / 180)
    ) +
    annotate(
        geom = "text",
        label   = 450,
        x = -10 + c(450, 450) * cos((c(6, 84) * pi) / 180),
        y = -10 + c(450, 450) * sin((c(6, 84) * pi) / 180)
    ) +
    annotate(
        geom = "text",
        label   = 500,
        x = -10 + c(500, 500) * cos((c(6, 84) * pi) / 180),
        y = -10 + c(500, 500) * sin((c(6, 84) * pi) / 180)
    ) +
    annotate(
        geom = "text",
        label   = 550,
        x = -10 + c(550, 550) * cos((c(6, 84) * pi) / 180),
        y = -10 + c(550, 550) * sin((c(6, 84) * pi) / 180)
    ) +
    annotate(
        geom = "point",
        x = c(rogers.centre.df$y[2], 0),
        y = c(0, rogers.centre.df$x[length(rogers.centre.df$x) - 1]),
        shape = 21,
        size = 2,
        color = "black",
        fill = "white"
    ) +
    
    geom_spoke(
        data = rogers_centre_data,
        aes(angle = horiz_angle,
            radius = -hit_distance_sc),
        color = "gray30",
        size = .375,
        alpha = .475) +
    
    geom_point(
        data = rogers_centre_data,
        aes(fill = launch_speed),
        shape = 21,
        color = "black",
        size = 2.75,
        alpha = 1
    ) +
    
    ## Field percentages ##
    
    annotate(
        geom = "text",
        label = rogers_centre_data %>% 
            with(., prop.table(table(horiz_angle_3)) * 100) %>%
            round(digits = 1) %>%
            paste(., "%", sep = "") %>%
            c(), 
        x = c(525, 475, 475) * cos((c(45, 75, 15) * pi) / 180),
        y = c(525, 475, 475) * sin((c(45, 75, 15) * pi) / 180),
        size = 6,
        fontface = "bold"
    ) +
    
    # scale_fill_gradientn(
    #     colors = c("gray80", "green2", "yellow", "#FF0000", "#FF4DA7"),
    #     name = "Launch\nSpeed",
    #     limits = c(90, 115),
    #     breaks = c(90, 95, 100, 105, 110, 115)
    # ) +
    
    scale_fill_viridis_c(
        name = "Launch\nSpeed",
        option = "C",
        limits = c(90, 120),
        breaks = c(90, 95, 100, 105, 110, 115, 120)
    ) +
    
    guides(
        fill = guide_colorbar(
            ticks = TRUE,
            label = TRUE,
            order = 2,
            title.position = "top",
            label.position = "left"),
        color = guide_colorbar(
            ticks = TRUE,
            label = TRUE,
            order = 1,
            title.position = "top",
            label.position = "left")
    ) +
    
    theme(
        axis.text = element_blank(),
        axis.title = element_blank(),
        axis.line = element_blank(),
        axis.ticks = element_blank(),
        panel.border = element_rect(color = "black", linetype = "solid", size = .5, fill = NA),
        strip.background = element_blank(),
        strip.text = element_text(size = 9),
        legend.key.size = unit(.75, "cm"),
        legend.key.height = unit(1.675, "cm"),
        legend.position = "left",
        legend.title = element_text(size = 15),
        legend.text = element_text(size = 10),
        panel.background = element_rect(fill = "ivory4"), 
        panel.grid = element_blank(),
        legend.title.align = 1
    )


## Saving plots ####

ggsave(
    plot = diamond.plot.field.RC,
    "plots/Rogers_Centre_Home_Runs_By_Field.png",
    width = 10,
    height = 9,
    dpi = 250,
    scale = 1
)

ggsave(
    plot = diamond.plot.field.RC,
    glue("plots/Rogers_Centre_Home_Runs_By_Field_{most_recent_day}.png"),
    width = 10,
    height = 9,
    dpi = 250,
    scale = 1
)


# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
# #
# #                             ---- THIS IS THE END! ----
# #
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #