Vibration displacement extraction based on an auto-tuning Kalman smoother from GNSS

The Global Navigation Satellite System (GNSS) is a widely used positioning technique. Computational efficiency is crucial to applications such as real-time GNSS positioning and GNSS network data processing. Many researchers have made great efforts to address this problem by means such as parameter elimination or satellite selection. However, parameter estimation is rarely discussed when analyzing GNSS algorithm efficiency. In addition, most studies on Kalman filter (KF) efficiency commonly have defects, such as neglecting application-specified optimization and limiting specific hardware platforms in the conclusion. The former reduces the practicality of the solution, because applications that need such analyses on filters are often optimized, and the latter reduces its generality because of differences between platforms. In this paper, the computational cost enhancement of replacing the conventional KF with the information filter (IF) is tested considering GNSS application-oriented optimization conditions and hardware platform differences. First, optimization conditions are abstracted from GNSS data-processing scenarios. Then, a thorough analysis is carried out on the computational cost of the filters, considering hardware–platform differences. Finally, a case of GNSS dynamic differencing positioning is studied. The simulation shows that the IF is slightly faster for precise point positioning and much faster for the code-based single-difference GNSS (SDGNSS) with the constant velocity (CV) model than the conventional KF, but is not a good substitute for the conventional KF in the other algorithms mentioned. The real test shows that the IF is about 50% faster than the conventional KF handling code-based SDGNSS with the CV model. Also, the information filter is theoretically equivalent to and can produce results that are consistent with the Kalman filter. Our conclusions can be used as a reference for GNSS applications that need high process speed or real-time capability.

Global Navigation Satellite System (GNSS) plays an increasingly important role in all walks of life. In order to meet the demands of different users, it is crucial to establish corresponding correct mathematical models, especially in the field of precise positioning and navigation. However, due to the spatiotemporal complexity of and limited knowledge on GNSS errors, some residual errors would inevitably remain even after being corrected with differencing and linear combination, empirical model correction, and traditional parameterization. These residual observation errors are referred to as unmodeled errors. In fact, the unmodeled errors have adverse impacts on high-precision GNSS positioning. However, most existing studies mainly focus on handling the part of systematic errors that can be adequately modeled and then simply ignore the part of unmodeled errors that may actually exist. To make a breakthrough in the precision and reliability of GNSS applications currently, this study focuses on the theory and method for processing the GNSS unmodeled errors based on the mathematical model compensation, including resilient functional model adjustment and resilient stochastic model optimization. Specifically, according to the significance and properties of unmodeled errors, the unmodeled-error-ignored model, unmodeled-error-corrected model, unmodeled-error-fixed model, unmodeled-error-float model, and unmodeled-error-weighted model are proposed, where the approaches of significance testing, geometry-free/based/fixed models, hemispherical/hierarchy maps, composite stochastic model, multi-epoch partial parameterization, and inequality and equality constraints are adopted. Ultimately, an unmodeled error processing flow that can adaptively adjust as the external conditions change is proposed, then one can obtain high-precision GNSS positioning solutions.

High accurate global navigation satellite systems (GNSS) require to correct a signal delay caused by the troposphere. The delay can be estimated along with other unknowns or introduced from external models. We assess the impact of the recently developed augmentation tropospheric model on real-time kinematic precise point positioning (PPP). The model is based on numerical weather forecast and thus reflects the actual state of weather conditions. Using the G-Nut/Geb software, we processed GNSS and meteorological data collected during the experiment using a hot-air balloon flying up to an altitude of 2000 m. We studied the impacts of random walk noise setting of zenith total delay (ZTD) on estimated parameters and the mutual correlations, the use of external tropospheric corrections, the use of data from a single or dual GNSS constellation and the use of Kalman filter and backward smoothing processing methods. We observed a significant negative correlation of the estimated rover height and ZTD which depends on constraining ZTD estimates. Such correlation caused a degraded performance of both parameters when estimated simultaneously, in particular for a single GNSS constellation. The impact of ZTD constraining reached up to 50-cm differences in the rover height. Introducing external tropospheric corrections improved the PPP solution regarding: (1) shortened convergence, (2) better overall robustness, particularly, in case of degraded satellite geometry, (3) less adjusted parameters with lower correlations. The numerical weather model-driven PPP resulted in 9–12- and 5–6-cm uncertainties in the rover altitude using the Kalman filter and the backward smoothing, respectively. Compared to standard PPP, it indicates better performance by a factor of 1–2 depending on the availability of GNSS constellations, the troposphere constraining and the processing strategy.

During the past few decades, the presence of global navigation satellite systems (GNSSs) such as GPS, GLONASS, Beidou and Galileo has facilitated positioning, navigation and timing (PNT) for various outdoor applications. With the rapid increase in the number of orbiting satellites per GNSS, enhancements in the satellite-based augmentation systems (SBASs) such as EGNOS and WAAS, as well as commissioning new GNSS constellations, the PNT capabilities are maximized to reach new frontiers. Additionally, the recent developments in precise point positioning (PPP) and real time kinematic (RTK) algorithms have provided more feasibility to carrier-phase precision positioning solutions up to the third-dimensional localization. With the rapid growth of internet of things (IoT) applications, seamless navigation becomes very crucial for numerous PNT dependent applications especially in sensitive fields such as safety and industrial applications. Throughout the years, GNSSs have maintained sufficiently acceptable performance in PNT, in RTK and PPP applications however GNSS experienced major challenges in some complicated signal environments. In many scenarios, GNSS signal suffers deterioration due to multipath fading and attenuation in densely obscured environments that comprise stout obstructions. Recently, there has been a growing demand e.g. in the autonomous-things domain in adopting reliable systems that accurately estimate position, velocity and time (PVT) observables. Such demand in many applications also facilitates the retrieval of information about the six degrees of freedom (6-DOF - x, y, z, roll, pitch, and heading) movements of the target anchors. Numerous modern applications are regarded as beneficiaries of precise PNT solutions such as the unmanned aerial vehicles (UAV), the automatic guided vehicles (AGV) and the intelligent transportation system (ITS). Hence, multi-sensor fusion technology has become very vital in seamless navigation systems owing to its complementary capabilities to GNSSs. Fusion-based positioning in multi-sensor technology comprises the use of multiple sensors measurements for further refinement in addition to the primary GNSS, which results in high precision and less erroneous localization. Inertial navigation systems (INSs) and their inertial measurement units (IMUs) are the most commonly used technologies for augmenting GNSS in multi-sensor integrated systems. In this article, we survey the most recent literature on multi-sensor GNSS technology for seamless navigation. We provide an overall perspective for the advantages, the challenges and the recent developments of the fusion-based GNSS navigation realm as well as analyze the gap between scientific advances and commercial offerings. INS/GNSS and IMU/GNSS systems have proven to be very reliable in GNSS-denied environments where satellite signal degradation is at its peak, that is why both integrated systems are very abundant in the relevant literature. In addition, the light detection and ranging (LiDAR) systems are widely adopted in the literature for its capability to provide 6-DOF to mobile vehicles and autonomous robots. LiDARs are very accurate systems however they are not suitable for low-cost positioning due to the expensive initial costs. Moreover, several other techniques from the radio frequency (RF) spectrum are utilized as multi-sensor systems such as cellular networks, WiFi, ultra-wideband (UWB) and Bluetooth. The cellular-based systems are very suitable for outdoor navigation applications while WiFi-based, UWB-based and Bluetooth-based systems are efficient in indoor positioning systems (IPS). However, to achieve reliable PVT estimations in multi-sensor GNSS navigation, optimal algorithms should be developed to mitigate the estimation errors resulting from non-line-of-sight (NLOS) GNSS situations. Examples of the most commonly used algorithms for trilateration-based positioning are Kalman filters, weighted least square (WLS), particle filters (PF) and many other hybrid algorithms by mixing one or more algorithms together. In this paper, the reviewed articles under study and comparison are presented by highlighting their motivation, the methodology of implementation, the modelling utilized and the performed experiments. Then they are assessed with respect to the published results focusing on achieved accuracy, robustness and overall implementation cost-benefits as performance metrics. Our summarizing survey assesses the most promising, highly ranked and recent articles that comprise insights into the future of GNSS technology with multi-sensor fusion technique.

For precise positioning in arctic regions with Global Navigation Satellite Systems (GNSS), such as GPS and Glonass, there are a number of issues to be resolved or taken into account. The most important ones are satellite geometry, ionospheric effects and distribution of correction data. Precise Point Positioning (PPP) is currently the standard for precise offshore positioning. PPP requires a sparse global network of GNSS reference stations to estimate precise satellite orbit and clock parameters in real-time. These parameters are transmitted to users, who can compute their position with an accuracy of 0.1 m using code and carrier observations. The main disadvantage of PPP is that in general it takes a long time to converge to this accuracy, about 30–45 minutes. GNSS carrier observations are similar to pseudo ranges, but have an additional bias which needs to be estimated. Under certain conditions, these biases are integer. This property is used in what is generally referred to as Real-Time Kinematic (RTK) positioning, a technique which uses corrections generated at reference stations to make sure these ambiguities are integer. Integer Ambiguity Resolution (IAR) results in RTK accuracies of 0.01–0.03 m or better and short convergence times, ranging from instantaneous to several minutes for long inter-station distances. The disadvantage is the limited area of applicability, because RTK in general does not use precise orbits and clocks. Current research focuses on the combination of PPP and RTK techniques. The main goal is to achieve cm accuracy with a sparse network of reference stations. Another goal is to reduce convergence time. Applications are in the field of precise real-time offshore positioning in general, and subsidence monitoring and tidal estimation in particular. Results will be presented from standard PPP and PPP combined with IAR for arctic regions and from a test bed in the Gulf of Mexico. These results, although preliminary, indicate centimeter level accuracy is indeed possible, also in the arctic. We will also look at further improvements, mainly in convergence time, using all operational satellites and signals from current and future navigation systems.

This paper addresses the Kalman filter divergence problem for Global Navigation Satellite System (GNSS) in Precise Point Positioning (PPP). Continuous improvement of precision of GNSS PPP makes it widely used in various fields, however the Kalman filter divergence resolutions for PPP signal attenuating and satellite missing condition are still insufficient. Kalman filter is a recursive, linear unbiased, minimum variance method. Compared with the standard Kalman filter, the dynamic Kalman filter with attenuation factor in this paper makes up for the Kalman filter divergence problem. The effectiveness of dynamic Kalman Filter is demonstrated via matlab simulation.

Currently Global Navigation Satellite Systems (GNSS) are the primary source for the determination of absolute position, velocity and time (PVT) for merchant vessel navigation. On a ship bridge a number of systems like Electronic Chart and Information Systems (ECDIS), automatic track control systems and the Automatic Identification System (AIS) are relying on that PVT provision. Nevertheless it is well known that GNSS performance can strongly degrade due to space weather events, jamming and spoofing. Especially the availability and usage of low cost jammers, so called personal privacy devices (PPD), lead to the question on how a continuous onboard provision of PVT data can be realized in the presence of such a PPD jammer. Generally three possible solutions for that challenge are contemplated: i) the usage of a terrestrial backup system, ii) a jamming resistant GNSS receiver or iii) the integration of GNSS with other onboard self-contained navigation sensors like, e.g. speed log, gyrocompass and inertial sensors (Inertial Measurement Unit - IMU). The paper will focus on option (iii) and addresses the question whether it is possible to bridge real maritime GNSS Jamming attacks by a Kalman Filter (KF) using GNSS, inertial and Doppler velocity log (DVL). The main challenge for the KF here is, that a free inertial integration needs to be performed after a time period with already degraded GNSS signal reception. As KF are known to be optimum estimators only under the assumption of known statistics of the measurements, the question arises, whether the pseudoranges under the influence of jamming still follow a certain statistics. In order answer this question an additional lab experiment was performed and analyzed. Concretely the sensor fusion was realized by using both loosely and tightly coupled architectures which have been implemented using an Unscented Kalman Filter (UKF). The fusion algorithm is formulated in ECEF coordinate frame with quaternion attitude parameterization, where the heading observability is addressed by fusing the baseline measurements from the GNSS compass. In order to evaluate the performance of the proposed sensor fusion scheme under the influence of a real maritime jamming attack we performed an experiment in our civilian maritime jamming testbed in the Baltic, which has been established together with the German Federal Network Agency approximately 6 nmi north of the peninsula Dars. Here a PPD jammer was mounted on the monkey deck of an anchored vessel. A test vessel was maneuvering around the jammer. As a test vessel the multipurpose research and diving vessel Baltic Diver II (length 29 m, beam 6.7 m and draught 2.8 m) was used, which was equipped with three dual frequency GNSS receivers, a fiber-optic gyroscope IMU and a Doppler velocity log. The obtained results confirm that real maritime jamming attacks using PPD Jammers can be bridged by the proposed tightly coupled GNSS-IMU-DVL integration scheme resulting in horizontal position errors of less than 10m only.

This study proposes a new set of processing procedures based on the strain model and the Kalman filter (SM-Kalman) to obtain high-precision three-dimensional surface deformation time series from interferometric synthetic aperture radar (InSAR) and global navigation satellite system (GNSS) data. Implementing the Kalman filter requires the establishment of state and observation equations. In the time domain, the state equation is generated by fitting the pre-existing deformation time series based on a deformation model containing linear and seasonal terms. In the space domain, the observation equation is established with the assistance of the strain model to realize the spatial combination of InSAR and GNSS observation data at each moment. Benefiting from the application of the Kalman filter, InSAR and GNSS data at different moments can be synchronized. The time and measurement update steps are performed dynamically to generate a 3-D deformation time series with high precision and a high resolution in the temporal and spatial domains. Sentinel-1 SAR and GNSS datasets in the Los Angeles area are used to verify the effectiveness of the proposed method. The datasets include twenty-seven ascending track SAR images, thirty-four descending track SAR images and the daily time series of forty-eight GNSS stations from January 2016 to November 2018. The experimental result demonstrates that the proposed SM-Kalman method can produce high-precision deformation results at the millimeter level and provide two types of 3-D deformation time series with the same temporal resolution as InSAR or GNSS observations according to the needs of users. The new method achieves a high degree of temporal and spatial fusion of GNSS and InSAR data.

Global Navigation Satellite Systems (GNSS) has been the primary positioning solution for Unmanned Aerial Vehicles (UAVs) due to their worldwide coverage, high precision and lightweight receivers. However, GNSS is prone to electromagnetic interference and malicious assaults, including jamming or spoofing because of its low signal-to-noise ratio (SNR). To ensure the continuity and protection of UAV operations, using redundant navigation systems is essential. In recent years, the phased array radio system (PARS) has established itself as a local navigation solution. PARS is robust towards malicious assaults because of an much higher SNR than GNSS regarding directed and encrypted transmission. An essential factor of PARS is that the orientation of the radio antenna at a ground station needs to be precisely determined to obtain the correct positioning of UAVs. This paper presents a method for extending a previously proposed calibration algorithm to estimate the ground antenna orientation with an inertial navigation system (INS) aided by redundant positioning sensors (GNSS, PARS or barometer) using a multiplicative extended Kalman filter (MEKF) so that the calibration can be activated during flights whenever GNSS is available. In other words, the proposed navigation system is essentially an aided-INS which switches between two modes depending on the availability of GNSS: calibration and GNSS aiding mode when GNSS is available (Mode 1) and PARS and barometer aiding mode when GNSS is unavailable (Mode 2). Considering that the navigation system needs to include the effect of Earth’s curvature for a long-distance flight, PARS horizontal measurement and the barometer measurement were treated independently, and the navigation equations were propagated in Earth Centred Earth Fixed (ECEF) frame. The independent treatment of barometer measurement, and the propagation in ECEF frame were also beneficial when using multiple ground antennas to have a common reference point and reference frame. The proposed method was validated using data (Inertial Measurement Unit (IMU), GNSS, PARS, Pixhawk autopilot (including barometer) measurements) collected during a field test. In the validation, GNSS was made available at the middle of the flight and the calibration mode was activated for 200s. The proposed navigation system successfully estimated the precise orientation of multiple ground antennas and the navigation solutions were verified using GNSS and Pixhawk autopilot solutions as ground truth.

In complex environments, signals are inevitably subject to phenomena such as reflection, refraction, diffraction, and obstruction, which result in significant unmodeled errors like colored noise, residual systematic errors, and other special outliers. Consequently, the achievement of high-precise and high-reliable global navigation satellite systems (GNSS) precise point positioning (PPP) is not a readily solvable problem in complex environments, especially in varying terrain topography. This paper proposed the resilient inequality constrained GNSS kinematic PPP method considering the terrain topography to improve the abnormal positioning results caused by unmodeled errors. Specifically, the proposed method is composed of the adaptive inequality constraint with dynamic penalty function and the timing-varying inequality considering the terrain topography. Two representative experiments including one set of designed data and three sets of daily measured data were conducted. The results show that the proposed method can improve the positioning results resulting from the unmodeled errors while preserving the trend of the original data. Typically, the proposed method decreases the standard deviations by 2.47, 0.25, and 1.46 cm in the U direction of the three real datasets, respectively. Consequently, the proposed method exhibits prospects in precision and reliability for complex environments.

PPP-RTK which takes full advantages of both Real-Time Kinematic (RTK) and Precise Point Positioning (PPP), is able to provide centimeter-level positioning accuracy with rapid integer Ambiguity Resolution (AR). In recent years, with the development of BeiDou Navigation Satellite System (BDS) and Galileo navigation satellite system (Galileo) as well as the modernization of Global Positioning System (GPS) and GLObal NAvigation Satellite System (GLONASS), more than 140 Global Navigation Satellite System (GNSS) satellites are available. Particularly, the new-generation GNSS satellites are capable of transmitting signals on three or more frequencies. Multi-GNSS and multi-frequency observations become available and can be used to enhance the performance of PPP-RTK. In this contribution, we develop a multi-GNSS and multi-frequency PPP-RTK model, which uses all the available GNSS observations, and comprehensively evaluate its performance in urban environments from the perspectives of positioning accuracy, convergence and fixing percentage. In this method, the precise atmospheric corrections are derived from the multi-frequency and multi-GNSS observations of a regional network, and then disseminated to users to achieve PPP rapid AR. Furthermore, a cascade ambiguity fixing strategy using Extra‐Wide‐Lane (EWL), Wide-Lane (WL) and L1 ambiguities is employed to improve the performance of ambiguity fixing in the urban environments. Vehicle experiments in different scenarios such as suburbs, overpasses, and tunnels are conducted to validate the proposed method. In suburbs, an accuracy of within 2 cm in the horizontal direction and 4 cm in the vertical direction, with the fixing percentage of 93.7% can be achieved. Compared to the GPS-only solution, the positioning accuracy is improved by 87.6%. In urban environments where signals are interrupted frequently, a fast ambiguity re-fixing can be achieved within 5 s. Moreover, multi-frequency GNSS signals can further improve the positioning performance of PPP-RTK, particularly in the case of small amount of observations. These results demonstrate that the multi-frequency and multi-GNSS PPP-RTK is a promising tool for supporting precise vehicle navigation.

Ground-based GNSS ZTD/IWV estimation system for numerical weather prediction in challenging weather conditions

Structural health monitoring (SHM) is vital for ensuring the service safety of aging bridges. As one of the most advanced sensing techniques, Global Navigation Satellite Systems (GNSS) could capture massive spatiotemporal information for effective bridge structural health monitoring (BSHM). Unfortunately, GNSS measurements often contain outliers due to various factors (e.g., severe weather conditions, multipath effects, etc.). All such outliers could jeopardize the accuracy and reliability of BSHM significantly. Previous studies have examined the feasibility of integrating the conventional multi-rate Kalman filter (MKF) with an adaptive algorithm in the data processing processes to ensure BSHM accuracy. However, frequent parameter adjustments are still needed in tedious data processing processes. This study proposed an outlier detection method using a Nelder-Mead simplex robust multi-rate Kalman filter (RMKF) for supporting trustworthy BSHM using GNSS and accelerometer. In the end, the authors have validated the proposed method using the monitoring data collected at the Wilford Bridge in the UK. Results showed that the accuracy of the total dynamic vibration displacement time series has been improved by 21% compared with the results using the conventional MKF approach. The authors envision that the proposed method will shed light on reliable and explainable data processing policy and trustworthy BSHM.

Acousto-optic devices using lead molybdate and related compounds: Bell Telephone Labs Inc US patent 3617931 (1971) filed 5 May 1969

Cycle slip detection and repair based on the unmodeled-error-constrained geometry-free combining geometry-based models for a single-frequency receiver