转自https://www.cnblogs.com/herenzhiming/articles/6526741.html  https://blog.csdn.net/qaz_wz/article/details/79052246

  • 算法原理

索贝尔算子(Sobel operator)主要用作边缘检测,在技术上,它是一离散性差分算子,用来运算图像亮度函数的灰度之近似值。在图像的任何一点使用此算子,将会产生对应的灰度矢量或是其法矢量

Sobel卷积因子为:

该算子包含两组3x3的矩阵,分别为横向及纵向,将之与图像作平面卷积,即可分别得出横向及纵向的亮度差分近似值。如果以A代表原始图像,Gx及Gy分别代表经横向及纵向边缘检测的图像灰度值,其公式如下:

具体计算如下:

Gx = (-1)*f(x-1, y-1) + 0*f(x,y-1) + 1*f(x+1,y-1)

+(-2)*f(x-1,y) + 0*f(x,y)+2*f(x+1,y)

+(-1)*f(x-1,y+1) + 0*f(x,y+1) + 1*f(x+1,y+1)

= [f(x+1,y-1)+2*f(x+1,y)+f(x+1,y+1)]-[f(x-1,y-1)+2*f(x-1,y)+f(x-1,y+1)]

Gy =1* f(x-1, y-1) + 2*f(x,y-1)+ 1*f(x+1,y-1)

+0*f(x-1,y) 0*f(x,y) + 0*f(x+1,y)

+(-1)*f(x-1,y+1) + (-2)*f(x,y+1) + (-1)*f(x+1, y+1)

= [f(x-1,y-1) + 2f(x,y-1) + f(x+1,y-1)]-[f(x-1, y+1) + 2*f(x,y+1)+f(x+1,y+1)]

其中f(a,b), 表示图像(a,b)点的灰度值;

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

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

如果梯度G大于某一阀值则认为该点(x,y)为边缘点。

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

Sobel算子根据像素点上下、左右邻点灰度加权差,在边缘处达到极值这一现象检测边缘。对噪声具有平滑作用,提供较为精确的边缘方向信息,边缘定位精度不够高。当对精度要求不是很高时,是一种较为常用的边缘检测方法。

附带知识:

普利维特算子(Prewitt operate)

除sobel边缘检测外 还有Prewitt算子, 它的卷积因子如下:

其他计算 和sobel差不多;

Prewitt算子利用像素点上下、左右邻点灰度差,在边缘处达到极值检测边缘。对噪声具有平滑作用,定位精度不够高。

罗伯茨交叉边缘检测(Roberts Cross operator

卷积因子如下:

灰度公式为:

近似公式为:

具体计算如下:

G(x,y)=abs(f(x,y)-f(x+1,y+1))+abs(f(x,y+1)-f(x+1,y))

灰度方向 计算公式为:

Roberts算子采用对角线方向相邻两像素之差近似梯度幅值检测边缘。检测水平和垂直边缘的效果好于斜向边缘,定位精度高,对噪声敏感

其他边缘检测技术:

参考文章:

http://homepages.inf.ed.ac.uk/rbf/HIPR2/featops.htm

http://homepages.inf.ed.ac.uk/rbf/HIPR2/sobel.htm

  • OpenCV实现

Sobel函数使用扩展的 Sobel 算子,来计算一阶、二阶、三阶或混合图像差分。

void Sobel (  
InputArray src,//输入图  
 OutputArray dst,//输出图  
 int ddepth,//输出图像的深度  
 int dx,  
 int dy,  
 int ksize=3,  
 double scale=1,  
 double delta=0,  
 int borderType=BORDER_DEFAULT );

列表内容
第一个参数,InputArray 类型的src,为输入图像,填Mat类型即可。
第二个参数,OutputArray类型的dst,即目标图像,函数的输出参数,需要和源图片有一样的尺寸和类型。
第三个参数,int类型的ddepth,输出图像的深度,支持如下src.depth()和ddepth的组合: 
若src.depth() = CV_8U, 取ddepth =-1/CV_16S/CV_32F/CV_64F 
若src.depth() = CV_16U/CV_16S, 取ddepth =-1/CV_32F/CV_64F 
若src.depth() = CV_32F, 取ddepth =-1/CV_32F/CV_64F 
若src.depth() = CV_64F, 取ddepth = -1/CV_64F
第四个参数,int类型dx,x 方向上的差分阶数。
第五个参数,int类型dy,y方向上的差分阶数。
第六个参数,int类型ksize,有默认值3,表示Sobel核的大小;必须取1,3,5或7。
第七个参数,double类型的scale,计算导数值时可选的缩放因子,默认值是1,表示默认情况下是没有应用缩放的。我们可以在文档中查阅getDerivKernels的相关介绍,来得到这个参数的更多信息。
第八个参数,double类型的delta,表示在结果存入目标图(第二个参数dst)之前可选的delta值,有默认值0。
第九个参数, int类型的borderType,我们的老朋友了(万年是最后一个参数),边界模式,默认值为BORDER_DEFAULT。这个参数可以在官方文档中borderInterpolate处得到更详细的信息。

调用Sobel函数

#include <opencv2/opencv.hpp>
#include <iostream>

using namespace cv;
using namespace std;

const char* path = "C:/Users/Administrator/Desktop/opencv/1.jpg";
const char* load_win = "load image";
const char* gray_win = "gray imgae";

int main(int argc, char** argv)
{
    Mat src, dst;
    src = imread(path);
    if (src.empty())
    {
        cout << "could not load..." << endl;
        return -1;
    }
    imshow(load_win, src);

GaussianBlur(src, dst, Size(3, 3), 0, 0);
    Mat src_gray;
    cvtColor(dst, src_gray, CV_BGR2GRAY);
    imshow(gray_win, src_gray);

Mat x_grad, y_grad;
    Sobel(src_gray, x_grad, CV_16S, 1, 0, 3);
    Sobel(src_gray, y_grad, CV_16S, 0, 1, 3);
    convertScaleAbs(x_grad, x_grad);
    convertScaleAbs(y_grad, y_grad);
    imshow("x_grad", x_grad);
    imshow("y_grad", y_grad);

waitKey(0);
    return 0;
}

Sobel边缘检测算法及OpenCV函数实现相关推荐

  1. matlab差分算子的灰度图像边缘检测,灰度图像的 Sobel 边缘检测算法的 HDL实现(一)...

    1.1 边缘检测算法介绍 所谓边缘是指其周围像素灰度急剧变化的那些象素的集合,它是图像最基本的特征.边缘存在于目标.背景和区域之间,所以,它是图像分割所依赖的最重要的依据.由于边缘是位置的标志,对灰度 ...

  2. Sobel边缘检测算法verilog实现及仿真

    实验verilog语言对sobel边缘检测算法进行设计 实验modelsim仿真工具进行仿真,程序和仿真截图如下图所示: 1.程序截图: 2.仿真截图: module edge_judge ( clk ...

  3. Python:实现sobel边缘检测算法(附完整源码)

    Python:实现sobel边缘检测算法 # coding=gbk import cv2 img = cv2.imread("1.jpg", 0) x = cv2.Sobel(im ...

  4. 基于MATLAB的Sobel边缘检测算法实现

    图像边缘就是图像灰度值突变的地方,也就是图像在该部分的像素值变化速度非常之快,就比如在坐标轴上一条曲线有刚开始的平滑突然来个大转弯,在变化出的导数非常大. Sobel算子主要用作边缘检测,它是一离散型 ...

  5. 【OpenCV】 - 图像分割之分水岭算法,watershed()函数的输出,对marker和image的改变

    一.背景 最近在学分水岭算法的opencv函数watershed()时,对函数执行完后image和marker的变化一无所知.懵懵懂懂. 于是便结合网上资料和自己现身说法,给大家分享一下[waters ...

  6. 【算法随记一】Canny边缘检测算法实现和优化分析。

    以前的博文大部分都写的非常详细,有很多分析过程,不过写起来确实很累人,一般一篇好的文章要整理个三四天,但是,时间越来越紧张,后续的一些算法可能就以随记的方式,把实现过程的一些比较容易出错和有价值的细节 ...

  7. canny边缘检测算法 opencv_OpenCV-Python Canny边缘检测 | 十九

    目标 在本章中,我们将学习 Canny边缘检测的概念 OpenCV函数: cv.Canny() 理论 Canny Edge Detection是一种流行的边缘检测算法.它由John F. Canny发 ...

  8. 边缘检测算法及各自优缺点

    边缘检测算法及各自优缺点 边缘提取其实也是一种滤波,不同的算子有不同的提取效果.比较常用的方法有三种,Sobel算子,Laplacian算子,Canny算子. Sobel算子检测方法对灰度渐变和噪声较 ...

  9. Matlab与FPGA图像处理系列——基于FPGA的实时边缘检测系统设计,sobel边缘检测流水线实现

    注:下载链接的资源是图片存 ROM 后读取进行 Sobel 检测显示在 VGA上,可供参考. 摘要:本文设计了一种基于 FPGA 的实时边缘检测系统,使用OV5640 摄像头模块获取实时的视频图像数据 ...

最新文章

  1. 小程序订单点击不同页面_小程序跳转页面参数丢失
  2. 黑马训练营自学笔记(03)
  3. linux内核编译及添加系统调用(hdu)_浅谈关于Linux内核write系统调用操作的原子性
  4. SQL Server (MSSQLSERVER) 启动又停止
  5. 变量、属性、函数、方法总结
  6. Vue之ElementUI导航菜单
  7. [MAC] 小技巧– 取消屏幕缩放功能,以免不小心误触
  8. opencv 基本绘图函数
  9. android蓝牙服务端设置,低功耗蓝牙BLE外围模式(peripheral)-使用BLE作为服务端
  10. python曲线拟合预测_python曲线拟合
  11. 计算机组成原理课程试题,计算机组成原理课程复习考试试题及答案B.doc
  12. 普林斯顿邓嘉学生亲述:一定要博士学位?不,我本科生也能在大厂当应用科学家...
  13. Rust更适合经验较少的程序员?
  14. 亚马逊6页纸开会方法!
  15. unity Font字体替换
  16. 以数字驱动未来,望美实业携手华为开启数字化转型新征程
  17. 字符串转json对象
  18. 2017阿里校招内推面试回忆
  19. python读取sav文件_在Python中读取SPSS(.sav)文件时,获取“title already used as a name or title”错误...
  20. 【牛投客】:让我们再来看看牛投客的战略网络

热门文章

  1. 程序员和美工是否可共存?
  2. 人生感悟:生活磨练有时也是一种财富
  3. [转载] Python-科赫雪花(科克曲线)
  4. 页面固定定位超出一屏
  5. P1160 队列安排 洛谷
  6. 【opencv入门篇】 10个程序快速上手opencv【上】
  7. iOS: AFNetworking手动配置(iOS7.1, AF2.2.4)
  8. 微信公众平台消息接口开发之校验签名与消息响应合并
  9. 杭州恒生数米基金网招聘1-3年本科.NET软件工程师
  10. 编写Oracle简单的存储过程