Search

There is an increasing demand for source analysis of neonatal EEG, but currently there is inadequate knowl- edge about i) the spatial patterning of neonatal scalp EEG and hence ii) the number of electrodes needed to capture neonatal EEG in full spatial detail. This study addresses these issues by using a very high density (2.5 mm interelectrode spacing) linear electrode array to assess the spatial power spectrum, by using a high density (64 electrodes) EEG cap to assess the spatial extent of the common oscillatory bouts in the neo- natal EEG and by using a neonatal size spherical head model to assess the effects of source depth and skull conductivities on the spatial frequency spectrum.The linear array recordings show that the spatial power spectrum decays rapidly until about 0.5–0.8 cycles per centimeter. The dense array EEG recordings show that the amplitude of oscillatory events decays within 4–6 cm to the level of global background activity, and that the higher frequencies (12–20 Hz) show the most rapid spatial decline in amplitude. Simulation with spherical head model showed that realistic variation in skull conductivity and source depths can both introduce orders of magnitude difference in the spatial fre- quency of the scalp EEG.Calculation of spatial Nyquist frequencies from the spatial power spectra suggests that an interelectrode dis- tance of about 6–10 mm would suffice to capture the full spatial texture of the raw EEG signal at the neonatal scalp without spatial aliasing or under-sampling. The spatial decay of oscillatory events suggests that a full representation of their spatial characteristics requires an interelectrode distance of 10–20 mm.The findings show that the conventional way of recording neonatal EEG with about 10 electrodes ignores most spatial EEG content, that increasing the electrode density is necessary to improve neonatal EEG source localization and information extraction, and that prospective source models will need to carefully consider the neonatally relevant ranges of tissue conductivities and source depths when source localizing cortical ac- tivity in neonates.

Electrical dipoles oriented perpendicular to the cortical surface are the primary source of the scalp EEGs and MEGs. Thus, in particular, gyri and sulci structures on the cortical surface have a definite possibility to influence the EEGs and MEGs. This was examined by comparing the spatial power spectral density (PSD) of the upper portion of the human cortex in MRI slices to that of simulated scalp EEGs and MEGs. The electrical activity was modeled with 2,650 dipolar sources oriented normal to the local cortical surface. The resulting scalp potentials were calculated with a finite element model of the head constructed from 51 segmented sagittal MR images. The PSD was computed after taking the fast Fourier transform of scalp potentials. The PSD of the cortical contour in each slice was also computed. The PSD was then averaged over all the slices. This was done for sagittal and coronal view both. The PSD of EEG and MEG showed two broad peaks, one from 0.05 to 0.22 cycles/cm (wavelength 20–4.545 cm) and the other from 0.22 to 1.2 cycles/cm (wavelength 4.545–0.834 cm). The PSD of the cortex showed a broad peak from 0.08 to 0.32 cycles/cm (wavelength 12.5–3.125 cm) and other two peaks within the range of 0.32 to 0.9 cycles/cm (wavelength 3.125–1.11 cm). These peaks are definitely due to the gyri structures and associated larger patterns on the cortical surface. Smaller peaks in the range of 1–3 cycles/cm were also observed which are possibly due to sulci structures. These results suggest that the spatial information was present in the EEG and MEG at the spatial frequencies of gyri. This also implies that the practical Nyquist frequency for sampling scalp EEGs should be 3.0 cycles/cm and an optimal interelectrode spacing of about 3 mm is needed for extraction of cortical patterns from scalp EEGs in humans.

Electrical dipoles oriented perpendicular to the cortical surface are the primary source of the scalp EEGs and MEGs. Thus, in particular, gyri and sulci structures on the cortical surface have a definite possibility to influence the EEGs and MEGs. This was examined by comparing the spatial power spectral density (PSD) of the upper portion of the human cortex in MRI slices to that of simulated scalp EEGs and MEGs. The electrical activity was modeled with 2,650 dipolar sources oriented normal to the local cortical surface. The resulting scalp potentials were calculated with a finite element model of the head constructed from 51 segmented sagittal MR images. The PSD was computed after taking the fast Fourier transform of scalp potentials. The PSD of the cortical contour in each slice was also computed. The PSD was then averaged over all the slices. This was done for sagittal and coronal view both. The PSD of EEG and MEG showed two broad peaks, one from 0.05 to 0.22 cycles/cm (wavelength 20–4.545 cm) and the other from 0.22 to 1.2 cycles/cm (wavelength 4.545–0.834 cm). The PSD of the cortex showed a broad peak from 0.08 to 0.32 cycles/cm (wavelength 12.5–3.125 cm) and other two peaks within the range of 0.32 to 0.9 cycles/cm (wavelength 3.125–1.11 cm). These peaks are definitely due to the gyri structures and associated larger patterns on the cortical surface. Smaller peaks in the range of 1–3 cycles/ cm were also observed which are possibly due to sulci structures. These results suggest that the spatial information was present in the EEG and MEG at the spatial frequencies of gyri. This also implies that the practical Nyquist frequency for sampling scalp EEGs should be 3.0 cycles/cm and an optimal interelectrode spacing of about 3 mm is needed for extraction of cortical patterns from scalp EEGs in humans.

Our objective was to examine if the high-density, 256 channel, scalp interictal EEG data can be used for localizing the epilepsy areas in patients. This was done by examining the long-range temporal correlations (LRTC) of EEGs and also that of the phase synchronization index (SI) of EEGs. It was found that the LRTC of scalp SI plots were better in localizing the seizure areas as compared with the LRTC of EEGs alone. The EEG data of one minute duration was filtered in the low Gamma band of 30-50 Hz. A detrended fluctuation analysis (DFA) was used to find LRTC of the scalp EEG data. Contour plots were constructed using a montage of the layout of 256 electrode positions. The SI was computed after taking Hilbert transform of the EEG data. The SI between a pair of channel was inferred from a statistical tendency to maintain a nearly constant phase difference over a given period of time even though the analytic phase of each channel may change markedly during that time frame. The SI for each electrode was averaged over with the nearby six electrodes. LRTC of the SI was computed and spatial plots were made. It was found that the LRTC of SI was highest at the location of the epileptic sites. A similar pattern was not found in the LRTC of EEGs. This provides a noninvasive way to localize seizure areas from scalp EEG data.

Phase slips arise from state transitions of the coordinated activity of cortical neurons which can be extracted from the EEG data. The phase slip rates (PSRs) were studied from the high-density (256 channel) EEG data, sampled at 16.384 kHz, of five adult subjects during covert visual object naming tasks. Artifact-free data from 29 trials were averaged for each subject. The analysis was performed to look for phase slips in the theta (4-7 Hz), alpha (7-12 Hz), beta (12-30 Hz), and low gamma (30-49 Hz) bands. The phase was calculated with the Hilbert transform, then unwrapped and detrended to look for phase slip rates in a 1.0 ms wide stepping window with a step size of 0.06 ms. The spatiotemporal plots of the PSRs were made by using a montage layout of 256 equidistant electrode positions. The spatiotemporal profiles of EEG and PSRs during the stimulus and the first second of the post-stimulus period were examined in detail to study the visual evoked potentials and different stages of visual object recognition in the visual, language, and memory areas. It was found that the activity areas of PSRs were different as compared with EEG activity areas during the stimulus and post-stimulus periods. Different stages of the insight moments during the covert object naming tasks were examined from PSRs and it was found to be about 512 ± 21 ms for the 'Eureka' moment. Overall, these results indicate that information about the cortical phase transitions can be derived from the measured EEG data and can be used in a complementary fashion to study the cognitive behavior of the brain. [ABSTRACT FROM AUTHOR]

Functional brain imaging and source localization based on the scalp's potential field require a solution to an ill-posed inverse problem with many solutions. This makes it necessary to incorporate a priori knowledge in order to select a particular solution. A computational challenge for some subject-specific head models is that many inverse algorithms require a comprehensive sampling of the candidate source space at the desired resolution. In this study, we present an algorithm that can accurately reconstruct details of localized source activity from a sparse sampling of the candidate source space. Forward computations are minimized through an adaptive procedure that increases source resolution as the spatial extent is reduced. With this algorithm, we were able to compute inverses using only 6% to 11% of the full resolution lead-field, with a localization accuracy that was not significantly different than an exhaustive search through a fully-sampled source space. The technique is, therefore, applicable for use with anatomically-realistic, subject-specific forward models for applications with spatially concentrated source activity. Index Terms--EEG, finite element method, FOCUSS, head model, inverse method, LORETA, source localization.

Functional brain imaging and source localization based on the scalp's potential field requires a solution to the inverse electrostatic problem. This is an underdetermined problem with many solutions. Minimum norm and regularization methods involving the norm are often used, but generally give solutions in which current is widely distributed. One method for reducing the spatial distribution of a solution is to apply it iteratively within the bounds of a shrinking ellipsoid. This paper compares the performance of this approach with an exhaustive search at various noise levels using a numeric simulation of the electroencephalogram in a realistic conductor model. The results show that inverting a single dipolar source with a location accuracy comparable to an exhaustive search requires in the range of 5 to 10 dB higher signal-to-noise ratio. Index Terms--EEG, finite element, head model, inverse method, minimum norm, regularization, source localization.

The current dipole is a widely used source model in forward and inverse electroencephalography and magnetoencephalography applications. Analytic solutions to the governing field equations have been developed for several approximations of the human head using ideal dipoles as the source model. Numeric approaches such as the finite-element and finite-difference methods have become popular because they allow the use of anatomically realistic head models and the increased computational power that they require has become readily available. Although numeric methods can represent more realistic domains, the sources in such models are an approximation of the ideal dipole. In this paper, we examine several methods for representing dipole sources in finite-element models and compare the resulting surface potentials and external magnetic field with those obtained from analytic solutions using ideal dipoles. Index Terms--Dipole, EEG, finite-element method, MEG, spherical head model.

Data from simulations of the anterior myocardial infarction (AMI) and inferior myocardial infarction (IMI) are presented. One infarct located in the anterior section of the left ventricle and a second one in the inferior wall of the left ventricle were modeled. A high-resolution finite element model of a heart and torso was used in this study. Differences in the normal and infarcted fields were computed. Our data suggest that the infarcted region contribution to the total magnetic field can be accounted for by an equivalent current dipole. It might also be possible to detect an infarct from these difference fields constructed for different cases of myocardial infarction. More simulations are needed to determine the relations between infarct sizes and locations and magnetic fields. These relations might then be used to detect various cases of myocardial infarction. Index Terms - Biomagnetic fields, biomagnetism, infarct simulations, magnetic fields of infarcts, MCG, torso models.

Magnetocardiograms (MCG's) simulated with a high-resolution heart-torso model of an adult subject were compared with measured MCG's acquired from the same individual. An exact match of the measured and simulated MCG's was not found due to the uncertainties in tissue conductivities and cardiac source positions. However, general features of the measured MCG's were reasonably represented by the simulated data for most, but not all of the channels. This suggests that the model accounts for the most important mechanisms underlying the genesis of MCG's and may be useful for cardiac magnetic field modeling under normal and diseased states. MCG's were simulated with a realistic finite-element heart-torso model constructed from segmented magnetic resonance images with 19 different tissue types identified. A finite-element model was developed from the segmented images. The model consists of 2.51 million brick-shaped elements and 2.58 million nodes, and has a voxel resolution of 1.56 x 1.56 x 3 mm. Current distributions inside the torso and the magnetic fields and MCG's at the gradiometer coil locations were computed. MCG's were measured with a Philips twin Dewar first-order gradiometer SQUID-system consisting of 31 channels in one tank and 19 channels in the other. Index Terms - Biomagnetic fields, biomagnetism, MCG, MCG simulations, torso fields, torso models.

Modeling in magnetoencephalography (MEG) and electroencephalography (EEG) requires knowledge of the in vivo tissue resistivities of the head. The aim of this paper is to examine the influence of tissue resistivity changes on the neuromagnetic field and the electric scalp potential. A high-resolution finite element method (FEM) model (452 162 elements, 2-mm resolution) of the human head with 13 different tissue types is employed for this purpose. Our main finding was that the magnetic fields are sensitive to changes in the tissue resistivity in the vicinity of the source. In comparison, the electric surface potentials are sensitive to changes in the tissue resistivity in the vicinity of the source and in the vicinity of the position of the electrodes. The magnitude (strength) of magnetic fields and electric surface potentials is strongly influenced by tissue resistivity changes, while the topography is not as strongly influenced. Therefore, an accurate modeling of magnetic field and electric potential strength requires accurate knowledge of tissue resistivities, while for source localization procedures this knowledge might not be a necessity. Index Terms - Biological tissues, conductivity, electroencephalography, finite element methods, magnetoencephalography.

Resolution of biomagnetic images using the technique of the alternating projections is proposed. Our image reconstruction procedure is divided in two steps. First, the biomagnetic inverse problem is solved by use of the projection theorem to reconstruct an initial image of the current distribution from a given magnetic field profile. Although the current distribution thus obtained has poor resolution, it can resemble the original shape of the current distribution. The second step improves the resolution of the reconstructed image by using the method of alternating projections. The procedure assumes that images can be represented by line like elements and involves finding the line like elements based on the initial image and projecting back onto the original solution space. Simulation studies were performed on a set of parallel conductors and a shape of the conductors in the form of letters, UWB @. All conductors were of line like thickness. Restored images closely resemble the original shape of the conductors.

The reconstruction of planar and three-dimensional current distributions from measured biomagnetic signals is a new field of research, known as biomagnetic computed tomography. This noninvasive imaging technique promises to provide precise, millimeter-sized resolution images of the electrical currents in tissues or organs. We performed simulation studies on phantom models of electrical sources. As a first step towards the development of an imaging algorithm, we addressed a simplified problem to identify the shape and direction of current flow in a planar surface. The problem was formulated by identifying a space in which the image was to be reconstructed. The space was segmented into a grid. Each grid space represented a current element. The magnetic field at a sampling point due to the current elements was computed using the Biot-Savart law. Since there were many more current elements than sample points, the problem was undetermined and had an uncountable number of solutions. The projection theorem was used to define an analytic solution for the magnitude and orientation of the current elements in the grid space. The solution required the inversion of large matrices in double precision. Such arrays were preprocessed on a mainframe computer, which permitted them to be rendered on any workstation. The accuracy of the image was determined by comparing it with the known location of the sources. Our results show that shape of the filamentary current flow can be imaged with our techniques. The resolution of images based on the sampling of the field, number of voxels in the reconstruction space, and noise is also analyzed.

Abstract Background The magnetoencephalograms (MEGs) are mainly due to the source currents. However, there is a significant contribution to MEGs from the volume currents. The structure of the anatomical surfaces, e.g., gray and white matter, could severely influence the flow of volume currents in a head model. This, in turn, will also influence the MEGs and the inverse source localizations. This was examined in detail with three different human head models. Methods Three finite element head models constructed from segmented MR images of an adult male subject were used for this study. These models were: (1) Model 1: full model with eleven tissues that included detailed structure of the scalp, hard and soft skull bone, CSF, gray and white matter and other prominent tissues, (2) the Model 2 was derived from the Model 1 in which the conductivity of gray matter was set equal to the white matter, i.e., a ten tissuetype model, (3) the Model 3 consisted of scalp, hard skull bone, CSF, gray and white matter, i.e., a five tissue-type model. The lead fields and MEGs due to dipolar sources in the motor cortex were computed for all three models. The dipolar sources were oriented normal to the cortical surface and had a dipole moment of 100 μA meter. The inverse source localizations were performed with an exhaustive search pattern in the motor cortex area. A set of 100 trial inverse runs was made covering the 3 cm cube motor cortex area in a random fashion. The Model 1 was used as a reference model. Results The reference model (Model 1), as expected, performed best in localizing the sources in the motor cortex area. The Model 3 performed the worst. The mean source localization errors (MLEs) of the Model 3 were larger than the Model 1 or 2. The contour plots of the magnetic fields on top of the head were also different for all three models. The magnetic fields due to source currents were larger in magnitude as compared to the magnetic fields of volume currents. Discussion These results indicate that the complexity of head models strongly influences the MEGs and the inverse source localizations. A more complex head model performs better in inverse source localizations as compared to a model with lesser tissue surfaces.

Abstract Background The structure of the anatomical surfaces, e.g., CSF and gray and white matter, could severely influence the flow of volume currents in a head model. This, in turn, will also influence the scalp potentials and the inverse source localizations. This was examined in detail with four different human head models. Methods Four finite element head models constructed from segmented MR images of an adult male subject were used for this study. These models were: (1) Model 1: full model with eleven tissues that included detailed structure of the scalp, hard and soft skull bone, CSF, gray and white matter and other prominent tissues, (2) the Model 2 was derived from the Model 1 in which the conductivity of gray matter was set equal to the white matter, i.e., a ten tissue-type model, (3) the Model 3 was derived from the Model 1 in which the conductivities of gray matter and CSF were set equal to the white matter, i.e., a nine tissue-type model, (4) the Model 4 consisted of scalp, hard skull bone, CSF, gray and white matter, i.e., a five tissue-type model. How model complexity influences the EEG source localizations was also studied with the above four finite element models of the head. The lead fields and scalp potentials due to dipolar sources in the motor cortex were computed for all four models. The inverse source localizations were performed with an exhaustive search pattern in the motor cortex area. The inverse analysis was performed by adding uncorrelated Gaussian noise to the scalp potentials to achieve a signal to noise ratio (SNR) of -10 to 30 dB. The Model 1 was used as a reference model. Results The reference model, as expected, performed the best. The Model 3, which did not have the CSF layer, performed the worst. The mean source localization errors (MLEs) of the Model 3 were larger than the Model 1 or 2. The scalp potentials were also most affected by the lack of CSF geometry in the Model 3. The MLEs for the Model 4 were also larger than the Model 1 and 2. The Model 4 and the Model 3 had similar MLEs in the SNR range of -10 dB to 0 dB. However, in the SNR range of 5 dB to 30 dB, the Model 4 has lower MLEs as compared with the Model 3. Discussion These results indicate that the complexity of head models strongly influences the scalp potentials and the inverse source localizations. A more complex head model performs better in inverse source localizations as compared to a model with lesser tissue surfaces. The CSF layer plays an important role in modifying the scalp potentials and also influences the inverse source localizations. In summary, for best results one needs to have highly heterogeneous models of the head for accurate simulations of scalp potentials and for inverse source localizations.

Abstract Background Our objective was to quantify and compare the extent of synchronization of the spatial-temporal myometrial activity over the human uterus before and during a contraction using transabdominal magnetomyographic (MMG) recordings. Synchronization can be an important indicator for the quantification of uterine contractions. Methods The spatialtermporal myometrial activity recordings were performed using a 151-channel noninvasive magnetic sensor system called SARA. This device covers the entire pregnant abdomen and records the magnetic field corresponding to the electrical activity generated in the uterine myometrium. The data was collected at 250 samples/sec and was resampled with 25 samples/sec and then filtered in the band of 0.1–0.2 Hz to study the primary magnetic activity of the uterus related to contractions. The synchronization between a channel pair was computed. It was inferred from a statistical tendency to maintain a nearly constant phase difference over a given period of time even though the analytic phase of each channel may change markedly during that time frame. The analytic phase was computed after taking Hilbert transform of the magnetic field data. The process was applied on the pairs of magnetic field traces (240 sec length) with a stepping window of 20 sec duration which is long enough to cover two cycle of the lowest frequency of interest (0.1 Hz). The analysis was repeated by stepping the window at 10 sec intervals. The spatial patterns of the synchronization indices covering the anterior transabdominal area were computed. For this, regional coil-pairs were used. For a given coil, the coil pairs were constructed with the surrounding six coils. The synchronization indices were computed for each coil pair, averaged over the 21 coil-pairs and then assigned as the synchronization index to that particular coil. This procedure was tested on six pregnant subjects at the gestational age between 29 and 40 weeks admitted to the hospital for contractions. The RMS magnetic field for each coil was also computed. Results The results show that the spatial patterns of the synchronization indices change and follow the periodic pattern of the uterine contraction cycle. Spatial patterns of synchronization indices and the RMS magnetic fields show similarities in few window frames and also show large differences in few other windows. For six subjects, the average synchronization indices were: 0.346 ± 0.068 for the quiescent baseline period and 0.545 ± 0.022 at the peak of the contraction. Discussion These results show that synchronization indices and their spatial distributions depict uterine contractions and relaxations.