在频率域中直接生成滤波器

掌握在频率域中直接生成滤波器的方法实验内容1.用 GLPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果2.用 GHPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果3.用 BHPF 对图像高频强调滤波,n=2,观察滤波结果

3-3在频率域中直接生成滤波器

一、 实验目的

掌握在频率域中直接生成滤波器的方法

二、 实验内容

1.用 GLPF 对图像滤波， D0 值为填充图像宽度的 5%,观察滤波结果

2.用 GHPF 对图像滤波， D0 值为填充图像宽度的 5%,观察滤波结果

3.用 BHPF 对图像高频强调滤波，n=2,观察滤波结果

三、 实验步骤

1. 低通滤波器 dftfilt dftuv

1)读入图像

clc

clear

f = imread('lena.jpg');

2）傅立叶频谱图像

F1 = fft2(f); % 注意 F1 和 下面 F 的区别 T=log(1+abs(fftshift(F1)));

T= mat2gray(T);

Imshow(T, []);

title('傅立叶频谱图像')

PQ = paddedsize(size(f));

[U V] = dftuv(PQ(1),PQ(2));

D0 = 0.05*PQ(2);

F = fft2(f,PQ(1),PQ(2));

z=log(1+abs(fftshift(F)));

掌握在频率域中直接生成滤波器的方法实验内容1.用 GLPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果2.用 GHPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果3.用 BHPF 对图像高频强调滤波,n=2,观察滤波结果

z=mat2gray(z);

imshow(z,[]);

title('傅立叶频谱图像')

3)高斯低通滤波器频谱图像:

H = exp(-(U.^2+V.^2)/(2*(D0^2)));

imshow (fftshift(H),[])

title('高斯低通滤波器频谱图像')

g = dftfilt(f,H); %

imshow (g,[])

title('高斯低通处理后图像')

掌握在频率域中直接生成滤波器的方法实验内容1.用 GLPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果2.用 GHPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果3.用 BHPF 对图像高频强调滤波,n=2,观察滤波结果

2. 绘制线框图 mesh lpfilter 低通：

clc

clear

H = fftshift(lpfilter('gaussian',500,500,50));

mesh(H(1:10:500,1:10:500))

axis([0 50 0 50 0 1])

imshow (H,[])

掌握在频率域中直接生成滤波器的方法实验内容1.用 GLPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果2.用 GHPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果3.用 BHPF 对图像高频强调滤波,n=2,观察滤波结果

3. 高通滤波器 hpfilter 高通

clc

clear

H = fftshift(hpfilter('ideal',500,500,100)); % 半径是100

mesh(H(1:10:500,1:10:500))

axis([0 50 0 50 0 1])

colormap([0 0 0])

axis off

grid off

imshow(H,[])

4. 高通滤波

clc clear

掌握在频率域中直接生成滤波器的方法实验内容1.用 GLPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果2.用 GHPF 对图像滤波, D0 值为填充图像宽度的 5%,观察滤波结果3.用 BHPF 对图像高频强调滤波,n=2,观察滤波结果

f = imread('lena.jpg');

PQ = paddedsize(size(f));

D0 = 0.05*PQ(1); % 半径是 D0

H = hpfilter('gaussian',PQ(1),PQ(2),D0);

g = dftfilt(f,H);

imshow (g,[])

title('高斯高通滤波后的图像')

四、 实验总结

通过本次实验，我掌握了在频率域中直接生成滤波器的方法，了解了图像的傅里叶变换，了解了高斯低通滤波、高通滤波的基本原理。

本文来源：https://www.bwwdw.com/article/z921.html