CN115127473B - 一种大口径光学元件面形检测装置及检测图像畸变校正方法 - Google Patents
一种大口径光学元件面形检测装置及检测图像畸变校正方法 Download PDFInfo
- Publication number
- CN115127473B CN115127473B CN202210610659.3A CN202210610659A CN115127473B CN 115127473 B CN115127473 B CN 115127473B CN 202210610659 A CN202210610659 A CN 202210610659A CN 115127473 B CN115127473 B CN 115127473B
- Authority
- CN
- China
- Prior art keywords
- point
- optical element
- vector
- surface shape
- polynomial
- 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
- 230000003287 optical effect Effects 0.000 title claims abstract description 91
- 238000001514 detection method Methods 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012937 correction Methods 0.000 title claims abstract description 30
- 239000013598 vector Substances 0.000 claims abstract description 85
- 238000013507 mapping Methods 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 49
- 238000005259 measurement Methods 0.000 claims description 30
- 239000011159 matrix material Substances 0.000 claims description 21
- 230000005540 biological transmission Effects 0.000 claims description 18
- 210000001747 pupil Anatomy 0.000 claims description 15
- 239000003550 marker Substances 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 9
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 abstract description 8
- 238000012360 testing method Methods 0.000 abstract description 4
- 238000013461 design Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
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)
- Holo Graphy (AREA)
Abstract
高精度的非球面镜面形检测是一项关键技术问题,为解决测试结果中存在的畸变问题,本发明基于此提供了一种大口径光学元件面形检测装置及检测图像畸变校正方法,本发明根据待测光学元件的形状选择合适的多项式基函数并生成与之对应的正交矢量多项式集,在光学镜面上粘贴按一定大小及形状分布的标记点,使用生成的正交矢量多项式集求取镜面标记点坐标和标记点图像坐标之间的映射关系,进而将这种映射关系应用于畸变面形上每一个像素点,并通过特定的插值处理,完成镜面面形的畸变校正。本发明提供的方法可广泛应用于大口径光学元件的面形检测中,而畸变校正方法可用于更加普遍的各种成像系统中。
Description
技术领域
本发明属于光学检测领域,具体涉及一种大口径光学元件面形检测装置及检测图像畸变校正方法。
背景技术
大口径光学镜面的检测与制备是当前硬科技的前沿项目之一,无疑对空间远距离成像,激光武器聚焦等一系列应用具有重要的应用价值。而随着镜面口径的增大,对光学检测也提出了更多的要求。
非球面镜在大口径光学系统中的应用广泛,非球面镜的面形误差会直接影响光学系统的成像及后端光学系统的探测精度,所以高精度的非球面镜面形检测是一项关键技术问题。
常用的面形检测方法,如平面镜自准直检测法、子孔径拼接检测法、零位补偿检测法,计算全息检测法等。由于计算全息图原则上能够产生任意形状的波前,非常适合于非球面镜尤其是大口径非球面镜的在位检测。另一方面随着激光直写等技术的发展,最小刻线间隔达微米级的计算全息图的制造业变得容易实现,这使得计算全息检测方法得到了广泛的应用。在使用计算全息图检测离轴非球面时,镜面坐标和干涉条纹图之间存在非线性畸变,这种非线性畸变可以通过减小计算全息图到非球面镜面之间的距离来减小,然而,此时计算图的尺寸也相应的增大,从而对计算全息图的制造精度提出了更高的要求,所以要在计算全息图的尺寸和非线性畸变之间做出权衡。考虑到计算全息的制作成本、加工精度等问题,制作的计算全息图的尺寸通常在某一范围内,这使得通常的离轴非球面面形测试系统中存在一定的畸变。这种非线性畸变使得Zernike多项式代表的像差项在畸变区域上存在耦合,从而影响镜面面形数据的分析,此外,畸变导致的面形数据所处的位置拉伸或压缩也使得进一步的面形修磨加工变得困难。所以高精度的畸变校正也变得非常必要。
传统的畸变校正方法借助于光学设计模型来完成,镜面坐标和计算全息平面之间的坐标对应关系通过光线追迹计算,然而,这种方法未考虑干涉仪镜头和内部光学部件所产生的畸变,其次,理想的光学设计模型和实际的测试系统之间存在差异,这使得这种方法给出的畸变校正结果精度不高。
发明内容
为了解决大口径非球面光学元件面形检测过程中存在的畸变问题,本发明提供了一种大口径光学元件面形检测装置及检测图像畸变校正方法。
本发明的发明构思是:
设计一种计算全息来完成大口径非球面光学元件的测量,其包含三个功能区,分别为透射式测量全息、反射式对准全息以及透射式标记全息。
其中,透射式测量全息完成将从非球面反射过来的光束转化为球面波,由干涉仪进行探测。反射式对准全息完成计算全息和干涉仪之间的对准,而透射式标记全息完成计算全息和待测大口径非球面光学元件之间的对准。
为解决畸变问题,本发明根据待测大口径非球面光学元件的形状选择合适的多项式基函数及其与之对应的正交矢量多项式集。在待测大口径非球面光学元件的光学镜面上制作一定分布大小及形状的标记点,并使用坐标测量设备这些标记点的坐标。干涉仪采集到的干涉条纹图中包含有这些标记点的信息,通过使用生成的正交矢量多项式集拟合镜面坐标和条纹图中标记点坐标之间的映射关系,将这种映射关系应用于畸变面形上每一个像素点,并通过特定的插值处理,从而完成镜面面形的畸变校正。
本发明的技术解决方案如下:
该大口径光学元件面形检测及图像畸变校正方法,包括以下步骤:
1】于待进行面形检测的大口径光学元件上制作标记点;
2】对标记点进行测量;
2.1】点光源发出的光束经计算全息后入射至待进行面形检测的大口径光学元件上,且点光源发出的光束全覆盖大口径光学元件表面,光束的方向和待进行面形检测的大口径光学元件上每一点处的法线方向重合;
2.2】大口径光学元件对入射光进行反射;
2.3】反射光于点光源位置处形成干涉条纹,调整干涉条纹至最稀疏时,测量各标记点映射于干涉条纹处的位置坐标;
3】根据映射关系对测量值进行面形矫正;
3.1】计算各标记点位置坐标与各标记点映射于干涉条纹处位置坐标的S矢量多项式和T矢量多项式集;
3.2】根据计算结果完成面形矫正。
进一步地,所述步骤3.1具体是:
3.1.1】定义在光瞳区域Σ内的两个矢量函数A(x,y)和B(x,y)的内积为
其中Sex是光瞳面积;假设在特定光瞳区域上使用的正交基函数为F(x,y),F(x,y)的梯度矢量多项式集可计算为:
其中,分别为F(x,y)在x方向和y方向的导数。/>为x方向和y方向的单位矢量。
3.1.2】对矢量多项式集进行正交化处理,得到S矢量多项式,其具体计算过程如下:计算矢量多项式集/>与其本身的内积CFF:
对内积CFF进行Cholesky分解得到中间变量Q,其满足:
QTQ=CFF, (4)
求取转换矩阵M:
M=(QT)-1 (5)
使用转换矩阵对多项式集进行转换,得到S矢量多项式:
由于S矢量多项式是由梯度矢量多项式集正交化得出的,因此将其看作是一标量函数φ的梯度,即
3.1.3】求取标量多项式φ的旋度来获取矢量多项式T:
其中为第i个标量多项式φi在x方向和y方向的导数,/>为x方向和y方向的单位矢量。
在求取各标记点在干涉条纹上的质心位置坐标之后,建立如下映射关系:
其中,是畸变图中以像素为单位的点;
3.1.4】使用正交的矢量Zernike多项式S和T来拟合映射矩阵M;
假设使用的S矢量多项式索引为S_ind=[4~J],T矢量多项式索引为T_ind=[4~L],标记点个数为n;假设在一个非畸变点p上,S矢量多项式的第j项在x和y方向的分量分别为Sxp,j和Syp,j,其中j∈S_ind,同样的,T矢量多项式的第j项在x和y方向的分量分别为Txp,j和Typ,j,其中j∈T_ind,则映射关系可以由下式给出:
其中s4~sJ,t4~tJ为拟合系数,使用最小二乘法计算公式(10)计算出映射矩阵M(s2…s13t4…t13);
3.1.5】通过映射矩阵M得出每一个基准点对应的畸变像素点坐标;
假设一基准点映射到畸变像素点后的坐标为p,在该处的面形数据由W(p)表示,而在该点最邻近的四个点的面形已知,假设分别为Wlu,Wru,Wld和Wrd。p点离左上方最近的整数点在行和列两方向的距离分别为α和β,那么W(p)可以由如下公式给出的插值方法给出:
Wp=(1-α)·(1-β)·Wlu+α·(1-β)·Wru+(1-α)·β·Wld+α·β·Wrd (11)。
进一步地,所述步骤3】中,假设在特定光瞳区域上使用的正交基函数为F(x,y),典型的函数为圆形区域上的Zemike多项式,方形区域上的Legendre多项式及Z-square多项式。
进一步地,所述步骤2.1】实施前,先调整点光源与计算全息之间的位置,具体是:点光源发出的光束经计算全息后部分反射回点光源位置处,并形成干涉条纹,调整干涉条纹至最稀疏后开始步骤2.1】。
进一步地,所述步骤2.1】中,点光源发出的光束全覆盖大口径光学元件表面具体是:
于计算全息上设置至少三个标记全息,点光源发出的光束经标记全息后投射在大口径光学元件表面边缘位置处。
或:一种大口径光学元件面形检测及图像畸变校正方法,包括以下步骤:
1】搭建检测光路
1.1】建立计算全息,所述计算全息上设置有用于测量大口径光学元件非球面面形的透射式测量全息,实现计算全息和点光源之间对准的反射式对准全息,以及用于计算全息和大口径光学元件非球面面形对准的透射式标记全息;
1.2】调整点光源和计算全息之间的位置,使得由点光源发出的光束经过对准全息反射,反射后的光束在干涉仪上的条纹最稀疏;
1.3】固定点光源和计算全息的相对位置姿态不变,调整待进行面形检测的大口径光学元件的位置姿态,使得由点光源发出的光束经过透射式测量全息后形成的图像刚好位于非球面光学元件的孔径边缘,此时测量全息4对应的干涉区域上出现干涉条纹;
1.4】精调点光源和计算全息的位置姿态,使得测量全息对应的干涉条纹最稀疏;
2】对标记点进行测量
2.1】于待进行面形检测的大口径光学元件非球面面形上设置多个标记点,测量各标记点映射于干涉条纹处的位置坐标;
3】根据映射关系对测量值进行面形矫正;
3.1】计算各标记点位置坐标与各标记点映射于干涉条纹处位置坐标的S矢量多项式和T矢量多项式集;
3.2】根据计算结果完成面形矫正。
上述步骤中的细节优选算法和前一种方法相同。
本发明的优点在于:
1、与使用光学设计模型进行光线追迹从而实现畸变校正的方法相比,本发明提供的大口径光学元件面形检测装置及检测图像畸变校正方法考虑了干涉仪镜头及干涉仪内部光学元件存在的畸变,消除了光学设计的模型和实际测试系统不一致引入的精度损失。
2、本发明提供的大口径光学元件面形检测装置及检测图像畸变校正方法针对不同形状的光学元件,采用相应的正交多项式基函数来描述元件面形,使得像差分析更加合理可行。
3、使用本发明提供的大口径光学元件面形检测装置及检测图像畸变校正方法产生的正交矢量多项式是完备的,可以描述任意形状和大小分布的畸变,具有很强的适用性。
4、本发明提供的大口径光学元件面形检测装置及检测图像畸变校正方法非常适合于解决大口径光学元件的面形测量中存在的畸变问题,应用前景广泛,值得被大力推广。
附图说明
图1是本发明提供的一种使用计算全息测量非球面光学元件面形的测量示意图。
图2是设计的一种计算全息示意图。
图3是在大口径光学元件上粘贴的标记点分布以及其对的应像方点阵示意图。
图4是获取的面形图像及校正图像。
附图标记说明:
1-点光源,2-计算全息,3-非球面光学元件,4-测量全息,5-对准全息,6-标记全息;7-标记点位置分布,8-标记点成像点阵,9-标记点位置,10-标记点成像位置,11-畸变面形图像,12-校正面形图像。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明,显然所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明的保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
以下结合附图及具体实施例对本发明进行详细说明:
如图1所示,干涉仪发出的光束汇聚于一点,形成点光源1,点光源1发出的光束传输到计算全息2表面,计算全息2对由点光源1发出的球面波进行波面的变换,当光束传输到非球面光学元件3时,光束的方向和非球面光学元件3上每一点处的法线方向重合,光束经由非球面光学元件3反射,再次经过计算全息2,到达点光源1后传输到干涉仪内部发生干涉,从而测量出面形。
为在实际测量中方便调节和操作,本发明设计了一种计算全息图如图2所示。图2中的计算全息分为三个区域,分别为透射式测量全息4,反射式对准全息5以及透射式标记全息6;这不同的全息在不同的区域上刻蚀有不同的刻线,在非球面测量时具有不同的功能。测量全息4用于测量非球面面形,对准全息5用于实现计算全息2和点光源1之间的对准,而标记全息6用于计算全息2和非球面光学元件3之间的对准。
具体地讲,测量全息4用于测量非球面面形,从点光源1发出的光束,经过测量全息4后传播到非球面光学元件3表面时,光束方向与非球面光学元件3表面法线方向一致;点光源1发出的光线入射到反射式对准全息5的光线全覆盖,且经反射式对准全息5反射回的光线在干涉仪上形成稀疏的直条纹;也就是说,通过对准全息5用于实现计算全息2和点光源1之间的对准。
点光源1发出的光线经透射式测量全息4透射出的光线对非球面光学元件3全覆盖;
点光源1发出的光线经透射式标记全息6透射出的光线投射在非球面光学元件3的孔径边缘位置处,透射式标记全息6为至少3个,且经过透射式标记全息6之后投射在非球面光学元件3的孔径边缘的不同位置处。一般情况下,投射出的图像为清晰的十字丝图像,这些该十字丝图像用于计算全息2和非球面光学元件3之间的对准。
使用图1给定的方法测量非球面光学元件3的面形,具体包含以下步骤:
1】搭建检测光路
1.1】建立计算全息,所述计算全息上设置有用于测量大口径光学元件非球面面形的透射式测量全息,实现计算全息和点光源之间对准的反射式对准全息,以及用于计算全息和大口径光学元件非球面面形对准的透射式标记全息;
1.2】调整点光源和计算全息之间的位置,使得由点光源发出的光束经过对准全息反射,反射后的光束在干涉仪上的条纹最稀疏;
1.3】固定点光源和计算全息的相对位置姿态不变,调整待进行面形检测的大口径光学元件的位置姿态,使得由点光源发出的光束经过透射式测量全息后形成的图像刚好位于非球面光学元件的孔径边缘,此时测量全息4对应的干涉区域上出现干涉条纹;
1.4】精调点光源和计算全息的位置姿态,使得测量全息对应的干涉条纹最稀疏;
图1所示的测量结构通常会导致面形测量结果中存在有图像畸变,为实现畸变校正,本发明在非球面光学元件3上粘贴一系列标记点,用以实现镜面坐标和图像坐标的映射。图3中,左图给出粘贴于非球面光学元件3上的标记点位置分布7以及一个标记点位置9。右图为标记点在干涉仪靶面上的成像点阵8以及标记点位置9所对应的标记点成像点阵10。具体如步骤2】对标记点进行测量:
2.1】于待进行面形检测的大口径光学元件非球面面形上设置多个标记点,测量各标记点映射于干涉条纹处的位置坐标;
3】根据映射关系对测量值进行面形矫正;
3.1】计算各标记点位置坐标与各标记点映射于干涉条纹处位置坐标的S矢量多项式和T矢量多项式集;
3.2】根据计算结果完成面形矫正。
当然,根据实际情况,也可以考虑制作标记点后再一并进行校准及测量:
1】于待进行面形检测的大口径光学元件上制作标记点;
2】对标记点进行测量;
2.1】点光源发出的光束经计算全息后入射至待进行面形检测的大口径光学元件上,且点光源发出的光束全覆盖大口径光学元件表面,光束的方向和待进行面形检测的大口径光学元件上每一点处的法线方向重合;
2.2】大口径光学元件对入射光进行反射;
2.3】反射光于点光源位置处形成干涉条纹,调整干涉条纹至最稀疏时,测量各标记点映射于干涉条纹处的位置坐标;
3】根据映射关系对测量值进行面形矫正;
3.1】计算各标记点位置坐标与各标记点映射于干涉条纹处位置坐标的S矢量多项式和T矢量多项式集;
3.2】根据计算结果完成面形矫正。
为实现标记点位置分布7以及标记点成像点阵10之间准确的映射关系,须产生正交的S矢量多项式和T矢量多项式集,其具体方法为:
定义在光瞳区域Σ内的两个矢量函数A(x,y)和B(x,y)的内积为:
其中Sex是光瞳面积。
假设在特定光瞳区域上使用的正交基函数为F(x,y),典型的函数为圆形区域上的Zernike多项式,方形区域上的Legendre多项式及Z-square多项式等。F(x,y)的梯度矢量多项式集可计算为:
对矢量多项式集进行正交化处理,得到S矢量多项式,其具体计算过程如下:
(1)计算矢量多项式集与其本身的内积CFF:
(2)对内积CFF进行Cholesky分解得到中间变量Q,其满足:
QTQ=CFF, (4)
(3)求取转换矩阵M:
M=(QT)-1, (5)
(4)使用转换矩阵对多项式集进行转换,得到S矢量多项式:
由于S矢量多项式是由梯度矢量多项式集正交化得出的,所以可看作是一标量函数φ的梯度,即
为了获取完备的矢量多项式集,求取标量多项式φ的旋度来得到矢量多项式T:
在求取这些标记点在干涉仪探测器图像上的质心位置坐标之后,建立如下映射关系:
其中,是畸变图中以像素为单位的点,假设其坐标为(x′,y′),/>是镜面上以mm为单位的,假设其坐标为(x,y)。映射关系M则可以由已知的/>和集合拟合得出。使用正交的矢量Zemike多项式S和T来拟合映射矩阵M。假设使用的S矢量多项式索引为S_ind=[4~J],T矢量多项式索引为T_ind=[4~L],标记点个数为n。假设在一个非畸变点p上,S矢量多项式的第j项在x和y方向的分量分别为Sxp,j和Syp,j,其中j∈S_ind,同样的,T矢量多项式的第j项在x和y方向的分量分别为Txp,j和Typ,j,其中j∈T_ind。则映射关系可以由下式给出:
其中s4~sJ,t4~tJ为拟合系数,使用最小二乘法计算公式(10),便能计算出映射矩阵M(s2…s13 t4…t13)。
在得到映射矩阵M之后,便能得出每一个基准点对应的畸变像素点坐标。假设一基准点映射到畸变像素点后的坐标为p,在该处的面形数据由W(p)表示,而在该点最邻近的四个点的面形已知,假设分别为Wlu,Wru,Wld和Wrd。p点离左上方最近的整数点在行和列两方向的距离分别为α和β,那么W(p)可以由如下公式给出的插值方法给出:
Wp=(1-α)·(1-β)·Wlu+α·(1-β)·Wru+(1-α)·β·Wld+α·β·Wrd (11)
图4为使用以上处理方法获取的结果,包括使用干涉仪测量出的非球面光学元件3的畸变面形图像11以及校正面形图像12。
Claims (6)
1.一种大口径光学元件面形检测及图像畸变校正方法,其特征在于:
1】于待进行面形检测的大口径光学元件上制作标记点;
2】对标记点进行测量;
2.1】点光源发出的光束经计算全息后入射至待进行面形检测的大口径光学元件上,且点光源发出的光束全覆盖大口径光学元件表面,光束的方向和待进行面形检测的大口径光学元件上每一点处的法线方向重合;
2.2】大口径光学元件对入射光进行反射;
2.3】反射光于点光源位置处形成干涉条纹,调整干涉条纹至最稀疏时,测量各标记点映射于干涉条纹处的位置坐标;
3】根据映射关系对测量值进行面形矫正;
3.1】计算各标记点位置坐标与各标记点映射于干涉条纹处位置坐标的S矢量多项式和T矢量多项式集,具体为:
3.1.1】定义在光瞳区域Σ内的两个矢量函数A(x,y)和B(x,y)的内积为
其中Sex是光瞳面积;假设在特定光瞳区域上使用的正交基函数为F(x,y),F(x,y)的梯度矢量多项式集可计算为:
其中,分别为F(x,y)在x方向和y方向的导数,/>为x方向和y方向的单位矢量;
3.1.2】对矢量多项式集进行正交化处理,得到S矢量多项式,其具体计算过程如下:计算矢量多项式集/>与其本身的内积CFF:
对内积CFF进行Cholesky分解得到中间变量Q,其满足:
QTQ=CFF, (4)
求取转换矩阵M:
M=(QT)-1 (5)
使用转换矩阵对多项式集进行转换,得到S矢量多项式:
由于S矢量多项式是由梯度矢量多项式集正交化得出的,因此将其看作是一标量函数φ的梯度,即
3.1.3】求取标量多项式φ的旋度来获取矢量多项式T:
在求取各标记点在干涉条纹上的质心位置坐标之后,建立如下映射关系:
其中,是畸变图中以像素为单位的点;
3.1.4】使用正交的矢量Zernike多项式S和T来拟合映射矩阵M;
假设使用的S矢量多项式索引为S_ind=[4~J],T矢量多项式索引为T_ind=[4~L],标记点个数为n;假设在一个非畸变点p上,S矢量多项式的第j项在x和y方向的分量分别为Sxp,j和Syp,j,其中j∈S_ind,同样的,T矢量多项式的第j项在x和y方向的分量分别为Txp,j和Typ,j,其中j∈T_ind,则映射关系可以由下式给出:
其中s4~sJ,t4~tJ为拟合系数,使用最小二乘法计算公式(10)计算出映射矩阵M(s2 …s13 t4 … t13);
3.1.5】通过映射矩阵M得出每一个基准点对应的畸变像素点坐标;
假设一基准点映射到畸变像素点后的坐标为p,在该处的面形数据由W(p)表示,而在该点最邻近的四个点的面形已知,假设分别为Wlu,Wru,Wld和Wrd,p点离左上方最近的整数点在行和列两方向的距离分别为α和β,那么W(p)可以由如下公式给出的插值方法给出:
Wp=(1-α)·(1-β)·Wlu+α·(1-β)Wru+(1-α)·β·Wld+α·β·Wrd (11);
3.2】根据计算结果完成面形矫正。
2.根据权利要求1所述的大口径光学元件面形检测及图像畸变校正方法,其特征在于:所述步骤3】中,假设在特定光瞳区域上使用的正交基函数为F(x,y),典型的函数为圆形区域上的Zernike多项式,方形区域上的Legendre多项式及Z-square多项式。
3.根据权利要求2所述的大口径光学元件面形检测及图像畸变校正方法,其特征在于:所述步骤2.1】实施前,先调整点光源与计算全息之间的位置,具体是:点光源发出的光束经计算全息后部分反射回点光源位置处,并形成干涉条纹,调整干涉条纹至最稀疏后开始步骤2.1】。
4.根据权利要求3所述的大口径光学元件面形检测及图像畸变校正方法,其特征在于:所述步骤2.1】中,点光源发出的光束全覆盖大口径光学元件表面具体是:
于计算全息上设置至少三个标记全息,点光源发出的光束经标记全息后投射在大口径光学元件表面边缘位置处。
5.一种大口径光学元件面形检测及图像畸变校正方法,其特征在于:
1】搭建检测光路
1.1】建立计算全息,所述计算全息上设置有用于测量大口径光学元件非球面面形的透射式测量全息,实现计算全息和点光源之间对准的反射式对准全息,以及用于计算全息和大口径光学元件非球面面形对准的透射式标记全息;
1.2】调整点光源和计算全息之间的位置,使得由点光源发出的光束经过对准全息反射,反射后的光束在干涉仪上的条纹最稀疏;
1.3】固定点光源和计算全息的相对位置姿态不变,调整待进行面形检测的大口径光学元件的位置姿态,使得由点光源发出的光束经过透射式测量全息后形成的图像刚好位于非球面光学元件的孔径边缘,此时测量全息4对应的干涉区域上出现干涉条纹;
1.4】精调点光源和计算全息的位置姿态,使得测量全息对应的干涉条纹最稀疏;
2】对标记点进行测量
2.1】于待进行面形检测的大口径光学元件非球面面形上设置多个标记点,测量各标记点映射于干涉条纹处的位置坐标;
3】根据映射关系对测量值进行面形矫正;
3.1】计算各标记点位置坐标与各标记点映射于干涉条纹处位置坐标的S矢量多项式和T矢量多项式集,具体为:
3.1.1】定义在光瞳区域Σ内的两个矢量函数A(x,y)和B(x,y)的内积为
其中Sex是光瞳面积;假设在特定光瞳区域上使用的正交基函数为F(x,y),F(x,y)的梯度矢量多项式集可计算为:
其中,分别为F(x,y)在x方向和y方向的导数,/>为x方向和y方向的单位矢量;
3.1.2】对矢量多项式集进行正交化处理,得到S矢量多项式,其具体计算过程如下:计算矢量多项式集/>与其本身的内积CFF:
对内积CFF进行Cholesky分解得到中间变量Q,其满足:
QTQ=CFF, (4)
求取转换矩阵M:
M=(QT)-1 (5)
使用转换矩阵对多项式集进行转换,得到S矢量多项式:
由于S矢量多项式是由梯度矢量多项式集正交化得出的,因此将其看作是一标量函数φ的梯度,即
3.1.3】求取标量多项式φ的旋度来获取矢量多项式T,其中 为第i个标量多项式φi在x方向和y方向的导数,/>为x方向和y方向的单位矢量;
在求取各标记点在干涉条纹上的质心位置坐标之后,建立如下映射关系:
其中,是畸变图中以像素为单位的点;
3.1.4】使用正交的矢量Zernike多项式S和T来拟合映射矩阵M;
假设使用的S矢量多项式索引为S_ind=[4~J],T矢量多项式索引为T_ind=[4~L],标记点个数为n;假设在一个非畸变点p上,S矢量多项式的第j项在x和y方向的分量分别为Sxp,j和Syp,j,其中j∈S_ind,同样的,T矢量多项式的第j项在x和y方向的分量分别为Txp,j和Typ,j,其中j∈T_ind,则映射关系可以由下式给出:
其中s4~sJ,t4~tJ为拟合系数,使用最小二乘法计算公式(10)计算出映射矩阵M(s2 …s13 t4 … t13);
3.1.5】通过映射矩阵M得出每一个基准点对应的畸变像素点坐标;
假设一基准点映射到畸变像素点后的坐标为p,在该处的面形数据由W(p)表示,而在该点最邻近的四个点的面形已知,假设分别为Wlu,Wru,Wld和Wrd,p点离左上方最近的整数点在行和列两方向的距离分别为α和β,那么W(p)可以由如下公式给出的插值方法给出:
Wp=(1-α)·(1-β)Wlu+α·(1-β)·Wru+(1-α)·β·Wld+α·βWrd (11);
3.2】根据计算结果完成面形矫正。
6.根据权利要求5所述的大口径光学元件面形检测及图像畸变校正方法,其特征在于:所述步骤3】中,假设在特定光瞳区域上使用的正交基函数为F(x,y),典型的函数为圆形区域上的Zernike多项式,方形区域上的Legendre多项式及Z-square多项式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210610659.3A CN115127473B (zh) | 2022-05-31 | 2022-05-31 | 一种大口径光学元件面形检测装置及检测图像畸变校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210610659.3A CN115127473B (zh) | 2022-05-31 | 2022-05-31 | 一种大口径光学元件面形检测装置及检测图像畸变校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115127473A CN115127473A (zh) | 2022-09-30 |
CN115127473B true CN115127473B (zh) | 2023-12-08 |
Family
ID=83378967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210610659.3A Active CN115127473B (zh) | 2022-05-31 | 2022-05-31 | 一种大口径光学元件面形检测装置及检测图像畸变校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115127473B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117346687B (zh) * | 2023-12-04 | 2024-02-13 | 中国科学院长春光学精密机械与物理研究所 | 一种修正干涉仪测量反射镜面形误差数据点的方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005156434A (ja) * | 2003-11-27 | 2005-06-16 | Fujinon Corp | 計算機ホログラムを用いた光波干渉測定方法およびこれを用いた干渉計装置 |
CN101261183A (zh) * | 2008-04-15 | 2008-09-10 | 中国科学院光电技术研究所 | 一种大口径非球面镜检测系统 |
CN101270975A (zh) * | 2008-05-16 | 2008-09-24 | 湖南大学 | 一种用于非球面面形测量的组合式干涉装置 |
CN111912534A (zh) * | 2020-07-28 | 2020-11-10 | 中国科学院西安光学精密机械研究所 | 一种光学系统波前测量装置及方法 |
CN112504117A (zh) * | 2020-12-01 | 2021-03-16 | 中国人民解放军国防科技大学 | 光学面形补偿干涉测量中的投影畸变校正方法、系统及介质 |
-
2022
- 2022-05-31 CN CN202210610659.3A patent/CN115127473B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005156434A (ja) * | 2003-11-27 | 2005-06-16 | Fujinon Corp | 計算機ホログラムを用いた光波干渉測定方法およびこれを用いた干渉計装置 |
CN101261183A (zh) * | 2008-04-15 | 2008-09-10 | 中国科学院光电技术研究所 | 一种大口径非球面镜检测系统 |
CN101270975A (zh) * | 2008-05-16 | 2008-09-24 | 湖南大学 | 一种用于非球面面形测量的组合式干涉装置 |
CN111912534A (zh) * | 2020-07-28 | 2020-11-10 | 中国科学院西安光学精密机械研究所 | 一种光学系统波前测量装置及方法 |
CN112504117A (zh) * | 2020-12-01 | 2021-03-16 | 中国人民解放军国防科技大学 | 光学面形补偿干涉测量中的投影畸变校正方法、系统及介质 |
Non-Patent Citations (2)
Title |
---|
Guang-ming Dai.Nonrecursive determination of orthonormalpolynomials with matrix formulation.OPTICS LETTERS.2007,第32卷(第1期),74-76. * |
Ke-Wei E.Distortion mapping correction in the AIMS primary mirror testing by a computer-generated hologram.Research in Astronomy and Astrophysics.2021,第21卷(第7期),165-1-165-8. * |
Also Published As
Publication number | Publication date |
---|---|
CN115127473A (zh) | 2022-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102997863B (zh) | 一种全口径光学非球面面形误差直接检测系统 | |
CN106840027B (zh) | 光学自由曲面的像散补偿型干涉检测装置与检测方法 | |
CN110726381B (zh) | 一种光学自由曲面全频段像差检测系统及检测方法 | |
US5485275A (en) | Apparatus and method for measuring the error of an apparatus which measure a cylindrical shape using an interferometer | |
CN108489421A (zh) | 一种条纹投影检测平面元件面形方法及装置 | |
WO2018000942A1 (zh) | 一种柱面及柱面汇聚镜的检测方法及装置 | |
CN104655053B (zh) | 基于针孔式点衍射干涉仪球面镜曲率半径测量方法 | |
CN110702009A (zh) | 一种基于逆向Hartmann计算机辅助法的三维测量系统 | |
CN113776455B (zh) | 一种离轴非球面反射镜零位补偿检测非线性误差校正方法 | |
US6008904A (en) | Apparatus and methods for detecting and correcting distortion of interference fringes | |
CN103776388A (zh) | 非球面非零位干涉检测中部分补偿透镜对准方法 | |
CN115127473B (zh) | 一种大口径光学元件面形检测装置及检测图像畸变校正方法 | |
CN108061514B (zh) | 一种利用轴向扫描光干涉法检测非球面的动态建模方法 | |
CN108955565B (zh) | 自由曲面干涉仪中自适应零位补偿器空间距离自标定方法 | |
CN113946041A (zh) | 一种折反式卡塞格林望远镜系统及其偏振像差校正方法 | |
CN115166932A (zh) | 一种大口径长焦距离轴光学系统的光轴调整方法 | |
CN102589472B (zh) | 在球面面形干涉检测中高精度消除调整误差的方法 | |
CN103134443B (zh) | 一种大口径大径厚比反射镜面形自准直检测装置及方法 | |
CN111189386B (zh) | 离轴抛物面反射镜干涉测量面形投影畸变的修正方法 | |
CN108267094B (zh) | 一种基于旋转cgh的非圆柱面干涉拼接测量系统及方法 | |
CN115901192A (zh) | 具有实时对准功能的光学系统波前拼接检测方法及装置 | |
CN113091638B (zh) | 一种标定干涉法测量面形中的回程误差的装置及方法 | |
Chen et al. | Fiducial free correction of mapping distortion in null test of aspheres and freeforms | |
CN115164775A (zh) | 大口径凸非球面反射镜面形检测装置及其检测方法 | |
CN111207679B (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 |