【opencv练习20 - Sobel算子 一阶导数】

/*****************************************************
    测试程序 【Sobel边缘 一阶导数】      
    时间：2016年8月26日
******************************************************/
int main(void)
{

  Mat src, src_gray;
  Mat grad;
  const char* window_name = "Sobel Demo - Simple Edge Detector";
  int scale = 1;
  int delta = 0;
  int ddepth = CV_16S;

  // Load an image
  src = imread("YY01.jpg");


  //【1、高斯去噪】
  GaussianBlur( src, src, Size(3,3), 0, 0, BORDER_DEFAULT );        
  //【2、灰度化】
  cvtColor( src, src_gray, COLOR_RGB2GRAY );

  namedWindow( window_name, WINDOW_AUTOSIZE );

  //【3、生成梯度——X，Y】Generate grad_x and grad_y
  Mat grad_x, grad_y;
  Mat abs_grad_x, abs_grad_y;

  /// Gradient X【4、计算X梯度】
  //Scharr( src_gray, grad_x, ddepth, 1, 0, scale, delta, BORDER_DEFAULT );
  Sobel( src_gray, grad_x, ddepth, 1, 0, 3, scale, delta, BORDER_DEFAULT );
  convertScaleAbs( grad_x, abs_grad_x );

  /// Gradient Y【5、计算Y梯度】
  //Scharr( src_gray, grad_y, ddepth, 0, 1, scale, delta, BORDER_DEFAULT );
  Sobel( src_gray, grad_y, ddepth, 0, 1, 3, scale, delta, BORDER_DEFAULT );
  convertScaleAbs( grad_y, abs_grad_y );

  /// Total Gradient (approximate)【6、XY梯度合成】
  addWeighted( abs_grad_x, 0.5, abs_grad_y, 0.5, 0, grad );

  imshow( window_name, grad );
  imshow( "Disp YY（1）", src_gray );

  waitKey(0);
  return 0;
}