R语言【raster】——calc()：使用公式从另一个Raster*对象计算新Raster*对象的值

Package raster version 3.6-27

---

Description

使用公式从另一个Raster*对象计算新Raster*对象的值。

如果x是一个RasterLayer，那么fun就是一个典型的函数，它可以将单个向量作为输入，并返回相同长度的值向量（例如sqrt）。如果x是一个RasterStack或RasterBrick, fun应该在一个值向量上操作(（每个单元一个向量）。如果fun返回单个值（如sum）， calc返回一个RasterLayer;如果fun返回多个数字（如fun=quantile）， calc返回一个RasterBrick。

在许多情况下，可以用微积分实现的，也可以用更直观的“栅格代数”符号来实现（参见算术方法）。例如，r <- r * 2而不是r <- calc(r, fun=function(x){x * 2}，或者 r <- sum(s)而不是r <- calc(s, fun=sum)。然而，在大型数据集上使用复杂公式时，calc应该更快。使用calc可以设置输出文件名和文件类型首选项。

---

Usage

calc(x, fun, filename='', na.rm, forcefun=FALSE, forceapply=FALSE, ...)

---

Arguments

参数【x】：Raster* 对象。

参数【fun】：函数。

参数【filename】：字符。输出文件名（可选）

参数【na.rm】：删除NA值，如果“fun”支持（仅在将多层光栅对象汇总为光栅层时相关）。

参数【forecefun】：逻辑值。强制calc不使用fun运行；用于歧义函数和调试。

参数【foreceapply】：逻辑值。强制calc使用fun运行；用于歧义函数和调试。

参数【...】：其余传递给writeRaster的参数。

---

Details

有些函数的意图可能是不明确的。思考以下代码：

library(raster)

r <- raster(volcano)

calc(r, function(x) x * 1:10)

在本例中，单元格值以矢量化的方式相乘，并返回一个单层，其中第一个单元格与1相乘，第二个单元格与2相乘，第11个单元格与1相乘，以此类推。但也许意图是创建10个新层(x*1, x*2，…)？这可以通过使用参数forceapply=TRUE来实现：

calc(r, function(x) x * 1:10, forceapply=TRUE)

---

Value

一个 Raster* 对象。

---

Note

对于大型对象，calc将逐块计算值。这意味着要使fun的结果正确，它不应该依赖于立即访问_all_值。例如，要通过减去Raster*对象的平均值（对于每一层）来缩放它的值，对于Raster对象x，你不应该这样做：

calc(x, function(x)scale(x, scale=FALSE))

因为每个块的平均值可能是不同的。不如这样做：

m <- cellStats(x, 'mean')

x - m