基本边缘检测算子—Roberts

#include<opencv2/opencv.hpp>
#include<iostream>
using namespace std;
using namespace cv;
Mat Robets(Mat &SrcImage)
{
	Mat DstImage = SrcImage.clone();
	int nRows = SrcImage.rows;
	int nCols = SrcImage.cols;
	for (int i = 0; i < nRows - 1; ++i)
	{
		for (int j = 0; j < nCols - 1; ++j)
		{
			int t1 = (SrcImage.at<uchar>(i, j) - SrcImage.at<uchar>(i + 1, j + 1))*
				(SrcImage.at<uchar>(i, j) - SrcImage.at<uchar>(i + 1, j + 1));
			int t2 = (SrcImage.at<uchar>(i + 1, j) - SrcImage.at<uchar>(i, j + 1))*
				(SrcImage.at<uchar>(i + 1, j) - SrcImage.at<uchar>(i, j + 1));
			DstImage.at<uchar>(i, j) = (uchar)(abs(t1)+ abs(t2));
			
		}
	}
	return DstImage;
}
int main()
{
	Mat SrcImage = imread("22.jpg");
	namedWindow("原图", WINDOW_NORMAL);
	imshow("原图", SrcImage);
	Mat GrayImage;
	cvtColor(SrcImage, GrayImage, COLOR_RGB2GRAY);
	Mat Image = Robets(GrayImage);
	namedWindow("edge", WINDOW_NORMAL);
	imshow("edge", Image);
	waitKey(0);
	return 0;

}

此图为用绝对值近似梯度幅值

*此图为用正常公式计算梯度幅值

其中注意代码for (int i = 0; i < nRows - 1; ++i)
{
for (int j = 0; j < nCols - 1; ++j)
要求i j 小于 ncols-1 否则用at访问像素时 指针越界
roberts算子是利用局部差分寻找边缘的一种算子
代码参考《opencv图像处理编程实例》