1. Introduction
With the rapid development of fifth-generation (5G) communication, autonomous driving, the Internet of Things (IoT) and other information technologies, the era of the Internet of Everything is coming. The demands for real-time precise and reliable relative positioning services in multiple moving equipment are increasing, especially in application areas requiring short distances and high speeds, such as the formation of aircrafts, vehicles, and ships. This brings a higher challenge to real-time high-precision positioning services. Multiple moving equipment usually form a joint platform system, and their data are exchanged to realize self-organized movement, i.e., collaborative work like bees. The accurate and reliable real-time position and attitude information of each device are the key for the safe and effective running of the moving platform system. Besides, to ensure a safe distance and cooperative operation between two devices, only the synchronous precise relative position and attitude information are sufficient. Although precise point positioning (PPP) technology can obtain centimeter-level absolute positioning results with only one receiver, the long initialization time and the dependency on additional network communications and augmentation products in real-time processing limit its application [
1,
2,
3]. However, real-time kinematic (RTK) positioning technology, which utilizes differential information from a reference station, can realize fast and easily implemented real-time centimeter-level relative positioning; thus, it is widely employed in numerous real-time high-precision applications [
4].
According to the move status of the reference station, the RTK technique can be divided into two modes, i.e., the “static-to-kinematic” mode, where a static reference station with known coordinates is required to transmit the differential data to moving devices, and the “kinematic-to-kinematic” mode, where both stations at the baseline are moving for real-time relative positioning. In the majority of contributions on kinematic positioning for a joint moving platform system, the “static-to-kinematic” mode is applied, and each moving device acquires centimeter-level positioning results relative to the nearby static reference station once the integer ambiguity is correctly fixed. Through the coordinate transfer, the precise relative position between any two moving devices in this platform system can be determined. However, this strategy relies heavily on the static reference station. The unexpected interruption of the communication link or data processing in the static reference station will paralyze the entire moving platform system. In the meantime, since the reference station is static, the coverage area of the sent differential data is limited, and the performance of the ambiguity initialization and positioning accuracy will decrease if the moving platform system is away from the static reference station. Furthermore, the real-time relative positioning of each device in the moving platform system requires large data transmission and communication ability in the static reference station if there are lots of moving devices, which means a higher cost for equipment at the static reference station. Alternatively, the other “kinematic-to-kinematic” mode has been gradually attracting more and more attention [
5].
In the “kinematic-to-kinematic” mode, any moving device in the platform system can be selected as the reference station for real-time positioning relative to the other moving devices. In open-source software RTKLIB version 2.4.2, Takasu and Yasuda realized the kinematic-to-kinematic relative positioning in the module’s “moving baseline” [
6]. Single-point positioning (SPP) was used first to obtain the approximate coordinates of the reference moving station; then, Kalman filtering was applied to estimate the position, velocity, and single-differenced ambiguity, and finally, the baseline solution was achieved once the ambiguity was fixed [
7]. However, the actual performance of the positioning was not inspiring due to the unsatisfactory success rate of the ambiguity resolution.
It should be noted that, although this software is very popular and lots of the literature on kinematic data processing with RTKLIB can be found, unfortunately, most of them only utilized the “PPP kinematic” or “static-to-kinematic” modes, and works in the “kinematic-to-kinematic” mode can hardly be found. For example, Berber et al. analyzed the results from other open-source online services together with RTKLIB in the PPP mode [
8], Gurturk et al. analyzed the results of flight data from RTKLIB in both the PPP with ambiguity resolution and postprocessing kinematic methods [
9], and Wu et al. analyzed the positioning performance of smartphone data based on RTKLIB with, also, the PPP algorithm [
10]. Hemant et al. utilized RTKLIB, along with two stations for kinematic positioning, but the base station was static; thus, the “static-to-kinematic” mode was applied [
11].
Besides the results from RTKLIB, Luo constructed the optimal Delaunay triangulation in the platform system and selected the optimal baseline as the constraint condition to accelerate the ambiguity convergence [
12]. Based on this approach, MultiKIN software was developed for relative positioning in a joint moving platform system. Dong et al. (2020) considered a data delay during the movement of a reference station and presented an asynchronous RTK (ARTK)/time-differenced carrier phase (TDCP) and inertial navigation system (INS)-integrated method, and a decimeter-to-centimeter-level higher-rate (≥100 Hz) relative positioning was obtained [
13].
It was obvious that all the estimated coordinates of the moving devices during the whole observation were compatible and consistent in the “static and kinematic” mode, since all their coordinates referred to the only static reference station with a fixed coordinate, which means the position datum of the entire moving platform system was fixed to the static reference station and would not change over time. In this way, all devices in the entire moving platform system could achieve their consecutive trajectory by epochs [
14,
15]. However, in the “kinematic-to-kinematic” mode, a simple epoch-independent SPP cannot reflect the precise position change of the moving reference station due to the accuracy limitation of the SPP. Under these circumstances, although centimeter-level relative positioning can be reached between the moving reference station and the moving rover station in each epoch, their absolute coordinates over time will have frequent jumps along the true moving direction, which is caused by the positioning errors of the SPP. Even if Kalman filtering is applied to smooth the coordinates, the accuracy of the positioning results achieved by the pseudo-range-only positioning is still unsatisfactory [
16]. In other words, the position datum in the moving platform system is very unstable during the entire movement, which leads to an inconsistency between the positioning results and the actual motion trajectory of each moving device in the platform system. Although this problem can be solved by the precise real-time PPP, the longtime ambiguity initialization and the requirements of the additional GPS/BDS state space representation (SSR) corrections described above make it unpractical to apply this technique in “kinematic-to-kinematic” positioning [
17,
18].
Considering that the position variations of the moving reference station can be accurately determined if observables with high sampling rates (≥1 Hz) are processed, the stability of the position datum, which refers to the accuracy of the position variations of the moving reference station, can be maintained. Usually, three different methods can obtain the position variations, including position difference (PD), Doppler measurements, and TDCP [
19]. The PD method differences two consecutive position solutions and is very easy to implement, but the accuracy of the determined position variations relies much on the position accuracy, and the change of the satellite number will introduce a discontinuity of the position. The Doppler measurements produced by the receiver satellite relative motion enable an average velocity accuracy of around 2 cm/s, and the position variations can only be calculated by multiplying the epoch interval. However, the velocity of the moving device may vary during the epoch interval, and the product cannot accurately reflect the actual position variation. Therefore, when the moving devices accelerate rapidly, the accuracy cannot be guaranteed [
20]. On the contrary, with the TDCP method, differences in the carrier phase observables between consecutive epochs can directly lead to position variations of the moving devices [
21]. Meanwhile, for data with a sampling rate no lower than 1Hz, the TDCP method almost removes the time-correlated errors, such as the clock and ephemeris errors in the satellite segment and the atmospheric delays, as well as eliminates the integer ambiguity [
22]. In this way, centimeter-level position variations can be achieved [
23]. As a result, the precise positions variation obtained by the TDCP method in kinematic relative positioning can be used to transfer the coordinates of the moving reference station along with the time and, thus, ensure the stability of the position datum.
On the other hand, we should realize that, in environments with severe signal obstructions such as urban canyon areas, the accuracy, reliability, and continuity of Global Navigation Satellite System (GNSS)-only real-time relative positioning cannot be guaranteed. Usually, the GNSS system is integrated with the inertial navigation system (INS) to improve the navigation and positioning performance [
24,
25]. Besides, considering the cost of the traditional inertial measurement unit (IMU) equipment, the low-cost micro-electro-mechanical system (MEMS) becomes an ideal alternative in the civil consumer market due to its small size, high rate, and good dynamic performance, and it is more and more popular in precise kinematic positioning applications such as unmanned aerial vehicles (UAVs), automatic driving, etc. [
26,
27,
28].
In view of the current instability of position datum in “kinematic-to-kinematic” relative positioning in a platform system with multiple moving devices, in this paper, we will propose an integrated method in which a loose combination (LC) between the GPS/BDS TDCP method and the INS is employed to maintain the stability of the position datum in the moving platform system. The consistency between the positioning results and the actual motion trajectory further ensures the safety distance and collaborative operation between the moving devices. The performance of the proposed method will be analyzed and compared with the tight combined (TC) RTK/INS results [
29,
30].
2. Materials and Methods
2.1. SPP
In a simple GPS or BDS SPP model, the code observation equation of any satellite is expressed as:
where
refers to the code observation.
is the geometric distance between the satellite and the receiver.
c is the speed of light.
and
are the clock error and the code hardware delay, whose subscript
r and superscript
s are the receiver and the satellite, respectively.
Orb denotes the orbital error, and
and
are the ionospheric delay and the tropospheric delay, respectively.
indicates the multipath effect and the measurement noise.
For at least four tracked satellites, the linearized observation equation in the matrix could be written as follows:
where
P represents the observation vector of the code,
A is the coefficient matrix, and
A = (
B,
E), where
B is the direction matrix from the receiver to the satellites, and
E is a column vector with all elements being 1. The parameters to be estimated
X include the coordinate corrections and the clock error of the receiver, i.e.,
. Herein, the hardware delays of the receiver and the satellite are absorbed into the receiver clock and cannot be separated. The remaining matrix
L is the calculated ranges between the receiver and the satellites, which include the initial geometry distances between the receiver and the satellites, the modeled satellite clock errors, and corrections of the ionospheric delays and tropospheric delays by empirical models such as the Klobuchar model [
31] and the Saastamoinen model [
32], respectively.
Furthermore, in a combined GPS/BDS model, the intersystem code bias (ISCB) in different systems should be considered, and the estimators in Equation (2) are modified to
where the superscripts G and C denote the GPS system and BDS system, respectively. As a result, the least-squares (LS) estimation can only be employed for more than five tracked satellites, and the specific observation equation can be expressed as
As for the stochastic model, all observations are regarded as independent, and their variances are determined by a simple satellite elevation-related model, expressed as
where
is the a priori code noise of 1 m for both GPS and BDS, and
E is the satellite elevation. Besides, the weights of different systems are determined separately and formed as a diagonal matrix. Therefore, the weight matrix of the combined GPS/BDS system
PW is expressed as
2.2. Kalman-Filtered SPP (KFSPP)
For a quick-moving device, the only position is not enough to fully describe its motion state, and the velocity is an important supplement. By acquiring the velocity information, one can establish an approximate kinematic equation for the moving device, and its state of motion can be precisely depicted. For this reason, we can estimate the velocity of the moving device by using Doppler measurements.
The user satellite relative motion leads to changes in the phase and frequency, or the Doppler shift. Scaling the Doppler shift by the carrier wavelength, the code changing rate observation can be obtained, whose linearized form can be written as follows:
where
is the wavelength of the carrier phase,
is the Doppler shift value,
is the coefficient matrix,
and
are the velocities of the receiver and the satellite, and
is the receiver clock drift.
Since the motion state of the moving device can be described, the Kalman filtering (KF) technique can then be used to determine its position and velocity together with the observation equations of the code and code changing rate. Provided that the velocity of the moving device is constant between consecutive epochs, its motion state can then be expressed as
where
k represents the epoch.
is the vector of the device motion state, and
.
is the transition matrix between the consecutive epochs. According to the kinematical equation for a moving object with a constant velocity, the transition matrix can be expressed as
where
I is a 5 × 5 identity matrix, and
is the time interval.
is the process noise vector whose covariance matrix is determined as follows:
where
is a 3 × 3 diagonal matrix representing the power spectral density (PSD) of the acceleration, and
and
are 2 × 2 diagonal matrices denoting the PSD of the clock bias and clock drift. Herein, all the diagonal elements in
, , and
are set to an empirical value of 0.05. Besides, the variance of the code and the Doppler observations are set to 4.5 and 0.18, respectively.
Using at least five simultaneous code and Doppler measurements, the position, velocity, clock offset, and clock drift of the moving device could be obtained by the KF technique. Compared with the simple SPP technique, Doppler observations are involved, and Kalman filtering is employed in KFSPP; thus, better accuracy of the positioning is expected.
2.3. TDCP
In the above KFSPP, the velocity is determined, and the kinematic equation is utilized to improve the accuracy of the position. However, the assumption that the velocity of the moving device is constant between consecutive epochs is probably unrealistic, and the velocity of the device may vary from one epoch to another. Alternatively, the TDCP observation can provide another way to directly determine the position variations between adjacent epochs.
The TDCP observation is derived from the original phase observation, which can also provide receiver satellite distance estimations, and the original phase observation equation in meters can be expressed as [
33]:
where λ is the wavelength,
is the carrier phase in cycles,
is the phase hardware delay,
N is the integer ambiguity, and
includes the multipath and the phase noise. The TDCP observation is then obtained by differencing the carrier phase observations of the same satellite at two successive epochs:
t and
t+1, which can be expressed as:
where
represents the time-differencing operation. Assuming that no cycle slips occur between these two epochs, in the TDCP observation, the unknown integer ambiguity is totally canceled, because it is a constant. Most of the other common errors such as
,
,
,
,
and
vary slowly with the time and can also be almost eliminated in a very short time interval [
34]. Therefore, the TDCP observation equation can be simplified to
Just like Equation (2), the linearized TDCP observation equation in the matrix can be obtained if at least 5 satellites are continuously tracked at two consecutive epochs for the combined GPS/BDS system. Then, the LS estimation can be applied to obtain the position variations
and the clock drift of the moving device. The position of the moving device at the current epoch can be easily derived from the position at the previous epoch as
Since carrier-phase observations with millimeter-level noise are used, there is a promising probability of centimeter-level position variations being obtained. In the meantime, the carrier phase also brings possible cycle slips. To handle this problem, the TurboEdit method [
35] was first utilized for data preprocessing. Then, the posterior mean error of the unit weight was checked after a standard least-square estimation. If this value is remarkably increased at one epoch, it means the cycle slip occurs in this epoch, and a robust least-square estimation was thereby applied to detect the specific satellites with cycle slips. At last, additional parameters were added in the observation equation to resolve these cycle slips. In this way, the possible cycle slip is detected and repaired.
2.4. TDCP/INS
Despite the high precision of carrier-phase observations in the TDCP method, the signal is susceptible to the environment around the station. In order to ensure a high reliability and availability of the moving reference position, in this study, the INS data will be collected and integrated with the GPS/BDS data, and an improved integrated method will be proposed. Since this method combines the results from the TDCP method and the INS data, we call this method the TDCP/INS method.
In the TDCP/INS method, the precise position variations are determined via the TDCP method, and the position at the current epoch is updated through Equation (13). Once the updated position is achieved, a loose combination integration (LC) between the position information and the velocity/acceleration information from the INS data can be employed. In this way, the stability and reliability of the position results in complex urban environments can be further improved.
In the INS error model, only the stochastic biases of the gyros and accelerometers are to be considered and are modeled as the first-order Gauss-Markov process. As a result, a 15-state INS dynamic model, including the three-dimensional attitude
, velocity
, position
, gyros bias
and acceleration bias
is established. The form of its state equation is the same as that in Equation (7), but the parameter to be estimated is modified to
. The definitions of the transition matrix
and the process noise vector
are referred to in Reference [
23], where the variances of the position in three different directions were all set as 9 cm
2, considering the centimeter-level accuracy of the position variations.
The observation equation of the TDCP/INS model is constructed based on the updated position results determined by the combined GPS/BDS system and the corresponding information predicted by INS data and can be expressed as:
where
is the observation vector, expressed as
where
and
represent the position and velocity vectors obtained by GPS/BDS and INS, respectively.
is the coefficient matrix, and
is the observation noise vectors. Based on the established state equation and the observation equation, the KF technique is then applied in integrated GPS/BDS/INS positioning to obtain the attitude, velocity, and position information.
Compared with the TDCP alone method, the INS data further ensures the reliability of the position results in complex urban environments. In this way, the stability of the position datum can be maintained.
2.5. Method Comparison
To demonstrate and compare the effects of the positioning results from different methods, we drew a sketch map showing the position information of one moving device along with the time, as seen in
Figure 1. The positions obtained from different methods and the corresponding trajectories are displayed with dots and lines in different colors.
Specifically, the moving reference position at the beginning epoch t can only be obtained by the SPP method; thus, there is inevitably an obvious deviation between the calculated position and the true position caused by the positioning error of the SPP. This situation occurs for all methods, including the SPP, KFSPP, TDCP, and TDCP/INS.
In the subsequent epochs, the position variations are obtained by the TDCP alone or the TDCP/INS method, and the moving position at the next epoch is updated from Equation (13) and further improved through the integration with the INS data. Since the determined position variations are accurate and the improved position is reliable, the final positions will continuously keep almost the same bias from the true ones. Therefore, the trajectory from the TDCP alone or TDCP/INS method is approximately a translation of the actual trajectory, as is shown with the green dots and line in
Figure 1.
On the contrary, the positions obtained by the SPP method have random errors of 3–10 m. Even if the KF technique is applied in the KFSPP method, the random errors are hardly reduced to below 1 m. Therefore, the position results obtained by these two methods show frequently jumps around the true positions at all the epochs, as we can see with the blue dots and line in
Figure 1.
Although the trajectory by the TDCP alone or TDCP/INS method has a systematic bias from the actual trajectory, the bias is very stable along the timeline. For a platform system with multiple moving devices, the constant systematic bias does not matter, and the stability of the position datum can be maintained during the whole motion period.
3. Results
In this section, a vehicle field experiment is introduced when the GPS/BDS and INS data were collected. The satellite availability of the GPS and BDS and the position dilution of precision (PDOP) of the moving device during the observation are analyzed to reflect the observation environment. Then, a quantitative analysis is made, and the numerical results of the position variations acquired from the SPP, KFSPP, TDCP alone, and TDCP/INS are evaluated. Finally, a qualitative analysis is made to visually assess and compare the consistency of the motion trajectory obtained from different methods in three typical urban driving scenarios, which include going along a large, curved road, turning at a crossroad, and passing under an urban viaduct.
3.1. Data Collection
In order to analyze and assess the performance of the TDCP/INS method, nearly one-hour’s vehicle field experiment was carried out on November 30th, 2020 in an urban environment of Wuhan, China. A K708 GNSS original equipment manufacturer (OEM) board manufactured by ComNav Technology Ltd. in China was used to collect the GPS/BDS data, and an ADIS16470 MEMS module manufactured by Analog Devices Inc. was used to synchronously collect the IMU data. The specifications of these two products and their configurations in use are listed in
Table 1 and
Table 2, respectively. Since only the low-cost MEMS is applied, the noise of its observations is relatively larger than other expensive INS modules. On the other hand, no blunder detection or noise elimination was employed in this experiment, which may have limited the accuracy of the final results.
It should be noted that the sampling rate of the GNSS receiver was set to 1 Hz, while that of the MEMS module was set to 100 Hz, and only the results at synchronous epochs will be compared for the three methods. More details about the models and strategies applied in the data processing of the TDCP/INS method are shown in
Table 3.
The GNSS OEM board was connected with an antenna via the cables, and then, it was integrated with the MEMS module. Two sets of the integrated devices were placed on the roof rack of a sport utility vehicle. One was used in this study, and the other two were used for another purpose that was outside the scope of this study. The experimental vehicle and the installed equipment are displayed in
Figure 2.
To obtain accurate reference positions of the moving vehicle, one static reference station was set up in the Science Park of Wuhan University. The high-precision vehicle positioning and attitude determination system M39 developed by Wuhan MAP Space Time Navigation Technology Co., LTD was used, and the centimeter-level TC RTK/INS results were regarded as the “true” positions and velocities [
36]. It should be noted that the proposed loosely combined TDCP/INS method belongs to the position domain and can directly obtain absolute positions without the information of the base station, while the tightly combined RTK/INS method belongs to the observation domain, and the achieved relative positions can only be changed to the absolute ones by the known coordinates of the base station.
The actual complete driving trajectory of the vehicle experiment from the TC RTK/INS method is shown in
Figure 3, and those from the SPP, KFSPP, TDCP, and TDCP/INS methods are shown in
Figure 4. From the actual trajectory, we can see that the vehicle turned at several corners at the beginning and then went straight along a long road. Additionally, the trajectories from the four other methods were similar to the actual one; this was due to the limited plotting scale on the maps.
Meanwhile, the vehicle velocities from TC RTK/INS in the east (VE), north (VN), and upward (VU) directions along the timeline are shown in
Figure 5. The velocity results also showed various scenarios of motion states during the nearly one-hour field test, which included frequent accelerations and decelerations.
3.2. Satellite Availability and PDOP
The satellite availability and the PDOP value are two important indicators when we evaluate the service performance of the system, and they have a significant impact on the final positioning accuracy. If more satellites are tracked and smaller PDOP values are achieved, the positioning results will be more accurate.
The tracked GPS, BDS-2, and BDS-3 satellites during the nearly one-hour field test are shown in
Figure 6 when the cut-off elevation angle was set to 15°. Generally speaking, the observation conditions of the test were pretty good, because approximately 8 GPS satellites, 11 BDS-2 satellites, and 10 BDS-3 satellites could be tracked. On the other hand, we also found frequent signal interruptions for several satellites. This may be due to the changes or the interferences of the environment around the vehicle, such as the tall trees beside the road or the viaduct on the road. In these situations, the number of tracked satellites rapidly decreases, and the positioning performance probably gets worse.
Furthermore, we calculated the number of visible satellites for the single and combined systems and their PDOP values at every epoch, as displayed in
Figure 7. It can be seen that, compared with the only GPS system, the number of available satellites significantly increased after combining with the BDS system, and the PDOP values decreased accordingly. In the meantime, we also found that the sudden decrease of the satellite numbers in several epochs immediately caused sharp increases of the PDOP values. This implies that the satellite number and the PDOP value were negatively correlated.
The subsequent statistical results showed that the average numbers of the GPS, BDS, and GPS/BDS satellites were 7.5, 17.1, and 23.8, and the corresponding average PDOP values were 2.21, 1.71, and 1.36, respectively. Compared with a single GPS, the average PDOP values of the BDS and the GPS+BDS combined system increased by 22.6% and 38.5%. That is the reason why we selected the combined GPS/BDS system for data processing.
3.3. Evaluation of Position Variation
In order to analyze and assess the performance of the TDCP/INS method in the stability of the position datum, the position results from TC RTK/INS were used as the reference. Meanwhile, the position variations between adjacent epochs were also obtained in the SPP, KFSPP, and TDCP methods for comparison. The position variations along the timeline in the SPP, KFSPP, TDCP, TDCP/INS, and the referenced TC RTK/INS methods are displayed in
Figure 8.
From
Figure 8, we can see that the position variations obtained from the TC RTK/INS method ranged within 15 m in the east/north direction and 1 m in the upward direction, and these results were regarded as the reference values. The results obtained from the TDCP/INS method were much similar compared to the reference, which means that the TDCP/INS method also accurately described the actual motion state of the moving devices. Comparatively, although the TDCP-only method achieved similar position variations along with the time, there were obviously many more outliers during the whole experiment, especially in the upward direction. This means that additional INS can remarkably improve the reliability of the solution in severe environments.
On the contrary, the results from the SPP method showed lots of outliers larger than 20 m, which means a poor accuracy was achieved. The results from the KFSPP method were slightly better, but there were also a few outliers at about 10 min before the end of the field test. These could be easily distinguished from the referenced results by the TC RTK/INS method.
Furthermore, the position variations obtained from the SPP, KFSPP, TDCP, and TDCP/INS methods were differenced with the reference values, and the biases of their position variations during the field test were acquired and are shown in
Figure 9. The bias results in the figure directly reflected the accuracies of the position variations obtained from the SPP, KFSPP, TDCP, and TDCP/INS methods.
Generally speaking, the SPP method achieved the worst accuracy, and the KFSPP method had a moderate improvement. Comparatively, the TDCP and TDCP/INS methods achieved more accurate results. As expected, almost 100 times better results were obtained from the TDCP/INS method than those from the SPP method. Due to the lack of INS, there were many more outliers in the TDCP results than the TDCP/INS results, and its robustness decreased. The biases from the TDCP/INS method were almost limited within only 0.05 m, and this performance was comparable to that in the popular RTK positioning.
For a quantitative analysis, the standard deviation (STD) and root mean square (RMS) values of the position variations obtained from the four methods were calculated and are listed in
Table 4. From the table, we can find that, compared with the results from the TC RTK/INS method, the biases of the position variations in both the SPP and the KFSPP were at the meter level. Surprisingly, the biases of the position variations in the TDCP method were still at the meter level, with only a slight decrease compared to the KFSPP method. This should be attributed to the poor quality of the GNSS signals in severe environments, which lead to a lot of outliers. With the help of the INS data, the biases in the TDCP/INS method were dramatically reduced to the centimeter level. The STD values and the RMS values were much similar for all three methods, which implies the results from different methods are consistent regardless of their noise differences.
Specifically, compared with the worst results from the SPP method, the accuracies of the position variations obtained from the KFSPP and TDCP methods had considerable improvements of about 32.1%, 6.9%, and 29.8% and 39.5%, 26.2%, and 43.7% in the east/north/upward directions, respectively, and that, from the TDCP/INS method, improved by nearly 100 times. The STD and RMS values of the TDCP/INS method reached the centimeter level, and the accuracies were 1.9 cm, 2.9 cm, and 3.1 cm in the east/north/upward directions. As a result, the TDCP/INS method obtained the most accurate position variations and, thus, the most stable moving reference positions.
3.4. Analysis of the Trajectory Consistency
After the quantitative analysis of the accuracy of the position variations from the SPP, KFSPP, and TDCP/INS methods, we further made a qualitative analysis of the trajectory consistencies from these methods. The ground trajectory of the vehicle obtained from the TC RTK/INS method was used as the benchmark trajectory, and the position results obtained from the SPP, KFSPP, TDCP, and TDCP/INS methods were compared with the benchmarks to analyze their consistencies.
For the above purpose, three typical urban driving scenarios, which included going along a large, curved road, turning at a crossroad, and passing under an urban viaduct, were selected for the discussion. The comparative vehicle trajectories obtained from the different methods with respect to the benchmarks from the TC RTK/INS method in the above three scenarios are demonstrated from
Figure 10,
Figure 11 and
Figure 12, where the red pins, the blue pins, the green pins, the orange pins, and the yellow pins represent the position results from the TC RTK/INS method, the SPP method, the KFSPP method, the TDCP method, and the TDCP/INS method, respectively.
From
Figure 10 and
Figure 11, we can see that the ground trajectory with the red pins obtained from the TC RTK/INS method was very smooth when the vehicle went along the curved road and turned at the crossroad. Additionally, the density of the red pins also reflected the driving velocity of the vehicle, since the pins in the figures are marked in equal time intervals of 1 s. Therefore, it clearly reflects in
Figure 10 that the vehicle went along the curved road at an almost-constant speed, while we found obvious deceleration and acceleration of the vehicle before and after it turned at the cross in
Figure 11. Note that traffic keeps on the right in China. From this perspective, the position results from the TC RTK/INS method can be reliably regarded as the “true” trajectory.
We first investigated the blue pins in
Figure 10a, which denoted the results from the SPP method, and found obvious jumps along the true trajectory when the vehicle drove from the east–south direction into the curved section. In
Figure 11a, the blue pins changed from the left side of the true trajectory to its right side when the vehicle turned at the cross. All these results imply that inaccurate results were obtained from the SPP method.
The green pins, which denoted the results from the KFSPP method, had a moderate improvement compared with the blue pins. Few jumps could be found for the green pins in
Figure 10b, but the changes from the left side of the true trajectory to its right side were still observed when the vehicle turned at the cross in
Figure 11b, just like what the blue pins showed. This indicated that the accuracy of the KFSPP method was still unsatisfactory.
In contrast, the orange and yellow pins, which denoted the results from the TDCP and TDCP/INS methods, were very consistent with the red pins, despite a fixed offset between them in
Figure 10c,d and
Figure 11c,d. We also understand that the offset was caused by the SPP result in the very beginning epoch and had no effect on the trajectory consistency. Additionally, the orange and yellow pins had almost the same densities as that for the red pins, which demonstrated that the TDCP and TDCP/INS methods could achieve almost the same levels of performance as the benchmarks from the TC RTK/INS method in open environments.
The position results when the vehicle passed under a viaduct demonstrated the performances of different methods when the signals were seriously interrupted, as can be seen in
Figure 12. Poor positioning results were still witnessed for the SPP method and the KFSPP method, since the trajectories formed by the blue pins in
Figure 12a and the green pins in
Figure 12b swayed left and right along the true trajectory formed by the red pins. In contrast to the results from the former two scenarios, in the trajectory formed by the orange pins in
Figure 12c, there were also jumps around the true trajectory, just like the SPP and KFSPP methods. Additionally, there were result interruptions when the vehicle just passed under the viaduct for all four methods. This was caused by the loss of satellite signals and could be confirmed by the number of visible satellites in
Figure 7. The limited satellite number made positioning impossible in that epoch. Fortunately, the position results of the TDCP/INS method were right immediately after the vehicle passed the viaduct, and its trajectory was almost in parallel with the true trajectory.
Since the position results from the TDCP/INS method in a single moving device were consistent with the ground true trajectory, the kinematic positioning results of the entire moving platform system were also consistent with the ground true trajectory. From this perspective, the TDCP/INS method realized the stable position datum and obtained a consistent moving trajectory with the true one without dependency on additional static reference stations.
4. Discussion
In the “static and dynamic” mode for a moving platform system, the coordinate of the reference station is fixed as the position datum of the entire platform system. Therefore, no matter how the platform system moves, once the relative positions are determined by the kinematic relative positioning, the position of every device in the platform system can be referred to the static position datum, and the positioning results of all devices must be consistent with their true motion trajectories.
However, in the “kinematic-to-kinematic” mode, the inaccurate coordinates of the moving reference device obtained through the noisy SPP or KFSPP methods will cause a stochastic swaying of the entire platform system, along with the time, despite that the relative positions of each device in the system can be accurately described. The key issue lies in the inconsistency between the selected position datum of the moving platform system and its true position. Since accurate relative position can be easily obtained, this inconsistency will eventually lead to inconsistency between the positioning results of the entire moving platform system and the actual motion trajectory. As a consequence, the safe distance and coordinated operation between each device in the moving platform system cannot be reliably guaranteed.
The alternative solution proposed in this paper chose to determine the position variations between adjacent epochs rather than obtain the direct position information. The TDCP observations were applied to obtain accurate position variations, and then, the position of the moving device could be precisely updated. The INS data was further integrated with the GNSS-derived position to enhance the reliability of the updated position. This indirect method demonstrated a good performance in the field test.
Once the position datum of the moving platform system became stable, without dependency on additional static reference stations, a novel operational pattern for a platform system with multiple moving devices was expected. In this pattern, any moving devices in the platform system could be regarded as the position datum for itself and acquire the relative positions of the nearby devices. It was clearly shown that this pattern was decentralized, and the breakdown of any devices in the system had no effect on the normal operation of other devices and the entire platform system. This is very different from the current pattern, in which the system operation relies a lot on the static reference station near the platform. From this perspective, the expected pattern was more robust in the system level.
5. Conclusions
In the conventional relative positioning in a platform system with multiple moving devices, the positioning results obtained by the SPP or KFSPP methods were inconsistent with the actual motion trajectory due to the instability of the position datum. To solve this problem, we proposed the TDCP/INS method to transfer and update the position datum of the moving platform via the position variations from the TDCP observations and the INS data. The positioning results of the moving platform were consistent with the actual motion trajectory.
A field test using a driving vehicle was carried out, and the performance of the proposed TDCP/INS method was analyzed and compared with the SPP method, KFSPP method, TDCP method, and the referenced TC RTK/INS method. The experiment results showed that there were large biases of the position variations for the SPP and KFSPP methods. The accuracies of the calculated position variations for the SPP method were 1.148 m, 1.635 m, and 2.499 m in the east/north/upward directions, which, for the KFSPP and TDCP methods, had improvements of 32.1%, 6.9%, and 29.8% and 39.5%, 26.2%, and 43.7%, and which, for the TDCP/INS method, improved by nearly 100 times and reached the centimeter level.
A further qualitative analysis demonstrated the perfect consistency between the position results from the TDCP/INS method and the reference positions in three typical urban driving scenarios, including going along a large, curved road, turning at a crossroad, and passing under an urban viaduct. The motion states of the vehicle, such as the acceleration and deceleration, were also reflected with the TDCP/INS method. On the contrary, there were frequent jumps and sways along the driving direction, and obvious outliers of the positioning results were detected with the SPP, KFSPP, and TDCP methods. The comparative analysis indicated an outstanding positioning performance from the proposed TDCP/INS method, which could achieve a stable position datum in the platform system with multiple moving devices.
In a future work, a platform system with UAVs will be established, and a corresponding field test will be carried out. Additionally, quality control algorithms for the INS data, especially the blunder detection and noise elimination, are planned to be employed to further improve the stability and reliability of the positioning results.