如何在R的sf包中高效转换多层嵌套GeometryCollection为MultiPolygon

这是Stack Overflow问题《How to simplify a MULTISURFACE, CURVEPOLYGON, or COMPOUNDCURVE geometry?》的延伸问题，场景如下：

• 把MultiSurface转成GeometryCollection后，用lapply能处理单层嵌套的CURVEPOLYGON，成功转成MultiPolygon
• 但碰到多层嵌套的GeometryCollection（比如里面的CURVEPOLYGON还包含COMPOUNDCURVE）时，直接用lapply会报错
• 现在需要找到能高效转换这两种结构为MultiPolygon的方法，同时确认是不是得根据嵌套结构调整lapply的用法

可复现代码示例 #

library(sf)

# 构造多层嵌套的几何示例
# 先做一个包含复合曲线的CurvePolygon
compound_curve <- st_compound_curve(list(
  st_linestring(rbind(c(0,0), c(1,0))),
  st_circularstring(rbind(c(1,0), c(1,1), c(0,1)))
))
curve_polygon <- st_curve_polygon(list(compound_curve))
# 包装成MultiSurface再转成GeometryCollection
multi_surface <- st_multisurface(list(curve_polygon))
geom_collection <- st_geometrycollection(list(multi_surface))

# 尝试单层转换（会报错）
tryCatch({
  converted <- lapply(geom_collection, function(x) st_cast(x, "MULTIPOLYGON"))
}, error = function(e) {
  cat("转换报错：", e$message, "\n")
})

解决方案：用递归函数处理嵌套结构 #

单层lapply只能处理一层嵌套，碰到多层结构就不行了，得用递归函数遍历所有层级的几何元素，根据类型针对性转换：

# 递归转换函数
recursive_cast <- function(geom) {
  if (inherits(geom, "sfc_GEOMETRYCOLLECTION")) {
    # 遍历集合内每个元素递归转换，再合并成MultiPolygon
    converted_parts <- lapply(geom, recursive_cast)
    do.call(st_multipolygon, lapply(unlist(converted_parts, recursive = FALSE), st_geometry))
  } else if (inherits(geom, "sfc_CURVEPOLYGON")) {
    # CurvePolygon直接转MultiPolygon
    st_cast(geom, "MULTIPOLYGON")
  } else if (inherits(geom, "sfc_COMPOUNDCURVE")) {
    # 先把复合曲线转成LineString，再构建Polygon和MultiPolygon
    line <- st_cast(geom, "LINESTRING")
    polygon <- st_polygon(list(line))
    st_multipolygon(list(polygon))
  } else if (inherits(geom, "sfc_MULTISURFACE")) {
    # MultiSurface先转成GeometryCollection再递归处理
    recursive_cast(st_cast(geom, "GEOMETRYCOLLECTION"))
  } else {
    # 已经是目标类型或其他兼容类型，直接返回
    geom
  }
}

# 测试转换
result <- recursive_cast(geom_collection)
print(result)

关键说明 #

• 递归函数会自动识别所有嵌套层级，不管是单层还是多层结构都能处理，不用手动调整lapply的嵌套深度
• 针对不同几何类型做了适配：复合曲线先转成直线串再构建多边形，CurvePolygon直接转换，MultiSurface先转集合再递归
• 最终输出统一为MultiPolygon类型，满足后续处理需求

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