基于小波变换的遥感影像融合与评价

1、    基于小波变换的遥感影像融合与评价        林 卉1， 杜培军2 时间:2008年09月08日     字 体: 大 中 小        关键词:        摘要：关键词： 小波变换 影像融合 统计参数 影像评价1 基于小波变换的图像融合算法在现实中，图像中的物

2、体是以不同的尺度出现的。以一条边缘为例，它可能是一条从黑到白的陡峭边缘，也可能是一条跨越相当距离的缓变边缘。图像表示和分析中的多分辨率方法就是基于这种考虑。Mallat在Burt和Adelson图像分解与重构金字塔算法(即高斯-拉普拉斯金字塔算法)的启发下，基于小波变换的多分辨率分析，提出了Mallat算法。小波变换等价于用低通滤波器和高通滤波器对图像信号进行滤波处理。若Hr,Gr和Hc,Gc分别表示镜像共轭滤波器H和G分别作用在行和列上，这样小波变换可以简单表示如下：对于一幅二维图像，式(4)中算子Hr,Hc是一个低通滤波器，它过滤出图像的低频分量，因此Cj+1显示了Cj低频分量，即图像的低

3、频部分。算子Hr,Gc相当于二维低通滤波器，它对列作平滑，检测行的差异，因此D1j+1显示Cj的垂直方向的高频分量，即图像的水平边缘。算子GrHc对行作平滑，检测列的差异，因此D2j+1显示Cj的水平方向的高频分量，即图像的垂直边缘。算子GrGc是两个方向的高频滤波，检测对角边缘。由此看出，对一副影像进行小波变换，就是将其分解到不同频率下的不同特征域上。把待融合的原始图像进行小波变换后，图像被分解到不同的频率区域上，图像的融合处理就要在每一频率段上分别采用不同的算法进行图像的融合。对于图像可以进行层小波变换，每一层的小波变换只需对上一层的小波变换后的低频分量进行变换。这样就形成了小波变换的金字

4、塔结构。由于图像融合的最高层需要对数据进行选取或均衡，因此最高层的低频部分所用的融合算子是融合图像的细节取舍的最关键一步，对图像高频算子的合理选取可以起到增强图像边缘、突出边缘的作用。一般，最高层的低频部分采用比较算子，高频部分采用简单的加权算子，其他各层均采用简单的加权算子。2 小波特性分析与图像融合规则2.1融合方法的理论根据通过数值分布统计,信源图像和经过小波分解后的子图像具有如下特性:(1)原图像中区域的数据变化幅度与在子图像中相应区域的数据变化幅度一致；(2)对于同一目标或物体的不同信源图像,其低频图像相应区域的数据值相同或相近,而高频子图像却有显著差别。小波变换的上述特性,为有效融

5、合方法的选择提供了理论依据。2.2 图像融合规则及融合因子采用武汉市某地区的两幅TM遥感图像进行融合实验。设A、B分别为两幅原始图像,F为融合后的图像。在图像融合过程中,融合规则及融合算子的选择是图像融合中至今尚未解决的问题。(1)融合规则一：对于边缘分量,即小波分解中的高频分量LHi,HLi,HHi,取两幅图像相应系数矩阵中对应项的最大值(i=1,2,3,，N);对于低频分量LL,由于这部分对恢复图像质量影响很大,采用：F(j,k)=(A(j,k)+K×B(j,k)×-A(j,k)-K×B(j,k)×计算。其中K、为加权因子。前半部分(A(j,k)+K

6、×B(j,k)×表示取两幅图像的加权均值,影响融合后图像的能量,对融合后图像的高度起决定作用;后半部分A(j,k)-K×B(j,k)×表示取两幅图像的加权差值,包含两幅图像的模糊信息。因子K调节两幅图像的占优比例,使两幅亮度不同的图像达到均衡。随着因子的增大,图像加亮;随着因子增大,图像的边缘加强。对于不同的图像,适当调整K、,可以消减模糊边缘并确保在消减时不过度地丧失边缘信息。对其他高频成分,取两组系数的最大值,可以得到最强的边缘信息,从而得到质量良好的输出图像。(2)融合规则二：分别计算两幅图像小波分解后的高频子图像上每个像素的局部平均梯度,以将像素

7、的局部平均梯度作为准则,确定融合后的高频子图像上的像素值。设A(x,y)和B(x,y)分别为可见光和红外图像,不同分辨率高频子图像分别为： Ajk(x,y)和Bjk(x,y)(k=1,2,3； j为尺度参数)不同分辨率的高频子图像的梯度图像为GAjk(x,y)和GBjk(x,y),若GAjk(x,y)GBjk(x,y),则小波分解后的高频信息为： 3 融合结果的定量评价，其中pi为灰度等于i的像素数与总的像素数的比。通过对图像信息熵的比较可以对比出图像的细节表现能力。熵的大小，反应了图像携带的信息量的多少。融合图像的熵值越大，说明融合图像携带的信息量越大。如果图像中所有灰度级出现概率趋于相等，

8、则包含的信息量趋于最大。(3)图像交叉熵直接反映了两幅图像对应像素的差异，可以用来确定各种融合效果的优劣，交叉熵值越小，说明融合结果图像与两幅原始图像的差异越小。令P=p1,p2,p3pm，Q=q1,q2,q3qm,则P与Q的交叉熵定义为：CEN=-piln(pi/qi)，设CEN1和CEN2为融合前两幅影像与融合结果图像的交叉熵，则融合图像与两幅原始图像的总体交叉熵为：。下面用n表示图像大小，K表示多光谱图像中的光谱分量的个数，k表示第k个光谱分量，分别为第k个光谱分量在原始图像与融合图像(i,j)点的颜色灰度值。(4)扭曲程度。图像光谱扭曲程度直接反映了多光谱图像的光谱失真程度。第k个光谱

9、分量的光谱扭曲定义为：(5)偏差指数。偏差指数用来比较融合影像和低分辨率多光谱影像偏离程度。第k个光谱分量的偏差指数定义为：(6)相关系数。图像的相关系数反映了两幅图像的相关程度。通过比较融合增强前后的图像相关系数可以看出多光谱图像的光谱信息改变程度。第k个光谱分量融合前后的相关系数定义为：4 实验结果与分析实验采用的两幅源图像为武汉某地区在2003年所摄取的同一地区TM4、TM5波段遥感影像，大小均为289×289，如图1所示。它们在TM七个波段的影像中的信息熵和标准方差是最大的，影像最清晰，质量最好。从视觉上看，TM4影像光线较暗，影像很模糊，对地物的轮廓反映较差，尤其对河中的孤

10、岛表达不完整。而TM5影像比TM4影像要清晰，地物的轮廓反映较突出，河中孤岛的轮廓较清晰，信息量很丰富，图像质量要比TM4影像好。表1中的均值、方差、信息熵的定量统计数值也说明了这点。对这两幅影像分别采用基于像元的加权融合、基于主分量变换(K-L)的融合、小波变换规则一和规则二融合，融合后的影像如图2所示。从融合结果图可以看出，无论是空间分辨率，还是图像的清晰程度，这几种融合方法的结果图都比两幅原始影像提高很多。图2(a)中影像模糊，地物的轮廓不突出，光线很暗，对河中孤岛的反映很清晰，这种简单的处理会引起信息的失真，造成边缘模糊。图2(b)中影像比图2(a)清晰，光线要亮，对细节部分的体现更明

11、显，尤其对湖中养鱼网箱的体现，但很难识别河中的孤岛。图2(c)中影像的信息有较大提高，光线较明亮，影像中的地物清晰可见，道路、网箱、孤岛的轮廓非常清晰，容易辨识。图2(d)中的影像视觉效果最佳，整幅影像的地物清晰，层次分明，信息量丰富，孤岛、道路和湖中的网箱非常清晰。利用原始影像与融合影像的亮度分量统计得到均值、方差、信息熵、交叉熵。对于光谱信息，利用影像的RGB三个光谱分量统计得到各光谱分量的扭曲程度、偏差指数和相关系数，并计算了平均扭曲程度、平均偏差指数和平均相关系数。通过对两幅原始影像和融合后的影像进行统计，其性能统计参数如表1所示。从表1可看出，这四种融合方法中反映影像亮度信息的均值是

12、递增的；反映空间信息细节的参数，如方差和信息熵呈现递增的趋势，交叉熵呈递减的趋势；反映光谱信息的参数如各个光谱的扭曲程度与偏差指数是递增的，而各个光谱的相关系数递减，平均扭曲程度、平均偏差指数和平均相关系数也有类似的变化趋势。在这几种融合方法中，最后一种方法得到的图像交叉熵最小，熵值最大，使得图像的信息量增加，与原图像差异最小，图像微小细节反差和纹理变化特征的表达能力大大增强，图像的清晰度得到了很大幅度提高。这些统计参数的变化趋势与各种融合方法的理论分析结果一致，与人的视觉效果的变化一致。对于多光谱遥感影像的融合处理，若用一般的图像融合方法和规则或者融合算子选择不当，就很难获得满意的融合效果。

13、通过对几种融合规则的融合结果比较，提出了规则二的融合方法，其算法简单，稳定性好，图像增强效果明显。实验结果表明，该融合规则取得了良好的效果。同时，对比融合图像统计参数表明，它们的变化趋势与目视结果是一致的。说明这些参数可以作为多光谱影像融合的客观效果评价准则，统计特性法是一种客观、实用的分析与评价方法。参考文献1 Jo Yew Tham,Shen Luan Lee,et al.A general approach for analysis and application of discrete multi-wavelet transformsJ.IEEE Trans on Signal Pro

14、cessing,2000;48(2):4574642 Pohl C ,Van Garderen J L.Multisensor Image Fusion in Re-mote Sensing:Concepts,Methods, Application J.Int J Remote Sensing,1998;19(5):3623673 Zhong Zhang ,Blum R S A.Categorization of multiscale-de-compositon-based image fusion schemes with a performance study for a digital camera applicationJ.Proceeding of the IEEE,1999;87(8):130