数字图像处理实验3冈萨雷斯..docxVIP

实验三图像的傅里叶变换 一、实验目的 1、掌握二维DFT变换及其物理意义2、掌握二维DFT变换的MATLAB程 序3、空间滤波及频域滤波 二、实验内容 1、利用MATLAB实现数字图像的傅里叶变换 读入并显示图Fig0316(3(third_from_top.tif,作该图的二维FFT变换F，将其 直流分量移到频谱中心F1，计算其实部RR、虚部II，用两种方法计算及幅值 A1=abs(F1 和 A2=sqrt(RR.A2+II.A2，分别显示 A1 和 A2，并加以比较。 2、近似冲击函数二维傅里叶变换 A=zeros(99,99; A(49:51,49:51=1; 作A的二维傅里叶变换B，将B直流分量移到频谱中心B1，分别用函数 imshow 和 mesh 显示 A 和 B1 模的对数(log(1+abs(B1) 频域图像 空域图像 1 0.8 10 86 4 2 0.6 PAGE # PAGE # 30 0.4 0.2 0 50 40 30 20 10 10 20 30 40 50 0 50 40 30 20 10 40 50 10 3、空间滤波与频域滤波 将图Fig0504(a(gaussia n-n oise.tif (f )分别进行空间与频域滤波。空间滤波： 用fspecial产生9*9，标准差为2的高斯滤波器w，用函数imfilter对f进行空 间滤波，得到滤 波后的图像fi1。频域滤波 将上述高斯滤波器w用函数freqz2获得其256*256频域形式W，用W在频 域滤波图像f (注意W的直流分量在频谱中心)，得到 Ff ,求其傅里叶逆变换 (ifft2)，得到滤波后的图像fi2。 用函数imshow显示w，用函数 mesh显示 W，用imshow显示f，fi1， Ff的 模的对数，fi2。比 较fi1与fi2。 1 0.8 PAGE # PAGE # 40 0.6 0.4 0.2 0 60 二、实验要求 在屏幕上显示出各 个图像。 程序： i=imread(Fig0316(3(third_from _top.tif; figure,imshow(i; title(Fig0316(3(third_from_top. tif; f=fft2(i; 40 20 20 0 0 f仁 fftshift(f; rr=real(f1; ii=imag(f1; a仁abs(f1; a2=sqrt(rr.A2+iiA2; aa1= log(1+a1; aa2=log(1+a2; figure,imshow(aa1,[]; title(a1 abs; figure,imshow(aa2,[]; title(a2 sqrt; A=zeros(99,99; A(49:51,49:51=1; B=fft2(A; B1= fftshift(B; figure,imshow(A; title( 冲击函数 imsow; figure,mesh(A; title(冲击函数 mesh; figure,imshow(log(1+abs(B1; title(冲击函数频域模 mesh; figure,mesh(log(1+abs(B1 title(冲击函数频域模 mesh; g=imread(Fig0504(a(gaussia n-n oise.tif; G=fft2(g; w=fspecial(gaussia n,9,2; fi1=imfilter(g,w; WHf「eqz2(w256259 W1ll.3shift(w? FlrG*w1? =h2H「ea_(i32(F+ figu「pimshow(w 八 m-e(w 八 figurpmesh(w? m-e(-w 八 figu「pimshow(g 八 m-e(-Fig0504(a(gaussian—noise.Hf(g「figure jmshow(f= 口 K)}「凹?藩二 =hgueimshow(=h2 口 m_e(酋ao藩-