Next Article in Journal
A Sequential Color Correction Approach for Texture Mapping of 3D Meshes
Next Article in Special Issue
Effects of Motor Task Difficulty on Postural Control Complexity during Dual Tasks in Young Adults: A Nonlinear Approach
Previous Article in Journal
Remote Photoplethysmography with a High-Speed Camera Reveals Temporal and Amplitude Differences between Glabrous and Non-Glabrous Skin
Previous Article in Special Issue
Accuracy of Inertial Measurement Units When Applied to the Countermovement Jump of Track and Field Athletes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Gait Feature Assessment on Two-Dimensional Body Axis Projection Planes Converted from Three-Dimensional Coordinates Estimated with a Deep Learning Smartphone App

1
Department of Neurosurgery, Nagoya City University Graduate School of Medical Science, Nagoya 467-8601, Japan
2
Normal Pressure Hydrocephalus Center, Rakuwakai Otowa Hospital, Kyoto 607-8062, Japan
3
Interfaculty Initiative in Information Studies/Institute of Industrial Science, The University of Tokyo, Tokyo 113-8654, Japan
4
Digital Standard Co., Ltd., Osaka 536-0013, Japan
5
Division of Neurology and Clinical Neuroscience, Department of Internal Medicine III, Yamagata University School of Medicine, Yamagata 990-9585, Japan
6
Human Augmentation Research Center, National Institute of Advanced Industrial Science and Technology (AIST), Kashiwa II Campus, University of Tokyo, Kashiwa 277-0882, Japan
7
Shin-Aikai Spine Center, Katano Hospital, Katano 576-0043, Japan
8
School of Medicine, Shiga University of Medical Science, Otsu 520-2192, Japan
9
Department of Informatics, Faculty of Engineering, Yamagata University, Yamagata 992-8510, Japan
10
Rakuwa Villa Ilios, Rakuwakai Healthcare System, Kyoto 607-8062, Japan
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(2), 617; https://doi.org/10.3390/s23020617
Submission received: 5 December 2022 / Revised: 29 December 2022 / Accepted: 4 January 2023 / Published: 5 January 2023

Abstract

:
To assess pathological gaits quantitatively, three-dimensional coordinates estimated with a deep learning model were converted into body axis plane projections. First, 15 healthy volunteers performed four gait patterns; that is, normal, shuffling, short-stepped, and wide-based gaits, with the Three-Dimensional Pose Tracker for Gait Test (TDPT-GT) application. Second, gaits of 47 patients with idiopathic normal pressure hydrocephalus (iNPH) and 92 healthy elderly individuals in the Takahata cohort were assessed with the TDPT-GT. Two-dimensional relative coordinates were calculated from the three-dimensional coordinates by projecting the sagittal, coronal, and axial planes. Indices of the two-dimensional relative coordinates associated with a pathological gait were comprehensively explored. The candidate indices for the shuffling gait were the angle range of the hip joint < 30° and relative vertical amplitude of the heel < 0.1 on the sagittal projection plane. For the short-stepped gait, the angle range of the knee joint < 45° on the sagittal projection plane was a candidate index. The candidate index for the wide-based gait was the leg outward shift > 0.1 on the axial projection plane. In conclusion, the two-dimensional coordinates on the body axis projection planes calculated from the 3D relative coordinates estimated by the TDPT-GT application enabled the quantification of pathological gait features.

Graphical Abstract

1. Introduction

A pathological gait is a gait that deviates from a normal gait, and its features are characterized as shuffling, short-stepped, freezing, wide-based, festination, hemiplegic, spastic, ataxic, and instability. [1,2,3,4,5,6]. Although videorecording has been clinically recommended when evaluating pathological gaits [2,3,4,5], these gait features have been evaluated subjectively without any standardized rating systems. The qualitative or semi-quantitative subjective assessments are simple and widely used in clinical practice, including medical record input, but their validity and reliability are generally low [1,2]. Therefore, gait analysis studies to assess movement disorders require quantitative values or indices using some type of instrumentation; for example, a reduced stride length and diminished step height were reported to be typical spatiotemporal and kinematic characteristics of a pathological gait in Parkinson’s disease and idiopathic normal pressure hydrocephalus (iNPH) [6,7,8,9,10,11,12], asymmetry of step length, and leg swing velocity were typical for a hemiplegic gait [13,14], and coefficients of variation of stride length and step width were typical for an ataxic gait and instability [15]. The three-dimensional (3D) acceleration of the trunk has been reported to be a useful parameter for detecting pathological gait features. In our previous study, however, 3D trunk acceleration during straight walking could not clearly distinguish the six gait features of shuffling, short-stepped, freezing, wide-based, festination, and instability in patients with iNPH [16]. Therefore, we reaffirmed the need to analyze the spatiotemporal and kinematic parameters of gait; for instance, the joint angle, stride length, and step width using a motion capture technology. An optical 3D motion capture system with multipoint cameras and numerous body markers is a precise method for gait analysis and has become the gold standard [17,18,19,20]. However, it has not been clinically applied because it takes time and effort to train technicians and carry out measurements and analysis, and the equipment is expensive. Therefore, we developed a novel iOS application named Three-Dimensional Pose Tracker for Gait Test (TDPT-GT) (https://digital-standard.com/tdptios_en/, accessed on 1 December 2022) in which the 3D motion of the full body was simply tracked from two-dimensional (2D) video images captured by a single iPhone camera [21]. This application can estimate the 3D relative coordinates of 24 key points in a full body-pose during movement without time-consuming complicated preparation. Furthermore, we developed a new method for assessing pathological gaits more easily by converting 3D relative coordinates into 2D coordinates projected on the sagittal, coronal, and axial planes with respect to the body axis. In this study, we explored the novel distinctive indices that express the essence of each pathological gait feature from 2D relative coordinates projected on the sagittal and axial planes.

2. Materials and Methods

2.1. Ethical Approvals

The study design and protocol of this study were approved by the ethics committee for human research at our institute (IRB number: R2019-337). All volunteers and patients participated in this study after providing written informed consent. The study design was prospective and observational. This study was conducted according to the approved guidelines of the Declaration of Helsinki.

2.2. Study Population

From November 2020 to August 2022, 15 healthy medical students, medical staff, and their families, aged 20 years or older, were recruited by open recruitment. The inclusion criteria for this study were independence in their daily lives, without any problems in walking, normal cognitive function, and no history of neurodegenerative disorders, severe head injury, or cerebrovascular disease. The volunteers walked normally and acted out pathological gait features, such as a short-stepped gait, short-stepped shuffling wide-based gait, and short-stepped shuffling closed gait. A short-stepped gait was defined as a reduced stride length, and a shuffling gait had a diminished step height. In addition, we acquired two additional datasets to assess the generalizability of the candidate parameters and indices. The first dataset contained data from 47 patients diagnosed with iNPH, and the second dataset contained data from 92 healthy volunteers aged 60 years or older who participated in the Takahata population-based cohort study [22,23,24]. The inclusion criteria for all participants, including iNPH patients, were independence in their ability to walk without any assistance and the ability to make independent decisions. These patients participated in the study on their own initiative. The patients with iNPH were subjectively evaluated for the presence or absence of the following six pathological gait features: a shuffling gait, a short-stepped gait, a wide-based gait, a freezing gait, a spastic gait, and instability [16]. The simulated pathological gaits of the healthy volunteers were quite extreme, compared with those of the patients with iNPH.

2.3. Data Acquisition of Estimated Three-Dimensional Relative Coordinates

The walking movements of participants while walking in a circle with a diameter of 1 m for two laps, clockwise and counterclockwise, were recorded from head to toe in the videoframes of the TDPT-GT application. The TDPT-GT application worked precisely under the following conditions. The lines of the subject’s body and legs were as easy to identify as possible. The recording took place in a bright space with a uniform background where no mirror was present and no one other than the subject appeared in the video frame. Details of the technology of the TDPT-GT application and data acquisition were described in our prior publication [21]. In brief, the TDPT-GT application is an artificial intelligence (AI)-based motion capture system using a 3D heatmap estimation that could measure the 3D relative coordinates of the human body at 30 fps with 448 × 448 pixels and RGB color using an iPhone camera without any markers for motion capture. The TDPT-GT application estimates the 3D relative coordinates of the following 24 key points: the center of the body (navel), nose, left and right ears, eyes, shoulders, elbows, wrists, thumbs, middle fingers, hips, knees, heels, and toes, based on the center of the body, not the absolute coordinates that defined the ground. The x- and y-axis coordinates on the screen were converted directly into the three-dimensional x- and y-axis coordinates, and the center point of depth on the screen was always set at the subject’s navel. The raw and smoothed coordinates were applied to the low-pass filter, and the AI scores, which represented the certainty and probability of coordinate estimations, were automatically calculated and saved as a csv file format in the iPhone. The head center was defined as the center coordinates of both ears, and the neck center was defined as the center coordinates of both shoulders. For relative 3D coordinates, the length of the upper body was calculated as the length from the navel center to the neck center and the length from the neck center to the head center, and this measurement was set as 1. The leg length was calculated as the total length from the center to the hip joint, from the hip joint to the knee, and from the knee to the heel, and this measurement was set as 1. The reliability of the 3D relative coordinates estimated by the TDPT-GT application was evaluated and compared with 3D optical motion capture systems, described in a previous study [21].

2.4. Projection of Relative Coordinates on the Body Axis Planes

As shown in Figure 1, 2D relative coordinates of the 24 key points were calculated from 3D relative coordinates with an AI score of 0.7 or higher and projected onto the sagittal, coronal (frontal), and axial planes based on the body axis. The normal vector to a plane passing through the navel center and right and left shoulder points was defined as U F and the normal vector to a plane passing through the navel center and right and left hip joints was defined as D F . The f for the forward direction of the body axis was calculated as follows: ( U F + D F ) / 2 . A sagittal plane was defined as a plane composed of f and c from the navel center to the center between the bilateral hip joints. A coronal (frontal) plane was defined as a plane composed of the c and n that was orthogonal to the two vectors of f and c . An axial plane was defined as a plane composed of f and n ; that is, f · c = f · n = n · c = 0. Each 2D coordinate projected on the sagittal, coronal, and axial planes of each 3D relative coordinate was the intersection of the projected plane and the normal line with respect to the projected plane. In addition, the 2D relative coordinates of the bilateral hip joints, knees, heels, and toes were corrected by the total leg length as 1 on the projected plane, and those of the upper body points were also corrected by the total length of the upper body as 1.
Estimated three-dimensional coordinates of 24 full body key points were converted to the two-dimensional relative coordinates by the projection onto the sagittal, coronal, and axial planes based on the body axis. First, the forward vector f relative to the body axis from the center and both shoulder points (light brown) and hip joints (light salmon) was determined. Next, the downward vector c relative to the body axis from the center and both hip joints was determined.

2.5. Data Processing for Gait Analysis

The chronological changes in the 2D relative coordinates excluding the upper limbs were plotted on the sagittal, coronal (frontal), and axial planes, as shown in Figure 2. Additionally, 75% tolerance ellipses of all plots for each key point were drawn on the projection planes. The elliptical area, center coordinate of the ellipses, tilt angles of major and minor axes, and the x-axis and y-axis lengths of the 75% tolerance ellipses were calculated at each 2D relative coordinate on the projection planes as follow: dataEllipse (x coordinates of key point, y coordinates of key point, levels = 0.75). The center coordinate of the ellipse was calculated as cov.wt (cbind (x coordinates of key point, y coordinates of key point)), and the elliptical area was calculated as the length of the major axis of the ellipse × the length of the minor axis of the ellipse × π (pi).
On the sagittal projection plane (Figure 2, lower left), an angle range of the hip joint was defined as an angle between two vectors from the center coordinate of the 75% tolerance ellipse for the hip joint to two coordinates of the major axis for the knee ellipse, and an angle range of the knee joint was defined as an angle between two vectors from the knee ellipse center to two coordinates of the major axis for the heel ellipse. A relative vertical amplitude of the heel was defined as the y-axis width of the heel plots.
On the axial projection plane (Figure 2, lower right), a foot angle [6], i.e., outward rotation of the feet, was defined as an angle between two vectors from the right and left heel ellipse centers to the right and left toe ellipse centers, respectively. A toe outward shift was defined as a lateral deviation of the right and left x-coordinates of toe ellipse centers from those of the heel ellipse centers, a heel outward shift was defined as a lateral deviation of the right and left x-coordinates of the heel ellipse centers from those of the hip ellipse centers, and a leg outward shift was defined as the sum of the toe and heel outward shifts. The indices of the 2D relative coordinates projected onto the coronal (frontal) plane were not investigated in this study because the bending direction of the knee joint or ankle joint was sagittal, and the foot angle or outward shift was easier to assess on the axial plane rather than on the coronal plane.

2.6. Statistical Analysis

After detecting appropriate indices for each pathological gait feature by the first comprehensive parameter search in the gait dataset of 15 healthy volunteers, optimal cutoff values of the candidate indices were investigated based on the distribution of the measurements and indices in the iNPH patients with gait disturbances and healthy individuals aged 60 years or older in the Takahata cohort study. In addition, the detectability of the shuffling gait, short-stepped gait, and wide-based gait at the setting cutoff values was evaluated with the area under the receiver-operating characteristic curve (AUC), sensitivity, specificity, and odds ratio. The statistical significance was assumed at a probability (p) value of <0.05. All missing data points were treated as deficit data that did not affect other variables. Statistical analyses were performed using R (version 4.1.1; The R Foundation for Statistical Computing; http://www.R-project.org accessed on 1 December 2022).

3. Results

3.1. Clinical Characteristics

Fifteen healthy volunteers (mean age 39.1 ± 20.1 years; range, 22–78 years; nine males and six females; height 148–194 cm) were assessed for a normal gait and the acting pathological gait features as follows: short-stepped shuffling wide-based gait, short-stepped shuffling closed gait, and short-stepped gait, according to the TDPT-GT application. As the second step, the gait of 47 patients diagnosed with iNPH (mean age 77.3 ± 6.3 years; range, 61–87 years; 32 males and 15 females) and 92 healthy elderly individuals (mean age 73.0 ± 6.3 years; range, 60–91 years; 36 males and 56 females) in the Takahata cohort study were assessed using the TDPT-GT. Table 1 summarizes the participants’ characteristics.

3.2. Two-Dimensional Relative Coordinates on the Sagittal and Axial Projection Planes

The movement trajectories (plot distribution) and 75% tolerance ellipse parameters for the 2D relative coordinates projected onto the sagittal and axial planes were compared between the normal gaits and acting pathological gaits in the first 15 healthy volunteer group (Figure 3).
Compared to the normal gait, the range of the back-and-forth leg swing was apparently smaller in the short-stepped shuffling gait, regardless of whether the legs were opened or closed. In the short-stepped but non-shuffling gait (Figure 3b), the range of the back-and-forth swing of the lower legs, i.e., the heels and toes, was also smaller because of the smaller stride length, but the knees moved markedly in an upward and forward direction.

3.2.1. Sagittal Projection Plane

The shuffling gait could be clearly distinguished by the angle range of the hip joint, angle range of the knee joint, and relative vertical amplitude of the heel (Figure 4, Figure 5 and Figure 6), whereas the short-stepped gait could be clearly distinguished only by the angle range of the knee joint (Figure 5). As the optimal cutoff values, <30° of the angle range of the hip joint for detecting a shuffling gait (Figure 4) and <45° of the angle range of the knee joint for detecting a shuffling gait and short-stepped gait (Figure 5) were proposed. As the optimal cutoff value for detecting the shuffling gait, we proposed <0.1 of the relative vertical amplitude of the heel with the total leg length set as 1 on the sagittal projection plane (Figure 6).

3.2.2. Axial Projection Plane

On the axial projection plane, the foot angle between the right and left vectors from the heel to toes could not distinguish the wide-based and normal gaits (Figure 7).
The wide-based gait could be clearly distinguished by the heel outward shift, toe outward shift, and leg outward shift, which was the sum of the heel and toe outward shift (Figure 8, Figure 9 and Figure 10). As the optimal cutoff values for a wide-based gait, we proposed >−0.08 of the heel outward shift (Figure 8), >0.18 of the toe outward shift (Figure 9), and >0.1 of the leg outward shift, with the total leg length set as 1 on the axial projection plane (Figure 10).

3.3. Discrimination of Pathological Gaits by Indices on 2D Projection Planes Using the TDPT-GT App

Although we subjectively evaluated shuffling and short-stepped and wide-based gaits in 47 patients with iNPH, as shown in Table 1, these gait types were not used for investigating indices of 2D relative coordinates projected onto the sagittal or axial planes of the body axis, because the validity of subjective assessments of gait features is known to be low [2] and the shuffling, short-stepped, and wide-based gait are closely related to each other [16]. Instead, the AUC (95% CIs), sensitivity, specificity, and odds ratio (95% CIs) at the setting cutoff values were evaluated (Table 2). Of them, <45° of the angle range of the knee joint for detecting a short-stepped gait had the highest AUC (81.7, 95% CIs = 77.4–86.0) and the highest odds ratio (10.55, 95% CIs = 5.63–21.12).

4. Discussion

To establish generalized indices for pathological gaits that can be commonly used, it is necessary to first extract the feature quantities of pathological gaits using a simple method and then compare the distribution of them between healthy and affected subjects to determine the optimal cutoff values. Prior to this study, we developed a novel markerless monocular motion capture smartphone application that could obtain the moving 3D relative coordinates of the entire human body based on simple 2D video recordings [21]. In this study, we suggested that the angle range of the hip joint be <30° and the relative vertical amplitude of the heel < 0.1 as novel indices for detecting a shuffling gait, and the angle range of the knee joint < 45° as a useful index for detecting a shuffling gait and short-stepped gait on the sagittal projection plane calculated from the 3D relative coordinates, estimated by the TDPT-GT application. In addition, to detect a wide-based gait, a leg outward shift > 0.1 on the axial projection plane was proposed as a useful index rather than the foot angle. The most characteristic gait feature specific to iNPH to differentiate it from Parkinson’s disease was suggested to be an opening foot angle, i.e., an outward rotation of the feet [6]. The reason why our result was inconsistent with previous reports could be that our task included walking around a circle with a diameter of 1 m in order to hold the smartphone close to the examinee. Using the TDPT-GT application with our method, it is possible to assess gait in a confined space, such as a private clinic or one’s own room. However, toes are likely to be open with a normal gait and the foot angle is not a useful index for detecting a wide-based gait with our measurement method.
In general, a person’s normal gait changes with age. Compared to younger adults, older adults walk slowly with a shorter and wider step due to decreased forward propulsion [19,20,25,26,27]. Therefore, we needed to evaluate the gait of patients with iNPH and age-matched controls using the same method to determine candidate indices and their cutoff values for detecting pathological gait features. Furthermore, we needed to evaluate simulated pathological gaits in volunteers, separate from the case-controlled group, to ensure the identification of the pathological gait features of the shuffling, short-stepped, and wide-based gaits because these gait patterns are closely related to each other and could not be clearly distinguished in the real patients with iNPH and older adults. For example, a wide-based gait reduces the stride length, resulting in a short-stepped gait and a shuffling gait. Conversely, a shuffling gait with large steps or a short-stepped gait with rising legs is not a normal situation when walking.
The greatest advantage of this study is that any clinician or researcher can easily evaluate a pathological gait using a simple smartphone application. Optical 3D motion capture systems with multipoint cameras and markers for motion capture are prohibitively expensive and may be unsuitable for general use. In contrast, our proposed gait assessment can be applied easily without considerable amounts of time and money. Furthermore, in our system, large amounts of human motion data can be processed within a few milliseconds. This enables the application to perform in real time to support capturing the 3D pose, analyzing gait, and encouraging physical exercise. Consequently, these measurements can be used in multicenter collaborative studies to assess the severity of the pathological gait and changes in gait disturbances using a monocular smartphone camera.
This study has some limitations that warrant discussion. First, the 2D relative coordinates projected onto the sagittal or axial planes of the body axis based on the navel as the center without the ground information could not be compared with the previous conventional gait research, which was usually assessed by separating the stance phase from the swing phase or measuring the distance from the floor surface. However, the gait cycle time, stride length, step frequency, and angle range of the joint can be estimated easily using our method. Second, this system could track the body pose of one person only, and the iPhone camera must always show the entire body. Therefore, it may not be usable in a cluttered environment, such as a group health check. In order to implement our simple gait assessment in society using the TDPT-GT application, it is necessary to accumulate experience in more collaborate studies and clarify its advantages and disadvantages in our proposed gait analysis without specifying the floor surface. The development of objective measurements would greatly help with assessing the efficacy of therapeutic interventions, including medications, surgery, and the rehabilitation for any movement disorders.

5. Conclusions

The 2D coordinates projected on the sagittal, coronal, and axial planes with respect to the body axis were calculated from the 3D relative coordinates estimated by our AI-based iPhone application. Using the TDPT-GT application, anyone can easily evaluate various pathological gait patterns quantitatively, which previously required a large-scale 3D motion analysis system. The candidate indices for the shuffling gait were the angle range of the hip joint < 30°, the angle range of the knee joint < 45°, and relative vertical amplitude of the heel < 0.1 on the sagittal projection plane. The candidate index for the short-stepped gait was the angle range of the knee joint < 45° on the sagittal projection plane. The candidate index for the wide-based gait was not the foot angle but the leg outward shift > 0.1 on the axial projection plane. In the future, this AI-based application will quantify pathological gaits, making it possible for non-specialist individuals to easily detect fall risks and gait disorders at an early stage, which will help to prevent worse conditions that would require nursing care.

Author Contributions

Conceptualization, S.Y. and Y.A.; methodology, S.Y. and Y.A.; software, Y.A.; validation, S.Y., Y.A. and Y.K.; formal analysis, S.Y., Y.A., T.K. and C.I.; investigation, S.Y., S.U., T.K. and C.I.; resources, Y.A.; data curation, S.Y., Y.A., T.K., K.M. and T.F.; writing (original draft preparation), S.Y.; writing (review and editing), S.Y.; visualization, S.Y. and Y.A.; supervision, M.T., M.H., M.I., Y.O. and M.M.; project administration, S.Y.; funding acquisition, S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by research grants from the Japan Society for the Promotion of Science, KAKENHI for 3 years, beginning in 2021 (grant number: 21K09098); from the G-7 Scholarship Foundation in 2020; from the Taiju Life Social Welfare Foundation in 2020 and 2022; and from the Osaka Gas Group Welfare Foundation in 2022. The funders had no effect or involvement in the writing of this article. The funders played no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

The study was conducted according to the Declaration of Helsinki and approved by the Ethics Committee for Human Research of Shiga University of Medical Science (protocol code: R2019-337 and date of approval: 17 April 2020).

Informed Consent Statement

Written informed consent was obtained from all subjects involved in this study.

Data Availability Statement

Data generated or analyzed during the study are available from the corresponding author upon reasonable request.

Acknowledgments

We would like to thank the rehabilitation staff at the Shiga University of Medical Science, Yamagata University Hospital, Katano Hospital, Tesseikai Neurosurgery Hospital, and Rakuwakai Otowa Hospital for their help in data acquisition. We thank Genjiro Katayama, Miki Yoshioka, and Digital Standard Co., Ltd. for allowing the use to use their application products.

Conflicts of Interest

Yukihiko Aoyagi was employed by Digital Standard Co., Ltd. The remaining authors declare that this research was conducted in the absence of any commercial or financial relationships that could be interpreted as a potential conflict of interest.

References

  1. Gor-García-Fogeda, M.D.; Cano de la Cuerda, R.; Carratalá Tejada, M.; Alguacil-Diego, I.M.; Molina-Rueda, F. Observational gait assessments in people with neurological disorders: A systematic review. Arch. Phys. Med. Rehabil. 2016, 97, 131–140. [Google Scholar] [CrossRef]
  2. Ishikawa, M.; Yamada, S.; Yamamoto, K. Agreement study on gait assessment using a video-assisted rating method in patients with idiopathic normal-pressure hydrocephalus. PLoS ONE 2019, 14, e0224202. [Google Scholar] [CrossRef]
  3. Marmarou, A.; Bergsneider, M.; Relkin, N.; Klinge, P.; Black, P.M. Development of guidelines for idiopathic normal-pressure hydrocephalus: Introduction. Neurosurgery 2005, 57, S1–S3. [Google Scholar] [CrossRef] [Green Version]
  4. Nakajima, M.; Yamada, S.; Miyajima, M.; Ishii, K.; Kuriyama, N.; Kazui, H.; Kanemoto, H.; Suehiro, T.; Yoshiyama, K.; Kameda, M.; et al. Guidelines for Management of Idiopathic Normal Pressure Hydrocephalus (Third Edition): Endorsed by the Japanese Society of Normal Pressure Hydrocephalus. Neurol. Med.-Chir. 2021, 61, 63–97. [Google Scholar] [CrossRef] [PubMed]
  5. Scully, A.E.; Lim, E.C.W.; Teow, P.P.; Tan, D.M.L. A systematic review of the diagnostic utility of simple tests of change after trial removal of cerebrospinal fluid in adults with normal pressure hydrocephalus. Clin. Rehabil. 2018, 32, 942–953. [Google Scholar] [CrossRef] [PubMed]
  6. Stolze, H.; Kuhtz-Buschbeck, J.P.; Drücke, H.; Jöhnk, K.; Diercks, C.; Palmié, S.; Mehdorn, H.M.; Illert, M.; Deuschl, G. Gait analysis in idiopathic normal pressure hydrocephalus--Which parameters respond to the CSF tap test? Clin. Neurophysiol. 2000, 111, 1678–1686. [Google Scholar] [CrossRef] [PubMed]
  7. Hausdorff, J.M. Gait dynamics, fractals and falls: Finding meaning in the stride-to-stride fluctuations of human walking. Hum. Mov. Sci. 2007, 26, 555–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Kitade, I.; Kitai, R.; Neishi, H.; Kikuta, K.I.; Shimada, S.; Matsumine, A. Relationship between gait parameters and MR imaging in idiopathic normal pressure hydrocephalus patients after shunt surgery. Gait Posture 2018, 61, 163–168. [Google Scholar] [CrossRef]
  9. Murray, M.P.; Sepic, S.B.; Gardner, G.M.; Downs, W.J. Walking patterns of men with parkinsonism. Am. J. Phys. Med. 1978, 57, 278–294. [Google Scholar]
  10. Panciani, P.P.; Migliorati, K.; Muratori, A.; Gelmini, M.; Padovani, A.; Fontanella, M. Computerized gait analysis with inertial sensor in the management of idiopathic normal pressure hydrocephalus. Eur. J. Phys. Rehabil. Med. 2018, 54, 724–729. [Google Scholar] [CrossRef]
  11. Williams, M.A.; Thomas, G.; de Lateur, B.; Imteyaz, H.; Rose, J.G.; Shore, W.S.; Kharkar, S.; Rigamonti, D. Objective assessment of gait in normal-pressure hydrocephalus. Am. J. Phys. Med. Rehabil. 2008, 87, 39–45. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, L.; Li, Y.; Xiong, F.; Zhang, W. Gait recognition using optical motion capture: A decision fusion based method. Sensors 2021, 21, 3496. [Google Scholar] [CrossRef] [PubMed]
  13. Hsu, A.L.; Tang, P.F.; Jan, M.H. Analysis of impairments influencing gait velocity and asymmetry of hemiplegic patients after mild to moderate stroke. Arch. Phys. Med. Rehabil. 2003, 84, 1185–1193. [Google Scholar] [CrossRef] [PubMed]
  14. Patterson, K.K.; Parafianowicz, I.; Danells, C.J.; Closson, V.; Verrier, M.C.; Staines, W.R.; Black, S.E.; McIlroy, W.E. Gait asymmetry in community-ambulating stroke survivors. Arch. Phys. Med. Rehabil. 2008, 89, 304–310. [Google Scholar] [CrossRef]
  15. Serrao, M.; Pierelli, F.; Ranavolo, A.; Draicchio, F.; Conte, C.; Don, R.; Di Fabio, R.; LeRose, M.; Padua, L.; Sandrini, G.; et al. Gait pattern in inherited cerebellar ataxias. Cerebellum 2012, 11, 194–211. [Google Scholar] [CrossRef]
  16. Yamada, S.; Aoyagi, Y.; Ishikawa, M.; Yamaguchi, M.; Yamamoto, K.; Nozaki, K. Gait assessment using three-dimensional acceleration of the trunk in idiopathic normal pressure hydrocephalus. Front. Aging Neurosci. 2021, 13, 653964. [Google Scholar] [CrossRef]
  17. Tsuchida, W.; Kobayashi, Y.; Inoue, K.; Horie, M.; Yoshihara, K.; Ooie, T. Kinematic characteristics during gait in frail older women identified by principal component analysis. Sci. Rep. 2022, 12, 1676. [Google Scholar] [CrossRef]
  18. Kobayashi, Y.; Hobara, H.; Matsushita, S.; Mochimaru, M. Key joint kinematic characteristics of the gait of fallers identified by principal component analysis. J. Biomech. 2014, 47, 2424–2429. [Google Scholar] [CrossRef]
  19. Kobayashi, Y.; Hobara, H.; Heldoorn, T.A.; Kouchi, M.; Mochimaru, M. Age-independent and age-dependent sex differences in gait pattern determined by principal component analysis. Gait Posture 2016, 46, 11–17. [Google Scholar] [CrossRef]
  20. Kobayashi, Y.; Ogata, T. Association between the gait pattern characteristics of older people and their two-step test scores. BMC Geriatr 2018, 18, 101. [Google Scholar] [CrossRef]
  21. Aoyagi, Y.; Yamada, S.; Ueda, S.; Iseki, C.; Kondo, T.; Mori, K.; Kobayashi, Y.; Fukami, T.; Hoshimaru, M.; Ishikawa, M.; et al. Development of smartphone application for markerless three-dimensional motion capture based on deep learning model. Sensors 2022, 22, 5282. [Google Scholar] [CrossRef] [PubMed]
  22. Iseki, C.; Kawanami, T.; Nagasawa, H.; Wada, M.; Koyama, S.; Kikuchi, K.; Arawaka, S.; Kurita, K.; Daimon, M.; Mori, E.; et al. Asymptomatic ventriculomegaly with features of idiopathic normal pressure hydrocephalus on MRI (AVIM) in the elderly: A prospective study in a Japanese population. J. Neurol. Sci. 2009, 277, 54–57. [Google Scholar] [CrossRef]
  23. Iseki, C.; Takahashi, Y.; Wada, M.; Kawanami, T.; Adachi, M.; Kato, T. Incidence of idiopathic normal pressure hydrocephalus (iNPH): A 10-year follow-up study of a rural community in Japan. J. Neurol. Sci. 2014, 339, 108–112. [Google Scholar] [CrossRef] [PubMed]
  24. Iseki, C.; Takahashi, Y.; Adachi, M.; Igari, R.; Sato, H.; Koyama, S.; Ishizawa, K.; Ohta, Y.; Kato, T. Prevalence and development of idiopathic normal pressure hydrocephalus: A 16-year longitudinal study in Japan. Acta Neurol. Scand. 2022, 146, 680–689. [Google Scholar] [CrossRef]
  25. Arvin, M.; Mazaheri, M.; Hoozemans, M.J.M.; Pijnappels, M.; Burger, B.J.; Verschueren, S.M.P.; van Dieën, J.H. Effects of narrow base gait on mediolateral balance control in young and older adults. J. Biomech. 2016, 49, 1264–1267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Hernández, A.; Silder, A.; Heiderscheit, B.C.; Thelen, D.G. Effect of age on center of mass motion during human walking. Gait Posture 2009, 30, 217–222. [Google Scholar] [CrossRef] [Green Version]
  27. Hurt, C.P.; Rosenblatt, N.; Crenshaw, J.R.; Grabiner, M.D. Variation in trunk kinematics influences variation in step width during treadmill walking by older and younger adults. Gait Posture 2010, 31, 461–464. [Google Scholar] [CrossRef]
Figure 1. Projection onto the sagittal, coronal, and axial planes based on the body axis.
Figure 1. Projection onto the sagittal, coronal, and axial planes based on the body axis.
Sensors 23 00617 g001
Figure 2. Two-dimensional relative coordinates projected onto the sagittal, coronal, and axial planes during two rounds of a 1-m diameter circular walk of a representative healthy volunteer. The chronological changes of the head (brown), neck (sky blue), right hip joint (orange), left hip joint (gold), right knee (red), left knee (yellow), right heel (blue), left heel (dark green), right toe (light green), and left toe (pale green) were plotted. The coordinates between joints, such as the hip and knee, knee and heel, and heel and toe, are connected by colored lines in the sagittal and coronal projection planes. In addition, 75% tolerance ellipses of all plots for each key point and its major axes and y-axis widths were drawn with each colored line on the left side of plots for the head and neck and the bilateral side of plots for legs. On the axial projection plane, plots and 75% tolerance ellipses for both legs and the horizontal and front-back range of motion for each key point were drawn with each colored line on the left, right, and bottom. The black arrows indicate vectors from the elliptical centers of the right and left hip joints to those of the knees.
Figure 2. Two-dimensional relative coordinates projected onto the sagittal, coronal, and axial planes during two rounds of a 1-m diameter circular walk of a representative healthy volunteer. The chronological changes of the head (brown), neck (sky blue), right hip joint (orange), left hip joint (gold), right knee (red), left knee (yellow), right heel (blue), left heel (dark green), right toe (light green), and left toe (pale green) were plotted. The coordinates between joints, such as the hip and knee, knee and heel, and heel and toe, are connected by colored lines in the sagittal and coronal projection planes. In addition, 75% tolerance ellipses of all plots for each key point and its major axes and y-axis widths were drawn with each colored line on the left side of plots for the head and neck and the bilateral side of plots for legs. On the axial projection plane, plots and 75% tolerance ellipses for both legs and the horizontal and front-back range of motion for each key point were drawn with each colored line on the left, right, and bottom. The black arrows indicate vectors from the elliptical centers of the right and left hip joints to those of the knees.
Sensors 23 00617 g002
Figure 3. Two-dimensional relative coordinates projected onto the sagittal and axial planes of the body axis. A representative healthy volunteer performed a normal gait (a), short-stepped gait (b), short-stepped shuffling wide-based gait (c), and short-stepped shuffling closed gait (d) while the TDPT-GT application was recording. The colored coordinates and lines of the head (brown), neck (sky blue), right hip joint (orange), left hip joint (gold), right knee (red), left knee (yellow), right heel (blue), left heel (dark green), right toe (light green), and left toe (pale green) were plotted on the sagittal and axial projection planes.
Figure 3. Two-dimensional relative coordinates projected onto the sagittal and axial planes of the body axis. A representative healthy volunteer performed a normal gait (a), short-stepped gait (b), short-stepped shuffling wide-based gait (c), and short-stepped shuffling closed gait (d) while the TDPT-GT application was recording. The colored coordinates and lines of the head (brown), neck (sky blue), right hip joint (orange), left hip joint (gold), right knee (red), left knee (yellow), right heel (blue), left heel (dark green), right toe (light green), and left toe (pale green) were plotted on the sagittal and axial projection planes.
Sensors 23 00617 g003
Figure 4. Distribution of angle ranges of the hip joint. The three violin plots in the upper row show the distribution of the angle range of the hip joint for the shuffling gait, short-stepped gait, and wide-based gait in the first 15 healthy volunteer group. The histogram in the lower row shows the distribution of the angle range of the hip joint for 47 patients with idiopathic normal pressure hydrocephalus (salmon pink) and 92 healthy elderly individuals aged 60 years or older in the Takahata cohort (blue emerald). The purple dotted line and number 30 indicate the proposed optimal cutoff value, and the green number in the violin plot indicates the median value.
Figure 4. Distribution of angle ranges of the hip joint. The three violin plots in the upper row show the distribution of the angle range of the hip joint for the shuffling gait, short-stepped gait, and wide-based gait in the first 15 healthy volunteer group. The histogram in the lower row shows the distribution of the angle range of the hip joint for 47 patients with idiopathic normal pressure hydrocephalus (salmon pink) and 92 healthy elderly individuals aged 60 years or older in the Takahata cohort (blue emerald). The purple dotted line and number 30 indicate the proposed optimal cutoff value, and the green number in the violin plot indicates the median value.
Sensors 23 00617 g004
Figure 5. Distribution of angle ranges of the knee joint.
Figure 5. Distribution of angle ranges of the knee joint.
Sensors 23 00617 g005
Figure 6. Distribution of relative vertical amplitudes of the heel.
Figure 6. Distribution of relative vertical amplitudes of the heel.
Sensors 23 00617 g006
Figure 7. Distribution of foot angles. The three violin plots in the upper row show the distribution of the foot angle for the shuffling gait, short-stepped gait, and wide-based gait in the first 15 healthy volunteer group. The histogram in the lower row shows the distribution of the angle range of the hip joint for 47 patients with idiopathic normal pressure hydrocephalus (salmon pink) and 92 healthy elderly individuals aged 60 years or older in the Takahata cohort (blue emerald). The purple dotted line and number 30 indicate the proposed optimal cutoff value, and the green number in the violin plot indicates the median value.
Figure 7. Distribution of foot angles. The three violin plots in the upper row show the distribution of the foot angle for the shuffling gait, short-stepped gait, and wide-based gait in the first 15 healthy volunteer group. The histogram in the lower row shows the distribution of the angle range of the hip joint for 47 patients with idiopathic normal pressure hydrocephalus (salmon pink) and 92 healthy elderly individuals aged 60 years or older in the Takahata cohort (blue emerald). The purple dotted line and number 30 indicate the proposed optimal cutoff value, and the green number in the violin plot indicates the median value.
Sensors 23 00617 g007
Figure 8. Distribution of heel outward shifts.
Figure 8. Distribution of heel outward shifts.
Sensors 23 00617 g008
Figure 9. Distribution of toe outward shifts.
Figure 9. Distribution of toe outward shifts.
Sensors 23 00617 g009
Figure 10. Distribution of leg outward shifts.
Figure 10. Distribution of leg outward shifts.
Sensors 23 00617 g010
Table 1. Clinical characteristics of three datasets.
Table 1. Clinical characteristics of three datasets.
Primary Healthy VolunteerTakahata CohortiNPH Patient
Total Number159247
Sex (male:female)9:636:5632:15
Mean ± SD 1 of age (years)39.1 ± 20.173.0 ± 6.377.3 ± 6.3
Range of age (years)22–7860–9161–87
Shuffling gait0026
Short-stepped gait0038
Wide-based gait0035
Freezing gait008
Spastic gait002
Instability0041
Alzheimer’s disease004
Fall history 01930
1 SD: standard deviation.
Table 2. Detectability of a shuffling gait, short-stepped gait, and wide-based gait.
Table 2. Detectability of a shuffling gait, short-stepped gait, and wide-based gait.
Cutoff ValueAUC * (95% CI) SensitivitySpecificityOD # (95% CI)
Shuffling gait
Angle range of the hip joint (°)3077.1 (72.0–82.2)69.0 74.3 6.39 (3.90–10.64)
Angle range of the knee joint (°)4578.6 (73.7–83.4)52.6 87.5 7.71 (4.54–13.28)
Relative vertical amplitude of the heel0.171.7 (66.3–77.2)59.5 76.0 4.64 (2.87–7.57)
Short-stepped gait
Angle range of the knee joint (°)4581.7 (77.4–86.0)43.0 93.4 10.55 (5.63–21.12)
Short-stepped gait
Heel outward shift−0.0878.1 (73.6–82.5)81.4 60.9 3.36 (2.16–5.27)
Toe outward shift0.1869.8 (64.6–75.0)70.2 58.8 6.77 (4.15–11.30)
Leg outward shift0.180.6 (76.4–84.8)81.4 63.8 7.65 (4.68–12.80)
* AUC: area under the receiver-operating characteristic curve. 95% CI: 95% confidence interval. # OD: odds ratio.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yamada, S.; Aoyagi, Y.; Iseki, C.; Kondo, T.; Kobayashi, Y.; Ueda, S.; Mori, K.; Fukami, T.; Tanikawa, M.; Mase, M.; et al. Quantitative Gait Feature Assessment on Two-Dimensional Body Axis Projection Planes Converted from Three-Dimensional Coordinates Estimated with a Deep Learning Smartphone App. Sensors 2023, 23, 617. https://doi.org/10.3390/s23020617

AMA Style

Yamada S, Aoyagi Y, Iseki C, Kondo T, Kobayashi Y, Ueda S, Mori K, Fukami T, Tanikawa M, Mase M, et al. Quantitative Gait Feature Assessment on Two-Dimensional Body Axis Projection Planes Converted from Three-Dimensional Coordinates Estimated with a Deep Learning Smartphone App. Sensors. 2023; 23(2):617. https://doi.org/10.3390/s23020617

Chicago/Turabian Style

Yamada, Shigeki, Yukihiko Aoyagi, Chifumi Iseki, Toshiyuki Kondo, Yoshiyuki Kobayashi, Shigeo Ueda, Keisuke Mori, Tadanori Fukami, Motoki Tanikawa, Mitsuhito Mase, and et al. 2023. "Quantitative Gait Feature Assessment on Two-Dimensional Body Axis Projection Planes Converted from Three-Dimensional Coordinates Estimated with a Deep Learning Smartphone App" Sensors 23, no. 2: 617. https://doi.org/10.3390/s23020617

APA Style

Yamada, S., Aoyagi, Y., Iseki, C., Kondo, T., Kobayashi, Y., Ueda, S., Mori, K., Fukami, T., Tanikawa, M., Mase, M., Hoshimaru, M., Ishikawa, M., & Ohta, Y. (2023). Quantitative Gait Feature Assessment on Two-Dimensional Body Axis Projection Planes Converted from Three-Dimensional Coordinates Estimated with a Deep Learning Smartphone App. Sensors, 23(2), 617. https://doi.org/10.3390/s23020617

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop