基于信息量遥感图像最佳波段选择

《基于信息量遥感图像最佳波段选择研究》简介开始：基于信息量遥感图像最佳波段选择研究摘要：本文介绍基于信息量最佳波段选择中的单波段信息特征量、相关性系数、熵和联合熵、最佳指数、协方差矩阵特征值、波段指数的计算方法及特点，以石家庄市的TM影像为例，分析单波段信息特征量，相关性，最佳指数法，找出了最佳波段组合为1、4、5。

关键词：信息量，最佳指数，波段指数，最佳波段组合。

1引言随着空间技术，数字图像处理技。

论文→ 电子论文》栏目，以上内容为2011-6-8 4:25:20》简单介绍，正文正式开始》》》基于信息量遥感图像最佳波段选择研究摘要：本文介绍基于信息量最佳波段选择中的单波段信息特征量、相关性系数、熵和联合熵、最佳指数、协方差矩阵特征值、波段指数的计算方法及特点，以石家庄市的TM影像为例，分析单波段信息特征量，相关性，最佳指数法，找出了最佳波段组合为1、4、5。

关键词：信息量，最佳指数，波段指数，最佳波段组合。

1 引言随着空间技术，数字图像处理技术和计算机技术的不断发展，遥感技术得到突飞猛进的发展多光谱和高光谱技术的出现，是21世纪遥感技术的发展前沿和当今世界遥感关注的焦点之一，多光谱遥感数据的最佳波段选取是遥感图像增强处理的关键部分，直接影响到目视解译和研究对象的信息提取。

目前遥感图像解译在相当的程度上仍依赖于目视解译．由于人眼对彩色比较敏感且分辨能力强，故应充分利用信息丰富的彩色合成图像进行目标判读．一般的数字图像处理系统都采用三色合成原理形成彩色图像，即在3个通道上安置3个波段图像，然后分别赋以红、绿、蓝色，叠合在一起形成彩色图像[1]。

因此，如何从遥感提供的多光谱数据中快速、准确选取最佳波段，以便于图像的目视解译和信息的有效提取，是遥感数字图像处理的关键问题之一。

本文是基于TM图像信息量的最佳波段选择。

通常选择最佳波段的原则有3点[2]：(1)所选的波段信息量要大；(2)波段间的相关性要小；(3)波段组合对所研究地物类型的光谱差异要大。

那些信息含量多、相关性小、地物光谱差异大、可分性好的波段组合就是最佳组合，据此，可以认为相关性较强的波段组合在一起不会是最佳组合，高光谱遥感数据波段间的存在着不同程度的信息量重复和冗余。

一般选择波段的主要依据是：波段辐射量的方差应尽可能大，因为方差的大小体现了所含信息的多少，但由于地物在各波段的辐射特性之间存在相关性，用3个方差最大的波段合成的效果并不一定能获得更多的信息。

当三者相关性很强时，各波段所包含的信息之间有可能出现大量的重复和冗余。

因此，选择三个波段进行组合时，必须同时考虑方差要大而相关性要小这样两个条件[3]。

2 最佳波段选择的理论模型[4] 目前应用比较广泛的选取方法有各波段信息量的比较、波段间相关性比较、最佳指数法（OIF）、各波段数据的信息熵和联合熵、协方差矩阵特征值法、波段指数法。

2.1单波段信息量的比较根据遥感图像各波段包含的信息量进行数值评价来选择波段是进行波段组合的第一步[5]。

通过分析，可以确定哪几部分或哪几个波段(即波段子集)包含信息量的多少。

各波段的标准差反映了图像各像元灰度值与平均值总的离散度，一定程度上反映了各波段的信息量，其值越大，所包含的信息量越大。

TM 图像各波段所包含的地物信息量，一般采用该波段图像覆盖的辐射量化级，即亮度值范围或称亮度差(最大亮度值最小亮度值)来衡量[6]。

标准差公式为 S= 某幅图像为Aij,(i =1，2，3…… ；j=l，2，3……)，其图像矩阵大小为MxN，Ao代表整幅图像的平均灰度值。

亮度差反映灰度值的变化程度，其大小等于最大亮度值减去最小亮度值，公式:frange(i,j)=fmax(i,j)fmin (i,j)。

均值向量则表征图像中地物的平均反射强度，f= 。

2.2波段间相关性比较信息量最大的三个波段组合的图像，不一定信息量最大，因为在多光谱影像间存在着谱间冗余和空间冗余，所以我们在基于信息量的最佳波段选择时，应该考虑到波段间的相关性。

相关系数描述了波段图像间的相关程度，反映了波段图像所包含信息的重叠程度.反映了的是各通道之间的相关性。

因为图像数据的标准差越大，所包含的信息量也越大，而波段间的相关系数越小，表明各波段图像数据的独立性越高，信息的冗余度越小。

多光谱影像各波段包含的地物信息量固然多寡不一，但基本上有一定的顺序。

各波段之问信息的重迭与分异程度也表现出明显的规律。

相关系数的大小反映的是各通道之间的相关性。

波段间的相关系数反映了两个波段间的信息重叠度，如果两个波段间的相关系数大，则说明它们的信息重叠度高，因此可将相关性大的两个通道合二为一或者取其中一个[7]。

波段问的相关系数公式如下：R=Sij2 /SijSjj 式中：Sij=COV(i,j)= (Yjk )其中i、j为1，2,3，4，5，7为波段数，Sij为第i波段与第j 波段的协方差；为第i波段的光谱灰度均值，为第i波段的第k个像元灰度值，Yjk 为第j 波段的第k个像元灰度值，为第j 波段的光谱灰度均值；k=l，2，3……，n为实验样区的像元数。

R越大，说明波段问信息重叠越大。

波段相关系数矩阵计算相对比较简单，一般的图像处理软件都有该功能。

相关系数分析成了波段选择的基础，相关系数的计算也是最佳指数法的法所必需的。

2.3 最佳指数法（OIF）因为图像数据的标准差越大所包含的信息量也越大，波段间的相关系数越小，表明各波段图像数据的独立性越高，信息的冗余度越小。

故也可采用美国查维茨提出的最佳指数(OIF)的概念。

OIF= /其中：Si为第i个波段的标准差，Sij为i、j两波段的相关系数。

OIF越大，则相应组合图像的信息量越大. 对OIF值进行由大到小排列，最大OIF值对应的波段为最佳波段组合。

此方法是目前最常用的波段选择方法。

计算方法简单，易于操作，且更接近于波段选择的原则[8].2.4 各波段数据的信息熵和联合熵一种既全面、效果又较好且较简便的方法。

信息熵定义来源于信息论和信号处理，它是对信息内容的客观评估[9]。

近年来在遥感数据融合[7][10]处理研究中，信息熵常作为对融合图像信息质量的客观评价指标。

联合熵表征了多个波段组合的信息量，联合熵最大的组合波段即“信息量最大”这个意义上的最佳波段。

但联合熵方法计算量相当大，一般不易操作。

计算时需要编写程序实现。

根据仙农信息论的原理[11]，一幅8bit表示的图像X的熵为：H(x)= (比特)式中为图像像素灰度值为i的概率。

同理， 23幅图像的联合熵分别为：H(x1,x2)=H(x1,x2,x3)=式中Pi1i2i3表示图像X1中像素灰度值为i1 ；图像X2中同名像素灰度值为i2 与图像X3中同名像素灰度值为i3的联合概率。

一般来说，H (x)和H (x1 ，x2 ，x3)越大，图像(或图像组)所含信息越丰富。

用熵和联合熵可以得出信息量大的波段及波段组合。

2.5 协方差矩阵特征值法在多波段遥感图像的波段组合选择中，章孝灿等、Charles [3，12]引提出了用计算Ⅳ维数据熵值原理进行优选最佳波段组合的方法，最后选出了包含信息量丰富的波段组合。

其公式为：H=由于要选取的是三个波段，所以式中Ms= 为三阶方阵。

其中对角线上的元素分别为各波段的方差，非对角线上的元素为两个不同波段的协方差。

因为图像的熵(信息量)H随着变量协方差矩阵行列式的值变化，所以计算任意三波段组合的协方差矩阵行列式的值并进行由大到小排序，就可选出包含信息量丰富的三波段组合。

由协方差矩阵进一步可计算波段间的相关系数矩阵，波段间相关系数由下式计算：Rij= 式中．Sij为波段间的协方差。

利用协方差矩阵特征值法可得到包含信息量丰富的波段，但是它是在完全服从正态分布的情况下得出的模型，选出的波段组合并非完全按照信息量大小排列。

另外这种模型的缺点是计算数据量大，处理速度慢。

2.6 波段指数法波段指数[7] 的定义：设 Rij为通道i与j 之间的相关系数，高光谱数据被分为k组，每组的波段数分别为n1,n2…….nk 。

定义波段指数Pt= ，其中Rt =Rw +Ra，Rw = (i≠j)。

式中Si为第i波段的标准差， Rw为第i波段与所在组内其他波段相关系数的绝对值之和的平均值，Ra 为第i波段与所在组外的其他波段之间的相关系数的绝对值之和。

从公式可以看出，它综合考虑了波段的信息量和相关性两个因素。

标准差越大，表明波段的离散程度越大，所含的信息量越丰富，而波段的总体相关系数的绝对值越小，表明通道数据的独立性越强，信息冗余度越小。

波段指数是从高光谱数据具有的特点出发，综合考虑了图像的信息量和相关性两个因素，并且是在分组成块的情况下设计的模型，而且还发现选出的最佳波段组合与由最佳指数法(考虑光谱特性后)得出的基本一致。

波段指数能综合地反映波段的信息含量和相关性两个因素，可作为波段选择的依据[13]。

因此，波段指数法对高光谱数据来说是一种较理想的最佳波段选择方法。

从波段指数公式可以看出，分子是关于标准差的，分母是关于相关系数的，这与最佳指数很相似，最佳指数的计算简单于波段指数，所以在波段选择时，最佳指数法用的更多。

3 基于信息量的最佳波段选择实验研究3.1 试验背景选择一景轨道号12434的TM影像，时间为1993年10月19日。

ERDAS IMAGINE中的AOI工具切出石家庄市影像。

所有实验在微机上进行，所用软件为ERDAS IMAGINE8.6。

由于TM6是与热有关的波段，所以在波段组合时不考虑，但一起计算了出来。

从信息量和信息冗余这两方面来考虑，最佳指数为最优方法。

各波段组合的最佳指数值同各波段的均方差之和成正比、与各波段间相关系数之和成反比。

满足最佳波段组合选择的基本原则。

又由于其简单的计算方法，所以本文以最佳指数法为例。

3.2 单波段信息量比较从ERDAS中LAYER INFOR 可以得到各波段的统计信息，表1。

表1 TM各波段光谱信息统计表项目TM1TM2TM3TM4TM5TM6TM7MAX1407496103250145254MIN481716961154亮度差9257809424430250均值向量64.81029.65935.32841.28760.613128.80534.436标准差4.7833.4456.0267.67613.9123.7439.638从表中可以看出（1）标准差大小顺序为TM6> TM1〉TM5>TM4>TM3>TM7> TM2（2）亮度差TM7> TM5〉TM4〉TM1>TM3>TM2> TM6(3) 标准差与亮度差大小顺序不一样（4）如果仅仅考虑单波段，TM1信息量最丰富，TM5波段相对也较丰富。