按槽从空间多边形对象中细分多边形

Question

我有一个“SpatialPolygons”对象，我想在@hole==FALSE处对多边形进行子集，并创建一个新对象。

使用下面的代码，一次可以预测得到1个多边形（如果@hole==TRUE，则没有），但我无法从对象中分割多个多边形。

newSP <- aDis[aDis@polygons[[1]]@Polygons[[1]]@hole == FALSE]

我认为我的问题在于我不知道如何“一般”引用第二个列表，即从列表（Polygon）中获取项目，其中@孔==FALSE。将第二个[[]]留空不起作用。

附言：抱歉，没有可重复的例子；我不确定如何以简单的方式重新创建此问题。

Answer 1

最简单的方法是使用< code > sappy :

hasHole <- sapply(
  aDis@polygons[[1]]@Polygons,
  slot,
  "hole"
)
aDis@polygons[[1]]@Polygons[!hasHole]

不过，作为最终用户，不得不求助于这种东西是不寻常的。我不熟悉这个特定的包——它的其他地方没有定义功能来帮助您访问这些插槽吗？