CN104915938B - 一种高密度椒盐噪声污染图像修复方法 - Google Patents
一种高密度椒盐噪声污染图像修复方法 Download PDFInfo
- Publication number
- CN104915938B CN104915938B CN201510377191.8A CN201510377191A CN104915938B CN 104915938 B CN104915938 B CN 104915938B CN 201510377191 A CN201510377191 A CN 201510377191A CN 104915938 B CN104915938 B CN 104915938B
- Authority
- CN
- China
- Prior art keywords
- pixel
- uncontaminated
- image
- pollution
- value
- 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.)
- Active
Links
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
为提高高密度椒盐噪声污染图像信噪比,本发明公开了一种高密度椒盐噪声污染图像修复方法。该方法将污染图像像素标记为污染像素和未污染像素,在图像修复过程中,未污染像素的像素值保持不变,污染像素的像素值由未污染像素在基于图像空间距离最小二乘回归模型下的估计值代替。本发明的优点是:可显著提高修复图像信噪比,增加修复图像的视觉可视信息,特别是在高密度(比如99%)椒盐噪声污染情况下,仍可得到图像基本形貌信息。
Description
技术领域
本发明涉及数字图像处理技术领域,特指一种高密度椒盐噪声污染图像修复方法。
背景技术
数字成像传感器CCD或CMOS等被广泛应用于工业、娱乐、民用等领域,在实际使用过程中,受制造缺陷、器件老化等因素影响,在获得的成像图像中存在着椒盐噪声污染。椒盐噪声常表现为稳定的极亮或极暗像素,以取值范围为0-255的灰度图像为例,像素值通常为255或0。
目前,为消除椒盐噪声污染,通常采用中值滤波对噪声图像进行处理。但是,传统的中值滤波,会将未污染像素错误地用邻域内像素的中值取代,而产生失真结果。为此,近期学者们提出了Switching Median Filter(SMF)。该方法的基本思想是:先在受污染图像中标记污染像素和未污染像素,在进行中值滤波时,只对污染像素进行更新。具体更新方法是,以污染像素邻域内未污染像素的中值取代污染像素,未污染像素不作更新。这样可保证未污染像素不被邻域内像素中值替代,具有更好保真度。但是,当椒盐噪声污染非常严重时,比如噪声密度高达90%以上(噪声密度是指:被污染像素占整幅图像像素的比例)时,SMF方法难以得到较理想的修复图像。
发明内容
为克服现有技术存在的缺陷,本发明提供一种高密度椒盐噪声污染图像修复方法。该方法将污染图像像素标记为污染像素和未污染像素,在图像修复过程中,未污染像素的像素值保持不变,污染像素的像素值由未污染像素在基于图像空间距离最小二乘回归模型下的估计值代替。
具体地,本发明一种高密度椒盐噪声污染图像修复方法,其步骤如下:
(1)遍历图像中所有像素,判定其像素值是否为0或255,若是,则将其标记为污染像素E,若不是,则标记为未污染像素T;
给定椒盐噪声污染图像I,I中任意像素可表示为p(u,v,q),u,v分别为图像空间坐标,q为图像像素值,q的取值范围为[0,255],将图像中像素值q为0或255的像素标记为污染像素,记作E={e1,…},取值为0-255之间的像素为未污染像素,记作T={t1,…};
(2)基于所有未污染像素T,计算基于图像空间距离的最小二乘回归模型的回归系数A*;
根据图像的马尔科夫性,未污染图像中的相邻像素像素值q存在强关联性,任意一个像素p,可由其邻近像素G={g1,...,gm}线性表示:
其中,qp是像素p的像素值,是第i个邻近像素的像素值,是像素p与第i个邻近像素的图像空间距离,wi是与第i个邻近像素的加权权重;
加权权重wi与邻近像素gi对应,可将两者乘积看作线性回归系数,那么式(1) 可重新写为下面的线性方程:
根据式(2),当线性回归系数ai已知时,利用污染像素与未污染像素之间的图像空间距离dj,即可估计出污染像素的像素值:
其中,aj是与第j个未污染像素tj对应的线性回归系数,是污染像素e与第j个未污染像素tj的图像空间距离。
假定未污染像素数量为k,那么回归系数A=[a1,...,ak]T可由这些未污染像素学习得到。k 个未污染像素t的线性回归方程可表示为:
式(4)写成矩阵形式为:
Q=DA (5)
其中,
Q和D已知,需要估计回归系数A。估计的最优回归系数A*应使现有的未污染像素的线性回归误差最小,即A*满足最小二乘条件:
式(6)的解为:
A*=(DDT)-1DTQ (7)
(3)遍历所有污染像素E,计算污染像素e与所有未污染像素间的空间距离d=[d1,...,dk],其中代入式(8)的回归预测模型,得到其估计值
(4)用步骤(3)得到的污染像素的像素估计值代替污染像素原像素值,未污染像素的像素值保持不变。
最小二乘回归是统计学中常用的统计分析方法,在该模型下,给定一组输入、输出数据,通过最小二乘约束,可以估计出预测模型的相关系数参数。本发明将在该框架下,在现有SMF 方法基础上,提出基于像素图像空间距离的最小二乘回归模型,用于估计污染像素的像素值。采用最小二乘回归模型的估计值代替SMF方法中的邻域内未污染像素的中值,用于污染像素的修复,可进一步提高修复图像的信噪比。
本发明公开的高密度椒盐噪声污染图像修复方法的优点是:可显著提高修复图像信噪比,增强修复图像的视觉可视信息,特别是在高密度(比如99%)椒盐噪声污染情况下,仍可得到图像基本形貌信息。以lena图像为例进行了仿真测试,与SMF算法修复效果进行对比的结果如图1所示,其中1是原始输入图像,2-4分别是椒盐噪声污染污染图像,噪声密度(噪声密度是指污染像素占所有像素的比率)依次为90%,95%,99%。5-7是SMF滤波在7*7搜索窗口大小情况下的修复结果,8-10是本发明方法修复结果。5-7的峰值信噪比(PSNR)为:24.47,22.11,12.31;8-10的峰值信噪比(PSNR)为:26.41,24.18,19.91。
附图说明
图1为分别采用本发明提供的方法、SMF方法进行高密度椒盐噪声去噪结果对比结果图;
图2为图像像素三维空间分布图;
图3为基于图像空间距离加权的图像像素值估计示意图;
图4为基于未污染像素加权表示的污染像素值估计示意图。
具体实施方式
为了使本发明的目的、技术方案及有益效果更加清楚明白,下面结合附图及实施例,对本发明进行进一步详细说明。应当注意,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
下面以8bits灰度图像为例,给出本发明的具体实施例。
给定椒盐噪声污染图像I,I中任意像素可表示为p(u,v,q),u,v分别为图像空间坐标,q为图像像素值,取值范围为[0,255]。图像椒盐噪声通常表现为像素值q为0或255的像素。因此,图像中像素值q为0或255的像素标记为污染像素,记作E={e1,…},取值为0-255之间的像素为未污染像素,记作T={t1,…}。如图2所示,14是图像横坐标u,15是图像纵坐标v, 16是图像像素值坐标q,11是图像I的三维曲面,12是未污染像素t,13是污染像素e,17是污染像素图像空间坐标,18是未污染像素图像空间坐标。本发明基本思想是:根据未污染像素12,训练基于图像空间距离的最小二乘回归模型,根据污染像素13的图像空间坐标,估计其像素值,对其进行修复。
根据图像的马尔科夫性,未污染图像中的相邻像素像素值q存在强关联性,任意一个像素p,可由其邻近像素G={g1,...,gm}线性表示:
其中,qp是像素p的像素值,是第i个邻近像素的像素值,是像素p与第i个邻近像素的图像空间距离(这里取欧式距离),wi是与第i个邻近像素的加权权重。
图3给出了一个例子,以像素19为例,它的像素值可由与之邻近的4个像素20-23的像素值加权表示,其中24是像素19与像素21之间的图像空间距离。
加权权重wi与邻近像素gi对应,可将两者乘积 看作 线性回归系数,那么式(1) 可重新写为 下面的线性方程:
根据式2,当线性回归系数ai已知时,利用污染像素与未污染像素之间的图像空间距离dj,即可估计出污染像素的像素值:
其中,aj是与第j个未污染像素tj对应的线性回归系数,是污染像素e与第j个未污染像素tj的图像空间距离。
图4给出了一个示例,25是污染像素,20-23是与之相邻的四个未污染像素,26是其中一个未污染像素21与污染像素25之间的图像空间距离,那么污染像素25的像素估计值,可根据25与邻近像素20-23的距离,通过式(3)进行加权表示。
假定未污染像素数量为k,那么回归系数A=[a1,...,ak]T可由这些未污染像素学习得到。k 个未污染像素t的线性回归方程可表示为:
式(4)写成矩阵形式为:
Q=DA (5)
其中,
Q和D已知,需要估计回归系数A。估计的最优回归系数A*应使现有的未污染像素的线性回归误差最小,即A*满足最小二乘条件:
式(6)的解为:
A*=(DDT)-1DTQ (7)
因为,D是像素间图像空间距离,因此本发明的算法叫作“基于图像空间距离的最小二乘回归方法”。得到线性回归系数A*后,给定某个污染像素e与未污染像素间的图像空间距离d=[d1,...,dk],即可通过式(8),估计出该污染像素e的像素值
上面,介绍了本发明的基本技术原理,下面给出具体操作步骤:
1)遍历图像中所有像素,判定其像素值是否为0或255,若是,则将其标记为污染像素 E,若不是,则标记为未污染像素T;
2)基于所有未污染像素T,按式(7),计算基于图像空间距离的最小二乘回归模型的回归系数A*;
3)遍历所有污染像素E,计算污染像素e与所有未污染像素间的空间距离d=[d1,...,dk],其中代入式(8)的回归预测模型,得到其估计值
4)用步骤3得到的污染像素的像素估计值代替污染像素原像素值,未污染像素的像素值保持不变。
Claims (1)
1.一种高密度椒盐噪声污染图像修复方法,其特征在于:包括以下步骤:
(1)遍历图像中所有像素,判定其像素值是否为0或255,若是,则将其标记为污染像素E,若不是,则标记为未污染像素T;
给定椒盐噪声污染图像I,I中任意像素可表示为p(u,v,q),u,v分别为图像空间坐标,q为图像像素值,q的取值范围为[0,255],将图像中像素值q为0或255的像素标记为污染像素,记作E={e1,…},取值为0-255之间的像素为未污染像素,记作T={t1,…};
(2)基于所有未污染像素T,计算基于图像空间距离的最小二乘回归模型的回归系数A*;
根据图像的马尔科夫性,未污染图像中的相邻像素像素值q存在强关联性,任意一个像素p,可由其邻近像素G={g1,...,gm}线性表示:
其中,qp是像素p的像素值,是第i个邻近像素的像素值,是像素p与第i个邻近像素的图像空间距离,wi是与第i个邻近像素的加权权重;
加权权重wi与邻近像素gi对应,可将两者乘积看作线性回归系数,那么式(1)可重新写为下面的线性方程:
根据式(2),当线性回归系数ai已知时,利用污染像素与未污染像素之间的图像空间距离dj,即可估计出污染像素的像素值:
其中,aj是与第j个未污染像素tj对应的线性回归系数,是污染像素e与第j个未污染像素tj的图像空间距离;
假定未污染像素数量为k,那么回归系数A=[a1,...,ak]T可由这些未污染像素学习得到;k个未污染像素t的线性回归方程可表示为:
式(4)写成矩阵形式为:
Q=DA(5)
其中,
Q和D已知,需要估计回归系数A;估计的最优回归系数A*应使现有的未污染像素的线性回归误差最小,即A*满足最小二乘条件:
式(6)的解为:
A*=(DDT)-1DTQ (7)
(3)遍历所有污染像素E,计算污染像素e与所有未污染像素间的空间距离d=[d1,...,dk],其中代入式(8)的回归预测模型,得到其估计值
(4)用步骤(3)得到的污染像素的像素估计值代替污染像素原像素值,未污染像素的像素值保持不变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510377191.8A CN104915938B (zh) | 2015-07-02 | 2015-07-02 | 一种高密度椒盐噪声污染图像修复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510377191.8A CN104915938B (zh) | 2015-07-02 | 2015-07-02 | 一种高密度椒盐噪声污染图像修复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104915938A CN104915938A (zh) | 2015-09-16 |
CN104915938B true CN104915938B (zh) | 2017-08-01 |
Family
ID=54084977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510377191.8A Active CN104915938B (zh) | 2015-07-02 | 2015-07-02 | 一种高密度椒盐噪声污染图像修复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104915938B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105719257B (zh) * | 2016-01-28 | 2018-08-03 | 河南师范大学 | 图像超高密度椒盐噪声的降除方法 |
CN106683054B (zh) * | 2016-12-08 | 2019-12-27 | 南京理工大学 | 基于图像抖动技术抑制椒盐噪声的方法 |
CN106851322A (zh) * | 2017-01-15 | 2017-06-13 | 四川精目科技有限公司 | 一种分块最小二乘回归高速相机压缩图像重建方法 |
CN106851321A (zh) * | 2017-01-15 | 2017-06-13 | 四川精目科技有限公司 | 一种最小二乘回归高速相机压缩图像重建方法 |
CN107067376A (zh) * | 2017-01-15 | 2017-08-18 | 四川精目科技有限公司 | 一种rbf插值椒盐噪声图像修复方法 |
CN108320269A (zh) * | 2017-01-18 | 2018-07-24 | 重庆邮电大学 | 一种高效消除高密度椒盐噪声的方法 |
CN107730457B (zh) * | 2017-08-28 | 2020-02-14 | 广东数相智能科技有限公司 | 一种图像补全方法、装置、电子设备及存储介质 |
CN111179190B (zh) * | 2019-12-16 | 2020-12-08 | 郑州轻工业大学 | 一种基于Criminisi的图像复原方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103761707A (zh) * | 2013-12-20 | 2014-04-30 | 浙江大学 | 一种快速高效去除图像椒盐噪声的均值滤波方法 |
CN104167005A (zh) * | 2014-07-07 | 2014-11-26 | 浙江大学 | 一种基于相似度函数的自适应加权椒盐噪声滤波方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101669819B1 (ko) * | 2010-11-11 | 2016-10-27 | 삼성전자주식회사 | 깊이 영상의 고정밀 복원을 위한 필터링 장치 및 방법 |
-
2015
- 2015-07-02 CN CN201510377191.8A patent/CN104915938B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103761707A (zh) * | 2013-12-20 | 2014-04-30 | 浙江大学 | 一种快速高效去除图像椒盐噪声的均值滤波方法 |
CN104167005A (zh) * | 2014-07-07 | 2014-11-26 | 浙江大学 | 一种基于相似度函数的自适应加权椒盐噪声滤波方法 |
Non-Patent Citations (2)
Title |
---|
Spatially Adaptive Total Variation Image Denoising under Salt and Pepper Noise;Renan Rojas等;《19th European Signal Processing Conference》;20110902;第278-282页 * |
基于噪声监测的高密度椒盐噪声自适应滤波算法;周军妮 等;《电视技术》;20141231;第38卷(第19期);第41-46页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104915938A (zh) | 2015-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104915938B (zh) | 一种高密度椒盐噪声污染图像修复方法 | |
CN107194872B (zh) | 基于内容感知深度学习网络的遥感图像超分辨率重建方法 | |
CN105208374B (zh) | 一种基于深度学习的无参考图像质量客观评价方法 | |
CN102722876B (zh) | 一种基于残差的图像超分辨率重建方法 | |
CN102883175B (zh) | 深度图提取、判断视频场景切换及深度图边缘优化方法 | |
CN104217404B (zh) | 雾霾天视频图像清晰化处理方法及其装置 | |
CN104902267B (zh) | 一种基于梯度信息的无参考图像质量评价方法 | |
CN107316295A (zh) | 一种基于深度神经网络的织物瑕疵检测方法 | |
CN111079734B (zh) | 铁路货车三角孔异物检测方法 | |
CN103500449B (zh) | 一种星上可见光遥感图像云检测方法 | |
CN102572499B (zh) | 基于小波变换多分辨率预测的无参考图像质量评价方法 | |
CN105657402A (zh) | 一种深度图恢复方法 | |
CN104303208A (zh) | 用于去除包含在视频中的雾的图像处理装置及其方法 | |
CN110070539A (zh) | 基于信息熵的图像质量评价方法 | |
CN106780727B (zh) | 一种车头检测模型重建方法及装置 | |
CN108022225A (zh) | 基于快速导向滤波改进的暗通道先验图像去雾算法 | |
TWI376648B (en) | Method and device for keeping image background by multiple gauss models | |
CN104504676A (zh) | 一种多视觉敏感特征相似度的全参考图像质量评价方法 | |
CN107389693A (zh) | 一种基于机器视觉的印刷品缺陷自动检测方法 | |
CN106341677A (zh) | 虚拟视点视频质量评价方法 | |
CN105654442A (zh) | 一种冲击噪声图像去噪方法 | |
CN104036493A (zh) | 一种基于多重分形谱的无参考图像质量评价方法 | |
CN116542932A (zh) | 一种基于改进的YOLOv5s的注塑件表面缺陷检测方法 | |
CN115619718A (zh) | 基于嵌入式边缘平台的钢板缺陷实时检测系统与方法 | |
CN102609689B (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 |