Next Article in Journal
An Autonomous Cooperative Navigation Approach for Multiple Unmanned Ground Vehicles in a Variable Communication Environment
Previous Article in Journal
A Comprehensive Literature Review on Volatile Memory Forensics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dual Kalman Filter Based on a Single Direction under Colored Measurement Noise for INS-Based Integrated Human Localization

1
School of Civil Engineering, Qingdao University of Technology, Qingdao 266520, China
2
Shandong Luqiao Group Company, Ltd., Jinan 250014, China
3
School of Electrical Engineering, University of Jinan, Jinan 250022, China
*
Author to whom correspondence should be addressed.
Electronics 2024, 13(15), 3027; https://doi.org/10.3390/electronics13153027
Submission received: 21 May 2024 / Revised: 16 July 2024 / Accepted: 26 July 2024 / Published: 31 July 2024

Abstract

:
For inertial-based integrated pedestrian navigation, the navigation environment might affect the positioning accuracy in different directions. Meanwhile, complex filtering algorithms can reduce computational efficiency. Therefore, one dual Kalman filter (KF) based on a single direction under a colored measurement noise (CMN) scheme is developed herein to improve the robustness and operational efficiency. The proposed method involves designing a data fusion model for the KF that integrates data from an inertial navigation system (INS) and ultrawideband (UWB). Subsequently, the INS/UWB integrated model-based KF under CMN (cKF) will be derived. Then, two sub-cKFs are proposed to fuse the data in the east and north directions, respectively. The empirical findings highlight the superior performance of the proposed approach over the KF for position estimation accuracy and runtime reduction, demonstrating its effectiveness.

1. Introduction

Pedestrian positioning has extensive applications across diverse industries. In intelligent transportation, precise positioning can help urban planners better understand the flow of people, optimizing road design and traffic management [1]. In security monitoring, pedestrian positioning can help monitors respond to security incidents quicker and more accurately [2]. In healthcare, pedestrian positioning can help monitor patient activity, provide better care and management, and more [3,4]. However, with increasing diversity of application areas and the complexity of positioning environments, higher demands are being placed on pedestrian positioning technology. Higher positioning accuracy, faster response time, lower power consumption, and better ability to adapt to various environments are required [5,6]. Therefore, the development and innovation of pedestrian positioning technologies have become critical to meet the growing demands and promote the development and progress of various industries.
Many approaches have been proposed to achieve human localization. Among them, inertial navigation systems (INSs) have long been a cornerstone of navigation and localization, offering continuous position updates by integrating measurements from inertial sensors. However, the performance of an INS degrades over time due to sensor errors, drift, and noise accumulation [7,8]. Therefore, researchers have proposed integrating INSs with other sensors, such as global positioning system (GPS) receivers, to improve localization accuracy [9,10]. This paper investigates data fusion with an INS using ultrawideband (UWB). With technological advancements and large-scale production in recent years, the implementation cost of UWB technology has been gradually reduced. The low power consumption, high-precision positioning, and strong penetration ability are more prominent than in other positioning sensors. Therefore, this sensor is selected for data fusion [11,12,13,14]. This combination approach has critical applications in many fields [15,16]. For example, Pan and Jia et al. proposed a tight combination of a UWB single-base station and an INS for smart vehicle localization through the data fusion of the UWB single-base station distance and angle measurements for the narrow and long space environment [17]. In indoor localization, Kang and Duan et al. used UWB normal range measurements to correct the predicted position of the employed INS for improving the problem of non-line-of-sight interference affecting localization accuracy [18]. In hospitals, wheelchair or pedestrian high accuracy positioning and continuous detection, the user’s safety, and ensuring normal workflow are crucial. In hospital-oriented outdoor scenarios, complex environments, such as centralized high-rise buildings and trees, frequently lead to satellite signal loss, resulting in discontinuous positioning information. Therefore, the INS system is used for positioning in the case of satellite signal loss to ensure the continuous detection of pedestrians in hospitals [19,20,21,22]. Note that the Inertial Measurement Units (IMUs) typically operate at significantly higher sampling frequencies compared to GPS and/or UWB systems. This discrepancy makes the study of localization using UWB and IMU particularly relevant to the academic community [23].
The diversity of localization techniques also means that the use of data fusion filters is more demanding to improve localization accuracy. The Kalman filter (KF) is a widely used sensor fusion technique in integrated navigation systems. It combines noisy sensor measurements with a dynamic model of the system to estimate its true state [24,25,26]. Researchers have proposed many types of KFs for different models and application scenarios. To reduce the computational load of nonlinear navigation systems, Zhao introduced a hybrid method combining the cubature KF (CKF) and the extended KF (EKF). This approach employs a linear filtering method for estimating the linear state and a nonlinear filtering method for estimating the nonlinear state [27]. The hybrid filtering method retains the CKF’s ability to manage nonlinear problems while reducing the system’s computation time. In bearing-only tracking (BOT) systems, Sun and Davies et al. proposed the adaptive kernel KF (AKKF), which substantially reduces the number of particles and improves the performance of the algorithm compared to the unseasoned Kalman filter (UKF), particle filter (PF), and Gaussian particle filter (GPF) [28]. The creation of hybrid KF and neural network models is a considerable development in state estimation, and this hybrid model has high flexibility in its use, either alternatively or in a mixture [29]. The hybrid model can be divided into two types: the state–space model, where the KF neural network trains the equations or parameters, and the neural network model, whose parameters are updated using the KF [30,31]. This study shows that the hybrid model performs better in accuracy and flexibility than a single KF or neural network model [32,33]. The traditional KF assumes that the measurement noise is Gaussian white noise; however, when applied in real-world environments, it typically exhibits colored characteristics, meaning that the noise at one point correlates with the noise at adjacent points [34,35]. This outcome violates the white noise assumption and can lead to suboptimal fusion performance. To address this limitation, the KF filter is proposed to be improved for the colored measurement noise (CMN) model using the backward Euler (BE) method, and the effectiveness of this method for CMN removal is also demonstrated; however, the system operation speed is not considered [36].
In order to improve the localization accuracy and reduce the computational load, this paper proposes a unidirectional dual Kalman-based filter for the sensor localization problem under CMN by using matrix segmentation. First, the data from the sensor is processed, and the system is modeled. Then, the state vector is divided into two states in a single direction, and one KF is used to estimate the state in this direction. Another KF is used to estimate the state in the other direction while considering the effect of cross-covariance on these two subfilters. Through experimental validation in real-world environments, the results show that the proposed dual KF has better localization accuracy and stability under colored noise than traditional methods, and there is an improvement in the system operation speed.
This study has made notable progress in the following aspects:
  • First, a single-direction equation of state and an observation equation are established based on the characteristics of inertial sensors and UWB measurements.
  • Then, for the influence of the navigation environment on the positioning accuracy in different directions, a dual KF based on a single direction under CMN is proposed to effectively address the positioning issue in the presence of CMN.
  • Finally, the feasibility and effectiveness of the proposed method in practical applications are proved through actual measurements, and a comparison is made regarding the operation speed of the system. The experimental results show that the algorithm proposed herein is valuable for improving the system running time.
The remaining sections of this paper are organized as follows: first, a model for the unidirectional INS/UWB-based navigation system is introduced. Then, the unidirectional dual KF-based method under CMN is described in detail. Next, experimental tests are conducted to compare the positioning accuracy and system runtime for KF, dual KF, and dual cKF. Finally, conclusions are given.

2. Dual KF based on a Single Direction under CMN

This section outlines the derivation of the dual KF under CMN, focusing on a single direction. The derivation is based on the model presented in Section 2. First, we will derive the improved data fusion model under CMN and its KF. Then, we split the filter into two subfilters according to different directions.

2.1. KF under CMN

In this article, we consider a system composed of Equations (1) and (2).
s n = A s n 1 + μ n ,
z n = H n s n + a η n 1 + v n ,
where s n is the state vector, n is the time index, A is system matrix, μ n N ( 0 , Q ) is the system noise, and Q is its covariance, a is the color measurement noise factor, η n 1 is the color measurement noise, and v n N ( 0 , R ) is the measurement noise, and R is its covariance.
We reduce the effect of the CMN by introducing one new measurement y n , which can be computed using Equation (2).
y n = z n a z n 1 ,
With Equation (1), we build one new measurement equation, which is shown in Equation (2).
y n = H n U n s n + U n μ n + γ n = D n s n + γ ¯ n ,
where U n = a H n A , D n = H n U n is the observation matrix, and  γ ¯ n = U n μ n + γ n is the measurement noise.
R ¯ n = E γ ¯ n γ ¯ n T = U n Q U n T + R n = U n Φ n + R n ,
where Φ n = Q U n T . Now, with the model (1) (3), we can list the method of the KF under CMN as Algorithm 1. With the new measurement Equation (3), we compute the new equation of state.
Algorithm 1: Kalman filter (KF) under CMN based on the model (27) (4) when μ n and γ ¯ n are correlated
Electronics 13 03027 i001
Equation (1) can be rewritten as (6).
s n = A s n 1 + μ n + β n [ y n ( D n s n + γ ¯ n ) ] = ( I β n D n ) A F n s n 1 + β n y n w n + ( I β n D n ) μ n β n γ ¯ n v n = F n s n 1 + w n + v n ,
where w n is the input, and  v n N ( 0 , Θ n ) is the new system noise, which can be computed using Equation (7). The algorithm for the KF under CMN is presented as Algorithm 2.
Θ n = E [ ( I β n D n ) μ n β n γ ¯ n ] [ ( I β n D n ) μ n β n γ ¯ n ] T = ( I β n H n ) Q ( I β n H n ) T + β n R β n T ,
Algorithm 2: Kalman filter (KF) under CMN based on the model (27) (4) when μ n and γ ¯ n are not correlated
Electronics 13 03027 i002

2.2. Dual KF under CMN Based on a Single Direction

This study employs the dual KF. For the direction, we rewrite the state vector as s n = Δ x n Δ v x n | Δ y n Δ v y n T = s x n | s y n T , and we rewrite the following matrices: U n = U x n U y n . Here, U x n = ( a H x n ) / A x n , U y n = ( a H y n ) / A y n , D n = D x n D y n . Here, D x n = H x n U x n , D y n = H y n U y n , and Φ n = Φ x n Φ y n T . Here, Φ x n = Q x n U x n T , Φ y n = Q y n U y n T ,
H n = x n I x n 1 x n I x n 1 2 + y n I y n 1 2 1 2 0 y n I y n i x n I x n 1 2 + y n I y n 1 2 1 2 0 x n I x n 2 x n I x n 2 2 + y n I y n 2 2 1 2 0 y n I y n i x n I x n 2 2 + y n I y n 2 2 1 2 0 x n I x n 3 x n I x n 3 2 + y n I y n 3 2 1 2 0 y n I y n i x n I x n 3 2 + y n I y n 3 2 1 2 0 x n I x n 4 x n I x n 4 2 + y n I y n 4 2 1 2 0 y n I y n i x n I x n i 2 + y n I y n i 2 1 2 0 = H x n | H y n ,
Thus, we can rewrite Equation (32) so that
z n = H x n s x n + H y n s y n + a η n 1 + v n ,
and R ¯ n can be computed as
R ¯ n = U x n Φ x n + U y n Φ y n + R n ,
When v n and γ ¯ n is uncorrelated, we can obtain
β n = Φ n ( H n Φ n + R ) 1 = Φ x n ( H x n Φ x n + R ) 1 0 0 Φ x n ( H x n Φ x n + R ) 1 = β x n 0 0 β y n ,
Θ n = ( I β x n H x n ) Q x ( I β x n H x n ) T 0 0 ( I β y n H y n ) Q y ( I β y n H y n ) T = Θ x n 0 0 Θ y n ,
where β x n = Φ x n ( H x n Φ x n + R ) 1 , β y n = Φ y n ( H y n Φ y n + R ) 1 ,
F n = ( I β x n D x n ) A x 0 0 ( I β y n D y n ) A y = F x n 0 0 F y n ,
In this work, we set K ¯ n = K ¯ x n K ¯ y n , P ^ n = P ^ 11 n P ^ 12 n P ^ 21 n P ^ 22 n , thus, we can compute the following equations.
P ^ 11 n = F x n P ^ 11 n F x n T + Θ x n ,
P ^ 12 n = F x n P ^ 12 n F y n T ,
P ^ 21 n = F y n P ^ 21 n F x n T ,
P ^ 22 n = F y n P ^ 22 n F y n T + Θ y n ,
Based on the equations above, we can obtain the measurement update steps based on the direction, which are listed in Equations (18) and (19).
s x n = s x n + K ¯ x n ( y n D x n x n D y n y n ) ,
s y n = s y n + K ¯ y n ( y n D x n x n D y n y n ) ,
The KF’s gain under CMN can be computed using the formulas given in Equations (20) and (21).
K ¯ x n = P ^ 11 n D x n T + P ^ 12 n D y n T G ¯ n 1 ,
K ¯ y n = P ^ 21 n D x n T + P ^ 22 n D y n T G ¯ n 1 ,
 where
G n = D x n P ^ 11 n D x n T + D y n P ^ 21 n D x n T + D x n P ^ 12 n D y n T + D y n P ^ 22 n D y n T + U x n Φ x n + U y n Φ y n + R ,
We can obtain the following equations:
P ^ 11 n = P ^ 11 n ( K ¯ x n D x n P ^ 11 n + K ¯ x n D y n P ^ 21 n ) ,
P ^ 12 n = P ^ 12 n ( K ¯ x n D x n P ^ 12 n + K ¯ x n D y n P ^ 22 n ) ,
P ^ 21 n = P ^ 21 n ( K ¯ y n D x n P ^ 11 n + K ¯ y n D y n P ^ 21 n ) ,
P ^ 22 n = P ^ 22 n ( K ¯ y n D x n P ^ 12 n + K ¯ y n D y n P ^ 22 n ) ,
   The algorithm for the KF under CMN based on a single direction is presented as Algorithm 3. Figure 1 illustrates the configuration of the KF under CMN based on a single direction. The figure shows that in the proposed method, the state vector is split into two subfilters based on direction. Each subfilter estimates its own one-step prediction and measurement update with the distance difference between the INS and UWB. Additionally, the one-step prediction from one subfilter is used to assist with the measurement update of the other. Finally, the outputs of both subfilters are merged to calculate the final output.
Algorithm 3: Kalman filter (KF) under the colored measurement noise (CMN) based on a single direction.
Electronics 13 03027 i003

3. UWB-Assisted INS-Based Human Localization Using a Dual KF Based on a Single Direction under CMN

This section focuses on using dual KFs for UWB-assisted INS human positioning under CMN. First, the scheme for the UWB-assisted INS-based human localization will be designed. Second, a data fusion model will be derived based on the characteristics of the measured data.

3.1. Scheme of the Dual KF under CMN

Herein, we employ UWB to assist with INS human localization. Figure 2 illustrates the configuration of the dual KF designed for INS-based integrated human localization under CMN based on a single direction. The foot-mounted inertial measurement unit (IMU) measures the human position, and the INS-based distances r n I , 1 , 2 , 3 , 4 between the human and the UWB reference nodes (RNs) are computed. Furthermore, the UWB-based distances r n U , 1 , 2 , 3 , 4 between the human and UWB RNs are computed. It should be pointed out that the IMUs cannot directly measure position, it uses acceleration and gyroscope data to calculate the velocity in the navigation frame and further calculate the position. Then, the difference of Δ r n I , 1 , 2 , 3 , 4 = r n I , 1 , 2 , 3 , 4 r n U , 1 , 2 , 3 , 4 is computed, serving as an input to the dual KF under CMN. The filter produces an error estimate of the INS position L n I . Then, the difference of L n I Δ L ^ n I is computed as the final output.

3.2. The Data Fusion Model for the Dual KF under CMN

This subsection focuses on deriving the data fusion model for the proposed dual KF under CMN. Equation (27) presents the state equation utilized in this study.
Δ x n Δ v x n Δ y n Δ v y n s n = A x 0 0 A y A Δ x n 1 Δ v x n 1 Δ y n 1 Δ v y n 1 s n 1 + μ n ,
A x = A y = 1 Δ n 0 1 ,
where ( Δ x n , Δ y n ) is the position error at the time index n, ( Δ v x n , Δ v y n ) is the velocity error at the time index n, Δ n is the sample time, μ n N ( 0 , Q n ) is the system noise, and Q n = Q Δ x , n 0 0 Q Δ y , n is its covariance.
For the measurement equation, the distance r n I , i , i [ 1 , 4 ] measured by the INS can be computed using Equation (29).
r n I , i = x n I x n i 2 + y n I y n i 2 1 2 , i 1 , 4 ,
where ( x n I , y n I ) is the position measured by the INS at the time index n, and ( x n i , y n i ) , i [ 1 , 4 ] is the i t h UWB RN’s position. Perform a first-order Taylor expansion on the above equation to obtain
r n I , i = r n R , i + x n I x n i x n I x n i 2 + y n I y n i 2 1 2 Δ x n + y n I y n i x n I x n i 2 + y n I y n i 2 1 2 Δ y n , i 1 , 4 ,
where r n R , i is the real position at the time index n. Then, we can obtain the following equation:
r n I , i r n U , i = x n I x n i x n I x n i 2 + y n I y n i 2 1 2 Δ x n + y n I y n i x n I x n i 2 + y n I y n i 2 1 2 Δ y n + η n i , i 1 , 4 ,
where η n i represents the CMN for the i t h distance at time index n. It can be expressed by η n i = a i η n 1 i + v n i . Here, v n i N ( 0 , R ) is the Gaussian noise. Thus, we obtained the measurement equation employed in this study.
r n I , 1 r n U , 1 r n I , 2 r n U , 2 r n I , 3 r n U , 3 r n I , 4 r n U , 4 z n = x n I x n 1 x n I x n 1 2 + y n I y n 1 2 1 2 0 y n I y n i x n I x n 1 2 + y n I y n 1 2 1 2 0 x n I x n 2 x n I x n 2 2 + y n I y n 2 2 1 2 0 y n I y n i x n I x n 2 2 + y n I y n 2 2 1 2 0 x n I x n 3 x n I x n 3 2 + y n I y n 3 2 1 2 0 y n I y n i x n I x n 3 2 + y n I y n 3 2 1 2 0 x n I x n 4 x n I x n 4 2 + y n I y n 4 2 1 2 0 y n I y n i x n I x n i 2 + y n I y n i 2 1 2 0 H n Δ x n Δ v x n Δ y n Δ v y n s n + a η n 1 1 η n 1 2 η n 1 3 η n 1 4 η n 1 + v n 1 v n 2 v n 3 v n 4 v n , ,

4. Discussion

In this section, a self-built experiment is conducted to analyze the localization accuracy and algorithm runtime of the proposed algorithm. First, the initial stage of the experiment will be introduced. Then, the proposed dual KF based on a single direction under CMN will be investigated for localization performance.

4.1. Initial Experimental Stage

In this study, the real-world test is conducted near the University of Jinan’s No. 1 building, as depicted in Figure 3. The figure shows that we employ four UWB RNs placed around the experimental site. The target pedestrian is equipped with a foot-mounted IMU, real-time kinematic (RTK), UWB blind node (BN), and one computer. In order to obtain the 2D human’s position, the layout of UWB base stations needs to ensure that the target pedestrian can obtain the distance of at least three UWB stations while walking; the parameter of the UWB localization system used in this work is listed in Table 1. In this article, we ensure that the distance information from all the UWB stations can be obtained during pedestrian walking. Here, the foot-mounted IMU measures the r n I , 1 , 2 , 3 , 4 , the UWB BN measures the r n U , 1 , 2 , 3 , 4 , and the RTK is the reference standard for the experiment; its parameter is listed in Table 2. All sensor data are input to the computer. The test structure is shown in Figure 4. Note that the inertial sensors need to be calibrated before the test, including static and dynamic calibration. Static calibration is performed without any movement by recording the output offset value of the sensor and correcting it to the zero position. Dynamic calibration criteria are used to calibrate sensors during specific actions or movements. By performing actions or movements under known conditions and comparing the sensor output, the offset error of the sensor can be accurately calculated and corrected. In this work, we use the software provided by the sensor to compensate before the experiment begins directly.

4.2. Localization Accuracy Investigation

In this subsection, we performed two real tests using the test system proposed in the previous section. This study employs the KF and dual KF as the comparison filters.
Figure 5 shows the reference path and those measured using the INS + zero velocity update (ZUPT), UWB, KF, dual KF, and the dual cKF in Test 1. In this figure, we employ the black dashed line to denote the reference path, the gray solid line is used to denote the INS with ZUPT method, the green solid line is used to denote the UWB-only solution, the cyan dashed line is employed to denote the path measured by the KF, the blue dashed line represents the trajectory obtained from the dual KF, while the proposed dual KF based on a single direction under CMN (denoted as dual cKF) is depicted by the red dashed line. Meanwhile, we use pink to show the positions of the UWB RNs, and the start and finish points are also shown in this figure.
From this figure, we can easily find that while ZUPT effectively reduces position drift errors, there is a noticeable error accumulation over time. The UWB system demonstrates stability in determining the target human’s position, as observed in its solution when compared to INS. All the filters’ solutions can provide trajectories similar to the reference trajectory, which are better than the UWB’s solution. The solutions of the KF and dual KF methods are similar to those in Test 1. The proposed dual cKF method’s path can provide a stable solution.
Figure 6 and Figure 7 display the reference value and the position in the east and north directions measured by the INS + ZUPT, UWB, KF, dual KF, and the dual cKF in Test 1. In this figure, we employ the black dashed line to denote the reference path, the gray solid line is used to denote the INS with the ZUPT method, the green dashed line is used to denote the UWB-only solution, the cyan solid line is employed to denote the position measured by the KF, the solid blue line represents the trajectory derived from the dual KF method, while the proposed dual cKF method is indicated by the solid red line. From the Figure 6, we can see that the difference between the value of the INS+ZUPT and the reference value is significant. When compared with the INS+ZUPT’s solution, the UWB and three filter solutions are closer to the reference value. Here, the Dual KF’s position is similar to the KF, and its value is slightly closer to the reference value, which shows the effectiveness of the dual filter model in reducing the localization error. Compared with the dual KF, the dual cKF’s curve is smoother. Although the values of the cKF algorithm are far from the reference value compared to the KF algorithm in some time periods, most of the time, its estimated values are closer to the reference value than KF and dual KF algorithms, which shows the effectiveness. In the north direction, from Figure 7, we can see that dual filter methods’ solutions are closer to the reference value, and the proposed method can provide a closer solution than the other methods, demonstrating the efficacy of the proposed approach.
Figure 8 shows the UWB distance errors and the position errors of the methods in Test 1. Here, the distance i , i [ 1 , 4 ] means the distance error between the target human and the ith UWB RN. From the figure, we can infer that the proposed method has the best performance. Figure 9 shows the position root-mean-square errors (RMSEs) in the east and north directions measured by UWB, the KF, the dual KF, and the dual cKF in Test 1. In this figure, we employ the green solid line to denote the UWB solution’s RMSEs, the cyan solid line is employed to denote the position measured by the dual KF, the solution obtained from the KF method is represented by the solid blue line, while the proposed dual cKF method is depicted by the solid red line.
Figure 9a shows that all the filters are effective in reducing the localization error when compared with the UWB’s solution in the east direction. When compared with the KF, the dual KF can further improve the localization error, which shows the dual model’s effectiveness. To the north direction, from Figure 9b, we can see that the KF’s position RMSE in the north direction is smaller than that of the UWB at the last time index. The dual KF provides the north position with a smaller RMSE than the KF method does at the last time index. The proposed dual cKF method’s north position RMSE is the smallest at the last time index, highlighting the effectiveness of the proposed approach.
Figure 10 illustrates the cumulative distribution function (CDF) of position error measured by the UWB, KF, dual KF, and dual cKF in Test 1. In this figure, we employ the black solid line to denote the UWB solution’s position error, the blue solid line is employed to denote the position measured by the dual KF, the solid green line represents the trajectory derived from the KF method, while the proposed dual cKF method is indicated by the solid red line. This figure illustrates that the proposed method achieves the smallest position error compared to other methods at the 0.9 point. This result underscores the effectiveness of the proposed approach in reducing localization errors.
Table 3 shows the RMSEs of the INS, UWB, KF, dual KF, and the dual cKF in Test 1. From the table, we can see that the UWB’s solution has better localization performance, and its mean localization error is 0.185 m. Both the KF and dual KF methods decrease the localization error from 0.185 m to 0.154 m and 0.140 m, respectively. In contrast, the proposed dual cKF achieves a mean localization error of 0.131 m, representing a reduction in localization error of approximately 29.19% compared to the UWB solution.
Table 4 shows the average running time of the KF, dual KF, and dual cKF in Test 1. Note that the sub-KF 1 and 2 are the subfilters of the dual KF, and the sub-cKF 1 and 2 are the subfilters of the dual cKF. The table shows that the average running time of the subfilters of the dual KF and dual cKF is lower than that of the KF method. We must emphasize that in actual operation, the two subfilters of the dual KF and dual cKF run in parallel; therefore, the running time of the dual KF and dual cKF can take the larger of the two subfilters. Therefore, the proposed method effectively reduces the running time.
In Test 2, the reference path and those for the INS + ZUPT, UWB, KF, dual KF, and dual cKF measurements are shown in Figure 11. In this figure, we employ the black dashed line to denote the reference path, the gray solid line is used to denote the INS with the ZUPT method, the green solid line is used to denote the UWB-only solution, the cyan dashed line is employed to denote the path measured by the KF, the trajectory from the dual KF method is shown as the blue dashed line, and the proposed dual cKF method is depicted by the red dashed line. Meanwhile, we use pink to show the positions of the UWB RNs, and the start and finish points are also shown.
Similar to Test 1, the figures show a noticeable accumulation of errors in the INS’s solution. All the filters’ and the UWB’s solutions can provide positions similar to the reference value. Compared with other methods, the proposed dual cKF method can provide a solution closer to the reference path, demonstrating the effectiveness of this approach.
Figure 12 and Figure 13 show the positions of the INS + ZUPT, UWB, KF, dual KF, and dual cKF measured in the east and north directions in Test 2 and the reference values. The figures show different trajectories in Test 2. Figure 14 shows the UWB distance errors and the position errors of the methods in Test 2. From the figure, we can infer that the proposed method has the best performance. The black dashed line indicates the reference path, while the gray solid line shows the INS with the ZUPT method, the green dashed line represents the UWB-only solution, the cyan solid line depicts the position estimated by the KF, the blue solid line illustrates the trajectory from the dual KF method, and the red solid line represents the trajectory from the proposed dual cKF method. From Figure 12, we can see that the INS’s east position has some big errors from the time index of 100 to 250 and from the time index of 500 to 600. Moreover, all the filters’ solutions can provide positions similar to the reference value, which is better than the UWB’s solution. Compared to the other methods, the proposed dual cKF method demonstrates its effectiveness by providing a solution that closely matches the reference path. From Figure 13, we can see that all the filters’ solutions can provide positions similar to the reference value, which is better than the UWB’s solution. Compared to the other methods, the proposed dual cKF method demonstrates its effectiveness by providing a solution that closely matches the reference path.
Figure 15 shows the position RMSEs in the east and north directions for UWB, the KF, the dual KF, and the dual cKF measured in Test 2. The green solid line represents the RMSEs of the UWB solution, the cyan solid line indicates the position measured by the dual KF, the solution derived from the KF method is shown by the blue solid line, and the proposed dual cKF method is depicted by the red solid line. From the figures, we can see that the KF’s position RMSE in the east direction is smaller than the UWB’s value at the last time index. The dual KF provides the east position with a smaller RMSE than the KF method does at the last time index. The proposed dual cKF method’s east position RMSE is the smallest at the last time index, demonstrating the efficacy of the proposed method. In the north direction, all the filter position RMSEs are similar but smaller than the UWB solutions. The proposed dual cKF method has the smallest localization error, although its value is slightly lower than that of the KF and dual KF.
Figure 16 shows the position error CDF measured by UWB, the KF, the dual KF, and the dual cKF in Test 2. The black solid line denotes the UWB solution’s position error, the blue solid line indicates the position measured by the dual KF, the solution derived from the KF method is represented by the solid green line, while the proposed dual cKF method is indicated by the solid red line. From this figure, it is evident that the proposed method exhibits the smallest position error compared to the other methods at the 0.9 point, indicating its effectiveness in reducing localization errors. However, its value is only slightly lower than the KF and dual KF.
Table 5 shows the RMSEs of the INS, UWB, KF, dual KF, and dual cKF in Test 2. The table shows that the UWB’s solution has better localization performance, with a mean localization error of 0.264 m. The localization error was reduced from 0.264 m to 0.248 m and 0.246 m when using the KF and dual KF methods, respectively. The proposed dual cKF has a mean localization error of 0.242 m, the smallest of all filters.
Table 6 lists the average running times of the KF, dual KF, and dual cKF. The table shows that the average running time of the dual KF and dual cKF subfilters is lower than that of the KF method. We must emphasize that in an actual operation, the two subfilters of the dual KF and dual cKF run in parallel; therefore, the running time of the dual KF and dual cKF can take the larger of the two subfilters. Therefore, the proposed method effectively reduces the running time.

5. Conclusions

This study addresses the issue of varying positioning accuracy in different directions due to environmental factors in inertial-based integrated pedestrian navigation by introducing a dual KF based on a single direction under a CMN scheme. In the proposed method, the data fusion model for the KF fusing the data of INS and UWB was designed initially. Subsequently, the dual INS/UWB integrated model-based cKF was designed, which employed two sub-cKFs to fuse the data independently in the east and north directions. Then, the real test was employed to verify the proposed method’s performance. The empirical results indicate that the proposed approach outperforms the KF in terms of position estimation accuracy and runtime reduction, highlighting its effectiveness.

Author Contributions

Conceptualization, Q.W., Y.X. and M.S.; methodology, Q.W., Y.X. and R.Y.; software, R.Y. and K.L.; validation, R.Y., K.L., J.M. and Y.X.; formal analysis, R.Y. and K.L.; investigation, Q.W., R.Y. and K.L.; resources, Q.W., R.Y., J.M. and M.S.; data curation, Q.W., R.Y. and M.S.; writing—original, Y.X. and R.Y.; draft preparation, Y.X. and M.S.; writing—review and editing, Q.W., Y.X. and R.Y.; visualization, Q.W., Y.X. and R.Y.; supervision, M.S.; project administration, Y.X. and M.S.; funding acquisition, Y.X. and M.S. All authors have reviewed and consented to the published version of the manuscript.

Funding

This research received support from Grant JNSX2023024 from the Strategic Project for Integrated Development between the City of Jinan and Universities.

Data Availability Statement

The article includes the data.

Conflicts of Interest

Author Qingdong Wu was employed by the company Shandong Luqiao Group Company Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

This manuscript uses the following abbreviations:
BEBackward Euler
BNBlind Node
BOTBearing-only Tracking
CDFCumulative Distribution Function
CKFCubature Kalman Filter
CMNColored Measurement Noise
AKKFAdaptive Kernel Kalman Filter
EKFExtended Kalman Filter
GPFGaussian Particle Filter
IMUInertial Measurement Unit
GPSGlobal Positioning System
KFKalman Filter
PFParticle Filter
RTKReal-time Kinematic
INSInertial Navigation System
RMSERoot-Mean-Square Error
RNsReference Nodes
UKFUnseasoned Kalman Filter
UWBUltrawideband
ZUPTZero Velocity UPdaTe

References

  1. Weng, D.; Chen, W.; Ji, S.; Wang, J. Intelligent urban positioning using smartphone-based GNSS and pedestrian network. IEEE Internet Things J. 2024, 11, 22537–22549. [Google Scholar] [CrossRef]
  2. Wang, J.; Liu, J.; Xu, X.; Yu, Z.; Li, Z. A yaw correction method for pedestrian positioning using two low-cost MIMUs. Measurement 2023, 217, 112992. [Google Scholar] [CrossRef]
  3. Qi, L.; Liu, Y.; Gao, C.; Feng, T.; Yu, Y. Hybrid pedestrian positioning system using wearable inertial sensors and ultrasonic ranging. Def. Technol. 2024, 33, 327–338. [Google Scholar] [CrossRef]
  4. Ji, M.; Liu, J.; Xu, X.; Lu, Z. The improved 3D pedestrian positioning system based on foot-mounted inertial sensor. IEEE Sens. J. 2020, 21, 25051–25060. [Google Scholar] [CrossRef]
  5. Xu, Y.; Shmaliy, Y.S.; Ahn, C.K.; Chen, X.; Guo, H.; Zhuang, Y. Blind robust multi-horizon EFIR filter for tightly integrating INS and UWB. IEEE Sens. J. 2021, 21, 23037–23045. [Google Scholar] [CrossRef]
  6. Cao, B.; Jiang, C.; Fan, S.; Zhang, H.; Liu, W. Improving the localization accuracy and robustness of a UWB system using VB-CSRUKF and RTS in harsh underground NLOS environments. IEEE Internet Things J. 2024, 11, 22790–22802. [Google Scholar] [CrossRef]
  7. Iosa, M.; Picerno, P.; Paolucci, S.; Morone, G. Wearable inertial sensors for human movement analysis. Expert Rev. Med. Devices 2016, 13, 641–659. [Google Scholar] [CrossRef]
  8. Shahid, R.; Baloch, A.; Tahir, H.; Ullah, A. Odometry and inertial sensor-based localization of a snake robot. In Proceedings of the 2023 International Conference on Robotics and Automation in Industry (ICRAI), Peshawar, Pakistan, 3–5 March 2023; pp. 1–6. [Google Scholar]
  9. Basso, M.; Galanti, M.; Innocenti, G.; Miceli, D. Triggered INS/GNSS data fusion algorithms for enhanced pedestrian navigation system. IEEE Sens. J. 2020, 20, 7447–7459. [Google Scholar] [CrossRef]
  10. Chen, J.; Zhou, B.; Bao, S.; Liu, X.; Gu, Z.; Li, L.; Zhao, Y.; Zhu, J.; Li, Q. A data-driven inertial navigation/Bluetooth fusion algorithm for indoor localization. IEEE Sens. J. 2021, 22, 5288–5301. [Google Scholar] [CrossRef]
  11. Wu, Q.; Li, C.; Shen, T.; Xu, Y. Improved adaptive iterated extended kalman filter for GNSS/INS/UWB-integrated fixed-point positioning. Comput. Model. Eng. Sci. 2023, 134, 1761–1772. [Google Scholar]
  12. Li, X.; Wu, Z.; Shen, Z.; Xu, Z.; Li, X.; Li, S.; Han, J. An indoor and outdoor seamless positioning system for low-cost UGV using PPP/INS/UWB tightly coupled integration. IEEE Sens. J. 2023, 23, 24895–24906. [Google Scholar] [CrossRef]
  13. Xu, Y.; Wan, D.; Shmaliy, Y.S.; Chen, X.; Shen, T.; Bi, S. Dual free-size LS-SVM assisted maximum correntropy kalman filtering for seamless INS-based integrated drone localization. IEEE Trans. Ind. Electron. 2024, 71, 9845–9854. [Google Scholar] [CrossRef]
  14. Feng, D.; Wang, C.; He, C.; Zhuang, Y.; Xia, X.G. Kalman-filter-based integration of IMU and UWB for high-accuracy indoor positioning and navigation. IEEE Internet Things J. 2020, 7, 3133–3146. [Google Scholar] [CrossRef]
  15. Xu, Y.; Li, Y.; Ahn, C.K.; Chen, X. Seamless indoor pedestrian tracking by fusing INS and UWB measurements via LS-SVM assisted UFIR filter. Neurocomputing 2020, 388, 301–308. [Google Scholar] [CrossRef]
  16. Cheng, L.; Zhao, F.; Zhao, P.; Guan, J. UWB/INS fusion positioning algorithm based on generalized probability data association for indoor vehicle. IEEE Trans. Intell. Veh. 2023, 9, 446–458. [Google Scholar] [CrossRef]
  17. Pan, S.; Jia, F.; Wei, J.; Chen, Q.; Wu, P.; Baoguo, Y. Intelligent vehicle positioning method based on UWB single basestation system / INS tight combination. J. Chin. Inert. Technol. 2023, 31, 1196–1202. [Google Scholar]
  18. Kang, H.; Duan, Y.; Ling, Y. Indoor positioning system based on tight coupling of UWB/INS. J. Xi’an Univ. Posts Telecommun. 2023, 28, 12–20. [Google Scholar]
  19. Chen, K.H.; Chen, P.C.; Liu, K.C.; Chan, C.T. Wearable sensor-based rehabilitation exercise assessment for knee osteoarthritis. Sensors 2015, 15, 4193–4211. [Google Scholar] [CrossRef]
  20. Cicirelli, G.; Impedovo, D.; Dentamaro, V.; Marani, R.; Pirlo, G.; D’Orazio, T.R. Human gait analysis in neurodegenerative diseases: A review. IEEE J. Biomed. Health Inform. 2021, 26, 229–242. [Google Scholar] [CrossRef]
  21. Wong, W.Y.; Wong, M.S.; Lo, K.H. Clinical applications of sensors for human posture and movement analysis: A review. Prosthetics Orthot. Int. 2007, 31, 62–75. [Google Scholar] [CrossRef]
  22. Zhang, Y.; Wang, K.; Jiang, J.; Tan, Q. Research on intraoperative organ motion tracking method based on fusion of inertial and electromagnetic navigation. IEEE Access 2021, 9, 49069–49081. [Google Scholar] [CrossRef]
  23. D’Ippolito, F.; Garraffa, G.; Sferlazza, A.; Zaccarian, L. A hybrid observer for localization from noisy inertial data and sporadic position measurements. Nonlinear Anal. Hybrid Syst. 2023, 49, 101360. [Google Scholar] [CrossRef]
  24. Hu, G.; Xu, L.; Gao, B.; Chang, L.; Zhong, Y. Robust unscented Kalman filter based decentralized multi-sensor information fusion for INS/GNSS/CNS integration in hypersonic vehicle navigation. IEEE Trans. Instrum. Meas. 2023, 72, 1–11. [Google Scholar]
  25. Xu, J.; Yang, G.; Sun, Y.; Picek, S. A multi-sensor information fusion method based on factor graph for integrated navigation system. IEEE Access 2021, 9, 12044–12054. [Google Scholar] [CrossRef]
  26. Xu, Y.; Shmaliy, Y.S.; Bi, S.; Chen, X.; Zhuang, Y. Extended kalman/UFIR filters for UWB-based indoor robot localization under time-varying colored measurement noise. IEEE Internet Things J. 2023, 10, 15632–15641. [Google Scholar] [CrossRef]
  27. Zhao, Y. Cubature + extended hybrid kalman filtering method and its application in PPP/IMU tightly coupled navigation systems. IEEE Sens. J. 2015, 15, 6973–6985. [Google Scholar] [CrossRef]
  28. Sun, M.; Davies, M.E.; Hopgood, J.R.; Proudler, I. Adaptive kernel Kalman filter multi-sensor fusion. In Proceedings of the 2021 IEEE 24th International Conference on Information Fusion (FUSION), Sun City, South Africa, 1–4 November 2021; pp. 1–8. [Google Scholar]
  29. Li, D.; Jia, X.; Zhao, J. A novel hybrid fusion algorithm for low-cost GPS/INS integrated navigation system during GPS outages. IEEE Access 2020, 8, 53984–53996. [Google Scholar] [CrossRef]
  30. Bi, S.; Ma, L.; Shen, T.; Xu, Y.; Li, F. Neural network assisted Kalman filter for INS/UWB integrated seamless quadrotor localization. PeerJ Comput. Sci. 2021, 7, e630. [Google Scholar] [CrossRef]
  31. Chen, C.; Pan, X. Deep learning for inertial positioning: A survey. IEEE Trans. Intell. Transp. Syst. 2024, 1–18. [Google Scholar] [CrossRef]
  32. Revach, G.; Shlezinger, N.; Ni, X.; Escoriza, A.L.; Van Sloun, R.J.; Eldar, Y.C. KalmanNet: Neural network aided Kalman filtering for partially known dynamics. IEEE Trans. Signal Process. 2022, 70, 1532–1547. [Google Scholar] [CrossRef]
  33. Sung, S.; Kim, H.; Jung, J.I. Accurate indoor positioning for UWB-based personal devices using deep learning. IEEE Access 2023, 11, 20095–20113. [Google Scholar] [CrossRef]
  34. Xu, Y.; Shmaliy, Y.S.; Shen, T.; Chen, D.; Sun, M.; Zhuang, Y. INS/UWB-based quadrotor localization under colored measurement noise. IEEE Sens. J. 2020, 21, 6384–6392. [Google Scholar] [CrossRef]
  35. Xu, Y.; Cao, J.; Shmaliy, Y.S.; Zhuang, Y. Distributed kalman filter for UWB/INS integrated pedestrian localization under colored measurement noise. Satell. Navig. 2021, 2, 1–10. [Google Scholar] [CrossRef]
  36. Shmaliy, Y.S.; Zhao, S.; Ahn, C.K. Kalman and UFIR state estimation with coloured measurement noise using backward Euler method. IET Signal Process. 2020, 14, 64–71. [Google Scholar] [CrossRef]
Figure 1. Structure of the Kalman filter (KF) under CMN based on a single direction.
Figure 1. Structure of the Kalman filter (KF) under CMN based on a single direction.
Electronics 13 03027 g001
Figure 2. Structure of the dual Kalman filter (KF) based on a single direction under colored measurement noise (CMN) for inertial navigation system (INS)-based integrated human localization.
Figure 2. Structure of the dual Kalman filter (KF) based on a single direction under colored measurement noise (CMN) for inertial navigation system (INS)-based integrated human localization.
Electronics 13 03027 g002
Figure 3. Actual testing scenarios.
Figure 3. Actual testing scenarios.
Electronics 13 03027 g003
Figure 4. The test structure.
Figure 4. The test structure.
Electronics 13 03027 g004
Figure 5. Reference path and the paths measured using the INS + ZUPT, UWB, KF, dual KF, and the dual cKF in Test 1.
Figure 5. Reference path and the paths measured using the INS + ZUPT, UWB, KF, dual KF, and the dual cKF in Test 1.
Electronics 13 03027 g005
Figure 6. Reference value and the position in the east direction measured using the INS + ZUPT, UWB, KF, dual KF, and the dual cKF in Test 1.
Figure 6. Reference value and the position in the east direction measured using the INS + ZUPT, UWB, KF, dual KF, and the dual cKF in Test 1.
Electronics 13 03027 g006
Figure 7. Reference value and the position in the north direction measured using the INS + ZUPT, UWB, KF, dual KF, and the dual cKF in Test 1.
Figure 7. Reference value and the position in the north direction measured using the INS + ZUPT, UWB, KF, dual KF, and the dual cKF in Test 1.
Electronics 13 03027 g007
Figure 8. UWB distances errors and the position errors of the methods in Test 1.
Figure 8. UWB distances errors and the position errors of the methods in Test 1.
Electronics 13 03027 g008
Figure 9. Position root-mean-square errors (RMSEs) in the east and north directions measured by UWB, the KF, the dual KF, and the dual cKF in Test 1: (a) east direction, (b) north direction.
Figure 9. Position root-mean-square errors (RMSEs) in the east and north directions measured by UWB, the KF, the dual KF, and the dual cKF in Test 1: (a) east direction, (b) north direction.
Electronics 13 03027 g009
Figure 10. Position error cumulative distribution function (CDF) measured by UWB, the KF, the dual KF, and the dual cKF in Test 1.
Figure 10. Position error cumulative distribution function (CDF) measured by UWB, the KF, the dual KF, and the dual cKF in Test 1.
Electronics 13 03027 g010
Figure 11. Reference path and those measured by the INS + ZUPT, UWB, KF, dual KF, and dual cKF in Test 2.
Figure 11. Reference path and those measured by the INS + ZUPT, UWB, KF, dual KF, and dual cKF in Test 2.
Electronics 13 03027 g011
Figure 12. Reference value and the position in the east direction measured by the INS + ZUPT, UWB, KF, dual KF, and dual cKF in Test 2.
Figure 12. Reference value and the position in the east direction measured by the INS + ZUPT, UWB, KF, dual KF, and dual cKF in Test 2.
Electronics 13 03027 g012
Figure 13. Reference value and the position in the north direction measured by the INS + ZUPT, UWB, KF, dual KF, and dual cKF in Test 2.
Figure 13. Reference value and the position in the north direction measured by the INS + ZUPT, UWB, KF, dual KF, and dual cKF in Test 2.
Electronics 13 03027 g013
Figure 14. UWB distances error and the position errors of the methods in Test 2.
Figure 14. UWB distances error and the position errors of the methods in Test 2.
Electronics 13 03027 g014
Figure 15. Position root-mean-square errors (RMSEs) in the east and north directions measured by UWB, the KF, the dual KF, and the dual cKF in Test 2: (a) east direction, (b) north direction.
Figure 15. Position root-mean-square errors (RMSEs) in the east and north directions measured by UWB, the KF, the dual KF, and the dual cKF in Test 2: (a) east direction, (b) north direction.
Electronics 13 03027 g015
Figure 16. Position error CDF measured by UWB, the KF, the dual KF, and the dual cKF in Test 2.
Figure 16. Position error CDF measured by UWB, the KF, the dual KF, and the dual cKF in Test 2.
Electronics 13 03027 g016
Table 1. The parameters of the UWB localization system used in the test.
Table 1. The parameters of the UWB localization system used in the test.
ParameterValue
Data availability frequency1 s
Distance error0.25 m
Localization error0.35 m
Table 2. The parameters of the R-T-K used in the test.
Table 2. The parameters of the R-T-K used in the test.
ParameterValue
Data availability frequency1 s
R-T-K accuracy (plane) ± 8 + 1 × 10 6 × D mm
R-T-K accuracy (height) ± 15 + 1 × 10 6 × D mm
size119 mm × 119 mm × 85 mm
Weigh0.73 kg
Operation temperature−45 ℃ ∼ +75 ℃
Table 3. Root-mean-square errors (RMSEs) of the INS, UWB, KF, dual KF, and dual cKF in Test 1.
Table 3. Root-mean-square errors (RMSEs) of the INS, UWB, KF, dual KF, and dual cKF in Test 1.
MethodsEast Direction (m)North Direction (m)Mean (m)
INS0.2070.9810.153
UWB0.2200.1490.185
KF0.1630.1440.154
dual KF0.1490.1300.140
dual cKF0.1430.1190.131
Table 4. The average running time of the KF, dual KF, and dual cKF in Test 1.
Table 4. The average running time of the KF, dual KF, and dual cKF in Test 1.
MethodsRunning Time (ms)
KF0.16
Sub-KF 10.06
Sub-KF 20.05
Sub-cKF 10.15
Sub-cKF 20.09
Table 5. RMSEs of the INS, UWB, KF, dual KF, and dual cKF in Test 2.
Table 5. RMSEs of the INS, UWB, KF, dual KF, and dual cKF in Test 2.
MethodsEast Direction (m)North Direction (m)Mean (m)
INS0.5810.6810.631
UWB0.2900.2370.264
KF0.2830.2120.248
dual KF0.2810.2110.246
dual cKF0.2740.2100.242
Table 6. Average running time of the KF, dual KF, and dual cKF in Test 2.
Table 6. Average running time of the KF, dual KF, and dual cKF in Test 2.
MethodsRunning Time (ms)
KF0.10
Sub-KF 10.04
Sub-KF 20.04
Sub-cKF 10.08
Sub-cKF 20.06
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, Q.; Yang, R.; Liu, K.; Xu, Y.; Miao, J.; Sun, M. Dual Kalman Filter Based on a Single Direction under Colored Measurement Noise for INS-Based Integrated Human Localization. Electronics 2024, 13, 3027. https://doi.org/10.3390/electronics13153027

AMA Style

Wu Q, Yang R, Liu K, Xu Y, Miao J, Sun M. Dual Kalman Filter Based on a Single Direction under Colored Measurement Noise for INS-Based Integrated Human Localization. Electronics. 2024; 13(15):3027. https://doi.org/10.3390/electronics13153027

Chicago/Turabian Style

Wu, Qingdong, Ruohan Yang, Kaixin Liu, Yuan Xu, Jijun Miao, and Mingxu Sun. 2024. "Dual Kalman Filter Based on a Single Direction under Colored Measurement Noise for INS-Based Integrated Human Localization" Electronics 13, no. 15: 3027. https://doi.org/10.3390/electronics13153027

APA Style

Wu, Q., Yang, R., Liu, K., Xu, Y., Miao, J., & Sun, M. (2024). Dual Kalman Filter Based on a Single Direction under Colored Measurement Noise for INS-Based Integrated Human Localization. Electronics, 13(15), 3027. https://doi.org/10.3390/electronics13153027

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