GADM-Maps 跨国比较图

Question

也许由于我对 R 比较陌生，我在 http://www 上使用 gadm-Mapfiles 时遇到问题.gadm.org/。

我尝试绘制几个国家的地图并将它们相互比较（使用不同的颜色）。

这就是我所做的

library('sp')
##
load(url('http://biogeo.ucdavis.edu/data/gadm2/R/ARG_adm0.RData')) 
# loads an Object "gadm" with shape of Argentinia
arg <- gadm # is there a more convenient way to do this in one line?
load(url('http://biogeo.ucdavis.edu/data/gadm2/R/CHL_adm0.RData'))
# loads an Object "gadm" with shape of Chile
chl <-gadm
load(url('http://biogeo.ucdavis.edu/data/gadm2/R/BOL_adm0.RData'))
# loads an Object "gadm" with shape of Bolivia
bol <- gadm
##
spplot(c(arg, chl, bol))
# output: unable to find an inherited method for function "spplot", for signature "list"

这是我的问题：（

1. 这个问题可能是由我的新手引起的）是否有更方便的方法来加载形状文件？我发现一直重命名 gadm-Object 是非常愚蠢的。也许甚至有一种方法，R 只下载一次数据，然后将它们存储在工作区/本地某个地方？
2. 我如何说服 R 将所有这些国家绘制在一张地图上？

先感谢您！

[编辑]

一些不错的功能 在 Gavin Simpson 的帮助下，我能够创建一些不错的函数，将整个地图合并减少到一行：

## you will need the sp-package
library('sp')

## load a file from GADM (you just have to specify the countries "special part" of the file name, like "ARG" for Argentina. Optionally you can specify which level you want to have
loadGADM <- function (fileName, level = 0, ...) {
    load(url(paste("http://biogeo.ucdavis.edu/data/gadm2/R/", fileName, "_adm", level, ".RData", sep     = "")))
    gadm
}

## the maps objects get a prefix (like "ARG_" for Argentina)
changeGADMPrefix <- function (GADM, prefix) {
    GADM <- spChFIDs(GADM, paste(prefix, row.names(GADM), sep = "_"))
    GADM
}

## load file and change prefix
loadChangePrefix <- function (fileName, level = 0, ...) {
    theFile <- loadGADM(fileName, level)
    theFile <- changeGADMPrefix(theFile, fileName)
    theFile
}

## this function creates a SpatialPolygonsDataFrame that contains all maps you specify in "fileNames".
## E.g.: 
## spdf <- getCountries(c("ARG","BOL","CHL"))
## plot(spdf) # should draw a map with Brasil, Argentina and Chile on it.
getCountries <- function (fileNames, level = 0, ...) {
    polygon <- sapply(fileNames, loadChangePrefix, level)
    polyMap <- do.call("rbind", polygon)
    polyMap
}

当您找到此页面时，请务必阅读以下答案： https://stackoverflow.com/a/33264548/263589

Answer 1

最简单的解决方案是

library(raster)
bol <- getData('GADM', country='BOL', level=1)

plot(bol)

将 R 数据添加到 GADM 网站以支持此功能。另请注意，文件格式已更改，因此本页描述的其他功能不再起作用。

将不同国家结合起来

per <- getData('GADM', country='PER', level=1)
bp <- bind(bol, per)
plot(bp)

Answer 2

对于问题 1，这是 R，因此您可以推出自己的 load() 函数来执行您想要的操作，例如：

loadGADM <- function(file, ...) {
    load(file, ...)
    gadm
}

并将其用作

> ls()
character(0)
> loadGADM <- function(file, ...) {
+     load(file, ...)
+     gadm
+ }
> arg <- loadGADM(url('http://gadm.org/data/rda/ARG_adm0.RData'))
> ls()
[1] "arg"      "loadGADM"

：当您知道对象已加载时，这是一个本地解决方案将被称为 gadm - 您可以改进该函数以不需要它，例如：

loadGADM <- function(file, ...) {
    f <- load(file, ...)
    get(f)
}

它可以工作，因为 load() 返回已加载对象名称的字符串。

对于问题 2，您需要将三个 sp 对象 rbind() 放在一起，而不是连接它们。但是，这对这些对象不起作用，并且多边形 ID 不是唯一的：

> sa <- rbind(arg, chl, bol)
Error in validObject(res) : 
  invalid class "SpatialPolygons" object: non-unique Polygons ID slot values

我正在解决此问题，如果我找到解决方法，将会更新。 解决方案是更改多边形使用 spChFIDs() 获取 ID 槽值。在这里，我们将 "arg_" 等附加到对象的行名称中，这样它们并不都是唯一的：

arg <- spChFIDs(arg, paste("arg", row.names(arg), sep = "_"))
chl <- spChFIDs(chl, paste("chl", row.names(chl), sep = "_"))
bol <- spChFIDs(bol, paste("bol", row.names(bol), sep = "_"))
sa <- rbind(arg, chl, bol)

然后我们可以绘制组合的 sp 对象：

plot(sa) ## beware might take a long time to plot...

Answer 3

我可以提供一些小的改进。

首先，如果您想绘制根据某些标准对国家/地区进行着色的地图图形，则必须在相应的 GADM 文件中加载该值，例如从另一个 data.frame 加载该值。我是这样做的：

## load a file from GADM (you just have to specify the countries "special part" of the file name, like "ARG" for Argentina. Optionally you can specify which level you want to have

loadGADM <- function (fileName, level = 0, ...) {
    load(paste("K:/Rdata/gadm/", fileName, "_adm", level, ".RData", sep = ""))
    gadm$coef <- land$coef[land[["abr"]]==fileName]
    gadm
}

如果您绘制较大区域的地图，您的地图将会扭曲。使用 rgdal 包（及其依赖项）应用地图投影：

newproj <- CRS("+proj=wintri ellps=WGS84 +lon_0=15E")
spdf.wintri <- spTransform(spdf, newproj)

为了绘制属性 coef，sp 包提供了一个很好的方法方法：

karte <- spplot(spdf.wintri,"coef",
        xlim=c(-13,46),ylim=c(33,72),
        col.regions = rainbow(100, start = 2/6, end = 4/6),
        main="Country Dummies")