matlab图像处理低通滤波器,基于matlab数字图像处理之低通滤波器

基于Matlab数字图像处理实践报告

实践一:理想低通滤波器、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('图像傅里叶变换取对数所得频谱');

[a,b]=size(s);

a0=round(a/2);

b0=round(b/2);

d=10;

for i=1:a

for j=1:b

distance=sqrt((i-a0)^2+(j-b0)^2);

if distance<=d h=1;

else h=0;

end;

s(i,j)=h*s(i,j);

end;

end;

s=uint8(real(ifft2(ifftshift(s))));

subplot(222),

imshow(s);

title('低通滤波所得图像');

I=imread('couple.bmp');

Hd=ones(size(I));

Hd(r>0.2)=0;

figure

surf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong');%画三维曲面(色)图