# Map census data
#'
#' @param census_data, dataframe of census data for specified variable of interest
#' @param conus_sf, sf of conus states outline
#' @param leg_title, character string for legend title
#' @param outfile_path, outfile path for pngs
#' @param viz_config_df `data.frame` width, height, counties outline color, conus outline color, background color, font nam, and font size
#' @param viz_config_pal `data.frame` assign colors for postively and negatively correlated dimensions for census maps
#' @param width, set figure width dimension
#' @param height, set figure height dimension
#' @param percent_leg if else statement where if TRUE, apply 0-100 legend, otherwise retain 0 - max of variable name
#' @param dollar_leg if else statement where if TRUE, apply $ to the estimate column to display cost in USD on map
#' @param font_size, set font size
#' @param barheight, set colorbar bar height
#' @param barwidth, set colorbar bar width
#' @param lim_vals, set legend limits
#' @param break_vals, set legend breaks
#' @param low_ramp_col, assign color to low end of color ramp for legend
plot_census_map <- function(census_data, conus_sf, leg_title, outfile_path, var,
percent_leg, viz_config_df, viz_config_pal, width, height,
font_size, barwidth, barheight, lim_vals, break_vals,
dollar_leg, low_ramp_col){
font_legend <- viz_config_df$load_font
font_add_google(font_legend)
showtext_opts(dpi = 300, regular.wt = 200, bold.wt = 700)
showtext_auto(enable = TRUE)
census_map <- census_data |>
ggplot(aes(fill = .data[[var]])) +
geom_sf(color = viz_config_df$counties_outline_col,
linewidth = 0.05) +
geom_sf(data = conus_sf,
fill = NA,
color = viz_config_df$conus_outline_col,
linewidth = 0.5,
linetype = "solid") +
theme_void() +
theme(text = element_text(family = viz_config_df$font_legend, size = font_size),
legend.margin = margin(t = 5, b = 5), legend.position = 'bottom', legend.title.align = 0.5) +
guides(fill = guide_colorbar(
title.position = "top",
title.theme = element_text(face = 'bold', family = viz_config_df$font_legend, size = font_size),
direction = "horizontal",
position = "bottom",
barwidth = barwidth,
barheight = barheight
))
if (percent_leg == FALSE) {
census_map <- census_map +
scale_fill_gradientn(
colors = colorRampPalette(c(low_ramp_col, viz_config_pal))(100),
name = leg_title,
limits = lim_vals,
labels = if (dollar_leg) scales::label_dollar() else scales::comma,
na.value = "#F5F5F5"
)
} else {
census_map <- census_map +
scale_fill_gradientn(
colors = colorRampPalette(c("#F8F9FF", viz_config_pal))(100),
name = leg_title,
limits = lim_vals, # c(0, 100),
breaks = break_vals, # c(0, 25, 50, 75, 100),
labels = if (dollar_leg) scales::label_dollar() else function(x) paste0(x, "%"),
na.value="#F5F5F5"
)
}
background_color = "white"
plot_margin = 0.025
# cowplot to get map sizes larger
canvas <- grid::rectGrob(
x = 0, y = 0,
width = width, height = height,
gp = grid::gpar(fill = background_color, alpha = 1, col = background_color
)
)
census_legend <- get_plot_component(census_map, 'guide-box', return_all = TRUE)
# compose final plot
final_map <- ggdraw(ylim = c(0,1),
xlim = c(0,1)) +
# White background
draw_grob(canvas,
x = 0, y = 1,
height = height, width = width,
hjust = 0, vjust = 1) +
# Add main plot
draw_plot(census_map + theme(legend.position="none"),
x = -0.01,
y = 0.08,
height = 0.98,
width = (1-plot_margin)*1.03) +
# Add legend
draw_plot(census_legend[[3]],
x = 0.48,
y = 0.02,
height = 0.09 ,
width = 0.1 - plot_margin)
ggsave(outfile_path, final_map, width = width, height = height, dpi = viz_config_df$dpi, bg = viz_config_df$bg_col, units = "in")
return(outfile_path)
}
# Map raster data
#'
#' @param raster_data .tif file of raster data for specified variable of interest
#' # notes: update colorpalette according to Elmera's request
plot_raster <- function(in_raster, conus_sf, conus_inner, outfile_path,
viz_config_df, viz_config_pal, width, height,
font_size, barwidth, barheight, counties_sf,
low_ramp_col){
# in_raster = p2_pop_density_processed, conus_sf = p2_conus_sf_proj, conus_inner = p2_conus_inner, viz_config_pal = p0_viz_config_pal$living_conditions
# width = p0_viz_config_df$width_desktop, height = p0_viz_config_df$height_desktop, font_size = p0_viz_config_df$font_size_desktop, barwidth = 20, barheight = 1,
#outfile_path = "3_visualize/out/pop_density_rast_2020_en.png", low_ramp_col = "#eef0ff",
raster_data <- rast(in_raster)
# plot
font_legend <- viz_config_df$load_font
font_add_google(font_legend)
showtext_opts(dpi = 300, regular.wt = 200, bold.wt = 700)
showtext_auto(enable = TRUE)
# doesn't work with cowplot::get_plot_component()
legend_label <- expression(paste(bold("Population per 1 "), bold(km^2), bold(", 2020")))
(pop_map <- ggplot() +
#geom_spatraster(data = pop_usa_proj)+
geom_sf(data = conus_sf, fill = low_ramp_col, color = low_ramp_col, linewidth = 0.5) +
geom_spatraster(data = raster_data)+ #x = x, y = y,
geom_sf(data = conus_inner, fill = NA, color = viz_config_df$counties_outline_col, linewidth = 0.5) +
geom_sf(data = counties_sf, fill = NA, color = viz_config_df$counties_outline_col, linewidth = 0.1) +
scale_fill_gradientn(
colors = colorRampPalette(c(low_ramp_col, viz_config_pal))(100),
name = legend_label,
labels = c(1, 10, 100, '1k', '10k', '100k'),
na.value=NA
)+
labs(fill = legend_label) +
theme(plot.background = element_rect(fill = "white", color = "white"),
panel.background = element_rect(fill = "white", color = "white"),
axis.title = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
panel.grid = element_blank(),
legend.background = element_blank(),
legend.direction = "horizontal",
legend.position = "bottom",
legend.title.align = 0.5) +
guides(fill = guide_colorbar(
title.position = "top",
title.theme = element_text(face = 'bold', family = font_legend, size = font_size),
direction = "horizontal",
position = "bottom",
barwidth = barwidth,
barheight = barheight
))
)
background_color = "white"
plot_margin = 0.025
# cowplot to get map sizes larger
canvas <- grid::rectGrob(
x = 0, y = 0,
width = 6, height = 6,
gp = grid::gpar(fill = background_color, alpha = 1, col = background_color
)
)
pop_map_legend <- get_plot_component(pop_map, 'guide-box', return_all = TRUE)
# compose final plot
final_map <- ggdraw(ylim = c(0,1),
xlim = c(0,1)) +
# White background
draw_grob(canvas,
x = 0, y = 1,
height = 6, width = 6,
hjust = 0, vjust = 1) +
# Add main plot
draw_plot(pop_map + theme(legend.position="none"),
x = -0.01,
y = 0.08,
height = 0.98,
width = (1-plot_margin)*1.03) +
# Add legend
draw_plot(pop_map_legend[[3]],
x = 0.48,
y = 0.02,
height = 0.09 ,
width = 0.1 - plot_margin)
ggsave(outfile_path, final_map, width = width, height = height, dpi = viz_config_df$dpi, bg = viz_config_df$bg_col, units = "in")
}
plot_imp_surf_raster <- function(in_raster, conus_sf, conus_inner, conus_col, outfile_path,
viz_config_df, viz_config_pal, width, height,
font_size, barwidth, barheight, counties_sf,
low_ramp_col){
# in_raster = p2_imp_surf_processed, conus_sf = p2_conus_sf_proj, conus_col = "gray90", viz_config_pal = p0_viz_config_pal$living_conditions
# width = p0_viz_config_df$width_desktop, height = p0_viz_config_df$height_desktop, font_size = p0_viz_config_df$font_size_desktop, barwidth = 20, barheight = 1,
#outfile_path = "3_visualize/out/imp_surface_rast_2020_en.png", low_ramp_col = "white"
raster_data <- rast(in_raster)
# plot
font_legend <- viz_config_df$load_font
font_add_google(font_legend)
showtext_opts(dpi = 300, regular.wt = 200, bold.wt = 700)
showtext_auto(enable = TRUE)
(pop_map <- ggplot() +
geom_sf(data = conus_sf, fill = NA, color = conus_col, linewidth = 0.9) +
geom_spatraster(data = raster_data)+
geom_sf(data = conus_inner, fill = NA, color = conus_col, linewidth = 0.5) +
geom_sf(data = counties_sf, fill = NA, color = conus_col, linewidth = 0.1) +
scale_fill_gradientn(
colors = colorRampPalette(c(low_ramp_col, viz_config_pal))(100),
name = "Percent impervious surface, 2022",
limits = c(0, 100),
breaks = c(0, 25, 50, 75, 100),
labels = c("0%", "25%", "50%", "75%", "100%"),
na.value=NA
)+
labs(fill = "% Impervious Surface, 2022") +
theme(plot.background = element_rect(fill = "white", color = "white"),
panel.background = element_rect(fill = "white", color = "white"),
axis.title = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
panel.grid = element_blank(),
legend.background = element_blank(),
legend.direction = "horizontal",
legend.position = "bottom",
legend.title.align = 0.5) +
guides(fill = guide_colorbar(
title.position = "top",
title.theme = element_text(face = 'bold', family = font_legend, size = font_size),
direction = "horizontal",
position = "bottom",
barwidth = barwidth,
barheight = barheight
))
)
background_color = "white"
plot_margin = 0.025
# cowplot to get map sizes larger
canvas <- grid::rectGrob(
x = 0, y = 0,
width = 6, height = 6,
gp = grid::gpar(fill = background_color, alpha = 1, col = background_color
)
)
pop_map_legend <- get_plot_component(pop_map, 'guide-box', return_all = TRUE)
# compose final plot
final_map <- ggdraw(ylim = c(0,1),
xlim = c(0,1)) +
# White background
draw_grob(canvas,
x = 0, y = 1,
height = 6, width = 6,
hjust = 0, vjust = 1) +
# Add main plot
draw_plot(pop_map + theme(legend.position="none"),
x = -0.01,
y = 0.08,
height = 0.98,
width = (1-plot_margin)*1.03) +
# Add legend
draw_plot(pop_map_legend[[3]],
x = 0.48,
y = 0.02,
height = 0.09 ,
width = 0.1 - plot_margin)
ggsave(outfile_path, final_map, width = width, height = height, dpi = viz_config_df$dpi, bg = viz_config_df$bg_col, units = "in")
}