这里做的是RH（土壤相对湿度）和PA、SPI、MI、MCI之间的关系： 结果： 看下R方，0.221，说明预测变量PA、SPI、MI、MCI能够解释因变量RH 22.1%，证明RH的波动只有22.1%是由于PA、SPI、MI、MCI造成的，一般统计学要求30%以上是可以接受的。实际数据达到10%就算 Read More
几何校正我后来选择了两个直接运行能出结果的软件，为啥不用envi来直接做，是因为太复杂了，之前也尝试了很多方法，但是要输入数据范围啥的，最终直接选择了用软件。 这里借助了2个工具能达到2种几何校正的结果。 第一个软件： 用cmd进行：先cd到数据和软件在的文件夹内，然后运行命令。 命令格式: FY_ Read More
想去除month列里的1，2，3，4，10，11，12月对应的行 留下5，6，7，8，9月 #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: deletestaion.py @time: 2023/09/22 @de Read More
isin 函数的使用 返回的结果是根据从isin函数传入的列表(list)，筛选出与列表中包含的数值或字符串相同的数据记录, 用法有点类似sql中的"in" 原始数据一堆站点，想筛选目标站点 #!usr/bin/env python # -*- coding:utf-8 -*- """ @autho Read More
目前有10℃以上的降水距平百分率 按照表格这个等级进行划分 原表： 代码如下： #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: jisuan.py @time: 2023/08/25 @desc: """ impo Read More
前面各指标的算法在我的博客中已经有介绍，下面开始计算各指标的求和情况 原始数据： #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: calculateMCI.py @time: 2023/08/24 @desc: " Read More
=VLOOKUP(B2&E2&F2,CHOOSE({1,2},$Q$2:$Q$7027&$T$2:$T$7027&$U$2:$U$7027,$V$2:$V$7027),2,FALSE) 此公式是数组公式，需要以三键输入：CTRL + SHIFT + ENTER。 这么做是为了方便后面算MI。 上面的 Read More
目标是想计算相对湿润度指数MI 得先计算PET： PET的方法可以采用Thornthwaite方法和FAO Penman-Monteith方法，但是彭曼公式涉及到辐射，需要辐射数据比较麻烦，所以采用Thornthwaite方法 原始数据： 处理后的数据： 代码如下： #!usr/bin/env py Read More
目标计算的是SPIW60为近60天的标准化权重降水指数 1、先计算好WAP，上一个帖子已经算好 2、再计算SPI2(WAP) 代码如下： #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: getWAP.py @tim Read More
之前计算过单站点的WAP，现在计算多站点的WAP 数据样子直接参考上几次的文章 直接上代码： #!usr/bin/env python # -*- coding:utf-8 -*- """ @author: Su @file: getWAP.py @time: 2023/07/27 @desc: " Read More