图象变换正交变换傅立叶变换.pptx

《图象变换正交变换傅立叶变换.pptx》由会员分享，可在线阅读，更多相关《图象变换正交变换傅立叶变换.pptx（58页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、12023年2月10日一、正交变换第1页/共58页22023年2月10日连续函数集合的正交性正交函数集合当C C=1=1时，称集合为归一化正交函数集合第2页/共58页32023年2月10日正交函数集合的完备性若f(x)是定义在t0和t0+T区间的实值信号，平方可积。可以表示为：对任意小的0，存在充分大的N，其中，则称函数U集合是完备的。意味着f(x)可以由无穷级数来表示第3页/共58页42023年2月10日离散情况n个正交向量 当C=1时，称归一化正交。即每一个相量为单位相量。第4页/共58页52023年2月10日满足上式的基相量组成矩阵：则一定满足：第5页/共58页62023年2月10日一维

2、正交变换对于一向量f，用上述正交矩阵进行运算：g=Af若要恢复f，则:以上过程称为正交变换。我们把原为A-1可以用AT来代替的A阵称为正交矩阵。第6页/共58页72023年2月10日二维正交变换 NN二维函数可以类似于一维正变换核反变换核l显然，这两个变换核应该满足正交性和完备性。第7页/共58页82023年2月10日二、傅立叶变换第8页/共58页92023年2月10日傅立叶变换傅立叶变换域也称为频域变换，它把图像从图像空间变换到频率空间。将原定义在图像空间的图像以某种形式转换（正变换）到另外一些空间，并利用在这些空间的特有性质方便地进行一定的加工，最后再转换回图像空间（反变换或逆变换）以得到

3、所需要的效果。第9页/共58页102023年2月10日1.连续函数的傅立叶变换 若把一个一维输入信号作一维傅立叶变换，该信号就被变换到频域上的一个信号，即得到了构成该输入信号的频谱，频谱反映了该输入信号由哪些频率构成。当一个一维信号f(x)满足狄里赫莱条件，即f(x)（1）具有有限个间断点；（2）具有有限个极值点；（3）绝对可积。则其傅立叶变换对（正变换和逆变换）一定存在。第10页/共58页112023年2月10日一维傅立叶变换的定义f(x)为连续可积函数，其傅立叶变换定义为：其反变换为：式中：，x称为时域变量，u为频域变量。通常傅立叶变换为复数形式F(u)=R(u)+jI(u)幅度谱：相位谱

4、：第11页/共58页122023年2月10日变换分析的直观说明 把一个信号的波形分解为许多不同频率正弦波之和。第12页/共58页132023年2月10日一维傅立叶变换举例方波信号：经过傅立叶变换后：第13页/共58页142023年2月10日一维离散傅立叶变换(DFT)一维离散傅立叶变换公式为：逆变换为：数学上建立傅立叶变换的f(x)是连续的模拟信号，而计算机处理的是离散的数字信号，同时数学上用无穷大概念，而计算机只能进行有限次计算。通常就将这种受限的傅立叶变换称为离散傅立叶变换（DFT）。第14页/共58页152023年2月10日由欧拉公式可知可得 可见，离散序列的傅立叶变换仍然是一个离散的序

5、列，每一个u对应的傅立叶变换结果是所有输入序列f(x)的加权和。每个f(x)都乘以不同频率的正弦和余弦值。第15页/共58页162023年2月10日二维离散傅立叶变换对于二维傅立叶变换，由一维推广而来，其离散形式为：逆变换为：幅谱(频谱)、相谱：第16页/共58页172023年2月10日二维傅立叶变换举例对于二维方波信号傅立叶变换为：幅度：第17页/共58页182023年2月10日例：函数在以原点为中心的一个正方形内为正值常数，而在其它地方为零。傅立叶频谱幅度的灰度图显示。第18页/共58页192023年2月10日二维傅立叶变换的性质 1.分离性 二维二维傅立叶变换可由连续两次运用一维傅立叶变

6、换来实现。第19页/共58页202023年2月10日 由可分离性可知，一个二维傅立叶变换可分解为两步进行，其中每一步都是一个一维傅立叶变换。先对f(x,y)按列进行傅立叶变换得到F(x,v)，再对F(x,v)按行进行傅立叶变换，便可得到f(x,y)的傅立叶变换结果，如图所示。显然对f(x,y)先按行进行离散傅立叶变换，再按列进行离散傅立叶变换也是可行的。第20页/共58页212023年2月10日2.平移性 将f(x,y)f(x,y)与一个指数项相乘就相当于把其变换后的频域中心移动到新的位置。F(u,v)F(u,v)与一个指数项相乘就相当于把其反变换后的空域中心移动到新的位置。对f(x,y)f(

7、x,y)的平移不影响其傅立叶变换的幅值。第21页/共58页222023年2月10日3.周期性和共扼对称性 如果f(x,y)f(x,y)是实函数，则它的傅立叶变换具有共扼对称性：F F*(u,v)(u,v)为F(u,v)F(u,v)的复共扼。假定傅立叶变换和反变换均以N为周期。第22页/共58页232023年2月10日 由旋转不变性可知，如果时域中离散函数旋转0角度，则在变换域中该离散傅立叶变换函数也将旋转同样的角度0。离散傅立叶变换的旋转不变性如图3-3所示。图3-3 离散傅立叶变换的旋转不变性（a）原始图像；（b）原图像的傅立叶频谱；（c）旋转45后的图像；（d）图像旋转后的傅立叶频谱(a)

8、(b)(d)(c)4.旋转性质 借助极坐标变换 x=rcosx=rcos，y=rsiny=rsin，u=w cosu=w cos，v=w v=w sinsin，将f(x,y)f(x,y)和F(u,v)F(u,v)转换为f(r,)f(r,)和F(w,)F(w,)。第23页/共58页242023年2月10日傅立叶变换和反变换对加法满足分配律，但对乘法则不满足。5.分配律 6.尺度变换(缩放)比例性质第24页/共58页252023年2月10日将u=v=0u=v=0代入正变换式，可以得到：7 7.平均值 2 2个函数的卷积定义为：8.卷积平均值计算计算步骤:折叠位移相乘积分（求和）第25页/共58页图

9、像变换(二)快速傅立叶变换、离散余弦变换第26页/共58页272023年2月10日普通傅立叶变换：完成全部DFT运算的计算量与N2成正比。特别是当N较大时，其运算时间将迅速增长，以至于无法容忍。为此，研究离散傅立叶变换的快速算法（Fast Fourier Transform，FFT）非常必要。研究快速傅立叶变换的必要性第27页/共58页 快速离散傅立叶变换一种称为逐次加倍法的快速傅立叶变换算法（FFT），它是1965年Cooley和Tukey首先提出的。算法时间复杂度为Nlog2N。当N很大时计算量可以大大减少。第28页/共58页292023年2月10日记称为旋转因子。则有：单位圆表示：快速傅

10、立叶变换(3-1)第29页/共58页302023年2月10日式中，由Wux构成的矩阵称为W阵或系数矩阵。一维离散傅立叶变换（DFT）用矩阵的形式表示为：第30页/共58页312023年2月10日 W的定义表达式W=e-j2N，由欧拉公式知系数W是以N为周期的。这样，W阵中很多系数就是相同的，且由于W的对称性，即 因此可进一步减少计算工作量。例如，对于N=4，W阵为 W4W0，W6W2，W9W1；W3W1，W2W0第31页/共58页322023年2月10日 可见N=4的W阵中只需计算W0和W1两个系数即可。说明W阵的系数有许多计算工作是重复的，如果把一个离散序列分解成若干短序列，并充分利用旋转因

11、子W的周期性和对称性来计算离散傅立叶变换，便可以简化运算过程，这就是FFT的基本思想。设N为2的正整数次幂，即 如令M为正整数，且 N=2M 第32页/共58页332023年2月10日将式N=2M代入式（3-1），离散傅立叶变换可改写成如下形式：由旋转因子W的定义可知，因此式(3-2)变为 现定义(3-2)(3-3)(3-4)(3-5)第33页/共58页342023年2月10日于是式(3-3)变为(3-6)进一步考虑W的对称性和周期性可知 和，于是(3-7)由此，可将一个N点的离散傅立叶变换分解成两个N2短序列的离散傅立叶变换，即分解为偶数和奇数序列的离散傅立叶变换Fe(u)和Fo(u)。第3

12、4页/共58页352023年2月10日 在此，以计算N=8的DFT为例，此时n=3，M=4。由式(3-6)和式(3-7)可得(3-8)第35页/共58页362023年2月10日 式（3-8）中，u取07时的F(u)、Fe(u)和Fo(u)的关系可用图3.1描述。左方的两个节点为输入节点，代表输入数值；右方两个节点为输出节点，表示输入数值的叠加，运算由左向右进行。线旁的W18和W18为加权系数，定义由F(1)、F(5)、Fe(1)和Fo(1)所构成的结构为蝶形运算单元,其表示的运算为(3-9)图3.1 蝶形运算单元第36页/共58页372023年2月10日 由于Fe(u)和Fo(u)都是4点的D

13、FT，因此，如果对它们再按照奇偶进行分组，则有第37页/共58页382023年2月10日图3.2 4点DFT分解为2点DFT的蝶形流程图 第38页/共58页392023年2月10日二维快速傅立叶变换的Matlab实现简单图像及其傅立叶变换Eg3.4d=zeros(32,32);d(13:20,13:20)=1;figure(1);imshow(d,notruesize);D=fft2(d);figure(2);imshow(abs(D),-1 5,notruesize);第39页/共58页402023年2月10日第40页/共58页412023年2月10日二维快速傅立叶变换的Matlab实现eg

14、3.3figure(1);load imdemos saturn2;imshow(saturn2);figure(2);S=fftshift(fft2(saturn2);imshow(log(abs(S),);第41页/共58页422023年2月10日频域变换的一般表达式 1、可分离变换 二维傅立叶变换可用通用的关系式来表示：(3-10)(3-11)式中：x,u=0,1,2,M1；y,v=0,1,2,N1；g(x,y,u,v)和h(x,y,u,v)分别称为正向变换核和反向变换核。第42页/共58页432023年2月10日如果 g(x,y,u,v)=g1(x,u)g2（y,v）(3-12)h(x

15、,y,u,v)=h1(x,u)h2(y,v)(3-13)则称正、反变换核是可分离的。如果g1和g2，h1和h2在函数形式上一样，则称该变换核是对称的。二维傅立叶变换对是一个特殊情况，它们的核为 可分离对称第43页/共58页442023年2月10日 对于图像变换，只要其变换核是可分离的，就可用两次一维变换来实现。如果先对f(x,y)的每一列进行一维变换得到F(y,u)，再沿F(y,u)每一行取一维变换得到F(u,v)，和先对f(x,y)的每一行进行一维变换得到F(x,v)，再沿F(x,v)每一列取一维变换得到F(u,v)其最终结果是一样的。该结论对反变换核也适用。第44页/共58页452023年

16、2月10日3.4 离散余弦变换第45页/共58页462023年2月10日3.4.1 一维离散余弦变换 一维DCT的变换核定义为 式中，x,u=0,1,2,N1；(3-47)(3-48)一维DCT定义如下：设f(x)|x=0,1,N-1为离散的信号列。(3-49)式中，u,x=0,1,2,N1。第46页/共58页472023年2月10日将变换式展开整理后，可以写成矩阵的形式，即 F=Gf（3-50）其中(3-51)第47页/共58页482023年2月10日 一维DCT的逆变换IDCT定义为(3-52)式中，x,u=0,1,2,N1。可见一维DCT的逆变换核与正变换核是相同的。第48页/共58页4

17、92023年2月10日3.4.2 二维离散余弦变换将一维DCT的定义推广到二维DCT。其正变换核为 (3-53)式中，C(u)和C(v)的定义同式（3-48）；x,u=0,1,2,M1；y,v=0,1,2,N1。第49页/共58页502023年2月10日3.4.2 二维离散余弦变换二维DCT定义如下：设f(x,y)为MN的数字图像矩阵，则(3-54)式中：x,u=0,1,2,M1；y,v=0,1,2,N1。第50页/共58页512023年2月10日 二维DCT逆变换定义如下：(3-55)式中：x,u=0,1,2,M1；y,v=0,1,2,N1。类似一维矩阵形式的DCT，可以写出二维DCT的矩阵

18、形式如下：F=GfGT (3-56)第51页/共58页522023年2月10日 同时，由式(3-55)和式(3-54)可知二维DCT的逆变换核与正变换核相同，且是可分离的，即(3-57)式中：C(u)和C(v)的定义同式（3-48）；x,u=0,1,2,M1；y,v=0,1,2,N1。第52页/共58页532023年2月10日 通常根据可分离性，二维DCT可用两次一维DCT来完成，其算法流程与DFT类似，即 第53页/共58页542023年2月10日3.4.3 快速离散余弦变换 离散余弦变换的计算量相当大，在实用中非常不方便，也需要研究相应的快速算法。目前已有多种快速DCT（FCT），在此介绍

19、一种由FFT的思路发展起来的FCT。首先，将f(x)延拓为 x=0,1,2,N-1x=N，N+1,2N-1(3-59)按照一维DCT的定义，fe(x)的DCT为(3-60)第54页/共58页552023年2月10日式中，Re表示取复数的实部。第55页/共58页562023年2月10日 由于 为fe(x)的2N点DFT。因此，在作DCT时，可把长度为N的f(x)的长度延拓为2N点的序列fe(x)，然后对fe(x)作DFT，最后取DFT的实部便可得到DCT的结果。同理对于离散余弦逆变换IDCT，可首先将F(u)延拓为u=0,1,2,N-1u=N,N+1,2N-1(3-62)第56页/共58页572023年2月10日由式(3-52)可得，DCT的IDCT为（3-63）由式（3-63）可见，IDCT可由 的2N点的IDFT来实现。第57页/共58页第三章 图像变换582023年2月10日感谢您的观看。第58页/共58页