Standardized Precipitation Index (SPI)

[bling-P]

Hi, thanks for the great work on this project!

I’m currently looking into extreme heat and drought indices. As extreme heat and drought events have become more frequent in recent years, heat–drought indices such as the Standardized Precipitation Index (SPI) are becoming increasingly important.

In Python, I’ve found package such as climate_indices, gma, and standard-precip, but they all seem to only handle 1-D data. This means that for multiple stations I have to loop over each one, and for gridded (lat/lon) data the computation becomes very slow and resource-intensive.

Is there a more efficient or vectorized way to compute these indices for multiple stations or gridded datasets, without having to loop through everything?

Thanks in advance for any suggestions!