一种光学干涉条纹快速处理方法
【专利摘要】本发明通过二维连续小波分析方法，针对数字干涉图，设计基于可编程逻辑器件的高速干涉条纹相位提取算法。本发明利用二维快速傅立叶变换实现二维连续小波变换（2D-CWT)，给出干涉条纹相位提取快速实现算法的FPGA实现，在实现中能对2D-CWT的缩放因子和旋转因子等参数的设置，以及对不同母小波（核函数）的选择。本发明对基于二维连续小波变换的光学干涉条纹相位提取技术具有重要意义，可为基于干涉图分析的光学测量技术在工业领域的应用提供有效可靠的分析方法。
【专利说明】一种光学干涉条纹快速处理方法
【技术领域】
[0001] 本发明属于光学测量和动态特性分析领域，涉及一种光学干涉条纹快速处理实现 方法。
【背景技术】
[0002] 在光学测量领域中经常需要对获取的光学千涉图（或称为干涉条纹）进行研究分 析。根据干涉方法的不同，常用的千涉图一般包括莫尔条纹图（或称为云纹、叠栅条纹），全 息干涉图，电子散斑干涉图，光栅直接投影的调制条纹图等多种。光学条纹分析技术即利用 干涉图的强度信息，通过信号分析理论与方法试图分析条纹的振幅和相位，从而提取出干 涉图中蕴含的物理信息，例如物体轮廓、形变/位移、应变/应力、机械特性（热、力分布） 等等。由于每种干涉图中条纹的振幅和相位的大小及其变化包含了待研究物体的信息，所 以干涉条纹分析，尤其是相位提取技术在光学检测技术中占有十分重要的地位。
[0003] 通过采用CCD摄像机和数字图像处理技术，现代数字化光学干涉条纹分析技术已 经能够实现相位提取的自动化分析。但是在需要对引起相位变化的瞬态变量进行动态测量 领域中，例如在电子封装的热可靠性分析、焊点失效特性分析、动态测量中的应变/应力分 析、电子器件和机械零部件振动测试等领域，往往只能获取到单幅干涉图，这时采用传统的 移相干涉技术分析干涉条纹就不适用了。因此，在测量物体动态特性领域中，要实现高速干 涉条纹分析，必须对单幅干涉图进行高精度和高速度的相位分析。
[0004] 其中连续小波变换（Continuous Wavelet Transform--CWT)技术，由于时频分析 窗口可变，具有获取局部信息的能力，或者称为多尺度分析能力，在分析含富频率成分信号 时具有得天独厚的优势，逐渐成为光学测量领域中的重要技术之一。但是由于干涉图是二 维信号且干扰噪声较强，因此使用一维连续小波变换（1D-CWT)分析千涉条纹时抗干扰能 力较差。因此已逐渐被二维连续小波变换（2D-CWT)技术所取代。与传统的单幅干涉图分 析技术相比，由于具有更好的抗噪性能，更灵活的时频分析窗，2D-CWT技术在ESPI等光学 条纹干涉分析领域成为近些年来的研究热点。运用2D-CWT技术进行高速千涉条纹分析， 需要解决的关键问题主要是运算速度问题和噪声干扰等。通过改变缩放因子和旋转因子， 2D-CWT具有可变的时频分析窗，从本质上来说是为分析千涉条纹构造了一组参数可变的匹 配滤波器，这一特点极大地改善了干涉场的信噪比，使得相位测量精度得到大幅度的提高， 测量结果更为可靠。因此，2D-CWT技术非常适合用于千涉图的分析上。另外，由于2D-CWT 技术仅需要一幅干涉图即可完成高精度的相位测量，利用这一优势，可以对引起相位变化 的瞬态变量进行测量和深入分析，例如动态测量中的应变/应力分析等。
[0005] 但是，对于单幅512X512像素的干涉图，分析的时间至少要200?400ms，如果提 高分析精度或图像像素数，甚至需要10--30分钟。在动态测量中，例如振动测量，当振动 频率为1kHz,处理一幅千涉图的时间最多只有lms，PC机上实现的算法对高速实时分析场 合明显不适用。因此，运算速度问题极大地限制了基于2D-CWT的光学干涉条纹分析技术在 工业生产领域的推广和应用。目前，对2D-CWT技术的研究普遍集中在非实时测量领域上， 通过PC机对采集后的千涉图进行分析。为了解决高速干涉图分析问题，本发明通过深入研 究2D-CWT的原理和千涉图的本质特征，提出将2D-CWT用于高速测量领域的方法，给出干涉 条纹相位提取快速实现算法及其FPGA实现。
【发明内容】
[0006] 本发明的目的在于提供一种2D-CWT的光学干涉条纹快速处理方法，解决了目前 2D-CWT的光学干涉条纹处理方法速度过慢，不适用于高速实时分析场合的问题。
[0007] 本发明的技术方案按照以下步骤进行：
[0008] 步骤1.通过千涉仪和振动平台，对待测对象进行测试，并利用高速摄像机采集干 涉图，输入至图像数据缓存模块；图像数据缓存模块，用于缓存至少2帧高速采集的干涉图 数据，通过先进先出队列实现，图像数据缓存模块包括干涉图输入接口，用于连接高速摄像 设备以读取干涉图；图像大小设置接口，用于设置干涉图的尺寸以调整缓存区每帧数据的 大小；图像缓存设置接口，用于指定缓存区的帧数和缓存的方式；
[0009] 步骤2.设置2D-CWT参数设置模块，包括母小波设置接口，用于指定参与后续计算 的小波核函数类型；缩放因子和旋转因子设置接口，用于输入参与后续计算的缩放因子序 列的起始值和步进大小，以及旋转因子序列的起始值和步进大小；
[0010] 步骤3 :二维小波核函数频谱模块和二维图像频谱模块对步骤1中图像数据缓存 模块输入的干涉图进行2D-CWT变换；
[0011] 二维小波核函数频谱模块，用于根据步骤2中2D-CWT参数设置模块输入的母小波 类型，计算母小波的频谱；根据输入的缩放因子序列的起始值和步进大小，计算完整的缩放 因子序列；根据输入的旋转因子序列的起始值和步进大小，计算完整的旋转因子序列；
[0012] 二维图像频谱模块，用于对缓存处理后的千涉图进行快速傅立叶变换，得到千涉 图的频谱；
[0013] 在频域中进行2D-CWT变换的公式如下：
[0014] W(u^：0) = ￡3Ι(ω)#*(^^ι))^?2ω
[0015] 其中I，V，u，s，θ，χ，ω分别代表千涉图，小波核函数，位移因子，缩放因子，旋转 因子，二维时域空间坐标和频域空间坐标；r 0是标准2X2旋转矩阵；符号*代表复共辄， 代表傅立叶变换；R2代表二维实数空间，通过傅立叶变换，完成2D-CWT变换；
[0016] 步骤4.小波脊提取模块，用于计算经过步骤3中2D-CWT变换后的小波脊数据； 根据步骤3中缩放因子序列和旋转因子序列，取每一对缩放因子和旋转因子，计算I (X)的 2D - CWT变换时，每一次变换都有一对参数因子（s，Θ)与之相对应，在完成所有（s，Θ)对 的计算后，需要比较每一次2D-CWT变换后的结果，进而得到小波脊序列，以便对高速采样 的干涉条纹图进行相位提取，并在提取过程中实现脊点参数设置模块；
[0017] 步骤5.千涉图分析结果生成模块，用于根据步骤4中的小波脊序列生成干涉图相 位数据和幅度数据；小波脊数据是复数，利用其实部和虚部数据计算出幅度值和相位值，从 而得到与千涉图对应的相位图和幅度图。
[0018] 本发明的有益效果是给出千涉条纹相位提取快速实现算法的FPGA实现。
【专利附图】
【附图说明】
[0019] 图1是本发明一种光学千涉条纹快速处理实现方法计算模块示意图。
【具体实施方式】
[0020] 下面结合【具体实施方式】对本发明进行详细说明。
[0021] 本发明通过千涉仪和高速图像采集系统，在FPGA上实现2D-CWT相位提取算法，以 完成光学干涉条纹相位的高速实时提取。本发明实施例如图 1所示，包括图像数据缓存模 块、2D-CWT参数设置模块、二维图像频谱计算模块、二维小波核函数频谱计算模块、小波脊 提取算法模块和干涉图分析结果生成模块。
[0022] 步骤1.图像数据缓存模块包括千涉图输入接口，用于连接高速摄像设备以读取 干涉图；图像大小设置接口，用于设置千涉图的尺寸以调整缓存区每帧数据的大小；图像 缓存设置接口，用于指定缓存区的帧数和缓存的方式。
[0023] 通过干涉仪和振动平台，对待测对象进行测试，并利用高速摄像机采集干涉图，输 入至图像数据缓存模块；图像数据缓存模块，用于缓存至少2帧髙速采集的干涉图数据，通 过先进先出（FIFO)队列实现。
[0024] 步骤2. 2D-CWT参数设置模块，包括母小波设置接口，用于指定参与后续计算的小 波核函数类型；缩放因子和旋转因子设置接口，用于输入参与后续计算的缩放因子序列的 起始值和步进大小，以及旋转因子序列的起始值和步进大小。
[0025] 步骤3 :二维小波核函数频谱模块和二维图像频谱模块对步骤1中图像数据缓存 模块输入的干涉图进行2D-CWT变换；
[0026] 二维小波核函数频谱模块，用于根据步骤2中2D-CWT参数设置模块输入的母小波 类型，计算母小波的频谱；根据输入的缩放因子序列的起始值和步进大小，计算完整的缩放 因子序列；根据输入的旋转因子序列的起始值和步进大小，计算完整的旋转因子序列。
[0027] 二维图像频谱模块，用于对缓存处理后的干涉图进行快速傅立叶变换，得到干涉 图的频谱。
[0028] 本发明中，在频域中进行2D-CWT变换的公式如下：
[酬 r(u，0): = 现d2_
[0030] 其中I，nr，U，s，θ，χ，ω分别代表干涉图，小波核函数，位移因子，缩放因子，旋转 因子，二维时域空间坐标和频域空间坐标；r 0是标准2X2旋转矩阵；符号*代表复共轭，"" 代表傅立叶变换；R2代表二维实数空间。其中，要计算I (X)的2D-CWT，通过计算I (X)的频 谱％3)和Ψ(χ)的频谱_，然后计算其乘积與纏)的反傅立叶变换即可。因此，在 FPGA中，通过快速傅立叶变换，例如基-2算法或者基-4算法，就能完成2D-CWT变换。
[0031] 步骤4.小波脊提取模块，用于计算经过步骤3中2D-CWT变换后的小波脊数据。 根据步骤3中缩放因子序列和旋转因子序列，取每一对缩放因子和旋转因子，计算I (X)的 2D-CWT变换时，每一次变换都有一对参数因子（s，Θ)与之相对应，在完成所有（s，Θ)对 的计算后，需要比较每一次2D-CWT变换后的结果，进而得到小波脊序列，以便对高速采样 的干涉条纹图进行相位提取，并在提取过程中实现脊点参数设置模块。
[0032]步骤5.干涉图分析结果生成模块，用于根据步骤4中的小波脊序列生成干涉图相 位数据和幅度数据。小波脊数据是复数，利用其实部和虚部数据计算出幅度值和相位值，从 而得到与干涉图对应的相位图和幅度图。
[0033]与现有技术相比，本发明一种2D-CWT的光学干涉条纹快速处理方法的优点为：1、 利用二维快速傅立叶变换（2D-FFT)实现二维连续小波变换（2D-CWT) ;2、给出干涉条纹相 位提取快速实现算法的FPGA实现，在实现中能对2D-CWT的缩放因子和旋转因子等参数的 设置，以及对不同母小波（核函数）的选择；3、为高速实时的干涉相位提取系统提供了有效 的实现手段。
[0034]以上所述仅是对本发明的较佳实施方式而已，并非对本发明作任何形式上的限 制，凡是依据本发明的技术实质对以上实施方式所做的任何简单修改，等同变化与修饰，均 属于本发明技术方案的范围内。
【权利要求】
1. 一种光学千涉条纹快速处理实现方法，其特征在于按照以下步骤进行： 步骤1.通过千涉仪和振动平台，对待测对象进行测试，并利用高速摄像机采集干涉 图，输入至图像数据缓存模块；图像数据缓存模块，用于缓存至少2帧高速采集的干涉图数 据，通过先进先出队列实现，图像数据缓存模块包括千涉图输入接口，用于连接高速摄像设 备以读取千涉图；图像大小设置接口，用于设置千涉图的尺寸以调整缓存区每帧数据的大 小；图像缓存设置接口，用于指定缓存区的帧数和缓存的方式； 步骤2.设置2D-CWT参数设置模块，包括母小波设置接口，用于指定参与后续计算的小 波核函数类型；缩放因子和旋转因子设置接口，用于输入参与后续计算的缩放因子序列的 起始值和步进大小，以及旋转因子序列的起始值和步进大小； 步骤3 :二维小波核函数频谱模块和二维图像频谱模块对步骤1中图像数据缓存模块 输入的干涉图进行2D-CWT变换； 二维小波核函数频谱模块，用于根据步骤2中2D-CWT参数设置模块输入的母小波类 型,计算母小波的频谱；根据输入的缩放因子序列的起始值和步进大小，计算完整的缩放因 子序列；根据输入的旋转因子序列的起始值和步进大小，计算完整的旋转因子序列； 二维图像频谱模块，用于对缓存处理后的干涉图进行快速傅立叶变换，得到干涉图的 频谱； 在频域中进行2D-CWT变换的公式
如下： 其中I，Ψ，u，s，Θ，X，ω分别代表干涉图，小波核函数，位移因子，缩放因子，旋转因 子，二维时域空间坐标和频域空间坐标；r0是标准2X2旋转矩阵；符号*代表复共轭/代 表傅立叶变换；R 2代表二维实数空间，通过傅立叶变换，完成2D-CWT变换； 步骤4.小波脊提取模块，用于计算经过步骤3中2D-CWT变换后的小波脊数据；根据步 骤3中缩放因子序列和旋转因子序列，取每一对缩放因子和旋转因子，计算I (X)的2D-CWT 变换时，每一次变换都有一对参数因子（s，Θ)与之相对应，在完成所有（s，Θ)对的计算 后，需要比较每一次2D-CWT变换后的结果，进而得到小波脊序列，以便对高速采样的干涉 条纹图进行相位提取，并在提取过程中实现脊点参数设置模块； 步骤5.干涉图分析结果生成模块，用于根据步骤4中的小波脊序列生成千涉图相位数 据和幅度数据；小波脊数据是复数，利用其实部和虚部数据计算出幅度值和相位值，从而得 到与干涉图对应的相位图和幅度图。
【文档编号】G01J9/02GK104215340SQ201410442467
【公开日】2014年12月17日 申请日期:2014年9月2日 优先权日:2014年9月2日
【发明者】马峻, 王新强, 黄新, 徐翠锋 申请人:桂林电子科技大学