如何找到三个栅格的平均值?

How can I find the mean of three rasters?

我有三个数据框,我正试图为其找到逐个单元格的平均值。

r1<-raster('a.tif')
r2<-raster('b.tif')
r3<-raster('c.tif')

但是,这样做会给我以下错误

q<-mean(r1,r2,r3)

q<-(r1+r2+r3)/3

错误

Error in .local(.Object, ...) : options(warn) not set
Warning message:
closing unused connection 4 .....

这是一个奇怪的错误消息。如果您重新启动 R 而没有 加载旧的工作区(可能是陈旧的),这种情况通常会消失。如果这是正在发生的事情,请使用 unlink(".RData"),不保存退出 R 并重新开始。

回答你的旁白问题,是的,将它们堆叠起来要容易得多。例如

f <- list.files(pattern='tif$')
s <- stack(f)
x <- sum(s)