add year to population density map

3_visualize/src/plot_utils.R

@@ -129,7 +129,7 @@ plot_raster <- function(in_raster, conus_sf, conus_inner, outfile_path,
   showtext_auto(enable = TRUE)
   
   # doesn't work with cowplot::get_plot_component()
-  legend_label <- expression(bold("Population per 1 km"^2)) #bquote('Population per 1'~km^2)
+  legend_label <- expression(paste(bold("Population per 1 "), bold(km^2), bold(", 2020")))
   
   (pop_map <- ggplot() +
       #geom_spatraster(data = pop_usa_proj)+