Cartographier un département de France

Voici comment réaliser simplement deux cartes sur une même fenêtre pour localiser un site.

Exemple de cartes de France et de Dordogne réalisées avec le logiciel R project

1. Installer les fonctions de cartographies nécessaires à la réalisation de cartes avec R

install.packages("maptools");install.packages("rgdal");install.packages("maps");

library("maptools");library(rgdal);library(maps)

install.packages("raster");library(raster) # pour scalebar

2. Charger la carte vectorielle de France avec contours des départements - cette carte est disponible ici (source du fichier vectoriel)

# Level 2 permet d'aller au niveau des départements

adm_fr <- getData('GADM', country='FRA', level=2)

3. Déclarer le département que vous désirez colorer et les coordonnées GPS (longitude, latitude) du point à afficher

cat("Liste des départements : \n")

adm_fr$NAME_2  # Afficher les noms de départements.

# Département concernée

departement <- "Dordogne"  # A MODIFIER EN CONSERVANT L'ORTHOGRAPHE INDIQUEE

# Ville concernee

longitude <-0.719049 # A MODIFIER

latitude <- 45.187834 # A MODIFIER

4. Tracer la carte

# Afficher les cartes dans une fenêtre

layout(matrix(1:2,1,2)) # pour afficher 2 cartes

# AFFICHER LA CARTE DE FRANCE AVEC LE DÉPARTEMENT COLORE

indice = which(adm_fr$NAME_2==departement)

plot(adm_fr,col=c(rep("white",indice-1),"red",rep("white",length(adm_fr$NAME_2)-indice)),

main=departement, cex = 0.5)

# BARRE D'ECHELLE

# d = distance en km, xy = c(long, lat) de position de la barre

cat("Cliquer là où vous voulez la barre d'échelle");l = locator(n=1)

scalebar(d = 200, xy = c(l$x,l$y), type = "bar", below = "km",lwd = 4, divs = 2, col = "black", cex = 0.75, lonlat = T)

# AFFICHER LA CARTE DE DORDOGNE EN IDENTIFIANT PRÉCISÉMENT UN SITE

plot(adm_fr[indice,],col="orange",main="")

points(longitude ,latitude ,pch=16,col="red",cex=3) # cex : taille du point

# BARRE D'ECHELLE

cat("Cliquer là où vous voulez la barre d'échelle");l = locator(n=1)

scalebar(d = 20, xy = c(l$x,l$y), type = "bar", below = "km",lwd = 4, divs = 2, col = "black", cex = 0.75, lonlat = T)

Mise à jour de la page septembre 2018