Extraction of event-related signals from multichannel bioelectrical measurements

Localization accuracy of a common beamformer for the comparison of two conditions

The linearly constrained minimum variance beamformer is frequently used to reconstruct sources underpinning neuromagnetic recordings. When reconstructions must be compared across conditions, it is considered good practice to use a single, "common" beamformer estimated from all the data at once. This is to ensure that differences between conditions are not ascribable to differences in beamformer weights. Here, we investigate the localization accuracy of such a common beamformer. Based on theoretical derivations, we first show that the common beamformer leads to localization errors in source reconstruction. We then turn to simulations in which we attempt to reconstruct a (genuine) source in a first condition, while considering a second condition in which there is an (interfering) source elsewhere in the brain. We estimate maps of mislocalization and assess statistically the difference between "standard" and "common" beamformers. We complement our findings with an application to experimental MEG data. The results show that the common beamformer may yield significant mislocalization. Specifically, the common beamformer may force the genuine source to be reconstructed closer to the interfering source than it really is. As the same applies to the reconstruction of the interfering source, both sources are pulled closer together than they are. This observation was further illustrated in experimental data. Thus, although the common beamformer allows for the comparison of conditions, in some circumstances it introduces localization inaccuracies. We recommend alternative approaches to the general problem of comparing conditions.

SCTICA: Sub-packet constrained temporal ICA method for fMRI data analysis

Independent component analysis (ICA) has become a widely used method for functional magnetic resonance imaging (fMRI) data analysis. However, spatial ICA usually performs better than temporal ICA with regard to the stability and accuracy of functional connectivity detection, and temporal ICA is often not feasible when it is applied to the analysis of real fMRI data of the whole brain because of the excessive spatial dimensions. In this paper, to overcome these problems, we propose a sub-packet constrained temporal ICA (SCTICA) method to take advantage of the a priori information using a multi-objective optimization framework with the Newton iterative algorithm. Moreover, a splitting strategy is presented to improve the feasibility of the temporal ICA for whole brain fMRI data analysis. The experimental results of real data show that the splitting strategy improved the ability of the temporal ICA to analyze whole brain fMRI data. Furthermore, the experimental results also demonstrated that the proposed SCTICA method can not only improve the stability of the temporal ICA, but can also improve the functional connectivity detection ability compared with the classical ICA and ICA with a priori information methods. In brief, the proposed SCTICA method overcomes the problem that prevents temporal ICA from being applied to fMRI data of the whole brain, and the functional connectivity detection performance is greatly improved compared with that of traditional methods.

An improved multi-objective optimization-based CICA method with data-driver temporal reference for group fMRI data analysis

Group independent component analysis (GICA) has been successfully applied to study multi-subject functional magnetic resonance imaging (fMRI) data, and the group independent component (GIC) represents the commonality of all subjects in the group. However, some studies show that the performance of GICA can be improved by incorporating a priori information, which is not always considered when looking for GICs in existing GICA methods. In this paper, we propose an improved multi-objective optimization-based constrained independent component analysis (CICA) method to take advantage of the temporal a priori information extracted from all subjects in the group by incorporating it into the computational process of GICA for group fMRI data analysis. The experimental results of simulated and real data show that the activated regions and the time course detected by the improved CICA method are more accurate in some sense. Moreover, the GIC computed by the improved CICA method has a higher correlation with the corresponding independent component of each subject in the group, which means that the improved CICA method with the temporal a priori information extracted from the group can better reflect the commonality of the subjects. These results demonstrate that the improved CICA method has its own advantages in fMRI data analysis.

A novel fMRI group data analysis method based on data-driven reference extracting from group subjects

Group-independent component analysis (GICA) is a well-established blind source separation technique that has been widely applied to study multi-subject functional magnetic resonance imaging (fMRI) data. The group-independent components (GICs) represent the commonness of all of the subjects in the group. Similar to independent component analysis on the single-subject level, the performance of GICA can be improved for multi-subject fMRI data analysis by incorporating a priori information; however, a priori information is not always considered while looking for GICs in existing GICA methods, especially when no obvious or specific knowledge about an unknown group is available. In this paper, we present a novel method to extract the group intrinsic reference from all of the subjects of the group and then incorporate it into the GICA extraction procedure. Comparison experiments between FastICA and GICA with intrinsic reference (GICA-IR) are implemented on the group level with regard to the simulated, hybrid and real fMRI data. The experimental results show that the GICs computed by GICA-IR have a higher correlation with the corresponding independent component of each subject in the group, and the accuracy of activation regions detected by GICA-IR was also improved. These results have demonstrated the advantages of the GICA-IR method, which can better reflect the commonness of the subjects in the group.

MEG source imaging method using fast L1 minimum-norm and its applications to signals with brain noise and human resting-state source amplitude images

The present study developed a fast MEG source imaging technique based on Fast Vector-based Spatio-Temporal Analysis using a L1-minimum-norm (Fast-VESTAL) and then used the method to obtain the source amplitude images of resting-state magnetoencephalography (MEG) signals for different frequency bands. The Fast-VESTAL technique consists of two steps. First, L1-minimum-norm MEG source images were obtained for the dominant spatial modes of sensor-waveform covariance matrix. Next, accurate source time-courses with millisecond temporal resolution were obtained using an inverse operator constructed from the spatial source images of Step 1. Using simulations, Fast-VESTAL's performance was assessed for its 1) ability to localize multiple correlated sources; 2) ability to faithfully recover source time-courses; 3) robustness to different SNR conditions including SNR with negative dB levels; 4) capability to handle correlated brain noise; and 5) statistical maps of MEG source images. An objective pre-whitening method was also developed and integrated with Fast-VESTAL to remove correlated brain noise. Fast-VESTAL's performance was then examined in the analysis of human median-nerve MEG responses. The results demonstrated that this method easily distinguished sources in the entire somatosensory network. Next, Fast-VESTAL was applied to obtain the first whole-head MEG source-amplitude images from resting-state signals in 41 healthy control subjects, for all standard frequency bands. Comparisons between resting-state MEG sources images and known neurophysiology were provided. Additionally, in simulations and cases with MEG human responses, the results obtained from using conventional beamformer technique were compared with those from Fast-VESTAL, which highlighted the beamformer's problems of signal leaking and distorted source time-courses.

Independent complexity patterns in single neuron activity induced by static magnetic field

We applied a combination of fractal analysis and Independent Component Analysis (ICA) method to detect the sources of fractal complexity in snail Br neuron activity induced by static magnetic field of 2.7 mT. The fractal complexity of Br neuron activity was analyzed before (Control), during (MF), and after (AMF) exposure to the static magnetic field in six experimental animals. We estimated the fractal dimension (FD) of electrophysiological signals using Higuchi's algorithm, and empirical FD distributions. By using the Principal Component Analysis (PCA) and FastICA algorithm we determined the number of components, and defined the statistically independent components (ICs) in the fractal complexity of signal waveforms. We have isolated two independent components of the empirical FD distributions for each of three groups of data by using FastICA algorithm. ICs represent the sources of fractal waveforms complexity of Br neuron activity in particular experimental conditions. Our main results have shown that there could be two opposite intrinsic mechanisms in single snail Br neuron response to static magnetic field stimulation. We named identified ICs that correspond to those mechanisms - the component of plasticity and the component of elasticity. We have shown that combination of fractal analysis with ICA method could be very useful for the decomposition and identification of the sources of fractal complexity of bursting neuronal activity waveforms.

ICA decomposition of EEG signal for fMRI processing in epilepsy

In this study, we introduce a new approach to process simultaneous Electroencephalography and functional Magnetic Resonance Imaging (EEG-fMRI) data in epilepsy. The method is based on the decomposition of the EEG signal using independent component analysis (ICA) and the usage of the relevant components' time courses to define the event related model necessary to find the regions exhibiting fMRI signal changes related to interictal activity. This approach achieves a natural data-driven differentiation of the role of distinct types of interictal activity with different amplitudes and durations in the epileptogenic process. Agreement between the conventional method and this new approach was obtained in 6 out of 9 patients that had interictal activity inside the scanner. In all cases, the maximum Z-score was greater in the fMRI studies based on ICA component method and the extent of activation was increased in 5 out of the 6 cases in which overlap was found. Furthermore, the three cases where an agreement was not found were those in which no significant activation was found at all using the conventional approach.

Comparison of saccade-related EEG signal with saccade-related independent component

We have been research saccade-related EEG signals in order to predict beginning of saccade by EEG signal. We have already detected saccade-related EEG signal by ensemble averaging and saccade-related independent components (ICs) by independent component analysis (ICA). However, features of saccade-related EEG signals and saccade-related ICs were not compared. In this paper, saccade-related EEG signals and saccade-related ICs were compared in the point of the latency between starting time of a saccade and time when a saccade-related EEG signal or an IC has maximum value and in the point of the peak scale where a saccade-related EEG signal or an IC has maximum value.

Mapping scalp topographies of rhythmic EEG activity using temporal decorrelation based constrained ICA

Independent component analysis (ICA) methods are being increasingly applied to the analysis of electromagnetic (EM) brain signals. However, these powerful techniques still generally require subjective a posteriori analysis in order to visualise neurophysiologically meaningful components in the outputs. Standard implementations of ICA are restrictive mainly due to the square mixing assumption (i.e., as many sources as measurement channels) - this is especially so with large multichannel recordings. There are many instances in neurophysiological analysis where there is strong a priori information about the signals being sought; as in tracking the changing scalp topographies of rhythmic activities. Through constraining the ICA solution it is possible to extract signals that are statistically independent, yet which are similar to some reference signal which incorporates the a priori information. We demonstrate this method on a multichannel recording of an epileptiform electroencephalogram (EEG), where we automate the repeated simultaneous extraction of both rhythmic seizure activity, as well as alpha-band activity, over an epoch of EEG. Subjective analysis of the results shows scalp topographies with realistic spatial distributions which conform to our neurophysiologic expectations. This work shows that constraining ICA can be a very useful technique, especially in automated systems and we demonstrate that this can be successfully applied to EM brain signal analysis.

Independent component analysis of high-resolution imaging data identifies distinct functional domains

In the vertebrate brain external stimuli are often represented in distinct functional domains distributed across the cortical surface. Fast imaging techniques used to measure patterns of population activity record movies with many pixels and many frames, i.e., data sets with high dimensionality. Here we demonstrate that principal component analysis (PCA) followed by spatial independent component analysis (sICA), can be exploited to reduce the dimensionality of data sets recorded in the olfactory bulb and the somatosensory cortex of mice as well as the visual cortex of monkeys, without loosing the stimulus-specific responses. Different neuronal populations are separated based on their stimulus-specific spatiotemporal activation. Both, spatial and temporal response characteristics can be objectively obtained, simultaneously. In the olfactory bulb, groups of glomeruli with different response latencies can be identified. This is shown for recordings of olfactory receptor neuron input measured with a calcium-sensitive axon tracer and for network dynamics measured with the voltage-sensitive dye RH 1838. In the somatosensory cortex, barrels responding to the stimulation of single whiskers can be automatically detected. In the visual cortex orientation columns can be extracted. In all cases artifacts due to movement, heartbeat or respiration were separated from the functional signal by sICA and could be removed from the data set. sICA following PCA is therefore a powerful technique for data compression, unbiased analysis and dissection of imaging data of population activity, collected with high spatial and temporal resolution.

A Novel Extraction Method of Fetal Electrocardiogram From the Composite Abdominal Signal

In contrast to the ultrasonic measurement of fetal heart motion, the fetal electrocardiogram (ECG) provides clinically significant information concerning the electrophysiological state of a fetus. In this paper, a novel method for extracting the fetal ECG from abdominal composite signals is proposed. This method consists of the cancellation of the mother's ECG and blind source separation with the reference signal (BSSR). The cancellation of the mother's ECG component was performed by subtracting the linear combination of mutually orthogonal projections of the heart vector. The BSSR is a fixed-point algorithm, the Lagrange function of which includes the higher order cross-correlation between the extracted signal and the reference signal as the cost term rather than a constraint. This realizes the convexity of the Lagrange function in a simple form, which guarantees the convergence of the algorithm. By practical application, the proposed method has been shown to be able to extract the P and T waves in addition to the R wave. The reliability and accuracy of the proposed method was confirmed by comparing the extracted signals with the directly recorded ECG at the second stage of labor. The gestational age-dependency of the physiological parameters of the extracted fetal ECG also coincided well with that of the magnetocardiogram, which proves the clinical applicability of the proposed method.

Vector-based spatial–temporal minimum L1-norm solution for MEG

Minimum L1-norm solutions have been used by many investigators to analyze MEG responses because they provide high spatial resolution images. However, conventional minimum L1-norm approaches suffer from instability in spatial construction, and poor smoothness of the reconstructed source time-courses. Activity commonly "jumps" from one grid point to (usually) the neighboring grid points. Equivalently, the time-course of one specific grid point can show substantial "spiky-looking" discontinuity. In the present study, we present a new vector-based spatial-temporal analysis using a L1-minimum-norm (VESTAL). This approach is based on a principle of MEG physics: the magnetic waveforms in sensor-space are linear functions of the source time-courses in the imaging-space. Our computer simulations showed that VESTAL provides good reconstruction of the source amplitude and orientation, with high stability and resolution in both the spatial and temporal domains. "Spiky-looking" discontinuity was not observed in the source time-courses. Importantly, the simulations also showed that VESTAL can resolve sources that are 100% correlated. We then examined the performance of VESTAL in the analysis of human median-nerve MEG responses. The results demonstrated that this method easily distinguishes sources very spatially close to each other, including individual primary somatosensory areas (BA 1, 2, 3b), primary motor area (BA 4), and other regions in the somatosensory system (e.g., BA 5, 7, SII, SMA, and temporal-parietal junction) with high temporal stability and resolution. VESTAL's potential for obtaining information on source extent was also examined.

BCI Competition 2003—Data Set IIb: Enhancing P300 Wave Detection Using ICA-Based Subspace Projections for BCI Applications

An algorithm based on independent component analysis (ICA) is introduced for P300 detection. After ICA decomposition, P300-related independent components are selected according to the a priori knowledge of P300 spatio-temporal pattern, and clear P300 peak is reconstructed by back projection of ICA. Applied to the dataset IIb of BCI Competition 2003, the algorithm achieved an accuracy of 100% in P300 detection within five repetitions.

Application of independent component analysis to magnetic resonance imaging for enhancing the contrast of gray and white matter

An application of independent component analysis (ICA) was attempted to develop a method of processing magnetic resonance (MR) images to extract physiologically independent components representing tissue relaxation times and achieve improved visualization of normal and pathologic structures. Anatomical T1-weighted, T2-weighted and proton density images were obtained from 10 normal subjects, 3 patients with brain tumors and 1 patient with multiple sclerosis. The data sets were analyzed using ICA based on the learning rule of Bell and Sejnowski after prewhitening operations. The three independent components obtained from the three original data sets corresponded to (1) short T1 components representing myelin of white matter and lipids, (2) relatively short T1 components representing gray matter and (3) long T2 components representing free water. The involvement of gray or white matter in brain tumor cases and the demyelination in the case of multiple sclerosis were enhanced and visualized in independent component images. ICA can potentially achieve separation of tissues with different relaxation characteristics and generate new contrast images of gray and white matter. With the proper choice of contrast for the original images, ICA may be useful not only for extracting subtle or hidden changes but also for preprocessing transformation before clustering and segmenting the structure of the human brain.

Temporally constrained ICA: an application to artifact rejection in electromagnetic brain signal analysis

Independent component analysis (ICA) is a technique which extracts statistically independent components from a set of measured signals. The technique enjoys numerous applications in biomedical signal analysis in the literature, especially in the analysis of electromagnetic (EM) brain signals. Standard implementations of ICA are restrictive mainly due to the square mixing assumption-for signal recordings which have large numbers of channels, the large number of resulting extracted sources makes the subsequent analysis laborious and highly subjective. There are many instances in neurophysiological analysis where there is strong a priori information about the signals being sought; temporally constrained ICA (cICA) can extract signals that are statistically independent, yet which are constrained to be similar to some reference signal which can incorporate such a priori information. We demonstrate this method on a synthetic dataset and on a number of artifactual waveforms identified in multichannel recordings of EEG and MEG. cICA repeatedly converges to the desired component within a few iterations and subjective analysis shows the waveforms to be of the expected morphologies and with realistic spatial distributions. This paper shows that cICA can be applied with great success to EM brain signal analysis, with an initial application in automating artifact extraction in EEG and MEG.

Extraction of specific signals with temporal structure

In this work we develop a very simple batch learning algorithm for semiblind extraction of a desired source signal with temporal structure from linear mixtures. Although we use the concept of sequential blind extraction of sources and independent component analysis, we do not carry out the extraction in a completely blind manner; neither do we assume that sources are statistically independent. In fact, we show that the a priori information about the autocorrelation function of primary sources can be used to extract the desired signals (sources of interest) from their linear mixtures. Extensive computer simulations and real data application experiments confirm the validity and high performance of the proposed algorithm.

Imaging Brain Dynamics Using Independent Component Analysis

The analysis of electroencephalographic (EEG) and magnetoencephalographic (MEG) recordings is important both for basic brain research and for medical diagnosis and treatment. Independent component analysis (ICA) is an effective method for removing artifacts and separating sources of the brain signals from these recordings. A similar approach is proving useful for analyzing functional magnetic resonance brain imaging (fMRI) data. In this paper, we outline the assumptions underlying ICA and demonstrate its application to a variety of electrical and hemodynamic recordings from the human brain.

Independence: a new criterion for the analysis of the electromagnetic fields in the global brain?

The impressive increase in the understanding of some basic processing in the human brain has recently led to the formulation of efficient computational methods, which when applied in the design of better signal processing tools, provides a deeper and clearer view to study the functioning of the human brain. The recently developed independent component analysis (ICA) has been shown to be an efficient tool for artifact identification and extraction from electroencephalographic and magnetoencephalographic recordings. In addition, ICA has been applied to the analysis of brain signals evoked by sensory stimuli. Extensions of the basic ICA methodology have also been employed to reveal otherwise hidden information. This paper reviews our recent results in this field.

Independent component approach to the analysis of EEG and MEG recordings

Multichannel recordings of the electromagnetic fields emerging from neural currents in the brain generate large amounts of data. Suitable feature extraction methods are, therefore, useful to facilitate the representation and interpretation of the data. Recently developed independent component analysis (ICA) has been shown to be an efficient tool for artifact identification and extraction from electroencephalographic (EEG) and magnetoencephalographic (MEG) recordings. In addition, ICA has been applied to the analysis of brain signals evoked by sensory stimuli. This paper reviews our recent results in this field.