meg_file = 'sbj_test_UR.meg.mat'; % load channel position [pos, channel_info] = vb_load_channel(meg_file); % load data [data, channel_info2, time_info] =vb_load_meg_data(meg_file); % average timewindow between 80msec to 120msec trial_number = 1; time_ix = intersect(find(time_info.time > 0.08), find(time_info.time<0.12)); data_average = mean(data(:, time_ix, trial_number), 2); % extract axial channel ch_axial = find(channel_info.Type==2); % plot topography figure; vb_plot_sensor_2d(pos(ch_axial, :), data_average(ch_axial)); axis equal; axis off; colorbar % add facial contour if you want. vb_plot_sensor_2d_head_plot_add(gca);