Come posso tracciare una mappa dei continenti con R?

Question

Esistono molte soluzioni per tracciare mappe a livello di paese, ma nel mio caso voglio stampare statistiche a livello di continente.

L'unica cosa che mi viene in mente è utilizzare le mappe a livello di paese e utilizzare un elenco di paesi per ogni continente, ma mi chiedevo se esiste una soluzione semplice per questo tipo di mappe. Per implementare la mia idea sarebbe come questo:

## produce the world map 
map() 
## list of countries per continent 
SA <- c("argentina", "bolivia", "brazil", "chile", "colombia", "ecuador", "guyana", "paraguay", "peru", "suriname", "uruguay", "venezuela") 
map(regions = SA, fill=TRUE, add=TRUE) 

Answer 1

rworldmap ha funzioni per la stampa o l'aggregazione dei dati a livello regionale tra i continenti.

Un inizio semplice che dovrebbe produrre la trama di seguito:

library(rworldmap) 
#get coarse resolution world from rworldmap 
sPDF <- getMap() 
#mapCountries using the 'continent' attribute 
mapCountryData(sPDF, nameColumnToPlot='continent') 

O per il modello 7 continenti:

mapCountryData(sPDF, nameColumnToPlot='REGION') 

ad aggregare i propri dati da paese a guardare livello regionale a:

?mapByRegion 

Answer 2

library(sp) #Load your libraries 
library(maptools) 
#Download the continents shapefile 
download.file("http://baruch.cuny.edu/geoportal/data/esri/world/continent.zip", 
       "cont.zip") 
#Unzip it 
unzip("cont.zip") 
#Load it 
cont <- readShapeSpatial("continent.shp") 
#Plot it 
plot(cont, 
    col=c("white","black","grey50","red","blue","orange","green","yellow")) 
#Or any other combination of 8 colors 

Answer 3

Facendo seguito a @ risposta di Andy, è possibile unire i poligoni di campagna all'interno di ogni continente in questo modo:

library(rworldmap) 
library(rgeos) 
library(maptools) 
library(cleangeo) ## For clgeo_Clean() 

sPDF <- getMap() 
sPDF <- clgeo_Clean(sPDF) ## Needed to fix up some non-closed polygons 
cont <- 
    sapply(levels(sPDF$continent), 
      FUN = function(i) { 
       ## Merge polygons within a continent 
       poly <- gUnionCascaded(subset(sPDF, continent==i)) 
       ## Give each polygon a unique ID 
       poly <- spChFIDs(poly, i) 
       ## Make SPDF from SpatialPolygons object 
       SpatialPolygonsDataFrame(poly, 
             data.frame(continent=i, row.names=i)) 
      }, 
      USE.NAMES=TRUE) 

## Bind the 6 continent-level SPDFs into a single SPDF 
cont <- Reduce(spRbind, cont) 

## Plot to check that it worked 
plot(cont, col=heat.colors(nrow(cont))) 

## Check that it worked by looking at the SPDF's data.frame 
## (to which you can add attributes you really want to plot on) 
data.frame(cont) 
#     continent 
# Africa    Africa 
# Antarctica  Antarctica 
# Australia   Australia 
# Eurasia    Eurasia 
# North America North America 
# South America South America