0

我有一个与这个问题类似的设置。然而，就我而言，我有一个sfdata.frame ，其中包含不同几何类型的混合，有些是这样POLYGON的，有些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)

就像上面提到的问题一样，我希望只保留那些GEOMETRYCOLLECTION二维的部分，因为最终我对这些几何形状的区域感兴趣并且st_area()不会在GEOMETRYCOLLECTIONs 上工作。

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

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

不起作用，因为st_cast()叶子GEOMETRYCOLLECTION未受影响。

4

2 回答 2

2

你试过了吗

st_collection_extract(ii)

? 你说的st_area不从事收藏是什么意思？我懂了

 > st_area(ii)
[1] 0.625 7.500
于 2017-11-13T19:57:03.237 回答
0

我不确定这是最优雅的解决方案，但它有效：

split(ii, 1:nrow(ii)) %>% 
  purrr::map(~ st_cast(.x) %>% 
                 filter(st_dimension(.) == 2)) %>% 
  do.call(rbind, .)
于 2017-11-07T20:49:37.473 回答