Group analysis: Subject coregistration

Authors: Francois Tadel

In the typical Brainstorm workflow, you process each subject individually: you have access to the MRI of the subject, you extract the cortex envelope from it using the software of your choice (FreeSurfer, BrainVISA, BrainSuite or CIVET), you import the recordings and estimate the sources for each subject separately.

After importing and pre-processing all the individual information, you will probably be interested in calculating some group statistics, to extract the significant features across all the subjects in your study. in order to do this at the source level, all the source files must be using the same support cortex source (= the same source space). The problem you typically face is that the shape and number of vertices of the cortex surface varies across subject: one source index doesn't mean the same thing for all the subjects, hence it's impossible to compare them directly.

The solution is to reproject the sources of all the subjects on the default anatomy. The reprojected versions share the same source space, the cortex from ICBM152 anatomy, so they can be averaged or used to calculate any other statistics across subjects. Different methods are available in Brainstorm depending on the software you are using for the MRI segmentation. The recommended option is to use FreeSurfer and benefit from the very efficient registration system this software offers.

Using FreeSurfer registered spheres

The FreeSurfer software registers all the subjects to the FSAverage atlas using a spherical representation of the cortex. Each hemisphere of the subject's brain is inflated to a sphere, which is then deformed to match the curvature map of the equivalent sphere in the FSAverage subject. The result of this registration is saved in the FreeSurfer output folder in surf/lh.sphere.reg and surf/rh.sphere.reg. When you import the FreeSurfer anatomy folder in Brainstorm, it reads these files and save the position of the vertices in the field ".Reg.Sphere.Vertices" of the cortex, white, and mid-gray surfaces.

We illustrate this method using the dataset from the tutorial MEG median nerve (CTF), but it is very easy to transpose to the protocol TutorialIntroduction you created in the introduction tutorials.

  1. Select the protocol TutorialRaw, go the anatomy view

  2. Right-click on the high-resolution cortex surface (Subject01/cortex_300000V) > MRI registration > View FreeSurfer sphere. These two spheres represent the inflated left and the right hemispheres, registered on the FSAverage subject. The darker stripes represent the bottom of the sulci: the map is exactly the same as the one that is displayed on the cortex surface when you click on the "Sulci" button in the Surface tab.

  3. Do the same for the low-resolution cortex surface (Subject01/cortex_300000V). This is the downsampled version of the same thing. You can click on the "Edge" button on the surface tab to see the same triangulation as you would see when displaying the low-resolution cortex surface.


  4. Now do the same for the cortex of the default anatomy. All these spheres are registered to the FSaverage subject (the vertices of the spheres are moved to match the curvature maps of the FSAverage sphere). Hence they are also registered together: we can use them to reproject the source calculated on the individual anatomies to the default anatomy.


  5. Now switch to the functional part of the protocol. Right-click on the average sources for condition Left > Cortical activations > Display on spheres. Double-click on the file to display the regular brain at the same time. Change the current time to a meaningful event (35ms for instance).


  6. To project this source file to the default anatomy ICBM152, right-click on the Left/Average/MN source file > Project sources > Default anatomy > cortex_15000V. The source values are re-interpolated from the registered spheres of Subject01 to the spheres of ICBM152 using the Shepard’s method (weighted combination of a few nearest neighbors)


  7. It creates a new subject "Group analysis" to store the projected source maps:


  8. Double-click on this new file, and display it on the registered spheres as well:


Projection on high-resolution cortex

You can also use the same method to re-project the sources on the high-resolution cortex, to produce nicer figures. The process is the same: right-click on Left/Average/MN source file > Project sources > Subject01 > cortex_300000V. It creates a new file in the same condition.

Important notes:

It doesn't work

If you have started your analysis before these options were made available in Brainstorm you would probably not be able to apply these operations to your own data.

Registration with BrainSuite

BrainSuite also provides an accurate registration method to the Brainstorm anatomy templates (Colin27, ICBM152). The registration procedure works in the same way as the FreeSurfer-based registration illustrated in this tutorial.

To use this registration, you need to use a template anatomy processed with BrainSuite.
Right-click on Default anatomy > Use template > Colin27_BrainSuite or ICBM152_BrainSuite.

A cheap alternative

If you did not process your MRI files with FreeSurfer, or if you are working with an older database, you have another option available to re-project your subject's source maps on the default anatomy. The results are much less reliable, but the method is quite robust and can be used for all sorts of purposes. We can illustrate this method using the protocol TutorialCTF, that you created while reading the introduction tutorials.

Running this method is done exactly in the same way as with the FreeSurfer spheres, using the "Project sources" menu. There is no additional option, Brainstorm would just automatically choose the best method available for the source projection. If the registered spheres are available in the field ".Reg.Sphere.Vertices", it uses them, if not it uses this method.

  1. Right-click on Subject01/Right/MN > Project sources > Default anatomy > Cortex_15000V


  2. The source cortex (Subject01) and the destination cortex (ICBM152) are split in two hemispheres, and then for each hemisphere:
    • Alignment along the anterior commissure/posterior commissure axis
    • Smooth intensely to keep only the main featuresof the surfaces
    • Deformation of the individual surface to match the ICBM152 surface with an iterative closest point algorithm (ICP). The following figures show the left hemispheres before and after the ICP deformation:

      projectBefore.gif projectAfter.gif

  3. Interpolation of the source amplitudes from the subject surface to the ICBM152 surface using the Shepard’s method (weighted combination of a few nearest neighbors)
  4. The interpolation matrix subject->Colin is computed only once and then saved in the subject's cortex surface file, in the field "tess2tess_interp"

  5. You can see that a new subject and a new condition were created in the database, that contain the sources from Subject01/Right/ERF projected on the default anatomy:


  6. Display at the same time the original source file and the projected one to see how it looks like (left: Subject01, right: Colin)


  7. Delete the "Group analysis" subject.

Warped brains

It is important to note that pseudo-individual anatomies that were generated by deforming the ICBM152 brain using the digitized head points do not need to be re-projected on the default anatomy.

If you have two subjects for which you generated pseudo anatomies using the popup menu Digitized head points > Warp > Deform default anatomy to fit these points, you can directly use their source maps for computing any statistics or difference. You don't need to use this menu "Project sources".


Volume source models

Source files using volume head models contain a list of grid points that are not necessarily positioned on the cortex surface. We cannot use the surface-based atlases from FreeSurfer or BrainSuite. Instead, we rely on the MNI transformation computed using SPM functions.

The procedure to project these source files on a template is different: You need to define a grid at the template level, then project this grid in the subject space and estimate the sources for this grid.

Additional discussions on the forum

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

Tutorials/CoregisterSubjects (last edited 2016-08-02 17:21:32 by FrancoisTadel)