Next Article in Journal
Effect of Loading Rate on Tensile and Failure Behavior of Concrete
Previous Article in Journal
High-Temperature Hydrogen Sensing Performance of Ni-Doped TiO2 Prepared by Co-Precipitation Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of IMU Sensor Placement for the Measurement of Lower Limb Joint Kinematics

1
Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20993, USA
2
Division of Clinical Evidence and Analysis 2, Office of Clinical Evidence and Analysis, Office of Product Evaluation and Quality, Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20993, USA
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(21), 5993; https://doi.org/10.3390/s20215993
Submission received: 2 September 2020 / Revised: 25 September 2020 / Accepted: 15 October 2020 / Published: 22 October 2020
(This article belongs to the Section Wearables)

Abstract

:
There is an increased interest in using wearable inertial measurement units (IMUs) in clinical contexts for the diagnosis and rehabilitation of gait pathologies. Despite this interest, there is a lack of research regarding optimal sensor placement when measuring joint kinematics and few studies which examine functionally relevant motions other than straight level walking. The goal of this clinical measurement research study was to investigate how the location of IMU sensors on the lower body impact the accuracy of IMU-based hip, knee, and ankle angular kinematics. IMUs were placed on 11 different locations on the body to measure lower limb joint angles in seven participants performing the timed-up-and-go (TUG) test. Angles were determined using different combinations of IMUs and the TUG was segmented into different functional movements. Mean bias and root mean square error values were computed using generalized estimating equations comparing IMU-derived angles to a reference optical motion capture system. Bias and RMSE values vary with the sensor position. This effect is partially dependent on the functional movement analyzed and the joint angle measured. However, certain combinations of sensors produce lower bias and RMSE more often than others. The data presented here can inform clinicians and researchers of placement of IMUs on the body that will produce lower error when measuring joint kinematics for multiple functionally relevant motions. Optimization of IMU-based kinematic measurements is important because of increased interest in the use of IMUs to inform diagnose and rehabilitation in clinical settings and at home.

1. Introduction

The importance of gait in clinical evaluation is well established. Measurement of human gait has the potential to aid clinicians in making diagnoses, targeting areas for rehabilitation, informing approaches for orthopedic surgery, and tracking rehabilitation progress [1]. While direct observation of human gait is easily done, the quantification of gait is more difficult but important to effectively track changes over time or compare gait to other clinical populations [2]. Optical motion capture (MOCAP) is a common means of quantitively measuring human gait. However, MOCAP is vulnerable to marker occlusion and limited capture volume, confining data collection to a dedicated laboratory space [3]. This prevents the monitoring and evaluation of patients in more realistic environments. There is increased interest in the use of inertial measurement units (IMUs) to measure gait kinematics due to their portability and immunity to MOCAP-specific issues such as occlusion and limited capture volume. As such, there is considerable research on IMU-based gait kinematics to study healthy and pathological gait, with most studies concluding IMU-based joint angle calculations are comparable to MOCAP systems [3,4,5,6,7].
The literature supports the notion that IMUs offer a reasonable alternative to MOCAP when collecting data outside of the laboratory. While this suggests publications related to IMU kinematic measurements in “real-world” contexts might increase, a recent review article indicates most studies assess gait on a treadmill or walkway with few studies examining common functional movements besides straight overground walking [3]. Examining functional movements such as sit-to-stand and turning are important to clinicians because they aid in assessing lower extremity strength and balance [8,9]. If IMU-derived gait kinematics are expected to accurately describe common daily activities, then they need to be proven using more functionally relevant motions. The anatomical placement of sensors is also variable and oftentimes poorly documented [3]. A few publications optimized sensor placement when measuring gait parameters such as stance/swing percent and cadence [10,11] and evaluated sensitivity of pose estimation accuracy to IMU sensor placements during single leg squats [12]. However, to our knowledge, no studies optimize IMU sensor placement for measuring lower limb joint angles during every day movements such as walking, turning, sitting, and standing.
While the utility of IMUs in clinics is recognized, these knowledge gaps may prevent more widespread adoption. Therefore, the aim of this study was to address the two gaps identified above by examining the impact of IMU location on lower limb joint angle accuracy while participants perform multiple functionally relevant movements. Participants performed the Timed-Up-And-Go (TUG) test given its extensive validation in clinical populations and the variety of functional movements performed [13,14,15]. Eleven sensor positions on the torso, thigh, shank, and foot were chosen based on commonly reported locations in the literature. All possible combinations of IMU sensors were used to calculate hip, knee, and ankle angles in the sagittal plane and compared to a reference MOCAP system. The output of this study could inform clinicians and researchers of IMU sensor locations that will produce lower error when measuring joint kinematics for multiple functionally relevant motions.

2. Materials and Methods

2.1. Participants

A convenience sample of seven participants was recruited (4 male/3 female; 26.0 ± 4.0 years of age). The inclusion criteria were as follows: greater than 18 years of age and no self-declared gait impediments or abnormalities Exclusion criteria were individuals younger than 18 years of age or individuals with gait impairments. All participants provided written informed consent prior to participation. The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the U.S. FDA Institutional Review Board (No. 2019-CDRH-002). Participants were asked to wear flat, close-toed shoes. Two-strap Velcro sandals were provided to participants who did not have proper footwear.

2.2. IMU System and Sensor Placements

Xsens MTw Awinda IMUs were used (Xsens, Enschede, The Netherlands). Data capture was linked to an Xsens base station and processed in the MTw Workstation to obtain angular velocities, accelerations, and quaternions. Due to a capture settings error, some data were collected at 40 Hz while the remaining data were captured at 60 Hz. The human 46.1 profile filter was used. A Vicon MOCAP system ran concurrently (Vicon Motion Systems Ltd., Oxford, UK). Regardless of the sampling frequency, a general-purpose output configuration file was created in Vicon to generate an output synchronization signal, similar to a TTL signal (Vicon MX System: Vicon MX Hardware System Reference R1.6). That signal was received by the XSens MTw Awinda receiver station to trigger the start of capture of both systems simultaneously.
A brief literature review informed IMU placements. The body segment and location of sensors placed on the lower body, along with their literature sources, are in Table 1. Recommendations from Xsens tutorials (https://tutorial.xsens.com/) are also included and denoted with an asterisk in the source column. Eleven anatomical locations were selected. Sensors were placed on the right leg (where applicable). Table 2 includes detailed sensor location descriptions. Figure 1 shows sensor placements on a subject.

2.3. MOCAP System and Marker Placement

IMU accuracy was evaluated against a reference MOCAP system as in previous studies [7,16,17,18]. The MOCAP system (Vicon) consisted of eight B10 Bonita and four Vero v1.3 optical cameras. The cameras sampled at 100 Hz or 120 Hz and their positions were optimized to the targeted capture volume. The system was calibrated prior to data collection according to the manufacturer’s specifications. The Vicon Plug-in-Gait (PiG) lower body model was used to analyze movement at the hip, knee, and ankle joints. Fourteen reflective markers were placed on the pelvis and legs of participants. Markers were placed on the anterior superior iliac spines and posterior superior iliac spines of the pelvis, lateral side of the thighs and shanks; flexion–extension axis of the knees, heel, lateral malleolus; and over the second metatarsal head on the mid-foot side of the equinus break. Subject-specific measurements of body mass, height, ankle width, knee width, and leg length were included in the lower body model. Ankle width was defined as the medio-lateral distance across the malleoli; knee width was defined as the medio-lateral width of the knee across the line of the knee axis; and leg length was measured between the anterior superior iliac spine markers and the medial malleolus. Figure 1 shows the front, side, and back views of a participant with the markers in place. Further details on PiG have been previously described and are available on Vicon’s website (https://docs.vicon.com/display/Nexus26/Full+body+modeling+with+Plug-in+Gait).

2.4. Functional Task–TUG Test

Participants performed the Timed-Up-and-Go (TUG) task. The TUG required participants to stand from a chair, walk 3 m, turn around, walk back, and sit down. A standard height chair was used without assistive devices or arm rests. Participants performed three TUG trials at minimum, with two participants performing five. In total, 25 TUG trials were analyzed. Data capture began with subjects standing for IMU calibration purposes.

2.5. IMU Joint Angle Calculations

There are several different approaches to sensor-to-body IMU calibration. No study has compared outcomes of these to date so there is no consensus on the best approach [37]. However, all approaches generally produce kinematic profiles that align well with reference systems such as MOCAP [37]. A calibration similar to the one described by Palermo et al. was used to define calibration vectors describing body segment orientations in terms of sensor coordinate systems and is briefly described here [18]. Participants assumed two poses for static calibration (standing upright and sitting while leaning back with outstretched legs) to determine vectors in the sagittal plane using gravitational acceleration. Cross products between these vectors defined medial–lateral (ML) axis while the standing vector alone defined superior–inferior (SI) axis for each sensor. SI and ML cross products defined anterior–posterior (AP) axis. A final cross product between the AP and SI axes redefined the ML axis to ensure orthogonality. The SI, ML, and AP vectors were normalized. Gravity vectors were averaged over multiple frames during a separate calibration trial for the seated pose. The standing pose vectors were defined during frame one of the TUG data capture. Body segments laid in the sagittal plane during both poses.
Quaternions describing sensor orientation were produced using Xsens software. These quaternions rotated the calibration vectors to determine body segment orientations in terms of sensor coordinate systems for each frame of the TUG. To accomplish this, the calibration vectors were converted into quaternion format (Equation (1)), where Cx, Cy, and Cz are the three components in x, y, and z, and rotated using quaternion conjugation (Equation (2)), where q is sensor quaternion and p is calibration unit vector in quaternion format. The resultant series of quaternions were converted back to unit vectors Equation (1).
[ C x   C y   C z ]   =   [   0   C x   C y   C z ]
q p q 1   =   p
These orientation vectors needed to be transformed to a global coordinate system before angle calculations could occur. Shank and thigh segment axes were assumed parallel during frame one when participants stood upright. This facilitated transformation by defining an instant when the otherwise unrelated sensor coordinate systems were known to be coincident. Direction cosine matrices (DCMs) were constructed between the frame one vectors and a common right-handed coordinate system [1,0,0], [0,1,0], and [0,0,1]. These DCMs rotated unit vectors at each frame to this global coordinate system. DCMs can be constructed as shown in Equation (3) where a, b, and c are the normalized calibration vectors, [x, y, z] is a vector in terms of a sensor coordinate system, and [x′, y′, z′] is the same vector in terms of the global coordinate system.
[ x y z ]   =   [ a 1 b 1 c 1 a 2 b 2 c 2 a 3 b 3 c 3 ] [ x y z ]
Flexion/extension angles were calculated between IMUs placed above and below each joint. DCMs were constructed at each frame to project the SI vector of the segment below a joint into the coordinate system of the segment above a joint. Flexion/extension angles were then calculated using Equation (4) for all three joints and sensor combinations (hip: 8 combinations; knee: 12 combinations; ankle: 6 combinations).
A 1   =   arctan ( C S I , 1 C S I , 3 )
where CSI,1 is the anterior–posterior component of projection and CSI,3 = superior–inferior component of projection.

2.6. MOCAP Joint Angle Calculations

The Vicon PiG model calculated hip flexion/extension, knee flexion/extension, and ankle dorsiflexion/plantarflexion angles. Hip flexion/extension is calculated between the pelvis AP axis, and a projection of the thigh AP axis into the sagittal plane of the pelvis. This plane is perpendicular to an axis passing transversely through the pelvis at the hip joint center. Knee flexion/extension is calculated between the thigh AP axis and the projection of the shank AP axis into the plane perpendicular to the knee flexion axis. Ankle dorsiflexion/plantarflexion is taken between the shank AP axis and the projection of the axis formed by the heel and toe markers into the sagittal plane of the foot. More information on PiG angle calculations is on Vicon’s website (https://docs.vicon.com/display/Nexus26/Full+body+modeling+with+Plug-in+Gait).

2.7. Data Analysis

Each trial of the TUG was segmented manually in the Vicon acquisition software to assess error during specific movements. Since 25 TUG trials were performed by seven subjects, each segment has 25 possible samples to include in the analysis. Walk Pass 1 and Walk Pass 2 were deemed similar and combined during data processing. In the tables presented in the Results Section 3, the row labeled “Walk (1 and 2)” represents the combination of these segments, yielding 50 possible samples for this analysis. The segment definitions are defined below in Table 3:
Error between angles determined by IMU combinations and MOCAP were calculated for all TUG segments for each iteration of TUG. The error was presented as mean bias and root mean square error (RMSE). These metrics were calculated for each segment of every TUG iteration. The mean biases for a single TUG were calculated as an average of the difference between IMU angle measurements and the MOCAP for each TUG iteration. Similarly, the RMSEs were calculated as the square root of the average of squared difference between IMU angle measurements and the MOCAP for each iteration of a given segment of TUG. The mean bias and RMSE for each iteration were compared among the IMU sensor combinations using generalized estimating equations (GEE) with exchangeable correlation to account for the clustering of biases or RMSEs for different TUG iterations within the same participants. The statistical analysis was conducted using SAS 9.4 (SAS Institute Inc., Cary, NC, USA).
Over the course of data collection, data from six sensors were excluded due to faulty recordings across all subjects and trials. Tables showing results of the study indicate the number of samples (n) for a given segment and subjects (subs) used to calculate each metric.

3. Results

Bias and RMSE values between joint angles measured by IMUs and the reference MOCAP system are presented in tabular format (Table 4, Table 5 and Table 6). Rows specify the averaged bias and RMSE over all TUG trials for each TUG segment; columns specify sensor combinations. Cells in the table provide the least squares means of bias and RMSE and their corresponding standard errors, as well as p-values. For bias, the p-value is calculated from a hypothesis test to determine whether the mean bias is different from zero. Since a p-value < 0.05 indicates the bias is significantly different from zero and therefore the measurement is biased, measurements for which p-value > 0.05 are highlighted in grey to emphasize those measurements which are unbiased. For RMSE, the p-value is calculated from the comparison between the RMSE and the lowest RMSE for each combination of sensors of a given joint and TUG segment. A p-value < 0.05 indicates there is statistically significant difference between this RMSE and the lowest RMSE for the given angle and TUG segment. Lower RMSEs indicate better performance; therefore, measurements for which p-value > 0.05 were highlighted in grey to emphasize measurements with values similar to the lowest RMSE.
Table 4 shows the bias and RMSE values for ankle flexion for each segment of the TUG. The lowest absolute bias over all segments was 0.02° (1.11°) (least squares mean (standard error)) obtained during Turn 1 using the LLS/Heel sensor combination. The largest absolute bias of 2.35° (0.54°) was obtained during Turn 2 using the MLS/DFoot sensor combination. Over all segments, the lowest RMSE values of 2.33° (0.57°) and 2.58° (0.61°) were obtained during the Sit-to-Stand and Stand-to-Sit movements, respectively, using the MLS/DFoot IMU sensor combination. The highest RMSE value for ankle flexion was 5.49° (0.91°) during the Turn 2 segment with the LLS/Heel sensor combination.
Table 5 shows the bias and RMSE values for knee flexion for each segment of the TUG. The lowest absolute bias over all segments was 0.08° ± 1.01° obtained during Turn 1 using the MLS/LAT sensor combination. The largest absolute bias values over all segments was −7.06° (1.52°) and −7.00° (2.08°) obtained using the Shin/LLT sensor combination during Sit-to-Stand and Stand-to-Sit movements, respectively. The Shin/LLT sensor combination also generated the lowest and highest RMSE values of 4.62° (0.63°) during Turn 2 and 9.31° (1.59°) during Sit-to-Stand, respectively.
For hip flexion, the lowest absolute bias over all segments was 0.13° (2.24°) obtained during Turn 2 using the Sacrum/MLT sensor combination; the largest absolute bias was 16.30° (6.43°) obtained during Stand-to-Sit using the L4-L5/LPT sensor combination (Table 6). The lowest RMSE value of 4.35° (0.64°) was obtained during Walk using the L4-L5/LLT sensor combination while the highest RMSE value of 21.46° (4.52°) was obtained during Stand-to-Sit using the L4-L5/LPT sensor combination.

4. Discussion

This study evaluated the impact of IMU sensor location on lower limb joint angle accuracy and bias during the TUG. A brief literature review pointed to common sensor locations that were incorporated into an experimental protocol comparing joint angles derived from combinations of IMU sensors to a reference MOCAP system. The results have several implications for research in motion analysis and clinical implementation of IMU-based joint kinematics.
There are several different approaches to calibration of the sensor coordinate system and calculating joint angles based on IMU sensor data. In comparing RMSE values to previous studies investigating joint angles measured with IMUs and MOCAP systems, our results are similar and even surpass the accuracy reported for other approaches, yielding confidence in the validity of our kinematic calculation. The average level-walking RMSE values across all sensor combinations in this study were 3.90°, 6.35°, and 5.93°, for the ankle, knee, and hip, respectively. These values are close to, and even surpass, the ankle, knee, and hip RMSE values presented by Tadano et al. (9.75°, 7.88°, and 10.14°, respectively) and Dorschky et al. (4.60°, 5.30°, and 8.70°, respectively) [30,38]. Other studies have reported lower RMSE values at the same joints [39,40]. A comparison of kinematic calculation approaches was not within the scope of this paper, but sufficient detail on the methodology was provided such that this approach can be replicated and applied.
Many studies have investigated the effectiveness of using IMUs for diagnosing and monitoring diseases in clinics and at home, where cost and space requirements make MOCAP infeasible [12,14,22,23,25,26,27,28]. The interest in using IMUs for clinical diagnosis and rehabilitation makes the question of their accuracy paramount. Optimization of sensor placement for measuring kinematics is one step towards making this technology a reliable tool for clinicians. Our experimental and analysis approaches allow for the determination of the best sensor placement based on task being performed and the joint being evaluated. For example, if the Sit-to-Stand movement were being investigated across the ankle, knee, and hip joints, the results presented in Table 4, Table 5 and Table 6 suggest the DFoot, MLS, LAT, and Sacrum sensors would produce unbiased, accurate results across all joints. Similarly, these data can be used to inform sensor placement for a specific joint. For example, unbiased and accurate measurements for knee flexion can be derived using the MLS/LAT, MLS/MLT, or LLS/LAT sensor combinations (Table 5). Sensor combinations that are highly biased and inaccurate were also identified in this study. The L4-L5/LPT sensor combination for hip flexion produced some of the highest bias and RMSE values across all functional movements of the TUG test (Table 6).
Determining an overall best configuration of sensors across all functional movements and joints would be useful. Based on the data, a configuration of sensors including the sacrum, lower anterior thigh, lower lateral shank, and heel locations will generally perform well. However, this may not be the best approach in all cases. While the current study expands on previous work by investigating multiple functional movements of interest, recommendations are based on data from healthy individuals and may have limited applicability to clinical populations with impaired gait. Further research in determining ideal sensor locations for clinical populations is needed. The analysis in this study was also limited to sagittal plane kinematics to enable comparisons between previous studies. Although the sagittal plane kinematics are more commonly reported in IMU gait studies, meaningful clinical information can be derived from frontal and transverse plane kinematics. Additional analysis is needed to assess accuracy in other planes of motion to determine the optimal sensor placements for all DOFs during specific functional activities. We also acknowledge that our sample size was small, although significant differences in the sensor combinations were still identified. To increase confidence in the generalizability of these results to a healthy, unimpaired population, additional participants should be included in the analysis.
The IMU coordinate systems are also slightly offset from those used in the PiG model. IMUs inherently cannot measure absolute location relative to anatomical landmarks meaning there are inherent differences in how both systems align segments [9]. The assumption is made that the axes of the sensors are parallel when the participant is standing still and that the joint angles are zero. This provides a good estimate of segment orientation but does not perfectly align axes with Vicon which uses its own calibration procedure. The direct attachment of reflective markers to the IMU sensors would allow for direct alignment of the IMU and Vicon coordinate systems and a more technical comparison of the kinematic models, but there is value in understanding the IMU kinematic output independent of Vicon since this more accurately reflects how the IMUs will be used in real-world scenarios.
In conclusion, this study compared accuracy and bias of lower limb joint kinematics derived from various IMU sensor locations to a reference MOCAP system. The findings can be used to inform the wearable sensors community of anatomical locations that are less prone to error when measuring joint angles for specific lower limb tasks. The findings here are of interest in clinical research, diagnosis, and rehabilitation by making IMU technology more reliable for measuring gait kinematics. It can be suggested that the sacrum, lower anterior thigh, middle lateral shank, and heel sensors will produce relatively low error for the three joints as compared to the reference system. However, the combination of sensor locations used should ultimately be driven by the motions and/or joints of interest.
The mention of commercial products, their sources, or their use in connection with material reported herein is not to be construed as either an actual or implied endorsement of such products by the Department of Health and Human Services.
The data and code used for the current analysis are available online at: https://github.com/dbp-osel/IMU-Sensor-Placement-Optimization.

Author Contributions

Conceptualization (K.K.); Methodology (K.K., W.N.); Software (W.N., W.W.); Validation (W.N., W.W.); Formal Analysis (W.N., W.W.); Investigation (K.K., W.N.); Resources (K.K.); Data curation (K.K., W.N.); Writing—original draft (W.N.); Writing—reviewing and editing (K.K., W.N., W.W.); Visualization (K.K., W.N.); Supervision (K.K.); Funding acquisition (K.K.). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the FDA Critical Path Initiative, and base funding from the Division of Biomedical Physics and Division of Applied Mechanics (FDA). The research was also supported in part by an appointment to the Research Participation Program at the U.S. FDA administered by the Oak Ridge Institute for Science and Education (ORISE) through an interagency agreement between the U.S. Department of Energy and FDA. The funders had no role in data collection and analysis, decision to publish, or preparation of the manuscript.

Acknowledgments

The authors would like to thank Stephen Cain from the University of Michigan and Melissa Goodman from the Mayo Clinic for their valuable insights regarding the use and calibration of body worn IMUs.

Conflicts of Interest

The authors have no conflict to disclose.

References

  1. Hodgins, D. The importance of measuring human gait. Med. Device Technol. 2008, 19, 44–47. [Google Scholar]
  2. Muro-de-la-Herran, A.; Garcia-Zapirain, B.; Mendez-Zorrilla, A. Gait analysis methods: An overview of wearable and non-wearable systems, highlighting clinical applications. Sensors 2014, 14, 3362–3394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Weygers, I.; Kok, M.; Konings, M.; Hallez, H.; De Vroey, H.; Claeys, K. Inertial sensor-based lower limb joint kinematics: A methodological systematic review. Sensors 2020, 20, 673. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Reinfelder, S.; Hauer, R.; Barth, J.; Klucken, J.; Eskofier, B.M. Timed Up-and-Go phase segmentation in Parkinson’s disease patients using unobtrusive inertial sensors. In Proceedings of the 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Milan, Italy, 25–29 August 2015; pp. 5171–5174. [Google Scholar]
  5. Maqbool, H.F.; Husman, M.A.B.; Awad, M.I.; Abouhossein, A.; Mehryar, P.; Iqbal, N.; Dehghani-Sanij, A.A. Real-time gait event detection for lower limb amputees using a single wearable sensor. In Proceedings of the 2016 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Orlando, FL, USA, 16–20 August 2016; pp. 5067–5070. [Google Scholar]
  6. Picerno, P. 25 years of lower limb joint kinematics by using inertial and magnetic sensors: A review of methodological approaches. Gait Posture 2017, 51, 239–246. [Google Scholar] [CrossRef]
  7. Bolink, S.; Naisas, H.; Senden, R.; Essers, H.; Heyligers, I.; Meijer, K.; Grimm, B. Validity of an inertial measurement unit to assess pelvic orientation angles during gait, sit–stand transfers and step-up transfers: Comparison with an optoelectronic motion capture system. Med. Eng. Phys. 2016, 38, 225–231. [Google Scholar] [CrossRef]
  8. Janssen, W.G.; Bussmann, H.B.; Stam, H.J. Determinants of the sit-to-stand movement: A review. Phys. Ther. 2002, 82, 866–879. [Google Scholar] [CrossRef]
  9. Adusumilli, G.; Lancia, S.; Levasseur, V.A.; Amblee, V.; Orchard, M.; Wagner, J.M.; Naismith, R.T. Turning is an important marker of balance confidence and walking limitation in persons with multiple sclerosis. PLoS ONE 2018, 13, e0198178. [Google Scholar] [CrossRef]
  10. Panebianco, G.P.; Bisi, M.C.; Stagni, R.; Fantozzi, S. Analysis of the performance of 17 algorithms from a systematic review: Influence of sensor position, analysed variable and computational approach in gait timing estimation from IMU measurements. Gait Posture 2018, 66, 76–82. [Google Scholar] [CrossRef]
  11. Anwary, A.R.; Yu, H.; Vassallo, M. Optimal foot location for placing wearable IMU sensors and automatic feature extraction for gait analysis. IEEE Sens. J. 2018, 18, 2555–2567. [Google Scholar] [CrossRef]
  12. Kianifar, R.; Joukov, V.; Lee, A.; Raina, S.; Kulić, D. Inertial measurement unit-based pose estimation: Analyzing and reducing sensitivity to sensor placement and body measures. J. Rehabil. Assist. Technol. Eng. 2019, 6. [Google Scholar] [CrossRef]
  13. Hafsteinsdóttir, T.B.; Rensink, M.; Schuurmans, M. Clinimetric properties of the timed up and go test for patients with stroke: A systematic review. Top. Stroke Rehabil. 2014, 21, 197–210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Mollinedo, I.; Ma Cancela, J. Evaluation of the psychometric properties and clinical applications of the Timed Up and Go test in Parkinson disease: A systematic review. J. Exerc. Rehabil. 2020, 16, 302–312. [Google Scholar] [CrossRef]
  15. Long, J.; Cai, T.; Huang, X.; Zhou, Y.; Kuang, J.; Wu, L. Reference value for the TUGT in healthy older people: A systematic review and meta-analysis. Geriatr. Nurs. 2020, 41, 325–330. [Google Scholar] [CrossRef] [PubMed]
  16. Kong, W.; Sessa, S.; Cosentino, S.; Zecca, M.; Saito, K.; Wang, C.; Imtiaz, U.; Lin, Z.; Bartolomeo, L.; Ishii, H. Development of a real-time IMU-based motion capture system for gait rehabilitation. In Proceedings of the 2013 IEEE International Conference on Robotics and Biomimetics (ROBIO), Shenzhen, China, 12–14 December 2013; pp. 2100–2105. [Google Scholar]
  17. Kwakkel, S.; Godha, S.; Lachapelle, G. Foot and Ankle Kinematics during Gait Using Foot Mounted Inertial System; ION NTM: San Diego, CA, USA, 2007. [Google Scholar]
  18. Palermo, E.; Rossi, S.; Marini, F.; Patanè, F.; Cappa, P. Experimental evaluation of accuracy and repeatability of a novel body-to-sensor calibration procedure for inertial sensor-based gait analysis. Measurement 2014, 52, 145–155. [Google Scholar] [CrossRef]
  19. Laudanski, A.; Brouwer, B.; Li, Q. Measurement of lower limb joint kinematics using inertial sensors during stair ascent and descent in healthy older adults and stroke survivors. J. Healthc. Eng. 2013, 4, 555–576. [Google Scholar] [CrossRef] [Green Version]
  20. Barrois, R.; Gregory, T.; Oudre, L.; Moreau, T.; Truong, C.; Pulini, A.A.; Vienne, A.; Labourdette, C.; Vayatis, N.; Buffat, S. An automated recording method in clinical consultation to rate the limp in lower limb osteoarthritis. PLoS ONE 2016, 11, e0164975. [Google Scholar] [CrossRef]
  21. Spain, R.; George, R.S.; Salarian, A.; Mancini, M.; Wagner, J.; Horak, F.; Bourdette, D. Body-worn motion sensors detect balance and gait deficits in people with multiple sclerosis who have normal walking speed. Gait Posture 2012, 35, 573–578. [Google Scholar] [CrossRef] [Green Version]
  22. Mancini, M.; Chiari, L.; Holmstrom, L.; Salarian, A.; Horak, F.B. Validity and reliability of an IMU-based method to detect APAs prior to gait initiation. Gait Posture 2016, 43, 125–131. [Google Scholar] [CrossRef] [Green Version]
  23. Esser, P.; Dawes, H.; Collett, J.; Feltham, M.G.; Howells, K. Assessment of spatio-temporal gait parameters using inertial measurement units in neurological populations. Gait Posture 2011, 34, 558–560. [Google Scholar] [CrossRef]
  24. Esser, P.; Dawes, H.; Collett, J.; Howells, K. IMU: Inertial sensing of vertical CoM movement. J. Biomech. 2009, 42, 1578–1581. [Google Scholar] [CrossRef]
  25. Doheny, E.P.; McGrath, D.; Greene, B.R.; Walsh, L.; McKeown, D.; Cunningham, C.; Crosby, L.; Kenny, R.A.; Caulfield, B. Displacement of centre of mass during quiet standing assessed using accelerometry in older fallers and non-fallers. In Proceedings of the 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, San Diego, CA, USA, 28 August–1 September 2012; pp. 3300–3303. [Google Scholar]
  26. Vargas-Valencia, L.S.; Elias, A.; Rocon, E.; Bastos-Filho, T.; Frizera, A. An IMU-to-body alignment method applied to human gait analysis. Sensors 2016, 16, 2090. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Bourgeois, A.B.; Mariani, B.; Aminian, K.; Zambelli, P.; Newman, C. Spatio-temporal gait analysis in children with cerebral palsy using, foot-worn inertial sensors. Gait Posture 2014, 39, 436–442. [Google Scholar] [CrossRef] [PubMed]
  28. Guo, Y.; Wu, D.; Liu, G.; Zhao, G.; Huang, B.; Wang, L. A low-cost body inertial-sensing network for practical gait discrimination of hemiplegia patients. Telemed. e-Health 2012, 18, 748–754. [Google Scholar] [CrossRef] [PubMed]
  29. Hsu, Y.-L.; Chung, P.-C.; Wang, W.-H.; Pai, M.-C.; Wang, C.-Y.; Lin, C.-W.; Wu, H.-L.; Wang, J.-S. Gait and balance analysis for patients with Alzheimer’s disease using an inertial-sensor-based wearable instrument. IEEE J. Biomed. Health Inform. 2014, 18, 1822–1830. [Google Scholar] [CrossRef]
  30. Tadano, S.; Takeda, R.; Miyagawa, H. Three dimensional gait analysis using wearable acceleration and gyro sensors based on quaternion calculations. Sensors 2013, 13, 9321–9343. [Google Scholar] [CrossRef]
  31. Scapellato, S.; Cavallo, F.; Martelloni, C.; Sabatini, A.M. In-use calibration of body-mounted gyroscopes for applications in gait analysis. Sens. Actuators A Phys. 2005, 123, 418–422. [Google Scholar] [CrossRef]
  32. Khan, A.; Biddiss, E. Musical Stairs: A motivational therapy tool for children with disabilities featuring automated detection of stair-climbing gait events via inertial sensors. Med. Eng. Phys. 2017, 40, 95–102. [Google Scholar] [CrossRef]
  33. Lau, H.; Tong, K. The reliability of using accelerometer and gyroscope for gait event identification on persons with dropped foot. Gait Posture 2008, 27, 248–257. [Google Scholar] [CrossRef]
  34. Rebula, J.R.; Ojeda, L.V.; Adamczyk, P.G.; Kuo, A.D. Measurement of foot placement and its variability with inertial sensors. Gait Posture 2013, 38, 974–980. [Google Scholar] [CrossRef] [Green Version]
  35. Rampp, A.; Barth, J.; Schülein, S.; Gaßmann, K.-G.; Klucken, J.; Eskofier, B.M. Inertial sensor-based stride parameter calculation from gait sequences in geriatric patients. IEEE Trans. Biomed. Eng. 2014, 62, 1089–1097. [Google Scholar] [CrossRef]
  36. Sijobert, B.; Denys, J.; Coste, C.A.; Geny, C. IMU based detection of freezing of gait and festination in parkinson’s disease. In Proceedings of the 2014 IEEE 19th International Functional Electrical Stimulation Society Annual Conference (IFESS), Kuala Lumpur, Malaysia, 7–19 September 2014. [Google Scholar]
  37. Pacher, L.; Chatellier, C.; Vauzelle, R.; Fradet, L. Sensor-to-segment calibration methodologies for lower-body kinematic analysis with inertial sensors: A systematic review. Sensors 2020, 20, 3322. [Google Scholar] [CrossRef] [PubMed]
  38. Dorschky, E.; Nitschke, M.; Seifer, A.-K.; van den Bogert, A.J.; Eskofier, B.M. Estimation of gait kinematics and kinetics from inertial sensor data using optimal control of musculoskeletal models. J. Biomech. 2019, 95, 109278. [Google Scholar] [CrossRef] [PubMed]
  39. Takeda, R.; Tadano, S.; Natorigawa, A.; Todoh, M.; Yoshinari, S. Gait posture estimation using wearable acceleration and gyro sensors. J. Biomech. 2009, 42, 2486–2494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Ohtaki, Y.; Sagawa, K.; Inooka, H. A method for gait analysis in a daily living environment by body-mounted instruments. JSME Int. J. Ser. C Mech. Syst. Mach. Elem. Manuf. 2001, 44, 1125–1132. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Participant set-up. Depiction of the location of each IMU sensor (orange) and Plug-In-Gait body model marker locations for: (A) the front view; (B) the side view; and (C) the back view. Orientations of the superior–inferior (SI), medial–lateral (ML), and anterior–posterior (AP) axes are also shown in (A).
Figure 1. Participant set-up. Depiction of the location of each IMU sensor (orange) and Plug-In-Gait body model marker locations for: (A) the front view; (B) the side view; and (C) the back view. Orientations of the superior–inferior (SI), medial–lateral (ML), and anterior–posterior (AP) axes are also shown in (A).
Sensors 20 05993 g001
Table 1. Description of IMU sensors placements found in the literature and their sources. Recommendations from Xsens tutorials (https://tutorial.xsens.com/) are also included and denoted with a * in the source column.
Table 1. Description of IMU sensors placements found in the literature and their sources. Recommendations from Xsens tutorials (https://tutorial.xsens.com/) are also included and denoted with a * in the source column.
Body SegmentLocationNumber of SourcesSources
PelvisL4-L58Laudanski 2013 [19], Panebianco 2018 [10], Barrois 2016 [20], Spain 2012 [21], Mancini 2016 [22], Esser 2011 [23], Esser 2009 [24], Doheny 2012 [25]
Sacrum2 *Vargas-Valencia 2016 [26]
FootDorsal foot13 *Laudanski 2013 [19], Panebianco 2018 [10], Barrois 2016 [20], Vargas-Valencia 2016 [26], Bourgeois 2014 [27], Guo 2012 [28], Hsu 2014 [29], Tadano 2013 [30], Kong 2013 [16], Scapellato 2005 [31], Kwakkel 2007 [17], Anwary 2018 [11]
Heel5Kwakkel 2007 [17], Anwary 2018 [11], Khan 2017 [32], Lau 2008 [33], Rebula 2013 [34]
Lateral, below lateral malleolus3Anwary 2018 [11], Rampp 2014 [35], Reinfelder 2015 [4]
ShankLateral mid-shank2Laudanski 2013 [19], Kong 2013 [16]
Flat surface of shin bone1 *
Lateral, just above lateral malleolus4Panebianco 2018 [10], Vargas-Valencia 2016, Guo 2012 [28], Sijobert 2014 [36]
Anterior4Spain 2012 [21], Tadano 2013 [30], Kwakkel 2007 [17], Maqbool 2016 [5]
Tibial tuberosity1Lau 2008 [33]
ThighLateral mid-thigh3 *Laudanski 2013 [19], Kong 2013 [16]
Lateral near knee2Vargas-Valencia 2016 [26], Guo 2012 [28]
Anterior near knee2Tadano 2013 [30], Lau 2008 [33]
Table 2. Description of the IMU sensor locations on the lower body.
Table 2. Description of the IMU sensor locations on the lower body.
TorsoL4-L5L4/L5 lumbar spine
Sacrumon the sacrum
ThighLATLower Anterior Thigh: anterior thigh, 5 cm above the knee joint axis.
MLTMiddle Lateral Thigh: lateral thigh, halfway between the hip and knee joints.
LLTLower Lateral Thigh: lateral thigh, 5 cm above the knee joint axis.
LPTLower Posterior Thigh: posterior thigh, 5 cm above the knee joint axis.
ShankShinShin Bone: hard surface of tibial bone, below the knee and above the thickest part of the calf.
MLSMiddle Lateral Shank: lateral shank, halfway between the knee and ankle.
LLSLower Lateral Shank: lateral shank, 5 cm above the lateral malleolus.
FootHeelAdhered to heel on the back of participant’s shoe
DFootDorsal Foot: under the tongue of the participant’s shoe, approximately over the distal end of the third and fourth metatarsal bones.
Table 3. Description of TUG segmentation.
Table 3. Description of TUG segmentation.
Sit-to-StandStarts when participant begins to lean forward; ends with first heel strike
Walk Pass 1Starts with first heel strike; ends with final toe off prior to participant turning
Turn 1Starts with final toe off prior to starting turn; ends with first heel strike out of turn
Walk Pass 2Starts with first heel strike out of turn; ends with final toe off prior to participant turning
Turn 2Starts with final toe off prior to starting turn; ends when participant begins to bend knees to sit
Stand-to-SitStarts when participant begins to bend knees to sit; ends when participant is seated upright
Table 4. Ankle flexion bias and RMSE [least squares mean (standard error)] for each sensor combination versus gold standard MOCAP system.
Table 4. Ankle flexion bias and RMSE [least squares mean (standard error)] for each sensor combination versus gold standard MOCAP system.
Sit to Stand ShinMLSLLS
Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)
Bias−1.98 (0.81)−1.00 (0.93)−1.44 (0.62)−0.40 (0.85)−1.30 (0.64)0.09 (0.96)
p value *0.0150.2800.0210.6390.0420.927
RMSE3.18 (0.66)2.86 (0.50)2.41 (0.48)2.33 (0.57)2.72 (0.25)2.69 (0.58)
p value *0.3320.2710.866n/a0.3510.265
Stand to Sit ShinMLSLLS
Heel (n = 16 from 5 subs)DFoot (n = 17 from 6 subs)Heel (n = 16 from 5 subs)DFoot (n = 17 from 6 subs)Heel (n = 16 from 5 subs)DFoot (n = 17 from 6 subs)
Bias−1.53 (1.19)−1.87 (1.25)−0.76 (1.13)−0.83 (1.03)−0.56 (1.41)−0.43 (1.19)
p value *0.2000.1350.4980.4180.6910.717
RMSE3.80 (0.49)3.73 (0.68)2.71 (0.71)2.58 (0.61)3.31 (0.49)2.89 (0.57)
p value #0.0260.1140.538n/a0.0090.444
Turn 1 ShinMLSLLS
Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)
Bias−0.84 (0.60)−1.10 (1.09)−1.68 (0.39)−1.96 (0.63)0.02 (1.11)−1.58 (0.66)
p value *0.1670.316< 0.0010.0020.9890.017
RMSE4.51 (0.66)4.56 (0.66)3.94 (0.45)4.13 (0.68)4.54 (0.71)4.21 (0.53)
p value *0.1240.070n/a0.7620.3110.647
Turn 2 ShinMLSLLS
Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)Heel (n = 19 from 5 subs)DFoot (n = 20 from 6 subs)
Bias0.79 (0.88)−1.56 (0.73)−1.09 (0.56)−2.35 (0.54)−1.49 (1.48)−1.76 (0.75)
p value *0.3730.0330.052< 0.0010.3140.020
RMSE3.93 (0.68)4.71 (0.70)3.85 (0.58)4.58 (0.77)5.49 (0.91)4.74 (0.59)
p value #0.6260.016n/a0.0950.0770.030
Walk (1 and 2) ShinMLSLLS
Heel (n = 38 from 5 subs)DFoot (n = 40 from 6 subs)Heel (n = 38 from 5 subs)DFoot (n = 40 from 6 subs)Heel (n = 38 from 5 subs)DFoot (n = 40 from 6 subs)
Bias0.97 (0.84)−0.39 (0.90)−0.16 (0.65)−1.10 (0.53)1.03(0.73)−0.05 (0.52)
p value *0.2460.6620.8010.0380.1570.927
RMSE4.10 (0.56)4.44 (0.46)3.40 (0.32)3.90 (0.32)3.95 (0.51)3.59 (0.31)
p value #0.0350.005n/a0.2380.1180.532
* Bias p-values were calculated from hypothesis test to determine whether the mean bias is different from 0. # RMSE p-values were calculated from the comparison between each RMSE and the least RMSE among different combinations. The number of samples (n) and number of subjects (subs) included in each analysis are indicated for each sensor combination.
Table 5. Knee flexion bias and RMSE [least squares mean (standard error)] for each sensor combination versus gold standard MOCAP system.
Table 5. Knee flexion bias and RMSE [least squares mean (standard error)] for each sensor combination versus gold standard MOCAP system.
Sit to Stand ShinMLSLLS
LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)
Bias−2.49 (1.21)−3.19 (1.07)−7.06 (1.52)−0.87 (1.53)−2.14 (1.36)−2.75 (1.17)−6.58 (1.62)−0.51 (1.68)−1.91 (1.21)−2.62 (1.39)−6.46 (1.68)−0.24 (1.58)
p value *0.0390.003<0.0010.5690.1140.019<0.0010.7630.1140.060<0.0010.881
RMSE5.51 (0.86)5.58 (1.04)9.31 (1.59)5.05 (1.29)5.23 (0.93)5.23 (1.08)8.64 (1.66)4.93 (1.42)4.82 (0.79)5.37 (1.15)8.32 (1.67)4.93 (1.23)
p value #0.0130.307<0.0010.8170.1340.5580.0010.923n/a0.4860.0030.911
Stand to Sit ShinMLSLLS
LAT (n = 21 from 7 subs)MLT (n = 21 from 7 subs)LLT (n = 21 from 7 subs)LPT (n = 21 from 7 subs)LAT (n = 21 from 7 subs)MLT (n = 21 from 7 subs)LLT (n = 21 from 7 subs)LPT (n = 21 from 7 subs)LAT (n = 21 from 7 subs)MLT (n = 21 from 7 subs)LLT (n = 21 from 7 subs)LPT (n = 21 from 7 subs)
Bias−1.27 (1.57)−1.34 (1.41)−7.00 (2.08)1.33 (1.41)−0.76 (1.70)−0.63 (1.61)−6.27 (2.23)1.89 (1.63)−0.75 (1.68)−0.69 (1.89)−6.35 (2.36)1.91 (1.62)
p value *0.4170.3430.0010.3440.6550.6950.0050.2460.6560.7150.0070.240
RMSE6.16 (0.54)5.13 (0.78)8.82 (2.07)5.67 (0.88)5.79 (0.69)4.74 (0.93)7.78 (2.2)5.75 (1.06)5.34 (0.72)4.78 (1.23)7.27 (2.49)5.56 (0.99)
p value #0.0720.2620.0100.5120.148n/a0.0470.5380.1780.9350.1400.621
Turn 1 ShinMLSLLS
LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)
Bias1.19 (1.06)1.62 (1.25)−0.53 (1.22)3.44 (0.76)0.08 (1.01)0.61 (1.12)−1.59 (1.15)2.3 (0.95)1.19 (0.67)1.76 (0.70)−0.49 (0.50)3.42 (0.78)
p value *0.2650.1950.664<0.0010.9340.5880.1680.0150.0760.0120.326<0.001
RMSE6.26 (0.56)6.81 (0.79)6.31 (0.86)6.76 (0.71)6.23 (0.60)6.64 (0.66)6.29 (0.90)6.32 (0.68)6.37 (0.64)6.81 (0.72)6.13 (0.77)6.94 (0.69)
p value #0.7650.3270.4890.4050.8140.4280.4930.8400.6110.292n/a0.217
Turn 2 ShinMLSLLS
LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)
Bias1.77 (0.71)2.61 (0.93)0.10 (0.56)4.12 (1.05)0.11 (0.92)1.03 (0.94)−1.51 (0.90)2.41 (0.88)−0.52 (1.62)0.5 (1.46)−2.2 (1.63)1.76 (1.29)
p value *0.0130.0050.858<0.0010.9040.2720.0940.0060.7490.7300.1780.173
RMSE5.4 (0.54)4.84 (1.01)4.62 (0.63)5.68 (0.88)5.48 (0.58)4.63 (0.82)4.96 (0.92)4.98 (0.68)6.37(0.89)6.01 (1.17)5.96 (1.33)5.93 (0.86)
p value #0.1500.756n/a0.0140.0190.9790.4810.506<0.0010.0240.1530.013
Walk (1 and 2) ShinMLSLLS
LAT (n = 50 from 7 subs)MLT (n = 50 from 7 subs)LLT (n = 50 from 7 subs)LPT (n = 50 from 7 subs)LAT (n = 50 from 7 subs)MLT (n = 50 from 7 subs)LLT (n = 50 from 7 subs)LPT (n = 50 from 7 subs)LAT (n = 50 from 7 subs)MLT (n = 50 from 7 subs)LLT (n = 50 from 7 subs)LPT (n = 50 from 7 subs)
Bias−0.45 (0.79) 0.5731.65 (0.76)−0.80 (0.67) 0.2342.76 (0.5)−1.43 (1.07) 0.1810.71 (0.94) 0.451−1.75 (0.97) 0.0711.75 (0.94) 0.062−0.41(1.05) 0.6971.78 (0.73)−0.79 (0.89) 0.3762.73 (0.86)
p value *0.030<0.0010.0150.002
RMSE6.02 (0.34)6.48 (0.37)6.02 (0.56)6.05 (0.33)6.37 (0.50)6.53 (0.44)6.30 (0.66)5.83 (0.41)6.43 (0.41)7.24 (0.41)6.38 (0.59)6.51 (0.42)
p value #0.7170.2260.7850.2760.4280.2300.566n/a0.3190.0140.4680.001
* Bias p-values were calculated from hypothesis test to determine whether the mean bias is different from 0. # RMSE p-values were calculated from the comparison between each RMSE and the least RMSE among different combinations. The number of samples (n) and number of subjects (subs) included in each analysis are indicated for each sensor combination.
Table 6. Hip flexion bias and RMSE [least squares mean (standard error)] for each sensor combination versus gold standard MOCAP system.
Table 6. Hip flexion bias and RMSE [least squares mean (standard error)] for each sensor combination versus gold standard MOCAP system.
Sit to Stand SacrumL4–L5
LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)
Bias0.41 (2.78)−0.33 (2.79)−4.11 (2.97)2.13 (2.85)9.50 (5.00)8.76 (4.27)4.98 (4.50)11.21 (5.31)
p value *0.8820.9070.1670.4550.0580.0400.2680.035
RMSE6.95 (1.53)6.76 (1.49)7.03 (1.41)8.44 (1.39)15.59 (3.70)14.39 (3.13)12.69 (2.57)17.89 (3.76)
p value #0.759n/a0.8390.1240.0450.0430.0670.014
Stand to Sit SacrumL4–L5
LAT (n = 21 from 7 subs)MLT (n = 21 from 7 subs)LLT (n = 21 from 7 subs)LPT (n = 21 from 7 subs)LAT (n = 21 from 7 subs)MLT (n = 21 from 7 subs)LLT (n = 21 from 7 subs)LPT (n = 21 from 7 subs)
Bias1.55 (3.48)1.37 (3.50)−4.23 (3.92)4.33 (3.29)13.53 (5.75)13.39 (5.14)7.75 (5.17)16.30 (6.43)
p value *0.6570.6950.2800.1890.0190.0090.1340.011
RMSE7.18 (2.01)6.95 (2.04)7.32 (1.40)9.31 (1.53)18.28 (4.06)17.31 (3.71)14.41 (2.80)21.46 (4.52)
p value #0.790n/a0.8630.0360.0260.0250.0540.012
Turn 1 SacrumL4–L5
LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)
Bias−1.12 (1.65)−0.38 (2.20)−2.34 (2.06)1.38 (1.91)2.16 (1.75)2.84 (1.62)0.91 (1.55)4.65 (2.17)
p value *0.4970.8620.2550.4710.2190.0810.5580.032
RMSE5.36 (0.58)6.30 (0.95)5.34 (1.18)6.06 (0.74)6.22 (0.67)6.48 (0.77)5.20 (0.63)7.93 (1.11)
p value #0.8760.3840.9170.4890.1200.010n/a0.002
Turn 2 SacrumL4–L5
LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)LAT (n = 25 from 7 subs)MLT (n = 25 from 7 subs)LLT (n = 25 from 7 subs)LPT (n = 25 from 7 subs)
Bias−1.00 (1.50)0.13 (2.24)−2.31 (1.82)1.57 (1.93)2.93 (1.94)4.04 (2.14)1.59 (1.93)5.50 (2.71)
p value *0.5050.9540.2040.4160.1310.0590.4100.042
RMSE4.47 (0.64)5.21 (1.44)4.74 (0.87)5.34 (0.95)6.95 (0.99)7.14 (1.28)5.91 (1.04)8.75 (1.68)
p value #n/a0.3990.6150.0950.0780.0700.3160.039
Walk (1 and 2) SacrumL4–L5
LAT (n = 50 from 7 subs)MLT (n = 50 from 7 subs)LLT (n = 50 from 7 subs)LPT (n = 50 from 7 subs)LAT (n = 50 from 7 subs)MLT (n = 50 from 7 subs)LLT (n = 50 from 7 subs)LPT (n = 50 from 7 subs)
Bias−3.08 (1.65)−0.89 (1.97)−3.21 (1.85)0.14 (1.83)−0.29 (1.57)1.86 (1.53)−0.47 (1.45)2.92 (2.01)
p value *0.0620.6520.0820.9390.8520.2250.7470.147
RMSE5.95 (0.90)6.65 (0.80)5.37 (1.08)5.78 (0.87)5.74 (0.64)6.46 (0.64)4.35 (0.64)7.10 (0.76)
p value #0.1710.0250.4630.201<0.0010.002n/a0.003
* Bias p-values were calculated from hypothesis test to determine whether the mean bias is different from 0. # RMSE p-values were calculated from the comparison between each RMSE and the least RMSE among different combinations. The number of samples (n) and number of subjects (subs) included in each analysis are indicated for each sensor combination.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Niswander, W.; Wang, W.; Kontson, K. Optimization of IMU Sensor Placement for the Measurement of Lower Limb Joint Kinematics. Sensors 2020, 20, 5993. https://doi.org/10.3390/s20215993

AMA Style

Niswander W, Wang W, Kontson K. Optimization of IMU Sensor Placement for the Measurement of Lower Limb Joint Kinematics. Sensors. 2020; 20(21):5993. https://doi.org/10.3390/s20215993

Chicago/Turabian Style

Niswander, Wesley, Wei Wang, and Kimberly Kontson. 2020. "Optimization of IMU Sensor Placement for the Measurement of Lower Limb Joint Kinematics" Sensors 20, no. 21: 5993. https://doi.org/10.3390/s20215993

APA Style

Niswander, W., Wang, W., & Kontson, K. (2020). Optimization of IMU Sensor Placement for the Measurement of Lower Limb Joint Kinematics. Sensors, 20(21), 5993. https://doi.org/10.3390/s20215993

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