Displacement responses of a bridge as a result of external loadings provide crucial information regarding structural integrity and current conditions. Due to the relative characteristic of displacement, the conventional measurement approach requires reference points to firmly install the transducers, while the points are often unavailable for bridges. In this paper, a displacement estimation approach using Kalman filter-based data fusion is proposed to provide a practical means for displacement measurement. The proposed method enables accurate displacement estimation by optimally utilizing acceleration and strain in combination that have high availability and are free from reference points for sensor installation. The Kalman filter is formulated using a state-space model representing the double integration of acceleration and model-based strain-displacement relationship. The validation of the proposed method is conducted successfully by a numerical simulation and a field experiment, which shows the efficacy and accuracy of the proposed approach in bridge displacement measurement.
Ministry of Oceans and Fisheries201101711. Introduction
Structural health monitoring (SHM) has been considered as an essential procedure for the sustainability of bridges which are the national assets supporting reliable and efficient operation of our society. Technology advances in SHM research have expanded beyond theoretical development and laboratory-scale verification to real-world applications [1–3]. SHM is generally composed of three processes: (1) acquisition of structural responses and environmental factors, (2) assessment of structural condition based on the measurement, and (3) maintenance action when required based on the assessment. Because the first task builds the foundation of the whole SHM procedure, data acquisition is a key component for successful SHM. Before acquiring the structural responses, the types of physical quantities to be monitored are determined priorly with considering necessity and availability of responses. The availability generally includes the accuracy, accessibility, convenience, and cost in the measurement.
Acceleration and strain are the most widely measured responses for the SHM due to their high availability. Acceleration is generally sensitive in measuring dynamic characteristics of a large structure with high convenience in the installation, and strain can be directly converted into current stress condition resisted by the instrumented structural member. Though acceleration and strain possess the majority of the measured responses for the SHM, other responses are often required in practice with specific purposes, including displacement [4–6], inclination [7], structural impedance [8], and pH [9].
Displacement is an intuitive response generated by external forces applied to a structure. Displacement has been used as a practical measurement that indicates structural stability and soundness. Indeed, a variety of modern design codes adopt limit displacement levels under given loadings to assure structural safety. For example, Korea and the US limit maximum deflections of vehicular bridges up to 1/800 of the bridge span length [10, 11] and Canada uses the maximum deflection and the first natural frequency for the evaluation of pedestrian bridges [12]. The practical usage of displacement in the design has been extended to various efforts to use it in the SHM of full-scale civil engineering structures [13–15].
Despite the intuitional characteristic and high demands in practice, displacement has not been as popular as acceleration and strain. Displacement is a relative measurement whose transducer must be placed between a measurement point and a stationary reference point. The sensor installation on bridges is often challenging because bridges are generally constructed over a road, a river, a swamp, or a sea, where a proper reference point is difficult to find. Though a number of noncontact transducers have been developed to resolve the issue, such as global positioning system (GPS) [16–18], LASER-based transducer [19, 20], and computer vision system [21–24], they are still unpopular due to the high cost and synchronization issues when densely deployed.
A noticeable alternative is indirect estimation approaches that convert other structural responses to displacement based on their physical relationships. The convertible responses include acceleration [25–27], strain [28–31], and inclination [32]. Knowing the fact that these responses can be measured without the reference point, the indirect estimation approaches have quite strong potential to be widely used in field applications. However, the indirect methods generally suffer from the low accuracy due to numerical errors and/or imperfect physical relationships. For example, the conversion of acceleration to displacement involves large low-frequency drift errors occurring in a numerical integration [25]; the conversion of strain exhibits errors due to an assumed mapping model between strain and displacement [31]. Thus, the indirect estimation approach has not been commonly accepted in the SHM field yet.
The accuracy enhancement can be achieved by employing multimetric (i.e., a heterogeneous mixture of multiple responses) data based on the data fusion model, as recently reported in [33–36]. The synergic usage of multimetric data minimizes the drawback of each sole response as observed in the indirect displacement estimation methods in the literature. Roberts et al. reported that the conversion of acceleration could be enhanced by referencing displacement measured by a high-precision GPS in the Kalman filtering [37]. Another effort for employing multimetric data is a data fusion method that uses acceleration and strain in combination [38–42]. The method has showed significantly improved accuracy, though a low-frequency noise in strain measurement still brings undesired errors in the estimation. In the aerospace engineering field, the Kalman filter-based data fusion has been actively studied for the accuracy enhancement of vehicle attitude and position assessment [43, 44], autonomous navigation of unmanned vehicles [45, 46], and indoor positioning [47, 48]. Many of them are available on the market by the company such as Trimble Navigation Limited [48] and SBG Systems [49].
In this paper, a Kalman filter-based indirect estimation method is proposed for convenient but accurate measurement of displacement on bridges. The proposed method uses both acceleration and strain, which have the highest availability on bridges without the reference point issue. Furthermore, acceleration and strain have sensitivity in relatively high- and low-frequency range, respectively; the fusion of acceleration and strain accurately estimates displacement in a broad frequency region by taking the synergic effect. The Kalman filter has been selected as the data fusion model due to its high reliability and versatility. A mathematical formulation of the Kalman filter for the displacement estimation is presented based on the double integration of acceleration and modal-mapping of strain. The performance and efficacy of the proposed method are validated in a numerical simulation using a simply supported beam model and further investigated in a field experiment on a single-span prestressed concrete bridge.
2. Formulation for Displacement Estimation Using Kalman Filter
This section provides a mathematical formulation for the proposed displacement estimation approach using the Kalman filter to use acceleration and strain measurements in combination. The fundamental idea of the proposed method is to optimally extract necessary information for displacement estimation from two different measurements, which the sole use of each measurement is not able to do. In general, acceleration responses have rich information in the high-frequency region, whereas strains are sensitive in the low-frequency and static regions. The numerical double integration of an acceleration measurement yields the corresponding displacement at the same position with accurate dynamic components, while resulting in the large low-frequency drift. This error can be compensated by introducing the strain-displacement relationship, which generally does not involve the numerical integration. Thus, selectively utilizing information contained in each measurement can lead to accurate displacement estimation in all frequency regions.
The Kalman filter for displacement estimation needs to be carefully designed to enable such optimal selective filtering. In this study, the formulation in [33] to fuse acceleration and displacement sampled at different rates is modified to accommodate acceleration and strain measurements. The state-space model for the acceleration-displacement relationship can be written using the definition of acceleration, which is the second derivative of displacement with respect to time, as(1)x˙tx¨t=0100xtx˙t+01at+01wt,ut=10xtx˙t+vt,where xt,a(t), u(t), w(t), and v(t) are state variable, acceleration, displacement, and measurement noises associated with acceleration and displacement, respectively. w(t) and v(t) are random processes with the Gaussian distribution with covariance of q and r. Because all measured signals are discrete, a state-space model in the discrete-time domain is desired. The discrete version of (1) can be written as(2)xk=1dt01xk-1+dt22dtak+wk-1=Axk-1+Bak+wk-1,uk=10xk+vk=Hxk+vk,where ak, uk are the measured acceleration and displacement, respectively, and wk and vk are the corresponding noise processes. The noise processes have the covariance matrices defined as(3)Q=qdt33dt22dt22dt,R=rdt,where dt is the sampling time.
The discrete-time state-space model in (2) and (3) is then utilized to develop the Kalman filter for displacement estimation. The measured acceleration is used to project the state ahead via the transition matrix A as(4)x^k-=Ax^k-1+Bak,Pk-=APk-1AT+Q,where x^k- and x^k are respective prior and updated state estimates and Pk- and Pk are respective prior and updated error covariance matrices. Subsequently, the Kalman gain Kk is computed as(5)Kk=Pk-HTHPk-HT+R-1.The measured strain is then employed to update the estimate with the Kalman gain obtained in (5): (6)x^k=x^k-+Kkαuk-Hx^k-,Pk=1-KkHPk-,where uk is the displacement converted from measured strain data and α is an error correction factor for uk [38]. The updated state yields the estimated displacement as(7)u^k=10x^k.This Kalman filter formulation for the displacement estimation is based on the state-space model that describes the definition of acceleration in the time updated in (4) and the strain-displacement relationship in the measurement update in (6). While the acceleration definition as the second derivative of displacement is clearly seen in (4), the strain-displacement relationship needs to be determined depending on structures and measurement locations.
The conversion relationship from strain to displacement can be obtained from a numerical model of a structure [28–31]. Displacement and strain responses can be written in the modal coordinates as(8)u~t=∑iϕiqit=Φq~tε~t=∑iψiqit=Ψq~t,where u~ and ε~ are respective displacement and strain vectors; ϕi and ψi are respective displacement mode shapes of displacement and strain in the ith mode; Φ and Ψ are respective mode shape matrices of displacement and strain; qi is the ith modal coordinate; and q~ is the modal coordinate vector. The conversion relationship is obtained as(9)u~=ΦΨ†ε~,where † denotes the matrix pseudoinverse. Given a finite element model of a bridge, the conversion relationship in (9) is readily available. One particular and common type of bridges is a single-span bridge with homogeneous sectional properties and simply supported boundary conditions. The mode shapes can be reasonably assumed to be sinusoidal functions, which allow (9) to be written as [31](10)u~=L2yπ2Φsinπz1L⋯r2sinrπz1L⋮⋱⋮sinπznL⋯r2sinrπznL†ε~,where y is the location of the neutral axis, r is the maximum number of natural modes being used, z1,…,zn are the locations of strain measurements, and L is the beam length. To compensate the inevitable discrepancy between physical structures and the numerical model such as the finite element model and the simple-beam model, the error correction factor α is employed. More detailed information regarding α can be found in [38].
The proposed displacement estimation method has a distinct advantage in accurate displacement estimation in both low- and high-frequency regions with a minimized noise floor. The displacement contained in the state estimate in the time update of (4) inherits the frequency content of the measured acceleration, which has rich information in the high-frequency region. Estimated from the acceleration, the displacement has a low-frequency noise floor, whereas the large low-frequency drift is involved. The measurement update in (6) is to compensate this low-frequency drift error by employing the strain-displacement relationship. Thus, the resulting displacement of the Kalman filter shall be accurate in all frequency regions. In the following numerical example, the performance of the proposed approach will be verified regarding these aspects of accuracy in both time and frequency domains.
3. Numerical Validation
To validate the proposed displacement estimation method, a beam model with a moving load shown in Figure 1 is considered. The beam is modeled as simply supported with 16 Euler beam elements, and its property is tabulated in Table 1. The moving load is introduced to simulate truck loading which is often used in bridge testing. The load moving at the speed of 0.1 m/s is designed to have a dynamic variance in time using a normal distribution (mean: 10 N, standard deviation: 3 N). Time histories of strain and acceleration responses from the beam under the moving load are simulated at nodes 5, 9, and 13 shown in Figure 1. 5% and 10% RMS noise signals are introduced to acceleration and strain, respectively. Displacement responses at the nodes are also directly calculated from the beam model to verify the accuracy of the estimation using acceleration and strain data.
Beam properties.
Properties
Values
Length
1.6 m
Thickness
4 cm
Width
15 cm
Elastic modulus
206 GPa
Mass density
7850 kg/m3
Beam model.
The performance of the proposed method is evaluated by comparing the estimated displacements with those obtained by the strain-only method using the strain-displacement relationship [28, 31] and the reference displacements. As the simply supported beam is considered in this numerical example, (10) is utilized to obtain the strain-displacement relationship for the strain-only and the proposed methods. A model error is introduced by intentionally using 2.6 cm as the location of the neutral axis (i.e., y of (10)), of which the true location is 2 cm. As previously described, the error correction factor α to compensate this model error [38] is used in the estimation process.
The comparison for displacements at nodes 5 and 9 is shown in Figures 2 and 3, respectively. Figures 2(a) and 3(a) show that the displacements obtained by the strain-only and the proposed method apparently agree well with the reference displacements. However, the close-up views (i.e., Figures 2(b) and 3(b)) and the frequency domain data obtained (i.e., Figures 2(c) and 3(c)) obviously show that the proposed method exhibits better performance in both the low- and high-frequency regions. The close-up views (i.e., Figures 2(b) and 3(b)) from 9.5 to 11.5 seconds show that the proposed method estimates the amplitude of the displacement (related to low-frequency accuracy) more accurately than the strain-based method. This is because of the intentional introduction of erroneous neutral axis to simulate the possible model error to develop the strain-displacement relationship. In the proposed method, the error correction factor α, which is only available when both acceleration and strain are measured together, is used to calibrate the amplitude of displacement, whereas the sole strain measurement does not allow calculating α. Meanwhile, the frequency domain data obtained by fast Fourier transform (i.e., Figures 2(c) and 3(c)) shows that the noise floor (related to high-frequency accuracy) of the proposed method is lower than that of the strain-only method.
Estimated and reference displacement at node 5: (a) time domain data, (b) close-up view of (a), and (c) frequency domain data.
Estimated and reference displacement at node 9: (a) time domain data, (b) close-up view of (a), and (c) frequency domain data.
The quantitative estimation error, defined in terms of ur and uε as (11), is calculated for the obtained displacements shown in Figures 2 and 3:(11)Eerr%=σur-uεσur×100.The estimation errors Eerr of the proposed method are 7.4% and 6.0% for the displacements at nodes 5 and 9, respectively. These error levels can be considered to be reasonably low and agree with the graphical observations from Figures 2 and 3.
4. Experimental Validation on a Test Bridge
The experimental validation of the proposed method is carried out on the prestressed concrete bridge shown in Figure 4. The bridge has a single span with four prestressed concrete girders, and its span length is 11 m. A truck of 28.93 ton is used to load the bridge. The truck runs on the bridge with two different speeds: (1) 5 km/h and (2) 15 km/h. The bridge is instrumented to measure the reference displacement as well as acceleration and strain for the proposed method. The reference displacements are measured by two laser displacement sensors at Locations B and C marked in Figure 4. For the proposed method, three strain gauges are installed at Locations A, C, and D and two accelerometers are placed at the same locations where the displacements are to be estimated and compared with the reference displacements (i.e., Locations B and C).
Prestressed concrete bridge under the truck loading and sensor configuration.
The displacements at Locations B and C are estimated by the proposed method. The neutral axis of the bridge is assumed to be 0.25 m for (10), whose error is to be compensated by the calibration factor α as described in Section 2. The displacements estimated at two locations by the proposed method are compared with the references measured by the laser displacement sensors in Figures 5 and 6 for both truck speeds. The estimated displacements show a good correlation with the reference displacements for both locations and both truck speeds. The errors are quantified between 7.6 and 8.3% by (11) as tabulated in Table 2. As to the location, the error is smaller at the center of the bridge (i.e., Location C), which seems substantial considering larger amplitude of the displacement at Location C. The errors are not seriously affected by the different truck speeds, which may infer different loading patterns. The errors are slightly larger than those obtained in the numerical simulation (i.e., 7.4% and 6.0%). Considering that the experiment is carried out on a large structure, the errors evidently show the good performance of the proposed method. Therefore, the proposed method can obtain structural displacement with high accuracy using highly available responses (i.e., acceleration and strain), especially at the bridges where the conventional displacement transducers are hard to be instrumented.
Error measurements at Locations B and C.
Truck speed
Eest (%)
Location B
Location C
5 km/h
8.33
7.85
15 km/h
8.00
7.63
Comparison of estimated and reference displacements (truck speed: 5 km/h): (a) displacements at Location B and (b) displacements at Location C.
Comparison of estimated and reference displacements (truck speed: 15 km/h): (a) displacements at Location B and (b) displacements at Location C.
5. Conclusions
In this paper, a Kalman filter-based approach using multimetric responses (i.e., acceleration and strain) that have reference-free nature was proposed for indirect displacement estimation of bridges. The formulation of the Kalman filter was constructed using a state-space model representing the double integration of acceleration and model-based strain-displacement relationship. The validation was carried out successfully by a numerical simulation and a field experiment. The results of the validation can be summarized as follows:
In the numerical simulation, the proposed method estimated the nonzero-mean displacements at two locations with good agreement to the reference displacements in both time and frequency domains. The quantitative estimation errors at two locations are 7.4 and 6.0%.
In the experimental test on the prestressed concrete bridge, the proposed method estimated the displacements excited by a truck with the quantitative estimation error between 7.6 and 8.3%. The slightly larger errors than those from the numerical simulation evidently showed good performance of the proposed method on real full-scale bridges.
The errors were not affected by the truck speeds (representing loading patterns) but amplitude of the displacement.
The proposed method is expected to be advantageous for bridge structures where the reference points to stably locate conventional displacement transducers are unavailable.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This research was a part of the project titled “Development of Active-Controlled Tidal Stream Generation Technology” funded by the Ministry of Oceans and Fisheries, Korea (20110171).
WangH.TaoT.LiA.ZhangY.Structural health monitoring system for Sutong cable-stayed bridge201618317334ChenC.-C.WuW.-H.LiuC.-Y.LaiG.Damage detection of a cable-stayed bridge based on the variation of stay cable forces eliminating environmental temperature effects201617685988010.12989/sss.2016.17.6.859SeoJ.HatfieldG.KimnJ.-H.Probabilistic structural integrity evaluation of a highway steel bridge under unknown trucks201616572StephenG. A.BrownjohnJ. M. W.TaylorC. A.Measurements of static and dynamic displacement from visual monitoring of the Humber Bridge199315319720810.1016/0141-0296(93)90054-82-s2.0-0027334412LawS. S.WardH. S.ShiG. B.ChenR. Z.WaldronP.TaylorC.Dynamic assessment of bridge load-carrying capacities. I1995121347848710.1061/(asce)0733-9445(1995)121:3(478)2-s2.0-0029275933YiJ.-H.ChoS.KooK.-Y.YunC.-B.KimJ.-T.LeeC.-G.LeeW.-T.Structural performance evaluation of a steel-plate girder bridge using ambient acceleration measurements20073328129810.12989/sss.2007.3.3.2812-s2.0-34249697681BurdetO.ZanellaJ. L.Automatic monitoring of bridges using electronic inclinometers16Proceedings of the IABSE Congress Report, 16th Congress of IABSE2000Lucerne, Switzerland1574158110.2749/222137900796314284ParkG.CudneyH. H.InmanD. J.Impedance-based health monitoring of civil structural components20006415316010.1061/(ASCE)1076-0342(2000)6:4(153)2-s2.0-0034529418DuffóG. S.FarinaS. B.Development of an embeddable sensor to monitor the corrosion process of new and existing reinforced concrete structures20092382746275110.1016/j.conbuildmat.2009.04.0012-s2.0-67349261723American Association of State Highway and Transportation Officials (AASHTO)2008Washington, DC, USAAmerican Association of State Highway and Transportation Officials (AASHTO)Ministry of Land, Transport, and Maritime (MLTM) Affairs (2012) Korea highway bridge design code (LRFD), 201219781stOntario Ministry of Transportation and CommunicationsFaulknerB. C.BartonF.BaberT. T.McKeelW. T.Determination of bridge response using acceleration data1996Charlottesville, Va, USAVirgina Transportation Research CouncilSumitoroS.MatsuiY.KonoM.OkamotoT.FujiiK.Long span bridge health monitoring system in Japan4337Health Monitoring and Management of Civil Infrastructure Systems2001Newport Beach, Calif, USAInternational Society for Optics and Photonics517524Proceedings of SPIE10.1117/12.435628MiyataT.YamadaH.KatsuchiH.KitagawaM.Full-scale measurement of Akashi-Kaikyo Bridge during typhoon20029012-151517152710.1016/s0167-6105(02)00267-22-s2.0-0036905359FujinoY.MurataM.OkanoS.TakeguchiM.Monitoring system of the Akashi Kaikyo Bridge and displacement measurement using GPS3995Nondestructive Evaluation of Highways, Utilities, and Pipelines IVJune 2000San Diego, Calif, USAInternational Society for Optics and Photonics229236Proceedings of SPIE10.1117/12.387814Kijewski-CorreaT.KochlyM.Monitoring the wind-induced response of tall buildings: GPS performance and the issue of multipath effects2007959–111176119810.1016/j.jweia.2007.02.0022-s2.0-34548189037JoH.SimS.-H.TatkowskiA.SpencerB. F.Jr.NelsonM. E.Feasibility of displacement monitoring using low-cost GPS receivers20132091240125410.1002/stc.15322-s2.0-84880770300NassifH. H.GindyM.DavisJ.Comparison of laser Doppler vibrometer with contact sensors for monitoring bridge deflection and vibration200538321321810.1016/j.ndteint.2004.06.0122-s2.0-12344308501BaiR.OstachowiczW.RadzienskiM.CaoM.Vibrational damage detection using fractal surface singularities with noncontact laser measurement201622112569258110.1177/1077546314548088OlaszekP.Investigation of the dynamic characteristic of bridge structures using a computer vision method199925322723610.1016/s0263-2241(99)00006-82-s2.0-0033114821YeX. W.DongC. Z.LiuT.Image-based structural dynamic displacement measurement using different multi-object tracking algorithms201617693595610.12989/sss.2016.17.6.935WahbehA. M.CaffreyJ. P.MasriS. F.A vision-based approach for the direct measurement of displacements in vibrating systems200312578579410.1088/0964-1726/12/5/0162-s2.0-0142124230LeeJ. J.ShinozukaM.A vision-based system for remote sensing of bridge displacement200639542543110.1016/j.ndteint.2005.12.0032-s2.0-33748261010ParkK.-T.KimS.-H.ParkH.-S.LeeK.-W.The determination of bridge displacement using measured acceleration200527337137810.1016/j.engstruct.2004.10.0132-s2.0-13444263518LeeH. S.HongY. H.ParkH. W.Design of an FIR filter for the displacement reconstruction using measured acceleration in low-frequency dominant structures201082440343410.1002/nme.27692-s2.0-77949933602ParkJ. W.SimS. H.JungH. J.Development of a wireless displacement measurement system using acceleration responses20131378377839210.3390/s130708377KangL.-H.KimD.-K.HanJ.-H.Estimation of dynamic structural displacements using fiber Bragg grating strain sensors2007305353454210.1016/j.jsv.2007.04.0372-s2.0-34447531830RappS.KangL.-H.HanJ.-H.MuellerU. C.BaierH.Displacement field estimation for a two-dimensional structure using fiber Bragg grating sensors200918202500610.1088/0964-1726/18/2/0250062-s2.0-65349173020ChangS.-J.KimN.-S.Estimation of displacement response from FBG strain sensors using empirical mode decomposition technique201252657358910.1007/s11340-011-9522-z2-s2.0-84861480918ShinS.LeeS.-U.KimY.KimN.-S.Estimation of bridge displacement responses using FBG sensors and theoretical mode shapes201242222924510.12989/sem.2012.42.2.2292-s2.0-84860253307HouX.YangX.HuangQ.Using inclinometers to measure bridge deflection200510556456910.1061/(ASCE)1084-0702(2005)10:5(564)2-s2.0-24944535369SmythA.WuM.Multi-rate Kalman filtering for the data fusion of displacement and acceleration response measurements in dynamic system monitoring200721270672310.1016/j.ymssp.2006.03.0052-s2.0-33750495377SimS. H.SpencerB. F.Jr.NagayamaT.Multimetric sensing for structural damage detection20101371223010.1061/(asce)em.1943-7889.00001992-s2.0-78650228731KimJ.KimK.SohnH.Autonomous dynamic displacement estimation from data fusion of acceleration and intermittent displacement measurements2014421-219420510.1016/j.ymssp.2013.09.0142-s2.0-84887116396SungS. H.ParkJ. W.NagayamaT.JungH. J.A multi-scale sensing and diagnosis system combining accelerometers and gyroscopes for bridge health monitoring201423101500510.1088/0964-1726/23/1/0150052-s2.0-84891382400RobertsG. W.MengX.DodsonA. H.Integrating a global positioning system and accelerometers to monitor the deflection of bridges20041302657210.1061/(ASCE)0733-9453(2004)130:2(65)2-s2.0-2342510243ParkJ.-W.SimS.-H.JungH.-J.Displacement estimation using multimetric data fusion20131861675168210.1109/TMECH.2013.22751872-s2.0-84898599019ParkJ.-W.SimS.-H.JungH.-J.Wireless displacement sensing system for bridges using multi-sensor fusion201423404502210.1088/0964-1726/23/4/0450222-s2.0-84896370088ChoS.SimS.-H.ParkJ.-W.LeeJ.Extension of indirect displacement estimation method using acceleration and strain to various types of beam structures201414469971810.12989/sss.2014.14.4.6992-s2.0-84908086071ChoS.YunC.-B.SimS.-H.Displacement estimation of bridge structures using data fusion of acceleration and strain measurement incorporating finite element model2015153-464566310.12989/sss.2015.15.3.6452-s2.0-84925745900SimS.Estimation of flexibility matrix of beam structures using multisensor fusion2016126064Gebre-EgziabherD.HaywardR. C.PowellJ. D.Design of multi-sensor attitude determination systems200440262764910.1109/TAES.2004.13100102-s2.0-4444319524SchneiderJ.ElingC.KlingbeilL.KuhlmannH.ForstnerW.StachnissC.Fast and effective online pose estimation and mapping for UAVsProceedings of the IEEE International Conference on Robotics and Automation (ICRA '16)May 2016Stockholm, SwedenIEEE4784479110.1109/icra.2016.7487682SunR.HanK.HuJ.WangY.HuM.OchiengW. Y.Integrated solution for anomalous driving detection based on BeiDou/GPS/IMU measurements20166919320710.1016/j.trc.2016.06.006ZhaoY.Performance evaluation of Cubature Kalman filter in a GPS/IMU tightly-coupled navigation system2016119677910.1016/j.sigpro.2015.07.0142-s2.0-84939442945ChenZ.ZouH.JiangH.ZhuQ.SohY. C.XieL.Fusion of WiFi, smartphone sensors and landmarks using the kalman filter for indoor localization201515171573210.3390/s1501007152-s2.0-84920518240ChungH.-Y.HouC.-C.ChenY.-S.Indoor intelligent mobile robot localization using fuzzy compensation and kalman filter to fuse the data of gyroscope and magnetometer201562106436644710.1109/tie.2015.24166922-s2.0-84941584734http://www.sbg-systems.com