[发明专利]基于特征指标的矿物成分遥感信息提取方法及装置

权利要求书

1.基于特征指标的矿物成分遥感信息提取方法，其特征在于包括如下步骤：

A、对遥感图像进行预处理；

B、对经过步骤A处理的遥感图像的像元的特征指标进行计算，所述特征指标包括：

波段相关斜率SL：SL=（B_i-B_j）/（B_j-B_m）；

波段斜率指数SI：SI=e^((Bi-Bm)/Bi)；

波段相关吸收深度指数SD：SD=e^{((Bi-Bk)/(Bi-Bm))}；

波段相关多项式指数SP：SP=e^{((Bi-Bm)/Bi+(Bj-Bk)/Bm)}；

其中，B_i为所述像元的反射率或DN值，对应的波段为样本波段光谱图中反射率或DN值的最大值对应的波段；B_j为所述像元的反射率或DN值，对应的波段为样本波段光谱图中反射率或DN值的第二大值对应的波段；B_k为所述像元的反射率或DN值，对应的波段为样本波段光谱图中反射率或DN值的第三大值对应的波段；B_m为所述像元的反射率或DN值，对应的波段为样本波段光谱图中反射率或DN值的最小值对应的波段；

C、建立回归方程：

Con_ij=a×SL_ij+b×SI_ij+c×SD_ij+d×SP_ij+e

其中，Con_ij为像元处的矿物成分含量，SL_ij为所述像元的波段相关斜率，SI_ij为所述像元的波段斜率指数，SD_ij为所述像元的波段相关吸收深度指数，SP_ij为所述像元的波段相关多项式指数，a、b、c、d为系数，e为回归分析截距；

D、根据实际样品反推出所述回归方程中a、b、c、d和e的值；

E、根据步骤C得到的回归方程计算像元处的矿物成分含量。

2.根据权利要求1所述基于特征指标的矿物相关遥感信息提取方法，其特征在于：还包括步骤F：基于步骤E得到的所述像元处的矿物成分含量形成彩色图像。

3.根据权利要求2所述基于特征指标的矿物相关遥感信息提取方法，其特征在于：在步骤C中还包括对所述回归方程进行回归检验的步骤。

4.根据权利要求1至3之一所述基于特征指标的矿物相关遥感信息提取方法，其特征在于：步骤A所述预处理包括如下步骤：

选取成像时间为春季或秋季的遥感图像作为原始数据；对属于ASTER数据的所述原始数据选择6、3和1波段组合形成第一优化数据，或对属于ETM数据的所述原始数据选择7、4和3波段组合形成第一优化数据。

5.根据权利要求4所述基于特征指标的矿物相关遥感信息提取方法，其特征在于：步骤A所述预处理还包括如下步骤：在所述第一优化数据中选取研究区域：分割第一优化数据，判断分割区内每个波段的直方图是否为正态分布，如果是正态分布，则将该分割区作为研究区域；如果不是正态分布，则将所述分割区划分为多个次级区域，使得所述次级区域内的每个光谱波段的直方图为正态分布，并将所述次级区域作为所述研究区域。

6.根据权利要求5所述基于特征指标的矿物相关遥感信息提取方法，其特征在于：步骤A所述预处理还包括如下步骤：将所述原始数据中不包含全部波段数据的边界区域的去除步骤。