Sobel kernel 简介

本文介绍索伯算子(Sobel operator),一种在图像处理和计算机视觉中用于边缘检测的重要工具。文章详细解释了如何使用索伯算子进行横向和纵向的边缘检测,并通过数学公式展示了梯度大小和方向的计算方法。

索伯算子(Sobel operator)是图像处理中的算子之一,有时又称为索伯-费德曼算子索贝滤波器,在影像处理电脑视觉领域中常被用来做边缘检测。Sobel 算子是一个离散微分算子 (discrete differentiation operator)

如果以{\displaystyle \mathbf {A} }\mathbf{A}代表原始图像,{\displaystyle \mathbf {G_{x}} }{\mathbf  {G_{x}}}及{\displaystyle \mathbf {G_{y}} }{\mathbf  {G_{y}}}分别代表经横向及纵向边缘检测的图像,其公式如下:

{\displaystyle \mathbf {G_{x}} ={\begin{bmatrix}+1&0&-1\\+2&0&-2\\+1&0&-1\end{bmatrix}}*\mathbf {A} \quad {\mbox{and}}\quad \mathbf {G_{y}} ={\begin{bmatrix}+1&+2&+1\\0&0&0\\-1&-2&-1\end{bmatrix}}*\mathbf {A} }{\displaystyle \mathbf {G_{x}} ={\begin{bmatrix}+1&0&-1\\+2&0&-2\\+1&0&-1\end{bmatrix}}*\mathbf {A} \quad {\mbox{and}}\quad \mathbf {G_{y}} ={\begin{bmatrix}+1&+2&+1\\0&0&0\\-1&-2&-1\end{bmatrix}}*\mathbf {A} }

图像的每一个像素的横向及纵向梯度近似值可用以下的公式结合,来计算梯度的大小。

{\displaystyle \mathbf {G} ={\sqrt {\mathbf {G_{x}} ^{2}+\mathbf {G_{y}} ^{2}}}}{\mathbf  {G}}={\sqrt  {​{\mathbf  {G_{x}}}^{2}+{\mathbf  {G_{y}}}^{2}}}

然后可用以下公式计算梯度方向。

{\displaystyle \mathbf {\Theta } =\operatorname {arctan} \left({\mathbf {G_{y}} \over \mathbf {G_{x}} }\right)}{\mathbf  {\Theta }}=\operatorname {arctan}\left({​{\mathbf  {G_{y}}} \over {\mathbf  {G_{x}}}}\right)

以纵向边缘为例,如果角度{\displaystyle \Theta }\Theta等于零,即代表图像的纵向边缘右方较亮;若为 {\displaystyle \pi }\pi,则左方较亮。

 

例子:

 

参考:

维基、百度

https://blog.youkuaiyun.com/qq_37124237/article/details/82183177

module ov7670_driver ( input wire clk, // 系统时钟 input wire rst_n, // 复位信号,低电平有效 // OV7670 相关信号 output reg pwdn, // 电源控制信号 output reg reset, // 复位信号 output reg scl, // SCCB 时钟信号 inout scd, // SCCB 数据信号 output reg vcm_en, // 闪光灯控制信号 output reg ov7670_pclk, // 像素时钟信号 output reg ov7670_vsync, // 垂直同步信号 output reg ov7670_href, // 行同步信号 input wire [7:0] ov7670_dout // 数据输出信号 ); // 初始化状态 reg init_state; // SCCB 通信控制 reg scb_wr_en; reg scb_rd_en; reg [7:0] scb_wr_data; reg [6:0] scb_reg_addr; reg scb_busy; // 数据采集控制 reg capture_en; reg [15:0] line_count; reg [15:0] pixel_count; // 控制信号生成 reg pclk_div; assign ov7670_pclk = pclk_div; // 图像处理相关 reg [7:0] image_buffer [0:307199]; // 假设分辨率为 640x480 reg [7:0] processed_image_buffer [0:307199]; reg [7:0] edge_buffer [0:307199]; reg [15:0] object_width; reg [15:0] object_height; reg [15:0] object_depth; reg processing_image; reg [7:0] Sobel_kernel_x [0:2][0:2]; reg [7:0] Sobel_kernel_y [0:2][0:2]; // 初始化 Sobel 算子 initial begin Sobel_kernel_x[0][0] = -1; Sobel_kernel_x[0][1] = 0; Sobel_kernel_x[0][2] = 1; Sobel_kernel_x[1][0] = -2; Sobel_kernel_x[1][1] = 0; Sobel_kernel_x[1][2] = 2; Sobel_kernel_x[2][0] = -1; Sobel_kernel_x[2][1] = 0; Sobel_kernel_x[2][2] = 1; Sobel_kernel_y[0][0] = 1; Sobel_kernel_y[0][1] = 2; Sobel_kernel_y[0][2] = 1; Sobel_kernel_y[1][0] = 0; Sobel_kernel_y[1][1] = 0; Sobel_kernel_y[1][2] = 0; Sobel_kernel_y[2][0] = -1; Sobel_kernel_y[2][1] = -2; Sobel_kernel_y[2][2] = -1; end // 初始化过程 always @(posedge clk or negedge rst_n) begin if (!rst_n) begin // 初始化代码保持不变 end else begin // 初始化 OV7670 代码保持不变 // 图像处理部分 if (processing_image) begin // 调用图像处理任务 process_image; end end end // 图像处理任务 task process_image; reg [15:0] i, j; reg [15:0] width, height, depth; begin // 边缘检测 for (i = 1; i < 479; i = i + 1) begin for (j = 1; j < 639; j = j + 1) begin // 计算梯度,代码省略具体实现 // 根据梯度确定边缘 end end // 物体识别和尺寸测量 // 代码省略具体实现,通常会使用连通域分析等算法 // 假设识别出物体后得到宽度、高度 object_width = width; object_height = height; // 假设深度通过其他方法得到 // object_depth = depth; end endtask // 其他模块代码保持不变 endmodule分析一下这串代码的功能
07-09
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值