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

说明书

技术领域

本发明涉及一种用于黄钾铁矾信息提取的高光谱影像处理方法， 特别是涉及一种降低其他地物或噪声的影响，减少了处理的数据量的 用于黄钾铁矾信息提取的高光谱影像处理方法。

背景技术

当前的高光谱遥感影像的黄钾铁矾信息提取方法主要是光谱全 波段匹配或是部分连续波段的光谱匹配，具体算法有光谱角、混合解 调滤波等，由于地表的物质组成很少是由单一矿物组成的，这些方法 在信息提取的过程易受其他地物波谱或噪声的影响，提取信息精度相 对较低。其次现有的光谱提取方法人工操作步骤多，增加了人为判断 误差。第三是高光谱数据波段多，数据量大，现有的方法处理时间长， 降低了数据处理的速度和应用规模。因此，如何在黄钾铁矾信息提取 的过程中减少其他地物或噪声的影响、人工操作步骤和处理数据量， 成为当前高光谱遥感影像处理的前沿之一。

发明内容

本发明要解决的技术问题是提供一种降低其他地物或噪声的影 响，减少了处理的数据量的用于黄钾铁矾信息提取的高光谱影像处理 方法。通过对高光谱影像重采样，提取特定波段，进行一系列判断和 计算，计算出影像范围内的不同区域黄钾铁矾的丰度值。

为解决上述技术问题，本发明一种用于黄钾铁矾信息提取的高光 谱影像处理方法，依次包括：

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

第二步，对地面反射率的影像数据进行重采样，提取波段在 980nm，1355nm，1415nm，1460nm，1550nm，1595nm，1790nm， 2015nm，2075nm，2120nm，2225nm，2255nm，2330nm，2405nm 的影像；

第三步，当980nm波段影像的像元值小于1355nm波段影像相对 应的像元值，判断该像元满足条件一；获得1355nm波段影像减去 980nm波段影像的结果影像H1；

第四步，当1415nm波段影像的像元值大于1460nm波段影像相对 应的像元值，该像元满足条件二；获得1415nm波段影像减去1460nm 波段影像的结果影像H2；

第五步，当1550nm波段影像的像元值小于1595nm波段影像相 对应的像元值，判断该像元满足条件三；获得1595nm波段影像减去 1550nm波段影像的结果影像H3；

第六步，当1790nm波段影像的像元值大于2015nm波段影像相 对应的像元值，判断该像元满足条件四；获得1790nm波段影像减去 2015nm波段影像的结果影像H4；

第七步，当2075nm波段影像的像元值大于2120nm波段影像相 对应的像元值，判断该像元满足条件五；获得2075nm波段影像减去 2120nm波段影像的结果影像H5；

第八步，当2120nm波段影像的像元值小于2225nm波段影像相 对应的像元值，判断该像元满足条件六；获得2225nm波段影像减去 2120nm波段影像的结果影像H6；

第九步，当2225nm波段影像的像元值大于2255nm波段影像相 对应的像元值，判断该像元满足条件七；获得2225nm波段影像减去 2255nm波段影像的结果影像H7；

第十步，当2255nm波段影像的像元值小于2330nm波段影像相 对应的像元值，判断该像元满足条件八；获得2330nm波段影像减去 2255nm波段影像的结果影像H8；

第十一步，当2330nm波段影像的像元值大于2405nm波段影像 相对应的像元值，判断该像元满足条件九；获得2330nm波段影像减 去2405nm波段影像的结果影像H9；

第十二步，选择同时满足条件一至条件九的像元，获得上述像元 范围内的H1、H2、H3、H4、H5、H6、H7、H8、H9相加的和H。

本发明本方法只用了14个波段，相对于高光谱影像SASI全波段 101个波段，需要处理的数据量减少了86％，并且由于是计算机自动 一步提取，减少了主成分变换、端元波谱的选择等操作步骤，运算速 度可以提高了7倍以上。由于去除了大部分对信息提取关系不大的波 段，减少其他物质或噪声对其光谱的干扰，提高了信息提取的精度。 对高光谱影像数据中黄钾铁矾信息的快速提取具有较好的作用和意 义。

具体实施方式