Co-registration of EEG and MRI data using matching of spline interpolated and MRI-segmented reconstructions of the scalp surface

Variability of EEG electrode positions and their underlying brain regions: visualizing gel artifacts from a simultaneous EEG-fMRI dataset

INTRODUCTION

We investigated the between-subject variability of EEG (electroencephalography) electrode placement from a simultaneously recorded EEG-fMRI (functional magnetic resonance imaging) dataset.

METHODS

Neuro-navigation software was used to localize electrode positions, made possible by the gel artifacts present in the structural magnetic resonance images. To assess variation in the brain regions directly underneath electrodes we used MNI coordinates, their associated Brodmann areas, and labels from the Harvard-Oxford Cortical Atlas. We outline this relatively simple pipeline with accompanying analysis code.

RESULTS

In a sample of 20 participants, the mean standard deviation of electrode placement was 3.94 mm in x, 5.55 mm in y, and 7.17 mm in z, with the largest variation in parietal and occipital electrodes. In addition, the brain regions covered by electrode pairs were not always consistent; for example, the mean location of electrode PO7 was mapped to BA18 (secondary visual cortex), whereas PO8 was closer to BA19 (visual association cortex). Further, electrode C1 was mapped to BA4 (primary motor cortex), whereas C2 was closer to BA6 (premotor cortex).

CONCLUSIONS

Overall, the results emphasize the variation in electrode positioning that can be found even in a fixed cap. This may be particularly important to consider when using EEG positioning systems to inform non-invasive neurostimulation.

Semi-Automated and Direct Localization and Labeling of EEG Electrodes Using MR Structural Images for Simultaneous fMRI-EEG

Electroencephalography (EEG) source reconstruction estimates spatial information from the brain’s electrical activity acquired using EEG. This method requires accurate identification of the EEG electrodes in a three-dimensional (3D) space and involves spatial localization and labeling of EEG electrodes. Here, we propose a new approach to tackle this two-step problem based on the simultaneous acquisition of EEG and magnetic resonance imaging (MRI). For the step of spatial localization of electrodes, we extract the electrode coordinates from the curvature of the protrusions formed in the high-resolution T1-weighted brain scans. In the next step, we assign labels to each electrode based on the distinguishing feature of the electrode’s distance profile in relation to other electrodes. We then compare the subject’s electrode data with template-based models of prelabeled distance profiles of correctly labeled subjects. Based on this approach, we could localize EEG electrodes in 26 head models with over 90% accuracy in the 3D localization of electrodes. Next, we performed electrode labeling of the subjects’ data with progressive improvements in accuracy: with ∼58% accuracy based on a single EEG-template, with ∼71% accuracy based on 3 EEG-templates, and with ∼76% accuracy using 5 EEG-templates. The proposed semi-automated method provides a simple alternative for the rapid localization and labeling of electrodes without the requirement of any additional equipment than what is already used in an EEG-fMRI setup.

Impact of different registration methods in MEG source analysis

For this study the impact of different co-registration procedures on MEG source localization of somatosensory evoked fields was evaluated. Two different co-registration procedures were used to calculate the transformation matrix which specifies how to align the MRI data to the MEG head coordinate system. In order to depict the differences, caused by the method, the Euclidian distance between the reconstructed sources was noted. It was shown that, erroneous MRI and MEG data co- registration effects source localization results. Most dipoles are located more posterior and superior when the more advanced registration procedure was applied. In conclusion the results show, that an iterative matching procedure allows an accurate knowledge of the MEG gradiometer sensor position relative to the head which is crucial to correctly reconstruct neuronal activity derived from MEG measurements.

The Age-ility Project (Phase 1): Structural and functional imaging and electrophysiological data repository

Our understanding of the complex interplay between structural and functional organisation of brain networks is being advanced by the development of novel multi-modal analyses approaches. The Age-ility Project (Phase 1) data repository offers open access to structural MRI, diffusion MRI, and resting-state fMRI scans, as well as resting-state EEG recorded from the same community participants (n=131, 15-35 y, 66 male). Raw imaging and electrophysiological data as well as essential demographics are made available via the NITRC website. All data have been reviewed for artifacts using a rigorous quality control protocol and detailed case notes are provided.

Approximate average head models for EEG source imaging

We examine the performance of approximate models (AM) of the head in solving the EEG inverse problem. The AM are needed when the individual's MRI is not available. We simulate the electric potential distribution generated by cortical sources for a large sample of 305 subjects, and solve the inverse problem with AM. Statistical comparisons are carried out with the distribution of the localization errors. We propose several new AM. These are the average of many individual realistic MRI-based models, such as surface-based models or lead fields. We demonstrate that the lead fields of the AM should be calculated considering source moments not constrained to be normal to the cortex. We also show that the imperfect anatomical correspondence between all cortices is the most important cause of localization errors. Our average models perform better than a random individual model or the usual average model in the MNI space. We also show that a classification based on race and gender or head size before averaging does not significantly improve the results. Our average models are slightly better than an existing AM with shape guided by measured individual electrode positions, and have the advantage of not requiring such measurements. Among the studied models, the Average Lead Field seems the most convenient tool in large and systematical clinical and research studies demanding EEG source localization, when MRI are unavailable. This AM does not need a strict alignment between head models, and can therefore be easily achieved for any type of head modeling approach.

Validation of a method for coregistering scalp recording locations with 3D structural MR images

A common problem in brain imaging is how to most appropriately coregister anatomical and functional data sets into a common space. For surface-based recordings such as the event related optical signal (EROS), near-infrared spectroscopy (NIRS), event-related potentials (ERPs), and magnetoencephalography (MEG), alignment is typically done using either (1) a landmark-based method involving placement of surface markers that can be detected in both modalities; or (2) surface-fitting alignment that samples many points on the surface of the head in the functional space and aligns those points to the surface of the anatomical image. Here we compare these two approaches and advocate a combination of the two in order to optimize coregistration of EROS and NIRS data with structural magnetic resonance images (sMRI). Digitized 3D sensor locations obtained with a Polhemus digitizer can be effectively coregistered with sMRI using fiducial alignment as an initial guess followed by a Marquardt-Levenberg least-squares rigid-body transform (df = 6) to match the surfaces. Additional scaling parameters (df = 3) and point-by-point surface constraints can also be employed to further improve fitting. These alignment procedures place the lower-bound residual error at 1.3 +/- 0.1 mm (micro +/- s) and the upper-bound target registration error at 4.4 +/- 0.6 mm (micro +/- s). The dependence of such errors on scalp segmentation, number of registration points, and initial guess is also investigated. By optimizing alignment techniques, anatomical localization of surface recordings can be improved in individual subjects.

EEG to MRI registration based on global and local similarities of MRI intensity distributions

In this paper, a novel method for EEG to MRI registration is proposed. Initial registration is achieved by extracting and matching symmetry planes of MRI and EEG data, followed by iterative registration based on minimizing a cost function. Comparison of the intensity distributions of the whole MR image and MRI voxels around a head surface point yields global similarities, while the comparison of intensity distributions of MRI voxels around corresponding EEG points, which reflects the head's sagittal symmetry, yields local similarities. Therefore, when the EEG points are registered to the MR image, maximal global and local similarities should be obtained. The cost function, incorporating global and local similarities, was the sum of Kullback-Leibler divergences between corresponding intensity distributions. The proposed method was evaluated on clinical MRI data with simulated EEG data, yielding mean registration error of 0.48 +/- 0.33 mm, while with real EEG data an average root-mean-square point-to-surface error of 2.27 +/- 0.02 mm was obtained.

Improving the finite element forward model of the human head by warping using elastic deformation

As the use of realistic geometry in the forward model of electrical impedance tomography (EIT) of brain function appears to improve image reconstruction, the generation of patient-specific finite element meshes has been the subject of much recent work. This paper presents a more rapid method of generating more geometrically accurate finite element meshes of the human head by warping existing meshes such that the surface boundary beneath the electrodes closely matches that of the subject with minimal degradation to the quality of the mesh. Pre-existing meshes of spheres and adult head models incorporating key internal anatomical features are warped, using elastic deformation, to match a phantom latex tank incorporating a real skull. The algorithm is described and tests are carried out to optimize the key parameters to ensure minimal degradation of mesh quality and distortion of internal features. Results show that the algorithm operating with the optimum parameters produces meshes of sound quality and could represent an important step in the timely and productive creation of forward models in clinical applications.

Segmentation of skull and scalp in 3-D human MRI using mathematical morphology

We present a new technique for segmentation of skull and scalp in T(1)-weighted magnetic resonance images (MRIs) of the human head. Our method uses mathematical morphological operations to generate realistic models of the skull, scalp, and brain that are suitable for electroencephalography (EEG) and magnetoencephalography (MEG) source modeling. We first segment the brain using our Brain Surface Extractor algorithm; using this, we can ensure that the brain does not intersect our skull segmentation. We next generate a scalp mask using a combination of thresholding and mathematical morphology. We use the scalp mask in our skull segmentation procedure, as it allows us to automatically exclude background voxels with intensities similar to those of the skull. We find the inner and outer skull boundaries using thresholding and morphological operations. Finally, we mask the results with the scalp and brain volumes to ensure closed and nonintersecting skull boundaries. Visual evaluation indicated accurate segmentations of the cranium at a gross anatomical level (other than small holes in the zygomatic bone in eight subjects) in all 44 MRI volumes processed when run using default settings. In a quantitative comparison with coregistered CT images as a gold standard, MRI skull segmentation accuracy, as measured using the Dice coefficient, was found to be similar to that which would be obtained using CT imagery with a registration error of 2-3 mm.

Individual differences in brain activity during visuo-spatial processing assessed by slow cortical potentials and LORETA

Using slow-cortical potentials (SCPs), Vitouch et al. demonstrated that subjects with low ability to solve a complex visuo-spatial imagery task show higher activity in occipital, parietal and frontal cortex during task processing than subjects with high ability. This finding has been interpreted in the sense of the so-called "neural efficiency" hypothesis, which assumes that the central nervous system of individuals with higher intellectual abilities is functioning in a more efficient way than the one of individuals with lower abilities. Using a higher spatial resolution of SCP recordings, and by employing the source localization method of LORETA (low-resolution electromagnetic tomography), we investigated this hypothesis by performing an extended replication of Vitouch et al.'s study. SCPs during processing of a visuo-spatial imagery task were recorded in pre-selected subjects with either high or low abilities in solving the imagery task. Topographic and LORETA analyses of SCPs revealed that a distributed network of extrastriate occipital, superior parietal, temporal, medial frontal and prefrontal areas was active during task solving. This network is well in line with former studies of the functional neuroanatomy of visuo-spatial imagery. Contrary to our expectations, however, the results of Vitouch et al. as well as of other studies supporting the neural efficiency hypothesis could not be confirmed since no difference in brain activity between groups was observed. This inconsistency between studies might be due to differing task processing strategies. While subjects with high abilities in the Vitouch et al. study seemed to use a visuo-perceptual task solving approach, all other subjects relied upon a visuo-motor task processing strategy.

Registration and real-time visualization of transcranial magnetic stimulation with 3-D MR images

This paper describes a method for registering and visualizing in real-time the results of transcranial magnetic stimulations (TMS) in physical space on the corresponding anatomical locations in MR images of the brain. The method proceeds in three main steps. First, the patient scalp is digitized in physical space with a magnetic-field digitizer, following a specific digitization pattern. Second, a registration process minimizes the mean square distance between those points and a segmented scalp surface extracted from the magnetic resonance image. Following this registration, the physician can follow the change in coil position in real-time through the visualization interface and adjust the coil position to the desired anatomical location. Third, amplitude of motor evoked potentials can be projected onto the segmented brain in order to create functional brain maps. The registration has subpixel accuracy in a study with simulated data, while we obtain a point to surface root-mean-square error of 1.17+/-0.38 mm in a 24 subject study.

Reconstruction of 3-D Head Geometry From Digitized Point Sets: An Evaluation Study

In this paper, we evaluate different methods to estimate patient-specific scalp, skull, and brain surfaces from a set of digitized points from the target's scalp surface. The reconstruction problem is treated as a registration problem: An a priori surface model, consisting of the scalp, skull, and brain surfaces, is registered to the digitized surface points. The surface model is generated from segmented magnetic resonance (MR) volume images. We study both affine and free-form deformation (FFD) registration, the use of average models, the averaging of individual registration results, a model selection procedure, and statistical deformation models. The registration algorithms are mainly previously published, and the objective of this paper is to evaluate these methods in this particular application with sparse data. The main interest of this paper is to generate geometric head models for biomedical applications, such as electroencephalography and magnetoencephalographic. However, the methods can also be applied to other anatomical regions and to other application areas. The methods were validated using 15 MR volume images, from which the scalp, skull, and brain were manually segmented. The best results were achieved by averaging the results of the FFD registrations of the database: the mean distance from the manually segmented target surface to a deformed a priori model surface for the studied anatomical objects was 1.68-2.08 mm, depending on the point set used. The results support the use of the evaluated methods for the reconstruction of geometric models in applications with sparse data.