激光损伤非制冷微测辐射热计的实验研究

周冰, 贺宣, 刘贺雄, 李秉璇, 张俊祥

周冰, 贺宣, 刘贺雄, 李秉璇, 张俊祥. 激光损伤非制冷微测辐射热计的实验研究[J]. 红外技术, 2021, 43(1): 1-7.
引用本文: 周冰, 贺宣, 刘贺雄, 李秉璇, 张俊祥. 激光损伤非制冷微测辐射热计的实验研究[J]. 红外技术, 2021, 43(1): 1-7.
ZHOU Bing, HE Xuan, LIU Hexiong, LI Bingxuan, ZHANG Junxiang. Experimental Study on Laser-induced Damage to an Uncooled Microbolometer[J]. Infrared Technology , 2021, 43(1): 1-7.
Citation: ZHOU Bing, HE Xuan, LIU Hexiong, LI Bingxuan, ZHANG Junxiang. Experimental Study on Laser-induced Damage to an Uncooled Microbolometer[J]. Infrared Technology , 2021, 43(1): 1-7.

激光损伤非制冷微测辐射热计的实验研究

详细信息
    作者简介:

    周冰(1976-),男,博士,副教授,硕士生导师,主要从事光电对抗方面的研究。E-mail:zhbgxgc@163.com

  • 中图分类号: TN977

Experimental Study on Laser-induced Damage to an Uncooled Microbolometer

  • 摘要: 在市场需求驱动下,非制冷红外测辐射热计阵列在多个领域正逐步应用,对其激光干扰的研究也被提上了日程。本文以多晶硅探测器为例,结合非制冷微测辐射热计阵列的构造和工作原理,分析了激光辐照作用下的温度响应;通过10.6 μm连续激光的辐照实验,得到了像元阵列进入不同损伤状态对应的激光功率范围,分析了激光损伤的热效应机制,指出干扰面积大于光斑面积是由于热量的“倒灌”;将战场上激光传输的影响简单等效为斩波调制,研究了调制频率和占空比对激光干扰效果的影响。经分析得出,在调制频率100~500 Hz、占空比0.1~0.5的范围内,连续激光对非制冷微测辐射热计在传输等效调制频率较低、占空比较大的条件下能得到更好的干扰效果。
    Abstract: Driven by market demand, uncooled infrared bolometer arrays are gradually becoming popular in several fields, and their laser interference has become a major research topic. In this study, based on the structure and working principles of an uncooled microbolometer array, temperature response under laser irradiation is analyzed. An experiment of 10.6 μm continuous laser irradiation of an uncooled polysilicon detector shows that the pixel array enters different damage states. The laser power range and thermal effect mechanism of laser damage are also analyzed. Previous studies have shown that the interference area is larger than the spot area as a result of "heat inversion" and that the influence of laser transmission on the battlefield is equivalent to chopping modulation. T he effects of modulation frequency and duty cycle are also studied. Analytical results show that a continuous laser can achieve a better interference effect on an uncooled microbolometer with a low transmission modulation frequency and large duty ratio.
  • 夜间高分辨率图像被广泛应用于安防安保监控系统、航天航空军事制导、对地检测卫星遥感以及汽车自动驾驶等领域。然而,有雾的夜间图像经增强处理后呈现出对比度低,细节模糊等特征,不易辨识。目前最广泛使用的基于暗通道先验理论(dark channel prior, DCP)的去雾算法[1]对白天降质图像进行去雾时具有较好的效果,但对夜间降质图像进行去雾后,图像在人工光源区域存在严重的颜色失真[2],暗区域的亮度不够[3]

    目前,基于DCP算法的夜间图像去雾算法的研究已经取得了较大的进展:陈志恒提出了一种自适应全局亮度补偿算法[4],自适应调节夜间有雾图像的全局亮度,提高了夜间降质图像待去雾区域的可视度,但对微弱光源的夜间图像去雾效果有待改进;皮燕燕将夜间有雾图片分为结构层和纹理层,利用DCP算法得到结构层的无雾反射光图像,对纹理层进行增强后得到夜间去雾图像[5],但对图像天空部分的去雾存在失真;王柳哲针对夜晚环境的光照不均匀会影响暗通道的大气光估计,采用中值滤波与加权正则化模型对环境光进行优化的方法来估计大气光值[6],改进了去雾图像的色偏现象,但是不能很好地去除光晕。

    针对现有算法存在的图像颜色失真、微弱光源图像去雾效果差、光晕不能很好去除等问题,本文提出了采用暗态点光源模型的夜间有雾图像去雾算法,针对传统图像透射率分布估计不准导致颜色失真的问题,通过暗通道可信度权值因子对粗透射率分布进行修正,再利用边窗导向滤波对粗透射率分布进行细化。实验结果表明,本文算法去雾效果好、在微弱光源处还原性好,有效减少了还原图像的光晕与颜色失真。

    He通过统计发现,晴朗天气图像的暗通道图像中至少86%以上的灰度值集中分布在[0, 16]的灰度空间范围之内,且在无天空区域图像的R、G、B通道中,至少存在一个通道的灰度值是趋近于0的[1]。He将这一结果定义为暗通道先验,用式(1)进行表示:

    $$ {J_{{\text{dark}}}}(x, y) = \mathop {\min }\limits_{c \in \{ R, G, B\} } \left( {\mathop {\min }\limits_{y \in \mathit{\Omega } (x, y)} \left( {{J_{\text{c}}}(x, y)} \right)} \right) $$ (1)

    式中:Jdark(x, y)表示无雾图像的暗通道图像;Jc(x, y)表示无雾图像的RGB通道中的某一通道;Ω(x, y)表示滤波的窗口大小。

    通过计算大气光值A和透射率分布t″(x, y)解得无雾图像J(x, y)的表达式如式(2)所示:

    $$ J(x, y) = \frac{{I(x, y) - A}}{{t''(x, y)}} + A $$ (2)

    然而,当透射率分布t″(x, y)趋近于零时,会导致还原后的去雾图像J(x, y)包含大量的噪声,表现出复原的图像整体偏白。为避免该类情况的发生,在实际处理中,给透射率分布设置一个下限t0,复原后的降质图像可以表示为式(3):

    $$ J(x, y) = \frac{{I(x, y) - A}}{{\max \left[ {t''(x, y), {t_0}} \right]}} + A $$ (3)

    式中:t0一般取作0.1,相当于去雾图像仍然保留了少量的雾层。

    DCP算法对白天图像进行处理时具有较好的去雾效果,去雾结果颜色真实自然。但对夜间图像进行去雾时,处理后的图像在人工光源区域产生严重的颜色失真,暗区域的亮度不够导致细节信息缺失,光源区域的纹理受损严重。

    针对夜间有雾图像的纹理细节保持、色彩保持等问题,本文在DCP算法的基础上,提出了一种采用暗态点光源模型的夜间有雾图像去雾算法。算法实现过程如图 1所示,主要由夜间环境光分布估计、透射率分布估计以及颜色校正3个模块组成,具体过程如下:

    图  1  采用暗态点光源模型的夜间图像去雾算法流程
    Figure  1.  Process of nighttime image dehazing algorithm using dark point light source model

    利用最小值滤波算法、自适应伽马变换算法建立暗态点光源模型,并利用最大值滤波算法获取夜间降质图像的亮度分布。使用暗态点光源模型、亮度分布以及联合双边滤波算法估计夜间降质图像的环境光分布[7]

    通过颜色空间转换、HSV通道分离,对降质图像的亮度通道和饱和度通道分别进行最大值滤波、最小值滤波处理,利用联合双边滤波对其进行细化处理,得到暗通道可信度权值因子。然后,建立扇形模型,利用高斯均值函数、边界约束、均值不等关系,获取“伪”去雾图像,再结合夜间环境光分布,将获取的暗通道可信度权值因子和伪去雾图像作为修正参数,改进透射率分布的求取公式,获取夜间降质图像的透射率分布。最后利用边窗导向滤波细化纹理,得到精细化透射率分布,再对降质图像进行还原后获取去雾图像。

    根据输入图像G通道的最大灰度值,计算输入图像的景深图像,将其与输入图像B通道、R通道的最大灰度值做减法,求取景深差分图像DDF。再利用φθ两个校正参数分别对图像R通道和B通道进行色偏校正,再将处理后的3个通道进行融合,得到颜色校正的图像。

    在计算夜晚降质图像的透射率分布时,容易引起透射率分布t′(x, y)的估计值偏小[8]。为提高透射率分布估算结果的准确度,减少去雾图像的颜色失真,本文通过暗通道可信度权值因子ω和“伪”去雾图像$\tilde J(x, y)$对图像粗透射率分布进行修正,同时利用边窗导向滤波对其进行细化,得到细透射率分布t″(x, y)。

    暗通道可信度权值因子获取流程如图 2所示,将夜间降质图像I(x, y)转换至HSV颜色空间,并进行通道分离得到饱和度通道S(x, y)、亮度通道V(x, y)图像,分别对其进行最大值滤波、最小值滤波,得到Smax(x, y)和Vmin(x, y)。

    图  2  暗通道可信度权值因子获取流程
    Figure  2.  Dark channel credibility weight factor acquisition process

    Smax(x, y)和Vmin(x, y)图像自身边缘信息损失严重,需用联合双边滤波对其进行纹理细化,得到细化结果Smin′(x, y)和Vmin′(x, y),处理过程如式(4)所示:

    $$ \left\{ \begin{array}{l} {{S'}_{\max }}\left( {x, y} \right) = {\rm{CroBilFilter}}\left( {{S_{\max }}\left( {x, y} \right)} \right) \\ {{V'}_{\min }}\left( {x, y} \right) = {\rm{CroBilFilter}}\left( {{V_{\min }}\left( {x, y} \right)} \right) \end{array} \right. $$ (4)

    此时以细化后Smax′(x, y)和Vmin′(x, y)为处理对象,对暗通道可信度权值因子ω进行求取,求取过程如式(5)、(6)所示。

    $$ \left\{ \begin{array}{l} {\omega _1} = {\left( {1 + {k_1} \cdot \exp \left( {{{S'}_{\max }}\left( {x, y} \right)} \right)} \right)^{ - 1}} \\ {\omega _2} = {\left( {1 + {k_2} \cdot \exp \left( {\frac{{{{V'}_{\min }}\left( {x, y} \right)}}{{A(x, y)}} - 1} \right)} \right)^{ - 1}} \end{array} \right. $$ (5)
    $$ \omega=\omega_1 \cdot \omega_2 $$ (6)

    式中:A(x, y)表示夜间环境光分布;ω1表示饱和度系数;ω2表示亮度系数;ω表示暗通道可信度权值因子。暗通道可信度权值因子ω会随着Smax′(x, y)的增加、Vmin′(x, y)的减少而出现大幅度下降;k1k2表示sigmoid函数的缩放系数,当k1=8、k2=20时,具备较好的实验效果。

    通过构建一个扇形模型,对夜间无雾图像的暗通道图像作初始估计;设定一种自适应高斯均值衰减函数用来估计扇形阴影区域的最小值和最大值;通过几何均值不等关系进行逼近拟合,得到最佳的夜间无雾图像暗通道图,将其称作夜间伪去雾图像,记作$ {\tilde J_{{\rm{night}}}}(x, y) $。

    本文通过边缘最值约定对min(Jc(x, y))的范围进行约束,构建如图 3所示的扇形模型,由图可知:min(Jc(x, y))的最优解存在于min(Icmin(x, y))与min(Icmax(x, y))之间的阴影区域中。

    图  3  扇形模型几何表示
    Figure  3.  Sector model geometric representation

    采用不等式逼近对阴影区域内的最优解min(Jc(x, y))进行求取。采用较为剧烈的衰减ρ2处理上边界min(Icmax(x, y)),对下边界min(Icmax(x, y))采用较为平和的衰减ρ1对其进行估值,以防止逼近过程中产生溢出现象。得到扇形区域两侧边界的估计值如式(7)所示:

    $$ \left\{ \begin{array}{l} \min \left( {{I_{{\text{cmin}}}}(x, y)} \right) = {\rho _1}\min \left( {{I_{{\text{cmin}}}}(x, y)} \right) \\ \min \left( {{I_{{\text{cmax}}}}(x, y)} \right) = {\rho _2}\min \left( {{I_{{\text{cmax}}}}(x, y)} \right) \end{array} \right. $$ (7)

    在得到阴影区域两侧边界的估计值min(Icmin(x, y))和min(Icmax(x, y))后,根据式(8)通过不等式形式进行逼近,从而得到最佳min(Jc(x, y)),本文称作夜间伪去雾图像,并记作${\tilde J_{{\rm{night}}}}(x, y)$,如式(8)所示:

    $$ \begin{array}{l} {{\tilde J}_{{\rm{night}}}}(x, y) = \min \left( {{J_{\text{c}}}(x, y)} \right) \\ \;\; \;\;\; \;\; \;\; \;\; \;\; \;\;= \sqrt {\min \left( {{J_{{\text{cmin}}}}(x, y)} \right) \cdot \min \left( {{J_{{\text{cmax}}}}(x, y)} \right)} \end{array} $$ (8)

    结合夜间环境光分布A(x, y),将暗通道可信度权值因子ω和夜间伪去雾图像${\tilde J_{{\rm{night}}}}(x, y)$作为修正参数,对透射率分布求取公式进行改进,得到式(9):

    $$ t'(x, y) = \frac{{1 - \psi \cdot \left( {\frac{{I\left( {x, y} \right)}}{{A(x, y)}}} \right)}}{{1 - \omega \cdot \left( {\frac{{{{\tilde J}_{{\rm{night}}}}(x, y)}}{{A(x, y)}}} \right)}} $$ (9)

    式中:t′(x, y)为粗透射率分布;ω为暗通道可信度权值因子;ψ为雾层保留参数,用作去雾图像中少量雾气的保留,以保证图像的视觉真实性,对于夜间雾图的处理中,一般取作0.65。

    为减小去雾后出现的光晕现象,本文根据文献[9]提出的边窗导向滤波技术对粗透射率分布进行纹理细化处理,以提升图像边缘信息的保持能力。通过对边窗的滤波核函数F进行加权操作,并将其应用在{L, R, U, D, NW, NE, SW, SE}8个窗口类型上,从输出的8个结果中,找出与输入最接近的边窗类型,将其滤波结果作为最终的边窗滤波结果,得到精细化透射率分布t″(x, y)。

    在得到夜间环境光分布A(x, y)以及细透射率分布t″(x, y)以后,对式(3)进行改进得到式(10):

    $$ J(x, y) = \frac{{I(x, y) - A(x, y)}}{{\max \left[ {t''(x, y), {t_0}} \right]}} + A(x, y) $$ (10)

    为验证本文算法的去雾效果,本文将从主观和客观两个评价角度进行实验对比。

    本文以重度雾霾下的城市夜景实拍等夜间降质图像建立降质图库,随机选取了4幅具有低照度夜间特征代表性的降质图像,基于低照度夜间增强算法对其处理之后,分别采用基于暗通道先验理论去雾算法与本文算法对其进行去雾处理,得到结果如图 4所示。

    图  4  本文算法的去雾结果对比
    Figure  4.  Comparison of dehazing results of the algorithm in this paper

    可以看到,采用基于暗通道先验理论的去雾图像在人工光源区域有严重的光晕效应,而本文的去雾算法的去雾结果中图像中的雾气得到有效去除,去雾后的图像颜色自然,符合人眼视觉观感。

    使用图像信息熵、平均梯度和对比度对经过文献[1]、[6]、[10]中算法与本文算法的处理结果进行客观评价。统计其客观评价结果如表 1~表 3所示,其中,表 1为图像信息熵的统计结果,表 2为图像平均梯度的统计结果,表 3为图像对比度的统计结果。表中,算法1为传统DCP算法[1],算法2为基于多光源模型与暗通道先验的去雾算法[6],算法3为基于暗通道补偿与大气光值改进的去雾算法[10]

    表  1  图像信息熵统计结果
    Table  1.  Image information entropy
    Comparative object Original image Enhanced image Algorithm 1 Algorithm 2 Algorithm 3 Method of this paper
    Fig. 1 6.87 7.15 7.28 7.46 9.04 9.77
    Fig. 2 6.74 7.03 7.26 7.2 7.21 7.61
    Fig. 3 7.14 7.95 7.29 8.44 11.17 13.81
    Fig. 4 6.57 7.28 7.68 8.06 8.98 9.74
    下载: 导出CSV 
    | 显示表格
    表  2  图像平均梯度统计结果
    Table  2.  Image average gradient
    Comparative object Original image Enhanced image Algorithm 1 Algorithm 2 Algorithm 3 Method of this paper
    Fig. 1 4.21 6.59 8.13 10.51 7.19 13.61
    Fig. 2 1.56 2.56 3.56 4.42 7.83 8.04
    Fig. 3 9.77 13.01 14.41 17.79 19.05 22.11
    Fig. 4 5.53 8.92 9.09 14.22 11.95 17.16
    下载: 导出CSV 
    | 显示表格
    表  3  图像对比度统计结果
    Table  3.  Image contrast
    Comparative object Original image Enhanced image Algorithm 1 Algorithm 2 Algorithm 3 Method of this paper
    Fig. 1 61.05 202.34 247.98 365.75 446.84 558.29
    Fig. 2 14.23 23.58 28.97 50.72 47.31 51.79
    Fig. 3 229.85 493.12 581.56 707.65 943.99 976.51
    Fig. 4 82.88 196.79 231.62 320.78 337.73 441.38
    下载: 导出CSV 
    | 显示表格

    可以看出,经过本文算法处理过后的图像的质量普遍有较为明显的改善,主要表现在图像的信息熵值、平均梯度以及对比度均高于对比算法,说明处理后的夜间有雾图像所包含的信息量更为丰富,在纹理细节上损失较小、图像清晰度高,同时夜间降质图像的明暗对比度得到较好地拉伸。

    为解决基于暗通道先验理论的去雾算法在处理夜间有雾图像时细节信息缺失、光源区域的纹理受损严重的问题,本文提出采用暗态点光源模型的夜间有雾图像去雾算法。同时通过暗通道可信度权值因子和“伪”去雾图像对图像粗透射率分布进行修正,再利用边窗导向滤波对其进行细化,提高了透射率分布估算结果的准确度。实验结果表明,本文算法处理后的图像在微弱光源处还原性好,有效减少了还原图像的光晕与颜色失真且在去雾前经过视见度增强处理,降质图像的亮度得到一定程度的保留,可视性较好。

  • 图  1   桥式结构的像元阵列

    Figure  1.   Cell array of bridge structure

    图  2   非制冷微测辐射热计纵向分布

    Figure  2.   Longitudinal distribution of uncooled microbolometers

    图  3   TEC结构图

    Figure  3.   TEC structure diagram

    图  4   实验系统布局

    Figure  4.   Experimental system layout

    图  5   斩波调制检测系统

    Figure  5.   Chopper modulation and detection system

    图  6   实验系统实景图

    Figure  6.   Actual appearance of the experimental platform

    图  7   点饱和现象的热成像

    Figure  7.   Thermal imaging of point saturation

    图  8   中度干扰现象的热成像

    Figure  8.   Thermal imaging of moderate saturation

    图  9   重度饱和现象的热成像

    Figure  9.   Thermal imaging of severe saturation

    图  10   点饱合现象的热成像

    Figure  10.   Thermal imaging of point damage

    图  11   斩波调制激光辐照导致的光斑闪烁现象(斩波调制频率100 Hz,占空比0.5)

    Figure  11.   Spot flashing caused by chopping modulated laser irradiation

    图  12   探测器灰度均值与激光调制频率的曲线拟合

    Figure  12.   Curve fitting of detector gray mean value and laser modulation frequency

    图  13   探测器灰度均值与激光调制占空比的曲线拟合

    Figure  13.   Curve fitting of detector gray mean value and laser modulation duty ratio

    表  1   激光调制频率变化时的闪烁光斑灰度均值

    Table  1   Gray mean of scintillation spots under changing laser modulation frequency

    Modulation frequency f/Hz Duty cycle δ Gray mean spot
    100 0.1 28.1
    200 0.1 26.0
    300 0.1 25.5
    400 0.1 25.2
    500 0.1 25.0
    下载: 导出CSV

    表  2   激光调制占空比变化时的光斑灰度闪烁均值

    Table  2   Gray mean values of scintillation spots under changing laser modulation duty ratio

    Modulation frequency f/Hz Duty cycle δ Gray mean spot
    100 0.1 28.3
    100 0.2 53.1
    100 0.3 72.7
    100 0.4 88.6
    100 0.5 102.9
    下载: 导出CSV
  • [1]

    Wood R A. Uncooled thermal imaging with monolithic silicon focal planes[J]. Proc. of SPIE on Infrared Technology XIX, 1993, 202: 322-329.

    [2]

    Kruse P. Uncooled thermal imaging, arrays, systems, and applications[J]. Proc. of SPIE, 2001, 3238: 203-211.

    [3] 孔令彬, 易新建, 王典洪, 等. 320×240长波非致冷微测辐射热计红外热象仪的研制[J].光子学报, 2002, 31(5): 596-599.

    KONG L B, YI X J, WANG D H, et al. Development of 320×240 long-wave non-cooling microbolometer infrared thermal imager[J]. Acta Photonica Sinica, 2002, 31(5): 596-599.

    [4]

    Kohin M, Butler N R. Performance limits of uncooled VOx microbolometer focal plane arrays[C]//Proc. of SPIE, 2004, 5406: 447-453.

    [5] 杨风暴.红外物理与技术[M].北京:电子工业出版社, 2014.

    YANG F B. Infrared Physics and Technology[M]. Beijing: Publishing House of Electronics Industry, 2014.

    [6]

    Fieque B, Tissot J L, Trouilleau C, et al. Uncooled microbolometer detector: recent developments at ULIS[J]. Opto-Electronics Review, 2006, 14(3): 187-191.

    [7] 孙恋君.微测辐射热计阵列性能及信号处理电路研究[D].南京: 南京理工大学, 2008.

    SUN L J. Research on Performance and Signal Processing Circuit of Microbolometer Array[D]. Nanjing: Nanjing University of Science and Technology, 2008.

    [8] 陈振林, 孙中泉.半导体制冷器原理与应用[J].微电子技术, 1999(5): 63-65.

    CHEN Z L, SUN Z Q. Principle and application of semiconductor refrigerator[J]. Microelectronic Technology, 1999(5): 63-65.

    [9] 刘华军, 李来风.半导体热电制冷材料的研究进展[J].低温工程, 2004(1): 35-41.

    LIU H Y, LI L F. Research progress of semiconductor thermoelectric refrigeration materials[J]. Low Temperature Engineering, 2004(1): 35-41.

    [10]

    ZHANG Y, Ashe J, LUO Y, et al. Semiconductor lasers and applications[C]//Proc. of SPIE, 2002, 4913: 177-183.

    [11] 孙承纬, 陆启生, 范正修.激光辐照效应[M].北京:国防工业出版社, 2002.

    SUN C W, LU Q S, FAN Z X, et al. Laser Irradiation Effect[M]. Beijing: National Defense Industry Press, 2002.

  • 期刊类型引用(5)

    1. 王怡雯,郑伟,邢成龙. 方向显著性特征分布的加权散度红外弱小目标增强检测. 空间科学学报. 2025(01): 215-225 . 百度学术
    2. 杜妮妮,单凯东,王建超. HRformer:基于多级回归Transformer网络的红外小目标检测. 红外技术. 2024(02): 199-207 . 本站查看
    3. 杜妮妮,单凯东,卫莎莎. LPformer:基于拉普拉斯金字塔多级Transformer的红外小目标检测. 红外技术. 2023(06): 630-638 . 本站查看
    4. 蔺素珍,张海松,禄晓飞,李大威,李毅. RBNSM:一种复杂背景下红外弱小目标检测新方法. 红外技术. 2022(07): 667-675 . 本站查看
    5. 李常宝. 基于滑窗的自适应红外小目标检测方法. 光学技术. 2022(04): 506-512 . 百度学术

    其他类型引用(6)

图(13)  /  表(2)
计量
  • 文章访问数:  2100
  • HTML全文浏览量:  219
  • PDF下载量:  89
  • 被引次数: 11
出版历程
  • 收稿日期:  2019-11-24
  • 修回日期:  2020-10-05
  • 刊出日期:  2021-01-19

目录

/

返回文章
返回