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

  • 126
  • 0
  • 约5.11千字
  • 约 4页
  • 2020-10-08 发布于浙江
  • 举报
图像的二维傅里叶变换和频谱 一、实验目的 通过本实验使学生掌握使用MATLAB 进行二维傅里叶变换的方法,加深对二维傅里叶变换的理解 和图像频谱的理解。 二、实验原理 本实验是基于数字图像处理课程中的二维傅里叶变换理论来设计的。 本实验的准备知识:第四章 频域图像增强中的一维傅里叶变换和二维傅里叶变换,频 域图像增 强的步骤,频域滤波器。 实验用到的基本函数: 一维傅里叶变换函数: fft, 一维傅里叶反变换函数:ifft 频谱搬移函数:fftshift 二维傅里叶变换函数:fft2 二维傅里叶反变换函数:ifft2 绘图函数:imshow, mesh 【说明,如对上述函数的使用方法有疑问,请先用help 命令查询。建议先用help 命令查询器应 用方法,再做具体实验内容。】 例:计算图像 f 的频谱并显示 F=fft2(f); S=abs(F); %求幅度 imshow(S,[]);%显示图像幅度频谱 Fc=fftshift(F); %将图像频谱原点移动到中心显示 imshow(abs(Fc)); 三、实验内容 (一)一维傅里叶变换的实现和分析 1、生成一个一维向量,x=[1 2 3 4 5 6 7 8]; 计算该向量的傅里叶变换,并由傅里叶变换求 反变换,验证结果。 2 在时间域中将x 乘以(-1)n,计算其傅里叶变换,实现傅里叶变换的平移性质使用 fftshift 函 数,实现频谱的平移。 (二)二维傅里叶变换的实现和分析 产生如图所示图象 f1(x,y) (64×64 大小,中间亮条宽 16,高 40,居中,暗处=0,亮处=255), 用 MATLAB 中的 fft2 函数求其傅里叶变换,要求: 1、同屏显示原图f1 和FFT(f1)的幅度谱图; 2、若令 f2(x,y)= (-1)x+y f1(x,y),重复过程 1,比较二者幅度谱的异同,简述理由; 3、若将 f2(x,y)顺时针旋转 90 度得到 f3(x,y),试显示 FFT(f3)的幅度谱,并与 FFT(f2)的幅度谱进行比较。 (三)任意图像的频谱显示任意图像的频谱显示 1、读入图像lenagray.tif,计算该图像的频谱,并将频谱原点移到中心位置显示。 2、读入图像rice.tif,计算该图像的频谱,并将频谱原点移到中心位置显示。 四、实验步骤 (一)一维傅里叶变换的实现和分析 1、程序: x=[1 2 3 4 5 6 7 8]; F=fft(x) F=ifft(F) 运行结果: 1 F = 36.0000 -4.0000 + 9.6569i -4.0000 + 4.0000i -4.0000 + 1.6569i -4.0000 -4.0000 - 1.6569i -4.0000 - 4.0000i -4.0000 - 9.6569i F = 1 2 3 4 5 6 7 8 2 、程序: x=[1 2 3 4 5 6 7 8]; F=fft(x) for i=1:10 y=(-1).^(i-1); end x1=x.*y F1=fft(x1) 运行结果: F = 36.0000 -4.0000 + 9.6569i -4.0000 + 4.0000i -4.0000 + 1.6569i -4.0000 -4.0000 - 1.6569i -4.0000 - 4.0000i -4.0000 - 9.6569i x1 = -1 -2 -3 -4 -5 -6 -7 -8 F1 = -36.0000 4.0000 - 9.6569i 4.0000 - 4.0000i 4.0000 - 1.6569i 4.0000 4.0000 + 1.6569i

文档评论(0)

1亿VIP精品文档

相关文档