如何在R中的SpatialPolygons对象中将多边形移动到其他列表?

时间:2019-06-22 21:54:42

标签: r polygon spatial sp

我有Large SpatialPolygons个对象rtu,其中包含161个多边形。在第一列表中有146 Polygons。我将一个多边形从同一对象(索引为92的列表)移至该列表(尝试进行追加)。最后,删除该行。

r = rtu@polygons[[92]]@Polygons
rp <- Polygons(r, 1)
ps = c(rtu@polygons[[1]]@Polygons, rp@Polygons)
rtu@polygons[[1]]@Polygons <- ps
rtu@polygons[[92]] <- NULL

此时,一切看起来都还不错。结构现在有147 Polygons,并且绘图正常。但是当我进一步尝试时:

over(rtu[1], rtu[2])

这会引发错误:

> ov = over(rtu[1], rtu[2])
Error in RGEOSBinPredFunc(spgeom1, spgeom2, byid, func) : 
  lengths of comment and Polygons slot differ

我想我的转换有时是错误的,但是在哪里呢?我猜想有一些功能可以通过这种转换来更好。请帮助我找到此错误。

0 个答案:

没有答案