<div dir="ltr"><div style="font-size:14px">Hi, Tzvetan.</div><div style="font-size:14px"><br></div><div style="font-size:14px">Thank you for your first response.</div><div style="font-size:14px"><br></div><div style="font-size:14px">The weird output that you asked in first response is first attached figure 1.</div><div style="font-size:14px">I thought It was weird because the sensor level output is so different form the source level output.</div><div style="font-size:14px"><br></div><div style="font-size:14px">I revised the script regarding your first advise and following tutorial. The output is second figure.</div><div style="font-size:14px">(<a href="http://www.fieldtriptoolbox.org/example/read_neuromag_mri_and_create_single-subject_grids_in_individual_head_space_that_are_all_aligned_in_mni_space" target="_blank">http://www.fieldtriptoolbox.org/example/read_neuromag_mri_and_create_single-subject_grids_in_individual_head_space_that_are_all_aligned_in_mni_space</a>)</div><div style="font-size:14px">It is better than first output. However,<b> the sensor level output correspond with the source level output, exactly.</b></div><div style="font-size:14px"><br></div><div style="font-size:14px">Is it anything that need to revise? </div><div style="font-size:14px">If the data is plotted at the source and sensor level, What the difference between two plots can emerge in general?</div><div style="font-size:14px"><br></div><div style="font-size:14px">I really appreciate for your help in advance.</div><div style="font-size:14px"><br></div><div style="font-size:14px">Best regards,<br></div><div style="font-size:14px">MK CHOE</div><div style="font-size:14px"><br></div><div style="font-size:14px">P.S. The revised script is as follows.</div><div style="font-size:14px"><br></div><div style="font-size:14px">%% Calculating the cross spectral density matrix </div><div style="font-size:14px">cfg = [];</div><div style="font-size:14px">cfg.output = 'powandcsd';</div><div style="font-size:14px">cfg.method = 'mtmfft';</div><div style="font-size:14px">cfg.foilim = [10 10]; % alpha [8 12]</div><div style="font-size:14px">cfg.tapsmofrq = 2;</div><div style="font-size:14px">cfg.channel = {'MEG'};</div><div style="font-size:14px">cfg.keeptrials = 'yes';</div><div style="font-size:14px">fft_rest = ft_freqanalysis(cfg, data_clean);</div><div style="font-size:14px"><br></div><div style="font-size:14px">%% Computing the leadfield matrices</div><div style="font-size:14px">cfg = [];</div><div style="font-size:14px">cfg.grad = fft_rest.grad;</div><div style="font-size:14px">cfg.vol = hdm;</div><div style="font-size:14px">cfg.grid = grid;</div><div style="font-size:14px">cfg.keeptrials = 'yes';</div><div style="font-size:14px">cfg.normalize = 'yes'; </div><div style="font-size:14px">cfg.channel = {'MEG'};</div><div style="font-size:14px">rest_grid = ft_prepare_leadfield(cfg);</div><div style="font-size:14px"><br></div><div style="font-size:14px">%% Source analysis </div><div style="font-size:14px">cfg = [];</div><div style="font-size:14px">cfg.frequency = 'all';</div><div style="font-size:14px">cfg.grad = fft_rest.grad;</div><div style="font-size:14px">cfg.method = 'dics';</div><div style="font-size:14px">cfg.grid = rest_grid;</div><div style="font-size:14px">cfg.vol = hdm;</div><div style="font-size:14px">cfg.keeptrials = 'yes';</div><div style="font-size:14px">cfg.keepfilter = 'yes';</div><div style="font-size:14px">cfg.keepcsd = 'yes';</div><div style="font-size:14px">cfg.feedback = 'textbar';</div><div style="font-size:14px">cfg.dics.lambda = 0.1</div><div style="font-size:14px">cfg.dics.projectnoise = 'yes';</div><div style="font-size:14px">rest_source = ft_sourceanalysis(cfg, fft_rest);</div><div style="font-size:14px"><br></div><div style="font-size:14px">%% sourcedescriptive</div><div style="font-size:14px">cfg = [];</div><div style="font-size:14px">rest_source.dim = rest_grid.dim;</div><div style="font-size:14px">descri_source = ft_sourcedescriptives(cfg, rest_source);</div><div style="font-size:14px">descri_source.pos = t1grid.pos;</div><div style="font-size:14px"><br></div><div style="font-size:14px">%% sourceinterpolating_nai</div><div style="font-size:14px">cfg = [];</div><div style="font-size:14px">cfg.parameter = 'nai';</div><div style="font-size:14px">cfg.voxelcoord = 'no';</div><div style="font-size:14px">cfg.interpmethod = 'linear';</div><div style="font-size:14px">sdint = ft_sourceinterpolate(cfg, descri_source, norm_mri);</div><div style="font-size:14px"><br></div><div style="font-size:14px">%% plot fft results </div><div style="font-size:14px">cfgp = [];</div><div style="font-size:14px">cfgp.layout = 'neuromag306mag_rev.lay';</div><div style="font-size:14px">cfgp.interactive='yes';</div><div style="font-size:14px">colorbar;</div><div style="font-size:14px">ft_topoplotER(cfgp, fft_rest);</div><div style="font-size:14px"><br></div><div style="font-size:14px">%% source plot _nai</div><div style="font-size:14px">cfg = [];</div><div style="font-size:14px">cfg.method = 'slice';</div><div style="font-size:14px">cfg.funparameter = 'nai';</div><div style="font-size:14px">cfg.surfdownsample = 2;</div><div style="font-size:14px">cfg.projmethod = 'project';</div><div style="font-size:14px">cfg.funcolormap = 'jet';</div><div style="font-size:14px">ft_sourceplot(cfg, sdint);</div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><font color="#888888">Mi Kyung CHOE</font></div><span style="color:rgb(136,136,136)">Department of Brain & Cognitive Sciences, </span><span style="color:rgb(136,136,136)">Seoul National University </span><div><font color="#888888">Human Brain Function Laboratory<br></font><div><font color="#888888">E-mail : <a href="mailto:cmk0803@hbf.re.kr" target="_blank">cmk0803@hbf.re.kr</a>, <a href="mailto:cmk0803@meg.re.kr" target="_blank">cmk0803@meg.re.kr</a></font></div></div></div></div>
</div>