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

Question

我正在使用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部分才能获得边框。

我非常感谢你的帮助！

谢谢！

Answer 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)

这会给你以下等级：