SOBEL算子

《SOBEL算子》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、Sobel算子及改进算法一、Sobel边缘检测算子在讨论边缘算子之前，首先给出一些术语的定义：(1)边缘点：图像中具有坐标[i，j]，且处在强度显著变化的位置上的点。(2)边缘段：对应于边缘点坐标[i，j]及其方位，边缘的方位可能是梯度角。(3)边缘检测器：从图像中抽取边缘(边缘点和边缘段)集合的算法。(4)边缘跟踪：一个用来确定轮廊的图像(指滤波后的图像)搜索过程。二、Sobel算子的基本原理采用3×3邻域可以避免在像素之间内插点上计算梯度。Sobel算子也是一种梯度幅值，即：其中的偏导数用下式计算：Sx=(a2+ca3+a4)一(a0+ca7+a6)Sy=(a

2、0+cal+a2)一(a6+ca5+a4)其中常数c=2。和其他的梯度算子一样，Sx和Sy可用卷积模板来实现(见图2)。121000-1-2-1-101-202-101SxSy图2卷积模板这一算子把重点放在接近于模板中心的像素点图2和图3表明了这一算子的作用。Sobel算子是边缘检中最常用的算子之一。a0ala2a7[i，j]a3a6a5a4图3用于说明Sobel算子的邻域像素点标记算法的基本原理：由于图像边缘附近的亮度变化比较大，所以可以把那些在邻域内灰度超过某个值的像素点当作边缘点。算法的主要步骤：1)分别将2个方向模板沿着图像从一个像素移动到另一个像素，并将

3、像素的中心与图像中的某个像素位置重合；2)将模板内的系数与其图像上相对应的像素值相乘；3)将所有相乘的值相加；4)将2个卷积的最大值。赋给图像中对应模板中心位置的像素，作为该像素新的灰度值；5)选取合适的阈值TH，若新像素灰度值≥TH，则判读该像素点为图像边缘点。三、Sobel算子图像边缘检测的MATLB程序实现及仿真>>f=imread('peppers.png');>>f=rgb2gray(f);>>f=im2double(f);%使用垂直Sobel算子，自动选择阈值>>[VSFATThreshold]=edge(f,'sobel','vertical');%

4、边缘检测>>figure,imshow(f),title('原始图像');%显示原始图像>>figure,imshow(VSFAT),title('垂直图像边缘检测');%显示边缘检测图像>>s45=[-2-10;-101;012];%使用指定45度角Sobel算子滤波器，指定阈值>>FST45=imfilter(f,s45,'replicate');>>SFST45=SFST45>=Threshold;>>figure,imshow(SFST45),title('45度角图像边缘检测');%显示边缘检测图像>>原始图像垂直边缘检测图像45度角边缘检测图像本实验中

5、使用Sobel算子在3个方向进行了图像边缘检测，从程序运行结果可以看出，45度角Sobel算子生成的边缘检测图像呈现出浮雕效果，水平和垂直Sobel算子检测出的边缘多于单个方向上检测出的边缘。Sobel算子利用像素的上、下、左、右邻域的灰度加权算法，根据在边缘点处达到极值这一原理进行边缘检测。该方法不但产生较好的检测效果，而且对噪声具有平滑作用，可以提供较为精确的边缘方向信息。但是，在抗噪声好的同时也存在检测到伪边缘，定位精度不高的缺点。四、Sobel的改进算法1、方向模板的改变因为图像目标的边缘是灰度不连续性的反映，边缘的种类可以粗略地区分为以下2类：1)阶跃性

6、边缘，它两边的像素点灰度值显著不同。2)屋顶状边缘，它位于灰度值从增加到减少的变化转折点。为了能够更准确地描述出图像缘点，减少噪声对检测结果的影响，提高算子的抗噪能力，重新构造了4个5×5大小的模板，模板中各个位置的权重是由该位置到中心点的距离以及该位置在模板中所在的方位决定的，等距离的点，具有相同的权重，如图4和图5所示。在本文的算法中，我们选用有最高输出模板所对应边缘梯度值来作为像元的边缘梯度强度。2、边缘的细化处理图像的边缘中包含了图像的大量信息。但是由上述模板计算出来的值仅仅是图像对应的梯度图像，边缘出现在梯度值较大的地方，若想提取图像的边缘必须设定阈值将

7、其二值化。图4边缘检测的4个方向x水平方向y垂直方向23-3-234-4-366-6-634-4-223-3-22363234643-3-4-6-4-3-2-3-6-3-2-6-2-3-2-2-6-42-3-442-2462232-245°方向135°方向-2-3-2-62-4-6-234-4-3264-26232图5由一阶梯度算子得到图像的梯度图像。一般来说，其边缘较粗。若直接对梯度图像设定阈值进行二值化，很难找到合适的阈值，使得检测出来的边缘达到要求，这不便于边缘连接与边缘特征提取等后期处理。因此在对图像梯度图像进行二值化前，有必要对所检测出来的梯度边缘进行细

8、化处理。在