Next Article in Journal
Methylation Genome-Wide Profiling in Lowly and Highly Efficient Somatic Cell Nuclear Transfer in Pigs
Previous Article in Journal
Special Issue on Metal-Based Composite Materials: Preparation, Structure, Properties and Applications
Previous Article in Special Issue
Evaluations and Measurements of the Occurrence of Maxilla and Palatine Bone Asymmetry Based on 3D Printed Stereolithographic Models in Elderly Edentulous People
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Stride Length, Foot Clearance, and Foot Progression Angle Using UWB Sensors

1
Center for Bionics, Korea Institute of Science and Technology, Seoul 02792, Republic of Korea
2
Department of Mechanical Engineering, Korea University, Seoul 02841, Republic of Korea
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2023, 13(8), 4801; https://doi.org/10.3390/app13084801
Submission received: 28 February 2023 / Revised: 6 April 2023 / Accepted: 10 April 2023 / Published: 11 April 2023
(This article belongs to the Special Issue Sports Science, Medicine and Rehabilitation)

Abstract

:
Stride length (SL), foot clearance (FC), and foot progression angle (FPA) are the key parameters for diagnosing gait disorders. This study used the distance data between two feet measured by ultra-wideband (UWB) sensors installed on shoes and proposed a method for estimating the three gait parameters. Here, a method of compensating the offset of the UWB sensor and estimating the distances between a base sensor installed on one foot during the stance phase and three UWB sensors on the other during the swing phase was applied. Foot trajectory was acquired in a gait experiment with ten healthy adults walking on a treadmill. The results were compared with those obtained using a motion capture system (MCS). The UWBs sensor displayed average errors of 45.84 mm, 7.60 mm, and 2.82° for SL, FC, and FPA, respectively, compared with the MCS. A similar accuracy level was achieved in a previous study that used an inertial measurement unit (IMU). Thus, these results suggest that UWB sensors can be extensively applied to sensor systems used to analyze mobile gait systems.

1. Introduction

Gait analysis assesses and measures the temporal and spatial characteristics of the body during motion. Gait characteristics data can be used as diagnostic parameters for various adult diseases, resulting in numerous clinical researchers focusing on developing technologies for measuring gait characteristics [1,2]. Among gait parameters, stride length (SL), gait speed (GS), and foot clearance (FC) are used to determine the risk of falling [3,4,5]. In addition, the foot progression angle (FPA) could be applied to examine the characteristics of the load applied to the knee joint while assessing degenerative arthritis [6,7,8,9].
Clinicians have reported that the gait asymmetry problem caused by slow GS [10] or low FC [11,12] could be improved through visual biofeedback. In addition, a study that fed the FPA to a subject walking on a treadmill confirmed that the subject improved the FPA on their own [13]. However, these studies were inapplicable in their daily life because they were mostly conducted in experimental environments that used motion capture systems (MCSs). Therefore, numerous studies have analyzed gait parameters using wearable devices. An inertial measurement unit (IMU) is commonly used for gait analysis with wearable devices. The SL [14,15,16], GS [17,18], gait phase [19,20,21], joint angle [22,23], and FC [24,25,26] can be estimated using these devices. The IMU can also be used simultaneously with an ultrasonic sensor to measure the step length and stride width [27].
The disadvantage of using an IMU is the drift phenomenon, which is the occurrence of a cumulative error during speed acquisition by integrating acceleration data. Therefore, researchers conducting mobile gait analysis have improved the accuracy by applying the zero-velocity update (ZUPT) method [28], which corrects the sensor measurement by constraining the speed to zero if the sensor is determined to be at rest owing to the set threshold value. However, although ZUPT may be effective under circumstances of a short activity time and repeated state of rest, it does not function when: (1) the stance phase only occurs for a very short period, such as while fast walking or running; and (2) the sensor used to detect zero velocity is affected by noise and the drift is not adequately corrected [29,30,31].
Ultra-wideband (UWB) is an alternative to the IMU with the abovementioned problems. UWB sensors determine distance using the time-of-flight (TOF) between two UWB sensors based on the reception and transmission of radio signals, making it highly advantageous compared with IMU because of its capability to measure distances directly. UWB technology has mainly been used to track locations within indoor spaces using fixed UWB anchor [32,33]. However, as the power consumption of UWB is reduced [34] and the measurement accuracy is increased [35] owing to recent technological developments, the utilization of UWB for gait parameter analysis, such as tracking the joint angle [36,37] and measuring the SL [38], is increasing.
The limitations of previous studies using UWB sensors are as follows: (1) since the signal receiver is installed in a fixed location, such as a wall or tripod, it can only operate within a designated place; and (2) gait analysis is possible only in the double stance phase.
To overcome these limitations, this study proposes a UWB sensor system (Figure 1) that tracks foot trajectory using four UWB sensors installed on each foot to analyze gait characteristics. First, the trilateration method was applied, which involved simultaneously applying the distance data acquired from each UWB sensor to determine the overlapping location and estimate the actual distance between the sensors. However, the distance data provided by the UWB sensors had errors owing to variations in the antenna direction. Therefore, the variation in the strength of the received sensor signals was considered to compensate for this error. The SL, FC, and FPA data for clinical assessments were determined by conducting a gait experiment using the developed UWB sensor system. Simultaneously, an MCS was applied to consider the data as the gold standard and review the measurement accuracy of the UWB sensor system.

2. Materials and Methods

2.1. Hardware Description

As shown in Figure 2a, the equipment used in this experiment is a UWB system, comprising UWB sensors attached to each shoe to measure the distance (collected on a personal computer (PC)), an MCS that measured the location data of the reflection marker attached to the shoe, and a device that synchronized the data of the two systems.
The UWB system included a UWB sensor (model: DWM1000, Qorvo, Greensboro, NC, USA) and an Arduino Nano. The UWB sensor was installed on an in-house printed circuit board (PCB) and fabricated to have a total size of 43 mm (width) × 20 mm (length) × 15 mm (height) and a data acquisition speed of 100 Hz. The trilateration method involves attaching three anchors and a tag to a shoe, as shown in Figure 2b, using eight UWB sensors (four in each shoe). The UWB sensors were placed between the two shoes facing each other so that wireless signals could be sent and received without interruption. As illustrated in Figure 2c, the tag on each foot measures the distance between the three anchors located on the opposite foot. The data were transmitted to and stored on a PC connected to a wire and used in the analysis. Gait data were measured using the MCS (Vicon Nexus MX, Oxford Metric, Oxford, UK) installed in the gait experiment space. The data measured by MCS were considered as the gold standard and were applied to evaluate the trajectory estimated by the UWB system. A reflection marker for motion capture was attached to the PCB for the UWB sensors, as shown in Figure 2a,b. The location of the reflective marker was placed at a site that does not affect the antenna signal of the UWB sensor. Meanwhile, a push switch (5 V) for triggers was attached to the right-foot tag, and its signal was used to synchronize the UWB data obtained from the gait experiment with the motion capture data before comparing them (Figure 2a).
The gait experiment proceeded with the participant walking on a treadmill (M-gait, Motek, Amsterdam, The Netherlands) wearing shoes equipped with UWB sensors and MCS markers.

2.2. Participants and Test Method

The participants in the gait experiment were 10 healthy adults (5 males and five females) in their 20 s and 30 s (Table 1). They were devoid of neurological diseases and gait disorders and were capable of normal gait. The experimental protocol was approved by the Institutional Review Board (IRB) of the Korea Institute of Science and Technology (KIST). All participants provided written informed consent for participation in the study.
Each participant installed a UWB system and MCS marker on a shoe that fit their size and participated in the experiment. The participants ascended the treadmill and stood waiting while gazing toward the front until they received a cue from the test performer. When they received the signal, they walked at designated speeds (2, 3, and 4 km/h) for 3 min and stopped the gait according to the stopping speed. The conductor first operated the trigger for data synchronization, then started the treadmill at the designated speed, along with the cue sign. When the walking time exceeded three minutes, the conductor stopped the trigger, declared the end of the test to the subjects, and stopped the treadmill. The first five steps after the trigger operation and the final five steps before the treadmill stopped were excluded from the collected gait experiment data during data analysis to increase accuracy.

2.3. Calculation of Gait Parameters Using UWB

This section describes the method for estimating gait parameters using a UWB sensor. The gait characteristic analysis comprises three steps, as shown in Figure 3a. The three gait parameters, SL, CL, and FPA (see Figure 3b,c), were used in the analysis.

2.3.1. UWB Sensor Offset Compensation

The UWB sensor used in this study performs distance measurements using the two-way ranging (TWR) method [39]. In this method, (1) the TOF, as the time required to transmit and receive wireless signals between the two directions of the anchor (as the fixed antenna) and tag (as the mobile antenna) is determined; and (2) the distance between the two sensors is measured. According to the UWB sensor datasheet (model: DWM1000), the precision of the sensor resolution was estimated to be 10 cm. However, the resolution of the oscillator mounted on the UWB was 16 ps [38], and according to the latest research, it has been proven that the UWB sensor can measure more stably at a short distance than the proposed precision. However, this sensor can experience signal delays according to variations in the surrounding temperature and antenna direction. Therefore, most studies to resolve this issue used the offset compensation method provided by Qorvo Inc. (Greensbor, NC, USA) [40].
The method used by previous researchers involved the variation in the distance between the anchor and tag terminals; measurement of the distance, signal strength, and temperature; and generation of an offset compensation function using the signal strength and temperature used to measure the distance. This method may improve the accuracy of the UWB sensor according to offset compensation in a laboratory environment. However, if the compensation method assumes that the anchor and tag face each other, it is infeasible to compensate for the offset according to the variation in the direction of the antenna. In reality, the offset varies according to the variation in the relative angle between the two UWB sensors [41], and the antenna angle varies periodically while walking, which causes the occurrence of a residual error after offset compensation.
This research used support vector regression (SVR) to compensate for the residual error generated by the variation in the antenna angle according to the gait cycle. Therefore, before the gait characteristics experiment, the following procedure was conducted for the distance data-offset compensation.
  • Before the treadmill gait experiment, the participants wore shoes equipped with UWB sensors and walked on the treadmill for 5 min. In this experiment, the signal strength and distance data measured by the UWB sensor according to the gait cycle before offset compensation was measured, and the location of the motion capture marker next to the UWB sensor was acquired.
  • The function that compensates for the offset of the sensor consists of SVR, which is composed of the signal strength and distance data measured using the UWB sensor, and the target variable Y, which is composed of the motion capture marker data. The composed dataset D is expressed by Equation (1). The UWB distance data used in the dataset were obtained after applying a third-order Savitzky–Golay filter [38].
    x i , y i i = 1 ι , x i R N , y i 1 ,   1
  • To distinguish different classes using input data, the target variable should be mapped to the domain of a higher dimension, and the hyperplane is expressed as a shape space vector with mapping ψ to the shape space Z. The linear decision function composed of the weight vector and the bias vector pair ω ,   b is as follows:
    f x = w T x + b
  • Weight and bias can be estimated by minimizing the Euclidean norm. The minimization equation composed of cost function C and slack variable ξ is as follows:
    Minimize   1 2 w 2 + C i = 1 n ξ i + ξ i *
    subject   to   y i w · z i + b 1 ξ i
  • Equation (3) solves this issue using positive Lagrange multipliers. The corresponding equations are as follows:
    Maximize   W α = i = 1 l α l 1 2 i = 1 l j = 1 l α i α j y i y j z i · z j
    subject   to   i = 1 l y i α i = 0 ,   0 α i C ,   i = 1 ,   2 ,   ,   l
  • Finally, the kernel function was used to map the data space to the shape space. In addition, a linear kernel is used to compensate for the offset of the data acquired from the UWB sensor. Previous research results [42] can be referred to for further details on the SVR.

2.3.2. Foot Trajectory Estimation

The UWB sensor, after compensation, was installed on the shoe, as shown in Figure 2. The location coordinates of the tag were estimated using the gait experiment. In general, the location coordinates are estimated from the relative location of the tag, based on fixed anchor coordinates. However, the coordinates of the anchor cannot be fixed in mobile-type smart insoles. Therefore, for each step in this study, the location of the foot tag in the swing phase was calculated based on the anchor on the foot in the stance phase.
To determine the location of the foot, the coordinate system was first defined based on the anchor located on the big toe ( A n c h o r t o e ) and that located at the heel of the foot ( A n c h o r h e e l ) in the stance phase. Here, the location coordinates of the tag, defined as (x, y, z), were A n c h o r h e e l as (0, 0, 0), and that of the anchor located at the three-fourths point of the shoe length ( A n c h o r q u a t e r ) as (0, y Q , z Q ), and that of X as (0, y T , 0). Then, the distance between each anchor and tag can be expressed as follows:
d T 2 = x 2 + y y T 2 + z 2
d Q 2 = x 2 + y y Q 2 + z z Q 2
d H 2 = x 2 + y 2 + z 2
The above equations can be evaluated to find the tag coordinates (x, y, z), which are expressed as follows:
y = d H 2 d T 2 + y T 2 2 y T
z = d H 2 d Q 2 + y Q 2 + z Q 2 2 y y Q 2 z Q
x = ± d H 2 y 2 z 2
In Equation (10), the positive and negative signs are obtained when the location coordinates of the right and left feet are estimated according to the reference coordinate.

2.3.3. Stride Length Estimation

SL refers to the linear distance between the coordinate at the time of starting the swing and at the ending of swing from the foot trajectory measured from one stride (Figure 3b). This is determined as follows:
SL = x S 1 x S 2 2 + y S 1 y S 2 2 + z S 1 z S 2 2
x S 1 , y S 1 , z S 1 = argmax y i , x S 2 , y S 2 , z S 2 = argmin y i , i = 1 , 2 , 3 , ,   n

2.3.4. Foot Clearance Estimation

FC is the value with the largest magnitude among the foot trajectory z data measured from one stride.
F C = max z i , i = 1 ,   2 ,   3 ,   ,   n

2.3.5. Foot Progression Angle Estimation

The FPA is calculated as the sum of the angle between the foot direction vector and d H vector and the walking direction vector, which can be obtained from the coordinates of the start and end points of the swing (Figure 4), and d H vector. This is calculated as follows:
θ w = cos 1 S T n 2 + d H n + 1 2 d H n 2 2 S T n d H n + 1
θ F = cos 1 y T 2 + d H n + 1 2 d T n + 1 2 2 y T d H n + 1
FPA = θ w + θ T π

3. Results

3.1. Sensor Offset Compensation

Examples of the motion capture data acquired from the pre-processing experiment to compensate for the UWB sensor data are shown in Figure 5. In Figure 5, the blue line represents the motion capture data used as the gold standard, the black dotted line represents the UWB sensor data, and the red line represents the UWB data that underwent offset compensation using SVR. As shown in Figure 5, compensating for the UWB data using SVR reduces the error. This value was similar to that of the MCS after compensation.
The offset compensation method applied in this research was used to compensate for the offset of the data acquired from six anchors, namely, the A n c h o r h e e l , A n c h o r q u a t e r , and A n c h o r t o e of the right and left foot. These showed mean square errors (MSE) of 10.27, 24.77, 9.27, 8.32, 10.62, and 9.91 mm, respectively. Thus, the average MSE was 12.19 mm.

3.2. Results for Gait Parameters

Figure 6 shows the UWB data acquired from participant #3 of the gait experiment and the 3D trajectory of the right foot determined using Equations (8)–(10). The foot trajectory estimated using UWB is shown in Figure 6a. Here, a closed curve is drawn using the combined trajectory of the foot when the tag is in the swing and stance phases. Herein, after the heel strike point, when the foot is in a swing state and higher than the opposite foot, the z-direction location value reverses as the opposite foot lifts from the ground. Therefore, the foot’s trajectory at a wholesome swing phase could be illustrated only if the data with the z value of the anchor of the opposite foot is higher than the z value data of the tag location installed on the stepping foot in the stance phase. This is illustrated in Figure 6b.
Figure 7 shows the Bland–Altman plot that assesses SL, FC, and FPA using the 3D trajectory data of the foot during the swing phase. It shows the data acquired from all participants. The x-axis represents the average of the values measured by the MSC (or optical reference) and the UWB system, and the y-axis represents the difference between the two measured values (MCS (or optical reference)—UWB).
It can be seen that the SL data in Figure 7a is divided into two groups, and the average SL values of each group are 719 mm and 874 mm, respectively. The data near SL = 719 mm, were mainly obtained from a group of female participants with heights < 170 cm. The data near SL = 874 mm were mainly from a mixed group of participants with heights ≥ 170 cm. Figure 7 shows that the obtained data were mostly within 95% confidence interval. Therefore, the foot trajectory data estimated using UWB could be considered to be in good agreement with the data measured using MCS.
The data in Figure 7 were analyzed and are presented in Table 2. The differences in SL at speeds of 2, 3, and 4 km/h were 45.13, 46.99, and 45.40 mm, respectively. Furthermore, the average difference in the entire gait speed section was 45.84 mm. The differences in FC at speeds of 2, 3, and 4 km/h were 7.80, 7.45, and 7.54 mm, respectively, and the average difference for the entire gait speed section was 7.60 mm. The differences in the FPA for the same speeds were 2.53°, 2.86°, and 3.12°, respectively, and the average difference for the entire gait speed section was 2.82°. In Figure 8, the average SL, FC, and FPA were measured for each gait speed. Moreover, a 95% confidence interval was identified. Overall, the difference in SL was trivially impacted by the variation in gait speed. In addition, the difference in FC was relatively large when the gait speed was low and tended to converge to a certain value as the gait speed increased. The difference in FPA increased continuously as the gait speed increased, although the amount of variation was marginal.
The estimated gait parameters from UWB and MCS were compared: they show good correlation coefficients (r2 > 0.7) of 0.898 in SL, 0.761 in CL, and 0.743 in FPA.

4. Discussion

Table 3 compares the SL, FC, and FPA data obtained in this study (expressed as KIST) with data acquired in previous studies that used an IMU.
The SL estimation error was 45.84 mm in this study, compared to 54 mm [43] and (−) 27 mm [44] in previous studies. However, two previous studies applied IMU sensors to the shank for gait experiments. Therefore, the measurement error for the SL value acquired from this research was considered to be in a range similar to that of the previous research that used an IMU.
FC estimation error was 7.60 mm in this study, compared to 4.2 mm/2.7 mm [45], and 8.6 mm [46] in previous studies. In a study by Fan et al. [45], the FS estimation error varied according to the attachment location of the IMU. The heel clearance and toe clearance measurement errors were 7.0 mm and 2.7 mm, respectively, when the IMU was located at the toes, and 4.2 mm and 13.1 mm, respectively, when located at the heel. This phenomenon was similar to that observed by Huang et al. [46]. These results were attributed to the identical foot experiencing different accelerations according to the part within the three-dimensional movement of the foot. Therefore, the detected foot trajectory was considered to differ according to the installation location of the IMU on the foot. This was also a significant error factor when tracking the foot trajectory using the IMU. In addition, the FC measurement error in this study was within a range similar to that used in previous studies that used an IMU.
The FPA estimation error was 2.82° in this study, compared to 2.6° [31] and 2.5° [29] in previous studies. In a study by Wouda et al. [31], the FPA measurement error was observed to increase to 5.22° when the sensor location was misarranged during the positioning of the coordinate system of the IMU attached to the foot. In summary, the error that occurred during the process of applying ZUPT and arranging the sensor coordinate system in the gait experiment using the IMU sensors was considered to have a significant impact. The measurement error of the FPA acquired from this study differed from the results of previous studies, which used an IMU of approximately 0.22°–0.32°.
Overall, the SL, FC, and FPA values estimated using UWB were similar to those estimated using IMU or displayed marginal differences.
There are two likely sources of error in the methodology with UWB sensors used in this research:
First, in the case of offset removal using SVR, the signal strength information is used to indirectly assess the variation in the relative angle of the anchor and tag attached to different feet. In this case, the relative angle and signal strength of the antenna do not display a linear relationship [40], and the signal strength can vary from the external impact of gadgets, such as Bluetooth devices. Therefore, to address this issue, it is necessary to compensate for the offset more accurately using a continuous gait phase that directly relates the varying antenna angle during walking to the ankle angle.
Second, an error occurred in the hypothesis applied while calculating the foot trajectory using the trilateration method. An anchor with four fixed points is required to determine the foot trajectory using the trilateration method. However, the locations of the two feet in this study were assumed to not intersect in the medial–lateral direction (feet are not entangled), and the number of anchors was reduced to three. In addition, the three locations of the anchor in the stance phase were assumed to be in a fixed coordinate system on the ground before the trajectory of the opposite foot was tracked. However, in reality, the repeated occurrence of a situation wherein one foot is at the heel contact point and the opposite foot is at the heel lift-off point, and the anchor is not entirely fixed to the ground, results in errors in the results of the characteristic analysis.

5. Conclusions

This study estimated the three gait parameters SL, FC, and FPA through gait experiments using shoes installed with four UWB sensors on the right and left feet. The accuracy of the results was similar to that of the studies that applied to the IMU.
The UWB sensor used in this study contains measurement error occurrence factors that originate from compensating data or the setting of the coordinate system for identifying the UWB sensor location. However, it enables direct determination of the SL. Furthermore, although it displays a measurement error similar to that in previous studies using IMU sensors, the drift that occurred with the IMU sensors was absent. Therefore, the applicability of the UWB sensor for gait characteristic analysis was considered high.
In future research, based on the results of this study, the gait phase detection method using an insole device could be applied to develop a more accurate foot trajectory tracking technology and conduct gait disorder diagnosis on actual patients.

Author Contributions

Conceptualization, J.S.P. and C.H.K.; methodology, J.S.P. and C.H.K.; software, J.S.P.; validation, J.S.P. and C.H.K.; formal analysis, J.S.P. and C.H.K.; investigation, J.S.P. and C.H.K.; data curation, J.S.P. and B.L.; writing—original draft preparation, J.S.P.; writing—review and editing, S.P. and C.H.K.; supervision, C.H.K.; project administration, C.H.K.; funding acquisition, C.H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board of the Korea Institute of Science and Technology (approval no. KIST-HR-001 in 2022).

Informed Consent Statement

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

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available owing to ethical concerns (because they were obtained in a clinical trial).

Acknowledgments

This research was funded by the Korea Institute of Science and Technology (KIST) Institutional Program (project nos. 2E31642 and 2E32341).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pirker, W.; Katzenschlager, R. Gait disorders in adults and the elderly: A clinical guide. Wien. Klin. Wochenschr. 2017, 129, 81–95. [Google Scholar] [CrossRef] [PubMed]
  2. Rubino, F.A. Gait disorders. Neurologist 2002, 8, 254–262. [Google Scholar] [CrossRef] [PubMed]
  3. Maki, B.E. Gait changes in older adults: Predictors of falls or indicators of fear? J. Am. Geriatr. Soc. 1997, 45, 313–320. [Google Scholar] [CrossRef]
  4. Ng, K.D.; Mehdizadeh, S.; Iaboni, A.; Mansfield, A.; Flint, A.; Taati, B. Measuring gait variables using computer vision to assess mobility and fall risk in older adults with dementia. IEEE J. Transl. Eng. Health Med. 2020, 8, 2100609. [Google Scholar] [CrossRef] [PubMed]
  5. Barrett, R.S.; Mills, P.M.; Begg, R.K. A systematic review of the effect of ageing and falls history on minimum foot clearance characteristics during level walking. Gait Posture 2010, 32, 429–435. [Google Scholar] [CrossRef] [PubMed]
  6. Simic, M.; Wrigley, T.V.; Hinman, R.S.; Hunt, M.A.; Bennell, K.L. Altering foot progression angle in people with medial knee osteoarthritis: The effects of varying toe-in and toe-out angles are mediated by pain and malalignment. Osteoarthr. Cartil. 2013, 21, 1272–1280. [Google Scholar] [CrossRef] [PubMed]
  7. Rutherford, D.J.; Hubley-Kozey, C.L.; Deluzio, K.J.; Stanish, W.D.; Dunbar, M. Foot progression angle and the knee adduction moment: A cross-sectional investigation in knee osteoarthritis. Osteoarthr. Cartil. 2008, 16, 883–889. [Google Scholar] [CrossRef]
  8. Wang, S.; Mo, S.; Chung, R.C.K.; Shull, P.B.; Ribeiro, D.C.; Cheung, R.T.H. How foot progression angle affects knee adduction moment and angular impulse in patients with and without medial knee osteoarthritis: A meta-analysis. Arthritis Care Res. 2021, 73, 1763–1776. [Google Scholar] [CrossRef]
  9. Jamari, J.; Ammarullah, M.I.; Santoso, G.; Sugiharto, S.; Supriyono, T.; Permana, M.S.; van der Heide, E. Adopted walking condition for computational simulation approach on bearing of hip joint prosthesis: Review over the past 30 years. Heliyon 2022, 8, e12050. [Google Scholar] [CrossRef]
  10. Karatsidis, A.; Richards, R.E.; Konrath, J.M.; van den Noort, J.C.; Schepers, H.M.; Bellusci, G.; Harlaar, J.; Veltink, P.H. Validation of wearable visual feedback for retraining foot progression angle using inertial sensors and an augmented reality headset. J. Neuroeng. Rehabil. 2018, 15, 78. [Google Scholar] [CrossRef]
  11. Begg, R.; Galea, M.P.; James, L.; Sparrow, W.A.T.; Levinger, P.; Khan, F.; Said, C.M. Real-time foot clearance biofeedback to assist gait rehabilitation following stroke: A randomized controlled trial protocol. Trials 2019, 20, 317. [Google Scholar] [CrossRef] [PubMed]
  12. Jacob, S.A. Real-Time Minimum Foot Clearance Estimation Using a Wearable Biofeedback System to Prevent Trip-Related Falls. Master’s Thesis, University of Toronto, Toronto, ON, Canada, 2022. [Google Scholar]
  13. Kaźmierczak, K.; Wareńczak-Pawlicka, A.; Miedzyblocki, M.; Lisiński, P. Effect of treadmill training with visual biofeedback on selected gait parameters in subacute hemiparetic stroke patients. Int. J. Environ. Res. Public Health 2022, 19, 16925. [Google Scholar] [CrossRef] [PubMed]
  14. Sui, J.D.; Chang, T.S. IMU based deep stride length estimation with self-supervised learning. IEEE Sens. J. 2021, 21, 7380–7387. [Google Scholar] [CrossRef]
  15. Singh, P.; Esposito, M.; Barrons, Z.; Clermont, C.A.; Wannop, J.; Stefanyshyn, D. Measuring gait velocity and stride length with an ultrawide bandwidth local positioning system and an inertial measurement unit. Sensors 2021, 21, 2896. [Google Scholar] [CrossRef]
  16. Fusca, M.; Negrini, F.; Perego, P.; Magoni, L.; Molteni, F.; Andreoni, G. Validation of a wearable IMU system for gait analysis: Protocol and application to a new system. Appl. Sci. 2018, 8, 1167. [Google Scholar] [CrossRef]
  17. Wang, L.; Sun, Y.; Li, Q.; Liu, T.; Yi, J. IMU-based gait normalcy index calculation for clinical evaluation of impaired gait. IEEE J. Biomed. Health Inform. 2021, 25, 3–12. [Google Scholar] [CrossRef]
  18. Yang, S.; Laudanski, A.; Li, Q. Inertial sensors in estimating walking speed and inclination: An evaluation of sensor error models. Med. Biol. Eng. Comput. 2012, 50, 383–393. [Google Scholar] [CrossRef]
  19. Han, Y.C.; Wong, K.I.; Murray, I. Gait phase detection for normal and abnormal gaits using IMU. IEEE Sens. J. 2019, 19, 3439–3448. [Google Scholar] [CrossRef]
  20. Perez-Ibarra, J.C.; Siqueira, A.A.; Krebs, H.I. Real-time identification of gait events in impaired subjects using a single-IMU foot-mount device. IEEE Sens. J. 2019, 20, 2616–2624. [Google Scholar] [CrossRef]
  21. 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]
  22. Fan, B.; Li, Q.; Tan, T.; Kang, P.; Shull, P.B. Effects of IMU sensor-to-segment misalignment and orientation error on 3D knee joint angle estimation. IEEE Sens. J. 2021, 22, 2543–2552. [Google Scholar] [CrossRef]
  23. Zjhajehzadeh, S.; Park, E.J. A novel biomechanical model-aided IMU/UWB fusion for magnetometer free lower body motion capture. IEEE Trans. Syst. Man Cybern. 2016, 47, 927–938. [Google Scholar] [CrossRef]
  24. Benoussaad, M.; Sijobert, B.; Mombaur, K.; Coste, C.A. Robust foot clearance estimation based on the integration of foot-mounted IMU acceleration data. Sensors 2015, 16, 12. [Google Scholar] [CrossRef] [PubMed]
  25. Arami, A.; Raymond, N.S.; Aminian, K. An Accurate Wearable Foot Clearance Estimation System: Toward a Real-Time Measurement System. IEEE Sens. J. 2017, 17, 2542–2549. [Google Scholar] [CrossRef]
  26. Delfi, G.; Al Bochi, A.; Dutta, T. A scoping review on minimum foot clearance measurement: Sensing modalities. Int. J. Environ. Res. Public Health 2021, 18, 10848. [Google Scholar] [CrossRef]
  27. Weenk, D.; Roetenberg, D.; van Beijnum, B.J.J.; Hermens, H.J.; Veltink, P.H. Ambulatory estimation of relative foot positions by fusing ultrasound and inertial sensor data. IEEE Trans. Neural. Syst. Rehabil. Eng. 2014, 23, 817–826. [Google Scholar] [CrossRef] [PubMed]
  28. Sabatini, A.M. Quaternion-based strap-down integration method for applications of inertial sensing to gait analysis. Med. Biol. Eng. Comput. 2005, 43, 94–101. [Google Scholar] [CrossRef]
  29. Huang, Y.; Jirattigalachote, W.; Cutkosky, M.R.; Zhu, X.; Shull, P.B. Novel foot progression angle algorithm estimation via foot-worn, magneto-inertial sensing. IEEE Trans. Biomed. Eng. 2016, 63, 2278–2285. [Google Scholar] [CrossRef]
  30. De Vries, W.H.K.; Veeger, H.E.J.; Baten, C.T.M.; Van Der Helm, F.C.T. Magnetic distortion in motion labs, implications for validating inertial magnetic sensors. Gait Posture 2009, 29, 535–541. [Google Scholar] [CrossRef]
  31. Wouda, F.J.; Jaspar, S.L.J.O.; Harlaar, J.; van Beijnum, B.F.; Veltink, P.H. Foot progression angle estimation using a single foot-worn inertial sensor. J. Neuroeng. Rehabil. 2021, 18, 37. [Google Scholar] [CrossRef]
  32. Elsanhoury, M.; Makela, P.; Koljonen, J.; Valisuo, P.; Shamsuzzoha, A.; Mantere, T.; Elmusrati, M.; Kuusniemi, H. Precision Positioning for Smart Logistics Using Ultra-Wideband Technology-Based Indoor Navigation: A Review. IEEE Access 2022, 10, 44413–44445. [Google Scholar] [CrossRef]
  33. Zhang, Y.; Tan, X.; Zhao, C. UWB/INS integrated pedestrian positioning for robust indoor environments. IEEE Sens. J. 2020, 20, 14401–14409. [Google Scholar] [CrossRef]
  34. Gaetano, D.; McEvoy, P.; Ammann, M.J.; Brannigan, C.; Keating, L.; Horgan, F. On-body fidelity factor for footwear antennas over different ground materials. In Proceedings of the 8th European Conference on Antennas and Propagation (EuCAP 2014), The Hague, The Netherlands, 6–11 April 2014; pp. 1410–1414. [Google Scholar] [CrossRef]
  35. MacGougan, G.; O’Keefe, K.; Klukas, R. Ultra-wideband ranging precision and accuracy. Meas. Sci. Technol. 2009, 20, 095105. [Google Scholar] [CrossRef]
  36. Shaban, H.A.; Abou El-Nasr, M.; Buehrer, R.M. Toward a highly accurate ambulatory system for clinical gait analysis via UWB radios. IEEE Trans. Inf. Technol. Biomed. 2010, 14, 284–291. [Google Scholar] [CrossRef]
  37. Qi, Y.; Soh, C.B.; Gunawan, E.; Low, K.S.; Maskooki, A. A novel approach to joint flexion/extension angles measurement based on wearable UWB radios. IEEE J. Biomed. Health Inform. 2014, 18, 300–308. [Google Scholar] [CrossRef] [PubMed]
  38. Anderson, B.; Shi, M.; Tan, V.Y.F.; Wang, Y. Mobile gait analysis using foot-mounted UWB sensors. Proc. ACM Interact. Mob. Wearable Ubiquitous Technol. 2019, 3, 1–22. [Google Scholar] [CrossRef]
  39. Zhang, C.; Bao, X.; Wei, Q.; Ma, Q.; Yang, Y.; Wang, Q. A Kalman filter for UWB positioning in LOS/NLOS scenario. In Proceedings of the 4th International Conference on Ubiquitous Positioning, Indoor Navigation and Location Based Services (UPINLBS), Shanghai, China, 2–4 November 2016. [Google Scholar]
  40. Smaoui, N. Improving Angle and Time Estimation for Concurrent Ultra-Wideband Localization through Transmitter-Side Techniques. Ph.D. Dissertation, University of Houston, Houston, TX, USA, 2021. [Google Scholar]
  41. Su, T.; Ling, H. On modeling mutual coupling in antenna arrays using the coupling matrix. Microw. Opt. Technol. Lett. 2001, 28, 231–237. [Google Scholar] [CrossRef]
  42. Smola, A.J.; Scholkopf, B. Learning with Kernels; GMD-Forschungszentrum Informationstechnik: Berlin, Germany, 1998; Volume 4. [Google Scholar]
  43. Uchitomi, H.; Hirobe, Y.; Miyake, Y. Three-dimensional continuous gait trajectory estimation using single shank-worn inertial measurement units and clinical walk test application. Sci. Rep. 2022, 12, 5368. [Google Scholar] [CrossRef]
  44. Hori, K.; Mao, Y.; Ono, Y.; Ora, H.; Hirobe, Y.; Sawada, H.; Inaba, A.; Orimo, S.; Miyake, Y. Inertial Measurement Unit-Based Estimation of Foot Trajectory for Clinical Gait Analysis. Front. Physiol. 2019, 10, 1530. [Google Scholar] [CrossRef]
  45. Fan, B.; Li, Q.; Liu, T. Accurate foot clearance estimation during level and uneven ground walking using inertial sensors. Meas. Sci. Technol. 2020, 31, 055106. [Google Scholar] [CrossRef]
  46. Huang, C.; Fukushi, K.; Wang, Z.; Nihey, F.; Kajitani, H.; Nakahara, K. An algorithm for real time minimum toe clearance estimation from signal of in-shoe motion sensor. In Proceedings of the 43rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Jalisco, Mexico, 26 July 2021; pp. 6775–6778. [Google Scholar] [CrossRef]
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 humans or property resulting from any concepts, methods, instructions, or products referred to in the content.
Figure 1. Schematic illustration of the UWB sensor system to analyze gait characteristics.
Figure 1. Schematic illustration of the UWB sensor system to analyze gait characteristics.
Applsci 13 04801 g001
Figure 2. Test setup: (a) Data acquisition system; (b) UWB and MCS marker installation; (c) Coordinate system for the UWB sensors.
Figure 2. Test setup: (a) Data acquisition system; (b) UWB and MCS marker installation; (c) Coordinate system for the UWB sensors.
Applsci 13 04801 g002
Figure 3. Gait analysis and its parameters: (a) Flowchart of the gait parameter estimation algorithm; (b) Gait parameters: foot clearance and stride length; (c) Gait parameter: foot progression angle.
Figure 3. Gait analysis and its parameters: (a) Flowchart of the gait parameter estimation algorithm; (b) Gait parameters: foot clearance and stride length; (c) Gait parameter: foot progression angle.
Applsci 13 04801 g003
Figure 4. Schematic of foot progression angle definition.
Figure 4. Schematic of foot progression angle definition.
Applsci 13 04801 g004
Figure 5. Comparison of the UWB data with motion capture system data. Black-dotted line: UWB raw data with the Savitzky–Golay filter, blue-solid line: data obtained using motion capture system, red-solid line: UWB data offset compensated using SVR.
Figure 5. Comparison of the UWB data with motion capture system data. Black-dotted line: UWB raw data with the Savitzky–Golay filter, blue-solid line: data obtained using motion capture system, red-solid line: UWB data offset compensated using SVR.
Applsci 13 04801 g005
Figure 6. Foot trajectory estimation using UWB sensors: (a) All the foot trajectory results; (b) Only the results estimated in the swing phase.
Figure 6. Foot trajectory estimation using UWB sensors: (a) All the foot trajectory results; (b) Only the results estimated in the swing phase.
Applsci 13 04801 g006
Figure 7. Bland–Altman graphs for gait parameter values determined using UWB and motion capture system: (a) Stride length; (b) Foot clearance; (c) Foot progression angle.
Figure 7. Bland–Altman graphs for gait parameter values determined using UWB and motion capture system: (a) Stride length; (b) Foot clearance; (c) Foot progression angle.
Applsci 13 04801 g007
Figure 8. Box plot results for difference between motion capture and UWB approach: (a) Difference in stride length; (b) Difference in foot clearance; (c) Difference in foot progression angle test-results.
Figure 8. Box plot results for difference between motion capture and UWB approach: (a) Difference in stride length; (b) Difference in foot clearance; (c) Difference in foot progression angle test-results.
Applsci 13 04801 g008
Table 1. Subject information.
Table 1. Subject information.
Subject GroupMaleFemale
Number of participants55
Age (years)32.2 ± 3.728.0 ± 4.5
Height (m)1.75 ± 0.041.61 ± 0.06
Weight (kg)71.0 ± 6.149.6 ± 6.0
Table 2. Means and standard deviations of differences in gait parameters between motion capture and UWB approach.
Table 2. Means and standard deviations of differences in gait parameters between motion capture and UWB approach.
SubjectStride Length (mm)Foot Clearance (mm)Foot Progression Angle (°)
2.0 km/h3.0 km/h4.0 km/h2.0 km/h3.0 km/h4.0 km/h2.0 km/h3.0 km/h4.0 km/h
Sb0140.67 ± 24.6134.63 ± 25.8726.48 ± 46.728.46 ± 1.848.17 ± 1.327.89 ± 1.521.34 ± 3.262.10 ± 2.161.50 ± 1.98
Sb0243.09 ± 23.1447.90 ± 24.7851.82 ± 19.818.23 ± 1.717.82 ± 2.118.63 ± 1.461.58 ± 1.871.54 ± 1.012.95 ± 1.33
Sb0346.48 ± 22.7747.84 ± 25.7250.01 ± 20.368.22 ± 1.766.28 ± 1.926.65 ± 1.923.12 ± 1.323.26 ± 1.924.09 ± 0.95
Sb0443.94 ± 25.2647.13 ± 39.0840.92 ± 13.358.91 ± 2.068.41 ± 1.827.24 ± 1.892.42 ± 0.243.31 ± 0.353.24 ± 0.17
Sb0545.17 ± 33.0447.49 ± 26.5746.14 ± 25.286.04 ± 2.007.10 ± 2.776.45 ± 4.882.00 ± 0.322.14 ± 0.412.19 ± 0.31
Sb0646.10 ± 34.5955.94 ± 23.9952.23 ± 19.178.68 ± 1.257.72 ± 0.957.88 ± 0.942.68 ± 1.242.60 ± 0.783.06 ± 0.89
Sb0744.23 ± 21.5038.96 ± 15.7029.77 ± 43.565.63 ± 1.455.74 ± 1.636.28 ± 1.673.22 ± 0.213.14 ± 1.183.31 ± 1.17
Sb0845.58 ± 10.6648.94 ± 12.6152.47 ± 11.218.40 ± 0.828.27 ± 1.308.51 ± 1.354.16 ± 0.434.23 ± 1.424.70 ± 0.33
Sb0946.68 ± 10.7846.39 ± 18.7944.26 ± 20.327.82 ± 2.716.33 ± 3.407.43 ± 1.091.80 ± 1.992.60 ± 1.412.00 ± 1.14
Sb1047.97 ± 19.0552.13 ± 19.8049.87 ± 8.098.07 ± 2.018.78 ± 3.347.71 ± 2.152.49 ± 0.543.17 ± 0.682.78 ± 0.43
Average45.13 ± 21.9946.99 ± 22.3445.40 ± 24.837.80 ± 1.887.45 ± 3.177.54 ± 2.272.53 ± 1.542.86 ± 1.333.12 ± 1.31
45.84 ± 23.107.60 ± 2.502.82 ± 1.42
Table 3. Mean gait parameter difference of motion capture compared with other sensors.
Table 3. Mean gait parameter difference of motion capture compared with other sensors.
Gait ParameterResearcherDifference of Motion Capture
Stride lengthKIST45.84 mm
Uchitomi et al. [43]54 mm
Hori et al. [44]−27 mm
Foot clearanceKIST7.60 mm
Fan et al. [45]4.2 mm (Heel)/13.1 mm (Toe)
Huang et al. [46]7.0 mm/2.7 mm (Toe)
Foot progression angleKIST2.82°
Wouda et al. [31]2.6°
Huang et al. [29]2.5°
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

Park, J.S.; Lee, B.; Park, S.; Kim, C.H. Estimation of Stride Length, Foot Clearance, and Foot Progression Angle Using UWB Sensors. Appl. Sci. 2023, 13, 4801. https://doi.org/10.3390/app13084801

AMA Style

Park JS, Lee B, Park S, Kim CH. Estimation of Stride Length, Foot Clearance, and Foot Progression Angle Using UWB Sensors. Applied Sciences. 2023; 13(8):4801. https://doi.org/10.3390/app13084801

Chicago/Turabian Style

Park, Ji Su, Bohyun Lee, Shinsuk Park, and Choong Hyun Kim. 2023. "Estimation of Stride Length, Foot Clearance, and Foot Progression Angle Using UWB Sensors" Applied Sciences 13, no. 8: 4801. https://doi.org/10.3390/app13084801

APA Style

Park, J. S., Lee, B., Park, S., & Kim, C. H. (2023). Estimation of Stride Length, Foot Clearance, and Foot Progression Angle Using UWB Sensors. Applied Sciences, 13(8), 4801. https://doi.org/10.3390/app13084801

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