我怎样才能找到三个栅格的平均值?

时间:2015-11-13 02:28:30

标签: r spatial raster

我有三个数据帧,我试图找到一个逐个单元格的意思。

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 .....

1 个答案:

答案 0 :(得分:2)

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

回答你的问题,是的,堆叠它们要容易得多。 E.g。

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