勿忘初心做自己 发表于 2024-8-6 15:15:15

图像处置处罚 -- Sobel滤波器的实现原理与使用案例

Sobel滤波器

概述

Sobel滤波器是一种边沿检测方法,用于图像处置处罚和计算机视觉范畴。它通过计算图像灰度值的梯度来检测边沿。Sobel滤波器结合了高斯平滑和微分利用,以减少噪声并加强边沿检测结果。
实现原理

Sobel滤波器通过使用两个3x3卷积核(也称为掩模)来计算图像灰度值的水平和垂直梯度。分别称为                                              G                            x                                       G_x                  Gx​ 和                                              G                            y                                       G_y                  Gy​。
水平梯度核                                              G                            x                                       G_x                  Gx​

                                                    G                               x                                    =                                       [                                                                                                   −                                              1                                                                                                            0                                                                                             1                                                                                                                                 −                                              2                                                                                                            0                                                                                             2                                                                                                                                 −                                              1                                                                                                            0                                                                                             1                                                                                 ]                                          G_x = \begin{bmatrix} -1 & 0 & 1 \\ -2 & 0 & 2 \\ -1 & 0 & 1 \end{bmatrix}                     Gx​=               ​−1−2−1​000​121​               ​
垂直梯度核                                              G                            y                                       G_y                  Gy​

                                                    G                               y                                    =                                       [                                                                                                   −                                              1                                                                                                                            −                                              2                                                                                                                            −                                              1                                                                                                                                  0                                                                                             0                                                                                             0                                                                                                                   1                                                                                             2                                                                                             1                                                                                 ]                                          G_y = \begin{bmatrix} -1 & -2 & -1 \\ 0 & 0 & 0 \\ 1 & 2 & 1 \end{bmatrix}                     Gy​=               ​−101​−202​−101​               ​
梯度计算

对于每个像素                                    I                         (                         x                         ,                         y                         )                              I(x, y)                  I(x,y),应用这些卷积核以获得水平和垂直方向上的梯度值:
                                                    G                               x                                    (                            x                            ,                            y                            )                            =                                       ∑                                           i                                  =                                  −                                  1                                          1                                                 ∑                                           j                                  =                                  −                                  1                                          1                                    I                            (                            x                            +                            i                            ,                            y                            +                            j                            )                            ⋅                                       G                               x                                    (                            i                            +                            1                            ,                            j                            +                            1                            )                                  G_x(x, y) = \sum_{i=-1}^{1} \sum_{j=-1}^{1} I(x+i, y+j) \cdot G_x(i+1, j+1)                     Gx​(x,y)=i=−1∑1​j=−1∑1​I(x+i,y+j)⋅Gx​(i+1,j+1)
                                                    G                               y                                    (                            x                            ,                            y                            )                            =                                       ∑                                           i                                  =                                  −                                  1                                          1                                                 ∑                                           j                                  =                                  −                                  1                                          1                                    I                            (                            x                            +                            i                            ,                            y                            +                            j                            )                            ⋅                                       G                               y                                    (                            i                            +                            1                            ,                            j                            +                            1                            )                                  G_y(x, y) = \sum_{i=-1}^{1} \sum_{j=-1}^{1} I(x+i, y+j) \cdot G_y(i+1, j+1)                     Gy​(x,y)=i=−1∑1​j=−1∑1​I(x+i,y+j)⋅Gy​(i+1,j+1)
梯度幅值

然后,计算梯度幅值(也称为梯度强度):
                                       G                            =                                                                G                                     x                                     2                                              +                                             G                                     y                                     2                                                                   G = \sqrt{G_x^2 + G_y^2}                     G=Gx2​+Gy2​             ​
为了便于计算,也可以使用近似计算梯度幅值:
                                       G                            ≈                            ∣                                       G                               x                                    ∣                            +                            ∣                                       G                               y                                    ∣                                  G \approx |G_x| + |G_y|                     G≈∣Gx​∣+∣Gy​∣
梯度方向

                                       θ                            =                            arctan                            ⁡                                       (                                                      G                                     y                                                         G                                     x                                                      )                                          \theta = \arctan\left(\frac{G_y}{G_x}\right)                     θ=arctan(Gx​Gy​​)
使用场景


[*]边沿检测:Sobel滤波器常用于检测图像中的边沿,如在计算机视觉和图像处置处罚使命中的对象轮廓提取。
[*]特性提取:在图像处置处罚的前期阶段,边沿信息可以作为特性,用于后续的图像分析和识别使命。
[*]图像加强:通过突出表现图像中的边沿,可以加强图像的视觉结果,应用于图像加强和视觉体系中。
[*]目标检测与识别:在自动驾驶、机器人视觉和安防监控等范畴,通过边沿检测获取目标物体的轮廓信息。
代码实现

以下是使用Python和OpenCV实现Sobel滤波器的示例代码:
import cv2
import numpy as np
from matplotlib import pyplot as plt

# 读取图像并转换为灰度图
image = cv2.imread('image.jpg', cv2.IMREAD_GRAYSCALE)

# 使用OpenCV的Sobel函数计算梯度
sobel_x = cv2.Sobel(image, cv2.CV_64F, 1, 0, ksize=3)
sobel_y = cv2.Sobel(image, cv2.CV_64F, 0, 1, ksize=3)

# 计算梯度幅值
sobel = np.hypot(sobel_x, sobel_y)
sobel = np.uint8(sobel / np.max(sobel) * 255)

# 显示结果
plt.figure(figsize=(10, 5))
plt.subplot(1, 2, 1), plt.title('Original Image'), plt.imshow(image, cmap='gray')
plt.subplot(1, 2, 2), plt.title('Sobel Filtered Image'), plt.imshow(sobel, cmap='gray')
plt.show()

免责声明:如果侵犯了您的权益,请联系站长,我们会及时删除侵权内容,谢谢合作!更多信息从访问主页:qidao123.com:ToB企服之家,中国第一个企服评测及商务社交产业平台。
页: [1]
查看完整版本: 图像处置处罚 -- Sobel滤波器的实现原理与使用案例