- 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
您可能关注的文档
最近下载
- (民主生活会)2025年度班子成员相互批评意见清单+对照检查查摆问题清单+批评与自我批评意见建议.docx VIP
- 基于单片机的智能衣柜控制系统设计.docx VIP
- 原神家具负荷表及计算器说明书(多功能小鹏).docx VIP
- 2025 年大学轨道交通信号与控制(信号系统)试题及答案.doc VIP
- 生活中的化学.ppt VIP
- XX村新任支部书记任职发言稿.docx VIP
- 痛苦的缺席——试论苏童的《米》.pdf VIP
- 2024高端装备制造业技术升级改造合同.docx VIP
- 603A使用说明书.doc VIP
- 江苏图集-JS2016Q5 住宅厨卫组合变压式耐火排烟气道系统 江苏省工程建设企业标准设计图集.pdf VIP
原创力文档

文档评论(0)