CN104200434B - 一种基于噪声方差估计的非局部均值图像去噪方法 - Google Patents
一种基于噪声方差估计的非局部均值图像去噪方法 Download PDFInfo
- Publication number
- CN104200434B CN104200434B CN201410431699.7A CN201410431699A CN104200434B CN 104200434 B CN104200434 B CN 104200434B CN 201410431699 A CN201410431699 A CN 201410431699A CN 104200434 B CN104200434 B CN 104200434B
- Authority
- CN
- China
- Prior art keywords
- noise
- image
- noise image
- estimation
- extension
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000009499 grossing Methods 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 4
- 230000003362 replicative effect Effects 0.000 claims description 6
- 235000013399 edible fruits Nutrition 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 4
- 238000001914 filtration Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Landscapes
- Image Processing (AREA)
Abstract
本发明属于数字图像处理技术领域,具体涉及一种应用于图像去噪并作为后续目标识别的预处理的基于噪声方差估计的非局部均值图像去噪方法。本发明包括:输入噪声图像,获取噪声图像尺寸;生成一个与噪声图像相同尺寸的零矩阵;对噪声图像边缘进行对称扩展;估计噪声方差,确定全局平滑参数;遍历噪声图像中每个像素,计算权值;利用非局部均值算法计算去噪图像。本发明中基于噪声方差估计的非局部均值图像去噪方法,可以显著提高噪声图像清晰度,并且去噪后能更清晰地保留边缘和细节信息。
Description
技术领域
本发明属于数字图像处理技术领域,具体涉及一种应用于图像去噪并作为后续目标识别的预处理的基于噪声方差估计的非局部均值图像去噪方法。
背景技术
图像去噪作为图像处理中的最为基础与核心的技术,是保证后续图像处理顺利实现的前提,又被称作图像滤波。其最终目的是改善实际图像受到噪声干扰引起的质量下降,通过各种技术手段的应用,有效提高图像的视觉质量及信噪比,更好的还原出图像的本质信息,作为重要的预处理手段为后续操作做好准备。
图像中存在的噪声一般可被认为高斯白噪声,传统去噪方法大致可分为基于空域和变换域的两大类。这些方法在去除噪声的同时常常会引起细节与边缘的缺失。实际上图像中的每个像素并非孤立存在,与其周围像素一起构成图像中的几何结构,以目标像素为中心的领域称为图像块。而且处于不同位置的像素表现出很强的相关性,即图像的自相似性。
Buades等人基于自然图像中冗余信息的存在,利用图像中的自相似性提出了非局部均值算法。其核心思想是以目标像素为中心选择合适的图像块,然后在整幅图像中搜索具有相似结构的图像块,以各图像块间相似度对于目标像素进行调整。但该方法中存在噪声图像的参数选取问题,包括图像块大小的选择与噪声方差的估计。传统的经验估计方法不能很好地获得准确的噪声方差,而过于复杂地估计方法又会导致整体计算量的加大。
发明内容
本发明的目的在于针对现有非局部均值图像去噪方法中存在的不完善之处,提出了一种通过噪声估计以及非局部均值滤波方法的结合,有效保持图像边缘和目标特性的同时实现图像快速去噪,为后续处理和分析提供准确的目标特性和边缘保持度的基于噪声方差估计的非局部均值图像去噪方法。
本发明的目的是这样实现的:
(1)输入噪声图像,获取噪声图像尺寸W×H;
(2)生成一个与噪声图像相同尺寸的零矩阵;
(3)对噪声图像边缘进行对称扩展,经扩展后噪声图像尺寸由W×H变为(W+4)×(H+4);
(4)估计噪声方差,确定全局平滑参数h;
(5)遍历噪声图像中每个像素,计算权值w(x,y);
(6)利用非局部均值算法计算去噪图像NL[I](x)。
对噪声图像边缘进行对称扩展具体为:
(3.1)复制噪声图像第一行与第H行的灰度值,并分别添加至噪声图像图像的顶部与底部,则噪声图像尺寸变为W×(H+2);
(3.2)复制原噪声图像的第二行与第H-1行的灰度值,并分别添加至扩展后噪声图像图像的顶部与底部,则噪声图像尺寸变为W×(H+4);
(3.3)复制扩展后噪声图像第一列与第W列的灰度值,并分别添加至扩展后噪声图像图像的左侧与右侧,则噪声图像尺寸变为(W+2)×(H+4);
(3.4)复制扩展后噪声图像第二列与第W-1列的灰度值,并分别添加至扩展后噪声图像图像的左侧与右侧,最终噪声图像尺寸变为(W+4)×(H+4)。
确定全局平滑参数h具体为:
(4.1)生成噪声估计算子N,
两个方向元素分别为:
噪声估计算子N为:
(4.2)将噪声图像与噪声估计算子进行二维卷积操作;
(4.3)计算获得噪声方差,并将噪声标准差确定为全局平滑参数h:
噪声估计算子具有零均值,每个像素位置的噪声标准差为则噪声估计算子方差为
I(x)*N表示噪声估计算子与图像进行卷积运算,W×H的图像与3×3的算子卷积后的结果尺寸为(W+2)×(H+2),其中真正属于原图像的范围为(W-2)×(H-2),噪声方差估计值为:
本发明的有益效果在于:本发明中基于噪声方差估计的非局部均值图像去噪方法,可以显著提高噪声图像清晰度,并且去噪后能更清晰地保留边缘和细节信息。
附图说明
图1为本发明实施例的图像去噪方法示范性流程图;
图2为本发明实施例使用的测试图像;
图3为本发明实施例图2加入标准差20的高斯噪声后的噪声图像;
图4为利用小波硬阈值去噪方法,维纳滤波方法以及本发明方法对图3进行去噪的结果图。
具体实施方式
下面结合附图对本发明做进一步描述。
噪声检测应当对图像中边缘具有敏感度,因此在估计算子的选取上应考虑到不同的方向性,假设两个方向元素分别为:
噪声估计算子N通过下式获得
噪声估计算子具有零均值,假设每个像素位置的噪声标准差为则噪声估计算子方差为
I(x)*N表示噪声估计算子与图像进行卷积运算。W×H的图像与3×3的算子卷积后的结果尺寸为(W+2)×(H+2)。考虑到图像边缘的扩展,其中真正属于原图像的范围为(W-2)×(H-2),因此噪声方差估计值为
实现本发明的技术方案是以非局部均值滤波为工具,采用噪声估计算子对原始噪声图像进行噪声方差估计,然后计算估计所得的噪声标准差作为全局平滑参数h,利用非局部均值完成图像去噪,最终获得去除噪声的图像。
本发明的目的是这样实现的:
(1.1)读取噪声图像
(1.2)通过噪声估计算子计算图像噪声的方差,获得全局平滑参数h
(1.3)选取合适的搜索窗口尺寸、相似性窗口尺寸、利用全局平滑参数h完成图像非局部均值去噪
本发明的核心技术内容在于图像噪声的方差估计以及非局部均值滤波算法的综合应用。应用本发明所去噪的图像能为用户带来功能更好的视觉体验且利于进一步的处理操作。
假设f(x)表示原始无噪声图像,图像的尺寸为W×H像素。在独立同分布的加性高斯零均值噪声的作用下,噪声图像I(x)可以表示为
I(x)=f(x)+n(x),x∈X (3)
其中x=(x1,x2)表示像素的空间坐标位置,X表示离散图像域。
下面结合附图1和具体实施例对本发明作进一步详细说明。
对于给定的噪声图像I(x),图像中每一个像素的输出估计值由整幅图像中所有像素的加权平均决定
有着相似灰度级的像素在平均计算中被分配更大的权值w(x,y)。权值函数为指数形式,两像素点间欧氏距离越小,则权值越大。在搜索区域中的每个像素点y,二范数表示灰度值向量i(x)和i(y)间的欧氏距离。
其中Z(x)为归一化常数。h为全局平滑参数,它决定指数函数的衰减速度并影响滤波的程度。权值w(x,y)的确定依赖于像素x与像素y间的相似度,且满足0≤w(x,y)≤1以及∑yw(x,y)=1。通过遍历计算图像中所有像素得到去噪结果。
步骤1,输入含有噪声的图像后,获取噪声图像尺寸W×H;
步骤2,生成一个与噪声图像相同尺寸的零矩阵;
步骤3,对于噪声图像边缘进行对称扩展,经扩展后噪声图像尺寸由W×H变为(W+4)×(H+4);
对于噪声图像边缘对称扩展,具体按如下步骤进行:
1.复制噪声图像第一行与第H行的灰度值,并分别添加至噪声图像图像的顶部与底部,则噪声图像尺寸变为W×(H+2);
2.复制原噪声图像的第二行与第H-1行的灰度值,并分别添加至扩展后噪声图像图像的顶部与底部,则噪声图像尺寸变为W×(H+4);
3.复制扩展后噪声图像第一列与第W列的灰度值,并分别添加至扩展后噪声图像图像的左侧与右侧,则噪声图像尺寸变为(W+2)×(H+4);
4.复制扩展后噪声图像第二列与第W-1列的灰度值,并分别添加至扩展后噪声图像图像的左侧与右侧,最终噪声图像尺寸变为(W+4)×(H+4)。
以一幅原始尺寸为2×2图像为例进行上述扩展:
步骤4,确定全局平滑参数h。
在本发明方法中选取噪声标准差的估计值,即
对于噪声方差的估计,具体按如下步骤进行:
1.依据式(1)生成噪声估计算子N;
2.将噪声图像与噪声估计算子进行二维卷积操作;
3.利用式(2)计算获得方差,同时可确定全局平滑参数h。
步骤5,遍历图像中每个像素,计算权值w(x,y);
步骤6,计算每个像素的输出估计值NL[I](x),输出去噪后图像。
本发明的效果可通过以下仿真实验进行证实:
实验条件:实验中使用的原始图像如图4所示,其中图2(a)为测试图像lena,图2(b)为测试图像cameraman,尺寸均为256×256像素,色彩深度为8位,即灰度值范围为0~255。测试图中包含丰富的边缘轮廓及细节。图3为加入噪声标准差为20的噪声图像。
实验内容:
(1)在上述实验条件下,分别使用小波硬阈值滤波,维纳滤波方法,以及本发明方法对图3进行去噪。结果如图4所示
(2)表1为利用本发明噪声方差估计方法所得的方差估计结果。σn和分别表示实际加入噪声标准差和估计所得标准差,δ表示估计的相对误差。11组数据的平均相对误差2.2%
(3)本发明采用峰值信噪比PSNR(Peak Signal-to-Noise Ratio)和结构相似度SSIM(Structural Similarity Index Measurement)作为去噪效果的评价指标。
PSNR定义为
其中f(i,j),分别为原始图像和去噪后图像对应位置的灰度值。实验结果中PSNR的单位为dB,数值越大表明图像去噪效果越好,图像的失真越少SSIM定义为
其中μf(x)和表示原始图像和去噪后图像的均值,σf(x)和表示原始图像和去噪后图像的标准差,表示原始图像和去噪后图像的协方差。C1,C2为给定常数。SSIM从图像组成的角度衡量两幅图像的相似度,其取值范围为0~1,实验结果中SSIM越接近1则表明去噪后图像越接近原始图像,图像的失真越少
对图2中的测试图像分别加入噪声标准差为10,20,30的高斯加性白噪声,对上述三种方法进行比较,各方法的去噪结果PSNR值与SSIM值分别如表2和表3所示。
表1本发明噪声估计结果
表2不同去噪方法实验结果PSNR对比
表3不同去噪方法实验结果SSIM对比
从表2中可以看出,本发明方法可以将PSNR提升1~3dB。表3的结果表明相比于其他方法,本发明方法的结果更接近于原始图像
图1为根据本发明实施例的图像去噪方法的示范性流程图。
如图1所示主要步骤为:
1.读取噪声图像
2.通过噪声估计算子计算图像噪声的方差
3.选取合适的搜索窗口尺寸,利用估计所得方差确定全局平滑参数并计算权值
4.输出去噪后图像
图2为本发明实施例的使用的测试图像。图2(a)为256×256大小的8bit lena图像,图2(b)为256×256大小的8bit cameraman图像。
图3为本发明实施例的使用的噪声图像。图3(a)为加入标准差20的高斯噪声后的lena图像,图3(b)加入标准差20的高斯噪声后的cameraman图像。
图4(a)~(b)为小波硬阈值去噪的结果示意图,图4(c)~(d)为维纳滤波的结果示意图,图4(e)~(f)为本发明方法结果示意图,噪声环境为标准差20的高斯噪声,通过结果对比可以看出,本发明方法更好地去除了噪声图像中的噪声点,提升了对比度,边缘保持效果更好。
综上所述,本发明中基于噪声方差估计的非局部均值图像去噪方法,可以显著提高噪声图像清晰度,并且去噪后能比较清晰地保留边缘和细节信息。本发明方法对于后续的图像分析处理有很大的帮助,根据不同的应用领域,本发明经过适当的修改同样可以适用于其他相关的图像处理领域。
最后应说明的是,以上所述实施例仅用以描述本发明的技术方案,而不是对本技术方法进行限制,本发明在应用上可以延伸为其他的修改、变化、应用和实施例,并且因此认为所有这样的修改、变化、应用、实施例都在本发明的精神和教导范围内。
Claims (1)
1.一种基于噪声方差估计的非局部均值图像去噪方法,包括如下步骤:
(1)输入噪声图像,获取噪声图像尺寸W×H;
(2)生成一个与噪声图像相同尺寸的零矩阵;
(3)对噪声图像边缘进行对称扩展,经扩展后噪声图像尺寸由W×H变为(W+4)×(H+4);
(4)估计噪声方差,确定全局平滑参数h;
(5)遍历噪声图像中每个像素,计算权值w(x,y);
(6)利用非局部均值算法计算去噪图像NL[I](x);
所述对噪声图像边缘进行对称扩展具体为:
(3.1)复制噪声图像第一行与第H行的灰度值,并分别添加至噪声图像图像的顶部与底部,则噪声图像尺寸变为W×(H+2);
(3.2)复制原噪声图像的第二行与第H-1行的灰度值,并分别添加至扩展后噪声图像图像的顶部与底部,则噪声图像尺寸变为W×(H+4);
(3.3)复制扩展后噪声图像第一列与第W列的灰度值,并分别添加至扩展后噪声图像图像的左侧与右侧,则噪声图像尺寸变为(W+2)×(H+4);
(3.4)复制扩展后噪声图像第二列与第W-1列的灰度值,并分别添加至扩展后噪声图像图像的左侧与右侧,最终噪声图像尺寸变为(W+4)×(H+4);
所述的确定全局平滑参数h具体为:
(4.1)生成噪声估计算子N,
两个方向元素分别为:
根据两个方向元素计算获得噪声估计算子N为:
(4.2)将噪声图像与噪声估计算子进行二维卷积操作;
(4.3)计算获得噪声方差,并将噪声标准差确定为全局平滑参数h:
噪声估计算子具有零均值,每个像素位置的噪声标准差为则噪声估计算子方差为
I(x)*N表示噪声估计算子与图像进行卷积运算,W×H的图像与3×3的算子卷积后的结果尺寸为(W+2)×(H+2),其中真正属于原图像的范围为(W-2)×(H-2),噪声方差估计值为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410431699.7A CN104200434B (zh) | 2014-08-28 | 2014-08-28 | 一种基于噪声方差估计的非局部均值图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410431699.7A CN104200434B (zh) | 2014-08-28 | 2014-08-28 | 一种基于噪声方差估计的非局部均值图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104200434A CN104200434A (zh) | 2014-12-10 |
CN104200434B true CN104200434B (zh) | 2017-04-12 |
Family
ID=52085720
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410431699.7A Expired - Fee Related CN104200434B (zh) | 2014-08-28 | 2014-08-28 | 一种基于噪声方差估计的非局部均值图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104200434B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104574305A (zh) * | 2014-12-25 | 2015-04-29 | 深圳市一体太赫兹科技有限公司 | 一种毫米波图像处理方法及系统 |
CN106023093A (zh) * | 2016-05-09 | 2016-10-12 | 吉林大学 | 一种基于改进图像块匹配的非局部均值图像去噪方法 |
CN108389173B (zh) * | 2018-03-24 | 2021-08-31 | 厦门大学嘉庚学院 | 一种基于机会成本的参数优化方法 |
CN113269687A (zh) * | 2021-05-17 | 2021-08-17 | 武汉博宇光电系统有限责任公司 | 一种基于参数优化的非局部平均红外图像去噪方法 |
CN116563136B (zh) * | 2023-03-27 | 2024-03-15 | 山东科技大学 | 一种侧扫声呐图像去噪方法、系统、存储介质和执行系统 |
CN117528274B (zh) * | 2023-11-13 | 2024-05-14 | 北京信息科技大学 | 基于深度学习的红外成像非均匀噪声提取及优化方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6718068B1 (en) * | 2000-03-10 | 2004-04-06 | Eastman Kodak Company | Noise reduction method utilizing statistical weighting, apparatus, and program for digital image processing |
CN101308573A (zh) * | 2008-06-30 | 2008-11-19 | 北京中星微电子有限公司 | 一种消除噪声的方法和装置 |
CN101655977A (zh) * | 2009-09-17 | 2010-02-24 | 浙江大学 | 基于差分图像检测、多窗口滤波去除图像脉冲噪声的方法 |
CN101661611A (zh) * | 2009-09-25 | 2010-03-03 | 西安电子科技大学 | 基于贝叶斯非局部均值滤波器的实现方法 |
CN103020919A (zh) * | 2013-01-09 | 2013-04-03 | 西安电子科技大学 | 基于非局部Lee的极化SAR相干斑噪声抑制方法 |
CN103020925A (zh) * | 2011-09-27 | 2013-04-03 | 上海微电子装备有限公司 | 模板图像的获取方法 |
CN103077506A (zh) * | 2013-03-06 | 2013-05-01 | 西安电子科技大学 | 结合局部和非局部的自适应图像去噪方法 |
-
2014
- 2014-08-28 CN CN201410431699.7A patent/CN104200434B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6718068B1 (en) * | 2000-03-10 | 2004-04-06 | Eastman Kodak Company | Noise reduction method utilizing statistical weighting, apparatus, and program for digital image processing |
CN101308573A (zh) * | 2008-06-30 | 2008-11-19 | 北京中星微电子有限公司 | 一种消除噪声的方法和装置 |
CN101655977A (zh) * | 2009-09-17 | 2010-02-24 | 浙江大学 | 基于差分图像检测、多窗口滤波去除图像脉冲噪声的方法 |
CN101661611A (zh) * | 2009-09-25 | 2010-03-03 | 西安电子科技大学 | 基于贝叶斯非局部均值滤波器的实现方法 |
CN103020925A (zh) * | 2011-09-27 | 2013-04-03 | 上海微电子装备有限公司 | 模板图像的获取方法 |
CN103020919A (zh) * | 2013-01-09 | 2013-04-03 | 西安电子科技大学 | 基于非局部Lee的极化SAR相干斑噪声抑制方法 |
CN103077506A (zh) * | 2013-03-06 | 2013-05-01 | 西安电子科技大学 | 结合局部和非局部的自适应图像去噪方法 |
Non-Patent Citations (1)
Title |
---|
NSCT De-noising Algorithm with Image Partition and;Yue Liu et al;《Journal of Information & Computational Science》;20121115;第4075-4082页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104200434A (zh) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | An experimental-based review of image enhancement and image restoration methods for underwater imaging | |
CN104200434B (zh) | 一种基于噪声方差估计的非局部均值图像去噪方法 | |
CN109767439B (zh) | 一种自适应窗口的多尺度差异与双边滤波的目标检测方法 | |
CN112819772B (zh) | 一种高精度快速图形检测识别方法 | |
CN107301661A (zh) | 基于边缘点特征的高分辨率遥感图像配准方法 | |
CN103440644B (zh) | 一种基于最小描述长度的多尺度图像弱边缘检测方法 | |
CN108846844B (zh) | 一种基于海天线的海面目标检测方法 | |
CN108932699B (zh) | 基于变换域的三维匹配调和滤波图像去噪方法 | |
CN102930576A (zh) | 一种基于特征流的抽象线条画生成方法 | |
CN103839234B (zh) | 一种基于可控核的双几何非局部均值图像去噪方法 | |
CN111783583B (zh) | 基于非局部均值算法的sar图像相干斑抑制方法 | |
CN103679173A (zh) | 图像显著区域检测方法 | |
CN112801141B (zh) | 基于模板匹配和孪生神经网络优化的异源图像匹配方法 | |
CN108257153B (zh) | 一种基于方向梯度统计特征的目标跟踪方法 | |
CN104599288A (zh) | 一种基于肤色模板的特征跟踪方法及装置 | |
CN113344810A (zh) | 基于动态数据分布的图像增强方法 | |
CN110766657A (zh) | 一种激光干扰图像质量评价方法 | |
CN113822352A (zh) | 基于多特征融合的红外弱小目标检测方法 | |
CN107944497A (zh) | 基于主成分分析的图像块相似性度量方法 | |
CN103268482A (zh) | 一种低复杂度的手势提取和手势深度获取方法 | |
CN103903228A (zh) | 一种基于hwd变换的非局部图像去噪方法 | |
Gupta et al. | A noise robust edge detector for color images using hilbert transform | |
CN105303538B (zh) | 一种基于nsct和pca的高斯噪声方差估计方法 | |
CN104217430A (zh) | 基于l1正则化的图像显著性检测方法 | |
CN107945119B (zh) | 基于拜耳模式的图像内相关噪声估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170412 |
|
CF01 | Termination of patent right due to non-payment of annual fee |