网站如何被百度快速收录,网站建设 投标书,哈尔滨网站建设托管,网店推广营销方案1 边缘检测介绍
图像边缘检测技术是图像处理和计算机视觉等领域最基本的问题#xff0c;也是经典的技术难题之一。如何快速、精确地提取图像边缘信息#xff0c;一直是国内外的研究热点#xff0c;同时边缘的检测也是图像处理中的一个难题。早期的经典算法包括边缘算子方法…
1 边缘检测介绍
图像边缘检测技术是图像处理和计算机视觉等领域最基本的问题也是经典的技术难题之一。如何快速、精确地提取图像边缘信息一直是国内外的研究热点同时边缘的检测也是图像处理中的一个难题。早期的经典算法包括边缘算子方法、曲面拟合的方法、模板匹配方法、阈值法等。
近年来随着数学理论与人工智能技术的发展出现了许多新的边缘检测方法如Roberts、Laplacan、Canny等图像的边缘检测方法。这些方法的应用对于高水平的特征提取、特征描述、目标识别和图像理解有重大的影响。然而在成像处理的过程中投影、混合、失真和噪声等会导致图像模糊和变形这使得人们一直致力于构造具有良好特性的边缘检测算子。
1.1 什么是边缘检测
边缘检测是图像处理和计算机视觉中的基本问题边缘检测的目的是标识数字图像中亮度变化明显的点。图像属性中的显著变化通常反映了属性的重要事件和变化包括深度不连续、表面方向不连续、物质属性变化和场景照明变化。边缘检测特征是提取中的一个研究领域。图像边缘检测大幅度地减少了数据量并且剔除了可以认为不相关的信息保留了图像重要的结构属性。 1.2 边缘检测的方法
人类视觉系统认识目标的过程分为两步首先把图像边缘与背景分离出来然后到图像的细节辨认出图像的轮廓。计算机视觉正是模仿人类视觉的过程。
因此在检测物体边缘时先对轮廓点进行粗略检测然后通过链接规则把原来检测到的轮廓点连接起来同时检测和连接遗漏的边界点及去除虚假的边界点。图像的边缘是图像的重要特征是计算机视觉、模式识别等的基础因此边缘检测是图像处理中一个重要的环节。然而边缘检测是图像处理中的一个难题因为实际景物图像的边缘往往是各种类型的边缘及它们模糊化后结果的组合且实际图像信号存在噪声。噪声和边缘都属于高频信号很难用频带做取舍。
边缘是指图像周围像素灰度有阶跃变化或屋顶状变化的像素集合存在于目标与背景、目标与目标、区域与区域、基元与基元之间。边缘具有方向和幅度两个特征沿边缘走向像素值变化比较平缓垂直于边缘走向像素值变化比较剧烈可能呈现阶跃状也可能呈现斜坡状。因此边缘可以分为两种 一种为阶跃性边缘两边的像素灰度值有着明显的不同 另一种为屋顶状边缘位于灰度值从增加到减少的变化转折点。
对于阶跃性边缘二阶方向导数在边缘处呈零交叉对于屋顶状边缘二阶方向导数在边缘处取极值。有许多方法可以用于边缘检测绝大部分可以划分为两类基于搜索的一类和基于零穿越的一类。 基于搜索通过寻找图像一阶导数中的最大值来检测边界然后利用计算结果估计边缘的局部方向通常采用梯度的方向并利用此方向找到局部梯度模的最大值代表算法是Sobel算子和Scharr算子。 基于零穿越通过寻找图像二阶导数零穿越来寻找边界代表算法是Laplacian算子。 1.3 典型算子比较
算子优缺点Roberts对具有陡峭的低噪声的图像处理效果较好但利用Roberts算子提取边缘的结果是边缘比较粗因此边缘定位不是很准确Sobel对灰度渐变和噪声较多的图像处理效果比较好Sobel算子对边缘定位比较准确Kirsch对灰度渐变和噪声较多的图像处理效果较好Prewitt 对灰度渐变和噪声较多的图像处理效果较好Laplacian对图像中的阶跃性边缘点定位准确对噪声非常敏感丢失一部分边缘的方向信息造成一些不连续的检测边缘LoGLoG算子经常出现双边缘像素边界而且该检测算法对噪声比较敏感所以很少用LoG算子检测边缘而是用来判断边缘像素是位于图像的明区还是暗区Canny此方法不容易受噪声的干扰能够检测到真正的弱边缘。在edge函数中最有效的边缘检测方法是Canny方法。该方法的优点在于使用两种不同的阈值分别检测强边缘和弱边缘并且仅当弱边缘和强边缘相连时才将弱边缘包含在输出图像中。因此这种方法不容易被噪声”填充“更容易检测出真正的弱边缘。
2 使用opencv的Canny算子实现边缘检测
Canny边缘检测是一种经典的边缘检测算法由John F.在1986年提出。它被广泛应用于计算机视觉和图像处理领域是一种多阶段的边缘检测算法能够有效地检测图像中的边缘并抑制噪声得到清晰准确的边缘信息并且对噪声具有一定的鲁棒性。Canny被认为是最优的边缘检测算法。
2.1 检测原理
2.1.1 应用高斯滤波去除图像噪声
由于图像边缘非常容易受到噪声的干扰因此为了避免检测到错误的边缘信息通常需要对图像进行滤波以去除噪声。滤波的目的是平滑一些纹理较弱的非边缘区域以便得到更准确的边缘。在实际处理过程中通常采用高斯滤波去除图像中的噪声。
图 10-1 演示了使用高斯滤波器 T 对原始图像 O 中像素值为 226 的像素点进行滤波得到该点在滤波结果图像 D 内的值的过程。 在滤波过程中我们通过滤波器对像素点周围的像素计算加权平均值获取最终滤波结果。对于高斯滤波器 T越临近中心的点权值越大。在图 10-1 中对图像 O 中像素值为 226 的像素点使用滤波器 T 进行滤波的计算过程及结果为
结果 156×(197×125×1106×2156×1159×1149×140×3107×45×371×1163×2198×4226×8223×4156×2222×137×368×4193×3157×142×172×1250×241×175×1) 138 当然高斯滤波器高斯核并不是固定的例如它还可以是 滤波器的大小也是可变的高斯核的大小对于边缘检测的效果具有很重要的作用。滤波器的核越大边缘信息对于噪声的敏感度就越低。不过核越大边缘检测的定位错误也会随之增加。通常来说一个 5×5 的核能够满足大多数的情况。
2.1.2 计算梯度幅值和方向
前面一节讲了如何计算图像梯度的幅度。在这里我们关注梯度的方向梯度的方向与边缘的方向是垂直的。
边缘检测算子返回水平方向的Gx和垂直方向的Gy。梯度的幅度和方向用角度值表示为 如果某个像素点是边缘则其梯度方向总是垂直与边缘垂直。梯度方向被归为四类垂直水平和两个对角线方向。
因此在计算梯度时我们会得到梯度的幅度和角度代表梯度的方向两个值。
图 10-2 展示了梯度的表示法。其中每一个梯度包含幅度和角度两个不同的值。为了方便观察这里使用了可视化表示方法。例如左上角顶点的值“2↑”实际上表示的是一个二元数对“(2, 90)”表示梯度的幅度为 2角度为 90°。
2.1.3 非极大值抑制
在获得了梯度的幅度和方向后遍历图像中的像素点去除所有非边缘的点。
在具体实现时逐一遍历像素点判断当前像素点是否是周围像素点中具有相同梯度方向的最大值并根据判断结果决定是否抑制该点。
通过以上描述可知该步骤是边缘细化的过程。针对每一个像素点 如果该点是正/负梯度方向上的局部最大值则保留该点。 如果不是则抑制该点归零。
在图 10-3 中A、B、C 三点具有相同的方向梯度方向垂直于边缘。判断这三个点是否为各自的局部最大值如果是则保留该点否则抑制该点归零。
经过比较判断可知A 点具有最大的局部值所以保留 A 点称为边缘其余两点B和 C被抑制归零。
在图 10-4 中黑色背景的点都是向上方向梯度水平边缘的局部最大值。因此这些点会被保留其余点被抑制处理为 0。这意味着这些黑色背景的点最终会被处理为边缘点而其他点都被处理为非边缘点。
“正/负梯度方向上”是指相反方向的梯度方向。例如在图 10-5 中黑色背景的像素点都是垂直方向梯度向上、向下方向上即水平边缘的局部最大值。这些点最终会被处理为边缘点。 经过上述处理后对于同一个方向的若干个边缘点基本上仅保留了一个因此实现了边缘细化的目的。
2.1.4 应用双阈值确定边缘
完成上述步骤后图像内的强边缘已经在当前获取的边缘图像内。但是一些虚边缘可能也在边缘图像内。这些虚边缘可能是真实图像产生的也可能是由于噪声所产生的。对于后者必须将其剔除。
设置两个阈值其中一个为高阈值 maxVal另一个为低阈值 minVal。根据当前边缘像素的梯度值指的是梯度幅度下同与这两个阈值之间的关系判断边缘的属性。具体步骤为 如果当前边缘像素的梯度值大于或等于 maxVal则将当前边缘像素标记为强边缘。 如果当前边缘像素的梯度值介于 maxVal 与 minVal 之间则将当前边缘像素标记为虚边缘需要保留。 如果当前边缘像素的梯度值小于或等于 minVal则抑制当前边缘像素。
在上述过程中我们得到了虚边缘需要对其做进一步处理。一般通过判断虚边缘与强边 缘是否连接来确定虚边缘到底属于哪种情况。通常情况下如果一个虚边缘 与强边缘连接则将该边缘处理为边缘。 与强边缘无连接则该边缘为弱边缘将其抑制。
在图 10-6 中左图显示的是三个边缘信息右图是对边缘信息进行分类的示意图具体划分如下 A 点的梯度值值大于 maxVal因此 A 是强边缘。 B 和 C 点的梯度值介于 maxVal 和 minVal 之间因此 B、C 是虚边缘。 D 点的梯度值小于 minVal因此 D 被抑制抛弃。 图 10-7 显示了对图 10-6 中的虚边缘 B 和 C 的处理结果。其中 B 点的梯度值介于 maxVal 和 minVal 之间是虚边缘但该点与强边缘不相连故将其抛弃。 C 点的梯度值介于 maxVal 和 minVal 之间是虚边缘但该点与强边缘 A 相连故将其保留。 注意高阈值 maxVal 和低阈值 minVal 不是固定的需要针对不同的图像进行定义。
10-8 给出了一个 Canny 边缘检测的效果图。 2.1.5 检测边缘连接
Canny边缘检测的边缘连接是指将非极大值抑制得到的边缘点连接成连续的边缘线。这个过程是Canny边缘检测算法中的最后一步它的目的是去除由于非极大值抑制产生的间断的边缘点从而得到更加准确和连续的边缘检测结果。
Canny边缘检测的边缘连接过程 遍历梯度幅值图像首先遍历经过非极大值抑制后的梯度幅值图像即只有边缘上的像素点的梯度幅值被保留其他像素点的梯度幅值为零。 标记边缘点对于每个强边缘像素点梯度幅值大于高阈值将其标记为边缘点。强边缘像素点是图像中梯度值较大的像素点它们代表了图像中明显的边缘。 边缘连接对于每个边缘点的相邻像素点如果其梯度幅值大于低阈值并且没有被标记为边缘点则将其标记为弱边缘点并递归地进行边缘连接。这一步的目的是将与强边缘像素点相邻的弱边缘像素点连接到边缘线上。 递归连接在边缘连接过程中如果某个弱边缘像素点被标记为边缘点则会继续检查该像素点的相邻像素点以便将所有与强边缘像素点相邻的弱边缘像素点连接到边缘线上。 结束条件边缘连接的递归过程会一直进行直到所有与强边缘像素点相邻的弱边缘像素点都被标记为边缘点没有更多的像素点可以连接。 非边缘点处理经过边缘连接后所有未被标记为边缘点的像素点视为非边缘点并抑制其梯度幅值为零从而得到最终的边缘图像。
2.2 Canny函数原型
edges cv.Canny( image, threshold1, threshold2[, apertureSize[, L2gradient]])
参数 edges 为计算得到的边缘图像。 image 为 8 位输入图像。 threshold1 表示处理过程中的第一个阈值。 threshold2 表示处理过程中的第二个阈值。 apertureSize 表示 Sobel 算子的孔径大小。 L2gradient 为计算图像梯度幅度gradient magnitude的标识。其默认值为 False。如果为 True则使用更精确的 L2 范数进行计算即两个方向的导数的平方和再开方否则使用 L1 范数直接将两个方向导数的绝对值相加。 2.3 检测代码
import cv2 as cv
import numpy as np
from matplotlib import pyplot as plt# 1 读取图像
img cv.imread(../data/cat_dog.jpg)# 2 将图片由BGR转为RGB
img cv.cvtColor(img, cv.COLOR_BGR2RGB)# 2 边缘检测canny算子阈值128-200低于128的像素点认为是边缘
# 高于200的像素点认为是边缘中间值的像素点如果与边缘点相连则认为是边缘
edge cv.Canny(img, 128, 200)# 3 图像显示
plt.figure(figsize(10, 8), dpi100)
plt.subplot(121), plt.imshow(img), plt.title(original)
plt.xticks([]), plt.yticks([])
plt.subplot(122), plt.imshow(edge, cmapplt.cm.gray), plt.title(Canny)
plt.xticks([]), plt.yticks([])
plt.show()
运行代码显示
