FY-3E微波湿度计发射前定标的非线性系数在星上定标中的适用性分析

王振占; 肖雨伟; 张升伟; 何杰颖; 谷松岩

doi:10.11834/jrs.20221447

国产遥感卫星历史数据再定标共性技术 | 浏览量 : 756 下载量: 740 CSCD: 2 更多指标

R-PDF
PDF
导出
分享给微信好友或者朋友圈

分享
收藏
专辑

FY-3E微波湿度计发射前定标的非线性系数在星上定标中的适用性分析
Analysis on applicability of nonlinearity coefficients derived from prelaunch calibration tests to onboard calibration of Microwave Humidity Sounder （MWHS-Ⅲ） on FY-3E satellite
- 角色： First author 第一作者
- 机构：
  中国科学院国家空间科学中心中国科学院微波遥感技术重点实验室, 北京 100190
- 邮箱： wangzhenzhan@mirslab.cn
- 简介：王振占，研究方向为微波遥感定标、定量反演与应用技术。E-mail： wangzhenzhan@mirslab.cn
王振占
1 ，
- 角色： Corresponding author 通讯作者
- 机构：
  中国科学院国家空间科学中心中国科学院微波遥感技术重点实验室, 北京 100190
  中国科学院大学, 北京 100049
- 邮箱： 15611650126@163.com
- 简介：肖雨伟，研究方向为微波遥感定标及应用技术。E-mail： 15611650126@163.com
肖雨伟
1 2 ，
- 机构：
  中国科学院国家空间科学中心中国科学院微波遥感技术重点实验室, 北京 100190
张升伟
1 ，
- 机构：
  中国科学院国家空间科学中心中国科学院微波遥感技术重点实验室, 北京 100190
何杰颖
1 ，
- 机构：
  国家卫星气象中心(国家空间天气监测预警中心) 中国气象局, 北京 100081
  中国遥感卫星辐射测量和定标重点开放实验室中国气象局, 北京 100081
  许健民气象卫星创新中心, 北京 100081
谷松岩
3 4 5
2023年27卷第10期页码：2327-2336
纸质出版日期： 2023-10-07 ，
DOI： 10.11834/jrs.20221447

扫描看全文

引用

阅读全文PDF
王振占，肖雨伟，张升伟，何杰颖，谷松岩.2023.FY-3E微波湿度计发射前定标的非线性系数在星上定标中的适用性分析.遥感学报，27（10）： 2327-2336

Wang Z Z，Xiao Y W，Zhang S W，He J Y and Gu S Y. 2023. Analysis on applicability of nonlinearity coefficients derived from prelaunch calibration tests to onboard calibration of Microwave Humidity Sounder （MWHS-Ⅲ） on FY-3E satellite. National Remote Sensing Bulletin， 27（10）：2327-2336
王振占，肖雨伟，张升伟，何杰颖，谷松岩.2023.FY-3E微波湿度计发射前定标的非线性系数在星上定标中的适用性分析.遥感学报，27（10）： 2327-2336 DOI： 10.11834/jrs.20221447.

Wang Z Z，Xiao Y W，Zhang S W，He J Y and Gu S Y. 2023. Analysis on applicability of nonlinearity coefficients derived from prelaunch calibration tests to onboard calibration of Microwave Humidity Sounder （MWHS-Ⅲ） on FY-3E satellite. National Remote Sensing Bulletin， 27（10）：2327-2336 DOI： 10.11834/jrs.20221447.

FY-3E微波湿度计发射前定标的非线性系数在星上定标中的适用性分析

王振占 ; 肖雨伟 ; 张升伟 ; 何杰颖 ; 谷松岩 ;

PDF下载

摘要

星载微波湿度计发射前定标和在轨定标的最大不同是冷源的使用：在轨定标的冷源是宇宙冷空背景温度，通常为2.7 K，而发射前定标的冷源通常控制在80—100 K。这样通过发射前定标得到的非线性系数是否可以用于在轨定标一直是一个非常关心和值得研究的问题。本文通过在FY-3E微波湿度计发射前热真空定标中除了传统需要的冷热、变温定标源之外，另外设计了一个独立的位于180—200 K附近的定标源。利用该验证源，本文提出了三点定标公式，基于冷空不变性并通过外推到冷空亮温的模拟观测计算非线性系数。随后通过比较传统和新型非线性系数推导方法的结果，并将两组系数带入同一组数据中分析两组发射前定标系数用于在轨定标的适用性。结论表明，通过对传统方法和新方法的计算结果比较，发现将两组非线性方法应用于同一数据集的结果存在差异，说明传统方法得到的发射前系数存在误差。并且该方法简化了发射前标定过程，提高了非线性系数用于星上标定的精度。

译

Abstract

Nonlinearity is one of the most important errors in calibrating a microwave radiometer. It cannot be corrected in orbit through direct measurements, and its correction usually uses a set of coefficients derived from prelaunch calibration tests, namely, thermal vacuum test. The main difference between prelaunch and on-board calibration of a microwave radiometer is that of the cold calibration target used in two-point calibration. On board calibration usually uses the cold sky, which is also called cosmic background temperature of 2.7 K, as the lower end in the response of a microwave radiometer, whereas in the prelaunch calibration, a blackbody at a temperature between 80 K and 100 K is used in the calibration practice. Therefore, the applicability of the nonlinear coefficients derived from the prelaunch calibration test for on-board calibration is a long-standing concern that has never been validated. In this study, a special program of prelaunch calibration was designed, adding a validation blackbody target, which works at ambient temperature of 180—200 K in the thermal vacuum chamber, in addition to the traditional cold blackbody target, a warm blackbody target, and a varied-temperature blackbody target. This validation target provides an opportunity to compute nonlinear coefficients directly through three-point calibration and to evaluate the errors in applying the nonlinear coefficients by the traditional calibration method. Furthermore, a novel method is invented through extrapolating the measurements to the cold sky temperature by the equation of three-point calibration during a scan. It serves as a virtual on-board calibration target to calibrate the radiometer and to produce nonlinear coefficients combined with three other calibration targets. The nonlinear coefficients of the traditional and novel methods are compared, and the two groups of nonlinearities are applied to the same data sets. Thus, the errors in applying the nonlinear coefficients by the traditional calibration method can be evaluated. The results show a deviation between the two group of nonlinear coefficients; the greater absolute value of the nonlinear coefficients indicates greater error. The nonlinear coefficients of the traditional and novel methods are compared, and discrepancies are found in applying the two groups of nonlinearities to the same data sets. The finding implies that errors exist in the prelaunch coefficients when using the traditional method. The novel method simplifies the procedure of prelaunch calibration and improves the accuracy of nonlinear coefficients using for on-board calibration.

译

关键词

微波湿度计定标; FY-3微波湿度计; 发射前定标; 在轨定标; 非线性系数

译

Keywords

Microwave radiometer calibration; FY-3E Microwave Humidity Sounder; Prelaunch calibration; Nonlinearity; On-board calibration

译

1　引言

星载微波湿度计是用来测量地球表面和大气微波辐射的高灵敏度接收机，在系统设计上通常使其工作在系统响应的线性区域，可以通过两点定标来获取这个线性方程，进而实现对地观测数据的定标。然而在实际的微波湿度计系统中，微波器件的响应特性不可能做到理想的线性，必然会存在实际测量值偏离线性方程对应的理论值，这个偏移量通常称为系统的非线性误差，通常把在观测范围内最大偏离的亮温或者辐亮度称为湿度计的非线性误差。非线性误差是系统响应特性决定的，只有在准确获得系统的响应特性之后，才能确定系统的非线性误差。由于微波湿度计的观测物理量是辐亮度，或者亮温，因此通常非线性系数的单位也是辐亮度或者亮温，用亮温表示的系统非线性误差定义为（

Wang等，2011）：

译

$Δ T_{B} = T_{B 0} - T_{B M}$

(1)

式中，T_B0为目标（热真空实验中通常是指用于地物表面辐射模拟的变温定标源）的真实亮温，T_BM为微波湿度计实际测量的亮温。当不考虑微波湿度计天线方向图影响的时候（实际情况天线方向图的影响需要考虑，这个过程称为天线方向图校正或者定标，这里目的是为了获得非线性系数，而在发射前定标中由于定标源紧邻天线，而且定标源后面足够大，所以天线方向图的影响可以忽略），星载微波湿度计测量的亮温T_BM等于基于两点定标方法得到的天线温度T_A（

Weng等，2013）：

译

$T_{A} = a_{1} V_{A} + b_{1}$

(2)

式中， $a_{1}$ 、 $b_{1}$ 为定标方程的系数， $V_{A}$ 为湿度计观测目标的输出电压。定标系数分别表示为

译

$a_{1} = \frac{T_{H} - T_{C}}{V_{H} - V_{C}}$

(3)

$b_{1} = \frac{V_{H} T_{C} - V_{C} T_{H}}{V_{H} - V_{C}}$

(4)

式中， $V_{H}$ 、 $V_{C}$ 分别为热源电压和冷源电压，单位为伏特或计数值； $T_{H}$ 、 $T_{C}$ 分别是热源和冷源的温度，单位为K。研究表明，两点定标得到的 $T_{B M}$ 并不能表征完整的星载微波湿度计系统，因为该系统还存在非线性成分，而这个非线性成分可以用二次项足够表示出来。微波湿度计的非线性误差与系统的工作状态（如仪器温度、接收机增益）、观测目标的动态范围有关，可以表示为如下形式：

译

$Δ T_{B} = u ρ_{A C} ρ_{A H} (T_{H} - T_{C})^{2}$

(5)

$u$ 为湿度计系统的非线性系数，此外，

译

$ρ_{A C} = \frac{V_{A} - V_{C}}{V_{H} - V_{C}}$

(6)

$ρ_{A H} = \frac{V_{A} - V_{H}}{V_{H} - V_{C}}$

(7)

微波湿度计在轨运行过程中，虽然可以实时获得测量数据，但是 $T_{B 0}$ 无法已知，因此在轨过程中的 $u$ 是未知的，一般基于发射前定标的数据计算 $u$ 。

译

星载微波湿度计发射前定标通常是指热真空定标，即在真空罐中利用液氮致冷的定标源作为冷源，工作在300 K附近的常温源作为定标的热源，在95—335 K变温的定标源，来计算微波湿度计在不同工作温度下的非线性系数（

Saunders等，1995；Mo，1996，2007；Kim等，2014；He等，2015；Wang等，2019），确定星载微波湿度计在不同的工作环境下非线性系数

$u$ 的变化。而在轨运行时星上定标所需非线性系数无法实时计算，是根据湿度计实际工作温度，对发射前热真空定标的非线性系数进行温度插值得到。

译

发射前定标和在轨星上定标的最大不同是冷源的使用：发射前定标的冷源通常控制在80—100 K（微波湿度计通常为95 K，后面以这个为例），而在轨星上定标的冷源是宇宙冷空背景温度2.7 K。从式（5）可知，非线性系数 $u$ 的计算除了与输出电压相关，还与冷热源亮温值有关，所以采用热真空定标中的冷热源和变温源计算的非线性系数不一定适合在轨星上定标非线性修正。所以本文通过分析式（1）并在FY-3E微波湿度计（MWHTS-FY-3E）热真空试验过程中加入独立第3个定标基准，来分析冷空条件下非线性系数的适用性。

译

2　研究方法或原理

虽然微波湿度计系统是一个线性系统，但是对于在较大亮温范围内还是存在较大的非线性，非线性的大小取决于湿度计系统的具体技术指标，由式（5）可知系统非线性误差 $Δ T_{B}$ 是关于 $V_{A}$ 的二次项，根据MWHTS仪器性能和定标特点，因此我们可以假设湿度计系统响应是一个二次曲线，定标方程可以表示为（https：//trs.jpl.nasa.gov/handle/2014/14259［2021-07-08］）：

译

$T_{A^{'}} = a V_{A}^{2} + b V_{A} + c$

(8)

式中， $a$ 、 $b$ 、 $c$ 为定标系数。这种表达方式有助于直接确定输出电压与入射的天线温度之间的关系。从式（8）可以看出，微波湿度计的准确定标需3个有效观测点： $(V_{C}, T_{C})$ 、 $(V_{M}, T_{M})$ 、 $(V_{H}, T_{H})$ 分别代表定标的低温、中间温度和高温点的电压和天线温度。利用这3个点，可以得到二次项系数为

译

$a = \frac{\frac{T_{H} - T_{M}}{V_{H} - V_{M}} - \frac{T_{M} - T_{C}}{V_{M} - V_{C}}}{V_{H} - V_{C}}$

(9)

$b$ 、 $c$ 通过将3个有效观测点代入式（8）解方程即可得到，这里不再列出。另外，结合两点定标公式（式（2））和微波湿度计非线性系数定义式（5）可得：

译

$\begin{array}{l} T_{A^{'}} = T_{A} + △ T_{B} = a_{1} V_{A} + b_{1} + u ρ_{A C} ρ_{A H} (T_{H} - T_{C})^{2} = \\ u {(\frac{T_{H} - T_{C}}{V_{H} - V_{C}})}^{2} {V_{A}}^{2} + [\frac{T_{H} - T_{C}}{V_{H} - V_{C}} - u {(\frac{T_{H} - T_{C}}{V_{H} - V_{C}})}^{2} (V_{H} + V_{C})] V_{A} + \\ u {(\frac{T_{H} - T_{C}}{V_{H} - V_{C}})}^{2} V_{C} V_{H} + \frac{V_{H} T_{C} - V_{C} T_{H}}{V_{H} - V_{C}} \end{array}$

(10)

因此系数 $a$ 还可以表示为

译

$a = u {(\frac{T_{H} - T_{C}}{V_{H} - V_{C}})}^{2}$

(11)

可以证明：对于湿度计定标测量数据的误差精度要求来说，上述两个定义在二次曲线单侧（对称轴的一侧）相同的范围内的表达结果是相同的。结合式（10）和（11），可以得到：

译

$u = \frac{T_{H} - T_{M}}{{(T_{H} - T_{C})}^{2}} \frac{V_{H} - V_{C}}{V_{H} - V_{M}} - \frac{T_{M} - T_{C}}{{(T_{H} - T_{C})}^{2}} \frac{V_{H} - V_{C}}{V_{M} - V_{C}}$

(12)

从式（12）可知，如果存在一个目标亮温 $T_{M}^{'}$ ，其输出电压位于 $V_{C}$ 、 $V_{H}$ 中间，则有：

译

$\frac{V_{H} - V_{C}}{V_{H} - V_{M}} = \frac{V_{H} - V_{C}}{V_{M} - V_{C}} = 2$

(13)

因此：

译

$u = \frac{2 (T_{H} + T_{C} - 2 T_{M}^{'})}{{(T_{H} - T_{C})}^{2}}$

(14)

式（14）表示当寻找到输出的电压位于冷源和热源电压的中间值时，非线性系数与冷热源输出电压无关，只与冷热源温度和中间电压对应的定标源的温度有关。式（14）的另外含义表明系统的非线性系数与冷热源的温度有关，发射前定标和发射后在轨定标的非线性系数是不同的。在冷空背景下， $T_{C}$ 的值很小，对 $u$ 的结果值作用几乎为零，当忽略 $T_{C}$ 的贡献，式（14）可变为

译

$u \approx \frac{2}{T_{H}} - \frac{4 T_{M}^{'}}{T_{H}^{2}} = \frac{2}{T_{H}} (1 - \frac{2 T_{M}^{'}}{T_{H}})$

(15)

这个公式是非线性系数的一个近似表达式，可以用来估算在轨非线性系数的大小量级。非线性的大小取决于 $2 T_{M}$ 和 $T_{H}$ 的相对大小，进而决定 $u$ 的符号。当 $2 T_{M} = T_{H}$ 的时候，系统的非线性系数 $u$ 近似为0。

译

另外式（12）直接提供了一种非线性系数的计算方法，通过3个定标点就可以计算对应范围的非线性系数。由于发射前热真空环境和在轨条件的差异，通常的发射前真空定标无法直接获取冷空定标的非线性系数。但是当增加一个中间亮温基准的时，就可以通过式（8）实时计算冷空背景 $T_{B C}$ 所对应的湿度计输出电压 $V_{B C}$ ，进而通过式（12）计算出在轨利用冷空作为定标低端的非线性系数。这种方法有望减小发射前热真空定标获得非线性的复杂度和成本。

译

3　发射前确定系统非线性系数的实验方案及其数据分析方法

3.1　实验方案

试验设计的方案如图1所示，真空罐内的热沉温度保持在100 K左右（

王振占等，2013）。除了微波湿度计系统自带的热源（位于微波湿度计的环境温度）以外，试验还包括3个定标源：变温源、冷源和验证源。变温源温度在95—335 K范围变化，在试验中步进温度为15 K，用来模拟地球表面亮度温度变化；冷源的温度稳定在95 K附近，使用液氮系统冷却；验证源不进行温控，只是采取被动包扎的方式，利用铂电阻实测其辐射体温度，最终的温度恒定在某个平衡温度。图1给出微波湿度计及定标源装配示意图，以天底位置（49，50像元之间）为0度参考，各定标源所在角度和观测位置如表1所示。

译

图1 热真空定标微波湿度计及定标源装配示意图

Fig. 1 Installing structure of MWHTS and three targets onboard the fixed test platform during T/V test

下载: 原图 | 高精图 | 低精图

表1 各定标源所在角度和观测位置

Table 1 The angle and observation position of each calibration source

定标源	所在角度	所在位置	观测像元
验证源	37.4°	第83—84像元	第68—98个
变温源	17.6°	第33—34像元	第18—48个
冷源	72.9°	冷空位置	冷空视场

下载: 导出CSV

下载: 下载表格图片

试验过程中，仪器温度的变化通过温控系统调整微波湿度计的底板温度在0 ℃—30 ℃变化范围来实现，工作温度设定为5 ℃、12 ℃、18 ℃、25 ℃。通过分析计算可以获取湿度计非线性系数及其误差、辐射灵敏度、线性度、辐射测量精度等主要技术指标，同时验证源可以独立验证非线性系数的定标结果。

译

3.2　冷空不变性及虚拟冷空观测的概念

MWHTS-FY-3E采用宇宙背景温度作为低温参考点，结合卫星上微波湿度计自带的定标黑体，在每个扫描周期2.667 s内产生实时两点定标方程，对这个扫描周期内98个对地观测数据进行定标，获得其对应的天线温度，进而经过天线方向图校正，获得场景亮温（

Wang等，2019）。对于MWHTS-FY-3E的15个通道，冷空的理论亮温见表2。这个亮温只与频率相关，与观测时间是无关的。因此对于一个系统指标确定的微波湿度计而言，其冷空亮温是不变的，我们称为冷空不变性。冷空不变性定标原理就是利用冷空的理论亮温是不随时间的变化而变化，在发射前利用这个亮温来模拟系统在轨输出电压，进而确定在轨定标所需要的非线性系数。冷空背景亮温可以表示为（Kramm和Mölders，2009）：

译

$T_{B C} = 2 h c^{2} (\frac{υ^{3}}{e^{\frac{h c υ}{k T_{c o s}}} - 1})$

(16)

式中，Planck常数 $h = 6.626 \times 10^{- 34} J \cdot S$ ；玻尔兹曼常数k=1.381×10^-23 J/K；光速c=2.9979×10¹⁰ cm/s；波数 $υ$ 为 $\frac{c}{f}$ ，单位为cm^-1；频率 $f$ 单位Hz；T_cos=2.7 K。FY-3E的MWHS的4个中心频率的冷空亮温如表1所示。

译

表2 MWHTS冷空理论亮温值

Table 2 Cold space Theoretical brightness temperature of MWHTS

通道中心频率/GHz	等效冷空亮温/K
89.0	3.24
118.75	3.63
166.0	4.42
183.31	4.75

下载: 导出CSV

下载: 下载表格图片

根据第2节的原理，在一个扫描周期内，利用3个不同亮温点的观测数据就可以得到一条二次定标曲线，这个曲线就是当前周期内系统的真实响应曲线。利用这个曲线，可以外推得到湿度计在当前观测周期（由于辐射计系统增益和接收机噪声温度的变化，只能在一个定标周期内的外推才能有效）和动态范围内的任意一个亮温对应的观测电压，这个过程称为虚拟观测。基于这个虚拟观测概念，可以很容易产生当前扫描周期内冷空的虚拟观测电压。在虚拟过程中，冷空是实际在轨的冷空，是不变的，而其电压是虚拟的，是通过实测三点定标的方程实时外推的结果。实时三点定标曲线代表了当前扫描周期内稳定的响应特性，在其动态范围内都是成立的，因此利用这个曲线获得的当前扫描周期的冷空亮温和电压（ $T_{B C}$ 、 $V_{B C}$ ）也是成立的。用这对冷空亮温和电压，结合其他的定标源亮温及电压就可以计算在轨时的非线性系数以及其他性能参数。作为例子，图2给出1和10通道（89 GHz）在5 ℃时变温源从95—330 K变化（步长15 K）时观测4个定标源并外推冷空的电压比较。

译

下载: | 高精图 | 低精图

图2 通道1和10在5 ℃时观测4个定标源和外推冷空电压及其亮温的分布比较

Fig. 2 The voltage and brightness temperature distributions of four calibration sources observed and cold space extrapolated of ch1， 10 at 5 ℃

在图2中，每一个点都是200包的平均值，代表了一个变温源温度点全部的有效观测。可以看出对于一个定标曲线，可以实时产生冷空的观测值。通道1输出的冷空电压值较为一致，而通道10的输出冷空电压在变温源95 K时与其他温度点不一致，这个点是整个定标试验的第一个点，系统还没有完全处于稳定状态。图2中外推得到的冷空虚拟观测电压也是200包数据的平均值，它的准确度只受热源、验证源和冷源亮温准确度的影响。由于虚拟冷空电压是通过外推产生的，因此必须在一个扫描周期内的二次曲线外推得到的冷空电压才是有效的。从这个意义上说，发射前定标源起到基准传递的作用，获取准确的非线性系数的关键是各个定标源的准确度。

译

4　系统非线性系数的确定结果及分析

4.1　系统非线性系数的确定方法及其适用性分析

在获得冷空观测电压 $V_{B C}$ 后，结合宇宙背景亮温 $T_{B C}$ 和其他几个定标源的亮温及其电压，将定标源的偏差去除，就可以利用式（12）计算系统的非线性系数，进而分析不同冷源对于系统非线性系数的影响。为了直观，本文以亮温（单位K）作为辐射量的单位，非线性系数u的单位是“1/K”。

译

根据前面的分析，这里非线性系数的确定方法包括两种：

译

一种是常规定标方法确定的非线性系数，是指通过消除定标源偏差以后，计算变温源亮温和两点定标的亮温差，然后利用式（5）确定的非线性系数（

Wang等，2019），包含了冷源、热源和变温源的亮温及其输出电压；另一种是利用式（15）的特殊点确定非线性系数的方法，包含冷源和热源，以及特殊点。本试验的特殊点可以来自验证源，或者变温源。详见下述分析。由于验证源的温度不进行控制，其输出电压刚好位于冷热源电压之间，因此特殊点不一定总是存在。参见图3通道1的定标结果，图3（a）是电压，图3（b）是对应的亮度。在上图中，红色点线是冷热源中值的电压，其不但与变温源电压有交点，还与验证源电压有交点；蓝色点线是虚拟冷空和热源中值电压，其只与变温源有交点。因此验证源只能确定冷热源条件下的特殊点方法的非线性系数，而不能产生虚拟冷空的非线性系数。

译

下载: | 高精图 | 低精图

图3 在5 ℃升温循环（5U）通道1的各个定标源实测电压、辐亮度及虚拟冷空电压及其辐亮度在17个温度点（横轴为17个温度点对应的时间序列）的比较

Fig. 3 Comparison of the measured voltage， radiance of each calibration target， virtual cold space voltage and its radiance at 17 temperature points （Horizontal axis is the corresponding time series of 17 temperature points） of ch1 at 5U

为了比较使用冷端定标基准对非线性系数的影响，我们把使用冷源得到的非线性系数称为实测非线性系数，把使用虚拟冷空产生的非线性系数称为虚拟非线性系数。因此上述两系数的差异在于计算过程中使用的冷端定标源的不同。前者使用95 K的定标黑体，后者使用虚拟冷空观测值及其亮温。

译

图4给出通道1在95 K冷源定标和2.7 K冷空定标时，采用常规定标和变温源特殊点定标4种情况下计算的实测非线性系数（对应图中冷源95 K）和虚拟非线性系数（对应图中虚拟冷空）的比较。对比图中四条折线可以看出：在实测的非线性系数和虚拟的非线性系数存在一定的差异，也就是说发射前采用95 K冷源定标获得的非线性系数与在轨采用冷空的非线性系数是不同的。由于验证源在虚拟冷空条件下没有中间温度点存在，因此图中省略了特殊点方法获得的非线性系数的比较。其它各个通道的结果总体上类似，但差异的大小略有不同，由于篇幅关系，这里不再赘述。

译

图4 通道1在95 K冷源定标和2.7 K冷空定标时，采用两种定标方法的非线性系数的比较（常规定标和变温源特殊点定标得到的实测非线性系数对应图中冷源95 K，虚拟非线性系数对应图中虚拟冷空）

Fig. 4 Comparison of nonlinear coefficients of the two calibration methods in 95 K cold target calibration and 2.7 K cold space calibration of ch1 （The measured nonlinear coefficient obtained from the conventional calibration and the special point calibration of the earth target corresponds to the cold target 95 K and the virtual nonlinear coefficient corresponds to the virtual cold space in the figure）

下载: 原图 | 高精图 | 低精图

下面以通道1为例，分析在不同非线性系数下不同方法定标误差之间的差异。图5（a）给出两种方法的实测非线性系数误差；图5（b）给出虚拟非线性系数的误差；图5（c）给出将常规定标方法得到的非线性系数用于在轨条件下的定标误差。从图5可以明显看出，传统的实测非线性定标对于发射前测量数据来讲误差很小，基本可以忽略（<0.1 K）。如果冷源采用虚拟冷空，其定标误差就会略微增大，总体位于0.1—0.2 K。但是把传统的实测非线性系数代入虚拟冷空情况下的观测数据，其误差明显增大，整体达到±1 K，同时结果中仍然存在非线性成分，而且出现单调增加的趋势，这说明这个非线性系数不适用于在轨情况下非线性误差的校正。这是试验使用验证源验证传统方法的非线性系数适用性的目的。

译

下载: | 高精图 | 低精图

图5 不同非线性系数下不同方法定标误差之间差异的比较

Fig. 5 Comparison of calibration errors of different methods under different nonlinear coefficients

在图4中，每个通道的非线性系数是变温源在110—330 K范围内16个点非线性系数的平均值（去掉了第一个观测点，此时定标环境温度尚未稳定）。为了分析非线性系数变化产生的定标误差，图6给出两种情况下利用变温源的实测值和虚拟观测值获得的非线性系数均值及其标准差引起的定标误差随着变温源温度的变化。从图中可以看出，非线性系数波动导致的定标误差基本上在0.2 K以内，发射前实测数据的定标误差略小于在轨的误差，除了330 K情况下的误差较大，这与定标过程中的温度稳定性相关，这里不再深入分析。

译

下载: | 高精图 | 低精图

图6 在5 ℃升温循环（5U）通道1的非线性系数波动引起的定标误差随着变温源温度的变化

Fig. 6 The calibration error caused by nonlinearity coefficient fluctuation varies with the temperature of earth targets of ch1 at 5U

4.2　系统非线性系数的确定方法及其适用性分析

无论是以前利用冷热源和变温源获得的非线性系数，还是现在直接三点定标得到的非线性系数，定标源的亮温的准确性是这种方法使用的关键和难点，也是获得的 $u$ 是否有效的决定性因素。我们知道，定标源亮温的准确度无法直接测量，只能通过部分测试结合计算得到。而定标源的发射率、反射率、温度准确度、温度梯度等因素是制约准确度的关键参数。

译

验证源由于处于温度平衡状态没有温度梯度，因此其偏差假设为零。冷源和热源的温度梯度是存在的。热真空罐的温度约100 K，接近冷源温度95 K，二者温差较小，但是真空罐中存在湿度计和定标热源等高温热辐射体，所以可能会导致冷源表面的温度梯度较大。热源处于微波湿度计内部，直接受到外界温度辐射的影响和间接受到由变温源升降温带来的环境温度变化的影响不大，因此温度梯度也较小。变温源是处于不断加热和制冷控制的，其温度梯度受到自身温度和外界环境温度之差的影响而不断变化，所以温度梯度最大的是变温源。

译

传统的方法不修正变温源的亮温，而是修正冷热源相对于变温源的冷热偏差。如果变温源自身亮温不存在非线性，那么两种方法是等价的。如果变温源自身亮温存在非线性，使用传统的非线性计算方法就会把这个非线性传递给湿度计的非线性。为了避免引入变温源的非线性问题，本文首先基于三点绝对定标，对变温源的亮温进行修正。图7和图8给出利用三点定标方法湿度计实测的变温源亮温相对于原始亮温（利用物理温度和发射率计算的亮温。本文使用的4个定标源的发射率>0.9992，因此发射率导致的误差可以忽略，主要来自温度梯度的误差）的差（后者减去前者）。实测数据包括4条测量曲线3组AGC下全部17个点的温度采集数据，图中每个温度点的数据有180组，一组数据包括200包有效观测数据，包含了每个点定标过程中低频两个通道和高频两个通道单独定标的温度数据。图中结果也证明了变温源温度存在误差，而且这些误差是稳定的，与频率的相关性较小。同时从图中可以看出变温源亮温确实存在非线性，因此必须事先去除，以得到准确的接收机非线性。

译

下载: | 高精图 | 低精图

图7 原始变温源亮温减去真实变温源亮温差的实测结果

Fig. 7 The actual result of the difference between the original earth target and the real earth target

图8 原始变温源亮温减去真实变温源亮温差的统计结果

Fig. 8 The statistical result of the difference between the original earth target and the real earth target

下载: 原图 | 高精图 | 低精图

下面简单分析非线性系数对于定标源的敏感性。根据式（14）可以得到：

译

$\frac{\partial u}{\partial T_{H}} = \frac{8 T_{M} - 2 T_{H} - 6 T_{C}}{{(T_{H} - T_{C})}^{3}}$

(17)

$\frac{\partial u}{\partial T_{C}} = \frac{- 8 T_{M} + 2 T_{C} + 6 T_{H}}{{(T_{H} - T_{C})}^{3}}$

(18)

$\frac{\partial u}{\partial T_{M}} = \frac{- 4}{{(T_{H} - T_{C})}^{2}}$

(19)

从式（19）可以看出，非线性系数受到 $T_{M}$ 的影响只与 $T_{H}$ 和 $T_{C}$ 的范围有关，与观测的电压无关。表2给出基于前文的观测数据得到的通道1发射前和发射后非线性系数对于定标源亮温的敏感性。可见虚拟冷空情况下，对热源的敏感性增强，对冷空的敏感性减弱，对特殊点的敏感性变化较小。但是从数据的量级来看，特殊点的影响要远大于冷源和热源，因此在轨定标中，如果根据中间点的亮温来确定在轨非线性系数，则需要较高的亮温精度，否则可能对于非线性系数带来较大的误差。

译

表3 通道1发射前和发射后非线性系数对于定标源亮温的敏感性比较

Table 3 The comparison of sensitivity of nonlinear coefficients to the brightness temperature of the calibration source before and after launch in ch1

状态	热源	冷源	特殊点
发射前	$8.19 \times 10^{- 5}$	$4.05 \times 10^{- 5}$	$- 1.22 \times 10^{- 4}$
在轨	$9.08 \times 10^{- 5}$	$3.13 \times 10^{- 5}$	$- 1.22 \times 10^{- 4}$

下载: 导出CSV

下载: 下载表格图片

5　结论

本文围绕微波湿度计系统非线性系数的确定方法及其适用性开展研究，通过首次在热真空定标中引入独立的验证源，基于微波湿度计在一个或多个定标周期内系统处于稳态响应的特性，提出了基于三点定标的方法来确定在轨非线性系数确定方法、虚拟冷空定标和虚拟观测的概念，同时提出了基于特殊点的非线性系数计算方法。在试验中，这个特殊点可以来自变温源，也可以来自验证源。文中的特殊点定标是三点定标的一个特例，其提供了一种计算非线性系数的特殊方法，其对于分析非线性系数的敏感性、简化定标过程，以及在轨定标中如何确定系统的非线性系数具有一定的理论和指导意义。

译

随后以FY-3E的通道1和10为例，比较了各种方法的非线性系数结果，以及定标误差，并对传统发射前非线性系数的在轨应用误差进行了计算。

译

文章最后对于非线性系数对于定标源亮温的敏感性进行了计算，结果表明：在轨冷空的非线性系数敏感性相对发射前使用95 K附近冷源的敏感性减弱，而在轨热源的敏感性相对增强。非线性系数对于特殊点（冷热源电压中值）亮温的敏感性要远大于冷源和热源，因此在轨定标中，如果根据中间点的亮温来确定在轨非线性系数，则这个中间点需要较高的亮温精度，否则可能对于非线性系数带来较大的误差。

译

致谢

致谢:感谢FY-3E真空定标团队的辛苦工作并提供基础数据。

参考文献（References）

He J Y, Zhang S W and Wang Z Z. 2015. Advanced microwave atmospheric sounder (AMAS) channel specifications and T/V calibration results on FY-3C satellite. IEEE Transactions on Geoscience and Remote Sensing, 53(1): 481-493 [DOI: 10.1109/TGRS.2014.2324173] [百度学术]

Kim E, Lyu C H J, Anderson K, Leslie R V and Blackwell W J. 2014. S-NPP ATMS instrument prelaunch and on-orbit performance evaluation. Journal of Geophysical Research: Atmospheres, 119(9): 5653-5670 [DOI: 10.1002/2013JD020483] [百度学术]

Kramm G and Mölders N. 2009. Planck’s blackbody radiation law: presentation in different domains and determination of the related dimensional. arXiv: 0901.1863 [DOI: 10.48550/arXiv.0901.1863] [百度学术]

Mo T. 1996. Prelaunch calibration of the advanced microwave sounding unit-A for NOAA-K. IEEE Transactions on Microwave Theory and Techniques, 44(8): 1460-1469 [DOI: 10.1109/22.536029] [百度学术]

Mo T. 2007. Postlaunch calibration of the NOAA-18 advanced microwave sounding unit-A. IEEE Transactions on Geoscience and Remote Sensing, 45(7): 1928-1937 [DOI: 10.1109/TGRS.2007.897451] [百度学术]

Saunders R W, Hewison T J, Stringer S J and Atkinson N C. 1995. The radiometric characterization of AMSU-B. IEEE Transactions on Microwave Theory and Techniques, 43(4): 760-771 [DOI: 10.1109/22.375222.] [百度学术]

Wang Z Z, Li J, Zhang S W and Li Y. 2011. Prelaunch calibration of microwave humidity sounder on China’s FY-3A meteorological satellite. IEEE Geoscience and Remote Sensing Letters, 8(1): 29-33 [DOI: 10.1109/LGRS.2010.2050676] [百度学术]

Wang Z Z, Li J Y, He J Y, Zhang S W, Gu S Y, Li Y, Guo Y and He B Y. 2019. Performance analysis of microwave humidity and temperature sounder onboard the FY-3D satellite from prelaunch multiangle calibration data in thermal/vacuum test. IEEE Transactions on Geoscience and Remote Sensing, 57(3): 1664-1683 [DOI: 10.1109/TGRS.2018.2868324] [百度学术]

Wang Z Z, Zhang S W, Li J, Li Y and Wu Q. 2013. Thermal/vacuum calibration of microwave humidity sounder on FY-3B satellite. Strategic Study of CAE, 15(10): 33-46, 53 [百度学术]

王振占, 张升伟, 李靖, 李芸, 吴琼. 2013. FY-3B卫星微波湿度计热真空定标方法和结果分析. 中国工程科学, 15(10): 33-46, 53 [DOI: 10.3969/j.issn.1009-1742.2013.10.005] [百度学术]

Weng F Z, Zou X L, Sun N H, Yang H, Tian M, Blackwell W J, Wang X, Lin L and Anderson K. 2013. Calibration of Suomi national polar-orbiting partnership advanced technology microwave sounder. Journal of Geophysical Research: Atmospheres, 118(19): 11187-11200 [DOI: 10.1002/jgrd.50840] [百度学术]

摘要

Abstract

Nonlinearity is one of the most important errors in calibrating a microwave radiometer. It cannot be corrected in orbit through direct measurements

and its correction usually uses a set of coefficients derived from prelaunch calibration tests

namely

thermal vacuum test. The main difference between prelaunch and on-board calibration of a microwave radiometer is that of the cold calibration target used in two-point calibration. On board calibration usually uses the cold sky

which is also called cosmic background temperature of 2.7 K

as the lower end in the response of a microwave radiometer

whereas in the prelaunch calibration

a blackbody at a temperature between 80 K and 100 K is used in the calibration practice. Therefore

the applicability of the nonlinear coefficients derived from the prelaunch calibration test for on-board calibration is a long-standing concern that has never been validated. In this study

a special program of prelaunch calibration was designed

adding a validation blackbody target

which works at ambient temperature of 180—200 K in the thermal vacuum chamber

in addition to the traditional cold blackbody target

a warm blackbody target

and a varied-temperature blackbody target. This validation target provides an opportunity to compute nonlinear coefficients directly through three-point calibration and to evaluate the errors in applying the nonlinear coefficients by the traditional calibration method. Furthermore

a novel method is invented through extrapolating the measurements to the cold sky temperature by the equation of three-point calibration during a scan. It serves as a virtual on-board calibration target to calibrate the radiometer and to produce nonlinear coefficients combined with three other calibration targets. The nonlinear coefficients of the traditional and novel methods are compared

and the two groups of nonlinearities are applied to the same data sets. Thus

the errors in applying the nonlinear coefficients by the traditional calibration method can be evaluated. The results show a deviation between the two group of nonlinear coefficients; the greater absolute value of the nonlinear coefficients indicates greater error. The nonlinear coefficients of the traditional and novel methods are compared

and discrepancies are found in applying the two groups of nonlinearities to the same data sets. The finding implies that errors exist in the prelaunch coefficients when using the traditional method. The novel method simplifies the procedure of prelaunch calibration and improves the accuracy of nonlinear coefficients using for on-board calibration.

关键词

微波湿度计定标FY-3微波湿度计发射前定标在轨定标非线性系数

Keywords

Microwave radiometer calibrationFY-3E Microwave Humidity SounderPrelaunch calibrationNonlinearityOn-board calibration

references

Mo T. 1996. Prelaunch calibration of the advanced microwave sounding unit-A for NOAA-K. IEEE Transactions on Microwave Theory and Techniques, 44(8): 1460-1469 [DOI: 10.1109/22.536029http://dx.doi.org/10.1109/22.536029]

Wang Z Z, Zhang S W, Li J, Li Y and Wu Q. 2013. Thermal/vacuum calibration of microwave humidity sounder on FY-3B satellite. Strategic Study of CAE, 15(10): 33-46, 53

王振占, 张升伟, 李靖, 李芸, 吴琼. 2013. FY-3B卫星微波湿度计热真空定标方法和结果分析. 中国工程科学, 15(10): 33-46, 53 [DOI: 10.3969/j.issn.1009-1742.2013.10.005http://dx.doi.org/10.3969/j.issn.1009-1742.2013.10.005]

文章被引用时，请邮件提醒。

提交

暂无数据