在R中，我如何加入和子集SpatialPolygonsDataFrame？

Question

我试图找出在R中如何执行(在地理信息系统中很容易)操作的方法。

让我们以spdep包中的多边形数据集为例

library("spdep")
c <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1])
plot(c)

我已经找到了可以使用subset选择逻辑语句的多边形。例如：

cc <-子集(c，c@data$POLYID<5)图(Cc)

现在，假设我有另一个数据框架，我想加入到我的空间数据：

POLYID=1:9
TO.LINK =101:109
link.data <- data.frame(POLYID=POLYID, TO.LINK=TO.LINK)

使用这两个数据集，如何获得两个空间数据帧：

1. 首先，由在第二个数据帧中有ID的多边形组成。
2. 第二，由第二数据帧中不存在的相反集多边形组成.

我怎么能说到这一点？

Answer 1

这可能会成功的。首先，您需要相关的ID。

myIDs <- link.data$POLYID

然后，如您所指出的那样，使用subset：

subset(c, POLYID %in% myIDs)
subset(c, !(POLYID %in% myIDs))

请注意，这假设您的第一个dataframe，c，也有一个相关的列名为POLYID。