一种改进小波阈值函数的图像去噪方法研究
- 格式:pdf
- 大小:2.59 MB
- 文档页数:6
第17卷 第4期2019年12月南京工程学院学报(自然科学版)JournalofNanjingInstituteofTechnology(NaturalScienceEdition)Vol.17ꎬNo.4Dec.ꎬ2019
doi:10.13960/j.issn.1672-2558.2019.04.008投稿网址:http://xb.njit.edu.cn
一种改进小波阈值函数的图像去噪方法研究
周 峡1ꎬ徐善顶2
(1.南京工程学院经济与管理学院ꎬ江苏南京211167ꎻ
2.南京工程学院数理部ꎬ江苏南京211167)
摘 要:性质良好的小波阈值去噪方法在图像去噪领域中应用十分广泛.考虑到传统硬、软阈值函数存在不连续和
恒定偏差的缺陷ꎬ基于均方根插值阈值函数引入一种改进的阈值函数图像去噪方法.该方法通过调节参数ꎬ能够更
好地适应不同图像不同噪声的去噪需求ꎬ且对于图像的边缘细节信息能够处理得更好.试验结果表明ꎬ改进阈值函
数图像去噪方法较之现有方法具有更小的均方误差、更高的峰值信噪比、更高的结构相似性和更小的梯度幅值相
似度偏差ꎬ去噪效果更优.
关键词:图像去噪ꎻ小波阈值ꎻ均方根插值ꎻ改进阈值函数
中图分类号:TN911.73
收稿日期:2019-07-23ꎻ修回日期:2019-09-29基金项目:南京工程学院校级科研基金项目(QKJ201804ꎬCKJB201606)作者简介:周峡ꎬ硕士ꎬ工程师ꎬ研究方向为数字图像处理、信息管理与数据应用.E ̄mail:zxxzy11@gmail.com
引文格式:周峡ꎬ徐善顶.一种改进小波阈值函数的图像去噪方法研究[J].南京工程学院学报(自然科学版)ꎬ2019ꎬ17(4):44-49. 图像在形成、记录、传输和显示等过程中ꎬ由于
传感器性能不稳定、数据传输网络故障、人为或自
然环境带来的干扰和破坏等因素ꎬ常常导致采集到
的图像有噪声产生ꎬ去除图像中的噪声ꎬ一直是图
像处理领域中非常有必要的一项工作.近几十年
来ꎬ随着人们对小波的研究不断深入ꎬ小波作为一
种功能强大的分析工具ꎬ已在图像去噪领域中得到
了非常广泛的应用.其中ꎬ小波阈值去噪方法由于
实现简单并且可以达到较好的去噪效果ꎬ已成为图
像去噪应用的主流技术.小波阈值去噪法的关键在
于阈值函数的设计ꎬ其在一定程度上关系到图像去
噪的质量.在现有的研究中ꎬ已经有很多阈值函数
的设计方法ꎬ如传统的硬阈值函数和软阈值函
数[1-2].但硬阈值函数不连续ꎬ去噪过程中会产生
振荡ꎻ软阈值函数虽然连续ꎬ但经处理得到的图像
比较平滑ꎬ会造成图像边缘细节的模糊失真现
象[3].针对软、硬阈值函数的不足ꎬ不同的改进阈
值函数算法相继被提出[4-10]ꎬ文献[4]提出的均方根插值阈值函数就是其中极具代表性的一种方法ꎬ
均方根插值阈值函数既能降低硬阈值函数法产生
的振荡ꎬ又能减小软阈值函数法引入的恒定偏差ꎬ
但是其阈值范围内直接置零ꎬ会损失图像部分有用
信息.本文针对均方根插值阈值函数存在的不足ꎬ
提出一种改进的阈值函数图像去噪方法ꎬ使去噪形
式更灵活、图像信息损失更少ꎬ通过试验证明改进
方法较传统硬、软阈值函数和均方根插值阈值函数
的去噪效果更好ꎬ图像信息保留更完整.
1 小波阈值去噪法
1.1 去噪原理及步骤
利用小波变换进行图像去噪ꎬ其本质是要找到
图像本身与干扰噪声在小波域的可分离性[11].含
噪声的图像信息在经各尺度的小波变换后ꎬ图像的
小波系数绝对值要大于噪声的小波系数绝对值ꎬ通
过选择合适的阈值ꎬ使图像对应的小波系数予以保
留ꎬ而噪声对应的小波系数尽量缩小为零ꎬ
从而达第17卷第4期周 峡ꎬ等:一种改进小波阈值函数的图像去噪方法研究
到有效去除噪声的目的.
小波阈值图像去噪法的基本步骤[12]为:
1)小波分解ꎬ选择合适的小波基和分解层数ꎬ
对待处理的含噪图像进行小波分解ꎻ
2)阈值量化ꎬ选择合适的阈值函数和阈值ꎬ对
分解得到的小波系数的高频部分进行量化处理ꎻ
3)小波重构ꎬ通过重构处理后的小波系数得
到去噪后的图像.
1.2 传统的小波阈值函数
传统的小波阈值函数设计方式有硬阈值函数
和软阈值函数.
硬阈值函数的表达式为:
^ωjꎬk=ωjꎬkꎬ|ωjꎬk|≥T
0ꎬ|ωjꎬk|
式中:ωjꎬk、^ωjꎬk分别为阈值量化前、量化后的小波系
数ꎻT为阈值.当|ωjꎬk|≥T时ꎬ保持该小波系数不
变ꎻ当|ωjꎬk|
见ꎬ硬阈值函数在阈值处不连续ꎬ这就导致小波重
构时容易出现伪吉布斯现象[1].
软阈值函数的表达式为:
^ωjꎬk=sgn(ωjꎬk)(|ωjꎬk|-T)ꎬ|ωjꎬk|≥T
0ꎬ|ωjꎬk|
式中:ωjꎬk、^ωjꎬk分别为阈值量化前、量化后的小波系
数ꎻT为阈值.软阈值函数在阈值处连续ꎬ但当
|ωjꎬk|≥T时ꎬ量化处理后的小波系数与真实图像的
小波系数之间存在恒定偏差ꎬ这势必会给重构图像
带来不可避免的偏差[2].
1.3 均方根插值阈值函数
均方根插值阈值函数的表达式为:
^ωjꎬk=sgn(ωjꎬk)ω2jꎬk+(|ωjꎬk|-T)22ꎬ|ωjꎬk|≥T
0ꎬ|ωjꎬk|
îíïï
ïï(3)
ω2jꎬk+(|ωjꎬk|-T)22总介于|ωjꎬk|和|ωjꎬk|-T
之间ꎬ由于用软阈值函数量化处理后的^ωjꎬk绝对值总要比ωjꎬk小T(|ωjꎬk|≥T时)ꎬ因此要尽量减小这
种偏差ꎬ但若把这种偏差减小为0ꎬ恰好是硬阈值函
数的情况ꎬ也不是最好的办法ꎬ因为ωjꎬk本身就是由
图像对应的小波系数和噪声对应的小波系数组成
的ꎬ为了使|^ωjꎬk-ωjꎬk|更小ꎬ令|^ωjꎬk|的取值为|ωjꎬk|
和|ωjꎬk|-T的均方根插值ꎬ可使处理后的小波系数
^ωjꎬk更接近ωjꎬk[4]ꎬ然而此阈值函数也是非连续的.
1.4 改进的阈值函数综合前面所述各种阈值函数的优点与不足ꎬ本
文引入调节参数μ和以下函数对均方根插值阈值
函数进行改进.
f(x)=xe-1x2ꎬx≠0
0ꎬx=0{(4)
改进的阈值函数为:
^ωjꎬk=sgn(ωjꎬk)μω2jꎬk+(1-μ)(|ωjꎬk|-T)2ꎬ
|ωjꎬk|≥T
μωjꎬke1T2-1ω2jꎬkꎬ|ωjꎬk|
îíïï
ïï
(5)
式中:ωjꎬk和^ωjꎬk分别为阈值量化前和量化后的小
波系数ꎻT为阈值ꎬμ为调节参数ꎬ0≤μ≤1ꎬ将μ的
取值做适当的调节ꎬ即可取得更好的去噪效果.当
μ=0时ꎬ式(5)即为软阈值函数ꎻ当μ=1时ꎬ对于
|ωjꎬk|≥T部分ꎬ式(5)与硬阈值函数一致.由于较
小的小波系数也是由部分图像信息和噪声构成的ꎬ
如果直接将其置零ꎬ就会使有用的图像信息产生一
定的损失.本文对较小的小波系数进行处理ꎬ保留
其中对应的部分有用图像信息ꎬ能在去噪的同时更
好地保留图像的边缘细节信息ꎬ并且优化该函数在
|ωjꎬk|=T时的连续性ꎬ在图像重建时避免了一些
附加振荡.当ωjꎬk∈[-1ꎬ1]ꎬT=0.5ꎬμ=0.2
时ꎬ
四种阈值函数的函数示意图如图1所示.
2 仿真试验研究
本文应用Matlab软件编程实现所有程序ꎬ选
取Peppers、Lena和Lion三种不同类型的图像作为54南京工程学院学报(自然科学版)2019年12月
图1 小波阈值函数示意图
试验对象ꎬ并在这三幅图像中加入均值为0ꎬ方差
分别为0.005、0.01和0.02的高斯噪声ꎬ用sym4小波对其进行分解ꎬ分解层数为3层ꎬ采用统一阈
值T=σ2log(M×N)去噪ꎬ其中ꎬM×N为图像
的像素数ꎬσ为噪声的标准差.
2.1 试验参数的选择
从表达式(5)可以看出ꎬ改进算法的去噪效果
依赖于参数μ.考虑到图像的多样性ꎬ本文在对参
数μ的选择上仍需改进ꎬ这里通过试验来对参数进
行合理的选择ꎬ以峰值信噪比(peaksignalnoise
ratioꎬPSNR)为客观评价指标ꎬPNSR值越大ꎬ表明
去噪效果越好.图2显示了各图像在不同μ值下去
噪的PSNR值.
图2 不同μ值下的PSNR值
从图2可以看出ꎬ对于Peppers、Lena和Lion
图像ꎬ当μ值分别在0.6、0.5、0.7处ꎬ去噪图像的
PSNR值达到最大ꎬ表明此时本文算法的去噪效果
最好.
2.2 算法有效性验证
为了对比分析改进阈值法的优越性ꎬ对硬阈值
法、软阈值法、均方根插值阈值法和改进阈值法四
种小波阈值去噪算法进行试验ꎬ去噪结果如图3至
图5所示.
从图3到图5可以看出ꎬ经过改进阈值法处理
之后ꎬ加噪图像中的大部分噪声得到了去除ꎬ且相
对于硬阈值法、软阈值法及均方根插值阈值法而
言ꎬ改进阈值法得到的图像具有更高的清晰度ꎬ对
于图像的边缘细节信息也保留的更好ꎬ达到了更好的去噪效果.对于不同图像含有不同强度的噪声ꎬ
各阈值函数法的去噪能力也不同ꎬ而改进阈值法可
以通过调节μ来使其具有更好的灵活性ꎬ以适应不
同的情况.
为了客观评价各阈值函数的去噪效果ꎬ本文采
用的评价指标为:均方误差(meansquareerrorꎬ
MSE)、峰值信噪比、结构相似性测度(structual
similarityꎬSSIM)和梯度幅值相似度偏差(gradient
magnitudesimilaritydeviationꎬGMSD)[13-14].根据图
像去噪的要求ꎬ期望改进阈值法处理后的图像具有
较小的均方误差、较大的峰值信噪比、较大的结构
相似性测度及较小的梯度幅值相似度偏差ꎬ从而说
明去噪后的图像更接近原始图像ꎬ也即去噪效果更
好.结果见表
1至表3.64