Next Article in Journal
Microsystems for Cell Cultures
Next Article in Special Issue
The Relevance of Calibration in Machine Learning-Based Hypertension Risk Assessment Combining Photoplethysmography and Electrocardiography
Previous Article in Journal
Review of Label-Free Monitoring of Bacteria: From Challenging Practical Applications to Basic Research Perspectives
Previous Article in Special Issue
Explainable Artificial Intelligence and Wearable Sensor-Based Gait Analysis to Identify Patients with Osteopenia and Sarcopenia in Daily Life
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proof of Concept in Artificial-Intelligence-Based Wearable Gait Monitoring for Parkinson’s Disease Management Optimization

by
Robert Radu Ileșan
1,2,
Claudia-Georgiana Cordoș
3,
Laura-Ioana Mihăilă
3,
Radu Fleșar
4,
Ana-Sorina Popescu
1,
Lăcrămioara Perju-Dumbravă
1 and
Paul Faragó
3,*
1
Department of Neurology and Pediatric Neurology, Faculty of Medicine, University of Medicine and Pharmacy “Iuliu Hatieganu” Cluj-Napoca, 400012 Cluj-Napoca, Romania
2
Clinic of Oral and Cranio-Maxillofacial Surgery, University Hospital Basel, CH-4031 Basel, Switzerland
3
Bases of Electronics Department, Faculty of Electronics, Telecommunications and Information Technology, Technical University of Cluj-Napoca, 400114 Cluj-Napoca, Romania
4
Computer Science, Faculty of Mathematics and Computer Science, West University of Timișoara, 300223 Timișoara, Romania
*
Author to whom correspondence should be addressed.
Biosensors 2022, 12(4), 189; https://doi.org/10.3390/bios12040189
Submission received: 31 January 2022 / Revised: 10 March 2022 / Accepted: 17 March 2022 / Published: 23 March 2022
(This article belongs to the Special Issue Intelligent Biosignal Processing in Wearable and Implantable Sensors)

Abstract

:
Parkinson’s disease (PD) is the second most common progressive neurodegenerative disorder, affecting 6.2 million patients and causing disability and decreased quality of life. The research is oriented nowadays toward artificial intelligence (AI)-based wearables for early diagnosis and long-term PD monitoring. Our primary objective is the monitoring and assessment of gait in PD patients. We propose a wearable physiograph for qualitative and quantitative gait assessment, which performs bilateral tracking of the foot biomechanics and unilateral tracking of arm balance. Gait patterns are assessed by means of correlation. The surface plot of a correlation coefficient matrix, generated from the recorded signals, is classified using convolutional neural networks into physiological or PD-specific gait. The novelty is given by the proposed AI-based decisional support procedure for gait assessment. A proof of concept of the proposed physiograph is validated in a clinical environment on five patients and five healthy controls, proving to be a feasible solution for ubiquitous gait monitoring and assessment in PD. PD management demonstrates the complexity of the human body. A platform empowering multidisciplinary, AI-evidence-based decision support assessments for optimal dosing between drug and non-drug therapy could lay the foundation for affordable precision medicine.

1. Introduction

More than 200 years ago, in 1817, Dr. James Parkinson published a scientific work entitled “The Essay on Shaking Palsy” [1] and with it the foundation of the disease that bears his name. After two centuries, we still struggle to understand and treat neurodegenerative diseases, such as Parkinson’s disease (PD), which are growing exponentially, especially in industrialized regions, and given that no one is immune to them, specialists are concerned about “the Parkinson pandemic” [2].
In 2017, we started the PDxOne project, intending to optimize the management of PD by offering sustainable automated or semiautomated solutions.

1.1. Considerations on Parkinson’s Disease

Parkinson’s disease (PD) is a common progressive neurodegenerative disorder that can cause significant disability and decreased quality of life [3]. PD is today more present than ever, being the second most common neurodegenerative disease after Alzheimer’s disease and affecting 0.3% of the population [4]. It is estimated that there are 6.2 million people diagnosed with PD disease worldwide, and the disease caused the death of 117,000 people in 2015 alone [5]. The American Parkinson Disease Association estimates that there are already 1 million people with PD living in the U.S. alone and over ten million worldwide, giving some researchers the evidence of a Parkinson pandemic [2]. While studies are divided on the prevalence by gender, the affinity of the pathology for the aging tissue, particularly neural tissue, is much clearer. Statistically, PD occurs in people over 60 years old, affecting 1% of this age group and increasing to 4% for people over 80 [4]. Unfortunately, PD can also occur in younger people under 50, known as Young Onset Parkinson’s disease (YOPD). Studies show that 5–10% of patients diagnosed with PD are between 20 and 50 years old [6]. It is becoming even more complex with the literature showing cases of patients diagnosed with PD who were younger than 20—even some rare cases with patients under ten years of age, and the first symptoms appeared as early as two years of age [7]. This form is known as the juvenile form of PD and was first described in 1875. It has a substantial genetic component that can be diagnosed today with genetic testing [7].
In 2014, Mary Ann Thenganalt et al. performed a systematic review of articles cited in PubMed between 1980 and 2013. They concluded that PD could be divided into several subtypes, the most representative being tremor-dominant and postural instability gait difficulty form (PIGD) [8]. The clinical presentation of PD can take the form of many symptoms, from which the easiest to notice are the motor impairments [9]. These symptoms are caused by the degeneration of the dopaminergic neurons located in the substantia nigra from the ventral midbrain [10]. According to the Movement Disorder Society (MDS), the clinical diagnosis of PD is based on the presence of bradykinesia, along with either rest tremor or rigidity [11].
Motor impairment in PD has long been the focus of researchers, with significant advances being made in diagnostic accuracy, implementation, improvement of more accurate assessment scales, and better management of therapeutic strategies [12]. Unfortunately, at the moment of the diagnosis, a significant number of neurons that produce dopamine are already dysfunctional [6]. During the 15–20 years before the onset of motor symptoms, the patients experience a phase called “prodromal PD”, during which the neurodegeneration starts and progresses [13]. It has been proved that olfaction impairment, constipation, depression, rapid eye movement (REM) and sleep behavior disorder (RBD) can be present in the prodromal period of PD [14].
Although the focus of research has been on motor symptoms, clinical studies have shown that non-motor symptoms in PD, such as depression, pain, psychosis and sleep disturbances, should be regarded as equally important when analyzed using quality-of-life questionnaires as well as economic and health indicators [12]. Therefore, we are dealing with a pathology that is increasingly present in our lives, that will double its number globally by 2030 [15], and that will have a significant socio-economic impact. Socially, PD patients become isolated, stigmatized [16], and even discriminated, with severe implications for the course of the disease and the clinical picture. To diagnose PD in a patient means to put a verdict that irreversibly alters their lifestyle and that of their families. Once diagnosed, the patient will also undergo a lifelong treatment, which aims not to cure but to improve or stagnate the symptoms’ evolution. Studies, therefore, show a mortality rate double that of the healthy population, which increases and presents a more aggressive clinical picture in patients with YOPD or juvenile PD [17]. Therefore, the aspect of early onset raises many questions about our current ability to understand the pathophysiology of PD and treatment errors that may result in individual and global socio-economic consequences.
Indeed, medical assistance for patients with PD is a major drain on the healthcare budget. This financial amount is given by the complexity of PD, which affects the patient on several levels at once. To understand the real-life implications that affect every one of us and why we should spend the necessary resources, a quantification of socio-economic indicators is necessary. One figure we could start with is the financial effort that countries in Europe make to treat patients with PD, i.e., almost EUR 14 billion per year [18]. Interestingly, this figure also represents the amount the USA is spending per year treating PD patients, namely USD 14.4 billion [19]. The average annual cost per patient for PD in Germany is EUR 20,095 [20]. On average, direct costs represent 65.5% and indirect costs 34.5% [21]. Of the total direct costs (EUR 13,158), EUR 3526 is spent on medication, which is also the largest expenditure, and EUR 3789 is spent annually on hospitalization and home care costs [21]. Costs for home care by the family amount to 20% of the direct costs [21]. The same study shows a directly proportional relation between costs and disease progression, i.e., EUR 18,660 annual costs for stage 1–2, increase to EUR 31,660 annual costs for stage 2–5 (according to Hoehn and Yahr) [21]. Annual costs for PD differ quite a lot between European countries. While in Russia, EUR 5240 is spent per patient per year [21], in England, the annual costs for a patient with advanced stage PD (3–5 according to Hoehn and Yahr) can reach EUR 72,277 [22]. Attention is required, especially regarding the increase in treatment costs if patients are taken out of their environment, out of their home, and moved to a nursing home. In this case, studies show a 500% increase in costs for PD treatment [23]. All these costs are strictly related to PD, but a patient may also have other associated pathologies aggravated by PD and vice versa. Moreover, the financial impact on society is difficult to quantify because it should be taken on a patient-by-patient basis.
A staggering EUR 798 billion is spent annually at the European level to treat brain diseases, according to a study carried out in 2010 [19]. This amount is almost four times higher than Romania’s GDP. It shows the importance of continuous research to provide a sustainable medical system for Europe and beyond in a growing and aging population that wants to maintain its standard of living in old age.
These figures motivate the ongoing research on wearables for the early detection of PD-specific symptomatology and prediagnosis of PD in incipient stages, as well as long-term monitoring of the disease in a ubiquitous healthcare environment which provides intelligent decision support algorithms for assessment and patient-specific treatment plans in PD.

1.2. Related Work—Wearables in PD Monitoring

For exemplification, Boroojerdi et al. report on the employment of the NIMBLE wearable biosensor patches, composed of an accelerometer and an electromyography (EMG) sensor, for motor evaluation in PD [24]. As for another example, Jauhiainen et al. report on the employment of a Movesense sensor and a Forciot insole to observe walking patterns in PD [25]. Phan et al. report on the use of BioKin devices to assess daily tasks: pointing, pouring, walking, and walking around a chair [26]. Lonini et al. report on employing BioStamplRC flexible wearable sensors, consisting of a tri-axial accelerometer and gyroscope, to record motion data [27].
Continuous long-term monitoring of motor symptoms in PD using inertial sensors is described by Borzì et al., aiming for the identification of bradykinesia and FOG [28], or by Powers et al., aiming for the identification of tremors and dyskinesia [29]. The employment of built-in smartphone sensors, with dedicated smartphone applications, is described by Heijmans et al. in [30] or by Motolese et al. in [31], for the remote monitoring of the PD patients during daily activities.

1.3. Related Work—AI-Based Decisional Support in PD Assessment

Some noteworthy examples regarding AI-based decisional support are presented as follows: Lonini et al. report on the employment of Random Forest classifiers to identify bradykinesia and tremor [27]. Random Forest for classification in PD was also studied by Aich et al. in [32], along with support vector machine, K-Nearest Neighbor, and Naïve Bayes. PD-specific symptomatology detection and classification using convolutional neural networks (CNN) was reported by Taewoong et al., who assessed daily activities based on 3D acceleration and angular velocity data measured with a Microsoft Band 2 [33]. Further employment of CNN was reported by Lonini et al. for the identification of bradykinesia and tremor [27], and by Steinmetzer et al. for the arm oscillation monitored with Mbientlab portable Motion Rectangle sensor bracelets under a Timed Up and Go (TUG) test scenario [34].

1.4. This Work

The PDxOne research project desires to develop and use the latest technologies to collect and, above all, interpret medical data. With a world population reaching 8 billion people and with today’s medical requirements and demands, these tasks simply cannot be hand-operated anymore. Given the sheer medical data volume that must be collected every day, the request for economical and sustainable solutions forces healthcare systems to embrace a worldwide digitalized implementation.
Thus, our work is placed in the context of today’s demand for ubiquitous monitoring and intelligent decisional support in healthcare. We target to develop, at the end of our project, a small-size wearable and portable monitoring system for patients diagnosed with PD, aiming for long-term quantitative and qualitative assessment of the PD symptomatology in a continuous fashion, and intended for AI-based decisional support in the discrimination of the pathology and formulation of dedicated treatment plans, which constitutes a novelty in the field.
This article targets the monitoring and assessment of gait in patients diagnosed with PD. Gait disorders are a hallmark of the condition and are associated with a loss of independence and an increased risk of falls. Disturbances of the gait, even if hardly noticeable, are described from the earliest stages of the disease [35] and include shuffling gate, shortened stride length, reduced overall velocity, and increased stance phase (up to doubling), along with reduced or absent arm swing, reduced trunk rotation, and decreased amplitude of motion in the hips, knees, and ankles [3]. In advanced stages, gait disorders often become increasingly complex, including motor blocks, festination, and imbalance [36].
Multiple studies have also been conducted for the early detection of motor deficient behaviors to apply proper therapeutic interventions, which are proved to slow down the motor dysfunction and maintain functional independency (in patients with preserved cognitive function) [3,37,38,39]. Symptoms such as dyskinesia, which is induced by therapy and manifests as involuntary movement of any body parts, appear in advanced stages of PD [9]. The symptomatic therapy for the classic motor features is usually satisfactory, but antiparkinsonian therapy that does not induce motor complications is still needed [40].
The solution proposed for ubiquitous gait monitoring and AI-based decisional support in gait assessment is envisioned in the shape of a wearable physiograph. The proposed physiograph performs bilateral tracking of the foot biomechanics assessed by means of plantar pressure distribution and lower-limb EMG, in correlation to upper limb balance, which is evaluated by means of arm balance magnitude of acceleration (MA) and variation of acceleration (VA). The recorded signals are transmitted over a Bluetooth radio link to a mobile device, e.g., smartphone or tablet. They are uploaded and stored into an online database and made available for future access, either in real-time or offline, for processing and interpretation.
The proposed physiograph enables both qualitative and quantitative assessment of gait. As such, we perform gait evaluation based on biomechanical parameters, expressed in terms of arm balance, heel strike, and lift-off, and temporal parameters, expressed in terms of cadence, single support, double support, single support to double support ratio, and stride time variability. Next, we evaluate the physiological interdependencies involved during the gait cycle by applying the cross-correlation function to each recorded signal pair. We illustrate that PD-specific gait is identifiable based on the evaluated gait assessment parameters following the evaluation results. Consequently, the biomechanical and temporal parameters and the cross-correlation results are applicable as inputs to an expert system for identifying and discriminating PD-specific gait pathology.
The novelty of the proposed physiograph consists of the underlying AI-based decisional support procedure for gait assessment. We generate a correlation coefficient matrix from the gait monitoring signals to visually represent the gait pattern. Gait assessment using the biomechanical and temporal parameters and the cross-correlation function is contained in the correlation coefficient matrix. Then, we apply the surface plot of the correlation coefficient matrix to a convolutional neural network (CNN) for gait classification.
A proof of concept of the proposed physiograph with AI-based decisional support is validated in the clinical environment on a group of ten subjects consisting of five PD patients and five healthy controls. As such, the proposed solution provides a feasible method for AI-based support for gait monitoring and assessment in a ubiquitous healthcare environment.

2. Materials and Methods

This paper proposes a wearable miniature physiograph with AI-based decisional support for gait monitoring and assessment in PD. Gait evaluation is performed in accordance with the Unified Parkinson’s Disease Rating Scale (UPDRS)—motor subscale, and the Movement Disorder Society UPDRS (MDS-UPDRS) [41,42].
The proof of concept of the proposed wearable gait monitoring physiograph was tested extensively in the laboratory and validated indoors in the clinical environment, with a study group consisting of five patients diagnosed with PD and five healthy controls. The PD group includes three males and two females. The healthy control group includes four males and one female. The healthy controls do not have any previously diagnosed neurodegenerative disorder or podiatric condition.
All procedures performed in this study involving human participants were following the ethical standards of the institutional and/or national research committee. Informed consent was obtained from all individual participants involved in the study.

2.1. The Proposed Physiograph for Gait Monitroing in PD

The proposed gait monitoring physiograph is presented in the block diagram from Figure 1a and the practical realization from Figure 1b.
The proposed wearable physiograph is developed around an ATmega2560 microcontroller (µC), which reads six Aidong IMS C20B thin-film resistive pressure sensors and four EMG channels over the analog ports and a LSM9DS0 module over the I2C interface. Signal acquisition is performed with synchronized sampling, with an fs = 100 Hz sampling frequency and an on-chip 10-bit analog-to-digital converter (ADC). Under this setup, the proposed physiograph performs bilateral tracking of the foot biomechanics through the plantar pressure progression pattern, lower-limb muscular activation, and unilateral monitoring of the arm balance.
The µC development board is attached to a Velcro strip and is worn around the user’s waist, as illustrated in Figure 2.
Foot biomechanics is assessed using three pressure sensors and two EMG channels, clustered into a foot biomechanics assessment module [43]. Two such modules are considered for bilateral monitoring.
Operation of the foot biomechanics assessment module is described as follows: bilateral tracking of the plantar pressure progression pattern during the gait cycle is performed using three pressure sensors, attached onto an insole below the toe (FSR0), metatarsal arch (FSR1), and heel area (FSR2), respectively, following the center of pressure (COP) progression line, as illustrated in Figure 3.
The sensors are deployed into a resistive divider topology with a R = 1 MΩ resistance, as illustrated in Figure 4a, and operate as force sense resistors (FSR) with the sensor resistance value derived as
F S R = V F S R V D D V F S R · R ,
where VDD = 5 V is the supply voltage and VFSR is the FSR voltage drop. The sensor resistance (kΩ) can be converted to mass (kg) according to the mass vs. resistance characteristics provided in the sensor datasheet and plotted in blue in Figure 4b. Mass can further be converted to pressure (kg/cm2) by dividing the mass to the sensor area.
In this work, we target gait pattern assessment rather than podiatric assessment. As such, the FSR resistance derived with (1) is sufficient to indicate the application of plantar pressure. In addition to (1), we have changed the polarity of the FSR signal,
F S R = max ( F S R ) F S R ,
to have the “HIGH” signal level indicating pressure, and the “LOW” signal level indicating absence of pressure. This also changes the mass vs. resistance characteristics as plotted in red in Figure 4b.
Tracking of the lower-limb muscular activation pattern during the gait cycle is performed with two EMG channels which acquire the EMG of the Tibialis anterior (TA) and Gastrocnemius medialis (GM) muscles. Off-the-shelf MikroElektronika EMG Click boards were used for each EMG channel analog front end (AFE), respectively. Wet Ag/AgCl electrodes were employed for EMG acquisition. Electrode placement is illustrated in Figure 5, with the active electrodes (white and red) placed onto the TA and GM muscles and the reference electrodes (black) placed onto the lateral and medial malleolus, respectively.
A 10× AFE gain is set from the on-board potentiometer to accommodate the prescribed 1 µV–10 mV EMG amplitude range, accounting for 2 mV motor unit action potential (MUAP) amplitude of the healthy muscle, 0.5 mV MUAP amplitude for primary muscular disease, as well as 10 mV MUAP amplitude of intramuscular sprouting and chronic partial denervation [43,44,45,46]. On-board filtering is performed with three analog filter stages: two high-pass filters with the cutoff frequencies set to 1.6 Hz and 0.16 Hz, respectively, and a low-pass filter with the cutoff frequency set to 60 Hz.
After acquisition, EMG signal processing accounts for averaging with an 8-sample rectangular window with 50% overlap, a 4th order Butterworth approximation high-pass filter with fcL = 0.5 Hz to suppress the DC component, and then a 4th order Butterworth approximation low-pass filter with fcH = 10 Hz. To be noted is that the low-pass frequency of 10 Hz was considered as we were interested in the identification rather than evaluation of muscular activity [47].
The accelerometer from a LSM9DS0 module is employed to perform arm balance monitoring. The sensor is attached to the patient’s right-hand wrist using a Velcro strip and, if necessary, tightened with an adhesive band, as illustrated in Figure 6.
The accelerometer was configured for a 2G acceleration range, and the sensor data were read using the Adafruit LSM9DS0 library. Accelerometer signal processing assumes averaging with an 8-sample rectangular window with 50% overlap, a 4th order Butterworth approximation high-pass filter with fcL = 0.5 Hz to suppress the DC component standing for the accelerometer initial position [48], and then a 4th order Butterworth approximation low-pass filter with fcH = 30 Hz. These filter specifications cover the targeted 1 Hz–10 Hz frequency range of gait-related informational content (most relevant information is available up to 4 Hz) [49], as well as the 4 Hz–6Hz frequency range of tremor [50]. Additionally, the low-pass filter suppresses higher frequency components due to, for example, vibrations as well as noise. The raw signals on the three axes are then converted to acceleration and expressed in m/s2.
Two metrics are employed for arm motion tracking based on the dynamic acceleration, defined as follows: the magnitude of acceleration (MA) is determined by applying the Pythagorean theorem to the readings on the three axes, respectively [51], according to equation:
M A = x 2 + y 2 + z 2 ,
and defines the absolute acceleration value. The variation of acceleration (VA) has the average of the past readings subtracted from each axis, respectively [47], as defined by equation:
V A ( k ) = ( x ( k ) x a v g ( k ) ) 2 + ( y ( k ) y a v g ( k ) ) 2 + ( z ( k ) z a v g ( k ) ) 2 ,
where k is the current index and
x a v g ( k ) = 1 k 1 · i = 1 k 1 x i ,
y a v g ( k ) = 1 k 1 · i = 1 k 1 y i ,
z a v g ( k ) = 1 k 1 · i = 1 k 1 z i ,
are the average of the past readings, i.e., up to index k − 1 on each axis, respectively.
The proposed physiograph is aimed at long-term monitoring in a ubiquitous healthcare environment. Wireless connectivity is achieved by deploying the proposed physiograph with a HC-05 Bluetooth module. After acquisition, the raw data are sent over UART to the HC-05 module and transferred to an Android mobile terminal, e.g., a smartphone or tablet. The mobile terminal collects the user data from the physiograph over Bluetooth and assembles it into JavaScript Object Notation (JSON) files. The JSON files are sent one by one over the Internet to the server via REST API and stored in an online database. From this database, the signals are available for later retrieval to a desktop computer. A diagram of the application is presented in Figure 7.
An example of a message sent to the API via the mobile phone is provided in Figure 8. The data will be stored in an online database. An example of the database content in Figure 9 illustrates that the data dictionary holds all the data acquired with the gait monitoring physiograph and is made available for later retrieval onto a laptop or personal computer for processing and interpretation.

2.2. Gait Assessment—Correlation

The study group was instructed to undertake steady-state walking, at a pace of their own choice, and walk around the room. One or two walking trials were performed to make sure that the users were comfortable with the wearable devices and that they understood the requirements of the exercise. The subsequent steady-state walking activity was then recorded for the proposed gait monitoring and assessment procedure. Photographs taken during the trials of gait assessment are illustrated in Figure 10.
A physiological gait cycle is considered from one heel strike to the next heel strike of the same foot and consists of a stance and a swing phase, respectively [52,53]. The stance phase is further split into:
  • Heel strike—with pressure applied onto the heel area;
  • Support (or foot flat);
  • Midstance;
  • Heel off—with pressure applied onto the metatarsal arch and hallux;
  • Toe off—with pressure applied only onto the hallux before lift-off [43].
The physiological gait cycle is associated with arm balance, which describes a forward sway during the stance phase (from heel strike to lift-off) and a backwards sway during the swing phase (from lift-off to the next heel strike) [43]. This is visible on the arm balance MA waveform which describes a U-shaped pattern. The MA maxima account for the arm sway direction changes, constituting a good indicator for the stance phase initiation and ending.
The gait pattern in PD differs from the physiological gait. The literature describes a flat-foot strike for the PD gait pattern, or toe-to-heel plantar pressure progression in more advanced stages [54], associated with reduced lifting of the foot after lift-off [55] and limited or no arm balance along the gait cycle [56]. As such, the MA waveform exhibits a larger number of peaks, corresponding to the oscillations of the body’s center of mass and tremor. In this case, the VA waveform exhibits a larger variability corresponding to the increased number of MA peaks.
Due to the large number of peaks, the MA waveform cannot be employed to provide indication regarding stance initiation and ending in PD. In this work, we have rather employed the plantar pressures for stance identification. Plantar pressure detection was performed by comparing the FSR value to an empirical threshold level FSRth computed as a fraction of the FSR signals. Accordingly, one stance phase ranges from the first to the last occurrence of plantar pressure, regardless of which pressure point, as illustrated in Figure 11. The stance phases identified in this manner account for the signal frames applied for cross-correlation in the gait assessment procedure described further on.
Physiological gait assumes a precisely defined interdependency between arm balance, plantar pressure, and lower-limb muscular activation. We evaluate signal interdependency in the time domain using the cross-correlation function given by
R s i g 1 , s i g 2 ( m ) = { n = 0 N m + 1 s i g 1 ( n + m ) · s i g 2 ( n ) ,   m 0 R s i g 1 , s i g 2 ( m ) ,   m < 0 ,
where sig1 and sig2 are the signal frames being correlated, m is the cross-correlation index, and N is the frame length [57]. The cross-correlation function defined in (8) provides a measure of the similarity between the two signals sig1 and sig2 as a function of m. As such, cross-correlation maxima in the origin account for the identification of signal interdependencies. Provided the cross-correlation maxima are situated outside the origin, the index of cross-correlation peaks accounts for the displacement between the signals.
First, we assess whether the arm balance MA peaks, corresponding to arm balance initiation and ending, are synchronous with the stance initiation and ending determined from the plantar pressure progression pattern. This should be the case for physiological gait. Next, we evaluate the cross-correlation functions for:
  • arm balance MA vs. lower-limb muscular activation signals, i.e., TA and GM respectively,
  • lower limb muscular activation signals vs. FSR signals respectively.
As we move forward, we employ the correlation coefficient matrix to quantify the interdependency between either signal pair. A generic correlation coefficient matrix for M signals is expressed as:
R = ( 1 ρ ( s i g 1 , s i g 2 ) ρ ( s i g 1 , s i g M ) ρ ( s i g 2 , s i g 1 ) 1 ρ ( s i g 2 , s i g M ) ρ ( s i g M , s i g 1 ) ρ ( s i g M , s i g 2 ) 1 ) ,
where
ρ ( s i g i , s i g j   ) = 1 N 1 n = 1 N ( s i g i ( n ) μ s i g i σ s i g i ) ( s i g j ( n ) μ s i g j σ s i g j ) ,   i , j = 1 , M ¯ ,
is the Pearson correlation coefficient of two signals sigi and sigj; µ and σ are the mean and standard deviation of the signals indicated in the signal subscripts, respectively; and N is the signal length [58,59]. We aim to generate the correlation coefficient matrix for the signal space
S I G = [ M A ,   V A , T A l e f t , G M l e f t , F S R 0 , l e f t , F S R 1 , l e f t F S R 2 , l e f t , T A r i g h t , G M r i g h t , F S R 0 , r i g h t , F S R 1 , r i g h t F S R 2 , r i g h t ] ,
consisting of the arm balance MA and VA (determined from the accelerometer signals), lower-limb EMG and plantar pressures. To be noted is that the main diagonal of the cross-correlation matrix consists of unity elements and accounts for the fact that the signals are correlated to themselves.
In contrast to (8), the definition of the correlation coefficients given in (10) does not account for the displacement between the signals, but only provides a quantification for signal similarity. Physiological delays originating from the biomechanical processes involved during gait, e.g., arm balance initiated before or after heel strike, arm balance terminated before or after heel strike, etc., which are determined using (8) as shifting of the cross-correlation peak form the origin, are missed using the correlation coefficients in (10).
To address the displacement of the signals in between one another and visualize them on the correlation coefficient plot, we have generated 10 shifted versions of the signal frame. Consequently, we extended the signal space to 120 signals which are to be correlated, resulting in a 120 × 120 correlation coefficient matrix. A 10 × 10 section from this matrix illustrates the interdependency between either shifted versions of the signals, rather than the signals themselves. Then, the largest coefficient value, accounting for the best similarity, determines the lag between the signals.

2.3. Gait Assessment—AI-Based Decisional Support

AI-based decisional support for the identification of PD gait pattern is implemented in this work using convolutional neural networks (CNN). The CNN is a deep learning algorithm that takes an image as input, assigns importance to the features in the image, and can differentiate them from each other. Thus, this type of network has the ability to extract local features based on the convolution operation between the original bidimensional data and certain series of the convolution kernels. The preprocessing required in a CNN is much lower compared to other classification algorithms, and such networks are used in applications for image recognition.
One of the benefits of deep learning is the ability to generalize and to learn massive amounts of data. Good network generalization capacities are obtained by accounting for the relationship between the size of the learning database and the complexity of the network architecture. The higher this ratio, the better the network performance on the test dataset. Furthermore, a big advantage of CNN networks is the weight sharing feature, which reduces the number of trainable network parameters and in turn helps the network to enhance generalization and avoid overfitting.
In this study, we used several architectures such as MobileNet, EfficientNetB0, and Xception. MobileNet is a CNN architecture model used for image classification and mobile vision. The advantage of this network is the very low computing power to apply transfer learning, because the model is based on depthwise separable convolution that has the effect of reducing the calculations and the size of the model. MobileNet uses 3 × 3 depthwise separable convolutions, using 8 to 9 times fewer calculations than standard convolutions with only a small reduction in accuracy. Counting the deep convolutions as separate layers, MobileNet has 28 layers [60]. The EfficientNet model is based on the uniform scaling of the network width, depth, and resolution to improve performance. This network has been extended to a family of deep learning architectures with very good accuracy and efficiency [61]. Xception is built on two main points: depthwise separable convolution, i.e., a depthwise convolution followed by a pointwise convolution, and shortcuts between convolution blocks.
The surface plot of the correlation coefficient matrix is saved as a jpeg image and is applied to the CNN for classification into physiological and pathological gait. The flowchart of the proposed solution is shown in Figure 12.
The parameters used to train the models are listed in Table 1. A very important parameter is the learning rate, which was chosen to be 0.05 for all models and has the role of controlling the model in response to the estimated error each time the model weights are updated. To reduce the nonlinearity of the output, the Softmax output layer activation function is used for all models. This function determines the type of predictions that the model can make. At the same time, the loss is the prediction error of the network, and the loss function has the role of determining the error. In the proposed binary classification system, the binary cross-entropy compares the predicted probability of the model with the actual result, which can be 0 or 1.
Models were trained on the graphic processing unit (GPU) in Google Colab using Keras. The motivation for Keras is ease of use and extension as neural layers, cost functions, optimizers, initialization schemes, and activation functions. As such, the activation functions are standalone modules that can be combined to create new models defined in Python. Keras offers scalability because it can run on tensor processing units (TPU) or large groups of GPUs, and the model can be exported to run in the browser or on a mobile device.
We have generated a total number of 236 images of correlation coefficient matrix plots, corresponding to the 10-subject database. This constitutes the data set for the CNN which aims to discriminate the walking pattern between physiological gait and PD. Since the data set is small, an augmentation was performed. The Adam optimizer was used to optimize the neural network. The RMSProp optimizer was used in the optimization of the EfficientNetB0 convolutional neural network. Model training was performed with 150 epochs for the MobileNet, EfficientNetB0, and Xception models. The batch size for MobileNet is 32, for EfficientNet it is 64, and for Xception it is 128. To evaluate the performance of each model, the data set was divided as follows: 70% for the training set, 10% for the validation set, and 20% for the test set.

3. Results

3.1. Gait Assessment—Time Domain

A section consisting of three gait cycles acquired during the walking trials is illustrated in Figure 13 for one healthy control and three PD patients. The arm balance MA and bilateral EMG of the TA and GM as well as the FSR resistance values, respectively, are plotted. Red markers are placed to indicate the MA maxima.
Figure 13a exhibits the physiological gait pattern (Healthy control 1). Indeed, the peaks of the arm balance MA correspond to the initiation and ending of the stance phases, and the plantar pressures follow the physiological heel→metatarsal arch→hallux progression pattern. Moreover, TA activity can be observed simultaneously with pressure under the heel area and GM activity simultaneously with pressure under the hallux.
In contrast to Figure 13a, the three-cycle gait section plotted in Figure 13b illustrates the gait pattern of a patient with PD (i.e., Patient 5). The walking pattern exhibits bilateral flat-foot strike, which is typical in incipient and mid-stage PD. Muscular activation is present in a direct link to the plantar pressure points. The magnitude of acceleration also exhibits the U-shaped variation due to arm balance, although it is more pronounced during the stance than the swing phase of the gait cycle. However, a series of small-amplitude local peaks originated by tremor are also visible.
The three-cycle gait section recorded on two other patients with PD (i.e., Patient 1 and 3) is plotted in Figure 13c,d. Both patients exhibited flat-foot strike. What stands out, however, on the MA waveform is the absence of arm balance, in which case the peaks are originated by the tremor.
One thing that stands out in Figure 13d for Patient 3 is an asymmetry between the left and right foot during gait. The left foot exhibits a flat-foot strike, typical for PD, although with reduced plantar pressure on the metatarsal arch. The right foot on the other hand exhibits a plantar progression pattern which has toe pressure exerted until the next heel strike with a slight overlap. This accounts for the fact that the patient does not lift the right foot during the sway phase, but rather pulls the right foot with continuous contact between the toe and the ground. According to the FSR plots, the stance phase is terminated when pressure is reduced under the metatarsal arch and tow, i.e., the spike on FSR0, and the forefoot is only lifted after heel contact.
Tremor on the accelerometer signals can be also assessed on the VA waveform. Regarding arm balance, the gait cycle in PD exhibits limited to no arm balance whatsoever. This is clearly visible on the signals acquired from the accelerometers. The arm balance MA of a healthy control, plotted in Figure 14a, exhibits signal peaks (indicated with a red triangle) which correspond to the stance phase initiation and termination, respectively, and a U-shaped waveform corresponding to arm sway. Consequently, the arm balance VA, plotted in Figure 14a, exhibits peaks during the sway.
In contrast, the arm balance MA of a patient with PD exhibits a larger number of peaks, as illustrated in Figure 14b. Although some MA peaks with larger amplitudes, originated by the oscillation of the body center of mass during gait, can be identified as stance initiation and termination (see Figure 13 for a clear identification), most peaks are rather tremor-induced local peaks. As such, the VA waveform reveals a larger variation.
The time-domain plots further enable the assessment of the biomechanical and temporal parameters, used for motor evaluation in accordance with the UPDRS and MDS-UPDRS rating scales for Parkinson’s Disease.
Heel strike accounts for initial contact with the ground at the stance initiation phase under a physiological gait pattern. Heel strike identification is performed by following the plantar pressure progression pattern. The identification of a heel strike is provided by a heel→metatarsal arch→toe pressure progression detection. For illustration, the FSR signals during gait initiation, with heel strike and flat-foot strike, are plotted in Figure 15.
Lift-off accounts for the fact that the foot is lifted from the ground during the swing phase, and consequently stance phase with the opposite foot. Accordingly, no plantar pressure should be recorded during the swing phase, provided the foot is fully lifted from the ground, accounting for the identification of lift-off. For illustration, the FSR signals during sway, with lift-off and without lift-off, are plotted in Figure 16.
The biomechanical parameters of the test group are listed in Table 2.
The biomechanical parameters of gait, assessed based on the signals recorded with the proposed physiograph, point out the asymmetry between the stance phases with the left and right foot, respectively. The patients exhibit flat-foot strike, which is typical for PD in incipient and middle stages. Only one instance of absent lift-off was observed for Patient 3. Arm balance was identifiable for Patients 1, 4, and 5 who exhibit a small-magnitude tremor.
The healthy control group was not previously diagnosed with any podiatric condition. As illustrated in Table 2, the healthy control group exhibits a physiological gait pattern, except for Healthy Control 2 who exhibits flat-foot strike.
The temporal gait analysis parameters are determined from the time-domain plots of the gait monitoring signals as follows [38,62]:
  • Cadence, i.e., the number of steps per minute, is determined as the number of gait cycles in a 10 s segment’s gait signals and extrapolated to one minute.
  • Single support, i.e., pressure is exerted by only one foot and not the other, determined from a 10 s segment of the gait signals, and expressed in percentages.
  • Double support, i.e., pressure is exerted by both feet during double limb support, determined from a 10 s segment of the gait signals, and expressed in percentages.
  • Ratio of the single support to double support (s/d) is the ratio for single support and double support.
  • Stride time variability is expressed as a coefficient of variation (CoV) computed from the mean (µstride) and standard deviation (σstride) of the strides over a 10 s segment of the gait signals [62,63,64], as given in the following equation
    C o V = σ s t r i d e μ s t r i d e · 100 ,
    and is expressed in percentages.
To be noted is that a 10 s interval is sufficient for the assessment considering that the walking trials assume steady-state walking. The temporal parameters of the test group are listed in Table 3.
What stands out in both Healthy control group and PD group is that the cadence is considerably smaller than the nominal 60 steps/min [52,53]. Both healthy controls and PD patients expressed that their pace was reduced because of physiograph wiring.
The PD group exhibits reduced percentage of single support and S/D in comparison to the healthy control group, while the percentage of double support is not significantly different. On the other hand, stride time variability is increased for the PD group in contrast to the healthy control group.
What stands out for the PD group is that Patient 4 exhibits a very small percentage of single support with the left leg. This is because Patient 4 has a very small right leg lift-off during the swing phase of the gait cycle.

3.2. Gait Assessment—Cross-Correlation

As discussed in Section 2, a physiological gait pattern assumes some specifically defined interdependencies between the plantar pressure, lower-limb muscular activation, and arm balance, assessed using the cross-correlation function. Illustration of the cross-correlation between the signals acquired with the proposed physiograph is presented next.
The arm balance magnitude of acceleration, the TA and GM activation, and the FSR signals for a stance-phase, acquired on a healthy control for both the left foot and the right foot, are plotted in Figure 17. Red markers indicate stance initiation and termination, respectively.
The cross-correlation of the arm balance MA to the lower-limb muscular activation is plotted in Figure 18. As illustrated for the left foot, RMA,EMG-TA exhibits a peak to the left of the origin, indicating a delay of 70 ms of the TA activation vs. the MA stance initiation peak. Similarly, RMA,EMG-GM exhibits a peak to the left of the origin, indicating an advancement of 70 ms of the GM activation vs. the MA stance termination peak. Nevertheless, both delays are accounted for as physiological. A similar reasoning can be formulated for the right foot, indicating a delay of 20 ms of the TA vs. the MA stance initiation peak and an advancement of 110 ms of the GM activation vs. the MA stance termination peak.
Next, the cross-correlation of the TA activation signal to the FSR signals is plotted in Figure 19. As illustrated, REMG-TA,FSR2 exhibits a cross-correlation peak in the origin, indicating a direct interdependency between TA contraction and plantar pressure under the heel area. In contrast, the peaks of REMG-TA,FSR1 and REMG-TA,FSR0 are shifted from the origin, indicating that plantar pressure under the metatarsal arch and the hallux are not directly linked to the TA. The lag indicates the displacement, i.e., delay, between heel strike (TA contraction), heel off (FSR1 and FSR0), and toe off (only FSR0). This reasoning holds for both feet.
A similar reasoning is formulated regarding the cross-correlation of the GM activation signal to the FSR signals plotted in Figure 20. As illustrated, REMG-TA,FSR0 exhibits a cross-correlation peak in the origin, indicating a direct interdependency between GM contraction and plantar pressure under the hallux. In contrast, the peaks of REMG-TA,FSR2 and REMG-TA,FSR1 are shifted from the origin indicating that plantar pressure under the heel area and metatarsal arch are not directly linked to the TA. The lag indicates the displacement, i.e., delay, between heel strike (FSR2), heel off (FSR1 and FSR0), and toe off (GM activation). This reasoning holds for both feet.
In contrast, the arm balance MA, the TA and GM activation, and the plantar pressures for a stance-phase acquired under a bilateral monitoring scenario from a PD patient are plotted in Figure 21.
Two aspects stand out in Figure 21. On one hand, the MA peaks are not in a direct correspondence with the onset and termination of stance. With this consideration, the cross-correlation of the MA to either one of the lower-limb muscles becomes irrelevant, as the cross-correlation peaks indicate correspondence of the muscular activation to MA maxima originating from tremor rather than arm balance. On the other hand, the gait pattern exhibits flat foot strike for the left leg, identified by having pressure exerted simultaneously by the heel, metatarsal arch, and hallux, respectively, rather than the physiological heel→metatarsal arch→toe progression. This latter remark results in cross-correlation maxima in the origin, or close to the origin, for all signal pairs, as plotted in Figure 22 and Figure 23.
The surface plot of the correlation coefficient matrix defined in (9) was employed in this work to visualize the biomechanical and temporal parameters of gait. Consider, for exemplification, the surface plot of the correlation coefficient matrix illustrated in Figure 24. Each 10 × 10 matrix section, delimited by the squares, displays 100 correlation coefficients for the shifted versions of the signal pairs indicated in the matrix header. Yellow corresponds to large coefficient values, i.e., good similarity between the signals, whereas blue corresponds to small coefficient values. As such, the largest correlation coefficient values, i.e., the brightest colors in the correlation coefficient matrix surface plot, yield the displacement between the signals given in the matrix header. As shown, the displacement determined from the correlation coefficient matrix is equal to the one determined using the cross-correlation function in Figure 19.
For illustration, the surface plot of the correlation coefficient matrix computed for the healthy control and the three PD patients from Figure 13 is illustrated in Figure 25. Interpretation of the correlation coefficient matrices follows in Section 4—Discussion.

3.3. Gait Assessment—AI-Based Decisional Support

AI-based decisional support for the discrimination of the PD-specific gait pattern was implemented using CNNs. The CNN takes the jpeg of the correlation coefficient matrix surface plot to classify the gait into physiological or PD pattern.
The CNN model performances were evaluated based on accuracy (Acc), sensitivity (Se), specificity (Sp), and precision, i.e., positive predicted values (PPV). These are defined using the true positive (TP), false positive (FP), true negative (TN), and false negative (FN) counts, respectively, as indicated in Equations (13)–(16). TP and TN count the number of correct classifications of Parkinsonian and physiological walking, respectively, whereas FP and FN count the number of incorrect classifications. According to Equations (13)–(16), Acc shows the probability of correct classification, Se evaluates the model’s ability to identify the true positive samples, Sp evaluates the model’s ability to identify the true negative samples, and PPV indicates the probability of the model to correctly classify a sample as positive. Finally, the error is then computed as the difference between the model training result and the desired result.
A c c = T N + T P T P + T N + F P + F N ,
S e = T P T P + F N ,  
S p = T N T N + F P ,
P P V = T P T P + F P
At the end of the training and data validation process, the performance metrics were calculated and listed in Table 4, suggesting that this algorithm for training a CNN network achieves efficient identification. As indicated, the best results were obtained on the MobileNet model with 95% accuracy, 90% sensitivity, and 95% precision. The best results were obtained on the MobileNet model.
The best performance metrics which we achieved in this work are listed in Table 5 for further comparison with the classification performances achieved by some classifiers reported in the literature, including Random Forest and Support Vector Machines (SVM).

4. Discussion

The work described in this article is developed in the context of the PDxOne research project, which aims to implement AI-based support for the collection and interpretation of medical data in PD, in the framework of ubiquitous healthcare.
The application described in this paper targets gait monitoring and assessment, based on the foot biomechanics, i.e., plantar pressure and lower-limb EMG, in correlation to upper-limb balance. To evaluate the gait problems that characterize PD, clinicians use semiquantitative rating scales such as the unified Parkinson’s disease rating scale (UPDRS) [66] or the movement disorders society unified Parkinson’s disease rating scale (MDS-UPDRS) [42]. Objective gait evaluation was performed in this work accordingly. For qualitative gait assessment, we have evaluated biomechanical parameters expressed in terms of arm balance, heel strike, and foot lift-off, with the results listed in Table 2. For quantitative gait assessment, we have evaluated temporal parameters expressed in terms of cadence, single support, double support, single support to double support ratio, and stride time variability, with the results listed in Table 3.
Gait impairment is evolving throughout the progression of the disease, and the patterns of gait disturbances that are detected can differ from early to mild/moderate and advanced stages of PD [67], but the relationship between gait features and disease progression is not completely explained [68].
In PD, non-motor symptoms such as anxiety, depression, and cognitive impairment develop along with the motor symptoms, influencing the subjects’ ability to perform motor tasks [9]. In the presence of such non-motor symptoms, motor tasks performed under trial conditions with a device attached to the body become a real challenge to the patient. Indeed, having the physiograph modules attached to the body produces an unusual sensorial stimulation to the patient. On the other hand, the presentation to the doctor’s office or the medical laboratory is a stress factor itself for many patients, which strikes the emotional component. Consequently, we expected the gait analysis results of the study group to be influenced: smaller stride length, slower gait velocity, and smaller activity motor unit recruitment, although gait is an activity of daily living. Indeed, the results reported in Table 3 illustrate that both PD patients and healthy controls exhibit a smaller cadence compared to the nominal 60 steps/min [43,52,53]. As a natural consequence of the reduced cadence, Table 3 also shows longer single-support and shorter double-support durations compared to the nominal 30% of the gait cycle [43,52,53].
Gait physiology was further presented in a visual manner using the correlation coefficient matrix defined in Section 2.2. Identification of the gait pattern on the correlation coefficient matrix is discussed as follows: physiological gait, as monitored with the proposed physiograph and illustrated in the correlation coefficient matrix surface plot from Figure 25a, accounts for a heel→metatarsal arch→hallux plantar pressure progression pattern in direct correlation to the lower-limb muscular activation. The plantar pressure progression pattern is visualized in the plantar pressure correlation sections by the shift of the in-section yellow diagonal toward the matrix main diagonal. Shifting of the in-section diagonal away from the matrix main diagonal would account for an inverse plantar pressure progression pattern, namely hallux→metatarsal arch→heel.
Activation of the TA accounts for eccentric contraction during heel strike and initial double limb support [43,52,69], which is visualized by the yellow diagonal in the TA-FSR correlation sections, for both left and right feet, respectively. Next, activation of the TA accounts for concentric contraction during the swing phase [43,52,69], which is visualized by the yellow areas in the top left and bottom right, i.e., dark blue diagonal, in the TA-FSR correlation sections for opposing feet.
Activation of the GM describes eccentric contraction during midstance and concentric contraction during heel off and toe off [43,52,69], which is visualized in the GM-FSR correlation sections for both left and right foot, respectively. The GM is inactive during the swing phase, which is visualized in the GM-FSR correlation sections for opposing feet.
Physiological gait assumes a complete arm balance during the gait cycle, with the forward swing accounting for right foot stance and the backward swing accounting for left foot stance. As such, the arm balance is fully correlated to the lower-limb EMG and plantar pressure signals, as visualized in the MA-EMG and MA-FSR correlation sections. To be noted is that only the right arm was considered for assessment in the present work, which explains the yellow diagonal in the correlation sections of the MA with the opposite, i.e., left foot, and the yellow corner areas (dark blue diagonal) in the correlation sections of the MA with the same side, i.e., right foot. Indeed, the forward swing of the right arm−left foot stance produces a larger magnitude of the U-shaped MA signal, in contrast to the backward arm swing−right foot stance, which produces a smaller magnitude (see Figure 14). The VA signal on the other hand exhibits peaks during midstance, thus being correlated to the MA, which explains the yellow diagonal on the MA-VA correlation section. Furthermore, the VA is correlated to shifted versions of the EMG and FSR signals, which explains the yellow areas in the VA-EMG and VA-FSR correlation sections.
Parkinsonian gait is clearly distinguishable from the physiological gait. One of the most representative but non-specific early features of Parkinsonian gait is reduced speed [67]. It has been demonstrated that early PD subjects exhibit a reduced amplitude of arm swing and smoothness of locomotion, as well as increased interlimb asymmetry, all of these being more specific to PD and often the first motor symptoms [70].
Such features of the Parkinsonian gait pattern are identifiable on the correlation coefficient matrix. Some gait features attributable to PD, which we were able to identify and assess during the clinical test of the proposed physiograph, are presented as follows: the correlation coefficient matrix surface plot from Figure 25b corresponds to Patient 5, whose walking pattern described in Table 2 consists of flat-foot strike, bilateral lift-off, presence of arm balance, and small tremor. In contrast to the healthy control, the flat-foot strike is visualized in the left-foot FSR−FSR correlation sections as a yellow main diagonal. Plantar pressure is applied simultaneously to all sensors during flat foot strike, and consequently the FSR signals are correlated to one another (see Figure 13b). The patient keeps a regular lower-limb muscular response during the stance phases, visualized by the TA-FSR and GM-FSR correlation sections. The arm balance is also present, visualized by the MA-EMG and MA-FSR correlation sections. Tremor, although existent, is small in magnitude and consequently allows for the visualization of the yellow diagonals in the VA-EMG and VA-FSR correlation sections.
The correlation coefficient matrix surface plot from Figure 25c corresponds to Patient 1, whose walking pattern described in Table 2 exhibits bilateral flat-foot strike, bilateral lift-off, presence of arm balance, and large tremor. The correlation coefficient matrix follows the same pattern regarding foot biomechanics as for Patient 5. Regarding tremor, however, Figure 13c shows that the tremor and balance magnitudes in the MA and VA signals are comparable. Consequently, both MA and VA signals are uncorrelated to the EMG and FSR signals, respectively. This is visualized by the correlation sections of both MA and VA to the other monitoring sections, which exhibit a rather uniform coloring.
A different walking pattern was identified for Patient 3, with the correlation coefficient matrix surface plot from Figure 25d, who according to Table 2 exhibits absent heel strike and absent lift-off for the right foot. In this regard, the patient pulls the right foot during the swing phase of the gait cycle. This is visualized on the correlation coefficient plot by the right-foot FSR−FSR correlation sections and the left FSR−right FSR correlation sections, which deviate from the yellow diagonal pattern.
As the disease progresses bilaterally, the asymmetry might decrease, and movement becomes more bradykinetic [68,71]. At the same time, along with the neurodegeneration, the movement of the limbs becomes more impaired, and the patients develop shuffling steps with the increased need of double-limb support [72,73]. The further decline in gait is also caused by the postural changes, which are altering the kinematics of the gait, as is the case of the stooped posture [74].
The changes in gait worsen and motor fluctuations, dyskinesias, and freezing of gait become frequent and are accompanied by reduced balance and postural control, all of these exposing the patient to a severe risk of falling. More than that, the decline in the motor capacity can lead, in advanced stages, to the need of wheelchair use [75,76,77]. We expect to see these changes in the correlation coefficient matrix to facilitate gait pattern interpretation. This could be a game changer for stage-related personalized medicine therapy options.
The proposed physiograph is also envisioned as a portable monitoring solution. For this purpose, the recorded signals are transmitted over a Bluetooth radio link to a mobile device, e.g., smartphone or tablet, and stored in an online database for future retrieval. As such, the proposed solution is applicable for gait monitoring in PD outside the hospital environment. In this scenario, the monitoring protocol can be extended to include the assessment and quantification of influences exercised by daily activities—both domestic and in the community [36], environmental demands [78], environmental manipulation [37], and dual/multiple tasking [36,37].
The proposed physiograph is in a proof-of-concept phase, which was validated in clinical environment. As such, it was to be expected that the performance of both patients and healthy controls was influenced as reported, due to their subjective perception on the worn-in discomfort. There are two reasons for this. Firstly, the proposed physiograph is still in proof-of-concept phase. In accordance with this, the subjects complained of worn-in discomfort and discomfort from wires interfering with their mobility. Secondly, the gait measurements were conducted in a relatively small room, which likely inhibited the automaticity of walking at a normal pace [79]. The clinical environment created conditions for directed attention during the walking trial. As such, the subjects are aware of the motor task and concentrate on performing the task, which may have led to the abovementioned deviations from standard.
The solution to overcome such issues is the repetition of the motor task for a prolonged period (several days) as well as removing the patient from the clinical environment and placing them into their usual living environment, both inside and outside the home. This is enabled by the wearability and portability features of the proposed physiograph.
In continuation of our research, we aim to optimize the system design toward miniaturization, targeting a wireless body area network topology. One step further, the external environment brings multiple sensory information that can divert the patient’s attention from the motor task. Thus, monitoring in an external environment provides multiple functional pieces of information with high value in the selection of the therapeutic and recovery interventions.
The solution proposed in this article was validated in a clinical environment for specific applicability in PD. Nevertheless, we envision that it can be easily extrapolated to further neurodegenerative conditions with gait affection.
On the other hand, PD does not exclude the presence of age-specific concurrent diseases such as (degenerative) osteoarthritis and/or cardiovascular diseases and their sequelae; comorbidities are common among patients with PD [80]. For example, one of the study group patients (Patient 3) also presents post-AVC hemiplegia with mild upper arm sequelae. In addition, inherent age-induced changes which affect joint function, e.g., sarcopenia, loss of proprioception and balance, and increased joint laxity, cannot be neglected [80]. All are associated with loss of independence and an increased incidence of falls.
Along with the clinical, biological, and imaging exam, the functional exam is also very important (including the functional exam performed with our device) for diagnostic purposes, as well as for monitoring the pharmacological disease management and rehabilitation outcome. It can be also used to guide kinetic trainings, to impose gait rhythms and velocities, and even to aid in dual tasks.
Along with the paraclinical examinations, the clinical testing, including functional performance, is key for diagnostic purposes. Part of this examination could be performed with our device and could be useful to provide guidance in kinetic trainings, to impose gait rhythms and velocities, and to aid in dual tasks. With further improvements, the pharmacological monitoring, disease management, and rehabilitation outcomes are in reach.

5. Conclusions

PD clearly demonstrates the complexity of the human body and the difficulty in choosing the appropriate treatment. In this paper, we focused on one of the main motor symptoms in PD, the gait. To develop a foundation for sustainable decision-support medical devices, we aim to include other important motor symptoms, as well as the non-motor ones, in future studies. Non-motor symptoms (e.g., cognitive impairment, depression, anxiety, sleep disorders, pain, and other autonomic disturbances) correlate with advanced age and disease severity [81], so they can be considered suggestive for the prognosis of disease progression. Some of these non-motor symptoms may appear much earlier in the course of disease progression [12]. AI-based assessment of these data could raise the importance of new technologies in prediagnosis and digital monitoring of patients. By pre-diagnostic technologies, we refer to those technologies that have the probabilistic ability to provide data interpretation at a high level of certitude for a given diagnosis, often before it is identified by a specialist without computational aid. In our experience, an additional interesting feature of AI systems is that they can identify non-specific symptoms, provided enough data and solid input are available. Regardless of the future development and promising results, the final diagnosis should be reserved, in our opinion, exclusively for the human medical practitioner, following established ethical/legal rules and evidence-based practice.
Back in 2017 when the PDxOne started, the main background idea was to find solutions for PD management optimization by bridging the gap between medicine and computational engineering. We used this approach in our study. Taking into consideration everything mentioned above, we see the necessity for an increased commitment in multidisciplinary treatment by interdisciplinary teams. In our opinion, this interdisciplinary team should include, at least, the following specialists in addition to the neurologist: general practitioner and dentist, geriatrician, trained medical caregiver for patients with PD, speech therapist, occupational therapist, nutritionist, psychologist/psychotherapist ideally specialized in patients with PD, and pharmacologists. To the abovementioned list, we would add new kinds of professionals that have medical degrees such as in a biomedical field (e.g., bioengineer, bioprinting expert, etc.) combined with knowledge of AI and augmented intelligence, as well as lifestyle strategists (advise patients with their health data), telemedicine, and/or health data analyst/biostatistician. Considering that by 2030 we will lack over 18 million healthcare workers worldwide [82], these future professions will be an asset in bridging interdisciplinary activities, they will provide sustainability in healthcare, and they could serve as good choice for internationally trained medical professionals toward alternative career pathways [83] because of labor migration in the era of globalization. We see today this new professional category in its early phase, with clinics having incorporated bioengineers, deep learning developers, bioprinting experts, and so on on their teams. New research projects are already hard to imagine without the congregation of such a team. As with most advancements, this one will also be led by necessity. In the time of Big Data, such an interdisciplinary team could provide personalized treatment for hospitalized patients and those at home. Perhaps the most important aspect of this interdisciplinary team is the medical data, and how interdisciplinary teams could correlate assessments to identify common patterns for optimizing patient-specific treatments. A good line of action would be the development of an international research platform, with standardized parameters, for automated or semiautomated data input. The assessment algorithms for this platform would benefit from the interdisciplinary approach and the fact that coding and treatment could be performed “at the patient’s bed”. A possible output would be an evidence-based decision support software for optimal dosing between drug and non-drug treatments. With such an effort, adequate treatment could be provided for this complex and multi-faceted disease, as PD is scripted [18]. Such an approach, once established, could offer individualized treatment plans and monitoring programs that would further improve PD management, and with it, socio-economic implications. The ultimate goal of our research community should be to offer affordable precision medicine for everyone.

Author Contributions

Conceptualization, R.R.I. and P.F.; methodology, R.R.I. and P.F.; software, C.-G.C., L.-I.M., R.F. and P.F.; validation, R.R.I., L.P.-D. and P.F.; formal analysis, P.F.; investigation, R.R.I. and P.F.; resources, R.R.I., A.-S.P. and L.P.-D.; data curation, P.F.; writing—original draft preparation, R.R.I., C.-G.C., L.-I.M., R.F., A.-S.P. and P.F.; writing—review and editing, R.R.I., L.P.-D. and P.F.; visualization, C.-G.C., L.-I.M., R.F. and P.F.; supervision, R.R.I. and P.F.; project administration, P.F.; funding acquisition, P.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by ELYSEUM S.R.L., grant number 21738/30.07.2021.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the University of Medicine and Pharmacy “Iuliu Hatieganu” Cluj-Napoca, Romania (protocol code 86 and date of approval 1 February 2018).

Informed Consent Statement

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

Acknowledgments

The authors would like to thank Florian Thieringer from the Clinic of Oral and Cranio-Maxillofacial Surgery, University Hospital Basel, Basel, Switzerland, for his valuable support in the final review of the paper. The authors would also like to thank Martin Douglas-Jones from the Department of Maxillofacial and Oral Surgery, Frere Provincial Hospital, Eastern Cape, South Africa, for his valuable support in the final review of the paper. The authors would also like to thank Monica-Adriana Farago from the Psycho-Neurosciences and Rehabilitation Department, Faculty of Medicine and Pharmacy, University of Oradea, Oradea, Romania, for her valuable support in the final review of the paper.

Conflicts of Interest

Ileșan Robert Radu is also CEO of Elyseum SRL, the company which funded the grant for the present research. As a physician, maxillofacial surgeon, researcher, and professional, he declares that the outcome of this paper was not influenced from a conflict-of-interest point of view. We hereby declare that the conduit of the research activity is straight and follows the international educational and research values, Dr. med. Dr. med. dent. Ileșan Robert Radu.

References

  1. Parkinson, J. An essay on the shaking palsy. 1817. J. Neuropsychiatry Clin. Neurosci. 2002, 14, 223–236. [Google Scholar] [CrossRef] [PubMed]
  2. Dorsey, E.R.; Sherer, T.; Okun, M.S.; Bloem, B.R. The Emerging Evidence of the Parkinson Pandemic. J. Parkinson’s Dis. 2018, 8, S3–S8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Svehlík, M.; Zwick, E.B.; Steinwender, G.; Linhart, W.E.; Schwingenschuh, P.; Katschnig, P.; Ott, E.; Enzinger, C. Gait analysis in patients with Parkinson’s disease off dopaminergic therapy. Arch. Phys. Med. Rehabil. 2009, 90, 1880–1886. [Google Scholar] [CrossRef] [PubMed]
  4. De Lau, L.M.; Breteler, M.M. Epidemiology of Parkinson’s disease. Lancet. Neurol. 2006, 5, 525–535. [Google Scholar] [CrossRef]
  5. Lipton, R.B.; Schwedt, T.J.; Friedman, B.W. GBD 2015 Disease and Injury Incidence and Prevalence, Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990–2015: A systematic analysis for the Global Burden of Disease Study 2015. Lancet 2016, 388, 1545–1602. [Google Scholar]
  6. Samii, A.; Nutt, J.G.; Ransom, B.R. Parkinson’s Disease. Lancet 2004, 363, 1783–1793. [Google Scholar] [CrossRef] [Green Version]
  7. Kreutzer, J.S.; John, D.; Bruce, C. Encyclopedia of Clinical Neuropsychology; Springer: New York, NY, USA, 2011; pp. 1380–1383. [Google Scholar]
  8. Thenganatt, M.A.; Jankovic, J. Parkinson disease subtypes. JAMA Neurol. 2014, 71, 499–504. [Google Scholar] [CrossRef]
  9. Jankovic, J. Parkinson’s disease: Clinical features and diagnosis. J. Neurol. Neurosurg. Psychiatry 2008, 79, 368–376. [Google Scholar] [CrossRef] [Green Version]
  10. Triarhou, L.C. Dopamine and Parkinson’s Disease. In Madame Curie Bioscience Database; Landes Bioscience: Austin, TX, USA, 2013. Available online: https://www.ncbi.nlm.nih.gov/books/NBK6271/ (accessed on 19 March 2022).
  11. Postuma, R.B.; Berg, D.; Stern, M.; Poewe, W.; Olanow, C.W.; Oertel, W.; Obeso, J.; Marek, K.; Litvan, I.; Lang, A.E.; et al. MDS clinical diagnostic criteria for Parkinson’s disease. Mov. Disord. Off. J. Mov. Disord. Soc. 2015, 30, 1591–1601. [Google Scholar] [CrossRef]
  12. Chaudhuri, K.R.; Healy, D.G.; Schapira, A.H. Non-motor symptoms of Parkinson’s disease: Diagnosis and management. Lancet. Neurol. 2006, 5, 235–245. [Google Scholar] [CrossRef]
  13. Mantri, S.; Morley, J.F.; Siderowf, A.D. The importance of preclinical diagnostics in Parkinson disease. Parkinsonism Relat. Disord. 2019, 64, 20–28. [Google Scholar] [CrossRef] [PubMed]
  14. Postuma, R.B.; Aarsland, D.; Barone, P.; Burn, D.J.; Hawkes, C.H.; Oertel, W.; Ziemssen, T. Identifying prodromal Parkinson’s disease: Pre-motor disorders in Parkinson’s disease. Mov. Disord. Off. J. Mov. Disord. Soc. 2012, 27, 617–626. [Google Scholar] [CrossRef] [PubMed]
  15. Dorsey, E.R.; Constantinescu, R.; Thompson, J.P.; Biglan, K.M.; Holloway, R.G.; Kieburtz, K.; Marshall, F.J.; Ravina, B.M.; Schifitto, G.; Siderowf, A.; et al. Projected number of people with Parkinson disease in the most populous nations, 2005 through 2030. Neurology 2007, 68, 384–386. [Google Scholar] [CrossRef] [PubMed]
  16. Tickle, D.L.; Zebrowitz, L.A.; Ma, H.I. Culture, gender and health care stigma: Practitioners’ response to facial masking experienced by people with Parkinson’s disease. Soc. Sci. Med. 2011, 73, 95–102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Poewe, W. The natural history of Parkinson’s disease. J. Neurol. 2006, 253, VII2–VII6. [Google Scholar] [CrossRef] [PubMed]
  18. Gustavsson, A. Cost of disorders of the brain in Europe 2010. Eur. Neuropsychopharmacol. 2011, 21, 718–779. [Google Scholar] [CrossRef] [Green Version]
  19. Kowal, S.L.; Stacey, L. The current and projected economic burden of Parkinson’s disease in the United States. Mov. Disord. 2013, 28, 311–318. [Google Scholar] [CrossRef]
  20. Winter, Y. Longitudinal study of the socioeconomic burden of Parkinson’s disease in Germany. Eur. J. Neurol. 2010, 17, 1156–1163. [Google Scholar] [CrossRef]
  21. Kosten von Parkinson Nach Land Sehr Unterschiedlich. Available online: https://healthcare-in-europe.com/de/news/kosten-von-parkinson-nach-land-sehr-unterschiedlich.html (accessed on 24 June 2010).
  22. Findley, L.J. The economic burden of advanced Parkinson’s disease: An analysis of a UK patient dataset. J. Med. Econ. 2011, 14, 130–139. [Google Scholar] [CrossRef]
  23. Findley, L. Direct economic impact of Parkinson’s disease: A research survey in the United Kingdom. Mov. Disord. 2003, 18, 1139–1145. [Google Scholar] [CrossRef]
  24. Boroojerdi, B.; Ghaffari, R.; Mahadevan, N.; Markowitz, M.; Melton, K.; Morey, B.; Sheth, N. Clinical feasibility of a wearable, conformable sensor patch to monitor motor symptoms in Parkinson’s disease. Parkinsonism Relat. Disord. 2019, 61, 70–76. [Google Scholar] [CrossRef] [PubMed]
  25. Jauhiainen, M.; Puustinen, J.; Mehrang, S.; Ruokolainen, J.; Holm, A.; Vehkaoja, A.; Nieminen, H. Identification of motor symptoms related to Parkinson disease using motion-tracking sensors at home (KÄVELI): Protocol for an observational case-control study. JMIR Res. Protoc. 2019, 8, 12808. [Google Scholar] [CrossRef] [PubMed]
  26. Phan, D.; Horne, M.; Pathirana, P.N.; Farzanehfar, P. Measurement of axial rigidity and postural instability using wearable sensors. Sensors 2018, 18, 495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lonini, L.; Dai, A.; Shawen, N.; Simuni, T.; Poon, C.; Shimanovich, L.; Jayaraman, A. Wearable sensors for Parkinson’s disease: Which data are worth collecting for training symptom detection models. NPJ Digit. Med. 2018, 1, 64. [Google Scholar] [CrossRef] [PubMed]
  28. Luigi, B.; Marilena, V.; Gabriella, O.; Carlo, A.A.; Margherita, F.; Mario, G.R.; Alberto, R.; Maurizio, Z.; Leonardo, L. Home monitoring of motor fluctuations in Parkinson’s disease patients. J. Reliab. Intell. Environ. 2019, 5, 145–162. [Google Scholar]
  29. Rob, P.; Maryam, E.A.; Edith, M.A.; Sara, K.; Irida, M. Smartwatch inertial sensors continuously monitor real-world motor fluctuations in Parkinson’s disease. Sci. Transl. Med. 2021, 13, 7865. [Google Scholar] [CrossRef]
  30. Margot, H.; Jeroen, G.V.H.; Christian, H.; Jos, A.; An, S.; Mark, L.K.; Pieter, L.K. Monitoring Parkinson’s disease symptoms during daily life: A feasibility study. NPJ Parkinson’s Dis. 2019, 5, 21. [Google Scholar]
  31. Francesco, M.; Alessandro, M.; Fiorella, P.; Mariagrazia, R.; Fioravante, C.; Keren, K.; Alit, S.I.; Ziv, Y.; Vincenzo, D.L.; Massimo, M. Parkinson’s Disease Remote Patient Monitoring During the COVID-19 Lockdown. Front. Neurol. 2020, 11, 567413. [Google Scholar]
  32. Satyabrata, A.; Jinyoung, Y.; Sabyasachi, C.; Pyari, M.P. A Supervised Machine Learning Approach to Detect the on/off State in Parkinson’s Disease Using Wearable Based Gait Signals. Diagnostics 2020, 10, 421. [Google Scholar] [CrossRef]
  33. Um, T.T.; Terry, T. Parkinson’s disease assessment from a wrist-worn wearable sensor in free-living conditions: Deep ensemble learning and visualization. arXiv 2018, arXiv:1808.02870. [Google Scholar]
  34. Tobias, S.; Michele, M.; Ingrid, B.; Carlos, M.T. Analysis and Classification of Motor Dysfunctions in Arm Swing in Parkinson’s Disease. Electronics 2019, 8, 1471. [Google Scholar] [CrossRef] [Green Version]
  35. Pirker, W.; Katzenschlager, R. Gait disorders in adults and the elderly: A clinical guide. Wien. Klin. Wochenschr. 2017, 129, 81–95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Muthukrishnan, N.; Abbas, J.J.; Shill, H.A.; Krishnamurthi, N. Cueing Paradigms to Improve Gait and Posture in Parkinson’s Disease: A Narrative Review. Sensors 2019, 19, 5468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Ebersbach, G.; Moreau, C.; Gandor, F.; Defebvre, L.; Devos, D. Clinical syndromes: Parkinsonian gait. Mov. Disord. 2013, 28, 1552–1559. [Google Scholar] [CrossRef]
  38. Seong, B.C.; Kun, W.P.; Dae, H.L.; Se, J.K.; Joon, S.Y. Gait Analysis in Patients with Parkinson’s Disease: Relationship to Clinical Features and Freezing. J. Mov. Disord. 2008, 1, 59–64. [Google Scholar]
  39. Sofuwa, O.; Nieuwboer, A.; Desloovere, K.; Willems, A.M.; Chavret, F.; Jonkers, I. Quantitative gait analysis in Parkinson’s disease: Comparison with a healthy control group. Arch. Phys. Med. Rehabil. 2005, 86, 1007–1020. [Google Scholar] [CrossRef] [Green Version]
  40. Bovolenta, T.M.; de Azevedo Silva, S.; Saba, R.A.; Borges, V.; Ferraz, H.B.; Felicio, A.C. Average annual cost of Parkinson’s disease in São Paulo, Brazil, with a focus on disease-related motor symptoms. Clin. Interv. Aging 2017, 12, 2095–2108. [Google Scholar] [CrossRef] [Green Version]
  41. Hauser, R.A.; Lyons, K.E.; Pahwa, R. The UPDRS-8: A brief clinical assessment scale for Parkinson’s disease. Int. J. Neurosci. 2012, 122, 333–340. [Google Scholar] [CrossRef] [Green Version]
  42. Goetz, C.G.; Tilley, B.C.; Shaftman, S.R.; Stebbins, G.T.; Fahn, S.; Martinez-Martin, P.; Poewe, W.; Sampaio, C.; Stern, M.B.; Dodel, R. Movement Disorder Society-sponsored revision of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS): Scale presentation and clinimetric testing results. Mov. Disord. 2008, 23, 2129–2170. [Google Scholar] [CrossRef]
  43. Faragó, P.; Grama, L.; Farago, M.A.; Hintea, S. A Novel Wearable Foot and Ankle Monitoring System for the Assessment of Gait Biomechanics. Appl. Sci. 2021, 11, 268. [Google Scholar] [CrossRef]
  44. Narwaria, R.P.; Verma, S.; Singhal, P.K. Removal of Baseline Wander and Power Line Interference from ECG Signal—A Survey Approach. Int. J. Electron. Eng. 2011, 3, 107–111. [Google Scholar]
  45. Hung, C.C.; Halonen, I.; Ismail, M.; Porra, V. Micropower CMOS Gm-C Filters for Speech Signal Processing. In Proceedings of the 1997 IEEE International Symposium on Circuits and Systems, Hong Kong, China, 9–12 June 1997. [Google Scholar]
  46. Mills, K.R. The basics of electromyography. J. Neurol. Neurosurg. Psychiatry 2005, 76, 32–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Faragό, P.; Groza, R.; Hintea, S. High Precision Activity Tracker Based on the Correlation of Accelerometer and EMG Data. In Proceedings of the 42nd International Conference on Telecommunications and Signal Processing (TSP), Budapest, Hungary, 1–3 July 2019. [Google Scholar]
  48. San-Segundo, R.; Zhang, A.; Cebulla, A. Parkinson’s Disease Tremor Detection in the Wild Using Wearable Accelerometers. Sensors 2020, 20, 5817. [Google Scholar] [CrossRef] [PubMed]
  49. Fridolfsson, J.; Börjesson, M.; Arvidsson, D. A Biomechanical Re-Examination of Physical Activity Measurement with Accelerometers. Sensors 2018, 18, 3399. [Google Scholar] [CrossRef] [Green Version]
  50. Lee, H.J.; Lee, W.W.; Kim, S.K.; Park, H.; Jeon, H.S.; Kim, H.B.; Park, K.S. Tremor frequency characteristics in Parkinson’s disease under resting-state and stress-state conditions. J. Neurol. Sci. 2016, 362, 272–277. [Google Scholar] [CrossRef]
  51. Lin, W.Y.; Verma, V.K.; Lee, M.Y.; Lai, C.S. Activity Monitoring with a Wrist-Worn, Accelerometer-Based Device. Micromachines 2018, 9, 450. [Google Scholar] [CrossRef] [Green Version]
  52. Karadsheh, M. Orthobullets—Gait Cycle. Available online: https://www.orthobullets.com/foot-and-ankle/7001/gait-cycle (accessed on 23 December 2021).
  53. Alamdari, A.; Krovi, V.N. A Review of Computational Musculoskeletal Analysis of Human Lower Extremities. In Human Modeling for Bio-Inspired Robotics: Mechanical Engineering in Assistive Technologies; Ueda, J., Kurita, Y., Eds.; Academic Press: San Diego, CA, USA, 2017; pp. 37–73. [Google Scholar]
  54. Hughes, J.R.; Bowes, S.G.; Leeman, A.L.; O’Neill, C.J.; Deshmukh, A.A.; Nicholson, P.W. Parkinsonian abnormality of foot strike: A phenomenon of ageing and/or one responsive to levodopa therapy? Br. J. Clin. Pharmacol. 1990, 29, 179–186. [Google Scholar] [CrossRef] [Green Version]
  55. 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]
  56. O’Sullivan, S.O. Parkinson’s Disease: Physical Therapy Intervention. In Physical Rehabilitation, 5th ed.; O’Sullivan, S.B., Schmitz, T.J., Eds.; EA Davis Company: Philadelphia, PA, USA, 2007; pp. 853–893. [Google Scholar]
  57. Marple, S.L., Jr. Digital Spectral Analysis, 2nd ed.; Dover Publications: Mineola, NY, USA, 2019; ISBN 978-048-6780-52-8. [Google Scholar]
  58. Park, Y.; Chung, W. A Novel EEG Correlation Coefficient Feature Extraction Approach Based on Demixing EEG Channel Pairs for Cognitive Task Classification. IEEE Access 2020, 8, 87422–87433. [Google Scholar] [CrossRef]
  59. Diwaker, S.; Gupta, S.C.; Gupta, N.K. Classification of EEG Signal using Correlation Coefficient among Channels as Features Extraction Method. Indian J. Sci. Technol. 2016, 9, 1–7. [Google Scholar] [CrossRef]
  60. Howard, A.G.; Zhu, M.; Chen, B.; Kalenichenko, D.; Wang, W.; Weyand, T.; Adam, H. Mobilenets: Efficient convolutional neural networks for mobile vision applications. arXiv 2017, arXiv:1704.04861. [Google Scholar]
  61. Tan, M.; Le, Q. EfficientNet: Rethinking Model Scaling for Convolutional Neural Networks. In Proceedings of the 36th International Conference on Machine Learning, Long Beach, CA, USA, 9–15 June 2019; pp. 6105–6114. [Google Scholar]
  62. Moon, Y.; Wajda, D.A.; Motl, R.W.; Sosnoff, J.J. Stride-Time Variability and Fall Risk in Persons with Multiple Sclerosis. Mult. Scler. Int. 2015, 2015, 964790. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Beauchet, O.; Annweiler, C.; Lecordroch, Y. Walking speed-related changes in stride time variability: Effects of decreased speed. J. Neuroeng. Rehabil. 2009, 6, 32. [Google Scholar] [CrossRef]
  64. Beauchet, O.; Freiberger, E.; Annweiler, C.; Kressig, R.W.; Herrmann, F.R.; Allali, G. Test-retest reliability of stride time variability while dual tasking in healthy and demented adults with frontotemporal degeneration. J. Neuroeng. Rehabil. 2011, 8, 37. [Google Scholar] [CrossRef] [PubMed]
  65. Jeon, H.-S.; Han, J.; Yi, W.-J.; Jeon, B.; Park, K.S. Classification of Parkinson gait and normal gait using Spatial-Temporal Image of Plantar pressure. In Proceedings of the 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Vancouver, BC, Canada, 20–25 August 2008; pp. 4672–4675. [Google Scholar] [CrossRef]
  66. Fahn, S.; Elton, R.L. Recent Developments in Parkinson’s Disease; Fahn, S., Marsden, C.D., Goldstein, M., Calne, D.B., Eds.; Macmillan Health Care Information: Florham Park, NJ, USA, 1987; pp. 153–163. [Google Scholar]
  67. Mirelman, A.; Bonato, P.; Camicioli, R.; Ellis, T.D.; Giladi, N.; Hamilton, J.L.; Hass, C.J.; Hausdorff, J.M.; Pelosin, E.; Almeida, Q.J. Gait impairments in Parkinson’s disease. Lancet Neurol. 2019, 18, 697–708. [Google Scholar] [CrossRef]
  68. Albani, G.; Cimolin, V.; Fasano, A.; Trotti, C.; Galli, M.; Mauro, A. “Masters and servants” in parkinsonian gait: A three-dimensional analysis of biomechanical changes sensitive to disease progression. Funct. Neurol. 2014, 29, 99–105. [Google Scholar] [PubMed]
  69. Honeine, J.L.; Schiepati, M.; Gagey, O.; Do, M.C. The Functional Role of the Triceps Surae Muscle during Human Locomotion. PLoS ONE 2013, 8, 52943. [Google Scholar] [CrossRef] [Green Version]
  70. Sterling, N.W.; Cusumano, J.P.; Shaham, N. Dopaminergic modulation of arm swing during gait among Parkinson’s disease patients. J. Parkinson’s Dis. 2015, 5, 141–191. [Google Scholar] [CrossRef] [Green Version]
  71. Galna, B.; Lord, S.; Burn, D.J.; Rochester, L. Progression of gait dysfunction in incident Parkinson’s disease: Impact of medication and phenotype. Mov. Disord. 2015, 30, 359–426. [Google Scholar] [CrossRef]
  72. Mirelman, A.; Bernad-Elazari, H.; Thaler, A. Arm swing as a potential new prodromal marker of Parkinson’s disease. Mov. Disord. 2016, 31, 1527–1561. [Google Scholar] [CrossRef]
  73. Kwon, K.Y.; Kim, M.; Lee, S.M.; Kang, S.H.; Lee, H.M.; Koh, S.B. Is reduced arm and leg swing in Parkinson’s disease associated with rigidity or bradykinesia? J. Neurol. Sci. 2014, 341, 32–35. [Google Scholar] [CrossRef] [PubMed]
  74. Pistacchi, M.; Gioulis, M.; Sanson, F. Gait analysis and clinical correlations in early Parkinson’s disease. Funct. Neurol. 2017, 32, 28–34. [Google Scholar] [CrossRef] [PubMed]
  75. Rezvanian, S.; Lockhart, T.; Frames, C.; Soangra, R.; Lieberman, A. Motor subtypes of Parkinson’s disease can be identified by frequency component of postural stability. Sensors 2018, 18, 1102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Okuma, Y.; Silva de Lima, A.L.; Fukae, J.; Bloem, B.R.; Snijders, A.H. A prospective study of falls in relation to freezing of gait and response fluctuations in Parkinson’s disease. Parkinsonism Relat. Disord. 2018, 46, 30–35. [Google Scholar] [CrossRef]
  77. Yusupov, E.; Chen, D.; Krishnamachari, B. Medication use and falls: Applying Beers criteria to medication review in Parkinson’s disease. SAGE Open Med. 2017, 5, 2050312117743673. [Google Scholar] [CrossRef] [Green Version]
  78. Barbosa, A.F.; Chen, J.; Freitag, F.; Valente, D.; Souza, C.O.; Voos, M.C.; Chien, H.F. Gait, posture and cognition in Parkinson’s disease. Dement. Neuropsychol. 2016, 10, 280–286. [Google Scholar] [CrossRef]
  79. Clark, D.J. Automaticity of walking: Functional significance, mechanisms, measurement and rehabilitation strategies. Front. Hum. Neurosci. 2015, 9, 246. [Google Scholar] [CrossRef] [Green Version]
  80. Shane Anderson, A.; Loeser, R.F. Why is osteoarthritis an age-related disease? Best Pract. Res. Clin. Rheumatol. 2010, 24, 15–26. [Google Scholar] [CrossRef] [Green Version]
  81. Zhong, L.L.; Song, Y.Q.; Cao, H.; Ju, K.J.; Yu, L. The non-motor symptoms of Parkinson’s disease of different motor types in early stage. Eur. Rev. Med. Pharmacol. Sci. 2017, 21, 5745–5750. [Google Scholar] [CrossRef]
  82. Limb, M. World will lack 18 million health workers by 2030 without adequate investment, warns un. BMJ 2016, 354, i5169. [Google Scholar] [CrossRef] [Green Version]
  83. Turin, T.C.; Chowdhury, N.; Ekpekurede, M. Alternative career pathways for international medical graduates towards job market integration: A literature review. Int. J. Med. Educ. 2021, 12, 45–63. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The proposed wearable physiograph for gait monitoring in PD: (a) block diagram and (b) practical realization.
Figure 1. The proposed wearable physiograph for gait monitoring in PD: (a) block diagram and (b) practical realization.
Biosensors 12 00189 g001aBiosensors 12 00189 g001b
Figure 2. Illustration of the waist-worn ATmega2560 microcontroller development board, attached to a Velcro strip, which reads the peripherals involved in the bilateral monitoring of foot biomechanics and unilateral tracking of the arm balance.
Figure 2. Illustration of the waist-worn ATmega2560 microcontroller development board, attached to a Velcro strip, which reads the peripherals involved in the bilateral monitoring of foot biomechanics and unilateral tracking of the arm balance.
Biosensors 12 00189 g002
Figure 3. Illustration of the resistive pressure sensor placement onto the shoe insoles, under the toe, metatarsal arch, and heel area, respectively, for bilateral tracking of the plantar pressure progression pattern along the center of pressure progression line.
Figure 3. Illustration of the resistive pressure sensor placement onto the shoe insoles, under the toe, metatarsal arch, and heel area, respectively, for bilateral tracking of the plantar pressure progression pattern along the center of pressure progression line.
Biosensors 12 00189 g003
Figure 4. Deployment of the resistive pressure sensors for operation as force sense resistors: (a) schematics, (b) mass vs. resistance conversion characteristics.
Figure 4. Deployment of the resistive pressure sensors for operation as force sense resistors: (a) schematics, (b) mass vs. resistance conversion characteristics.
Biosensors 12 00189 g004
Figure 5. Illustration of the electrode placement for EMG acquisition of the Tibialis anterior and Gastrocnemius medialis, with the reference electrodes placed onto the lateral and medial malleolus, respectively.
Figure 5. Illustration of the electrode placement for EMG acquisition of the Tibialis anterior and Gastrocnemius medialis, with the reference electrodes placed onto the lateral and medial malleolus, respectively.
Biosensors 12 00189 g005
Figure 6. Illustration of the attachment of the LSM9DS0 accelerometer to the user’s right wrist, for arm balance tracking during the gait cycle.
Figure 6. Illustration of the attachment of the LSM9DS0 accelerometer to the user’s right wrist, for arm balance tracking during the gait cycle.
Biosensors 12 00189 g006
Figure 7. Diagram of the implemented solution to have the signals recorded from the users wearing the gait monitoring physiograph transmitted to a mobile terminal over a Bluetooth radio link and stored into an online database for future retrieval.
Figure 7. Diagram of the implemented solution to have the signals recorded from the users wearing the gait monitoring physiograph transmitted to a mobile terminal over a Bluetooth radio link and stored into an online database for future retrieval.
Biosensors 12 00189 g007
Figure 8. Example of a message sent to the API via the mobile phone.
Figure 8. Example of a message sent to the API via the mobile phone.
Biosensors 12 00189 g008
Figure 9. The online database content, resembling the data acquired with gait monitoring physiograph.
Figure 9. The online database content, resembling the data acquired with gait monitoring physiograph.
Biosensors 12 00189 g009
Figure 10. Photographs taken during the trials of gait assessment.
Figure 10. Photographs taken during the trials of gait assessment.
Biosensors 12 00189 g010
Figure 11. Illustration of the procedure for stance identification: the stance phase lasts from the first to the last occurrence of plantar pressure, identified by having the FSR signals compared to an empirical threshold level.
Figure 11. Illustration of the procedure for stance identification: the stance phase lasts from the first to the last occurrence of plantar pressure, identified by having the FSR signals compared to an empirical threshold level.
Biosensors 12 00189 g011
Figure 12. Flowchart of the CNN employed for the classification of the gait pattern into physiological and PD classes.
Figure 12. Flowchart of the CNN employed for the classification of the gait pattern into physiological and PD classes.
Biosensors 12 00189 g012
Figure 13. Three-cycle gait section plot of the arm balance MA, bilateral EMG of the TA and GM, and FSR signals, with the red triangles indicating the MA maxima, acquired with the proposed physiograph during walking trials on (a) healthy control, (b) PD Patient 5, (c) PD Patient 1, and (d) PD Patient 3.
Figure 13. Three-cycle gait section plot of the arm balance MA, bilateral EMG of the TA and GM, and FSR signals, with the red triangles indicating the MA maxima, acquired with the proposed physiograph during walking trials on (a) healthy control, (b) PD Patient 5, (c) PD Patient 1, and (d) PD Patient 3.
Biosensors 12 00189 g013
Figure 14. The arm balance MA and VA waveforms during the walking trials, with the red triangles indicating the MA maxima. (a) Healthy control—the MA exhibits a U-shaped variation and the VA accounts for smaller variation, and (b) PD patient—the MA exhibits a larger number of local peaks with reduced amplitude and the VA accounts for larger variation.
Figure 14. The arm balance MA and VA waveforms during the walking trials, with the red triangles indicating the MA maxima. (a) Healthy control—the MA exhibits a U-shaped variation and the VA accounts for smaller variation, and (b) PD patient—the MA exhibits a larger number of local peaks with reduced amplitude and the VA accounts for larger variation.
Biosensors 12 00189 g014
Figure 15. The FSR waveforms (FSR2—heel, FSR1—metatarsal arch, FSR0—toe) during the walking trials for (a) healthy control—exhibits the physiological heel→metatarsal arch→toe pressure progression pattern, and (b) Patient 5—exhibits flat-foot strike as pressure is applied simultaneously on all three FSRs.
Figure 15. The FSR waveforms (FSR2—heel, FSR1—metatarsal arch, FSR0—toe) during the walking trials for (a) healthy control—exhibits the physiological heel→metatarsal arch→toe pressure progression pattern, and (b) Patient 5—exhibits flat-foot strike as pressure is applied simultaneously on all three FSRs.
Biosensors 12 00189 g015
Figure 16. The FSR waveforms (FSR2—heel, FSR1—metatarsal arch, FSR0—toe) during the walking trials for (a) healthy control—exhibits a section when no pressure is applied, thus indicating lift-off, and (b) Patient 3—exhibits no section when no pressure is applied, indicating that the patient does not lift the foot from the ground but rather pulls the foot during swing.
Figure 16. The FSR waveforms (FSR2—heel, FSR1—metatarsal arch, FSR0—toe) during the walking trials for (a) healthy control—exhibits a section when no pressure is applied, thus indicating lift-off, and (b) Patient 3—exhibits no section when no pressure is applied, indicating that the patient does not lift the foot from the ground but rather pulls the foot during swing.
Biosensors 12 00189 g016
Figure 17. One-stance gait section plot of the arm balance MA, EMG of the TA and GM, and FSR signals for a healthy control, with the red triangles indicating the MA maxima: (a) left foot and (b) right foot.
Figure 17. One-stance gait section plot of the arm balance MA, EMG of the TA and GM, and FSR signals for a healthy control, with the red triangles indicating the MA maxima: (a) left foot and (b) right foot.
Biosensors 12 00189 g017
Figure 18. Cross-correlation of the MA and EMG signals plotted in Figure 17 for the (a) left foot and (b) right foot.
Figure 18. Cross-correlation of the MA and EMG signals plotted in Figure 17 for the (a) left foot and (b) right foot.
Biosensors 12 00189 g018
Figure 19. Cross-correlation of the TA and FSR signals plotted in Figure 17 for the (a) left foot and (b) right foot.
Figure 19. Cross-correlation of the TA and FSR signals plotted in Figure 17 for the (a) left foot and (b) right foot.
Biosensors 12 00189 g019
Figure 20. Cross-correlation of the GM and FSR signals plotted in Figure 17 for the (a) left foot and (b) right foot.
Figure 20. Cross-correlation of the GM and FSR signals plotted in Figure 17 for the (a) left foot and (b) right foot.
Biosensors 12 00189 g020
Figure 21. One-stance gait section plot of the arm balance MA, EMG of the TA and GM, and FSR signals for a PD patient, with the red triangles indicating the MA maxima: (a) left foot and (b) right foot.
Figure 21. One-stance gait section plot of the arm balance MA, EMG of the TA and GM, and FSR signals for a PD patient, with the red triangles indicating the MA maxima: (a) left foot and (b) right foot.
Biosensors 12 00189 g021
Figure 22. Cross-correlation of the TA and FSR signals plotted in Figure 21 for the (a) left foot and (b) right foot.
Figure 22. Cross-correlation of the TA and FSR signals plotted in Figure 21 for the (a) left foot and (b) right foot.
Biosensors 12 00189 g022
Figure 23. Cross-correlation of the GM and FSR signals plotted in Figure 21 for the (a) left foot and (b) right foot.
Figure 23. Cross-correlation of the GM and FSR signals plotted in Figure 21 for the (a) left foot and (b) right foot.
Biosensors 12 00189 g023
Figure 24. The surface plot of a correlation coefficient matrix split in 12 × 12 sections corresponding to the monitored signal pairs, with a 10 × 10 section providing the quantification of interdependency on shifted versions of the corresponding signals.
Figure 24. The surface plot of a correlation coefficient matrix split in 12 × 12 sections corresponding to the monitored signal pairs, with a 10 × 10 section providing the quantification of interdependency on shifted versions of the corresponding signals.
Biosensors 12 00189 g024
Figure 25. The correlation coefficient matrix surface plots, generated from the signals acquired with the proposed physiograph during the walking trials, for (a) healthy control, (b) Patient 5, (c) Patient 1, and (d) Patient 3.
Figure 25. The correlation coefficient matrix surface plots, generated from the signals acquired with the proposed physiograph during the walking trials, for (a) healthy control, (b) Patient 5, (c) Patient 1, and (d) Patient 3.
Biosensors 12 00189 g025
Table 1. CNN training parameters.
Table 1. CNN training parameters.
CNN ParametersMobileNetEfficientNetB0Xception
Image dimension160 × 160 × 3160 × 160 × 3160 × 160 × 3
Learning rate0.050.050.05
Epochs150150150
Batch size3264128
OptimizerAdamAdamAdam
Loss functionBinary Cross EntropyBinary Cross EntropyBinary Cross Entropy
Output layer activation functionSoftmaxSoftmaxSoftmax
Table 2. Biomechanical parameters of the test group.
Table 2. Biomechanical parameters of the test group.
Test GroupArm Balance TremorHeel StrikeLift-Off
Healthy control 1PresentAbsentLeft: presentLeft: present
Right: presentRight: present
Healthy control 2PresentAbsentLeft: absent (flat-foot strike)Left: present
Right: absent (flat-foot strike)Right: present
Healthy control 3PresentAbsentLeft: presentLeft: present
Right: presentRight: present
Healthy control 4PresentAbsentLeft: presentLeft: present
Right: presentRight: present
Healthy control 5PresentAbsentLeft: presentLeft: present
Right: presentRight: present
Patient 1PresentLarge tremorLeft: absent (flat-foot strike)Left: present
Right: absent (flat-foot strike)Right: present
Patient 2AbsentLarge tremorLeft: absent (flat-foot strike)Left: present
Right: presentRight: present
Patient 3AbsentLarge tremorLeft: absent (flat-foot strike)Left: present
Right: absentRight: absent
Patient 4PresentSmall tremorLeft: presentLeft: present
Right: presentRight: present
Patient 5PresentSmall tremorLeft: absent (flat-foot strike)Left: present
Right: presentRight: present
Table 3. Temporal parameters of the test group.
Table 3. Temporal parameters of the test group.
Cadence (Steps/min)Single Support (%)Double Support (%)S/DStride Time Variability (%)
Healthy Control 162Left: 41.814.9Left: 2.83.8
Right: 43.3Right: 2.9
Healthy Control 227Left: 36.229.8Left: 1.29.8
Right: 34Right: 1.1
Healthy Control 324Left: 35.923.4Left: 1.54.8
Right: 40.8Right: 1.7
Healthy Control 456Left: 38.517.1Left: 2.23.1
Right: 44.3Right: 2.6
Healthy Control 555Left: 3919.7Left: 25.7
Right: 41Right: 2.1
Statistics for healthy control group44.8 ± 17.9Left: 38.3 ± 2.421 ± 5.9Left: 1.9 ± 0.65.4 ± 2.6
Right: 40.7 ± 4Right: 2.1 ± 0.7
Patient 139Left: 38.425.6Left: 1.56.3
Right: 36Right: 1.4
Patient 244Left: 29.131.6Left: 0.93.4
Right: 39.2Right: 1.2
Patient 348Left: 4720.6Left: 2.39
Right: 32.3Right: 1.6
Patient 432Left: 48.58.9Left: 5.413.1
Right: 42.5Right: 4.8
Patient 524Left: 3537.5Left: 0.98.2
Right: 27.4Right: 0.7
Statistics for PD group37.4 ± 9.6Left: 39.6 ± 8.124.9 ± 10.9Left: 2.2 ± 1.98 ± 3.6
Right: 35.5 ± 5.9Right: 1.9 ± 1.6
Table 4. Performance metrics of the neural network for gait pattern discrimination.
Table 4. Performance metrics of the neural network for gait pattern discrimination.
Performance MetricsCNN Model
MobilNetEfficientNetB0Xception
Accuracy0.950.850.85
Error0.190.390.25
Sensitivity0.900.850.9
Specificity0.960.800.73
Precision0.950.850.74
Table 5. Classification performances of the proposed work in comparison to other solutions reported in the literature.
Table 5. Classification performances of the proposed work in comparison to other solutions reported in the literature.
Performance MetricsClassifier
This Work[32][33][34][65]
ClassifierCNN
MobilNet
Random ForestCNNCNNSVM
Accuracy0.950.960.630.8560.917
Error0.19----
Sensitivity0.900.982-0.880.752
Specificity0.960.96-0.840.99
Precision0.950.972-0.86-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ileșan, R.R.; Cordoș, C.-G.; Mihăilă, L.-I.; Fleșar, R.; Popescu, A.-S.; Perju-Dumbravă, L.; Faragó, P. Proof of Concept in Artificial-Intelligence-Based Wearable Gait Monitoring for Parkinson’s Disease Management Optimization. Biosensors 2022, 12, 189. https://doi.org/10.3390/bios12040189

AMA Style

Ileșan RR, Cordoș C-G, Mihăilă L-I, Fleșar R, Popescu A-S, Perju-Dumbravă L, Faragó P. Proof of Concept in Artificial-Intelligence-Based Wearable Gait Monitoring for Parkinson’s Disease Management Optimization. Biosensors. 2022; 12(4):189. https://doi.org/10.3390/bios12040189

Chicago/Turabian Style

Ileșan, Robert Radu, Claudia-Georgiana Cordoș, Laura-Ioana Mihăilă, Radu Fleșar, Ana-Sorina Popescu, Lăcrămioara Perju-Dumbravă, and Paul Faragó. 2022. "Proof of Concept in Artificial-Intelligence-Based Wearable Gait Monitoring for Parkinson’s Disease Management Optimization" Biosensors 12, no. 4: 189. https://doi.org/10.3390/bios12040189

APA Style

Ileșan, R. R., Cordoș, C. -G., Mihăilă, L. -I., Fleșar, R., Popescu, A. -S., Perju-Dumbravă, L., & Faragó, P. (2022). Proof of Concept in Artificial-Intelligence-Based Wearable Gait Monitoring for Parkinson’s Disease Management Optimization. Biosensors, 12(4), 189. https://doi.org/10.3390/bios12040189

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