不同高度不同大气模式下红外图像仿真

李金萍, 蔡友林, 刘其涛

李金萍, 蔡友林, 刘其涛. 不同高度不同大气模式下红外图像仿真[J]. 红外技术, 2024, 46(10): 1218-1223.
引用本文: 李金萍, 蔡友林, 刘其涛. 不同高度不同大气模式下红外图像仿真[J]. 红外技术, 2024, 46(10): 1218-1223.
LI Jinping, CAI Youlin, LIU Qitao. Simulation of Infrared Images at Different Heights and Atmospheric Modes[J]. Infrared Technology , 2024, 46(10): 1218-1223.
Citation: LI Jinping, CAI Youlin, LIU Qitao. Simulation of Infrared Images at Different Heights and Atmospheric Modes[J]. Infrared Technology , 2024, 46(10): 1218-1223.

不同高度不同大气模式下红外图像仿真

基金项目: 

江西省放射性地学大数据技术工程实验室开放基金 JELRGBDT202103

江西省核地学数据科学与系统工程技术研究中心开放基金 JETRCNGDSS201806

详细信息
    作者简介:

    李金萍(1972-),女,副教授,硕士,主要从事卫星遥感图像处理,红外图像仿真,目标特征提取与识别,数字图像处理算法研究等。E-mail: 379896761@qq.com

  • 中图分类号: TP721

Simulation of Infrared Images at Different Heights and Atmospheric Modes

  • 摘要:

    本文以飞机平台的红外图像为数据源,研究了不同大气模式下成像高度对红外仿真图像的影响,推导了随热像仪成像高度变化的辐射亮度计算公式。分别采用MODTRAN和LOWTRAN两种大气辐射传输模型进行不同高度的红外图像大气修正,并对所得仿真数据进行了分析比较得出:在高度低于20 km时,同种大气模式下两种模型仿真数据接近,误差不大,都很理想;在高于20 km时,LOWTRAN模型仿真数据参考价值不大,MODTRAN模型结果更接近实际情况、更实用。

    Abstract:

    The atmosphere and imaging height directly affect the infrared image quality of thermal imagers. Using the infrared image of an aircraft platform as the data source, infrared simulation images at different heights and atmospheric modes were studied, and a simulation calculation formula for the height of the thermal imager was derived. By selecting different atmospheric modes, examining MODTRAN and LOWTRAN atmospheric radiation transmission models for infrared image atmospheric correction at different heights, and analyzing and comparing the obtained simulation data, it is concluded that when the distance is less than 20 km, the simulation data of the two models in the same atmospheric mode are close and that the error is small, which is ideal. When the distance was greater than 20 km, the reference value of the LOWTRAN model simulation data was small, and the results of the MODTRAN model were closer to the actual situation and therefore more practical.

  • 碲镉汞(mercury cadmium telluride,MCT)薄膜具有带隙在整个红外波段连续可调、高吸收系数、高量子效率、低暗电流等优点,碲镉汞红外探测器技术成为制冷型高性能红外探测器发展近60年来最主流,也是工程化最成功的技术路线[1-2]。目前红外探测器正向着第三代探测器工程化及第四代探测器机理探索的方向发展,主要涵盖大面阵、双/多色、近室温工作温度、雪崩光电二极管(avalanche photon diode,APD)、集成化、智能化以及多维探测等方向[3]。碲锌镉(cadmium zinc telluride,CZT)衬底与碲镉汞薄膜有着最小的晶格失配,相比于其他衬底,可实现更低缺陷密度、更高晶体质量碲镉汞薄膜材料的制备,为开发低暗电流、高量子效率红外焦平面探测器提供了有利保障。

    超大尺寸碲镉汞薄膜制备是亿像素规模红外探测器的研制基础,提升宇航等领域对空间目标和对地目标的大视场、超高分辨率监测能力,同时实现千万像素高性能红外探测器的批产,满足军用和民用领域的批量装备需求。近年来,红外探测器各主流机构在该领域投入了大量时间和精力,推动了大尺寸碲镉汞薄膜材料技术的快速发展,如图 1所示[4-6]。国外方面,美国Raytheon Vision Systems公司采用分子束外延(molecular beam epitaxy,MBE)技术已成功制备了80 mm×80 mm CZT衬底基MCT薄膜[7];Teledyne Image sensors公司也采用MBE技术实现了70 mm×75 mm CZT衬底基MCT薄膜的生长,厚度标准差为0.06 μm,组分标准差为0.0004,位错腐蚀坑密度(etch pit density, EPD)<5×105 cm-2 [8];法国Lynred(Sofradir)采用液相外延(liquid phase epitaxy, LPE)技术实现了72 mm×73 mm(CZT衬底)MCT薄膜的外延,厚度高差为0.12 μm,波长差为0.02 μm(50 mm×50 mm)X射线双峰回摆曲线半峰宽(full width at half maxima, FWHM)为32 arcsecs[9];美国BAE具备了60 mm×80 mm及以下尺寸CZT基MCT薄膜的批量生产能力[10]。国内方面,中电11所报道了80 mm×80 mm MCT材料技术,EPD<4×104 cm-2;上海技物所采用富碲垂直液相外延,在CZT衬底上外延50 mm×50 mm MCT薄膜,组分均方差为0.0004,厚度均方差为0.4 μm,FWHM<30 arcsecs,EPD<1×105 cm-2 [11];高德红外公司采用富碲垂直液相外延技术将MCT材料尺寸拓展到了40 mm×50 mm,组分均匀性为0.00065[12]

    图  1  国内外研发碲镉汞薄膜进展报道[4-12]
    Figure  1.  Research progress on the development of HgCdTe thin films[4-12]

    自2019年以来,昆明物理研究所在大尺寸MCT外延材料方面不断取得突破,陆续实现了40 mm×50 mm、50 mm×60 mm以及70 mm×75 mm的碲镉汞材料的研制及批产能力建设[8]。在下一代更高性能红外探测器的研制及装备应用需求牵引下,昆明物理研究所实现了高质量100 mm×100 mm大尺寸短波、中波MCT材料的制备。

    CZT与MCT具有相同的闪锌矿结构,通过调节其Zn组分可实现与MCT材料的完全晶格匹配,成为高质量MCT薄膜研制的首选衬底,有利支撑了高性能红外探测器的研发和生产。然而,实现高质量大尺寸CZT晶锭生长的难度极大,因存在层错能低、临界切应力小、热导率低等情况,导致CZT晶锭中极易出现多晶、孪晶/层错、位错等缺陷;同时由于生长过程中的Zn组分分凝效应,CZT晶体轴向Zn组分分布不均匀,只有实现定向生长并通过垂直轴向切割才能实现高Zn组分均匀性大尺寸衬底的制备。

    昆明物理研究所为实现大尺寸、Zn组份均匀、低缺陷密度CZT衬底的制备,采用籽晶质量提升、温场优化等技术,基于垂直梯度凝固法实现ϕ150 mm CZT单晶的定向生长,如图 2(a)所示。在晶体生长过程中,通过组份偏析抑制等技术实现缺陷密度的有效抑制,生长的晶体EPD≤4×104 cm-2,沉积相尺寸小于5 μm、密度小于5×103 cm-2图 2(b))。上述制备的ϕ150 mm碲锌镉晶锭经切割、磨抛等工艺后,获得了100 mm×100 mm大尺寸高质量外延级碲锌镉衬底(图 2(c)),Zn组分分布极差达到±0.2%(图 3(a))。通过化学腐蚀的方法,在衬底磨抛工艺中始终保证应力平衡状态,同时完全消除损伤层残余应力对面形的影响,成功使得100 mm×100 mm碲锌镉衬底的峰谷(peak to valley, PV)值达到10 μm,如图 3(b)所示,在后续亿像素红外探测器研制中降低了器件制备、芯片与读出电路互连工艺难度。

    图  2  100 mm×100 mm碲锌镉衬底研制情况:(a)ϕ150 mm碲锌镉晶锭;(b) ϕ150 mm碲锌镉晶锭局部沉积相分布;(c)100 mm×100 mm碲锌镉衬底
    Figure  2.  The development status of 100 mm×100 mm CZT substrate: (a) ϕ150 mm cadmium zinc telluride (CZT) ingots; (b)Typical precipitate density of ϕ150 mm CZT ingots; (c) 100 mm×100 mm CZT substrate
    图  3  100 mm×100 mm碲锌镉衬底Zn组分分布情况与((a))100 mm×100 mm碲锌镉衬底PV值((b))
    Figure  3.  Zn mole fraction distribution over a 100 mm×100 mm CZT substrate((a)) and peak to valley of a 100 mm×100 mm CZT substrate((b))

    大面积水平液相外延技术制备MCT薄膜的技术难点主要有两个。首先是薄膜厚度和组分均匀性控制难度大,可能会阻碍高响应均匀性红外探测器组件的研制。随着材料尺寸的增加,母液在降温生长过程中由于热传导速率的差异会产生更大的温度梯度,薄膜面内生长速率不均匀,从而引起薄膜横向厚度和组分的不均匀;严重时可能会直接产生表面波纹,影响探测器响应均匀性。

    其次是大面积材料缺陷控制,MCT薄膜缺陷类型主要包含位错缺陷和表面缺陷。位错缺陷主要和衬底有关,通过有效抑制衬底本身的位错密度,并筛选合适的衬底Zn组分并优化外延生长程序,MCT薄膜的位错密度能够得到有效抑制。表面缺陷主要包含两种,一种与衬底相关的表面缺陷,另一种是由于外延石墨舟结构不合理或温场均匀性、降温速率等参数不合适导致的,具体形成机理及抑制方法在其他文献中已经进行阐述[2, 8]

    昆明物理研究所针对以上问题,对石墨舟结构进行了重新设计,优化了石墨舟结构,提高散热及温场均匀性;另外,通过对设备温度补偿进行调节,保证生长过程中母液处于一个优异的恒温区间,减小母液自身的温度梯度。在此基础上,采用富碲水平推舟液相外延技术成功实现了100 mm×100 mm MCT材料的生长。图 4(a)是在CZT衬底(111)B面上外延生长所得的MCT材料,该材料表面光亮、核心区域无母液残留及其他宏观材料缺陷,材料表面质量优异。采用Olympus 3D数码显微镜对该薄膜进行了表征,在整个100 mm×100 mm范围内的薄膜表面未出现表面宏观特征波纹,同时材料表面无明显微观缺陷(图 4(b))。

    图  4  100 mm×100 mm碲镉汞薄膜:(a) 照片;(b) 微观形貌
    Figure  4.  100 mm×100 mm MCT films: (a) Photo and (b) Microstructure morphology

    使用XRD对薄膜晶体衍射半峰宽(FWHM)进行了测试,薄膜半峰宽均值≤33 arcsecs(图 5(a))。使用Chen腐蚀液对材料进行了腐蚀并进行了EPD测试,密度在5×104 cm-2以下(图 5(b))。

    图  5  100 mm×100 mm碲镉汞薄膜半峰宽和EPD:(a) 薄膜半峰宽分布(每一个半峰宽测试点区域为1 mm×1 mm);(b) 薄膜局部EPD测试结果
    Figure  5.  FWHM distribution and typical EPD of the 100 mm×100 mm MCT films: (a) FWHM distribution(the area of each test point is 1 mm×1 mm); (b) Typical EPD of the 100 mm×100 mm MCT films

    用傅里叶红外透过显微镜对室温薄膜组分进行了测试,测试结果如图 6所示。在100 mm×100 mm范围内,材料厚度分布呈现中间薄、四周偏高的盆地型分布。短波碲镉汞薄膜厚度均值为6.1 μm,波长均值为2.91 μm,厚度极差±1.24 μm,波长极差±0.025 μm,组分极差为±0.0022;中波薄膜的厚度均值为7.9 μm,波长均值为4.22 μm,厚度极差±1.25 μm,波长极差±0.06 μm,组分极差为±0.0031。

    图  6  100 mm×100 mm碲镉汞薄膜波长厚度分布情况:(a) 短波碲镉汞薄膜厚度分布;(b) 短波碲镉汞薄膜波长分布;(c) 中波碲镉汞薄膜厚度分布;(d) 中波碲镉汞薄膜波长分布
    Figure  6.  Thickness and wavelength distribution of 100 mm×100 mm infrared MCT films: (a) Thickness distribution of the short wavelength infrared MCT films; (b) Wavelength distribution of the short wavelength infrared MCT films; (c) Thickness distribution of the middle wavelength infrared MCT films; (d) Wavelength distribution of the middle wavelength infrared MCT films

    为验证100 mm×100 mm碲镉汞薄膜性能,采用昆明物理研究所成熟的1280×1024@15 μm中波器件工艺进行流片,在单片100 mm×100 mm薄膜上一次性制备20个1280×1024@15 μm碲镉汞中波红外焦平面芯片,大幅提升流片效率,管芯排布如图 7(a)所示。制备的典型器件性能良好,该器件响应图及盲元分布图如图 7(b)(c)所示,器件总盲元数仅275个。

    图  7  100 mm×100 mm碲镉汞薄膜中波1280×1024@15 μm流片芯片排布情况及红外焦平面性能测试结果:(a) 一次流片碲镉汞焦平面芯片排布;(b) 1280×1024(15 μm)规格焦平面探测器信号响应图(77 K);(c) 1280×1024(15 μm)规格焦平面探测器盲元分布(77 K)
    Figure  7.  Chip arrangement and performance of 1280×1024@15 μm MWIR focal plane arrays of (a)Chip arrangement on one 100 mm×100 mm MCT film; (b) Thermal response (77 K) of 1280×1024 (15 μm) focal plane array; (c) Spatial pixel operability map (77 K) of 1280×1024 (15 μm) focal plane array

    通过改进型垂直梯度凝固法,突破了ϕ150 mm碲锌镉晶锭定向生长,实现了100 mm×100 mm碲锌镉衬底的制备,EPD≤4×104 cm-2,沉积相尺寸小于5 μm、密度小于5×103 cm-2。采用富碲水平推舟液相外延技术实现了100 mm×100 mm大尺寸短波、中波碲镉汞材料的制备,其中短波碲镉汞薄膜厚度极差±1.24 μm,波长极差±0.025 μm,组分极差±0.0022;中波薄膜厚度极差±1.25 μm,波长极差±0.06 μm,组分极差±0.0031,是目前国际上报道的尺寸最大的碲锌镉基碲镉汞薄膜材料。在单片100 mm×100 mm薄膜上实现了20个1280×1024@15 μm碲镉汞中波红外焦平面芯片制备,芯片性能良好,为后续亿像素级红外焦平面探测器的研制奠定了基础。

  • 图  1   MODTRAN大气传输模型

    Figure  1.   MODTRAN atmospheric transport model

    图  2   LOWTRAN大气传输模型

    Figure  2.   LOWTRAN atmospheric transport model

    图  3   仿真图像比较

    Figure  3.   Simulation image comparison

    表  1   20 km高度的平台仿真结果

    Table  1   Simulation results of 20 km platform

    Airborne image radiation luminance data/(W·cm-2·sr-2) Simulation platform image radiation brightness data /(W·cm-2·sr-2)
    Mid-latitude winter Tropical atmosphere
    (MODTRAN) (LOWTRAN) (MODTRAN) (LOWTRAN)
    33.62 28.92 28.35 25.8 25.41
    33.64 28.93 28.37 25.83 25.43
    33.62 28.92 28.35 25.80 25.41
    33.64 28.93 28.37 25.83 25.43
    33.95 29.21 28.62 26.09 25.68
    33.9 29.15 28.58 26.04 25.54
    33.93 29.2 28.6 26.07 25.67
    33.95 29.21 28.62 26.09 25.68
    33.06 28.61 28.33 25.52 25.46
    33.04 28.59 28.31 25.47 25.43
    33.11 28.69 28.37 25.64 25.6
    33.12 28.71 28.39 25.65 25.62
    Relative error 1.67% 1.12%
    下载: 导出CSV

    表  2   100 km高度的平台仿真结果

    Table  2   Simulation results of 100 km platform

    Airborne image radiation luminance data/(W·cm-2·sr-2) Simulation platform image radiation brightness data /(W·cm-2·sr-2)
    Mid-latitude winter Tropical atmosphere
    (MODTRAN) (LOWTRAN) (MODTRAN) (LOWTRAN)
    33.62 25.26 28.33 22.81 25.4
    33.64 25.61 28.36 22.83 25.43
    33.62 25.26 28.33 22.81 25.4
    33.64 25.61 28.36 22.83 25.43
    33.95 26.58 28.59 23.09 25.66
    33.9 26.55 28.54 23.04 25.54
    33.93 26.57 28.57 23.07 25.65
    33.95 26.58 28.59 23.09 25.66
    33.06 24.67 28.32 22.18 25.46
    33.04 24.61 28.30 22.16 25.43
    33.11 24.70 28.38 22.34 25.61
    33.12 24.73 28.41 22.35 25.63
    Relative error 10.10% 11.00%
    下载: 导出CSV
  • [1] 肖甫, 吴慧中, 肖亮, 等. 地面坦克目标红外热成像物理模型研究[J]. 系统仿真学报, 2005(11): 2577-2579, 2585. DOI: 10.3969/j.issn.1004-731X.2005.11.002

    XIAO Fu, WU Huizhong, XIAO Liang, et al. Research on infrared imaging model of land tank target[J]. Journal of System Simulation, 2005(11): 2577-2579, 2585. DOI: 10.3969/j.issn.1004-731X.2005.11.002

    [2] 姚涛, 李一凡. 场景红外成像仿真原理和应用[J]. 计算机仿真, 2004, 21(1): 96-98, 135. DOI: 10.3969/j.issn.1006-9348.2004.01.035

    YAO Tao, LI Yifan. Principle and application of infrared imaging simulation for scene[J]. Computer Simulation, 2004, 21(1): 96-98, 135. DOI: 10.3969/j.issn.1006-9348.2004.01.035

    [3] 韩涛, 朱学光. 一种计算大气透过率的方法[J]. 红外技术, 2002, 24(6): 51-53. DOI: 10.3969/j.issn.1001-8891.2002.06.012

    HAN Tao, ZHU Xueguang. A method of calculating the atmosphere transmissibility[J]. Infrared Technology, 2002, 24(6): 51-53. DOI: 10.3969/j.issn.1001-8891.2002.06.012

    [4] 王兵学, 张启衡, 陈昌彬, 等. 凝视型红外捜索跟踪系统的作用距离模型[J]. 光电工程, 2004, 31(7): 8-11.

    WANG Bingxue, ZHANG Qiheng, CHEN Changbin, et al. A mathematical model for operating range of a staring IR search and track system[J]. Opto-Electronic Engineering, 2004, 31(7): 8-11.

    [5] 贾庆莲, 乔彦峰, 邓文渊. 周视搜索系统对点目标的作用距离分析[J]. 光学学报, 2009, 29(4): 937-943.

    JIA Qinglian, QIAO Yanfeng, DENG Wenyuan. Analysis for point- target detection range of panoramic searching system[J]. Acta Optica Sinica, 2009, 29(4): 937-943.

    [6] 张乐, 梁冬明, 姚梅, 等. 红外搜索跟踪系统作用距离等效折算[J]. 红外与激光工程, 2013, 42(1): 26-30.

    ZHANG Le, LIANG Dongming, YAO Mei, et al. Equivalent calculation of operating range of IRST[J]. Infrared and Laser Engineering. 2013, 42(1): 26-30.

    [7]

    ZHAO Y, WU P, SUN W. Calculation of infrared system operating distance by spectral bisection method[J]. Infrared Physics & Technology, 2014, 63(3): 198-203.

    [8] 江飞虹, 吴文贤, 石岩. 机载热红外图像特征的提取与分析[J]. 烟台大学学报(自然科学与工程版版), 2007(4): 295-298.

    JIANG Feihong, WU Wenxian, SHI Yan. Features extraction and analysis from airborne thermal infrared image[J]. Journal of Yantai University (Natural Science and Engineering Edition), 2007(4): 295-298.

    [9] 李金萍. 机载红外图像转换生成卫星仿真图像的大气修正方法研究[J]. 激光与红外, 2023, 53(1): 75-78. DOI: 10.3969/j.issn.1001-5078.2023.01.011

    LI Jinping. Research on atmospheric correction method of airborne infrared image conversion to satellite simulation image[J]. Laser & Infrared, 2023, 53(1): 75-78. DOI: 10.3969/j.issn.1001-5078.2023.01.011

    [10] 刘其涛. 大气对红外热像仪成像仿真的影响[J]. 红外, 2006(1): 1-4.

    LIU Qitao. Influence of atmosphere on simulation of infrared thermal-imager imaging[J]. Infrared, 2006(1): 1-4.

    [11] 刘美玲, 明冬萍. 遥感地学应用实验教程[M]. 北京: 科学出版社, 2018.

    LIU Meiling, MING Dongping. Experimental Course on Remote Sensing Geoscience Applications[M]. Beijing: Science Press, 2018.

    [12] 田杰. 目标红外辐射特征与建模仿真技术研究[D]. 南京: 南京理工大学, 2016.

    TIAN Jie. Research on Target Infrared Radiation Characteristics and Modeling and Simulation Technology[D]. Nanjing: Nanjing University of Science and Technology, 2016.

    [13] 孟凡斌, 郑丽. 基于LOWTRAN7的红外大气透过率计算方法[J]. 光电技术应用, 2009, 24(3): 29-32, 66. DOI: 10.3969/j.issn.1673-1255.2009.03.008

    MENG Fanbin, ZHENG Li. LOWTRAN7-based calculation method of IR transmittance in the atmosphere[J]. Electro-Optic Technology Application, 2009, 24(3): 29-32, 66. DOI: 10.3969/j.issn.1673-1255.2009.03.008

图(3)  /  表(2)
计量
  • 文章访问数:  87
  • HTML全文浏览量:  16
  • PDF下载量:  39
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-10-29
  • 修回日期:  2023-11-25
  • 刊出日期:  2024-10-19

目录

/

返回文章
返回