Abstract
Analyzing human gait with inertial sensors provides valuable insights into a wide range of health impairments, including many musculoskeletal and neurological diseases. A representative and reliable assessment of gait requires continuous monitoring over long periods and ideally takes place in the subjects’ habitual environment (real-world). An inconsistent sensor wearing position can affect gait characterization and influence clinical study results, thus clinical study protocols are typically highly proscriptive, instructing all participants to wear the sensor in a uniform manner. This restrictive approach improves data quality but reduces overall adherence. In this work, we analyze the impact of altering the sensor wearing position around the waist on sensor signal and step detection. We demonstrate that an asymmetrically worn sensor leads to additional odd-harmonic frequency components in the frequency spectrum. We propose a robust solution for step detection based on autocorrelation to overcome sensor position variation (sensitivity = 0.99, precision = 0.99). The proposed solution reduces the impact of inconsistent sensor positioning on gait characterization in clinical studies, thus providing more flexibility to protocol implementation and more freedom to participants to wear the sensor in the position most comfortable to them. This work is a first step towards truly position-agnostic gait assessment in clinical settings.
Introduction
Various neurological and musculoskeletal diseases result in motor impairments that negatively affect human locomotion [1]. Changes in human gait can be monitored through a variety of parameters including walking speed [2, 3], step length [4, 5], or cadence [6]. Besides the analysis of human gait in special gait laboratories, inertial measurement units can be used to monitor gait in free-living, real-world environments [7]. Small body-mounted sensors represent a promising alternative to specialized gait analysis facilities, enabling unobtrusive long-term monitoring in the subject’s habitual environment. Therefore, these sensors are increasingly used in clinical studies [8, 9], even though the analysis of data derived from such wearable sensors remains challenging [10]. A critical human factor for the reliability of sensor-based gait analysis is the sensor wearing position, as it strongly impacts the recorded signal and the level of comfort during wearing. A sensor position close to the body’s center of mass, for example a belt buckle [11], allows for an easier transformation and interpretation of the inertial measurement unit signal when compared to other positions, such as the wrist, which has more dynamic and unpredictable movements [12]. The wrist is generally more comfortable than the ankle or foot, which has been shown to be most accurate for measuring gait parameters [13]. A good compromise of signal periodicity, stability, and wearing comfort might be a sensor position around the waist [11].
In a real-world environment, the correct sensor attachment is dependent upon the subjects’ commitment, discipline, and functional ability to comply with the provided instructions, which is in turn dependent on the burden such as mounting complexity and the level of comfort to adhere to wearing instructions over a long period and possibly on the type of clothing (e.g., when worn on a belt around the waist). When a subject wears a sensor over multiple weeks, it is likely that they unintentionally change the sensor wearing position from day to day or deliberately change the position to improve the comfort of wearing or for other reasons. Medically relevant changes in a subject’s gait over time could therefore be confounded by systematic as well as short-term frequent changes in wearing position.
A non-central sensor attachment around the waist leads to a different sensor displacement during left and right steps and introduces additional frequencies to the spectrum, in particular odd harmonic frequencies that are also observed during asymmetric gait [14]. How severely the sensor position affects gait analysis is influenced by the underlying algorithm and in particular impacts frequency-based methods [15] that describe gait as a quasiperiodic signal in which successive steps generate a repeating pattern over time [16]. The frequency spectrum changes when the sensor wearing position is altered [17]. In the patient’s real world, it would be desirable to allow flexibility in sensor wearing positions to increase the subject’s compliance and wearing comfort while achieving reliable analysis of the data to avoid confounding of clinical and technical changes.
In this study, we systematically investigated how changing sensor position around the waist affects the recorded signal, and how these changes impact frequency-based step detection which is fundamental for continuous gait monitoring. Based on harmonic frequency analysis we explain the differences in sensor position and explore an autocorrelation algorithm for the reliable estimation of steps for different sensor positions. This work is a first step in sensor position agnostic assessment of gait characteristics, enabling greater flexibility in protocol implementation and collection of more robust results from performance and behavioral gait assessments.
Materials and Methods
Data Acquisition
Twenty healthy volunteers walked a 20-m distance in 3 different walking styles, each with up to 3 repetitions (see Table 1 for study demographics). To analyze the algorithm performance during different walking styles we included straight and turn walks at a self-selected speed and straight walks at low speed. The study had no specific inclusion criteria other than “fit for work” and the first 20 participants were selected.
A marked start and finish line predetermined the path length. The subject started and stopped walking 2 m behind the line to reduce the impact of acceleration and deceleration. The duration of the walk was measured with a stopwatch to ensure a walking speed difference between self-selected and slow walk styles. The tasks were performed indoors on a hard surface and with no specific shoe wear. Each subject was equipped with 11 time-synchronized inertial sensors (Physilog 5, sampling frequency 128 Hz, GaitUp SA, Lausanne, Switzerland): 1 × right shoe, 1 × left shoe, 5 × hip, 1 × right pocket, 1 × collar, 1 × upper arm, 1 × wrist. The sensors were attached by the subject with instructions from the study tester. The SREF sensor was placed with particular care to ensure its correct positioning, while for the other sensors we allowed slight individual changes as expected in a real-world scenario. For this work, we restrict our analysis to sensors worn around the waist (Fig. 1), but data from all positions are available for download [18].
SREF represents the sensor position for which the step detection algorithm was originally designed [19] and to which we compare the tri-axial acceleration data from the other sensors. One subject had to be excluded due to a wrong sensor attachment. In total, 147 unique walks were recorded. No participant stated any health-related issues or gait impairments. Detailed data processing is described in the online supplementary materials (for all online suppl. material, see www.karger.com/doi/10.1159/000511611).
Biomechanics and Harmonic Frequencies during Gait
A non-central sensor wearing position alters the acceleration signal [20, 21] and impacts the harmonic frequency properties of the corresponding spectrum. A sensor attached close to the center of mass (SREF) experiences approximately the same displacement (ds) during a right or a left step in symmetric gait. The non-central sensor attachment around the waist (e.g., S90) leads to a different sensor displacement during left and right steps (Fig. 2a), introducing additional odd harmonic frequencies of the underlying stride frequency (Fig. 2b). Odd harmonic frequencies alternatingly increase and decrease the signal magnitude of right and left steps (see online suppl. materials).
Frequency Analysis Based on Harmonic Theory
The gait signal of each walk was initially filtered and a fast Fourier transformation (FFT) was applied on every axis to decompose the signal into its frequency content. The first 4 harmonic frequency components of the underlying stride frequency were extracted for every axis. The influence of the sensor position on the harmonic frequency components of all axes was analyzed with a principle component analysis.
The harmonic ratio (HR) is a measure of signal symmetry [14]. Here, HR is calculated as the sum of the even harmonics divided by the sum of odd harmonic frequency components per walk and axis.
Algorithms for Gait Analysis
Starting from our previously published FFT-based gait algorithm [19] we developed an autocorrelation-based step estimation with the aim of overcoming problems with frequencies introduced by asymmetric sensor wearing positions and compared performance between the 2 algorithms. The FFT-based approach utilizes the dominant frequency of the 3D signal for step detection. The autocorrelation approach works similarly but uses the fundamental (lowest) frequency of all 3 axes. Detailed descriptions of these algorithms can be found in the online supplementary materials.
Algorithm Performance Assessment
The step detection was evaluated based on the already validated SREF sensor [19]. For this purpose, the time stamps of the detected steps were compared between SREF and the other sensors per walk with a maximum delay of 0.3 s accounting for synchronization imprecision and overall resulting in best performance. An SREF step was matched by identifying the corresponding step with the lowest temporal difference for the other positions. True positive (TP) steps are those matching between the 2 sensors, false positives (FP) are steps in the sensor not matched by SREF steps, and false negatives (FN) are steps in SREF that are not matched by the other sensors. We used a critical success index (CSI; calculated as TP divided by the sum of TP, FP, and FN) as the performance measure.
Results
Impact of Wear Position on Harmonic Frequency Components
Principal component analysis was applied on the first 4 harmonic frequency components per axis and for each walk (shown in Fig. 3). Sensors attached centrally (SREF, S0, S180) show strong overlap, while the centroids of sensors attached on the outside of the hip (S45 and S90) are in a separate group and show higher spread than the central positions. The third harmonic frequency (odd harmonic) of the vertical and the anterior-posterior axes have the highest magnitude and orthogonal direction and explain most of the variance between the signals of the sensor positions. The high spread in S45 and S90 is likely due to larger differences in sensor placement between participants and subject-specific gait characteristics impacting the outer positions more than the central positions.
Figure 4 shows the distribution of the inverse HR for all subjects and sensors for the vertical axis. The HR and thus the asymmetry increases with attachment distance from the center of mass (note that the position on the back at 180° is also near the center of mass). However, the variability of the HR and the overlap in positions is to some extend subject dependent (e.g., subjects 3 and 8). The HR of the anterior-posterior axis shows a similar pattern whereas the medio-lateral axis has no clear pattern and shows a higher variance (see online suppl. figures).
Assessment of Step Detection Accuracy
Figure 5 shows the step detection performance per sensor and walk type for the FFT-based and our autocorrelation-based algorithm. As the sensor is positioned further away from the belt buckle position, the step detection accuracy of the FFT algorithm significantly decreases (green boxes). The FFT assumes a specific axis rotation, in our case a center of mass positioning, that can only be estimated with limited certainty in a real-world study. To demonstrate the impact of axis rotation, we included results obtained from the FFT method on “rotated” data, which has been preprocessed to reorientate the signal to the expected “belt buckle” center of mass orientation. Rotating the data substantially improves the performance of the FFT algorithm at all positions except S180 (blue boxes), but precise sensor data rotation is only possible if the wearing position is a priori known. Therefore, in a real-world setting the determination of the precise rotation is challenging. The autocorrelation method performs highly accurately for all positions without requiring a priori data rotation (red boxes). The performance of both methods is slightly lower during slow walks. The performance of straight walks and turn walks at self-selected speed were almost identical and are therefore pooled in Figure 5 (also see online suppl. Table 1). Overall, the autocorrelation achieves a CSI of 0.99 (across walks) compared to 0.84 for the position-specific FFT and 0.96 for the a priori rotated FFT. The average sensitivity and precision across walks and asymmetric sensor positions were both >0.99 for the autocorrelation and almost 0.99 for the rotated FFT, compared to 0.92 and 0.9, respectively, for the FFT.
Discussion
Our study shows that altering the sensor wearing position systematically affects step detection of frequency-based algorithms. We propose an autocorrelation-based approach, which eliminates the need for sensor data rotation, to determine the underlying gait frequency more reliably. This approach robustly detects steps for different sensor positions around the waist and allows more flexibility in clinical study protocols to monitor gait continuously with increased analyzable wear time in real-world studies.
The analysis of the harmonic gait frequencies has revealed clear differences between the various sensor positions and shows that the dominant signal frequency does not always represent the gait frequency for asymmetric sensor positions. With increasing distance from the center of mass the magnitude of odd harmonic frequencies increases too, explaining most of the asymmetry in the signal. However, the source of asymmetry can originate from a non-central wearing position as well as from gait asymmetry, for example due to pathological gait, such as limping. Therefore, gait asymmetry is confounded with asymmetry in sensor location, and we recommend a central attachment when analyzing gait asymmetries. Note, none of the study participants had any known gait impairment, and this is not representative of an elderly population, for example. In this work we could not analyze asymmetric gait together with sensor displacement [22]. Also, algorithm parameters were optimized based on the walks of a few subjects, and we cannot exclude possible overfitting effects. An additional validation study, possibly including pathological gait, is needed to investigate the generalizability.
We have not yet characterized the data from the additional wearing positions from our study (feet, wrist, pocket, upper trunk). These sensor positions also generate a characteristic harmonic frequency pattern during human gait [23, 24]. It is unclear whether the autocorrelation can robustly delineate steps from these sensor positions.
Our results are based on a study conducted under controlled conditions and were performed in a laboratory setting. Measuring gait in the real world poses several problems, such as unknown context and distinguishing gait periods from other activities to avoid FPs [10, 25]. In a long-term, real-world setting, sensor position can change and may impact data analysis and study conclusions. It is therefore necessary to develop robust analytical methods that cope with flexible low-burden study protocols. The Mobilise-D consortium [28] aims to develop robust digital mobility endpoints collected in the patient’s real world over the upcoming years.
Acknowledgements
We thank our colleagues Sarah Hemsley, Ruth Frommherz, and Roland Huber for their help in setting up the study, and Jason Laramie, Anna Kostikova, and Scott Kennedy for their support. We would also like to thank all study participants.
Statement of Ethics
Data were collected as part of a single site, non-interventional study carried out at the Novartis Institutes for Biomedical Research (NIBR) with Novartis employees. The Ethics Committee involved in the study granted approval and written consent was given by every participant.
Conflict of Interest Statement
Ar.M., A.K., F.K., M.U., and B.M.E. are members of the IMI Mobilise-D consortium. I.C. is an employee of and holds stock options in Evidation Health. He has received payment for lecturing on Digital Health at the ETH Zurich and FHNW Muttenz. He is an Editorial Board Member at Karger Digital Biomarkers and a founding member of the Digital Medicine Society. H.H., Am.M., and Ar.M are employees of Novartis and hold stock in Novartis. M.T. is a Novartis employee and a student at University Erlangen-Nürnberg. B.M.E. reports grants outside the submitted work from Adidas AG, Agaplesion gAG, and Bosch Sensortec GmbH. He received compensation from lecturing for AbbVie Deutschland GmbH & Co. KG and Agaplesion gAG. He is a co-founder and co-owner of Portabiles GmbH and Portabiles HealthCare Technologies GmbH, and co-inventor of gait analysis patent application EP 16174268.9. F.K., M.U., and A.K. received funding from IMI Mobilise-D (grant agreement 820820).
Funding Sources
This work was sponsored by Novartis Pharma AG and received no direct funding.
Author Contributions
M.T., Am.M, H.H., F.K., and Ar.M. contributed to the design of data collection and performed analysis and interpretation. I.C. ran the data acquisition. M.U., F.K., and A.K. contributed ideas to the autocorrelation. Ar.M., M.T., and I.C. drafted the manuscript. H.H., F.K., A.K., M.U., Am.M., and B.M.E. revised the manuscript.