最近,我需要计算 12 个世界气候栅格图层的平均值。我们有两种方法可以做到。第一个更直接:
mean.layer <- mean(L1, L2,......,L12) # L1 means the first layer
或者
mean.layer <- (L1+L2+......+L12)/12
另一个对我来说是新的:
layer.stack <- stack(L1,L2,......,L12)
mean.layer <- calc(layer.stack, mean, na.rm = T)
有人可以解释使用calc
和stack
代替mean
直接使用函数的优点吗?据我所知,我们可以直接以相同的分辨率和扩展名操作栅格数据。
2021.7.10 编辑。我重写了第二种方法来纠正一些错误输入。
请始终包含一些示例数据
library(raster)
b <- stack(system.file("external/rlogo.grd", package="raster"))
这两个语句是等价的
x <- mean(b)
y <- calc(b, mean)
但是calc
有一个文件名参数,因此您可以一步将结果保存到磁盘。
calc
在处理大型栅格和更复杂的函数时尤其有利。例如
z1 <- calc(b, function(i) 100 - sqrt(mean(i + 10)))
相当于这个
z2 <- 100 - sqrt(mean(b + 10))
但是后者可能需要写 4 个临时文件来存储这些值,而前者只需要 1 个这样的文件。
你应该避免像这样获得平均值的方法
mean.layer <- (L1+L2+......+L12)/12
编写起来很麻烦,容易出错,而且不能缩放(想象一下为 1200 个栅格做这件事!)。
我不知道你在哪里找到了calc
你的问题中的公式,但这没有意义,也行不通。
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句