(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 (43)申请公布日 (21)申请号 201710118180.7
(22)申请日 2017.03.01
(71)申请人 武汉科技大学
地址 430000 湖北省武汉市青山区平大道
(72)发明人 伍世虔 何松 陈鹏
(74)专利代理机构 北京科亿知识产权代理事务
所(普通合伙) 11350
代理人 傅海鹏
(51)Int.Cl.
G06T 5/00(2006.01)
(54)发明名称
法
(57)摘要
本发明公开了一种基于Hessian矩阵的图像
噪点检测及去噪方法，针对椒盐噪声利用噪声
点，边缘点和平滑区域点的Hessain矩阵的特征
值的不同特点来检测出噪声点，对检测出来的点
以其为中心的3*3窗口中利用中值滤波的思想，
用窗口的中心值来代替噪声点，对于其他点不做
任何处理；同时在算法评价上提出噪点检测率的
概念；由于在检测噪点时可能会将平滑区域点或
边缘点当成噪点，故设定第三个判定条件来进一
步提高检测准确率和去噪效果；而对于噪声密度
很大的图像本发明通过多次迭代也能够得到较
好的效果；本发明具有更好去噪效果，同时保留
了更多的图像边缘和细节信息；在噪声密度很大
的情况下相比于中值滤波也有较好的效果。权利要求书3页 说明书8页 附图2页CN 107038688 A 2017.08.11
C N 107038688
A
1.一种基于Hessian矩阵的图像噪点检测及去噪方法，其特征在于，包括以下步骤：
1)建立噪声模型
1.1)建立椒盐噪声的理论模型为式(1)：
其中，φ是图像像素点受椒盐噪声污染的概率，f(x，y)、g(x，y)分别是原图像和被污染后图像像素点的灰度值；由式(1)可知，椒盐噪声只会污染图像中的部分像素点，而其它像素的灰度值则保持不变。
1.2)利用噪声点的Hessian矩阵特点检测出噪声点，分析图像L采用式(2)的方法将其在的邻域内进行泰勒展开：
其中H o，s分别表示在尺度s下处的梯度向量和Hessian矩阵，Hessian矩阵由二维
图像L的偏导数组成，如式(3)：
1.3)根据线性尺度空间理论，尺度空间图像偏导数为原图像与高斯函数偏导数的卷积，方程式为式(4)：
其中G(x,y,s)为二维高斯函数，s为标准差，定义如式(5)：
1.4)假设二维Hessian矩阵H的两个特征值为λ1、λ2，由于L xy＝L yx，即H为实对称矩阵，因此Hessian的两个特征值由式(6)和式(7)计算得出：
其中，M＝(L xx+L yy)/2 (8)；
2)根据椒盐噪声和图像边缘分析其Hessian矩阵特征值的不同进行去噪，通过建立多元函数二阶导数组成的Hessian矩阵，计算Hessian矩阵特征值，分析并确定图像最小曲率的方向；
3)利用Hessian矩阵特征值判定边缘点和噪声点；
4)利用中值滤波原理将数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替，让周围像素值接近真实值，消除孤立的噪声点，数学表达式如式(10)：
其中，f(x，y)为滤波输出值，为窗口W的中心值，代表窗口中的所有点；g(x,y)代表点(x,y)处的灰度值；
5)算法分析，将步骤3)中判定出的噪声点，按式(11)、式(12)的判定条件判定出椒盐噪声点：
(λ1＞T1)&&(λ2＞T1) (11)；
λ1×λ2＞0 (12)；
以中值滤波采用3*3的窗口，定义第三个判定条件提高噪声点的检测准确率和去噪效果，如式(13)：
其中，λ1、λ2为Hessian矩阵的两个特征值，T1,T2为判定阈值；
通过式(11)-(13)的判定条件筛选得到检测到的椒盐噪声点，对检测到的噪声点运用中值滤波原理，以椒盐噪声点为中心选用3*3窗口，用窗口中所有点的中值来代替噪声点；
6)算法评价，采用均方误差和峰值信噪比对处理后的图像进行算法评价，均方误差值MSE越小，处理后的图像质量越高，峰值信噪比PSNR值越大，处理后的图像视觉效果越好；
尺寸为R×C的图像峰值信噪比算法如式(14)：
其中，f(x，y)和f(x，y)分别代表原始图像和处理后的图像；
7)噪声点检测评价，采用式(15)的算法计算噪声点检测率：
M＝φ×I
其中，I为原图像像素点总数，为噪声点污染概率，M为噪声点总数，N为处理后噪声点总
数，N的值由式噪声理论模型得知处理后的图像中0和255的点近似认为未处理的噪声污染点，原图像中存在0和255的点，则噪声检测率的实验结果偏小于实际值的，噪声点检测率的值越大说明检测出的噪点越多，去噪效果越好。
2.根据权利要求1所述的基于Hessian矩阵的图像噪点检测及去噪方法，其特征在于，所述步骤3)中，边缘点和噪声点的判定依据点的类型与特征值的关系为：边缘点-L、H；盐噪声-H-、H-；胡椒噪声-H+、H+；平滑区域点-L、L，其中L＝Low，H＝High代表特征值大小，±代表特征值的符号。
3.根据权利要求1所述的基于Hessian矩阵的图像噪点检测及去噪方法，其特征在于，所述步骤4)中，中值滤波窗口的大小为3*3、5*5或7*7。
4.根据权利要求1所述的基于Hessian矩阵的图像噪点检测及去噪方法，其特征在于，所述步骤5)中，当椒盐噪声点密度大于0.1时，采用多次迭代处理椒盐噪声点，当处理后的图像峰值信噪比小于上次迭代处理结果则终止迭代，实现去噪。
基于Hessian矩阵的图像噪点检测及去噪方法技术领域
[0001]本发明涉及图像处理技术领域，具体是一种基于Hessian矩阵的图像噪点检测及去噪方法。
背景技术
[0002]图像是人们获取外界信息的主要来源，在数字图像中，噪声主要来源于图像的获取和传输过程，噪声不仅影响着图像的主观视觉效果，且其他图像处理操作的好坏也依赖于前期噪声处理的好坏，因此图像去噪在图像预处理中是非常重要的，是图像处理中的一个重要环节。
[0003]图像在编码和传输过程中影响着图像质量的主要因素是脉冲噪声，脉冲噪声的去除得到了广泛的研究和快速的发展，其中效果较好的为中值滤波，但传统的中值滤波是将整幅图像的全部像素点用邻域窗口的中心值代替，虽然具备较好的去噪效果，但是对图像的全部像素点都进行了处理，故而模糊了图像的边缘、细线等重要细节信息；同时当脉冲噪声密度很大时，中值滤波处理噪声效果则达不到所需的要求；为了解决这些问题，近年来已经提出了各种各样改进中值滤波算法，如加权中值滤波算法、开关中值滤波算法、自适应中值滤波算法、minmax算法、极值中值滤波算法等，但都或多或少的存在缺陷，基于上述原因，需要提出一种新的噪点检测和去噪方法来满足使用需求。
发明内容
[0004]本发明的目的在于提供一种在去噪点同时保护图像边缘和细节、提高峰值信噪比、降低均方误差的基于Hessian矩阵的图像噪点检测及去噪方法，以解决上述背景技术中提出的问题。
[0005]为实现上述目的，本发明提供如下技术方案：
[0006]一种基于Hessian矩阵的图像噪点检测及去噪方法，包括以下步骤：
[0007]1)建立噪声模型
[0008] 1.1)建立椒盐噪声的理论模型为式(1)：
[0009]
[0010]其中，φ是图像像素点受椒盐噪声污染的概率，f(x，y)、g(x，y)分别是原图像和被污染后图像像素点的灰度值；由式(1)可知，椒盐噪声只会污染图像中的部分像素点，而其它像素的灰度值则保持不变。
[0011] 1.2)利用噪声点的Hessian矩阵特点检测出噪声点，分析图像L采用式(2)的方法将其在的邻域内
进行泰勒展开：
说 明 书
1/8页CN 107038688 A
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议