梯度,sobel算子的理解

本文转于:https://blog.csdn.net/u014020344/article/details/50196635 若用侵权 请联系删除 ?

在这一讲中我们来学习一下opencv中最基本的边缘检测的知识,首先我们来介绍一下图像梯度

1.图像梯度

可以把图像看成二维离散函数,图像梯度其实就是这个二维离散函数的求导:
图像梯度:

梯度算子 Gradient operators

函数f(x,y)在(x,y)处的梯度为一个向量:

计算这个向量的大小为:

近似为:

梯度的方向角为:

可能上面的数学表达式大家觉得还是过于抽象,那么我们将图片看成是二维的离散函数,那么图像的梯度其实就是这个二维离散函数的求导(偏导数)

图像梯度: G(x,y) = dx i + dy j;
dx(i,j) = I(i+1,j) - I(i,j);
dy(i,j) = I(i,j+1) - I(i,j);
其中,I是图像像素的值(如:RGB值),(i,j)为像素的坐标。
图像梯度一般也可以用中值差分:
dx(i,j) = [I(i+1,j) - I(i-1,j)]/2;
dy(i,j) = [I(i,j+1) - I(i,j-1)]/2;

图像边缘一般都是通过对图像进行梯度运算来实现的。

图像梯度的最重要性质是,梯度的方向在图像灰度最大变化率上,它恰好可以反映出图像边缘上的灰度变化 
上面说的是简单的梯度定义,其实还有更多更复杂的梯度公式。

2.边缘检测

关于边缘的细节知识,请参考博主 http://blog.csdn.net/xiaowei_cqu/article/details/7829481

下面让我们来看一下边缘检测中常用的一些算子:

sobel算子边缘检测算法

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

Sobel卷积因子为:

具体计算如下:

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)点的灰度值;

 

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

其实在边缘检测中还会用到很多算子,如普利维特算子(Prewitt operate);罗伯茨交叉边缘检测(Roberts Cross operator);拉普拉斯算子等等。

2.opencv中Sobel算子实现边缘检测效果图

先上效果图:

image

image

当然,opencv中其实上面的代码完全可以用

void cvSobel( const CvArr* src, CvArr* dst, int xorder, int yorder, int aperture_size=3 );

这个函数来代替,相应的参数大家可以去百度,在这就不一一介绍了,本文主要针对sobel算子的检测原理来进行说明。

代码如下


 
 
  1. #include "cv.h"
  2. #include "highgui.h"
  3. #include "math.h"
  4. #include "stdio.h"
  5. #include "malloc.h"
  6. IplImage *image; //声明IplImage指针
  7. int height, width;
  8. CvScalar s;
  9. int sobel_y[ 9] = { 1, 2, 1, 0, 0, 0, -1, -2, -1 }; //y方向sobel算子
  10. int sobel_x[ 9] = { 1, 0, -1, 2, 0, -2, 1, 0, -1 }; //x方向sobel算子
  11. void sobel()
  12. {
  13. int i, j, k;
  14. int grayx = 0, grayy = 0, gray;
  15. int *data;
  16. int a[ 9];
  17. data = ( int *) malloc(height*width* sizeof( int));
  18. for (i = 0; i<height; i++)
  19. {
  20. for (j = 0; j<width; j++)
  21. {
  22. s = cvGet2D(image, i, j);
  23. gray = ( int)s.val[ 0];
  24. data[i*width + j] = gray;
  25. }
  26. }
  27. for (i = 1; i<height - 1; i++)
  28. {
  29. for (j = 1; j<width - 1; j++)
  30. {
  31. grayx = 0;
  32. grayy = 0;
  33. s = cvGet2D(image, i, j);
  34. a[ 0] = data[width*(i - 1) + j - 1];
  35. a[ 1] = data[width*(i - 1) + j];
  36. a[ 2] = data[width*(i - 1) + j + 1];
  37. a[ 3] = data[width*i + j - 1];
  38. a[ 4] = data[width*i + j];
  39. a[ 5] = data[width*i + j + 1];
  40. a[ 6] = data[width*(i + 1) + j - 1];
  41. a[ 7] = data[width*(i + 1) + j];
  42. a[ 8] = data[width*(i + 1) + j + 1];
  43. for (k = 0; k< 9; k++)
  44. grayy += a[k] * sobel_y[k];
  45. for (k = 0; k< 9; k++)
  46. grayx += a[k] * sobel_x[k];
  47. s.val[ 0] = ( abs(grayx) + abs(grayy)); //此处需要取绝对值
  48. cvSet2D(image, i, j, s);
  49. }
  50. }
  51. free(data);
  52. }
  53. int main(int argc, char** argv)
  54. {
  55. image = cvLoadImage( "image1.jpg", 0);
  56. IplImage* image1 = cvLoadImage( "image1.jpg", 1);
  57. height = image->height;
  58. width = image->width;
  59. sobel();
  60. cvNamedWindow( "Image", 1); //创建窗口
  61. cvNamedWindow( "Sobel", 1); //创建窗口
  62. cvShowImage( "Image", image1); //显示图像
  63. cvShowImage( "Sobel", image); //显示图像
  64. cvWaitKey( 0); //等待按键
  65. cvDestroyWindow( "Image"); //销毁窗口
  66. cvReleaseImage(&image); //释放图像
  67. return 0;
  68. }


参考文章:

http://blog.csdn.net/xiaowei_cqu/article/details/7829481

  • 1
    点赞
  • 15
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值