Optimal Design of Wide Angle Diffractive Optical Element
-
摘要: 为进一步研究入射角度的增大对衍射光学元件(diffractive optical element, DOE)衍射效率及微结构高度等参数的影响,分析了入射角度和周期宽度对带宽积分平均衍射效率的影响。基于扩展标量衍射理论,建立了DOE的微结构高度与入射角度和周期宽度的数学模型,提出了工作在一定入射角度范围内,基于复合带宽积分平均衍射效率(comprehensive polychromatic integral diffraction efficiency, CPIDE)最大化实现设计波长和微结构高度等结构参数的优化设计方法。以工作在红外波段的DOE为例进行分析。结果表明:当相对周期宽度为20,入射角度范围为0°~40°时,该DOE的CPIDE为94.15%,微结构高度为1.3396 μm。该设计方法可以实现广角DOE的优化设计。Abstract: The influence of the incident angle on the diffraction efficiency and microstructure height of the diffractive optical element (DOE) was analyzed to further study the influence of the incident angle and period width on the polychromatic integral diffraction efficiency (PIDE). Based on the extended scalar diffraction theory (ESDT), a mathematical model of the relationship among the microstructure height, incident angle, and period width of the DOE was established. An optimal design method for structural parameters, such as the design wavelength and microstructure height, was proposed based on maximizing the comprehensive PIDE (CPIDE) within a certain range of incident angles. A DOE operating within the infrared waveband was considered as an example. The results indicate that when the relative period width is 20 and the incidence angle range is 0° to 40°, the CPIDE of the DOE is 94.15%, and the microstructure height is 1.3396 μm. This design method can realize the optimal design of a wide-angle DOE.
-
0. 引言
当前在核电监测、无人机遥感和军事等领域都经常采用红外成像系统,且都要求有一个大的视场,以获取更多的有效消息。但是大视场会带来成像畸变的问题,成像畸变导致像点与物点定位关系不准确,影响成像效果。想要利用这些具有严重畸变的图像就需要对其进行畸变校正。使用传统光学的办法进行畸变校正会导致以下问题:
1)光学系统尺寸增大,变得复杂,不利于很多使用场景下的小型化要求;
2)由于光线的透过率会随透镜的数量增加而降低,红外波段的光学镜头受限于使用材料,透过率衰减更为突出。
因此,通常都采用图像处理的办法来校正畸变,可以避免上述问题,且校正效果优良。
对相机进行畸变校正可以分为两步,第一步就是获取畸变参数,称之为相机标定;第二步是根据标定的畸变参数对畸变图像还原[1-2]。在进行标定时,通常利用标定板对相机进行标定[3],这种方法过程比较繁琐,且对标定板精度要求较高,而红外标定板由于成像原理导致加工更为困难。因此,更适合采用自标定的方法进行红外相机标定。自标定方法可以通过两幅图像中特征点对的坐标关系,来标定畸变参数[4]。常用的自标定方法又以非线性迭代的方法为主,但是非线性迭代的方法涉及到相机内参、畸变参数初值的选择,且对初值敏感。除此之外还有非迭代的方法求解相机的畸变参数。
文献[5]通过实验仿真结果证明了摄像机参数关联的问题,提出可以通过其他算法来替代非线性迭代从而获得需要的参数。文献[6]提出了九点非迭代优化算法对相机进行标定,但这种方法对图像质量要求高,不能直接适用于红外相机。
本文使用了非迭代的方法对大视场红外相机进行标定。在不知道场景信息和相机参数的前提下,用一种改进算子的SURF算法对畸变图像进行特征点检测与匹配,使用除DM模型和核密度估计算法对相机进行标定,并采用基于边缘保持的插值方法对畸变图像进行校正。
1. 成像原理和畸变定义
1.1 成像原理
相机的成像过程涉及到坐标系的转换,经典成像过程就是在4个坐标系中投影映射的过程,如图 1所示。所遵循的模型为针孔成像模型,物点的物像之间的理想关系是线性关系,也称为线性投影模型。
大视场相机为了使用有限大的成像器件获取更大的视场,使用的成像模型是非线性模型,常见的有等距投影模型、等区域投影模型、正射投影模型、球面投影模型。这些模型投影半径与视场角之间的投影关系如图 2所示,编号1~5的曲线分别是线性投影模型、球面投影、等距投影、等区域投影、正射投影。
1.2 畸变定义
相机的光学畸变是成像过程中的非线性几何失真导致的,这种畸变是由多种类型的畸变组合起来的,通常包括径向畸变(Radial distortion)、切向畸变(Tangential distortion)和薄棱镜畸变(Thin prism distortion)等这几种类型[7]。
径向畸变是光学畸变中对成像效果影响最大的,表现为成像点沿着径向方向有向内或向外的偏移,其中向内偏移为负畸变,向外偏移为正畸变。假设物面为正方形,则正畸变表现为枕形畸变,负畸变表现为桶形畸变,如图 3所示。
根据相机模型,径向畸变可用下式表示为:
$$ \left\{ {\begin{array}{*{20}{c}} {{{\hat x}_{\text{d}}} = x\left[ {1 + {k_1}{r^2} + {k_2}{r^4} + {k_3}{r^6} + \cdots } \right]} \\ {{{\hat y}_{\text{d}}} = y\left[ {1 + {k_1}{r^2} + {k_2}{r^4} + {k_3}{r^6} + \cdots } \right]} \end{array}} \right. $$ (1) 上式表示径向畸变和半径r的偶次方正相关,其中k1决定了图像中心位置畸变较小,k2决定边缘区域畸变量较大。对于鱼眼相机一般需要3阶以上进行描述。
除径向畸变以外,根据相机的内部构造及安装工艺水平可知,不能完全共线的光学系统的镜头光轴产生了偏心畸变;根据Brown-Conrady[8]提出的模型,这种畸变由径向分量和切向分量共同产生,它使成像过程中的实际的成像点在径向和切向都产生微小偏移,将偏心畸变分解到平面坐标系的x、y轴,可得数学表达为:
$$ \left\{ {\begin{array}{*{20}{c}} {{{\hat x}_{\text{d}}} = {p_1}\left( {3{x^2} + {y^2}} \right) + 2{p_2}xy} \\ {{{\hat y}_{\text{d}}} = {p_2}\left( {{x^2} + 3{y^2}} \right) + 2{p_1}xy} \end{array}} \right. $$ (2) 式中:p1和p2是为泰勒展开后的切向畸变系数;(x, y)是理想坐标点。
由于偏心畸变的系数p1、p2远小于泰勒展开的径向畸变的系数k1、k2,因此偏心畸变在数量上要比径向畸变小得多,在进行校正时通常只考虑径向畸变的影响。
2. 参数标定与畸变校正
2.1 径向畸变模型
畸变图像的畸变是多种畸变共同作用的结果,但对相机的成像效果影响最大的是径向畸变,因此进行畸变校正时,仅考虑径向畸变。由于镜头的实际成像模型较为复杂,一般都采用多项式模型来近似。
假设无畸变图像上一点坐标为xu=(xu, yu)T,在畸变图像上对应的点坐标为xd=(xd, yd)T,多项式模型可表示为:
$$ \boldsymbol{x}_{\mathrm{u}}-\boldsymbol{e}=\left(\boldsymbol{x}_{\mathrm{d}}-\boldsymbol{e}\right) L\left(r_{\mathrm{d}}, k\right) $$ (3) 式中:e=(x0, y0)T为畸变中心坐标:
$$\begin{gathered} L\left(r_{\mathrm{d}}, k\right)=1+k_1 r_{\mathrm{d}}^2+k_2 r_{\mathrm{d}}^4+k_3 r_{\mathrm{d}}^6+\cdots \\ r_{\mathrm{d}}=\sqrt{\left(x_{\mathrm{d}}-x_0\right)^2+\left(y_{\mathrm{d}}-y_0\right)^2} \end{gathered} $$ 上述多项式模型(Polynomial model, PM)称为乘法模型,在描述较大的畸变时需要使用高阶PM模型,但是参数过多会导致求解困难。Fitzgibbon提出了除法模型,在参数较少的情况下,DM模型具有更强大的畸变描述的能力[9],且便于求解。DM模型可表示为:
$$ \boldsymbol{x}_{\mathrm{u}}-\boldsymbol{e}=\frac{\left(\boldsymbol{x}_{\mathrm{d}}-\boldsymbol{e}\right)}{L\left(r_{\mathrm{d}}, k\right)}$$ (4) 本文使用单参数的DM模型,即:
$$ \boldsymbol{x}_{\mathrm{u}}-\boldsymbol{e}=\frac{\boldsymbol{x}_{\mathrm{d}}-\boldsymbol{e}}{1+k r_{\mathrm{d}}^2}$$ (5) 2.2 特征点的检测与匹配
本文中的自标定算法利用的是红外图像中的特征点满足的几何约束关系,所以首先要提取出红外图像中的特征点并能将它们正确匹配。因为图像存在较大畸变,所以需要选择一种适用于具有较大畸变的红外图像的方法。
红外图像相比于可见光图像存在分辨率低、对比度低、信噪比低、信息量少的特点,因此针对可见光的特征点检测方法不能直接用于红外图像,需要针对红外图像进行改进。
SURF算法是对SIFT算法的一种改进,提高了运算速度和鲁棒性。经典SURF算法的特点是:①SURF特征是图像的局部特征,它对旋转、尺度缩放、亮度变化保持不变性,对视角变化、仿射变换、噪声也保持一定程度的稳定性;②独特性好;③多量性;④可扩展性强。但是也存在一些不足:①特征点精度不高;②提取的特征点容易存在错误检测,导致特征点信息的冗余,稳定性差[10]。在对红外图像进行特征点提取时,这些问题更为凸显。因此,本文在使用SURF算法进行特征点提取时,使用了一种改进的SURF算子,用来对特征点进行提纯,剔除一些错误的特征点。具体做法是引入局部二维熵来进一步刻画特征点的独特性,通过对比特征点的局部二维熵与阈值的关系来实现特征点的提纯[11]。
局部二维熵的定义为:
$$ H = - \sum\limits_{i = 0}^{255} {\frac{{f\left( {i,j} \right)}}{{{N^2}}}} {\log _2}\frac{{f\left( {i,j} \right)}}{{{N^2}}} $$ (6) 式中:i表示像素的灰度值;j表示邻域灰度均值;f(i, j)表示(i, j)出现的频数;N代表尺度。
改进的SURF特征点检测过程如图 4。
在对检测到的特征点进行提纯后,我们使用双向匹配算法对特征点进行匹配。
2.3 相机的标定
由对极几何关系和基础矩阵我们可以知道,同一场景两幅无畸变图像中的一对成像点xu=(xu, yu)T,xu′=(xu′, yu′)T,设畸变中心e为坐标原点,则两者之间满足:
$$ \boldsymbol{x}_{\mathrm{u}}{ }^{\prime} \boldsymbol{F} \boldsymbol{x}_{\mathrm{u}}=0 $$ (7) 式中:F是基本矩阵[12],将式(5)代入式(7)中可得到:
$$ {\left[ {\frac{{{{x'}_{\text{d}}}}}{{1 + k{{r'}_{\text{d}}^2}}}} \right]^{\text{T}}}F\left[ {\frac{{{x_{\text{d}}}}}{{1 + kr_{\text{d}}^{\text{2}}}}} \right] = 0 $$ (8) 将图像上两点坐标用齐次坐标表示,进行坐标还原后在方程的两边均乘以(1+krd2),不改变等式的恒等性。使用直积(Kronecker积)把基础矩阵F拉直,则式(8)变为如下形式:
$$ \left(x_{\mathrm{d}}^{\prime}, y_{\mathrm{d}}^{\prime}, 1+k r_{\mathrm{d}}^{\prime 2}\right) \otimes\left(x_{\mathrm{d}}, y_{\mathrm{d}}, 1+k r_{\mathrm{d}}^2\right) \vec{F}=0 $$ (9) 式(9)表示了两幅畸变图像中一对成像点之间的关系。当两幅图像中9对成像点堆积起来时,我们就可以得到一个关于基础矩阵拉直后向量f的齐次线性方程组:
$$\boldsymbol{M}\left(x, x^{\prime}, k\right) \boldsymbol{f}=0 $$ (10) 式中:矩阵M称为测量矩阵;x和x′表示9对成像点坐标的堆积;f表示基础矩阵F的拉直向量。为方便书写,以下M(x, x′, k)记为M。
由基础矩阵的性质可以知道,等式(10)的解必须是一个奇异矩阵,于是可以得到如下两个等式:
$$ \operatorname{det}(\boldsymbol{M})=0 $$ (11) $$ \operatorname{det}(\operatorname{Mtx}[\operatorname{Ker}[\boldsymbol{M}]])=0 $$ (12) 式中:Ker[]表示零空间算子;Mtx[]表示将向量转为矩阵的运算。在式(11)中,仅包含有一个未知数k,于是可以得到一个关于k的6阶多项式方程,求解出k。
虽然理论上有9对匹配点就能求得k值,但是由于噪声与误差的影响,求得的k值会非常不稳定,可以通过使用多组匹配点进行重复计算来提高结果的稳定性。
通过实验发现,在进行多组实验时,求得的k值始终包围着真正的根,于是可以采用核密度估计的方法来估计根值。
这种方法就是使用核密度估计去寻找最大峰值分布的概率,最大峰值对应的位置就是我们要找的最优根值。假设对随机变量进行核密度估计:
$$ \hat{f}_h(x)=\sum\limits_{i=1}^n \frac{K\left(x_i-x\right)}{h} $$ (13) 式中:K()表示核函数;h表示带宽。高阶方程求得的根作为随机变量。选择带宽为0.3的高斯核来估计概率密度函数,基于对实根的概率密度估计,很容易得到最优根。
2.4 畸变校正
图像畸变校正的过程就是根据相机的畸变模型,将畸变图像恢复成无畸变图像的过程。这个过程需要将无畸变图像的每个像素坐标映射到畸变图像上,但通常情况下,得到的坐标映射结果不是一个整数,因此需要对坐标像素进行插值运算。
在进行坐标变换时,我们使用逆向映射,通过计算无畸变图像任意像素点经映射后在畸变图像中的坐标位置进行校正,这样可以减少计算量。
最常使用的插值方法包括:临近点插值,双线性插值,双三次插值。临近点插值计算量最小,但是效果不好,会存在明显的灰度不连续性。双线性插值在边缘处过渡较为自然,但是也会导致损失图像高频细节,使轮廓变得模糊。双三次插值效果最好,但是计算量非常大,不利于图像畸变校正的实际应用。因此提出一种基于边缘保持的灰度插值方法。基于边缘保持的灰度插值方法可以通过计算插值区域的局部标准差,判断插值区域属于平滑区域或者边缘区域。在图像平滑区域使用双线性插值,减少计算量的同时也不会导致边缘模糊;在边缘区域使用双三次插值,保证图像细节和边缘。
基于边缘保持的插值算法主要思路是:首先确定阈值T,计算当前像素区域标准差;如果当前区域像素标准差小于T,则表明它在平滑区域,当前像素选择用邻域加权替代;如果当前区域像素标准差大于T,则表明它在边缘区域,属于边缘点,此时选择用最大梯度方向的像素点代替当前像素。实现流程如图 5所示。
3. 实验结果及分析
为了验证理论算法的有效性,本文使用大视场红外相机进行了实验。使用的相机为中波红外相机,探测器响应波段为3.3~4.3 μm,系统输出图像分辨率为1280×1024,图像为16 bit,水平视场角约50°,最大畸变约40%。该系统的光学畸变设计结果如图 6所示,可以清楚地看到畸变率随视场角的增大而快速增大,符合2.1节中径向畸变模型的预测结果。本文的实验算法均在同一台电脑上进行,电脑配置为:CPU Intel i7 9750,内存16 GB,操作系统为Windows 10旗舰版;运行环境:Matlab 2020a。
使用大视场红外相机对室外场景进行拍摄,对同一场景从不同位置或角度拍摄两张图像,并对图像进行简单预处理,便于后续利用,如图 7所示。
使用改进算子的SURF算法对两张红外图像进行特征点检测,并使用双向匹配算法进行特征点匹配,得到的结果如图 8所示。从图 8可以看到,使用改进算子的SURF算法可以自动获得特征匹配点对,且准确率较高。
根据2.3节的方法标定相机的畸变参数,我们默认采用文献中常用的方法,设畸变中心坐标为图像中心,即:
$$ \boldsymbol{e} = {\left( {\frac{{{\text{IMW}} + 1}}{2},\frac{{{\text{IMW}} + 1}}{2}} \right)^{\text{T}}} $$ (14) 将匹配到的坐标点每9对堆叠在一起作为一组,可以得到一个关于畸变参数k的6阶方程,并使用Matlab自带的solve函数求解方程的实根。将多组匹配点对求解出的畸变参数的实根,作为样本集,并使用基于核密度估计的方法求取最优解。在进行这一步时,可以将图像坐标归一化,有利于提高畸变参数的估计精度。基于单参数DM模型的畸变参数k值范围在-2~2之间。
设置带宽为0.3,得到实根的概率密度分布曲线如图 9所示,由图 9可知求得的最优解为k=-0.8725。
对畸变图像进行校正时,使用后向映射的方法,需要用ru来表示rd,即:
$$ {r_{\text{d}}} = \frac{{\sqrt {1 - 4kr_{\text{u}}^2} }}{{2kr_{\text{u}}^{}}} $$ (15) 然后使用基于边缘保持的插值方法对校正后坐标不是整数的像素进行灰度插值运算,最终得到校正后的图像。从图 10可以看出,图像保持了良好的边缘特征,没有出现插值导致的锯齿和模糊现象。而且建筑的边缘基本都变成了直线,这意味着原图像中的畸变得到了校正。
图 11对比了校正前后图像像高的变化。校正前,像高与视场角的关系符合鱼眼投影模型,导致了畸变;校正后,像高与视场角的关系符合线性投影模型,消除了畸变。
为了验证这种方法对于不同畸变程度的图像的有效性,本文使用畸变程度不同的图像进行了实验。图 12中的红外图像畸变率小于10%,用上述方法进行畸变校正后得到图 13,从图像中可以看出,房顶边缘处的弧线基本变为直线,图像边缘处的楼房畸变消失,校正效果较为理想。
本文的方法对畸变程度不同的红外图像都实现了畸变校正,且效果较为理想。以上结果说明本文提出的方法对大视场红外图像的畸变校正具有可行性和普适性,能有效处理各种程度的径向畸变。受限于不同畸变程度图像数据的不足,因此,本文未能通过具体实验确定此畸变校正方法的校正上限,但是通过参考可见光图像畸变校正和像差校正的相关文献[13-14],在不考虑特征点检测与匹配误差影响的情况下,该方法可以做到对视场约180°的鱼眼镜头进行畸变校正。
4. 总结
本文通过输入同一相机拍摄的两幅相同场景的图像,实现了大视场红外相机的畸变自动校正,为大视场红外相机的畸变校正提供了一种可行的方案。该方法不需要人工参与,且不需要预先知道相机本身的参数信息,即可完成畸变校正。但是仍然存在一些不足,在标定畸变参数时,需要的特征点对数量较多,否则会导致求解的最优解陷入局部最优。因此在进行图片选取时,需要两幅图像尽可能包含较多的相同场景,便于进行特征点检测时,两幅图像中能够尽可能地获取到更多的特征点对。除此之外,本文的算法在校正畸变图像时,导致了原有的视场有所减小,这些都需要进行后续的优化和修正。
-
表 1 微结构高度与周期宽度的关系
Table 1 Relationship between microstructure height and period width
Incident angle/(°) Period width /λ 5 10 20 30 ∞ 0 1.3758 1.3734 1.3728 1.3727 1.3726 20 1.3579 1.3463 1.3411 1.3395 1.3365 40 1.2689 1.2511 1.2430 1.2404 1.2355 60 1.1266 1.1076 1.0988 1.0960 1.0907 表 2 基于带宽积分平均衍射效率最大化确定的结构参数
Table 2 Structural parameters determined by maximum PIDE
Parameters Incident angle/° 0 20 40 60 Maximum PIDE/% 94.47 94.47 94.48 94.50 Design wavelength/μm 1.7399 1.7298 1.7220 1.7182 Microstructure height/μm 1.3724 1.3360 1.2350 1.0904 表 3 基于带宽积分平均衍射效率最大化确定的结构参数
Table 3 Structural parameters determined by maximum CPIDE
Parameters Incident angle range/° 0-20 0-40 0-60 Design wavelength/μm 1.73 1.72 1.72 Design angle/° 8.25 16 24 Microstructure height/μm 1.3615 1.3396 1.3142 CPIDE/% 94.45 94.15 92.67 -
[1] Greisukh G I, Ezhov E G, Zakharov O A, Diffractive microstructures of zoom lenses for visible and near-infrared ranges based on novel optical plastics[J]. Journal of Optical Technology, 2022, 89(3): 127-131. DOI: 10.1364/JOT.89.000127
[2] 田晓航, 薛常喜. 小F数红外双波段无热化折衍摄远物镜设计[J]. 光学学报, 2022, 42(14): 181-187. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202214023.htm TIAN Xiaohang, XUE Changxi. Athermalization design of small f-number refractive-diffractive telephoto objective lens in infrared dual-band[J]. Acta Optica Sinica, 2022, 42(14): 181-187. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202214023.htm
[3] 段慧慧, 杨艳芳, 何英, 等. 4π聚焦系统中衍射光学元件对聚焦场多光球结构的影响[J]. 光学学报, 2021, 41(20): 174-179. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202120021.htm DUAN Huihui, YANG Yanfang, HE Ying, et al. Influence of diffractive optical elements on multiple spherical spots in a 4π focusing system[J]. Acta Optica Sinica, 2021, 41(20): 174-179. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202120021.htm
[4] 张博, 崔庆丰, 朴明旭, 等. 双波段多层衍射光学元件的基底材料选择方法研究及其在变焦系统中的应用[J]. 光学学报, 2020, 40(6): 0605001. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202006005.htm ZHANG Bo, CUI Qingfeng, PIAO Mingxu, et al. Substrate material selection method for dual-band multilayer diffractive optical elements and its application in the zoom system[J]. Acta Optica Sinica, 2020, 40(6): 0605001. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202006005.htm
[5] Shima G G, David Fischer, Stefan Sinzinger. Multifocal multi-value phase zone plate for 3D focusing[J]. Applied Optics, 2019, 58(32): 8943-8949. DOI: 10.1364/AO.58.008943
[6] XUE Changxi, CUI Qingfeng. Design of multilayer diffractive optical elements with polychromatic integral diffraction efficiency[J]. Optics Letters, 2010, 35(7): 986-988. DOI: 10.1364/OL.35.000986
[7] YANG Liangliang, LIU Chenglin, LI Shengqiang. Optimal design of depth-scaling error for multilayer diffractive optical elements with oblique incidence[J]. Applied Optics, 2017, 56(15): 4532-4536. DOI: 10.1364/AO.56.004532
[8] GAO Long, To Suet, YANG Hongfang, et al. Effect of assembling errors on the diffraction efficiency for multilayer diffractive optical elements[J]. Applied Optics, 2014, 53: 7341-7347. DOI: 10.1364/AO.53.007341
[9] Laborde V, Loicq J, Habraken S. Modeling infrared behavior of multilayer diffractive optical elements using Fourier optics[J]. Applied Optics, 2021, 60(7): 2037-2045. DOI: 10.1364/AO.414082
[10] YANG Hongfang, XUE Changxi. Sensitivity of diffraction efficiency to period width errors for multilayer diffractive optical elements[J]. Applied Optics, 2018, 57(14): 855-860.
[11] PANG H, YIN S Y, DENG Q L, et al. A novel method for the design of diffractive optical elements based on the Rayleigh–Sommerfeld integral[J]. Optics and Lasers in Engineering, 2015, 70: 38-44. DOI: 10.1016/j.optlaseng.2015.02.007
[12] Noponen E, Turunen J. Binary high-frequency-carrier diffractive optical elements: electromagnetic theory[J]. Journal of the Optical Society of America A, 1994, 11(3): 1097-1109. DOI: 10.1364/JOSAA.11.001097
[13] Swanson G J. Binary Optics Technology: Theoretical Limits on the Diffraction Efficiency of Multilevel Diffractive Optical Elements[R/OL]. M. I. T. Technical Report, 1991, https://api.semanticscholar.org/CorpusID:26445902.
[14] Greisukh G I, Danilov V A, Ezhov E G, et al. Comparison of electromagnetic and scalar methods for evaluation of efficiency of diffractive lenses for wide spectral bandwidth[J]. Optics Communications, 2015, 338: 54-57. DOI: 10.1016/j.optcom.2014.10.037
[15] HUO Furong, WANG Wensheng, XUE Changxi. Limits of scalar diffraction theory for multilayer diffractive optical elements[J]. Optik, 2016, 127: 5688-5694. DOI: 10.1016/j.ijleo.2016.03.062