Next Article in Journal
LiReD: A Light-Weight Real-Time Fault Detection System for Edge Computing Using LSTM Recurrent Neural Networks
Previous Article in Journal
A Novel Eddy Current Testing Error Compensation Technique Based on Mamdani-Type Fuzzy Coupled Differential and Absolute Probes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real Time Precise Relative Positioning with Moving Multiple Reference Receivers

College of Automation, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(7), 2109; https://doi.org/10.3390/s18072109
Submission received: 11 May 2018 / Revised: 23 June 2018 / Accepted: 27 June 2018 / Published: 30 June 2018
(This article belongs to the Section Remote Sensors)

Abstract

:
The stationary reference receiver with precisely known coordinates is difficult to establish in some special real-time relative positioning applications. To improve the relative position estimation accuracy and the reliability simultaneously for the RTK without a precisely known reference receiver, multiple receivers mounted on a moving platform are used as the base station. A code and phase measurement fusion model is proposed to reduce the communication burden and generate measurements at any virtual position where it is inconvenient to install the GPS receiver. To keep the integer property of the ambiguity of fused phase measurements, the RTK method with the moving reference receivers is proposed by implementing the integer ambiguity transformation and error absorption strategy based on the known geometry of multiple receivers. Static and kinematic experiments were carried out to evaluate the performance of the proposed relative positioning method. When compared with the single-receiver solution, static results have shown that the proposed method can improve position accuracy by 15.9% and 15.7% for the horizontal and the vertical component, respectively. The kinematic results have shown that the proposed method can achieve position accuracy enhancement by 26.9% for the vertical component.

1. Introduction

Real Time Kinematic (RTK) has been proven to be an efficient and reliable technology for high-precision positioning over the past decades [1]. By resolving integer ambiguity, centimeter-level position estimation accuracy can be achieved in both real-time and post-processing modes [2,3]. Since the RTK outputs the relative position between the base station and the rover, the accurate absolute position of the rover needs to use the coordinates of a stationary reference station. Many applications require the relative position, such as the positioning of cranes on ships, aerial docking of spacecraft, formation flying, relative positioning of vehicles, fleet management, and deformation monitoring of large moving vehicles [4,5]. In such applications, it is difficult to establish a stationary reference station with accurately known coordinates. The moving reference receiver cannot give its precise position as a priori, and furthermore, as it is in motion, the measurements can be easily corrupted by, for example, multipaths [6,7]. The traditional RTK technique cannot be applied directly to moving reference receiver without precise position [4,8].
The research on moving reference receivers was started in the 1990s [9]. Binning studied the single difference positioning method for satellite to satellite navigation using simulated data [10]. Kawano used dynamic to dynamic technology based on pseudo difference for on-orbit satellite docking [11]. The idea of using multiple receivers to improve position estimation accuracy or to determine attitude was explored by several research groups [8,12,13,14,15]. Since carrier phase integer ambiguity resolution plays an important role in the precise positioning and attitude determination, there are many documented studies focusing on ambiguity resolution in the ambiguity domain [16,17,18,19]. These studies incorporate known baseline constraints in the ambiguity search algorithms to get efficient integer search strategies. Another class of studies mainly enhances the performance of float solutions of ambiguities in the positioning domain. Since multiple receivers mounted on the platform in an array of known geometry increases redundancy [20], these studies use multiple receiver configurations to strengthen the observational model [12,21,22,23]. In these studies, multiple receivers are utilized without considering the communication burden, and the baseline constraints are used to strengthen the position estimation model to improve the relative positioning performance using a stationary base station. Since multiple receiver configuration improves positioning performance efficiently, it can be also used in the base station. For fixed base stations, the technique using multiple reference receivers is known as the Network RTK (NRTK). In the NRTK technique, multiple reference receivers are stationary with precisely known coordinates, and the inter-station baseline distances in the current NRTK implementations are typically restricted to around 50 km or less [24,25]. Therefore, the high-precision network-derived atmospheric corrections in both spatial and temporal domain can be computed to correct the rover measurement error when ambiguity between reference stations are fixed [26]. However, these technologies cannot be applied in the RTK with multiple moving reference receivers. Most of the applications mentioned above are practically implemented in air or sea. They encounter three problems. First, land-based GPS stationary stations cannot be used directly, because the users are far away from land. Second, the reference receivers used are usually in motion. Third, the accurately absolute coordinates of reference receivers are difficult to determine. Therefore, modifications are required for moving RTK with multiple moving reference receivers. Luo proposed the ‘Multi-Kin’ method to process GPS observations from multiple (more than three) moving platforms simultaneously, making full use of constraints through the multiplicity of platforms to improve the ambiguity resolution [5]. The ‘Multi-Kin’ method is only suitable for applications that have multiple moving platforms, each of which is equipped with one receiver. Schrader used multiple GPS receivers and supporting electronic components to improve the accuracy of GPS data [27]. Trinklein employed truncated mean and moving average filters of the latitudes and longitudes to post process multiple moving inexpensive GPS receiver data, and yielded approximately one meter accuracy for the relative position vector [14].
In some critical real-time systems, both measurements and positioning performance can be deteriorated severely under challenging conditions, such as electromagnetic interference. In this case, multiple moving reference receiver configuration as the base station is needed to improve the situation. As we use all the measurements from multiple reference receivers, the communication burden between the reference receivers and the user has to be considered. To improve relative positioning performance without affecting the real time performance of positioning, a real time code and phase measurement fusion model is proposed. The model makes full use of known antenna geometry, and similar ionospheric and tropospheric delays to fuse the original measurements from different reference receivers on a moving platform. Moreover, to avoid destroying the integer property of ambiguity of fused phase measurements, both the integer ambiguity transformation and error absorption are applied. Errors that may be introduced in the process of generating fused measurements are evaluated. The static experiment and kinematic experiment were formed to assess the performance of proposed methodologies.
The contribution is organized as follows. In Section 2, the details of the model are outlined, including the original un-differenced phase and code observation models, measurements fusion model, integer ambiguity re-parameterization, and stochastic error model. In Section 3, experiments and analyses of the methodology are given. In Section 4, the conclusion remarks are presented.

2. Methodology

2.1. Functional Model

Considering the satellite and receiver clock biases and various propagation delays, the original un-differenced phase, and code observation models on frequency j can be expressed as [28,29,30]:
P r , j = R r + μ j I r + T r + d t r , j d t , j + ε P r , j ,
λ j ϕ r , j = R r μ j I r + T r + λ j N r , j + δ t r , j δ t , j + ε ϕ r , j ,
where the subscripts r and j denotes the receiver and frequency fj, which are used to emphasize the receiver-specific and frequency-specific term, respectively. Pr,j = [Pr,j1, …, Pr,js]T and φr,j = [ϕr,j1, …, ϕr,js]T are the un-differenced code measurement in meters and phase measurement in unit of cycle, and the superscripts s denote the satellite. Rr = [Rr1, …, Rrs]T is geometric distance; Ir = [Ir1, …, Irs]T is the ionosphere delay on frequency f1 with uj = f12/fj2; Tr = [Tr1, …, Trs]T is tropospheric delay; dtr,j and dt,j = [dt,1, …, dt,s]T are the receiver and satellite clock offset for code respectively, which includes clock error and frequency-dependent code biases; δtr,j and δt,j = [δt,1, …, δt,s]T are the receiver and satellite clock offset for phase respectively, which includes clock error and frequency-dependent phase biases; λj is the phase wavelength and Nr,j = [Nr,j1, …, Nr,js]T is un-differenced integer ambiguity; εPr,j is the receiver code noise in meters, including multipath error in codes; εϕr,j is the receiver phase noise in meters, including phase multipath errors.

2.2. Measurement Fusion Model

The measurement fusion model can combine the code and phase measurements from different reference receivers to generate fused measurements at any virtual locations within the network of reference receivers. Multiple reference receivers can be deployed in different locations of the common moving platform, where there are minimal electromagnetic interferences and multipath effects. Constructing fused measurements in real time requires an accurate knowledge of the baseline lengths between the reference receivers and the position of the reference point. Figure 1 schematically depicts the measurements transferred from multiple receivers, taking four receivers as an example. We choose the point R as the location of a virtual receiver and transfer the code and phase measurements from the reference receivers to this reference point. Point A, B, C, D are reference receivers, and yr = [Pr, φr]T denotes the code and phase measurements from the reference receivers.
To generate fused measurements at the virtual point R, we need to convert the measurements from different receivers, as shown by the red dotted lines in Figure 1. As various reference receivers are installed on the same moving platform, the distances between each receiver and the point R are relatively short. Therefore, their satellite line-of-sight vectors are assumed to be parallel with each other. The projection of the baseline vector on the satellite observation vector is equal to the change of geometric distance between satellite and receiver. Taking advantage of the close distance between reference receivers, we can assume that the atmospheric delays of the satellite signals at all receivers (A, B, C, and D) are similar. Therefore, the main difference of measurement between the reference point R and each receiver is the difference of their geometric distance between satellite and receiver caused by the baseline vector. For receiver r, the transferred measurements at the point R can be given as:
y R , r = y r + Δ R R , r + ε y R , r ,
Δ R R , r = e r , e b R r , e ,
where yR,r = [PR,r, λϕR,r]T, PR,r = [PR,r,1, …, PR,r,n]T and ϕR,r = [ϕR,r,1, …, ϕR,r,n]T are the transferred measurements from receiver r. λ = [λ1, …, λn]T is the phase wavelength of different frequency. εPR,r and εyR,r = [εPR,r, εϕR,r]T is the transferred measurements noise in meters. ΔRR,r is the changed geometric distance between the reference point R and receiver r. er,e is the unit observation vectors of receiver r; the subscripts e denote the WGS-84 ECEF coordinate frame; bRr,e is baseline between the receiver r and the point R.
Equations (1)–(4) can be used to generate the transferred measurements from receiver r. The fused code and phase measurements on point R can be expressed as:
y R = r = 1 m α r y R , r ,
where yR = [PR, ϕR]T, αr is weighting coefficients, and m denotes the number of the reference receivers.

2.3. Integer Ambiguity Re-Parameterization

For the code measurement, Equation (5) provides the final fused code measurements of point R. However, for the fused phase measurement, the introduction of weighting coefficients αr will destroy the ambiguity’s integer property. From Equations (2), (3) and (5) the final fused phase measurements of satellite s at the point R can be generated as:
λ j ϕ R , j s = ρ R , j s + λ j N R , j s + r = 1 m α r δ t r , j δ t , j s + ε ϕ R , j ,
N R , j s = r = 1 m α r N r , j s ,
ρ R , j s = r = 1 m α r ( R r s + Δ R R , r s μ j I r s + T r s ) ,
ε ϕ R , j = r = 1 m ( α r ε ϕ r , j ) 2 .
In Equation (6), N R , j s is not an integer. Therefore, we need to re-parameterize the combined integer ambiguity. The multiple reference receiver configuration makes it possible to form ambiguity transfer. In this approach, the ambiguities between moving reference receivers are used. A main reference receiver needs to be chosen from all the reference receivers. Without the accurately known coordinates of the reference receivers, we use the baseline length constraints to strengthen the double-differenced (DD) phase and code measurement equations. The DD integer ambiguity among the reference receivers can be obtained reliably. Using the DD integer ambiguity, N R , j s in Equation (6) can be expressed as:
N R , j s = r = 1 m α r N 1 r , j 1 s + N 1 , j s + r = 1 m α r ( N r , j 1 N 1 , j 1 )
where the subscript 1 denotes the main reference receiver and the superscripts 1 denotes the main reference satellite. The term α r ( N r , j 1 N 1 , j 1 ) is only affected by the receiver for a fixed frequency. It has the same properties as the receiver clock offset, therefore it can be absorbed by the receiver clock offset as:
δ t R , j = r = 1 m α r δ t r , j + λ j r = 1 m α r ( N r , j 1 N 1 , j 1 )
The DD integer ambiguity between reference receivers which is N 1 r , j 1 s in Equation (10) can be calculated based on the known body-frame geometry of multiple reference receivers [22]. Therefore, the final fused phase measurements to point R can be given as:
λ j ϕ ˜ R , j s = ρ R , j s + λ j N 1 , j s + δ t R , j δ t , j s + ε ϕ R , j ,
where N 1 r , j 1 s represents the ambiguity in the fused phase measurements. It can be found that the combined integer ambiguity is not only immune to the coefficient αr as 1/m, but also preserves the integer property of the integer ambiguity.

2.4. Stochastic Error Model

Since the stochastic model plays an important role in resolving the integer ambiguity [31,32], we need to characterize the stochastic error model of the fused measurements. Based on Equations (1) and (2), we define the observation vector from each receiver yr = [Pr,1, …, Pr,n, ϕr,1, …, ϕr,n]T, where the subscripts n represents the number of frequency. Qr, the stochastic model of yr, can be specified by:
Q r = Q m Q w ,
where Qm = blkdiag(Qp, Qϕ) stands for the precision of phase and code measurement. Qp = diag(σp12, …, σpn2) and Qϕ = diag(σϕ12, …, σϕn2), where σpj2 and σϕj2 are the variance of the un-differenced phase and code measurements. It is assumed that the precision of phase and code measurements are unique for both frequencies. Qw is the weight matrix of the un-differenced measurements. During the construction of the fused GNSS measurements, accurate knowledge of the baseline length relative to WGS-84 ECEF coordinate frame between each receiver and the virtual point R is required. If the baseline length in ECEF frame is unknown, it needs to be calculated by baseline length in the body frame and attitude in the navigation frame. Therefore, along with the GPS measurement errors of each receiver, the baseline length errors between each receiver and point R are of great importance in the generation of fused measurements. The baseline lengths in the body coordinate frame can be determined from a platform survey or other methods which are transferred into the WGS-84 ECEF coordinate frame with the knowledge of the attitude and position of the platform. Therefore, the calculation of baseline length in WGS-84 ECEF coordinate frame is coupled with uncertainties due to baseline length errors in body frame and attitude errors of the platform. An expression for baseline length b ^ R r , e can be simply written as
b ^ R r , e = b R r , e + δ b R r , e ,
in which b ^ R r , e is the estimated value of baseline length used in Equation (4), which is corrupted by estimation error δbRr,e. With the baseline length in body coordinate frame, b ^ R r , e can be obtained as
b ^ R r , e = T ^ n e T ^ b n b ^ R r , b ,
where the subscripts n and b denotes the navigation coordinate frame and the body coordinate frame respectively. The superscript and subscripts e denote the WGS-84 ECEF coordinate frame. T ^ n e is the transformation matrix from navigation frame to WGS-84 ECEF frame and T ^ b n is the transformation matrix from body frame to navigation frame. The hat notation indicates that the terms have nominal errors.
Due to the platform attitude errors, T ^ b n can be expressed as
T ^ b n = T b n + C T b n ,
where C is the skew symmetric matrix caused by attitude errors and determined by the rotation vector of attitude errors in the navigation frame. Compared with attitude errors of the platform, the errors of the transformation matrix caused by platform position estimation error are negligible, i.e., T ^ n e = T n e .
Therefore, the baseline vector errors δbRr,e can be expressed as follows:
δ b R r , e = T n e T b n δ b R r , b + T n e C T b n b R r , b + T n e C T b n δ b R r , b ,
in which δbRr,e contains three parts. The first part TneTbnδbRr,b in the right hand represents baseline length survey errors, which are mainly affected by the survey methods. The second part TneCTbnbRr,b is determined by the attitude errors of the platform The third part couples the survey errors with attitude errors. In this paper, we assume that the attitude and baseline length survey errors are small enough, so that the third term can be neglected.
The covariance for baseline vector errors δbRr,e can be equated as
Q b = E ( δ b R r , e δ b R r , e T ) = R D ( δ b R r , b ) R T + M D ( δ e ) M T ,
M = T n e [ 0 b R r , n ( 3 ) b R r , n ( 2 ) b R r , n ( 3 ) 0 b R r , n ( 1 ) b R r , n ( 2 ) b R r , n ( 1 ) 0 ] ,
b R r , n = T b n b R r , b ,
where E and D are the expectation and covariance operators, respectively, and R = T n e T b n . The first term in the right hand is caused by baseline length errors, while the second term is caused by attitude errors. Equation (18) can be used to assess the stochastic model of the baseline errors. A characterization of baseline errors will benefit ambiguity resolution using the fused measurements. From Equations (14)–(20), it is clear that the term δbRr,e is caused by baseline errors in body frame and the rotation vector of attitude errors in the navigation frame. Therefore, the transferred measurement yR,r in Equation (3) includes the translating error caused by δbRr,e in translating process. The weighting coefficients αr in Equation (5) is set to 1/m, which means the fused measurements are generated by calculating the mean of the transferred measurements. It is noted that when making an average operation, errors and variances are pulled down at the same time. Therefore, the measurement fusion model is helpful in diminishing the influence of the individual measurement error.

3. Experiment and Analysis

We have verified the proposed method by conducting the static and kinematic experiments. We extensively compared the ambiguity resolution and the positioning accuracy performance of the single receiver-antenna (SA) solution and multiple moving receiver-antenna (MA) solution. The SA solution means using a single antenna and receiver as the base station. Both experiments use dual-frequency code and phase GPS measurements. The standard deviation of un-differenced code and phase measurement is 30 cm and 3 mm respectively, by keeping measurement type weighting. It has been widely proven that the LAMBDA algorithm has the best success rate among the different integer aperture estimation methods [33]. We therefore used the LAMBDA for ambiguity resolution. Since the proposed method can mitigate the effect of measurement noise from the combination of multiple reference receivers, an optimistic threshold is preferable to prove the edge of our proposed method. The threshold for ratio test is set up to 1/2 [34,35]. The ambiguity success rates are computed by comparing the single-epoch estimated ambiguities to the reference ambiguities obtained from the whole span of data. To evaluate position estimation accuracy, the horizontal and vertical position error (HPE/VPE) are:
H P E = Δ x E 2 + Δ x N 2 ,
V P E = | Δ x U | ,
where ΔxE, ΔxN, ΔxU is the relative position estimation error of east, north and up component respectively. To quantify relative position estimation precision, assuming Rayleigh and Gaussian distribution of horizontal and vertical error, Equations (23) and (24) are used to compute the 95% horizontal and vertical error bounds, where operator std() denotes standard deviation operator [36].
σ H , 95 % = 2 × s t d ( Δ x E ) 2 + s t d ( Δ x N ) 2
σ V , 95 % = 1.96 × s t d ( Δ x U )

3.1. Static Experiment

In the experiment, the performance of real-time relative positioning using MA solution are compared with commonly used SA solutions. The data from Curtin GNSS Research Centre was used as multiple reference receiver data in the static experiment. There are four stations marked CUTB0, CUTC0, CUT00, CUTA0 on the roof of the building of Curtin University in Perth, Australia, whose locations are precisely known. Each station is equipped with a TRM59800.00-SCIS antenna (Trimble Inc., Sunnyvale, CA, USA) which is connected with TRIMBLE NETR9 receiver (Trimble Inc., Sunnyvale, CA, USA). We assumed that entire roof would act as a base station; the receivers of CUTB0, CUTC0, CUT00 are used as reference receivers. We selected the location of station CUTA0 as the virtual point. Given the coordinates of each antenna in ECEF frame, we calculated the baseline length between the multiple reference receiver antennas and the virtual point. The distribution of CUTB0, CUTC0, CUT00, CUTA0 and the calculated baseline lengths are shown in Figure 2. The details of the data used including baseline length between each reference receiver and the virtual point, the time of data collection, and the sampling rate are given in Table 1.
To simulate the moving base station scenario, we assume that the absolute position of CUT00, CUTB0, CUTC0, CUTA0 are unknown, given the baseline lengths. We set CUT00 as the rover. Figure 3 shows the number of the visible GPS satellites and the positional dilution of precision (PDOP) value for baseline between rover and the base station during the experiment, which helps to evaluate the observing conditions. The cutoff elevation angle is 5 degree. The number of visible GPS satellites varies from 7 to 13, while the PDOP values were smaller than 2 at most times.
With the baseline lengths shown in Table 1, the fused code and phase measurements at the virtual point can be calculated as the measurements of base station. We can get the GPS double difference (DD) relative position estimation between the virtual point and the rover. The estimation accuracy and precision are shown in Figure 4 and Table 2.
The experiment initially demonstrates the feasibility of the multiple moving receiver data fusion method. Considering all receivers of the base station and rover used in this experiment are installed at the top of the building in an open field of vision, the number of GPS satellites received and satellite geometry is appropriate. In addition, the baseline between the reference point and rover can be considered a short baseline. Therefore, in Figure 4, the difference in position estimation accuracy obtained by using MA and the SA is not obvious. In Table 2, the position estimation precision using MA and the SA is at the same level, and the ambiguity resolution success rates are all 100%.
To test the relative positioning performance with a longer baseline distance, we use PERT as the rover instead of CUT00. PERT is a permanent GNSS reference station in Australia in the IGS network, which is equipped with a TRM59800.00 antenna and TRIMBLE NETR9 receiver. The distance between PERT and the virtual point is about 22.4 km. To mitigate the effect of DD measurement residual errors, we increased the elevation angle to 10 degrees for the long baseline test. The visible GPS satellites and the PDOP value for baseline between rover and the base station are shown in Figure 5, with a cutoff elevation angle of 10 degrees. During the experiment, 5 to 13 satellites were observed at the rover station. The PDOP values were oscillating near the value of 2 most of the time.
With the same fused code and phase measurements at the virtual point, the relative position estimation accuracy and precision are shown in Figure 6 and Table 3.
The relative position estimation errors in Figure 6 are increased compared with Figure 4 for both horizontal and vertical components. The position accuracy of MA solution is comparable with the SA solution. When compared with the SA, it can found that the ambiguity resolution success rate of MA is increased by 3.06%. Furthermore, the relative position estimation precision can be reduced by 15.9% for the horizontal component and 15.7% for the vertical component, compared with SA solution shown in Table 3. The relative position estimation accuracy decreases mainly due to the de-correlation of the ionospheric and tropospheric delays with growing rover to reference point distance. MA solution can obtain better quality measurements and then have better ambiguity resolution success rate and position estimation precision compared with the SA solution. This is mainly because the MA solution pulls down the errors and generates higher precision fused measurements than the individual measurement used in SA solution.
The fused measurements obtained using baseline length in ECEF frame directly are not affected by the platform baseline errors and the attitude errors. However, if the baseline length in ECEF frame is unknown, it needs to be calculated by baseline length in body frame and attitude in navigation frame. In this case, the uncertainty of fused measurements is influenced by platform baseline errors and attitude errors, as shown in Equation (15). We use the static experiment data to evaluate the effect of attitude errors and baseline errors on relative position estimation results. As all the reference receivers are static to simplify simulation, the body frame and the navigation frame are assumed to be coincident. Different error models including baseline errors and attitude errors are used in the simulation.
b ^ R r , b is estimated as the baseline vector in the body frame, which can be obtained by platform survey or estimation in real time. Any estimation process for baseline vector cannot have a perfect solution. Therefore, there will be a residual uncertainty of b ^ R r , b . We assume the uncertainty of baseline has a standard deviation of 1 cm [37] in the simulation. Figure 7 and Table 4 shows the baseline-induced extra position estimation errors with fused measurements. In Figure 7, the cloud of points represents the position estimation errors caused by baseline errors. The red ellipses represent 95% of position estimation errors. From Figure 7 and Table 4, it is seen that baseline uncertainty with a standard deviation of 1 cm causes centimeter-level accuracy changes in position estimation accuracy.
For attitude errors, it is also difficult to specify values as the expected errors in the simulation. This is because the attitude errors mainly depend on grade of attitude sensors. To show the effect of attitude errors on the final position estimation result, we assume that the standard deviation of attitude errors is approximately 0.05° to 0.3° in the simulation [38,39,40,41].
Figure 8 and Table 5 show the relative position estimation errors in east, north, up component, and position estimation precision caused by attitude errors. It is clear that if the attitude standard deviation is small enough, for example smaller than 0.05°, the attitude-induced relative position estimation errors will change in millimeter-level compared with the results without attitude error. Figure 9 and Table 5 show the position estimation error due to the combined effect of baseline errors and attitude errors. The baseline standard deviation is set to 1 cm, while the attitude standard deviation ranges from 0.05° to 0.3°. If the baseline standard deviation is restricted to 1 cm, the combined baseline and attitude-induced position estimation errors have far less difference compared with attitude-induced position estimation errors. This means that the baseline length error with the standard deviation of 1 cm has little effect on the position estimation errors compared with the attitude errors. These figures and tables show the influence of the attitude errors and baseline errors on the position estimation results, thereby providing a reference for the accuracy of baseline and attitude.

3.2. Kinematic Experiment

The kinematic experiment was conducted on the Songhua River of Harbin, China. The receiver antennas are mounted on the roof of the ship, as shown by Figure 10. Three reference receivers (BDM683, Unicorecomm, Beijing, China) are placed in the cabin and connect with two CHCNAV A230GRB antennas and a Novatel GPS-703-GGG antenna, marked as A1, A2, A3, respectively. We use the post-processing software to obtain the coordinate of antennas A1~A3 as the benchmark. We set the receivers RCV1, RCV2, RCV3, which are connected with three antennas A1, A2, A3 respectively, as the reference receivers and set the location of antennas A3 as the reference point (R). The details of the data are given in Table 6. The cutoff elevation angle is set up to 5 degrees.
In the kinematic experiment, the ship sailed on the river for about 39 min. The horizontal trajectory is shown in Figure 11, and the velocity of the ship are shown in the right part of Figure 11. The ship moved irregularly during the experiment, and the three-dimensional velocity ranged from 0 m/s to 3 m/s.
We use the fused measurements at virtual R as the measurements of the base station. The fused measurements are obtained by the measurements from three moving receivers RCV1, RCV2, RCV3 which are connected with three antennas A1, A2, A3 respectively. One receiver located at Harbin Engineering University, China, is used as the rover. The distance between the rover and start point in Figure 11 is about 4.1 km. During the kinematic experiment, the number of visible GPS satellites and PDOP values for baseline between rover and the base station are shown as Figure 12. The visible GPS satellites number varied from 6 to 7, while the PDOP values were lager than 2. With the fused code and phase measurements at virtual point, the GPS DD relative position estimation accuracy and precision are shown in Figure 13 and Table 7.
Figure 13 gives the final relative position estimation accuracy of the kinematic experiment, from which we can see that the position estimation accuracy using MA and the SA is at the same level after ambiguity resolution. However, the ambiguity resolution success rate using MA is increased by 1.33%, which is obtained with a cutoff elevation angle of 5; the position error can be reduced by 26.9% for vertical component when compared with SA solution, as shown in Table 7. Compared with SA solution, the MA solution can reduce the errors of the measurement. Multipath is the dominant error in short-baseline differential GPS systems. The observation data at low elevation angle is more severely affected by multipath noise. MA can mitigate multipath noise to some extent, so it got a better ambiguity resolution rate than SA.

4. Conclusions

In this paper, we use moving multiple reference receivers-based RTK for the real time precise relative positioning. To reduce the communication burden and improve the relative positioning accuracy, we have developed a code and phase measurements fusion method to combine GPS measurements from multiple receivers mounted on a moving platform for real time relative positioning. The proposed method takes advantage of known antenna geometry, integer ambiguity transformation, and error absorption, as well as similar ionospheric and tropospheric delays, to generate the fused measurements.
The presented static and kinematic experiment results were based on single epoch double difference positioning with a moving base station. In the static experiment, when the baseline is increased to 22.4 km, the ambiguity resolution success rate is increased by 3.06%, and the relative position estimation precision is reduced by 15.9% for the horizontal component and by 15.7% for the the vertical component compared with the SA solution. In the kinematic experiment, the ambiguity resolution success rate is increased by 1.33%, and the relative position estimation precision is reduced by 26.9% for the vertical component, compared with the SA solution. Along with the relative positioning performance test, we also evaluated the baseline-induced, attitude-induced, combined baseline and attitude-induced relative position estimation errors in the MA solution with different error models in the simulation. The results showed that if the attitude standard deviation is smaller than 0.05° and the baseline standard deviation is restricted to 1 cm, the 95% horizontal and vertical relative position estimation precision bounds are less than 0.0169 m and 0.0137 m respectively.
The proposed approach is suitable for applications which have multipath noise and require high real time performance. Since the measurements fusion model can convert the measurements of receivers distributed in different locations on the platform to measurements at any point on the platform and obtain the relative position between the point and the user, it can be used in applications that need to determine the relative position of a given point where it is inconvenient to install the GPS receiver.

Author Contributions

H.L., S.G. and L.L. conceived the idea and wrote the paper; H.L. designed the experiments with L.L. and L.Z.; C.J. performed the experiments.

Funding

This research is supported by the National Natural Science Funds of China (Nos. 61773132, 61633008), the Fundamental Research Funds for Central Universities (No. HEUCFP201768) and the Post-Doctoral Scientific Research Foundation, Heilongjiang Province (No. LBH-Q15033).

Acknowledgments

The GPS data sets collected for static experiment from Curtin GNSS Research Centre were used in this study, which is acknowledged.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. He, H.; Li, J.; Yang, Y.; Xu, J.; Guo, H.; Wang, A. Performance assessment of single-and dual-frequency BeiDou/GPS single-epoch kinematic positioning. GPS Solut. 2014, 18, 393–403. [Google Scholar] [CrossRef]
  2. He, X.; Zhang, X.; Tang, L.; Liu, W. Instantaneous Real-Time Kinematic Decimeter-Level Positioning with BeiDou Triple-Frequency Signals over Medium Baselines. Sensors 2016, 16, 1. [Google Scholar] [CrossRef] [PubMed]
  3. Nadarajah, N.; Teunissen, P.J.G.; Sleewaegen, J.-M.; Montenbruck, O. The mixed-receiver BeiDou inter-satellite-type bias and its impact on RTK positioning. GPS Solut. 2015, 19, 357–368. [Google Scholar] [CrossRef]
  4. Maorong, G. Dynamic to Dynamic GPS Real Time Differential Positioning. Geotech. Investig. Surv. 1998, 4, 55–57. [Google Scholar]
  5. Luo, N.; Lachapelle, G. Relative positioning of multiple moving platforms using GPS. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 936–948. [Google Scholar]
  6. Lachapelle, G.; Sun, H.; Cannon, M.E.; Lu, G. Precise aircraft-to-aircraft positioning using a multiple receiver configuration. Can. Aeronaut. Space J. 1994, 40, 74–80. [Google Scholar]
  7. Ahn, Y.W.; Kim, D.; Dare, P.; Langley, R.B. Long baseline GPS RTK performance in a marine environment using NWP ray-tracing technique under varying tropospheric conditions. In Proceedings of the ION GNSS, Long Beach, CA, USA, 13–16 September 2005; pp. 2092–2103. [Google Scholar]
  8. Luo, N. Centimetre-level relative positioning of multiple moving platforms using ambiguity constraints. In Proceedings of the ION GPS2000 Conference, Salt Lake City, UT, USA, 19–22 September 2000; pp. 113–122. [Google Scholar]
  9. Hermann, B.R. Kinematic on-the-Fly GPS Positioning Relative to a Moving Reference. Navigation 1995, 42, 487–501. [Google Scholar] [CrossRef]
  10. Binning, P.W. Absolute and Relative Satellite to Satel Lite Navigation Using GPS. Ph.D. Thesis, Department of Aerospace Engineering Sciences, University of Colorado, Boulder, CO, USA, 1997; pp. 55–61. [Google Scholar]
  11. Kawano, I.; Mokuno, M.; Kasai, T.; Suzuki, T. First Autonomous Rendezvous Using Relative GPS Navigation by ETS-VII. Navig. J. Inst. Navig. 2001, 48, 49–56. [Google Scholar] [CrossRef]
  12. Nair, S. A Multiple Antenna Global Positioning System Configuration for Enhanced Performance. Master’s Thesis, Ohio University, Athens, OH, USA, 2004. [Google Scholar]
  13. Renga, A.; Fasano, G.; Simonetti, A.; Accardo, D.; Grassi, M.; Senatore, R. Fusion of Multi-Antenna Carrier Phase Differential GPS and Inertial Measurements for Performance Evaluation of High Accuracy Integrated Aircraft Navigation Systems. In Proceedings of the AIAA Infotech at Aerospace Conference and Exhibit, Garden Grove, CA, USA, 19–21 June 2012; Volume 51, p. 262. [Google Scholar]
  14. Trinklein, E.; Parker, G. Combining multiple GPS receivers to enhance relative distance measurements. In Proceedings of the 2013 IEEE Sensors Applications Symposium, Galveston, TX, USA, 19–21 February 2013. [Google Scholar]
  15. Dorn, M.; Lesjak, R.; Wieser, M. Improvement of the standard GNSS/IMU solution for UAVs using already existing redundant navigation sensors. In Proceedings of the IEEE 2016 European Navigation Conference (ENC), Helsinki, Finland, 30 May–2 June 2016; pp. 1–9. [Google Scholar]
  16. Teunissen, P.J.G. A general multivariate formulation of the multi-antenna GNSS attitude determination problem Artif. Satellites 2007, 42, 97–111. [Google Scholar]
  17. Giorgi, G.; Teunissen, P.; Buist, P. A search and shrink approach for the baseline constrained LAMBDA method: Experimental results. In Proceedings of the International Symposium GPS/GNSS, Curran Associates, Tokyo, Japan, 11–14 November 2008. [Google Scholar]
  18. Teunissen, P.J.G. Integer least-squares theory for the GNSS compass. J. Geodesy 2010, 84, 433–447. [Google Scholar] [CrossRef] [Green Version]
  19. Giorgi, G.; Teunissen, P.J.G.; Verhagen, S.; Buist, P. Instantaneous ambiguity resolution in global-navigation satellite-system-based attitude determination applications: A multivariate constrained approach. J. Guid. Control Dyn. 2012, 35, 51–67. [Google Scholar] [CrossRef]
  20. Ray, J.; Cannon, M.E.; Fenton, P. GPS code and carrier multipath mitigation using a multiantenna system. IEEE Trans. Aerosp. Electron. Syst. 2001, 37, 183–195. [Google Scholar] [CrossRef]
  21. Paziewski, J. Precise GNSS single epoch positioning with multiple receiver configuration for medium-lengthbaselines: Methodology and performance analysis. Meas. Sci. Technol. 2015, 26, 035002. [Google Scholar] [CrossRef]
  22. He, K.; Xu, T.; Förste, C.; Petrovic, S.; Barthelmes, F.; Jiang, N.; Flechtner, F. GNSS precise kinematic positioning for multiple kinematic stations based on a priori distance constraints. Sensors 2016, 16, 470. [Google Scholar] [CrossRef] [PubMed]
  23. He, K.; Xu, G.; Xu, T.; Flechtner, F. GNSS navigation and positioning for the GEOHALO experiment in Italy. GPS Solut. 2016, 20, 215–224. [Google Scholar] [CrossRef]
  24. Wielgosz, P.; Kashani, I.; Grejner-Brzezinska, D. Analysis of long-range network RTK during severe ionospheric storm. J. Geodesy 2005, 79, 524–531. [Google Scholar] [CrossRef]
  25. Li, B.; Shen, Y.; Feng, Y.; Gao, W.; Yang, L. GNSS ambiguity resolution with controllable failure rate for long baseline network RTK. J. Geodesy 2014, 88, 99–112. [Google Scholar] [CrossRef]
  26. Odolinski, R. Temporal correlation for network RTK positioning. GPS Solut. 2011, 16, 147–155. [Google Scholar] [CrossRef]
  27. Schrader, D.K.; Min, B.-C.; Matson, E.T.; Dietz, J.E. Combining multiple, inexpensive GPS receivers to improve accuracy and reliability. In Proceedings of the Sensors Applications Symposium, Brescia, Italy, 7–9 February 2013. [Google Scholar]
  28. Teunissen, P.J.G.; Odijk, D.; Zhang, B. PPP-RTK: Results of CORS Network-Based PPP with Integer Ambiguity Resolution. J. Aeronaut. Astronaut. Aviat. Ser. A 2010, 42, 223–229. [Google Scholar]
  29. Teunissen, P.J.G.; Kleusberg, A. GPS for Geodesy, 2nd ed.; Springer-Verlag: Berlin, Germany, 1998; Volume 60, pp. 271–318. [Google Scholar]
  30. Zhang, B.; Teunissen, P.J.G.; Odijk, D. A Novel Un-differenced PPP-RTK Concept. J. Navig. 2011, 64, S180–S191. [Google Scholar] [CrossRef]
  31. Teunissen, P.J.G. The success rate and precision of GPS ambiguities. J. Geodesy 2000, 74, 321–326. [Google Scholar] [CrossRef]
  32. Amiri-Simkooei, A.R.; Jazaeri, S.; Zangeneh-Nejad, F.; Asgari, J. Role of stochastic model on GPS integer ambiguity resolution success rate. GPS Solut. 2016, 20, 51–61. [Google Scholar] [CrossRef]
  33. Teunissen, P.J.G. The least aquares ambiguity decorrelation adjustment: A method for fast GPS integer estimation. J. Geodesy 1995, 70, 65–82. [Google Scholar] [CrossRef]
  34. Petovello, M.G.; O’Keefe, K.; Lachapelle, G.; Cannon, M.E. Measuring aircraft carrier flexure in support of autonomous aircraft landings. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 523–535. [Google Scholar] [CrossRef]
  35. Wei, M.; Schwarz, K.P. Fast ambiguity resolution using an integer nonlinear programming method. In Proceedings of the ION GPS, Palm Springs, CA, USA, 12–15 September 1995; pp. 1101–1110. [Google Scholar]
  36. Kube, F.; Bischof, C.; Alpers, P.; Wallat, C.; Schön, S. A virtual receiver concept and its application to curved aircraft-landing procedures and advanced LEO positioning. GPS Solut. 2018, 22, 1–12. [Google Scholar] [CrossRef]
  37. Zhang, X.H.; He, X.Y. Performance analysis of triple-frequency ambiguity resolution with BeiDou observations. GPS Solut. 2015, 20, 1–13. [Google Scholar] [CrossRef]
  38. Zhang, X.; Wu, M.; Liu, W. Receiver time misalignment correction for GPS-based attitude determination. J. Navig. 2015, 68, 646–664. [Google Scholar] [CrossRef]
  39. Hu, P.; Gao, Z.; She, Y.; Cai, L.; Han, F. Shipborne heading determination and error compensation based on a dynamic baseline. GPS Solut. 2015, 19, 403–410. [Google Scholar] [CrossRef]
  40. Cohen, C.E. Attitude Determination Using GPS. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 1992. [Google Scholar]
  41. Vander Kuylen, L.; Nemry, P.; Simsky, A.; Lorga, J.F.M. Comparison of attitude performance for multi-antenna receivers. Eur. J. Navig. 2006, 4, 1–9. [Google Scholar]
Figure 1. Multiple reference receiver configuration architecture.
Figure 1. Multiple reference receiver configuration architecture.
Sensors 18 02109 g001
Figure 2. Distribution of multiple receivers. The (Left) panel is the location of receivers. The (Right) panel is the baseline between each receiver and point R.
Figure 2. Distribution of multiple receivers. The (Left) panel is the location of receivers. The (Right) panel is the baseline between each receiver and point R.
Sensors 18 02109 g002
Figure 3. Visible GPS satellites number and PDOP value.
Figure 3. Visible GPS satellites number and PDOP value.
Sensors 18 02109 g003
Figure 4. Relative position estimation errors using MA and SA. The (Left) panel is horizontal position estimation error. The (Right) panel is vertical position estimation error.
Figure 4. Relative position estimation errors using MA and SA. The (Left) panel is horizontal position estimation error. The (Right) panel is vertical position estimation error.
Sensors 18 02109 g004
Figure 5. Visible GPS satellites number and PDOP value.
Figure 5. Visible GPS satellites number and PDOP value.
Sensors 18 02109 g005
Figure 6. Relative position estimation errors using MA and SA. The (Left) panel is horizontal position estimation error. The (Right) panel is vertical position estimation error.
Figure 6. Relative position estimation errors using MA and SA. The (Left) panel is horizontal position estimation error. The (Right) panel is vertical position estimation error.
Sensors 18 02109 g006
Figure 7. Baseline-induced relative position estimation errors. The (Left) panel is the errors in east and north direction. The (Middle) panel is the errors in east and up direction. The (Right) panel is the errors in north and up direction.
Figure 7. Baseline-induced relative position estimation errors. The (Left) panel is the errors in east and north direction. The (Middle) panel is the errors in east and up direction. The (Right) panel is the errors in north and up direction.
Sensors 18 02109 g007
Figure 8. Attitude-induced position estimation errors.
Figure 8. Attitude-induced position estimation errors.
Sensors 18 02109 g008
Figure 9. Combined error-induced position estimation errors.
Figure 9. Combined error-induced position estimation errors.
Sensors 18 02109 g009
Figure 10. Distribution of multiple receivers and antennas. The (Left) panel is mounted receivers and antennas. The (Right) panel is baseline between the antennas and the selected virtual point.
Figure 10. Distribution of multiple receivers and antennas. The (Left) panel is mounted receivers and antennas. The (Right) panel is baseline between the antennas and the selected virtual point.
Sensors 18 02109 g010
Figure 11. Ship’s trajectory and velocity during the kinematic experiment. The (Left) panel is trajectory. The (Right) panel is velocity.
Figure 11. Ship’s trajectory and velocity during the kinematic experiment. The (Left) panel is trajectory. The (Right) panel is velocity.
Sensors 18 02109 g011
Figure 12. Visible GPS satellites number and PDOP value.
Figure 12. Visible GPS satellites number and PDOP value.
Sensors 18 02109 g012
Figure 13. Relative position estimation errors using MA and SA. The (Left) panel is horizontal position estimation error. The (Right) panel is vertical position estimation error.
Figure 13. Relative position estimation errors using MA and SA. The (Left) panel is horizontal position estimation error. The (Right) panel is vertical position estimation error.
Sensors 18 02109 g013
Table 1. Details of data in static experiment.
Table 1. Details of data in static experiment.
StationDateRover/BaseBaseline (to R)IntervalDuration
CUT0019 April 2017Base 18.42 m30 s24 h
CUTB019 April 2017Base 29.47 m30 s24 h
CUTC019 April 2017Base 35.75 m30 s24 h
CUTA019 April 2017the virtual point (R)0 m30 s24 h
Table 2. Relative position estimation precision and ambiguity success rate.
Table 2. Relative position estimation precision and ambiguity success rate.
UserσH,95% (cm)σV,95% (cm)Ambiguity Success Rate
MA0.610.68100%
SA0.750.81100%
Table 3. Relative position estimation precision and ambiguity success rate.
Table 3. Relative position estimation precision and ambiguity success rate.
UserσH,95% (cm)σV,95% (cm)Ambiguity Success Rate
MA2.903.9396.98%
SA3.454.6693.92%
Table 4. Baseline–induced relative position estimation precision and ambiguity success rate.
Table 4. Baseline–induced relative position estimation precision and ambiguity success rate.
Standard DeviationσH,95% (cm)σV,95% (cm)Ambiguity Success Rate
1 cm1.410.97100%
Table 5. Attitude-induced and combined error-induced relative position estimation precision.
Table 5. Attitude-induced and combined error-induced relative position estimation precision.
Standard DeviationAttitude-Induced OnlyCombined Error-Induced
σH,95% (cm)σV,95% (cm)σH,95% (cm)σV,95% (cm)
0.05°0.940.941.691.37
0.1°1.881.842.322.08
0.15°2.792.783.112.95
0.2°3.923.724.123.84
0.25°4.714.724.894.80
0.3°5.805.625.845.56
Table 6. Details of data in kinematic experiment.
Table 6. Details of data in kinematic experiment.
StationDateUser/BaseBaseline (to R)IntervalDuration
RCV113 November 2017Base 13.22 m1 s38.78 min
RCV213 November 2017Base 22.94 m1 s38.78 min
RCV313 November 2017Base 30 m1 s38.78 min
Table 7. Relative position estimation precision and ambiguity success rate.
Table 7. Relative position estimation precision and ambiguity success rate.
UserσH,95% (cm)σV,95% (cm)Ambiguity Success Rate
MA0.451.28100%
SA0.501.7598.67%

Share and Cite

MDPI and ACS Style

Li, H.; Gao, S.; Li, L.; Jia, C.; Zhao, L. Real Time Precise Relative Positioning with Moving Multiple Reference Receivers. Sensors 2018, 18, 2109. https://doi.org/10.3390/s18072109

AMA Style

Li H, Gao S, Li L, Jia C, Zhao L. Real Time Precise Relative Positioning with Moving Multiple Reference Receivers. Sensors. 2018; 18(7):2109. https://doi.org/10.3390/s18072109

Chicago/Turabian Style

Li, Hui, Shuang Gao, Liang Li, Chun Jia, and Lin Zhao. 2018. "Real Time Precise Relative Positioning with Moving Multiple Reference Receivers" Sensors 18, no. 7: 2109. https://doi.org/10.3390/s18072109

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

Article Metrics

Back to TopTop