第3章图像变换.ppt
《第3章图像变换.ppt》由会员分享,可在线阅读,更多相关《第3章图像变换.ppt(90页珍藏版)》请在三一文库上搜索。
1、第3章 图像变换,傅立叶变换(FFT) 离散余弦变换(DCT) K-L变换 小波变换,问题的提出,为达到某种目的将原始图象变换映射到另一个空间上,使得图象的某些特征得以突出,以便于后面的处理和识别。 一般变换后的图象,大部分能量都分布于低频谱段,这对以后图象的压缩、传输都比较有利。使得运算次数减少,节省时间,傅立叶变换FFT,线性系统,对于一般线性系统,往往是用时间作为参数来描述的,表示为一维(t)系统。在图像处理中是用空间作为参数来描述的,通常表示为二维(x,y)系统。输入函数f(x,y)表示原始图像,输出函数g(x,y)表示经处理后的图像,线性系统可看作是一种映射,它反映了各种线性的图像处
2、理方法,其输入和输出的关系表示为:,二维线性系统,叠加原理,点源和狄拉克函数,一幅图像可以看成由无穷多极小的象素组成,每一个象素都可以看作为一个点源,因此,一幅图像也可以看成由无穷多点源所组成 数学上,点源可以用狄拉克函数来表示,二维函数定义为 满足,狄拉克函数性质,(1) 函数为偶函数,即,(2)位移性,或用卷积符号*表示为,(3)可分性,(4)乘积性,(5)筛选性,当且仅当0时,(6)指数函数,卷积,假设f(x)(x=0,1,A-1)以及g(x)(x=0,1,C-1)是两个有限离散函数,其线性卷积为,对于图像二维函数的卷积,则,相关,2个函数的相关定义为,其中f*(i)为f(i)的复共轭,
3、一维连续傅立叶变换,定义及基本概念,设f(x)为x的函数,如果f(x)满足下面的狄里赫莱条件:,(1)具有有限个间断点(2)具有有限个极值点(3)绝对可积,则有下列式成立:,x为时域变量,u为频域变量,(1),(2),如果令w2u,则有,称为傅立叶变换对,函数f(x)的傅立叶变换一般是一个复量,它可以用下式表示:,实例,傅立叶变换,其相位为,傅立叶谱,周期函数的傅立叶变换,傅立叶级数来表示,结论,1、只要满足一定条件,连续函数就可以进行傅立叶变换,实际上这个条件在工程应用中很容易满足。 2、连续非周期函数的傅立叶谱是连续的非周期函数,连续的周期函数的傅立叶谱是离散的非周期函数。,二维连续傅立叶
4、变换,如果二维函数f(x,y)满足狄里赫莱条件,则将有下面的傅立叶变换对存在:,与一维傅立叶变换类似,二维傅立叶变换的傅立叶谱和相位谱为:,例:求如图所示的函数的傅立叶谱,f(x,y)函数,其傅立叶变换为:,其傅立叶谱为:,离散傅立叶变换,由于实际问题的时间或空间函数的区间是有限的,或者是频谱有截止频率 离散傅立叶变换(Discrete Fourier Transform简称DFT)在数字信号处理和数字图像处理中应用十分广泛,它建立了离散时域和离散频域之间的联系,一维离散傅立叶变换,二维离散傅立叶变换,由于图像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度。 如:大面积的沙
5、漠在图像中是一片灰度变化缓慢的区域,对应的频率值很低;而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域,对应的频率值较高。 傅立叶变换在实际中的物理意义,设f是一个能量有限的模拟信号,则其傅立叶变换就表示f的谱。 从纯粹的数学意义上看,傅立叶变换是将一个函数转换为一系列周期函数来处理的。 从物理效果看,傅立叶变换是将图像从空间域转换到频率域,其逆变换是将图像从频率域转换到空间域。 换句话说,傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数,傅立叶逆变换是将图像的频率分布函数变换为灰度分布函数。,一个MN大小的二维函数f(x,y),其离散傅立叶变换对为,在数字图
6、像处理中,图像一般取样为方形矩阵,即NN,则其傅立叶变换及其逆变换为,(1)可分性,从上式可以看出,一个二维傅立叶变换可用二次一维傅立叶变换来实现,傅立叶变换的性质,其意义:一个二维傅立叶变换或反变换都可以分解为二步进行,其中每一步都是一个一维傅立叶变换或反变换,行变换,列变换,二维傅立叶变换分离成两个一维变换,(2)平移性,在空域中,图像原点平移到(x0,y0)时,其对应的频谱F(u,v)要乘上一个负的指数项,在频域中,原点平移到(u0,v0)时,其对应的f(x,y)要乘上一个正的指数项,也就是说,当空域中f(x,y)产生移动时,在频域中只发生相移,而傅立叶变换的幅值不变,反之,当频域中F(
7、u,v)产生移动时,相应的f(x,y)在空域中也只发生相移,而幅值不变,平移,(3)周期性和共轭对程性,周期性可表示为,该式表明,如果空间域函数f(x,y)旋转0角度后,相应的傅立叶变换F(u,v)在频域中也旋转同一0角,反之,F(u,v)在频域中旋转0角,其反变换f(x,y)在空间域中也旋转0角,(5)分配性和比例性,傅立叶变换的分配性表明,傅立叶变换和反变换对于加法可以分配,而对乘法则不行,即,傅立叶变换的比例性表明,对于二个标量a和b,有,在空间比例尺度的展宽,相应于频域中比例尺度的压缩,其幅值也减少为原来的,(6)平均值性质,定义二维离散函数的平均值为,将u=v=0代入二维离散傅立叶公
8、式,可得,比较上面两式,可看出,(7)卷积定理,卷积定理和相关定理都是研究两个函数的傅立叶变换之间的关系,这构成了空间域和频域之间的基本关系,对于两个二维连续函数f(x,y)和g(x,y)的卷积定义为,指出傅立叶变换一个主要好处:与其在一个域中作不直观的,和难懂的卷积,不如在另外一个域中作乘法,可以达到相同的效果,它表明两个二维连续函数在空间域中的卷积可用求其相应的两个傅立叶变换乘积的逆变换而得。 反之,在频域中的卷积可用空间域中乘积的傅立叶变换而得,应用卷积定理明显的好处是避免了直接计算卷积的麻烦,它只需要先算出各自的频谱,然后相乘,再求其逆变换,即可得到卷积。,离散傅立叶变换和逆变换都是周
9、期函数,为了防止卷积后产生交叠误差,需要对离散的二维函数的定义域加以扩展。 利用一维函数来说明函数定义域的扩展,设f(x)(x=0,1. A-1)和g(x)(x=0,1,.C-1)是两个有限的离散函数,也就是说f(x)定义域为0xA-1,g(x)的定义域为0xC-1。则其线性卷积为,如果利用卷积定理来计算该卷积,则相当于把f(x)和g(x)分别以A和C为周期进行周期延拓,因此,将上式改写为(一个周期) 求得结果将不是需要的z(x),而是周期循环的,为了避免循环卷积,可以对原被卷积函数补零。由于卷积结果的长度为N=A+C-1,因此,可以把两个被卷积的函数的长度扩展到N,并在原函数定义区间外的部分
10、补零,即取,同样地,二维离散卷积计算时,也必须对被卷积函数进行延拓和补零,如果被卷积函数f(x,y)和g(x,y)地大小分别为AB和CD,则延拓后地函数为:,(8)相关定理,对于二维连续函数f(x,y)和g(x,y)的相关定义为,相关定理可表示为,在离散情况下,与离散卷积一样,需要用增补零的方法扩充f(x,y)和g(x,y)为fe(x,y)和ge(x,y),并根据卷积定理相类似的方法来选取M和N,以避免在相关函数周期内产生交叠误差。,快速傅立叶变换(FFT),离散傅立叶变换已成为数字信号处理的重要工具,然而,它的计算量达,运算时间长,在某种程度上却限制了它的使用范围。 快速算法大大提高了运算速
11、度,在某些应用场合已能作实时处理,并且应用在控制系统中 快速傅立叶变换不是一种新的变换,它是离散傅立叶变换的一种算法,它是在分析离散傅立叶变换中的多余运算的基础上,进而消除这些重复工作的思想指导下得到的,其原理:,对于一个有限长序列f(x)(0x N-1),它的傅立叶变换由下式表示:,令,将正变换(1)展开得到:,矩阵表示为:,从上式可以看出,要得到每一个频率分量,需进行N次乘法和N-1次加法运算。要完成整个变换需要N2次乘法和N(N-1)次加法运算。当序列较长时,必然要花费大量的时间,观察上面的系数矩阵,发现Wmn是以N为周期的,即,当N=8时,其周期性如图所示,由于,所以当N=8时,可得:
12、,可见,离散傅立叶变换中的乘法运算有许多重复内容。1965年库利-图基提出原始的N点序列依次分解成一系列短序列,然后,求出这些短序列的离散傅立叶变换,以此来减少乘法运算,例如,设:,由此,离散傅立叶变换可写成下面的形式:,因为:,所以:,F1(u)和F2(u)分别是f1(x)和f2(x)的N/2点的傅立叶变换,由于F1(u)和F2(u)均是以N/2为周期,所以,这说明当mN/2时,上式也是成立的,因此下式成立:,由上面的分析可见,一个N点的离散傅立叶变换可由两个N/2点的傅立叶变换得到,离散傅立叶变换的计算时间主要由乘法决定,分解后所需的乘法次数大为减少。第一项为(N/2)2次,第二项为(N/
13、2)2N次,总共为2(N/2)2N次运算即可完成,而原来需要N2次运算,可见分解后的乘法计算次数减少了近一半。,当N为2的整数幂时,则上式中的F1(u)和F2(u)还可以再分成两个更短的序列,因此计算时间会更短。 由此可见,利用WNm的周期性和分解运算,从而减少乘法运算次数是实现快速运算的关键,实例,离散余弦变换,离散余弦变换(Discrete Cosine Transform-简称DCT)是傅里叶变换的一种特殊情况。在傅里叶级数展开式中,被展开的函数是实偶函数时,其傅里叶级数中只包含余弦项,称之为余弦变换 DCT计算复杂性适中,又具有可分离特性,还有快速算法,所以被广泛地用在图象数据压缩编码
14、算法中,如JPEG、MPEG-1、MPEG-2及H.261等压缩编码国际标准都采用了离散余弦变换编码算法 其变换核是为实数的余弦函数,因而DCT的计算速度比DFT快得多 对于具有一阶马尔可夫过程的随机信号,DCT是K-L变换的最好近似,已被广泛应用到图像压缩编码、语音信号处理等领域,一维离散余弦变换,一维离散余弦反变换,二维离散余弦变换,二维离散反余弦变换,如果令N=4,由一维解析式定义可得如下展开式:,写成矩阵形式:,F(u)=Af(x),同理,可得到反变换展开形式:,写成矩阵形式:,f(x)=ATF(u),二维离散余弦变换为:,F(u,v)=Af(x,y)AT f(x,y)=ATF(u,v
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 图像 变换
链接地址:https://www.31doc.com/p-2526522.html