海洋学报

观测量	最小值	最大值
20 Hz测距数目	10	−
20 Hz测距标准差/mm	0	200
轨道高度−测距值/mm	−130 000	100 000
模型干对流层校正/mm	−2 500	−1 900
辐射计湿对流层校正/mm	−500	−1
Ku波段电离层校正/mm	−400	40
Ku波段海况偏差/mm	−500	0
潮汐校正/mm	−5 000	5 000
固体潮校正/mm	−1 000	1 000
极潮校正/mm	−150	150
Ku波段有效波高/mm	0	11 000
Ku波段向后散射系数/dB	7	30
风速/(m·s⁻¹)	−0	30
Ku波段偏指向角平方/（°）²	−0.2	0.64
后向散射系数标准差/dB	−	1
后向散射系数观测数	10	−

观测量

最小值

最大值

20 Hz测距数目

−

20 Hz测距标准差/mm

200

轨道高度−测距值/mm

−130 000

100 000

模型干对流层校正/mm

−2 500

−1 900

辐射计湿对流层校正/mm

−500

−1

Ku波段电离层校正/mm

−400

Ku波段海况偏差/mm

−500

潮汐校正/mm

−5 000

5 000

固体潮校正/mm

−1 000

1 000

极潮校正/mm

−150

150

Ku波段有效波高/mm

11 000

Ku波段向后散射系数/dB

风速/(m·s⁻¹)

−0

Ku波段偏指向角平方/（°）²

−0.2

0.64

后向散射系数标准差/dB

−

后向散射系数观测数

−

观测量	最小值	最大值
20 Hz测距数目	10	−
20 Hz测距标准差/mm	0	200
轨道高度−测距值/mm	−130 000	100 000
模型干对流层校正/mm	−2 500	−1 900
辐射计湿对流层校正/mm	−500	−1
Ku波段电离层校正/mm	−400	40
Ku波段海况偏差/mm	−500	0
潮汐校正/mm	−5 000	5 000
固体潮校正/mm	−1 000	1 000
极潮校正/mm	−150	150
Ku波段有效波高/mm	0	11 000
Ku波段向后散射系数/dB	7	30
风速/(m·s⁻¹)	−0	30
Ku波段偏指向角平方/（°）²	−0.2	0.64
后向散射系数标准差/dB	−	1
后向散射系数观测数	10	−

观测量

最小值

最大值

20 Hz测距数目

−

20 Hz测距标准差/mm

200

轨道高度−测距值/mm

−130 000

100 000

模型干对流层校正/mm

−2 500

−1 900

辐射计湿对流层校正/mm

−500

−1

Ku波段电离层校正/mm

−400

Ku波段海况偏差/mm

−500

潮汐校正/mm

−5 000

5 000

固体潮校正/mm

−1 000

1 000

极潮校正/mm

−150

150

Ku波段有效波高/mm

11 000

Ku波段向后散射系数/dB

风速/(m·s⁻¹)

−0

Ku波段偏指向角平方/（°）²

−0.2

0.64

后向散射系数标准差/dB

−

后向散射系数观测数

−

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	14 746	14 699	14 312
平均值/cm	0.03	0.03	0.03
标准差/cm	6.63	6.89	6.17

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

14 746

14 699

14 312

平均值/cm

0.03

标准差/cm

6.63

6.89

6.17

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	14 746	14 699	14 312
平均值/cm	0.03	0.03	0.03
标准差/cm	6.63	6.89	6.17

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

14 746

14 699

14 312

平均值/cm

0.03

标准差/cm

6.63

6.89

6.17

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	5 005	5 045	4 936
平均值/cm	0.01	0	0.01
标准差/cm	4.60	5.1	4.21

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

5 005

5 045

4 936

平均值/cm

0.01

标准差/cm

4.60

5.1

4.21

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	5 005	5 045	4 936
平均值/cm	0.01	0	0.01
标准差/cm	4.60	5.1	4.21

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

5 005

5 045

4 936

平均值/cm

0.01

标准差/cm

4.60

5.1

4.21

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	22 639	22 610	22 233
平均值/cm	−0.01	−0.01	0
标准差/cm	6.63	5.45	4.70

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

22 639

22 610

22 233

平均值/cm

−0.01

标准差/cm

6.63

5.45

4.70

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	22 639	22 610	22 233
平均值/cm	−0.01	−0.01	0
标准差/cm	6.63	5.45	4.70

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

22 639

22 610

22 233

平均值/cm

−0.01

标准差/cm

6.63

5.45

4.70

	HY-2C和HY-2B	HY-2C和Jason-3
交叉点数（SLA<20 cm）	15 470	15 799
平均值/cm	−0.47	−0.3
标准差/cm	5.32	5.32

HY-2C和HY-2B

HY-2C和Jason-3

交叉点数（SLA<20 cm）

15 470

15 799

平均值/cm

−0.47

−0.3

标准差/cm

5.32

	HY-2C和HY-2B	HY-2C和Jason-3
交叉点数（SLA<20 cm）	15 470	15 799
平均值/cm	−0.47	−0.3
标准差/cm	5.32	5.32

HY-2C和HY-2B

HY-2C和Jason-3

交叉点数（SLA<20 cm）

15 470

15 799

平均值/cm

−0.47

−0.3

标准差/cm

5.32

HY-2C卫星海面高度数据的质量分析

PDF下载

孟菊 ¹ , 于方杰 ¹^,²^,^* , 庄志远 ¹ , 齐娟娟 ¹ , 陈戈 ¹^,²

海洋学报 | 论文 2022,44(10): 173-181

收起

海洋学报 | 论文 2022, 44(10): 173-181

HY-2C卫星海面高度数据的质量分析

全屏

孟菊¹, 于方杰¹^,²^,^*, 庄志远¹, 齐娟娟¹, 陈戈¹^,²

作者信息

1.中国海洋大学信息科学与工程学部海洋技术学院，山东青岛 266100

2.青岛海洋科学与技术试点国家实验室区域海洋学与数值模拟实验室，山东青岛 266237

孟菊（1996-），女，山西省晋城市人，主要从事海洋卫星遥感研究。E-mail：2530511502@qq.com

通讯作者:

于方杰，男，山东省青岛市人，教授，主要从事卫星−无人移动平台协同海洋观测技术研究。E-mail：yufangjie@ouc.edu.cn

Quality analysis of sea level height data of HY-2C satellite

Ju Meng¹, Fangjie Yu¹^,²^,^*, Zhiyuan Zhuang¹, Juanjuan Qi¹, Ge Chen¹^,²

Affiliations

1. Department of Marine Technology, Faculty of Information Science and Engineering, Ocean University of China, Qingdao 266100, China

2. Laboratory for Regional Oceanography and Numerical Modeling, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China

出版时间: 2022-10-01 doi: 10.12284/hyxb2022168

文章导航

摘要

收起

卫星海面高度数据对于监测全球海面高度具有重要的意义，所以卫星高度数据的定标和检验变得至关重要。海洋二号C（HY-2C）卫星是继海洋二号B卫星后的第二颗业务卫星，于2020年成功发射升空。然而，目前对HY-2C卫星高度计的数据质量了解甚少，所以对HY-2C卫星的海面高度数据进行质量分析具有重要的意义。本文以同期观测的HY-2B卫星和Jason-3卫星的地球物理数据（GDR）为参考，对HY-2C卫星遥感地球物理数据（SGDR）中的海面高度数据进行质量分析。结果显示，在星星交叉定标中使用3种常见的交叉定标插值方法对HY-2C卫星的海面高度异常数据进行自交叉点分析时，HY-2C卫星海面高度异常数据质量分析的结果不同。其中使用三次样条插值方法进行质量分析的结果最优，得到海平面高度异常差的平均值为0.03 cm，标准差为6.17 cm。此外，对HY-2C卫星和HY-2B卫星互交叉点海面高度异常差异的平均值为−0.47 cm，标准差为5.32 cm；HY-2C卫星SGDR与Jason-3卫星GDR的海面高度异常数据进行互交叉点分析，得到海平面高度异常差的平均值为−0.3 cm，标准差为5.32 cm，这些数据表明HY-2C卫星的测高精度与HY-2B卫星、Jason-3卫星一致。因此HY-2C高度计产品数据质量稳定，能满足海洋应用和科学研究的需要。

关键词

HY-2C卫星 / 高度计 / 海面高度 / 交叉点分析

Abstract

收起

As an important part of monitoring oceanic phenomenon, it is essential to calibrate altimeter sea level height data. In 2020, the Haiyang-2C (HY-2C) satellite was sent to orbit as a follow-up mission of the HY-2B satellite. However, the quality of HY-2C sea level height is not fully known yet. In this paper, the quality and accuracy of HY-2C sensor geophysical data records (SGDR) of the latest version were assessed through the comparison with the HY-2B geophysical data records (GDR) and Jason-3 GDR measured at the same time. According to the results: self-crossover points analysis is carried out by star cross calibration of HY-2C satellite. The quality analysis results of sea level anomaly data of HY-2C SGDR are different. Through the comparison results, the optimal cubic spline interpolation method has the best quality analysis performance, with mean value is 0.03 cm and the standard deviation is 6.17 cm. The mean and standard deviation of sea level anomaly at the dual-crossing calibration of HY-2C and HY-2B satellites are −0.47 cm and 5.32 cm. The mean and standard deviation of sea level anomaly at the intersection of HY-2C and Jason-3 satellites are −0.3 cm and 5.32 cm. These results show that HY-2C has the same precision as HY-2B satellite and Jason-3 satellite. Data quality of HY-2C satellite is enough to meet the requirements of marine scientific research and application.

Key words

HY-2C satellite / altimeter / sea level height / crossover analysis

引用本文

孟菊, 于方杰, 庄志远, 齐娟娟, 陈戈. HY-2C卫星海面高度数据的质量分析. 海洋学报, 2022 , 44 (10) : 173 -181 . DOI: 10.12284/hyxb2022168

Ju Meng, Fangjie Yu, Zhiyuan Zhuang, Juanjuan Qi, Ge Chen. Quality analysis of sea level height data of HY-2C satellite[J]. Haiyang Xuebao, 2022 , 44 (10) : 173 -181 . DOI: 10.12284/hyxb2022168

正文

收起

1 引言

收起

海洋二号C（HY-2C）卫星是我国第三颗海洋动力环境卫星，该卫星集主、被动微波遥感于一体，属于我国海洋系列遥感卫星，具有高精度测轨、定轨能力与全天候、全天时、全球探测能力。卫星的主要使命是监测和调查海洋环境，获得海面高度等多种海洋动力环境参数，直接为灾害性海况预警预报提供实测数据，为海洋防灾减灾、海洋权益维护、海洋资源开发、海洋环境保护、海洋科学研究以及国防建设等提供支撑服务。

相对定标是国内外学者主要采用的一种数据质量和精度评估方式，其优点是可以实现在全球范围内对数据进行分析，而不存在观测设备高精度要求和地域限制。通过对星间数据的交叉点分析可探测卫星系统异常和地域性误差特征，评估卫星的测量精度，为实现卫星系统改进和有效应用卫星数据提供定量依据^[1-5]。相对定标分为单颗卫星升降轨道的自交叉和不同卫星升降轨道之间的互交叉。卫星自交叉点分析是高度计数据质量评估的重要指示，此方法可用来评价高度计数据的稳定性。互交叉点分析广泛应用于利用校准后的卫星雷达高度计数据对新发射的卫星雷达高度计数据进行评估^[6-7]。已有Jason-2、Jason-3、HY-2A、HY-2B等卫星通过相对交叉定标对海面高度进行了评估^{[1, 3, 6, 8-10]}。用交叉定标获取交叉点的插值方法不尽相同，Jiang等^[6]利用邻近观测点插值的方法获取交叉点，在时间窗口为3 d、纬度小于60°、水深大于200 m、交叉点的高度异常差小于10 cm的条件下对HY-2A卫星进行评估，结果显示HY-2A和Jason-2卫星的交叉点海平面高度异常差值的平均值为4.19 cm，标准差达到4.98 cm，HY-2A卫星成功达到国际同类高度计的测高水平；刘亚龙^[1]利用了三次样条插值的方法获取交叉点，在时间窗口为3 d、纬度小于50°、水深大于1 000 m、海面变化小于20 cm的条件下对HY-2A卫星进行评估，结果显示HY-2A和Jason-2卫星测量的海平面高度差值平均值为0.26 cm，均方根误差为6.61 cm，交叉定标精度达到国际同类高度计精度水平，表明 HY-2A卫星高度计观测海面高度稳定可靠；刘治中等^[9]在Jason-2与Jason-3卫星高度计编队飞行阶段利用最邻近方法获取交叉点，结果显示Jason-3和Jason-2卫星高度计的交叉点海平面高度差值的标准偏差仅在0.06 cm左右。以上为3种比较典型的获取交叉点的方法。而国内对于HY-2C卫星的数据质量分析的文章非常少，目前有Wang等^[11]通过利用浮标、欧洲中期天气预报中心模型和其他散射计风速参数与HY-2C卫星的风速参数进行对比，对HY-2C卫星风速参数的准确性进行了评估。最近Kong等^[12]提出了一种直接利用星载多普勒轨道和卫星无线电定位（DORIS）数据来验证HY-2C卫星轨道精度的新方法。Wang等^[13]提出了一种通过深度学习采集宽幅有效波高的方法。截至目前，还未见到有关HY-2C卫星高度计数据质量的报告。所以实现对HY-2C卫星海面高度精度分析至关重要，有利于实现后续对于HY-2C卫星在海洋观测上的应用。

本文主要对HY-2C卫星高度计遥感地球物理数据（Sensor Geophysical Data Records, SGDR）中的海面高度数据进行全球质量分析。本文以Jason-3卫星和HY-2B卫星同期观测的地球物理数据（Geophysical Data Records, GDR）为参考，对HY-2C卫星的SGDR进行数据完整性和有效性分析，分析HY-2C卫星数据质量；并且使用3种常见的交叉定标插值方法对HY-2B、Jason-3、HY-2C卫星雷达高度计SGDR的海面高度异常数据进行自交叉点分析，探究不同的插值方法对卫星高度计海面高度数据质量分析结果的影响；最后通过分别与HY-2B卫星、Jason-3卫星GDR的海面高度异常数据进行互交叉点计算，完成对HY-2C卫星海面高度数据的质量分析。

2 数据

收起

2.1 数据源

HY-2C卫星高度计SGDR发布于国家卫星海洋应用中心，从2020年9月24日开始发布。用户可以从中国海洋卫星数据服务系统（https://osdds.nsoas.org.cn/）下载，本文使用了HY-2C卫星雷达高度计第13～21周期的SGDR，时间范围是2021年1月21日至4月20日。SGDR是利用MOE（Medium Accuracy Orbit Ephemeris）定轨数据和波形重构等方法得到的未经校正的数据产品，数据中除了主要包括有效波高、海面风速、海面高度及用于计算海面高度所需的相关校正参数，还包含波形数据。因为HY-2B、Jason-3卫星高度计已经在现场经过校准，具有较高的数据质量，所以本文利用HY-2B、Jason-3卫星高度计采集的同期GDR作为对比数据，对HY-2C卫星高度计的海面高度数据进行质量分析。

目前HY-2C卫星只发布了MOE精密定轨的SGDR。从SGDR到GDR的制作流程缺少一个精密定轨校正，因此，利用HY-2B、Jason-3卫星对HY-2C卫星的海面高度数据定标结果可能会因为数据产品的不同存在系统偏差。

2.2 数据处理

HY-2C卫星和HY-2B、Jason-3卫星的海平面高度（Sea Surface Height, SSH）和海平面高度异常（Sea Level Anomaly, SLA）采用下面的公式计算：

(1)

$ \begin{array}{l}{{\rm{SSH}}}={\rm{Altitude}}- ({{{{\rm{Altimeter}}}}}\;{{{{\rm{Range}}}}}+\sum {{{{\rm{Corrections}}}}}) \text{，}\end{array} $

(2)

$ \mathrm{S}\mathrm{L}\mathrm{A}=\mathrm{S}\mathrm{S}\mathrm{H}-\mathrm{M}\mathrm{S}\mathrm{S} \text{，} $

式中，${\rm{Altitude}}$表示轨道的高度；${\rm{Altimeter\;Range}}$表示高度计测距值；$\sum {\rm{Corrections}}$表示干对流层延迟、湿对流层延迟、模型电离层延迟、逆大气压改正、海潮、极潮、固体潮和海况偏差改正项之和；MSS为在一个较长的时期之内剔除了年际、半年际、季节性和其他周期性的海面高度变化信息之后平均的海面高度。本文中HY-2C卫星的SGDR、HY-2B卫星的GDR、Jason-3卫星的GDR中使用的MSS均来自MSS CLS2015模型（https://www.aviso.altimetry.fr/en/data/products/auxiliary-products/mss.html）。该模型是由法国卫星中心利用1993−2012年卫星测高数据中的海平面高度数据解算得到的。SLA可以很好地捕捉海洋变化特征，利用SLA来检验卫星高度计的高度测量能力要比用海平面高度来检验更为准确^{[4-6, 9]}。所以本文用SLA来检验卫星高度计的精度。

2.3 数据完整性与有效性分析

通过数据编辑标准可以剔除掉高度计测量中精度较低的值，得到有效的观测结果^{[7, 10, 14]}。对不同的参数，编辑标准按照最大阈值和最小阈值的方式定义，阈值范围外的参数都需要被剔除。为了提高数据的时空相关性，对HY-2C卫星采用和HY-2B、Jason-3卫星一致的数据编辑策略，本文采用的是Jason-3卫星用户手册中二级产品（O）（I）GDR 的编辑标准。主要编辑标准见表1。

为了准确分析HY-2C卫星高度计在海洋中的测量精度，我们除了通过表1对数据筛选，还进行了以下条件的限制：

（1）为了排除海冰的影响，将纬度限制在南北纬50°之间^[14]。

（2）为了排除近海的影响，将水的深度限制在大于1 000 m的海域^[14]。

（3）为了减少时间变化对海面高度的影响，交叉点处的最大时间间隔为3 d^[14]。

（4）为了剔除一些异常值，卫星交叉点处的SLA之差应该大于20 cm^[6-7,14]。

2.3.1 数据完整性分析

通过分析卫星轨迹数量来检验卫星高度计的完整性。理论上HY-2C卫星1个完整周期有274个轨迹，但实际上由于卫星测量模式的调整或其他原因（地面接收站问题或数据传输异常，也可能是由于地面控制等^{[7, 15]}），导致HY-2C卫星SGDR存在轨迹缺失现象。图1a为HY-2C卫星每周期的轨道数据缺失情况，其中第21周期轨道数据缺失可能是因为数据传输异常，导致出现数据冗余。HY-2B卫星1个完整周期有386条轨迹，图1b是HY-2B卫星每周期轨道数据缺失情况，在和HY-2C卫星相同的时间段内，HY-2B卫星和HY-2C卫星一样存在着数据缺失和冗余现象。图1c是Jason-3卫星每周期的轨道数据缺失情况，Jason-3卫星的GDR 1个完整周期有254个轨迹，在和HY-2C卫星相同的时间段内，Jaosn-3卫星的数据相对稳定，卫星轨道数据未产生缺失和冗余现象。

2.3.2 数据有效性分析

计算SLA时，按照表1的编辑标准对测高数据进行编辑，可以将不符合条件的观测数据滤除。数据的编辑比例指的是不符合约束条件被滤除的数据比例，能够反映出卫星测高数据的质量，编辑比例越低，有效性越好^{[8, 16]}。

根据表1给定的阈值对卫星数据产品进行约束，从而得到HY-2C卫星和Jason-3卫星的数据编辑比例，如图2a所示，在没有经过海深和纬度约束的情况下，HY-2C卫星数据编辑比例的均值为30.7%，标准差为6.62%；HY-2B卫星数据编辑比例的均值为39.9%，标准差为5.29%；Jason-3卫星数据编辑比例的均值为34.17%，标准差为0.60%。在海深1000 m以下和纬度在南北纬50°之间的条件下，得到HY-2C卫星、HY-2B卫星、Jason-3卫星数据编辑比例（图2b），HY-2C卫星数据编辑比例的均值为14.02%，标准差为8.18%；HY-2B卫星数据编辑比例的均值为12.20%，标准差为7.64%；Jason-3卫星数据编辑比例的均值为2.44%，标准差为0.39%。从全球范围来看，HY-2C卫星的有效数据点更多。在海深1 000 m以下和纬度在南北纬50°之间的限制条件下，Jason-3卫星数据有效点更多。并且HY-2C卫星、HY-2B卫星和Jason-3卫星不符合编辑标准的数据产品多集中在高纬度和浅海地区。

3 研究方法

收起

3.1 确定交叉点精确的位置

卫星测高技术发展早期通常只对独立的一颗卫星进行处理，因此交叉点被认为只是升弧段和降弧段的交点^[17-18]。随着时代的发展，卫星高度计的数量越来越多，多种卫星数据的联合处理成为趋向。对于卫星倾角不同的情况，还要考虑到升弧段与升弧段、降弧段与降弧段也可能存在交叉点。因此本文在此基础上做了改进，不再对轨道的升降进行判断，而是直接判断卫星轨道地面轨迹是否有交叉点。首先直接输入两个卫星轨道，根据轨道的两个端点坐标拟合二元一次线性方程组，公式为

(3)

$ Ax+B=y \text{，} $

式中，A是拟合得到的卫星轨道的系数；B是拟合得到的卫星轨道的截距；x是经度坐标；y是纬度坐标。本文中设定卫星数据的经度范围为0°～360°，其中0°～180°代表东经，180°～360°代表西经；纬度范围为−90°～90°，正值代表北纬，负值代表南纬。通过拟合得到两个卫星轨道的二元一次方程，联立方程组判断是否有解。如果没有解，进一步判断是否是由卫星跳变引起（如果卫星轨迹穿过本初子午线，其星下点经度在0°（360°）附近会产生跳变），判断卫星轨道地面轨迹两个端点经度差的绝对值是否大于180°，如果大于180°，通过减去360°的方法使经度变化连续。然后再一次联立方程组进行判断是否有解，有解的情况下我们认为两个轨道相交。在这个条件下，遍历求出两颗卫星轨道星下点之间距离在1°范围内的星下点。由于两颗卫星轨道星下点之间距离在1°的范围内的星下点数较少，轨道也较短，因此可以视作直线。根据公式（3）输入两个端点坐标，求出两个轨道地面轨迹的二元一次方程，联立方程组得到近似交叉点位置（在1°范围内拟合的直线求出来的交叉点肯定和实际交叉点存在一定距离）。在近似交叉点位置的0.4°、0.6°、0.8°范围内分别再对轨道点拟合直线方程组求交叉点，选取在近似交叉点范围内最近的距离拟合出来的交叉点作为精确交叉点。通过卫星数据产品计算得到的HY-2C卫星两条轨道之间的距离在1°范围内的星下点以及精确交叉点位置如图3所示。

3.2 确定交叉点海面高度异常

在求取交叉点精确位置之后，为了比较在星星交叉定标中使用不同的插值方法是否会影响到卫星高度计的海面高度数据质量结果，本文用3种典型的交叉定标插值方法将SLA插值到交叉点精确位置。

方法一（临近点插值）：在两颗卫星轨道星下点之间距离在1°范围内的所有星下点中选择距离交叉点最近的两个星下点。如图4所示，每条轨道选择距离交叉点最近的两个星下点，然后通过式（4）计算交叉点海面高度异常。

(4)

$ H=ha+\left(hb-ha\right)\times \frac{lat-lata}{latb-lata} \text{，} $

式中，H是不同轨道交叉点插值得到的SLA；ha是纬度较小星下点的SLA；hb是纬度较大星下点的SLA；lat是交叉点的纬度值；lata是纬度较小星下点的纬度值；latb是纬度较大星下点的纬度值。

方法二（最近点插值）：在两颗卫星轨道星下点之间距离在1°范围内的所有星下点中选择距离交叉点最近的一个星下点，如图5所示，将最近星下点的高度异常作为精确交叉点位置的海面高度异常。

方法三（三次样条插值）：为了通过三次样条函数得到更好的插值效果，要保证距离1°范围内的两颗卫星轨道交叉点两侧至少存在4个点以上，如图中蓝色的卫星星下点（图6 ）。然后通过三次样条插值把选取的星下点的海面高度异常插值到精确交叉点位置。

卫星轨道的自交叉是评估和确认测高精度稳定性的重要手段。为了了解在星星交叉定标中使用不同的插值方法是否会影响卫星高度计的海面高度数据质量结果，本文在交叉定标中使用了3种典型插值方法对HY-2C卫星、HY-2B卫星、Jason-3卫星的数据质量进行自交叉点分析。为了减少时间变化对海面高度的影响，卫星交叉点处的最大时间间隔限制为3 d。采用交叉点海面高度异常差的平均值和标准差来对海面高度数据质量进行评估。理想情况下，交叉点海面高度异常差的平均值将接近0，标准差将非常小^{[3, 5]}。这种评估方法已经被用于TOPEX/Poseidon、Jason-1、Jason-2、Jason-3、Sentinel-3A、SARAL和HY-2A等卫星高度计高度产品的精度评估^{[6-10, 12, 15]}。表2至表4是HY-2C卫星、HY-2B卫星和Jason-3卫星进行自交叉点分析的结果。

上述3种交叉定标得到的交叉点数量相近，HY-2C卫星、HY-2B卫星和Jason-3卫星使用三次样条插值可以得到最优的结果。结果表明，HY-2C卫星、HY-2B卫星和Jason-3卫星在星星交叉定标中使用不同的插值方法会影响卫星海面高度数据质量评估的结果，并且本文经过对3种插值方法的比较后得出三次样条插值法相较于前两种方法更适合对卫星数据质量进行分析。为了避免不同的插值方法对HY-2C卫星SGDR产品的海面高度数据质量分析结果造成影响，在后续对卫星的海面高度数据质量进行分析时，将使用三次样条插值获取交叉点。

4 结果与分析

收起

4.1 HY-2C自交叉点分析

通过分析交叉点海面高度差值可以实现对高度计海面高度质量的评估，表2结果显示，对HY-2C卫星高度计SGDR中海面高度异常数据进行自交叉点分析的平均值为0.03 cm，标准差为6.17 cm。这证明了HY-2C卫星高度计观测性能的稳定性。

为了更好地对HY-2C卫星海面高度数据质量进行评估，选取从2021年1月21日至4月20日相同时间段内的HY-2B卫星数据和Jason-3卫星数据进行周期内海平面高度异常差异平均值和标准差的比较，如图7所示。

如图7a所示，HY-2C卫星各个观测周期自交叉点海面高度异常差的平均值在0 cm上下波动，标准差在6 cm上下波动。图7b和图7c是在和HY-2C卫星相同的数据筛选标准和限制条件下得出的结果。图7b显示，HY-2B卫星各个观测周期自交叉点海面高度异常差的平均值在0 cm上下波动，标准差在4 cm上下波动。图7c显示，Jason-3卫星各个观测周期的交叉点海面高度异常差的平均值在0 cm上下波动，标准差在5 cm上下波动。HY-2C卫星与HY-2B卫星、Jason-3卫星的自交叉定标结果相近，卫星间海面高度异常差异自交叉结果的均值在0.02 cm以内，标准差在2 cm以内，这些结果证实HY-2C卫星高度计测高能力与HY-2B卫星、Jason-3卫星的相近。

4.2 HY-2C卫星与HY-2B卫星和Jason-3卫星互交叉点比较分析

对HY-2C卫星和HY-2B卫星、Jason-3卫星的海面高度异常数据进行互交叉点分析的方法和自交叉点分析的计算方法一样。HY-2C卫星和HY-2B卫星、Jason-3卫星互交叉点海面高度异常数据统计结果如表5所示，HY-2C卫星和HY-2B卫星海面高度异常差异的平均值为−0.47 cm，标准差为5.32 cm。HY-2C卫星和Jason-3卫星海面高度异常差异的平均值为−0.3 cm，标准差为5.32 cm。图8a显示了HY-2C卫星和HY-2B卫星在不同观测周期内互交叉点海面高度异常差异的平均值和标准差的变化，平均值在−1～0 cm之间波动，标准差在5 cm上下波动。图8b显示了HY-2C卫星和Jason-3卫星在不同观测周期内互交叉点海面高度异常差异的平均值和标准差的变化，平均值在−1～1 cm之间波动，标准差在5 cm上下波动。结果非常接近Jia等^[7]提供的HY-2B卫星和Jaosn-3卫星互交叉点精度水平。

5 结论

收起

本文使用了HY-2B卫星、Jason-3卫星同期观测数据对我国HY-2C卫星雷达高度计SGDR产品第13～21周期的海面高度数据进行质量分析。HY-2C卫星在全球范围内的有效数据点较多，不符合观测条件的数据产品多集中在高纬度和浅海地区。使用临近点插值、最近点插值、三次样条插值法分别对3颗卫星的海面高度异常数据进行了自交叉点分析，从均值和标准差来看，在星星交叉定标中使用三次样条插值法对HY-2C卫星、HY-2B卫星、Jason-3卫星的定标结果最好。这说明在星星交叉定标中使用合适的插值方法会影响卫星高度计的海面高度数据质量分析结果，此外，HY-2C卫星和HY-2B卫星互交叉点海面高度异常差异的平均值为−0.47 cm，标准差为5.32 cm；HY-2C卫星和Jason-3卫星互交叉点海面高度异常差异的平均值为−0.3 cm，标准差为5.32 cm。这些数据表明HY-2C卫星的测高精度接近Jason-3卫星的数据，说明HY-2C卫星的测高精度与HY-2B卫星和Jason-3卫星一致，均具有较高的精度。本文研究结果还可为HY-2C卫星SGDR的应用提供依据，并为GDR的生产提供参考。

基金

收起

“观澜二号”(GL-2）敏捷编队高度计小卫星工程项目（2021WHZZB1602)
“观澜号”海洋科学卫星校飞试验项目(2021WHZZB1705)
基于人工智能的孪生海洋关键技术与应用示范项目(2022QNLM050301-4)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

刘亚龙. HY-2雷达高度计海面高度定标技术研究[D]. 青岛: 中国海洋大学, 2014.

Liu Yalong. Caliration technology for HY-2 radar altimeter sea surface height[D]. Qingdao: Ocean University of China, 2014.

汪栋. HY-2卫星高度计仪器误差评估与校正研究[D]. 青岛: 中国海洋大学, 2013.

Wang Dong. The research for evaluations and corrections of HY-2 satellite altimeter instrument error[D]. Qingdao: Ocean University of China, 2013.

Sepulveda H H, Queffeulou P, Ardhuin F. Assessment of SARAL/AltiKa wave height measurements relative to buoy, Jason-2, and Cryosat-2 data[J]. Marine Geodesy, 2015, 38(S1): 449−465.

Dorandeu J, Ablain M, Faugère Y, et al. Jason-1 global statistical evaluation and performance assessment: calibration and cross-calibration results[J]. Marine Geodesy, 2004, 27(3/4): 345−372.

Ablain M, Philipps S, Picot N, et al. Jason-2 global statistical assessment and cross-calibration with Jason-1[J]. Marine Geodesy, 2010, 33(S1): 162−185.

Jiang Xingwei, Jia Yongjun, Zhang Youguang. Measurement analyses and evaluations of sea-level heights using the HY-2A satellite’s radar altimeter[J]. Acta Oceanologica Sinica, 2019, 38(11): 134−139.

Jia Yongjun, Yang Jungang, Lin Mingsen, et al. Global assessments of the HY-2B measurements and cross-calibrations with Jason-3[J]. Remote Sensing, 2020, 12(15): 2470.

杨磊, 周兴华, 林明森, 等. HY-2A卫星雷达高度计全球IGDR数据质量评估[J]. 地球物理学进展, 2016, 31(2): 629−636.

Yang Lei, Zhou Xinghua, Lin Mingsen, et al. Global statistical assessment of HY-2A altimeter IGDR data[J]. Progress in Geophysics, 2016, 31(2): 629−636.

刘治中, 杨俊钢, 张杰, 等. 基于与Jason-2数据比对的Jason-3卫星高度计全球数据质量评估[J]. 海洋学报, 2020, 42(3): 129−139.

Liu Zhizhong, Yang Jungang, Zhang Jie, et al. Jason-3 global statistical assessment based on Jason-2[J]. Haiyang Xuebao, 2020, 42(3): 129−139.

Ji Min, Ma Huadong, Li Gang, et al. Quality evaluation of sea level height data of Jason-3 satellite based on Jason-2[C]//Proceedings of the 7th International Conference on Geoinformatics in Sustainable Ecosystem and Society. Guangzhou: Springer, 2019: 103−121.

Wang Zhixiong, Zou Juhong, Stoffelen A, et al. Scatterometer sea surface wind product validation for HY-2C[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14: 6156−6164.

Kong Qiaoli, Guo Jinyun, Sun Yu, et al. Centimeter-level precise orbit determination for the HY-2A satellite using DORIS and SLR tracking data[J]. Acta Geophysica, 2017, 65(1): 1−12.

Wang Jichao, Yu Ting, Deng Fangyu, et al. Acquisition of the wide swath significant wave height from HY-2C through deep learning[J]. Remote Sensing, 2021, 13(21): 4425.

彭海龙, 林明森, 穆博, 等. HY-2A卫星雷达高度计数据的全球统计评价及质量分析[J]. 海洋学报, 2015, 37(7): 54−66.

Peng Hailong, Lin Mingsen, Mu Bo, et al. Global statistical evaluation and performance analysis of HY-2A satellite radar altimeter data[J]. Haiyang Xuebao, 2015, 37(7): 54−66.

Zaouche G, Perbos J, Lafon T, et al. OSTM/Jason-2: assessment of the system performances (ocean surface topography mission: OSTM)[J]. Marine Geodesy, 2010, 33(S1): 26−52.

Wang Jinbo, Lee L F, Qiu Bo, et al. An observing system simulation experiment for the calibration and validation of the surface water ocean topography sea surface height measurement using in situ platforms[J]. Journal of Atmospheric and Oceanic Technology, 2018, 35(2): 281−297.

叶小敏, 林明森, 宋庆君. 基于现场观测数据的卫星雷达高度计海面风速和有效波高真实性检验方法研究[J]. 遥感技术与应用, 2014, 29(1): 26−32.

Ye Xiaomin, Lin Mingsen, Song Qingjun. Research on validation methods of satellite radar altimeter retrieved sea surface wind speed and significant wave height using in-situ data[J]. Remote Sensing Technology and Application, 2014, 29(1): 26−32.

杨磊, 周兴华, 彭海龙, 等. 基于Jason-2的Saral/AltiKa高度计全球统计评估与交叉定标[J]. 海洋科学进展, 2014, 32(4): 482−490.

Yang Lei, Zhou Xinghua, Peng Hailong, et al. Global assessment and cross-calibration of Saral/AltiKa based on Jason-2 altimeter[J]. Advances in Marine Science, 2014, 32(4): 482−490.

2022年第44卷第10期

PDF下载

175

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2022168

接收时间：2022-05-04
首发时间：2026-02-01
出版时间：2022-10-01

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2022-05-04
修回日期：2022-06-06

基金

“观澜二号”(GL-2）敏捷编队高度计小卫星工程项目（2021WHZZB1602)

“观澜号”海洋科学卫星校飞试验项目(2021WHZZB1705)

基于人工智能的孪生海洋关键技术与应用示范项目(2022QNLM050301-4)

作者信息

1.中国海洋大学信息科学与工程学部海洋技术学院，山东青岛 266100

2.青岛海洋科学与技术试点国家实验室区域海洋学与数值模拟实验室，山东青岛 266237

通讯作者:

于方杰，男，山东省青岛市人，教授，主要从事卫星−无人移动平台协同海洋观测技术研究。E-mail：yufangjie@ouc.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2022168

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

观测量	最小值	最大值
20 Hz测距数目	10	−
20 Hz测距标准差/mm	0	200
轨道高度−测距值/mm	−130 000	100 000
模型干对流层校正/mm	−2 500	−1 900
辐射计湿对流层校正/mm	−500	−1
Ku波段电离层校正/mm	−400	40
Ku波段海况偏差/mm	−500	0
潮汐校正/mm	−5 000	5 000
固体潮校正/mm	−1 000	1 000
极潮校正/mm	−150	150
Ku波段有效波高/mm	0	11 000
Ku波段向后散射系数/dB	7	30
风速/(m·s⁻¹)	−0	30
Ku波段偏指向角平方/（°）²	−0.2	0.64
后向散射系数标准差/dB	−	1
后向散射系数观测数	10	−

观测量

最小值

最大值

20 Hz测距数目

−

20 Hz测距标准差/mm

200

轨道高度−测距值/mm

−130 000

100 000

模型干对流层校正/mm

−2 500

−1 900

辐射计湿对流层校正/mm

−500

−1

Ku波段电离层校正/mm

−400

Ku波段海况偏差/mm

−500

潮汐校正/mm

−5 000

5 000

固体潮校正/mm

−1 000

1 000

极潮校正/mm

−150

150

Ku波段有效波高/mm

11 000

Ku波段向后散射系数/dB

风速/(m·s⁻¹)

−0

Ku波段偏指向角平方/（°）²

−0.2

0.64

后向散射系数标准差/dB

−

后向散射系数观测数

−

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	14 746	14 699	14 312
平均值/cm	0.03	0.03	0.03
标准差/cm	6.63	6.89	6.17

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

14 746

14 699

14 312

平均值/cm

0.03

标准差/cm

6.63

6.89

6.17

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	5 005	5 045	4 936
平均值/cm	0.01	0	0.01
标准差/cm	4.60	5.1	4.21

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

5 005

5 045

4 936

平均值/cm

0.01

标准差/cm

4.60

5.1

4.21

	临近点插值	最近点插值	三次样条插值
交叉点数量（SLA<20 cm）	22 639	22 610	22 233
平均值/cm	−0.01	−0.01	0
标准差/cm	6.63	5.45	4.70

临近点插值

最近点插值

三次样条插值

交叉点数量（SLA<20 cm）

22 639

22 610

22 233

平均值/cm

−0.01

标准差/cm

6.63

5.45

4.70

	HY-2C和HY-2B	HY-2C和Jason-3
交叉点数（SLA<20 cm）	15 470	15 799
平均值/cm	−0.47	−0.3
标准差/cm	5.32	5.32

HY-2C和HY-2B

HY-2C和Jason-3

交叉点数（SLA<20 cm）

15 470

15 799

平均值/cm

−0.47

−0.3

标准差/cm

5.32