图像处理Sobel 算子

Sobel算子是一种常用的图像处理算法,用于边缘检测。它利用了图像中灰度值的变化来识别边缘的位置。在OpenCV中,Sobel算子通常用于图像梯度计算,其中包括水平方向和垂直方向的梯度。这两个方向的梯度合并后可以得到图像的边缘信息。

原理:

Sobel算子的原理基于图像中的灰度变化。它使用了一个小的卷积核(3x3的矩阵)对图像进行卷积操作,以检测图像中像素值的梯度。对于水平方向和垂直方向的梯度,Sobel算子使用以下两个核:

水平方向的Sobel算子核:
SobelX = [ − 1 0 1 − 2 0 2 − 1 0 1 ] \text{SobelX} = \begin{bmatrix} -1 & 0 & 1 \\ -2 & 0 & 2 \\ -1 & 0 & 1 \end{bmatrix} SobelX= −1−2−1000121

垂直方向的Sobel算子核:
SobelY = [ − 1 − 2 − 1 0 0 0 1 2 1 ] \text{SobelY} = \begin{bmatrix} -1 & -2 & -1 \\ 0 & 0 & 0 \\ 1 & 2 & 1 \end{bmatrix} SobelY= −101−202−101

这两个核分别对图像进行水平和垂直方向的卷积运算,得到了水平方向和垂直方向上的梯度。

作用和适用场景:

  • 边缘检测:Sobel算子可以帮助识别图像中的边缘,因为边缘处的灰度值通常会发生剧烈变化。
  • 特征提取:梯度信息可以用于图像的特征提取,例如在计算机视觉任务中用于对象识别、图像分割等。

公式:

Sobel算子在水平方向(x轴)的应用:
G x = SobelX ∗ 图像 G_x = \text{SobelX} * \text{图像} Gx=SobelX∗图像

Sobel算子在垂直方向(y轴)的应用:
G y = SobelY ∗ 图像 G_y = \text{SobelY} * \text{图像} Gy=SobelY∗图像

最终梯度图像的合并:
梯度图像 = G x 2 + G y 2 \text{梯度图像} = \sqrt{{G_x}^2 + {G_y}^2} 梯度图像=Gx2+Gy2

代码示例:

python 复制代码
import cv2
import numpy as np

def show_images(image):
    cv2.namedWindow('image',cv2.WINDOW_KEEPRATIO)
    cv2.imshow('image',image)
    cv2.waitKey()
    cv2.destroyAllWindows()

def Sobel(image):
    # 使用Sobel算子计算水平方向和垂直方向的梯度
    sobel_x = cv2.Sobel(img, cv2.CV_64F, 1, 0, ksize=3)
    sobel_y = cv2.Sobel(img, cv2.CV_64F, 0, 1, ksize=3)
    # 计算梯度的幅值
    gradient_img = np.sqrt(sobel_x ** 2 + sobel_y ** 2)

    return sobel_x,sobel_y,gradient_img

if __name__ == '__main__':
    # 读取图像
    img = cv2.imread('cat-dog.png', flags=0)
    re_img=Sobel(img)
    # 将四张图像连接成一个大图像
    top_row = np.hstack((img, re_img[0]))
    bottom_row = np.hstack((re_img[1], re_img[2]))
    combined_img = np.vstack((top_row, bottom_row))

    show_images(combined_img)
相关推荐
心情好的小球藻26 分钟前
Python应用进阶DAY9--类型注解Type Hinting
开发语言·python
都叫我大帅哥27 分钟前
LangChain加载HTML内容全攻略:从入门到精通
python·langchain
惜.己37 分钟前
使用python读取json数据,简单的处理成元组数组
开发语言·python·测试工具·json
都叫我大帅哥2 小时前
Python的Optional:让你的代码优雅处理“空值”危机
python
曾几何时`4 小时前
基于python和neo4j构建知识图谱医药问答系统
python·知识图谱·neo4j
写写闲篇儿6 小时前
Python+MongoDB高效开发组合
linux·python·mongodb
杭州杭州杭州7 小时前
Python笔记
开发语言·笔记·python
路人蛃9 小时前
通过国内扣子(Coze)搭建智能体并接入discord机器人
人工智能·python·ubuntu·ai·aigc·个人开发
qiqiqi(^_×)9 小时前
卡在“pycharm正在创建帮助程序目录”
ide·python·pycharm
Ching·10 小时前
esp32使用ESP-IDF在Linux下的升级步骤,和遇到的坑Traceback (most recent call last):,及解决
linux·python·esp32·esp_idf升级