[发明专利]一种用于多卫星传感器的积雪覆盖度混合像元分解方法

权利要求书

1.一种用于多卫星传感器的积雪覆盖度混合像元分解方法，其特征在于，包括：

S1、基于预先设定的多个遥感卫星传感器，获取相应的影像数据；

所述影像数据包括：多光谱地表反射率数据、成像几何以及掩膜数据；

S2、基于所述多光谱地表反射率数据和预先设定的影像端元提取规则进行多种地类的端元提取，获取地类的端元集合；

所述地类的端元集合包括多个地类的端元；

当提取失败则从预先获取的辅助数据中调取地物的地面测量光谱曲线，并根据预先获取的遥感卫星传感器的光谱响应函数，计算匹配传感器波段光谱响应的端元，获取地面测量端元库；

所述预先获取辅助数据包括积雪、植被、土壤光谱曲线，以及光谱响应数据；

S3、基于所述地类的端元集合，获取典型地类端元库；

S4、基于所述典型地类端元库和地面测量端元库，进行多端元光谱混合分析获取积雪覆盖度；

S5、计算卫星像元积雪覆盖度产品。

2.根据权利要求1所述的方法，其特征在于，

所述影像数据包括：MODIS影像数据、Landsat TM/ETM+/OLI影像数据、Sentinel-2MSI影像数据、FY-4A AGRI影像数据以及Himawari-8AHI影像数据。

3.根据权利要求2所述的方法，其特征在于，

所述步骤S2中预先设定的影像端元提取规则为：

所述MODIS影像数据、Sentinel-2/MSI影像数据、Landsat TM/ETM+/OLI影像数据中的积雪端元提取规则为：NDVI-0.035且NDSI0.75且R_0.550.7；植被端元提取规则为：NDVI0.7且NDSI-0.4；土壤和岩石端元的提取规则为：0NDVI0.15且NDSI-0.4；水体端元提取规则为：NDWI＞0.2且R_0.86＜0.2；

Himawari-8/AHI影像数据和FY-4A/AGRI影像数据中的积雪端元提取规则为：NDVI-0.03且NDSI0.74且R_0.550.5；植被端元提取规则为：NDVI0.65且NDSI-0.4；土壤和岩石端元的提取规则为：0NDVI0.15且NDSI-0.4。

4.根据权利要求3所述的方法，其特征在于，所述步骤S3包括：

基于所述地类的端元集合，对所述地类的端元集合中的影像端元进行筛选处理，获取典型地类端元库。

5.根据权利要求4所述的方法，其特征在于，所述步骤S3具体包括：

S31、采用公式(1)计算每一条端元光谱的矢量长度；

其中，||r||为端元光谱的矢量长度；r_k为第k波段的反射率；N为波段数目；

S32、将积雪、植被、土壤地物每一类端元内的所有矢量长度进行升序排列，获取积雪、植被、土壤地物每一类端元相应的序列；

S33、针对所述积雪、植被、土壤地物每一类端元相应的序列，将每一类内端元按照等间隔分割为n个子集；

其中，所述间隔为||R_i||；

n为预先设定值；||R||_max像元最大反射率；||R||_min像元最小反射率；

S34、对每一个子集内的所有光谱取中值或均值，作为所述子集的典型端元；

S35、基于所述每一子集的典型端元，获取典型地类端元库；

所述典型地类端元库库中包括：积雪类地物的n个典型端元、植被类地物的n个典型端元以及土壤类地物的n个典型端元。