使用R语言terra包extract函数提取结果为空的问题求助

核心问题：缓冲区创建顺序错误 #

你当前是先在原矢量坐标系下创建1km缓冲区，再投影到栅格坐标系。如果原矢量的坐标系是地理坐标系（单位为度），buffer(1000)会生成1000度的缓冲区，这显然不是你要的1km，最终导致缓冲区和栅格完全不重叠，返回空值/NA。

修复步骤 #

1. 调整缓冲区创建流程：先投影，再缓冲 #

先将原始矢量数据投影到栅格的坐标系（CDL数据通常采用UTM投影，单位为米），再创建1km缓冲区：

# 读取矢量
CLU_MN <- vect("~/clu_public_a_mn001.shp")
# 先投影到栅格的坐标系
CLU_MN_proj <- terra::project(CLU_MN, crs(ST_CDL))
# 再创建1km缓冲区（此时单位为米，符合需求）
CLU_MN_buffer <- terra::buffer(CLU_MN_proj, 1000)

2. 验证空间重叠 #

执行以下代码确认缓冲区与栅格是否重叠：

# 绘制栅格边界
plot(ext(ST_CDL), main="栅格与缓冲区范围")
# 叠加缓冲区
lines(CLU_MN_buffer, col="red")

如果红色缓冲区不在栅格范围内，说明矢量数据本身与栅格的空间范围不匹配，需要检查矢量的AOI是否对应明尼苏达州27001县。

3. 正确使用extract函数 #

• 若要获取每个缓冲区覆盖的所有栅格值（原始明细），直接调用extract即可：

# 返回每个缓冲区对应的栅格值列表
extract_result <- terra::extract(ST_CDL, CLU_MN_buffer)
# 查看第一个缓冲区的栅格值
extract_result[[1]]

• 若要统计每个缓冲区的土地利用类型频次，使用fun=table即可（无需weights=TRUE，若需加权占比可改用fun=sum结合权重参数）：

# 统计每个缓冲区的土地利用类型频次
type_count <- terra::extract(ST_CDL, CLU_MN_buffer, fun=table)
# 转换为易读的数据框
type_count_df <- do.call(rbind, lapply(type_count, as.data.frame))

4. 排查其他可能问题 #

• 检查栅格数据是否加载正确：ST_CDL的ncell()和ext()应返回合理值
• 验证缓冲区有效性：is.valid(CLU_MN_buffer)确保没有无效几何

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