Bogdan Mijović

Source Separation From Single-Channel Recordings by Combining Empirical-Mode Decomposition and Independent Component Analysis

In biomedical signal processing, it is often the case that many sources are mixed into the measured signal. The goal is usually to analyze one or several of them separately. In the case of multichannel measurements, several blind source separation techniques are available for decomposing the signal into its components [e.g., independent component analysis (ICA)]. However, only a few techniques have been reported for analyses of single-channel recordings. Examples are single-channel ICA (SCICA) and wavelet-ICA (WICA), which all have certain limitations. In this paper, we propose a new method for a single-channel signal decomposition. This method combines empirical-mode decomposition with ICA. We compare the separation performance of our algorithm with SCICA and WICA through simulations, and we show that our method outperforms the other two, especially for high noise-to-signal ratios. The performance of the new algorithm was also demonstrated in two real-life applications.

The "why" and "how" of JointICA: results from a visual detection task

Since several years, neuroscience research started to focus on multimodal approaches. One such multimodal approach is the combination of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI). However, no standard integration procedure has been established so far. One promising data-driven approach consists of a joint decomposition of event-related potentials (ERPs) and fMRI maps derived from the response to a particular stimulus. Such an algorithm (joint independent component analysis or JointICA) has recently been proposed by Calhoun et al. (2006). This method provides sources with both a fine spatial and temporal resolution, and has shown to provide meaningful results. However, the algorithms performance has not been fully characterized yet, and no procedure has been proposed to assess the quality of the decomposition. In this paper, we therefore try to answer why and how JointICA works. We show the performance of the algorithm on data obtained in a visual detection task, and compare the performance for EEG recorded simultaneously with fMRI data and for EEG recorded in a separate session (outside the scanner room). We perform several analyses in order to set the necessary conditions that lead to a sound decomposition, and to give additional insights for exploration in future studies. In that respect, we show how the algorithm behaves when different EEG electrodes are used and we test the robustness with respect to the number of subjects in the study. The performance of the algorithm in all the experiments is validated based on results from previous studies.

The dynamics of contour integration: A simultaneous EEG-fMRI study

To study the dynamics of contour integration in the human brain, we simultaneously acquired EEG and fMRI data while participants were engaged in a passive viewing task. The stimuli were Gabor arrays with some Gabor elements positioned on the contour of an embedded shape, in three conditions: with local and global structure (perfect contour alignment), with global structure only (orthogonal orientations interrupting the alignment), or without contour. By applying JointICA to the EEG and fMRI responses of the subjects, new insights could be obtained that cannot be derived from unimodal recordings. In particular, only in the global structure condition, an ERP peak around 300ms was identified that involved a loop from LOC to the early visual areas. This component can be interpreted as being related to the verification of the consistency of the different local elements with the globally defined shape, which is necessary when perfect local-to-global alignment is absent. By modifying JointICA, a quantitative comparison of brain regions and the time-course of their interplay were obtained between different conditions. More generally, we provide additional support for the presence of feedback loops from higher areas to lower level sensory regions.

ICA extracts epileptic sources from fMRI in EEG-negative patients: a retrospective validation study.

Simultaneous EEG-fMRI has proven to be useful in localizing interictal epileptic activity. However, the applicability of traditional GLM-based analysis is limited as interictal spikes are often not seen on the EEG inside the scanner. Therefore, we aim at extracting epileptic activity purely from the fMRI time series using independent component analysis (ICA). To our knowledge, we show for the first time that ICA can find sources related to epileptic activity in patients where no interictal spikes were recorded in the EEG. The epileptic components were identified retrospectively based on the known localization of the ictal onset zone (IOZ). We demonstrate that the selected components truly correspond to epileptic activity, as sources extracted from patients resemble significantly better the IOZ than sources found in healthy controls. Furthermore, we show that the epileptic components in patients with and without spikes recorded inside the scanner resemble the IOZ in the same degree. We conclude that ICA of fMRI has the potential to extend the applicability of EEG-fMRI for presurgical evaluation in epilepsy.

Single trial ERP reading based on parallel factor analysis

The extraction of task-related single trial ERP features has recently gained much interest, in particular in simultaneous EEG-fMRI applications. In this study, a specific decomposition known as parallel factor analysis (PARAFAC) was used, in order to retrieve the task-related activity from the raw signals. Using visual detection task data, acquired in normal circumstances and simultaneously with fMRI, differences between distinct task-related conditions can be captured in the trial signatures of specific PARAFAC components when applied to ERP data arranged in Channels × Time × Trials arrays, but the signatures did not correlate with the fMRI data. Despite the need for parameter tuning and careful preprocessing, the approach is shown to be successful, especially when prior knowledge about the expected ERPs is incorporated.

A new and fast approach towards sEMG decomposition

The decomposition of high-density surface EMG (HD-sEMG) interference patterns into the contribution of motor units is still a challenging task. We introduce a new, fast solution to this problem. The method uses a data-driven approach for selecting a set of electrodes to enable discrimination of present motor unit action potentials (MUAPs). Then, using shapes detected on these channels, the hierarchical clustering algorithm as reported by Quian Quiroga et al. (Neural Comput 16:1661–1687, 2004) is extended for multichannel data in order to obtain the motor unit action potential (MUAP) signatures. After this first step, more motor unit firings are obtained using the extracted signatures by a novel demixing technique. In this demixing stage, we propose a time-efficient solution for the general convolutive system that models the motor unit firings on the HD-sEMG grid. We constrain this system by using the extracted signatures as prior knowledge and reconstruct the firing patterns in a computationally efficient way. The algorithm performance is successfully verified on simulated data containing up to 20 different MUAP signatures. Moreover, we tested the method on real low contraction recordings from the lateral vastus leg muscle by comparing the algorithm’s output to the results obtained by manual analysis of the data from two independent trained operators. The proposed method showed to perform about equally successful as the operators.

Synergistic control of forearm based on accelerometer data and artificial neural networks

In the present study, we modeled a reaching task as a two-link mechanism. The upper arm and forearm motion trajectories during vertical arm movements were estimated from the measured angular accelerations with dual-axis accelerometers. A data set of reaching synergies from able-bodied individuals was used to train a radial basis function artificial neural network with upper arm/forearm tangential angular accelerations. The trained radial basis function artificial neural network for the specific movements predicted forearm motion from new upper arm trajectories with high correlation (mean, 0.9149-0.941). For all other movements, prediction was low (range, 0.0316-0.8302). Results suggest that the proposed algorithm is successful in generalization over similar motions and subjects. Such networks may be used as a high-level controller that could predict forearm kinematics from voluntary movements of the upper arm. This methodology is suitable for restoring the upper limb functions of individuals with motor disabilities of the forearm, but not of the upper arm. The developed control paradigm is applicable to upper-limb orthotic systems employing functional electrical stimulation. The proposed approach is of great significance particularly for humans with spinal cord injuries in a free-living environment. The implication of a measurement system with dual-axis accelerometers, developed for this study, is further seen in the evaluation of movement during the course of rehabilitation. For this purpose, training-related changes in synergies apparent from movement kinematics during rehabilitation would characterize the extent and the course of recovery. As such, a simple system using this methodology is of particular importance for stroke patients. The results underlie the important issue of upper-limb coordination.

Bayesian model selection of template forward models for EEG source reconstruction.

Several EEG source reconstruction techniques have been proposed to identify the generating neuronal sources of electrical activity measured on the scalp. The solution of these techniques depends directly on the accuracy of the forward model that is inverted. Recently, a parametric empirical Bayesian (PEB) framework for distributed source reconstruction in EEG/MEG was introduced and implemented in the Statistical Parametric Mapping (SPM) software. The framework allows us to compare different forward modeling approaches, using real data, instead of using more traditional simulated data from an assumed true forward model. In the absence of a subject specific MR image, a 3-layered boundary element method (BEM) template head model is currently used including a scalp, skull and brain compartment. In this study, we introduced volumetric template head models based on the finite difference method (FDM). We constructed a FDM head model equivalent to the BEM model and an extended FDM model including CSF. These models were compared within the context of three different types of source priors related to the type of inversion used in the PEB framework: independent and identically distributed (IID) sources, equivalent to classical minimum norm approaches, coherence (COH) priors similar to methods such as LORETA, and multiple sparse priors (MSP). The resulting models were compared based on ERP data of 20 subjects using Bayesian model selection for group studies. The reconstructed activity was also compared with the findings of previous studies using functional magnetic resonance imaging. We found very strong evidence in favor of the extended FDM head model with CSF and assuming MSP. These results suggest that the use of realistic volumetric forward models can improve PEB EEG source reconstruction.

Decoupling between fundamental frequency and energy envelope of neonate cries

BACKGROUND The presence of decoupling, i.e. the absence of coupling between fundamental frequency variation and energy envelope during phonetic crying, and its extent, reflects the degree of maturation of the central nervous system. AIM We hereby wanted to assess the existence and extent of decoupling in term neonates (neurodevelopmental relevance) and whether an association between decoupling and clinical pain expression could be unveiled (clinical relevance). STUDY DESIGN To assess decoupling in healthy term neonates during procedural pain, newborns were videotaped and crying was recorded during venous blood sampling. Besides acoustic analysis, pain expression was quantified based on the Modified Behavioral Pain Scale (MBPS). SUBJECTS 47 healthy term neonates underwent venous blood puncture at the 3rd day of life. OUTCOME MEASURES Beside the MBPS score, the correlation coefficients were calculated between the fundamental frequency variation and energy envelope of the cries. RESULTS Based on data collected in 47 healthy term neonates, correlation coefficients varied between 0.20 and 0.68. The degree of decoupling displayed extensive variability between the neonates and also in different cry bouts in a crying sequence within an individual neonate. A negative association was found between MBPS value and decoupling (r(2)=-0.12), the same as for the intra-subject variability although less extensive (r(2)=-0.02). CONCLUSION Decoupling only relates weakly with the amount of distress in 3day old newborns, even though a great intra-subject variability is present. This study suggests that there is no evidence of extensive decoupling as the newborn still has to fully develop the control of larynx and abdominal muscles.

Combining EMD with ICA for extracting independent sources from single channel and two-channel data

Blind Source Separation (BSS) techniques are frequently needed in the processing of biomedical signals. This need comes from the fact that these signals are often composed of many different sources, which are mixed in the measured signal. However, we are usually only interested in examining one or a limited set of sources of interest separately. A variety of algorithms exist for separating multichannel mixtures into its independent sources (e.g. different Independent Component Analysis (ICA) techniques). These techniques only work if the number of channels is larger than, or equal to the number of sources present in the signal. On the other hand, only a few algorithms have been reported for the analysis of single channel sources, or other mixtures where the number of sources is higher than the number of channels. In this work we show a new technique which combines Empirical Mode Decomposition (EMD) and Independent Component Analysis (ICA). We will show that this technique is capable in separating independent sources when the number of these sources is higher than the number of channels available. We show the performance in single channel and two-channel biosignal processing.