小波阈值图像降噪及MATLAB仿真.docVIP

小波阈值图像降噪及MATLAB仿真

PAGE

PAGE9

小波阈值图像降噪及MATLAB仿真

图像信号在生成和传输过程中常常因受到各种噪声的干扰和影响而使图像降质，这对后续图像的处理(如分割、压缩和图像理解等)将产生不利影响。为了抑制噪声，改善图像质量，便于更高层次的处理，必须对图像进行降噪预处理。

小波降噪的方法有多种，如利用小波分解与重构的方法滤波降噪、利用小波变换模极大值的方法去噪、利用信号小波变换后空域相关性进行信噪分离、非线性小波阈值方法去噪、平移不变量小波降噪法，以及多小波降噪等等。归结起来主要有三类:模极大值检测法、阈值降噪法和屏蔽(相关)降噪法。其中最常用的就是阈值法去噪，其基本思想就是利用图像小波分解后，各个子带图像的不同特性选取不同的阈值，从而达到较好的降噪目的。而且，小波变换本身是一种线形变换，因而对于类似于高斯噪声的效果较好。线性运算往往还会造成边缘模糊，小波分析技术正因其独特的时频局部化特性在图像信号和噪声信号的区分以及有效去除噪声并保留有用信息等方面较之传统的降噪具有明显的优势，且在降噪的同时实现了图像一定程度的压缩和边缘特征的提取。所以小波降噪具有无可比拟的优越性。小波降噪主要优点有：

低熵性，小波系数的稀疏分布，使得图象变换后的熵降低；

多分辨率，由于采用了多分辨率的方法，所以可以非常好地刻画信号的非平稳特征，如边缘、尖峰、断点等；

去相关性，因为小波变换可以对信号进行去相关，且噪声在变换后有白化趋势，所以小波频域比时域更利于降噪；

选基灵活性，由于小波变换可以灵活选择变换基，从而对不同应用场合、不同的研究对象，可以选用不同的小波函数，以获得最佳的效果。

一、阈值去噪法简述

1992年，斯坦福大学的DonohoDL和Johnstone教授提出一种具有良好的统计优化特性的降噪方法，称作“WaveletShrinkage”（即阈值收缩法）。

该方法的主要思想是：基于图像和噪声在经小波变换后具有不同的统计特性：图像本身的能量对应着幅值较大的小波系数,主要集中在高频（）；噪声能量则对应着幅值较小的小波系数，并分散在小波变换后的所有系数中。根据该特征，设置一个阈值门限，认为大于该阈值的小波系数的主要成份为有用的信号，给予收缩后保留；小于该阈值的小波系数，主要成份为噪声，予以剔除，这样就可以达到降噪的目的。

降噪时，通常认为低通系数含有大量的图像能量，一般不作处理，只对剩余三个高通部分进行处理。因此，一次阈值降噪并不能完全去除噪声，还需要对未作处理的低频部分（）再次进行小波分解和阈值去噪，直到实际图像与估计图像的偏差达到最小值。

但是，随着分解和降噪次数的增加，小波系数中的噪声能量越来越少，并且趋于分散，降噪的效果将逐渐降低。一般来说，进行3-4层小波分解和降噪就可以达到满意的降噪效果。

小波阈值降噪的基本思路是：

（1）先对含噪信号做小波变换，得到一组小波系数；

（2）通过对进行阈值处理，得到估计系数，使得与两者的差值尽可能小；

（3）利用进行小波重构，得到估计信号即为降噪后的信号。

Donoho提出了一种非常简洁的方法对小波系数进行估计。对连续做几次小波分解后，有空间分布不均匀信号各尺度上小波系数在某些特定位置有较大的值，这些点对应于原始信号的奇变位置和重要信息，而其他大部分位置的较小；对于白噪声，它对应的小波系数在每个尺度上的分不都是均匀的，并随尺度的增加，系数的幅值减小。因此，通常的降噪办法是寻找一个合适的数作为阈值（门限），把低于λ的小波函数（主要由信号引起），设为零，而对于高于的小波函数（主要由信号引起），则予以保留或进行收缩，从而得到估计小波系数，它可理解为基本由信号引起的，然后对进行重构，就可以重构原始信号。

估计小波系数的方法如下，取：

其语法格式为：

[THR,SORH,KEEPAPP,CRIT]=ddencmp（IN1，IN2，X）

[THR,SORH,KEEPAPP]=ddencmp（IN1，‘wv’，X）

[THR,SORH,KEEPAPP,CRIT]=ddencmp（IN1，‘wp’，X）

2阈值降噪

MATLAB中实现阈值降噪的函数有wden、wdencmp、wpdencmp、wthresh、wpthcoef和wthcoef2。这里主要介绍函数wdencmp。

其语法格式为：

[XC,CXC,LXC,PERF0,PERFL2]=wdencmp(gbl,X,wname,N,THR,SORH,KEEPAPP)

[XC,CXC,LXC,PERF0,PERFL2]=wdencmp(lvd,X,wname,N,THR,SORH)

[XC,CXC,LXC,PERF0,PERFL2]=wdencmp(lvd,C,L,wname,N,THR,SORH)

函