Calculation of Parameters for Long Wave Infrared FPA Detectors Applied in Low-temperature Background
-
摘要: 对于在低温背景条件下应用的红外焦平面探测器,在其设计阶段首先需要确认的是探测器在低温背景下的主要性能参数是否满足应用需求。通过单元探测器的基本理论公式,从性能测试的角度,推导出了焦平面探测器主要性能指标的理论计算公式,包括峰值响应率、噪声、峰值探测率、动态范围等。提出了红外焦平面探测器主要性能参数计算的设计流程,并通过实例对长波红外焦平面探测器在低温背景下的性能参数进行了计算、设计及验证,实测结果与理论计算符合较好。表明了理论公式及设计流程具有较好的实用性,可为红外焦平面探测器的应用设计提供参考。Abstract: For infrared FPA detectors applied to a low-temperature background, the first aspect that must be confirmed in the design stage is whether the performance of the detectors at low temperatures meets the application requirements. Using the basic theoretical formula of the unit detector, from the perspective of testing, the main performance parameters of the theoretical calculation formula for FPA detectors were deduced, including the peak responsivity, noise, peak detectivity, and dynamic range. In addition, flowchart of the calculation design for the detectors was proposed. The main parameters of the long-wave infrared FPA detector at low temperatures were designed, calculated, and verified using an example, and the results are in good agreement with the theoretical calculations. Thus, the theoretical formula and design calculation flow are practical and can provide a reference for designing infrared FPA detectors.
-
0. 引言
目前,类似于微光夜视仪的直视型微光装备可靠性试验主要通过“拷机”即静态工作过程来实现。试验时,参试人员间歇性地观察被试品显示端图像,通过人工判读被试品目镜视场的方法来确定被试品是否处于正常工作状态[1]。此方法存在以下不足:一是人工试验鉴定效率低且易发生漏检。可靠性试验具有工作量大、重复性高的特点。由于设备批量大,采用人工方法需要多名检测人员每天长时间工作、连续几个月以上才能完成,耗时耗力。而且,试验中人员间隔对设备进行检查,检查间隙内出现的偶发故障则很难被发现。加之检查很大程度取决于检查人员技术水平和检查经验,主观性大,极易发生漏检。二是人工试验鉴定缺乏有效的故障数据存储处理,不利于技术的持续改进。以人工检测为主的试验方法主要依靠肉眼识别,对出现的故障进行简短的文字描述与记载。这样的检测记录方式存在两个问题:一方面,人工检测记录在装备出现故障时难以完备具体、定性定量地记录故障对装备各个部件的影响,在提供故障检测报告时缺乏有力佐证,且由于故障出现的随机性,难以复现;另一方面,人工检测记录的方式缺乏有效的检测数据存储,无法形成故障数据样本集,这为技术改进增加了难度。
综上所述,现有直视型装备的可靠性试验所存在的问题归根结底是检测手段无法满足日益增长的检测需求,因此亟需一种智能的、高效的替代人工检查的技术[2-3]。从图像中提取特征获得场景信息是区域监视、目标识别等领域的重要手段。图像质量评价方法可通过衡量图像变化程度来实现视场异常检测[4-5]。传统的图像质量评价方法包括利用图像像素偏离量度量图像下降质量的方法[6],如均方误差度量(Mean-Square Error,MSE)、峰值信噪比(Peak Signal to Noise Ratio,PSNR)等;计算图像像素灰度幅度变化(等价为图像边缘信息)衡量图像质量的方法[7];利用图像灰度分布计算信息熵的方法[8]等。基于小波变换[9]与结构相似性(Structural SIMilarity,SSIM)理论[10-12]的评估方法与人类视觉系统类似,其出现使图像质量评价的准确性有了较大提升。
本文通过设计转接环连接相机与直视型微光装备进行图像的采集获取,采用结构相似性算法进行异常图像的自动检测,从而实现对直视型微光装备在可靠性试验中的故障诊断。
1. 直视型微光装备图像采集技术
基于SSIM算法的直视型微光装备可靠性试验故障诊断通过在被试品目镜上采用连接环固定图像采集设备的方法来获取图像,连续采集记录在被试品可靠性试验过程中目镜里的图像信息。通过设定图像模板,计算被试品监测图像与模板图像的相似度,识别被试品工作过程中的故障,标记、记录并报警,实现被试品工作过程中的故障自动监测。
1.1 被试品与图像采集设备的可靠连接
被试品与图像采集设备间的可靠连接是稳定采集图像的前提。文中利用直视型微光装备目镜上预留的眼罩螺纹与工业相机镜头上的滤镜螺纹进行可靠连接。按照被试品目镜上螺纹规格和相机镜头滤镜螺纹规格设计加工连接环,一端与被试品连接,另一端与视频采集设备镜头连接,即可实现被试品与图像采集设备间的可靠连接。采集图像时,通过手动调节被试品视度、相机焦距和被试品焦距来实现被试品图像的高清采集。在图像采集相机和镜头选定后,适配该相机的转接环尺寸也随即确定。实际应用中不同被试品眼罩转接环尺寸可能有所差异,为保证转接环能够适用多型装备,可预先按照常用眼罩螺纹规格设计加工系列与相机连接的转接环,如图 1所示的M21×0.7-M25.5×0.5,使用时,根据被试品眼罩螺纹规格,选用相匹配的转接环。
1.2 图像采集设备的选择
图像采集设备应尽量满足以下要求:
1)图像采集的效果应尽可能模拟或超过人眼观察效果,保证软件处理的准确性;
2)相机与被试品间的连接不影响被试品的正常操作和整体平衡,相机重量应尽量轻,体积尽量小;
3)实际试验中可能存在多路图像同时采集的需求,为便于后期图像处理和管理,要求相机接口以Ethernet为主。
文中选择高分辨率彩色相机(型号A3600MU60)作为图像采集设备,配置12mm 1:2.0 1/1.8″镜头。表 1展示了其基本参数。
表 1 图像采集设备基本参数Table 1. Basic parameters of image acquisition equipmentParameter Value Resolution 1920×1200 Maximum frame rate 54fps@1920×1200 Camera type Color camera Pixel size 4.8μm×4.8μm Exposure time 59μs-10s Data interface Gigabit Ethernet Power supply Voltage range 5-15 V,Support PoE power supply Power consumption 2.9 W@12VD Lens interface C-Mout 2. 图像异常诊断技术
2.1 SSIM算法
SSIM算法利用图像像素间的相关性,结合亮度l、对比度c以及两个复杂结构信号间的结构变化s对图像质量进行客观量化评估。
定义图像亮度为:
$${\mu _x} = \frac{1}{N}\sum\limits_{i = 1}^N {{x_i}} $$ 式中:xi为图像第i个像素点的灰度值;N为图像的像素总个数。
设x为模板图像信号,y为待测图像信号。采用平均亮度(μx, μy)计算亮度相似性,定义为:
$$l\left( {x,y} \right) = \frac{{2{\mu _x}{\mu _y} + {C_1}}}{{\mu _x^2 + \mu _y^2 + {C_1}}}$$ 采用标准差(σx, σy)来估计对比度相似性,定义为:
$$c\left( {x,y} \right) = \frac{{2{\sigma _x}{\sigma _y} + {C_2}}}{{\sigma _x^2 + \sigma _y^2{\rm{ + }}{C_2}}}$$ 采用$\left( {\frac{{{\mu _i} - {\mu _x}}}{{{\sigma _x}}},\frac{{{y_i} - {\mu _y}}}{{{\sigma _y}}}} \right)$来计算结构相似度s(x, y),定义为:
$$s(x,y) = \frac{{2{\sigma _{xy}} + {C_3}}}{{{\sigma _x} + {\sigma _y} + {C_3}}}$$ 式中:σxy为参考图像信号x和失真图像信号y的相关系数,计算为:
$${\sigma _{xy}} = \frac{1}{{N - 1}}\sum\limits_{i = 1}^N {\left( {{x_i} - {\mu _x}} \right)\left( {{y_i} - {\mu _y}} \right)} $$ 上述公式中的C1、C2、C3为引入常数,避免μx2+μy2或σx2+σy2值接近0时引起公式无意义。
在实际操作过程中,通常引入参数α, β, γ>0调整亮度、对比度和结构性信息的权重,确定较为合适的SSIM值,定义如下:
SSIM=[l(x, y)]α[c(x, y)]β[s(x, y)]γ
在评价被测图像质量时,SSIM算法用一定尺寸的窗口沿图像逐像素地从左上角到右下角移动,每次计算两幅图像中窗口对应子图像的SSIM值。得到的所有子图像的SSIM值的平均即为该待测图像相较于模板图像的相似度,记为MSIMM(X, Y)表示为:
$${\rm{MSIMM}}(X,Y) = \sum\limits_{j = 1}^M {{\rm{SSIM}}\left( {{X_j},{Y_j}} \right)} $$ 式中:M为子图像块的数量。
2.2 基于结构相似性的图像异常诊断软件设计
异常诊断软件主要根据SSIM算法计算得到的待测图像与模板图像的相似度判定待测设备是否出现故障。软件运行情况如图 2,诊断流程主要包括:确定模板图像、选定相似度判决阈值、SSIM算法计算相似度、实时检测筛选异常图像并对异常信息做出警告与记录。其主要步骤如下:
1)确定模板图像:对于直视型微光装备可靠性试验而言,通常在一个大的、较为黑暗条件下的工房中或晚上照度较小的野外进行。SSIM算法需提前确定好模板图像,即提前将直视型微光装备的视场固定。模板图像作为判定被试装备异常的基础,在选取过程中需严格控制相关外部条件。在视场中心放置一个特征明显、占据直视型微光装备目镜视场1/3~2/3视场的物体,同时控制现场外部光线不出现明显变化,保证直视型微光装备空间位置不变,点击软件拍照功能,存储为模板;
2)相似度判决阈值确定:正式试验前,录制一段被试品目镜视场的视频,期间通过人为开关机操作和遮挡目镜镜头等方式模拟故障,将视频逐帧图像与选定的模板进行比较,得到系列相似度值,将其与故障现象对照,确定出故障出现的最大相似度值,进而作为相似度判决阈值;
3)对于筛选出的异常图像,需发出警告信息并记录序号和时间信息,与图像一起存入告警文件夹。
本文基于MATLAB GUI设计开发直视型微光装备可靠性自动诊断软件,界面如图 2所示。
3. 工程实践与结果分析
根据1.2的描述选择图像采集设备,按照图 3所示的方式通过转接环连接图像采集设备与直视型微光装备目镜。被试微光夜视仪开机,通过相机参数调节,使拍摄到的图像满足观察和软件处理需求。本文实践中依据测试实际情况设置相机采样频率为1 Hz。
3.1 模拟故障识别率检验试验
设置模板图像为设备正常工作时相机传回的监视图像,如图 4。
直视型微光装备可靠性试验主要检测一定工作时间内装备是否正常工作,可能出现的故障有:设备非正常关机、图像质量非正常变化等[3, 13]。故障发生时,监视器中的图像主要表现为:视场出现黑斑、十字刻度消失、黑屏。
采集过程中,通过人为开关机操作和遮挡被试品目镜镜头等方式模拟微光夜视仪可靠性试验过程中出现故障的情况,可覆盖故障时视场变化情况。共获得93幅图像,其中故障图像15幅,如图 5(下标为对应图像在模拟试验获取的图像序列中的编号),构成本文装备可靠性试验工作图像数据集。
利用SSIM算法进行图像相似度计算,采集的监视图像与模板图像相似度计算结果如图 6所示;根据图像序号挑选出相似度小于设定阈值(根据SSIM分布特征和历史经验数据提前确定,本文取0.5)的图像,用“*”表示,这些图像正是数据集中模拟所得的故障图像。故障检测率达到100%且未出现虚警。
该试验证明,通过合理地设置相似度阈值,SSIM算法能够准确地诊断出可靠性试验过程中,被试装备图像的异常情况。
3.2 室内照度变化对被试品故障诊断的影响
直视型微光装备的可靠性试验一般在暗室条件下进行。由于装备对环境照度变化十分敏感,需要保持室内环境照度相对稳定。实际试验中,物体移动、光源发生变化等都会引起室内环境照度不同程度变化[14-15]。图 7(a)显示了某次试验中模拟人员走动、开关门和使用手机时暗室环境的照度变化曲线,图 7(b)展示了采集的对应图像与模板图像的SSIM值。其中,用“*”标记的点为与模板图像相似度小于0.5的图像。可以看到,微小的照度变化引起了被试品图像的显著变化。图 8(下标为对应图像在模拟试验获取的图像序列中的编号)展示了部分相似度小于0.5也被错误识别为异常的图像。
可见,室内照度变化会显著影响SSIM算法图像诊断的准确性。在实际可靠性试验过程中,应该尽可能保证室内照度的稳定,从而提高算法诊断的准确性。
4. 结束语
本文针对传统的直视型微光装备可靠性试验故障诊断中的不足,讨论了一种基于机器视觉的自动故障诊断方法。该方法通过设计转接环连接工业相机与被试装备目镜,自动采集传输监视图像,采用SSIM算法计算测试图像与模板的相似度从而判断设备是否发生异常。实践表明该方法很好地符合了人工视觉观察判断的结果,可为实现装备可靠性自动故障检测提供有力支撑,同时指出试验环境照度变化对试验结果的显著影响。
-
表 1 长波红外探测器低温背景下主要性能参数计算结果
Table 1 Calculation results of main parameters of long wave infrared FPA detector in low temperature background
Capacitance/pF Dark
current/nAHalf well integration time/ms Response signal/V Background noise of T0 /mV Dark current noise /mV Total noise/mV Peak responsivity/(V·W-1) Peak detectivity/
(cmHz1/2/W)Dynamic range/dB 1.0 0.5 0.57 0.919 0.441 0.213 0.538 2.73E+09 4.53E+11 74.9 1.5 0.5 0.85 0.919 0.360 0.174 0.458 2.74E+09 4.34E+11 76.3 2.0 0.5 1.13 0.919 0.312 0.151 0.412 2.74E+09 4.18E+11 77.2 2.5 0.5 1.42 0.919 0.279 0.135 0.382 2.74E+09 4.03E+11 77.9 3.0 0.5 1.70 0.919 0.255 0.123 0.361 2.74E+09 3.90E+11 78.4 3.5 0.5 1.98 0.919 0.236 0.114 0.344 2.74E+09 3.78E+11 78.8 4.0 0.5 2.27 0.919 0.221 0.106 0.332 2.74E+09 3.67E+11 79.1 1.0 2 0.36 0.587 0.353 0.340 0.539 1.75E+09 3.62E+11 74.9 1.5 2 0.54 0.587 0.288 0.278 0.458 1.75E+09 3.47E+11 76.3 2.0 2 0.72 0.587 0.249 0.241 0.412 1.75E+09 3.34E+11 77.2 2.5 2 0.90 0.587 0.223 0.215 0.382 1.75E+09 3.22E+11 77.9 3.0 2 1.09 0.587 0.203 0.196 0.361 1.75E+09 3.12E+11 78.4 3.5 2 1.27 0.587 0.188 0.182 0.344 1.75E+09 3.02E+11 78.8 4.0 2 1.45 0.587 0.176 0.170 0.332 1.75E+09 2.94E+11 79.1 表 2 典型长波红外探测器主要性能参数测试结果
Table 2 Main test results of typical long wave infrared FPA detector
Blackbody temperature/K Half well
integration time/msPeak responsivity/(V·W-1) Peak detectivity/(cmHz1/2/W) Dynamic range/dB 200-215 1.5 2.30E+09 3.21E+11 78.2 293-308 0.17 2.73E+08 1.24E+11 79.0 -
[1] 丁瑞军, 梁平治, 唐红兰, 等. 红外焦平面阵列参数测试方法GB/T 17444-2013[S]. 北京: 中国标准出版社, 2014. DING Ruijun, LIANG Pingzhi, TANG Honglan, et al. Measuring methods for parameters of infrared focal plane arrays GB/T 17444-2013[S]. Beijing: China Standard Press, 2014.
[2] Campos T. Test Bench for Infrared Detectors[C]//Proc. of SPIE, 2005, 5640: 183-192.
[3] 王群, 朱牡丹, 赵亮. 红外探测器参数测试研究[J]. 红外技术, 2006, 28(10): 599-601. DOI: 10.3969/j.issn.1001-8891.2006.10.010 WANG Qun, ZHU Mudan, ZHAO Liang. Parameter test of infrared detectors[J]. Infrared Technology, 2006, 28(10): 599-601. DOI: 10.3969/j.issn.1001-8891.2006.10.010
[4] 白丕绩, 姚立斌. 第三代红外焦平面探测器读出电路[J]. 红外技术, 2015, 37(2): 89-96. DOI: 10.11846/j.issn.1001_8891.201502001 BAI Piji, YAO Libin. Read out integrated circuit for third-generation infrared focal plane detector[J]. Infrared Technology, 2015, 37(2): 89-96. DOI: 10.11846/j.issn.1001_8891.201502001
[5] 翟永成, 丁瑞军. 320×256大电荷容量的长波红外读出电路结构设计[J]. 红外与激光工程, 2016, 45(9): 0904003-1-6. DOI: 10.3788/IRLA201645.0904003 ZHAI Yongcheng, DING Ruijun. 320×256 LW IRFPA ROIC with large charge capacity[J]. Infrared and Laser Engineering, 2016, 45(9): 0904003-1-6. DOI: 10.3788/IRLA201645.0904003
[6] 景松, 杨波, 黄张成, 等. 超高增益低噪声红外焦平面读出电路研究[J]. 红外技术, 2019, 41(12): 1117-1123. http://hwjs.nvir.cn/article/id/hwjs201912005 JING Song, YANG Bo, HUANG Zhangcheng, et al. Study on high-gain and low noise infrared focal plane readout circuit[J]. Infrared Technology, 2019, 41(12): 1117-1123. http://hwjs.nvir.cn/article/id/hwjs201912005
[7] Bertrand Terrier, Anne Delannoy, Philippe Chorier. LWIR and VLWIR detectors development at SOFRADIR[C]// Proc. of SPIE, 2010, 7826: 78261J-1-12.
[8] Willardson R K, Albert C Beer. Semiconductors and Semimetals, V18 Mercury Cadmium Telluride[M]. New York: ACADEMIC PRESS, A Subsidiary of Harcourt Brace Jovanovich, Publishers, 1981.
[9] Vincent J D. Fundamentals of Infrared Detector Operation and Testing[M]. New York: A Wiley-Interscience Publication, John Wiley & Sons, 1990.
[10] 吴宗凡, 柳美琳, 张绍举, 等. 红外与微光技术[M]. 北京: 国防工业出版社, 1998. WU Zongfan, LIU Meilin, ZHAO Shaoju, et al. Infrared and Night Vision Technology[M]. Beijing: National Defence of Industry Press, 1998.
[11] 褚君浩. 窄禁带半导体物理学[M]. 北京: 科学出版社, 2005. CHU Junhao. Physics of Narrow Bandgap Semiconductor[M]. Beijing: Science Press, 2005.
[12] 吴东生. 光子探测器的光子噪声[J]. 长春光学精密机械学院学报, 1983(3): 42-46. https://www.cnki.com.cn/Article/CJFDTOTAL-CGJM198303005.htm WU Dongshen. Photon noise of photon detector[J]. Journal of Changchun Institute of Optics and Precision Mechanics, 1983(3): 42-46. https://www.cnki.com.cn/Article/CJFDTOTAL-CGJM198303005.htm
[13] 徐丽娜, 董杰, 戴立群, 等. 红外探测器探测率的研究与分析[J]. 激光与红外, 2018, 48(4): 497-502. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201804016.htm XU Lina, DONG Jie, DAI Liqun, et al. Study and analysis of the detectivity of infrared detector[J]. Laser & Infrared, 2018, 48(4): 497-502. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201804016.htm
[14] Antoni Rogalski. Infrared Detectors[M]. Second Edition, Boca Raton: CRC Press, 2011.
[15] Shkedy L, Armon E, Avnon E, et al. Hot MWIR detector with 5 μm pitch[C]//Proc. of SPIE, 2021, 11741: 117410W1-11.
[16] 毛京湘, 庄继胜, 陈国华. 红外焦平面探测器性能测量中g因子的计算[J]. 红外, 2008, 29(4): 9-11. https://www.cnki.com.cn/Article/CJFDTOTAL-HWAI200804004.htm MAO Jingxiang, ZHUANG Jisheng, CHENG Guohua. Calculation of factor g in measurement of parameters of IRFPA detectors[J]. Infrared, 2008, 29(4): 9-11. https://www.cnki.com.cn/Article/CJFDTOTAL-HWAI200804004.htm