仅保留sf数据中的二维分量。框架

我的设置与中类似 this question . 然而,就我而言,我有一个 sf 数据包含不同几何体类型混合的框架,一些 POLYGON s、 一些 GEOMETRYCOLLECTION 是这样的:

a <- st_polygon(list(cbind(c(0,0,7.5,7.5,0),c(0,-1,-1,0,0))))
b <- st_polygon(list(cbind(c(0,1,2,3,4,5,6,7,7,0),c(1,0,.5,0,0,0.5,-0.5,-0.5,1,1))))
i <- st_intersection(a,b)

a1 <- st_sf(a=1, geom = st_sfc(i))
a2 <- st_sf(a=2, geom = st_sfc(a))

ii <- rbind(a1, a2)

正如上面提到的问题,我只想保留 几何集合 st_area() 无法继续 几何集合 s

然而,在混合几何的情况下,给出了答案

st_cast(ii)[which(st_is(st_cast(ii), c("POLYGON", "MULTIPOLYGON"))),]

不起作用 st_cast() 离开 几何集合 未触及。