Cartographier un département de France

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

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

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

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

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

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

# 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)