目标计算的是SPIW60为近60天的标准化权重降水指数 1、先计算好WAP,上一个帖子已经算好 2、再计算SPI2(WAP) 代码如下: #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: getWAP.py @tim Read More
posted @ 2023-07-27 16:37 秋刀鱼CCC Views(336) Comments(0) Diggs(0) Edit
之前计算过单站点的WAP,现在计算多站点的WAP 数据样子直接参考上几次的文章 直接上代码: #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: getWAP.py @time: 2023/07/27 @desc: " Read More
posted @ 2023-07-27 16:11 秋刀鱼CCC Views(222) Comments(0) Diggs(0) Edit
前面提过如何计算单站点多年的标准化降水指数SPI,但是在现实中不仅只计算一个站点的标准化降水指数,有时要计算多个站点的。 原始数据是11个站点1961年至2022年每个月的降水值: 编程思路是: 1、先需要进行站点分类: 此处的难点在于需要写个for循环,根据站点名称遍历数据 for i in na Read More
posted @ 2023-07-26 16:00 秋刀鱼CCC Views(649) Comments(2) Diggs(0) Edit
原excel: 目标将year和month合并: year、month里放的1961等是数字,合并日期的时候需要把它们变成字符串再合并,采用.astype(str) #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: Read More
posted @ 2023-07-19 09:13 秋刀鱼CCC Views(275) Comments(0) Diggs(0) Edit
上一节计算了每月的平均降水量,现在要计算降水距平,算法是把每年每个月的降水量减去每月的平均降水量 数据分为两个excel 方法是把name和month作为key,pd.merge,再新增一列,两列做差 这两张表用两个dataframe存,把原excel的两个rain区分一下,变成rainA和rain Read More
posted @ 2023-07-03 16:52 秋刀鱼CCC Views(275) Comments(0) Diggs(0) Edit
#!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: calculatejsjp.py @time: 2023/06/09 @desc: """ import pandas as pd # 读取数据 Data = pd. Read More
posted @ 2023-06-27 17:07 秋刀鱼CCC Views(152) Comments(0) Diggs(0) Edit
原始数据跟上一篇一样 计算公式如图: 公式解析: 实则就是当前的降水量一列,前一个月的降水量a的二次方×p2一列,这两列错位相加。 #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: calculateWAP.py @ Read More
posted @ 2023-05-17 16:44 秋刀鱼CCC Views(279) Comments(0) Diggs(0) Edit
采用干旱气象指标的气象干旱综合指数MCI时因公式要涉及标准化降水指数SPI,故需要先进行SPI的计算: SPI的计算依赖python里的包gma,安装方法参考: https://blog.csdn.net/weixin_42155937/article/details/121023363 原始数据是 Read More
posted @ 2023-05-16 11:19 秋刀鱼CCC Views(1683) Comments(6) Diggs(0) Edit
解决每一个遇到的异常,建议大家遵循如下三点: 1. 不慌张,细看信息,定位错误。看清楚报的错误信息,并定位发生错误的地方。 2. 百度并查看十个相关帖子。将异常类信息进行百度,至少查看十个以上的相关帖子。 3. 以上两步仍然无法解决,找老师和同学协助解决 Read More
posted @ 2023-05-15 16:03 秋刀鱼CCC Views(49) Comments(0) Diggs(0) Edit
python版本为3.7的安装不了gma包,得安装一个python3.8的虚拟环境,再尝试安装gma包。 Read More
posted @ 2023-05-15 16:02 秋刀鱼CCC Views(102) Comments(0) Diggs(0) Edit