我想计算 R 中每月降雨量栅格的标准化降水指数 (SPI)。我在 SO（R 标准化降水指数 .nc 文件）中看到了一个与此相关的问题，但它没有任何答案。所以，我提出一个新问题。我从以下来源下载了每月降雨量数据 https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.03/cruts.1905011326.v4.03/pre/
我已经使用raster包在 R 中阅读了它。
library(raster)
file <- "cru_ts4.03.1901.2018.pre.dat.nc"
rain <- brick(file)
我知道它可以使用包的calc功能来完成，raster并且对于 SPI 计算，R 中有许多包可用（例如spi，、、SPEI）precintcon。我能够计算点数据的 SPI。但我无法使用该calc函数进行栅格数据的 SPI 计算。现在我的问题是如何在不同的时间尺度（1、2、3、6、12、24）上计算这个栅格上的 SPI？