mpib_memoreye/gaze_heatmaps/plotting_gaze_heatmap_v4.m

clear all;

 mainfolder=''

%% Loading data of interest

% load([mainfolder gaze_heatmaps/stim1.mat'])
%load([mainfolder 'gaze_heatmaps/stim2.mat'])
%load([mainfolder 'gaze_heatmaps/delay1_stim1.mat'])
%load([mainfolder 'gaze_heatmaps/delay1_stim2.mat'])
%load([mainfolder 'gaze_heatmaps/delay2_stim1.mat'])
load([mainfolder 'gaze_heatmaps/delay2_stim2.mat'])


%% Options
ave_across_object=0; % averaging acrros object? yes= 1; no =0;
plot_avg_point=0;
c=[0 0.00025]; %colorbars limits

%% Preparing for plots

hhmap_grand=squeeze(nanmean(squeeze(nanmean(hhmap,1)),2));
avg_coord_rot_grand=squeeze(nanmean(nanmean(avg_coord_rot,1),3));

%%
allobjectc={'01 lamppost','02 gazebo','01 nightstand',...
    '02 table','01 lighthouse','02 radio',...
    '02 crown','01 outdoor chair','01 candelabra'};

figure('Position',([200 900 800 780]))

counter=1;
pixtovd=0.0217; %pixel to visual degrees
windowsize=400; %pixel

colormap('hot')

for i=1:9 % per object
    subplot(3,3,counter)
    
    ctr=round(size(hhmap_grand,2)/2,-1)
    ctrw=round((windowsize)/2)
    
    p1=(ctr-ctrw)
    p2=ctr+ctrw
    
    if p1==0
        p1=1
    end
    
    if ave_across_object
        toplot=squeeze(nanmean(hhmap_grand(:,p1:p2,p1:p2),1));
    else
        toplot=squeeze((hhmap_grand(i,p1:p2,p1:p2)));
    end
    
    toplot_gausfil = imgaussfilt(toplot,2.5);
    
    sizegaz=size(toplot_gausfil);
    
    j=getitemimg(i,sizegaz,2);
        imagesc(j);
    
    hold on;
    
    if plot_avg_point
        plot(avg_coord_rot_grand(i,1),avg_coord_rot_grand(i,2), 'y.','MarkerSize',20);
        hold on
    end
    
    imagesc(toplot_gausfil,'AlphaData', .6,c)
    
    
    hold on
    plot([ctrw],[ctrw], 'k.','MarkerSize',20);
    hold off
    
    % for visual degrees
    
    vd=round((-ctrw:1:ctrw).*0.0217,3);
    for ci=1:5
        n=[-4 -2 0 2 4];
        [val,idx(ci)]=min(abs(vd-n(ci)));
    end
    % minVal=vd(idx)
        
    space=windowsize/3;
    
    set(gca,'XTick',(idx), 'XTickLabel',n,'YTick', (idx), 'YTickLabel',n,'FontName', 'Arial','FontSize',16 )
    
    %   for pixels
    
    %    space=windowsize/5;
    %    set(gca,'XTick',([1:space:windowsize windowsize]), 'XTickLabel', (-ctrw:space:ctrw),'YTick', ([1:space:windowsize windowsize]), 'YTickLabel', (-ctrw:space:ctrw))
    %
    counter=counter+1;
    ylabel('° vis. angle');
    xlabel('° vis. angle');
    
end