[go: up one dir, main page]

Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Kinematic Validation of a Multi-Kinect v2 Instrumented 10-Meter Walkway for Quantitative Gait Assessments

  • Daphne J. Geerse ,

    d.j.geerse@vu.nl

    Affiliations MOVE Research Institute Amsterdam, Department of Human Movement Sciences, VU University Amsterdam, Amsterdam, The Netherlands, Department of Neurology, Leiden University Medical Center, Leiden, The Netherlands

  • Bert H. Coolen,

    Affiliation MOVE Research Institute Amsterdam, Department of Human Movement Sciences, VU University Amsterdam, Amsterdam, The Netherlands

  • Melvyn Roerdink

    Affiliation MOVE Research Institute Amsterdam, Department of Human Movement Sciences, VU University Amsterdam, Amsterdam, The Netherlands

Abstract

Walking ability is frequently assessed with the 10-meter walking test (10MWT), which may be instrumented with multiple Kinect v2 sensors to complement the typical stopwatch-based time to walk 10 meters with quantitative gait information derived from Kinect’s 3D body point’s time series. The current study aimed to evaluate a multi-Kinect v2 set-up for quantitative gait assessments during the 10MWT against a gold-standard motion-registration system by determining between-systems agreement for body point’s time series, spatiotemporal gait parameters and the time to walk 10 meters. To this end, the 10MWT was conducted at comfortable and maximum walking speed, while 3D full-body kinematics was concurrently recorded with the multi-Kinect v2 set-up and the Optotrak motion-registration system (i.e., the gold standard). Between-systems agreement for body point’s time series was assessed with the intraclass correlation coefficient (ICC). Between-systems agreement was similarly determined for the gait parameters’ walking speed, cadence, step length, stride length, step width, step time, stride time (all obtained for the intermediate 6 meters) and the time to walk 10 meters, complemented by Bland-Altman’s bias and limits of agreement. Body point’s time series agreed well between the motion-registration systems, particularly so for body points in motion. For both comfortable and maximum walking speeds, the between-systems agreement for the time to walk 10 meters and all gait parameters except step width was high (ICC ≥ 0.888), with negligible biases and narrow limits of agreement. Hence, body point’s time series and gait parameters obtained with a multi-Kinect v2 set-up match well with those derived with a gold standard in 3D measurement accuracy. Future studies are recommended to test the clinical utility of the multi-Kinect v2 set-up to automate 10MWT assessments, thereby complementing the time to walk 10 meters with reliable spatiotemporal gait parameters obtained objectively in a quick, unobtrusive and patient-friendly manner.

Introduction

Walking speed is associated with falls [13], adverse events [4,5] and life expectancy [6] in older adults. A standardized clinical test often used to assess walking speed is the 10-meter walking test (10MWT). However, the 10MWT only provides a single performance measure (i.e., walking speed derived from the time to walk 10 meters), reflecting just one aspect of walking ability. To yield a more comprehensive evaluation of walking ability, quantitative gait assessments (e.g., step length, cadence and step width) may be conducted using high-end motion-registration systems. Yet, even the best motion-registration systems yield limitations when conducting quantitative gait assessments in clinical settings (e.g., costs, patient-preparation time, calibration procedures, marker occlusion, and delays in availability of results [7]).

A promising motion-registration system to instrument the 10MWT is the Microsoft Kinect sensor, a RGB-D camera that was launched in 2011 in combination with a Software Development Kit for 3D human-pose estimation, originating from the gaming industry [8]. The development of 3D human-pose estimation software, using a large and highly varied training dataset of paired depth images and ground truth body parts to train very deep decision forests for efficient and accurate body part recognition [8], was a major undertaking by Microsoft. It successfully eliminated the need for markers and calibration procedures, thereby enabling fast and patient-friendly 3D full-body motion registration (Fig 1). This motion-registration system has gained enormous interest from developers and scientists in the context of assessment and rehabilitation of balance, posture and gait (e.g., [918]), since it allows for motion registration in a quick and affordable manner. Recently, the second generation of the Kinect sensor has been introduced. Key differences with the previous Kinect v1 sensor are that the Kinect v2 sensor is a time-of-flight camera with an increased resolution of the depth image, a wider field of view and improved body point tracking [19], possibly leading to improved results.

thumbnail
Fig 1. Body points derived with the human-pose estimation software of Kinect v1.

(A) RGB image and (B) depth image with the corresponding body points derived with the human-pose estimation software of Kinect v1.

https://doi.org/10.1371/journal.pone.0139913.g001

Several studies have demonstrated that spatiotemporal gait parameters can be validly obtained using a single Kinect v1 sensor [9,11,13,14,17], and recently also for a single Kinect v2 sensor [15]. However, these studies only analyzed a few steps since accurate body point tracking with the Kinect sensor is only possible between 0.8 and 4.0 meters from the Kinect v1 sensor and between 0.5 and 4.5 meters from the Kinect v2 sensor due to the limited field of view and poorer depth-image quality at greater distances. One way to cover a larger volume, such as the walkway of the 10MWT, is to use multiple spatially and temporally integrated Kinect sensors. Hereby measurement volume may be increased, while preserving good quality depth images for accurate body point tracking. This supposedly allows for the parametrization of a large number of steps during walking from high quality 3D body point’s time series. In view of Kinect’s v2 higher resolution depth images, improved body point tracking and enlarged area for accurate body point tracking, the current study will explore the potential of a multi-Kinect v2 set-up for instrumenting the 10MWT.

The objective of this study is to determine the usability of a multi-Kinect v2 set-up to quantitatively assess gait during the 10MWT. Because the multi-Kinect v2 set-up has not yet been validated for 3D full-body motion registration, its performance will be compared to a gold standard in 3D measurement accuracy (i.e., the Optotrak active-marker 3D optical tracking system, Northern Digital Inc., Waterloo, Canada). The between-systems agreement will be examined for raw data (i.e., body point’s time series) and spatiotemporal gait parameters (e.g., step length, cadence and step width). In addition, the between-systems agreement for the performance measure of the 10MWT (i.e., time to walk 10 meters) will be assessed between the multi-Kinect v2 set-up, the Optotrak motion-registration system (i.e., the gold-standard reference) and the stopwatch (i.e., the clinical standard).

Methods

Subjects

A heterogeneous group of 21 healthy subjects in terms of gender (11 males, 10 females), age (mean [range]: 30.2 [19–63] years), height (176.1 [158–190] cm) and weight (70.5 [53–83] kg) took part in this experiment. Subjects did not have any medical condition that would influence walking.

Ethics statement

The current study was approved by the ethics committee of the Department of Human Movement Sciences (VU University Amsterdam, Amsterdam). All subjects provided written informed consent prior to participation. The subjects in Fig 1 have given written informed consent, as outlined in the PLOS consent form, to publish this photograph.

Experimental set-up and procedure

Full-body kinematics was recorded with four spatially and temporally integrated Microsoft Kinect v2 sensors and the Optotrak system (Northern Digital Inc., Waterloo, Canada). The multi-Kinect v2 set-up is displayed in Fig 2. The four Kinect v2 sensors were positioned on tripods alongside a walkway of 10 by 0.5 meters at a height of 0.75 meters. The sensors were placed 0.5 meters from the left border of the walkway with an angle of 70 degrees relative to the walkway direction. The first sensor was positioned at 4 meters from the start of the walkway. The other three sensors were placed at inter-sensor distances of 2.5 meters. In addition, five Optotrak cameras (i.e., a combination of two Optotrak 3020 and three Optotrak Certus cameras, which are all compatible with each other) were positioned around the walkway to cover the same area as the multi-Kinect v2 set-up. The so-obtained Optotrak set-up ensured sub-millimeter accuracy throughout the 10-meter walkway. The coordinate systems of the multi-Kinect v2 set-up and the Optotrak system were aligned using a spatial calibration grid.

The Kinect for Windows Software Development Kit (SDK 2.0, www.microsoft.com) provides, with a sampling rate of 30 Hz, the 3D positions of 25 body points (Fig 3B). These body points are: head, neck, spine shoulder, spine mid, spine base and right and left shoulder, elbow, wrist, hand, thumb, hand tip, hip, knee, ankle and foot. For motion registration with the Optotrak system (Northern Digital Inc., Waterloo, Canada, using First Principles data acquisition software with a sampling rate of 60 Hz), subjects were asked to wear tight-fitting shorts and a t-shirt to limit clothing-related marker occlusion. Smart Marker Rigid Bodies (Northern Digital Inc., Waterloo, Canada) were attached to the head, upper arms, forearms, lower abdomen, upper legs, lower legs and feet (Fig 3A), allowing for 6 degrees of freedom tracking of body segments. In addition, 30 anatomical landmarks were digitized using a 3-marker digitizing probe to define various body point positions (so-called virtual markers) on abovementioned body segments. Smart markers were also placed on the sternum, hands and feet. The body points represented by Optotrak’s virtual markers and/or smart markers were selected to closely match Kinect’s body points (see S1 Table), although sometimes arbitrary positional differences between the body point’s time series of the two motion-registration systems could not be prevented because 1) the exact definitions of the body points given by the human-pose estimation algorithms of Kinect v2 are not known and 2) virtual markers and smart markers are by definition positioned at the contours of the body while Kinect v2 body points are typically estimated within the body. For example, the smart marker representing Kinect’s spine shoulder was placed on the sternum (see S1 Table), which deviates in AP direction from the within-body spine shoulder given by the human-pose estimation algorithm of Kinect v2, thus resulting in a between-systems positional mismatch. Positions of the neck, spine mid, thumbs and hand tips body points were not tracked with the Optotrak system due to the limited number of available smart markers, rendering a total of 19 out of aforementioned 25 body points eligible for a between-systems agreement analysis (as specified in S1 Table).

thumbnail
Fig 3. Body point determination with the Optotrak and Kinect v2 systems.

(A) Subject with all markers of the Optotrak system; (B) Same subject with body points derived with the human-pose estimation algorithm of Kinect v2.

https://doi.org/10.1371/journal.pone.0139913.g003

Before conducting the experiment, the quality of depth image of the subject was checked since some textiles are known to corrupt the infrared radiation emitted by the previous Kinect v1 sensor, making human-pose estimation less accurate [17]. No problems were encountered with clothing of the subjects, possibly owing to the improved properties of the Kinect v2 sensor. Subsequently, subjects performed the 10MWT at two different walking speeds, namely comfortable walking speed (CWS) and maximum walking speed (MWS). Both conditions were performed three times in a fixed order (i.e., three times CWS followed by three times MWS). Subjects were instructed to start walking at the fourth, high-pitched beep of a standardized auditory start command (i.e., three low-pitched beeps followed by one high-pitched beep) and to continue walking until they had fully crossed the finish line. The standardized auditory start command was synchronized with the multi-Kinect v2 set-up. Synchronization between the two motion-registration systems was achieved by a synchronization movement (i.e., ab- and adduction of both arms) that participants performed prior the auditory start command of each trial. Motion registration started before the synchronization movement and ended well after the subject had passed the 10-meter line. Time to walk 10 meters (i.e., from final beep onset until the moment that the most forward ankle passed the 10-meter line, according to the recommendations of Graham et al. [20]) was determined using a stopwatch. A video showing body point’s time series simultaneously for both measurement systems during the 10MWT is available in S1 Video. This video also includes the synchronization movement and the standardized auditory start command.

Data pre-processing

The 3D positional data of body points were first pre-processed per Kinect sensor separately. Inferred body points (i.e., when a body point was not visible due to for example occlusion, Kinect’s human-pose estimation software inferred its position) were considered as missing values. Moreover, since the sampling frequency of the Kinect system is not constant (i.e., apart from 20 outliers in inter-sample intervals for multiple subjects but confined to one Kinect sensor, the remaining inter-sample intervals ranged from 32 to 34 ms), the body point’s time series were linearly interpolated using Kinect’s timestamps to ensure a constant sampling frequency of 30 Hz, without filling in the parts with missing values. Data points not adhering to the requirements for valid human-pose estimation (e.g., minimum of 15 tracked body points out of the 25 body points, tracked data points for the head and at least one foot and no outliers in segment lengths) were removed from the time series. Subsequently, data of the four Kinect sensors were combined by taking for each sample the 3D positions of the body points of a validly estimated human pose. If, for a given sample, more than one sensor contained valid human-pose data, the associated body point’s 3D positions were averaged for that specific sample. Optotrak data were down-sampled to 30 Hz. Subsequently, the cross-covariance and time lag were determined for paired time series in the mediolateral (ML) and vertical (V) direction of the elbows, wrists and hands during the synchronization movement. These time series were first interpolated with a spline algorithm in case of missing data. The median of the time lags was used to temporally align the time series of the two motion-registration systems. Time-synchronized 3D body point’s time series of both systems are presented as supplementary material, starting from final beep onset until the moment that for both systems the most forward ankle passed the 10-meter line (see S1 Data). Body point’s time series with more than 50 percent of missing values were excluded from further analyses. No time series were excluded for the multi-Kinect v2 set-up, whereas 17 out of 2394 time series were excluded for Optotrak, including two time series of the ankles from which gait parameters were derived. The missing values of the remaining data were interpolated with a spline algorithm. The so-obtained time series were used for assessing the between-systems agreement in body point’s time series (see Data analysis) and for the quantification of several gait parameters, as specified in the next paragraph.

Several gait parameters were calculated from the body point’s time series, separately for both measurement systems. The following spatiotemporal gait parameters were all determined for the intermediate 6 meters (i.e., from the 2-meter to the 8-meter line), reducing the effect of gait acceleration and deceleration on the gait parameters [21]. Walking speed (in cm/s) was defined as the distance travelled between the 2-meter and 8-meter line on the walkway divided by the time, using the data of the spine shoulder. For the other gait parameters, estimates of foot contact and foot off were required, stemming from respectively the maxima and minima of the anterior-posterior (AP) time series of the ankles relative to that of the spine base [22] (Fig 4A and 4C). For spatial gait parameters, first right and left step locations were determined, defined as the median value of the right and left ankle position in the AP and ML direction during the respective single-support stance phases (i.e., between foot off and foot contact of the contralateral foot). Based on these AP and ML step locations, various spatial gait parameters were determined. Step length (in cm) was calculated as the AP difference of consecutive step locations (Fig 4D). Stride length (in cm) was calculated as the AP difference of consecutive ipsilateral step locations. Moreover, step width (in cm) was estimated by taking the absolute ML difference of consecutive step locations. Cadence (in steps/min) was calculated from the number of steps in the time interval between the first and last estimate of foot contact. Step time (in s) was calculated as the time interval between two consecutive instants of foot contact (Fig 4D). Consequently, stride time (in s) was calculated as the time interval between two consecutive ipsilateral instants of foot contact. For step length, stride length, step width, step time and stride time, median values within the 6-meter window were used as outcome measures per trial since Baldewijns et al. [9] demonstrated superior agreement between registration systems on a per walk basis.

thumbnail
Fig 4. Overview of the analysis of spatiotemporal gait parameters.

Analyses for CWS (panels A and B) and MWS (panels C and D) conditions are based on AP displacement data of the right (gray lines) and left (black lines) ankles as a function of time for the multi-Kinect v2 set-up. AP ankle time series relative to the spine base (panels A and C) were used to estimate instants of foot contact (black dots) and foot off (gray dots) for each step. Step location was defined as the median value of the AP ankle time series during the single-support stance phase (i.e., the horizontal plateaus delimited by foot off and foot contact events of the contralateral foot in panels B and D). Vertical bars represent the four beep onsets of the auditory start command. The shaded area in panels B and D represent the 6-meter window from which spatiotemporal gait parameters were derived. Dashed lines in panels B and D schematically define the time to walk 10 meters (T10), step time (ST) and step length (SL).

https://doi.org/10.1371/journal.pone.0139913.g004

The performance measure of the 10MWT, that is the time to walk 10 meters (in s), was defined as the time from final beep onset until the moment that the most forward ankle passed the 10-meter line (Fig 4B). For comparison with the stopwatch score, serving as the clinical reference, the time to walk 10 meters was also determined from data of the multi-Kinect v2 set-up and the Optotrak system, the latter serving as the gold-standard reference.

Data analysis

First, the between-systems agreement was calculated for the body point’s time series from final beep onset until the moment that the most forward ankle passed the 10-meter line. For the AP direction, the trend was removed using a bidirectional, second-order Butterworth high-pass filter (cutoff frequency of 0.5 Hz) to reduce the effect of a large within-subject variation (increasing from 0 to 10 meter) on the agreement statistic, which would become arbitrarily high [23]. The agreement between the time series of the two motion-registration systems was calculated for each body point in the AP, ML and V direction by means of the intraclass correlation coefficient for consistency (ICC(C,1); [24]). We selected ICC(C,1) in view of abovementioned somewhat arbitrary between-systems mismatches in body point’s time series (see S1 Table). The average ICC(C,1) was constructed over all trials per system, body point and direction for each subject. From these values, the average ICC(C,1) over subjects was calculated for each system, body point and direction, including confidence intervals.

Second, the between-systems agreement for spatiotemporal gait parameters was calculated. Spatiotemporal gait parameters were based on specific within-system time series’ features (e.g., minima or maxima, consecutive step locations) and hence less susceptible to arbitrary systematic between-systems positional differences in body point’s time series. Therefore, the ICC for absolute agreement (ICC(A,1); [24]) was selected. The agreement in the time to walk 10 meters obtained with the multi-Kinect v2 set-up, the Optotrak system (gold standard) and a stopwatch (clinical standard) was also assessed using ICC(A,1).

In line with Cicchetti [25], we regard ICC values above 0.60 as good and ICC values above 0.75 as excellent. ICC(A,1) values were complemented by mean differences and precision values obtained with a Bland-Altman analysis (i.e., the bias and the limits of agreement, respectively; [26]). Since large differences were expected between CWS and MWS conditions for all gait parameters, leading to large within-subject variation that would arbitrarily inflate the between-systems agreement [23], the agreement for gait parameters and time to walk 10 meters was analyzed separately for both conditions. In line with Flansbjer et al. [27], the average time to walk 10 meters was constructed over the three trials per condition per subject. For the spatiotemporal gait parameters the average was hence also constructed over the three trials per condition per subject. For each condition, at least two trials had to be valid (i.e., less than 50 percent of missing values and, for the time to walk 10 meters, data around the 10-meter line and no error in pressing the stopwatch) in order to compute the average over the trials. This resulted in the exclusion of one subject for further analysis of the between-systems agreement for the time to walk 10 meters for the MWS condition.

Results

Agreement between body point’s time series

The agreement (ICC(C,1)) between the body point’s time series of the multi-Kinect v2 set-up and the gold-standard Optotrak motion-registration system for all 19 matched body points in AP (detrended), ML and V directions are listed in Table 1. Apart from the hips, there was a good to excellent agreement in body point’s time series between the two motion-registration systems in the AP direction. Furthermore, all gait parameters were derived from time series with high (i.e., ML time series of the right ankle) or excellent levels of agreement (all other time series), as highlighted in Table 1 (bold values). Fig 5 shows an example of a part of the AP (detrended) and ML time series of the right and left ankle for the multi-Kinect v2 set-up and the Optotrak system during a CWS trial with corresponding ICC(C,1) values (as well as ICC(A,1) values to illustrate the effect of a systematic between-systems mismatch in body point’s time series on ICC values).

thumbnail
Fig 5. Representative ankle time series of the multi-Kinect v2 set-up and the Optotrak system.

Multi-Kinect v2 (gray lines) and Optotrak (black lines) time series of the right (panels A and C) and left (panels B and D) ankle in the AP (detrended) and ML direction for a part of a CWS trial, including between-systems agreement assessed with the intraclass correlation coefficient for consistency (ICC(C,1)) and absolute agreement (ICC(A,1)).

https://doi.org/10.1371/journal.pone.0139913.g005

thumbnail
Table 1. Agreement statistics for body point’s time series.

https://doi.org/10.1371/journal.pone.0139913.t001

Agreement of spatiotemporal gait parameters

The agreement statistics of the spatiotemporal gait parameters are presented in Table 2. Apart from step width, the between-systems agreement for spatiotemporal gait parameters was excellent for CWS (ICC(A,1) ≥ 0.888) and MWS (ICC(A,1) ≥ 0.951) conditions. This was supported by relatively small biases and narrow limits of agreement (Table 2). Step width showed a good between-systems agreement (CWS: 0.646, MWS: 0.705) with proportionally higher biases and wider limits of agreement (Table 2). Bland-Altman plots for spatiotemporal gait parameters are available in S1 Fig.

thumbnail
Table 2. Agreement statistics for spatiotemporal gait parameters.

https://doi.org/10.1371/journal.pone.0139913.t002

Agreement of time to walk 10 meters

Mean values of the time to walk 10 meters for CWS and MWS conditions are presented in Fig 6. There was a high level of agreement between the measurement systems according to the ICC(A,1) for both conditions. For the multi-Kinect v2 set-up and the Optotrak system, ICC(A,1) values were excellent for CWS (ICC(A,1) = 0.998) and MWS (ICC(A,1) = 0.999), with biases being smaller than one sample (CWS: -0.01 s, MWS: -0.01 s) and narrow limits of agreement (CWS: [-0.11 0.09] s, MWS: [-0.07 0.06] s). The comparison between the multi-Kinect v2 set-up and the stopwatch also revealed excellent ICC(A,1) values (CWS: 0.988, MWS: 0.989), but biases were greater (CWS: -0.09 s, MWS: -0.08 s) and limits of agreement wider (CWS: [-0.23 0.05] s, MWS: [-0.21 0.06] s). The same was true for the comparison between the Optotrak system and the stopwatch: excellent ICC(A,1) values (CWS: 0.987, MWS: 0.990) but biases were approximately two samples (CWS: -0.08 s, MWS: -0.07 s) and limits of agreement were again wider (CWS: [-0.26 0.11] s, MWS: [-0.21 0.07] s).

thumbnail
Fig 6. Time to walk 10 meters for CWS and MWS conditions.

Bars represent average time to walk 10 meters for the multi-Kinect v2 set-up (gray bars), the Optotrak motion-registration system as the gold-standard reference (black bars) and the stopwatch as the clinical standard (white bars).

https://doi.org/10.1371/journal.pone.0139913.g006

Discussion

In the current study, we evaluated a multi-Kinect v2 set-up for quantitative gait assessment during the 10MWT by determining between-systems agreement for body point’s time series, for spatiotemporal gait parameters and for the time to walk 10 meters. Performance of the multi-Kinect v2 set-up was compared to the Optotrak system (i.e., the gold-standard reference) to validate 3D full-body kinematical data of the just-released Kinect v2 sensor. We observed a good to excellent agreement between the two motion-registration systems for raw data (i.e., relevant body point’s time series), spatiotemporal gait parameters and the time to walk 10 meters.

To the best of our knowledge, this study is the first to statistically compare unfiltered body point’s time series stemming from a multi-Kinect v2 set-up to a gold-standard reference. Covering the entire measurement volume with a marker-based motion-registration system was quite difficult and required many cameras to avoid marker occlusion. In fact, the number of excluded body point’s time series due to excessive missing values was substantially larger for the marker-based gold standard in 3D measurement accuracy (17 excluded time series, average percentage of missing values was 6.8%) than for the multiple-Kinect v2 set-up (no excluded time series, average percentage of missing values was 5.0%). For the remaining 2377 time series, ICC(C,1) values were generally exceeding 0.60 for all directions, indicating a good to excellent between-systems agreement. Nevertheless, some time series only demonstrated a poor to fair between-systems agreement, especially time series exhibiting a small range of motion. Note that the ICC is constructed using models that assume equal variance between two variables [24]. With a small range of motion (i.e., with low signal power and hence low true within-system variation), the noisier Kinect v2 data may have caused the error-variances of the two motion-registration systems to differ, with consequently a lower between-systems agreement. This is supported by results of a previous study [28], showing that larger movements of Parkinson’s disease patients were better tracked by a Kinect v1 sensor than smaller movements. Thus, as long as body points are moving (i.e., high signal power), the resultant time series of Kinect v2 match well with those stemming from a gold standard in 3D measurement accuracy. Furthermore, low-pass filtering time series may also increase the between-systems agreement.

In the current study, all spatiotemporal gait parameters were derived from body point’s time series with high (for the ML time series of the right ankle) or excellent levels of agreement (for all other time series; see Table 1, bold values). This resulted in excellent between-systems agreement (high ICC(A,1) values) of the from these time series derived spatiotemporal gait parameters walking speed, cadence, step length, stride length, step time and stride time. These spatiotemporal gait parameters can be accurately obtained with the multi-Kinect v2 set-up, as testified by negligible biases and narrow limits of agreement (Table 2). Step width was the only gait parameter that demonstrated good instead of excellent absolute agreement (Table 2). The deviant findings for step width may be due to systematic within-subject differences in ML ankle position time series between the two motion-registration systems. An example of such a systematic positional difference is presented in Fig 5. The left ML ankle position obtained with the multi-Kinect v2 set-up was about 3 to 4 centimeters more lateral compared to Optotrak’s left ML ankle position (Fig 5D) while the right ML ankle positions matched well between the two systems (Fig 5C), resulting in a substantial bias of 3.6 cm in step width for this specific subject. This systematic between-systems mismatch for the left ML ankle position was confirmed by a clear difference between ICC values for consistency and absolute agreement (ICC(C,1) = 0.830, ICC(A,1) = 0.405; Fig 5D), whereas for the right ML ankle positions the ICC values were similar (ICC(C,1) = 0.818, ICC(A,1) = 0.783; Fig 5C). Note that this positional mismatch in ankle time series was not consistent among subjects in terms of its size, sign and side, which may explain the relatively larger between-subjects variation in the between-systems difference for step width (i.e., relatively wider limits of agreement in Table 2).

Kitsunezaki et al. [29] also assessed the possibility of instrumenting the 10MWT with multiple Kinect sensors. Specifically, they used two temporally integrated Kinect v1 sensors that were positioned at the 2-meter and 8-meter lines of a 10-meter walkway to determine the walking time of the intermediate 6 meters of the 10MWT. The mean difference in walking times obtained with the clinical standard (i.e., stopwatch) and the two Kinect v1 sensors was 0.15 seconds, which led the authors to conclude that a Kinect-based assessment was acceptable for practical use [29]. In the current study we quantified the time to walk 10 meters with a multi-Kinect v2 set-up, a gold-standard motion-registration system and a stopwatch. Despite examining walking time over a greater walking distance than Kitsunezaki et al. [29], we found smaller differences between the three measurement systems (≤ 0.09 s), especially between the multi-Kinect v2 set-up and the gold-standard motion-registration system (0.01 s). Noteworthy is that the agreement between these two motion-registration systems–in terms of ICC(A,1), biases and limits of agreement–was better than the agreement of either one with the clinical standard (i.e., stopwatch). To put these findings in perspective, the between-systems differences in the time to walk 10 meters were about 30 to 300 times smaller than the within-system differences between CWS and MWS conditions. Moreover, the meaningful change in walking speed of 5 cm/s according to Perera et al. [30] is at least twice as large as the between-systems differences in walking speed observed in the current study (i.e., after transforming the time to walk 10 meters to walking speed, ≤ 2.5 cm/s).

A multi-Kinect v2 set-up, such as the one described in the current study, may in practice be employed to automate the assessment of the 10MWT. An advantage of this set-up is that the 10MWT and quantitative gait assessment can be conducted simultaneously to reduce the time needed for a comprehensive assessment of walking ability. This could be beneficial for clinical applications, especially in view of our observation that the set-up can provide reliable estimates of the time to walk 10 meters and commonly used spatiotemporal gait parameters in a very quick, unobtrusive and patient-friendly manner. Other advantages of the Kinect v2 sensor are that 3D positional data of 25 body points (of up to six persons!) are tracked and available in real time, without markers, and not requiring time-consuming pre-registration calibration and post-registration labeling/tracking. Considering these assets, one may consider a multi-Kinect v2 set-up as a serious alternative for quantitative gait assessments.

A limitation of the multi-Kinect v2 set-up is the relatively low sampling frequency of 30 Hz. Although a good agreement between the multi-Kinect v2 set-up and the Optotrak system was found for almost all outcome measures of the current study, other outcome measures of interest may require higher sampling rates (e.g., the analysis of stride-to-stride fluctuations in stride times [31]). Another limitation of the study was that the between-systems agreement was only assessed for healthy subjects. Before implementing the multi-Kinect v2 walkway in the clinic, gait parameters for the patient groups of interest should be validated first. Moreover, one can imagine that in a clinical context an accompanying person such as a therapist wants to walk along with a patient for safety reasons. Because 3D positional data of body points of up to six persons can be tracked with a Kinect v2 sensor, each being allocated with a unique body identification number, it is important to ensure the correct allocation of data to a specific person when tracking multiple persons with multiple Kinects (e.g., using minimization of 3D positional data when moving from one camera’s field of view to another). Therefore, gait parameters need to be validated in various patient groups both with and without an accompanying person. As in healthy controls, good human-pose estimation is to be expected for patients. Clark et al. [32], for example, recently concluded that gait parameters of stroke patients derived from Kinect v1 data were highly reliable and could provide valuable additional information for gait analysis alongside the 10WMT. They stated that their findings provide support for implementing Kinect-based gait assessments in clinical settings [32]. With the development and validation of the multi-Kinect v2 instrumented 10-meter walkway, the current study may help pave the way to fulfill that premise.

Conclusion

Body point’s time series obtained with a multi-Kinect v2 set-up match well with those derived with a gold standard in 3D measurement accuracy, particularly so for body points in motion. The excellent absolute agreements with the gold standard observed for time to walk 10 meters, walking speed, cadence, step length, stride length, step time and stride time emphasize that those parameters can be reliably obtained with the multi-Kinect v2 set-up. Future studies are recommended to test the clinical utility of the multi-Kinect v2 set-up to automate 10MWT assessments, thereby complementing the time to walk 10 meters with reliable spatiotemporal gait parameters obtained objectively in a quick, unobtrusive and patient-friendly manner.

Supporting Information

S1 Data. Body point’s time series in the AP, ML and V direction for the multi-Kinect v2 set-up and the Optotrak system.

https://doi.org/10.1371/journal.pone.0139913.s001

(ZIP)

S1 Fig. Bland-Altman plots for the spatiotemporal gait parameters for CWS and MWS conditions.

https://doi.org/10.1371/journal.pone.0139913.s002

(PDF)

S1 Table. Overview of the body points obtained with the multi-Kinect v2 set-up and the Optotrak system.

https://doi.org/10.1371/journal.pone.0139913.s003

(PDF)

S1 Video. Body point’s time series obtained with the multi-Kinect v2 set-up and the Optotrak system of a single representative trial during the CWS condition of the 10MWT.

https://doi.org/10.1371/journal.pone.0139913.s004

(MP4)

Acknowledgments

We are grateful to Microsoft for accepting MR and HC as members of the Kinect for Windows developers program, allowing us to work with the Kinect v2 sensor and SDK before its commercial release. We would like to acknowledge Aukje de Vrijer from NDI for her invaluable help with setting up the Optotrak system and introducing us to the First Principles software.

Author Contributions

Conceived and designed the experiments: DJG BHC MR. Performed the experiments: DJG. Analyzed the data: DJG MR. Contributed reagents/materials/analysis tools: BHC MR. Wrote the paper: DJG MR. Acquisition of funding: MR.

References

  1. 1. Clemson L, Kendig H, Mackenzie L, Browning C. Predictors of injurious falls and fear of falling differ: an 11-year longitudinal study of incident events in older people. J Aging Health. 2015;27(2):239–256. pmid:25117181
  2. 2. Quach L, Galica AM, Jones RN, Procter-Gray E, Manor B, Hannan MT, et al. The nonlinear relationship between gait speed and falls: the maintenance of balance, independent living, intellect, and zest in the elderly of Boston study. J Am Geriatr Soc. 2011;59(6):1069–1073. pmid:21649615
  3. 3. Verghese J, Holtzer R, Lipton RB, Wang C. Quantitative gait markers and incident fall risk in older adults. J Gerontol. 2009;64A(8):896–901.
  4. 4. Montero-Odasso M, Schapira M, Soriano ER, Varela M, Kaplan R, Camera LA, et al. Gait velocity as a single predictor of adverse events in healthy seniors aged 75 years and older. J Gerontol. 2005;60A(10):1304–1309.
  5. 5. Van Kan GA, Rolland Y, Andrieu S, Bauer J, Beauchet O, Bonnefoy M, et al. Gait speed at usual pace as a predictor of adverse outcomes in community-dwelling older people an International Academy on Nutrition and Aging (IANA) Task Force. J Nutr Health Aging. 2009;13(10):881–889. pmid:19924348
  6. 6. Studenski S, Perera S, Patel K, Rosano C, Faulkner K, Inzitari M, et al. Gait speed and survival in older adults. JAMA. 2011;305(1):50–58. pmid:21205966
  7. 7. Simon SR. Quantification of human motion: gait analysis—benefits and limitations to its application to clinical problems. J Biomech. 2004;37(12):1869–1880. pmid:15519595
  8. 8. Shotton J, Fitzgibbon A, Cook M, Sharp T, Finocchio M, Moore R, et al. Real-time human pose recognition in parts from single depth images. Proceeding of the IEEE Conference on Computer Vision and Pattern Recognition; 2011 June 20–25; Providence, United States. IEEE; 2011. p. 1297–1304.
  9. 9. Baldewijns G, Verheyden G, Vanrumste B, Croonenborghs T. Validation of the Kinect for gait analysis using the GAITRite walkway. Proceeding of the 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society; 2014 Aug 26–30; Chicago, United States. IEEE; 2014. p. 5920–5923.
  10. 10. Clark RA, Pua Y, Fortin K, Ritchie C, Webster KE, Denehy L, et al. Validity of the Microsoft Kinect for assessment of postural control. Gait Posture. 2012;36(3):372–377. pmid:22633015
  11. 11. Clark RA, Bower KJ, Mentiplay BF, Paterson K, Pua Y. Concurrent validity of the Microsoft Kinect for assessment of spatiotemporal gait variables. J Biomech. 2013;46(15):2722–2725. pmid:24016679
  12. 12. Clark RA, Pua YH, Oliveira CC, Bower KJ, Thilarajah S, McGaw R, Hasanki K, Mentiplay BF. Reliability and concurrent validity of the Microsoft Xbox One Kinect for assessment of standing balance and postural control. Gait Posture. 2015;42(2):210–213. pmid:26009500
  13. 13. Dubois A, Charpillet F. A gait analysis method based on a depth camera for fall prevention. Proceeding of the 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society; 2014 Aug 26–30; Chicago, United States. IEEE; 2014. p. 4515–4518.
  14. 14. Gabel M, Gilad-Bachrach R, Renshaw E, Schuster A. Full-body Gait Analysis with Kinect. Proceeding of the Annual International Conference of the IEEE on Engineering in Medicine and Biology Society; 2012 Aug 28 –Sep 1; San Diego, United States. IEEE; 2012. p. 1964–1967.
  15. 15. Mentiplay BF, Perraton LG, Bower KJ, Pua YH, McGaw R, Heywood S, Clark RA. Gait assessment using the Microsoft Xbox One Kinect: Concurrent validity and inter-day reliability of spatiotemporal and kinematic variables. J Biomech. 2015;48(10):2166–2170. pmid:26065332
  16. 16. Olesh EV, Yakovenko S, Gritsenko V. Automated assessment of upper extremity movement impairment due to stroke. PLOS ONE. 2014;9(8):e104487. pmid:25100036
  17. 17. Stone EE, Skubic M. Evaluation of an Inexpensive Depth Camera for In-Home Gait Assessment. Proceeding of the 5th International Conference on Pervasive Computing Technologies for Healthcare; 2011 May 23–26; Dublin, Ireland. IEEE; 2011. p. 71–77.
  18. 18. Sun T, Lee C. An impact study of the design of exergaming parameters on body intensity from objective and gameplay-based player experience perspectives, based on balance training exergame. PLOS ONE. 2013;8(7):e69471. pmid:23922716
  19. 19. Kinect for Windows features [Internet]. Cited 23 March 2015. Available from: http://www.microsoft.com/en-us/kinectforwindows/meetkinect/features.aspx.
  20. 20. Graham JE, Ostir GV, Fisher SR, Ottenbacher KJ. Assessing walking speed in clinical research: a systematic review. J Eval Clin Pract. 2008;14(4):552–562. pmid:18462283
  21. 21. Salbach NM, O'Brien KK, Brooks D, Irvin E, Martino R, Takhar P, et al. Reference values for standardized tests of walking speed and distance: a systematic review. Gait Posture. 2015;41(2):341–360. pmid:25542397
  22. 22. Zeni JA, Richards JG, Higginson JS. Two simple methods for determining gait events during treadmill and overground walking using kinematic data. Gait Posture. 2008;27(4):710–714. pmid:17723303
  23. 23. Bland JM, Altman DG. A note on the use of the intraclass correlation coefficient in the evaluation of agreement between two methods of measurement. Comput Biol Med. 1990;20(5):337–340. pmid:2257734
  24. 24. McGraw KO, Wong SP. Forming inferences about some intraclass correlation coefficients. Psychol Methods. 1996;1(1):30–46.
  25. 25. Cicchetti DV. Guidelines, criteria, and rules of the thumb for evaluating normed and standardized assessment instruments in psychology. Psychol Assess. 1994;6(4):284–290.
  26. 26. Bland JM, Altman DG. Measuring agreement in method comparison studies. Stat Methods Med Res. 1999;8(2):135–160. pmid:10501650
  27. 27. Flansbjer UB, Holmbäck AM, Downham D, Patten C, Lexell J. Reliability of gait performance tests in men and women with hemiparesis after stroke. J Rehabil Med. 2005;37(2):75–82. pmid:15788341
  28. 28. Galna B, Barry G, Jackson D, Mhiripiri D, Olivier P, Rochester L. Accuracy of the Microsoft Kinect sensor for measuring movement in people with Parkinson's disease. Gait Posture. 2014;39(4):1062–1068. pmid:24560691
  29. 29. Kitsunezaki N, Adachi E, Masuda T, Mizusawa J. KINECT Applications for The Physical Rehabilitation. Proceeding of the IEEE International Symposium on Medical Measurement and Applications; 2013 May 4–5; Gatineau, Canada. IEEE; 2013. p. 294–299.
  30. 30. Perera S, Mody SH, Woodman RC, Studenski SA. Meaningful change and responsiveness in common physical performance measures in older adults. J Am Geriatr Soc. 2006;54(5):743–749. pmid:16696738
  31. 31. Roerdink M, Daffertshofer A, Marmelat V, Beek PJ. How to sync to the beat of a persistent fractal metronome without falling off the treadmill? PLoS ONE. 2015;10(7):e0134148. pmid:26230254
  32. 32. Clark RA, Vernon S, Mentiplay BF, Miller KJ, McGinley JL, Pua YH, et al. Instrumenting gait assessment using the Kinect in people living with stroke: reliability and association with balance tests. J Neuroeng Rehabil. 2015;12(1):15.