数字图像处理领域平均滤波中值滤波

1、东华大学实验报告课程 数字图像处理 名称 数字图像变换 实验名称: 邻域平均法(box模板)和中值滤波处理一、 实验目的图像变换是数字图像处理中的一种综合变换，如直方图变换、几何变换等。通过本实验，使得学生掌握两种变换的程序实现方法。二、 实验任务 请设计程序，分别用邻域平均法，其模板为：和中值滤波法对testnoise图像进行去噪处理（中值滤波的模板的大小也设为3×3）。三、实验环境本实验在Windows平台上进行，对内存及cpu主频无特别要求，使用VC或者MINGW（gcc）编译器均可。四、设计思路介绍代码的框架结构、所用的数据结构、各个类的介绍（类的功能、类中方法的功能、类的成

2、员变量的作用）、各方法间的关系试验要求中以给出大致的编程思路和源代码以及代码注释，只有黑框部分需要自己填写。在此不进行赘述。五、具体实现 实现设计思路中定义的所有的数据类型，对每个操作给出实际算法。对主程序和其他模块也都需要写出实际算法。 注意：源代码中要加上注释。代码：(红色为重点代码)<邻域平均法>(3*3)/*-利用第一次实验课提供的 dhc.h 和 dhc.c文件以获取位图的高 宽 以及从文件头到实际的位图数据的偏移字节数，从而实现对位图实际数据的操作。-*/#include <stdio.h>#include <stdlib.h>#include

3、<memory.h>#include "hdr.h"/*-定义结构指针-*/struct bmphdr *hdr; /定义用于直方图变量unsigned char *bitmap,*count,*new_color;/*-main()函数编写-*/int main() /定义整数 i， j 用于函数循环时的，nr_pixels为图像中像素的个数int i, j ,nr_pixels,nr_w,nr_h;/定义两个文件指针分别用于提取原图像的数据和生成直方图均衡化后的图像FILE *fp, *fpnew;/定义主函数的参数包括：输入的位图文件名和输出的位图文件名,

4、此处内容可以不要，在DOS下执行命令的时候再临时输入也可，为了方便演示，我这里直接把函数的参数确定了。/argc=3;/argv1="test.bmp"/argv2="testzf.bmp"/参数输入出错显示/*if (argc != 3) printf("please input the name of input and out bitmap filesn");exit(1);*/ 获取位图文件相关信息/hdr = get_header(argv1);hdr = get_header("testnoise.bmp&quo

5、t;);if (!hdr) exit(1);/以二进制可读方式打开输入位图文件fp = fopen("testnoise.bmp", "rb");if (!fp) printf("File open error!n");exit(1);/ 文件指针指向数据区域fseek(fp, hdr->offset, SEEK_SET);/计算位图像素的个数nr_pixels = hdr->width * hdr->height;nr_w = hdr->width;nr_h = hdr->height;bitmap =

6、 malloc(nr_pixels);new_color = malloc(nr_pixels);count = malloc(nr_w+2)*(+nr_h+2);/读取位图数据到bitmap中fread(bitmap, nr_pixels, 1, fp);fclose(fp);/因为图像边缘无法使用邻域平均，所以根据邻近颜色填补图像的周围一圈，存入count数组中/中心图像存入countfor(i=nr_w+3;i<(nr_w+2)*(+nr_h+2)-nr_w-3;i+)j=i/(nr_w+2);if(i%(nr_w+2)!=0&&(i+1)%(nr_w+2)!=0)

7、 counti=bitmapi-nr_w-1-2*j;/填补第一排for(i=1;i<nr_w+1;i+)counti=bitmapi-1;/填补最后一排for(i=1;i<nr_w+1;i+)count(nr_w+2)*(nr_h+1)+i=bitmapnr_w*(nr_h-1)+i-1;/填补左边一排for(i=0;i<nr_h+3;i+)counti*(nr_w+2)=counti*(nr_w+2)+1;/填补右边一排for(i=0;i<nr_h+3;i+)count(i+1)*(nr_w+2)-1=count(i+1)*(nr_w+2)-2;/邻域平均3*3fo

8、r(j=nr_w+3,i=0;j<(nr_w+2)*(+nr_h+2)-nr_w-3;j+)if(j%(nr_w+2)!=0&&(j+1)%(nr_w+2)!=0) new_colori=(countj+countj-1+countj+1+countj-nr_w-2+countj-1-nr_w-2+countj+1-nr_w-2+countj+nr_w+2+countj-1+nr_w+2+ countj+1+nr_w+2)/9,i+;/结果存入bitmap中for (i = 0; i < nr_pixels; i+)bitmapi=new_colori;/ 打开一个

9、以输出文件名命名的文件,设为可写的二进制形式fpnew = fopen("test_lynoise.bmp", "wb+");/由于位图文件的头部信息并没有因直方图均衡化而改变，因此输出图像的头部信息从原位图文件中拷贝即可：fwrite(hdr->signature, 2, 1, fpnew);fwrite(&hdr->size, 4, 1, fpnew);fwrite(hdr->reserved, 4, 1, fpnew);fwrite(&hdr->offset, 4, 1, fpnew);fwrite(&

10、;hdr->hdr_size, 4, 1, fpnew);fwrite(&hdr->width, 4, 1, fpnew);fwrite(&hdr->height, 4, 1, fpnew);fwrite(&hdr->nr_planes, 2, 1, fpnew);fwrite(&hdr->bits_per_pixel, 2, 1, fpnew);fwrite(&hdr->compress_type, 4, 1, fpnew);fwrite(&hdr->data_size, 4, 1, fpnew);fw

11、rite(&hdr->resol_hori, 4, 1, fpnew);fwrite(&hdr->resol_vert, 4, 1, fpnew);fwrite(&hdr->nr_colors, 4, 1, fpnew);fwrite(&hdr->important_color, 4, 1, fpnew);if (hdr->offset > 54)fwrite(hdr->info, (hdr->offset - 54), 1, fpnew);/直方图均衡化的数据(bitmap)赋值fwrite(bitmap, nr

12、_pixels, 1, fpnew);/关闭fclose(fpnew);/释放内存（优化程序必需）free(hdr);free(bitmap);free(new_color);free(count);return 0;<中值滤波>/*-利用第一次实验课提供的 dhc.h 和 dhc.c文件以获取位图的高 宽 以及从文件头到实际的位图数据的偏移字节数，从而实现对位图实际数据的操作。-*/#include <stdio.h>#include <stdlib.h>#include <memory.h>#include "hdr.h"

13、/*-定义结构指针-*/struct bmphdr *hdr; /定义用于直方图变量unsigned char *bitmap,*count,*new_color;/*-main()函数编写-*/int main() /定义整数 i， j 用于函数循环时的，nr_pixels为图像中像素的个数int i, j, m, n, nr_pixels,nr_w,nr_h,temp,t9;/定义两个文件指针分别用于提取原图像的数据和生成直方图均衡化后的图像FILE *fp, *fpnew;/定义主函数的参数包括：输入的位图文件名和输出的位图文件名,此处内容可以不要，在DOS下执行命令的时候再临时输入也可

14、，为了方便演示，我这里直接把函数的参数确定了。/argc=3;/argv1="test.bmp"/argv2="testzf.bmp"/参数输入出错显示/*if (argc != 3) printf("please input the name of input and out bitmap filesn");exit(1);*/ 获取位图文件相关信息/hdr = get_header(argv1);hdr = get_header("testnoise.bmp");if (!hdr) exit(1);/以二进制可

15、读方式打开输入位图文件fp = fopen("test.bmp", "rb");if (!fp) printf("File open error!n");exit(1);/ 文件指针指向数据区域fseek(fp, hdr->offset, SEEK_SET);/计算位图像素的个数nr_pixels = hdr->width * hdr->height;nr_w = hdr->width;nr_h = hdr->height;bitmap = malloc(nr_pixels);new_color = ma

16、lloc(nr_pixels);count = malloc(nr_w+2)*(+nr_h+2);/读取位图数据到bitmap中fread(bitmap, nr_pixels, 1, fp);fclose(fp);/因为图像边缘无法使用邻域平均，所以根据邻近颜色填补图像的周围一圈，存入count数组中/中心图像存入countfor(i=nr_w+3;i<(nr_w+2)*(+nr_h+2)-nr_w-3;i+)j=i/(nr_w+2);if(i%(nr_w+2)!=0&&(i+1)%(nr_w+2)!=0) counti=bitmapi-nr_w-1-2*j;/填补第一排

17、for(i=1;i<nr_w+1;i+)counti=bitmapi-1;/填补最后一排for(i=1;i<nr_w+1;i+)count(nr_w+2)*(nr_h+1)+i=bitmapnr_w*(nr_h-1)+i-1;/填补左边一排for(i=0;i<nr_h+3;i+)counti*(nr_w+2)=counti*(nr_w+2)+1;/填补右边一排for(i=0;i<nr_h+3;i+)count(i+1)*(nr_w+2)-1=count(i+1)*(nr_w+2)-2;/中值平均3*3for(j=nr_w+3,i=0;j<(nr_w+2)*(+nr

18、_h+2)-nr_w-3;j+)if(j%(nr_w+2)!=0&&(j+1)%(nr_w+2)!=0) t0=countj;t1=countj-1;t2=countj+1;t3=countj-nr_w-2;t4=countj-1-nr_w-2;t5=countj+1-nr_w-2;t6=countj+nr_w+2;t7=countj-1+nr_w+2;t8=countj+1+nr_w+2;for(m=0;m<9;m+)for(n=0;n<9-m;n+)if(tn>tn+1)temp=tn;tn=tn+1;tn+1=temp;new_colori=t4;i+;

19、/结果存入bitmap中for (i = 0; i < nr_pixels; i+)bitmapi=new_colori;/ 打开一个以输出文件名命名的文件,设为可写的二进制形式fpnew = fopen("test_zznoise.bmp", "wb+");/由于位图文件的头部信息并没有因直方图均衡化而改变，因此输出图像的头部信息从原位图文件中拷贝即可：fwrite(hdr->signature, 2, 1, fpnew);fwrite(&hdr->size, 4, 1, fpnew);fwrite(hdr->reser

20、ved, 4, 1, fpnew);fwrite(&hdr->offset, 4, 1, fpnew);fwrite(&hdr->hdr_size, 4, 1, fpnew);fwrite(&hdr->width, 4, 1, fpnew);fwrite(&hdr->height, 4, 1, fpnew);fwrite(&hdr->nr_planes, 2, 1, fpnew);fwrite(&hdr->bits_per_pixel, 2, 1, fpnew);fwrite(&hdr->compress_type, 4, 1, fpnew);fwrite(&hdr->data_size, 4, 1, fpnew);fwrite(&hdr->resol_hori, 4, 1, fpnew);fwrite(&hdr->resol_vert, 4, 1, fpnew);fwrite(&hdr->nr_colors, 4, 1, fpnew);fwrite(&hdr->important_color, 4, 1, fpnew);if (hdr