标准化降水蒸散发指数 SPEI 算法及代码（Java）

前言

关于SPEI

SPEI：Standardized Precipitation Evapotranspiration Index 标准化降水蒸散发指数。该指数是使用降水和气温数据计算得到的表征干湿状态的指数，可以计算多种时间尺度（1，3，12，24月等等），SPEI值越大越湿润，值越小越干旱，在干旱研究中应用广泛。它是从SPI Standardized Precipitation Index 归一化降水指数进一步发展而来，同时纳入蒸散发的影响，在气温有显著变化趋势的地区更加适用，尤其适合长时间尺度的研究。
关于SPEI发展背景及详细描述，可以查阅文献[1]

SPEI算法

这里尽量详细讲解一下SPEI的计算过程。
首先，SPEI指数，使用一种简单的方法计算潜在蒸散发PET，计算公式如下：

$$PET=16K(\frac{10T}{I}{)}^m$$
其中，$T$是月均温，单位摄氏度；$I$是热量指数，该指数通过12个月的月度指数相加得到，计算公式如下：
$$I=\sum _{mon=1}^{12}i=\sum _{mon=1}^{12}(\frac{T_{mon,ave}}{5}{)}^{1.514}$$

$T_{mon,ave}$为某月研究时期内多年均值，如一月多年均值为$T_{1,ave}$。$m$是一个取决于热量指数的系数：
$$m=6.75\times 10^{-7}{I}^3-7.71\times 10^{-5}{I}^2+1.79\times 10^{-2}I+0.492$$
系数$K$ 是通过 该月天数$NDM$ 和 最大日照时数$N$ 计算得到的：
$$K=($$