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

基于小波变换的图像去噪方法的研究

（附送程序，见上传者“我的文档”）

摘要

图像在生成和传输的过程中会受到各种噪声的干扰,对信息的处理、传输和存储造成极大的影响。寻求一种既能有效地减小噪声,又能很好地保留图像边缘信息的方法,是人们一直追求的目标。小波分析是局部化时频分析,它用时域和频域联合表示信号的特征,是分析非平稳信号的有力工具。它通过伸缩、平移等运算功能对信号进行多尺度细化分析,能有效地从信号中提取信息。

本文对基于小波变换的图像去噪方法进行了研究分析, 详细介绍了其去噪原理和算法,分析了去噪过程中参数的选取问题,并给出了一些选取依据;详细介绍了小波系数相关性去噪方法的原理和算法。最后对均值滤波、中值滤波和维纳滤波方法在高斯噪声下进行了分析比较,并给出了仿真实验结果，结果证明小波去噪十分有效，其结果好于其它3种滤波。

关键词：小波变换，图像去噪，阈值，阈值函数

1．引言

数字图像在我们日常生活中起着非常重要的作用，它与我们的日常生活息息相关，例如在卫星、电视、核磁共振、计算机视觉、地球信息系统以及天文学中应用非常广泛。但是一般情况下采集到的数字图像是含有噪声的。噪声[1]可以理解为“妨碍人们感觉器官对所接收的信源信息理解的因素”。图像在生成和传输的过程中灰受到各种噪声的干扰，对信号的处理、传输和存储造成极大的影响。数字图像之所以含有噪声这是因为在图像的采集、获取、编码和传输的过程中，所有的图像均不同程度地被可见或不可见的噪声“污染”。对于这种“污染”，如果信噪比(SNR)低于一定水平，就会影响图像场景内容的表示，直接导致图像质量的下降。除了视觉质量上下降外，噪声还可能掩盖一些重要的图像细节，使图像的熵增大，从而对于图像数据的有效压缩起到了一定的妨碍作用。对于图像在采集、获取过程造成的“污染”，我们虽然尽量提高硬件设备以获取质量更高的图像，但图像传感器的截止频率总是有一定的，受硬件水平和价格的限制，且图像在编码和传输过程中造成的“污染”，必需采取有效的降噪技术才能提高图像的质量。

对图像进行去噪最初主要是在空域内进行的，图像空域去噪方法很多，主要是通过各种滤波器对图像进行去噪。例如均值滤波器、顺序统计滤波器、维纳滤波器等。为了进一步提高去噪的效果，在变换域中进行降噪处理成为有效的方法，图像变换域去噪就是对图像进行某一种变换，然后将图像从时域变换到变换域中，再对变换域中的图像变换系数按照某种方法进行处理，最后再对处理后的系数按照某种方法进行反变换，这样就实现了将图像去除图像噪声的目的。将图像从时域转换到变换域的变换方法很多，例如傅立叶变换、小波变换等等。

小波变换是在短时傅立叶变换的基础上发展起来的一种新型的变换方法。小波变换具有多分辨率分析的特点，在时域、频域都具有较强的表征信号局部特征的能力，因此基于小波分析的图像去噪技术已成为图像去噪的一个重要方法。本文采用小波阈值去噪的方法，从去噪的效果上比较了多种去噪方法的优劣，实验证明小波去噪在图像噪声处理中起到很好的效果。

2．小波变换概述

2.1 小波变化去噪技术研究现状

上个世纪八十年代Mallet 提出了 MRA(Multi_Resolution Analysis)，并首先把小波理论运用于信号和图像的分解与重构，利用小波变换模极大值原理进行信号的奇异性检测，提出了交替投影算法用于信号重构，为小波变换用于图像处理奠定了基础[1]。后来，人们根据信号与噪声在小波变换下模极大值在各尺度上的不同传播特性，提出了基于模极大值去噪的基本思想。1992年，Donoho和Johnstone提出了“小波收缩”，它较传统的去噪方法效率更高。“小波收缩”被 Donoho和Johnstone证明是在极小化极大风险中最优的去噪方法，但在这种方法中最重要的就是确定阈值。1995年，Stanford大学的学者D.L.Donoho和I.M.Johnstone提出了通过对小波系数进行非线性阈值处理来降低信号中的噪声[2]。从这之后的小波去噪方法也就转移到从阈值函数的选择或最优小波基的选择出发来提高去噪的效果。影响比较大的方法有以下这么几种：Eero P.Semoncelli和Edward H.Adelson提出的基于最大后验概率的贝叶斯估计准则确定小波阈值的方法[3]；Elwood T.Olsen等在处理断层图像时提出了三种基于小波相位的去噪方法：边缘跟踪法、局部相位方差阈值法以及尺度相位变动阈值法；学者Kozaitis结合小波变换和高阶统计量的特点提出了基于高阶统计量的小波阈值去噪方法[4]；G.P.Nason等利用原图像和小波变换域中图像的相关性用GCV(general cross- validation)法对图像进行去噪；Hang.X和Woolsey等人提出结合维纳滤波器和小波阈值的方法对信号进行去噪处理[5]，Vasily Strela等人将一类新的特性良好的小波(约束对)应用于图像去噪的方法[6]；同时，在19世纪60年代发展的隐马尔科夫模型(Hidden Markov Model)，是通过对小波系数建立模型以得到不同的系数处理方法；后又有人提出了双变量模型方法[7]，它是利用观察相邻尺度间父系数与子系数的统计联合分布来选择一种与之匹配的二维概率密度函数。这些方法均取得了较好的效果，对小波去噪的理论和应用奠定了一定的基础。

小波具有低墒性、多分辨率、去相关性、选基灵活性等特点。

2.2 连续小波变换

1. 小波基函数

所谓小波(Wavelet)，即存在于一个较小区域的波。小波函数的数学定义是：

??w?满足： 设??t?为一平方可积函数，即??t??L2?R?，若其傅立叶变换?C?????w?w2Rdw??（2.1）

时，则称??t?为一个基本小波或小波母函数，并称上式是小波函数的可容许条件。

根据小波函数的定义，小波函数一般在时域具有紧支集或近似紧支集，即函数的非零值定义域具有有限的范围，这即所谓“小”的特点;另一方面，根据可容许性条件可知??w?w?0?0，即直流分量为零，因此小波又具有正负交替的波动性。

将小波母函数??t?进行伸缩和平移，设其伸缩因子(亦称尺度因子)为a，平移因子为b，并记平移伸缩后的函数为?a,b?t?，则:

???a,b?t??a??t?a,a,b?R;a?0（2.2）

?12并称?a,b?t?为参数a和b小波基函数。由于a和b均取连续变换的值，因此又称为连续小波基函数，它们是由同一母函数??t?经伸缩和平移后得到的一组函数系列。

2. 连续小波变换

将L2?R?空间的任意函数f?t?在小波基下进行展开，称其为函数f?t?的连续小波变换CWT，变换式为：

WTf?a,b???f,?a,b??1a?f?t?????dt（2.3）

Rt?ba当小波的容许性条件成立时，其逆变换为：

f?t??其中C?????w?w21C????da2??a?????b?WTf?a,b????t?adb（2.4）

Rdw??为??t?的容许性条件。

另外，在小波变换过程中必须保持能量成比例，即:

?Rdaa2?WTf?a,b?db?C??f?x?dx（2.5）

2RR2由CWT的定义可知，小波变换和傅立叶变换一样，也是一种积分变换，其中

可见小波变换对函数f?t?在小波基上的展开具有多分WTf?a,b?为小波变换系数。

辨率的特性，这种特性正是通过缩放因子a和平移因子b来得到的。

一个一维函数f?t?的连续小波变换是一双变量的函数，变量比f?t?多一个，因此称连续小波变换是超完备的，因为它要求的存储量和它代表的信息量都显著增加了。对于变量超过一个的函数来说，这个变换的维数也将增加。

若f?t?是一个二维函数，则它的连续小波变换是：

WTf?a,bx,by??1a??????????f?t1,t2???x?bxa,y?bya?dxdy（2.6）

y其中，bx，by表示在两个维度上的平移，二维连续小波逆变换为：

f?x,y??1C????0????daa3??????WTf?a,bx,by???x?bxa,y?bya?dbdb（2.7）

x同样的方法可以推广到两个或两个以上的变量函数上。

2.3 离散小波变换

计算机中的图像信息是以离散信号形式存放的，所以需要将连续小波变换离散化。而最基本的离散化方法就是二进制离散，一般将这种经过离散化的小波及其变换叫做二进小波和二进变换。需要注意的是这里的离散化都是针对连续的尺度因子a和连续平移因子b的，而不是针对时间t的。这儿限制尺度因子a总是正数。

1. 尺度与位移的离散化

对连续小波基函数?a,b?t?尺度因子a和平移因子b进行离散化可以得到离散小波变换WTf?a,b?，从而减少小波变换系数的冗余度。在离散化时通常对尺度因

mm子a和平移因子b按幂级数进行离散化，即取a?a0（m为整数，a0?1,但,b?b0一般都假定a0?1），得到离散小波函数为：

?m,n?t??其对应系数为：

1a0??mt?na0b0ma0??1a0?m??a0t?nb0?（2.8）

Cm,n??f?t?,?m,n???f?t??m,n?t?dt（2.9）

????2. 二进制小波变换

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库基于小波变换的图像去噪方法研究报告附MATLAB程序在线全文阅读。