[FieldTrip] Adding MRI Overlay to ft_sourceplot Images

Julian Keil julian.keil at gmail.com
Mon Jul 9 08:01:33 CEST 2018


Hi Lindsey,

What is "grid" bases on? Is this the individual mri or a standard mri?
In any case, you can use ft_sourceinterpolate to interpolate to a standard
MRI.
Then, you can set the anatomy parameter in ft_spurceplot. Check the details
of both functions for the relevant cfg settings.

Hope this helps,

Cheers,
Julian

Tate, Lindsey R. <lindseyrtate at ou.edu> schrieb am Mo. 9. Juli 2018 um 07:35:

> Hello fieldtrip community,
>
>
> I have a project with a lot of source analysis of EEG data. I'm using
> ft_sourceanalysis and ft_sourceplot. I would like to add standard (e.g.,
> Colin 27) MRI slices to these images as an overlay (underlay?) for
> reference. Is there a way to do this? How do I align the MRI with my source
> space?
>
>
> Example code for ft_sourceanalysis and ft_sourceplot are provided below my
> signature in case you need the reference.
>
>
> Thank you,
>
>
> Lindsey Tate
> PhD Candidate
> University of Oklahoma
> lindseyrtate at ou.edu
>
>
> ///
>
> %baseline period source analysis
>
> cfg = [];
> cfg.elec = elec ;
> cfg.method = 'dics';
> cfg.frequency = 25 ;
> cfg.grid = grid;
> cfg.headmodel = headmodel;
> cfg.dics.projectnoise = 'yes';
> cfg.dics.lambda = 0;
> cfg.latency = -0.7490 ;
> [presource25] = ft_sourceanalysis(cfg, prefreq25) %from ft_freqanalysis
>
> %same is done for a post-stimulus period, called pmsource25; code omitted
> for brevity
>
> %difference between pre- and post-stimulus periods
> source25diff = pmsource25 ;
> source25diff.avg.pow = (pmsource25.avg.pow - presource25.avg.pow) ./
> presource25.avg.pow;
>
> %interpolate the difference between pre- and post-stimulus periods
> cfg = [];
> cfg.downsample = 2;
> cfg.parameter = 'avg.pow';
> source25DiffInt = ft_sourceinterpolate(cfg, source25diff, grid);
>
> %plot
> cfg = [];
> cfg.method = 'slice';
> cfg.funparameter = 'avg.pow';
> cfg.maskparameter = cfg.funparameter;
> cfg.funcolorlim = [-1 2];
> cfg.opacitylim = [-5 1];
> cfg.opacitymap = 'rampup';
> ft_sourceplot(cfg, source25DiffInt);
>
> _______________________________________________
> fieldtrip mailing list
> https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> https://doi.org/10.1371/journal.pcbi.1002202
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20180709/9caa4f57/attachment-0001.html>


More information about the fieldtrip mailing list