2021年图像的傅里叶变换实验报告.docxVIP

计算机科学与技术系 实 验 报 告 专业名称 计算机科学与技术课程名称 数字图像处理 项目名称 Matlab 语言、 图像傅里叶变换班 级 14 计科 2 班 学 号 姓 名 卢爱胜 同组人员 张佳佳、 王世兜、 张跃文试验日期 .11.30 一、 试验目与要求: （简述此次试验要求达成目, 包含到相关知识点, 试验具体要求。） 试验目: 了解图像变换意义和手段; 熟悉傅立叶变换基础性质; 熟练掌握 FFT 变换方法及应用; 经过试验了解二维频谱分布特点; 经过本试验掌握利用 MATLAB 编程实现数字图像傅立叶变换。 评价人眼对图像幅频特征和相频特征敏感度。 试验要求: 应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析一个有力工具 , 它能够定量地分析诸如数字化系统、 采样点、 电子放大器、 卷积滤波器、 噪音和显示点等作用。经过试验培养这项技能, 将有利于处理大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术人来说, 把时间用在学习和掌握博里叶变换上是很有必需。 二、 试验内容 （依据此次试验项目具体任务和要求, 完成相关内容, 可包含: 试验目、算法原理、 试验仪器、 设备选型及连线图、 算法描述或步骤图、 源代码、 试验运行步骤、 关键技术分析、 测试数据与试验结果、 其她 ） 傅立叶（Fourier）变换定义 对于二维信号, 二维 Fourier 变换定义为: F (u,v) f(x,y)e j2 (ux uy)dxdy 逆变换: f(x,y) F (u,v)ej2 (ux uy)dudv 二维离散傅立叶变换为: F (m ,n) 1 N 1 N 1 f(i,k)e j2 (m i n k ) 逆变换: N N N i 0 k 0 f(i,k) 1 N 1 N 1 F (m ,n)ej2 (m i n k ) N N N m 0 n 0 图像傅立叶变换与一维信号傅立叶变换变换一样 , 有快速算法, 具体参见参考书目, 相关傅立叶变换快速算法程序不难找到。实际上, 现在有实现傅立叶变换芯片, 能够实时实现傅立叶变换。 利用 MATLAB 软件实现数字图像傅立叶变换程序: I=imread(‘原图像名.gif’);%读入原图像文件imshow(I); %显示原图像 fftI=fft2(I); %二维离散傅立叶变换sfftI=fftshift(fftI); %直流分量移到频谱中心RR=real(sfftI); %取傅立叶变换实部 II=imag(sfftI); %取傅立叶变换虚部A=sqrt(RR.^2+II.^2); %计算频谱幅值 A=（A-min(min(A))）/(max(max(A))-min(min(A)))*225 %归一化figure; %设定窗口 imshow(A); %显示原图像频谱 四、 源代码 clc;clear all I=imread(Fig0707(a)(Original).tif.tif); %读入原图像文件imshow(I); %显示原图像 title(原始图像) fftI=fft2(I); %二维离散傅立叶变换sfftI=fftshift(fftI); %直流分量移到频谱中心RR=real(sfftI); %取傅立叶变换实部 II=imag(sfftI); %取傅立叶变换虚部 A=sqrt(RR.^2+II.^2); %计算频谱幅值 A=(A-min(min(A)))/(max(max(A))-min(min(A)))*225; %归一化figure; %设定窗口 imshow(A); %显示原图像频谱title(原始图像频谱) f1=ifft2(A); %用 Fourier 系数幅度进行 Fourier 反变换f2=ifft2(angle(fftI)); %用 Fourier 系数相位进行 Fourier 反变换; figure subplot 121;imshow(f1,[]) title(幅度进行 Fourier 反变换) subplot 122;imshow(f2,[]) title(相位进行 Fourier 反变换) 五、 试验结果及分析 原始图像 原始图像的频谱 幅度进行Fourier反 变换 相位进行Fourier反 变换 试验分析: 此次试验研究了相关傅里叶算法方面知识, 将傅里叶变换应用在图像处理上, 让我学习到了傅里叶算法方面知识, 实践才是成长好道路。 得分（百分制）