图像滤波处理_第1页
图像滤波处理_第2页
图像滤波处理_第3页
图像滤波处理_第4页
图像滤波处理_第5页
已阅读5页,还剩3页未读 继续免费阅读

下载本文档

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

文档简介

1、图像平滑滤波学号:A01314063 姓名:陈文鹏1.实验目的:分别使用空域法、频域法实现图像平滑滤波,分析空域滤波器和频域滤波器之间的转化过程2.实验方法和步骤:1.空域法图像平滑滤波:方法:邻域平均法步骤:1.imread()读取图像。2. rgb2gray()函数将原图转换为灰度图。3. imnoise()函数为原图加高斯白噪声。4. double()将加了噪声的图转换为双精度型。5.创建3*3模板,用模板对噪声图邻域平均,Image()显示处理后的图像。6.创建5*5模板,用模板对噪声图邻域平均,Image()显示处理后的图像。7.创建7*7模板,用模板对噪声图邻域平均,Image()

2、显示处理后的图像。其源程序代码如下:clear; clc; close all; I=imread('E:cwp.jpg'); %读取图像f=rgb2gray(I);%将原图转换为灰度图subplot(2,3,1);%分割2*3个窗口。取第一个窗口,下面在第一个窗口处显示图像imshow (f);%显示灰度图colormap(gray); title('原始图')%给显示的图像命名为“原始图”J=imnoise(f,'gaussian',0,0.09);%给原始图加入参数为0.09的高斯白噪声subplot(2,3,2);%分割2*3个窗口。取第

3、二个窗口imshow(J);%显示加了高斯白噪声的图像title('噪声图')%命名为“噪声图”m n=size(f); %获取灰度图的大小f=double(f);%转换f为双精度型c=1/9*1 1 1;1 1 1;1 1 1; %3*3模板for i=1:m for j=1:n L=f(i:i,j:j).*c; %求点积 G(i,j)=sum(sum(L);%求和 end endsubplot(2,3,3);image(G);%取第三个窗口title('3*3模板')%命名为“3*3模板”c=1/25*1 1 1 1 1;1 1 1 1 1;1 1 1 1

4、1; %5*5模板for i=1:m for j=1:n L=f(i:i,j:j).*c; %求点积 G(i,j)=sum(sum(L);%求和 end end subplot(2,3,4);image(G);%在第四个窗口中显示图像title('5*5模板')%命名为“5*5模板”c=1/49*1 1 1 1 1 1 1;1 1 1 1 1 1 1;1 1 1 1 1 1 1; %7*7模板for i=1:m for j=1:n L=f(i:i,j:j).*c; %求点积 G(i,j)=sum(sum(L);%求和 end end subplot(2,3,5);image(G

5、);%在第五个窗口中显示图像title('7*7模板')%命名为“7*7模板”end2.频域法图像平滑滤波:方法:二阶布特沃斯低通滤波法步骤:1.imread()读取图像。2.Size()函数获取图像矩阵大小。3.将矩阵转换为双精度型,再将原图转换为灰度图。4.wgn()函数产生高斯白噪声并加入原来图像中。5.fft2()对图像进行傅立叶变换。6.巴特沃斯低通滤波器传递函数与图像傅立叶变换相乘。7.ifft2()函数进行逆傅立叶变换。8.image()函数显示图像。其源程序代码如下:clear; clc; close all; I=imread('E:cwp.jpg&#

6、39;); %读取图像m n p=size(I);%获取图像矩阵大小I=double(I);%将原来的图像矩阵转换为双精度型的I=I(1:m,1:n,1).*0.3+I(1:m,1:n,2).*0.51+I(1:m,1:n,1).*0.11;%将原图转换为灰度图subplot(2,2,1);image(I);colormap(gray); %分割2*2个窗口。取第一个窗口,显示灰度图title('原始图');%命名为“原始图”Noise=wgn(m,n,25);%产生25dBm的高斯白噪声New = Noise + I;%将高斯白噪声与原始信号叠加subplot(2,2,2);

7、image(New);%在第二个窗口中显示加了噪声的图像title('噪声图')%命名为“噪声图”g=fft2(New); % 傅立叶变换g=fftshift(g); % 转换数据矩阵M,N=size(g); %获取矩阵大小nn=2; % 定义二阶巴特沃斯(Butterworth)低通滤波器d0=20; %截止频率为20 m=fix(M/2); n=fix(N/2);%圆整函数for i=1:M for j=1:N d=sqrt(i-m)2+(j-n)2); h=1/(1+0.414*(d/d0)(2*nn); % 计算低通滤波器传递函数 result(i,j)=h*g(i,j

8、);%传递函数与信号傅立叶变换函数相乘 end end result=ifftshift(result); J2=ifft2(result); %逆傅立叶变换J3=uint8(real(J2);%将输出图像转换为8位数据存储subplot(2,2,3);image(J3); %在第三个窗口中显示滤波后的图像title('截止频率20 ') ;%命名为“截止频率20”nn=2; % 定义二阶巴特沃斯(Butterworth)低通滤波器d0=100; %截止频率为100 m=fix(M/2); n=fix(N/2);%圆整函数for i=1:M for j=1:N d=sqrt(i

9、-m)2+(j-n)2); h=1/(1+0.414*(d/d0)(2*nn); % 计算低通滤波器传递函数 result(i,j)=h*g(i,j); %传递函数与信号傅立叶变换函数相乘 end end result=ifftshift(result); J2=ifft2(result);%二维逆傅立叶变换J3=uint8(real(J2); %转换为8位存储的图像subplot(2,2,4);image(J3); %在第四个窗口中显示绿波后的图像title('截止频率100 ')%命名为“截止频率100”end3.实验结果与分析:空域法图像平滑滤波处理结果如下:由噪声图与滤波后的图像对比可看出,邻域平均法对抑制噪声有明显的效果,但随着邻域的加大,就是随着模板的加

温馨提示

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

评论

0/150

提交评论