原理:
代码:
#include <opencv2\opencv.hpp>
using namespace std;
using namespace cv;
int main()
{
Mat img = imread("C:\\Users\\Yang\\Desktop\\1.bmp", 0);
imshow("原图",img);
int rows = img.rows;
int cols = img.cols;
Mat Dx = img.clone();//用于存放水平方向微分算子的运算结果
Mat Dy = img.clone();//用于存放垂直方向微分算子的运算结果
//卷积
//注意这里的i一定要从1开始,到ows-1结束。
for (int i = 1; i < rows-1; i++)
{
for (int j = 1; j < cols-1; j++)
{
Dx.at<uchar>(i, j) = img.at<uchar>(i + 1, j - 1) - img.at<uchar>(i - 1, j - 1)
+ 2 * (img.at<uchar>(i + 1, j) - img.at<uchar>(i - 1, j))
+ img.at<uchar>(i + 1, j + 1) - img.at<uchar>(i - 1, j + 1);
Dy.at<uchar>(i, j) = img.at<uchar>(i - 1, j + 1) - img.at<uchar>(i - 1, j - 1)
+ 2 * (img.at<uchar>(i, j + 1) - img.at<uchar>(i, j - 1))
+ img.at<uchar>(i + 1, j + 1) - img.at<uchar>(i + 1, j - 1);
}
}
Mat newImg = img.clone();
for (int i = 0; i < rows-1; i++)
{
for (int j = 0; j < cols-1; j++)
{
newImg.at<uchar>(i, j) = abs(Dx.at<uchar>(i, j)) + abs(Dy.at <uchar>(i, j));
}
}
imshow("效果图", newImg);
waitKey(0);
return 0;
}