在 R 中绘制具有不同层的地图答案

【问题标题】：Plot Maps with different layers in R在 R 中绘制具有不同层的地图
【发布时间】：2014-09-11 12:34:23
【问题描述】：

我正在使用 spplot 在 R 中为德国地区绘制一些数据。 - 效果很好。

为了更容易理解/阅读地图或定位地点，我现在想添加联邦州的边界。地图来自German Borders

我使用 readShapeSpatial 加载数据，然后使用 spplot with layer 来绘制地图。但是，我不确定

krs <- readShapeSpatial("VG250_Kreise.shp")
bdl <- readShapeSpatial("VG250_Bundeslaender.shp")

krs@data$n <- cut(rnorm(402, 0, 1), breaks=5)
seq.red <- rev(sequential_hcl(n=nlevels(krs@data$n), h=0))
spplot(krs, "n", col=grey(.9), col.regions=seq.red, pretty=TRUE)+
layer(sp.polygons(bdl))

我不确定我需要添加 bdl 的哪一部分才能仅获得边框。

非常感谢您的帮助！

谢谢！

【问题讨论】：

注意：不要使用readShapeSpatial，使用rgdal 包和readOGR 函数——这将读入shapefile中被readShapeSpatial忽略的任何投影信息。

标签： r maps spatial

【解决方案1】：

您可以使用sp.layout 参数进行过度绘图，请参阅this relevant question。因此，在您的情况下，您将在代码中添加以下内容：

bunder<-list("sp.polygons", bdl)
spplot(krs, "n", col=grey(.9), col.regions=seq.red, pretty=TRUE,
       sp.layout=bunder)

这将为您提供以下等值线：

【讨论】：