0

我正在使用rworldmap创建世界地图来显示数据。就我的数据而言，小岛国是相当重要的，但它们太小了，以至于无法看到整个世界。所以我想在岛国的位置有一个小的可见点/圆圈，以适当类别的颜色代替。我环顾四周，但没有找到解决方案。当然，我可以为每个岛国手动执行此操作。我的问题是，有没有更快的解决方案？

作为一个例子，这里是用于绘制数据（在本例中为生物多样性）的标准示例，该示例随rworldmap软件包一起提供。我怎么能在这张地图上显示岛国？任何建议将不胜感激！

library(rworldmap)
mapCountryData()
data("countryExData",envir=environment(),package="rworldmap")
sPDF <- joinCountryData2Map(countryExData
              , joinCode = "ISO3"
              , nameJoinColumn = "ISO3V10"
              )
mapCountryData( sPDF
              , nameColumnToPlot="BIODIVERSITY" 
              )
4

1 回答 1

0

感谢您lmo为我指明了正确的方向，我现在找到了解决方案。我将在此处提供答案，包括一个可重现的示例。我还认为我必须解释一个我没有说清楚的问题。我希望绘制世界上所有国家以获取世界地图。然后，我想为岛国添加额外的点以使其可见。我找到的解决方案是首先绘制正常的世界地图，然后在第一个绘图中添加第二张地图，但仅适用于所建议的选定岛国lmo。

首先，这里是再次生成基本图的代码

par(mar=c(0,0,1,0))
data("countryExData",envir=environment(),package="rworldmap")
sPDF <- joinCountryData2Map(countryExData, joinCode = "ISO3", nameJoinColumn = "ISO3V10")
mapCountryData( sPDF, nameColumnToPlot="BIODIVERSITY", catMethod="fixedWidth")

在这里，这些岛屿是不可见的，它们太小了。现在我必须只为这些岛国生成数据集。它们不在数据集中，所以为了演示我想要做什么，我选择了几个岛国并重用了其他数据集中的一些数据。所以这里的数据不是实际数据。

mapdata2 <- data.frame(ccode=c("ATG", "COM", "CPV", "DMA", "FJI", "FSM",
     "GRD", "KIR", "KNA", "LCA", "MDV", "MHL", "MUS", "NRU","PLW", "SLB", 
     "STP", "SYC", "TON", "TUV", "VCT", "VUT", "WSM"), 
     biodiv=countryExData$BIODIVERSITY[1:23], size=1)

最后需要大小变量，以便所有岛点具有相同的大小。现在我可以使用mapBubbles命令将所选岛国的点添加到原始绘图中，如下所示。

sPDF2 <- joinCountryData2Map(mapdata2, joinCode = "ISO3", nameJoinColumn = "ccode")
mapBubbles(sPDF2,nameZSize='size', nameZColour='biodiv',add = T, 
     addColourLegend = F, addLegend = F, pch=21, symbolSize=.22,catMethod="fixedWidth")

这产生了下面的地图。可能有更好或更快的方法，但我现在对结果很满意。当然，现在可以更改颜色、符号大小等...

在此处输入图像描述

于 2016-09-27T10:06:35.640 回答