基于小波变换的图像去噪方法研究报告附MATLAB程序

基于小波变换的图像去噪方法的研究（附送程序，见上传者“我的文档”）摘要图像在生成和传输的过程中会受到各种噪声的干扰,对信息的处理、传输和存储造成极大的影响。寻求一种既能有效地减小噪声,又能很好地保留图像边缘信息的方法,是人们一直追求的目标。小波分析是局部化时频分析,它用时域和频域联合表示信号的特征,是分析非平稳信号的有力工具。它通过伸缩、平移等运算功能对信号进行多尺度细化分析,能有效地从信号中提取信息。本文对基于小波变换的图像去噪方法进行了研究分析,详细介绍了其去噪原理和算法,分析了去噪过程中参数的选取问题,并给出了一些选取依据;详细介绍了小波系数相关性去噪方法的原理和算法。最后对均值滤波、中值滤波和维纳滤波方法在高斯噪声下进行了分析比较,并给出了仿真实验结果，结果证明小波去噪十分有效，其结果好于其它3种滤波。关键词：小波变换，图像去噪，阈值，阈值函数1．引言数字图像在我们日常生活中起着非常重要的作用，它与我们的日常生活息息相关，例如在卫星、电视、核磁共振、计算机视觉、地球信息系统以及天文学中应用非常广泛。但是一般情况下采集到的数字图像是含有噪声的。噪声[1]可以理解为“妨碍人们感觉器官对所接收的信源信息理解的因素”。图像在生成和传输的过程中灰受到各种噪声的干扰，对信号的处理、传输和存储造成极大的影响。数字图像之所以含有噪声这是因为在图像的采集、获取、编码和传输的过程中，所有的图像均不同程度地被可见或不可见的噪声“污染”。对于这种“污染”，如果信噪比(SNR)低于一定水平，就会影响图像场景内容的表示，直接导致图像质量的下降。除了视觉质量上下降外，噪声还可能掩盖一些重要的图像细节，使图像的熵增大，从而对于图像数据的有效压缩起到了一定的妨碍作用。对于图像在采集、获取过程造成的“污染”，我们虽然尽量提高硬件设备以获取质量更高的图像，但图像传感器的截止频率总是有一定的，受硬件水平和价格的限制，且图像在编码和传输过程中造成的“污染”，必需采取有效的降噪技术才能提高图像的质量。对图像进行去噪最初主要是在空域内进行的，图像空域去噪方法很多，主要是通过各种滤波器对图像进行去噪。例如均值滤波器、顺序统计滤波器、维纳滤波器等。为了进一步提高去噪的效果，在变换域中进行降噪处理成为有效的方法，图像变换域去噪就是对图像进行某一种变换，然后将图像从时域变换到变换域中，再对变换域中的图像变换系数按照某种方法进行处理，最后再对处理后的系数按照某种方法进行反变换，这样就实现了将图像去除图像噪声的目的。将图像从时域转换到变换域的变换方法很多，例如傅立叶变换、小波变换等等。小波变换是在短时傅立叶变换的基础上发展起来的一种新型的变换方法。小波变换具有多分辨率分析的特点，在时域、频域都具有较强的表征信号局部特征的能力，因此基于小波分析的图像去噪技术已成为图像去噪的一个重要方法。本文采用小波阈值去噪的方法，从去噪的效果上比较了多种去噪方法的优劣，实验证明小波去噪在图像噪声处理中起到很好的效果。2．小波变换概述2.1小波变化去噪技术研究现状上个世纪八十年代Mallet提出了MRA(Multi_ResolutionAnalysis)，并首先把小波理论运用于信号和图像的分解与重构，利用小波变换模极大值原理进行信号的奇异性检测，提出了交替投影算法用于信号重构，为小波变换用于图像处理奠定了基础[1]。后来，人们根据信号与噪声在小波变换下模极大值在各尺度上的不同传播特性，提出了基于模极大值去噪的基本思想。1992年，Donoho和Johnstone提出了“小波收缩”，它较传统的去噪方法效率更高。“小波收缩”被Donoho和Johnstone证明是在极小化极大风险中最优的去噪方法，但在这种方法中最重要的就是确定阈值。1995年，Stanford大学的学者D.L.Donoho和I.M.Johnstone提出了通过对小波系数进行非线性阈值处理来降低信号中的噪声[2]。从这之后的小波去噪方法也就转移到从阈值函数的选择或最优小波基的选择出发来提高去噪的效果。影响比较大的方法有以下这么几种：EeroP.Semoncelli和EdwardH.Adelson提出的基于最大后验概率的贝叶斯估计准则确定小波阈值的方法[3]；ElwoodT.Olsen等在处理断层图像时提出了三种基于小波相位的去噪方法：边缘跟踪法、局部相位方差阈值法以及尺度相位变动阈值法；学者Kozaitis结合小波变换和高阶统计量的特点提出了基于高阶统计量的小波阈值去噪方法[4]；G.P.Nason等利用原图像和小波变换域中图像的相关性用GCV(generalcross-validation)法对图像进行去噪；Hang.X和Woolsey等人提出结合维纳滤波器和小波阈值的方法对信号进行去噪处理[5]，VasilyStrela等人将一类新的特性良好的小波(约束对)应用于图像去噪的方法[6]；同时，在19世纪60年代发展的隐马尔科夫模型(HiddenMarkovModel)，是通过对小波系数建立模型以得到不同的系数处理方法；后又有人提出了双变量模型方法[7]，它是利用观察相邻尺度间父系数与子系数的统计联合分布来选择一种与之匹配的二维概率密度函数。这些方法均取得了较好的效果，对小波去噪的理论和应用奠定了一定的基础。小波具有低墒性、多分辨率、去相关性、选基灵活性等特点。2.2连续小波变换1.小波基函数所谓小波(Wavelet)，即存在于一个较小区域的波。小波函数的数学定义是：设t为一平方可积函数，即RLt2，若其傅立叶变换wˆ满足：dwCRww2（2.1）时，则称t为一个基本小波或小波母函数，并称上式是小波函数的可容许条件。根据小波函数的定义，小波函数一般在时域具有紧支集或近似紧支集，即函数的非零值定义域具有有限的范围，这即所谓“小”的特点;另一方面，根据可容许性条件可知00ww，即直流分量为零，因此小波又具有正负交替的波动性。将小波母函数t进行伸缩和平移，设其伸缩因子(亦称尺度因子)为a，平移因子为b，并记平移伸缩后的函数为tba,，则:0;,,21,aRbaatatba（2.2）并称tba,为参数a和b小波基函数。由于a和b均取连续变换的值，因此又称为连续小波基函数，它们是由同一母函数t经伸缩和平移后得到的一组函数系列。2.连续小波变换将RL2空间的任意函数tf在小波基下进行展开，称其为函数tf的连续小波变换CWT，变换式为：dttffbaWTabtRabaf1,,,（2.3）当小波的容许性条件成立时，其逆变换为：dbbaWTtfabtfadaC,21（2.4）其中dwCRww2为t的容许性条件。另外，在小波变换过程中必须保持能量成比例，即:RRfRadadxxfCdbbaWT22,2（2.5）由CWT的定义可知，小波变换和傅立叶变换一样，也是一种积分变换，其中baWTf,为小波变换系数。可见小波变换对函数tf在小波基上的展开具有多分辨率的特性，这种特性正是通过缩放因子a和平移因子b来得到的。一个一维函数tf的连续小波变换是一双变量的函数，变量比tf多一个，因此称连续小波变换是超完备的，因为它要求的存储量和它代表的信息量都显著增加了。对于变量超过一个的函数来说，这个变换的维数也将增加。若tf是一个二维函数，则它的连续小波变换是：dxdyttfbbaWTabyabxayxfyx,,,,211（2.6）其中，xb，yb表示在两个维度上的平移，二维连续小波逆变换为：yxabyabxyxfadaCdbdbbbaWTyxfyx,,,,013（2.7）同样的方法可以推广到两个或两个以上的变量函数上。2.3离散小波变换计算机中的图像信息是以离散信号形式存放的，所以需要将连续小波变换离散化。而最基本的离散化方法就是二进制离散，一般将这种经过离散化的小波及其变换叫做二进小波和二进变换。需要注意的是这里的离散化都是针对连续的尺度因子a和连续平移因子b的，而不是针对时间t的。这儿限制尺度因子a总是正数。1.尺度与位移的离散化对连续小波基函数tba,尺度因子a和平移因子b进行离散化可以得到离散小波变换baWTf,，从而减少小波变换系数的冗余度。在离散化时通常对尺度因子a和平移因子b按幂级数进行离散化，即取mmbbaa00,（m为整数，,10a但一般都假定10a），得到离散小波函数为：0011,00000nbtatmaabnatanmmm（2.8）其对应系数为：dtttftfCnmnmnm,,,,（2.9）2.二进制小波变换二进小波变换是一种特殊的离散小波变换，特别地令参数20a，10b，则有ntmnmm222,。该二进尺度分解的原理在二十世纪三十年代由Littlewood和Paley在数学上进行了研究证明。离散小波变换为：dtttfnmnmWTnmf,,,（2.10）离散二进小波变换为：dtttfnmnmWTnmf,,,（2.11）二维离散小波变换：我们考虑二维尺度函数是可分离的情况，也就是：2121,xxxx（2.12）设ix是与ix对应的一维小波函数，则有：21211,xxxx（2.13）21212,xxxx（2.14）21213,xxxx（2.15）以上三式就建立了二维小波变换的基础。3.图像的小波变换及Mallat算法图像是二维信号，二维多分辨率分析与一维类似，但空间)(2RL变成)*(2RRL，一维中引入的尺度函数)(x变为),(yx。设ZjjV是)(2RL的一个多分辨率分析，则可以证明，张量空间ZjjV2：jjjVVV2构成)*(2RRL的一个多分辨率分析，并且二维多分辨率分析ZjjV2的二维尺度函数),(yx为)()(),(yxyx(2.16)式中：)(x是ZjjV尺度函数(一维)。式(2.28)说明了二维尺度函数的可分离性。对于每一个Zj,函数系}),()()(),({2,,,,Zmnyxyxmjnjmnj构成ZjjV2的规范正交基，这里下标j，n，m的含义是：)0(),2()2(2),(,,jmynxyxjjjmnj（2.17）我们将ZjjV2称为)*(2RRL的可分离多分辨率分析。因)(x、)(y都是低通的尺度函数，所以ZjjV2平滑的低通空间。如果)(x是一维多分辨率分析ZjjV的正交小波基，则二维多分辨率分析的三个小波函数为：)()(),()()(),()()(),(321yxyxyxyxyxyx（2.18）对于每一个Zj，它们的整数平移系为：)()(),()()(),()()(),(,,,,3,,,,2,,,,1yxyxyxyxyxyxmjnjmnjmjnjmnjmjnjmnj（2.19）注意这里的上标只是索引而不是指数。它们构成了ZjjW}{2的规范正交基。因此以上的三个正交基中都至少包含一个带通的)(x或)(y，所以它们都是带通的。也就是说这三部分反应的都是细节信息。具体来说，函数系}2,2(2{)},({,,mynxyxjjjmnj，3,2,1,0j是)*(2RRL的正交归一基，其中均为整数，=1,2,3分别对应于水平、垂直和对角三个方向。对于任一二维图像信号)*(),(2RRLyxf，在分辨率j2下有：2,,332,,222,,112,,),(,),