Recomendadtion for working with raster data

[AmandaBatlle]

Hello,

I have a time series of 30 years of Andalucia Region. What would be the best/most efficinet way to generate SPI in raster format?
I saw the example using 3-d array input:
data(cruts4)
dim(cruts4)
spei_12 <- spei(cruts4, 12)

dim(spei_12$fitted)
[1] 1440 2 3

But I'm not sure how the multidimensional SPEI output can be slices and exported in a format that can be converted back to raster.

Thank you very much.

[sbegueria]

Hi Amanda, I don’t have example code to do that right away, but if you want to export slices of the 3D array to make maps (to use then in a GIS, for instance), I recommend you have a look at the sf and terra packages.