Next Article in Journal
Construction of a Fluxgate Magnetic Gradiometer for Integration with an Unmanned Aircraft System
Next Article in Special Issue
Robust TOA-Based UAS Navigation under Model Mismatch in GNSS-Denied Harsh Environments
Previous Article in Journal
Model-Assisted Bird Monitoring Based on Remotely Sensed Ecosystem Functioning and Atlas Data
Previous Article in Special Issue
On the Recursive Joint Position and Attitude Determination in Multi-Antenna GNSS Platforms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Design a TDCP-Smoothed GNSS/Odometer Integration Scheme with Vehicular-Motion Constraint and Robust Regression

1
Department of Geomatics, National Cheng Kung University, Tainan 701, Taiwan
2
Interdisciplinary Division of Aeronautical and Aviation Engineering, The Hong Kong Polytechnic University, Kowloon, Hong Kong, China
3
Department of Land Economics, National Chengchi University, Taipei 116, Taiwan
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(16), 2550; https://doi.org/10.3390/rs12162550
Submission received: 17 June 2020 / Revised: 18 July 2020 / Accepted: 3 August 2020 / Published: 7 August 2020

Abstract

:
Global navigation satellite system (GNSS) is widely regarded as the primary positioning solution for intelligent transport system (ITS) applications. However, its performance could degrade, due to signal outages and faulty-signal contamination, including multipath and non-line-of-sight reception. Considering the limitation of the performance and computation loads in mass-produced automotive products, this research investigates the methods for enhancing GNSS-based solutions without significantly increasing the cost for vehicular navigation system. In this study, the measurement technique of the odometer in modern vehicle designs is selected to integrate the GNSS information, without using an inertial navigation system. Three techniques are implemented to improve positioning accuracy; (a) Time-differenced carrier phase (TDCP) based filter: A state-augmented extended Kalman filter is designed to incorporate TDCP measurements for maximizing the effectiveness of phase-smoothing; (b) odometer-aided constraints: The aiding measurement from odometer utilizing forward speed with the lateral constraint enhances the state estimation; the information based on vehicular motion, comprising the zero-velocity constraint, fault detection and exclusion, and dead reckoning, maintains the stability of the positioning solution; (c) robust regression: A weighted-least-square based robust regression as a measurement-quality assessment is applied to adjust the weightings of the measurements adaptively. Experimental results in a GNSS-challenging environment indicate that, based on the single-point-positioning mode with an automotive-grade receiver, the combination of the proposed methods presented a root-mean-square error of 2.51 m, 3.63 m, 1.63 m, and 1.95 m for the horizontal, vertical, forward, and lateral directions, with improvements of 35.1%, 49.6%, 45.3%, and 21.1%, respectively. The statistical analysis exhibits 97.3% state estimation result in the horizontal direction for the percentage of epochs that had errors of less than 5 m, presenting that after the intervention of proposed methods, the positioning performance can fulfill the requirements for road level applications.

Graphical Abstract

1. Introduction

To realize safe intelligent transportation systems (ITS), the demand for land vehicle navigation has been rapidly increasing. Global Navigation Satellite System (GNSS) is commonly used to provide global three-dimensional (3D) navigation with good long-term accuracy. It is currently the principal technique for providing globally referenced positioning without other aiding information. In kinematic positioning, code measurement is preferred for position estimation because it does not contain integer ambiguities, and meter-level accuracy can be expected from the single-point positioning (SPP) technique in open-sky environments [1]. However, unsatisfactory navigation could occur because of its low precision and the multipath effect [2]. Further, range errors due to the multipath effect could be introduced up to half a code chip (150 m for GPS C/A code) because of the code tracking error [3]. On the other hand, the maximum tracking error for carrier-phase measurements due to multipath interference from a reflected signal of the same amplitude as the direct signal is a quarter of a wavelength (e.g., 4.76 cm for GPS L1) [4]. The carrier-smoothed-code method is generally used for dynamic navigation applications, which combines a series of delta-phase and pseudorange into a single noise-reduced measurement. Here, a pseudorange provides absolute position information and an incremental carrier-phase measurement provides displacement information [5,6]. The advantage of using this method is the high spatial variation in multipath errors by averaging out most of the code multipath errors without significantly increasing the cost and computation load [4]. However, it is difficult to maintain pre-supposition of this phase-smoothing method in the range domain for continuous tracking of the same satellites in urban areas. Every time a signal outage occurs, all of the smoothed pseudo-range information is lost, and the accuracy of the combined measurements reverts to the nominal unsmoothed level [4,7]. On the other hand, the position-domain phase-smoothing method utilizes time-differenced carrier phase (TDCP) measurements as a constraint of the current position relative to the previous position. Theoretically, this epoch-to-epoch position change can be determined to the level of TDCP accuracy through the observables directly related to the position increment [7]. In addition to the enhancement of redundancy, the advantage of such information in the position domain over the range domain is that carrier-phase measurements from the same satellites are required, only since the immediately previous time epoch and not several previous epochs [5,7], which increases the availability of phase-smoothing. Nonetheless, a problem with the phase-smoothing methods including both domains is the frequent occurrence of cycle slip. Limited to the accuracy of cycle slip detection, especially for single-frequency GNSS receivers, the ambiguities emerging due to undetected loss of signals will be introduced into the incremental carrier-phase information and further influence the performance of estimation. Notably, phase-smoothing does not mitigate the effects of non-line-of-sight (NLOS) reception, in which the signal is received only via reflections. The range measurement error, caused by the NLOS reception, is equal to the additional path delay, which is the difference between the length of the path taken by the reflected signal, and the blocked-direct path, between the satellite and the receiver antenna [4,8]. In this situation, because the code and the carrier are affected in the same manner, the phase-smoothing techniques, including TDCP, cannot improve the position estimation. In these GNSS-degraded environments, the combination of multipath/NLOS contamination significantly challenges the reliability of GNSS positioning.
Over the past years, many techniques have been developed to address this contamination problem. In the antenna-design technique, based on polarization discrimination, a polarization-sensitive antenna can reduce the magnitude of signal tracking errors caused by the multipath effect. The usage of carrier-power-to-noise-density ratio measurements based on this antenna could enhance the performance of NLOS detection, except for doubly reflected signals between the transmitter and receiver [4]. However, due to the cost and size limitations, this technique is less effective in the patch-antenna product for land-vehicle applications [4,9,10]. Other designs, such as using choke-ring antennas and multiple antennas, are expensive approaches as well [4,9]. In the post-receiver technique, consistency checking is the conventional approach based on the fact that a multipath/NLOS-contaminated measurement is inconsistent with normal measurements [4,11,12,13]. This method can improve the GNSS performance under open environments with the emergence of multi-constellation GNSS, in which the number of visible satellites can be between 15 and 30, especially in the Asia Pacific region where is easy to receive GPS, GLONASS, Beidou, Galileo and QZSS signals at the same time [13,14]. However, in the meantime, less satellites would be received due to signal blockage in urban areas, while more reflected signals can be received in such environments [15]. Similar to the concept of receiver autonomous integrity monitoring (RAIM), the consistency checking fails in GNSS-challenging environments where a large proportion of the signals are NLOS-contaminated [13]. Although, using the random sample consensus (RANSAC) method can improve its performance in case of multiple outliers [16], it is a time-intensive process as it involves an iterative selection of the observation subset. Compared with the previous fault detection and exclusion (FDE) based methods, robust regression as a quality assessment method adaptively adjusts the weightings of the measurements [17,18]. This method is based on the minimization of cost function depending on the least-squares residuals, which represent an approximation of the measurement errors [19]. The magnitude of the residuals can be used to verify the consistency between the measurements and reweight each measurement accordingly. On the basis of this concept, robust regression has the capability of more efficient computation when multiple faulty measurements emerge, compared with RAIM using the RANSAC. Reference [19] further extends this approach and implements in extended Kalman filter (EKF) to enhance the performance when outliers exist. Nonetheless, the performance of reweighting method could be unsatisfied when abnormal measurements are in majority, because the resulting residuals are significantly influenced by contaminated signals. In the multi-sensor aiding technique, some of the methods utilize additional information from other sensors, such as LiDAR, camera, or 3D-building models to detect the NLOS signals [20,21,22,23,24,25,26]. However, it cannot discriminate multipath-contaminated signals from the received group, and suffers from cost and computational challenges, thereby, limiting their application in current consumer products. One of the most popular methods in this technique is integrating with aiding sensors to improve the natural GNSS limitations (i.e., satellite availability and spatial geometrical strength), thereby, enhancing the accuracy and continuity for standalone GNSS. Among these sensors, the integration of inertial navigation system (INS) and GNSS is the most common scheme applied in navigation applications due to its complementarity [27,28,29,30,31]. Nonetheless, it is well-known that the error of low-cost INS grows rapidly, due to large uncertainties in the sensor output, including having high drift rates, being vulnerable to non-linear error behavior, and sensitivity to environmental effects, such as temperature variations and shocks [31,32,33]. There are some methods for improving the accuracy of the low-cost INS, the most important being stochastic error modeling [31,32,33]. However, for low-cost INS, it is difficult to find an adequate error parameter due to statistically inconsistent in modeling stochastic errors (e.g., Allan variance method) [33]. In addition, the quality of the initial alignment is also another factor influencing the performance [31,33], which is hard to ensure in civil navigation application. Also, this GNSS/INS combination does not always satisfy the product design due to the high power consumption, memory storage, and increased cost. Another popular sensor for aiding navigation is the odometer, which is a self-contained velocity sensor available as a standard component in antilock braking systems (ABS). It generates digital pulses for each revolution of the wheel and allows estimation of the distance covered by a vehicle [34]. For modern land-vehicle designs, this information can be obtained through On-board Diagnostics (OBD II), a standardized electrical data interface that communicates with a vehicle’s computer. In urban areas, it generally provides more accurate and frequent velocity measurements than GNSS [35,36]. However, velocity measurements from the odometer are obtained in the vehicle frame (v-frame), and not the navigation frame, such as earth frame (e-frame). Therefore, such velocity information cannot be directly used in the absolute coordinate system. Another drawback is the position drift phenomenon due to error accumulation of the measuring distance. Therefore, performance of only the navigation task over a long time is not adequate.
It is clear that the current automotive GNSS solutions cannot easily meet the tradeoff between cost and performance in GNSS-challenging environments. Therefore, this study presents a GNSS/odometer integration scheme with odometer-aided constraints and robust regression for navigation applications. A state-augmented EKF-based GNSS estimator incorporated with the TDCP measurements is implemented for the online-smoothing of the state estimation in the position domain. The v-frame velocity information obtained from the in-vehicle odometer through the OBD II interface is integrated with the proposed vehicular-motion constraints, fault detection and exclusion (FDE), and dead reckoning (DR) aiding to enhance the performance without increasing the cost. Moreover, robust regression as an adaptive weighting method is combined with odometer-aided measurement, providing more reliable weighting indicators for state estimation in a scenario with increased multi-false measurements, without significantly burdening the navigation system with increased computation. This research investigates the feasibility of a cost-effective onboard solution for land-vehicle navigation applications without integrating INS and complex statistical models.

2. Methods

Figure 1 illustrates the GNSS/odometer integration scheme with weighted least squares (WLS)-based robust regression, which includes the measurements implemented in this paper. In the proposed design, GNSS acts as a major sensor, comprising pseudorange, pseudorange rate, and TDCP measurements, for providing absolute positioning state (position and velocity). Whereas, the odometer with the forward speed acts as the constraint in the position and velocity domain. Those measurements will be processed by the WLS-based robust regression to reweight the scale indicators adaptively, and the state-augmented EKF integrates all of the inputs for the state estimation. In the following, the design of state-augmented EKF corresponding to delta measurement (TDCP) for estimating the position change is introduced in Section 2.1. Section 2.2 describes the measurement model, including conventional measurements and TDCP. Section 2.3 introduces the measurement model for the odometer, and the corresponding designs of vehicular-motion constraint and FDE mechanism in this study are also included. The robust regression based on the odometer aiding method proposed in this study is in Section 2.4. The process flow comprised of all of the proposed methods will be illustrated afterward.

2.1. State-Augmented EKF Design

In this study, the position and velocity error states of GNSS positioning are estimated using EKF. Conventionally, the EKF design presents only the system error at a particular time. Whereas, a delta phase or delta position represents the integrated velocity over time [5]. To sufficiently use such delta measurement for estimating the change in the receiver position, the conventional system model should be modified and extended to include the estimated position and receiver clock bias from the last epoch as follows [5,37]:
δ x ^ = [ δ r u e ( t k ) δ v u e ( t k ) δ r u e ( t k 1 ) δ b u ( t k ) δ b d ( t k ) δ b u ( t k 1 ) ] 1 × 12 T
Φ p 1 v p 0 t = [ Φ p 1 v 0 6 × 3 0 6 × 3 [ I 0 3 × 3 ] 3 × 6 0 3 × 3 0 3 × 3 0 3 × 6 0 3 × 3 [ Φ t 0 [ 1 0 ] 1 × 2 0 ] 3 × 3 ] 12 × 12
Q p 1 v p 0 t = [ Q p 1 v 0 6 × 3 0 6 × 3 0 3 × 6 0 3 × 3 0 3 × 3 0 3 × 6 0 3 × 3 [ Q t 0 1 × 1 0 1 × 2 0 1 × 1 ] 3 × 3 ] 12 × 12
with
Φ p 1 v = [ 1 0 0 Δ t 0 0 0 1 0 0 Δ t 0 0 0 1 0 0 Δ t 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 ] 6 × 6 ,   I = [ 1 0 0 0 1 0 0 0 1 ] 3 × 3 Q p 1 v = [ q p Δ t + 1 3 q v Δ t 3 0 0 1 2 q v Δ t 2 0 0 0 q p Δ t + 1 3 q v Δ t 3 0 0 1 2 q v Δ t 2 0 0 0 q p Δ t + 1 3 q v Δ t 3 0 0 1 2 q v Δ t 2 1 2 q v Δ t 2 0 0 q v Δ t 0 0 0 1 2 q v Δ t 2 0 0 q v Δ t 0 0 0 1 2 q v Δ t 2 0 0 q v Δ t ] 6 × 6
where δ x ^ is the error state vector including the position errors   δ r u e , velocity error δ v u e , receiver clock bias δ b u , and receiver clock drift δ b d ; t k and t k 1 denote the corresponding epochs of the error states; each of δ r u e and δ v u e is a 3 × 1 vector, where the superscript e denotes a vector in the e-frame (earth–centered-earth-fixed coordinate system); Φ p 1 v and Q p 1 v denote the transition matrix ( t k to t k + 1 epoch) and the corresponding process noise matrix under assumed constant velocity conditions, for the current time; q p and q v are the spectral densities of the above process noise sequences; Φ t and Q t denote the transition matrix and process noise matrix for the receiver clock state; Δ t is the time interval between adjacent epochs; Φ p 1 v p 0 t and Q p 1 v p 0 t are the augmented models from Φ p 1 v and Q p 1 v to incorporate with the delta-measurement update; I is the identity matrix, corresponding to each component of the current position error state. The subscript   p 1 ,   v , and p 0 denotes that the matrix contains the components corresponding to the current position, velocity, and previous position error state. This design both supports the previously defined dynamic system for the random walk model, and transfers the current position error vector to the previous spot in the error state vector during the propagation. In this design, the delta observables are directly related to the position increment as a constraint of the current position, relative to the previous position. Based on the phase-smoothing approach in the position domain, this EKF design can maximize the effectiveness of such measurements without assumptions about the dynamics of the vehicle [5] and control the growth of position errors. The correspondingly adopted measurement model in this study is described in the following subsection.

2.2. Measurement Model for Pseudorange, Doppler Shift, and TDCP

The models for GNSS measurements, including pseudorange, Doppler shift, and TDCP, which are related to the state of a vehicle are considered in the SPP mode in this study. The models of pseudorange and pseudorange rate measurement are well-documented in the literature, which can be found in [38]. The matrix form of the fundamental observations corresponding to the designed EKF can be expressed as,
H c o d e = [ e i ( t k ) 0 n × 3 0 n × 3 1 n × 1 0 n × 1 0 n × 1 ] n × 12
H d o p p l e r = [ 0 n × 3 e i ( t k ) 0 n × 3 0 n × 1 1 n × 1 0 n × 1 ] n × 12
H = [ H c o d e H d o p p l e r ] 2 n × 12
δ z = [ ρ ˜ i ρ i ρ ˙ ˜ i ρ ˙ i ] 2 n × 1 = H δ x ^ + ε 2 n × 1
R c o d e / d o p p l e r = σ 0 2 [ 1 s i n 2 ( e l 1 ) 0 0 1 s i n 2 ( e l i ) ] n × n
with
e i ( t k ) = 1 ρ i ( t k ) [ ( r u , x e ( t k ) r s , x e , i ( t k ) ) ( r u , y e ( t k ) r s , y e , i ( t k ) ) ( r u , z e ( t k ) r s , z e , i ( t k ) ) ] 1 × 3
where H c o d e , H d o p p l e r , ρ ˜ i , ρ i , ρ ˙ ˜ i , ρ ˙ i , R c o d e / d o p p l e r , σ 0 2 , e l i , and e i ( t k ) denote the design matrices of the pseudorange and pseudorange rate, the observed pseudorange, geometric range between the receiver and the satellite, the observed pseudorange rate, geometric pseudorange rate, weighting matrix of the pseudorange/pseudorange rate measurement, a priori variance factor, the elevation angle of the ith satellite with respect to the rover, and LOS unit vector from the satellite to the receiver at t k , respectively. The design of the weighting matrix of pseudorange rate is similar to R c o d e , but with a different σ 0 2 .
The TDCP measurements, which are the delta-phase measurements from the same satellites between consecutive epochs, are conventionally used to estimate the average velocity [39,40,41,42,43]. In other words, the difference between the previous and current positions with mm-level noise can be observed through such phase differences. It could be partially observable if fewer than four satellites are continuously available in the absence of cycle slip. The advantage of this approach over conventional techniques is that carrier-phase measurements from the same satellites are merely available since the previous time epoch, rather than several seconds, to average out the code multipath error [5,37]. The original measurement equation of the carrier phase is expressed as follows:
λ Φ i = ρ i + c ( δ b u δ b s i + δ d e p h i δ d i o n o i + δ d t r o p i ) + λ N i + η i
where λ , Φ i , and N i represent the wavelength, measured carrier phase in cycles, and integer ambiguity, respectively. The term η i includes a multipath/NLOS-contaminated source and receiver noise for the carrier phase. The difference between the carrier-phase measurements at consecutive epochs t k and t k 1 is described as follows:
λ Δ Φ i = λ · [ Φ i ( t k ) Φ i ( t k 1 ) ] = Δ ρ i + c ( Δ δ b u + Δ δ b s i + Δ δ d e p h i Δ δ d i o n o i + Δ δ d t r o p i ) + Δ η i
where Δ indicates the differencing operation. Under the condition of no cycle slip, the integer ambiguity is eliminated by time differencing. Other common error sources, which include satellite clock bias, ephemeris error, tropospheric error, and ionospheric error, vary slowly within a sampling, and are, Δ δ b s i ,   Δ δ d e p h i , Δ δ d i o n o i , and Δ δ d t r o p i are negligibly small. Notably, the transition between the adjacent ephemeris sets causes the discontinuity of Δ δ d e p h i , especially for Δ δ b s i , when different sets are respectively selected between consecutive epochs for TDCP measurements [43]. Selecting the same ephemeris can prevent this problem during the transition. In transforming (10) into a relation of receiver position change, Δ ρ i has to be recombined as [43],
Δ ρ i = Δ D Δ g [ e i ( t k ) × Δ r u ]
with
Δ g = [ e i ( t k ) × r u e ( t k 1 ) ] [ e i ( t k 1 ) × r u e ( t k 1 ) ] , Δ D = [ e i ( t k ) × r s e , i ( t k ) ] [ e i ( t k 1 ) × r s e , i ( t k 1 ) ] ,  
where r u e and r s e , i are the receiver position vector and satellite position vector, respectively, including the x, y, and z components in the e-frame. The term Δg reflects the orientation change of the LOS vector in the relative satellite–receiver geometry, and ΔD represents a change in the range and is proportional to the average Doppler shift caused by the relative satellite–receiver motion along the LOS vector [44]. The derivation of (11) is presented in detail in [39,43]. By substituting (11) into (10), the TDCP measurement equation becomes (12), and the corresponding measurement model is expressed by (13) and (14). The weighting matrix of the TDCP measurements, R TDCP , is also similar to the design of (17) with a different σ 0 2 .
λ Δ Φ i Δ D + Δ g = [ e i ( t j ) × Δ r u e ] + c Δ δ b u + Δ η i
H TDCP = [ e i ( t k ) 0 1 × 3 e i ( t k 1 ) 1 0 1 ] m × 12
δ z TDCP = [ Δ Φ i Δ ρ i ] m × 1
The effectiveness of TDCP is affected by cycle slip. The cycle-slip detection method applied in this study is “Doppler integration” for a single-frequency receiver. This method uses an observed Doppler shift to predict the growth in a carrier phase observation as an indicator of a cycle slip, whose fundamental equation is presented in [38]. Nonetheless, due to the vehicular dynamics and the environment, the signal tracking quality of Doppler shift strongly affects the accuracy of detection. When a cycle slip occurs with false detection, new ambiguities are introduced in delta-phase measurements. On the other hand, the path delay of the carrier phase measurements due to NLOS reception could be similar to the pseudorange error, as mentioned previously in [8]. These errors could be introduced into carrier phase ambiguities and further degrade GNSS positioning if the phase-smoothing approach is applied. In this situation, standalone GNSS positioning cannot satisfy the requirement of accuracy and reliability in vehicular navigation. This study proposes the aiding from odometer velocity via OBD II to enhance the robustness of positioning in case of occurrence of such a situation or frequent cycle-slip. The corresponding aiding model based on the proposed EKF design is introduced in the following subsection.

2.3. Measurement Model for Odometer Observation with Vehicular-Motion Constraint

The odometer-derived velocity indicates that unless a vehicle is off the ground or slides on the ground, the velocity of the vehicle in the plane perpendicular to the forward direction is nearly zero [44,45]. Because an odometer provides the forward speed of a vehicle, the measurement is expressed as follows:
v v = [ v 0 0 ] T
where v v denotes the velocity measurement in the v-frame, and the lateral and vertical speeds are presumed to be zero. Because (15) need not be obtained in the e-frame, it cannot be used directly in the state estimation. In this study, based on the fact that v is the observation as the resultant speed of x, y, and z components for the velocity state in (1), the corresponding measurement model can be formulated as,
H o d o = [ 0 1 × 3 ( v 3 d ) 1 2 v x e ( v 3 d ) 1 2 v y e ( v 3 d ) 1 2 v z e 0 1 × 3 0 1 × 1 0 1 × 1 0 1 × 1 ] 1 × 12
δ z o d o = [ v v 3 d ] 1 × 1 = H o d o δ x ^ + η o d o
R o d o = [ σ o d o 2 ] 1 × 1
with
v 3 d = ( v x e ) 2 + ( v y e ) 2 + ( v z e ) 2
where H o d o is the design matrix for the odometer; η o d o comprises the measurement noise and abnormal sensing such as sliding on the ground; v x e , v y e , v z e are the initial GNSS velocity at the current epoch in the three directions; v 3 d is the initial resultant speed. The determination of σ o d o 2 can be given by a pre-calibrated process (e.g., comparison between the resultant speed of the reference system and the OBD II odometer). It is worth mentioning that this model only constrains the combined velocity for an estimation that is non-directional, and could face the local minimum problem during state estimation. Due to the lateral velocity of a land vehicle is nearly zero under normal conditions, it can act as an additional aiding source to enhance such state estimation. The relationship between the vehicle frame and the local-level frame (l-frame) in the horizontal direction is depicted in Figure 2, and the corresponding equation given by,
H l a t e r a l = [ 0 1 × 3 A ( t k 1 ) B ( t k 1 ) C ( t k 1 ) 0 1 × 3 0 1 × 1 0 1 × 1 0 1 × 1 ] 1 × 12
δ z l a t e r a l = [ 0 v l a t e r a l ] 1 × 1 = H l a t e r a l δ x ^ + η l a t e r a l
R l a t e r a l = [ σ l a t e r a l 2 ] 1 × 1
with
v l a t e r a l = A ( t k 1 ) × v x e ( t k ) + B ( t k 1 ) × v y e ( t k ) + C ( t k 1 ) × v z e ( t k ) , A = c o s ( θ ) × s i n ( Λ ) + s i n ( θ ) × s i n ( φ ) × c o s ( Λ ) , B = c o s ( θ ) × c o s ( Λ ) + s i n ( θ ) × s i n ( φ ) × s i n ( Λ ) ,   C = s i n ( θ ) × cos ( φ )
where v l a t e r a l is the lateral velocity derived from the velocity state estimation; η l a t e r a l is the combination of velocity estimation error and heading error; H l a t e r a l is the design matrix of the lateral constraint based on the transformation of the e-frame, the l-frame, and vehicle frame; θ , φ , and Λ denote the vehicular heading, latitude, and longitude angles at t k 1 epoch. The determination of σ l a t e r a l 2 is based on the uncertainty of the estimated heading, which is an empirical coefficient in this study. θ is derived from the velocity state transformed into the l-frame and is thus vulnerable to the quality of GNSS measurement and significant vehicular motion. To secure the performance of the lateral constraint, the moving average filter (MAF) and the detection of heading condition are designed before (19)–(21) update EKF, which is formulated as,
θ = θ ¯ i 1 + θ i θ i 1 γ ,   | θ i θ ¯ i 1 | d θ t h r e   &   β > γ
with
θ ¯ = { 1 n i = 1 β θ i , β = γ i = 1 β θ i , β < γ
where θ ¯ i 1 , θ i , d θ t h r e , γ , and β denote the averaged heading angle at i 1 epoch, heading angle derived from the estimated velocity at the current epoch, threshold of heading change between the current heading and the previously averaged angle, window size, and the number of θ i used in MAF. The implementation of the proposed lateral velocity constraint in (19)–(22) further enhances the directivity of smoothing between adjacent epochs, because a priori information is used. In this design, when the heading change is less than d θ t h r e , and β exceeds γ , the lateral constraint will use the averaged heading based on the outcome of (22); otherwise, the lateral constraint will be closed until the counts of qualified θ i reach γ . The intervention of d θ t h r e is used to assure that the designed constraint will not be activated during significant vehicular motion, such as a turning maneuver. Notably, the calculation of angular distance in (22) must consider the vibration of angular conversion from the estimated velocity between adjacent epochs, such as the case of | θ i θ ¯ i 1 | > 180° or < −180°, averting a false output/determination in (22).
The advantage of the odometer-based methods mentioned above is that they provide additional measurements (constraints) and improve the state estimation result under the contamination of multipath/NLOS reception. However, if a limited number of tracked satellites with bad geometry exists and false GNSS measurements update EKF, unsatisfactory estimation could still emerge. Using the similar concept of (19)–(22), odometer-based FDE can be implemented. The velocity, based on an odometer, is a forward speed in the v-frame, as given in (15), following the forward direction of the vehicle. If the result of the position state estimation from GNSS presents an abnormal motion in the lateral, the estimation result is dominated by unhealthy signals. Such an FDE method can be formulated as follows:
Δ d l a t e r a l = A ( t k 1 ) × ( r u ,   x e ( t k ) r u , x e ( t k 1 ) ) + B ( t k 1 ) × ( r u ,   y e ( t k ) r u ,   y e ( t k 1 ) ) + C ( t k 1 ) × ( r u , z e ( t k ) r u , z e ( t k 1 ) ) , δ x ^ { f e e d b a c k ,   Δ d l a t e r a l s × v ( t k 1 ) f e e d b a c k   r e j e c t e d , Δ d l a t e r a l > s × v ( t k 1 )
where Δ d l a t e r a l , s , and v k 1 indicate the displacement between adjacent epochs in the lateral direction of the vehicle, scale coefficient, and instantaneous velocity from an odometer at the last GNSS epoch, respectively. As shown in Figure 2, the concept of (23) is based on the fact that, theoretically, for a land vehicle, the displacement in the lateral direction between consecutive GNSS epochs will not surpass the resultant speed in the forward. While, s is a variable coefficient depending on different vehicular conditions, such as low or high dynamics, coping with various scenarios in urban environments. If an unsatisfactory state estimation is detected and excluded by (23) or a GNSS outage occurs, the information, including odometer-based velocity and smoothed-heading from the last epoch can provide the DR solution. Although this solution provides only a 2-D solution with a constant height state, which is vulnerable to variations in height and heading, such design can bridge the gap in positioning and maintain the estimation accuracy over a short time.
Another additional aid from an odometer is the zero-velocity detection. When a vehicle is stationary, the velocity state estimation is nearly zero. However, in the densely packed urban areas, NLOS/multipath-contaminated reception becomes more severe for the duration that the vehicle is stationary. Such continuous faulty-signal reception not only degrades the velocity estimation but also misleads the position state due to the path delay. These measurements falsely update the EKF during the static motion, causing the EKF to diverge. To address this situation, the use of an odometer can provide precise constraints by assigning the velocity to be zero in any direction, thereby preventing error growth. This zero-velocity constraint is expressed as,
H o d o = [ 0 3 × 3 I 0 3 × 3 0 3 × 1 0 3 × 1 0 3 × 1 ] 3 × 12
δ z o d o = [ 0 3 × 1 v u e ] 3 × 1
where the identity matrix I corresponds to each component of the velocity state in (1). Moreover, during the activation of zero-velocity constraint, GNSS measurements will not update EKF and fix the position state and heading output, thereby preventing an erroneous position change. This position-fix design can be formulated in (26), where v t h r e denote the velocity threshold of determining the stationary motion. This mechanism can maintain stability for the state and its covariance matrix of EKF, and the accuracy of the solution. The flowchart depicting the utilization of zero-velocity constraints is presented next (Figure 3):
r u e ( t k ) { f i x   p o s i t i o n   e q u a l   t o   r u e ( t k 1 ) ,   v ( t k ) v t h r e p o s i t i o n   s t a t e   e s t i m a t i o n , v ( t k ) > v t h r e

2.4. Robust Regression with Odometer Aiding

WLS estimation is more susceptible to outliers compared with KF-based estimation. Robust regression, based on WLS residuals analysis, can adjust the weight measurements probably affected by multipath/NLOS reception whose corresponding residuals are large [18,19]. Such a process is a generalization of the maximum likelihood estimator (M-estimator), and is applied in this study. Considering that the state estimation for WLS is an iterative process (until the convergence is reached), the weightings depending on the residuals are adaptively reweighted in each iteration. The weighting function of this robust regression in each iteration can be expressed as follows [18,19]:
( w i ) j + 1 = { ( w i ) j ,   ( w i ) j α | ( r i ) j σ ^ 0 | , | ( r i ) j | σ ^ 0 α | ( r i ) j | σ ^ 0 > α
where ( w i ) j denotes the diagonal element of the weighting matrix corresponding to the ith satellite measurement at the jth iteration; σ ^ 0 is the a posteriori scale factor given by median absolute deviation due to the spread of the residuals;   α is the tuning constant and is set to 1.345 based on [46]; α / | ( r i ) j / σ ^ 0 | is the weighting scale. Using (27), the residuals are considered as an approximation of the measurement errors, and the results until WLS convergence are implemented to adaptively adjust the measurement covariance matrix used in EKF estimation [19]. A potential limitation of this method emerges if abnormal measurements are in majority, as mentioned previously. To improve the performance of robust regression under the emergence of multiple faulty measurements, this study proposes the use of accumulated distance based on the odometer information as a positional-displacement constraint to the WLS filter during the robust regression operation. When the vehicle does not change its heading orientation significantly, the moving distance from the wheel is similar to the difference between the corresponding position states of the adjacent epochs. Considering the error accumulation, due to the non-linear acceleration and the driving path challenges this assumption, such distance information is used only as an aided constraint for obtaining a weighting scale in (27), rather than as a direct measurement update in EKF estimation. The wheel information from the OBD II interface is transformed into instantaneous velocity; therefore, the accumulated distance can be calculated as follows,
Δ d 3 d = 1 ξ κ = 1 ξ v κ
where Δ d 3 d , ξ , and v κ indicate the averaged displacement between consecutive GNSS epochs, odometer rating between GNSS consecutive epochs, and instantaneous velocity measurement through the OBD II interface. Similar to (16)–(18), this measurement model of position constraint for the WLS filter can be formulated as,
δ x ^ W L S = [ δ r u e ( t k ) δ b u ( t k ) δ r u e ( t k 1 ) b ( t k 1 ) ] 1 × 8 T
H W L S = [ r 3 d 1 2 Δ r x e r 3 d 1 2 Δ r y e r 3 d 1 2 Δ r z e 0 r 3 d 1 2 Δ r x e r 3 d 1 2 Δ r y e r 3 d 1 2 Δ r z e 0 ] 1 × 8
δ z 3 d = [ Δ d 3 d r 3 d ] 1 × 1 = H W L S δ x ^ W L S + η Δ d
with
r 3 d = ( Δ r x e ) 2 + ( Δ r y e ) 2 + ( Δ r z e ) 2
where H W L S , δ x ^ W L S , η Δ d , and r 3 d are the design matrix of averaged displacement, error state, including position and receiver clock bias for WLS filter, combined errors comprising sliding, heading change and noise, and initial distance derived from the change of initial rover position between consecutive epochs; Δ r x e , Δ r y e , and Δ r z e denote the difference between the initial position vectors in the three directions of the e-frame. The measurement model in the WLS filter is similar to the design in EKF, as mentioned previously, and is not differentiated here. The benefits of odometer aiding for robust regression, include the improvement in observability of geometry, which combines the ground-based (odometer) and space-based (GNSS) sources, and the enhanced connection of relative-position change in addition to the TDCP measurements. Those contributions enhance WLS filter to find more accurate results; therefore, a priori weighting indicators in (27) will be more reliable to implement in EKF estimation.
The flowchart depicting the combination of the proposed methods in this study is shown in Figure 3. The data flow from process 1 to process 6 aims to determine the condition of combining TDCP measurements to update EKF (process 6) according to the result of receiving carrier-phase measurements (process 3) and cycle slip detection (process 4). After odometer measurements are introduced, the priority is the detection of stationary motion (process 8), since the strong intervention from zero velocity constraint updates EKF without any GNSS measurements (process 9 and 13) and fix the position output (process 14), ending the whole process flow for the current epoch. While the vehicle is not stationary, odometer-aided robust regression process based on WLS will be conducted to obtain weighting scales (process 11), and the corresponding GNSS measurements, odometer-aided velocity constraint (based on the resultant speed, process 10), and lateral velocity constraint (process 12) will be used to update EKF accordingly (process 13). After the position and velocity states have been estimated, and the heading (derived from the velocity states) does not change significantly (process 15), odometer-based FDE in (23) will be activated (process 16). Once an unsatisfactory estimation is detected, the odometer-based DR solution will bridge the gap during the estimation (process 17). Otherwise, MAF will be activated (process 18) and produce the averaged heading for other processes (process 12, 15, and 17) for the next epoch.

3. Experiment

In the experiment, the reference system uses a high tactical-grade INS (NovAtel SPAN® LCI, Canada) with a multi-band reception antenna (NovAtel GPS-GGG-703-HV), and the evaluated system is an automotive-grade Allystar EVK-2024 module with a single-band active patch antenna (Allystar AGR 6301, China); these sensors are mounted on the same vehicle platform (Figure 4). Table 1 lists the specifications of the reference system. Table 2 shows the performance of the OBD II odometer applied in this study, which is evaluated from the other individual tests in an urban area. To evaluate the performance with the proposed methods, an experimental route (red trajectory) is designed under a GNSS-challenging environment (Tainan city), as shown in Figure 5. The referenced position solution from tactical-grade INS is based on post-smoothing processing with full constellations and frequencies receiving, whereas the evaluated EVK-2024 module uses two satellite constellations (GPS, Beidou) with single-frequency (L1 band) reception in the SPP mode. All of the results are synchronized with GPS time, and the relationship between the reference system and the EVK-2024 antenna (as shown in Figure 4) are measured by the surveying tapes. The performance analysis can be implemented after the referenced results are projected to the evaluated antenna.
The total experiment duration is approximately 40 min. The tested results are based on forward-filter estimation with the proposed EKF design and without any post-smoothing processing. The scenarios A–E are selected to discuss the contribution of the proposed methods, where are multipath/NLOS-contaminated environments, and their locations are represented by white blocks in Figure 5. The corresponding street-view images of these areas are shown in Figure 6a–f. Moreover, as shown in Figure 7a,b, respectively, the position dilution of precision (PDOP) increases, and the number of visible satellites significantly decreases after the vehicle approaches urbanized areas, although two constellations are used. The observability of healthy satellite signals is frequently unsatisfactory during the operation. For analysis, the implemented methods are classified into proposed odometer-aided constraints and robust regression, i.e., method A, and method B, respectively. The tested results include code-based measurements only (conventional GNSS), code/TDCP measurements (TDCP-aided GNSS), code/TDCP measurements with odometer-aided constraints (TDCP-aided GNSS with method A), code/TDCP measurements with robust regression (TDCP-aided GNSS with method B), and code/TDCP measurements with both methods (TDCP-aided GNSS with both methods). The evaluated solutions are colored in cyan, orange, blue, green, and black, in order for the following analysis for the trajectory and the line chart of errors. Moreover, to evaluate the performance of these solutions during a GNSS outage (area I of scenario E, Figure 6e), the solutions without method A will use assumptions, such as constant velocity from the last epoch to bridge the gap during GNSS positioning. Whereas, the others with method A perform odometer-aided DR in this experiment. For the error analysis, the results from the reference system at the first epoch in this experiment is selected as the origin of the l-frame, and all of the estimated results are transformed from the e-frame into the l-frame accordingly in the following.

4. Results and Discussion

In this section, the analysis of the proposed methods will be categorized into separate subsections. Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 show the trajectory and the corresponding positioning errors in the E, N, and U directions of the l-frame for the discussed scenarios. In addition, the statistical analyses are shown in Table 3, Table 4, Table 5 and Table 6 accordingly.

4.1. Performance Analysis of TDCP-Aided GNSS

The contribution of TDCP aiding could be mainly found in a multipath/NLOS-contaminated environment. It can be shown in Figure 8 (scenario A, time 520–625 s), the TDCP-aided GNSS solution produced more accurate results than the conventional GNSS solution. Pseudorange measurements, suffering from such contamination, presented unsatisfactory estimation with approximately 20 m of maximum error in the horizontal direction and 22 m in the vertical direction. Whereas, TDCP-aided GNSS showed approximately 12 m and 15 m of maximum error, mitigating some of the influence from unhealthy signals. Moreover, similar performance was also found in Figure 11 and Figure 12 (scenarios D and E, during time 1550–1690 s and 1965–2030 s), where the TDCP-aided GNSS solutions were able to provide relatively smoothed results when the trajectory of conventional GNSS was misled into the building area.
Most of the time, the use of TDCP measurements with state-augmented EKF lessened the multipath effect; however, it is similar to other phase-smoothing approaches that lead to unsatisfactory results. For example, in Figure 9 (scenario B), the TDCP-aided GNSS solution (orange line) exhibited deteriorated position estimation, where the trajectory was located at the building. Its position error was almost the same as the unsmoothed level, where the E-direction position error reached up to 12.53 m at time 815 s. The cause was attributed to the emergence of frequent cycle-slip due to the signal path delay, and the remaining code-based measurements were significantly affected by the multipath/NLOS effect. Due to the locations of the surrounding buildings, degraded signals mainly came from the east-west direction, causing position errors in the E-direction analysis in scenario B. Another specific situation in urban areas is the duration of the stationary motion, where the abnormal signals were continuously received due to the reflection. The results are presented in Figure 8, Figure 9 and Figure 10 (scenario A–C, during time 550–610 s, 835–863 s, and 907–954 s), respectively (highlighted by yellow circles). The degraded estimation could be obviously observed, particularly in the vertical direction. Taking scenario B as an example, despite using the TDCP measurements, the maximum error was still up to 15 m in height, and the maximum error is 13 m in the E-direction, which is similar to the conventional GNSS. The reason is that the multipath/NLOS reception and the limited satellite geometry result in a deviated estimation in these challenging environments.
In addition to the influence of multipath/NLOS signals, the potential problem of using the phase-smoothing approach would be the occurrence of cycle-slip detection. When unstable tracking qualities frequently occur for Doppler-shift measurements, the cycle-slip detection generates false outputs, and range errors due to the path delay are intervened, further degrading the state estimation. This situation can be found in Figure 9 and Figure 10 (scenarios B and C, during time 872–875 s and 975 s), where TDCP-based GNSS presented a degraded position in the horizontal direction. Even worse, this degradation of state estimation could become severe in a frequent GNSS-outage environment. In Figure 12 (area I of scenario E, during time 1945–1954 s), the TDCP-based solution exhibited a deteriorated trajectory, due to the use of false information from the last epoch, where the maximum errors are 11.75 m and 5.94 m in the horizontal and vertical directions, respectively, during this period. Without additional constraint and aiding, such degraded situations are frequently presented for GNSS positioning in urban environments.

4.2. Performance Analysis of TDCP-Aided GNSS with Method A

Under an open-sky condition (e.g., PDOP value was less than four, such as Figure 11 (scenario D)), the results of TDCP-aided GNSS with/without method A were acceptable. It indicates that phase-smoothing is a critical factor in obtaining a reliable result. However, in constrained environments, due to enhanced positional connection between the adjacent epochs, TDCP-aided GNSS plus method A exhibited a more reliable outcome (blue and black lines) compared with the others without method A. For instance, in Figure 9 (scenario B, time 815 s), when the phase-smoothing approach was unavailable, the use of method A was able to improve the position state estimation, reducing the E-direction error by around 3.60 m when comparing with TDCP-based GNSS. Moreover, for the situation of deteriorated results due to the degraded cycle-slip detection (during time 872–875 s and at time 975 s), in these signal-degraded environments, the positioning errors produced with method A mitigated from 7.82 to 3.58 m (at time 873 s) and 5.53 to 2.38 m (at time 975 s) in the horizontal direction, being able to reach the SPP-nominal level. A similar case could be found in Figure 10 (scenario C, during time 963–969 s).
With the help of using method A, the benefit of the proposed zero-velocity constraint with position-fix design in (24)–(26) was significantly presented when the vehicle was stationary. Due to continuously faulty signal reception, the major contribution of method A could include the precise detection for vehicle speed, assuring the timing of stoppage, the duration time, and preventing the estimated position from presenting a deteriorated change. As shown in Figure 8, Figure 9, Figure 10 and Figure 11 (scenarios A–D, during time 550–610 s, 835–863 s, 907–954 s, and 1636–1663 s, indicated by yellow circles), the solution with method A yielded a more accurate outcome. The aid of method A maintains the position output that is the same as the position state before transiting to the static, averting the use of unhealthy measurements. Table 3 shows the related analysis for the stationary condition in the experiment. The position accuracy in this situation improved because of method A, and the position accuracy in 3D attained 4.87 m with an improvement of 41.7%, demonstrating that the zero-velocity constraint with position-fix design successfully enhanced the performance in the challenging circumstance.
For method A, the odometer provides only a relative constraint in positional displacement, and this study uses the constraint design in (16)–(22) in the velocity domain to improve this constraint, as mentioned earlier. However, this constraint is along the forward direction of the vehicle, meaning that its aid is less influential in the lateral direction, especially for the case where the estimated heading is degraded. Table 4 shows that method A improves the positional accuracy in the forward (RMSE 1.87 m), but has limited effectiveness in the lateral (RMSE 2.10 m, similar to the TDCP-based solution without method A). A representative situation was found in Figure 8 (scenario A, during time 519–547 s), where the solution with method A still presented a deteriorated trajectory misled to the oncoming road (indicated by the white dotted line). In addition, there is still a potential risk on the position-fix design for method A. In Figure 8 (scenario A, during time 550–610 s), the height error based on the scheme with method A presented an error of 12.12 m during the static motion. In this representative case, because the result of position estimation had insufficient time to revert to the stable state before transiting to the stationary, this degraded result was selected to be the position-fix output. Although, the velocity error was corrected back to the stationary state, and averted to influence the state estimation further as the vehicle started to move, this erroneous position result was still unable to adjust until the end of stationary. Moreover, even if method A was implemented, an unsatisfactory result was still observed in the previously discussed GNSS-outage environment (area I of scenario E, during time 1945–1954 s). The proposed odometer-aided FDE did not detect the false result, because the lateral errors during time 1947–1954 s did not surpass the threshold. Specifically, the effectiveness of the detection is vulnerable to the quality of the initial position, such as an inaccuracy result at time 1946 s. Although the result was relatively stable than the GNSS solutions without method A in this area (cyan and orange result), the maximum error presented an unsatisfactory outcome (9.94 m and 5.36 m in the horizontal and vertical directions) due to the limited effectiveness of method A.

4.3. Performance Analysis of TDCP-Aided GNSS with Method B

The contribution of method B is to adjust the measurement weights adaptively. In this case, the influence of inconsistent measurements can be reduced, and the estimated position becomes closer to the designated path in comparison with TDCP-based GNSS. For instance, In Figure 8 and Figure 11 (scenario A and D, during time 525–540 s and 877–884 s), the solution with method B (green and black lines) showed superior robustness during the operation in challenging environments, whereas the solutions without method B misdirected the vehicular path towards an oncoming-traffic road. In Figure 9 (scenario B, time 815 s), the positional distortion due to the multipath/NLOS contamination was improved after the intervention of method B, even when the TDCP-aided approach was inactive.
However, unsatisfactory estimation was still found in static situations. Although the performance of the TDCP-based solution with method B in the stationary condition was superior to the typical solutions without it (cyan and orange lines), the present finding exhibits that the remaining error was still up to 13.47 m in the horizontal direction (scenario B, time 853 s) and 10.62 m in the vertical direction (scenario C, time 932 s) for the harsh cases. This situation was attributed to the increase in the number of abnormal signals for the stationary condition. Another discussed case was in the GNSS-outage environment, as shown in Figure 12 (area I of scenario E, during time 1945–1954 s). In this area, the aid of method B was not obvious, where the positional error reached a maximum of up to 8.95 m and 5.43 m in the horizontal and the vertical directions. That was because method B failed to find appropriate scales for reweighting. Taking the situation of time 1948 for example (the corresponding information of measurements as shown in Table 5), it is seen that the range errors of five observations were over 10 m, even reaching 37 m. In this situation, for method B, residuals could not adequately reflect the approximation of measurement errors due to the influence of the contaminated signals, and some false measurements such as G31, G32, and C10 listed in Table 5 still significantly influenced the estimation.

4.4. Performance Analysis of TDCP-Aided GNSS with Both Methods

The proposed combination of method A and method B presents complementary attributes of both the methods and outperforms the standalone schemes of each aiding approach. For a multipath/NLOS-contaminated environment, the improvement of the absolute position estimation is mainly contributed by method B. This situation can be obviously found in Figure 8 and Figure 11 (scenario A and D), where the combined result showed that the trend of the black line was similar to that of the green line. Another contribution can also be observed under the stationary condition when method A and method B are adopted at the same time. The previous situations can be used as examples. In Figure 8 (scenario A, during time 550–610 s), the estimation of the combined methods reduced the error by 2 m in the horizontal direction and 4 m in the vertical direction, compared with those of method A only (blue line). In Figure 9, Figure 10 and Figure 11 (scenarios B, C, and D, during time 835–863 s, 907–954 s, and 1636–1663 s), this combined scheme mitigated the limitation of method B and prevented the related error growth, and the maximum of the remaining error for these areas was reduced by approximately 8.93 m, 5.97 m, and 2.91 m in the horizontal direction, and 2.71 m, 8.04 m, and 3.99 m in the vertical direction, respectively. According to the previous results, it can be found that: (1) method B provides improved position before the dynamical motion transits to stationary, and therefore, method A can provide a superior position-fix result. (2) method A effectively averts the degradation of continuous reception of abnormal signals (where method B is limited) and maintains the performance under stationary motion. According to the statistical analysis (Table 3), the TDCP-based GNSS with both methods showed an RMSE of 4.06 m in 3D with an improvement of 51.4% under the static condition. The maximum error in 3D improved from 37.24 m (conventional solution) to 8.88 m, which is significantly better than the performance based on method A only (12.45 m).
The other significant advantage of the combined method is in the frequent GNSS-outage environment. In Figure 12 (area I of scenario E, during time 1945–1954 s), the unsatisfactory performance was improved by 63.2% and 82.3% at most in the horizontal and vertical directions in contrast with method B. The main reasons are the following mechanisms from both methods in this area: (1) odometer-aided robust regression, providing a more reliable scale for reweighting and more accurate position state (compared with method B only) before driving into the GNSS-degraded environment (epochs before time 1948 s); (2) odometer-based FDE, averting unsatisfactory position state estimation due to the majority of faulty signals; (3) odometer-aided DR, bridging the duration of GNSS outage. The combined performance significantly increased the availability of satisfactory results in such a harsh situation. Although the degradation of state estimation still emerged during time 1963–1967 s, where the vehicle was performing a turning motion, and the auxiliary constraint and FDE in (19)–(23) were deactivated due to the continuous heading change, the remaining aiding methods still provided better performance compared with others, as shown in the trajectory (Figure 12).
Based on this proposed solution, the overall analysis in Table 6 presents an RMSE of 2.51 m and 3.63 m in the horizontal and vertical directions with improvements of 35.1%, and 49.6%, respectively. The indicators of the percentage of epochs, including errors of less than 3 m (narrow-road accuracy) and 5 m (wide-road accuracy), were 86.1% and 97.3% in the horizontal direction and 63.9% and 84.4% in the vertical direction for these two levels, respectively. The corresponding analysis of the empirical cumulative distribution function (CDF) is shown in Figure 13. An analysis of vehicular tracking in Table 4 shows that the proposed combination scheme reduced the forward error and lateral error with an RMSE of 1.63 m and 1.95 m, corresponding to improvements of 45.3% and 21.1%. Statistically, this finding shows that, by using a low-cost automotive-grade receiver, the proposed GNSS/odometer scheme and the combined methods without any GNSS correction services can satisfy the accuracy requirement for locating a vehicle’s position on a designated road in GNSS-challenging environments.

5. Conclusions

This study investigated GNSS/odometer integration with TDCP-designed EKF by utilizing vehicular motion constraint and robust regression for an automotive-grade receiver in an urban environment. The findings show that the phase-smoothing technique improved the state estimation and smoothed out some multipath effects after the TDCP updated the EKF state, but was subject to cycle-slip and NLOS reception. The degradation was apparent during the stationary motion. After the odometer-based aiding (method A) intervened, the aided constraints based on the resultant velocity and zero-assumption velocity of the lateral significantly enhanced the performance, when the state estimation reverted to the unsmoothed level due to a cycle slip. Moreover, the additional aiding from method A, including the zero-velocity constraint and the position-fix design, effectively averted the contamination of abnormal signals when the vehicle was stationary. However, the results show that method A had limited effectiveness in the lateral direction and presented a risk of fixing the position at an unsatisfactory level. When robust regression (method B) was added, the combination of method A and method B not only improved the accuracy in the lateral direction and the stationary condition, but also apparently enhanced the contribution of the odometer-aided FDE and DR, which effectively avoided erroneous outputs from the degradation of state estimation and performed seamless positioning, with adequate performance over a short period during a GNSS outage.
The experimental results show that the GNSS/odometer integration with the proposed methods can satisfy the demand of guiding a vehicle on designated roads with a single-frequency SPP mode in an urban area, without INS aiding. However, the performance of the proposed methods significantly depends on the empirical coefficients by the tuning process, including (21)–(23). An adequate setting for general cases is a time-consuming process. The auxiliary mechanism can include machine learning approaches, such as the random forest algorithm, for improving this situation in the future. Moreover, to mitigate the limitations of the proposed odometer-aided FDE and DR during a heading-change motion, such as a turning maneuver, the aid of z-channel gyro is considerable, and must be investigated in the future.

Author Contributions

Conceptualization, K.-W.C. and Y.-H.L.; methodology, Y.-H.L.; validation, Y.-H.L.; formal analysis, Y.-H.L., L.-T.H., and F.-Y.C.; data curation, Y.-H.L.; writing—original draft preparation, Y.-H.L.; writing—review and editing, K.-W.C., L.-T.H., and F.-Y.C.; visualization, Y.-H.L.; supervision, K.-W.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Interior under the grant number 109CCL013C.

Acknowledgments

The authors would like to acknowledge the financial support and the project under the Ministry of Interior (MOI). The authors thank the journal editor and anonymous reviewers for their constructive comments on this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, Z.; Li, B. Optimal Doppler-aided smoothing strategy for GNSS navigation. GPS Solut. 2017, 21, 197–210. [Google Scholar] [CrossRef]
  2. Kuusniemi, H.; Lachapelle, G.; Takala, J.H. Position and velocity reliability testing in degraded GPS signal environments. GPS Solut. 2017, 8, 226–237. [Google Scholar] [CrossRef]
  3. Groves, P.D. GNSS: User Equipment Processing and Errors. In Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems, 2nd ed.; Artech House: Norwood, MA, USA, 2015; pp. 401–407. [Google Scholar]
  4. Groves, P.D.; Jiang, Z.; Rudi, M.; Strode, P. A Portfolio Approach to NLOS and Multipath Mitigation in Dense Urban Areas; The Institute of Navigation: Manassas, VA, USA, 2013. [Google Scholar]
  5. Ford, T.J.; Hamilton, J. A new positioning filter: Phase smoothing in the position domain. Navigation 2003, 50, 65–78. [Google Scholar] [CrossRef]
  6. Lee, H.K.; Rizos, C.; Jee, G.I. Position domain filtering and range domain filtering for carrier-smoothed-code DGNSS: An analytical comparison. IEE Proc.-Radar Sonar Navig. 2005, 152, 271–276. [Google Scholar] [CrossRef] [Green Version]
  7. Soon, B.K.; Scheding, S.; Lee, H.K.; Lee, H.K.; Durrant-Whyte, H. An approach to aid INS using time-differenced GPS carrier phase (TDCP) measurements. GPS Solut. 2008, 12, 261–271. [Google Scholar] [CrossRef]
  8. Hsu, L.T. Analysis and modeling GPS NLOS effect in highly urbanized area. GPS Solut. 2018, 22, 7. [Google Scholar] [CrossRef] [Green Version]
  9. Braasch, M.S. Multipath Effects. In Global positioning System: Theory and Applications Volume I; Parkinson, B.W., Spilker, J.J., Jr., Eds.; AIAA: Washington, DC, USA, 1996; pp. 547–568. [Google Scholar]
  10. Leisten, O.; Knobe, V. Optimizing small antennas for body-loading applications. GPS World 2012, 23, 40–44. [Google Scholar]
  11. Blanch, J.; Walter, T.; Enge, P. Fast multiple fault exclusion with a large number of measurements. In Proceedings of the 2015 International Technical Meeting of The Institute of Navigation, Dana Point, CA, USA, 26–28 January 2015. [Google Scholar]
  12. Iwase, T.; Suzuki, N.; Watanabe, Y. Estimation and exclusion of multipath range error for robust positioning. GPS Solut. 2013, 17, 53–62. [Google Scholar] [CrossRef]
  13. Hsu, L.T.; Tokura, H.; Kubo, N.; Gu, Y.; Kamijo, S. Multiple faulty GNSS measurement exclusion based on consistency check in urban canyons. IEEE Sens. J. 2017, 17, 1909–1917. [Google Scholar] [CrossRef] [Green Version]
  14. Sato, K.; Tateshita, H.; Wakabayashi, Y. Asia Oceania Multi-GNSS Demonstration Campaign. In XXV FIG Congress; GIM International: Kuala Lumpur, Malaysia, 2014. [Google Scholar]
  15. Gao, H.; Groves, P.D. Environmental context detection for adaptive navigation using GNSS measurements from a smartphone. Navig. J. Inst. Navig. 2018, 65, 99–116. [Google Scholar] [CrossRef]
  16. Castaldo, G.; Angrisano, A.; Gaglione, S.; Troisi, S. P-RANSAC: An integrity monitoring approach for GNSS signal degraded scenario. Int. J. Navig. Obs. 2014, 2014, 1–11. [Google Scholar] [CrossRef] [Green Version]
  17. Wieser, A.; Brunner, F.K. Short static GPS sessions: Robust estimation results. GPS Solut. 2002, 5, 70–79. [Google Scholar] [CrossRef]
  18. Gaglione, S.; Innac, A.; Carbone, S.P.; Troisi, S.; Angrisano, A. Robust estimation methods applied to GPS in harsh environments. In Proceedings of the 2017 IEEE European Navigation Conference (ENC), Lausanne, Switzerland, 9–12 May 2017; pp. 14–25. [Google Scholar]
  19. Gaglione, S.; Angrisano, A.; Crocetto, N. Robust Kalman Filter applied to GNSS positioning in harsh environment. In Proceedings of the 2019 IEEE European Navigation Conference (ENC), Warsaw, Poland, 9–12 April 2019; pp. 1–6. [Google Scholar]
  20. Groves, P.D.; Jiang, Z.; Wang, L.; Ziebart, M.K. Intelligent urban positioning using multi-constellation GNSS with 3D mapping and NLOS signal detection. In Proceedings of the ION ITM, Nashville, TN, USA, 17–21 September 2012; pp. 458–472. [Google Scholar]
  21. Ziedan, N.I. Urban positioning accuracy enhancement utilizing 3D buildings model and accelerated ray tracing algorithm. In Proceedings of the 30th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS+), Portland, OR, USA, 25–29 September 2017. [Google Scholar]
  22. Hsu, L.T.; Gu, Y.; Kamijo, S. NLOS correction/exclusion for GNSS measurement using RAIM and city building models. Sensors 2015, 15, 17329–17349. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wang, L.; Groves, P.D.; Ziebart, M.K. Smartphone shadow matching for better cross-street GNSS positioning in urban environments. J. Navig. 2015, 68, 411–433. [Google Scholar] [CrossRef] [Green Version]
  24. Peyraud, S.; Bétaille, D.; Renault, S.; Ortiz, M.; Mougel, F.; Meizel, D.; Peyret, F. About non-line-of-sight satellite detection and exclusion in a 3D map-aided localization algorithm. Sensors 2013, 13, 829–847. [Google Scholar] [CrossRef] [Green Version]
  25. Ng, H.-F.; Zhang, G.; Hsu, L.-T. A Computation Effective Range-based 3D Mapping Aided GNSS using Skymask. J. Navig. 2020, 1, 1–21. [Google Scholar] [CrossRef]
  26. Wen, W.; Zhang, G.; Hsu, L.T. Correcting NLOS by 3D LiDAR and Building Height to improve GNSS Single Point Positioning. Navigation 2019, 66, 705–718. [Google Scholar] [CrossRef]
  27. Gao, Z.; Shen, W.; Zhang, H.; Niu, X.; Ge, M. Real-time kinematic positioning of INS tightly aided multi-GNSS ionospheric constrained PPP. Sci. Rep. 2016, 6, 1–16. [Google Scholar] [CrossRef] [Green Version]
  28. El-Sheimy, N. The Development of VISAT: A Mobile Survey System for GIS Applications; University of Calgary: Calgary, AB, Canada, 2016. [Google Scholar]
  29. Shin, E.H. Accuarcy Improvement of Low Cost INS/GPS for Land Applications. Graduate Studies; University of Calgary: Calgary, AB, Canada, 2001. [Google Scholar]
  30. Sharaf, R.; Noureldin, A.; Osman, A.; El-Sheimy, N. Online INS/GPS integration with a radial basis function neural network. IEEE Aerosp. Electron. Syst. Mag. 2005, 20, 8–14. [Google Scholar] [CrossRef]
  31. El-Sheimy, N.; Youssef, A. Inertial sensors technologies for navigation applications: State of the art and future trends. Satell. Navig. 2020, 1, 2. [Google Scholar] [CrossRef] [Green Version]
  32. Shin, E.H. Estimation Techniques for Low-Cost Inertial Navigation. Ph.D. Thesis, University of Calgary, Calgary, Alberta, Canada, May 2005. [Google Scholar]
  33. El-Sheimy, N.; Hou, H.; Niu, X. Analysis and modeling of inertial sensors using Allan variance. IEEE Trans. Instrum. Meas. 2007, 57, 140–149. [Google Scholar] [CrossRef]
  34. El-Sheimy, N.; Chiang, K.W.; Noureldin, A. The utilization of artificial neural networks for multisensor system integration in navigation and positioning instruments. IEEE Trans. Instrum. Meas. 2006, 55, 1606–1615. [Google Scholar] [CrossRef]
  35. Lahrech, A.; Boucher, C.; Noyer, J.C. Fusion of GPS and odometer measurements for map-based vehicle navigation. In Proceedings of the 2004 IEEE International Conference on Industrial Technology (ICIT’04), Hammamet, Tunisia, 8–10 December 2004; Volume 2, pp. 944–948. [Google Scholar]
  36. Seo, J.; Lee, H.K.; Lee, J.G.; Park, C.G. Lever arm compensation for GPS/INS/odometer integrated system. Int. J. Control Autom. Syst. 2006, 4, 247–254. [Google Scholar]
  37. Mosavi, M.R.; Azad, M.S.; EmamGholipour, I. Position estimation in single-frequency GPS receivers using Kalman filter with pseudo-range and carrier phase measurements. Wirel. Pers. Commun. 2013, 72, 2563–2576. [Google Scholar] [CrossRef]
  38. Misra, P.; Enge, P. Global Positioning System: Signals, Measurements and Performance, 2nd ed.; Ganga-Jamuna Press: Lincoln, MA, USA, 2006; Volume 206. [Google Scholar]
  39. van Graas, F.; Soloviev, A. Precise Velocity Estimation Using a Stand-Alone GPS Receiver. Navigation 2004, 51, 283–292. [Google Scholar] [CrossRef]
  40. Serrano, L.; Kim, D.; Langley, R.B.; Itani, K.; Ueno, M. A GPS velocity sensor: How accurate can it be?—A first look. In Proceedings of the ION NTM, San Diego, CA, USA, 26–28 January 2004; pp. 875–885. [Google Scholar]
  41. Wendel, J.; Meister, O.; Monikes, R.; Trommer, G.F. Time-differenced carrier phase measurements for tightly coupled GPS/INS integration. In Proceedings of the 2006 IEEE/ION Position, Location, and Navigation Symposium, San Diego, CA, USA, 25–27 April 2006; pp. 54–60. [Google Scholar]
  42. Ding, W.; Wang, J. Precise velocity estimation with a stand-alone GPS receiver. J. Navig. 2001, 64, 311–325. [Google Scholar] [CrossRef] [Green Version]
  43. Freda, P.; Angrisano, A.; Gaglione, S.; Troisi, S. Time-differenced carrier phases technique for precise GNSS velocity estimation. GPS Solut. 2015, 19, 335–341. [Google Scholar] [CrossRef]
  44. Aftatah, M.; Lahrech, A.; Abounada, A. Fusion of GPS/INS/Odometer measurements for land vehicle navigation with GPS outage. In Proceedings of the 2016 IEEE 2nd International Conference on Cloud Computing Technologies and Applications (CloudTech), Marrakech, Morocco, 24–26 May 2016; pp. 48–55. [Google Scholar]
  45. Georgy, J.; Noureldin, A.; Bayoumi, M. Mixture particle filter for low cost INS/odometer/GPS integration in land vehicles. In Proceedings of the VTC Spring 2009-IEEE 69th Vehicular Technology Conference, Barcelona, Spain, 26–29 April 2009; pp. 1–5. [Google Scholar]
  46. Knight, N.L.; Wang, J. A comparison of outlier detection procedures and robust estimation methods in GPS positioning. J. Navig. 2009, 62, 699–709. [Google Scholar] [CrossRef]
  47. Stephenson, S. Automotive Applications of High Precision GNSS. Ph.D. Thesis, University of Nottingham, Nottingham, UK, June 2016. [Google Scholar]
  48. Basnayake, C.; Williams, T.; Alves, P.; Lachapelle, G. Can GNSS Drive V2X? GPS World 2010, 21, 35–43. [Google Scholar]
  49. National Highway Traffic Safety Administration and Department of Transportation. Federal Motor Vehicle Safety Standards; V2V Communications, Docket No. NHTSA–2016–0126, RIN 2127–AL55; Federal Register: Washington, DC, USA, 2017.
  50. Reid, T.G.; Houts, S.E.; Cammarata, R.; Mills, G.; Agarwal, S.; Vora, A.; Pandey, G. Localization requirements for autonomous vehicles. arXiv 2019, arXiv:1906.01061. [Google Scholar] [CrossRef] [Green Version]
Figure 1. GNSS/Odometer integration scheme with WLS-based robust regression.
Figure 1. GNSS/Odometer integration scheme with WLS-based robust regression.
Remotesensing 12 02550 g001
Figure 2. Relationship between the vehicle frame and the global frame in the horizontal plane.
Figure 2. Relationship between the vehicle frame and the global frame in the horizontal plane.
Remotesensing 12 02550 g002
Figure 3. Process flow of the designed scheme with the proposed methods.
Figure 3. Process flow of the designed scheme with the proposed methods.
Remotesensing 12 02550 g003
Figure 4. (a) NovAtel SPAN® LCI and GPS-GGG-703-HV antenna; (b) Allystar EVK-2040 and AGR 6301 antenna; (c) Experimental platform.
Figure 4. (a) NovAtel SPAN® LCI and GPS-GGG-703-HV antenna; (b) Allystar EVK-2040 and AGR 6301 antenna; (c) Experimental platform.
Remotesensing 12 02550 g004
Figure 5. Experimental trajectory.
Figure 5. Experimental trajectory.
Remotesensing 12 02550 g005
Figure 6. (ad) Street views of the experimental route for the evaluated scenarios A–D, (e) area I of scenario E, and (f) other areas of scenario E.
Figure 6. (ad) Street views of the experimental route for the evaluated scenarios A–D, (e) area I of scenario E, and (f) other areas of scenario E.
Remotesensing 12 02550 g006
Figure 7. (a) PDOP and (b) number of satellites.
Figure 7. (a) PDOP and (b) number of satellites.
Remotesensing 12 02550 g007aRemotesensing 12 02550 g007b
Figure 8. (a) Evaluated trajectory and (b) error analysis for scenario A.
Figure 8. (a) Evaluated trajectory and (b) error analysis for scenario A.
Remotesensing 12 02550 g008aRemotesensing 12 02550 g008b
Figure 9. (a) Evaluated trajectory and (b) error analysis for scenario B.
Figure 9. (a) Evaluated trajectory and (b) error analysis for scenario B.
Remotesensing 12 02550 g009aRemotesensing 12 02550 g009b
Figure 10. (a) Evaluated trajectory and (b) error analysis for scenario C.
Figure 10. (a) Evaluated trajectory and (b) error analysis for scenario C.
Remotesensing 12 02550 g010aRemotesensing 12 02550 g010b
Figure 11. (a) Evaluated trajectory and (b) error analysis for scenario D.
Figure 11. (a) Evaluated trajectory and (b) error analysis for scenario D.
Remotesensing 12 02550 g011
Figure 12. (a) Evaluated trajectory and (b) error analysis for scenario E.
Figure 12. (a) Evaluated trajectory and (b) error analysis for scenario E.
Remotesensing 12 02550 g012
Figure 13. CDF for the different GNSS solutions; (a) horizontal position error; and (b) vertical position error. (The scope of the x-axis is limited within 15 m for better viewing.)
Figure 13. CDF for the different GNSS solutions; (a) horizontal position error; and (b) vertical position error. (The scope of the x-axis is limited within 15 m for better viewing.)
Remotesensing 12 02550 g013
Table 1. Specifications of NovAtel SPAN® LCI.
Table 1. Specifications of NovAtel SPAN® LCI.
AccelerometerGyroscope
Bias instability 0.1   mg 0.05 ° / hr
Random walk noise 0.0588   m / s / Hz 0.012 ° / hr
Table 2. Performance of OBD II Odometer.
Table 2. Performance of OBD II Odometer.
Unit: m/sOverall 1Straight 2Turning 3
Mean0.0530.0680.116
Max1.2820.9291.043
STD0.1410.1320.205
RMSE0.1510.1480.235
1 “Overall” indicates all of the data for evaluation. 2 “Straight” indicates the evaluated data without turning and stationary motion. 3 “Turning” indicates the evaluated data based on turning motion only.
Table 3. Performance analysis for the stationary condition.
Table 3. Performance analysis for the stationary condition.
Unit: m3D
Conventional GNSSTDCP-Aided GNSSTDCP-Aided GNSS +
Method A
TDCP-Aided GNSS +
Method B
TDCP-Aided GNSS + Method A Method B
Mean5.844.754.164.803.65
Max37.2419.4212.4515.008.88
RMSE8.355.684.875.574.06
RMSE
Improvement
-32.0%41.7%33.3%51.4%
Table 4. Performance analysis for GNSS solution in vehicular path tracking.
Table 4. Performance analysis for GNSS solution in vehicular path tracking.
Unit: mForward/Lateral
Conventional GNSSTDCP-Aided GNSSTDCP-Aided GNSS +
Method A
TDCP-Aided GNSS +
Method B
TDCP-Aided GNSS + Method A Method B
Mean−0.41/−0.470.39/−0.54−0.49/−0.61−0.47/−0.55−0.37/−0.49
Max25.95/17.2114.92/10.0512.6/9.3015.86/10.6010.75/12.75
RMSE2.98/2.472.21/2.071.87/2.102.09/2.161.63/1.95
RMSE
Improvement
-25.8%/16.2%37.2%/15.0%29.9%/12.6%45.3%/21.1%
Table 5. Code measurement information in WLS-based filter at time 1948.
Table 5. Code measurement information in WLS-based filter at time 1948.
Satellite 1Range Measurement 2True RangeRange Error 3Residual 4Scale 4
G1422,492,867.5622,492,830.4037.1626.962.77
G2520,580,188.5620,580,187.610.953.781.00
G3121,729,429.2521,729,415.4913.76−1.261.00
G3221,213,735.1221,213,711.4823.6518.871.94
C635,967,621.8835,967,615.376.515.391.00
C736,142,589.4936,142,586.113.38−3.681.00
C936,559,641.7536,559,610.2831.4724.932.56
C1037,151,565.4737,151,546.0019.477.931.00
1 “G” and “C” indicate GPS and Beidou satellite system, respectively. 2 Range measurements were corrected by the process of broadcast ephemeris. Moreover, considering that WLS residuals are derived from the estimated states (including the receiver clock bias), range measurements in this table are also adjusted by the estimated receiver bias for consistency. 3 Range error means the difference between range measurements and true ranges. 4 The definitions correspond to (27).
Table 6. Performance analysis in position domain.
Table 6. Performance analysis in position domain.
Unit: mHorizontal/Vertical
Conventional GNSSTDCP-Aided GNSSTDCP-Aided GNSS +
Method A
TDCP-Aided GNSS +
Method B
TDCP-Aided GNSS + Method A Method B
Mean2.83/4.782.47/3.412.50/3.172.57/3.052.28/2.88
Max27.66/34.1815.20/14.9513.15/17.4716.82/14.9714.61/12.67
RMSE3.87/7.203.01/4.552.79/4.162.98/4.112.51/3.63
RMSE
Improvement
-22.3%/36.8%28.0%/42.0%23.0%/42.9%35.1%/49.6%
Which Road1(3 m)272.0%/42.4%75.4%/49.3%81.8%/60.4%76.7%/58.1%86.1%/63.9%
Which Road1(5 m)289.5%/65.8%93.5%/75.4%96.3%/82.9%94.4%/79.5%97.3%/84.4%
1 The accuracy category refers to [47,48,49,50]. 2 Percentage of epochs that had errors of less than 3 m or 5 m.

Share and Cite

MDPI and ACS Style

Chiang, K.-W.; Li, Y.-H.; Hsu, L.-T.; Chu, F.-Y. The Design a TDCP-Smoothed GNSS/Odometer Integration Scheme with Vehicular-Motion Constraint and Robust Regression. Remote Sens. 2020, 12, 2550. https://doi.org/10.3390/rs12162550

AMA Style

Chiang K-W, Li Y-H, Hsu L-T, Chu F-Y. The Design a TDCP-Smoothed GNSS/Odometer Integration Scheme with Vehicular-Motion Constraint and Robust Regression. Remote Sensing. 2020; 12(16):2550. https://doi.org/10.3390/rs12162550

Chicago/Turabian Style

Chiang, Kai-Wei, Yu-Hua Li, Li-Ta Hsu, and Feng-Yu Chu. 2020. "The Design a TDCP-Smoothed GNSS/Odometer Integration Scheme with Vehicular-Motion Constraint and Robust Regression" Remote Sensing 12, no. 16: 2550. https://doi.org/10.3390/rs12162550

APA Style

Chiang, K. -W., Li, Y. -H., Hsu, L. -T., & Chu, F. -Y. (2020). The Design a TDCP-Smoothed GNSS/Odometer Integration Scheme with Vehicular-Motion Constraint and Robust Regression. Remote Sensing, 12(16), 2550. https://doi.org/10.3390/rs12162550

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop