- 1Biomechanics and Technical Aids Unit, National Hospital for Paraplegics, Toledo, Spain
- 2Biorobotics Group, CAR-Centre of Automation and Robotics, CSIC-Spanish National Research Council, Madrid, Spain
- 3International Doctoral School, Rey Juan Carlos University, Madrid, Spain
When assessing gait analysis outcomes for clinical use, it is indispensable to use an accurate system ensuring a minimal measurement error. Inertial Measurement Units (IMUs) are a versatile motion capture system to evaluate gait kinematics during out-of-lab activities and technology-assisted rehabilitation therapies. However, IMUs are susceptible to distortions, offset and drifting. Therefore, it is important to have a validated instrumentation and recording protocol to ensure the reliability of the measurements, to differentiate therapy effects from system-induced errors. A protocol was carried out to validate the accuracy of gait kinematic assessment with IMUs based on the similarity of the waveform of concurrent signals captured by this system and by a photogrammetry reference system. A gait database of 32 healthy subjects was registered synchronously with both devices. The validation process involved two steps: 1) a preliminary similarity assessment using the Pearson correlation coefficient, and 2) a similarity assessment in terms of correlation, displacement and gain by estimating the offset between signals, the difference between the registered range of motion (∆ROM), the root mean square error (RMSE) and the interprotocol coefficient of multiple correlation (CMCP). Besides, the CMCP was recomputed after removing the offset between signals (CMCPoff). The correlation was strong (r > 0.75) for both limbs for hip flexion/extension, hip adduction/abduction, knee flexion/extension and ankle dorsal/plantar flexion. These joint movements were studied in the second part of the analysis. The ∆ROM values obtained were smaller than 6°, being negligible relative to the minimally clinically important difference (MCID) estimated for unaffected limbs, and the RMSE values were under 10°. The offset for hips and ankles in the sagittal plane reached -9° and -8°, respectively, whereas hips adduction/abduction and knees flexion/extension were around 1°. According to the CMCP, the kinematic pattern of hip flexion/extension (CMCP > 0.90) and adduction/abduction (CMCP > 0.75), knee flexion/extension (CMCP > 0.95) and ankle dorsi/plantar flexion (CMCP > 0.90) were equivalent when captured by each system synchronously. However, after offset correction, only hip flexion/extension (CMCPoff = 1), hip adduction/abduction (CMCPoff > 0.85) and knee flexion/extension (CMCPoff > 0.95) satisfied the conditions to be considered similar.
1 Introduction
Biomechanics is considered an important tool to assess gait during rehabilitation therapies, since it allows the quantitative analysis of human walking, gait features recognition, and its use for diagnostic purposes and treatment planning (Growney et al., 1997; Yavuzer et al., 2008; Xu et al., 2024). However, in technology-assisted rehabilitation therapies, evaluating the immediate biomechanical effects of using technology becomes a challenge. Multiple devices allow assessing biomechanics such as optoelectronic systems based on 3D photogrammetry and Inertial Measurement Units (IMUs) to evaluate kinematics, or force platforms in the case of kinetic analysis. Currently, in neurorehabilitation, the gold standard to assess therapy and intervention outcomes are photogrammetry systems, either with active or passive markers, as they allow in-depth analysis of gait kinematics due to their high accuracy (Eichelberger et al., 2016). Nevertheless, these systems have some drawbacks related to the quality of the recordings, which may be influenced by the number of cameras, the occlusion of markers, the time spent in the execution of the test due to instrumentation, or the expertise of the evaluator (Cano de la Cuerda and Collado Vázquez, 2012; Hassani et al., 2022). In accordance with these limitations, photogrammetry systems are not an option when performing tests to evaluate kinematics in out-of-lab environments, to study the immediate effect of rehabilitation technologies, or when therapies involve multiple devices. Therefore, as a more versatile alternative, IMUs motion capture (mocap) systems are used to evaluate gait kinematics.
IMUs are small and light motion sensor devices based on micro-electro-mechanical technology that estimate the orientation of a body segment to which they are attached from the inertial forces experienced by that segment. The orientation of the IMU is expressed with respect to a fixed coordinate system based on the magnetic north and Earth’s gravitational force. In this way, no specially equipped laboratories are necessary to use them. This makes the system portable and useable both outdoors and during technology-assisted therapies (Ferrari et al., 2010b; Francisco and Tejada, 2020; Hassani et al., 2022).
However, registration will be valid as long as the magnetic field is not distorted (de Vries et al., 2009). Besides, IMUs also suffer from drifting biases, a type of cumulative noise in their measurements that hinders an integration-based analysis to estimate kinematics (Sabatini, Ligorio, and Mannini, 2015; Dorschky et al., 2019). In recent years, multiple approaches based on filtering and global optimization have been proposed to cope with sensor noise and drift to correctly estimate the relative position and orientation of each body segment. For example, incorporating an extended Kalman filter in the sensory fusion process to obtain the corrected orientation of each sensor and segment, or applying a Gaussian distribution to model accelerometer noise and gyroscope bias (Roetenberg et al., 2005; Kok, Hol, and Schön, 2014). Nevertheless, their implementation does not always achieve drift-free estimation of joint angles. In addition, unlike photogrammetry systems, it does not consider the anthropometric measures of the users, there is no standardized placement zone, and its angular measurements always start from zero regardless of the initial posture. These lead to greater sources of error, causing offset and not recording the real joint range of motion (ROM) in their registrations, therefore biomechanical constraints are included in previous studies (Kok, Hol, and Schön, 2014; Dorschky et al., 2019). All these factors make the system susceptible to distortions, so it is important to design and validate an instrumentation and recording protocol able to control these aspects.
Reliability of gait parameters with minimal measurement error is an important consideration in the clinical use of quantitative gait analysis outcomes (Yavuzer et al., 2008), therefore gait analysis requires an accurate, reproducible and precise measurement system (Hassani et al., 2022). It is important to investigate whether a variation between measurements is a therapy effect or is solely due to variation in registrations (Hammer and Lindmark, 2003). Significant information will be lost if recording errors mask gait impairments. Therefore, before using IMUs to measure gait kinematics, it is indispensable to validate the technology with a reliable reference system. Common approaches are to use an optoelectronic mocap system as a reference due to its proven precision (Ferrari et al., 2010b; Kim and Nussbaum, 2013; Zhang et al., 2013; Schiefer et al., 2014).
For all these reasons, the aim of this study was to present a validation protocol to evaluate the accuracy of gait kinematic assessment with IMUs based on the similarity of the waveform of the signals captured by this system and those captured by a photogrammetry system synchronously during gait tests. For this purpose, concurrent measurements were taken with IMUs and the reference photogrammetry system in individuals without gait disorders, since they have a repeatable gait pattern that allows comparing the equivalence of the waveforms taken by both devices and determining the reliability of the IMUs system (Kadaba et al., 1989).
2 Materials and methods
2.1 Participants
A gait database of 32 healthy adult subjects was gathered, volunteers were between 20 and 63 y. o. (33.64 ± 12.44) and 71.88% were females. The detailed demographic and clinical characteristics of the sample are presented in Table 1. The data were collected between June and November 2023. Every individual underwent a barefoot walking test recorded simultaneously with IMUs and the photogrammetry system. The dataset contains the kinematic gait information of the hip, knee, and ankle joints in the three planes of motion: sagittal, frontal, and transversal (Blanco-Coloma et al., 2024).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-t001.jpg)
Table 1. Demographic and anthropometric characteristics of the 32 subjects that make up the healthy gait database for validation.
All subjects were informed of the purpose of the study, the possibility of withdrawing from the same, and signed an informed consent for gait analysis. The study protocol was approved by the local bioethics committee (Clinical Research Ethics Committee at University Hospital Complex of Toledo, CEIC-CHTO-NO 1006 of 26 April of 2023 and NO 949 of 25 January of 2023).
2.2 Experimental procedure and data acquisition
Each subject was instrumented with 8 IMUs of the Tech-MCS V3 mocap system (Technaid S.L., Spain), and with 23 passive markers of Vicon photogrammetry system (Vicon Motion System, Oxford, United Kingdom), following the Plug-in Gait marker set model (Motion Capture System, 2017; Plug-in Gait Reference Guide, 2020). Two additional markers were placed on the medial condyles to adjust femur rotation during processing (Figure 1). For each subject a maximum of 15 captures were registered.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g001.jpg)
Figure 1. IMUs and the photogrammetry system Vicon marker set instrumentation for synchronous captures recording.
For IMUs instrumentation, a standardized model was developed according to the anthropometric measures of each subject. Thigh sensors were placed in the upper third of the segment, displaced frontally 5 cm from the vertical line formed by the trochanter and the lateral condyle of the knee. Tibia sensors were placed in the upper third of the segment slightly displaced towards the inner part, laterally touching the tibial spine. Ankle sensors were placed in the instep, the lumbar sensor at L4, and the chest sensor on the sternum, close to the clavicle. The axes orientation of the sensors was also defined and equal between subjects. The instrumentation with IMUs and the photogrammetry marker set is shown in Figure 1.
The calibration position of the system was fixed to avoid introducing offsets in the captures. This position was as follows: arms extended in a T-shape, trunk and legs extended and opened to the width of the hips, keeping the ankles completely aligned with the tibia segment in a neutral 0° position.
2.3 Data analysis
2.3.1 Signal processing
Gait kinematics recorded with the photogrammetry system were preprocessed with Nexus 2.10.3 software (Vicon Motion Systems, Oxford, United Kingdom). Standard processing operations combined with the anthropometric data entered in the system allowed estimating the position of the joint centers and, subsequently, obtaining the kinematic trajectories of each joint angle of the lower extremities in the three planes of motion, resulting in a personalized and accurate gait analysis. Likewise, captures recorded with the IMUs were transformed from quaternions to Euler angles with the Tech MCS software (Technaid S.L., Spain). Afterward, the trials registered with each system were exported and further processed in MATLAB_R2021b software (The MathWorks, Inc., Natick, Massachusetts, United States). The IMUs signals were smoothed with the Saviztky-Golay filter of order three and with frame length of 21 samples. The photogrammetry signals were captured with a sampling frequency of 100 Hz; therefore, these were resampled to the IMUs sampling frequency, 50 Hz, and aligned with the corresponding IMUs signal. Henceforth, in this article, the signals recorded by the IMUs are referred to as I(t) and those recorded by Vicon V(t). Then, gait cycles were extracted for each pair of trials, obtaining a pair of waveforms for each gait cycle [I(t), V(t)]. A total of 268 synchronous gait cycles [I(t), V(t)] were recovered and analyzed in the three planes of motion per limb.
2.3.2 Waveform similarity assessment
To validate the robustness and feasibility of the IMUs configuration, the similarity and variability of the waveforms of the extracted gait cycles was studied.
To determine whether these two mocap devices were interchangeable (i.e., equivalent) for measuring kinematics, the evaluation of similar waveforms was performed. To assess the similarity between I(t) and V(t) in terms of correlation, displacement, and gain, four parameters were calculated for each [I(t), V(t)] and each joint-angle: their Pearson correlation coefficient (r), the offset between I(t) and V(t), the difference between the registered range of motion (∆ROM) and the root mean square error (RMSE) (Ferrari et al., 2010b). The formulas used to calculate each are presented in Equations 1–3.
In addition, the adjusted variation of the within-day coefficient of multiple correlation (CMC) of Kadaba, named as the interprotocol CMC (CMCP), was calculated. This parameter assesses the repeatability of kinematics removing all other sources of “gait-cycle-to-gait-cycle” variability: 1) the biological variability of the subject’s lower limb kinematics, 2) the variability in the spread of soft tissue artefact in the lower limb kinematics and 3) the variability in the performance of the measurement system (Mcginley et al., 2009; Ferrari et al., 2010a; Ferrari et al., 2010b). Given that each I(t) can be compared only with its synchronous V(t), the aim of this new CMCP statistic formulation is to assess the similarity of waveform (joint angles) acquired synchronously with different protocols and different measurement systems, within each gait cycle, when the effect of the media on waveform similarity is the only parameter of interest. It considers the magnitude of the waveform data and provides a value between 0 and 1, with a value of 1 indicating perfect similarity between two waveforms, [I(t), V(t)]. The formulation shown in Equation 4 was used to evaluate the interprotocol similarity, the CMCP. Suppose that for a subject and a joint angle, the kinematics are measured synchronously through P protocols, in G gait cycles. Consequently, P waveforms are available for each gth gait cycle, one per protocol, each of Fg frames (Ferrari et al., 2010a).
Where
If within each gait cycle, the variability of the P waveforms around their mean waveform is less than the variance around their overall mean, the CMCP approaches one. Otherwise, the CMCP tends to zero or even turns into a complex number. This happens, for instance, when the ROM of the P waveforms is comparable to the phase difference (offset) among them (Kadaba et al., 1989; Ferrari et al., 2010a). To interpret the CMCP and r values obtained, the following ranges were considered: poor (0–0.65), moderate (0.65–0.75), good (0.75–0.8), very good (0.85–0.95) and excellent (0.95–1) (Yavuzer et al., 2008; Ferrari et al., 2010b).
The CMCP considers the overall effect of the offset, r, and gain between waveforms, but it has limitations when recording gait curves with low ROM, resulting in complex values, whose interpretation is not agreed upon, nor evident in the formula breaking down (Røislien et al., 2012). Unlike other studies that drew conclusions from validation by focusing mainly on this parameter, this protocol established a stage-by-stage analysis to identify the joint angles that have equivalent waveforms with both systems by analyzing each parameter independently and assessing how they influence the CMCP value.
The first part of the waveform similarity assessment between the signals I(t) and V(t) was based on the calculation of their correlation. The analysis continued for those joint angles that meet the following requirement:
- Condition 1: The median value of r should follow a strong positive tendency (>0.7) (Kotu and Deshpande, 2019). For this reason, it was established that in the sagittal plane, the median value of r should belong at least to the very good range (0.85–0.95). In the frontal and transversal plane, the median value of r should belong at least to a good range (0.75–0.85).
Next, for those joint angles that met condition 1, the remaining parameters related to displacement and gain were calculated: the offset, the ∆ROM and the RMSE; as well as the CMCP. The signals I(t) and V(t) were totally equivalent, and therefore the recording systems completely substitutable, for those movement planes that satisfied the following condition:
- Condition 2: In the sagittal plane, the median values of CMCP should belong to the excellent range (0.95–1). In the frontal and transversal planes, the CMCP median values should belong at least to the very good range (0.85–0.95).
Conditions had higher acceptance thresholds for the sagittal plane because there is evidence that it is the plane with most reliable and repeatable kinematics, especially for the hip and knee (Kadaba et al., 1989; Mcginley et al., 2009). Furthermore, photogrammetry, the gold standard for kinematic assessment, has higher inter-trial, intra- and inter-evaluator precision and reliability in the sagittal plane than in the other two (Fonseca et al., 2023). Thus, the acceptance thresholds were lower for frontal and transversal motion planes.
The values of all parameters are presented with box-and-whisker plots as well as in terms of median and whisker range for each limb and each joint angle, allowing the variability and dispersion of the recorded data to be studied. Median and whiskers were used because all the parameters did not follow a normal distribution for all joint angles; normality was tested with Lilliefors.
Besides, for those planes that satisfied condition 2, the CMCP was recomputed after zeroing the offset (CMCPoff) for each couple [I(t), V(t)] to measure the effect of the displacement on the similarity. The offset was corrected by subtracting from each gait cycle I(t) the offset between it and its partner V(t), that was previously calculated with Equation 1 (Kadaba et al., 1989).
At last, a visualization of the kinematics recorded is given for every joint angle before correcting the offset, including those that did not satisfy condition 1, and after removing the offset, for those that were analyzed in condition 2.
3 Results
3.1 Waveform similarity assessment: condition 1
The results of the r parameter calculated for each pair of gait cycles [I(t), V(t)] are shown in Figure 2. It shows the distribution of r for each joint angle considering all acquired gait cycles. The figure shows two boxplots, one for each leg, with data from 268 pairs of gait cycles [I(t), V(t)] each. In total, each boxplot contains 268*9 = 2412 values (9 joint angles). Data are displayed in boxplots because not all parameters follow a normal distribution. The results of the normality assessment performed for each parameter with the Lilliefors test are given in Table 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g002.jpg)
Figure 2. Box-and-whisker plot for r regarding the comparison between the 9 joint angles acquired with IMUs and the photogrammetry system. The median value for each joint angle is reported.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-t002.jpg)
Table 2. Results of the Lilliefors normality test computed for the parameters estimated for each pair of gait cycles [I(t), V(t)] and each joint angle.
For the sagittal plane, both hips had a median r value of 0.99 and both knees 0.96, showing an excellent correlation. Thus, there was a positive direct relationship between the gait cycles of each couple [I(t), V(t)], showing a similar waveform (Kotu and Deshpande, 2019). For the ankle dorsi/plantar flexion the median r values decreased slightly but were also within the very good correlation range, being 0.86 for the right limb and 0.87 for the left one. The three joints of both limbs complied condition 1.
In the frontal plane both hips had a strong correlation, showing a good correlation for the right hip adduction/abduction movement (r = 0.75) and a very good correlation for the left one (r = 0.86), satisfying condition 1, despite having more scattered data and more outliers. However, in the frontal plane for the knees (right: r = 0.22, left: r = 0.32) and ankles (right: r = 0.09, left: r = 0.08) the r medians presented low values with a high dispersion.
Finally, the hip rotation in both limbs presented a poor correlation, being the median values of r 0.15 and −0.07 for the right and left limbs, respectively. Equally, the knees (right: r = −0.41, left: r = −0.35) and the ankles (right: r = −0.41, left: r = −0.30) in the transversal plane showed no similarity between [I(t), V(t)] gait cycles and a lot of dispersion in the data.
3.2 Waveform similarity assessment: condition 2
The results obtained for the ∆ROM, offset, RMSE and CMCP estimations are shown for those joint angles that satisfied condition 1: hip flexion/extension, hip adduction/abduction, knee flexion/extension and ankle dorsal/plantar flexion.
When the resulting value of the ∆ROM is negative, it means that the range measured by the IMUs was smaller than that measured by the photogrammetry system, while if it is positive, the opposite was true. As seen in Figure 3, in the sagittal plane, the median ∆ROM values for the hips were −2.10° in the right limb and 1.71° in the left one, showing slightly more dispersion and more outliers. In the case of the knees, the ROM registered by the IMUs was higher than the ROM registered by the photogrammetry for both limbs (right: ∆ROM = 5.14°, left: ∆ROM = 6.86°). For the ankles in the sagittal plane (right: ∆ROM = -1.14°, left: ∆ROM = 0.65°) and the hip in the frontal plane (right: ∆ROM = 1.33°, left: ∆ROM = 0.32°) smaller differences were observed.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g003.jpg)
Figure 3. Box-and-whisker plot for ∆ROM regarding the comparison between the 4 joint angles that satisfy condition 1: hip flexion/extension, hip adduction/abduction, knee flexion/extension and ankle dorsi/plantar flexion. The median value for each joint angle is reported.
The offset box-plot graph is shown in Figure 4. According to the offset formula, when the resulting value is negative, it means that the IMUs measurement had a negative offset (lower values) with respect to the photogrammetry system, and if it is positive, the opposite happened. For the hip, the sagittal plane had a median offset value of −9.06° for the right limb and −7.98° for the left one. However, the offset in the frontal plane was lower (right: 0.18°, left: −0.75°). The knee flexion/extension also presented a small offset for both limbs, being the median values −1.02° and 1.59° for right and left, respectively. Besides, both ankles showed a negative offset of -8°. In this case, the hips and ankles in the sagittal plane had more data dispersion and outliers.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g004.jpg)
Figure 4. Box-and-whisker plot for offset regarding the comparison between the 4 joint angles that satisfy condition 1: hip flexion/extension, hip adduction/abduction, knee flexion/extension and ankle dorsi/plantar flexion. The median value for each joint angle is reported.
Analyzing the RMSE, depicted in the Figure 5, both hips and ankles in the sagittal plane presented more scattered data and the median RMSE values were around 10°. The knees presented a median error of 7°, and the hip adduction/abduction movement showed an error of 6.28° in the right leg and 4.26° in the left leg.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g005.jpg)
Figure 5. Box-and-whisker plot for RMSE regarding the comparison between the 4 joint angles that satisfy condition 1: hip flexion/extension, hip adduction/abduction, knee flexion/extension and ankle dorsi/plantar flexion. The median value for each joint angle is reported.
In terms of the CMCP calculation (Figure 6), whose formula shown in Equation 4 included the whole effect of the offset, r, and ∆ROM, the knee flexion/extension movement was the only one that reached an excellent median value of 0.99 for both limbs with hardly any data dispersion, satisfying condition 2. For the hip in the sagittal plane, the left one also fulfilled condition 2, displaying an excellent median CMCP value of 0.96. The right hip also presented a very good CMCP of 0.91 between pairs of [I(t), V(t)]. For the frontal plane, the right hip had a good CMCP value of 0.78, whereas the left one reached a very good coefficient value of 0.85. Lastly, both ankles presented median CMCP values around 0.9, showing a very good similarity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g006.jpg)
Figure 6. Box-and-whisker plot for CMCP regarding the comparison between the 4 joint angles that satisfy condition 1: hip flexion/extension, hip adduction/abduction, knee flexion/extension and ankle dorsi/plantar flexion. The median value for each joint angle is reported.
The effects of correcting the offset are shown in Figure 7. Once the offsets between corresponding waveforms were removed, the CMCPoff values improved for every joint angle, except for the ankles dorsi/plantar flexion (right: CMCPoff = 0.74, left: CMCPoff = 0.64), whose data dispersion also increased considerably. For both limbs in the sagittal plane the hips and knees flexion/extension described an excellent similarity of 1 and 0.99, respectively. The dispersion of the data enhanced significantly after removing the offset for the hips in the sagittal plane. The hips CMCPoff also increased in the frontal plane, becoming 0.86 for the right limb and 0.91 for the left one, both in the very good similarity range. Therefore, when the offset was corrected, condition 2 was fulfilled by hip flexion/extension and adduction/abduction, and by knee flexion/extension for both limbs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g007.jpg)
Figure 7. Box-and-whisker plot for CMCPoff regarding the comparison between the 4 joint angles that satisfy condition 1: hip flexion/extension, hip adduction/abduction, knee flexion/extension and ankle dorsi/plantar flexion. The median value for each joint angle is reported.
In addition, the kinematic pattern of each pair of gait cycles [I(t), V(t)] for the 9 joint angle is shown in Figures 8, 9, containing kinematic data of the right and left limbs, respectively. Lastly, the kinematic pattern with the offset arrangement is displayed only for those joint angles that have been analyzed in condition 2. Figure 10 shows the offset correction for the right limb joints and Figure 11 for the left ones.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g008.jpg)
Figure 8. Right limb kinematic pattern of the gait cycles registered with IMUs (yellow) and the photogrammetry system (purple) in the 9 joint angles.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g009.jpg)
Figure 9. Left limb kinematic pattern of the gait cycles registered with IMUs (yellow) and the photogrammetry system (purple) in the 9 joint angles.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g010.jpg)
Figure 10. Right limb kinematic pattern after offset zeroing of the gait cycles registered with IMUs (yellow) and the photogrammetry system (purple) in the 4 joint angles that satisfy condition 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1449698/fbioe-12-1449698-HTML/image_m/fbioe-12-1449698-g011.jpg)
Figure 11. Left limb kinematic pattern after offset zeroing of the gait cycles registered with IMUs (yellow) and the photogrammetry system (purple) in the 4 joint angles that satisfy condition 1.
4 Discussion
The kinematic pattern of the hip flexion/extension and adduction/abduction, the knee flexion/extension and the ankle dorsi/plantar flexion were equivalent when captured by each mocap system synchronously. However, when the offset was corrected, only the hip flexion/extension, the hip adduction/abduction and the knee flexion/extension satisfied the conditions to be considered similar with high confidence.
Although the CMCP included the effect of every estimated parameter, it was considered necessary to study individually the variability, trend, and dispersion of each of them to properly evaluate the accuracy of the synchronous pairwise measurements [I(t), V(t)]. Likewise, it was decided to study each parameter individually to identify the root cause of the differences in the measurements registered by each mocap system, since some of them could be controlled by adjusting the IMUs instrumentation model or the registration protocol, thus improving the precision of the movement recorded with the IMUs.
Initially, given the limitations observed in CMC in previous papers and in this work, it was decided to perform an initial analysis focused uniquely on the assessment of the correlation-centered waveform similarity (Growney et al., 1997; Steinwender et al., 2000; Røislien et al., 2012). If the waveforms of each pair of gait cycles [I(t), V(t)] did not maintain a linear dependence relationship, that is, they did not change in the same way, it was assumed that the kinematic pattern was not similar, and the in-depth analysis was not performed.
As stated in the results, the only joint angles that showed a sufficiently strong correlation (r > 0.75) to meet condition 1 were the movements of all joints (hips, knees, and ankles) in the sagittal plane and the hips in the frontal plane (Figure 2). In addition, these joint boxes showed less data dispersion and thus less variability in the recordings. Furthermore, this result is confirmed by the graphs in Figures 8, 9, which show the kinematic pattern of all the pairs of gait cycles registered with the IMUs and the photogrammetry system for all joint angles. It can be observed that all the signals in the sagittal plane and the hips in the frontal plane follow the same waveform behavior with both systems. Some outliers are also observed in the dorsi/plantar flexion of the right ankle, showing a registration error by the IMUs. The hypothesis of the origin of this error is the distortions of the magnetometers and the drift accumulated during the captures.
Besides, it can also be seen that the IMUs and photogrammetry signals of the remaining joint angles did not have the same waveform, confirmed by their poor correlation (r < 0.65). In other studies it has been demonstrated that the kinematics recorded in the sagittal plane in adult gait are the most repeatable within the same day and between days, while the repeatability of gait kinematic pattern in the other planes is much more variable for the same subject even within the same day, mainly in the transverse plane for every joint and knee valgus/varus (Kadaba et al., 1989; Growney et al., 1997; Besier et al., 2003; Mcginley et al., 2009). In this study, this can be appreciated especially in hip rotation and knee valgus/varus movements (Figures 8, 9), where the average kinematic pattern recorded with photogrammetry does not describe a common trend, showing more data variability, and does not correspond to that defined in the theory (van der Linden, 2011; Duarte et al., 2018; Francisco and Tejada, 2020). Additionally, there is evidence that, despite their accuracy, photogrammetry systems are less reliable in the transverse plane and knee varus/valgus during gait analysis (Fonseca et al., 2023). In fact, most studies report the highest errors in those planes (Growney et al., 1997; Kim and Nussbaum, 2013; Kok, Hol, and Schön, 2014). So, adjusting the accuracy of other less robust devices in these planes based on the similarity assessment could not be as reliable.
The second part of the analysis focused on the four planes that satisfied condition 1. The position of each IMU sensor was defined so that it would record the ROM as completely as possible. Looking at Figure 3, the ∆ROM recorded for hips flexion/extension, hips adduction/abduction and ankles dorsi/plantar flexion barely reached 2°. These values are closed to those obtained in similar IMUs and photogrammetry systems validation studies (Picerno, Cereatti, and Cappozzo, 2008; Ferrari et al., 2010b) and, in addition, fall within the ranges considered despicable according to the minimal clinically important differences (MCID) estimated for unaffected limbs (Guzik et al., 2021). The knee flexion/extension difference values reached 6.8° in the left limb, higher than those reported in the mentioned validation studies. However, this is still a negligible value with respect to the MCID estimated in clinical practice for unaffected limbs in other articles, which almost reaches 7° (Guzik et al., 2020). It should also be noted that previous assessment studies include samples of one to four subjects, whereas in this work the sample was increased to 32 users, therefore, more variability was captured in the dataset.
Data from the articles reporting errors reveal that most of the studies and gait variables present errors between 2° and 5° for every joint angle, although few of them reach the sample size presented here. The lower RMSE obtained in this validation, in Figure 5, was the one corresponding to the hips adduction/abduction, whereas in the sagittal plane the knees reached 7° and the hips and the ankles reached 10°. However, there are also studies reporting errors between 5° and 10° in the sagittal plane during clinical assessment in the lower limb (McDowell et al., 2000; Fosang et al., 2003). In this study, therefore, these errors can be considered negligible.
As explained in the methodology, the mode of computing CMCP is significantly influenced by joint ROM. Previous studies note that joints with a large ROM tend to record high CMCP and, conversely, joints with a low ROM tend to show lower reliability (Growney et al., 1997; Steinwender et al., 2000; Mcginley et al., 2009; Røislien et al., 2012). This limitation was also shown in this validation (Figure 6), because the lowest CMCP values were reported for the joint angles with smaller ROMs evaluated which were hip adduction/abduction (CMCP = 0.8) and the ankle dorsi/plantar flexion (CMCP = 0.9), even though the measured error and waveform correlation values obtained were good. However, in this study the reliability value increased in the frontal plane because most studies report minimum reliability indices of 0.7. This could be thanks to the standardized instrumentation proposed for the IMUs, defined following the anthropometry measurements to enhance reliability. For the hip and knee in the sagittal plane, the best repeatability and reliability values were obtained (CMCP > 0.9), as reported in the other articles, confirming the feasibility of recording precise kinematics of this plane with IMUs (Mcginley et al., 2009; Ferrari et al., 2010b; Robert-Lachaine et al., 2017).
On the other hand, the CMCPoff was computed to assess how the offset affects the calculation of the repeatability (Figure 7). Once removed, the CMCPoff of the hips in the sagittal plane became excellent, this movement was the one with more offset registered (right: offset = −9.06°; left: offset = −7.98°), and in the frontal plane it increased to 0.9, demonstrating that reliability was directly influenced by the offset. The dispersion of the data, and therefore its variability, improved significantly. Meanwhile, the knees were not disturbed by the offset, whose value is 1°, since they maintained an excellent reliability of 0.99. On the contrary, the CMCPoff of the ankles decreased once the offset is removed. As the CMCP includes the effect of the r, the offset and ∆ROM, if the CMCPoff did not improve when correcting the offset, the root of the problem is another of these parameters. However, as previously stated, the values of ∆ROM, RMSE and r were good. In addition, it was found that most of the values obtained for the [I(t), V(t)] couples were complex, thus, the hypothesis is that the CMCPoff estimation was limited in this movement because of the small ROM of this joint angle, resulting in the formula breaking down (Røislien et al., 2012). As there is no consensus for interpreting these complex values, the acceptance of ankle flexion/extension as reliable was centered on all the other parameters studied. For the CMCP with offset included, this plane showed a value of 0.92 and 0.91 for right and left extremities respectively, showing a very good reliability range, despite not fulfilling condition 2.
According to the IMUs instrumentation protocol, it is known that the offset is caused due to two main reasons: the difference in degrees introduced due to anatomical morphology and the difference in degrees introduced during calibration. In both cases these introduced degrees remained constant during all the captures for the same subject. Likewise, both hips and ankles showed a constant offset for all the subjects, as can be seen in Figures 4, 8, 9, and in the kinematic patterns whose offsets were corrected, shown in Figures 10, 11. Considering that it could be arranged in post-processing to achieve a satisfactory reliability for gait recordings with IMUs, it was agreed that values obtained with CMCPoff satisfy the requirements established in the methodology.
For experiments in which the IMUs mocap system is used to register healthy adults kinematic gait data, it is proposed to perform the offset correction in post-processing by subtracting to the whole signal the median offset value obtained in this study for each joint angle, shown in Figure 4, or by measuring with a goniometer the offsets of each joint angle intrinsic to anatomical morphology that could be observed before starting the recording.
5 Conclusion
The kinematic pattern of the hip and the knee in the sagittal plane and the hip in the frontal plane satisfies condition 2 showing an excellent and a very good similarity, respectively, between the waveforms captured with the IMUs and those captured with the photogrammetry system synchronously. Therefore, it is concluded that the movements registered by the IMUs in those planes are completely reliable with offset included. However, the reliability increases if the offset correction is performed, especially for the hip adduction/abduction movement. On the other hand, the ankle dorsi/plantar flexion shows a very good range of similarity between the IMUs and photogrammetry records. Although it does not fulfill condition 2, it can be used to obtain measurements of this joint angle with very good reliability.
It is also concluded that standardizing the instrumentation of the IMUs to the anthropometric measurements of each subject favors the complete recording of the ROM of each joint. However, this is not enough to adjust the offset control, since it also depends on a good calibration position in which the joint angles are aligned at 0°.
Data availability statement
The raw data registered for this study can be found in Zenodo, an open research data repository, at this link https://doi.org/10.5281/zenodo.11242602.
Ethics statement
The studies involving humans were approved by the Clinical Research Ethics Committee at University Hospital Complex of Toledo, CEIC-CHTO-NO 1006 of 26 April of 2023 and NO 949 of 25 January of 2023. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
LB-C: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing. LG-G: Investigation, Software, Supervision, Writing–review and editing. IS-A: Investigation, Supervision, Writing–review and editing. ST-Á: Data curation, Investigation, Supervision, Writing–review and editing. PM-H: Data curation, Investigation, Supervision, Writing–review and editing. SG-E: Data curation, Investigation, Supervision, Writing–review and editing. ÁG-A: Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Visualization, Writing–review and editing. DH-V: Conceptualization, Data curation, Funding acquisition, Investigation, Project administration, Resources, Software, Visualization, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research is part of grant PID2021-124111OB- C33, funded by MCIN/AEI/10.13039/501100011033 and by ERDF A way of making Europe.
Acknowledgments
We acknowledge all the healthy volunteers who accepted voluntarily to be part of this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Besier, T. F., Sturnieks, D. L., Alderson, J. A., and Lloyd, D. G. (2003). Repeatability of gait data using a functional hip joint Centre and a mean helical knee Axis. J. Biomechanics 36 (8), 1159–1168. doi:10.1016/S0021-9290(03)00087-3
Blanco-Coloma, L., García-González, L., Sinovas-Alonso, I., Herrera-Valenzuela, D., and Gil-Agudo, Á. (2024). Synchronously recorded gait kinematic data with inertial measurement units and a photogrammetry system for a validation assessment. Zenodo 2024. doi:10.5281/ZENODO.11242602
Cano de la Cuerda, R., and Collado Vázquez, S. (2012). “Neurorrehabilitación. Métodos Específicos de Valoración y Tratamiento,” in Neurorrehabilitación. Métodos Específicos de Valoración y Tratamiento. 1st ed. (Madrid: byEditorial Médica Panamericana SA), 1, 161–172. doi:10.1016/J.RH.2012.02.003
de Vries, W. H. K., Veeger, H. E. J., Baten, C. T. M., and van der Helm, F. C. T. (2009). Magnetic distortion in motion labs, implications for validating inertial magnetic sensors. Gait and Posture 29 (4), 535–541. doi:10.1016/J.GAITPOST.2008.12.004
Dorschky, E., Nitschke, M., Seifer, A. K., van den Bogert, A. J., and Eskofier, B. M. (2019). Estimation of gait kinematics and kinetics from inertial sensor data using optimal control of musculoskeletal models. J. Biomechanics 95 (October), 109278. doi:10.1016/J.JBIOMECH.2019.07.022
Duarte, M., Fukuchi, C. A., and Fukuchi, R. K. (2018). A public dataset of overground and treadmill walking kinematics and kinetics in healthy individuals. doi:10.7717/peerj.4640
Eichelberger, P., Ferraro, M., Minder, U., Denton, T., Blasimann, A., Krause, F., et al. (2016). Analysis of accuracy in optical motion capture – a protocol for laboratory setup evaluation. J. Biomechanics 49 (10), 2085–2088. doi:10.1016/J.JBIOMECH.2016.05.007
Ferrari, A., Cutti, A. G., and Cappello, A. (2010a). A new formulation of the coefficient of multiple correlation to assess the similarity of waveforms measured synchronously by different motion analysis protocols. Gait and Posture 31, 540–542. doi:10.1016/j.gaitpost.2010.02.009
Ferrari, A., Cutti, A. G., Garofalo, P., Raggi, M., Heijboer, M., Cappello, A., et al. (2010b). First in vivo assessment of ‘“Outwalk”’: a novel protocol for clinical gait analysis based on inertial and magnetic sensors. Med. Biol. Eng. Comput. 48, 1–15. doi:10.1007/s11517-009-0544-y
Fonseca, M., Gasparutto, X., Grouvel, G., Bonnefoy-Mazure, A., Dumas, R., and Armand, S. (2023). Evaluation of lower limb and pelvic marker placement precision among different evaluators and its impact on gait kinematics computed with the conventional gait model. Gait and Posture 104, 22–30. doi:10.1016/j.gaitpost.2023.05.028
Fosang, A. L., Galea, M. P., McCoy, A. T., Reddihough, D. S., and Story, I. (2003). Measures of muscle and joint performance in the lower limb of children with cerebral palsy. Dev. Med. and Child Neurology 45 (10), 664–670. doi:10.1111/J.1469-8749.2003.TB00868.X
Francisco, M. R., and Tejada, M. C. (2020). La marcha humana: biomecánica, evaluación y patología. Madrid: Editorial Médica Panamericana S.A.
Growney, E., Meglan, D., Johnson, M., Cahalan, T., and An, K.-N. (1997). Repeated measures of adult normal walking using a video tracking system. Gait Posture 6, 147–162. doi:10.1016/s0966-6362(97)01114-4
Guzik, A., Drużbicki, M., Perenc, L., Wolan-Nieroda, A., Turolla, A., and Kiper, P. (2021). Establishing the minimal clinically important differences for sagittal hip range of motion in chronic stroke patients. Front. Neurology 12 (September), 700190. doi:10.3389/FNEUR.2021.700190
Guzik, A., Drużbicki, M., Wolan-Nieroda, A., Turolla, A., and Kiper, P. (2020). Estimating minimal clinically important differences for knee range of motion after stroke. J. Clin. Med. 9 (10), 3305–3314. doi:10.3390/JCM9103305
Hammer, A., and Lindmark, B. (2003). Test-retest intra-rater reliability of grip force in patients with stroke. J. Rehabilitation Med. 35 (4), 189–194. doi:10.1080/16501970306132
Hassani, R. H., Willi, R., Rauter, G., Bolliger, M., and Seel, T. (2022). Validation of non-restrictive inertial gait and kinematic analysis of incomplete spinal cord injured patients in clinical settings. Submitt. Sensors, 1–12. doi:10.3390/s22114237
Kadaba, M. P., Ramakrishnan, H. K., Wootten, M. E., Gainey, J., Gorton, G., and Cochran, G. V. B. (1989). Repeatability of kinematic, kinetic, and electromyographic data in normal adult gait. J. Orthop. Res. 7 (6), 849–860. doi:10.1002/JOR.1100070611
Kim, S., and Nussbaum, M. A. (2013). Performance evaluation of a wearable inertial motion capture system for capturing physical exposures during manual material handling tasks. Ergonomics 56 (2), 314–326. doi:10.1080/00140139.2012.742932
Kok, M., Hol, J. D., and Schön, T. B. (2014). An optimization-based approach to human body motion capture using inertial sensors. IFAC Proc. Vol. (IFAC-PapersOnline) 19, 79–85. doi:10.3182/20140824-6-za-1003.02252
McDowell, B. C., Vera, H., Nurse, A., Weston, T., and Baker, R. (2000). The variability of goniometric measurements in ambulatory children with spastic cerebral palsy. Gait Posture 12 (2), 114–121. doi:10.1016/S0966-6362(00)00068-0
Mcginley, J. L., Baker, R., Wolfe, R., and Morris, M. E. (2009). The reliability of three-dimensional kinematic gait measurements: a systematic review. Gait and Posture 29, 360–369. doi:10.1016/j.gaitpost.2008.09.003
Motion Capture System (2017). Motion capture system, technaid - inertial motion capture. Available at: https://www.technaid.com/products/motion-capture-system-tech-mcs-hub-imu/.
Picerno, P., Cereatti, A., and Cappozzo, A. (2008). Joint kinematics estimate using wearable inertial and magnetic sensing modules. Gait and Posture 28 (4), 588–595. doi:10.1016/J.GAITPOST.2008.04.003
Plug-in Gait Reference Guide (2020). Plug-in gait reference Guide. Available at: https://docs.vicon.com/display/Nexus213/Plug-in+Gait+Reference+Guide.
Robert-Lachaine, X., Mecheri, H., Larue, C., and Plamondon, A. (2017). Validation of inertial measurement units with an optoelectronic system for whole-body motion analysis. Med. Biol. Eng. Comput. 55, 609–619. doi:10.1007/s11517-016-1537-2
Roetenberg, D., Luinge, H. J., Baten, C. T. M., and Veltink, P. H. (2005). Compensation of magnetic disturbances improves inertial and magnetic sensing of human body segment orientation. IEEE Trans. Neural Syst. Rehabilitation Eng. 13 (3), 395–405. doi:10.1109/TNSRE.2005.847353
Røislien, J., Skare, A. O., and Rennie, L. (2012). Evaluating the properties of the coefficient of multiple correlation (CMC) for kinematic gait data. J. Biomechanics 45 (11), 2014–2018. doi:10.1016/j.jbiomech.2012.05.014
Sabatini, A. M., Ligorio, G., and Mannini, A. (2015). Fourier-based integration of quasi-periodic gait accelerations for drift-free displacement estimation using inertial sensors. doi:10.1186/s12938-015-0103-8
Schiefer, C., Ellegast, R. P., Hermanns, I., Kraus, T., Ochsmann, E., Larue, C., et al. (2014). Optimization of inertial sensor-based motion capturing for magnetically distorted field applications. J. Biomechanical Eng. 136 (12), 121008. doi:10.1115/1.4028822
Steinwender, G., Saraph, V., Scheiber, S., Zwick, E. B., Uitz, C., and Hackl, K. (2000). Intrasubject repeatability of gait analysis data in normal and spastic children. Clin. Biomech. 15 (2), 134–139. doi:10.1016/S0268-0033(99)00057-1
van der Linden, M. (2011). “Gait analysis,”in Normal and pathological function. Editors J. Perry, and J. M. Burnfield 2nd Ed. (Physiotherapy journal: Slack Inc., Physiotherapy), 97, 180. 576 Pages, ISBN 978-1-55642r-R766-4. doi:10.1016/J.PHYSIO.2010.05.007
Xu, D., Zhou, H., Quan, W., Jiang, X., Liang, M., Li, S., et al. (2024). A new method proposed for realizing human gait pattern recognition: inspirations for the application of Sports and clinical gait analysis. Gait and Posture 107 (January), 293–305. doi:10.1016/J.GAITPOST.2023.10.019
Yavuzer, G., Ken, O. ¨ Z. O. ¨, Elhan, A., and Stam, H. J. (2008). Repeatability of lower limb three-dimensional kinematics in patients with stroke. Gait and Posture 27, 31–15. doi:10.1016/j.gaitpost.2006.12.016
Keywords: three-dimensional (3D) kinematic gait data, inertial measurement units (IMUs), photogrammetry, waveform similarity assessment, interprotocol coefficient of multiple correlation (CMCP), feasibility, repeatability
Citation: Blanco-Coloma L, García-González L, Sinovas-Alonso I, Torio-Álvarez S, Martos-Hernández P, González-Expósito S, Gil-Agudo Á and Herrera-Valenzuela D (2024) Validation of inertial measurement units based on waveform similarity assessment against a photogrammetry system for gait kinematic analysis. Front. Bioeng. Biotechnol. 12:1449698. doi: 10.3389/fbioe.2024.1449698
Received: 15 June 2024; Accepted: 30 July 2024;
Published: 12 August 2024.
Edited by:
Elvira Padua, Università telematica San Raffaele, ItalyCopyright © 2024 Blanco-Coloma, García-González, Sinovas-Alonso, Torio-Álvarez, Martos-Hernández, González-Expósito, Gil-Agudo and Herrera-Valenzuela. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Laura Blanco-Coloma, laurablancocoloma@gmail.com
†These authors share senior authorship