4.1 Point sur une carte

Il est plutôt simple de mettre en évidence un point sur une carte lorsqu’on connait les coordonnées GPS. Prenons un exemple avec le pavillon Pouliot de l’Université Laval.

Tout d’abord, on garde les coordonnées du pavillon dans l’objet pouliot.

pouliot <- c(latitude = 46.779495, longitude = -71.275398)

Ensuite, on stocke la carte dans l’objet carte_pouliot. À l’aide de la fonction leaflet, on modifie la largeur et la hauteur par défaut.

carte_pouliot <- leaflet(width = 910, height = 240)

Avec la fonction addTiles, on ajoute les tuiles de fond provenant de OpenStreetMap.

carte_pouliot <- addTiles(carte_pouliot)    

La fonction addMarkers permet d’ajouter du texte en glissant le curseur sur le point.

carte_pouliot <- addMarkers(
  map = carte_pouliot,
  lng = pouliot["longitude"],
  lat = pouliot["latitude"],
  label = "Université Laval, Pavillon Pouliot"  # Texte en glissant le curseur
)

Voici une autre façon d’ajouter du texte avec la fonction addPopups. Il s’agit d’un popup qui disparaît lorsqu’on clique quelque part dans la carte.

carte_pouliot <- addPopups(
  map = carte_pouliot,
  lng = pouliot["longitude"],
  lat = pouliot["latitude"],
  popup = "Université Laval, Pavillon Pouliot",  # Texte qui apparaît au début
  options = popupOptions(closeButton = FALSE)
)

On ajoute aussi un bouton permettant de mesurer des distances et des superficies. La fonction addMeasure ajoute un bouton avec une équerre dessus. On choisit le mètre et le mètre carré pour les unités de mesures.

carte_pouliot <- addMeasure(
  map = carte_pouliot,
  position = "bottomleft",
  primaryLengthUnit = "meters",
  primaryAreaUnit = "sqmeters",
  activeColor = "#FF0000",
  completedColor = "#7D4479"
)

Finalement, on affiche la carte.

carte_pouliot

4.2 Grouper des points par couleur

Il est parfois pratique de montrer plusieurs points sur une même carte. L’exemple qui suit a été produit avec le data frame quakes du package datasets. Ces données représentent la position géographique, la magnitude et la profondeur de plusieurs séismes près des îles Fidji. Une variable contient aussi le nombre de stations sismiques qui ont rapporté chaque tremblement de terre. La colonne region a été ajoutée selon la longitude.

Dans cet exemple, on peut obtenir d’autres informations sur les points. En glissant le curseur, on obtient la magnitude du séisme et en cliquant on obtient les coordonnées GPS.

Tout d’abord, on appelle le package datasets et on ajoute la variable.

library(datasets)

quakes$region <- factor(quakes$long >= 175, labels = c("Ouest", "Est"))

On ajoute ensuite la couleur en fonction de la région avec la fonction colorFactor.

quake_couleur <- colorFactor(c("blue", "red"), domain = c("Est", "Ouest"))

On stocke ensuite la carte dans l’objet carte_quake.

carte_quake <- leaflet()

On utilise la fonction addTiles pour le fond de carte.

carte_quake <- addTiles(carte_quake)

Avec la fonction addCircleMarkers, on ajoute des marqueurs en forme de cercle. On ajoute également un popup qui contient la magnitude en glissant le curseur et les coordonnées en cliquant.

carte_quake <- addCircleMarkers(
  map = carte_quake,
  lng = quakes$long,
  lat = quakes$lat,
  popup = paste("Latitude:", quakes$lat, ", longitude:", quakes$long),
  label = paste("Magnitude", quakes$mag),
  color = quake_couleur(quakes$region),
  radius = quakes$depth/100  # Rayon des points. On divise par 100, sinon les points sont trop gros.
)

La fonction addProviderTiles permet de changer le fond de la carte.

carte_quake <- addProviderTiles(
  map = carte_quake,
  provider = providers$Stamen.Watercolor
)

On ajoute ensuite la légende des couleurs avec la fonction addLegend.

carte_quake <- addLegend(
  map = carte_quake,
  position = "topright",
  pal = quake_couleur,
  values = quakes$region,
  title = "Région",
  opacity = 1
)

Finalement, on affiche le résultat.

carte_quake

Il existe plusieurs fonds de carte. Pour avoir la liste complète disponible, il suffit d’exécuter la commande suivante:

names(providers)

4.3 Groupes avec nombres variables

Lorsqu’il y a plusieurs points à montrer à la fois, il est possible de faire une carte avec des groupes qui varient selon le niveau de zoom.

On stocke la carte dans l’objet carte_groupe.

carte_groupe <- leaflet()

On ajoute ensuite le fond de carte OpenStreetMap.

carte_groupe <- addTiles(carte_groupe)

Avec la fonction addMarkers, on ajoute des points qui sont groupés selon le niveau de zoom. On ajoute également le nombre de stations ayant détecté le séisme lorsqu’on clique sur un point.

carte_groupe <- addMarkers(
  map = carte_groupe,
  lng = quakes$long,
  lat = quakes$lat,
  popup = paste(quakes$stations, "stations ont détecté le séisme"),
  clusterOptions = markerClusterOptions()
)

On change le fond de la carte avec la fonction addProviderTiles.

carte_groupe <- addProviderTiles(
  map = carte_groupe, 
  provider = providers$Stamen.TonerLite
)

On ajoute un bouton qui permet de revenir au zoom initial. Le niveau utilisé est 4 et le bouton est en forme de planète.

carte_groupe <- addEasyButton(
  map = carte_groupe,
  button = easyButton(
    icon = "fa-globe",
    title = "Zoom initial",
    onClick = JS("function(btn, map){ map.setZoom(4); }")
  )
)

On affiche la carte.

carte_groupe