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

[练习]对傅立叶变换后图像空间域与频率域中垂直现象的研究 对傅立叶变换后图像空间域与频率域中垂直现象的研究 【摘要】本文就图像二维傅立叶变换常令人感到困惑的问题进行深入的讨论，并对傅立叶变换后原图和频谱图出现的垂直现象进行分析，同时给出数学证明。 【关键词】傅立叶变换;频谱图;垂直现象;图像 1 引言 傅立叶变换是线性系统分析的一个有力工具，是信号处理中最重要、应用最广泛的变换，但遗憾的是很多人可能还是不习惯在频域中思考问题，尤其是面对图像时，空间域、频率域都是二维的，更是对空域、频域的对应关系不甚了解。如何理解傅立叶变换后的频谱图，为何对不同图像进行傅立叶变换后产生的频谱图往往会出现水平和垂直的“亮线”，而对一些规则图像进行傅立叶变换后得到的频谱图往往会在与原图垂直的方向上出现较亮的点或线(如垂直条纹图像，傅立叶变换后会在水平方向上出现一些较亮的点)。这些问题困扰着每个初学者，更让初学者困惑的是，几乎市面上所有有关数字图像处理的书籍都没有给出详细的解释。下面就围绕傅立叶变换以及图像处理中一些经常让人困惑的问题进行深入的讨论。 [1]2一维傅立叶变换中的问题讨论 一维傅立叶变换是二维变换的基础，下面就先对一些相关的、基本的但又很重要的概念进行讨论。 (1)实信号的傅立叶变换 就信号处理来说，大家所关心的都是实信号，所以单独对它进行讨论，可以简化工作。 下面是大家所熟知的傅立叶变换公式 , Fftitdt()()exp(2),,,,, ,,, 从定义式不难推出 ,*Fftitdt()()exp(2),,,, ,,, *FF()(),,,, 所以可得出结论:。这说明实函数的傅立叶变换是实部为偶函数，虚部为奇函数，也就是说:求某一实函数的傅立叶变换时，它的幅度谱总是关于原点对称的，而相位谱左右两边只是差一个负号，即左右互为复共轭。由此可见，就实际应用来讲，无 论哪边的频谱都是完备的，并且负频率本身也不具有什么意义，但是当用更为通用的数学方法去对物理过程建模时，保留负频率部分会使分析更加容易。 (2)正弦分量的分解 大家知道，任何满足狄里赫利条件的信号都可以通过傅立叶变换 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 征为一组正(余)弦信号的和或积分，而由图1可以从频域得出这一结论: 由于任何实函数的傅立叶变换都是偶函数，所以对于任意一个实信号，将它的傅ft()立叶变换在频率域进行抽样时，都能得到无数的抽样脉冲对(见图2)，也即将这些脉F(), 冲对累加起来就得到了，而当脉冲宽度趋于零时，每个脉冲对正好是某个频率的余弦F(), 信号的傅立叶变换，这从频率域的角度验证了任何一个实信号都可以看作是由若干个正(余)弦信号以及相应的幅度所组成。 图2 实信号及傅立叶变换 ft()F(), 通过以上的讨论可得出两个简单的结论:(1)实信号的频谱是对称的;(2)信号在时域和频域中是相互对应的，总是能把实信号看作是由若干不同频率、振幅的正弦波组成。 3.图像的二维傅立叶变换 3.1二维离散傅立叶变换的定义 图像经数字化处理后，可以用二维离散信号fmn(,)表示。对于二维离散信号{(,)|0,1,,1;0,1,,1}fmnmMnN,,,,？？，其离散傅立叶变换定义为: MN,,111munv (1),,，,Fuvfmnj(,)(,)exp(2()),,MNMNmn,,00 式中，称为空间频率。反变换定义为 uMvN,,,,0,1,,1;0,1,,1？？ MN,,111munv (2),，,fmnFuvj(,)(,)exp(2()),,MNMNmn,,00 式中。 mMnN,,,,0,1,,1;0,1,,1？？ NN，MN, 在图像处理时，一般选取图像块为的方阵，即取，这时二维离散傅立叶变换和反变换式: NN,,111munv为 (3),,，,Fuvfmnj(,)(,)exp(2()),,NNNmn,,00 NN,,111munv及 (4),，,fmnFuvj(,)(,)exp(2()),,NNNmn,,00 NN，在(3)(4)两式中，。本文都是选取的图像进行讨论的。uvmnN,,,0,1,,1,,？ [2]3.2频谱图的理解 由(3)式可知图像经傅立叶变换后，往往得到的是复数形式。要直接表示结果就必须用到两幅图像:一幅表示实部，一幅表示虚部。这样表示十分不方便，同时也没有得到有用的信息，因此引入变换结果的模作为值在频谱图中表示出来，以灰度的明暗代表模的大小。 作为典型的二维信号，图像的频率相应地也是二维的。其分别对应着图像的像素值在两个相互垂直的方向上变化的情况(如图3.2.1)。根据对一维离散信号频谱的分析结果可 fffff,,,0,,0 知，频谱在三点处的频率分别为为图像信号的最uNN,,0,2,1021NcNc, 高截止频率。将其推广到二维，则在图像频谱图中所有沿方向的频率值变化情况与一维u ffff,,,0,,0相同，也有。同理，在方向上也有相同的结果。因此，在频谱图v021NcN, 四角处沿u和v方向的频率分量均为0，在频谱图中(0,0), (1,0), (0,1), (1,1)NNNN,,,, f心点(2,2)NN处沿u和v方向的频率分量均为最大值。 c N/2N-10 v AB N/2(N/2,N/2)(N/2,N-1) CD N-1(N-1,N/2)(N-1,N-1) u 图3.2.1 图像频谱示意图 由于图像中的大部分能量集中在低频分量上，因此频谱图中四角部分的幅度值较大。然而，在实际的图像频谱分析过程中，由于低频分两区域较小，并且分散在四角，因此不利于对其进行分析。此时，可以根据图像频谱的周期性和共轭对称性对频谱图坐标进行移位，将所有低频分量集中在频谱图中心，同时高频分量分散在四周。在具体实现频谱的移位过程中，将图3.2.1中的区域A和区域D对换位置，将区域B和区域C对换位置即可。 .2.2所以。 进行移位之后大大增加了图像频谱的可读性，如图3 (a)原图像 (b)移位前的频谱 (c)移位后的频谱 图3.2.2 简单矩形图像坐标移动前后的频谱图 .3频谱图的垂直现象及单条直线段垂直现象3 对图像进行傅立叶变换后所得到的移位后的频谱图往往会在垂直方向和水平方向上出现两条“亮线”(如图3.2.2)，这个现象常常困扰初学者，为什么会出现这两条亮线，什么时候会出现这两条亮线,要很完整的回答这些问题很不容易，因为输入图像可以有很多种，但产生的频谱图大多都有这两条亮线。因此应该选取一些特殊的图像进行研究，接下来就选取:单条直线段，fringe patterns图像进行讨论。如图上方的为原图，下方为傅立叶变换并后未移位的频谱图: (a) sin(x+y) (b) x = a (c) x + y = N-1 图3.3 一些特殊图像及其傅立叶变换后的频谱图 [3] 图3.3(a)是fringe patterns sin(x+y)的图像，可以看到图像变换后的频谱图中出 [4]现两个点，这两点分别与原点及(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，同样出现了垂直现象。 上述对一些比较特殊的图像进行了傅立叶变换后中都出现了原图像与频谱图存在一定垂直关系的现象。对这样的垂直现象的研究能较好的理解频谱图和二维傅立叶变换，接下来就给出垂直现象的数学证明和物理意义的解释。 4.垂直现象的数学证明 本节给出只由一条直线段组成的图像的数学证明，在4.1中给出对图像内任意一条水平的直线段的证明，在4.2节中给出任意一条斜线段的证明。 因为频谱图中的灰度代表傅立叶变换的模值，因此证明变换前后出现垂直现象，实际上就是求傅立叶变换后模将会在什么位置或者说什么方向上出现最大值。即讨论求模公式: 1222NNNN,,,,,,1111，，，，uxvyuxvy,,(,)(,)cos(2())(,)sin(2())(1),,？？？？？Fuvfxyfxy,，，，,,,,,,,,,,NNNNxyxy,,,,0000，，，，,,,, 的最大值将出现在的什么方向上。 uv, 4.1水平直线段垂直现象的证明 设输入图像为一个的二维矩阵，图像内任意一条水平直线段定义如下:(1,1)NN,, 设图像内任意一条强度为宽度为一个像素的水平线为:1 fxyfxy(,)1(,)0,,,, 且有 ,,xaaNyNxa,,,,,, 0,1,21 ; 0,1,21？？,, 带入(1)式可得: 22NN11,,，，，，uavyuavy2Fuv(,)cos(2())sin(2()),，，，,,,,,,,,NNNNyy00,,，，，， 22NN11,,，，，，22222222uavyuavyuavyuavy,,,,,,,, (coscos sinsin)(sincos cossin),,，，,,,,,,NNNNNNNNyy00,,，，，， NN,,112222,,,,uavyuavy2222 (c,os)(cos)(sin)(sin)，,,NNNNyy,,00 NN,,112222uavyuavy,,,,2222 (cos)(sin)(sin)(cos)，，,,NNNNyy,,00 NN,,1122vyvy22,, (cos)(sin)..................,，.........................................................(2) ,,NNyy,,00 由(2)式可以看出，模的值只可能出现在的一列或者多列，也就是说在变换后的频v 谱图里亮线可能出现在垂直的方向上，这就很好的证明了图3.3(b)中出现的垂直现象。但 v,0是(2)式并不能说明模值只出现在这一列，因此还要对(2)式进行进一步的讨论。 N,1N 如果把(2)式中的连加看成是对积分的取样，并把扩充到，则可以得到如下的(3)式，并对其计算可得以下的等式: 22NN22,,vyvy，，，， cos()sin()..................................................(3),，dydy,,,,,,00NN，，，， NNNN2222vyvyvyvy,,,, cos()cos() sin()sin(),，dydydydy,,,,0000NNNN NNNN2222,,,,vyvtvyvt,， cos()cos() sin()sin()dydtdydt ,,,,0000NNNN NN2222vyvtvyvt,,,, cos()cos()sin()sin(),，dydt,,00NNNN NN2v, cos(()),,tydydt,,00N NNNv2,Nv2,)) sin() dyydy， sin((,,Ny,,002vN2vN,,22NN cos(2)1 cos(2)1........................................(4),,,,vv,,,,,,22(2)(2)vv,, 由(4)式可知:当时，当时vFuvvFuvN,,,,0 (,)0 0 (,) ,由(1)式可知:与实验相符 4.2单条斜线段垂直现象的数学证明 x-yl推广到平面上幅度值相同的任意斜线斜线的方程 fxycfxy(,) (,)0,,,,可定义为:为常量不为水平线abcl,, ,,yaxbyaxb,，,，,, 01,,,xN,假设与正方形区域交与其中lxyxyxxN(,),(,)01,,,,,11221201,,,yN, 则可知在区间内才有值，这里令是为了后面讨论方便xxxfxyxx,,,(,)1212 22NNNN,,,,1111，，，，uxvyuxvy2Fuvfxy(,)(,)cos(2(,,，，，))(,)sin(2())fxy,,,,,,,,,NNNNxyxy,,,,0000，，，， 22xx22，，，，uxvaxbuxvaxb()()，， cos(2())sin(2()),，，，cc,,,,,,,,NNNNxxxx,,，，，，11 22xx22，，，，2()2xuvavb，2()2,,xuvavb，,, cos(),，c，，csin(),,,,,,NNNNx,,xxx，，，，11 22xx2()22()2xuvavbxuvavb，，22，，，，,,,, cos()sin(),，，，cdxcdx,,,,,,xx 11NNNN，，，， xx2()22()2,,,,xuvavbtuvavb，，22 cos()cos(),，，cdxcdt,,xx11NNNN xx2()22()2xuvavbtuvavb,,,,，，22 sin()sin()，，，cdxcdt,,xx11NNNN 2()2,,tuvavb，22xx2()2xuvavb,,，2cos()，dxdt, cos(),，c,NNxx11NNxx22()22()2xuvavbtuvavb,,,,，，22 sin()sin()，，，cdxdt,,xx11NNNN xx22()22()2xuvavbtuvavb,,,,，，22 cos()cos() ,，，c,,xx11NNNN2()22()2,,,,xuvavbtuvavb，， sin()sin()，，dxdt ，NNNN xx2222(),uva， cos[()] ................................................ (5),,cxtdxdt,,xx11N 2()uva，,设则:,, Nx2x22 式,x(5)cos[()],,cxtdxdt,1,x12cx2 sin()sin(),,,,,,,xxxxdx,,21, ,x1 2c,, 22cos(),,,,xx212, ,()()uvaxx，,21设则:,,N 2sin,22 (6)()式,,cxx212, 易证得: 22sinsin,, 1 lim1 ,,,,当趋近于零的时候出现最大值也就是:,22,, 0,, ()()uvaxx，,21, 0()0,,,,,时出现最大值，这里有，因此:xxxx1221N 1 当时试出现最大值uvavu，,,,,0(6) 时出现最大值a对比直线方程可知:变换后在平面上lyaxbuv,， Fourier-幅度最大值将出现在与直线垂直的方向上l 5.垂直现象的物理意义 5.1图像傅立叶变换的物理意义 由于空间是三维的，图像是二维的，因此空间中物体在另一个维度上的关系就由梯度来表示，这样我们就可以通过观察图像得知物体在三维空间中的对应关系。而像素之间的梯度关系在傅立叶变换后表现为频率，也就是说图像的频率是表征图像中灰度变换剧烈成都的指标。 NN， 假设为一个的图像则它的离散二维傅立叶变换可定义为: fxy(,) uxvy,,11NN,,，2()j1NNFuvfxye,(,)(,) ,,N,,00xy 其傅立叶反变换为: uxvyNN,,11,j2()，1NNfxyFuve,(,)(,) ,,Nuv,,00 对上述的傅立叶正变换做一些变形可得到下面的式子: uxvyNN,,11，，,,,,jj2()2()11NNFuvfxyee(,)(,), ,,,,NN,,00xy，， 可以看到经过变形后，二维的傅立叶变换可以通过两次一维傅立叶变换得到。因此对图像进行傅立叶变换可以先对行(水平方向)进行一维的傅立叶变换，得出的结果再对列(垂直方向)进行一维的傅立叶变换。这样就是使得图像的傅立叶变换实际上把频率分成水平分量和垂直分量，即分量和分量。用这个思想再去分析图3.3b可知:图3.3b是一幅uv 只包括单条水平线的简单图像，显然它只在垂直方向上有灰度跳变(垂直方向是梯度的方向)，所以从它的频谱图中只能看到垂直分量。这样就能很好的解释变换前后所出现的垂直现象。 5.2从公式角度解释垂直现象的物理意义 连续二维傅立叶变换可以表示为: ,jxuyv2(),,，Fuvfxyedxdy(,)(,), ,,,, 22uv令,，, uv,, (,) (,) (,) (,)nuvrxyxyxy,,表示沿着方向的单位向量;表示在平面上沿着方向的向量 ,, xuyv,,j2(),,，jxuyvjrn2()2()，,,,,,,eee,,则有: 代回原式可得: ,,,jrn2( ),,,,Fuvfxyedxdy(,)(,), ,,,, ,,( )rn,内积表示空域内点在向量方向上的投影。这样就把图像向相应的所指向(,)uv(,)uv的方向投影。这就能很好的解释fringe patterns经过傅立叶变换后，在频谱图中最亮的点是沿着fringe patterns变化的方向。如图3.3a的fringe patterns:，当把sin()xy，图像向= (或者说是)方向投影并以灰度值为纵坐标时，可以近似得到图xy,uvfxy(,) 22,,，uv5.2.1中的右图，可以看出在=这个方向上有固定频率物理意义就uvfxy(,)000 是:输入函数的所有点在= (或)方向的投影所得到的函数只有一个频率xy,uvsin()xy， 为的正弦分量，也就是说所有的点在=方向上只能在上有值。这样就会产生,(,)uvuv000 “叠加”的效果，我们就可以在-平面上看见两个亮点。 uv 图5.2.1将图像投影到的方向上 (,)uv 现在就能解释为什么对图像进行傅立叶变换后常常会出现水平和垂直的“亮线”，如图5.2.2的右图是左图的频谱图，由前面的结论可知，频谱图中的水平亮线是由原图中背景的垂直条纹产生，而垂直亮线则主要是由原图中的镜框及一些灰度规律变化的水平线产生的。 图5.2.2原图像与傅立叶变换后的频谱图 6. 总结 初级经济法重点总结下载党员个人总结TXt高中句型全总结.doc高中句型全总结.doc理论力学知识点总结pdf 通过对单条直线段的简单图像的分析和证明，可以从数学角度和物理意义两个方面解释傅立叶变换后原图像与生成频谱图之间出现的垂直现象，从而能更好的理解傅立叶变换在图像处理中应用的意义，同时解决了在利用傅立叶变换对图像进行处理过程中的一些疑惑。 参考文献 [1] 张若岚，刘劲松.图像信号的频域理解[J].数字电视与数字视频.2002, 1002-8692 [2] 马平编著.数字图像处理和压缩[M].2007 [3] S. De Nicola a, *, P. Ferraro . A two-dimensional fast Fourier transform method for measuring the inclination angle of parallel fringe patterns[J]. Optics & Laser Technology 30 (1998) 167?173 [4] D.G. Bailey, Dept. of Phys., Massey Univ., Palmerston North, New Zealand . Detecting regular patterns using frequency domain self-filtering[J]. Image Processing, 1997. Proceedings., International Conference on Download from IEEE