Instantaneous 3D EEG Signal Analysis Based on Empirical Mode Decomposition and the Hilbert–Huang Transform Applied to Depth of Anaesthesia

Depth of anaesthesia (DoA) is an important measure for assessing the degree to which the central nervous system of a patient is depressed by a general anaesthetic agent, depending on the potency and concentration with which anaesthesia is administered during surgery. We can monitor the DoA by observing the patient’s electroencephalography (EEG) signals during the surgical procedure. Typically high frequency EEG signals indicates the patient is conscious, while low frequency signals mean the patient is in a general anaesthetic state. If the anaesthetist is able to observe the instantaneous frequency changes of the patient’s EEG signals during surgery this can help to better regulate and monitor DoA, reducing surgical and post-operative risks. This paper describes an approach towards the development of a 3D real-time visualization application which can show the instantaneous frequency and instantaneous amplitude of EEG simultaneously by using empirical mode decomposition (EMD) and the Hilbert–Huang transform (HHT). HHT uses the EMD method to decompose a signal into so-called intrinsic mode functions (IMFs). The Hilbert spectral analysis method is then used to obtain instantaneous frequency data. The HHT provides a new method of analyzing non-stationary and nonlinear time series data. We investigate this approach by analyzing EEG data collected from patients undergoing surgical procedures. The results show that the EEG differences between three distinct surgical stages computed by using sample entropy (SampEn) are consistent with the expected differences between these stages based on the bispectral index (BIS), which has been shown to be quantifiable measure of the effect of anaesthetics on the central nervous system. Also, the proposed filtering approach is more effective compared to the standard filtering method in filtering out signal noise resulting in more consistent results than those provided by the BIS. The proposed approach is therefore able to distinguish between key operational stages related to DoA, which is consistent with the clinical observations. SampEn can also be viewed as a useful index for evaluating and monitoring the DoA of a patient when used in combination with this approach

Heterogeneous data fusion for brain psychology applications

This thesis aims to apply Empirical Mode Decomposition (EMD), Multiscale Entropy (MSE), and collaborative adaptive filters for the monitoring of different brain consciousness states. Both block based and online approaches are investigated, and a possible extension to the monitoring and identification of Electromyograph (EMG) states is provided. Firstly, EMD is employed as a multiscale time-frequency data driven tool to decompose a signal into a number of band-limited oscillatory components; its data driven nature makes EMD an ideal candidate for the analysis of nonlinear and non-stationary data. This methodology is further extended to process multichannel real world data, by making use of recent theoretical advances in complex and multivariate EMD. It is shown that this can be used to robustly measure higher order features in multichannel recordings to robustly indicate ‘QBD’. In the next stage, analysis is performed in an information theory setting on multiple scales in time, using MSE. This enables an insight into the complexity of real world recordings. The results of the MSE analysis and the corresponding statistical analysis show a clear difference in MSE between the patients in different brain consciousness states. Finally, an online method for the assessment of the underlying signal nature is studied. This method is based on a collaborative adaptive filtering approach, and is shown to be able to approximately quantify the degree of signal nonlinearity, sparsity, and non-circularity relative to the constituent subfilters. To further illustrate the usefulness of the proposed data driven multiscale signal processing methodology, the final case study considers a human-robot interface based on a multichannel EMG analysis. A preliminary analysis shows that the same methodology as that applied to the analysis of brain cognitive states gives robust and accurate results. The analysis, simulations, and the scope of applications presented suggest great potential of the proposed multiscale data processing framework for feature extraction in multichannel data analysis. Directions for future work include further development of real-time feature map approaches and their use across brain-computer and brain-machine interface applications

Energy weighting method and its application to fault diagnosis of rolling bearing

Feature extraction of vibration signal is the key factor of machine fault diagnosis. This paper proposes a novel method of capturing shock energy based on multi-scale weight evaluation of high definition time-frequency map. Specifically, the proposed method is conducted by the following steps. First, ensemble empirical mode decomposition (EEMD) preprocessor-based Hilbert-Huang Transform (HHT) is applied to the raw signal for high definition time-frequency map acquisition. Second, an original algorithm named multi-scale binary spectrum was applied to extract impulsive energy features with different frequency characteristics. Then weights of energy can be calculated by dimensionality reduction of each binary spectrum and merged by summation after blank processing. Finally, power spectrum of compound weight of energy can reveal characteristic frequency corresponding to specific fault. In this method, the key point is enhancement of frequency spectrum using higher dimension details. The process of multi-scale binarization and weight summation were aligned and the effectiveness is verified by simulated signal processing and actual case of train bearing experiment

BREAKDOWN OF CAUSALITY AND CORTICAL DOWNSTATE WITHIN THE SLEEPING BRAIN

Theoretically, consciousness depends on the brain\u2019s ability to engage in complex activity patterns that are, at once, distributed among interacting cortical areas (integrated) and differentiated in space and time (information-rich). In a recent series of experiments the electroencephalographic response to a direct cortical stimulation in humans was recorded during wakefulness and non-rapid eyes movement sleep (NREM) by means of a combination of transcranial magnetic stimulation (TMS) and high-density electroencephalogram (hd-EEG). TMS/hd-EEG measurements showed that, while during wakefulness the brain is able to sustain long-range specific patterns of activation, during NREM sleep, when consciousness fades, this ability is lost: the thalamocortical system, despite being active and reactive, either breaks down in causally independent modules (producing a local slow wave), or it bursts into an explosive and non-specific response (producing a global EEG slow wave). We hypothesize that, like spontaneous sleep slow waves, the slow waves triggered by TMS during deep sleep are due to bistability between periods of hyperpolarized down-state in cortical neurons, and periods of activation (up-state). In this condition, the inescapable occurrence of a silent, down-state after an initial activation could impair the ability of thalamocortical circuits to sustain long-range, differentiated patterns of activation, a theoretical requisite for consciousness. According to animal experiments the extracellular signature of the downstate is a transient suppression of high frequency (20Hz) oscillations \u2013 that is followed by a loss of both PLF and PLV, in spite of restored levels of neuronal activity. These results point to bistability as the underlying critical mechanism that prevents the emergence of complex interactions in human thalamocortical networks when consciousness is lost during NREM sleep. This finding is particularly relevant because a similar mechanism may play a role in other conditions where loss of consciousness is paralleled by the appearance of spontaneous (or TMS evoked) slow waves such as some kind of anesthesia and in brain injured subjects

Integrated Real-Time Control And Processing Systems For Multi-Channel Near-Infrared Spectroscopy Based Brain Computer Interfaces

This thesis outlines approaches to improve the signal processing and anal- ysis of Near-infrared spectroscopy (NIRS) based brain-computer interfaces (BCI). These approaches were developed in conjunction with the implemen- tation of a new customized exible multi-channel NIRS based BCI hardware system (Soraghan, 2010). Using a comparable functional imaging modality the assumptions on which NIRS-BCI have been reassessed, with regard to cognitive task selection, active area locations and lateralized motor cortex activation separability. This dissertation will also present methods that have been implemented to allow reduced hardware requirements in future NIRS-BCI development. We will also examine the sources of homeostatic physiological interference and present new approaches for analysis and at- tenuation within a real-time NIRS-BCI paradigm

Neonatal seizure detection based on single-channel EEG: instrumentation and algorithms

Seizure activity in the perinatal period, which constitutes the most common neurological emergency in the neonate, can cause brain disorders later in life or even death depending on their severity. This issue remains unsolved to date, despite the several attempts in tackling it using numerous methods. Therefore, a method is still needed that can enable neonatal cerebral activity monitoring to identify those at risk. Currently, electroencephalography (EEG) and amplitude-integrated EEG (aEEG) have been exploited for the identification of seizures in neonates, however both lack automation. EEG and aEEG are mainly visually analysed, requiring a specific skill set and as a result the presence of an expert on a 24/7 basis, which is not feasible. Additionally, EEG devices employed in neonatal intensive care units (NICU) are mainly designed around adults, meaning that their design specifications are not neonate specific, including their size due to multi-channel requirement in adults - adults minimum requirement is ≥ 32 channels, while gold standard in neonatal is equal to 10; they are bulky and occupy significant space in NICU. This thesis addresses the challenge of reliably, efficiently and effectively detecting seizures in the neonatal brain in a fully automated manner. Two novel instruments and two novel neonatal seizure detection algorithms (SDAs) are presented. The first instrument, named PANACEA, is a high-performance, wireless, wearable and portable multi-instrument, able to record neonatal EEG, as well as a plethora of (bio)signals. This device despite its high-performance characteristics and ability to record EEG, is mostly suggested to be used for the concurrent monitoring of other vital biosignals, such as electrocardiogram (ECG) and respiration, which provide vital information about a neonate's medical condition. The two aforementioned biosignals constitute two of the most important artefacts in the EEG and their concurrent acquisition benefit the SDA by providing information to an artefact removal algorithm. The second instrument, called neoEEG Board, is an ultra-low noise, wireless, portable and high precision neonatal EEG recording instrument. It is able to detect and record minute signals (< 10 nVp) enabling cerebral activity monitoring even from lower layers in the cortex. The neoEEG Board accommodates 8 inputs each one equipped with a patent-pending tunable filter topology, which allows passband formation based on the application. Both the PANACEA and the neoEEG Board are able to host low- to middle-complexity SDAs and they can operate continuously for at least 8 hours on 3-AA batteries. Along with PANACEA and the neoEEG Board, two novel neonatal SDAs have been developed. The first one, termed G prime-smoothed (G ́_s), is an on-line, automated, patient-specific, single-feature and single-channel EEG based SDA. The G ́_s SDA, is enabled by the invention of a novel feature, termed G prime (G ́) and can be characterised as an energy operator. The trace that the G ́_s creates, can also be used as a visualisation tool because of its distinct change at a presence of a seizure. Finally, the second SDA is machine learning (ML)-based and uses numerous features and a support vector machine (SVM) classifier. It can be characterised as automated, on-line and patient-independent, and similarly to G ́_s it makes use of a single-channel EEG. The proposed neonatal SDA introduces the use of the Hilbert-Huang transforms (HHT) in the field of neonatal seizure detection. The HHT analyses the non-linear and non-stationary EEG signal providing information for the signal as it evolves. Through the use of HHT novel features, such as the per intrinsic mode function (IMF) (0-3 Hz) sub-band power, were also employed. Detection rates of this novel neonatal SDA is comparable to multi-channel SDAs.Open Acces

Connectivity Analysis in EEG Data: A Tutorial Review of the State of the Art and Emerging Trends

Understanding how different areas of the human brain communicate with each other is a crucial issue in neuroscience. The concepts of structural, functional and effective connectivity have been widely exploited to describe the human connectome, consisting of brain networks, their structural connections and functional interactions. Despite high-spatial-resolution imaging techniques such as functional magnetic resonance imaging (fMRI) being widely used to map this complex network of multiple interactions, electroencephalographic (EEG) recordings claim high temporal resolution and are thus perfectly suitable to describe either spatially distributed and temporally dynamic patterns of neural activation and connectivity. In this work, we provide a technical account and a categorization of the most-used data-driven approaches to assess brain-functional connectivity, intended as the study of the statistical dependencies between the recorded EEG signals. Different pairwise and multivariate, as well as directed and non-directed connectivity metrics are discussed with a pros-cons approach, in the time, frequency, and information-theoretic domains. The establishment of conceptual and mathematical relationships between metrics from these three frameworks, and the discussion of novel methodological approaches, will allow the reader to go deep into the problem of inferring functional connectivity in complex networks. Furthermore, emerging trends for the description of extended forms of connectivity (e.g., high-order interactions) are also discussed, along with graph-theory tools exploring the topological properties of the network of connections provided by the proposed metrics. Applications to EEG data are reviewed. In addition, the importance of source localization, and the impacts of signal acquisition and pre-processing techniques (e.g., filtering, source localization, and artifact rejection) on the connectivity estimates are recognized and discussed. By going through this review, the reader could delve deeply into the entire process of EEG pre-processing and analysis for the study of brain functional connectivity and learning, thereby exploiting novel methodologies and approaches to the problem of inferring connectivity within complex networks

Design of a Simulator for Neonatal Multichannel EEG: Application to Time-Frequency Approaches for Automatic Artifact Removal and Seizure Detection

The electroencephalogram (EEG) is used to noninvasively monitor brain activities; it is the most utilized tool to detect abnormalities such as seizures. In recent studies, detection of neonatal EEG seizures has been automated to assist neurophysiologists in diagnosing EEG as manual detection is time consuming and subjective; however it still lacks the necessary robustness that is required for clinical implementation. Moreover, as EEG is intended to record the cerebral activities, extra-cerebral activities external to the brain are also recorded; these are called “artifacts” and can seriously degrade the accuracy of seizure detection. Seizures are one of the most common neurologic problems managed by hospitals occurring in 0.1%-0.5% livebirths. Neonates with seizures are at higher risk for mortality and are reported to be 55-70 times more likely to have severe cerebral-palsy. Therefore, early and accurate detection of neonatal seizures is important to prevent long-term neurological damage. Several attempts in modelling the neonatal EEG and artifacts have been done, but most did not consider the multichannel case. Furthermore, these models were used to test artifact or seizure detection separately, but not together. This study aims to design synthetic models that generate clean or corrupted multichannel EEG to test the accuracy of available artifact and seizure detection algorithms in a controlled environment. In this thesis, synthetic neonatal EEG model is constructed by using; single-channel EEG simulators, head model, 21-electrodes, and propagation equations, to produce clean multichannel EEG. Furthermore, neonatal EEG artifact model is designed using synthetic signals to corrupt EEG waveforms. After that, an automated EEG artifact detection and removal system is designed in both time and time-frequency domains. Artifact detection is optimised and removal performance is evaluated. Finally, an automated seizure detection technique is developed, utilising fused and extended multichannel features along a cross-validated SVM classifier. Results show that the synthetic EEG model mimics real neonatal EEG with 0.62 average correlation, and corrupted-EEG can degrade seizure detection average accuracy from 100% to 70.9%. They also show that using artifact detection and removal enhances the average accuracy to 89.6%, and utilising the extended features enhances it to 97.4% and strengthened its robustness.لمراقبة ورصد أنشطة واشارات المخ، دون الحاجة لأي عملیات (EEG) یستخدم الرسم أو التخطیط الكھربائي للدماغ للدماغجراحیة، وھي تعد الأداة الأكثر استخداما في الكشف عن أي شذوذأو نوبات غیر طبیعیة مثل نوبات الصرع. وقد أظھرت دراسات حدیثة، أن الكشف الآلي لنوبات حدیثي الولادة، ساعد علماء الفسیولوجیا العصبیة في تشخیص الاشارات الدماغیة بشكل أكبر من الكشف الیدوي، حیث أن الكشف الیدوي یحتاج إلى وقت وجھد أكبر وھوذو فعالیة أقل بكثیر، إلا أنھ لا یزال یفتقر إلى المتانة الضروریة والمطلوبة للتطبیق السریري.علاوة على ذلك؛ فكما یقوم الرسم الكھربائي بتسجیل الأنشطة والإشارات الدماغیة الداخلیة، فھو یسجل أیضا أي نشاط أو اشارات خارجیة، مما یؤدي إلى -(artifacts) :حدوث خلل في مدى دقة وفعالیة الكشف عن النوبات الدماغیة الداخلیة، ویطلق على تلك الاشارات مسمى (نتاج صنعي) . 0.5٪ولادة حدیثة في -٪تعد نوبات الصرع من أكثر المشكلات العصبیة انتشارا،ً وھي تصیب ما یقارب 0.1المستشفیات. حیث أن حدیثي الولادة المصابین بنوبات الصرع ھم أكثر عرضة للوفاة، وكما تشیر التقاریر الى أنھم 70مرة أكثر. لذا یعد الكشف المبكر والدقیق للنوبات الدماغیة -معرضین للإصابة بالشلل الدماغي الشدید بما یقارب 55لحدیثي الولادة مھم جدا لمنع الضرر العصبي على المدى الطویل. لقد تم القیام بالعدید من المحاولات التي كانتتھدف الى تصمیم نموذج التخطیط الكھربائي والنتاج الصنعي لدماغ حدیثي الولادة, إلا أن معظمھا لم یعر أي اھتمام الى قضیة تعدد القنوات. إضافة الى ذلك, استخدمت ھذه النماذج , كل على حدة, أو نوبات الصرع. تھدف ھذه الدراسة الى تصمیم نماذج مصطنعة من شأنھا (artifact) لإختبار كاشفات النتاج الصنعيأن تولد اشارات دماغیة متعددة القنوات سلیمة أو معطلة وذلك لفحص مدى دقة فعالیة خوارزمیات الكشف عن نوبات ضمن بیئة یمكن السیطرة علیھا. (artifact) الصرع و النتاج الصنعي في ھذه الأطروحة, یتكون نموذج الرسم الكھربائي المصطنع لحدیثي الولادة من : قناة محاكاة واحده للرسم الكھربائي, نموذج رأس, 21قطب كھربائي و معادلات إنتشار. حیث تھدف جمیعھا لإنتاج إشاراة سلیمة متعدده القنوات للتخطیط عن طریق استخدام اشارات مصطنعة (artifact) الكھربائي للدماغ.علاوة على ذلك, لقد تم تصمیم نموذجالنتاج الصنعيفي نطاقالوقت و (artifact) لإتلاف الرسم الكھربائي للدماغ. بعد ذلك تم انشاء برنامج لكشف و إزالةالنتاج الصناعينطاقالوقت و التردد المشترك. تم تحسین برنامج الكشف النتاج الصناعيالى ابعد ما یمكن بینما تمت عملیة تقییم أداء الإزالة. وفي الختام تم التمكن من تطویر تقنیة الكشف الآلي عن نوبات الصرع, وذلك بتوظیف صفات مدمجة و صفات الذي تم التأكد من صحتھ. (SVM) جدیدة للقنوات المتعددة لإستخدامھا للمصنفلقد أظھرت النتائج أن نموذج الرسم الكھربائي المصطنع لحدیثي الولادة یحاكي الرسمالكھربائي الحقیقي لحدیثي الولادة بمتوسط ترابط 0.62, و أنالرسم الكھربائي المتضرر للدماغ قد یؤدي الى حدوث ھبوطفي مدى دقة متوسط الكشف عن نوبات الصرع من 100%الى 70.9%. وقد أشارت أیضا الى أن استخدام الكشف والإزالة عن النتاج الصنعي (artifact) یؤدي الى تحسن مستوى الدقة الى نسبة 89.6 %, وأن توظیف الصفات الجدیدة للقنوات المتعددة یزید من تحسنھا لتصل الى نسبة 94.4 % مما یعمل على دعم متانتھا

Analysis of pattern dynamics for a nonlinear model of the human cortex via bifurcation theories

This thesis examines the bifurcations, i.e., the emergent behaviours, for the Waikato cortical model under the influence of the gap-junction inhibitory diffusion D₂ (identified as the Turing bifurcation parameter) and the time-to-peak for hyperpolarising GABA response γi (i.e., inhibitory rate-constant, identified as the Hopf bifurcation parameter). The cortical model simplifies the entire cortex to a cylindrical macrocolumn (∼ 1 mm³) containing ∼ 10⁵ neurons (85% excitatory, 15% inhibitory) communicating via both chemical and electrical (gap-junction) synapses. The linear stability analysis of the model equations predict the emergence of a Turing instability (in which separated areas of the cortex become activated) when gap-junction diffusivity is increased above a critical level. In addition, a Hopf bifurcation (oscillation) occurs when the inhibitory rate-constant is sufficiently small. Nonlinear interaction between these instabilities leads to spontaneous cortical patterns of neuronal activities evolving in space and time. Such model dynamics of delicately balanced interplay between Turing and Hopf instabilities may be of direct relevance to clinically observed brain dynamics such as epileptic seizure EEG spikes, deep-sleep slow-wave oscillations and cognitive gamma-waves. The relationship between the modelled brain patterns and model equations can normally be inferred from the eigenvalue dispersion curve, i.e., linear stability analysis. Sometimes we experienced mismatches between the linear stability analysis and the formed cortical patterns, which hampers us in identifying the type of instability corresponding to the emergent patterns. In this thesis, I investigate the pattern-forming mechanism of the Waikato cortical model to better understand the model nonlinearities. I first study the pattern dynamics via analysis of a simple pattern-forming system, the Brusselator model, which has a similar model structure and bifurcation phenomena as the cortical model. I apply both linear and nonlinear perturbation methods to analyse the near-bifurcation behaviour of the Brusselator in order to precisely capture the dominant mode that contributes the most to the final formed-patterns. My nonlinear analysis of the Brusselator model yields Ginzburg-Landau type amplitude equations that describe the dynamics of the most unstable mode, i.e., the dominant mode, in the vicinity of a bifurcation point. The amplitude equations at a Turing point unfold three characteristic spatial structures: honeycomb Hπ, stripes, and reentrant honeycomb H₀. A codimension-2 Turing–Hopf point (CTHP) predicts three mixed instabilities: stable Turing–Hopf (TH), chaotic TH, and bistable TH. The amplitude equations precisely determine the bifurcation conditions for these instabilities and explain the pattern-competition mechanism once the bifurcation parameters cross the thresholds, whilst driving the system into a nonlinear region where the linear stability analysis may not be applicable. Then, I apply the bifurcation theories to the cortical model for its pattern predictions. Analogous to the Brusselator model, I find cortical Turing pattens in Hπ, stripes and H₀ spatial structures. Moreover, I develop the amplitude equations for the cortical model, with which I derive the envelope frequency for the beating-waves of a stable TH mode; and propose ideas regarding emergence of the cortical chaotic mode. Apart from these pattern dynamics that the cortical model shares with the Brusselator system, the cortical model also exhibits “eye-blinking” TH patterns latticed in hexagons with localised oscillations. Although we have not found biological significance of these model pattens, the developed bifurcation theories and investigated pattern-forming mechanism may enrich our modelling strategies and help us to further improve model performance. In the last chapter of this thesis, I introduce a Turing–Hopf mechanism for the anaesthetic slow-waves, and predict a coherence drop of such slow-waves with the induction of propofol anaesthesia. To test this hypothesis, I developed an EEG coherence analysing algorithm, EEG coherence, to automatically examine the clinical EEG recordings across multiple subjects. The result shows significantly decreased coherence along the fronto-occipital axis, and increased coherence along the left- and right-temporal axis. As the Waikato cortical model is spatially homogenous, i.e., there are no explicit front-to-back or right-to-left directions, it is unable to produce different coherence changes for different regions. It appears that the Waikato cortical model best represents the cortical dynamics in the frontal region. The theory of pattern dynamics suggests that a mode transition from wave–Turing–wave to Turing–wave–Turing introduces pattern coherence changes in both positive and negative directions. Thus, a further modelling improvement may be the introduction of a cortical bistable mode where Turing and wave coexist