java 图像梯度检测_opencv学习笔记(六)---图像梯度

图像梯度的算法有很多方法：sabel算子，scharr算子，laplacian算子，sanny边缘检测(下个随笔)。。。

下面是我的一些理解：

sabel算子：

sobel算子主要用于获得数字图像的一阶梯度，常见的应用和物理意义是边缘检测。

函数：

Python: cv2.Sobel(src, ddepth, dx, dy[, dst[, ksize[, scale[, delta[, borderType]]]]]) → dst  (参数就不一一说了，常用的就那几个，其他默认即可)

Parameters:

src – Source image.

dst – Destination image of the same size and the same number of channels as src .

ddepth – Destination image depth.

xorder – Order of the derivative x.

yorder – Order of the derivative y.

ksize – Size of the extended Sobel kernel. It must be 1, 3, 5, or 7.

scale – Optional scale factor for the computed derivative values. By default, no scaling is applied. See getDerivKernels() for details.

delta – Optional delta value that is added to the results prior to storing them in dst .

borderType – Pixel extrapolation method. See borderInterpolate() for details.

原理

算子使用两个33的矩阵(图1)算子使用两个33的矩阵(图1)去和原始图片作卷积，分别得到横向Ｇ(ｘ)和纵向Ｇ(ｙ)的梯度值，如果梯度值大于某一个阈值，则认为该点为边缘点

Gx方向的相关模板：

Gy方向的相关模板：

具体计算如下：

图像的每一个像素的横向及纵向灰度值通过以下公式结合，来计算该点灰度的大小：

通常，为了提高效率使用不开平方的近似值：

#sobels算子

img = cv.imread("E:/pictures/lena.jpg",cv.IMREAD_UNCHANGED)

sobelx= cv.Sobel(img,cv.CV_64F,1,0) #cv.CV_64F将像素值转换为double型，不然计算后为负值的像素会被截断为0

sobelx = cv.convertScaleAbs(sobelx) #转换为uint8类型(x方向)

sobely = cv.Sobel(img,cv.CV_64F,0,1)

sobely = cv.convertScaleAbs(sobely)

sobelxy11 = cv.Sobel(img,cv.CV_64F,1,1) #直接x，y方向一起计算(效果不好，应分开计算，再求权重和)

sobelxy11 = cv.convertScaleAbs(sobelxy11)

sobelxy = cv.addWeighted(sobelx,0.5,sobely,0.5,0) #图像权重和

cv.imshow("orginal",img) #dst = cv.addWidget(src1,alpha,src2,beta,gamma)

cv.imshow("sobelx",sobelx) #src1 图一 alpha->图一的权重 src2->图二 beta->图二的权重 gamma->修正值

cv.imshow("sobely",sobely) #dst = src1*alpha+src2*beta+gamma

cv.imshow("sobelxy",sobelxy)

cv.imshow("sobelxy11",sobelxy11)

cv.waitKey()

cv.destroyAllWindows()

用sobel函数同时对想x，y方向检测的效果并不好，一般不用。

scharr算子：

函数：

Python: cv2.Scharr(src, ddepth, dx, dy[, dst[, scale[, delta[, borderType]]]]) → dst

Parameters:

src – Source image.

dst – Destination image of the same size and the same number of channels as src .

ddepth – Destination image depth.

xorder – Order of the derivative x.

yorder – Order of the derivative y.

scale – Optional scale factor for the computed derivative values. By default, no scaling is applied. See getDerivKernels() for details.

delta – Optional delta value that is added to the results prior to storing them in dst .

borderType – Pixel extrapolation method. See borderInterpolate() for details.

scharr算子是对sabel算子的增强，可以看到上图中很多细小的边缘都没检测到，那么scharr算子就是解决这个问题的，它比sabel算子更精确，速度和复杂程度却一样，只是因为用的核不一样

是scharr 的卷积核，他的原理和sabel算子一样。

#scharr算子 scharr算子是对sabel算子的增强 scharr算子等价于ksize=-1的sabel算子

img = cv.imread("E:/pictures/lena.jpg",cv.IMREAD_GRAYSCALE)

scharrx= cv.Scharr(img,cv.CV_64F,1,0) #scharrx算子要满足dx>=0&&dy>=0&&dx+dy=1

scharrx = cv.convertScaleAbs(scharrx)

scharry = cv.Scharr(img,cv.CV_64F,0,1)

scharry = cv.convertScaleAbs(scharry)

scharrxy = cv.addWeighted(scharrx,0.5,scharry,0.5,0)

#sabel算子和 scharr算子的比较

sobelx= cv.Sobel(img,cv.CV_64F,1,0)

sobelx = cv.convertScaleAbs(sobelx)

sobely = cv.Sobel(img,cv.CV_64F,0,1)

sobely = cv.convertScaleAbs(sobely)

sobelxy = cv.addWeighted(sobelx,0.5,sobely,0.5,0)

cv.imshow("orginal",img)

cv.imshow("sobelxy",sobelxy)

cv.imshow("scharrxy",scharrxy)

cv.waitKey()

cv.destroyAllWindows()

laplacian算子：

Laplace算子和Sobel算子一样，属于空间锐化滤波操作，只不过是用的二阶微分，看官网的一些解释：

The function calculates the Laplacian of the source image by adding up the second x and y derivatives calculated using the Sobel operator:

This is done when ksize > 1 . When ksize == 1 , the Laplacian is computed by filtering the image with the following

aperture:

cv.Laplace(src, dst, ksize=3) → None

Parameters:

src – Source image.

dst – Destination image of the same size and the same number of channels as src .

ddepth – Desired depth of the destination image.

ksize – Aperture size used to compute the second-derivative filters. See getDerivKernels() for details. The size must be positive and odd.

scale – Optional scale factor for the computed Laplacian values. By default, no scaling is applied. See getDerivKernels() for details.

delta – Optional delta value that is added to the results prior to storing them in dst .

borderType – Pixel extrapolation method. See borderInterpolate() for details.

#拉普拉斯算子

img = cv.imread("E:/pictures/erode1.jpg",cv.IMREAD_GRAYSCALE)

r = cv.Laplacian(img,cv.CV_64F)

r = cv.convertScaleAbs(r) #关键代码就这两行，拉普拉斯算子不用再求x，y的权重和，因为这个函数都计算好了

cv.imshow("orginal",img)

cv.imshow("laplacian",r)

cv.waitKey()

cv.destroyAllWindows()

Canny边缘检测：

Python: cv2.Canny(image, threshold1, threshold2[, edges[, apertureSize[, L2gradient]]]) → edges

threshold1和threshold2是两个阈值，越小检测效果越好

import cv2 as cv

import numpy as np

img = cv.imread("E:/pictures/lena.jpg")

result1 = cv.Canny(img,50,100)

result2 = cv.Canny(img,100,200)

cv.imshow("orginal",img)

cv.imshow("result1",result1)

cv.imshow("result2",result2)

cv.waitKey()

cv.destroyAllWindows()