[发明专利]一种秸秆覆盖和人为扰动改进RULSE-C因子的方法

权利要求书

1.一种秸秆覆盖和人为扰动改进RULSE-C因子的方法，其特征在于，包括以下步骤：

S1：获取土壤数据，并对其进行预处理，其中，土壤数据包括：卫星遥感数据、秸秆覆盖度野外实测数据、气象观测数据、水文数据、土壤栅格数据和再分析数据；

S2：构建归一化耕作指数来反演作物秸秆覆盖度，以提高秸秆覆盖度的精度；

S3：根据步骤S1和S2的数据计算降雨侵蚀力因子、土壤可蚀性因子、地形因子、植被覆盖与管理因子、土壤保持因子；其中，植被覆盖与管理因子分为非耕地的植被覆盖与管理因子和耕地的植被覆盖与管理因子，耕地的植被覆盖与管理因子通过作物生育期和休耕期的覆盖因子以及耕地管理因子来计算；

S4：根据RUSLE模型计算年土壤侵蚀模数。

2.如权利要求1所述的秸秆覆盖和人为扰动改进RULSE-C因子的方法，其特征在于，步骤S3中，耕地的植被覆盖与管理因子C_cropland的计算公式如下：

C_cropland＝C_crop·C_CRC+C_disturbance

式中：C_cropland为耕地的植被覆盖与管理因子；C_crop和C_CRC分别为作物生育期和休耕期的覆盖因子；C_disturbance是反映人为扰动强度的耕地管理因子。

3.如权利要求2所述的秸秆覆盖和人为扰动改进RULSE-C因子的方法，其特征在于，步骤S3中，作物生育期的覆盖因子C_crop，计算公式如下：

式中：C_crop为耕地的作物生育期的覆盖因子；c为植被覆盖度；NDVI为归一化植被指数，取值范围为–1至1；NDVI_max和NDVI_min分别是累积分布概率为5％和95％时的NDVI值。

4.如权利要求2所述的秸秆覆盖和人为扰动改进RULSE-C因子的方法，其特征在于，步骤S3中，休耕期的覆盖因子C_CRC的计算公式为：

C_CRC＝e^(-0.03·CRC)

式中，C_CRC为休耕期的覆盖因子的值；CRC为基于步骤S2计算得到的作物秸秆覆盖度。

5.如权利要求2所述的秸秆覆盖和人为扰动改进RULSE-C因子的方法，其特征在于，步骤S3中，耕地管理因子C_disturbance的计算公式如下：

式中：C_disturbance为表征人为扰动强度的耕地管理因子，取值范围为0至0.20；SS为表征人为扰动因子的县域综合得分，(SS–minSS)/(maxSS–minSS)为归一化后的SS值，maxSS和minSS分别为SS最大值和最小值。

6.如权利要求5所述的秸秆覆盖和人为扰动改进RULSE-C因子的方法，其特征在于，基于主成分分析法计算的县域综合得分SS，其计算公式如下：

SS＝F₁·VCR

式中，F₁为PCA值；VCR为PCA的方差贡献率，为90.50％。

7.如权利要求2所述的秸秆覆盖和人为扰动改进RULSE-C因子的方法，其特征在于，步骤S4中，年土壤侵蚀模数的计算公式为：

A＝R·K·LS·C·P

式中：A为年土壤侵蚀模数；R为降雨侵蚀力因子；K为土壤可蚀性因子；LS为地形因子；C为植被覆盖度和管理因子；P为土壤保持因子。