数字图像处理实验报告四.pdf

实验四图像频域滤波技术

姓名学号专业年级

刘真志0807030424自动化08级

实验目的：

(1)了解图像增强的目的及意义，加深对图像增强的感性认识，巩固所学的图像增强的理论

知识和相关算法。

(2)熟练掌握频域平滑、锐化滤波的实现方法。

实验重点：频域平滑、锐化滤波的实现方法

实验难点：频域平滑、锐化滤波的实现方法

实验内容：

(1)选择一幅图像，利用高斯低通滤波器对此图像进行平滑滤波，观察不同滤波器的截止半

径对平滑效果的影响。

f=imread(mohu.jpg);

subplot(321);

imshow(f,[]);

title(原图);

F=fftshift(fft2(f,size(f,1),size(f,2)));

u=1:size(F,1);

v=1:size(F,2);

[V,U]=meshgrid(v,u);

D=sqrt((U-(floor(size(F,1)/2)+1)).^2+(V-(floor(size(F,2)/2)+1)).^2);

H=zeros(size(f));

D0=5;

H=exp(-D.^2./(2*D0.^2));

G=F.*H;

G1=ifft2(ifftshift(G));

subplot(322);

imshow(real(G1),[]);

title(半径小于5的高斯低通滤波器);

D0=15;

H=exp(-D.^2./(2*D0.^2));

G=F.*H;

G1=ifft2(ifftshift(G));

subplot(323);

imshow(real(G1),[]);

title(半径小于15的高斯低通滤波器);

D0=30;

H=exp(-D.^2./(2*D0.^2));

G=F.*H;

G1=ifft2(ifftshift(G));

subplot(324);

imshow(real(G1),[]);

title(半径小于30的高斯低通滤波器);

D0=80;

H=exp(-D.^2./(2*D0.^2));

G=F.*H;

G1=ifft2(ifftshift(G));

subplot(325);

imshow(real(G1),[]);

title(半径小于80的高斯低通滤波器);

D0=230;

H=exp(-D.^2./(2*D0.^2));

G=F.*H;

G1=ifft2(ifftshift(G));

subplot(326);

imshow(real(G1),[]);

title(半径小于230的高斯低通滤波器);

(2)选择一幅边缘较模糊的图像，利用高斯高通滤波器对此图像进行增强滤波，以突出图像

中的高频细节信息。

f=imread(mohu.jpg);

subplot(321);

imshow(f,[]);

title(原图);

F=fftshift(fft2(f,size(f,1),size(f,2)));

u=1:size(F,1);

v=1:size(F,2);

[V,U]=meshgrid(v,u);

D=sqrt((U-(floor(size(F,1)/2)+1)).^2+(V-(floor(size(F,2)/2)+1)).^2);

H=zeros(size(f));

D0=5;

H=1-exp(-D.^2./(2*D0.^2));

G=F.*H;

G1=ifft2(ifftshift(G));

subplot(322);

imshow(real(G1),[]);

title(半径小于5的高斯低通滤波器);

D0=15;

H=1-exp(-D.^2./(2*D0.^2));

G=F.*H;

G1=ifft2(ifftshift(G));

subplot(323);

imshow(real(G1),[]);

title(半径小于15的高斯低通滤波器);

D0=30;

H=1-exp(-D.^2./(2*D0.^2));

G=F.*H;

G1=ifft2(ifftshift(G));

subplot(324);

imshow(real(G1),[]);

title(半径小于30的高斯低通滤波器);

D0=80;