|
Size: 5582
Comment:
|
Size: 7931
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| = Group studies: Subject co-registration = | = Group studies: Subject coregistration = |
| Line 14: | Line 14: |
| 1. 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. 1. Do the same for the low-resolution cortex surface (Subject01/cortex_300000V). You see the downsampled version of the same thing. 1. Now do the same for the high-resolution cortex of the default anatomy. |
1. 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. 1. 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.<<BR>><<BR>> {{attachment:sphereSubj.gif||height="173",width="538"}} 1. 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. <<BR>><<BR>> {{attachment:sphereDefault.gif||height="173",width="538"}} 1. 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).<<BR>><<BR>> {{attachment:sourcesSubj.gif||height="209",width="561"}} 1. 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) <<BR>><<BR>> {{attachment:projectRaw.gif||height="226",width="394"}} 1. It creates a new subject "Group analysis" to store the projected source maps:<<BR>><<BR>> {{attachment:projectTreeRaw.gif||height="203",width="386"}} 1. Double-click on this new file, and display it on the registered spheres as well:<<BR>><<BR>> {{attachment:sourcesDefault.gif||height="209",width="561"}} 1. 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.<<BR>><<BR>>{{attachment:projectTreeHires.gif}}<<BR>>{{attachment:sourcesHires.gif||height="209",width="561"}} 1. Note that 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. |
| Line 18: | Line 24: |
| The best available option is selected automatically == Those menus don't show up == |
== It doesn't work == |
| Line 24: | Line 28: |
| The best available option is selected automatically |
|
| Line 32: | Line 38: |
| 1. Interpolation of the source amplitudes from the subject surface to the Colin27 surface using Shepard’s method (weighted combination of a few nearest neighbors) | 1. 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) |
| Line 42: | Line 48: |
| 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". | 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". |
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 each hemisphere. The hemisphere of the subject's brain are inflated to a sphere, then this sphere is 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.
We can illustrate this method using the protocol TutorialRaw, that you created while reading the introduction tutorials.
Select the protocol TutorialRaw, go the anatomy view
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.
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.
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.
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).
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)
It creates a new subject "Group analysis" to store the projected source maps:
Double-click on this new file, and display it on the registered spheres as well:
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.
- Note that 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.
It doesn't work
If you have started your analysis before those options were made available in Brainstorm
A cheap alternative
The best available option is selected automatically
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.
Right-click on Subject01/Right/MN > Project sources > Default anatomy > Cortex_15000V
- 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 intensly 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:
- 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)
The interpolation matrix subject->Colin is computed only once and then saved in the subject's cortex surface file, in the field "tess2tess_interp"
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:
Display at the same time the original source file and the projected one to see how it looks like (left: Subject01, right: Colin)
- 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".
Feedback