高斯滤波

#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