Voici un script R possible pour le TP 2 de ce document :
Activité 1 : Manipulation des cartes1
1- Télécharger le shapefile du Cameroun2
library(rgdal) # pour lire les shapefiles library(sf) # pour manipuler les objets
spatiaux library(sp) # pour afficher les cartes library(tmap) # pour faire des cartes
thématiques
Changer le répertoire de travail vers le
dossier contenant le shapefile
setwd(“C:/Users/…/gadm36_CMR_shp”)
Lire le shapefile du Cameroun sans
subdivisions
cameroun <- readOGR(dsn = “.”, layer = “gadm36_CMR_0”)
2- Représenter la carte du Cameroun sans
les subdivisions3
plot(cameroun, main = “Carte du Cameroun”)
Lire le shapefile du Cameroun avec les
subdivisions de niveau 1 (régions)
cameroun_reg <- readOGR(dsn = “.”, layer = “gadm36_CMR_1”)
Représenter la carte du Cameroun avec les
régions
plot(cameroun_reg, main = “Carte du Cameroun par région”)
Lire le shapefile du Cameroun avec les
subdivisions de niveau 2 (départements)
cameroun_dep <- readOGR(dsn = “.”, layer = “gadm36_CMR_2”)
Représenter la carte du Cameroun avec les
départements
plot(cameroun_dep, main = “Carte du Cameroun par département”)
Lire le shapefile du Cameroun avec les
subdivisions de niveau 3 (arrondissements)
cameroun_arr <- readOGR(dsn = “.”, layer = “gadm36_CMR_3”)
Représenter la carte du Cameroun avec les
arrondissements
plot(cameroun_arr, main = “Carte du Cameroun par arrondissement”)
3- Mettre en exergue une région de votre
choix (Centre)4
centre <- subset(cameroun_reg, NAME_1 == “Centre”) plot(cameroun_reg, col =
“grey”) plot(centre, col = “blue”, add = TRUE, main = “Carte du Centre”)
Mettre en exergue les départements de
votre choix (Kadey, Lom et Djerem, Mefou
et Afamba)4
kadey <- subset(cameroun_dep, NAME_2 == “Kadey”) lom_djerem <-
subset(cameroun_dep, NAME_2 == “Lom-et-Djerem”) mefou_afamba <-
subset(cameroun_dep, NAME_2 == “Méfou-et-Afamba”) plot(cameroun_dep, col
= “grey”) plot(kadey, col = “red”, add = TRUE) plot(lom_djerem, col = “green”,
add = TRUE) plot(mefou_afamba, col = “yellow”, add = TRUE, main = “Carte des
départements choisis”)
4- Représenter le département de la Kadey
individuellement5
plot(kadey, main = “Carte du Kadey”)
Représenter les départements qui lui sont
limitrophes5
voisins <- gTouches(kadey, cameroun_dep, byid = TRUE, returnDense = TRUE)
voisins <- [Link](voisins) voisins <- subset(voisins, V1 == TRUE) voisins
<- cameroun_dep[rownames(voisins), ] plot(cameroun_dep, col = “grey”)
plot(kadey, col = “red”, add = TRUE) plot(voisins, col = “green”, add = TRUE,
main = “Carte du Kadey et ses voisins”)
5- Matérialiser le centroïde du centre6
centre_centroide <- gCentroid(centre, byid = TRUE) plot(cameroun_reg, col =
“grey”) plot(centre, col = “blue”, add = TRUE) plot(centre_centroide, pch = 19, col
= “red”, add = TRUE, main = “Carte du centre avec son centroïde”)
Matérialiser le centroïde de la Kadey6
kadey_centroide <- gCentroid(kadey, byid = TRUE) plot(cameroun_dep, col =
“grey”) plot(kadey, col = “red”, add = TRUE) plot(kadey_centroide, pch = 19, col
= “blue”, add = TRUE, main = “Carte du Kadey avec son centroïde”)
6- Matérialiser une zone tampon de 300km
autour du centroïde de la Kadey7
kadey_buffer <- gBuffer(kadey_centroide, width = 300000, byid = TRUE)
plot(cameroun_dep, col = “grey”) plot(kadey, col = “red”, add = TRUE)
plot(kadey_buffer, col = “blue”, add = TRUE, main = “Carte du Kadey avec une
zone tampon de 300km”)7
7- Récupérer sur une carte isolée la zone
tampon8
plot(kadey_buffer, main = “Carte de la zone tampon de 300km autour du Kadey”)7
8- Refaire les questions 1-7 en considérant
la carte de votre pays et votre
arrondissement de naissance9
A vous de jouer !
