%%
clear all

datafile='VisualSearch';
load(datafile);

image_path='.\DispImg\';

imageNo = 03222;
subNumber = 6;

TrialIndex=find(data(subNumber).images==imageNo);

mode = 'LR';
algorithm = 'IDT';

EyeMove=data(subNumber).gaze{TrialIndex};

mode = 'L';
durationThreshold = 150;
dispersionThreshold = 50;

%%

switch upper(mode)
    case 'L'
        eyeData = EyeMove(:,1:3);
    case 'R'
        eyeData = [EyeMove(:,1) EyeMove(:, 4:5)];
end


%%

d_t_EyeData = diff(eyeData(:,1));
plot(1:length(d_t_EyeData),d_t_EyeData);

unique(d_t_EyeData);

%%

figure;
subplot(3,2,1)
title('diff in time');
plot(d_t_EyeData);
xlim([1 length(d_t_EyeData)]);

subplot(3,2,2)
title('diff in time');
plot(d_t_EyeData);
xlim([1 length(d_t_EyeData)]);

subplot(3,2,[3 5])
title('diff in x and y');
d_x_EyeData = diff(eyeData(:,2));
plot(1:length(d_x_EyeData),d_x_EyeData);
hold on;
d_y_EyeData = diff(eyeData(:,3));
plot(1:length(d_y_EyeData),d_y_EyeData);
legend({'diff in x','diff in y'},'Location','SouthWest');
xlim([1 length(d_t_EyeData)]);

subplot(3,2,[4 6])
title('x and y');
plot(eyeData(:,2));
hold on;
plot(eyeData(:,3));
legend({'x','y'},'Location','SouthWest');
xlim([1 length(d_t_EyeData)]);

%%

cs_EyeData = cumsum(d_t_EyeData);
window = eyeData(cs_EyeData<durationThreshold,:)

window_dispersion = calcDispersion(window)

window_dispersion<dispersionThreshold