Time Frequency (complex) and projection on default anatomy

Francois · October 29, 2019, 9:30am

The scouts are defined on the default anatomy and comes from FreeSurfer/HCP. The default anatomy I'm using is a fsaverage on which I imported the HCP-MM atlas from Glasser et al., 2016 (A multi-modal parcellation of human cerebral cortex | Nature). Thus everything (registered spheres etc..) are available. However, I don't have the atlas in the individual surfaces.

Everything should work then!
Project your scouts of interest from the template to the individual subjects, and do a ROI-based analysis. Then you don't need to any additional interpolation for the group analysis.

I computed it using the same approach as in the Matti's MNE manual (morphing) and used it to switch from individual to group space by matrix multiplication. Is it computed differently in Brainstorm?

This is described in this tutorial:
https://neuroimage.usc.edu/brainstorm/Tutorials/CoregisterSubjects#Using_FreeSurfer_registered_spheres

If you need additional details, refer to the code of tess_interp_tess2tess.m:

github.com

brainstorm-tools/brainstorm3/blob/master/toolbox/anatomy/tess_interp_tess2tess.m#L318


      
          iScoutDest = find(strcmpi(sScoutSrc.Label, {sScoutStructDest.Label}));

          if isempty(iScoutDest)

              disp(['PROJECT> Warning: Structure not found in destination surface: ' sScoutSrc.Label ]);

              continue;

          end

          sScoutDest = sScoutStructDest(iScoutDest);

          % Is it a cortex region

          isCortexL = ismember(sScoutSrc.Label, {'lh', '01_Lhemi L', 'Cortex L'});

          isCortexR = ismember(sScoutSrc.Label, {'rh', '01_Rhemi R', 'Cortex R'});

          

          % ===== USE FREESURFER SPHERES =====

          % Interpolate using the sphere and the Shepard's algorithm

          if isCortexL && isFreeSurfer

              Wmat(sScoutDest.Vertices, sScoutSrc.Vertices) = bst_shepards(vertSphLdest, vertSphLsrc, nbNeighbors, 0);

          elseif isCortexR && isFreeSurfer

              Wmat(sScoutDest.Vertices, sScoutSrc.Vertices) = bst_shepards(vertSphRdest, vertSphRsrc, nbNeighbors, 0);

              

          % ===== USE BRAINSUITE SQUARES =====

          % Interpolate using the Brainsuite squares and the Shepard's algorithm

          elseif isCortexL && isBrainSuite

              % Interpolation: Subject => BrainSuiteAtlas1