基于matlab数字图像处理之低通滤波器Word文档下载推荐.doc

1、imshow（abs（s）,）;图像傅里叶变换所得频谱subplot（224）,imshow（log（abs（s）,）;图像傅里叶变换取对数所得频谱 a,b=size（s）;a0=round（a/2）;b0=round（b/2）;d=10;for i=1:afor j=1:b distance=sqrt（i-a0）2+（j-b0）2）; if distance0.2）=0;figuresurf（Hd,Facecolor,interpEdgecolornoneFacelightingphong% 画三维曲面（色）图1.1.2理想低通滤波器实践结果截图：1.2.1 Butterworth低通滤波器

2、实践代码：clear all;I1=imread（girl.bmpsubplot（221）,imshow（I1）;xlabel（a）原始图像f=double（I1）;%数据类型转换g=fft2（f）;%图像傅里叶转换g=fftshift（g）;%傅里叶变换平移F2=log（abs（g）;%对傅里叶变换结果取绝对值，然后取对数subplot（222）,imshow（F2,InitialMagnificationfit%将计算后的矩阵用图像表示colormap（jet）;%设置色彩索引图colorbar %显示色彩索引条（b）原始图像的傅里叶变换图像N1,N2=size（g）;%傅里叶变换图像尺寸

3、n=2;%参数赋初始值d0=5;n1=fix（N1/2）;%数据圆整n2=fix（N2/2）;N1 %遍历图像像素N2 d=sqrt（i-n1）2+（j-n2）2）;if d=0h=0;elseh=1/（1+（d/d0）（2*n）;endresult（i,j）=h*g（i,j）;% 图像矩阵计算处理F3=log（abs（result）; %对傅里叶变换结果取绝对值，然后取对数subplot（223）,imshow（F3,（c）滤波后的傅里叶变换图像）result=ifftshift（result）;X2=ifft2（result）;X3=uint8（real（X2）;subplot（224）,

4、imshow（X3）（d）Butterworth低通滤波图像f1,f2=freqspace（size（I1）,meshgridD=0.3;r=f1.2+f2.2;n=4;size（I1,1） for j=1:size（I1,2） t=r（i,j）/（D*D）; Hd（i,j）=1/（tn+1）; end1.2.2 Butterworth低通滤波器实践结果截图：1.3.1 高斯低通滤波器实践代码：IA=imread（f1,f2=freqspace（size（IA）,D=100/size（IA,1）;Hd=ones（size（IA）;size（IA,1）size（IA,2） Hd（i,j）=exp（-t）;Y=fft2（double（IA）;Y=fftshift（Y）;Ya=Y.*Hd;Ya=ifftshift（Ya）;Ia=real（ifft2（Ya）;subplot（2,2,1）,imshow（uint8（IA）;subplot（2,2,2）,imshow（uint8（Ia）;高斯低通滤波处理 % 画三维曲面（色）图1.3.2 高斯低通滤波器实践结果截图：