C实现Sobel算子边缘检测

 

计算机中的图像是由点阵形式保存的。

通常，边缘检测是检测图像灰度值的不连续性，可以使用Sobel算子检测一阶导数，从而发现边缘。

Sobel算子是3*3矩阵:

分别对原图像进行滤波，分别得到水平和垂直边缘。

对矩阵使用滤波，要用到卷积运算。实现的方法有很多种，例如matlab的imfilter和conv2，opencv的filter2D等，甚至可以它们封装的Sobel进行检测，我在文章C简单实现矩阵相关和卷积中已经实现了一遍，这一次利用自己实现的卷积方法对图像分别进行水平和垂直的边缘检测。

回顾一下两个矩阵的卷积：

与之前实现稍微不同的是，为方便计算，传递的参数模板矩阵w变成了已知的二维数组，类型也变为int型

var **conv2(var **img, int w[3][3], int rows, int cols, int bm=3, int bn=3) {
	int n1 = rows + bm - 1;
	int n2 = cols + bn - 1;
	var **result = (var**)malloc(n1 * sizeof(var*));
	for (int i = 0; i < n1; i++) {
		result[i] = (var*)malloc(n2 * sizeof(var*));
		for (int j = 0; j < n2; j++) {
			int sum = 0;
			for (int m = 0; m < bm; m++) {
				for (int n = 0; n < bn; n+