<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div><br></div><div>Allright then,</div><div><br><blockquote type="cite"><div dir="ltr"><div>Sorry for my continuous questions, but i am new in this world, the final objective of this is try of generate my own sources to generate my own simulate eeg data, and next aplicate the source analysis to found the same position of the dipole</div></div></blockquote>here is what I would do using the standard head model and standard electrodes provided with fieldtrip. Just exchange the head model with your desired head model and make sure that the units match and the alignment of elec’s and vol is correct. </div><div>good luck</div><div>tzvetan</div><div><br></div><div><br></div><div><br></div><div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(178, 69, 243);"><span style="color: #000000">load (</span>'/Users/zez/Documents/matlabtools/fieldtrip-20140518/template/headmodel/standard_bem'<span style="color: #000000">);</span></div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(178, 69, 243);"><span style="color: #000000">load (</span>'/Users/zez/Documents/matlabtools/fieldtrip-20140518/template/headmodel/standard_mri'<span style="color: #000000">);</span></div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(178, 69, 243);"><span style="color: #000000">elec = ft_read_sens(</span>'/Users/zez/Documents/matlabtools/fieldtrip-20140518/template/electrode/standard_1020.elc'<span style="color: #000000">);</span></div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);">%%</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">ft_plot_mesh(vol.bnd(1), <span style="color: #b245f3">'facecolor'</span>,[0.2 0.2 0.2], <span style="color: #b245f3">'facealpha'</span>, 0.3, <span style="color: #b245f3">'edgecolor'</span>, [1 1 1], <span style="color: #b245f3">'edgealpha'</span>, 0.05);</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">hold <span style="color: #b245f3">on</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">ft_plot_mesh(vol.bnd(2),<span style="color: #b245f3">'edgecolor'</span>,<span style="color: #b245f3">'none'</span>,<span style="color: #b245f3">'facealpha'</span>,0.4);</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">hold <span style="color: #b245f3">on</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">ft_plot_mesh(vol.bnd(3),<span style="color: #b245f3">'edgecolor'</span>,<span style="color: #b245f3">'none'</span>,<span style="color: #b245f3">'facecolor'</span>,[0.4 0.6 0.4]);</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">hold <span style="color: #b245f3">on</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">ft_plot_sens(elec,<span style="color: #b245f3">'style'</span>, <span style="color: #b245f3">'sk'</span>);</div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);">%%</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg=[];</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.method = <span style="color: #b245f3">'ortho'</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">figure;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">ft_sourceplot(cfg,mri);</div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);">%%</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg=[];</div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);"><span style="color: #000000">cfg.dip.pos = [55 -7 5];</span>%dipole located in the right heschl gyrus</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.dip.mom = [0 0 -1];</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.dip.frequency = 10;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.dip.phase     = 2;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.dip.amplitude = 1;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.ntrials       = 20;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.triallength   = 2;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.fsample = 100;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.vol     = vol;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.elec    = elec;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.relnoise = 0.1;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">data  = ft_dipolesimulation(cfg);</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">data.elec = elec;</div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);">%%</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg              = [];</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.output       = <span style="color: #b245f3">'pow'</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.channel      = <span style="color: #b245f3">'EEG'</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.method       = <span style="color: #b245f3">'mtmconvol'</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.taper        = <span style="color: #b245f3">'hanning'</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.foi          = 2:2:30;                         <span style="color: #25992d">% analysis 2 to 30 Hz in steps of 2 Hz </span></div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.t_ftimwin    = ones(length(cfg.foi),1).*0.5;   <span style="color: #25992d">% length of time window = 0.5 sec</span></div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);"><span style="color: #000000">cfg.toi          = 0:0.05:2;                  </span>% time window "slides" from 0 to 2 sec in steps of 0.05 sec (50 ms)</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">TFRhann = ft_freqanalysis(cfg, data);</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">TLK = ft_timelockanalysis([],data);</div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);">%% plot the simulated signal in time and time-freq space</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg = [];   </div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(178, 69, 243);"><span style="color: #000000">cfg.channel = {</span>'FC2'<span style="color: #000000">, </span>'FC4'<span style="color: #000000">, </span>'FC6'<span style="color: #000000">, </span>'FT8'<span style="color: #000000">, </span>'C2'<span style="color: #000000">, </span>'C4'<span style="color: #000000">, </span>'C6'<span style="color: #000000">, </span>'T8'<span style="color: #000000">, </span>'T4'<span style="color: #000000">};</span></div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.renderer = <span style="color: #b245f3">'zbuffer'</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.xlim = [0.5 1.5];</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">figure </div><div style="margin: 0px; font-size: 8px; font-family: Courier;">subplot(2,2,1);ft_singleplotTFR(cfg, TFRhann);</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">subplot(2,2,2);ft_singleplotER(cfg,TLK);</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.channel = {<span style="color: #b245f3">'EEG'</span>};</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">subplot(2,2,3);ft_topoplotTFR(cfg, TFRhann);</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">subplot(2,2,4);ft_topoplotER(cfg,TLK);</div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);">%%</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg = [];</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.latency = [0.5 1.5];  <span style="color: #25992d">% specify latency window </span></div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.numdipoles = 1;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.vol = vol;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.grid.resolution = 10;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.grid.unit = <span style="color: #b245f3">'mm'</span>;</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">dip = ft_dipolefitting(cfg, TLK);</div><div style="margin: 0px; font-size: 8px; font-family: Courier; color: rgb(37, 153, 45);">%% plot the location of the dipole</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg = [];</div><div style="margin: 0px; font-size: 8px; font-family: Courier;">cfg.location = dip.dip.pos;  </div><div style="margin: 0px; font-size: 8px; font-family: Courier;">figure; ft_sourceplot(cfg, mri)</div><div><br></div><blockquote type="cite"><div dir="ltr"><div><br></div><div><br></div><div>best </div><div><br></div>-- <br><div dir="ltr"><div style="text-align:center"><font color="#20124d"><b>Gamaliel Huerta</b></font></div><div style="text-align:center"><font color="#20124d"><b>Ingeniería Civil Biomédica</b></font></div><div style="text-align:center"><font color="#20124d"><b>Universidad de Valparaíso</b></font></div></div>
</div>
_______________________________________________<br>fieldtrip mailing list<br><a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</blockquote></div><br></body></html>