sobel算子

VS2017&OPENCV4.0

#include <opencv2/opencv.hpp>
#include <iostream>
#include <math.h>
#include <opencv2\imgproc\types_c.h>

using namespace cv;
int main(int argc, char** argv) {
	Mat src, dst;
	src = imread("D:\\test.jpg");
	if (!src.data) {
		printf("could not load image...\n");
		return -1;
	}

	char INPUT_TITLE[] = "input image";
	char OUTPUT_TITLE[] = "sobel-demo";
	namedWindow(INPUT_TITLE, WINDOW_AUTOSIZE);
	namedWindow(OUTPUT_TITLE, WINDOW_AUTOSIZE);
	imshow(INPUT_TITLE, src);

	Mat gray_src;
	GaussianBlur(src, dst, Size(3, 3), 0, 0);//高斯平滑
	cvtColor(dst, gray_src, CV_BGR2GRAY);//先转化为灰度图像
	imshow("gray image", gray_src);

	/*cv::Sobel (
	InputArray Src // 输入图像
		OutputArray dst// 输出图像，大小与输入图像一致
		int depth // 输出图像深度. 
		Int dx.  // X方向，几阶导数
		int dy // Y方向，几阶导数. 
		int ksize, SOBEL算子kernel大小，必须是1、3、5、7、
		double scale = 1
		double delta = 0
		int borderType = BORDER_DEFAULT
		)
		*/
	Mat xgrad, ygrad;
	Sobel(gray_src, xgrad, CV_16S, 1, 0, 3);//sobel的x方向的梯度
	Sobel(gray_src, ygrad, CV_16S, 0, 1, 3);//sobel的y方向的梯度
	//如果是函数Scharr,把对应的函数改变即可
	convertScaleAbs(xgrad, xgrad);//计算图像A的像素绝对值，输出到图像B
	convertScaleAbs(ygrad, ygrad);
	imshow("xgrad", xgrad);
	imshow("ygrad", ygrad);

	//下面计算x和y方向的梯度图像
	Mat xygrad = Mat(xgrad.size(), xgrad.type());
	printf("type : %d\n", xgrad.type());
	int width = xgrad.cols;
	int height = ygrad.rows;
	for (int row = 0; row < height; row++) {
		for (int col = 0; col < width; col++) {
			int xg = xgrad.at<uchar>(row, col);
			int yg = ygrad.at<uchar>(row, col);
			int xy = xg + yg;
			xygrad.at<uchar>(row, col) = saturate_cast<uchar>(xy);//限制xy 像素的范围
		}
	}

	//addWeighted(xgrad, 0.5, ygrad, 0.5, 0, xygrad);//xy权衡值为0.5的图像
	imshow(OUTPUT_TITLE, xygrad);

	waitKey(0);
	return 0;
}