mv-chp06边缘检测

*,北京邮电大学自动化,*,边缘检测,Edge Detection,第,6,章,2024/10/1,1,北京邮电大学自动化,主讲课本：,1,、,贾云得，机器视觉，科学出版社，2000年,2024/10/1,2,北京邮电大学自动化,深度阅读：,1,、,D.A.Forsyth,J.Ponce,Computer Vision-A modern approach,，清华出版社，,2004,2,、,Milan,Sonka,等著，艾海舟等译，,图像处理、分析与机器视觉,，人民邮电出版社，,2003,年,2024/10/1,3,北京邮电大学自动化,边缘(,edge),是指图像局部强度变化最显著的部分边缘主要存在于目标与目标、目标与背景、区域与区域(包括不同色彩)之间，,图像分割、纹理特征和形状特征等图像分析的重要基础,图像强度的不连续可分为：,(1)阶跃不连续，即图像强度在不连续处的两边的像素灰度值有着显著的差异；,(2)线条不连续，即图像强度突然从一个值变化到另一个值，保持一较小行程后又回到原来的值,1、边缘检测的基本定义,2024/10/1,4,北京邮电大学自动化,两种常见的边缘一阶导数和二阶导数示意图,(,a),阶跃函数 (,b),线条函数,理论曲线,实际曲线,一阶导数,二阶导数,2024/10/1,5,北京邮电大学自动化,边缘点,(Edge point),：在亮度显著变化的位置上的点,边缘段,(Edge segment),：对应于边缘点坐标及其方位,边缘检测器,(Edge detector),：从图像中抽取边缘集合的算法,轮廓,(Boundary),：边缘列表或一条表示边缘列表的拟合曲线,边缘连接,(Edge linking),：从无序边缘表形成有序边缘表的过程,边缘跟踪,(Edge tracking),：一个用来确定轮廊的图像搜索过程,2、术语定义,2024/10/1,6,北京邮电大学自动化,3、梯度,梯度是一阶导数的二维等效式，定义为矢量,(1)向量的方向就是函数增大时的最大变化率方向；,(2)梯度的幅值和方向:,2024/10/1,7,北京邮电大学自动化,用差分来近似梯度：,用简单卷积模板表示:,注意：与课本不同！,2024/10/1,8,北京邮电大学自动化,采用上面公式计算的梯度近似值,Gx,和,Gy,并不位于同一位置，,Gx,实际上是内插点,i,+1/2,j,处的梯度近似值，,Gy,是内插点,i,j,+1/2,处的梯度近似值由于这个缘故，人们常常使用 一阶差分模板（而不用 或 模板）来求,x,和,y,的偏导数。,2024/10/1,9,北京邮电大学自动化,基本步骤：,（,1,）滤波：,边缘检测算法主要是基于图像强度的一阶和二阶导数，但导数的计算对噪声很敏感，因此必须使用滤波器来改善与噪声有关的边缘检测器的性能需要指出，大多数滤波器在降低噪声的同时也导致了边缘强度的损失，因此，增强边缘和降低噪声之间需要折衷,4、边缘检测算法,2024/10/1,10,北京邮电大学自动化,（,2,）增强：增强边缘的基础是确定图像各点邻域强度的变化值增强算法可以将邻域（或局部）强度值有显著变化的点突显出来边缘增强一般是通过计算梯度幅值来完成的,2024/10/1,11,北京邮电大学自动化,（,3,）检测：在图像中有许多点的梯度幅值比较大，而这些点在特定的应用领域中并不都是边缘，所以应该用某种方法来确定哪些点是边缘点最简单的边缘检测判据是梯度幅值阈值判据,2024/10/1,12,北京邮电大学自动化,常用边缘检测器,(1),Roberts,算子,用卷积模板表示：,梯度交叉算子,梯度幅值计算近似方法：,2024/10/1,13,北京邮电大学自动化,(2)Sobel,算子,梯度幅值:,其中的偏导数用下式计算：,c=2,用卷积模板来实现,请注意：,Sobel,算子把重点放在接近于模板中心的像素点,Sobel,算子是边缘检测器中最常用的算子之一,2024/10/1,14,北京邮电大学自动化,(3)Prewitt,算子,与,Sobel,算子的方程完全一样，但,c=1,该算子没有把重点放在接近模板中心的 像素点,2024/10/1,15,北京邮电大学自动化,(4)各种算法的比较,图,6.2,用于边缘检测的测试图像,(a),原始图像,(b)7x7,高斯滤波的图像,2024/10/1,16,北京邮电大学自动化,图,6.3,各种边缘检测器对未经滤波的图像进行边缘检测的比较,(a)Roberts,交叉算子,(b),Sobel,算子,(c)Prewitt,算子,2024/10/1,17,北京邮电大学自动化,图,6.4,各种边缘检测器对滤波后的图像进行边缘检测的结果,Roberts,交叉算子,Sobel,算子,Prewitt,算子,2024/10/1,18,北京邮电大学自动化,5、,二阶微分算子,如果所求的一阶导数高于某一阈值，则确定该点为边缘点这样做会导致检测的边缘点太多一种更好的方法就是求梯度局部最大值对应的点，并认定它们是边缘点，通过去除一阶导数中的非局部最大值，可以检测出更精确的边缘一阶导数的局部最大值对应着二阶导数的零交叉点这样，通过找图像强度的二阶导数的零交叉点就能找到边缘点,阈值,边缘太宽,零交叉点,2024/10/1,19,北京邮电大学自动化,拉普拉斯算子,拉普拉斯算子是二阶导数的二维等效式:,这一近似式是以点,i+1,j,为中心的用,i-1,替换,得到以,i,j,为中心的,:,2024/10/1,20,北京邮电大学自动化,用算子表示：,希望邻域中心点具有更大的权值,2024/10/1,21,北京邮电大学自动化,拉普拉斯的二阶方向导数算子在机器视觉中并不常用，因为任何包含有二阶导数的算子比只包含有一阶导数的算子更易受噪声的影响甚至一阶导数很小的局部峰值也能导致二阶导数过零点为了避免噪声的影响，必须采用特别有效的滤波方法下面我们讨论高斯滤波与二阶导数相结合的边缘检测方法,拉普拉斯算子的缺点,2024/10/1,22,北京邮电大学自动化,6、,LoG,算法,基本特征：,平滑滤波器是高斯滤波器,增强步骤采用二阶导数(二维拉普拉斯函数,),边缘检测判据是二阶导数零交叉点并对应一阶导数的较大峰值,使用线性内插方法在子像素分辨率水平上估计边缘的位置,Marr,和,Hildreth,将高斯滤波和拉普拉斯边缘检测结合在一起，形成,LoG（Laplacian,of Gaussian）,算法，也称之为拉普拉斯高斯算法,2024/10/1,23,北京邮电大学自动化,LoG,算子的输出是通过卷积运算得到的,根据卷积求导法有,:,其中：,两种方法在数学上是等价的：,1,求图像与高斯滤波器卷积，再求卷积的拉普拉斯变换,2,求高斯滤波器的拉普拉斯变换，再求与图像的卷积,该函数称为,:,LoG,2024/10/1,24,北京邮电大学自动化,一维和二维高斯函数的拉普拉斯变换图的翻转图，其中 =2,称之为墨西哥草帽算子,LoG,函数,2024/10/1,25,北京邮电大学自动化,5,X5,拉普拉斯高斯模板,2024/10/1,26,北京邮电大学自动化,拉普拉斯高斯边缘检测结果,2024/10/1,27,北京邮电大学自动化,7、,Canny,边缘检测器,检测阶跃边缘的基本思想是在图像中找出具有局部最大梯度幅值的像素点检测阶跃边缘的大部分工作集中在寻找能够用于实际图像的梯度数字逼近图像梯度逼近必须满足两个要求：,(1),逼近必须能够抑制噪声效应，,(2),必须尽量精确地确定边缘的位置抑制噪声和边缘精确定位是无法同时得到满足的，也就是说，边缘检测算法通过图像平滑算子去除了噪声，但却增加了边缘定位的不确定性；反过来，若提高边缘检测算子对边缘的敏感性，同时也提高了对噪声的敏感性有一种线性算子可以在抗噪声干扰和精确定位之间提供最佳折衷方案，它就是高斯函数的一阶导数。,2024/10/1,28,北京邮电大学自动化,2）使用,2x2,一阶有限差分计算局部梯度的两个阵列,P,与,Q：,1）求图像与高斯平滑滤波器卷积:,实际上可以采用前面所说的,Roberts,、,Prewitt,和,Sobel,算子！,2024/10/1,29,北京邮电大学自动化,3）幅值和方位角:,注意：与课本不同！,2024/10/1,30,北京邮电大学自动化,4）非极大值抑制（,Non-Maxima Suppression,，,NMS,）,：,幅值图像阵列,M,i,j,的值越大，其对应的图像梯度值也越大，但这还不足以确定边缘，因为这里仅仅把图像快速变化的问题转化成求幅值阵列,M,i,j,的局部最大值问题为确定边缘，必须细化幅值图像中的屋脊带（,Ridge,），即只保留幅值局部变化最大的点，即追踪所有脊的顶部，并将所有不在脊的顶部的像素设为零，以便在输出中给出一条细线这一过程叫非极大值抑制（,Non-Maxima Suppression,，,NMS,），它会生成细化的边缘,2024/10/1,31,北京邮电大学自动化,*方向角：,非极大值抑制通过抑制梯度线上所有非屋脊峰值的幅值来细化,M,i,j,中的梯度幅值屋脊这一算法首先将梯度角,i,j,的变化范围减小到圆周的四个扇区之一，如图所示。四个扇区的标号为,0,到,3,，对应着,3x3,邻域内元素的四种可能组合，任何通过邻域中心的点必通过其中一个扇区,2024/10/1,32,北京邮电大学自动化,*幅值：,算法使用一个,3x3,邻域作用于幅值阵列,M,i,j,的所有点在每一点上，邻域的中心像素,M,i,j,与沿着梯度线的两个元素进行比较，其中梯度线是由邻域的中心点处的扇区值,i,j,给出的如果在邻域中心点处的幅值,M,i,j,不比沿梯度线方向上的两个相邻点幅值大，则,M,i,j,赋值为零这一过程可以把,M,i,j,宽屋脊带细化成只有一个像素点宽在非极大值抑制过程中，保留了屋脊的高度值,2024/10/1,33,北京邮电大学自动化,*幅值：,表示非极大值抑制过程,N,i,j,中的非零值对应着图像强度阶跃变化处的对比度尽管在边缘检测的第一步对图像进行了平滑，但非极大值抑制幅值图像,N,i,j,仍会包含许多由噪声和细纹理引起的假边缘段,2024/10/1,34,北京邮电大学自动化,5）阈值化,*将低于阈值的所有值赋零值，得到图像的边缘阵列,*阈值,太低,假边缘（假阳性）,.,*阈值,太高,部分轮廊丢失（假阴性）.,*选用两个阈值:更有效的阈值方案,基本思想：,取高低两个阈值,t1,t2,作用在幅值图,Ni,j，t1=2t2;,得到两个边缘图，高阈值和低阈值边缘图,;,连接高阈值边缘图，出现断点时，在低阈值边缘图,中的8邻点域搜寻边缘点。,2024/10/1,35,北京邮电大学自动化,算法6.1,Canny,边缘检测,(1),用高斯滤波器平滑图像；,(2),用一阶偏导有限差分计算梯度幅值和方向；,(3),对梯度幅值应用非极大值抑制；,(4),用双阈值算法检测和连接边缘,2024/10/1,36,北京邮电大学自动化,7,X7,高斯滤波模板,13X13,高斯滤波模板,2024/10/1,37,北京邮电大学自动化,8,、基于函数逼近的边缘检测,获得图像对应连续函数，基于函数进行检测,图,6.11,连续图像强度函数的图形表示,2024/10/1,38,北京邮电大学自动化,小面模型（,facet model,）,现在的任务是从数字图像的灰度值重构连续函数对复杂的图像来说，连续强度函数可能包含,x,和,y,的超阶幂方，从而使得重构原始函数十分困难，因此，我们将采用简单的分段函数来建立图像的模型这样，任务就变成了重构每一个分段解析函数，或小面（,Facets,）。换言之，就是试图找到在每一像素点的邻域内最能逼近该邻域强度值的简单函数，这种逼近叫小面模型（,Facet model,）。每一个像素对应一个局部函数。,2024/10/1,39,北京邮电大学自动化,图像局部函数例（三次多项式）,函数逼近：,用最小二乘法等方法计算函数参数,2024/10/1,40,北京邮电大学自动化,小面模型边缘检测器