苏 悦,拜晓锋,2,党小刚,2,冯丹青,2,程宏昌,2,李周奎,2,韩 坤,2
(1.昆明物理研究所,云南 昆明,650223;2.微光夜视技术重点实验室,陕西 西安 710065)
目前针对强光照射条件下目标背景对比度的研究较为广泛。然而,已有研究结果表明,随着环境照度的降低,若光照的色温条件发生变化[1],则其目标背景对比度也会发生变化,因此强光照射条件下的目标背景对比度不再适用,或者说会产生较大偏差,而关于夜间目标背景对比度对视距的影响研究却鲜有报道。对于微光夜视系统而言,其典型应用环境多为星光条件,即10-3lx 环境。像增强器在微光领域对图像信息倍增作用被广泛使用,所以采取用微光像增强器和CCD 通过光纤光锥或透镜进行耦合[2],利用像增强器的微光下的增强性能,使得CCD 在低照度环境下工作,并且可以输出稳定图像。
本文利用微光像增强器的高亮度增益放大特性,对星光条件下的目标背景对比度进行测量及验证。针对像增强器亮度增益特性表现出在一定范围内输入照度和输出亮度为线性变化的现象,对线性变化范围内不同亮度增益下的照度-灰度关系以及照度-目标背景对比度关系进行研究。由此对亮度增益在微光环境下以微光像增强器为核心的目标背景对比度的影响展开讨论。满足微光夜视整机视距评估时对夜间星光环境下基于国产三代ICCD 的目标背景对比度特性的准确测量需求。
由微光环境下物体辐射光传递过程,利用目标对比度理论公式,以国产ICCD 作为测试器件,建立ICCD 目标对比度测试数学模型,图1所示为微光环境对比度测试示意图[3-4]。
图1 对比度测试示意图Fig.1 Schematic diagram of contrast test
利用灰度反应目标与背景调制对比度的公式为[5]:
式中:LM为探测环境的目标亮度值;LB为探测目标的亮度值;gB为探测环境的背景灰度值;gM为探测目标的灰度值;C为调制对比度。
测试系统中的光信息传递分为3 个主要过程:一是光源光束在目标处反射并入射到探测仪器中;二是在ICCD 中通过像增强器光放大;三是通过耦合的CCD 器件将光信号转化为灰度信息。
第一个传输过程主要影响因素为在光传播过程中大气的衰减损耗以及物镜的衰减损耗τobj,光源在物体表面发生反射,反射光成像到像增强器输入面上。光谱反射率ρ表征物体的颜色信息,表现了颜色材质特征,是物体的固有属性,目标探测时只与物体本身的相关,无关于环境照度[6]。目标和背景的光谱反射率不同形成差异。经过物镜成像,反射光信息成像到像增强器输入面上的辐射照度为[3]:
式中:EK为目标反射光到达物镜表面的辐射照度;D为物镜孔径大小;f为物镜的焦距;LA为目标反射光到达物镜表面的辐射亮度。
第二个传输过程中,光信息在ICCD 中通过像增强器光放大。如图2所示,三代微光像增强器ABC曲线在最大输出亮度Em(点D所对应位置)和最小输出亮度ES(图2点A所对应位置)区间内表现出近似线性变化,即亮度增益为定值。C点(B点)处于线性区域内,由于接近上限点D(下限点A),输出信号出现上限饱和(下限截断)的现象[3]。当输入照度大于最大输出照度或小于最小输出照度时呈现出非线性变化。线性区间内亮度增益表达式为:
图2 三代微光像增强器的ABC 曲线Fig.2 ABC curve of the third generation low light level image intensifier
第三个传输过程中,通过耦合的CCD 器件将光信号转化为灰度信息。像增强器输出辐射亮度为光锥输入辐射亮度。考虑光纤光锥耦合对光信息带来的损耗以及CCD 工作时的光-电-光变化过程,光纤光锥的有效透过率T为下式[7-9]:
式中:Φout光纤光锥小端面输出的光通量;Φin为光纤光锥大端面的输入辐射通量;Dsmall光纤光锥小端面直径;Dlarge光纤光锥大端面直径;TA光纤光锥材料的透过率;TR光锥光纤反射损耗后的透过率;KC为填充效率。
CCD 接收光纤光锥输出光信息,将光信息转换为灰度图像,整个输出图像的灰度分布,是由数个单个像素共同组成的。每个像素接收到的辐射照度Ee为:
式中:s光敏面为光敏面面积;τ粘合剂为粘合剂损耗;s为像素面积。
根据电荷转移原理,可以得到单个像素的电荷量Q为式(6):
式中:η为光电转换效率;e为电子电荷量(1e=1.602×10-19C);h为普朗克常数(h=6.62607015×10-34J·s);ν为辐射电磁波的频率。
综合前文讨论光学传递链的影响因素得到输出图像上对应的灰度值为:
式中:o为输出过程中的直流偏置;σ为输出过程中存在的随机噪声[10-11]。
由于CCD 的像素不均匀性,每一个像素的灰度响应、光电转换效率和像素面积不尽相同。CCD 上任意位置(i,j)处灰度值可表示为:
当选取目标位置为(m,n),背景位置为(u,v)时,对比度可表示为下式:
式中:ρM为目标的光谱反射率;ρB为背景的光谱反射率;νM为目标辐射电磁波频率;νB为背景辐射电磁波频率;EM为目标的输入照度;EB为背景的输入照度;亮度增益G在非线性变化区间表现为以输入照度为自变量的函数。
式(9)中可以看出,当照度在线性变化区间(ES<E<Em)内时,对比度变化不受亮度增益影响,当照度大于像增强器的最大输入照度或低于最小输入照度时,对比度变化与亮度增益有关。
目标对比度测试系统中,主要采用高信噪比、高分辨率的CCD 相机、国产三代微光像增强器、光纤光锥、照度计、经过标定的靶板以及计算机等。光源采用2856 K 卤钨灯来模拟夜天光环境;选用照度范围可覆盖1×10-1lx~1×10-4lx 范围的宽量程微弱光照度计进行测试;靶标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所示。
图3 对比度测试系统示意图Fig.3 Schematic diagram of contrast test system
由式(1)所示,实际测试中图像的对比度由目标灰度与背景灰度计算所得,通过光学信息传递所接收得到灰度图像。以亮度增益为12000 cd/(m2·lx)条件为例,可以看出对于相同靶标在不同灰度时,靶标上的黑白条纹的灰度值均表现为随着照度的增加而递增,图4为照度和灰度对应图线。从图4可以看出,当照度大于一定值时,照度-目标灰度呈现非线性变化关系,根据这种现象对G=12000 cd/(m2·lx)进行线性区间拟合。
图4 G=12000 cd/(m2·lx)时照度与目标灰度的关系Fig.4 Relationship between illuminance and target gray level when 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。
图5 照度-目标灰度在10-3 lx 量级拟合直线Fig.5 Illuminance-target gray fitted line in the order of 10-3lx
像增强器与CCD 作为ICCD 主要组成部分,非线性变化出现原因有主要由两方面:像增强器的输出亮度非线性区间以及CCD 的动态范围。其一,对于像增强器输入照度达到最大输入照度值后,空间上出现电荷量饱和,并且由于荧光屏发光能力的有限性和光电阴极疲劳的原因造成像增强器的输入照度-输出亮度之间的非线性变化[12],此时,由光电阴极输入照度和荧光屏输出亮度表征的亮度增益公式无法适用[13]。其二,像素大小以及噪声影响决定了CCD 的动态范围,表现在输出图像灰度和照度之间的非线性变化关系。在对比度计算中,灰度的非线性变化影响对比度变化。
建立目标背景对比度数学模型如图6(a)所示,平面投影如图6(b)所示,可以看出,在照度-灰度所对应的线性区间内,对比度的实际变化量与理想变化量接近。当照度值过低时,对比度出现阶跃变化;当照度值过高时,对比度曲线偏离理想变化曲线。因此对于实际测试中,输入照度过低和过高时,目标对比度均不能正确表征图像的对比度信息。
图6 目标背景对比度数学模型Fig.6 Mathematical model of target background contrast
图 7 中(a)(b)(c)分别为亮度增益为 12000 cd/(m2·lx)、8300 cd/(m2·lx)、6000 cd/(m2·lx)的情况下照度和对比度的对应关系。可以看出曲线均变化平缓,照度和对比度未表现出明显相关性。
图7 不同照度对应目标背景对比度的变化关系Fig.7 Relationship between different illuminance and target background contrast
根据图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 为高度相关。
图8 3 组亮度增益样本的目标背景对比度Pearson 相关系数热力图Fig.8 Target background contrast Pearson correlation of three groups of brightness gain samples
本实验结果说明不同亮度增益条件下对比度表现出中度以上的正相关性,呈现较为良好的相似性。结合式(9)可以看出,在当前的试验条件下,在照度为[2×10-3,8×10-3](lx)区间内,对比度大小变化与亮度增益未表现出明显相关性;当照度小于2×10-3lx时,亮度增益为6000 cd/(m2·lx)条件下采集图像灰度值过低,目标背景对比度数值出现阶跃变化,不能正确表征图像目标与背景的差异程度,亮度增益为6000 cd/(m2·lx)与12000 cd/(m2·lx)、8300 cd/(m2·lx)两组试验结果的Pearson 相关系数表现出低于-0.70 负相关性。当照度大于8×10-3lx 时,在当前试验条件限制下,可得到8300 cd/(m2·lx)、6000 cd/(m2·lx)在[8×10-3,3×10-2](lx)区间内Pearson 相关系数为-0.37,表现出弱负相关性。线性区间外,对比度表现出与输入照度的相关性,亮度增益数值影响目标背景对比度数值。
本文从目标对比度测试原理出发,讨论亮度增益对以微光像增强器为核心的目标背景对比度测试的影响,搭建测试平台,进行暗室实验,对获得实验数据进行整理分析。通过实验得到不同亮度增益下照度-灰度、照度-目标背景对比度的变化特性以及不同亮度增益条件下目标背景对比度之间的相关性。分析结果表明: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 目标背景对比度测试的准确测量提供支撑。