Python的gma包计算标准化降水指数SPI
采用干旱气象指标的气象干旱综合指数MCI时因公式要涉及标准化降水指数SPI,故需要先进行SPI的计算:
SPI的计算依赖python里的包gma,安装方法参考:
https://blog.csdn.net/weixin_42155937/article/details/121023363
原始数据是每个月的降水数据:
python计算SPI代码如下:
#!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: calculateSPI.py @time: 2023/05/16 @desc: """ import gma import pandas as pd # 读取数据 Data = pd.read_excel('processdata/datas/rainmonth.xlsx',sheet_name='50525') # 看一下 Excel 文件前 5 行内容 a = Data.head(6) print(a) # 提取 rain 用于 SPI 运算 rain = Data['rain'].values # 分别计算1个月、3个月、5个月尺度的 SPI 数据 SPI1 = gma.climet.SPI(rain) SPI3 = gma.climet.SPI(rain, Scale = 3) SPI5 = gma.climet.SPI(rain, Scale = 5) # 将结果保存到文件 OUT = pd.DataFrame([SPI1, SPI3, SPI5], index = ['SPI1','SPI3','SPI5']).T OUT.to_excel(r'processdata/datas/SPI.xlsx', index = False)
结果:
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 【.NET】调用本地 Deepseek 模型
· CSnakes vs Python.NET:高效嵌入与灵活互通的跨语言方案对比
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· Plotly.NET 一个为 .NET 打造的强大开源交互式图表库