#include<iostream>
#include<opencv2/opencv.hpp>
#include<opencv2/core/core.hpp>
#include<opencv2/imgproc/imgproc.hpp>
#include<opencv2/highgui/highgui.hpp>
using namespace std;
using namespace cv;
void myGaussianBlur(Mat SrcImage, Mat DstImage, int besarkernal, double delta)
{
int ker = besarkernal / 2;//定义高斯核半径
Mat Gaussian = Mat_<double>(besarkernal, besarkernal);//定义高斯核
double alpha = 1 / (2 * (22 / 7)*delta*delta);
//生成高斯核函数
for (int i = -ker; i < ker; i++)
{
for (int j = -ker; j < ker; j++)
{
Gaussian.at<double>(i + ker, j + ker) = exp(-(i*i + j*j) / 2 * delta*delta)*alpha;
}
}
DstImage = SrcImage.clone();
double pixel;//定义图像像素初始值
//遍历像素
for (int i = ker; i < SrcImage.rows - ker; ++i)
{
for (int j = ker; j < SrcImage.cols - ker; ++j)
{
pixel = 0;
for (int n = -ker; n < ker; ++n)
{
for (int k = -ker; k < ker; ++k)
{
pixel = SrcImage.at<uchar>(i + n, j + k)*Gaussian.at<double>(n + ker, k + ker);
}
}
DstImage.at<uchar>(i - ker, j - ker) = pixel;
}
}
namedWindow("滤波后图片", WINDOW_NORMAL);
imshow("滤波后图片", DstImage);
waitKey(0);
}
int main()
{
Mat image = imread("1.png");
namedWindow("原图", WINDOW_NORMAL);
imshow("原图", image);
Mat GrayImage;
cvtColor(image, GrayImage, COLOR_RGB2GRAY);
namedWindow("灰度图", WINDOW_NORMAL);
imshow("灰度图", GrayImage);
int k = 3;
double b = 0.8;
Mat out;
myGaussianBlur(GrayImage, out, k, b);
return 0;
}
本代码参考《Opencv图像处理编程实例》
高斯滤波详解参考 https://lps-683.iteye.com/blog/2251180
编程中遇到Mat.at问题参考https://blog.csdn.net/weixin_43213208/article/details/84617183