SPI

gma.climet.Index.SPI(PRE, Axis = None, Scale = 1, Periodicity = 12, Distribution = 'Gamma', FitMethod = 'MLE', Calibration = None) 1.0.10 +

---

功能：【标准化降水指数】。计算标准化降水指数（Standardized Precipitation Index）。

参数：

 PRE: array。降水量（mm）。

可选参数：

  Axis = int。计算轴。如果不设置（None），多维数据会将所有数据展开到一维计算。

  Scale = int。时间尺度。默认为 1。例如：1月、3月或其他。

  Periodicity = int。周期。默认为 12。

关于周期

周期定义了参与拟合数据的分组方式，例如：

输入月数据，Periodicity = 12 时，每个月份（共 12 组数据）之间互不干扰，独立拟合，更适合针对月份独立分析；

输入月数据，Periodicity = 1 时，所有月份（共 1 组）同时参与拟合，相互影响，更适合分析月份之间的差异；

  Distribution = str1.1.1 +。用于内部拟合/变换计算的分布类型。默认为 Gamma。

支持的分布类型

'Gamma'：伽马分布；

'LogLogistic'：对数逻辑斯蒂分布；

'Pearson3'：皮尔逊 III 分布。

  FitMethod = str 2.0.4 +。用于内部拟合/变换计算的参数估计方法。默认为 MLE。

支持的分布类型

'MLE'：最大似然估计；

'LMoment'：L-矩估计（PWD，概率加权矩）；

'LMoment2'：L-矩估计。

  Calibration = list||slice||None2.0.4 +。参与内部参数拟合运算数据的周期。默认（None）为全部数据。

例如

全部数据：Calibration = None

第1~10个周期：Calibration = [0, 10]

在第1~50个周期中，每隔 2 周期选一组数据：Calibration = slice(0, 50, 2)

返回：array。

参考文献：

 McKee T B, Doesken N J, Kleis J. The relationship of drought frequency and duration to time scales[R]. Eighth Conf. on Applied Climatology, Anaheim, CA: American Meteor Society, 1993.

---

示例 ：

from gma import climet

基于 Excel 表数据（下载 示例数据）

from gma import io

ELSXLayer = io.ReadVector('PRE_ET0.xlsx')
Data = ELSXLayer.ToDataFrame()

PRE = Data['PRE'].values

# 分别计算1个月、3个月、6个月、12个月、24个月、60个月尺度的 SPI 指数
SPI1 = climet.Index.SPI(PRE)
SPI3 = climet.Index.SPI(PRE, Scale = 3)
SPI6 = climet.Index.SPI(PRE, Scale = 6)
SPI12 = climet.Index.SPI(PRE, Scale = 12)
SPI24 = climet.Index.SPI(PRE, Scale = 24)
SPI60 = climet.Index.SPI(PRE, Scale = 60)

不同尺度 SPI 结果

基于栅格数据（下载 示例数据）

import numpy as np
from gma import io
# 读取数据集
PRESet = io.ReadRaster('PRE_Luoyang_1981-2020.tif')
PRE = PRESet.ToArray()
PRE[PRE == PRESet.NoData] = np.nan
# 读取的数据为三维数据（波段，行，列），第一维为时间序列（月数据）。因此按照轴 0 来计算
SPI1 = climet.Index.SPI(PRE, Axis = 0)
SPI3 = climet.Index.SPI(PRE, Axis = 0, Scale = 3)
SPI6 = climet.Index.SPI(PRE, Axis = 0, Scale = 6)
SPI12 = climet.Index.SPI(PRE, Axis = 0, Scale = 12)
SPI24 = climet.Index.SPI(PRE, Axis = 0, Scale = 24)
SPI60 = climet.Index.SPI(PRE, Axis = 0, Scale = 60)
# 存储计算结果
S = [1,3,6,12,24,60]
for i in S:
	# 保存所有结果中的非全 nan 波段。即：去除时间尺度累积时序列前无效的波段。
    io.SaveArrayAsRaster(eval(f'SPI{i}')[i-1:],
                         fr'.\1981-2020_SPI{i}.tif', 
                         Projection = PRESet.Projection,
                         Transform = PRESet.GeoTransform, 
                         DataType = 'Float32', 
                         NoData = np.nan)  

最后一个月（2020年12月）计算结果