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