Question

在土壤制图的框架中，我需要总结未确定的数量的栅格。我尝试使用'raster'包和'do.call'函数来完成它。但是，如果'sum'函数可以总计多个栅格，则使用do.call执行相同的操作会导致错误。我做错了什么？

library(raster)

r1 <- raster(ncol=10, nrow=10)   # dataset for test
values(r1) <- runif(ncell(r1))
r2 <- raster(ncol=10, nrow=10)
values(r2) <- runif(ncell(r2))
r3 <- raster(ncol=10, nrow=10)
values(r3) <- runif(ncell(r3))

sum(r1,r2,r3)    # works nice

do.call(sum,list(r1,r2,r3))

##Erreur dans as.character(sys.call()[[1L]]) : 
##cannot coerce type 'builtin' to vector of type 'character'

谢谢你的帮助，

弗朗索瓦

Answer 1

您可以使用Reduce和+来计算列表中的总和：

Reduce("+",list(r1,r2,r3))
class       : RasterLayer 
dimensions  : 10, 10, 100  (nrow, ncol, ncell)
resolution  : 36, 18  (x, y)
extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 
data source : in memory
names       : layer 
values      : 0.4278222, 2.476625  (min, max)

至于为什么你的原始命令不起作用，这有点令人困惑。将函数名称作为字符提供似乎有效：

do.call("sum",list(r1,r2,r3))
class       : RasterLayer 
dimensions  : 10, 10, 100  (nrow, ncol, ncell)
resolution  : 36, 18  (x, y)
extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 
data source : in memory
names       : layer 
values      : 0.4278222, 2.476625  (min, max)

但在其他情况下不需要这样做：

do.call(sum,list(1,2,3))
[1] 6

Answer 2

我不知道为什么这不起作用（如James所指出的那样，没有关于sum的引号），也许这是与“sum”成为S4 Summary组泛型的成员相关的bug（或特征）;其他成员如“max”和“prod”具有相同的行为。

无论哪种方式，而不是

 do.call("sum", list(r1,r2,r3))

你也可以

 sum(stack(r1,r2,r3))

或者如果您已经有一个列表

 sum(stack(list(r1,r2,r3)))

使用do.call和raster包对未确定数量的栅格求和

2 个答案: