Opencv-Sobel滤波器

Sobel用来求一阶导数。
Sobel滤波器原理看视觉笔记（二）
此例求了x和y方向的一阶导数。

#include<iostream>  
#include <opencv2/core/core.hpp>  
#include <opencv2/highgui/highgui.hpp>  
#include <opencv2/imgproc/imgproc.hpp>

using namespace std;
using namespace cv;

int scale = 1;
int delta = 0;
int ddepth = CV_16S;

int main()
{
    // 读入一张图片  
    Mat image = imread("car.png");
    //高斯滤波降噪
    GaussianBlur(image, image, Size(3, 3), 0, 0, BORDER_DEFAULT);
    //彩色图像转灰度图像
    cvtColor(image, image, COLOR_BGR2GRAY);
    Mat grad_x, grad_y;
    Mat abs_grad_x, abs_grad_y;
    /*
    Sobel( src_gray, grad_x, ddepth, 1, 0, 3, scale, delta, BORDER_DEFAULT );
    src_gray: 在本例中为输入图像，元素类型 CV_8U
    grad_x/grad_y: 输出图像.
    ddepth: 输出图像的深度，设定为 CV_16S 避免外溢。
    x_order: x 方向求导的阶数。
    y_order: y 方向求导的阶数。
    scale, delta 和 BORDER_DEFAULT: 使用默认值
    */
    //求 X方向梯度，竖直边缘
    Sobel(image, grad_x, ddepth, 1, 0, 3, scale, delta, BORDER_DEFAULT);
    //求 Y方向梯度，水平边缘
    Sobel(image, grad_y, ddepth, 0, 1, 3, scale, delta, BORDER_DEFAULT);
    //将中间结果转换到 CV_8U:
    convertScaleAbs(grad_x, abs_grad_x);
    imshow("x", abs_grad_x);
    convertScaleAbs(grad_y, abs_grad_y);
    imshow("y", abs_grad_y);
    //将两个方向的梯度相加来求取近似 
    addWeighted(abs_grad_x, 0.5, abs_grad_y, 0.5, 0, image);

    imshow("sobel", image);

    // 等待6000 ms后窗口自动关闭  
    waitKey(60000);
}