terra::union处理零范围SpatVector异常及优雅合并方案咨询

问题现象 #

当使用terra::union合并两个正常的SpatVector对象时，结果符合预期：

a = terra::vect(cbind(c(0, 0, 1, 1, 0), c(0, 1, 1, 0, 0)), type = "polygons", crs = "local")
terra::union(a, a)
all.equal(a, terra::union(a, a)) # 返回TRUE

但合并零范围SpatVector（比如所有顶点重合的多边形）时，结果会出现额外维度，导致与原对象不一致：

b = terra::vect(cbind(rep(0, 5), rep(0, 5)), type = "polygons", crs = "local")
terra::union(b, b)
all.equal(b, terra::union(b, b)) # 返回FALSE（要素长度不同）

原因分析 #

零范围几何体（无有效面积/边界的多边形）属于拓扑上的“空要素”，terra::union的核心逻辑是基于拓扑边界的合并计算。对于正常几何体，union会合并重叠区域并保留唯一要素；但零范围对象没有可参与拓扑计算的有效边界，合并时会被视为独立的空要素重复添加，最终导致结果的要素长度超出原对象。

解决方案 #

针对你的使用场景（合并多个SpatVector，零范围对象不产生任何影响），可以通过先过滤零范围对象，再执行合并的方式实现优雅处理。

步骤1：编写过滤空要素的辅助函数 #

先实现一个函数，筛选出SpatVector中的非空（有效范围）要素：

filter_non_empty <- function(x) {
  # 筛选出非空的要素（零范围要素会被标记为empty）
  non_empty_mask <- !terra::is.empty(x)
  x[non_empty_mask, ]
}

步骤2：安全合并函数 #

基于过滤函数，实现一个忽略零范围对象的安全union逻辑：

safe_union <- function(x, y) {
  x_clean <- filter_non_empty(x)
  y_clean <- filter_non_empty(y)
  
  # 如果其中一个对象过滤后为空，直接返回另一个非空对象（或空对象）
  if (nrow(x_clean) == 0) return(y_clean)
  if (nrow(y_clean) == 0) return(x_clean)
  
  # 对有效对象执行正常union
  terra::union(x_clean, y_clean)
}

步骤3：测试验证 #

# 测试1：正常对象 + 零范围对象
all.equal(a, safe_union(a, b)) # 返回TRUE，符合预期
# 测试2：零范围对象 + 零范围对象
safe_union(b, b) # 返回空的SpatVector，无多余要素
# 测试3：多个对象合并
c = terra::vect(cbind(c(1,1,2,2,1), c(0,1,1,0,0)), type = "polygons", crs = "local")
result <- safe_union(safe_union(a, b), c)
# 结果为a和c的union，b被忽略

如果需要合并多个SpatVector列表，还可以扩展为批量处理：

batch_safe_union <- function(vec_list) {
  # 过滤所有空对象
  cleaned_list <- lapply(vec_list, filter_non_empty)
  cleaned_list <- Filter(function(x) nrow(x) > 0, cleaned_list)
  
  # 批量合并
  if (length(cleaned_list) == 0) {
    return(terra::vect())
  }
  Reduce(safe_union, cleaned_list)
}

# 批量测试
batch_result <- batch_safe_union(list(a, b, c))

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