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