分数傅里叶变换

1、分数傅里叶变换分数傅里叶定义:分数傅里叶变换 的物理意义即做傅里叶变换.次，其中不一定要为整数（比傅里叶变换更加广泛）；通过分数傅里叶变换 之后，图像或信号便会同时拥有时域与频域两者的特征。1.1 （维基百科） 第一种定义：第二种定义:X申(it)=1 jcot(/2?r厂1沖i.J1.2从数学上分数傅立叶变换定义了积分形式:A. C. Me bride和F. H Ken在1987不给出了 V Narnia s的分数傅里叶变换 的枳井形式讥 具体地说，对信号空间用）中的任何信号/（/），它的分数 傅里叶变换（时卩）可以写成积分形式0/）2） = / *（卩；叮）由（2-5）| i(cota-2

2、v/csc(a 4- r cot ( al) | p 2/i(2-6)其积分核是r町）exp斤（p;以）=*公式中各记号的含义是a =(2-7)式中小整数.p分数傅里叶变换的级次可取任何实数.Wigner分布函数相空间定义的分数傅立叶变换A.W.Lohmann在1993年利用傅里叶变换相当于在 Wigner分布函数相空间中角度为 n 12的旋转这一性质，说明分数傅里叶变换在 Wigner分布函数空间中相当于角度是 pn 12的 旋转，这里，p是分数傅里叶变换的级次。伙I此* AAVLohmaiin 义杲次是p的分数傅里叶变换（产丁）（列为It中，矩阵眉/）疑时 频相平面x-r l.ffi度为（

3、pn/2）的旋转矩阵砒）=Zcosu JwinI-jh mnd4D（5）对恳也）进行正立归一代得到叫仏）可比采用了两种不同正交映 射的方法，并根捌采用的正交化方法不同而分别命名OPA方法和GSA方 法.（6）再由暫甘）构造U彳讪I |产 J 叫T；VmQg（2-41）7）构造少尸3-2A- A rnod? I(2-42)（8）最后由和U得到离散分数傅立叶变换矩阵尸户=UD（,Ur（2-43）离散分数傅立叶算子尸尸无论p为任恿实数，均满足特征方程（2-27）.井 且满足阶数可加性尸尸戸=严。信号（左）的离散分数傅立叶变换兀仗）通过如I：公式计算：（2-44）与连绩的分数傅立叶变换相似，信号”約也

4、可通过我逆变换恢复： x（k=Kf，X（2-45）基于正交投影的离散化算法 对连续分数傅立叶变换的特征函数进行离散化近似和正交投影， 得到一组与 Hermite-Gaussian函数形状相似的离散傅立叶变换矩阵的正交化离散Hermite特征向量。然后，仿照连续分数傅立叶变换的核函数谱分解表达式，构造了离散分数傅立叶变换矩阵。此算法适合用来计算连续的分数傅里叶变换。基于chirp分解的离散化算法，将分数傅立叶变换分解为信号的卷积形式后，直接离散化，利用FFT来计算分数傅立叶变换。图像的分数傅里叶变换对图像进行分数傅里叶变换分析的目的是确定图像经过分数傅里叶变换后的特性表现， 主要包含分数傅里叶变

5、换对图像能量分布和频率分布影响两方面的内容。其中能量分布表现分数傅里叶变换图像的能量聚积性与分数变换阶数的关系，能量聚集性强烈地依赖于其接近于傅里叶变换的程度；频率分布表现在分数傅里叶变换的相位函数包含了图像的纹理频率信 息，变换阶数不同，相位函数所含的图像边缘高频信息也不相同。图像经过某种二维离散变换之后的能量分布体现了图像的变换特征。图像分数傅里叶变换域的能量分布特点是：能量向中心区域聚集性。(1) 当分数阶次p由小变大时，由相位函数恢复的图像呈现出图像边缘轮廓变得越来越清晰，这类似于原始图像经历了不同截止频率的高通滤波器。当p较小时对应于截止频率较低的高通滤波器，低频成份浮现出来，图像边

6、缘模糊；当p较大时，对应于截止频率较高的高通滤波器，大部分低频成份被滤掉，图像边缘比较清晰，FRFT逐渐向FT退化。(2) 当变换阶数p由小变大时，仅由幅度函数恢复的图像越来越接近原图像的背景，这类似于原图像经历了不同截止频率的低通滤波器。p较小时，对应于截止频率较高的低通滤波器，高频分量残留较多，能清晰看到原图像的轮廓；p较大时，对应于截止频率较低的低通滤波器，大部分高频分量被滤出只显现原图像背景。(3) 当变换阶数p为其它值时，由FRFT相位函数和幅度函数所恢复的图像既包含了原 图像的背景信息又包含了原图像的纹理频率信息。由此可以推论这类似于原图像经历了FRFT的时频滤波，也即将时频平面旋

7、转某一角度后再进行滤波。假如频域滤波器截止频率 和带宽固定，当旋转的角度不同(阶数不同)时，时间轴和频率轴上的投影不同，所以频域滤波器输出的频率成分也不同。 这表现在恢复的图像上即为相位函数和幅度函数包含的频率成 份随阶数而变化。当变换阶数较小时，由图像的FRFT的幅度函数和相位函数恢复的图像都显示出很强的 图像信息，体现出了较强的空域特性；当变换阶数逐渐增大时，图像的 FRFT的幅度函数恢 复的图像所包含的原图像的空域特征逐渐减弱直至消失，相位函数恢复的图像包含原图像的边缘纹理特征逐渐增强。当FRFT的变换阶数增大到一定程度时，其幅度和相位特征越来越接近FT域即频域特征。这些结论有力地体现了

8、FRFT域的空-频双域特征。采用分数傅立叶变换的图像边缘提取方法对图像作连续小级数的分数傅立叶变换，相当于对图像作连续的微小变换，当分数级次很小时，肉眼几乎看不出与原图的区别， 当级次略有增加，图像边缘与原图有了明显的区别， 当继续缓慢增加分数级次时，图像与原图明显不同。 通过分析可以看出， 图像中对比度低的区域随级次变化缓慢，对比度高的区域(即图像边缘)随级数变化快。由此，取不同级数的 分数傅立叶变换后的图像减去原图像，即可得到图像的边缘。 不同的分数级数对应不同的形变，选取不同级数变换后的图像相减，即可提取不同尺度的边缘。图像分数阶Fourier变换的幅度和相位信息假设 F(k,h )是二

9、维图像 f(x, y)的二维 Fourier 变换 F (k,h )= FT2D f (x, y )( 8)我们可以把F k,h分解成幅度部分和相位部分，即(9)F (k,h )=F (k,h j_|P(k,h ) = A(k,h )P(k,h )其中A k,h二F k,h为幅度函数，P k,h二F k,hjAk,h为相位函数,结论：1. 当变换阶数P由小变大时，仅由相位函数恢复的图像， 显现原图像的边缘越来越清晰，这类似于原图像经历了不同截止频率的高通滤波器。p较小时(0.01)对应于截止频率较低的高通滤波器，低频成份浮现出来，使提取的边缘模糊，如图4(b)所示；p较大时(0.8)，对应于截

10、止频率较高的高通滤波器，大部分低频成份被滤出，提取的边缘较清晰，如图5(d)所示，此时FRFT基本退化为FT。2. 同理,当变换阶数P由小变大时，仅由幅度函数恢复的图像越来越接近原图像的背景，这类似于原图像经历了不同截止频率的低通滤波器。p较小时(0.01)，对应于截止频率较高的低通滤波器，高频分量残留较多，还能清晰看到原图像的轮廓，如图4(a)所示；p较大时(0.8)，对应于截止频率较低的低通滤波器，大部分高频分量被滤出，此时仅显现原图像的背景，如图6(c)所示。3. 当变换阶数P为其它任意值时，由FEFT相位函数和幅度函数所恢复的图像既包含了原图像的背景又包含了原图像的纹理，如图4(b)、

11、图4( c)、图5( a)、图5( b)所示。针对这种情况，我们可以推论这类似于原图像经历了FRFT的时频滤波，也即将时频平面旋转某一角度后再进行滤波。假如频域滤波器截止频率和带宽固定，当旋转的角度不同(阶数不同)时，时间轴和频率轴上的投影不同，所以频域滤波器输出的频率成分也不同。表现在恢复的图像上即为相位函数和幅度函数包含的频率成份随阶数而变。分数傅里叶变换(采用 chirp信号的方法)a从0.1-1过程中，相应的分数傅里叶变换。 从图中可以发现，两个区域信息的相互转换。相应的幅值信息(a从0.1-1)相应的相位信息（即：f/fl，其中f为相应分数傅里叶变化后，输出的信号 ）相应的幅值信息进行逆变换根据相位信息进行相应的逆变换，可以发现图像的边缘信息，消除了图像的背景信息。总结通过分数傅里叶变换和其逆变换可以找到一些图像的主体边缘信息， 但是有些图像边缘 的采集并不是很好， 对于背景信息的剔除不是很好。 对于区域的选择不是很准确， 它主要是 根据灰度的变换率，而不是根据像素去识别边界。傅里叶变换主要是进行像素变化率的识别， 即高频和低频的识别和分类。 因此容易受到 其他