基于matlab數(shù)字圖像處理之低通濾波器

1、實踐一:理想低通濾波器、Butterworth低通濾波器、高斯低通濾波器1.1.1理想低通濾波器實踐代碼： I=imread(couple.bmp);%I=rgb2gray(I); subplot(221),imshow(I);title(原圖像);s=fftshift(fft2(I);subplot(223),imshow(abs(s),);title(圖像傅里葉變換所得頻譜);subplot(224),imshow(log(abs(s),); title(圖像傅里葉變換取對數(shù)所得頻譜); a,b=size(s);a0=round(a/2);b0=round(b/2);d=10;for i=

2、1:afor j=1:b distance=sqrt(i-a0)2+(j-b0)2); if distance0.2)=0;figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong);% 畫三維曲面（色）圖1.1.2理想低通濾波器實踐結(jié)果截圖：1.2.1 Butterworth低通濾波器實踐代碼：clear all;I1=imread(girl.bmp);subplot(221),imshow(I1);xlabel(a)原始圖像);f=double(I1);%數(shù)據(jù)類型轉(zhuǎn)換g=fft2(f);%圖像傅里葉轉(zhuǎn)換g=fftshi

3、ft(g);%傅里葉變換平移F2=log(abs(g);%對傅里葉變換結(jié)果取絕對值，然后取對數(shù)subplot(222),imshow(F2,InitialMagnification,fit);%將計算后的矩陣用圖像表示colormap(jet);%設(shè)置色彩索引圖colorbar %顯示色彩索引條xlabel(b)原始圖像的傅里葉變換圖像);N1,N2=size(g);%傅里葉變換圖像尺寸n=2;%參數(shù)賦初始值d0=5;n1=fix(N1/2);%數(shù)據(jù)圓整n2=fix(N2/2);%數(shù)據(jù)圓整for i=1:N1 %遍歷圖像像素for j=1:N2 d=sqrt(i-n1)2+(j-n2)2);i

4、f d=0h=0;elseh=1/(1+(d/d0)(2*n);endresult(i,j)=h*g(i,j);% 圖像矩陣計算處理endendF3=log(abs(result); %對傅里葉變換結(jié)果取絕對值，然后取對數(shù)subplot(223),imshow(F3,InitialMagnification,fit);colormap(jet);%設(shè)置色彩索引圖colorbar %顯示色彩索引條xlabel(c)濾波后的傅里葉變換圖像)result=ifftshift(result);X2=ifft2(result);X3=uint8(real(X2);subplot(224),imshow(

5、X3)xlabel(d)Butterworth低通濾波圖像);I1=imread(couple.bmp);f1,f2=freqspace(size(I1),meshgrid);D=0.3;r=f1.2+f2.2;n=4;for i=1:size(I1,1) for j=1:size(I1,2) t=r(i,j)/(D*D); Hd(i,j)=1/(tn+1); endendsurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong);% 畫三維曲面（色）圖1.2.2 Butterworth低通濾波器實踐結(jié)果截圖：1.3.1 高斯低通濾波

6、器實踐代碼：IA=imread(girl.bmp);f1,f2=freqspace(size(IA),meshgrid);D=100/size(IA,1);r=f1.2+f2.2;Hd=ones(size(IA);for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=exp(-t); endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);title(原圖像);subplot(2,2,2),imshow(uint8(Ia);title(高斯低通濾波處理);f