《图像滤波研究.doc》由会员分享，可在线阅读，更多相关《图像滤波研究.doc（21页珍藏版）》请在金锄头文库上搜索。
1、太原理工大学现代科技学院 数字图像处理 课程设计 设计名称 图像滤波研究 专业班级 通信0903 学 号 姓 名 指导教师 太原理工大学现代科技学院 课程设计装订线专业班级 通信0903 学号 成绩 装订线 一、 目的与要求 提高分析问题、解决问题的能力，进一步巩固数字图像处理系统中的基本原理与方法。熟悉掌握一门计算机语言，可以进行数字图像的应用处理的开发设计。开发工具可以选用matlab、VC+、 VB、 C#等。本人选用matlab作为编程开发工具。二、 设计的内容刚获得的图像有很多的噪声，这是由平时的工作和环境引起的，图像增强是减弱噪音，增强对比度。想得到比较干净清晰的图像并不是容易的事
2、情。为这个目标而为处理图像所涉及的操作是设计一个适合、匹配的滤波器和恰当的阈值。这就是所谓的图像滤波研究。本设计就是用图像平均、平滑滤波、中值滤波、理想低通滤波方法、巴特沃斯滤波方法五种方法来实现图像滤波处理。三、 总体方案设计首先找一幅图像，然后对图像加噪声，获得加噪图像。然后分别选用以下五种方法进行图像滤波处理的研究。用图像平均的方法消除噪声并计算信噪比的改善。用平滑滤波方法消除噪声并计算信噪比的改善用中值滤波方法消除噪声并计算信噪比的改善用理想低通滤波方法消除噪声并计算信噪比的改善用巴特沃斯低通滤波方法消除噪声并计算信噪比的改善 装订线 然后改变图像的加噪类型，通过分析比较各种滤波结果，
3、总结各种滤波方法的适用范围。在该设计中，所采用的噪声类型主要有：高斯噪声，椒盐噪声。四、 各个功能模块的主要实现程序以及测试和调试1、 用图像平均的方法消除噪声并计算信噪比的改善。M文件：I=imread(circuit.tif);J=imnoise(I,gaussian,0,0.01);subplot(231),imshow(I);%title(yuanshituxiang)；subplot(232),imshow(J);%title(添加椒盐噪声的图像)；k1=filter2(fspecial(average,3),J);%进行3*3模板平滑滤波k2=filter2(fspecial(av
4、erage,5),J);%进行5*5模板平滑滤波k3=filter2(fspecial(average,7),J);%进行7*7模板平滑滤波k4=filter2(fspecial(average,9),J);%进行9*9模板平滑滤波subplot(233),imshow(uint8(k1);%title(3*3模板平滑滤波)；subplot(234),imshow(uint8(k2);%title(5*5模板平滑滤波)；subplot(235),imshow(uint8(k3);%title(7*7模板平滑滤波)；subplot(236),imshow(uint8(k4);%title(9*9
5、模板平滑滤波)；k=8;fmax=2k-1 装订线 s=0;d=0;for i=1:M for j=1:N y=double(I(i,j)-double(J(i,j); s=s+double(y2); endendPSNR1=10*LOG10(N*M*fmax2/s);i=0;j=0;for i=1:M for j=1:N y=double(I(i,j)-double(Q(i,j); d=d+double(y2); endendPSNR2=10*log10(N*M*fmax2/d);Y=PSNR2-PSNR1;Y该liangxujiao.m文件运行后的结果 装订线峰值信噪比PSNRY =11.
6、7799 装订线 2、用平滑滤波方法消除噪声并计算信噪比的改善M文件编写I=imread(circuit.tif);J=imnoise(I,gaussian,0,0.01);subplot(231),imshow(I);%title(yuanshituxiang)；subplot(232),imshow(J);%title(添加椒盐噪声的图像)；k1=filter2(fspecial(average,3),J);%进行3*3模板平滑滤波k2=filter2(fspecial(average,5),J);%进行5*5模板平滑滤波k3=filter2(fspecial(average,7),J);
7、%进行7*7模板平滑滤波k4=filter2(fspecial(average,9),J);%进行9*9模板平滑滤波subplot(233),imshow(uint8(k1);%title(3*3模板平滑滤波)；Q1=uint8(k1);subplot(234),imshow(uint8(k2);%title(5*5模板平滑滤波)；Q2=uint8(k2);subplot(235),imshow(uint8(k3);%title(7*7模板平滑滤波)；Q3=uint8(k2);subplot(236),imshow(uint8(k4);%title(9*9模板平滑滤波)；Q4=uint8(k2
8、);k=8;fmax=2k-1;s=0; 装订线 d=0;for i=1:M for j=1:N y=double(I(i,j)-double(J(i,j); s=s+double(y2); endendPSNR1=10*LOG10(N*M*fmax2/s);i=0;j=0;for i=1:M for j=1:N y=double(I(i,j)-double(Q1(i,j); d=d+double(y2); endendPSNR2=10*log10(N*M*fmax2/d);Y1=PSNR2-PSNR1;Y1i=0;j=0;for i=1:M 装订线 for j=1:N y=double(I(i,j)-double(Q(i,j); d=d+double(y2); endendPSNR3=10*log10(N*M*fmax2/d);Y2=PSNR3-PSNR1;Y2 i=0;j=0;for i=1:M for j=1:N y=double(I(i,j)-double(Q(i,j); d=d+double(y2); endendPSNR4=10*log10(N*M*fmax2/d);Y3=PSNR4-PSNR1;Y3i=0;j=0;for i=1:M 装订线 for j=1:N y=double(I(i,j)-double(Q(i,j);
