高斯滤波器的c++,matlab实现

最新推荐文章于 2025-06-30 14:24:44 发布

只取一勺

最新推荐文章于 2025-06-30 14:24:44 发布

阅读量1.1k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： C++ 文章标签：高斯滤波器

本文链接：https://blog.youkuaiyun.com/qq_41736190/article/details/83514482

C++ 专栏收录该内容

2 篇文章

订阅专栏

本文介绍了一个使用C++实现的高斯模糊矩阵生成算法。该算法首先创建一个动态二维数组，然后根据高斯分布公式计算每个元素的权重，并进行归一化处理，最后输出归一化后的矩阵。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

c++版本

#include <iostream>
using namespace std;
int gaussian(double **gaus, const int size, const double sigma);
int main()
{
	int size = 4;
	double **gaus = new double *[size];
	for (int i = 0; i < size; i++) {
		gaus[i] = new double[size];  //生成动态矩阵
	}
	gaussian(gaus,4, 1);
	return 0;
	 
}
int gaussian(double **gaus, const int size, const double sigma)
{
	const double PI = 4*atan(1.0);
	double radius;
	double center = (size+1) / 2;
	double sum = 0;
	for (int i = 0; i < size; i++)
	{
		for (int j = 0; j < size; j++) {
			radius = ((i + 1 - center)*(i + 1 - center) + (j + 1 - center)*(j + 1 - center));
			gaus[i][j] = (1/(2 * PI *sigma*sigma))* exp(-(radius / (2 * sigma*sigma)));
			sum += gaus[i][j];
		}
	}
	cout << sum << endl;
	for (int i = 0; i < size; i++) {
		for (int j = 0; j < size; j++) {
			gaus[i][j] /= sum;
			cout << gaus[i][j] << "	  ";
		}
		cout << endl;
	}
	return 0;
}

matlab版本

#include <iostream>
using namespace std;
int gaussian(double **gaus, const int size, const double sigma);
int main()
{
	int size = 4;
	double **gaus = new double *[size];
	for (int i = 0; i < size; i++) {
		gaus[i] = new double[size];  //生成动态矩阵
	}
	gaussian(gaus,4, 1);
	return 0;
	 
}
int gaussian(double **gaus, const int size, const double sigma)
{
	const double PI = 4*atan(1.0);
	double radius;
	double center = (size+1) / 2;
	double sum = 0;
	for (int i = 0; i < size; i++)
	{
		for (int j = 0; j < size; j++) {
			radius = ((i + 1 - center)*(i + 1 - center) + (j + 1 - center)*(j + 1 - center));
			gaus[i][j] = (1/(2 * PI *sigma*sigma))* exp(-(radius / (2 * sigma*sigma)));
			sum += gaus[i][j];
		}
	}
	cout << sum << endl;
	for (int i = 0; i < size; i++) {
		for (int j = 0; j < size; j++) {
			gaus[i][j] /= sum;
			cout << gaus[i][j] << "	  ";
		}
		cout << endl;
	}
	return 0;
}