一种土壤侵蚀地形因子计算方法、系统及存储介质技术方案

本发明专利技术公开了一种土壤侵蚀地形因子计算方法、系统及存储介质，包括数据输入，获取并输入计算土壤侵蚀地形因子所需的数据；数据处理，包括填洼方法设置、算法选择、阈值设置及模型计算，首先选择填洼方法，选择计算土壤侵蚀因子的算法，然后设置算法所需的阈值参数，之后通过土壤侵蚀模型计算；结果输出，输出土壤侵蚀地形因子数据和图像。具有精确的计算出土壤侵蚀地形因子，获得土壤侵蚀地形因子的数据和图像的技术效果。和图像的技术效果。和图像的技术效果。

【技术实现步骤摘要】
一种土壤侵蚀地形因子计算方法、系统及存储介质


[0001]本专利技术涉及软件算法模型和土壤侵蚀地形因子计算领域，具体而言，涉及一种土壤侵蚀地形因子计算方法、系统及存储介质。

技术介绍

[0002]目前，地形因子具体分为坡长因子和坡度因子，坡长因子是指某一坡面土壤流失量与坡长为22.13m、其它条件一致的坡面产生土壤流失量的比值。坡度因子是指某一坡度土壤流失量与坡度为5.13
°
，其它条件一致的坡面产生土壤流失量的比值。地形因子计算的基础数据为1:10000电子等高线图。
[0003]综上所述，存在如下技术问题：
[0004]如何快速精确的计算出土壤侵蚀地形因子，获得土壤侵蚀地形因子的数据和图像。

技术实现思路

[0005]本专利技术的主要目的在于提供一种土壤侵蚀地形因子计算方法、系统及存储介质，以解决
技术介绍
中至少一个技术问题。
[0006]为了实现上述目的，根据本专利技术的一个方面，提供了一种土壤侵蚀地形因子计算方法，包括：
[0007]数据输入，获取并输入计算土壤侵蚀地形因子所需的数据；
[0008]数据处理，包括填洼方法设置、算法选择、阈值设置及模型计算，首先选择填洼方法，选择计算土壤侵蚀因子的算法，然后设置算法所需的阈值参数，之后通过土壤侵蚀模型计算；
[0009]结果输出，输出土壤侵蚀地形因子数据和图像。
[0010]优选的，所述填洼方法的算法包括MV算法和快速算法:选择MV算法或快速算法，以及设置限制因子，其中MV算法为将洼地底点栅格高程升高至其流向出口的高程值，从而使水流能够顺畅地流向出口的方法；快速算法为利用优先队列快速填洼的方法；限制因子即填洼的高程大小与当前数据中的最高点的乘积。
[0011]优选的，所述MV算法和快速算法中的设置包括：坡度计算设置、流向设置、流量计算设置、坡向计算设置和坡长因子计算方法设置，其中坡长因子计算方法包括坡长法和汇流面积法，其中坡向计算设置包括简单算法、南北走向算法和东西走向算法，其中流向计算设置包括流向汇集和快速汇集，其中流向设置中包括D6算法和快速算法，所述坡度计算设置中包括三阶差分。
[0012]优选的，所述阈值设置中设置的参数包括：坡度、坡长、回流量、去短枝、缓坡中断因子和陡坡中断因子。
[0013]优选的，填洼方法采用扫描窗口法进行，先找到每个平地栅格或是洼地底点栅格，即潜在出流点，然后找到最低的潜在出流点，比较其和洼地栅格的高程，如果出流点高程高
于洼地栅格，那么洼地是一个凹地，否则是一个平地。
[0014]优选的，所述坡度计算和坡向计算采用拟合曲面法对填洼处理后的DEM模型进行坡度和坡向计算，拟合曲面采用二次曲面。
[0015]优选的，采用最大坡降算法计算水流流向，即栅格坡度的最佳代表值是以之为中心3
×
3窗口内其周围八个方向坡度最大值，水流方向与最大坡降一致，流向模型采用D8流向模型。
[0016]优选的，使用计算出的每个单元格的流向生成流向网络，流向网络是一种以单元格为节点，以水流方向为边的有向图，每个单元格都可以通过流向网络找到其紧邻的下游单元格，从下游到上游逐一遍历流向网络，累加每个单元格的面积，计算出每个单元格的汇流累积量，汇流累积量表示流向该单元格的所有水流的总和，其中最下游的单元格的汇流累积量为0。
[0017]根据本专利技术的另一方面提供了一种土壤侵蚀地形因子计算系统，包括：
[0018]数据输入模块，用于获取并输入计算土壤侵蚀地形因子所需的数据；
[0019]数据处理模块，用于设置算法和设置阈值及模型计算，用于选择或自定义计算土壤侵蚀因子的算法，用于设置算法所需的阈值参数，用于土壤侵蚀模型计算；
[0020]结果输出模块，用于输出土壤侵蚀地形因子数据和图像。
[0021]根据本专利技术的另一方面提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现土壤侵蚀地形因子计算方法的步骤。
[0022]应用本专利技术的技术方案，具有如下技术效果：
[0023]通过获取并输入计算土壤侵蚀地形因子所需的数据；然后数据处理，包括填洼方法设置、算法选择、阈值设置及模型计算，首先选择填洼方法，选择计算土壤侵蚀因子的算法，然后设置算法所需的阈值参数，之后通过土壤侵蚀模型计算；最后结果输出，输出土壤侵蚀地形因子数据和图像。具有精确的计算出土壤侵蚀地形因子，获得土壤侵蚀地形因子的数据和图像的技术效果。
附图说明
[0024]构成本申请的一部分的说明书附图用来提供对本专利技术的进一步理解，本专利技术的示意性实施例及其说明用于解释本专利技术，并不构成对本专利技术的不当限定。在附图中：
[0025]图1示出了根据本专利技术的土壤侵蚀地形因子计算方法的流程示意图；
[0026]图2示出了图1中的土壤侵蚀地形因子计算方法的计算逻辑图；
[0027]图3示出了图1中的土壤侵蚀地形因子计算方法的拟合窗体示意图；
[0028]图4示出了图1中的土壤侵蚀地形因子计算方法的D8流向模型图；
[0029]图5示出了图1中的土壤侵蚀地形因子计算方法的数据表和输出图像图；
[0030]图6示出了图1中的土壤侵蚀地形因子计算方法的DEM图；
[0031]图7示出了图1中的土壤侵蚀地形因子计算方法的坡度图像图；
[0032]图8示出了图1中的土壤侵蚀地形因子计算方法的坡长图像图；
[0033]图9示出了图1中的土壤侵蚀地形因子计算方法的坡向图像图；
[0034]图10示出了图1中的土壤侵蚀地形因子计算方法的坡长因子图像图；
[0035]图11示出了图1中的土壤侵蚀地形因子计算方法的流向图像图；
[0036]图12示出了图1中的土壤侵蚀地形因子计算方法的坡度因子图像图；
[0037]图13示出了图1中的土壤侵蚀地形因子计算方法的流量图像图；
[0038]图14示出了图1中的土壤侵蚀地形因子计算方法的数据图。
具体实施方式
[0039]需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本专利技术。
[0040]如图1至图14所示，本专利技术实施例提供了一种土壤侵蚀地形因子计算方法，包括数据输入，获取并输入计算土壤侵蚀地形因子所需的数据；数据处理，包括填洼方法设置、算法选择、阈值设置及模型计算，首先选择填洼方法，选择计算土壤侵蚀因子的算法，然后设置算法所需的阈值参数，之后通过土壤侵蚀模型计算；结果输出，输出土壤侵蚀地形因子数据和图像。可以提出快速生成需要的数据同时提供用户自由成图的操作，生成的数据为.t i f，可以将软件生成的数据直接用于其他软件进行后续操作。
[0041]本实施例中，数据输入，获取并输入计算土壤侵蚀地形因子所需的数据；数据处理，包括填洼方法设置、算法选择、阈值设置及模型计算，首先本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种土壤侵蚀地形因子计算方法，其特征在于，包括：数据输入，获取并输入计算土壤侵蚀地形因子所需的数据；数据处理，包括填洼方法设置、算法选择、阈值设置及模型计算，首先选择填洼方法，选择计算土壤侵蚀因子的算法，然后设置算法所需的阈值参数，之后通过土壤侵蚀模型计算；结果输出，输出土壤侵蚀地形因子数据和图像。2.如权利要求1所述的土壤侵蚀地形因子计算方法，其特征在于，所述填洼方法的算法包括MV算法和快速算法:选择MV算法或快速算法，以及设置限制因子，其中MV算法为将洼地底点栅格高程升高至其流向出口的高程值，从而使水流能够顺畅地流向出口的方法；快速算法为利用优先队列快速填洼的方法；限制因子即填洼的高程大小与当前数据中的最高点的乘积。3.如权利要求2所述的土壤侵蚀地形因子计算方法，其特征在于，所述MV算法和快速算法中的设置包括：坡度计算设置、流向设置、流量计算设置、坡向计算设置和坡长因子计算方法设置，其中坡长因子计算方法包括坡长法和汇流面积法，其中坡向计算设置包括简单算法、南北走向算法和东西走向算法，其中流向计算设置包括流向汇集和快速汇集，其中流向设置中包括D6算法和快速算法，所述坡度计算设置中包括三阶差分。4.如权利要求1所述的土壤侵蚀地形因子计算方法，其特征在于，所述阈值设置中设置的参数包括：坡度、坡长、回流量、去短枝、缓坡中断因子和陡坡中断因子。5.如权利要求1所述的土壤侵蚀地形因子计算方法，其特征在于，填洼方法采用扫描窗口法进行，先找到每个平地栅格或是洼地底点栅格，即潜在出流点，然后找到最低的潜在出流点，比较其和洼地栅格...

【专利技术属性】
技术研发人员：陈志元，陈国坤，黄燕，赵晶晶，胡春田，全诗雨，
申请(专利权)人：昆明理工大学，
类型：发明
国别省市：