专利名称：：适用于高光谱遥感图像的自动波段选择方法
技术领域：
：本发明涉及一种适用于高光谱遥感图像的自动波段选择方法，通过应用此方法，可以将信息丰富而与其它波段相关性小的波段选择出来，为后续的地物精确分类和目标探测提供信息丰富而波段数少的数据源。在各类军用或者民用的遥感信息处理系统、数字城市空间信息系统、航空航天机载高光谱成像仪等领域中均可有广泛的应用。技术背景目前对于高光谱遥感图像的处理方法主要集中在数值分析处理上，如大气校正、信息提取、压縮、目标探测与分类等。而高光谱遥感图像降维方法的研究是做好后续处理的一个关键步骤，降维方式的正确选取与使用，对于发展和完善那些针对高光谱海量数据和丰富信息等特点的算法和软件有极大的好处。现有的降维方法，可以分为两类一类是基于变换的方法，如主成分分析、正交子空间投影、正则分析和离散小波变换等。基于变换的降维方式的优点是可以经过若干变换直接将高维数据降低到几维甚至--维，降维速度快；其缺点是对图像进行了变换，改变了图像原有的特性。另一类是基于非变换的，如波段选择，数据源划分等。基于非变换的降维方式是在考察图像整体特点之后对图像进行选择和划分的，它克服了变换法改变图像特性的缺点，因此更有利于保持图像的原有特性。目前的波段选择方法，一般都是对于多光谱图像(4-7个波段)进行组合选优的。最初的波段选择主要依赖于人工目视解译，这个工作需要由有经验的专家来进行。但是随着遥感图像数字化的不断发展，依赖于专家目视解译的方法越来越显示出它的局限性。例如，高光谱波段数目通常达到100个以上，人工判读的过程会比较慢，影响了影像的实时处理；对一个新的图像数据进行选择容易出现人为的判断失误等。目前，有一些方法对多光谱或者高光谱遥感图像进行波段选择。如利用JM距离和一般地表覆盖类型的光谱集合的专家系统来进行波段选择。然而通过进一步的研究发现该方法不仅需要对地表覆盖类型和成像时间的关系有着深入的了解，而且要发展有效的JM距离选择策略。这些要求对于高光谱波段选择，显然是不适用的。另外还有一种基于类间可分性的最佳组合波段选择方法，但它所涉及的波段选择研究都是对于特定的波段组合进行的。在实际应用中，由于高光谱图像维数过多，如果从220个波段中选择出100个波段，那么其组合数为《2°。°，如果要想成功的辨识某些地物，比如湖泊，就要考虑到波段的排列问题，那么可能的波段排列总数为尸。对于上述所有可能的波段组合，再求其相应的Bhattacharyya距离、离散度、协方差矩阵等，计算量将会十分巨大。高光谱遥感图像各波段间存在高相关性和高冗余度，通过选择或者提取有用的波段或者特征可以大大地降低计算量，从信息处理实效性的角度来看，波段选择是有必要的。并且很多实际处理的例子也表明，并非波段越多越好，在高光谱图像预处理时去除有些受噪声影响的波段对于后继处理也是大大有益的。因此，对于数字化的高光谱图像进行准确的波段选择，必须发展有效而简便的识别方法。
发明内容本发明的目的在于针对现有技术的上述不足，提供一种适用于高光谱遥感图像的自动波段选择(Automaticbandselection)降维的新方法，应用此方法可以将信息含量丰富并且与其他波段相关性小的波段选择出来。该方法简单有效，计算量小。不仅能降低后继处理的计算量，还可以将信息集中到几个波段中，有利于目标探测和分类。自动波段选择(Automaticbandselection)方法充分考虑了各波段的空间相关性和光谱间相关性，对求出的各个波段指数进行由大到小的排列，系统根据设定的阈值自动选择需要的波段。自动波段选择方法依据的原则是(1)所选择的波段信息量要大；(2)所选择的波段与其他波段的相关性要小。本发明的方法包括如下具体步骤1)计算高光谱遥感图像各个波段的局部标准差，记为a,，/=1，一,"，"为高光谱遥感图像的波段总数；a,的具体计算公式如公式(l):<formula>formulaseeoriginaldocumentpage5</formula>其中，m、iv分别是图像的行、列像素总数，yxx,力是第/个波段图像(x,力处的像素值，",是第；个波段图像的像素平均值。计算局部标准差的目的是考察待目标探测区域的信息富集程度，比应用整个图像的方差更为准确。2)计算高光谱遥感图像各个波段之间的相关系数。(lSKn),^具体计算公式如公式(2):<formula>formulaseeoriginaldocumentpage5</formula>这里乂(x,力为波段/在(;c,力处的像素值，乂+1(1,力为波段/+1在(^，力处的像素值，w,为第个波段的像素平均值，",+,为第/+1个波段的像素平均值。所有r,形成一个相关系数矩阵/;3)将当前波段与其前后波段的相关系数相加并取平均值，艮P:d,,+':.,+1)/2(5)其中，rM,是第/个波段与前一个波段的相关系数，r,,+|是第/个波段与后一个波段的相关系数，相关系数越小，两个波段数据之间的独立性越高。4)计算波段指数/mfec，具体计算公式为/mfex,=(6)5)通过执行以上歩骤，计算机可自动计算出每个波段所对应的指数值。本发明的方法最终波段数目的选择包含两种选择方法。方法一根据实际目标探测和分类的需要以及计算机快速处理图像的能力选择指数排在前面的"("s")个波段，"为高光谱遥感图像的波段总数。方法二根据实际目标探测和分类的需要以及计算机快速处理图像的能力选择指数大于某个阈值的前A(A:^")个波段，w为高光谱遥感图像的波段总数。确定了阈值之后，指数大于该阈值的波段就被选择出来。本发明的特色在于1)充分利用了高光谱遥感数据的空间信息和光谱间的信息；2)可计算机自动执行，适合于星载或者机载高光谱成像仪实吋处理数据；3)算法简单有效，计算量小，最终波段数目选择方法简单灵活。4)应用自动波段选择方法，后继目标探测和分类的计算量可大大降低。附图1为220波段的AVIRIS高光谱遥感图像的第10波段图像，该图像取自1992年6月拍摄的美国印第安纳州西北部印第安遥感试验区的一部分。具体实施方式为了更好地理解本发明的技术方案，以下结合附图对本发明的实施方式做进一步描述。附图为220波段的AVIRIS高光谱遥感图像的第IO波段图像，该图像取自1992年6月拍摄的美国印第安纳州西北部印第安遥感试验区的一部分。去掉受噪声和水汽吸收较明显的谱段，从原始的220个波段中选取了180个作为研究对象。它包含了农作物和森林植被的混合区。在应用了上述图像之后，本发明的方法各部分的具体实施细节如下1)计算高光谱遥感图像各个波段的局部标准差，记为(T,，/=1，一,"，n为高光谱遥感图像的波段总数180;cx,的具体计算公式如公式(7):1MW丄其中，M、iV分别是图像的行、列像素总数，/是第/个波段图像(x，力处的像素值，",是第Z个波段图像的像素平均值。2)计算高光谱遥感图像各个波段之间的相关系数。(lS/《")，r,具体计算公式如公式(8):其中，<formula>formulaseeoriginaldocumentpage7</formula>(8)(9)(10)这里乂(x，力为波段z'在(;c，力处的像素值，/;+1(;^)为波段/+1在(1，力处的像素值，w,为第/个波段的像素平均值，",+1为第Z+l个波段的像素平均值。所有。形成一个相关系数矩阵及；通过自动波段选择方法获得的180波段的部分相关系数矩阵如表1，可以看出相邻波段的相关性一般比较大，随着波段距离的增大相关系数相应地减小。表1谱带B28-B39的相关系数矩阵<table>tableseeoriginaldocumentpage7</column></row><table>3)将当前波段与其前后波段的相关系数相加并取平均值，艮口;，，=(r'-U+。,w)/2(11)其中，rM,,是第/个波段与前一个波段的相关系数，是第/个波段与后一个波段的相关系数，相关系数越小，两个波段数据之间的独立性越高。4)计算波段指数/^fec，具体计算公式为/"麵=~^(12)用自动波段选择方法得出的指数见表2，其中/n&x由大到小排列，后面是对应的波段号(该波段号为在180个波段集中的序号)。表2/"A;c及其对应的波段号(前18个)序<table>tableseeoriginaldocumentpage8</column></row><table>5)最终波段数目的选择。本发明的方法最终波段数目的选择包含两种选择方法。方法一根据实际目标探测和分类的需要以及计算机快速处理图像的能力选择指数排在前面的减"S")个波段，《为高光谱遥感图像的波段总数。方法二根据实际目标探测和分类的需要以及计算机快速处理图像的能力选择指数大于某个阈值的前A(A^W)个波段，M为高光谱遥感图像的波段总数。确定了阈值之后，指数大于该阈值的波段就被选择出来。通过执行以上步骤，计算机可自动将所需的波段选择出来。权利要求1.一种适用于高光谱遥感图像的自动波段选择方法，适用于高光谱遥感图像的降维处理，其特征在于包含如下具体步骤1)计算高光谱遥感图像各个波段的局部标准差，记为σi，i＝1，…，n，n为高光谱遥感图像的波段总数，σi的具体计算公式如公式(1)2.根据权利要求l所述的高光谱遥感图像的自动波段选择方法，其特征在于最终波段数目的选择可以根据实际目标探测和分类的需要以及计算机快速处理图像的能力选择指数排在前面的减"^n)个波段，"为高光谱遥感图像的波段总数。3.根据权利要求l所述的高光谱遥感图像的自动波段选择方法，其特征在于最终波段数目的选择还可以根据实际目标探测和分类的需要以及计算机快速处理图像的能力选择指数大于某个阈值的前"A:《")个波段，"为高光谱遥感图像的波段总数。确定了阈值之后，指数大于该阈值的波段就被选择出来。全文摘要一种适用于高光谱遥感图像的自动波段选择方法，有利于地物精确分类和目标探测。该方法充分考虑了高光谱图像各波段的空间相关性和光谱间相关性，使得信息丰富并且与其他波段相关性小的波段被选择出来，降低了后继处理的计算量，适合于高光谱成像仪实时处理数据。最终波段选择的方法有两种一种是系统根据设定的阈值选择进行目标探测与分类的波段数目，阈值的设定可以根据后续的处理需要自动地选择适合的指数值，确定了阈值之后，指数大于该阈值的波段就被选择出来；另一种方法就是选择波段指数排在前面的d(d≤n)个波段，n为高光谱遥感图像的波段总数。该方法可通过计算机自动执行。文档编号G06K9/36GK101131734SQ20071011780公开日2008年2月27日申请日期2007年6月25日优先权日2007年6月25日发明者刘春红,赵慧洁申请人:北京航空航天大学