[发明专利]黄土高原中尺度流域年侵蚀产沙模型计算方法

权利要求书

1.一种黄土高原中尺度流域年侵蚀产沙模型计算方法，其特征在于，其具体过程为；

步骤1，把流域划分为若干单元，每个单元近似的代表一个坡面；

步骤2，利用坡面土壤侵蚀试验观测研究成果和认识，基于坡面土壤侵蚀试验观测数据，建立各因子与侵蚀量之间关系，即单元模型；

步骤3，用所述单元模型完成对每个单元侵蚀量的计算；

步骤4，利用空间分析和空间统计方法，通过聚合运算，完成对流域土壤侵蚀总体状况的评价和分析；

上述步骤2中，获取各因子的计算过程如下：

步骤21，计算站点气候因子R1；

步骤22，计算气候因子计算插值R2；

步骤23，计算植物措施因子B，用土地利用图和植被盖度图，计算某区域植物措施因子B；

步骤24，计算工程措施因子E，利用统计数据计算各区域的工程措施因子E值；

步骤25，计算耕作措施因子T，耕作措施因子值T根据不同坡度条件下等高耕作减少土壤流失来确定。

2.根据权利要求1所述的黄土高原中尺度流域年侵蚀产沙模型计算方法，其特征在于，在上述步骤21中，计算站点气候因子R1包括一个或多个站点的降雨侵蚀力，其根据相应站点的降雨量来计算；

计算公式如下述(1)和(2)所示：

F=(ΣPi2)/P---(1)]]>

R=αF^β    (2)

式中，P为年平均降雨量(mm)；Pi为第i月的平均降雨量(mm)；R为多年平均降雨侵蚀力(MJ·mm·hm^-2·h^-1·a^-1)；α=0.3589，β=1.9462；F指数mm大小与年平均雨量P的季节分布有关，取值范围在P·12-1～P之间。

3.根据权利要求1或2所述的黄土高原中尺度流域年侵蚀产沙模型计算方法，其特征在于，在上述步骤23中计算植物措施因子B的具体过程为：

步骤231，获取土地利用图；

步骤232，计算植被盖度图；

步骤233，计算年平均植被盖度；

步骤234，计算植物措施因子B。

4.根据权利要求3所述的黄土高原中尺度流域年侵蚀产沙模型计算方法，其特征在于，在上述步骤232中，利用N DVI根据下述公式(3)计算植被盖度，

F=100×float(NDVI)-NDVIminNDVImax-NDVImin---(3)]]>

式中，F为植被盖度，NDVI为归一化植被指数，NDVImin和NDVImax分别为无植被地区的N DVI值和植被良好覆盖地区的NDVI值；在现场测试基础上，结合遥感图像典型地类的采样，确定最高和最低植被盖度的NDVI阈值；将采样的结果与对应的NDVI值进行统计比较，当NDVI值>≥NDVImax时，盖度可视为100％，当NDVI值＜NDVImin时，盖度可视为0％。