opencv017-Sobel算子

l卷积应用-图像边缘提取

l相关API

l代码演示

卷积应用-图像边缘提取

\

卷积应用-图像边缘提取

l边缘是什么 – 是像素值发生跃迁的地方，是图像的显著特征之一，在图像特征提取、对象检测、模式识别等方面都有重要的作用。

l如何捕捉/提取边缘 – 对图像求它的一阶导

  delta =  f(x) – f(x-1), delta越大，说明像素在X方向变化越大，边缘信号越强，

l我已经忘记啦，不要担心，用Sobel算子就好！卷积操作！

Sobel算子

l是离散微分算子（discrete differentiation operator），用来计算图像灰度的近似梯度

lSoble算子功能集合高斯平滑和微分求导

l又被称为一阶微分算子，求导算子，在水平和垂直两个方向上求导，得到图像X方法与Y方向梯度图像

Sobel算子

水平梯度:      垂直梯度

最终图像梯度

l求取导数的近似值，kernel=3时不是很准确，OpenCV使用改进版本Scharr函数，算子如下：

API说明cv::Sobel

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

)

API说明cv::Scharr

cv::Scharr (

InputArray Src // 输入图像

OutputArray dst// 输出图像，大小与输入图像一致

int depth // 输出图像深度.

Int dx.  // X方向，几阶导数

int dy // Y方向，几阶导数.

double scale  = 1

double delta = 0

int borderType = BORDER_DEFAULT

)

其它API

lGaussianBlur( src, dst, Size(3,3), 0, 0, BORDER_DEFAULT );

lcvtColor( src,  gray, COLOR_RGB2GRAY );

laddWeighted( A, 0.5,B, 0.5, 0, AB);

lconvertScaleAbs(A, B)// 计算图像A的像素绝对值，输出到图像B

防止数据溢出[0~255]

//高斯平滑-灰度处理-求梯度x和y-振幅图像
#include <opencv2/opencv.hpp>
#include <iostream>
using namespace std;
using namespace cv;

Mat src, dst, gray_src;

int main(int agrc, char** agrv) {

	src = imread("C:/Users/liyangxian/Desktop/bjl/nm3.jpg");
	if (!src.data) {
		printf("no load..\n");
		return -1;
	}
	const char* input_win = "input";
	const char* out_put = "out_put";
	namedWindow(input_win, CV_WINDOW_AUTOSIZE);
	namedWindow(out_put, CV_WINDOW_AUTOSIZE);
	imshow(input_win, src);
	GaussianBlur(src, dst, Size(3, 3), 0, 0);
	cvtColor(dst, gray_src, CV_BGR2GRAY);

	Mat xgrad, ygrad;
	Sobel(gray_src, xgrad, CV_16S, 1, 0, 3);
	Sobel(gray_src, ygrad, CV_16S, 0, 1, 3);
	convertScaleAbs(xgrad, xgrad);
	convertScaleAbs(ygrad, ygrad);
	imshow("xgrad", xgrad);
	imshow("ygrad", ygrad);

	Mat xygrad=Mat(xgrad.size(),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);
		}
	}
	//addWeighted(xgrad, 0.5, ygrad, 0.5, 0,xygrad);
	imshow(out_put, xygrad);
	waitKey(0);
	return 0;
}