光学学报, 2018, 38 (11): 1122002, 网络出版: 2019-05-09   

宽波段光学系统三级光谱的校正方法 下载: 682次

Calibration Method of Tertiary Spectrum for Broadband Optical Systems
作者单位
长春理工大学光电工程学院, 吉林 长春 130022
摘要
提出了一种子波段缝合光学玻璃组合的选择方法,通过最小二乘法求解复消色差方程,得到了能有效校正三级光谱的最优化玻璃组合。基于该组合设计了一个宽波段长焦距光学系统,该系统的波长范围为0.45~1.014 μm,焦距为400 mm,F数为6,全视场角为10°。系统共有3种光学玻璃,7片球面透镜,调制传递函数(MTF)接近衍射极限。研究结果表明,所提方法可以在复消色差光学系统中实现三级光谱残差的校正。
Abstract
A selection method for the sub-band stitching optical glass combination is proposed. The apochromatic equation is solved by the least square method, and the optimal glass combination which can effectively correct the tertiary spectrum is obtained. A broadband long-focus optical system is designed based on this optical glass combination with a wavelength range of 0.45-1.014 μm, a focal length of 400 mm, an F number of 6, and a full field of view of 10°. This system consists of three kinds of optical glasses and seven spherical lenses, and the modulation transfer function (MTF) of this system is closed to the diffraction limit. The research results show that the calibration of the tertiary-spectral residuals in the apochromatic optical system can be achieved with this proposed method.

1 引言

在消色差系统中,残余色差以二级光谱的形式存在,当对成像质量要求更高时,需要对二级光谱进行校正[1-2]。已校正二级光谱的系统称为复消色差光学系统。在需要更进一步提高光学系统的分辨率时,需要对复消色差光学系统中的三级光谱进行校正。

校正色差可利用基于几何图形和数学方法求解出的最佳玻璃组合[3-7]。这些方法有助于由两种或三种玻璃组成的简单系统的二级光谱的校正。但是,当光学系统变得更复杂时,波段加宽,所使用的光学玻璃种类增多,这些方法不再适用。此外,还有研究者提出使用基因优化算法[8]来优化玻璃选择。尽管这种方法可以得到非常好的结果,但计算量大,即使对于简单的光学系统也需要对大量的数据进行处理。

为了能够快速选择复消色差,得到可对三级光谱进行校正的最佳光学玻璃组合,本文提出了一种子波段缝合的方法,对设计波长进行拆分和缝合,并针对拆分和缝合后产生的子波段进行复消色差以选择最佳的光学玻璃组合。

2 光学玻璃选择方法

单透镜与消色差系统的初级轴向色差曲线与二级光谱曲线分别如图1(a)、(b)所示,其中F为焦点。可以看出,消色差光学系统的轴向色差曲线呈U形,并具有一个拐点,可以使消色差波长对应相同的焦面位置。图1(c)、(d)所示为复消色差透镜系统的典型轴向色差曲线,可以看出,复消色差光学系统的轴向色差曲线有两种形式:一种为S形,该曲线具有两个拐点,可以使得两个边缘波长与中心波长对应同样的焦平面位置,两边缘波长与中心波长之间的波长对应的焦面位置相对复消色差焦平面位置的距离即为三级光谱;另一种为W形,该曲线具有三个拐点,故其对应的三级光谱残余量要小于S形的[9]

校正三级光谱就是在指定的光谱范围内对4个波长校正轴向色差,使得4个波长对应的焦面位置重合[10]。对于图1(d)所示的这种情况,将焦平面向左移动一点就可以产生4个焦点,因此能够产生W形轴向色差曲线的复消色差光学系统,其三级光谱已经得到了校正。

图 1. 光谱。(a)初级色差;(b)二级;(c)三级,S形;(d) 三级,W形

Fig. 1. Spectra. (a) Primary; (b) secondary; (c) tertiary, S-shaped; (d) tertiary, W-shaped

下载图片 查看所有图片

为了能够产生W形的轴向色差曲线,除对边缘波长和中心内波长进行复消色差外,还要考虑其余波长。为能引入其余波长,采用子波段缝合型复消色差校正方法。

2.1 子波段缝合光学玻璃组合的选择方法

j种玻璃构成的薄透镜模型如图2所示,为了使该系统能够对某一波段实现复消色差[11],有

i=1jϕi=1,(1)i=1jϕi1νi=0,(2)i=1jϕiPiνi=0,(3)

式中ϕi为第i个透镜的光焦度;Pi为第i个透镜所使用的材料的相对部分色散系数;νi为第i个透镜所使用材料的色散系数。为了能够实现对多色焦点偏移曲线的控制,需要对整个设计波段进行拆分再缝合,过程如图3所示。设这个波段中含有m个特征波长(λ1,λ2,…,λmλ12<m),相邻的三个特征波长构成一个子波段区间,并将其作为缝合的基本元素。一级缝合是在基础子波段λi-1λiλi+1上引入一个相邻的特征波长λi+2,将新引入的波长作为边缘波长,以原有的中心波长λi和边缘波长λi+1分别作为新生成的范围更大的子波段区间各自的中心波长,生成两个新的子波段区间作为一级缝合的结果,即λi-1λiλi+2λi-1λi+1λi+2。例如,以λ1λ2λ3所构成的子波段作为基础子波段,一级缝合在其基础上加入λ4,一级缝合产生的结果是以λ1λ4作为边缘波长,以λ2λ3作为中心波长的两个子波段区间,即λ1λ2λ4λ1λ3λ4。二级缝合则是在一级缝合的基础之上进一步引入一个相邻的特征波长,并将其作为新的子波段区间的边缘波长,再重复缝合的过程。以此类推直至包含m个特征波长的m级缝合,覆盖整个设计波段。

图 2. 薄透镜系统示意图

Fig. 2. Schematic of thin lens system

下载图片 查看所有图片

图 3. 波段的拆分与缝合

Fig. 3. Splitting and stitching of frequency bands

下载图片 查看所有图片

对全部的光学玻璃重新计算相对于每个基础子波段区间以及缝合产生的新子波段区的相对部分色散P和色散系数ν,同时计算 1νPν的值,分别记为Vξ。则由(2)、(3)式可以得到

Vλ1λ2λ31Vλ1λ2λ32Vλ1λ2λ33Vλ1λ2λ3jVλ1λ2λ41Vλ1λ2λ42Vλ1λ2λ43Vλ1λ2λ4jVλ1λ3λ41Vλ1λ3λ42Vλ1λ3λ43Vλ1λ3λ4jVλm-2λm-1λm1Vλm-2λm-1λm2Vλm-2λm-1λm3Vλm-2λm-1λmjϕ1ϕ2ϕ3ϕj=0000,(4)ξλ1λ2λ31ξλ1λ2λ32ξλ1λ2λ33ξλ1λ2λ3jξλ1λ2λ41ξλ1λ2λ42ξλ1λ2λ43ξλ1λ2λ4jξλ1λ3λ41ξλ1λ3λ42ξλ1λ3λ43ξλ1λ3λ4jξλm-2λm-1λm1ξλm-2λm-1λm2ξλm-2λm-1λm3ξλm-2λm-1λmjϕ1ϕ2ϕ3ϕj=0000,(5)

式中各参数 ξλ1λ2λ3=Pλ1λ2λ3νλ1λ2λ3, ξλ1λ2λ4=Pλ1λ2λ4νλ1λ2λ4, ξλ1λ3λ4=Pλ1λ3λ4νλ1λ3λ4,…, ξλm-2λm-1λm=Pλm-2λm-1λmνλm-2λm-1λm, Vλ1λ2λ3=1νλ1λ2λ3, Vλ1λ2λ4=1νλ1λ2λ4, Vλ1λ3λ4=1νλ1λ3λ4,…, Vλm-2λm-1λm=1νλm-2λm-1λm;Vξ的下角标代表计算波段,上角标代表不同的玻璃。

将(4)、(5)式中由Vξ构成的矩阵分别记为 V¯ξ¯,故(4)、(5)式可合写为

 V¯ξ¯·Φ¯=0,(6)

式中 Φ¯=[ ϕ1ϕ2ϕj]T;T表示矩阵的转置,下同;0为零矩阵。

故(1)式可写为

A¯·Φ¯=1,(7)

式中 A¯=[ 111]。

将(6)、(7)式联立可以得到

 A¯V¯ξ¯·Φ¯=a¯,(8)

式中 a¯=[ 100]T。(8)式即为j块光学玻璃对m(m>j)个波长的复消色差条件。

对所设计波段的拆分和缝合使得方程的数量多于未知数的数量,故(8)式无解。在这种情况下使用最小二乘解作为方程的近似解。令矩阵参量

H¯= A¯V¯ξ¯,(9)

则有

Φ¯^=(H¯T·H¯)-1·H¯T·a¯,(10)

由(10)式可以求解出j种薄透镜的光焦度近似解,这种近似解并不能严格满足m个主要波长上的复消色差条件,即不能使(8)式两边严格相等。此时计算各个子波段的色差残余值,将 Φ¯^代入(6)式中,令矩阵

S¯= V¯ξ¯·Φ¯^,(11)

则方差值D= S¯D表示各个子波长的焦距对于各自中心波长焦距的偏离程度,该值可以用来评估该组光学玻璃的色差校正能力。如果D小于所设计的光学系统的焦深,则可认为该组光学玻璃可消除该系统内产生的色差。

除了对色差的校正能力的考量,还希望玻璃所承担的光焦度在一个合理的范围内。单个透镜如果承担很大光焦度,其表面曲率会很大,会给单色像差的校正和加工带来诸多问题。同时希望该组光学玻璃可产生较小的高级球差和色球差。

2.2 边界条件的计算

对于已知光焦度的一个简单薄透镜,其球差[12]可以表示为

L=-2β/(yϕ2),(12)

式中y为半口孔径;ϕ为光焦度;β为角弥散,其表达式为

β=y3ϕ3[n2-(2n+1)K+(n+2)K2/n]/4(n-1)2,(13)

式中n为折射率;K用来表示透镜的形状,其表达式为

K=r2/(r2-r1),(14)

式中r1r2分别为透镜的前、后表面半径。

由(12)~(14)式可知,单薄透镜的球差是关于折射率和形状的函数。

单薄透镜的球差最小形状为

r2/(r2-r1)=K=n(2n+1)/2(n+2)(15)

此时,球差为只与折射率有关的函数:

L=-n(4n-1)8(n-1)2(n+2)y2ϕ(16)

此形状K下的色球差为

ΔLFC=LF-LC,(17)

式中长波时的球差LC=-2βC/(yϕC2);短波时的球差LF=-2βF/(yϕF2),其中βCβF分别为长波和短波的角弥散,ϕCϕF分别为长波和短波的光焦度。

根据高级球差估计式可知,此时的高级球差[13]

ΔL'Sn=L0.707-Lm2,(18)

式中L0.707为0.707孔径处的球差;Lm为边缘孔径处的球差。

通过计算可以获得j个薄透镜组中每个透镜的光焦度 Φ¯^、多色焦点偏离程度D、高级球差ΔL'Sn以及色球差ΔLFC。合适组合中单个透镜的光焦度不应过大,残余的色差值应小于系统焦深,高级球差和色球差应尽量小。

2.3 玻璃组合的选择步骤

要明确光学系统的使用波段和有效焦距(EFL)等设计参数,再根据经济性等因素对玻璃库进行初步的筛选,剔除不常见玻璃或成本较高的玻璃。将筛选后的玻璃库中的玻璃按照每一种组合中含有j种光学玻璃进行组合,并将所有可能的组合进行编号。对设计波段进行拆分和缝合,对每一个子波段计算每一种光学玻璃的Pν值,利用(11)式得到光焦度分配 Φ¯^,令 F¯=i=1jΦ¯^i, F¯表示各个元件的光焦度绝对值之和。若 F¯值过大,则认为该组合的光焦度分配不合理。对光焦度分配合理的组合计算D值,将D值乘以有效焦距(EFL)判断该乘积是否小于光学系统焦深,若是则将结果保留。对所有可能的光学玻璃组合进行计算,再对满足 Φ¯^D值要求的组合计算ΔL'Sn与ΔLFC。选择ΔL'Sn与ΔLFC大小合适的组合进行光学系统设计,该过程如图4所示,其中N为循环次数, F number为相对孔径倒数。

图 4. 玻璃组合选择方法的流程图

Fig. 4. Flow chart of glass combination selection methods

下载图片 查看所有图片

3 设计实例

根据设计流程完成一个宽波段长焦相机镜头的优化设计,在设计过程中只使用球面镜。选择的设计实例参数见表1

表 1. 光学系统要求

Table 1. Specifications for optical system

ParameterRequirement
EFL /mm400
Field of view /(°)±5
F number6
Spectral waveband /μm0.45-1.014
MTFClose to diffraction limit

查看所有表

从给定波段中选择5个特征波长:0.45, 0.5461, 0.6438,0.852,1.014 μm,将0.6438 μm作为EFL的中心波长。计算所采用的光学玻璃库为从成都光明(CDGM)玻璃库中选取的综合性能优且常年在库的光学玻璃,计算光学玻璃在所有的子波段区间上的Vξ。将材料库中的材料每三个一组进行排列组合,利用(10)式计算每一组光学玻璃中每一种材料的归一化光焦度Φ1Φ2Φ3。计算归一化光焦度绝对值之和 F¯,通过 F¯值的大小可以判断该组透镜的光焦度是否合理,记录 F¯值小于10的玻璃组合,将所使用的光学玻璃牌号写入表2中。当测试完所有的玻璃组合,对满足 F¯值小于10条件的玻璃组合按照(11)式计算残差值,再由残差值求出方差值D。由于方差D值是按照归一化的光焦度进行计算的,在与光学系统的焦深进行比较时需要乘以EFL。D值乘以光学系统的EFL得出的值为d。再利用(17)、(18)式计算ΔLFC与ΔL'Sn的值,结果见表2

表 2. 光学玻璃组合的参数

Table 2. Optical glass combination parameters

NMark 1Mark 2Mark 3Φ1Φ2Φ3F¯d /(10-2 mm)ΔL'Sn/10-4ΔLFC /10-4
70H-FK61H-BAK7H-ZF623.70-2.790.16.593.072.335.89
74H-FK61H-BAK7H-ZF883.70-2.770.086.554.082.346.81
34H-FK61H-ZF62H-BAK84.230.14-3.387.754.402.49-12.5
45H-FK61H-ZF52H-BAK84.270.18-3.467.914.442.49-17.8
119H-FK61H-ZF88H-BAK84.220.11-3.347.676.102.5020.0

查看所有表

表2中的数据按照d值由小到大排列。根据设计要求,计算出该光学系统的焦深为±0.046 mm,前4组的d值都小于焦深,再比较ΔL'Sn与ΔLFC。第1组的ΔL'Sn与ΔLFC明显比其他几组的都小,故使用第1组光学玻璃组合来完成设计。

光学系统图如图5所示,采用7片球面镜片、3种玻璃,其中第1、2、7片使用了牌号为H-FK61的光学玻璃,第3、4、6片使用了牌号为H-BAK7的光学玻璃,第5片使用了牌号为H-ZF62的光学玻璃。

图 5. 光学系统结构示意图

Fig. 5. Structural diagram of optical system

下载图片 查看所有图片

光学系统的成像性能如图6所示,其中RMS代表均方根,T表示子午方向,R表示弧矢方向。由图6(a)知,各个视场中的弥散斑尺寸与艾里斑大小十分接近,这与图6(c)中MTF接近衍射极限相对应,说明该光学系统的成像质量良好。由图6(b)可知,该系统的最大畸变小于0.25%。

图 6. 光学系统的成像性能。(a)点列图;(b)畸变;(c) MTF

Fig. 6. Imaging performances of optical system. (a) Spot diagram; (b) distortion; (c) MTF

下载图片 查看所有图片

图7所示为最佳像点位置随着波长变化的曲线图,最佳像点的最大偏移量为-0.03 mm,小于光学系统本身的焦深±0.046 mm。该曲线图呈W形,与零点出现了四处交点,说明该系统使得四种波长的最佳像点位置重合在一起,实现了三级光谱的校正。

图 7. 多色焦点漂移曲线

Fig. 7. Chromatic focal shift

下载图片 查看所有图片

4 结论

提出了一种针对宽波段长焦距光学系统的光学玻璃选择的优化计算方法。该方法可以快速得出符合三极光谱校正的最优光学玻璃组合,并给出了设计实例。研究结果表明,通过所提计算方法得到的光学玻璃组合实现了0.45~1.014 μm波段的复消色差,并且能够实现W形的多色焦点偏移曲线,即对三级光谱进行了校正,且各个视场的成像质量良好,MTF接近衍射极限。

参考文献

[1] 李玉瑶, 唐勇, 林昭珩, 等. 基于大面阵CCD的复消色差航空相机物镜设计[J]. 光学学报, 2009, 29(11): 3183-3185.

    Li Y Y, Tang Y, Lin Z H, et al. Design of apochromatism objective lens of aerial camera based on large frame array CCD[J]. Acta Optica Sinica, 2009, 29(11): 3183-3185.

[2] 伍雁雄, 吴洪波, 张继真, 等. 亚秒级甚高精度星相机光学系统设计[J]. 中国激光, 2015, 42(7): 0716001.

    Wu Y X, Wu H B, Zhang J Z, et al. Optical system design of star camera with high precision better than second level[J]. Chinese Journal of Lasers, 2015, 42(7): 0716001.

[3] Rayces J L, Rosete-Aguilar M. Selection of glasses for achromatic doublets with reduced secondary color[J]. Proceedings of SPIE, 2000, 4093: 36-46.

[4] Sun W S, Chu C H, Tien C L. Well-chosen method for an optimal design of doublet lens design[J]. Optics Express, 2009, 17(3): 1414-1428.

[5] 曲锐, 邓键, 彭晓乐, 等. 0.4~1.7 μm宽波段大相对孔径光学系统设计[J]. 光学学报, 2015, 35(8): 0822007.

    Qu R, Deng J, Peng X L, et al. 0.4-1.7 μm wideband fast F-number optical system design[J]. Acta Optica Sinica, 2015, 35(8): 0822007.

[6] de Albuquerque B F, Sasian J, de Sousa F L, et al. . Method of glass selection for color correction in optical system design[J]. Optics Express, 2012, 20(13): 13592-13611.

[7] Sigler R D. Glass selection for airspaced apochromats using the Buchdahl dispersion equation[J]. Applied Optics, 1986, 25(23): 4311.

[8] Fang Y C, Tsai C M, Macdonald J, et al. Eliminating chromatic aberration in Gauss-type lens design using a novel genetic algorithm[J]. Applied Optics, 2007, 46(13): 2401-2410.

[9] Maxwell J. Tertiary-spectrum manipulation in apochromats[J]. Applied Optics, 1992, 31(13): 2194-2198.

[10] Greǐsukh G I, Ezhov E G, Kazin S V, et al. Diffraction-refraction corrector of the tertiary spectrum[J]. Journal of Optical Technology, 2010, 77(9): 542-547.

[11] MalacaraD, MalacaraZ. Handbook of optical design[M]. 3rd ed. New York: CRC Press, 2013: 65- 87.

[12] 史密斯. 现代光学工程[M]. 周海宪, 程云芳, 译. 第4版. 北京: 化学工业出版社, 2011: 167-168, 409- 411.

    Smith WJ. Modern optical engineering[M]. Zhou H X, Cheng Y F, Transl. 4th ed. Beijing: Chemical Industry Press, 2011: 167-168, 409- 411.

[13] 袁旭沧. 光学设计[M]. 北京: 北京理工大学出版社, 1988: 103- 105.

    Yuan XC. Optical design[M]. Beijing: Beijing Institute of Technology Press, 1988: 103- 105.

王靖锋, 崔庆丰. 宽波段光学系统三级光谱的校正方法[J]. 光学学报, 2018, 38(11): 1122002. Jingfeng Wang, Qingfeng Cui. Calibration Method of Tertiary Spectrum for Broadband Optical Systems[J]. Acta Optica Sinica, 2018, 38(11): 1122002.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!