Estimation of the electric conductivity from scalp measurements: feasibility and application to source localization

Contrast enhancement in EIT imaging of the brain

We consider electrical impedance tomography (EIT) imaging of the brain. The brain is surrounded by the poorly conducting skull which has low conductivity compared to the brain. The skull layer causes a partial shielding effect which leads to weak sensitivity for the imaging of the brain tissue. In this paper we propose an approach based on the Bayesian approximation error approach, to enhance the contrast in brain imaging. With this approach, both the (uninteresting) geometry and the conductivity of the skull are embedded in the approximation error statistics, which leads to a computationally efficient algorithm that is able to detect features such as internal haemorrhage with significantly increased sensitivity and specificity. We evaluate the approach with simulations and phantom data.

Development of volume conductor and source models to localize epileptic foci

SUMMARY

There is increasing interest in mapping and source reconstruction from electrocorticoencephalographic (ECoG) grid data and comparison to surface EEG evaluations of epileptic patients. ECoG mapping onto three-dimensional renderings of the individual cortical anatomy derived from magnetic resonance images and computed tomography (CT) is performed after coregistration of anatomical and functional coordinate systems. Source reconstructions from ECoG and EEG are compared using different source models and realistically shaped volume conductor models. Realistically shaped volume conductor models for EEG source reconstruction are a prerequisite for improved localization accuracy. Individual boundary element method (BEM) models derived from MRI represent the "gold standard" and can approximate isotropic homogeneous head compartments and thus give an improved description of the head shape as compared with classical oversimplifying spherical shell models. Anisotropic volume conduction properties of the bone layer or the white matter fibers can be described by the finite element method (FEM); unfortunately, these models require a huge computational effort and are thus not used in daily applications. To avoid this computational effort, head models derived from an averaged MRI dataset can be used. Highly refined models with a large number of nodes and thus better numerical accuracy can be used by this approach, because the setup is performed only once and the decomposed models or precomputed leadfield matrices are saved for later application. Individual image data are not at all needed, if an overlay of the reconstruction results with the anatomy is not desired. With precomputed leadfield matrices and linear interpolation techniques, at least standardized BEM and FEM volume conductor models derived from averaged MRI datasets can achieve the same computational speed as analytical spherical models. The smoothed cortical envelope is used as a realistically shaped single-shell volume conductor model for ECoG source reconstruction, whereas three-compartment BEM-models are required for EEG. The authors describe how to localize ECoG-grid electrode positions and how to segment the cortical surface from coregistered magnetic resonance and CT images. Landmark-based coregistration is performed using common fiducials in both image modalities. Another more promising automatic approach is based on mutual three-dimensional volume gray-level information. The ECoG electrode positions can be retrieved from three-dimensional CT slices manually using cursors in thresholded images with depth information. Alternatively, the smoothed envelope of the cortical surface segmented from the MRI is used to semiautomatically determine the grid electrode positions by marking the four corners and measuring distances along the smoothed surface. With extended source patches for cortically constrained scans and current density reconstructions, results from ECoG and surface EEG data were compared. Single equivalent dipoles were used to explain the EEG far fields, and results were compared with the original current density distributions explaining the ECoG data. The authors studied the performance of spherical and realistically shaped BEM volume conductor models for EEG and ECoG source reconstruction in spherical and nonspherical parts of the head with simulations and measured epileptic spike data. Only small differences between spherical and realistically shaped models were found in the spherical parts of the head, whereas realistically shaped models are superior to spherical approximations in both single-shell ECoG and three-shell EEG cases in the nonspherical parts, such as the temporal lobe areas. The ECoG near field is more complicated to interpret than the surface EEG far field and cannot be explained in general by simple equivalent dipoles. However, from simulations with realistically shaped volume conductor models and cortically constrained source models, the authors studied how the bone and skin layer act as spatial low pass filters that smooth and simplify the surface EEG maps generated by much more complicated-looking source configurations derived from measured ECoG data.

Use of the isolated problem approach for multi-compartment BEM models of electro-magnetic source imaging

The isolated problem approach (IPA) is a method used in the boundary element method (BEM) to overcome numerical inaccuracies caused by the high-conductivity difference in the skull and the brain tissues in the head. Hämäläinen and Sarvas (1989 IEEE Trans. Biomed. Eng. 36 165-71) described how the source terms can be updated to overcome these inaccuracies for a three-layer head model. Meijs et al (1989 IEEE Trans. Biomed. Eng. 36 1038-49) derived the integral equations for the general case where there are an arbitrary number of layers inside the skull. However, the IPA is used in the literature only for three-layer head models. Studies that use complex boundary element head models that investigate the inhomogeneities in the brain or model the cerebrospinal fluid (CSF) do not make use of the IPA. In this study, the generalized formulation of the IPA for multi-layer models is presented in terms of integral equations. The discretized version of these equations are presented in two different forms. In a previous study (Akalin-Acar and Gençer 2004 Phys. Med. Biol. 49 5011-28), we derived formulations to calculate the electroencephalography and magnetoencephalography transfer matrices assuming a single layer in the skull. In this study, the transfer matrix formulations are updated to incorporate the generalized IPA. The effects of the IPA are investigated on the accuracy of spherical and realistic models when the CSF layer and a tumour tissue are included in the model. It is observed that, in the spherical model, for a radial dipole 1 mm close to the brain surface, the relative difference measure (RDM*) drops from 1.88 to 0.03 when IPA is used. For the realistic model, the inclusion of the CSF layer does not change the field pattern significantly. However, the inclusion of an inhomogeneity changes the field pattern by 25% for a dipole oriented towards the inhomogeneity. The effect of the IPA is also investigated when there is an inhomogeneity in the brain. In addition to a considerable change in the scale of the potentials, the field pattern also changes by 15%. The computation times are presented for the multi-layer realistic head model.

Evaluation of different cortical source localization methods using simulated and experimental EEG data

Different cortical source localization methods have been developed to directly link the scalp potentials with the cortical activities. Up to now, these methods are the only possible solution to noninvasively investigate cortical activities with both high spatial and time resolutions. However, the application of these methods is hindered by the fact that they have not been rigorously evaluated nor compared. In this paper, the performances of several source localization methods (moving dipoles, minimum Lp norm, and low resolution tomography (LRT) with Lp norm, p equal to 1, 1.5, and 2) were evaluated by using simulated scalp EEG data, scalp somatosensory evoked potentials (SEPs), and upper limb motor-related potentials (MRPs) obtained on human subjects (all with 163 scalp electrodes). By using simulated EEG data, we first evaluated the source localization ability of the above methods quantitatively. Subsequently, the performance of the various methods was evaluated qualitatively by using experimental SEPs and MRPs. Our results show that the overall LRT Lp norm method with p equal to 1 has a better source localization ability than any of the other investigated methods and provides physiologically meaningful reconstruction results. Our evaluation results provide useful information for choosing cortical source localization approaches for future EEG/MEG studies.

Effects of skull thickness, anisotropy, and inhomogeneity on forward EEG/ERP computations using a spherical three-dimensional resistor mesh model

Bone thickness, anisotropy, and inhomogeneity have been reported to induce important variations in electroencephalogram (EEG) scalp potentials. To study this effect, we used an original three-dimensional (3-D) resistor mesh model described in spherical coordinates, consisting of 67,464 elements and 22,105 nodes arranged in 36 different concentric layers. After validation of the model by comparison with the analytic solution, potential variations induced by geometric and electrical skull modifications were investigated at the surface in the dipole plane and along the dipole axis, for several eccentricities and bone thicknesses. The resistor mesh permits one to obtain various configurations, as local modifications are introduced very easily. This has allowed several head models to be designed to study the effects of skull properties (thickness, anisotropy, and heterogeneity) on scalp surface potentials. Results show a decrease of potentials in bone, depending on bone thickness, and a very small decrease through the scalp layer. Nevertheless, similar scalp potentials can be obtained using either a thick scalp layer and a thin skull layer, and vice versa. It is thus important to take into account skull and scalp thicknesses, because the drop of potential in bone depends on both. The use of three different layers for skull instead of one leads to small differences in potential values and patterns. In contrast, the introduction of a hole in the skull highly increases the maximum potential value (by a factor of 11.5 in our case), because of the absence of potential drop in the corresponding volume. The inverse solution without any a priori knowledge indicates that the model with the hole gives the largest errors in both position and dipolar moment. Our results indicate that the resistor mesh model can be used as a robust and user-friendly simulation tool in EEG or event-related potentials. It makes it possible to build up real head models directly from anatomic magnetic resonance imaging without tessellation, and is able to take into account head heterogeneities very simply by changing volume elements conductivity. Hum. Brain Mapping 21:84-95, 2004.

Current constrained voltage scaled reconstruction (CCVSR) algorithm for MR-EIT and its performance with different probing current patterns

Conventional injected-current electrical impedance tomography (EIT) and magnetic resonance imaging (MRI) techniques can be combined to reconstruct high resolution true conductivity images. The magnetic flux density distribution generated by the internal current density distribution is extracted from MR phase images. This information is used to form a fine detailed conductivity image using an Ohm's law based update equation. The reconstructed conductivity image is assumed to differ from the true image by a scale factor. EIT surface potential measurements are then used to scale the reconstructed image in order to find the true conductivity values. This process is iterated until a stopping criterion is met. Several simulations are carried out for opposite and cosine current injection patterns to select the best current injection pattern for a 2D thorax model. The contrast resolution and accuracy of the proposed algorithm are also studied. In all simulation studies, realistic noise models for voltage and magnetic flux density measurements are used. It is shown that, in contrast to the conventional EIT techniques, the proposed method has the capability of reconstructing conductivity images with uniform and high spatial resolution. The spatial resolution is limited by the larger element size of the finite element mesh and twice the magnetic resonance image pixel size.

Interictal spike EEG source analysis in hypothalamic hamartoma epilepsy

OBJECTIVE

The epilepsy associated with the hypothalamic hamartomas constitutes a syndrome with peculiar seizures, usually refractory to medical therapy, mild cognitive delay, behavioural problems and multifocal spike activity in the scalp electroencephalogram (EEG). The cortical origin of spikes has been widely assumed but not specifically demonstrated.

METHODS

We present results of a source analysis of interictal spikes from 4 patients (age 2-25 years) with epilepsy and hypothalamic hamartoma, using EEG scalp recordings (32 electrodes) and realistic boundary element models constructed from volumetric magnetic resonance imaging (MRIs). Multifocal spike activity was the most common finding, distributed mainly over the frontal and temporal lobes. A spike classification based on scalp topography was done and averaging within each class performed to improve the signal to noise ratio. Single moving dipole models were used, as well as the Rap-MUSIC algorithm.

RESULTS

All spikes with good signal to noise ratio were best explained by initial deep sources in the neighbourhood of the hamartoma, with late sources located in the cortex. Not a single patient could have his spike activity explained by a combination of cortical sources.

CONCLUSIONS

Overall, the results demonstrate a consistent origin of spike activity in the subcortical region in the neighbourhood of the hamartoma, with late spread to cortical areas.

Weighted regularization in electrical impedance tomography with applications to acute cerebral stroke

We apply electrical impedance tomography to detect and localize brain impedance changes associated with stroke. Forward solutions are computed using the finite-element method in two dimensions. We assume that baseline conductivity values are known for the major head tissues, and focus on changes in the brain compartment only. We use singular-value decomposition (SVD) to show that different impedance measurement patterns, which are theoretically equivalent by the reciprocity theorem, have different sensitivities to the brain compartment in the presence of measurement noise. The inverse problem is solved in part by standard means, using iterated SVD, and regularizing by truncation. To improve regularization we introduce a weighting scheme which normalizes the sensitivity matrix for voxels at different depths. This increases the number of linearly independent components which contribute to the solution, and forces the different measurement patterns to have similar sensitivity. When applied to stroke, this weighted regularization improves image quality overall.

A standardized boundary element method volume conductor model

OBJECTIVES

We used a 3-compartment boundary element method (BEM) model from an averaged magnetic resonance image (MRI) data set (Montreal Neurological Institute) in order to provide simple access to realistically shaped volume conductor models for source reconstruction, as compared to individually derived models. The electrode positions were transformed into the model's coordinate system, and the best fit dipole results were transformed back to the original coordinate system. The localization accuracy of the new approach was tested in a comparison with simulated data and with individual BEM models of epileptic spike data from several patients.

METHODS

The standard BEM model consisted of a total of 4770 nodes, which describe the smoothed cortical envelope, the outside of the skull, and the outside of the skin. The electrode positions were transformed to the model coordinate system by using 3-5 fiducials (nasion, left and right preauricular points, vertex, and inion). The transformation consisted of an averaged scaling factor and a rigid transformation (translation and rotation). The potential values at the transformed electrode positions were calculated by linear interpolation from the stored transfer matrix of the outer BEM compartment triangle net. After source reconstruction the best fit dipole results were transformed back into the original coordinate system by applying the inverse of the first transformation matrix.

RESULTS

Test-dipoles at random locations and with random orientations inside of a highly refined reference BEM model were used to simulate noise-free data. Source reconstruction results using a spherical and the standardized BEM volume conductor model were compared to the known dipole positions. Spherical head models resulted in mislocation errors at the base of the brain. The standardized BEM model was applied to averaged and unaveraged epileptic spike data from 7 patients. Source reconstruction results were compared to those achieved by 3 spherical shell models and individual BEM models derived from the individual MRI data sets. Similar errors to that evident with simulations were noted with spherical head models. Standardized and individualized BEM models were comparable.

CONCLUSIONS

This new approach to head modeling performed significantly better than a simple spherical shell approximation, especially in basal brain areas, including the temporal lobe. By using a standardized head for the BEM setup, it offered an easier and faster access to realistically shaped volume conductor models as compared to deriving specific models from individual 3-dimensional MRI data.

Modeling of post-surgical brain and skull defects in the EEG inverse problem with the boundary element method

OBJECTIVES

In order to obtain accurate EEG inverse solutions in patients subjected to surgery, we have studied the feasibility and influence of incorporating brain and skull defects in realistic head models.

METHODS

We first measured the conductivity of the methacrylate used for cranioplasty. Then, we designed realistic boundary element method head models with a skull burr hole, a methacrylate plug or a temporal-lobe resection. We simulated the potentials that would be produced at 71 electrode locations (10/10 system) by dipoles located near the defects. Then, we fitted dipoles on these potentials using a defect-free head model. We also ran simulations in a noisy situation and with higher skull and cerebrospinal fluid (CSF) conductivity.

RESULTS

The largest errors were found for burr holes, with a localization error up to 20 mm for a radial dipole located 30 mm below the hole and an amplification factor of 8. Methacrylate plugs lead to errors up to 5 mm and 0.5; the resection only lead to errors of 2 mm and 1.3. Results obtained with noise were consistent with those obtained without noise. Doubling the skull conductivity led to errors that were reduced by 10-20%, while doubling CSF conductivity increased the errors by up to 31%.

CONCLUSIONS

We have shown that it is important to incorporate skull defects in realistic head models when sources are located near the defects and precision is sought. Brain cavities of the size of a typical anterior temporal lobe resection may be omitted without a significant impact on dipole localization.

Boundary element method volume conductor models for EEG source reconstruction

OBJECTIVES

The boundary element method (BEM) approximates the different compartments of volume conductor models by closed triangle meshes with a limited number of nodes. The shielding effect of the weakly conducting skull layer of the human head leads to decreasing potential gradients from the inside to the outside. Thus, there may be an optimum distribution of nodes to the compartments for a given number of nodes corresponding to a fixed computational effort, resulting in improved accuracy as compared to standard uniform distributions.

METHODS

Spherical and realistically shaped surfaces are approximated by 500, 1000, 2000, and 3000 nodes, each leading to BEM models with 1500-9000 nodes in total. Electrodes are placed on extended 10/20-system positions. Potential distributions of test-dipoles at 4000 random positions within the innermost compartment are calculated. Dipoles are then fitted using 192 different models to find the optimum node distribution.

RESULTS

Fitted dipole positions for all BEM models are evaluated to show the dependency of the averaged and maximum localization errors on their node distributions. Dipoles close to the innermost boundary exhibit the largest localization errors, which mainly depend on the refinement of this compartment's triangle mesh.

CONCLUSIONS

More than 500 nodes per compartment are needed for reliable BEM models. For a state-of-the-art model consisting of 6000 nodes overall, the best model consists of 3000, 2000, and 1000 nodes from the inside to the outside.