[R] Faire une carte par points avec R

Voici ci-dessous le code R pour réaliser une carte en points à partir d’un fond de carte et d’une donnée quantitative absolue.

Screenshot from 2014-03-04 09:58:58

# VARIABLES
mygeom<-“nuts3_2006.shp”
mydata<-“data_nuts3.csv”
myvar<-“pop_t_2008”
nb<-100 # Un point vaut nb units
library(maptools)
setwd(“/home/nlambert/Documents/R/carto”)
list.files()fdc <- readShapeSpatial(mygeom)
donnees<-read.csv(mydata ,header=TRUE,sep=”\t”,dec=”,”,encoding=”latin1″,)
colnames(fdc@data)[1]<-“id”
names(donnees)[1]<-“id”
fdc@data <- data.frame(fdc@data, donnees[match(fdc@data[,”id”], donnees[,”id”]),])
fdc@data<-data.frame(fdc@data$id,fdc@data[,myvar])
colnames(fdc@data)<-c(“id”,”var”)
fdc@data$var[is.na(fdc@data$var)] <- 0
fdc@data$nbdots <-as.integer(fdc@data$var/nb)
dots<-dotsInPolys(fdc, as.integer(fdc@data$nbdots),f = “random”)
plot(fdc,col=”#CCCCCC”,border=”white”)
plot(dots, pch=20, coll=”black”,add=TRUE,cex=0.01)

Nicolas Lambert

Ingénieur de recherche CNRS en sciences de l'information géographique. Membre de l'UMS RIATE et du réseau MIGREUROP / CNRS research engineer in geographical information sciences. Member of UMS RIATE and the MIGREUROP network.

Vous aimerez aussi...

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *

Ce site utilise Akismet pour réduire les indésirables. En savoir plus sur comment les données de vos commentaires sont utilisées.

Rechercher dans OpenEdition Search

Vous allez être redirigé vers OpenEdition Search