（OpenCV+Python）六种算子 Sobel Scharr 拉普拉斯 Laplacian Roberts Prewitt.docx

《（OpenCV+Python）六种算子 Sobel Scharr 拉普拉斯 Laplacian Roberts Prewitt.docx》由会员分享，可在线阅读，更多相关《（OpenCV+Python）六种算子 Sobel Scharr 拉普拉斯 Laplacian Roberts Prewitt.docx（16页珍藏版）》请在第一文库网上搜索。

1、我们在教程顶帽与黑帽操作中留了一个小彩蛋一一形态学的梯度问题， 通常情况下，它被用于提取图像的轮廓。今天我们来了解图像边缘的另一种方法, 它将比形态学梯度更有效，适用范围也更广。Sobel算子前面的例子，已经接触到了图像卷积运算。最重要的卷积运算之一是用于计 算图像的导数（或近似导数）。为什么图像中导数的计算很重要，看下面边缘检 测的例子：很容易观察到上面图像中像素灰度值变化没有规律。一种比较好的描述这种 变化的方法是采用导数。其中梯度剧烈变化的地方代表图像灰度值变化强烈的地 方，也就是边缘。为了更好的说明，以1维图像（也就是图像的1行）为例。边缘出现在灰度 值跳变的地方，如下图所示：如果对上

2、面的1维图像求导数，得到下图，可以很明显的看到边缘所在的位 置。从上面的解释，我们可以设置一个阈值，根据局部像素变化强烈程度获取图 像边缘。Sobel算子是一个离散微分算子，计算得到的是图像梯度的近似值。SobeI算 子结合了高斯平滑和微分。假设输入图像是L ,核大小为3,通过下面运算分别计算水平方向和垂直方 向的微分：a.水平方向：b.垂直方向：-1一 2 -1Gy =000* I.+1+2 +1.具体运算为：Gx = (-1)f(-1,y-1) 0f(xty-1) 开(x+1,yT)+(-2)依-7M + O(x,y)+2f(x+1,y)(-1)40rzy*T; + (x,y*1) + f

3、(*1,y1)=EXry-I)+2 枢沙依力-*旗XlyHf(Xly+1)Gy =f(x-1,y-1) 2f(fy1)+ 开(X+1,yT)+Qf(x- 1,y) Cf(,y) + 0*f(x+1,y)+(-i)*f(x-,y+i) (-2)*f(,yi) * (-*f(.y+=If(X-Ly) + 2f(,y-i) + f(+.y-i)J-(f(-i,y+i) 2*f(x,y+i)+f(xi,y+1)结合上面结果可以计算出图像中一个点的近似梯度：G = Gl + Gl或者表示为：G=Gl + IGyl需要注意的是，当核的大小为3时，也就是上面所示的SObel核可能会产生 明显的误差（毕竟So

4、bel只是微分的近似值）。在C）PenCV-PythOn中，使用Sobel的算子的函数原型如下：dst = cv2.Sobel（src, ddepth, dx, dy, dst, ksize, scale, delta, borderTypeJ）（1）前四个是必须的参数：第一个参数是需要处理的图像。第二个参数是图像的深度，-1表示采用的是与原图像相同的深度。目标图像 的深度必须大于等于原图像的深度。dx和dy表示的是求导的阶数，0表示这个方向上没有求导，一般为0、1、 2o（2）其后是可选的参数：dst不用解释了。ksize是SObel算子的大小,必须为1、3、5、7。scale是缩放导数的比

5、例常数，默认情况下没有伸缩系数。delta是一个可选的增量，将会加到最终的dst中，同样，默认情况下没有额 外的值加到dst中。bordefype是判断图像边界的模式。这个参数默认值为cv2.BORDER_DEFAULTo现在我们来看实验图像：实战代码：import cv2import numpy as npimg = cv2.imread(,pie.png,0)x=cv2.Sobel(img,cv2.CV_l 6S, 1,0)cv2.imshow(x,x)cv2.waitKey(0)cv2.destroyAI lWindows()我们发现没有输出，原因在于图像格式的问题，在SObel函数的第

6、二个参数 这里使用了 cv2.CV_16So因为OPenCV文档中对Sobel算子的介绍中有这么一 句：“in the case of 8-bit input images it will result in truncated derivatives,o BR Sobel 函数求完导数后会有负值，还有会大于255的值。而原图像是uint8,即8位无 符号数，所以SObel建立的图像位数不够，会有截断。因此要使用16位有符号 的数据类型，即cv2.CV6S。在经过处理后，别忘了用ConVertSCaIeAbSO函数将其转回原来的Uint8形式。 否则将无法显示图像，而只是一副灰色的窗口。COn

7、VertSCaIeAbS()的原型为：dst = cv2.convertScaleAbs(src, dst, alpha, beta)其中可选参数alpha是伸缩系数,beta是加到结果上的一个值。结果返回uint8 类型的图片。由于Sobel算子是在两个方向计算的，最后还需要用cv2.addWeighted。.)函 数将其组合起来。其函数原型为：dst = cv2.addWeighted(src 1, alpha, src2, beta, gamma, dst, dtypeJ) 当然，这个函数我们在前面已经讲述过在这里就不多讲述了。 我们再来看代码：import cv2import nump

8、y as npimg = cv2.imread(,pie.png,0)x=cv2.Sobel(img,cv2.CV_ 16S, 1,0)y = cv2.Sobel(img,cv2.C V_ 16S,0,1 )absX=cv2.convertScaleAbs(x)absY =cv2.convertScaleAbs(y)dst=cv2.addWeighted(absX,0.5,absY,0.5,0)cv2.imshow(x,absX)cv2.imshow(,y,ab sY)cv2.imshow(,res,dst)cv2.waitKey(0)cv2.destroyAll Windows()先分别来看

9、下，x, y两个方向上的：再看看最终相加在一起的:由于Sobel算子是滤波算子的形式，用于提取边缘，可以利用快速卷积函数, 简单有效，因此应用广泛。美中不足的是，Sobel算子并没有将图像的主体与背 景严格地区分开来，即SobeI算子没有严格地模拟人的视觉生理特征，所以提取 的图像轮廓有时并不能令人满意。Scharr算子ScharrO函数提供了比标准SObeI函数更精确的计算结果。它使用了下面的 核：-3 O +3-3 -10 -3-10 O +10和Gy =OOO-3 O +3+3 +10 +3除了卷积核与SobeI不同，在其余方面它与SObel基本一致，我们来看它的 函数原型：dst =

10、cv2.Scharr(src, ddepth, dx, dy, dst, ksize, scale, delta, borderTypel) 参数就不再介绍了，与SObel是完全一致的，我们来看看代码演示：import cv2import numpy as npimg = cv2.imread(pie.png,0)x =cv2.Scharr(img,cv2.C V_ 16S, 1,0)y = cv2.Scharr(img,cv2.CV_ 16S,0,1 )absX =cv2.convertScaleAbs(x)absY = cv2.convertScaleAbs(y)dst =cv2.addW

11、eighted(absX,0.5,absY,0.5,0)cv2.imshow(,res,dst)cv2.waitKey(0)cv2. destroyAHWindows()本次教程所讲述的Sobel算子与Scharr算子在图形的边缘检测方面有些缺 陷，想必大家也看到了，不过在下面我们将会讲述这个问题。拉普拉斯算子我们在上面教程中的例子学习了使用Sobel边缘检测。原理是利用边缘区域 像素值的跳变。通过求一阶导数，可以使边缘值最大化。如下图所示：那么，如果求二阶导数会得到什么呢?可以观察到二阶导数为O的地方。因此,可以利用该方法获取图像中的边缘。 然而，需要注意的是二级导数为O的不只出现在边缘地方

12、，还可能是一些无意义 的位置，根据需要通过滤波处理该情况。二阶微分现在我们来讨论二阶微分，它是拉普拉斯算子的基础，与微积分中定义的微 分略有不同，数字图像中处理的是离散的值，因此对于一维函数的一阶微分的基 本定义是差值：空=(x+l)-(x)J +y-l可以用多种方式将其表示为数字形式。对于一个3*3的区域，一般情况下被 推荐最多的形式是：?/ =第 + 3 = /(z 1. !/) + f(x - 1,!/) + /(I,!/ + 1) f(x,!/ - 1) - 4(z, y)实现上式的滤波器模板为：O101-41010_我们可以发现，拉普拉斯算子不需要向Sobel算子那样分别对X, y方

13、向进 行处理，它可以直接处理，现在我们来看看OPenCV中的拉普拉斯算子的函数原 型：dst = cv2.Laplacian(src, ddepth, dst, ksize, scale, delta, borderType)(1)前两个是必须的参数：第一个参数是需要处理的图像。第二个参数是图像的深度，-1表示采用的是与原图像相同的深度。目标图像 的深度必须大于等于原图像的深度。(2)其后是可选的参数：dst不用解释了。ksize是算子的大小，必须为1、3、5、7o默认为1。SCale是缩放导数的比例常数，默认情况下没有伸缩系数。delta是一个可选的增量，将会加到最终的dst中，同样，默认情

14、况下没有额 外的值加到dst中。borderType是判断图像边界的模式。这个参数默认值为cv2.BORDER_DEFAULTo我们来看代码：import cv2import numpy as npimg = cv2.imread(,pie.png,)dst =cv2.Laplacian(img,cv2.CV_16S,ksize=3)dst =cv2.convertScaleAbs(dst)cv2.imshow(uimg,img)cv2.imshow(nres,dst)cv2.waitKey(0)cv2.destroyAllWindows()l现在可以拿这个结果对比上一个教程的结果了，我们发现，这个结果要比上 一个教程的结果好的多，对于边缘检测没有大的偏差。然而事实上，这只是对于简单的图像而言，而对于一幅复杂的图像，那么边 缘提取就