16
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Zeffiro user interface for electromagnetic brain imaging: a GPU accelerated FEM tool for forward and inverse computations in Matlab

      Preprint
      , ,

      Read this article at

      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          This article introduces the Zeffiro interface (ZI) for brain imaging. ZI is an attempt to provide a simple and accessible and multimodal open source platform for FEM based forward and inverse computations. Zeffiro is Italian for a 'gentle breeze'. ZI has been implemented for the Matlab environment. It allows one to (1) generate the FE mesh for a given multi-layer head model, (2) form a basis for computing a LF matrix and inverting a given set of measurements, and (3) to visualize the reconstructions on the FE mesh. To tackle the issue of the high computational cost, ZI uses graphics processing unit (GPU) acceleration for each of the processing stages (1)-(3). In its current configuration, ZI includes forward solvers for E/MEG and linearized EIT LF matrices and and a set of inverse solvers utilizing the hierarchical Bayesian model (HBM). In this paper, we describe the mathematics behind the ZI, describe the principle operations and usage of ZI, concentrating especially on EEG/MEG, and introduce some central points for the developer perspective. We also give examples of computing times, and analyze numerically how the inversion parameters should be selected, when inverting brain activity with HBM.

          Related collections

          Most cited references17

          • Record: found
          • Abstract: not found
          • Article: not found

          Influence of transcortical d-c currents on cortical neuronal activity.

            Bookmark
            • Record: found
            • Abstract: found
            • Article: not found

            Electrode models for electric current computed tomography.

            This paper develops a mathematical model for the physical properties of electrodes suitable for use in electric current computed tomography (ECCT). The model includes the effects of discretization, shunt, and contact impedance. The complete model was validated by experiment. Bath resistivities of 284.0, 139.7, 62.3, 29.5 omega.cm were studied. Values of "effective" contact impedance zeta used in the numerical approximations were 58.0, 35.0, 15.0, and 7.5 omega.cm2, respectively. Agreement between the calculated and experimentally measured values was excellent throughout the range of bath conductivities studied. It is desirable in electrical impedance imaging systems to model the observed voltages to the same precision as they are measured in order to be able to make the highest resolution reconstructions of the internal conductivity that the measurement precision allows. The complete electrode model, which includes the effects of discretization of the current pattern, the shunt effect due to the highly conductive electrode material, and the effect of an "effective" contact impedance, allows calculation of the voltages due to any current pattern applied to a homogeneous resistivity field.
              Bookmark
              • Record: found
              • Abstract: found
              • Article: not found

              A guideline for head volume conductor modeling in EEG and MEG.

              For accurate EEG/MEG source analysis it is necessary to model the head volume conductor as realistic as possible. This includes the distinction of the different conductive compartments in the human head. In this study, we investigated the influence of modeling/not modeling the conductive compartments skull spongiosa, skull compacta, cerebrospinal fluid (CSF), gray matter, and white matter and of the inclusion of white matter anisotropy on the EEG/MEG forward solution. Therefore, we created a highly realistic 6-compartment head model with white matter anisotropy and used a state-of-the-art finite element approach. Starting from a 3-compartment scenario (skin, skull, and brain), we subsequently refined our head model by distinguishing one further of the above-mentioned compartments. For each of the generated five head models, we measured the effect on the signal topography and signal magnitude both in relation to a highly resolved reference model and to the model generated in the previous refinement step. We evaluated the results of these simulations using a variety of visualization methods, allowing us to gain a general overview of effect strength, of the most important source parameters triggering these effects, and of the most affected brain regions. Thereby, starting from the 3-compartment approach, we identified the most important additional refinement steps in head volume conductor modeling. We were able to show that the inclusion of the highly conductive CSF compartment, whose conductivity value is well known, has the strongest influence on both signal topography and magnitude in both modalities. We found the effect of gray/white matter distinction to be nearly as big as that of the CSF inclusion, and for both of these steps we identified a clear pattern in the spatial distribution of effects. In comparison to these two steps, the introduction of white matter anisotropy led to a clearly weaker, but still strong, effect. Finally, the distinction between skull spongiosa and compacta caused the weakest effects in both modalities when using an optimized conductivity value for the homogenized compartment. We conclude that it is highly recommendable to include the CSF and distinguish between gray and white matter in head volume conductor modeling. Especially for the MEG, the modeling of skull spongiosa and compacta might be neglected due to the weak effects; the simplification of not modeling white matter anisotropy is admissible considering the complexity and current limitations of the underlying modeling approach.
                Bookmark

                Author and article information

                Journal
                19 November 2018
                Article
                1811.07717
                4db8cc0b-5851-4fe5-b913-1b2b7d8e430a

                http://arxiv.org/licenses/nonexclusive-distrib/1.0/

                History
                Custom metadata
                65L60, 65N30, 34A55, 34M50, 92C50
                cs.MS cs.CE math.AP

                Analysis,Mathematical software,Applied computer science
                Analysis, Mathematical software, Applied computer science

                Comments

                Comment on this article