小波变换论文Word下载.docx

1、以Matlab为平台，进行了基于小波变换的图像融合实验，并分析了实验结果。关键词：小波变换 图像压缩 图像去噪 图像融合ABSTRACTThe paper give a brief introduction of wavelet transforms basic conception and analysis the applications of wavelet transform in image compression, image denoising and image fusion. Then it introduces some algorithms about image pro

2、sessing. Finally, give a experiment of image fusion based on wavelet transform, which is programmed in Matlab platform, and analyze the experimental results.Key words: Wavelet transform Image compression Image denoising Image fusion第1章引言当从时域中观察一个信号时，得到的信息是信号随着时间的变化，其幅度的起起伏伏。但是，如果更进一步想研究起伏速度较快或较慢的部分，

3、就不太容易从时域中信号的波形直接得到所需的信息。因此，需要将时域中的信号转换到频域中分析。传统的转换方式是利用傅立叶变换，然而，傅立叶变换潜在的假设了信号是平稳信号。所谓的平稳信号就是信号的规律不随时间的变化而改变，而现实生活中的信号往往是非平稳信号和平稳信号交织在一起的。另一方面，用傅立叶变换提取信号的频谱需要利用信号全部时域的信息，也就无法通过傅立叶分析来刻画时域信号的局部特性。为解决傅立叶变换的不足，Gabor提出在傅立叶变换中加入高斯窗函数，将窗函数沿时间轴挪移，得到一系列包含时间信息的傅立叶变换结果，从而能同时分析信号的时间信息和频率信息。根据Heisenberg的测不准原理，窗口傅

4、立叶变换对信号的时间定位和频率定位能力是相互矛盾的，时间分辨率和频率分辨率不可能同时提高，而且变换窗口没有自适应性，只适于分析所有特征尺度大致相同的信号，不适于分析多尺度信号和突变过程。由此，引入了小波变换。顾名思义，“小波”就是小的波形。所谓“小”是指它具有衰减性，而称之为“波”则是指它的波动性，其振幅呈正负相间的震荡形式。傅立叶分析是将信号分解成一系列不同频率的正弦波的叠加，同样小波分析是将信号分解为一系列小波函数的叠加，而这些小波函数都是由一个母小波函数经过平移和尺度伸缩得来的。小波分析优于傅立叶分析的地方是，它在时域和频域同时具有良好的局部化性质，且具有多分辨分析的特点。它是一种窗口大

5、小可以改变的分析方法，可以改变其时间窗和频率窗，根据高频和低频的不同，可以使时间频率窗变窄或变宽，即：在低频部分时具有较高的频率分辨率和较低的时间分辨率，在高频部分时具有较低的频率分布率和较高的时间分辨率，非常适合于加带、瞬态、反常现象的探测正常信号中并展示其成分，所以被誉为分析信号的显微镜。第2章小波变换的基本概念2.1连续小波变换给定基本小波函数，信号f（t）的连续小波变换定义为： （a0,bR） （2.1）式（2.1）也可以表示为，它可以看做是求函数f（t）在的各尺度平移信号上的投影。其给出了f（t）的一种多尺度表示，a代表尺度因子，=称为小波。若a1，则函数具有伸展作用，a1时，函数具

6、有收缩作用。因此随着参数a的变化，就有可能实现窗口大小自适应变化，当信号频率增高时，时窗宽度变窄，而频窗宽度增大，有利于提高时域分辨率，反之亦然。参数b反映的位移，把基本小波（母小波）的函数作位移后，再在不同尺度下与待分析信号作内积，就可以得到一个小波序列。然而，小波函数的选择既不是唯一的，也不是任意的，它应满足以下几个条件：1）定义域应是紧支撑的，即在一个很小的区间外，函数为零，也就是函数应有速降特性。2）函数平均值为零，而函数只有在t轴上取值有正有负才能保证均值为零，所以函数应有震荡性。在实际应用尤其是信号处理以及图像处理的应用中，变换只是一种简化问题、处理问题的有效手段，最终目的需要回到

7、原问题的求解，因此，还要保证连续小波变换存在逆变换。对于所有的f（t）、 （R），连续小波逆变换由式（2.2）给出：f（t）= （2.2）2.2离散小波变换连续小波变换往往用于理论分析，在实际应用中，需要将连续小波加以离散化才能应用于计算机中。因此，需要对尺寸参数a和平移参数b进行离散化处理，设定 ，则： （2.3）定义相应的离散小波变换为： （2.4）其逆变换为： （2.5）2.3二维小波变换一维小波变换是将一维原始信号分别经过低通滤波和高通滤波以及二元下抽样得到信号的低频部分L和高频部分H。根据Mallat算法，二维小波变换可以用一系列的一维小波变换得到。对一幅m行n列的图像，二维小波变换

8、的过程是先对图像的每一行做一维小波变换，得到L和H两个对半部分；然后对得到的LH图像（仍是m行n列）的每一列做一维小波变换。这样经过一级小波变换后的图像就可以分为LL，HL，LH，HH四个部分。而二级、三级以至更高级的二维小波变换则是对上一级小波变换后图像的LL部分再进行一级二维小波变换，是一个递归过程，从而得到塔式结构的图像，如下图所示：一个图像经过小波分解后，可以得到一系列不同分辨率的子图像，不同分辨率的子图像对应的频率也不同。高分辨率（即高频）子图像上大部分点的数值都接近于0，分辨率越高，这种现象越明显。2.4多分辨分析多分辨分析又称为多尺度分析，是建立在函数空间概念上的理论。其主要思想

9、是将L （ R ）分解为一串具有不同分辨率的子空间序列，该子空间序列的极限就是L （ R ），然后将L （ R ）中的 f 函数描述为具有一系列近似函数的逼近极限，其中每一个近似函数都是 f 函数在不同分辨率子空间上的投影。通过这些投影可以分析和研究f 函数在不同分辨率子空间上的形态和特征。类似于人的视觉系统，多分辨率分析在各尺度上可以由粗到精地逐级观察目标。利用多分辨率分析能构造L （ R ）的标准小波正交基。第3章小波变换在图像处理中的应用利用小波变换进行图像处理均需要有三个基本的步骤，先通过小波变换将图像变换到频域中，再对小波系数进行相应的处理，最后对处理后的小波系数进行重构，得到处理后

10、的图像。针对不同的应用，具体的小波系数处理过程会有所不同。3.1小波变换在图像压缩中的应用经过小波变换很容易得到图像的低频部分和高频部分，图像的大部分能量一般集中于低频部分，而高频部分则反应图像的细节。因此，最简单的压缩方式，就是在小波重构时，将高频系数置0，也可根据图像将局部区域的高频系数置0，或者根据阈值来决定高频系数的取舍。这种方法图像细节损失严重，压缩后的图像模糊，图像信息损失严重。由于遥感数据的海量增加、保存的重要性及其应用对精度保证的要求，因而遥感数据的压缩必须满足信息保持压缩及能达到较高的压缩比这两个条件时才有意义。小波变换能提供原始图像的多尺度分解表示，而且各分辨率层上的子图像

11、具有不同的频率特征和不同的方向取向，从而可对这些信息表示进行相应的编码。图像信号的统计特征表明大幅值的系数往往集中于低频区内，这样就可以给那些小幅值系数分配很小的比特数，甚至可以不传输或存储，得到很高的压缩比和很小的失真度。通常认为遥感数据中存在有空间冗余和谱间冗余，前者表现为同一类地物相邻像素间存在的空间相关性，后者表现为相邻波段同一位置的像素之间存在相关性。因此有学者提出用遥感图像零树压缩编码算法和基于整数小波变换的压缩等方法来去除遥感数据中的冗余，达到更高的压缩比和更小的失真度。遥感图像零树压缩编码算法的原理是子代中某一小波系数不重要的时候,便认为其所有后代的小波系数都是不重要的,即形成

12、所谓的零树。通过对零树进行编码得到压缩后的遥感影像。这种算法假设不重要的小波系数其后代也不重要，事实上，对于图像边缘往往不满足假设条件，因此存在压缩影像损失边缘信息的情况。为解决这一问题，有许多人提出了改进算法，如Said 和Pearlman 提出的SPIHT （Set Partitio2ning In Hierarchical Trees） 算法。基于整数小波变换的压缩方法将小波变换压缩技术中的零树编码推广到高光谱图像压缩中，其基本原理是利用高光谱图像的结构相关性，对多幅小波系数图像，只构造一幅有效图（共享有效图） 来确定多幅小波系数图像中非零值的位置，通过各波段零树相“与”得到共享零树，同

13、时去除空间冗余和谱间冗余。还可以先进行K- L 变换再进行共享有效图的小波变换编码，进一步去除谱间冗余,提高压缩效率。3.2小波变换在图像去噪中的应用在图像中噪声一般表现为高频信息，并且与图像中其他信息呈弱相关或者不相关。所以消噪过程一般可按以下方法进行处理。首先对图像进行小波分解，选择小波并确定分解层次为N，则噪声部分通常包含在高频中。然后对小波分解的高频系数进行门限阈值量化处理。最后根据小波分解的第N层低频系数和经过量化后的1N层高频系数进行小波重构，达到消除噪声的目的，即抑制图像的噪声。基于小波分析的去噪方法大体有小波萎缩法、投影方法、相关方法等几种。小波萎缩法的出发点是较大的小波系数通

14、常反映景物信号，而较小的系数通常反映噪声。因此可以设定阈值，过滤掉噪声。阈值函数的设定决定了降噪的效果，广泛使用的阈值函数有硬阈值函数和软阈值函数。此外，还可以通过判断系数被噪声污染的程度，并为这种程度引入各种度量方法（例如概率和隶属度等） ，进而确定萎缩的比例，来实现去噪的目的。投影方法的原理是将带噪信号以一种迭代的方式，投影到逐步缩小的空间，由于最后的空间能更好地体现原信号的特点，所以投影方法也能够有效地区分噪声和信号。其中典型的算法是MatchingPursuits 法，通过指定一族小波或波函数，并将带噪声信号向此函数进行投影，接着又对残差投影，并循环反复，直到残差符合一定限值。相关方法

15、是基于信号在各层相应位置上的小波系数间往往相关、而噪声的小波系数则具有弱相关或不相关的特点来去噪。3.3小波变换在图像融合中的应用在图像融合中，小波变换的目的是将原始图像分别分解到一系列频率通道中,利用分解后的金字塔或树结构对不同分解层、不同频带进行融合处理,这种方法有助于将来自不同图像的感兴趣的细节融合在一起。针对不同类型的图像,学者们提出了各种有效的融合规则如取系数绝对值较大法、加权平均法、消除高噪声法、高低频混合双阈值法等。取系数绝对值较大法适合高频成分较丰富，亮度、对比度较高的原图像，否则在融合图像中只保留一幅图像的特征，其他的特征被覆盖；融合图像中基本保留原图像的特征，图像对比度与原

16、图像基本相同。小波变换的实际作用是对信号解相关，并将信号的全部信息集中到一部分具有大幅值的小波系数中。这些大的小波系数含有的能量远比小系数含有的能量大，从而在信号的重构中，大的系数比小的系数更重要。加权平均法的权重系数可调，适用范围广，可消除部分噪声，原图像信息损失较少，但会造成图像对比度的下降，需要进行图像灰度增强。消除高频噪声法的高频噪声可以基本消除，融合图像对比度较高，原图像特征可较好地保留在融合图像中，但在消除高频噪声的同时，损失了部分高频信息。双阈值法适于原图像中一幅图像的灰度分布均衡，高频成分较多。双阈值可选，增加了算法的实用性，但选择阈值时要考虑原图像灰度分布的特点，否则有可能出

17、现边缘跳跃的现象。以上是基于单个像素的融合规则，其在融合处理时表现出对边缘的高度敏感性，使得在预处理时要求图像是严格对准的，否则处理结果将不尽人意，这就加大了预处理的难度。而基于区域的融合规则由于考虑了与相邻像素间的相关性，降低了对边缘的敏感性，所以具有更加广泛的适用性。基于区域特征的融合规则主要包括基于梯度的方法、基于局域方差的方法、基于局域能量的方法等。3.4其他应用此外小波变换还能应用于图像拼接与镶嵌、特征提取与分类、图像复原、图像插值、纹理分析和边缘检测等方面。第4章基于小波变换的图像融合实验根据文献中的各种图像融合方法，我利用Matlab做了一个简单的图像融合实验。4.1算法原理1）

18、将待融合的两幅影像进行小波变换2）对于高频小波系数采用基于像素点小波系数绝对值取较大的规则进行融合。3）进行小波重构，得到融合影像。4.2实验结果4.3实验分析可以看到融合后的影像基本上只保证了第二幅影像的特征，但是在亮度上比第二幅影像更亮，而趋近于第一幅影像。第一幅影像的信息只有很少一部分反映在融合图像中，融合效果好像并不是太理想。这个可能与图像的选择有关，也与算法的设计有关。总之，通过实验基本了解了利用小波变换进行图像融合的过程，有了比较直观的体会。便于以后进行深入研究。参考文献1 贾永红.数字图像处理M. 武汉大学出版社, 20042 林宏裔, 孔亮. 在 MATLAB 语言环境下基于小波变换的图像处理J. 华北科技学院学报, 2003, 5（2）: 60-64.3 刘金亨. 基于小波变换的遥感图像处理研究D. 重庆大学, 2010.4 葛伟. 基于小波的图像处理技术的实现与探讨J. 现代计算机 （专业版）, 2002 （5）.5 何国金. 小波变换在遥感图像处理中的应用综述J. 遥感信息, 1999 （1）: 14-17.6 徐建达, 王洪华. 基于 IHS 变换和小波变换的遥感影像融合J. 测绘学院学报, 2002, 19（3）: 198-202.7 晁锐, 张科, 李言俊. 一种基于小波变换的图像融合算法J. 电子学报, 2004, 32（5）: 750-753.