SPEI的python代码

### SPEI Python 实现

为了计算标准化降水蒸散发（SPEI），可以利用 `climate_indices` 庌库，这是一个专门用于气候指数分析的强大工具[^2]。

首先，确保已经安装了必要的库：

pip install spei-python climate-indices

下面是一个完整的代码示例来展示如何使用这些库来进行 SPEI 的计算:

import pandas as pd
from climate_indices import compute, indices
from datetime import date

# 假设有一个CSV文件包含了降水量和潜在蒸发量的数据
data = pd.read_csv('path_to_your_data.csv')

precipitation_series = data['Precipitation'].values  # 提取降水量序列
pet_series = data['PotentialEvapotranspiration'].values  # 提取潜在蒸发量序列
dates = [date(year=year, month=month, day=1) for year in range(1980, 2010) for month in range(1, 13)]  # 创建日期列表

# 计算SPEI 使用时间尺度为6个月
spei_6m = compute.spei(precips_mm=precipitation_series,
                       pet_mm=pet_series,
                       scale=6,
                       periodicity='monthly',
                       dates=dates)

print(spei_6m)

这段代码读入了一个 CSV 文件中的数据，并调用了 `compute.spei()` 函数来执行实际的 SPEI 计算。这里假设输入数据是以毫米为单位记录的月度累积降水量 (`precip`) 和潜在蒸发量 (`pet`) 。此外还指定了一个时间窗口大小 (scale)，在这个例子中设置成了六个月的时间跨度。