数字图像处理-实现理想、布特沃斯以及高斯低通滤波器

《数字图像处理-实现理想、布特沃斯以及高斯低通滤波器》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、02•理想低通滤波器04•高斯低通滤波器01•原图・・.aIIIIIIIIIt数字图像处理03作业1•实现理想、布特沃斯以及高斯低通滤波器【实验结果】03•布特沃斯低通滤波器图1如图1所示，02、03、04分别是01表示的原图在d二15时经过理想滤波器、2阶布特沃斯滤波器以及高斯滤波器之后形成的图像，由图可知，理想滤波器有振铃现象发生，平滑效果较为粗糙；2阶布特沃斯滤波器几乎没有振铃现象，图像较理想滤波器而言平滑效果也比较好；高斯滤波器没有振铃现象，其平滑效果是三者中最好的。【代码】cIearaIIIMO二imread('Fig0441(a)(cha

2、racterstestpattern).tif');subplot(2,2,1)imshow(IMO)title(f01-原图’)%理想低通滤波器d二15;IM1=fftshift(fft2(IMO));[m,n]=size(IM1);a1=round(m/2);b1二round(n/2);fori=1:mforj二1:ndistance=sqrt((i-a1)"2+(j-b1)"2);ifdistance〈二dh二1;eIseh二0;end;IM1(i,j)二h*IM1(i,j);end;end;IM1=uint8(real(ifft2(iffts

3、hift(IM1))));subplot(2,2,2)imshow(IMI)title(*02-理想低通滤波器')%2阶的布特沃斯低通滤波器k二2;IM2=fftshift(fft2(IMO));[m,n]二size(IM2);a2二floor(m/2);b2二floor(n/2);fori=1:mforj=1:ndistance二sqrt((i-a2)"2+(j-b2)八2);h二1/(1+(distance/d厂(2*k));IM2(i,j)二h*IM2(i,j);end;end;IM2=ifftshift(IM2);IM2=uint8(real

4、(ifft2(IM2)));subplot(2,2,3)imshow(IM2)titleC03-布特沃斯低通滤波器,)%高斯低通滤波器IMtemp1=im2double(IMO);[m,n]二size(IMO);fori=1:mforj二1:nI毗emp1(i,j)二(T厂(i+j)*IMtemp1(i,j);end;end;IMtemp2=fft2(IMtempl);[mfn]=size(lMtemp2);foru=1:mforv二1:nD(u,v)二[(u-m/2厂2+(v-n/2)^2]^0.5;factc=-D(u,v)"2/(2*d"2);

5、H(u,v)二exp(facte);G(u,v)二H(u,v)*IMtemp2(u,v);end;end;IMtemp3=ifft2(G);fori=1:mforj=1:nIM3(i,j)二(-1厂(i+j)*IMtemp3(i,j);end;end;IM3二real(IM3);subplot(2,2,4)imshow(IM3)title(*04-高斯低通滤波器')