Next Article in Journal
Application of Moire Profilometry in Three-Dimensional Profile Reconstruction of Key Parts in Railway
Next Article in Special Issue
Human Movement Quality Assessment Using Sensor Technologies in Recreational and Professional Sports: A Scoping Review
Previous Article in Journal
Industry 4.0 Technologies Applied to the Rail Transportation Industry: A Systematic Review
Previous Article in Special Issue
Neurofeedback Therapy for Sensory Over-Responsiveness—A Feasibility Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Accurate Prediction of Knee Angles during Open-Chain Rehabilitation Exercises Using a Wearable Array of Nanocomposite Stretch Sensors

1
Department of Mechanical Engineering, Brigham Young University, Provo, UT 84602, USA
2
Department of Exercise Science, Brigham Young University, Provo, UT 84602, USA
3
Department of Statistics, Brigham Young University, Provo, UT 84602, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2022, 22(7), 2499; https://doi.org/10.3390/s22072499
Submission received: 24 February 2022 / Revised: 20 March 2022 / Accepted: 21 March 2022 / Published: 24 March 2022
(This article belongs to the Special Issue Biomedical Sensing for Human Motion Monitoring)

Abstract

:
In this work, a knee sleeve is presented for application in physical therapy applications relating to knee rehabilitation. The device is instrumented with sixteen piezoresistive sensors to measure knee angles during exercise, and can support at-home rehabilitation methods. The development of the device is presented. Testing was performed on eighteen subjects, and knee angles were predicted using a machine learning regressor. Subject-specific and device-specific models are analyzed and presented. Subject-specific models average root mean square errors of 7.6 and 1.8 degrees for flexion/extension and internal/external rotation, respectively. Device-specific models average root mean square errors of 12.6 and 3.5 degrees for flexion/extension and internal/external rotation, respectively. The device presented in this work proved to be a repeatable, reusable, low-cost device that can adequately model the knee’s flexion/extension and internal/external rotation angles for rehabilitation purposes.

1. Introduction

The rate of invasive knee surgeries is rising. The number of total knee arthroplasty (TKA) procedures has steadily increased over the past three decades [1,2]. Future projections estimate that the number will increase to 1.16–3.48 million annually by 2030 [2,3], resulting in about 90 procedures per orthopedic surgeon per year in the United States [4]. ACL reconstructions are also expected to increase [5]. There is a proportional need for post-operative physical therapy. However, recent estimates show a decreasing ratio of physical therapists (PTs) per unit of population [6,7], partly due to burnout and career dissatisfaction [8] among PTs. The imbalance between the rising number of patients requiring rehabilitation and a decreasing workforce is expected to grow.
Home-based rehabilitation methods have been developed to accommodate recent trends relating to available PTs and the number of patients requiring therapy. Compared to inpatient or outpatient rehabilitation, home-based interventions are less expensive [9,10,11,12] and more accessible [12,13]. One critique of at-home rehabilitation, however, is decreased regimen adherence that is usually seen in unguided or unsupervised rehabilitation [14]. However, important advantages of at-home rehabilitation, as described by patients, include ease of use of many currently available systems, decreased stress, and “not having to travel when we are in pain” [13].
One solution to increase rehabilitation availability and patient accountability is sensor-based technology. Such technology could potentially maintain quality of rehabilitative care despite increasing patient-to-PT ratios. Sensor-based technology could also give the PT quantitative feedback on the patient’s progress while performing at-home exercises. However, sensor-based rehabilitation also has drawbacks: battery limitations, sensor drift over time [15], device cost, and mobility. Furthermore, depending on the specific approach, and its implementation, limitations might outweigh benefit and restrict widespread use of sensor technology in physical rehabilitation. For example, if a sensor-based device exhibits a large amount of drift, the error contained within the data may be large enough to render the system unusable. This work presents a cost-efficient, low-power device to facilitate passively supervised at-home rehabilitation via a knee sleeve instrumented with an array of piezoresistive sensors.
High-deflection strain gauges are desirable for biomechanical measurement applications because of their large strain capacity and sensitivity to changes in strain [16,17]. There are many types of high-deflection strain gauges used in biomechanics including liquid metal [18,19,20], gel-based [21,22], polymer optical fiber [23,24,25], and piezoresistive sensors. The most common filler materials used in high-deflection strain gauges are derivatives of carbon [26,27,28,29,30], due to the favorable mechanical properties, high electrical conductivity, and high thermal conductivity of carbon-based polymers.
A notable biomechanical application of polymer-based high-deflection strain gauges is smart textiles, garments ranging from braces to shirts that are instrumented with sensors to measure forces, positions, or pressure in the body [31]. For example, Shyr et al. [32] estimated knee and elbow angles using a woven elastic textile. Gholami et al. [33] developed a pair of leg tights capable of measuring knee flexion and extension with a high degree of repeatability via machine learning regressors. These and other examples demonstrate the potential for wearable sensors in TKA rehabilitation.
Although knee flexion and extension are essential for successful TKA rehabilitation, the knee should not be considered a single degree of freedom joint. It has been well established that transverse plane (internal/external of the tibia) rotation occurs naturally in gait [34] and is a conjunct (or linked) movement necessary for full range of motion [35]. All mentions of internal/external in this work are concerning the tibia unless otherwise noted. Post-operative tibia rotation correlates with post-operative flexion/extension [36], but may not necessarily be correctly rehabilitated due to prosthesis misalignment [37] or alterations in the patient’s gait patterns [36]. Recovery of normal ranges of rotation of the tibia is especially imperative for flexion while bearing weight [38]. Inadequate internal rotation of the tibia is a common cause of chronic pain after TKA [37]. Thus, for a complete recovery of the knee’s pre-operative range of motion with decreased pain, flexion/extension and internal/external rotation should be simultaneously monitored and improved during rehabilitation.
This work utilizes a triphasic silicone/nickel composite [39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55] with two different nickel filler materials—nickel nanostrands (NiNs, (18.76 wt%); Figure 1a) and nickel-coated carbon fiber (NCCF, (3.13 wt%), Figure 1b). The resulting sensors are inversely piezoresistive (negative correlation between strain and electrical resistance). A single triphasic silicone/nickel sensor was previously used to measure knee flexion/extension [54], but an array of these sensors can measure additional degrees of freedom in the knee. We hypothesize that, by placing piezoresistive strain gauges in a garment across locations of the anterior knee, where the skin experiences the greatest variance in strain, we can accurately model and predict two rotational degrees of freedom: flexion/extension and tibia/femur rotation, during open-chain knee flexion.
The device presented in this work was used to measure two rotational degrees of freedom (flexion/extension and internal/external rotation) of the knee in healthy subjects to determine the accuracy and repeatability of a 16-sensor, skin-mounted instrumented sleeve. Future application is intended to occur in physical therapy applications during TKA rehabilitation. The technology represents a cost-efficient instrument that can eventually be used to enhance at-home rehabilitation for TKA patients by providing a quantifiable measurement of two rotational degrees of freedom. Future work will incorporate these measurements into a real-time biofeedback system for improving the quality of at-home physical therapy exercises [56].

2. Materials and Methods

2.1. Data Collection

The sensor system for this work, shown in Figure 2, was mounted onto an off-the-shelf 88% copper nylon, 12% spandex knee sleeve fitted so as not to squeeze the knee joint, yet tight enough to not migrate during exercise. The positioning of an array of piezoresistive sensors was developed based on previous work [57]. The fitting of the device was devised with the guidance of an orthopedic surgeon and experienced PT. The instrumented knee sleeve conforms to the skin of the knee without adhering to the skin and can be directly donned and doffed by the user. A custom printed circuit board was developed to read and simultaneously log data from all sixteen sensors. The sensors were connected with a multiplexer into a Wheatstone bridge circuit that was excited with a 5-Volt square wave at 300 Hz. The printed circuit board was externally powered and logged the data to an accompanying computer through a serial connection.
The sensor component materials, nickel nanostrands (18.76 wt%), nickel-coated carbon fibers (3.13 wt%), and Ecoflex 00-30 silicone (Smooth-On Inc., Macungie, PA, USA), were mixed and the resultant slurry was pressed into aluminum molds of three different lengths. Five sensors were 21.0 mm long, ten were 31.2 mm, and one was 42.0 mm in length. All sensors were 0.8 mm thick and 5.0 mm wide. Copper wires were embedded into the ends of every sensor by placing 32-gauge wires within the slurry material before curing at 190 °F for 90 min. Then, a thin layer of Libra gloss silicone from Zodiac (Kennesaw, GA, USA) was spread onto a nylon/spandex knee sleeve over the desired location for each sensor and cured at 190 °F for 30 min. This silicone layer created a base under each sensor that prevented the Ecoflex silicone sensor adhesive from wicking into the nylon during the curing process. Finally, sixteen sensors were adhered to the silicone footings and cured at 190 °F for 30 min. To connect the sensors to the circuit board described previously, 36-gauge copper magnet wire was embroidered into a separate piece of rayon spandex cloth and secured to the sleeve with stitch points between the sensors. By ensuring the additional material fit more loosely than the knee sleeve, the layer of wires did not inhibit the motion of the knee sleeve on a subject’s skin. The approximate manufacturing cost of this sensing system was $7–8 plus electronics (roughly $20) and the cost of an off-the-shelf knee sleeve.
Previous work identified positions around the knee that are sensitive to knee motion and informed the layout of the sensors in the device presented here [57]. The sensor array was designed to capture redundant information in the array due to overlaps in sensor orientation and positioning. Because PT goals are oriented strictly around knee flexion/extension and internal/external rotation, the present work only included these movements, but future work may include knee translation, valgus/varus motion, and patella movement.
To develop a predictive model for the sensor array and validate the array performance, an optical tracking system was used in conjunction with the sensor array. Twenty one markers were placed on the right leg and pelvis to measure knee kinematics (Figure 3). Motion capture data were collected at 100 Hz, and sensor data were captured at 18.75 Hz/sensor. The same computer captured both sets of data so the two sets could be synchronized with respect to time using the computer’s local timestamp. Knee angles were calculated from the marker data using a generic lower-body model in Visual3D (C-Motion, Germantown, MD).
Eighteen healthy subjects (9 male, 9 female) between the ages of 20 and 47 (27.4 ± 7.1 years, mean ± SD), of varying body sizes (body mass indexes ranged from 18.3 to 33.9, with a mean of 23.9 ± 4.1) participated in this study, to determine the accuracy and repeatability of the multi-sensor system. Information about the cohort can be found in Table 1. Exclusion criteria for these subjects were past knee or hip surgeries or major non-surgical knee injuries within the past two years. Subjects performed open-chain knee flexion (OCKF), an exercise that replicates the range of motion of movements seen during post-operative TKA rehabilitation. OCKF was continually performed at an instructed pace of “two seconds up and two seconds down” for 60 s.

2.2. Modeling and Predictions

The collected data were used to produce subject-specific models and predict knee angles. A subject-specific model is defined in this work as a model that was trained on the data from a single subject and predicts knee angles of the same subject. Each model was trained and tested with ten-fold cross-validation from the data collected from one subject in the previous section. First, high-frequency noise was removed from the data using a Butterworth filter, with cutoff frequencies of 10 Hz for the marker motion data and 2 Hz for the strain sensor data. Kinematic experimental data are commonly filtered with a Butterworth filter with a cut-off frequency between 3 and 10 Hz [58]. The lower cut-off frequency of 2 Hz was used to filter the strain sensor data because the signal-to-noise ratio of the instrumented knee sleeve was higher than the motion capture data. A cut-off frequency that was double the frequency of the captured motion was used. Although it is lower than common cut-off frequencies, 2 Hz was experimentally found to provide the greatest increase in the strain sensor data signal-to-noise ratio. Once filtered, the motion capture data were downsampled to match the capture rate of the instrumented knee sleeve.
Unfiltered data from the sixteen sensors contained 1- and 2-point spikes that were more than triple the magnitude of the rest of the signal. Static electricity between the nylon–spandex knee sleeve material and exposed soldered points around the sensors likely caused these transient spikes. After these outlying spikes were eliminated from the sensor data using a Hampel filter (threshold = 2 SD), a Savitzky-Golay smoothing filter using a centered window smoothed the sensor data. At this point, 16 additional variables were created from the derivatives of each sensor’s data with time. Derivatives were calculated using finite differences between consecutive data points. Lastly, the mean values were subtracted from the sensor and derivative signals, and the data were scaled to unit variance before predicting knee angles. The model’s explanatory variables were the scaled sensor outputs and their derivatives. Figure 4 shows the post-processing pipeline.
Eighteen machine learning models were evaluated as candidates for mapping the composite sensor responses to the knee angles. These models include linear regression, neural networks, random forest regressors, adaptive boosted random forests, etc. All included models are listed in Table A1 and Table A2. The root mean square error (RMSE) and coefficient of determination (R2) were used to assess model accuracy. The model with the lowest mean RMSE over all subjects was selected, excluding models with excessively large outlier RMSEs. Removing models with outlying RMSEs decreased the likelihood of the chosen model failing for out-of-sample participants. To further test the validity of the instrumented knee sleeve, an analysis of bias and variability in predictions was performed to evaluate the predictive model’s agreement with the knee angles calculated from motion capture data.
Predictive modeling of data were completed using scikit-learn [59]. The results from ten separate models trained on different portions of a subject’s data were averaged to determine the robustness of the model’s predictive capability. Each model was trained on 90% of a subject’s data and tested on the remaining 10% of data. The training and testing data were separated in a way to preserve the time series component in the data. Multivariate machine learning regressor hyperparameters were optimized using the entire dataset from a subject using grid search and multivariate Bayesian optimization approaches. The average RMSE and R2 values from the ten unique models trained and tested on different portions of a subject’s data determined each model’s predictive power. Three instrumented knee sleeves were tested; therefore, both subject-specific and device-specific models were analyzed.

3. Results

3.1. Data Collection

Subjects performed exercises in similar ranges of motion with respect to knee angle; a complete record of knee angle ranges is found in Figure 5. Walking trials included knee flexion angles between 2.9 ± 3 . 9 and 58.7 ± 5 . 1 (mean ± standard deviation), and OCKF trials included knee flexion angles between 2.0 ± 4 . 8 and 97.5 ± 12 . 2 . Walking trials included internal/external rotation angles between 10.6 ± 7 and 6.9 ± 12 . 7 , and OCKF trials included internal/external rotation angles between 6.5 ± 8 . 3 and 9.1 ± 8 . 3 .
Data from motion capture indicated ranges of knee flexion/extension angles between 2.0 ± 4 . 8 and 97.5 ± 12 . 2 and internal/external rotation angles between 6.5 ± 8 . 3 and 9.1 ± 8 . 3 . All participants exhibited similar ranges of motion in flexion/extension and internal/external rotation and were within normal ranges of motion (see Figure 5). Participants completed an average of 14 movements in 60 s.

3.2. Modeling and Predictions

Figure 6 shows representative sensor data and its calculated derivative that was used as explanatory variables in the predictive models. Table A1 and Table A2 summarize the results from all 18 models. Adaptive boosting of a random forest regressor (RFR) predicted flexion/extension angles and internal/external rotation most accurately. The accuracy of the model’s predictions is indicated by both the R2 and RMSE values. Figure 7 shows the performance of each model in the two rotational degrees of freedom of interest. The average coefficient of determination (R2) of flexion/extension and internal/external rotation for the final subject-specific models were 0.940 and 0.731, respectively. The average RMSE across all subjects was 7.6° for flexion/extension and 1.8° for internal/external rotation. By including the derivatives of the sensor outputs in the model, the average RMSE improvement in the flexion/extension and internal/external rotation degrees of freedom were −16.9% and −9.3%, respectively. The average improvement of the R2 values in the flexion/extension and internal/external rotation degrees of freedom were 1.7% and 24.7%, respectively. All reported values are from models that included the signal’s derivative. Table 2 shows the accuracy of the device-specific models.
Figure 8 shows the dependency of the model’s bias and standard deviation as a function of angle of flexion/extension (Figure 8a,b) and internal/external rotation (Figure 8c,d). Figure 8 demonstrate that there is the most bias at the extreme of the model’s flexion/extension angle, yet the smallest standard deviation. The same holds for the model’s internal/external rotation. The largest standard deviation occurs in the middle of the range of knee angles (for both degrees of freedom), where the average bias is the smallest.
Using the presented pipeline, shown graphically in Figure 4, the resulting predictions and actual knee angles in Figure 7 were calculated to create the subject-specific models. The models correctly predicted the frequency of a motion and reasonably captured the magnitude of the extrema during OCKF. However, the magnitude of the models’ underestimation at the extrema is consistent (as shown in the bias plots in Figure 8a,c). Furthermore, the predictive model’s largest errors occurred at angles less than or equal to 0° flexion (with a straightened knee). This agrees with the bias results for all subjects shown in Figure 8.
The residuals of the predictive models’ bias showed a strong linear correlation with the actual knee angle, see Figure 8. Adding a linear regression model of the residuals to correct the bias was analyzed but proved to be statistically insignificant and was not included in these results. A paired t-test, between the models, including the linear regression of the residuals’ bias and those without, resulted in p-values greater than or equal to 0.191 for the R2 values and RMSE of both knee angles.

4. Discussion

An instrumented knee sleeve was developed, presented, and modeled using an adaptively boosted RFR model. The presented knee sleeve is a low-cost alternative to measuring knee kinematics outside a laboratory setting with comparable accuracy to current technologies. Unfortunately, during testing on two subjects, connectivity was lost between the sensors and the printed circuit board due to severed wires. Thus, the results from these two subjects resulted in failed models and were excluded from the predictive model analysis.
The predictive models’ most significant bias occurred at angles less than or equal to 0° flexion (an extended knee). This can be partially due to the laxity in the knee sleeve once the garment began to slip. After full flexion, the sleeve exhibited some buckling when the leg extended again. This buckling also occurred in several strain gauges near the patella and the insertion of the vastus lateralis, likely leading the model to believe the knee was in flexion (when in actuality, it was hyperextended), thus underestimating the peaks of the knee angle. The predictive models exhibited the most significant deviation in the middle of the range of motion. This is most likely caused by the high angular velocities at these angles and related to the inherent hysteresis caused by the polymer-based piezoresistive sensors’ viscoelastic nature. This has been seen previously in [54,60], and an attempt to compensate for this effect post hoc can be seen in [23]. Similar to all polymer-based wearable sensors, the sensors used in the current work exhibited viscoelastic effects such as creep, drift, and hysteresis. The magnitude of these effects was not found to affect the accuracy of the current analysis, and additional work characterizing these important sensor properties is currently underway. The impact of these errors is similar to previous IMUs [61,62,63,64,65] and high deflection strain gauges [54,66,67] used to measure joint angles. By including the derivatives of the sensor outputs in the model, the average improvement of the RMSE in the flexion/extension and internal/external rotation degrees of freedom were −16.9% and −9.3%, respectively. The average increase of the R2 values in the flexion/extension and internal/external rotation degrees of freedom were 1.7% and 24.7%, respectively. Accounting for differences in angular velocities significantly decreased the RMSE of the models developed in this work; accounting for angular velocities in optical fiber strain gauges results in similar improvements [68]. All results from the predictive models are for within-subject predictions, and the prediction of angles for previously unmeasured subjects would still have to be addressed.
Numerous studies have looked at methods to maximize the efficacy of at-home rehabilitation with biofeedback systems, namely teleconferencing with a PT [69,70,71,72], inertial measurement units (IMUs) [73,74], Nintendo Wii based systems/virtual reality (VR) [75,76,77,78], or PT guidance with robotic assistance [79]. Lockdowns associated with the COVID-19 pandemic have accelerated the development of devices for at-home rehabilitation [80]. However, telerehabilitation—rehabilitation that is still under the active guidance of a PT via video or telephone conferencing—continues to be the norm [15,69,71,72,81]. Telerehabilitation does little to relieve the increasing demand for TKA rehabilitation due to its reliance on one-to-one patient-to-therapist structure. Affordable IMUs are portable solutions that work for brief periods. However, issues with sensor drift [82] and calibration sensitivity [83] keep this solution from being widely adopted. Nintendo Wii-based systems are another promising solution. However, Nintendo Wii and VR-based systems address only stationary rehabilitation exercises, namely lateral weight shifting, multidirectional balance, and postural control [84]. Lastly, robotic assistance devices are the most accurate of these alternatives. However, they are also the most costly. Most robotic assistance measures are employed for inpatient use (as in [79]) but may be engineered for at-home use in the future. The sensor system presented in this work may be integrated into the current technologies to enhance their effectiveness or address their shortcomings. For example, the presented device may enhance telerehabilitation by relieving the workload of PTs. By combining the two approaches, most rehabilitation may be completed without direct PT observance. The presented device also allows the rehabilitation of many patients simultaneously, and a PT may intervene when necessary. This approach can increase the patient-to-therapist ratio without unduly increasing a PTs workload, yet still provide quantifiable feedback on a patient’s progress.
As noted in the Results section, our device could not fully capture the magnitude of the peaks and valleys in the flexion/extension or internal/external rotation degrees of freedom for all knee flexions for all subjects. Similar to other polymer-based devices [24], our device underestimated the extremes. The KneeHapp device developed by Haladjian et al. in [85] is a similar device that exhibited errors of similar magnitude to this work (4.82 ± 3.92°) for the flexion/extension angle. Other technologies, namely IMU-based devices (such as the one seen in [65]), overestimate the extremes. Additional filtering techniques are out of the scope of this work. However, they could be used in future works to mitigate this error since it was fairly consistent for each subject; i.e., the amount of underestimation was of a similar proportion to the actual peak angle for each knee flexion performed within each subject. The device presented in this work represents a cost-efficient instrument that can enhance at-home rehabilitation for TKA patients by providing a quantifiable measurement of two rotational degrees of freedom.
The present work focuses on the development and subject-specific modeling. However, the longer-term goals of the project are to generate device-specific models, and eventually a universal model that can be applied to all subjects. These goals will likely require a more inclusive incorporation of the viscoelastic characteristics of the sensors, as well as the inclusion of subject-specific demographic and anthropometric information.
We do note three limitations to this study. First, motion capture data were downsampled to equal the sampling rate of the instrumented knee sleeve resulting in some interpolation. Although data were interpolated with a quadratic fit with continuous derivatives, this step may have caused minute errors in the flexion/extension or internal/external rotation angles referenced as ground truth throughout the data collection. Second, the presented device was not adhered to the skin’s surface. Therefore, motion of the device relative to the skin’s surface within the data are possible and may add errors to the data if the device was not positioned correctly or was loosely fitted. This may be especially relevant to closed chain movement where the potential for movement of the device is amplified. Similar limitations have also been observed previously by Gholami et al. [33]. Third, there is debate around the accuracy of motion capture methods when measuring tibia internal/external rotation. Merriaux et al. found the three-dimensional positioning error of a motion capture system to be 2 mm during dynamic motion [86]. When using motion capture to measure internal/external rotation, Keizer and Otten concluded that rotations less than 1.70 degrees should be taken with caution due to the sum of errors with the camera system [87]. Additional errors can be introduced from improper marker placement on a limb’s anatomical landmarks, as well as from skin artifacts [88,89]. To minimize errors related to improper placement or motion relative to the skin, the same researcher placed markers on all participants to eliminate inter-researcher placement error. The researcher also used double-sided adhesive to attach the markers to all participants. Other techniques to more confidently measure bone motion may be used to mitigate this potential source of error. However, the sensor prediction data shown here reflected the motion capture data, despite the debate of the accuracy of the technique used.
In conclusion, the instrumented knee sleeve presented in this work accurately predicted two rotational degrees of freedom within the knee during OCKF. The viscoelastic nature of the piezoresistive sensors used to predict the knee angle was apparent and a source of error in predicting two rotational degrees of freedom at high angular velocities. Nevertheless, the instrumented knee sleeve presented in this work proved to be a repeatable, reusable, low-cost device that can adequately model the flexion/extension and internal/external rotation angles in the knee.

Author Contributions

Conceptualization, D.S.W., D.T.F., M.K.S. and A.E.B.; methodology, D.S.W. and A.E.B.; software, D.S.W.; validation, D.S.W.; formal analysis, D.S.W., W.F.C., J.G. and A.S.; investigation, D.S.W., K.J., A.C., H.L. and H.D.; resources, M.K.S. and U.H.M.; data curation, D.S.W., K.J. and A.C.; writing—original draft preparation, D.S.W. and A.C.; writing—review and editing, K.J., J.G., A.S., D.T.F., M.K.S., U.H.M., W.F.C. and A.E.B.; visualization, D.S.W.; supervision, A.E.B.; project administration, A.E.B.; funding acquisition, D.T.F., M.K.S., U.H.M. and A.E.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation grant number CBET-1901845.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Brigham Young University (X2020-239, date of approval 17 May 2021).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. RMSE and R2 values for all models tested on flexion\slash extension. The failed column shows the number of subjects predicted with negative R2 values from each model. All models were optimized using Bayesian or grid search optimization. Data presented in this work were calculated with the adaptive boosted regression (RFR) model.
Table A1. RMSE and R2 values for all models tested on flexion\slash extension. The failed column shows the number of subjects predicted with negative R2 values from each model. All models were optimized using Bayesian or grid search optimization. Data presented in this work were calculated with the adaptive boosted regression (RFR) model.
Model NameRMSER2Failed
MeanSDMaxMinMeanSDMaxMin
Linear Regression14.183.6823.789.450.790.160.930.346
Gradient Boosting Regression14.7810.2733.964.540.720.370.990.091
Adaptive Boosted Regression (DTR)7.752.5114.344.110.940.040.990.851
RFR8.342.6213.544.510.930.050.980.831
Multilayer Perceptron Regression14.505.3226.588.780.670.330.930.014
Ridge14.083.5623.279.440.800.150.930.386
Bayesian Ridge14.093.5523.249.450.800.150.930.396
Elastic Net16.266.0729.409.300.690.280.930.023
PLS Regression17.305.1528.0710.150.690.220.920.134
K-Neighbors Regression10.042.4414.925.380.900.050.980.811
SGD Regression15.354.7124.839.700.740.210.920.184
Lasso14.563.3222.309.420.800.080.930.634
SVR16.337.3229.318.860.700.290.930.086
NuSVR27.237.2335.8515.120.310.290.770.0311
Bagging Regression (RFR)24.495.5330.0212.790.490.190.890.151
Linear Regression (PCA)15.254.9428.4710.060.760.160.880.231
RFR (PCA)8.242.6313.714.530.930.040.980.851
Adaptive Boosted Regression (RFR)7.572.5513.584.190.940.040.990.861
Table A2. RMSE and R2 values for all models tested on internal\slash external rotation. The failed column shows the number of subjects predicted with negative R2 values from each model. All models were optimized using Bayesian or grid search optimization. Data presented in this work were calculated with the adaptive boosted regression (RFR) model.
Table A2. RMSE and R2 values for all models tested on internal\slash external rotation. The failed column shows the number of subjects predicted with negative R2 values from each model. All models were optimized using Bayesian or grid search optimization. Data presented in this work were calculated with the adaptive boosted regression (RFR) model.
Model NameRMSER2Failed
MeanSDMaxMinMeanSDMaxMin
Linear Regression2.771.145.241.630.650.130.790.389
Gradient Boosting Regression3.272.219.311.250.380.380.960.014
Adaptive Boosted Regression (DTR)1.850.482.971.090.710.210.970.222
RFR1.890.462.921.230.700.230.970.222
Multilayer Perceptron Regression2.320.552.911.280.740.120.930.486
Ridge2.751.125.171.620.650.130.790.399
Bayesian Ridge2.741.125.171.600.660.130.790.419
Elastic Net2.651.054.701.360.580.240.780.077
PLS Regression2.791.125.101.350.550.220.770.127
K-Neighbors Regression2.160.754.161.270.650.240.920.192
SGD Regression3.482.178.861.560.540.280.800.047
Lasso3.102.018.041.420.580.230.780.057
SVR2.391.325.271.220.620.270.880.199
NuSVR3.041.265.061.830.460.270.740.0411
Bagging Regression (RFR)4.052.449.921.330.360.230.710.035
Linear Regression (PCA)2.971.355.741.380.500.240.710.014
RFR (PCA)1.860.493.021.210.710.210.970.262
Adaptive Boosted Regression (RFR)1.800.462.861.110.730.200.970.291

Abbreviations

The following abbreviations are used in this manuscript:
OCKFopen chain knee flexion
RFRrandom forest regressor
DTRdecision tree regressor
SGDstochastic gradient descent
PLSpartial least squares
SVRsupport vector regression
PCAprincipal component analysis
VRvirtual reality
PTphysical therapist

References

  1. Kurtz, S.M.; Ong, K.L.; Lau, E.; Bozic, K.J. Impact of the economic downturn on total joint replacement demand in the United States: Updated projections to 2021. J. Bone Jt. Surg. 2014, 96, 624–630. [Google Scholar]
  2. Inacio, M.C.; Paxton, E.W.; Graves, S.E.; Namba, R.S.; Nemes, S. Projected increase in total knee arthroplasty in the United States—An alternative projection model. Osteoarthr. Cartil. 2017, 25, 1797–1803. [Google Scholar] [CrossRef] [Green Version]
  3. Kurtz, S.; Ong, K.; Lau, E.; Mowat, F.; Halpern, M. Projections of Primary and Revision Hip and Knee Arthroplasty in the United States from 2005 to 2030. J. Bone Jt. Surg. 2007, 89, 780–785. [Google Scholar] [CrossRef]
  4. Sloan, M.; Premkumar, A.; Sheth, N.P. Future Demand for Total Joint Arthroplasty Drives Renewed Interest in Arthroplasty Fellowship. HSS J. 2020, 16, 210–215. [Google Scholar] [CrossRef]
  5. Herzog, M.M.; Marshall, S.W.; Lund, J.L.; Pate, V.; Mack, C.D.; Spang, J.T. Trends in Incidence of ACL Reconstruction and Concomitant Procedures Among Commercially Insured Individuals in the United States, 2002–2014. Sport. Health 2018, 10, 523–531. [Google Scholar] [CrossRef]
  6. Landry, M.D.; Ricketts, T.C.; Verrier, M.C. The precarious supply of physical therapists across Canada: Exploring national trends in health human resources (1991 to 2005). Hum. Resour. Health 2007, 5, 23. [Google Scholar] [CrossRef] [Green Version]
  7. Jesus, T.S.; Landry, M.D.; Hoenig, H.; Dussault, G.; Koh, G.C.; Fronteira, I. Is Physical Rehabilitation Need Associated With the Rehabilitation Workforce Supply? An Ecological Study Across 35 High-Income Countries. Int. J. Health Policy Manag. 2020, 2017, 1–9. [Google Scholar] [CrossRef]
  8. Lo, K.; Curtis, H.; Keating, J.L.; Bearman, M. Physiotherapy clinical educators’ perceptions of student fitness to practise. BMC Med. Educ. 2017, 17, 16. [Google Scholar] [CrossRef] [Green Version]
  9. Lavernia, C.J.; D’Apuzzo, M.R.; Hernandez, V.H.; Lee, D.J.; Rossi, M.D. Postdischarge Costs in Arthroplasty Surgery. J. Arthroplast. 2006, 21, 144–150. [Google Scholar] [CrossRef]
  10. Tran, V.; Doctor, D.; Nelson, C. The cost of postoperative home health extended care rehabilitation facilities for total joint replacement patients. In Proceedings of the American Academy of Orthopaedic Surgeons Annual Meeting, Anaheim, CA, USA, 4–8 February 1999. [Google Scholar]
  11. Wang, W.L.; Rondon, A.J.; Tan, T.L.; Wilsman, J.; Purtill, J.J. Self-Directed Home Exercises vs Outpatient Physical Therapy After Total Knee Arthroplasty: Value and Outcomes Following a Protocol Change. J. Arthroplast. 2019, 34, 2388–2391. [Google Scholar] [CrossRef]
  12. Joice, M.G.; Bhowmick, S.; Amanatullah, D.F. Perioperative physiotherapy in total knee arthroplasty. Orthopedics 2017, 40, e765–e773. [Google Scholar] [CrossRef] [Green Version]
  13. Kairy, D.; Tousignant, M.; Leclerc, N.; Côté, A.M.; Levasseur, M. The patient’s perspective of in-home telerehabilitation physiotherapy services following total knee arthroplasty. Int. J. Environ. Res. Public Health 2013, 10, 3998–4011. [Google Scholar] [CrossRef]
  14. Oatis, C.A.; Li, W.; Dirusso, J.M.; Hoover, M.J.; Johnston, K.K.; Butz, M.K.; Phillips, A.L.; Nanovic, K.M.; Cummings, E.C.; Rosal, M.C.; et al. Variations in Delivery and Exercise Content of Physical Therapy Rehabilitation Following Total Knee Replacement Surgery: A Cross-Sectional Observation Study HHS Public Access. Int. J. Phys. Med. Rehabil. 2014, 5. [Google Scholar] [CrossRef] [Green Version]
  15. Naeemabadi, M.; Fazlali, H.; Najafi, S.; Dinesen, B.; Hansen, J. Telerehabilitation for Patients With Knee Osteoarthritis: A Focused Review of Technologies and Teleservices. JMIR Biomed. Eng. 2020, 5, e16991. [Google Scholar] [CrossRef]
  16. Yee, M.J.; Mubarak, N.M.; Abdullah, E.C.; Khalid, M.; Walvekar, R.; Karri, R.R.; Nizamuddin, S.; Numan, A. Carbon Nanomaterials Based Films for Strain Sensing Application—A Review. Nano-Struct. Nano-Objects 2019, 18, 100312. [Google Scholar]
  17. Lu, Y.; Biswas, M.C.; Guo, Z.; Jeon, J.W.; Wujcik, E.K. Recent developments in bio-monitoring via advanced polymer nanocomposite-based wearable strain sensors. Biosens. Bioelectron. 2019, 123, 167–177. [Google Scholar] [CrossRef]
  18. Abd, M.A.; Al-Saidi, M.; Lin, M.; Liddle, G.; Mondal, K.; Engeberg, E.D. Surface Feature Recognition and Grasped Object Slip Prevention with a Liquid Metal Tactile Sensor for a Prosthetic Hand. In Proceedings of the IEEE RAS and EMBS International Conference on Biomedical Robotics and Biomechatronics, New York, NY, USA, 29 November–1 December 2020; pp. 1174–1179. [Google Scholar] [CrossRef]
  19. Abd, M.A.; Paul, R.; Aravelli, A.; Bai, O.; Lagos, L.; Lin, M.; Engeberg, E.D. Hierarchical Tactile Sensation Integration from Prosthetic Fingertips Enables Multi-Texture Surface Recognition. Sensors 2021, 21, 4324. [Google Scholar] [CrossRef]
  20. Park, Y.L.; Chen, B.R.; Wood, R.J. Design and fabrication of soft artificial skin using embedded microchannels and liquid conductors. IEEE Sens. J. 2012, 12, 2711–2718. [Google Scholar] [CrossRef]
  21. Ying, B.; Chen, R.Z.; Zuo, R.; Li, J.; Liu, X. An Anti-Freezing, Ambient-Stable and Highly Stretchable Ionic Skin with Strong Surface Adhesion for Wearable Sensing and Soft Robotics. Adv. Funct. Mater. 2021, 31, 2104665. [Google Scholar]
  22. Kazanskiy, N.L.; Butt, M.A.; Khonina, S.N. Recent Advances in Wearable Optical Sensor Automation Powered by Battery versus Skin-like Battery-Free Devices for Personal Healthcare—A Review. Nanomaterials 2022, 12, 334. [Google Scholar]
  23. Leal-Junior, A.G.; Frizera-Neto, A.; Pontes, M.J.; Botelho, T.R. Hysteresis compensation technique applied to polymer optical fiber curvature sensor for lower limb exoskeletons. Meas. Sci. Technol. 2017, 28, 125103. [Google Scholar]
  24. Rezende, A.; Alves, C.; Marques, I.; Silva, M.A.; Naves, E. Polymer optical fiber goniometer: A new portable, low cost and reliable sensor for joint analysis. Sensors 2018, 18, 4293. [Google Scholar] [CrossRef] [Green Version]
  25. Donno, M.; Palange, E.; Di Nicola, F.; Member, S.; Bucci, G.; Ciancetta, F. A New Flexible Optical Fiber Goniometer for Dynamic Angular Measurements: Application to Human Joint Movement Monitoring. IEEE Trans. Instrum. Meas. 2008, 57, 1614–1620. [Google Scholar] [CrossRef]
  26. Liu, Q.; Chen, J.; Li, Y.; Shi, G. High-Performance Strain Sensors with Fish-Scale-Like Graphene-Sensing Layers for Full-Range Detection of Human Motions. ACS Nano 2016, 10, 7901–7906. [Google Scholar] [CrossRef]
  27. Gholami, M.; Rezaei, A.; Cuthbert, T.J.; Napier, C.; Menon, C. Lower Body Kinematics Monitoring in Running Using Fabric-Based Wearable Sensors and Deep Convolutional Neural Networks. Sensors 2019, 19, 5325. [Google Scholar] [CrossRef] [Green Version]
  28. Li, Y.Q.; Zhu, W.B.; Yu, X.G.; Huang, P.; Fu, S.Y.; Hu, N.; Liao, K. Multifunctional Wearable Device Based on Flexible and Conductive Carbon Sponge/Polydimethylsiloxane Composite. ACS Appl. Mater. Interfaces 2016, 8, 33189–33196. [Google Scholar] [CrossRef]
  29. Amjadi, M.; Yoon, Y.J.; Park, I. Ultra-stretchable and skin-mountable strain sensors using carbon nanotubes–Ecoflex nanocomposites. Nanotechnology 2015, 26, 744–747. [Google Scholar] [CrossRef]
  30. Luhrs, C.C.; Daskam, C.D.; Gonzalez, E.; Phillip, J. Fabrication of a low density carbon fiber foam and its characterization as a strain gauge. Materials 2014, 7, 3699–3714. [Google Scholar] [CrossRef] [Green Version]
  31. Singha, K.; Kumar, J.; Pandit, P. Recent advancements in wearable & smart textiles: An overview. Mater. Today Proc. 2019, 16, 1518–1523. [Google Scholar] [CrossRef]
  32. Shyr, T.W.; Shie, J.W.; Jiang, C.H.; Li, J.J. A textile-based wearable sensing device designed for monitoring the flexion angle of elbow and knee movements. Sensors 2014, 14, 4050–4059. [Google Scholar] [CrossRef] [Green Version]
  33. Gholami, M.; Ejupi, A.; Rezaei, A.; Ferrone, A.; Menon, C. Estimation of Knee Joint Angle Using a Fabric-Based Strain Sensor and Machine Learning: A Preliminary Investigation. In Proceedings of the 2018 7th IEEE International Conference on Biomedical Robotics and Biomechatronics (Biorob), Enschede, The Netherlands, 26–29 August 2018. [Google Scholar] [CrossRef]
  34. Levens, A.; Inman, V.T.; Blosser, J. Transverse rotation of the segments of the lower extremity in locomotion. J. Bone Jt. Surg. 1948, 30, 859–872. [Google Scholar]
  35. Lützner, J.; Kirschner, S.; Günther, K.P.; Harman, M.K. Patients with no functional improvement after total knee arthroplasty show different kinematics. Int. Orthop. 2012, 36, 1841. [Google Scholar] [CrossRef] [Green Version]
  36. Matsuzaki, T.; Matsumoto, T.; Muratsu, H.; Kubo, S.; Matsushita, T.; Kawakami, Y.; Ishida, K.; Oka, S.; Kuroda, R.; Kurosaka, M. Kinematic factors affecting postoperative knee flexion after cruciate-retaining total knee arthroplasty. Int. Orthop. 2013, 37, 803. [Google Scholar] [CrossRef] [Green Version]
  37. Nicoll, D.; Rowley, D.I. Internal rotational error of the tibial component is a major cause of pain after total knee replacement. J. Bone Jt. Surg. Ser. B 2010, 92, 1238–1244. [Google Scholar] [CrossRef]
  38. Meccia, B.; Komistek, R.D.; Mahfouz, M.; Dennis, D. Abnormal axial rotations in TKA contribute to reduced weightbearing flexion knee. Clin. Orthop. Relat. Res. 2014, 472, 248–253. [Google Scholar] [CrossRef]
  39. Gardner, C.; Johnson, O.; Fullwood, D.; Hansen, G.; Adams, B. Piezoresistive Effect in Nickel Nanostrand-Polymer Composites. In Proceedings of the TMS 2009 Annual Meeting & Exhibition, San Francisco, CA, USA, 15–19 February 2009. [Google Scholar]
  40. Johnson, O.; Gardner, C.; Fullwood, D.; Adams, B.; Hansen, G.; Kalidindi, S. Textures of Dispersion of Nickel Nanostrand Composites, and Modeling of Piezoresistive Behavior. In Proceedings of the Materials Science & Technology, Pittsburgh, PA, USA, 25–29 October 2009. [Google Scholar]
  41. Calkins, T.B.; Fullwood, D.T.; Ghosh, S.; Hyatt, T.B.; Johnson, O.K.; Hansen, N.; Hansen, G. Applications For A Nano-composite High Displacement Strain Gauge. In Proceedings of the 42nd ISTC, Salt Lake City, UT, USA, 11–14 October 2010. [Google Scholar]
  42. Johnson, O.K.; Gardner, C.J.; Fullwood, D.T.; Adams, B.L.; Hansen, N.; Hansen, G. The colossal piezoresistive effect in nickel nanostrand polymer composites and a quantum tunneling model. Comput. Mater. Contin. 2010, 15, 87–111. [Google Scholar] [CrossRef]
  43. Johnson, O.; Mara, N.; Kaschner, G.; Mason, T.; Fullwood, D.; Adams, B. Multi-scale Model for the Extreme Piezoresistivity in Silicone/Nickel “Nanostrand”/Nickel Coated Carbon Fiber Nanocomposite. Metall. Mater. Trans. A 2010, 42, 3898–3906. [Google Scholar] [CrossRef]
  44. Johnson, O.K.; Kaschner, G.C.; Mason, T.A.; Fullwood, D.T.; Hansen, G. Optimization of nickel nanocomposite for large strain sensing applications. Sens. Actuators A Phys. 2011, 166, 40–47. [Google Scholar] [CrossRef]
  45. Johnson, T.M.; Fullwood, D.T.; Hansen, G. Strain monitoring of carbon fiber composite via embedded nickel nano-particles. Compos. Part B Eng. 2012, 43, 1155–1163. [Google Scholar]
  46. Koecher, M.C.; Pande, J.H.; Merkley, S.S.; Henderson, S.; Fullwood, D.T. Evaluation of nickel nanostrands as strain sensors in CFRP. In Proceedings of the SAMPE 2012 Conference, Baltimore, MD, USA, 21–24 May 2012. [Google Scholar]
  47. Remington, T.D.; Merrell, A.; Stolworthy, D.; McArthur, D.; Fullwood, D.; Bowden, A.; Hansen, N. Biomechanical applications of nano-composite strain gauges. In Proceedings of the SAMPE 2013 Conference, Long Beach, CA, USA, 6–9 June 2012. [Google Scholar]
  48. Koecher, M.; Yeager, J.D.; Park, T.; Fullwood, D.; Colton, J.S.; Mara, N.; Hansen, N. Characterization of nickel nanostrand nanocomposites through dielectric spectroscopy and nanoindentation. Polym. Eng. Sci. 2013, 53, 2666–2673. [Google Scholar]
  49. Remington, D. Biomechanical Applications and Modeling of Quantum Nano-Composite Strain Gauges. Ph.D. Thesis, Brigham Young University, Provo, UT, USA, 2014. [Google Scholar]
  50. Adam Bilodeau, R.; Fullwood, D.T.; Colton, J.S.; Yeager, J.D.; Bowden, A.E.; Park, T. Evolution of nano-junctions in piezoresistive nanostrand composites. Compos. Part B Eng. 2015, 72, 45–52. [Google Scholar] [CrossRef] [Green Version]
  51. Hansen, N.; Adams, D.O.; Fullwood, D.T. Evaluation and development of electrical conductivity models for nickel nanostrand polymer composites. Polym. Eng. Sci. 2015, 55, 549–557. [Google Scholar]
  52. Koecher, M.C.; Pande, J.H.; Merkley, S.; Henderson, S.; Fullwood, D.T.; Bowden, A.E. Piezoresistive in-situ strain sensing of composite laminate structures. Compos. Part B Eng. 2015, 69, 534–541. [Google Scholar]
  53. Baradoy, D.A. Composition Based Modaling of Silicone Nano-Composite Strain Gauges. Ph.D. Thesis, Brigham Young University, Provo, UT, USA, 2015. [Google Scholar]
  54. Martineau, A.D. Estimation of Knee Kinematics Using Non-Monotonic Nanocomposite High-Deflection Strain Gauges. Ph.D. Thesis, Brigham Young University, Provo, UT, USA, 2018. [Google Scholar]
  55. Clayton, M. Modeling Piezoresistive Effects in Flexible Sensors. Ph.D. Thesis, Brigham Young University, Provo, UT, USA, 2019. [Google Scholar]
  56. Mitchell, U.H.; Lee, H.; Dennis, H.E.; Seeley, M.K. Quality of knee strengthening exercises performed at home deteriorates after one week. BMC Musculoskelet. Disord. 2022, 23, 164. [Google Scholar]
  57. Wood, D.S.; Jensen, K.; Fullwood, D.; Seeley, M.; Mitchell, U.; Christensen, W.; Bowden, A. General Method of Measuring Skin Strain of the Anterior Knee using Linear Strain Theory and Application. CMBBE, 2021; under review. [Google Scholar]
  58. Winter, D.A. Biomechanics and Motor Control of Human Movement; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  59. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  60. Amjadi, M.; Kyung, K.U.; Park, I.; Sitti, M. Stretchable, Skin-Mountable, and Wearable Strain Sensors and Their Potential Applications: A Review. Adv. Funct. Mater. 2016, 26, 1678–1698. [Google Scholar]
  61. Bakhshi, S.; Mahoor, M.H.; Davidson, B.S. Development of a body joint angle measurement system using IMU sensors. In Proceedings of the 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Boston, MA, USA, 30 August–3 September 2011; pp. 6923–6926. [Google Scholar] [CrossRef]
  62. Cooper, G.; Sheret, I.; McMillian, L.; Siliverdis, K.; Sha, N.; Hodgins, D.; Kenney, L.; Howard, D. Inertial sensor-based knee flexion/extension angle estimation. J. Biomech. 2009, 42, 2678–2685. [Google Scholar] [CrossRef] [Green Version]
  63. Fennema, M.C.; Bloomfield, R.A.; Lanting, B.A.; Birmingham, T.B.; Teeter, M.G. Repeatability of measuring knee flexion angles with wearable inertial sensors. Knee 2019, 26, 97–105. [Google Scholar] [CrossRef] [Green Version]
  64. Watanabe, T.; Saito, H. Tests of wireless wearable sensor system in joint angle measurement of lower limbs. In Proceedings of the 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Boston, MA, USA, 30 August–3 September 2011; pp. 5469–5472. [Google Scholar] [CrossRef]
  65. Ajdaroski, M.; Tadakala, R.; Nichols, L.; Esquivel, A. Validation of a Device to Measure Knee Joint Angles for a Dynamic Movement. Sensors 2020, 20, 1747. [Google Scholar] [CrossRef] [Green Version]
  66. Mohamed, A.A.; Baba, J.; Beyea, J.; Landry, J.; Sexton, A.; McGibbon, C.A. Comparison of strain-gage and fiber-optic goniometry for measuring knee kinematics during activities of daily living and exercise. J. Biomech. Eng. 2012, 134, 084502. [Google Scholar] [CrossRef] [PubMed]
  67. Nakamoto, H.; Yamaji, T.; Hirata, I.; Ootaka, H.; Kobayashi, F. Joint angle measurement by stretchable strain sensor. J. Ambient. Intell. Humaniz. Comput. 2018, 9, 1–6. [Google Scholar] [CrossRef]
  68. Avellar, L.; Leal-Junior, A.; Marques, C.; Frizera, A. Performance Analysis of a Lower Limb Multi Joint Angle Sensor Using CYTOP Fiber: Influence of Light Source Wavelength and Angular Velocity Compensation. Sensors 2020, 20, 326. [Google Scholar] [CrossRef] [Green Version]
  69. Tousignant, M.; Moffet, H.; Biossy, P.; Corriveau, H.; Cabana, F.; Marquis, F. A randomized controlled trial of home telerehabilitation for post-knee arthroplasty. J. Telemed. Telecare 2011, 17, 195–198. [Google Scholar] [PubMed]
  70. Moffet, H.; Tousignant, M.; Nadeau, S.; Mérette, C.; Boissy, P.; Corriveau, H.; Marquis, F.F.F.; Cabana, F.F.; Ranger, P.; Belzile, É.L.; et al. In-home telerehabilitation compared with faceto-face rehabilitation after total knee arthroplasty: A noninferiority randomized controlled trial. J. Bone Jt. Surg. Am. Vol. 2015, 97, 1129–1141. [Google Scholar] [CrossRef] [Green Version]
  71. Xie, S.H.; Wang, Q.; Wang, L.Q.; Wang, L.; Song, K.P.; He, C.Q. Effect of internet-based rehabilitation programs on improvement of pain and physical function in patients with knee osteoarthritis: Systematic review and meta-analysis of randomized controlled trials. J. Med. Internet Res. 2021, 23, e21542. [Google Scholar] [CrossRef]
  72. Jansson, M.M.; Rantala, A.; Miettunen, J.; Pikkarainen, M.; Puhto, A.P.; Pikkarainen, M. The effects and safety of telerehabilitation in patients with lower-limb joint replacement: A systematic review and narrative synthesis. J. Telemed. Telecare 2022, 28, 96–114. [Google Scholar] [CrossRef]
  73. Correia, F.D.; Nogueira, A.; Magalhães, I.; Guimarães, J.; Moreira, M.; Barradas, I.; Teixeira, L.; Tulha, J.; Seabra, R.; Lains, J.; et al. Home-based Rehabilitation with A Novel Digital Biofeedback System versus Conventional In-person Rehabilitation after Total Knee Replacement: A feasibility study. Sci. Rep. 2018, 8, 11299. [Google Scholar] [CrossRef] [Green Version]
  74. Piqueras, M.; Marco, E.; Coll, M.; Escalada, F.; Ballester, A.; Cinca, C.; Belmonte, R.; Muniesa, J.M. Effectiveness of an interactive virtual telerehabilitation system in patients after total knee arthroplasty: A randomized controlled trial. J. Rehabil. Med. 2013, 45, 392–396. [Google Scholar] [CrossRef] [Green Version]
  75. Ficklscherer, A.; Stapf, J.; Meissner, K.M.; Niethammer, T.; Lahner, M.; Wagenhäuser, M.; Müller, P.E.; Pietschmann, M.F. Testing the feasibility and safety of the Nintendo Wii gaming console in orthopedic rehabilitation: A pilot randomized controlled study. Arch. Med. Sci. 2016, 12, 1273–1278. [Google Scholar] [CrossRef] [Green Version]
  76. Fung, V.; Ho, A.; Shaffer, J.; Chung, E.; Gomez, M. Use of Nintendo Wii Fit TM in the rehabilitation of outpatients following total knee replacement: A preliminary randomised controlled trial. Physiotherapy 2012, 98, 183–188. [Google Scholar] [CrossRef] [PubMed]
  77. Christiansen, C.L.; Bade, M.J.; Davidson, B.S.; Dayton, M.R.; Stevens-Lapsley, J.E. Effects of weight-bearing biofeedback training on functional movement patterns following total knee arthroplasty: A randomized controlled trial. J. Orthop. Sport. Phys. Ther. 2015, 45, 647–655. [Google Scholar] [CrossRef] [Green Version]
  78. Jin, C.; Feng, Y.; Ni, Y.; Shan, Z. Virtual reality intervention in postoperative rehabilitation after total knee arthroplasty: A prospective and randomized controlled clinical trial. Int. J. Clin. Exp. Med. 2018, 11, 6119–6124. [Google Scholar]
  79. Li, J.; Wu, T.; Xu, Z.; Gu, X. A pilot study of post-total knee replacement gait rehabilitation using lower limbs robot-assisted training system. Eur. J. Orthop. Surg. Traumatol. 2014, 24, 203–208. [Google Scholar] [CrossRef] [PubMed]
  80. Lim, M.A.; Pranata, R. Teleorthopedic: A Promising Option During and After the Coronavirus Disease 2019 (COVID-19) Pandemic. Front. Surg. 2020, 7, 62. [Google Scholar] [CrossRef]
  81. Eichler, S.; Rabe, S.; Salzwedel, A.; Müller, S.; Stoll, J.; Tilgner, N.; John, M.; Wegscheider, K.; Mayer, F.; Völler, H. Effectiveness of an interactive telerehabilitation system with home-based exercise training in patients after total hip or knee replacement: Study protocol for a multicenter, superiority, no-blinded randomized controlled trial. Trials 2017, 18, 438. [Google Scholar] [CrossRef] [Green Version]
  82. Zhao, Y.; Horemuz, M.; Sjöberg, L.E. Stochastic modelling and analysis of IMU sensor errors. Arch. Fotogram. Kartogr. Teledetekcji 2011, 22, 437–449. [Google Scholar]
  83. Walchko, K.J.; Mason, P.A. Inertial navigation. In Proceedings of the 2002 Florida Conference on Recent Advances in Robotics, Gainesville, FL, USA, 10–11 May 2002; pp. 1–9. [Google Scholar]
  84. Byra, J.; Czernicki, K. The Effectiveness of Virtual Reality Rehabilitation in Patients with Knee and Hip Osteoarthritis. J. Clin. Med. 2020, 9, 2639. [Google Scholar] [CrossRef]
  85. Haladjian, J.; Bredies, K.; Brugge, B. KneeHapp textile: A smart textile system for rehabilitation of knee injuries. In Proceedings of the 2018 IEEE 15th International Conference on Wearable and Implantable Body Sensor Networks, Las Vegas, NV, USA, 4–7 March 2018. [Google Scholar] [CrossRef]
  86. Merriaux, P.; Dupuis, Y.; Boutteau, R.; Vasseur, P.; Savatier, X. A Study of Vicon System Positioning Performance. Sensors 2017, 17, 1591. [Google Scholar] [CrossRef]
  87. Keizer, M.N.; Otten, E. Technical note: Sensitivity analysis of the SCoRE and SARA methods for determining rotational axes during tibiofemoral movements using optical motion capture. J. Exp. Orthop. 2020, 7, 6. [Google Scholar] [CrossRef] [Green Version]
  88. Taylor, W.R.; Ehrig, R.M.; Duda, G.N.; Schell, H.; Seebeck, P.; Heller, M.O. On the influence of soft tissue coverage in the determination of bone kinematics using skin markers. J. Orthop. Res. 2005, 23, 726–734. [Google Scholar] [CrossRef] [PubMed]
  89. Akbarshahi, M.; Schache, A.G.; Fernandez, J.W.; Baker, R.; Banks, S.; Pandy, M.G. Non-invasive assessment of soft-tissue artifact and its effect on knee joint kinematics during functional activity. J. Biomech. 2010, 43, 1292–1301. [Google Scholar] [CrossRef] [PubMed]
Figure 1. SEM images of the two electrically conductive filler materials. (a) SEM of NiNs. (b) SEM of NCCFs.
Figure 1. SEM images of the two electrically conductive filler materials. (a) SEM of NiNs. (b) SEM of NCCFs.
Sensors 22 02499 g001
Figure 2. (a) shows the placement of the instrumented knee sleeve and markers used in the data collection portion of the study. The placement of 16 sensors within the knee sleeve are shown in yellow in (b). (a) Marker/sleeve placement. (b) Sensor placement.
Figure 2. (a) shows the placement of the instrumented knee sleeve and markers used in the data collection portion of the study. The placement of 16 sensors within the knee sleeve are shown in yellow in (b). (a) Marker/sleeve placement. (b) Sensor placement.
Sensors 22 02499 g002
Figure 3. (a) Demonstration of open-chain knee flexion. (b) The captured motion with marker placement.
Figure 3. (a) Demonstration of open-chain knee flexion. (b) The captured motion with marker placement.
Sensors 22 02499 g003
Figure 4. Flow chart of the pipeline used to predict knee angles.
Figure 4. Flow chart of the pipeline used to predict knee angles.
Sensors 22 02499 g004
Figure 5. Knee angle ranges of all subjects of two rotational degrees of freedom. The red line represents the mean of all subjects. The blue lines represent the extrema. Positive angles correspond to knee flexion and internal rotation of the tibia. (a) Flexion/extension angle. (b) Tibia internal/external rotation.
Figure 5. Knee angle ranges of all subjects of two rotational degrees of freedom. The red line represents the mean of all subjects. The blue lines represent the extrema. Positive angles correspond to knee flexion and internal rotation of the tibia. (a) Flexion/extension angle. (b) Tibia internal/external rotation.
Sensors 22 02499 g005
Figure 6. Representative sensor data and its derivative from 30 s of open chain knee flexion.
Figure 6. Representative sensor data and its derivative from 30 s of open chain knee flexion.
Sensors 22 02499 g006
Figure 7. (a,b) show the distribution of the best model’s RMSE and R2 value in both degrees of freedom for subject-specific models. Red bars denote the overall average and whiskers denote standard deviations of the dataset. Samples of the models’ outputs (red) and actual knee angles (blue) are also shown. Positive angles correspond to knee flexion and internal rotation. (c,d) show representative data from one trial. (a) RMSE. (b) R2. (c) Predicted (red) versus actual (blue) flexion/extension knee angle. (d) Predicted (red) versus actual (blue)internal/external rotation angle.
Figure 7. (a,b) show the distribution of the best model’s RMSE and R2 value in both degrees of freedom for subject-specific models. Red bars denote the overall average and whiskers denote standard deviations of the dataset. Samples of the models’ outputs (red) and actual knee angles (blue) are also shown. Positive angles correspond to knee flexion and internal rotation. (c,d) show representative data from one trial. (a) RMSE. (b) R2. (c) Predicted (red) versus actual (blue) flexion/extension knee angle. (d) Predicted (red) versus actual (blue)internal/external rotation angle.
Sensors 22 02499 g007
Figure 8. Distributions of all subjects’ bias and standard deviations in both degrees of freedom. Blue dots represent the local average across all subjects within 2 degree ranges. Lines represent the average of all subjects.
Figure 8. Distributions of all subjects’ bias and standard deviations in both degrees of freedom. Blue dots represent the local average across all subjects within 2 degree ranges. Lines represent the average of all subjects.
Sensors 22 02499 g008
Table 1. Biometric information of the cohort of subjects.
Table 1. Biometric information of the cohort of subjects.
GenderAgeHeight (cm)Weight (kg)BMI
M301859327.2
M29189101.728.5
M33179.587.627.2
M2918174.322.7
M2118372.821.7
M4717971.322.3
F22161.558.422.4
M24182.565.919.8
F4017474.824.7
F2518077.824.0
F2416554.219.9
F211686021.3
F2116859.120.9
F26162.7548.518.3
M20187.571.320.3
F31163.568.825.7
M2817793.829.9
F2316693.533.9
Table 2. RMSE and R 2 values of the device-specific models.
Table 2. RMSE and R 2 values of the device-specific models.
DeviceFlexion/ExtensionInt./Ext. RotationNumber of Participants
RMSE R 2 RMSE R 2
A18.5250.7235.5280.27114
B5.7440.9752.7140.8623
C13.5770.8592.2850.6761
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wood, D.S.; Jensen, K.; Crane, A.; Lee, H.; Dennis, H.; Gladwell, J.; Shurtz, A.; Fullwood, D.T.; Seeley, M.K.; Mitchell, U.H.; et al. Accurate Prediction of Knee Angles during Open-Chain Rehabilitation Exercises Using a Wearable Array of Nanocomposite Stretch Sensors. Sensors 2022, 22, 2499. https://doi.org/10.3390/s22072499

AMA Style

Wood DS, Jensen K, Crane A, Lee H, Dennis H, Gladwell J, Shurtz A, Fullwood DT, Seeley MK, Mitchell UH, et al. Accurate Prediction of Knee Angles during Open-Chain Rehabilitation Exercises Using a Wearable Array of Nanocomposite Stretch Sensors. Sensors. 2022; 22(7):2499. https://doi.org/10.3390/s22072499

Chicago/Turabian Style

Wood, David S., Kurt Jensen, Allison Crane, Hyunwook Lee, Hayden Dennis, Joshua Gladwell, Anne Shurtz, David T. Fullwood, Matthew K. Seeley, Ulrike H. Mitchell, and et al. 2022. "Accurate Prediction of Knee Angles during Open-Chain Rehabilitation Exercises Using a Wearable Array of Nanocomposite Stretch Sensors" Sensors 22, no. 7: 2499. https://doi.org/10.3390/s22072499

APA Style

Wood, D. S., Jensen, K., Crane, A., Lee, H., Dennis, H., Gladwell, J., Shurtz, A., Fullwood, D. T., Seeley, M. K., Mitchell, U. H., Christensen, W. F., & Bowden, A. E. (2022). Accurate Prediction of Knee Angles during Open-Chain Rehabilitation Exercises Using a Wearable Array of Nanocomposite Stretch Sensors. Sensors, 22(7), 2499. https://doi.org/10.3390/s22072499

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop