高光谱观测卫星高精时间同步方案设计与应用.pdf

摘要：为了满足星上高光谱相机、全谱段光谱成像仪、星敏感器等时间精度敏感仪器的高精度时间同步需求，高光谱观测卫星上设计了可适应多种精度需求的时间同步方案。对于系统时间，综合通过星地时差集中校时、均匀校时、GPS总线校时、GPS硬件秒脉冲校时等校时方法进行时间修正；对于终端用户，分别给出直接采用广播的进行本地时间修正的时间同步方法、采用广播及内部秒脉冲的时间同步方法、采用广播授时和不同秒脉冲校时方案的地面时间重建算法。利用在轨数据，对卫星星敏感器、全谱段光谱成像仪及可见短波红外高光谱相机的曝光时刻时间重建方法进行了应用说明，并对不同方法的时统误差进行了分析。关键词：高光谱观测卫星；时间同步；GPS；秒脉冲中图分类号：V448文献标识码：A0 引言高光谱遥感技术具有观测通道多、光谱分辨率高等特点，是遥感技术的重要发展方向1-2。高光谱遥感技术已经被广泛应用于环境综合监测、地矿调查、大气成分探测等诸多领域3-10。高光谱观测卫星（GF-5）卫星是我国第1颗高光谱综合观测卫星，设计运行于太阳同步轨道，轨道高度705 km，主要用于获取从紫外到长波红外谱段的高光谱分辨率遥感数据产品，是我国高分专项的重要组成部分，是体现国家高光谱分辨率对地观测能力的重要标志11-12。GF-5配置了可见短波红外高光谱相机、全谱段光谱成像仪、大气环境红外甚高光谱分辨率探测仪、大气主要温室气体监测仪、大气痕量气体差分吸收光谱仪和大气气溶胶多角度偏振探测仪。这6台载荷，具有高光谱、太阳掩星、耀斑观测、红外、偏振、多角度等多种观测手段13。为满足卫星高精度姿态控制、有效载荷高图像定位精度和配准精度等需求，要求卫星具有高精度时间同步精度14-15。本文设计了适用于GF-5的高精度时间同步系统方案，综合通过星地时差集中校时、GPS总线校时、GPS秒脉冲校时、校匀校时、总线广播授时等多种校时方案，同时适应一般精度需求和高精度时间同步应用需求。针对星敏感器、可见短波红外高光谱相机及全谱段光谱成像仪这3类对曝光时刻时间同步精度要求最高的远程终端，设计了3种高精度时间同步方案，并利用在轨遥测数据进行验证。1 卫星时间同步总体方案数管计算机（Central Terminal Unit，CTU）的时间是整星的系统时钟，卫星的姿态数据、遥感数据、广播消息和广播时间等均采用系统时间进行描述。系统时间采用UTC时间基准，以48位0.1 ms累计数作为星上的计时方式，数管周期性通过1553B总线将当前系统时钟广播给各终端（Remote Terminal，RT）用户使用。为了获得高精度的秒内时间精度和高稳定的时钟周期，数管计算机采用高稳恒温晶振触发产生16位长度、分辨率为0.1 ms的时钟计数供软件读取。数管每一节拍开始时，通过读取该时钟计数，与上一节拍读取的计数值获取差值T，累加到系统软件时钟上。数管通过1553B总线进行广播的频率为1 Hz，在广播前，提前把广播产生的时延量进行扣除。1553B总线时间广播的准确度受软件处理时延和总线通信传输时延估算精度和稳定度的影响，仅用广播时间进行校时的用户存在ms级的误差。时间误差将会对星敏感器的姿态测量和高分辨率相机的图像定位带来影响。表1给出不同时间精度对星敏感器及高分辨率相机的影响分析。可看出：对于星敏感器，当时间误差在5 ms以内时，星敏的指向误差约小于1，可以满足姿态解算需求；对于高分辨率相机，当时间误差在5 ms左右时，对定位造成的误差就达37.5 m，远不能满足高分辨率图像定位的要求。收稿日期：2022-07-08；修回日期：2022-12-27作者简介：温渊（1987），男，高工，硕士；李云端（1976），男，研究员，硕士（通信作者）。海 军 航 空 大 学 学 报海 军 航 空 大 学 学 报第38卷表1 不同时间精度对星敏感器及高分辨率相机的影响分析Tab.1 Influence of different time accuracytowards star tracker and camera序号1234567时间误差/ms0.11234510星敏指向误差/()0.021 60.216 00.432 00.648 00.864 01.080 02.160 0载荷定位误差/m0.757.5015.0022.5030.0037.5075.00我国现有的常规遥感卫星校时方案中：针对时间精度需求在ms级的一般精度需求的终端，通过数管1553B广播等手段进行时间同步16；针对时间精度需求在s级的终端（如高分辨率相机），可通过高精度时间基准单元进行高精度时间同步17。GF-5同时存在了多种一般精度需求的终端及多种高精度需求的终端，为了简化系统的设计，区别于其他高分辨率遥感卫星，在不引入额外的时间基准单元的情况下：对于一般精度需求的终端，仅采用1553B广播同步的时统方案；对于高精度需求的终端，设计通过数管广播结合GPS秒脉冲的时统方案，以满足高精度时间同步需求18。GF-5时间同步架构示意图如图1所示。图1 GF-5时间同步架构示意图Fig.1 Time synchronization framework of GF-5各终端用户根据自身需求，通过数管广播获取1553B接收时刻的时间，根据不同的精度需求，通过内部的计数器等获取进一步细化终端时间。对于终端分系统下位机以外的其他单机，可以采用内部秒脉冲等校时方式，如星敏感器由于没有直接连接到1553B总线上，设计通过内部秒脉冲与姿轨控计算机进行时间同步，以避免姿轨控计算机至星敏感器之间再引入更大的时间误差。已知秒脉冲本身的时间精度优于1 s，对时间精度要求较高的可见短波红外高光谱相机及全谱段光谱成像仪设计采用GPS秒脉冲进行时间同步，通过1553B获取整秒时间，通过内部计数器结合秒脉冲信号获取秒内时间。2 具有多级故障重构功能的系统校时方案数管计算机时间是整星时间的基准，其时间精度直接影响整星各终端的时间。数管计算机虽然采用了高稳恒温晶振，但其精度和稳定度仍有限，长期积累也会发生时钟漂移。为此，卫星设计采用软件时钟，通过具有多级故障重构功能的校时方案进行校时，优先级从高到低依次为GPS秒脉冲校时、GPS总线校时、高稳晶振均匀校时、CPU计时器均匀校时及软钟校时等，确保在GPS秒脉冲故障、GPS通信故障、高稳晶振故障、CPU计时器故障等情况下，卫星系统时钟误差仍可控。为了进一步提高可靠性，系统时钟在内存中同时维护3份。2.1 GPS秒脉冲校时GPS 接收机每 1 s 通过 RS422 输出脉冲宽度为1 ms的秒脉冲信号，信号形式为宽度1 ms的负脉冲，时间精度可达1 s。CTU在接收到GPS秒脉冲信号后，自动锁存时钟计数器的值，软件周期性计算对应GPS 秒脉冲时刻的数管系统时间，通过取模运算（Mod）获取整秒的余数，即与整秒的偏差值，用其修正本地时钟。数管GPS秒脉冲校时的流程如图2所示。实际应用中，需要考虑到16位内部计数器的进位，采用规范化处理的公式计算数管时间到秒脉冲时刻的间隔dT，计算公式如式（1）所示：dT=mod()TPulse-TLast+5 000,216-5 000。（1）式（1）中：TPulse为秒脉冲时刻锁存计数值；TLast为上一次数管时间；mod()为求余运算。数管递推到秒脉冲时刻在秒内的偏差TMod可由mod()T+dT,10 000计算。但由于MOD公式的特性，计算得到的TMod结果在数值上可能为靠近0的正数（图2所示的TMod1），或者为在数值上接近但小于10000的数（图2所示的TMod2），须对2种情况进行规范化处理，计算公式如式（2）所示：TMod=mod()T+dT+5 000,10 000-5 000。（2）结合式（1）（2），得到有修正系统时间T的计算公式：T=T-TMod=T-modT+mod()TPulse-TLast+5 000,216,10 000-5 000。（3）1 5 5 31 5 5 31 5 5 31 5 5 3 BBBB 总线总线总线总线G P S接收机数管计算机A C E全谱段光谱成像仪获取G P S时间输出G P S时间可见短波红外高光谱相机秒脉冲获取数管时间获取数管时间秒脉冲其他4 台大气载荷获取数管时间G P S-C T U校时数据姿轨控计算机星敏感器获取数管时间内部秒脉冲R T S 信号 248第3期温渊，等：高光谱观测卫星高精时间同步方案设计与应用图2 CTU秒脉冲校时流程图Fig.2 CTU pulse per second time adjusting scheme图 2 中：dT为数管时间到秒脉冲时刻的间隔；TPulse为秒脉冲时刻锁存计数值；TLast为上一次数管时间；TMod为数管递推到秒脉冲时刻的整秒偏差，修正后系统时间T=T-TMod。通过对在轨的遥测数据进行分析可知，在轨秒脉冲时差测量受限于时钟量化精度为0.1 ms，误差稳定在0.1 ms以内，如图3所示。图3 CTU在轨GPS秒脉冲时差测量（校正）值Fig.3 CTU in-orbit GPS pulse persecond time adjusting scheme2.2 GPS总线校时当GPS秒脉冲信号故障时，数管计算机自动降级到GPS总线校时。数管计算机每秒通过1553B总线向各终端用户广播已修正软件时延和通信时延的当前系统时钟。GPS接收机收到数管广播时钟后，在接收机内部进行GPS-CTU时差的计算，把计算差值通过总线消息返回给数管计算机。数管计算机收到该消息后，在判断正确与否及地面授权的情况下，将该时差修正到当前系统时间中。GPS总线校时的优先级仅次于GPS秒脉冲校时，GPS总线时差在轨实测精度接近时间量化精度，在-0.10.2 ms，如图4所示。图4 CTU在轨GPS总线时差测量精度Fig.4 CTU in-orbit GPS bus timing gap measuring accuracy2.3 高稳晶振均匀校时在默认的情况下，数管计算机采用精度达10-8、稳定度达10-9的高稳恒温晶振产生的100 s的时钟计数作为星上时间的基准。地面根据星上下传的遥测帧中的数管时钟与地面时钟进行比对，扣除设备固定时延和轨道空间时延，可计算星地时差。地面根据一段时间的星地时差变化情况，评估生成集中校时注数对星上时钟进行整体修正.对时差的线性变化部分，生成均匀校时注数，即周期性地进行时钟拨快或拨慢1个时间单位，从而实现对星上时钟的漂移情况修正。根据时钟的精度和稳定度，使用该方法，每天1次星地时差测量与校正，可以保证卫星星地时差优于5 ms的要求。利用高稳晶振进行星地时差测量后，进行集中及均匀校时的校时流程如图5所示。2.4 CPU计时器均匀校时星上设计了计时器切换的硬件结构，可通过地面遥控注数选择采用高稳晶振计时器或CPU晶振计时器。CPU自带的晶振计时器精度不如高稳晶振计时器，默认选择高稳晶振计时器。通过更频繁地星地时0.150.100.050-0.05-0.10-0.15秒脉冲时差测量值/ms07-01 07-02 07-03 07-04 07-05 07-06 07-07 07-08 07-09 07-10日期0.250.200.150.100.050-0.05-0.10-0.15GPS时间与CTU时间差/ms07-01 07-02 07-03 07-04 07-05 07-06 07-07 07-08 07-09 07-10日期 249海 军 航 空 大 学 学 报海 军 航 空 大 学 学 报第38卷差测量分析及集中校时，在高稳晶振失效后，也可以保证卫星的星地时差优于 5 ms。星上高稳晶振和CPU晶振的硬件切换结构，如图6所示。图5 高稳晶振均匀校时流程Fig.5 Time Synchronization scheme of high stabilitycrystal homogeneous adjustment图6 星上高稳晶振和CPU晶振的硬件切换结构Fig.6 Hardware switch structure of high stabilitycrystal and CPU crystal2.5 软钟校对星上设计是在硬件计数器均不可用的情况下，采用的1种软件节拍累加计时手段，即每个0.5 s软件周期，直接给系统时间加上0.5 s，再结合均匀校时手段，实现具有一定精度的最小配置、最简时间生成手段。3 适应不同精度需求的终端时间同步方案GF-5 的终端时间同步方案基于 1553B 广播及GPS秒脉冲，不采用时间基准单元。一般精度需求下直接采用1553B总线时间广播结合内部计数器获得本地时间，终端内部的其他单机再设计内部秒脉冲等方式进行时间同步；高精度需求下，设计通过1553B广播结合GPS秒脉冲进行时间获取。下面分别给出星上4种典型的时间获取方式。3.1 基于1553B中断的时间同步方案姿轨控计算机、低分辨率观测载荷等终端用户设计采用基于1553B中断的时间同步方案。终端在本地维护1个本地时间，收到1553B广播中断后，对本地时间的偏差进行计算，以便进行时间误差修正。终端本地具有1个时钟计数器，当终端下位机软件每收到1次数管广播（中断）时，可先通过判断前后2次数管广播时钟的差值等状态，判断数管时钟是否在正常范围。如果正常，在1553B时钟广播触发的中断内，计算终端系统时钟与1553B广播时钟的钟差T，作为校时量使用。在终端下位机软件的时间周期任务开始时，将终端系统钟加上钟差T获得校正后的时间。基于1553B中断的时间同步方案如图7所示。图7 基于1553B中断的终端时钟方案Fig.7 Time synchronization schemebased on 1553B interruption图 7 中：TRT为 终 端 上 一 拍 的 时 间；TRT=TRT+mod()Nnow-Nlast,M n+T；Tctu为终端收到的数管时间，在广播时已经对运算及总线传输时延进行了修正；Nlast为终端上一拍的内部计数值；Nint为收到数管时刻的内部计数值。需要考虑到内部计数器的溢出影响，广播接收时刻的时差T由式（4）计算：T=Tctu-TRT+mod()Nint-Nlast,M n。（4）式（4）中：n为内部计数器时间量化精度；M为内部计数器发生进位的数值。在每个终端时间周期的开始，进行新终端时间的维护，由式（5）计算：TRT=TRT+mod()Nnow-Nlast,M n+T=Tctu+mod()Nnow-Nlast,M-mod()Nint-Nlast,M n。该方案允许终端进行整体校时，以修正与1553B广播的时间延迟。可以对数管时间进行合法性、均匀性判断，增加可靠性措施。3.2 基于1553B及内部秒脉冲的时间同步方案卫星星敏感器采用基于1553B及内部秒脉冲的时间同步方案。姿轨控计算机通过1553B进行校时，与星敏感器通过 RS422 进行连接，采用内部秒脉冲（Pulse Per Second，PPS）信号进行时间同步，采用RTS（Request to Send，RTS）信号提出数据获取请求。当姿轨控计算机发送PPS信号时，星敏感器锁存对应时刻的本地钟TPPS，单位ms。当姿轨控计算机置低RTS信号时，星敏感器将锁存该时刻的本地钟TRTS，单位ms，并发送最新解算姿态数据包，在数据包中给出曝光时（5）250第3期温渊，等：高光谱观测卫星高精时间同步方案设计与应用刻TExp至最近1次PPS之间的曝光时差TEP。由于RTS为0.25 s一次，PPS为1 s一次，两者的周期不一致，有可能存在RTS信号与曝光时刻TExp之间有PPS信号更新的情况。在计算曝光时刻TExp时间时，需要根据TRTS与TPPS的时差()TRP=TRTS-TPPS进行判断。基于内部秒脉冲的星敏感器各类时间示意图，如图8所示。图8 基于内部秒脉冲的星敏感器各类时间示意图Fig.8 Time schematic of star sensor based on internal PPS根据逻辑关系有星敏感器曝光时刻TExp由式（6）计算：TPPS+TEP,TEPTRPTpps-1 000+TEP,TEPTRP。（6）3.3 可见短波红外高光谱相机高可靠时间同步方案可见短波红外高光谱相机空间分辨率为30 m。为了满足图像定位精度和配准精度的需求，设计采用GPS秒脉冲进行校时，当接收到秒脉冲时，载荷微秒计数器清零。为了适应发生秒脉冲异常、广播异常等可能性，设计有由电控箱维护的内部100 s计数和头部电路维护的内部s计数，分别由时间广播和GPS秒脉冲触发清零操作。另设计“本地缓存时钟1”和“本地缓存时钟2”对时间进行缓存，进一步提高可靠性。可见短波红外高光谱相机的时间码的传递流程如图9所示。图9 可见短波红外高光谱相机时间传递流程Fig.9 AHSI timing flow diagram其中，“本地缓存时钟1”“本地缓存时钟2”和“本地时钟”的定义如下。1）本地缓存时钟1：在可见短波相机内部维护的1个载荷实时钟，用于缓存最新获取的星上时间。在与“本地缓存时钟2”比较后，判断是否需要更新本地时钟，若需要，则可将“本地缓存时钟1”赋值给“本地缓存时钟2”。2）本地缓存时钟2：在可见短波相机内部维护的1个载荷实时钟，用于缓存即将需要更新的本地时钟，“本地缓存时钟2”在成像脉冲触发后延迟赋值到本地时钟。3）本地时钟：可见短波相机内部维护的1个载荷实时钟，在成像脉冲触发时将每帧的成像时刻打入辅助数据包。为了适应数管广播及秒脉冲可能存在故障后的应用，对数管广播及秒脉冲的计数值进行逻辑判断，根据正常范围进行处理。星上GPS秒脉冲设计的精度为1 s。当GPS秒脉冲正常时，微秒计数器计数值范围在1 000 0001以内，超出则认为GPS秒脉冲失效；若星上100 s计数广播正确，则本地缓存时钟1和2 中的时间之差在-5 ms,5 ms内，若监测到超出范围，则可认为星上100 s计数广播失效。3.4 全谱段光谱成像仪简化的时间同步方案全谱段光谱成像仪最高空间分辨率为20 m。为了满足图像定位和通道配准的需求，引入了GPS秒脉冲进行高精度校时。全谱段光谱成像仪管理控制盒通过1553B总线接收数管广播时间码，通过RS-422接口接收GPS接收机发送的秒脉冲。管理控制盒在收到时间码和秒脉冲后不做处理，直接向信号处理盒转 251海 军 航 空 大 学 学 报海 军 航 空 大 学 学 报第38卷发时间码和秒脉冲信号。信号处理盒完成图像辅助数据的编码，包括广播时间码、行同步对应内部计时锁存值、秒脉冲对应内部计时锁存值的记录和打包。全谱段光谱成像仪数管广播时间码和GPS秒脉冲的传递流程如图10所示。图10 全谱段光谱成像仪时间传递流程Fig.10 VIMI timing flow diagram基于该硬件方案，原设计入轨后，采用直接时间码重建方法，即使用广播时间的整秒，加上成像时内部计数器与秒脉冲时刻的差值对应的秒内时差，t=Tn+1+mod()NH-Nn+1,M100。（7）式（7）中：Tn+1为上一个整秒计数；NH为行同步计数器数值；Nn+1为上一个秒脉冲对应计数器数值；M为计数器的计数数量。实际上，由于广播更新时刻与秒脉冲下降沿更新时刻完全不相关，二者的秒长存在差异，如果直接采用广播时钟的整秒，不管是直接取整数部分还是通过四舍五入取整数，均有可能产生整秒值的1 s模糊。为了解决该问题，利用数管广播时间更新时的成像行对应的内部计数，以及秒脉冲更新时刻的内部计数之间的关系，严格重新计算数管广播时刻的精确时间，四舍五入计算秒脉冲更新时刻的准确整秒时间。全谱段光谱成像仪秒脉冲校时实现方案中各计数之间的时序关系，如图11所示。图11 全谱段光谱成像仪秒脉冲校时实现方案原理图Fig.11 Schematic diagram of realizationscheme of VIMI pulse per second time correction4 时间精度分析下面分别给出各时间重建算法的地面重建的精度分析结果。4.1 数管高稳晶振在GPS失效后的时间精度分析已知高稳晶振的时钟频率f为5 MHz。经测试，数 管 计 算 机 A 机 的 高 稳 晶 振 的 稳 定 度f为1.210-9，用t表示待求的漂移时间（守时误差）。t时间内，总晶振脉冲个数为ft，当存在频率偏差时，根据脉冲个数一致，有时钟的漂移方程如下：ft=()t+t()f+ff。（8）两边约去f，则上式可展开为：t=t+t+tf+tf。忽略二阶小量ft，漂移时间（守时误差）可简化为：t=t f。（9）当 门 限t取 0.1 ms 时，可 计 算 出t为8.33104s，约等于0.965 d，即当GPS秒脉冲失效时，时钟误差在1 d左右才漂移出1个时间码字0.1 ms。4.2 星敏感器时间精度分析卫星姿态测量误差包括失重后的重力释放误差、在轨热变化误差、测量低频误差、测量噪声等，在入轨后难以直接评价，地面多通过高分辨率载荷的地标匹配等对姿态绝对测量误差进行评价，通过星敏的夹角法等对姿态的测量噪声进行评价。图12给出进行时间码重建后，在轨星敏感器连续3轨的夹角噪声情况。对该星敏的测量噪声求3倍标准差有3噪声为5.2，而1 ms的时间误差就达到3.6，由此推断出星敏感器的时间重建没有出现跳变等错误。卫星地面系统的在图像定位过程中，对星敏事后处理精度分析结果如表2所示。G P S 接收机管理控制盒信号处理盒图像辅助数据数管分系统秒脉冲广播时间码G P S 秒脉冲（光耦隔离）广播时间码（辅助数据）谱段行计数行同步对应内部计数器数值G P S 秒脉冲对应内部计数器数值广播时间码全谱段光谱成像仪分系统n 1n+2 时刻时刻N时刻秒脉冲本地计数器行同步行同步计时值秒脉冲计时数据+1nN+2nNHNHT1nT+2nT+252第3期温渊，等：高光谱观测卫星高精时间同步方案设计与应用图12 2台星敏感器的夹角噪声情况Fig.12 Noise of the included angle of two star sensors表2 在轨星敏感器随机误差情况统计Tab.2 Random error statistics of star sensors in orbit数据编号数据1数据2数据3星敏组合模式ABACBCABACBCABACBC最大值/()4.8436.6994.4944.7515.8965.0916.5896.9076.813最小值/()-4.198-7.506-7.332-4.464-8.126-4.196-4.885-7.510-5.093平均值/()0.0000.0000.0000.0000.0000.0000.0000.0000.000中误差/()1.9633.2012.2651.7972.8232.1652.0952.9752.299从表2的结果可知，卫星星敏感器的夹角中误差优于3.2，星敏时间重建方法可满足使用的要求。经过热变形标定，图像定位精度从约200 m提升至优于45 m。4.3 可见短波红外高光谱相机时间精度分析利用文中给出的可见短波红外高光谱相机高可靠时间同步方法，任意取1天（2018年11月9日）的光照期遥感数据进行处理，得到成像帧间时差曲线如图13所示。图13 可见短波红外高光谱相机成像帧间时差Fig.13 AHSI timing gaps of adjacent frame从图13可看出，可见短波红外高光谱相机的相邻成像帧间时差的波动值小于6 s，没有出现整秒的不确定情况，时间码的重建精度可满足图像定位及配准要求。4.4 全谱段光谱成像仪时间精度分析利用全谱段光谱成像仪时间同步方法，对卫星第3 505圈（2019年1月4日）的全谱段光谱成像仪遥感数据进行分析，从广播时间递推至下一秒脉冲时间的修正量及秒脉冲时刻数管时间的小数部分，如图14所示。a）秒脉冲时刻距离广播更新修正量a）Time modified from broadcast to PPSb）秒脉冲时刻广播时间递推的小数部分b）Fractional part of PPS from broadcast time图14 从广播时间递推的秒脉冲时刻的修正值及秒内余数Fig.14 Broadcast-derived PPSmodified time and fractional part全谱段光谱成像仪时间码重建后，相邻成像帧间的时差如图15所示。可见，不会出现原直接时间码重建方法出现的整秒跳动情况。从图15可看出，全谱段光谱成像仪的相邻成像帧间时差的波动值小于3 s，时间码的重建精度可满足图像定位及配准要求。图15 全谱段光谱成像仪相邻成像帧间时差Fig.15 VIMI timing gaps of adjacent frame0.970.950.930.910.890.87秒脉冲距离广播修正量/s14181121161201241281累计秒脉冲/次6420-2-4-6夹角/()123456星敏曝光/次1040.004 3520.004 3510.004 3500.004 3490.004 3480.004 3470.004 346相邻帧帧间时差/s14 0018 00112 00116 00120 001累计成像/帧累计秒脉冲/次141811211612012412810.990.980.970.960.950.940.930.920.910.90秒脉冲广播递推间尾数/s相邻帧帧间时差/s0.005 8680.005 8670.005 8660.005 8650.005 8640.005 86315 00110 00115 00120 00125 00130 001累计成像/帧 253海 军 航 空 大 学 学 报海 军 航 空 大 学 学 报第38卷图16给出采用本文提出时间码重建方法的2019年2月25日第4 263圈及2019年5月1日第5 209圈前后帧帧间时差分析结果。可以看出：采用本文的时间码重建方法，时间码不存在跳变；而采用最早提出的基于广播时间整秒+秒脉冲秒内计数的直接时间码重建方法，存在时间码跳变。a）2019年2月25日第4 263圈时间码a）Timecode for lap 4 263 on February 25,2019b）2019年5月1日第5 209圈时间码b）Timecode for lap 5 209 on May 1,2019图16 全谱段光谱成像仪不同时间的相邻成像帧间时差Fig.16 VIMI timing gaps of adjacent framein different imaging times图17给出采用直接时间码重建方法的2019年2月25日第4 263圈及2019年5月1日第5 209圈前后帧帧间时差分析结果。a）2019年2月25日第4 263圈时间码a）Timecode for lap 4 263 on February 25,2019b）2019年5月1日第5 209圈时间码b）Timecode for lap 5 209 on May 1,2019图17 全谱段光谱成像仪不同时间的相邻成像帧间时差（直接法）Fig.17 VIMI timing gaps of adjacent frame in differentimaging times by direct timing method可以看出，采用原来的直接秒脉冲加秒内计数法，在第4 263圈的成像中，时间码存在来回跳变，在第5 209圈的成像中，第1帧时间偏差1 s，后续时间整体偏差1 s。5 结论本文设计了适用于GF-5的高精度时间同步系统方案，对于系统时间，综合采用集中校时、GPS总线校时、GPS秒脉冲校时、校匀校时等多种校时方案。对于终端用户，给出适应一般精度需求的基于1553B中断的时间同步方案、基于1553B及内部秒脉冲的时间同步方案、适应高精度应用需求的可见短波红外高光谱相机高可靠时间同步方案、全谱段光谱成像仪简化的时间同步方案以及推荐的标准秒脉冲应用时间同步方案，可用作遥感卫星时统方案设计的参考。经过评估，在GPS秒脉冲失效后，数管计算机的高稳晶振可以使系统时钟误差在1 d左右漂移出0.1 ms。利用在轨遥感下传数据对星敏感器、可见短波红外高光谱相机及全谱段光谱成像仪曝光时刻的时间重建精度进行了评估。分析表明：星敏感器、可见短波红外高光谱相机、全谱段光谱成像仪这3台仪器的时统方案可满足高精度图像定位及配准的要求。参考文献：1童庆禧,张兵,郑兰芬.高光谱遥感:原理、技术与应用M.北京:高等教育出版社,2006:1-5.TONG QINGXI,ZHANG BING,ZHENG YIFEN.Hy-perspectral remote sensingM.Beijing:Higher Educa-5.867 55.867 05.866 55.866 05.865 55.865 05.864 55.864 05.863 5相邻帧前后帧的时间差/s10-3024681012104成像帧数0.40.20.0-0.2-0.4-0.6-0.8-1.0012345成像帧数104相邻帧前后帧的时间差/s1041.510.50-0.5-1直接相邻帧前后帧的时间差/s024681012成像帧数5.867 55.867 05.866 55.866 05.865 55.865 05.864 55.864 05.863 5相邻帧前后帧的时间差/s012345104成像帧数10-3 254第3期温渊，等：高光谱观测卫星高精时间同步方案设计与应用tion Press,2006:1-5.(in Chinese)2童庆禧,张兵,张立福.中国高光谱遥感的前沿进展J.遥感学报,2016,20(5):689-707.TONG QINGXI,ZHANG BING,ZHANG LIFU.Cur-rent progress of hyperspectral remote sensing in ChinaJ.Journal of Remote Sensing,2016,20(5):689-707.(inChinese)3赵少华,刘思含,吴迪,等.“高分五号”卫星生态环境领域应用前景J.航天返回与遥感,2018,39(3):115-120.ZHAO SHAOHUA,LIU SIHAN,WU DI,et al.Potentialapplication analysis of GF-5 satellite in ecology and envi-ronment fieldsJ.Spacecraft Recovery&Remote Sens-ing,2018,39(3):115-120.(in Chinese)4王锦锦,李真,朱玉玲.高光谱影像在海洋环境监测中的应用J.卫星应用,2019(8):36-40.WANG JINJIN,LI ZHEN,ZHU YULING.Applicationof hyperspectral imagery in marine environmental moni-toring J.Satellite Application,2019(8):36-40.(in Chi-nese)5甘甫平,牟兴林.国土资源业务卫星应用及需求J.卫星应用,2017(11):37-40.GAN FUPING,MOU XINGLIN.Application and de-mand of Land and resources business satelliteJ.SatelliteApplication,2017(11):37-40.(in Chinese)6魏红艳,肖晨超,宋庆君,等.资源一号02D卫星在轨测试自然资源典型应用J.航天器工程,2020,29(6):148-154.WEI HONGYAN,XIAO CHENCHAO,SONG QING-JUN,et al.Natural resources typical applications in on-or-bit test of ZY-1-02D satelliteJ.Spacecraft Engineering,2020,29(6):148-154.(in Chinese)7孙雨,刘家军,赵英俊,等.基于GF-5高光谱数据的蚀变矿物填图及地质应用:以甘肃省瓜州县花牛山地区为例J.中国地质,2022,49(2):558-574.SUN YU,LIU JIAJUN,ZHAO YINGJUN,et al.Altera-tion mineral mapping based on the GF-5 hyperspectral da-ta and its geological application:an example of the Hua-niushan area in Guazhou County of Gansu ProvinceJ.Geology in China,2022,49(2):558-574.(in Chinese)8INGMANN P,VEIHELMANN B,LANGEN J,et al.Re-quirements for the GMES Atmosphere Service and ESAsimplementation concept:sentinels-4/-5 and-5pJ.Re-mote Sensing of Environment,2012,120:58-69.9陈良富,尚华哲,范萌,等.高分五号卫星大气参数探测综述J.遥感学报,2021,25(9):1917-1931.CHEN LIANGFU,SHANG HUAZHE,FAN MENG,etal.Mission overview of the GF-5 satellite for atmospher-ic parameter monitoringJ.National Remote SensingBulletin,2021,25(9):1917-1931.(in Chinese)10 张兴赢,孟晓阳,周敏强,等.卫星高光谱大气CO2探测精度验证研究进展J.气候变化研究进展,2018,14(6):602-612.ZHANGXINGYING,MENGXIAOYANG,ZHOUMINQIANG,et al.Review of the validation of atmo-spheric CO2from satellite hyper spectral remote sensingJ.Climate Change Research,2018,14(6):602-612.(inChinese)11 孙允珠,蒋光伟,李云端,等.高光谱观测卫星及应用前景J.上海航天,2017,34(3):1-13.SUN YUNZHU,JIANG GUANGWEI,LI YUNDUAN,et al.Hyper-spectral observation satellite and its applica-tion prospectsJ.Aerospace Shanghai,2017,34(3):1-13.(in Chinese)12 孙允珠,蒋光伟,李云端,等.“高分五号”卫星概况及应用前景展望J.航天返回与遥感,2018,39(3):1-13.SUN YUNZHU,JIANG GUANGWEI,LI YUNDUAN,et al.GF-5 satellite:overview and application prospectsJ.Spacecraft Recovery&Remo