Canny边缘检测是一种经典的边缘检测算法,由John F.
Canny在1986年提出。它被广泛应用于计算机视觉和图像处理领域,是一种多阶段的边缘检测算法,能够有效地检测图像中的边缘并抑制噪声。
Canny边缘检测的主要步骤如下:
噪声抑制:首先,通过使用高斯滤波器对图像进行平滑处理,以去除图像中的噪声。高斯滤波器可以有效地平滑图像,同时保持边缘的细节。
计算梯度幅值和方向:使用Sobel算子计算图像中每个像素点的水平和垂直方向的梯度值。然后,根据梯度值计算每个像素点的梯度幅值和方向。
非极大值抑制:在计算得到的梯度幅值图像上进行非极大值抑制。这一步的目的是将边缘宽度变窄,使得边缘更加细化和明确。
双阈值处理:根据设定的高阈值和低阈值,将梯度幅值图像中的像素点分为强边缘、弱边缘和非边缘三类。通常选择高阈值和低阈值使得强边缘像素点的梯度幅值大于高阈值,非边缘像素点的梯度幅值小于低阈值,而弱边缘像素点的梯度幅值处于高阈值和低阈值之间。
边缘连接:最后,通过连接强边缘像素点和与之相邻的弱边缘像素点,得到完整的边缘图像。
Canny边缘检测算法通过多个阶段的处理,能够得到清晰准确的边缘信息,并且对噪声具有一定的鲁棒性。因此,它在图像处理和计算机视觉中得到广泛应用,特别是在要求高精度边缘检测的场景中。
Canny边缘检测在图像处理和计算机视觉领域有许多应用场景,下面列举了一些常见的应用场景:
物体检测与目标定位:Canny边缘检测能够帮助检测图像中物体的边缘,从而实现目标检测和定位。在计算机视觉任务中,这对于目标识别、物体追踪和目标定位等是至关重要的。
图像分割:Canny边缘检测可以在图像中检测出物体和背景之间的边缘,有助于将图像分割成不同的区域,使得图像处理更加高效和准确。
视觉导航与SLAM:在机器人视觉导航和同时定位与地图构建(SLAM)中,Canny边缘检测有助于提取环境中的地标和边缘特征,用于机器人的定位和导航。
图像增强:Canny边缘检测可以突出图像中的边缘特征,使得图像在可视化和分析上更加清晰明了,从而用于图像增强和美化。
图像匹配与对准:Canny边缘检测能够提取图像中的特征点,用于图像匹配和图像对准,常用于图像拼接、图像融合等应用。
视觉检测与安全:Canny边缘检测在视觉检测和安全领域也有应用,例如边缘检测在视频监控中用于检测异常行为,或者在车辆驾驶辅助系统中用于车道检测和车辆识别。
医学影像处理:Canny边缘检测在医学影像处理中广泛应用,用于检测器官边缘、病变区域等,辅助医生进行疾病诊断和治疗。
总的来说,Canny边缘检测在图像处理和计算机视觉的许多领域都扮演着重要的角色,它是一种经典且有效的边缘检测算法,被广泛应用于实际场景中。
由于图像边缘非常容易受到噪声的干扰,因此为了避免检测到错误的边缘信息,通常需要对图像进行滤波以去除噪声。滤波的目的是平滑一些纹理较弱的非边缘区域,以便得到更准确的边缘。在实际处理过程中,通常采用高斯滤波去除图像中的噪声。
图 10-1 演示了使用高斯滤波器 T 对原始图像 O 中像素值为 226 的像素点进行滤波,得到该点在滤波结果图像 D 内的值的过程。
在滤波过程中,我们通过滤波器对像素点周围的像素计算加权平均值,获取最终滤波结果。对于高斯滤波器 T,越临近中心的点,权值越大。在图 10-1 中,对图像 O 中像素值为 226 的像素点,使用滤波器 T 进行滤波的计算过程及结果为:
结果 = 156×(197×1+25×1+106×2+156×1+159×1
+149×1+40×3+107×4+5×3+71×1
+163×2+198×4+226×8+223×4+156×2
+222×1+37×3+68×4+193×3+157×1
+42×1+72×1+250×2+41×1+75×1)
= 138
当然,高斯滤波器(高斯核)并不是固定的,例如它还可以是:
滤波器的大小也是可变的,高斯核的大小对于边缘检测的效果具有很重要的作用。滤波器的核越大,边缘信息对于噪声的敏感度就越低。不过,核越大,边缘检测的定位错误也会随之增加。通常来说,一个 5×5 的核能够满足大多数的情况。
前面一节讲了如何计算图像梯度的幅度。在这里,我们关注梯度的方向,梯度的方向与边缘的方向是垂直的。
边缘检测算子返回水平方向的Gx和垂直方向的Gy。梯度的幅度和方向(用角度值表示)为:
式中,atan2(•)表示具有两个参数的 arctan 函数。
梯度的方向总是与边缘垂直的,通常就近取值为水平(左、右)、垂直(上、下)、对角线(右上、左上、左下、右下)等 8 个不同的方向。
因此,在计算梯度时,我们会得到梯度的幅度和角度(代表梯度的方向)两个值。
图 10-2 展示了梯度的表示法。其中,每一个梯度包含幅度和角度两个不同的值。为了方便观察,这里使用了可视化表示方法。例如,左上角顶点的值“2↑”实际上表示的是一个二元数对“(2, 90)”,表示梯度的幅度为 2,角度为 90°。
在获得了梯度的幅度和方向后,遍历图像中的像素点,去除所有非边缘的点。
在具体实现时,逐一遍历像素点,判断当前像素点是否是周围像素点中具有相同梯度方向的最大值,并根据判断结果决定是否抑制该点。
通过以上描述可知,该步骤是边缘细化的过程。针对每一个像素点:
在图 10-3 中,A、B、C 三点具有相同的方向(梯度方向垂直于边缘)。判断这三个点是否为各自的局部最大值:如果是,则保留该点;否则,抑制该点(归零)。
经过比较判断可知,A 点具有最大的局部值,所以保留 A 点(称为边缘),其余两点(B和 C)被抑制(归零)。
在图 10-4 中,黑色背景的点都是向上方向梯度(水平边缘)的局部最大值。因此,这些点会被保留;其余点被抑制(处理为 0)。这意味着,这些黑色背景的点最终会被处理为边缘点,而其他点都被处理为非边缘点。
“正/负梯度方向上”是指相反方向的梯度方向。例如,在图 10-5 中,黑色背景的像素点都是垂直方向梯度(向上、向下)方向上(即水平边缘)的局部最大值。这些点最终会被处理为边缘点。
经过上述处理后,对于同一个方向的若干个边缘点,基本上仅保留了一个,因此实现了边缘细化的目的。
完成上述步骤后,图像内的强边缘已经在当前获取的边缘图像内。但是,一些虚边缘可能也在边缘图像内。这些虚边缘可能是真实图像产生的,也可能是由于噪声所产生的。对于后者,必须将其剔除。
设置两个阈值,其中一个为高阈值 maxVal,另一个为低阈值 minVal。根据当前边缘像素的梯度值(指的是梯度幅度,下同)与这两个阈值之间的关系,判断边缘的属性。具体步骤为:
(1)如果当前边缘像素的梯度值大于或等于 maxVal,则将当前边缘像素标记为强边缘。
(2)如果当前边缘像素的梯度值介于 maxVal 与 minVal 之间,则将当前边缘像素标记为虚
边缘(需要保留)。
(3)如果当前边缘像素的梯度值小于或等于 minVal,则抑制当前边缘像素。
在上述过程中,我们得到了虚边缘,需要对其做进一步处理。一般通过判断虚边缘与强边
缘是否连接,来确定虚边缘到底属于哪种情况。通常情况下,如果一个虚边缘:
在图 10-6 中,左图显示的是三个边缘信息,右图是对边缘信息进行分类的示意图,具体划分如下:
图 10-7 显示了对图 10-6 中的虚边缘 B 和 C 的处理结果。其中:
注意,高阈值 maxVal 和低阈值 minVal 不是固定的,需要针对不同的图像进行定义。
10-8 给出了一个 Canny 边缘检测的效果图。
Canny边缘检测的边缘连接是指将非极大值抑制得到的边缘点连接成连续的边缘线。这个过程是Canny边缘检测算法中的最后一步,它的目的是去除由于非极大值抑制产生的间断的边缘点,从而得到更加准确和连续的边缘检测结果。
下面详细讲解Canny边缘检测的边缘连接过程:
遍历梯度幅值图像:首先,遍历经过非极大值抑制后的梯度幅值图像,即只有边缘上的像素点的梯度幅值被保留,其他像素点的梯度幅值为零。
标记边缘点:对于每个强边缘像素点(梯度幅值大于高阈值),将其标记为边缘点。强边缘像素点是图像中梯度值较大的像素点,它们代表了图像中明显的边缘。
边缘连接:对于每个边缘点的相邻像素点,如果其梯度幅值大于低阈值,并且没有被标记为边缘点,则将其标记为弱边缘点,并递归地进行边缘连接。这一步的目的是将与强边缘像素点相邻的弱边缘像素点连接到边缘线上。
递归连接:在边缘连接过程中,如果某个弱边缘像素点被标记为边缘点,则会继续检查该像素点的相邻像素点,以便将所有与强边缘像素点相邻的弱边缘像素点连接到边缘线上。
结束条件:边缘连接的递归过程会一直进行,直到所有与强边缘像素点相邻的弱边缘像素点都被标记为边缘点,没有更多的像素点可以连接。
非边缘点处理:经过边缘连接后,所有未被标记为边缘点的像素点视为非边缘点,并抑制其梯度幅值为零,从而得到最终的边缘图像。
OpenCV 提供了函数 cv2.Canny()来实现 Canny 边缘检测,其语法形式如下:
edges = cv.Canny( image, threshold1, threshold2[, apertureSize[, L2gradient]])
式中:
代码如下:
import cv2
#读取图像,灰度模式
o=cv2.imread("feibu.png",cv2.IMREAD_GRAYSCALE)
#边缘检测,canny算子,阈值128-200,低于128的像素点认为是边缘,高于200的像素点认为是边缘,中间值的像素点如果与边缘点相连则认为是边缘
r1=cv2.Canny(o,128,200)
cv2.imshow("original",o)
cv2.imshow("result1",r1)
cv2.waitKey()
cv2.destroyAllWindows()
运行效果: