# Alcuni comandi per facilitare la costruzione di mappe delle # province e regioni italiane in R, attraverso il "package maps". # Si tenga presente che i dati di partenza si riferiscono al 1989, # e quindi le province di piĆ¹ recente istituzione non appaiono. # Comandi provati con "maps" versione 2.0-40. # A.Azzalini, 05 luglio 2004 (modificato 22.04.2005, 10.12.2008) # #------------------------- library(maps) data(italyMapEnv) #----- # esempio 1 it <- map('italy', fill = TRUE, col="grey96") map.axes() map.cities(country="Italy", capitals=1, col=2) map.cities(country="Italy", capitals=0, col=3) #------------ # definizioni delle regioni # it <- map('italy', fill = TRUE, col="grey96") val.d.aosta <- it$names[15] piemonte <- it$names[c(6,16,24,32,35,40)] lombardia <- it$names[c(4,8,9,10,12,20,25,28,29)] liguria <- it$names[c(44,51,46,42)] trentino_alto.adige <- it$names[c(1,5)] veneto <- it$names[c(2,11,14,17,18,23, 27, 30, 31, 33)] friuli_venezia.giulia <- it$names[c(3,7,13,22,19,21)] emilia_romagna <- it$names[c(34,36,37,38,39,41,43,47)] toscana <- it$names[c(45,48,49,50,53,54,55,57,59,63,66,62,69,73,74)] marche <- it$names[c(52,56,60,61)] umbria <- it$names[c(58,64)] abruzzo <- it$names[c(65,70,71,72)] molise <- it$names[c(76,79)] lazio <- it$names[c(67,68,75,78,80,92)] campania<- it$names[c(81,82,84,88,91,96,99,100,101)] basilicata <- it$names[c(89,97)] puglia <- it$names[c(77,83,94,98,102)] calabria <- it$names[c(104,106,111)] sicilia <- it$names[c(109,110,112:133)] sardegna <- it$names[c(85,86,87,90,93,95,103,105,107)] #----- # esempio 2 (eseguire questo dopo le assegnazioni precedenti) # map('italy', regions=veneto, fill = TRUE, col="gray90") map('italy', regions="Treviso", fill = TRUE, col="yellow", add=TRUE) map('italy', regions="Belluno", fill = TRUE, col="tan2", add=TRUE) map.cities(country="Italy", capitals=0, col=3) map('italy', regions=c(trentino_alto.adige, veneto, friuli_venezia.giulia), fill = TRUE, col="gray90") #----- # esempio 3 it <- map('italy', fill = TRUE, col="grey96") cat("selezionare una provincia con il pulsante del topo..\n") name <- identify(map("italy",add=TRUE)) print(name) print(which(it$names==name)) map('italy', regions=name, fill=TRUE, col=4, add=TRUE)