使用R/terra处理大型SpatVector时intersect/crop返回Large list问题

背景与目标 #

• 资深ArcMap用户，熟练掌握R及数据处理，刚接触Terra库，需自动化ArcGIS工作流完成多场景计算
• 数据集：约2000种澳大利亚（含海洋）物种的多部件多边形分布数据，分解后共15万个单个多边形，分辨率为1km或10km
• 核心目标：计算指定网格每个单元的物种相关指标（如物种数量、物种覆盖面积等），输出包含单元内物种分布信息的data.frame

遇到的问题 #

• 使用rasterize()处理聚合数据时，无法得到正确的物种计数（推测因大量微小多边形导致，参考相关issue未解决）
• 改用intersect()方案时，未得到预期的SpatVector，反而返回包含数千元素的Large list；对数据执行crop()操作时也出现相同问题
• 已尝试多种数据修复操作（空间子集化、移除问题物种、分解/聚合多边形、修复几何等），但crop或intersect步骤仍返回Large list

示例代码 #

基础流程代码 #

library(terra)
library(dplyr)

# 加载示例数据
p <- vect(system.file("ex/lux.shp", package="terra"))
v <- p
cell_size <- 0.1 # 单位：十进制度

# 创建网格单元
r <- rast(v, res=cell_size)
# 为每个单元格分配ID并命名
values(r) <- 1:ncell(r)
names(r) <- "CELLID"
# 将栅格网格转换为多边形网格
z <- as.polygons(r)
# 与物种数据执行相交操作
u <- intersect(z,v)

后续分析代码 #

# 按单元格ID和物种名称聚合，得到每个单元格的唯一物种列表
pa <- aggregate(u, by=c("CELLID","NAME_1"))

数据说明 #

实际数据为澳大利亚环境部发布的物种分布数据集

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