第五章图象正交变换PPT讲稿.ppt
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《第五章图象正交变换PPT讲稿.ppt》由会员分享,可在线阅读,更多相关《第五章图象正交变换PPT讲稿.ppt(129页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第五章图象正交变换第1页,共129页,编辑于2022年,星期三什么是图象变换什么是图象变换将图象看成是线性叠加系统将图象看成是线性叠加系统图象在空域上相关性很强图象在空域上相关性很强图象变换是将图象从空域变换到其它域图象变换是将图象从空域变换到其它域如频域的数学变换如频域的数学变换常用的变换:傅立叶变换、离散余弦变常用的变换:傅立叶变换、离散余弦变换、沃尔什变换、离散换、沃尔什变换、离散K-L变换、小波变变换、小波变换等换等第2页,共129页,编辑于2022年,星期三连续函数集合的正交性连续函数集合的正交性正交函数集合正交函数集合当当C=1时,时,称集合为归一化正交函数集合称集合为归一化正交函
2、数集合 第3页,共129页,编辑于2022年,星期三正交函数集合的完备性正交函数集合的完备性若若f(x)是定义在是定义在t0和和t0+T区间的实值信号区间的实值信号,平,平方可积。可以表示为:方可积。可以表示为:对任意小的对任意小的0,存在充分大的,存在充分大的N,其中其中,则称函数则称函数U集合是完备的。集合是完备的。第4页,共129页,编辑于2022年,星期三离散情况离散情况n个正交向量个正交向量当当C=1时,时,称归一化正交称归一化正交 第5页,共129页,编辑于2022年,星期三满足:满足:第6页,共129页,编辑于2022年,星期三一维正交变换一维正交变换对于一向量对于一向量f,用上
3、述正交矩阵进行,用上述正交矩阵进行运算:运算:g=Af若要恢复若要恢复f,则,则以上过程称为以上过程称为正交变换正交变换。第7页,共129页,编辑于2022年,星期三酉变换酉变换若若A为复数矩阵,正交的条件为:为复数矩阵,正交的条件为:其中其中A*为为A的复数共轭矩阵,满足这个的复数共轭矩阵,满足这个条件的矩阵为条件的矩阵为酉矩阵酉矩阵。对于任意向量。对于任意向量f的运算称为的运算称为酉变换酉变换:第8页,共129页,编辑于2022年,星期三酉变换、正交变换与信号分析酉变换、正交变换与信号分析正交变换是酉变换的特例正交变换是酉变换的特例它们都可以用于信号分析它们都可以用于信号分析用于信号分析的
4、基函数集合和正交矩阵用于信号分析的基函数集合和正交矩阵都应满足正交性和完备性都应满足正交性和完备性第9页,共129页,编辑于2022年,星期三二维酉变换二维酉变换 NN二维函数可以类似于一维用正交序列展开和恢二维函数可以类似于一维用正交序列展开和恢复复正变换核正变换核反变换核反变换核第10页,共129页,编辑于2022年,星期三变换核的可分离性变换核的可分离性其中其中au(x),u=0,1,N-1,bv(y),v=0,1,N-1为为一维完备正交基向量的集合。用矩阵表示:一维完备正交基向量的集合。用矩阵表示:A=a(u,x),B=b(v,y)通常选择通常选择A=B。第11页,共129页,编辑于2
5、022年,星期三二维酉变换二维酉变换A=B时,二维酉变换正变换表示为时,二维酉变换正变换表示为用矩阵表示:用矩阵表示:F=AfAT类似的,对于类似的,对于MN的二维函数的二维函数f(x,y)第12页,共129页,编辑于2022年,星期三基图象基图象反变换反变换看成是基图象看成是基图象F(u,v)权因子权因子图象图象f(x,y)可以用可以用N2个基图象的加权和来表示个基图象的加权和来表示第13页,共129页,编辑于2022年,星期三酉变换的性质酉变换的性质1.酉矩阵是正交阵酉矩阵是正交阵2.AA*T=A*TA=INN3.2.A为酉阵,则为酉阵,则A-1和和AT都是酉阵都是酉阵4.3.酉变换是能量
6、保持的变换酉变换是能量保持的变换5.对于一维酉变换对于一维酉变换F=Af,有有|F|=|f|6.二维情况下,则有:二维情况下,则有:第14页,共129页,编辑于2022年,星期三酉变换的性质酉变换的性质(2)设设f(x,y)的均值和协方差为的均值和协方差为f和和f4.均值和方差均值和方差则则F(u,v)的均值为:的均值为:则则F(u,v)的协方差为:的协方差为:第15页,共129页,编辑于2022年,星期三酉变换的性质酉变换的性质(3)5.其他性质:其他性质:(1)A为酉阵,则其行列式值为酉阵,则其行列式值|A|=1(2)若若a为向量,则作酉变换后向量为向量,则作酉变换后向量模保持不变:模保持
7、不变:b=Aa,则,则|b|=|a|。第16页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 5.1.1 一维傅立叶变换一维傅立叶变换1.一维连续函数的傅立叶变换(一维连续函数的傅立叶变换(FT)定义:若函数定义:若函数f(x)满足条件:)满足条件:1)具有有限个间断点;)具有有限个间断点;2)具有有限个极值点;)具有有限个极值点;3)绝对可积,)绝对可积,则把变换称为:则把变换称为:傅立叶正变换:傅立叶正变换:傅立叶反变换:傅立叶反变换:傅立叶变换对:傅立叶变换对:F(u)f(x)第17页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 2.一维离散傅立
8、叶变换(一维离散傅立叶变换(DFT)傅立叶正变换:傅立叶正变换:傅立叶反变换:傅立叶反变换:运算量为运算量为N*N次复数相乘和次复数相乘和N*(N-1)次复数相加次复数相加第18页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 2.一维离散傅立叶变换(一维离散傅立叶变换(DFT)实序列的实序列的FT:第19页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 目的:(1)用矩阵乘法的程序进行FT;(2)理论推导用。一维一维DFT的矩阵表示的矩阵表示根据定义:令:则:展开:第20页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 一维一维DF
9、T的矩阵表示的矩阵表示当当N=4时:时:第21页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 一维一维DFT的矩阵表示的矩阵表示例:例:第22页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 系数 和 具有 周期 性和对称性 傅立叶系数的性质:傅立叶系数的性质:第23页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 时间抽取基时间抽取基2 FFT算法算法 要求N为2的幂,即 ,M为正整数将输入时间序列按奇、偶抽取两个N/2点的FT3.快速傅立叶变换(快速傅立叶变换(FFT)第24页,共129页,编辑于2022年,星期三5.1 傅立叶
10、变换傅立叶变换 时间抽取基时间抽取基2 FFT算法算法对于uN/2-1的点,有以下规则给出:N=84点DFT4点DFT第25页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 时间抽取基时间抽取基2 FFT算法算法蝶形运算蝶形运算(u=0,1,N/2-1)第26页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 时间抽取基时间抽取基2 FFT算法算法用两个2点DFT代替4点DFT:第27页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 时间抽取基时间抽取基2 FFT算法算法2点DFT2点DFT2点DFT2点DFT输入倒序排列第28页,共1
11、29页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 时间抽取基时间抽取基2 FFT算法算法输入倒序排列第29页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 标号二进制表示 按位倒序的二进制表示倒序后标号0123456700000101001110010111011100010001011000110101111104261537时间抽取基时间抽取基2 FFT算法算法第30页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 4.用计算机实现快速傅里叶变换用计算机实现快速傅里叶变换利用利用FFT蝶形流程图算法在计算机上实现快速付里叶变换必须解决如下
12、蝶形流程图算法在计算机上实现快速付里叶变换必须解决如下问题:迭代次数问题:迭代次数r的确定、对偶节点的计算、加权系数的确定、对偶节点的计算、加权系数 的计算、重新的计算、重新排序问题。排序问题。(1)迭代次数迭代次数r的确定的确定 由蝶形流程图可见,迭代次数与由蝶形流程图可见,迭代次数与N有关。有关。r值可由下式确定:值可由下式确定:式中式中N是变换序列的长度,对于基是变换序列的长度,对于基2的蝶形程图,的蝶形程图,N是是2的整的整数次幂。例如,序列长度为数次幂。例如,序列长度为8则要三次迭代,序列长度为则要三次迭代,序列长度为16时就要时就要4次迭代。次迭代。第31页,共129页,编辑于20
13、22年,星期三5.1 傅立叶变换傅立叶变换 4.用计算机实现快速傅里叶变换用计算机实现快速傅里叶变换(2)对偶节点的计算对偶节点的计算 在流程图中把标有在流程图中把标有xi(k)的点称为节点。其中下标的点称为节点。其中下标i为列数,也就是为列数,也就是第几次迭代,例如,第几次迭代,例如,x1(k)则说明它是第一次迭代的结果。则说明它是第一次迭代的结果。k代表代表流程图中的行数,也就是序列的序号数。其中每节点的值均是流程图中的行数,也就是序列的序号数。其中每节点的值均是用前一用前一 节点对计算得来的。例如,节点对计算得来的。例如,x1(0)和和x1(4)均是均是x(0)和和x(4)计计算得来的。
14、在蝶形流程图中,把具有相同来源的一对节点叫做对偶节点。算得来的。在蝶形流程图中,把具有相同来源的一对节点叫做对偶节点。如:如:x1(0)和和x1(4)就是一对对偶节点,因为它们均来源于就是一对对偶节点,因为它们均来源于x(0)和和x(4)。对偶节点的计算也就是求出在每次迭代小对偶节点的间隔。由。对偶节点的计算也就是求出在每次迭代小对偶节点的间隔。由流程图可见,第一次迭代的间距为计流程图可见,第一次迭代的间距为计N2,第二次迭代的节距为,第二次迭代的节距为N4,第三次迭代的节距为人,第三次迭代的节距为人N23等等。由以上分析可得到对等等。由以上分析可得到对偶节点的计算方法。如果某一节点为偶节点的
15、计算方法。如果某一节点为xi(k),那么它的对偶节点为,那么它的对偶节点为第32页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 4.用计算机实现快速傅里叶变换用计算机实现快速傅里叶变换(3)加权系数加权系数 的计算的计算 的计算主要是确定的计算主要是确定p值。值。p值可用下述方法求得:值可用下述方法求得:把把k值写成值写成r位的二进制数位的二进制数(k序列的序号数,序列的序号数,r是迭代次数是迭代次数);把这个二进制数右移把这个二进制数右移r-1位,并把左边的空位补零位,并把左边的空位补零(结果仍为结果仍为r位位)把这个右移后的二进制数倒序把这个右移后的二进制数倒序把倒序
16、后的二进制数转换维十进制数就得到把倒序后的二进制数转换维十进制数就得到p值。值。第33页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 4.用计算机实现快速傅里叶变换用计算机实现快速傅里叶变换(4)重新排序重新排序 由蝶形流程图可见,如果序列由蝶形流程图可见,如果序列x(n)是按顺序排列的,经过蝶式是按顺序排列的,经过蝶式运算后,其变换序列是倒序排列的;反之,如果运算后,其变换序列是倒序排列的;反之,如果x(n)是倒序排列是倒序排列的,那么输出就是顺序排列的。因此,为了便于输出的,那么输出就是顺序排列的。因此,为了便于输出使用,最好加入倒序程序以保证使用,最好加入倒序程序以
17、保证x(n)与它的变换系数与它的变换系数X(m)的对应关的对应关系。系。第34页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 5.如何提高如何提高FFT的速度?的速度?(1)减少乘法次数;(2)基4、基8算法;(3)实数FFT;(4)硬件实现(DSP芯片,FFT集成块)第35页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 6.FFT举例举例w2w1w3-1-1-1-1F(u)-1-1-1-1w2w2-1-1-1-11001-1001111-1-1-1202000001001f(x)其中:第36页,共129页,编辑于2022年,星期三5.1 傅立叶变换
18、傅立叶变换 6.FFT举例举例21+j01-j000010101000-101111-11-1-j-1 j第37页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 6.FFT举例举例F(u)=2,0,幅度谱:幅度谱图:0121234567uF(u)第38页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 5.1.2 二维傅立叶变换二维傅立叶变换1.二维连续函数傅立叶变换(二维连续函数傅立叶变换(2D FT)定义定义:若若f(x,y)是连续图象函数是连续图象函数正变换正变换:反变换反变换:变换对变换对:第39页,共129页,编辑于2022年,星期三5.1 傅立
19、叶变换傅立叶变换 2.幅度谱、相位谱、能量谱幅度谱、相位谱、能量谱一般一般F(u,v)是复函数是复函数,即即:幅度谱幅度谱:相位谱相位谱:能量谱能量谱:第40页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 二维连续傅立叶变换举例:二维连续傅立叶变换举例:函数函数 f(x,y)=A,0,其他其他 求求F(u,v)。f(x,y)xy0XYAXY(0,0)图象屏幕显示图象屏幕显示第41页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 ;代入函数;分离变量;查积分表;欧拉公式第42页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 两个SIN
20、C函数的乘积幅度谱:u=0,v方向v=0,u方向谱图(截面图)第43页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 幅度谱:幅度谱的屏幕显示:u=0,v方向第44页,共129页,编辑于2022年,星期三Fourier 变换示意图变换示意图第45页,共129页,编辑于2022年,星期三Fourier 变换示意图第46页,共129页,编辑于2022年,星期三Fourier变换的频率特性 第47页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 5.1.2 二维傅立叶变换二维傅立叶变换2.二维离散傅立叶变换(二维离散傅立叶变换(2D DFT)定义定义:若若f(
21、x,y)是离散图象函数是离散图象函数正变换正变换:反变换反变换:第48页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 二维离散傅立叶变换的矩阵表示二维离散傅立叶变换的矩阵表示令:正变换:反变换:(忽略1/N)第49页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 根据可分离性:令:忽略1/N第50页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 FT:IFT:(忽略1/N)第51页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 5.1.3 二维傅立叶变换的性质二维傅立叶变换的性质1.可分离性可分离性正变换第52页
22、,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 同样,反变换也具有可分离性同样,反变换也具有可分离性第53页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 利用二维傅立叶变换的可分离性,可将二维利用二维傅立叶变换的可分离性,可将二维DFT转化成一维转化成一维DFT计计算。即,先在算。即,先在x(或(或y)方向进行一维)方向进行一维DFT,再在,再在y(或(或x)方向进行)方向进行一维一维DFT:第一步:第二步:第54页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 二维离散傅立叶变换过程图示:二维离散傅立叶变换过程图示:第一步:第二步
23、:f(x,y)=F(u,y)=在在x方向逐行进行一维方向逐行进行一维FT在在y方向逐列进行一维方向逐列进行一维FT1/NF(u,v)=第55页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 二维离散傅立叶变换举例二维离散傅立叶变换举例例例1:x方向FFT11001111201-j1+j21-j01+j-1-1-1-1w1y方向FFT1/4第56页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 二维离散傅立叶变换举例二维离散傅立叶变换举例1111220040004000-1-1-1-1w1第57页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立
24、叶变换 二维离散傅立叶变换举例二维离散傅立叶变换举例例例2:二维傅立叶反变换:二维傅立叶反变换用下式求反变换,与正变换使用同一流程:第58页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 二维离散傅立叶变换举例二维离散傅立叶变换举例例例2:二维傅立叶反变换:二维傅立叶反变换第59页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 二维离散傅立叶变换举例二维离散傅立叶变换举例例例2:二维傅立叶反变换:二维傅立叶反变换逐列FFT逐行FFT1/4逐像素求共轭第60页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 5.1.3 二维傅立叶变换的性
25、质二维傅立叶变换的性质2.平移性平移性FT则:则:相当于相当于F(u,v)的坐标原点)的坐标原点移到(移到(u0,v0)点)点第61页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 即:移中性移中性同理:第62页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 移中性移中性移中性的用途:移中性的用途:图象作傅立叶变换时,若采用以下公式变换,则变换后主要能量图象作傅立叶变换时,若采用以下公式变换,则变换后主要能量(低频分量)集中在频率平面的中心。(低频分量)集中在频率平面的中心。第63页,共129页,编辑于2022年,星期三5.1 傅立叶变换傅立叶变换 移中
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第五 图象 正交 变换 PPT 讲稿
![提示](https://www.taowenge.com/images/bang_tan.gif)
限制150内