使用 ifelse 在 r 中叠加栅格以满足三种不同的条件

Question

我正在尝试使用栅格运行叠加函数，其中我想在每个像元（使用 == 和 > 或 < 运算符）满足所有 3 个不同条件并生成单个栅格作为输出。

使用 & 运算符运行ifelse似乎以从左到右的线性方式查看条件 - 如果满足前两个条件，那么它将生成 if 条件作为输出，而不管第三个条件如何。 此处不能使用&&，因为结果未矢量化。

您可以在下面的示例中看到这一点，其中对于生成的栅格，很明显它没有评估所有三个参数。第一个线索是，即使第三个栅格中的某些值为 NA，它也会产生结果。

我想我可以通过首先检查 x 和 y 的条件，然后用这个结果分别检查 z 的条件来获得我想要的结果，但我希望能够在一个函数中完成这一切（似乎这应该是可能的，至少）。

希望有人能为我指出正确的方向。

library(raster)
fn <- system.file("external/test.grd", package="raster")
  s <- stack(fn, fn,fn)
  #Create grids
  s[[1]] <- round(runif(ncell(s), 1, 2))
  s[[2]] <- round(runif(ncell(s), 1, 2))
  s[[3]] <- round(runif(ncell(s), 1, 2))
  #convert some values in s[[3]] to NA
  s[[3]][s[[3]] == 1]<- NA
  #run overlay function
  result.rast <- overlay(s[[1]], s[[2]], s[[3]], fun =   
  function(x,y,z) { 
        ifelse( x == 2 & y == 1 & z ==2, 1, 0) 
  } )

Answer 1

我没有看到第三个条件未被使用的证据。NA 值是一种特例。请参阅函数f2，了解您可以执行的一些操作。

使用

较小的栅格更容易看到

发生了什么

library(raster)
set.seed(0)
r <- raster(ncol=10, nrow=10, xmn=0, xmx=10, ymn=0, ymx=10)
r1 <- setValues(r, round(runif(ncell(r), 1, 2)))
r2 <- setValues(r, round(runif(ncell(r), 1, 2)))
r3 <- setValues(r, round(runif(ncell(r), 1, 2)))
r3[r3 == 1] <- NA
s <- stack(r1, r2, r3)

res1 <- overlay(s, fun =   
  function(x,y,z) { 
        ifelse( x == 2 & y == 1 & z ==2, 1, 0) 
  } )

#A more complex function, that keeps NAs
f2 <- function(x,y,z) {
    a <- rep(0, length(x))
    a[x == 2 & y == 1 & z ==2] <- 1
    a[is.na(x) | is.na(y) | is.na(z)] <- NA
    a 
} 
res2 <- overlay(s, fun = f2)

Answer 2

你试过stackApply吗？

您还可以将栅格堆栈的每个图层用作矢量。下面是一个示例（不过，这可能是引用 rasterStack 中的像元的更好方法）

tt <- raster(ncol=4,nrow=5)
tt[] <- 1
tts <- stack(tt,tt,tt)
tts[[1]][4,2]<-NA
# now the condition
tt2 <- (tts[[1]] == 1 & tts[[2]] == 1 & tts[[3]] == 1)
plot(tt2)