第 19卷第 2期四川理工学院学报 (自然科学版 V ol . 19 No. 2JOURNAL OF SICHUAN UNIVERSITY OF2006年 4月 SCIENCE & ENGINEERING (NATURAL SCIENCE EDITION Apr . 2006文章编号:1673-1549(2006 02-0008-04基于小波变换与阀值收缩法的图像增强去噪高飞,杨平先,孙兴波(四川理工学院电子与信息工程系,四川自贡 643000摘要:提出了一种基于小波变换与阀值收缩法的图像增强去噪方法。
图像经过小波分解后可以得到一系列不同尺度上的子带图像, 在不同尺度的子带图像上进行基于阈值收缩滤波的细节系数增强, 再进行小波重构,即可得到增强后的图像。
该方法可以有效地去除噪声,增强图像的平均梯度,改善图像的视觉效果。
关键词:图像增强;小波变换;去噪;阀值收缩中图分类号:TP391 文献标识码:A前言小波变换是传统傅里叶变换的继承和发展, 由于小波的多分辨率分析具有良好的空间域和频率域局部化特性, 对高频采用逐渐精细的时域或空域步长, 可以聚焦到分析对象的任意细节, 因此特别适合于图像信号这一类非平稳信源的处理,已成为一种信号/ 图像处理的新手段。
目前,小波分析已被成功地应用于信号处理、图象处理、语音与图像编码、语音识别与合成、多尺度边缘提取和重建、分形及数字电视等科学领域 [1]。
图像增强是图像处理中一个非常重要的研究领域,已经有许多非常成熟和有效的方法如直方图均衡、高通滤波、反掩模锐化法等,但是这些传统的图像增强方法都存在着不足,如噪声放大、有时可能引入新的噪声结构等。
目前已经有许多关于小波变换在图像处理方面的应用研究, 取得了非常不错的效果。
针对传统图像增强中存在的一些问题,如增强噪声、丢失细节等,本文提出了一种基于阈值收缩法[2]的小波图像增强方法, 实验结果表明该方法能较好地解决图像增强中的噪声放大的问题, 并能非线性地增强图像的细节信息,保持图像的边缘特征,改善图像的视觉效果,是一种很有效的方法。
1 小波变换小波变换的基本思想是用一族函数去表示或逼近一信号, 这一族函数称为小波函数系。
它是通过一小波母函数的伸缩和平移产生其“子波”来构成的,用其变换系数描述原来的信号 [3]。
设相应的尺度函数为 (x ϕ,小波函数为(x ψ,二维尺度函数 , (y x ϕ,是可分离的,即: ( ( , (y x y xϕϕϕ=,即可以构造 3个二维基本小波函数:( ( , (1y x y x ψϕψ=, ( ( , (2y x y x ϕψψ=, ( ( , (3y x y x ψψψ=那么,二维小波基可以通过以下伸缩平移实现:2, 2(2 , (, , n y m x y x j j i j in m j −−=−−−ψψ 3, 2, 1, , , =∈i Z n m j这样,一个二维图像信号 , (y x f 在尺度 j2下的平滑成分(低频分量可用二维序列 , (n m D j 表示为: , ( , ( , (, , y x y x f n m D n m j j ϕ=细节成分可以表示为: , ( , ( , (1, , 1y x y x f n m C n m j j ψ= , (, ( , (2, , 2y x y x f n m C n m j j ψ= , (, ( , (3, , 3y x y x f n m C n m j j ψ=收稿日期:2005-10-31基金项目:四川省教育厅青年基金项目资助(2004B018作者简介:高飞(1978- ,男,河南南阳人,2004级硕士研究生,研究方向为数据融合、图像处理及模式识别。
第 19卷第2期高飞等:基于小波变换与阀值收缩法的图像增强去噪 9二维小波变换的重构公式可用下式表示:((( , ( , (, , , 1, , , , y x n m D y x n m C y x f J j i m n n m j j i n m j m n i j∑∑∑∑∑∑=+=如果引入增强系数 i j W , ,则重构公式为:((( , ( , (, , , 1, , , , , y x n m D y x n m C W y x f J j i m n n m j j i n m j m n i j i j ∑∑∑∑∑∑=+=一幅 N N ×的图像经过一层二维 Mallat 快速算法实现的小波分解后,将得到四个大小均为 22N N ×的子带图像, 即一个逼近信号(对应着 jD ,水平和垂直方向均为低频分量和 3个细节信号 (对应着 321, , j j j C C C ,水平和垂直方向至少有一个高频分量 , 从上面的分析可以知道, 图像信号的小波分解实质上就是把图像信号分解成不同频带范围内的图像分量,每一层小波分解都将待分解图象分解成四个子带图像:LL (水平与垂直方向皆为低频成分、 LH (水平低频、垂直高频、 HL (水平高频、垂直低频、 HH (水平与垂直方向皆为高频成分 ,因此,可以采用不同的方法来增强不同频率范围内图像的细节分量,突出不同尺度的细节,从而改善图像的视觉效果,如本文采用的阀值收缩法。
在实际应用中,可以根据噪声水平和感兴趣的细节所处的尺度,选用不同的阀值和增强系数对分解后的图像进行重构。
2 小波图像增强去噪方法研究2.1 算法思路现实中采集来的图像往往含有大量噪声, 这些噪声主要分布在小波变换域的小尺度小波系数上, 而这些细节系数也包含了大量的图像细节信息,传统的方法对噪声考虑不足,只是简单的增强细节信号, 存在着噪声放大的问题, 因此这里提出利用小波阀值收缩法去噪。
它的主要理论依据是, 小波变换特别是正交小波变换具有很强的去数据相关性,它能够使图像的能量在小波域集中在一些大的小波系数中; 而噪声的能量却分布于整个小波域内, 因此, 经小波分解后, 图像的小波系数幅值要大于噪声的系数幅值,可以认为幅值较大的小波系数一般以图像信号为主,而幅值较小的系数在很大程度上是噪声。
小波阀值收缩法增强去噪的具体处理过程为:将含噪图像在各尺度上进行小波分解, 保留大尺度低分辨率下的全部小波系数; 对于各尺度高分辨率下的小波系数, 可以设定一个阀值, 幅值低于该阀值的小波系数置为 0, 高于该阀值的小波系数或者完全保留, 或者做相应的“收缩” 处理。
结合上面的分析, 设定增强系数, 对每层得到的小波系数进行增强变换, 对经过上述处理的小波系数进行逆小波变换, 即对图像进行重构,可得到处理结果即增强后的图像。
2.2 阀值的选择采用阀值收缩法进行图像增强去噪, 最重要的一步是在小波域上对小波系数进行阈值操作, 阈值选择恰当与否直接影响到算法的有效性, 文献[4-6]给出了下列阀值公式N log 2σλ=, N 为图像的像素点数, σ为噪声的标准差。
如果已知图像的噪声,则从统计意义上可以知道各层小波变换子带的噪声的上界,并可以该上界为阈值来进行滤波。
在实际应用中,噪声的标准差一般是未知的,通常需要进行估计, 由于噪声经过小波变换后, 其能量大部分都集中在 HH 子图上, 因此也可以选取第一层小波分解系数的 HH 子图来估计σ,估计出来的方差要经过调整,例如将3σ[7]作为噪声方差代入上述公式,用HH 子图估计噪声的方差的通常做法就是用每层 HH 子图的方差作为该层的噪声方差,如果图像有较平缓的区域, 用 HH 子图中对应的该区域来做估计, 将能得到噪声的一个很好的估计。
由于噪声主要集中在最高分辨级1−J ,所以我们也可以利用小波系数{}1, 12, 2, 1, −−=J k J k W " 估计噪声标准差,如取 6745. ~, 12, 11k J k W median J −=−=" σ[2]。
2.3 基于小波变换与阀值收缩的图像增强算法算法描述:(1计算含噪声图像的正交小波变换。
对于像素点数为 N 的含噪图像 X ,不妨设JN 2=,利用10 四川理工学院学报(自然科学版 2006年 4月正交小波变换的 Mallat 快速算法获得分辨率(J L L <≤0下的尺度系数 {}L k L k V 2, 2, 1, , " =,及各分辨率下的小波系数{}j k j k J L L j W 2, , 1, 1, , 1, , , " " =−+=,其中尺度系数和小波系数共 N 个。
处理边界时, 常采用周期延拓方法。
(2对小波系数进行非线性阀值处理。
为保持图像的整体形状不变,保留所有的低频系数L k L k V 2, 2, 1, , " =。
取阀值N log σλ=[4], σ为图像噪声的标准差(度量噪声的强弱 ,对每个小波系数,采用软阀值和硬阀值方法进行处理:软阀值: ⎪⎩⎪⎨⎧+−=λλk j k j k j W W W , , , 0~ λλλ−≤<≥k j k j k j W W , , , 即将含噪声的图像的小波系数与所选定的阀值λ进行比较,大于阀值的点收缩为该点值与阀值的差值;小于阀值相反数的点收缩为该点值与阀值的和;幅值小于等于阀值的点变为零。
硬阀值: ⎩⎨⎧=0~, , k j k j W W λλ<≥k j kj , , 即把含噪图像的小波系数的绝对值与所选定的阀值λ进行比较,小于等于阀值的点变为零,大于阀值的点保持不变。
(3结合上面对小波系数的处理,设增强系数 i v W , ,其中, v 为尺度系数, i 取值1 ,2 ,3 , 分别代表 HH 、 HL 、 LH 子带图像,对每层得到的小波系数进行增强变换:k j i v k j W W W , , , ~~⋅=(4进行逆小波变换。
由所有低频尺度系数,以及经由阀值处理后的小波系数做逆小波变换进行重构,得到恢复的原始信号的估计值。
3 实验结果分析利用 matlab6.5在 PC 机上仿真试验,我们将原始图像人为地加上高斯白噪声,利用本文介绍的算法进行降噪增强处理,使用的是具有 8阶消失矩的近似对称小波sym8及正交周期小波变换,图像大小是 256*256, 高斯白噪声方差为 30, 阈值操作采用软阈值操作方法, 阈值采用文献[5]提出的公式和文献[2]提出的噪声估计方法。
并对没有使用阀值操作的算法也进行了实验,我们还对这两种方法进行了对 比,结果如图 1所示:图 1中(a 是原始图像, (b 是人为的加上高斯白噪声以后的图像, (c 是没有使用阀值操作的 算法产生的结果, (d 是采用软阀值操作的算法产生的结果。
从上图可以看出,没有使用阀值的算法的 细节信息损失很大,边缘特征不明显,不够清晰,而用软阀值法不但噪声几乎完全得到控制,而且较好 地保留了图像信号的边缘特征, 增强了图像细节信息, 视觉效果好, 处理结果明显优于没有使用阀值的 算法。