SPI计算方法原理

SPI计算公式

由于降水量的分布一般不是正态分布，而是一种偏态分布。所以在进行降水分析和干旱监测、评估中，采用 Gamma 分布概率来描述降水量的变化。

标准化降水指标（简称 SPI）就是在计算出某时段内降水量的 Gamma 分布概率后，再进行正态标准化处理，最终用标准化降水累积频率分布来划分干旱等级。

标准化降水指数（简称 SPI）的计算步骤为：

1) 假设某时段降水量为随机变量 x，则其 Gamma 分布的概率密度函数如式（1）：

                 （1）

其中:

β>0，γ>0 分别为尺度和形状参数，β和γ可用极大似然估计方法求得，如式(2)，式(3):

                       （2）

                                      （3）

其中

                    （4）

式中:

x_i——为降水量资料样本;

\bar{x}—为降水量气候平均值。

确定概率密度函数中的参数后，对于某一年的降水量x_0，可求出随机变量x小于x_0事件的概率为:

 F(x < x_0) = \int_0^{x_0} f(x) dx                （5）

利用数值积分可以计算用式(1)代人式(5)后的事件概率近似估计值。

b）降水量为0时的事件概率由式（6）估计:

F(x = 0) = \frac{m}{n} 

式中：

m——降水量为0的样本数；

n——总样本数。



c）对Γ分布概率进行正态标准化处理，即将式（5）、式（6）求得的概率值代人标准化正态分
布函数，即:

F(x < x_0) = \frac{1}{\sqrt{2\pi}} \int_{0}^{x_0} e^{-z^2/2} dzF(x<x0​)=2π​1​∫0x0​​e−z2/2dz

对式(7)进行近似求解可得:

 Z = S \left( t - \frac{(c_2t + c_1)t + c_0}{(d_3t + d_2)t + d_1t + 1.0} \right) 

其中：t = \sqrt{\ln\left(\frac{1}{F^2}\right)} 

F为式(5)或式(6)求得的概率；并当F>0.5时,F值取1.0-F,S =1；当F≤0.5时，S=-1。

C_0=2.515 517;

C_1=0.802 853;

C_2=0.010 328;

d_1=1.432 788;

d_2=0.189 269;

d_3=0.001 308。
 

STI计算公式

1. 估计Gamma分布参数：

   - 平均温度

   - 参数

   - 估计的形状参数

   - 估计的尺度参数

  

2. 计算累积概率函数：

   - 使用估计的形状参数 和尺度参数 计算Gamma分布的累积概率函数。

3. 标准化处理：

   - 使用标准正态分布的累积概率函数计算z-score。

   - 根据 z-score 计算STI值。

这些公式是计算标准化温度指数（STI）的关键步骤。

标准化温度指数（The Standardized Temperature Index，STI）是一种简单的、易计
算的温度指数。由于温度异常可以假设为正态分布，因此一年中每一天的温度都符合
正态分布。对一年中每一天的温度进行正态分布拟合，计算出对应的概率分布，
再进行正态标准化处理。根据下面列出的累积分布函数G（x）计算STI 如下所示：

G(x) = \frac{1}{\sigma\sqrt{2\pi}} \int_{-\infty}^{x} \exp\left(-\frac{(x - \mu)2}{2\sigma^2}\right) dx

STI = \phi^{-1}(q)

其中，x 是温度时间序列。x 和σ 分别是均值和标准差参数。q 是累积概率，𝜑−1是逆标
准正态分布函数。

高温事件使用STI进行表征,STI是根据SPI进行定义,其计算方法一致