基于土地覆被条件与矢量多边形提取SOC均值的技术咨询

1. 矢量范围差异的影响及对齐方法 #

• 影响：若矢量与栅格范围不匹配，提取操作只会处理两者重叠区域，矢量超出栅格的部分因无对应栅格数据会被忽略；栅格超出矢量的部分不影响结果，但会增加不必要的计算量。
• 对齐步骤：
  1. 先确认矢量与栅格的坐标参考系（CRS）完全一致，不一致时用project(poly, crs(raster))转换矢量CRS。
  2. 用栅格的范围裁剪矢量，确保两者范围匹配：

     poly <- crop(poly, ext(raster_object))
     

2. 矢量是否需要转为栅格？ #

不需要，直接保留SpatVector格式即可：

• 转栅格会引入多边形边界的像素化误差，导致提取结果不准确；
• terra包的extract()函数支持直接基于矢量多边形提取栅格统计值，效率和精度都更优。

3. 实现草地像素平均SOC提取的代码 #

假设你的土地覆被栅格中，草地对应的数值为2（请替换为实际值），完整代码如下：

library(terra)

# 加载数据（替换为你的文件路径）
land_cover <- rast("land_cover_30m.tif")
soc_resampled <- rast("soc_30m.tif")
polygons <- vect("your_polygons.shp")

# 统一CRS（若矢量与栅格CRS不同）
if (!identical(crs(polygons), crs(land_cover))) {
  polygons <- project(polygons, crs(land_cover))
}

# 对齐矢量与栅格范围
polygons <- crop(polygons, ext(land_cover))

# 创建仅保留草地区域的SOC栅格（非草地设为NA）
soc_grass <- mask(soc_resampled, land_cover == 2)

# 提取每个多边形内草地的平均SOC
mean_soc <- extract(soc_grass, polygons, fun = mean, na.rm = TRUE)

# 将结果合并到矢量属性表
polygons$mean_soc_grass <- mean_soc[, 2]

# 导出结果（可选）
writeVector(polygons, "polygons_with_grass_soc.shp")

代码说明： #

• mask()函数：把非草地对应的SOC像素设为NA，后续计算均值时会自动忽略这些无效值；
• extract()的fun=mean指定计算均值，na.rm=TRUE确保忽略NA值（即非草地像素）；
• 最终结果会添加到矢量的属性列mean_soc_grass中，可直接查看或导出为矢量文件。

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