|
|
|
|
@ -0,0 +1,67 @@
|
|
|
|
|
#include <stdio.h>
|
|
|
|
|
#include <time.h>
|
|
|
|
|
|
|
|
|
|
// 定义GaussianBlur函数
|
|
|
|
|
void applyGaussianBlur(float src[5][5], float dst[5][5], int h, int w, float kernel[3][3]) {
|
|
|
|
|
int i, j, m, n;
|
|
|
|
|
float sum;
|
|
|
|
|
|
|
|
|
|
// 对图像进行高斯模糊
|
|
|
|
|
for (i = 1; i < h - 1; i++) { // 遍历除边界以外的所有像素
|
|
|
|
|
for (j = 1; j < w - 1; j++) {
|
|
|
|
|
sum = 0.0;
|
|
|
|
|
// 使用3×3高斯核对周围像素进行加权平均
|
|
|
|
|
for (m = -1; m <= 1; m++) {
|
|
|
|
|
for (n = -1; n <= 1; n++) {
|
|
|
|
|
sum += src[i + m][j + n] * kernel[m + 1][n + 1];
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
// 将计算结果存储到目标图像
|
|
|
|
|
dst[i][j] = sum;
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
int main() {
|
|
|
|
|
// 示例输入图像矩阵(5x5)
|
|
|
|
|
float src[5][5] = {
|
|
|
|
|
{1, 1, 1, 1, 1},
|
|
|
|
|
{1, 2, 2, 2, 1},
|
|
|
|
|
{1, 2, 4, 2, 1},
|
|
|
|
|
{1, 2, 2, 2, 1},
|
|
|
|
|
{1, 1, 1, 1, 1}
|
|
|
|
|
};
|
|
|
|
|
|
|
|
|
|
// 3x3高斯模糊核
|
|
|
|
|
float kernel[3][3] = {
|
|
|
|
|
{1.0 / 16.0, 2.0 / 16.0, 1.0 / 16.0},
|
|
|
|
|
{2.0 / 16.0, 4.0 / 16.0, 2.0 / 16.0},
|
|
|
|
|
{1.0 / 16.0, 2.0 / 16.0, 1.0 / 16.0}
|
|
|
|
|
};
|
|
|
|
|
|
|
|
|
|
// 输出图像矩阵
|
|
|
|
|
float dst[5][5] = {0};
|
|
|
|
|
|
|
|
|
|
clock_t start = clock();
|
|
|
|
|
|
|
|
|
|
// 调用GaussianBlur函数
|
|
|
|
|
applyGaussianBlur(src, dst, 5, 5, kernel);
|
|
|
|
|
|
|
|
|
|
clock_t end = clock();
|
|
|
|
|
|
|
|
|
|
// 输出结果
|
|
|
|
|
int i,j;
|
|
|
|
|
printf("Blurred Image:\n");
|
|
|
|
|
for (i = 0; i < 5; i++) {
|
|
|
|
|
for ( j = 0; j < 5; j++) {
|
|
|
|
|
printf("%.2f ", dst[i][j]);
|
|
|
|
|
}
|
|
|
|
|
printf("\n");
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
// 计算并输出运行时间
|
|
|
|
|
printf("Execution time: %.6f seconds\n", (double)(end - start) / CLOCKS_PER_SEC);
|
|
|
|
|
|
|
|
|
|
return 0;
|
|
|
|
|
}
|
|
|
|
|
|