图像变换 - sobel算子(cvSobel)

首先,我们来开一下计算机是如何检测边缘的。以灰度图像为例,它的理论基础是这样的,如果出现一个边缘,那么图像的灰度就会有一定的变化,为了方便假设由黑渐变为白代表一个边界,那么对其灰度分析,在边缘的灰度函数就是一个一次函数y=kx,对其求一阶导数就是其斜率k,就是说边缘的一阶导数是一个常数,而由于非边缘的一阶导数为零,这样通过求一阶导数就能初步判断图像的边缘了。通常是X方向和Y方向的导数,也就是梯度。理论上计算机就是通过这种方式来获得图像的边缘。

但是,具体应用到图像中你会发现这个导数是求不了的,因为没一个准确的函数让你去求导,而且计算机在求解析解要比求数值解麻烦得多,所以就想到了一种替代的方式来求导数。就是用一个3×3的窗口来对图像进行近似求导。拿对X方向求导为例,某一点的导数为第三行的元素之和减去第一行元素之和,这样就求得了某一点的近似导数。其实也很好理解为什么它就近似代表导数,导数就代表一个变化率,从第一行变为第三行,灰度值相减,当然就是一个变化率了。这就是所谓的Prewitt算子。这样近似X方向导数就求出来了。Y方向导数与X方向导数求法相似,只不过是用第三列元素之和减去第一列元素之和。X方向和Y方向导数有了,那么梯度也就出来了。这样就可以找出一幅图中的边缘了。 还有一个问题,由于求的是3×3中心点的导数,所以给第二列加了一个权重,它的权重为2,第一列和第三列的权重为1,好了,这就是Sobel算子了。相比Prewitt算子,Sobel

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
首先,需要导入必要的库和读取灰度图像。这里使用Python的OpenCV进行处理。 ```python import cv2 import numpy as np # 读取灰度图像 img = cv2.imread('image.jpg', 0) ``` 接下来,分别应用log算子sobel算子、prewitt算子和roberts算子进行边缘检测。 1. Log算子 Log算子是一种基于拉普拉斯变换的边缘检测算法,它可以增强图像的高频部分,使得边缘更加明显。 ```python # 定义log算子 log_kernel = np.array([[0, 0, -1, 0, 0], [0, -1, -2, -1, 0], [-1, -2, 16, -2, -1], [0, -1, -2, -1, 0], [0, 0, -1, 0, 0]]) # 应用log算子 log_img = cv2.filter2D(img, -1, log_kernel) # 显示边缘检测结果 cv2.imshow('Log Edge Detection', log_img) cv2.waitKey(0) ``` 2. Sobel算子 Sobel算子是一种基于图像梯度的边缘检测算法,它可以检测出图像中的水平和竖直边缘。 ```python # 定义sobel算子 sobel_x = np.array([[-1, 0, 1], [-2, 0, 2], [-1, 0, 1]]) sobel_y = np.array([[-1, -2, -1], [0, 0, 0], [1, 2, 1]]) # 应用sobel算子 sobel_x_img = cv2.filter2D(img, -1, sobel_x) sobel_y_img = cv2.filter2D(img, -1, sobel_y) sobel_img = cv2.bitwise_or(sobel_x_img, sobel_y_img) # 显示边缘检测结果 cv2.imshow('Sobel Edge Detection', sobel_img) cv2.waitKey(0) ``` 3. Prewitt算子 Prewitt算子也是一种基于图像梯度的边缘检测算法,它与Sobel算子类似,但是使用了不同的卷积核。 ```python # 定义prewitt算子 prewitt_x = np.array([[-1, 0, 1], [-1, 0, 1], [-1, 0, 1]]) prewitt_y = np.array([[-1, -1, -1], [0, 0, 0], [1, 1, 1]]) # 应用prewitt算子 prewitt_x_img = cv2.filter2D(img, -1, prewitt_x) prewitt_y_img = cv2.filter2D(img, -1, prewitt_y) prewitt_img = cv2.bitwise_or(prewitt_x_img, prewitt_y_img) # 显示边缘检测结果 cv2.imshow('Prewitt Edge Detection', prewitt_img) cv2.waitKey(0) ``` 4. Roberts算子 Roberts算子是一种基于图像梯度的边缘检测算法,它使用了不同的卷积核计算图像中的边缘。 ```python # 定义roberts算子 roberts_x = np.array([[0, 0, 0], [0, 1, 0], [0, 0, -1]]) roberts_y = np.array([[0, 0, 0], [0, 0, 1], [0, -1, 0]]) # 应用roberts算子 roberts_x_img = cv2.filter2D(img, -1, roberts_x) roberts_y_img = cv2.filter2D(img, -1, roberts_y) roberts_img = cv2.bitwise_or(roberts_x_img, roberts_y_img) # 显示边缘检测结果 cv2.imshow('Roberts Edge Detection', roberts_img) cv2.waitKey(0) ``` 最后,需要释放窗口和内存。 ```python # 释放窗口和内存 cv2.destroyAllWindows() ``` 完整代码如下:

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值