CN104034277B - 一种双波长微纳结构相位测量方法 - Google Patents
一种双波长微纳结构相位测量方法 Download PDFInfo
- Publication number
- CN104034277B CN104034277B CN201410139316.9A CN201410139316A CN104034277B CN 104034277 B CN104034277 B CN 104034277B CN 201410139316 A CN201410139316 A CN 201410139316A CN 104034277 B CN104034277 B CN 104034277B
- Authority
- CN
- China
- Prior art keywords
- wavelength
- phase
- measurement
- height
- nano structure
- 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
Landscapes
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明公开一种双波长微纳结构相位测量方法,通过使用两个不同波长的激光器作为光源,分别获取相位包裹图,利用同一位置(像素)上的两个波长已知的固定关系来确定该位置(像素)的相位跳变系数,重构出物体的相位。由于物体不同位置(像素)的相位信息是独立获得的,不依赖相邻像素之间的关系,可以用于测量具有孤立非连续物体的相位信息,可以有效避免解包裹算法中由于路径依赖的原因产生误差积累传播的现象,不仅可用于连续光滑物体,也可以用于孤立非连续(存在突变)物体的相位测量,同时,本发明方法大大拓展了有效测量突变物体的量程。本发明测量方法简单、测量范围大、精度高、测量实时高效快捷。
Description
技术领域
本发明涉及一种应用于微纳结构物体形貌的干涉测量方法,属于光电检测技术的领域。
背景技术
干涉测量方法是获取物体形貌(相位)微纳结构的常用方法,通过携带了物体相位信息的物光与参考光之间形成的干涉图,使用相位提取算法,可以获取物体的三维信息。由于使用激光作为计量手段,可以实现纳米精度的测量。通常使用的干涉测量方法,结合使用相移技术,其测量精度可以达到百分之一波长的精度。然而由于所获取的相位使用反正切函数来表达,这类方法获取的相位值分布在[-π,π]的相位包裹中,当待测物体的高度超过一个波长时,需要对包裹相位进行解包裹以获取真实的相位分布。从20世纪70年代开始,国内外经过大量的研究,已经出现了不少于40种的相位解包裹算法,这些算法主要分为两大类:最小范数法和路径跟踪算法,由于这两类算法依靠相邻像素之间的关系来实现解包裹,为了满足像素间的采样定律,待测物体必须是光滑连续的,使得这些方法难于应用于工程中常见的具有孤立非连续(存在突变)物体的测量。本发明方法给出了一种可以克服现有方法不足的新方法,通过使用两个不同波长的激光器作为光源,分别获取相位包裹图,利用同一位置(像素)上的两个波长已知的固定关系来确定该位置(像素)的相位跳变系数,重构出物体的相位。由于物体不同位置(像素)的相位信息是独立获得的,不依赖相邻像素之间的关系,可以用于测量具有孤立非连续物体的相位信息,可以有效避免解包裹算法中由于路径依赖的原因产生误差积累传播的现象,不仅可用于连续光滑物体,也可以用于孤立非连续(存在突变)物体的相位测量,同时,本发明方法大大拓展了有效测量突变物体的量程。
发明内容
针对目前各种干涉测量解相位包裹方法的局限,本发明提供一种可以测量孤立非连续(存在突变)物体的干涉测量方法,扩大了干涉测量的量程,突破了传统相位解包裹的限制。
在干涉测量中,待测物体的高度超过测量波长(量程)时,可以表示为:
h(x,y)=n(x,y)λ+φ(x,y) (1)
其中λ为测量波长;n(x,y)是当待测物体高度超出波长时跳变的次数,为非负整数;为高度包裹,为包裹相位,在测量中,我们可由干涉光路或者干涉仪获取的干涉图通过相位提取算法获得包裹相位因此,要获得完整的待测物体的高度信息,还需要确定跳变系数n(x,y)。
从(1)式我们可以知道当使用两个波长来测量同一待测物体时,有:
我们将(2)中上下两式相减时得到的高度包裹差分z(x,y)为:
z(x,y)=φ1(x,y)-φ2(x,y)=n2(x,y)λ2-n1(x,y)λ1 (3)
由于测量使用的波长是已知的,我们可以确定波长差分函数:
δ(n1,n2)=n2λ2-n1λ1,(n1,n2为非负整数) (4)
对比(3)、(4)式,我们可以发现,高度包裹差分z(x,y)刚好与波长差分函数δ(n1,n2)一致,据此,我们可以从波长差分函数中找到实际测量时的跳变系数,从而确定待测物体的高度。
本发明中,n1或n2是由δ(n1,n2)决定的,也就是说,本发明方法成立的条件是:只要能保证有且仅有一对n1,n2,使得δ(n1,n2)的值是唯一的,就能根据δ(n1,n2)的值确定n1或n2的值。以下证明该方法成立的条件:假设存在另一对非负整数m1,m2,使δ(m1,m2)=δ(n1,n2),即:
m2λ2-m1λ1=n2λ2-n1λ1 (5)
从(5)式我们可以得到:
考虑λ1,λ2存在最大公约数k,可得:
当m1,m2存在时,不能保证有且仅有一对n1,n2,使得δ(n1,n2)的值是唯一的,此时本发明失效,即超出本发明的量程,将出现跳变模糊。为满足本发明方法成立的条件,必须保证m1,m2不存在,此时有:
考虑高度包裹φ1,φ2的最大值为λ1,λ2,由此n1,n2确定的高度为:
由(10)中确定的高度即为本发明的有效量程,在这个有效量程内,可以保证有且仅有一对n1,n2,使得δ(n1,n2)的值是唯一的。从(10)式中也可以看出,有效量程与两个波长的最大公约数成反比,如果在选择激光波长的时候,使两波长的最大公约数为1,就能获得最大量程。
实现本发明方法包括以下步骤:
步骤一:使用干涉光路或者干涉仪分别获取两种波长λ1、λ2作为计量波长时的干涉条纹图;使用已知的相位提取算法获取两种波长的包裹相位图并转换成高度包裹图φ1(x,y)、φ2(x,y),其中
步骤二:获取高度包裹差分图z(x,y)=φ1(x,y)-φ2(x,y);
步骤三:确定波长差分函数δ(n1,n2)=n2λ2-n1λ1,其中,n1,n2为非负整数;
步骤四:将步骤二获得的高度包裹差分图中的点z(x,y)与波长差分δ(n1,n2)逐一比较,找出最相近的波长差分函数值δ(n1,n2),获取n1(x,y)或n2(x,y);
步骤五:则物体的高度可以表示为h(x,y)=n1(x,y)λ1+φ1(x,y)或者h(x,y)=n2(x,y)λ2+φ2(x,y);
本发明利用不同波长的高度包裹差分与已知波长差分之间一一对应的关系,解决了干涉测量中相位跳变的问题。与现有技术相比,本发明有以下的优点:
(1)确定跳变系数不需要借鉴相邻像素之间的关系,克服了传统相位解包裹算法只能用于测量光滑连续的物体,可以实现孤立非连续(存在突变)物体的微结构三维轮廓测量。
(2)大大扩展了干涉测量的量程,可以解决当物体存在数百倍波长突变时的相位测量问题。
(3)提出了在双波长干涉测量中有效量程的限制条件和确定方法,可以根据物体突变的范围灵活选择波长以改变有效量程,实现精确测量。
(4)本发明方法的实施过程简单,重构物体微纳形貌的速度快,精度高。
附图说明
图1为本发明方法采用的双波长相移同轴迈克尔逊干涉测量系统示意图。
图2为本发明方法在使用常见相位提取方法获得相位后,重构物体相位(高度)的流程图。
图3a为使用532纳米波长为光源时获取到的凹面镜的高度包裹灰度图。
图3b为使用632.8纳米波长为光源时获取到的凹面镜的高度包裹灰度图。
图3c为将两波长高度包裹相减后得到的两波长高度包裹差分灰度图。
图3d为使用本发明方法获取的凹面镜的高度分布图。
图3e为使用传统解包裹方法获取的凹面镜的高度分布图。
图4a为随机生成的孤立待测方柱图。
图4b为使用532纳米波长为光源时获取到的方柱的高度包裹灰度图。
图4c为使用632.8纳米波长为光源时获取到的方柱的高度包裹灰度图。
图4d为将两波长高度包裹相减后得到的两波长高度包裹差分灰度图。
图4e为使用本发明方法获取的孤立方柱的高度分布图。
图5a为模拟生成的待测斜坡图。
图5b为使用532纳米波长为光源时获取到的斜坡的高度包裹灰度图。
图5c为使用632.8纳米波长为光源时获取到的斜坡的高度包裹灰度图。
图5d为将两波长高度包裹相减后得到的两波长高度包裹差分灰度图。
图5e为使用本发明方法获取的斜坡的高度分布图。
图5f为第70行处的斜坡的高度分布图。
图5g为使用本发明方法获取的第70行处的斜坡高度分布图。
具体实施方式
下面将结合附图和具体实施例对本发明作进一步的解释说明。
实施例一
本实施例采用了如图1所示的迈克尔逊干涉光路构成的同轴相移的干涉测量光路,所测物体为常用的光学元件凹面镜。测量系统包含光源、压电陶瓷微位移器(PZT)、单色图像传感器(CCD)、计算机(PC)。其中,光源由波长为632.8nm的He-Ne激光器和波长为532nm的固体激光器组成。压电陶瓷微位移器为电子工业集团26研究所生产的型号为2D020的微位移器,由计算机控制形成相移子系统。图像传感器为台湾敏通公司生产的型号为MTV-1082CB黑白工业摄像机,通过图像采集卡将图像数据存入计算机,从而形成图像采集子系统。
在实施本测量方法时,在两个光源和分束镜BS1之间分别插入可变光衰减器ND1和ND2,用来调节各光源的强度满足测量要求,调节分束镜BS1,使两不同波长的光束重叠形成一束共路光。共路光经过分束镜BS2后形成两束:物光和参考光。物光照亮待测物体samples并被待测物体调制后沿原路返回,另一路参考光穿过BS2后被安装在微位移器(PZT)上的反射镜M反射后也沿原路返回,返回的物光和参考光在分束镜BS2处叠加形成干涉条纹,被图像传感器CCD记录下来形成干涉条纹图。
光路调整好以后,按照图2所示的流程,执行以下步骤:
步骤一:仅使用波长为λ1=532nm的光源,在单波长的情况照明的情况下,通过电脑控制微位移器带动反射镜M产生λ1/4的相移,同时用计算机控制图像传感器CCD每当完成λ1/4相移时采集干涉图,共采集四幅相移干涉图I1,I2,I3,I4,使用四步相移算法提取出λ1对应的包裹相位保持光路和待测物体不变,挡住光源λ1,仅使用波长为λ2=632.8nm的光源,通过电脑控制微位移器带动反射镜M产生λ2/4的相移,同时控制图像传感器CCD每当完成λ2/4相移时采集干涉图,同样采集四幅相移干涉图I'1,I'2,I'3,I'4,使用四步相移算法提取出λ2对应的包裹相位将两种波长对应的包裹相位转换成高度包裹φ1(x,y)、φ2(x,y),其中高度包裹图分别如图3a和图3b所示;
步骤二:获取高度包裹差分图z(x,y)=φ1(x,y)-φ2(x,y),转换成灰度图如图3c所示;
步骤三:确定波长差分函数δ(n1,n2)=n2λ2-n1λ1,其中,n1,n2为非负整数;
步骤四:将步骤二获得的高度包裹差分图中的点z(x,y)与波长差分函数δ(n1,n2)逐一比较,找出最相近的波长差分值δ(n1,n2),确定n1(x,y)或n2(x,y);
步骤五:则物体的高度可以表示为h(x,y)=n1(x,y)λ1+φ1(x,y)或者h(x,y)=n2(x,y)λ2+φ2(x,y),重构的物体的三维微形貌如图3d所示。
本实施例选择凹面镜作为待测物体,是因为这种光滑连续物体可以用传统常见的依靠横向相邻像素之间的关系解相位包裹,用这种方法得到的待测物体信息,来衡量本发明方法的正确性。选取λ1=632.8nm波长对应的包裹相位进行解包裹,重构的物体的三维微形貌如图3e所示,从图3d与图3e的对比可以发现,这两种方法重构的三维微形貌是一致的。这两种方法解包裹的原理完全不同,但是最后可以获得相同的技术效果,因此,可以证明本发明方法的正确性。
实施例二
为了进一步证明本方法可以用于孤立非连续(存在突变)物体的相位测量。本实施例采用模拟仿真的方法,将本发明方法运用于测量孤立的突变的方柱,如图4a所示。这些方柱的高度是在0-60000纳米的范围内随机生成的,也就是说存在0-60000纳米的随机突变,它们相互之间不存在任何关联。实施步骤同样如图2所示流程图。这些方柱在532纳米波长的光源下获取的高度包裹灰度图如图4b所示,在632.8纳米波长光源下获取的高度包裹灰度图如图4c所示。将两波长高度包裹图相减,得到高度包裹差分灰度图如图4d所示,通过与波长差分函数匹配,找到对应的跳变系数,重构的方柱如图4e所示,可见本发明方法能够准确重构出孤立非连续(存在突变)物体的相位测量。
实施例三
本发明方法大大扩展了干涉精密测量的量程,根据本发明方法,如果采用的测量光源波长之间存在一个最大公约数k,那么本方法可以测量达到的量程为:(λ2/k+1)λ1或(λ1/k+1)λ2。当测量光源波长为632.8纳米和532纳米时,其最大公约数为56,若使用532纳米波长重构,则量程可以达到60648纳米;若使用632.8纳米波长重构,则量程可以达到60748.8纳米。为了证明本发明方法的量程,本实施例仍然采用模拟的方法,将本发明方法用于测量0-60648纳米线性增高的斜坡,斜坡的两侧同时存在突变,待测斜坡如图5a所示。测量中使用波长为532纳米的光源来重构斜坡。图5b和图5c分辨展示了波长为532和632.8纳米波长下的高度包裹灰度图。图5d展示了两波长高度包裹相减后的差分图的灰度图。重构后的斜坡如图5e所示。图5f与图5g分别表示了待测物体与测量后的物体在第70行处的高度分布,可见,测量结果与待测物体的高度是一致的。
以上的实施例中,很好地验证了本发明方法所披露的不仅可以用于测量光滑缓变的物体,还可以用于测量孤立非连续(存在突变)物体的相位(高度),同时,大大拓展了有效测量突变物体的量程。
Claims (5)
1.一种双波长微纳结构相位测量方法,其特征在于包含以下步骤:
步骤一:使用干涉光路或者干涉仪分别获取两种波长λ1、λ2作为计量波长时的干涉条纹图;使用已知的相位提取算法获取两种波长的包裹相位图并转换成高度包裹图φ1(x,y)、φ2(x,y),其中
步骤二:获取高度包裹差分图z(x,y)=φ1(x,y)-φ2(x,y);
步骤三:确定波长差分函数δ(n1,n2)=n2λ2-n1λ1,其中,n1,n2分别为使用两个波长测量且待测物体高度超出测量波长时的跳变次数,为非负整数;
步骤四:将步骤二获得的高度包裹差分图中的点z(x,y)与波长差分函数δ(n1,n2)逐一比较,找出最相近的波长差分函数值δ(n1,n2),获取n1(x,y)或n2(x,y);
步骤五:则物体的高度表示为h(x,y)=n1(x,y)λ1+φ1(x,y)或者h(x,y)=n2(x,y)λ2+φ2(x,y)。
2.根据权利要求1所述的一种双波长微纳结构相位测量方法,其特征在于步骤一中获取两种波长的干涉图时保持待测物体和干涉光路不变的情况下同时或先后采集干涉条纹图。
3.根据权利要求1或2所述的一种双波长微纳结构相位测量方法,其特征在于所述步骤一中的包裹相位图由四步相移算法获得。
4.根据权利要求3所述的一种双波长微纳结构相位测量方法,其特征在于所述两种波长是不成倍频关系的两种波长。
5.根据权利要求4所述的一种双波长微纳结构相位测量方法,其特征在于所述干涉条纹图使用同轴相移迈克尔逊干涉光路获取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410139316.9A CN104034277B (zh) | 2014-04-08 | 2014-04-08 | 一种双波长微纳结构相位测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410139316.9A CN104034277B (zh) | 2014-04-08 | 2014-04-08 | 一种双波长微纳结构相位测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104034277A CN104034277A (zh) | 2014-09-10 |
CN104034277B true CN104034277B (zh) | 2016-08-17 |
Family
ID=51465144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410139316.9A Active CN104034277B (zh) | 2014-04-08 | 2014-04-08 | 一种双波长微纳结构相位测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104034277B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105300276B (zh) * | 2015-11-14 | 2017-11-10 | 华南师范大学 | 一种双波长单曝光干涉测量方法及系统 |
CN106643474B (zh) * | 2016-09-23 | 2019-03-08 | 华南师范大学 | 一种归一化正交多步相移相位测量方法 |
CN106568392B (zh) * | 2016-10-19 | 2019-10-01 | 江苏大学 | 一种双波长相位显微成像系统和方法、以及对应相位恢复方法 |
CN106767523B (zh) * | 2016-11-17 | 2019-03-08 | 南方科技大学 | 一种提高相位精度的方法及装置 |
CN108132026B (zh) * | 2018-01-24 | 2024-02-27 | 赵智亮 | 半导体中红外可见光双波长透射式干涉测试装置 |
CN108917641B (zh) * | 2018-05-15 | 2020-08-04 | 广东工业大学 | 基于激光器波数合成的样件内部轮廓检测方法及系统 |
CN111159251A (zh) * | 2019-12-19 | 2020-05-15 | 青岛聚好联科技有限公司 | 一种异常数据的确定方法及装置 |
CN113155058A (zh) * | 2021-04-29 | 2021-07-23 | 范金坪 | 一种基于空域载频相移的多波长动态相位测量装置及方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001343221A (ja) * | 2000-06-02 | 2001-12-14 | Fujitsu Ltd | 表面計測装置および表面計測方法 |
CN101375332A (zh) * | 2005-09-26 | 2009-02-25 | 皇家飞利浦电子股份有限公司 | 光学补偿器、光学元件、光学扫描头和光学扫描设备 |
CN103630086B (zh) * | 2013-11-13 | 2017-03-15 | 华南师范大学 | 一种基于单色ccd的双波长同时相移干涉测量方法 |
-
2014
- 2014-04-08 CN CN201410139316.9A patent/CN104034277B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104034277A (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104034277B (zh) | 一种双波长微纳结构相位测量方法 | |
CN105180892B (zh) | 一种飞秒激光频率梳脉冲啁啾干涉测距方法及测距系统 | |
Su et al. | Fourier transform profilometry:: a review | |
EP2606311B1 (en) | Apparatus and method for measuring distance | |
CN103630086B (zh) | 一种基于单色ccd的双波长同时相移干涉测量方法 | |
CN102788558B (zh) | 散斑相关和散斑干涉相结合的三维变形测量系统及方法 | |
CN102221342B (zh) | 一种时域多波长外差散斑干涉测量物体变形的方法 | |
CN101629813B (zh) | 基于计算全息的自由曲面三维形貌的测量方法 | |
Kumar et al. | Two-wavelength micro-interferometry for 3-D surface profiling | |
CN106568392A (zh) | 一种双波长相位显微成像系统和方法、以及对应相位恢复方法 | |
CN102620665A (zh) | 使用白光扫描干涉测量法进行膜厚度测量的方法和装置 | |
Du et al. | Iterative two-step temporal phase-unwrapping applied to high sensitivity three-dimensional profilometry | |
CN102538714A (zh) | 一种平面高精度平行度的检测装置 | |
Chen et al. | Phase-unwrapping algorithm for the measurement of 3D object | |
CN106247950A (zh) | 基于广义相移数字全息的微小位移测量方法 | |
Kimbrough et al. | The spatial frequency response and resolution limitations of pixelated mask spatial carrier based phase shifting interferometry | |
Reichold et al. | 2500-Channel single-shot areal profilometer using hyperspectral interferometry with a pinhole array | |
CN103267485A (zh) | 一种点衍射三维绝对位移测量方法 | |
Davila et al. | Measurement of sub-surface delaminations in carbon fibre composites using high-speed phase-shifted speckle interferometry and temporal phase unwrapping | |
Genovese et al. | A comprehensive ESPI based system for combined measurement of shape and deformation of electronic components | |
CN106482633B (zh) | 一种基于π/4相移的多光束干涉相位提取方法 | |
Puga et al. | General model to predict and correct errors in phase map interpretation and measurement for out-of-plane ESPI interferometers | |
CN104330034A (zh) | 台阶表面高度测量的单光栅干涉仪及台阶表面高度测量方法 | |
Cordero et al. | Uncertainty evaluation of displacements measured by electronic speckle-pattern interferometry | |
Cordero et al. | Strain maps obtained by phase-shifting interferometry: An uncertainty analysis |
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 |