r - r-raster中overlay和ifelse函数之间的冲突

Question

我正在运行一个代码并在使用 ifelse 时发现了一个奇怪的行为overlay 中的函数在光栅包中。简而言之，如果每个栅格的前5个值是NA，函数就会报错.
为什么会这样？
下面是一个简短的代码，它模拟了我在使用 R 3.2.3 和光栅版本 2.5-2 时发现的问题，以及我正在考虑使用的一些临时解决方案。

谢谢

卡洛斯·阿尔贝罗

library(raster)
cob1d <- raster(matrix(1,nr=6,nc=6))
cob1 <- cob1d; cob2 <- cob1d; cob3 <- cob1d
overlay(cob1, cob2, cob3, fun=function(x1, x2, x3)  ifelse(x1 > 0, x1 + x2 + x3,  x3))

# class       : RasterLayer 
# dimensions  : 6, 6, 36  (nrow, ncol, ncell)
# resolution  : 0.1666667, 0.1666667  (x, y)
# extent      : 0, 1, 0, 1  (xmin, xmax, ymin, ymax)
# coord. ref. : NA 
# data source : in memory
# names       : layer 
# values      : 3, 3  (min, max)


# Changing the first 5 values...

cob1[1:5] <- NA; cob2[1:5] <- NA; cob3[1:5] <- NA
overlay(cob1, cob2, cob3, fun=function(x1, x2, x3)  (x1 + x2 + x3))

给出了相同的结果...

# but if I use `ifelse`, there is a problem:

overlay(cob1, cob2, cob3, fun=function(x1, x2, x3)  ifelse(x1 > 0, x1 + x2 + x3,  x3))
# Error in ifelse(x1, x1 + x2 + x3, x3) : 
#  argument "x2" is missing, with no default 

# Another way to solve it is adding a useless extra variable without `NA`.

cob4 <- cob1d
overlay(cob1, cob2, cob3, cob4, fun=function(x1, x2, x3, x4)  ifelse(x1 > 0, x1 + x2 + x3,  x3))
# same result as before...
# class       : RasterLayer 
# dimensions  : 6, 6, 36  (nrow, ncol, ncell)
# resolution  : 0.1666667, 0.1666667  (x, y)
# extent      : 0, 1, 0, 1  (xmin, xmax, ymin, ymax)
# coord. ref. : NA 
# data source : in memory
# names       : layer 
# values      : 3, 3  (min, max)

# or just avoiding the use of the `ifelse` function...

overlay(cob1, cob2, cob3, cob4, fun=function(x1, x2, x3, x4)  (x1 > 0) * (x1 + x2 + x3) + (x1 <= 0)*x3)

Answer 1

这是一个有趣的案例。 overlay使用前 5 个单元格来确定如何处理数据(通过 apply 或通过 do.call)。第一个测试是看是否apply可以使用。使用的函数可以通过 apply 运行如果所有值都是 NA

f <- function(x1, x2, x3)  ifelse(x1 > 0, x1 + x2 + x3,  x3)
m <- matrix(NA, 5, 3)
apply(m, 1, f)

基于该测试， apply用于所有单元格。然而， apply当并非所有值都是 NA 时，此函数失败:

m[1] <- 1
apply(m, 1, f)

这失败了，因为只有一个第一个参数 x (不是 f 要求的三个参数)。

而相反的情况很常见(当仅提供 NA 值时函数失败)，这种情况很少见。
raster:::calc有论据 forcefun避免使用 apply ，但这在 overlay 中不可用.我已将其添加到 raster 的 future 版本中(版本 >= 2.5-4)，这样你就可以:

f <- function(x1, x2, x3)  ifelse(x1 > 0, x1 + x2 + x3,  x3)
overlay(cob1, cob2, cob3, fun=f, forcefun=TRUE)

在任何情况下，您都可以使用 calc作为解决方法:

s <- stack(cob1, cob2, cob3)
r <- calc(s, fun=function(x)  ifelse(x[1] > 0, x[1] + x[2] + x[3],  x[3]))

或另一种选择，基于您的替代功能:

r <- (cob1 > 0) * (cob1 + cob2 + cob3) + (cob1 <= 0)*cob3