A Performance Improvement Method for Low-Cost Land Vehicle GPS/MEMS-INS Attitude Determination

Global positioning system (GPS) technology is well suited for attitude determination. However, in land vehicle application, low-cost single frequency GPS receivers which have low measurement quality are often used, and external factors such as multipath and low satellite visibility in the densely built-up urban environment further degrade the quality of the GPS measurements. Due to the low-quality receivers used and the challenging urban environment, the success rate of the single epoch ambiguity resolution for dynamic attitude determination is usually quite low. In this paper, a micro-electro-mechanical system (MEMS)—inertial navigation system (INS)-aided ambiguity resolution method is proposed to improve the GPS attitude determination performance, which is particularly suitable for land vehicle attitude determination. First, the INS calculated baseline vector is augmented with the GPS carrier phase and code measurements. This improves the ambiguity dilution of precision (ADOP), resulting in better quality of the unconstrained float solution. Second, the undesirable float solutions caused by large measurement errors are further filtered and replaced using the INS-aided ambiguity function method (AFM). The fixed solutions are then obtained by the constrained least squares ambiguity decorrelation (CLAMBDA) algorithm. Finally, the GPS/MEMS-INS integration is realized by the use of a Kalman filter. Theoretical analysis of the ADOP is given and experimental results demonstrate that our proposed method can significantly improve the quality of the float ambiguity solution, leading to high success rate and better accuracy of attitude determination.


Introduction
In low-cost land vehicle navigation applications, not only are accurate position and velocity information required, but the navigation system also needs to provide accurate attitude [1][2][3][4]. Making use of the GPS signal for attitude determination has several advantages such as its small size, low cost, lack of cumulative errors and high accuracy. However, GPS attitude determination also has its own drawbacks. It is susceptible to the external environment and unstable in dynamic applications [4]. The traditional attitude determination method is to use INS. The drawback of INS is that the errors are accumulated over time especially for the low-cost MEMS-INS. If the GPS and MEMS-INS are combined effectively, the reliability and accuracy of attitude determination can be improved [5][6][7].
When a multiple-antenna GPS system is used for attitude determination, the challenge is how to determine the carrier phase integer ambiguity quickly, accurately and reliably. However, due to the loss of lock and noise disturbance, cycle slips often occur in land vehicle applications. Therefore, the single epoch method is usually adopted for dynamic attitude determination in such applications. This is because attitude determination in single epoch is insensitive to cycle slips. In addition, due to the low cost requirement, single frequency GPS receivers are widely used in this type of application, but the carrier phase and code (pseudorange) measurement qualities are both very low for low-cost GPS receivers, especially the code measurement [8]. Moreover, GPS signal is often blocked, attenuated or contaminated by multipath signals in urban areas. As pointed out in [9], these factors contribute to low ambiguity success rate, so the key to high accuracy GPS/MEMS-INS attitude determination is to improve the GPS ambiguity success rate.
The ideas of using inertial sensors to improve GPS ambiguity success rate have been proposed recently [10][11][12][13][14][15]. These methods adopt MEMS inertial sensors to aid different ambiguity resolution methods. They make use of the MEMS-INS attitude information to reduce the integer ambiguity search region, thus improve the rapidity and reliability of ambiguity resolution. For example, the attitude information of rate gyros was employed in [10] to reduce the ambiguity search region to a small cube to improve the ambiguity resolution process which was based on the least squares ambiguity search technique. On the other hand, Zhu et al. [11] applied rate-gyro-constraints to filter the candidates in the ambiguity search stage, which can speed up the initialization of attitude parameters under dynamic circumstances. Eling [13] presented an instantaneous GNSS/MEMS attitude determination method which used AFM aided by MEMS to perform the single epoch ambiguity resolution. Roth et al. [14,15] explored a method that combined low-cost inertial and magnetic field sensors with a GNSS compass to provide a multi-sensor attitude system for portable, small-sized launcher applications. The method can improve the availability and reliability of pure GNSS attitude determination by using an extension of the LAMBDA method which accounts for baseline length and attitude constraints.
LAMBDA is popular for ambiguity resolution, since it is known to maximize the ambiguity success rate [9,16]. It has been used to obtain precise baseline estimations that are then used to extract the platform attitude [17][18][19][20]. However, the LAMBDA method does not include the prior knowledge of the length of the baseline, which is usually known in GPS attitude determination problems [9,21]. Therefore, the CLAMBDA method [22,23] is proposed for attitude determination by integrating the nonlinear baseline constraint into the ambiguity objective function. It further improves the success rate of attitude determination, especially for single frequency single epoch cases [24,25], but the CLAMBDA method is also challenged for the single frequency single epoch case in low-cost land vehicle applications, since the unconstrained float solution is usually of very low quality due to the poor quality GPS measurements, especially when the GPS signal is blocked or contaminated by multipath signals in urban areas, which results in decreased success rates [21,26]. Therefore, for such application, the quality of the float solution should be improved.
In this paper, a new method to improve the attitude determination performance by improving the quality of float ambiguity solution is proposed. First, the INS calculated baseline vector is augmented with GPS carrier phase and code measurements. The quality improvement of the float ambiguity solution is verified by theoretical analysis. Second, we use ADOP combined with the INS derived reference value to assess the quality of the float ambiguity solution. The undesirable float solution is replaced with the float solution obtained through the INS-aided AFM, where the INS attitude information is used to optimize the AFM search region. Moreover, some constraints are adopted for the selection of the correct float solution. Then, the fixed solution is obtained by the CLAMBDA method, and the GPS attitude measurement can be calculated. Finally, the GPS/MEMS-INS integrated filter is used to estimate the navigation errors and sensor errors of MEMS-INS, which can effectively improve the quality of the INS attitude measurement for the estimation of the float ambiguity solution.
The rest of this paper is organized as follows: Section 2 gives an overview of the CLAMBDA method. Section 3 describes our improved attitude determination method for single frequency and single epoch. Section 4 introduces the GPS/MEMS-INS integrated filter algorithm. The experimental results of the new method are shown in Section 5. The final conclusions of this paper are given in Section 6.

GPS Attitude Determination with the CLAMBDA
For the GPS compass system, two antennas are utilized to provide the observability of yaw (or heading) and pitch (or elevation). The GPS compass model with the baseline constraint is presented as [9]: where y is the given GPS data vector, and a and b are the ambiguity vector and baseline vector of order n and 3, respectively. E( )  and D( )  denote the expectation and dispersion operators, respectively.
A and B are the given design matrices that link the data vector to the unknown ambiguity vector and baseline vector, respectively. The variance matrix of y is given by the positive definite matrix y Q .
l is the baseline length, which is a priori given information.
If we apply the least squares (LS) estimation principle to solve for the unknown ambiguity vector a and baseline vector b, we need to solve the minimization problem: y Aa Bb (2) Based on the orthogonal decomposition, the minimization problem of (2) can be formulated as: where ê is the LS residual vector. â Q and ˆ( ) The CLAMBDA method can minimize the following objective function to obtain the fixed ambiguity solution: Equation (4) has no analytic solution, and it should be solved using an efficient searching method. The detailed description of the searching method used here can be found in [23], which based on the "expansion" approach. After we obtain the fixed ambiguity solution, we can also obtain the baseline vector with respect to the local East-North-Up frame, which is used to calculate yaw and pitch.

Improved Attitude Determination Method for Single Frequency and Single Epoch
Our method is based on the orthogonal transformation model, which is a numerically stable approach [27]. For 1 m n   satellites in view, the orthogonal transformation model [28] of single difference carrier phase and pseudorange can be described as follows: 2 1 ,~( ,2 ) the same distribution because orthogonal transformation will not change the statistical properties of white noise. The first step of attitude determination with CLAMBDA is to obtain the float ambiguity solution â (which is the so-called unconstrained float ambiguity solution [21] in this paper) and its vc-matrix â Q by the LS method [29]. â Q contains all the information necessary to infer the quality of the float ambiguity solution. It can be seen that the smaller the â Q , the higher the quality of the float solution. In order to capture the main characteristics of â Q , the ADOP [30] is defined as: ADOP equals the geometric mean of the conditional standard deviations of â . Therefore, ADOP can be used for the quality assessment of the float solution. However, for the single frequency single epoch GPS attitude determination, if the quality of GPS measurements is very low, especially when GPS signal is blocked or contaminated by multipath signals, the float solution is usually of low quality. Although the CLAMBDA method can maximize the ambiguity success rate of attitude determination, the success rate is also not very high due to the low-quality float solution caused by poor quality GPS measurements. Therefore, we augment GPS measurements with MEMS-INS measurement to improve the quality of float solution. This is described in detail as follows.

Float Solution by GPS/INS Augmented Measurements
Assuming the expression of the INS calculated baseline vector in the local East-North-Up frame is , and the expression of baseline vector in the body frame is For the short baseline, the following equation is obtained:  I  I  I  I  I  I  I  I  I  I  I  I  n  b  I  I  I  I  I  I  I  I  I  I  I  I   I  I  I  I where n b C is the INS attitude matrix, which is orthogonal. , ,    are the yaw, pitch and roll of INS, respectively. Assuming the vehicle is a rigid body, b b is unchanged, which can be obtained by accurate measurement. In order to unify the symbols and consider the measurement noise, we adjust Equation (8) as: 2 3 ,~( , ) where ε represents the measurement noise vector of the INS calculated baseline vector in the local frame and 2 I  is the corresponding noise variance. Actually, if the baseline vector is calculated by INS attitude, the entries of the measurement noise vector are correlated. Therefore, the measurement noise variance in Equation (10) is not diagonal. The detailed analysis is as follows: the main baseline vector in the local level frame is expressed as: cos sin cos cos sin When the baseline length is fixed, the baseline vector can be expressed as the nonlinear function of yaw and pitch [31,32]. The nonlinear equations are linearized as:  IE  I  I  I  I  I  IN  I  I  I  I  I  IU  I   I  I  I  I  I  I   I  I  I  I  I  I   I  I where the given Taylor point of expansion is at the INS attitude , I I   . The second and higher order terms are neglected. The expression of the vc-matrix of the measurement noise vector is obtained as:  I  I   I  I   I   I  I   I   EE  I  I  I  I   NN  I  I  I  I   UU  I   EN  NE  I  I  I  I  I  I   EU  UE  I  I  I  To simplify the calculation, the upper bounding approach [33,34] is adopted for the decorrelation of the measurement noise vector. A "more positive definite" diagonal matrix is selected as the upper bound of the vc-matrix of Equation (13). According to the positive definite matrix theory [33], the upper bound of the vc-matrix of Equation (13) should satisfy: Therefore, 2 I  should be large enough to make sure above matrix is a positive definite matrix. It can be seen from Equation (13) . Equation (15) can be expressed as follow: . The state estimate obtained by the LS method is: The vc-matrix of state estimate is: It is easy to prove that the quality of the float solution with the INS calculated baseline vector augmentation is higher than the quality of the float solution without the augmentation. Let: We can obtain the vc-matrix of state estimate: For the float solution without the augmentation, the vc-matrix of state estimate is: It also can be seen from Equations (15), (20) and (21) that more accurate ambiguity float solution can be obtained as long as the number of tracked satellites is not less than two satellites, which is more important for the land vehicle attitude determination.

Quality Assessment of Float Solution
Theoretical analysis results of the ADOP in Section 3.1 demonstrate that better quality of the float solution can be obtained by the GPS/INS augmented measurements. However, low-cost GPS receivers may produce abnormal measurement errors (especially pseudorange errors), which are mainly caused by signal attenuation or multipath in the urban area [35]. This results in part of the float solutions still having large deviations.
In order to reduce the effect of large outliers in pseudorange on the float ambiguity solution, a simple quality control (QC) method is adopted here, which is based on the residual chi-square test [6,10]. As the float solution is calculated by LSQ method, the key issue is how to construct the residual vector. Here, we use the INS calculated baseline vector to calculate the predicted psedorange measurements, and obtain the residual vector as: where the variance of the residual vector is: Then we can perform the residual chi-square test [6,10] on psedorange measurements. In order to detect the large outliers in pseudorange measurements correctly, 2 I  should be conservatively selected to describe the error of the INS calculated baseline vector.
However, undetected measurement errors may still result in some low-quality float ambiguity solutions. In order to achieve high success rate, these low-quality float solutions need to be filtered and improved. So the quality of the float solution estimated by Equation (17) should be assessed first.
Since the attitude errors of MEMS-INS are estimated and corrected by the GPS/MEMS-INS integrated filter in real time, which will be introduced in Section 4, the errors usually are not very large [13,20]. Thus the float solution calculated by MEMS-INS attitude can be regarded as the reference. The error range of the float solution is related with the GPS/INS augmented measurements. So the ADOP calculated by the augmented measurement equations can be used for the quantitative description of the error range. Therefore, the range of the float solution can be determined by the MEMS-INS attitude and the ADOP.
First, the INS attitude calculated baseline vector bI in the local frame can be obtained by Equation (8). Then according to the carrier phase measurement Equation (5), the INS calculated float ambiguity solution ˆI a can be solved as follows: Second, the ADOP can be calculated with Equations (7) and (18). Then, according to the 3 principle, the float solution lies within the 3 bound with a high probability, thus the range of each entry of DD float ambiguity solution â is chosen as: If one entry of the float solution estimated by Equation (17) does not satisfy Equation (24), the float solution is marked as low-quality.

Float Solution with INS-Aided AFM
For the low-quality float solution filtered by Equation (24), we will replace it with the float solution resolved by the INS-aided AFM. The basic idea of AFM is that for the correct attitude, the value of the adaptive function should theoretically have the value 1 as the maximum, since the ambiguity is integer. The maximum of the adaptive function ideally leads to the correct ambiguity [36]. AFM does not use pseudorange measurement to calculate the float ambiguity solution, so it rejects the influence of bad pseudorange, which is an advantage of the AFM method. MEMS-INS attitude can be used to determine approximate attitude to reduce the size of the search space [13]. It can improve the computational efficiency and reliability of the ambiguity resolution. When searching in the space constrained by INS attitude, the correct ambiguity is always contained, so a high-quality float solution can be obtained.
Given the pitch k  , yaw k  and baseline length l , the main baseline vector bk of the vehicle can be calculated by Equation (11), then the float solution ˆk a can be calculated by Equation (23).
Accordingly, the adaptive function of the float solution can be expressed as: The search space is determined by MEMS-INS attitude. , However, despite the search space reduction, several candidates of similar magnitudes may still exist in the adaptive function. Therefore, in order to identify the most likely correct candidate, we need the validation procedure with some constrain conditions. Before the validation, for each float solution candidate, the fixed solution is obtained through rounding. The corresponding baseline vector ˆk b can be estimated by LS method, then the attitude ˆ( , ) k k   can be calculated.
The validation procedure is divided into three consecutive steps: (1) Baseline length verification For the correct float solution candidate, the error between the estimated and known baseline length is very small, so if the estimated baseline length satisfies , the candidate will be sent to next step. When the fixed solution and the estimated baseline vector of each candidate are substituted back to Equation (5), the corresponding residual can be obtained. The correct candidate will make the residual reach the minimum. Thus the remaining candidates are sorted according to the ascending order of their residuals. The float solution of the maximum value of adaptive function is selected from the first k candidates. The k can be set according to the number of remaining candidates. If the validation procedure returns empty, increase the l  , and repeat the whole validation procedure from step (1) to step (3).

The Implementation of the Proposed Method
The flow diagram of the method is shown in Figure 1. First, the INS calculated baseline vector is augmented with the GPS carrier phase and code measurements. Second, according to the result of quality assessment, the undesirable ambiguity float solution is replaced with the float solution obtained through INS aided AFM, where the INS attitude is used to reduce the AFM search region. Then, the fixed solution is obtained by the CLAMBDA method, which is used to calculate GPS attitude measurement. Finally, in order to effectively control the quality of the INS attitude measurement for the estimation of the float ambiguity solution, the GPS/INS integrated filter is used to estimate the navigation errors and sensor errors of MEMS-INS, which will be introduced in the next section. The proposed method can significantly improve the quality of the float ambiguity solution. It can not only improve the success rate of attitude determination, but also the computational efficiency. Since when the float solution is accurate, the search time of the fixed solution by CLAMBDA is always very short.

GPS/MEMS-INS Integrated Filter
After the attitude determination using MEMS-INS aided multiple-antenna GPS, GPS/MEMS-INS integration is also necessary. The integrated Kalman filter estimates the navigation errors and sensor errors using the GPS pseudorange, pseudorange rate and attitude measurements.

State Equations of the Integrated Filter
The states of integrated filter are given by: 17 where t R , m R and 0 R are the tangential radius, the meridional radius and the mean radius of the Earth, respectively; the symbols 0 g and  are local gravity and the Earth's rotation rate, respectively. The gyro, accelerometer, and the receiver frequency errors are modeled as first-order Markov processes, and   , f  and tr  are the corresponding correlation coefficients.

Measurement Equations of the Integrated Filter
The measurements of the integrated filter are the difference of pseudorange, pseudorange rate and attitude between GPS and MEMS-INS, which are given by:  I  G  I  G  I  G  I  G  I  G  a  I  G   I  G  Im Gm Im Gm The measurement equation of pseudorange is described by [37]:  cos cos cos sin sin The measurement equation of attitude is described by: Then the measurement equations of the integrated filter can be obtained: Based on the discrete forms of Equations (27) and (32), we can implement Kalman filter algorithm to estimate the MEMS-INS errors in real time. The Kalman filter algorithm is composed of prediction and measurement update [6]. Using the error estimates of the integrated filter to correct MEMS-INS errors, we can obtain the final attitude determination results. After feedback correction to the MEMS-INS navigation processing, MEMS-INS can provide better attitude measurement for the estimation of the float ambiguity solution.

Experiment Test Results
Field tests in an urban area were conducted to verify the performance of the GPS/MEMS-INS integrated attitude determination system. The test system is shown in Figure 2. It consists of a MEMS-INS and the GPS attitude determination system developed by our research group. The GPS attitude determination system includes three receivers, an interface unit, a navigation processing unit, and a display unit. The interface unit includes a serial port for MEMS-INS. The system can work in real-time processing mode or data collection mode.
The receivers used here are the SUPERSTAR II, which is the NovAtel 12 channels single frequency GPS receiver with 1 Hz output. Its code measurement precision is 0.75 m RMS. The difference carrier phase measurement precision is 0.01 m RMS. The MEMS-INS is the XW-IMU5220 from Beijing Starneto Technology Co. Ltd. (Beijing, China). The bias stability of gyro and accelerometer are 0.02°/s and 8 mg, respectively. The output frequency is 100 Hz.

Test Settings
The experimental setup of the attitude determination system is shown in Figure 3. Two GPS baselines were mounted on the roof of the car. The main baseline which is used to determine the vehicle pitch and yaw comprises two antennas aligned in the driving direction. The auxiliary baseline which is used to determine the roll is at right-angled to the main baseline. The main baseline length is 1.0 m and the auxiliary baseline length is 0.9 m. The attitude determination system was mounted in the car trunk, and the measurement axes of MEMS-INS were corresponded with the baseline directions. The experimental site is on the east side of DaTun Road, Beijing Olympic Park, which is shown in Figure 4. The car moves along a long and narrow rectangular block for about five laps and both ends of the rectangle block are arc-shaped. The experimental time is about ten minutes. The distribution, number and PDOP of actual visible GPS satellites are shown in Figure 5. The number of visible satellites changes frequently due to the blockage of surrounding buildings.  In order to further test the performance of the attitude determination system, three satellites are removed to create a GPS challenged environment, i.e., three satellites which have the lowest elevation angles are regarded as invisible. Figure 6 shows the distribution, number and PDOP of the visible GPS satellites after three satellites with lowest elevation angles removed. In this GPS challenged environment, the scenario with less than four visible satellites frequently occurs.

Test Results
The results of the proposed MEMS-INS aided GPS attitude determination method are compared with the unaided method. The performance improvements are verified through examing the quality of float solution, success rate and attitude accuracy. In the following, the default environment is the actual GPS environment shown in Figure 5, and the default baseline is the main baseline. It can be seen that the quality of the float solution is improved (i.e., smaller ADOP) when the INS calculated baseline vector is augmented with the GPS carrier phase and code measurements. Since the precisions of the INS calculated baseline vector and GPS carrier phase are higher than that of the pseudorange, the contribution of GPS pseudorange is decreased after using the INS calculated baseline vector augmentation. The improvement is especially significant when additional three satellites with lowest elevation angles are removed. In other words, with the INS calculated baseline vector augmentation, the ADOP is less affected by the number of visible satellites. This can be seen clearly in Figure 8. It demonstrates that high-quality float solution can be obtained by using GPS/INS augmented measurements, even if the number of visible satellites is less than four. Since the float solution cannot be calculated using only GPS measurements when the number of visible satellites is less than four, the ADOP value is set to 45 for this situation to draw the figures. We next investigate the float solutions of our proposed method. In the test runs, satellite No. 2 and No. 5 were always visible. Satellite No. 2 had the highest elevation angle and it was selected as the reference satellite. Figure 9 shows the DD ambiguity float solution of satellite No. 5. From Figure 10b, it can been seen that 100% of the replaced float solutions are inside the range determined by Equation (24). Figure 11 shows the attitude determination results corresponding to the float solutions with and without the low-quality float solutions replacement. It clearly indicates that when the low-quality float solutions are used, the calculated attitudes are always wrong. The quality improvement of float solution is an important factor for successful attitude determination with CLAMBDA. In the experiment, since the car moves on a flat road, the pitches of the two baselines can be considered constant in a certain allowable error range, here we use their initial values as the references and set 5° as the error range. The comparison of the attitudes calculated by the two attitude determination methods is shown in Figure 12.
From Figure 12, it can be seen that the attitude calculated by the MEMS-INS aided GPS method is more accurate and stable than that of the unaided GPS method. The yaw calculated by the MEMS-INS aided GPS method clearly shows the movement of the car. The lines approximately vertical to the angle axis represent the movement along a straight line, the two changing parts of the curve before 150 s show the starting and 90° turn to the rectangular block, and the dramatically changing parts of the curve after 150 s show the 180° turn of the car for the 5 laps. However, it can not be obtained by the attitude calculated by the unaided GPS method, because unsuccessful attitude determination often occurs in the unaided single frequency single epoch case, which results in large attitude error. Figure 13 shows the coordinates of baseline in the local frame of the MEMS-INS aided GPS method, where N, E, and U denote the north, east and up coordinates of baseline, respectively; L denotes the length of baseline. The maximum error of the baseline length is 0.05 m, which indicates that accurate baseline coordinates can also be obtained by the aided method. Then, the success rates of the MEMS-INS aided GPS and the unaided GPS attitude determination methods are given in Table 1. The results can be expected after the float solution quality verification.
Since external factors such as multipath and low satellite visibility in urban environment further degrade the measurement quality of the GPS receivers, the success rate of the unaided GPS attitude determination is not very high due to the low-quality float solution. The success rate of the proposed MEMS-INS aided GPS method is above 98%, which is much higher than that of the unaided method.   Figures 14 and 15 further show the influence of the number of visible satellites on the attitudes calculated by the two attitude determination methods, respectively. When the number of visible satellites is less than four, the unaided GPS method cannot determine the attitude, so the pitch and roll are set to 90°, and yaw is set to 0° for drawing the figures. The corresponding success rates are listed in Table 2.
As shown in Figures 14 and 15 and Table 2, the success rate of the unaided GPS method obviously decreases to below 50% after the three satellites with lowest elevation angles are removed, but for the MEMS-INS aided GPS method, the success rate slightly decreases about 1%, which in general is still very high (above 97%). Its advantage is very evident, since the attitude determination is also successful at times of two visible satellites.   From the above results, it can be seen that the success rate of GPS/MEMS-INS attitude determination is usually unable to achieve 100%, since the maximum success rate depends on the quality of carrier phase measurement and the performance of MEMS-INS in the land vehicle application. When the integer ambiguity resolution is unsuccessful, correct attitude can not be obtained, and the attitude error may be large. It is no doubt that if the abnormal GPS attitude measurement is used in the integrated filter, the large error will degrade the overall attitude accuracy of the integrated system. To avoid this problem, a simple quality control (QC) method is adopted in the GPS/MEMS-INS integration, which is based on the residual chi-square test [6,10].  Figure 16, whose initial value is assigned by GPS attitude determination. The standard deviations of pitch and roll are shown in Table 3.  As compared with Figure 12, Figure 16 demonstrates that the attitude accuracy is improved through the integration of GPS and MEMS-INS. It also demonstrates that the attitude results of the aided GPS/MEMS-INS integration are much better than that of the unaided integration. Since the success rate of MEMS-INS aided GPS attitude determination is much higher, the attitude errors of MEMS-INS can be corrected effectively by the GPS/MEMS-INS integration. It can be seen from Table 3 that the attitude accuracy improvement by the aided GPS/MEMS-INS integration is significant, especially for a GPS challenged environment.

Conclusions
In order to improve the performance of the low-cost GPS/MEMS-INS attitude determination system used in land vehicles, a new integrated attitude determination method is presented. The core issue of the method is how to improve the success rate of the single frequency single epoch GPS ambiguity resolution with low-quality measurements. We adopt the GPS/INS measurements augmentation, the quality assessment of float solution and the undesirable float solution replacement by INS aided AFM to improve the quality of the float solution. Then, the CLAMBDA method is used to obtain the fixed solution, which can maximize the ambiguity success rate for attitude determination. Finally, the GPS/MEMS-INS integrated filter is designed to increase the attitude accuracy. Field test results in urban area demonstrate that our proposed method can significantly improve the quality of the float ambiguity solution, the success rate, and the accuracy of attitude determination especially for GPS challenged environment. The success rate increases to above 97%. Compared with the unaided GPS/MEMS-INS integration, the attitude accuracy is improved at least 35%. The improved low-cost GPS/MEMS-INS attitude determination system can offer a superior performance and efficiently fulfill the task in the land vehicle application.