第七章图像分割（2）.pptVIP

第七章 图像分割（2） 7.4 图像的边缘检测 边缘检测是基于灰度不连续性进行的分割方法。在第4章中，我们讨论了用差分、梯度、拉普拉斯算子及各种高通滤波处理方法对图像边缘进行增强 ，只要再进行一次门限化的处理，便可以将边缘增强的方法用于边缘检测。 7.4.1 梯度算子 第四章介绍了梯度对应于一阶导数，相应的梯度算子就对应于一阶导数算子。对于一个连续函数f(x,y) ，其在(x,y)处的梯度定义如下： 在实际应用中常常采用小型模板，然后利用卷积运算来近似，Gx和Gy各自使用一个模板。常用的有一下这些： 1. Roberts算子 2. Prewitt算子 3. Sobel算子 通过算子检测后，还需作二值处理从而找到边界点。 后面分别给出了利用这三个算子进行边缘检测的不同效果。 完成上述检测的实现程序如下： I = imread(blood1.tif); imshow(I); BW1 = edge(I,roberts); %进行Roberts算子边缘检测，门限值采用matlab默认值 BW2 = edge(I,prewitt); %进行Prewitt算子边缘检测，门限值采用matlab默认值 BW3 = edge(I,sobel); %进行Sobel算子边缘检测，门限值采用matlab默认值 figure,imshow(BW1,[]); figure,imshow(BW2,[]); figure,imshow(BW3,[]); 7.4.2 拉普拉斯算子 拉普拉斯（Laplacian）是一种二阶导数算子，也是借助模板来实现的。这里对模板有一些基本要求：模板中心的系数为正，其余相邻系数为负，且所有的系数之和为零。 % H = FSPECIAL(log,HSIZE,SIGMA) returns a rotationally symmetric Laplacian of Gaussian filter of size HSIZE with standard deviation SIGMA (positive). HSIZE can be a vector specifying the number of rows and columns in H or a scalar, in which case H is a square matrix. %The default HSIZE is [5 5], the default SIGMA is 0.5. H = fspecial(log) H = 0.0448 0.0468 0.0564 0.0468 0.0448 0.0468 0.3167 0.7146 0.3167 0.0468 0.0564 0.7146 -4.9048 0.7146 0.0564 0.0468 0.3167 0.7146 0.3167 0.0468 0.0448 0.0468 0.0564 0.0468 0.0448 实现的MATLAB程序如下： I = imread(blood1.tif); imshow(I); BW4 = edge(I,log); figure,imshow(BW4,[]); 7.4.3 方向算子 方向算子是利用一组模板对图像中的同一像素求卷积，然后选取其中最大的值作为边缘强度，而将与之对应的方向作为边缘方向。它相对于梯度算子的优点是不仅仅只考虑水平和垂直方向，还可以检测其他方向上的边缘，但计算量将大大增加。 常用的有8方向Kirsch（3×3）模板，如图7.7所示，方向间的夹角为45o。 7.4.4 Canny边缘检测算子 Canny的主要工作是推导了最优边缘检测算子。他考核边缘检测算子的指标是：低误判率，即尽可能少地把边缘点误认为是非边缘点；高定位精度，即准确地把边缘点定位在灰度变化最大的像素上；抑制虚假边缘。 判断一个像素是否为边缘点的条件为： （1）像素(i, j)的边缘强度大于沿梯度方向的两个相邻像素的边缘强度； （2）与该像素梯度方向上相邻两点的方向差小于45o； （3）以该像素为中心的3×3邻域中的边缘强度的极大值小于某个阈值 。 I=imread(rice.png); imshow(I); BW1=edge(I,log); BW2=edge(I,canny); figure,imshow(BW1); figure,imshow(BW2); 结论 Canny算子的检测比较优越，它可以减少小模板检测中边缘中断，有利于得到较完整的边缘。 另