Optical system design of four-channel low light level polarization imaging
-
摘要: 为了提高微光成像质量,实现针对动态目标或变化场景的实时性成像的需求,基于偏振成像原理,并结合孔径分割技术设计了微光偏振实时成像光学系统。光学系统采用共口径四通道阵列结构,将4个待测偏振态分成4个独立的成像通道,通过子孔径成像镜组对探测器靶面进行四象限分割成像,每个偏振通道通过放置起偏状态不同的偏振片获取目标不同偏振态的强度图像,从而实现实时偏振成像。光学系统的设计焦距为100 mm、系统整体F数为1.2、工作波长范围为0.4 μm~0.85 μm、最低工作照度为1×10-3 lx、可输出1 080 pixel全高清图像。系统光学总长为167.5 mm,单通道镜头的MTF值在40 lp/mm处全视场均高于0.57。Abstract: In order to improve the quality of low-light-level(LLL) imaging and achieve the requirement for real-time imaging for dynamic targets or changing scenes, a LLL polarization real-time imaging optical system was designed based on the polarization imaging principle and the division of aperture technique. To achieve simultaneously polarization imaging, the optical system uses an array of common-caliber four-channel configuration. This structure divides four polarization states to be measured into four independent imaging channels, and four-quadrant segmentation imaging is performed on the target surface of the detector through a sub-aperture imaging group set. Each polarized channel obtains an intensity image with different polarization states by placing polarizing plates with different polarization states. The focal length of the optical system is 100 mm, the system whole F number is 1.2, the operating wavelength is 0.4 μm~0.85 μm, the minimum working illumination is 1×10-3 lx, and the 1 080 pixel full high-definition (HD image) can be outputted. The total optical length of the system is 167.5 mm, and the modulation transfer function (MTF) value of the single-channel is higher than 0.57 at 40 lp/mm in full field of view.
-
Keywords:
- polarization imaging /
- low light level imaging /
- division of aperture /
- four-channel
-
引言
在光照环境极其微弱的情况下,光电设备对目标的侦查能力受到空间、时间光谱及灵敏度等因素的限制导致目标与背景对比度降低,成像质量下降而无法完成侦查探测任务。如何改进微光成像光学系统,加强其探测能力和提升成像质量,依然是一个备受关注的研究方向[11]。偏振成像技术可以探测目标表面的偏振信息,因此可以提高光电设备在相同场景下对目标的探测能力[12, 15]。最初,偏振成像系统多为分时型,在采集目标时,需要多次成像,系统引起的误差较大,虚假偏振信息过多。为了能够同时获取动态场景全部偏振信息,实时成像系统应运而生,解决了分时系统存在的问题。Taylor等人[1]最早设计了一套微光电视光学系统,完成了在空域和时域对偏振态的测量;Sean Moultrie等人[2]基于孔径分割技术,结合尼康变焦镜头,设计了子孔径成像镜头组,组合实现了可用于微光条件下的可见光偏振成像系统, 虽然实现了实时偏振成像,但镜头并不是完全匹配微光设计,光谱范围未涵盖微光波段,镜头成像分辨率低。2014年,陈振跃等[8]研究了基于高灵敏度CCD的孔径分割偏振成像系统,实验表明能够提高成像对比度,但系统同样是结合已有变焦镜头,仅设计子通道成像镜头并配以激光辅助照明,缺点是光谱范围太窄且属于主动夜视。2015年,茹志兵等[6]设计了基于三代像增强器和高动态范围数字CCD的分时偏振成像的微光成像系统,设计了相应大口径成像物镜和中继光学耦合透镜, 存在的问题是系统为分时成像系统,且利用三代像增强器,系统整体笨重,偏振探测灵敏度低。近年国外研制出分焦平面的偏振传感器[5, 13],单个像元对应单个偏振通道,3~4个像元及微偏振片组成一个偏振成像单元,成像效率高,但目前这种技术对器件加工工艺水平要求较高。
综上所述,根据已有报道, 在微光偏振成像领域,国内完全自主设计的微光偏振实时成像光学系统较少。本文基于孔径分割技术[7],设计了共口径四通道微光偏振实时成像光学系统,可实现在最低照度1×10-3 lx环境下实时偏振成像。
1 偏振成像系统原理
设计共口径四通道微光偏振实时成像光学系统,其光学系统结构如图 1所示。
通过将同轴微光光学系统在光阑处进行离轴化,使整个光学系统形成4个独立的成像通道,结合偏振成像原理,通过放置4个完全相同的成像透镜组,并在成像光路中放置4个起偏状态不同(0°、45°、90°、135°)的线偏振片进行偏振成像,通过每个成像通道实时获取同一目标不同偏振方向的强度图像。
偏振信息可用Stokes矢量描述,Stokes矢量用4个具有强度量纲的参量来表征偏振状态[13],且可被探测器直接探测,其形式为
$$ \mathit{\boldsymbol{S = }}\left[ \begin{array}{l} {S_0}\\ {S_1}\\ {S_2}\\ {S_3} \end{array} \right] = \left[ \begin{array}{l} I\\ Q\\ U\\ V \end{array} \right] = \left[ \begin{array}{l} E_x^2 + E_y^2\\ E_x^2 - E_y^2\\ 2{E_x}{E_y}\cos \delta \\ 2{E_x}{E_y}\sin \delta \end{array} \right] $$ (1) 式中:I表示成像系统接受到的光波总强度;Q表示水平、垂直2个偏振方向的强度差;U表示2个对角线方向的强度差;V表示左、右旋圆偏振分量的强度差。确定了光束的Stokes矢量之后,偏振度(degree of polarization,DoP)和偏振角(angle of polarization,AoP)可由式(2)计算得到:
$$ \left\{ \begin{array}{l} {\rm{DoP}} = \frac{{\sqrt {{Q^2} + {U^2} + {V^2}} }}{I}\\ {\rm{AoP}} = \frac{1}{2}\arctan \left[ {\frac{U}{Q}} \right] \end{array} \right. $$ (2) 当V=0时,偏振度退化成线偏振度(degree of linear polarization,DoLP)为
$$ {\rm{DoLP}} = \frac{{\sqrt {{{\rm{Q}}^2} + {{\rm{U}}^2}} }}{{\rm{I}}} $$ (3) 光学元件或系统对入射斯托克斯向量的影响可以用4×4穆勒矩阵来描述,其定义如下:
$$ M = \left[ {\begin{array}{*{20}{l}} {{M_{11}}}&{{M_{12}}}&{{M_{13}}}&{{M_{14}}}\\ {{M_{21}}}&{{M_{22}}}&{{M_{23}}}&{{M_{24}}}\\ {{M_{31}}}&{{M_{32}}}&{{M_{33}}}&{{M_{34}}}\\ {{M_{41}}}&{{M_{42}}}&{{M_{43}}}&{{M_{44}}} \end{array}} \right] $$ (4) 入射光学系统的入射辐射量可用斯托克斯矢量Si=[Ii Qi Ui Vi]T表示,经过光学系统后的出射辐射量可以用斯托克斯矢量So=[Io Qo Uo Vo ]T表示。两者直接可以通过穆勒矩阵建立连接,即:
$$ {\mathit{\boldsymbol{S}}_o} = \mathit{\boldsymbol{M}} \cdot {\mathit{\boldsymbol{S}}_i} $$ (5) 则成像设备接收到的强度为
$$ I = {M_{11}}{I_i} + {M_{12}}{Q_i} + {M_{13}}{U_i} + {M_{14}}{V_i} $$ (6) 通过4次改变成像系统的偏振状态,可得到类似于公式(6)的4个公式。它们中的任意2个彼此独立,则可以从4个强度方程计算斯托克斯参数。一般情况下,自然大气背景及目标物对太阳入射的偏振效应中,圆偏振分量极少,故V=0,除V之外所有的斯托克斯参数可以从3个强度方程中计算出来。
对于多个通道实时偏振成像来讲,探测成像示意图如图 2所示。
图 2中FP为分割的视场通道,L为成像透镜组,M表示偏振组件、FPA表示面阵探测器焦面,DoAPi表示分孔径通道。由于横向2个通道视图中重叠,因此图 2中用3个通道来代替显示。则探测偏振强度可写成
$$ X = A \cdot {{{\pmb{S}}}_{{\rm{in}}}} = {\left\{ {\left[ {S_D^{\left( {\rm{1}} \right)}} \right]\left[ {S_D^{\left( N \right)}} \right]} \right\}^{\rm{T}}} \cdot {{{\pmb{S}}}_{{\rm{in}}}} $$ (7) 式中:Sin是未知的Stokes矢量;SD(i)是第i个通道测得的Stokes矢量;N是向量X中的强度测量值的数量。
2 光学系统设计
微光条件下,为了使光电设备得到高质量的强度图像,要求所使用的探测器必须具有高灵敏度、高动态范围以及高增益。同时,由于分孔径成像通过分割探测器靶面的方式会带来分辨率损失问题,因此需选择高分辨率大靶面探测器来提供大的成像分辨率,以保证单通道成像分辨率。
微光机芯采用国产化CIS2003芯片,能够在星光条件下(1×10-3 lx)连续输出清晰真彩色1 080 pixel全高清图像,适应实时观察及低照度的应用需求,由于芯片成像靶面不含有微透镜阵列,因此在系统设计中,不需要实现光路远心要求。
2.1 系统设计参数
系统工作环境最低照度为1×10-3 lx,F数为1.2,中心波长为0.6 μm,另根据所选择的探测器及相关要求最终确定所设计的系统指标参数如表 1所示。
表 1 系统设计参数Table 1. System design parameters参数 指标 工作波段/μm 0.4~0.85 F数 1.2 FOV/(°) ≥16 探测器尺寸/pixel 1 920 H×1 080 V 单个像元尺寸/μm 13 最低工作环境照度/lx 1×10-3 探测器的奈奎斯特频率Nn的计算方法为
$$ {N_n} = \frac{1}{{2\varphi }} = \frac{1}{{2 \times 13 \times {{10}^{ - 3}}}} = 38.46\;\;{\rm{lp/mm}} $$ (8) 通过上式计算,结合微光探测特性,同时为保证分孔径后单个成像通道分辨率,取40 lp/mm作为评价光学系统光学传递函数的最大参考空间截止频率值。
根据选择的CCD探测器横纵像元数及单个像元尺寸,可得到光学系统对角线尺寸为
$$ \mathit{L}{\rm{ = }}\sqrt {{M^2} + {N^2}} d $$ (9) 式中:M、N分别表示探测器的横、纵像素数。探测器像元数为1 920 H×1 080 V pixel,像元大小为13 μm,由公式(6)可计算出像面对角线约为28.63 mm。
系统由4个完全相同的分孔径成像镜组成像于探测器靶面,根据靶面大小对CCD进行四象限分割,探测器靶面分割及子孔径示意图如图 3所示。故单个偏振通道成像所占用CCD外接圆的对角线尺寸约为14.315 mm。
2.2 初始结构搭建
为扩大系统相对孔径,在双高斯结构的基础上,增加两片透镜搭建同轴偏振成像系统,得到初始结构如图 4所示, 系统参数如表 2所示。
表 2 同轴系统参数Table 2. Coaxial system design parameters透镜序号 物距l/mm 像距l′/mm 光焦度φ 第1透镜 -120.179 140 102.696 777 0.008 581 第2透镜 -73.894 095 63.438 309 0.013 409 第3透镜 35.527 512 -32.845 421 -0.030 539 第4透镜 37.433 792 -37.245 051 -0.027 925 第5透镜 -38.743 023 35.426 071 0.023 144 第6透镜 234.072 358 -216.712 40 -0.004 471 第7透镜 -55.746 440 55.837 441 0.016 247 第8透镜 -532.191 340 371.169 284 0.002 226 为实现孔径分割同时保证系统整体体积相对较小,通过向后移动光阑位置到第6透镜前面,初步离轴光阑面设置并简单优化设计,得到共口径偏振成像单通道光学初始结构如图 5所示。
共口径部分由一片单透镜及2个双胶合透镜组成,子孔径成像组由3个单透镜组成。表 3给出光学系统单通道初始结构的各项初级像差,由表可得初始结构球差、彗差、像散均较大,其他初级像在0.01左右,系统需进一步完善设计。
表 3 初始结构初级像差Table 3. Primary aberrations of initial structure球差 彗差 像散 场曲 畸变 1.130 290 -0.038 128 -0.010 110 0.065 645 0.013 133 2.3 设计结果
系统由前置共口径物镜组,孔径分割成像镜头组及后置偏振组件构成。设计得到具体的光学结构如图 6所示,四通道渲染光学结构如图 7所示。
系统实现单个成像通道设计指标为像高14 mm,系统焦距100 mm,4个通道除偏振组件外透镜各参数完全一致,成像引入的差异最小,保证后续信息处理及图像处理阶段的方便、快速性等特点。得到的单个通道凝视视场像高大小如表 4所示。
表 4 各通道凝视视场像高Table 4. Image height in gaze field of each channel偏振通道 设计像高/mm 实际半像高/mm 0° 14 6.955 935 45° 14 6.955 932 90° 14 6.956 289 135° 14 6.955 209 由于系统存在一定量的光学畸变,实际像高小于14 mm,各个通道像高大小与设计值最大差值小于0.005 mm。
四通道成像将探测器划分为四象限成像,系统光线追迹示意图如图 8所示,整个系统成像完全对称。
3 像质评价
针对微光夜视系统,像差及像质评价主要考察传递函数、衍射能量集中度、点列斑大小以及相对照度。系统MTF曲线如图 9所示,在奈奎斯特频率为40 lp/mm处所有视场均不小于0.6,成像质量良好。
系统的单个像元大小内衍射能量集中度曲线如图 10所示。从图 10中可以看出,系统在单个探测器像元范围内的衍射能量集中度大于0.9,能量集中度较高。
微光成像光学系统本身工作环境光照环境差,因此对系统相对照度有一定要求,图 11为单通道光学系统相对照度曲线图。从图 11中可以看出,系统单通道相对照度各个视场均大于75%,可以正常成像,且像质良好。
系统探测器单个像元尺寸大小为13 μm,最终设计结果各通道点列斑最大尺寸均小于3倍的像元大小且一致性较好,具体参数如表 5所示。
表 5 各个通道点列斑大小Table 5. Spot size for each channel point通道 像高H×V/mm (0×0)
0视场(1.75×3.0)
0.5视场(2.5×4.2)
0.7视场(2.9×4.8)
0.8视场(3.5×6.0)
全视场0°偏振通道 均方根半径/μm 5.736 5.348 5.246 5.172 5.087 最大半径/μm 12.719 9.481 9.814 9.627 10.181 45°偏振通道 均方根半径/μm 5.084 5.244 5.345 5.427 5.732 最大半径/μm 10.183 9.81 9.489 9.302 12.654 90°偏振通道 均方根半径/μm 5.749 7.232 8.478 9.359 11.255 最大半径/μm 12.901 18.159 19.493 22.32 26.952 135°偏振通道 均方根半径/μm 5.758 5.492 5.734 5.861 6.307 最大半径/μm 13.018 9.958 11.811 12.581 14.19 光学系统的相对孔径直接决定系统聚光能力。但由于微光光电设备普遍工作在低照度环境,因此无法单一地以相对孔径大小来评价设备聚光能力的强弱,因此国内外用T数来作为另一个衡量标准。T数是针对视距和分辨能力之间转换的一个评判标准,需要根据光学系统对分辨率的要求结合物镜焦距及通过口径的大小,最后得出T数数值。
根据T数定义:
$$ T=\frac{{{f}'}}{D\centerdot \sqrt{\tau }} $$ (10) 式中:f′为物镜焦距;D为物镜有效孔径;Τ为透过率。
可计算出本微光偏振成像光学系统T约为1.342。国军标要求焦距在100 mm以下、相对口径为1/1.296的最低T数值不小于1.25,因此该偏振成像系统的设计满足实际使用要求。
4 成像质量公差分析
光学系统在实际加工和光学装调过程中会存在厚度、间隔、面型及元件的倾斜和偏心等误差,如果误差很大,会严重影响光学系统的成像性能,因此一个完善的设计必须按照目前的主流加工工艺水平,对光学系统进行公差分析,并进行合理的公差分配使系统成像性能达到最佳。
公差分析的方法主要有Rms spot radius法、Rms wavefront法、Diff.MTF法。其中MTF是一种直观全面的像质评价方法,因此利用MTF法,考察系统MTF曲线的下降程度。
传统光学系统会在镜片间装上隔圈,通过修切隔圈保证间距公差,同轴度[14]靠镜筒结构精度保证。本文四通道孔径分光偏振成像系统由于四路完全对称,要保证后端探测器上获取的四幅图像完全一致,则子孔径成像组镜片性能要求完全一致。由于本系统设计思想在于光阑处离轴分孔径成像,因此偏心及倾斜公差对于加工及装调尤为重要,所以必须在设计完成阶段需要对子孔径成像组的加工公差进行分析,给出明确的公差容限范围。另根据实际工程经验,同时需要保证光学系统在公差容限内,系统传递函数全视场在40 lp/mm处大于0.3。
结合设计经验, 对系统进行初始公差分配,表 6为初步光学系统统一公差分配表。
表 6 光学零件/装配统一公差分配表Table 6. Uniform tolerance distribution for optical parts/optical assemblies公差项 TFRN
/friTTH
I/mmTEDX
/mmTEDY
/mmTETX
/(°)TETY
/(°)TSDX
/mmTSDY
/mmTIRX
/mmTIRY
/mmTIRR
/mmTIND 数值 5 0.05 0.05 0.05 0.05 0.05 0.03 0.03 0.05 0.05 0.2 0.001 表 6中TFRN为光圈公差,TIND为折射率公差。
运用光学设计软件Zemax的公差分析功能,设定系统焦面到最后一面透镜的距离为补偿量,按照公差分配表进行衍射MTF均值分析,得到光学系统单通道MTF分布如图 12所示。
从图 12中可以看出在初始分配的公差容限下,系统MTF下降量在0.3左右,传递函数下降严重将导致光学系统最终像质变差,不符合实际工程要求。为确保实际光学系统在公差容限内传递函数在40 lp/mm处大于0.3,因此按照分析结果对部分影响严重的面型及元件间倾斜偏心重新进行公差分配,分配结果如表 7所示,未重新分配面按照表 6公差分配。
表 7 光学零件/装配最终公差分配表Table 7. Final tolerance distribution for optical parts/optical assemblies面序号 初次TIRX(Y)
分配收紧TIRX(Y)
公差初次TFRN
分配收紧TFRN
公差初次TSDX(Y)
分配收紧TSDX(Y)
公差10~12 0.05 0.03 5 4 0.05 0.02 12~15 0.05 0.025 5 3 0.05 0.03 3~5 0.05 0.04 5 4 0.05 0.03 再次运用Zemax的公差分析功能,进行Monte-Carlo分析,得到光学系统单通道MTF分布如图 13所示。
通过公差二次分配,从图 13中可看到最差MTF下降到0.4。具体传递函数分布概率如表 8所示。
表 8 公差分析MTF值Table 8. Tolerance analysis MTF value概率 MTF值 90%≥ 0.570 994 51 80%≥ 0.604 671 56 50%≥ 0.643 934 52 20%≥ 0.671 457 94 10%≥ 0.673 026 90 从表中可以看出在分配的公差容限下,系统MTF下降量在0.1左右,有90%的概率光学系统MTF值在0.57之上,80%的概率MTF值高于0.6046,因此公差分配是合理的,所设计的共口径四通道微光实时偏振成像光学系统符合实际加工水平,满足设计及实际使用要求。
5 结论
本文根据光的偏振特性,结合微光夜视成像原理,按照偏振成像技术的具体指标要求,设计共口径四通道微光偏振实时成像光学系统。系统主要由共口径物镜和4个偏心分孔径成像镜头组成,且各个通道均成像质量较好,能够通过一次曝光实时获取动态目标的同一时刻的4个偏振分量的强度图像,通过偏振解算便可以得到完整的Stokes信息。通过对系统进行像质评价,并进行成像质量公差分析,结果显示设计的光学系统MTF值在40 lp/mm处全视场高于0.57。后续将进行加工制造,并进一步研究分焦平面高分辨率实时偏振探测光学系统。
-
表 1 系统设计参数
Table 1 System design parameters
参数 指标 工作波段/μm 0.4~0.85 F数 1.2 FOV/(°) ≥16 探测器尺寸/pixel 1 920 H×1 080 V 单个像元尺寸/μm 13 最低工作环境照度/lx 1×10-3 表 2 同轴系统参数
Table 2 Coaxial system design parameters
透镜序号 物距l/mm 像距l′/mm 光焦度φ 第1透镜 -120.179 140 102.696 777 0.008 581 第2透镜 -73.894 095 63.438 309 0.013 409 第3透镜 35.527 512 -32.845 421 -0.030 539 第4透镜 37.433 792 -37.245 051 -0.027 925 第5透镜 -38.743 023 35.426 071 0.023 144 第6透镜 234.072 358 -216.712 40 -0.004 471 第7透镜 -55.746 440 55.837 441 0.016 247 第8透镜 -532.191 340 371.169 284 0.002 226 表 3 初始结构初级像差
Table 3 Primary aberrations of initial structure
球差 彗差 像散 场曲 畸变 1.130 290 -0.038 128 -0.010 110 0.065 645 0.013 133 表 4 各通道凝视视场像高
Table 4 Image height in gaze field of each channel
偏振通道 设计像高/mm 实际半像高/mm 0° 14 6.955 935 45° 14 6.955 932 90° 14 6.956 289 135° 14 6.955 209 表 5 各个通道点列斑大小
Table 5 Spot size for each channel point
通道 像高H×V/mm (0×0)
0视场(1.75×3.0)
0.5视场(2.5×4.2)
0.7视场(2.9×4.8)
0.8视场(3.5×6.0)
全视场0°偏振通道 均方根半径/μm 5.736 5.348 5.246 5.172 5.087 最大半径/μm 12.719 9.481 9.814 9.627 10.181 45°偏振通道 均方根半径/μm 5.084 5.244 5.345 5.427 5.732 最大半径/μm 10.183 9.81 9.489 9.302 12.654 90°偏振通道 均方根半径/μm 5.749 7.232 8.478 9.359 11.255 最大半径/μm 12.901 18.159 19.493 22.32 26.952 135°偏振通道 均方根半径/μm 5.758 5.492 5.734 5.861 6.307 最大半径/μm 13.018 9.958 11.811 12.581 14.19 表 6 光学零件/装配统一公差分配表
Table 6 Uniform tolerance distribution for optical parts/optical assemblies
公差项 TFRN
/friTTH
I/mmTEDX
/mmTEDY
/mmTETX
/(°)TETY
/(°)TSDX
/mmTSDY
/mmTIRX
/mmTIRY
/mmTIRR
/mmTIND 数值 5 0.05 0.05 0.05 0.05 0.05 0.03 0.03 0.05 0.05 0.2 0.001 表 7 光学零件/装配最终公差分配表
Table 7 Final tolerance distribution for optical parts/optical assemblies
面序号 初次TIRX(Y)
分配收紧TIRX(Y)
公差初次TFRN
分配收紧TFRN
公差初次TSDX(Y)
分配收紧TSDX(Y)
公差10~12 0.05 0.03 5 4 0.05 0.02 12~15 0.05 0.025 5 3 0.05 0.03 3~5 0.05 0.04 5 4 0.05 0.03 表 8 公差分析MTF值
Table 8 Tolerance analysis MTF value
概率 MTF值 90%≥ 0.570 994 51 80%≥ 0.604 671 56 50%≥ 0.643 934 52 20%≥ 0.671 457 94 10%≥ 0.673 026 90 -
[1] TAYLOR M J. A wide field, low light level TV system to measure the state of polarisation of light[J]. Journal of Physics E Scientific Instruments, 1981, 14(7):865. doi: 10.1088/0022-3735/14/7/020
[2] MOULTRIE S, ROCHE M, LOMPADO A, et al. Design of a dual use imager incorporating polarimetric capabilities[J]. SPIE :The International Society for Optical Engineering, 2007, 6682:66820B-1-10. http://cn.bing.com/academic/profile?id=4cc540264a1dab259d39e574ad9c4040&encoded=0&v=paper_preview&mkt=zh-cn
[3] FENG B, SHI Z L, LIU H Z, et al. Calibration method for equivalent extinction ratio of polarized pixel in integrated micropolarizer array camera[J]. SPIE, 2017, 10605: 106051B. http://cn.bing.com/academic/profile?id=ae5ade4a37ce0891841fc9f4bcefc14e&encoded=0&v=paper_preview&mkt=zh-cn
[4] LENIN A S, TYO J S. Design of channeled partial Mueller matrix polarimeters[J]. Journal of the Optical Society of America A-Optics Image Science and Vision, 2016, 33(6):1060. doi: 10.1364/JOSAA.33.001060
[5] TYO J S. Hybrid division of aperture/division of a focal-plane polarimeter for real-time polarization imagery without an instantaneous field-of-view error[J]. Optics Letters, 2006, 31(20):2984-2989. doi: 10.1364/OL.31.002984
[6] 茹志兵, 刘冰, 李双全, 等.基于微光像增强器的偏振成像系统设计与实验[J].应用光学, 2015, 36(3):435-441. doi: 10.5768/JAO201536.0304002 RU Zhibing, LIU Bing, LI Shuangquan, et al. Design and experiment for polarization imaging system of low-light-level image intensifier[J]. Journal of Applied Optics, 2015, 36(3):435-441. doi: 10.5768/JAO201536.0304002
[7] 贺虎成, 季轶群, 周建康, 等.偏心分孔径偏振成像光学系统的设计[J].光学学报, 2013, 33(6):279-284. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gxxb201306040 HE Hucheng, JI Yiqun, ZHOU Jiankang, et al. Optical design of decentered aperture-divided polarization imaging system[J]. Acta Optica Sinica, 2013, 33(6):279-284. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gxxb201306040
[8] 陈振跃, 王霞, 马斌, 等.微光偏振成像系统设计及实验[J].光子学报, 2014, 43(4):411003-04. http://d.old.wanfangdata.com.cn/Periodical/yygx201503018 CHEN Zhenyue, WANG Xia, MA Bin, et al. Low light level polarization imaging system design and experiment[J]. Acta Photonica Sinica, 2014, 43(4):411003-04. http://d.old.wanfangdata.com.cn/Periodical/yygx201503018
[9] TYO J S, GOLDSTEIN D L, CHENAULT D B, et al. Review of passive imaging polarimetry for remote sensing applications.[J]. Applied Optics, 2015, 45(22):5453-5469. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ae70a057befbcc9f36d2a6df7301c1ec
[10] 谭欣明.关于微光物镜T数系列数值的确定[J].云光技术, 1998, 30(5):7-13. http://www.cqvip.com/qk/92686X/199805/3308799.html TAN Xinming. About T parameter confirm ation of lens in low-light level night vision device[J]. Journal of Yunnan Optics, 1998, 30(5) :7-13. http://www.cqvip.com/qk/92686X/199805/3308799.html
[11] 刘宇.微光成像探测优化理论与高性能系统技术研究[D].南京: 南京理工大学, 2009. http://cdmd.cnki.com.cn/Article/CDMD-10288-2010097954.htm LIU Yu. Research on low-light-level imaging detection optimization theory and high performance system[D]. Nanjing: Nanjing University of Science & Technology, 2009. http://cdmd.cnki.com.cn/Article/CDMD-10288-2010097954.htm
[12] 高明, 王川, 刘钧, 等.四通道阵列偏振成像自适应光学系统设计[J].激光与红外, 2016, 46(6):731-736. doi: 10.3969/j.issn.1001-5078.2016.06.017 GAO Ming, WANG Chuan, LIU Jun, et al. Design of four-channel array polarized imaging adaptive optical system[J]. Laser & Infrared, 2016, 46(6):731-736. doi: 10.3969/j.issn.1001-5078.2016.06.017
[13] ALENIN A S, VAUGHN I J, SCOTT T J. Optimal bandwidth micropolarizer arrays[J]. Optics Letters, 2017, 42(3):458. doi: 10.1364/OL.42.000458
[14] 王春艳, 姜会林, 王陆, 等.用于火控动态性能测试的多光轴光学系统[J].兵工学报, 2011, 32(6):746-751. http://d.old.wanfangdata.com.cn/Periodical/bgxb201106019 WANG Chunyan, JIANG Huilin, WANG Lu, et al. The multi-axial optical system study for fire control dynamic performance testing[J]. Acta Armamentarii, 2011, 32(6):746-751. http://d.old.wanfangdata.com.cn/Periodical/bgxb201106019
[15] 黄诗喆, 田继文, 薛庆增.红外偏振成像系统光学设计[J].应用光学, 2018, 39(3): 316-320. doi: 10.5768/JAO201839.0301003 HUANG Shizhe, TIAN Jiwen, XUE Qingzeng. Optical design of infrared polarization imaging system[J]. Journal of Applied Optics, 2018, 39(3): 316-320. doi: 10.5768/JAO201839.0301003