ET0处理杂记
发布网友
发布时间：2024-10-01 15:10
共1个回答
热心网友
时间：2024-10-17 19:13
吐哈地区
精度评价：气象站点数据验证
将吐哈地区1955-2020年逐月均值导出一幅图
【目标】依据吐哈地区11个气象站点shp文件，提取1955-2020年逐月/年的蒸散发栅格数据上对应的值
【解决】主要参考前两个博文，然后根据报错内容进行略微的修改，第三个因为没能提取出shp失败了。
将下面代码复制粘贴过去，根据自己数据的路径修改进行修改即可
【试错】前面试了两次，因为是直接粘过去，第四行的arcpy竟然变成了 a rcpy，还有就是需要将空行清除及路径的设置错误，导致结果没跑出来
打开生成的txt文件，将结果粘贴到excel；选中刚才粘进来的那一列数据，点击数据→分列→分列→下一步→分割符号勾选逗号→完成；将shp文件里的站点信息粘进来，大功告成，结果如下：
其余数据同理
与ET0产品对比
利用FileZilla3.60.1从国家青藏高原科学数据中心上下载1990-2020年中国1km逐月潜在蒸散发数据集
参考 R语言学习记录：nc格式转tif - 哔哩哔哩 (bilibili.com)，这里用的法三，将下载的nc文件转成tif
批量转tif结果如下，1990-2020年，共31年，372个月
接下来，利用QGIS与研究区矢量文件批量裁剪，转为真值
首先加载研究区矢量文件
菜单栏→栅格→提取→按掩膜图层裁剪栅格→以批处理方式运行
输入图层→从一个目录中添加全部文件（根据自己的情况更改）；目标CRS→工程CRS EPSG：4326→向下填充；无数据值设为0；保留输入栅格的分辨率；设置输出路径；运行
转换为真值，批量*0.1
菜单栏→栅格处理→工具箱→GDAL→栅格杂项→栅格计算器→以批处理方式运行
处理MODIS16A2潜在蒸散发
参考PML下载
批量修改文件名
批量合成年ET0批量获取均值
ArcPy批量获取并输出一个文件内所有文件的均值，主要参考以下博文
MAE平均绝对误差+RMSE均方根误差
Hargreaves-Samani计算1955-2020年吐哈地区逐月ET0
月尺度合成年尺度 Monthly scale combined with annual scale
根据气象学上对季节的划分，3–5月为春季，6–8月为夏季，9–11月为秋季，12月至翌年2月为冬季。
ET0计算方法大致分为温度法、辐射法、综合法和蒸发皿法，如公认的Penman-Monteith[4,8,9]、最佳候选Thornthwaite[10]、应用广泛的Hargreaves-Samani[11,12]、人工神经网络ANN[13]等。
吐哈地区ET0在年月尺度上的变化与均温相一致(Fig 4)。年尺度上，吐哈地区潜在蒸散量介于2519.92~2735.93 mm/year，多年平均ET0为2597.77mm,变化幅度较大，于1984年达到最低点，在整个时段内呈现缓慢上升的趋势。
月尺度上，吐哈地区ET0介于24.18-448.14mm，呈单峰型，峰值范围位于5-8月，最低值出现在最冷月，这与水热分布时间相似。其中7月温度最高，但因本月降水增多，空气中的相对湿度增大，ET0也会随之减少。
处理1960年-2020年全国1km逐月降水数据集
【来源】 数据分享 | 中国1公里分辨率月降水数据集（1960-2020年）
【下载】 pan.baidu.com/s/1x5WayN...
【趋势分析】 python进行栅格一元回归分析（slope），一学就会！（不用改代码参数，只换文件地址）
热力图也可以玩这种花样，学到了
1.学硕： 全国学术型硕士研究生学生资助调查问卷
2.专硕： 全国专业型硕士研究生学生资助调查问卷
作物蒸散发-计算作物需水量指导手册（符号与缩写列表）
Remote Sensing专刊 | 作物系数和蒸散发的遥感估算
