CN100451535C - 移相干涉图像的信息处理方法 - Google Patents
移相干涉图像的信息处理方法 Download PDFInfo
- Publication number
- CN100451535C CN100451535C CNB2006101474829A CN200610147482A CN100451535C CN 100451535 C CN100451535 C CN 100451535C CN B2006101474829 A CNB2006101474829 A CN B2006101474829A CN 200610147482 A CN200610147482 A CN 200610147482A CN 100451535 C CN100451535 C CN 100451535C
- Authority
- CN
- China
- Prior art keywords
- phase
- point
- picture element
- value
- fiduciary level
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 230000010365 information processing Effects 0.000 title claims abstract description 10
- 238000003672 processing method Methods 0.000 title claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000000605 extraction Methods 0.000 claims abstract description 16
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 10
- 230000010363 phase shift Effects 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 11
- 238000004458 analytical method Methods 0.000 abstract description 7
- 230000003287 optical effect Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 210000001747 pupil Anatomy 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000001900 immune effect Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Landscapes
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明涉及一种用于移相干涉的信息自动化处理系统及其处理方法。其信息处理系统至少包括干涉系统,图像采集系统,计算机分析处理系统和输出系统。其信息处理的方法步骤:首先对获得的数字干涉图做向量滤波;用对误差不敏感的组合相位提取方法提取相位,得到包裹相位;再对上述包裹相位应用可靠度相位展开方法,恢复相位信息;最后通过获取的相位分布,可进一步得到待测的参量。本发明充分考虑了误差对测量结果的影响,改善了测量精确度。
Description
技术领域
本发明属于光学测量领域,特别是一种涉及激光干涉仪的移相干涉图像的信息处理方法。
背景技术
光学干涉计量是以光的波长为计量单位的一种高精度、高灵敏度的计量测试方法,它在光学测量中有广泛的应用。干涉计量检测的结果是通过干涉图的识别与分析实现的,它具有多信息量,非接触,高精度等特点。
从双光束干涉原理知道,干涉条纹是干涉场中光程差相同点的轨迹。干涉场强度的空间变化率主要取决于相位的变化。在移相干涉仪面形测量中,干涉条纹分析就是要从一帧或者多帧干涉图中恢复出波面的原始特征。
采用主动相位调制得到的光强分布I可以通写为:
I(x,y,t)=a(x,y)+b(x,y)cos[φ(x,y)+δ(t)](1)
式中:(x,y)为干涉场坐标,a是背景光强,b是干涉图调制度,φ是待测相位,δ(t)是相位调制项。
现有技术[1]CN 1584533A的相位提取,采用传统的四步法,即连续帧之间相位步距为π/2,选取δ(t)分别在0,π/2,π和3π/2处的光强为I1,I2,I3和I4,则所求相位解为:
相位提取是移相干涉测量中的重要技术,相位提取的精度直接决定了仪器的精度,但相位的值很容易受到系统自身移相器的误差和不稳定因素的影响,如移相器的误差和系统的振动、温度和气流等不稳定因素,系统误差是不可避免的。移相上存在的非线性误差及其它测量误差,在测量信号中以谐波形式存在。现有技术[1]的相位提取法,对误差都是比较敏感的。在移相中的轻微的误差将导致最后的相位测量结果的很大变化。现有技术[2](Danie Malacara,Interferogram Analysis forOptical Testing,Second Edition,Taylor & Francis Group,CRC Press,2005:275~282)中的四步法和五步法相位提取都是对失谐误差很敏感的,造成相位提取的精确度低。
由(2)式可知,得到的相位信息均以反正切表示。由于反正切函数的固有性质,干涉相位图像中的相位值被限制在[-π,π)之间,即相位值是在限定范围内呈阶跃的不连续的分布。需要将多个截断相位的区域拼接展开成连续相位,以反映被测面形真实的相位变化。在已公开的现有技术[3]CN 1673666A中涉及到的信息处理方法,没有全面考虑残差点对结果的影响,其测量精度仍有待进一步提高。
发明内容
本发明的目的在于克服上述现有技术的不足,提供一种移相干涉图像的信息处理系统及其处理方法,本发明应充分考虑了误差对测量结果的影响,改善测量精确度,以实现高精度的检测。
本发明的技术方案如下:
一种移相干涉图像信息处理系统,由一台移相干涉仪、CCD摄像机、图像采集卡,图像监视器、具有显示器和打印机的计算机组成,该移相干涉仪对待测元件检测产生的干涉图像由CCD摄像机摄取,经图像采集卡模数转换后,送图像监视器显示,或送计算机进行数据处理。
一种移相干涉图像信息处理方法,该方法的步骤如下:
(1)调整移相干涉仪对待测元件检测产生的干涉图像用CCD摄像机摄取N帧干涉图像,经图像采集卡模数转换后,送入计算机进行数据处理;
(2)根据测量精度的需要,进行向量滤波,以减少系统随机噪声对干涉图的影响,提高信噪比;
(3)用组合相位提取法在N帧干涉图像的光强分布图提取相位,形成包裹相位图;
(4)对所述的包裹相位图用可靠度相位展开法展开相位;
(5)计算相位梯度分布值;在相位分布的基础上进行泽尼克多项式拟合,或作其它进一步的计算。
所述的向量滤波的具体操作过程如下:
(2)分别对余弦和正弦分量进行平均滤波处理:
式中:T1和T2分别是滤波后的余弦和正弦值,m为滤波窗口的尺寸;
(3)通过逆映射得到滤波后的干涉图像相位值φ:
φ=arctg(T1/T2)。
所述的滤波窗口的尺寸的选取原则是:如果噪声干扰比较严重,可以选择大一点的窗口,否则要尽量选择小一些的窗口以便保留更多的原始图像的信息,通常滤波窗口的尺寸为3×3。
所述的组合相位提取方法是:
对连续以π/2为步距的五帧图,对应地把相位π/4,3π/4,5π/4,7π/4,9π/4处的光强为I1,I2,I3,I4,I5,相位提取算法为:
然后,通过对上式的反正切变换,求得包裹的相位;
所述的可靠度相位展开法的具体步骤如下:
(1)计算每个象素点的可靠度:
对N帧移相干涉图,第n帧的干涉图中坐标为(x,y)的像素点(x,y)的光强In,则包裹相位图中象素点(x,y)的可靠度为:
(2)选可靠度最高的象素点为起始点,该起始点的相邻的象素点按可靠度高低排序;
(3)在上述相邻的象素点中选择可靠度较高的象素点,以起始点为参照点展开该象素点的相位:
相位展开时遵循如下原则:若该象素点与参照点的包裹相位之差在[-π,π)内,该点的展开相位值为其包裹相位值;若该点与参照点的包裹相位差超过π,认为有2π正跳变,则该点的展开相位值在原来基础上减2π;反之,若该点与参照点的包裹相位差小于-π,则该点的展开相位值在原来基础上加2π;
(4)再以已经展开相位的象素点和起始点为参照区,将该参照区周围的还未进行相位展开的邻近像素点按可靠度高低排队,选取可靠度较高的象素点并以参照区中与该象素点最相邻的象素点的相位为参照,按第(3)步的原则进行相位展开;
(5)重复上述第(4)步,直至所有象素点的相位全都展开为止。
本发明的效果如下:
根据信号提取在频域上的评价函数方法,其具体方法参见在先技术(DanieMalacara,Interferogram Analysi s for Optical Testing,Second Edition,Taylor& Francis Group,CRC Press,2005:190~224)。
本发明采用的组合相位提取方法的形成是这样的:
对连续以π/2为步距获得五帧图,对应地把相位π/4,3π/4,5π/4,7π/4,9π/4处的光强分两组A(I1,I2,I3,I4)和B(I2,I3,I4,I5),
我们实施一种新的组合,该组合是将上述两式的分子相加作分子,分母相加作分母而构成:
通过上式(4)的反正切变换求得包裹的相位φ。
对本发明组合相位提取方法与已有的传统的相位提取方法仿真分析,结果如图2所示,图中横坐标是归一化频率的相位误差考察,纵坐标为所提取的相位误差。图中分别表示了现有技术[1]和现有技术[2]相位提取的四步法,五步法相位提取与本发明的相位提取方法的误差敏感度对比图。由图可见,本发明组合相位提取方法对失谐误差不敏感,有很好的免疫特性,应用本发明组合相位提取方法提取相位有助于更好提高检测的准确率。
本发明应用了向量滤波方法,实践表明向量滤波方法在消除噪声、减少残差点的同时还保持了相位图的清晰度和图像的分辨率。
本发明的可靠度相位展开法,可有效避免误差的传递,有利于提高检测的精度。模块化的软件结构设计,易于编程实现。
附图说明
图1是本发明移相干涉信息处理系统组成示意图
图2是本发明相位提取方法与在先技术方法误差敏感度对比图
具体实施方式
下面结合附图和实施例对本发明作进一步说明,但不应以此限制本发明的保护范围。
先请参阅图1,图1是本发明移相干涉信息处理系统的组成示意图。由图可见,本发明移相干涉图像信息处理系统,由一台移相干涉仪2、CCD摄像机3、图像采集卡4,图像监视器5、具有显示器6和打印机8的计算机7组成,该移相干涉仪2对待测元件1检测产生的干涉图像由CCD摄像机3摄取,经图像采集卡4模数转换后,送图像监视器5显示,送计算机7进行数据处理。
所述移相干涉仪是能实现移相的菲索干涉仪,对待测元件1进行检测。用图像采集卡4将电荷藕合器件(Charge-Coupled Devices,CCD)摄像机3上的光学干涉图像的模拟信号,经过模数转换,使干涉图像数字化,随后便可借助计算机7对所采集到的干涉图像进行各种必要的处理和分析。在该系统中还包括图像监视器5。在专用的图像监视器上显示干涉图像,供检测者观察,调节干涉图像使条纹达到测量要求的最佳状态。
所述的计算机7是分析处理系统,是整个检测系统的控制中枢,计算机7除了用于完成对干涉图像进行分析与计算的工作外,还担负着对整个检测系统实时控制任务,如对CCD图像摄取和图像采集动作的控制等。
用于交互键盘和显示器6和打印机8是输出系统,把测量结果进行输出的设备。当系统完成整个检测任务之后,所得到有关物理信息,以文字、图形的形式输出到相应设备。
本发明移相干涉图像信息处理方法,其特征在于该方法的步骤如下:
(1)通过CCD摄像机3获得移相干涉仪2上的图像,经图像采集卡4模数转换后送如计算机7处理;
(2)对获得的数字干涉图应用向量滤波法做预处理。该向量滤波方法的实现步骤是:对干涉图中的每一个象素点都取相应的余弦和正弦值,即将象素的光强值映射到一个二维向量空间,两个分量分别为余弦和正弦分量,得到变换的余弦和正弦图像。再选择经典的平均滤波法分别对余弦和正弦图像做滤波处理。最后通过逆映射便可得到滤波后的干涉图像光强值。该向量滤波法在消除噪声、减少残差点的同时还保持了相位图的清晰度和图像的分辨率,更有利与提高检测精度。
(3)在干涉图光强分布图中应用本发明对误差不敏感的组合相位提取方法提取相位;
(4)对获得的上述包裹相位采用可靠度相位展开法进行展开,恢复原来的相位信息。
(5)计算相位梯度分布值;在相位分布的基础上进行泽尼克多项式拟合,或作其它进一步的计算。
获取相位分布函数φ(x,y),并以此为基础得到波前梯度分布值。
波前梯度在位相分布函数φ(x,y)的基础上得到,基本的方法是如下
一维波前梯度
二维波前梯度
从梯度的数学定义出发可以得到波前相位梯度的计算式,梯度表示的是一个矢量场沿着一个方向的变化。通过计算得到的波前梯度分布图形上,可以直观的得到被检测元件的波前变化情况,梯度大,分布密集,表示波前位相畸变大,反之,表示位相畸变小。
基于双光束干涉原理,波前函数W(x,y)与相位的关系式为 W(ρ,θ)是在极坐标下的波前表示。波前函数一定可以表达成为一个完备的基底函数的线性组合,或一个线性无关的基底函数的线性组合。这里应用泽尼克多项式进行线形拟合:
式中:Ar为线性组合系数,Ur为泽尼克多项式,L为选择的多项式级数。
其它检测结果:
光学元件质量用在整个瞳面上的峰谷值和瞳孔面内光学位相的均方根偏差来说明。峰谷差P-V值:即通过波前干涉所反映出来的被检测元件表面最高突出部分(峰)与最低凹陷部分(谷)之间的差值:
PV=max(Wi)-min(Wi) (11)
均方根偏差用RMS表示,对采样点总数为N的均方根偏差值,可由下式得到:
本发明充分考虑了误差对测量结果的影响,改善了测量精确度。
Claims (2)
1、一种移相干涉图像信息处理方法,其特征在于该方法包括如下步骤:
(1)移相干涉仪(2)对待测元件(1)检测产生的干涉图像用CCD摄像机(3)摄取N帧干涉图像,经图像采集卡(4)模数转换后,送入计算机(7)进行数据处理;
(2)用组合相位提取法在所述的N帧干涉图像提取相位,形成包裹相位图:
对连续以π/2为步距获得五帧图,对应地把相位π/4,3π/4,5π/4,7π/4,9π/4处的光强为I1,I2,I3,I4,I5,用下列公式计算:
然后,通过对上式的反正切变换,求得包裹的相位;
(3)根据测量精度的需要,进行向量滤波:
②分别对余弦和正弦分量进行平均滤波处理:
式中:T1和T2分别是滤波后的余弦和正弦值,m为滤波窗口的尺寸;
(4)对所述的包裹相位图用可靠度相位展开法展开相位;
(5)计算相位梯度分布值;在相位分布的基础上进行泽尼克多项式拟合。
2、根据权利要求1所述的移相干涉图像信息处理方法,其特征在于所述的可靠度相位展开法的具体步骤如下:
(1)计算每个象素点的可靠度:
对N帧移相干涉图,第n帧的干涉图中坐标为(x,y)的像素点(x,y)的光强In,则包裹相位图中象素点(x,y)的可靠度为:
(2)选可靠度最高的象素点为起始点,该起始点的相邻的象素点按可靠度高低排序;
(3)在上述相邻的象素点中选择可靠度较高的象素点,以起始点为参照点展开该象素点的相位:
相位展开时遵循如下原则:若该象素点与参照点的包裹相位之差在[-π,π)内,该点的展开相位值为其包裹相位值;若该点与参照点的包裹相位差超过π,认为有2π正跳变,则该点的展开相位值在原来基础上减2π;反之,若该点与参照点的包裹相位差小于-π,则该点的展开相位值在原来基础上加2π;
(4)再以已经展开相位的象素点和起始点为参照区,将该参照区周围的还未进行相位展开的邻近像素点按可靠度高低排队,选取可靠度较高的象素点并以参照区中与该象素点最相邻的象素点的相位为参照,按第(3)步的原则进行相位展开;
(5)重复上述第(4)步,直至所有象素点的相位全都展开为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101474829A CN100451535C (zh) | 2006-12-19 | 2006-12-19 | 移相干涉图像的信息处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101474829A CN100451535C (zh) | 2006-12-19 | 2006-12-19 | 移相干涉图像的信息处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1975321A CN1975321A (zh) | 2007-06-06 |
CN100451535C true CN100451535C (zh) | 2009-01-14 |
Family
ID=38125559
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006101474829A Expired - Fee Related CN100451535C (zh) | 2006-12-19 | 2006-12-19 | 移相干涉图像的信息处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100451535C (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101877124A (zh) * | 2009-04-28 | 2010-11-03 | 北京捷科惠康科技有限公司 | 一种医学图像的滤波方法及系统 |
CN102175332B (zh) * | 2011-01-21 | 2012-06-06 | 南京理工大学 | 一种从含有移相误差的干涉图中恢复相位的方法 |
CN103778650A (zh) * | 2013-10-25 | 2014-05-07 | 南京农业大学 | 主成分分析相位恢复算法 |
CN104330027B (zh) * | 2014-11-18 | 2017-04-12 | 哈尔滨工业大学 | 一种基于误差互补修正的移相干涉测量位相提取方法 |
CN106482633B (zh) * | 2015-08-24 | 2019-01-18 | 南京理工大学 | 一种基于π/4相移的多光束干涉相位提取方法 |
CN106649922B (zh) * | 2016-09-18 | 2020-05-08 | 中国科学院长春光学精密机械与物理研究所 | 前处理接口程序的光机集成分析方法及镜面面形优化方法 |
CN107680063B (zh) * | 2017-10-23 | 2021-03-12 | 西华大学 | 一种直接数字化图像的增强方法 |
CN108170636B (zh) * | 2018-01-11 | 2021-04-06 | 苏州科技大学 | 一种从线性载频移相干涉图中恢复出真实相位的方法 |
CN111461295B (zh) * | 2020-03-20 | 2022-08-16 | 南京理工大学 | 基于多尺度生成对抗神经网络的单帧条纹分析方法 |
CN113048876A (zh) * | 2021-03-16 | 2021-06-29 | 上海乾曜光学科技有限公司 | 一种应用于移相式激光干涉仪的振动检测预处理方法 |
CN115183697B (zh) * | 2022-07-18 | 2024-07-16 | 西安交通大学 | 一种基于干涉光强信息的相位提取方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1096873A (zh) * | 1993-06-26 | 1994-12-28 | 南京理工大学 | 高精度大孔径移相式数字平面干涉仪 |
WO2004072629A1 (en) * | 2003-02-17 | 2004-08-26 | Nanyang Technological University | System and method for inspection of silicon wafers |
CN1584533A (zh) * | 2004-06-02 | 2005-02-23 | 中国科学院上海光学精密机械研究所 | 光纤连接器端面几何参数测量仪 |
-
2006
- 2006-12-19 CN CNB2006101474829A patent/CN100451535C/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1096873A (zh) * | 1993-06-26 | 1994-12-28 | 南京理工大学 | 高精度大孔径移相式数字平面干涉仪 |
WO2004072629A1 (en) * | 2003-02-17 | 2004-08-26 | Nanyang Technological University | System and method for inspection of silicon wafers |
CN1584533A (zh) * | 2004-06-02 | 2005-02-23 | 中国科学院上海光学精密机械研究所 | 光纤连接器端面几何参数测量仪 |
Non-Patent Citations (2)
Title |
---|
在模拟干涉相位图中噪声的滤波研究. 宋芳,苏显渝.四川大学学报(自然科学版),第43卷第1期. 2006 * |
数字相移干涉术中图像处理和波前再现误差的分析及校正算法的研究. 刘青.山东大学博士学位论文. 2005 * |
Also Published As
Publication number | Publication date |
---|---|
CN1975321A (zh) | 2007-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100451535C (zh) | 移相干涉图像的信息处理方法 | |
Huang et al. | Review of phase measuring deflectometry | |
CN109253708B (zh) | 一种基于深度学习的条纹投影时间相位展开方法 | |
Zappa et al. | Comparison of eight unwrapping algorithms applied to Fourier-transform profilometry | |
Hu et al. | Absolute three-dimensional micro surface profile measurement based on a Greenough-type stereomicroscope | |
CN108955575B (zh) | 一种基于单幅干涉条纹高精度地恢复波面的方法 | |
CN109870129A (zh) | 一种基于相位偏折原理的晶圆表面粗糙度检测装置 | |
Xu et al. | Enhancement of measurement accuracy of optical stereo deflectometry based on imaging model analysis | |
CN111307063A (zh) | 一种消除单幅干涉条纹波面恢复中的符号模糊问题的方法 | |
CN111537203A (zh) | 一种利用相位测量偏折术的光学透镜波前测量方法 | |
Nguyen et al. | Real-time 3D measurement of freeform surfaces by dynamic deflectometry based on diagonal spatial carrier-frequency pattern projection | |
Xiong et al. | Evaluation of the pitch deviation of a linear scale based on a self-calibration method with a Fizeau interferometer | |
CN111561877A (zh) | 一种基于点衍射干涉仪的可变分辨率相位解包裹方法 | |
CN106482664B (zh) | 一种基于圆载频莫尔条纹理论的合成波长相位提取方法 | |
Chang et al. | Phase-measuring profilometry using sinusoidal grating | |
CN116559179B (zh) | 一种反光表面形貌与缺陷检测方法及其系统 | |
Wang | Development and application of computer-aided fringe analysis | |
Miao et al. | Analysis of phase distortion in phase-shifted fringe projection | |
CN108489422B (zh) | 一种变频相移最小二乘迭代叠加面形分离的方法 | |
CN103398666B (zh) | 一种用于双层周期性微结构的层间错位测试方法 | |
Bertani et al. | High-resolution optical topography applied to ancient painting diagnostics | |
Zhu et al. | Analysis and optimization of shielding in measuring helical gear tooth flank by laser interferometry | |
Wang et al. | High-precision technique for in-situ testing of the PZT scanner based on fringe analysis | |
Yang et al. | A radial-shearing interference system of testing laser-pulse wavefront distortion and the original wavefront reconstructing | |
Morimoto et al. | Theory and application of sampling moiré method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090114 Termination date: 20111219 |