CN109239913A - 基于平板相位元件的成像光学系统的设计方法 - Google Patents

基于平板相位元件的成像光学系统的设计方法 Download PDF

Info

Publication number
CN109239913A
CN109239913A CN201810922636.XA CN201810922636A CN109239913A CN 109239913 A CN109239913 A CN 109239913A CN 201810922636 A CN201810922636 A CN 201810922636A CN 109239913 A CN109239913 A CN 109239913A
Authority
CN
China
Prior art keywords
phase
point
imaging optical
face
optical system
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
Application number
CN201810922636.XA
Other languages
English (en)
Other versions
CN109239913B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201810922636.XA priority Critical patent/CN109239913B/zh
Publication of CN109239913A publication Critical patent/CN109239913A/zh
Application granted granted Critical
Publication of CN109239913B publication Critical patent/CN109239913B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012Optical design, e.g. procedures, algorithms, optimisation routines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Studio Devices (AREA)

Abstract

本发明涉及一种基于平板相位元件的成像光学系统的设计方法,包括:建立初始系统;选取特征光线;依次逐点初步构建每个相位面的相位函数,得到一第一基于平板相位元件的成像光学系统;将所述第一基于平板相位元件的成像光学系统再次作为初始系统,经过多次迭代,得到一第二基于平板相位元件的成像光学系统;将所述第二基于平板相位元件的成像光学系统作为初始系统,进行优化,得到最终系统。

Description

基于平板相位元件的成像光学系统的设计方法
技术领域
本发明涉及光学设计领域,尤其涉及一种基于平板相位元件的成像光学系统的设计方法。
背景技术
成像光学系统在日常生活、科学研究、医疗、工业、国防等方面有非常重要的应用。随着科学技术的不断发展,人们对成像系统的要求越来越高,非球面和自由曲面等高级曲面也被逐步引入到成像光学系统设计中。一般来说,成像系统是由传统几何曲面组成的。但是此类系统有一些固有的缺陷:首先,传统基于几何曲面的系统经常使用大块的光学材料或镜面基底材料来实现复杂的光学面形,从而造成系统重量和体积较大;第二,对于紧凑型光学系统设计,光线经常需要在系统中实现大角度折转,从而可能对应于极度非常规的几何面形;第三,采用非对称结构的系统的装调难度极大。而使用非对称的自由曲面会进一步加大装调难度。
相位元件可以通过另一种方式来调控光线。而平板相位元件有重要的应用价值。相对于传统几何曲面元件,使用平板相位元件可以有效降低系统的重量和体积;此外,采用平板相位元件可以大大降低系统的装调难度。常见的平板元件有两类:一类是基于衍射面的平板相位元件,另一类是基于超颖表面的相位元件。
设计由平板相位元件组成的成像系统即是求解得到每一个元件的相位函数。现有的一大类相位求解方法是数值化的相位恢复方法,例如G-S方法、快速傅里叶变换方法、模拟退火法等。但上述方法一般应用于光束整形等非成像系统设计。一般来说,成像系统设计希望得到解析形式的相位函数。因此针对有解析形式相位函数的包含相位面的系统设计一般是先找到一个好的初始结构,然后进行计算机辅助优化。然而对于由平板相位元件组成的成像系统设计,尤其是结构非对称的系统设计,可用的初始结构几乎没有,导致设计非常困难。 Mendes-Lopes等人提出了一种改进的SMS方法,可以根据成像要求实现衍射面的逐点求解。但该方法只适用于使用一到二个相位面控制两到三个波面,且系统限于回转对称。而成像光学系统设计的基本要求是要控制多视场多孔径的光线,且能够实现多曲面以及离轴非对称光学系统的设计。综上,目前尚无有效的设计基于平板相位元件的离轴非对称成像光学系统的方法
发明内容
综上所述,确有必要提供一种基于平板相位元件的成像光学系统的设计方法。
一种基于平板相位元件的成像光学系统设计方法,其包括以下步骤:
步骤S1,建立一初始系统,该初始系统包括多个初始的几何平面,每个平面无相位函数。
步骤S2,在全视场和全孔径范围内选取K条特征光线Ri(i=1,2…K),根据物像关系求解所述K条特征光线在像面处的理想像点。
步骤S3,以初始系统为起点,依次逐点初步构建每个相位面的相位函数,得到一第一基于平板相位元件的成像光学系统
步骤S4,将所述第一基于平板相位元件的成像光学系统再次作为初始系统,经过多次迭代,得到一第二基于平板相位元件的成像光学系统;
步骤S5,将所述第二基于平板相位元件的成像光学系统作为初始系统,进行优化,得到最终系统。
相较于现有技术,本发明提供的基于平板相位元件的成像光学系统设计方法,可以仅使用由简单几何平面组成的初始系统作为设计的输入,根据给定的物像关系或光线映射关系,实现成像系统的快速设计。设计得到的系统可以作为后续优化的良好初始结构。本发明提供的方法克服了现有的基于平板相位元件的光学系统设计方法中初始结构匮乏的问题,可以考虑不同视场不同孔径的光线进行设计,且可以设计离轴非对称系统中的多个平板元件。本发明提供的方法也可以作为一种快速的相位恢复方法。
附图说明
图1为本发明提供的求解特征数据点时特征光线和系统中平板元件的交点示意图。
图2为本发明提供的初始构建相位函数的方法示意图。
具体实施方式
下面根据说明书附图并结合具体实施例对本发明的技术方案进一步详细表述。
本发明提供一种基于平板相位元件的成像光学系统设计方法,其包括以下步骤:
步骤S1,建立一初始系统,该初始系统包括多个初始的几何平面,即每个平面无相位函数。该初始系统中的一个初始平面对应待设计成像光学系统中一个相位面;初始平面应该大致处于曲面最终期望的位置处。可以基本按照要求控制光线的走向。
步骤S2,在全视场和全孔径范围内选取K条特征光线Ri(i=1,2…K),根据物像关系求解所述K条特征光线在像面处的理想像点。
步骤S3,以初始系统为起点,依次逐点初步构建每个相位面的相位函数,得到一第一基于平板相位元件的成像光学系统
步骤S4,将所述第一基于平板相位元件的成像光学系统再次作为初始系统,经过多次迭代,得到一第二基于平板相位元件的成像光学系统;
步骤S5,将所述第二基于平板相位元件的成像光学系统作为初始系统,进行优化,得到最终系统。
步骤S2中,所述TR条特征光线应在全视场和全孔径范围内选取的选取,例如:针对一般的采用圆形孔径的系统,可以根据需求选取TF个视场,并将每个视场的孔径在角度方向上分成TA等份,沿着每个角度的半径方向取 TP个不同的孔径位置,每一个孔径位置对应一条特征光线,这样一共选取了 TR=TF×TA×TP条对应不同视场不同孔径位置的特征光线。
特征光线Ri(i=1,2…TR)与一个待求的相位面Ω在实际三维空间中的交点定义为特征点Pi;相位面的相位函数或者函数图像可认为处于一个三维空间中,定义为相位函数空间。特征点Pi在相位函数空间中的对应的投影点定义为特征点Pi,φ。Pi和Pi,φ有相同的x和y方向上的局部坐标。为了得到一个待求的相位面Ω上特征光线Ri(i=1,2…TR)的所有特征数据点Pi处的相位值(i=1,2…TR),将借助特征光线与待求的相位元件Ω的前一个相位面Ω'及后一个相位面Ω”的交点。请参阅图1,在求解每条特征光线Ri(i=1,2… TR)对应的待求的相位元件Ω上的特征数据点Pi(i=1,2…TR)时,将该特征光线Ri与前一个相位面Ω'的交点定义为该特征光线的起点Si,特征光线Ri与后一个相位面Ω”的交点定义为该特征光线的终点Ei。当待设计的系统与特征光线确定后,该特征光线Ri的起点Si是确定的,且易于通过光线追迹得到。在理想状态下,特征光线Ri从Ω'上的Si射出后,经过Ω上的Pi,交于Ω”上的Ei,并最终交目标面于其理想像点Ii,ideal。如果Ω”是像面,特征光线的终点就是其理想像点Ii,ideal;如果在待求相位面Ω和目标面之间还有其他相位面(共有q个),可将特征光线Ri与q个相位面的交点定义为 E* i(ξ=1,2…q).E* i在各自相位面局部坐标系中的坐标为(x* i,ξ,y* i,ξ,0)。若每个相位面的相位函数分别为φξ,从E* i,ξ到E* i,ξ+1的折射率为nξ。从Pi到Ii,ideal的光程可写为
其中n和n'分别为待求相位面Ω前后的折射率,m是衍射级次(如果相位面为超颖表面,则m=1),λ为波长。根据费马原理,光程OPL对x* i,ξ和y* i,ξ (ξ=1…q)偏微分等于零。由此可以得到E* i,1的坐标。该点即可定义为特征光线的终点Ei。定义ri=(Pi-Si)/|Pi-Si|=[ri,x,ri,y,ri,z]T和ri'=(Ei-Pi)/|Ei-Pi|=[ri,x', ri,y',ri,z']T分别为特征光线Ri在Pi处的单位入射和出射矢量;Ni是Pi处实际三维空间中的单位法向量。考虑系统关于子午面对称的常见情况。由于相位面是平面,同一个相位平面上每一点处的法向量为Ni=[0,sinα,cosα]T,矩阵为从全局坐标系到局部坐标系的转换矩阵。α为相位面相对于x轴全局的旋转角度。根据广义的斯涅耳定律
可以求得相位函数在Pi处的梯度根据微分几何知识,在相位函数空间中,相位函数三维面形在Pi,φ处的法向为
请一并参阅图2。步骤S3中,每个相位面的相位函数的构建方法包括以下步骤:
步骤S31,取定一第一个相位面(j=1)。
步骤S32,在构建第j个相位面时,其他所有相位面的位置以及相位函数均保持当前状态不变。请参阅图2(a),取定一第一条特征光线R1与所述第j 个相位面对应的初始平面在实际三维空间中的交点为特征数据点P1;在相位函数空间中的投影特征数据点为P1,φ。设P1,φ(或P1)处的相位值为φ(P1)。在得到所有特征点处的相位后,将各点处的相位整体增加或降低,从而使相位函数顶点处的相位为0。
步骤S33,请参阅图2(b),在得到第i个特征数据点Pi以及其在相位函数空间中的投影点Pi,φ后,为了使光线能够交像面与理想像点Ii,ideal,可以得出相位函数空间中Pi,φ处的法向Ni,φ以及切平面。
步骤S34,请参阅图2(b),计算从Pi到剩余未使用的特征光线和相位面在实际三维空间中交点的距离。从这些交点中找到和Pi距离最短的交点,定义为Pi+1,对应的光线定义为Ri+1
步骤S35,请参阅图2(c),计算从Pi+1到Pη(η=1…i)的距离。找到距离 Pi+1最近的交点,定义为Qi。其在相位函数空间中的投影点为Qi,φ。在Qi,φ的切平面上找到与Pi+1点xy坐标相同的点,将该点定义为Pi+1,φ。此时将i增加 1。
步骤S36,重复步骤S33到S35,直到计算得到第j个相位面上的所有特征数据点。通过拟合可以得到所述相位函数的表达式,请参阅图2(d)。这里采用综合考虑数据点坐标和法向的拟合方法。此时得到了第j个相位面的相位函数。将j增加1。
步骤S37,重复步骤S32到S36,直到得到所有相位面的相位函数。所述的相位函数的数学形式可以为XY多项式、Zernike多项式、径向基函数的叠加函数,非均匀有理B样条等。
步骤S4中,将所述第一基于平板相位元件的成像光学系统再次作为初始系统,经过多次迭代,得到一第二基于平板相位元件的成像光学系统。包括以下步骤:
步骤S41,取定一第一个相位面(j=1)。
步骤S42,在设计第j个相位面时,其他所有相位面的位置以及相位函数均保持当前状态不变。提取所有特征光线和当前第j个相位面的交点Pi (i=1,2…TR)以及在相位函数空间的投影点Pi,φ(i=1,2…TR)。为了使光线能够交像面与理想像点Ii,ideal(i=1,2…TR),可以得出相位函数空间中Pi,φ处的法向Ni,φ(i=1,2…TR)。通过拟合可以得到所述相位函数的表达式。这里采用综合考虑数据点坐标和法向的拟合方法。此时得到了第j个相位面的相位函数。将j增加1。
步骤S43,重复步骤S42,直到得到所有相位面的相位函数。将此系统作为新的初始系统。
步骤S44,步骤a至步骤c是完整的一轮迭代。重复步骤a至步骤c,直到像质收敛或者达到规定的迭代步数。此时得到一第二基于平板相位元件的成像光学系统;
步骤S5,将所述第二基于平板相位元件的成像光学系统作为初始结构,进行优化,得到最终系统。优化可以在商用光学设计软件中进行。
本发明提供的基于平板相位元件的成像光学系统设计方法具有以下优点:其一,可以仅使用由简单纯几何平面组成的系统作为设计的输入,在给定预期的物像关系或光线映射关系,实现满足一定结构与成像要求的基于平板相位元件的成像光学系统的自动设计,克服了初始结构极度匮乏的问题。其二,设计过程中考虑了来自不同视场不同孔径位置的光线,符合成像系统设计的一般要求。其三、此方法可以用于设计有多个相位面的系统设计,且可以设计非常规的离轴非对称结构,适用范围广,且步骤简单,计算速度快。
另外,本领域技术人员还可在本发明精神内作其它变化,当然这些依据本发明精神所作的变化,都应包含在本发明所要求保护的范围内。

Claims (5)

1.一种基于平板相位元件的成像光学系统设计方法,其包括以下步骤:
步骤S1,建立一初始系统,该初始系统包括多个初始的几何平面,即每个平面无相位函数。该初始系统中的一个初始曲面对应待设计平板相位成像光学系统中一个相位面;初始平面应该大致处于曲面最终期望的位置处。可以基本按照要求控制光线的走向。
步骤S2,在全视场和全孔径范围内选取特征光线,根据物像关系求解所述特征光线在像面处的理想像点。
步骤S3,以初始系统为起点,依次逐点初步构建每个相位面的相位函数,得到一第一基于平板相位元件的成像光学系统
步骤S4,将所述第一基于平板相位元件的成像光学系统再次作为初始系统,经过多次迭代,得到一第二基于平板相位元件的成像光学系统;
步骤S5,将所述第二基于平板相位元件的成像光学系统作为初始系统,进行优化,得到最终系统。
2.如权利要求1所述的基于平板相位元件的成像光学系统设计方法,其特征在于,求解得到的相位函数为解析形式,相位函数的拟合过程同时考虑了数据点的坐标与法向。
3.如权利要求1所述的基于平板相位元件的成像光学系统设计方法,其特征在于,设计过程中逐点求解每条光线处的相位函数处的相位值,并考虑了来自不同视场不同孔径位置的光线。
4.如权利要求1所述的基于平板相位元件的成像光学系统设计方法,其特征在于,初步构建每个相位面相位函数时包括以下步骤:
步骤S31,取定一第一个相位面(j=1)。
步骤S32,在构建第j个相位面时,其他所有相位面的位置以及相位函数均保持当前状态不变。取定一第一条特征光线R1与所述第j个相位面对应的初始平面在实际三维空间中的交点为特征数据点P1;在相位函数空间中的投影特征数据点为P1,φ。设P1,φ(或P1)处的相位值为φ(P1)。在得到所有特征点处的相位后,将各点处的相位整体增加或降低,从而使相位函数顶点处的相位为0。
步骤S33,在得到第i个特征数据点Pi以及其在相位函数空间中的投影点Pi,φ后,为了使光线能够交像面与理想像点Ii,ideal,可以得出相位函数空间中Pi,φ处的法向Ni,φ以及切平面。
步骤S34,计算从Pi到剩余未使用的特征光线和相位面在实际三维空间中交点的距离。从这些交点中找到和Pi距离最短的交点,定义为Pi+1,对应的光线定义为Ri+1
步骤S35,计算从Pi+1到Pη(η=1…i)的距离。找到距离Pi+1最近的交点,定义为Qi。其在相位函数空间中的投影点为Qi,φ。在Qi,φ的切平面上找到与Pi+1点xy坐标相同的点,将该点定义为Pi+1,φ。此时将i增加1。
步骤S36,重复步骤S33到S35,直到计算得到第j个相位面上的所有特征数据点。通过拟合可以得到所述相位函数的表达式。这里采用综合考虑数据点坐标和法向的拟合方法。此时得到了第j个相位面的相位函数。将j增加1。
步骤S37,重复步骤S32到S36,直到得到所有相位面的相位函数。
5.如权利要求1所述的基于平板相位元件的成像光学系统设计方法,其特征在于,进行每个相位面相位函数迭代时包括以下步骤:
步骤S41,取定一第一个相位面(j=1)。
步骤S42,在设计第j个相位面时,其他所有相位面的位置以及相位函数均保持当前状态不变。取所有特征光线和当前第j个相位面的交点Pi(i=1,2…TR)以及在相位函数空间的投影点Pi,φ(i=1,2…TR)。为了使光线能够交像面与理想像点Ii,ideal(i=1,2…TR),可以得出相位函数空间中Pi,φ处的法向Ni,φ(i=1,2…TR)。通过拟合可以得到所述相位函数的表达式。这里采用综合考虑数据点坐标和法向的拟合方法。此时得到了第j个相位面的相位函数。将j增加1。
步骤S43,重复步骤S42,直到得到所有相位面的相位函数。将此系统作为新的初始系统。
步骤S44,重复步骤a至步骤c,直到像质收敛或者达到规定的迭代步数。此时得到一第二基于平板相位元件的成像光学系统。
CN201810922636.XA 2018-08-14 2018-08-14 基于平板相位元件的成像光学系统的设计方法 Active CN109239913B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810922636.XA CN109239913B (zh) 2018-08-14 2018-08-14 基于平板相位元件的成像光学系统的设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810922636.XA CN109239913B (zh) 2018-08-14 2018-08-14 基于平板相位元件的成像光学系统的设计方法

Publications (2)

Publication Number Publication Date
CN109239913A true CN109239913A (zh) 2019-01-18
CN109239913B CN109239913B (zh) 2020-11-10

Family

ID=65071301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810922636.XA Active CN109239913B (zh) 2018-08-14 2018-08-14 基于平板相位元件的成像光学系统的设计方法

Country Status (1)

Country Link
CN (1) CN109239913B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110989075A (zh) * 2019-12-30 2020-04-10 中国科学技术大学 基于光纤单元的焦面板非平面阵列孔均分制造方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040125453A1 (en) * 2002-07-03 2004-07-01 Mcaulay Alastair D. Inverse computation of phase masks for two layer diffractive optic element system
US20100042363A1 (en) * 2008-08-08 2010-02-18 Lpi-Europe, S.L. Imaging optics designed by the simultaneous multiple surface method
CN104062757A (zh) * 2014-06-30 2014-09-24 太原理工大学 一种用于相控阵多光束三维扫描的相位分布设计方法
CN104898274A (zh) * 2014-03-05 2015-09-09 清华大学 应用于光学系统的自由曲面的设计方法
CN106569332A (zh) * 2016-11-21 2017-04-19 成都精密光学工程研究中心 连续型螺旋相位板设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040125453A1 (en) * 2002-07-03 2004-07-01 Mcaulay Alastair D. Inverse computation of phase masks for two layer diffractive optic element system
US20100042363A1 (en) * 2008-08-08 2010-02-18 Lpi-Europe, S.L. Imaging optics designed by the simultaneous multiple surface method
CN104898274A (zh) * 2014-03-05 2015-09-09 清华大学 应用于光学系统的自由曲面的设计方法
CN104062757A (zh) * 2014-06-30 2014-09-24 太原理工大学 一种用于相控阵多光束三维扫描的相位分布设计方法
CN106569332A (zh) * 2016-11-21 2017-04-19 成都精密光学工程研究中心 连续型螺旋相位板设计方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110989075A (zh) * 2019-12-30 2020-04-10 中国科学技术大学 基于光纤单元的焦面板非平面阵列孔均分制造方法
CN110989075B (zh) * 2019-12-30 2021-04-23 中国科学技术大学 基于光纤单元的焦面板非平面阵列孔均分制造方法

Also Published As

Publication number Publication date
CN109239913B (zh) 2020-11-10

Similar Documents

Publication Publication Date Title
US11573942B2 (en) System and method for determining exact location results using hash encoding of multi-dimensioned data
CN108227184B (zh) 自由曲面成像光学系统的设计方法
Lin et al. Constructing iterative non-uniform B-spline curve and surface to fit data points
CN108122277B (zh) 一种建模方法及装置
US10210289B2 (en) Method for designing three-dimensional freeform surface to be used in design of the reflective mirrors
CN105988214A (zh) 具有实出瞳的自由曲面离轴成像系统的设计方法
CN104898274A (zh) 应用于光学系统的自由曲面的设计方法
CN111487766A (zh) 自由曲面离轴三反成像系统的设计方法
Jowers et al. Shape detection with vision: implementing shape grammars in conceptual design
Yuan et al. High-accuracy phase demodulation method compatible to closed fringes in a single-frame interferogram based on deep learning
CN109239913A (zh) 基于平板相位元件的成像光学系统的设计方法
Jin et al. Weakly-supervised single-view dense 3D point cloud reconstruction via differentiable renderer
Han et al. Curved LCD based deflectometry method for specular surface measurement
Rinaldo et al. Inverse and approximation problem for two-dimensional fractal sets
Söderkvist Introductory overview of surface reconstruction methods
Loos et al. Modeling of surfaces with fair reflection line pattern
JP2012215855A (ja) 非中央カタディオプトリック光学系において投影を求める方法およびシステム
Chen et al. Rotation-invariant completion network
Wang et al. B-spline freeform surface tailoring for prescribed irradiance based on differentiable ray-tracing
Li et al. Neural gradient learning and optimization for oriented point normal estimation
Zhang Free-form surface fitting for precision coordinate metrology
Chen et al. 3D symmetry detection by a single image and geometric transformation
Zhong et al. PSR-GAT: Arbitrary point cloud super-resolution using graph attention networks
CN116740300B (zh) 一种基于多模态的素体与纹理融合家具模型重建方法
Mukherjee An error reduced and uniform parameter approximation in fitting of B-spline curves to data points

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