[FieldTrip] Using an fMRI cluster as an ROI

Julian Keil julian.keil at gmail.com
Mon Jul 25 14:03:26 CEST 2016


Hi Ross,

one way to attempt this, is to create a 3D grid on your fMRI data, so that each voxel is a grid point.
Then you can compute the distance between grid points from the fMRI grid and your virtual electrode grid and find the nearest neighbours.
This sould work for the whole brain as well as for predefined ROIs.

I've attached a function I use for similar problems. I'm sure there are more elegant ways but this works for me. It won't do exactly what you want to do, but maybe some bits and pieces can help you along.

Good luck,

Julian




Am 25.07.2016 um 13:47 schrieb Ross Wilson:

> Hello all!
> 
> I am attempting to use an fMRI cluster as an ROI to look at virtual electrode timecourses post-beamforming.
> 
> I have been using the aal atlas to find virtual electrodes of maximum change between baseline and stimulus period within, for example, the motor cortex afer beamforming. However, the virtual electrodes found are not necessarily in a position similar to that of the fMRI cluster and I’d like to be more specific than that if possible.
> 
> The template grid I have been using for beamforming is created from the T1 (code as below, in MNI space: …fieldtrip\external/spm8/templates/T1.nii). Is it possible to transform an fMRI cluster (in MNI space) onto the template grid so that the data from specific grid points relating to that cluster can be extracted post-beamforming??
> 
> Any help with this would be great!
> 
> Thank you,
> 
> Ross
> %%%tempale grid creation
> T1 = ft_read_mri('F:\fieldtrip\external/spm8/templates/T1.nii');
> T1.coordsys = 'spm';
> 
> cfg          = [];
> template_seg = ft_volumesegment(cfg, T1);
> 
> cfg          = [];
> cfg.method   = 'singleshell';
> template_vol = ft_prepare_headmodel(cfg, template_seg);
> template_vol = ft_convert_units(template_vol, 'cm');
> 
> cfg = [];
> cfg.grid.xgrid  = -20:0.5:20;
> cfg.grid.ygrid  = -20:0.5:20;
> cfg.grid.zgrid  = -20:0.5:20;
> cfg.grid.unit   = 'cm';
> cfg.grid.tight  = 'yes';
> cfg.inwardshift = 0.15;
> cfg.vol        = template_vol;
> template_grid  = ft_prepare_sourcemodel(cfg);
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> https://mailman.science.ru.nl/mailman/listinfo/fieldtrip

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20160725/1c223b0a/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: vt_make_roifield.m
Type: application/octet-stream
Size: 6646 bytes
Desc: not available
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20160725/1c223b0a/attachment.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20160725/1c223b0a/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 495 bytes
Desc: Message signed with OpenPGP using GPGMail
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20160725/1c223b0a/attachment.sig>


More information about the fieldtrip mailing list