#include<stdio.h>
#include<time.h>
void applyGaussianBlur(float src[5][5], float dst[5][5], int h, int w, float kernel[3][3])
{
	for (int i = 1; i < h - 1; i++)
	{
		for (int j = 1; j < w - 1; j++)
		{
			float sum = 0.0;
			for (int ki = -1; ki <= 1; ki++)
			{
				for (int kj = -1; kj <= 1; kj++)
				{
					sum += src[i + ki][j + kj] * kernel[ki + 1][kj + 1];
				}
			}
			dst[i][j] = sum;
		}
	}
}
int main()
{
	float src[5][5] = {{1,2,3,4,5},{6,7,8,9,10},{11,12,13,14,15},{16,17,18,19,20},{21,22,23,24,25}};
	float dst[5][5] = { 0 };
	float kernel[3][3] ={{1/16.0,2/16.0,1/16.0},{2/16.0,4/16.0,2/16.0},{1/16.0,2/16.0,1/16.0}};
	clock_t start = clock();
	applyGaussianBlur(src, dst, 5, 5, kernel);
	clock_t end = clock();
	printf("模糊后的图像矩阵:\n");
	for (int i = 0; i <= 5; i++)
	{
		for (int j = 0; j < 5; j++)
		{
			printf("%.2f", dst[i][j]);
		}
		printf("\n");
	}
	double time_taken = ((double)(end - start)) / CLOCKS_PER_SEC;
	printf("运行时间:%f 秒\n", time_taken);
	return 0;
}