我正在运行代码，并在光栅包中使用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)
这是一个有趣的案例。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（版本> = 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
本文收集自互联网，转载请注明来源。
如有侵权，请联系 [email protected] 删除。
我来说两句