Search

The analysis of simultaneous EEG and fMRI data is generally based on the extraction of regressors of interest from the EEG, which are correlated to the fMRI data in a general linear model setting. In more advanced approaches, the spatial information of EEG is also exploited by assuming underlying dipole models. In this study, we present a semi automatic and efficient method to determine electrode positions from electrode gel artifacts, facilitating the integration of EEG and fMRI in future EEG/fMRI data models. In order to visualize all electrode artifacts simultaneously in a single view, a surface rendering of the structural MRI is made using a skin triangular mesh model as reference surface, which is expanded to a "pancake view". Then the electrodes are determined with a simple mouse click for each electrode. Using the geometry of the skin surface and its transformation to the pancake view, the 3D coordinates of the electrodes are reconstructed in the MRI coordinate frame. The electrode labels are attached to the electrode positions by fitting a template grid of the electrode cap in which the labels are known. The correspondence problem between template and sample electrodes is solved by minimizing a cost function over rotations, shifts and scalings of the template grid. The crucial step here is to use the solution of the so-called "Hungarian algorithm" as a cost function, which makes it possible to identify the electrode artifacts in arbitrary order. The template electrode grid has to be constructed only once for each cap configuration. In our implementation of this method, the whole procedure can be performed within 15 min including import of MRI, surface reconstruction and transformation, electrode identification and fitting to template. The method is robust in the sense that an electrode template created for one subject can be used without identification errors for another subject for whom the same EEG cap was used. Furthermore, the method appears to be robust against spurious or missing artifacts. We therefore consider the proposed method as a useful and reliable tool within the larger toolbox required for the analysis of co-registered EEG/fMRI data. © 2011 Elsevier Inc.

The simultaneous recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can be used to localize interictal epileptiform discharges (IEDs). Previous studies have reported varying degrees of concordance of EEG-fMRI with electroclinical findings. The aim of the present study is to evaluate to what extent this variability is determined by the analytical strategy or by the properties of the EEG data. For that purpose, 42 IED sets obtained in 29 patients with epilepsy were reanalyzed using a finite impulse response approach, which estimates the hemodynamic response function (HRF) from the data and allows non-causal effects. Cardiac effects were treated as additional confounders in the model. This approach was compared to the classical approach assuming a fixed HRF for each voxel in the brain. The performance of each method was assessed by comparing the fMRI results to the EEG focus. The flexible model revealed more significantly activated voxels, which resulted in more activated brain regions concordant with the EEG focus (26 vs. 16). Correction for cardiac effects improved the results in 7 out of the 42 data sets. Furthermore, design theory for event-related experiments was applied in order to determine the influence of the number of IEDs and their temporal distribution on the success of an experiment. It appeared that this success is highly dependent upon the number of IEDs present during the recording and less on their temporal spacing. We conclude that the outcome of EEG-fMRI can be improved by using an optimized analytical strategy, but also depends on the number of IEDs occurring during the recording.

Fitting an equivalent current dipole (ECD) to magnetoencephalography (MEG) data is a widely used method for analyzing MEG data, especially in clinical applications. However, as is well known, the iterative fitting algorithms routinely employed can converge on a local minima far from the best fit. While a host of fitting algorithms have been proposed to remedy this problem, the simple exhaustive search optimization algorithm has not been considered in the literature — possibly because it is assumed to be too slow for routine use. Taking advantage of the speed of modern computers, it is demonstrated that using exhaustive search to fit the parameters of a single ECD to 151 MEG channels yields both robust and reproducible fits within reasonable computation times. On a high end personal computer, with a 2 mm grid spacing over the head, exhaustive search can fit 450 ECD's per hour with only an additional 5 min spent on calculating the lead fields. As the algorithm is fully automatic, it is very useful when the goal is to detect ECD's in a large dataset. Also, an interpreter can be provided with parameters from the ECD that fits for all latencies after stimulus, allowing a more informed decision of which latency to use for further analysis.

In six healthy adults we examined the sources underlying P1 and N2 of the motion VEP. For this purpose we acquired, in addition to the VEP, MRI images and patterns of regional cerebral blood flow with SPECT for three of the subjects. With the same motion stimulus we also examined the spatial distribution of N2 in children. In both adults and children left and right half-field responses were compared. It was found that N2 is generated by extrastriate activity and that motion stimuli are not equivalently processed in the two cerebral hemispheres. In adults, N2 dominates in one hemisphere irrespective of the visual half-field used for stimulation whereas children show an ipsilateral maximum for N2 upon half-field presentation.