Next Article in Journal
Sol-Gel Thin Films for Plasmonic Gas Sensors
Previous Article in Journal
Nonintrusive Finger-Vein Recognition System Using NIR Image Sensor and Accuracy Analyses According to Various Factors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Medium to Long Range Kinematic GPS Positioning with Position-Velocity-Acceleration Model Using Multiple Reference Stations

1
Department of Geoinformatics Engineering, Kyungil University, 50 Gamasilgil, Kyeongsan, Gyeongbuk 712-701, Korea
2
Division of IoT and Robotics Convergence Research, Daegu Gyeongbuk Institute of Science & Technology, 333 Techno jungang-daero, Hyeonpung-myeon, Dalseong-gun, Dagegu 711-873, Korea
3
Department of Geoinformatics, University of Seoul, 163 Seoulsiripdaero, Dongdaemun-gu, Seoul 130-743, Korea
*
Author to whom correspondence should be addressed.
Sensors 2015, 15(7), 16895-16909; https://doi.org/10.3390/s150716895
Submission received: 12 May 2015 / Revised: 9 July 2015 / Accepted: 9 July 2015 / Published: 13 July 2015
(This article belongs to the Section Remote Sensors)

Abstract

:
In order to obtain precise kinematic global positioning systems (GPS) in medium to large scale networks, the atmospheric effects from tropospheric and ionospheric delays need to be properly modeled and estimated. It is also preferable to use multiple reference stations to improve the reliability of the solutions. In this study, GPS kinematic positioning algorithms are developed for the medium to large-scale network based on the position-velocity-acceleration model. Hence, the algorithm can perform even in cases where the near-constant velocity assumption does not hold. In addition, the estimated kinematic accelerations can be used for the airborne gravimetry. The proposed algorithms are implemented using Kalman filter and are applied to the in situ airborne GPS data. The performance of the proposed algorithms is validated by analyzing and comparing the results with those from reference values. The results show that reliable and comparable solutions in both position and kinematic acceleration levels can be obtained using the proposed algorithms.

1. Introduction

For the past two decades, global positioning systems (GPS) has been widely used for precise positioning in various engineering and scientific fields. While its application areas in the early stage were limited to the determination of a position at a static environment, more and more applications which require kinematic positioning have been developed. For the kinematic GPS positioning, either single or multi-reference stations can be used for the data processing. The multi-reference station approach generally offers better positioning accuracy and reliability [1,2]. Therefore, the multi-reference station approach is recommended when precise kinematic positioning is required. However, many distance-dependent errors, such as atmospheric and ionospheric delays, are not fully removed through the double-differencing technique when the baseline length increases. Thus, these errors should be properly modeled for the medium to long range kinematic positioning [3,4,5]. GPS kinematic positioning is usually performed by using the Kalman filter, and so called Position-Velocity (PV) model can be adopted when the object is assumed to move with a nearly constant velocity. In addition, Position-Velocity-Acceleration (PVA) model can be introduced for the cases where the near-constant velocity assumption does not hold [6]. This means that the PVA model can provide not only the flexibility in estimation procedure but also the position, velocity, and kinematic acceleration of the moving object simultaneously.
It is notable that one of the popular applications utilizing the GPS derived acceleration is the airborne gravimetry. The airborne gravity survey is performed by measuring the kinematic accelerations of the aircraft using GPS. Then, the gravity information is extracted using the well-known equation as shown in Equation (1). The Equation (1) expressed in inertial frame shows that the kinematic accelerations consist of the gravitational acceleration and the specific force sensed by an accelerometer [7,8];
x ¨ = g + a
where x ¨ is kinematic acceleration; g is gravitational acceleration; a is specific force.
In general, the GPS measurements are processed in relative positioning mode to obtain accurate and reliable positioning results in 3-D space. Then, the kinematic accelerations are computed by taking second-order time-derivatives of positions [9,10]. This means that two steps of data processing are required to get the kinematic accelerations from the positions, which decrease the efficiency in a computational point of view. To overcome the drawback of this method, one can also use the PVA model approach so that the position, velocity, and acceleration of the aircraft can be determined at the same time. An efficient PVA model for Kalman filter is proposed and demonstrated using the measurements from the electronic tacheometer, i.e., geodetic kinematic measurements [11]. However, it should be noted that the experiments are conducted in an indoor environment, and the measurement used in the experiments are not from the GPS.
In this paper, the PVA model to directly determine not only the positions but also the kinematic accelerations is proposed, and the validation of the algorithm is demonstrated with in situ GPS measurements. The proposed algorithm is designed for the medium to long baseline scenario so that it can perform properly even in a survey for a large area under unfavorable surveying condition. The atmospheric effects such as tropospheric and ionospheric delays are modeled in the Kalman filter for the medium to long baseline scenario. In addition, multi-baseline approach is adopted for the improvement of the accuracy and reliability of the solutions. The algorithms are developed and numerical test are performed to validate the proposed algorithms by analyzing the positioning results and comparing the kinematic accelerations with those from the second-order time-derivatives of positions.

2. Methodology

The Kalman filter is an optimal recursive estimator which uses a system’s dynamics model and sequential measurements to estimate the system’s unknown states in a minimum variance sense. The Equation (2) shows the discrete Kalman filter equations categorized into two components, i.e., the prediction of state vector through the state-transition matrix and the update of the state vector using the measurements [12,13].
x k = Φ k 1 x k 1 + w k 1 w k ~ N ( 0 , Q k ) z k = H k x k + v k v k ~ N ( 0 , R k )
where subscript k indicates the epoch; x k is state vector; Φ k is state-transition matrix; w k is process noise vector; Q k is covariance matrix for the processes noise; z k is measurement vector; H k is design matrix; v k is measurement noise vector; R k is covariance matrix for the measurement noise.

2.1. State Vector and System Equation

The objective of the proposed algorithm is to determine the kinematic accelerations of the aircraft in unfavorable conditions. Therefore, the algorithm is designed for the medium to long baseline scenario implementing multi-baseline data processing. This means that the proper modeling for the atmospheric effects such as tropospheric and ionospheric delays is required and, as a consequence, the zenith wet delay (ZWD) residuals and the double-differenced (DD) ionospheric delays are included in the state vector as shown in Equation (3). The system state vector is composed of positions of reference stations ( x i _ p o s ), position-velocity-acceleration of aircraft ( x j _ p v a ), zenith wet delay residual ( Δ Z W D ), ambiguities for L1 and wide-lane ( x a m b ), and DD ionospheric delay ( x i o n ) as follows:
x k = [ x i y i z i x i _ p o s x j x ˙ j x ¨ j y j y ˙ j y ¨ j z j z ˙ j z ¨ j x j _ p v a Δ Z W D i Δ Z W D j x Z W D N i j , 1 k N i j , w k x a m b I i j k / f 1 2 x i o n ] T
where subscript i and j denote reference station and the moving aircraft, respectively; superscript k and are the satellite indices; x , y , z are the position components; x ˙ , y ˙ , z ˙ and x ¨ , y ¨ , z ¨ are velocities and accelerations of the aircraft; Δ Z W D is the total zenith tropospheric delay residual; N i j , 1 k and N i j , w k ( = N i j , 1 k N i j , 2 k ) are the DD ambiguities for L1 and wide-lane, respectively; I i j k / f 1 2 is the DD ionospheric delay.
The corresponding state-transition matrix can be expressed by:
Φ = [ Φ i _ p o s 0 0 0 0 0 Φ j _ p v a 0 0 0 0 0 Φ Δ Z W D 0 0 0 0 0 Φ a m b 0 0 0 0 0 Φ i o n ]
The coordinates of reference stations are included in the state vector and estimated with an assumption of random constant stochastic process, so the transition matrix for those are given as in Equation (5).
Φ i _ p o s = [ 1 0 0 0 1 0 0 0 1 ]
The state-transition matrix for position-velocity-acceleration components of the aircraft can be derived with the assumption that the kinematic acceleration is Gaussian white noise process [6,11,14]:
Φ j _ p v a = [ 1 Δ t ( Δ t ) 2 2 0 0 0 0 0 0 0 1 Δ t 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 Δ t ( Δ t ) 2 2 0 0 0 0 0 0 0 1 Δ t 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 Δ t ( Δ t ) 2 2 0 0 0 0 0 0 0 1 Δ t 0 0 0 0 0 0 0 0 1 ]
For the modeling of tropospheric effect, either random walk or first-order Gauss-Markov process is used because either process is known to be sufficient to describe the stochastic property of tropospheric effect [4]. Here, the zenith wet delay residual, i.e., increment with respect to a priori value, is assumed to be random walk stochastic process and the transition matrix is:
Φ Δ Z W D = [ 1 0 0 1 ]
The remaining atmospheric effect, i.e., DD ionospheric delays, is modeled with first-order Gauss-Markov process. The first-order Gauss-Markov process can describe a large number of physical processes with reasonable accuracy [15]. However, it should be noted that the DD ionospheric delays are correlated with not only time but also baseline length. Hence, the DD ionospheric delays are modeled with the first-order Gauss-Markov process which incorporates the effect of both the time and baseline length changes as shown by [3].
Φ i o n = e ( Δ t / T + | δ | / D )
where T is first-order correlation time for the DD ionosphere; | δ | is baseline length change; D = 1500 km is the first-order correlation distance.
The DD ambiguities for L1 and wide-lane measurements can be modeled as a random constant process as the case of reference stations’ coordinates. It should be noted that the ambiguity resolution is not attempted because it is an issue outside the scope of this study in long range static or kinematic positioning.
For the Kalman filtering, the process noise covariance matrix should be properly modeled as it confines the variability of the state estimates. The discrete form of covariance matrices for the reference stations’ coordinates and the position-velocity-acceleration of aircraft can be written as shown in Equations (9) and (10), respectively [6,11].
Q i _ p o s = [ 0 0 0 0 0 0 0 0 0 ]
Q j _ p v a = q a c c [ ( Δ t ) 5 20 ( Δ t ) 4 8 ( Δ t ) 3 6 0 0 0 0 0 0 ( Δ t ) 4 8 ( Δ t ) 3 3 ( Δ t ) 2 2 0 0 0 0 0 0 ( Δ t ) 3 6 ( Δ t ) 2 2 Δ t 0 0 0 0 0 0 0 0 0 ( Δ t ) 5 20 ( Δ t ) 4 8 ( Δ t ) 3 6 0 0 0 0 0 0 ( Δ t ) 4 8 ( Δ t ) 3 3 ( Δ t ) 2 2 0 0 0 0 0 0 ( Δ t ) 3 6 ( Δ t ) 2 2 Δ t 0 0 0 0 0 0 0 0 0 ( Δ t ) 5 20 ( Δ t ) 4 8 ( Δ t ) 3 6 0 0 0 0 0 0 ( Δ t ) 4 8 ( Δ t ) 3 3 ( Δ t ) 2 2 0 0 0 0 0 0 ( Δ t ) 3 6 ( Δ t ) 2 2 Δ t ]
The initial spectral density, q a c c for the acceleration states can be determined using the initial kinematic solutions by computing the standard deviation of the time-differenced accelerations [14]. In this study, q a c c is set to 1 m2/s4/Hz for the data processing through the analysis of the spectral density from the estimated acceleration data.
As mentioned already, the zenith wet delay residual is modeled with a random walk process so that the corresponding process noise covariance matrix can be expressed as following:
Q Δ Z W D = q Δ Z W D Δ t
It should be noted that the spectral density, q Δ Z W D may be assigned with different values for reference stations and aircraft because the behavior of the zenith wet delay residuals are expected to be different from each other. The spectral density for zenith wet delay, q Δ Z W D are determined by using the empirical auto-covariance functions method proposed by [4]. Then, 3 × 10−10 m2/Hz and 5 × 10−8 m2/Hz are chosen for the spectral densities for reference stations and aircraft, respectively. The process noise covariance matrix for DD ambiguities for L1 and wide-lane is modeled with zero matrix. The covariance matrix for DD ionospheric delay is provided by [3] as shown in Equation (12).
Q i o n = σ 2 ( 1 e 2 ( | Δ t | / T + | δ | / D ) )
where T = 74 min is correlation time, D = 1500 km is correlation distance; σ 2 = 2.0 m2 is variance.

2.2. GPS Measurement Equations

The DD GPS measurements are widely used for the relative positioning due to the fact that various systematic errors can be eliminated by differencing operation [16]. In other words, many error sources in the GPS measurements, for example, receiver and satellite clock errors, and atmospheric effects can be removed or significantly reduced in the DD measurements. However, the spatially correlated atmospheric errors should be modeled as the spatial separation between the two receivers increases. The Equation (13) shows the DD measurement equations for both L1 and L2 frequencies [17,18].
Φ i j , 1 k l = ρ i j k l + T i j k l I i j k l f 1 2 + λ 1 N i j , 1 k l + ε i j , 1 k l Φ i j , 2 k l = ρ i j k l + T i j k l I i j k l f 2 2 + λ 2 N i j , 2 k l + ε i j , 2 k l P i j , 1 k l = ρ i j k l + T i j k l + I i j k l f 1 2 + e i j , 1 k l P i j , 1 k l = ρ i j k l + T i j k l + I i j k l f 2 2 + e i j , 2 k l
where, subscript i and j indicate receivers while superscript k and denote satellites; ( ) i j k is double-differencing operator defined as ( ) i j k = ( ) i k ( ) j k ( ) i + ( ) j ; Φ 1 , Φ 2 are phase pseudo-range measurements on L1 and L2, respectively; P 1 , P 2 are code pseudo-range measurements on L1 and L2, respectively; ρ is the geometric distance between the receiver and satellite; T is tropospheric delay, f 1 and f 2 are carrier frequencies of L1 and L2, respectively; I / f 1 2 and I / f 2 2 are ionospheric delays for L1 and L2, respectively; λ 1 and λ 2 are wavelengths of L1 and L2, respectively; N 1 and N 2 are ambiguities associated with L1 and L2, respectively. ε 1 , ε 2 , e 1 , e 2 are measurement noises.
The DD GPS measurement equations shown in Equation (13) are nonlinear forms; thus, the linearization step is required. After linearization of Equation (13), the design matrix, H has the matrix form as shown in Equation (14).
H = [ h x i h y i h z i h x j 0 0 h y j 0 0 h z j 0 0 h Δ Z W D i h Δ Z W D j λ 1 0 1 h x i h y i h z i h x j 0 0 h y j 0 0 h z j 0 0 h Δ Z W D i h Δ Z W D j λ 2 λ 2 f 1 2 / f 2 2 h x i h y i h z i h x j 0 0 h y j 0 0 h z j 0 0 h Δ Z W D i h Δ Z W D j 0 0 1 h x i h y i h z i h x j 0 0 h y j 0 0 h z j 0 0 h Δ Z W D i h Δ Z W D j 0 0 f 1 2 / f 2 2 ]
As shown in Equation (14), the components of the design matrix can be categorized into partial derivatives with respect to the coordinates for i (reference station) and j (aircraft), ZWD residual, DD ambiguities, and DD ionospheric delays.
Hence, the components corresponding to the coordinates of reference station and rover in the design matrix can be expressed as following:
h x i = x k x i , o ρ i , 0 k + x x i , o ρ i , 0 , h x j = x k x j , o ρ j , 0 k x x j , o ρ j , 0 h y i = y k y i , o ρ i , 0 k + y x i , o ρ i , 0 , h y j = y k y j , o ρ j , 0 k y y j , o ρ j , 0 h z i = z k z i , o ρ i , 0 k + z z i , o ρ i , 0 , h z j = z k z j , o ρ j , 0 k z z j , o ρ j , 0
The tropospheric delay consists of dry and wet components. Approximately 90% of the total zenith delay (TZD) is from the dry component, and about 10% is from the wet component [16]. It is well known that the dry component can be modeled with high accuracy while the wet component is less predictable. Therefore, the ZWD residual with respect to a priori value is estimated in this approach. Equation (16) shows the tropospheric delays between the receiver and satellite.
T i k = Z D D i m i , D k + ( Z W D i + Δ Z W D i ) m i , W k T i = Z D D i m i , D + ( Z W D i + Δ Z W D i ) m i , W T j k = Z D D j m j , D k + ( Z W D j + Δ Z W D j ) m j , W k T j = Z D D j m j , D + ( Z W D j + Δ Z W D j ) m j , W
where Z D D is modeled zenith dry delay; Z W D is modeled zenith wet delay; Δ Z W D is ZWD residual; m is mapping function; subscript D and W indicate “dry” and “wet”, respectively.
The modeled Z D D , Z W D and their mapping functions are computed by using the Saastamoinen model [19]. The tropospheric delays in Equation (16) can be used to make the DD forms and the components of the design matrix correspond to the Δ Z W D is as following:
h Δ Z W D i = m ( z i k ) m ( z i ) h Δ Z W D j = m ( z j k ) + m ( z j )
To obtain the covariance matrix for the DD measurements, R k , the error propagation law is applied to the covariance matrix for the un-difference (UD) measurements. The magnitudes of the noises for code and phase pseudo-range observations are set to decimetre-level and milimetre-level, respectively [14].

3. Results and Discussion

An airborne gravity survey was conducted in South Korea to develop a new geoid model in 2009. A Cessna Grand Caravan was used for the survey and it was flown with the help of autopilot at a constant altitude, i.e., 10,000 feet, in order to get as smooth a flight as possible. The GPS data were collected from both the GPS receiver on board and six Continuously Operating Reference Station (CORS) with 1 s data interval. Figure 1 shows one of the trajectories of aircraft (blue line) and the locations of reference CORS used to demonstrate the proposed algorithm. The trajectories of the aircraft started and ended at the Gimpo (GIMP) airport and the total length of the aircraft trajectory is about 1200 km.
Figure 1. Ground track of aircraft and locations of CORS.
Figure 1. Ground track of aircraft and locations of CORS.
Sensors 15 16895 g001
The coordinates of the reference stations in ITRF2005 frame are determined using the BERNESE GPS data processing software [20], and the results are shown in Table 1. The coordinates of reference stations are tightly constrained by assigning very small variance values (i.e., 1 × 10−10 m2) to the initial covariance matrix.
For the medium to long baseline data processing scenario, the SHAO reference station is included in the data processing. Hence, the diameter of the CORS GPS network reaches up to approximately 900 km. The GPS data was collected on 11 January 2009 and the data span was 4 h and 18 min.
Table 1. ITRF2000 Cartesian coordinates for continuously operating reference station (CORS) (unit: m).
Table 1. ITRF2000 Cartesian coordinates for continuously operating reference station (CORS) (unit: m).
Station NameXYZ
GIMP−3031,894.8284054,598.9883866,287.717
SEOS−3042,060.3694111,978.7573797,578.729
JUNJ−3124,886.9194126,580.5263714,170.141
KWNJ−3134,404.4854173,081.8273654,100.961
CHJU−3168,622.3054277,489.5653501,650.004
SHAO−2831,733.4444675,665.9533275,369.395
As described in Section 2, the kinematic accelerations together with positions and velocities are estimated through the Kalman filter approach. The multi-baseline data processing is performed with radial type of network configuration. To evaluate the filter efficiency, the convergence of the covariance matrix Pk is examined. Figure 2a presents the variations of the trace of Pk with respect to time. As shown in Figure 2a, the trace of Pk converges rapidly with respect to time, which indicates the reliability of the filter. It is also notable that two peaks in Figure 2a correspond to the epochs at which new satellites are observed. Then, new states for DD ionospheric delays and ambiguities are included in state vector, and predefined variance values, i.e., (0.5 m)2 for DD ionospheric delay and 1 × 1010 for DD ambiguity, are assigned to the new states.
Figure 2. Convergence of the covariance matrix, Pk: (a) whole trajectory; (b) beginning part of the trejectory.
Figure 2. Convergence of the covariance matrix, Pk: (a) whole trajectory; (b) beginning part of the trejectory.
Sensors 15 16895 g002
Figure 3a–c present the estimated positions, velocities and accelerations of the aircraft with respect to time, respectively. The estimated states are the results from forward-backward smoothing.
Figure 3. (a) Positions; (b) velocities; (c) accelerations.
Figure 3. (a) Positions; (b) velocities; (c) accelerations.
Sensors 15 16895 g003
As shown in Figure 3b, the constant velocities in most of the data span for airborne gravity survey are retained. The peaks observed in Figure 3c correspond to the aircraft maneuvering for takeoff, turnaround, and landing. Table 2 shows the statistics of the estimated accelerations which ranges from about −4.5 m/s2 to 4.5 m/s2.
Table 2. Statistical characteristics of estimated accelerations (unit: m/s2).
Table 2. Statistical characteristics of estimated accelerations (unit: m/s2).
x ¨ y ¨ z ¨
Min.−3.92−4.18−4.62
Max.3.754.344.49
Std.0.530.510.52
The atmospheric effects are estimated in forms of ZWD residuals and DD ionospheric delays, respectively. As explained in Section 2, the ZWD residuals are estimated at each station including the aircraft and then final TZDs are computed by adding the modeled values to the estimated ZWD residuals. Figure 4a shows the final TZDs obtained from the data processing. As shown in Figure 4a, no significant variations of TZDs for the reference stations are observed while abrupt changes in the magnitude of values can be seen in the TZD estimated from the aircraft. This is caused by the fact that the modeled values are highly depended on the altitude of the aircraft. The DD ionospheric delays are also computed for each of the baseline and satellite pairs, and one example of the results for SHAO-AIRP baseline is shown in Figure 4b. Each continuous line corresponds to the DD ionospheric delay for each pair of GPS satellites. The variation of the DD ionospheric delays ranges from about −1 m to 1 m.
Figure 4. (a) Total zenith delays (TZDs); (b) DD ionospheric delays for SHAO-AIRP baseline.
Figure 4. (a) Total zenith delays (TZDs); (b) DD ionospheric delays for SHAO-AIRP baseline.
Sensors 15 16895 g004
The next step is to evaluate the quality of the estimated kinematic accelerations once all the solutions are obtained through the proposed algorithm. However, there is a limitation on direct comparison because the reference values for kinematic accelerations are not available. Therefore, the coordinates of aircraft are computed using the same GPS datasets with the commercial software called Waypoint®, and the computed coordinates are used as reference values for the evaluation. Figure 5 presents the differences between the estimated coordinates using the proposed method and the reference values from the Waypoint software. The differences between the two solutions range from about −0.15 m to 0.1 m as shown in Table 3. From the results, one can conclude that the filter based on the proposed algorithm is working properly.
Figure 5. Differences between the estimated coordinates and reference values.
Figure 5. Differences between the estimated coordinates and reference values.
Sensors 15 16895 g005
Table 3. Statistical characteristics of the differences between the estimated coordinates and the reference values (unit: m).
Table 3. Statistical characteristics of the differences between the estimated coordinates and the reference values (unit: m).
Δ x Δ y Δ z
Min.−0.09−0.17−0.17
Max.0.110.110.02
Std.0.040.060.04
Then, the second-order time-derivatives of determined positions are applied to obtain the kinematic accelerations. It should be mentioned that taking the time-derivatives is performed after the data fitting with B-spline function. In the final step, the differences between the two kinematic accelerations are computed as shown in Figure 6. This procedure is applied to only 1 h data span, i.e., 1–2 h (UTC) because no significant dynamics of the aircraft is observed. This also indicates that the actual gravity survey is performed during that time span.
Table 4 shows the statistical characteristics of the differences between the estimated kinematic accelerations and the reference values. The differences between the two solutions range from −0.15 m/s2 to 0.16 m/s2 and the standard deviation of the differences is about 0.03 m/s2. The mean value of the differences is 10−5 m/s2 level for each component.
From the results, one can state that the comparable kinematic accelerations for airborne gravity survey can be obtained by using the proposed algorithmIt should be mentioned that the kinematic acceleration for airborne gravimetry needs the accuracy at least at the level of 10−4 m/s2 to detect the gravity signal at the level of 10 mGal. In airborne gravimetry, however, a smoother is applied to the acceleration to extract meaningful gravity signal. Although the extraction of the gravity signal is out of the scope in this study, a B-spline smoother with window size of 120 s are applied to both accelerations and verified that the differences reside at the level of 10−4 m/s2. Therefore, one can conclude that the proposed algorithm is reliable for the application of kinematic accelerations. The validation for the proposed algorithm in this study, however, is performed with the dataset of one trajectory collected under the relatively low dynamic conditions. Therefore, it might be necessary to process more dataset under the high dynamic environments in the future.
Figure 6. Differences between the estimated kinematic accelerations and reference values.
Figure 6. Differences between the estimated kinematic accelerations and reference values.
Sensors 15 16895 g006
Table 4. Statistical characteristics of the differences between the estimated kinematic accelerations and the reference values (unit: m/s2).
Table 4. Statistical characteristics of the differences between the estimated kinematic accelerations and the reference values (unit: m/s2).
Δ x ¨ Δ y ¨ Δ z ¨
Min.−0.15−0.13−0.14
Max.0.130.160.13
Std.0.030.030.02

4. Conclusions

In this study, we proposed the kinematic GPS positioning algorithm using multiple reference stations for medium to large scale networks. The Kalman filter with PVA dynamic model is used for the kinematic positioning so that the kinematic acceleration information can be obtained simultaneously. For the long range kinematic applications, the tropospheric and ionospheric delays are modeled with random walk and first-order Gauss-Markov process models, respectively. The algorithm is implemented and tested with in situ airborne GPS data collected at 10,000 ft altitude with 1 s data interval. For the validation of the proposed algorithm, the positioning results are compared with the reference values first. The computed differences between the estimated and reference values range from −0.17 m to 0.11 m, which indicates the validity of the proposed algorithm. The estimated kinematic accelerations from the proposed algorithms are also compared with those from the second-order time-derivatives of the reference positions. The computed differences between the two solutions range from −0.15 m/s2 to 0.16 m/s2 and the standard deviation of the differences is about 0.03 m/s2. From these results, we can conclude that comparable kinematic accelerations, which can be used for airborne gravimetry, are obtainable using the proposed algorithm.

Acknowledgments

This work was supported by the DGIST R&D Program of the Ministry of Science, ICT & Future Planning of Korea (15-IT-01).

Author Contributions

Chang-Ki Hong carried out the development and implementation of the proposed algorithms. Chi Ho Park and Joong-Hee Han performed numerical test and data analysis. Jay Hyoun Kwon conducted all of the research and wrote the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cannon, M.E.; Lachapelle, G.; Alves, P.; Fortes, L.P.; Townsend, B. GPS RTK positioning using a regional reference network: Theory and results. In Proceedings of the 5th GNSS international symposium, Seville, Spain, 8–11 May 2001.
  2. Wanninger, L. Virtual Reference Stations for Centimeter-Level Kinematic Positioning. In Proceedings of the 15th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 2002), Portland, OR, USA, 24–27 September 2002; pp. 1400–1407.
  3. Yang, M.; Tang, C.H.; Yu, T.T. Development and assessment of a medium-range real-time kinematic GPS algorithm using an ionospheric information filter. Earth Planets Space 2000, 52, 783–788. [Google Scholar] [CrossRef]
  4. Schüler, T. On Ground-Based GPS Tropospheric Delay Estimation; Univ. der Bundeswehr: München, Germany, 2001. [Google Scholar]
  5. Dekkiche, H.; Kahlouche, S.; Abbas, H. Differential ionosphere modelling for single-reference long-baseline GPS kinematic positioning. Earth Planets Space 2010, 62, 915–922. [Google Scholar] [CrossRef]
  6. Brown, R.G.; Hwang, P.Y.C. Introduction to Random Signals and Applied Kalman Filtering, 3rd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1997. [Google Scholar]
  7. Jekeli, C.; Garcia, R. GPS phase accelerations for moving-base vector gravimetry. J. Geod. 1997, 71, 630–639. [Google Scholar] [CrossRef]
  8. Jekeli, C. Inertial Navigation Systems with Geodetic Applications; Walter de Gruyter: New York, NY, USA, 2001. [Google Scholar]
  9. Li, X. Moving Base INS/GPS Vector Gravimetry on a Land Vehicle. Ph.D. Thesis, The Ohio State University, Columbus, OH, USA, 2007. [Google Scholar]
  10. Kwon, J.H. Airborne Vector Gravimetry Using GPS/INS; Dept. of Geodetic Science and Surveying, The Ohio State University: Columbus, OH, USA, 2000. [Google Scholar]
  11. Bogatin, S.; Foppe, K.; Wasmeier, P.; Wunderlich, T.A.; Schäfer, T.; Kogoj, D. Evaluation of linear Kalman filter processing geodetic kinematic measurements. Measurement 2008, 41, 561–578. [Google Scholar] [CrossRef]
  12. Gelb, A. Applied Optimal Estimation; The M.I.T. Press: London, UK, 1974. [Google Scholar]
  13. Simon, D. Optimal State Estimation: Kalman, H Infinity, and Nonlinear Approaches; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2006. [Google Scholar]
  14. Petovello, M.G. Real-Time Integration of a Tactical-Grade IMU and GPS for High-Accuracy Positioning and Navigation. Ph.D. Thesis, Dept. of Geomatics Engineering, University of Calgary, Calgary, AB, Canada, 2003. [Google Scholar]
  15. Shin, E.H. Estimation Techniques for Low-Cost Inertial Navigation. Ph.D. Thesis, Dept. of Geomatics Engineering, University of Calgary, Calgary, AB, Canada, 2005. [Google Scholar]
  16. Hofmann-Wellenhof, B.; Lichtenegger, H.; Collins, J. GPS Theory and Practice, 5th ed.; Springer: Wien, Austria, 2001. [Google Scholar]
  17. Grejner-Brzezinska, D.A. Analysis of GPS Data Processing Techniques: In Search of Optimized Strategy of Orbit and Earth Rotation Parameter Recovery; Dept. of Geodetic Science and Surveying, The Ohio State University: Columbus, OH, USA, 1996. [Google Scholar]
  18. Goad, C.C.; Yang, M.A. New Approach to Precision Airborne GPS Positioning for Photogrammetry. Photogramm. Eng. Remote Sens. 1997, 63, 1067–1077. [Google Scholar]
  19. Saastamoinen, J. Atmospheric correction for the troposphere and stratosphere in radio ranging satellites. Geophys. Monogr. Ser. 1972, 15, 247–251. [Google Scholar]
  20. Dach, R.; Hugentobler, U.; Fridez, P.; Meindl, M. Bernese GPS Software Version 5.0; Astron. Institute, Univ. Berne: Berne, Swizerland, 2001. [Google Scholar]

Share and Cite

MDPI and ACS Style

Hong, C.-K.; Park, C.H.; Han, J.-h.; Kwon, J.H. Medium to Long Range Kinematic GPS Positioning with Position-Velocity-Acceleration Model Using Multiple Reference Stations. Sensors 2015, 15, 16895-16909. https://doi.org/10.3390/s150716895

AMA Style

Hong C-K, Park CH, Han J-h, Kwon JH. Medium to Long Range Kinematic GPS Positioning with Position-Velocity-Acceleration Model Using Multiple Reference Stations. Sensors. 2015; 15(7):16895-16909. https://doi.org/10.3390/s150716895

Chicago/Turabian Style

Hong, Chang-Ki, Chi Ho Park, Joong-hee Han, and Jay Hyoun Kwon. 2015. "Medium to Long Range Kinematic GPS Positioning with Position-Velocity-Acceleration Model Using Multiple Reference Stations" Sensors 15, no. 7: 16895-16909. https://doi.org/10.3390/s150716895

APA Style

Hong, C. -K., Park, C. H., Han, J. -h., & Kwon, J. H. (2015). Medium to Long Range Kinematic GPS Positioning with Position-Velocity-Acceleration Model Using Multiple Reference Stations. Sensors, 15(7), 16895-16909. https://doi.org/10.3390/s150716895

Article Metrics

Back to TopTop