21  Import de fonds de carte

21.0.1 Fonds de cartes raster

Le package maptiles (Giraud, 2025b) permet de télécharger et d’afficher des fonds de cartes raster. La fonction get_tiles() permet de télécharger des fonds de cartes OSM au format SpatRaster du package terra.
Dans cette exemple nous utilisons le package mapsf pour afficher la carte.
Les rendus sont meilleurs si les données en entrée utilisent le même système de coordonnées que les tuiles (EPSG:3857).

library(sf)
library(maptiles)
library(mapsf)

com_raw <- st_read("data/lot.gpkg", layer = "communes", quiet = TRUE)
com <- st_transform(com_raw, 3857)

# Récupération d'un fond de carte OSM
osm_tiles <- get_tiles(x = com, zoom = 9, crop = TRUE)

mf_theme(mar = c(0,0,0,0))
mf_raster(osm_tiles)
mf_map(com, border = "grey20", col = NA, lwd = .7, add = TRUE)
mf_credits(get_credit("OpenStreetMap"), bg = "#ffffff80")

De nombreux styles de tuiles sont disponibles avec le package. En voici quelque uns:

Certains styles ne contiennent que des labels et peuvent être utilisés en complément de données vectorielles:

# Récupération d'un fond de carte OSM des labels
mf_theme(mar = c(0,0,0,0))
osm_labels <- get_tiles(x = com, provider = "CartoDB.PositronOnlyLabels")
mf_map(com, col = 'ivory', border = 'ivory3')
mf_raster(osm_labels, add = T)
mf_credits(get_credit("CartoDB.PositronOnlyLabels"), bg = "#ffffff80")

Pour certains styles, ceux fournis par Stadia ou Thunderforest par exemple, vous aurez besoin d’une clef d’API. Vous devez vous inscrire sur le site Web de ces fournisseurs pour obtenir une clef.

21.0.2 Fonds de cartes vectoriels

La base de données OpenStreetMap propose de nombreuses couches géographiques très détaillées à l’échelle mondiale. Pour obtenir des informations synthétiques pouvant être utilisées à des fins cartographiques, les couches doivent être sélectionnées, simplifiées, fusionnées ou modifiées. Le package maposm (Giraud, 2025a) fournit un ensemble de couches composites construites à partir d’extractions effectuées avec le package osmdata.

Les couches proposées sont les suivantes :

  • les zones urbaines
  • les bâtiments
  • les espaces verts
  • les routes principales
  • les outes secondaires
  • les voies ferrées
  • les étendues d’eau

Le package n’est pas sur le CRAN. Pour l’installer :

install.packages("maposm", repos = "https://riatelab.r-universe.dev")

om_get() est utilisé pour télécharger, fusionner et simplifier les couches géographiques. om_map() peut être utilisé pour cartographier les couches.

library(maposm)
coords <- com_raw[com_raw$NOM_COM == "Cahors", ] |>
  st_geometry() |>
  st_centroid() |>
  st_transform("EPSG:4326") |>
  st_coordinates() |>
  as.vector()

cahors <- om_get(x = coords, r = 2000)
#> Getting urban areas: 1.548 sec elapsed
#> Getting buildings: 6.939 sec elapsed
#> Getting green areas: 1.141 sec elapsed
#> Getting roads: 0.68 sec elapsed
#> Getting streets: 1.867 sec elapsed
#> Getting railways: 0.475 sec elapsed
#> Getting water bodies: 2.208 sec elapsed
om_map(x = cahors, title = "Cahors", theme = "light")

NoteSite web de maposm

Get Map Layers from OpenStreetMap