(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211470082.7 (22)申请日 2022.11.23 (71)申请人 水利部交通 运输部国家能源局南京 水利科学研究院 地址 210029 江苏省南京市 鼓楼区广州路 223号 (72)发明人 雷少华　尤俊坚　金秋　耿韧　 徐春　周璐　沈珊珊　 (74)专利代理 机构 北京弘权知识产权代理有限 公司 11363 专利代理师 逯长明　温瑞鑫 (51)Int.Cl. G06F 17/11(2006.01) G06F 30/20(2020.01) G06F 119/14(2020.01) (54)发明名称 一种土壤侵蚀强度估算方法、 系统及 存储介 质 (57)摘要 本申请涉及土壤侵蚀分析技术领域， 具体而 言， 涉及一种土壤侵蚀强度估算方法、 系统及存 储介质， 一定程度上可以解决计算获得的土壤侵 蚀强度与实际情况存在较大的误差的问题。 所述 土壤侵蚀强度估算方法包括： 获取年降雨侵蚀 力， 并根据卫星降雨量数据对年降雨侵蚀力进行 拟合优化， 输出年降雨侵蚀力的拟合优化计算公 式， 其中， 所述年降雨侵蚀力是基于气象站逐日 降雨数据计算得到的； 使用年降雨侵蚀力的拟合 优化计算公式对年降雨量卫星栅格数据进行采 样， 获取并输出优化后的年降雨侵蚀力； 获取土 地利用类型， 并根据土壤流失方程， 将优化后的 年降雨侵蚀力带入至土壤流失方程， 获取并输出 土壤侵蚀强度。 权利要求书2页 说明书7页 附图4页 CN 115510377 A 2022.12.23 CN 115510377 A 1.一种土壤侵蚀强度估算方法， 其特 征在于， 所述方法包括： 获取年降雨侵蚀力， 并根据卫星降雨量数据对年降雨侵蚀力进行拟合优化， 输出年降 雨侵蚀力的拟合优化计算公式， 其中， 所述年降雨侵蚀力是基于气象站逐日降雨数据计算 得到的； 使用年降雨侵蚀力的拟合优化计算公式对年降雨量卫星栅格数据进行采样， 获取并输 出优化后的年降雨侵蚀力； 获取土地利用类型， 并根据土壤流失方程， 将优化后的年降雨侵蚀力带入至土壤流失 方程， 获取并输出土壤侵蚀强度。 2.如权利要求1所述土壤侵蚀强度估算方法， 其特征在于， 使用年降雨侵蚀力的拟合优 化计算公式对年降雨量卫星栅格数据进行采样， 获取并输出优化后的年降雨侵 蚀力的过程 中， 所述方法包括： 获取空间分辨率为1  km的逐月降雨量栅格数据， 所述空间分辨率1  km的逐月降雨量栅 格数据为 根据全球气候数据集 通过Delta空间降尺度生成； 将空间分辨率为1  km的逐月降雨量栅格数据按照年度进行累加获得年降雨量格栅数 据； 基于年降雨侵蚀力的拟合优化计算公式， 对年降雨量格栅数据进行采样并生成空间分 辨率为2 m的年降雨侵蚀力。 3.如权利要求1所述土壤侵蚀强度估算方法， 其特征在于， 获取年降雨侵蚀力的过程 中， 使用如下的公式： 其中， 为年降雨侵蚀力； 为第i个半月的降雨侵蚀力， i取1~24； 为该年 第y个侵蚀性日雨量； 为回归系数， 当时间为冷季 （10~12月， 1~4月） 时， =0.3101， 当降 雨时间为暖季 （5~9月） ， =0.3937； y为该年第i个半月内降雨量大于等于10  mm的侵蚀性降 雨总日数； m为半月内侵蚀性降雨日数； 为第i个半月平均降雨侵蚀力 占年降 雨侵蚀力 的权重。 4.如权利要求1所述土壤侵蚀强度估算方法， 其特征在于， 根据卫星降雨量数据对年降 雨侵蚀力进行拟合优化的过程中， 所述方法包括： 获取卫星降雨量数据， 使用卫星降雨量数据对获得的年降雨侵蚀力进行拟合， 输出年 降雨侵蚀力对于卫星降雨 量数据的幂函数； 使用决定系数以及平均绝对百分比误差对年降雨侵蚀力对于卫星降雨量数据的幂函权　利　要　求　书 1/2 页 2 CN 115510377 A 2数进行判断； 输出拟合效果 最优的年降雨侵蚀力对于卫星降雨 量数据的幂函数。 5.如权利要求1所述土壤侵蚀强度估算方法， 其特征在于， 获取土地利用类型的步骤 中， 所述方法包括： 获取第一影像数据和第二影像数据， 并融合第一影像数据和第二影像数据， 输出融合 后的第三影像数据， 其中， 所述第一影像数据为空间分辨率为30  m的多光谱RGB 波段影像数 据， 所述第二影 像数据为空间分辨 率为15m的全色波段影 像数据； 获取第四影像数据， 融合第 三影像数据和第四影像数据， 输出融合后的第五影像数据， 其中， 第四影 像数据为空间分辨 率为2m的多光谱RGB波段影 像数据； 对第五影 像数据进行区域分割， 获取并输出土地利用类型。 6.如权利要求5所述土壤侵蚀强度估算方法， 其特征在于， 融合第 一影像数据和第 二影 像数据， 输出融合后的第三影 像数据的步骤中， 所述方法包括： 使用至少两种影像数据融合方法对第 一影像数据和第 二影像数据进行融合， 并输出各 自对应的第三影 像数据； 通过对两种影像数据融合方法生成的第 三影像数据进行比较， 选取包含信 息量较多的 第三影像数据作为 最优的第三影 像数据输出； 其中， 所述信息量包括： 均值、 标准差相关系数、 光谱扭曲度、 平均梯度以及信息熵。 7.一种土壤侵蚀强度估算系统， 其特 征在于， 所述土壤侵蚀强度估算系统包括： 降雨侵蚀力计算单元， 获取年降雨侵蚀力， 并根据卫星降雨量数据对年降雨侵蚀力进 行拟合优化， 输出年降雨侵 蚀力的拟合优化计算 公式， 其中， 所述年降雨侵 蚀力是基于气象 站逐日降雨数据计算得到的； 降雨侵蚀力优化单元， 用于使用年降雨侵蚀力的拟合优化计算公式对年降雨量卫星栅 格数据进行采样， 获取并输出优化后的年降雨侵蚀力； 土壤侵蚀强度获取单元， 用于获取土地利用类型， 并根据土壤流失方程， 将优化后的年 降雨侵蚀力带入至 土壤流失方程， 获取并输出土壤侵蚀强度。 8.一种存储介质， 其特征在于， 存储有计算机可执行指令， 所述计算机可执行指令被处 理器执行时执行根据权利要求1 ‑6中任一项所述的土壤侵蚀强度估算方法。权　利　要　求　书 2/2 页 3 CN 115510377 A 3