CN114088010B - 一种干涉三维形貌解算方法 - Google Patents

一种干涉三维形貌解算方法 Download PDF

Info

Publication number
CN114088010B
CN114088010B CN202111370114.1A CN202111370114A CN114088010B CN 114088010 B CN114088010 B CN 114088010B CN 202111370114 A CN202111370114 A CN 202111370114A CN 114088010 B CN114088010 B CN 114088010B
Authority
CN
China
Prior art keywords
zero
point
envelope curve
extreme point
order
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
Application number
CN202111370114.1A
Other languages
English (en)
Other versions
CN114088010A (zh
Inventor
姚东
梁瀚钢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN202111370114.1A priority Critical patent/CN114088010B/zh
Publication of CN114088010A publication Critical patent/CN114088010A/zh
Application granted granted Critical
Publication of CN114088010B publication Critical patent/CN114088010B/zh
Priority to PCT/CN2022/132783 priority patent/WO2023088409A1/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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)
  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供了一种干涉三维形貌解算方法，包括：S1、对干涉仪进行标定，获得干涉仪的顶峰空间偏移量Δl；S2、利用干涉仪获取待测目标上每个像素采集点P(xi，yj)的时序干涉图阵列；时序干涉图阵列为N张具有干涉条纹的干涉图像；S3、获取时序干涉图阵列的零级包络曲线以及潜在零级包络曲线；S4、利用零级包络曲线得到零级包络曲线的零级极值点Z0(xi，yj)；S5、根据零级极值点Z0(xi，yj)求解待测目标的高程位置Zd(xi，yj)。本发明在空域进行，避免了傅里叶变换、卷积等频域大计算量过程，提高了算法的实时性；通过对子包络曲线拟合进行具有较高的拟合精度，可用于低反射率物体表面的三维形貌解算。

Description

一种干涉三维形貌解算方法
技术领域
本发明属于表面貌形测量技术领域，具体涉及一种干涉三维形貌解算方法。
背景技术
现有技术中已有众多学者围绕白光干涉图解算问题开展了一系列研究工作，将这些研究工作可以总结为三个主要技术发展方向：
1、直接求解法，具体包括插值法、移相法、空间频域法；
2、加权平均法，具体包括重心法、相干相关法；
3、包络曲线拟合法，具体包括多项式拟合法、高斯拟合法、傅里叶变换法、希尔伯特变换法、小波变换法、基于采样定理的包络曲线函数计算法。
在面对完美的信号时，上述每一种方法均会得出近于理想的解，但当光谱退化回射信号或者强度衰减回射信号与标准反射镜的回射信号生成干涉图时，将会发生干涉信号包络曲线展宽、信号对比度下降等问题，同时夹杂多级干涉杂光信号，将产生更为劣化的待处理干涉图(如图2a、图2b所示)。利用上述方法均难以获取高还原精度的目标表面三维形貌。
发明内容
本发明为了解决现有技术中的缺陷，提出了一种干涉三维形貌解算方法，由于多级叠加光场是仪器内部产生的，具有确定性；相同表面零级干涉条纹的顶峰具有唯一性；以及“光谱退化效应”不影响包络曲线峰值相对位置。因此，尽管在“包络曲线展宽”、“信号对比度下降”等问题的影响下，被测目标的干涉强度空间分布依然具有稳定性，通过获取零级干涉条纹位置Z0，将Z0与标定状态下顶峰空间偏移量Δl做和便可求得待测目标的真实位置Zd。
为实现上述目的，本发明采用以下具体技术方案：
一种干涉三维形貌解算方法，包括步骤：
S1、利用标准镜对干涉仪进行标定，获得干涉仪的顶峰空间偏移量Δl；
S2、利用干涉仪获取待测目标上每个像素采集点P(xi,yj)的时序干涉图阵列；时序干涉图阵列为N张具有干涉条纹的干涉图像；
S3、获取时序干涉图阵列的零级包络曲线以及潜在零级包络曲线；
S4、利用零级包络曲线得到零级包络曲线的零级极值点Z0(xi,yj)；
S5、根据零级极值点Z0(xi,yj)求解待测目标的高程位置Zd(xi,yj)：
Zd(xi,yj)＝Δl(xi,yj)+Z0(xi,yj) (1)。
优选地，标准镜为理想标准平面反射镜、球面镜或非球面镜；
标准镜与待测目标具有相同的光谱反射特性。
优选地，步骤S3包括以下步骤：
S31、判断N个位于时序干涉图阵列中同一像素位置(xi,yj)处的采样点P(n,xi,yj)是否为极值点Pe(n,xi,yj)，或为过渡点Pt(n,xi,yj)， n∈[1,M]，M为预设值；
S32、根据极值点Pe(n,xi,yj)得到时序干涉图阵列的零级包络曲线以及潜在零级包络曲线。
优选地，步骤S31包括以下步骤：
S311、对N张干涉图像上的各个采样点P(n,xi,yj)进行曝光，获得各个采样点P(n,xi,yj)在干涉图像上的强度值I(n,xi,yj)；
S312、根据位于时序干涉图阵列中同一位置的N个采样点 P(n,xi,yj)所在的干涉图像上的强度值I(n,xi,yj)，构建数据立方体S；数据立方体S的构建模型如下：
S(n,i,j)＝[I(n+1,xi,yj)-I(n,xi,yj)]×[I(n,xi,yj)-I(n-1,xi,yj)] (2)；
其中，
I(n,xi,yj)表示采样点P(n,xi,yj)在时序干涉图阵列中第n张干涉图像上的强度值；
I(n+1,xi,yj)表示采样点P(n+1,xi,yj)在时序干涉图阵列中第n+1 张干涉图像上的强度值；
I(n-1,xi,yj)表示采样点P(n-1,xi,yj)在时序干涉图阵列中第n-1 张干涉图像上的强度值；
S313、判断采样点P(n,xi,yj)是否为时序干涉图阵列的极值点：
若S(n,i,j)≤0，则采样点P(n,xi,yj)为极值点Pe(n,xi,yj)；
若S(n,i,j)>0，则采样点P(n,xi,yj)为过渡点Pt(n,xi,yj)。
优选地，步骤S32的方法为：
将每个像素采集点P(xi,yj)对应在时序干涉图阵列上的极值点排序，选取至少三个强度值按降序排列的第一极值点Pe(na,xi,yj)、第二极值点Pe(nb,xi,yj)和第三极值点Pe(nc,xi,yj)，存在如下关系：
Ie(na,xi,yj)＜Ie(nb,xi,yj)；
Ie(nb,xi,yj)＞Ie(nc,xi,yj)；
其中，Ie(na,xi,yj)为第一极值点Pe(na,xi,yj)的强度值，
Ie(nb,xi,yj)为第二极值点Pe(nb,xi,yj)的强度值，
Ie(nc,xi,yj)为第三极值点Pe(nc,xi,yj)的强度值；
则第二极值点Pe(nb,xi,yj)位于零级包络曲线；
第一极值点Pe(na,xi,yj)和第三极值点Pe(nc,xi,yj)位于潜在零级包络曲线。
优选地，步骤S4包括以下步骤：
S41、根据位于零级包络曲线的极值点Pe(nk,xi,yj)及过渡点 Pt(nk,xi,yj)获取零级包络曲线的零级极值点Z0(xi,yj)。
优选地，零级极值点Z0(xi,yj)的获取方法为：
选取至少两个位于零级包络曲线的极值点Pe(nk,xi,yj)附近的第一过渡点Pt(nk+1,xi,yj)和第二过渡点Pt(nk-1,xi,yj)，利用拟合的方法，对零级包络曲线的极值点Pe(nk,xi,yj)、第一过渡点Pt(nk+1,xi,yj)和第二过渡点Pt(nk-1,xi,yj)进行拟合，得到零级极值点Z0(xi,yj)。
优选地，拟合的方法包括拟合法、重心法或插值法中的至少一种。
优选地，步骤S4还包括以下步骤：
S42、拟合潜在零级包络曲线的零级极值点Z′0(xi,yj)，利用先验条件对错误的零级极值点Z0(xi,yj)进行替换；
先验条件为：待测目标的连续平面测量数据无高程跳变错位。
优选地，步骤S1包括以下步骤：
S11、利用已知面型的标准镜对干涉仪进行标定，得到的标准镜的标准面型分布
Figure BDA0003362033780000041
S12、通过利用干涉仪获取标准镜上每个像素采集点的时序干涉图阵列，计算获得标准镜的变形分布Z(xi,yj)；
S13、计算干涉仪内部由于光程不匹配所引起的顶峰空间偏移量Δl(xi,yj)：
Figure BDA0003362033780000042
优选地，干涉仪为宽谱段白光干涉仪或窄谱段激光干涉仪。
本发明能够取得以下技术效果：
1、本发明在空间域中进行，避免了傅里叶变换、卷积等频域大计算量过程，提高了算法的实时性。
2、本发明通过对子包络曲线拟合，具有较高的拟合精度，可以应用于低反射率表面的三维形貌解算。
3、本发明适用于低反射率样品，在工业测量、科学研究等实际工况下具有前景。
附图说明
图1是本发明一个实施例的一种干涉三维形貌解算方法的流程图；
图2a是背景技术中利用迈克尔逊白光干涉系统测量获得的干涉信号示意图
图2b是背景技术中的相同系统由于光谱退化效应而导致干涉图发生改变的干涉信号示意图；
图3是本发明一个实施例的零级干涉条纹极大值点Z0、顶峰位置 Zd和标定状态下顶峰空间位置偏移量Δl的相对位置关系图；
图4a是本发明一个实施例的台阶处三维形貌解算图；
图4b是图4a的重心截面图；
图5a是利用傅里叶变换法获得的台阶处三维形貌解算图；
图5b是图5a的重心截面图；
图6a是利用重心法获得的台阶处三维形貌解算图；
图6b是图6a的重心截面图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及具体实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，而不构成对本发明的限制。
本发明的目的是提供一种干涉三维形貌解算方法，下面将对本发明提供的一种干涉三维形貌解算方法，通过具体实施例来进行详细说明。
在常规白光干涉图中，由于“多重干涉叠加效应”、“光谱退化效应”等因素影响，外包络曲线的顶峰值Zd和零级干涉条纹(子包络曲线)的中心位置Z0往往不重合，存在如图3所示的位置关系。因此，本发明提出的干涉三维形貌解算方法通过精确求解出零级干涉条纹中心所在位置Z0和标定状态下顶峰位置偏移Δl，便可获得待测目标的真实位置Zd。
参照图1示出的干涉三维形貌解算方法的流程，该方法包括以下步骤：
S1、利用标准镜对干涉仪进行标定，获得干涉仪的顶峰空间偏移量Δl，存在如下表达式：
Figure BDA0003362033780000051
其中，
Figure BDA0003362033780000052
为经干涉仪得到的标准镜的标准面型分布；
Z(xi,yj)为利用本发明的一种干涉三维形貌解算方法得到的标准镜的变形分布。
S2、利用干涉仪获取待测目标上每个像素采集点P(xi,yj)的时序干涉图阵列。
在本发明的一个优选实施例中，以待测目标为台阶为例，利用完成标定的白光干涉仪获取台阶的采集域内每一个像素采集点P(xi,yj) 的时序干涉图阵列。时序干涉图阵列为N张具有干涉条纹的干涉图像，干涉图像上的每一个像素点具有“空间-强度”的变化信息。
S3、获取时序干涉图阵列的零级包络曲线以及潜在零级包络曲线，具体包括以下步骤：
S31、判断N个位于时序干涉图阵列中同一像素位置(xi,yj)处的采样点P(n,xi,yj)是否为极值点Pe(n,xi,yj)，或为过渡点Pt(n,xi,yj)， n∈[1,M]，M为预设值；
在台阶的采集域内，每一个像素采集点P(xi,yj)都对应一个获取的台阶时序干涉图阵列中的采样点P(n,xi,yj)，判断该采样点P(n,xi,yj) 是极值点Pe(n,xi,yj)还是过渡点Pt(n,xi,yj)。
步骤S31具体包括以下步骤：
S311、对N张干涉图像上的各个采样点P(n,xi,yj)进行曝光，获得每一个采样点P(n,xi,yj)在干涉图像上的强度值I(n,xi,yj)；
其中，i,j代表位置，n代表采样时序，
n∈[1,M]，M为预设值，M的取值根据待测目标的高程值与采样间隔设定。
S312、根据位于时序干涉图阵列中同一位置的N个采样点 P(n,xi,yj)所在的干涉图像上的强度值I(n,xi,yj)，构建数据立方体S。
以任意一个采集点(x1,y1)为例，因时序干涉图阵列由N个干涉图像组成，因此，采集点(x1,y1)对应在时序干涉图阵列中各个干涉图像的采样点为P(1,x1,y1)、P(2,x1,y1)、…、P(n,x1,y1)、P(n+1,x1,y1)、…、 P(N,x1,y1)。对应得到的强度值为I(1,x1,y1)、I(2,x1,y1)、…、I(N,x1,y1)。
同理，获得同一位置(xi,yj)，不同时序的所有采样点P(n,xi,yj)在干涉图像上的强度值I(n,xi,yj)，建立如下的数据立方体S模型：
S(n,i,j)＝[I(n+1,xi,yj)-I(n,xi,yj)]×[I(n,xi,yj)-I(n-1,xi,yj)] (2)；
其中，I(n,xi,yj)表示采样点P(n,xi,yj)在时序干涉图阵列中第n张干涉图像上的强度值；
I(n+1,xi,yj)表示采样点P(n+1,xi,yj)在时序干涉图阵列中第n+1 张干涉图像上的强度值；
I(n-1,xi,yj)表示采样点P(n-1,xi,yj)在时序干涉图阵列中第n-1 张干涉图像上的强度值。
S313、判断采样点P(n,xi,yj)为时序干涉图阵列的极值点，否则为过渡点：
若S(n,i,j)≤0，则采样点P(n,xi,yj)为极值点Pe(n,xi,yj)；
若S(n,i,j)>0，则采样点P(n,xi,yj)为过渡点Pt(n,xi,yj)。
S32、根据极值点Pe(n,xi,yj)得到时序干涉图阵列的零级包络曲线以及潜在零级包络曲线。
在本发明的一个优选实施例中，将得到的每个像素采集点P(xi,yj) 对应在时序干涉图阵列上的极值点排序，选取强度值按降序排列的第一极值点Pe(na,xi,yj)、第二极值点Pe(nb,xi,yj)和第三极值点 Pe(nc,xi,yj)。
即存在第一极值点Pe(na,xi,yj)的强度值Ie(na,xi,yj)＜第二极值点Pe(nb,xi,yj)的强度值Ie(nb,xi,yj)；同时第二极值点Pe(nb,xi,yj)的强度值Ie(nb,xi,yj)＞第三极值点Pe(nc,xi,yj)的强度值Ie(nc,xi,yj)。
则第二极值点Pe(nb,xi,yj)所在的子包络曲线为零级包络曲线；
第一极值点Pe(na,xi,yj)和第三极值点Pe(nc,xi,yj)所在的子包络曲线为潜在零级包络曲线。
对于选取任意前奇数个强度值最大的极值点，中间强度值对应的极值点所在的子包络曲线为零级包络曲线，其他极值点位于潜在零级包络曲线。
在本发明的另一个实施例中，还可以选取两个强度值最大的第四极值点P(nf,xi,yj)和第五极值点P(ng,xi,yj)，对应如下的强度关系： Ie(nf,xi,yj)＜Ie(ng,xi,yj)；
则强度值较大的极值点P(ng,xi,yj)所在的子包络曲线为零级包络曲线；
极值点P(nf,xi,yj)所在的子包络曲线为潜在零级包络曲线。
对于选取任意前偶数个强度值最大的极值点，强度值最大的极值点所在的子包络曲线为零级包络曲线，其他极值点位于潜在零级包络曲线。
对于如图2b所示的多层膜系的情况，可根据实际情况选择各个外包络组中强度极大值的极值点所在的子包络曲线为零级包络曲线。
S4、利用零级包络曲线得到零级包络曲线的零级极值点Z0(xi,yj)。
步骤S4包括以下步骤：
S41，根据位于零级包络曲线的极值点Pe(nk,xi,yj)及过渡点 Pt(nk,xi,yj)通过拟合的方法获取零级包络曲线的零级极值点Z0(xi,yj)。
在本发明的一个优选实施例中，选择第二极值点Pe(nb,xi,yj)及其附近的四个过渡点Pe(nb-2,xi,yj)、Pe(nb-1,xi,yj)、Pe(nb+1,xi,yj)和 Pe(nb+2,xi,yj)共同拟合出零级包络曲线的零级极值点Z0(xi,yj)。
进一步的，选则相对拟合精度较高的拟合法对第二极值点 Pe(nb,xi,yj)、第一过渡点Pe(nb-2,xi,yj)、第二过渡点Pe(nb-1,xi,yj)、第三过渡点Pe(nb+1,xi,yj)和第四过渡点Pe(nb+2,xi,yj)进行拟合，通过计算获得连续的数学曲线，获得零级极值点Z0(xi,yj)。拟合法为现有技术，不在此赘述。
在本发明的另一个实施例中，还可以选择包括但不限于重心法或插值法的方法对上述极值点和过渡点进行拟合，得到零级极值点 Z0(xi,yj)。
S5、根据零级极值点Z0(xi,yj)求解出待测目标的高程位置Zd(xi,yj)：
Zd(xi,yj)＝Δl(xi,yj)+Z0(xi,yj) (1)。
在实际零级包络曲线的极大值点Z0(xi,yj)求解过程中，可能发生解算错位现象，即误把潜在零级包络曲线代替零级包络曲线重新生成测量数据，因此步骤S4还包括：
S42、拟合潜在零级包络曲线的零级极值点Z′0(xi,yj)，利用先验条件对错误的零级极值点Z0(xi,yj)进行替换；
先验条件可根据实际情况设定，可以设定为：待测目标的连续平面测量数据无高程跳变错位。
具体的，按照步骤S5中的方法，对位于潜在零级包络曲线的第一极值点Pe(na,xi,yj)和第三极值点Pe(nc,xi,yj)分别进行拟合，获取潜在零级包络曲线的零级极值点Z′0(xi,yj)和Z″0(xi,yj)。
即对第一极值点Pe(na,xi,yj)以及第一极值点附近的四个过渡点 Pe(na-2,xi,yj)、Pe(na-1,xi,yj)、Pe(na+1,xi,yj)和Pe(na+2,xi,yj)进行拟合，得到潜在零级包络曲线的零级极值点Z′0(xi,yj)；
对第三极值点Pe(nc,xi,yj)以及第三极值点附近的四个过渡点 Pe(nc-2,xi,yj)、Pe(nc-1,xi,yj)、Pe(nc+1,xi,yj)和Pe(nc+2,xi,yj)进行拟合，得到潜在零级包络曲线的零级极值点Z″0(xi,yj)。
假设已知待测台阶的表面为平滑表面，而得到的解算后的三维貌形数据发生了断崖式跳变，则考虑用潜在零级包络曲线的零级极值点 Z′0(xi,yj)或Z″0(xi,yj)替代错误的零级极值点Z0(xi,yj)，利用式(1)重新求解出待测目标的高程位置Zd(xi,yj)。
因此，本发明提出的干涉三维形貌解算方法，一方面没有使用高斯函数进行滤波，降低了计算量；
一方面不需要对全部的采样点进行滤波，提高了算法的实时性；同时利用了子包络曲线具有更大的陡度的特点，使得算法的拟合精度高。
另一方面，本发明同样适用于窄谱段激光光源的干涉仪，如法布里-珀罗干涉仪、迈克尔逊干涉仪、马赫曾德干涉仪、塞格纳克干涉仪、斐索干涉仪等。
在本发明的一个优选实施例中，步骤S1中对干涉仪进行标定，获得干涉仪的顶峰空间偏移量Δl，具体包括以下步骤：
S11、利用已知面型的标准镜对干涉仪进行标定，计算经干涉仪获得的标准镜的标准面型分布
Figure BDA0003362033780000091
在本发明的另一个优选实施例中，标准镜选取理想标准平面反射镜，利用已知面型的具有理想平面的标准平面反射镜对宽谱段白光干涉仪进行标定。
由于标准镜的反射信号强，对比度高，因此可利用常规方法，如小波变换法计算获取标准平面反射镜的面型分布
Figure BDA0003362033780000101
将其视为标准值。
在本发明的另一个实施例中，标准镜还可选取已知面型的球面镜或非球面镜；通过选择与待测目标具有相同的光谱反射特性的标准镜，可以使用于标定的干涉图与实际测量的干涉图具有更好的相似度，有利于提高测量准确度。
S12、通过利用干涉仪获取标准镜上每个像素采集点的时序干涉图阵列，计算获得标准镜的变形分布Z(xi,yj)。
具体包括以下步骤：
S121、获取标准镜的时序干涉图阵列的零级包络曲线；
S122、利用标准镜的时序干涉图阵列的零级包络曲线得到标准镜的零级包络曲线的极值点和过渡点；
S123、利用标准镜的零级包络曲线的极值点和过渡点获取标准镜上每个像素采集点的零级包络曲线的极点高程位置；
S124、根据高程位置计算获得标准镜的变形分布Z(xi,yj)。
步骤S12中标准镜的零级包络曲线的极值点和过渡点的获取方法与上述实施例中台阶的零级包络曲线的极值点和过渡点的获取方法相同，在此不再赘述。
S13、计算干涉仪内部由于光程不匹配所引起的顶峰空间偏移量Δl(xi,yj)：
Figure BDA0003362033780000102
图4a、图5a和图6a分别示出了利用本发明的方法、傅里叶变换法以及重心法获取同一台阶处三维形貌解算图,图4b、图5b和图6b 分别示出了对应的重心截面图。
从图中可以看出，利用本发明计算获取台阶上表面面型平整度RMS值为0.038μm，下表面面型平整度RMS值为0.009μm，高度为 1.755μm，与标准值1.761μm的差值为6nm。
相比于传统的傅里叶变换法(上表面RMS值为0.045μm，下表面RMS值为0.015μm)，上表面RMS值减小了15.6％，下表面RMS 值减小了40％；
相比于传统的重心法(上表面RMS值为0.045μm，下表面RMS 值为0.017μm)，上表面RMS值减小了15.6％，下表面RMS值减小了47％。
通过上述数据可以发现，本发明的方法能够表现出更为准确的平面特征，测量获取的高程测量误差位于±0.01μm区间(标准台阶给定误差带范围)，具有非常高的可信度。
另一方面，台阶样品以硅为基底，经过化学刻蚀而制成，硅材质偏黑色，具有相对较低的光谱反射率，导致信号的信噪比相对较低，因此解算效果并不理想。但利用本发明提出的方法，即可以大幅提高表面面型精度(RMS值)，又使得高度差在标定范围内。可应用于低反射率样品，具有应用于工业测量、科学研究等实际工况的前景。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
以上本发明的具体实施方式，并不构成对本发明保护范围的限定。任何根据本发明的技术构思所作出的各种其他相应的改变与变形，均应包含在本发明权利要求的保护范围内。

Claims (11)

1.一种干涉三维形貌解算方法，其特征在于，包括步骤：
S1、利用标准镜对干涉仪进行标定，获得所述干涉仪的顶峰空间偏移量Δl；
S2、利用所述干涉仪获取待测目标上每个像素采集点P(xi，yj)的时序干涉图阵列；所述时序干涉图阵列为N张具有干涉条纹的干涉图像；
S3、获取所述时序干涉图阵列的零级子包络曲线以及潜在零级子包络曲线；
S4、利用所述零级子包络曲线得到所述零级子包络曲线的零级极值点Z0(xi，yj)；
S5、根据所述零级极值点Z0(xi，yj)求解待测目标的高程位置Zd(xi，yj)：
Zd(xi，yj)＝Δl(xi，yj)+Z0(xi，yj) (1)。
2.根据权利要求1所述的干涉三维形貌解算方法，其特征在于，所述标准镜为理想标准平面反射镜、球面镜或非球面镜；
所述标准镜与待测目标具有相同的光谱反射特性。
3.根据权利要求1所述的干涉三维形貌解算方法，其特征在于，所述步骤S3包括以下步骤：
S31、判断N个位于所述时序干涉图阵列中同一像素位置(xi，yj)处的采样点P(n，xi，yj)是否为极值点Pe(n，xi，yj)，或为过渡点Pt(n，xi，yj)，n∈[1，M]，M为预设值；
S32、根据所述极值点Pe(n，xi，yj)得到所述时序干涉图阵列的零级子包络曲线以及潜在零级子包络曲线。
4.根据权利要求3所述的干涉三维形貌解算方法，其特征在于，所述步骤S31包括以下步骤：
S311、对N张干涉图像上的各个采样点P(n，xi，yj)进行曝光，获得各个采样点P(n，xi，yj)在所述干涉图像上的强度值I(n，xi，yj)；
S312、根据位于所述时序干涉图阵列中同一位置的N个所述采样点P(n，xi，yj)所在的所述干涉图像上的强度值I(n，xi，yj)，构建数据立方体S；
所述数据立方体S的构建模型如下：
S(n，i，j)＝[I(n+1，xi，yj)-I(n，xi，yj)]×[I(n，xi，yj)-I(n-1，xi，yj)] (2)；
其中，
I(n，xi，yj)表示采样点P(n，xi，yj)在所述时序干涉图阵列中第n张干涉图像上的强度值；
I(n+1，xi，yj)表示采样点P(n+1，xi，yj)在所述时序干涉图阵列中第n+1张干涉图像上的强度值；
I(n-1，xi，yj)表示采样点P(n-1，xi，yj)在所述时序干涉图阵列中第n-1张干涉图像上的强度值；
S313、判断所述采样点P(n，xi，yj)是否为所述时序干涉图阵列的极值点：
若S(n，i，j)≤0，则所述采样点P(n，xi，yj)为所述极值点Pe(n，xi，yj)；
若S(n，i，j)＞0，则所述采样点P(n，xi，yj)为过渡点Pt(n，xi，yj)。
5.根据权利要求3所述的干涉三维形貌解算方法，其特征在于，所述步骤S32的方法为：
将每个像素采集点P(xi，yj)对应在所述时序干涉图阵列上的所述极值点排序，选取至少三个强度值按降序排列的第一极值点Pe(na，xi，yj)、第二极值点Pe(nb，xi，yj)和第三极值点Pe(nc，xi，yj)，存在如下关系：
Ie(na，xi，yj)＜Ie(nb，xi，yj)；
Ie(nb，xi，yj)＞Ie(nc，xi，yj)；
其中，Ie(na，xi，yj)为所述第一极值点Pe(na，xi，yj)的强度值，
Ie(nb，xi，yj)为所述第二极值点Pe(nb，xi，yj)的强度值，
Ie(nc，xi，yj)为所述第三极值点Pe(nc，xi，yj)的强度值；
则所述第二极值点Pe(nb，xi，yj)位于所述零级子包络曲线；
所述第一极值点Pe(na，xi，yj)和所述第三极值点Pe(nc，xi，yj)位于所述潜在零级子包络曲线。
6.根据权利要求1所述的干涉三维形貌解算方法，其特征在于，所述步骤S4包括以下步骤：
S41、根据位于零级子包络曲线的极值点Pe(nk，xi，yj)及过渡点Pt(nk，xi，yj)获取零级子包络曲线的零级极值点Z0(xi，yj)。
7.根据权利要求6所述的干涉三维形貌解算方法，其特征在于，所述零级极值点Z0(xi，yj)的获取方法为：
选取至少两个位于所述零级子包络曲线的极值点Pe(nk，xi，yj)附近的第一过渡点Pt(nk+1，xi，yj)和第二过渡点Pt(nk-1，xi，yj)，利用拟合的方法，对所述零级子包络曲线的极值点Pe(nk，xi，yj)、所述第一过渡点Pt(nk+1，xi，yj)和所述第二过渡点Pt(nk-1，xi，yj)进行拟合，得到所述零级极值点Z0(xi，yj)。
8.根据权利要求7所述的干涉三维形貌解算方法，其特征在于，所述拟合的方法包括拟合法、重心法或插值法中的至少一种。
9.根据权利要求1所述的干涉三维形貌解算方法，其特征在于，所述步骤S4还包括以下步骤：
S42、拟合所述潜在零级子包络曲线的零级极值点Z′0(xi，yj)，利用先验条件对错误的零级极值点Z0(xi，yj)进行替换；
所述先验条件为：待测目标的连续平面测量数据无高程跳变错位。
10.根据权利要求1所述的干涉三维形貌解算方法，其特征在于，所述步骤S1包括以下步骤：
S11、利用已知面型的标准镜对所述干涉仪进行标定，得到的所述标准镜的标准面型分布
Figure FDA0003742688550000031
S12、通过利用所述干涉仪获取所述标准镜上每个像素采集点的时序干涉图阵列，计算获得所述标准镜的变形分布Z(xi，yj)；
S13、计算所述干涉仪内部由于光程不匹配所引起的顶峰空间偏移量Δl(xi，yj)：
Figure FDA0003742688550000032
11.根据权利要求1或10所述的干涉三维形貌解算方法，其特征在于，所述干涉仪为宽谱段白光干涉仪或窄谱段激光干涉仪。
CN202111370114.1A 2021-11-18 2021-11-18 一种干涉三维形貌解算方法 Active CN114088010B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202111370114.1A CN114088010B (zh) 2021-11-18 2021-11-18 一种干涉三维形貌解算方法
PCT/CN2022/132783 WO2023088409A1 (zh) 2021-11-18 2022-11-18 一种干涉三维形貌解算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111370114.1A CN114088010B (zh) 2021-11-18 2021-11-18 一种干涉三维形貌解算方法

Publications (2)

Publication Number Publication Date
CN114088010A CN114088010A (zh) 2022-02-25
CN114088010B true CN114088010B (zh) 2022-08-19

Family

ID=80301745

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111370114.1A Active CN114088010B (zh) 2021-11-18 2021-11-18 一种干涉三维形貌解算方法

Country Status (2)

Country Link
CN (1) CN114088010B (zh)
WO (1) WO2023088409A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114088010B (zh) * 2021-11-18 2022-08-19 中国科学院长春光学精密机械与物理研究所 一种干涉三维形貌解算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102916339A (zh) * 2012-04-19 2013-02-06 广东工业大学 一种基于干涉条纹图像相位反馈的激光控制器装置及其方法
CN109163672A (zh) * 2018-08-03 2019-01-08 广东工业大学 一种基于白光干涉零光程差位置拾取算法的微观形貌测量方法
CN110488292A (zh) * 2019-08-23 2019-11-22 长沙天仪空间科技研究院有限公司 一种基于卫星编队的遥感系统
CN113446930A (zh) * 2020-03-26 2021-09-28 山东大学 一种基于希尔伯特变换的白光干涉信号非均匀采样的校正方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1939581A1 (en) * 2006-12-27 2008-07-02 Heliotis AG Apparatus for the contact-less, interferometric determination of surface height profiles and depth scattering profiles
JP4908273B2 (ja) * 2007-03-15 2012-04-04 富士通株式会社 表面形状計測装置及び表面形状計測方法
CN108759709B (zh) * 2018-03-15 2020-03-27 北京航空航天大学 一种适用于表面形貌检测的白光干涉三维重建方法
CN113091634B (zh) * 2021-03-01 2022-09-13 南京理工大学 一种适用于白光扫描干涉的微观形貌快速测量方法
CN113066083B (zh) * 2021-04-25 2022-08-02 青岛海信医疗设备股份有限公司 流体的多普勒参数确定方法和电子设备
CN114088010B (zh) * 2021-11-18 2022-08-19 中国科学院长春光学精密机械与物理研究所 一种干涉三维形貌解算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102916339A (zh) * 2012-04-19 2013-02-06 广东工业大学 一种基于干涉条纹图像相位反馈的激光控制器装置及其方法
CN109163672A (zh) * 2018-08-03 2019-01-08 广东工业大学 一种基于白光干涉零光程差位置拾取算法的微观形貌测量方法
CN110488292A (zh) * 2019-08-23 2019-11-22 长沙天仪空间科技研究院有限公司 一种基于卫星编队的遥感系统
CN113446930A (zh) * 2020-03-26 2021-09-28 山东大学 一种基于希尔伯特变换的白光干涉信号非均匀采样的校正方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于白光干涉术的微器件三维表面形貌恢复研究;申路等;《佛山科学技术学院学报(自然科学版)》;20120731;第30卷(第4期);第36页倒数第2段至第38页最后一段 *

Also Published As

Publication number Publication date
CN114088010A (zh) 2022-02-25
WO2023088409A1 (zh) 2023-05-25

Similar Documents

Publication Publication Date Title
Hu et al. Absolute three-dimensional micro surface profile measurement based on a Greenough-type stereomicroscope
JPH02500217A (ja) 拡張範囲モアレ等高線作成法
TWI646303B (zh) 決定裝置之儀器轉移函數之方法以及系統
CN108759709A (zh) 一种适用于表面形貌检测的白光干涉三维重建方法
WO2013134966A1 (zh) 用于非球面测量的波长扫描干涉仪及其应用方法
CN109737892B (zh) 基于区域定位拟合算法的数字莫尔移相干涉面形测量方法
CN114088010B (zh) 一种干涉三维形貌解算方法
CN112525070B (zh) 基于非均匀采样校正的抗振动白光干涉测量方法
Lehmann et al. Reduction of chromatic aberration influences in vertical scanning white-light interferometry
CN106123805B (zh) 基于白光扫描干涉的镀膜器件三维形貌测量方法
CN104330027B (zh) 一种基于误差互补修正的移相干涉测量位相提取方法
CN102620680A (zh) 一种三平面绝对测量光学面形的检测装置及方法
CN111121661B (zh) 用于光滑表面形貌测量的窄带非单色光n+1幅相移测试算法
US7158914B2 (en) Precision surface measurement
Dai et al. A novel material measure for characterising two-dimensional instrument transfer functions of areal surface topography measuring instruments
WO2023124867A1 (zh) 基于伪Wigner-Ville分布的白光扫描干涉三维重建方法
Fortmeier et al. Comparison of form measurement results for optical aspheres and freeform surfaces
You et al. Minimization of spectral phase errors in spectrally resolved white light interferometry by the iterative least-squared phase-shifting method
CN103398666B (zh) 一种用于双层周期性微结构的层间错位测试方法
US9702829B1 (en) Systems and methods for wafer surface feature detection and quantification
Lowitzsch et al. Vision and modeling of specular surfaces
US8692999B1 (en) Crosstalk cancellation for a simultaneous phase shifting interferometer
WO2001007864A2 (en) Method and apparatus for performing calibration-based absolute shape measurement using phase-shifting projected fringe profilometry
Qian et al. Calibration of phase shift from two fringe patterns
Chen et al. Dual phase-shifting moire projection with tunable high contrast fringes for three-dimensional microscopic surface profilometry

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