青岛大学专业课程设计院系: 自动化学院专业: 电子信息工程班级: 08级电子信息工程3班学生姓名: 刘法指导教师: 王汉萍庄晓东日期: 2011年12月23日题目：图像边缘检测方法的研究与实现一、边缘检测以及相关概念1．1边缘，边缘检测的介绍边缘(edge)是指图像局部强度变化最显著的部分．边缘主要存在于目标与目标、目标与背景、区域与区域(包括不同色彩)之间，是图像分割、纹理特征和形状特征等图像分析的重要基础．图像分析和理解的第一步常常是边缘检测(edge detection)．边缘检测是指使用数学方法提取图像像元中具有亮度值（灰度）空间方向梯度大的边、线特征的过程。
在讨论边缘算子之前，首先给出一些术语的定义：边缘点：图像中具有坐标]i且处在强度显著变化的位置上的点．,[j边缘段：对应于边缘点坐标]i及其方位 ，边缘的方位可能是梯度角．,[j边缘检测器：从图像中抽取边缘(边缘点和边缘段)集合的算法．轮廓：边缘列表，或是一条表示边缘列表的拟合曲线．边缘连接：从无序边缘表形成有序边缘表的过程．习惯上边缘的表示采用顺时针方向序．边缘跟踪：一个用来确定轮廊的图像（指滤波后的图像）搜索过程．边缘点的坐标可以是边缘位置像素点的行、列整数标号，也可以在子像素分辨率水平上表示．边缘坐标可以在原始图像坐标系上表示，但大多数情况下是在边缘检测滤波器的输出图像的坐标系上表示，因为滤波过程可能导致图像坐标平移或缩放．边缘段可以用像素点尺寸大小的小线段定义，或用具有方位属性的一个点定义．请注意，在实际中，边缘点和边缘段都被称为边缘．边缘连接和边缘跟踪之间的区别在于：边缘连接是把边缘检测器产生的无序边缘集作为输入，输出一个有序边缘集；边缘跟踪则是将一幅图像作为输入，输出一个有序边缘集．另外，边缘检测使用局部信息来决定边缘，而边缘跟踪使用整个图像信息来决定一个像素点是不是边缘．1．2 边缘检测算子边缘检测是图像特征提取的重要技术之一, 边缘常常意味着一个区域的终结和另一个区域的开始. 图像的边缘包含了物体形状的重要信息,它不仅在分析图像时大幅度地减少了要处理的信息量,而且还保护了目标的边界结构. 因此,边缘检测可以看做是处理许多复杂问题的关键.边缘检测的实质是采用某种算法来提取出图像中对对象与背景间的交界线。
图像灰度的变化情况可以用图像灰度分布的梯度来反映，因此可以用局部图像微分技术来获取边缘检测算子。
经典的边缘检测方法是对原始图像中的像素的某个邻域来构造边缘检测算子。
以下是对几种经典的边缘检测算子进行理论分析，并对各自的性能特点做出比较和评价。
边缘检测的原理是：由于微分算子具有突出灰度变化的作用，对图像进行微分运算，在图像边缘处其灰度变化较大，故该处微分计算值教高，可将这些微分值作为相应点的边缘强度，通过阈值判别来提取边缘点，即如果微分值大于阈值，则为边缘点。
Roberts,Sobel,Prewwit是基于一阶导数的边缘检测算子，图像的边缘检测是通过2*2或者3*3模板作为核与该图像中的每个像素点做卷积和运算，然后选取合适的阈值以提取边缘。
Laplace边缘检测算子是基于二阶导数的边缘检测算子，该算子对噪声敏感。
Laplace算子的改进方式是先对图像进行平滑处理，然后再应用二阶导数的边缘检测算子，其代表是拉普拉斯高斯（LOG）算子。
前边介绍的边缘检测算法是基于微分方法的，其依据是图像的边缘对应一阶导数的极大值点和二阶导数过零点。
Canny算子是另外一类边缘检测算子，它不是通过微分算子检测边缘，而是在满足一定约束条件下推导出的边缘检测最优化算子。
1．3 边缘检测算法对于边缘的检测常常借助于空域微分算子进行,通过将其模板与图像卷积完成。
两个具有不同灰度值的相邻区域之间总存在灰度边缘。
灰度边缘是灰度值不连续(或突变) 的结果,这种不连续常可利用求一阶和二阶导数方便地检测到。
已有的局部技术边缘检测方法,主要有一次微分(Sobel 算子、Robert s 算子等) 、二次微分(拉普拉斯算子等)。
这些边缘检测器对边缘灰度值过渡比较尖锐且噪声较小等不太复杂的图像,大多数提取算法均可以取得较好的效果。
但对于边缘复杂、采光不均匀的图像来说,则效果不太理想。
主要表现为边缘模糊、边缘非单像素宽、弱边缘丢失和整体边缘的不连续等方面。
用算子检测图像边缘的方法是用小区域模板对图像进行处理,即采用卷积核作为掩模模板在图像中依次移动,完成图像中每个像素点同模板的卷积运算,最终输出的边缘幅度结果可以检测出图像的边缘。
卷积运算是一种邻域运算。
图像处理认为:某一点像素的结果不但和本像素灰度有关,而且和其邻域点值有关。
运用模板在图像上依此对每一个像素进行卷积, 即模板上每一个点的值与其在图像上当前位置对应的像素点值相乘后再相加,得出的值就是该点处理后的新值。
边缘检测算法有如下四个步骤：滤波：边缘检测算法主要是基于图像强度的一阶和二阶导数，但导数的计算对噪声很敏感，因此必须使用滤波器来改善与噪声有关的边缘检测器的性能。
需要指出，大多数滤波器在降低噪声的同时也导致了边缘强度的损失，因此，增强边缘和降低噪声之间需要折衷。
增强：增强边缘的基础是确定图像各点邻域强度的变化值。
增强算法可以将邻域（或局部）强度值有显著变化的点突显出来。
边缘增强一般是通过计算梯度幅值来完成的。
检测：在图像中有许多点的梯度幅值比较大，而这些点在特定的应用领域中并不都是边缘，所以应该用某种方法来确定哪些点是边缘点。
最简单的边缘检测判据是梯度幅值阈值判据。
定位：如果某一应用场合要求确定边缘位置，则边缘的位置可在子像素分辨率上来估计，边缘的方位也可以被估计出来。
二、经典边缘检测算子的理论分析和比较2．1 Roberts 算子Roberts 算子是一种斜向偏差分的梯度计算方法, 梯度的大小代表边缘的强度, 梯度的方向与边缘走向垂直。
该算子定位精度高, 但容易丢失部分边缘。
因为没进行平滑处理, 不具有抑制噪声的能力。
用该算子处理边缘陡峭度高且噪声小的图像效果较佳。
Roberts 交叉算子为梯度幅值计算提供了一种简单的近似方法： ]1,[],1[]1,1[],[],[+-++++-=j i f j i f j i f j i f j i G (1．1)用卷积模板，上式变成：y x G G j i G +=],[ (1．2)其中x G 和y G 由下面的模板计算：G x =-1001G y =-0110(1．3) 同前面的22⨯梯度算子一样，差分值将在内插点[/,/]i j ++1212处计算．Roberts 算子是该点连续梯度的近似值，而不是所预期的点[,]i j 处的近似值。
2．2 Sobel 算子、Prewitt 算子Sobel 算子和P rew it t 算子都是一阶的微分算子,都是先对图像进行平滑处理, 虽然两者都是加权平均滤波, 但是前者邻域的像素对当前像素产生的影响不是等价的, 距离不同的像素具有不同的权值, 对算子结果产生的影响也不同。
这两种算子对噪声都有一定的抑制作用, 但不能完全排除检测结果中出现虚假边缘的情况。
这两者对灰度渐变低噪声的图像有较好的检测效果, 但是对于混合多复杂噪声的图像处理效果就不理想了。
正如前面所讲，采用33⨯邻域可以避免在像素之间内插点上计算梯度．考虑一下上图中所示的点[,]i j 周围点的排列。
Sobel 算子也是一种梯度幅值，M s s x y=+22 (1．4) 其中的偏导数用下式计算：)()()()(456210670432a ca a a ca a s a ca a a ca a s y x ++-++=++-++= (1．5)其中常数2=c和其他的梯度算子一样，s x 和s y 可用卷积模板来实现： 图1．1请注意这一算子把重点放在接近于模板中心的像素点。
Sobel 算子是边缘检测器中最常用的算子之一。
图1．2图1．2 用于说明Sobel 算子和Prewitt 算子的邻域像素点标记Prewitt 算子与Sobel 算子的方程完全一样，只是常量c=1。
所以101101101---=x s 111000111---=y s ．6) 请注意，与Sobel 算子不同，这一算子没有把重点放在接近模板中心的像素点。
2．3拉普拉斯算子平滑过的阶跃边缘二阶导数是一个在边缘点处过零的函数。
拉普拉斯算子是二阶导数的二维等效式。
函数),(y x f 的拉普拉斯算子公式为∇=+22222f f x f y ∂∂∂∂ (1．7) 使用差分方程对x 和y 方向上的二阶偏导数近似如下：],[])1,[2]2,[(],[]1,[]),[]1,[(22j i f j i f j i f xj i f x j i f xj i f j i f x G x f x ++-+=-+=-+==∂∂∂∂∂∂∂∂∂∂ (1．8)这一近似式是以点]1,[+j i 为中心的。
用1-j 替换j ，得到=22xf ∂∂]1,[]),[2]1,[(-+-+j i f j i f j i f (1．9) 它是以点[,]i j 为中心的二阶偏导数的理想近似式，类似地，=22yf ∂∂],1[]),[2],1[(j i f j i f j i f -+-+ (1．10) 把这两个式子合并为一个算子，就成为下面能用来近似拉普拉斯算子的模板：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-≈∇010******* （1．11)当拉普拉斯算子输出出现过零点时就表明有边缘存在，其中忽略无意义的过零点(均匀零区)。
原则上，过零点的位置精度可以通过线性内插方法精确到子像素分辨率，不过由于噪声，结果可能不会很精确。
2．4 LoG 算法正如上面所提到的，利用图像强度二阶导数的零交叉点来求边缘点的算法对噪声十分敏感，所以，希望在边缘增强前滤除噪声。
为此，Marr 和Hildreth[146]将高斯滤波和拉普拉斯边缘检测结合在一起，形成LoG （Laplacian of Gaussian, LoG ）算法，也称之为拉普拉斯高斯算法。
LoG 边缘检测器的基本特征是：1． 平滑滤波器是高斯滤波器；2． 增强步骤采用二阶导数(二维拉普拉斯函数)；3． 边缘检测判据是二阶导数零交叉点并对应一阶导数的较大峰值；4． 使用线性内插方法在子像素分辨率水平上估计边缘的位置。
LoG 算子的输出h x y (,)是通过卷积运算得到的：h x y g x y f x y (,)[(,)(,)]=∇*2 (1．12) 根据卷积求导法有h x y g x y f x y (,)[(,)](,)=∇*2 (1．13) 其中：∇=+-⎛⎝ ⎫⎭⎪-+2222422222g x y x y e x y (,)σσσ (1．14)称之为墨西哥草帽算子．2．5 Canny 边缘检测器检测阶跃边缘的基本思想是在图像中找出具有局部最大梯度幅值的像素点。
