《频域滤波实验指导书》由会员分享，可在线阅读，更多相关《频域滤波实验指导书（3页珍藏版）》请在装配图网上搜索。
1、频域的图像滤波1.实验目的(1)学习MATLAB的相关使用方法。(2) 了解如何给图像添加噪声。(3) 学习如何在频域对图形进行滤波处理。2.实验原理(1) 陷波滤波器指的是一种可以在某一个频率点迅速衰减输入信号，以达到阻碍此频率信号通过的滤波效果。从通过信号的频率范围的角度讲，陷波滤波器属于带阻滤波器的一种，只是他的阻带非常狭窄。既然陷波滤波器属于带阻滤波器,它的阶数必须是二阶(含二阶)以上.(2)使用快速傅立叶变换(fft)和快速傅立叶逆变换(IFFT)。而fftshift的作用是对fft的输出进行重新排列，将零频分量移到频谱的中心。调用格式：y=fftshift(x)当x为向量时，fft
2、shift(x)直接将x中左右两半交换而产生y。当x为矩阵时，fftshift(x)直接将x中左右、上下进行交换而产生y。3.实验内容(1)对原图像添加已知频率的噪声。(2) 对含噪声的图像进行频域陷波滤波。(3) 原图像有用成分一般位于低频部分，可对图像添加纹理噪声，形成噪声图像；分析纹理的频率，设计陷波滤波器滤除噪声。(4) 实验方法：对沿x轴方向的波纹加性噪声进行陷波滤波。参考程序如下：%1生成波纹噪声图像img=double(rgb2gray(imread(peppers.b)m);pfigure;imshow(img,);sizec=size(img);w=0.4*2*pi;%噪声的
3、数字频率N=2*pi/w;%噪声每一周期的采样点数img_noice=img+20*ones(sizec(l),l)*sin(w*l:sizec(2);figure;imshow(img_noice,);%图像频谱F0=fft2(img);FO=fftshift(FO);figure;imshow(log(abs(FO),);F=fft2(img_noice);F=fftshift(F);figure;imshow(log(abs(F),);%2.设计理想陷波滤波器H=ones(sizec(1),sizec(2);%图像中心点xO=sizec(l)/2+1;yO=sizec(2)/2+l;%噪声所处频率点(x,y)x=x0;y=yO-round(sizec(2)/N);H(x,y-3:y+3)=0;H(x,(yO-y)+yO-3:(yO-y)+yO+3)=O;%3.滤波结果I=ifftshift(F.*H);imgl=ifft2(I);figure;imshow(imgl,);(b)波纹噪声图像(a)原图(c)原图频谱(d)含噪声的图像频谱(e)陷波滤波结果4.思考题仝.廿工口皆AArr卞丄t厶匕；7亠垃3乙(1)对参考程序绐岀功能注释。对实验结果进行分析。中？该实验可以应用到哪些实际问题
频域滤波实验指导书
