GPGPU-aided ensemble empirical-mode decomposition for EEG analysis during anesthesia

GPU Implementation of the Improved CEEMDAN Algorithm for Fast and Efficient EEG Time-Frequency Analysis

Time-frequency analysis of EEG data is a key step in exploring the internal activities of the human brain. Studying oscillations is an important part of the analysis, as they are thought to provide the underlying mechanism for communication between neural assemblies. Traditional methods of analysis, such as Short-Time FFT and Wavelet Transforms, are not ideal for this task due to the time-frequency uncertainty principle and their reliance on predefined basis functions. Empirical Mode Decomposition and its variants are more suited to this task as they are able to extract the instantaneous frequency and phase information but are too time consuming for practical use. Our aim was to design and develop a massively parallel and performance-optimized GPU implementation of the Improved Complete Ensemble EMD with the Adaptive Noise (CEEMDAN) algorithm that significantly reduces the computational time (from hours to seconds) of such analysis. The resulting GPU program, which is publicly available, was validated against a MATLAB reference implementation and reached over a 260× speedup for actual EEG measurement data, and provided predicted speedups in the range of 3000-8300× for longer measurements when sufficient memory was available. The significance of our research is that this implementation can enable researchers to perform EMD-based EEG analysis routinely, even for high-density EEG measurements. The program is suitable for execution on desktop, cloud, and supercomputer systems and can be the starting point for future large-scale multi-GPU implementations.

A Novel Method Based on Combination of Independent Component Analysis and Ensemble Empirical Mode Decomposition for Removing Electrooculogram Artifacts From Multichannel Electroencephalogram Signals

Electrooculogram (EOG) is one of common artifacts in recorded electroencephalogram (EEG) signals. Many existing methods including independent component analysis (ICA) and wavelet transform were applied to eliminate EOG artifacts but ignored the possible impact of the nature of EEG signal. Therefore, the removal of EOG artifacts still faces a major challenge in EEG research. In this paper, the ensemble empirical mode decomposition (EEMD) and ICA algorithms were combined to propose a novel EEMD-based ICA method (EICA) for removing EOG artifacts from multichannel EEG signals. First, the ICA method was used to decompose original EEG signals into multiple independent components (ICs), and the EOG-related ICs were automatically identified through the kurtosis method. Then, by performing the EEMD algorithm on EOG-related ICs, the intrinsic mode functions (IMFs) linked to EOG were discriminated and eliminated. Finally, artifact-free IMFs were projected to obtain the ICs without EOG artifacts, and the clean EEG signals were ultimately reconstructed by the inversion of ICA. Both EOGs correction from simulated EEG signals and real EEG data were studied, which verified that the proposed method could achieve an improved performance in EOG artifacts rejection. By comparing with other existing approaches, the EICA obtained the optimal performance with the highest increase in signal-to-noise ratio and decrease in root mean square error and correlation coefficient after EOG artifacts removal, which demonstrated that the proposed method could more effectively eliminate blink artifacts from multichannel EEG signals with less error influence. This study provided a novel promising method to eliminate EOG artifacts with high performance, which is of great importance for EEG signals processing and analysis.

Monitoring the level of hypnosis using a hierarchical SVM system

Monitoring level of hypnosis is a major ongoing challenge for anesthetists to reduce anesthetic drug consumption, avoiding intraoperative awareness and prolonged recovery. This paper proposes a novel automated method for accurate assessing of the level of hypnosis with sevoflurane in 17 patients using the electroencephalogram signal. In this method, a set of distinctive features and a hierarchical classification structure based on support vector machine (SVM) methods, is proposed to discriminate the four levels of anesthesia (awake, light, general and deep states). The first stage of the hierarchical SVM structure identifies the awake state by extracting Shannon Permutation Entropy, Detrended Fluctuation Analysis and frequency features. Then deep state is identified by extracting the sample entropy feature; and finally light and general states are identified by extracting the three mentioned features of the first step. The accuracy of the proposed method of analyzing the brain activity during anesthesia is 94.11%; which was better than previous studies and also a commercial monitoring system (Response Entropy Index).

A Novel GPU-Based Acceleration Algorithm for a Longwave Radiative Transfer Model

Graphics processing unit (GPU)-based computing for climate system models is a longstanding research area of interest. The rapid radiative transfer model for general circulation models (RRTMG), a popular atmospheric radiative transfer model, can calculate atmospheric radiative fluxes and heating rates. However, the RRTMG has a high calculation time, so it is urgent to study its GPU-based efficient acceleration algorithm to enable large-scale and long-term climatic simulations. To improve the calculative efficiency of radiation transfer, this paper proposes a GPU-based acceleration algorithm for the RRTMG longwave radiation scheme (RRTMG_LW). The algorithm concept is accelerating the RRTMG_LW in the g- p o i n t dimension. After implementing the algorithm in CUDA Fortran, the G-RRTMG_LW was developed. The experimental results indicated that the algorithm was effective. In the case without I/O transfer, the G-RRTMG_LW on one K40 GPU obtained a speedup of 30.98× over the baseline performance on one single Intel Xeon E5-2680 CPU core. When compared to its counterpart running on 10 CPU cores of an Intel Xeon E5-2680 v2, the G-RRTMG_LW on one K20 GPU in the case without I/O transfer achieved a speedup of 2.35×.

Complexity of Wake Electroencephalography Correlates With Slow Wave Activity After Sleep Onset

Sleep electroencephalography (EEG) provides an opportunity to study sleep scientifically, whose chaotic, dynamic, complex, and dissipative nature implies that non-linear approaches could uncover some mechanism of sleep. Based on well-established complexity theories, one hypothesis in sleep medicine is that lower complexity of brain waves at pre-sleep state can facilitate sleep initiation and further improve sleep quality. However, this has never been studied with solid data. In this study, EEG collected from healthy subjects was used to investigate the association between pre-sleep EEG complexity and sleep quality. Multiscale entropy analysis (MSE) was applied to pre-sleep EEG signals recorded immediately after light-off (while subjects were awake) for measuring the complexities of brain dynamics by a proposed index, CI₁₋₃₀. Slow wave activity (SWA) in sleep, which is commonly used as an indicator of sleep depth or sleep intensity, was quantified based on two methods, traditional Fast Fourier transform (FFT) and ensemble empirical mode decomposition (EEMD). The associations between wake EEG complexity, sleep latency, and SWA in sleep were evaluated. Our results demonstrated that lower complexity before sleep onset is associated with decreased sleep latency, indicating a potential facilitating role of reduced pre-sleep complexity in the wake-sleep transition. In addition, the proposed EEMD-based method revealed an association between wake complexity and quantified SWA in the beginning of sleep (90 min after sleep onset). Complexity metric could thus be considered as a potential indicator for sleep interventions, and further studies are encouraged to examine the application of EEG complexity before sleep onset in populations with difficulty in sleep initiation. Further studies may also examine the mechanisms of the causal relationships between pre-sleep brain complexity and SWA, or conduct comparisons between normal and pathological conditions.

Localization of Active Brain Sources From EEG Signals Using Empirical Mode Decomposition: A Comparative Study

The localization of active brain sources from Electroencephalogram (EEG) is a useful method in clinical applications, such as the study of localized epilepsy, evoked-related-potentials, and attention deficit/hyperactivity disorder. The distributed-source model is a common method to estimate neural activity in the brain. The location and amplitude of each active source are estimated by solving the inverse problem by regularization or using Bayesian methods with spatio-temporal constraints. Frequency and spatio-temporal constraints improve the quality of the reconstructed neural activity. However, separation into frequency bands is beneficial when the relevant information is in specific sub-bands. We improved frequency-band identification and preserved good temporal resolution using EEG pre-processing techniques with good frequency band separation and temporal resolution properties. The identified frequency bands were included as constraints in the solution of the inverse problem by decomposing the EEG signals into frequency bands through various methods that offer good frequency and temporal resolution, such as empirical mode decomposition (EMD) and wavelet transform (WT). We present a comparative analysis of the accuracy of brain-source reconstruction using these techniques. The accuracy of the spatial reconstruction was assessed using the Wasserstein metric for real and simulated signals. We approached the mode-mixing problem, inherent to EMD, by exploring three variants of EMD: masking EMD, Ensemble-EMD (EEMD), and multivariate EMD (MEMD). The results of the spatio-temporal brain source reconstruction using these techniques show that masking EMD and MEMD can largely mitigate the mode-mixing problem and achieve a good spatio-temporal reconstruction of the active sources. Masking EMD and EEMD achieved better reconstruction than standard EMD, Multiple Sparse Priors, or wavelet packet decomposition when EMD was used as a pre-processing tool for the spatial reconstruction (averaged over time) of the brain sources. The spatial resolution obtained using all three EMD variants was substantially better than the use of EMD alone, as the mode-mixing problem was mitigated, particularly with masking EMD and EEMD. These findings encourage further exploration into the use of EMD-based pre-processing, the mode-mixing problem, and its impact on the accuracy of brain source activity reconstruction.

Design and Evaluation of a Real Time Physiological Signals Acquisition System Implemented in Multi-Operating Rooms for Anesthesia

With critical importance of medical healthcare, there exist urgent needs for in-depth medical studies that can access and analyze specific physiological signals to provide theoretical support for practical clinical care. As a consequence, obtaining the valuable medical data with minimal cost and impacts on hospital work comes as the first concern of researchers. Anesthesia plays a widely recognized role in surgeries, which attracts people to undertake relevant research. In this paper, a real-time physiological medical signal data acquisition system (PMSDA) for the multi-operating room applications is proposed with high universality of the hospital practical settings and research requirements. By utilizing a wireless communication approach, it provides an easily accessible network platform for collection of physiological medical signals such as photoplethysmogram (PPG), electrocardiograph (ECG) and electroencephalogram (EEG) during the surgery. In addition, the raw data is stored on a server for safe backup and further analysis of depth of anesthesia (DoA). Results show that the PMSDA exhibits robust, high quality performance and efficiently reduces costs compared to previously manual methods and allows seamless integration into hospital environment, independent of its routine work. Overall, it provides a pragmatic and flexible surgery-data acquisition system model with low impact and resource cost applicable to research in critical and practical medical circumstances.

Monitoring the Depth of Anesthesia Using a New Adaptive Neurofuzzy System

Accurate and noninvasive monitoring of the depth of anesthesia (DoA) is highly desirable. Since the anesthetic drugs act mainly on the central nervous system, the analysis of brain activity using electroencephalogram (EEG) is very useful. This paper proposes a novel automated method for assessing the DoA using EEG. First, 11 features including spectral, fractal, and entropy are extracted from EEG signal and then, by applying an algorithm according to exhaustive search of all subsets of features, a combination of the best features (Beta-index, sample entropy, shannon permutation entropy, and detrended fluctuation analysis) is selected. Accordingly, we feed these extracted features to a new neurofuzzy classification algorithm, adaptive neurofuzzy inference system with linguistic hedges (ANFIS-LH). This structure can successfully model systems with nonlinear relationships between input and output, and also classify overlapped classes accurately. ANFIS-LH, which is based on modified classical fuzzy rules, reduces the effects of the insignificant features in input space, which causes overlapping and modifies the output layer structure. The presented method classifies EEG data into awake, light, general, and deep states during anesthesia with sevoflurane in 17 patients. Its accuracy is 92% compared to a commercial monitoring system (response entropy index) successfully. Moreover, this method reaches the classification accuracy of 93% to categorize EEG signal to awake and general anesthesia states by another database of propofol and volatile anesthesia in 50 patients. To sum up, this method is potentially applicable to a new real-time monitoring system to help the anesthesiologist with continuous assessment of DoA quickly and accurately.

Instantaneous phase difference analysis between thoracic and abdominal movement signals based on complementary ensemble empirical mode decomposition

BACKGROUND

Thoracoabdominal asynchrony is often adopted to discriminate respiratory diseases in clinics. Conventionally, Lissajous figure analysis is the most frequently used estimation of the phase difference in thoracoabdominal asynchrony. However, the temporal resolution of the produced results is low and the estimation error increases when the signals are not sinusoidal. Other previous studies have reported time-domain procedures with the use of band-pass filters for phase-angle estimation. Nevertheless, the band-pass filters need calibration for phase delay elimination.

METHODS

To improve the estimation, we propose a novel method (named as instantaneous phase difference) that is based on complementary ensemble empirical mode decomposition for estimating the instantaneous phase relation between measured thoracic wall movement and abdominal wall movement. To validate the proposed method, experiments on simulated time series and human-subject respiratory data with two breathing types (i.e., thoracic breathing and abdominal breathing) were conducted. Latest version of Lissajous figure analysis and automatic phase estimation procedure were compared.

RESULTS

The simulation results show that the standard deviations of the proposed method were lower than those of two other conventional methods. The proposed method performed more accurately than the two conventional methods. For the human-subject respiratory data, the results of the proposed method are in line with those in the literature, and the correlation analysis result reveals that they were positively correlated with the results generated by the two conventional methods. Furthermore, the standard deviation of the proposed method was also the smallest.

CONCLUSIONS

To summarize, this study proposes a novel method for estimating instantaneous phase differences. According to the findings from both the simulation and human-subject data, our approach was demonstrated to be effective. The method offers the following advantages: (1) improves the temporal resolution, (2) does not introduce a phase delay, (3) works with non-sinusoidal signals, (4) provides quantitative phase estimation without estimating the embedded frequency of breathing signals, and (5) works without calibrated measurements. The results demonstrate a higher temporal resolution of the phase difference estimation for the evaluation of thoracoabdominal asynchrony.

Theta Oscillation Reveals the Temporal Involvement of Different Attentional Networks in Contingent Reorienting

In the visual world, rapidly reorienting to relevant objects outside the focus of attention is vital for survival. This ability from the interaction between goal-directed and stimulus-driven attentional control is termed contingent reorienting. Neuroimaging studies have demonstrated activations of the ventral and dorsal attentional networks (DANs) which exhibit right hemisphere dominance, but the temporal dynamics of the attentional networks still remain unclear. The present study used event-related potential (ERP) to index the locus of spatial attention and Hilbert-Huang transform (HHT) to acquire the time-frequency information during contingent reorienting. The ERP results showed contingent reorienting induced significant N2pc on both hemispheres. In contrast, our time-frequency analysis found further that, unlike the N2pc, theta oscillation during contingent reorienting differed between hemispheres and experimental sessions. The inter-trial coherence (ITC) of the theta oscillation demonstrated that the two sides of the attentional networks became phase-locked to contingent reorienting at different stages. The left attentional networks were associated with contingent reorienting in the first experimental session whereas the bilateral attentional networks play a more important role in this process in the subsequent session. This phase-locked information suggests a dynamic temporal evolution of the involvement of different attentional networks in contingent reorienting and a potential role of the left ventral network in the first session.

EEG entropy measures in anesthesia

HIGHLIGHTS

► Twelve entropy indices were systematically compared in monitoring depth of anesthesia and detecting burst suppression.► Renyi permutation entropy performed best in tracking EEG changes associated with different anesthesia states.► Approximate Entropy and Sample Entropy performed best in detecting burst suppression.

OBJECTIVE

Entropy algorithms have been widely used in analyzing EEG signals during anesthesia. However, a systematic comparison of these entropy algorithms in assessing anesthesia drugs' effect is lacking. In this study, we compare the capability of 12 entropy indices for monitoring depth of anesthesia (DoA) and detecting the burst suppression pattern (BSP), in anesthesia induced by GABAergic agents.

METHODS

Twelve indices were investigated, namely Response Entropy (RE) and State entropy (SE), three wavelet entropy (WE) measures [Shannon WE (SWE), Tsallis WE (TWE), and Renyi WE (RWE)], Hilbert-Huang spectral entropy (HHSE), approximate entropy (ApEn), sample entropy (SampEn), Fuzzy entropy, and three permutation entropy (PE) measures [Shannon PE (SPE), Tsallis PE (TPE) and Renyi PE (RPE)]. Two EEG data sets from sevoflurane-induced and isoflurane-induced anesthesia respectively were selected to assess the capability of each entropy index in DoA monitoring and BSP detection. To validate the effectiveness of these entropy algorithms, pharmacokinetic/pharmacodynamic (PK/PD) modeling and prediction probability (Pk) analysis were applied. The multifractal detrended fluctuation analysis (MDFA) as a non-entropy measure was compared.

RESULTS

All the entropy and MDFA indices could track the changes in EEG pattern during different anesthesia states. Three PE measures outperformed the other entropy indices, with less baseline variability, higher coefficient of determination (R (2)) and prediction probability, and RPE performed best; ApEn and SampEn discriminated BSP best. Additionally, these entropy measures showed an advantage in computation efficiency compared with MDFA.

CONCLUSION

Each entropy index has its advantages and disadvantages in estimating DoA. Overall, it is suggested that the RPE index was a superior measure. Investigating the advantages and disadvantages of these entropy indices could help improve current clinical indices for monitoring DoA.

A novel dynamic update framework for epileptic seizure prediction

Epileptic seizure prediction is a difficult problem in clinical applications, and it has the potential to significantly improve the patients' daily lives whose seizures cannot be controlled by either drugs or surgery. However, most current studies of epileptic seizure prediction focus on high sensitivity and low false-positive rate only and lack the flexibility for a variety of epileptic seizures and patients' physical conditions. Therefore, a novel dynamic update framework for epileptic seizure prediction is proposed in this paper. In this framework, two basic sample pools are constructed and updated dynamically. Furthermore, the prediction model can be updated to be the most appropriate one for the prediction of seizures' arrival. Mahalanobis distance is introduced in this part to solve the problem of side information, measuring the distance between two data sets. In addition, a multichannel feature extraction method based on Hilbert-Huang transform and extreme learning machine is utilized to extract the features of a patient's preseizure state against the normal state. At last, a dynamic update epileptic seizure prediction system is built up. Simulations on Freiburg database show that the proposed system has a better performance than the one without update. The research of this paper is significantly helpful for clinical applications, especially for the exploitation of online portable devices.

EOG artifact correction from EEG recording using stationary subspace analysis and empirical mode decomposition

Ocular contamination of EEG data is an important and very common problem in the diagnosis of neurobiological events. An effective approach is proposed in this paper to remove ocular artifacts from the raw EEG recording. First, it conducts the blind source separation on the raw EEG recording by the stationary subspace analysis, which can concentrate artifacts in fewer components than the representative blind source separation methods. Next, to recover the neural information that has leaked into the artifactual components, the adaptive signal decomposition technique EMD is applied to denoise the components. Finally, the artifact-only components are projected back to be subtracted from EEG signals to get the clean EEG data. The experimental results on both the artificially contaminated EEG data and publicly available real EEG data have demonstrated the effectiveness of the proposed method, in particular for the cases where limited number of electrodes are used for the recording, as well as when the artifact contaminated signal is highly non-stationary and the underlying sources cannot be assumed to be independent or uncorrelated.

Measuring center of pressure signals to quantify human balance using multivariate multiscale entropy by designing a force platform

To assess the improvement of human body balance, a low cost and portable measuring device of center of pressure (COP), known as center of pressure and complexity monitoring system (CPCMS), has been developed for data logging and analysis. In order to prove that the system can estimate the different magnitude of different sways in comparison with the commercial Advanced Mechanical Technology Incorporation (AMTI) system, four sway tests have been developed (i.e., eyes open, eyes closed, eyes open with water pad, and eyes closed with water pad) to produce different sway displacements. Firstly, static and dynamic tests were conducted to investigate the feasibility of the system. Then, correlation tests of the CPCMS and AMTI systems have been compared with four sway tests. The results are within the acceptable range. Furthermore, multivariate empirical mode decomposition (MEMD) and enhanced multivariate multiscale entropy (MMSE) analysis methods have been used to analyze COP data reported by the CPCMS and compare it with the AMTI system. The improvements of the CPCMS are 35% to 70% (open eyes test) and 60% to 70% (eyes closed test) with and without water pad. The AMTI system has shown an improvement of 40% to 80% (open eyes test) and 65% to 75% (closed eyes test). The results indicate that the CPCMS system can achieve similar results to the commercial product so it can determine the balance.

Global Synchronization Measurement of Multivariate Neural Signals with Massively Parallel Nonlinear Interdependence Analysis

The estimation of synchronization amongst multiple brain regions is a critical issue in understanding brain functions. There is a lack of an appropriate approach which is capable of 1) measuring the direction and strength of synchronization of activities of multiple brain regions, and 2) adapting to the quickly increasing sizes and scales of neural signals. Nonlinear Interdependence (NLI) analysis is an effective method for measuring synchronization direction and strength of bivariate neural signal. However, the method currently does not directly apply in handling multivariate signal. Its application in practice has also long been largely hampered by the ultra-high complexity of NLI algorithms. Aiming at these problems, this study 1) extends the conventional NLI to quantify the global synchronization of multivariate neural signals, and 2) develops a parallelized NLI method with general-purpose computing on the graphics processing unit (GPGPU), namely, G-NLI. The approach performs synchronization measurement in a massively parallel manner. The G-NLI has improved the runtime performance by more than 1000 times comparing to the original sequential NLI. Meanwhile, the G-NLI was employed to analyze 10-channel local field potential (LFP) recordings from a patient suffering from temporal lobe epilepsy. The results demonstrate that the proposed G-NLI method can support real-time global synchronization measurement and it could be successful in localization of epileptic focus.

Consciousness and depth of anesthesia assessment based on Bayesian analysis of EEG signals

This study applies Bayesian techniques to analyze EEG signals for the assessment of the consciousness and depth of anesthesia (DoA). This method takes the limiting large-sample normal distribution as posterior inferences to implement the Bayesian paradigm. The maximum a posterior (MAP) is applied to denoise the wavelet coefficients based on a shrinkage function. When the anesthesia states change from awake to light, moderate, and deep anesthesia, the MAP values increase gradually. Based on these changes, a new function B(DoA) is designed to assess the DoA. The new proposed method is evaluated using anesthetized EEG recordings and BIS data from 25 patients. The Bland-Alman plot is used to verify the agreement of B(DoA) and the popular BIS index. A correlation between B(DoA) and BIS was measured using prediction probability P(K). In order to estimate the accuracy of DoA, the effect of sample n and variance τ on the maximum posterior probability is studied. The results show that the new index accurately estimates the patient's hypnotic states. Compared with the BIS index in some cases, the B(DoA) index can estimate the patient's hypnotic state in the case of poor signal quality.

Hilbert-Huang transformation-based time-frequency analysis methods in biomedical signal applications

Hilbert-Huang transformation, wavelet transformation, and Fourier transformation are the principal time-frequency analysis methods. These transformations can be used to discuss the frequency characteristics of linear and stationary signals, the time-frequency features of linear and non-stationary signals, the time-frequency features of non-linear and non-stationary signals, respectively. The Hilbert-Huang transformation is a combination of empirical mode decomposition and Hilbert spectral analysis. The empirical mode decomposition uses the characteristics of signals to adaptively decompose them to several intrinsic mode functions. Hilbert transforms are then used to transform the intrinsic mode functions into instantaneous frequencies, to obtain the signal's time-frequency-energy distributions and features. Hilbert-Huang transformation-based time-frequency analysis can be applied to natural physical signals such as earthquake waves, winds, ocean acoustic signals, mechanical diagnosis signals, and biomedical signals. In previous studies, we examined Hilbert-Huang transformation-based time-frequency analysis of the electroencephalogram FPI signals of clinical alcoholics, and 'sharp I' wave-based Hilbert-Huang transformation time-frequency features. In this paper, we discuss the application of Hilbert-Huang transformation-based time-frequency analysis to biomedical signals, such as electroencephalogram, electrocardiogram signals, electrogastrogram recordings, and speech signals.