图像的二维傅里叶变换和频谱==

1、图像的二维傅里叶变换和频谱一、实验目的通过本实验使学生掌握使用MATLAB进行二维傅里叶变换的方法，加深对二维傅里叶变换的理解和图像频谱的理解。二、实验原理本实验是基于数字图像处理课程中的二维傅里叶变换理论来设计的。本实验的准备知识：第四章频域图像增强中的一维傅里叶变换和二维傅里叶变换，频域图像增强的步骤，频域滤波器。实验用到的基本函数：一维傅里叶变换函数：fft,一维傅里叶反变换函数：ifft频谱搬移函数：fftsh辻t二维傅里叶变换函数：fft2二维傅里叶反变换函数：ifft2绘图函数：imshow,mesh【说明，如对上述函数的使用方法有疑问，请先用help命令查询。建议先用help命令

2、查询器应用方法，再做具体实验内容。】例：计算图像f的频谱并显示F=fft2(f);S=abs(F);%求幅度imshow(S,)：力显示图像幅度频谱Fc二fftsh辻t(F);%将图像频谱原点移动到中心显示imshow(abs(Fc);三、实验内容(一) 一维傅里叶变换的实现和分析1、生成一个一维向量，x=l2345678;计算该向量的傅里叶变换，并由傅里叶变换求反变换，验证结果。2在时间域中将x乘以(-l)n,计算其傅里叶变换，实现傅里叶变换的平移性质使用fftshift函数，实现频谱的平移。(二) 二维傅里叶变换的实现和分析产生如图所示图象fl(x,y)(64X64大小，中间亮条宽16,高

3、40,居中，暗处二0,亮处二255),用MATLAB中的fft2函数求其傅里叶变换，要求：1、同屏显示原图fl和FFT(fl)的幅度谱图；2、若令f2(x,y)=(-1)x+yfl(x,y),重复过程1,比较二者幅度谱的异同，简述理由；3、若将f2(x,y)顺时针旋转90度得到f3(x,y),试显示FFT(f3)的幅度谱,并与FFT(f2)的幅度谱进行比较。(三) 任意图像的频谱显示任意图像的频谱显示1、读入图像1enagray.t辻，计算该图像的频谱，并将频谱原点移到中心位置显示。2、读入图像rice.tif,计算该图像的频谱，并将频谱原点移到中心位置显示。四、实验步骤(一)一维傅里叶变换的

4、实现和分析1、程序：x=l2345678F=fft(x)F=ifit(F)运行结果：CombandVindoColojine1through636.0000-4.0000+9.65691-4.0000+q.0000i70000-*1.6%阪-4.0000-4.0000-4.0Cl00i-4.0000-9.6569iF=123456?82、程序：x=l2345678F=fft(x)fbii=l:10endxl=x.*yFl=fft(xl)运行结果：andovColunnoItHrouek6nnnn-dadoo4d.ascdi.-<i.naooa.nnaat-4.nnaai.eersoi-d

5、.onnnUqX'utuwiT-tHrca><KD-4.DUUU-QODUU3-4QDDU-O.CJ5DS1.-dFl-000000U0000-95691a0000-QOOOOiQO000-16G091q0000Colunnt?LHrouch87400004-aoooo：iaodooos>gijgpi.3、程序:clcx=l2345678F=fft(x)Fa=fftshift(x)Fb=fftsluft(F)运行结果：Co»WindowP=CoLunnsIthroueh636.0000-a.0000f9.65691-4.0000±4.00001-4

6、.0000十1.6569a-4.0000-4.0000-1656加CoLunn?7throueh8-4.0000-1.00001-*3.0000-9.656J>iFb-ColtmrxKIiKrouiK6-4.00004DDOO-1.6560a-4.0000-40OODi-4.ODDO-9.6B60aM6O00(-4-OODO-9.6569aCoLux%rff74hrou(h8-4.000040OODIDDOO*1.(-)1、程序：clccleai-x=zeros(64.64);x(32-20:32+20,32-8:32+8)=255;subplot(1,2,1)±ushow(x

7、);titled原图幅度谱图);F=fft2(x);subplot(l,2,2)±ushow(log(abs(F),);title。傅里叶变换的幅度谱图，);运行结果：x=zeios(64.64);/?:程2x(3220:32+20,32-8:32+8)=255;subplot(2,2,1),iinshow(x);htle(，原图x幅度谱图?；F=fft2(x);subplot(2,2,2),imshow(log(abs(F),);傅里叶变换的幅度谱图?；fori=l:64forj=l:64endendxl=h.*xFl=fft2(xl);subplot(2,2,3),imshow(

8、x1);titleCxl幅度谱图subplot(2,2,4)?imshow(log(abs(Fl),);title('xl傅里叶变换的幅度谱图运行结果：xiwaifiHB受妙的幅啖遷他3、程序：x=zeios(64.64);x(3220:32+20,32-8:32+8)=255;subplot(3,2,l),imshow(x);htle(，原图x幅度谱图?；F=fft2(x);subplot(3,2,2),iinshow(log(abs(F),);title('x傅里叶变换的幅度谱图J;fori=l:64forj=l:64endendxl=h.*xFl=fft2(xl);sub

9、plot(3,2,3),imshow(xl);titleCxl幅度谱图);subplot(3,2,4)，imshow(log(abs(F1),);title('xl幅度谱图，);x2=muotate(xL90);subplot(3,2,5),iinshow(x2);titleCxl旋转90度的x2幅度谱图，)；F2=fft2(x2);subplot(3,2,6),imshow(log(abs(F2),);title('x2傅里叶变换的幅度谱图?；运行结果：原耋>:唱度追團x傅里叶变换的馄度浩耋乂耳冒度ifi耋x2溥里叶变换的Cg度淹耋（三）任意图像的频谱显示程序：i=mu

10、ead(,D:iiiiagelena.bmp')订=subplot(2,3J)iiiishow(i)title('lena原图像*)subplot(2,3,2)imshow(log(abs(i1)，)title('lena频谱图？subplot(2,3,3)mishow(log(abs(i2),)title('lena频谱原点移到中心位置Jill=iiiiiead(,D：iniagerice.png,)subplot(2,3,4)iiiishow(ill)title(Tice原图像?subplot(2,3,5)imshow(log(abs(i22)J)titleCrice频谱图？subplot(2,3,6)i33=fftsliift(i22)imshow(log(abs(i33)?)titleCnce频谱原点移到中心位置?运行结果：leg原田僚38烦暗Elens荻诺床点移到中心燈E：“花预涪国nm频谱加点移到中心位逬五、实验思考题图像频谱有何特点？低频分量和高频分量在图像频谱中是怎样分布的？(1) 频谱图，四个