Group studies: Subject coregistration

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 Colin27 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 those files and save the position of the vertices in the field ".Reg.Sphere.Vertices" of the cortex, white, and mid-gray surfaces.

We can illustrate this method using the protocol TutorialRaw, that you created while reading 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. Those 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 those 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 Colin27, 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 Colin27 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:

• You probably don't need 300,000 vertices on the cortex to get figures nice enough for your publications. You can downsample the high-resolution to something around 90,000 vertices, and then re-project the sources on this mid-resolution brain.

• To get a smoother display instead of the all those little dots due to the Shepards interpolation, you can run the process Sources > Spatial smoothing on the high-resolution cortical maps.

It doesn't work

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

• You don't see the menus described above: update Brainstorm.
• You get the error "There is no registered sphere available for this surface" for the subject cortex:

  • Re-import the anatomy of the subject from the FreeSurfer folder after making sure that the files surf/?h.sphere.reg are not missing. Then delete all the head models and sources results, and re-calculate them.

  • Warning: Don't do this if you already have a lot of results for this subject, it would probably make all the source results unusable.
• You get this error for the default anatomy:

  • You have to update the default anatomy for this protocol. Right-click on the folder (Default anatomy) > Use template > Colin27

  • Warning: Don't do this if you already have projected all your results on the default anatomy, or if you have some subjects that use the default anatomy for the source estimation. It would make all the results unusable for sure.
• When you use the menu "Project sources", you see that it's using the second method described below (you see figures with those smoothed red and gray surfaces): it means that either the subject or the default anatomy doesn't have the registered spheres saved in its structure. You should also observe one or the other of the above error, fix accordingly.

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 chose automatically 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 (Colin27) 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 Colin27 surface with an iterative closest point algorithm (ICP). The following figures show the left hemispheres before and after the ICP deformation:
     
     

3. Interpolation of the source amplitudes from the subject surface to the Colin27 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 Colin27 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".

Additional discussions on the forum

• Spatial smoothing of sources: http://neuroimage.usc.edu/forums/showthread.php?1409

Feedback