gma 教程：计算相对湿润度指数（RMI）

安装 gma：pip install gma （依赖的 gdal 需自行安装）

本文基于：gma 2.0.8，Python 3.10

本文用到数据请从 gma 网站获取：https://gma.luosgeo.com/UserGuide/climet/Index/SPI.html 。

公式原理请参考 <参考文献/文档> ，本文不做重复解释。

RMI 函数简介

gma.climet.Index.RMI(PRE, PET, Axis = None, Scale = 1)

---

功能：【相对湿润度指数】。计算相对湿润度指数。

参数：

 PRE: array。降水量（mm）。

 PET: array。潜在蒸散量（mm）。

可选参数：

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

  Scale = int。时间尺度。默认为 1。

返回：array。

参考：

 https://image.data.cma.cn/staticDry/article/20170927151824947.pdf

---

示例：

from gma import climet, io
from gma import io

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

PRE = Data['PRE'].values
PET = Data['ET0'].values

# 分别计算1个月、3个月、6个月、12个月、24个月、60个月尺度的 RMI 数据
RMI1 = climet.Index.RMI(PRE, PET)
RMI3 = climet.Index.RMI(PRE, PET, Scale = 3)
RMI6 = climet.Index.RMI(PRE, PET, Scale = 6)
RMI12 = climet.Index.RMI(PRE, PET, Scale = 12)
RMI24 = climet.Index.RMI(PRE, PET, Scale = 24)
RMI60 = climet.Index.RMI(PRE, PET, Scale = 60)

不同尺度 RMI 结果