Chapitre 12 Créer des cartes pour le web
12.1 Les cartes Leaflet
{leaflet} est une bibliothèque javascript de cartographie en ligne.
R permet de produire des cartes en exploitant cette bibliothèque.
La limitation de {leaflet} est qu’il ne permet de visualiser que des données en projection WGS84.
Dans ce chapitre, nous utiliserons les packages suivants
# CRAN
library(htmltools)
library(htmlwidgets)
library(leaflet)
library(sf)
library(tidyverse)
# Github
library(COGiter) # remotes::install_github("MaelTheuliere/COGiter")
library(ggiraph)
# Gitlab-forge
library(mapfactory) # remotes::install_gitlab("dreal-pdl/csd/mapfactory", host = "gitlab-forge.din.developpement-durable.gouv.fr")Les données utilisées sont issues de {COGiter}.
12.1.1 Carte choroplète
Ci-dessous un exemple avancé de carte choroplète avec {leaflet}.
- Définition des données de départ :
pop2015_dep_geo <- pop2015 %>%
# mise à jour du COG et agrégation à l'échelle départementale
cogifier(code_commune = DEPCOM, communes = FALSE, epci = FALSE, departements = TRUE, regions = FALSE, metro = FALSE) %>%
inner_join(COGiter::departements_geo, ., by = c("DEP" = "CodeZone")) %>%
mutate(densite = pop2015 / as.numeric(AREA) * 1000000) %>%
filter(substr(DEP, 1, 2) != "97")
glimpse(pop2015_dep_geo)Rows: 96
Columns: 9
$ DEP <fct> 01, 02, 03, 04, 05, 06, 07, 08, 09, 10, 11, 12, 13, 14,…
$ AREA [m^2] 5762580000 [m^2], 7361650000 [m^2], 7340110000 [m^2], 6…
$ geometry <GEOMETRY [m]> POLYGON ((904448.6 6504924,..., POLYGON ((7421…
$ TypeZone <fct> Départements, Départements, Départements, Départements,…
$ Zone <fct> Ain, Aisne, Allier, Alpes-de-Haute-Provence, Hautes-Alp…
$ pop2015 <dbl> 631877, 538659, 341613, 161799, 140916, 1082440, 324209…
$ pop2015_a_part <dbl> 17135, 13870, 10013, 4836, 5144, 15116, 9572, 7860, 540…
$ pop2015_totale <dbl> 649012, 552529, 351626, 166635, 146060, 1097556, 333781…
$ densite <dbl> 109.65175, 73.17096, 46.54058, 23.36373, 25.49205, 251.…- Transformation de la projection car
leafletne connait que le WGS 84 :
pop2015_dep_geo <- st_transform(pop2015_dep_geo, crs=("+proj=longlat +datum=WGS84 +no_defs"))- Discrétisation de la variable d’intérêt :
bornes <- quantile(pop2015_dep_geo$densite, na.rm = TRUE)- Création d’une palette de couleurs associée :
{leaflet} intègre une fonction colorBin() qui permet d’associer à un vecteur de valeurs numériques un vecteur de couleurs en fonction d’une palette et d’un vecteur de bornes.
pal <- colorBin("YlOrRd", domain = pop2015_dep_geo$densite, bins = bornes)- Création d’étiquettes popover, à afficher en surbrillance au passage de la souris sur la carte :
popover <- sprintf("<strong>%s</strong><br/>%.1f habitants par km² en 2015",
pop2015_dep_geo$Zone, pop2015_dep_geo$densite) %>%
lapply(htmltools::HTML)- Réalisation de la carte :
La fonction addPolygons() permet de créer la carte choroplète.
La fonction addProviderTiles() permet de rajouter un fond cartographique parmis les couches fournies par leaflet.
leaflet(pop2015_dep_geo) %>%
addProviderTiles(providers$CartoDB.Positron) %>%
addPolygons(data = pop2015_dep_geo,
fillColor = ~pal(densite),
weight = 2,
opacity = 1,
color = "white",
dashArray = "3",
fillOpacity = 0.7,
highlight = highlightOptions(
weight = 5,
color = "#666",
dashArray = "",
fillOpacity = 0.7,
bringToFront = TRUE),
label = popover,
labelOptions = labelOptions(
style = list("font-weight" = "normal", padding = "3px 8px"),
textsize = "15px",
direction = "auto"))12.1.2 Carte à ronds proportionnels
Pour réaliser des ronds proportionnels, il va faut partir des centroïdes de nos polygones. On cherche cette fois à représenter les populations 2015.
- Création de la table en entrée :
dep_centr <- st_centroid(pop2015_dep_geo)- Discrétisation de la variable d’intérêt :
bornes <- quantile(dep_centr$pop2015 / 1000, na.rm = TRUE)- Création d’une palette de couleurs associée :
pal <- colorBin("YlOrRd", domain = dep_centr$pop2015 / 1000, bins = bornes)- Création d’un label ad-hoc à afficher en surbrillance au passage de la souris sur la carte :
popover <- sprintf("<strong>%s</strong><br/>%.1f milliers habitants en 2015",
dep_centr$Zone, dep_centr$pop2015 / 1000) %>%
lapply(htmltools::HTML)- Création de la carte :
carte_rond_proportionnel <- leaflet(dep_centr) %>%
addProviderTiles(providers$CartoDB.Positron) %>%
addCircles(data = dep_centr,
fillColor = ~pal(pop2015 / 1000),
radius = ~2500*log(pop2015),
weight = 2,
opacity = 1,
color = "gray",
dashArray = "3",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 5,
color = "#666",
dashArray = "",
fillOpacity = 0.98,
bringToFront = TRUE),
label = popover,
labelOptions = labelOptions(
style = list("font-weight" = "normal", padding = "3px 8px"),
textsize = "15px",
direction = "auto")) %>%
addLegend("bottomright", pal = pal, values = ~pop2015,
title = "Nombre d'habitants en 2015 - INSEE - RP",
opacity = 1)
carte_rond_proportionnel12.2 Les cartes ggplot
Les cartes créées par {ggplot2}, comme tous les autres graphiques produits avec ggplot d’ailleurs, peuvent être rendues dynamiques grâce à {ggiraph} développé par David Gohel.
Ce package permet l’ajout d’info-bulles, des effets de survol et d’ajouter des actions JavaScript aux graphiques. L’interactivité est s’applique aux géométries ggplot, aux légendes et aux éléments de thème, grâce aux paramètres esthétiques suivants :
tooltip: info-bulles à afficher lorsque la souris se trouve sur les éléments ;
onclick: fonction JavaScript à exécuter lorsque les éléments sont cliqués ;
data_id: id à associer aux éléments (utilisé pour les actions de survol et de clic).
Les choses que vous devez savoir pour créer un graphique interactif :
- Au lieu d’utiliser
geom_sf(), utilisezgeom_sf_interactive(), au lieu d’utilisergeom_sf_label(), utilisezgeom_sf_label_interactive()… Fournissez au moins une des esthétiquestooltip,data_idandonclickpour créer des éléments interactifs.
- Appelez la fonction girafe avec l’objet ggplot afin que le graphique soit traduit en un graphique web interactif.
On reprend notre carte ggplot normande en ronds proportionnels :
carte_ronds_prop_dyn <- ggplot(data = pop2015_epci_normands_geo) +
geom_sf_interactive(aes(tooltip = Zone)) +
geom_sf_interactive(data = st_point_on_surface(pop2015_epci_normands_geo),
mapping = aes(size = pop2015_totale, color = poids_pop_a_part,
tooltip = paste0(Zone, "\n", format_fr_nb(pop2015_totale, dec = 0, big_mark = " "), " hab.\ndont ",
format_fr_nb(poids_pop_a_part, dec = 1, big_mark = " "), " % comptés à part"))) +
labs(
title = "Population par EPCI en Normandie",
subtitle = "En 2015", x = "", y = "",
size = "Nombre d'habitants",
color = "Poids population\ncomptée à part (%)",
caption = "Source : INSEE - RP"
) +
theme_gouv_map(plot_title_size = 20, subtitle_size = 12) +
scale_color_gouv_continuous(palette = "pal_gouv_i", reverse = TRUE) +
scale_size_continuous(labels = function(x) format(x, big.mark = " ", scientific = FALSE)) +
annotation_north_arrow(location = "bl", height = unit(0.7, "cm"), width = unit(0.7, "cm"))
girafe(ggobj = carte_ronds_prop_dyn)12.3 Les cartes {mapfactory}
{mapfactory} s’appuie sur {ggiraph} pour rendre ses cartes dynamiques.
Les arguments de la fonction creer_carte() en rapport avec l’interactivité sont
interactive
TRUE si on veut une carte html interactive.
popover
Chaîné de texte désignant le contenu de l’info-bulles. Par défaut, L’étiquette qui s’affiche au survol contient “nom du territoire : valeur de l’indicateur et suffixe”, obtenu a partir de “.data\(nom_ter : .data\)valeur”.
indicateur_popover
Un nom d’indicateur présent dans data (sans quote) a mobiliser pour le popover. Utiliser .data$valeur_popover au niveau du parametre popover pour afficher les valeurs de cet indicateur, par exemple .data$nom_ter : .data$valeur_popover unites.
creer_carte(data = cogifier(pop2015), code_region = "52", indicateur = pop2015, espace = " ",
interactive = TRUE, type_viz = "proportionnel", maillage = "Départements")