Prediction grid

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)
}