Targeted density electrode placement achieves high concordance with traditional high-density EEG for electrical source imaging in epilepsy.

OBJECTIVE: High-density (HD) electroencephalography (EEG) is increasingly used in presurgical epilepsy evaluation, but it is demanding in time and resources. To overcome these issues, we compared EEG source imaging (ESI) solutions with a targeted density and HD-EEG montage. METHODS: HD-EEGs from patients undergoing presurgical evaluation were analyzed. A low-density recording was created by selecting the 25 electrodes of a standard montage from the 83 electrodes of the HD-EEG and adding 8-11 electrodes around the electrode with the highest amplitude interictal epileptiform discharges. The ESI solution from this "targeted" montage was compared to that from the HD-EEG using the distance between peak vertices, sublobar concordance and a qualitative similarity measure. RESULTS: Fifty-eight foci of forty-three patients were included. The median distance between the peak vertices of the two montages was 13.2 mm, irrespective of focus' location. Tangential generators (n = 5/58) showed a higher distance than radial generators (p = 0.04). We found sublobar concordance in 54/58 of the foci (93%). Map similarity, assessed by an epileptologist, had a median score of 4/5. CONCLUSIONS: ESI solutions obtained from a targeted density montage show high concordance with those calculated from HD-EEG. SIGNIFICANCE: Requiring significantly fewer electrodes, targeted density EEG allows obtaining similar ESI solutions as traditional HD-EEG montage.

Neurophysiology, Neuropsychology, and Epilepsy, in 2022: Hills We Have Climbed and Hills Ahead. Neurophysiology in epilepsy.

Since the discovery of the human electroencephalogram (EEG), neurophysiology techniques have become indispensable tools in our armamentarium to localize epileptic seizures. New signal analysis techniques and the prospects of artificial intelligence and big data will offer unprecedented opportunities to further advance the field in the near future, ultimately resulting in improved quality of life for many patients with drug-resistant epilepsy. This article summarizes selected presentations from Day 1 of the two-day symposium "Neurophysiology, Neuropsychology, Epilepsy, 2022: Hills We Have Climbed and the Hills Ahead". Day 1 was dedicated to highlighting and honoring the work of Dr. Jean Gotman, a pioneer in EEG, intracranial EEG, simultaneous EEG/ functional magnetic resonance imaging, and signal analysis of epilepsy. The program focused on two main research directions of Dr. Gotman, and was dedicated to "High-frequency oscillations, a new biomarker of epilepsy" and "Probing the epileptic focus from inside and outside". All talks were presented by colleagues and former trainees of Dr. Gotman. The extended summaries provide an overview of historical and current work in the neurophysiology of epilepsy with emphasis on novel EEG biomarkers of epilepsy and source imaging and concluded with an outlook on the future of epilepsy research, and what is needed to bring the field to the next level.

Deconvolution of hemodynamic responses along the cortical surface using personalized functional near infrared spectroscopy.

In functional near infrared spectroscopy (fNIRS), deconvolution analysis of oxy and deoxy-hemoglobin concentration changes allows estimating specific hemodynamic response functions (HRF) elicited by neuronal activity, taking advantage of the fNIRS excellent temporal resolution. Diffuse optical tomography (DOT) is also becoming the new standard reconstruction procedure as it is more accurate than the modified Beer Lambert law approach at the sensor level. The objective of this study was to assess the relevance of HRF deconvolution after DOT constrained along the cortical surface. We used local personalized fNIRS montages which consists in optimizing the position of fNIRS optodes to ensure maximal sensitivity to subject specific target brain regions. We carefully evaluated the accuracy of deconvolution when applied after DOT, using realistic simulations involving several HRF models at different signal to noise ratio (SNR) levels and on real data related to motor and visual tasks in healthy subjects and from spontaneous pathological activity in one patient with epilepsy. We demonstrated that DOT followed by deconvolution was able to accurately recover a large variability of HRFs over a large range of SNRs. We found good performances of deconvolution analysis for SNR levels usually encountered in our applications and we were able to reconstruct accurately the temporal dynamics of HRFs in real conditions.

Differences in MEG and EEG power-law scaling explained by a coupling between spatial coherence and frequency: a simulation study.

Electrophysiological signals (electroencephalography, EEG, and magnetoencephalography, MEG), as many natural processes, exhibit scale-invariance properties resulting in a power-law (1/f) spectrum. Interestingly, EEG and MEG differ in their slopes, which could be explained by several mechanisms, including non-resistive properties of tissues. Our goal in the present study is to estimate the impact of space/frequency structure of source signals as a putative mechanism to explain spectral scaling properties of neuroimaging signals. We performed simulations based on the summed contribution of cortical patches with different sizes (ranging from 0.4 to 104.2 cm2). Small patches were attributed signals of high frequencies, whereas large patches were associated with signals of low frequencies, on a logarithmic scale. The tested parameters included i) the space/frequency structure (range of patch sizes and frequencies) and ii) the amplitude factor c parametrizing the spatial scale ratios. We found that the space/frequency structure may cause differences between EEG and MEG scale-free spectra that are compatible with real data findings reported in previous studies. We also found that below a certain spatial scale, there were no more differences between EEG and MEG, suggesting a limit for the resolution of both methods.Our work provides an explanation of experimental findings. This does not rule out other mechanisms for differences between EEG and MEG, but suggests an important role of spatio-temporal structure of neural dynamics. This can help the analysis and interpretation of power-law measures in EEG and MEG, and we believe our results can also impact computational modeling of brain dynamics, where different local connectivity structures could be used at different frequencies.

Optimal positioning of optodes on the scalp for personalized functional near-infrared spectroscopy investigations.

BACKGROUND: Application of functional Near InfraRed Spectroscopy (fNIRS) in neurology is still limited as a good optical coupling and optimized optode coverage of specific brain regions remains challenging, notably for prolonged monitoring. METHODS: We propose to evaluate a new procedure allowing accurate investigation of specific brain regions. The procedure consists in: (i) A priori maximization of spatial sensitivity of fNIRS measurements targeting specific brain regions, while reducing the number of applied optodes in order to decrease installation time and improve subject comfort. (ii) Utilization of a 3D neuronavigation device and usage of collodion to glue optodes on the scalp, ensuring good optical contact for prolonged investigations. (iii) Local reconstruction of the hemodynamic activity along the cortical surface using inverse modelling. RESULTS: Using realistic simulations, we demonstrated that maps derived from optimal montage acquisitions showed, after reconstruction, spatial resolution only slightly lower to that of ultra high density montages while significantly reducing the number of optodes. The optimal montages provided overall good quantitative accuracy especially at the peak of the spatially reconstructed map. We also evaluated real motor responses in two healthy subjects and obtained reproducible motor responses over different sessions. COMPARISON WITH EXISTING METHODS: We are among the first to propose a mathematical optimization strategy, allowing high sensitivity measurements. CONCLUSIONS: Our results support that using personalized optimal montages should allow to conduct accurate fNIRS studies in clinical settings and realistic lifestyle conditions.

Comparison of the spatial resolution of source imaging techniques in high-density EEG and MEG.

BACKGROUND: The present study aims at evaluating and comparing electrical and magnetic distributed source imaging methods applied to high-density Electroencephalography (hdEEG) and Magnetoencephalography (MEG) data. We used resolution matrices to characterize spatial resolution properties of Minimum Norm Estimate (MNE), dynamic Statistical Parametric Mapping (dSPM), standardized Low-Resolution Electromagnetic Tomography (sLORETA) and coherent Maximum Entropy on the Mean (cMEM, an entropy-based technique). The resolution matrix provides information of the Point Spread Functions (PSF) and of the Crosstalk functions (CT), this latter being also called source leakage, as it reflects the influence of a source on its neighbors. METHODS: The spatial resolution of the inverse operators was first evaluated theoretically and then with real data acquired using electrical median nerve stimulation on five healthy participants. We evaluated the Dipole Localization Error (DLE) and the Spatial Dispersion (SD) of each PSF and CT map. RESULTS: cMEM showed the smallest spatial spread (SD) for both PSF and CT maps, whereas localization errors (DLE) were similar for all methods. Whereas cMEM SD values were lower in MEG compared to hdEEG, the other methods slightly favored hdEEG over MEG. In real data, cMEM provided similar localization error and significantly less spatial spread than other methods for both MEG and hdEEG. Whereas both MEG and hdEEG provided very accurate localizations, all the source imaging methods actually performed better in MEG compared to hdEEG according to all evaluation metrics, probably due to the higher signal-to-noise ratio of the data in MEG. CONCLUSION: Our overall results show that all investigated methods provide similar localization errors, suggesting very accurate localization for both MEG and hdEEG when similar number of sensors are considered for both modalities. Intrinsic properties of source imaging methods as well as their behavior for well-controlled tasks, suggest an overall better performance of cMEM in regards to spatial resolution and spatial leakage for both hdEEG and MEG. This indicates that cMEM would be a good candidate for studying source localization of focal and extended generators as well as functional connectivity studies.

Zoomed MRI Guided by Combined EEG/MEG Source Analysis: A Multimodal Approach for Optimizing Presurgical Epilepsy Work-up and its Application in a Multi-focal Epilepsy Patient Case Study.

In recent years, the use of source analysis based on electroencephalography (EEG) and magnetoencephalography (MEG) has gained considerable attention in presurgical epilepsy diagnosis. However, in many cases the source analysis alone is not used to tailor surgery unless the findings are confirmed by lesions, such as, e.g., cortical malformations in MRI. For many patients, the histology of tissue resected from MRI negative epilepsy shows small lesions, which indicates the need for more sensitive MR sequences. In this paper, we describe a technique to maximize the synergy between combined EEG/MEG (EMEG) source analysis and high resolution MRI. The procedure has three main steps: (1) construction of a detailed and calibrated finite element head model that considers the variation of individual skull conductivities and white matter anisotropy, (2) EMEG source analysis performed on averaged interictal epileptic discharges (IED), (3) high resolution (0.5 mm) zoomed MR imaging, limited to small areas centered at the EMEG source locations. The proposed new diagnosis procedure was then applied in a particularly challenging case of an epilepsy patient: EMEG analysis at the peak of the IED coincided with a right frontal focal cortical dysplasia (FCD), which had been detected at standard 1 mm resolution MRI. Of higher interest, zoomed MR imaging (applying parallel transmission, 'ZOOMit') guided by EMEG at the spike onset revealed a second, fairly subtle, FCD in the left fronto-central region. The evaluation revealed that this second FCD, which had not been detectable with standard 1 mm resolution, was the trigger of the seizures.

Complex patterns of spatially extended generators of epileptic activity: Comparison of source localization methods cMEM and 4-ExSo-MUSIC on high resolution EEG and MEG data.

Electric Source Imaging (ESI) and Magnetic Source Imaging (MSI) of EEG and MEG signals are widely used to determine the origin of interictal epileptic discharges during the pre-surgical evaluation of patients with epilepsy. Epileptic discharges are detectable on EEG/MEG scalp recordings only when associated with a spatially extended cortical generator of several square centimeters, therefore it is essential to assess the ability of source localization methods to recover such spatial extent. In this study we evaluated two source localization methods that have been developed for localizing spatially extended sources using EEG/MEG data: coherent Maximum Entropy on the Mean (cMEM) and 4th order Extended Source Multiple Signal Classification (4-ExSo-MUSIC). In order to propose a fair comparison of the performances of the two methods in MEG versus EEG, this study considered realistic simulations of simultaneous EEG/MEG acquisitions taking into account an equivalent number of channels in EEG (257 electrodes) and MEG (275 sensors), involving a biophysical computational neural mass model of neuronal discharges and realistically shaped head models. cMEM and 4-ExSo-MUSIC were evaluated for their sensitivity to localize complex patterns of epileptic discharges which includes (a) different locations and spatial extents of multiple synchronous sources, and (b) propagation patterns exhibited by epileptic discharges. Performance of the source localization methods was assessed using a detection accuracy index (Area Under receiver operating characteristic Curve, AUC) and a Spatial Dispersion (SD) metric. Finally, we also presented two examples illustrating the performance of cMEM and 4-ExSo-MUSIC on clinical data recorded using high resolution EEG and MEG. When simulating single sources at different locations, both 4-ExSo-MUSIC and cMEM exhibited excellent performance (median AUC significantly larger than 0.8 for EEG and MEG), whereas, only for EEG, 4-ExSo-MUSIC showed significantly larger AUC values than cMEM. On the other hand, cMEM showed significantly lower SD values than 4-ExSo-MUSIC for both EEG and MEG. When assessing the impact of the source spatial extent, both methods provided consistent and reliable detection accuracy for a wide range of source spatial extents (source sizes ranging from 3 to 20cm2 for MEG and 3 to 30cm2 for EEG). For both EEG and MEG, 4-ExSo-MUSIC localized single source of large signal-to-noise ratio better than cMEM. In the presence of two synchronous sources, cMEM was able to distinguish well the two sources (their location and spatial extent), while 4-ExSo-MUSIC only retrieved one of them. cMEM was able to detect the spatio-temporal propagation patterns of two synchronous activities while 4-ExSo-MUSIC favored the strongest source activity. Overall, in the context of localizing sources of epileptic discharges from EEG and MEG data, 4-ExSo-MUSIC and cMEM were found accurately sensitive to the location and spatial extent of the sources, with some complementarities. Therefore, they are both eligible for application on clinical data.

Wavelet-based localization of oscillatory sources from magnetoencephalography data.

Transient brain oscillatory activities recorded with Eelectroencephalography (EEG) or magnetoencephalography (MEG) are characteristic features in physiological and pathological processes. This study is aimed at describing, evaluating, and illustrating with clinical data a new method for localizing the sources of oscillatory cortical activity recorded by MEG. The method combines time-frequency representation and an entropic regularization technique in a common framework, assuming that brain activity is sparse in time and space. Spatial sparsity relies on the assumption that brain activity is organized among cortical parcels. Sparsity in time is achieved by transposing the inverse problem in the wavelet representation, for both data and sources. We propose an estimator of the wavelet coefficients of the sources based on the maximum entropy on the mean (MEM) principle. The full dynamics of the sources is obtained from the inverse wavelet transform, and principal component analysis of the reconstructed time courses is applied to extract oscillatory components. This methodology is evaluated using realistic simulations of single-trial signals, combining fast and sudden discharges (spike) along with bursts of oscillating activity. The method is finally illustrated with a clinical application using MEG data acquired on a patient with a right orbitofrontal epilepsy.

Detection of hemodynamic responses to epileptic activity using simultaneous Electro-EncephaloGraphy (EEG)/Near Infra Red Spectroscopy (NIRS) acquisitions.

Simultaneous recordings of Electro-EncephaloGraphy (EEG) with Near InfraRed Spectroscopy (NIRS) allow measuring hemodynamic changes (changes in the concentration of oxy- and deoxyhemoglobin) at the time of epileptic discharges detected on scalp EEG. Two NIRS detection methods based on the General Linear Model (GLM) respectively in the time domain and in the time-frequency domain are investigated in this study using realistic simulations of spontaneous interictal epileptic activity. We evaluated the sensitivity at different Signal to Noise Ratios (SNR), the effect of either a large or a small number of discharges and the impact of model misspecification (e.g. omission or false detection of epileptic discharges). We also explored the effect on the estimation of key parameters, which set the model order. Simulations showed that both methods become inaccurate in lower SNR conditions, leading to many false positive detections. However, the time-frequency estimator showed better performance than the time-domain one. Key parameters for each algorithm were identified and results suggest to model confounds in the GLM differently for oxy- and deoxyhemoglobin. We also demonstrated that an inaccurate marking of epileptic events has a small impact on the detection statistics whereas an inaccurate specification of the hemodynamic response function delay decreases drastically the detection abilities. Finally, we illustrated the two methods on clinical EEG/NIRS data of one patient with focal epilepsy, showing an increase of regional Cerebral Blood Volume (rCBV) spatially concordant with the presumed epileptogenic focus.

Concordance between distributed EEG source localization and simultaneous EEG-fMRI studies of epileptic spikes.

In order to analyze where epileptic spikes are generated, we assessed the level of concordance between EEG source localization using distributed source models and simultaneous EEG-fMRI which measures the hemodynamic correlates of EEG activity. Data to be compared were first estimated on the same cortical surface and two comparison strategies were used: (1) MEM-concordance: a comparison between EEG sources localized with the Maximum Entropy on the Mean (MEM) method and fMRI clusters showing a significant hemodynamic response. Minimal geodesic distances between local extrema and overlap measurements between spatial extents of EEG sources and fMRI clusters were used to quantify MEM-concordance. (2) fMRI-relevance: estimation of the fMRI-relevance index alpha quantifying if sources located in an fMRI cluster could explain some scalp EEG data, when this fMRI cluster was used to constrain the EEG inverse problem. Combining MEM-concordance and fMRI-relevance (alpha) indexes, each fMRI cluster showing a significant hemodynamic response (p<0.05 corrected) was classified according to its concordance with EEG data. Nine patients with focal epilepsy who underwent EEG-fMRI examination followed by EEG recording outside the scanner were selected for this study. Among the 62 fMRI clusters analyzed (7 patients), 15 (24%) found in 6 patients were highly concordant with EEG according to both MEM-concordance and fMRI-relevance. EEG concordance was found for 5 clusters (8%) according to alpha only, suggesting sources missed by the MEM. No concordance with EEG was found for 30 clusters (48%) and for 10 clusters (16%) alpha was significantly negative, suggesting EEG-fMRI discordance. We proposed two complementary strategies to assess and classify EEG-fMRI concordance. We showed that for most patients, part of the hemodynamic response to spikes was highly concordant with EEG sources, whereas other fMRI clusters in response to the same spikes were found distant or discordant with EEG sources.

Anatomically informed interpolation of fMRI data on the cortical surface.

Analyzing functional magnetic resonance imaging (fMRI) data restricted to the cortical surface is of particular interest for two reasons: (1) to increase detection sensitivity using anatomical constraints and (2) to compare or use fMRI results in the context of source localization from magneto/electro-encephalography (MEEG) data, which requires data to be projected on the same spatial support. Designing an optimal scheme to interpolate fMRI raw data or resulting activation maps on the cortical surface relies on a trade-off between choosing large enough interpolation kernels, because of the distributed nature of the hemodynamic response, and avoiding mixing data issued from different anatomical structures. We propose an original method that automatically adjusts the level of such a trade-off, by defining interpolation kernels around each vertex of the cortical surface using a geodesic Voronoï diagram. This Voronoï-based interpolation method was evaluated using simulated fMRI activation maps, manually generated on an anatomical MRI, and compared with a more standard approach where interpolation kernels were defined as local spheres of radius r=3 or 5 mm. Several validation parameters were considered: the spatial resolution of the simulated activation map, the spatial resolution of the cortical mesh, the level of anatomical/functional data misregistration and the location of the vertices within the gray matter ribbon. Using an activation map at the spatial resolution of standard fMRI data, robustness to misregistration errors was observed for both methods, whereas only the Voronoï-based approach was insensitive to the position of the vertices within the gray matter ribbon.

Widespread and intense BOLD changes during brief focal electrographic seizures.

BACKGROUND: Combined recording of EEG and fMRI has shown changes in blood oxygenation level dependent (BOLD) signal during focal interictal epileptic spikes. Due to difficult assessment of seizures inside the scanner little is known about BOLD changes during seizures. OBJECTIVES: To describe BOLD changes related to brief focal electrographic seizures in a patient with right temporo-parietal gray matter nodular heterotopia. METHODS: The patient underwent two EEG-fMRI sessions during which several focal seizures were recorded. EEG was acquired continuously during scanning and seizure timing was used for statistical analysis. Functional maps were thresholded to disclose positive (activation) and negative (deactivation) BOLD changes. RESULTS: Twenty-five focal electrographic seizures were analyzed, consisting of runs of polyspikes lasting 2 to 6 s in the right temporal region. Activation included a large volume, involving the heterotopia and the abnormal temporo-parietal cortex overlying the nodule, with a clear maximum over the angular gyrus. Deactivation was bilateral and maximum in the occipital regions. The hemodynamic response function showed a return to baseline of the BOLD signal 30 s after seizure end. CONCLUSIONS: The brief focal seizures resulted in high amplitude and widespread blood oxygenation level dependent (BOLD) responses taking 30 s to return to baseline. This suggests that such brief events could have important behavioral consequences despite absent overt manifestations. A clear focal BOLD peak was found at some distance from the main EEG discharge, raising the possibility that the seizure could have started in a region that did not generate a visible EEG change despite its superficial location.

Evaluation of EEG localization methods using realistic simulations of interictal spikes.

Performing an accurate localization of sources of interictal spikes from EEG scalp measurements is of particular interest during the presurgical investigation of epilepsy. The purpose of this paper is to study the ability of six distributed source localization methods to recover extended sources of activated cortex. Due to the frequent lack of a gold standard to evaluate source localization methods, our evaluation was performed in a controlled environment using realistic simulations of EEG interictal spikes, involving several anatomical locations with several spatial extents. Simulated data were corrupted by physiological EEG noise. Simulations involving pairs of sources with the same amplitude were also studied. In addition to standard validation criteria (e.g., geodesic distance or mean square error), we proposed an original criterion dedicated to assess detection accuracy, based on receiver operating characteristic (ROC) analysis. Six source localization methods were evaluated: the minimum norm, the minimum norm weighted by multivariate source prelocalization (MSP), cortical LORETA with or without additional minimum norm regularization, and two derivations of the maximum entropy on the mean (MEM) approach. Results showed that LORETA-based and MEM-based methods were able to accurately recover sources of different spatial extents, with the exception of sources in temporo-mesial and fronto-mesial regions. Several spurious sources were generated by those methods, however, whereas methods using the MSP always located very accurately the maximum of activity but not its spatial extent. These findings suggest that one should always take into account the results from different localization methods when analyzing real interictal spikes.

Generalized epileptic discharges show thalamocortical activation and suspension of the default state of the brain.

Our objective was to evaluate the brain regions showing increased and decreased metabolism in patients at the time of generalized bursts of epileptic discharges in order to understand their mechanism of generation and effect on brain function. By recording the electroencephalogram during the functional MRI, changes in the blood oxygenation level-dependent signal were obtained in response to epileptic discharges observed in the electroencephalogram of 15 patients with idiopathic generalized epilepsy. A group analysis was performed to determine the regions of positive (activation) and negative (deactivation) blood oxygenation level-dependent responses that were common to the patients. Activations were found bilaterally and symmetrically in the thalamus, mesial midfrontal region, insulae, and midline and bilateral cerebellum and on the borders of the lateral ventricles. Deactivations were bilateral and symmetrical in the anterior frontal and parietal regions and in the posterior cingulate gyri and were seen in the left posterior temporal region. Activations in thalamus and midfrontal regions confirm known involvement of these regions in the generation or spread of generalized epileptic discharges. Involvement of the insulae in generalized discharges had not previously been described. Cerebellar activation is not believed to reflect the generation of discharges. Deactivations in frontal and parietal regions remarkably followed the pattern of the default state of brain function. Thalamocortical activation and suspension of the default state may combine to cause the actual state of reduced responsiveness observed in patients during spike-and-wave discharges. This brief lapse of responsiveness may therefore not result only from the epileptic discharge but also from its effect on normal brain function.

Detection of inter-hemispheric asymmetries of brain perfusion in SPECT.

Technetium-99m HMPAO and technetium-99m ECD single photon emission computed tomography (SPECT) imaging is commonly used to highlight brain regions with altered perfusion. It is particularly useful in the investigation of intractable partial epilepsy. However, SPECT suffers from poor spatial resolution that makes interpretation difficult. In this context, we propose an unsupervised voxel neighbourhood based method to assist the detection of significant functional inter-hemispheric asymmetries in brain SPECT, using anatomical information from MRI. For each MRI voxel, the anatomically homologous voxel in the contralateral hemisphere is identified. Both homologous voxel coordinates are then mapped into the SPECT volume using SPECT-MRI registration. Neighbourhoods are then defined around each SPECT voxel and compared to obtain a volume of inter-hemispheric differences. A volume including only the statistically significant inter-hemispheric differences is deduced from this volume using a non-parametric approach. The method was validated using realistic analytical simulated SPECT data including known asymmetries (in size and amplitude) as ground truth (gold standard). Detection performance was assessed using an ROC (receiver operating characteristic) approach based on the measures of the overlap between known and detected asymmetries. Validation with computer-simulated data demonstrates the ability to detect asymmetric zones with relatively small extension and amplitude. The registration of these detected functional asymmetries on the MRI enables good anatomical localization to be achieved.

A methodology for generating normal and pathological brain perfusion SPECT images for evaluation of MRI/SPECT fusion methods: application in epilepsy.

Quantitative evaluation of brain MRI/SPECT fusion methods for normal and in particular pathological datasets is difficult, due to the frequent lack of relevant ground truth. We propose a methodology to generate MRI and SPECT datasets dedicated to the evaluation of MRI/SPECT fusion methods and illustrate the method when dealing with ictal SPECT. The method consists in generating normal or pathological SPECT data perfectly aligned with a high-resolution 3D T1-weighted MRI using realistic Monte Carlo simulations that closely reproduce the response of a SPECT imaging system. Anatomical input data for the SPECT simulations are obtained from this 3D T1-weighted MRI, while functional input data result from an inter-individual analysis of anatomically standardized SPECT data. The method makes it possible to control the 'brain perfusion' function by proposing a theoretical model of brain perfusion from measurements performed on real SPECT images. Our method provides an absolute gold standard for assessing MRI/SPECT registration method accuracy since, by construction, the SPECT data are perfectly registered with the MRI data. The proposed methodology has been applied to create a theoretical model of normal brain perfusion and ictal brain perfusion characteristic of mesial temporal lobe epilepsy. To approach realistic and unbiased perfusion models, real SPECT data were corrected for uniform attenuation, scatter and partial volume effect. An anatomic standardization was used to account for anatomic variability between subjects. Realistic simulations of normal and ictal SPECT deduced from these perfusion models are presented. The comparison of real and simulated SPECT images showed relative differences in regional activity concentration of less than 20% in most anatomical structures, for both normal and ictal data, suggesting realistic models of perfusion distributions for evaluation purposes. Inter-hemispheric asymmetry coefficients measured on simulated data were found within the range of asymmetry coefficients measured on corresponding real data. The features of the proposed approach are compared with those of other methods previously described to obtain datasets appropriate for the assessment of fusion methods.

A data fusion environment for multimodal and multi-informational neuronavigation.

OBJECTIVE: Part of the planning and performance of neurosurgery consists of determining target areas, areas to be avoided, landmark areas, and trajectories, all of which are components of the surgical script. Nowadays, neurosurgeons have access to multimodal medical imaging to support the definition of the surgical script. The purpose of this paper is to present a software environment developed by the authors that allows full multimodal and multi-informational planning as well as neuronavigation for epilepsy and tumor surgery. MATERIALS AND METHODS: We have developed a data fusion environment dedicated to neuronavigation around the Surgical Microscope Neuronavigator system (Carl Zeiss, Oberkochen, Germany). This environment includes registration, segmentation, 3D visualization, and interaction-applied tools. It provides the neuronavigation system with the multimodal information involved in the definition of the surgical script: lesional areas, sulci, ventricles segmented from magnetic resonance imaging (MRI), vessels segmented from magnetic resonance angiography (MRA), functional areas from magneto-encephalography (MEG), and functional magnetic resonance imaging (fMRI) for somatosensory, motor, or language activation. These data are considered to be relevant for the performance of the surgical procedure. The definition of each entity results from the same procedure: registration to the anatomical MRI data set (defined as the reference data set), segmentation, fused 3D display, selection of the relevant entities for the surgical step, encoding in 3D surface-based representation, and storage of the 3D surfaces in a file recognized by the neuronavigation software (STP 3.4, Leibinger; Freiburg, Germany). RESULTS: Multimodal neuronavigation is illustrated with two clinical cases for which multimodal information was introduced into the neuronavigation system. Lesional areas were used to define and follow the surgical path, sulci and vessels helped identify the anatomical environment of the surgical field, and, finally, MEG and fMRI functional information helped determine the position of functional high-risk areas. CONCLUSION: In this short evaluation, the ability to access preoperative multi-functional and anatomical data within the neuronavigation system was a valuable support for the surgical procedure.