- 博客(14)
- 收藏
- 关注
RSS订阅原创 Python：SPEI指数的计算
将nc文件修改为想要格式：参考sample_precip.nc和sample_tavg.nc(1）根据ERA5小时数据求月总降水和月平均气温在python中(2) 在nco中修改nc文件的格式ncatted -a units,precip,c,c,millimeter in.nc #给precip变量添加单位属性为“millimeter”ncatted -a units,tavg,c,c,K in.nc #给tavg变量添加单位属性为“K”
2023-08-09 15:12:05
1649
1
原创 WINDOWS系统下NCO安装与使用
ncatted -a units,precip,c,c,millimeter in.nc #给precip变量添加单位属性为“millimeter”ncatted -O units,precip,o,c,"millimeter" in.nc #修改precip的单位为millimeter。ncatted -a units,tavg,c,c,K in.nc #给tavg变量添加单位属性为“K”ncatted -a ,precip,d,, in,nc #删除precip变量的所有属性。
2023-08-09 14:57:13
934
1
原创 nc时间object格式转化为datetime64格式
data['time'] = data.indexes['time'].to_datetimeindex() ## 在原xarray数据中修改时间变量格式。time = data.indexes['time'].to_datetimeindex() # 转化格式--->只有时间变量。time_new = time+datetime.timedelta(hours=8) #得到往后推移了8h后的时间。data = xr.open_dataset(r"") # 输入文件。###后面可以做时间匹配。
2023-08-08 17:20:39
500
1
原创 Linux系统cdo对nc数据重采样
在终端输入sudo apt-get install cdo。安装成功后输入cdo -V 来检验是否成功安装以及安装版本。输入cdo griddes input.nc > grid。主要修改x间隔、y间隔、x初始、y初始、x格网数、y格网数。在终端中输入sudo apt install vim。输入 cdo griddes .nc。2.下载vim命令，可以编辑文件。grid为修改过后的网格文件。3.查看原始nc文件的网格。6.利用新的网格进行重采样。4.导出原始数据的网格。
2023-08-07 16:26:18
783
2
原创 基于灾损曲线的综合风险与指标体系综合风险验证
（2）计算空间上灾害的重现期，重现期越大，代表灾害程度越严重，发生可能性越小。计算过程：得到灾害各危险性要素37年间的内蒙古全域平均值，构建单变量边际分布，选择最优copula，三种灾害均为Clayton Copula函数，计算联合重现期，计算联合概率。（4）加载到arcgis中，在工具箱中打开Band Collection Statistics工具，计算空间相关系数。（3）得到空间上的联合重现期后，在python中根据灾损曲线分别计算经济损失（风险）、人口损失（风险）。1. 计算基于灾损曲线的综合风险。
2023-08-05 14:15:37
195
1
原创 不同重现期下单/多灾害综合风险计算
加权综合：在ArcGIS栅格计算器输入“Float("15年重现期下危险性_CON") * 0.52 + Float("暴雨-层次分析法_敏感性1_CON") * 0.19 + Float("暴雨-层次分析法_脆弱性1_CON") * 0.19 +(1- Float("暴雨-层次分析法_防灾减灾能力_CON")) * 0.10”；工具箱-栅格计算器，输入公式：（栅格最大值-栅格）/（栅格最大值-栅格最小值）在栅格计算器中根据权重加权综合三种灾害风险得到多灾害的综合风险图。计算步骤：（以暴雨灾害为例）
2023-08-01 11:14:19
121
1
原创 路网密度计算
1.利用【数据管理——投影与转换工具——矢量——投影】给只有地理坐标系的两个图层定义投影，注意：地理坐标系转换为投影坐标系用【投影】工具，而不是【定义投影】工具，投影类型为WGS1984 UTM-49N。2.创建渔网：【数据管理——sampling——create fishnet】，定义分辨率为0.099°（0.1°）。4.将路网按渔网进行分割：【Analysis Tools——Overlay——Identity】5.按照fish的id进行字段合并：【geoprocessing——dissolve】
2023-06-20 10:09:03
1522
原创 ARCGIS批量标准化栅格数据（最大最小值标准化）
env.workspace = "E:/历史时期综合风险评估/暴雨洪涝灾害风险评估/敏感性数据/"out_path = "E:/历史时期综合风险评估/暴雨洪涝灾害风险评估/敏感性标准化数据/"#设置读取的栅格类型，可以参考ArcGIS的官方文档，下面的是读取目录下所有类型的栅格文件。... print("最大值:" + str(maxValue))... print("最小值:" + str(minValue))#这个是工作空间目录，待处理的所有影像都在这个文件型数据库中。
2023-06-15 19:30:22
1107
原创 matlab 趋势面分析（二）
三次：coef3=[4953.2 -8389.4 -47.091 -3171 -159.15 33308 -304.95 10888 -37105 -11688]一次：coef1 = [-4.237070308572415e+02;一次：coef1 = [-1.329694198065215e+02;二次：coef2 = [-2528.7 -740.27 -2.5713 -39.158 2732.9 16.35]
2023-05-31 22:39:27
281
1
原创 matlab 趋势面分析 记录
num = xlsread("G:\Data\灾损曲线\趋势面分析-数据.xlsx",'B21:N37');ylabel('暴雨强度（mm/d）');
2023-05-31 17:39:36
1113
1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人