摘要
风云三号E星(FY-3E)搭载的高光谱大气探测仪(HIRAS-II)能够实现大气的垂直探测,具有高光谱、高灵敏度、高精度的特点。仪器在轨之后由于仪器衰减和环境变化的原因产生非线性响应,影响在轨定标精度。针对非线性响应的问题,提出了一种基于带内光谱的非线性校正方法。首先基于带外低频光谱的非线性特征求解非线性校正系数,将此系数作为初值输入到辐射定标模型中,以星上测量的黑体带内光谱与理想光谱的偏差为目标函数,通过迭代优化非线性校正系数。通过辐射定标实验得出,校正后的黑体亮温偏差明显低于未校正和基于带外光谱的校正方法。将HIRAS-II的观测数据与IASI进行交叉比对并计算平均亮温偏差和偏差绝对值,经过带内校正法非线性校正后的亮温平均偏差为-0.13 K,优于带外校正方法。
风云三号E星(FY-3E)是全球首颗民用晨昏极轨气象卫星,于2021年7月5日在酒泉卫星发射中心成功发射。风云三号E星经过在轨测试和校准后投入使用,与风云三号上午星C星和风云三号下午星D星组网运行,可每6小时为数值天气预报提供一次全球覆盖的卫星资料,有效补充了6小时同化时间窗内卫星观测资料的空
美国、欧洲和我国都在发展星载红外高光谱大气探测仪,美国于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,都采用了干涉式分光技术。红外高光谱大气探测仪为高精度仪器,在实验室中经过了严格定标和检测。是在卫星发射后,受到外界因素的干扰,会造成仪器的性能相比于试验结果发生一定偏差,因此需要进行在轨定标。仪器自身光学系统和探测器材质、光学系统设计、电子学性能、仪器发射和在轨运行状态的变化和元器件衰减等因素的影响,造成了仪器入瞳辐射值和仪器输出值的非线性响应,如果依然采用线性模型进行近似,势必造成辐射定标误差。
基于目前的研究结果显示,热红外通道普遍存在不同程度的非线性响应,需要引入二次项或更高项定标系数,从而达到高定标精度。因此,实现红外通道定标的非线性处理算法,通过定标数据拟合出热红外通道的辐射定标的非线性项系数是需要解决的技术关键。国内外学者对探测器的非线性响应进行了大量研究,矫正方法可分为硬件校正法和软件校正法,其中硬件方法主要通过补偿电路来校正非线性响应,软件方法主要通过构建干涉数据的非线性响应模型,获取非线性响应系数进而校正探测器产生的光谱畸
本文基于上述研究结果,提出了一种基于带内(in-band)光谱的非线性校正方法。以带外校正法求解的非线性校正系数为初值,再通过带内光谱辐射定标模型进行迭代优化,取得了较好的非线性校正效果。
HIRAS-II采用了干涉式分光技术,测量数据为干涉图。当探测器的输出信号与入射能量呈非线性关系时,需要进行非线性校正。

图1 非线性响应示意图
Fig. 1 Nonlinear response diagram
非线性关系一般可以由如下关系式表示:
, | (1) |
其中为入射辐射强度,为系数,为光生导带电子浓度。理想的干涉图光强与入射辐射强度成正比,实测干涉图光强与光生导带电子浓度成正比,所以理想干涉图光强和实测干涉图光强具备如下关系:
, | (2) |
其中表示探测器输出的直流电压,该电压由探测器光电流和暗电流产生,表示理想直流电压。表示理想的干涉图光强,表示实测的干涉图光强,表示二次项系数,表示三次项系数。由于三阶项及以上在探测器中的影响表现不显著,可以省略其非线性影
, | (3) |
将
, | (4) |
通过傅里叶变换到频谱域,得到:
, | (5) |
其中为理想光谱,为实测光谱。非线性引起的光谱畸变由缩放因子和二次项产生,需要根据光谱特征求解二次项系数。
利用普朗克(Max Planck)定律可以得到理想黑体的光谱辐射能量。再通过仪器的光谱响应函数(Spectral Response Function,SRF)对理想光谱进行采样。光谱响应函数表达了传感器在每个波长处接收的单色辐亮度与入射的单色辐亮度的比值。可以模拟出仪器采集的理想光谱信号值Spc。
, | (6) |
其中E为入瞳处的黑体辐射能量。为了得到
. | (7) |
将干涉图I求平方然后再进行傅里叶变换,得到光谱自卷积信号Spc2:
. | (8) |
以长波为例,设置黑体温度为300 K,光谱分辨率为0.625波数。长波光谱有效范围650~1136波数,得到的模拟光谱辐射信号见

图2 模拟光谱信号
Fig. 2 Simulated spectral signal
将二次项光谱(长波)与原始光谱(长波)合并,得到模拟的非线性特征光谱(长波)和理想光谱(长波)的对比图见

图3 模拟非线性数据信号
Fig.3 Simulated nonlinear data signal
从模拟结果来看,光谱自卷积信号出现在低频(0~500 c
根据以上模拟结果可以得出,在带外低频范围内,理想光谱信号值应该为0,所以
, | (9) |
进而求解:
, | (10) |
令
, | (11) |
则
. | (12) |
直流电压V需要通过光谱信号进行估算。理论上,直流电压是由仪器的背景发射和探测器的暗电流产生的,可以利用观测冷空的光谱信号进行估计,即。
, | (13) |
为冷空的观测光谱,n为通道数,G为仪器观测引入的增益。当仪器观测黑体或地面目标时,前置放大器输出的直流电压会产生增量或,则观测黑体或地面目标时的直流量或可以表示为
, | (14) |
. | (15) |
观测黑体和目标的电压增量可以通过观测光谱计算:
, | (16) |
. | (17) |
得到直流量V后,可对
辐射定标一般采用线性模型:
, | (18) |
为探测器第i通道对应的辐亮度值,为探测器数字量化后输出的光谱幅度值,为增益系数,为辐亮度偏差。
在轨辐射定标一般采用两点法,利用冷空和变温黑体某一温度点对应的幅度值、和辐亮度、求解某一观测目标的辐亮度。则在轨辐射定标方程如下:
, | (19) |
其中Re[]代表取实部,代表地面目标的观测数据幅度值,、的理想辐亮度值可以利用普朗克定律求得。
辐射定标模型对采集数据进行线性校正,无法校正非线性误差,因此在进行辐射定标之前要先进行非线性校正。
利用带外光谱进行非线性校正具有一定的不确定性,会导致校正效果并不理想原因有以下两点:1)在
在求解出非线性校正系数后,需要代入
, | (20) |
代入辐射定标方程
. | (21) |
利用方程(21)求解非线性系数,需要已知目标的辐亮度。在轨定标中,将星上黑体变温后的黑体的观测辐度值和理想辐亮度作为目标辐亮度代入方程式中,以辐射偏差最小为目标求解非线性校正系数。目标函数为
. | (22) |
基于目标函数进行迭代求解非线性系数。为了提升计算效率,将带外校正法求解的非线性系数作为迭代初值进行运算。
实验数据为FY-3E/HIRAS-II的在轨测量数据。根据设计,HIRAS-II 仪器的常规对地观测模式下扫描镜跨轨扫描工作,每5个扫描行观测32个驻留视场,包括28个连续对地目标、2个冷空和2个星上黑体观测目标(这些测量既包含干涉仪正向摆扫测量,也包含反向摆扫测量)。每个驻留视场称为 FOR(Field of Regard),包含 3×3排列的9个探元,每个探元称为 FOV(Field of View

图4 实测光谱数:(a) 长波;(b) 中波1;(c) 中波2
Fig.4 Measured spectral data:(a) LW;(b) MW1;(c) MW2
从光谱图像中可以看出,长波、中波1具有明显的非线性特征,中波2没有明显的非线性特征,所以本文进行非线性校正的主要研究对象是长波红外和中波1红外光谱。选取FY-3E/HIRAS-II 2021年11月4日在轨运行的部分数据用于研究,黑体处于降温状态,见

图5 黑体温度变化曲线
Fig.5 Black body temperature curve
将降温数据分成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;
从
利用未经过非线性校正光谱数据和经过带内、带外非线性校正方法得到的数据进行对比实验。将三种数据分别进行辐射定标得到的黑体光谱辐射亮温,计算不同温度点黑体的辐射亮温偏差,得到的结果见

(a)

(b)
图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

(a)

(b)
图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
从对比结果来看带内校正法校正后黑体不同温度的辐射亮温偏差长波在1 K以内,中波1在2 K以内,明显小于未校正和带外校正方法的结果。
对于线性系统,光谱响应度P在不同温度点应该具有一致性。基于此原理可以通过计算光谱响应度来验证非线性校正的效
, | (23) |
其中代表观测目标的光谱,代表了对冷空间观测光谱,代表目标理想光谱,代表冷空间理想光谱。将黑体作为观测目标代入公式,分别计算校正前和校正后在不同温度点的光谱响应度。从

图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
为了检验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的交叉匹配数据(见

图11 交叉匹配点的地理位置
Fig. 11 The geographic location of the cross-match points
非线性校正情况 | 未校正 | 带外校正法 | 带内校正法 |
---|---|---|---|
平均亮温偏差的均值/k | 4.60 | -0.82 | -0.13 |
平均亮温偏差的绝对均值/k | 4.81 | 2.94 | 1.63 |

图12 与IASI对地目标观测的平均辐射亮温偏差
Fig.12 The average radiant brightness temperature bias of ground targets compared with IASI
从HIRAS-II和IASI交叉匹配数据的平均亮温偏差可以看出,经过带内校正的辐射亮温偏差在长波和中波1光谱范围在5 K以内,大部分分布在2 K内,平均亮温偏差的均值为-0.13,绝对均值为1.63,优于带外校正法和没有校正的光谱。说明经过带内校正法非线性校正的观测数据与IASI的匹配精度优于带外校正法和没有经过非线性校正的光谱,此方法能够提升对地面目标的观测精度。
本文针对HIRAS-II在轨定标数据,提出了一种基于带内光谱的非线性校正方法,即带内校正法。通过实验得到以下结论:
(1)FY-3E/HIRAS-II的在轨定标测试中的黑体的长波、中波1光谱数据具有较明显的非线性特征;
(2)利用本文方法计算的非线性系数在不同温况下具有较高的一致性。利用辐射定标后的光谱计算光谱响应度得出,经过本文方法非线性校正后的不同温度黑体的观测数据的光谱响应度具有较好的一致性;
(3)将本文方法校正后的数据进行辐射定标,在不同温度下,得到的黑体光谱辐射亮温偏差长波在1 K以内,中波1在2 K以内,明显小于未校正和带外校正法的结果;
(4)将HIRAS-II辐射定标后的在轨数据与IASI进行交叉比对。在满足时间、空间和观测几何一致性的条件下,经过本文方法非线性校正的目标观测亮温与IASI的匹配度明显优于带外校正法和没有经过非线性校正的光谱,提升了高光谱数据的对地观测精度。
References
Zhang P, Hu X Q, Lu Q F, et al. FY-3E: The first operational meteorological satellite mission in an early morning orbit[J]. Adv. Atmos. Sci., 2022,39(1), 1-8, https://doi.org/10.1007/s00376-021-1304-7 [百度学术]
Zuo Fenghua, Hu Xiuqing, Wang Wia, et al. Positioning and calibration accuracy evaluation of HIRAS-II by FY-3E imager on same platform [J]. Acta Optica Sinica. 左丰华,胡秀清,王霞,等.FY-3E同平台成像仪对HIRAS-II定位与定标精度评估[J].光学学报,2022,42(24):240-248. [百度学术]
Sun Yongfeng, Xu Liang, Shen Xianchun, et al. High-order nonlinear response correction method for infrared radiation detector [J]. Acta Physica Sinica.孙永丰,徐亮,沈先春,等.红外光谱辐射计探测器高阶非线性响应校正方法[J].物理学报,2021,70(06):159-166. [百度学术]
Algorithm theoretical basis document for the cross track infrared sounder (CrIS) sensor data records (SDR). Goddard Space Flight Center. http: // www. star. nesdis. noaa. gov/jpss /documents/ATBD/D0001-M01-S01-002_JPSS_ATBD_CrIS-SDR_C.pdf.2014 [百度学术]
Fiedler L,Newman S, Bakan S2005Appl. Opt. 44 5332 [百度学术]
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 [百度学术]
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.杨敏珠,邹曜璞,张磊,等.干涉数据中非线性的校正及其对辐射定标的影响[J].中国激光,2017,44(01):272-278. [百度学术]
Yang Tianhang, Gu Mingjian, Shao Chunyuan,et al. Nonlinearity correction of FY-3E HIRAS-II in pre-launch thermal vacuum calibration tests [J]. Journal of Infrared and Millimeter Waves.杨天杭,顾明剑,邵春沅,等.FY-3E/HIRAS-Ⅱ发射前热真空定标试验非线性校正[J].红外与毫米波学报,2022,41(03):597-607. [百度学术]