留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

星载遥感短波红外高速高光谱成像仪坏像元识别

刘永征 杜剑 安秦宇宁 杨凡超 张昕 李洪波

刘永征, 杜剑, 安秦宇宁, 杨凡超, 张昕, 李洪波. 星载遥感短波红外高速高光谱成像仪坏像元识别[J]. 红外与激光工程, 2023, 52(2): 20220308. doi: 10.3788/IRLA20220308
引用本文: 刘永征, 杜剑, 安秦宇宁, 杨凡超, 张昕, 李洪波. 星载遥感短波红外高速高光谱成像仪坏像元识别[J]. 红外与激光工程, 2023, 52(2): 20220308. doi: 10.3788/IRLA20220308
Liu Yongzheng, Du Jian, An Qinyuning, Yang Fanchao, Zhang Xin, Li Hongbo. Identification of bad pixels in shortwave infrared high-speed hyperspectral imager for spaceborne remote sensing[J]. Infrared and Laser Engineering, 2023, 52(2): 20220308. doi: 10.3788/IRLA20220308
Citation: Liu Yongzheng, Du Jian, An Qinyuning, Yang Fanchao, Zhang Xin, Li Hongbo. Identification of bad pixels in shortwave infrared high-speed hyperspectral imager for spaceborne remote sensing[J]. Infrared and Laser Engineering, 2023, 52(2): 20220308. doi: 10.3788/IRLA20220308

星载遥感短波红外高速高光谱成像仪坏像元识别

doi: 10.3788/IRLA20220308
详细信息
    作者简介:

    刘永征,男,高级工程师,硕士,主要从事星载高速高光谱成像仪电子学设计方面的研究

  • 中图分类号: TH741

Identification of bad pixels in shortwave infrared high-speed hyperspectral imager for spaceborne remote sensing

  • 摘要: 为改善干涉成像短波红外高速高光谱成像仪的坏像元对复原光谱的影响,利用高光谱成像仪测试流程建立了坏像元识别模板,以提高坏像元识别效率。首先,按照高光谱成像仪测试流程设置增益模板和帧频模板并采集图像数据,依据正常像元增益响应设定合理判定阈值Th1,识别不同增益下异常像元并记录对应坐标值;再依据正常像元帧频响应灰度值设定合理判定阈值Th2,识别不同帧频下异常像元并记录坐标值。最后,对比增益模板和帧频模板判定的异常像元,融合确定坏像元。实验结果表明基于增益模板和帧频模板的识别方法在不增加设备研制测试成本的同时有效识别出短波红外高光谱成像仪探测器的坏像元,为可靠识别短波红外高光谱成像仪坏像元提供了一种低成本、高效可靠的新方法,提高了干涉成像高光谱成像仪光谱反演准确性。
  • 图  1  高分辨率高光谱成像仪原理示意图

    Figure  1.  Schematic diagram of the high-resolution hyperspectral imager

    图  2  短波红外谱段LASIS图像

    Figure  2.  LASIS image in shortwave infrared waveband

    图  3  干涉成像高光谱图像数据的处理流程

    Figure  3.  Function diagram of the procedure of interferometer imaging hyperspectral data

    图  4  模板法检测坏像元技术流程

    Figure  4.  Technology flow chart of detecting bad pixels according our temples

    图  5  增益设计模板

    Figure  5.  Gain design template

    图  6  增益4和增益5下的亮点坏像元

    Figure  6.  Bad pixels in image on Gain4 and Gain5

    图  7  短波红外图像中各像元的帧频响应值

    Figure  7.  The frame rate response value of each pixel in the shortwave infrared image

    图  8  制冷型1000×256大面阵高速短波红外探测器

    Figure  8.  Outline drawing of 1 000×256 refrigerated large area array shortwave infrared detector

    图  9  基于增益模板和帧频模板的坏像元识别流程图

    Figure  9.  Flow chart of bad pixel identification based on gain template and frame rate template

    图  10  高速高光谱成像仪不同增益下亮元及暗元检测

    Figure  10.  Bright pixels and dark pixels recognition for high speed hyperspectral imager at different gains

    图  11  坏像元对大孔径干涉成像短波红外光谱图像干涉图的影响

    Figure  11.  The influence of bad pixels on the short wave infrared spectral image of large aperture interference imaging

    表  1  采用短波红外图像三种方法检测坏像元的附加用时对比

    Table  1.   Comparing additional times of three bad pixels detection methods using in SWIR images

    Image data#Time of M1/sTime of M2/sTime of M3/s
    Data12.2839.57312.796
    Data22.2659.59112.889
    Data32.2369.59912.989
    Data42.2679.57712.799
    Data52.2739.59412.896
    下载: 导出CSV
  • [1] Chen Jianjun, Cui Jicheng, Liu Jianan, et al. Correction technology of HgCdTe short-wave infrared focal plane arrays [J]. Acta Optica Sinica, 2019, 39(2): 0204001. (in Chinese) doi:  10.3788/AOS201939.0204001
    [2] Li Yingying, Zhao Yongchao, Wu Hao, et al. Detection and correction of bad pixels in large aperture static imaging spectrometer [J]. Spacecraft Recovery & Remote Sensing, 2015, 36(5): 76-82. (in Chinese)
    [3] Yu Yi, Chang Songtao, Wang Min, et al. Fast non-uniformity correction for high dynamic infrared radiometric system [J]. Optics and Precision Engineering, 2015, 23(7): 2315-2319. (in Chinese)
    [4] Huang Fuyu, He Yongqiang, Zhou Bing, et al. A new method of bad pixels detection based on SWIR FPA [J]. Infrared Technology, 2010, 32(11): 611-614. (in Chinese)
    [5] Li Qian, Yuan Yan, Cao Jun. Bad pixel correction in temporally and spatially modulated fourier transform imaging spectrometer [J]. Spacecraft Recovery & Remote Sensing, 2016, 37(6): 76-84. (in Chinese)
    [6] Wang Rongbin, Zeng Chao, Jiang Wan, et al. Terra MODIS band 5th stripe noise detection and correction using MAP-based algorithm [J]. Infrared and Laser Engineering, 2013, 42(1): 273-277. (in Chinese)
    [7] Huang Fuyu, He Yongqiang, Zhou Bing, et al. A new method of bad pixels detection based on SWIR FPA [J]. Infrared Technology, 2010, 32(10): 611-614. (in Chinese)
    [8] Hu Mingpeng, Ma Dongmei, Liu Zhixiang, et al. Study of IRFPA non-effective pixel discrimination [J]. Semiconductor Optoelectronics, 2008, 29(2): 184-187. (in Chinese)
    [9] Zhang Jinyue, Lv Junpeng, Ni Zhenhua. Highly sensitive infrared detector based on a two-dimensional heterojunction [J]. Chinese Optics, 2021, 14(1): 87-96. (in Chinese)
    [10] Li Yanjun, Cui Ruiqing, Zhao Guilan, et al. A new method of bad pixel detection and compensation of IRFPA [J]. Laser & Infrared, 2007, 37(1): 51-52. (in Chinese)
    [11] Lai Rui, Liu Shangqian, Zhou Huixin, et al. A new method of bad pixel detection and compensation of IRFPA [J]. Semiconductor Optoelectronics, 2005, 26(3): 51-52. (in Chinese)
    [12] Liang Chao, Ma Tianxiang. Design of infrared imaging nonuniformity correction system based on black body calibration [J]. Chinese Optics, 2016, 9(3): 385-392. (in Chinese)
    [13] Zhang Changxing, Liu Chengyu, Qi Hongxing, et al. Blind pixel detection algorithm using spectral matching for thermal infrared hyperspectral imager [J]. Infrared and Laser Engineering, 2020, 49(1): 0104002. (in Chinese)
    [14] Zhu Shuangshuang, Wu Yang, Zou Peng, et al. Environmental adaptability analysis and performance optimization of space-borne infrared detector assembly [J]. Infrared and Laser Engineering, 2020, 49(2): 0204001. (in Chinese)
    [15] Huang Liang, Li Mingxuan, Lv Hengyi, et al. Large-range and high-resolution temperature measurement system for satellite-borne infrared detector [J]. Optics and Precision Engineering, 2019, 27(11): 2315-2319. (in Chinese)
  • [1] 张前程, 钟胜, 吕劲松, 李显成.  大视场全天时星敏感器光学系统设计 . 红外与激光工程, 2023, 52(3): 20220583-1-20220583-9. doi: 10.3788/IRLA20220583
    [2] 汪鸿祎, 陶文刚, 陆逸凡, 张永刚, 黄松垒, 方家熊.  基于可编程开窗IP核的低功耗读出电路研究 . 红外与激光工程, 2023, 52(12): 20230241-1-20230241-9. doi: 10.3788/IRLA20230241
    [3] 史衍丽, 李云雪, 白容, 刘辰, 叶海峰, 黄润宇, 侯泽鹏, 马旭, 赵伟林, 张家鑫, 王伟, 付全.  短波红外单光子探测器的发展(特邀) . 红外与激光工程, 2023, 52(3): 20220908-1-20220908-16. doi: 10.3788/IRLA20220908
    [4] 曾祥伟, 李亚红, 张燕, 褚金奎.  烟雾环境中532 nm/1550 nm线偏振光前向传输实验研究 . 红外与激光工程, 2022, 51(3): 20210157-1-20210157-6. doi: 10.3788/IRLA20210157
    [5] 马旭, 李云雪, 黄润宇, 叶海峰, 侯泽鹏, 史衍丽.  短波红外探测器的发展与应用(特邀) . 红外与激光工程, 2022, 51(1): 20210897-1-20210897-12. doi: 10.3788/IRLA20210897
    [6] 于春蕾, 龚海梅, 李雪, 黄松垒, 杨波, 朱宪亮, 邵秀梅, 李淘, 顾溢.  2560×2048元短波红外InGaAs焦平面探测器(特邀) . 红外与激光工程, 2022, 51(3): 20210941-1-20210941-10. doi: 10.3788/IRLA20210941
    [7] 张刘, 李博楠, 卢勇男, 邹阳阳, 王泰雷.  基于Offner凸面光栅星载CO2成像光谱仪光学系统设计 . 红外与激光工程, 2022, 51(7): 20220431-1-20220431-9. doi: 10.3788/IRLA20220431
    [8] 赵慧洁, 张晓元, 贾国瑞, 邱显斐, 翟亮.  运动补偿下双通道星载高光谱成像仪图像配准 . 红外与激光工程, 2021, 50(3): 20211022-1-20211022-9. doi: 10.3788/IRLA20211022
    [9] 李春来, 吕刚, 袁立银, 王跃明, 金健, 徐艳, 刘成玉, 何志平, 王建宇.  机载热红外高光谱成像仪的光谱性能测试与初步应用 . 红外与激光工程, 2020, 49(5): 20190117-20190117-7. doi: 10.3788/IRLA20190117
    [10] 张长兴, 刘成玉, 亓洪兴, 张东, 蔡能斌.  热红外高光谱成像仪光谱匹配盲元检测算法 . 红外与激光工程, 2020, 49(1): 0104002-0104002(7). doi: 10.3788/IRLA202049.0104002
    [11] 李雪, 邵秀梅, 李淘, 程吉凤, 黄张成, 黄松垒, 杨波, 顾溢, 马英杰, 龚海梅, 方家熊.  短波红外InGaAs焦平面探测器研究进展 . 红外与激光工程, 2020, 49(1): 0103006-0103006(8). doi: 10.3788/IRLA202049.0103006
    [12] 刘永征, 孔亮, 刘学斌, 陈小来, 刘文龙, 张昕.  星载大孔径静态干涉高光谱成像仪高速成像电路设计 . 红外与激光工程, 2020, 49(S1): 20200096-20200096. doi: 10.3788/IRLA20200096
    [13] 唐绍凡, 鲁之君, 王伟刚, 李欢.  航天高光谱成像仪简述(特邀) . 红外与激光工程, 2019, 48(3): 303003-0303003(9). doi: 10.3788/IRLA201948.0303003
    [14] 吴俊政, 倪维平, 严卫东, 张晗.  圆周阵列太赫兹干涉成像中目标场景仿真 . 红外与激光工程, 2019, 48(1): 125004-0125004(8). doi: 10.3788/IRLA201948.0125004
    [15] 于磊, 陈素娟, 陈结祥, 薛辉.  精准农业观测高数值孔径短波红外成像光谱仪光学系统 . 红外与激光工程, 2018, 47(12): 1218007-1218007(7). doi: 10.3788/IRLA201847.1218007
    [16] 王建宇, 李春来, 王跃明, 吕刚, 袁立银, 金健, 陈小文, 谢峰.  热红外高光谱成像仪的灵敏度模型与系统研制 . 红外与激光工程, 2017, 46(1): 102001-0102001(7). doi: 10.3788/IRLA201746.0102001
    [17] 吴俊政, 严卫东, 倪维平, 张晗.  圆周阵列太赫兹干涉成像仿真 . 红外与激光工程, 2017, 46(8): 825002-0825002(8). doi: 10.3788/IRLA201746.0825002
    [18] 巩盾, 王红.  空间高光谱成像仪的光学设计 . 红外与激光工程, 2014, 43(2): 541-545.
    [19] 李永强, 郭永祥, 王静怡.  基于中阶梯光栅的光谱定标装置设计方法 . 红外与激光工程, 2014, 43(1): 208-211.
    [20] 闵敏, 张勇, 胡秀清, 董立新, 戎志国.  FY-3A中分辨率光谱成像仪红外通道辐射定标的场地评估 . 红外与激光工程, 2012, 41(8): 1995-2001.
  • 加载中
图(11) / 表(1)
计量
  • 文章访问数:  141
  • HTML全文浏览量:  53
  • PDF下载量:  55
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-05
  • 修回日期:  2022-06-16
  • 刊出日期:  2023-02-25

星载遥感短波红外高速高光谱成像仪坏像元识别

doi: 10.3788/IRLA20220308
    作者简介:

    刘永征,男,高级工程师,硕士,主要从事星载高速高光谱成像仪电子学设计方面的研究

  • 中图分类号: TH741

摘要: 为改善干涉成像短波红外高速高光谱成像仪的坏像元对复原光谱的影响,利用高光谱成像仪测试流程建立了坏像元识别模板,以提高坏像元识别效率。首先,按照高光谱成像仪测试流程设置增益模板和帧频模板并采集图像数据,依据正常像元增益响应设定合理判定阈值Th1,识别不同增益下异常像元并记录对应坐标值;再依据正常像元帧频响应灰度值设定合理判定阈值Th2,识别不同帧频下异常像元并记录坐标值。最后,对比增益模板和帧频模板判定的异常像元,融合确定坏像元。实验结果表明基于增益模板和帧频模板的识别方法在不增加设备研制测试成本的同时有效识别出短波红外高光谱成像仪探测器的坏像元,为可靠识别短波红外高光谱成像仪坏像元提供了一种低成本、高效可靠的新方法,提高了干涉成像高光谱成像仪光谱反演准确性。

English Abstract

    • 随着我国高光谱成像遥感持续发展,高分辨率高光谱成像仪愈加显得重要。高分辨率短波红外高光谱成像仪,特别是大孔径静态干涉光谱成像仪(Large Aperture Static Imaging Spectrometry,LASIS)在植物覆盖、海岸线检测、矿石分布统计、自然灾害预警、灾后评估、食品及药品检测、军事敏感目标识别等方面发挥着重要作用。

      碲镉汞(Mercury Cadmium Telluride,HgCdTe)为材料的短波红外探测器具有量子效率高、热激发暗电流小的优点,是LASIS采用的主流器件[1]。但其普遍存在坏像元和响应不均匀的问题,影响了光谱复原准确性,需要进行校正。坏像元识别是校正的前提,国内外关于坏像元识别的研究较多,例如3σ阈值法、图像差分域识别[2]、光谱角β阈值法[3]、滑动窗口法[4]等,但对于LASIS坏像元检测识别及其工程应用的研究极少。

      作者结合多年来的LASIS高光谱成像仪研制经验,提出了适用于LASIS高光谱成像仪的基于模板的坏像元检测方法,多个航天项目也验证了本方法具有方法简捷、效果可靠的优点。

    • 基于大孔径干涉成像光谱技术的高光谱成像仪系统,集成了时间调制和空间调制的优点,克服了高辐射通量与高稳定性之间的矛盾,实现了较高的高分辨高光谱成像性能。LASIS成像数据在目标的二维空间影像上同时叠加了包含其光谱信息的干涉图。基于LASIS的高分辨率高光谱成像仪原理如图1所示,短波红外谱段LASIS获取的实际图像如图2所示,要提取出所关注的光谱信息,数据处理的标准流程如图3所示。

      高光谱图像与普通相机图像对被摄目标关注的侧重点不同,它更关注的是从图像数据所提取出的目标点的光谱信息,故对存在的坏点更为敏感。因此,在提取目标点光谱信息前,必须进行高光谱图像数据的坏像元检测和校正。

      高分辨率高光谱成像仪图像数据量巨大,为确保卫星数传通道的可靠性,高光谱图像数据下传前都进行过图像数据压缩。要利用高光谱图像数据,处理的流程如图3所示,分别是:图像解压缩、坏像元检测、坏像元均匀性校正、抽成大孔径空间调制谱(Large Spatially Modulated Interference Spectral, LSMIS)图像、去除暗电平、干涉图直流去除、切趾、相位修正,傅里叶逆变换,最终获得光谱数据曲线。

      图  1  高分辨率高光谱成像仪原理示意图

      Figure 1.  Schematic diagram of the high-resolution hyperspectral imager

      图  2  短波红外谱段LASIS图像

      Figure 2.  LASIS image in shortwave infrared waveband

      图  3  干涉成像高光谱图像数据的处理流程

      Figure 3.  Function diagram of the procedure of interferometer imaging hyperspectral data

      可以看到,坏像元检测是整个高光谱成像仪图像数据进行应用的最前端步骤,对光谱复原影响最为直接。

      高光谱图像数据从最初的LASIS图像到最终的输出光谱,将所有预处理环节省略掉,数据处理可概括为LASIS图像选取某列干涉曲线到光谱曲线的逆向傅里叶变换[5]

      传统的短波红外成像系统坏像元检测方法是计算整个成像系统均匀光照下的平均响应,将那些响应值小于平均响应值1/10的像元识别为坏像元,并归入死像元类;再计算整个成像系统所有像元的噪声电压,将噪声电压大于10倍的整个探测器像元数平均后的平均噪声电压的像元也识别为坏像元,并归为过热像元类[6]

      目前有研究人员在传统坏像元检测方法的基础上,从应用自动识别算法的角度进行了研究性尝试[7-11],取得了一定效果,但往往仅能在实验室使用,处于科研探索的阶段,并不能应用在成熟稳定的业务化的高光谱图像数据应用中。

      作者根据长期在高光谱成像仪工程项目研发中积累的经验基础上,结合大孔径干涉光谱成像遥感应用流程的特点,提出了一套实用有效的坏像元识别及分类方法,有效解决了多个星载遥感用的干涉成像高光谱成像仪的坏像元识别和分类问题。

      从电子学实现的角度看高光谱图像数据,红外探测器光电转换获取的模拟视频信号,首先经过模拟前端进行信号调理和放大,放大倍数以G=1.4NN= −2、−1、0、1、2)为常规设置,目的是为了匹配不同太阳光高度角照射下及不同反射率时进入星载成像光谱仪光学系统的辐照能量。另外,星载遥感系统往往有变化轨道以适应不同分辨率和幅宽的具体需求,这也要求高光谱成像仪帧频可调以匹配速高比。

      针对航天遥感使用中的变化模拟视频放大倍数(即模拟视频的增益)和变化帧频的使用特点,笔者设计了一种依据设计模板来识别坏像元的方法。操作流程如图4所示。

      根据高光谱成像仪响应辐照亮度范围所对应确定的成像电路需设置的增益级数,设计最小均匀照度(太阳模拟器辐照度设置不使高光谱成像仪饱和)下,依次增大各级增益时的成像电路响应的图像,计算各档增益下1000×256面阵各像元相应响应的平均值。

      确认是否在模板1划定范围内?对不在划定范围内的灰度值进行标记其对应像元,并进行分类。

      根据高光谱成像仪需要响应不同轨道高度下所应摄影的辐照亮度范围所对应确定的成像电路的帧频设置,设计最小帧频时,确定辐射亮度下,依次增大帧频时的成像电路响应的图像,计算各级帧频下1000×256面阵各像元相应给定辐照亮度响应的平均值。

      确认是否在模板2划定范围内?对不在划定范围内的灰度值进行标记其对应的像元,并进行分类。

      图  4  模板法检测坏像元技术流程

      Figure 4.  Technology flow chart of detecting bad pixels according our temples

    • 航天遥感为获取高的空间分辨率和宽的对地遥感幅宽,所选用短波红外探测器图像面阵为1000×256面阵,像元尺寸为30 μm×30 μm,最高帧频为390 fps。

      高光谱成像仪设计的5档增益中,对应的曝光时间控制范围从2.6~5.6 ms。在星载系统中,卫星平台根据观测目标,及时修正高光谱成像仪的增益,故此,增益设置指令模块是发射前例行且必须检查的检测和演练项目。设计增益模板来对高光谱成像仪短波红外谱段的坏像元进行识别和分类,没有增加高光谱成像仪研制过程的额外工作量和检查项目及流程。其过程为依次设定高光谱成像仪工作在5档增益下,并相应记录对应的图像数据;计算各档增益下的图像多帧平均灰度值,作为该增益下对应的正常响应;根据各档增益下图像响应值,设计增益响应模板,作为界定非正常响应的灰度响应界限;按照模板,逐档增益判读异常像元。将响应值未落入增益模板内的像元,记录其位置,识别为坏像元。

      图5中以最简单的增益模板,高光谱成像仪设计为从增益1依次增大到最高档增益,图像响应值从平均值989~3800,在此过程中,计算出正常图像平均响应服从响应曲线L6,考虑测量误差及图像面阵的不一致性容限,模板刻画为曲线L5 (相对L6,DN值低30)和L7 (相对L6,DN值高30)内为正常响应像元(L5L7容限值10 DN);按照增益1到增益5增大响应时,实际响应未落入模板内的即为异常像元,标记异常像元的响应曲线,并按其特征进行分类。

      图  5  增益设计模板

      Figure 5.  Gain design template

      图  6  增益4和增益5下的亮点坏像元

      Figure 6.  Bad pixels in image on Gain4 and Gain5

      高光谱成像仪的短波红外图像中,出现1类异常像元(简称第1类)为响应在各档增益下整体偏低,且其灰度值不随增益增大而增大(固定为暗背景值),其响应特征曲线为L1,定义为暗点像元,简称暗点;图像中还出现了1类像元(记作第2类),随增益增大时整体响应偏低,但灰度值在增益增大时会略微增大,其响应特征曲线为L2,称之为弱响应像元;图像中也大量出现了一类像元(记为第3类),其表现为在前三档增益下,其响应值高于正常像元响应值,在后两档增益下,其响应值高于正常响应值,响应特征曲线表现为L3,称之为非线性响应像元;图像中第4类像元,其在各档增益下,灰度值均很高,且不随积分时间增大而变化,其响应特征曲线表现为L4,现象比较典型,如图6所示,简单描述为常亮,称之为亮点。绝大多数像元响应灰度都值服从随着增益增大而增大,响应特征曲线可以被L5L6L7三条曲线所表征,这些都是正常像元。

      以往坏点去除研究的不足主要是坏点的漏判和多判比较严重,漏判是指对探测器坏点识别不全面,多判是指将探测器的部分非坏像元识别为坏像元,导致出现误判和漏判现象的主要原因是坏像元检测算法不适用。

    • 星载高光谱成像仪要满足在轨变轨条件下任务需求,需测试不同帧频下的成像能力,故帧频响应能力是必测项目,作者充分利用装星及发射前必测项目设计了帧频检测模板来识别坏像元。

      在星载遥感高速高光谱成像仪中,设计短波红外高光谱成像仪在1000×256图像幅面下依次工作帧频为230、260、290、320、350、380 fps,采集并记录连续帧图像数据;计算连续帧图像中各个像元的不同帧频时响应灰度值。

      在此设计模板下,短波红外光谱通道(0.8~3.0 μm)图像数据量高达1000×256×12 bit×390 fps =1.14 Gbps。在每个工作帧频下,计算每个像元灰度值是指计算该帧频下采集到图像的每个像元的灰度响应。以帧频数值为横坐标,以灰度值为纵坐标,对所有像元进行画图,即可得到成像仪各像元的帧频响应分布,如图7所示。

      图  7  短波红外图像中各像元的帧频响应值

      Figure 7.  The frame rate response value of each pixel in the shortwave infrared image

    • 碲镉汞短波红外焦面探测器是干涉成像高光谱成像仪的核心器件,技术水平与制作工艺的限制使得HgCdTe短波红外焦面探测器仍普遍存在坏像元和感应非均匀性的像元,且坏像元的数量会随着探测器开关机次数的增加而增加,因此对HgCdTe短波红外焦面探测器进行坏像元识别非常必要。

      通过前述的两种识别模板,对HgCdTe短波红外焦面探测器的所有像元遍历模板涉及的所有测试点响应后,可有效识别出所有坏像元,通过对其进行分类,可以为采取有效的坏像元校正方法提供有效的参考依据。

      通过对碲镉汞短波红外焦面探测器的异常响应像元的响应特征进行分析,可以将碲镉汞短波红外焦面探测器的坏像元分为以下几类:

      (1)暗点像元,其对两类检测模板的增益和帧频变化时的输入光照能量变化没有响应上的变化,图像上表现为常暗,灰度值在暗背景灰度值附近,图像灰度值响应曲线如模板1(增益模板)和模板2(帧频模板)中的L1曲线。

      (2)弱响应像元,其对两类检测模板的增益和帧频变化时的输入光照能量变化有响应,但响应均较大多数正常像元的灰度值偏暗,灰度值总是低于大多数像元,图像灰度值响应曲线如模板1和模板2中的L2曲线。

      (3)非线性响应像元,其对两类检测模板中的增益和帧频变化时输入的能量变化有响应,且其比大多数正常像元的响应差异小。只有仔细比对完整个响应范围才能发现其与大多数线性响应的像元有响应特征上的差异,表现为输入能量不强时,其响应比正常像元的响应值高;输入能量强时,其响应却比正常的像元响应值低,其响应特征曲线如模板1和模板2中的曲线L3所示。

      (4)亮点像元,其对两类检测模板中的增益和帧频变化没有响应上的变化,上电即常亮。其响应特征表现为模板1和模板2中的曲线L4所示。

      HgCdTe短波红外焦面探测器的多数像元均在模板1和模板2的测试中表现出线性变化,其在图像上灰度值的特征表现为曲线L6所示。因为像元响应值的分布误差,会表现为部分像元响应特征在L6附近,以曲线L5为响应下限边缘,曲线L7为响应上限边缘,也就是说大多数在L5L7划定边界的模板内。

      坏像元修正和非均匀性修正都是红外探测器投入高光谱遥感应用前必须的修正过程,坏像元识别和分类是进行修正的前提。通过坏像元校正将坏像元的错误输出值校正为合理的输出值,通过非均匀性校正来校正像元的响应曲线,使他们对同样强度的红外辐射具有趋于正常像元的响应值。

    • 在高光谱成像仪电子学设计技术流程中,筛选短波红外探测器阶段在组装之前,对短波红外高光谱成像仪的坏像元进行识别后,然后按其特征进行分类,这样为有效的修正提供了便利。短波红外高光谱成像仪采用定制大面阵制冷型探测器(外观如图8所示)。

      图  8  制冷型1000×256大面阵高速短波红外探测器

      Figure 8.  Outline drawing of 1 000×256 refrigerated large area array shortwave infrared detector

      按照短波红外高光谱成像仪的设计工作模式,测试增益模板中其不同增益(G表示Gain,分别取值 G1G2G3G4G5)下的图像,统计面阵图像中的坏像元数量;再按照帧频模板中工作模式,测试帧频变化范围内图像的响应值,特别保留选取模板帧频为180、220、260、300、340、380 fps的图像,复核增益模板。对两次模板统计出的坏像元,按类分别筛查其面阵坐标位置,生成坏像元坐标文件。图9是设计模板检测坏像元的流程图。

      暗背景下采集图像,识别亮点像元如图10所示。

      通过增益模板的检测,可检测出短波红外高光谱成像仪图像中亮点像元、暗点像元、弱响应像元、非线性响应像元[12]

      通过对星载干涉成像短波红外高光谱成像仪进行亮点及暗点像元检测,在12 bit量化下,在增益4下,图像响应灰度值在2500附近时,检测结果最为准确。

      图  9  基于增益模板和帧频模板的坏像元识别流程图

      Figure 9.  Flow chart of bad pixel identification based on gain template and frame rate template

      一般用来表示坏像元检测准确率的参数有两个,分别是:漏检测系数(Detection Missed Factor):定义为漏检坏像元数占比实际坏像元数的百分比。虚假检测系数(Detection False Factor):定义为虚假检测坏像元数占比实际坏像元数的百分比。

      $$ {{{N}}}_{\text {miss }}=\left(1-\frac{b}{d+h+w}\right) {\text{%}} $$ (1)

      式中:b代表所采用方法检测到的坏像元;d表示实际没有任何响应的暗点像元;h表示响应值过高的亮点像元;w表示有响应但响应辐照不随正常像元随曝光时间增加呈现线性响应的像元。

      $$ {N}_{\mathrm{err}}=\left(\frac{a}{M \times N-d-h-w}\right) {\text{%}} $$ (2)

      式中:a代表被所用方法检测出的坏像元;d表示实际没有任何响应的暗点像元;h表示响应值过高的亮点像元;w表示对光照有响应但响应辐照不随正常帧长像元呈现线性响应的像元;分母中M表示图像面阵行数, N表示图像列数[13],在文中项目中M=256,N=1000。

      图  10  高速高光谱成像仪不同增益下亮元及暗元检测

      Figure 10.  Bright pixels and dark pixels recognition for high speed hyperspectral imager at different gains

      为了比对本方法的便捷性,通过对比同一图像数据处理工程师编写实现的三个方法在对相同的5组图像数据进行坏像元识别时所附加消耗的时长(见表1,运算用计算处理平台相同),考核处理防范的便捷性。方法1:记作M1,代表文中提出的模板法;方法2:记作M2,同一软件设计师参考标准流程编写的“滑动窗口法”;方法3:记作M3,同一设计师参考标准流程编写“图像差分域识别法”。附加时间指的是,软件不识别短波红外图像用时与读入并遍历显示一遍对应图像数据进行坏像元识别所额外占用的附加时间,所增加的额外附加时间。

      表 1  采用短波红外图像三种方法检测坏像元的附加用时对比

      Table 1.  Comparing additional times of three bad pixels detection methods using in SWIR images

      Image data#Time of M1/sTime of M2/sTime of M3/s
      Data12.2839.57312.796
      Data22.2659.59112.889
      Data32.2369.59912.989
      Data42.2679.57712.799
      Data52.2739.59412.896
    • 对按照设定模板识别并分类的坏像元中的暗点像元和亮点像元进行了修正,以验证坏像元识别及分类的有效性。

      按照之前筛选检测生成的坏像元坐标文件,对图像中暗点像元和亮点像元进行了修正,修正方法选取工程中常用的拟合插值替代法[14-15]。大孔径干涉成像高光谱图像在修正坏像元时沿干涉曲线方向(列)进行拟合插值替代,采用分段二次函数拟合的方法。

      修正前各列典型的坏像元对本列干涉图曲线影响见图11 (a)~(c),修正后的LASIS图像结果见图11 (b)。图11 (a)是增益3下识别的第486列的两个坏像元(标记为1#、2#),(c)是第605列的坏像元修正后的列干涉图。测试结果及工程应用表明:文中提出的坏像元检测方法,具有实用性强、工程开支成本低、方法便捷等优点。

      图  11  坏像元对大孔径干涉成像短波红外光谱图像干涉图的影响

      Figure 11.  The influence of bad pixels on the short wave infrared spectral image of large aperture interference imaging

    • 大面阵制冷型碲镉汞短波红外焦面探测器在星载遥感干涉成像高光谱成像仪中的广泛应用,使得其坏像元对光谱复原的影响问题变得日益显著。文中结合大孔径干涉成像高速高光谱成像仪研制的工程化流程,提出了一种新的低成本的增益模板+帧频模板的坏像元检测方法,经过广泛验证,结果表明该方法效率高成本低,为有效矫正坏像元提供了便利,使得从干涉图到光谱数据立方体的提取更高效准确。目前作者已在具体遥感高光谱成像项目中实际使用该方法,方法的工程化应用大大提升了短波红外干涉成像高光谱成像仪遥感业务的可靠性,为遥感用干涉成像高光谱成像仪的规范化提供了参考。

参考文献 (15)

目录

    /

    返回文章
    返回