[发明专利]一种用于石膏信息提取的高光谱影像处理方法

权利要求书

1.一种用于石膏信息提取的高光谱影像处理方法，其特征在于：包括下述步骤

步骤一，预处理

获取高光谱影像；对高光谱影像的进行预处理，进行大气校正，获取地面反射率的影像数据；

步骤二：采样

对步骤一得到的影像数据进行采样，对波段在1145nm，1175nm，1325nm，1445nm，1460nm，1700nm，1745nm，1820nm，1940nm，2060nm,2210nm，2330nm，2405nm的图像采样，并依次记录为b1～b13，每次采样得到的都是一幅灰度图，图中每一个像元的值都是其灰度值，即b1为波段1145nm的采样数据，b2为波段1175nm的采样数据，以此类推；

步骤三：判断

进行下面一系列判断，并记录结果

a1＝(b1大于b2)；

a2＝(b3大于b4)；

a3＝(b4小于b5)；

a4＝(b6大于b7)；

a5＝(b7小于b8)；

a6＝(b8大于b9)；

a7＝(b9小于b10)；

a8＝(b10大于b11)；

a9＝(b12大于b13)；

上述判断是针对每次判断图像的相应像元进行判断的，当判断结果为“是”时，记录判断结果为1，否则记录结果为0，

本步骤结束后得到a1～a9，共9个矩阵；

步骤四：计算用下述公式进行计算b14

b14＝b1+b3+b5+b6+2*b8+2*b10+b12-b2-2*b4-2*b7-2*b9-b11-b13

所述的*表示相乘；

用下面公式计算a0

a0＝a1*a2*a3*a4*a5*a6*a7*a8*a9*b14

a0为石膏信息的丰度图；

上述所有计算均为相应像元计算，即使用不同图像的相应像元计算。