图像的频域滤波

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,*,第一题：图像的频域滤波,汇报人：王宪会,第四组：徐德来 、孟凡超、王利辉、王宪会 、王晓光、田伟,1,题目要求,自选图像进行傅里叶变换，在频域分别去除不同大小的高频和低频区域（置零即可），比较空域图像效果 。,2,解题思路,载入图像,傅里叶变换,频域滤波,傅里叶反变换,输出图像,3,解题过程,利用matlab自带的函数FFT2和IFFT2完成实验,FFT2 Two-dimensional discrete Fourier Transform.,FFT2(X) returns the two-dimensional Fourier transform of matrix X. If X is a vector, the result will have the same orientation.,IFFT2 Two-dimensional inverse discrete Fourier transform.,IFFT2(F) returns the two-dimensional inverse Fourier transform of matrix F. If F is a vector, the result will have the same orientation.,4,程序清单,clear;,clc;,f=zeros(512,512);,f(226:286,230:276)=1;,subplot(221);,imshow(f,);,title(原始图像),F=fft2(f); %对图像进行快速傅里叶变换,S=abs(F);,subplot(222),imshow(S,) %显示幅度谱,title(幅度谱（四角）),Fc=fftshift(F); %把频谱坐标原点由左上角移至屏幕中央,5,subplot(223),Fd=abs(Fc);,imshow(Fd,),ratio=max(Fd(:)/min(Fd(:),%ratio = 2.3306e+007,动态范围太大， 显示器无法正常显示,title(幅度谱（中央）),S2=log(1+abs(Fc);,subplot(224),imshow(S2,),title(以对数方式显示频谱),figure(2);,subplot(221);,W=ifft2(F);,imshow(W);,title(原始图像),6,W=F;,W(1:20,1:20)=0;,W=ifft2(W);,subplot(222),imshow(W);,title(左上角置0);,W=F;,W(1:20,492:512)=0;,W=ifft2(W);,subplot(223),imshow(W);,title(右上角置0);,W=F;,W(492:512,1:20)=0;,W=ifft2(W);,subplot(224),imshow(W);,title(左下角置0);,7,实验结果,8,9,实验中遇到的问题和解决方法,%领域平均法是指用某点领域的灰度值平均值来代替该点的灰度值，算法简单，处理速度快，但是在衰减噪声的同时也会使图像产生模糊。,%中值滤波法对噪声的抑制很好，但是不太适合点，线，尖顶等细节较多的图像 。,10,