Influence of Brightness Gain on the Object-Background Contrast of an Image Intensifier
-
摘要: 亮度增益为微光像增强器的重要参数,影响微光夜视环境下以微光像增强器为核心的目标背景对比度测试结果精度。本文基于对比度测试原理,通过对不同亮度增益的ICCD进行目标背景对比度测试,分析得到亮度增益与目标背景对比度的变化关系及影响因素。测试结果表明:亮度增益分别为12000 cd/(m2·lx)、8300 cd/(m2·lx)、6000 cd/(m2·lx)时,环境照度在[2×10-3, 8×10-3](lx)区间内,目标背景对比度受亮度增益影响较小。照度在区间外时,受亮度增益影响,目标背景对比度测试结果可能失真。本文研究结果为提高ICCD目标背景对比度测试精度提供数据支持。Abstract: Brightness gain is an important parameter of low-light image intensifiers. It affects the accuracy of the object-background contrast in low-light night vision environments. Based on the principle of contrast tests, this study analyzes the relationship between the brightness gain and the object-background contrast, as well as influencing factors, via the object-background contrast test on an ICCD with different brightness gains. The results show that when the brightness gain is 12000, 8300, or 6000 cd/(m2·lx) and the ambient illuminance is within the range of [2×10-3, 8×10-3] (lx), the object-background contrast is minimally affected by the brightness gain. However, when the illuminance is outside this range, the object-background contrast test result may be distorted due to the influence of the brightness gain. The results of this study can guide the improvement of accuracy of ICCD object-background contrast tests.
-
Keywords:
- image intensifier /
- ICCD /
- object-background contrast /
- brightness gain
-
0. 引言
目前针对强光照射条件下目标背景对比度的研究较为广泛。然而,已有研究结果表明,随着环境照度的降低,若光照的色温条件发生变化[1],则其目标背景对比度也会发生变化,因此强光照射条件下的目标背景对比度不再适用,或者说会产生较大偏差,而关于夜间目标背景对比度对视距的影响研究却鲜有报道。对于微光夜视系统而言,其典型应用环境多为星光条件,即10-3 lx环境。像增强器在微光领域对图像信息倍增作用被广泛使用,所以采取用微光像增强器和CCD通过光纤光锥或透镜进行耦合[2],利用像增强器的微光下的增强性能,使得CCD在低照度环境下工作,并且可以输出稳定图像。
本文利用微光像增强器的高亮度增益放大特性,对星光条件下的目标背景对比度进行测量及验证。针对像增强器亮度增益特性表现出在一定范围内输入照度和输出亮度为线性变化的现象,对线性变化范围内不同亮度增益下的照度-灰度关系以及照度-目标背景对比度关系进行研究。由此对亮度增益在微光环境下以微光像增强器为核心的目标背景对比度的影响展开讨论。满足微光夜视整机视距评估时对夜间星光环境下基于国产三代ICCD的目标背景对比度特性的准确测量需求。
1. 测试原理和方法
1.1 测试原理
由微光环境下物体辐射光传递过程,利用目标对比度理论公式,以国产ICCD作为测试器件,建立ICCD目标对比度测试数学模型,图 1所示为微光环境对比度测试示意图[3-4]。
利用灰度反应目标与背景调制对比度的公式为[5]:
$$ C{\text{ = }}\left| {\frac{{{L_{\text{M}}} - {L_{\text{B}}}}}{{{L_{\text{M}}}{\text{ + }}{L_{\text{B}}}}}} \right|{\text{ = }}\left| {\frac{{{g_{\text{M}}} - {g_{\text{B}}}}}{{{g_{\text{M}}}{\text{ + }}{g_{\text{B}}}}}} \right| $$ (1) 式中:LM为探测环境的目标亮度值;LB为探测目标的亮度值;gB为探测环境的背景灰度值;gM为探测目标的灰度值;C为调制对比度。
测试系统中的光信息传递分为3个主要过程:一是光源光束在目标处反射并入射到探测仪器中;二是在ICCD中通过像增强器光放大;三是通过耦合的CCD器件将光信号转化为灰度信息。
第一个传输过程主要影响因素为在光传播过程中大气的衰减损耗以及物镜的衰减损耗τobj,光源在物体表面发生反射,反射光成像到像增强器输入面上。光谱反射率ρ表征物体的颜色信息,表现了颜色材质特征,是物体的固有属性,目标探测时只与物体本身的相关,无关于环境照度[6]。目标和背景的光谱反射率不同形成差异。经过物镜成像,反射光信息成像到像增强器输入面上的辐射照度为[3]:
$$ {E_\rm K} = \frac{\pi }{4}{(\frac{D}{f})^2}{\tau _{\rm obj}}{L_\rm A} $$ (2) 式中:EK为目标反射光到达物镜表面的辐射照度;D为物镜孔径大小;f为物镜的焦距;LA为目标反射光到达物镜表面的辐射亮度。
第二个传输过程中,光信息在ICCD中通过像增强器光放大。如图 2所示,三代微光像增强器ABC曲线在最大输出亮度Em(点D所对应位置)和最小输出亮度ES(图 2点A所对应位置)区间内表现出近似线性变化,即亮度增益为定值。C点(B点)处于线性区域内,由于接近上限点D(下限点A),输出信号出现上限饱和(下限截断)的现象[3]。当输入照度大于最大输出照度或小于最小输出照度时呈现出非线性变化。线性区间内亮度增益表达式为:
$$ G = \frac{{{L_\rm A}}}{{{E_\rm K}}} $$ (3) 第三个传输过程中,通过耦合的CCD器件将光信号转化为灰度信息。像增强器输出辐射亮度为光锥输入辐射亮度。考虑光纤光锥耦合对光信息带来的损耗以及CCD工作时的光-电-光变化过程,光纤光锥的有效透过率T为下式[7-9]:
$$ T = \frac{{{\mathit{\Phi} _{\rm out}}}}{{{\mathit{\Phi} _{\rm in}}}} = {\left| {\frac{{{D_{\rm small}}}}{{{D_{\rm large}}}}} \right|^2}{T_\rm R}{T_\rm A}{K_\rm C} $$ (4) 式中:Φout光纤光锥小端面输出的光通量;Φin为光纤光锥大端面的输入辐射通量;Dsmall光纤光锥小端面直径;Dlarge光纤光锥大端面直径;TA光纤光锥材料的透过率;TR光锥光纤反射损耗后的透过率;KC为填充效率。
CCD接收光纤光锥输出光信息,将光信息转换为灰度图像,整个输出图像的灰度分布,是由数个单个像素共同组成的。每个像素接收到的辐射照度Ee为:
$$ {E}_{\rm e}=\frac{{\mathit{\Phi} }_{\rm out}s{\tau }_{粘合剂}}{{s}_{光敏面}} $$ (5) 式中:s光敏面为光敏面面积;τ粘合剂为粘合剂损耗;s为像素面积。
根据电荷转移原理,可以得到单个像素的电荷量Q为式(6):
$$ Q = \frac{{\eta es}}{{h\nu }}\int {{E_\rm e}{\rm d}t} $$ (6) 式中:η为光电转换效率;e为电子电荷量(1e=1.602×10-19 C);h为普朗克常数(h=6.62607015×10-34 J·s);ν为辐射电磁波的频率。
综合前文讨论光学传递链的影响因素得到输出图像上对应的灰度值为:
$$ g=Q+o+σ $$ (7) 式中:o为输出过程中的直流偏置;σ为输出过程中存在的随机噪声[10-11]。
由于CCD的像素不均匀性,每一个像素的灰度响应、光电转换效率和像素面积不尽相同。CCD上任意位置(i, j)处灰度值可表示为:
$$ g\left[ {i,j} \right] = \frac{{\eta \left[ {i,j} \right]e{s^2}\left[ {i,j} \right]}}{{h\nu }}{E_\rm e}t + o + \sigma $$ (8) 当选取目标位置为(m, n),背景位置为(u, v)时,对比度可表示为下式:
$$ C= \left\{\begin{array}{l}\left|\frac{\frac{\eta \left[m,n\right]{s}^{2}[m,n]}{{\nu }_{\rm m}}{\rho }_{\rm M}^{}-\frac{\eta \left[u,v\right]{s}^{2}[u,v]}{{\nu }_{\rm m}}{\rho }_{ B}^{}}{\frac{\eta \left[m,n\right]{s}^{2}[m,n]}{{\nu }_{\rm M}}{\rho }_{\rm M}+\frac{\eta \left[u,v\right]{s}^{2}[u,v]}{{\nu }_{\rm B}}{\rho }_{\rm B}+2o+2\sigma }\right|\quad {E}_{\rm s}<E<{E}_{\rm m}\\ \left|\frac{\frac{\eta \left[m,n\right]{s}^{2}[m,n]}{{\nu }_{\rm M}}G({E}_{\rm M}){\rho }_{\rm M}-\frac{\eta [u,v]{s}^{2}[u,v]}{{\nu }_{\rm B}}G({E}_{\rm B}){\rho }_{\rm B}}{\frac{\eta \left[m,n\right]{s}^{2}[m,n]}{{\nu }_{\rm M}}{\rm G}({E}_{\rm M}){\rho }_{\rm M}+\frac{\eta [u,v]{s}^{2}[u,v]}{{\nu }_{\rm B}}G({E}_{\rm B}){\rho }_{\rm B}+2o+2\sigma }\right|\quad E<{E}_{\rm s}或E>{E}_{\rm m}\end{array}\right. $$ (9) 式中:ρM为目标的光谱反射率;ρB为背景的光谱反射率;νM为目标辐射电磁波频率;νB为背景辐射电磁波频率;EM为目标的输入照度;EB为背景的输入照度;亮度增益G在非线性变化区间表现为以输入照度为自变量的函数。
式(9)中可以看出,当照度在线性变化区间(ES<E<Em)内时,对比度变化不受亮度增益影响,当照度大于像增强器的最大输入照度或低于最小输入照度时,对比度变化与亮度增益有关。
1.2 测试方法
目标对比度测试系统中,主要采用高信噪比、高分辨率的CCD相机、国产三代微光像增强器、光纤光锥、照度计、经过标定的靶板以及计算机等。光源采用2856 K卤钨灯来模拟夜天光环境;选用照度范围可覆盖1×10-1 lx~1×10-4 lx范围的宽量程微弱光照度计进行测试;靶标1#、2#、3#、4#、5#采用对比度依次为100%、90%、75%、60%以及45%;ICCD器件选用国产三代像增强器以及光纤光锥作为连接CCD与像增强器的中继元件;日本KOWA公司生产的物镜,焦距值f为50 mm,F数为1.4。
在实验室中,该系统选用的三代像增强器亮度增益设置分别为12000 cd/(m2·lx)、8300 cd/(m2·lx)、6000 cd/(m2·lx),在暗室环境进行目标背景对比度测试,对比度测试系统的示意图如图 3所示。
2. 实验结果及分析
由式(1)所示,实际测试中图像的对比度由目标灰度与背景灰度计算所得,通过光学信息传递所接收得到灰度图像。以亮度增益为12000 cd/(m2·lx)条件为例,可以看出对于相同靶标在不同灰度时,靶标上的黑白条纹的灰度值均表现为随着照度的增加而递增,图 4为照度和灰度对应图线。从图 4可以看出,当照度大于一定值时,照度-目标灰度呈现非线性变化关系,根据这种现象对G=12000 cd/(m2·lx)进行线性区间拟合。
其中,线性拟合的程度由决定系数R2决定,系数越接近1,则线性度越好。如图 5(a)所示为G=12000 cd/(m2·lx)的线性拟合结果,当照度区间为[0.002, 0.008](单位:lx)变化区间时,线性拟合的决定系数R2大于0.9273,具有良好的拟合程度,因此认为线性关系可信。以同样方法对G=8300 cd/(m2·lx)进行线性拟合,拟合结果如图 5(b)所示,当照度区间为[0.002, 0.008](lx)变化区间时,决定系数大于0.9306;当G=6000 cd/(m2·lx)进行线性拟合,拟合结果如图 5(c)所示,当照度区间为[1.8×10-3, 1×10-2](lx)变化区间时,决定系数大于0.9244。
像增强器与CCD作为ICCD主要组成部分,非线性变化出现原因有主要由两方面:像增强器的输出亮度非线性区间以及CCD的动态范围。其一,对于像增强器输入照度达到最大输入照度值后,空间上出现电荷量饱和,并且由于荧光屏发光能力的有限性和光电阴极疲劳的原因造成像增强器的输入照度-输出亮度之间的非线性变化[12],此时,由光电阴极输入照度和荧光屏输出亮度表征的亮度增益公式无法适用[13]。其二,像素大小以及噪声影响决定了CCD的动态范围,表现在输出图像灰度和照度之间的非线性变化关系。在对比度计算中,灰度的非线性变化影响对比度变化。
建立目标背景对比度数学模型如图 6(a)所示,平面投影如图 6(b)所示,可以看出,在照度-灰度所对应的线性区间内,对比度的实际变化量与理想变化量接近。当照度值过低时,对比度出现阶跃变化;当照度值过高时,对比度曲线偏离理想变化曲线。因此对于实际测试中,输入照度过低和过高时,目标对比度均不能正确表征图像的对比度信息。
图 7中(a)(b)(c)分别为亮度增益为12000 cd/(m2·lx)、8300 cd/(m2·lx)、6000 cd/(m2·lx)的情况下照度和对比度的对应关系。可以看出曲线均变化平缓,照度和对比度未表现出明显相关性。
根据图 7中不同照度和对比度的变化趋势,在线性区间内对不同增益条件下的目标背景对比度进行比较,亮度增益分别为12000 cd/(m2·lx)、8300 cd/(m2·lx)x、6000 cd/(m2·lx)时,得到图 8所示的在3组不同亮度增益实验样本在输入照度[2×10-3,8×10-3](lx)条件下的Pearson相关系数图。Pearson相关系数值在-1~1之间,当数值越接近1时,说明各组实验结果表现出越强的正相关性,当数值越接近0时,说明各组实验结果越独立。其中,相关系数在0.1~0.4表现为弱相关,0.40~0.69表现为中度相关,0.70~0.89为高度相关。
本实验结果说明不同亮度增益条件下对比度表现出中度以上的正相关性,呈现较为良好的相似性。结合式(9)可以看出,在当前的试验条件下,在照度为[2×10-3,8×10-3](lx)区间内,对比度大小变化与亮度增益未表现出明显相关性;当照度小于2×10-3 lx时,亮度增益为6000 cd/(m2·lx)条件下采集图像灰度值过低,目标背景对比度数值出现阶跃变化,不能正确表征图像目标与背景的差异程度,亮度增益为6000 cd/(m2·lx)与12000 cd/(m2·lx)、8300 cd/(m2·lx)两组试验结果的Pearson相关系数表现出低于-0.70负相关性。当照度大于8×10-3 lx时,在当前试验条件限制下,可得到8300 cd/(m2·lx)、6000 cd/(m2·lx)在[8×10-3,3×10-2](lx)区间内Pearson相关系数为-0.37,表现出弱负相关性。线性区间外,对比度表现出与输入照度的相关性,亮度增益数值影响目标背景对比度数值。
3. 结论
本文从目标对比度测试原理出发,讨论亮度增益对以微光像增强器为核心的目标背景对比度测试的影响,搭建测试平台,进行暗室实验,对获得实验数据进行整理分析。通过实验得到不同亮度增益下照度-灰度、照度-目标背景对比度的变化特性以及不同亮度增益条件下目标背景对比度之间的相关性。分析结果表明:1)对比度大小受光谱反射率、CCD器件量子效率、像素面积以及辐射电磁波频率等影响综合决定;2)基于当前实验条件,亮度增益分别为12000 cd/(m2·lx)、8300 cd/(m2·lx)、6000 cd/(m2·lx)的测试表现为,在照度为[2×10-3, 8×10-3](lx)区间内对比度变化与亮度增益未表现出明显相关性;3)当照度过大和过小时,受像增强器亮度增益的非线性变化以及CCD的动态范围和噪声影响,测试结果失真。该结果为国产ICCD目标背景对比度测试的准确测量提供支撑。
-
-
[1] 刘必鎏, 时家明, 赵大鹏, 等. 目标与背景的对比度分析及其应用[J]. 航天电子对抗, 2007, 24(3): 48-51. https://www.cnki.com.cn/Article/CJFDTOTAL-HTDZ200803019.htm LIU Biliu, SHI Jiaming, ZHAO Dapeng, et al. Analysis and application of object and background[J]. Aerospace Electronic Warfare, 2007, 24(3): 48-51. https://www.cnki.com.cn/Article/CJFDTOTAL-HTDZ200803019.htm
[2] 韩坤, 姚泽, 乔凯, 等. 微光ICCD相机的动态调制传递函数理论模型研究[J]. 红外技术, 2020, 42(3): 294-299. http://hwjs.nvir.cn/article/id/hwjs202003013 HAN Kun, YAO Ze, QIAO Kai, et al. Theoretical model of dynamic MTF of low-light-level ICCD[J]. Infrared Technology, 2020, 42(3): 294-299. http://hwjs.nvir.cn/article/id/hwjs202003013
[3] 向世明, 倪国强. 光电子成像器件原理[M]. 北京: 国防工业出版社, 1999. XIANG Shiming, NI Guoqiang. The Principe of Photoelectonic Imaging Devices[M]. Beijing: National Defense Industry Press, 1999.
[4] 陈旭浪, 党小刚, 郭欣达, 等. 基于三代微光ICCD成像装置的目标对比度影响因素测试分析[J]. 应用光学, 2020, 41(6): 1241-1246. https://www.cnki.com.cn/Article/CJFDTOTAL-YYGX202006017.htm CHEN Xulang, DANG Xiaogang, GUO Xinda, et al. Test and analysis of factors influencing target contrast ratio based on 3rd generation LLL ICCD imaging device[J]. Journal of Applied Optics, 2020, 41(6): 1241-1246. https://www.cnki.com.cn/Article/CJFDTOTAL-YYGX202006017.htm
[5] 田园. CCD图像对比度测试方法研究[D]. 长春: 长春理工大学, 2008. TIAN Yuan. Researching on the Testing Method of CCD Picture Contrast[D]. Changchun: Changchun University of Science and Technology.
[6] 周建勋, 王利平, 徐志良. 像增强型夜视系统的图像信息传递[J]. 红外与毫米波学报, 1997, 16(4): 5. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYH704.008.htm ZHOU Jianxun, WANG Liping, XU Zhiliang. Image information transfer of image intensified night vision system[J]. Journal of Infrared & Millimeter Waves, 1997, 16(4): 5. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYH704.008.htm
[7] 汪丽, 田维坚, 王耀祥, 等. 光锥与CCD耦合器件调制传递函数的测试方法[J]. 光子学报, 2005, 34(4): 3. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB20050400X.htm WANG Li, TIAN Weijian, WANG Yaoxiang, et al. Measurement method for modulation transfer function of fiber optic taper coupled with CCD[J]. Acta Photonica Sinica, 2005, 34(4): 3. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB20050400X.htm
[8] 王耀祥, 田维坚, 黄琨, 等. 光锥与CCD耦合效率的理论分析[J]. 光子学报, 2004, 33(3): 4. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB200403016.htm WANG Yaoxiang, TIAN Weijian, HUANG Kun, et al. Theoretical analysis of the coupling efficient between fiber taper and CCD[J]. Acta Photonica Sinica, 2004, 33(3): 4. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB200403016.htm
[9] 程万胜, 赵杰, 菜鹤皋. CCD像素响应非均匀的校正方法[J]. 光学精密工程, 2008, 16(2): 103-108. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM200802024.htm CHENG Wansheng, ZHAO Jie, CAI Hegao. Correction method for pixel response nonuniformity of CCD[J]. Optics and Precision Engineering, 2008, 16(2): 103-108. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM200802024.htm
[10] 吴厚德, 许文海. 多输出CCD接缝的校正[J]. 光学精密工程, 2013, 21(2): 454-461. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201302031.htm WU Houde, XU Wenhai. Correction of seams for multiple output CCD[J]. Optics and Precision Engineering, 2013, 21(2): 454-461. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201302031.htm
[11] 袁鹏程. 基于像素光子转换特性的图传感器校正方法研究[D]. 哈尔滨: 哈尔滨工程大学, 2019. YUAN Pengcheng. Research on Image Sensor Correction Method Based on Pixel Transfer Characteristics[D]. Harbin: Harbin Engineering University, 2019.
[12] 林祖伦, 王小菊. 光电成像导论[M]. 北京: 国防工业出版社, 2019 LIN Zulun, WANG Xiaoju. Introduction to Optoelectronic Imaging[M]. Beijing: National Defense Industry Press, 2019.
[13] 阴浩, 刘广荣, 金伟其, 等. 像增强型CCD成像系统的分辨力分析[J]. 光子学报, 2010, 39(S1): 96-100. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB2010S1023.htm YIN Hao, LIU Guangrong, JIN Qiwei, et al. Resolution analysis of image enhanced CCD imaging system[J]. Acta Photonica Sinica, 2010, 39(S1): 96-100. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB2010S1023.htm