第六章 频率图像增强.pptVIP

第六章 频率图像增强 * * 6.1低通滤波法 g(x,y)=h(x,y)*f(x,y) (6.1) G(u,v)=H(u,v)F(u,v) (6.2) g(x,y)=F -1(H(u,v)F(u,v)) (6.3) 理想低通滤波器 (u0,v0)是低频中心，D0是滤波器的截止频率。 H(u,v) ||(u,v)|| O D0 a=imread(c:\1.jpg); b=a(:,:,1); f=fft2(b); ff=f; h=50; xy=imfinfo(c:\1.jpg); x=xy.Width/2; y=xy.Height/2; %(x,y)为图象的中心 ff(x-h:x+h,y-h:y+h)=0; fff=ifft2(ff); imshow(abs(fff),[0,255]); 例1 Butterworth低通滤波器 (u0,v0)是低频中心，D0是滤波器的截止频率。 n阶Butterworth低通滤波器的传递函数由下式定义： 其中 n=2 n=3 H(u,v) D(u,v)/ D0 O 理想低通滤波器在数学上定义得很清楚，在计算机模拟中也可以实现，但不能用实际的电子器件实现，是“非物理”的理想滤波器。Butterworth滤波器是物理上可以实现的。令 B为低通能量百分比。 6.2高通滤波器 g(x,y)=h(x,y)*f(x,y) (6.4) G(u,v)=H(u,v)F(u,v) (6.5) g(x,y)=F -1(H(u,v)F(u,v)) (6.6) 理想高通滤波器 (u0,v0)是低频中心，D0是滤波器的截止频率。 H(u,v) ||(u,v)|| O D0 例2 a=imread(a.jpg); f=dct2(a); h=400; xy=imfinfo(a.jpg); x=xy.Width; y=xy.Height; f(1:y-h,1:x-h-200)=0; imshow(log(1+abs(f)),[0 5],notruesize); %fff=idct2(f); %imshow(5*abs(fff),[0 256],notruesize); Butterworth高通滤波器 (u0,v0)是低频中心，D0是滤波器的截止频率。 n阶Butterworth高通滤波器的传递函数由下式定义： 其中 n=3 H(u,v) D(u,v)/ D0 O 前提：阶数n 相同，D0相等 b=imread(t.gif); %b=(a(:,:,1)+a(:,:,2)+a(:,:,3))/3; colormap(gray); imagesc(b); c=size(b); x=c(1); y=c(2); for i=1:x for j=1:y d(i,j)=double(b(i,j)); end; end; e=fft2(d); f=e; （仅供参考） 例3 n=-4; d00=x/5; d0=d00*d00; for i=1:x for j=1:y f(i,j)=1/(1+(((i-x/2)^2+(j-y/2)^2)/d0)^n); end; end; g=f.*e; h=ifft2(g); imagesc(abs(h),[0,50]); xlabel(n=-4, d0=100);