CN108759709B - 一种适用于表面形貌检测的白光干涉三维重建方法 - Google Patents
一种适用于表面形貌检测的白光干涉三维重建方法 Download PDFInfo
- Publication number
- CN108759709B CN108759709B CN201810212750.3A CN201810212750A CN108759709B CN 108759709 B CN108759709 B CN 108759709B CN 201810212750 A CN201810212750 A CN 201810212750A CN 108759709 B CN108759709 B CN 108759709B
- Authority
- CN
- China
- Prior art keywords
- phase
- optical path
- path difference
- zero optical
- interference
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2441—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明公开了一种适用于表面形貌检测的白光干涉三维重建方法,主要包括强度信息和相位信息,强度信息的获得是通过高斯函数对干涉曲线快速生成其相应的包络曲线,对包络曲线进行函数拟合,连续函数的峰值所在的位置即为零光程差位置。相位信息基于Carré算法解算得到。综合两者结果得到高度信息,三维形貌得以重建。本发明适用于各种样品的表面形貌测量,所求的零光程差位置为连续函数的最大值,比离散采样点的最大值更为准确,并且不受制于移相器采样频率波动的影响。降低了对测试环境的要求,在有外界干扰的情况(振动,光源中心波长波动,相机噪声,被测样品表面反射系数低)下依然可以得到准确的测试结果,其鲁棒性大大提高。
Description
技术领域
本发明涉及一种适用于表面形貌检测的白光干涉三维重建方法,属于表面形貌测量技术领域。
背景技术
表面是指机械加工、表面处理等工艺形成的曲面,直接影响其功能特性和使用,如摩擦磨损、抗腐蚀性能、导电性、导热性等,同时也是评价加工工艺的重要方面。表面形貌测量是制造产品质量表征的重要环节和质量控制与功能可靠保证的前提。随着科学研究和先进制造技术的发展,产品表面形貌与结构越来越多样化,由传统工程表面,发展到越来越多科学表面,如生物表面、IC表面、MEMS表面、能源器件表面、信息存储表面、结构光学表面、功能表面等,由平滑表面发展到结构性表面,由简单几何表面发展到自由曲面,由毫米微米尺度形貌发展到纳米尺度形貌,精度要求也越来越高。表面形貌的迅速发展对其检测技术提出了新的要求。
表面形貌测量在制造领域有着举足轻重的地位,统计表明有10%的制造零部件故障源于表面形貌因素。由于光学测量法的非接触和快速测量特性,使得它在应用中得到了飞速的发展。白光干涉测量法是利用白光的低相干特性,将物体表面三维形貌信息反映到干涉信号上,通过对干涉图像的分析可以恢复出待测样品表面三维形貌。白光干涉测量法最大的优点是测量精度高、速度快,但其对测量环境的要求也较高。
在干涉测量的过程中,通过相移器来改变参考光束和测试光束之间的光程差,同时采集一系列干涉图,对干涉图进行解算得到表面各点的高度信息,三维形貌得以重现。而在实际检测中,由于机械振动,光源波长漂移,移相器的线性误差及相机散粒噪声等环境干扰会引入误差,同时被测样品的表面形貌不同,均会导致测量结果偏差,鲁棒性差。因此设计一种能够抑制干扰,对被测样品无选择性,鲁棒性好的完整重建算法是必要的。
发明内容
本发明的目的是为了解决上述问题,提出一种适用于表面形貌检测的白光干涉三维重建方法,在理想环境和带有干扰的环境中均可解算出准确的结果。
本发明由强度信息p与相位信息组成高度信息以重建表面形貌,其中,带有干涉条纹的图像作为方法处理的原始数据,移相器带动参考镜作单方向等步长的移动,使得干涉条纹划过整个视野,以实现对被测物的干涉扫描;在该过程中,移相器每完成一个步距的移动,相机采集一幅图像并存储;计算机通过提取一系列干涉图像的相同位置的单个像素的像素值获得一组干涉曲线数据;所述的强度信息,通过定位零光程差的位置获得;其过程为:通过高斯函数生成干涉曲线的包络曲线,定位包络曲线的峰值位置作为零光程差的位置;所述的相位信息,由解算零光程差处的光强值得到;其过程为:基于强度信息,取零光程差处附近的四个光强值,由Carré法解算出对应的相位值;结合强度信息和相位信息获得单像素点的高度信息,遍历所有像素点获得对应的高度,从而重建整个表面的三维形貌。
具体的,本发明的一种适用于表面形貌检测的白光干涉三维重建方法,包括以下几个步骤:
步骤一:借助移相器和相机,实现原始数据采集工作。
步骤二:对原始数据进行处理,获取零光程差位置信息。
步骤三:对数据进一步处理,获取相位信息。
步骤四:获取每一点的零光程差位置和相位信息,重现三维形貌。
本发明的优点在于:
本发明能够抑制由机械振动,光源波长漂移,移相器线性误差,相机散粒噪声以及样品表面反射率低引入的误差,对被测样品无选择性,在理想情况下和带有干扰的情况下三维重现结果均准确稳定,鲁棒性好。
附图说明
图1为本发明白光干涉测量系统原理图;
图2为本发明方法流程图;
图3为本发明步骤二和步骤三方法流程图;
图4为本发明重现的三维形貌。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
如图1所示的是白光干涉测量系统的原理示意图,光源经过柯勒照明系统整形后经由分光棱镜分为两束光,一束为照射在被测样品的物光,另一束为照射在参考镜的参考光。两束光由分光棱镜到达相机焦平面处合束并发生干涉。移相器单方向移动带动干涉条纹划过整个图像。相机记录并存储图像作为原始数据。经由方法解算得到每一点的高度值,从而完成整个表面的三维重建。
其方法的关键在于准确的找到零光程差的位置即光强值最大的零光程差的位置以及相应的相位信息。基于高精度和高稳定性的综合考量,本发明提供了一种适用于表面形貌检测的白光干涉三维重建方法,如图2所示。首先通过移相器和相机的配合,记录下原始数据。对原始数据进行处理,提取零光程差位置和相位校正信息。对图像中的每一点对进行上述操作,获得每一点的高度信息,从而重建样品的表面形貌。
步骤一:原始数据采集;
假设移相器步距为s,扫描范围为d。移相器带动参考镜以步距s作单方向移动,相机在每个步距s完成后采集一幅图像。保证干涉条纹划过整个视野。
步骤二:零光程差位置定位;
流程图如图3所示,包含以下两个步骤:
(1)生成包络曲线;
首先对干涉曲线进行低通滤波降噪处理。由下式(1)所示的高斯函数s(τ)调制降噪后的干涉曲线,生成包络曲线。生成的包络曲线与干涉曲线具有相同的峰值位置。
其中:τ是坐标变量,σs为标准差,两者均为实数常数。
这样生成的包络曲线对干涉曲线有一定的容错能力,当干涉曲线发生一定的变型时,依然可以生成的正确的包络曲线。对外界干扰有很强的抑制效果。
(2)获取零光程差位置;
对生成的包络曲线用(2)式所示的二次函数p(x)进行拟合,得到函数表达式,并在区间内求取函数的极值,如(3)式所示,即为零光程差的位置p。
p(x)=A+Bx+Cx2. (2)
p=-B/2C (3)
其中:A、B、C是二次函数的常量系数。
步骤三:相位校正;
如图3所示,基于Carré相移算法对提取到的光强极大值进行相位校正,其基本原理为取零光程差邻域内的4个条纹值进行解算,过程如下:
其中,I′为背景光强值,I1,I2,I3,I4分别为相对于所定位的零光程差领域内的条纹光强值前3帧、前1帧、后1帧、后3帧对应的光强值。α为初始相位,V是条纹对比度。φ为相位校正后的信息。
解得
步骤四:三维形貌重现;
通过对零光程差和相位的解算,获取被测物单点的相对高度:
其中,h为相对高度,N为光强极大值对应的图片帧数,p为步骤二中求解的零光程差位置,λ为光源中心波长,φ为上步解算后的相位信息。
依据步骤二、三对图像中的每一点提取相应的零光程差位置p和相位φ,从而重现样品表面的三维形貌,如图4所示。
本发明方法适用于各种样品的表面形貌测量,对被测样品无选择性。本发明方法所求的零光程差位置为连续函数的最大值,比离散采样点的最大值更为准确,并且不受至于移相器采样频率波动的影响。本发明方法降低了对测试环境的要求,在有外界干扰的情况(振动,光源中心波长波动,相机噪声,被测样品表面反射系数低)下依然可以得到准确的测试结果,其鲁棒性大大提高。
Claims (1)
1.一种适用于表面形貌检测的白光干涉三维重建方法,首先通过移相器和相机的配合,记录样品原始数据,对原始数据进行处理,提取零光程差位置和相位校正信息,对图像中的每一点对进行上述操作,获得每一点的高度信息,重建样品的表面形貌;
其特征在于,具体包括:
步骤一:原始数据采集;
假设移相器步距为s,扫描范围为d,移相器带动参考镜以步距s作单方向移动,相机在每个步距s完成后采集一幅样品图像,确保干涉条纹划过整个视野;
步骤二:零光程差位置定位,包含以下两个步骤:
(1)生成包络曲线;
首先对干涉曲线进行低通滤波降噪处理,通过高斯函数s(τ)调制降噪后的干涉曲线,生成包络曲线,生成的包络曲线与干涉曲线具有相同的峰值位置:
其中:τ是坐标变量,σs为标准差;
(2)获取零光程差位置;
对生成的包络曲线采用二次函数p(x)进行拟合,得到函数表达式,并在区间内求取函数的极值,即为零光程差的位置p;
p(x)=A+Bx+Cx2. (2)
p=-B/2C (3)
其中:A、B、C是二次函数的常量系数;
步骤三:相位校正;
采用Carré相移算法对提取到的光强极大值进行相位校正,得到相位φ;
采用Carré相移算法对提取到的光强极大值进行相位校正,其基本原理为取零光程差邻域内的4个条纹值进行解算,过程如下:
其中,I'为背景光强值,I1,I2,I3,I4分别为相对于所定位的零光程差领域内的条纹光强值前3帧、前1帧、后1帧、后3帧对应的光强值;α为初始相位,V是条纹对比度;φ为相位校正后的信息;
解得
步骤四:三维形貌重现;
通过对零光程差和相位的解算,获取被测物单点的相对高度:
其中,h为相对高度,N为光强极大值对应的图片帧数,λ为光源中心波长;
依据步骤二、三对图像中的每一点提取相应的零光程差位置p和相位φ,最终重现样品表面的三维形貌。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810212750.3A CN108759709B (zh) | 2018-03-15 | 2018-03-15 | 一种适用于表面形貌检测的白光干涉三维重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810212750.3A CN108759709B (zh) | 2018-03-15 | 2018-03-15 | 一种适用于表面形貌检测的白光干涉三维重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108759709A CN108759709A (zh) | 2018-11-06 |
CN108759709B true CN108759709B (zh) | 2020-03-27 |
Family
ID=63980368
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810212750.3A Active CN108759709B (zh) | 2018-03-15 | 2018-03-15 | 一种适用于表面形貌检测的白光干涉三维重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108759709B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109358334B (zh) * | 2018-11-09 | 2023-03-14 | 复旦大学 | 基于部分相干光的纳米位移台精密定位方法 |
CN110223384B (zh) * | 2019-06-27 | 2022-12-27 | 广东工业大学 | 一种白光干涉三维形貌重建方法、装置、系统及存储介质 |
CN110307805B (zh) * | 2019-07-31 | 2020-11-03 | 北京航空航天大学 | 一种用于表面三维形貌测量的白光干涉系统 |
CN111121661B (zh) * | 2019-12-17 | 2021-09-17 | 南京理工大学 | 用于光滑表面形貌测量的窄带非单色光n+1幅相移测试算法 |
CN112525070B (zh) * | 2020-11-18 | 2022-04-01 | 南京理工大学 | 基于非均匀采样校正的抗振动白光干涉测量方法 |
CN113091634B (zh) * | 2021-03-01 | 2022-09-13 | 南京理工大学 | 一种适用于白光扫描干涉的微观形貌快速测量方法 |
CN114088010B (zh) * | 2021-11-18 | 2022-08-19 | 中国科学院长春光学精密机械与物理研究所 | 一种干涉三维形貌解算方法 |
CN114322833B (zh) * | 2021-12-31 | 2022-09-06 | 中国科学院长春光学精密机械与物理研究所 | 基于伪Wigner-Ville分布的白光扫描干涉三维重建方法 |
CN115046469B (zh) * | 2022-05-20 | 2023-05-02 | 浙江大学 | 一种面向光纤白光干涉的干涉条纹包络提取方法 |
CN115325963B (zh) * | 2022-10-11 | 2023-01-13 | 三代光学科技(天津)有限公司 | 一种晶圆表面三维形貌测量装置及其测量方法 |
CN116929572A (zh) * | 2023-05-06 | 2023-10-24 | 深圳市中图仪器股份有限公司 | 白光光源中心波长的测量装置 |
CN116972771B (zh) * | 2023-09-25 | 2023-12-12 | 板石智能科技(深圳)有限公司 | 白光扫描干涉三维重建方法、装置、电子设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59131110A (ja) * | 1983-08-20 | 1984-07-27 | Nitto Seiko Co Ltd | ねじ山成形不良検出装置 |
US5471303A (en) * | 1994-04-29 | 1995-11-28 | Wyko Corporation | Combination of white-light scanning and phase-shifting interferometry for surface profile measurements |
CN1952594A (zh) * | 2005-10-20 | 2007-04-25 | 致茂电子股份有限公司 | 形貌测量方法及其测量装置 |
CN103791853A (zh) * | 2014-01-20 | 2014-05-14 | 天津大学 | 基于彩色条纹信息处理的微结构测量系统及测量方法 |
CN105571517A (zh) * | 2016-01-19 | 2016-05-11 | 北京航空航天大学 | 一种适用于光纤端面检测的改进型相干峰解调方法 |
CN106767496A (zh) * | 2016-11-18 | 2017-05-31 | 中国科学院光电技术研究所 | 一种结合相移干涉与垂直扫描干涉的3d形貌恢复方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9377292B2 (en) * | 2013-08-06 | 2016-06-28 | Zygo Corporation | Interferometry employing refractive index dispersion broadening of interference signals |
-
2018
- 2018-03-15 CN CN201810212750.3A patent/CN108759709B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59131110A (ja) * | 1983-08-20 | 1984-07-27 | Nitto Seiko Co Ltd | ねじ山成形不良検出装置 |
US5471303A (en) * | 1994-04-29 | 1995-11-28 | Wyko Corporation | Combination of white-light scanning and phase-shifting interferometry for surface profile measurements |
CN1952594A (zh) * | 2005-10-20 | 2007-04-25 | 致茂电子股份有限公司 | 形貌测量方法及其测量装置 |
CN103791853A (zh) * | 2014-01-20 | 2014-05-14 | 天津大学 | 基于彩色条纹信息处理的微结构测量系统及测量方法 |
CN105571517A (zh) * | 2016-01-19 | 2016-05-11 | 北京航空航天大学 | 一种适用于光纤端面检测的改进型相干峰解调方法 |
CN106767496A (zh) * | 2016-11-18 | 2017-05-31 | 中国科学院光电技术研究所 | 一种结合相移干涉与垂直扫描干涉的3d形貌恢复方法 |
Non-Patent Citations (3)
Title |
---|
垂直扫描白光干涉表面三维形貌测量系统;戴蓉;《光学技术》;20060731;第32卷(第4期);全文 * |
干涉图零光程差位置的确定方法研究;冯绚;《第34届中国气象学会年会》;20170927;全文 * |
微纳特征结构及表面形貌表征方法的研究;简黎;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20170415(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108759709A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108759709B (zh) | 一种适用于表面形貌检测的白光干涉三维重建方法 | |
CN109163672B (zh) | 一种基于白光干涉零光程差位置拾取算法的微观形貌测量方法 | |
CN110779464B (zh) | 一种时域频域联合分析宽光谱相干测量方法及系统 | |
CN113091634B (zh) | 一种适用于白光扫描干涉的微观形貌快速测量方法 | |
WO2013134966A1 (zh) | 用于非球面测量的波长扫描干涉仪及其应用方法 | |
CN109269438A (zh) | 一种用于多层复杂微纳结构检测的结构光照明显微测量方法 | |
Hu et al. | Microscopic 3D measurement of shiny surfaces based on a multi-frequency phase-shifting scheme | |
CN112525070B (zh) | 基于非均匀采样校正的抗振动白光干涉测量方法 | |
US8553231B2 (en) | Method and apparatus for determining the height of a number of spatial positions on a sample defining a profile of a surface through white light interferometry | |
CN105571517B (zh) | 一种适用于光纤端面检测的改进型相干峰解调方法 | |
CN112097678B (zh) | 基于频率盲估计的多表面面形测量方法 | |
CN102818786A (zh) | 正弦相位调制的并行复频域光学相干层析成像系统和方法 | |
KR20180089910A (ko) | 간섭계의 광학 성능을 최적화하기 위한 방법 및 장치 | |
CN111272099A (zh) | 一种面向航空发动机叶片表面三维形貌的面结构光精密检测系统 | |
CN101881603A (zh) | 横向扫描干涉量测方法与系统 | |
CN113295385B (zh) | 一种镜头内部形变的弱相干光学检测方法与系统 | |
Quan et al. | Determination of surface contour by temporal analysis of shadow moiré fringes | |
CN114440789A (zh) | 旋转体速度、距离和三维形貌同步干涉测量方法及系统 | |
CN109341571A (zh) | 一种双波长同步干涉的表面形貌测量装置和方法 | |
CN111964586A (zh) | 一种基于随机噪声修正的白光干涉信号处理方法 | |
CN114065650B (zh) | 基于深度学习的材料裂纹尖端多尺度应变场测量跟踪方法 | |
KR20200115932A (ko) | Wlsi의 간섭무늬 신호를 이용한 샘플 표면 높이 측정 장치 | |
Fu | Low-frequency vibration measurement by temporal analysis of projected fringe patterns | |
KR20210157969A (ko) | 복합막 단층 측정을 위한 영상 시스템 및 그 방법 | |
CN114279984B (zh) | 用于研究牙科复合材料聚合收缩动力学的非接触式测量系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |