我有三个数据帧,我试图找到一个逐个单元格的意思。
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 .....
答案 0 :(得分:2)
这是一个奇怪的错误消息。如果在没有加载旧工作区(可能是陈旧的)的情况下重新启动R ,这种情况通常会消失。如果这是正在使用的unlink(".RData")
,请退出R而不保存,然后重新开始。
回答你的问题,是的,堆叠它们要容易得多。 E.g。
f <- list.files(pattern='tif$')
s <- stack(f)
x <- sum(s)