Rで大陸地図をプロットするにはどうすればよいですか？

Question

国レベルでプロットマップを作成するソリューションはたくさんありますが、私の場合は大陸レベルで統計を印刷したいと考えています。

唯一のことは、国レベルのマップを使用して各大陸の国のリストを使用することですが、この種のマップには単純な解決策があるかどうかは疑問でした。私のアイデアを実現するために、それはこのようになります：

## 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大陸を含む地域レベルにデータをプロットするか、集約するための機能を持っています。以下のプロットを作成する必要があり

シンプルなスタート：

：

mapCountryData(sPDF, nameColumnToPlot='REGION') 

がで地域レベルの外観に国から独自のデータを集計する：

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

または7つの大陸モデルの

?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

アンディの答え@フォローアップ、あなたはとてもように各大陸内の国のポリゴンをマージすることができます：大陸シェイプファイルの

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