- 博客(12)
- 资源 (31)
- 收藏
- 关注
RSS订阅原创 python绘制雷达图
在以前的研究中喜欢做各类因子与其中一个因子的相关系数，放在表格里面。但当时想到能用图片直观表达出来看哪个因子相关更高，于是做了雷达图，例如下面的图片，绘制的是各个因子和城市内涝相关性的绝对值。
2024-10-03 20:13:46
221
原创 python 计算纳什系数，KGE，MAE，RMSE等指标
水文计算的时候常用NSE纳什系数、KGE这两个指标计算模型模拟的准确性，比较喜欢用python来单独对模拟值、观测值进行对比，分享一下简单的代码。‘./test.txt’ 第一列为观测值，第二列为模拟值。
2024-10-03 20:06:27
323
原创 ERA5降雨（Total Precipitation）全球逐小时数据下载
数据量实在太多，ERA5全球数据需要一个月一个月逐个下载，所以也参考了网上的大佬们的代码，改编了一下，这里分享一下每个月下载的代码，因为传输过程随时可能中断，所以最好用个try反复一下，但本人的python 还没精进到那个地步，不会用-- 以后再说吧。不过前提是注册了账号，以及在 vim $HOME/.cdsapirc 目录里面有自己的账号和API，详细过程可以参照这位的过程。今年暴雨也挺多的，下了ERA5所有数据想进行分析极端暴雨，过程中遇到很多问题。保存后才能用批量下载的代码正常下载。
2024-09-29 13:55:06
512
原创 数据处理，读取每一天逐三小时asc文件、逐日asc文件，将数据上下（南北）颠倒过来
由于以前没怎么用C++自带的处理时间的函数，习惯了用自己写的get_nexttime函数，用string 类型的年月日时，获取增加3小时的年月日时，如果时获取一天以后的，则直接循环8次。当然也可以用C++自带的，本人太懒了。最近处理CMIP6数据发生了许多问题，由于python读取nc文件写完了之后存储成了南北颠倒的asc文件，这里需要再将文件的南北倒过来。
2024-09-29 13:33:19
443
原创 读取CMIP6的nc文件，转换成每日一个asc格式文件
由于模型方便读取降雨数据，由之前处理好的标准化CMIP6 nc文件读取每日的降雨，改写成行列标准化的asc文件。
2024-09-28 11:03:36
343
原创 ArgGIS数据存储asc转tif压缩节省空间，linux gdal_translate 命令
由于是Linux环境，所以采用gdal_translate命令，先把规范行列的数据写成AAIGrid格式，也就是asc格式，假如初始文件时 Prec_2010010100.asc，浮点类型，想要得到压缩的存储文件Prec_2010010100.tif， 运行命令。最开始的数据格式采用二进制Bin格式，为了数据规范读取，输出，但是后来发现Bin格式极大浪费了空间。不过在降雨数据存储上，可以将一位或者两位小数点的数据乘以10或者100，因为降雨数据日最大数值范围不会超过1000，可以用Int类型。
2024-09-28 10:24:31
456
原创 ERA5 Land 降雨（Total precipitation）和径流（Runoff）数据
每日的00时是前一天的累积量，比如 2010年5月2日的TP_2010050200 是 2010年5月1日00时至23时的总累计量。所以当需要逐小时数据，需要把从每天的时次减去前一个时次，从01时开始，2010050100不用变，02时则需要将 2010050200-2010050100，03时则 2010050300 - 2010050200，依次类推至下一日的00时，2010050300-2010050223，然后 2010050301不变，循环如此。ERA5 逐小时数据遇到的问题。
2024-09-27 13:01:31
629
4
原创 CMIP6数据处理
根据grid信息双曲线性插值命令，这里假设原nc文件为pr_day_MRI-ESM2-0_ssp585_r1i1p1f1_gn_20780622-20991231.nc， 输出文件名为MRI-ESM2-0_ssp585_20780622-20991231.nc。这样调整好了，用python再读取每个数据的全球空间范围的降雨数值，空间分辨率1X1°，由于单位不是mm/天，把每个数值乘以86400，得到每一天的数值。下载好了，要处理时间分辨率和空间分辨率，推荐用cdo，首先要pip install cdo。
2024-09-27 11:38:07
486
原创 生成有箭头的流向Line.shp 并通过ArcMap绘制带箭头的图，根据D8算法，输入Flow Direction
简单说明：代码C++，配置GDAL环境（必须），ArcMap附一个配置GDAL环境的教程，感谢这位大佬VS2017编译配置GDAL库输入流向数据的 asc文件，附网盘测试链接链接：https://pan.baidu.com/s/1-Q35OktO8dOPbjdkI1XC8w提取码：p2ap代码如下：#include <iostream>#include <fstream>#include <string>#include <string.h>
2022-04-01 19:12:00
1413
原创 C++ 调用 SWMM模型.swmm5.dll
现在啥都要讲究自动化，一年前在深入学习了SWMM模型后经历了一段技术困扰，特别是如何自动运行SWMM模型生成报告并作处理最终变成TIFF等能用GIS展示的结果。那会从自动处理降雨(变成inpfile里面的格式)、输入模型、运行模型、输出结果、调用GDAL库生成shp，TIFF等等…纪念一下如何调用 swmm5.dll,也忘记是从哪里学来的了，在这里感谢网上分享技术、代码的大佬们！//run swmm model.char inpfile[]=“D:\\New.inp”;char out_rpt[]=
2020-08-22 03:30:47
1404
2
黑龙江省DEM（30米分辨率）
2022-04-08
海南省DEM（30米分辨率）
2022-04-03
江苏省DEM（30米分辨率）
2022-04-03
ERA5 再分析数据例子，范围覆盖全球，降雨，逐日数据，浮点格式，GeoTIFF文件类型
2024-09-27
广东Landcover数据(10m分辨率)
2022-04-12
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人