check_denoised_meg

 Check denoised MEG data

 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function check_denoised_meg(p)
0002 % Check denoised MEG data
0003 %
0004 % Copyright (C) 2011, ATR All Rights Reserved.
0005 % License : New BSD License(see VBMEG_LICENSE.txt)
0006 
0007 % Make directory to save figures
0008 save_dir = fullfile(p.fig_root, mfilename);
0009 if exist(save_dir, 'dir') ~= 7
0010     vb_mkdir(save_dir);
0011 end
0012 
0013 % Set files
0014 file=1;
0015 original_file = fullfile(p.proj_root, p.meg_dirname, p.load_dirname, [p.file_list{file} '.meg.mat']);
0016 denoised_file = fullfile(p.proj_root, p.meg_dirname, p.denoise_dirname, [p.file_list{file} '.meg.mat']);
0017 
0018 % Load original and denoised MEG data
0019 [original_data, ~, time_info] = vb_load_meg_data(original_file);
0020 denoised_data = vb_load_meg_data(denoised_file);
0021 
0022 % Calculate power spectrum
0023 p_o = vb_fftpower(original_data, time_info.sample_frequency);
0024 [p_d, f] = vb_fftpower(denoised_data, time_info.sample_frequency);
0025 
0026 % Plot power spectrum
0027 h = figure;
0028 plot(f, log([mean(p_o,1)', mean(p_d,1)']))
0029 xlim([0 100])
0030 xlabel('Frequency [Hz]')
0031 ylabel('log(power)')
0032 title('Power spectrum')
0033 legend('Original MEG', 'Denoised MEG')
0034 vb_savefig_as_shown(h, fullfile(save_dir, p.sub))
0035 
0036