课程设计-小波-图像融合

1、目录 TOC o 1-5 h z 摘要11、设计目的与意义22、题目分析 33、设计原理64、总体设计 65、算法设计与功能描述 76、测试结果与分析 107、设计总结11 8、设计体会11参考文献 121 / 13摘要小波变换是一种新的变换分析方法，它继承和发展了短时傅立叶 变换局部化的思想，同时又克服了窗口大小不随频率变化等缺点， 能 够提供一个随频率改变的时间一频率窗口， 是进行信号时频分析和处 理的理想工具。它的主要特点是通过变换能够充分突出问题某些方面 的特征，因此，小波变换在许多领域都得到了成功的应用，特别是小 波变换的离散数字算法已被广泛用于许多问题的变换研究中。从此， 小波变换

2、越来越引进人们的重视，其应用领域来越来越广泛。数据融合是80年代形成和发展起来的一种自动化信息综合处理 技术，它将来自多传感器或多源的信息和数据进行综合处理，从而得出更为准确可信的结论，它充分利用多源数据的互补性和计算机 的高速运算与智能来提高结果信息的质量。图像融合是数据融合技术 在数字图像处理方面的一个应用。高效的图像融合方法可以根据需要 综合处理多源通道的信息，从而有效地提高了图像信息的利用率、 系 统对目标探测识别地可靠性及系统的自动化程度。 本文着重讨论了基 于小波变换的图像融合。关键词：图像融合，小波变换2 / 131设计目的与意义通常地，图像融合是指将来自不同探测器的图像进行合并

3、，以得到一个更为完整的图片或场景。图像融合的主要目的是通过对多幅 图间的冗余数据的处理来提高图像的可靠性，通过对多幅图像间的互补信息的处理来提高图像的清晰度。 高效的图像融合方法可以根据 需要综合处理多源通道的信息，从而有效地提高了图像信息的利用 率、系统对目标探测识别地可靠性及系统的自动化程度。 其目的是将 单一传感器的多波段信息或不同类传感器所提供的信息加以综合，消除多传感器信息之间可能存在的冗余和矛盾， 以增强影像中信息透明 度，改善解译的精度、可靠性以及使用率，以形成对目标的清晰、完 整、准确的信息描述。 图像融合从抽象层次上分为：像素级、特征级 和决策级图像融合。本论文主要研究像素级

4、图像融合 ，研究重点是基 于小波变换的图像融合。由于人的视网膜是在不同的频道中进行处理 因而基于小波变换的融合方法可以获得与人的视觉特性更接近的融 合效果。小波变换将原图像分解成一系列具有不同空间分辨率和频域 特性的子图像，反应了原始图像的局部特征变化，在多个分解层、多 个频带上进行融合。通过小波变换能更好的对图像进行融合， 得到更 好的效果。3 / 132题目分析用小波变换来进行图像融合。图像融合从抽象层次上分为：像素级、特征级和决策级图像融合。 本论文主要研究像素级图像融合，研究重点是基于小波变换的图像融 合。小波(wavelet)是什么在有限时间范围内变化且其平均值为零的数学函数具有有限

5、的 持续时间和突变的频率和振幅在有限的时间范围内，它的平均值等于零小波分析/小波变换：变换目的是获得时间和频率域之间的相互关系小波变换：对一个函数在空间和时间上进行局部化的一种数学变换通过平 移母小波(mother wavelet)获得信号的时间信息通过缩放母小波的 宽度(或称尺度)获得信号的频率特性对母小波的平移和缩放操作是为计算小波的系数，这些系数代表局部信号和小波之间的相互关系对比傅立叶变换：提供了频率域的信息，但丢失了时间域的局部化信息小波分析中常用的三个基本概念：连续小波变换离散小波变换小波重构4 / 13(一)连续小波变换所谓小波(wavelet)是由满足条件:的解析函数经过平移、

6、缩放得到的正交函数族小波变换(WT Wavelet Transform )是用小波函数族ya,b (t)按不同尺度对函数f(t)?L 2 (R)进行的一种线性分解运算：y a J对应的逆变换为：W产m a (f)G cladb小波变换有如下性质：(1)小波变换是一个满足能量守恒方程的线形运算，它把一个 信号分解成对空间和尺度(即时间和频率)的独立贡献，同时又不失 原信号所包含的信息；(2)小波变换相当于一个具有放大、缩小和平移等功能的数学 显微镜，通过检查不同放大倍数下信号的变化来研究其动态特性；(3)小波变换不一定要求是正交的，小波基不惟一。小波函数 系的时宽-带宽积很小，且在时间和频率轴上

7、都很集中，即展开系数 的能量很集中；(4)小波变换巧妙地利用了非均匀的分辨率，较好地解决了时 间和频率分辨率的矛盾；在低频段用高的频率分辨率和低的时间分辨 率(宽的分析窗口)，而在高频段则用低的频率分辨率和高的时间分5 / 13 辨率(窄的分析窗口)，这与时变信号的特征一致；(5)小波变换将信号分解为在对数坐标中具有相同大小频带的 集合，这种以非线形的对数方式而不是以线形方式处理频率的方法对 时变信号具有明显的优越性；(6)小波变换是稳定的，是一个信号的冗余表示。由于 a、b是 连续变化的，相邻分析窗的绝大部分是相互重叠的，相关性很强；(7)小波变换同傅立叶变换一样，具有统一性和相似性，其正

8、反变换具有完美的对称性。小波变换具有基于卷积和QMF勺塔形快速 算法。(二)离散二进小波变换在实际应用中，常常要把连续小波变换离散化。若对连续小波变 换w| (a, b )的伸缩因子a和b进行采样，选取a=2-j , b=2-j kb 0, 则可得到离散的二进小波变换；火品，上式,)=67甲3工-点九)% (瓦/左)=%=卮匚人工)以之x 柩尔这里j, k? Z ,采样率bo 0.由于离散二进小波变换是对连续小波变换的伸缩因子和平移因 子按一定规则采样而得到的，因此，连续小波变换所具有的性质，离 散二进小波变换一般仍具备。(三)小波重构重构概念把分解的系数还原成原始信号的过程叫做小波重构(wa

9、veletreconstruction) 或合成(synthesis),数学上叫做逆离散小波变换 (inverse discrete wavelet transform , IDWT)(四)Mallat算法Mallat算法是便于计算机软件和硬件实现的快速离散算法。这 是Mallat在Burt和Adelson的图像分解和重构的塔式算法的启发 下，根据多分辨率框架提出的算法。此算法在小波分析中的地位相当 于FFT在经典傅立叶分析的地位。按Mallat算法，我们可以把函数f(x)分解为不同频率通道的成 分，并把每一频率通道的成分按相位进行分解，频率越高，相位划分 越细，频率越低，相位划分越粗。Mal

10、lat算法完全是离散的，便于数值计算。6 / 133设计原理小波变换是一种新的变换分析方法，它继承和发展了短时傅立叶 变换局部化的思想，同时又克服了窗口大小不随频率变化等缺点， 能 够提供一个随频率改变的时间一频率窗口， 是进行信号时频分析和处 理的理想工具。小波变换是一个时间和频域的局域变换因而能有效地 从信号中提取信息，通过伸缩和平移等运算功能对函数或信号进行多 尺度细化分析。而通过小波变换处理的图像进行融合， 可得到原有图 像的特性，而不是简单的直接对图像表面进行融合。选着基于小波变换的方法来进行图像融合，提高图像融合的质 量，得到更好的融合效果。4总体设计对待融合的图像用sym4小波基

11、进行分解，然后再对分解图像进行融合，最后对所得到的图像进行逆变换， 就可得到基于小波变换的融合图像。下面是基于小波变换的图像融合流程图1:7 / 135算法设计与功能描述(1)分别对两幅原始图像进行分解首先对两幅原始图像进行分解，得到两组小波分解系数 TOC o 1-5 h z M 1N 1I 1 zr z x za(x, y)= fA(x,y) (x, y)x 0 y 0M 1N 1J 2a(x, y)= j1A(x, y) (x, y)x 0 y 0M 1N 1Ja(x, y)=(x, y) (x, y)x 0 y 0M 1N 1j1 B(x,y)= fB(x,y) (x, y) x 0

12、y 0M 1 N 1j2b(x, y)= j1b(x, y) (x, y) x 0 y 0 M 1 N 1Jb(x, y)= j1b(x, y) (x, y) x 0 y 0得到高频分量和低频分量J a(x, y) = DaJ,AaJ和Jb(x, y)=DBJ,ABJ甘士吉怖公鸟.c Jr J J Jlc Jr J J J,/、干同以刀 JU.DaA , A , A ,DbB , B , B ,记aL a AJ分别为水平方向，垂直方向和对角方向的高频分 量。低频分量：AaJ=SaJ, AbJ=SbJ(2)对分解两幅图像A、B所得的分量进行重构高频分量的确定：首先计算两幅图像在三个方向上的局部能

13、量：8 / 13123)1,2,3)一 ，、 j ,123)1,2,3)Ej,A(m,n)W (m,n)DA (m m,n n)(m M,n N一 ，、.j，、-2,Ej,B(m,n) W (m ,n )DB (m m ,n n )( m M ,n NnWj,maQj,A 叫 minDj,min,Ej,A jDj,叫 minDj,A 叫 maQj,B，Ej,A5其中 Wj,max Wj,min 1.低频分量的确定：处理后的低频分量为Aj，FsA,|sA | |sB |, sA或sB,isA | |sB | sB ,| sA | | sB |.Aj，F(3)图像重构M NfF(x,y)(ADj,

14、F) (x,y)x 0 y 09 / 13MATLA篁现，代码如下：X1=imread(5.jpg);%读取图片 5X2=imread(6.jpg);%读取图片 6%版图片的矩阵行列一致，数字数组的尺寸必须匹配二进制阵列if ndims(X1)=3 % 计算图像X1的维数X3=rgb2gray(X1);% 转换成灰度图 elseX3=X1;endif ndims(X2)=3 % 计算图像X2的维数X4=rgb2gray(X2);% 转换成灰度图 elseX4=X2;end subplot(221) imshow(X3),title(原图像 1);% 显示 X3图像subplot(222)ims

15、how(X4),title(原图像 2);% 显示 X4图像subplot(223);imshow(X3+X4)/2),title(直接进行融合图像结果);%显示直接融合的图像X3=double(X3);%转换成双精度浮点图像X4=double(X4);%行小波变换C1,L1=wavedec2(X3,2,sym4);% 对图像X3用wname、波基函数实现2层分解 C2,L2=wavedec2(X4,2,sym4);% 对图像X4用wname、波基函数实现2层分解 %混像融合 C=(C1+C2)*0.5;X=waverec2(C,L1,sym4);% 多尺度二维小波重构X=uint8(X);%

16、转换为无符号8为整数型图像subplot(224)imshow(X),title(基于小波变换图像融合结果)%显示基于小波变换的图像融合结果10 / 136测试结果与分析运行程序得到如下结果:通过融合结果可以看出基于小波变换的图像融合比直接进行图像融合效果要好很多。基于小波变换的融合图像弥补了 2幅原图像不 同的缺陷，得到完整的清晰图像，采用小波分解融合的方法不会产生 明显的丢失信息现象。而直接进行融合所得的图像灰度值改变与原图像不同。11 / 137设计总结本次设计通过小波变换将原图像分解成一系列具有不同空间 分辨率和频域特性的子图像，反应了原始图像的局部特征变化，在 多个分解层、多个频带上

17、进行融合从而得到较好的融合效果。 通过图 像融合，我们可以看到比较清晰地图像，以及互补原图像的缺点。8设计体会图像融合的重要性，将来自多传感器或多源的信息和数据进行综 合处理，从而得出更为准确可信的结论。图像融合的主要目的是通过 对多幅图间的冗余数据的处理来提高图像的可靠性，通过对多幅图像间的互补信息的处理来提高图像的清晰度。只有更细致的融合才能 得到更完美的图像效果。12 / 13参考文献1陈武凡.小波分析及其在图像处理中的应用.北京：科学出版社. 2002.2林福宗.小波与小波变换.清华大学计算机科学与技术系智能技术与系统国家重点实验室.2001-9-25.3徐建华，王洪华.基于 HIS变换和小波变换的遥感图像融合.信息工程大学测绘学 院.2002 .4晁锐，张科，李言俊.一种基于小波变换的图像融合算法J.电子学报.May. 2004,Vol 32,No 5: 157-159.5 Mallat S. Atheory for multi-resolution signaldecomposition : the wavelet representation J. IEEE Trans . Patt . Anal . Machine Intell .Vol.7,pp.674-693,1989.6Martin Wavelet and Filter