Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Sagittal plane knee kinematics can be measured during activities of daily living following total knee arthroplasty with two IMU

  • Bradley M. Cornish ,

    Roles Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

    Bradley.cornish@griffith.edu.au

    Affiliations Griffith Centre of Biomedical and Rehabilitation Engineering, Menzies Health Institute Queensland, Griffith University, Gold Coast, Queensland, Australia, School of Health Sciences and Social Work, Griffith University, Gold Coast, Queensland, Australia

  • Laura E. Diamond,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliations Griffith Centre of Biomedical and Rehabilitation Engineering, Menzies Health Institute Queensland, Griffith University, Gold Coast, Queensland, Australia, School of Health Sciences and Social Work, Griffith University, Gold Coast, Queensland, Australia

  • David John Saxby,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliations Griffith Centre of Biomedical and Rehabilitation Engineering, Menzies Health Institute Queensland, Griffith University, Gold Coast, Queensland, Australia, School of Health Sciences and Social Work, Griffith University, Gold Coast, Queensland, Australia

  • David G. Lloyd,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliations Griffith Centre of Biomedical and Rehabilitation Engineering, Menzies Health Institute Queensland, Griffith University, Gold Coast, Queensland, Australia, School of Health Sciences and Social Work, Griffith University, Gold Coast, Queensland, Australia

  • Beichen Shi,

    Roles Data curation, Investigation, Methodology, Writing – review & editing

    Affiliations Griffith Centre of Biomedical and Rehabilitation Engineering, Menzies Health Institute Queensland, Griffith University, Gold Coast, Queensland, Australia, School of Health Sciences and Social Work, Griffith University, Gold Coast, Queensland, Australia

  • Jenna Lyon,

    Roles Conceptualization, Funding acquisition, Project administration, Resources, Writing – review & editing

    Affiliation Stryker Corporation, Kalamazoo, Michigan, Unites States of America

  • Kevin Abbruzzese,

    Roles Conceptualization, Funding acquisition, Methodology, Resources, Writing – review & editing

    Affiliation Stryker Corporation, Kalamazoo, Michigan, Unites States of America

  • Price Gallie,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Writing – review & editing

    Affiliation Coast Orthopaedics and Sports Medicine, Gold Coast, Queensland, Australia

  • Jayishni Maharaj

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

    Affiliations Griffith Centre of Biomedical and Rehabilitation Engineering, Menzies Health Institute Queensland, Griffith University, Gold Coast, Queensland, Australia, School of Health Sciences and Social Work, Griffith University, Gold Coast, Queensland, Australia

Abstract

Knee function is rarely measured objectively during functional tasks following total knee arthroplasty. Inertial measurement units (IMU) can measure knee kinematics and range of motion (ROM) during dynamic activities and offer an easy-to-use system for knee function assessment post total knee arthroplasty. However, IMU must be validated against gold standard three-dimensional optical motion capture systems (OMC) across a range of tasks if they are to see widespread uptake. We computed knee rotations and ROM from commercial IMU sensor measurements during walking, squatting, sit-to-stand, stair ascent, and stair descent in 21 patients one-year post total knee arthroplasty using two methods: direct computation using segment orientations (r_IMU), and an IMU-driven iCloud-based interactive lower limb model (m_IMU). This cross-sectional study compared computed knee angles and ROM to a gold-standard OMC and inverse kinematics method using Pearson’s correlation coefficient (R) and root-mean-square-differences (RMSD). The r_IMU and m_IMU methods estimated sagittal plane knee angles with excellent correlation (>0.95) compared to OMC for walking, squatting, sit-to-stand, and stair-ascent, and very good correlation (>0.90) for stair descent. For squatting, sit-to-stand, and walking, the mean RMSD for r_IMU and m_IMU compared to OMC were <4 degrees, < 5 degrees, and <6 degrees, respectively but higher for stair ascent and descent (~12 degrees). Frontal and transverse plane knee kinematics estimated using r_IMU and m_IMU showed poor to moderate correlation compared to OMC. There were no differences in ROM measurements during squatting, sit-to-stand, and walking across the two methods. Thus, IMUs can measure sagittal plane knee angles and ROM with high accuracy for a variety of tasks and may be a useful in-clinic tool for objective assessment of knee function following total knee arthroplasty.

1. Introduction

Total knee arthroplasty (TKA) is an effective surgical treatment for knee osteoarthritis (OA) that reduces pain, lessens functional impairments, and improves quality of life [1, 2]. Nevertheless, patient outcomes, including post-operative complication rates, knee function, and patient satisfaction can vary due to differences in surgical techniques, prosthetic design, and implant alignment [3, 4]. Recently, robotically assisted Total Knee Arthroplasty (TKA) has improved alignment of prosthetic implants, optimised ligament balancing, increased knee range of motion (ROM), and reduced surgical times [1, 5, 6]. The use of robots has led to lower levels of post-operative pain and analgesic use, greater improvements in pain at 6- and 12-months post-operation, lower episode-of-care costs, and greater patient satisfaction compared to standard surgical approaches [68]. However, accurately assessing post-operative functional outcomes remains difficult.

Post-operative evaluations include clinical assessments of function and ROM, patient reported outcome measures (PROMS), and radiographic assessment of implant position [7]. Though recommended, the use of quantitative functional-based measures is less common due to timely and costly procedures, and consequently, are often inferred from PROMS [9]. However, there is poor concurrent validity between PROMS and functional measures in people with TKA [9]. Patient reports may relate more to pain than performance in functional tests, are confounded by the function and strength of the uninvolved limb, and fail to detect changes in the affected limb post TKA [10]. Additionally, PROMS commonly exhibit a ceiling effect in which they do not reliably discriminate between moderate and excellent outcomes, as clustering occurs at the upper limit of scoring [11].

Quantitative measures of knee kinematics during walking and other activities of daily living (ADL) may better describe knee function and disease status following knee arthroplasty. Prior to surgery, individuals with knee osteoarthritis (OA) typically display less knee excursion during the loading phase of gait and when climbing stairs, with greater knee adduction angle throughout gait when compared to age-matched controls [12, 13]. Following knee arthroplasty, knee mechanics partially restore [14, 15], although functional deficits often remain when compared to healthy controls [16]. Post TKA, individuals typically show reduced knee excursion (i.e., lower functional range of motion) and reduced external knee flexion moments [17, 18], that may be attributed to reduced knee extensor strength [19]. Evaluating knee biomechanics post-operatively during TKA recovery may better describe post-operative knee function, objectively inform rehabilitation decisions, and lend itself to improved surgical techniques. In the case of robotically assisted TKA, objective data such as knee kinematics may be analysed with respect to surgical information such as implant alignment, implant type, ligament balancing, and surgical duration, which can be logged automatically by the robot systems [6].

Current methods for measuring knee kinematics are time-consuming, costly, require specialised equipment such as three-dimensional optical motion capture (OMC), and operator expertise [20]. Optical motion capture is largely unportable and involves lengthy setup and participant preparation that make it impractical for clinical or home use. Wearable sensors, such as inertial measurement units (IMU), may provide an affordable, portable, and practical solution to movement analysis outside the laboratory [2123] that may inform clinical practice [24]. Inertial measurement units combine multiple sensors; typically an accelerometer, a gyroscope, and a magnetometer, into one device that measures orientation, rotations, and accelerations [25]. For individuals post total hip arthroplasty, IMU have been shown to detect improvements in gait following surgery and provide objective data to complement PROMs [26]. Additionally, a recent review suggests that gait parameters and joint kinematics measured using wearable sensors offer greater insight to functional recovery than traditional functional outcome measures such as the timed-up-and go, and static ROM tests [27].

Inertial measurement units have measured joint angles during functional tasks with root-mean-square-difference (RMSD) less than 5 degrees compared to OMC [22, 2830], which were not different to OMC for measurement of knee ROM during walking in individuals post total hip arthroplasty [31]. However, few studies have validated their use during multiple ADL [32] in a clinical population. For those following TKA, IMUs have shown to measure knee flexion ROM with an average error of roughly 4 degrees during clinical tests [33] but was not assessed during functional tasks and ADL that may indicate a patient’s knee function. Evaluating the accuracy of IMUs across multiple ADL would provide support for their use as a clinical measurement tool and may identify tasks during which IMU measurements are accurate and acceptable for objective assessments of knee function in individuals post TKA.

The aim of this study was to compare knee kinematics using IMU sensors and OMC in individuals 1-year post robotically assisted TKA during walking, squatting, sit-to-stand, stair ascent, and stair descent tasks. Knee kinematics from IMU were calculated using an interactive lower limb model (m_imu) and from raw quaternions (r_imu) measured using commercial IMU sensors (Notch Interfaces Inc. New York, USA). We hypothesised that knee angles for the m_IMU would have lower RMSD and higher correlation (R) than the r_IMU when compared to OMC. Furthermore, we hypothesised that RMSD in the sagittal plane would be lower than out-of-plane angles for all activities. We hypothesised that the sagittal plane knee ROM calculated using m_IMU would have smaller differences t compared to OMC than r_IMU.

2. Methods

2.1 Participants

Twenty-one individuals (9F/12M; mean ± standard deviation; age: 69 yrs ± 8; height: 1.71 m ± 0.11; body mass: 93.89 ± 24 kg) one year post robotically assisted TKA were recruited for this study between August and December 2020. All surgeries were performed by a single surgeon (PG), using the MAKO system (Stryker, Michigan, USA). Participants were included if they were twelve months post Mako assisted TKA, were older than 18 years, and self-reported their ability to walk for 30 minutes continuously. Participants were excluded if they were unable to walk independently; had any other serious lower limb injury in the past twelve months; had muscle degenerative disorders or motor impairments; had a body mass index greater than 40. Ethical approval was obtained from the Griffith University Research Ethics Committee and participants provided their written informed consent prior to testing.

2.2 Experimental protocol

Participants were asked to complete five tasks, in a fixed order: (i) walking on level-ground; (ii) bilateral squatting to max depth with arms crossed over chest; (iii) sit-to-stand from chair (44cm) with arms folded over chest and without use of arm rests; (iv) stair ascent; and (v) stair descent. Participants were given time to familiarise with each task, with all movements completed at a self-selected pace. Each participant was asked to perform 3 practice repetitions followed by a short break (1 to 2min) and 5 recorded repetitions of each task, with breaks (10 to 20min) between tasks to reduce fatigue. To allow for automatic calculation of joint kinematics, the Notch sensor system was calibrated according to manufacturer instructions by instructing the participants to stand stationary in a neutral reference posture, prior to each activity. For scaling of the OpenSim model, marker positions were collected by the OMC system from the participant in a static relaxed standing trial.

2.3 Experimental setup

A 10-camera Vicon OMC system (Vicon, Oxford, UK) was used to collect marker trajectories (200Hz) from 17 retro-reflective markers placed on the pelvis, thigh, shank, and foot of the affected limb (Fig 1). The marker set consisted of a three-marker cluster on the thigh and shank segments and anatomical markers placed over the right and left anterior-superior iliac spines; right and left posterior-superior iliac spines; medial and lateral femoral epicondyles; medial and lateral malleoli; posterior aspect of the foot and the first and fifth metatarsal heads. Ground reaction forces (1000Hz) were recorded for level overground-ground walking, squatting, and sit-to-stand activities using 3 force plates (510 mm × 465 mm, AMTI, USA) embedded into the ground and arranged in series in the centre of a walkway, approximately 10-metres long. For stair-ascent and stair-descent tasks, ground reaction forces (1000Hz) were recorded using a force-instrumented staircase with 3 force-plate embedded within the first three stairs (Bertec Corporation, Columbus, OH).

thumbnail
Fig 1. Participant preparation.

A. Schematic showing the placement IMUs on body segments based according to manufacturer recommendations B. Images illustrate the methodological setup used to simultaneously measure knee kinematics from inertial measurement units (IMUs) placed on the thigh, shank, and foot and an optical motion capture (OMC) system.

https://doi.org/10.1371/journal.pone.0297899.g001

Segment rotations and orientations (200Hz) were recorded simultaneously to the OMC and ground reaction forces using IMUs, which contained an accelerometer, a gyroscope, and a magnetometer (Notch Interfaces Inc. New York, USA). The IMUs were affixed to the pelvis, thigh, shank, and foot of the affected limb using double sided tape and Velcro straps, following manufacturer instructions with devices placed on the anterior portion of each segment oriented upwards (Fig 1). IMU data were collected simultaneously with the OMC system.

2.4 Data processing

Knee joint angles and ROM were calculated for each task using three methods: (i) OMC driven inverse kinematics in Opensim (OMC); (ii) direct calculation of relative angles from raw quaternion outputs of IMU (r_IMU); and by (iii) a proprietary IMU-driven kinematic model (Wearnotch, Notch Interfaces,Inc., NJ, USA) (m_IMU) (Fig 2). Ground reaction forces were used solely to identify gait events as well as the beginning and end of the squat and sit-to-stand movements. Marker trajectories and ground reaction forces were filtered using a second-order low-pass Butterworth filter with a cut-off frequency of 6 Hz prior to model scaling and inverse kinematics analysis. Marker positions from each participant’s static relaxed standing trials were used to linearly scale a generic 3-segment (thigh, shank, foot) musculoskeletal model with 3 degrees of freedom (sagittal, frontal, and transverse rotations) at both the knee and ankle in OpenSim v4.1 [34, 35]. Scaled models for each participant were subsequently used for all inverse kinematics in the OMC analysis.

thumbnail
Fig 2. Schematic of the data processing pipeline.

Marker trajectories from optical motion capture and accelerations and rotations captured from the Notch IMU system (m_IMU) were recorded simultaneously during 5 activities of daily living. Tri-planar knee angles of were calculated using inverse kinematics analysis in OpenSim using marker trajectories (OMC) and from direct calculation from quaternions (r_IMU). Only sagittal plane knee angles were calculated for IMU system with embedded kinematic model (m_IMU) as depicted in the bottom right-hand plot.

https://doi.org/10.1371/journal.pone.0297899.g002

For the r_IMU method, the three-dimensional relative motions between the thigh and shank were computed using the raw quaternions exported from the IMUs on-board sensor fusion algorithm. These quaternions were registered to the model during the calibration pose to determine the sensor to segment-to-segment alignment. Quaternions were converted to rotation matrixes, which were used to calculate the change in angle between the shank and thigh segments, prior to conversion to Euler angles using a XYZ sequence. For the OMC, marker trajectories were used in OpenSim inverse kinematics to compute flexion-extension, adduction-abduction, and internal-external rotation of the knee. For r_IMU method, joint angles were calculated as the rotation from thigh orientation to the shank orientation. The Notch IMU system automatically filtered and optimised raw data to calculate segment orientations. The Notch system computed sagittal knee joint angles only, using the measured orientations from the thigh and shank IMUs and Notch proprietary software. Thus, only the r_IMU method was analysed and compared to OMC for frontal and transverse plane motion.

A single representative trial (clear heel-strike/toe-off on force plates, no use of safety rail, clear start and end to squat and sit-to-stand movement, and no marker occlusion) for each participant was analysed for each task. Kinematic data from each method were time aligned using cross-correlation and spline interpolated to 101 time-points. To reduce errors associated with incorrect posture during calibration (i.e., standing with flexed knee), the relative angles were compared by subtracting the mean angle of each trial from all values, for each system. ROM was calculated as the maximum angle minus the minimum angle for each movement cycle.

2.5 Statistical analysis

Knee angles across the movement cycle and ROM using r_IMU and m_IMU methods were compared to OMC using Pearson’s correlation coefficient (R) and RMSD. R scores range between 0 and 1 and were stratified as: poor similarity (0 to 0.60); moderate (0.60 to 0.75); good (0.75 to 0.85); very good (0.85 to 0.95); and excellent (0.95 to 1) [36]. Sagittal plane knee ROM was calculated for each system and compared using a repeated measures ANOVA (P < 0.05). For significant differences, post-hoc paired t-tests with Bonferroni correction were performed. Pearson’s correlation coefficients and RMSD were computed in MATLAB v2022a, whilst the repeated measures ANOVA was conducted using SPSS v28.

Results

3.1 Correlation and RMSD–sagittal plane

In the sagittal plane the m_IMU and r_IMU methods estimated knee angles across the entire movement with excellent correlation (>0.95) compared to OMC method for walking, squatting, sit-to-stand, and stair-ascent, and very good correlation (>0.90) for stair descent (Table 1, Fig 3). The task with the highest correlation and lowest RMSD for the m_IMU was squatting (R = 1.00±0.01, RMSD = 3.42±1.53), followed by sit-to-stand (R = 0.98±0.04, RMSD = 4.19±1.91), walking (R = 0.93±0.09, RMSD = 6.16±2.94), stair ascent (R = 0.88±0.11, RMSD = 11.95±5.30), and stair descent (R = 0.88±0.13, RMSD = 11.81±8.38). Similarly, the r_IMU method demonstrated highest correlation and lowest RMSD for squatting (R = 0.99±0.01, RMSD = 3.94±2.24), followed by sit-to-stand (R = 0.98±0.04, RMSD = 4.57±1.90), walking (R = 0.93±0.08, RMSD = 6.19±2.90), stair ascent (R = 0.88±0.11, RMSD = 12.06±5.21), and stair descent (R = 0.88±0.13, RMSD = 11.81±8.38).

thumbnail
Fig 3. Comparison of IMU measured knee kinematics to optical motion capture.

Group mean (± 95% confidence interval) time-series normalised knee kinematics during walking (n = 15), squatting (n = 17), sit-to-stand (n = 17), stair ascent (n = 12), and stair descent (n = 11) from OMC (grey dashed lines), r_IMU (green solid lines), and the m_IMU (blue solid lines) methods. Positive values indicate flexion, adduction, and internal rotation.

https://doi.org/10.1371/journal.pone.0297899.g003

thumbnail
Table 1. Group mean (± standard deviation) Pearson’s correlation coefficient (R) and root-mean-square-difference (RMSD, in degrees) of the m_IMU and r_IMU calculated knee kinematics compared to the OMC.

https://doi.org/10.1371/journal.pone.0297899.t001

3.2 Correlation and RMSD–frontal and transverse planes

In the frontal plane, the r_IMU estimated knee kinematics with poor (R < 0.60 for walking, sit-to-stand, stair ascent, and stair descent) to moderate (R = 0.6 to 0.75 for squatting) correlation with OMC values (Table 1, Fig 3). The highest R (0.69±0.23) and lowest RMSD (5.89±63.80) for the frontal plane were found in the squatting and sit-to-stand tasks, respectively. In the transverse plane, the r_IMU estimated knee kinematics with poor (R < 0.60 for walking, stair ascent, and stair descent), moderate (R = 0.6 to 0.75 for sit-to-stand) and good (R = 0.75 to 0.85 for squatting) correlation with OMC values (Table 1, Fig 3). For the transverse plane, the highest R (0.78±0.15) was found in the squatting task whilst the lowest RMSD (5.19±2.94) was found in walking. The errors observed in the frontal and transverse planes are high considering the low ROM observed in these planes (Fig 3).

3.3 Range of motion–sagittal plane

Between the three methods, there was no significant difference in sagittal plane ROM for walking, squatting, and sit-to-stand tasks (Table 2). For stair ascent, a repeated measures ANOVA found a significant difference between the OMC (88.56±9.72), m_IMU (84.93±12.19), and r_IMU (82.42±11.70) (F = 5.37, p = 0.03). For stair descent, there was a significant difference between OMC (88.90±9.63), m_IMU (83.46±9.76), and r_IMU (81.72±9.35) (F = 4.95, p = 0.04). However, post-hoc comparison with Bonferroni correction (*p<0.017) showed no significant difference between methods for both stair ascent and stair descent.

thumbnail
Table 2. Group mean (± standard deviation) knee sagittal plane range of motion (ROM) and mean difference between methods (95% confidence interval) calculated using optical motion capture (OMC), model IMU (m_IMU) and raw IMU (r_IMU) data. r_IMU.

https://doi.org/10.1371/journal.pone.0297899.t002

Discussion

Inertial measurement units provide an easy-to-use and objective tool to assess knee function and recovery following TKA. In participants 1-year post robotically assisted TKA, we found IMU methods (m_IMU and r_IMU) could measure sagittal plane knee angles with very good to excellent correlation across five activities of daily living (walking, squatting, sit-to-stand, stair ascent, and stair descent) compared to OMC, in agreement with our first hypothesis. For example, walking showed excellent correlation and mean RMSD of less than 6.5 degrees compared to OMC. The highest correlation and lowest RMSD were observed during squatting and sit-to-stand, suggesting that IMU-based estimations better match OMC-based estimations during closed-chain, low velocity movements. Frontal and transverse plane knee kinematics calculated using IMU methods, should however be interpreted with caution, as they typically demonstrated lower correlations and higher RMSD relative to their total ROM compared to OMC across all tasks. Differences between the three methods were identified in ROM for stair ascent and stair descent tasks only, although pairwise comparison failed to identify significant difference between comparisons, suggesting differences are not meaningful. Therefore, our findings support the use of IMU as a clinical tool for measurement of sagittal plane knee kinematics and ROM during activities of daily living in individuals following TKA. The implementation of functional assessments should, however, consider tasks with excellent correlation and low error compared to gold standard measurement such as squatting, sit-to-stand, and walking. IMU offer an easy-to-use, rapid, and cost-effective method for sagittal plane knee kinematics analysis in the clinic, outdoors, or in the home where OMC is impractical.

While minimum differences were found between the r_IMU and m_IMU methods, an IMU system with an underlying kinematic model (m_IMU) may be better suited to clinical applications. IMU systems with an underlying kinematic model and user-friendly software simplify data acquisition (~ 5 mins) and overcome many of the challenges of recording biomechanical data outside a research laboratory [37]. In addition, the real-time kinematic model provides instant feedback on knee function that can be used to assess improvements, identify tasks with suboptimal biomechanics, and guide rehabilitation exercises. Knee flexion ROM, an important clinical measure indicative of recovery following TKA [38], is typically measured passively or actively using goniometers [39]. Assessment during functional tasks, like walking, offers clinically relevant insight of knee function following TKA [40, 41]. A potential application is the use of IMU at scale to elucidate differences in functional outcomes due to surgical techniques, surgery duration, implant positioning, and type of implant, and may be enhanced after robotically assisted TKA during which surgical data may be automatically logged.

These results are consistent with previous studies that estimated sagittal plane knee angles during walking with RMSD between 4–7 degrees [30, 42, 43]. The differences observed may be attributed to several factors including gyroscopic drift for IMUs [43], a non-zero knee angle during the calibration pose, and movement of the soft tissue to which the sensors and markers are fixed [44]. The magnitude of soft-tissue artefact is affected by the location of markers/sensor placement, with studies suggesting the thigh is particularly susceptible to translational and rotational errors [45, 46]. Soft tissue artefact is also dependent on individual characteristics [46] and may be greater in participants with high levels of adipose tissue. Additionally, soft-tissue artefact may affect IMU and OMC methods differently, with a recent study showing a IMU method over-estimated sagittal plane knee angles during stance, and underestimated them during swing [47]. Characterising the differing effects of soft-tissue artefact on data collection systems will allow for the development of correction methods that will improve IMU-based predictions. Finally, although OMC is considered the gold-standard for non-invasive human motion capture, error in the OMC method should also be considered as a source of difference, with a comparison of five marker-based OMC protocols showing a mean absolute variance of greater than 3 degrees between protocols [47]. Additionally, the mean RMSD in the flexion-extension axis for both the m_IMU and r_IMU methods when compared to OMC were within the reported uncertainty for OMC-based kinematic measurements for all tasks except stair ascent and stair descent [48].

There are limitations in our study that should be considered. The r_IMU method measured frontal and transverse plane kinematics with poor to moderate correlation to the OMC method. The disagreement in the frontal and transverse planes may be attributed to sensitivity to marker placement protocols [47], the high ratio of soft-tissue movement to joint motion in the frontal and transverse planes, and the definition of the flexion-extension axis [49]. Errors in the flexion-extension axis alignment during calibration can lead to angular crosstalk, where motion is incorrectly represented in the transverse and frontal planes. While our results support the use of IMU for the measurement of sagittal plane knee kinematics across a range of ADL, we did not assess reliability during long recordings, for which orientation estimates may drift, and error can increase [42]. Additionally, we used embedded force plates (ground and instrumented staircase) to identify gait events which may affect IMU orientations due to electromagnetic interference [50]. Greater angular crosstalk due to more time spent in knee flexion compared to walking, and greater magnetic electromagnetic disturbance of the largely exposed instrumented staircase likely contributed to poorer results in the stair ascent and descent tasks. Further, several participants were excluded from analysis due to marker occlusion, particularly in stair-related tasks caused by the required safety railings. There may be differences in r_IMU and m_IMU range of motion compared to OMC that were not detected due to our small sample size. Nevertheless, we believe the high correlation and low RMSD in the flexion-extension axis supports the use of IMUs for measurement of sagittal plane knee kinematics during ADL. To ensure participants could complete multiple tasks in a single session, participant inclusion criteria included the ability to walk for 30 minutes continuously and unaided. Excluding participants who could not walk for 30 minutes continuously likely biased our study to individuals with good TKA outcomes who may have different knee kinematics to those with poor TKA outcomes. Future studies may build on our findings by assessing the reliability of IMU-based measurements of kinematics and ROM during ADL in the real-world (e.g., walking outside, functional tasks in the home), assess the ability for IMU to detect clinically meaningful changes in knee kinematics, and prospectively evaluate changes in knee kinematics post TKA over the long term. Additionally, future studies may improve results in the frontal and transverse planes by utilising functional calibration tasks and using biplanar radiography to assess the frontal and transverse accuracy of IMU measured knee kinematics [51]. Estimation of knee kinematics may also be improved by combining IMUs with marker-less motion capture [52], body-worn cameras [50] or machine learning techniques [53, 54], which have shown promising results in predicting human movement kinematics.

Conclusions

Our results demonstrate that IMU may accurately quantify sagittal plane knee motion across various ADL that may complement subjective patient reported outcome measures throughout the episode of care for TKA patients. Both IMU methods have excellent agreement with gold-standard OMC methods for sagittal plane knee kinematics during squatting, sit-to-stand, and walking. Frontal and transverse plane kinematics measured with IMU showed poor-to-moderate correlation with OMC and should be interpreted with caution. Objective measures of knee kinematics may provide insight into the effectiveness of surgical techniques and can be measured easily and cost-effectively using IMU. Future studies should determine if knee kinematics from IMU measurements can discriminate between excellent, moderate, and poor outcomes, and to what extent these differences are associated with surgical techniques.

Supporting information

S1 Checklist. STROBE statement—checklist of items that should be included in reports of observational studies.

https://doi.org/10.1371/journal.pone.0297899.s001

(DOCX)

S1 Dataset. Time-series knee angles measured by OMC and IMU.

https://doi.org/10.1371/journal.pone.0297899.s002

(CSV)

Acknowledgments

The authors thank all participants for their time and contribution. The authors also thank Erin Bickell for her assistance with patient recruitment and data collection for this study.

References

  1. 1. Batailler C, Fernandez A, Swan J, Servien E, Haddad FS, Catani F, et al. MAKO CT-based robotic arm-assisted system is a reliable procedure for total knee arthroplasty: a systematic review. Knee Surg Sports Traumatol Arthrosc. 2021;29(11):3585–98. pmid:32975626
  2. 2. Kim MS, Koh IJ, Choi YJ, Lee JY, In Y. Differences in Patient-Reported Outcomes Between Unicompartmental and Total Knee Arthroplasties: A Propensity Score-Matched Analysis. J Arthroplasty. 2017;32(5):1453–9. pmid:27979407
  3. 3. Kahlenberg CA, Nwachukwu BU, McLawhorn AS, Cross MB, Cornell CN, Padgett DE. Patient Satisfaction after Total Knee Replacement: A Systematic Review. HSS Journal®. 2018;14(2):192–201. pmid:29983663
  4. 4. Matsuda S, Kawahara S, Okazaki K, Tashiro Y, Iwamoto Y. Postoperative Alignment and ROM Affect Patient Satisfaction After TKA. Clin Orthop Relat Res. 2013;471(1). pmid:22903282
  5. 5. Bell SW, Anthony I, Jones B, MacLean A, Rowe P, Blyth M. Improved Accuracy of Component Positioning with Robotic-Assisted Unicompartmental Knee Arthroplasty: Data from a Prospective, Randomized Controlled Study. J Bone Joint Surg. 2016;98(8). pmid:27098321
  6. 6. Elliott J, Shatrov J, Fritsch B, Parker D. Robotic-assisted knee arthroplasty: an evolution in progress. A concise review of the available systems and the data supporting them. Arch Orthop Trauma Surg. 2021;141(12):2099–117. pmid:34491411
  7. 7. Agarwal N, To K, McDonnell S, Khan W. Clinical and Radiological Outcomes in Robotic-Assisted Total Knee Arthroplasty: A Systematic Review and Meta-Analysis. J Arthroplasty. 2020;35(11):3393–409.e2. pmid:32234326
  8. 8. Cool CL, Jacofsky DJ, Seeger KA, Sodhi N, Mont MA. A 90-day episode-of-care cost analysis of robotic-arm assisted total knee arthroplasty. J Comp Eff Res 2019;8(5):327–36. Epub 2019/01/29. pmid:30686022.
  9. 9. Mizner RL, Petterson SC, Clements KE, Zeni JA Jr., Irrgang JJ, Snyder-Mackler L. Measuring functional improvement after total knee arthroplasty requires both performance-based and patient-report assessments: a longitudinal analysis of outcomes. J Arthroplasty. 2011;26(5):728–37. Epub 2010/09/20. pmid:20851566.
  10. 10. Jacobs CA, Christensen CP. Correlations between knee society function scores and functional force measures. Clin Orthop Relat Res. 2009;467(9):2414–9. Epub 2009/04/03. pmid:19343460.
  11. 11. Eckhard L, Munir S, Wood D, Talbot S, Brighton R, Walter B, et al. The ceiling effects of patient reported outcome measures for total knee arthroplasty. Orthop Traumatol-Sur. 2021;107(3):102758. https://doi.org/10.1016/j.otsr.2020.102758.
  12. 12. Childs JD, Sparto PJ, Fitzgerald GK, Bizzini M, Irrgang JJ. Alterations in lower extremity movement and muscle activation patterns in individuals with knee osteoarthritis. Clin Biomech. 2004;19(1):44–9. pmid:14659929
  13. 13. SchmitT LC, Rudolph KS. Influences on knee movement strategies during walking in persons with medial knee osteoarthritis. Arthritis Care Res. 2007;57(6):1018–26. pmid:17665469
  14. 14. McClelland JA, Webster KE, Feller JA. Gait analysis of patients following total knee replacement: a systematic review. Knee. 2007;14(4):253–63. pmid:17531493
  15. 15. Sosdian L, Dobson F, Wrigley TV, Paterson K, Bennell K, Dowsey M, et al. Longitudinal changes in knee kinematics and moments following knee arthroplasty: a systematic review. Knee. 2014;21(6):994–1008. Epub 20141012. pmid:25311517.
  16. 16. Naili JE, Iversen MD, Esbjörnsson A-C, Hedström M, Schwartz MH, Häger CK, et al. Deficits in functional performance and gait one year after total knee arthroplasty despite improved self-reported function. Knee Surg Sports Traumatol Arthrosc. 2017;25(11):3378–86. pmid:27436196
  17. 17. Levinger P, Menz HB, Morrow AD, Feller JA, Bartlett JR, Bergman NR. Lower Limb Biomechanics in Individuals With Knee Osteoarthritis Before and After Total Knee Arthroplasty Surgery. J Arthroplasty. 2013;28(6):994–9. pmid:23523504
  18. 18. Risitano S, Cacciola G, Capella M, Bosco F, Giustra F, Fusini F, et al. Comparison between gaits after a medial pivot and posterior stabilized primary total knee arthroplasty: a systematic review of the literature. Arthroplasty. 2023;5(1):15. pmid:36927464
  19. 19. Christensen JC, Capin JJ, Hinrichs LA, Aljehani M, Stevens-Lapsley JE, Zeni JA. Gait mechanics are influenced by quadriceps strength, age, and sex after total knee arthroplasty. J Orthop Res. 2021;39(7):1523–32. pmid:33034899
  20. 20. Lloyd D. The future of in-field sports biomechanics: wearables plus modelling compute real-time in vivo tissue loading to prevent and repair musculoskeletal injuries. Sports Biomech. 2021:1–29. pmid:34496728
  21. 21. Mundt M, Johnson WR, Potthast W, Markert B, Mian A, Alderson J. A Comparison of Three Neural Network Approaches for Estimating Joint Angles and Moments from Inertial Measurement Units. Sensors. 2021;21(13). pmid:34283080
  22. 22. Robert-Lachaine X, Mecheri H, Larue C, Plamondon A. Validation of inertial measurement units with an optoelectronic system for whole-body motion analysis. Med Biol Eng Comput. 2017;55(4):609–19. pmid:27379397
  23. 23. 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 (Basel). 2020;20(3). Epub 2020/01/30. pmid:31991862; PubMed Central PMCID: PMC7038336.
  24. 24. Diamond LE, Grant T, Uhlrich SD. Osteoarthritis year in review 2023: biomechanics. Osteoarthr Cartil. 2023. pmid:38043858
  25. 25. Adesida Y, Papi E, McGregor AH. Exploring the Role of Wearable Technology in Sport Kinematics and Kinetics: A Systematic Review. Sensors. 2019;19(7). pmid:30987014
  26. 26. Kaufmann M, Nüesch C, Clauss M, Pagenstert G, Eckardt A, Ilchmann T, et al. Functional assessment of total hip arthroplasty using inertial measurement units: Improvement in gait kinematics and association with patient-reported outcome measures. J Orthop Res. n/a(n/a). https://doi.org/10.1002/jor.25421.
  27. 27. Follis S, Chen Z, Mishra S, Howe CL, Toosizadeh N, Dohm M. Comparison of wearable sensor to traditional methods in functional outcome measures: A systematic review. J Orthop Res. 2021;39(10):2093–102. pmid:33300119
  28. 28. Teufl W, Miezal M, Taetz B, Fröhlich M, Bleser G. Validity of inertial sensor based 3D joint kinematics of static and dynamic sport and physiotherapy specific movements. PLOS ONE. 2019;14(2):e0213064. pmid:30817787
  29. 29. Seel T, Raisch J, Schauer T. IMU-based joint angle measurement for gait analysis. Sensors (Basel, Switzerland). 2014;14(4):6891–909. pmid:24743160.
  30. 30. Niswander W, Wang W, Kontson K. Optimization of IMU Sensor Placement for the Measurement of Lower Limb Joint Kinematics. Sensors (Basel, Switzerland). 2020;20(21):5993. pmid:33105876.
  31. 31. Zügner R, Tranberg R, Timperley J, Hodgins D, Mohaddes M, Kärrholm J. Validation of inertial measurement units with optical tracking system in patients operated with Total hip arthroplasty. BMC Musculoskelet Disord. 2019;20(1):52. pmid:30727979
  32. 32. Al-Amri M, Nicholas K, Button K, Sparkes V, Sheeran L, Davies JL. Inertial Measurement Units for Clinical Movement Analysis: Reliability and Concurrent Validity. Sensors (Basel, Switzerland). 2018;18(3):719. pmid:29495600.
  33. 33. Antunes R, Jacob P, Meyer A, Conditt MA, Roche MW, Verstraete MA. Accuracy of Measuring Knee Flexion after TKA through Wearable IMU Sensors. J Funct Morphol Kinesiol. 2021;6(3). Epub 2021/07/22. pmid:34287303; PubMed Central PMCID: PMC8293382.
  34. 34. Seth A, Hicks JL, Uchida TK, Habib A, Dembia CL, Dunne JJ, et al. OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. PLOS Computl Biol. 2018;14(7):e1006223. pmid:30048444
  35. 35. Delp SL, Anderson FC, Arnold AS, Loan P, Habib A, John CT, et al. OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans Biomed Eng. 2007;54(11):1940–50. Epub 2007/11/21. pmid:18018689.
  36. 36. Di Marco R, Scalona E, Pacilli A, Cappa P, Mazzà C, Rossi S. How to choose and interpret similarity indices to quantify the variability in gait joint kinematics. Int Biomech. 2018;5(1):1–8.
  37. 37. Rana M, Mittal V. Wearable Sensors for Real-Time Kinematics Analysis in Sports: A Review. IEEE Sensors Journal. 2021;21(2):1187–207.
  38. 38. Oka T, Wada O, Asai T, Maruno H, Mizuno K. Importance of knee flexion range of motion during the acute phase after total knee arthroplasty. Phys Ther Res. 2020;23(2):143–8. Epub 2021/01/26. pmid:33489651; PubMed Central PMCID: PMC7814201.
  39. 39. Lavernia C D’Apuzzo M, Rossi MD, Lee . Accuracy of Knee Range of Motion Assessment After Total Knee Arthroplasty. J Arthroplasty. 2008;23(6, Supplement):85–91. pmid:18722308
  40. 40. McClelland JA, Webster KE, Feller JA, Menz HB. Knee kinematics during walking at different speeds in people who have undergone total knee replacement. The Knee. 2011;18(3):151–5. pmid:20537900
  41. 41. Yoshida Y, Mizner RL, Ramsey DK, Snyder-Mackler L. Examining outcomes from total knee arthroplasty and the relationship between quadriceps strength and knee function over time. Clin Biomech. 2008;23(3):320–8. pmid:18060669
  42. 42. Al Borno M O’Day J, Ibarra V, Dunne J, Seth A, Habib A, et al. OpenSense: An open-source toolbox for inertial-measurement-unit-based measurement of lower extremity kinematics over long durations. J Neuroeng Rehabilitation. 2022;19(1):22. pmid:35184727
  43. 43. Bailey CA, Uchida TK, Nantel J, Graham RB. Validity and Sensitivity of an Inertial Measurement Unit-Driven Biomechanical Model of Motor Variability for Gait. Sensors (Basel, Switzerland). 2021;21(22):7690. pmid:34833766.
  44. 44. Johnson CD, Outerleys J, Tenforde AS, Davis IS. A comparison of attachment methods of skin mounted inertial measurement units on tibial accelerations. J Biomech. 2020;113:110118. pmid:33197691
  45. 45. Garling EH, Kaptein BL, Mertens B, Barendregt W, Veeger HEJ, Nelissen RGHH, et al. Soft-tissue artefact assessment during step-up using fluoroscopy and skin-mounted markers. J Biomech. 2007;40:S18–S24. pmid:17462655
  46. 46. Peters A, Galna B, Sangeux M, Morris M, Baker R. Quantification of soft tissue artifact in lower limb human motion analysis: A systematic review. Gait Posture. 2010;31(1):1–8. pmid:19853455
  47. 47. Ferrari A, Benedetti MG, Pavan E, Frigo C, Bettinelli D, Rabuffetti M, et al. Quantitative comparison of five current protocols in gait analysis. Gait Posture. 2008;28(2):207–16. pmid:18206374
  48. 48. Cappozzo A, Catani F, Leardini A, Benedetti MG, Croce UD. Position and orientation in space of bones during movement: experimental artefacts. Clin Biomech (Bristol, Avon). 1996;11(2):90–100. Epub 1996/03/01. pmid:11415604.
  49. 49. Schache AG, Baker R, Lamoreux LW. Defining the knee joint flexion–extension axis for purposes of quantitative gait analysis: An evaluation of methods. Gait Posture. 2006;24(1):100–9. pmid:16191481
  50. 50. Drobnjakovic F, Douangpaseuth JB, Gadea C, Haider M, Ionescu D, Ionescu B, et al., editors. Fusing data from inertial measurement units and a 3D camera for body tracking. 2018 IEEE International Instrumentation and Measurement Technology Conference (I2MTC); 2018 14–17 May 2018.
  51. 51. Miranda DL, Rainbow MJ, Crisco JJ, Fleming BC. Kinematic differences between optical motion capture and biplanar videoradiography during a jump-cut maneuver. J Biomech. 2013;46(3):567–73. Epub 20121022. pmid:23084785; PubMed Central PMCID: PMC3551998.
  52. 52. Halilaj E, Shin S, Rapp E, Xiang D. American society of biomechanics early career achievement award 2020: Toward portable and modular biomechanics labs: How video and IMU fusion will change gait analysis. J Biomech. 2021;129:110650. pmid:34644610
  53. 53. Rapp E, Shin S, Thomsen W, Ferber R, Halilaj E. Estimation of kinematics from inertial measurement units using a combined deep learning and optimization framework. J Biomech. 2021;116:110229. pmid:33485143
  54. 54. Burton WS, Myers CA, Rullkoetter PJ. Machine learning for rapid estimation of lower extremity muscle and joint loading during activities of daily living. J Biomech. 2021;123:110439. pmid:34004394