如何使用两个RasterStack进行计算？

Question

我有两个大型Rasterstack（每个6000个）。我想做一些简单的计算，比如

         NewRasterStack = RasterStack_1 + RasterStack2*(-3)

您能帮助我建议使用代码进行此计算吗？所有数据都使用R raster包处理，所以我正在寻找栅格计算功能。

提前谢谢！

Answer 1

如果堆栈具有相同的范围和分辨率，则应该起作用：

NewRasterStack = RasterStack_1 + RasterStack2*(-3)

这可能会更有效率：

nr <- overlay(RasterStack_1, RasterStack2, fun=function(x,y) x - 3 * y)