CN103852030A - 用于倾斜波面非零位干涉***的自由曲面面形重构方法 - Google Patents
用于倾斜波面非零位干涉***的自由曲面面形重构方法 Download PDFInfo
- Publication number
- CN103852030A CN103852030A CN201410098162.3A CN201410098162A CN103852030A CN 103852030 A CN103852030 A CN 103852030A CN 201410098162 A CN201410098162 A CN 201410098162A CN 103852030 A CN103852030 A CN 103852030A
- Authority
- CN
- China
- Prior art keywords
- measured
- point
- opl
- centerdot
- corrugated
- 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.)
- Granted
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明公开了一种用于倾斜波面非零位干涉***的自由曲面面形重构方法,该方法基于微扰理论及哈密尔顿特征函数对倾斜波面非零位干涉***测试光路的分析,提出了一种回程误差逆向消除算法,同时将计算生成虚拟波面技术应用于倾斜波面非零位干涉***,解决了通过采集的阵列干涉图如何高精度反演待测件面形的问题。在求解实际待测件面形过程中,提出了一种迭代算法用于面形的最终求解。
Description
技术领域
本发明属于光学精密测试领域,具体涉及一种用于倾斜波面非零位干涉***的自由曲面面形重构方法,可以从倾斜波面非零位干涉***获得的干涉条纹高精度的计算和重构自由曲面的面形。
背景技术
光学自由曲面元件在校正各类像差的同时还能简化光学***的结构,使光学***的性能得到显著提升,因此在照明领域、显示领域、生物医学领域等被大量使用。但其过多的“自由度”,使得高精度加工和检测该类元件变得非常困难,制约了它在成像领域的应用推广。自由曲面光学元件的面形精密检测问题已经成为世界光学领域的一个研究热点。2007年德国斯图加特大学Osten教授团队首先提出倾斜波面非零位干涉法(Tilted-Wave-Interferometer,TWI)针对大梯度非球面元件面形检测提出了一种给测试波前预载倾斜量的非零位干涉检测方法。该方法已经成功应用于表面梯度变化达到±14°的非球面面形测量,测量精度优于λ/30(RMS值)。但是由于倾斜波面干涉法是针对大梯度非球面面形检测提出的,要将其应用于光学自由曲面面形检测,需要解决众多关键技术。
倾斜波面非零位干涉***由于采用了轴外球面波点源产生多束倾斜波面,该干涉***在光路形式上是不满足零位干涉条件的,造成了干涉条纹的相位与待测件面形偏差之间的关系不再那么直接。同时,非零位干涉光路在测量过程中还会给干涉条纹引入较大的回程误差。因此,从干涉条纹中如何高精度的计算和重构自由曲面面形是研究倾斜波面非零位干涉***的核心问题。
目前,国际上只有德国斯图加特大学针对倾斜波面非零位干涉***提出了待测件面形的重构方法。该方法的基本思想就是研究倾斜波面非零位干涉***采集到的干涉条纹的相位与待测件面形表达式的系数之间的数学关系,建立数学模型。每次测量时,首先计算得到干涉图的相位,然后通过建立的数学模型,得到待测件面形表达多项式的实际系数,然后再通过多项式系数设计值与实际值的差值计算待测件的面形偏差。该方法的前提基础是待测件的面形可以用多项式精确表述。由于该方法是针对大梯度非球面测量提出的,而非球面面形都可以用多项式进行精确表述,所以在非球面测量中可以完全适用。但是面对光学自由曲面,特别是那些无法用多项式精确描述的曲面,该方法的使用就受到了限制。
如果能够像传统干涉***一样,从倾斜波面非零位干涉***采集到的干涉条纹的相位直接计算重构出待测件面形,那么就能适合自由曲面的面形测量。
发明内容
本发明的目的在于从倾斜波面非零位干涉***获得的干涉条纹中高精度的计算和重构被测自由曲面面形。
技术方案为:
第一步、根据矢量光学理论建立非旋转对称空间的倾斜波面非零位干涉***的光线追迹模型,生成虚拟波面,获得待测件是理想面形时倾斜波面干涉***中CCD采集到的波面Wnom,使用倾斜波面非零位干涉***测量得到实际待测件时CCD采集到的干涉波面Wrea;
第二步、根据微扰理论及哈密尔顿特征函数,得到理想干涉***中,面形偏差引起的入射至待测面的光线长度变化加上由待测面反射的光线长度变化δOPL(∈)的数学表达式,即:
δOPL(∈)=Wrea-Wrom (1)
第三步、根据回程误差逆向消除算法,消除在实际倾斜波面非零位干涉***中存在的回程误差:
3-1)δOPL(E理论表达式为:
3-2)进行待测件面形重构,相对于CCD的各个像素坐标,其中用像素坐标矩阵表示,的光线在待测件表面的交点B的坐标,用坐标矩阵BC表示,相对于CCD的各个像素坐标的待测件表面B点的入射光线和反射光线的方向余弦,分别用矢量矩阵Inc和Ref表示,即:
式中M表示CCD的横向像素总数,N表示CCD的纵向像素总数,ij表示CCD上像素的位置序号。
3-3)采用迭代算法,求解代表实际待测件面形的矩阵A:
3-3-1)、已知B点坐标矩阵入射至待测件B点的入射光线矢量矩阵从待测件B点反射的反射光线矢量矩阵入射光线与干涉***焦面的交点OC的坐标矩阵光线的光程矩阵待测件为理想面形时的干涉波面实际待测件的干涉波面以及入射光线的光程式中Q表示入射光线的序号,Q的初值为1;
第四步、求得每个点源负责测试的待测件对应区域的面形数据矩阵A,最后将得到的序列矩阵Am×n,其中m×n为点源数量,根据坐标位置进行拼接得到全口径待测件的面形。
在第三步中,回程误差逆向消除算法使理想情况和实际情况下的反射光线到达CCD的同一位置,即式(1)的波面点对点相减具有物理意义。
附图说明
图1传统方法计算分析干涉***中面形偏差引起的δOPL(∈)的过程(1表示理想面形,2表示实际面形,3表示焦平面)
图2计算分析干涉***中面形偏差引起的δOPL(∈)的新方法示意图(未考虑***误差;1表示理想面形,2表示实际面形,3表示焦平面,4表示入射球面波)
图3计算分析干涉***中面形偏差引起的δOPL(∈)的新方法示意图(考虑***误差;1表示理想面形,2表示实际面形,3表示焦平面,5表示入射波面)
图4用于倾斜波面非零位干涉***的自由曲面面形重构算法流程图
具体实施方式
下面结合附图对本发明做进一步详细描述。
根据矢量光学理论,建立非旋转对称空间的倾斜波面非零位干涉***的光线追迹模型,运用计算生成虚拟波面技术生成待测件是理想面形时倾斜波面非零位干涉***中CCD采集到的波面Wnom。使用倾斜波面非零位干涉***测量得到实际待测件时CCD采集到的干涉波面Wrea;根据微扰理论及哈密尔顿特征函数分析倾斜波面非零位干涉***,可得到结论:待测面为理想面形时的干涉***的测试光路的光程与待测面为实际面形时的测试光路的光程相比,它们的差仅仅是由入射至实际待测面的光线长度以及由待测面反射的光线长度与理想情况下的这两种光线长度的差引起的,与光线在干涉***内部的路径是否发生改变无关。所以有式子:
δOPL(∈)=Wrea-Wnom (1)
式中δOPL(∈)表示面形偏差引起的入射至待测面的光线长度变化加上由待测面反射的光线长度变化。
结合图1,描述了传统方法计算分析干涉***中面形偏差引起的δOPL(∈)的过程。当待测件为理想面形时,入射光线入射至待测件表面,交点为A,由待测件反射的光线为当待测件为实际面形时,入射光线入射至待测件表面,交点为B,由待测件反射的光线为由图1的几何关系可知,δOPL(∈)可由下式决定:
图1所展示的传统分析方法对于零位干涉***是非常适合的,并且通过式子(1)以及式子(2)可以精确的计算出待测件的面形偏差。针对非零位干涉***的回程误差,本发明提出了一种回程误差逆向消除算法计算分析倾斜波面非零位干涉***的δOPL(∈)。
结合图2,当待测件是理想面形时,干涉***出射的光线入射至待测件,在B点发生反射,反射光线为反射回干涉***中。同样,当待测件为实际面形时(即存在面形偏差)时,干涉***的出射光线为其入射至待测件在点A发生反射,为反射回干涉***的光线。为了使式子(1)的波面点对点相减具有物理意义,那么理想情况和实际情况下的反射光线必须到达CCD的同一位置,即光线与光线必须为同一条反射光线。根据图2中几何关系,我们可以得到δOPL(∈)的表达式为:
式中点O为干涉***出射波面的会聚点即干涉测试光路的焦点。
将式子(4)代入式子(3),则δOPL(∈)可由下式决定:
结合图3,实际倾斜波面非零位干涉***出射的波阵面具有很大的畸变,导致入射光线和与干涉***测试光路的焦面的交点不会是同一个点,即出射波阵面不会会聚于O点一点。入射光线与干涉***测试光路的焦面相交于点O1,而入射光线与干涉***测试光路的焦面相交于点O2,因此式子(5)就变化为:
下面根据式子(1)和式子(6)来重构出实际待测件的面形。
因为倾斜波面非零位干涉***各器件光学及结构参数已知,可得:光线的光程和光线的光程入射光线与干涉***测试光路的焦面的交点的坐标矩阵:相对于CCD的各个像素坐标(用像素坐标矩阵表示)的光线在待测件表面的交点B的坐标,用坐标矩阵BC;相对于CCD的各个像素坐标的待测件表面B点的入射光线和反射光线的方向余弦,分别用矢量矩阵Inc和Ref表示,即:
式中M表示CCD的横向像素总数,N表示CCD的纵向像素总数,ij表示CCD上像素的位置序号。
在由式子(1)、(6)、(7)、(8)和(9)联立的方程式中,含有两个未知量,即:与矩阵A相关的量和
结合图4,根据上述情况,本发明提出了一种迭代算法来最终求解代表实际待测件面形的矩阵A。算法过程如下:
1、已知B点坐标矩阵入射至待测件B点的入射光线矢量矩阵从待测件B点反射的反射光线矢量矩阵入射光线与干涉***焦面的交点OC的坐标矩阵光线的光程矩阵待测件为理想面形时的干涉波面实际待测件的干涉波面以及入射光线的光程式中Q表示入射光线的序号,Q的初值为1。
求得每个点源负责测试的待测件对应区域的面形数据矩阵A,最后将得到的序列矩阵Am×n(m×n为点源数量)根据坐标位置进行拼接就能得到全口径待测件的面形。
实施例1
下面将通过仿真实验来验证该算法的有效性及精度。
仿真验证实验的基本思想是:利用计算生成虚拟波面技术得到一个已知面形的待测件的理想干涉波面Wnom。然后利用随机函数给这个待测件的面形叠加一个人为的面形误差,并记录该面形误差数据ΔWerror。再利用计算生成虚拟波面技术得到已经叠加了面形误差的待测件干涉波阵面Wrea,将它作为干涉***测试实际待测件得到的波面。接着利用本发明提出的可消除回程误差的面形重构算法,获得待测件的面形偏差ΔW′error。最后将波面ΔWerror和ΔW′error进行点对点相减,得到的残差波面的PV值和RMS值就可以验证本发明提出的光学自由曲面面形重构及评价方法的准确性。
仿真验证实验过程如下:利用一个口径为Φ60mm、最佳拟合球面曲率半径为143mm的凸抛物面作为待测件。
当待测件面形为理想面形时,通过计算生成虚拟波面技术计算得到的干涉***一个点光源的干涉波面。测试待测件Φ40mm口径的区域,该区域内的面形PV值为5.157λ、RMS值为1.452λ。
利用随机函数给该凸抛物面的全口径叠加了一个PV值为1um的面形偏差数据。在该凸抛物面Φ40mm被测区域内人为叠加的面形偏差数据波面,其PV值为0.321λ、RMS值为0.066λ。
利用计算生成虚拟波面技术对叠加过面形偏差数据的凸抛物面再次进行仿真计算,得到干涉波面。该波面的PV值为5.445λ、RMS值为1.458λ。
针对上述的波面Wnom和Wrea,利用本发明提出的可消除回程误差的面形重构算法,计算得到了该待测件的面形偏差数据波面,该波面PV值为0.3461λ、RMS值为0.0655λ。
可知计算得到的面形偏差数据波面与已知的面形偏差数据波面的PV值和RMS值都十分接近。另外,将计算得到的面形偏差数据波面与已知的面形偏差数据波面进行点对点相减,得到残差波面PV值为0.04λ、RMS值为0.003λ,这也验证了本发明提出的光学自由曲面面形重构算法的精度优于λ/20(PV值)。
Claims (2)
1.一种用于倾斜波面非零位干涉***的自由曲面面形重构方法,其特征在于,方法步骤如下:
第一步、根据矢量光学理论建立非旋转对称空间的倾斜波面非零位干涉***的光线追迹模型,生成虚拟波面,获得待测件是理想面形时倾斜波面干涉***中CCD采集到的波面Wnom,使用倾斜波面非零位干涉***测量得到实际待测件时CCD采集到的干涉波面Wrea;
第二步、根据微扰理论及哈密尔顿特征函数,得到理想干涉***中,面形偏差引起的入射至待测面的光线长度变化加上由待测面反射的光线长度变化δOPL(∈)的数学表达式,即:
δOPL(∈)=Wrea-Wrom (1)
第三步、根据回程误差逆向消除算法,消除在实际倾斜波面非零位干涉***中存在的回程误差:
3-1)δOPL(∈)理论表达式为:
3-2)进行待测件面形重构,相对于CCD的各个像素坐标,其中用像素坐标矩阵表示,的光线在待测件表面的交点B的坐标,用坐标矩阵BC表示,相对于CCD的各个像素坐标的待测件表面B点的入射光线和反射光线的方向余弦,分别用矢量矩阵Inc和Ref表示,即:
式中M表示CCD的横向像素总数,N表示CCD的纵向像素总数,ij表示CCD上像素的位置序号。
3-3)采用迭代算法,求解代表实际待测件面形的矩阵A:
3-3-1)、已知B点坐标矩阵入射至待测件B点的入射光线矢量矩阵从待测件B点反射的反射光线矢量矩阵入射光线与干涉***焦面的交点OC的坐标矩阵光线的光程矩阵待测件为理想面形时的干涉波面实际待测件的干涉波面以及入射光线的光程式中Q表示入射光线的序号,Q的初值为1;
第四步、求得每个点源负责测试的待测件对应区域的面形数据矩阵A,最后将得到的序列矩阵Am×n,其中m×n为点源数量,根据坐标位置进行拼接得到全口径待测件的面形。
2.根据权利要求1所述用于倾斜波面非零位干涉***的自由曲面面形重构方法,其特征在于:第三步,回程误差逆向消除算法使理想情况和实际情况下的反射光线到达CCD的同一位置,即式(1)的波面点对点相减具有物理意义。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410098162.3A CN103852030B (zh) | 2014-03-17 | 2014-03-17 | 用于倾斜波面非零位干涉***的自由曲面面形重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410098162.3A CN103852030B (zh) | 2014-03-17 | 2014-03-17 | 用于倾斜波面非零位干涉***的自由曲面面形重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103852030A true CN103852030A (zh) | 2014-06-11 |
CN103852030B CN103852030B (zh) | 2016-04-27 |
Family
ID=50859976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410098162.3A Active CN103852030B (zh) | 2014-03-17 | 2014-03-17 | 用于倾斜波面非零位干涉***的自由曲面面形重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103852030B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104251672A (zh) * | 2014-10-13 | 2014-12-31 | 南京理工大学 | 非零位干涉***中自由曲面待测件空间姿态调整方法 |
CN105783780A (zh) * | 2016-04-29 | 2016-07-20 | 浙江大学 | 一种自由曲面非常规子孔径拼接干涉检测装置与方法 |
CN108362202A (zh) * | 2018-02-02 | 2018-08-03 | 苏州科技大学 | 倾斜波面干涉测量非球面中的参数确定方法 |
CN111308593A (zh) * | 2020-03-30 | 2020-06-19 | 南京理工大学 | 用于手机镜头的超薄型反射式自由曲面设计方法 |
CN113091638A (zh) * | 2021-03-25 | 2021-07-09 | 中国科学院光电技术研究所 | 一种标定干涉法测量面形中的回程误差的装置及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070195084A1 (en) * | 2006-02-03 | 2007-08-23 | Chi-Fai Cheung | Methods and apparatus of aligning surfaces |
CN102607454A (zh) * | 2011-02-24 | 2012-07-25 | 南京理工大学 | 光学自由曲面干涉检测装置 |
CN102997864A (zh) * | 2012-12-17 | 2013-03-27 | 北京理工大学 | 一种大口径光学非球面镜检测*** |
CN103575229A (zh) * | 2012-07-27 | 2014-02-12 | 南京理工大学 | 共光路型多重倾斜波面补偿非零位干涉测量装置 |
-
2014
- 2014-03-17 CN CN201410098162.3A patent/CN103852030B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070195084A1 (en) * | 2006-02-03 | 2007-08-23 | Chi-Fai Cheung | Methods and apparatus of aligning surfaces |
CN102607454A (zh) * | 2011-02-24 | 2012-07-25 | 南京理工大学 | 光学自由曲面干涉检测装置 |
CN103575229A (zh) * | 2012-07-27 | 2014-02-12 | 南京理工大学 | 共光路型多重倾斜波面补偿非零位干涉测量装置 |
CN102997864A (zh) * | 2012-12-17 | 2013-03-27 | 北京理工大学 | 一种大口径光学非球面镜检测*** |
Non-Patent Citations (1)
Title |
---|
沈华等: "基于点源阵列的自由曲面非零位干涉检测***设计方法", 《光学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104251672A (zh) * | 2014-10-13 | 2014-12-31 | 南京理工大学 | 非零位干涉***中自由曲面待测件空间姿态调整方法 |
CN104251672B (zh) * | 2014-10-13 | 2016-09-21 | 南京理工大学 | 非零位干涉***中自由曲面待测件空间姿态调整方法 |
CN105783780A (zh) * | 2016-04-29 | 2016-07-20 | 浙江大学 | 一种自由曲面非常规子孔径拼接干涉检测装置与方法 |
CN108362202A (zh) * | 2018-02-02 | 2018-08-03 | 苏州科技大学 | 倾斜波面干涉测量非球面中的参数确定方法 |
CN108362202B (zh) * | 2018-02-02 | 2020-09-25 | 苏州科技大学 | 倾斜波面干涉测量非球面中的参数确定方法 |
CN111308593A (zh) * | 2020-03-30 | 2020-06-19 | 南京理工大学 | 用于手机镜头的超薄型反射式自由曲面设计方法 |
CN113091638A (zh) * | 2021-03-25 | 2021-07-09 | 中国科学院光电技术研究所 | 一种标定干涉法测量面形中的回程误差的装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103852030B (zh) | 2016-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Review of phase measuring deflectometry | |
Hu et al. | Absolute three-dimensional micro surface profile measurement based on a Greenough-type stereomicroscope | |
CN103852030B (zh) | 用于倾斜波面非零位干涉***的自由曲面面形重构方法 | |
CN101957182B (zh) | 一种大口径高陡度光学镜面在线测量*** | |
Su et al. | SCOTS: a reverse Hartmann test with high dynamic range for Giant Magellan Telescope primary mirror segments | |
CN102749143B (zh) | 一种提高夏克-哈特曼波前传感器测量精度的波前重构方法 | |
Zhou et al. | Improved method for rapid shape recovery of large specular surfaces based on phase measuring deflectometry | |
Wang et al. | Systematic error control for deflectometry with iterative reconstruction | |
CN103471521B (zh) | 快速、准确的光学非球面的实时检测方法 | |
CN105318847A (zh) | 基于***建模的非球面非零位环形子孔径拼接方法 | |
CN104165755A (zh) | 光栅剪切波像差检测干涉仪及其检测方法 | |
CN103674493A (zh) | 光栅剪切干涉仪波像差检测的***误差的消除方法 | |
Sun et al. | High-accuracy simultaneous phase extraction and unwrapping method for single interferogram based on convolutional neural network | |
Nguyen et al. | Real-time 3D measurement of freeform surfaces by dynamic deflectometry based on diagonal spatial carrier-frequency pattern projection | |
CN101504276A (zh) | 大口径光学镜面Ronchi光栅定量检测仪 | |
Chen et al. | Bayesian uncertainty evaluation of stitching interferometry for cylindrical surface | |
Rhee et al. | Improved iterative method for wavefront reconstruction from derivatives in grid geometry | |
CN116105633A (zh) | 自由曲面光学镜片检测方法 | |
JP6223368B2 (ja) | 光学要素の幾何学的構造を測定する方法及びツール | |
Fortmeier et al. | Steps towards traceability for an asphere interferometer | |
Wang et al. | Vision ray model based stereo deflectometry for the measurement of the specular surface | |
CN205719273U (zh) | 动态高分辨率波前测量装置 | |
Fan et al. | A deep learning-based approach to solve the height-slope ambiguity in phase measuring deflectometry | |
CN113432731A (zh) | 一种光栅横向剪切干涉波前重建过程中的补偿方法 | |
Liu et al. | Double transmission-mediums based geometric phase analysis for determining the two surface profiles of transparent object |
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 |