[FieldTrip] EEG Electrode / Volume Alignment Problem

Julian Keil julian.keil at gmail.com
Fri Dec 20 10:07:42 CET 2013


Hi Stephan,

Ah, I forgot that the electrode definition is added to the EEG-structure.
Unfortunately, neither removing this nor setting all .elec-fields to the aligned elec_new helps.

Happy end-of-the-year wishes,

Julian

Am 19.12.2013 um 20:44 schrieb smoratti at psi.ucm.es:

> 
> I thought of changing the dat4beam.elec to newelec. Does this help?
> 
> Stephan
> 
> ________________________________________________________
> Stephan Moratti, PhD
> 
> see also: http://web.me.com/smoratti/
> 
> Universidad Complutense de Madrid
> Facultad de Psicología
> Departamento de Psicología Básica I
> Campus de Somosaguas
> 28223 Pozuelo de Alarcón (Madrid)
> Spain
> 
> and
> 
> Center for Biomedical Technology
> Laboratory for Cognitive and Computational Neuroscience
> Parque Científico y Tecnológico de la Universidad Politecnica de Madrid
> Campus Montegancedo
> 28223 Pozuelo de Alarcón (Madrid)
> Spain
> 
> 
> email: smoratti at psi.ucm.es
> Tel.:    +34 679219982
> 
> El 19/12/2013, a las 17:08, Julian Keil escribió:
> 
>> Hi Stephan,
>> 
>> thanks for your quick reply.
>> 
>> If I change the cfg.elec definition in the call to the beamformer, the source projection follows the electrodes, i.e. when I use the misaligned electrodes, the source is in the visual cortex, and if I use the aligned electrodes, the source is left temporal.
>> 
>> If I leave the info for cfg.elec out, ft_sourceanalysis takes the electrode info from the lead field.
>> 
>> Am I maybe missing some step during the call to ft_prepare_leadfield? Maybe some rotation of the volume or something like this?
>> 
>> Thanks a lot once more.
>> 
>> Julian
>> 
>> Am 19.12.2013 um 16:35 schrieb smoratti at psi.ucm.es:
>> 
>>> 
>>> hmm… in theory every thing should be fine as the lead field is calculated with aligned electrodes to bnd_s. However, try to change the electrode positions as well in dat4beam. Does this help?
>>> 
>>> best,
>>> 
>>> Stephan
>>> 
>>> ________________________________________________________
>>> Stephan Moratti, PhD
>>> 
>>> see also: http://web.me.com/smoratti/
>>> 
>>> Universidad Complutense de Madrid
>>> Facultad de Psicología
>>> Departamento de Psicología Básica I
>>> Campus de Somosaguas
>>> 28223 Pozuelo de Alarcón (Madrid)
>>> Spain
>>> 
>>> and
>>> 
>>> Center for Biomedical Technology
>>> Laboratory for Cognitive and Computational Neuroscience
>>> Parque Científico y Tecnológico de la Universidad Politecnica de Madrid
>>> Campus Montegancedo
>>> 28223 Pozuelo de Alarcón (Madrid)
>>> Spain
>>> 
>>> 
>>> email: smoratti at psi.ucm.es
>>> Tel.:    +34 679219982
>>> 
>>> El 19/12/2013, a las 15:38, Julian Keil escribió:
>>> 
>>>> Dear all,
>>>> 
>>>> I have a problem of which I can't really find the source.
>>>> 
>>>> When I use the standard BEM-Model and the standard MRI together with our vendor-provided EEG electrode positions, I can nicely and reliably localize effects in source space (e.g. the visual N1, see attached source plot).
>>>> 
>>>> However, when I use an individual MRI (DICOM, imported via ft_read_mri, realigned interactively, resliced and segmented) to build a BEM-Model (Dipoli-Method), realign the electrodes to the new volume (using the interactive mode in ft_electroderealign) and subsequently build an individual lead field, the effect is shifted away from visual cortex and it appears as if the source and the MRI are misaligned.
>>>> 
>>>> The odd thing is, that if I leave the alignment step out, i.e. use the misaligned electrodes to build the lead field, the effect at least moves back to the visual cortex.
>>>> 
>>>> I hope the attached comparison between aligned and misaligned electrodes makes this a bit clearer.
>>>> 
>>>> Does anyone have an idea at which point this problem might occur?
>>>> 
>>>> Thanks a lot.
>>>> 
>>>> Julian
>>>> <EEG_Vol_Alignment.pdf>
>>>> <N1_standard_brain.png>
>>>> 
>>>> P.S.: Below are the basic steps I use:
>>>> 
>>>> %% 2.1. Build individual BEM-Model
>>>> % First Realign Volume. Set Nas LPA, RPA, and positive z-values
>>>> 
>>>> cfg=[];
>>>> cfg.method='interactive';
>>>> cfg.coordsys = 'ctf';
>>>> 
>>>> mri_r = ft_volumerealign(cfg,mri);
>>>> 
>>>> cfg=[];
>>>> 
>>>> mri_rs = ft_volumereslice(cfg,mri_r);
>>>> 
>>>> %% Segment the individual volume
>>>> cfg=[];
>>>> cfg.output={'brain' 'scalp' 'skull'};
>>>> cfg.brainsmooth = 5;
>>>> cfg.scalpsmooth = 5;
>>>> cfg.brainthreshold = .5;
>>>> cfg.scalpthreshold = .1;
>>>> mri_s = ft_volumesegment(cfg,mri_rs);
>>>> 
>>>> %% Create Headmodel
>>>> cfg=[];
>>>> cfg.interactive = 'no';
>>>> cfg.numvertices = 1000;
>>>> 
>>>> bnd_s = ft_prepare_mesh(cfg,mri_s);
>>>> 
>>>> cfg=[];
>>>> cfg.method='dipoli'; %
>>>> cfg.conductivity=[0.3300 0.0041 0.3300];
>>>> 
>>>> vol_s = ft_prepare_headmodel(cfg,bnd_s);
>>>> 
>>>> cfg = [];
>>>> cfg.mri            = mri_s;
>>>> grid               = ft_prepare_sourcemodel(cfg);
>>>> 
>>>> grid.inside = find(grid.brain==1)';
>>>> grid.outside = find(grid.brain==0)';
>>>> grid.pos = grid.pos*10; % Set right dimensions
>>>> 
>>>> %% Align Electrodes
>>>> cfg=[];
>>>> cfg.method='interactive';
>>>> cfg.elec=elec_128;
>>>> cfg.headshape=vol_s.bnd(1);
>>>> 
>>>> elec_new=ft_electroderealign(cfg);
>>>> 
>>>> %% 4. Make the Leadfield.
>>>> elec_new.type='eeg';
>>>>  
>>>> cfg=[];
>>>> cfg.channel = 1:126;
>>>> cfg.elec=elec_new;
>>>> cfg.vol=vol_s;
>>>> cfg.grid=grid;
>>>> 
>>>> lf = ft_prepare_leadfield(cfg);
>>>> 
>>>> %% BEAMFORMER 
>>>> 
>>>> cfg=[];
>>>> cfg.channel=1:126;
>>>> cfg.method='lcmv'; %Time-Domain-Beamformer
>>>> cfg.grid=lf;
>>>> cfg.vol=vol_s; 
>>>> cfg.elec=elec_new; 
>>>> cfg.reducerank  = 3; 
>>>> cfg.keepfilter='yes'; 
>>>> 
>>>> beamfilter = ft_sourceanalysis(cfg, dat4beam);
>>>> 
>>>> %% Keep the cfg, just add the filter.
>>>> 
>>>> cfg.grid.filter = beamfilter.avg.filter; % Now we'll use the filter defined above
>>>> 
>>>> lcmv_bl = ft_sourceanalysis(cfg, bl_dat);
>>>> lcmv_act = ft_sourceanalysis(cfg, act_dat);
>>>> 
>>>> 
>>>> 
>>>> ********************
>>>> Dr. Julian Keil
>>>> 
>>>> AG Multisensorische Integration
>>>> Psychiatrische Universitätsklinik
>>>> der Charité im St. Hedwig-Krankenhaus
>>>> Große Hamburger Straße 5-11, Raum E 307
>>>> 10115 Berlin
>>>> 
>>>> Telefon: +49-30-2311-1879
>>>> Fax:        +49-30-2311-2209 
>>>> http://psy-ccm.charite.de/forschung/bildgebung/ag_multisensorische_integration
>>>> 
>>>> _______________________________________________
>>>> fieldtrip mailing list
>>>> fieldtrip at donders.ru.nl
>>>> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>>> 
>>> _______________________________________________
>>> fieldtrip mailing list
>>> fieldtrip at donders.ru.nl
>>> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>> 
>> _______________________________________________
>> fieldtrip mailing list
>> fieldtrip at donders.ru.nl
>> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> 
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20131220/6a4c2371/attachment-0002.html>


More information about the fieldtrip mailing list