《成像原理论》word版.docx

成像原理论文三维形貌成像摘 要本文研究了一种还原物体表面三维形貌的技术-傅里叶变换轮廓术。傅里叶变换轮廓术只需要一幅干涉图即可对整个范围进行分析，其测量系统简约，分辨率高，计算机处理数据快，计算精度高，且具有很强的非接触测量优势，使其得到广泛的应用。在实际傅里叶变换轮廓术测量中，获取的条纹图扩展的零频分量对傅里叶变换轮廓术的测量精度和测量范围有很大影响，甚至妨碍正确三维面形的恢复。相移技术常被用来消除零频分量对测量的影响，但需要在测量系统中安装精密相移装置，并需要采集两帧具有相位差的条纹图。传统傅里叶变换轮廓术中，完成精密相移需要较长的时间，影响了傅里叶变换轮廓术测量方法的实时性。本文采用双色正弦条纹投影来实现从一帧条纹图中消除零频对傅里叶变换轮廓术测量的影响。该方法同传统的相移方法相比， 不需要相移装置，测量系统简单，并且能真正实现高速测量。最后介绍了基于线阵CCD的FTP测量技术的应用，分别设计了静态物体测量系统和振动过程测量系统，并且分析了系统中潜在误差来源。关键词：傅里叶变换轮廓术；三维形貌测量；双色正弦条纹投影1 三维形貌成像1.1 光学三维成像简介光学三维成像就是指用光学的手段获得物体三维空间信息的方法和技术，目前主要指获得物体表面三维空间形状信息的方法和技术。随着计算机技术、信息技术的迅速发展，极大地改变了传统的光学计量技术。光学计量初期所采用的感光胶片记录方式已为固态摄像机技术所取代，高性能的微型计算机和图像处理系统使光学图像的计算机辅助分析技术迅速发展，这些信息的获取和处理技术上的进步又给光学传感和计量方法上的革新和发展以新的活力，使新的三维传感和计量方法不断涌现。为此，国际光学学会1994年以信息光学的前沿为主题的年会上，首次将光学三维传感列为信息光学前沿七个主要领域和方向之一。物体三维信息的数字化是目前重要的研究课题之一，在机器视觉、实物仿形、工业自动检测、产品外观质量检查、生物及医学检测等领域中具有重要意义，日益受到人们的重视。相信看过电影侏罗纪公园的人一定会对影片中那些活灵活现的大恐龙感到疑惑，对片中大恐龙穿墙过壁、以假乱真的场面记忆犹新，这种影视特技到底是怎么做出来的呢？这就是三维传感技术的魅力。在特技制作中，技术人员先雕刻好一个恐龙的模型，然后用光学三维传感方法得到恐龙的三维彩色数字模型，再用三维动画软件使其做出各种动作，并完成与背景、人物的合成，最后才形成我们看到的惊心动魄的场面。另外，尤其在某些应用场合常常要求非损伤式、高精度、高速度的获取物体的表面形状信息，例如沙体形态三维测量、井下模式识别、产品的自由曲面检测、医学诊断中对人脸和牙齿等生物器官的三维造型等。因此，对非接触、快速、高精度的三维面形测量技术的研究更为重要。近些年来，随着半导体激光器、数字投影仪、电荷耦合器件、高分辨率数码相机、高速摄像机等高性能光电设备的出现，以及基于计算机的数字图像处理技术的广泛应用，集光、机、电、计算机技术于一体的光学三维传感技术得到了很大发展。其中，基于结构光照明的主动三维传感方法由于具有快速、非接触、高精度、易于实现自动化等测量特点，日益受到人们的重视。在各种主动三维面形测量方法中，傅里叶变换轮廓术由于只需获取一幅条纹图就可计算出物体的三维面形信息，特别适合实时动态测量，因此受到众多研究者的关注。1.2 光学三维面形测量方法概述获取三维面形信息的基本方法可以分为两大类：被动三维传感和主动三维传感。被动三维传感采用非结构照明方式，从一个或多个摄像系统获取的二维图像中确定距离信息，形成三维面形数据。从一个摄像系统获取的二维图像中确定距离信息时，人们必须依赖对于物体形态、光照条件等的先验知识。如果这些知识不完整，对距离的计算可能产生错误。从两个或多个摄像系统获取的不同视觉方向的二维图像中，通过相关或匹配等运算可以重建物体的三维面形。双摄像机的传感系统如图1.1所示，它与人眼双目立体视觉的原理相似。从两个或多个摄像系统获取的不同视觉方向的二维图像中确定距离信息，常常要求大量的数据运算。当被测目标的结构信息过分简单或过分复杂，以及被测目标上各点反射率没有明显差异时，这种计算变得更加困难。因此，被动三维传感的方法常常用于对三维目标的识别、理解，以及用于位置、形态分析。这种方法的系统构成比较简单，在无法采用结构照明的时侯更具有独特的优点。随着计算技术的发展，运算速度已不再是一个主要的限制因素。在机器视觉领域已经广泛地应用被动三维传感技术。图1.1 被动三维传感：双摄象机系统主动三维传感采用结构照明方式。由于三维面形对结构光场的空间或时间调制，可以从携带有三维面形信息的观察光场中解调得到三维面形数据。由于这种方法具有较高的测量精度，因此大多数以三维面形测量为目的的三维传感系统都采用主动三维传感方式。作为结构照明所采用的光源，原则上可以采用激光光源和普通光源。激光具有亮度高、方向性和单色性好、易于实现强度调制等优点，所以在很多应用领域常常采用以激光为光源的三维传感系统。采用白光光源的结构照明方式具有噪声低、结构简单的优点，特别是在面结构照明的三维传感系统中受到愈来愈多的重视。1.2.1 激光扫描法激光扫描三维轮廓术是目前应用最广，技术最成熟的三维测量技术。根据测量物体的要求，激光光束可以是点、线，还可以是环。其中以线结构光使用的多，称之为光刀。图1.2是线结构光测距原理示意图。在待测定距离的场景前放置光源和相机。光源可绕一垂直方向的轴旋转，对景物进行扫描以覆盖整个待测定的场景。在扫描过程中，相机不断从场景获取图像。光源从场景的一端扫描到另一端，摄像机采集到的是图像序列。图1.2 线结构光测量示意图结构光扫描法几乎可以适应任何物体的外形轮廓测量。由于扫描激光的强度是二值化的，并且光强较高，所以不会受到物体表面纹理的干扰。它的缺点是需要一个扫描装置，而且测量速度慢，不能对动态物体进行实时测量。每当激光移动一个角度就要获取一幅图象。在整个测量过程中要处理大量的图象。而每幅图象中有用的信息只占整个图象的一小部分，造成大量的信息冗余。 1.2.2 相位测量轮廓术相位测量轮廓术（Phase Measuring Profilometry，简称PMP）是借鉴相移干涉术的相位计算方法发展而来的。它的基本思想就是采用正弦条纹投影和数字相移技术，通过获得全场条纹的空间信息和一个条纹周期内相移条纹的时序信息，重建物体的三维面形。PMP测量系统光路图如图1.3所示，图中光栅为由相移装置移动的投影光栅，成像设备采用CCD，投影系统的光轴P1P2与成像系统的光轴I1I2交于参考平面上的点O，它们的夹角为，当一个正弦光栅（或罗奇光栅）被投影到三维漫反射物体表面上时，成像系统获取的由物体高度所调制的变形光栅像为： I(x,y)=R(x,y)A(x,y)+B(x,y)cos(x,y) (l.1)其中R(x,y)为物体表面不均匀反射率，A(x,y)表示背景光强，B(x,y)/ A(x,y)是条纹的对比度（条纹的对比度是一个等于或小于1的数）。相位函数(x,y)包含受到被测物体高度调制而引起的变形信息。 图1.3 PMP测量系统光路图直接从强度分布确定相位分布(x,y)是比较困难的，相移算法却提供了一种精确测定相位分布的手段。以四步相移算法为例，即相位移动的增量为/2，所产生的四个变形光栅图像为： I1x,y=Ax,y+Bx,ycosx,y I2(x,y)=A(x,y)+B(x,y)cosx,y+/2 I3(x,y)=A(x,y)+B(x,y)cos(x,y)+ I4(x,y)=A(x,y)+B(x,y)cos(x,y)+3/2 (1.2)由以上方程组可得： (x,y)=I4(x,y)-I2(x,y)/ I1x,y-I3(x,y) (1.3)对于更为普遍的N步相移算法，当投影的正弦光栅被移动其周期的l/N时，条纹图的相位被移动2/N产生新的相位强度，依次相移N-1次，可以从N帧相移条纹图中计算出相位函数： (x,y)=tan-1n=1NIn(x,y)sin(2n/N)/n=1NIn(x,y)cos(2n/N) (1.4)式中In表示由CCD获取的第n次变形光栅像强度分布。计算所得到的相位(x,y)被截断在反三角函数主值区(-, )内，要对其进行相位展开，以得到连续的相位分布。再通过相位和高度的对应关系计算出物体表面的高度分布。相移技术的优点是对相位测量的精度可以达到几十分之一到几百分之一等效波长（条纹周期）；同时对条纹背景、对比度和噪声的变化不敏感；计算得到相位值是一个均匀分布的正交网格上的点的测值，测点与探测器阵列或图像处理板上的阵列一一对应，有利于进一步的信号处理，实现自动的三维面形测量。1.2.3 傅里叶变换轮廓术1983年，M.Takeda和K.Muloh将傅里叶变换用于三维物体面形测量，提出了从频域来获取物体三维面形的方法，即傅里叶变换轮廓术（Fourier Transform Profilometry，简称FTP），由于其只用一帧或两帧条纹图就可以恢复物体的三维面形，具有很强的实时性，因此在物体动态三维面形测量中得到了广泛的应用。这种方法以正弦光栅（或罗奇光栅）产生的结构光场投影到待测物体表面，得到被三维物体面形分布调制的变形条纹光场，成像系统将此变形条纹光场成像于面阵探测器上，然后利用计算机对图像的强度进行傅里叶变换、频谱滤波和逆傅里叶变换等处理，求解出物体的高度分布信息，得到物体的三维面型分布。该方法具有单帧获取、全场分析、高分辨率和易于处理动态过程等特点而受到人们的广泛关注。三维传感的方法很多，以上介绍的是最基本几种。随着计算机技术、光学和光电子技术的飞速发展，新的光学三维传感和计量方法也不断进步。1.3 傅里叶变换轮廓术研究现状和发展趋势FTP与时间飞行法(TOF)、激光逐点扫描法、光切法等光学三维轮廓术相比，具有全场、快速的特点，不需要扫描装置；同PMP比较，FTP只需要一幅或两幅条纹图就可以求得物面的三维信息，省去了复杂的相移装置，设备简单，测量速度快，在实时和动态测量领域具有广泛的应用前景。因此，FTP成为当前光学三维面形测量技术研究的重点，众多研究者对FTP进行了深入的研究和进一步改进。在提高FTP测量精度方面，如：用准正弦投影和相移技术消除零频对携带有用信息的基频影响，减小甚至消除频谱混叠；针对计算机处理的是离散信息这一前提，为了完善FTP的理论分析，研究了抽样对FTP测量精度的影响，并给出了抽样频率的选择依据；用二维汉宁滤波窗口进行加权滤波来减小噪声的影响；根据条纹频率与相位可靠度之间的关系，采用数字加权滤波建立调制度函数作为衡量相位可靠性的模板，指导截断相位沿一条最佳路径展开消除条纹不连续对测量的影响；采用多角度彩色条纹投影解决阴影问题；采用双频光栅和双色光栅的快速FTP等；研究了频谱泄漏和CCD量化误差对基于条纹投影的FTP测量精度的影响。在FTP的应用方面：将FTP轮廓术用于高度不连续性或表面孤立复杂物体的三维面形测量；将FTP用于实时动态或360度三维面形测量。 目前FTP的测量精度还有待于进一步提高，在实际应用中还存在若干困难，比如不完善相位图无法展开的问题，频谱混叠的影响等。未来的研究重点是进一步提高FTP的测量精度，克服频谱混叠，不断优化频域滤波和相位展开算法，尽量减少测量过程中各种干扰因素，以满足相关领域对三维面形数据快速、自动、实时获取的要求。本文中主要利用相移技术消除零频对携带有用信息的基频影响，由于传统的相移技术需要两帧条纹图及复杂的实验装置，我们考虑用双色正弦条纹投影光栅来解决以上问题，可以简化实验装置并且只需一帧条纹图。2 傅里叶变换轮廓术2.1 傅里叶变换定义及条件利用公式： fx=-W(f)(x;f)df (2.1)对非周期函数g(x)作正交展开，即把函数看作复指数函数在整个连续的频率区间上的积分和： g(x)=-G(f)ej2fxdf (2.2)式中： G(f)=-g(x)e-j2fxdx (2.3)式(2.2)和式(2.3)这两个积分即傅里叶积分，G(f)称为 g(x)的傅里叶变换或频谱。若g(x)表示某空间域的物理量，G(f)则是该物理量在频率域的表示形式。G(f)的作用类似于傅里叶系数Cn，即作为各种频率成分的权重因子，描述各复指数分量的相对幅值和相移。当G(f)是复函数时，可以表示为： G(f)=A(f)ej(f) (2.4)式中，A(f)=G(f)，是g(x)的振幅频谱；(f)是g(x)的位相频谱。非周期函数的频谱不是离散的，而是频率连续或分段连续的函数。傅里叶变换的意义：所有适当加权的各种频率的复指数分量叠加起来就得到原函数g(x)，称它为G(f)的傅里叶逆变换。g(x)和G(f)构成傅里叶变换对，常简记为： G(f)=Fg(x) , g(x)= F-1G(f) (2.5)二维傅里叶变换只是一维傅里叶变换的推广： gx,y=-G(fx,fy)expj2fxx+fyydfxdfy=F-1G(fx,fy) (2.6) Gfx,fy=-gx,yexp-j2fxx+fyydxdy= Fgx,y (2.7)傅里叶变换存在的条件：(1) gx,y在整个x y平面绝对可积，即：-|g(x,y)|dxdy0 (2.25)图2.4 正弦光栅投影时变形条纹的一维频谱图其中：(fQ(f+f0,y)min代表基频分量Q(f+f0,y)的最小值。所以基频最大可以扩展到频域轴上的零点，所以正弦投影光栅的FTP测量范围限制为： h(x,y)xmaxl/d (2.26)其中：l/d=1/tg （是投影光轴与成像光轴之间的夹角）。我们通常称式(2.26)为FTP测量系统的结构条件。另外，正弦光栅投影的傅里叶变换频谱，除了保证同一个周期内正负基频分量之间不混叠外，还要保证基频分量与相邻周期的频率分量不混叠。即：要防止抽样不当引起的频谱混叠，就必须保证Q(f-f0)的频谱同相邻周期Q*(f+f0-mf0)的频谱相分离。相邻周期间的频谱分布如图2.5所示： 图2.5 相邻两周期间的频谱分布示意图当m取得足够大时，相邻的“频谱岛”是分离的；当m取得太小，不满足抽样定理时，相邻的“频谱岛”就会发生混叠。设抽样频率为光栅频率f0的m倍，利用调制信号的瞬时频率概念可以得到： f0+12xmaxmf0-(f0+12xmax) (2.27) 结合式(2.17)式进一步化简得： hxmax(m-22)ld (2.28) 结合式(2.28)和同一个周期内的频谱分离条件式(2.26)可得： (m-22)ldld (2.29)故由式(2.28)可得：m4 (2.30)上式(2.30)表明，使用正弦光栅投影的FTP，一个条纹周期至少需要4个抽样值才能保证周期间频谱不发生混叠。所以说为了保证FTP测量精度，一方面必须满足抽样条件；另一方面被测物体的高度变化必须满足测量系统的结构条件。3 双色正弦条纹投影原理传统傅里叶变换轮廓术中，采用正弦光栅投影和相移技术获得2幅变形条纹图像，然后可以消除零频，留下有用基频分量。但是完成精密相移需要较长的时间，并且需要复杂的实验装置，既影响了傅里叶变换轮廓术测量方法的实时性，又增加了实验的难度。本文采用双色正弦条纹投影来实现从一帧条纹图中消除零频对傅里叶变换轮廓术测量的影响。该方法同传统的相移方法相比，不需要相移装置，测量系统简单，并且能真正实现高速测量。3.1 双色正弦条纹投影FTP测量流程传统FTP投影黑白光栅条纹在被测物体表面上，双色正弦条纹投影则采用红、绿、蓝三基色中的两种颜色相互间具有相位差的正弦条纹，通过RGB图像模式组合成双色条纹投影到被测物体表面。根据RGB图像模式，一幅彩色条纹强度分布可用如下方程表示： I(x,y,k)=kIk(x,y),(x=1n;y=1m;k=1,2,3) (3.1)其中Ik为一幅单色条纹的强度分布，k=1,2,3分别对应红、绿、蓝三种颜色。因此利用RGB色彩分离技术，可以从采集到的一幅双色条纹中分离出两个颜色分量的条纹，得到具有相位差的两幅正弦条纹图。通过对这两幅条纹图进行背景和对比度的矫正，使他们具有相同的背景和对比度，然后将矫正后的两幅条纹图相减可以得到消零频的变形条纹，相加可以得到条纹图中包含的背景分量。如果在实际测量中，将物体置于黑色背景中，那么通过对矫正后的两幅条纹图相加生成的背景分量进行二值化，便可以得到反映物到轮廓的二元模板。图3.1给出了红绿双色正弦条纹投影FTP的测量流程图。投影并获取基于RGB的双色条纹提取各个颜色分量的单色条纹红色条纹绿色条纹对条纹进行背景和对比度校正物体轮廓二元模板生成去零频的条纹二维FFT、频域滤波、二维IFFT调制度信息截断相位分布可靠性模板反映物体表面高度变化的连续分布相位三维物体的面形信息相加相减指导相位展开图3.1 双色正弦条纹投影FTP流程框图3.2 消零频原理若投影一幅红绿双色条纹，则红绿双色条纹的光场强度可表示为： I(x,y)=ag+bgcos(2f0x)+ ar+brcos(2f0x+) (3.2) ar和ag分别是红、绿条纹的背景分量，br和bg分别是红、绿条纹的对比度。为了简单，通常使ar=ag，br=bg。将此双色条纹投影到被测物体表面后，于投影光轴成一定夹角的方向上，获取的一幅双色变形条纹图像可表示为： Grg(x,y)=Rg(x,y) ag+bgcos(2f0x+(x,y) +Rr(x,y) ar+brcos(2f0x+(x,y)+) (3.3) Rr(x,y)、Rg(x,y) 分别是物体对不同颜色投影光的反射率分布。从式(3.3)所示的双色条纹中提取红、绿两个颜色分量，就可以得到具有相移的两幅条纹图： gg(x,y)=Rg(x,y) ag+ Rg(x,y)bgcos(2f0x+(x,y) (3.4) gr(x,y)=Rr(x,y)ar+Rr(x,y)brcos(2f0x+(x,y)+) (3.5)由于Rr(x,y)Rg(x,y)，当ar= ag，br=bg，红、绿条纹的背景和对比度仍不相等，需要进行矫正，使它们相等。如果把条纹图看作是条纹灰度在空间上的一种随机分布，当Rr(x,y)和Rr(x,y)变化比较缓慢的时候，可以认为条纹的背景分布趋近于条纹图的均值，条纹的对比度分布趋近于条纹图的均方差，因此可以通过调整红、绿条纹的均值和均方差相等，来使他们具有近似相等的背景和对比度，采用如下公式进行矫正(例如将红色条纹的背景和对比度矫正为绿色条纹的背景和对比度）： gr(x,y)=gr(x,y)-mrdg/dr+mg (3.6)其中mr、mg分别代表红、绿条纹灰度的均值，dr、dg分别代表红、绿条纹灰度的均方差。gr(x,y)为调整后的红色条纹灰度。理想情况下，通过以上矫正便得到了两组具有相同背景和对比度的，相互间具有相移的两幅条纹图： gr(x,y)=A(x,y)+B(x,y)cos(2f0x+(x,y)+) (3.7) gg(x,y)=A(x,y)+B(x,y)cos(2f0x+(x,y) (3.8)其中：A(x,y)= Rg(x,y) ag代表红、绿条纹具有相同的零频背景， B(x,y)= Rg(x,y) bg代表红、绿条纹具有相同的对比度。 于是利用式(3.7)、式(3.8)两组条纹相加可以生成条纹的背景分布图，利用式(3.7)、式(3.8)两组条纹相减可以消去零频。3.3 相位展开Mx,y=1 0 通常情况下，条纹的背景分布中并不能分辨出物体的轮廓，为了得到反映物体轮廓的二元模板，需要将被测物体置于黑色背景中。此时利用式(3.7)、式(3.8)两组条纹相加在物体轮廓区域内得到的是物体的较亮的背景，在物体轮廓之外得到的是暗背景。这样，条纹的背景图中物体轮廓内区域的灰度与周围黑色背景的灰度呈现出明显的差别，再利用最大方差准则确定灰度阈值进行图像分割（即通过选择阈值使得条纹图被此阈值分割的两部分区域的平均灰度与整幅条纹图的平均灰度的差别最大，区域间的方差作为描述这种差别的有效参数），将大于此阈值的点灰度设置为1，小于此阈值的点灰度设置为零，从而得到一个反映物体轮廓的二元模板，用一个二元矩阵表示即：M(x,y)=1 0 物体轮廓以内的点 物体轮廓以外的点 (3.9)这种条纹相加提取物体轮廓模板的方法，解决了物体轮廓内区域的暗条纹与周围暗背景无法区分的问题。利用式(3.7)、式(3.8)两组条纹相减，可以将零频分量减去。但是由于黑色背景的存在，将导致条纹的背景在物体轮廓处发生突变，如图3.2所示的一幅黑色背景存在下采集到的变形条纹的剖面图，其中虚线a是物体区域内的条纹灰度均值，实线b是整幅条纹的均值，两者有较大差别，因而物体区域内变形条纹的背景分布不再趋近于整幅条纹的均值，这样在应用式(3.6)矫正变形条纹的背景和对比度时，就会产生比较大的误差，影响减零频的效果。为了提高消零频的效果，可以利用式(3.9)得到的二元模板将物体轮廓区域内的红、绿条纹提取出来，然后利用式(3.6)对提取出来的条纹进行第二次校正后，相减消去零频，得到条纹g(x,y)。图3.2 一列真实条纹剖面图通过上述方法得到的条纹图g(x,y)，只在物体轮廓内区域有数据，因而条纹所包含的像素点数不够多，在对条纹数据进行离散傅里叶变换时，频域的分辨率不高。同时条纹在空域内的突然截止也会造成物体轮廓的边缘处产生误差。为了解决以上两个问题，需要在物体轮廓外扩充条纹数据。我们在保证其他条件不变的情况下，采集一幅无黑色背景的红绿参考条纹，根据式(3.6)对采集的红、绿参考条纹进行矫正，相减消去零频后，得到条纹g0(x,y)。利用g0(x,y)结合式(3.9)生成的物体轮廓二元模板便可以完成对条纹g(x,y)的扩充，其具体方法如式(3.10)所示，生成新的条纹g(x,y)如下： g(x,y)=g(x,y)+ g0(x,y)M(x,y) (3.10)M(x,y)=1 0 其中M(x,y)由Mx,y取反得到即： 物体轮廓以外的点 物体轮廓以内的点 (3.111)从式(3.10)可以看出，新生成的条纹g0(x,y)有两部分组成，一部分是原来的条纹数据 g(x,y)，一部分是由模板M(x,y)从条纹g0(x,y)中提取的。因此，使用这种方法扩充条纹，增加了参与离散傅里叶变换的总的采样点数，从而提高了频域分辨率，同时抑制了因条纹锐截止产生的边缘误差，当生成的条纹在连接处有断裂时，还可以进一步采用迭代算法弥补。对条纹g(x,y)做二维傅里叶变换，用二维汉宁窗滤出基频，做逆傅里叶变换，得到： g(x,y)= B(x,y)expj2f0x+(x,y) (3.12)对参考条纹g0(x,y)作相同运算得到： g0(x,y)=B0(x,y)expj2f0x+0(x,y) (3.13)(x,y)=(x,y)- 0(x,y)=arctan Img(x,y)g0*(x,y)Reg(x,y)g0*(x,y)由式(3.12)和式(3.13)得到由物体高度引起的相位差：(3.14) 调制度函数：m(x,y)=g(x,y) (3.15)Im表示取函数的虚部，Re表示取函数的实部，“*”表示复共轭。将调制度和物体轮廓二元模板结合起来形成可靠性模板指导截断在三角函数主值区域的(x,y)进行相位展开。利用物体轮廓的二元模板生成新的调制度函数如下： m(x,y)= m(x,y)M(x,y) (3.16)新的调制度函数将物体轮廓以外的点的调制度设为零，当此调制度函数作为模板指导相位展开时，一方面通过对调制度排序，优先展开调制度高的地方，可以避免条纹断裂、阴影、物体高度突变造成调制度过低的点对相位展开的影响；一方面在相位展开中可以忽略那些调制度为零的点，从而只对物体轮廓以内的有效点进行展开，节省了运算时间，避免了物体轮廓以外的相位不可靠点引起的相位展开错误。得到连续的自然位相(x,y)后，根据相位和高度的映射关系，得到物体的高度分布函数： h(x,y)=(x,y) K (3.17)其中K为系统参数，当在远心光路条件下，K=l/2f0d 。4基于线阵CCD的FTP测量技术的应用4.1 测量系统设计基于线阵CCD的FTP测量系统包括投影系统、成像系统和被测物三大部分。计算机生成双色正弦光栅通过投影仪成像，在被测物表面形成光栅的变形条纹。但由于线阵CCD只能采集一行光栅条纹图的特殊性，成像系统中线阵CCD必须垂直于被测物表面光栅分布方向，测量系统结构如图4.1所示。图4.1 测量系统结构示意图 lp为投影仪到被测物的水平距离，lc为线阵CCD相机到被测物的水平距离，d1为投影仪到线阵CCD相机之间的水平距离，d2为投影仪到线阵CCD相机之间的垂直距离。系统参数为lp=128mm，lc= 86mm，d1= 39mm，d2= 40mm。4.2 静态物体测量线阵CCD只能对单独的一行表面条纹图成像，得到完整的被测物表面条纹分布则需要通过扫描技术来实现。针对静止物体表面三维形貌的测量，设计基于线阵CCD的FTP移动扫描测量系统。根据前面所述的FTP测量原理，设计的实验系统如图4.2所示。计算机编写出数字光栅图像，由投影仪投射到物体表面上，形成分布在表面上的光栅条纹。线阵CCD采集一行条纹，将釆集到的数据传输给计算机，然后步进电机导轨沿X方向移动1mm，再进行下一行的条纹采集，如此重复直至完成整个面的扫描测量。通过移动扫描采集到被测物表面每一行的条纹图后，得到完整的被测物表面变形条纹图。然后，逐行进行一维傅里叶变换、提取基频、一维傅里叶逆变换，获得每一行的截断相位，用基于传统方法相位展开进行相位展开。最后，根据标定结果恢复每一行的高度分布，从而得到完整的被测物表面高度分布。图4.2 移动扫描测量装置结构示意图4.3 振动过程测量基于线阵CCD的FTP振动过程测量技术的测量对象是处于受迫振动过程中的物体。线阵CCD对被测物表面一行的测量，可视作对被测物截面进行的测量。测量过程中，连续采集处于振动过程中的被测物表面条纹图，并以时间为轴，显示每一时刻恢复的振动状态，从而恢复被测物在连续采集时间段内的振动过程。根据前文所述的FTP测量原理，设计的实验装置图如图4.3所示。计算机编写出数字光栅图像，由投影仪投射到物体表面上，形成分布在表面上的光栅条纹。将被测物固定振动源上，给定振动频率，待振动达到稳定状态后，由线阵CCD进行连续采集，完成条纹图采集后进行数据处理。连续采集到振动过程中被测物截面条纹图后，以时间为轴排列得到完整的被测物截面振动过程条纹图。然后，逐行进行一维傅里叶变换、提取基频、一维傅里叶逆变换，获得每一行的截断相位，用基于传统方法相位展开进行相位展幵。最后，根据标定结果恢复每一行的高度分布，从而得到连续采集时间内的被测物截面振动过程。实验系统结构设计图如图4.3所示：图4.3 振动过程测量装置结构示意图4.4 误差分析测量系统中产生误差的因素有多种，既有来自于系统的误差，也有基于原理本身的误差，同时还存在来自外部的误差。4.4.1光学系统引起的误差对于实际的光学系统来说，像差是不可避免的。光栅条纹图质量不够好的原因之一是系统中的存在像差。像差大小可以由一维傅里叶变换之后的频域的二次分量的强弱