一种新的单幅条纹图的相位解调方法.doc

 一种新的单幅条纹图的相位解调方法 3 杨福俊1 何小元1 王蕴珊2 (1 东南大学 工程力学系 , 江苏 南京 210096) (2 山东大学 物理与微电子学院 , 山东 济南 250061) 摘 要 基于非载波条纹图的条纹灰度分布的极值准确定位 ,提出了一种新的由单幅条纹图解调 相位方法 1 该方法首先准确定位每个条纹灰度的极值包括最大和最小 ,再基于极值图将条纹的强 度值分布线性地变换到21 和 1 之间 ,最后利用 a rco s 反余弦取代传统的 at a n2 反正切算法求出相 位 ,结合极值图就能快速准确地解调出含有封闭条纹的条纹图相位 1 本方法仅用一幅条纹图 ,不需 要传统的条纹定级和对分数级条纹插值与拟合 ,自动判别条纹相位变化的波谷 ,准确获得条纹相位 分布 1 这一方法使光测法研究动态和瞬态问题成为可能 ,而不需引入载波法调制条纹 1 文中通过 实例说明了本方法的处理过程 1 关键词 相位解调 ; 条纹分析 ; 图像处理 ; 自适应加权滤波 中图分类号 文献标识码 O348 . 1 ; T P301 A 条纹相位的 (至少是局域的) 连续性以及相应的条纹 图灰度分布也应该是连续分布的假设 ,利用迭代法 求解出相位 1 由于条纹相位和条纹图的连续性的要 求 ,这样在迭代求解相位的过程中 ,不规则的随机噪 声同时也被滤除或减少其对解调出的相位值的影 响 ;类似还有二维非线性回归法[ 12 ] , 基于遗传算法 的多项式参数法[ 13 ,14 ] ,多层神经网络法[ 15 ] , 这些方 法都类似于相位跟踪法 ,但后面的这几种方法都要 求条纹图的信噪比比较高 ,采用迭代法求解相位 ,如 果噪声影响大 ,将导致最终结果不能收敛到实际值 , 且求解时间都很长 ,一般中小图像采用 Pe nti u m III 以上的计算机处理需要一小时左右 1 (二) 座标变换 法[ 16 ] ,采用座标变换将直角坐标下的封闭条纹转换 到极座标下的非封闭条纹 , 再利用 Fo urie r 变换法 解调出相位 ,再将解调结果的反演回原坐标 1 这种 方法虽然处理很快 ,但只能处理含有一个中心的封 闭条纹 ,而且封闭条纹中心处的解调相位误差大 1 本文针对上述现有方法的不足 ,通过非载波条 纹图的条纹灰度分布的极值准确定位 ,提出了一种 新的由单幅条纹图解调相位方法 ,文中同时应用局 域灰度梯度信息对条纹进行加权滤波处理 ,获得了 较理想的结果 1 0 引言 除通过图像相关运算等方法外 ,光测力学的数 据 多是通过采集处理光学干涉条纹图来获得结果的1 条纹图的条纹相位与诸如三维形貌 、位移 、应变 、振 动以及温度梯度等待测物理量密切相关 1 二维条纹 图的条纹相位解调方法有相移法[ 1 ] 、基于载波条纹 调制的 Fo urier 变换法[ 2～5 ] 、同步探测法及 锁相 法 等[ 6 ,7 ]1 这些方法中 ,采用相移法解调出的相位精度 最高 ,但要求在相移过程中保持测试系统的机械稳 定 ,避免大气扰动 ,同时一般要求采 3 幅以上的条纹 图 1 显见相移法不能用于研究瞬态或动态问题 1 对 于不含有封闭条纹的非载波条纹图 , Fo urier 变换法 通过分 离 提 取 单 侧 频 谱 实 现 相 位 解 调[ 3 ] 1 另 外 , Fo urie r 变换法通过载波调制条纹 ,只需一幅条纹图 就可以解调出相位 1 光学计量中的所谓载波条纹通 常是加进有一定斜率的线性载波调制非单调的条纹 相位分布 ,这样可以用 Fo urie r 变换法研究动态或 瞬态问题 1 然而据 N yqui st 采样定律 ,引进的载波 条纹的频率必须是被调制条纹最高频率的两倍 ,因 此除非被调制条纹的频谱很窄 ,空间载波法才能使 用 1 基于上述原因 ,国内外学者开始致力于单幅非 载波条纹图 ,尤其是含有封闭条纹的相位解调研究1 目前已开发研究出解调单幅条纹相位的方法有两 类 : ( 一 ) 规 则 化 相 位 跟 踪 法 ( Regula rized p ha se2 t racki ng) 和 Re gula rizatio n 法[ 8～11 ] , 两 者 都 是 基 于 1 1 . 1 原理 条纹灰度值的分布与条纹相位关系 对于大小为 M ×N , 带有噪声的干涉条纹图上 任意点 r ( i , j) 灰度值分布 I ( r) 一般可以写成 I ( r) = a ( r) + b( r) co s ( r) 3 自然科学基金资助项目 ( 、国家高技术研究发展 计划“863”资助项目 (2002A A404140) Tel :025 E