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

1、實(shí)踐二 : 理想高通濾波器、 Butterworth 高通濾波器、高斯高通濾波器 2.1.1 理想高通濾波器實(shí)踐代碼： I=imread(); subplot(221),imshow(I); title( 原圖像 ); s=fftshift(fft2(I); subplot(223), imshow(abs(s),); title( 圖像傅里葉變換所得頻譜 ); subplot(224), imshow(log(abs(s),); title( 圖像傅里葉變換取對(duì)數(shù)所得頻譜 ); a,b=size(s); a0=round(a/2); b0=round(b/2); d=10; p=;q=; f

2、ori=1:a forj=1:b dista nce=sqrt(i-aOF2+(j-bOF2); ifdistance=dh=0; elseh=1; end; s(i,j)=(p+q*h)*s(i,j); end; end; s=uint8(real(ifft2(ifftshift(s); subplot(222), imshow(s);title( 高通濾波所得圖像 ); I=imread(); f1,f2=freqspace(size(I),meshgrid); Hd=ones(size(I); r=sqrt(f1.A2+f2.A2); Hd(r=O; figure surf(Hd,Fac

3、ecolor,interp,Edgecolor,none,Facelighting,phong);% 維曲面(色)圖 2.1.2 理想高通濾波器實(shí)踐結(jié)果截圖： 2.2.1Butterworth 高通濾波器實(shí)踐代碼： I1=imread(); subplot(121),imshow(I1);畫 title( 原始圖像 ); f=double(I1); g=fft2(f); g=fftshift(g); N1,N2=size(g); n=2; d0=5; n1=fix(N1/2); n2=fix(N2/2); fori=1:N1 forj=1:N2 d=sqrt(i-n1F2+(j-n2)T);

4、ifd=0 h=0; else h=1/(1+(dO/d)A(2* n); end result(i,j)=h*g(i,j); end end result=ifftshift(result); X2=ifft2(result); X3=uint8(real(X2); subplot(122),imshow(X3); title(Butterworth 高通濾波 ); I1=imread(); f1,f2=freqspace(size(I1),meshgrid); D=; r=f1.A2+f2.A2; n=4; fori=1:size(I1,1) forj=1:size(I1,2) t=(D*

5、D)/r(i,j); Hd(i,j)=1/(tAn+1); end end figure surf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong);% 維曲面(色)圖 2.2.2Butterworth 高通濾波器實(shí)踐結(jié)果截圖： 2.3.1 高斯高通濾波器實(shí)踐代碼： clearall IA=imread(); f1,f2=freqspace(size(IA),meshgrid); %D=100/size(IA,1); D=; r=f1.A2+f2.A2; fori=1:size(IA,1) forj=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1-exp(-t); end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; 畫Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(uint8(IA);title( 原始圖像 ); subplot(2,2,2),imshow(uint8(Ia);title( 高斯高通濾波 )