数字图像处理实验报告材料(邻域平均法和中值滤波法

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

2、功能、类的成员变量的作用）、各方法间的关系写。在此 不进行赘述。五、具体实现实现设计思路中定义的所有的数据类型，对每个操作给出实际算法 对主程序和其他模块也都需要写出实际算法。代码：(3*3)#include #include #include #include hdr.h /* 定义结构指针 */struct bmphdr *hdr;/ 定义用于直方图变量unsigned char *bitmap,*count,*new_color; /*main() 函 数编*/int main() / 定义整数 i， j 用于函数循环时的， nr_pixels 为图像中像素的个 数int i, j ,n

3、r_pixels,nr_w,nr_h; / 定义两个文件指针分别用于提取原图 的数据和生成直方图均衡化后的图像FILE *fp, *fpnew; / 定义主函数的参数包括：输入的位图文件名和 输出的位图文件名 ,此处内容可以不要， 在 DOS 下执行命令的时候再 临时输入也可， 为了方便演示， 我这里直接把函数的参数确定了。 / argc=3; /argv1=test.bmp; / argv2=testzf.bmp; / 参数输入出错显示/* if (argc != 3) printf(please input the name of input and out bitmap filesn);

4、exit(1);*/ 获取位图文件相关信息 /hdr = get_header(argv1);hdr = get_header(testnoise.bmp);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;n

5、r_h = hdr-height;bitmap = 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 数组中/ 中心图像存入 count for(i=nr_w+3;i(nr_w+2)*(+nr_h+2)-nr_w-3;i+) j=i/(nr_w+2);if(i%(nr_w+

6、2)!=0&(i+1)%(nr_w+2)!=0) counti=bitmapi-nr _w-1-2*j;/ 填补第一排 for(i=1;inr_w+1;i+) counti=bitmapi-1; / 填补最后一排 for(i=1;inr_w+1;i+) count(nr_w+2)*(nr_h+1)+i=bitmapnr_w*(nr_h-1)+i-1; / 填补左边一排for(i=0;inr_h+3;i+) counti*(nr_w+2)=counti*(nr_w+2)+1; / 填补右边一排for(i=0;inr_h+3;i+) count(i+1)*(nr_w+2)-1=count(i+1)

7、*(nr_w+2)-2;/ 邻域平均 3*3for( 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 signature, 2, 1, fpnew);fwri

8、te(&hdr-size, 4, 1, fpnew); fwrite(hdr-reserved, 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_

9、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-offset 54) fwrite(hdr-info, (hdr-offset - 54), 1, fpnew); /直 方图 均衡 化的数据 (bitmap

10、) 赋值 fwrite(bitmap, nr_pixels, 1, fpnew);/ 关闭 fclose(fpnew);/ 释放内存(优化程序必需) free(hdr);free(bitmap);free(new_color);free(count); return 0;#include #include #include #include hdr.h/* 定义结构指针 */struct bmphdr *hdr;/ 定义用于直方图变量 unsigned char *bitmap,*count,*new_color;/*main() 函数编写 */int main() / 定义整数 i ， j

11、用于函数循环时的， nr_pixels 为图像中像素的 个数int i, j, m, n, nr_pixels,nr_w,nr_h,temp,t9; / 定义两个文件指针 分别用于提取原图像的数据和生成直方图均衡化后的图像FILE *fp, *fpnew; / 定义主函数的参数包括：输入的位图文件名和 输出的位图文件名 ,此处内容可以不要， 在 DOS 下执行命令的时候再 临时输入也可 ./ argc=3;/ argv1=test.bmp;/ argv2=testzf.bmp;/ 参数输入出错显示/* if (argc != 3) printf(please input the name of

12、 input and out bitmap filesn );exit(1);*/ 获取位图文件相关信息 /hdr = get_header(argv1);hdr = get_header(testnoise.bmp);if (!hdr) exit(1); / 以二进制可读方式打开输入位图文件 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

13、-height;nr_w = hdr-width;nr_h = hdr-height; bitmap = 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 数组中/ 中心图像存入 count for(i=nr_w+3;i(nr_w+2)*(+nr_h+2)-nr_w-

14、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; / 填补第一排for(i=1;inr_w+1;i+) counti=bitmapi-1; / 填补最后一排for(i=1;inr_w+1;i+) count(nr_w+2)*(nr_h+1)+i=bitmapnr_w*(nr_h-1)+i-1; / 填补左边一排for(i=0;inr_h+3;i+) counti*(nr_w+2)=counti*(nr_w+2)+1; / 填补右边一排 for(i=0;inr_h+3;i+) co

15、unt(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_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;m9;m+)

16、for(n=0;ntn+1) temp=tn;tn=tn+1;tn+1=temp;new_colori=t4;i+; / 结果存入 bitmap 中for (i = 0; i signature, 2, 1, fpnew);fwrite(&hdr-size, 4, 1, fpnew);fwrite(hdr-reserved, 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);