小波变换图像去噪MATLAB实现.docx
《小波变换图像去噪MATLAB实现.docx》由会员分享,可在线阅读,更多相关《小波变换图像去噪MATLAB实现.docx(10页珍藏版)》请在冰豆网上搜索。
小波变换图像去噪MATLAB实现
基于小波图像去噪的MATLAB实现
一、论文背景数字图像处理(DigitalImage
Processing,DIP)是指用计算机辅助技术对图像信
号进行处理的过程。
数字图像处理最早出现于
20世
纪50年代,随着过去几十年来计算机、网络技术和通信的快速开展,为信号处理这个学科领域的开展奠定了根底,使得DIP技术成为信息技术中最重要的
学科分支之一。
在现实生活中,DIP应用十分广泛,医疗、艺术、军事、航天等图像处理影响着人类生活
和工作的各个方面。
然而,在图像的采集、获取、
编码和传输的过程中,都存在不同程度被各种噪声所
“污染〞的现象。
如果图像被污染得比拟严重,噪声会
变成可见的颗粒形状,导致图像质量的严重下降。
根
据研究说明,当一张图像信噪比(SNR)低于
14.2dB时,图像分割的误检率就高于0.5%,而参
数估计的误差高于0.6%。
通过一些卓有成效的噪声处理技术后,尽可能地去除图像噪声,我们在从图像
中获取信息时就更容易,有利于进一步的对图像进行
如特征提取、信号检测和图像压缩等处理。
小波变换
处理应用于图像去噪外,在其他图像处理领域都有着
十分广泛的应用。
本论文以小波变换作为分析工具处
理图像噪声,研究数字图像的滤波去噪问题,以提高
图像质量。
二、课题原理1.小波根本原理在
数学上,小波定义为对给定函数局部化的新领域,小
波可由一个定义在有限区域的函数来构
造,称为母小波,〔motherwavelet〕或者叫做基
本
x
x
{x}
小波。
一组小波基函数,,可
以通过缩放和平移根本小波
来生成:
a,b
1xb
1
〔〕
x
a,baa()其中,a为进行缩放的缩放参数,
反映特定基函数的宽度,b为进行平移的j平移参数,指定沿x轴平移的位置。
当a=2和b=ia的情况下,
一维小波基函数序列定义为:
j
j
x22x1
〔2
〕
2i,j
〔〕
x其中,i为平移参数,j为缩放因
子,函数
fx
以小波为基的连续小
〔〕
x
波
变换定义为函数fx和的内积:
a,b
x1b〔3〕Wxf,f(x)()dx
a,ba,baa与时域函数对应,在频域上那么有:
j
可以看出,当|a|
度增大,而且的
〔4〕xaea()a,b
减小时,时域宽度减小,而频域宽
ax,b窗口中心向|ω|增大方向移动。
这说明连续小波的局部是变化的,在高频时分辨率高,在低频时分辨率低,这便是它优于经典傅里叶变换的地方。
总体说来,小波变换具有更好的时频
窗口特性。
2.图像去噪综述所谓噪声,就是指阻碍人的视觉或相关传感器对图像信息进行理解或分析
的各种因素。
通常噪声是不可预测的随机信号。
由于噪声影响图像的输入、采集、处理以及输出的各个环节,尤其是图像输入、采集中的噪声必然影响图像处理全过程乃至最终结果,因此抑制噪声已成
为图像处理中极其重要的一个步骤。
依据噪声对图像的影响,可将噪声分为加性噪声和乘性噪声两大类。
由于乘性噪声可以通过变换当加性噪声来处理,
因此我们一般重点研究加性噪声。
设,,,f(xy)力为理想图像,n(xy)力为噪声,实际输入图像为为g(xy),那么加性噪声可表示为:
,,,,g(xy)=f(xy)+n(xy)
(5〕,其中,n(xy)和图像光强大小无关。
,,图像
去噪的目的就是从所得到的降质图像以g(xy)中尽可能地去除噪声n(x,y),从而复原理想图像f(xy)。
图像去噪就是为了尽量减少图像的均方误差,提高图
像的信噪比,从而尽可能多地保存图像的特征信息。
图像去噪分为时域去噪和频域去噪两种。
传统图像去噪方法如维纳滤波、中值滤波等都属于时域去噪
方法。
而采用傅里叶变换去噪那么属于频域去噪。
这
些方法去噪的依据是一致的,即噪声和有用信号在
频域的不同分布。
我们知道,有用信号主要分布于
图像的低频区域,噪声主要分布在图像的高频区域,
但图像的细节信息也分布在高频区域。
这样在去除
高频区域噪声的同时,难免使图像的一些细节也变
得模糊,这就是图像去噪的一个两难问题。
因此如
何构造一种既能降低图像噪声,又能保存图像细节特征的去噪方法成为图像去噪研究的一个重大课题。
3.小波阈值去噪法3.1小波变换去噪的过程小波去噪是小波变换较为成功的一类应用,其去噪的根本
思路可用框图3-1
来概括,即带噪信号经过预处理,然后利用小波变换把信号分解到各尺度中,在每一尺度下把属于噪声的小波系数去掉,保存并增强属于信号的小波系
数,最后再经过小波逆变换恢复检测信号。
带噪小波分尺
度逆小波恢复图像
分解去噪变换图像
图3-1
小波去噪框图因此,
利用小波变换在去除噪声时,可提取并保存对视觉起主要作用的边缘信息。
而传统的傅立叶变换去噪方法在去除噪声和边沿保持上存在着矛盾,原因是傅立叶变换方法在时域不能局部化,难以检测到局域突变信号,在去除噪声的同时,也损失了图像边
沿信息。
由此可见,与傅立叶变换去噪方法相比,
小波变换去噪方法具有明显的优越性。
3.2小波阈值去噪的根本方法3.2.1阈值去噪原理Donoho提出的小波阈值去噪方法的根本思想是当w小于某个临界
阈值jk,时,认为这时的小波系数主要是由噪声引起的,予以舍弃。
当w大于这个临界jk,阈值时,认为这时的小波系数主要是由信号引起,那么就把这一
局部的w直接jk,保存下来(硬阈值方法),或者按某一个固定量向零收缩(软阈值方法),然后用新的小波系数进行小波重构得到去噪后的信号。
此方法可通
过以下三个步骤实现:
(1)先对含噪声信号f(t)做小波变换,得到一组小波分解系数w。
jk,
(2)通过对分解得到的小波系数w得出估计小波系数使进行阈值处理,wjkj,k,得w-u,尽可能的小。
jkjk,,(3)利用估计小波系数进行小波重构,得到估计信号了,即为
去噪f(t)wj,k之后的信号。
需要说明的是,在小波阈值去噪法中,最重要的是闭值函数和闲值的选取。
3.2.2阈值函数的选取阈值函数关系着重构信号的连续性和精度,对小波去噪的效果有很大影响。
目前,
阈值的选择主要分硬阈值和软阈值两种处理方式。
其中,软阈值处理是将信号的绝对值与阈值进行比
较,当数据的绝对值小于或等于阈值时,令其为零;
大于阈值的数据点那么向零收缩,变为该点值与阈值
之差。
而硬阈值处理是将信号的绝对值阈值进行比
较,小于或等于阈值的点变为零,大于阈值的点不
变。
但硬阈值函数的不连续性使消噪后的信号仍然
含有明显的噪声;采用软阈值方法虽然连续性好,但估计小波系数与含噪信号的小波系数之间存在恒定的偏差,当噪声
信号很不规那么时显得过于光滑。
4、基于小波变换的图像分解与重构,二维离散小波主要解决二维多分辨率分析问题,如一幅二维离散图像{c(mn)},二小波可以将它分解为各层各个分辨率上的近似分量
cAj,水平方向细节分量cHj,垂直方向细节分量cVj,对角线方向细节分量cDj,其二层小波图像分解过程
如图4-1所示:
图4-1小波图像分解过程
图4-2小波图像分解过程其二层小波图像重构过程正好与此相反如图4-2所示,基于小波变换的图像处理,是通过对图像分解过程中所产生的近似分量与细节分量系数的调整,使重构图像满足特定条件,而实现图像处理。
三、程序实现图像消噪常用的图像去噪方法是小波阈值去噪法,它是一种实现简单而效果较好的去噪方法,阈值去噪方法的思想很简单,就是对小波分解后的各层稀疏模大于和小于某阈值的系数
分别进行处理,然后利用处理后的小波系数重构出去噪后的图像。
在阈值去噪中,阈值函数表达了对小波分解稀疏的不同处理策略以及不同的估计方法,常用的阈值函数有硬阈值和软阈值函数,硬阈值函数可以很好的保存图像边缘等局部特征,但图像会出现伪吉布斯效应,等视觉失真现象,而软阈值处理相对较平稳,但可能会出现边缘模糊等失真现象,为此人们又提出了半软阈值函数。
小波阈值去噪方法处理阈值的选取,另一个关键因素是阈值的具体估计,如果阈值
太小,去噪后的图像仍然存在噪声,相反如果阈值太大,重要图像特
征又将被滤掉,引起偏差。
从直观上讲,对给定的小
波系数,噪声越大,阈值就越大。
图像信号的小波去噪步骤与一维信号的去噪步骤完全相同,只使用二维小波2分析工具代替了一维小波分析工具,如果用固定阈值形式,那么选择的阈值用m代替了一维信号中的n。
这三步是:
1〕二维信号的小波分解。
选择一个小波和小波分解的层次N,然后计算信号S到第
N层的分解。
2〕对高频系数进行阈值量化,对于从一到N的每一层,选择一个阈值,并对这一层的高频系数进行软阈值化处理。
3〕二维小波的重构,根据小波分解的第N层的低频系数和经过修改的从第
一层到第N层的高频系数,来计算二维信号的小波
重构。
下面就通过具体实例来说明利用小波分析进行
图像去噪的问题。
对给定图像进行去噪的二维小波
去噪源程序:
clear;
%清理工作空间
load
wbarb;
%装载原始图像
subplot(221);
%
新建窗口
image(X);
%显示图像
colormap(map);
%设置色彩索引图title('原始图像
');
%设置图像标题axissquare;
%设置显示
比例,生成含噪图像并图示
init=2055615866;
%初
始值
randn('seed',init);
%随机值
XX=X+8*randn(size(X));
%
添加随机噪声
subplot(222);
%新建窗口
image(XX);
%
显示图像colormap(map);
%
设置色彩索引图
title('
含噪图像');
%
设置图像标题
axis
square;
%用小波函数coif2对图像XX进行2层
[c,l]=wavedec2(XX,2,'coif2');%分解
n=[1,2];
%
设置尺度向量p=[10.28,24.08];%设置阈值向量,对
高频小波系数进行
阈%nc=wthcoef2('h',c,l,n,p,'s');%nc=wthcoef2('v',c,l,
n,p,'s');
重构
X1=waverec2(nc,l,'coif2');
subplot(223);
%图像的二维小波%新建窗口
image(X1);
%显示图像
colormap(map);%设置色彩索引图title('第
一次消噪后的图像');%设置图像标题axissquare;%
设置显示比例,再次对高频小波系数进行阈值处
理%mc=wthcoef2('h',nc,l,n,p,'s');mc=wthcoef2('v',nc,l,
n,p,'s');%mc=wthcoef2('d',nc,l,n,p,'s');
X2=waverec2(mc,l,'coif2');%图像的二维小波重构
subplot(224);
%新建窗口image(X2);
%
显示图像colormap(map);
%设置色彩索引图
title('第二次消噪后的图像');
%设置图像标题
axis
square;
%设置显示比例程序运行结果:
原始图像含噪
图像
第一次消噪
后的图像第二次消噪后的图像
图5-1去噪前后图像比拟上图中
几幅图像,可见第一次去早滤除了大局部的高频噪,
但与原图比拟,依然有不少的高频噪声,第二次去噪
在第一次的去噪根底上,再次滤除高频噪声,去噪效
果较好,但图像的质量比原图稍差。
六、总结随着信息时代计算机的日益普及,人们对
数字图像的质量要求越来越高。
但是数字图像在采集
和传输过程中,难免会受到噪声的污染,这不仅不符合人们的视觉效果,而且也不利于图像的进一步处理。
因此,图像去噪具有很强的理论意义和应用价值。
图像消噪是信号处理中的一个经典问题,传统的消噪方
法多采用平均或线性方法进行,但是其消噪效果不
好,随着小波理论的不断完善,它以自身良好的时频
特性在图像消噪领域受到越来越多的关注,文中将以
MATLAB为平台介绍以小波变换去除图像噪声的根本
方法。
采用基于小波变换的图像增强技术可以通
过对低频分解系数进行增强处理,对高频分解系数进
行衰减处理到达图像增强的作用。
小波阈值去噪方法
是小波去噪领域使用较多的方法,因为其理论相比照
较成熟,而且去噪效果也比拟好。
它是集图像去噪和增强为一体的优秀处理方法,被广泛的应用于图像处理中,大量实践也证明该算法优于其他增强技术。