Tutorial 20: Head modelling

Authors: Francois Tadel, Elizabeth Bock, John C Mosher, Richard Leahy, Sylvain Baillet

The following tutorials describe how brain activity can be estimated from the MEG/EEG recordings we have processed so far. This step consists in solving two separate modeling problems: the modeling of the electromagnetic properties of the head and of the sensor array (a.k.a. head model or forward model) and the estimation of the brain sources which have produced the data. That second step is known as source modeling or solving an inverse problem. It requires that the forward modeling of the head tissues and sensor instrumentation is completed first. This tutorial explains how to compute a head model for the subject of the auditory oddball experiment. As far as source modeling is concerned, we focus the tutorial on linear estimates of distributed source models, which are popular and physiologically plausible approaches (there is no dipole fitting available in the software): you may want to refer to other sources for a complete description of source imaging with MEG and EEG.

Advanced

Why estimating sources?

Reconstructing the activity of the brain from MEG or EEG recordings involves several sophisticated steps. Athough Brainstorm simplifies the procedures required, it is important to understand whether source modeling is essential to answer the neuroscience question which brought you to collect data in the first place.

If one of your primary objectives is to identify and map the regions of the brain involved in a specific stimulus response or behaviour, source estimation can help address this aspect. Empirical interpretations of sensor topographies can inform where about brain generators might be located: which hemisphere, what broad aspect of the anatomy (e.g., frontal vs. posterior regions). Source estimation methods improves anatomical resolution with respect to the interpreration of sensor patterns. Spatial resolution in MEG and EEG depends on source depth and orientation and overall SNR: still, one can expect to be able to source activations within the millimeter to centimeter range, especially relatively, when contrasts between conditions are implemented in the study design.

Source mapping is a form of spatial deconvolution of sensor data. In EEG in particular, scalp topographies are very smooth and it is common to that contributions from distant brain regions overalp over large clusters of electrodes. Moving to the source space can help separating the active regions.

More specifically in MEG, source maps can be a great assest to alleviate some specific issues with that modality. In MEG, the head of the subject is not fixed. Hence sensor topographies depend on the actual position of the subject under the rigid helmet. Therefore, between two acquisition runs, or between subjects with different head shapes and sizes, the same MEG sensors may pick up signals from different parts of the brain. This problem does not exist in EEG, where electrode montages are attached to the head and arranged according to standard positions.

Another important point to consider when interpreting MEG sensor maps is that MEG manufacturers use different types of sensor technology (e.g., magnetometers vs. gradiometers; axial vs. tangential gradiometers, etc. yielding different physical measures). This is not an issue with EEG, with essentially one sensor type (electrodes, dry or active, all measuring Volts): two EEG electrode caps with same number of electrodes only differ in how well they ensure good contact with the scalp. Working in source space alleviates all these aspects.

Nevertheless, if your neuroscience question can be solved by measuring signal latencies over broad regions, or other aspects which do not depend crucially on anatomical localization (such as global signal properties integrated over all or clusters of sensors), source modeling might not be a requirement. To sort out this question wil influence the time and hardware requirements to complete your data analysis (source analysis multiplies the needs in terms of disk storage and computational specifications).

Advanced

The origins of MEG/EEG signals

To better understand how the forward model is elaborated, we need to have at least a basic understabnding of the physiological origins of MEG/EEG signals. Note that, as always when dealing with modeling, we need to deal with various degrees of approximation.

Overall, it is assumed that most of the MEG/EEG signals are generated by postsynaptic activity of ensembles of cortical pyramidal neurons of the cerebral cortex. The reason is essentially in the morphology and mass effect of these cells, which present elongated shapes and are oriented perpendicularly to the cortical surface. Mass effects of close-to-simulatenous changes in post-synaptic potentials pyramidal neural assemblies add up in time and space. These effects can conveniently be modeled at a mesoscopic spatial scale with electric dipoles distributed along the cortical mantle (green arrows in figure below). Note that there is growing evidence that MEG and EEG are also sensitive to deeper, cortical and subcortical structures, including brain nuclei and the cerebellum, where pyramidal cells are rare or absent. Brainstorm features advanced models that can include as an option, these structures. The emphasis in the present tutorial is on pyramidal cell assemblies for simplicity.

The primary and volume currents generated by these dipoles create differences in electrical potentials and magnetic fields that can be detected outside the head. They can be measured with electrodes placed on the skin (EEG, with respect to a reference) or very sensitive magnetic detectors (MEG).

Advanced

Source space

Dipole fitting vs distributed models

MEG/EEG source estimation consists in modeling brain activity with current dipoles. A current dipole is a convenient model equivalent to the net electrophysiological activity of local assemblies of neurons. Two main approaches have been expplored for source MEG/EEG estimation: dipole fitting methods - where the position and amplitude of a one to a few equivalent current dipoles (ECD) are estimated over relativelyshort time windows - and distributed models - where the location (and typically, the orientation) of a large number dipoles is fixed. The dipoles sample a spatial grid covering the entire brain volume or the cortical surface - only their amplitude is estimated at each time point.

Equivalent dipole fitting approaches are quite straightforward and can be adequate when the number of brain regions expected to be active is small. Therefore, it is most adequate for responses at early latencies. They cannot generalize to capture complex dynamics over extended period of time (epochs) and the associated estimation techniques are quite sensitive to intial conditions (how many dipoles to fit? Where does the serch start? etc.). Our strategy in Brainstorm is to promote distributed source models, which are less user-sensitive, can genralize to all experimental conditions, and yield time-resolved image volumes that can be processed in many different, powerful ways (group statistics, spatial segmentation, use of regions fo interest, correspondence with fMRI, etc.)

Location constraints

When opting for distributed source models, the positions and orientations of the elementary dipoles that will define the "voxel" grid of the source images produced need to be defined. This set of dipoles is called the source space. By default, Brainstorm constrains the source space to the cortex, where signal-to-noise and sensitivity is maximum in MEG/EEG. Note however that more complete models that include subcortical structures and the cerebellum are available in Brainstorm. Therefore, one decision you need to make before proceeding with source imaging is whether more complete source spaces are required to answer your neuroscience question.

For this tutorial, we use the simple approach where a current dipole is automatically assigned to each of the vertices of the cortical surface we have imported so far (see the nodes in the grey mesh in the leftmost image below). Also in this turorial, we elected not to contrain the dipole orientations to be perpendicular to the cortical surface. This choice can be motivated by the fact that you have no access to the individual MRI of your participant and thefore resort to a generic template, or that you find the cortical orientation constraint to be poorly unjustified.

In this context, three orthogonal dipoles are assigned to each vertex of the cortex. This triplet forms a basis to account for local currents flowing in any arbitrary direction, not just orthogonal to the cortex.

When importing the anatomy of the subject, we downsampled the cortex surface to 15,000 vertices. This will correspond to a source space of 45,000 dipoles. We will compute a forward model that connects the activity of these 45,000 dipoles with the 275 MEG sensors we have in this dataset.

This default number of 15,000 vertices is empirical. Over the years, our experience seemed to show that it represents a good balance between the representation of the brain circumvolutions, the surface sampling and the amount of data that is generated. Using less vertices makes it difficult to preserve the shape of the brain, using more vertices produces more data without adding to the spatial resolution of the method and may lead to computational memory issues.

Orientation constraints

Additionally, we can impose constraints of orientation on the dipoles, to match the physiological observation that the pyramidal cells are mostly organized perpendicularly to the cortex surface. This has the advantage of limiting the number of dipoles to 15,000 (one per vertex) and making the results much easier to display and process. However, this constraint is most of the time too strong and distorts the reconstruction. This orientation constraint is an option of the inverse model and will be discussed in the following introduction tutorials.

Fully unconstrained

The spatial constraint of imposing all the dipoles to be on the cortical surface might also be too restrictive in some cases, because our model is then not able to correctly represent the activity in deeper brain structures or in the cerebellum. Therefore we also offer an option to use the entire brain volume as our source space (the green dots below represent dipoles locations in volume model). This produces results that can be better or worse depending on the data, but in all the cases much more difficult to review. Volume and mixed head volumes are discussed in the advanced tutorials about source modeling.

Advanced

Forward problem

The first step of the source reconstruction consists in computing a model that explains how the electric currents or the magnetic fields flow from the electric generators in the brain (source space) through the different tissues of the head (mostly brain, CSF, skull and skin), to finally reach the sensors.

Available methods for MEG forward modeling

Computation

The forward models are related with the anatomy of the subject and the description of the sensors, therefore the menus associated to its computation are attached to the channel file.

Repeat the same operation for the other file. We have two different acquisition runs with two different relative position of the head and the sensors, therefore we need to compute two different head models.

Advanced

Database explorer

Additional considerations about the management of the head model files.

Advanced

On the hard drive

Right-click on any head model > File > View file contents:

Structure of the head model files

Gain matrix

References

Additional documentation








Feedback: Comments, bug reports, suggestions, questions
Email address (if you expect an answer):


Tutorials/HeadModel (last edited 2016-01-25 22:37:24 by FrancoisTadel)