Next Article in Journal
Advancing Learning Assignments in Remote Sensing of the Environment Through Simulation Games
Previous Article in Journal
An Atmospheric Correction Method over Bright and Stable Surfaces for Moderate to High Spatial-Resolution Optical Remotely Sensed Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Kalman Filtering Based on Chi-square Increment and Its Application

1
Key Laboratory of Geographic Information Science, Ministry of Education, East China Normal University, Shanghai 200241, China
2
Engineering Center of SHMEC for Space Information and GNSS, School of Communication and Electronic Engineering, East China Normal University, Shanghai 200241, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(4), 732; https://doi.org/10.3390/rs12040732
Submission received: 29 November 2019 / Revised: 17 January 2020 / Accepted: 20 February 2020 / Published: 22 February 2020
(This article belongs to the Section Urban Remote Sensing)

Abstract

:
In Global Navigation Satellite System (GNSS) positioning, gross errors seriously limit the validity of Kalman filtering and make the final positioning solutions untrustworthy. Thus, the detection and correction of gross errors have become indispensable parts of Kalman filtering. Starting by defining an incremental Chi-square method of recursive least squares, this paper extends this definition to Kalman filtering to detect gross errors, explains its nature and its relation with the currently adopted Chi-square variables of Kalman filtering in model and data spaces, and compares them with the predictive residual statistics. Two robust Kalman filtering models based on an incremental Chi-square method (CI-RKF) were established, and the one with a better incremental Chi-square component was selected based on a static accuracy evaluation experiment. We applied the selected robust model to the GNSS positioning and the GNSS/inertial measurement unit (IMU)/visual odometry (VO) integrated navigation experiment in an occluded urban area at the East China Normal University. We compared the results for conventional Kalman filtering (CKF) with a robust Kalman filtering constructed using predictive residual statistics and an Institute of Geodesy and Geophysics (IGGШ) weight factor, abbreviated as “PRS-IGG-RKF”. The results show that the overall accuracy of CI-RKF in GNSS positioning was improved by 22.68%, 54.33%, and 72.45% in the static experiment, and 12.30%, 7.50%, and 16.15% in the kinematic experiment. The integrated navigation results indicate that the CI-RKF fusion method increased the system positioning accuracy by 66.73%, 59.59%, and 59.62% in one of the severe occlusion areas, and 42.04%, 59.04%, and 52.41% in the other.

1. Introduction

Kalman filtering is a recursive algorithm widely used in positioning and navigation. It dynamically updates the state variables of a system and recursively estimates the state variables by assigning proper weights to the observations and state predictions. The weights of conventional Kalman filtering are defined based on the observation signal-to-noise ratio and satellite elevation angles [1]. When the real observations do not satisfy the mathematical model or the statistical properties of the measurement noise and the dynamic state process (in particular, when signals undergo severe interference from multipath, no-line-of-sight (NLOS) and gross errors), such a recursive process of conventional Kalman filtering not only generates a poor solution for this epoch, but also propagates poor quality solutions for subsequent epochs. To remedy such a vulnerability, scholars have proposed various methods to mitigate the impact of these abnormal situations, such as function model compensation [2], adaptive Kalman filtering [3], and robust Kalman filtering [4]. Both the function model compensation and adaptive Kalman filtering methods primarily examine the reasonableness of the underlying mathematical model and noise statistics. When a violation is detected, these methods revise the original inaccurate dynamic model and process noise statistics to improve the performance of the recursive filtering. These methods usually consume extra computational resources to detect abnormal observations and gross errors, and hence are not effective or conducive for civil real-time positioning in cases where the impact of gross errors is overwhelming. On the other hand, robust Kalman filtering is very effective for gross error detection. This paper focuses on the detection and elimination of abnormal observations and gross errors of Kalman filtering, in order to ensure its stability and reliability. Hence, for the rest of the paper our attention is confined to robust Kalman filtering.
Current robust algorithms are mainly composed of two parts: gross error detection and robust estimation. The main idea of these robust Kalman filtering algorithms is to use learning statistics for judging errors [5] to detect gross errors in observations, then eliminate them or assign proper reduced weights to minimize their contamination [6]. Commonly used algorithms such as the Huber scheme [7], Institute of Geodesy and Geophysics (IGGIII) scheme [8], two-factor equivalent weights [9], robust Bayesian estimation [10], Danish function [11], and so on all rely on the learning statistics from the discrepancy between the predictive residuals and posterior residuals of observations [12]. A commonly adopted approach is to construct a Chi-square variable with a Chi-square test to detect gross errors and monitor system failure. One type of Chi-square variable is based on the departure between the propagated state vector and the state vector estimated from Kalman filtering; its corresponding statistical test is called the state Chi-square test (SCST) [13,14]. Another type of Chi-square variable is based on the residuals between observations and predicted values from the estimated state; its corresponding statistic test is called the innovation Chi-square test (ICST) [15]. In ordinary least squares (LS), the Chi-square variable is defined based on the Sum of Squares of Errors (SSE) [16]. In particular, the increment of Chi-square ( Δ χ 2 ) between two adjacent epochs is sensitive to abnormal observations and gross errors at this epoch [17]. The concept of the increment of Δ χ 2 also can be extended to Kalman filtering [18]. However, the connection between the Δ χ 2 extending from the recursive LS (RLS) and Chi-square variables of Kalman filtering that are defined directly in the model and data spaces has not yet been discussed. In this paper, we reveal that the Δ χ 2 extending from the RLS is equivalent to the ICST and, hence, that its contributors are from both model and data spaces. We then deploy the Δ χ 2 -based robust Kalman filtering in real observation scenarios by defining the corresponding robust weight function.
The emerging integrated positioning of multi-constellation satellite navigation systems enlarges spatio-temporal coverage while enhancing positioning accuracy and stability. Meanwhile, it increases the probability of observation anomalies and blunders [19]. In occluded urban areas, the ratio of abnormal gross errors becomes much higher due to complex multipath environments and limited sky visibility [20]. Thus, the gross error detection and on-the-spot contamination reduction appear critical for real-time Global Navigation Satellite System (GNSS) positioning, an inertial measurement unit (IMU), and a visual odometry (VO)-assisted GNSS (GNSS/IMU/VO) integrated navigation system [21]. We tested our Δ χ 2 -based robust Kalman filtering scheme in real-time static and kinematic experiments of both GNSS positioning and GNSS/IMU/VO integrated navigation, analyzed its positioning accuracy improvement, and compared its performance with the robust Kalman filtering based on predictive residual statistics. The observations were obtained from instrumental measurements, and their contents were related to their practical applications, such as the observation results from various sensors [22]. In this paper, the observations represent the pseudo-range measurements from the receiver in GNSS positioning, and refer to the positioning parameters obtained by GPS, IMU, and VO in integrated navigation. The validity of this model verifies that Δ χ 2 -based robust Kalman filtering does improve the accuracy and stability of real-time kinematic positioning, and has an advantage over the other schemes of Kalman filtering that are compared in the paper.

2. Methodology

2.1. Gross Error Detection

The SSE is the sum of squared deviations of observed values from predicted values. It stems from the regression analysis and is used to measure discrepancies between observed data and an estimation model. It is also an important statistic in LS positioning [16], which is equal (in numerals) to the weighted sum of the squares of errors (WSSE) [23]. Since the LS positioning solution minimizes the WSSE components [24], the magnitude of the SSE of the LS reflects the degree of consistency between the observed values, and can be used to evaluate the fitting goodness of the estimated parameters to the actual observations.
Differently from the ordinary LS, where the model prediction uses full observations, the recursive LS (RLS) builds its model based on a subset of observations (called as “past history”). It has been shown that such derived recursive residuals are very sensitive to such a point that the model is not supported by real data [25]. Later, it was suggested that the recursive residuals were more sensitive to data outliers than the residuals from the ordinary LS [26]. Scientists have proposed the routine calculation of recursive residuals in an RLS analysis to diagnose both model and data outliers [27].
The observation equation at epoch k is expressed in forms of vector and matrix:
z k = h k x k + ν k
where z k , h k , ν k are the observation vector, design matrix, and noise vector measurements, respectively, and their subscript k represents the kth epoch; the observation error covariance matrix at this epoch is r k ; x k denotes the state vector at epoch k. In our positioning applications, we assume that observations and the error covariance matrix at each epoch are independent, and that the design matrix has a full ranking. After this is established, the solution and its covariance matrix can be expressed in a recursive form [28]:
x ^ k = x ^ k 1 + C x ^ k 1 h k T ( h k C x ^ k 1 h k T + r k ) 1 ( z k h k x ^ k 1 ) C x ^ k = C x ^ k 1 C x ^ k 1 h k T ( h k C x ^ k 1 h k T + r k ) 1 h k C x ^ k 1
where the symbol ^ is the estimated solution vector and the superscript T denotes matrix transposal. The estimated solution covariance matrix at epoch k is represented by C.
The SSE (Chi-square) of such a recursive LS from epoch 1 to epoch k is defined as
χ k 2 = i = 1 k ( z i h i x ^ k ) r i 1 ( z i h i x ^ k )
The increment of Chi-square at epoch k is defined as
Δ χ k 2 = χ k 2 χ k 1 2
It was proven that Δ χ k 2 can also be expressed in a recursive form [18]:
Δ χ k 2 = Δ x ^ k T C x ^ k 1 1 Δ x ^ k + ( z k h k x ^ k ) T r k 1 ( z k h k x ^ k ) Δ x ^ k = x ^ k x ^ k 1
This concept of Δ χ k 2 from RLS can be extended to Kalman filtering. Suppose that the state equation and observation equation of Kalman filtering at epoch k are [29,30]
x k = A k x k 1 + B k u k 1 + Γ k ω k
z k = h k x k + v k
where subscript k still represents the epoch index; x k is the state vector; u k is the control function; z k and h k represent the observation vector and design matrix; A k ,   B k are the transition coefficient matrix of the state and the control function transition; Γ k is the input matrix of the state disturbance noise; ω k ~ N ( 0 , q ) is the state disturbance noise; v k ~ N ( 0 , r ) is the observation noise. Suppose that ω k ,   v k are independent of each other and obey a normal distribution of white noise.
The state transition and update equations are as follows:
x ^ k , k 1 = A k x ^ k 1 + B k u k 1 C x ^ k , k 1 = A k C x ^ k 1 A k T + Γ k q k Γ k T K k = C x ^ k , k 1 h k T ( h k C x ^ k , k 1 h k T + r k ) 1 x ^ k = x ^ k , k 1 + K k ( z k h k x ^ k , k 1 ) C x ^ k = ( I K k h k ) C x ^ k , k 1
where x ^ k represents the estimated state variable at epoch k and x ^ k , k 1 denotes the predicted state vector at epoch k from epoch k − 1; I is the identity matrix and K k is the Kalman gain matrix at epoch k. According to the Δ χ k 2 of the recursive LS, the Δ χ k 2 of the k epoch in Kalman filtering is defined as [31]
Δ χ k 2 = Δ x ^ k T C x ^ k , k 1 1 Δ x ^ k + ( z k h k x ^ k ) T r k 1 ( z k h k x ^ k )
where Δ x ^ k = x ^ k x ^ k , k 1 . It can be seen from the above equation that the Δ χ k 2 of the k epoch is composed of two parts, with the first part being the increment caused by the change of fitting solution, and the second part being the increment caused by the observation residual of the k epoch.
The Δ χ k 2 derived from Equation (9) acts as an indicator of gross error detection. Although it demonstrates its nature clearly, it is inconvenient in practical use because it requires an updated state variable x ^ k . Once a gross error occurs, it indicates that the solution and its covariance matrix of the Kalman filtering have already been distorted and can no longer be used recursively to get the correct solution of the next epoch [32]. In order to obtain a normal solution, the solution and its covariance matrix of the previous epoch must first be saved. If no gross error is found, the solution and covariance matrix will be updated. Otherwise, it is necessary to remove the gross error of the k epoch and restore the previously saved uncontaminated solution and its covariance matrix, then continue the operation. Such a procedure not only occupies more computational memory, but also reduces the operation speed greatly, which is not conducive to real-time positioning. To obtain an estimated value of Δ χ k 2 before updating the solution, we express the Δ χ k 2 and the observation residuals at k epoch by the form of their propagated values from the k − 1 epoch:
Δ x ^ k = C x ^ k , k 1 h k T ( h k C x ^ k , k 1 h k T + r k ) 1 ( z k h k x ^ k , k 1 )
z k h k x ^ k = r k ( h k C x ^ k , k 1 h k T + r k ) 1 ( z k h k x ^ k , k 1 )
Bring Equations (10) and (11) into Equation (9):
Δ χ k 2 = ( z k h k x ^ k , k 1 ) T ( h k C x ^ k , k 1 h k T + r k ) 1 ( z k h k x ^ k , k 1 )
Equation (12) is equivalent to Equation (9), but Equation (9) uses the updated x ^ k , while Equation (12) uses the propagated x ^ k , k 1 and C x ^ k , k 1 . Therefore, Formula (12) can judge whether the epoch is normal or not through the prediction stage only. If the raw observation causes the Δ χ k 2 anomaly, it will be regarded as a gross error, which can be further processed by eliminating data or reducing the Kalman filtering weight to get the normal positioning results.
For comparison, the SCST-adopted Chi-square variable is defined as [14]
Δ χ k 2 = Δ x ^ k T ( C x ^ k , k 1 C x ^ k ) 1 Δ x ^ k
which is different from our Equation (9) and is confined to state variables only in a parameter space. Meanwhile the ICST-adopted Chi-square variable is defined as [15]
Δ χ k 2 = ( z k h k x ^ k , k 1 ) T ( h k C x ^ k , k 1 h k T + r k ) 1 ( z k h k x ^ k , k 1 )
which is exactly the same as Equation (12). Thus, the nature of the ICST-adopted Chi-square variable is equivalent to the extension of the RLS-defined Δ χ k 2 , which actually covers the variance caused by the new data in both the data and model spaces according to Equation (9). The matrix h k C x ^ k , k 1 h k T + r k is called the covariance matrix of innovation, and is introduced by the new observation.
The learning statistics for judging errors, such as predictive residual statistics, are express as [33,34,35]
Δ V ^ k = ( v ^ k T v ^ k t r ( C v ^ k ) ) 1 2
where v ^ k = h k x ^ k , k 1 z k is the predictive residual vector, and C v ^ k = h k C x ^ k , k 1 h k T + r k is the covariance matrix of innovation. Equation (14) can be converted to
Δ χ ^ k 2 = v ^ k T C ν ^ k 1 v ^ k
A comparison shows that Equation (15) is a simplified version of Equation (16), which uses a trace of the covariance matrix of innovation to replace the full covariance matrix. Hence, Equation (16) is more rigorous, as all correlations among the predictive residuals are considered. Therefore, compared with learning statistics, Δ χ 2 reflects the impact of gross errors more accurately.

2.2. Robust Kalman Filtering Based on Chi-Square Increment

Ideally, z k h k x ^ k , k 1 should obey the normal distribution with a zero mean value [36]:
z k h k x ^ k , k 1 ~ N ( 0 , h k C x ^ k , k 1 h k T + r k )
When a gross error occurs, the predictive residual no longer obeys the above normal distribution. The Chi-square test for Δ χ k 2 provides useful statistics to identify gross errors in observation.
Method 1 is used to test the components of the vector z k h k x ^ k , k 1 totally:
( z k h k x ^ k , k 1 ) T ( h k C x ^ k , k 1 h k T + r k ) 1 ( z k h k x ^ k , k 1 ) ~ χ 2 ( n )
Set the following hypothesis test:
{ H 0 : Δ χ k 2 ~ χ 2 ( n , 0 ) H 1 : Δ χ k 2 ~ χ 2 ( n , λ )
where H0 denotes that there is no gross error in the system and H1 denotes that there is a gross error in the system; n is the degree of freedom of Chi-square distribution; and λ is a centralized parameter matrix.
For observations with gross errors, robust estimation is usually achieved by enlarging the observation covariance noise matrix [37,38,39]. According to the above hypothesis, the robust factors are constructed with Δ χ k 2 :
β k = { 1 | Δ χ k 2 χ α 2 ( n ) | < c 0 | Δ χ k 2 χ α 2 ( n ) | c 0 < | Δ χ k 2 χ α 2 ( n ) | < c 1 | Δ χ k 2 χ α 2 ( n ) | 2 | Δ χ k 2 χ α 2 ( n ) | > c 1
where α is the confidence level, and c 0 ,   c 1 are two pre-defined constants. Then, the observation covariance matrix is updated to r ¯ k = β k r k .
Method 2 is used to test the components of the vector z k h k x ^ k , k 1 individually.
From Equation (18), it can be deduced that an individual Δ χ k , i 2 should obey the Chi-square distribution with a degree of freedom equal to 1:
( z k h k x ^ k , k 1 ) i T ( h k C x ^ k , k 1 h k T + r k ) i , i 1 ( z k h k x ^ k , k 1 ) i ~ χ 2 ( 1 )
where the subscript i denotes the i-th component. Then, the hypothesis test is as follows:
{ H 0 : Δ χ k , i 2 ~ χ 2 ( 1 , 0 ) H 1 : Δ χ k , i 2 ~ χ 2 ( 1 , λ )
Similarly, the robust factors can be constructed as
β i = { 1 | Δ χ k , i 2 χ α 2 ( 1 ) | < c 0 | Δ χ k , i 2 χ α 2 ( 1 ) | c 0 < | Δ χ k , i 2 χ α 2 ( 1 ) | < c 1 | Δ χ k , i 2 χ α 2 ( 1 ) | 2 | Δ χ k , i 2 χ α 2 ( 1 ) | > c 1
The observation covariance matrix is updated to r ¯ k , i = β i r k , i ,   r ¯ k = [ r ¯ k , 1 r ¯ k , i r ¯ k , m ] .
The updating process of robust Kalman filtering is
r ¯ k = C x ^ k , k 1 h k T ( h k C x ^ k , k 1 h k T + r ¯ k ) 1 x ^ k = x ^ k , k 1 + K ¯ k ( z k h k x ^ k , k 1 ) C x ^ x = ( I K ¯ K h k ) C x ^ k , k 1

2.3. Mathematical Model of CI-RKF in GNSS Positioning

The observation equation of the pseudo-range [39] is
ρ k i = r k i + δ τ k δ t k i + I k i + T k i + ε ρ i
where subscript k still denotes the epoch and superscript i is the observation index or satellite index at each epoch; ρ is the pseudo-range measurement; r is the geometric distance between the satellite i and the receiver; δ τ k is the receiver error; I is the ionosphere delay; T is the atmosphere delay; and ε ρ is the pseudo-range measurement noise. In order to establish the pseudo-range positioning equation, Equation (25) can be simplified as the following:
r k i + δ τ k = ρ ¯ k i ε ρ i
where ρ ¯ k i = ρ k i + δ t k i I k i T k i is the pseudo-range measurement value after correction.
The initial position x ^ 0 is determined by LS positioning [40], Γ = e y e ( 8 ) , and B = z e r o s ( 8 ) , where eye indicates the unit diagonal matrix, zeros indicates the zero matrix, and the number in brackets is the dimension of the matrix. In the static experiment of GNSS positioning, the state transition coefficient matrix is defined as A = [ e y e ( 4 ) d i a g ( [ 0 0 0 T s ] ) z e r o s ( 4 , 4 ) d i a g ( [ 0 0 0 T s ] ) ] , where T s is the sampling interval. In the kinematic experiment, A = [ e y e ( 4 ) T s e y e ( 4 ) z e r o s ( 4 , 4 ) e y e ( 4 ) ] [41].
In Figure 1, all calculations except Δ χ k 2 are the routine Kalman filtering formula. For the gross error checking, the only increased computational burden is the Δ χ k 2 calculation, in which all matrix values are ready from routine analysis and the dimension of the matrix inversion is the observation number at epoch k. Such an additional computation is tolerable for real-time positioning. In our GNSS positioning static experiment, the central processing unit (CPU) time is increased from 90.643 to 93.284 seconds (only a 2.91% hike), then increased from 4.431 to 4.561 seconds in the GNSS positioning kinematic experiment (only a 2.96% hike) due to the gross checking.

2.4. Mathematical Model of CI-RKF in GNSS/IMU/VO Integrated Positioning

An IMU is a commonly used auxiliary navigation device [42]. It has the advantages of being independent from external information, freedom from electromagnetic interference, and good short-term accuracy and stability [43]. However, an IMU requires sporadic external positioning information to set up its initial position and calibrate its solutions in order to avoid cumulative positioning errors. VO is a key module of lidar sensor navigation technology based on simultaneous localization and mapping (SLAM) [44]. It locates its position autonomously through the surrounding environment, independent of satellite signals [45,46]. Due to the image blur, lens blocking, and other issues, it is difficult to perform positioning applications using VO technology on its own. As a result, the combination of an IMU and VO can improve the reliability of the system. An IMU and VO provide the relative location, while GNSS provides the absolute location. This system can provide short-term and relatively accurate position information by using a VO/IMU when the GNSS signal is temporarily missing, and for this reason the GNSS/IMU/VIO has become a trend in current integrated navigation.
The construction of a combined GNSS and VO/IMU is based on the loose coupling method [47] in which the VO and IMU are tightly coupled (the combined process can be found in [48,49]). The velocity information from the VO/IMU and the position information of GNSS are used as inputs [50], and the system positions are estimated via filtering [51]. The system coupling mode is shown in Figure 2.
The state equation and the observation equation of the system are
X k = A X k 1 + Q
Z k = H X k + R
where A = [ I 3 3 T s × I 3 3 0 3 3 I 3 3 ] is the state transition matrix, and I is the identity matrix; H = [ I 6 6 ] is the design matrix; X k = [ x k y k z k v x v y v z ] T is the matrix of six variables to be estimated, x k , y k , and z k are the positional variables; v x , v y , and v z are the velocity variables; Q is the state noise; and R is the observation noise.
The system state prediction and state estimation equations are
X ^ k , k 1 = A X ^ k 1
X ^ k = X ^ k , k 1 + K k ( Z k H k X ^ k , k 1 )
The system status update process is as shown in Equation (24).

3. Data and Experiments

Experiment 1: A corridor between Building A and Building B of the Estuary and Coastal Research Institute in East China Normal University was selected as our platform, as shown in Figure 3. The experimental data were collected by a U-blox9 receiver. The sampling frequency was 1 Hz, and total data covered 5522 epochs. The positioning results from the Trimble R8 were taken as the reference solutions.
Experiment 2: A route was planned around the campus as the kinematic experiment environment, with the occlusion caused mainly by the surrounding trees and buildings. The experimental equipment was placed on the roof of a car to receive the satellite signal (Figure 4). The data were also collected by both U-blox9 and Trimble BD982 receivers with a 1-Hz acquisition frequency. The solutions were derived from the data of the Trimble BD982 receiver and were solved using Inertial Explorer 8.60, a high-precision post-processing software that served as the ground truth. In order to reduce the influence of instability in the attitude and route, the vehicle tried to keep running in the middle of road at a near-constant speed.
Experiment 3: Selecting part of the route from Experiment 2, the vehicle moved around the building- dense area of the campus at a relatively constant speed. A HUAWEI HonorV9 mobile phone was placed in the carriage to collect the raw GPS observations through software named rinexOn [52], and the receiving frequency was set to 1 Hz. The IMU (MIT-G-710 of xsen series) and VO (ZED Stereo Camera) were fixed on the roof of the car to collect the data, and the sampling frequency and frame rate were 100 Hz and 30 fps, respectively. At the same time, the Trimble BD982 receiver was placed on the roof with a receiving frequency of 1 Hz, and its data were processed by Inertial Explorer 8.60. The vehicle’s experimental devices are shown in Figure 5.

4. Results and Discussion

Two robust methods were constructed by a Chi-square increment in Section 2. In this section, we first compare their performances based on the data collected through Experiment 1, select the better one, and record it as robust Kalman filtering based on a Chi-square increment (CI-RKF). The results are shown in Section 4.1. Next, the CI-RKF is compared with conventional Kalman filtering (CKF) and the robust Kalman filtering constructed using predictive residual statistics and IGGШ weight factor (PRS-IGG-RKF) based on the data from Experiments 1 and 2, respectively. The results are discussed in Section 4.2 and Section 4.3. Finally, the three schemes are applied to the GNSS/IMU/VO integrated navigation, and the results are obtained by comparing the fusion positioning accuracy in the shelter zone based on Experiment 3, which is shown in Section 4.4.

4.1. Comparison of Chi-Square Increment Robust Methods

The experiment results of two robust incremental Chi-Square methods are shown in Figure 6. The α is set to 0.15 in both methods, and c 0 ,   c 1 is set to 3.7 and 2.3 in methods 1 and 2, respectively.
It can be seen in Figure 6 that the deviations between the final positioning results and the true values of Method 2 were less than those of Method 1, and the convergence speed of Method 2 was faster than that of Method 1, meaning the detection index in Method 2 was more sensitive to the gross error. Because the vector z k H k x ^ k , k 1 was tested as a whole in Method 1, the entire observed noise covariance matrix will be enlarged if the sum of errors per epoch is greater than the detection threshold, and the magnification of each covariance will also be the same. However, in Method 2, the components in vector z k H k x ^ k , k 1 were tested individually. If a satellite fails, the corresponding noise variance term will be amplified, and the magnification of each variance varies according to the deviation degree of the residual. The final positioning accuracies of Method 2 were 21.02%, 60.69%, and 85.57% higher than that of Method 1 in the east, north and up (E/N/U) direction. Thus, Method 2 was selected for discussion in the following section.

4.2. Comparison of GNSS Positioning Schemes in a Static Experiment

In order to evaluate the quality of the raw observations, the experimental data were solved by the CKF using the incremental Chi-square method from Section 2.1. The change of the Chi-square increment in the static experiment is shown in Figure 7a. The normalized Chi-square increment was obtained by dividing the Chi-square increment by the total number of satellites per unit epoch, as shown in Figure 7b. It can be seen from Figure 8 that there are gross errors in the observation at the 500–600 and the 2000–2500 epochs. In Figure 3, the experimental equipment was placed in the area similar to the urban canyon [53], so that as the satellites were moving, some signals (called no-line-of sight (NLOS) signals) could only be received by reflections off the building [53], and this caused the outline (the huge leaps in 500–600 and the 2000–2500 epochs) in the incremental Chi-square. The delayed effect and the calculation model of the NLOS are discussed in [54], and here we used the incremental Chi-square to detect them, mitigating the influence of this kind of gross error via the robust Kalman filtering in this section.
In PRS-IGG-RKF, the constants c 0 ,   c 1 were set to 1 and 5, respectively. In CI-RKF, the constants c 0 , c 1 were set to 2 and 3, respectively, and α was set to 0.15. Figure 8 illustrates the time series of positioning errors from the three schemes in the E/N/U direction, and Table 1 shows their position errors at the final epoch and the improved accuracy of the other two schemes relative to the CKF.
In Figure 8, the positioning errors from the three schemes coincided during the first 500 epochs, as no gross errors occurred (Figure 7). When the outliers emerged in the 500–600 and 2000–2500 epochs, the positioning results of the PRS-IGG-RKF fluctuated more significantly at these points, and the positioning results of the CI-RKF were relatively smooth as the numbers of iterations increased. The errors of the PRS-IGG-RKF tended to converge after about 4000 epochs, while those of the CI-RKF tended to converge after about 3500 epochs. Finally, the positioning errors of CI-RKF in the E/N/U direction were lower than those of the PRS-IGG-RKF. Both the PRS-IGG-RKF and CI-RKF were able to improve positioning accuracy, but the final positioning accuracies of CI-RKF were 22.68%, 54.33%, and 72.45% higher than those of the PRS-IGG-RKF in the E/N/U direction. The static experiment using GNSS positioning confirms that the Chi-square increment could detect outliers effectively, and the CI-RKF could resist the influence of gross errors in real time by enlarging the observation covariance matrix, providing a better performance compared with the PRS-IGG-RKF.

4.3. Comparison of GNSS Positioning Schemes in the Kinematic Experiment

The incremental Chi-square and its normalized form of the kinematic data from Experiment 2 are shown in Figure 9. It can be seen that gross errors occurred in the intervals within the 100–150 and the 200–300 epochs.
In the PRS-IGG-RKF, the constants c 0 ,   c 1 were 1.5 and 5, respectively; in the CI-RKF, the constants c 0 ,   c 1 were 1 and 4, respectively, and α was set to 0.15. The positioning errors in the E/N/U direction are shown in Figure 10.
During the abnormal epochs (Figure 9), the filtering solutions of CKF fluctuated noticeably, and there were large deviations between the positioning results and the real values. The maximum error of the CKF was about 15 meters at the 100–150 epoch, and 10 meters in the 200–300 epoch. The accuracy of the subsequent solutions was influenced by these outliers. Although there was continuous GNSS navigation, it can be seen that gross errors still had a significant impact on positioning in the area with poor observation. However, the overall performance of the PRS-IGG-RKF and CI-RKF were relatively stable, especially at these epochs, which greatly reduced the errors and caused them to fluctuate near zero in the east and north direction. Table 2 shows the positioning errors of the total region and occluded areas from three schemes. For the total region, the performances of both the PRS-IGG-RKF and CI-RKF were better than of the CKF, and the positioning accuracies of the CI-RKF were about 12.3%, 7.5%, and 16.15% is higher in the E/N/U direction than those of the PRS-IGG-RKF. For the occlusion area, the positioning accuracies of the CI-RKF were 12.11% and 8.99% higher than those of the PRS-IGG-RKF in the east and north directions, but they underperformed in the up direction. Overall, the CI-RKF performed better, but more attention should be paid to the accuracy of the horizontal direction.

4.4. Integrated Navigation Experiment in an Occluded Urban Area

In heavily occluded areas, crowded buildings limit the sky visibility and increase interference from multipath and NLOS signals. The GNSS/IMU/VO integrated navigation system is commonly deployed to improve accuracy and stability of the positioning, while filling in GNSS observation gaps caused by eclipsed signals and huge outliers (in particular for low-cost GNSS receivers). The critical issue then is how to detect and identify outliers effectively when the true trajectory remains unknown. In this experiment, the robust Kalman filtering was an effective method, and three schemes were applied to a GNSS/IMU/VO integrated navigation system. Figure 11 compares the trajectories from three schemes. Fusion combined the GNSS raw observations with the positioning results from the VO/IMU data in a loosely coupled way via the CKF. The robust fusion trajectories of the integrated navigation system were obtained by the PRS-IGG-RKF (Robust Fusion 1) or CI-RKF (Robust Fusion 2). In this part, the robust weighting factors of the CI-RKF were adjusted slightly, according to the IMU’s ability to provide positioning information without GNSS signals over a short time. The GPS information from the mobile phone will be disused, as it provided very poor observations (that is, the corresponding observation noise variance was amplified to infinity).
The results from different Kalman filtering schemes in occluded areas were also projected into Google Earth, as shown in Figure 12. The respective position errors in the E/N/U directions of shelter zones 1 and 2 were calculated (Table 3).
Figure 11 indicates that the results of mobile phone positioning by GPS only had a large deviation, even interruption. In the GNSS/IMU/VO integrated navigation, when the receiver had no sufficient GPS satellite signals, or the observed GPS signals contained severe gross errors (and hence needed to be eliminated), the IMU/VO was able to assist navigation within a short interval, so that users could obtain a sustained positioning trajectory with good stability. The upper left and lower right corner of Figure 11a show the areas where the GPS-only solutions deviated seriously, and the integrated navigation trajectory of GNSS/IMU/VO also fluctuated significantly. This shows that in areas with serious occlusion or where GPS signals experienced gross errors, the CKF integrated navigation was not completely immune and the positioning results still suffered. Figure 11b shows the PRS-IGG-RKF results for data fusion, and that it improved positioning to a certain extent and mitigates *but not completely eliminated) the influence of gross errors. Figure 11c demonstrates the ability of the CI-RKF results to fuse with GNSS/IMU/VO integrated navigation data. This robust fusion method generated a smoother trajectory than in Figure 11a,b, and was closer to the ground truth. Figure 12 zooms the solution trajectories in on the two sheltered areas, displaying that the trajectory of the CI-RKF approach was significantly better than the other two schemes. Table 3 shows the statistics of the position errors for the two occluded areas, and that accuracy was greatly improved when the CI-RKF was applied.

5. Conclusions

This paper derives an incremental Chi-square using Kalman filtering by extending the definition of RLS and comparing its expression with currently adopted Chi-square variables. The comparison indicates that the definitions of Δ χ k 2 from the extension of RLS and ICST are equivalent, and reveals that Δ χ k 2 actually comes from both the change of fitting solutions and the residuals of new observations. In this paper, two methods were proposed for detecting anomalous observations via an incremental Chi-square. The experiment results show that the second method (called the CI-RKF), which provided an adaptive magnification to each term of the observation covariance, was superior. We applied the CI-RKF to the environment with an obvious NLOS effect, and its final positioning accuracies were 28.72%, 95.94%, and 95.82% higher than the CKF, and 22.68%, 54.33%, and 72.45% higher than the PRS-IGG-RKF, respectively. The static experiment indicates that the CI-RKF could detect the outlines and obtain satisfactory positioning solutions. Subsequently, we applied the CI-RKF to two kinematic experiments. First, the CI-RKF was tested using the GNSS positioning model, and the results show that its accuracies were 28.72%, 95.94%, and 95.82% higher than the CKF, and 22.68%, 54.33%, and 72.45% higher than the PRS-IGG-RKF. The overall positioning accuracies were 41.22%, 54.65%, and 18.61% higher than the CKF, and 12.3%, 7.5%, and 16.15% higher than the PRS-IGG-RKF. In next experiment, we replaced the U-blox by a mobile phone with an internal low-cost GNSS antenna, whose observation quality was relatively poorer than the U-blox and easily influenced by surrounding environment. We use its output as the satellite signal for an integrated navigation combining an IMU and VO. Compared with the positioning results obtained from the CKF fusion method, the positioning accuracy of the CI-RKF fusion method increased by 66.73%, 59.59%, and 59.62% in shelter zone 1, and by 42.04%, 59.04%, and 52.41% in shelter zone 2. When compared with the PRS-IGG-RKF fusion method, the CI-RKF was 25.22%, 30.66%, 30.66% for shelter zone 1, and 4.19%, 6.11%, and 16.13% higher for shelter zone 2. These kinematic experiments illustrate that the CI-RKF is able to improve the positioning accuracy of the GNSS and the integrated navigation system, while also enhancing accuracy and stability for real-time positioning in occluded urban areas.
The robust Kalman filtering approach using Δ χ k 2 is general. This paper demonstrates its usage in GNSS time of arrival (TOA) positioning. It can also be used in other positioning approaches, such as time difference of arrival (TDOA), angle of arrival (AOA), received signal strength (RSS), among others.

Author Contributions

Conceptualization, D.D. and W.C.; methodology, B.L. and T.X.; software, B.L and T.X.; validation, W.C., D.D. and B.L.; formal analysis, D.D. and W.C.; investigation, B.L. and W.C.; resources, B.L. and Y.P.; data curation, Y.P. and Z.W.; writing—original draft preparation, B.L.; writing—review and editing, B.L.; visualization, B.L. and Z.W.; supervision, C.Y.; project administration, W.C. and M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is sponsored by the National Key R&D Program of China (No. 2017YFE0100700), the National Natural Science Foundation of China (Nos. 41771475), and the Fund of Director of Key Laboratory of Geographic Information Science (Ministry of Education), East China Normal University (Grant No. KLGIS2017C01).

Acknowledgments

We are very grateful to the editor and anonymous reviewers for their constructive comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Z.; Li, B.; Shen, Y.; Gao, Y.; Wang, M.M. Site-Specific Unmodeled Error Mitigation for GNSS Positioning in Urban Environments Using a Real-Time Adaptive Weighting Model. Remote Sens. 2018, 10, 7. [Google Scholar] [CrossRef] [Green Version]
  2. Teunissen, P.J.G. Quality Control in Integrated Navigation Systems. IEEE Aerosp. Electron. Syst. Mag. 1990, 5, 35–41. [Google Scholar] [CrossRef]
  3. Geng, Y.; Wang, J. Adaptive Estimation of Multiple Fading Factors in Kalman Filter for Navigation Applications. GPS Solut. 2008, 12, 273–279. [Google Scholar] [CrossRef]
  4. Yang, Y. Robust Kalman Filter for Dynamic Systems. J. Pla Inst. Surv. Mapp. 1997, 14, 79–84. [Google Scholar]
  5. Yang, Y.; Ren, X.; Xu, Y. Main Progress of Adaptively Robust Filter with Applications in Navigation. J. Navig. 2013, 1, 9–15. [Google Scholar]
  6. Wang, Q.; Xu, T.; Xu, G. Application of Combining Method of Outlier Detection and Robust Estimation to GPS Kinematic Relative Positioning. Geomat. Inf. Sci. Wuhan Univ. 2011, 36, 476–480. [Google Scholar]
  7. Wang, H.; Li, H.; Zhang, W.; Zuo, J.; Wang, H. Derivative-Free Huber-Kalman Smoothing Based on Alternating Minimization. Signal Process. 2019, 163, 115–122. [Google Scholar] [CrossRef]
  8. Yang, Y.; Song, L.; Xu, T. Robust Estimator for Correlated Observations Based on Bifactor Equivalent Weights. J. Geod. 2002, 76, 353–358. [Google Scholar] [CrossRef]
  9. Yang, Y.; Cheng, M.; Shum, C.; Tapley, B. Robust Estimation of Systematic Errors of Satellite Laser Range. J. Geod. 1999, 73, 345–349. [Google Scholar] [CrossRef]
  10. Yang, Y. Robust Bayesian Estimation. J. Geod. 1991, 65, 145–150. [Google Scholar]
  11. Wu, H.; Chen, S.; Zhang, Y.; Yang, B. Robust Bearings-only Tracking Algorithm Using Structured Total Least Squares-based Kalman Filter. Automatika 2015, 56, 275–280. [Google Scholar] [CrossRef]
  12. Koch, K.; Yang, Y. Robust Kalman Filter for Rank Deficient Observation Models. J. Geod. 1998, 72, 436–441. [Google Scholar] [CrossRef]
  13. Brumback, B.; Srinath, M. A Fault-Tolerant Multisensor Navigation System Design. IEEE Trans. Aerosp. Electron. Syst. 1987, 6, 738–755. [Google Scholar] [CrossRef]
  14. Ren, D. Failure Detection of Dynamical Systems with the State Chi-square Test. J. Guid. Control. Dyn. 1994, 17, 271–277. [Google Scholar]
  15. Groves, P. Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems; Artech House: Norwood, MA, USA, 2007; pp. 458–459. [Google Scholar]
  16. He, Y.; Bilgic, A. Iterative Least Squares Method for Global Positioning System. ARS 2011, 9, 203–208. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, R.; Zhao, C.; Zhang, M.; Sun, P.; Du, X. The Analysis of Adjacent Epoch Error Related Robust Kalman Filtering Algorithm. J. Geo Inf. Sci. 2015, 17, 1506–1510. [Google Scholar]
  18. Dong, D.; Herring, T.A.; King, R.W. Estimating Regional Deformation from a Combination of Space and Terrestrial Geodetic Data. J. Geod. 1998, 72, 200–214. [Google Scholar] [CrossRef]
  19. Angrisano, A.; Gioia, C.; Gaglione, S.; del Core, G. GNSS Reliability Testing in Signal-Degraded Scenario. Int. J. Navig. Obs. 2013, 2013, 1–12. [Google Scholar] [CrossRef]
  20. Zhang, G.; Wen, W.; Hsu, L. Rectification of GNSS-based Collaborative Positioning Using 3D Building Models in Urban Areas. GPS Solut. 2019, 23, 83. [Google Scholar] [CrossRef]
  21. Mohamed, A.; Schwarz, K. Adaptive Kalman Filtering for INS/GPS. J. Geod. 1999, 73, 193–203. [Google Scholar] [CrossRef]
  22. Jo, B.; Baloch, Z. Internet of Things-based Arduino Intelligent Monitoring and Cluster Analysis of Seasonal Variation in Physicochemical Parameters of Jungnangcheon, an Urban Stream. Water 2017, 9, 220. [Google Scholar] [CrossRef] [Green Version]
  23. Hatch, R.; Sharpe, R.; Yang, Y. An Innovative Algorithm for Carrier-Phase Navigation. In Proceedings of the International Technical Meeting of the Satellite Division of the Institute of Navigation, Long Beach, CA, USA, 21–24 September 2004; pp. 1431–1437. [Google Scholar]
  24. Dag, T.; Arsan, T. Received Signal Strength Based Least Squares Lateration Algorithm for Indoor Localization. Comput. Electr. Eng. 2018, 66, 114–126. [Google Scholar] [CrossRef]
  25. Brown, R.; Durbin, J.; Evans, J. Techniques for Testing the Constancy of Regression Relationships over Time. J. R. Stat. Soc. 1975, B37, 149–192. [Google Scholar] [CrossRef]
  26. Belsley, D.; Kuh, E.; Welsch, R. Regression Diagnostics; John Wiley: New York, NY, USA, 1980. [Google Scholar]
  27. Galpin, J.; Hawkins, D. The Use of Recursive Residuals in Checking Model Fit in Linear Regression. Am. Stat. 1984, 38, 94–105. [Google Scholar]
  28. Lawrence, R.; Paul, N. The Kalman Filter: An Introduction to the Mathematics of Linear Least Mean Square Recursive Estimation. Int. J. Math. Educ. Sci. Technol. 1986, 17, 347–366. [Google Scholar]
  29. Kluga, A.; Kluga, J. Dynamic Data Processing with Kalman Filter. Elektron. Elektrotech. 2011, 5, 33–36. [Google Scholar]
  30. Einicke, G.; White, L. Robust Extended Kalman Filtering. IEEE Signal Process. 1999, 47, 2596–2599. [Google Scholar] [CrossRef]
  31. Dong, D.; Chen, J.; Wang, J. Principals of GNSS High Precision Positioning; Science Press: Beijing, China, 2018; pp. 243–248. [Google Scholar]
  32. Xie, Y.; Wang, J.; Liang, Y.; Yu, R. Robust Kalman Filter as a Chemometric Method for Analytical Data Processing. Anal. Chim. Acta 1992, 269, 307–316. [Google Scholar] [CrossRef]
  33. Yang, Y.; Gao, W. An Optimal Adaptive Kalman Filter. J. Geod. 2006, 80, 177–183. [Google Scholar] [CrossRef]
  34. Guo, F.; Zhang, X. Adaptive Robust Kalman Filtering for Precise Point Positioning. Meas. Sci. Technol. 2014, 25, 105011. [Google Scholar] [CrossRef]
  35. Yang, Y.; Gao, W. A New Learning Statistic for Adaptive Filter Based on Predicted Residuals. Prog. Nat. Sci. 2006, 16, 833–837. [Google Scholar]
  36. Mehra, R.; Peschon, J. An Innovations Approach to Fault Detection and Diagnosis in Dynamic Systems. Automatica 1971, 7, 637–640. [Google Scholar] [CrossRef]
  37. Ding, W.; Wang, J.; Rizos, C.; Kinlyside, D. Improving Adaptive Kalman Estimation in GPS/INS Integration. J. Navig. 2007, 60, 517–529. [Google Scholar] [CrossRef] [Green Version]
  38. Yang, Y.; He, H.; Xu, G. Adaptively Robust Filtering for Kinematic Geodetic Positioning. J. Geod. 2001, 75, 109–116. [Google Scholar] [CrossRef]
  39. Huang, J.; Liu, P.; Lu, M.; Xu, B.; Rao, Z.; Wan, Y.; Lu, C.; Zhang, H. Minimum Mean Square Error Estimator for Mobile Location Using the Pseudo-Range TOA Measurements; ICISMME: Chongqing, China, 2015; p. 6. [Google Scholar]
  40. RTKLIB: An Open Source Program Package for GNSS Positioning. Available online: http://www.rtklib.com/ (accessed on 23 January 2017).
  41. Cheng, C.; Tourneret, J.; Quan, P.; Vincent, C. Detecting, Estimating and Correcting Multipath Biases Affecting GNSS Signals Using a Marginalized Likelihood Ratio-based Method. Signal Process. 2015, 118, 221–234. [Google Scholar] [CrossRef] [Green Version]
  42. Chun, T.; Guo, N.; Backen, S.; Akos, D. Monocular Camera/IMU/GNSS Integration for Ground Vehicle Navigation in Challenging GNSS Environments. Sensors 2012, 12, 3162–3185. [Google Scholar]
  43. Caron, F.; Duflos, E.; Pomorski, D.; Vanheeghe, P. GPS/IMU Data Fusion Using Multisensor Kalman Filtering: Introduction of Contextual Aspects. Inf. Fusion 2006, 7, 221–230. [Google Scholar] [CrossRef]
  44. Davison, A.; Reid, I.; Molton, N.; Stasse, O. MonoSLAM: Real-Time Single Camera SLAM. IEEE Pattern Anal. 2007, 29, 1052–1067. [Google Scholar] [CrossRef] [Green Version]
  45. David, N.; Naroditsky, O.; Bergen, J. Visual Odometry for Ground Vehicle Applications. J. Field Robot. 2006, 23, 3–20. [Google Scholar]
  46. Durrant-Whyte, H.; Bailey, T. Simultaneous Location and Mapping: Part I. IEEE Robot. Autom. Mag. 2006, 13, 99–110. [Google Scholar] [CrossRef] [Green Version]
  47. Luan, Y.; Fan, L. Research on GPS/INS Integrated Navigation system. Microcomput. Inf. 2009, 34, 203–204. [Google Scholar]
  48. Liu, Y.; Chen, Z.; Zheng, W.; Wang, H.; Liu, J. Monocular Visual-Inertial Slam: Continuous Preintegration and Reliable Initialization. Sensors 2017, 17, 2613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Kim, J.; Lyou, J.; Kwak, H. Vision Coupled GPS/INS Scheme for Helicopter Navigation. J. Mech. Sci. Technol. 2010, 24, 489–496. [Google Scholar] [CrossRef]
  50. Liu, F. Tightly Coupled Integration of GNSS/INS/Stereo Vision/Map Matching System for Land Vehicle Navigation; UCGE Report; University of Calgary: Calgary, AB, Canada, 2018. [Google Scholar]
  51. Liu, Y.; Fan, X.; Chen, L.; WU, J.; Li, L.; Ding, D. An Innovative Information Fusion Method with Adaptive Kalman Filter for Integrated INS/GPS Navigation of Autonomous Vehicles. Mech. Syst. Signal Process. 2018, 100, 605–616. [Google Scholar] [CrossRef] [Green Version]
  52. NSL. Launches a New Free Android App as Part of FLAMINGO—Discover rinexON. Available online: https://www.flamingognss.com/rinexon (accessed on 10 September 2018).
  53. Groves, P. Shadow Matching: A new GNSS Positioning Technique for Urban Canyons. J. Navig. 2011, 64, 417–430. [Google Scholar] [CrossRef]
  54. Betaille, D.; Peyret, F.; Ortiz, M.; Miquel, S.; Fontenay, L. A new Modeling Based on Urban trenched to Improve GNSS positioning Quality of Service in Cities. IEEE Intell. Transp. Syst. Mag. 2013, 5, 59–70. [Google Scholar] [CrossRef]
Figure 1. The calculation process of the incremental Chi-square (CI-RKF) using Global Navigation Satellite System (GNSS) positioning.
Figure 1. The calculation process of the incremental Chi-square (CI-RKF) using Global Navigation Satellite System (GNSS) positioning.
Remotesensing 12 00732 g001
Figure 2. GNSS/IMU/VO coupling model; P represents the position, and V represents the velocity.
Figure 2. GNSS/IMU/VO coupling model; P represents the position, and V represents the velocity.
Remotesensing 12 00732 g002
Figure 3. The setup of the static experiment using GNSS positioning; (a) overview of the test platform (marked by the line of dashes) in Google Earth; (b) observation environment of the static experiment.
Figure 3. The setup of the static experiment using GNSS positioning; (a) overview of the test platform (marked by the line of dashes) in Google Earth; (b) observation environment of the static experiment.
Remotesensing 12 00732 g003
Figure 4. The setup of the kinematic experiment using GNSS positioning; (a) trajectory of the vehicle; (b) onboard equipment for the kinematic experiment.
Figure 4. The setup of the kinematic experiment using GNSS positioning; (a) trajectory of the vehicle; (b) onboard equipment for the kinematic experiment.
Remotesensing 12 00732 g004
Figure 5. The setup of the kinematic experiment using integrated positioning; (a) integrated navigation vehicle experiment; (b) Huawei Honor V9 placed in the carriage to receive satellite signals.
Figure 5. The setup of the kinematic experiment using integrated positioning; (a) integrated navigation vehicle experiment; (b) Huawei Honor V9 placed in the carriage to receive satellite signals.
Remotesensing 12 00732 g005
Figure 6. Comparison of robust incremental Chi-square methods.
Figure 6. Comparison of robust incremental Chi-square methods.
Remotesensing 12 00732 g006
Figure 7. Chi-square increment in the static experiment: (a) Chi-square increment; (b) normalized Chi-square increment.
Figure 7. Chi-square increment in the static experiment: (a) Chi-square increment; (b) normalized Chi-square increment.
Remotesensing 12 00732 g007
Figure 8. Position errors in the E/N/U direction in the static experiment.
Figure 8. Position errors in the E/N/U direction in the static experiment.
Remotesensing 12 00732 g008
Figure 9. Chi-square increment in the GNSS kinematic experiment: (a) Chi-square increment; (b) normalized Chi-square increment.
Figure 9. Chi-square increment in the GNSS kinematic experiment: (a) Chi-square increment; (b) normalized Chi-square increment.
Remotesensing 12 00732 g009
Figure 10. Position errors in the E/N/U direction in the kinematic experiment.
Figure 10. Position errors in the E/N/U direction in the kinematic experiment.
Remotesensing 12 00732 g010
Figure 11. Comparison of the trajectories from three schemes: (a) the positioning results of the CKF fusion method, with the black trajectory as the true value, the green trajectory as the GPS positioning results, and the red as the result of the fusion method; (b) the positioning results of the PRS-IGG-RKF fusion method, with the same color representations as (a); (c) the positioning results of the CI-RKF fusion method, with the same color representations as (a).
Figure 11. Comparison of the trajectories from three schemes: (a) the positioning results of the CKF fusion method, with the black trajectory as the true value, the green trajectory as the GPS positioning results, and the red as the result of the fusion method; (b) the positioning results of the PRS-IGG-RKF fusion method, with the same color representations as (a); (c) the positioning results of the CI-RKF fusion method, with the same color representations as (a).
Remotesensing 12 00732 g011
Figure 12. Performance of different Kalman filtering positioning results in shelter zones: (a) shelter zone 1, with the true trajectory represented by yellow points, the positioning results of the CKF fusion method in the integrated navigation system represented by green points, the positioning results of the CI-RKF fusion method in the integrated navigation system represented by red points, and the positioning results of the PRS-IGG-RKF fusion method in the integrated navigation system represented by orange points; the numbers are the points’ index; and (b) shelter zone 2, with the same color and number representations as (a).
Figure 12. Performance of different Kalman filtering positioning results in shelter zones: (a) shelter zone 1, with the true trajectory represented by yellow points, the positioning results of the CKF fusion method in the integrated navigation system represented by green points, the positioning results of the CI-RKF fusion method in the integrated navigation system represented by red points, and the positioning results of the PRS-IGG-RKF fusion method in the integrated navigation system represented by orange points; the numbers are the points’ index; and (b) shelter zone 2, with the same color and number representations as (a).
Remotesensing 12 00732 g012
Table 1. Position errors at the final epoch in the E/N/U direction in the static experiment.
Table 1. Position errors at the final epoch in the E/N/U direction in the static experiment.
DirectionPosition Error (m)Improve (%)
CKFPRS-IGG-RKFCI-RKFPRS-IGG-RKFCI-RKF
East−3.012−2.830−2.1476.0428.72
North1.3050.762−0.05341.6195.94
Up3.2262.4720.13523.3795.82
Table 2. The RMS of position errors in the E/N/U direction in the kinematic experiment.
Table 2. The RMS of position errors in the E/N/U direction in the kinematic experiment.
Position Error (m)Improve (%)
CKFPRS-IGG-RKFCI-RKFPRS-IGG-RKFCI-RKF
OverallEast2.3341.6591.37228.9241.22
North2.8601.5111.29747.1554.65
Up4.5474.4353.7012.4618.61
Shelter ZonesEast3.5852.3901.95633.3345.44
North3.8882.0241.67547.9456.93
Up6.4283.3113.68748.5042.64
Table 3. The RMS of position errors in the E/N/U direction for shelter zones in the integrated navigation experiment.
Table 3. The RMS of position errors in the E/N/U direction for shelter zones in the integrated navigation experiment.
Position Error (m)Improve (%)
FusionRobust Fusion 1Robust Fusion 2Robust Fusion 1Robust Fusion 2
Shelter Zone 1East4.9392.8891.64341.5166.73
North2.0291.4420.82028.9359.59
Up2.3651.6800.95528.9659.62
Shelter Zone 2East25.78216.03014.94437.8542.04
North12.9466.0945.30352.9359.04
Up15.0899.6157.18136.2852.41

Share and Cite

MDPI and ACS Style

Li, B.; Chen, W.; Peng, Y.; Dong, D.; Wang, Z.; Xiao, T.; Yu, C.; Liu, M. Robust Kalman Filtering Based on Chi-square Increment and Its Application. Remote Sens. 2020, 12, 732. https://doi.org/10.3390/rs12040732

AMA Style

Li B, Chen W, Peng Y, Dong D, Wang Z, Xiao T, Yu C, Liu M. Robust Kalman Filtering Based on Chi-square Increment and Its Application. Remote Sensing. 2020; 12(4):732. https://doi.org/10.3390/rs12040732

Chicago/Turabian Style

Li, Bo, Wen Chen, Yu Peng, Danan Dong, Zhiren Wang, Tingting Xiao, Chao Yu, and Min Liu. 2020. "Robust Kalman Filtering Based on Chi-square Increment and Its Application" Remote Sensing 12, no. 4: 732. https://doi.org/10.3390/rs12040732

APA Style

Li, B., Chen, W., Peng, Y., Dong, D., Wang, Z., Xiao, T., Yu, C., & Liu, M. (2020). Robust Kalman Filtering Based on Chi-square Increment and Its Application. Remote Sensing, 12(4), 732. https://doi.org/10.3390/rs12040732

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