当前位置:文档之家› (完整版)图像的二维傅里叶变换和频谱==

(完整版)图像的二维傅里叶变换和频谱==

(完整版)图像的二维傅里叶变换和频谱==
(完整版)图像的二维傅里叶变换和频谱==

一、实验目的

通过本实验使学生掌握使用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)

运行结果:

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)

运行结果:

3、程序:

clc

x=[1 2 3 4 5 6 7 8]

F=fft(x)

Fa=fftshift(x)

Fb=fftshift(F)

运行结果:

(二)

1、 程序:

clc

clear

x=zeros(64,64);

x(32-20:32+20,32-8:32+8)=255;

subplot(1,2,1),imshow(x);

title('原图幅度谱图');

F=fft2(x);

subplot(1,2,2),imshow(log(abs(F)),[]);

title('傅里叶变换的幅度谱图');

运行结果:

2、程序:

原图幅度谱

图傅里叶变换的幅度谱图

x=zeros(64,64);

x(32-20:32+20,32-8:32+8)=255;

subplot(2,2,1),imshow(x);

title('原图x 幅度谱图');

F=fft2(x);

subplot(2,2,2),imshow(log(abs(F)),[]);

title('x 傅里叶变换的幅度谱图');

for i=1:64

for j=1:64

h(i,j)=(-1).^(i-1+j-1);

end

end

x1=h.*x

F1=fft2(x1);

subplot(2,2,3),imshow(x1);

title('x1幅度谱图');

subplot(2,2,4),imshow(log(abs(F1)),[]);

title('x1傅里叶变换的幅度谱图');

运行结果:

3、程序:

x=zeros(64,64);

x(32-20:32+20,32-8:32+8)=255;

subplot(3,2,1),imshow(x);

title('原图x 幅度谱图');

F=fft2(x);

subplot(3,2,2),imshow(log(abs(F)),[]);

title('x 傅里叶变换的幅度谱图');

for i=1:64

for j=1:64

h(i,j)=(-1).^(i-1+j-1);

end

end

x1=h.*x

F1=fft2(x1);

subplot(3,2,3),imshow(x1);

原图x 幅度谱

图x 傅里叶变换的幅度谱

x1幅度谱

图x1傅里叶变换的幅度谱图

title('x1幅度谱图');

subplot(3,2,4),imshow(log(abs(F1)),[]);

title('x1幅度谱图');

x2=imrotate(x1,90);

subplot(3,2,5),imshow(x2);

title('x1旋转90度的x2幅度谱图');

F2=fft2(x2);

subplot(3,2,6),imshow(log(abs(F2)),[]);

title('x2傅里叶变换的幅度谱图');

运行结果:

(三)任意图像的频谱显示

程序:

i=imread('D:\image\lena.bmp')

i1=fft2(i)

subplot(2,3,1)

imshow(i)

title('lena 原图像')

subplot(2,3,2)

imshow(log(abs(i1)),[])

title('lena 频谱图')

subplot(2,3,3)

i2=fftshift(i1)

imshow(log(abs(i2)),[])

title('lena 频谱原点移到中心位置')

i11=imread('D:\image\rice.png')

i22=fft2(i11)

subplot(2,3,4)

imshow(i11)

title('rice 原图像')

subplot(2,3,5)

imshow(log(abs(i22)),[])

title('rice 频谱图')

subplot(2,3,6)

x1旋转90度的x2幅度谱

图x2傅里叶变换的幅度谱

原图x 幅度谱

图x 傅里叶变换的幅度谱

x1幅度谱

图x1幅度谱图

i33=fftshift(i22)

imshow(log(abs(i33)),[])

title('rice频谱原点移到中心位置')

运行结果:

五、实验思考题

图像频谱有何特点?低频分量和高频分量在图像频谱中是怎样分布的?

(1)频谱图,四个角对应低频成分,中央部分对应高频成分;图像亮条的平移影响频谱的分布,但当频谱搬移到中心时,图像亮条的平移后频谱图是相同的。图像旋转,频谱也会旋转,并且角度相同。频谱具有平移特性,可分离性。

(2)图像的高低频是对图像各个位置之间强度变化的一种度量方法。低频分量:主要对整副图像的强度的综合度量。高频分量:主要是对图像边缘和轮廓的度量。

图像的傅立叶变换与频域滤波

实验四 图像的傅立叶变换与频域滤波 一、 实验目的 1了解图像变换的意义和手段; 2熟悉傅里叶变换的基本性质; 3熟练掌握FFT 方法的应用; 4通过实验了解二维频谱的分布特点; 5通过本实验掌握利用MATLAB 编程实现数字图像的傅立叶变换。 6、掌握怎样利用傅立叶变换进行频域滤波 7、掌握频域滤波的概念及方法 8、熟练掌握频域空间的各类滤波器 9、利用MATLAB 程序进行频域滤波 二、 实验原理 1应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具,它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能,将有助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说,把时间用在学习和掌握博里叶变换上是很有必要的。 2傅立叶(Fourier )变换的定义 对于二维信号,二维Fourier 变换定义为 : ??∞ ∞ -+-==dxdy e y x f v u F y x f F vy ux j )(2),(),()},({π

二维离散傅立叶变换为: ∑ ∑-=+--==10)(21 01 ),(),(N y N y u M x u j M x MN e y x f v u F π 图像的傅立叶变换与一维信号的傅立叶变换变换一样,有快速算法,具体参见参考书目,有关傅立叶变换的快速算法的程序不难找到。实际上,现在有实现傅立叶变换的芯片,可以实时实现傅立叶变换。 3利用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); %显示原图像的频谱 域滤波分为低通滤波和高通滤波两类,对应的滤波器分别为低通滤波器和 高通滤波器。频域低通过滤的基本思想: G(u,v)=F(u,v)H(u,v) F(u,v)是需要钝化图像的傅立叶变换形式,H(u,v)是选取的一个低通过滤

傅里叶变换图像压缩

傅里叶变换图像压缩

————————————————————————————————作者:————————————————————————————————日期:

DSP实验进度汇报 组员:汪张扬、任艳波、陈雪松、谢聪、沈旭 任务分配:汪张扬由于考G,上周没有任务,沈旭负责自制二值图像的处理,陈雪松和谢聪负责其他图片的处理,任艳波负责搜集图像压缩评价的相关材料 以下为简要概括: 读入图像进行傅里叶变换和压缩 原始程序: a=imread('d:\1.jpg');b=figure;imshow(a);title('原始图像'); F=fft2(a); F_mm=abs(F);figure;imshow(F);title('原始幅度谱'); Fshift=fftshift(F); F_m=abs(Fshift);figure;imshow(F_m);title('幅度谱'); F_p=angle(Fshift);figure;imshow(F_p);title('相位谱'); T=@fft2; B1=blkproc(a,[8 8],T);%将图像分块为8×8矩阵进行处理 figure; imshow(a); title('原始图像'); mask=[100 000 00 0 10 0 0 0 0 0 00 1 000 0 0 00 0 1 000 0 000 0 0000 0 000 0 1 0 0 0 0 0 000 1 0 00 0 0 00 01];%与该矩阵相乘去掉中间行,即高频部分 B2=blkproc(B1,[88],'P1*x',mask); fun=@ifft2; F3=blkproc(B2,[88],fun); F=mat2gray(F3); figure; imshow(F); title('压缩87.5%的图像'); 刚开始的原始图像:

MATLAB数字图像处理几何变换傅里叶变换

Matlab数字图像处理实验指导 实验目的: 通过实验,深入理解和掌握图像处理的基本技术,提高动手实践能力。 实验环境: Matlab变成 实验一图像的几何变换 实验内容:设计一个程序,能够实现图像的各种几何变换。 实验要求:读入图像,打开图像,实现图像的平移变换、比例缩放、转置变换、镜像变换、旋转变换等操作。 实验原理: 图像几何变换又称为图像空间变换,它将一幅图像中的坐标位置映射到另一幅图像中的新坐标位置。学习几何变换的关键就是要确定这种空间映射关系,以及映射过程中的变化参数。 几何变换不改变图像的像素值,只是在图像平面上进行像素的重新安排。一个几何变换需要两部分运算:首先是空间变换所需的运算,如平移、镜像和旋转等,需要用它来表示输出图像与输入图像之间的(像素)映射关系;此外,还需要使用灰度插值算法,因为按照这种变换关系进行计算,输出图像的像素可能被映射到输入图像的非整数坐标上。 设原图像f(x0,y0)经过几何变换产生的目标图像为g(x1,y1),则该空间变换(映射)关系可表示为: x1=s(x0,y0) y1=t(x0,y0) 其中,s(x0,y0)和t(x0,y0)为由f(x0,y0)到g(x1,y1)的坐标换变换函数。 一、图像平移 图像平移就是将图像中所有的点按照指定的平移量水平或者垂直移动。

二、图像镜像 镜像变换又分为水平镜像和垂直镜像。水平镜像即将图像左半部分和右半部分以图像竖直中轴线为中心轴进行对换;而竖直镜像则是将图像上半部分和下半部分以图像水平中轴线为中心轴进行对换。 三、图像转置 图像转置是将图像像素的x坐标和y坐标呼唤。图像的大小会随之改变——高度和宽度将呼唤。

图像傅里叶变换的物理意义

傅里叶变换在图像处理中的作用 图像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度。如:大面积的沙漠在图像中是一片灰度变化缓慢的区域,对应的频率值很低;而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域,对应的频率值较高。傅立叶变换在实际中有非常明显的物理意义,设f是一个能量有限的模拟信号,则其傅立叶变换就表示f的谱。从纯粹的数学意义上看,傅立叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看,傅立叶变换是将图像从空间域转换到频率域,其逆变换是将图像从频率域转换到空间域。换句话说,傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数,傅立叶逆变换是将图像的频率分布函数变换为灰度分布函数 傅立叶变换以前,图像(未压缩的位图)是由对在连续空间(现实空间)上的采样得到一系列点的集合,我们习惯用一个二维矩阵表示空间上各点,则图像可由z=f(x,y)来表示。由于空间是三维的,图像是二维的,因此空间中物体在另一个维度上的关系就由梯度来表示,这样我们可以通过观察图像得知物体在三维空间中的对应关系。为什么要提梯度?因为实际上对图像进行二维傅立叶变换得到频谱图,就是图像梯度的分布图,当然频谱图上的各点与图像上各点并不存在一一对应的关系,即使在不移频的情况下也是没有。傅立叶频谱图上我们看到的明暗不一的亮点,实际上图像上某一点与邻域点差异的强弱,即梯度的大小,也即该点的频率的大小(可以这么理解,图像中的低频部分指低梯度的点,高频部分相反)。一般来讲,梯度大则该点的亮度强,否则该点亮度弱。这样通过观察傅立叶变换后的频谱图,也叫功率图,我们首先就可以看出,图像的能量分布,如果频谱图中暗的点数更多,那么实际图像是比较柔和的(因为各点与邻域差异都不大,梯度相对较小),反之,如果频谱图中亮的点数多,那么实际图像一定是尖锐的,边界分明且边界两边像素差异较大的。对频谱移频到原点以后,可以看出图像的频率分布是以原点为圆心,对称分布的。将频谱移频到圆心除了可以清晰地看出图像频率分布以外,还有一个好处,它可以分离出有周期性规律的干扰信号,比如正弦干扰,一副带有正弦干扰,移频到原点的频谱图上可以看出除了中心以外还存在以某一点为中心,对称分布的亮点集合,这个集合就是干扰噪音产生的,这时可以很直观的通过在该位置放置带阻滤波器消除干扰 注: 1、图像经过二维傅立叶变换后,其变换系数矩阵表明: 若变换矩阵Fn原点设在中心,其频谱能量集中分布在变换系数短阵的中心附近(图中阴影区)。若所用的二维傅立叶变换矩阵Fn的原点设在左上角,那么图像信号能量将集中在系数矩阵的四个角上。这是由二维傅立叶变换本身性质决定的。同时也表明一股图像能量集中低频区域。 2 、变换之后的图像在原点平移之前四角是低频,最亮,平移之后中间部分是低频,最亮,亮度大说明低频的能量大(幅角比较大) 傅立叶变换在图像处理中有非常非常的作用。因为不仅傅立叶分析涉及图像处理的很多方面,傅立叶的改进算法, 比如离散余弦变换,gabor与小波在图像处理中也有重要的分量。 印象中,傅立叶变换在图像处理以下几个话题都有重要作用: 1.图像增强与图像去噪 绝大部分噪音都是图像的高频分量,通过低通滤波器来滤除高频——噪声; 边缘也是图像的高频分量,可以通过添加高频分量来增强原始图像的边缘; 2.图像分割之边缘检测 提取图像高频分量

图像的傅立叶变换和边缘提取

实验五图像的傅立叶变换和边缘提取 兰州大学信息学院 0级通信工程一班赵军伟 第一部分图像的傅立叶变换 一、实验目的 1.了解图像变换的意义和手段; 2. 熟悉傅里叶变换的基本性质; 3. 熟练掌握FFT的方法及应用; 4. 通过实验了解二维频谱的分布特点; 5. 通过本实验掌握利用MATLAB编程实现数字图像的傅立叶变换。 二、实验原理 1.应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具,它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能,将有助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说,把时间用在学习和掌握博里叶变换上是很有必要的。b5E2RGbCAP 2.傅立叶

对于二维信号,二维Fourier变换定义为: 二维离散傅立叶变换为: 三、实验步骤 1.打开计算机,安装和启动MATLAB程序;程序组中“work”文件夹中应有待处理的图像文件; 2.利用MatLab工具箱中的函数编制FFT频谱显示的函数。 3. a>调入、显示三张不同的图像; b>对这三幅图像做FFT并利用自编的函数显示其频谱。 c>讨论不同的图像内容与FFT频谱之间的对应关系。 4.记录和整理实验报告。 四、实验仪器 1计算机, MATLAB软件; 3移动式存储器<软盘、U盘等)。 4记录用的笔、纸。 五、实验结果及程序 1.程序

I1=imread('F:\MATLAB学习\实验\picture\LENA.TIF'>。 %读入原图像文件p1EanqFDPw I2=imread('F:\MATLAB学习\实验\picture\cell.tif'>。 %读入原图像文件DXDiTa9E3d I3=imread('cameraman.tif'>。 %读入原图像文件 subplot(3,2,1>。imshow(I1>。 %显示原图像 fftI1=fft2(I1>。 %二维离散傅立叶变换sfftI1=fftshift(fftI1>。 %直流分量移到频谱中心 RR1=real(sfftI1>。 %取傅立叶变换的实部II1=imag(sfftI1>。 %取傅立叶变换的虚部A1=sqrt(RR1.^2+II1.^2>。 %计算频谱幅值 A1=(A1-min(min(A1>>>/(max(max(A1>>-min(min(A1>>>*225。%归一化RTCrpUDGiT subplot(3,2,2>。imshow(A1>。 %显示原图像的频谱 subplot(3,2,3>。imshow(I2>。 %显示原图像 fftI2=fft2(I2>。 %二维离散傅立叶变换sfftI2=fftshift(fftI2>。 %直流分量移到频谱中心RR2=real(sfftI2>。 %取傅立叶变换的实部II2=imag(sfftI2>。 %取傅立叶变换的虚部A2=sqrt(RR2.^2+II2.^2>。 %计算频谱幅值

数字图像的傅里叶变换

数字图像的傅里叶变换 一. 课程设计目的 (1)了解图像变换的意义和手段 (2)熟悉傅里叶变换的基本性质 (3)热练掌握FFT的方法反应用 (4)通过本实验掌握利用MATLAB编程实现数字图像的傅里叶变换 二.课程设计要求 (1)熟悉并掌握傅立叶变换 (2)了解傅立叶变换在图像处理中的应用 (3)通过实验了解二维频谱的分布特点 (4)用MATLAB实现傅立叶变换仿真 三.设计思路 1.相关知识原理 (1)应用傅里叶变换进行数字图像处理 数字图像处理(digital image processing)是用计算机对图像信息进行处理的一门技术,使利用计算机对图像进行各种处理的技术和方法。 20世纪20年代,图像处理首次得到应用。20世纪60年代中期,随电子计算机的发展得到普遍应用。60年代末,图像处理技术不断完善,逐渐成为一个新兴的学科。利用数字图像处理主要是为了修改图形,改善图像质量,或是从图像中提起有效信息,还有利用数字图像处理可以对图像进行体积压缩,便于传输和保存。数字图像处理主要研究以下内容:傅立叶变换、小波变换等各种图像变换;对图像进行编码和压缩;采用各种方法对图像进行复原和增强;对图像进行分割、描述和识别等。随着技术的发展,数字图像处理主要应用于通讯技术、宇宙探索遥感技术和生物工程等领域。 傅里叶变换在数字图像处理中广泛用于频谱分析,傅里叶变换是线性系统分析的一个有力工具,它使我们能够定量地分析诸如数字化系统,采样点,电子放大器,卷积滤波器,噪声,显示点等地作用(效应)。傅里叶变换(FT)是数字图像处理技术的基础,其通过在时空域和频率域来回切换图像,对图像的信息特

二维傅里叶变换推倒及理解

2D 傅里叶变换理解心得 一、 目的 完整推倒2D 傅里叶变换公式,加深对2D 傅里叶变换公式的理解。 二、 内容 2维傅里叶变换,针对的信号函数是2维空间平面内的函数,2维傅里叶变换也有四种不同的形式。 1、 连续周期时域信号<---->非周期离散频谱。2D_CFS (,)XY f x y 表示2维周期连续信号,可以理解为一幅连续的图像信号(这里(,)XY f x y 可以为复数信号,但工程实践中常为实信号),(,)F k l 表示2维频谱信号,其中,k l 取-∞ +∞上的整数。 00000000002()2()00 00 2()2()0000 2()00 (,).(,).(,).1(,).,,-+X Y X Y j ku x lv y j ku x lv y XY XY X Y X Y j ku x lv y j ku x lv y X Y j ku x lv y XY f x y e dxdy f x y e dxdy F k l e e dxdy dxdy f x y e dxdy k l XY πππππ-+-++-+-+= = = ∞ ∞?? ?? ?????? 取上的实整数 其中X,Y 为(,)XY f x y 在x 坐标和y 坐标上各自的最小正周期。00,u v 表示在x 坐标和y 坐标上各自的基频率,这里有0011 ,u v X Y = =,,k l 取-∞+∞上的整数,对应不同的频率成分,(,) F k l 的图像为离散的,且在x 坐标和y 坐标上的频率间隔分别为0011 ,u v X Y = =。 002() (,)(,).,,-+j ku x lv y XY k l f x y F k l e x y π+∞ +∞ +=-∞=-∞ = ∞∞∑ ∑ 取上的实数 这里,(,)F k l 为复数。 所以得到2D_CFS (2维连续傅里叶级数) 00002() 002()(,).(,),,-+(,)(,).,,-+X Y j ku x lv y XY j ku x lv y XY k l f x y e dxdy F k l k l XY f x y F k l e x y ππ-++∞+∞ +=-∞=-∞????=∞∞???=∞∞?? ?? ∑∑取上的实整数 取上的实数

图像的二维傅里叶变换

图像傅立叶变换(二维傅立叶变换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变换。 h和k的范围在-N/2到N/2-1之间。 通常I(n,m)是实数,F(0,0)总是实数,并且F(h,k)具有对偶性。 如果写成指数形式,即: -------------------------------- 图像傅立叶变换的物理意义

图像的傅里叶变换

图像的傅里叶变换

成绩 实验论文 题目:图像的傅里叶变换 学生姓名:代朋车 学生学号:1114020207 系别:电气信息工程学院 专业:电子信息工程 年级:2011级 任课教师:沈晓波

电气信息工程学院制 2013年12月 图像的傅里叶变换 学生:代朋车 任课教师:沈晓波 电气信息工程学院 电子信息工程 1实验题目 图像的傅里叶变换 2实验对象 声音信号(自己的声音)/图像信号(自选)/自定义时域信号 3 实验任务 (1) 自行录制一段自己的语音或下载一段音乐,完成FFT 运算,具体参数自行规 定 基本要求:显示原语音信号图像、FFT 频谱,要求以上图像显示在同一个fig ,最终还原语音信号。 (2) 自行下载一幅图,格式不限,完成FFT 运算,具体参数可自行设定。 基本要求:显示原点在中心位置,画出幅度谱、相位谱能量或功率谱、虚部、实部,最终还原图像信号,要求以上图像显示在同一个fig 。 4 实验原理 4.1理论基础 连续傅里叶变换的定义:函数f(x)的一维连续傅里叶变换由下式定义: 2()()j ux F u f x e dx π∞ --∞=? 式中, 21J =-。 F(u)的傅里叶反变换定义为: 2()()j ux f x F u e du π∞ -∞==? 这里f(x)是实函数,它的傅里叶变换F(u)通常是复函数。F(u)的实部、虚部、振幅、

能量和相位分别表示如下:

所示。F(u)傅里叶频谱、相位谱和能量谱分别表示如下; F(u)=R(u)+jI(u) 傅里叶频谱 )()(R |F |22u I u u +=)( 相位谱 )(/)(arctan u R u I u =)(φ 能量谱 )()(|)(|222u I u R u F u +==E ) ( 同连续函数的傅里叶变换一样,离散函数的傅里叶变换也可推广到二维的情形,其二维离散傅里叶变换定义为: F(u,v)=∑∑==+1-N 0x 1-N 0y vy )/N (ux j2y)e F(x ,1πN 式中,u=0,1,2,…,N-1, v=0,1,2,…,N-1。二维离散傅里叶反变换定义为: F(x,y)=∑∑-=-=+101 0/)(2),(1N u N v N vy ux j e v u F N π 式中,x=0,1,2,…,N-1, y=0,1,2,…,N-1,u,v 是频率变量。与一维情况一样,二维函数的离散傅里叶频谱、相位谱和能量谱为: 傅里叶频谱 F(u,v)=),(),(22v u I v u R + 相位谱 ),(/),(arctan ,v u R v u I v u =)(φ 能量谱 ),(),(|),(|,222v u I v u R v u F v u +==E ) ( 4.2设计方案 4.2.1对音频进行FFT 运算 clear all; fs=8000; x1=wavread('安妮的仙境.wav'); t=(0:length(x1)-1)/8000; subplot(1,2,1); plot(t,x1) grid on;axis tight; title('原始语音信号'); xlabel('time(s)'); ylabel('幅度'); y1=fft(x1,2048); f=fs*(0:1023)/2048; %对信号做2048点FFT 变换 subplot(1,2,2);

图像傅里叶变换详解

图像傅里叶变换 冈萨雷斯版<图像处理>里面的解释非常形象:一个恰当的比喻是将傅里叶变换比作一个玻璃棱镜。棱镜是可以将光分解为不同颜色的物理仪器,每个成分的颜色由波长(或频率)来决定。傅里叶变换可以看作是数学上的棱镜,将函数基于频率分解为不同的成分。当我们考虑光时,讨论它的光谱或频率谱。同样, 傅立叶变换使我们能通过频率成分来分析一个函数。 Fourier theory讲的就是:任何信号(如图像信号)都可以表示成一系列正弦信号的叠加,在图像领域就是将图像brightness variation 作为正弦变量。比如下图的正弦模式可在单傅里叶中由三个分量编码:频率f、幅值A、相位γ这 三个value可以描述正弦图像中的所有信息。1.frequency frequency在空间域上可由亮度调节,例如左图的frequency比右图的frequency 低…… 2.幅值magnitude(amplitude)sin函数的幅值用于描述对比度,或者说是图像中最明和最暗的峰值之间的差。(一个负幅值表示一个对比逆转,即明暗交换。) 3.相位表示相对于原始波形,这个波形的偏移量(左or右)。=================================================================一个傅里叶变换编码是一系列正弦曲线的编码,他们的频率从0开始(即没有调整,相位为0,平均亮度处),到尼奎斯特频率(即数字图像中可被编码的最高频率,它和像素大小、resolution有关)。傅里叶变换同时将图像中所有频率进行编码:一个只包含一个频率f1的信号在频谱上横坐标f为f1的点处绘制一个单峰值,峰值高度等于对应的振幅amplitude,或者正弦曲线信号的高度。如下图所示。

图像的傅里叶变换

实验三 图像的傅里叶变换 一、 实验目的 1.了解图像变换的意义和手段; 2.掌握FFT 变换方法及应用; 3.通过实验了解二维频谱的分布特点; 4.通过本实验掌握利用MATLAB 编程实现数字图像的傅立叶变换。 二、 实验原理 1 应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具,它能够定量地分析诸如数字化系统、采样点、 电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能,将有助于解决 大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说,把时间用在 学习和掌握博里叶变换上是很有必要的。 2 傅立叶(Fourier )变换的定义 对于二维信号,二维Fourier 变换定义为: 2()(,)(,)j ux uy F u v f x y e dxdy π∞∞ -+-∞-∞= ?? 逆变换: 2()(,)(,)j ux uy f x y F u v e dudv π∞∞ +-∞-∞= ?? 二维离散傅立叶变换为: 11 2()00 1(,)(,)i k N N j m n N N i k F m n f i k e N π---+===∑∑ 逆变换: 11 2()00 1(,)(,)i k N N j m n N N m n f i k F m n e N π--+===∑∑ 三、 实验步骤及结果 步骤: 1将图像内容读入内存; 2用Fourier 变换算法,对图像作二维Fourier 变换; 3将其幅度谱进行搬移,在图像中心显示; 4用Fourier 系数的幅度进行Fourier 反变换;

5用Fourier系数的相位进行Fourier反变换; 6比较4、5的结果,评价人眼对图像幅频特性和相频特性的敏感度。 7记录和整理实验报告。 结果: 四、程序源代码 clear; I=imread('1.gif'); I=rgb2gray(I); subplot(3,3,1); imshow(I); title('1.gif'); E=fft2(double(I)); sfftI=fftshift(E); %正半轴部分和负半轴部分的图像分别关于各自的中心对称RR=real(sfftI); II=imag(sfftI); A=sqrt(RR.^2+II.^2); A=(A-min(min(A)))/(max(max(A))-min(min(A)))*225 ; subplot(3,3,2);

傅里叶变换及其在图像处理中的应用

傅里叶变换及其在数字图像处理中的应用 王家硕 学号:1252015 一、 Fourier 变换 1. 一维连续傅里叶变换 设 f (x)为x 的实变函数,如果f (x)满足下面的狄里赫莱条件: (1)具有有限个间隔点。 (2)具有有限个极点。 (3)绝对可积。 则 f (x )的傅里叶变换(Fourier Transformation ,FT )定义为: Fourier 正变换:dt e t f t f f F t j ? +∞ ∞ --==ωω)()]([)(; Fourier 逆变换:ωωπ ωd e f t F f t f t j ? ∞ +∞ ---= =)(21)]([)(1 , 式中:1-= j ,ω 为频域变量。 f (x )与F (w )构成傅里叶变换对,可以证明傅里叶变换对总是存在的。由于f (x )为实函数,则它的傅里叶变换F (w )通常是复函数,于是F (w )可写成 F (w ) = R (w ) + j I (w ) (1) 式中:R (w )和I (w )分别是F (w )的实部和虚部。公式1可表示为指数形式: 式中: F (w ) 为f (x )的傅里叶幅度谱,f (w )为f (x )的相位谱。 2. 二维连续傅里叶变换 如果二维函数f (x , y )是连续可积的,即∞

图像的傅立叶变换数字图像处理实验报告

图像的傅立叶变换数字图像处理实验报告

————————————————————————————————作者:————————————————————————————————日期:

计算机与信息工程学院验证性实验报告 专业:通信工程 年级/班级:2011级 2013—2014学年第一学期 课程名称 数字图像处理 指导教师 段新涛 本组成员 学号姓名 实验地点 计科楼111 实验时间 周五5-6节 项目名称 图像的傅立叶变换 实验类型 验证性 一、实验目的 1了解图像变换的意义和手段; 2熟悉傅立叶变换的基本性质; 3熟练掌握FFT 变换方法及应用; 4通过实验了解二维频谱的分布特点; 5通过本实验掌握利用MATLAB 编程实现数字图像的傅立叶变换。 6评价人眼对图像幅频特性和相频特性的敏感度。 二、实验原理 1 应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具,它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能,将有助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说,把时间用在学习和掌握博里叶变换上是很有必要的。 2 傅立叶(Fourier )变换的定义 对于二维信号,二维Fourier 变换定义为: 2()(,)(,)j ux uy F u v f x y e dxdy π∞∞ -+-∞-∞= ?? 逆变换: 2()(,)(,)j ux uy f x y F u v e dudv π∞∞ +-∞-∞= ?? 二维离散傅立叶变换为: 112()00 1(,)(,)i k N N j m n N N i k F m n f i k e N π---+===∑∑ 逆变换: 112()001(,)(,)i k N N j m n N N m n f i k F m n e N π--+===∑∑

傅立叶变换在图像处理中的作用

傅立叶变换在图像处理中的作用 摘要:本文首先简述了傅立叶变换的原理及应用领域,介绍了傅立叶变换在数字图象处理中的重要地位和应用,分析了其变换的数学原理和方法,特别着重的是二维傅立叶变换和FFT(快速傅立叶变换)的原理,然后介绍了Matlab 软件,分析了Matlab 的好处,及其在数字图像处理和傅立叶变换计算上的使用,编出程序实现了其变换功能,给出了应用于图象压缩和图像去噪的实例。 关键词: 图象处理 傅立叶变换 Matlab 正文 傅里叶变换是将时域信号分解为不同频率的正弦信号或余弦函数叠加之和。傅立叶变换是数字图像处理技术的基础,其通过在时空域和频率域来回切换图像,对图像的信息特征进行提取和分析,简化了计算工作量,被喻为描述图像信息的第二种语言,广泛应用于图像变换,图像编码与压缩,图像分割,图像重建等。因此,对涉及数字图像处理的工作者,深入研究和掌握傅立叶变换及其扩展形式的特性,是很有价值得。把傅立叶变换的理论通其物理解释相结合,将有助于解决大多数图像处理问题。傅里叶变换可分为连续傅里叶变换、离散傅里叶变换、快速傅里叶变换。 3.1.1 连续傅里叶变换 函数f(x)的傅里叶变换存在的条件是满足狄里赫莱条件,即: 1)具有有限个间断点; 2)具有有限个极值点; 3)绝对可积。 (1)一维连续傅里叶变换及反变换: 单变量连续函数f(x)的傅里叶变换F(u)定义为: dx e x f u F ux j ? ∞∞--=π2)()( 其中12-=j ,x 称为时域变量,u 为频率变量。 当给定F(u),通过傅里叶反变换可以得到f(x) du e u F x f ux j ?∞ ∞-=π2)()( (2)二维连续傅里叶变换及反变换: 二维连续函数f(x,y)的傅里叶变换F(u,v) 定义为:

5.图像的频域增强及傅里叶变换

5.图像的频域增强及傅里叶变换 傅立叶变换在图像处理中有非常非常的作用。因为不仅傅立叶分析涉及图像处理的很多方而,傅立叶的改进算法,比如离散余弦变换,gabor与小波在图像处理中也有重要的分量。 印象中,傅立叶变换在图像处理以下几个话题都有重要作用: 1.图像增强与图像去噪 绝大部分噪音都是图像的高频分呈:,通过低通滤波器来滤除髙频一一噪声;边缘也是图像的髙频分量,可以通过添加髙频分量来增强原始图像的边缘; 2?图像分割Z边缘检测 提取图像高频分量 3.图像特征提取: 形状特征:傅里叶描述子 纹理特征:直接通过傅里叶系数来汁算纹理特征 英他特征:将提取的特征值进行傅里叶变换来使特征具有平移、伸缩、旋转不变性 4.图像压缩 可以直接通过傅里叶系数来压缩数据:常用的离散余弦变换是傅立叶变换的实变换:傅立叶变换傅里叶变换是将时域信号分解为不同频率的正弦信号或余弦函数叠加之和。连续情况下要求原始信号在一个周期内满足绝对可积条件。离散情况下,傅里叶变换一左存在。冈萨雷斯版<图像处理>里而的解释非常形象:一个恰当的比喻是将傅里叶变换比作一个玻璃棱镜。棱镜是可以将光分解为不同颜色的物理仪器,每个成分的颜色由波长(或频率)来决泄。傅里叶变换可以看作是数学上的棱镜,将函数基于频率分解为不同的成分。当我们考虑光时, 讨论它的光谱或频率谱。同样,傅立叶变换使我们能通过频率成分来分析一个函数。 傅立叶变换有很多优良的性质。比如线性,对称性(可以用在计算信号的傅里叶变换里而);时移性:函数在时域中的时移,对应于其在频率域中附加产生的相移,而幅度频谱则保持不变;频移性:函数在时域中乘以』wt,可以使整个频谱搬移W U这个也叫调制左理,通讯里而信号的频分复用需要用到这个特性(将不同的信号调制到不同的频段上同时传输): 卷积泄理:时域卷积等于频域乘枳:时域乘积等于频域卷积(附加一个系数)。(图像处理里而这个是个重点)信号在频率域的表现在频域中,频率越大说明原始信号变化速度越快:频率越小说明原始信号越平缓。当频率为O时,表示直流信号,没有变化。因此,频率的大小反应了信号的变化快慢。高频分疑解释信号的突变部分,而低频分量决左信号的整体形象。 在图像处理中,频域反应了图像在空域灰度变化剧烈程度,也就是图像灰度的变化速度, 也就是图

5.图像的频域增强及傅里叶变换

5. 图像的频域增强及傅里叶变换 傅立叶变换在图像处理中有非常非常的作用。因为不仅傅立叶分析涉及图像处理的很多方面,傅立叶的改进算法,比如离散余弦变换,gabor与小波在图像处理中也有重要的分量。 印象中,傅立叶变换在图像处理以下几个话题都有重要作用: 1.图像增强与图像去噪 绝大部分噪音都是图像的高频分量,通过低通滤波器来滤除高频——噪声; 边缘也是图像的高频分量,可以通过添加高频分量来增强原始图像的边缘; 2.图像分割之边缘检测 提取图像高频分量 3.图像特征提取: 形状特征:傅里叶描述子 纹理特征:直接通过傅里叶系数来计算纹理特征 其他特征:将提取的特征值进行傅里叶变换来使特征具有平移、伸缩、旋转不变性 4.图像压缩 可以直接通过傅里叶系数来压缩数据;常用的离散余弦变换是傅立叶变换的实变换;傅立叶变换 傅里叶变换是将时域信号分解为不同频率的正弦信号或余弦函数叠加之和。连续情况下要求原始信号在一个周期内满足绝对可积条件。离散情况下,傅里叶变换一定存在。冈萨雷斯版<图像处理>里面的解释非常形象:一个恰当的比喻是将傅里叶变换比作一个玻璃棱镜。棱镜是可以将光分解为不同颜色的物理仪器,每个成分的颜色由波长(或频率)来决定。傅里叶变换可以看作是数学上的棱镜,将函数基于频率分解为不同的成分。当我们考虑光时,讨论它的光谱或频率谱。同样,傅立叶变换使我们能通过频率成分来分析一个函数。 傅立叶变换有很多优良的性质。比如线性,对称性(可以用在计算信号的傅里叶变换里面);时移性:函数在时域中的时移,对应于其在频率域中附加产生的相移,而幅度频谱则保持不变;频移性:函数在时域中乘以,可以使整个频谱搬移w。这个也叫调制定理,通讯里面信号的频分复用需要用到这个特性(将不同的信号调制到不同的频段上同时传输);卷积定理:时域卷积等于频域乘积;时域乘积等于频域卷积(附加一个系数)。(图像处理里面这个是个重点) 信号在频率域的表现 在频域中,频率越大说明原始信号变化速度越快;频率越小说明原始信号越平缓。当频

图像的傅里叶变换

一、 实验目的 1.了解图像变换的意义和手段; 2.掌握FFT 变换方法及应用; 3.通过实验了解二维频谱的分布特点; 4.通过本实验掌握利用MATLAB 编程实现数字图像的傅立叶变换。 二、 实验原理 1 应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具,它能够定量地分析诸如数字化系统、采样点、 电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能,将有助于解决 大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说,把时间用在 学习和掌握博里叶变换上是很有必要的。 2 傅立叶(Fourier )变换的定义 对于二维信号,二维Fourier 变换定义为: 2()(,)(,)j ux uy F u v f x y e dxdy π∞∞ -+-∞-∞= ?? 逆变换: 2()(,)(,)j ux uy f x y F u v e dudv π∞∞ +-∞-∞= ?? 二维离散傅立叶变换为: 11 2()00 1(,)(,)i k N N j m n N N i k F m n f i k e N π---+===∑∑ 逆变换: 11 2()00 1(,)(,)i k N N j m n N N m n f i k F m n e N π--+===∑∑ 三、 实验步骤及结果 步骤: 1将图像内容读入内存; 2用Fourier 变换算法,对图像作二维Fourier 变换; 3将其幅度谱进行搬移,在图像中心显示; 4用Fourier 系数的幅度进行Fourier 反变换; 5用Fourier 系数的相位进行Fourier 反变换; 6比较4、5的结果,评价人眼对图像幅频特性和相频特性的敏感度。 7记录和整理实验报告。 结果:

二维离散傅立叶变换

图像的二维离散傅立叶变换 一、实验目的 掌握图像的二维离散傅立叶变换以及性质 二、实验要求 1) 建立输入图像,在64?64的黑色图像矩阵的中心建立16?16的白色矩形图像点阵, 形成图像文件。对输入图像进行二维傅立叶变换,将原始图像及变换图像(三维、中心化)都显示于屏幕上。 2) 调整输入图像中白色矩形的位置,再进行变换,将原始图像及变换图像(三维、中 心化)都显示于屏幕上,比较变换结果。 3) 调整输入图像中白色矩形的尺寸(40?40,4?4),再进行变换,将原始图像及变 换图像(三维、中心化)都显示于屏幕上,比较变换结果。 三、实验仪器设备及软件 HP D538、MATLAB 四、实验原理 设),(y x f 是在空间域上等间隔采样得到的M ×N 的二维离散信号,x 和y 是离散实变量,u 和v 为离散频率变量,则二维离散傅里叶变换对一般地定义为 ∑∑-=-=+-=1010)],(2exp[),(1),(M x N y N yu M xu j y x f MN v u F π,1,0=u …,M-1;y=0,1,…N-1 ∑∑-=-=+=101 0)],( 2e x p [),(),(M x N y N uy M ux j v u F y x f π ,1,0=x …,M-1;y=0,1,…N-1 在图像处理中,有事为了讨论上的方便,取M=N ,这样二维离散傅里叶变换对就定义为,])(2exp[),(1),(1010∑∑-=-=+-=N x N y N yu xu j y x f N v u F π 1,0,=v u …,N-1 ,])(2exp[),(1),(1010∑∑-=-=+=N u N v N vy ux j v u F N y x f π 1,0,=y x ,…,N-1 其中,]/)(2exp[N yv xu j +-π是正变换核,]/)(2exp[N vy ux j +π是反变换核。 将二维离散傅里叶变换的频谱的平方定义为),(y x f 的功率谱,记为

图像处理之傅里叶变换matlab实现

傅里叶变换 一.实验内容: 1、傅里叶变换 二.实验目的: 1、理解傅里叶变换的原理 2、掌握傅里叶变换的性质 三.实验步骤: 1.首先构造一幅黑白二值图像,在128×128的黑色背景中心产生一个4×4的白 色方块,对其进行傅里叶变换;(Matlab 中用fft2实现2D 傅里叶变换) 2.把低频分量移到图象中心,而把高频分量移到四个角上;(方法有两种:其 一,在FT 以前对测试图象逐点加权(-1)^(i+j);其二,利用FFTSHIFT 函数); 3.利用图象增强中动态范围压缩的方法增强2DFT ;(Y =C*log (1+abs (X))); 4.构造一幅黑白二值图像,在128×128的黑色背景中令第32行至36行、第 32列至第36列的值为1(即产生一个4×4的白色方块),对其进行傅里叶变换; 5.将上图旋转300,再进行傅里叶变换 (imrotate ) 6.构造二幅黑白二值图像,在128×128的黑色背景中分别令第60行至68行、 第60列至第68列的值为1,第64行至65行、第64列至第65列的值为1产生两幅图像,分别对这两幅图像进行傅里叶变换 四、原理分析、技术讨论、回答问题 1、对于第二幅图像(第一步与第四步图像的比较),说明FOURIER 变换具有以下性质: )//(20000),(),(N vy M ux j e v u F y y x x f +-?--π 2、对于第三幅图像(第一步与第五步图像的比较),说明FOURIER 变换具有以下性质: θcos r x = θs i n r y = αωc o s =u αωs i n =v ),(),(00θαωθθ+?+F r f 3、对于第四幅图像(第一步与第六步图像的比较),说明FOURIER 变换具有以下性质: )/,/(||1 ),(b v a u F ab by ax f =

傅里叶变换

少将 图像的傅立叶变换,原始图像由N行N列构成,N必须是基2的,把这个N*N个包含图像的点称为实部,另外还需要N*N个点称为虚部,因为FFT是基于复数的,如下图所示: (//实数DFT将时域内的N个点变换为频域中两组各N/2+1个点(分别对应实部和虚部)) 计算图像傅立叶变换的过程很简单:首先对每一行做一维FFT,然后对每一列做一维FFT。具体来说,先对第0行的N个点做FFT(实部有值,虚部为0),将FFT输出的实部放回原来第0行的实部,FFT输出的虚部放回第0行的虚部,这样计算完全部行之后,图像的实部和虚部包含的是中间数据,然后用相同的办法进行列方向上的FFT变换,这样N*N的图像经过FFT得到一个N*N的频谱。 下面展示了一副图像的二维FFT变换:

频域中可以包含负值,图像中灰色表示0,黑色表示负值,白色表示正值。可以看到4个角上的黑色更黑,白色更白,表示其幅度更大,其实4个角上的系数表示的是图像的低频组成部分,而中心则是图像的高频组成部分。除此以外,FFT的系数显得杂乱无章,基本看不出什么。 将上述直角坐标转换为极坐标的形式,稍微比较容易理解一点,幅度中4个角上白色的区域表示幅度较大,而相位中高频和低频基本看不出什么区别来。 上述以一种不同的方法展示了图像频谱,它将低频部分平移到了频谱的中心(//MATLAB中实现函数fftshift)。这个其实很好理解,因为经2D-FFT的信号是离散图像,其2D-FFT的输

出就是周期信号,也就是将前面一张图周期性平铺,取了一张以低频为中心的图。将原点放在中心有很多好处,比如更加直观更符合周期性的原理,但在这节中还是以未平移之前的图来解释。 行N/2和列N/2将频域分成四块。对实部和幅度来说,右上角和左下角成镜像关系,左上角和右下角也是镜像关系;对虚部和相位来说,也是类似的,只是符号要取反(//共轭?),这种对称性和1维傅立叶变换是类似的,你可以往前看看。 为简单起见,先考虑4*4的像素,右边是其灰度值,对这些灰度值进行2维fft变换。 h和k的范围在-N/2到N/2-1之间。 通常I(n,m)是实数,F(0,0)总是实数(//直流分量),并且F(h,k)具有对偶性。 如果写成复数形式,即: ------------------------------------------------------------- 图像傅立叶变换的物理意义

相关主题
文本预览
相关文档 最新文档