Question

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

由于

Carlos Albero

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：{/ p>，则NA会失败并显示此功能

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

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

虽然相反的情况是常见的（当仅提供NA值时失败的函数），但这种情况很少见。

raster:::calc使用参数forcefun以避免使用apply，但overlay中无此功能。我已将其添加到raster的未来版本（版本＆gt; = 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

Answer 2

也许你必须使用每一个＆＃34; x＆＃34;在你的函数中有一个条件，例如

overlay(cob1, cob2, cob3, fun=function(x1, x2, x3)  ifelse(x1 > 0 & x2 >0 & x3 >0 , x1 + x2 + x3,  x3))

不仅是x1（x1>0）

的条件

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

2 个答案: