• 中国科学引文数据库(CSCD)
  • 中文科技期刊数据库
  • 中国核心期刊(遴选)数据库
  • 日本科学技术振兴机构数据库(JST)
  • 中国学术期刊(网络版)(CNKI)
  • 中国学术期刊综合评价数据库(CAJCED)
  • 中国超星期刊域出版平台

北斗/GNSS机载动态PPP性能分析

赵龙, 郭将, 栗广才, 罗锋

赵龙, 郭将, 栗广才, 罗锋. 北斗/GNSS机载动态PPP性能分析[J]. 全球定位系统, 2022, 47(5): 57-64. DOI: 10.12265/j.gnss.2022096
引用本文: 赵龙, 郭将, 栗广才, 罗锋. 北斗/GNSS机载动态PPP性能分析[J]. 全球定位系统, 2022, 47(5): 57-64. DOI: 10.12265/j.gnss.2022096
ZHAO Long, GUO Jiang, LI Guangcai, LUO Feng. Performance analysis of BDS/GNSS airborne kinematic PPP[J]. GNSS World of China, 2022, 47(5): 57-64. DOI: 10.12265/j.gnss.2022096
Citation: ZHAO Long, GUO Jiang, LI Guangcai, LUO Feng. Performance analysis of BDS/GNSS airborne kinematic PPP[J]. GNSS World of China, 2022, 47(5): 57-64. DOI: 10.12265/j.gnss.2022096

北斗/GNSS机载动态PPP性能分析

基金项目: 中国自然资源航空物探遥感中心青年创新基金(2020YFL02);中央高校基本科研业务费专项资金(2042022kf1035);武汉大学测绘遥感信息工程国家重点实验室专项科研经费资助.
详细信息
    作者简介:

    赵龙: (2000—),男,研究方向为GNSS高精度数据处理

    郭将: (1994—),男,博士,研究方向为高精度定位与数据处理

    栗广才: (1989—),男,博士,研究方向为低成本高精度定位

    罗锋: (1979—),男,工程师,研究方向为航空地球物理勘探

    通信作者:

    郭将 E-mail: guojiang@whu.edu.cn

  • 中图分类号: P228

Performance analysis of BDS/GNSS airborne kinematic PPP

  • 摘要: 针对基于机载动态场景下的实测数据,使用开源的PRIDE PPP-AR软件,采用伪距和载波相位观测值构建双频的消电离层组合(IF)进行了动态精密单点定位(PPP)实验,并对比了单北斗卫星导航系统(BDS)以及BDS/GNSS在进行机载大动态PPP定位方面的性能. 结果表明:多系统组合在卫星数、卫星几何构型以及位置精度衰减因子(PDOP)等方面均优于单系统,且对比单GPS而言,平面(东(E)、北(N)、高程(U))方向的定位精度分别提升了10%和12%;此外将开源软件PRIDE PPP-AR的解算结果与商业软件WayPoint进行了对比,结果表明前者的定位精度在E和U方向分别提升了46%和36%,N方向提升最多,提升了近2倍,因而PRIDE PPP-AR具备更高的动态解算精度与可靠性.
    Abstract: Based on the measured data in the airborne kinematic scene, this paper uses the open-source PRIDE PPP-AR software to construct a dual-frequency ionosphere-free combination (IF) using the pseudorange and carrier phase observations for kinematic precise point positioning (PPP) experiment, and it also compares the performance of a single BDS system and BDS/GNSS in airborne large kinematic positioning. The results show that the number of satellites, satellite geometry and position dilution of precision (PDOP) of the multi-system combination are superior. In a single system, the positioning accuracy in the plane (east (E), north (N) and up (U) )directions are improved by 10% and 12% respectively. In addition, this paper compares the solution results of the open-source software PRIDE PPP-AR with the commercial software WayPoint, the results show that the positioning accuracy of the former is improved by 46% and 36% in the E and U directions respectively, and the N direction is the most improved, nearly doubled. Therefore, PRIDE PPP-AR has higher solution accuracy and better solution performance.
  • 自上世纪美国的GPS建成以来,在相当长一段时间内卫星导航定位领域主要以GPS定位为主,但随着21世纪以来俄罗斯的GLONASS现代化计划的实施,使得GLONASS焕发新生;中国的北斗卫星导航系统(BDS)历经“三步走”战略也已经实现全球组网,面向全球用户提供服务;欧盟的Galileo系统也在持续建设和发展. 当前全球卫星导航系统(GNSS)已由原来的单GPS主导逐渐发展为多系统鼎足而立的阶段,多系统融合定位已经成为当前卫星导航定位技术领域的主要发展趋势之一[1-3]. GNSS多系统融合定位相比单系统定位而言,具有可见卫星数更多、卫星几何构型更好以及可以利用多频观测值组合成波长更长的新观测量,使得模糊度参数能更快被固定[4].

    精密单点定位(PPP)是指利用高精度的载波相位观测值以及卫星星历来实现高精度单点定位的技术,它只需要架设一台接收机,就能实现厘米级甚至毫米级定位精度,自其出现以来就有十分显著的应用优势[5]. 相比相对定位技术,它不受距离限制,作业方便;此外,由于避免差分运算对观测信息造成的损失,PPP技术也能应用于气象监测等方面. 由于上述特点,近年来PPP技术已经成为GNSS领域的热门研究方向,展现出了广阔的应用前景. 如今以国际GNSS服务(IGS)组织推动的多模GNSS实验跟踪网(MGEX)为契机,PPP进一步朝着多系统融合的方向发展[6].

    近年来,许多学者对多系统PPP技术进行了大量研究,文献[7]进行了多系统PPP实验,结果表明多系统在定位稳定性和精度上均优于单系统;文献[8]研究分析了四系统融合PPP的定位性能,实验表明,在单系统几何观测、几何构型不理想的区域,多系统融合能显著提高PPP的定位精度和收敛速度;文献[9]对比分析了 GPS、BDS、GPS+BDS 不同情况下静态与动态PPP精度,结果表明多系统状态下定位精度更高. 文献[10]对多系统的PPP性能进行了分析,结果发现BDS+GPS的定位精度优于其他组合.

    目前PPP技术的研究主要采用MGEX的测网数据进行多系统组合静态或者静态模拟动态实验,对于实际的动态场景尤其是航空动态数据解算中的多系统PPP的定位性能研究较少[11]. 因此,本文选用机载动态场景下的实测的数据,利用IGS提供的精密产品,实现了对单GPS、单BDS和BDS+GPS组合定位系统,BDS+GPS+GLONASS组合系统动态模型下的消电离层组合(IF)的PPP解算,通过对可视卫星数、位置精度衰减因子(PDOP)、验后残差以及平面和高程方向定位精度等方面对定位性能进行了评估和分析.

    PPP广泛使用伪距和载波相位观测值进行解算. 相对论效应、电离层延迟、对流层延迟、天线相位中心偏差(PCO)、相位中心变化(PCV)、潮汐效应均对最终的定位精度产生影响[12],考虑到上述误差改正后的观测方程为:

    $$ \begin{split}{L}_{{\rm{r,}}i}^{{\rm{s}},j}=&{\rho }_{{\rm{r}}}^{{\rm{s}},j}+c·({\rm{d}}{t}_{{\rm{r}}}-{\rm{d}}{t}^{{\rm{s}}})+\lambda ·({N}_{{\rm{r}},i}^{{\rm{s}},j}-{b}^{{\rm{s}},j}+\\& {b}_{{\rm{r,}}i})+\alpha ·{T}_{z}-{\beta }_{f}·{I}_{z}+{\varepsilon }_{\phi } ,\end{split} $$ (1)
    $$ \begin{split}{P}_{{\rm{r,}}i}^{{\rm{s}},j}=&{\rho }_{{\rm{r}}}^{{\rm{s}},j}+c·({\rm{d}}{t}_{{\rm{r}}}-{\rm{d}}{t}^{{\rm{s}}})+c·({{{d}}}_{{\rm{r}},i}^{{\rm{s}}}-{{{d}}}_{{\rm{r}}}^{{\rm{s}},j})+\\& \alpha ·{T}_{z}+{\beta }_{f}·{I}_{z}+{\varepsilon }_{p}.\end{split} $$ (2)

    式中:P、L分别为伪距和载波相位观测值,单位为m;上标j、s分别为卫星系统和卫星,下标r、i分别为接收机和测站;$\;\rho$为卫星与接收机之间的几何距离,单位为m;c为真空中的光速;${\rm{d}}{t_{\rm{r}}}$${\rm{d}}{t^{\rm{s}}}$分别为接收机端钟差和卫星钟差,单位为m;$\lambda $为载波的波长;$N_{{\rm{r}},i}^{{\rm{s}},j}$为整周模糊度;$b_{}^{{\rm{s}},j}$$b_{{\rm{r}},i}^{}$分别表示卫星端和接收机端相位硬件延迟;$d_{}^{{\rm{s}},j}$${{d}}_{{\rm{r}},i}^{}$分别表示卫星端和接收机端的伪距硬件延迟偏差,单位为m;$\alpha $为对流层投影系数;${T_z}$为天顶对流层延迟;${\;\beta _f}$为电离层对应投影系数;${I_z}$为天顶电离层延迟;${\varepsilon _\phi }$${\varepsilon _p}$分别为载波相位与伪距观测噪声,单位为m.

    对于多频的伪距和载波相位观测值而言,为了消除电离层延迟的影响,可以通过构造无电离层组合,此时观测方程为[13]

    $$ P_{{\rm{IF}},i,{\rm{r}}}^{{\rm{s}},j} = \frac{{f_1^2 \cdot {P_1} - f_2^2 \cdot {P_2}}}{{f_1^2 - f_2^2}}, $$ (3)
    $$ L_{{\rm{IF}},i,{\rm{r}}}^{{\rm{s}},j} = \frac{{f_1^2 \cdot {L_1} - f_2^2 \cdot {L_2}}}{{f_1^2 - f_2^2}}.$$ (4)

    式中: ${f_1}$${f_2}$分别为两个载波的频率;${P_1}$${P_2}$分别为两个频率上的伪距观测值;${L_1}$$L{}_2$分别为两个频率上的相位观测值. 式(1)~(2)的观测方变为:

    $$ \begin{split}{P}_{{\rm{IF,r,}}i}^{{\rm{s}},j}=&{\rho }_{{\rm{r}}}^{{\rm{s}},j}+c\cdot \left({\rm{d}}{t}_{{\rm{r}}}-{\rm{d}}{t}^{{\rm{s}}}\right)+\alpha \cdot {T}_{Z}+ \beta \cdot {I}_{Z}+\\&c\cdot \left({{{d}}}_{i,{\rm{r}}}^{{\rm{s}},j}-{{{d}}}_{{\rm{j}}}^{{\rm{s}},j}\right)+{\varepsilon }_{P}\cdot ({P}_{{\rm{IF}}}),\end{split} $$ (5)
    $$ \begin{split} L_{{\rm{IF,r,}}i}^{{\rm{s}},j} =& \rho _{\rm{r}}^{{\rm{s}},j} + c({\rm{d}}{t_{\rm{r}}} - {\rm{d}}{t^{\rm{s}}}) + \lambda \cdot (N_{{\rm{IF}},i}^{{\rm{s}},j} - b_{{i}}^{{\rm{s}},j}+ \\& b_{i,{\rm{r}}}^{{\rm{s}},j}) + \alpha \cdot {T_Z} + {\beta _f} \cdot {I_Z} + {\varepsilon _\phi }({L_{{\rm{IF}}}}),i . \end{split} $$ (6)

    在实际数据解算中,文献[14-15]指出接收机钟差、卫星钟差会分别吸收接收机伪距硬件延迟、卫星伪距硬件延迟;文献[16-17]说明了相位硬件延迟与整周模糊度参数的相关性较强,PPP解算过程中卫星端和接收机端的相位延迟会被模糊度参数吸收,使得模糊度参数失去整数特性,式(5)~(6)就变为:

    $$ P_{{\rm{IF,r}},i}^{{\rm{s}},j} = \rho _{\rm{r}}^{{\rm{s}},j} + c\bar {\rm{d}}{t_{\rm{r}}} + \alpha \cdot {T_Z} + \varepsilon _{\rm{r}}^{\rm{s}}({P_{{\rm{IF}}}}), $$ (7)
    $$ L_{{\rm{IF,r}},i}^{{\rm{s}},j} = \rho _{\rm{r}}^{{\rm{s}},j} + c\bar {\rm{d}}{t_{\rm{r}}} + \alpha \cdot {T_Z} + \varepsilon _{\rm{r}}^{\rm{s}}({P_{{\rm{IF}}}}) + \bar N_{{\rm{r,IF}}}^{\rm{s}}. $$ (8)

    其中:

    $$ c\bar {\rm{d}}{t_{\rm{r}}} = \delta {t_{\rm{r}}} + {{d}}_{i,{\rm{r}}}^{{\rm{s}},j}, $$ (9)
    $$ \bar N_{{\rm{IF,r,}}i}^{{\rm{s}},j} = N_{{\rm{IF,r,}}i}^{{\rm{s}},j} + c(b_{{\rm{IF,r,}}i}^{{\rm{s}},j} + b_{{\rm{IF,}}i}^{{\rm{s}},j}). $$ (10)

    本文进行PPP参数解算时采用的是最小二乘估计,误差方程为

    $$ \mathop {\boldsymbol{V}}\limits_{1 \times 1} = \mathop {\boldsymbol{B}}\limits_{1 \times n} \mathop {\boldsymbol{X}}\limits_{n \times 1} - \mathop {\boldsymbol{L}}\limits_{1 \times 1} . $$ (11)

    式中:$\mathop {\boldsymbol{L}}\limits_{1 \times 1}$为常数矩阵;$\mathop {\boldsymbol{B}}\limits_{1 \times n}$为系数矩阵. 根据平差论中加权残差平方和最小的原理

    $$ {{\boldsymbol{V}}^{\rm{T}}}{\boldsymbol{PV}} = \min . $$ (12)

    $$ {({\boldsymbol{BPX}} - {\boldsymbol{L}})^{\rm{T}}}({\boldsymbol{BPX}} - {\boldsymbol{L}}) = \min. $$ (13)

    法方程可表示为

    $$ {{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{PBX}} = {{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{PL}}. $$ (14)

    则参数解为

    $$ {\boldsymbol{X}} = {({{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{PB}})^{ - 1}}{{\boldsymbol{B}}^{\rm{T}}}{\boldsymbol{PL}}. $$ (15)

    估计的精度为:

    $$ {{\boldsymbol{D}}_{ii}} = {\sigma _0}\sqrt {{{\boldsymbol{Q}}_{ii}}} , $$ (16)
    $$ {\sigma _0} = \sqrt {\frac{{{{\boldsymbol{V}}^{{\rm{T}}}}{\boldsymbol{PV}}}}{r}} . $$ (17)

    式中: ${{\boldsymbol{D}}_{ii}}$为方差矩阵,反映了精度信息;${\sigma _0}$为验前单位权中误差;${{\boldsymbol{Q}}_{ii}}$为协因数矩阵;1为观测数;$r$为多余观测量.

    本实验为在海南省进行的动态机载实验,飞行时间自2021年3月27日6点29分起,结束于当日11点50分(UTC时间),飞行轨迹如图1所示. 最终获得到的实测数据主要包括基准站数据和机载数据.

    图  1  飞机飞行轨迹图

    基准站数据由基准站1与基准站2数据组成. 基准站1的坐标已知,基准站2坐标未知,通过与基准站1联测的方式可以获得基准站2的坐标. 基准站1的数据为2 Hz的单GPS数据,其搭载的接收机类型为Spectra Precision PROFLEX800,数据起止时间为2021年3月27日6点14分至12点24分(UTC时间),基准站2数据情况与基准站1相同.

    机载数据由移动站搭载NovAtel Guiscard接收观测得到,数据为2 Hz的 GPS/BDS/GLONASS的高频实测数据.

    其中各系统包含频率分别为GPS的L1/L2; GLONASS的L1/L2 以及北斗二号(BDS-2)的B1I/B2I. 采集时段从2021年3月27日 6点29分起至11点50分(UTC时间).

    采用表1中的数据处理策略,解算时位置当作随机游走参数来估计,每历元以标准单点定位解作为初值;接收机参数视作噪声,大小给定为$9\;000\;{\text{m}}/\sqrt {\rm{s}} $的白噪声参数进行估计,其中s表示时间单位s;由于飞机速度很大,其所处对流层环境变化剧烈,因此本研究采用了每历元单独估计对流层参数的策略,随机游走噪声给定为 $ 2\;{\text{cm}}/\sqrt {\text{s}} $.

    表  1  数据处理策略
    类别参数策略
    观测值卫星系统与频率GPS:L1/L2/L5
    GLONASS:L1/L2
    BDS:B1I/B2I
    截止高度角
    采样率2 Hz
    DCB产品武汉大学WUM
    快速OSB产品
    ERP产品武汉大学WUM
    快速ERP产品
    轨道与钟差武汉大学WUM
    快速精密钟差轨道产品
    改正模型卫星端/接收机端PCO/PCVIgs14.atx
    对流层模型改正+随机游走
    估计[18]
    电离层无电离层组合
    参数估计位置当作随机游走
    接收机钟差当作白噪声
    模糊度当作常数估计
    下载: 导出CSV 
    | 显示表格

    在数据处理时,为了获得参考真值,首先使用RTKLIB软件对参考站和流动站做了双差解算,固定率为99.3%,这证明以此基线结果作为真值是可行可靠的. 由于此次机载测量实验中北斗三号(BDS-3)缺少双频数据,无法构建双频的无电离层组合,因而本文中所使用的BDS数据均是BDS-2的数据.

    使用PRIDE PPP-AR软件处理动态数据,PRIDE PPP-AR是由武汉大学PRIDE课题组开发的一款可对多系统 GNSS 数据进行后处理,能应用于大地测量学、重力测量、摄影测量、地震监测等方面研究的开源软件. 软件主要经过数据准备和预处理、最小二乘估计以及位置信息输出等步骤,输入精密星历文件和观测值文件最终解算出高精度单点定位坐标,软件处理具体流程如图2所示

    图  2  PRIDE PPP-AR数据处理流程图

    图3可以看出,GPS的卫星分布较均匀,在各个方向均有分布,几何构型也是其中最好的,BDS的卫星多数分布在西南侧,东北侧的卫星分布较少,空间几何构型较差,GLONASS的卫星分布较为均匀,但卫星数较少,因此分布较为稀疏,BDS+GPS+GLONASS的卫星分布更加均匀,高程(U)方向分布的卫星数也更多.

    图  3  GPS、BDS和GLONASS的卫星分布图

    图4可知,GPS+BDS-2+GLONASS的可视卫星数约为20~30颗,飞行过程中约稳定在25颗,相比单GPS、单BDS-2在亚太地区拥有更多得可用卫星,平均10~15颗,GLONASS在观测的开始和结束阶段可视卫星7~9颗,但观测过程大部分时间段内可视卫星数基本为零. 多系统的卫星数显著多于单系统的卫星数,更加适合在动态情况部分卫星被环境遮挡时保持PPP解算的准确估计和稳定持续解算.

    图  4  可视卫星数的时间序列

    PDOP值反映了定位的质量情况,其值越小,表明其卫星处于较好的几何分布状态,定位精度更好. 由图5可知,单BDS-2的PDOP最大,其值在2附近波动,单GPS其PDOP值在1.5波动. GPS+BDS+GLONASS的PDOP值与GPS+BDS相差不大,在1附近波动,主要是因为GLONASS在飞机起飞后其卫星数量急剧减少,推测是由于接收机未能有效跟踪GLONASS卫星,因而实际的定位主要由BDS-2和GPS决定. 相比单系统而言,飞机飞行过程中的大部分时刻多系统的PDOP降低了0.5~1倍,有望获得更高的定位精度.

    图  5  三系统的PDOP值和可用卫星数时间序列

    结合表2的定位精度情况,对于卫星数和几何构型更好的多系统PPP,其PDOP更小,卫星分布更均匀,可见卫星数更多、定位精度更高,表明更多的卫星数以及更好的空间几何构型有助于提高机载动态情况下PPP解算的精度.

    表  2  不同组合下的定位精度 cm
    卫星系统GPSBDSGPS/BDSGPS/BDS/GLONASS
    E 9.1 8.0 7.8 8.1
    N 5.7 5.7 5.4 5.4
    U15.416.713.613.4
    下载: 导出CSV 
    | 显示表格

    观测值残差是由观测值减去改正值后以及参数估计值后的残余部分,包括机载测量过程中来自卫星端、信号传播路径中、接收机端等带来的噪声以及其他未被精确模型化的误差,观测值残差大小可以作为评价一个模型好坏以及观测质量的重要标准.

    定位的精度与观测条件、数据质量直接相关、越精确的观测数据,其定位解算的结果就会相应越好. 残差的噪声量级以及包含的周期信号可以用来分析观测值的数据质量. 虽然PPP中载波相位的权重远大于伪距,但如果伪距的质量太差依然会影响对位置、模糊度参数的准确估计. 图6绘制了GPS、BDS、GLONASS三个系统的载波相位、伪距的残差以及对应的高度角图. 从图6中可以看出,载波相位测量值的残差在±0.1 cm内波动,伪距测量值在±5 m内波动;并且可以看到伪距和载波相位的残差值在飞机起飞和降落阶段波动较大,在飞行过程中保持稳定,主要原因是飞机起飞和降落时环境变化较大,测量误差较大. 从图中也可以看出,BDS-2的卫星高度角大部分在30°以上,且50°以上占比较大,GPS的卫星高度角集中在10°~40°,50°以上占比较小,这与BDS-2在亚太地区的可视卫星数较多有关,BDS-2相对更大的卫星高度角测量会使得多路径造成的误差更小.

    图  6  观测值残差以及卫星高度角

    根据实测数据计算出GPS、BDS、GLONASS的载波相位测量值残差分别是1.45 cm、1.54 cm和1.31 cm,伪距测量值残差分别为1.118 m、1.181 m和1.891 m,三个卫星导航系统的载波相位观测值残差相近,相差10~20 mm,但对于伪距观测值残差,GPS最小,GLONASS最大, GPS与BDS伪距测量残差相近,在1.15 m附近,GLONASS较之降低了约60%. 从测量残差来看,BDS的观测值残差与GPS接近,均低于GLONASS的观测值残差,因此从残差来看,BDS和GPS的观测质量相近,GLONASS的观测质量较差.

    图7为单GPS、单BDS-2、BDS-2+GPS以及BDS-2+GPS+GLONASS的定位精度序列,表2为不同组合系统下的平均定位精度. 定位序列在飞机起飞和降落时有明显的抖动,在U方向尤其明显,推测是对流层模型剧烈变化导致结果估计不稳定. 从表2中可以看出,无论是单GPS、单BDS还是组合系统GPS+BDS以及GPS+BDS+GLONASS,其动态PPP在东(E)方向和北(N)方向的定位精度都要远优于U方向的定位精度,一方面是因为处于平面方向的卫星分布较均匀,卫星空间几何构型较好,而只有U方向有卫星分布,几何构型较差;另一方面,由于电离层延迟、对流层延迟等误差改正项主要是投影到U方向,对于定位的精度影响较大,从表2 中也可以看出,E方向和N方向的平均精度高于U方向60%~90%.

    图  7  不同组合下的定位精度

    同时单BDS的PPP在E方向和N方向较单GPS精度更高,GPS+BDS的PPP在E、N和U方向上均优于单GPS和单BDS,其中在U方向提升更为明显. 我们也可以发现GPS+BDS+GLONASS与GPS+BDS的定位表现相近,主要是由于GLONASS卫星在飞机起飞后可视卫星数目急剧减少,对结果的贡献很小. 这与图4中可用卫星数也可以看出来总卫星数在飞行阶段与GPS+BDS的基本一致. 因此多系统PPP相比单系统PPP而言,能提高定位的精度.

    图8中我们将PRIDE PPP-AR的解算结果同WayPoint差分解做了比较.

    图  8  WayPoint软件基线解

    可以看到在飞机起飞后图中水平方向和U方向表现出明显的“阶跃”,推测是由于模糊度固定错误导致的位置跳变. 最终,定位结果在三个方向定位精度分别达到了13.3 cm、15.4 cm、18.6 cm,WayPoint的解算精度相比PRIDE PPP-AR而言,GPS+BDS的PPP在E方向和U方向分别降低了46%和36%,N方向下降最多,降低了近2倍. 因此PRIDE PPP-AR软件的PPP解算结果显著优于WayPoint软件的解算结果.

    本文通过机载NovAtel GPScard接收机在海南省所进行的动态PPP实验,由于GLONASS在飞机起飞后卫星数急剧减少,对解算结果贡献比较小,因而本文主要是利用PRIDE PPP-AR软件进行了单GPS、单BDS、GPS+BDS以及BDS+GPS+GLONASS的动态PPP解算. 结果表明:

    1)相比GPS,单BDS-2在亚太地区拥有更多的可视卫星,平均10~15颗,因此具有更好的几何构型. 从定位结果来看,单BDS-2动态PPP定位精度在E、N、U三个方向分别为8.0 cm、5.7 cm、16.7 cm,相比单GPS结果在E方向提高了12%. 如果引入BDS-3卫星,定位性能会得到进一步提升.

    2)机载动态情况下多系统GNSS的PPP在实际解算中优势更为明显,其可视卫星数更多,卫星几何构型更好;E、N、U方向的定位精度都高于单系统,U方向提升更加明显,在U方向上,GPS+BDS比单GPS、单BDS精度提升了12%和19%,这对于高程方向更加敏感的应用场景如航空重力测量会更有帮助.

    3)本文利用PRIDE PPP-AR软件的动态PPP解算结果与WayPoint软件差分解对比,结果表明PRIDE PPP-AR解算结果精度相比前者在E和U方向分别提高了46%和36%,N方向提高了近两倍,因此其解算结果显著优于WayPoint软件的解算结果,且前者完全开源,而后者为商业软件,因此PRIDE PPP-AR软件相比WayPoint软件具有更加显著的解算性能与应用能力.

    致谢:感谢武汉大学PRIDE Lab 提供的开源软件PRIDE PPP-AR的支持.感谢中国自然资源航空物探遥感中心提供的机载观测数据.

  • 图  1   飞机飞行轨迹图

    图  2   PRIDE PPP-AR数据处理流程图

    图  3   GPS、BDS和GLONASS的卫星分布图

    图  4   可视卫星数的时间序列

    图  5   三系统的PDOP值和可用卫星数时间序列

    图  6   观测值残差以及卫星高度角

    图  7   不同组合下的定位精度

    图  8   WayPoint软件基线解

    表  1   数据处理策略

    类别参数策略
    观测值卫星系统与频率GPS:L1/L2/L5
    GLONASS:L1/L2
    BDS:B1I/B2I
    截止高度角
    采样率2 Hz
    DCB产品武汉大学WUM
    快速OSB产品
    ERP产品武汉大学WUM
    快速ERP产品
    轨道与钟差武汉大学WUM
    快速精密钟差轨道产品
    改正模型卫星端/接收机端PCO/PCVIgs14.atx
    对流层模型改正+随机游走
    估计[18]
    电离层无电离层组合
    参数估计位置当作随机游走
    接收机钟差当作白噪声
    模糊度当作常数估计
    下载: 导出CSV

    表  2   不同组合下的定位精度 cm

    卫星系统GPSBDSGPS/BDSGPS/BDS/GLONASS
    E 9.1 8.0 7.8 8.1
    N 5.7 5.7 5.4 5.4
    U15.416.713.613.4
    下载: 导出CSV
  • [1] 宁津生, 姚宜斌, 张小红. 全球导航卫星系统发展综述[J]. 导航定位学报, 2013, 1(1): 3-8. DOI: 10.3969/j.issn.2095-4999.2013.01.002
    [2] 季宇虹, 王让会. 全球导航定位系统GNSS的技术与应用[J]. 全球定位系统, 2010, 35(5): 69-75. DOI: 10.3969/j.issn.1008-9268.2010.05.016
    [3] 刘基余. GNSS全球导航卫星系统的新发展[J]. 遥测遥控, 2007, 28(4): 1-6. DOI: 10.3969/j.issn.2095-1000.2007.04.001
    [4] 耿江辉, 常华, 郭将, 等. 面向城市复杂环境的3种多频多系统GNSS单点高精度定位方法及性能分析[J]. 测绘学报, 2020, 49(1): 1-13.
    [5] 张小红, 李星星, 李盼. GNSS精密单点定位技术及应用进展[J]. 测绘学报, 2017, 46(10): 1399-1407. DOI: 10.11947/j.AGCS.2017.20170327
    [6] 屈利忠, 赵齐乐, 郭靖, 等. GNSS融合动态精密单点定位性能分析[J]. 大地测量与地球动力学, 2016, 36(4): 298-302. DOI: 10.14075/j.jgg.2016.04.004
    [7] 苗维凯, 陈旭升, 刘洋洋. GNSS多系统组合PPP解算方法与成果分析[J]. 大地测量与地球动力学, 2019, 39(5): 521-527. DOI: 10.14075/j.jgg.2019.05.016
    [8] 任晓东, 张柯柯, 李星星, 等. BeiDou、Galileo、GLONASS、GPS多系统融合精密单点[J]. 测绘学报, 2015, 44(12): 1307-1313, 1339.
    [9] 祝会忠, 杨添宇, 赵洪涛, 等. GNSS多系统精密单点定位方法与性能分析[J]. 测绘科学, 2020, 45(12): 1-7,21. DOI: 10.16251/j.cnki.1009-2307.2020.12.001
    [10] 尹海博, 郭杭, 罗孝文. BDS/GNSS融合精密单点定位性能分析[J]. 全球定位系统, 2021, 46(3): 66-71. DOI: 10.12265/j.gnss.2020121501
    [11] 高永刚, 吴观烨, 郭金运, 等. GNSS多系统组合动态精密单点定位性能研究[J]. 应用基础与工程科学学报, 2022, 30(2): 501-518. DOI: 10.16058/j.issn.1005-0930.2022.02.020
    [12] 夏宇飞. INS辅助的多系统GNSS单频精密单点定位算法研究[D]. 上海: 华东师范大学, 2018.
    [13] 杨添宇. GNSS精密单点定位算法研究及性能分析[D]. 阜新: 辽宁工程技术大学, 2021.
    [14] 胡丽乐, 潘林. GPS/GLONASS/BDS/Galileo组合精密单点定位[J]. 导航定位学报, 2017, 5(1): 86-90. DOI: 10.16547/j.cnki.10-1096.20170118
    [15] 丁赫, 孙付平, 李亚萍, 等. BDS/GPS/GLONASS组合精密单点定位模型及性能分析[J]. 大地测量与地球动力学, 2016, 36(4): 303-307. DOI: 10.14075/j.jgg.2016.04.005
    [16] 李盼. GNSS 精密单点定位模糊度快速固定技术和方法研究[D]. 武汉: 武汉大学, 2016.
    [17] 屈利忠. 北斗/GNSS实时精密定位服务系统关键算法研究与实现[D]. 武汉: 武汉大学, 2017.
    [18]

    SAASTAMOINEN J. Contributions to the theory of atmospheric refraction[J]. Bulletin géodésique, 1972: 279-298.

  • 期刊类型引用(1)

    1. 易卿武,廖桂生,蔚保国,王彬彬,肖恭伟. 精密单点定位零基准非差模糊度固定方法. 空间电子技术. 2023(03): 39-45 . 百度学术

    其他类型引用(2)

图(8)  /  表(2)
计量
  • 文章访问数:  549
  • HTML全文浏览量:  95
  • PDF下载量:  63
  • 被引次数: 3
出版历程
  • 收稿日期:  2022-05-29
  • 网络出版日期:  2022-09-25
  • 刊出日期:  2022-11-04

目录

/

返回文章
返回
x 关闭 永久关闭