数字图像处理实验报告(邻域平均法和中值滤波法)_第1页
数字图像处理实验报告(邻域平均法和中值滤波法)_第2页
数字图像处理实验报告(邻域平均法和中值滤波法)_第3页
数字图像处理实验报告(邻域平均法和中值滤波法)_第4页
数字图像处理实验报告(邻域平均法和中值滤波法)_第5页
已阅读5页,还剩9页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、精品数字图像处理实验报告班级:姓名:学号:日期:感谢下载载邻域平均法和中值滤波处理一、实验目的图像变换是数字图像处理中的一种综合变换,如直方图变换、几何变换等。通过本实验,使得学生掌握两种变换的程序实现方法。二、实验任务请设计程序,分别用邻域平均法,其模板为:和中值滤波法对testnoise图像进行去噪处理(中值滤波的模板的大小也设为3X3)c三、实验环境本实验在Windows平台上进行,对内存及cpu主频无特别要求,使用VC或者MINGW(gcc)编译器均可。四、设计思路介绍代码的框架结构、所用的数据结构、各个类的介绍(类的功能、类中方法的功能、类的成员变量的作用)、各方法间的关系写。在此不

2、进行赘述。五、具体实现实现设计思路中定义的所有的数据类型,对每个操作给出实际算法。对主程序和其他模块也都需要写出实际算法。代码:<邻域平均法>(3*3)#include<stdio.h>#include<stdlib.h>#include<memory.h>#include"hdr.h"/*定义结构指针*/structbmphdr*hdr;/定义用于直方图变量unsignedchar*bitmap,*count,*new_color;/*main()函数编*/intmain()/定义整数i,j用于函数循环时的,nr_pixel

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

4、apfilesn");exit(1);*/获取位图文件相关信息/hdr=get_header(argv1);hdr=get_header("testnoise.bmp");if(!hdr)exit(1);/以二进制可读方式打开输入位图文件fp=fopen("testnoise.bmp","rb");if(!fp)printf("Fileopenerror!n");exit(1);/文件指针指向数据区域fseek(fp,hdr->offset,SEEK_SET);/计算位图像素的个数nr_pixels

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

6、-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;/填补第一排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&

7、lt;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_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+;/结果存入

8、bitmap中for(i=0;i<nr_pixels;i+;bitmapi=new_colori;/打开一个以输出文件名命名的文件,设为可写的二进制形式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(

9、&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->da

10、ta_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)赋值fwrite(bitmap,nr

11、_pixels,1,fpnew);/关闭fclose(fpnew);/释放内存(优化程序必需)free(hdr);free(bitmap);free(new_color);free(count);return0;<中值滤波>#include<stdio.h>#include<stdlib.h>#include<memory.h>#include"hdr.h"/*定义结构指针*/structbmphdr*hdr;/定义用于直方图变量unsignedchar*bitmap,*count,*new_color;/*main()函数编

12、写*/intmain()/定义整数i,j用于函数循环时的,nr_pixels为图像中像素的个数inti,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("pleasei

13、nputthenameofinputandoutbitmapfilesn");exit(1);*/获取位图文件相关信息/hdr=get_header(argv1);hdr=get_header("testnoise.bmp");if(!hdr)exit(1);/以二进制可读方式打开输入位图文件fp=fopen("test.bmp","rb");if(!fp)printf("Fileopenerror!n");exit(1);/文件指针指向数据区域fseek(fp,hdr->offset,SEEK_S

14、ET);/计算位图像素的个数nr_pixels=hdr->width*hdr->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数组中/中心图像存入countfor(i=nr_w+3;

15、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;/填补第一排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

16、+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_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

17、+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+;/结果存入bitmap中for(i=0;i<nr_pixels;i+)bitmapi=new_colori;/打开一个以输出文件名命名的文件,设为可写的二进制形式fpnew=fopen("test_zznoise.bmp","wb+");/由于位图文件的头部信息并没有因直方图均衡化而改变,因此输出图像的头部信

18、息从原位图文件中拷贝即可: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(&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. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论