2.1- Ouvrir un fond de carte vectoriel et colorer sélectivement des aires (régions, communes...)

plot(dordogne)

# Ajouter une coloration rouge pour la 516ème commune

plot (dordogne[516,],col="red",add=T) # add = T pour superposer à la carte précédente

# Ajouter une coloration verte pour la commune ayant le code_insee 24035

plot (dordogne[which(dordogne$code_insee==24035),],col="green",add=T)

plot(dordogne,col="#CCFFCC")

# Colorer la commune 01 de la liste en bleu

plot (dordogne[1,],col="blue",add=T)

# Récupérer les coordonnées moyennes de la commune colorée en bleu

x <- mean(coordinates(dordogne@polygons[[1]]@Polygons[[1]])[,1])

y <- mean(coordinates(dordogne@polygons[[1]]@Polygons[[1]])[,2])

x ; y

# Tracer le cercle autour du centre de gravité de la commune

points(x,y,col="red",cex=5,lwd=3)

Je télécharge une base de données qui contient plein de noms de villes avec des caractères spéciaux (é,ê en français - ou par exemple ici, des noms grecs). Cela va me poser tout de suite de nombreux problèmes (fichier lourd, caractères spéciaux qui s'affichent mal, données à trier...).

library("raster")

adm_gr <- getData('GADM', country='GRC', level=1)

J'affiche la carte :

plot(adm_gr,xlim=c(22.4,22.5),ylim=c(36.38,36.55),lwd=2,bg="#B6F3EC",col="#F2F3B6")

Je charge les noms de lieu :

install.packages(c("maps","maptools","rgdal"))

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

# le fichier à télécharger des noms de ville se retrouve ici.

places <- readOGR(dsn = ".", layer = "places",p4s=CRS("+proj=longlat"))

# Décrire le fichier

names(places) # on voit que les noms des villes apparaissent dans la sous-rubrique name

On peut donc les ajouter avec text()

text(places,places$name)

Et là, les problèmes commencent ! C'est long à afficher et c'est moche !

Problème 1, les noms sont abherrants car caractères imbuvables.

Encoding(levels(places$name)) <- "latin1"

Encoding(levels(places$name)) <- "UTF-8"

plot(adm_gr,xlim=c(22.4,22.5),ylim=c(36.38,36.55),lwd=2,bg="#B6F3EC",col="#F2F3B6")

text(places,places$name)

Maintenant, ça marche mais c'est toujours lent à afficher (car il y a tous les noms du pays : il faut extraire ce qui nous intéresse)

Problème 2, alléger le fichier nom en n'affichant que ceux visibles sur la carte

Je dois extraire les coordonnées du fichier pour extraire les villes qui m'intérèssent, problème ici ce sont des id osm, autrement dit des identifiants Open Street Map. Je dois donc convertir ces id osm avec la fonction coordinates().

lon <- coordinates(places)[,1]

lat <- coordinates(places)[,2]

id <- which(lon>=22.4&lon<22.6&lat>=36.38&lat<=36.55)

noms_sur_carte <- places[id,]

plot(adm_gr,xlim=c(22.4,22.5),ylim=c(36.38,36.55),lwd=2,bg="#B6F3EC",col="#F2F3B6")

text(noms_sur_carte,noms_sur_carte$name)

Problème 3, n'afficher que les principaux patelins

names(noms_sur_cartes)

noms_sur_carte$type

plot(adm_gr,xlim=c(22.4,22.5),ylim=c(36.38,36.55),lwd=2,bg="#B6F3EC",col="#F2F3B6")

text(noms_sur_carte[which(noms_sur_carte$type=="village"),],noms_sur_carte$name[noms_sur_carte$type=="village"],cex=1.1,las=1)

# Et on peut même ajouter les rivières

riviere <- readOGR(dsn = ".", layer = "waterways", p4s=CRS("+proj=longlat")) # cliquer pour télécharger

plot(rivieres,add=T,col="#177CED") # et on ajoute le tracé des rivières