[发明专利]基于Hessian矩阵的图像噪点检测及去噪方法

权利要求书

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的两个特征值为λ₁、λ₂，由于L_xy＝L_yx,即H为实对称矩阵，因此Hessian的两个特征值由式(6)和式(7)计算得出：

其中，M＝(L_xx+L_yy)/2(8)；

1

2)根据椒盐噪声和图像边缘分析其Hessian矩阵特征值的不同进行去噪，通过建立多元函数二阶导数组成的Hessian矩阵，计算Hessian矩阵特征值，分析并确定图像最小曲率的方向；

3)利用Hessian矩阵特征值判定边缘点和噪声点；

4)利用中值滤波原理将数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替，让周围像素值接近真实值，消除孤立的噪声点，数学表达式如式(10)：

其中，f(x，y)为滤波输出值，为窗口W的中心值，代表窗口中的所有点；g(x,y)代表点(x,y)处的灰度值；

5)算法分析，将步骤3)中判定出的噪声点，按式(11)、式(12)的判定条件判定出椒盐噪声点：

(λ₁＞T₁)&&(λ₂＞T1)(11)；

λ₁×λ₂＞0(12)；

以中值滤波采用3*3的窗口，定义第三个判定条件提高噪声点的检测准确率和去噪效果，如式(13)：

其中，λ₁、λ₂为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代表特征值大小，±代表特征值的符号。