IMU / Magnetometer / Barometer / Mass-Flow Sensor Integrated Indoor Quadrotor UAV Localization with Robust Velocity Updates

Velocity updates have been proven to be important for constraining motion-sensor-based dead-reckoning (DR) solutions in indoor unmanned aerial vehicle (UAV) applications. The forward velocity from a mass flow sensor and the lateral and vertical non-holonomic constraints (NHC) can be utilized for three-dimensional (3D) velocity updates. However, it is observed that (a) the quadrotor UAV may have a vertical velocity trend when it is controlled to move horizontally; (b) the quadrotor may have a pitch angle when moving horizontally; and (c) the mass flow sensor may suffer from sensor errors, especially the scale factor error. Such phenomenons degrade the performance of velocity updates. Thus, this paper presents a multi-sensor integrated localization system that has more effective sensor interactions. Specifically, (a) the barometer data are utilized to detect height changes and thus determine the weight of vertical velocity update; (b) the pitch angle from the inertial measurement unit (IMU) and magnetometer data fusion is used to set the weight of forward velocity update; and (c) an extra mass flow sensor calibration module is introduced. Indoor flight tests have indicated the effectiveness of the proposed sensor interaction strategies in enhancing indoor quadrotor DR solutions, which can also be used for detecting outliers in external localization technologies such as ultrasonics.


Introduction
Unmanned aerial vehicles (UAV) have shown great potential in civilian applications such as indoor/outdoor mapping [1], target tracking [2], victim searching [3], and industrial inspection [4].For these applications, a key is the real-time estimation of UAV navigation states (i.e., position, velocity, and attitude).Although the integration of data from global navigation satellite systems (GNSS), real-time kinematics (RTK), or precise point positioning (PPP) and inertial measurement units (IMU) have been successfully commercialized to provide accurate (i.e., decimeter and centimeter level location accuracy for dynamic and static applications, respectively) location solutions in outdoor areas [5], reliable indoor UAV localization is a challenge due to the degradation of GNSS signals.
To alleviate this issue, researchers have presented various systems and approaches.Table 1 lists part of the existing works from the years 2016 to 2018.The used sensors and algorithms are shown as well as their test areas and location accuracies.
Table 1.Selected systems and methods for indoor unmanned aerial vehicles (UAV) localization from the years 2016 to 2018.

•
Different types of sensors typically provide various localization accuracies and meanwhile have different costs and coverage areas.Thus, there is a trade-off between performance and cost/coverage.

•
High-precision wireless technologies (e.g., UWB and ultrasonic) can provide high localization accuracy (e.g., decimeter or even centimeter level).However, although the prices for low-cost commercial UWB and ultrasonic development kits have been reduced to the hundreds of dollars level, such systems have limited ranges (e.g., 30 m between nodes and anchors).Thus, other technologies are required to bridge their signal outages in wide-area applications.Meanwhile, for wireless ranging systems, there are inherent issues such as signal obstruction and multipath [36].Thus, other technologies are needed to ensure localization reliability and integrity.

•
Cameras and lidars can also provide high location accuracy when loop closures are correctly detected.Furthermore, some previous issues, such as heavy computational load, are being eliminated because of the development of modern processors and wireless data transmission technologies.However, the performance of vision-based localization systems is highly dependent on whether the measured features are distinct in space and stable over time.
For database matching, any inconsistency between the measured data and the database may cause mismatches [37].For mobile mapping, it is possible to add updates and loop closures to control errors.However, real-world localization conditions are complex and unpredictable; thus, it is difficult to maintain accuracy in challenging environments (e.g., areas with glass or solid-color walls).Therefore, external technologies may be needed to bridge such task periods as well as detect the outliers in vision sensor measurements.

•
Dead-reckoning (DR) solutions from IMUs have been widely used to bridge other localization technologies' signal outages and integrate with them to provide smoother and more robust solutions [38].However, traditional navigation-or tactical-grade IMUs are heavy and costly and thus are not suitable for consumer-level UAVs.Micro-electro-mechanical systems (MEMS) IMUs are light and low-cost, which have made them suitable for low-cost indoor localization.However, low-cost MEMS IMUs suffer from significant run-to-run biases and thermal drifts [39], which are issues inherent to MEMS sensors.Therefore, standalone IMU-based DR solutions will drift over time.Magnetometer measurements can be used to derive an absolute heading update.However, the indoor magnetic declination angle becomes unknown, which makes the magnetometer heading unreliable [40].Thus, it is still important to implement periodical updates to correct DR solutions.

•
Vehicle motion model updates can be used to enhance the navigation system observability [41], especially when there are significant vehicle dynamics (e.g., accelerating or turning).Sensors such as the mass flow and Hall-effect sensors can measure the forward velocity.Meanwhile, it is assumed that the lateral and vertical velocity components are zeroes plus noises when the UAV is being controlled to move horizontally, i.e., the non-holonomic constraint (NHC) [42].Accordingly, 3D velocity updates can be applied.Furthermore, there are other updates, such as the zero velocity update (ZUPT) and zero angular rate update (ZARU) when the UAV is hovering in a quasi-static mode [43].These updates are effective when the actual UAV motion meets the assumption.However, in contrast to land vehicles that are constrained by the road surface, UAVs may suffer from vertical velocity passively during task periods, which degrades the NHC performance.Meanwhile, UAVs may have a pitch angle when moving horizontally, which pollutes the forward velocity measurements.Therefore, some updates are needed to better use the velocity updates.
This research focuses on using low-cost sensors to provide a self-contained DR solution, so as to bridge the signal outages and resist outliers in high-precision localization solutions.Therefore, the above high-precision wireless and vision sensors are not investigated.Since IMU-based DR solutions drift quickly over time, the magnetometers, mass flow sensor, and barometer are also used.Compared to the existing works, better sensor interactions are utilized to enhance the localization solution.The main contributions of this paper are as follows: 1.
Velocity updates have been proven to be effective in constraining DR errors.However, it is observed that the quadrotor UAV may have vertical velocity even when it is controlled to move horizontally.Therefore, the barometer data are utilized to detect height changes and thus determine the weight for the vertical velocity update.

2.
According to the fact that the quadrotor may have a pitch angle when moving horizontally, the pitch angle, which is obtained from IMU and magnetometer data fusion, is used to set the weight of the forward velocity update.

3.
It is observed that the mass flow sensor may suffer from significant sensor errors, especially the scale factor error.Thus, a specific mass flow sensor calibration module is introduced.
This paper is organized as follows.Section 2 illustrates the methodology.Section 3 describes the experimental verification, and Section 4 draws the conclusions.The IMU data are used for predicting the navigation states through the inertial navigation system (INS) mechanization and constructing the EKF system model.The magnetometer data are calibrated and utilized to compute the magnetometer heading, which is further used to build a heading update in the MSL EKF.The mass flow sensor data are calibrated and used in the velocity update in the EKF.The barometer data are used to detect height changes, which is in turn used for quality control (QC) on velocity updates.Meanwhile, the horizontal angles from IMU data are used for QC on velocity updates.The ultrasonic data are used to provide absolute position updates for the MSL EKF, so as to provide a reference location solution.The position data from the MSL EKF are also used to detect the outliers in ultrasonic position solutions.The following subsections separately describe the inertial navigation system (INS)-based attitude, velocity, and position prediction, the magnetometer heading update, the velocity update, the position update, and the EKF computation.

EKF System Model
Inertial navigation is a DR technique in which an IMU is tracked relative to its initial navigation states (i.e., attitude, velocity, and position) from the IMU body frame (i.e., b-frame) to the navigation frame (i.e., n-frame).The INS mechanization processes angular rates and specific forces (or angular and velocity increments) from gyros and accelerometers in the IMU for navigation-state prediction.Refer to Reference [44] for details about the INS mechanization.The predicted navigation states are also used to construct the MSL EKF system model.The INS error model [45] is applied in the continuous system model as ẋms = F ms x ms + w ms where x ms and F ms are the state vector and the dynamics matrix, respectively.w ms is the system noise vector.The elements in the vectors and matrix are where the states δp n , δv n , ψ, b g , and b a are the vectors of position errors, velocity errors, attitude errors, gyro biases, and accelerometer biases, respectively; C n b is the direction cosine matrix (DCM) from the b-frame to the n-frame.f n is the specific force vector projected to the n-frame.The sign [l×] denotes the cross-product (skew-symmetric) form of the 3D vector l = l 1 l 2 l 3 T .w g and w a are noises in gyro and accelerometer readings, respectively; τ bg and τ ba denote the correlation time for b g and b a , and w bg and w ba are the gyro and accelerometer bias driving noises.0 3×3 and I 3×3 are the 3D zero matrix and identity matrix, respectively.

Magnetometer Heading Update
In the INS mechanization, the horizontal attitude (i.e., roll and pitch) errors can be controlled by accelerometers, while the heading error may grow as a result of the weak observability of the heading angle and vertical gyro bias [41].Thus, magnetometers are utilized to provide absolute heading updates.In indoor environments, the local magnetic field may be disturbed by man-made infrastructures.Therefore, calibration is required to use the magnetometer as a reliable source of heading.Refer to References [46,47] for magnetometer calibration and magnetometer heading computation, respectively.The calculated magnetometer heading are also fused with INS data in the attitude and heading reference system (AHRS) algorithm [48] to obtain the AHRS heading updates.The obtained heading is used to build the heading update model in the MSL EKF.The corresponding measurement model can be written as where z ψ,k , H ψ,k , and v ψ,k are the observation vector, design matrix, and measurement noise vector for heading update at time t k .These vectors and matrix can be described as where ψ m,k is the magnetometer-or AHRS-derived heading at time t k , and ψ in,k is the INS-predicted heading.

Velocity Update
It has been revealed in previous research that using a velocity update can enhance the heading observability when there are significant vehicle dynamics [42].As illustrated in Section 1, this paper uses a velocity update and presents multiple updates for QC of velocity updates.

Velocity Update for Multi-Sensor Localization EKF
In this research, 3D velocity updates are applied.The forward velocity is measured by the mass flow sensor, while the lateral and vertical velocity components are set at zeroes plus noises.The corresponding MSL EKF measurement model can be written as where z vc,k , H vc,k , and v vc,k are the observation vector, design matrix, and measurement noise vector for velocity update at time t k .According to Reference [49], z vc and H vc can be written as where ṽb = v f 0 0 T and v f is the velocity from mass flow sensor.

Mass Flow Sensor Calibration
When the reference velocity form external technologies (e.g., ultrasonic) is available, the mass flow sensor can be calibrated.The velocity error model is ṽ where ṽ is the reference velocity.b f and s f are the bias and scale factor of the mass flow sensor, respectively.Thus, the mass flow sensor calibration model is where z f c , H f c , and x f c are the observation vector, design matrix, and state vector for mass flow sensor calibration, and where N f c is the number of mass flow sensor measurements.The least squares method can be used to estimate x f c by The estimated mass flow sensor bias and scale factor errors are used to compensate for velocity errors before using the velocity update.

Availability for the Velocity Update
A QC mechanism is used to improve the robustness of velocity measurements.The basic ideas for QC include (1) the weight for the vertical velocity update is lowered when the vehicle has a significant height change and (2) the weight for the forward velocity update is decreased when the vehicle pitch angle is large.For (1), the barometer data are used for height change detection.With the model in Reference [50], the barometer-measured air pressure can be converted to the barometer height as (16) where h b,k is the barometer height at time t k , p k and p 0 are respectively the measured air pressure and the sea level reference pressure.The p 0 value is set at 101,325 Pa for calculation.
The height data are typically noisy.Thus, a smoother is used as where n bs is the smoother window size.c 1 to c n bs are the coefficients that meet the condition ∑ n bs i=1 c i = 1 and c 1 ≥ ... ≥ c n bs .This smoother causes a lag of n bs 2 data epochs.Such a lag is acceptable for indoor UAVs because UAV sensors typically have a high data rate (e.g., over 50 Hz).A time lag of within 0.1 s will occur if n bs is set at 10.
The smoothed barometer height is used for tracking height changes through a threshold-based piecewise model.To use the method, a scoring value, α hc , is computed and compared with the corresponding threshold values.Specifically, when α hc ≤ T hc,1 (i.e., the quasi-static-height mode), the vertical velocity measurement noise covariance is set at σ 2 vv .When T hc,1 < α hc ≤ T hc,2 (i.e., the low-height-change mode), the vertical velocity measurement noise covariance is set at α hc T hc,1 When α hc > T hc,2 (i.e., the high-height-change mode), the vertical velocity measurement noise covariance is set at a large value σ 2 vmax .In this situation, the vertical velocity update will not contribute to the solution.To use the method, the standard deviation (STD) of the n h smoothed height data epochs is computed as the scoring value α hc .The threshold values are set based on training data.
For (2), a similar threshold-based piecewise model is used by comparing the real-time pitch angle to the preset threshold.Specifically, when θ ≤ T θ,1 (i.e., the horizontal mode), the forward velocity measurement noise covariance is set at σ 2 v f .When T θ,1 < θ ≤ T θ,2 (i.e., the low-pitch-angle mode), the forward velocity measurement noise covariance is set at When θ > T θ,2 (i.e., the high-pitch-angle mode), the forward velocity measurement noise covariance is set at a large value σ 2 vmax .

Position Update
Although the integration of data from the IMU, magnetometers, barometer, and mass flow sensor can provide a short-term accurate DR solution, the solution will drift over time when an absolute update is not available.The drifts of the navigation solution will occur especially when low-cost MEMS sensors are used, since these sensors are susceptible to significant run-to-run biases and thermal drifts.Thus, to obtain a long-term accurate navigation solution, at least one type of the absolute updates is required.In this research, the device-anchor ranges from ultrasonic sensors are used to integrate with MSL EKF, so as to obtain references for the localization solution.
Ultrasonic sensors can provide accurate (i.e., centimeter-level) ranges in environments that have clear line-of-sight between the device and anchors.However, the ranging accuracy may be degraded by outliers that may occur as a result of obstructions between the device and anchors and the multipath effect.This section describes the method for localization using ultrasonic ranges, the position measurement model for the MSL EKF, and the method for removing outliers in ultrasonic ranges.

Ultrasonic Multilateration
For 3D localization, the model for the range between the device and the i-th anchor is Therefore, the multilateration model is where x r = δx δy δz represents the vector of device location errors, z r is the observation vector for multilateration, H r is the design matrix for multilateration, and The state vector x r is estimated by where x r = x r y r z r T and x r , y r , and z r are the device coordinates along east, north, and up directions.

Position Update for Multi-Sensor Localization EKF
The ultrasound-derived position solutions are used to build the MSL EKF position update model as where z pc,k , H pc,k , and v pc,k are the observation vector, design matrix, and measurement noise vector for position update at time t k , and where xk ỹk zk T is the multilateration solution using ultrasonic ranges at time t k .χ k β k h k T are the elements in p n .χ re f β re f h re f T is the global position (i.e., latitude, longitude, and ellipsoidal height) of the original point for the local coordinate frame (i.e., x, y, and z).

Ultrasonic Position Outlier Detection
To detect outliers in ultrasonic positioning solutions, statistical testing is implemented on the innovations of the MSL EKF.The innovation sequence and its covariance matrix can be calculated by [51] where P ms,k is the state covariance matrix in the MSL EKF and R pc,k is the position measurement covariance matrix.The assumption for outlier detection is where represents the normal distribution with mean of c 1 and covariance of c 2 .If the hypothesis is rejected, the measurements corresponding to the outliers are removed.
With the system model and measurement models, the EKF predicts the states and then obtains updates from noisy measurements.Refer to Reference [45] for details of EKF computation.

Test Description
To verify the proposed MSL method, indoor flight tests were conducted with a 3DR Solo quadrotor [52].Figure 2a illustrates the test environment (20 m × 20 m) and Figure 2b shows the quadrotor and sensors.Five Marvelmind ultrasonic (abbreviated as US) beacons [53] were used, including one fixed on the quadrotor and four installed on four static leveling pillars.The height of the four static beacons was 4 m.The quadrotor was equipped with an InvenSense MPU6000 MEMS-based IMU [54], a Honeywell HMC 5983 magnetometer (abbreviated as Mag) triad [55], a TE MS5611 barometer (abbreviated as Baro) [56], a Sensirion SFM3000 mass flow sensor (abbreviated as Flow) [57], and a Marvelmind US beacon.The data rates for IMU, Mag, Baro, Flow, and US were 50, 100, 100, 100, and 100 Hz, respectively.A LattePanda 4 GB/64 GB Windows 10 single board computer was used for data collection and sensor fusion computation [58].
Four flying tests with various quadrotor trajectories were conducted.The trajectories are shown by blue lines in Figure 3.Each trajectory lasted for five to ten minutes.The locations of four US anchors are shown by red pins in Figure 3.The reference trajectories were obtained by post-processing of fusing US with INS, Mag, Baro, and Flow data.The used US system can provide a centimeter-level ranging accuracy in line-of-sight environments [53].

Impact of Velocity Solutions
In this subsection, the impact of the velocity update is tested.Meanwhile, other factors, such as the detection of height changes and the pitch angle, are investigated.Finally, the AHRS/INS/Velocity integrated solutions that use various velocity strategies are evaluated.

Velocity Solutions (Mass Flow-Based)
To investigate the effect of mass flow sensor calibration, Figure 4a illustrates the raw and low-pass filtered mass flow data, as well as the reference 1D velocity from ultrasonic solutions.There was a difference between the filtered and reference data, which indicated the existence of mass flow sensor errors.Such data were used as training data for mass flow sensor calibration.
The estimated mass flow sensor bias and scale factor values were further used to compensate the mass flow data in the tests.The corresponding 1D velocity time series are shown in Figure 4b.Compared to the filtered velocity, the calibrated velocity was closer to the reference.This phenomenon indicates the effectiveness of mass flow sensor calibration.

Height-Change Detection (Barometer-Based)
The barometer height was utilized to detect the time periods that did not have a significant height change.These time periods were important for using the vertical velocity update.Figure 5 illustrates the raw and smoothed barometer heights, as well as the indicator for the time periods that did not have a significant height change.The time periods indicated by the cyan dots that had a legend of Flag-H.
The data in Figure 5 indicate that even when the quadrotor is controlled to move horizontally, it may have height changes.With the height-change detection technique, it is possible to process the barometer data to obtain the time periods during which the quadrotor had only horizontal movements.
According to training data, the n h value was set at 50 (i.e., barometer data epochs in half a second), the height-change threshold values T hc,1 and T hc,2 were set at 0.02 m and 0.06 m, respectively.The vertical velocity covariance values σ 2 vv and σ 2 vmax were set at (0.1 m/s) 2 and (100 m/s) 2 , respectively.

Impact of Pitch Angle on Velocity
In contrast to land vehicles, quadrotors may have significant horizontal angles, especially the pitch angle, during the flying process.Figure 6a illustrates the roll and pitch angle during a test.The pitch angle reached 15 deg. Figure 6b shows the theoretical relationship between the pitch angle and the forward velocity scale factor error.A scale factor error of around 4 % may be introduced by a pitch angle of 15 deg.
According to the training data, the pitch-angle threshold values T θ,1 and T θ,2 were set at 10 deg and 30 deg, respectively.The forward velocity covariance value σ 2 v f was set at (0.3 m/s) 2 .    2 illustrates the location error statistics, including the mean, root mean squares (RMS), 80% and 95% quantile errors, and the maximum value.Compared to AHRS/INS/Flow(Raw), the mean AHRS/INS/Vel(Raw) location error was reduced from 20.2 m to 15.9 m, with an accuracy improvement of 21.3%.This phenomenon indicates the benefits of using 3D velocity instead of 1D.
When using the calibrated mass flow data, the location error mean value in AHRS/INS/Vel(Cali) was reduced to 10.6 m, which was 33.3% lower than that in AHRS/INS/Vel(Raw).Such accuracy improvement is significant, which suggest calibration for the mass flow sensor.The AHRS/INS/ Vel(Cali,QC) strategy further reduced the mean location error from 10.6 m to 9.4 m, with an accuracy improvement of around 10%.

Use of Ultrasonic Positioning
The solutions in Section 3.2 indicate the effectiveness of using more reliable velocity solutions in enhancing AHRS/INS/Velocity integrated navigation.On the other hand, it is shown that it was challenging to obtain a long-term DR solution with such an AHRS/INS/Velocity integrated system on a quadrotor.Thus, external updates are still needed.Figure 9 shows the US position solution that was obtained through the method in Section 2.4.1.Figure 9a demonstrates the time series of the ranges from the quadrotor to four anchors.Figure 9b,c shows the 2D location solution and the east and north positions, and Figure 9d illustrates the height solution.
The US system can provide a centimeter-level ranging accuracy [53]; thus, the multilateration solution was expect to be generally at a centimeter to decimeter level, which was at least one level more accurate than the DR position accuracy (sub-meter to meter level).One issue for the US positioning system is that there were outliers in range measurements, which caused outliers in the US position solution in Figure 9a.

AHRS/INS/Velocity/Ultrasonic Integrated Solution
To mitigate the effect of outliers, the US positioning solution was utilized as the position update for the MSL EKF by following the method in Section 2.4.2.The outliers were detected and removed by using the approach in Section 2.4.3.Figure 10a-d demonstrate the 2D locations, east and north positions, 1D velocity magnitude, and the east and north velocity components from the MSL EKF, respectively.Compared to the solution in Figure 9, the position solutions became smoother, with all position outliers removed.The MSL EKF position solutions were used as the reference trajectories to evaluate the MSL solution during US signal outages in the next subsection.Figure 11a-d shows the 2D position results by processing four sets of test data with US outages of 5, 10, 30, and 60 s, respectively.The red dots indicates the solutions during US outage time periods.Figure 12a illustrates the east and north positions with and without US outages, and Figure 12b shows the corresponding location errors, which were the differences between the results with and without US outages.Thus, in Figure 12b, the location errors during the time periods that had US position updates were not used when computing the location error statistics in this subsection.Figure 13a illustrates the CDF of 2D location errors when using US outage time periods of 5, 10, 15, 20, 30, and 60 s on all test trajectories.Only the location errors during US outage periods were used for calculation.Figure 13b shows the corresponding statistics.
According to Figure 13b, the mean 2D location errors were 0.2, 0.6, 1.0, 1.3, 1.8, and 4.3 m when there were US outages of 5, 10, 15, 20, 30, and 60 s, respectively.This solution indicated that the AHRS/INS/Velocity integrated system generally provided a localization accuracy (in mean value) of approximately 1.0 m and 2.0 m when localizing using AHRS/INS/Velocity integration for 15 and 30 s, respectively.Such location accuracy is acceptable because one-meter accuracy for 15 s is enough for resisting many signal interference and outages for a commercial US system.Particularly, the location accuracy reached 0.2 and 0.6 m when localizing using AHRS/INS/Velocity integration for 5 and 10 s, respectively.Such accuracy was promising for self-contained quadrotor localization without using high-precision wireless or vision localization technologies.

Figure 1
Figure 1 illustrates the system diagram for the proposed multi-sensor integrated localization (MSL) method.The blue and red boxes indicate the inputs and outputs, respectively.The purple boxes represent the algorithm modules, while the green boxes indicate the prediction and update information for the MSL extended Kalman filter (EKF).

Figure 1 .
Figure 1.Diagram for the proposed multi-sensor integrated localization method.

Figure 4 .
Figure 4. (a) Raw, filtered, and calibrated mass flow training and (b) testing data.

Figure 5 .
Figure 5. Barometer data and height-change detection solution.

Figure 6 .
Figure 6.(a) Horizontal angles and (b) theoretical relation between pitch and forward velocity scale factor error. 3.2.4.AHRS/INS/Velocity Integrated Solutions with Various Velocity Strategies Figure 7 shows the position solutions from the following strategies: • AHRS/INS: integration of AHRS heading and INS mechanization, without using any velocity update.• AHRS/INS/Flow(Raw): using raw mass flow sensor data (i.e., 1D velocity) as the update in the MSL EKF.• AHRS/INS/Vel(Raw): using raw mass flow sensor data and NHC for 3D velocity updates in the MSL EKF.• AHRS/INS/Vel(Cali): using calibrated mass flow sensor data and NHC (i.e., 3D velocity) in the MSL EKF.• AHRS/INS/Vel(Cali,QC): using mass flow sensor data that were calibrated and had QC based on height-change and pitch-angle detection, as well as NHC (i.e., 3D velocity) in the MSL EKF.

Figure 7
Figure 7 indicates the importance of using proper velocity updates on motion-sensor-based quadrotor localization.Figure 8a,b shows the location error time series of the strategies that had velocity updates in one test, and the cumulative distribution function (CDF) of location errors in all four tests.Table2illustrates the location error statistics, including the mean, root mean squares (RMS), 80% and 95% quantile errors, and the maximum value.

Figure 8 .
Figure 8.(a) AHRS/INS/Velocity integrated location errors with various velocity strategies and (b) their cumulative distribution function (CDF).

3. 3 . 3 .
AHRS/INS/Velocity Integrated Solution during US Outages To focus on AHRS/INS/Velocity localization, the US positions were cut off to generate US outages.The US outage time length was set at various values (e.g., 5, 10, 15, 20, 30, and 60 s), so as to investigate the AHRS/INS/Velocity solution for different time periods.The AHRS/INS/Velocity system will be valuable in engineering practices if it can provide a reliable solution during US position outages.

Figures 11 and 12
indicate that the location errors drifted over time during US outage periods.The maximum drifts reached 1.3, 3.0, 7.2, and 17.7 m in the selected test data that had US outages of 5, 10, 30, and 60 s, respectively.

Figure 11 .
Figure 11.2D position results by processing four test data with US outages of (a) 5 s, (b) 10 s, (c) 30 s, and (d) 60 s, respectively.

Figure 13 .
Figure 13.(a) CDF and (b) statistics of 2D location errors during ultrasonic outage periods.
This paper has investigated the integration of low-cost IMU, magnetometer, barometer, and mass flow sensors for quadrotor UAV localization.Multiple indoor flying tests were conducted using ultrasonic ranging measurements to compute the reference trajectories.The introduction of forward velocity from the mass flow sensor improved the AHRS/INS-based DR location accuracy (in mean value) by 95.0%, and the use of 3D velocity updates further enhanced the location accuracy by 21.3%.Furthermore, the calibration of mass flow sensor improved the location accuracy by 33.3%, and the sensor interaction strategies further enhanced the location accuracy by 11.3%.The proposed AHRS/INS/Velocity integrated approach generally provided a localization accuracy of 0.2, 0.6, 1.0, 1.3, 1.8, and 4.3 m when localizing for 5, 10, 15, 20, 30 and 60 s, respectively.Such DR accuracy (1.0 m for 15 s and 1.8 m for 30 s) was promising to bridge the signal outages from high-precision localization technologies and to resist outliers in the high-precision localization data.

Table 2 .
Statistics of AHRS/INS/Velocity integrated location errors, with various velocity strategies.