构建与学术的桥梁 拉近与权威的距离
复制链接,关注《测绘学报》抖音!
【测绘学报的个人主页】长按复制此条消息,长按复制打开抖音查看TA的更多作品##7NsBSynuc88##[抖音口令]
本文内容来源于《测绘学报》2020年第9期,审图号GS(2020)4848号。
Galileo三频非组合精密定轨模型及精度评估
李星星, 黄健德,袁勇强,李婕,刘城伯,朱艺婷
武汉大学测绘学院, 湖北 武汉 430079
基金项目:国家自然科学基金(41774030;41974027);湖北省自然科学基金(2018CFA081);武汉市科技局应用基础前沿专项(2019010701011395);中德科学中心-中德国际合作项目(M-0054)
摘要:不断丰富的多频信号为GNSS精密数据处理带来了新的机遇与挑战。本文首先推导了适用于多频非组合(UC)观测值的GNSS卫星精密定轨模型,并给出了多频UC模糊度的双差约束策略。在此基础上,本文基于全球分布的150个MGEX测站的观测数据,对UC模型和无电离层组合(IF)模型分别使用E1/E5a、E1/E5b和E1/E5a/E5b观测值进行了Galileo卫星精密定轨。采用与外部精密产品对比、轨道边界不连续性比较和卫星激光测距(SLR)检核等方法评估了不同策略的定轨精度。结果表明:双频情况下,本文提出的UC模型与目前常用的IF模型定轨精度基本一致,两者定轨结果的1DRMS差异在1 mm以内,钟差STD差异在0.01 ns以内,SLR残差差异在2 mm以内。使用E1/E5a/E5b观测值后,UC模型和IF模型的浮点解精度相较于使用E1/E5b观测值的结果有1~2 mm的改善。
关键词:Galileo 非组合模型 无电离层组合 精密定轨

引文格式:李星星, 黄健德, 袁勇强, 等. Galileo三频非组合精密定轨模型及精度评估. 测绘学报,2020,49(9):1120-1130. DOI: 10.11947/j.AGCS.2020.20200320.
阅读全文:http://xb.sinomaps.com/article/2020/1001-1595/2020-9-1120.htm
全文概述
Galileo是由欧盟主持和建设的全球卫星导航系统,旨在提供高精度的导航、定位和授时服务[1]。Galileo于2002年3月开始建设,在2016年12月宣布具备初始服务能力,并计划于2020年全面建成[1]。目前Galileo共有包括4颗在轨试验(in-orbit validation,IOV)卫星和22颗全工作能力(full operational capability,FOC)卫星在内的26颗卫星[2]。稳定可靠的精密轨道产品是Galileo高精度应用的前提。现今,包括欧洲定轨中心(Centre for Orbit Determination in Europe, CODE)、欧空局(European Space Agency, ESA)、德国地学中心(German Research Centre for Geosciences, GFZ)和武汉大学(Wuhan University, WHU)在内的多家MGEX分析中心为Galileo提供精密轨道和钟差产品。上述机构均采用E1/E5a双频无电离层(ionosphere free, IF)组合观测值进行卫星精密定轨[4-6]。
除E1和E5a信号外,Galileo还提供了E5b、E5ab和E6等信号[7]。丰富的卫星频率信号为精密数据处理提供了新的机遇和挑战。非组合(uncombined, UC)模型具有表达形式简单统一且易扩展,观测信息利用充分等优势[8],逐渐成为多频GNSS精密数据处理领域新的研究热门[9]。在精密定位领域,文献[10]提出过一种非组合GPS精密定位方法。文献[11—12]等对非组合精密定位过程中的函数和参数模型进行了深入研究。文献[13]总结了多频GNSS非差非组合精密单点定位模型。文献[14]提出了北斗系统三频非组合和无电离层观测值组合的钟差基准约束准则。文献[15]通过非差非组合模型发展了电离层模型参数化方法。在精密定轨领域,UC模型的相关研究相较于精密定位领域更少。文献[16]基于非组合观测值进行了低轨卫星精密定轨。文献[17—18]曾对非组合精密定轨策略进行了初步的研究。文献[19]在他的研究中论证了双频UC模型和IF模型在GPS卫星精密定轨中的等价性。然而,很少有研究学者使用双频甚至三频信号通过UC模型来进行Galileo卫星精密定轨。本文建立了Galileo三频非组合精密定轨模型,并评估了使用双频和三频的UC和IF模型的Galileo卫星精密定轨结果。
1 Galileo多频定轨模型及策略
1.1 多频非组合定轨模型
测站r和卫星s在某一个历元接收到的卫星频率i的伪距观测量和载波相位观测量可以用如下观测方程表示

(1)
式中,Pi,rs和Li,rs分别表示以m为单位的伪距观测量和载波观测量;ρrs表示站星几何距离;Mrs表示对流层投影系数;Tr表示天顶对流层湿延迟;I1,rs表示f1的电离层延迟一阶项;γis表示第i频率与fi相关的电离层延迟放大因子,可以表示为fi2/f12;c为光速;dti,r和dtis分别表示接收机和卫星钟差;di,r和dis分别表示接收机和卫星端的伪距硬件延迟;bi,r和bis分别表示接收机和卫星端的相位硬件延迟;λis为第i频率的波长;Ni,rs表示以周为单位的整周模糊度;εi,rs和ξi,rs则是伪距观测量和载波观测量的信号噪声。
诸如天线相位中心偏差(phase center offset, PCO)和天线相位中心变化(phase center variation, PCV)、测站天顶对流层干延迟、相对论效应、潮汐改正(包括固体潮、海潮和极潮)和天线相位缠绕效应等需要的误差可通过对应模型正确改正[20-23]。精密定轨中,对式(1)线性化处理后可表示为

(2)
式中,αrs=(αx,αy,αz)为站星几何距离的一阶偏导数;t0为初始历元;t为当前历元;ϕ(t0,t)为初始历元到当前历元的状态转移矩阵;X0s为初始历元坐标;Xr为接收机的位置坐标。由于MGEX分析中心产品通常采用IF模型进行估计,在多频数据处理过程中,为了与MGEX分析中心产品保持一致,需要采用IF基准消除方程秩亏。假定以频率i和j的IF组合为基准,此时观测方程式可以表示为[24]

(3)
式中,dXrs=αrs(ϕ(t0,t)·X0s-Xr)+Mrs·Tr为非频率相关参数;αi,j=fi2/(fi2-fj2)和βi,j=-fj2/(fi2-fj2)为IF组合系数;Br,i,js=c(αi,jdi,r+βi,jdj,r)+c(αi,jdis+βi,jdjs)为引入的硬件延迟改正数;dtr、dts分别为IF组合接收机钟差和卫星钟差。由于钟差基准中所引入的偏差项和卫星及接收机相关,并且伪距和相位观测值中都存在该参数,因此在最小二乘平差准则中,其值能被电离层参数吸收,此时电离层参数可表示为I1,rs=I1,rs+βi,j(di,r-dj,r-dis+djs)。对于频率k,将以频率i和j的无电离层组合基准代入频率k的观测方程,可以得到式(4)

(4)
式中

(5)
从式(5)中可以看到,当以频率i和频率j的双频IF组合为基准时,频率i和频率j的硬件延迟偏差不需要额外估计,其值可以被卫星钟差、接收机钟差及电离层延迟参数吸收。
1.2 双差模糊度固定
在精密定轨过程中,非差UC模糊度和IF模糊度因为吸收了卫星硬件延迟和接收机硬件延迟而不具备整周特性。通过恢复模糊度的整周特性可以得到模糊度固定解并显著提高卫星精密定轨精度。考虑到卫星精密定轨需要在一个地面网络下进行,因此可借助双差策略实现模糊度的固定[25]。对于两个卫星和两个测站,通过线性组合可以将4个UC模糊度组成双差模糊度,此时模糊度中包含的卫星端和接收机端偏差项均被消除。对于IF模糊度,常用式(6)固定

(6)
式中,NNL为双差窄巷(narrow lane, NL)模糊度;NIF和NWL分别表示双差IF和宽巷(wide line,WL)整数模糊度;σNL和σIF分别为双差NL和双差IF模糊度误差。双差WL模糊度可以借助双差的MW(melbourne-wübbena)组合观测量表示[26-28]。在进行三频IF模糊度固定时,可以先对f1和f2频率的IF模糊度进行双差约束,再对f1和f3频率的IF模糊度进行双差约束,从而固定所有模糊度。此时模糊度的双差约束表达式为

(7)
UC模糊度与IF模糊度固定策略基本一致,首先借助UC模糊度得到双差WL整数模糊度,然后基于双差WL和双差UC模糊度得到双差NL整数模糊度。双差UC模糊度约束方程如下所示

(8)
式中,N(r1,r2), (L1,L2)(s1,s2)为f1和f2频率的双差UC模糊度;Ω(L1,L2)为f1和f2频率的双差UC模糊度的协方差阵;α12、β12为f1和f2频率的IF组合系数。最后,对于f1和f2频率,其UC模糊度可以表示为

(9)
在进行三频UC模糊度固定时,可以先对f1和f2频率的UC模糊度进行双差约束,再对f1和f3频率的UC模糊度进行双差约束,从而固定所有模糊度。此时模糊度的双差约束表达式为

(10)
3个频率的模糊度可以表示为

(11)
1.3 数据处理策略
本文利用全球分布的150个MGEX测站在2019年年积日080-120共41天的Galileo观测数据来验证Galileo多频UC模型的定轨效果。在选取的这段时间内,Galileo共有24颗可用卫星,所有测站均可接收到Galileo的E1、E5a和E5b的观测数据。具体测站的全球分布情况如图 1所示。

图 1 MGEX测站分布Fig. 1 Distribution of MGEX stations
图选项
在数据处理中,共使用6种不同的策略进行Galileo卫星的精密定轨和钟差估计。具体策略可见表 1。除了使用E1/E5b双频观测值进行的定轨策略,其他策略均选取E1/E5a的双频IF钟差作为钟差基准。定轨过程中,观测值的采样率为300 s,定轨弧长为24 h,观测值采用高度角加权,截止高度角为7°。电离层采用斜延迟电离层模型,且仅估计E1频率的电离层延迟。接收机端的E1和E5a观测值的天线相位中心改正由GPS的L1和L2观测值对应的模型改正代替,E5b观测值的天线相位中心改正保持和E5a观测值一致。卫星姿态模型选取名义模型,光压模型为ECOM2模型[29]。观测模型、动力学模型和参数估计等的详细策略见表 2。
表 1 Galileo精密定轨的数据选择策略Tab. 1 Data selection strategy for Galileo precise orbit determination
频率组合方式钟差基准策略标识E1/E5aIFIF E1/E5aIF-1E1/E5aUCIF E1/E5aUC-1E1/E5bIFIF E1/E5bIF-2E1/E5bUCIF E1/E5bUC-2E1/E5a/E5bIFIF E1/E5aIF-3E1/E5a/E5bUCIF E1/E5aUC-3表选项
表 2 Galileo精密定轨的处理策略Tab. 2 Processing strategy for Galileo precise orbit determination
模型项目处理策略观测模型观测值E1/E5a/E5b定轨采样率300 s定轨弧长24 h截止高度角7°观测值权比高度角加权。高度角E大于30°时为1,否则为2·sin (E)天线相位改正卫星端:按照IGS14天线文件改正接收机端:E1、E5a/E5b分别由GPS L1和L2的改正数代替相位缠绕采用模型改正对流层延迟对流层干延迟使用Saastamoinen模型计算,通过GMF投影函数将天顶延迟转换到观测值路径方向;对流层湿延迟采用分段常数估计,间隔为2 h动力学模型地球重力场EGM2008模型,12×12阶[30]N体引力根据JPL DE405星历计算地球潮汐形变海潮潮汐:FES2004模型[31]固体潮、极潮:IERS conventions 2010[32]相对论效应采用IERS conventions 2010的模型改正太阳光压ECOM2模型参数估计初始轨道状态卫星位置、速度和ECOM2力模型参数测站坐标IGS周解的测站坐标接收机钟差作为白噪声逐历元估计,频间码偏差作为常数每天估计一组卫星钟差作为白噪声逐历元估计相位模糊度每个弧段估计一个常数模糊度参数ERP参数作为常数估计表选项
2 试验结果与分析
2.1 精密产品比较结果
为了评估轨道和钟差的外符合精度,可以将不同策略的定轨和估钟结果与精密产品进行对比。本文选取CODE提供的精密轨道和钟差产品作为外部评估标准。研究表明,CODE提供的Galileo精密轨道产品与ESA、GFZ和WHU提供的精密轨道产品之间比较结果的1DRMS为16.2~28.3 mm,精密钟差产品之间比较结果的STD为0.09~0.24 ns,另外,CODE机构Galileo精密轨道产品的卫星激光测距(satellite laser ranging,SLR)检核残差的STD为30.0~32.0 mm[33]。
轨道方面,将每个定轨弧段解算的轨道与CODE精密轨道产品进行比较得到切向(along-track, A)、法向(cross-track, C)、径向(radial, R)及1D方向的差值序列,并计算得到差值的RMS值。钟差方面,对解算的钟差结果与CODE精密钟差产品进行比较,使用“二次差”的方法,扣除基准偏差,并计算得到每个定轨弧段的STD值[34]。
图 2给出了Galileo卫星使用双频观测值时,不同策略下定轨结果与CODE精密轨道产品比较结果的1DRMS平均值。其中UC-1、IF-1、UC-2和IF-2策略浮点解RMS平均值分别为44.6、44.7、48.2 mm和48.1 mm,UC-1与IF-1策略间轨道差异和UC-2与IF-2策略间轨道差异都小于1 mm;固定解RMS平均值分别为21.1、21.6、21.4 mm和22.0 mm,4者差异同样小于1 mm。上述结果表明,在使用相同观测值时,UC策略和IF策略的平均RMS差异在1 mm以内,两者精度相当。

图 2 双频UC和IF策略定轨结果Fig. 2 Orbit results of dual-frequency UC and IFmodel
图选项
图 3给出了双频情况下UC-1、UC-2、IF-1和IF-2这4种策略WL和NL模糊度的固定率。表 3进一步给出了上述4种策略下,WL和NL模糊度的平均固定率。从上述结果可以发现,UC-1和UC-2策略的WL和NL模糊度固定率均分别略低于IF-1和IF-2策略。其中,UC-1和UC-2策略的WL模糊度固定率分别为91.65%和92.54%,比IF-1和IF-2策略分别低0.17%和0.65%;UC-1和UC-2策略的NL模糊度固定率分别为82.23%和83.09%,也略低于IF-1和IF-2策略的82.93%和83.58%。

图 3 双频UC和IF策略精密定轨模糊度固定率Fig. 3 Ambiguity fixed rate of dual-frequency UC and IFmodel
图选项
表 3 双频UC和IF策略精密定轨模糊度固定率平均值Tab. 3 Average value of ambiguity fixed rate for dual-frequency UC and IF model
频率WL模糊度固定率表选项
图 4给出了双频和三频UC/IF策略定轨结果与CODE精密产品的互差结果。从图 4中可以发现,对于UC-1、UC-2和UC-3这3种策略,UC-2策略的浮点解精度最差,UC-1策略和UC-3策略的浮点解精度相当,3种策略的固定解精度相当。对于IF-1、IF-2和IF-3这3种策略,IF-2策略的浮点解精度最差,IF-1和IF-2策略的浮点解精度一致,3种策略的固定解精度相当。

图 4 双频和三频UC/IF策略定轨结果与CODE精密产品的互差结果Fig. 4 RMS values oforbit differences for dual-frequency and triple-frequency UC/IF model compared with CODE precision products
图选项
表 4总结了使用双频和三频观测值的UC和IF策略定轨结果与CODE精密产品的互差结果的平均值。可以发现,使用相同观测值时,UC和IF策略定轨结果的精度差异很小,两者在切向、法向与径向差异均在1 mm以内。使用不同观测值时,UC-2策略浮点解的1DRMS平均值为45.2 mm,比UC-1和UC-2策略分别大了3.1 mm和2.0 mm;IF-2策略浮点解的1DRMS平均值同样为45.2 mm,比IF-1和IF-2策略分别大了3.2 mm和2.1 mm。此外,UC-1、UC-2、UC-3、IF-1、IF-2和IF-3这6种策略固定解的1DRMS平均值分别为20.3、20.3、20.4、20.6、20.6 mm和20.4 mm,相互之间差异小于1 mm,可以认为精度相当。
表 4 双频和三频UC/IF策略精密定轨结果与CODE精密产品的互差结果Tab. 4 RMS values of orbit differences for dual-frequency and triple-frequency UC/IF model compared with CODE precision products
观测值UC策略表选项
图 5给出了双频和三频UC/IF策略下,Galileo卫星固定解定轨和估钟结果与CODE精密产品的互差结果。对于定轨结果,在使用相同观测值时,每颗卫星UC策略和IF策略的定轨精度大致相当。当使用不同观测值时,每颗卫星UC策略和IF策略定轨固定解结果精度基本一致。对于估钟结果,其中IF-1、IF-3、UC-1和UC-3这4种策略所有卫星的钟差STD平均值为0.090 ns、0.089 ns、0.087 ns和0.088 ns,相互间STD值差异小于0.01 ns,可以认为钟差精度基本一致。

图 5 双频和三频UC/IF策略精密定轨和估钟结果与CODE精密产品的互差结果Fig. 5 RMS values of clock differences for dual-frequency and triple-frequency UC/IF model compared with CODE precision products
图选项
2.2 轨道边界不连续性
除了和MGEX分析中心提供的精密产品比较外,轨道边界不连续性(orbit boundary discontinuities,OBD)同样可用于评估轨道精度。将24小时定轨结果外推1500 s(采样间隔为300 s时为5个历元),比较外推结果与次日定轨结果的重叠轨道弧段可以得到切向、法向、径向及1D方向的差值序列,并计算得到差值的RMS值。
图 6给出了双频和三频UC/IF策略定轨的OBD结果。从比较结果来看,使用相同观测值时,UC策略和IF策略OBD结果的量级基本一致。使用E1/E5b观测值时,UC策略和IF策略定轨浮点解的OBD结果都略差于使用E1/E5a和E1/E5a/E5b观测值时的结果。表 5进一步给出了双频和三频UC/IF策略定轨的OBD结果的RMS平均值。从表 5中可以发现,在观测值类型相同时,UC策略和IF策略的定轨结果在切向、法向和径向上的OBD结果差异均在2.5 mm以内。UC策略的浮点解中,UC-2策略的精度最差,UC-3策略次之;UC策略的固定解中,UC-2策略的精度最好,UC-3策略次之。对于IF策略,上述结论和UC策略类似。浮点解中,IF-2策略的精度最差,IF-3策略次之;固定解中,IF-2策略的精度最好,IF-3次之。

图 6 双频和三频UC/IF策略精密定轨OBD结果Fig. 6 RMS values of OBD for dual-frequency and triple-frequency UC/IF model compared with CODE precision products
图选项
表 5 双频和三频UC/IF策略精密定轨OBD结果RMS平均值Tab. 5 RMS values of OBD for dual-frequency and triple-frequency UC/IF model compared with CODE precision products
定轨结果策略切向法向径向1D浮点解UC-111.7310.205.639.55UC-213.2910.985.7410.49UC-312.0510.355.619.73IF-111.6910.165.659.52IF-213.3610.985.7710.52IF-312.1110.455.639.79固定解UC-18.377.184.196.81UC-27.896.784.116.46UC-38.016.884.066.53IF-18.327.174.256.80IF-28.086.844.256.59IF-38.147.114.146.68表选项
2.3 SLR检核结果
SLR是一种高精度的绝对距离观测值,因其不包含载波相位模糊度,并具有系统误差小、不受钟差和天线相位中心误差影响的优点,可用于对精密卫星轨道进行检核评估。所有Galileo卫星都配备了用于SLR检核的激光反射器阵列,并由国际激光测距服务(international laser ranging service, ILRS)组织的地面观测网络跟踪[35],因此本文还对不同策略定轨结果进行了SLR检核。
图 7给出了E01、E09、E13、E26卫星在UC-2、IF-2、UC-3和IF-3这4种策略下固定解定轨结果的SLR残差时序图。表 6进一步给出了不同策略下这4颗卫星固定解定轨结果SLR检核结果。需要注意的是,在统计过程中剔除了残差值大于1 m或超过3倍中误差的SLR检核激光点。通过比较可以发现,当观测值类型一致时,UC策略和IF策略定轨结果的SLR残差和STD值都基本相当。总体而言,所有策略的SLR残差平均值的绝对值都小于2 cm,不同策略的SLR残差平均值之间的差异在0.3 cm以内。使用三频观测值时,如果观测值组合策略和双频一致,其SLR残差平均值并没有显著小于双频观测值。以E09为例,其UC-1、UC-2和UC-3策略的SLR残差平均值分别为-2.6 cm、-2.5 cm和-2.6 cm,与IF-1、IF-2和IF-3策略的-2.7 cm、-2.5 cm和-2.5 cm相比,相互间差异都不大于0.2 cm。同样可以比较这几种策略SLR残差STD值,其中UC-1、UC-2和UC-3策略的SLR残差STD值分别为3.0 cm、3.2 cm和2.9 cm,与IF-1、IF-2和IF-3策略的3.0 cm、2.9 cm和3.0 cm相比,相互间差异同样都不大于0.3 cm。

图 7 双频和三频UC/IF策略定轨固定解SLR残差Fig. 7 SLR residuals of fixed solutions for dual-frequency and triple-frequency UC/IF model
图选项
表 6 双频和三频UC/IF策略定轨SLR检核结果Tab. 6 SLR results of fixed solutions for dual-frequency and triple-frequency UC/IF model
卫星PRNSolution#NPOffset/cmSTD/cmE01UC-1381-2.32.1UC-2381-2.22.1UC-3381-2.22.0IF-1381-2.32.1IF-3381-2.32.1IF-3381-2.32.1E09UC-1495-2.63.0UC-2495-2.53.2UC-3495-2.62.9IF-1495-2.73.0IF-2495-2.72.9IF-3495-2.53.0E13UC-1488-1.83.4UC-2488-1.53.7UC-3488-1.83.4IF-1488-1.83.4IF-2488-1.73.4IF-3488-1.73.3E26UC-1300-2.52.6UC-2300-2.23.0UC-3300-2.42.6IF-1300-2.52.5IF-2300-2.22.8IF-3300-2.52.5表选项
3 总结
本文推导了多频UC和IF定轨模型,给出了多频UC定轨模型的双差固定策略。在此基础上,利用全球分布的150个MGEX测站进行了不同策略的Galileo卫星精密定轨。使用与CODE产品比较、OBD和SLR检核这3种方法来评估Galileo定轨产品的精度。通过分析比较结果,得出如下结论:
(1) 在使用E1/E5a与E1/E5b双频观测值的情况下,观测值类型一致时,本文提出的UC策略和传统的IF策略的定轨精度基本相当,两者浮点解轨道精度差异和固定解轨道精度差异都不超过1 mm,浮点解钟差差异和固定解钟差差异都不超过0.01 ns。此外,在双频固定解中,UC-1和UC-2策略的宽巷和窄巷模糊度固定率都分别略低于IF-1和IF-2策略对应的模糊度固定率。
(2) 使用E1/E5a/E5b三频观测值时,浮点解定轨结果相较于双频观测值浮点解定轨结果略有提升。无论是UC策略还是IF策略,其固定解定轨结果的1D RMS和使用E1/E5a和E1/E5b观测值时定轨结果的RMS差异在1 mm以内。双频和三频条件下,UC策略和IF策略固定解的钟差差异在0.01 ns以内。
(3) 在观测值类型一致时,本文提出的UC策略和传统的IF策略的浮点解及固定解OBD结果在切向、法向和径向上的差异小于2.5 mm,不同策略间OBD的1D RMS差异不超过2 mm。
(4) 在观测值类型一致时,UC策略和IF策略的固定解定轨结果的SLR残差和STD值都基本相当。总体而言,所有策略的SLR检核残差小于20 mm,不同策略的SLR检核残差之间的差异不超过3 mm。
(5) Galileo多频非组合精密定轨精度目前仍需进一步改善, 包括E5b信号误差模型精化、UC模糊度双差固定等问题仍需要深入研究和分析。与Galileo类似,北斗三代卫星导航系统(BDS-3)共有B1、B3、B1C、B2a和B2b 5种信号,适用于BDS-3的多频非组合精密定轨模型同样需要进一步研究和评估。
作者简介

第一作者简介:李星星(1985-), 男, 教授, 研究方向为GNSS精密数据处理。E-mail:xxli@sgg.whu.edu.cn
《测绘学报(英文版)》(JGGS)专刊征稿:LiDAR数据处理论文推荐 | 施闯, 张卫星. 曹云昌, 等:基于北斗/GNSS的中国-中南半岛地区大气水汽气候特征及同降水的相关分析
论文推荐 | 张勤, 燕兴元,黄观文,解世超,曹钰:北斗卫星天线相位中心改正模型精化及对精密定轨和定位影响分析
资讯 | “2020年青年测绘地理信息科技创新人才奖”评选结果公告

权威 | 专业 | 学术 | 前沿
微信、抖音小视频投稿邮箱 | song_qi_fan@163.com
欢迎加入《测绘学报》作者QQ群: 751717395
进群请备注:姓名+单位+稿件编号
以上就是李星星抖音红人的全部内容,希望能够对大家有所帮助。
神秘海域手机版
动作格斗62.6M
下载
城市飞车手机版
赛车竞速56.7M
下载
数字华容道中文版
休闲益智81.3M
下载
hideonline中文版
休闲益智49.43MB
下载
逃生1安卓版
冒险解谜21.6M
下载
配音秀最新版本
休闲益智105.83MB
下载
巴士模拟器破解2021版
休闲益智768MB
下载
死亡效应2安卓汉化版
飞行射击1.06G
下载
暗影格斗2手机版
动作格斗147M
下载
艾希官方免费中文版
飞行射击364MB
下载