[FieldTrip] Interpolation for power values in surface source space.

[Park, Jinhan]

Hello,

I am Jinhan Park, a third-year Ph.D. Candidate at the University of Florida. I am currently working on the analysis of resting-state EEG data with eyes closed in the source space. For the data analysis, I employed EEGLAB for preprocessing, ICA, and dipole fitting with independent components. The processed datasets were then transformed into the FieldTrip format for source localization (eLORETA).

I created a source model using 'ft_prepare_leadfield' and obtained distributed dipoles' power values using 'ft_sourceanalysis'. With this data, I attempted to visualize a 3D brain figure using 'ft_sourceplot' with the surface method. However, I noticed that certain areas of the brain appeared empty, and despite their size being relatively small, it's causing some inconvenience. I've heard that this issue can potentially be addressed by interpolating the distributed values, but I'm uncertain about how to proceed.

Could you please let me know how to fix this problem?

By the way, I wrote down some codes I did.
% Source model
vol = load('-mat', EEG.dipfit.hdmfile);% Load headfile from EEGLAB data processed by DIFIT
cfg            = [];
cfg.elec       = dataPre.elec;% Channels' information
cfg.headmodel  = vol.vol;% Headmodel
cfg.resolution = 10;   % use a 3-D grid with a 1 cm resolution
cfg.unit       = 'mm';% Unit of voxels
cfg.channel    = { 'all' };% Channels of interest
cfg.normalize  ='yes';
sourcemodel = ft_prepare_leadfield(cfg);% Create a source model (leadfield)
%% Inverse solution
% eLORETA
cfg             = [];
cfg.method      = 'eloreta';% Method for source localization
cfg.sourcemodel = sourcemodel;% Put the leadfield
cfg.headmodel   = vol.vol;% Headmodel
source          = ft_sourceanalysis(cfg, freq_pow);% Conduct the inverse solution. freq_pow is the power values in the frequency domain.
source_paf = source;
source_paf.avg.pow = dominant_f;% Change original PSD values (1 to 40 Hz) to peak frequencies
cfg = [];
cfg.method         = 'surface';% 3D brain image
cfg.funparameter   = 'pow';% Represent "sourceoi.avg.pow"
cfg.maskparameter  = cfg.funparameter;
cfg.funcolormap    = 'jet';% Types of the colorbar
cfg.projmethod     = 'nearest';
cfg.surffile       = 'surface_white_both.mat';% Brain tamplate
cfg.surfdownsample = 10;% Downsample the voxels
ft_sourceplot(cfg, source_paf);% Plot the eLORETA figure
I would sincerely appreciate it if you could answer my question.

Best regards,
Jinhan




Jinhan Park (Hans) | Ph.D. candidate

Department of Applied Physiology and Kinesiology

University of Florida

Laboratory of Rehabilitation Neuroscience | lrnlab.org

Phone: +1 (352) 328-6605

Email: jinhan.park at ufl.edu


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20240117/d20efa5d/attachment.htm>