R/get_density_map.R

In bclroger/euparliament: Analysis of the EU Parliament elections of 2014 and 2019

# Function 6: plots density map

get_density_map <- function(data, fill_by, legend = "", title ="", gradmin = 0,
                            gradmax = 100, gradnum = 4, colors=c("#FFFFFF","#224d77")){
  #The arguments are:
  #data: dataframe with the polygon data and other data
  #fill_by: The column that with the desired info to be shown in the map
  #legend: string to be used as legend of the map
  #title: string to be used as title of map
  #colors: the color gradient of "fill_by" info. Default: "#FFFFFF","#224d77"
  #gradmin: lowest value shown in the gradient scale. Default: 0
  #gradmax: highest value shown in the gradient scale. Default: 100
  #grandnum: number of intervals shown in the gradient scale. Default: 4


  #Creating the breaks:
  gradbreak <- seq(from = gradmin, to = gradmax, by = ((gradmax - gradmin)/gradnum))
  #This creates a vector from the lowest value chosen to be shown in the gradient scale
  # to the highest chosen value with the desired number number of intervals. The vector
  # goes from gradmin to gradmax in a sequence of numbers increasing by equal numbers.
  # This vector is later used as an input to create the gradient scale.


  # Create world plot:
  p <- ggplot(data = data) +
    geom_sf(aes(fill = fill_by)) +
    labs(fill = legend) +
    scale_fill_gradient(low = colors[1], high = colors[2], breaks = gradbreak,
                        labels = gradbreak, limits = c(gradmin, gradmax), name = legend) +
    ggtitle(title) +
    theme(plot.title = element_text(hjust = 0.5, face = "bold")) +
    coord_sf(crs = "+init=epsg:3035", xlim = c(2500000, 6600000), ylim = c(5500000, 1200000), expand = FALSE)

  #A map is created using dataframe assigned to "data" plus the geometries of entries assigned to "fill_by"
  #shown using the color gradient assigned to "color" and using the text assigned to "legend"
  # The colors of the gradient are set by the values in argument vector "colors"
  # ggtitle creates the title of the map using the text from the argument title and the title is in bold
  # coord_sf sets an Azimuthalequal area projection for the map and the extent of the map
  # are set by xlim and ylim


  return(p)
}