Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Enhancing Speech Emotion Recognition Using Dual Feature Extraction Encoders
Next Article in Special Issue
Swimmers’ Effective Actions during the Backstroke Start Technique
Previous Article in Journal
Gaussian-Filtered High-Frequency-Feature Trained Optimized BiLSTM Network for Spoofed-Speech Classification
Previous Article in Special Issue
The Influence of Anthropometric Variables and Filtering Frequency on Center of Pressure Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Control Method for Gait Detection and Classification Devices with Inertial Measurement Unit

1
Division of Mechanical Engineering, (National) Korea Maritime and Ocean University, Busan 49112, Republic of Korea
2
Division of Biomedical Engineering, Konkuk University, Chungju 27478, Republic of Korea
3
BK21 Plus Research Institute of Biomedical Engineering, Konkuk University, Seoul 05029, Republic of Korea
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(14), 6638; https://doi.org/10.3390/s23146638
Submission received: 20 January 2023 / Revised: 20 July 2023 / Accepted: 21 July 2023 / Published: 24 July 2023

Abstract

:
Cueing and feedback training can be effective in maintaining or improving gait in individuals with Parkinson’s disease. We previously designed a rehabilitation assist device that can detect and classify a user’s gait at only the swing phase of the gait cycle, for the ease of data processing. In this study, we analyzed the impact of various factors in a gait detection algorithm on the gait detection and classification rate (GDCR). We collected acceleration and angular velocity data from 25 participants (1 male and 24 females with an average age of 62 ± 6 years) using our device and analyzed the data using statistical methods. Based on these results, we developed an adaptive GDCR control algorithm using several equations and functions. We tested the algorithm under various virtual exercise scenarios using two control methods, based on acceleration and angular velocity, and found that the acceleration threshold was more effective in controlling the GDCR (average Spearman correlation −0.9996, p < 0.001) than the gyroscopic threshold. Our adaptive control algorithm was more effective in maintaining the target GDCR than the other algorithms (p < 0.001) with an average error of 0.10, while other tested methods showed average errors of 0.16 and 0.28. This algorithm has good scalability and can be adapted for future gait detection and classification applications.

1. Introduction

1.1. Research Background and Related Literatures

Parkinson’s disease (PD) is a neurodegenerative disorder that is caused by a deficiency in dopamine. The main symptoms of PD can be classified as either non-motor or motor. Non-motor symptoms include things like depression and insomnia, while motor symptoms include gait disorders and freezing of gait. Music therapy [1], cueing training [2], and feedback training [3] are all methods that have been used to rehabilitate gait disorders in individuals with PD. These training methods can be effective, particularly when used in conjunction with medication or as part of a standalone rehabilitation program [4,5,6,7].
One form of rehabilitation that has been shown to be particularly effective for PD-related gait disorders is augmented feedback training. Many researchers have explored the use of augmented feedback training to improve gait and alleviate gait disorders [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. However, in order to effectively use these training methods, it is necessary to have a gait analysis system that can detect and evaluate a user’s gait. Gait detection is important for determining whether a person is walking or not, while gait classification is necessary for evaluating the quality of the detected gait. Augmented feedback training methods rely on this information to provide appropriate feedback, such as auditory, visual, or tactile signals.
Traditionally, the gold standard for quantitative gait analysis has been optoelectronic stereophotogrammetry. However, this method requires a large laboratory and a significant budget to support its operation and maintenance [25]. An alternative to optoelectronic stereophotogrammetry is the use of inertial measurement units (IMUs). IMUs are portable and easy to use, making them a popular choice for gait analysis in both daily life and rehabilitation settings [26,27]. The use of IMUs for gait analysis is now widespread.
The IMU gait analysis system was previously characterized by multi-point IMU (MIMU) systems [28], which were also validated through various research studies [29,30,31,32,33,34]. Although MIMU systems have demonstrated their superiority, issues with sensor adjustment, calibration, and consistency have limited their adoption in practice [35]. In recent years, single-point IMU (SIMU) devices have gained competitive advantages as medical devices, particularly in the assessment of the severity of Parkinson’s disease, fall risk, and step count [36,37,38,39,40,41].
According to a review of wearable gait detection devices [42], 23 journal papers have used SIMUs on the shank and 25 journal papers have used SIMUs on the foot in an attempt to make SIMUs a popular wearable system for gait analysis. These papers have primarily focused on identifying a user’s gait pattern using artificial intelligence (AI) and machine learning (ML) methods, such as artificial neural networks, extreme learning machines, convolutional neural networks, and long short-term memory [43,44,45,46,47], or rule-based algorithms to detect gait phases such as heel strike, toe-off, and swing phase [48,49,50,51,52].
However, the authors of this paper believed that a simpler gait analysis method was needed to use SIMUs in daily life. AI and ML require a large amount of data to set personal preferences before use, and rule-based gait phase detection algorithms are too comprehensive for our purposes. Therefore, we developed a device that detects and evaluates gait based on acceleration and angular velocity around the ankle joint at the swing phase of gait cycle [53]. This device uses a gait detection and classification (GDC) method to detect and classify a user’s gait based on various factors and the GDC frequency of the device user, known as the gait detection and classification rate (GDCR). However, we found that GDCR can only be roughly controlled and it is necessary to conduct further research to develop an adaptive GDCR control method.
The development of GDCR control methods is also an invaluable part of improving the previous GDC device and method. During feedback training, the results are divided into two categories: success (positive feedback) and failure (negative feedback). Usually, the frequency of provided positive and negative feedback is closely related to the user’s motivation systems in common [54,55,56,57]. Therefore, the adaptive GDCR control method is crucial for preparing for expected and unexpected events that can affect walking speed, such as medication treatment, accidental injury events, or the gradual degradation of motor control ability during rehabilitation, to maintain, enhance, or prevent the user’s motivation.

1.2. Summary of the Previous Research

This research was based on previous work [53] and involved the development of a device for detecting and classifying gait based on angular velocity around the ankle joint during the swing phase of a gait cycle as shown in Figure 1. The device used modified (out of anatomical axis, i.e., directly measured by the device) the medial-lateral axis angular velocity and the magnitude of acceleration on the Sagittal plane around the ankle joint. In what remains of this article, we will loosely use the terms “angular velocity around the ankle” and “ankle acceleration” to refer to the values measured by the IMU. As shown in Figure 1, the factors changed significantly during the swing phase of the gait cycle, so we focused on this phase for gait detection and classification.
The GDC method used a gyroscopic threshold (GT) to detect the user’s gait. If the z axis angular velocity of the user was recorded above the GT, the method would detect the user’s gait and set a “gait detection range.” The device also evaluated the user’s gait based on acceleration around the ankle joint using an Acceleration Threshold (AT). If the average of the magnitude of acceleration during the gait detection range was recorded above the AT, the user’s gait would be classified as a “good gait.” The difference between the AT and GT is that the AT is based on the “average value” of data during the gait detection range, while the GT is based on the “value” of the data.
In the previous research, the AT and GT were used to control the difficulty or ease of gait detection and classification. The GDC method and related devices were used to detect and classify the strength of a user’s gait using the AT and GT as factors. Figure 1 and Figure 2 provide more information on these definitions and how the GDC method works. The GDCR was found to be affected by factors such as walking speed and the presence of gait disorders. The GDCR was also a useful measure for observing and assessing gait conditions.
Twenty-five subjects participated in the assessments of the previous experiment and we collected the acceleration and angular velocity around the ankle joint during one gait cycle. Then, GDC algorithm was applied to the collected data and we observed relations between AT, GT, walking speed, and gait symptom to GDCR. The GDCR was inversely proportional to the AT, GT, and gait disorders, but directly proportional to the walking speed, as shown in Figure 3.

2. Materials and Methods

2.1. Experiments for Data Collection

Experimental data were collected to develop a GDCR control algorithm and analyze gait. A device that measures acceleration and angular velocity around the ankle joint was used to collect the data. The device was attached to a support band using a metal clip, and subjects were able to wear it themselves, although the experimenter adjusted the devices for better consistency of data. Before the experiment, the devices were calibrated by having subjects stand up straight for 30 s to measure acceleration in their default position. The default position coordinate data were used to calibrate the acceleration. The collected data were processed after the experiment and used to calculate the GDCR at different ATs and GTs. The system provided feedback at every GDC fault.
The subjects included 24 females and 1 male, with an average age of 62 years. None of the subjects had diseases that caused gait disorders. Each subject was asked to walk for one minute at 90 steps per minute cadence (45 times with each leg), with the experimenter providing auditory metronome sounds to indicate the exact timing for walking. Subjects walked three times per each case at treadmill speeds of 1.0, 1.5, 2.0, 2.5, and 3.0 km/h. Subjects took a one minute rest and then walked again, following the sequence walk–rest–walk–rest–walk. If there were more data than 90 steps, only the first 90 steps were used. The experimenter controlled the treadmill speed in ascending order (1 km/h–3 km/h).
The devices used were an MPU-6050 IMU sensor with an ATmega328P chip as the main control unit, and data were collected at a frequency of 30 Hz. 30 Hz was determined to prevent data shifting (phenomenon that some data was recorded at another data line) and data losing (phenomenon that losing IMU data) from low data processing speed. The accelerometers had a sensitivity of 16,384 sensor values per gravity acceleration (9.81 m per second squared (16,384 LSB/(9.81 m/s2)), and the gyroscope had a sensitivity of 131 sensor values per degree per second [131 LSB/(deg/s)]. Data generated by the subjects was transmitted via a 2.4 GHz Xbee wireless communication system. Figure 4 shows an image of the equipment and experimental setup.

2.2. Impact of the AT and GT on the GDCR

IBM’s SPSS software (SPSS 25) was used to perform statistical analysis to examine the impact of the walking speed, the left and right legs, and the AT and GT on the GDCR. Various statistical methods were used to determine the correlation or difference in the average between each element. Non-parametric tests, the independent sample Mann–Whitney U and Kruskal–Wallis tests, were used to assess the effect of walking speed and the left and right legs on the GDCR. Spearman’s correlation analysis was used to examine the correlation between the AT, GT, and GDCR. The relation between GT and GDCR came by comparing GT and GDCR according to GT at each walking speed. For example, the GT (15,000~25,000) and GDCR (0~1) according to the GT at each walking speed (1.0 km/h~3.0 km/h) was compared. The relation between AT and GDCR came by comparing AT and GDCR according to AT at each walking speed. For example, AT (1000~15,000) and GDCR (0~1) according to the AT at each walking speed (1.0 km/h~3.0 km/h) was compared. Pearson’s correlation analysis was not suitable for this analysis because the experimental dataset did not meet normality. Table 1 provides more details on the statistical analysis.
Twenty-five subjects participated in this research; however, some data on left leg and walking speed 1.0 km/h (n = 3) and 2.0 km/h (n = 1) were missed. We could compare the impact of walking speed after deleting missing data. However, we could not compute the statistical analysis for the relation between GDCR-GT and GDCR-AT since, for consistency, they can only be compared when we have data for every walking speed.
Number of subjects and samples are different, the number of samples were calculated from the number of subjects, factors, and flaws of data, and we used the number of samples as statistical analysis data. For example, the number of cases of GDCR-Leg relation was 25 on the left leg and 25 on the right leg, but data from 4 participants had flaws on left leg and walking speed 1.0 km/h (n = 3) and 2.0 km/h (n = 1). Therefore, the number of samples was 121 (cases * walking speeds—flaws, 25 × 5 − 4 = 121) on left and 125 (25 × 5 − 0 = 125) on the right leg. The dependent variable was GDCR, and the independent variable was the leg.
Also, the number of samples of the relation GDCR—walking speed was calculated as follows; the number of subjects × legs (2, left and right) − flaws = the number of samples. The number of samples at 1.0 km/h of walking speed was 47 (25 × 2 − 3 = 47), at 2.0 km/h of walking speed was 49 (25 × 2 − 1 = 49), 50 for other walking speed (25 × 2 − 0 = 50). The dependent variable was GDCR, and independent variable was walking speed.
The number of samples of the relation GDCR-AT was calculated as follows; the number of subjects without flaws of data × the number of AT points of the AT-GDCR curve (1901 points, from 1000 to 20,000, every 10 AT) ((25 − 4) × 1901 = 39,921). The dependent variable was GDCR, and the independent variable was AT. In the case of the relation GT-GDCR, the number of samples is calculated by the number of subjects without flaws of data − the number of GT points of the GT-GDCR curve (21 points, from 15,000 to 25,000, every 500 GT) ((25 − 4) × 21 = 441). Dependent variable was GDCR, independent variable was GT.

2.3. Algorithm

The device designed in previous research [53] and improved in this study provided GDC results for the users through various feedback types: visual, auditory, and vibration (somatosensory). The frequency of feedback given when a gait was classified as good is important for gait rehabilitation. The designed rehabilitation device provides feedback whenever a gait fails (when acceleration around the ankle joint cannot be higher than AT). Therefore, the GDCR and the frequency of the gait failure have an inverse relationship. The frequency may vary depending on the rehabilitation method or goal, but if it is too high (>30%) due to high difficulty, users may feel tired and overwhelmed by the device. On the other hand, users may feel that the rehabilitation is insufficient if the frequency is too low (<5%). Therefore, an important feature of the algorithm is its ability to maintain the GDCR, regardless of previous exercise history. Thus, an algorithm was developed that can set the GDCR to be as close as possible to the target GDCR, based on the average GDCR data from 25 subjects. The GDCR target of this algorithm was 0.9 (90% of classified gaits among detected gaits).
The relationship between the AT and the average GDCR of the experimental subjects can be represented as a sigmoid-shaped fitted curve. The function FG was calculated by multiplying the average GDCR from the experimental dataset by 1.38, as shown in Figure 5. The reason of multiplying 1.38 is that to match the range of two different functions: FG and other sigmoid shape functions. The value of 1.38 can be changed into others based on the value of the average GDCR at AT = 1000 (this can be seen in Figure 5). Therefore, the maximum value of FG was 1. In Equation (1), the maximum value of AT was 15,000, and this value was selected based on Figure 5, where it is seen that for high values of AT the GDCR approaches zero. Average GDCR at over 15,000 AT does not significantly affect GDCR. Also, AT at the equation started from 1000 for the same reason. Various sigmoid-shaped functions, such as logistic, error, algebraic, and arctangent functions, were then used to evaluate the similarity to FG. The average error and absolute value of the average differences of these functions with FG were used to determine the similarity. Equations (2)–(5) show the functions used. The average error was calculated using Equation (6). ‘The number of errors’ of Equation (6) was 1400, since the unit of AT was 10 (from 1000 to 15,000, every 10). Figure 6 shows the curves of the various functions and the average errors.
F G A T = N o r m a l i z e d   f i t t e d   c u r v e   ( 1000 A T 15,000 , 0 F G 1 )
f 1 x 1 = L o g i s t i c x 1 = 1 1 + e x 1 ( 5.6 < x 1 < 7 )
f 2 x 2 = E r r o r x 2 = 0.5 + 1 π 0 x 2 e t 2 d t ( 2.5 < x 2 < 2 )
f 3 x 3 = A l g e b r a i c x 3 = 0.5 x 3 2 1 + x 3 2 ( 16.5 < x 3 < 13 )
f 4 x 4 = A r c t a n g e n t x 4 = 0.5 1 π arctan π x 4 2 ( 7 < x 4 < 5.6 )
A v e r a g e   e r r o r n = ( F G ( x ) f n x n ) N u m b e r   o f   e r r o r s ( 1000 x 15,000 )
The logistic function with the least average error of 0.0131 was the function most similar to the FG among the various sigmoid-shaped functions. Therefore, the logistic function’s domain range was converted into the AT domain from 1000 to 15,000 using a trendline, as shown in Equation (7). Subsequently, the function G had converted the x1 domain, xAT based on logistic functions such as Equation (8).
x A T ( A T ) = 0.009 · A T 1000 10 + 7.009 ( 1000 A T 15,000 , 5.6 x A T 7 )
L o g i s t i c A T = G ( A T ) = 1 1 + e 0.009 · A T 1000 10 7.009 ( 1000 A T 15,000 )
To determine the features of an individual user, we developed an algorithm that personalizes the GDCR control. The algorithm is based on the average GDCR data from the 25 subjects and is designed to maintain the GDCR regardless of the user’s previous exercise history. The initial personalized coefficient, Cp1, was calculated as 1.11 times the GDCR at a walking speed of 2 km/h and an AT of 6000. This value was obtained by dividing the average GDCR of experimental data at AT = 1000 by the average GDCR of experimental data at AT = 6000.
For values of n greater than two, Cpn can be calculated using Equation (10). The personalized GDCR function, Gp, can be determined by AT to predict the expected GDCR of each user, as shown in Equation (11). Equation (12) describes the normalized GDCR, Gnor, which is a function derived from FG. Equation (13) describes the relationship between the AT function and the G function, and Equation (14) is the AT function that calculates the appropriate AT value to achieve the target GDCR, Gtarget.
C p 1 = a v e r a g e   G D C R   a t   A T = 1000 a v e r a g e   G D C R   a t   A T = 6000   actual   GDCR at   AT = 6000   and   2   km / h = 1.11   GDCR
Cp1 is the initial personalized coefficient. The average GDCR was from all of experimental subjects, and the actual GDCR would be from new experimental data that will be collected from experiments.
C p n = a v e r a g e   G D C R   a t   A T = 1000 · a c t u a l   G D C R E x p e c t e d   G D C R = 0.720 · a c t u a l   G D C R G n o r A T ( n 2 )
Cp2 is the personalized coefficient to control the algorithm. ‘Expected GDCR’ means expected GDCR with the designed logistic function at current AT and GT. Therefore, Gnor can be utilized for the purpose.
G p A T = C p n 1 + e 0.009 · A T 1000 10 7.009 1000 A T 15,000
G n o r ( A T ) = a v e r a g e   G D C R   a t   A T = 1000 1 + e 0.009 · A T 1000 10 7.009 ( 1000 A T 15,000 )
The form of Gp and Gnor was that of the selected sigmoid shape, the logistic function. The x domain of the logistic function changed into ‘0.009 × (AT − 1000)/10 − 7.009’ to utilize the range of AT as a x domain of Gp and Gnor equations. Also, each equation was multiplied by Cpn or ‘average GDCR at AT = 1000’ to change the range of value of the equation.
A T ( G t ) = G t 1 A T 1000 A T 15,000 , 0 G t 1
A T G t = 1111.11 · log C p G t 1 + 8787.77 ( 0 G t 1 )
AT(Gt) is the inverse function of Gt(AT) because Gt(AT) is defined as the target GDCR according to the current value of AT and its inverted function, AT(Gt) can mean ‘required AT to reach to Gt’. Gt can be determined from Gp. This will be described when we introduce the algorithm process.
To design the GDCR adaptive control algorithm using Equation (9) from Equation (14), we followed these steps: (1) Measure Cp1 at a walking speed of 2 km/h and AT = 6000; (2) Set Gtarget to 0.9 if Cp1 is greater than 0.9, or set Gtarget to 0.01 less than Cp1 if Cp1 is lower than 0.9; (3) Set a new AT at AT(Gtarget) (initial setting end); (4) Exercise with the new AT and collect the new GDCR; (5) Calculate the new Cpn; (6) Set Gtarget to 0.9 if Cpn is greater than 0.9, or set Gtarget to 0.01 less than Cpn if Cpn is lower than 0.9; and (7) Set a new AT with AT(Gtarget). Repeat these steps starting at step 4. Figure 7 shows the design process.
This algorithm can lead to confusion surrounding the following question: ‘Will the algorithm classify the first steps as bad and then change the AT value to classify the other steps as ‘good’ since it is an adaptive algorithm?’ However, this adaptive control algorithm is almost a daily-based algorithm. Rather than changing with every step, the idea is that over the workout of a day, for example, about 2 h of walking per day, the results are summarized to obtain an average GDCR and Cpn (as shown in Equation (10)) for the day, which is then fed into the algorithm to change for the next day’s or next time’s workout. In this sense, the goal of the algorithm is to keep the frequency of feedback at a similar level while increasing the difficulty so that it can be used for feedback training.
To evaluate the controllability of the adaptive control algorithm for GDCR, we designed various virtual walking-speed pattern scenarios based on experimental dataset (Walking speed-GDCR and its interpolations). These scenarios consisted of nine combinations of five walking-speed patterns: linear increase, linear decrease, random with upper and lower bounds, sudden increase, and sudden decrease, as shown in Figure 8.
The list below is description for each scenario:
  • The first scenario presents a linear increase in walking speed.
  • The second scenario presents a linear decrease in walking speed.
  • The third scenario presents walking speeds with punctuations within a range.
  • The fourth scenario presents an increasing trend in walking speed followed by a sudden decrease in the middle.
  • The fifth scenario presents a decreasing trend in walking speed followed by a sudden increase in the middle.
  • The sixth scenario presents a V-shaped trend in walking speed.
  • The seventh scenario presents a V-shaped trend in walking speed with a range of random punctuations in the middle.
  • The eighth scenario presents a sudden increase in ranged random punctuations after a linear decrease in walking speed.
  • The ninth scenario presents ranged random punctuations after a linear increase in walking speed.
Also, the list below is description for algorithm processing:
  • For an initial setting, measure the first Cp1 at 2 km/h of walking speed and 6000 AT.
  • If Cp1 is recorded as smaller than 0.9, set Gt as 0.01 below Cp1 (Cp1—0.01). Otherwise, if Cp1 was over 0.9, set Gt as 0.9.
  • Set new AT for later rehabilitation as AT (Gt);
  • Do exercise with set AT (Gt) and measure new GDCR as a result of exercise.
  • Calculate Cpn from the measured GDCR.
  • If Cpn was smaller than 0.9, set Gt as 0.01 smaller than Cpn (Cpn − 0.01). Otherwise, if Cpn was over 0.9, set Gt as 0.9.
  • Set new AT as AT (Gt)
  • Repeat process from 4.
The initial setting is from 1 to 3, and the adaptive control algorithm for further rehabilitation is from 4 to 8.
In addition, we prepared control methods to compare the accuracy of the adaptive control method with others;
  • Adaptative, designed method that was previously described;
  • Fixed, which simply fixes the AT at 6000;
  • Simple, which increases or decreases the AT by 100 from the initial value of ‘AT = 6000’ when the GDCR exceeds 0.95 or falls below 0.9, respectively.
These control methods were then applied to the exercise scenarios based on experimental data from subjects. This resulted in a total of 27 simulation results per subject, or 675 in total. The GDCR trends from the various scenarios and algorithms were statistically analyzed. The performance of the control methods was compared using the average error from the target GDCR to 0.9 in all scenarios. The statistical analysis included independent sample Kruskal–Wallis and Mann–Whitney U tests.
Simulation with each scenario was ran based on the number of algorithm cycle repetition. For example, if walking speeds change from one to three, as shown in Figure 9, the algorithm measure GDCR is based on the AT(Gtarget) at the walking speed = 1.0 km/h. Then, the algorithm return appropriate new AT value ‘AT(Gtarget)’ for the next algorithm repetition. The simulation applied new AT(Gtarget) to the new GDCR curve at a walking speed = 2.0 km/h and returned a new AT value. These processes repeated until simulation ended. Consequently, the average error to GDCR = 0.9 could be calculated.

2.4. Our Goal

Before presenting the results section, we would like to remind readers of what we obtained. Firstly, we calculated the average GDCR for both the left and right legs to determine whether there were any significant differences between the two. Secondly, we computed the average GDCR for different ascending walking speeds to observe the effect of walking speed on GDCR. Thirdly, we examined the impact of GT and AT on GDCR to determine which factor had a dominant effect on controlling GDCR. Finally, we investigated how the adaptive algorithm could sustain a constant GDCR as much as possible under various scenarios.

3. Results

3.1. Statistical Analysis Results of the Effects of the AT and GT on the GDCR

The statistical analysis in this section was conducted with IBM’s SPSS, and the results are based on how the program works. The average GDCR on the left and right legs exhibited a non-normal distribution (p < 0.001, Kolmogorov–Smirnov test), but it differed significantly (p < 0.001, n = 246, Mann–Whitney U test). The number 246 arises from considering each leg separately, summing 121 from the left and 125 from the right leg, as shown in Table 1. Between the left and right legs, at 0.49 and 0.45, respectively. Although the samples used to examine the correlation between walking speed and GDCR, throughout entire AT range (from 1000 to 15,000), were not normally distributed (p < 0.001, Kolmogorov–Smirnov test), the average GDCRs significantly increased (p < 0.001, n = 246, Kruskal–Wallis test) with walking speeds of 0.27, 0.45, 0.49, 0.52, and 0.59, respectively; 1.0 km/h on the left leg had 47 samples (number of cases × 2 − flaws, 25 × 2 − 3 = 47); 2.0 km/h on the left leg had 49 samples (25 × 2 − 1 = 49); and others had 50 samples (25 × 2 − 0 = 50). Therefore, the sum of samples was 246. These results are shown in Figure 10.
The dependent variable was GDCR for both, and independent variable were the GT and AT to calculate Spearman correlation coefficients. As a result, the GT was not a significant factor in controlling the GDCR. The Spearman correlation coefficients were −0.353, −0.210, −0.610, −0.109, and −0.056, respectively, and the walking speed were in ascending order of. GT (15,000~25,000) and GDCR according to GT (depending on walking speed) were compared. (p < 0.001, n = 441 for each walking speed as shown in Table 1, n = 2205 in sum). On the other hand, the AT was found to be a significant factor in controlling and reducing the GDCR. The Spearman correlation coefficients were −0.998 at 1 km/h and −1.000 for all other cases. AT (1000~15,000) and GDCR according to AT (depending on walking speed) were compared (p < 0.001 n = 39,921 for each walking speed as shown in Table 1, n = 199,605 in sum). The AT was the most significant factor in controlling the GDCR, and the adaptive control algorithm was designed based on the AT. Figure 11 and Figure 12 show the GDCR according to the GT and AT, respectively, and the corresponding Spearman correlation coefficients.

3.2. Algorithm Evaluation

In this section, the errors will not be computed from the entire range of AT but from a narrower range, dependent on the control algorithm’s actual performance. The statistical analysis showed that the adaptive control method was the best among three control algorithms (n = 675, p < 0.001, Kruskal–Wallis test). Error of the result means an absolute value of ‘actual GDCR—0.9.’ The adaptive control method differed significantly from the others (n = 675, p < 0.001, Mann–Whitney U test). The number of samples, n, from one control method and one scenario was 25. Therefore, each control method had 225 (the number of scenarios * the number of participants, 9 × 25 = 225) samples and 675 in sum. The average errors for the adaptive, fixed, and simple control methods were 0.10, 0.16, and 0.28, respectively. Figure 13 shows the results of this analysis.
Additionally, the average errors from the experimental dataset were observed for each control method and scenario. The adaptive control method had the best performance in scenario 3, a random scenario, with an average error of 0.06. The highest error was observed in scenario 4, with an average error of 0.14. The fixed control method also had the best results in scenario 3, with the lowest average error of 0.10. The worst performance of the fixed control method was observed in scenario 4, with an average error of 0.18. The simple control method had the lowest average error in scenario 7 (0.13) and the highest average error in scenario 1 (0.66). Figure 14 shows the average errors for all the control methods.
The goal of the adaptive algorithm is to maintain the GDCR close to 0.9; however, why were the average GDCRs for left and right legs 0.49 and 0.45, respectively, in Figure 10? These values came from the full range of AT (from 1000 to 15,000). Therefore, as shown in Figure 12a, the GDCR curve was close to 1.0 at a lower AT and close to 0 at a higher AT. On the contrary, the errors came from algorithm usage. Therefore, errors were calculated from the algorithm control results; consequently, GDCR was usually higher than 0.5 since the control algorithm was designed to sustain GDCR close to 0.9. In short, the GDCRs equal to 0.49 and 0.45 from each leg were from the entire AT range (from 1000 to 15,000), while errors from the GDCR control algorithm came from a narrower AT range.

4. Discussion and Conclusions

This study investigated the correlations between the GDCR and various elements used in GDC devices and developed an adaptive GDCR control algorithm to control the GDCR and feedback frequency during exercise using a set AT. The adaptive control algorithm improves the user experience by reducing mental fatigue associated with a higher feedback generation frequency and addresses user uncertainty about the feedback and the functionality of the device when the GDCR is high. However, some readers may still have some confusion from our research. To add more description for better understanding and to avoid confusion, we ask: why the ‘Steady walking speed scenario was not included into virtual scenarios?’ We did not need to include a scenario with constant walking speed. The main purpose of the designed algorithm is to observe how the algorithm can compensate fluctuating GDCR from different walking conditions. At the constant walking speed scenario (1.5 km/h), we expected the AT and GDCR to be invariant with a constant walking speed, and this is actually right. Moreover, we assumed that daily walking speeds would be subtly different in a real-world rehabilitation scenario, so the constant walking speed scenario was altered by including a randomized scenario, i.e., scenario 3.
Previous studies have collected and classified gait factors for diagnosing and rehabilitating patients with Parkinson’s disease, stroke, and other conditions [58,59,60,61,62,63]. Further experiments are needed to determine the effectiveness of the device developed in this study in improving gait symptoms. Other researchers have used IMU systems to quantitatively assess gait [64,65,66,67]. AI has also been applied to assess, evaluate, and extract gait features using IMU systems [43,44,45,46,47]. To the best of our knowledge, there has not been any similar research using an extremely simplified gait detection and assessment system with an associated control algorithm. The device developed in this study uses only one element, AT, to control the GDCR, and does not employ any machine learning methods or AI systems for gait analysis. Therefore, the GDC device and its adaptive control algorithm can roughly but easily detect changes in the user’s condition.
The developed algorithm has good scalability because it uses acceleration around the ankle joint, which can be measured by various devices. For example, the gait rehabilitation algorithm can be used with inertial motion capture systems commonly used in virtual reality production, allowing patients to rehabilitate in a virtual reality environment using a treadmill. The algorithm also has a simple design, making it reproducible and repairable. The device is lightweight, small, and wearable, and can be attached to an ankle brace or another support for ease of use. However, the device already adopted some components to generate a feedback signal needed to be advances. Therefore, additional advanced components such as better and smaller speakers for auditory feedback or lasers for visual feedback can be added to enhance the user experience and portability.
The results of this study confirmed the validity of the existing GDC algorithm after the statistical analysis of the research data. The device and algorithm in this study were designed based on this analysis. Further research is needed to confirm the medical effectiveness of the device and algorithm and to determine their potential for use in clinical settings. In conclusion, the GDC device and adaptive control algorithm developed in this study can provide a rough but easy detection of changes in the user’s condition and have the potential to improve gait symptoms and enhance the user experience.

Author Contributions

Conceptualization—data curation, validation, and writing (original draft, review, and editing), H.K., J.-W.K. and J.K.; formal analysis—investigation and visualization, H.K. and J.K.; project administration and supervision, J.K.; resources, J.-W.K.; software, H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the University Creative Asset Commercialization Support (BRIDEGE+) project from the Ministry of University Education (National Research Foundation of Korea).

Institutional Review Board Statement

Ethical review and approval were waived for this study because it belongs to the first class of the Korean medical device classification criteria, and the experiment was conducted with the consent of the subjects for research purposes, not commercial purposes.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patients to publish this paper.

Acknowledgments

We thank every experimental participant.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Parkinson’s Disease (PD); Multi-point IMU (MIMU); Single-point IMU (SIMU); Artificial Intelligence (AI); Machine Learning (ML); Gait Detection and Classification (GDC); Gait Detection and Classification Rate (GDCR); Gyroscopic Threshold (GT); Acceleration Threshold (AT).

References

  1. Pereira, A.P.S.; Marinho, V.; Gupta, D.; Magalhães, F.; Ayres, C.; Teixeira, S. Music Therapy and Dance as Gait Rehabilitation in Patients with Parkinson Disease: A Review of Evidence. J. Geriatr. Psychiatry Neurol. 2019, 32, 49–56. [Google Scholar] [CrossRef] [PubMed]
  2. Lim, I.; van Wegen, E.; de Goede, C.; Deutekom, M.; Nieuwboer, A.; Willems, A.; Jones, D.; Rochester, L.; Kwakkel, G. Effects of External Rhythmical Cueing on Gait in Patients with Parkinson’s Disease: A Systematic Review. Clin. Rehabil. 2005, 19, 695–713. [Google Scholar] [CrossRef] [PubMed]
  3. Kearney, E.; Shellikeri, S.; Martino, R.; Yunusova, Y. Augmented Visual Feedback-Aided Interventions for Motor Rehabilitation in Parkinson’s Disease: A Systematic Review. Disabil. Rehabil. 2019, 41, 995–1011. [Google Scholar] [CrossRef] [PubMed]
  4. Monteiro-Junior, R.S.; Cevada, T.; Oliveira, B.R.R.; Lattari, E.; Portugal, E.M.M.; Carvalho, A.; Deslandes, A.C. We Need to Move More: Neurobiological Hypotheses of Physical Exercise as a Treatment for Parkinson’s Disease. Med. Hypotheses 2015, 85, 537–541. [Google Scholar]
  5. van Duijnhoven, H.J.R.; Heeren, A.; Peters, M.A.M.; Veerbeek, J.M.; Kwakkel, G.; Geurts, A.C.H.; Weerdesteyn, V. Effects of Exercise Therapy on Balance Capacity in Chronic Stroke. Stroke 2016, 47, 2603–2610. [Google Scholar] [CrossRef] [Green Version]
  6. Tomlinson, C.L.; Patel, S.; Meek, C.; Herd, C.P.; Clarke, C.E.; Stowe, R.; Shah, L.; Sackley, C.M.; Deane, K.H.; Wheatley, K.; et al. Physiotherapy versus Placebo or No Intervention in Parkinson’s Disease. Cochrane Database Syst. Rev. 2013, CD002817. [Google Scholar] [CrossRef] [Green Version]
  7. Gunn, H.; Markevics, S.; Haas, B.; Marsden, J.; Freeman, J. Systematic Review: The Effectiveness of Interventions to Reduce Falls and Improve Balance in Adults with Multiple Sclerosis. Arch. Phys. Med. Rehabil. 2015, 96, 1898–1912. [Google Scholar] [CrossRef]
  8. Galea, J.M.; Mallia, E.; Rothwell, J.; Diedrichsen, J. The Dissociable Effects of Punishment and Reward on Motor Learning. Nat Neurosci 2015, 18, 597–602. [Google Scholar] [CrossRef] [Green Version]
  9. Trotter, A.B.; Inman, D.A. The Use of Positive Reinforcement in Physical Therapy. Phys. Ther. 1968, 48, 347–352. [Google Scholar] [CrossRef]
  10. Takeda, K.; Mani, H.; Hasegawa, N.; Sato, Y.; Tanaka, S.; Maejima, H.; Asaka, T. Adaptation Effects in Static Postural Control by Providing Simultaneous Visual Feedback of Center of Pressure and Center of Gravity. J. Physiol. Anthropol. 2017, 36, 31. [Google Scholar] [CrossRef] [Green Version]
  11. Mirelman, A.; Herman, T.; Nicolai, S.; Zijlstra, A.; Zijlstra, W.; Becker, C.; Chiari, L.; Hausdorff, J.M. Audio-Biofeedback Training for Posture and Balance in Patients with Parkinson’s Disease. J. NeuroEng. Rehabil. 2011, 8, 35. [Google Scholar] [CrossRef] [Green Version]
  12. Hasegawa, N.; Takeda, K.; Mancini, M.; King, L.A.; Horak, F.B.; Asaka, T. Differential Effects of Visual versus Auditory Bio-feedback Training for Voluntary Postural Sway. PLoS ONE 2020, 15, e0244583. [Google Scholar] [CrossRef] [PubMed]
  13. Ginis, P.; Nieuwboer, A.; Dorfman, M.; Ferrari, A.; Gazit, E.; Canning, C.G.; Rocchi, L.; Chiari, L.; Hausdorff, J.M.; Mirelman, A. Feasibility and Effects of Home-Based Smartphone-Delivered Automated Feedback Training for Gait in People with Parkinson’s Disease: A Pilot Randomized Controlled Trial. Park. Relat. Disord. 2016, 22, 28–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Dozza, M.; Chiari, L.; Chan, B.; Rocchi, L.; Horak, F.B.; Cappello, A. Influence of a Portable Audio-Biofeedback Device on Structural Properties of Postural Sway. J. NeuroEng. Rehabil. 2005, 2, 13. [Google Scholar] [CrossRef] [Green Version]
  15. Sienko, K.H.; Seidler, R.D.; Carender, W.J.; Goodworth, A.D.; Whitney, S.L.; Peterka, R.J. Potential Mechanisms of Sensory Augmentation Systems on Human Balance Control. Front. Neurol. 2018, 9, 944. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Mirelman, A.; Maidan, I.; Herman, T.; Deutsch, J.E.; Giladi, N.; Hausdorff, J.M. Virtual Reality for Gait Training: Can It Induce Motor Learning to Enhance Complex Walking and Reduce Fall Risk in Patients with Parkinson’s Disease? J. Gerontol. Ser. A 2011, 66A, 234–240. [Google Scholar] [CrossRef] [Green Version]
  17. Carpinella, I.; Cattaneo, D.; Bonora, G.; Bowman, T.; Martina, L.; Montesano, A.; Ferrarin, M. Wearable Sensor-Based Bio-feedback Training for Balance and Gait in Parkinson Disease: A Pilot Randomized Controlled Trial. Arch. Phys. Med. Rehabil. 2017, 98, 622–630.e3. [Google Scholar] [CrossRef] [Green Version]
  18. Dozza, M.; Chiari, L.; Peterka, R.J.; Wall, C.; Horak, F.B. What Is the Most Effective Type of Audio-Biofeedback for Postural Motor Learning? Gait Posture 2011, 34, 313–319. [Google Scholar] [CrossRef]
  19. Shen, X.; Mak, M.K.Y. Balance and Gait Training with Augmented Feedback Improves Balance Confidence in People with Parkinson’s Disease: A Randomized Controlled Trial. Neurorehabil. Neural Repair 2014, 28, 524–535. [Google Scholar] [CrossRef]
  20. Van den Heuvel, M.R.C.; Kwakkel, G.; Beek, P.J.; Berendse, H.W.; Daffertshofer, A.; van Wegen, E.E.H. Effects of Augmented Visual Feedback during Balance Training in Parkinson’s Disease: A Pilot Randomized Clinical Trial. Park. Relat. Disord. 2014, 20, 1352–1358. [Google Scholar] [CrossRef] [Green Version]
  21. Liao, Y.-Y.; Yang, Y.-R.; Wu, Y.-R.; Wang, R.-Y. Virtual Reality-Based Wii Fit Training in Improving Muscle Strength, Sensory Integration Ability, and Walking Abilities in Patients with Parkinson’s Disease: A Randomized Control Trial. Int. J. Gerontol. 2015, 9, 190–195. [Google Scholar] [CrossRef] [Green Version]
  22. Mhatre, P.V.; Vilares, I.; Stibb, S.M.; Albert, M.V.; Pickering, L.; Marciniak, C.M.; Kording, K.; Toledo, S. Wii Fit Balance Board Playing Improves Balance and Gait in Parkinson Disease. PM&R 2013, 5, 769–777. [Google Scholar]
  23. Baram, Y.; Aharon-Peretz, J.; Badarny, S.; Susel, Z.; Schlesinger, I. Closed-Loop Auditory Feedback for the Improvement of Gait in Patients with Parkinson’s Disease. J. Neurol. Sci. 2016, 363, 104–106. [Google Scholar] [CrossRef] [PubMed]
  24. Gooßes, M.; Saliger, J.; Folkerts, A.-K.; Nielsen, J.; Zierer, J.; Schmoll, P.; Niepold, A.; Colbach, L.; Leemhuis, J.; Engels, L.; et al. Feasibility of Music-Assisted Treadmill Training in Parkinson’s Disease Patients with and without Deep Brain Stimulation: Insights from an Ongoing Pilot Randomized Controlled Trial. Front. Neurol. 2020, 11, 790. [Google Scholar] [CrossRef] [PubMed]
  25. Cappozzo, A.; Della Croce, U.; Leardini, A.; Chiari, L. Human Movement Analysis Using Stereophotogrammetry: Part 1: Theoretical Background. Gait Posture 2005, 21, 186–196. [Google Scholar]
  26. Mobbs, R.J.; Phan, K.; Maharaj, M.; Rao, P.J. Physical Activity Measured with Accelerometer and Self-Rated Disability in Lumbar Spine Surgery: A Prospective Study. Glob. Spine J. 2016, 6, 459–464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Phan, K.; Mobbs, R.J. Long-Term Objective Physical Activity Measurements Using a Wireless Accelerometer Following Minimally Invasive Transforaminal Interbody Fusion Surgery. Asian Spine J. 2016, 10, 366–369. [Google Scholar] [CrossRef] [Green Version]
  28. Sprager, S.; Juric, M.B. Inertial Sensor-Based Gait Recognition: A Review. Sensors 2015, 15, 22089–22127. [Google Scholar] [CrossRef] [PubMed]
  29. Agostini, V.; Gastaldi, L.; Rosso, V.; Knaflitz, M.; Tadano, S. A Wearable Magneto-Inertial System for Gait Analysis (H-Gait): Validation on Normal Weight and Overweight/Obese Young Healthy Adults. Sensors 2017, 17, 2406. [Google Scholar] [CrossRef] [Green Version]
  30. Donath, L.; Faude, O.; Lichtenstein, E.; Pagenstert, G.; Nüesch, C.; Mündermann, A. Mobile Inertial Sensor Based Gait Analysis: Validity and Reliability of Spatiotemporal Gait Characteristics in Healthy Seniors. Gait Posture 2016, 49, 371–374. [Google Scholar] [CrossRef]
  31. Kluge, F.; Gaßner, H.; Hannink, J.; Pasluosta, C.; Klucken, J.; Eskofier, B.M. Towards Mobile Gait Analysis: Concurrent Validity and Test-Retest Reliability of an Inertial Measurement System for the Assessment of Spatio-Temporal Gait Parameters. Sensors 2017, 17, 1522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Maffiuletti, N.A.; Gorelick, M.; Kramers-de Quervain, I.; Bizzini, M.; Munzinger, J.P.; Tomasetti, S.; Stacoff, A. Concurrent Validity and Intrasession Reliability of the IDEEA Accelerometry System for the Quantification of Spatiotemporal Gait Parameters. Gait Posture 2008, 27, 160–163. [Google Scholar] [CrossRef] [PubMed]
  33. Orlowski, K.; Eckardt, F.; Herold, F.; Aye, N.; Edelmann-Nusser, J.; Witte, K. Examination of the Reliability of an Inertial Sensor-Based Gait Analysis System. Biomed. Eng./Biomed. Tech. 2017, 62, 615–622. [Google Scholar] [CrossRef]
  34. Rao, P.J.; Phan, K.; Maharaj, M.M.; Pelletier, M.H.; Walsh, W.R.; Mobbs, R.J. Accelerometers for Objective Evaluation of Physical Activity Following Spine Surgery. J. Clin. Neurosci. 2016, 26, 14–18. [Google Scholar] [CrossRef]
  35. Tunca, C.; Pehlivan, N.; Ak, N.; Arnrich, B.; Salur, G.; Ersoy, C. Inertial Sensor-Based Robust Gait Analysis in Non-Hospital Settings for Neurological Disorders. Sensors 2017, 17, 825. [Google Scholar] [CrossRef] [Green Version]
  36. Mobbs, R.J.; Katsinas, C.J.; Choy, W.J.; Rooke, K.; Maharaj, M. Objective Monitoring of Activity and Gait Velocity Using Wearable Accelerometer Following Lumbar Microdiscectomy to Detect Recurrent Disc Herniation. J. Spine Surg. 2018, 4, 792–797. [Google Scholar] [CrossRef]
  37. Demonceau, M.; Donneau, A.-F.; Croisier, J.-L.; Skawiniak, E.; Boutaayamou, M.; Maquet, D.; Garraux, G. Contribution of a Trunk Accelerometer System to the Characterization of Gait in Patients with Mild-to-Moderate Parkinson’s Disease. IEEE J. Biomed. Health Inform. 2015, 19, 1803–1808. [Google Scholar] [CrossRef]
  38. Van Schooten, K.S.; Pijnappels, M.; Rispens, S.M.; Elders, P.J.M.; Lips, P.; Daffertshofer, A.; Beek, P.J.; Dieën, J.H. van Daily-Life Gait Quality as Predictor of Falls in Older People: A 1-Year Prospective Cohort Study. PLoS ONE 2016, 11, e0158623. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Del Din, S.; Galna, B.; Godfrey, A.; Bekkers, E.M.J.; Pelosin, E.; Nieuwhof, F.; Mirelman, A.; Hausdorff, J.M.; Rochester, L. Analysis of Free-Living Gait in Older Adults with and Without Parkinson’s Disease and with and without a History of Falls: Identifying Generic and Disease-Specific Characteristics. J. Gerontol. Ser. A 2019, 74, 500–506. [Google Scholar] [CrossRef] [PubMed]
  40. Espinosa, H.G.; Thiel, D.V.; Sorell, M.; Rowlands, D. Can We Trust Inertial and Heart Rate Sensor Data from an APPLE Watch Device? Proceedings 2020, 49, 128. [Google Scholar]
  41. Mobbs, R.J.; Perring, J.; Raj, S.M.; Maharaj, M.; Yoong, N.K.M.; Sy, L.W.; Fonseka, R.D.; Natarajan, P.; Choy, W.J. Gait Metrics Analysis Utilizing Single-Point Inertial Measurement Units: A Systematic Review. Mhealth 2022, 8, 9. [Google Scholar] [CrossRef] [PubMed]
  42. Prasanth, H.; Caban, M.; Keller, U.; Courtine, G.; Ijspeert, A.; Vallery, H.; von Zitzewitz, J. Wearable Sensor-Based Real-Time Gait Detection: A Systematic Review. Sensors 2021, 21, 2727. [Google Scholar] [CrossRef]
  43. Semwal, V.B.; Gaud, N.; Lalwani, P.; Bijalwan, V.; Alok, A.K. Pattern Identification of Different Human Joints for Different Human Walking Styles Using Inertial Measurement Unit (IMU) Sensor. Artif. Intell. Rev. 2022, 55, 1149–1169. [Google Scholar] [CrossRef]
  44. Sarshar, M.; Polturi, S.; Schega, L. Gait Phase Estimation by Using LSTM in IMU-Based Gait Analysis—Proof of Concept. Sensors 2021, 21, 5749. [Google Scholar] [CrossRef]
  45. Trabassi, D.; Serrao, M.; Varrecchia, T.; Ranavolo, A.; Coppola, G.; De Icco, R.; Tassorelli, C.; Castiglia, S.F. Machine Learning Approach to Support the Detection of Parkinson’s Disease in IMU-Based Gait Analysis. Sensors 2022, 22, 3700. [Google Scholar] [CrossRef] [PubMed]
  46. Mundt, M.; Koeppe, A.; David, S.; Witter, T.; Bamer, F.; Potthast, W.; Markert, B. Estimation of Gait Mechanics Based on Simulated and Measured IMU Data Using an Artificial Neural Network. Front. Bioeng. Biotechnol. 2020, 8, 41. [Google Scholar] [CrossRef] [PubMed]
  47. Balakrishnan, A.; Medikonda, J.; Namboothiri, P.K.; Natarajan, M. Role of Wearable Sensors with Machine Learning Approaches in Gait Analysis for Parkinson’s Disease Assessment: A Review. Eng. Sci. 2022, 19, 5–19. [Google Scholar] [CrossRef]
  48. Pappas, I.P.; Popovic, M.R.; Keller, T.; Dietz, V.; Morari, M. A reliable gait phase detection system. IEEE Trans. Neural Syst. Rehabil. Eng. 2001, 9, 113–125. [Google Scholar] [CrossRef]
  49. Mazhar, O.; Bari, A.Z.; Faudzi, A.A.M. Real-time gait phase detection using wearable sensors. In Proceedings of the 2015 10th Asian Control Conference (ASCC), Kota Kinabalu, Malaysia, 31 May–3 June 2015; pp. 1–4. [Google Scholar]
  50. Li, J.; Zhou, X.; Li, C.; Li, W.; Zhang, H.; Gu, H. A Real-Time Gait Phase Detection Method for Prosthesis Control. In Assistive Robotics, Proceedings of the 18th International Conference on CLAWAR 2015, Hangzhou, China, 6–9 September 2015; National Natural Science Foundation of China: Beijing, China, 2016; pp. 577–584. [Google Scholar]
  51. Goršič, M.; Kamnik, R.; Ambrožič, L.; Vitiello, N.; Lefeber, D.; Pasquini, G.; Munih, M. Online phase detection using wearable sensors for walking with a robotic prosthesis. Sensors 2014, 14, 2776–2794. [Google Scholar] [CrossRef]
  52. Yang, J.; Chen, X.; Guo, H.; Zhang, Q. Implementation of omnidirectional lower limbs rehabilitation training robot. In Proceedings of the 2007 International Conference on Electrical Machines and Systems (ICEMS), Seoul, Republic of Korea, 8–11 October 2007; pp. 2033–2036. [Google Scholar]
  53. Kim, H.; Kim, J.-W.; Ko, J. Gait Disorder Detection and Classification Method Using Inertia Measurement Unit for Augmented Feedback Training in Wearable Devices. Sensors 2021, 21, 7676. [Google Scholar] [CrossRef]
  54. Karl, K.A.; O’Leary-Kelly, A.M.; Martocchio, J.J. The Impact of Feedback and Self-Efficacy on Performance in Training. J. Organ. Behav. 1993, 14, 379–394. [Google Scholar] [CrossRef]
  55. Burgers, C.; Eden, A.; van Engelenburg, M.D.; Buningh, S. How Feedback Boosts Motivation and Play in a Brain-Training Game. Comput. Hum. Behav. 2015, 48, 94–103. [Google Scholar] [CrossRef] [Green Version]
  56. Wilson, K.M.; Helton, W.S.; de Joux, N.R.; Head, J.R.; Weakley, J.J.S. Real-Time Quantitative Performance Feedback during Strength Exercise Improves Motivation, Competitiveness, Mood, and Performance. Proc. Hum. Factors Ergon. Soc. Annu. Meet. 2017, 61, 1546–1550. [Google Scholar] [CrossRef]
  57. Burke, J.W.; McNeill, M.D.J.; Charles, D.K.; Morrow, P.J.; Crosbie, J.H.; McDonough, S.M. Optimising Engagement for Stroke Rehabilitation Using Serious Games. Vis. Comput. 2009, 25, 1085. [Google Scholar] [CrossRef]
  58. Iosa, M.; Capodaglio, E.; Pelà, S.; Persechino, B.; Morone, G.; Antonucci, G.; Paolucci, S.; Panigazzi, M. Artificial Neural Network Analyzing Wearable Device Gait Data for Identifying Patients with Stroke Unable to Return to Work. Front. Neurol. 2021, 12, 650542. [Google Scholar] [PubMed]
  59. Mohan, D.M.; Khandoker, A.H.; Wasti, S.A.; Ismail Ibrahim Ismail Alali, S.; Jelinek, H.F.; Khalaf, K. Assessment Methods of Post-Stroke Gait: A Scoping Review of Technology-Driven Approaches to Gait Characterization and Analysis. Front. Neurol. 2021, 12, 650024. [Google Scholar] [CrossRef]
  60. Nadeau, S.; Betschart, M.; Bethoux, F. Gait Analysis for Poststroke Rehabilitation: The Relevance of Biomechanical Analysis and the Impact of Gait Speed. Phys. Med. Rehabil. Clin. 2013, 24, 265–276. [Google Scholar] [CrossRef] [PubMed]
  61. Scheffer, C.; Cloete, T. Inertial Motion Capture in Conjunction with an Artificial Neural Network Can Differentiate the Gait Patterns of Hemiparetic Stroke Patients Compared with Able-Bodied Counterparts. Comput. Methods Biomech. Biomed. Eng. 2012, 15, 285–294. [Google Scholar] [CrossRef]
  62. Lau, H.; Tong, K.; Zhu, H. Support Vector Machine for Classification of Walking Conditions of Persons after Stroke with Dropped Foot. Hum. Mov. Sci. 2009, 28, 504–514. [Google Scholar] [CrossRef]
  63. Zhou, Y.; Romijnders, R.; Hansen, C.; van Campen, J.; Maetzler, W.; Hortobágyi, T.; Lamoth, C.J.C. The Detection of Age Groups by Dynamic Gait Outcomes Using Machine Learning Approaches. Sci. Rep. 2020, 10, 4426. [Google Scholar] [CrossRef] [Green Version]
  64. Hutabarat, Y.; Owaki, D.; Hayashibe, M. Quantitative Gait Assessment with Feature-Rich Diversity Using Two IMU Sensors. IEEE Trans. Med. Robot. Bionics 2020, 2, 639–648. [Google Scholar] [CrossRef]
  65. Zhou, L.; Tunca, C.; Fischer, E.; Brahms, C.M.; Ersoy, C.; Granacher, U.; Arnrich, B. Validation of an IMU Gait Analysis Algorithm for Gait Monitoring in Daily Life Situations. In Proceedings of the 2020 42nd Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC), Montreal, QC, Canada, 20–24 July 2020; pp. 4229–4232. [Google Scholar]
  66. Gujarathi, T.; Bhole, K. Gait Analysis Using Imu Sensor. In Proceedings of the 2019 10th International Conference on Computing, Communication and Networking Technologies (ICCCNT), Kanpur, India, 6–8 July 2019; pp. 1–5. [Google Scholar]
  67. Anwary, A.R.; Yu, H.; Vassallo, M. An Automatic Gait Feature Extraction Method for Identifying Gait Asymmetry Using Wearable Sensors. Sensors 2018, 18, 676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Definition of elements used in previous research to detect and classify a user’s gait.
Figure 1. Definition of elements used in previous research to detect and classify a user’s gait.
Sensors 23 06638 g001
Figure 2. Figure illustrating the role of the AT and GT. (a) The case when a gait was not detected (b) The case when a gait was detected and classified as a good gait (c) The case when a gait was detected but not classified as a good gait.
Figure 2. Figure illustrating the role of the AT and GT. (a) The case when a gait was not detected (b) The case when a gait was detected and classified as a good gait (c) The case when a gait was detected but not classified as a good gait.
Sensors 23 06638 g002
Figure 3. Example results of previous research: (a) Results from non-PD subjects; and (b) Results from PD subjects.
Figure 3. Example results of previous research: (a) Results from non-PD subjects; and (b) Results from PD subjects.
Sensors 23 06638 g003
Figure 4. Equipment image and experimental setup: (a) equipment example; and (b) experimental setup.
Figure 4. Equipment image and experimental setup: (a) equipment example; and (b) experimental setup.
Sensors 23 06638 g004
Figure 5. Shape of average GDCR and FG. Solid line curve is the average GDCR from experimental dataset. Dashed line curve which has GDCR range from 0 to 1 is normalized fitted curve from average GDCR curve. Gray-colored area means the average GDCR between slowest (line below of the average GDCR, 1.0 km/h) and fastest (line above of the average GDCR, 3.0 km/h) walking speed.
Figure 5. Shape of average GDCR and FG. Solid line curve is the average GDCR from experimental dataset. Dashed line curve which has GDCR range from 0 to 1 is normalized fitted curve from average GDCR curve. Gray-colored area means the average GDCR between slowest (line below of the average GDCR, 1.0 km/h) and fastest (line above of the average GDCR, 3.0 km/h) walking speed.
Sensors 23 06638 g005
Figure 6. (a) FG from the experimental results when GT was 15,000 and various sigmoid functions; and (b) average error of sigmoid functions from FG.
Figure 6. (a) FG from the experimental results when GT was 15,000 and various sigmoid functions; and (b) average error of sigmoid functions from FG.
Sensors 23 06638 g006
Figure 7. Adaptive control algorithm flowchart.
Figure 7. Adaptive control algorithm flowchart.
Sensors 23 06638 g007
Figure 8. Various exercise scenarios and walking speeds. Each scenario is a combination of scenario elements: linear increase, linear decrease, random with upper and lower boundaries, sudden increase, and sudden decrease. The exercise procedure means exercise from start to end. #1: Linear increase. #2: Linear decrease. #3: Random with upper and lower boundaries. #4: Sudden decrease during linear increase. #5: Sudden increase during linear decrease. #6: Linear increase after linear decrease. #7: Linear decrease—random with upper and lower boundaries—linear increase. #8: Linear decrease—sudden increase—random with upper and lower boundaries. #9: Random with upper and lower boundaries after linear increase.
Figure 8. Various exercise scenarios and walking speeds. Each scenario is a combination of scenario elements: linear increase, linear decrease, random with upper and lower boundaries, sudden increase, and sudden decrease. The exercise procedure means exercise from start to end. #1: Linear increase. #2: Linear decrease. #3: Random with upper and lower boundaries. #4: Sudden decrease during linear increase. #5: Sudden increase during linear decrease. #6: Linear increase after linear decrease. #7: Linear decrease—random with upper and lower boundaries—linear increase. #8: Linear decrease—sudden increase—random with upper and lower boundaries. #9: Random with upper and lower boundaries after linear increase.
Sensors 23 06638 g008
Figure 9. Example of virtual simulation with scenarios. If the walking speed changed from 1.0 km/h, 2.0 km/h, and to 3.0 km/h, the example simulation in scenario one will follow the sequence as shown in this figure.
Figure 9. Example of virtual simulation with scenarios. If the walking speed changed from 1.0 km/h, 2.0 km/h, and to 3.0 km/h, the example simulation in scenario one will follow the sequence as shown in this figure.
Sensors 23 06638 g009
Figure 10. Average GDCR according to legs used and walking speeds. The range of AT was from 1000 to 15,000: (a) The relation between the GDCR and used legs (n = 246; *** p < 0.001); and (b) relation between the GDCR and walking speeds (n = 246; *** p < 0.001).
Figure 10. Average GDCR according to legs used and walking speeds. The range of AT was from 1000 to 15,000: (a) The relation between the GDCR and used legs (n = 246; *** p < 0.001); and (b) relation between the GDCR and walking speeds (n = 246; *** p < 0.001).
Sensors 23 06638 g010
Figure 11. Average GDCR according to the GT and walking speeds. Entire GT range (from 15,000 to 25,000): (a) shape of the GDCR curve according to the GT; and (b) relation between Spearman’s correlation coefficients and walking speeds (n = 441; *** p < 0.001 vs. GT, for each walking speed. In sum, n = 2205).
Figure 11. Average GDCR according to the GT and walking speeds. Entire GT range (from 15,000 to 25,000): (a) shape of the GDCR curve according to the GT; and (b) relation between Spearman’s correlation coefficients and walking speeds (n = 441; *** p < 0.001 vs. GT, for each walking speed. In sum, n = 2205).
Sensors 23 06638 g011
Figure 12. Average GDCR according to the AT and walking speeds: (a) shape of the GDCR curve according to the AT. AT more than 15,000 was trimmed in the figure since there were no significant changes in values; and (b) relation between Spearman’s correlation coefficients and walking speeds. (n = 39,921; *** p < 0.001 vs. AT, for each walking speed. In sum, n = 199,605).
Figure 12. Average GDCR according to the AT and walking speeds: (a) shape of the GDCR curve according to the AT. AT more than 15,000 was trimmed in the figure since there were no significant changes in values; and (b) relation between Spearman’s correlation coefficients and walking speeds. (n = 39,921; *** p < 0.001 vs. AT, for each walking speed. In sum, n = 199,605).
Sensors 23 06638 g012
Figure 13. The average error of each control method. (n = 225; *** p < 0.001, +++ p < 0.001, vs. adaptive, ++ p < 0.01, vs. adaptive). Errors from the GDCR control algorithm came from a narrower AT range than other results, e.g., those of Figure 10.
Figure 13. The average error of each control method. (n = 225; *** p < 0.001, +++ p < 0.001, vs. adaptive, ++ p < 0.01, vs. adaptive). Errors from the GDCR control algorithm came from a narrower AT range than other results, e.g., those of Figure 10.
Sensors 23 06638 g013
Figure 14. The average error of each control method on every scenario. The N of cases was 25 per each control method in one scenario. The N of cases in one scenario was 75 in total. Errors from GDCR control algorithm came from a narrower AT range, not the full range.
Figure 14. The average error of each control method on every scenario. The N of cases was 25 per each control method in one scenario. The N of cases in one scenario was 75 in total. Errors from GDCR control algorithm came from a narrower AT range, not the full range.
Sensors 23 06638 g014
Table 1. Statistical analyses were used to determine the impact of elements used in the GDC algorithm.
Table 1. Statistical analyses were used to determine the impact of elements used in the GDC algorithm.
RelationVariableNumber of SubjectsNumber of SamplesStatistical Analysis Method
GDCR-LegLeft leg25121Independent sample
Mann–Whitney U test
Right leg25125
GDCR-Walking speed1.0 km/h2547Independent sample
Kruskal–Wallis test
1.5 km/h2550
2.0 km/h2549
2.5 km/h2550
3.0 km/h2550
GDCR-GT1.0 km/h2539,921Spearman’s correlation analysis
1.5 km/h2539,921
2.0 km/h2539,921
2.5 km/h2539,921
3.0 km/h2539,921
GDCR-AT1.0 km/h25441
1.5 km/h25441
2.0 km/h25441
2.5 km/h25441
3.0 km/h25441
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, H.; Kim, J.-W.; Ko, J. Adaptive Control Method for Gait Detection and Classification Devices with Inertial Measurement Unit. Sensors 2023, 23, 6638. https://doi.org/10.3390/s23146638

AMA Style

Kim H, Kim J-W, Ko J. Adaptive Control Method for Gait Detection and Classification Devices with Inertial Measurement Unit. Sensors. 2023; 23(14):6638. https://doi.org/10.3390/s23146638

Chicago/Turabian Style

Kim, Hyeonjong, Ji-Won Kim, and Junghyuk Ko. 2023. "Adaptive Control Method for Gait Detection and Classification Devices with Inertial Measurement Unit" Sensors 23, no. 14: 6638. https://doi.org/10.3390/s23146638

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