Hi Martin,
Could I get your help again? I generated the below script from Brainstorm to cover some pre-processing steps for all my files. I trialed it on one raw file. The output looked fine but I got the below warning under 'process_resample' (see screenshot). The same message was listed numerous times. Could you advise if the code worked appropriately?
Also, I had originally tried to include ICA after the PSD step but the PSD and ICA steps don't seem to run together. I assume placing 'resample(256Hz)' after 'tag' under '% Process: Select data files in: /' will allow an ICA script to identify the filter file?
% Process: Band-pass:0.5Hz-30Hz
sFiles = bst_process('CallProcess', 'process_bandpass', sFiles, , ...
'sensortypes', 'EEG', ...
'highpass', 0.5, ...
'lowpass', 30, ...
'tranband', 0, ...
'attenuation', 'strict', ... % 60dB
'ver', '2019', ... % 2019
'mirror', 0, ...
'read_all', 1);
% Process: Resample: 256Hz
sFiles = bst_process('CallProcess', 'process_resample', sFiles, , ...
'freq', 256, ...
'read_all', 0);
% Process: Re-reference EEG
sFiles = bst_process('CallProcess', 'process_eegref', sFiles, , ...
'eegref', 'AVERAGE', ...
'sensortypes', 'EEG');
% Process: Power spectrum density (Welch)
sFiles = bst_process('CallProcess', 'process_psd', sFiles, , ...
'timewindow', [0, 2398.999756], ...
'win_length', 4, ...
'win_overlap', 50, ...
'sensortypes', 'EEG', ...
'win_std', 0, ...
'edit', struct(...
'Comment', 'Power', ...
'TimeBands', , ...
'Freqs', , ...
'ClusterFuncTime', 'none', ...
'Measure', 'power', ...
'Output', 'all', ...
'SaveKernel', 0));
% Save and display report
ReportFile = bst_report('Save', sFiles);
bst_report('Open', ReportFile);
% bst_report('Export', ReportFile, ExportDir);
