# Performance Analysis of Two-Dimensional Dead Reckoning Based on Vehicle Dynamic Sensors during GNSS Outages

Authors: Joong-hee Han; Chi-ho Park; Chang-Ki Hong; Jay Hyoun Kwon
Journal: Journal of Sensors
 (2017)

---

## Abstract

Recently, to improve safety and convenience in driving, numerous sensors are mounted on cars to operate advanced driver assistant systems. Among various sensors, vehicle dynamic sensors can measure the vehicle motions such as speed and rotational angular speed for dead reckoning, which can be applied to develop a land vehicle positioning system to overcome the weaknesses of the GNSS technique. In this paper, three land vehicle positioning algorithms that integrate GNSS with vehicle dynamic sensors including a wheel speed sensor (WSS), a yaw rate sensor (YRS), and a steering angle sensor (SAS) are implemented, and then a performance evaluation was conducted during GNSS outages. Using a loosely coupled strategy, three integration algorithms are designed, namely, GNSS/WSS, GNSS/WSS/YRS, and GNSS/WSS/YRS/SAS. The performance of the three types of integration algorithm is evaluated based on two data sets. The results indicate that both the GNSS/WSS/YRS integration and the GNSS/WSS/YRS/SAS integration could estimate the horizontal position with meter-level accuracy during 30-second GNSS outages. However, the GNSS/WSS integration would provide an unstable navigation solution during GNSS outages due to the accuracy limitation of the computed yaw rate using WSS.


---

## Body

## 1. Introduction

The vehicle positioning technique is a key component in car navigation to find and guide routes. Typical car navigation uses a low-cost GNSS (Global Navigation Satellite Systems) receiver that can provide the position and velocity of a vehicle with an accuracy that is appropriate for navigation. However, the major drawback of the GNSS based positioning technique is that the performance depends on the satellite signal reception environment. In particular, cars often move near buildings and tunnels where the GNSS signal environment is poor, and thus the GNSS based positioning technique cannot guarantee continuity and reliability of positioning [1, 2]. Therefore, in order to overcome this problem, the GNSS based technique is combined with a dead reckoning (DR) sensor, which does not rely on any external sources. A well-known DR sensor is the inertial navigation system (INS), which estimates the navigation solution by using the acceleration and angular rate measured by the inertial measurement unit (IMU) [1]. In particular, advances in MEMS (microelectromechanical system) technology have led to the development of a small-sized and low-cost INS [3], which is integrated with GNSS for a continuous and accurate solution in car navigation systems [3–5]. Although GNSS/MEMS-INS can fulfill the required performance for car navigation system, these sensors are too expensive for use in commercial navigation systems [6].Recently, several ADAS (Advanced Driver Assistance Systems) such as ABS (Antilock Brake System), ESC (Electronic Stability Control), and ACC (Adaptive Cruise Control) have been applied to most passenger cars for driver and passenger safety [7]. Since ADAS operation relies on information about the vehicle environment and dynamics, numerous vehicle dynamic sensors and surround sensors and actuators are built in a car [8]. Among the sensors, vehicle dynamic sensors such as wheel speed sensor (WSS), steering angle sensor (SAS), and yaw rate sensor (YRS) measure and provide the vehicle motions for two-dimensional DR navigation. Therefore, as vehicle dynamic sensors can be applied to overcome the limitations of the GNSS technique, various studies combining GNSS and vehicle dynamic sensors have been conducted. Bonnifait et al. [9] developed a localization system using GPS and WSS. They showed that a two-dimensional DR navigation based on WSS could estimate the horizontal positions with meter-level accuracy with GPS. However, this study has not been fully analyzed to evaluate the performance of WSS based DR navigation during GNSS outages. Iqbal [2] and Georgy et al. [10] studied the performance of two-dimensional DR navigation using the yaw rate measured by a gyro as well as velocity based on WSS in GNSS signal blockage situations. However, since these studies used the yaw rate measured from an IMU, they are limited with respect to verifying the performance of DR navigation based on vehicle dynamic sensors mounted on a car. The integration DR navigation based on a dynamic vehicle model using vehicle dynamic sensors with GNSS was studied [11, 12]. The dynamic model takes into account the vehicle dynamic characteristic such as lateral force and tire side slip angle that could estimate stable positioning results rather than the vehicle model based on the geometric relationship. However, to implement the dynamic vehicle model, the vehicle parameters such as front and rear slip angle, tire cornering stiffness, and angular speed of the wheels are required [3], and thus it is difficult to apply the above method to actual vehicle positioning.In this study, in order to overcome the limitations of the GNSS based positioning technique, three land vehicle positioning algorithms that integrate GNSS with vehicle dynamic sensors (WSS, YRS, and SAS) are implemented, and the performance evaluation was conducted in GNSS signal blockage situations. A description of the integration strategies as well as mathematical models of three types of the GNSS/vehicle dynamic sensors is presented in Section2. A performance evaluation of the three integration algorithms is provided in Section 3. The conclusion are given in Section 4.


## 2. GNSS/Vehicle Dynamic Sensor Integration

### 2.1. GNSS/Vehicle Dynamic Sensor Integration Strategies

DR is the process of estimating the current position based on the previous position using the velocity and the traveling direction measured by DR sensors. To calculate the horizontal position of the vehicle based on DR, the vehicle’s velocity and yaw information are required. To apply vehicle dynamic sensors for the two-dimensional DR navigation, the relationship of the vehicle dynamic sensors and navigation information is summarized in the Table1.Table 1
Relationship of the vehicle dynamic sensors and navigation information for DR.
Sensor name
Velocity
Yaw rate
Others
Wheel speed sensor (WSS)
Directly
Derived from the differential of left and right wheel speed
—
Steering angle sensor (SAS)
—
Derived from kinematic model
Can estimate the side slip by kinematic relationship
Yaw rate sensor (YRS)
DirectlyTo implement a two-dimensional DR navigation based on vehicle dynamic sensors, three different combination strategies are designed to ensure that the vehicle dynamic sensors’ information is not duplicated. The strategies for two-dimensional DR navigation based on vehicle dynamic sensors are summarized in Table2.Table 2
Description of combination strategies of vehicle dynamic sensors for DR.
Combination strategies
Description
WSS
(i) Nonholonomic constraints are applied in lateral direction(ii) WSS provides both the longitudinal velocity and yaw rate
WSS/YRS
(i) Nonholonomic constraints are applied in lateral direction(ii) WSS and YRS provides the longitudinal velocity and yaw rate, respectively
WSS/YRS/SAS
(i) The side slip angle is computed by SAS and side slip ratio(ii) YRS provides the yaw rate (iii) Speed from WSS and side slip angle compute longitudinal and lateral velocityIn this study, three integration algorithms are proposed by integrating GNSS and a two-dimensional DR navigation based on combination strategies of vehicle dynamic sensors using a loosely coupled mode. In addition, three integration algorithms are implemented through the extended Kalman filter. The output of integration algorithm rate is set to 50 Hz. The GNSS measurements are the position and the velocity calculated from the C/A code and the Doppler measurements at a GNSS receiver. The GNSS measurement update rate is set to 1 Hz.


### 2.2. GNSS/WSS Integration Algorithm

A block diagram of the GNSS/WSS integration algorithm is shown in Figure1. A two-dimensional DR based on WSS uses the speed measured at a pair of WSS of the rear wheels, and the procedure for this is as follows. When both the rear left and right WSS measure the wheel speed, the vehicle’s speed, vr, is calculated at the center of the rear wheel axle as follows [13]:(1)vr=vrL+vrR2,where vrL and vrR are the rear left and rear right WSS measurements, respectively.Figure 1
Block diagram of GNSS/WSS integration algorithm.To define the velocity in the body frame using the vehicle’s speed, we assume that the center of the body frame is set to the center of the rear wheel axle, the vehicle drives on a flat road, and no wheel slip occurs. In addition, assuming that the direction of the vehicle’s speed at the rear axle is equal to the longitudinal axis of the body frame, the velocities in the body frame are defined as follows:(2)vxvy=vr0,where vx and vy are the longitudinal and lateral velocity in the body frame, respectively.When a vehicle turns, the left and right wheel speeds are different. The speed measured by the individual WSS varies from along track while the vehicle is turning [1]. Therefore, the yaw rate could be calculated as the difference between the rear left wheel speed and the rear right wheel speed divided by the length between the rear wheels [13], as expressed in(3)ψ˙=vrL-vrRTr,where ψ˙ is the yaw rate of vehicle and Tr is the distance between the rear wheels.The yaw angle of the vehicle could be computed by numerical integration using the yaw at the previous time and the yaw rate, which is written as follows:(4)ψt=ψt-1+ψ˙Δt,where ψt and ψt-1 denote the yaw at time t and t-1, respectively, and Δt is the sampling time interval.The velocities in the navigation frame could be transformed by using(5)vnve=cos⁡ψ-sin⁡ψsin⁡ψcos⁡ψvxvy,where vn and ve denote the north velocity and the east velocity in the navigation frame, respectively, ψ is the yaw angle at the current time, and vx and vy are the longitudinal and lateral velocity in the body frame, respectively.In EKF-based GNSS/WSS integration, the state vector of the navigation error is composed of the latitude and longitude error, the north and east velocity error, and the yaw error. The state vector of the sensor error is composed of the average and the difference of the left and right WSS scale factor, defined as random constants. The white noise vector includes the white noise of WSS derived vehicle speed and the white noise of WSS derived yaw rate. The dynamic model for the GNSS/WSS integration is given in(6)δφ˙δλ˙δv˙nδv˙eδψ˙δs˙vδΔs˙v︸=δx˙WSS=001M+h00000vetan⁡φN+h0cos⁡φ001N+h0cos⁡φ0000000-vxsin⁡ψcos⁡ψ00000vxcos⁡ψsin⁡ψ000000ψ˙vxTr00000000000000︸=FWSSδφδλδvnδveδψδsvδΔsv︸=δxWSS+cos⁡ψM+h00sin⁡ψN+h0cos⁡φ00000010000︸=GWSSwvWSSwψ˙WSS︸=uWSS,where FWSS, δxWSS, GWSS, and uWSS are the dynamic matrix, the error state vector, the shaping matrix, and the white noise vector for GNSS/WSS integration, respectively, δφ and δλ are the latitude and longitude errors, respectively, δvn and δve are the north and east velocity errors, respectively, δψ is the yaw error, δsv and δΔsv are the average and the difference of the left and right WSS scale factors, respectively, M is the radius of curvature of the meridian, N is the radius of curvature of the prime vertical, h0 is the ellipsoidal height defined as a constant, vx and ve are the longitudinal and the east velocity, respectively, φ is the latitude, ψ is the yaw, Tr is the length between rear wheels, wvWSS is the white noise of WSS derived vehicle speed, and wψ˙WSS is the white noise of WSS derived yaw rate.The measurement model is generally expressed as follows:(7)z=Hδx+v,v~0,R,where z is the measurement vector, H is the design matrix, δx is the error state vector, v is the measurement error, and R is the variance-covariance matrix of the measurement error.In this study, the GNSS measurement vector consists of difference between the latitude, the longitude, the north velocity, and the east velocity estimated from DR based on vehicle dynamic sensors and GNSS receiver, as shown below:(8)zGNSS=φλvnveDR-φλvnveGNSS,where zGNSS is the GNSS measurement vector, the subscripts DR and GNSS denote the values computed by DR based on vehicle dynamic sensors and the values calculated from the GNSS receiver, respectively, φ and λ are the latitude and the longitude, respectively, and vn and ve are the north velocity and the east velocity, respectively.The design matrix and the variance-covariance matrix for GNSS measurements are shown, respectively, as follows:(9)HGNSS=I4×404×3,(10)RGNSS=σφ20000σλ20000σvn20000σve2,where HGNSS is the GNSS design matrix, RGNSS is the GNSS variance-covariance matrix, I4×4 is the 4×4 identity matrix, 04×3 is the 4×3 zero matrix, and σφ2, σλ2, σvn2, and σve2 are the variance of the estimated latitude, the longitude, the north velocity, and the east velocity from GNSS receiver, respectively.


### 2.3. GNSS/WSS/YRS Integration Algorithm

Figure2 shows a block diagram of the GNSS/WSS/YRS integration algorithm. Its structure is similar to that of the GNSS/WSS integration algorithm shown in Figure 1. The major difference between the DR mechanization based on WSS and the DR mechanization based on WSS/YRS is that the yaw rate value is measured from YRS. Therefore, in EKF-based GNSS/WSS/YRS integration, the state vector of the navigation error is the same as that in GNSS/YRS integration, and the state vector of the sensor error is composed of the average of the left and right WSS scale factor and the YRS bias, defined as random constants. The dynamic model for the GNSS/WSS/YRS integration is given in(11)δφ˙δλ˙δv˙nδv˙eδψ˙δs˙vδb˙YRS︸=δx˙WSS/YRS=001M+h00000vetan⁡φN+h0cos⁡φ001N+h0cos⁡φ0000000-vxsin⁡ψcos⁡ψ00000vxcos⁡ψsin⁡ψ0000000100000000000000︸=FWSS/YRSδφδλδvnδveδψδsvδbYRS︸=δxWSS/YRS+cos⁡ψM+h00sin⁡ψN+h0cos⁡φ00000010000︸=GWSS/YRSwvWSSwψ˙YRS︸=uWSS/YRS,where FWSS/YRS, δxWSS/YRS, GWSS/YRS, and uWSS/YRS are the dynamic matrix, the error state vector, the shaping matrix, and the white noise vector for GNSS/WSS/YRS integration, respectively,δbYRS is the YRS bias, and wψ˙YRS is the white noise of YRS.Figure 2
Block diagram of GNSS/WSS/YRS integration algorithm.Comparing the GNSS/WSS integration, the GNSS/WSS/YRS integration performs not only the GNSS measurement update but also the ZIHR (Zero Integrated Heading Rate) measurement update [3] to prevent the drift error of yaw when the vehicle is stationary. To update the error state vector by using GNSS based position and velocity, the GNSS measurement model is the same as that of GNSS/WSS integration, as shown in (8)–(10). When the vehicle stops, the measured speed by WSS is zero, and thus ZIHR can be applied to fix the yaw. The measurement vector for ZIHR is given by(12)zZIHR=ψk-ψk-1,where zZIHR is the ZIHR measurement vector, ψk is the estimated yaw at the current time, and ψk-1 is the estimated yaw at the previous time.The design matrix and the variance-covariance matrix for ZIHR are shown, respectively, as follows:(13)HZIHR=0000100,(14)RZIHR=σψ2,where HZIHR is the ZIHR design matrix, RZIHR is the ZIHR variance-covariance matrix, and σψ2 is the variance of the estimated yaw.


### 2.4. GNSS/WSS/YRS/SAS Integration Algorithm

Figure3 shows a block diagram of the GNSS/WSS/YRS/SAS integration algorithm. The idea of integrating SAS for DR is to compute the side slip angle by using the handle steering angle from SAS and then to employ the estimated side slip angle to set the longitudinal and lateral velocity from the WSS derived vehicle speed. In the case of the integration algorithm in Sections 2.2 and 2.3, a lateral constraint is applied to set the velocity in the body frame using WSS. However, in a land vehicle positioning system, since the lateral nonholonomic constraint is violated by a side slip during either cornering or change in road and tire conditions [14], the positioning accuracy would be degraded. Although the side slip angle is a very complicated phenomenon associated with road and tire conditions as well as high vehicle dynamics including fast driving and sharp turns [15], side slip always occurs during turning due to lateral tire deformation [16]. In order to define the lateral velocity by side slip in the two-dimensional DR based on WSS/YRS/SAS, the side slip angle model is defined in (15) with the assumption of a linear relationship between the side slip angle and the handle steering angle:(15)β=γ·α,where β is the side slip angle, γ is the side slip ratio, and α is the handle steering angle measured by SAS.Figure 3
Block diagram of GNSS/WSS/YRS/SAS integration algorithm.Considering the side slip angle, the velocities in the body frame are defined as follows:(16)vxvy=vrcos⁡βvrsin⁡β.In the two-dimensional DR based on WSS/YRS/SAS, the velocities in the navigation frame are given by (17)vnve=cos⁡ψ-sin⁡ψsin⁡ψcos⁡ψvxvy=vrcos⁡ψ+βvrsin⁡ψ+β.In EKF-based GNSS/WSS/YRS/SAS integration, the state vector of the navigation error is composed of the latitude and longitude error, the north and east velocity error, the yaw error, and the side slip angle error. The state vector of the sensor error consists of the average of the left and right WSS scale factor, the YRS bias, the SAS scale factor, and the SAS bias, defined as random constants. The dynamic model for the GNSS/WSS/YRS/SAS integration is expressed as follows:(18)δφ˙δλ˙δv˙nδv˙eδψ˙δβ˙δs˙vδb˙YRSδs˙SASδb˙SAS︸=δx˙WSS/YRS/SAS=001M+h00000000vetan⁡φN+h0cos⁡φ001N+h0cos⁡φ0000000000-vxsin⁡ψ-vycos⁡ψ-vxsin⁡ψ-vycos⁡ψcos⁡ψ+β0000000vxcos⁡ψ-vysin⁡ψvxcos⁡ψ-vysin⁡ψsin⁡ψ+β000000000010000000000β10000000000000000000000000000000000000000︸=FWSS/YRS/SASδφδλδvnδveδψδβδsvδbYRSδsSASδbSAS︸=δxWSS/YRS/SAS+cos⁡ψM+h000sin⁡ψN+h0cos⁡φ00000000010001000000000000︸=GWSS/YRS/SASwvWSSwψ˙YRSwβ˙SAS︸=uWSS/YRS/SAS,where FWSS/YRS/SAS, δxWSS/YRS/SAS, GWSS/YRS/SAS, and uWSS/YRS/SAS are the dynamic matrix, the error state vector, the shaping matrix, and the white noise vector for GNSS/WSS/YRS/SAS integration, respectively, δβ is the side slip angle error, δsSAS is the SAS scale factor, δbSAS is the SAS bias, and wβ˙YRS is the white noise of the side slip angle derived SAS.The measurement vector and the variance-covariance matrix for both the GNSS measurement model and the ZIHR measurement model are the same as that of GNSS/WSS/YRS integration. The design matrices for measurement models in the GNSS/WSS/YRS/SAS integration are expressed as follows:(19)HGNSS=I4×404×6,HZIHR=000010000000.


## 2.1. GNSS/Vehicle Dynamic Sensor Integration Strategies

DR is the process of estimating the current position based on the previous position using the velocity and the traveling direction measured by DR sensors. To calculate the horizontal position of the vehicle based on DR, the vehicle’s velocity and yaw information are required. To apply vehicle dynamic sensors for the two-dimensional DR navigation, the relationship of the vehicle dynamic sensors and navigation information is summarized in the Table1.Table 1
Relationship of the vehicle dynamic sensors and navigation information for DR.
Sensor name
Velocity
Yaw rate
Others
Wheel speed sensor (WSS)
Directly
Derived from the differential of left and right wheel speed
—
Steering angle sensor (SAS)
—
Derived from kinematic model
Can estimate the side slip by kinematic relationship
Yaw rate sensor (YRS)
DirectlyTo implement a two-dimensional DR navigation based on vehicle dynamic sensors, three different combination strategies are designed to ensure that the vehicle dynamic sensors’ information is not duplicated. The strategies for two-dimensional DR navigation based on vehicle dynamic sensors are summarized in Table2.Table 2
Description of combination strategies of vehicle dynamic sensors for DR.
Combination strategies
Description
WSS
(i) Nonholonomic constraints are applied in lateral direction(ii) WSS provides both the longitudinal velocity and yaw rate
WSS/YRS
(i) Nonholonomic constraints are applied in lateral direction(ii) WSS and YRS provides the longitudinal velocity and yaw rate, respectively
WSS/YRS/SAS
(i) The side slip angle is computed by SAS and side slip ratio(ii) YRS provides the yaw rate (iii) Speed from WSS and side slip angle compute longitudinal and lateral velocityIn this study, three integration algorithms are proposed by integrating GNSS and a two-dimensional DR navigation based on combination strategies of vehicle dynamic sensors using a loosely coupled mode. In addition, three integration algorithms are implemented through the extended Kalman filter. The output of integration algorithm rate is set to 50 Hz. The GNSS measurements are the position and the velocity calculated from the C/A code and the Doppler measurements at a GNSS receiver. The GNSS measurement update rate is set to 1 Hz.


## 2.2. GNSS/WSS Integration Algorithm

A block diagram of the GNSS/WSS integration algorithm is shown in Figure1. A two-dimensional DR based on WSS uses the speed measured at a pair of WSS of the rear wheels, and the procedure for this is as follows. When both the rear left and right WSS measure the wheel speed, the vehicle’s speed, vr, is calculated at the center of the rear wheel axle as follows [13]:(1)vr=vrL+vrR2,where vrL and vrR are the rear left and rear right WSS measurements, respectively.Figure 1
Block diagram of GNSS/WSS integration algorithm.To define the velocity in the body frame using the vehicle’s speed, we assume that the center of the body frame is set to the center of the rear wheel axle, the vehicle drives on a flat road, and no wheel slip occurs. In addition, assuming that the direction of the vehicle’s speed at the rear axle is equal to the longitudinal axis of the body frame, the velocities in the body frame are defined as follows:(2)vxvy=vr0,where vx and vy are the longitudinal and lateral velocity in the body frame, respectively.When a vehicle turns, the left and right wheel speeds are different. The speed measured by the individual WSS varies from along track while the vehicle is turning [1]. Therefore, the yaw rate could be calculated as the difference between the rear left wheel speed and the rear right wheel speed divided by the length between the rear wheels [13], as expressed in(3)ψ˙=vrL-vrRTr,where ψ˙ is the yaw rate of vehicle and Tr is the distance between the rear wheels.The yaw angle of the vehicle could be computed by numerical integration using the yaw at the previous time and the yaw rate, which is written as follows:(4)ψt=ψt-1+ψ˙Δt,where ψt and ψt-1 denote the yaw at time t and t-1, respectively, and Δt is the sampling time interval.The velocities in the navigation frame could be transformed by using(5)vnve=cos⁡ψ-sin⁡ψsin⁡ψcos⁡ψvxvy,where vn and ve denote the north velocity and the east velocity in the navigation frame, respectively, ψ is the yaw angle at the current time, and vx and vy are the longitudinal and lateral velocity in the body frame, respectively.In EKF-based GNSS/WSS integration, the state vector of the navigation error is composed of the latitude and longitude error, the north and east velocity error, and the yaw error. The state vector of the sensor error is composed of the average and the difference of the left and right WSS scale factor, defined as random constants. The white noise vector includes the white noise of WSS derived vehicle speed and the white noise of WSS derived yaw rate. The dynamic model for the GNSS/WSS integration is given in(6)δφ˙δλ˙δv˙nδv˙eδψ˙δs˙vδΔs˙v︸=δx˙WSS=001M+h00000vetan⁡φN+h0cos⁡φ001N+h0cos⁡φ0000000-vxsin⁡ψcos⁡ψ00000vxcos⁡ψsin⁡ψ000000ψ˙vxTr00000000000000︸=FWSSδφδλδvnδveδψδsvδΔsv︸=δxWSS+cos⁡ψM+h00sin⁡ψN+h0cos⁡φ00000010000︸=GWSSwvWSSwψ˙WSS︸=uWSS,where FWSS, δxWSS, GWSS, and uWSS are the dynamic matrix, the error state vector, the shaping matrix, and the white noise vector for GNSS/WSS integration, respectively, δφ and δλ are the latitude and longitude errors, respectively, δvn and δve are the north and east velocity errors, respectively, δψ is the yaw error, δsv and δΔsv are the average and the difference of the left and right WSS scale factors, respectively, M is the radius of curvature of the meridian, N is the radius of curvature of the prime vertical, h0 is the ellipsoidal height defined as a constant, vx and ve are the longitudinal and the east velocity, respectively, φ is the latitude, ψ is the yaw, Tr is the length between rear wheels, wvWSS is the white noise of WSS derived vehicle speed, and wψ˙WSS is the white noise of WSS derived yaw rate.The measurement model is generally expressed as follows:(7)z=Hδx+v,v~0,R,where z is the measurement vector, H is the design matrix, δx is the error state vector, v is the measurement error, and R is the variance-covariance matrix of the measurement error.In this study, the GNSS measurement vector consists of difference between the latitude, the longitude, the north velocity, and the east velocity estimated from DR based on vehicle dynamic sensors and GNSS receiver, as shown below:(8)zGNSS=φλvnveDR-φλvnveGNSS,where zGNSS is the GNSS measurement vector, the subscripts DR and GNSS denote the values computed by DR based on vehicle dynamic sensors and the values calculated from the GNSS receiver, respectively, φ and λ are the latitude and the longitude, respectively, and vn and ve are the north velocity and the east velocity, respectively.The design matrix and the variance-covariance matrix for GNSS measurements are shown, respectively, as follows:(9)HGNSS=I4×404×3,(10)RGNSS=σφ20000σλ20000σvn20000σve2,where HGNSS is the GNSS design matrix, RGNSS is the GNSS variance-covariance matrix, I4×4 is the 4×4 identity matrix, 04×3 is the 4×3 zero matrix, and σφ2, σλ2, σvn2, and σve2 are the variance of the estimated latitude, the longitude, the north velocity, and the east velocity from GNSS receiver, respectively.


## 2.3. GNSS/WSS/YRS Integration Algorithm

Figure2 shows a block diagram of the GNSS/WSS/YRS integration algorithm. Its structure is similar to that of the GNSS/WSS integration algorithm shown in Figure 1. The major difference between the DR mechanization based on WSS and the DR mechanization based on WSS/YRS is that the yaw rate value is measured from YRS. Therefore, in EKF-based GNSS/WSS/YRS integration, the state vector of the navigation error is the same as that in GNSS/YRS integration, and the state vector of the sensor error is composed of the average of the left and right WSS scale factor and the YRS bias, defined as random constants. The dynamic model for the GNSS/WSS/YRS integration is given in(11)δφ˙δλ˙δv˙nδv˙eδψ˙δs˙vδb˙YRS︸=δx˙WSS/YRS=001M+h00000vetan⁡φN+h0cos⁡φ001N+h0cos⁡φ0000000-vxsin⁡ψcos⁡ψ00000vxcos⁡ψsin⁡ψ0000000100000000000000︸=FWSS/YRSδφδλδvnδveδψδsvδbYRS︸=δxWSS/YRS+cos⁡ψM+h00sin⁡ψN+h0cos⁡φ00000010000︸=GWSS/YRSwvWSSwψ˙YRS︸=uWSS/YRS,where FWSS/YRS, δxWSS/YRS, GWSS/YRS, and uWSS/YRS are the dynamic matrix, the error state vector, the shaping matrix, and the white noise vector for GNSS/WSS/YRS integration, respectively,δbYRS is the YRS bias, and wψ˙YRS is the white noise of YRS.Figure 2
Block diagram of GNSS/WSS/YRS integration algorithm.Comparing the GNSS/WSS integration, the GNSS/WSS/YRS integration performs not only the GNSS measurement update but also the ZIHR (Zero Integrated Heading Rate) measurement update [3] to prevent the drift error of yaw when the vehicle is stationary. To update the error state vector by using GNSS based position and velocity, the GNSS measurement model is the same as that of GNSS/WSS integration, as shown in (8)–(10). When the vehicle stops, the measured speed by WSS is zero, and thus ZIHR can be applied to fix the yaw. The measurement vector for ZIHR is given by(12)zZIHR=ψk-ψk-1,where zZIHR is the ZIHR measurement vector, ψk is the estimated yaw at the current time, and ψk-1 is the estimated yaw at the previous time.The design matrix and the variance-covariance matrix for ZIHR are shown, respectively, as follows:(13)HZIHR=0000100,(14)RZIHR=σψ2,where HZIHR is the ZIHR design matrix, RZIHR is the ZIHR variance-covariance matrix, and σψ2 is the variance of the estimated yaw.


## 2.4. GNSS/WSS/YRS/SAS Integration Algorithm

Figure3 shows a block diagram of the GNSS/WSS/YRS/SAS integration algorithm. The idea of integrating SAS for DR is to compute the side slip angle by using the handle steering angle from SAS and then to employ the estimated side slip angle to set the longitudinal and lateral velocity from the WSS derived vehicle speed. In the case of the integration algorithm in Sections 2.2 and 2.3, a lateral constraint is applied to set the velocity in the body frame using WSS. However, in a land vehicle positioning system, since the lateral nonholonomic constraint is violated by a side slip during either cornering or change in road and tire conditions [14], the positioning accuracy would be degraded. Although the side slip angle is a very complicated phenomenon associated with road and tire conditions as well as high vehicle dynamics including fast driving and sharp turns [15], side slip always occurs during turning due to lateral tire deformation [16]. In order to define the lateral velocity by side slip in the two-dimensional DR based on WSS/YRS/SAS, the side slip angle model is defined in (15) with the assumption of a linear relationship between the side slip angle and the handle steering angle:(15)β=γ·α,where β is the side slip angle, γ is the side slip ratio, and α is the handle steering angle measured by SAS.Figure 3
Block diagram of GNSS/WSS/YRS/SAS integration algorithm.Considering the side slip angle, the velocities in the body frame are defined as follows:(16)vxvy=vrcos⁡βvrsin⁡β.In the two-dimensional DR based on WSS/YRS/SAS, the velocities in the navigation frame are given by (17)vnve=cos⁡ψ-sin⁡ψsin⁡ψcos⁡ψvxvy=vrcos⁡ψ+βvrsin⁡ψ+β.In EKF-based GNSS/WSS/YRS/SAS integration, the state vector of the navigation error is composed of the latitude and longitude error, the north and east velocity error, the yaw error, and the side slip angle error. The state vector of the sensor error consists of the average of the left and right WSS scale factor, the YRS bias, the SAS scale factor, and the SAS bias, defined as random constants. The dynamic model for the GNSS/WSS/YRS/SAS integration is expressed as follows:(18)δφ˙δλ˙δv˙nδv˙eδψ˙δβ˙δs˙vδb˙YRSδs˙SASδb˙SAS︸=δx˙WSS/YRS/SAS=001M+h00000000vetan⁡φN+h0cos⁡φ001N+h0cos⁡φ0000000000-vxsin⁡ψ-vycos⁡ψ-vxsin⁡ψ-vycos⁡ψcos⁡ψ+β0000000vxcos⁡ψ-vysin⁡ψvxcos⁡ψ-vysin⁡ψsin⁡ψ+β000000000010000000000β10000000000000000000000000000000000000000︸=FWSS/YRS/SASδφδλδvnδveδψδβδsvδbYRSδsSASδbSAS︸=δxWSS/YRS/SAS+cos⁡ψM+h000sin⁡ψN+h0cos⁡φ00000000010001000000000000︸=GWSS/YRS/SASwvWSSwψ˙YRSwβ˙SAS︸=uWSS/YRS/SAS,where FWSS/YRS/SAS, δxWSS/YRS/SAS, GWSS/YRS/SAS, and uWSS/YRS/SAS are the dynamic matrix, the error state vector, the shaping matrix, and the white noise vector for GNSS/WSS/YRS/SAS integration, respectively, δβ is the side slip angle error, δsSAS is the SAS scale factor, δbSAS is the SAS bias, and wβ˙YRS is the white noise of the side slip angle derived SAS.The measurement vector and the variance-covariance matrix for both the GNSS measurement model and the ZIHR measurement model are the same as that of GNSS/WSS/YRS integration. The design matrices for measurement models in the GNSS/WSS/YRS/SAS integration are expressed as follows:(19)HGNSS=I4×404×6,HZIHR=000010000000.


## 3. Performance Evaluation of GNSS/Vehicle Dynamic Sensor Integration Algorithm during GNSS Signal Blockages

### 3.1. Test Description

The test vehicle (Figure4) that is equipped with WSS, SAS, and YRS to operate ADAS is used for this study. The data from the vehicle dynamic sensors are acquired through a CAN (Controller Area Network) bus without installation of additional sensors. The specifications of the vehicle dynamic sensors are described in Table 3. The GNSS receiver used for acquiring the position and the velocity is a NoVATel DL-V3-GENERIC. The positioning method of the GNSS receiver was set to single point positioning using pseudorange and Doppler measurements from GPS and GLONASS. According to the specifications of the GNSS receiver, the precision of position and velocity are 1.5 m and 0.03 m/s. To evaluate the performance of the integration algorithms, the reference data were computed in postmission using the inertial, the GNSS, and DMI data obtained from Applanix’s POS LV 520. The precision of the horizontal position and yaw for reference data are 0.02 m and 0.015 deg, respectively.Table 3
Specifications of vehicle dynamic sensors.
Vehicle dynamic sensor
Output range
Resolution
Output rate
Wheel speed sensor
0~511.75 km/h
0.125 km/h
50 Hz
Yaw rate sensor
−40.95~40.95 deg/s
0.01 deg/s
100 Hz
Steering angle sensor
−3276.8~3276.6 deg
0.1 deg
100 HzFigure 4
Test vehicle.Two trajectories were used in the performance evaluation of the three GNSS/vehicle dynamic sensors integration algorithms during GNSS outages. An experiment for two trajectories data acquisition was conducted in the Daegu Technopolis area, Korea. These trajectories include various driving circumstances encountered during typical driving on urban roads with frequent stops, acceleration, deceleration, and speed bumps. However, since two trajectories were set in different paths in the same area, the two trajectories data were different in vehicle motion and driving condition. To compare the performance of three GNSS/Vehicle dynamic sensors integration algorithms during GNSS outages, fifteen GNSS outages of 30 seconds were simulated in each trajectory. The simulated GNSS outages covered a wide range of vehicle dynamics such as straight portions, turns, slopes, high speed, slow speeds, and jumping.


### 3.2. Performance Evaluation in the First Trajectory

Figure5 shows the reference data of the first trajectory with the simulated GNSS outages shown as a red line overlaid on the trajectory. The traveled distance of the first trajectory is about 11 km for 32 minutes. Figure 6 shows the velocity of the body frame and the attitude for reference data with the simulated GNSS outages shown as a red line. As can be seen in Figures 5 and 6, the GNSS outages (red line) were chosen by considering different dynamics motions.Figure 5
The first trajectory (red line: simulated GNSS outage).Figure 6
Velocity and attitude of the first trajectory (red line: simulated GNSS outage).
(a)
Velocity in body frame
(b)
AttitudeFirst, the drift error of the navigation solution during GNSS outages was analyzed in accordance with the integration strategies. Table4 compares the average RMS error of the horizontal position at the first and last epoch during GNSS outages by different integration strategies for DR navigation based on the combination of vehicle dynamic sensors. The corresponding results for horizontal velocity and yaw are summarized in Tables 5 and 6. During the GNSS outages, the error of GNSS/WSS integration drifts rapidly. However, the improvement of the horizontal position at the end of the 30-second GNSS outage is 24% using the YRS sensor as compared to the GNSS/WSS integration. Furthermore, the average RMS error of horizontal velocity at the end of the 30-second GNSS outage is improved by about 67%. This is due to the dependency of the limitation of WSS accuracy for yaw computation on factors such as wheel slipping, skidding, and slope. Comparing GNSS/WSS/YRS integration with GNSS/WSS/YRS/SAS integration, application of SAS based side slip angle in DR improves the positioning performance by about 20% during GNSS outages. Also, the average RMS error of horizontal velocity and yaw is reduced by 66% and 19%, respectively, from 0.06 m/s and 1.00 degrees for GNSS/WSS/YRS integration to 0.04 m/s and 0.81 degrees for GNSS/WSS/YRS/SAS integration. This indicates that the computation of the SAS based side slip angle could be employed to alleviate the effect of violating the lateral nonholonomic constraints.Table 4
Average RMS error of horizontal position during GNSS outages for first trajectory.
Integration strategies
Average of RMS error of horizontal position during GNSS outages [m]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
2.03
9.16
7.13
GNSS/WSS/YRS
1.77
7.18
5.41
GNSS/WSS/YRS/SAS
1.75
6.09
4.34Table 5
Average RMS error of horizontal velocity during GNSS outages for first trajectory.
Integration strategies
Average of RMS error of horizontal velocity during GNSS outages [m/s]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
0.25
0.43
0.18
GNSS/WSS/YRS
0.22
0.28
0.06
GNSS/WSS/YRS/SAS
0.22
0.26
0.04Table 6
Average RMS error of yaw during GNSS outages for first trajectory.
Integration strategies
Average of RMS error of yaw during GNSS outages [deg]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
1.07
2.61
1.54
GNSS/WSS/YRS
0.74
1.74
1.00
GNSS/WSS/YRS/SAS
0.60
1.41
0.81The maximum horizontal position error values during the fifteen outages are presented in Figure7. The cumulative relative frequency based on the maximum horizontal position error of 5 m and 10 m is summarized in Table 7. As expected, the integration strategy with all the vehicle dynamic sensors together (GPS/WSS/YRS/SAS) shows the overall best performance. Despite this, the maximum horizontal error calculated by the three GNSS/vehicle dynamic sensors integration algorithms was generally smaller than 10 m during the GNSS outage for 30 seconds. However, in the case of GNSS/WSS integration, the largest horizontal position error was observed for GNSS outage F. In this outage, the vehicle is traveling along a straight path, and the vehicle speed is greater than the other GNSS outages, as shown in both Figure 5 and Figure 6. Figure 8 shows the details of the horizontal position error, the longitudinal velocity error, and the yaw error during GNSS outages F. The longitudinal velocity error of three GNSS/Vehicle dynamic sensors integration algorithms were almost the same in the GNSS outages F. Despite this, the horizontal position error of the GNSS/WSS integration was higher than the other integrations due to the yaw error. As mentioned in Section 2, the major difference between GNSS/WSS integration and GNSS/WSS/YRS integration in DR mechanization is that the yaw computation in GNSS/WSS uses the calculated yaw rate based on the difference between the rear left wheel speed and the rear right wheel speed. If the vehicle is traveling along a straight path on flat road, the performance of GNSS/WSS integration is similar to the performance of GNSS/WSS/YRS integration because the yaw rate calculated by the rear wheels is zero in the ideal situation. However, since the wheel slipping and skidding in the left and right wheels could occurred differently due to the influence of either uneven road surface or vehicle dynamic force, the yaw rate calculated by the rear wheels is not zero in the real situation. This means that the accuracy of the computed yaw angle by using WSS could be significantly degraded. For this result, it is considered that WSS based DR has limited application to car positioning.Table 7
Cumulative relative frequency table of maximum horizontal position error for GNSS outages of first trajectory.
Integration strategies
Maximum horizontal position error
≤5 m
≤10 m
GNSS/WSS
27% (4 out of total 15 outages)
87% (13 out of total 15 outages)
GNSS/WSS/YRS
33% (5 out of total 15 outages)
93% (14 out of total 15 outages)
GNSS/WSS/YRS/SAS
53% (8 out of total 15 outages)
93% (14 out of total 15 outages)Figure 7
Maximum horizontal position error for GNSS outages of first trajectory.Figure 8
Details of the horizontal position error, the longitudinal velocity error, and the yaw error during GNSS outages F.


### 3.3. Performance Evaluation in the Second Trajectory

The test data for the second trajectory were acquired by driving a different path of 15.6 km for 39 min. The GNSS outages were selected by the same criteria. Figure9 shows the reference data of the second trajectory with the fifteen simulated GNSS outages. Both the velocity in the body frame and the attitude for reference data are shown in Figure 10. Compared with the first trajectory, the driving direction and the velocity for the second trajectory were changed more frequently.Figure 9
The second trajectory (red line: simulated GNSS outage).Figure 10
Velocity and attitude of the first trajectory (red line: simulated GNSS outage).
(a)
Velocity in body frame
(b)
AttitudeTable8 presents the average RMS error of the horizontal position at the first and last epoch during GNSS outages with respect to the integration strategies. Tables 9 and 10 summarize the average RMS error of the horizontal velocity and yaw at the first and last epoch during GNSS outages, respectively. The average RMS of the horizontal position of GNSS/WSS integration is significantly increased in comparison with the results of the first trajectory, as shown in Table 4. Despite this, the horizontal position RMS error of the GNSS/WSS integration of 13.62 m at the end of the 30-second GNSS outage was reduced to 4.22 m by using the WSS plus YRS, constituting a 69% improvement. The horizontal velocity drift error can be reduced from 0.81 m/s for GNSS/WSS integration to 0.20 m/s through the integration with YRS, and the yaw drift error is reduced from 5.23 degrees to 2.52 degrees. The horizontal velocity percentage improvement is 75% and the yaw percentage improvement is 52%. Therefore, using the yaw rate measured by YRS instead of that computed by WSS greatly enhances the DR accuracy during GNSS outages. The average of RMS horizontal position errors of the GNSS/WSS/YRS/SAS integration at the end of the 30-second GNSS outage is improved by about 7% compared to the GNSS/WSS/YRS solution. Furthermore, for the GNSS/WSS/YRS integration, the enhancement at the end of the 30-second GNSS outages is 15% for the average RMS error of horizontal velocity and 2% for the average RMS error of the yaw case as compared to the GNSS/WSS/YRS solution. Although the improvement gained from additional integration of the SAS is slight, using the SAS based side slip angle for DR can reduce the drift error.Table 8
Average RMS error of horizontal position during GNSS outages for second trajectory.
Integration strategies
Average of RMS error of horizontal position during GNSS outages [m]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
1.91
15.54
13.62
GNSS/WSS/YRS
1.65
5.86
4.22
GNSS/WSS/YRS/SAS
1.64
5.57
3.93Table 9
Average RMS error of horizontal velocity during GNSS outages for second trajectory.
Integration strategies
Average of RMS error of horizontal velocity during GNSS outages [m/s]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
0.16
0.97
0.81
GNSS/WSS/YRS
0.12
0.32
0.20
GNSS/WSS/YRS/SAS
0.12
0.28
0.17Table 10
Average RMS error of yaw during GNSS outages for second trajectory.
Integration strategies
Average of RMS error of yaw during GNSS outages [deg]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
0.52
5.75
5.23
GNSS/WSS/YRS
0.36
2.88
2.52
GNSS/WSS/YRS/SAS
0.34
2.80
2.46Figure11 shows the maximum horizontal position error during the fifteen outages. The cumulative relative frequency based on the maximum horizontal position error of 5 m and 10 m is summarized in Table 11. The maximum horizontal position error for GNSS/WSS integration is smaller than 10 m in nine out of the total 15 GNSS outages (60%), which is worse than the result of the first trajectory. It can be seen that the accuracy of the computed yaw rate using WSS is easily degraded by typical conditions encountered by a land vehicle, including uneven road, speed bumps, high speed, and sudden acceleration and deceleration. As seen in Figure 10, for the velocity and attitude of GNSS outages A, C, E, H, J, K, and L, where the maximum horizontal error of GNSS/WSS integration is higher than the other integrations, the accuracy of the computed yaw rate using WSS is easily degraded by typical driving conditions encountered by a land vehicle, including uneven road, speed bumps, high speed, and sudden acceleration, and deceleration. Therefore, it is estimated that the GNSS/WSS integration would provide an unstable navigation solution during GNSS outages. The maximum horizontal position error for GNSS/WSS/YRS integration as well as GNSS/WSS/YRS/SAS integration is smaller than 10 m in fourteen out of the total 15 GNSS outages (93%), which is equal to that of the first trajectory. Also, GNSS/WSS/YRS/SAS integration achieves the overall best performance.Table 11
Cumulative relative frequency table of maximum horizontal position error for GNSS outages of second trajectory.
Integration strategies
Maximum horizontal position error
≤5 m
≤10 m
GNSS/WSS
27% (4 out of total 15 outages)
60% (9 out of total 15 outages)
GNSS/WSS/YRS
53% (8 out of total 15 outages)
93% (14 out of total 15 outages)
GNSS/WSS/YRS/SAS
67% (10 out of total 15 outages)
93% (14 out of total 15 outages)Figure 11
Maximum horizontal position error for GNSS outages of second trajectory.


## 3.1. Test Description

The test vehicle (Figure4) that is equipped with WSS, SAS, and YRS to operate ADAS is used for this study. The data from the vehicle dynamic sensors are acquired through a CAN (Controller Area Network) bus without installation of additional sensors. The specifications of the vehicle dynamic sensors are described in Table 3. The GNSS receiver used for acquiring the position and the velocity is a NoVATel DL-V3-GENERIC. The positioning method of the GNSS receiver was set to single point positioning using pseudorange and Doppler measurements from GPS and GLONASS. According to the specifications of the GNSS receiver, the precision of position and velocity are 1.5 m and 0.03 m/s. To evaluate the performance of the integration algorithms, the reference data were computed in postmission using the inertial, the GNSS, and DMI data obtained from Applanix’s POS LV 520. The precision of the horizontal position and yaw for reference data are 0.02 m and 0.015 deg, respectively.Table 3
Specifications of vehicle dynamic sensors.
Vehicle dynamic sensor
Output range
Resolution
Output rate
Wheel speed sensor
0~511.75 km/h
0.125 km/h
50 Hz
Yaw rate sensor
−40.95~40.95 deg/s
0.01 deg/s
100 Hz
Steering angle sensor
−3276.8~3276.6 deg
0.1 deg
100 HzFigure 4
Test vehicle.Two trajectories were used in the performance evaluation of the three GNSS/vehicle dynamic sensors integration algorithms during GNSS outages. An experiment for two trajectories data acquisition was conducted in the Daegu Technopolis area, Korea. These trajectories include various driving circumstances encountered during typical driving on urban roads with frequent stops, acceleration, deceleration, and speed bumps. However, since two trajectories were set in different paths in the same area, the two trajectories data were different in vehicle motion and driving condition. To compare the performance of three GNSS/Vehicle dynamic sensors integration algorithms during GNSS outages, fifteen GNSS outages of 30 seconds were simulated in each trajectory. The simulated GNSS outages covered a wide range of vehicle dynamics such as straight portions, turns, slopes, high speed, slow speeds, and jumping.


## 3.2. Performance Evaluation in the First Trajectory

Figure5 shows the reference data of the first trajectory with the simulated GNSS outages shown as a red line overlaid on the trajectory. The traveled distance of the first trajectory is about 11 km for 32 minutes. Figure 6 shows the velocity of the body frame and the attitude for reference data with the simulated GNSS outages shown as a red line. As can be seen in Figures 5 and 6, the GNSS outages (red line) were chosen by considering different dynamics motions.Figure 5
The first trajectory (red line: simulated GNSS outage).Figure 6
Velocity and attitude of the first trajectory (red line: simulated GNSS outage).
(a)
Velocity in body frame
(b)
AttitudeFirst, the drift error of the navigation solution during GNSS outages was analyzed in accordance with the integration strategies. Table4 compares the average RMS error of the horizontal position at the first and last epoch during GNSS outages by different integration strategies for DR navigation based on the combination of vehicle dynamic sensors. The corresponding results for horizontal velocity and yaw are summarized in Tables 5 and 6. During the GNSS outages, the error of GNSS/WSS integration drifts rapidly. However, the improvement of the horizontal position at the end of the 30-second GNSS outage is 24% using the YRS sensor as compared to the GNSS/WSS integration. Furthermore, the average RMS error of horizontal velocity at the end of the 30-second GNSS outage is improved by about 67%. This is due to the dependency of the limitation of WSS accuracy for yaw computation on factors such as wheel slipping, skidding, and slope. Comparing GNSS/WSS/YRS integration with GNSS/WSS/YRS/SAS integration, application of SAS based side slip angle in DR improves the positioning performance by about 20% during GNSS outages. Also, the average RMS error of horizontal velocity and yaw is reduced by 66% and 19%, respectively, from 0.06 m/s and 1.00 degrees for GNSS/WSS/YRS integration to 0.04 m/s and 0.81 degrees for GNSS/WSS/YRS/SAS integration. This indicates that the computation of the SAS based side slip angle could be employed to alleviate the effect of violating the lateral nonholonomic constraints.Table 4
Average RMS error of horizontal position during GNSS outages for first trajectory.
Integration strategies
Average of RMS error of horizontal position during GNSS outages [m]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
2.03
9.16
7.13
GNSS/WSS/YRS
1.77
7.18
5.41
GNSS/WSS/YRS/SAS
1.75
6.09
4.34Table 5
Average RMS error of horizontal velocity during GNSS outages for first trajectory.
Integration strategies
Average of RMS error of horizontal velocity during GNSS outages [m/s]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
0.25
0.43
0.18
GNSS/WSS/YRS
0.22
0.28
0.06
GNSS/WSS/YRS/SAS
0.22
0.26
0.04Table 6
Average RMS error of yaw during GNSS outages for first trajectory.
Integration strategies
Average of RMS error of yaw during GNSS outages [deg]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
1.07
2.61
1.54
GNSS/WSS/YRS
0.74
1.74
1.00
GNSS/WSS/YRS/SAS
0.60
1.41
0.81The maximum horizontal position error values during the fifteen outages are presented in Figure7. The cumulative relative frequency based on the maximum horizontal position error of 5 m and 10 m is summarized in Table 7. As expected, the integration strategy with all the vehicle dynamic sensors together (GPS/WSS/YRS/SAS) shows the overall best performance. Despite this, the maximum horizontal error calculated by the three GNSS/vehicle dynamic sensors integration algorithms was generally smaller than 10 m during the GNSS outage for 30 seconds. However, in the case of GNSS/WSS integration, the largest horizontal position error was observed for GNSS outage F. In this outage, the vehicle is traveling along a straight path, and the vehicle speed is greater than the other GNSS outages, as shown in both Figure 5 and Figure 6. Figure 8 shows the details of the horizontal position error, the longitudinal velocity error, and the yaw error during GNSS outages F. The longitudinal velocity error of three GNSS/Vehicle dynamic sensors integration algorithms were almost the same in the GNSS outages F. Despite this, the horizontal position error of the GNSS/WSS integration was higher than the other integrations due to the yaw error. As mentioned in Section 2, the major difference between GNSS/WSS integration and GNSS/WSS/YRS integration in DR mechanization is that the yaw computation in GNSS/WSS uses the calculated yaw rate based on the difference between the rear left wheel speed and the rear right wheel speed. If the vehicle is traveling along a straight path on flat road, the performance of GNSS/WSS integration is similar to the performance of GNSS/WSS/YRS integration because the yaw rate calculated by the rear wheels is zero in the ideal situation. However, since the wheel slipping and skidding in the left and right wheels could occurred differently due to the influence of either uneven road surface or vehicle dynamic force, the yaw rate calculated by the rear wheels is not zero in the real situation. This means that the accuracy of the computed yaw angle by using WSS could be significantly degraded. For this result, it is considered that WSS based DR has limited application to car positioning.Table 7
Cumulative relative frequency table of maximum horizontal position error for GNSS outages of first trajectory.
Integration strategies
Maximum horizontal position error
≤5 m
≤10 m
GNSS/WSS
27% (4 out of total 15 outages)
87% (13 out of total 15 outages)
GNSS/WSS/YRS
33% (5 out of total 15 outages)
93% (14 out of total 15 outages)
GNSS/WSS/YRS/SAS
53% (8 out of total 15 outages)
93% (14 out of total 15 outages)Figure 7
Maximum horizontal position error for GNSS outages of first trajectory.Figure 8
Details of the horizontal position error, the longitudinal velocity error, and the yaw error during GNSS outages F.


## 3.3. Performance Evaluation in the Second Trajectory

The test data for the second trajectory were acquired by driving a different path of 15.6 km for 39 min. The GNSS outages were selected by the same criteria. Figure9 shows the reference data of the second trajectory with the fifteen simulated GNSS outages. Both the velocity in the body frame and the attitude for reference data are shown in Figure 10. Compared with the first trajectory, the driving direction and the velocity for the second trajectory were changed more frequently.Figure 9
The second trajectory (red line: simulated GNSS outage).Figure 10
Velocity and attitude of the first trajectory (red line: simulated GNSS outage).
(a)
Velocity in body frame
(b)
AttitudeTable8 presents the average RMS error of the horizontal position at the first and last epoch during GNSS outages with respect to the integration strategies. Tables 9 and 10 summarize the average RMS error of the horizontal velocity and yaw at the first and last epoch during GNSS outages, respectively. The average RMS of the horizontal position of GNSS/WSS integration is significantly increased in comparison with the results of the first trajectory, as shown in Table 4. Despite this, the horizontal position RMS error of the GNSS/WSS integration of 13.62 m at the end of the 30-second GNSS outage was reduced to 4.22 m by using the WSS plus YRS, constituting a 69% improvement. The horizontal velocity drift error can be reduced from 0.81 m/s for GNSS/WSS integration to 0.20 m/s through the integration with YRS, and the yaw drift error is reduced from 5.23 degrees to 2.52 degrees. The horizontal velocity percentage improvement is 75% and the yaw percentage improvement is 52%. Therefore, using the yaw rate measured by YRS instead of that computed by WSS greatly enhances the DR accuracy during GNSS outages. The average of RMS horizontal position errors of the GNSS/WSS/YRS/SAS integration at the end of the 30-second GNSS outage is improved by about 7% compared to the GNSS/WSS/YRS solution. Furthermore, for the GNSS/WSS/YRS integration, the enhancement at the end of the 30-second GNSS outages is 15% for the average RMS error of horizontal velocity and 2% for the average RMS error of the yaw case as compared to the GNSS/WSS/YRS solution. Although the improvement gained from additional integration of the SAS is slight, using the SAS based side slip angle for DR can reduce the drift error.Table 8
Average RMS error of horizontal position during GNSS outages for second trajectory.
Integration strategies
Average of RMS error of horizontal position during GNSS outages [m]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
1.91
15.54
13.62
GNSS/WSS/YRS
1.65
5.86
4.22
GNSS/WSS/YRS/SAS
1.64
5.57
3.93Table 9
Average RMS error of horizontal velocity during GNSS outages for second trajectory.
Integration strategies
Average of RMS error of horizontal velocity during GNSS outages [m/s]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
0.16
0.97
0.81
GNSS/WSS/YRS
0.12
0.32
0.20
GNSS/WSS/YRS/SAS
0.12
0.28
0.17Table 10
Average RMS error of yaw during GNSS outages for second trajectory.
Integration strategies
Average of RMS error of yaw during GNSS outages [deg]
First epoch
Last epoch
Difference (last − first)
GNSS/WSS
0.52
5.75
5.23
GNSS/WSS/YRS
0.36
2.88
2.52
GNSS/WSS/YRS/SAS
0.34
2.80
2.46Figure11 shows the maximum horizontal position error during the fifteen outages. The cumulative relative frequency based on the maximum horizontal position error of 5 m and 10 m is summarized in Table 11. The maximum horizontal position error for GNSS/WSS integration is smaller than 10 m in nine out of the total 15 GNSS outages (60%), which is worse than the result of the first trajectory. It can be seen that the accuracy of the computed yaw rate using WSS is easily degraded by typical conditions encountered by a land vehicle, including uneven road, speed bumps, high speed, and sudden acceleration and deceleration. As seen in Figure 10, for the velocity and attitude of GNSS outages A, C, E, H, J, K, and L, where the maximum horizontal error of GNSS/WSS integration is higher than the other integrations, the accuracy of the computed yaw rate using WSS is easily degraded by typical driving conditions encountered by a land vehicle, including uneven road, speed bumps, high speed, and sudden acceleration, and deceleration. Therefore, it is estimated that the GNSS/WSS integration would provide an unstable navigation solution during GNSS outages. The maximum horizontal position error for GNSS/WSS/YRS integration as well as GNSS/WSS/YRS/SAS integration is smaller than 10 m in fourteen out of the total 15 GNSS outages (93%), which is equal to that of the first trajectory. Also, GNSS/WSS/YRS/SAS integration achieves the overall best performance.Table 11
Cumulative relative frequency table of maximum horizontal position error for GNSS outages of second trajectory.
Integration strategies
Maximum horizontal position error
≤5 m
≤10 m
GNSS/WSS
27% (4 out of total 15 outages)
60% (9 out of total 15 outages)
GNSS/WSS/YRS
53% (8 out of total 15 outages)
93% (14 out of total 15 outages)
GNSS/WSS/YRS/SAS
67% (10 out of total 15 outages)
93% (14 out of total 15 outages)Figure 11
Maximum horizontal position error for GNSS outages of second trajectory.


## 4. Conclusions

This study presents the performance evaluation of a land vehicle positioning system encompassing GNSS combined with a two-dimensional DR based on vehicle dynamic sensors. To develop GNSS/vehicle dynamic sensor based positioning algorithms, vehicle dynamic sensors used WSS, YRS, and SAS, which were already installed in the test vehicle. Three two-dimensional DR mechanisms were designed to ensure that the vehicle dynamic sensors’ information is not duplicated. The GNSS/vehicle dynamic sensors integrations were implemented by EKF through a loosely coupled mode. The GNSS measurement are the position and the velocity calculated using the C/A code and the Doppler measurements from GPS and GLONASS at a GNSS receiver. The developed algorithms were tested on two trajectories acquired in various driving circumstances. A performance evaluation was conducted in fifteen simulated GNSS outages during 30 seconds for each trajectory. The results indicate that the integration algorithm with all the vehicle dynamic sensors together (GPS/WSS/YRS/SAS) provided the best performance. With respect to two trajectories, the maximum horizontal position error of both GNSS/WSS/YRS and GNSS/WSS/YRS/SAS integration was smaller than 10 m in 28 out of the total 30 GNSS outages. And the maximum horizontal positon error of GNSS/WSS integration was smaller than 10 m in 22 out of the total 30 GNSS outages. It is estimated that the GNSS/WSS integration would provide an unstable navigation solution during GNSS outages in comparison to both GNSS/WSS/YRS and GNSS/WSS/YRS/SAS integration since the accuracy of the computed yaw angle by using WSS could be significantly degraded due to frequent wheel slipping and skidding. Therefore, the proposed GNSS/vehicle dynamics sensor integrations excluding GNSS/WSS integration could be applied to an automotive navigation system with meter-level accuracy to overcome the limitations of the GNSS based positioning technique. Indeed, these algorithms might provide alternative solutions to the use of a low-cost MEMS-based IMU.



---

*Source: 9802610-2017-09-07.xml*