Attenuation of artifacts in EEG signals measured inside an MRI scanner using constrained independent component analysis

Artifact Reduction in Simultaneous EEG-fMRI: A Systematic Review of Methods and Contemporary Usage

Simultaneous electroencephalography-functional MRI (EEG-fMRI) is a technique that combines temporal (largely from EEG) and spatial (largely from fMRI) indicators of brain dynamics. It is useful for understanding neuronal activity during many different event types, including spontaneous epileptic discharges, the activity of sleep stages, and activity evoked by external stimuli and decision-making tasks. However, EEG recorded during fMRI is subject to imaging, pulse, environment and motion artifact, causing noise many times greater than the neuronal signals of interest. Therefore, artifact removal methods are essential to ensure that artifacts are accurately removed, and EEG of interest is retained. This paper presents a systematic review of methods for artifact reduction in simultaneous EEG-fMRI from literature published since 1998, and an additional systematic review of EEG-fMRI studies published since 2016. The aim of the first review is to distill the literature into clear guidelines for use of simultaneous EEG-fMRI artifact reduction methods, and the aim of the second review is to determine the prevalence of artifact reduction method use in contemporary studies. We find that there are many published artifact reduction techniques available, including hardware, model based, and data-driven methods, but there are few studies published that adequately compare these methods. In contrast, recent EEG-fMRI studies show overwhelming use of just one or two artifact reduction methods based on literature published 15–20 years ago, with newer methods rarely gaining use outside the group that developed them. Surprisingly, almost 15% of EEG-fMRI studies published since 2016 fail to adequately describe the methods of artifact reduction utilized. We recommend minimum standards for reporting artifact reduction techniques in simultaneous EEG-fMRI studies and suggest that more needs to be done to make new artifact reduction techniques more accessible for the researchers and clinicians using simultaneous EEG-fMRI.

Removal of BCG artefact from concurrent fMRI-EEG recordings based on EMD and PCA

BACKGROUND

Simultaneous electroencephalography (EEG) and functional magnetic resonance image (fMRI) acquisitions provide better insight into brain dynamics. Some artefacts due to simultaneous acquisition pose a threat to the quality of the data. One such problematic artefact is the ballistocardiogram (BCG) artefact.

METHODS

We developed a hybrid algorithm that combines features of empirical mode decomposition (EMD) with principal component analysis (PCA) to reduce the BCG artefact. The algorithm does not require extra electrocardiogram (ECG) or electrooculogram (EOG) recordings to extract the BCG artefact.

RESULTS

The method was tested with both simulated and real EEG data of 11 participants. From the simulated data, the similarity index between the extracted BCG and the simulated BCG showed the effectiveness of the proposed method in BCG removal. On the other hand, real data were recorded with two conditions, i.e. resting state (eyes closed dataset) and task influenced (event-related potentials (ERPs) dataset). Using qualitative (visual inspection) and quantitative (similarity index, improved normalized power spectrum (INPS) ratio, power spectrum, sample entropy (SE)) evaluation parameters, the assessment results showed that the proposed method can efficiently reduce the BCG artefact while preserving the neuronal signals.

COMPARISON WITH EXISTING METHODS

Compared with conventional methods, namely, average artefact subtraction (AAS), optimal basis set (OBS) and combined independent component analysis and principal component analysis (ICA-PCA), the statistical analyses of the results showed that the proposed method has better performance, and the differences were significant for all quantitative parameters except for the power and sample entropy.

CONCLUSIONS

The proposed method does not require any reference signal, prior information or assumption to extract the BCG artefact. It will be very useful in circumstances where the reference signal is not available.

Methods for pulse artefact reduction: experiences with EEG data recorded at 9.4 T static magnetic field

BACKGROUND

The feasibility of recording electroencephalography (EEG) at ultra-high static magnetic fields up to 9.4 T was recently demonstrated and is expected to be incorporated into functional magnetic resonance imaging (fMRI) studies at 9.4 T. Correction of the pulse artefact (PA) is a significant challenge since its amplitude is proportional to the strength of the magnetic field in which EEG is recorded.

NEW METHOD

We conducted a study in which different PA correction methods were applied to EEG data recorded inside a 9.4 T scanner in order to retrieve visual P100 and auditory P300 evoked potentials. We explored different PA reduction methods, including the optimal basis set (OBS) method as well as objective and subjective component rejection using independent component analysis (ICA).

RESULTS

ICA followed by objective rejection of components is optimal for retrieving visual P100 and auditory P300 from EEG data recorded inside the scanner.

COMPARISON WITH EXISTING METHODS

Previous studies suggest that OBS or OBS followed by ICA are optimal for retrieving evoked potentials at 3T. In our EEG data recorded at 9.4 T OBS performed alone was not fully optimal for the identification of evoked potentials. OBS followed by ICA was partially effective.

CONCLUSIONS

In this study ICA has been shown to be an important tool for correcting the PA in EEG data recorded at 9.4 T, particularly when automated rejection of components is performed.

Crosslinking EEG time-frequency decomposition and fMRI in error monitoring

Recent studies implicate a common response monitoring system, being active during erroneous and correct responses. Converging evidence from time-frequency decompositions of the response-related ERP revealed that evoked theta activity at fronto-central electrode positions differentiates correct from erroneous responses in simple tasks, but also in more complex tasks. However, up to now it is unclear how different electrophysiological parameters of error processing, especially at the level of neural oscillations are related, or predictive for BOLD signal changes reflecting error processing at a functional-neuroanatomical level. The present study aims to provide crosslinks between time domain information, time-frequency information, MRI BOLD signal and behavioral parameters in a task examining error monitoring due to mistakes in a mental rotation task. The results show that BOLD signal changes reflecting error processing on a functional-neuroanatomical level are best predicted by evoked oscillations in the theta frequency band. Although the fMRI results in this study account for an involvement of the anterior cingulate cortex, middle frontal gyrus, and the Insula in error processing, the correlation of evoked oscillations and BOLD signal was restricted to a coupling of evoked theta and anterior cingulate cortex BOLD activity. The current results indicate that although there is a distributed functional-neuroanatomical network mediating error processing, only distinct parts of this network seem to modulate electrophysiological properties of error monitoring.

Statistical feature extraction for artifact removal from concurrent fMRI-EEG recordings

We propose a set of algorithms for sequentially removing artifacts related to MRI gradient switching and cardiac pulsations from electroencephalography (EEG) data recorded during functional magnetic resonance imaging (fMRI). Special emphasis is directed upon the use of statistical metrics and methods for the extraction and selection of features that characterize gradient and pulse artifacts. To remove gradient artifacts, we use channel-wise filtering based on singular value decomposition (SVD). To remove pulse artifacts, we first decompose data into temporally independent components and then select a compact cluster of components that possess sustained high mutual information with the electrocardiogram (ECG). After the removal of these components, the time courses of remaining components are filtered by SVD to remove the temporal patterns phase-locked to the cardiac timing markers derived from the ECG. The filtered component time courses are then inversely transformed into multi-channel EEG time series free of pulse artifacts. Evaluation based on a large set of simultaneous EEG-fMRI data obtained during a variety of behavioral tasks, sensory stimulations and resting conditions showed excellent data quality and robust performance attainable with the proposed methods. These algorithms have been implemented as a Matlab-based toolbox made freely available for public access and research use.

Extraction of P300 using constrained independent component analysis

A brain computer interface (BCI) uses electrophysiological activities of the brain such as natural rhythms and evoked potentials to communicate with some external devices. P300 is a positive evoked potential (EP), elicited approximately 300 ms after an attended external stimulus. A P300-based BCI uses this evoked potential as a means of communication with the external devices. Until now this P300-based BCI has been rather slow, as it is difficult to detect a P300 response without averaging over a number of trials. Previously, independent component analysis (ICA) has been used in the extraction of P300. However, the drawback of ICA is that it extracts not only P300 but also non-P300 related components requiring a proper selection of P300 ICs by the system. In this study we propose an algorithm based on constrained independent component analysis (cICA) for P300 extraction which can extract only the relevant component by incorporating a priori information. A reference signal is generated as this a priori information of P300 and cICA is applied to extract the P300 related component. Then the extracted P300 IC is segmented, averaged, and classified into target and non-target events by means of a linear classifier. The method is fast, reliable, computationally inexpensive as compared to ICA and achieves an accuracy of 98.3% in the detection of P300.