Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

基于NTCM-BC模型的全球卫星导航系统单频电离层延迟修正

选择NTCM-BC模型作为单频电离层延迟修正模型,通过非线性最小二乘拟合的方法,利用提前一天预测的电离层图（COPG文件）,计算得到NTCM-BC模型修正系数;利用Klobuchar模型和IGS发布的GIM数据对NTCM-BC模型进行比较和分析.对太阳活动高、中、低年实测数据的分析结果表明:全球平均水平上,NTCM-BC模型的电离层延迟修正性能明显优于Klobuchar模型,NTCM-BC模型的TEC平均误差和均方根误差比Klobuchar模型分别下降了41%和30%;模型的TEC计算误差与太阳活动剧烈程度成正相关,即太阳活动高年模型误差较大,太阳活动低年误差相对较低.相较于磁静日,磁扰日期间Klobuchar模型和NCTM模型的误差均有一定程度的增加.此外,模型的电离层修正误差同时存在明显的纬度、季节和地方时差异.      

A real-time combined quality control method for GNSS precise positioning in harsh environments

Global Navigation Satellite System (GNSS) precise positioning can be significantly affected by severe multipath effects and outliers in harsh environments, and highly relies on quality control strategies. Previous studies mainly focus on the posterior residuals to check and exclude the outliers in GNSS observations, limited work emphasizes the combined quality control method considering both the prior and posterior knowledge simultaneously. This paper proposed a real-time combined quality control method to process the multipath effects and outliers in harsh environments simultaneously. Specifically, in the prior stage, a modified multipath processing strategy is proposed for both phase and code observations, then a modified detection, identification, and adaptation (DIA) method considering the maximum times of data snooping is studied in the posterior stage. Two dedicated experiments in real harsh environments were carried out to evaluate the performance of the proposed combined quality control method. For the static experiment, the proposed method exhibits smaller positioning errors, the best positioning accuracy, and the highest availability in this study. Specifically, the proposed method exhibits an improved percentage of 55.4 %, 56.3 %, and 59.7 % for positioning accuracy compared to those without the quality control method in the E, N, and U directions, respectively. Besides, the proposed method can further improve the performance of ambiguity resolution with an improved percentage of 32.2 %. For the kinematic experiment, the three-dimensional positioning accuracy of the proposed method is 0.577 m, which exhibits a 40.0 % improvement compared to those without the quality control method. Also, the proposed method exhibits better performance under relatively strong multipath effects. In this sense, the proposed real-time combined quality control method is highly appreciated in terms of positioning availability, accuracy, and ambiguity resolution for GNSS precise positioning, especially in harsh environments.     

On the estimation of regional covariance functions of TEC variations over Canada

Spatial and temporal variations of Total Electron Content (TEC) can affect GNSS high accuracy positioning. Enhanced estimation of ionospheric variations and their de-correlation can benefit differential and point positioning rapid solutions. Global and regional TEC maps can provide the overall state of ionopsheric variations in space and time domains within their accuracy limits. In this paper, these maps are exploited to retrieve ionospheric variations by means of variograms and their associated covariance functions of TEC residuals over Canadian region during geomagnetically quiet and disturbed conditions. A number of theoretical variogram functions are reviewed for modeling covariance of TEC residuals. The variogram modeling of residuals during a strong geomagnetic storm revealed variances of one order of magnitude larger compared to a rather quiet condition. Variogram models are also used in regional and local kriging interpolation experiments and their performances are evaluated. Global maps of TEC RMS by International GNSS Service and two of its analysis centres are also compared over the Canadian region during a two-year period. Realistic representation of regional variances using estimated variograms when compared to global ionospheric RMS maps are also presented.     

Rapid local ionosphere modeling based on Precise Point Positioning over Iran: A case study under 2014 solar maximum

Presently, the ionosphere effect is the main source of the error in the Global Positioning System (GPS) observations. This effect can largely be removed by using the two-frequency measurements, while to obtain the reasonable results in the single-frequency applications, an accurate ionosphere model is required. Since the global ionosphere models do not meet our needs everywhere, the local ionosphere models are developed. In this paper, a rapid local ionosphere model over Iran is presented. For this purpose, the GPS observations obtained from 40 GPS stations of the Iranian Permanent GPS Network (IPGN) and 16 other GPS stations around Iran have been used. The observations have been selected under 2014 solar maximum, from the days 058, 107, 188 and 271 of the year 2014 with different geomagnetic activities. Moreover, ionospheric observables based on the precise point positioning (PPP) have been applied to model the ionosphere. To represent our ionosphere model, the B-spline basis functions have been employed and the variance component estimation (VCE) method has been used to regularize the problem.To show the efficiency our PPP-derived local ionosphere model with respect to the International GNSS Service (IGS) global models, these models are applied on the single point positioning using single-frequency observations and their results are compared with the precise coordinates obtained from the double-differenced solution using dual-frequency observations. The results show that the 95th percentile of horizontal and vertical positioning errors of the single-frequency point positioning are about 3.1 and 13.6?m, respectively, when any ionosphere model are not applied. These values significantly improve when the ionosphere models are applied in the solutions. Applying CODE’s Rapid Global ionosphere map (CORG), improvements of 59% and 81% in horizontal and vertical components are observed. These values for the IGS Global ionosphere map (IGSG) are 70% and 82%, respectively. The best results are obtained from our local ionosphere model, where 84% and 87% improvements in horizontal and vertical components are observed. These results confirm the efficiency of our local ionosphere model over Iran with respect to the global models. As a by-product, the Differential Code Biases (DCBs) of the receivers are also estimated. In this line, we found that the intra-day variations of the receiver DCBs could be significant. Therefore, these variations must be taken into account for the precise ionosphere modeling.     

基于MLAMBDA的BDS/GPS单频RTK性能分析

GNSS RTK技术以其高精度、高效率、实时性的优点，被广泛应用于航空航天等领域.目前双频RTK技术已非常成熟并且应用较广.相比于双频，单频GNSS RTK在数据质量控制、定位误差处理等方面存在难点.因此单频RTK服务精度可能会受到限制，其定位性能有待研究.本文基于扩展卡尔曼滤波模型，通过MLAMBDA模糊度搜索方法和Ratio检验法，结合实测数据，对比分析BDS，GPS，BDS/GPS三种模式下的单频RTK定位性能.实验证明在静态场景下，三种模式的单频RTK定位精度都在厘米级，可满足高精度定位需求；动态场景下三种模式的模糊度固定率都在70%以上，可满足日常定位需求.在静态及动态应用场景下，北斗的模糊度固定率最高，模糊度解算所用时间短，能实现快速RTK定位.      

Effects of physical correlations on long-distance GPS positioning and zenith tropospheric delay estimates

The global positioning system (GPS) has become an essential tool for the high precision navigation and positioning. The quality of GPS positioning results mainly depends on the model’s formulations regarding GPS observations, including both a functional model, which describes the mathematical relationships between the GPS measurements and unknown parameters, and a stochastic model, which reflects the physical properties of the measurements. Over the past two decades, the functional models for GPS measurements have been investigated in considerable detail. However, the stochastic models of GPS observation data are simplified, assuming that all the GPS measurements have the same variance and are statistically independent. Such assumptions are unrealistic. Although a few studies of GPS stochastic models were performed, they are restricted to short baselines and short time session lengths. In this paper, the stochastic modeling for GPS long-baseline and zenith tropospheric delay (ZTD) estimates with a 24-h session is investigated using the residual-based and standard stochastic models. Results show that using the different stochastic modelling methods, the total differences can reach as much as 3–6 mm in the baseline component, especially in the height component, and 10 mm in the ZTD estimation. Any misspecification in the stochastic models will result in unreliable GPS baseline and ZTD estimations. Using the residual-based stochastic model, not only the precision of GPS baseline and ZTD estimation is obviously improved, but also the baseline and ZTD estimations are closer to the reference value.     

Analysis of GNSS clock prediction performance with different interrupt intervals and application to real-time kinematic precise point positioning

Continuous and timely real-time satellite orbit and clock products are mandatory for real-time precise point positioning (RT-PPP). Real-time high-precision satellite orbit and clock products should be predicted within a short time in case of communication delay or connection breakdown in practical applications. For prediction, historical data describing the characteristics of the real-time orbit and clock can be used as the basis for performing the prediction. When historical data are scarce, it is difficult for many existing models to perform precise predictions. In this paper, a linear regression model is used to predict clock products. Seven-day GeoForschungsZentrum (GFZ) final clock products sampled at 30 s are used to analyze the characteristics of GNSS clocks. It is shown that the linear regression model can be used as the prediction model for the satellite clock products. In addition, the accuracy of the clock prediction for different satellites are analyzed using historical data with different periods (such as 2 and 10 epochs). Experimental results show that the accuracy of the clock with the linear regression prediction model using historical data with 10 epochs is 1.0 ns within 900 s. This is higher accuracy than that achieved using historical data of 2 epochs. Finally, the performance analysis for real-time kinematic precise point positioning (PPP) is provided using GFZ final clock prediction results and state space representation (SSR) clock prediction results when communication delay or connection breakdown occur. Experimental results show that the positioning accuracy without prediction is better than that with prediction in general, whether using the final clock product or the SSR clock product. For the final clock product, the positioning accuracy in the north (N), east (E), and up (U) directions is better than 10.0 cm with all visible GNSS satellites with prediction. In comparison, the 3D positioning accuracy of N, E, and U directions with visible GNSS satellites whose prediction accuracy is better than 0.1 ns using historical data of 10 epochs is improved from 15.0 cm to 7.0 cm. For the SSR clock product, the positioning accuracy of N, E, and U directions is better than 12.0 cm with visible GNSS satellites with prediction. In comparison, the 3D positioning accuracy of N, E, and U directions with visible GNSS satellites whose prediction accuracy is better than 0.1 ns using historical data of 10 epochs is improved from 12.0 cm to 9.0 cm.     

The ionosphere,radio navigation,and global navigation satellite systems

This article is a review of Global Navigation Satellite Systems (GNSS) for space scientists who are interested in how GNSS signals and observables can be used to understand ionospheric dynamics and, conversely, how ionospheric dynamics affect the operational capabilities of GNSS receivers. The most common form of GNSS is the Global Positioning System (GPS); we will first review its operating principles and then present a discussion of errors, of which ionospheric propagation is the most significant. Methods and systems for mitigating errors will be introduced, along with a discussion of modernization plans for GPS and for entirely new systems such as Galileo. In the second half of this article the effects of the ionosphere on GPS signals will be examined in more detail, particularly ionospheric propagation, leading to a discussion of the relation of TEC to ranging errors. Next, the subject of scintillations will also be introduced and connected to the presence and scale sizes of irregularities. Scintillations will be examined as spatial and temporal structures. The method of measuring scintillation pattern drift and ionospheric velocity will be discussed. We conclude by examining ionospheric effects on GPS at midlatitudes.     

Testing a new multivariate GNSS carrier phase attitude determination method for remote sensing platforms

GNSS (Global Navigation Satellite Systems)-based attitude determination is an important field of study, since it is a valuable technique for the orientation estimation of remote sensing platforms. To achieve highly accurate angular estimates, the precise GNSS carrier phase observables must be employed. However, in order to take full advantage of the high precision, the unknown integer ambiguities of the carrier phase observables need to be resolved. This contribution presents a GNSS carrier phase-based attitude determination method that determines the integer ambiguities and attitude in an integral manner, thereby fully exploiting the known body geometry of the multi-antennae configuration. It is shown that this integral approach aids the ambiguity resolution process tremendously and strongly improves the capacity of fixing the correct set of integer ambiguities. In this contribution, the challenging scenario of single-epoch, single-frequency attitude determination is addressed. This guarantees a total independence from carrier phase slips and losses of lock, and it also does not require any a priori motion model for the platform. The method presented is a multivariate constrained version of the popular LAMBDA method and it is tested on data collected during an airborne remote sensing campaign.     

基于ORBCOMM卫星机会信号的定位技术

为摆脱对全球导航卫星系统（GNSS）的依赖，克服其有意或无意干扰情况下无法工作等问题，可采用机会信号（SOP）实现定位，低轨卫星机会信号具备信号功率高、覆盖性广及无需增建基础设施等优点。提出了利用轨道通信卫星（ORBCOMM）系统实现天基机会信号定位。通过对ORBCOMM卫星机会信号的通信体制进行深入研究，实现了利用ORBCOMM卫星机会信号获取多普勒测量信息，建立了瞬时多普勒定位及其几何精度因子的数学模型，并采用卫星TLE数据结合轨道预测模型获得的卫星轨道信息实现ORBCOMM卫星机会信号定位。实测结果表明:利用ORBCOMM卫星机会信号可实现精度优于140 m的定位。研究成果对基于天基机会信号定位技术的理论研究及应用具有重要意义。      

Comparative analysis of four different single-frequency PPP models on positioning performance and atmosphere delay retrieval

Single-frequency precise point positioning (SF-PPP) has attracted increasing attention due to its high precision and cost effectiveness. With various strategies to handle the dominant error, i.e., ionosphere delay, the ionosphere-float (IF), ionosphere-free-half (IFH), ionosphere-corrected (IC), and ionosphere-weighted (IW) SF-PPP models are certain to possess different characteristics and performance levels. This study is dedicated to assessing and comparing the four models from model characteristics, positioning performance, and atmosphere delay retrieval. The model comparison shows that IC and IW models are full-rank while IF and IFH models have a rank deficiency of size one that will result in biased estimations, which means the better solvability of IC and IW models. The experiments are carried out based on the 7-day Global Positioning System (GPS) observations collected at 57 global Multi-GNSS Experiment (MGEX) stations and Global Ionosphere Map (GIM) products. The results indicate that the IW model can accelerate SF-PPP convergence and achieve higher positioning accuracy compared to the other three SF-PPP models, especially in kinematic mode. With convergence criteria of 0.25 m in horizontal and 0.5 m in vertical, the east/north/up convergence times of IW model are 0.5/15.0/25.0 min and 0.5/16.0/36.5 min for static and kinematic modes, respectively. The IW model is able to achieve an instantaneous positioning accuracy of 0.28/0.35/0.75 m. In addition, a real kinematic test also demonstrates the best positioning solutions of IW model. Regarding troposphere delay retrieval, the IF, IFH, and IW models obtain a comparable daily accuracy of 3.0 cm on average, while the IC model achieves the worst accuracy of 8.0 cm. For precise ionosphere delay estimation, IW model only needs an average initialization time of 34.3 min, but a longer initialization time of 51.6 min is required for IF model. The daily precision of ionosphere delay estimation for IW model can reach up to 10.8 cm. At the present accuracy of GIM products, it is suggested that the IW model should be adopted for SF-PPP first due to its superior performance in positioning and atmosphere delay retrieval.     

NeQuick2 validation using multi-satellite and ground data

Global Navigation Satellite System’s (GNSS) positioning calculation is prone to ionospheric errors. Single frequency GNSS users receive ionospheric corrections through broadcast ionospheric models. Therefore, the accuracy of ionospheric models must be validated based on various geographic and geomagnetic conditions. In this work, an attempt is made to validate NeQuick2 electron density (Ne) using multiple sources of space-based and ground-based data at the Arabian Peninsula and for low solar activity conditions. These sources include space-based data from Swarm, DMSP and COSMIC-2 satellite constellations and ground-based data from GNSS receiver and the ionosonde. The period of this study is 1 year from October 2019 to September 2020. Analysis shows that the agreement between NeQuick2 and experimental Ne close to the peak density height depends on seasons and time of the day with the largest errors found in Autumn and during the daytime. NeQuick2 generally overestimates Ne during the daytime. During the early morning and evening hours, Ne estimates seem to be fairly accurate with slight underestimation in Winter and Spring. Estimation of slab thickness by NeQuick2 is found to be close to the values calculated using collocated ionosonde and GNSS receiver.     

Improving the estimation of zenith dry tropospheric delays using regional surface meteorological data

Global Navigation Satellite Systems (GNSS) are emerging as possible tools for remote sensing high-resolution atmospheric water vapour that improves weather forecasting through numerical weather prediction models. Nowadays, the GNSS-derived tropospheric zenith total delay (ZTD), comprising zenith dry delay (ZDD) and zenith wet delay (ZWD), is achievable with sub-centimetre accuracy. However, if no representative near-site meteorological information is available, the quality of the ZDD derived from tropospheric models is degraded, leading to inaccurate estimation of the water vapour component ZWD as difference between ZTD and ZDD. On the basis of freely accessible regional surface meteorological data, this paper proposes a height-dependent linear correction model for a priori ZDD. By applying the ordinary least-squares estimation (OLSE), bootstrapping (BOOT), and leave-one-out cross-validation (CROS) methods, the model parameters are estimated and analysed with respect to outlier detection. The model validation is carried out using GNSS stations with near-site meteorological measurements. The results verify the efficiency of the proposed ZDD correction model, showing a significant reduction in the mean bias from several centimetres to about 5 mm. The OLSE method enables a fast computation, while the CROS procedure allows for outlier detection. All the three methods produce consistent results after outlier elimination, which improves the regression quality by about 20% and the model accuracy by up to 30%.     

Towards high accuracy GNSS real-time positioning with smartphones

The possibility to access undifferenced and uncombined Global Navigation Satellite System (GNSS) measurements on smart devices with an Android operating system allows us to manage pseudorange and carrier-phase measurements to increase the accuracy of real-time positioning. The goal is to perform real-time kinematic network positioning with smartphones, evaluating the positioning accuracy regarding an external mass-market device. The positioning of Samsung Galaxy S8+ and Huawei P10 plus smartphones was performed using a dedicated tool developed by the authors, considering a continuous operating reference station (CORS) network with a mean inter-station distance of about 50?km. The same positioning technique was also applied to an external GNSS low-cost single-frequency receiver (u-blox EVK-M8T) to compare performance between the receiver and antenna embedded in the previous smartphones and this low-cost receiver coupled with a mass-market antenna (Garmin GA38). Attention was also focused on the phase ambiguity resolution, that it is still a challenging aspect for mass-market devices: even if the two smartphones provide slightly different results, the accuracy obtainable today is greater than 60?cm with a precision of few centimetres in real-time, if a CORS network is available. For real-time applications using portable devices, decimetre-level accuracy is sufficient for many applications, such as rapid mapping and search and rescue activities: these results will open new frontiers in terms of real-time positioning with portable low-cost devices.     

基于多约束因子图优化的无人车定位与建图方法

针对目前在特定场景下应用的低速无人车定位系统极度依赖全球导航卫星系统（GNSS）,存在定位精度不高、漂移误差大、受环境影响严重等问题,提出一种低成本、高精度的无人车定位与建图方法。该方法基于三维激光定位与建图（SLAM）技术。首先,使用点云主成分分析（PCA）实现基于特征匹配的激光里程计;其次,将GNSS位置信息、点云分割聚类得到的地平面和点云聚类特征作为位姿约束分别加入图优化框架,消除激光里程计的累积误差;最后,得到最优位姿和大规模场景的点云地图,以实现无人车的自主定位导航。利用包含大型户外城市街道环境的KITTI数据集对所提出的SLAM算法进行了评估,结果表明:系统在3km运动距离情况下定位偏差可控制在1.5 m以下,在局部精度和全局一致性方面均优于其他里程计系统,为无人车的定位提供了新思路。      

LEO enhanced Global Navigation Satellite System (LeGNSS) for real-time precise positioning services

Global Navigation Satellite System (GNSS) has been widely used in many geosciences areas with its Positioning, Navigation and Timing (PNT) service. However, GNSS still has its own bottleneck, such as the long initialization period of Precise Point Positioning (PPP) without dense reference network. Recently, the concept of PNTRC (Positioning, Navigation, Timing, Remote sensing and Communication) has been put forward, where Low Earth Orbit (LEO) satellite constellations are recruited to fulfill diverse missions. In navigation aspect, a number of selected LEO satellites can be equipped with a transmitter to transmit similar navigation signals to ground users, so that they can serve as GNSS satellites but with much faster geometric change to enhance GNSS capability, which is named as LEO constellation enhanced GNSS (LeGNSS). As a result, the initialization time of PPP is expected to be shortened to the level of a few minutes or even seconds depending on the number of the LEO satellites involved. In this article, we simulate all the relevant data from June 8th to 14th, 2014 and investigate the feasibility of LeGNSS with the concentration on the key issues in the whole data processing for providing real-time PPP service based on a system configuration with fourteen satellites of BeiDou Navigation Satellite System (BDS), twenty-four satellites of the Global Positioning System (GPS), and sixty-six satellites of the Iridium satellite constellations. At the server-end, Precise Orbit Determination (POD) and Precise Clock Estimation (PCE) with various operational modes are investigated using simulated observations. It is found out that GNSS POD with partial LEO satellites is the most practical mode of LeGNSS operation. At the user-end, the Geometry Dilution Of Precision (GDOP) and Signal-In-Space Ranging Error (SISRE) are calculated and assessed for different positioning schemes in order to demonstrate the performance of LeGNSS. Centimeter level SISRE can be achieved for LeGNSS.     

An iterative LiDAR DEM-aided algorithm for GNSS positioning in obstructed/rapidly undulating environments

This paper describes a new algorithm to aid stand-alone GNSS positioning in areas of bad signal reception using a Digital Elevation Model (DEM). Traditional Height-Aiding (HA) algorithms assume either a preset (fixed) value for the receiver elevation or rely on the elevation value that corresponds to the nearest available position fix. This may lead in erroneous receiver elevation estimates that, under circumstances, are inefficient to aid effectively GNSS positioning. In this study, the receiver elevation is updated at every iteration step of the navigation solution through dynamic interpolation of the elevation model. The algorithm, because of its ability to extract and fully exploit the elevation information derived from a digital model, it can prove particularly useful in forested areas with steep-sloped terrain. Extended test runs were undertaken to validate the correctness of the mathematical model and the feasibility of the algorithm and associated software. Particularly, analysis of a dataset acquired in a forested, rapidly undulating environment reveals significant average improvement in all performance metrics of positioning, namely the GNSS position availability (50%), accuracy (56%) and external reliability (86%) compared to the Standard Point Positioning (SPP) solution. Moreover, it was found that the method can cope successfully in marginal operating conditions with situations of bad satellite geometry and satellite signals affected by interference due to tree canopy.