[FieldTrip] dipole simulation on head model

gamaliel huerta urrea gamaliel.ghu at gmail.com
Tue Sep 16 19:20:36 CEST 2014


Hi all

I followed the tips above to simulate a dipole source on my own head
model generated from MRIs. However, using information from:
http://fieldtrip.fcdonders.nl/example/compute_forward_simulated_data_and_apply_a_dipole_fit

I could not simulate anything. Only one image is blank. The code I am
using is as follows:

cfg = [];
cfg.vol = vol; % See above
cfg.elec = elec; % See above
cfg.dip.pos = [0 0 0];
cfg.dip.mom = [2 2 2]; % Note, it Should be transposed
cfg.dip.frequency = 10;
cfg.ntrials = 10;
cfg.triallength = 1; % seconds
cfg.fsample = 25; % Hz
raw1 = ft_dipolesimulation (cfg);
avg1 = ft_timelockanalysis([], raw1);

plot (avg1.time, avg1.avg); % Plot the timecourse

where:

vol = ft_prepare_headmodel volumes generated from the segmentation of
my resonances.

and

ft_read_sens elec = ('standard_1020.elc');

to run the script the result is a blank image.

and Matlab tells me the following:

using headmodel specified in the configuration
using electrodes specified in the configuration
Determining source compartment (1)
projecting electrodes on skin surface
electrode and system combine combining transfer matrix
computing simulated data
computing simulated trial data for 10

RMS value of simulated data is NaN
the call to "ft_dipolesimulation" took 4 seconds
the input is raw Data with 97 channels and 10 trials
Warning: the data does not Contain a definition trial
>* In utilities \ private \ warning_once at 158
*   In utilities \ private \ fixsampleinfo at 66
   In ft_datatype_raw at 157
   In ft_checkdata at 224
   In ft_timelockanalysis at 105
   In codgen at 50
Warning: sampleinfo by reconstructing Assuming That the trials are
consecutive segments of
a continuous recording
>* In utilities \ private \ warning_once at 158
*   In utilities \ private \ fixsampleinfo at 79
   In ft_datatype_raw at 157
   In ft_checkdata at 224
   In ft_timelockanalysis at 105
   In codgen at 50
averaging trials
trial averaging 10 of 10

the call to "ft_timelockanalysis" took 1 seconds


I hope you can help


greetings

-- 
*Gamaliel Huerta*
*Ingeniería Civil Biomédica*
*Universidad de Valparaíso*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20140916/cf2e0fd5/attachment-0002.html>


More information about the fieldtrip mailing list