CN114820260B 基于耦合协调度的区域风水复合侵蚀风险评价方法 （西安理工大学）

段.基于耦合协调度的区域风水复合侵蚀风险本发明公开了一种基于耦合协调度的区域风水复合侵蚀风险评价方法，具体包括如下步骤：首先确定研究区域，收集研究区边界shp文NDVI三维耦合协调度的理论累积概率反映风水即可识别出研究区风水复合侵蚀风险变化趋势2步骤1，确定研究区域，收集研究区域边界shp文件步骤2，利用步骤1收集到的研究区气象站的逐日气象蚀气候侵蚀力，并通过插值和掩膜提取计算获得研究区域的逐月风蚀气候侵蚀力栅格数0m_2_1_2_1a_1_1_2_1S0_2_1S和RS0的计算公式如下：3a_2_1；a的计算公式如下：的风蚀气候侵蚀力栅格数据，使用黄河流域边界shp文件在ArcGIS软件中进行裁剪得到逐步骤3，利用步骤1收集到的研究区气象站4到的研究区逐月风蚀气候侵蚀力栅格数据、步骤3计算得到的研究区逐月降雨侵蚀力栅格平均值即可获得研究区若干年平均风水复合侵蚀风2.根据权利要求1所述的基于耦合协调度的区域风水复合侵蚀风险评价方法，其特征_1_1雨侵蚀力栅格数据，使用黄河流域边界shp文件在ArcGIS软件中进行裁剪得到逐月黄河流3.根据权利要求2所述的基于耦合协调度的区域风水复合侵蚀风险评价方法，其特征5ij为子系统第i个指标的第j个年份的归一化值；xij为子系统第i个指标的第j个计算NDVI子系统、风蚀气候侵蚀力子系统和降雨侵i为第i个子系统对总系统有序度的贡献；uij为第i个ij为第i个子系统中第j个指标的权重；64.根据权利要求3所述的基于耦合协调度的区域风水复合侵蚀风险评价方法，其特征（3675.根据权利要求4所述的基于耦合协调度的区域风水复合侵蚀风险评价方法，其特征利用Mann_Kendall趋势检验方法对步骤5计算得到的研究区每个栅格点耦合协调度逐年最优理论累积概率进行趋势计算，Mann_8[0002]土壤侵蚀是一种在全球广泛发生的自然灾害，严重影响人类社会的生产和生[0005]本发明所采用的技术方案是，基于耦合协调度的区域风水复合侵蚀风险评价方算得到的研究区逐月风蚀气候侵蚀力栅格数据、步骤3计算得到的研究区逐月降雨侵蚀力9均值即可获得研究区若干年平均风水复合侵蚀_1_2_1T为平均气温，a_1e_1_2_1S0_2_1a_2_1；(15)每月的风蚀气候侵蚀力栅格数据，使用黄河流域边界shp文件在ArcGIS软件中进行裁剪得_1的降雨侵蚀力栅格数据，使用黄河流域边界shp文件在ArcGIS软件中进行裁剪得到逐月黄雨侵蚀力子系统由每年12个月的月降雨侵蚀力指标构成，上述每个指标需要进行归一化，r阶线性矩估计量lr如下：[0096]步骤5.4，分别计算经验累积概率P经验和每种备用理论分布函数的理论累计概率[0100]利用Mann_Kendall趋势检验方法对步骤5计算得到的研究区每个栅格点耦合协调水平上是显著的，即可得到黄河流域风水复合侵蚀风险的变化趋势空间分布，Mann_[0105]本发明的有益效果是，本发明基于耦合协调度的区域风水复合侵蚀风险评价方[0106]图1是本发明基于耦合协调度的区域风水复合侵蚀风险评价方法中实现技术流程[0107]图2是本发明基于耦合协调度的区域风水复合侵蚀风险评价方法中黄河流域若干[0108]图3是本发明基于耦合协调度的区域风水复合侵蚀风险评价方法中黄河流域若干[0109]图4是本发明基于耦合协调度的区域风水复合侵蚀风险评价方法中黄河流域若干[0110]图5是本发明基于耦合协调度的区域风水复合侵蚀风险评价方法中黄河流域风险[0115]步骤2.1，首先根据步骤1收集到的研究区逐日气象数据，利用FAOPenman_Monteith方法计算每个气象站的逐月潜在蒸散发，FAOPenman_Monteith方法的具体计算_1_2_1T为平均气温，a_1e_1_2_1S0_2_1a_2_1；蚀气候侵蚀力栅格数据，使用研究区边界shp文件在ArcGIS软件中进行裁剪得到逐月研究_1到每月的降雨侵蚀力栅格数据，使用研究区边界shp文件在ArcGIS软件中进行裁剪得到逐算得到的研究区逐月风蚀气候侵蚀力栅格数据、步骤3计算得到的研究区逐月降雨侵蚀力[0158]步骤4的具体步骤为：步骤4.1[0163]式中，Gij表示指标数据xij对第j项指标的总贡献度；m为样本数(研究期年份长[0170]式中ui为第i个子系统对总系统有序度的贡献；uij为第i个子系统中第j个指标的T＝au1+bu2+cu3的r阶线性矩估计量lr如下：[0183]步骤5.4，分别计算经验累积概率P经验和每种备用理论分布函数的理论累计概率FAOPenman_Monteith方法的具体计算_1_2_1T为平均气温，a_1e_1_2_1S0_2_1一天中最高温度，℃；Tmin,K为一天中最低温度，℃；其中，RS和RS0的计算公式如下：a_1a的计算公式如下：;在ArcGIS软件中进行裁剪得到逐月黄河流域风蚀气候_1到每月的降雨侵蚀力栅格数据，使用黄河流域边界shp文件在ArcGIS软件中进行裁剪得到ij为子系统第i个指标的第j个年份的归ij为子系统第i个指标的第j个年份的值；由于风蚀气候侵蚀力和降雨侵蚀力对侵常常认为三个子系统的重要性相同，所以a＝b＝c＝1/3。黄河流域1982_2015年NDVI子系θ)；[0226]步骤5.4，分别计算经验累积概率P经验和每种备用理论分布函数的理论累