|
|
|
|
@ -0,0 +1,56 @@
|
|
|
|
|
#include <stdio.h>
|
|
|
|
|
#include <time.h>
|
|
|
|
|
|
|
|
|
|
// 应用3×3 GaussianBlur到给定的二维数组src上,并将结果存到二维数组dst中
|
|
|
|
|
void applyGaussianBlur(float src[][5], float dst[][5], int h, int w, float kernel[3][3]) {
|
|
|
|
|
// 遍历图像内部像素(忽略边界,只处理中心3×3区域对应的像素)
|
|
|
|
|
for (int i = 1; i < h - 1; i++) {
|
|
|
|
|
for (int j = 1; j < w - 1; j++) {
|
|
|
|
|
float sum = 0.0;
|
|
|
|
|
// 3×3 kernel卷积核加权平均计算
|
|
|
|
|
for (int m = 0; m < 3; m++) {
|
|
|
|
|
for (int n = 0; n < 3; n++) {
|
|
|
|
|
sum += kernel[m][n] * src[i - 1 + m][j - 1 + n];
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
dst[i][j] = sum;
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
int main() {
|
|
|
|
|
// 输入图像矩阵(这里示例一个5×5的矩阵)
|
|
|
|
|
float image[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}
|
|
|
|
|
};
|
|
|
|
|
// 常用的kernel_size=3的kernel
|
|
|
|
|
float kernel[3][3] = {
|
|
|
|
|
{1.0 / 16, 2.0 / 16, 1.0 / 16},
|
|
|
|
|
{2.0 / 16, 4.0 / 16, 2.0 / 16},
|
|
|
|
|
{1.0 / 16, 2.0 / 16, 1.0 / 16}
|
|
|
|
|
};
|
|
|
|
|
float result[5][5];
|
|
|
|
|
// 记录开始时间
|
|
|
|
|
clock_t start_time = clock();
|
|
|
|
|
applyGaussianBlur(image, result, 5, 5, kernel);
|
|
|
|
|
// 记录结束时间
|
|
|
|
|
clock_t end_time = clock();
|
|
|
|
|
double elapsed_time = (double)(end_time - start_time) / CLOCKS_PER_SEC;
|
|
|
|
|
|
|
|
|
|
// 输出dst[][]矩阵结果
|
|
|
|
|
printf("Gaussian Blur Result Matrix:\n");
|
|
|
|
|
for (int i = 0; i < 5; i++) {
|
|
|
|
|
for (int j = 0; j < 5; j++) {
|
|
|
|
|
printf("%.2f ", result[i][j]);
|
|
|
|
|
}
|
|
|
|
|
printf("\n");
|
|
|
|
|
}
|
|
|
|
|
// 输出运行时间
|
|
|
|
|
printf("Elapsed time: %.6f seconds\n", elapsed_time);
|
|
|
|
|
|
|
|
|
|
return 0;
|
|
|
|
|
}
|