多维度变换
在数学分析及应用中,“多维度变换”是用来分析讯号的二维或是多维的频率成分。
多维度傅立叶变换
其中一个常用的多维度变换就是傅立叶变换,是将一个讯号的表示式从时域/空域转换到频域。[1]离散域的多维度傅立叶变换可表示成下列式子:
其中F代表多维度傅立叶变换,m代表维度。将f定义成多维度的离散域讯号,则逆多维度傅立叶变换为:
连续域的多维度傅立叶变换可表示成下列式子:[1]
快速傅立叶变换(FFT)是一种用来计算离散傅立叶变换(DFT)和其逆变换的快速算法,快速傅立叶变换所得到的结果跟按照定义去算离散傅立叶变换的结果是一样的,但唯一的差别是快速傅立叶变换的速度快很多。(在舍入误差的存在下,很多快速傅立叶变换还比直接照定义算还更精准。)有很多种快速傅立叶变换,他们包含很广泛的数学运算,从简单的复数运算到数论和群论,详情可以看快速傅立叶变换。
多维度的离散傅立叶变换是离散域傅立叶变换的简单版本,其方法是在均匀间隔下的样本频率去估计其值。[2] N1 × N2 × ... NM离散傅立叶变换如下式:
其中0 ≤ Ki ≤ Ni − 1, i = 1, 2, ..., m。
多维度离散傅里叶变换对应的逆变换是:
其中0 ≤ n1, n2, ... , nm ≤ N(1, 2, ... , m) – 1。
多维度离散余弦变换
离散余弦变换被广泛的应用,像是资料压缩、特征萃取、影像重建等等。多维度离散余弦变换为:
其中 ki = 0, 1, ..., Ni − 1, i = 1, 2, ..., r.
应用
离散傅立叶变换和离散余弦变换常常被使用在讯号处理[3]和影像处理,也常被用来当作解偏微分方程式时更有效率的方法。离散傅立叶变换也可用在运算折积或是乘上很大的整数。下列只列出一些例子。
影像处理
![](http://images.weserv.nl/?url=//upload.wikimedia.org/wikipedia/commons/thumb/2/23/Dctjpeg.png/250px-Dctjpeg.png)
离散余弦变换被用在 JPEG 影像压缩、MJPEG、MPEG、DV和 Theora影片压缩上。压缩时使用NxN'格的二维的离散余弦变换(DCT-II)然后再被量化且用熵编码法编码,通常N为8,而DCT-II的运算就用在该格的每一行和每一排,结果会生成8x8的变换系数矩阵,其中(0,0)(左上角)的值是直流分量(频率为0),随着水平或垂直的编号增加,代表水平或垂直的空间频率增加,如右图所示。
在影像处理方面,利用二维的离散余弦变换可以分析并且描述非常规的图形加密方法,像是在二维图像平面中插入非可见的二进位制水印。[4]利用不同的方向,DCT-DWT混杂的转换也可以用来去除超声波影像的噪声。[5]三维的离散余弦变换可以被用来转换在使用水印影像迁入的影片资料或是三维影像资料。[6][7]
频谱分析
当使用离散傅立叶变换来做频谱分析时,{xn}的数列通常代表着从讯号 x(t)中在均匀的时间点做取样所得到的有限集合,这样将连续时间点经取样离散化后,也将原本的傅立叶变换转变成离散时间傅立叶变换(DTFT),通常也因此产生了混叠的失真。为了要最小化这种失真,选择适当的取样频率是重点(详情请看取样定理)。同样的,将一个非常长(或无限)的数列转变成一个容易处理的大小,会因此造成失真(Spectral leakage),选取一个适当的子数列长度是最小化这个问题的关键点。当资料量大于达到理想频率分辨率所需的适量时,标准的作法是使用多个DFT,例如产生频谱图的时候。如果所期望的结果是功率频谱而且有噪声或随机讯号出现在资料内的话,多个DFT的振幅平均值可以用来减少频谱的变异性,Welch方法和Bartlett方法就是这种技术。一般处理这种用来估计有噪声的讯号的功率频谱的方法就称为频谱估计。
其实会造成失真的主要源头就是DFT本身,因为DFT是将DTFT这种连续性的频域做离散取样的结果,可以利用提高DFT的频率分辨率来减缓这问题。
- 这种方法有时候也被认为是零填充,这是一种被用在快速傅立叶变换的一种特别应用。这种因为值为零的取样点而产生的乘法与加法比原本的FFT产生偏移还要没有效率。
- 如上面所言,失真(leakage)的问题对DTFT的频率分辨率造成了限制,因此会对透过提高频率分辨率的效益造成限制。
偏微分方程式
离散傅立叶变换时常被用来解偏微分方程式,其中DFT是被用来近似傅立叶级数,其优点在于将讯号延伸为复数指数函数einx,而这微分方程式的特征函数为:d/dx einx = in einx,因此,微分在傅立叶变换后的表示式下变得很简单,只要乘上i n (但是因为混叠的影响,n的选择不一定是唯一的,为了让这样的方法达到收敛,需要使用三角插值去选择这个值)。一个线性微分方程式且其系数为常数经由离散傅立叶变换后会变换成很容易解的代数等式,将解完等式的结果用逆离散傅立叶变换就能回到原来的时/空域表示式。
用快速傅立叶变换处理影像艺术面的分析
我们必须使用没有损害的方法去得到一些关于艺术稀有的资讯(从HVS的观点是着重于色度法以及空间资讯)。我们可以透过观察色彩变化或是测量表面一制性的变化来了解艺术,因为整个影像是非常大的,所以我们会使用一个双生的余弦窗去撷取影像:[8]
其中N是影像的维度, x, y是从影像中心(0,0)所扩展的座标(从0到N/2),可将空间频率表示成下式:[8]
“FFT”为快速傅立叶变换, f是空间频率。这种基于FFT的成像方法是一种诊断技术,用以确保文化艺术的长寿及稳定。这是一种简单、成本低且可用于博物馆又不影响日常使用的方法,但这种方法没办法定量的量测腐蚀速率。
参见
参考资料
- ^ 1.0 1.1 Smith,W. Handbook of Real-Time Fast Fourier Transforms:Algorithms to Product Testing, Wiley_IEEE Press, edition 1, pages 73–80, 1995
- ^ Dudgeon and Mersereau, Multidimensional Digital Signal Processing,2nd edition,1995
- ^ Tan Xiao, Shao-hai Hu, Yang Xiao. 2-D DFT-DWT Application to Multidimensional Signal Processing. ICSP2006 Proceedings, 2006 IEEE
- ^ Peter KULLAI, Pavol SABAKAI, JozefHUSKAI. Simple Possibilities of 2D DCT Application in Digital Monochrome Image Cryptography. Radioelektronika, 17th International Conference, IEEE, 2007, pp. 1–6
- ^ Xin-ling Wen, Yang Xiao. The 2-D Directional DCT-DWT Hybrid Transform and Its Application in Denoising Ultrasound Image. Signal Processing. ICSP 2008. 9th International Conference, Page(s): 946–949
- ^ Jinwei Wang, Shiguo Lian, Zhongxuan Liu, Zhen Ren, Yuewei Dai, Haila Wang. Image Watermarking Scheme Based on 3-D DCT.Industrial Electronics and Applications, 2006 1ST IEEE Conference, pp. 1–6
- ^ Jin Li, Moncef Gabbouj, Jarmo Takala, Hexin Chen. Direct 3-D DCT-to-DCT Resizing Algorithm for Video Coding. Image and Signal Processing and Analysis, 2009. ISPA 2009. Proceedings of 6th International Symposium pp. 105–110
- ^ 8.0 8.1 Angelini, E., Grassin, S. ; Piantanida, M. ; Corbellini, S. ; Ferraris, F. ; Neri, A. ; Parvis, M. FFT-based imaging processing for cultural heritage monitoring Instrumentation and Measurement Technology Conference (I2MTC), 2010 IEEE