基于相位偏差的实时PPP模糊度固定

采用CNES发布的实时相位偏差数据,实现包含模糊度固定的实时精密单点定位.对全球10个IGS测站10天观测数据进行RTPPP解算,分别统计模糊度首次固定时间和定位精度,结果显示利用实时相位偏差数据能在平均30min内实现模糊度首次固定,模糊度固定时水平位置误差由6cm迅速降低至2cm左右,三维位置误差由10cm迅速降低至5cm左右,同时RTPPP模糊度固定在3h观测内可保持水平3cm、三维5cm左右的定位精度.通过分析得出,基于相位偏差的RTPPP模糊度固定技术具有较高的定位精度和定位稳定性,能够快速实现cm级定位.     

Predicting atmospheric delays for rapid ambiguity resolution in precise point positioning

Integer ambiguity resolution in precise point positioning (PPP) can shorten the initialization and re-initialization time, and ambiguity-fixed PPP solutions are also more reliable and accurate than ambiguity-float PPP solutions. However, signal interruptions are unavoidable in practical applications, particularly while operating in urban areas. Such signal interruptions can cause discontinuity of carrier phase arc, which introduces new integer ambiguities. Usually it will take approximately 15 min of continuous tracking to a reasonable number of satellites to fix new integer ambiguities. In many applications, it is impractical for a PPP user to wait for such a long time for the re-initialization. In this paper, a method for rapid ambiguity fixing in PPP is developed to avoid such a long re-initialization time. Firstly, the atmospheric delays were estimated epoch by epoch from ambiguity-fixed PPP solutions before the data gap or cycle slip occurs. A random walk procedure is then applied to predict the atmospheric delays accurately over a short time span. The predicted atmospheric delays then can be used to correct the observations which suffer from signal interruptions. Finally, the new ambiguities can be fixed with a distinct WL-LX-L3 (here LX denotes either of L1, L2) cascade ambiguity resolution strategy. Comprehensive experiments have demonstrated that the proposed method and strategy can fix zero-difference integer ambiguities successfully with only a single-epoch observation immediately after a short data gap. This technique works even when all satellites are interrupted at the same time. The duration of data gap bridged by this technique could be possibly extended if a more precise atmospheric delay prediction is found or on-the-fly (OTF) technology is applied. Based on the proposed method, real-time PPP with integer ambiguity fixing becomes more feasible in practice.     

Towards PPP-RTK: Ambiguity resolution in real-time precise point positioning

Integer ambiguity resolution at a single station can be achieved by introducing predetermined uncalibrated phase delays (UPDs) into the float ambiguity estimates of precise point positioning (PPP). This integer resolution technique has the potential of leading to a PPP-RTK (real-time kinematic) model where PPP provides rapid convergence to a reliable centimeter-level positioning accuracy based on an RTK reference network. Nonetheless, implementing this model is technically subject to how rapidly we can fix wide-lane ambiguities, stabilize narrow-lane UPD estimates, and achieve the first ambiguity-fixed solution. To investigate these issues, we used 7 days of 1-Hz sampling GPS data at 91 stations across Europe. We find that at least 10 min of observations are required for most receiver types to reliably fix about 90% of wide-lane ambiguities corresponding to high elevations, and over 20 min to fix about 90% of those corresponding to low elevations. Moreover, several tens of minutes are usually required for a regional network before a narrow-lane UPD estimate stabilizes to an accuracy of far better than 0.1 cycles. Finally, for hourly data, ambiguity resolution can significantly improve the accuracy of epoch-wise position estimates from 13.7, 7.1 and 11.4 cm to 0.8, 0.9 and 2.5 cm for the East, North and Up components, respectively, but a few tens of minutes is required to achieve the first ambiguity-fixed solution. Therefore, from the timeliness aspect, our PPP-RTK model currently cannot satisfy the critical requirement of instantaneous precise positioning where ambiguity-fixed solutions have to be achieved within at most a few seconds. However, this model can still be potentially applied to some near-real-time remote sensing applications, such as the GPS meteorology.     

A new stochastic model considering satellite clock interpolation errors in precise point positioning

Precise clock products are typically interpolated based on the sampling interval of the observational data when they are used for in precise point positioning. However, due to the occurrence of white noise in atomic clocks, a residual component of such noise will inevitable reside within the observations when clock errors are interpolated, and such noise will affect the resolution of the positioning results. In this paper, which is based on a twenty-one-week analysis of the atomic clock noise characteristics of numerous satellites, a new stochastic observation model that considers satellite clock interpolation errors is proposed. First, the systematic error of each satellite in the IGR clock product was extracted using a wavelet de-noising method to obtain the empirical characteristics of atomic clock noise within each clock product. Then, based on those empirical characteristics, a stochastic observation model was structured that considered the satellite clock interpolation errors. Subsequently, the IGR and IGS clock products at different time intervals were used for experimental validation. A verification using 179 stations worldwide from the IGS showed that, compared with the conventional model, the convergence times using the stochastic model proposed in this study were respectively shortened by 4.8% and 4.0% when the IGR and IGS 300-s-interval clock products were used and by 19.1% and 19.4% when the 900-s-interval clock products were used. Furthermore, the disturbances during the initial phase of the calculation were also effectively improved.     

Precise orbit determination and baseline consistency assessment for Swarm constellation

Precise orbit determination (POD) and precise baseline determination (PBD) of Swarm satellites with 4 years of data are investigated. Ambiguity resolution (AR) plays a crucial role in achieving the best orbit accuracy. Swarm POD and PBD based on single difference (SD) AR and traditional double difference (DD) AR methods are explored separately. Swarm antenna phase center variation (PCV) corrections are developed to further improve the orbit determination accuracy. The code multipath of C1C, C1W and C2W observations is first evaluated and clear variations in code noise related to different receiver settings are observed. Carrier phase residuals of different time periods and different loop tracking settings of receiver are studied to explore the effect of ionospheric scintillation on POD. The reduction of residuals in the polar and geomagnetic equator regions confirms the positive impact of the updated carrier tracking loops (TLs) on POD performance. The SD AR orbits and orbits with float ambiguity (FA) are compared with the Swarm precise science orbits (PSOs). An average improvement of 27 %, 4 % and 16 % is gained in along-track, cross-track and radial directions by fixing the ambiguity to integer. For Swarm-A/B and Swarm-B/C formations, specific days are selected to perform the DD AR-based POD during which the average distance of the formation satellites is less than 5000 km. Satellite laser ranging (SLR) observations are employed to validate the performance of FA, SD AR and DD AR orbits. The consistency between the SD AR orbits and SLR data is at a level of 10 mm which shows an improvement of 25 % when comparing with the FA results. An SLR residuals reduction of 15 % is also achieved by the DD AR solution for the selected days. Precise relative navigation is also an essential aspect for spacecraft formation flying missions. The closure error method is proposed to evaluate the baseline precision in three dimensions. A baseline precision of 1–3 mm for Swarm-A/C formation and 3–5 mm for Swarm-A/B and Swarm-B/C satellite pairs is verified by both the consistency check and closure error method.     

Comprehensive outage compensation of real-time orbit and clock corrections with broadcast ephemeris for ambiguity-fixed precise point positioning

The main challenge in real-time precise point positioning (PPP) is that the data outages or large time lags in receiving precise orbit and clock corrections greatly degrade the continuity and real-time performance of PPP positioning. To solve this problem, instead of directly predicting orbit and clock corrections in previous researches, this paper presents an alternative approach of generating combined corrections including orbit error, satellite clock and receiver-related error with broadcast ephemeris. Using ambiguities and satellite fractional-cycle biases (FCBs) of previous epoch and the short-term predicted tropospheric delay through linear extrapolation model (LEM), combined corrections at current epoch are retrieved and weighted with multiple reference stations, and further broadcast to user for continuous enhanced positioning during outages of orbit and clock corrections. To validate the proposed method, two reference station network with different inter-station distance from National Geodetic Survey (NGS) network are used for experiments with six different time lags (i.e., 5 s, 10 s, 15 s, 30 s, 45 s and 60 s), and one set of data collected by unmanned aerial vehicle (UAV) is also used. The performance of LEM is investigated, and the troposphere prediction accuracy of low elevation (e.g., 10–20degrees) satellites has been improved by 44.1% to 79.0%. The average accuracy of combined corrections before and after LEM is used is improved by 12.5% to 77.3%. Without LEM, an accuracy of 2–3 cm can be maintained only in case of small time lags, while the accuracies with LEM are all better than 2 cm in case of different time lags. The performance of simulated kinematic PPP at user end is assessed in terms of positioning accuracy and epoch fix rate. In case of different time lags, after LEM is used, the average accuracy in horizontal direction is better than 3 cm, and the accuracy in up direction is better than 5 cm. At the same time, the epoch fix rate has also increased to varying degrees. The results of the UAV data show that in real kinematic environment, the proposed method can still maintain a positioning accuracy of several centimeters in case of 20 s time lag.     

GNSS global real-time augmentation positioning: Real-time precise satellite clock estimation,prototype system construction and performance analysis

Lots of ambiguities in un-differenced (UD) model lead to lower calculation efficiency, which isn’t appropriate for the high-frequency real-time GNSS clock estimation, like 1 Hz. Mixed differenced model fusing UD pseudo-range and epoch-differenced (ED) phase observations has been introduced into real-time clock estimation. In this contribution, we extend the mixed differenced model for realizing multi-GNSS real-time clock high-frequency updating and a rigorous comparison and analysis on same conditions are performed to achieve the best real-time clock estimation performance taking the efficiency, accuracy, consistency and reliability into consideration. Based on the multi-GNSS real-time data streams provided by multi-GNSS Experiment (MGEX) and Wuhan University, GPS + BeiDou + Galileo global real-time augmentation positioning prototype system is designed and constructed, including real-time precise orbit determination, real-time precise clock estimation, real-time Precise Point Positioning (RT-PPP) and real-time Standard Point Positioning (RT-SPP). The statistical analysis of the 6 h-predicted real-time orbits shows that the root mean square (RMS) in radial direction is about 1–5 cm for GPS, Beidou MEO and Galileo satellites and about 10 cm for Beidou GEO and IGSO satellites. Using the mixed differenced estimation model, the prototype system can realize high-efficient real-time satellite absolute clock estimation with no constant clock-bias and can be used for high-frequency augmentation message updating (such as 1 Hz). The real-time augmentation message signal-in-space ranging error (SISRE), a comprehensive accuracy of orbit and clock and effecting the users’ actual positioning performance, is introduced to evaluate and analyze the performance of GPS + BeiDou + Galileo global real-time augmentation positioning system. The statistical analysis of real-time augmentation message SISRE is about 4–7 cm for GPS, whlile 10 cm for Beidou IGSO/MEO, Galileo and about 30 cm for BeiDou GEO satellites. The real-time positioning results prove that the GPS + BeiDou + Galileo RT-PPP comparing to GPS-only can effectively accelerate convergence time by about 60%, improve the positioning accuracy by about 30% and obtain averaged RMS 4 cm in horizontal and 6 cm in vertical; additionally RT-SPP accuracy in the prototype system can realize positioning accuracy with about averaged RMS 1 m in horizontal and 1.5–2 m in vertical, which are improved by 60% and 70% to SPP based on broadcast ephemeris, respectively.     

RTK model and positioning performance analysis using Galileo four-frequency observations

This paper proposes a real-time kinematic (RTK) model that uses one common reference satellite for the Galileo system with four frequency observations. In the proposed model, the double-differenced (DD) pseudorange and carrier phase biases among the different frequencies are estimated as unknown parameters to recover the integer features of the DD ambiguities among the different frequencies for ambiguity resolution and precise positioning. Analysis results show that the E5a, E5b, and E5 frequencies have virtually the same performance in terms of the positioning accuracy, observation residuals, and ratio values of ambiguity resolution. However, the E1 frequency performs worse than the E5a, E5b, and E5 frequencies. The RTK results for the combination of multiple frequencies are much better than those for a single-frequency observation, the coordinates’ standard deviation is improved about 20–30%, and the ambiguity fix time is improved about 10%.     

Performance analysis of PPP ambiguity resolution with UPD products estimated from different scales of reference station networks

Integer ambiguity fixing with uncalibrated phase delay (UPD) products can significantly shorten the initialization time and improve the accuracy of precise point positioning (PPP). Since the tracking arcs of satellites and the behavior of atmospheric biases can be very different for the reference networks with different scales, the qualities of corresponding UPD products may be also various. The purpose of this paper is to comparatively investigate the influence of different scales of reference station networks on UPD estimation and user ambiguity resolution. Three reference station networks with global, wide-area and local scales are used to compute the UPD products and analyze their impact on the PPP-AR. The time-to-first-fix, the unfix rate and the incorrect fix rate of PPP-AR are analyzed. Moreover, in order to further shorten the convergence time for obtaining precise positioning, a modified partial ambiguity resolution (PAR) and corresponding validation strategy are presented. In this PAR method, the ambiguity subset is determined by removing the ambiguity one by one in the order of ascending elevations. Besides, for static positioning mode, a coordinate validation strategy is employed to enhance the reliability of the fixed coordinate. The experiment results show that UPD products computed by smaller station network are more accurate and lead to a better coordinate solution; the PAR method used in this paper can shorten the convergence time and the coordinate validation strategy can improve the availability of high precision positioning.     

Network based real-time precise point positioning

Network based real-time precise point positioning system includes two stages, i.e. real-time estimation of satellite clocks based on a reference network and real-time precise point positioning thereafter. In this paper, a satellite- and epoch-differenced approach, adopted from what is introduced by Han et al. (2001), is presented for the determination of satellite clocks and for the precise point positioning. One important refinement of our approach is the implementation of the robust clock estimation. A prototype software system is developed, and data from the European Reference Frame Permanent Network on September 19, 2009 is used to evaluate the approach. Results show that our approach is 3 times and 90 times faster than the epoch-difference approach and the zero-difference approach, respectively, which demonstrates a significant improvement in the computation efficiency. The RMS of the estimated clocks is at the level of 0.1 ns (3 cm) compared to the IGS final clocks. The clocks estimates are then applied to the precise point positioning in both kinematic and static mode. In static mode, the 2-h estimated coordinates have a mean accuracy of 3.08, 5.79, 6.32 cm in the North, East and Up directions. In kinematic mode, the mean kinematic coordinates accuracy is of 4.63, 5.82, 9.20 cm.     

北斗卫星导航系统实时定轨与钟差处理策略

目前鲜有对北斗卫星导航系统（BeiDouNavigationSatelliteSystem,BDS）实时精密定轨与钟差确定的研究,文章提出了BDS实时轨道与实时钟差处理策略,包括了观测与动力学模型、实时轨道与实时钟差处理流程与评估方法。尤其对于实时钟差,为了提高计算效率,联合使用两个独立并行的线程估计非差绝对钟差和历元间相对钟差。利用多模全球卫星导航系统试验(MGEX)与全球连续检测评估系统（iGMAS）实测数据进行了北斗实时轨道与钟差解算,BDS实时轨道径向平均精度对于GEO卫星优于20cm,对于IGSO与MEO一般优于10cm;钟差精度对于GEO卫星为0.5~4.5ns,对于IGSO／MEO为0.2~2.0ns。基于目前的轨道与钟差结果,实时精密单点定位（PrecisePointPositioning,PPP）结果可以达到分米量级。     

BDS/GPS对流层延迟估计及对精密单点定位的影响

针对对流层湿延迟参数估计及水平梯度处理策略对精密单点定位影响的问题, 以9个多系统实验网实测数据为基础, 设计了6种对流层延迟处理方案并进行实验. 实验结果表明, 观测时长在8 h时段内, BDS/GPS定位模式下, 对流层延迟的水平梯度处理不宜采用随机游走模型; 采用随机游走估算对流层湿延迟参数, 且不估算或使用分段常数处理水平梯度的策略相对占优. 动静态精密单点定位结果显示, BDS/GPS采用分段常数或随机游走进行对流层湿延迟参数估算和分段常数处理水平梯度的策略相对占优, 采用随机游走处理水平梯度的策略稍弱. 总体而言, BDS较GPS估算的对流层延迟及定位精度稍弱; 海洋周边6种方案解算结果差异较内陆站大. 静态精密单点定位模式下, 随着观测时长的增加, 6种方案定位误差趋势趋于一致, 误差大小基本在2 cm内, 观测时长大于12 h时, 随机游走估算对流层湿延迟参数和不考虑水平梯度的处理策略整体较优.     

GLONASS-based precise point positioning and performance analysis

Current precise point positioning (PPP) techniques are mainly based on GPS which has been extensively investigated. With the increase of available GLONASS satellites during its revitalization, GLONASS observations were increasingly integrated into GPS-based PPP. Now that GLONASS has reached its full constellation, there will be a wide interest in PPP systems based on only GLONASS since it provides a PPP implementation independent of GPS. An investigation of GLONASS-based PPP will also help the development of GPS and GLONASS combined PPP techniques for improved precision and reliability. This paper presents an observation model for GLONASS-based PPP in which the GLONASS hardware delay biases are addressed. In view of frequently changed frequency channel number (FCN) for GLONASS satellites, an algorithm has been developed to compute the FCN for GLONASS satellites using code and phase observations, which avoids the need to provide the GLONASS frequency channel information during data processing. The observation residuals from GLONASS-based PPP are analyzed and compared to those from GPS-based PPP. The performance of GLONASS-based PPP is assessed using data from 15 globally distributed stations.