在R语言中检测两幅栅格数据的变化区域问题求助

错误原因 #

你用了基础R的ifelse()函数直接处理栅格对象，而ifelse()仅支持向量/标量运算，无法识别raster包的栅格S4对象，导致which()函数收到非逻辑向量的参数，触发报错。

正确实现方法 #

以下是几种可行的解决方案（假设你已加载library(raster)）：

方法1：使用overlay()函数（推荐，专为多栅格运算设计） #

# 生成r3：目标区域设为1，其余为NA（可替换成0，按需求调整）
r3 <- overlay(r1, r2, fun = function(x, y) {
  ifelse(x == 1 & y != 1, 1, NA)
})

方法2：直接利用栅格的逻辑运算特性 #

栅格对象支持向量式逻辑运算，可先生成逻辑栅格再赋值：

# 先创建逻辑栅格：TRUE对应目标区域，FALSE对应其他
logic_r <- r1 == 1 & r2 != 1
# 将TRUE替换为标记值（比如0），FALSE替换为NA
r3 <- logic_r
r3[r3 == TRUE] <- 0
r3[r3 == FALSE] <- NA

方法3：使用calc()函数（需合并栅格栈） #

# 将r1和r2合并为栅格栈
rs <- stack(r1, r2)
# 对栈内栅格逐像元运算
r3 <- calc(rs, fun = function(x) {
  ifelse(x[1] == 1 & x[2] != 1, 0, NA)
})

注意事项 #

• 你已通过compareRaster(r1,r2)验证栅格属性匹配，这一步无需额外处理
• 标记值可按需调整：如果需要用0标记目标区域，把代码里的1替换成0即可

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