已阅读5页,还剩2页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
课程实验报告实验名称:实验3 二维离散Fourier变换实验目的:理解Fourier变换的意义和性质实验设备:微型计算机 编号: 软件 MATLAB 7.x实验原理: 1.选择Lenna图像,用fft2()这个函数对其进行离散Fourier变换。为离散图像的频谱。为幅度频谱(简称幅度谱),为相位谱。而幅度谱的数值范围很大,而显示器的动态范围为8比特,为了清楚地显示幅度谱,一般作对数处理,即 Matlab核心代码:A=fft2(I);%Fourier变换Ashft=fftshift(A);%移到中心Am=abs(A);%求幅度谱Am1=0.5*log10(1+Am);%频度谱Am2=0.5*log10(1+abs(Ashft);%把频度谱移到中间Aph=angle(A);%相位谱2.实验步骤:1. 选择Lenna图像,对其进行离散Fourier变换,求得幅度谱和相位谱,并显示幅度谱和相位谱。2. (1)利用完整的相位谱和幅度谱重构Lenna原图像,然后仅利用相位谱重构原图像,最后仅利用幅度谱重构原图像,比较3个实验结果;(2)选择两幅内容不同的图像Lenna和flowers,分别进行Fourier变换,交换二者的相位谱后求Fourier反变换,观察实验结果。3. 分析(1)、(2)的实验结果,说明图像Fourier相位谱的重要性。4. 读取文件名为text.tif的图像,从图中剪切字母“t”,利用相关定理进行模板匹配,找出图中所有t的位置。实验结果及分析1.选择一幅图像,对其进行离散Fourier变换,求得幅度谱和相位谱,并显示幅度谱和相位谱。2.图像的Fourier反变换:(1)利用完整的相位谱和幅度谱重构原图像,然后仅利用相位谱重构原图像,最后仅利用幅度谱重构原图像,比较3个实验结果;分析:(2)选择两幅内容不同的图像,分别进行Fourier变换,交换二者的相位谱后求Fourier反变换,观察实验结果。分析:(3)分析(1)、(2)的实验结果,说明图像Fourier相位谱的重要性。分析:4.读取文件名为text.tif的图像,从图中剪切字母“t”,利用相关定理进行模板匹配,找出图中所有t的位置。实验中的体会(如实验过程中遇到的问题及其解决的方法等)本次实验内容比较多一点,实验课堂上也完全得差不多,发现实验课上我做的速度还是挺快的,但是有时候总是去纠结一个结果。就像上次实验课,我检查了好几遍,发现我的代码是没有问题,也请老师和同学帮我看了,但是出不来效果,一直停留在这个问题,没心思做最后一题,可是回到宿舍,自己又敲了一次代码,结果一次就通过,也很顺利地完全了最后一题。这是第三次实验课了,我也比前两次课清晰一点我们究竟在做什么,可能是由于数学基础问题,有些原理还是不能理解得很透彻。虽然实验是知道要干什么,要怎么样做,但是对于背后的原理还是不清楚。不过我还是挺喜欢这个方向,挺喜欢这门课的,我会更加努力地学好这门课的。附录: (MATLAB程序)%实验3% EX1clear;close all;I=imread(E:study图形图像Lenna.tif);I=im2double(I);A=fft2(I);%Fourier变换Ashft=fftshift(A);%移动中心Am=abs(A);%求幅度谱Am1=0.5*log10(1+Am);%频度谱Am2=0.5*log10(1+abs(Ashft);%把频度谱移到中间Aph=angle(A);%相位谱figure;subplot(2,2,1),imshow(I);title(Lenna);subplot(2,2,2),imshow(Am1);title(FT magnitute);subplot(2,2,3),imshow(Aph);title(FT phase);subplot(2,2,4),imshow(Am2);title(FT shift); % EX2(1)I2=ifft2(Am.*exp(i*Aph);%利用幅度谱和相位谱做反变换I3=ifft2(Am);%利用幅度谱做反变换I4=ifft2(255*exp(i*Aph);%仅利用相位谱做反变换figure;subplot(2,2,1),imshow(I);title(Lenna);subplot(2,2,2),imshow(I2);title(Synthesis from the FT);subplot(2,2,3),imshow(I3);title(Magnitute-only Synthesis);subplot(2,2,4),imshow(I4);title(Phase-only Synthesis); % EX2(2)K=imread(E:study图形图像flowers.tif);J=rgb2gray(K);%将图片变成灰白J1=J(10+1:256,222+1:256);J=im2double(J1);B=fft2(J1);Bshft=fftshift(B);Bm=abs(B);Bm1=0.5*log10(1+Bm);Bm2=0.5*log10(1+abs(Bshft);Bph=angle(B);J2=ifft2(255.*exp(i*Bph);%利用Lenna的幅度谱和flowers的相位谱J3=ifft2(255.*exp(i*Aph);%利用flowers的幅度谱和Lenna的相位谱figure;subplot(2,2,1),imshow(I);title(Lenna);subplot(2,2,2),imshow(J1);title(flowers);subplot(2,2,3),imshow(J2);title(Lenna的幅度谱和flowers的相位谱);subplot(2,2,4),imshow(J3);title(flowers的幅度谱和Lenna的相位谱); % EX3I=imread(E:study图形图像text.tif);t=imcrop(I,145 30 6 20);t1=real(ifft2(fft2(I).*rot90(fft2(t,256,256),2);t2=max(t1(:);thresh=t2*0.9;
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024-2030年中国齿形皮带轮行业经营发展分析及未来发展战略规划报告
- 2024-2030年中国黄金饰品行业消费规模预测及前景经营效益分析报告
- 2024-2030年中国麦芽糖浆行业市场竞争格局及发展趋势与投资前景研究报告
- 2024-2030年中国高频头行业发展现状及投资策略分析研究报告
- 2024-2030年中国高纯度乙硼烷市场深度调查与发展前景预测分析报告
- 2024-2030年中国高端路由器行业发展潜力评估及市场前景预判报告
- 2024-2030年中国高熔点费托蜡行业供需态势及前景动态预测报告
- 2024-2030年中国高压清洗机行业经营策略及应用前景预测报告
- 2024-2030年中国骨移植和替代物行业市场发展趋势与前景展望战略研究报告
- 2024-2030年中国马铃薯行业供需趋势及投资风险研究报告
- 少年宫微机组活动记录表
- 《秋兴八首(其一)》课件34张
- 《无人机组装与调试》课件 第一章
- 水炮施工方案 消防水炮安装施工方案
- STA无痛麻醉仪课件
- 2023年南方出版社六年级信息技术下册教案
- 博腾变频器说明书
- 医院感染暴发应急预案演练脚本与医院感染暴发应急预案演练记录
- 建筑工地防台风安全安全措施专项检查表
- 阳光房用料报价单
- 鼠疫演练脚本
评论
0/150
提交评论