MPPC low-level-light imaging enhancement algorithm based on sub-window box filtering
-
摘要:
针对多像素光子计数器(MPPC)进行微光成像时,图像受光照不足和噪声影响出现的图像亮度低、对比度差、边缘模糊等问题,提出一种基于子窗口盒式滤波的自适应微光图像处理算法。为了减少算法运行时间的同时突出图像的边缘细节信息,利用子窗口盒式滤波器对图像进行分层得到基础层和细节层;对基础层图像采用自适应阈值直方图均衡化拉伸对比度,细节层图像采用自适应增益控制方式进行增强;根据基础层图像中有效灰度值个数占总灰度的比值自适应确定融合系数,将基础层图像与细节层图像融合得到增强后图像。通过微光实验平台设置3组不同照度的微光环境进行实验仿真,验证了本文算法在保持边缘信息和增强细节方面获得了更好的效果。实验结果表明本文算法在标准差、信息熵、平均梯度等客观评价方面优于改进前算法,提升了微光图像的成像效果。
Abstract:Aiming at the problems of low image brightness, poor contrast, and blurred edges caused by insufficient illumination and noise during low-level-light imaging for multi-pixel photon counter (MPPC), an adaptive low-level-light image processing algorithm based on sub-window box filtering was proposed. To reduce the algorithm running time while highlighting the edge detail information of the image, the sub-window box filter was used to layer the image to obtain the basic layer and detail layer. For the image of basic layer, the adaptive threshold histogram equalization was used to stretch the contrast, and the image of detail layer was enhanced by adaptive gain control method. The fusion coefficient was determined adaptively based on the ratio of the number of effective gray values to the total gray in the image of basic layer, and the image of basic layer was fused with the image of detail layer to obtain the enhanced image. Three sets of low-level-light environments with different illumination levels were set by the low-level-light experimental platform for experimental simulation, which verified that the algorithm obtained better results in maintaining edge information and enhancing details. Experimental results show that the proposed algorithm is superior to the previous algorithm in objective evaluation of standard deviation, information entropy, and average gradient, which improves the imaging effect of low-level-light image.
-
引言
目前对于高精度光学平面的加工与检测是科研领域的研究热点之一,高精度光学平面的加工必须有更高精度的检测技术作为支撑。光学干涉测量方法是检验高精度光学平面的有效手段,在干涉平面测量中,通常使用一个高精度平面(λ/20)作为参考镜,而检验一块高精度的平面需要更高精度的平面作为参考基准[1]。因此,测试精度受参考平面精度的制约,这就提出了平面绝对检测要求。平面绝对检测是通过消除干涉仪器的系统误差和参考平面误差对测量结果的影响,从而得到被测平面面形绝对值的测量。
三平面绝对检测以及它的扩展是目前平面绝对检测的主要手段之一。G. Schulz和J. Schwider[2]等人于1967年和1971年提出并发展了传统的三平面绝对检测,可得到平面上沿径向轴线上的绝对面形分布。1954年P. B. Keenan[3]提出了一种伪剪切干涉计量测试技术,通过两平面在某一固定位置的测量结果和被检平面在2个相互垂直方向移动后的测量结果相减后再累加求出被检平面的绝对面形分布。1984年B. S. Fritzz[4]提出了泽尼克多项式拟合法,利用泽尼克多项式特性将平面的表面面形分解为某些正交基函数,然后采用最小二乘法将这些基函数拟合成被检平面的绝对面形。1992年Chiayu Ai和J. C. Wyant[5-6]提出了奇偶函数法,这种方法将平面面形分解为偶奇、奇偶、偶偶、奇奇函数项,再根据6次测量的结果分别求出这4项的值。2001年K. R. Freischlad[7]提出了旋转剪切的方法,这种方法也是在传统三平面绝对检测的基础上提出的。
由于上述方法在实际应用中各有不同的缺点,因此它们还只是一种实验室方法[8]。基于奇偶函数的六步绝对检测算法[5-6]和仿真测量[9],在仿真测量中,旋转次数越多,其理论精度越高。然而,在实际的高精度测量中,旋转次数越多,引入的旋转角度误差和旋转偏心误差就越大。如$ \text{PV}\ll {1}/{20\lambda }\;$,多次旋转对绝对测量精度和重复性影响非常大。本文依据四步绝对检测算法[10],利用商用菲索型干涉仪搭建实验系统进行测量[11-12]。
1 测量原理
文中讨论的平面绝对检测基于菲索干涉系统,方法为奇偶函数四步旋转测量方法。在奇偶函数法中,函数F(x, y)可以用偶偶函数项Fee、奇偶函数项Foe、偶奇函数项Feo、奇奇函数项Foo:
$$ {F_{\left( {x, y} \right)}} = Fee + Foe + Feo + Foo $$ (1) 函数的奇偶项、偶奇项和偶偶项很容易精确得到,而奇奇项可以用傅立叶正弦级数第一项Foo, 2oddθ来进行近似,通过旋转90°的运算可以算出奇奇项分量中的第一项,其基本原理如图 1所示。
图 1所示的A(x, y)、B(x, y)、C(x, y)为测量过程中使用的3个平面,首先在每个平面上建立直角坐标系xoy,将这3个平面两两进行干涉,得到4组干涉测量的数据,测量步骤如图 2所示,其中平面B与平面C沿X轴翻转得到BX与CX,平面A逆时针旋转90°得到A90。
四步旋转测量方法的函数表达式为
$$ \begin{array}{l} {M_1} = A + {B^X},{M_2} = {A^{90}} + {B^X},{M_3} = A + {C^X},\\ {M_4} = B + {C^X} \end{array} $$ (2) 通过旋转90°的方式联立方程可解出平面A、B、C的奇偶项、偶奇项、偶偶项和奇奇项,因此通过4次组合测量就能够确定3个平面的绝对面形误差分布,即:
$$ \begin{array}{l} A = {A_{ee}} + {A_{oe}} + {A_{eo}} + {A_{{\rm{oo, 2odd \mathit{ θ} }}}}\\ B = {B_{ee}} + {B_{oe}} + {B_{eo}} + {B_{{\rm{oo, 2odd \mathit{ θ} }}}}\\ C = {C_{ee}} + {C_{oe}} + {C_{eo}} + {C_{{\rm{oo, 2odd \mathit{ θ} }}}} \end{array} $$ (3) 根据以上理论推导得出的奇偶函数平面绝对测量结果,须使第二步的90°旋转尽量精确,同时干涉系统的检测光学系统与参考镜及被测平面共光轴。
2 仿真实验与分析
在三面互检法的实验中,3个被测平面需要相互交换位置测量,且有一个被测平面需要旋转一定的角度,这就要求在互换、旋转和翻转后被测平面的空间坐标轴对准。实验中实际旋转角度与理论角度的差异、平面旋转中心与系统光学中心的差异都会影响测量结果的正确性。鉴于此,本文对绝对检测技术能够达到的理论精度和影响检测精度的因素进行了模拟分析,以实现平面面形绝对检测技术应用于实际时的优化配置。
2.1 仿真实验
首先,采用泽尼克多项式前36项构建原始平面A、B、C,把平面A、B、C按图 2所示步骤输入到计算机中进行模拟4次测量。分别得到4次测量结果M1、M2、M3、M4,并在将4次测量结果输入到计算机中得到平面A′、B′、C′,模拟的3个平面与计算出的3个平面的偏差,如图 3所示。其算法偏差的PV值在10-5 nm量级,RMS值在10-6 nm量级,达到高精度测量的理论要求。
2.2 旋转角度误差的影响
在实际测量中,实际旋转角度会与理论旋转角度有一定的差异,即旋转角度误差。在如图 2所示的第二步中,需要将平面A逆时针旋转90°,因此旋转角度误差的不同可能会导致不同的测量结果。使用泽尼克多项式前36项构建3个平面A、B、C,在仿真过程中,软件干涉图直径为400 pixel,平面A的实际旋转角度为理论旋转角度(90°)叠加一个旋转角度误差(1′),并对得到的测量结果两两相减,进行点对点作差处理,模拟结果如图 4所示。
由图 4可以看出,测量误差随着旋转角度误差的增大而增大,当旋转角度误差达到8′(0.13°)时,测量误差PV值为0.000 1 λ,与文献[9]中使用六步旋转测量法仿真实验,旋转角度误差为0.1°时对测量结果的影响相吻合。
2.3 旋转偏心误差的影响
在实际测量中,平面旋转中心与光学系统中心存在横向偏移,即旋转偏心误差。图 5为在X方向等量偏心(1/2 pixel)的情况下对绝对检测方法的测量误差进行模拟分析的结果。
由图 5可以看出,测量误差随着旋转偏心误差的增大而增大,当发生3 pixel的旋转偏心误差时,测量误差PV值为0.005 λ,此时,旋转偏心误差为系统的重要误差源。因此进行绝对检验时,必须进行精确的旋转中心定位。
3 实验研究
3.1 实验仪器与参数设计
在实际测量中,采用苏州慧利仪器生产的HOOL L9600干涉仪进行测量,如图 6所示。HOOL干涉仪是以移相干涉技术为原理的干涉仪,采用的是菲索干涉仪光路,它是一款高精度干涉测量仪器,其系统误差优于λ/500(λ=632.8 nm)。
为了减小环境的扰动,所有实验装置置于防震气浮光学平台上,并且由屏蔽罩笼罩,实验室中应尽可能减小人为等因素产生的扰动。
我们选取150 mm口径的标准平面镜B、C,以及150 mm口径的被测平面A进行测量。并将3个平面的x、y方向预先标定好,以便测量时确定平面的安放位置。
3.2 实验步骤
在本实验中,我们选被测平面A、标准平面镜B和C在HOOL干涉仪上进行测量。由于被测平面要旋转,将被测平面A放置在旋转台上,可以确保A逆时针旋转90°时的旋转角度误差小于2′, 即测量误差PV值小于0.000 05 λ。
对于旋转偏心误差,我们将一个可以随被测平面旋转的十字线装置安装在旋转台上,将被测平面在0°位置时十字线中心与被测平面在90°位置时十字线中心调整至与系统光轴重合,并将旋转中心的坐标输入到计算机中,可将旋转偏心误差控制在2 pixel(0.375 mm/pixel)以内,即测量误差PV值小于0.002 5 λ。
将3个平面与仪器放置在同一环境下进行等温处理后,按测量原理,实验步骤如下:
1) 在参考平面位置装上标准平面镜B,被测平面位置装上被测平面A;
2) 校准被测平面A的旋转中心;
3) 调整完毕后,测量并记录干涉结果M1;
4) 将被测平面A逆时针旋转90°,调整完毕后,测量并记录干涉结果M2;
5) 将标准平面镜B取出,换上标准平面镜C,并将被测平面A旋转至步骤1位置。调整完毕后,测量并记录干涉结果M3;
6) 将被测平面A取出,换上标准平面镜B,调整完毕后,测量并记录干涉结果M4;
7) 将4次测量得到的干涉结果M1、M2、M3、M4输入计算机进行计算,得到平面A、B、C的绝对形貌。
3.3 实验结果及分析
PV值通常被用于描述元件或系统的光学质量,描述一个波前的PV值可以理解为该波前的最高点到最低点的间距。PV值容易受到粗大误差的干扰,为了保证PV值测量的重复性,HOOL干涉仪以PV10代替PV,PV10即以波面上最高的10个点的均值减去最低的10个点的均值。PV、PV10分别用如下公式表示:
$$ \begin{array}{l} {\rm{PV}} = {W_{\max }} - {W_{\min }}\\ {\rm{P}}{{\rm{V}}_{10}} = \frac{{\sum\nolimits_k^{10} { = 1\omega \max , k} }}{{10}} - \frac{{\sum\nolimits_k^{10} { = 1\omega \min , k} }}{{10}} \end{array} $$ (4) 因此,在本文中,使用PV10和RMS这2个数字波面指标来评价实验数据的可重复性和实验检测精度。
按上述步骤进行测量,可得到各位置的测量结果M1、M2、M3、M4的面形分布, 并由此计算出被测平面A的绝对面形;各面形分布以及面形的PV10值、RMS值如图 7所示。
从图 7的这组结果可以看出,被测平面A的PV10值为0.041 5 λ,小于M1的0.054 λ和M2的0.044 2 λ;被测平面A的RMS值为0.008 7λ,小于M1的0.011 4 λ和M2的0.010 2λ。在本实验中对比说明这些数据,通过绝对检测的方法可以减少干涉仪系统误差的影响,提高被检面的检测精度。
由于空气扰动、环境震动的影响,平面面形绝对检测结果会在一定范围内浮动,而使用绝对检测的意义在于被检平面与参考平面精度相当,可剔除参考面面形的影响。所以在高精度检测中,检测结果的重复性是检测结果精度的前提。为了考察绝对检测的重复性,表 1记录了这次检测中10组数据的中间过程和最终结果的PV10值和RMS值,用标准差σ来表证检测的重复性:
$$ \sigma = \sqrt {\left( {1/N} \right){{\sum\nolimits_{{\rm{i}} = 1}^N {\left( {{{\rm{X}}_{\rm{i}}} - \mu } \right)} }^2}} $$ (5) 表 1 10次测量结果及标准差Table 1. Results of 10 measurements and standard deviation名目 M1 M2 M3 M4 A RMS/λ PV10/λ RMS/λ PV10/λ RMS/λ PV10/λ RMS/λ PV10/λ RMS/λ PV10/λ 1 0.011 4 0.054 0 0.010 2 0.044 2 0.012 3 0.054 3 0.012 4 0.049 7 0.008 7 0.041 5 2 0.011 3 0.053 8 0.009 8 0.043 2 0.012 4 0.054 8 0.012 6 0.050 0 0.008 8 0.041 1 3 0.011 4 0.054 0 0.010 0 0.043 7 0.012 7 0.054 5 0.012 6 0.050 5 0.008 6 0.041 2 4 0.011 3 0.053 8 0.010 0 0.044 3 0.012 3 0.053 4 0.012 5 0.050 5 0.008 7 0.042 4 5 0.011 4 0.054 0 0.010 1 0.044 5 0.012 3 0.053 2 0.012 6 0.050 4 0.008 8 0.042 0 6 0.011 4 0.053 6 0.010 1 0.045 0 0.012 3 0.053 2 0.012 6 0.050 8 0.008 8 0.042 1 7 0.011 3 0.054 6 0.010 0 0.044 2 0.012 6 0.054 7 0.012 6 0.051 1 0.008 8 0.041 4 8 0.011 6 0.054 2 0.009 9 0.043 3 0.012 4 0.053 3 0.012 5 0.050 9 0.008 7 0.040 8 9 0.011 5 0.054 5 0.010 0 0.044 7 0.012 4 0.053 8 0.012 5 0.050 0 0.008 8 0.041 8 10 0.011 5 0.054 8 0.010 0 0.044 7 0.012 4 0.053 3 0.012 5 0.050 1 0.008 9 0.041 9 平均值 0.011 4 0.054 1 0.010 0 0.044 1 0.012 4 0.053 8 0.012 5 0.050 4 0.008 7 0.041 6 标准偏差 0.000 1 0.000 3 0.000 1 0.000 6 0.000 1 0.000 6 0.000 1 0.000 4 0.000 1 0.000 5 从表 1中可以看出:1)绝对检测得到的结果A比传统的检测方法得到的结果(M1、M2)精度要高,尤其是RMS值精度;2)实验结果重复性较好,绝对检测得到结果的PV10值重复性在±λ/1 000范围内,RMS值重复性在±λ/7 000范围内。
上述的重复性仅仅是讨论的测量结果数字指标(即PV10、RMS)的重复性,为了讨论各个位置的整体面形分布的重复性,将同一被测平面两次独立的绝对检测结果进行点对点作差处理,用残留误差的大小反映整体面形分布的重复性。如图 8所示,其残差图PV10值为0.004 λ,RMS值为0.000 5 λ,所以对于各个位置,整体面形分布的重复性较好。
4 结论
本文根据现代高精度光学平面检测工业测量精度高、测量稳定可靠等要求,尤其是部分光学平面的表面面形精度远远超出现有商用干涉仪的面形检测水平,介绍了一种基于奇偶函数的绝对检测方法,分别就旋转角度误差及旋转偏心误差对测量误差的影响进行了仿真分析,仿真结果显示绝对检测技术对旋转偏心误差敏感。因此,在进行绝对检测时,应进行严格的旋转中心定位。并基于商用干涉仪搭建了实验系统进行验证,测得实际样品的绝对检测精度是PV10值为0.041 5λ,RMS值为0.008 7λ,小于常规干涉检测所得结果。实验结果表明:使用本文方法可以减少干涉仪器的系统误差和参考平面误差的影响,提高被检平面的检测精度。对测量重复性进行研究,实验结果显示其PV10值重复性在±λ/1 000范围内,RMS值重复性在±λ/7 000范围内;并对同一平面两次独立的绝对检测结果进行点对点作差处理,从而获得残差图,其残差图PV10值为0.004λ,RMS值为0.000 5λ,表明测量结果的重复性高,实现了高精度平面面形检测,对光学平面加工与检测具有一定的指导意义。
-
表 1 ${10^{ - 1}}\;{\rm{lx}}$微光环境下的图像增强评价指标
Table 1 Image enhancement evaluation indexes under ${10^{ - 1}}\;{\rm{lx}}$ low-level-light environment
算法 标准差 信息熵 平均梯度 A B C A B C A B C 原始图像 34.63 15.87 29.24 6.38 1.20 4.59 1.26 0.57 0.83 交叉子窗
口算法35.51 14.43 33.75 6.41 1.22 4.98 1.37 0.61 0.96 DDE算法 35.82 16.11 33.56 6.46 1.16 4.81 1.78 0.72 1.06 文献[4]算法 37.35 17.15 36.57 6.47 1.19 4.94 1.43 0.89 1.12 本文算法 39.58 20.56 40.36 6.49 0.77 5.19 2.27 1.15 1.46 表 2 ${10^{ - 3}}\;{\rm{lx}}$微光环境下的图像增强评价指标
Table 2 Image enhancement evaluation indexes under ${10^{ - 3}}\;{\rm{lx}}$ low-level-light environment
算法 标准差 信息熵 平均梯度 A B C A B C A B C 原始图像 21.29 16.71 19.42 5.85 1.81 6.04 0.77 0.69 0.64 交叉子窗
口算法22.39 17.24 19.57 6.04 1.73 6.06 0.87 0.72 0.68 DDE算法 24.67 18.08 19.65 5.87 1.66 5.98 1.06 0.80 0.71 文献[4]算法 26.89 18.70 19.82 6.01 1.63 6.01 0.98 0.87 0.78 本文算法 28.94 19.30 28.43 6.15 1.39 6.39 1.18 0.95 1.13 表 3 ${10^{ - 5}}\;{\rm{lx}}$微光环境下的图像增强评价指标
Table 3 Image enhancement evaluation indexes under ${10^{ - 5}}\;{\rm{lx}}$ low-level-light environment
算法 标准差 信息熵 平均梯度 A B C A B C A B C 原始图像 15.53 9.17 8.08 5.27 1.55 4.55 0.45 0.37 0.26 交叉子窗
口算法14.88 10.10 14.07 5.28 1.38 4.61 0.42 0.43 0.41 DDE算法 17.15 10.57 11.59 5.24 1.30 4.58 0.58 0.44 0.37 文献[4]算法 17.46 11.23 12.02 5.46 1.27 4.61 0.60 0.48 0.42 本文算法 16.57 12.48 15.16 5.21 1.32 3.89 0.63 0.57 0.51 表 4 4种算法运行时间对比
Table 4 Comparison of running time of four algorithms
s 算法 微光环境 $ {10^{ - 1}}\;{\rm{lx}} $ $ {10^{ - 3}}\;{\rm{lx}} $ $ {10^{ - 5}}\;{\rm{lx}} $ 交叉子窗口算法 10.371 9.114 6.947 文献[4]算法 9.210 7.062 3.763 本文算法 2.816 2.177 1.041 -
[1] 徐炜君, 刘国忠. 空间域和频域结合的图像增强技术及实现[J]. 中国测试,2009,35(4):52-54. XU Weijun, LIU Guozhong. Image enhancement technology combining spatial domain and frequency domain and its implementation[J]. China Measurement & Test,2009,35(4):52-54.
[2] 周克虎, 雷涛, 罗刚. 一种基于时域滤波的红外序列图像去噪算法[J]. 应用光学,2021,42(3):474-480. ZHOU Kehu, LEI Tao, LUO Gang. An infrared sequence image denoising algorithm based on time-domain filtering[J]. Journal of Applied Optics,2021,42(3):474-480.
[3] LIU Ning, ZHAO Dongxue. Detail enhancement for high-dynamic-range infrared images based on guided image filter[J]. Infrared Physics & Technology,2014,67:138-147.
[4] ZHOU B, LUO Y, YANG M, et al. An improved adaptive detail enhancement algorithm for infrared images based on guided image filter[J]. Journal of Modern Optics,2019,66(1):33-46. doi: 10.1080/09500340.2018.1511861
[5] 王炫, 尹丽菊, 高明亮, 等. 基于新符号函数与盲源分离的光子计数图像去噪方法[J]. 激光与光电子学进展. 2018, 55(10): 158-164. WANG Xuan, YIN Liju, GAO Mingliang, et al. Photon counting image denoising method based on new symbol function and blind source separation[J]. Laser & Optoelectronics Progress. 2018, 55(10): 158-164.
[6] 廖斌, 刘鸳鸳. 基于多尺度灰度变换的图像增强研究[J]. 量子电子学报,2015,32(5):550-554. LIAO Bin, LIU Yuanyuan. Research on image enhancement based on multi-scale gray transform[J]. Chinese Journal of Quantum Electronics,2015,32(5):550-554.
[7] YIN Hui, GONG Yuanhao, Qiu Guoping. Combined window filtering and its applications[J]. Multidimensional Systems and Signal Processing,2021,32(2):313-333.
[8] KITAGAWA G. The two-filter formula for smoothing and an implementation of the Gaussian-sum smoother[J]. Annals of the Institute of Statistical Mathematics,1994,46:605-623. doi: 10.1007/BF00773470
[9] 沈德海, 张龙昌, 鄂旭. 基于多子窗口的混合噪声滤波算法[J]. 计算机技术与发展,2015,25(6):69-72. SHEN Dehai, ZHANG Longchang, E Xu. Hybrid noise filtering algorithm based on multiple sub-windows[J]. Computer Technology and Development,2015,25(6):69-72.
[10] CÁTIA P, ANA O, CRISTINA J, et al. Automatic crackle detection algorithm based on fractal dimension and box filtering[J]. Procedia Computer Science,2015,64:705-712. doi: 10.1016/j.procs.2015.08.592
[11] 胡海波, 邱民仆. 微光图像对比度增强技术研究[J]. 测绘通报,2015,1:26-30. HU Haibo, QIU Minpu. Contrast enhancement of low light level image[J]. Bulletin of Surveying and Mapping,2015,1:26-30.
[12] 孙国栋, 徐昀, 徐亮. 不均匀光照和重影的仪表图像二值化方法[J]. 应用光学,2020,41(1):74-78. doi: 10.5768/JAO202041.0102002 SUN Guodong, XU Yun, XU Liang. Instrumentation image binarization method for inhomogeneous illumination and ghosting[J]. Journal of Applied Optics,2020,41(1):74-78. doi: 10.5768/JAO202041.0102002
[13] 汪伟, 许德海, 任明艺. 一种改进的红外图像自适应增强方法[J]. 红外与激光工程,2021,50(11):419-427. WANG Wei, XU Dehai, REN Mingyi. An improved adaptive enhancement method for infrared image[J]. Infrared and Laser Engineering,2021,50(11):419-427.
[14] 李佳, 李少娟, 段小虎, 等. 基于Retinex理论与概率非局部均值的红外图像增强方法[J]. 光子学报,2020,49(4):187-196. LI jia, LI Shaojuan, DUAN Xiaohu, et al. Infrared image enhancement method based on Retinex theory and probabilistic non-local mean[J]. Acta Photonica Sinica,2020,49(4):187-196.
[15] 汪子君, 罗渊贻, 蒋尚志, 等. 基于引导滤波的自适应红外图像增强改进算法[J]. 光谱学与光谱分析,2020,40(11):3463-3467. WANG Zijun, LUO Yuanyi, JIANG Shangzhi, el at. An improved adaptive infrared image enhancement algorithm based on guided filtering[J]. Spectroscopy and Spectral Analysis,2020,40(11):3463-3467.
-
期刊类型引用(5)
1. 周思羽,杨文奇,卢建华,宋立廷. 基于LESO的舰载机纵向着舰动态面抗饱和控制技术. 测控技术. 2023(03): 124-133 . 百度学术
2. 王国龙,白端元,王愚,刘云清. SOA在音圈电机控制系统中应用仿真研究. 计算机仿真. 2023(04): 277-281 . 百度学术
3. 谢娜,杜言鲁,姜世洲,巩全成,王明超,徐飞飞,王新伟,齐媛. 快速反射镜理论建模与仿真分析. 应用光学. 2023(03): 513-522+585 . 本站查看
4. 王惠林,宁飞,刘吉龙. 航空光电侦察技术发展研究. 应用光学. 2022(05): 825-832 . 本站查看
5. 姜世洲,杜言鲁,王惠林,王明超,巩全成,徐飞飞. 柔性铰链结构特性对快速反射镜性能影响分析. 应用光学. 2022(04): 626-634 . 本站查看
其他类型引用(4)