使用R中的spplot在图形上绘制多个shp文件

Question

我有3个shp文件，分别代表房子，房间和房子的床。我需要使用R在图表上绘制它们，以便它们彼此重叠。我知道在plot函数中，我可以使用line在现有绘图的顶部绘制新行，spplot中是否有相应的内容？感谢。

Answer 1

这是一种方法，使用 latticeExtra 包中的漂亮layer()函数：

# (1) Load required libraries
library(sp)
library(rgeos)        # For its readWKT() function
library(latticeExtra) # For layer()

# (2) Prepare some example data
sp1 = readWKT("POLYGON((0 0,1 0,1 1,0 1,0 0))")
sp2 = readWKT("POLYGON((0 1,0.5 1.5,1 1,0 1))")
sp3 = readWKT("POLYGON((0.5 0,0.5 0.5,0.75 0.5,0.75 0, 0.5 0))")

# spplot provides "Plot methods for spatial data with attributes",
# so at least the first object plotted needs a (dummy) data.frame attached to it.
spdf1 <- SpatialPolygonsDataFrame(sp1, data=data.frame(1), match.ID=1)

# (3) Plot several layers in a single panel
spplot(spdf1, xlim=c(-0.5, 2), ylim=c(-0.5, 2), 
       col.regions="grey90", colorkey=FALSE) +
layer(sp.polygons(sp2, fill="saddlebrown")) +
layer(sp.polygons(sp3, fill="yellow"))

或者，您可以通过spplot()的{​​{1}}参数获得相同的结果。 （指定sp.layout=可确保将“屋顶”和“门”绘制在灰色方块之上/之上，并将其作为first=FALSE的第一个参数。）

spplot()

Answer 2

您可以使用sp.layout中的spplot参数。或者，您可以使用ggplot2。一些示例代码（未经测试）：

library(ggplot2)
shp1_data.frame = fortify(shp1)
shp1_data.frame$id = "shp1"
shp2_data.frame = fortify(shp2)
shp2_data.frame$id = "shp2"
shp = rbind(shp1_data.frame, shp2_data.frame)

ggplot(aes(x = x, y = y, group = group, col = id), data = shp) + geom_path()

在ggplot2中，数据中的列链接到图中的图形比例。在这种情况下，x是x坐标，y是y坐标，group是data.frame shp中的一列，它指定一个点所属的多边形，{ {1}}是多边形的颜色。我使用的几何是col，它根据多边形输入data.frame绘制一系列线条。另一种方法是使用geom_path，它也支持填充多边形。