二维灰度图象的统计分析及FFT变换处理
《基础强化训练》设计报告题目：二维灰度图象的统计分析及FFT变换处理专业班级：学生姓名：指导教师：武汉理工大学信息工程学院2010年7月17日《基础强化训练》设计任务书学生姓名：专业班级：指导教师：工作单位：信息工程学院题目: 二维灰度图象的统计分析及FFT变换处理课程设计目的：①较全面了解常用的数据分析与处理原理及方法；②能够运用相关软件进行模拟分析；③掌握基本的文献检索和文献阅读的方法；④提高正确地撰写论文的基本能力。
课程设计内容和要求①采集一幅像素大于64*64黑白图像；②常规的数学统计数据处理：计算图象各象素点灰度值得均值、标准差、方差，并绘出灰度直方图；③采用[FFT（傅立叶变换）]对图像进行分析初始条件：①MATLAB软件。
②数字信号处理与图像处理基础知识。
时间安排：第18周周一：安排任务19～20周：仿真设计（鉴主13楼计算机实验室）第20周周六：完成（答辩，提交报告，演示）指导教师签名：年月日系主任（或责任教师）签名：年月日目录摘要........................................... 错误!未定义书签。
Abstract........................................ 错误!未定义书签。
1 数据采集...................................... 错误!未定义书签。
图像的选取.................................. 错误!未定义书签。
MATLAB读取方法............................. 错误!未定义书签。
编辑M文件.............................. 错误!未定义书签。
图像的读取.............................. 错误!未定义书签。
查看图像的格式.......................... 错误!未定义书签。
灰度值的获取............................ 错误!未定义书签。
2 数据统计处理.................................. 错误!未定义书签。
均值计算.................................... 错误!未定义书签。
原理介绍................................ 错误!未定义书签。
仿真结果................................ 错误!未定义书签。
标准差计算.................................. 错误!未定义书签。
原理介绍................................ 错误!未定义书签。
仿真结果................................ 错误!未定义书签。
方差计算.................................... 错误!未定义书签。
原理介绍................................ 错误!未定义书签。
仿真结果................................ 错误!未定义书签。
绘制灰度直方图.............................. 错误!未定义书签。
直方图定义.............................. 错误!未定义书签。
直方图的绘制............................ 错误!未定义书签。
3 快速傅立叶变换（FFT）......................... 错误!未定义书签。
原理介绍................................ 错误!未定义书签。
仿真结果.................................... 错误!未定义书签。
4 总结和心得.................................... 错误!未定义书签。
参考文献........................................ 错误!未定义书签。
摘要数字图像处理是一门新兴的综合性学科，它在交通、通信、工业、军事、医学、遥感等多领域被广泛应用。
伴随着计算机技术的发展，图像处理的应用领域还在不断扩大。
图像处理工程可分为图像处理、图像分析和图像理解三个层次,图像处理是比较低层的操作，它主要在图像的像素级上进行处理，而且处理的数据非常庞大，需要用到大量的工程数学知识。
Matlab集数值计算、符号运算及图形处理等强大功能于一体，是当今国际上公认的最优秀的科技应用软件之一。
它编写简单，具有强大的科学计算能力、可视化功能和开放式可扩展环境，因此在图像处理领域得到了广泛的应用。
本次练习就是Matlab环境下的一些最基本的图像处理操作，如图像的读取、显示、直方图统计、图像增强、傅里叶变换等处理，为将来迅速进入图像处理应用领域打下基础。
AbstractDigital image processing is a new comprehensive subject which has been widely used in many fields like transportation, communications, industrial, military, medical, remote the development of computer technology, image processing applications is still processing projects can be divided into three levels,image processing, image analysis and image processing operations are considered relatively low-level. It is mainly used on the pixel level and it handles very large it needs to use a lot of engineering mathematics.Matlab is good at numerical calculations, symbolic manipulation and graphics processing. It is internationally recognized as the best application of technology can be used for programming with its powerful scientific computing, visualization and open extensible it has been widely used in image processing.We will do the most basic image processing operations under the Matlab environment, such as image reading, displaying,histogram statistics,image enhancement and Fourier transform practice can help us prepare for do image processing quickly in the future.1 数据采集Matlab软件的图像采集工具箱提供了大量的函数用于采集图像和视频信号，该工具箱支持的硬件设备包括工业标准的PC图像采集卡和相应的设备。
所支持的硬件包括：Matrox和Data Translation公司提供的视频采集设备，同时还支持Windows平台下支持USB视频摄像头等设备。
Matlab可以支持很多很多图片的文件格式，例如﹡cur、﹡bmp、﹡hdf、﹡ico、﹡jpg、﹡tif、﹡gif、﹡pcx等。
用于读写图像数据的函数有imread（从图像格式的文件中读取数据）和imwrite（将数据写入图像格式的文件）。
这里我们主要是用到图像的读取。
图像的选取目前我们刚接触图像的处理，为了便于学习和降低难度，根据训练要求，我们采集一幅像素略大于64*64黑白图像如下：图选取的图片MATLAB读取方法编辑M文件打开Matlab软件，在File下选择New\M-File便可以编辑M文件。
M文件分为脚本式和函数式，我们可以利用M文件实现自己的算法即可以建立自己的函数工具箱。
例如读取图像我们可以编辑一个的文件来实现读取图像的功能。
在以后的处理过程中我们几乎每部分都要建立M文件，到时将具体列出其内容。
图像的读取imread函数用于读入各种图象文件，其一般的用法为：利用函数imread()可完成图形图像文件的读取，语法：A=imread(filename,fmt)fmt为图象的格式，filename为读取的图象文件（可以加上文件的路径）。
通常，读取的大多数图像均为8bit，当这些图像加载到内存中时，Matlab就将其存放在类uint8中。
例如本次练习中： I = imread（‘’）;利用函数imread函数实现图像文件的读取，然后运用imshow函数显示图像。
具体操作如下：将所要处理的图片存在D:\matlab\work里，若要将图片“”写入Matlab中。
需要先打开运行Matlab，然后在File里面新建M-File文件，再将读取图片M文件的程序如下：I= imread ('');figureimshow(I);title('西院的老树');在Matlab的command windows栏里运行read图读取后显示图片查看图像的格式图形图像文件信息的查询 imfinfo( ) 函数imfinfo( '')在Matlab的command windows栏里运行该函数，可查看该图像的格式>> imfinfo( '')ans =Filename: ''FileModDate: '25-May-2010 22:17:58'FileSize: 58954Format: 'jpg'FormatVersion: ''Width: 353Height: 250BitDepth: 8ColorType: 'grayscale'FormatSignature: ''NumberOfSamples: 1CodingMethod: 'Huffman'CodingProcess: 'Sequential'Comment: {'ACD Systems Digital Imaging'}灰度值的获取从上述文件的信息可以看到该文件的颜色类型ColorType: 'grayscale'，所以我们不需转换格式便可以查看该图片的灰度值。
图像的二维傅里叶变换
图像傅立叶变换（二维傅立叶变换fourier, 二维DFT, 2d-fft）的原理和物理意义图像傅立叶变换图像的傅立叶变换，原始图像由N行N列构成，N必须是基2的，把这个N*N个包含图像的点称为实部，另外还需要N*N个点称为虚部，因为FFT是基于复数的，如下图所示：计算图像傅立叶变换的过程很简单：首先对每一行做一维FFT，然后对每一列做一维FFT。
具体来说，先对第0行的N个点做FFT（实部有值，虚部为0），将FFT输出的实部放回原来第0行的实部，FFT输出的虚部放回第0行的虚部，这样计算完全部行之后，图像的实部和虚部包含的是中间数据，然后用相同的办法进行列方向上的FFT变换，这样N*N的图像经过FFT得到一个N*N的频谱。
下面展示了一副图像的二维FFT变换：频域中可以包含负值，图像中灰色表示0，黑色表示负值，白色表示正值。
可以看到4个角上的黑色更黑，白色更白，表示其幅度更大，其实4个角上的系数表示的是图像的低频组成部分，而中心则是图像的高频组成部分。
除此以外，FFT的系数显得杂乱无章，基本看不出什么。
将上述直角坐标转换为极坐标的形式，稍微比较容易理解一点，幅度中4个角上白色的区域表示幅度较大，而相位中高频和低频基本看不出什么区别来。
上述以一种不同的方法展示了图像频谱，它将低频部分平移到了频谱的中心。
这个其实很好理解，因为经2D-FFT的信号是离散图像，其2D-FFT的输出就是周期信号，也就是将前面一张图周期性平铺，取了一张以低频为中心的图。
将原点放在中心有很多好处，比如更加直观更符合周期性的原理，但在这节中还是以未平移之前的图来解释。
行N/2和列N/2将频域分成四块。
对实部和幅度来说，右上角和左下角成镜像关系，左上角和右下角也是镜像关系；对虚部和相位来说，也是类似的，只是符号要取反，这种对称性和1维傅立叶变换是类似的，你可以往前看看。
为简单起见，先考虑4*4的像素，右边是其灰度值，对这些灰度值进行2维fft变换。
维灰度图象的统计分析DWT变换处理matlab图像处理
《基础强化训练》设计报告题目:二维灰度图象的统计分析及DWT变换处理专业班级：学生姓名：指导教师：信息工程学院2010 年 7 月日《基础强化训练》设计任务书学生姓名：专业班级：指导教师：工作单位：题目: 二维灰度图象的统计分析及（FFT，DCT，DWT）变换处理课程设计目的：①较全面了解常用的数据分析与处理原理及方法；②能够运用相关软件进行模拟分析；③掌握基本的文献检索和文献阅读的方法；④提高正确地撰写论文的基本能力。
课程设计内容和要求①采集一幅像素大于64*64黑白图像；②常规的数学统计数据处理：计算图象各象素点灰度值得均值、标准差、方差，并绘出灰度直方图；③采用[FFT（傅立叶变换），DCT（离散余弦变化），DWT（离散小波变化）]对图像进行分析初始条件：①MATLAB软件。
②数字信号处理与图像处理基础知识。
时间安排：第18周周一：安排任务19～20周：仿真设计（鉴主13楼计算机实验室）第20周周六：完成（答辩，提交报告，演示）指导教师签名：年月日系主任（或责任教师）签名：年月日目录摘要 (1)Abstract (2)1 数据采集 (3)1.1图像的选取 (3)1.2 MA TLAB编辑 (3)1.2.1 编辑M文件 (3)1.2.2 图像的读取与显示 (3)2 数据统计处理 (6)2.1 均值计算 (6)2.1.1 原理介绍 (6)2.1.2 仿真结果 (6)2.2 标准差计算 (7)2.3 方差计算 (8)图2.5方差的计算程序 (8)图2.6 方差的结果显示 (9)2.4 灰度直方图 (9)3 离散小波变换 (11)3.1 原理介绍 (11)3.2 仿真结果 (12)图3.2结果显示 (13)4参考文 (14)5心得与体会 (15)摘要MATLAB是集数值计算，符号运算及图形处理等强大功能于一体的科学计算语言。
作为强大的科学计算平台，它几乎能够满足所有的计算需求。
MATLAB 本身就是功能强大的数据可视化工具，可以通过各种形式显示分析数据，例如灰度直方图、等高线、蒙太奇混合、像素分析、图层变换以及材质贴图等。
基于MATLAB的二维灰度图象的统计分析及FFT变换处理
目录摘要 (1)Abstract (3)1 数据采集 (4)1.1 选取图像 (4)1.2 MATLAB读取方法 (5)1.2.1 创建M文件 (5)1.2.2 读取图像 (5)1.2.3 查看图像格式 (6)1.2.4 灰度值的获取 (8)2 数据统计处理 (9)2.1各像素点灰度值的均值 (9)2.1.1算法原理及公式 (9)2.1.2 MATLAB源程序及运行结果 (9)2.2各像素点灰度值的方差 (11)2.2.1算法原理及公式 (11)2.2.2 MATLAB源程序及运行结果 (11)2.3各像素点灰度值的标准差 (12)2.3.1算法原理及公式 (12)2.3.2 MATLAB源程序及运行结果 (12)2.4灰度直方图 (12)3 傅里叶变换 (14)3.1算法原理及公式 (14)3.2 FFT及反FFT源程序及运行结果 (14)3.3对变换后的图像进行说明 (17)4 小结与体会 (18)5参考文献 (19)摘要MATLAB(Matrix Laboratory)是MathWorks公司于1989年推出的一套高性能的数值计算和可视化软件。
它集数值分析、矩阵运算、信号处理和图形于一体，构成了一个方便且界面友好的用户环境，可进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
这次基础强化训练就是运用MATLAB做基本的图像处理与数据分析，掌握基本的图像处理方法。
关键词：MATLAB，图像处理，数据分析AbstractMATLAB (Matrix Laboratory) is the MathWorks Inc. in 1989 on a high-performance numerical computation and visualization software.It combines numerical analysis, matrix computation, signal processing and graphics in one, form a convenient and user-friendly user environment, can be matrix operations Huizhi functions and data, implementation algorithms, create user interfaces, connect to other programming language procedures, mainly used in engineering calculations, control design, signal processing and communications, image processing, signal detection, financial modeling design and analysis and other fields.The basis of intensive training is to use MATLAB to do basic image processing and data analysis, a grasp of basic image processing method.Key words: MATLAB, Image Processing, Data Analysis1 数据采集1.1 选取图像从互联网上下载一幅图片，像素为874*590，符合要求。
二维傅里叶变换分离实现
二维傅里叶变换分离实现二维傅里叶变换（2D Fourier Transform）是信号处理中重要的工具之一，它能够将时域信号转换为频域信号，从而方便进行一系列的分析和处理。
在信号处理领域，这种技术被广泛应用于图像处理、声音处理、通信等各种领域。
要实现二维傅里叶变换的分离，需要以下几个步骤：1.读取图像：首先需要读取待处理的图像，并将其转换为灰度图像。
灰度图像是一个二维函数的离散表示，其中像素的灰度值代表了该像素点的亮度信息。
2.傅里叶变换：将灰度图像进行傅里叶变换，可以得到图像中各个频率的特征分量。
傅里叶变换可以通过快速傅里叶变换（FFT）算法来实现，该算法能够高效地计算离散傅里叶变换。
3.频域分析：对傅里叶变换得到的频域图像进行分析，可以得到图像中各个频率的特征分量。
这些特征分量可以表示为幅度和相位信息，幅度表示该频率的强度，相位表示该频率的相对位置。
4.分离处理：根据需要分离和处理特定的频率成分。
可以选择保留几个频率成分，或者将几个频率成分剔除，或者对几个频率成分进行增强等处理。
5.逆傅里叶变换：将处理后的频域图像进行逆傅里叶变换，得到经过处理后的图像。
逆傅里叶变换也可以通过FFT算法来实现。
下面给出一个简单的Python代码，实现了二维傅里叶变换的分离处理。
该代码使用了Python的NumPy库和OpenCV库。
```pythonimport cv2import numpy as np#读取图像并转换为灰度图像image = cv2.imread("image.jpg")gray_image = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY) #傅里叶变换（快速傅里叶变换）fft_image = np.fft.fftshift(np.fft.fft2(gray_image)) #频域分析magnitude_spectrum = 20 * np.log(np.abs(fft_image))#分离处理（保留部分频率成分）rows, cols = gray_image.shapecrow, ccol = int(rows / 2), int(cols / 2)mask = np.zeros((rows, cols), np.uint8)mask[crow-30:crow+30, ccol-30:ccol+30] = 1fft_image *= mask#逆傅里叶变换（快速逆傅里叶变换）ifft_image = np.fft.ifft2(np.fft.ifftshift(fft_image)) ifft_image = np.abs(ifft_image) # 取绝对值#显示结果cv2.imshow("Original Image", gray_image)cv2.imshow("Magnitude Spectrum", magnitude_spectrum)cv2.imshow("Processed Image", ifft_image.astype(np.uint8))cv2.waitKey(0)cv2.destroyAllWindows```在这个例子中，我们首先读取并转换图像为灰度图像。
二维灰度图象的统计分析 DWT变换处理 matlab图像处理word精品文档9页
《基础强化训练》设计报告题目:二维灰度图象的统计分析及DWT变换处理专业班级：学生姓名：指导教师：信息工程学院2019 年7 月日《基础强化训练》设计任务书学生姓名：专业班级：指导教师：工作单位：题目: 二维灰度图象的统计分析及（FFT，DCT，DWT）变换处理课程设计目的：①较全面了解常用的数据分析与处理原理及方法；②能够运用相关软件进行模拟分析；③掌握基本的文献检索和文献阅读的方法；④提高正确地撰写论文的基本能力。
课程设计内容和要求①采集一幅像素大于64*64黑白图像；②常规的数学统计数据处理：计算图象各象素点灰度值得均值、标准差、方差，并绘出灰度直方图；③采用[FFT（傅立叶变换），DCT（离散余弦变化），DWT（离散小波变化）]对图像进行分析初始条件：①MATLAB软件。
②数字信号处理与图像处理基础知识。
时间安排：第18周周一：安排任务19～20周：仿真设计（鉴主13楼计算机实验室）第20周周六：完成（答辩，提交报告，演示）指导教师签名：年月日系主任（或责任教师）签名：年月日目录摘要 (1)Abstract (1)1 数据采集 (1)1.1图像的选取 (2)1.2 MATLAB编辑 (2)1.2.1 编辑M文件 (2)1.2.2 图像的读取与显示 (2)2 数据统计处理 (3)2.1 均值计算 (3)2.1.1 原理介绍 (3)2.1.2 仿真结果 (3)2.2 标准差计算 (3)2.3 方差计算 (4)图2.5方差的计算程序 (4)图2.6 方差的结果显示 (4)2.4 灰度直方图 (4)3 离散小波变换 (4)3.1 原理介绍 (4)3.2 仿真结果 (5)图3.2结果显示 (5)4参考文 (5)5心得与体会 (5)摘要MATLAB是集数值计算，符号运算及图形处理等强大功能于一体的科学计算语言。
作为强大的科学计算平台，它几乎能够满足所有的计算需求。
MATLAB 本身就是功能强大的数据可视化工具，可以通过各种形式显示分析数据，例如灰度直方图、等高线、蒙太奇混合、像素分析、图层变换以及材质贴图等。
二维图像傅里叶变换
a 图
a 图的相位谱重 值谱结合，进行傅里叶反变换
a 图
b 图的幅值谱重构图
由此更加说明相 位谱较幅值谱更能 影响图像的轮廓。
a 图 的 幅 值 谱
b 图 的 相 位 谱
b 图的大体轮廓
b图的幅值谱与 a图的相位谱组合
a 图 的 相 位 谱
b 图 的 幅 值 谱
a 图的大体轮廓
由此可以说明相 位谱较幅值谱更能影 响图像的形状。通俗 的说，幅度决定图像 的强弱，相位决定图 像的频率。
先将幅值谱设为常数（这里 设为1），然后和图像原来的相 位谱结合，进行傅里叶反变换
图像的傅里叶 变换是图像在空域和 频域之间的变换
. 二维傅 . 里叶变换可 以由两步一 维傅里叶变 换来实现
幅度和频率哪个
更能影响图像的形状呢 请看如下试验
先 准 备 两 张 图 片
a 图
b 图
a 图 的 幅 值 谱
b 图 的 幅 值 谱
b 图的相位谱
a 图的相位谱
图a的幅值谱 和图b的相位谱重 新组合
二维灰度图象的统计分析及变换处理
《基础强化训练》报告题目：二维灰度图象的统计分析及变换处理专业班级：通信0701学生姓名：指导教师：刘新华武汉理工大学信息工程学院年月日课程设计任务书学生姓名：专业班级：通信0701指导教师：刘新华工作单位：信息工程学院题目: 二维灰度图象的统计分析及变换处理课程设计目的：①较全面了解常用的数据分析与处理原理及方法；②能够运用相关软件进行模拟分析；③掌握基本的文献检索和文献阅读的方法；④提高正确地撰写论文的基本能力。
课程设计内容和要求①采集一幅像素大于64*64灰度图像；②常规的数学统计数据处理：计算图象各象素点灰度值得均值、标准差、方差，并绘出灰度直方图；③采用[FFT（傅立叶变换），DCT（离散余弦变化），DWT（离散小波变化）]对图像进行分析初始条件：①MATLAB软件。
②数字信号处理与图像处理基础知识。
时间安排：第20周，安排任务（鉴3-206，6月30日）第21周，程序设计与计算（鉴主13楼计算机实验室）第22周，完成（答辩，提交报告，演示）指导教师签名：年月日系主任（或责任教师）签名：年月日目录摘要 (I)Abstract (II)1数据采集 (1)1.1图像的选取 (1)1.2图像的读取 (1)1.2.1编辑M文件 (1)1.2.2图像的读取 (1)1.2.3查看图片格式 (2)1.2.3.1判断图像真彩或是灰度图像 (2)1.2.3.2转变成灰度图像 (3)1.2.3.3判断并查看灰度图像的格式 (3)1.2.4灰度值的获取 (4)2数据统计处理 (5)2.1各像素点灰度值均值 (5)2.1.1原理及相关函数 (5)2.1.2程序执行和结果 (5)2.2像素点灰度值标准差 (6)2.3像素点灰度值方差 (6)2.4绘制灰度直方图 (6)2.4.1灰度直方图的定义 (6)2.4.2灰度直方图的绘制 (6)3傅里叶变换（FFT） (8)3.1二维离散傅里叶变换的定义 (8)3.2傅里叶变换的相关函数 (8)3.3灰度图像的傅里叶变换及其逆变换 (9)3.4傅里叶变换谱的图像分析和变换的意义 (10)4离散余弦变换(DCT) (11)4.1二维离散余弦变换的定义和应用 (11)4.2灰度图像的二维离散余弦变换 (11)5离散小波变换(DWT) (13)5.1离散小波变换的应用 (13)5.2灰度图像的离散小波变换 (13)6心得体会 (15)7参考文献 (16)摘要随着人类社会和科技的发展，人们对信息处理和交流的要求越来越高。
实验一 图像的二维离散傅立叶变换
实验一图像的二维离散傅立叶变换一、实验目的掌握图像的二维离散傅立叶变换以及性质二、实验要求1) 建立输入图像,在6464的黑色图像矩阵的中心建立1616的白色矩形图像点阵,形成图像文件。
对输入图像进行二维傅立叶变换,将原始图像及变换图像(三维、中心化)都显示于屏幕上。
2) 调整输入图像中白色矩形的位置,再进行变换,将原始图像及变换图像(三维、中心化)都显示于屏幕上,比较变换结果。
3) 调整输入图像中白色矩形的尺寸(4040,44),再进行变换,将原始图像及变换图像(三维、中心化)都显示于屏幕上,比较变换结果。
三、 实验仪器设备及软件HP D538、MATLAB四、 实验原理在二维情况下,定义 f(x,y)的傅立叶变换Ｆ(u,v) :2()2()(,)(,)(,)(,)j ux vy j ux vy F u v f x y e dxdy f x y F u v e dudv ππ∞∞-+-∞-∞∞∞+-∞-∞==⎰⎰⎰⎰它表明了空间频率成分与二维图像信号之间的相互关系对于我们要处理的实际二维图像,其傅氏变换一般就是在频率域上有界的,亦即有用成分总就是落在一定的频率域范围之内上述的频率域性质的依据在于:一就是图像中景物的复杂性具有一定的限度,其中大部分内容就是变化不大的区域完全像“雪花”点似的图像没有任何实际意义。
二就是人眼对空间复杂性(频率)的分辨率以及显示器的分辨能力都就是具有一定限度。
若实变量函数f(x)就是绝对可积的,即:且F(u)就是可积的,则傅立叶变换对一定存在。
(){}()()[](){}()()[]du ux j u F x f u F dx ux j x f u F ππ2exp 2exp ⎰⎰∞∞-∞∞-==-==1-F xf F如果f(x)考虑为实函数,它的傅立叶变换通常就是复数形式,即:()()()u jI u R u F +=也可表为:()()()u j e u F u F φ=若二变量函数f(x,y) 就是绝对可积的,即:且F(u,v)就是可积的,则傅立叶变换对一定存在。
对灰度图像作傅立叶变换并显示其频谱
实验项目 实验二、 对灰度图像作傅 指导教师 名 称 立叶变换并显示其频谱
二、 实验原理 离散傅立叶变换建立了离散空域和离散频域之间的联系，是图像处理的重要工具，广泛应用 于图像增强、复原、编码等方面 二维 DFT 的变换式及逆变换式为：
F (u , v ) 1 MN f ( x, y ) 1 MN
三、 程序流程
按要求生成所需 图像
对上述图像进行二 维离散傅立叶变换
将 0 频分量移到 图像中心位置
按公式 lg(1+|F(u,v)|)显示其 傅立叶频谱
四、程序代码 f=zeros(30,30); f(5:24,13:17)=1; F=fft2(f,256,256); F2=fftshift(F); F3=log(1+abs(F2)); subplot(1,2,1); imshow(f); title('原图象'); subplot(1,2,2); imshow(F3); title('傅立叶幅度谱');
M 1 N 1
j 2 (
ux M
vy N
)
f ( x , y )e
j 2 (
其中：u, x=0, 1, 2, …, M-1； v, y=0, 1, 2, …, N-1；
x0 y0 ux M vy N )
ቤተ መጻሕፍቲ ባይዱ
M 1 N 1
F (u , v )e
u 0 v0
MATLAB 中有函数 fft2 可直接对图像进行离散傅立叶变换 变换后可用 fftshit 将 0 频分量平移到图像的中心位置
实验报告
开课实验室： 物理学院 学院 课程 名称 教 师 评 语 一、 实验目的 掌握用 MATLAB 对图像进行傅立叶变换并显示其频谱的方法; 加深对图像傅立叶变换的理解 教师签名： 年 月 日 物理学院 年级、专业、班 09 电子 01 班 姓名 2012 年 3 月 7 日 成绩
图像处理与傅里叶变换原理与运用
图像处理与傅里叶变换1背景傅里叶变换是一个非常复杂的理论，我们在图像处理中集中关注于其傅里叶离散变换离散傅立叶变换(Discrete Fourier Transform) 。
1.1离散傅立叶变换图象是由灰度（RGB ）组成的二维离散数据矩阵，则对它进行傅立叶变换是离散的傅立叶变换。
对图像数据f(x,y)（x=0,1,… ,Ｍ-1； y=0,1,… ,N-1）。
则其离散傅立叶变换定义可表示为：式中，u=0,1,…, M-1;v= 0,1,…, N-1 其逆变换为式中，x=0,1,…, M-1;y= 0,1,…, N-1在图象处理中，一般总是选择方形数据，即Ｍ＝Ｎ影像f(x,y)的振幅谱或傅立叶频谱： 相位谱：能量谱（功率谱） )1(2exp ),(1),(101∑∑-=-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛+-=M x N y N vy M uxi y x f MNv u F π)2(2exp ),(1),(101∑∑-=-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛+=M u N v N vy M uxi v u F MNy x f π),(),(),(22v u I v u R v u F +=[]),(/),(),(v u R v u I arctg v u =ϕ),(),(),(),(222v u I v u R v u F v u E +==1.2快速傅里叶变化可分离性的优点是二维的傅立叶变换或逆变换由两个连续的一维傅立叶变换变换来实现，对于一个影像f(x,y)，可以先沿着其每一列求一维傅立叶变换，再对其每一行再求一维变换正变化逆变换由于二维的傅立叶变换具有可分离性，故只讨论一维快速傅立叶变换。
正变换 逆变换由于计算机进行运算的时间主要取决于所用的乘法的次数。
按照上式进行一维离散由空间域向频率域傅立叶变换时，对于N 个F∑∑∑∑-=-=-=-=⎥⎦⎤⎢⎣⎡⨯⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+=110101)(2exp ),(1)(2exp ),(1)(2exp ),(1),(N v N u N u N v N vy i v u F NN ux i v u F N N vy ux i v u F NNy x f πππ∑-=⎥⎦⎤⎢⎣⎡-=12exp )(1)(N x N ux i x f Nu F π∑∑∑∑-=-=-=-=⎥⎦⎤⎢⎣⎡-⨯⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡+-=11101)(2exp ),(1)(2exp ),(1)(2exp ),(1),(N y N x N x N y N vy i y x f NN ux i y x f NN vy ux i y x f NNv u F πππ∑-=⎥⎦⎤⎢⎣⎡=12exp )(1)(N u N ux i u F Nx f π（u)值，中的每一个都要进行N 次运算，运算时间与N 2成正比。