基于原始观测值的单频精密单点定位算法

研究了一种基于GPS原始观测值的单频PPP算法。该算法通过增加电离层延迟先验信息、空间和时间约束的虚拟观测方程,将电离层延迟当作未知参数与其他定位参数一并进行估计来高效修正电离层延迟误差。通过使用全球178个IGS站1d的实测数据对本算法的收敛速度、定位精度和电离层VTEC的精度进行检验与分析。结果表明,该算法的收敛速度和稳定性均得到了改善,其静态单频单天PPP解的精度可达2~3cm、模拟动态单频单天PPP解的精度可达2~3dm,并且单频PPP与双频PPP提取的电离层总电子含量平均偏差小于5个TECU,可作为一种附属定位产品使用。     

广域实时精密定位与时间服务系统

广域实时精密定位与时间服务已成为GNSS应用领域研究热点,目前国内外学者围绕其模型算法已展开大量的研究。本文重点论述广域实时精密定位与时间服务数据的处理方法和服务系统,给出了基于不同基准约束的卫星钟差解算数学模型,提出通过引入外接原子钟测站、标准时间源(UTC/BDT)等不同时间基准,构建卫星拟稳基准、外接原子钟跟踪站拟稳基准及标准时间源等约束下的钟差解算模型,分析了时间基准对精密单点定位和精密单点授时的影响。本文采用实时卫星轨道、钟差、相位偏差、电离层延迟等服务产品及跟踪站实时数据,验证了系统产品可靠性及终端定位与时间服务性能。实测结果表明:GPS轨道径向精度1.8 cm,钟差STD精度约0.05 ns;BDS-3轨道径向精度6.7 cm,钟差STD精度优于0.1 ns;GPS和BDS-2电离层改正精度分别为0.74 TECU与1.03 TECU。基于该产品实现了用户端PPP、PPP-RTK及PPT、PPT-RTK服务,满足了用户实时厘米级定位和优于0.5 ns的单站时间传递服务,当采用GPS+BDS-2 PPP-RTK解算时,平面收敛至5 cm约需要12 min。     

几种电离层模型对单频单点定位的影响分析

针对实时GNSS单频定位中电离层延迟改正问题,本文采用可用于实时GNSS单频定位的几种电离层模型对电离层延迟进行改正并分析其对GNSS单频单点定位性能的影响。其中,对单频SPP的电离层延迟采用模型直接进行改正,采用Klobuchar模型、CODE的预报产品c1pg、原国家测绘地理信息局的实时球谐电离层产品cosong和CODE事后产品codg计算的电离层精度依次提高;采用不同电离层模型作为电离层估计的先验约束进行单频PPP定位。结果表明:采用精度较好的电离层产品作为先验约束可加快单频PPP收敛。     

A method for improving uncalibrated phase delay estimation and ambiguity-fixing in real-time precise point positioning

In order to improve the performance of precise point positioning (PPP), this paper presents a new data processing scheme to shorten the convergence time and the observation time required for a reliable ambiguity-fixing. In the new scheme, L1 and L2 raw observations are used and the slant ionospheric delays are treated as unknown parameters. The empirical spatial and temporal constraints and the ionospheric delays derived from a real-time available ionospheric model are all considered as pseudo-observations into the estimation for strengthening the solution. Furthermore, we develop a real-time computational procedure for generating uncalibrated phase delays (UPDs) on L1 and L2 frequencies. The PPP solution is first carried out on all reference stations based on the proposed scheme, undifferenced float ambiguities on L1 and L2 frequencies can be directly obtained from the new scheme. The L1 and L2 UPDs are then generated and broadcasted to users in real-time. This data product and also the performance of the new PPP scheme are evaluated. Our results indicate that the new processing scheme considering ionospheric characteristics can reduce the convergence time by about 30 % for float kinematic solutions. The observation time for a reliable ambiguity-fixing is shortened by 25 % compared to that of the traditional ambiguity-fixed kinematic solution. When the new method is used for static reference stations, the observation time for ambiguity-fixing is about 10 min in static mode and only 5  min if the coordinates are fixed to well-known values.     

Multi-GNSS precise point positioning (MGPPP) using raw observations

A joint-processing model for multi-GNSS (GPS, GLONASS, BDS and GALILEO) precise point positioning (PPP) is proposed, in which raw code and phase observations are used. In the proposed model, inter-system biases (ISBs) and GLONASS code inter-frequency biases (IFBs) are carefully considered, among which GLONASS code IFBs are modeled as a linear function of frequency numbers. To get the full rank function model, the unknowns are re-parameterized and the estimable slant ionospheric delays and ISBs/IFBs are derived and estimated simultaneously. One month of data in April, 2015 from 32 stations of the International GNSS Service (IGS) Multi-GNSS Experiment (MGEX) tracking network have been used to validate the proposed model. Preliminary results show that RMS values of the positioning errors (with respect to external double-difference solutions) for static/kinematic solutions (four systems) are 6.2 mm/2.1 cm (north), 6.0 mm/2.2 cm (east) and 9.3 mm/4.9 cm (up). One-day stabilities of the estimated ISBs described by STD values are 0.36 and 0.38 ns, for GLONASS and BDS, respectively. Significant ISB jumps are identified between adjacent days for all stations, which are caused by the different satellite clock datums in different days and for different systems. Unlike ISBs, the estimated GLONASS code IFBs are quite stable for all stations, with an average STD of 0.04 ns over a month. Single-difference experiment of short baseline shows that PPP ionospheric delays are more precise than traditional leveling ionospheric delays.     

远海实时GPS潮汐的实时精密单点定位观测

针对传统事后精密单点定位技术的时间延迟问题,该文基于IGS RTS实时数据流产品,开展了实时精密单点定位技术在远海实时GPS验潮中的应用研究.对RTS改正的实时精密卫星轨道和钟差进行了精度验证和分析,给出了RT-PPP的数据处理策略以及实时GPS验潮的基本流程;组织和实施了渤海湾船载GPS验潮试验,以压力式验潮仪数据为参考,对远距离实时GPS潮汐测量结果进行了精度分析.结果表明:①以IGS最终卫星轨道和钟差产品为参考,RTS实时精密卫星轨道在X、y、Z方向的精度(RMS)均优于3 cm,卫星钟差的精度优于0.15 ns;②采用傅里叶低通滤波方法,消除波浪对潮汐观测的影响,进一步提取潮位信息.在忽略船体姿态改正的情况下,实时精密单点定位验潮相对于压力式验潮仪结果的最大偏差优于20 cm,RMS达到7.5 cm.     

Evaluation of three ionospheric delay computation methods for ground-based GNSS receivers

GNSS observables for ionospheric estimation are commonly based on carrier-to-code leveling (CCL) and precise point positioning (PPP) methods. The CCL method is a geometry-free method which uses carrier phase to level pseudorange observation for decreasing multipath error and observation noise. However, the ionospheric observable based on the CCL has been proven to be affected by leveling errors. The leveling errors are caused by pseudorange multipath and intraday variation of receiver DCB. To obtain more accurate ionospheric observable, the PPP method takes advantage of precise satellite-to-ground range for retrieving slant total electron content and is less affected by the leveling errors. Previous studies have only proven that the ionospheric observables extracted by the two methods are affected by the leveling errors. The influence on ionospheric observable by the pseudorange inter-receiver satellite bias (IRSB) of the receiver has not been taken into consideration. Also, the magnitude of the differences between the ionospheric observables extracted by the two methods has also not been given. In this work, three methods, namely, the CCL, the conventional ionospheric-free PPP method which uses the ionospheric-free Hatch–Melbourne–Wubbena (HMW) function, and the University of Calgary (UOFC) PPP method, are selected to analyze and compare the differences of ionospheric observables and the global ionospheric maps, using a large number of measured data from international GNSS service global stations. Experimental results show that the accuracy of ionospheric observables obtained by the three methods is not only related to the leveling error, but also pseudorange IRSB. The IRSB of the receiver exerts a major effect on the ionospheric observables obtained by the CCL method and a minor effect on the ionospheric observables obtained by the HMW and UOFC methods. The accuracies in the latter case are similar and superior to those obtained by the CCL. The differences of the ionospheric observables obtained by the CCL and UOFC methods, or the CCL and HMW methods, are at decimeter level, whereas the difference of the ionospheric observables obtained by the UOFC and HMW methods is at centimeter level. The UOFC method presented the highest single-frequency pseudorange positioning accuracy using estimated global ionospheric products, followed by the HMW and the CCL methods which presented the lowest positioning accuracy.     

单站区域电离层TEC建模及精度分析

为了分析单站区域电离层总电子含量(total electron content,TEC)模型的适用范围和精度,基于2～15阶次球谐函数,分别建立了欧洲区域16个单站区域电离层TEC模型,生成了区域格网TEC,并与欧洲定轨中心(Center for Orbit Determination in Europe,CODE)、...     

A preliminary evaluation of the performance of multiple ionospheric models in low- and mid-latitude regions of China in 2010–2011

Ionospheric delay is a dominant error source in Global Navigation Satellite System (GNSS). Single-frequency GNSS applications require ionospheric correction of signal delay caused by the charged particles in the earth’s ionosphere. The Chinese Beidou system is developing its own ionospheric model for single-frequency users. The number of single-frequency GNSS users and applications is expected to grow fast in the next years in China. Thus, developing an appropriate ionospheric model is crucially important for the Chinese Beidou system and worldwide single-frequency Beidou users. We study the performance of five globally accessible ionospheric models Global Ionospheric Map (GIM), International Reference Ionosphere (IRI), Parameterized Ionospheric Model (PIM), Klobuchar and NeQuick in low- and mid-latitude regions of China under mid-solar activity condition. Generally, all ionospheric models can reproduce the trend of diurnal ionosphere variations. It is found that all the models have better performances in mid-latitude than in low-latitude regions. When all the models are compared to the observed total electron content (TEC) data derived from GIM model, the IRI model (2012 version) has the best agreement with GIM model and the NeQuick has the poorest agreement. The RMS errors of the IRI model using the GIM TEC as reference truth are about 3.0–10.0 TECU in low-latitude regions and 3.0–8.0 TECU in mid-latitude regions, as observed during a period of 1 year with medium level of solar activity. When all the ionospheric models are ingested into single-frequency precise point positioning (PPP) to correct the ionospheric delays in GPS observations, the PIM model performs the best in both low and mid-latitudes in China. In mid-latitude, the daily single-frequency PPP accuracy using PIM model is ~10 cm in horizontal and ~20 cm in up direction. At low-latitude regions, the PPP error using PIM model is 10–20 cm in north, 30–40 cm in east and ~60 cm in up component. The single-frequency PPP solutions indicate that NeQuick model has the lowest accuracy among all the models in both low- and mid-latitude regions of China. This study suggests that the PIM model may be considered for single-frequency GNSS users in China to achieve a good positioning accuracy in both low- and mid-latitude regions.     

Joint estimation of vertical total electron content (VTEC) and satellite differential code biases (SDCBs) using low-cost receivers

Vertical total electron content (VTEC) parameters estimated using global navigation satellite system (GNSS) data are of great interest for ionosphere sensing. Satellite differential code biases (SDCBs) account for one source of error which, if left uncorrected, can deteriorate performance of positioning, timing and other applications. The customary approach to estimate VTEC along with SDCBs from dual-frequency GNSS data, hereinafter referred to as DF approach, consists of two sequential steps. The first step seeks to retrieve ionospheric observables through the carrier-to-code leveling technique. This observable, related to the slant total electron content (STEC) along the satellite–receiver line-of-sight, is biased also by the SDCBs and the receiver differential code biases (RDCBs). By means of thin-layer ionospheric model, in the second step one is able to isolate the VTEC, the SDCBs and the RDCBs from the ionospheric observables. In this work, we present a single-frequency (SF) approach, enabling the joint estimation of VTEC and SDCBs using low-cost receivers; this approach is also based on two steps and it differs from the DF approach only in the first step, where we turn to the precise point positioning technique to retrieve from the single-frequency GNSS data the ionospheric observables, interpreted as the combination of the STEC, the SDCBs and the biased receiver clocks at the pivot epoch. Our numerical analyses clarify how SF approach performs when being applied to GPS L1 data collected by a single receiver under both calm and disturbed ionospheric conditions. The daily time series of zenith VTEC estimates has an accuracy ranging from a few tenths of a TEC unit (TECU) to approximately 2 TECU. For 73–96% of GPS satellites in view, the daily estimates of SDCBs do not deviate, in absolute value, more than 1 ns from their ground truth values published by the Centre for Orbit Determination in Europe.     

GPS/BDS组合实时动态精密单点定位

针对常规模式下。单系统实时精密单点定位精度受接收机环境和可视卫星数量影响严重等问题,研究了GPS/BDS双系统实时精密单点定位,采用非差无电离层组合载波和伪距观测值,详细推论了Kalman滤波参数估计方法的基本原理,并利用其进行参数估计,最后通过IGS站和实测数据进行了实时PPP实验,实验表明:GPS/BDS双系统定位模式较GPS单系统有明显改善,在E、N、U方向收敛后RMS值分别达到0.125 m、0.117 m、0.289 m,较单系统在各方向分别改善了11.9%、18.1%、22.5%。证明了GPS/BDS实时PPP能够达到分米级到厘米级定位精度。     

Using DORIS measurements for modeling the vertical total electron content of the Earth’s ionosphere

The Doppler orbitography and radiopositioning integrated by satellite (DORIS) system was originally developed for precise orbit determination of low Earth orbiting (LEO) satellites. Beyond that, it is highly qualified for modeling the distribution of electrons within the Earth’s ionosphere. It measures with two frequencies in L-band with a relative frequency ratio close to 5. Since the terrestrial ground beacons are distributed quite homogeneously and several LEOs are equipped with modern receivers, a good applicability for global vertical total electron content (VTEC) modeling can be expected. This paper investigates the capability of DORIS dual-frequency phase observations for deriving VTEC and the contribution of these data to global VTEC modeling. The DORIS preprocessing is performed similar to commonly used global navigation satellite systems (GNSS) preprocessing. However, the absolute DORIS VTEC level is taken from global ionospheric maps (GIM) provided by the International GNSS Service (IGS) as the DORIS data contain no absolute information. DORIS-derived VTEC values show good consistency with IGS GIMs with a RMS between 2 and 3 total electron content units (TECU) depending on solar activity which can be reduced to less than 2 TECU when using only observations with elevation angles higher than \(50^\circ \) . The combination of DORIS VTEC with data from other space-geodetic measurement techniques improves the accuracy of global VTEC models significantly. If DORIS VTEC data is used to update IGS GIMs, an improvement of up to 12  % can be achieved. The accuracy directly beneath the DORIS satellites’ ground-tracks ranges between 1.5 and 3.5 TECU assuming a precision of 2.5 TECU for altimeter-derived VTEC values which have been used for validation purposes.     

估计接收机差分码偏差的GPS/BDS非组合精密单点定位模型

在传统多系统非差非组合精密单点定位（precise point positioning，PPP）模型中，电离层延迟会吸收部分接收机码硬件延迟，其估计值可能为负数。提出了一种估计接收机差分码偏差（differential code bias，DCB）参数的GPS（Global Positioning System）/BDS（BeiDou Navigation Satellite System）非组合PPP模型，将每个系统第1个频率上的接收机码硬件延迟约束为零，对接收机DCB进行参数估计，达到了分离电离层延迟和接收机码硬件延迟的目的，降低了接收机钟差和电离层延迟的相关程度。利用4个多星座实验（multi-GNSS experiment，MGEX）跟踪站的GPS/BDS数据进行了静态和动态PPP试验，结果表明，与不估计DCB参数的PPP模型相比，采用估计DCB参数PPP模型后，静态模式下定位精度和收敛速度平均提高了29.3%和29.8%，动态模式下定位精度和收敛速度平均提高了15.7%和21.6%。     

多系统非差非组合精密单点定位电离层延迟约束权阵的确定

非差非组合精密单点定位需要估计电离层延迟参数,采用电离层先验改正模型约束可以辅助电离层参数解算。针对先验电离层改正量与实际观测量之间权比关系难以确定的问题,本文提出一种电离层约束权因子搜索算法,采用权因子对先验电离层改正量的方差进行调整,根据验后残差加权平方和最小原则通过搜索找出较优的权因子,利用验后残差动态调整先验电离层改正量的方差从而达到改善定位结果的目的。采用8个MGEX跟踪站的GPS/BDS观测数据对该算法进行验证。静态结果表明：对比传统约束方法,采用搜索算法后平均三维定位精度由3.96 cm提高到3.40 cm,平均收敛时间由76.3 min缩短为59.9 min。     

Rapid re-convergences to ambiguity-fixed solutions in precise point positioning

Integer ambiguity resolution at a single receiver can be achieved if the fractional-cycle biases are separated from the ambiguity estimates in precise point positioning (PPP). Despite the improved positioning accuracy by such integer resolution, the convergence to an ambiguity-fixed solution normally requires a few tens of minutes. Even worse, these convergences can repeatedly occur on the occasion of loss of tracking locks for many satellites if an open sky-view is not constantly available, consequently totally destroying the practicability of real-time PPP. In this study, in case of such re-convergences, we develop a method in which ionospheric delays are precisely predicted to significantly accelerate the integer ambiguity resolution. The effectiveness of this method consists in two aspects: first, wide-lane ambiguities can be rapidly resolved using the ionosphere-corrected wide-lane measurements, instead of the noisy Melbourne–Wübbena combination measurements; second, narrow-lane ambiguity resolution can be accelerated under the tight constraints derived from the ionosphere-corrected unambiguous wide-lane measurements. In the test at 90 static stations suffering from simulated total loss of tracking locks, 93.3 and 95.0% of re-convergences to wide-lane and narrow-lane ambiguity resolutions can be achieved within five epochs of 1-Hz measurements, respectively, even though the time latency for the predicted ionospheric delays is up to 180 s. In the test at a mobile van moving in a GPS-adverse environment where satellite number significantly decreases and cycle slips frequently occur, only when the predicted ionospheric delays are applied can the rate of ambiguity-fixed epochs be dramatically improved from 7.7 to 93.6% of all epochs. Therefore, this method can potentially relieve the unrealistic requirement of a continuous open sky-view by most PPP applications and improve the practicability of real-time PPP.     

iGMAS全球电离层TEC格网产品精度分析

为了分析与评估国际GNSS监测评估系统(iGMAS)全球电离层TEC格网产品精度,该文基于iGMAS及IGS各电离层分析中心发布的全球电离层TEC格网产品,进行了精度比较分析,结果表明:iGMAS与IGS、CODE、JPL、ESOC、UPC等IGS电离层工作组发布的全球电离层TEC格网产品,在全球、不同纬度带和欧洲等不同区域均表现出较高的一致性和强相关性,互差为0~2.0 TECU;JPL分析中心GIM的内符合精度约为2.5 TECU,iGMAS、IGS、CODE、ESOC和UPC等分析中心GIM的内符合精度均小于1.5 TECU;在2~8 TECU的精度范围内,iGMAS全球电离层TEC格网产品的精度总体与IGS、CODE、JPL、ESOC、UPC等IGS电离层工作组的精度相当。     

Ionospheric and receiver DCB-constrained multi-GNSS single-frequency PPP integrated with MEMS inertial measurements

Single-frequency precise point positioning (SF-PPP) is a potential precise positioning technique due to the advantages of the high accuracy in positioning after convergence and the low cost in operation. However, there are still challenges limiting its applications at present, such as the long convergence time, the low reliability, and the poor satellite availability and continuity in kinematic applications. In recent years, the achievements in the dual-frequency PPP have confirmed that its performance can be significantly enhanced by employing the slant ionospheric delay and receiver differential code bias (DCB) constraint model, and the multi-constellation Global Navigation Satellite Systems (GNSS) data. Accordingly, we introduce the slant ionospheric delay and receiver DCB constraint model, and the multi-GNSS data in SF-PPP modular together. In order to further overcome the drawbacks of SF-PPP in terms of reliability, continuity, and accuracy in the signal easily blocking environments, the inertial measurements are also adopted in this paper. Finally, we form a new approach to tightly integrate the multi-GNSS single-frequency observations and inertial measurements together to ameliorate the performance of the ionospheric delay and receiver DCB-constrained SF-PPP. In such model, the inter-system bias between each two GNSS systems, the inter-frequency bias between each two GLONASS frequencies, the hardware errors of the inertial sensors, the slant ionospheric delays of each user-satellite pair, and the receiver DCB are estimated together with other parameters in a unique Kalman filter. To demonstrate its performance, the multi-GNSS and low-cost inertial data from a land-borne experiment are analyzed. The results indicate that visible positioning improvements in terms of accuracy, continuity, and reliability can be achieved in both open-sky and complex conditions while using the proposed model in this study compared to the conventional GPS SF-PPP.     

基于RT-PPP的低轨卫星实时高精度时间同步方法

实现低轨导航增强的关键前提是实现低轨星座的整网时间同步,本文针对低轨导航增强系统,提出了一种基于实时精密单点定位(RT-PPP)的低轨卫星高精度时间同步方法,以解决低轨星座实时高精度时间同步的问题. 本文分析了在处理过程中存在的各类误差,介绍了低轨卫星采用状态空间(SSR)改正信息通过精密单点定位(PPP)实现实时高精度时间同步方法的处理流程,将此方法应用于气象、电离层与气候星座观测系统(COSMIC)卫星实测数据的处理,并将该方法与采用广播星历伪距的方法以及事后精密星历的方法进行了比较分析. 结果表明:采用SSR改正信息PPP的方式对2颗COSMIC卫星进行GPS双频观测值的解算,得到的轨道误差的标准差在分米级,钟差误差标准差分别在2.4 ns和2.3 ns左右,可以达到纳秒级. 通过对不同方法解算的结果进行比较可以看出,采用SSR改正信息PPP的方法明显优于采用广播星历伪距方法的解算精度,且与事后精密星历PPP的方法解算精度相当.      

区域电离层建模的单频伪距定位验证

针对电离层单层模型无法满足单频用户定位精度要求的问题,该文建立了区域电离层斜路径模型和单星多项式函数模型。基于河北省区域CORS站的实验结果表明:该文建立的区域电离层模型的拟合精度比CODE中心的格网电离层模型提高了近80%。并将其应用于北斗、GPS单频伪距单点定位,得出北斗的单点定位高程和平面定位精度分别优于3和2m,GPS的单点定位高程和平面定位精度分别优于1.5和1m;相对于CODE,GPS和北斗在平面及高程方向的定位精度均提高了50%左右。结果证明,采用斜路径电离层模型和多项式函数模型,可以较好地反映区域电离层的精细结构。     

Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS

The calibration errors on experimental slant total electron content (TEC) determined with global positioning system (GPS) observations is revisited. Instead of the analysis of the calibration errors on the carrier phase leveled to code ionospheric observable, we focus on the accuracy analysis of the undifferenced ambiguity-fixed carrier phase ionospheric observable determined from a global distribution of permanent receivers. The results achieved are: (1) using data from an entire month within the last solar cycle maximum, the undifferenced ambiguity-fixed carrier phase ionospheric observable is found to be over one order of magnitude more accurate than the carrier phase leveled to code ionospheric observable and the raw code ionospheric observable. The observation error of the undifferenced ambiguity-fixed carrier phase ionospheric observable ranges from 0.05 to 0.11 total electron content unit (TECU) while that of the carrier phase leveled to code and the raw code ionospheric observable is from 0.65 to 1.65 and 3.14 to 7.48 TECU, respectively. (2) The time-varying receiver differential code bias (DCB), which presents clear day boundary discontinuity and intra-day variability pattern, contributes the most part of the observation error. This contribution is assessed by the short-term stability of the between-receiver DCB, which ranges from 0.06 to 0.17 TECU in a single day. (3) The remaining part of the observation errors presents a sidereal time cycle pattern, indicating the effects of the multipath. Further, the magnitude of the remaining part implies that the code multipath effects are much reduced. (4) The intra-day variation of the between-receiver DCB of the collocated stations suggests that estimating DCBs as a daily constant can have a mis-modeling error of at least several tenths of 1 TECU.