The Performance Analysis of a Real-Time Integrated INS/GPS Vehicle Navigation System with Abnormal GPS Measurement Elimination

The integration of an Inertial Navigation System (INS) and the Global Positioning System (GPS) is common in mobile mapping and navigation applications to seamlessly determine the position, velocity, and orientation of the mobile platform. In most INS/GPS integrated architectures, the GPS is considered to be an accurate reference with which to correct for the systematic errors of the inertial sensors, which are composed of biases, scale factors and drift. However, the GPS receiver may produce abnormal pseudo-range errors mainly caused by ionospheric delay, tropospheric delay and the multipath effect. These errors degrade the overall position accuracy of an integrated system that uses conventional INS/GPS integration strategies such as loosely coupled (LC) and tightly coupled (TC) schemes. Conventional tightly coupled INS/GPS integration schemes apply the Klobuchar model and the Hopfield model to reduce pseudo-range delays caused by ionospheric delay and tropospheric delay, respectively, but do not address the multipath problem. However, the multipath effect (from reflected GPS signals) affects the position error far more significantly in a consumer-grade GPS receiver than in an expensive, geodetic-grade GPS receiver. To avoid this problem, a new integrated INS/GPS architecture is proposed. The proposed method is described and applied in a real-time integrated system with two integration strategies, namely, loosely coupled and tightly coupled schemes, respectively. To verify the effectiveness of the proposed method, field tests with various scenarios are conducted and the results are compared with a reliable reference system.


Introduction
The past two decades have seen an increase in the use of positioning and navigation technologies in land vehicle applications. Applications of these technologies in land transportation are numerous, including automated car navigation, emergency assistance, fleet management, person/asset tracking, collision avoidance, environmental monitoring, and automotive assistance. The convergence of location, information management, and communication technologies has created a rapidly emerging market known as location-based services (LBS). LBS is a critical enabling technology that uses location as a filter or magnet to extract relevant information to provide value-added services such as location-aware billing, automated advertising services, and other location-based information sought by the user based on their location. Because of the importance of location information, the market, in turn, has pushed hard for the development of reliable vehicle navigation and guidance systems that provide not only location information but also route guidance and location-sensitive services. Virtually all modern land-vehicle navigation systems integrate two or more complimentary positioning technologies to provide the vehicle's position, velocity, and heading information in a seamless fashion. Typical candidates for these integrated navigation systems are the Global Positioning System (GPS) and Inertial Navigation Systems (INS). GPS can provide continuous, accurate positioning given clear lines of sight to more than four satellites. However, the accuracy and the availability of GPS-based vehicle navigation systems are subject to the open-sky conditions and degrade in the presence of signal blockage and reflected signals, as shown in Figures 1 and 2, respectively. The GPS trajectories shown in these figures were obtained from the real-time solutions of a consumer-grade GPS receiver, the AEK 4T from Ublox, and the reference trajectories were obtained from the post-processed, smoothed solutions using a tactical-grade integrated INS/GPS system aided by an odometer, the SPAN-CPT from NovAtel.
When GPS signals are unavailable, the INS can fill in the gaps to provide continuous navigation solutions (position, velocity, and attitude). An Inertial Measurement Unit (IMU) generally contains a set of inertial sensors, e.g., gyroscopes and accelerometers, to provide raw measurements including velocity changes (Δv) and orientation changes (Δθ) in the three directions of its body-fixed coordinate frame. The term "INS" usually refers to an IMU combined with an onboard computer that can provide navigation solutions in the chosen navigation frame directly in real time in addition to compensated raw measurements [1].
Although an integrated navigation system can work in GPS-denied environments, problems include the cost of the inertial sensors and the length of time that the GPS signals are unavailable, which affect its applicability. Tactical-grade or better inertial systems can achieve sufficient position accuracy and sustainability during long-duration GPS signal blockages [1,2]. For example, the high-end, expensive systems can provide less than 3 m real-time position accuracy with a GPS gap lasting one minute. However, the cost of the sophisticated inertial sensors is prohibitive for applications such as the primary navigation module for general land vehicles. For this reason, strap-down Micro-Electro-Mechanical Systems (MEMS) inertial sensors are preferred as the complementary component to GPS for general, seamless vehicle navigation applications. However, the position accuracy of these low-cost inertial sensors degrades rapidly with time when GPS signals are interrupted. The sustainability of an integrated INS/GPS system using currently available commercial MEMS inertial technology in typical GPS-denied environments is thus limited. However, the progress in MEMS inertial sensors has advanced rapidly. Thus, the inclusion of MEMS inertial sensors for general land-vehicle navigation has considerable potential in terms of cost and accuracy [3].  For INS/GPS integration applications, several architectures have been developed [4]. The most common integration scheme used today is the loosely coupled (LC) integration scheme. It is the simplest way of integrating the GPS processing engine into an integrated navigation system. The GPS processing engine calculates position fixes and velocities in the local level frame and then sends the solutions as measurement updates to the main Extended Kalman Filter (EKF). By comparing the navigation solutions provided by the INS mechanization with those provided by the GPS processing engine, the navigation states can be optimally estimated. As shown in Figure 3, the primary advantage of the LC architecture is the simplicity of its implementation; no advanced knowledge of processing GPS measurements is necessary. The disadvantage of this implementation is that the measurement update of the integrated navigation system is possible only when four or more satellites are in view [1,2,4]. The tightly coupled (TC) integration scheme uses a single EKF to integrate the GPS and IMU measurements. In the TC integration scheme, the GPS pseudo-range, Doppler signal, and carrier-phase measurements are processed directly in the main EKF, as shown in Figure 4. The primary advantage of this architecture is that the raw GPS measurements can be used to update the INS when fewer than four satellites are available. This scheme is particularly beneficial in problematic environments such as downtown city areas, where the reception of GPS satellite signals may be impeded by obstructions. However, according to Chiang and Huang [5], the EKF implemented with the conventional TC scheme may have serious problems related to the quality of the raw GPS measurements, especially when a low-cost, consumer-grade GPS receiver is used. Most of the expensive, geodetic-grade GPS receivers come with a sophisticated correlator and an antenna designed to manage reflected signals. These expensive GPS receivers usually fail to track a sufficient number of GPS satellites to provide a continuous trajectory when operating in GPS-degraded environments such as urban areas or forests because the reflected signals are actively eliminated [1]. Conversely, the low-cost, consumer-grade GPS receivers tend to provide a continuous although less accurate trajectory because they make use of all the available GPS signals, even the reflected ones that reduce the overall accuracy of the continuous trajectory, as shown in Figures 1 and 2. Conventional, EKF-based TC architectures are sensitive to the quality of the raw GPS measurements. This problem is usually observed in urban and suburban areas because of the reflected GPS signals [4,5].
Real-time INS/GPS systems using MEMS IMUs have been investigated with various computation platforms [6,7]. Liu [8] analyzed the performance of real-time, low-cost MEMS INS/GPS integration with additional constraints in land-based mobile mapping systems (MMSs). In the present study, a land-based, GPS-aided, MEMS inertial navigator is developed using a low-cost, consumer-grade GPS receiver with a modified version of a real-time, TC INS/GPS integration scheme to detect and eliminate abnormal GPS signals caused by multipath and other uncertainties to improve the stability and reliability of the system. In addition, tests using various scenarios are performed to evaluate the performance of the proposed system with actual data.

Configuration of Proposed Land-Based Global Positioning System Aided Inertial Navigator
The configuration of the proposed land-based, GPS-aided, MEMS inertial navigator and the real-time, modified TC INS/GPS integration scheme are described in detail in this section.

Inertial Navigation System (INS)/Global Positioning System (GPS) Integration
A modified TC INS/GPS integration is proposed in this research as described in the Figure 5. First, the outputs of the IMU, the angular rates sensed by gyroscopes and the specific forces sensed by accelerometers are processed by the INS mechanization to obtain navigation solutions, which are the position, the velocity, and the attitude in the navigation frame. The pseudo-range and Doppler measurements, the raw measurements from the GPS, are pre-processed to eliminate systematic errors such as the satellite clock, tropospheric, and ionospheric errors. Based on the navigation information estimated by the INS mechanization, GPS measurements are then checked to detect and eliminate abnormal measurements before the data-fusion engine is used. This process will be discussed in detail in Section 3. In addition, as a design of a closed-loop, the biases and the scale factors of the accelerometers and the gyroscopes are included in the state vector. Thus, during estimation, these errors are estimated to compensate the raw measurements coming from the IMU. The EKF is used as an estimator for data fusion to obtain estimated navigation solutions.

Additional Aiding Sources
In the presence of GPS outages, additional aids are essential to augment the system to maintain the stability and the accuracy of the navigation solutions. As illustrated in [8], the use of a Zero Velocity Update and a Zero Integrated Heading Rate (ZUPT/ZIHR) can significantly improve the accuracy of the navigation solutions during GPS signal blockages. ZUPT means the occasion stop of the system for short duration for estimating errors of the system and thus bounding the growth of inertial sensor errors. If vehicle stops, the velocity outputs in any directions should be zero. On the other hand, the ZIHR utilize the fact that the heading angle of the vehicle is not changed during the vehicle stops. However, the frequent stopping required for implementing ZUPT/ZIHR is not feasible in practice. In this study, only the Non-Holonomic Constraint (NHC) is applied as an aid during GPS signal outages.
As mentioned in [9], in a land vehicle platform, it is assumed that if the vehicle does not jump off the ground or slide sideways under normal conditions, so the velocities of the vehicle in the plane perpendicular to the forward direction are approximately zero. This assumption becomes a constraint condition for land-based navigation applications. In terms of implementation, the velocity components in the y and z directions in the body frame will be zero, as given in Equation (1) and illustrated in Figure 6: where the superscript (b) denotes the body frame To derive the measurement update equation for the KF, the velocity in the navigation frame is converted into the body frame: (2) where superscript (n) denotes the navigation frame and is the rotation matrix from the navigation frame to the body frame. For the EKF, in which the system model is derived from the error state model, the corresponding error equation to Equation (2) is derived as follows: is the velocity error vector in the body frame, is the velocity error vector in the navigation frame, is the cross-product form of the velocity vector in the navigation frame, and is the attitude error vector. The measurement equation can be constructed as follows: (4) where and are velocity noise in the y and z directions, respectively.
The NHC is an analytic correction; no additional sensor is required; therefore, it can be applied to any land-based integrated navigation system to improve the navigation accuracy. However, if the assumption of the vehicle behavior is violated, the NHC may cause more noise to the system. Normally, under an open sky, the GPS is more reliable than the NHC. Therefore, in the proposed system, the NHC is activated only when GPS signal outages take place. In addition, the update interval of the NHC is subject to change depending on the quality of the IMU: the higher the IMU quality, the longer the update interval of the NHC should be.

Model Design for Estimation
Models of the system and the measurements are required for data fusion using EKF [2]. With a high sampling rate and a seamless output, the INS is used to form the system model. GPS measurements and other external aids are used to build the measurement models.

System Model
The system model is derived based on the error model of a strap-down INS in a navigation frame. The "PSI" model is chosen for the integrated navigation system because of its simple attitude error dynamic equation. Details of the derivation of the "PSI" model can be found in [2]. The core of the system dynamic model using the PSI angle error model is expressed in time-continuous form as below: (5) where , and are the time derivative of the position, velocity and attitude error vector in the computer frame (the local level frame at the location erroneously estimated by the INS mechanization), respectively, is rotation matrix from b-frame to the n-frame, and are the specific force and angular rate error vector: where is the specific force, g is the acceleration of gravity, is the rotation rate vector of the navigation frame (n-frame) relative to the earth frame, is the rotation rate vector of the earth frame relative to the inertial frame expressed in computer frame, r e is the earth radius, h is the ellipsoid high. Equation (5) can be rewritten as follows: (7) where: (8) For the closed-loop architecture, the biases and the scale factors of the gyroscopes and the accelerometers are considered as unknown parameters. These parameters are included in the state vector to be estimated in every computational cycle of the KF. The models for these parameters are as follows: where , , and are biases and scale factors of the accelerometers and the gyroscopes, respectively, , , , and are the discrete-time sensor error model parameters that can be determined from the first term of Gauss-Markov processes. Let: ; Equation (9) can be rewritten as follows: Combining Equations (7) and (12) yields the following: Equation (13) can be rewritten as follows: Equation (14) is transformed into a discrete-time form based on [2]: (14) where is the state vector at time (epoch) k, is the state transition matrix from epoch k − 1 to k, and w k is the system noise. In KF-based, the system noise is assumed to have normal distribution with zero mean and variance Q k , . In fact, the system noise is not usually correctly estimated but can be modeled by the representation of the system noise model k . In theory, , however, because the is not directly estimated at every epoch k, the is normally fixed overtime and is modeled by IMU calibration based on Gauss-Markov processes. The IMU calibration process can be found in [10].

Measurement Model
For the LC scheme, the aiding measurements are position and velocity provided by the GPS receiver. For the EKF, the measurement model is as follows: ( 15) where are mapping matrices, is the position error vector, is the velocity vector expressed in the earth-centered, earth-fixed frame (ECEF or e-frame), and and are the position and velocity noise, respectively. Identical to the system noise, the measurement noises vector is not directly estimated, but is accounted in the EKF by a representation of the measurement noise model R, where R is modeled based on the GPS position and velocity uncertainty.
For the TC scheme, the pseudo-range and Doppler measurements are used to aid the INS mechanization. Pseudo-range is the erroneous rage from the GPS receiver to the satellite. The original pseudo-range measurement equation derived by [11] is shown as: where  is the pseudo-range, is the true range between the receiver and the satellite, is the satellite clock error, is the clock synchronization error, I  and T  are the delays associated with the transmission of the signal through the ionosphere and troposphere, respectively, c is the speed of light, and is pseudo-range noise.
In this research ionospheric, tropospheric and satellite clock errors are estimated before using pseudo-range in EKF. The Klobuchar model [12] is used to account for ionospheric error and Hofield model [13] is used to eliminate tropospheric error. The satellite clock error can be estimated based on information from ephemeris message. After error correction, the Equation (17) becomes: Doppler measurement (Doppler shift) is the difference between the frequency of the received signal (carrier phase received at the receiver) and the frequency at the source (carrier phase transmitted at the satellite) [11]. The equation of the Doppler measurement is expressed as: where f r is the received frequency from the receiver, f s is the original frequency transmitted by the satellite, is the change line of sight distance between the receiver and the satellite, is the signal wavelength ( , where c is the speed of light, f L1 is the GPS L1 carrier frequency), is the Doppler measurement noise. Account for the GPS receiver clock drift error, the Doppler measurement equation can be rewritten as: where is line of sight vector, is the receiver velocity, and is the satellite velocity, is the receiver clock drift error.
From Equations (18) and (20), with m satellites observed, the pseudo-range and Doppler measurements can be modeled for EKF as: where the superscript e denotes the e-frame, , and are the projected components of the line of sight to satellite k in the x, y, and z axes in the e-frame, respectively, and are the predictions of the pseudo-range and Doppler measurements by the INS with satellite k and determined as: From Equations (16), (21) and (22), after some transformations, the measurement model can be expressed in a form required by the EKF: In the TC INS/GPS integration scheme with GPS pseudo-range and Doppler phase measurements, the receiver clock error and the clock drift error are included in the system dynamic model: where and are the white noise of the receiver clock error and the clock drift error. The state vector becomes:

Data Fusion Strategies
This section first describes the use of the EKF to estimate the state variables and their corresponding covariances based on the system and measurement models. Subsequently, a strategy is proposed to eliminate abnormal GPS measurements.
Normally, the updated estimates are better than the estimates at the prediction step. Thus, the estimated biases and scale factors in the update steps are used to compensate the IMU measurements in subsequent epochs: where are the increments in the angle and the velocity at time k + 1, and are the increments in the angle and the velocity estimated at time k, respectively.
As discussed in Section 1, the accuracy of the GPS measurements is independent of time but dependent on the operating environment. This means that the GPS measurements may deteriorate at any time in a GPS-degraded environment, such as in an urban canyon or a tunnel. In contrast, the solutions from the INS are continuously provided in any environment and are accurate in the short term. Using this characteristic, the predicted states based on the INS mechanization are used to decide if the GPS measurements are sufficiently reliable to be used to perform measurement update in the EKF. This problem is considered for both the LC and TC schemes. If the position or the pseudo-range is biased because of abnormal signals, then the velocity and Doppler measurement are also biased. Therefore, only the position and the pseudo-range are considered for abnormal GPS measurement rejection in these schemes, respectively. Two methods of abnormal GPS measurement rejection are considered in this research, an analytic method and a statistical method.

Analytic Method for Rejecting Abnormal Global Positioning System Measurements
In the LC scheme, the difference vector that forms the aiding position error measurement vector in Equation (16)  In the modified TC mode, the difference between the measured GPS pseudo-range and the pseudo-range estimated by the INS is considered, i.e.,: An error criterion is as follows: where is the absolute value of : can be derived from Equation (23) Assuming that the Equation (42) becomes: The variable is the pseudo-range standard deviation. It can be modeled based on the satellite clock and ephemeris, the atmospheric propagation, and the multipath effect. For a single-frequency GPS receiver, a nominal value of the range error is estimated as follows [11]: where is the control-segment standard deviation, is the atmospheric propagation standard deviation, and is the uncertainty due to the multipath effect.

Statistical Method for Rejecting Abnormal Global Positioning System Measurements
One of the assumptions following Kalman filtering theory is that |z r | in Equation (37) and in Equation (40) have normal distributions with zero mean and associated standard deviation and , respectively. However, the output of the INS contains not only the white noise but also systematic errors such as the biases and scale factors. Therefore, the aforementioned statistical assumption is violated. Figure 7 shows the distribution of in a field test. In addition, the uncertainty of the INS and  therefore inno has a normal distribution with zero mean and a standard deviation . Figure 8 shows the histogram of a set of inno calculated from a field test data. Since is determined, the abnormal GPS measurements are detected and eliminated following the rules in Equation (37)  In some cases, only one or two GPS satellites containing range errors may cause drift in the GPS position and velocity. If the LC scheme is used, these aiding measurements are completely rejected, and the integrated system enters the stand-alone mode with the INS mechanization only. In contrast, if the modified TC scheme is used, only abnormal GPS measurements are eliminated while others continue to be used in data fusion for more accurate estimates. This effect is another advantage of the modified TC scheme compared with the LC and original TC schemes. Traditionally, users can set higher cut-off angles for a GPS receiver to reduce the effect of reflected signals in GPS-degraded regions, even with a low-cost GPS receiver. However, such a procedure may remove normal GPS measurements, thus leaving the low-cost TC integrated INS/GPS system unaided, which consequently degrades the overall position accuracy. However, the proposed algorithms determine the quality of the GPS measurements based on the detection of abnormal pseudo-range errors with the aid of the INS mechanization rather than on the elevations of the GPS satellites.

Software Design
Based on the proposed schemes, the software for processing raw measurements from GPS and the IMU was developed in the C++ programming language in the Microsoft Windows operating system. The major functions of the software include reading raw measurements from GPS receivers and various IMUs, processing the data with the EKF and smoothing in real-time, and post-processing with the LC and TC schemes.
For real-time processing, a multi-threaded program is designed to simultaneously implement multiple tasks, including reading data from sensors, real-time processing, logging data for post-processing, and displaying results. Figure 9 shows the architecture of the real-time processing module.
In the data reading engine, IMUs and GPS receivers are supported in the software. Depending on the format of the output data, functions for parsing and decoding the data are designed for a given sensor based on the specified output messages and the protocol specification.
Time synchronization is very important for an integrated system and for data fusion. In general, the incoming data from different sensors have different time counters, formats, and sampling rates. For integration, time synchronization is necessary. The effects and the methods of time synchronization have been previously investigated [8,14]. In this research, a hardware-based method is implemented. In a real-time system, because incoming data from the IMU and the GPS receiver are handled on the same computational platform (PC) at the same time, the system tags the INS data with the incoming GPS time. Because the sampling rate of the IMU is usually higher than that of the GPS, the time stamp for the incoming inertial data is interpolated based on the IMU sampling rate between two incoming GPS measurements. The times of the INS and GPS data are thus synchronized. Figure 10 shows the on-line time synchronization scheme implemented in this study. The time synchronized error depends on the sampling rate accuracy of IMU and GPS. In extra tests, the time synchronized error between Ublox GPS receiver and MIDG-II IMU (Robotics) is about 1 millisecond and between Propak-V3 GPS receiver (Novatel) with C-MIGITS IMU (BEI) is 0.3 millisecond. These errors are not significant to the overall performance of the system.

Experiments and Discussion
Field tests with various scenarios were conducted to verify the proposed strategies. In the first test, two integrated systems were set up to evaluate the effectiveness of the proposed method. The reference system comprised a high-end, tactical-grade IMU, a SPAN-LCI (NovAtel), and a dual-frequency, geodetic-grade GNSS receiver, a SPAN-SE (NovAtel). The specifications for the SPAN-LCI IMU are shown in Table 1. The reference data was generated by the reference system with its raw IMU measurements and GPS carrier-phase measurements processed in the differential mode with commercial software, Inertial Explorer (NovAtel), and sensor fusion was performed in the post-processed TC smoothing mode with an odometer. In general, the accuracy of the reference system was less than 10 cm, which is considered sufficient for this study. Accelerometer bias (milli-g) 1.0 Accelerometer scale factor (ppm) 250 The test system comprised a MEMS tactical-grade IMU, the C-MIGITS (BEI), with an integrated dual-frequency GPS receiver with Doppler measurement provided, Propak-V3 (Novatel). The specifications for the test IMU system are shown in Table 2. SPP technique is used to provide GPS measurements for the test. Both two integrated system were mounted on a mobile mapping van ( Figure 11) for the field tests.
In the first test, the test data sets were collected under various scenarios in urban and suburban areas in Kaohsiung, Taiwan. The test trajectory is shown in Figure 12. Abnormal GPS measurements were included in the test trajectory naturally and by simulation. Figure 13 shows an enlargement of the abnormal GPS measurement scenarios, and Figure 14 shows the number of satellites before and after abnormal GPS measurement rejection.   For the test scenario, three data fusion strategies, namely, pure LC, pure TC, and modified TC with abnormal GPS measurement rejection (TC-M), were implemented. The estimated position, velocity, and attitude were compared to the reference data for further analysis. Table 3 shows the statistics of root mean square errors (RMSEs), and Figure 15 illustrates the performance of the three algorithms in terms of position, velocity, and attitude errors. The analysis in the first test indicates that in good GPS signal environments, the performances of the three integrated strategies are comparable. However, in GPS-degraded environments, there are significant differences. As show in Figure 13, GPS measurements are frequently affected by the multipath effect or an insufficient number of satellites to derive navigation solutions in urban canyons. The GPS solutions are biased to large deviations from the reference trajectory. In these scenarios, the navigation solutions of the pure LC and pure TC schemes are strongly influenced because the abnormal GPS measurements are used to update in the EKF. In contrast, the modified TC scheme can detect and ignore abnormal GPS measurements; therefore, its performance is better than that of the pure LC and TC schemes. However, the navigation accuracy of the modified TC scheme slightly decreases in these scenarios because fewer GPS measurements remain to correct the INS mechanization after the bad measurements have been eliminated. Table 3 and Figure 15 show that, overall, the improvement of the pure TC scheme compared with the LC scheme is not significant, with the horizontal accuracy being worse. In GPS-noisy environments, the pure TC scheme is often worse than the pure LC scheme because of the effect of reflected GPS signals. The performance of the proposed modified TC with interference from abnormal GPS measurements is better than that of the LC and pure TC schemes. The improvement in the RMSE of position, velocity, and attitude is approximately 60%.
In the second test, a real-time integrated system with an LC + NHC algorithm was analyzed. The test system comprised a low-cost MEMS IMU, the MIDG-II (specifications shown in Table 4), and a single-frequency GPS receiver, the AEK-4T (Ublox) in SPP mode. The two devices were connected to a laptop via the RS232 communication protocol, as shown in Figure 16. The system was installed in the mobile mapping van for the field test. The longest GPS outage, 3 min, was generated by driving the car into an underground parking lot. The test results were compared with those of the reference system for performance analysis. Figure 17 shows the performance of the LC + NHC scheme, Figure 18 shows the performance of the pure LC scheme, and Table 5 shows the navigation accuracy of two integration strategies.    The preliminary results from the second test indicate that under open-sky conditions, the position accuracy of the proposed system mainly depends on the accuracy of GPS. The positional accuracy of the integrated system in that case is about 1 to 2 m. In GPS-rejected environments, the navigation solutions of the system completely rely on the output of the INS and aiding sources. With the given MEMS IMU in the second test, in the pure LC mode, the position error grows quickly over time if no GPS updates are available. The position error is approximately 25 m after a 15-s GPS outage. In a long-duration GPS outage (3 min, in this test), the position error drifts to a very large value because of the divergence in the estimator. In contrast, with the augmentation of NHC, for a 3-min GPS outage, the maximum position error is approximately 20 m and the overall position accuracy is approximately 10 m.
Generally speaking, the performance of the integrated system depends on the performance of the IMU and the GPS receiver, the operating environment, and the data fusion strategy. In open-sky conditions, the position accuracy mainly depends on the accuracy of the GPS. Therefore, improving the positioning technique and eliminating the GPS systematic errors are the keys to improving the overall navigation accuracy of the system. In GPS-degraded and -denied environments, the performance of an IMU with a data fusion strategy generally affects the final solutions. Because the goal of this research is to improve the performance of the system using a low-cost MEMS IMU, the data fusion strategy is most critical. As shown in the first test, the overall position accuracy improved by approximately 60% in a GPS-degraded environment by improving the estimation strategy. The position accuracy can reach 10 m with a 3-min GPS outage, as was shown in the second test using a low-cost MEMS IMU with the aid of the NHC in GPS-denied environments.

Conclusions
This study developed a real-time integrated Inertial Navigation System (INS)/Global Positioning System (GPS) navigation system for land vehicle applications. Issues related to GPS uncertainty and the systematic errors of an Inertial Measurement Unit (IMU) were investigated. A modified tightly coupled (TC) integration scheme was used to reject abnormal GPS measurements. Two Micro-Electro-Mechanical Systems (MEMS) IMUs were used for testing under various scenarios.
The preliminary results showed the effectiveness of the proposed strategies. Overall, for the INS/GPS integration mode, the improvement of the modified TC scheme over the standard loosely coupled (LC) and pure TC schemes is approximately 60% in terms of position uncertainty. The main advantage of the proposed system is that it can detect and eliminate abnormal GPS measurements, mainly caused by the multipath effect, thus improving the overall performance of the integrated system. In cases of long-duration GPS signal outages, with the aid of the Non-Holonomic Constraint (NHC), the integrated system using a low-cost MEMS IMU can provide reliable solutions for navigation applications.