<div dir="ltr">Hello All Fieldtrippers!<div><br></div><div>First i will tell you what my data is: I have epoched EEG data, where in each epoch participants task was to imagine the situations presented on screen. Each epoch is 5 seconds long, containing fixation cross and tag for imagination. </div><div><br></div><div>What i want: I want to perform source analysis on my preprocessed epoched data, then i want to parcellate sources to atlas, then i want to create a matrix where i will have reconstructed source data for each ROI from atlas. I need my data to be epoched (time locked, but not ERP), and not power, since i want to do multiple operations on this matrix. So the matrix will look like</div><div>ROI 1 signal value for 1st time point | signal value for 2nd time point .....</div><div>ROI 2 signal value for 1st time point | signal value for 2nd time point .....</div><div>ROI 3 signal value for 1st time point | signa value for 2nd time point .....</div><div>...</div><div><br></div><div>What i did: I've done preprocessing in eeglab, computed headmodel (im using template MRI from Fieldtrip repsitory), sourcemodel, leadfield and sourcemodel (i will paste code below), then i wanted to compute sourceinterpolate and sourceparcellate (i've used some tutorials from Your page), and then i wanted to create a matrix.</div><div><br></div><div>What is my problem: I faced first trouble after running sourceanalysis. I've check 'source.avg.filter' and 'source.avg.mom' field and is looks like this:<br>  {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {3×19 double}<br>    {3×19 double}<br>    {3×19 double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br>    {0×0  double}<br></div><div>Next, i've tried running sourceinterpolate on my source data and atlas, it worked without errors. I passed <span style="font-family:Menlo,Monaco,Consolas,"Courier New",monospace;font-size:10pt">cfg.parameter = </span><span style="font-family:Menlo,Monaco,Consolas,"Courier New",monospace;font-size:10pt;color:rgb(167,9,245)">'trial'</span><span style="font-family:Menlo,Monaco,Consolas,"Courier New",monospace;font-size:10pt">; </span>since i thought it would interpolate data as time locked epochs (trials). </div><div><br></div><div>Unfortunately, next function i wanted to use - sourceparcellate can not parcellate data like that, as i've seen it works only with 'pow' parameter, and i dont want power, since its like average power for ROI. <br>It only works like that:<br><div style="padding:30px"><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg = [];</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg.interpmethod = <span style="color:rgb(167,9,245)">'nearest'</span>;</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg.parameter = <span style="color:rgb(167,9,245)">'all'</span>; <span style="color:rgb(0,128,19)">%'avg.pow'</span></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg.atlas = atlas;</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg.keeptrials = <span style="color:rgb(167,9,245)">'yes'</span>;</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg.mri = mri_resliced;</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">source_int = ft_sourceinterpolate(cfg, source, atlas);</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg = [];</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg.parameter = <span style="color:rgb(167,9,245)">'all'</span>;</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg.method = <span style="color:rgb(167,9,245)">'mean'</span>;</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">cfg.atlas = atlas;</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace">source_parcel = ft_sourceparcellate(cfg, source_int,atlas);</div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small"><br></span></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small">but when i tried to create matrix on source_parcel data, i got a lot of NaNs:<br></span>    1.3660    1.3660    1.3660    1.3660    1.3660    1.3660    1.3660    1.3660<br>    0.4568    0.4568    0.4568    0.4568    0.4568    0.4568    0.4568    0.4568<br>    0.8530    0.8530    0.8530    0.8530    0.8530    0.8530    0.8530    0.8530<br>       NaN       NaN       NaN       NaN       NaN       NaN       NaN       NaN<br>       NaN       NaN       NaN       NaN       NaN       NaN       NaN       NaN<br>       NaN       NaN       NaN       NaN       NaN       NaN       NaN       NaN<br>    0.4486    0.4486    0.4486    0.4486    0.4486    0.4486    0.4486    0.4486<br>       NaN       NaN       NaN       NaN       NaN       NaN       NaN       NaN<br>       NaN       NaN       NaN       NaN       NaN       NaN       NaN       NaN<span style="font-family:Arial,Helvetica,sans-serif;font-size:small"><br></span></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><br></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small">So, in summary, i need some help with usage of sourceinterpolate and sourceparcellate on epoched data, so i will be able to create matrix containing signal, not one power value per ROI.</span><span style="font-family:Arial,Helvetica,sans-serif;font-size:small"> I haven't seen anyone do this on the data structure that i'm working with, so i feel a little bit lost by now.</span></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small">If you spot any mistakes prior to sourceinterpolate please let me know, i'm doing pipeline like that first time in my life so i know it can contain mistakes.</span></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small">I would appreciate any help.</span></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small"><br></span></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small">Thanks in advance and have a good day!</span></div><div style="font-size:10pt;font-family:Menlo,Monaco,Consolas,"Courier New",monospace"><span style="font-family:Arial,Helvetica,sans-serif;font-size:small">Monika<br><br>%%%%code<br></span><div style="padding:30px"><div style="font-size:10pt"><span style="white-space-collapse: preserve;">fid = fopen(electrode_file_path, <span style="color:rgb(167,9,245)">'r'</span>);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><span style="color:rgb(14,0,255)">if </span>fid == -1</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    error(<span style="color:rgb(167,9,245)">'Failed to open file: %s'</span>, electrode_file_path);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><span style="color:rgb(14,0,255)">end</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">file_content = textscan(fid, <span style="color:rgb(167,9,245)">'%s %s %f %f %f %f %f %f %f %f'</span>, <span style="color:rgb(167,9,245)">'HeaderLines'</span>, 1);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">fclose(fid);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">labels = file_content{2};</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">X = file_content{5};</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">Y = file_content{6};</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">Z = file_content{7};</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">elec = struct();</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">elec.label = labels;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">elec.pnt = [X, Y, Z];</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">n_elec = length(labels);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">elec.chantype = repmat({<span style="color:rgb(167,9,245)">'eeg'</span>}, n_elec, 1);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">elec.chanunit = repmat({<span style="color:rgb(167,9,245)">'V'</span>}, n_elec, 1);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">elec.unit = <span style="color:rgb(167,9,245)">'mm'</span>;</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><br></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.method = <span style="color:rgb(167,9,245)">'fiducial'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.coordsys = <span style="color:rgb(167,9,245)">'ctf'</span>; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">mri_realigned = ft_volumerealign(cfg, mri);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.method = <span style="color:rgb(167,9,245)">'flip'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.coordsys = <span style="color:rgb(167,9,245)">'ctf'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">mri_resliced = ft_volumereslice(cfg, mri_realigned);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg           = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.output    = {<span style="color:rgb(167,9,245)">'brain'</span>, <span style="color:rgb(167,9,245)">'skull'</span>, <span style="color:rgb(167,9,245)">'scalp'</span>};</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.coordsys = <span style="color:rgb(167,9,245)">'ctf'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">segmentedmri  = ft_volumesegment(cfg, mri_resliced);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><br></span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg        = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.method = <span style="color:rgb(167,9,245)">'hexahedral'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.coordsys = <span style="color:rgb(167,9,245)">'ctf'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">mesh = ft_prepare_mesh(cfg, segmentedmri);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.method = <span style="color:rgb(167,9,245)">'simbio'</span>; <span style="color:rgb(0,128,19)">%method used in artcile</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.conductivity = [0.33, 0.0125, 0.33]; <span style="color:rgb(0,128,19)">% brain skull scalp</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">headmodel = ft_prepare_headmodel(cfg, mesh);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.headmodel = headmodel;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.headmodel.type=<span style="color:rgb(167,9,245)">'simbio'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.xgrid  = <span style="color:rgb(167,9,245)">'auto'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.ygrid  = <span style="color:rgb(167,9,245)">'auto'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.zgrid  = <span style="color:rgb(167,9,245)">'auto'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.unit   = <span style="color:rgb(167,9,245)">'mm'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.tight  = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.inwardshift  = -1.5;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.resolution  = 6; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">sourcemodel   = ft_prepare_sourcemodel(cfg);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><br></span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">atlas = ft_read_atlas(fullfile(ftpath, </span><span style="font-size:10pt;color:rgb(167,9,245)">'template\atlas\aal\ROI_MNI_V4.nii'</span><span style="font-size:10pt">))</span><br></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">atlas = ft_convert_units(atlas,<span style="color:rgb(167,9,245)">'mm'</span>);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">sourcemodel = ft_convert_units(sourcemodel,<span style="color:rgb(167,9,245)">'mm'</span>);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">file_list = dir(fullfile(PATH, <span style="color:rgb(167,9,245)">'*.mat'</span>));</span></div><div style="font-size:10pt"><span style="font-size:10pt;color:rgb(14,0,255)">for </span><span style="font-size:10pt">idx = 1:length(file_list)</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    data_load = load(fullfile(PATH, file_list(idx).name));</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    EEG = data_load.EEG2;</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">    data = struct();</span><br></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">    data.label = {EEG.chanlocs.labels}';</span><br></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">    num_trials = EEG.trials;</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    num_points = EEG.pnts;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    srate = EEG.srate;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    time = linspace(EEG.xmin, EEG.xmax, num_points);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    data.time = repmat({time}, 1, num_trials);</span></div><div style="font-size:10pt"><span style="font-size:10pt">    data.trial = squeeze(num2cell(EEG.data, [1, 2]));</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    data.trial = data.trial';</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">    sample_info_start = repmat((0:num_trials-1)' * num_points, 1, 1);</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    sample_info_end = sample_info_start + num_points - 1;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    data.sampleinfo = [sample_info_start, sample_info_end];</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">   data.trialinfo = [EEG.epoch.eventtype]'; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    data.hdr = EEG;</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">    data.elec = elec;</span><br></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">    data.cfg = [];</span><br></div><div style="font-size:10pt"><span style="font-size:10pt"><br></span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">   cfg = [];</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.method = <span style="color:rgb(167,9,245)">'project'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.headshape = mesh;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.mri=mri_resliced;</span></div><div style="font-size:10pt"><span style="font-size:10pt">   elec_aligned = ft_electroderealign(cfg,elec);</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"> </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   elec_good = ft_prepare_vol_sens(headmodel, elec_aligned)</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   data.elec = elec_good;</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.trials = <span style="color:rgb(167,9,245)">'all'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.keeptrials = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.covariance = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   data.elec=data.elec.elec; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   timlock = ft_timelockanalysis(cfg,data);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt"><br></span></div><div style="font-size:10pt"><span style="font-size:10pt">   cfg = [];</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.elec      = timlock.elec;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.channel   = data.label;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.headmodel = headmodel; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.sourcemodel.pos    = sourcemodel.pos; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.grid = sourcemodel;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.channel = {<span style="color:rgb(167,9,245)">'EEG'</span>};</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   leadfield = ft_prepare_leadfield(cfg);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><br></span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">    cfg = [];</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    cfg.interpmethod = <span style="color:rgb(167,9,245)">'nearest'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    cfg.parameter = <span style="color:rgb(167,9,245)">'all'</span>;</span></div><div style="font-size:10pt"><span style="font-size:10pt">    cfg.atlas=atlas;</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    cfg.mri=mri_resliced;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    atlas_int = ft_sourceinterpolate(cfg, sourcemodel, atlas);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="font-size:10pt">  </span></div><div style="font-size:10pt"><span style="font-size:10pt">   cfg = [];</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.method       = <span style="color:rgb(167,9,245)">'lcmv'</span>; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.elec         = timlock.elec; <span style="color:rgb(0,128,19)">;</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.sourcemodel = sourcemodel;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.grid = sourcemodel;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.sourcemodel.leadfield = leadfield;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.headmodel    = headmodel;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.projectnoise = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.keepmom      = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.keepfilter   = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.keeptrials = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   source = ft_sourceanalysis(cfg, timlock);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><br></span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg=[];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   cfg.keeptrials = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">   source = ft_sourcedescriptives(cfg, source); </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><br></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><span style="color:rgb(14,0,255)">if </span>~isfield(source, <span style="color:rgb(167,9,245)">'filter'</span>)</span></div><div style="font-size:10pt"><span style="font-size:10pt">    source.filter = cell(size(source.pos, 1), 1);</span><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    <span style="color:rgb(14,0,255)">for </span>i = 1:length(source.filter)</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">        <span style="color:rgb(14,0,255)">if </span>source.inside(i)</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">            source.filter{i} = leadfield.leadfield{i} * source.avg.mom{i};</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">        <span style="color:rgb(14,0,255)">else</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">            source.filter{i} = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">        <span style="color:rgb(14,0,255)">end</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    <span style="color:rgb(14,0,255)">end</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><span style="color:rgb(14,0,255)">end</span></span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.interpmethod = <span style="color:rgb(167,9,245)">'nearest'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.parameter = <span style="color:rgb(167,9,245)">'trial'</span>; <span style="color:rgb(0,128,19)">%'avg.pow'</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.atlas = atlas;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.keeptrials = <span style="color:rgb(167,9,245)">'yes'</span>;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.mri = mri_resliced;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">source_int = ft_sourceinterpolate(cfg, source, atlas);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><br></span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg = [];</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.parameter = <span style="color:rgb(167,9,245)">'all'</span>; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.method = <span style="color:rgb(167,9,245)">'mean'</span>; </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">cfg.atlas = atlas;</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">source_parcel = ft_sourceparcellate(cfg, source_int,atlas);</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><br></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">%%big matrix creation works only for 'pow' parameter</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">num_regions = length(source_parcel.label); </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">num_epochs = size(EEG.data, 3);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><br></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">big_matrix = zeros(num_regions, num_epochs);</span></div><div style="font-size:10pt"></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><span style="color:rgb(14,0,255)">for </span>region_idx = 1:num_regions</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    <span style="color:rgb(14,0,255)">for </span>epoch_idx = 1:num_epochs</span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">        big_matrix(region_idx, epoch_idx) = source_parcel.pow(region_idx); </span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;">    <span style="color:rgb(14,0,255)">end</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><span style="color:rgb(14,0,255)">end</span></span></div><div style="font-size:10pt"><span style="white-space-collapse: preserve;"><span style="color:rgb(14,0,255)">end</span></span></div></div></div></div></div></div>