onsetFrameFiles = dir('*analysis*');
frameDataFiles = dir('*frameData.mat');

FigPos = [680 384 668 594];
for i=length(onsetFrameFiles)
    FigPos = [680 384 668 594];

    load(onsetFrameFiles(i).name);
    onsetFrame = DataMat(2,1);
    
    vidFile = dir(['*' onsetFrameFiles(i).name(21:32) '*avi']);
    
    movie_obj = VideoReader(vidFile.name); %#ok<TNMLP>
    frameRate = movie_obj.FrameRate;
    vidH = movie_obj.Height;
    vidW = movie_obj.Width;
    numFrames = onsetFrame+frameRate*30;
    
    firstFrame = onsetFrame-20*30;
    if firstFrame<1
        firstFrame=1;
    end
    
    
    totalFrames = numFrames-firstFrame;
    
   
    frame = (read(movie_obj, (firstFrame)));

    imshow(frame);
    axis([0 vidW 0 vidH])
    title('Set Screen Edges')
    screenEdges = ginput(4);
    
    
    
    Data.vidFile = vidFile;
    Data.onsetFrame = onsetFrame;
    Data.movie_obj = movie_obj;
    Data.numFrames = numFrames;
    Data.firstFrame = firstFrame;
    
    Data.screenEdges = screenEdges;
    
    screenEdges(end+1,:) = screenEdges(1,:);
    title('Find bird position')
    
    n = 1;
    figure(1)
    set(gcf,'Position', FigPos);
    axis([0 round(max(screenEdges(:,1))) 0 round(max(screenEdges(:,2)))])

    load(frameDataFiles(i-1).name);
    for k=firstFrame : numFrames

        frame = squeeze(frameData(:,:,:,n));

        imagesc(frame/255)
        hold on
        patch(screenEdges(:,1), screenEdges(:,2),'k','facecolor','w', 'facealpha', 0.7)
        
        title([num2str(k) '/' num2str(numFrames)])
        pos(n,:) = ginput(1);
        n = n+1;
        frameCount(n) = k;
        cla
    end
    
    Data.frameCount = frameCount;
    Data.pos = pos;
    saveName = [onsetFrameFiles(i).name(1:(end-13)) '_Data.mat'];
    save(saveName, 'Data')
    
    clearvars -except onsetFrameFiles i frameDataFiles
    
end