红外高光谱大气探测仪(FY-3E/HIRAS-II)在轨数据非线性校正方法【增强内容出版】
1 引言
风云三号E星(FY-3E)是全球首颗民用晨昏极轨气象卫星,于2021年7月5日在酒泉卫星发射中心成功发射。风云三号E星经过在轨测试和校准后投入使用,与风云三号上午星C星和风云三号下午星D星组网运行,可每6小时为数值天气预报提供一次全球覆盖的卫星资料,有效补充了6小时同化时间窗内卫星观测资料的空白。[1]FY-3E/HIRAS-II(High-spectral Infrared Atmospheric Sounder II)是我国风云系列极轨气象卫星中所搭载的第二个红外高光谱大气探测仪,具有高光谱分辨率、高灵敏度和高精度等特点,在大气温湿度廓线反演、数值天气预报、气候变化研究和大气痕量气体探测等应用上具有显著优势,红外高光谱数据是全球气象观测资料的重要组成部分[2]。
美国、欧洲和我国都在发展星载红外高光谱大气探测仪,美国于2002年5月4日发射了载于EOS-Aqua卫星的光栅式红外大气探测仪AIRS(Atomspheric Infrared Sounder),其通道数为2378个。欧洲气象局卫星应用组织MetOp-A/B/C系列星上的IASI(Infrared atmospheric sounder interferometer),其通道数为8461个。美国Suomi NPP和JPSS-1卫星的CrIS(cross-track infrared sounder),通道数为1305。我国自主研制的红外高光谱探测仪也实现了业务在轨运行,装载在风云三号D星的HIRAS,E星的HIRAS-II和装载在风云四号A星的GIIRS,都采用了干涉式分光技术。红外高光谱大气探测仪为高精度仪器,在实验室中经过了严格定标和检测。是在卫星发射后,受到外界因素的干扰,会造成仪器的性能相比于试验结果发生一定偏差,因此需要进行在轨定标。由于仪器自身光学系统和探测器材质、光学系统设计、电子学性能、仪器发射和在轨运行状态的变化和元器件衰减等因素的影响,造成了仪器入瞳辐射值和仪器输出值的非线性响应,如果依然采用线性模型进行近似,势必造成辐射定标误差。
基于目前的研究结果显示,热红外通道普遍存在不同程度的非线性响应,需要引入二次项或更高项定标系数,从而达到高定标精度。因此,实现红外通道定标的非线性处理算法,通过定标数据拟合出热红外通道的辐射定标的非线性项系数是需要解决的技术关键。国内外学者对探测器的非线性响应进行了大量研究,矫正方法可分为硬件校正法和软件校正法,其中硬件方法主要通过补偿电路来校正非线性响应,软件方法主要通过构建干涉数据的非线性响应模型,获取非线性响应系数进而校正探测器产生的光谱畸变。[3]例如美国极轨气象卫星的CrIS[4]以带外(out-of-band)低频光谱范围(50~300 cm-1)的光谱特征为依据,构建卷积方程求解非线性响应系数,校正光谱畸变,本文将此方法简称为带外校正法。欧洲气象局的机载大气测量(FAAM)中的红外干涉仪评价系统(ARIES)[5]选择以波数为50~500 cm-1和 2000~2500 cm-1两个特征区域为判断依据,通过交叉迭代二阶和三阶非线性响应系数对光谱误差进行拟合,校正光谱畸变[3]。此方法研究两个特征区域求解多次项系数,计算复杂度较高。Wu[6]等提出利用光谱响应度作为判断依据计算非线性校正系数,取得了较好辐射校正效果,但此类方法需要依赖大量的变温工况数据,计算量较大。
本文基于上述研究结果,提出了一种基于带内(in-band)光谱的非线性校正方法。以带外校正法求解的非线性校正系数为初值,再通过带内光谱辐射定标模型进行迭代优化,取得了较好的非线性校正效果。
1 高光谱探测仪的非线性
1.1 探测仪非线性响应原理
HIRAS-II采用了干涉式分光技术,测量数据为干涉图。当探测器的输出信号与入射能量呈非线性关系时,需要进行非线性校正。
非线性关系一般可以由如下关系式表示:
其中
其中
将
通过傅里叶变换到频谱域,得到:
其中
1.2 基于带外光谱的非线性校正方法
利用普朗克(Max Planck)定律可以得到理想黑体的光谱辐射能量。再通过仪器的光谱响应函数(Spectral Response Function,SRF)对理想光谱进行采样。光谱响应函数表达了传感器在每个波长处接收的单色辐亮度与入射的单色辐亮度的比值。可以模拟出仪器采集的理想光谱信号值Spc。
其中E为入瞳处的黑体辐射能量。为了得到
将干涉图I求平方然后再进行傅里叶变换,得到光谱自卷积信号Spc2:
以长波为例,设置黑体温度为300K,光谱分辨率为0.625波数。长波光谱有效范围650~1136波数,得到的模拟光谱辐射信号见
将二次项光谱(长波)与原始光谱(长波)合并,得到模拟的非线性特征光谱(长波)和理想光谱(长波)的对比图见
从模拟结果来看,光谱自卷积信号出现在低频(约0~500cm-1)和高频(约1500~2500cm-1)区域,都在长波的有效光谱范围以外,即带外。在带外高频区处会受到与非线性无关的失真信号所影响,例如二次反射谱(double-bounce),所以可以将带外低频区域作为非线性特征的研究区域。CrIS[4]带外低频的研究范围设为:50~300 cm-1,杨[7]等将50~500 cm-1作为研究区域。
根据以上模拟结果可以得出,在带外低频范围内,理想光谱信号值应该为0,所以
进而求解
令
则:
直流电压V需要通过光谱信号进行估算。理论上,直流电压是由仪器的背景发射和探测器的暗电流产生的,可以利用观测冷空的光谱信号进行估计,即
观测黑体和目标的电压增量可以通过观测光谱计算:
得到直流量V后,可对
2 基于带内光谱的非线性校正方法
2.1 辐射定标模型
辐射定标一般采用线性模型:
在轨辐射定标一般采用两点法,利用冷空和变温黑体某一温度点对应的幅度值
其中Re[]代表取实部,
辐射定标模型对采集数据进行线性校正,无法校正非线性误差,因此在进行辐射定标之前要先进行非线性校正。
2.2 带内校正法
利用带外光谱进行非线性校正具有一定的不确定性,会导致校正效果并不理想原因有以下两点:1)在
在求解出非线性校正系数
带入辐射定标方程
利用方程(21)求解非线性系数
基于目标函数进行迭代求解非线性系数。为了提升计算效率,将带外校正法求解的非线性系数作为迭代初值进行运算。
3 非线性校正结果与分析
3.1 在轨定标数据分析
实验数据为FY-3E/HIRAS-II的在轨测量数据。根据设计,HIRAS-II 仪器的常规对地观测模式下扫描镜跨轨扫描工作,每5个扫描行观测32个驻留视场,包括28个连续对地目标、2个冷空和2个星上黑体观测目标(这些测量既包含干涉仪正向摆扫测量,也包含反向摆扫测量)。每个驻留视场称为 FOR(Field of Regard),包含 3×3排列的9个探元,每个探元称为 FOV(Field of View)[2]。直接测量数据为干涉图,经过傅里叶变换得到光谱图像,见
从光谱图像中可以看出,长波、中波1具有明显的非线性特征,中波2没有明显的非线性特征,所以本文进行非线性校正的主要研究对象是长波红外和中波1红外光谱。选取FY-3E/HIRAS-II 2021年11月4日在轨运行的部分数据用于研究,黑体处于降温状态,见
将降温数据分成5个降温曲线,利用带外校正法和带内校正法(本文方法)分别求解非线性系数进行对比,非线性系数结果见
图 6. 非线性校正系数 (a) 长波,带内校正法,(b)长波:带外校正法,(c)中波1:带内校正法,(d) 中波1:带外校正法
Fig. 6. Nonlinear correction coefficient (a) LW: in-band correction method,(b) LW: out-of-band correction method,(c) MW1: in-band correction method,(d) MW1: out-of-band correction method;
从图像中可以看出,在同一个FOV中不同的降温曲线中两种方法求解的非线性系数具有较高的一致性,且两种方法在不同FOV的非线性系数具有相同的趋势。
利用未经过非线性校正光谱数据和经过带内、带外非线性校正方法得到的数据进行对比实验。将三种数据分别进行辐射定标得到的黑体光谱辐射亮温,计算不同温度点黑体的辐射亮温偏差,得到的结果见图
图 7. 不同探测元未校正和带内校正方法的平均亮温偏差 (a)长波,(b)中波1
Fig. 7. Average brightness temperature deviations of mediumwave 1 by uncorrected and in-band correction methods of different FOV (a) LW,(b)MW1
图 8. 不同探测元带内校正法和带外校正法的平均亮温偏差 (a)长波,(b)中波1
Fig. 8. Average brightness temperature deviations of mediumwave 1 by in-band and out-of-band correction methods of different FOV (a) LW,(b)MW1;
从对比结果来看带内校正法校正后黑体不同温度的辐射亮温偏差长波在1K以内,中波1在2K以内,明显小于未校正和带外校正方法的结果。
图 9. 光谱响应度对比 (a)长波校正前,(b)长波校正后,(c)中波1校正前,(d)中波1校正后
Fig. 9. Spectral responsivity comparison (a) longwave before correction,(b) longwave after correction,(c) mediumwave1 before correction,(d) mediumwave1 after correction
对于线性系统,光谱响应度P在不同温度点应该具有一致性。基于此原理可以通过计算光谱响应度来验证非线性校正的效果[8],计算公式如下:
其中
3.2 交叉比对
为了检验HIRAS-II的非线性校正的结果,可以将辐射定标后的数据与国际公认精度较高的Metop-C IASI的观测数据进行交叉比对。基于星下点交叉比对SNO方法,在满足时间、空间和观测几何一致性的条件下,匹配HIRAS-II和IASI的观测数据。从亮温对比图(
图 10. HIRAS-II与IASI亮温对比图像 (a)未非线性校正,(b)带外校正法,(c)带内校正法
Fig. 10. Bright-temperature comparation of HIRAS-II and IASI (a) no nonlinear correction,(b) out-of-band correction,(c) in-band correction
选取HIRAS-II与IASI在2021-11-02的交叉匹配数据(见
图 12. 与IASI对地目标观测的平均辐射亮温偏差
Fig. 12. The average radiant brightness temperature bias of ground targets compared with IASI
表 1. 平均亮温偏差
Table 1. Average brightness temperature deviation
|
从HIRAS-II和IASI交叉匹配数据的平均亮温偏差可以看出,经过带内校正的辐射亮温偏差在长波和中波1光谱范围在
4 结论
本文针对HIRAS-II在轨定标数据,提出了一种基于带内光谱的非线性校正方法,即带内校正法。通过实验得到以下结论:
(1)FY-3E/HIRAS-II的在轨定标测试中的黑体的长波、中波1光谱数据具有较明显的非线性特征;
(2)利用本文方法计算的非线性系数在不同温况下具有较高的一致性。利用辐射定标后的光谱计算光谱响应度得出,经过本文方法非线性校正后的不同温度黑体的观测数据的光谱响应度具有较好的一致性;
(3)将本文方法校正后的数据进行辐射定标,在不同温度下,得到的黑体光谱辐射亮温偏差长波在1K以内,中波1在2K以内,明显小于未校正和带外校正法的结果;
(4)将HIRAS-II辐射定标后的在轨数据与IASI进行交叉比对。在满足时间、空间和观测几何一致性的条件下,经过本文方法非线性校正的目标观测亮温与IASI的匹配度明显优于带外校正法和没有经过非线性校正的光谱,提升了高光谱数据的对地观测精度。
[2] 左丰华, 胡秀清, 王霞, 等. FY-3E同平台成像仪对HIRAS-II定位与定标精度评估[J]. 光学学报, 2022, 42(24): 240-248.
Hu Xiuqing, et al. Positioning and Calibration Accuracy Evaluation of HIRAS-II by FY-3E Imager on Same Platform[J]. Acta Optica Sinica, 2022, 42(24): 240-248.
[3] 孙永丰, 徐亮, 沈先春, 等. 红外光谱辐射计探测器高阶非线性响应校正方法[J]. 物理学报, 2021, 70(06): 159-166.
Sun Yongfeng, Xu Liang, Shen xianchun, et al. High-order nonlinear response correction method for infrared radiation detector[J]. Acta Physica Sinica, 2021, 70(06): 159-166.
[5] FiedlerL,NewmanS, Bakan S 2005 Appl. Opt. 445332
[6] Wu Chunqiang, Qi Chengli, Hu Xiuqing, et al. FY-3D HIRAS radiometric calibration and accuracy assessment[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 58(6): 3965-3976.
[7] 杨敏珠, 邹曜璞, 张磊, 等. 干涉数据中非线性的校正及其对辐射定标的影响[J]. 中国激光, 2017, 44(01): 272-278.
Yang Minzhu, Zou Yaopu, Zhang Lei, et al. Correction to Nonlinearity in Interferometric Data and Its Effect on Radiometric Calibration[J]. Chinese Journal of Lasers, 2017, 44(01): 272-278.
[8] 杨天杭, 顾明剑, 邵春沅, 等. FY-3E/HIRAS-Ⅱ发射前热真空定标试验非线性校正[J]. 红外与毫米波学报, 2022, 41(03): 597-607.
Yang Tianhang, Gu Mingjian, Shao Chunyuanet al. Nonlinearity correction of FY-3E HIRAS-II in pre-launch thermal vacuum calibration tests[J]. Journal of Infrared and Millimeter Waves, 2022, 41(03): 597-607.
Article Outline
黄硕, 顾明剑, 胡勇, 杨天杭, 邵春沅, 张春明. 红外高光谱大气探测仪(FY-3E/HIRAS-II)在轨数据非线性校正方法[J]. 红外与毫米波学报, 2024, 43(1): 98. Shuo HUANG, Ming-Jian GU, Yong HU, Tian-Hang YANG, Chun-Yuan SHAO, Chun-Ming ZHANG. Nonlinear response correction method for on-orbit data of FY-3E hyperspectral infrared atmospheric sounder II[J]. Journal of Infrared and Millimeter Waves, 2024, 43(1): 98.