使用R语言disagg命令匹配栅格分辨率时遇类型错误求解决方案

错误原因 #

disagg()函数的第二个参数需要传入整数型的拆分因子（fact），代表粗分辨率栅格的每个单元格要拆分成多少个细分辨率单元格。你直接传入了完整的栅格对象dataset2，导致S4类型对象无法转换为函数所需的数值类型，从而触发类型不兼容错误。

两种可行解决方案 #

方案一：计算拆分因子后使用disagg

先根据两个栅格的分辨率计算拆分倍数，取整后传入fact参数：

# 加载raster包（如果未加载）
library(raster)

# 计算x、y方向的拆分因子，取整为整数
fact <- round(res(dataset1) / res(dataset2))

# 执行拆分，指定最近邻方法
harmonised_data <- disagg(dataset1, fact = fact, method = "near", filename = "merged_data.tif")

注意：需确保两个栅格的空间范围（extent）完全一致，否则拆分后的结果可能与dataset2的范围不匹配。若范围不一致，需先使用crop()函数将dataset1裁剪到dataset2的范围。

方案二：用resample()直接匹配目标栅格（更推荐）

如果需要严格对齐dataset2的空间范围和分辨率，resample()函数更稳妥，它会自动处理范围和分辨率的匹配，且支持最近邻方法：

library(raster)

# 使用最近邻法（ngb = nearest neighbor）重采样到dataset2的规格
harmonised_data <- resample(dataset1, dataset2, method = "ngb", filename = "merged_data.tif")

额外提示 #

• 若你使用的是terra包的SpatRaster对象，语法类似：
  • 用disagg：fact <- round(terra::res(dataset1) / terra::res(dataset2))，再terra::disagg(dataset1, fact = fact, method = "near")
  • 用resample：terra::resample(dataset1, dataset2, method = "near")
• 输出文件名建议加上.tif等明确后缀，避免格式识别问题。

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