Abstract
Background
Despite its clinical relevance, there is a relative lack of research examining flexibility and stability based on the acceleration or speed of localized limb segments, such as the shank and thigh.
Objective
This study aimed to evaluate gait characteristics based on acceleration in the thigh and shank to identify differences between the affected and unaffected sides in stroke hemiplegic patients.
Method
Forty individuals with stroke-induced hemiplegia were assessed during a 5-meter walk using a 3D motion analysis system and Inertial Measurement Units (IMUs). Spatial-temporal and acceleration parameters were calculated.
Results
Significant differences were observed between the affected and unaffected sides in stance time, swing time, swing phase, and stance phase. In terms of acceleration, the mean acceleration in the anterior-posterior (AP) direction of the thigh and the mean value of the center of mass (CoM) in the AP direction differed significantly. These spatial-temporal findings were consistent with known characteristics of hemiplegic gait. A notable posterior shift of the thigh CoM on the affected side was identified, likely reflecting impaired propulsion and reduced stability.
Conclusion
The posterior displacement of the thigh CoM on the affected side may represent a compensatory mechanism to maintain balance during gait. Clinically, this posterior CoM shift could serve as a meaningful indicator of hemiplegic gait and a potential target for rehabilitation interventions aimed at restoring gait symmetry and improving functional mobility.
Introduction
Stroke occurs when a blood vessel in the brain is blocked or ruptures, and it is one of the leading causes of death. 1 The incidence of stroke has increased, especially as society has aged. Each year, approximately 15 million people suffer from stroke, and about 5 million of them experience disabilities. 2 When a stroke occurs, damage to peripheral nerves and reduced or abnormal signal transmission from the brain lead to impaired muscle activation. 3 Particularly, strokes caused by thrombosis or hemorrhage result in damage to parts of the brain, causing paralysis on one side of the body. The severity and pattern of hemiplegia vary depending on which part of the brain is affected. The main symptoms of hemiplegia include paralysis on one side of the body, decreased sensation, impaired motor ability, and sometimes speech and cognitive dysfunctions. Stroke also significantly impacts gait, one of the most essential functions for daily life, as the paralyzed leg cannot move properly, leading to difficulties in performing normal walking patterns. 4 This is due to abnormal muscle contractions or the occurrence of spasticity. 5 Typically, the characteristics of hemiplegic gait include reduced walking speed and step length, as well as the occurrence of circumduction gait and foot drop due to abnormal gait event ratios. 6 The characteristics of gait on the affected side differ significantly from the un-affected side, and the affected side also greatly influences the un-affected side. Therefore, an accurate diagnosis of the severity and characteristics of the affected side is crucial for gait pattern improvement and muscle strength enhancement through gait training during rehabilitation. Various previous studies have identified the characteristics of the affected side in hemiplegic gait. Some study found that during the swing phase of gait, the knee on the affected side was lifted higher due to compensatory movements, as measured by a force plate. 7 Another study confirmed asymmetry in the affected side of hemiplegic patients, showing reduced step length, walking speed, and ground contact time. 8
The method of evaluating affected side characteristics using spatial-temporal parameters is commonly used. 9 Additionally, various approaches such as comparing muscle activation through electromyography measurements and utilizing kinetic analysis with force plates have also been proposed.10,11 Studies comparing the gait of the un-affected and affected sides after stroke play a crucial role in understanding the impact of neurological damage on gait and in setting rehabilitation goals. 12 These studies analyze the differences in muscle strength, gait patterns, and balance between the affected and un-affected sides, identifying the asymmetrical gait characteristics of stroke patients and suggesting methods to improve them. The parameters used for kinematic analysis are based on gait events, defining repetitive sections of the gait cycle to calculate spatial-temporal variables. For such analyses, optical 3D motion analysis camera systems based on infrared cameras or gait analysis systems based on force plates are commonly employed. Most measurement systems are relatively expensive and have spatial limitations, making them somewhat difficult to use in general hospitals and similar settings. To address these drawbacks, relatively simple IMU sensors have been used more recently. These sensors are attached to the limbs based on the joints of the lower extremities, and the gyroscope calculates the tilting angles or derives gait events. 13 While traditional gait measurement systems can calculate overall spatial-temporal variables of gait, using IMU sensors allows not only the extraction of events and angle calculations but also the direct or indirect calculation of the speed and acceleration of the attached body parts. 14 This enables the identification of dynamic characteristics at each segment. There have been numerous previous studies using IMU sensors attached to the trunk to examine tilting angles and left-right asymmetry.15,16 However, there is a relative lack of research evaluating flexibility and stability based on the speed or acceleration of localized limb segments such as the shank and thigh. Measuring the acceleration of shank movement during gait plays an important role in the gait analysis of stroke patients. 17 The acceleration of the shank is controlled by the multi-joints of the knee and ankle and is adjusted by the coordinated activation of muscles such as the tibialis anterior and hamstrings.18,19 Additionally, the movement of the thigh is related to the body's balance and stability, and can be influenced by the gluteus muscles, iliopsoas, and quadriceps muscles, which affect the pelvic girdle and tilting. By understanding and predicting the rate of segmental speed change per unit time due to these factors, it can be used to assess the efficiency and safety of movement.
In this study, we aim to evaluate gait characteristics based on acceleration in the thigh and shank of stroke hemiplegic patients to identify the differences between the un-affected and affected sides at local positions. The purpose of this research is to derive new gait evaluation indicators based on the results and to facilitate the development of new gait rehabilitation strategies.
Materials and method
Participants
Outpatients or inpatients with hemiparesis due to stroke from Wonkwang University Gwangju Hospital (WKUGH) were selected. The inclusion criteria for participants were: 1) hemiplegia diagnosed with ischemic stroke within the last 10 years, 2) a functional ambulation category of level 3 or higher, and 3) low skin sensitivity to reduce the likelihood of issues arising from sensor attachment. A total of 40 participants were selected (Table 1). All participants were fully informed about the procedures prior to the experiment and signed an informed consent form. The study was conducted in accordance with the procedures approved by the Institutional Review Board (IRB: WKIRB 2022-07, June 29, 2022).
Data acquisition
Before the gait experiment, participants were thoroughly informed about the test protocol and practiced sufficiently to ensure they could walk as naturally as possible. The gait experiment was performed three times over a 5-meter distance at each participant's preferred walking speed in a straight line. 20 To collect kinematic data, an 8-camera optical infrared system for 3D motion analysis (BTS SMART DX, BTS Bioengineering Corp., Milanese, MI, Italy) was synchronized with a 4-channel inertial sensor system (Delsys Trigno Avanti wireless system, Natick, MA, USA). The sampling frequencies for the motion cameras and IMUs were set at 120 Hz and 1200 Hz, respectively, and the IMU's accelerometer and gyroscope data were resampled and synchronized with the motion data. A total of 18 reflective markers were attached according to the Helen Hayes lower limb set for the 3D motion analysis system, and 4 IMU sensors were attached to the front of the shank and thigh. 21
Gait measures
The average value of data measured three times while walking 5-meter in a straight line at the preferred walking speed was used. Motion analysis software was used to perform the calculations from 3D motion data (BTS Motion Analysis Lab., BTS Bioengineering Corp., Milanese, MI, Italy). The acceleration values of the shank and thigh on the un-affected side and the affected side measured with the IMU. The MATLAB R2023a software (Mathworks Inc., Natick, MA, USA) was used for noise removal and parameter calculation. The acceleration value in the measured three-axis acceleration values, Medial (+)-lateral (-) (ML), superior (+)-inferior (-) (SI), and Anterior (+)-Posterior (-) (AP), were used in the analysis. To remove noise from all acceleration values, a fourth-order zero-lag Butterworth bandpass filter was applied, and the cut-off was set to 0.5 to 5 Hz to preserve the typical gait signal that exists between 0.5 and 5 Hz.22,23 The identification of gait events and phases (stance and swing) was based on heel-strike and toe-off timings obtained from the 3D motion capture data. Specifically, the stance phase was defined from heel-strike to toe-off, and the swing phase from toe-off to the subsequent heel-strike. These events were automatically detected using the motion analysis software, which identifies gait cycle events based on the position and velocity of heel and toe markers, in line with established gait analysis protocols. Nine spatial-temporal parameters were derived from the data measured by the 3D motion analysis system: stride, stance and swing time, and gait events such as stance and swing phase ratio, mean velocity, cadence, and stride and step length. Variables derived from the IMU are basic parameters and acceleration movement parameters. The basic parameters are mean acceleration, which is the average of the measured acceleration values, maximum acceleration, which is the maximum value, and range acceleration, which is the range between the maximum and minimum values. The acceleration movement parameter is the average value representing the change in center of mass (CoM) calculated from the acceleration movement, which is the center of mass of the acceleration value itself, and the distance value obtained by integrating the acceleration value twice. The three basic parameters and the two acceleration movement parameters, total five variables are divided into ML, SI, and AP of the three axes, and there are 15 types, and since they are values measured from each thigh and shank, there are a total of 30 variables. In addition, in order to calculate the ellipse area for the CoM in the three-axis direction, the mean and covariance matrix of the data were calculated to calculate the eigenvector, and from this, the parameter settings and the ellipsoid based on the average value were expressed in the figure. This was calculated to confirm the relative tendency difference between the affected and un-affected side CoM in the thigh and shank.
The gait data of 40 hemiplegic patients were divided into two groups, un-affected and affected, based on the clinical diagnosis confirmed in advance. Paired T-test was performed for statistical comparison of all variables between un-affected and affected sides. In addition, Cohen's d was calculated to confirm the effect size.
Results
As a result of gait measurement, the mean velocity of all subjects was 0.78 ± 0.35 m/s, and the cadence was 98.78 ± 17.02 setps/min (mean ± standard deviation). In the results of spatial-temporal parameters, there were statistical differences in stance time, swing time, swing phase, and stance phase (Table 2).
Characteristics of subjects.
Characteristics of subjects.
Spatial-temporal parameter.
Paired T-test: **p < 0.01; Effect size: Cohen's d; value: Mean (Standard deviation)
Among the acceleration basic parameters of the thigh and shank, the mean acceleration and maximum acceleration, the variable in which statistical differences were confirmed between the un-affected and affected sides was the mean acceleration of AP-thigh. The acceleration in the posterior direction was higher on the affected side than on the un-affected side (Table 3).
Acceleration basic parameter.
Acc.: Acceleration; AP, anterior-posterior direction; ML, medial-lateral direction; SI, superior-inferior direction; Paired T-test: *p < 0.05; Effect size: Cohen's d; value: Mean (Standard deviation)
Acceleration movement parameter.
AP, anterior-posterior direction; ML, medial-lateral direction; SI, superior-inferior direction; Paired T-test: *p < 0.05; Effect size: Cohen's d; value: Mean (Standard deviation), CoM: Center of Mass
In Acceleration movement, which is the center of mass of the acceleration value itself, which is the acceleration center parameter of the thigh and shank, there was no statistical difference between the un-affected and affected sides. In CoM, which is calculated from the distance value obtained by integrating the acceleration value twice, there was a statistical difference in AP-thigh between the affected and un-affected sides (Table 4).
Figure 1 is the result of fitting the ellipse area with the CoM value calculated from the distance value, which is the integral value of the velocity calculated from the acceleration. Based on the transverse plane viewed from above the head, the area of the affected side of the thigh was larger than that of the un-affected side, and the CoM increased in the posterior and lateral directions. The difference in the area of the shank was not large, and the CoM of the affected side showed a tendency to slightly increase in the medial direction.

Ellipse area of affected and un-affected side center of mass.
In this study, in addition to the spatial-temporal parameters, which are commonly used variables in comparing the un-affected and affected sides of stroke hemiplegia gait, we conducted a basic study to derive parameters that can be evaluated at the local location of the patient's lower limb from acceleration data measured using a simplified sensor, IMU.
First, in the spatial-temporal parameters, which are common gait evaluation variables, the variables that could confirm significant differences between the un-affected and affected sides were stance and swing time, and stance and swing phase ratio. We confirmed the same trend as the shortening of the stance phase interval and time, which are characteristics of the affected side of stroke hemiplegia patients in general. This is the same as the result of gait asymmetry that occurs when the time to support the weight with the paralyzed leg becomes shorter due to muscle weakness and balance instability, as suggested in various previous studies.24,25 Since the subjects participating in this experiment showed the characteristics of common hemiplegic patients, it can be interpreted as the result of appropriate subject selection.
Second, the result of calculating the mean, max, and range of the acceleration raw data, which is the Acceleration basic parameter, a significant difference was confirmed between the affected side and the un-affected side only in the mean value of the Thigh in the AP direction. The affected side showed a larger acceleration value on average in the posterior direction. The section where the Thigh increases to a negative value in the opposite direction of the gait progression is the section from mid stance to terminal stance. It is thought that the characteristics of the affected side tend to be that the phase ratio and time of the mid stance tend to be shorter, and that it tends to end the mid stance with a larger acceleration on average. 26 However, since there was no statistical difference in the variable indicating the range of the maximum and minimum values of acceleration, additional study is necessary.
Third, the statistical comparison result of the acceleration movement parameter, the center of mass value defined as the average value of the distance calculated from the acceleration also showed a statistically significant difference between the affected side and the un-affected side only in the AP direction value of the Thigh, similar to the mean acceleration. The center of mass refers to the center of mass calculated from the distance value. When the center point for the three-axis value was considered, it increased in the posterior direction. This is related to the weakened muscle strength and difficulty maintaining balance due to paralysis. 27 Due to the influence of the weakened extensor muscle, the propulsive force to move the body forward on the affected side is insufficient, causing the center of mass to shift backward. 28 In addition, foot drop, etc. occur due to the decline in dorsiflexion function, which can also cause the center of mass to shift posteriorly. This can be explained by the results of a previous study. 29
Fourth, Figure 1 confirms the characteristics of the thigh and shank by the size of the ellipse area. Since the absolute numerical results of each subject are different depending on the characteristics such as severity and body length, the tendency of area movement was confirmed for relative comparison, and the thigh-shank coordination. 30 This is thought to be because the center of mass of the body is located in the pelvis during normal walking, so the change in the center of mass is more evident in the thigh, which is closer to it. 31 In addition, it is thought that the thigh, which induces the hip and knee angle, has a greater influence on movement than the shank, which is bound to the ankle.
The limitations of this study are as follows. First, the measured acceleration data was not verified based on the section set as a gait event. This is because the purpose of the result derivation method in this study was to compare the overall un-affected and affected sides using acceleration in the entire 5-m walking section rather than gait events. In a follow-up study, we plan to verify which walking phase the results actually occur in based on the data set as gait events. Second, there is a possibility of error in the values used by integrating the acceleration data of the IMU sensor. This is a shortcoming of most IMU sensors. However, it can be compensated to some extent by applying various filters. In future studies, it is thought that the shortcomings in data post-processing can be compensated for by applying more sophisticated filters. Additionally, given the possibility of Type I error, it is important to review the results by applying the Bonferroni correction. However, given the relatively small number of participants in this study, we believe that the results of the current statistical analysis should be interpreted at an exploratory level rather than generalizable.
Conclusion
In this study, rather than directly using acceleration data to identify the characteristics of the affected side of hemiplegic patients, we confirmed that gait evaluation is possible using the CoM value calculated from the distance value. In particular, the increase in the average CoM value in the posterior direction of the thigh, which reflects the characteristics of the affected side, is a gait characteristic of hemiplegic stroke patients, so it can be used as a criterion for gait rehabilitation treatment that aims to eliminate the difference between the un-affected side and the affected side. In future studies, we aim to develop more generalizable diagnostic criteria using the results targeting a larger number of stroke patients.
Footnotes
Acknowledgements
This research was supported by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number : RS-2020-KH088006)
Ethical approval
The studies involving human participants were reviewed and approved by the Ethics Committee of Wonkwang University Gwangju Hospital (WKUGH) (IRB No.: WKIRB 2022-07, June 29, 2022). The patients/participants provided their written informed consent to participate in this study.
Author contributions
SEO JS: writing—original draft, formal analysis, investigation, Lee S: writing—review and editing, and visualization, conceptualization and resources, supervision, and funding acquisition. All authors have read and agreed to the published version of the manuscript.
Funding
The author disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Korea Health Industry Development Institute, (grant number RS-2020-KH088006).
Conflicts of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Data availability statement
All relevant data is contained within the article: The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
