CN103208104A - 一种基于非局部理论的图像去噪方法 - Google Patents
一种基于非局部理论的图像去噪方法 Download PDFInfo
- Publication number
- CN103208104A CN103208104A CN2013101322436A CN201310132243A CN103208104A CN 103208104 A CN103208104 A CN 103208104A CN 2013101322436 A CN2013101322436 A CN 2013101322436A CN 201310132243 A CN201310132243 A CN 201310132243A CN 103208104 A CN103208104 A CN 103208104A
- Authority
- CN
- China
- Prior art keywords
- image
- details
- pixel
- edge
- positional information
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
基于非局部理论的图像去噪方法,包括如下步骤:输入含高斯白噪声图像u;设置算法中的参数;初步滤波得到u1;计算差图像uc1=u-u1;调用梯度算子(Candy算子、Roberts算子及Prewitt算子均可),检测出uc1的边缘ub1,并把它的位置当作所要提取的细节的位置信息,即检测出的边缘位置,也就是细节丢失的位置;利用步骤5中得到的细节的位置信息将u中丢失的细节ux1提取出来,包括位置信息所在像素点及其邻域的像素值,其他位置的像素值置零;将已经提取出的细节按比例叠加回图像u1;依据步骤2中的迭代参数ξ,ξ为一常数,然后利用步骤5中所调用的梯度算子所检测出的边缘图像矩阵ub1,计算矩阵ub1内所有像素的和θ,如θ≤ξ,则迭代停止,如θ≥ξ,则重复4~7步骤;迭代完成后,得到输出图像un。本发明的图像去噪新方法可以在滤除噪声的同时有效地保持图像的边缘细节信息。
Description
(一)技术领域
本发明属于图像处理技术中的图像去噪领域,具体涉及一种基于非局部理论的图像去噪方法。
(二)背景技术
数字图像在采集、转换以及传输时都会受到噪声的污染,因此,图像去噪是图像处理中最基本、最重要的问题,图像噪声的去除也一直是图像处理领域一个重要的课题。而广泛存在于实际数字图像中的噪声可以近似为高斯白噪声,所以去除高斯白噪声也成为去噪领域的重要部分。
自2005年,Buades等人提出了非局部均值算法(Nonlocal Means Filter,常缩写为NLM,见参考文献:Areview of image denoising algorithms with a new one,[J]multiscale model.Simul,2005,4:490-530),在空间域得到了较好的去噪效果,是目前最好的去噪算法之一,尤其是对高斯白噪声,图像去噪领域进入了一个新的阶段,其中的非局部理论已应用到图像处理的其他领域.然而非局部均值方法在去噪的同时会有部分纹理及边缘细节的丢失。非局部均值算法的应用有一定的局限性。本发明以非局部均值算法为基础,提出一种新算法,和NLM方法相比,在信噪比和边缘细节保持方面有明显提升。
(三)发明内容
本发明旨在克服非局部均值算法的不足,提供一种基于非局部理论的图像去噪方法,提高去噪效果。
本发明所述的基于非局部理论的图像去噪方法,包括如下步骤:
1、输入含高斯白噪声图像u;
2、设置算法中的参数:包括第一部分中的基本非局部均值算法(NLM)中的相关参数(搜索窗口大小、邻域窗口大小、滤波参数h)以及第二部分的迭代参数ξ的选取;
3、初步滤波:使用步骤2中选择的参数对输入图像u进行非局部均值滤波,得到的输出图像为丢失了部分细节的图像u1,作为下一步的基础;
4、计算差图像:uc1=u-u1;
5、调用梯度算子(Candy算子、Roberts算子及Prewitt算子均可),检测出uc1的边缘ub1,并把它的位置当作所要提取的边缘细节的位置信息,即检测出的边缘位置,也就是细节丢失的位置;
6、利用步骤5中得到的细节的位置信息将u中丢失的细节ux1提取出来,包括细节位置信息所在像素点及其邻域的像素值,其他位置的像素值置零;
7、将已经提取出的细节按比例叠加回图像u1:将细节ux1的每个像素点赋一个权值w(i,j),0≤w(i,j)≤1,(i,j为该像素点的横纵坐标),像素值为零的像素点权值为零,即w(i,j)=0;权值赋值原则为中心像素权值大于边缘像素权值,得到图像u2=w(i,j)ux1+(1-w(i,j))u1;
8、依据步骤2中的迭代参数ξ,ξ为一常数,然后利用步骤5中所调用的梯度算子所检测出的边缘图像矩阵ub1,计算矩阵ub1内所有像素的和θ,如θ≤ξ,则迭代停止,如θ≥ξ,则重复4~7步骤;
9、迭代完成后,得到输出图像un。
本发明的技术方案:本发明的算法结构分为两大部分,第一步是初步滤波部分,输入图像u,并选择合适的参数,应用非局部均值算法得到一个初步滤波后的图像,即丢失部分细节的图像u1。第二步细节添加部分,是通过模式识别技术将图像u1丢掉的细节逐步叠加融合回图像u1:计算差图像uc1=u-u1,识别出uc1中的边缘信息,用这些信息提取出图像u的细节ux1,然后将细节ux1叠加融合回图像u1,得到图像u2=ux1+u1,设置迭代参数ξ,当所提取的差图像边缘少于ξ时最后得到输出图像un。
本发明的优越性:将传统去噪方法的顺序逆转,用先得到初步滤波后的图像,然后添加回细节的思想来代替逐步迭代得到最优图像的方法。将图像丢掉的细节位置信息从差图像中识别出来并根据此信息得到图像的细节并融合回到初步滤波后的图像中去。可以在滤除噪声的同时有效地保持图像的边缘细节信息。
(四)附图说明
附图1为本发明的算法结构示意图;
附图2为本发明的迭代过程中细节提取的范围示意图;
附图3为本发明算法对lena图像进行仿真的结果。
附图4为本发明算法对onion图像进行仿真的结果。
附图5为本发明算法对cameraman图像进行仿真的结果。
(五)具体实施方式
本发明的算法实现步骤分为两大部分,第一部分是初步滤波部分;第二部分是迭代更新图像部分。具体步骤如下:
1、输入含高斯白噪声图像u;
2、设置算法中的一些参数:包括第一部分中的基本非局部均值算法中的相关参数(搜索窗口大小、邻域窗口大小、滤波参数h)以及第二部分的迭代参数ξ选取;
3、初步滤波:使用步骤2中选择的参数对输入图像u进行非局部均值滤波,得到的输出图像为丢失了部分细节的图像u1,作为下一步的基础;
4、计算差图像:uc1=u-u1;
5、调用边缘检测算子(Candy算子、Roberts算子及Prewitt算子均可),检测出uc1的边缘ub1,并把它的位置当作所要提取的边缘细节的位置信息(即检测出的边缘位置,也就是细节丢失的位置);
6、利用步骤5中得到的细节的位置信息将u中丢失的细节ux1提取出来,包括细节位置信息所在像素点及其邻域的像素值,其他位置的像素值置零。详细见附图2所示;
7、将已经提取出的细节按比例叠加回图像u1:将细节ux1的每个像素点赋一个权值w(i,j),0≤w(i,j)≤1,(i,j为该像素点的横纵坐标),像素值为零的像素点权值为零,即w(i,j)=0。权值赋值原则为中心像素权值大于边缘像素权值,得到图像u2=w(i,j)ux1+(1-w(i,j))u1;
8、依据步骤2中的迭代参数ξ,ξ为一常数,然后利用步骤5中所调用的梯度算子所检测出的边缘图像矩阵ub1,计算矩阵ub1内所有像素的和θ,如θ≤ξ,则迭代停止,如θ≥ξ,则重复4~7步骤;
9、迭代完成后,得到输出图像un。
本发明的效果通过以下实验证实:
一、实验条件:
使用matlab仿真软件对lena图像进行测试,输入图像所含高斯白噪声的标准差为σ=5,选取搜索窗口15×15,邻域窗口5×5,滤波参数h=1.35。未设置ξ,取4次迭代结果。
二、实验内容:
按照上面所述的实验步骤进行实验仿真,并将本发明算法与非局部均值算法进行比较。
三、实验结果:
实验结果见表1和附图3,附图4和附图5。实验结果表明本发明算法随着迭代次数的增多,信噪比其中信号能量 噪声能量 也逐渐提高,并有更好的边缘等细节保持性能。本发明算法比非局部均值算法有明显更好的去噪效果和更有效地保持图像的边缘细节信息。
表一去噪前后图像信噪比
以上所述,仅是本发明的较佳实施例而已,并不对本发明做形式上的限制,凡是依据本发明对以上实例所做的简单修改、等同变化与修饰,均仍属本发明技术方案的范围内。
Claims (1)
1.基于非局部理论的图像去噪方法,包括如下步骤:
步骤1、输入含高斯白噪声图像u;
步骤2、设置算法中的参数:包括第一部分中的基本非局部均值算法(NLM)中的相关参数(搜索窗口大小、邻域窗口大小、滤波参数h)以及第二部分的迭代参数ξ的选取;
步骤3、初步滤波:使用步骤2中选择的参数对输入图像u进行非局部均值滤波,得到的输出图像为丢失了部分细节的图像u1,作为下一步的基础;
步骤4、计算差图像:uc1=u-u1;
步骤5、调用梯度算子(Candy算子、Roberts算子及Prewitt算子均可),检测出uc1的边缘ub1,并把它的位置当作所要提取的边缘细节的位置信息,即检测出的边缘位置,也就是细节丢失的位置;
步骤6、利用步骤5中得到的细节的位置信息将u中丢失的细节ux1提取出来,包括细节位置信息所在像素点及其邻域的像素值,其他位置的像素值置零;
步骤7、将已经提取出的细节按比例叠加回图像u1:将细节ux1的每个像素点赋一个权值w(i,j),0≤w(i,j)≤1,(i,j为该像素点的横纵坐标),像素值为零的像素点权值为零,即w(i,j)=0;权值赋值原则为中心像素权值大于边缘像素权值,得到图像u2=w(i,j)ux1+(1-w(i,j))u1;
步骤8、依据步骤2中的迭代参数ξ,ξ为一常数,然后利用步骤5中所调用的梯度算子所检测出的边缘图像矩阵ub1,计算矩阵ub1内所有像素的和θ,如θ≤ξ,则迭代停止,如θ≥ξ,则重复4~7步骤;
步骤9、迭代完成后,得到输出图像un。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310132243.6A CN103208104B (zh) | 2013-04-16 | 2013-04-16 | 一种基于非局部理论的图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310132243.6A CN103208104B (zh) | 2013-04-16 | 2013-04-16 | 一种基于非局部理论的图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103208104A true CN103208104A (zh) | 2013-07-17 |
CN103208104B CN103208104B (zh) | 2015-10-07 |
Family
ID=48755321
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310132243.6A Active CN103208104B (zh) | 2013-04-16 | 2013-04-16 | 一种基于非局部理论的图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103208104B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103700076A (zh) * | 2013-12-26 | 2014-04-02 | 辽宁师范大学 | 基于非局部均值法的视频图像快速去噪方法 |
CN104036465A (zh) * | 2014-06-17 | 2014-09-10 | 南京邮电大学 | 基于边缘检测的自适应非局部均值去噪方法 |
CN106504207A (zh) * | 2016-10-24 | 2017-03-15 | 西南科技大学 | 一种图像数据处理方法 |
CN107248148A (zh) * | 2017-06-14 | 2017-10-13 | 上海晔芯电子科技有限公司 | 图像降噪方法及系统 |
CN110097518A (zh) * | 2019-04-28 | 2019-08-06 | 东软医疗系统股份有限公司 | 图像去噪方法、装置及终端设备 |
CN111402172A (zh) * | 2020-03-24 | 2020-07-10 | 湖南国科微电子股份有限公司 | 一种图像降噪方法、系统、设备及计算机可读存储介质 |
CN112184583A (zh) * | 2020-09-28 | 2021-01-05 | 成都微光集电科技有限公司 | 一种图像降噪方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012121280A1 (ja) * | 2011-03-09 | 2012-09-13 | 日本電信電話株式会社 | 映像符号化/復号方法,映像符号化/復号装置およびそのプログラム |
WO2012121269A1 (ja) * | 2011-03-09 | 2012-09-13 | 日本電信電話株式会社 | 画像処理方法,画像処理装置,映像符号化/復号方法,映像符号化/復号装置およびそれらのプログラム |
CN102800078A (zh) * | 2012-07-20 | 2012-11-28 | 西安电子科技大学 | 非局部均值图像细节修复方法 |
CN102842123A (zh) * | 2012-07-12 | 2012-12-26 | 南京理工大学 | 图像上采样中边缘清晰度提升的稀疏域残差补偿修正方法 |
CN103020918A (zh) * | 2013-01-09 | 2013-04-03 | 西安电子科技大学 | 基于形状自适应邻域均值的非局部均值去噪方法 |
-
2013
- 2013-04-16 CN CN201310132243.6A patent/CN103208104B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012121280A1 (ja) * | 2011-03-09 | 2012-09-13 | 日本電信電話株式会社 | 映像符号化/復号方法,映像符号化/復号装置およびそのプログラム |
WO2012121269A1 (ja) * | 2011-03-09 | 2012-09-13 | 日本電信電話株式会社 | 画像処理方法,画像処理装置,映像符号化/復号方法,映像符号化/復号装置およびそれらのプログラム |
CN102842123A (zh) * | 2012-07-12 | 2012-12-26 | 南京理工大学 | 图像上采样中边缘清晰度提升的稀疏域残差补偿修正方法 |
CN102800078A (zh) * | 2012-07-20 | 2012-11-28 | 西安电子科技大学 | 非局部均值图像细节修复方法 |
CN103020918A (zh) * | 2013-01-09 | 2013-04-03 | 西安电子科技大学 | 基于形状自适应邻域均值的非局部均值去噪方法 |
Non-Patent Citations (6)
Title |
---|
BO XIONG 等: "A Universal Denoising Framework With a New Impulse Detector and Nonlocal Means", 《IMAGE PROCESSING, IEEE TRANSACTIONS ON (VOLUME:21 , ISSUE: 4 )》, vol. 21, no. 4, 19 October 2011 (2011-10-19), pages 1663 - 1675, XP011491971, DOI: 10.1109/TIP.2011.2172804 * |
DENG SHI-WO 等: "Image denoising study on Nonlocal-Means Filter", 《IMAGE ANALYSIS AND SIGNAL PROCESSING (IASP), 2012 INTERNATIONAL CONFERENCE ON》, 9 November 2012 (2012-11-09), pages 1 - 4, XP032322126, DOI: 10.1109/IASP.2012.6424997 * |
HAO XU 等: "On the biased estimation of nonlocal means filter", 《MULTIMEDIA AND EXPO, 2008 IEEE INTERNATIONAL CONFERENCE ON》, 26 April 2008 (2008-04-26), pages 1149 - 1152 * |
单建华: "改进权值函数的非局部均值去噪算法", 《中国图象图形学报》, vol. 17, no. 10, 31 October 2012 (2012-10-31), pages 1227 - 1231 * |
易子麟 等: "基于非局部均值滤波的SAR图像去噪", 《电子与信息学报》, vol. 34, no. 4, 30 April 2012 (2012-04-30), pages 950 - 955 * |
杨学志 等: "基于非局部均值滤波的结构保持相干斑噪声抑制方法", 《中国图象图形学报》, vol. 14, no. 12, 31 December 2009 (2009-12-31), pages 2443 - 2450 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103700076B (zh) * | 2013-12-26 | 2016-09-14 | 辽宁师范大学 | 基于非局部均值法的视频图像快速去噪方法 |
CN103700076A (zh) * | 2013-12-26 | 2014-04-02 | 辽宁师范大学 | 基于非局部均值法的视频图像快速去噪方法 |
CN104036465A (zh) * | 2014-06-17 | 2014-09-10 | 南京邮电大学 | 基于边缘检测的自适应非局部均值去噪方法 |
CN106504207A (zh) * | 2016-10-24 | 2017-03-15 | 西南科技大学 | 一种图像数据处理方法 |
CN106504207B (zh) * | 2016-10-24 | 2019-08-02 | 西南科技大学 | 一种图像数据处理方法 |
CN107248148B (zh) * | 2017-06-14 | 2023-08-25 | 思特威(上海)电子科技股份有限公司 | 图像降噪方法及系统 |
CN107248148A (zh) * | 2017-06-14 | 2017-10-13 | 上海晔芯电子科技有限公司 | 图像降噪方法及系统 |
CN110097518A (zh) * | 2019-04-28 | 2019-08-06 | 东软医疗系统股份有限公司 | 图像去噪方法、装置及终端设备 |
CN110097518B (zh) * | 2019-04-28 | 2022-12-27 | 东软医疗系统股份有限公司 | 图像去噪方法、装置及终端设备 |
CN111402172B (zh) * | 2020-03-24 | 2023-08-22 | 湖南国科微电子股份有限公司 | 一种图像降噪方法、系统、设备及计算机可读存储介质 |
CN111402172A (zh) * | 2020-03-24 | 2020-07-10 | 湖南国科微电子股份有限公司 | 一种图像降噪方法、系统、设备及计算机可读存储介质 |
CN112184583A (zh) * | 2020-09-28 | 2021-01-05 | 成都微光集电科技有限公司 | 一种图像降噪方法及装置 |
CN112184583B (zh) * | 2020-09-28 | 2023-11-17 | 成都微光集电科技有限公司 | 一种图像降噪方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103208104B (zh) | 2015-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103208104B (zh) | 一种基于非局部理论的图像去噪方法 | |
CN107016642B (zh) | 用于对有噪输入图像进行分辨率上调的方法和装置 | |
Fan et al. | Median filtered image quality enhancement and anti-forensics via variational deconvolution | |
CN103247036A (zh) | 多曝光图像融合方法和装置 | |
CN103049906A (zh) | 一种图像深度提取方法 | |
CN103839234B (zh) | 一种基于可控核的双几何非局部均值图像去噪方法 | |
Kaur et al. | Comparative study of different edge detection techniques | |
CN103208097A (zh) | 图像多方向形态结构分组的主分量分析协同滤波方法 | |
CN102930508B (zh) | 基于图像残余信号的非局部均值图像去噪方法 | |
Shenbagavalli et al. | Classification of soil textures based on laws features extracted from preprocessing images on sequential and random windows | |
CN103337053A (zh) | 一种基于开关非局部全变分的椒盐噪声污染图像滤波方法 | |
Vijayakumar | Selective image enhancement and restoration for skin cancer identification | |
CN102663706A (zh) | 基于菱形模板的自适应加权均值滤波方法 | |
CN101655977B (zh) | 基于差分图像检测、多窗口滤波去除图像脉冲噪声的方法 | |
CN104778710A (zh) | 一种基于量子理论的形态学图像边缘检测方法 | |
CN105513014A (zh) | 一种多帧图像超分辨率重建方法及其重建系统 | |
CN110782385A (zh) | 一种基于深度学习的图像水印的去除方法 | |
CN105869133A (zh) | 一种基于非因果分数阶次微分的图像锐化方法 | |
CN104463811A (zh) | 基于能量泛函的图像平滑与锐化算法 | |
CN102118547A (zh) | 图像加权滤波的方法 | |
CN101859385A (zh) | 基于图像的局部模糊篡改盲检测方法 | |
CN101957912B (zh) | 基于多尺度各向异性扩散的人脸光照不变图像的获取方法 | |
CN103745442B (zh) | 基于非局部小波系数收缩的图像去噪方法 | |
CN105930811A (zh) | 基于图像处理的手掌纹理特征检测方法 | |
CN102903087B (zh) | 基于gpu编程的sar图像去噪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |