make_athens_plot <- function(spatraster, title) {
ggplot() +
geom_spatraster(data = spatraster) +
geom_sf(data = map, fill = NA) +
coord_sf(datum = "EPSG:2100") +
scale_fill_viridis_c(
option = "magma",
direction = -1,
na.value = "transparent"
) +
ggtitle(title)
}