【转】图像增强(频域低通)

 

图像增强(频域低通)

分类： Matlab学习与实验 2010-12-16 15:44  465人阅读  评论(0)  收藏  举报

fft 照片

1、理想低通滤波

I=imread('001.tif');%将要滤波的图像读到I中，当然要换成你的图片名字了

imshow(I);

[f1,f2]=freqspace(size(I),'meshgrid');

Hd=ones(size(I));

r=sqrt(f1.^2+f2.^2);

Hd(r>0.1)=0;

f=double(I);

Y=fft2(f);

Y=fftshift(Y);

Ya=Y.*Hd;

Ya=ifftshift(Ya);

Ya=ifft2(Ya);

J=uint8(Ya);

figure,imshow(J);

2、Butterworth低通滤波

I=imread('001.tif');%将要滤波的图像读到I中，当然要换成你的图片名字了

imshow(I);

[f1,f2]=freqspace(size(I),'meshgrid');

D=0.4;

n=4;

Hd=ones(size(I));

r=sqrt(f1.^2+f2.^2);

Hd(r>0.1)=0;

for i=1:size(I,1)

for j=1:size(I,2)

t=r(i,j)/(D*D);

Hd(i,j)=1/(t^n+1);

end

end

f=double(I);

Y=fft2(f);

Y=fftshift(Y);

Ya=Y.*Hd;

Ya=ifftshift(Ya);

Ya=ifft2(Ya);

J=uint8(Ya);

figure,imshow(J);

3、指数低通滤波器（高斯滤波器）

I=imread('Snow-Leopard-Prowl.jpg');%将要滤波的图像读到I中，当然要换成你的图片名字了

imshow(I);

[f1,f2]=freqspace(size(I),'meshgrid');

D=0.4;

n=2;

Hd=ones(size(I));

r=sqrt(f1.^2+f2.^2);

Hd(r>0.1)=0;

for i=1:size(I,1)

for j=1:size(I,2)

t=r(i,j)/(D*D);

Hd(i,j)=1/(t^n+1);

end

end

 

f=double(I);

Y=fft2(f);

Y=fftshift(Y);

Ya=Y.*Hd;

Ya=ifftshift(Ya);

Ya=ifft2(Ya);

J=uint8(Ya);

figure,imshow(J);

4、梯形低通滤波器

应用：

1、皮肤美容

2、照片补偿

3、去噪

4、横纹清除