-
如图1所示,基于Offner结构的中波红外编码孔径光谱成像仪主要由望远系统、双光路Offner光栅成像系统、中继成像系统、两个DMD以及制冷型中波红外探测器等几个部分组成[18]。望远系统作为前置成像单元,主要获取远处目标场景信息。DMD-1位于望远系统的焦面位置,同时也与双光路Offner光栅成像系统的物面重合,用来补偿DMD-2所产生的像面倾斜。DMD-2作为编码器件,配合双光路Offner光栅成像系统同时实现对目标场景信息的色散、编码、聚合等多个功能。中继成像系统将双光路Offner光谱成像系统与制冷型中波红外探测器进行光瞳衔接,实现冷光阑严格匹配。根据中波红外目标探测应用需求,编码孔径光谱成像仪的技术指标如表1所示。
图 1 基于Offner结构的编码孔径光谱成像系统光路图
Figure 1. Optical path diagram of Offner coded aperture spectral imaging system
表 1 编码孔径光谱成像仪技术指标
Table 1. Specifications of coded aperture spectral imager
Parameters Value Wavelength range/μm 3-5 F-number 3.0 Focal length/mm 240 Spatial resolution/mrad ≤0.1 Spectral resolution/nm 100 Detector array size 640×512 Detector pixel size/μm 15 DMD array size 1920×1080 DMD micromirror size/μm 10.8 双光路Offner光栅成像系统作为核心部分,主要由一个反射式凸面闪耀光栅和一个球面反射镜组成,其对应的技术指标要求如表2所示[18]。
表 2 双光路Offner光栅成像系统的技术指标
Table 2. Specifications of the dual-pass Offner spectroscopic sytem
Parameters Value Wavelength range/μm 3-5 Object plane size/mm2 9.6×7.68 Object space NA 0.164 Spectral resolution/nm 100 Dispersion distance/mm 2.4 Optical magification 1 -
凸面闪耀光栅的设计主要包括三个模块:宏观层面的传统光学系统设计;三维偏振光追迹算法;微观层面的凸面闪耀光栅微结构设计与优化。具体算法流程如图2所示。
首先,开展宏观层面的传统光学系统设计。以双光路Offner光谱成像系统的技术指标要求作为设计输入参数,如工作波长范围、物面尺寸、光谱分辨率、数值孔径、色散距离、放大倍率等;以物方远心和像面大小为约束条件,利用点列图和调制传递函数(MTF)来评估优化过程中系统的成像质量。通过传统的光学系统设计过程,可以获得具有良好成像质量、光谱分辨率和色散距离均满足要求的双光路Offner光谱成像系统。该优化设计过程对凸面闪耀光栅的设计具有两个作用:(1)直接获得了凸面闪耀光栅的曲率半径、衍射级次、光栅周期、有效口径等重要技术参数;(2)为三维偏振光追迹算法提供了准确的光学结构参数。
其次,采用三维偏振光追迹算法,对完成优化设计的双光路Offner光谱成像系统进行光线追迹,以获得不同视场的不同光线到达凸面闪耀光栅表面时的入射角、偏振角、方位角等入射光线信息。将这些信息作为微观层面的凸面闪耀光栅微结构设计的关键输入参数。因此,三维偏振光追迹算法是连接宏观层面的传统光学设计与微观层面的凸面闪耀光栅微结构设计的关键纽带。
最后,进行微观层面的凸面闪耀光栅微结构设计和优化。将光栅周期、衍射级次等由宏观层面的传统光学设计方法获得的参数,以及通过三维偏振光追迹算法获得的入射光线信息共同作为设计输入参数,通过严格耦合波分析方法(RCWA)精确计算工作波长范围内的平均衍射效率,结合粒子群优化算法(PSO)不断对凸面闪耀光栅的微结构参数进行迭代优化,以获得收敛的平均衍射效率值和最优的槽型参数。
-
传统的Offner光谱成像系统中凸面光栅主要实现色散分光,而在基于Offner结构的凝视型编码孔径光谱成像系统,凸面光栅同时承担色散分光和合光两个功能。此外,狭缝的取消以及编码器件DMD的引入导致双光路Offner光栅成像系统的光学设计与传统Offner光谱成像系统的设计产生一些区别。如图3(a)所示,假设凸面闪耀光栅的曲率半径为
${R_{\rm{G}}}$ ,其球心为${C_{\rm{G}}}$ 点。球面反射镜的曲率半径为${R_{\rm{M}}}$ ,其球心为${C_{\rm{M}}}$ 点。凸面闪耀光栅与球面反射镜的间隔为${L_{{\rm{GM}}}}$ ,物面尺寸为9.6 mm×7.68 mm。物面中心到球面反射镜的距离为${L_{\rm{O}}}$ ,物面中心的偏心量为$ \Delta H $ ,色散像平面的中心到球面反射镜的距离为${L_{\rm{D}}}$ ,色散距离为$ \delta $ 。为了实现图1所示的光路,DMD-1应该与双光路Offner光栅成像系统的物面重合,而DMD-2则位于色散像平面。如图3(b)所示,在YZ视图中,物平面和色散像平面均绕X轴倾斜12°,这样的设计使得色散分光光路与合光光路关于XY平面对称,合光之后像面与物面完全重合,有利于减小双光路Offner光栅成像系统的像散,特别是X方向的边缘视场。图 3 双光路Offner光栅成像系统的光学设计原理图
Figure 3. Optical design schematic of the dual-pass Offner spectroscopic system
双光路Offner光栅成像系统的初始结构参数应满足如下数学关系[18]:
$$ {R_{\rm{M}}} = 2{R_{\rm{G}}} $$ (1) $$ {R_{\rm{G}}} = \frac{{d\delta }}{{m\left( {{\lambda _2} - {\lambda _1}} \right)}} $$ (2) 式中:
$ d $ 为光栅周期;$ m $ 为光栅衍射级次;$ {\lambda _1} $ 和$ {\lambda _2} $ 分别表示工作波长的最小值和最大值。以表2所示的双光路Offner光谱成像系统的技术指标要求作为设计输入参数,以物方远心和像面大小为约束条件,利用点列图和MTF来评估优化过程中系统的成像质量,设计结果如表3所示。其中,凸面闪耀光栅的半径为120 mm,周期为100 μm,衍射级次为−1级,外形尺寸为102 mm×46 mm。表 3 双光路Offner光栅成像系统的设计结果
Table 3. Design result of the dual-pass Offner spectroscopic sysytem
Parameters Value ${R_{\rm{G}}}$/mm 120 ${R_{\rm{M}}}$/mm 240 ${L_{{\rm{GM}}} }$/mm 117.827 ${L_{\rm{O}}}$/mm 242.826 ${L_{\rm{D}}}$/mm 235.866 $ \Delta H $/mm 41 $ d $/μm 100 $ m $ -1 Length of the grating/mm 102 Width of the grating/mm 46 -
三维偏振光追迹是连接宏观层面的传统光学设计与微观层面的凸面闪耀光栅微结构设计的关键纽带。其功能为微观层面的凸面闪耀光栅微结构设计获取光栅表面光线的入射信息,包括入射角分布、偏振角分布和方位角分布。
三维偏振光追迹算法是传统二维琼斯算法在三维空间的扩展,假设入射光线和出射光线的偏振态分别为
${{\boldsymbol{E}}_{in}}$ 和${{\boldsymbol{E}}_{out}}$ ,光学系统的三维偏振光追迹矩阵为${{\boldsymbol{P}}_{total}}$ ,则相关算法如下[19-20]:$$ {{\boldsymbol{E}}_{out}} = {{\boldsymbol{P}}_{{\text{total}}}} \cdot {{\boldsymbol{E}}_{in}} $$ (3) $$ {{\boldsymbol{P}}_{total}} = \prod\limits_{q = N, - 1}^1 {{{\boldsymbol{P}}_q}} $$ (4) 式中:N为光学表面总数;
${{\boldsymbol{P}}_q}$ 为第$q$ 个光学界面的三维偏振光追迹矩阵,可由下式计算得到:$$ {{\boldsymbol{P}}_q} = {O_{out,q}} \cdot {J_q} \cdot O_{in,q}^{ - 1} $$ (5) $$ {O_{out,q}} = \left( {\begin{array}{*{20}{c}} {\hat s_{x,q}'}&{\hat p_{x,q}'}&{{{\hat k}_{x,q}}}\\ {\hat s_{y,q}'}&{\hat p_{y,q}'}&{{{\hat k}_{y,q}}}\\ {\hat s_{z,q}'}&{\hat p_{z,q}'}&{{{\hat k}_{z,q}}} \end{array}} \right) $$ (6) $$ {J_q} = \left( {\begin{array}{*{20}{c}} {{\alpha _{s,q}}}&0&0 \\ 0&{{\alpha _{p,q}}}&0 \\ 0&0&1 \end{array}} \right) $$ (7) $$ O_{in,q}^{ - 1} = \left( {\begin{array}{*{20}{c}} {{{\hat s}_{x,q}}}&{{{\hat s}_{y,q}}}&{{{\hat s}_{z,q}}} \\ {{{\hat p}_{x,q}}}&{{{\hat p}_{y,q}}}&{{{\hat p}_{z,q}}} \\ {{{\hat k}_{x,q - 1}}}&{{{\hat k}_{y,q - 1}}}&{{{\hat k}_{z,q - 1}}} \end{array}} \right) $$ (8) $$ {\hat s_q} = \dfrac{{{{\hat k}_{q - 1}} \times {{\hat k}_q}}}{{\left| {{{\hat k}_{q - 1}} \times {{\hat k}_q}} \right|}}\text{,}{\hat p_q} = {\hat k_{q - 1}} \times {\hat s_q}\text{,}\hat p_q' = {\hat k_q} \times {\hat s_q}\text{,} $$ (9) 式中:
${J_q}$ 为菲涅耳公式局部坐标系下的三维琼斯矩阵;${O_{out,q}}$ 和${O_{in,q}}$ 为三维坐标变换矩阵;${\alpha _{s,q}}$ 和${\alpha _{p,q}}$ 分别为第$q$ 个界面上S波和P波分量的振幅透射(反射)系数;$\hat s$ 和$\hat p$ 为菲涅耳公式的特征基底矢量;${\hat k_{q - 1}}$ 和${\hat k_q}$ 分别为第q个界面入射光线和出射光线的传播矢量。根据双光路Offner光栅成像系统的设计结果建立三维偏振光追迹模型,计算得到了入射到凸面闪耀光栅上的全部光线信息。如图4所示,图4(a)~(c)分别表示归一化视场为(0,−1)时凸面闪耀光栅入射光线的入射角分布、方位角分布和偏振角分布;图4(d)~(f)分别表示归一化视场为(1,1)时凸面闪耀光栅入射光线的入射角分布、方位角分布和偏振角分布。可见,对于同一视场,入射角分布比较集中,而方位角和偏振角的分布差异较大;而对于不同视场,入射角和方位角的差异较大,偏振角差异较小。
在归一化的二维面视场内进行均匀采样,依次计算每个视场下所有光线的入射角、方位角和偏振角,并统计全部视场下所有入射光线的信息,结果如图5所示。由图5(a)可知,入射角范围为16.35°~22.37°,且出现了21个明显的离散峰值。这是因为在归一化的二维面视场内采样数为21×21,对于双光路Offner光栅成像系统,在凸面闪耀光栅上的入射角对Z方向的视场变化不敏感,而对X方向的视场变化敏感,且每个视场的入射角分布都集中在各自峰值附近很小范围内。由图5(b)和图5(c)可知,方位角范围为58.51°~82.54°,偏振角范围为67.87°~89.53°,且方位角和偏振角的分布都较为离散。因此,传统的凸面闪耀光栅设计方法中,仅仅采用某一固定的入射角、方位角和偏振角来进行衍射效率的评估和槽型的优化设计具有明显的局限性。
-
严格耦合波分析方法是一种广泛应用于平面光栅槽型设计的电磁场理论,它将介质光栅当作周期调制的平面波导处理,通过麦克斯韦方程组和边界条件进行求解波导模式的本征方程[21],其获得的所有解都必须满足Hill微分方程[22]。由于曲率的存在,凸面闪耀光栅表面的介电常数函数并不是严格的周期函数,这打破了使用Hill微分方程的先决条件。文中利用三维偏振光追迹算法克服了这个冲突,将凸面闪耀光栅表面进行微元化,每一个微元表面根据菲涅耳公式建立自身的局部坐标系,并结合三维偏振光追迹算法获得该微元表面的入射光线信息。在局部坐标系中,任意微元表面都近似等效为平面光栅,满足Hill微分方程。
以中心波长作为闪耀波长
$ {\lambda _B} $ ,以能量占比最大的那一类光线的入射角作为入射角$ i $ ,根据闪耀光栅方程可得衍射角$ \theta $ :$$ d\left( {\sin i - sin\theta } \right) = m{\lambda _B} $$ (10) 进一步可得闪耀角的初始值
$ \gamma $ :$$ 2\gamma = \theta - i $$ (11) 假设凸面闪耀光栅的槽型为直角三角形,则可计算得到槽深的初始值
$ h $ :$$ h\left( {\tan \gamma + \dfrac{1}{{\tan \gamma }}} \right) = d $$ (12) 将光栅周期、衍射级次、波长、入射角、方位角、偏振角、闪耀角初始值、槽深初始值等作为设计输入参数,通过RCWA可计算出某一视场下任意入射光线信息对应的衍射效率
$ {\boldsymbol{\eta}} _{m,n}^{\lambda ,F} $ ,上标$ \lambda $ 表示工作波长,$ F $ 表示视场编号;下标m和n分别表示入射光线的光瞳采样序号。则工作波长范围${{\lambda _1} \sim {\lambda _2}}$ 内的平均衍射效率为:$$ {{\boldsymbol{\eta}} _G} = \dfrac{{\displaystyle\sum\limits_{\kappa = 1}^K {\displaystyle\sum\limits_{\lambda = {\lambda _1}}^{{\lambda _2}} {\displaystyle\sum\limits_{m = 1}^M {\displaystyle\sum\limits_{n = 1}^N {{\boldsymbol{\eta}} _{m,n}^{\lambda ,\kappa }} } } } }}{{K \times Y \times m \times n}} $$ (13) 式中:K为视场采样数;Y为工作波长采样数;M和N表示光瞳采样数。以工作波长范围内的平均衍射效率为优化目标函数,结合粒子群优化算法(PSO)不断对凸面闪耀光栅的微结构参数进行迭代优化。
$$ \left\{ {\begin{array}{*{20}{l}} {{\boldsymbol{P}} = \left[ {{{\boldsymbol{X}}^t}} \right],{{\boldsymbol{X}}^t} = \left[ {{\gamma _t},{h_t}} \right]} \\ {{{\boldsymbol{\eta }}_g} = F\left( {\left[ {{\boldsymbol{X}}_g^t} \right]} \right) = \left[ {{\boldsymbol{\eta}} _g^t} \right]} \\ {{{\boldsymbol{\eta}}_t} = F\left( {\left[ {{\boldsymbol{X}}_t^g} \right]} \right) = \left[ {{\boldsymbol{\eta}} _t^g} \right]} \end{array}} \right. $$ (14) 粒子群
$ {\boldsymbol{P}} $ 中每个单独的粒子记为$ {\boldsymbol{X}} $ ,其变量为闪耀角$ \gamma $ 和槽深$ h $ ,$ t $ 为粒子序号,$ g $ 为迭代序数。利用RCWA计算工作波长范围内平均衍射效率的算法表示为F。${{\boldsymbol{\eta}}_g}$ 表示第g次迭代中所有粒子衍射效率的集合,而${{\boldsymbol{\eta}}_t}$ 表示粒子群中第$ t $ 个粒子在所有迭代次数中衍射效率的集合,如公式(14)所示:$$ \left\{\begin{array}{c}{{\boldsymbol{\eta}} }_{g}^{best}=\mathrm{max}\left({\boldsymbol{\eta}}_{g}\right)=F\left({{\boldsymbol{X}}}_{g}^{best}\right)\\ {{\boldsymbol{\eta}}}_{t}^{best}=\mathrm{max}\left({\boldsymbol{\eta}}_{t}\right)=F\left({{\boldsymbol{X}}}_{t}^{best}\right)\end{array}\right. $$ (15) $$ \left\{ {\begin{array}{*{20}{l}} {\Delta _g^t = {\boldsymbol{X}}_{g - 1}^t - {\boldsymbol{X}}_g^t} \\ {\Delta _g^{best} = {\boldsymbol{X}}_g^t - {\boldsymbol{X}}_g^{best}} \\ {\Delta _t^{best} = {\boldsymbol{X}}_g^t - {\boldsymbol{X}}_t^{best}} \end{array}} \right. $$ (16) 如公式(15)~(16),在利用PSO进行优化过程中,有三种中间参量为关键参量:即当前迭代序数g中特定粒子t的衍射效率
${\boldsymbol{ \eta}} _g^t $ ,当前迭代序数g中所有粒子的最大衍射效率${\boldsymbol{ \eta}} _g^{best} $ ,以及特定粒子t在所有迭代次数中的最大衍射效率$ {\boldsymbol{\eta}} _t^{best} $ 。它们对应的变量$ {\boldsymbol{X}}_g^t $ 、$ {\boldsymbol{X}}_g^{best} $ 、$ {\boldsymbol{X}}_t^{best} $ 同时也被记录。此外,这三个变量与上一代的差值分别为$ \Delta _g^t $ 、$ \Delta _g^{best} $ 和$ \Delta _t^{best} $ ,分别与惯性系数$ {k_1} $ 、个体约束因子$ {k_2} $ 、粒子群约束因子$ {k_3} $ 相乘后相加创建下一代的变化。每个粒子$ {\boldsymbol{X}} $ 不断迭代,搜索优化目标函数的最大值,直到平均衍射效率收敛,如公式(17)所示:$$ \left\{ {\begin{array}{*{20}{l}} {\Delta _{g + 1}^t = {k_1}\Delta _g^t + {k_2}\Delta _g^{best} + {k_3}\Delta _t^{best}} \\ \begin{gathered} {\boldsymbol{X}}_{g + 1}^t = {\boldsymbol{X}}_g^t + \Delta _g^t \hfill \\ {{\boldsymbol{P}}_{g + 1}} = [{\boldsymbol{X}}_{g + 1}^t] \hfill \\ \end{gathered} \end{array}} \right. $$ (17) PSO优化算法的初始参数如表4所示,粒子数为20,最大迭代次数为100,惯性系数
$ {k_1} = 0.8 $ ,个体约束因子$ {k_2} = 0.5 $ ,粒子群约束因子$ {k_3} = 0.5 $ ,光谱通道数为20,基底材料为无氧铜,光栅材料为金。PSO优化算法的迭代过程及结果如图6所示,经过16次迭代优化后,平均衍射效率收敛于85.47%。优化后闪耀角为1.1895°,槽深为1.85 μm。表 4 粒子群优化算法的初始参数
Table 4. Initial Parameters of the PSO optimization algorithm
Parameters Value Number of particles 20 Maximum number of generations 100 Inertia factor $ {k_1} $ 0.8 Constriction factor of individual $ {k_2} $ 0.5 Constriction factor of swarm $ {k_3} $ 0.5 Core function RCWA Coating material Gold Substrate material Copper Dimension 2 Number of spectrum channels 20 Variables $ \gamma $ & $ h $ -
目前,凸面闪耀光栅的加工方法主要有电子束直写法[23]、全息离子束刻蚀法[24]、X射线光刻法[25]、机械刻划法[15,26]等。由于优化设计后凸面闪耀光栅的周期较大(100 μm),工作波长范围为3~5 μm,对凸面闪耀光栅的面型误差和表面粗糙度要求并不苛刻,且外形尺寸较大(102 mm×46 mm),因此采用基于单点金刚石车削工艺的机械刻划法[15]。
文中采用美国ATMTEK Precitech公司的Nanoform700 Ultra四轴超精密单点金刚石车床作为凸面闪耀光栅的加工设备。光栅基底材料为无氧铜,采用刀尖角为90°的单晶偏尖刀作为刀具,如图7(a)所示,对凸面光栅基底进行槽型的加工和制备,如图7(b)所示,获得的凸面闪耀光栅样品如图7(c)所示。加工后的凸面闪耀光栅镀制金膜,以在中波红外波段获得最高的衍射效率。
采用白光干涉仪对凸面闪耀光栅的槽型制造误差进行检测,如图8所示,光栅周期为99.945 μm,槽深为1.834 μm,闪耀角为1.1783°,光栅槽面基本平整,脊线较为尖锐,与理论槽型接近。
图 8 凸面闪耀光栅的槽型制造误差检测结果
Figure 8. Test result of groove fabrication error of the convex blazed grating
作为双光路Offner光栅成像系统的核心元件,凸面闪耀光栅同时承担色散分光和合光的功能,入射场景光波将在凸面闪耀光栅上发生两次衍射,因此,凸面闪耀光栅的加工面型误差将对衍射波前产生直接影响,从而降低光谱成像系统的成像质量。利用Zygo干涉仪对凸面闪耀光栅的0级衍射光进行干涉测试,获得了光栅表面的面型误差峰谷值PV为
$ 0.545{\lambda _{test}} $ ,均方根值RMS为$ 0.085{\lambda _{test}} $ ,其中测试波长$ {\lambda _{test}} $ 为632.8 nm,如图9所示。光栅表面产生的波前误差主要表现为像散,这主要是因为加工过程中凸面闪耀光栅沿两个正交方向的半径$ {R_{GC}} $ 和$ {R_{GX}} $ 不完全相等造成的。若以最小工作波长$ {\lambda _1} = $ 3 μm进行换算,可得PV为$ 0.115{\lambda _1} $ ,RMS为$ 0.018{\lambda _1} $ ,满足瑞利判据对波前误差PV小于$ {\lambda \mathord{\left/ {\vphantom {\lambda 4}} \right. } 4} $ 的要求。图 9 凸面闪耀光栅零级光波前形貌图
Figure 9. Topography of zero-order optical wavefront of the convex blazed grating
图10(a)为凸面闪耀光栅衍射效率测量装置原理图,该装置主要由氮化硅红外光源、单色仪、光功率计、球面反射镜以及凸面闪耀光栅组成,其中球面反射镜为双光路Offner光栅成像系统中的主镜,镀制金膜。氮化硅红外光源采用卓立汉光LSSiN-40,其光谱范围覆盖1~14 μm;单色仪采用卓立汉光Omni-
$ \lambda $ 200 i,焦距为200 mm,可在3~5 μm范围内进行波长扫描。首先,利用分光光度计测试获得金膜陪镀片在3~5 μm光谱范围内各波长采样点的反射率$ {r_\lambda } $ 。其次,将单色仪的出射狭缝换为圆孔,并将光功率计放置在单色仪出光口,记录单色仪进行波长扫描时各波长采样点的出射光功率$ {P_{0,\lambda }} $ 。最后,将单色仪出光口与双光路Offner光栅成像系统的物面重合,光功率计放置在Offner系统的色散像面上,再次记录单色仪进行波长扫描时各波长采样点的出射光功率$ {P_{G,\lambda }} $ 。则凸面闪耀光栅在各波长采样点的衍射效率为:$$ {{\boldsymbol{\eta}} _\lambda } = \dfrac{{{P_{G,\lambda }}}}{{r_\lambda ^2 \cdot {P_{0,\lambda }}}} $$ (18) 图10(b)为测试得到的凸面闪耀光栅衍射效率曲线,波长采样间隔为100 nm。可见最高衍射效率在波长3.6 μm附近达到93.46%,起始波长3 μm处衍射效率为77.65%,终止波长5 μm处衍射效率为71.18%,平均衍射效率为84.29%,与前文PSO优化结果较为吻合。
Development of convex blazed grating in coded aperture spectral imager
-
摘要: 针对一款基于DMD的光谱维编码Offner光谱成像仪对凸面闪耀光栅的性能要求,提出了一种凸面闪耀光栅的宏观-微观一体化优化设计方法,利用三维偏振光追迹算法有机融合了宏观层面的Offner光学系统设计与微观层面的凸面闪耀光栅槽型设计。介绍了编码孔径Offner光谱成像系统的组成和工作原理,并结合系统的使用要求设计了一款平均衍射效率为85.47%的中波红外凸面闪耀光栅。在此基础上,采用超精密单点金刚石车床成功制备了曲率半径为120 mm、周期为99.945 μm、闪耀角为1.1783°、槽深为1.834 μm的凸面闪耀光栅。测试结果表明,在3~5 μm光谱范围内,最大衍射效率为93.46%,平均衍射效率为84.29%,与理论设计值较为吻合,验证了凸面闪耀光栅设计方法的有效性。Abstract: Aiming at the performance requirements of a DMD-based Offner spectral imager encoding in spectral dimension for a convex blazed grating, a macro-micro integrated optimization design method for convex blazed gratings was proposed. The three-dimensional polarization ray tracing algorithm was used to organically integrate the optical design of the Offner system in macro-level and the groove design of the convex blazed grating in micro-level. The composition and working principle of the coded aperture Offner spectral imaging system were introduced, and a MWIR convex blazed grating with an average diffraction efficiency of 85.47% was designed according to the requirements of the system. On this basis, a convex blazed grating with curvature radius of 120 mm, grating period of 99.945 μm, blazed angle of 1.1783°, groove depth of 1.834 μm was successfully fabricated by using an ultra-precision single-point diamond lathe. The test results show that in the spectral range of 3-5 μm, the maximum diffraction efficiency is 93.46% and the average diffraction efficiency is 84.29%, which is in good agreement with the theoretical design value. Thus, the proposed design method of the convex blazed grating is verified to be effective and valuable.
-
Key words:
- spectral imaging /
- convex grating /
- diffraction efficiency /
- coded aperture
-
表 1 编码孔径光谱成像仪技术指标
Table 1. Specifications of coded aperture spectral imager
Parameters Value Wavelength range/μm 3-5 F-number 3.0 Focal length/mm 240 Spatial resolution/mrad ≤0.1 Spectral resolution/nm 100 Detector array size 640×512 Detector pixel size/μm 15 DMD array size 1920×1080 DMD micromirror size/μm 10.8 表 2 双光路Offner光栅成像系统的技术指标
Table 2. Specifications of the dual-pass Offner spectroscopic sytem
Parameters Value Wavelength range/μm 3-5 Object plane size/mm2 9.6×7.68 Object space NA 0.164 Spectral resolution/nm 100 Dispersion distance/mm 2.4 Optical magification 1 表 3 双光路Offner光栅成像系统的设计结果
Table 3. Design result of the dual-pass Offner spectroscopic sysytem
Parameters Value ${R_{\rm{G}}}$ /mm120 ${R_{\rm{M}}}$ /mm240 ${L_{{\rm{GM}}} }$ /mm117.827 ${L_{\rm{O}}}$ /mm242.826 ${L_{\rm{D}}}$ /mm235.866 $ \Delta H $ /mm41 $ d $ /μm100 $ m $ -1 Length of the grating/mm 102 Width of the grating/mm 46 表 4 粒子群优化算法的初始参数
Table 4. Initial Parameters of the PSO optimization algorithm
Parameters Value Number of particles 20 Maximum number of generations 100 Inertia factor $ {k_1} $ 0.8 Constriction factor of individual $ {k_2} $ 0.5 Constriction factor of swarm $ {k_3} $ 0.5 Core function RCWA Coating material Gold Substrate material Copper Dimension 2 Number of spectrum channels 20 Variables $ \gamma $ &$ h $ -
[1] Galvis L, Lau D, Ma X, et al. Coded aperture design in compressive spectral imaging based on side information [J]. Applied Optics, 2017, 56(22): 6332-6340. doi: 10.1364/AO.56.006332 [2] Wang L, Xiong Z, Gao D, et al. Dual-camera design for coded aperture snapshot spectral imaging [J]. Applied Optic, 2015, 54(4): 848-858. doi: 10.1364/AO.54.000848 [3] Manolakis D, Marden D, Shaw G A. Hyperspectral image processing for automatic target detection applications [J]. Lincoln Laboratory Journal, 2003, 14(1): 79-116. [4] Kamruzzaman M, Makino Y, Oshita S. Rapid and non-destructive detection of chicken adulteration in minced beef using visible near-infrared hyperspectral imaging and machine learning [J]. Journal of Food Engineering, 2016, 170(7): 8-15. [5] Adam E, Mutanga O, Rugege D. Multispectral and hyperspectral remote sensing for identification and mapping of wetland vegetation: A review [J]. Wetlands Ecology & Management, 2010, 18: 281-296. [6] Rueda H, Arguello H, Arce G R. DMD-based implementation of patterned optical filter arrays for compressive spectral imaging [J]. J Opt Soc Am A, 2015, 32: 80-89. doi: 10.1364/JOSAA.32.000080 [7] Wagadarikar A, John R, Willett R, et al. Single disperser design for coded aperture snapshot spectral imaging [J]. Applied Optics, 2008, 47: B44-B51. doi: 10.1364/AO.47.000B44 [8] Gehm M E, John R, Brady D J, et al. Single-shot compressive spectral imaging with a dual-disperser architecture [J]. Optics Express, 2007, 15: 14013-14027. doi: 10.1364/OE.15.014013 [9] Shao Xiujuan, Hu Bingliang, Yan Peng. Design of multi-spectral adaptive spectrometer based on Hadamard transform [J]. Infrared and Laser Engineering, 2010, 39(5): 963-966. (in Chinese) doi: 10.3969/j.issn.1007-2276.2010.05.039 [10] Wu Y, Mirza I O, Arce G R, et al. Development of a digital micromirror device based multishot snapshot spectral imaging system [J]. Optics Letters, 2011, 36: 2692-2694. doi: 10.1364/OL.36.002692 [11] Goldstein N, Vujkovic-Cvijin P, Fox M J, et al. DMD- based adaptive spectral imagers for hyperspectral imagery and direct detection of spectral signatures [C]//Proceedings of SPIE, 2009, 7210: 721008. [12] Vujkovic-Cvijin P, Goldstein N, Fox M J, et al. Adaptive spectral imager for space-based sensing [C]//Proceedings of SPIE, 2006, 6206: 62060X. [13] Zheng Zhizhong, Yang Zhong, Xiu Liancun. Development and application of shortwave infrared convex blazed grating with high diffraction efficiency [J]. Acta Optica Sinica, 2020, 40(12): 1205002. (in Chinese) doi: 10.3788/AOS202040.1205002 [14] Wang Dinglu, Zhang Daiwei, Xu Banglian, et al. Research progress in design and fabrication of convex grating [J]. Laser & Optoelectronics Progress, 2021, 58(11): 1100002. (in Chinese) [15] Xu D, Owen J D, Papa J C, et al. Design, fabrication, and testing of convex reflective diffraction gratings [J]. Optics Express, 2017, 25(13): 15252-15267. doi: 10.1364/OE.25.015252 [16] Liu Quan, Wu Jianhong, Guo Peiliang, et al. Fabrication of convex blazed grating with high diffraction efficiency [J]. Chinese Journal of Lasers, 2019, 46(3): 0313001. (in Chinese) doi: 10.3788/CJL201946.0313001 [17] Wang Haibin, Liu Quan, Wu Jianhong. Fabrication of convex blazed grating by Ar+ ion-beam etching [J]. Acta Optica Sinica, 2011, 31(4): 0405002. (in Chinese) doi: 10.3788/AOS201131.0405002 [18] Zhao Y, He W, Liu Z, et al. Optical design of an Offner coded aperture snapshot spectral imaging system based on dual-DMDs in the mid-wave infrared band [J]. Optics Express, 2021, 29(24): 39271-39283. doi: 10.1364/OE.444460 [19] Yun G, Crabtree K, Chipman R A. Three-dimensional polarization ray-tracing calculus I: Definition and diattenuation [J]. Applied Optics, 2011, 50(18): 2855-2865. doi: 10.1364/AO.50.002855 [20] Yun G, McClain S C, Chipman R A. Three- dimensional polarization ray-tracing calculus II: Retardance [J]. Applied Optics, 2011, 50(18): 2866-2874. doi: 10.1364/AO.50.002866 [21] Moharam M, Grann E B, Pommet D A, et al. Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings [J]. J Opt Soc Am A, 1995, 12: 1068-1076. doi: 10.1364/JOSAA.12.001068 [22] Moharam M, Pommet D A, Grann E B, et al. Stable implementation of the rigorous coupled-wave analysis for surfacerelief gratings: Enhanced transmittance matrix approach [J]. J Opt Soc Am A, 1995, 12: 1077-1086. doi: 10.1364/JOSAA.12.001077 [23] Wilson D W, Maker P D, Muller R E, et al. Recent advances in blazed grating fabrication by electron-beam lithography [C]//Proceedings of SPIE, 2003, 5173: 51730E. [24] Aoyagi Y, Namba S. Blazed ion-etched holographic gratings [J]. Optica Acta: International Journal of Optics, 1976, 23(9): 701-707. [25] Mouroulis P, Hartley F T, Wilson D W, et al. Blazed grating fabrication through gray-scale X-ray lithography [J]. Optics Express, 2003, 11(3): 270-281. doi: 10.1364/OE.11.000270 [26] Davies M A, Dutterer B S, Suleski T J, et al. Diamond machining of diffraction gratings for imaging spectrometers [J]. Precision Engineering, 2012, 36(2): 334-338. doi: 10.1016/j.precisioneng.2011.09.006