CFOSAT散射计几何定位验证与精度初步评估
1.
2.
3.
4.
Verification of CFOSAT Scatterometer Geolocation and Preliminary Evaluation of Geolocation Accuracy
1.
2.
3.
4.
收稿日期: 2021-01-25 修回日期: 2022-07-08
基金资助: |
|
Received: 2021-01-25 Revised: 2022-07-08
作者简介 About authors
云日升(1974-),男,内蒙古和林格尔人,副研究员,主要从事微波遥感系统仿真与数据预处理研究E⁃mail:
关键词:
Keywords:
本文引用格式
云日升, 刁宁辉, 崔海英, 高畅畅.
Yun Risheng, Diao Ninghui, Cui Haiying, Gao Changchang.
1 引 言
2018年10月29日8时,中法海洋卫星在酒泉卫星发射中心用长征二号丙运载火箭成功发射。中法海洋卫星(CFOSAT)是中法两国合作研制的首颗卫星,中方负责提供卫星平台、海风观测载荷以及发射测控,法方负责提供海浪观测载荷,卫星数据双方共享。CFOSAT可获取全球海面风场、海浪谱以及有效波高等海洋动力参数,有助于更好地认识和掌握海洋动力过程及其变化规律,为海洋预报提供海浪谱、海面风场等多参量的初始场信息,改进海况预报及同化模型,提高对热带风暴、风暴潮等灾害性海况预报的精度与时效[1]。通过对海浪、海面风场进行同步监测,可为海洋科学研究、全球气候变化提供实测数据并且积累长时间序列历史数据。
扇形波束旋转扫描散射计几何定位是CFOSAT散射计数据预处理的难点之一。与笔形波束扫描散射计固定俯仰角不同,其俯仰角有效范围约为26°—46°,同时Doppler频率随星下点经纬度、天线方位扫描角和俯仰角变化,使得全去斜处理后的回波频率同地面定位之间存在非线性对应关系,从而给观测波束条带几何定位带来了困难[5]。CFOSAT几何定位算法首先建立了脉冲回波频率与地面经纬度的基准对应模型,在此基础上对观测条带中心回波频率进行插值获得当前观测条带的地面经纬度。STK(Satellite Tool Kit)卫星工具包软件是美国Analytical Graphics公司开发的航天工业分析软件,STK的核心能力是产生卫星状态向量(位置和速度)和姿态数据、获取时间、遥感器覆盖分析、特定观测几何下目标定位等[6]。卫星散射计通常采用海岸线匹配的方法对几何定位算法进行验证和分析[7],与之相比较相比较,STK用于几何定位算法验证具有精度高、数据适应性更好等优势。但由于卫星微波遥感器系统体制和观测模式的特殊性,将STK用于几何定位验证时需要合理设计验证方案。
首先简要介绍了CFOSAT散射计及几何定位算法。基于CFOSAT散射计实际观测数据,根据CFOSAT散射计扇形波束扫描,地面足印内各条带对应的观测俯仰角不确定的特点,设计了CFOSAT散射计STK几何定位验证方案,采用STK对CFOSAT散射计几何定位原型算法、包含天线安装误差的完整工程算法进行了详细验证。由于天线在轨可能存在畸变等影响,实际定位精度要低于验证分析的结果。本文在算法验证的基础上,结合有源定标器实验数据对定位参数进行了校正,最终获得CFOSAT散射计几何定位偏差H极化约为5.91 km,V极化约为6.06 km。
2 CFOSAT散射计几何定位算法
CFOSAT散射计是国际上首个扇形波束旋转扫描散射计,工作在Ku波段,卫星轨道高度520 km,测绘带刈幅约1 000 km,地球表面日覆盖率约为90%。表1为CFOSAT散射计主要系统参数和轨道参数。
表1 CFOSAT散射计主要工作参数
Table1
系统参数 | 轨道参数 | ||
---|---|---|---|
载频/GHz 发射功率/W 脉冲宽度/ms 脉冲重复频率/Hz 俯仰角范围/° 极化方式 | 13.256 120 1.35 1.35 26—46 H/V | 半长轴/km 偏心率 近地点幅角/° 倾角/° 升交点赤经/° 降交点地方时 | 6 891.984 0.001 23 90 97.528 1 322.67 7:00 AM |
图1
对于CFOSAT扇形波束扫描散射计,Doppler频率随星下点经纬度、天线扫描方位角和俯仰角变化。特别是俯仰向观测范围约20°,对应的波束近端和远端的Doppler频率差约300 kHz左右[8],为了保持回波频率在零频附近,利用Doppler预补偿LUT(Look-Up Table)对发射频率进行频率预偏移。由于Doppler频率的影响,前视观测时回波信号的有效带宽被压缩。因此,系统对发射Chirp信号采取前视观测时变调频斜率符号的设计,这样可以使得天扫描前视观测和后视观测时频率分辨率基本一致。当后视观测时,回波频率为:
当前视观测时,回波频率为:
其中:
3 CFOSAT散射计几何定位验证
采用STK对CFOSAT散射计几何定位算法进行验证。图2为几何定位验证流程。作为商业仿真分析软件,STK暂未设定特定的外部星历和姿态数据导入接口。因此将CFOSAT某一天的TLE(Two Line Element)两行报文件和STK配置的姿态模型作为STK的输入,按照散射计观测时间、脉冲时序生成与卫星散射计相匹配的卫星星历和姿态数据。然后在确保STK和散射计在卫星星历、姿态数据、极化观测和天线转速等参数完全一致的基础上,分别获得考察时间点的STK和散射计的几何定位结果,进行验证分析。在验证层次上分别进行CFOSAT散射计几何定位原型算法验证和包含天线安装误差的完整几何定位算法验证;在验证方法上分别在给定时间段的固定俯仰角上进行多个脉冲统计比较验证和给定时刻CFOSAT散射计回波脉冲40个条带对应俯仰角上进行单个脉冲条带几何定位比较验证,后者需要先反推出散射计脉冲40个条带对应的俯仰角,然后再将俯仰角信息代回到STK中计算地面几何定位结果。验证过程还分别针对观测前视和后视、升轨和降轨等不同散射计工作条件进行了深入的验证分析。
图2
图2
利用STK对CFOSAT散射计几何定位验证流程
Fig.2
Verification process of CFOSAT scatterometer geolocation using STK
3.1 几何定位原型算法验证
首先利用STK对CFOSAT散射计几何定位原型算法进行验证。对CFOSAT原型算法的验证主要用以明确CFOSAT散射计条带几何定位处理中的星历和姿态数据的使用、坐标系转换算法、Doppler补偿算法和关于回波频率和固定入射角序列的几何定位基准模型算法。图3为几何定位原型算法的多脉冲验证(显示脉冲有抽取)。图3左为时间段(UTC):20181218T150000-20181218T151000的散射计若干脉冲几何定位验证,对应于卫星升轨;图3右为时间段(UTC):20181218T153000 -20181218T154000的散射计若干脉冲几何定位验证,对应于卫星降轨。不同极化的几何定位差异体现在图1天线坐标系到星体坐标系的转换上。对于CFOSAT散射计几何定位原型算法,不需要考虑不同极化情况下的几何定位验证。图3中STK定位结果是相对于俯仰角20°—50°,间隔5°的结果。由图3可知,散射计几何定位原型算法的结果与STK有很好的一致性。
将考察的升降轨两个时间段的各45 000个脉冲的定位偏差进行统计,结果如图4所示。与STK相比较,升轨与降轨时间段的几何定位偏差没有明显差异。同时在20°—50°俯仰角范围内也基本一致,平均定位偏差约36 m。
图3
图3
几何定位原型算法多脉冲验证
Fig.3
Multi-pulse verification of geolocation prototype algorithm
图4
图4
CFOSAT散射计几何定位原型算法定位偏差
Fig.4
Bias between CFOSAT scatterometer geolocation prototype algorithm and STK
3.2 完整几何定位算法验证
CFOSAT散射计采用两个相互平行的波导缝隙天线实现H/V极化的扫描观测,观测几何如图5所示。两个天线观测方位角设计为180°间隔,散射计交替发射并接收H极化和V极化脉冲。天线阵列坐标观测视向量首先通过天线安装矩阵变换到星体坐标系下。星体坐标系坐标原点在航天器质心,3个坐标轴分别沿航天器3个惯性主轴方向。卫星姿态决定了星体坐标系和轨道坐标系之间的关系。轨道坐标系坐标原点为航天器的质心,Z轴指向地心(反径向),X轴在轨道平面内垂直于Z轴,指向卫星运动方向,Y轴与X轴和Z轴构成右手正交坐标系统。散射计观测方位角是天线观测方向与卫星运行方向的夹角。
CFOSAT散射计H极化天线安装矩阵:
V极化天线安装矩阵:
图5
图5
CFOSAT散射计天线阵面观测几何
Fig.5
Observation geometry of CFOSAT scatterometer antenna
天线安装矩阵是星体坐标系统和散射计天线坐标系3个坐标轴之间的实测角度(单位:度),反映了天线在地面安装时的偏差。特别地,对于扇形波束旋转扫描散射计,不但需要考察天线在方位向上相对于星体坐标系方位角零位的固有偏差,而且需要考虑由于天线相对于旋转轴的安装偏差所引起的在扇形波束俯仰向上方位角偏差的分布。因此,安装矩阵揭示了由于天线安装偏差,观测方位角与实际方位角的偏差是方位向零位固有偏差和随俯仰角变化的方位角偏差二者的叠加。由安装矩阵计算得到H极化方位向零位固有偏差为-2.075 7°,V极化方位向零位固有偏差为-1.487 1°,两个天线方位角实际间隔为179.411 4°。图6给出了由安装矩阵得到的随俯仰角变化的天线方位角偏差值。对于特定观测脉冲,由于天线安装误差的存在,在整个观测俯仰角范围内,天线单元给出的方位角与实际方位角之间存在一定的偏差。特别是越接近波束近端,误差越大。如在20°俯仰角附近,天线方位角误差约±0.5度。因此,CFOSAT散射计完整的几何定位算法,需要对由天线安装误差带来的方位角固有零位偏差和随俯仰角变化的方位角偏差进行补偿,才能获得观测脉冲各条带的准确几何定位。
图6
图6
CFOSAT散射计天线方位角偏差随俯仰角变化
Fig.6
CFOSAT scatterometer antenna azimuth bias vs. elevation angle
表2 STK随俯仰角变化的方位角补偿值
Table2
STK俯仰角 /° | STK方位角补偿 H-极化/° | STK方位角补偿 V-极化/° |
---|---|---|
20 25 30 35 40 45 50 | 0.438 1 0.281 6 0.174 6 0.095 7 0.034 4 -0.015 4 -0.057 2 | -0.400 0 -0.247 3 -0.142 9 -0.065 9 -0.006 0 0.042 5 0.083 3 |
表2为对20°—50°天线俯仰角(5°间隔),分别对H极化和V极化给出的方位角补偿值。采用图2给出的多脉冲几何定位验证方法,散射计采用图1所示的完整几何定位处理。由于实际观测存在天线安装偏差,因此对STK的观测方位角进行固定方位角零位偏差补偿和表2中随俯仰角的方位角偏差补偿,从而对CFOSAT散射计完整几何定位算法进行验证。由于上节已对散射计几何定位原型算法进行了验证,这一验证过程可以验证CFOSAT扇形波束扫描散射计天线安装矩阵分析和使用是否正确,以及完整的坐标系变换和几何定位算法的准确性。多脉冲几何定位算法验证结果如图7所示。几何定位偏差与图4原型几何定位算法相比较结果相差不大,在输入星历、姿态、极化等几何定位参量一致的前提下,CFOSAT散射计完整几何定位算法与通过STK获得的几何定位结果相差约数十米。对于应用天线安装矩阵的CFOSAT散射计完整几何定位算法,几何定位偏差随俯仰角有增大的趋势。
图7
图7
CFOSAT散射计几何定位算法与STK定位结果偏差
Fig.7
Bias between CFOSAT scatterometer geometric positioning algorithm and STK results
由图可知,升轨和降轨条件下,不同极化和不同视向观测的几何定位结果差别不大,其平均偏差一致地在100 m以下,最大偏差接近150 m。几何定位偏差有随俯仰角略有增大的趋势。需要说明的是:所谓单脉冲验证,是首先计算散射计特定脉冲几何定位各条带的观测俯仰角,然后将各条带的实际俯仰角代回到STK中获得STK在此观测俯仰角下的几何定位结果。实际评估中,针对不同极化、观测视向和升降轨均要对考察时间段内的多个脉冲进行分析验证,并非只是在单一脉冲上的比较。
图8
图8
CFOSAT散射计条带几何定位与STK定位结果偏差
Fig.8
Bias between CFOSAT scatterometer slice geolocation algorithm and STK results
4 CFOSAT散射计数据几何定位精度初步评估
通过对CFOSAT散射计几何定位原型算法和完整几何定位算法的STK验证,表明CFOSAT散射计几何定位算法的正确性,对天线安装矩阵的使用也是合理的。同时,验证结果表明:CFOSAT散射计几何定位算法与STK几何定位算法偏差在100 m以内。为了确保验证的合理性,整个验证过程CFOSAT散射计几何定位算法输入的星历和姿态数据均采用STK仿真得到的数据,并非实际卫星给出的数据。同时,由于卫星在轨运行的复杂性和不确定性,散射计天线在轨可能存在畸变,其安装特性可能与地面发射前存在一定的差异。因此,在当前几何定位参数条件下,实际观测数据的几何定位误差远大于100 m。在验证了几何定位算法正确性的前提下,需要对CFOSAT散射计实际观测数据的几何定位参数进行优化校正,并对实际观测数据的定位精度进行评估。
采用2019年7月国家卫星海洋应用中心组织的在内蒙古新巴尔虎左旗嵯岗进行的有源定标器外定标实验相关数据对CFOSAT实际数据几何定位结果进行估计。有源定标器发射的点目标信号落在CFOSAT散射计数据特定脉冲和条带内。由于条带分辨率约10 km左右,因此从散射计数据中获得有源定标器的位置精度受到一定限制。排除数据缺失或异常、有源定标器信号落在边缘条带等特殊数据,从27轨包含有源定标器散射计数据中筛选出7轨数据用于几何定位算法精度评估和误差校正。以后向散射系数最大值条带位置为中心,选择若干相邻脉冲和条带进行二维数据拟合,搜索得到理论上更接近于有源定标器的真实位置。在此基础上,目前主要通过调整观测方位角参数,之后与有源定标器地面确定位置相比较,搜索最优定位结果对应的方位角。
表3 利用外定标数据估计散射计定位精度
Table3
轨道号 | 升/降轨 | H-极化/m | V-极化/m |
---|---|---|---|
3406 3428 3436 3497 3519 3625 3633 | D A D D A A D | 5 970.67 4 847.42 3 838.96 7 913.61 7 634.05 4 512.14 6 678.34 | 7 414.05 4 656.31 4 789.44 4 656.94 6 257.86 6 898.44 7 768.24 |
平均值 | 5 913.59 | 6 063.04 |
表3为利用有源定标器优化之后的7轨实际数据几何定位结果。其中“A”表示升轨,“D”表示降轨。第三列为H极化散射计几何定位结果与有源定标器地面实际位置之间的偏差值;第4列为V极化散射计几何定位结果与有源定标器地面实际位置之间的偏差值。H极化几何定位的平均偏差为5.91 km,V极化几何定位的平均偏差为6.06 km。由此可见,在STK对CFOSAT散射计几何定位算法验证之后,通过有源定标器数据的进一步校正,CFOSAT散射计实际数据几何定位精度在整体上约为6 km左右。
图9
图9
CFOSAT散射计后向散射系数数据海岸线匹配
Fig.9
Coastline matching of backscattering coefficient data of CFOSAT Scatterometer
5 结 论
中法海洋卫星(CFOSAT)是中法两国合作研制的首颗卫星,其上搭载的扇形波束扫描散射计是中方提供的唯一载荷,是我国首次从方案设计、载荷技术和数据处理全链条自主创新实现的卫星散射计。CFOSAT散射计科学数据主要应用于海面风场的实时观测,也可应用于极冰监测、植被分析等领域。本文介绍了CFOSAT扇形波束扫描散射计几何定位算法之后,主要采用STK软件对CFOSAT散射计几何定位算法进行了验证。利用STK进行CFOSAT散射计几何定位算法验证的基本流程是在确保CFOSAT几何定位算法模块和STK输入参数一致的条件下,对两者各自独立获得的地面定位结果进行比较。针对CFOSAT散射计扇形波束扫描观测体制,分别就扇形波束扫描散射计几何定位原型算法、包含天线安装误差的完整几何定位算法利用STK进行了验证。通过对CFOSAT散射计H极化和V极化、升轨和降轨、前视观测和后视观测数据的几何定位进行验证,其结果表明了CFOSAT扇形波束散射计几何定位算法的准确性,同时验证了扇形波束散射计天线安装矩阵使用的正确性。CFOSAT几何定位算法定位结果与STK定位结果平均偏差在100 m以内,不同极化、升降轨和扫描观测方位条件下差别不大,几何定位偏差有随俯仰角有增大的趋势。由于可能存在的天线在轨畸变等原因,将地面测试的天线安装矩阵应用于天线坐标系到星体坐标系的转换之后,实际数据的几何定位精度远大于100 m。结合有源定标器外定标实验数据,进一步校正了几何定位参数,给出了CFOSAT散射计几何定位精度的初步评估,得到CFOSAT散射计实际数据的几何定位精度约6 km左右。下一步需要对实际数据几何定位误差来源进行深入分析,并对几何定位算法进一步优化,有望获得更为准确的几何定位结果。
参考文献
SWIM: The first spaceborne wave scatterometer
[J]. ,
An analysis of a rotating, range-gated, fanbeam spaceborne scatterometer concept
[J]. ,
Orbit performances validation for CFOSAT scatterometer
[J]. ,
中法海洋卫星微波散射计在轨性能验证
,,
First results from the rotating fan beam scatterometer onboard CFOSAT
[J]. ,
The processing and simulation of the CFOSAT RFSCAT
[C]∥
Satellite orbit calculation based on SGP4 model
[J]. ,
基于SGP4模型的卫星轨道计算
,,
Estimation and correction of geolocation errors of the CFOSAT scatterometer using coastline backscatter coefficients
[J]. ,
Simulation and retrieval of wind of CFOSAT rotating-fan beam SCATterometer
[C]∥
CFOSAT Scatterometer Data Level-1 Processing and Preliminary Results
[C]∥
A perspective on the performance of the CFOSAT fotating fan-beam scatterometer
[J]. ,
/
〈 | 〉 |