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