胡朝龙,刘 磊,钱芸生,莫涌超,拜晓锋,石 峰
(1.南京理工大学 电子工程与光电技术学院,江苏 南京 210094;2.微光夜视技术重点实验室,陕西 西安 710065)
夜间低照度下的目标检测在军事领域一直有着重要意义。微光夜视技术是在夜天光条件下,光信息与电信息的互相转换、处理、加强、显示等过程实现的一门新的技术[1]。在微光夜视技术的发展过程中,微光夜视系统的性能评估一直是一个关键问题。
最早由ROSE A 提出了基于光量子噪声与阈值信噪比理论构建了理想的成像系统探测模型,将人眼观察视角、亮度、对比度三大要素和信噪比联系在一起[2]。1977 年,Richard J C 等结合像增强器和光电阴极的光谱响应建立了经典微光成像系统的理论视距方程[3]。1991 年,G E Giakoumakis 计算了常见的荧光粉和典型的光电阴极的光谱匹配系数,还给出了典型的白炽灯和He-Ne 激光以及人眼与荧光粉的匹配系数[4]。2008 年,张灿林等计算了两种背光照明电子倍增CCD 与景物反射光谱匹配系数[5],这两种CCD 都使用的是三代光电阴极材料。2018 年,刘乾坤计算了超二代光电阴极材料Super S251、Super S252和New S25与绿色草木反射光谱的光谱匹配系数[6],说明了光谱匹配系数对微光夜视系统的影响。2022 年,张琴等人计算了在三种夜天光环境下超二代和三代微光像增强器光阴极与A 光源的光谱匹配系数[7],并得到了这几种夜天辐射条件下的积分灵敏度转换系数。
上述研究中,光谱匹配系数可以反映出光电阴极的光谱响应与景物反射辐射光谱或光源辐射光谱的吻合程度,但仅限于光谱层面上进行比较,并不能直接反映出光电阴极探测景物的性能。本文通过比较超二代和三代微光夜视器件的光电阴极与景物的光谱匹配系数,验证了匹配系数的局限性。针对这一问题,本文结合光谱转换系数、背景与目标的初始对比度以及光电阴极的积分灵敏度,建立了光电阴极的评价因子方程,该评价因子公式可以反映出光电阴极探测景物的性能。传统的视距方程也可以反映出不同光电阴极的探测性能,但该方程计算参数非常复杂,评价因子方程简化了计算过程,在比较不同光电阴极在探测景物性能方面更有优势。
在微光夜视系统分析中,光谱分布的影响是一个重要的方面,良好的光谱匹配是保证和改善仪器性能的重要途径。从目标到人眼的一系列组元中,自然光源与人造光源、景物反射率、大气透过率、仪器的光电阴极直至人眼,都有各自不同的光谱分布。其中,光电阴极的光谱响应分布和目标的反射光谱分布之间的光谱匹配系数是一个重要的参数。
根据光电阴极的光谱响应特性曲线可以得到光谱响应率Sλ,将Sλ对其最大值Smax归一化,可得相对光谱响应率为[8]
此时,光电阴极面接收的是夜天光经过目标反射后的辐射,即:
式中:Pλ为目标反射辐射光谱分布;ρλ为目标的光谱反射系数,它随波长λ 而变化;Kλ为夜天光辐射光谱分布。将目标的反射辐射光谱分布Pλ对其最大值Pmax归一化,则得其相对光谱分布为
由此,定义光谱匹配系数[9]
由光谱匹配系数表达式(4)可知,它只与光电阴极的相对光谱响应率和目标的相对反射光谱分布有关。目标的相对反射光谱分布由式(3)可以获得,这里取 λ在400 nm~1 000 nm 范围内对其最大值的归一化。光电阴极的光谱响应可以通过实验曲线获得。本文讨论的是GaAs 和New S25光电阴极,New S25是美国夜视实验室研制的多碱光电阴极[6],属于超二代光电阴极;GaAs 是标准三代光电阴极[10]。它们的光谱响应曲线如图1 所示。
图1 两种光电阴极的光谱响应曲线Fig.1 Spectral response curves of two photocathodes
夜天光光谱成分相对复杂,包含了月光、星光、大气辉光等多种自然辐射源。满月晴朗天气条件的夜天光辐射光谱主要取决于月光,而星光和大气辉光则是无月天气条件的主要成分,但由于大气辉光主要集中在红外波段,即像增强器工作波段的夜天辐射主要取决于星光的光谱分布[11]。
图2 是满月光和晴朗星空下的夜天空光谱辐射曲线[12]。图3 是两种背景和三种目标的反射率曲线[13],背景环境是山石和水泥路面,目标分别为涤纶军装、军用篷布和军用吉普车。通过表达式(2)可以计算出夜天光经过景物反射后的辐射光谱分布。
图3 不同景物的反射率曲线Fig.3 Reflectance curves for different targets
通过公式(4)计算出不同夜天光环境下两种光电阴极材料与不同景物之间的光谱匹配系数。图4 是山石背景的反射光谱与光电阴极的匹配结果。图5 是涤纶军装反射光谱与光电阴极的匹配结果。图4、图5 中,曲线S1 和S2 分别表示New S25和GaAs 光电阴极相对光谱响应率,曲线P 表示景物在满月光和晴朗星光下的相对反射光谱分布,则光谱匹配系数α(S,t)即为相应的曲线S·P 与横坐标所包围的面积除以曲线P 与横坐标所包围的面积的结果。从表1 中可以看出:相同条件下,New S25与景物反射光谱的匹配系数都要大于GaAs 的匹配系数。原因是New S25材料的光谱响应的半宽度要远大于GaAs 光谱响应的半宽度,并且New S25的光谱响应曲线在所选的光谱范围内更加平滑,GaAs只有在500 nm~900 nm 范围内较为平滑。New S25与景物在满月光时的光谱匹配系数都要远大于在晴朗星空时的光谱匹配系数,这是因为New S25的光谱响应峰值波长离满月光下景物的反射光谱峰值波长更近,另外景物对月光的反射辐射曲线相比对晴朗星空的反射辐射曲线更加平坦。虽然在满月光时GaAs 与景物的匹配系数基本大于在晴朗星空时的匹配系数,但差距不大,这是由于GaAs 的光谱响应峰值波长在近红外范围内。
表1 光谱匹配系数计算结果Table 1 Calculation results of spectral matching factor
图4 两种光电阴极与山石的光谱匹配系数Fig.4 Calculation of spectral matching factors between two types of photocathodes and rock
图5 两种光电阴极与涤纶军装的光谱匹配系数Fig.5 Calculation of spectral matching factors between two types of photocathodes and polyester military uniforms
光谱匹配系数可以体现出光电阴极在波长范围内与景物的反射光谱的吻合程度,光谱匹配系数越大,两者之间的吻合程度越高。光电阴极的响应范围应与景物的光谱辐射波段尽量保持一致,如果匹配系数为0,则两者的响应范围完全没有重合,则该光电阴极对该光源没有探测能力。根据表1 中光谱匹配系数的计算结果,我们可以得到超二代光电阴极比三代光电阴极与景物反射光谱的吻合程度更好,但是并不能说明New S25在探测景物上的性能比GaAs 探测性能好。事实上,在微光成像系统中三代光电阴极比超二代光电阴极的性能好得多,这就体现出了光谱匹配系数的局限性,在同一代光电阴极材料中通过计算光谱匹配系数可以反映出它们性能的差异,因为同一代光电阴极的积分灵敏度差距并不大。可是不同代的光电阴极的积分灵敏度差异很大,三代光电阴极的积分灵敏度显著大于超二代光电阴极的积分灵敏度,因此,比较不同代的光电阴极与景物反射光谱的匹配系数并没有实际意义。想要直观地看出不同光电阴极探测景物的性能,需要结合光电阴极的其他参数提出一种新的评价因子方程。
文献[14]中微光像增强器的空间分辨率公式为
式中:C0是背景与目标的初始对比度;Cd是大气对比衰减系数;M(Ak)为全系统对应空间频率为Ak时的调制传递函数;φmin是阈值信噪比;fo′是物镜焦距;β(G)是亮度增益修正系数;D是物镜直径;τo是物镜的透过率;τd是大气透过率;E0是夜天空照度;SA是光电阴极的积分灵敏度;αλ是标准光源转换到景物反射辐射的光谱转换系数;t是系统累积时间;ε是目标长宽比;是背景和目标所对应的平均反射率;FΦ是噪声功率因子。
式(5)的参数中与光电阴极相关的有背景与目标的初始对比度C0、光电阴极的积分灵敏度SA和标准光源转换到景物反射辐射的光谱转换系数αλ。因此建立一个关于光电阴极的评价因子,需要结合3 个参数进行计算,定义光电阴极的评价因子δ为
转换系数 αλ有如下表示:
式中α(S,A)和α(e,A)分别为[14]
式中A(λ)表示标准光源的相对光谱辐射通量。
式(7)中α-(e,t)是光适应人眼与景物的光谱匹配系数(光适应人眼与目标光谱匹配系数α(e,t1)和光适应人眼与背景光谱匹配系数α(e,t2)的平均值):
式中V(λ)为光适应人眼的相对视见函数。
光电阴极对标准光源的积分灵敏度SA有如下表示[15]:
式中Sm为光电阴极的光谱响应峰值。
背景与目标的初始对比度C0有如下表示[15]:
公式(13)并不能准确地表达目标与背景的初始对比度,该表达式的数学意义为两条曲线与坐标轴的面积之差除以面积之和,这种方式并不能体现出两条曲线的对比度,可能会存在两条曲线不重合,但与坐标轴围成的面积却相等的情况。如图6 所示,如果利用公式(13)求图中两条曲线的对比度,则得到的结果应该是:
图6 两条假设的光谱曲线Fig.6 Two hypothetical spectral curves
此时得到的对比度是0,但是两条曲线并不重合,显然公式(13)并不能准确表示出不同曲线的对比度。因此本文对目标与背景的初始对比度公式进行了改进,即:
通过公式(15)求图6 中两条曲线的对比度,得到的结果是:
此时得到的结果可以表示两条曲线的对比度,体现出了两条曲线的差异性。在满月光时用表达式(13)计算GaAs 光电阴极作用于山石与涤纶军装的初始对比度为0.007 2,用改进后的表达式(15)计算得到的结果为0.220 1。说明初始对比度的计算结果明显变大,更符合实际情况。
将公式(7)、公式(12)、公式(15)带入表达式(6)中并进行化简,为方便计算可以去掉常数项,对评价因子没有影响。微光夜视器件中光电阴极对景物探测性能的评价因子表达式如下:
图7 人眼对光适应的相对视见函数Fig.7 Relative visual function of light adaptation to human eyes
根据公式(17)分别计算出满月光和晴朗星空条件下不同光电阴极与景物之间的评价因子,得到表2~表5 的计算结果。
表2 满月光条件下New S25 光电阴极的评价因子Table 2 Evaluation factors of New S25 photocathode under full moonlight conditions
表3 晴朗星空条件下NewS25 光电阴极的评价因子Table 3 Evaluation factors of New S25 photocathode under clear starlight conditions
表4 满月光条件下GaAs 光电阴极的评价因子Table 4 Evaluation factors of GaAs photocathode under full moonlight conditions
对比表2~表5 的数据可以发现大部分情况下GaAs 的评价因子都要大于New S25,这是因为三代光电阴极GaAs 的最大光谱响应要远远大于超二代光电阴极New S25的最大光谱响应,并且GaAs 的积分灵敏度也同样大于New S25的积分灵敏度。
通过结合背景与目标的初始对比度、光电阴极的积分灵敏度和光谱匹配系数构建的评价因子方程,考虑了在微光成像系统中所有光电阴极相关的参数,这使评价因子更具严谨性。相关研究表明,在不同夜天光环境下超二代的识别距离都明显逊色于三代[16]。通过计算不同光电阴极与景物反射光谱的评价因子,可以看出评价因子的大小符合光电阴极探测景物性能的规律,即三代光电阴极在探测景物的性能上要优于超二代光电阴极的探测性能。所以,GaAs 与景物反射光谱的评价因子大于New S25与景物反射光谱的评价因子是准确的,在微光成像系统中评价因子可以用来作为光电阴极探测性能的判断依据。
本文对两种光电阴极材料和不同景物进行了光谱匹配系数计算,改进了目标与背景的初始对比度公式,并在微光夜视系统下建立了一个光电阴极的评价因子,用来描述不同光电阴极材料对探测目标的性能。通过对比分析发现以下结论。
1)New S25光电阴极与景物在满月光下的光谱匹配系数都要大于晴朗星光下的情况。这是因为在满月光下景物反射光谱的峰值波长离New S25光谱响应峰值波长较近,而景物在晴朗星光下的反射光谱的峰值波长远大于该光电阴极的光谱响应峰值波长。GaAs 的光谱响应峰值波长在近红外范围内,因此GaAs 与景物的光谱匹配系数在两种夜天光条件下差距不大。
2)New S25光电阴极在任何一种情况下的光谱匹配系数都大于GaAs。这是因为New S25光谱响应在所选取的光谱范围内要平坦得多,而GaAs只在650 nm~900 nm 之间有较高的响应,其余部分下降很快。因此New S25光电阴极与景物反射光谱的匹配效果更好。不过三代光电阴极比超二代光电阴极在探测性能上更好,这反映出了光谱匹配系数的局限性。对不同代的光电阴极与景物反射光谱的匹配系数进行比较,并不能说明光谱匹配系数大的光电阴极探测性能就更好。
3)在相同条件下GaAs 的评价因子普遍大于New S25的评价因子,这是因为GaAs 光电阴极的积分灵敏度要远大于New S25的积分灵敏度,且GaAs 的峰值光谱灵敏度要大于New S25的峰值光谱灵敏度。这说明了三代光电阴极性能优于超二代光电阴极,验证了评价因子的可靠性。通过对比不同的光电阴极的评价因子,可以直接比较出光电阴极探测性能的优劣,并且比计算视距公式要简单得多。
微光夜视仪的发展必然离不开光电阴极材料的进步,给出了光电阴极的发展方向,未来的光电阴极的光谱响应应该有较大的半宽度,并且光谱响应应该向近红外波段靠近,可以与景物得到更好的匹配系数。另外,通过光电阴极的评价因子公式可以看出,想要提高评价因子,就需要光电阴极有更高的积分灵敏度和峰值光谱灵敏度。