小波分析论文

1、基于小波分析在图像处理中的应用摘要介绍了图像小波分析的基本理论和基于小波变换的分解与重构原理，利用小波变换对二 维图像进行分解，将原始图像分解成不同方向、不同频率成分的子图像。同时对含噪声图像 进行小波分解。通过选取适当的阈值，对小波分解系数进行阈值量化，再对高低频系数重构， 实现图像的去噪。最后运用 MATLA仿真平台进行仿真验证，仿真结果表明：利用小波分析对 图像进行压缩和去噪可以得到非常好的压缩效果和去噪效果。对工程应用具有一定的借鉴意 义。关键字：小波；图像压缩；图像去噪；1 引言小波分析 (Wavelet Analysis) 即小波变换是 80 年代中期发展起来的一门新兴的数学理论

2、和方法，它被认为是傅立叶分析方法的突破性进展，它具有许多优良的特性。小波变换的基 本思想类似于 Fourier 变换，就是用信号在一族基函数张成的空间上的投影表征该信号。经 典的 Fourier 变换把信号按三角正、余弦基展开，将任意函数表示为具有不同频率的谐波函 数的线性迭加，能较好地刻划信号的频率特性，但它在时空域上无任何分辨，不能作局部分 析，这在理论和应用上都带来了许多不便。小波分析优于傅立叶之处在于，小波分析在时域 和频域同时具有良好的局部化性质，因为小波函数是紧支集，而三角正、余弦的区间是无穷 区间，所以小波变换可以对高频成分采用逐渐精细的时域或空间域取代步长，从而可以聚焦 到对象

3、的任意细节。因此，小波变换被誉为分析信号的显微镜，傅立叶分析发展史上的一个 新的里程碑。小波分析的应用是与小波分析的理论研究紧密地结合在一起的。现在，它已经在科技信 息领域取得了令人瞩目的成就。现在，对性质随时间稳定不变的信号，处理的理想工具仍然 是傅立叶分析。但在实际应用中，绝大多数信号是非稳定的，小波分析正是适用于非稳定信 号的处理工具。图像处理是针对性很强的技术，根据不同应用、不同要求需要采用不同的处 理方法。采用的方法是综合各学科较先进的成果而成的，如数学、物理学、心理学、信号分 析学、计算机学、和系统工程等。计算机图像处理主要采用两大类方法：一类是空域中的处 理，即在图像空间中对图像

4、进行各种处理；另一类是把空间与图像经过变换，如傅立叶变换， 变到频率域，在频率域中进行各种处理，然后在变回到图像的空间域，形成处理后的图像。 图像处理是“信息处理”的一个方面，这一观点现在已经为人所熟知。它可以进一步细分为 多个研究方向：图片处理、图像处理、模式识别、景物分析、图像理解、光学处理等等。小 波分析用在图像处理方面，主要是用来进行图像压缩、图像去噪、图像增强(包括图像钝化和图像锐化 )、图像融合、图像分解。2 常用小波介绍2.1 Haar 小波A.Haar 于 1990 年提出一种正交函数系 , 定义如下：10 兰 XE1/2屮H丿-11/2 兰 xv1(2.1 )0其它这是一种最

5、简单的正交小波，即屮(t)屮(x -n)dx =0n = ±1,±2,(2.2)2.2 Daubechies (dbN)小波系该小波是Daubechies从两尺度方程系数h *出发设计出来的离散正交小波。一般简写为 dbN, N是小波的阶数。小波和尺度函数吁中的支撑区为 2N-1。的消失矩为除N= 1 外(Haar小波),dbN不具对称性即非线性相位；dbN没有显式表达式(除N= 1外)。但 ?N A的传递函数的模的平方有显式表达式。假设P(y)八CkNJkyk，其中，C"-1k为二项式的系”k=0数，则有-m0(«>) 2 =(cos2 竺)N

6、P(sin2?)(2.3)2 21 2N A其中 m°(co)=送 hkeJk°y2 k=02.3 Biorthogonal(biorNr.Nd )小波系Biorthogo nal函数系的主要特征体现在具有线性相位性，它主要应用在信号与图像的重 构中。通常的用法是采用一个函数进行分解，用另外一个小波函数进行重构。Biorthog onal函数系通常表示为biorNr.Nd的形式：Nr=1Nd=1，3,5Nr=2Nd=2，4,6, 8Nr=3Nd=1，3,5, 7, 9Nr=4Nd=4Nr=5Nd=5Nr=6Nd=8其中，r表示重构，d表示分解2.4 Coiflet (coi

7、fN )小波系coiflet 函数也是由Daubechies构造的一个小波函数，它具有coifN (N=1, 2，3，4，5) 这一系列，coiflet 具有比dbN更好的对称性。从支撑长度的角度看，coifN具有和db3N及 sym3N相同的支撑长度；从消失矩的数目来看，coifN具有和db2N及sym2N相同的消失矩数目。2.5 SymletsA (symN 小波系Symlets函数系是由Daubechies提出的近似对称的小波函数，它是对db函数的一种改进。Symlets函数系通常表示为symN(N=2 3, , ,8)的形式。2.6 Mexican Hat(mexh)小波2 2Mexi

8、can Hat 函数为空(x)=笫兀少4(1 x2)e/2(2.4)它是Gauss函数的二阶导数，因为它像墨西哥帽的截面，所以有时称这个函数为墨西哥 帽函数。墨西哥帽函数在时间域与频率域都有很好的局部化，并且满足：(x)dx=O(2.5)a由于它的尺度函数不存在，所以不具有正交性。2.7 Meyer 函数4 二<3_ 38 二 _Meyer小波函数？和尺度函数都是在频率域中进行定义的，是具有紧支撑的正交小波< co34兀< co3(2兀)"2e2 sin(工口(旦园1)22兀(2.6)0(国)=«(2巧"电曲2 cos(=u(2莎卜°)

9、22兀0其中，：(a)为构造Meyer小波的辅助函数，且有J 、-1/22兀(2.7)33(2巧?()= (2 二)/2cos( (f|-1)22兀03小波分析用于图像压缩3.1图像压缩概述通常所说的图像压缩主要指无损压缩(无失真)和有损压缩(有失真)两大类。所谓无 损压缩是指图像数据经压缩后可以完全得到复原，复原后的图像与原始图像完全一致。有损 压缩则是指经它处理的数据在基本保持原图像的特征的前提下，不可避免地要丢掉一部分原 始图像信息。图像能够进行压缩的主要原因是：(1) 原始图像信息存在着很大的冗余度，数据之间存在着相关性，如相邻像素之间色彩 的相关性等，消息中这些冗余信息将会产生额外的

10、编码。如果去掉冗余信息，就会减少消息 所占的空间。(2) 在美图系统的应用领域中，人眼作为图像信息的接收端，其视觉对于边缘急剧变化 不敏感(视觉掩盖效应)，以及人眼对图像的亮度信息敏感，而对颜色分辨率弱等，因此在高 压缩比的情况下，解压缩后的图像信号仍比较满意。基于上述两点，无论采用无损压缩还是有损压缩。只要损失的数据不太影响人眼主观接 受的效果，即可采用。一个图像作小波分解后，可得到一系列不同分辨率的子图像，不同分 辨率的子图像对应的频率是不相同的。高分辨率(即高频)子图像上大部分分点的数值都接 近于0，越是高频这种现象越明显。对一个图像来说，表现一个图像最主要的部分是低频部 分，所以一个最

11、简单的压缩方法是利用小波分解，去掉图像的高频部分而只保留低频部分。3.2主要调用命令“whoS'用于显示当前 MATLA工作空间的变量，而在命令窗口中输入data后，将显示该数据。变量查询函数 who与whos,作用都是列出在 matlab工作空间中已经驻留的变量名 清单，不同的是whos在给出驻留变量的同时，还给出他们的维数及性质。wavedec2是多尺度二维小波分解，调用格式为：C,L = wavedec2(X,N,'w name')即对信号X进行N尺度的小波分解，wname为所使用的小波名称。N为正整数。输出分 解结构包括行向量C,它包含计算出的小波变换系数及定义

12、了 C中系数的排列的记录矩阵L。图3.1图像压缩流程图4 小波分析用于图像去噪4.1 图像去噪概述噪声可以理解为妨碍人的视觉器官或系统传感器对所接收图像源进行理解或分析的各种 因素。一般噪声是不可预测的随机信号，它只能用概率统计的方法去认识， 。噪声对图像处理 十分重要，它影响图像处理的输入、采集、处理的各个环节以及输出结果的全过程。特别是 图像的输入、采集的噪声是个十分关键的问题，若输入伴有较大噪声，必然影响处理全过程 及输出结果。因此一个良好的图像处理系统，不论是模拟处理还是计算机处理无不把减少最 前一级的噪声作为主攻目标。去噪已成为图像处理中极其重要的步骤。对二维图像信号的去噪方法同样适

13、用于一维信号，尤其是对于几何图像更适合。二维模 型可以表述为s(i,j)=f( i,j)+ S e(i,j) i,j=0,1, ,m-1其中， e 是标准偏差不变的高斯白噪声。二维信号用二维小波分析的去噪步骤有 3 步：(1) 二维信号的小波分解。选择一个小波和小波分解的层次N,然后计算信号s到第N层的分解。( 2)对高频系数进行阈值量化。对于从 1 到 N 的每一层,选择一个阈值,并对这一层的 高频系数进行软阈值量化处理。(3) 二维小波的重构。根据小波分解的第N层的低频系数和经过修改的从第一层到第N层的各层高频系数计算二维信号的小波重构。在这 3 个步骤中,重点是如何选取阈值和阈值的量化。

14、4.2 主要调用命令dde ncmp的调用格式有以下三种：(1) THR,SORH,KEEPAPP,CRIT=ddencmp(IN1,IN2,X)(2) THR,SORH,KEEPAPP,CRIT=ddencmp(IN1,'wp',X)(3) THR,SORH,KEEPAPP,CRIT=ddencmp(IN1,'wv',X)函数ddencmp用于获取信号在消噪或压缩过程中的默认阈值。输入参数X为一维或二维信号；IN1取值为'den'或'cmp' , 'den'表示进行去噪，cmp'表示进行压缩；IN2取值为

15、'wv'或'wp'， wv表示选择小波，wp表示选择小波包。返回值 THR是返回的阈值；SORH是软阈值或硬阈值 选择参数；KEEPAP表示保存低频信号；CRIT是熵名(只在选择小波包时使用)。wdencmp用于一维或二维信号的消噪或压缩。其调用格式为：1. XC,CXC,LXC,PERF0,PERFL2=wdencmp('gbl',X,'wname',N,THR,SORH,KEEPAPP)2. XC,CXC,LXC,PERF0,PERFL2=wdencmp('lvd',X,'wname',N,TH

16、R,SORH)3. XC,CXC,LXC,PERF0,PERFL2=wdencmp('lvd',C,L,'wname',N,THR,SORH)wname是所用的小波函数，gbl(global的缩写)表示每层都采用同一个阈值进行处理，lvd表 示每层用不同的阈值进行处理，N表示小波分解的层数，THR为阈值向量，对于格式(2)(3) 每层都要求有一个阈值，因此阈值向量THR的长度为N, SOR表示选择软阈值还是硬阈值(分 别取为s'和h')，参数KEEPAP取值为1时，则低频系数不进行阈值量化处理，反之，则低频系数进行阈值量化。XC是消噪或压缩后的信

17、号，CXC, LXC是XC的小波分解结构,PERF0和PERFL2是恢复和压缩LA2的范数百分比。4.3程序流程图图4.1图像去噪流程图5运行结果5.1图像压缩结果压缩前图像X的大小：Name SizeBytes ClassX 256x256524288 double arrayGrand total is 65536 eleme nts usi ng 524288 bytes 第一次压缩图像的大小为：Name SizeBytes Classca1 135x135145800 double arrayGrand total is 18225 eleme nts usi ng 145800 by

18、tes 第二次压缩图像的大小为：Name S ize Bytes Classca2 75x7545000 double arrayGrand total is 5625 elements using 45000 bytes 运行结果如图所示：'原始图像分解后低频和高频信息100200300400500100 200 300 400 500 20曽 »40j601 j VWk i .80jIr100J QF120第一次压缩20 40 60 80 100120第二次压缩图4.1利用二维小波分析进行图像压缩图像对比如图所示。可以看出，第一次压缩提取的是原始图像中小波分解第一层的低频

19、 信息，此时压缩效果较好，压缩比较小：第二次压缩是提取第一层分解低频部分的低频部分 （即小波分解第二层的低频部分），其压缩比较大，压缩效果在视觉上也基本过的去。这是一 种最简单的压缩方法，只保留原始图像中低频信息，不经过其他处理即可获得较好的压缩效 果。在上面的例子中，我们还可以只提取小波分解第3、4、，层的低频信息。从理论上说，我们可以获得任意压缩比的压缩图像。MATLAB中实现图像压缩，还可利用现有的函数来实现。这种方法主要包括获取压缩阈 值和进行图像压缩两反面。实现获取压缩阈值的函数有ddencmp和wdcbm2两个，实现图像压缩的函数有 wdencmp、wpdencm 和 wthcoe

20、f2 三个。5.2图像去噪结果输出结果从图中5个图像的比较可以看出，Matlab中的ddencmp和wdencmpS数可以有 效地进行去噪处理原始图像X含噪声的图像Xnoise基于软阈值去噪后的图像5010015020050 100 150 20050 100 150 2005010015020050 100 150 200硬阈值去噪后的图像50100150默认阈值去噪后图像50 100 150 200图5.2图像去噪比较小波阈值法去噪：主要适用于信号中混有白噪声的情况。其优点是噪声几乎完全得到抑 制，且反映原始信号的特征尖峰点得到很好的保留。用软阈值法去噪可以使去噪信号是原始 信号的近似最优

21、估计，且估计信号至少和原始信号同样光滑而不会产生附加振荡。6结论本文对图像压缩和去噪的相关技术进行分析，并根据小波分析的原理解决图像压缩和去 噪的方法，详细论述了小波分析函数在图像压缩和去噪的应用，通过MATLA仿真结果表明，利用此方法进行图像的压缩，去噪具有比较好的效果。随着数字图像处理技术的发展，图像 压缩的应用已经深入到关系国计民生的许多领域。而基于小波分析变化的图像压缩的应用会 更显著。在今后的发展中，小波分析变换以其独特的分解重构算法，在图像处理应用中会越 来越广泛。参考文献1 张德丰.数字图像处理M.北京：人民邮电出版社，2009: 186-220 .2 王经民.小波分析.杨凌：西

22、北农林科技大学出版社，2004.43 高成.MATLAB波分析与应用M.北京：国防工业出版社，2007.4 张志勇，彭玉青.小波分析在图像处理中的应用J.江苏：常州工学院学报，2005,12刘贵忠，邸双亮.小波分析及其应用.西安：西安电子科技大学出版社，1992 周伟.MATLA小波分析高级技术M.西安：西安电子科技大学出版社，2005附加代码：喲像去噪；下面装载原始图像，X中含有被装载的图像load tire;subplot(221); image(X); colormap(map);title(' 原始图像 ');axis square%产生含噪图像init = 2055615866; randn('seed', init); %#ok<*RAND>x = X + 38*randn(size(X);%画出含噪声图像subplot(222); image(X); colormap(map);title(' 含噪声图像 ');axis square;%= %下面进行图像的消噪处理%用小波函数sym4寸x进行2层小波分