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

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

Info

Publication number
CN109239913B
CN109239913B CN201810922636.XA CN201810922636A CN109239913B CN 109239913 B CN109239913 B CN 109239913B CN 201810922636 A CN201810922636 A CN 201810922636A CN 109239913 B CN109239913 B CN 109239913B
Authority
CN
China
Prior art keywords
phase
plane
point
imaging optical
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.)
Active
Application number
CN201810922636.XA
Other languages
English (en)
Other versions
CN109239913A (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

Images

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方法、快速傅里叶变换方法、模拟退火法等。但上述方法一般应用于光束整形等非成像系统设计。一般来说,成像系统设计希望得到解析形式的相位函数。因此针对有解析形式相位函数的包含相位面的系统设计一般是先找到一个好的初始结构,然后进行计算机辅助优化。然而对于由平板相位元件组成的成像系统设计,尤其是结构非对称的系统设计,可用的初始结构几乎没有,导致设计非常困难。
Figure BDA0001764624060000011
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的光程可写为
Figure BDA0001764624060000041
其中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,矩阵
Figure BDA0001764624060000042
为从全局坐标系到局部坐标系的转换矩阵。α为相位面相对于x轴全局的旋转角度。根据广义的斯涅耳定律
Figure BDA0001764624060000043
可以求得相位函数在Pi处的梯度
Figure BDA0001764624060000044
根据微分几何知识,在相位函数空间中,相位函数三维面形在Pi,φ处的法向为
Figure BDA0001764624060000045
Figure BDA0001764624060000051
请一并参阅图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,重复步骤S41至步骤S43,直到像质收敛或者达到规定的迭代步数;此时得到一第二基于平板相位元件的成像光学系统。
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 CN109239913A (zh) 2019-01-18
CN109239913B true 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)

Families Citing this family (1)

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

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104062757A (zh) * 2014-06-30 2014-09-24 太原理工大学 一种用于相控阵多光束三维扫描的相位分布设计方法
CN104898274A (zh) * 2014-03-05 2015-09-09 清华大学 应用于光学系统的自由曲面的设计方法
CN106569332A (zh) * 2016-11-21 2017-04-19 成都精密光学工程研究中心 连续型螺旋相位板设计方法

Family Cites Families (2)

* 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
WO2010016930A2 (en) * 2008-08-08 2010-02-11 Lpi-Europe, S.L. Imaging optics designed by the simultaneous multiple surface method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104898274A (zh) * 2014-03-05 2015-09-09 清华大学 应用于光学系统的自由曲面的设计方法
CN104062757A (zh) * 2014-06-30 2014-09-24 太原理工大学 一种用于相控阵多光束三维扫描的相位分布设计方法
CN106569332A (zh) * 2016-11-21 2017-04-19 成都精密光学工程研究中心 连续型螺旋相位板设计方法

Also Published As

Publication number Publication date
CN109239913A (zh) 2019-01-18

Similar Documents

Publication Publication Date Title
CN108227184B (zh) 自由曲面成像光学系统的设计方法
Doskolovich et al. Designing double freeform surfaces for collimated beam shaping with optimal mass transportation and linear assignment problems
TWI557436B (zh) 自由曲面離軸光學系統的設計方法
EP2630541B1 (en) Methods and systems for creating free space reflective optical surfaces
TWI622813B (zh) 離軸非球面光學系統的設計方法
TWI499798B (zh) 應用於光學系統的自由曲面的設計方法
Deng et al. Interactive design exploration for constrained meshes
CN104516108A (zh) 自由曲面成像系统的设计方法
CN110764255A (zh) 波导显示亚波长衍射光栅的设计与优化方法
TWI667450B (zh) 光學系統自由曲面面形公差分佈的分析方法
Meyron et al. Light in power: a general and parameter-free algorithm for caustic design
CN111487766B (zh) 自由曲面离轴三反成像系统的设计方法
TW201629577A (zh) 三維自由曲面的設計方法
CN109239913B (zh) 基于平板相位元件的成像光学系统的设计方法
Rausch et al. Illumination design for extended sources based on phase space mapping
Bykov et al. Optimal mass transportation problem in the design of freeform optical elements generating far-field irradiance distributions for plane incident beam
Jansari et al. Adaptive enriched geometry independent field approximation for 2D time-harmonic acoustics
Petters Gravity’s action on light
US11294148B2 (en) Method for designing freeform surface imaging systems
US5222025A (en) Method of generating Fresnel mirrors suitable for use with image display systems
JP6348953B2 (ja) 自由空間反射光学表面を生成する方法及びシステム
Zhdanov et al. Problems of the realistic image synthesis in media with a gradient index of refraction
CN109613697B (zh) 一种反射式非球面光学系统设计方法
JP2012215855A (ja) 非中央カタディオプトリック光学系において投影を求める方法およびシステム
Wang et al. B-spline freeform surface tailoring for prescribed irradiance based on differentiable ray-tracing

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