[发明专利]土壤侵蚀植被覆盖-管理因子的遥感估算方法

权利要求书

1.一种土壤侵蚀植被覆盖-管理因子的遥感估算方法，其特征在于其实施步骤如下：

1)将遥感图像进行土地利用/覆盖类别划分得到多种地类；将所述地类进行分类，得到空间异质性小、季节变化小、年际变化小的第一组地类，空间异质性小、季节变化大、年际变化小的第二组地类，空间异质性大、季节变化小、年际变化大的第三组地类；

2)采用遥感分类法获取第一组地类和第二组地类中每一个地类的植被覆盖-管理因子；

3)将遥感图像进行线性光谱混合像元分解，得到绿色植被分量图、土壤分量图和阴影/水分量图；通过C＝F_土壤/(1+F_绿色植被+F_阴影)计算得到第三组地类的植被覆盖-管理因子，其中，F_土壤为土壤分量，F_绿色植被为绿色植被分量，F_阴影为阴影/水分量；

4)将各个地类的植被覆盖-管理因子进行综合得到遥感图像的植被覆盖-管理因子。

2.根据权利要求1所述的土壤侵蚀植被覆盖-管理因子的遥感估算方法，其特征在于：所述步骤4)中将遥感图像进行线性光谱混合像元分解的详细步骤包括：

A)采用主成分分析、最大噪声比变换或者奇异值分解对遥感图像进行数据降维；

B)选择端元模型，根据所述端元模型确定端元的数量和光谱；

C)根据所述端元的数量和光谱对数据降维后的遥感图像进行反演得到绿色植被分量图、土壤分量图和阴影/水分量图。

3.根据权利要求2所述的土壤侵蚀植被覆盖-管理因子的遥感估算方法，其特征在于：所述步骤C)中对数据降维后的遥感图像采用全约束最小二乘法进行反演，所述全约束最小二乘法的模型为其中，R_ib为第b波段第i像元的光谱反射率，N为端元数目，f_ki为对应于i像元的第k个端元所占分量值，C_ib为第b波段端元K的光谱反射率，ε_ib为该像素第b波段的分解残差；所述模型的约束条件为且0≤f_ki≤1，且所述模型通过误差ε_i来评价模型的优劣，误差ε_i的表达式为其中m为光谱波段数目。

4.根据权利要求1或2或3所述的土壤侵蚀植被覆盖-管理因子的遥感估算方法，其特征在于：所述多种地类包括密林地、建设用地、水体、水田、旱地、疏林地、园地、裸土地八种地类；所述第一组地类包括密林地、建设用地、水体；所述第二组地类包括水田、旱地；所述第三组地类包括疏林地、园地、裸土地。