[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
sensors-logo

Journal Browser

Journal Browser

Biosensors for Gait Measurements and Patient Rehabilitation

A special issue of Sensors (ISSN 1424-8220). This special issue belongs to the section "Biosensors".

Deadline for manuscript submissions: closed (30 April 2024) | Viewed by 22709

Special Issue Editors


E-Mail Website
Guest Editor
Department of Development and Regeneration, KU Leuven, 3000 Leuven, Belgium
Interests: gait; orthopedic surgery; knee; knee joint; kinematics; biomedical imaging; 3D motion analysis; biomechanics

E-Mail Website
Guest Editor
Human Movement Biomechanics Research Group, Department of Biomedical Kinesiology, KU Leuven, Leuven, Belgium
Interests: biomechanics; mechanobiology; cartilage; multi-scale modeling; bioreactor; motion capture; osteoarthritis; physics-based modeling; musculoskeletal loading

Special Issue Information

Dear Colleagues,

Musculoskeletal care is gradually but steadily evolving from a primary focus on structural aspects of the musculoskeletal system towards a more comprehensive approach that integrates and targets optimal function of the musculoskeletal system in every step of the care pathway. In other words, we are moving away from “targeting the best possible X-ray” and towards “targeting the best possible functional outcome”. Nevertheless, this evolution is inherently data-hungry, and its success relies on reliable, patient-friendly technology that allows accurately capturing the key biomechanical drivers for, and biomarkers of, an optimally functioning musculoskeletal system during daily-life activities, including gait as an example. Furthermore, the associated data are preferentially continuously captured in a patient’s natural setting, as recent findings suggest that collecting such data at discrete moments in a hospital or laboratory setting may not accurately reflect the patient’s true functional status.

Therefore, this Special Issue welcomes original research articles, experimental studies, and systematic reviews covering all types of biosensors that contribute to the above evolution and target capturing biomechanical data of the musculoskeletal system during gait and/or other daily-life motor tasks with the goal of better informing musculoskeletal care, from diagnosis, through conservative or surgical treatment and the associated rehabilitation, up to post-treatment follow-up.

Prof. Dr. Lennart Scheys
Prof. Dr. Ilse Jonkers
Guest Editors

Manuscript Submission Information

Manuscripts should be submitted online at www.mdpi.com by registering and logging in to this website. Once you are registered, click here to go to the submission form. Manuscripts can be submitted until the deadline. All submissions that pass pre-check are peer-reviewed. Accepted papers will be published continuously in the journal (as soon as accepted) and will be listed together on the special issue website. Research articles, review articles as well as short communications are invited. For planned papers, a title and short abstract (about 100 words) can be sent to the Editorial Office for announcement on this website.

Submitted manuscripts should not have been published previously, nor be under consideration for publication elsewhere (except conference proceedings papers). All manuscripts are thoroughly refereed through a single-blind peer-review process. A guide for authors and other relevant information for submission of manuscripts is available on the Instructions for Authors page. Sensors is an international peer-reviewed open access semimonthly journal published by MDPI.

Please visit the Instructions for Authors page before submitting a manuscript. The Article Processing Charge (APC) for publication in this open access journal is 2600 CHF (Swiss Francs). Submitted papers should be well formatted and use good English. Authors may use MDPI's English editing service prior to publication or during author revisions.

Benefits of Publishing in a Special Issue

  • Ease of navigation: Grouping papers by topic helps scholars navigate broad scope journals more efficiently.
  • Greater discoverability: Special Issues support the reach and impact of scientific research. Articles in Special Issues are more discoverable and cited more frequently.
  • Expansion of research network: Special Issues facilitate connections among authors, fostering scientific collaborations.
  • External promotion: Articles in Special Issues are often promoted through the journal's social media, increasing their visibility.
  • e-Book format: Special Issues with more than 10 articles can be published as dedicated e-books, ensuring wide and rapid dissemination.

Further information on MDPI's Special Issue polices can be found here.

Published Papers (7 papers)

Order results
Result details
Select all
Export citation of selected articles as:

Research

15 pages, 3498 KiB  
Article
Validation of Inertial-Measurement-Unit-Based Ex Vivo Knee Kinematics during a Loaded Squat before and after Reference-Frame-Orientation Optimisation
by Svenja Sagasser, Adrian Sauer, Christoph Thorwächter, Jana G. Weber, Allan Maas, Matthias Woiczinski, Thomas M. Grupp and Ariana Ortigas-Vásquez
Sensors 2024, 24(11), 3324; https://doi.org/10.3390/s24113324 - 23 May 2024
Cited by 1 | Viewed by 925
Abstract
Recently, inertial measurement units have been gaining popularity as a potential alternative to optical motion capture systems in the analysis of joint kinematics. In a previous study, the accuracy of knee joint angles calculated from inertial data and an extended Kalman filter and [...] Read more.
Recently, inertial measurement units have been gaining popularity as a potential alternative to optical motion capture systems in the analysis of joint kinematics. In a previous study, the accuracy of knee joint angles calculated from inertial data and an extended Kalman filter and smoother algorithm was tested using ground truth data originating from a joint simulator guided by fluoroscopy-based signals. Although high levels of accuracy were achieved, the experimental setup leveraged multiple iterations of the same movement pattern and an absence of soft tissue artefacts. Here, the algorithm is tested against an optical marker-based system in a more challenging setting, with single iterations of a loaded squat cycle simulated on seven cadaveric specimens on a force-controlled knee rig. Prior to the optimisation of local coordinate systems using the REference FRame Alignment MEthod (REFRAME) to account for the effect of differences in local reference frame orientation, root-mean-square errors between the kinematic signals of the inertial and optical systems were as high as 3.8° ± 3.5° for flexion/extension, 20.4° ± 10.0° for abduction/adduction and 8.6° ± 5.7° for external/internal rotation. After REFRAME implementation, however, average root-mean-square errors decreased to 0.9° ± 0.4° and to 1.5° ± 0.7° for abduction/adduction and for external/internal rotation, respectively, with a slight increase to 4.2° ± 3.6° for flexion/extension. While these results demonstrate promising potential in the approach’s ability to estimate knee joint angles during a single loaded squat cycle, they highlight the limiting effects that a reduced number of iterations and the lack of a reliable consistent reference pose inflicts on the sensor fusion algorithm’s performance. They similarly stress the importance of adapting underlying assumptions and correctly tuning filter parameters to ensure satisfactory performance. More importantly, our findings emphasise the notable impact that properly aligning reference-frame orientations before comparing joint kinematics can have on results and the conclusions derived from them. Full article
(This article belongs to the Special Issue Biosensors for Gait Measurements and Patient Rehabilitation)
Show Figures

Figure 1

Figure 1
<p>Cadaveric knee on the knee rig, side-view with optical markers (<b>left</b>) and from behind with IMU sensors (<b>right</b>).</p>
Full article ">Figure 2
<p>Raw knee joint angles, plotted over the entire squat movement, expressed as a percentage. The solid purple lines represent the angles estimated using inertial data. The dashed green lines represent the angles measured by the optical marker-based system. Each row represents one subject, while the columns represent flexion/extension, ab/adduction and ext/internal rotation (from <b>left</b> to <b>right</b>).</p>
Full article ">Figure 3
<p>REFRAME<sub>IMU→GOM</sub> knee joint angles, plotted over the entire squat movement, expressed as a percentage. The solid purple lines represent the angles estimated using inertial data. The dashed green lines represent the angles measured by the optical marker-based system. Each row represents one subject, while the columns represent flexion/extension, ab/adduction and ext/internal rotation (from <b>left</b> to <b>right</b>).</p>
Full article ">Figure 4
<p>REFRAME<sub>RMS</sub> knee joint angles, plotted over the entire squat movement, expressed as a percentage. The solid purple lines represent the angles estimated using inertial data. The dashed green lines represent the angles measured by the optical marker-based system. Each row represents one subject, while the columns represent flexion/extension, ab/adduction and ext/internal rotation (from <b>left</b> to <b>right</b>).</p>
Full article ">Figure 5
<p>Root-mean-square error (RMSE) comparison: average root-mean-square errors ± standard deviation, before optimisation (raw), after REFRAME<sub>IMU→GOM</sub> and REFRAME<sub>RMS</sub>: (<b>a</b>) flexion/extension, (<b>b</b>) abduction/adduction and (<b>c</b>) <span class="html-italic">external/</span>internal rotation. Statistical significance of differences (<span class="html-italic">p</span>-values) based on a paired <span class="html-italic">t</span>-test between raw and REFRAME<sub>IMU→GOM</sub> and raw and REFRAME<sub>RMS</sub>.</p>
Full article ">
12 pages, 1589 KiB  
Article
Four Days Are Enough to Provide a Reliable Daily Step Count in Mild to Moderate Parkinson’s Disease through a Commercial Smartwatch
by Edoardo Bianchini, Silvia Galli, Marika Alborghetti, Lanfranco De Carolis, Alessandro Zampogna, Clint Hansen, Nicolas Vuillerme, Antonio Suppa and Francesco E. Pontieri
Sensors 2023, 23(21), 8971; https://doi.org/10.3390/s23218971 - 4 Nov 2023
Cited by 4 | Viewed by 1663
Abstract
Daily steps could be a valuable indicator of real-world ambulation in Parkinson’s disease (PD). Nonetheless, no study to date has investigated the minimum number of days required to reliably estimate the average daily steps through commercial smartwatches in people with PD. Fifty-six patients [...] Read more.
Daily steps could be a valuable indicator of real-world ambulation in Parkinson’s disease (PD). Nonetheless, no study to date has investigated the minimum number of days required to reliably estimate the average daily steps through commercial smartwatches in people with PD. Fifty-six patients were monitored through a commercial smartwatch for 5 consecutive days. The total daily steps for each day was recorded and the average daily steps was calculated as well as the working and weekend days average steps. The intraclass correlation coefficient (ICC) (3,k), standard error of measurement (SEM), Bland–Altman statistics, and minimum detectable change (MDC) were used to evaluate the reliability of the step count for every combination of 2–5 days. The threshold for acceptability was set at an ICC ≥ 0.8 with a lower bound of CI 95% ≥ 0.75 and a SAM < 10%. ANOVA and Mann–Whitney tests were used to compare steps across the days and between the working and weekend days, respectively. Four days were needed to achieve an acceptable reliability (ICC range: 0.84–0.90; SAM range: 7.8–9.4%). In addition, daily steps did not significantly differ across the days and between the working and weekend days. These findings could support the use of step count as a walking activity index and could be relevant to developing monitoring, preventive, and rehabilitation strategies for people with PD. Full article
(This article belongs to the Special Issue Biosensors for Gait Measurements and Patient Rehabilitation)
Show Figures

Figure 1

Figure 1
<p>Boxplots showing daily steps for each monitoring day in PD patients. Thick line in the boxes indicates the median; lower and upper box limits indicate first (Q1) and third quartile (Q3), respectively; black vertical lines indicate lower and upper outliers boundaries calculated as Q1 − (1.5 × IQR) and Q3 + (1.5 × IQR), respectively. Red X indicates mean values for each disease stage. PD: Parkinson’s disease.</p>
Full article ">Figure 2
<p>Boxplots showing daily steps for working and weekend days in PD patients. The thick line in the boxes indicates the median; lower and upper box limits indicate first (Q1) and third quartile (Q3), respectively; black vertical lines indicate lower and upper outliers boundaries calculated as Q1 − (1.5 × IQR) and Q3 + (1.5 × IQR), respectively. Red X indicates mean values for each disease stage. PD: Parkinson’s disease.</p>
Full article ">Figure 3
<p>Graph showing the ICC(3,k) values across the different combinations of monitoring days. The top dark blue line represents 95% CI upper bound. The middle blue line represents the ICC (3,k) point estimate. The bottom light blue line represents 95% CI lower bound. Vertical brown lines were added to mark 2-day, 3-day, 4, day, and 5-day combinations. Dashed horizontal dark and light grey lines were added to mark the 0.8 threshold for the point estimate and the 0.75 threshold for 95% CI lower bound, respectively. CI: confidence interval; ICC: intraclass correlation coefficient. LB: lower bound; UB: upper bound.</p>
Full article ">
14 pages, 1713 KiB  
Article
Sensitivity of Model-Based Predictions of Post-TKA Kinematic Behavior to Residual Errors in Ultrasound-Based Knee Collateral Ligament Strain Assessment
by Félix Dandois, Orçun Taylan, Jacobus H. Müller and Lennart Scheys
Sensors 2023, 23(19), 8268; https://doi.org/10.3390/s23198268 - 6 Oct 2023
Viewed by 1249
Abstract
Ultrasound-based ligament strain estimation shows promise in non-invasively assessing knee joint collateral ligament behavior and improving ligament balancing procedures. However, the impact of ultrasound-based strain estimation residual errors on in-silico arthroplasty predictions remains unexplored. We investigated the sensitivity of post-arthroplasty kinematic predictions to [...] Read more.
Ultrasound-based ligament strain estimation shows promise in non-invasively assessing knee joint collateral ligament behavior and improving ligament balancing procedures. However, the impact of ultrasound-based strain estimation residual errors on in-silico arthroplasty predictions remains unexplored. We investigated the sensitivity of post-arthroplasty kinematic predictions to ultrasound-based strain estimation errors compared to clinical inaccuracies in implant positioning.Two cadaveric legs were submitted to active squatting, and specimen-specific rigid computer models were formulated. Mechanical properties of the ligament model were optimized to reproduce experimentally obtained tibiofemoral kinematics and loads with minimal error. Resulting remaining errors were comparable to the current state-of-the-art. Ultrasound-derived strain residual errors were then introduced by perturbing lateral collateral ligament (LCL) and medial collateral ligament (MCL) stiffness. Afterwards, the implant position was perturbed to match with the current clinical inaccuracies reported in the literature. Finally, the impact on simulated post-arthroplasty tibiofemoral kinematics was compared for both perturbation scenarios. Ultrasound-based errors minimally affected kinematic outcomes (mean differences < 0.73° in rotations, 0.1 mm in translations). Greatest differences occurred in external tibial rotations (−0.61° to 0.73° for MCL, −0.28° to 0.27° for LCL). Comparatively, changes in implant position had larger effects, with mean differences up to 1.95° in external tibial rotation and 0.7 mm in mediolateral translation. In conclusion, our study demonstrated that the ultrasound-based assessment of collateral ligament strains has the potential to enhance current computer-based pre-operative knee arthroplasty planning. Full article
(This article belongs to the Special Issue Biosensors for Gait Measurements and Patient Rehabilitation)
Show Figures

Figure 1

Figure 1
<p>Virtual knee-joint simulator. The femur is fixed to the hip assembly free to rotate in the sagittal plane and to translate along the coronal axis). The tibia is linked to the ankle assembly with a cylindric joint (allowing rotation in the coronal plane). The ankle assembly is free to rotate in all planes and translate along the sagittal axis. Tibiofemoral and patellofemoral joints are created using a contact function constraining the joints based on the bone and implant geometries. Additionally, the model is also constrained by active loading, i.e., quadriceps load between patella and hip assembly (3 bundles) and hamstrings load between tibia and hip assembly (1 bundle medial and 1 bundle lateral) as well as passive forces from soft tissues, i.e., patellar ligament between the tibia and the patella (3 bundles) and both collateral ligaments: MCL (2 bundles) and LCL (1 bundle).</p>
Full article ">Figure 2
<p>Tibiofemoral kinematics of both specimens (Specimen 1 on <b>top</b> and Specimen 2 on <b>bottom</b>) experimentally obtained (blue) and simulated with the developed model (red). The displayed kinematics are, from left to right, adduction, external tibial rotation, anterior–posterior (AP) translation, medio-lateral (ML) translation and inferior–superior (IS) translation.</p>
Full article ">Figure 3
<p>Comparisons of simulated adduction and external tibial rotation obtained after modification of ligaments stiffness (magenta and green) and results obtained after modification of implant internal rotation (blue and red) across the flexion range (dashed orange) with optimal kinematics (black). Parts (<b>A</b>,<b>B</b>) display comparisons of Specimen 1 for MCL and LCL, respectively. Part (<b>C</b>) displays comparisons of Specimen 2.</p>
Full article ">Figure 3 Cont.
<p>Comparisons of simulated adduction and external tibial rotation obtained after modification of ligaments stiffness (magenta and green) and results obtained after modification of implant internal rotation (blue and red) across the flexion range (dashed orange) with optimal kinematics (black). Parts (<b>A</b>,<b>B</b>) display comparisons of Specimen 1 for MCL and LCL, respectively. Part (<b>C</b>) displays comparisons of Specimen 2.</p>
Full article ">
15 pages, 1712 KiB  
Article
The Validity of Hawkin Dynamics Wireless Dual Force Plates for Measuring Countermovement Jump and Drop Jump Variables
by Andrew J. Badby, Peter D. Mundy, Paul Comfort, Jason P. Lake and John J. McMahon
Sensors 2023, 23(10), 4820; https://doi.org/10.3390/s23104820 - 17 May 2023
Cited by 21 | Viewed by 7621
Abstract
Force plate testing is becoming more commonplace in sport due to the advent of commercially available, portable, and affordable force plate systems (i.e., hardware and software). Following the validation of the Hawkin Dynamics Inc. (HD) proprietary software in recent literature, the aim of [...] Read more.
Force plate testing is becoming more commonplace in sport due to the advent of commercially available, portable, and affordable force plate systems (i.e., hardware and software). Following the validation of the Hawkin Dynamics Inc. (HD) proprietary software in recent literature, the aim of this study was to determine the concurrent validity of the HD wireless dual force plate hardware for assessing vertical jumps. During a single testing session, the HD force plates were placed directly atop two adjacent Advanced Mechanical Technology Inc. in-ground force plates (the “gold standard”) to simultaneously collect vertical ground reaction forces produced by 20 participants (27 ± 6 years, 85 ± 14 kg, 176.5 ± 9.23 cm) during the countermovement jump (CMJ) and drop jump (DJ) tests (1000 Hz). Agreement between force plate systems was determined via ordinary least products regression using bootstrapped 95% confidence intervals. No bias was present between the two force plate systems for any of the CMJ and DJ variables, except DJ peak braking force (proportional bias) and DJ peak braking power (fixed and proportional bias). The HD system may be considered a valid alternative to the industry gold standard for assessing vertical jumps because fixed or proportional bias was identified for none of the CMJ variables (n = 17) and only 2 out of 18 DJ variables. Full article
(This article belongs to the Special Issue Biosensors for Gait Measurements and Patient Rehabilitation)
Show Figures

Figure 1

Figure 1
<p>Example set-up for data collection (frontal plane).</p>
Full article ">Figure 2
<p>Example set-up for data collection (sagittal plane).</p>
Full article ">Figure 3
<p>A representative example of an original (<b>top</b>) and time-aligned (<b>bottom</b>) countermovement jump trial recorded by the AMTI (solid grey line) and HD (dotted black line) force plate systems. The bottom graph also illustrates the occurrence of key events. vGRF = vertical ground reaction force.</p>
Full article ">Figure 4
<p>A representative example of an original (<b>top</b>) and time-aligned (<b>bottom</b>) drop jump trial recorded by the AMTI (solid grey line) and HD (dotted black line) force plate systems. The bottom graph also illustrates the occurrence of key events. vGRF = vertical ground reaction force.</p>
Full article ">
12 pages, 1900 KiB  
Article
Peak Tibiofemoral Contact Forces Estimated Using IMU-Based Approaches Are Not Significantly Different from Motion Capture-Based Estimations in Patients with Knee Osteoarthritis
by Giacomo Di Raimondo, Miel Willems, Bryce Adrian Killen, Sara Havashinezhadian, Katia Turcot, Benedicte Vanwanseele and Ilse Jonkers
Sensors 2023, 23(9), 4484; https://doi.org/10.3390/s23094484 - 4 May 2023
Cited by 7 | Viewed by 3619
Abstract
Altered tibiofemoral contact forces represent a risk factor for osteoarthritis onset and progression, making optimization of the knee force distribution a target of treatment strategies. Musculoskeletal model-based simulations are a state-of-the-art method to estimate joint contact forces, but they typically require laboratory-based input [...] Read more.
Altered tibiofemoral contact forces represent a risk factor for osteoarthritis onset and progression, making optimization of the knee force distribution a target of treatment strategies. Musculoskeletal model-based simulations are a state-of-the-art method to estimate joint contact forces, but they typically require laboratory-based input and skilled operators. To overcome these limitations, ambulatory methods, relying on inertial measurement units, have been proposed to estimated ground reaction forces and, consequently, knee contact forces out-of-the-lab. This study proposes the use of a full inertial-capture-based musculoskeletal modelling workflow with an underlying probabilistic principal component analysis model trained on 1787 gait cycles in patients with knee osteoarthritis. As validation, five patients with knee osteoarthritis were instrumented with 17 inertial measurement units and 76 opto-reflective markers. Participants performed multiple overground walking trials while motion and inertial capture methods were synchronously recorded. Moderate to strong correlations were found for the inertial capture-based knee contact forces compared to motion capture with root mean square error between 0.15 and 0.40 of body weight. The results show that our workflow can inform and potentially assist clinical practitioners to monitor knee joint loading in physical therapy sessions and eventually assess long-term therapeutic effects in a clinical context. Full article
(This article belongs to the Special Issue Biosensors for Gait Measurements and Patient Rehabilitation)
Show Figures

Figure 1

Figure 1
<p>InCap-based workflow overview. (1) OpenIMUs—IMU calibration and inverse kinematics estimation; (2) PPCA model—estimation of GRFM; (3) ZMP—estimation of COP; (4) OpenSim JAM—KCF estimation and comparison.</p>
Full article ">Figure 2
<p>Measured (blue) vs. estimated (red) ground reaction forces (GRF) in BW (body weight) and ground reaction moments (GRM) in Nm/kg of the PPCA KOA population-based model (1787 gait cycles). Note that the illustrated comparison refers to the OA-involved knee joint. The shaded areas represent the standard deviation of the mean.</p>
Full article ">Figure 3
<p>Mean measured (blue solid) vs. mean estimated (orange solid) COP trajectories of one patient during walking (84 gait cycles) (COP trajectories in dashed—measured in blue and estimated in orange).</p>
Full article ">Figure 4
<p>MoCap-based (blue) vs. InCap-based (red) knee contact forces in BW (body weight) for the medial and lateral knee compartment of the OA-involved limb. The shaded areas represent the standard deviation of the mean.</p>
Full article ">
16 pages, 1606 KiB  
Article
Video-based Goniometer Applications for Measuring Knee Joint Angles during Walking in Neurological Patients: A Validity, Reliability and Usability Study
by Monica Parati, Matteo Gallotta, Beatrice De Maria, Annalisa Pirola, Matteo Morini, Luca Longoni, Emilia Ambrosini, Giorgio Ferriero and Simona Ferrante
Sensors 2023, 23(4), 2232; https://doi.org/10.3390/s23042232 - 16 Feb 2023
Cited by 2 | Viewed by 3691
Abstract
Easy-to-use evaluation of Range Of Motion (ROM) during walking is necessary to make decisions during neurological rehabilitation programs and during follow-up visits in clinical and remote settings. This study discussed goniometer applications (DrGoniometer and Angles - Video Goniometer) that measure knee joint ROM [...] Read more.
Easy-to-use evaluation of Range Of Motion (ROM) during walking is necessary to make decisions during neurological rehabilitation programs and during follow-up visits in clinical and remote settings. This study discussed goniometer applications (DrGoniometer and Angles - Video Goniometer) that measure knee joint ROM during walking through smartphone cameras. The primary aim of the study is to test the inter-rater and intra-rater reliability of the collected measurements as well as their concurrent validity with an electro-goniometer. The secondary aim is to evaluate the usability of the two mobile applications. A total of 22 patients with Parkinson’s disease (18 males, age 72 (8) years), 22 post-stroke patients (17 males, age 61 (13) years), and as many healthy volunteers (8 males, age 45 (5) years) underwent knee joint ROM evaluations during walking. Clinicians and inexperienced examiners used the two mobile applications to calculate the ROM, and then rated their perceived usability through the System Usability Scale (SUS). Intraclass correlation coefficients (ICC) and correlation coefficients (corr) were calculated. Both applications showed good reliability (ICC > 0.69) and validity (corr > 0.61), and acceptable usability (SUS > 68). Smartphone-based video goniometers could be used to assess the knee ROM during walking in neurological patients, because of their acceptable degree of reliability, validity and usability. Full article
(This article belongs to the Special Issue Biosensors for Gait Measurements and Patient Rehabilitation)
Show Figures

Figure 1

Figure 1
<p>Screenshot images of the current version of Angles (<b>a</b>) and the DrGoniometer apps (<b>b</b>), and a schematic representation of the experimental study design composed of a testing patient’s session (T0) and two assessment sessions (T1, T2). The first assessment session (T1) was completed by two experienced examiners (E1, E2) and two inexperienced examiners (I1, I2). The assessments were repeated at T2 by one experienced examiner (E1) and one inexperienced examiner (I1) (<b>c</b>).</p>
Full article ">Figure 2
<p>Bland-Altman plot of the knee range of motion of the three participant populations (SK: post-stroke patients, PD: patients with Parkinson’s disease, HC: healthy controls) collected by the experienced examiner (E1) using the Angles and the DrGoniometer applications. The blue line indicates the mean differences between the mobile app and electro-goniometer measures, and the red lines are the 95% lower and upper limits of agreement.</p>
Full article ">Figure 3
<p>Bland-Altman plot of the knee range of motion of the three participant populations (SK: post-stroke patients, PD: patients with Parkinson’s disease, HC: healthy controls) collected by the inexperienced examiner (I1) using the Angles and DrGoniometer applications. The blue line indicates the mean differences between the mobile app and electro-goniometer measures, and the red lines are the 95% lower and upper limits of agreement.</p>
Full article ">Figure 4
<p>Usability findings of the Angles app (<b>a</b>) and DrGoniometer app (<b>b</b>) evaluated with the system usability scale (SUS) by experienced and inexperienced raters, respectively, in yellow and blue. The mean score collected by experienced and inexperienced raters, respectively, in yellow and blue, was reported for the ten SUS items (1: totally disagree; 5: totally agree). The green stars indicated the best possible usability results.</p>
Full article ">
34 pages, 2005 KiB  
Article
Detecting Gait Events from Accelerations Using Reservoir Computing
by Laurent Chiasson-Poirier, Hananeh Younesian, Katia Turcot and Julien Sylvestre
Sensors 2022, 22(19), 7180; https://doi.org/10.3390/s22197180 - 21 Sep 2022
Cited by 7 | Viewed by 2555
Abstract
Segmenting the gait cycle into multiple phases using gait event detection (GED) is a well-researched subject with many accurate algorithms. However, the algorithms that are able to perform accurate and robust GED for real-life environments and physical diseases tend to be too complex [...] Read more.
Segmenting the gait cycle into multiple phases using gait event detection (GED) is a well-researched subject with many accurate algorithms. However, the algorithms that are able to perform accurate and robust GED for real-life environments and physical diseases tend to be too complex for their implementation on simple hardware systems limited in computing power and memory, such as those used in wearable devices. This study focuses on a numerical implementation of a reservoir computing (RC) algorithm called the echo state network (ESN) that is based on simple computational steps that are easy to implement on portable hardware systems for real-time detection. RC is a neural network method that is widely used for signal processing applications and uses a fast-training method based on a ridge regression adapted to the large quantity and variety of IMU data needed to use RC in various real-life environment GED. In this study, an ESN was used to perform offline GED with gait data from IMU and ground force sensors retrieved from three databases for a total of 28 healthy adults and 15 walking conditions. Our main finding is that despite its low complexity, ESN is robust for GED, with performance comparable to other state-of-the-art algorithms. Our results show the ESN is robust enough to obtain good detection results in all conditions if the algorithm is trained with variable data that match those conditions. The distribution of the mean absolute errors (MAE) between the detection times from the ESN and the force sensors were between 40 and 120 ms for 6 defined gait events (95th percentile). We compared our ESN with four different state-of-the-art algorithms from the literature. The ESN obtained a MAE not more than 10 ms above three other reference algorithms for normal walking indoor and outdoor conditions and yielded the 2nd lowest MAE and the 2nd highest true positive rate and specificity when applied to outdoor walking and running conditions. Our work opens the door to using the ESN as a GED for applications in wearable sensors for long-term patient monitoring. Full article
(This article belongs to the Special Issue Biosensors for Gait Measurements and Patient Rehabilitation)
Show Figures

Figure 1

Figure 1
<p>Procedure used to evaluate the performance of the ESN and the TKEO algorithms. Arrow boxes represent calculation steps and square boxes represent data. The experimental data were collected during walking tests as explained in <a href="#sec2dot1-sensors-22-07180" class="html-sec">Section 2.1</a> and preprocessed as explained in <a href="#sec2dot2-sensors-22-07180" class="html-sec">Section 2.2</a>. The input timeseries <math display="inline"><semantics> <msub> <mi mathvariant="bold">U</mi> <mi>k</mi> </msub> </semantics></math> is acceleration signals recorded on the participant left foot. The output timeseries <math display="inline"><semantics> <msubsup> <mi mathvariant="bold">Y</mi> <mrow> <mi>k</mi> </mrow> <mo>′</mo> </msubsup> </semantics></math> is the gait event targets identified with ground force sensors and formatted as a multi-channel binary signal. Each channel indicated the event time indices of one gait event class. The timeseries lists of the inputs <math display="inline"><semantics> <mrow> <mo>{</mo> <msub> <mi mathvariant="bold">U</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>…</mo> <mo>,</mo> <msub> <mi mathvariant="bold">U</mi> <mi>k</mi> </msub> <mo>,</mo> <mo>…</mo> <mo>,</mo> <msub> <mi mathvariant="bold">U</mi> <msub> <mi mathvariant="normal">N</mi> <mi>TS</mi> </msub> </msub> <mo>}</mo> </mrow> </semantics></math>, and the outputs <math display="inline"><semantics> <mrow> <mo>{</mo> <msubsup> <mi mathvariant="bold">Y</mi> <mn>1</mn> <mo>′</mo> </msubsup> <mo>,</mo> <mo>…</mo> <mo>,</mo> <msubsup> <mi mathvariant="bold">Y</mi> <mi>k</mi> <mo>′</mo> </msubsup> <mo>,</mo> <mo>…</mo> <mo>,</mo> <msubsup> <mi mathvariant="bold">Y</mi> <mrow> <msub> <mi mathvariant="normal">N</mi> <mi>TS</mi> </msub> </mrow> <mo>′</mo> </msubsup> <mo>}</mo> </mrow> </semantics></math> represent a subset of multiple input and output timeseries for each participant and walking condition considered in the procedure. The timeseries are separated in two groups, one for the training step of the ESN algorithm and one for the testing step of both ESN and TKEO algorithms. The ESN algorithm steps are represented by the red boxes and explained in <a href="#sec2dot3-sensors-22-07180" class="html-sec">Section 2.3</a>. The TKEO algorithm steps are represented by the blue boxes and explained in <a href="#sec2dot4-sensors-22-07180" class="html-sec">Section 2.4</a>. The green boxes evaluate the error between predicted and target events and establish performance criteria as presented in <a href="#sec2dot5-sensors-22-07180" class="html-sec">Section 2.5</a>.</p>
Full article ">Figure 2
<p>IMU and FSR positioning for the three databases. In CIRRIS tests, the IMU was placed outside the shoes using straps. For MAREA tests, the three IMU were attached with straps at the three positions shown, and the FSR were placed into the shoe soles. In UDS tests, the IMU and the FSR were placed on a sock inside the shoes. The orientations are defined as the antero-posterior (AP), the medio-lateral (ML) and the vertical (V) axis.</p>
Full article ">Figure 3
<p>Walking tracks of UDS walking tests. The subjects were walking continuously on the full arrow and coming back following the dotted arrow. T1 was a normal straight walk without obstacles. T2 was a walking circuit over small obstacles, including stepping on one 0.3 m wide by 1 m long by 0.2 m high hard platform (O1), stepping over 0.1 m high by 0.3 m square side platforms (O2) and stepping on 0.5 m wide by 1 m long by 4 cm thick exercise mats (O3). The participant stepped on O1 and O3 by putting their feet on it but had to pass over O2 without touching it. T3 was a walking circuit around 0.3 m wide by 1 m long platforms (O1) aligned in a row with 2 m distance between platforms.</p>
Full article ">Figure 4
<p>Event target identification process from foot pressure measurements for each three databases. Examples of foot pressure measurements, acceleration measurements, and target gait events identified in two gait cycles for each database are shown in subfigures (<b>b</b>–<b>d</b>). The heel strike (HS) is the first contact of the foot with the ground. The heel push (HP) is the maximum load response of the foot after the HS. The foot flat (FF) is the minimum load response of the foot in the middle of the stance phase. The heel off (HO) is the end of the heel contact on the ground in the second part of the stance phase. The toe push (TP) is the second maximum load response of the foot in the second part of the support phase. The toe off (TO) is the final contact of the foot on the ground. (<b>a</b>) Event identification; (<b>b</b>) MAREA database (acceleration of the foot position); (<b>c</b>) CIRRIS database; (<b>d</b>) UDS database. FSR1 is on the toe and FSR2 is on the heel.</p>
Full article ">Figure 5
<p>ESN algorithm representation considering the HS and TO gait event classes for the output prediction.</p>
Full article ">Figure 6
<p>Average (dots) and ±1 standard deviation (error bars) of the MAE for the values reported in reference [<a href="#B10-sensors-22-07180" class="html-bibr">10</a>] and our own implementation of the TKEO algorithm, under conditions of treadmill walking (TW), treadmill incline (TI) and indoor walking (IW).</p>
Full article ">Figure 7
<p>NRMSE and MAE as a function of the length of the training set (with the AP axis acceleration input) considering all events for each database. The errors were evaluated using the <math display="inline"><semantics> <msub> <mi>MAR</mi> <mi>all</mi> </msub> </semantics></math>, <math display="inline"><semantics> <msub> <mi>CIR</mi> <mi>all</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>UDS</mi> <mi>all</mi> </msub> </semantics></math> subsets. The three subsets were considered together in the set ALL, with the HS and TO gait event types only. The MAE and NRMSE are different for all databases, mostly because of the harder and higher number of gait event types detected in the CIRRIS and UDS databases than in the MAREA.</p>
Full article ">Figure 8
<p>MAE (ms) of the prediction with the ESN on HS and TO events as a function of the training/testing database combination considering all walking condition subsets <math display="inline"><semantics> <msub> <mi>MAR</mi> <mi>all</mi> </msub> </semantics></math>, <math display="inline"><semantics> <msub> <mi>CIR</mi> <mi>all</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>UDS</mi> <mi>all</mi> </msub> </semantics></math>.</p>
Full article ">Figure 9
<p>MAE, FPR, PPV and T1 score of the ESN as a function of the subsets of walking condition. (<b>a</b>) MAREA database; (<b>b</b>) CIRRIS database; (<b>c</b>) UDS database. The mean value of MEA is expressed as a green triangle and the median as a central line, while boxes are the 25th and 75th percentiles and the whiskers are the extensions of the last value of <math display="inline"><semantics> <mrow> <mi>MAE</mi> <msub> <mrow/> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> </mrow> </semantics></math> of the distribution over and under 1.5 times the interquartile range. The values of TPR, PPV and <span class="html-italic">T</span><math display="inline"><semantics> <msub> <mrow/> <mn>1</mn> </msub> </semantics></math> are plotted with three lines and are relative to the right-side axis. (Same for all boxplot).</p>
Full article ">Figure 10
<p>MAREA database results of MAE, TPV, PPV and T1 performance of the ESN for GED on subset <math display="inline"><semantics> <msub> <mi>MAR</mi> <mi>all</mi> </msub> </semantics></math>. The boxplots are related to the left axis and represent the MAE distribution over all the timeseries of the testing set. The lines, relative to the right axis, represent the value of the TPV, PPV and T1 score. The horizontal axis categories are relative to different input axis combinaisons. Each vertically aligned subplot represents the results for one event class. The three subfigures represent the result from different IMU positions. (<b>a</b>) Foot sensor position; (<b>b</b>) wrist sensor position; (<b>c</b>) waist sensor position.</p>
Full article ">Figure 11
<p>Same as <a href="#sensors-22-07180-f010" class="html-fig">Figure 10</a> for the CIRRIS database. The ESN for GED on subset <math display="inline"><semantics> <msub> <mi>CIR</mi> <mi>all</mi> </msub> </semantics></math> with foot IMU position.</p>
Full article ">Figure 12
<p>Same as <a href="#sensors-22-07180-f010" class="html-fig">Figure 10</a> for the UDS database. The ESN for GED on subset <math display="inline"><semantics> <msub> <mi>UDS</mi> <mi>all</mi> </msub> </semantics></math> with foot IMU position.</p>
Full article ">Figure 13
<p>Comparison of MEA, FPR, PPV and T1 score on the three databases for the ESN and the TKEO algorithms.</p>
Full article ">Figure 14
<p>Comparison of the TKEO and ESN algorithms to other algorithms in terms of the TPR, the <span class="html-italic">SPF</span> and the mean (central point) and ±1 standard deviation (error bar) of the MAE. The subsets <math display="inline"><semantics> <msub> <mi>MAR</mi> <mi>IW</mi> </msub> </semantics></math>, <math display="inline"><semantics> <msub> <mi>MAR</mi> <mi>OW</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>MAR</mi> <mi>OWR</mi> </msub> </semantics></math> of the MAREA database were used, as were the results of the Kh, Au and Ru algorithms from [<a href="#B15-sensors-22-07180" class="html-bibr">15</a>]. The performance criteria used here for TKEO and ESN were adjusted to match the criteria used for Kh, Au and Ru. A threshold of 40 ms was used to consider the events as TP. The MAE was evaluated with the error from the TP events only. The number of points in each timeseries was considered for a sampling rate of 128 Hz.</p>
Full article ">Figure A1
<p>Relationship between the CHARC metrics and the mean absolute error for all events. Exploration was done on the CIRRIS database with timeseries of 3 conditions: normal, right toe out and right side trunk leaning. Five different walking speeds and both subjects were selected (30 timeseries).</p>
Full article ">Figure A2
<p>Performance according to the temporal threshold for peak detection (the shaded areas represent the variability over all the datasets and event classes). The lower, middle and upper line of the shaded areas correspond, respectively, to the 10th, 50th and 90th percentiles of the metric.</p>
Full article ">Figure A3
<p>MAE, TPR, PPV and T1 comparison between the ESN and TKEO algorithms for HS and TO events for multiple walking condition subsets. (<b>a</b>) MAREA database; (<b>b</b>) CIRRIS database; (<b>c</b>) CIRRIS database (zoom on ESN results); (<b>d</b>) UDS database; (<b>e</b>) UDS database (zoom on ESN results).</p>
Full article ">Figure A4
<p>Examples of predictions on a timeseries for the UDS database for the ESN and the TKEO. Four gait cycles are represented. <span class="html-italic">U</span>(<span class="html-italic">n</span>) is the AP acceleration. The other subplots show the output signals computed from the TKEO (blue lines), <math display="inline"><semantics> <mi>ϕ</mi> </semantics></math>(<span class="html-italic">n</span>) and <math display="inline"><semantics> <mi>χ</mi> </semantics></math>(<span class="html-italic">n</span>), or the ESN (red lines), <math display="inline"><semantics> <msub> <mi>Y</mi> <mn>1</mn> </msub> </semantics></math>(<span class="html-italic">n</span>) and <math display="inline"><semantics> <msub> <mi>Y</mi> <mn>2</mn> </msub> </semantics></math>(<span class="html-italic">n</span>). The events predicted (dots) were identified from the peak finder algorithm applied on each output signal. All the events were correctly predicted by the ESN. For the TKEO, the HS of cycles 1–2 were correctly predicted. HS and TO events of cycles 3–4 were misdetected (×), and the TO of cycle 2 was not detected.</p>
Full article ">
Back to TopTop