如何将sfc_MULTIPOLYGON转换为sf/data.frame并计算面积？

sf::st_union()合并后返回的是单纯的几何对象（sfc类型），而dplyr::mutate()需要操作带数据框结构的sf对象，这就是报错的核心原因。只需把sfc对象转换为含data.frame属性的sf对象，就能正常计算面积。

具体步骤与代码示例 #

1. 模拟原始sf数据（仅作演示） #

library(sf)
library(dplyr)

# 创建带多个独立多边形的sf对象
original_sf <- st_sf(
  id = 1:3,
  geometry = st_sfc(
    st_polygon(list(rbind(c(0,0), c(1,0), c(1,1), c(0,1), c(0,0)))),
    st_polygon(list(rbind(c(1,1), c(2,1), c(2,2), c(1,2), c(1,1)))),
    st_polygon(list(rbind(c(0,1), c(1,1), c(1,2), c(0,2), c(0,1))))
  )
)

2. 合并后得到sfc_MULTIPOLYGON对象 #

union_sfc <- st_union(original_sf)
class(union_sfc) # 输出 "sfc_MULTIPOLYGON" "sfc"

3. 错误示例（直接用mutate会报错） #

# 此代码会触发错误：mutate需要data.frame/tibble，而非sfc对象
union_sfc %>% mutate(total_area = st_area(.))

4. 解决方案：转换为sf+data.frame结构 #

有两种常用方法：

方法1：用st_sf()直接包裹几何对象

union_sf <- st_sf(geometry = union_sfc)

方法2：结合tibble和st_as_sf()（支持添加额外字段）

union_sf <- tibble::tibble(geometry = union_sfc) %>% st_as_sf()

5. 正常计算面积 #

union_sf <- union_sf %>% mutate(total_area = st_area(geometry))

# 查看结果
print(union_sf)

关键说明 #

sf对象本质是带几何列的data.frame/tibble，而sfc只是单纯的几何向量集合。mutate()依赖数据框的结构来添加新字段，所以必须先完成类型转换。

内容的提问来源于stack exchange，提问作者tamarack