对傅立叶变换后图像空间域与频率域中垂直现象的研究.

1、对傅立叶变换后图像空间域与频率域中垂直现象的研究对傅立叶变换后图像空间域与频率域中垂直现象的研究【摘要】本文就图像二维傅立叶变换常令人感到困惑的问题进行深入的讨论，并对傅 立叶变换后原图和频谱图出现的垂直现象进行分析，同时给出数学证明。【关键词】傅立叶变换；频谱图；垂直现象；图像1引言傅立叶变换是线性系统分析的一个有力工具，是信号处理中最重要、应用最广泛的变 换，但遗憾的是很多人可能还是不习惯在频域中思考问题，尤其是面对图像时，空间域、 频率域都是二维的，更是对空域、频域的对应关系不甚了解。如何理解傅立叶变换后的频 谱图，为何对不同图像进行傅立叶变换后产生的频谱图往往会出现水平和垂直的“亮线”

2、 而对一些规则图像进行傅立叶变换后得到的频谱图往往会在与原图垂直的方向上出现较 亮的点或线(如垂直条纹图像，傅立叶变换后会在水平方向上出现一些较亮的点)。这些问题困扰着每个初学者，更让初学者困惑的是，几乎市面上所有有关数字图像处理的书籍都 没有给出详细的解释。下面就围绕傅立叶变换以及图像处理中一些经常让人困惑的问题进 行深入的讨论。2 一维傅立叶变换中的问题讨论17一维傅立叶变换是二维变换的基础，下面就先对一些相关的、基本的但又很重要的概 念进行讨论。(1)实信号的傅立叶变换就信号处理来说，大家所关心的都是实信号，所以单独对它进行讨论，可以简化工作。 下面是大家所熟知的傅立叶变换公式F()f(

3、t)exp( 2t)dt从定义式不难推出f (t)ex p(2t)dt所以可得出结论：F()F ()。这说明实函数的傅立叶变换是实部为偶函数，虚 部为奇函数，也就是说：求某一实函数的傅立叶变换时，它的幅度谱总是关于原点对称的, 而相位谱左右两边只是差一个负号，即左右互为复共轭。由此可见，就实际应用来讲，无 论哪边的频谱都是完备的，并且负频率本身也不具有什么意义，但是当用更为通用的数学 方法去对物理过程建模时，保留负频率部分会使分析更加容易。(2)正弦分量的分解大家知道，任何满足狄里赫利条件的信号都可以通过傅立叶变换表征为一组正(余)弦信号的和或积分，而由图1可以从频域得出这一结论：由于任何实函

4、数的傅立叶变换都是偶函数，所以对于任意一个实信号f(t)，将它的傅立叶变换F()在频率域进行抽样时，都能得到无数的抽样脉冲对 (见图2)，也即将这些脉 冲对累加起来就得到了 F()，而当脉冲宽度趋于零时，每个脉冲对正好是某个频率的余弦 信号的傅立叶变换，这从频率域的角度验证了任何一个实信号都可以看作是由若干个正 (余)弦信号以及相应的幅度所组成。)X图2实信号f (t)及傅立叶变换F ()(2)信号在时域通过以上的讨论可得出两个简单的结论：(1)实信号的频谱是对称的；和频域中是相互对应的，总是能把实信号看作是由若干不同频率、振幅的正弦波组成。3.图像的二维傅立叶变换3.1二维离散傅立叶变换的定

5、义图像经数字化处理后，可以用二维离散信号f(m,n)表示。对于二维离散信号 f (m, n)|m 0,1,L,M1;n0,1,L , N 1，其离散傅立叶变换定义为：M 1N 1F(u,v)一 Z C / mu nV/ 八f (m, n)exp( j2 ( )( 1)M N式中u 0,1丄,M1;v0,1丄,Nf(m, n)1，称为空间频率。反变换定义为1muF(u,v)exp(j2 (0n 0nv M N )式中m 0,1丄,M1;n0,1丄,N在图像处理时,变换和反变换式：一般选取图像块为N N的方阵，即取MN，这时二维离散傅立叶F(u,v)1 N 1N 1mu書 f (m, n)exp(

6、 j2 ( N m 0n 0Nnv亓)f (m, n) N 1 N 11匚 z、z s mu丁F(u,v)exp(j2 (N m 0 n 0N(4)在(3)(4)两式中，u,v,m, n 0,1丄,N 1。本文都是选取NN的图像进行讨论的。3.2频谱图的理解2要直接表示结果就必须由(3)式可知图像经傅立叶变换后，往往得到的是复数形式。用到两幅图像：一幅表示实部，一幅表示虚部。这样表示十分不方便，同时也没有得到有 用的信息，因此引入变换结果的模作为值在频谱图中表示出来，以灰度的明暗代表模的大 小。作为典型的二维信号，图像的频率相应地也是二维的。其分别对应着图像的像素值在 两个相互垂直的方向上变化

7、的情况(如图3.2.1 )。根据对一维离散信号频谱的分析结果可知，频谱在u 0,N/2, N 1三点处的频率分别为f0 0, fN2 fc,fN 1 0 fc为图像信号的最 高截止频率。将其推广到二维，则在图像频谱图中所有沿u方向的频率值变化情况与一维相同，也有f。 0, fN/2fc, fN 1 0。同理，在v方向上也有相同的结果。因此，在频谱图四角(0,0), (N 1,0), (0, N 1), (N 1,N1)处沿u和v方向的频率分量均为0,在频谱图中心点(N/2, N/2)处沿u和v方向的频率分量均为最大值fc。N-1AB1 1 /c(N/2,N/2)CDN/2N/2(N/2,N-1

8、)(N-1,N/2)(N-1,N-1)N-1v图3.2.1图像频谱示意图由于图像中的大部分能量集中在低频分量上，因此频谱图中四角部分的幅度值较大。 然而，在实际的图像频谱分析过程中，由于低频分两区域较小，并且分散在四角，因此不 利于对其进行分析。此时，可以根据图像频谱的周期性和共轭对称性对频谱图坐标进行移 位，将所有低频分量集中在频谱图中心，同时高频分量分散在四周。在具体实现频谱的移所以。位过程中，将图3.2.1中的区域A和区域D对换位置，将区域B和区域C对换位置即可。 进行移位之后大大增加了图像频谱的可读性，如图3.2.2XIf(a)原图像(b)移位前的频谱图3.2.2简单矩形图像坐标移动前

9、后的频谱图(c)移位后的频谱3.3频谱图的垂直现象及单条直线段垂直现象对图像进行傅立叶变换后所得到的移位后的频谱图往往会在垂直方向和水平方向上出 现两条“亮线”(如图3.2.2)，这个现象常常困扰初学者，为什么会出现这两条亮线，什 么时候会出现这两条亮线？要很完整的回答这些问题很不容易，因为输入图像可以有很多 种，但产生的频谱图大多都有这两条亮线。因此应该选取一些特殊的图像进行研究，接下 来就选取：单条直线段，fringe patterns图像进行讨论。立叶变换并后未移位的频谱图：如图上方的为原图，下方为傅lU(a) sin( x+y)(c) x + y = N-1zia6-6,2(b) x

10、= a图3.3 一些特殊图像及其傅立叶变换后的频谱图图3.3(a)是fringe patternssin(x+y)的图像，可以看到图像变换后的频谱图中出现两个点，这两点分别与原点及(N-1,N-1)的连线是沿着u = V的方向。而原图像在x+y=b （b为0到2N-2的正整数）的灰度值是相等的，这与u = v方向正好是垂直的。图（b） 中的垂直现象就更明显了，原图为一水平线，变换后的频谱图中在 V=0上出现一条亮线， 这与原图的方向正好垂直。图（C）为直线段x + y = N-1，变换后的频谱图为直线段u = v， 同样出现了垂直现象。上述对一些比较特殊的图像进行了傅立叶变换后中都出现了原图像

11、与频谱图存在一定 垂直关系的现象。对这样的垂直现象的研究能较好的理解频谱图和二维傅立叶变换，接下 来就给出垂直现象的数学证明和物理意义的解释。4.垂直现象的数学证明本节给出只由一条直线段组成的图像的数学证明，在4.1中给出对图像内任意一条水平的直线段的证明，在4.2节中给出任意一条斜线段的证明。N 1N 1f (x, y)sin(2x 0 y 012 2L L L L L (1)因为频谱图中的灰度代表傅立叶变换的模值，因此证明变换前后出现垂直现象，实际 上就是求傅立叶变换后模将会在什么位置或者说什么方向上出现最大值。即讨论求模公 式：2N 1N 1ux vyF(u,v)f (x, y)cos(

12、2 (- *)x 0 y 0NN的最大值将出现在u,v的什么方向上。4.1水平直线段垂直现象的证明设输入图像为一个（N 1,N1）的二维矩阵，图像内任意一条水平直线段定义如下:设图像内任意一条强度为 1宽度为一个像素的水平线为：f(x,y) 1x a a 0,1,2L1 ; y 0,1,2L N且有1f(x, y) 0x带入（1）式可得：2 N1cos（2y 0F(u,v)ua2vf)N 1si n(2y 0ua(uaN 1c,2 ua(coscosy0N,2 uaN1 2 vy、2 (cos) ( cos)N y 0 N(2 ua、2r 1 . 2 vy、2 (cos) ( sin)N y

13、0 NN 1 c2 . 2 vy 2)(sin) y 0 NN 1 c.2 vy、2(cos,y 0 N2 vyN22vy).2 ua . 2 sinsin N N(i 2 ua、2(N1 i 2(sin) ( sin N 1y 0(sin Nvy)2N y0 N,.2 uaN 12 vy、2(sin) ( cos)N y 0 N2 ua 2 vycos22 ua . 2 vy、cosSin)N N对傅立叶变换后图像空间域与频率域中垂直现象的研究29由(2)式可以看出，模的值只可能出现在 v的一列或者多列，也就是说在变换后的频谱图里亮线可能出现在垂直的方向上，这就很好的证明了图3.3(b)中出

14、现的垂直现象。但是(2)式并不能说明模值只出现在v 0这一列，因此还要对(2)式进行进一步的讨论。如果把(2)式中的连加看成是对积分的取样，并把N 1扩充到N，则可以得到如下的(3)式，并对其计算可得以下的等式：N 2 vycos()dy0 NN 2 vycos()dy0 N yN 2 vycos()dy0 NN 2 vya cos(0NN 2 vNcos(0NN 2 vt cos(0 N2 vt)cos()NN 2 vySin(0 N2 vy )dy2)dy)dt2sin(-NN 2 vy N 2 vy sin( )dy sin( )dy0 N y0 N yN 2 vy N 2 vt sin

15、()dy sin()dt0 N 0 Nvy 2 vty )sin()dydtN0 cos(u(t y)dydt si n( 2-v(NNN0 2 vN2(2 v)2y)dy:Fsi n(晋 y)dy0 2 v Ncos(2 v) 1N22 cos(2 v) 1 (2 v)2由(4)式可知：当v 0时，F(u,v) 0当v 0时，由(1)式可知：F(u,v) N与实验相符4.2单条斜线段垂直现象的数学证明推广到x-y平面上幅度值相同的任意斜线斜线 啲方程a,b,c为常量l不为水平线f(x,y) y ax可定义为：f(x,y)cy ax b00假设I与正方形区域X N 1、.亠y N 父与(为，y

16、1),(x2,y2)其中0则可知在区间X1 xF(u,v)N 1N 12X1 X2N 1X2内f (x, y)才有值，这里令X1 X2是为了后面讨论方便2X2ccos(2x X1f (X, y)cos(2 怜 3)2,ux(Nv(ax b)NN 1 N 12f (X, y)sin(2 (UX )X 0 y 0NN2v(ax b)X2x.csin(2 (UXX2 ccos(2x Xix(uNva)2 vb) N )n(十22 vb)N )x22 x(u va)X1 ccos-2-vb)dxNx22 x(u va)X1 csin(对傅立叶变换后图像空间域与频率域中垂直现象的研究15X22 x(u

17、va) 2 vb、，2 t(u va)2 vb、丄ccos( )dx ccos( )dtNX12 vbx2)dx csin(NX1N2 vb ,2 t(u va)cos(N2 vb)sin(Nvb (2)cos(x1c2Nx22 x(u va)csi n(XlX2Xlc2X2NX2, 2 x(u va)cos(X1X2x1Nx22 x(u va)x sin( 丿N2 x(u va)X2cos(X1N.2 x(u va) Sin(X2 X2cosX1 X1(5)式 c2 X2x12 cX2xiX2xicossinN2 (U va)(x(X(Xt) dxdtx2)sin2与 2 2 cos(u v

18、a)(X2 X1)贝y.式 c2(x22X1)2 芝(X2xi)NN2 t(U va) 2 vb)dtN2 t(u va)Nt(u va)2N ， Nvb2 t(u va)si n(t) dxdt(X x1) dx易证得：sin22 1.2 Sin lim 20 2(u va)(X2Xi)N2 vb )dxdtN2 vb )dxdtNvb)2Nvb )dxdt N当趋近于零的时候出现最大值也就是:0时出现最大值，这里有X1X2(x2 X,)0，因此:1v u时出现最大值a对比直线I方程y ax b 可知：Fourier变换后在u-v平面上 幅度最大值将出现在与直线I垂直的方向上当u va 0时

19、(6)试出现最大值(5)5.垂直现象的物理意义5.1图像傅立叶变换的物理意义由于空间是三维的，图像是二维的，因此空间中物体在另一个维度上的关系就由梯度 来表示，这样我们就可以通过观察图像得知物体在三维空间中的对应关系。而像素之间的 梯度关系在傅立叶变换后表现为频率，也就是说图像的频率是表征图像中灰度变换剧烈成 都的指标。假设f （x, y）为一个N N的图像则它的离散二维傅立叶变换可定义为：1 N 1 N 1- f （x,y）eN x 0 y 0F(u,v)ux vyj2 (卞-N )其傅立叶反变换为:f (x, y)1 N 1N 1j21 F(u,v)eN u 0 v 0对上述的傅立叶正变换

20、做一些变形可得到下面的式子:1 N 1 1 N 1F(u，v) 7Nx0 7Ny0f(x，y)ej2() e j2 即可以看到经过变形后，二维的傅立叶变换可以通过两次一维傅立叶变换得到。因此对图像 进行傅立叶变换可以先对行（水平方向）进行一维的傅立叶变换，得出的结果再对列（垂 直方向）进行一维的傅立叶变换。这样就是使得图像的傅立叶变换实际上把频率分成水平 分量和垂直分量，即u分量和v分量。用这个思想再去分析图 3.3b可知：图3.3b是一幅 只包括单条水平线的简单图像，显然它只在垂直方向上有灰度跳变（垂直方向是梯度的方 向），所以从它的频谱图中只能看到垂直分量。这样就能很好的解释变换前后所出现

21、的垂 直现象。5.2从公式角度解释垂直现象的物理意义连续二维傅立叶变换可以表示为：F(u,v)f(x,y)e j2 (xu yv)dxdy（,）表示沿着（u,v）方向的单位向量；r （x,y）表示在xy平面上沿着（x,y）方向的向量xu yvi2 (xu yv)j2 ( )j 2(r?n)则有:ej( y) eej代回原式可得：F(u,v)f(x,y)ej2 (r?n)dxdy(u,v)方向上的投影。这样就把图像向相应的(u,v)所指向 fringe P atterns经过傅立叶变换后，在频谱图中最亮的内积(rr ? n)表示空域内点在向量的方向投影。这就能很好的解释点是沿着fringe pa

22、tterns 变化的方向。如图3.3a的fringe patterns : sin(x y)，当把 图像向u=v (或者说是x y)方向投影并以f(x,y)灰度值为纵坐标时，可以近似得到图 5.2.1中的右图，可以看出f(x,y)在u=v这个方向上有固定频率0 Ju； v2物理意义就是:输入函数sin(x y)的所有点在u=v (或x y)方向的投影所得到的函数只有一个频率 为0的正弦分量，也就是说所有的点在u=V方向上只能在(u0,V0)上有值。这样就会产生“叠加”的效果，我们就可以在u-v平面上看见两个亮点。现在就能解释为什么对图像进行傅立叶变换后常常会出现水平和垂直的“亮线”，如图5.2