R est un langage et un environnement open source permettant de réaliser une variété de traitements statistiques et de représentations graphiques.

Son utilisation en géographie permet de simplifier les chaînes de traitements et d’inscrire pleinement les analyses dans une démarche de recherche reproductible.

Lors de cette séance nous aborderons les fonctionnalités spatiale du logiciel et plus particulièrement le package de cartographie thématique cartography.

Quelques manipulations de données

# Connaitre le répertoire de travail
getwd()
[1] "/data/prz/pi"
# changer le répertoire de travail setwd(dir = '/mon_chemin/')

# import d'un jeu de données, le format texte délimité (type csv) est
# fortement conseillé.
datanuts <- read.csv("data/datanuts.csv")

# Aperçu des données
View(datanuts)

# Calcul de quelques statistiques sur les variables
summary(datanuts)
       id         birth_2008       death_2008      gdppps1999    
 AT111  :   1   Min.   :   122   Min.   :  139   Min.   :   296  
 AT112  :   1   1st Qu.:  1217   1st Qu.: 1414   1st Qu.:  1844  
 AT113  :   1   Median :  2476   Median : 2550   Median :  3476  
 AT121  :   1   Mean   :  4722   Mean   : 3643   Mean   :  6440  
 AT122  :   1   3rd Qu.:  5448   3rd Qu.: 4638   3rd Qu.:  6666  
 AT123  :   1   Max.   :220464   Max.   :43827   Max.   :124610  
 (Other):1442                    NA's   :81                      
   gdppps2008        pop1999           pop2008        
 Min.   :   296   Min.   :  14873   Min.   :   15471  
 1st Qu.:  2643   1st Qu.: 137979   1st Qu.:  135413  
 Median :  5185   Median : 255661   Median :  255380  
 Mean   :  9516   Mean   : 388349   Mean   :  402646  
 3rd Qu.:  9764   3rd Qu.: 480745   3rd Qu.:  492828  
 Max.   :229459   Max.   :9692579   Max.   :12573836  
                                                      
# Affichage d'un graphique simple
plot(datanuts[, c("gdppps2008", "pop2008")])
title("GDP  - POP")

Import et manipulation de données spatiales

# chargement du package rgdal
library("rgdal")

# import d'un shapefile dans le dossier data
nuts3 <- readOGR(dsn = "data", layer = "nuts", verbose = TRUE)
OGR data source with driver: ESRI Shapefile 
Source: "data", layer: "nuts"
with 1448 features
It has 1 fields
# affichage du fond de carte
plot(nuts3)

# affichage plus personnalisé
plot(nuts3, col = "#DAE3E6", border = "#8A0641", lwd = 0.5)