R生长期天数计算中多个波段的栅格代数

Question

因此，我正在根据北美的每日最低和最高温度构建一个 30 年平均生长期天数栅格。输入数据由 60 个栅格、30 个 tmax 和 30 个 tmin 组成。这些栅格中的每一个都包含 365 个波段（每天一个波段）。 GDD 的计算相当简单，它是((tmax-tmin))/2)-10。问题是我似乎无法找到将所有波段和栅格一起计算的方法。

我最初的计划是将所有 tmin 栅格和所有 tmax 栅格叠加在一起，然后将它们相减。但在这种情况下，只计算第一个波段（第 1 天）。有没有什么办法可以不用一一提取所有波段来进行操作。

library(raster)

tmax <- stack("C:/gddtmax")
tmin <- stack("C:/gddtmin")
average <- ((tmax-tmin)/2)-10

Answer 1

你的例子对我来说很好。

library(raster)
f <- system.file("external/rlogo.grd", package="raster")
tmax <- brick(f) / 10
tmin <- tmax - 10

tmax 和tmin 各有三层

avg <- ((tmax-tmin)/2)-10

而avg也有三层

avg
#class       : RasterBrick 
#dimensions  : 77, 101, 7777, 3  (nrow, ncol, ncell, nlayers)
#resolution  : 1, 1  (x, y)
#extent      : 0, 101, 0, 77  (xmin, xmax, ymin, ymax)
#coord. ref. : +proj=merc +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
#data source : in memory
#names       : red, green, blue 
#min values  :  -5,    -5,   -5 
#max values  :  -5,    -5,   -5 

编辑 看起来问题是由于内存管理错误（在版本 2.7-3 中修复）。这里有两个解决方法。

使用overlay（它还允许您指定输出文件名）

avg <- overlay(tmin, tmax, fun=function(x, y){ ((y-x)/2)-10 })

或使用rasterOptions

rasterOptions(todisk=TRUE)
avg <- ((tmax-tmin)/2)-10
rasterOptions(todisk=FALSE)