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

说明书

技术领域

本发明属于信息提取方法，具体涉及一种用于伊利石信息提取的高光谱影像处理方法。

背景技术

通过对图像信息处理来提取图像区域内不同类型矿物的丰度，是目前勘探领域常用的方法。这种方法不需要工作人员直接到野外勘探、采样，可以提取不同类型矿物的丰度。当前的高光谱遥感影像的伊利石信息提取方法主要是光谱全波段匹配或是部分连续波段的光谱匹配，具体方法有光谱角、混合解调滤波、沙漏提取方法等，由于地表的物质组成很少是由单一矿物组成的，这些方法在信息提取的过程易受其他地物波谱或噪声的影响，提取信息精度相对较低。其次现有的光谱提取方法人工操作步骤多，增加了人为判断误差。第三是高光谱数据波段多，数据量大，现有的方法处理时间长，降低了数据处理的速度和应用规模。因此，如何在伊利石信息提取的过程中减少其他地物或噪声的影响、人工操作步骤和处理数据量，成为当前高光谱遥感影像处理的前沿之一。

发明内容

本发明针对现有技术的缺陷，提供一种用于伊利石信息提取的高光谱影像处理方法。

本发明是这样实现的：一种用于伊利石信息提取的高光谱影像处理方法，包括

步骤一：波段采样

对已有的地面反射率的影像数据进行采样，提取在1415nm、1430nm、1940nm、1970nm、2165nm、2195nm、2225nm、2240nm、2285nm、2345nm的影像，并依次记录为b1～b10，每次采样得到的都是一幅灰度图，图中每一个像元的值都是其灰度值，

步骤二：判断

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

a1＝(b1小于b2)；

a2＝(b3小于b4)；

a3＝(b5大于b6)；

a4＝(b7小于b8)；

a5＝(b9大于b10)；

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

步骤三：计算

用下述公式进行计算b11

b11＝b2+b4+b5+b8+b9-b1-b3-b6-b7-b10

所述的*表示相乘，

用下面公式计算a0，

a0＝a1*a2*a3*a4*a5*b11

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

本发明的效果是：本发明在对高光谱影像重采样，提取特定波段，进行一系列判断和计算，可以准确的计算出影像范围内的不同区域伊利石的丰度值。保留伊利石的光谱特征明显的波段，去除其他特征不明显的波段，从而在信息提取的过程突出伊利石的光谱特征，降低其他地物或噪声的影响，减少了处理的数据量，并可以用IDL程序达到用较少的人工操作实现最终结果信息提取的目的，提高了伊利石信息提取的精度和速度。

具体实施方式

下面结合实施例对本发明作进一步说明。一种用于伊利石信息提取的高光谱影像处理方法，包括

步骤一：波段采样

对已有的地面反射率的影像数据进行采样，提取在1415nm、1430nm、1940nm、1970nm、2165nm、2195nm、2225nm、2240nm、2285nm、2345nm的影像，并依次记录为b1～b10，即b1为波段1415nm的采样数据，b2为波段1430nm的采样数据，以此类推。每次采样得到的都是一幅灰度图，图中每一个像元的值都是其灰度值，即b1为一幅灰度图，b1图像的(1,1)点的值为灰度值，其余点依此类推，其余采样图也依次类推。

步骤二：判断

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

a1＝(b1小于b2)；

a2＝(b3小于b4)；

a3＝(b5大于b6)；

a4＝(b7小于b8)；

a5＝(b9大于b10)；

上述判断是针对每次判断图像的相应像元进行判断的，以a1＝(b1小于b2)为例，取b1图像的某像元(例如(1,1)点)的灰度值，与b2图像的相应像元(当b1图像取(1,1)点，则b2图像也必须取(1,1)点)的灰度值，然后根据判断规则“b1小于b2”判断，当判断结果为“是”时，记录判断结果为1，否则记录结果为0。因此当a1＝(b1小于b2)判断完毕时，得到的a1是与b1矩阵相同大小的矩阵，其中每个点的值是根据判断规则得到的判断结果(即每个点的值是0或1)。

其它判断也按照类似的规则进行。本步骤结束后得到a1～a5，共5个矩阵。