空间域图像增强：二阶导数图像边缘检测算子-Laplacian算子

0.综述

前面介绍了Sobel等几种一阶导数图像边缘检测算子，图像通过与这些算子的卷积计算可以提取出图像边缘信息。但是在利用一阶导数的边缘检测算子进行边缘检测时，有时会出现因检测到的边缘点过多而导致边缘(线)过粗的情况。

通过去除一阶导数中的非局部最大值就可以检测出更细的边缘，而一阶导数的局部最大值对应着二阶导数的零交叉点。所以通过找图像的二阶导数的零交叉点就能找到精确的边缘点。二阶导数可以通过快速的图像像素值强度的变化来检测图像边缘，其检测原理跟一阶导数有点类似，只是二阶导数计算的是图像 $x，y$ 两个方向的二阶导数。

1.Laplacian算子原理分析

与梯度方法类似，用于图像锐化和图像边缘检测的Laplacian算子的概念也是相同的。Laplacian算子是一种无方向性的二阶导数算子，其在点$(x+y)$处的拉普拉斯值定义为:
$${\nabla }^{2}f=\frac{{\partial }^{2}f}{\partial x^2}+\frac{{\partial }^{2}f}{\partial y^2}\text{\hspace{1em}(1)}$$
其中：
$$\frac{{\partial }^{2}f}{\partial x^2}=f(i+1,j)-2f(i,j)+f(i-1,j)\text{\hspace{1em}(2)}$$
$$\frac{{\partial }^{2}f}{\partial y^2}=f(i,j+1)-2f(i,j)+f(i,j-1)\text{\hspace{1em}(3)}$$

合并式$(2)$和式$(3)$可得到实现Laplacian运算的几种模板如下所示：

$$H_1=\left[\begin{array}{ccc}0& -1& 0\\ -1& 4& -1\\ 0& -1& 0\end{array}\right]$$
$$H_2=\left[\begin{array}{ccc}-1& -1& -1\\ -1& 8& -1\\ -1& -1& -1\end{array}\right]$$

Laplacian运算由于只需要用一个模板，所以计算量比较小。Laplacian算子对图像中的噪声比较敏感，并且由于通常产生的是2个像素宽的边缘,所以很少直接用于边缘检测，通常是在已知边缘像素后用于确定该像素是在图像的暗区那边还是在明区那边，一般先对图像进行一次高斯模糊处理，然后在使用Laplacian算子来提取图像边缘(LOG算子，由Marr提出，因此也称为Marr边缘检测算法)。但由于拉普拉斯算子是二阶算子,所以利用其过零点的特性可精确定位边缘。Laplacian二阶边缘检测算子具有各向同性和旋转不变性，是一个标量算子。

2.Marr边缘检测算法

在前面的内容中我们提到了Marr边缘检测算法，该算法克服了一般微分运算对于噪声敏感的缺点，利用能够反映人眼视觉特征的LOG算子对图像进行边缘检测，并结合二阶导数零交叉的性质对边缘进行定位。

典型的二维高斯函数的形式为：
$$G(x，y，\sigma )=-\frac{1}{2\pi {\sigma }^2}{e}^{-\frac{x^2+y^2}{2{\sigma }^2}}$$
其中 $\sigma$ 称为尺度因子，用于控制去噪效果，实验结果表明当 $\sigma =1$时图像去噪效果最好。

由上面的内容也可以知道，Marr算法分为两个主要过程：

• 首先利用二维高斯函数对图像进行低通滤波，也就是用二维高斯函数与图像进行卷积：
  $$g(x,y)=G(x,y,\sigma )\ast f(x,y)$$
  可以得到平滑后的图像$g_0(x,y)$，这一步可以消除图像空间中尺度小于$\sigma$的图像强度变化，从而去除部分噪声。
• 其次使用Laplacian算子对$g_0(x,y)$进行二阶导数运算：
  $${\nabla }^2{g}_0(x,y)={\nabla }^2(G(x,y,\sigma )\ast f(x,y))={\nabla }^{2}G(x,y,\sigma )\ast f(x,y)$$
  其中：
  $${\nabla }^{2}G(x,y,\sigma )=\frac{\partial G}{\partial x^2}+\frac{\partial G}{\partial y^2}=\frac{1}{\pi {\sigma }^4}(\frac{x^2+y^2}{2{\sigma }^2}-1)e^{-\frac{x^2+y^2}{2{\sigma }^2}}$$
  称为高斯-拉普拉斯算子(LOG)算子。由此可见，Marr算法将Gauss的平滑滤波器和Lapla的锐化滤波器结合了起来，也即把二维高斯函数与原图像$f(x,y)$进行卷积核用拉普拉斯算子对卷积结果进行二阶导数运算结合起来，先平滑噪声，在进行边缘检测。

3.OpenCV中的Laplacian算子API

void Laplacian( InputArray src, 
				OutputArray dst, 
				int ddepth,
                int ksize = 1, 
                double scale = 1, 
                double delta = 0,
                int borderType = BORDER_DEFAULT );

其中ddepth一般设为-1表示输入输出图像类型一致。ksize必须为奇数，默认为1，表示是四领域算子，大于1则是八领域算子。

4.代码实践Laplacian算子

Mat blur_result, lap_result;
GaussianBlur(input, blur_result, Size(3, 3), 0);
Laplacian(blur_result, lap_result, CV_32F, 1, 1.0, 127.0, 4);
convertScaleAbs(lap_result, lap_result);

imshow("blur result", blur_result);
imshow("lap result", lap_result);

运行结果如下：