Skip to main content

Indirect measurement of anterior-posterior ground reaction forces using a minimal set of wearable inertial sensors: from healthy to hemiparetic walking

Abstract

Background

The anterior-posterior ground reaction force (AP-GRF) and propulsion and braking point metrics derived from the AP-GRF time series are indicators of locomotor function across healthy and neurological diagnostic groups. In this paper, we describe the use of a minimal set of wearable inertial measurement units (IMUs) to indirectly measure the AP-GRFs generated during healthy and hemiparetic walking.

Methods

Ten healthy individuals and five individuals with chronic post-stroke hemiparesis completed a 6-minute walk test over a walking track instrumented with six forceplates while wearing three IMUs securely attached to the pelvis, thigh, and shank. Subject-specific models driven by IMU-measured thigh and shank angles and an estimate of body acceleration provided by the pelvis IMU were used to generate indirect estimates of the AP-GRF time series. Propulsion and braking point metrics (i.e., peaks, peak timings, and impulses) were extracted from the IMU-generated time series. Peaks and impulses were expressed as % bodyweight (%bw) and peak timing was expressed as % stance phase (%sp). A 75%-25% split of 6-minute walk test data was used to train and validate the models. Indirect estimates of the AP-GRF time series and point metrics were compared to direct measurements made by the forceplates.

Results

Indirect measurements of the AP-GRF time series approximated the direct measurements made by forceplates, with low error and high consistency in both the healthy (RMSE= 4.5%bw; R2= 0.93) and post-stroke (RMSE= 2.64%bw; R2= 0.90) cohorts. In the healthy cohort, the average errors between indirect and direct measurements of the peak propulsion magnitude, peak propulsion timing, and propulsion impulse point estimates were 2.37%bw, 0.67%sp, and 0.43%bw. In the post-stroke cohort, the average errors for these point estimates were 1.07%bw, 1.27%sp, and 0.31%bw. Average errors for the braking estimates were higher, but comparable.

Conclusions

Accurate estimates of AP-GRF metrics can be generated using three strategically mounted IMUs and subject-specific calibrations. This study advances the development of point-of-care diagnostic systems that can catalyze the routine assessment and management of propulsion and braking locomotor deficits during rehabilitation.

The neuromechanical processes underlying healthy bipedal locomotion are multi-factorial [13] and converge on locomotor patterns that are characteristically fast, efficient, and stable [1, 4]. An impaired ability to transition from step to step is a locomotor deficit common across diagnostic groups [513]. During the step-to-step transition of each gait cycle, a braking force is generated by the leading limb as it makes contact with the ground in front of the body. To efficiently accelerate the body into the next step, coordination of the timing and magnitude of the forward propulsion force generated by the trailing limb is required [1, 1416]. Moreover, to walk faster, healthy individuals symmetrically increase the magnitude of propulsion generated by each limb while maintaining the relative timing of the propulsion peak [15, 17, 18]. In individuals with impaired propulsion function, walking is often slow, metabolically expensive, and unstable [1922].

Laboratory equipment such as instrumented treadmills and forceplates are the gold standard in characterizing propulsion and braking function during healthy [23, 24] and impaired [5, 6, 9, 10, 20, 2527] walking by way of direct measurements of the anterior-posterior ground reaction forces (AP-GRFs) generated during walking and point metrics extracted from the AP-GRF time series (Fig. 1). For example, older adults are reported to generate up to 22% less peak propulsion (i.e., the peak of the anterior ground reaction force) compared to young adults [23, 24], and in people post-stroke, the propulsion generated by the paretic limb is up to 68% less than the non-paretic limb [9, 20, 26, 27]. Studies that have combined AP-GRF measurements with clinical evaluations have shown the clinical consequences of impaired propulsion function. Indeed, asymmetry in the propulsion impulses generated by the paretic and non-paretic limbs is correlated with hemiparetic severity [9, 28]. Moreover, deficits in propulsion function are highly related to walking speed [29] and long distance walking [30] after stroke—key determinants of community participation and perceived quality of life [19, 31, 32].

Fig. 1
figure 1

Anterior-posterior ground reaction force (AP-GRF) time series and salient propulsion and braking metrics

Despite the importance of propulsion to a functional bipedal gait, conventional rehabilitation efforts have, by and large, been unable to restore propulsion function after neurological injury or dysfunction. The development and study of interventions that target propulsion function is a highly active area of research [12, 3341]; however, the clinical translation of these experimental approaches is hindered by the limited access that rehabilitation clinicians have to the sophisticated instrumentation (i.e., forceplates and instrumented treadmills) and personnel with advanced training required to collect, analyze, and interpret ground reaction force data. Moreover, even in settings with access to a motion analysis laboratory, locomotor differences inherent to treadmill walking and the small collection footprint of most overground forceplate walkways limit ecological validity. Together, these limitations of the current state-of-the-art motivate the development of point-of-care propulsion diagnostic systems. The clinical management of locomotor propulsion deficits will remain untenable if the measurement instruments used to assess limb propulsion remain inaccessible.

Wearable sensors are a promising solution for this measurement gap. Indeed, wearable sensors have been used to extend gait measurements outside of the laboratory [4247] and a wide range of methods and sensors have proven effective in providing indirect measurements of the ground reaction forces generated during walking [4851]. These methods, however, have largely not been effective for the AP-GRFs and depend on assumptions of healthy, consistent walking patterns that may not translate to impaired locomotor patterns [51, 52]. Recent work has shown that inertial measurement units (IMUs) can be used to make measurements during healthy [44] and hemiparetic walking [53] that are highly correlated to key features of propulsion. The aims of this study were to extend this work by describing the use of a minimal set of IMUs to indirectly measure the AP-GRF generated during healthy and hemiparetic walking and provide estimates of: (i) the AP-GRF time series and (ii) salient propulsion and braking point metrics (i.e., peak magnitudes, peak timings, and impulses) extracted from the time series (see Fig. 1).

Methods

Participants

Ten healthy individuals (26±4 years, 171±10.8 cm, 68±17 kg) (Table 1, top) who were free of conditions that impaired walking ability (as per self-report and confirmation during the study visit) and five individuals with chronic post-stroke hemiparesis (58±15 years, 180±2.9 cm, 91±12 kg) (Table 1, bottom) were recruited to participate in this study. The inclusion criteria for study participants who were post-stroke consisted of being greater than six months post-stroke, having the ability to walk without the assistance of another individual, and presenting with observable gait deficits. Exclusion criteria included comorbidities other than stroke that impaired walking ability, resting heart rate outside the range of 40 to 100 beats per minute, resting blood pressure outside the range of 90/60 to 170/90 mmHg, inability to communicate with investigators, and pain in the lower limbs or spine. Individuals post-stroke were recruited from a research participant registry generated from the clinical programs at Boston University, referrals from local clinics and hospitals, and flyers distributed in and around Boston. All study procedures were approved by the Boston University Institutional Review Board and written informed consent was secured from all study participants.

Table 1 Study participant characteristics

Gait evaluation and data collection overview

Study participants completed a testing session that included a standing static trial, 10-meter walk test at a comfortable walking speed (CWS), and a 6-minute walk test with the instruction to cover as much distance as safely possible [54]. The standing static trial served as a reference for the orientation of the IMUs, the 10-meter walk test was used to measure usual walking speed and paretic propulsion symmetry (as described in [9]) to characterize study participants’ baseline walking function (see Table 1), and the 6-minute walk test provided the model training and validation datasets. All walking tests were performed around a 26.6m oval indoor track consisting of two 10m straightaways separated by 3.3m turns on either end. One of the 10m straightaways was instrumented with six forceplates (Bertec, Columbus, OH, USA) located level with the surrounding floor (Fig. 2a) to enable the collection of ground reaction forces (i.e., the reference standard) during the 6-minute walk test.

Fig. 2
figure 2

a IMU setup on a participant and walking track with forceplates mounted level with the surrounding floor. b IMU-measured pelvis acceleration and shank and thigh angles were used in combination with subject-specific estimation algorithms to generate indirect measurements of the anterior-posterior ground reaction force (AP-GRF). Moreover, the shank angle was used to identify the stance phase, defined as the period between heel strike and toe off of the target limb. c IMU-based AP-GRF (indirect) and forceplate measured AP-GRF (direct). Direct AP-GRF measurements were limited to forceplate strikes, thus one stride

Prior to testing, wireless inertial measurement units (IMUs, MTw Awinda, Xsens, Enschede, Netherlands) were securely attached to the posterior pelvis and laterally on the thigh and shank using fabric wraps (Fig. 2a). These three IMUs served as the primary sensor set evaluated in this study, with each IMU selected based on a biomechanics-based model of the propulsion and braking forces generated during walking. More specifically, when considered together, the shank and thigh sagittal angles measured by the IMUs attached to these segments provide information on the limb’s position relative to the body [44] and the acceleration measured by the pelvis IMU serves as a proxy for the body acceleration. Both the limb’s position relative to the body and the body’s acceleration during walking are highly correlated with propulsion and braking during walking [5557]. Each IMU was mounted such that one axis moves along the sagittal plane. Specifically, the pelvis IMU was placed below the fifth lumbar vertebrae on the sacrum and the thigh and shank IMUs were placed laterally, a third of the corresponding segment length away from the knee and ankle joints respectively. For the purposes of this study, two additional IMUs were attached to the contralateral thigh and shank to enable indirect measurement of each limb’s AP-GRF concurrently. The calibration routines and algorithm were implemented in the same way for both the healthy and post-stroke subgroups, with indirect measurement of each limb’s AP-GRF estimated using measurements of global acceleration from the pelvis IMU and segment angles from the respective limb’s thigh and shank segment IMUs. That is, using this approach, only three IMUs (i.e., pelvis, shank, and thigh) are necessary for AP-GRF estimates for one limb, whereas five IMUs (1 pelvis, 2x shank, and 2x thigh) are necessary for AP-GRF estimates from both limbs.

The IMU measurements of pelvis global acceleration and thigh and shank angles (Fig. 2b) were used to generate subject-specific estimates of the AP-GRF time series (Fig. 2c), from which key metrics were extracted (see “Analyses” section). Both IMU and forceplate signals were temporally synchronized using a synchronization pulse triggered at the start of data collection in both data collection systems [58]. Whereas our IMU approach provided indirect measurements of the AP-GRF for all strides (Fig. 2c), only strides with full direct measurement of the individual limb AP-GRF by the forceplates (i.e., no partial strikes) were used for comparative analyses. Because the six forceplates were located within the middle 4m of one of the 10m straightaways of the walking track, steps taken during turns were not included in these analyses. For eight of the healthy study participants, data from the right limb are reported. For the remaining two healthy study participants, physical drift of one of the IMUs on the body during the experiment (due to an insecure attachment) made the right limb data unusable. Thus, for these two individuals, data from the left limb are reported. For the post-stroke participants, data from both the paretic and non-paretic limbs are reported. The model training datasets consisted of an average 14±3 strides per healthy participant and 9±4 strides per post-stroke participant (Supplementary Table 1).

Data processing

IMU data were collected at 100Hz, with acceleration data filtered at 10Hz using a second order Butterworth filter. Forceplate data were collected at 2000Hz, filtered at 10 Hz using a second order Butterworth filter, and down sampled to match the IMU collection frequency. All gait data were segmented between initial contact and toe-off and time-normalized to 100 points to represent the stance phase of walking. When forceplate data were available, initial contact and toe-off events were defined using the vertical ground reaction force using a 30 N threshold. When forceplate data were not available, the maximum and minimum peaks in the IMU-measured shank angle were used, with initial contact defined by the maximum peak and toe-off defined by the minimum peak [43] (Fig. 2b).

Before modeling, the forceplate data were normalized to the body weight of the subject. The IMU orientations during the static standing trial served as the zero reference for the IMU orientation signals during walking [59, 60]. To compute the segment Euler angles, the quaternion vector relative to the static standing reference was found and subsequently rotated such that the roll axis of the IMU was aligned with the sagittal rotation axis of the segment [6163]. The thigh and shank angles were offset to start at zero degrees at each periodic heel strike to address any drift over time [59, 64]. We also bounded the thigh and shank angles to a sine function, as is common with inverse kinematics equations [65]. Together, the filtered pelvis acceleration signal and the sine of the thigh and shank angles form the basis for the modeling approach described.

Analyses

We analyzed the IMU data with two main goals. The first goal was to identify the IMU sensor set that provides the most accurate and reliable model-informed indirect measurement of the anterior-posterior ground reaction force (AP-GRF) time series. The second goal was to use this indirect measurement of the AP-GRF time series to estimate salient point metrics—i.e., the peaks, peak timings, and impulses of the anterior (propulsion) and posterior (braking) ground reaction force (see Fig. 1). A 75%-25% data split was used for model training and validation, respectively.

Indirect measurement of the anterior-posterior ground reaction force time series

Indirect measurements of the AP-GRF time series Fap,est were made by generating a subject-specific linear regression model using the IMU measurements of pelvis acceleration apelvis, the sine of the thigh angle in the sagittal plane sin(θthigh), the sine of the shank angle in the sagittal plane sin(θshank), and the interactions sin(θshankapelvis, sin(θthighapelvis, and sin(θthigh)× sin(θshank). Thus, we describe Fap,est as

$$ F_{a-p,est} = c \cdot x $$
(1)

where c is a vector [c1...c7] of subject-specific regression coefficients generated from the training dataset and x is a vector [apelvis, sin(θthigh), sin(θshank), sin(θshankapelvis, sin(θthighapelvis, sin(θthigh)× sin(θshank)] of IMU measured components. To evaluate the relative importance of each IMU, we recomputed Fap,est without each of the pelvis, thigh, and shank IMUs—i.e., by setting their respective components in Eq. 1 to zero. The overall fit (R2) and the root mean square error (RMSE) when comparing Fap,est from each model to the AP-GRF time series directly measured by the forceplates (i.e., Fap,act) is presented in Fig. 3.

Fig. 3
figure 3

Comparison of direct (i.e., forceplate-measured) and indirect (i.e., IMU-based) measurements of the anterior-posterior ground reaction force (AP-GRF). a The AP-GRF time series reconstruction enabled by the primary sensor set (i.e., pelvis, thigh, and shank IMUs) and with each of these IMUs removed is shown for exemplar participants from the healthy and post-stroke cohorts. b Root mean square error (RMSE) and consistency (R2) metrics for each sensor set are shown for each study subject

Indirect measurement of propulsion and braking point metrics

Estimates of the propulsion and braking peak magnitudes were defined initially as the maximum and minimum values of the Fap,est time series. If more than one peak was observed, the second peak was used. To improve the accuracy of these point estimates, we introduced an additional subject-specific model that leveraged highly consistent estimates of the peak propulsion and braking timings, which were recorded as a function of the stance phase (%sp) and described as Fpktim,est. More specifically, the values of apelvis, sin(θthigh), and sin(θshank) at Fpktim,est were identified from the Fap,est time series and used to create a new vector \(y = [1, a_{pelvis\_pk-tim}, \sin (\theta _{thigh\_pk-tim}), \sin (\theta _{shank\_pk-tim})]\), which was, in turn, used to describe Fpkmag,est as

$$ F_{pk-mag,est} = d \cdot y $$
(2)

where d is a vector [ d1... d4] consisting of updated regression coefficients. The interaction terms sin(θshankapelvis, sin(θthighapelvis, and sin(θthigh)× sin(θshank) were not included in Eq. 2 as they did not improve the model’s performance.

In preliminary work, we found that Fpktim,est consistently overestimated or underestimated the actual force-plate measured peak timings and was thus revised as Fpktim,revised and described as

$$ F_{revised} = e_{1}+ e_{2} \cdot F_{est} $$
(3)

The final point metrics of interest were the propulsion and braking impulses (Fimp,est). These were computed by summing all of the positive (i.e., for propulsion) and negative (i.e., for braking) values in each Fap,est cycle and dividing by the total number of points (i.e., 100) to yield the propulsion and braking impulses (%bw) per stride. Similar to the peak propulsion and braking timings, we found that Fimp,est consistently overestimated or underestimated the actual propulsion and braking impulses and was thus revised using Eq. 3.

Statistical analysis

All analyses were performed using custom MATLAB scripts (MATLAB, MathWorks, Natick, MA, USA). The regression model was fit using the fitlm() function in MATLAB, with R2 and root mean square error (RMSE) exported directly as a function output. Together, R2 and RMSE show the consistency and accuracy between the Fap,act and Fap,est time series.

For the point metrics of interest, in addition to computing the RMSE, the degree of absolute agreement among direct and indirect measurements was evaluated using two-way mixed effect, absolute agreement, single rater intraclass correlation coefficients (ICCs) [66, 68], with an alpha value of 0.05. ICC values were interpreted using the guidelines provided in [66] with an ICC above 0.90 considered to be excellent, 0.90 to 0.75 as good, 0.75 to 0.50 as moderate, and less than 0.50 as poor.

RMSEs and ICCs between direct and indirect measurements were computed for the training and validation datasets across all strides available from the healthy and post-stroke cohorts. Datapoints were identified as outliers if they were greater than three standard deviations from the mean. If an outlier was identified, the analysis was re-run with the datapoint removed. The results with and without the outlier are reported, and the outlier is shown in the plots.

Results

Measurement of the AP-GRF time series

The indirect measurement of the AP-GRF time series (Fap,est) strongly approximated the direct measurement of the time series made by the reference standard forceplates (Fap,act) in both the healthy (R2= 0.93, RMSE = 4.62%bw) and post-stroke (R2= 0.90, RMSE = 2.64%bw) cohorts. To assess the importance of each IMU used in the equation (Eq. 1), we recomputed R2 and RMSE with the removal of each IMU’s parameters (Fig. 3, Supplementary Table 1).

Removal of the thigh IMU parameters from Eq. 1 resulted in a substantial weakening of the AP-GRF reconstruction in both the healthy (ΔR2 = -0.40, % increase in RMSE = 155%) and post-stroke (ΔR2 = -0.34, % increase in RMSE = 125%) cohorts. Without the shank IMU parameters, the AP-GRF reconstruction was similarly weakened, but to a lesser degree, in both the healthy (ΔR2 = -0.13, % increase in RMSE = 66%) and post-stroke (ΔR2 = -0.13, % increase in RMSE = 55%) cohorts. Removal of the pelvis IMU parameters similarly weakened the AP-GRF reconstruction in both the healthy (ΔR2 = -0.06, % increase in RMSE = 31%) and post-stroke (ΔR2 = -0.07, % increase in RMSE = 33%) cohorts. Each IMU’s importance to the AP-GRF reconstruction varied across study participants (Fig. 3b, Supplementary Table 1).

Measurement of AP-GRF point metrics

Peak propulsion magnitude

In the training dataset, the magnitude of the propulsion peaks in the Fap,est time series strongly approximated the Fap,act time series in the healthy cohort (RMSE = 1.39%bw, ICC = 0.96) and for both the paretic (RMSE = 0.68%bw, ICC = 0.99) and non-paretic (RMSE = 0.97%bw, ICC = 0.99) limbs of the post-stroke cohort (Fig. 4a, training). In the validation dataset, the estimated propulsion peak magnitudes remained a strong approximation of the forceplate-measured propulsion peak magnitudes for the healthy cohort (RMSE = 2.37%bw, ICC = 0.87) and for the paretic (RMSE = 1.07%bw, ICC = 0.98) and non-paretic (RMSE = 1.10%bw, ICC = 0.99) limbs of the post-stroke cohort (Fig. 4a, validation).

Fig. 4
figure 4

Accuracy of propulsion point metrics. Comparison of direct (i.e., forceplate) and indirect (i.e., IMU) measurements of the a peak propulsion magnitude, b peak propulsion timing, and c propulsion impulse. RMSE and ICC means are reported for each training and validation analysis. A single statistical outlier—identified with a red outline—was identified for the post-stroke peak propulsion validation analysis. The findings with and without the outlier are presented, respectively, in parentheses and with an *. Note: high homogeneity in (b, Left), with multiple points overlaid. See Supplementary Table 2 for ICC 95% confidence interval

Peak propulsion timing

In the healthy cohort, the propulsion peak timings in the Fap,est time series strongly approximated the Fap,act time series for both the training (RMSE = 0.63%sp) and validation (RMSE = 0.67%sp) datasets (Fig. 4b). In the post-stroke cohort, similarly strong approximations were observed for both the paretic and non-paretic limb measurements in the training (paretic RMSE = 1.03%sp, non-paretic RMSE = 0.70%sp) and validation (paretic RMSE = 1.27%sp, non-paretic RMSE = 0.71%sp) datasets. High reliability was observed in the post-stroke training and validation datasets (ICCs >0.97); however, in the healthy datasets, highly homogeneous propulsion peak timings (i.e., between 84 to 89%sp) contributed to low reliability (ICCs <0.57), despite high agreement (see [66, 67]).

Propulsion impulse

In the training dataset, the propulsion impulses in the Fap,est time series strongly approximated those of the Fap,act time series in the healthy cohort (RMSE = 0.39%bw, ICC = 0.89) and for both the paretic (RMSE = 0.26%bw, ICC = 0.98) and non-paretic (RMSE = 0.31%bw, ICC = 0.99) limbs of the post-stroke cohort (Fig. 4c, training). In the validation dataset, the estimated propulsion impulses remained a strong approximation of the forceplate-measured propulsion impulses for the healthy cohort (RMSE = 0.43%bw, ICC = 0.86) and for the paretic (RMSE = 0.31%bw, ICC = 0.98) and non-paretic (RMSE = 0.33%bw, ICC = 0.98) limbs of the post-stroke cohort (Fig. 4c, validation).

Peak braking magnitude

In the training dataset, the magnitude of the braking peaks in the Fap,est time series strongly approximated the Fap,act time series in the healthy cohort (RMSE = 1.81%bw, ICC = 0.94) and for both the paretic (RMSE = 1.52%bw, ICC = 0.99) and non-paretic (RMSE = 1.09%bw, ICC = 0.95) limbs of the post-stroke cohort (Fig. 5a, training). In the validation dataset, the estimated braking peak magnitudes remained a strong approximation of the forceplate-measured braking peak magnitudes for the healthy cohort (RMSE = 2.72%bw, ICC = 0.88) and for the paretic (RMSE = 2.46%bw, ICC = 0.97) limbs of the post-stroke cohort (Fig. 5a, validation). The estimated braking peak magnitudes for the non-paretic limbs of the post-stroke cohort had a moderate approximation (RMSE = 2.98%bw, ICC = 0.64).

Fig. 5
figure 5

Accuracy of braking point metrics. Comparison of direct (i.e., forceplate) and indirect (i.e., IMU) measurements of the a peak braking magnitude, b peak braking timing, and c braking impulse. RMSE and ICC means are reported for each training and validation analysis. Note: high homogeneity in (b, Left), with multiple points overlaid. See Supplementary Table 2 for ICC 95% confidence interval

Peak braking timing

In the healthy cohort, the braking peak timings in the Fap,est time series strongly approximated the Fap,act time series for both the training (RMSE = 0.76%sp) and validation (RMSE = 1.37%sp) datasets (Fig. 5b). In the post-stroke cohort, similarly strong approximations were observed for both the paretic and non-paretic limb measurements made in the training (paretic RMSE = 0.73%sp, non-paretic RMSE = 0.78%sp) and validation (paretic RMSE = 1.21%sp, non-paretic RMSE = 1.52%sp) datasets. High reliability was observed in the post-stroke training and validation datasets (ICCs >0.96); however, in the healthy datasets, high homogeneity in the braking peak timings (i.e., between 16 to 22%sp) contributed to lower reliability (ICCs < 0.76), despite high agreement (see [66, 67]).

Braking impulse

In the training dataset, the braking impulses in the Fap,est time series strongly approximated those of the Fap,act time series in the healthy cohort (RMSE = 0.52%bw, ICC = 0.80) and for both the paretic (RMSE = 0.47%bw, ICC = 0.95) and non-paretic (RMSE = 0.40%bw, ICC = 0.95) limbs of the post-stroke cohort (Fig. 5c, training). In the validation dataset, the estimated braking impulses remained a strong approximation of the forceplate-measured braking impulses for the healthy cohort (RMSE = 0.68%bw, ICC = 0.66) and for the paretic (RMSE = 0.56%bw, ICC = 0.93) and non-paretic (RMSE = 0.55%bw, ICC = 0.90) limbs of the post-stroke cohort (Fig. 5c, validation).

Discussion

Laboratory-based instrumented treadmills and forceplates are the gold standard in the direct measurement of the AP-GRF generated during walking [65]. We present the use of a minimal set of wearable inertial sensors to provide, with a high level of accuracy and robustness, indirect measurements of the AP-GRF time series and salient propulsion and braking point metrics. This work advances the development of point-of-care measurement systems that can catalyze the routine assessment and management of propulsion and braking locomotor deficits. Indeed, the approach presented in this paper has near-term potential to overcome logistical hinderances to including overground AP-GRF assessments in laboratory-based movement research and has long-term potential to extend the measurement capabilities of instrumented treadmills and forceplates to clinical and real-world settings.

Despite reduced ecological validity, instrumented treadmills are often utilized instead of overground forceplates in laboratory-based movement research [41, 6973]. Indeed, instrumented treadmills facilitate the collection of a large number of consecutive samples with ease, whereas forceplate walkways have traditionally been limited to a small collection footprint of only several meters and require multiple trials to collect a small number of samples [65, 69]. Our sensor-based approach has the potential to extend the confined and limited length of forceplate walkways by providing indirect AP-GRF measurements for steps where a forceplate is not available, to replace partial forceplate strikes that are not usable, and to minimize threats to the validity of AP-GRF measurements that may arise when test subjects target the forceplates while walking (i.e., by only including steps without a forceplate in the analyses).

Toward a minimal sensor set for AP-GRF time series measurements

When compared to forceplate measurements, three IMUs strategically mounted to the pelvis, thigh, and shank were effective in producing highly accurate indirect measurements of the AP-GRF generated during overground walking by both healthy and post-stroke individuals (Fig. 3). The selection of this sensor set was based on a biomechanical framework linking the AP-GRF generated during walking with certain kinematic features that are readily measured by IMUs. That is, the angles measured by the thigh and shank IMUs were used as a proxy for the orientation of the limb relative to the body and the acceleration measured by the pelvis IMU was used as proxy for the body’s center of mass acceleration. These variables are highly related to propulsion and braking function [44, 5557] and our findings demonstrate the efficacy of this approach for reconstructing the AP-GRF time series.

We found that all three IMUs were required to produce the best estimates of the AP-GRF time series. Removing even one IMU had a substantial impact on estimation accuracy. Removal of the thigh IMU resulted in the greatest change across subjects: a 155% increase in the root mean square error and a 40% reduction in the model strength (R2) for healthy study participants, with a similar effect in the post-stroke cohort. Removing the pelvis or shank IMUs resulted in overall less reductions in accuracy; however, the effects varied across individual subjects (Fig. 3), indicating that including all three IMUs was most advantageous.

In recent work, a model-based approach to estimating the AP-GRF time series and other kinematic and kinetic variables during healthy walking used a 17-IMU sensor set to produce AP-GRF estimates with a root mean square error of 5.5%bw [48]. Our approach with only three IMUs resulted in lower error in healthy study participants (i.e., 4.6%bw) and was also shown to be highly effective for post-stroke hemiparetic gait, with error of 2.64%bw. Taken together with the availability and low cost of IMUs, these findings motivate future applications in both laboratory-based movement research and clinical practice—i.e., to complement and extend overground forceplate data collections in a research lab setting by providing indirect estimates of the AP-GRF time series every step, and providing clinicians with the access to these measurements that they need to inform locomotor diagnoses and treatments.

Accuracy of AP-GRF point metric measurements

Point metrics are often extracted from the AP-GRF time series to characterize deficits in locomotor function [55, 74, 75]. The minimal detectable change (MDC) for many of these point metrics has been computed for hemiparetic walking overground [75]. Comparing the magnitude of error in our estimates to the available overground MDCs is another approach to evaluate measurement accuracy. For the peak paretic propulsion estimates, the root mean square error of 1.07%bw that we observed was lower than the 1.80%bw MDC previously reported [75]. Similarly, for the propulsion impulse estimates, the root mean square error of 0.31%bw for the paretic propulsion impulse was lower than the 0.90%bw MDC previously reported [75]. For the non-paretic limb, the propulsion peak and impulse point estimates had similar error magnitudes as the paretic limb, but MDC references were not available for comparison. The paretic limb’s braking point metrics were observed to have slightly higher error than the propulsion point metrics; however, MDC references were also not available for comparison. Ultimately, the accuracy of our approach is highlighted by the excellent agreement (ICCs >0.90) between IMU and forceplate measurements for most of the point metrics of interest.

In addition to the peak and impulse point metrics, we were also interested in the accuracy of the timing of the propulsion and braking peaks. Across subjects, the root mean square error in the timing of the propulsion and braking peaks was less than 1.4%sp. Despite this low error, the ICC values for the propulsion and braking peak timings were poor to moderate in the healthy cohort (ICCs ranged from 0.33 to 0.57), whereas they were good to excellent in the post-stroke cohort (ICCs ranged from 0.79 to 0.97). The low ICCs for the peak timings observed in the healthy cohort is likely the result of reduced variability in the dataset [66, 67]. In the healthy cohort, peak propulsion timings ranged between 84 and 89%sp and peak braking timings ranged between 16 and 22%sp, whereas substantially more variability was observed for both point metrics in the post-stroke cohort (Figs. 4 and 5b, Supplementary Table 1).

Future considerations

An advantage of modeling the AP-GRF using a biomechanics-based equation is the potential to identify the relative contribution of each term, and how this relative contribution may change over time or as the result of intervention. Each of the terms included in the model reflects biomechanical processes related to the generation of the AP-GRF (i.e., segment/limb orientations and body acceleration) and examining changes in the relative contribution of each term—e.g., changes in standardized coefficients—would presumably reflect changes in these biomechanical processes. In contrast, model-free approaches centered on machine learning may provide more accurate estimates of the AP-GRF, but do not easily allow for examination of underlying biomechanical processes [49, 52]. Hybrid approaches that combine model and model-free terms are worth investigating. Furthermore, the extension of this work to other diagnostic groups with neuromotor impairments that alter the generation of anterior-posterior ground reaction forces (e.g., Multiple Sclerosis [5], Parkinson’s disease [25, 76], spinal cord injury [10], traumatic brain injury [11, 12]) would advance targeted rehabilitation approaches for these populations.

Limitations

There are inherent limitations to using IMU technology for human movement analysis. For example, in addition to drift in the component signals that must be accounted for, IMUs need to be securely mounted to the body segments to avoid physical drift or movement of the sensor relative to the body segment and minimize soft tissue artifacts. As described in the “Methods” section, data from the right limbs of two healthy individuals were not usable due to an insecure attachment of the sensor to the segment. Relatedly, a current requirement of our IMU-based approach is mounting the lower limb IMUs such that one of the IMU planes correspond with the sagittal plane of the segment. Failure to do so would result in erroneous measurements of the sagittal plane angles used in the equations. Future efforts that leverage techniques such as coordinate rotations [60, 62, 63, 77] may enable extension of this approach. Fundamentally, our approach depends on body and limb kinematics to produce a kinetic measurement. Thus, patient subsets and walking conditions characterized by higher locomotor variability or AP-GRFs outside of the range observed in this study (e.g., acute post-stoke patients) may not be suitable for this approach in its current form.

Conclusions

We show that indirect measurements of the anterior-posterior ground reaction forces generated during healthy and post-stroke walking can be made using three strategically-mounted inertial sensors in combination with subject-specific calibrations. This work has near-term potential to overcome logistical hinderances to including overground anterior-posterior ground reaction force assessments during laboratory-based movement research, and has long-term potential to catalyze the routine assessment and management of propulsion and braking deficits during locomotor rehabilitation. This foundational study is a step towards a fully wearable and autonomous point-of-care measurement system that can extend the measurement capabilities of laboratory-based instrumented treadmills and forceplates to ecological settings.

Availability of data and materials

Data will be made available upon reasonable request.

References

  1. Kuo AD, Donelan JM. Dynamic Principles of Gait and Their Clinical Implications. Phys Ther. 2010; 90(2):157–74. https://doi.org/10.2522/ptj.20090125.

    PubMed  PubMed Central  Google Scholar 

  2. Neptune RR, Kautz SA, Zajac FE. Contributions of the individual ankle plantar flexors to support, forward progression and swing initiation during walking. J Biomech. 2001; 34(11):1387–98. https://doi.org/10.1016/S0021-9290(01)00105-1.

    CAS  PubMed  Google Scholar 

  3. Ting LH, Chiel HJ, Trumbower RD, Allen JL, McKay JL, Hackney ME, Kesar TM. Neuromechanical principles underlying movement modularity and their implications for rehabilitation. 2016; 86(1):38–54. https://doi.org/10.1016/j.neuron.2015.02.042.Neuromechanical.

  4. Browne MG, Franz JR. Does dynamic stability govern propulsive force generation in human walking?R Soc Open Sci. 2017; 4(11). https://doi.org/10.1098/rsos.171673.

  5. Bethoux F, Bennett S. Evaluating Walking in Patients with Multiple Sclerosis. Int J MS Care. 2011; 13(1):4–14. https://doi.org/10.7224/1537-2073-13.1.4.

    PubMed  PubMed Central  Google Scholar 

  6. Morris ME, Huxham F, McGinley J, Dodd K, Iansek R. The biomechanics and motor control of gait in Parkinson disease. Clin Biomech. 2001; 16(6):459–70. https://doi.org/10.1016/S0268-0033(01)00035-3.

    CAS  Google Scholar 

  7. Mahon CE, Farris DJ, Sawicki GS, Lewek MD. Individual limb mechanical analysis of gait following stroke. J Biomech. 2015; 48(6):984–9. https://doi.org/10.1016/j.jbiomech.2015.02.006.

    PubMed  Google Scholar 

  8. Franz JR. The Age-Associated Reduction in Propulsive Power Generation in Walking. Exerc Sport Sci Rev. 2016; 44(4):129–36. https://doi.org/10.1249/JES.0000000000000086.

    PubMed  Google Scholar 

  9. Bowden MG, Balasubramanian CK, Neptune RR, Kautz SA. Anterior-Posterior Ground Reaction Forces as a Measure of Paretic Leg Contribution in Hemiparetic Walking. Stroke. 2006; 37(3):872–6. https://doi.org/10.1161/01.STR.0000204063.75779.8d.

    PubMed  Google Scholar 

  10. Peters DM, Thibaudier Y, Deffeyes JE, Baer GT, Hayes HB, Trumbower RD. Constraints on Stance-Phase Force Production during Overground Walking in Persons with Chronic Incomplete Spinal Cord Injury. J Neurotrauma. 2018; 35(3):467–77. https://doi.org/10.1089/neu.2017.5146.

    PubMed  PubMed Central  Google Scholar 

  11. Williams G, Morris ME, Schache A, McCrory PR. Incidence of Gait Abnormalities After Traumatic Brain Injury. Arch Phys Med Rehabil. 2009; 90(4):587–93. https://doi.org/10.1016/j.apmr.2008.10.013.

    PubMed  Google Scholar 

  12. Williams GP, Schache AG, Morris ME. Mobility After Traumatic Brain Injury. J Head Trauma Rehabil. 2013; 28(5):371–8. https://doi.org/10.1097/HTR.0b013e31824a1d40.

    PubMed  Google Scholar 

  13. Ivanenko YP, Cappellini G, Solopova IA, Grishin AA, MacLellan MJ, Poppele RE, Lacquaniti F. Plasticity and modular control of locomotor patterns in neurological disorders with motor deficits. Front Comput Neurosci. 2013; 7(SEP):1–11. https://doi.org/10.3389/fncom.2013.00123.

    Google Scholar 

  14. Donelan JM, Kram R, Kuo AD. Mechanical work for step-to-step transitions is a major determinant of the metabolic cost of human walking. J Exp Biol. 2002; 205(23):3717–27.

    PubMed  Google Scholar 

  15. Kuhman D, Hurt CP. The timing of locomotor propulsion in healthy adults walking at multiple speeds. Hum Mov Sci. 2019; 68(September):102524. https://doi.org/10.1016/j.humov.2019.102524.

    PubMed  Google Scholar 

  16. Zelik KE, Adamczyk PG. A unified perspective on ankle push-off in human walking. J Exp Biol. 2016; 219(23):3676–83. https://doi.org/10.1242/jeb.140376.

    PubMed  PubMed Central  Google Scholar 

  17. Lewek MD. The influence of body weight support on ankle mechanics during treadmill walking. J Biomech. 2011; 44(1):128–33. https://doi.org/10.1016/j.jbiomech.2010.08.037.

    PubMed  Google Scholar 

  18. Riley PO, Croce UD, Casey Kerrigan D. Propulsive adaptation to changing gait speed. J Biomech. 2001; 34(2):197–202. https://doi.org/10.1016/S0021-9290(00)00174-3.

    CAS  PubMed  Google Scholar 

  19. Combs SA, Van Puymbroeck M, Altenburger PA, Miller KK, Dierks TA, Schmid AA. Is walking faster or walking farther more important to persons with chronic stroke?,. Disabil Rehabil. 2013; 35(10):860–7. https://doi.org/10.3109/09638288.2012.717575.

    PubMed  Google Scholar 

  20. Farris D, Hampton A, Lewek MD, Sawicki GS. Revisiting the mechanics and energetics of walking in individuals with chronic hemiparesis following stroke: from individual limbs to lower limb joints. J NeuroEngineering Rehabil. 2015; 12(1):24. https://doi.org/10.1186/s12984-015-0012-x.

    Google Scholar 

  21. Reisman DS, Rudolph KS, Farquhar WB. Influence of speed on walking economy poststroke. Neurorehabil Neural Repair. 2009; 23(6):529–34. https://doi.org/10.1177/1545968308328732.

    PubMed  Google Scholar 

  22. Awad LN, Reisman DS, Pohlig RT, Binder-Macleod SA. Reducing The Cost of Transport and Increasing Walking Distance After Stroke: A Randomized Controlled Trial on Fast Locomotor Training Combined With Functional Electrical Stimulation. Neurorehabil Neural Repair. 2016; 30(7):661–670. https://doi.org/10.1177/1545968315619696.

    PubMed  Google Scholar 

  23. Franz JR, Kram R. How does age affect leg muscle activity/coactivity during uphill and downhill walking?Gait Posture. 2013; 37(3):378–84. https://doi.org/10.1016/j.gaitpost.2012.08.004.

    PubMed  Google Scholar 

  24. Franz JR, Maletis M, Kram R. Real-time feedback enhances forward propulsion during walking in old adults. Clin Biomech (Bristol, Avon). 2014; 29(1):68–74. https://doi.org/10.1016/j.clinbiomech.2013.10.018.

    Google Scholar 

  25. Martin M, Shinberg M, Kuchibhatla M, Ray L, Carollo JJ, Schenkman ML. Gait Initiation in Community-Dwelling Adults With Parkinson Disease: Comparison With Older and Younger Adults Without the Disease. Phys Ther. 2002; 82(6):566–77. https://doi.org/10.1093/ptj/82.6.566.

    PubMed  Google Scholar 

  26. Jonkers I, Delp S, Patten C. Capacity to increase walking speed is limited by impaired hip and ankle power generation in lower functioning persons post-stroke. Gait Posture. 2009; 29(1):129–37. https://doi.org/10.1016/j.gaitpost.2008.07.010.

    CAS  PubMed  Google Scholar 

  27. Turns LJ, Neptune RR, Kautz SA. Relationships Between Muscle Activity and Anteroposterior Ground Reaction Forces in Hemiparetic Walking. Arch Phys Med Rehabil. 2007; 88(9):1127–35. https://doi.org/10.1016/j.apmr.2007.05.027.

    PubMed  PubMed Central  Google Scholar 

  28. Brunnstrom S. Movement Therapy in Hemiplegia: A Neurophysiological Approach. New York: Med Dept Harper Row; 1970.

    Google Scholar 

  29. Hurt CP, Burgess JK, Brown DA. Limb contribution to increased self-selected walking speeds during body weight support in individuals poststroke. Gait Posture. 2015; 41(3):857–9. https://doi.org/10.1016/j.gaitpost.2015.02.004.

    PubMed  PubMed Central  Google Scholar 

  30. Awad LN, Binder-Macleod SA, Pohlig RT, Reisman DS. Paretic propulsion and trailing limb angle are key determinants of long-distance walking function after stroke. Neurorehabil Neural Repair. 2015; 29(6):499–508. https://doi.org/10.1177/1545968314554625.Paretic.

    PubMed  Google Scholar 

  31. Fulk GD, He Y, Boyne P, Dunning K. Predicting Home and Community Walking Activity Poststroke. Stroke. 2017; 48(2):406–11. https://doi.org/10.1161/STROKEAHA.116.015309.

    PubMed  Google Scholar 

  32. Fulk GD, Reynolds C, Mondal S, Deutsch JE. Predicting home and community walking activity in people with stroke. Arch Phys Med Rehabil. 2010; 91(10):1582–86. https://doi.org/10.1016/j.apmr.2010.07.005.

    PubMed  Google Scholar 

  33. Awad LN, Bae J, O’Donnell K, De Rossi SMM, Hendron K, Sloot LH, Kudzia P, Allen S, Holt KG, Ellis TD, Walsh CJ. A soft robotic exosuit improves walking in patients after stroke. Sci Transl Med. 2017; 9(400):9084. https://doi.org/10.1126/scitranslmed.aai9084.

    Google Scholar 

  34. Awad LN, Reisman DS, Kesar TM, Binder-Macleod SA. Targeting Paretic Propulsion to Improve Poststroke Walking Function: A Preliminary Study. Arch Phys Med Rehabil. 2014; 95(5):840–8. https://doi.org/10.1016/j.apmr.2013.12.012.

    PubMed  Google Scholar 

  35. Browne MG, Franz JR. Ankle power biofeedback attenuates the distal-to-proximal redistribution in older adults. Gait Posture. 2019; 71(October 2018):44–9. https://doi.org/10.1016/j.gaitpost.2019.04.011.

    PubMed  PubMed Central  Google Scholar 

  36. Kesar TM, Perumal R, Reisman DS, Jancosko A, Rudolph KS, Higginson JS, Binder-Macleod SA. Functional Electrical Stimulation of Ankle Plantarflexor and Dorsiflexor Muscles. Stroke. 2009; 40(12):3821–7. https://doi.org/10.1161/STROKEAHA.109.560375.

    PubMed  PubMed Central  Google Scholar 

  37. McCain EM, Dick TJM, Giest TN, Nuckols RW, Lewek MD, Saul KR, Sawicki GS. Mechanics and energetics of post-stroke walking aided by a powered ankle exoskeleton with speed-adaptive myoelectric control. J NeuroEngineering Rehabil. 2019; 16(1):57. https://doi.org/10.1186/s12984-019-0523-y.

    Google Scholar 

  38. Penke K, Scott K, Sinskey Y, Lewek MD. Propulsive Forces Applied to the Body’s Center of Mass Affect Metabolic Energetics Poststroke. Arch Phys Med Rehabil. 2019; 100(6):1068–75. https://doi.org/10.1016/j.apmr.2018.10.010.

    PubMed  Google Scholar 

  39. Phadke CP. Immediate Effects of a Single Inclined Treadmill Walking Session on Level Ground Walking in Individuals After Stroke. Am J Phys Med Rehabil. 2012; 91(4):337–45. https://doi.org/10.1097/PHM.0b013e31823cabe3.

    PubMed  Google Scholar 

  40. Takahashi KZ, Lewek MD, Sawicki GS. A neuromechanics-based powered ankle exoskeleton to assist walking post-stroke: a feasibility study. J NeuroEngineering Rehabil. 2015; 12(1):23. https://doi.org/10.1186/s12984-015-0015-7.

    Google Scholar 

  41. Genthe K, Schenck C, Eicholtz S, Zajac-Cox L, Wolf S, Kesar TM. Effects of real-time gait biofeedback on paretic propulsion and gait biomechanics in individuals post-stroke. Top Stroke Rehabil. 2018; 25(3):186–193.

    PubMed  PubMed Central  Google Scholar 

  42. Boutaayamou M, Schwartz C, Stamatakis J, Denoël V. Development and validation of an accelerometer-based method for quantifying gait events. 2015:1–8. https://doi.org/10.1016/j.medengphy.2015.01.001.

  43. Chia Bejarano N, Ambrosini E, Pedrocchi A, Ferrigno G, Monticone M, Ferrante S. A Novel Adaptive, Real-Time Algorithm to Detect Gait Events From Wearable Sensors. IEEE Trans Neural Syst Rehabil Eng. 2015; 23(3):413–22. https://doi.org/10.1109/TNSRE.2014.2337914.

    PubMed  Google Scholar 

  44. Miyazaki T, Kawada M, Nakai Y, Kiyama R, Yone K. Validity of Measurement for Trailing Limb Angle and Propulsion Force during Gait Using a Magnetic Inertial Measurement Unit. BioMed Res Int. 2019; 2019:1–8. https://doi.org/10.1155/2019/8123467.

    Google Scholar 

  45. Peruzzi A, Della Croce U, Cereatti A. Estimation of stride length in level walking using an inertial measurement unit attached to the foot: A validation of the zero velocity assumption during stance. J Biomech. 2011; 44(10):1991–4. https://doi.org/10.1016/j.jbiomech.2011.04.035.

    CAS  PubMed  Google Scholar 

  46. Seel T, Raisch J, Schauer T. IMU-Based Joint Angle Measurement for Gait Analysis. Sensors. 2014; 14(4):6891–909. https://doi.org/10.3390/s140406891.

    PubMed  PubMed Central  Google Scholar 

  47. Yang S, Zhang J-T, Novak AC, Brouwer B, Li Q. Estimation of spatio-temporal parameters for post-stroke hemiparetic gait using inertial sensors. Gait Posture. 2013; 37(3):354–8. https://doi.org/10.1016/j.gaitpost.2012.07.032.

    PubMed  Google Scholar 

  48. Karatsidis A, Jung M, Schepers HM, Bellusci G, de Zee M, Veltink PH, Andersen MS. Musculoskeletal model-based inverse dynamic analysis under ambulatory conditions using inertial motion capture. Med Eng Phys. 2019; 65:68–77. https://doi.org/10.1016/j.medengphy.2018.12.021.

    PubMed  Google Scholar 

  49. Lim, Kim, Park. Prediction of Lower Limb Kinetics and Kinematics during Walking by a Single IMU on the Lower Back Using Machine Learning. Sensors. 2019; 20(1):130. https://doi.org/10.3390/s20010130.

    PubMed Central  Google Scholar 

  50. Ryu HX, Park S. Estimation of unmeasured ground reaction force data based on the oscillatory characteristics of the center of mass during human walking. J Biomech. 2018; 71:135–43. https://doi.org/10.1016/j.jbiomech.2018.01.046.

    PubMed  Google Scholar 

  51. Shahabpoor E, Pavic A. Estimation of vertical walking ground reaction force in real-life environments using single IMU sensor. J Biomech. 2018; 79:181–90. https://doi.org/10.1016/j.jbiomech.2018.08.015.

    CAS  PubMed  Google Scholar 

  52. Ancillao A, Tedesco S, Barton J, O’Flynn B. Indirect Measurement of Ground Reaction Forces and Moments by Means of Wearable Inertial Sensors: A Systematic Review. Sensors. 2018; 18(8):2564. https://doi.org/10.3390/s18082564.

    PubMed Central  Google Scholar 

  53. Pieper NL, Lewek MD, Franz JR. Can shank acceleration provide a clinically feasible surrogate for individual limb propulsion during walking?J Biomech. 2020; 98:109449. https://doi.org/10.1016/j.jbiomech.2019.109449.

    PubMed  Google Scholar 

  54. Pohl PS, Duncan PW, Perera S, Liu W, Lai SM, Studenski S, Long J. Influence of stroke-related impairments on performance in 6-minute walk test. J Rehabil Res Dev. 2002; 39(4):439–44.

    PubMed  Google Scholar 

  55. Duclos NC, Duclos C, Nadeau S. Slow and faster post-stroke walkers have a different trunk progression and braking impulse during gait. Gait Posture. 2019; 68(July 2018):483–7. https://doi.org/10.1016/j.gaitpost.2018.12.037.

    PubMed  Google Scholar 

  56. Lewek MD, Sawicki GS. Trailing limb angle is a surrogate for propulsive limb forces during walking post-stroke. Clin Biomech. 2019; 67(May):115–8. https://doi.org/10.1016/j.clinbiomech.2019.05.011.

    Google Scholar 

  57. Neptune RR, Zajac FE, Kautz SA. Muscle mechanical work requirements during normal walking: The energetic cost of raising the body’s center-of-mass is significant. J Biomech. 2004; 37(6):817–25. https://doi.org/10.1016/j.jbiomech.2003.11.001.

    CAS  PubMed  Google Scholar 

  58. Xsens Technologies BV. MTi User Manual. Technical report. Enschede: Xsens; 2019.

    Google Scholar 

  59. Liu T, Inoue Y, Shibata K. Development of a wearable sensor system for quantitative gait analysis. Measurement. 2009; 42(7):978–88. https://doi.org/10.1016/j.measurement.2009.02.002.

    Google Scholar 

  60. Takeda R, Lisco G, Fujisawa T, Gastaldi L, Tohyama H, Tadano S. Drift Removal for Improving the Accuracy of Gait Parameters Using Wearable Sensor Systems. Sensors. 2014; 14(12):23230–47. https://doi.org/10.3390/s141223230.

    PubMed  PubMed Central  Google Scholar 

  61. Hamilton WR, Hamilton WE. Elements of Quaternions: Longmans, Green, & Company; 1866. https://books.google.com/books?id=GFYtAAAAYAAJ.

  62. Henderson M-D. Euler Angles, Quaternions, and Transformation matrics for space shuttle analysis. Technical report. Houston: McDonnell-Douglas Technical Services Co., Inc.; Astronautics Div; 1977.

    Google Scholar 

  63. Sabatini AM. Quaternion-based strap-down integration method for applications of inertial sensing to gait analysis. Med Biol Eng Comput. 2005; 43(1):94–101. https://doi.org/10.1007/BF02345128.

    CAS  PubMed  Google Scholar 

  64. Skog I, Handel P, Nilsson JO, Rantakokko J. Zero-Velocity Detection—An Algorithm Evaluation. IEEE Trans Biomed Eng. 2010; 57(11):2657–66. https://doi.org/10.1109/TBME.2010.2060723.

    Google Scholar 

  65. Winter DA. Biomechanics and Motor Control of Human Movement; 2009, pp. 1–367.

  66. Koo TK, Li MY. A Guideline of Selecting and Reporting Intraclass Correlation Coefficients for Reliability Research. J Chiropr Med. 2016; 15(2):155–63. https://doi.org/10.1016/j.jcm.2016.02.012.

    PubMed  PubMed Central  Google Scholar 

  67. Mehta S, Bastero-Caballero RF, Sun Y, Zhu R, Murphy DK, Hardas B, Koch G. Performance of intraclass correlation coefficient (ICC) as a reliability index under various distributions in scale reliability studies. Stat Med. 2018. https://doi.org/10.1002/sim.7679.

  68. McGraw KO, Wong SP. "Forming inferences about some intraclass correlations coefficients": Correction,. Psychol Methods. 1996; 1(4):390. https://doi.org/10.1037/1082-989X.1.4.390.

    Google Scholar 

  69. Chen G, Patten C, Kothari DH, Zajac FE. Gait differences between individuals with post-stroke hemiparesis and non-disabled controls at matched speeds. Gait Posture. 2005; 22(1):51–6. https://doi.org/10.1016/j.gaitpost.2004.06.009.

    PubMed  Google Scholar 

  70. Harris-Love ML, Forrester LW, Macko RF, Silver KHC, Smith GV. Hemiparetic Gait Parameters in Overground Versus Treadmill Walking. Neurorehabil Neural Repair. 2001; 15(2):105–12. https://doi.org/10.1177/154596830101500204.

    CAS  PubMed  Google Scholar 

  71. Harris-Love ML, Macko RF, Whitall J, Forrester LW. Improved Hemiparetic Muscle Activation in Treadmill versus Overground Walking. Neurorehabil Neural Repair. 2004; 18(3):154–60. https://doi.org/10.1177/0888439004267678.

    PubMed  Google Scholar 

  72. Kautz SA, Bowden MG, Clark DJ, Neptune RR. Comparison of Motor Control Deficits During Treadmill and Overground Walking Poststroke. Neurorehabil Neural Repair. 2011; 25(8):756–65. https://doi.org/10.1177/1545968311407515.

    PubMed  PubMed Central  Google Scholar 

  73. Kesar TM, Binder-Macleod SA, Hicks GE, Reisman DS. Minimal detectable change for gait variables collected during treadmill walking in individuals post-stroke. Gait Posture. 2011; 33(2):314–7. https://doi.org/10.1016/j.gaitpost.2010.11.024.

    PubMed  Google Scholar 

  74. Brouwer B, Parvataneni K, Olney SJ. A comparison of gait biomechanics and metabolic requirements of overground and treadmill walking in people with stroke. Clin Biomech. 2009; 24(9):729–34. https://doi.org/10.1016/j.clinbiomech.2009.07.004.

    Google Scholar 

  75. Campanini I, Merlo A. Reliabilty, smallest real difference and concurrent validity of indices computed from GRF components in gait of stroke patients. Gait Posture. 2009; 30(2):127–31. https://doi.org/10.1016/j.gaitpost.2009.03.011.

    CAS  PubMed  Google Scholar 

  76. Hass CJ, Buckley TA, Pitsikoulis C, Barthelemy EJ. Progressive resistance training improves gait initiation in individuals with Parkinson’s disease. Gait Posture. 2012; 35(4):669–73. https://doi.org/10.1016/j.gaitpost.2011.12.022.

    PubMed  Google Scholar 

  77. Ferraris F. Procedure for effortless in-field calibration of three-axis rate gyros and accelerometers. Sensors Mater. 1995; 7(5):311–30.

    Google Scholar 

Download references

Funding

The authors acknowledge funding support from the NIH (KL2TR001411) and AHA (18TPA34170171 and 18IPA34170487).

Author information

Authors and Affiliations

Authors

Contributions

D.A.R., A.M.A., C.J.W., S.D.R, and L.N.A. conceived the study concepts and experimental methods; D.A.R. and A.M.A collected the data; D.A.R., S.D.R., and L.N.A. conceived and implemented the algorithms; D.A.R. conducted the analyses; D.A.R., A.M.A, and L.N.A. prepared the manuscript; all authors revised and approved the final manuscript.

Corresponding author

Correspondence to Louis N. Awad.

Ethics declarations

Ethics approval and consent to participate

All study procedures were approved by the Boston University Institutional Review Board. Written informed consent was secured from all study participants.

Consent for publication

All study participants provided consent for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1

Supplementary Table 1: Changes in R2 and RMSE with removal of single IMUs. Supplementary Table 2: ICC 95% confidence intervals for propulsion and braking point metrics.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Revi, D.A., Alvarez, A.M., Walsh, C.J. et al. Indirect measurement of anterior-posterior ground reaction forces using a minimal set of wearable inertial sensors: from healthy to hemiparetic walking. J NeuroEngineering Rehabil 17, 82 (2020). https://doi.org/10.1186/s12984-020-00700-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12984-020-00700-7

Keywords