function varargout = tagtravel2(varargin)
% TAGTRAVEL2 M-file for tagtravel2.fig
%      TAGTRAVEL2, by itself, creates a new TAGTRAVEL2 or raises the existing
%      singleton*.
%
%      H = TAGTRAVEL2 returns the handle to a new TAGTRAVEL2 or the handle to
%      the existing singleton*.
%
%      TAGTRAVEL2('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in TAGTRAVEL2.M with the given input arguments.
%
%      TAGTRAVEL2('Property','Value',...) creates a new TAGTRAVEL2 or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before tagtravel2_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to tagtravel2_OpeningFcn via varargin.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help tagtravel2

% Last Modified by GUIDE v2.5 04-Nov-2008 14:54:05

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @tagtravel2_OpeningFcn, ...
                   'gui_OutputFcn',  @tagtravel2_OutputFcn, ...
                   'gui_LayoutFcn',  [] , ...
                   'gui_Callback',   []);
if nargin && ischar(varargin{1})
    gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
    gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before tagtravel2 is made visible.
function tagtravel2_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   command line arguments to tagtravel2 (see VARARGIN)

% Choose default command line output for tagtravel2
handles.output = hObject;

% Update handles structure
guidata(hObject, handles);

% UIWAIT makes tagtravel2 wait for user response (see UIRESUME)
% uiwait(handles.figure1);


% --- Outputs from this function are returned to the command line.
function varargout = tagtravel2_OutputFcn(hObject, eventdata, handles) 
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;


% --- Executes on button press in giddyup.
function giddyup_Callback(hObject, eventdata, handles)
% hObject    handle to giddyup (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

handles.start=str2num(get(handles.start_date,'string'));
handles.d=str2num(get(handles.num_days,'string'));
handles.window=str2num(get(handles.win_size,'string'));
handles.t=str2num(get(handles.refresh,'string'));
handles.tst=str2num(get(handles.tstep,'string'));

try
if isnan(handles.start)
   msgbox('You have not entered a valid number into the start date field.','Oops...','error')
   return;
elseif mod(handles.start,1)~=0
    msgbox('Start date must be a round number. Start date will be rounded down.','Oops...','warn')
    handles.start = handles.start-mod(handles.iters,1); %ceilings to the nearest 1000
end

end
%enter more of this error catching stuff later...

%__________________________________________________________________________
%%READ THESE IN...
%but first check to make sure they exist

if exist('moondat.txt','file')==0
    msgbox('The file "moondate.txt" is missing from the MixSIR directory!')
    return;
else
    load moondat.txt;    %creates a matrix called moondat   
end
if exist('tagdat.txt','file')==0
    msgbox('The file "tagdat.txt" is missing from the MixSIR directory!')
    return;
else
    load tagdat.txt; data=tagdat;		%creates a matrix called tagdat
end
if handles.draw_map == 1
        if exist('map.txt','file')==0
            msgbox('The file "map.txt" is missing from the MixSIR directory!')
        return;
        else
            load map.txt;    %creates a matrix called coast
            coast=map;
        end
end
%__________________________________________________________________________
% NOW PUT THE DATA INTO THE CORRECT FORMAT FOR WORK...

data = sortrows(data,4); %sort the data based on time

% pull a subset of data for processing based on the start date and the
% number of days to run based on user input
data = data(find(data(:,4)>handles.start & data(:,4)<handles.start+handles.d),:);

% get moon data only from dates near those used in the animation
moondat = moondat(find(moondat(:,1)>floor(min(data(:,4)))-10 & moondat(:,1)<floor(min(data(:,4)))+handles.d+10),:);

%make list of tag IDs
tag = unique(data(:,1));

%make list of hydro IDs
hyd = unique(data(:,2)); %assuming all hydros have unique Latitudes (BXS: this is dangerous... fix it)

n_tag = size(tag,1); %number of tags
n_hyd = size(hyd,1); %number of hydrophones

%CALCULATE AN APPROPRIATE FIGURE AXES BASED ON HYDRO POSITIONS
%furthest latitude distance (in lat coords) between hydrophones
 latsize=abs(min(data(:,2))-max(data(:,2)));
 %furthest longitude distance (in lat coords) between hydrophones
 lonsize=abs(min(data(:,3))-max(data(:,3)));
 %make axes 10% wider on either size of the min convex poly of hydrophones
 %(in both long and lat directions)
 axis_bounds= [min(data(:,3))-lonsize*.1 max(data(:,3))+lonsize*.1 ...
                min(data(:,2))-latsize*.1 max(data(:,2))+latsize*.1];
  %set date text position in top left of animation figure
 text_pos = [min(data(:,3))-lonsize*.05 max(data(:,2))+latsize*.05];

 
% NOW TURN ALL TIME-DATE RECORDS INTO HOUR STEPS
%the following line floors to the last hour (eg- 4:15pm, becomes 4pm)

tsteps = (floor(max(data(:,4)))-floor(min(data(:,4)))) / (handles.tst/24);%how many time steps?
tvec=(floor(min(data(:,4))):(handles.tst/24):floor(max(data(:,4))));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    map=colormap;
a1=(handles.moon_plot); %use subplot because it automatically writes over old plot in uipanel
a2=(handles.animate_plot);
    
for i=handles.window+1:tsteps-handles.window %loop over time steps (accounts for window)
    tic;
    
    %update (or initiate) the lunar phase plot-----------------------------
    cla(a1);
    plot(a1,moondat(:,1),moondat(:,2),'LineWidth',2) %plots lunar cycle as line
    set(a1,'XTick',[]) %clears the X axis tick values (serial dates not useful)
    hold (a1); %holds a1 axis
    plot(a1,[floor(tvec(i));floor(tvec(i))],[0;1],'r','LineWidth',4) %plots current date as red bar
    set(a1,'XTick',[]) %clears the X axis tick values (serial dates not useful)
    hold(a1); %unholds a1 axis
    
    
    %for each time step, make a map frame----------------------------------
    if handles.draw_map == 1
    cla(a2);
    plot(a2,coast(:,1),coast(:,2));
    axis(axis_bounds)
    hold on (a2);
    end

    
    for j = 1:n_tag
        clear temper & plt;
        
        %first identify each of the hydros the fish was heard at during the
        %time window
        temper = find(data(:,1) == tag(j,1) & data(:,4)>tvec(i-handles.window) & data(:,4)<tvec(i+handles.window));
        %now put the hydro lat longs into rows for plotting
        if isempty(temper)== 0 %only do if fish was heard
            tempdat =  data(temper,:);
            [B I J]=unique(tempdat(:,2),'first');
            I=sort(I);
            plt(:,1) = tempdat(I,3); 
            plt(:,2) = tempdat(I,2); 
            %now add 'jitter' by choosing random number between +-2%
            %of the lat and long window size
            plt(:,1) = plt(:,1) + ((rand(size(plt,1),1)-.5 ) *lonsize*1.2*.02); %1.2 is to add the 20% past the hydro MCP range
            plt(:,2) = plt(:,2) + ((rand(size(plt,1),1)-.5 ) *latsize*1.2*.02); %1.2 is to add the 20% past the hydro MCP range
            %now plot it
            plot(a2,plt(:,1),plt(:,2),'-rs','LineWidth',1,...
                'MarkerEdgeColor','k',...
                'MarkerFaceColor',map(j,:),...
                'MarkerSize',6);
            axis(axis_bounds); % set axis size
                
            hold on (a2);
            %MAKE LINES AND ARROW HEADS--------------------------------------------
            if size(plt,1)>1 %make an arrow end for each movement
                %this draws a tiny line back along the path so that only
                %the arrow head shows. To do this, I first needed to calc
                %the slope of the line.
                steps = size(plt,1);
                slope = (plt(steps,2) - plt(steps-1,2)) /(plt(steps,1) - plt(steps-1,1)); %rise over run.
                if plt(steps,1) < plt(steps-1,1)
                    arrow([plt(steps,1)+.000001,plt(steps,2)+slope*.000001],[plt(steps,1),plt(steps,2)],'FaceColor','r','Length',20);
                else
                    arrow([plt(steps,1)-.000001,plt(steps,2)+slope*-.000001],[plt(steps,1),plt(steps,2)],'FaceColor','r','Length',20);    
                end
            end
            %END ARROWHEADS----------------------------------------------------
  
        end
        
    end
        %MAKE DAY AND NIGHTS----------------------------------------------------
    %gray shading to axis area based on daylight
    if mod(tvec(i),1) < (5/24) || mod(tvec(i),1) > (22/24)
        whitebg(gcf,[.8,.8,.8]); %night hours are between 10pm and 5am
    elseif mod(tvec(i),1) < (7/24) || mod(tvec(i),1) > (21/24)
        whitebg(gcf,[.9,.9,.9]); %dusk hours are 8-10pm and 5-7am
    else
        whitebg(gcf,[1,1,1]); %daylight the rest of the time
    end
    %END DAYS ANDNIGHTS----------------------------------------------------

    text(text_pos(1),text_pos(2),datestr(tvec(i))); %puts the date and hour on the plot
    M(1)=getframe; %for now, this just refreshs screen
    while toc<handles.t;
        buy_time=1; %this is just a cpu eater to keep each loop step exactly t amount of time
    end
    hold off; %end loop, so allow a new fig to be built
end

% --- Executes on button press in draw_map.
function draw_map_Callback(hObject, eventdata, handles)
% hObject    handle to draw_map (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
handles.draw_map = get(hObject,'Value');
% Hint: get(hObject,'Value') returns toggle state of draw_map
guidata(hObject,handles);



function start_date_Callback(hObject, eventdata, handles)
% hObject    handle to start_date (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of start_date as text
%        str2double(get(hObject,'String')) returns contents of start_date as a double


% --- Executes during object creation, after setting all properties.
function start_date_CreateFcn(hObject, eventdata, handles)
% hObject    handle to start_date (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end



function num_days_Callback(hObject, eventdata, handles)
% hObject    handle to num_days (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of num_days as text
%        str2double(get(hObject,'String')) returns contents of num_days as a double


% --- Executes during object creation, after setting all properties.
function num_days_CreateFcn(hObject, eventdata, handles)
% hObject    handle to num_days (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end



function win_size_Callback(hObject, eventdata, handles)
% hObject    handle to win_size (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of win_size as text
%        str2double(get(hObject,'String')) returns contents of win_size as a double


% --- Executes during object creation, after setting all properties.
function win_size_CreateFcn(hObject, eventdata, handles)
% hObject    handle to win_size (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end



function refresh_Callback(hObject, eventdata, handles)
% hObject    handle to refresh (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of refresh as text
%        str2double(get(hObject,'String')) returns contents of refresh as a double


% --- Executes during object creation, after setting all properties.
function refresh_CreateFcn(hObject, eventdata, handles)
% hObject    handle to refresh (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end



function tstep_Callback(hObject, eventdata, handles)
% hObject    handle to tstep (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of tstep as text
%        str2double(get(hObject,'String')) returns contents of tstep as a double


% --- Executes during object creation, after setting all properties.
function tstep_CreateFcn(hObject, eventdata, handles)
% hObject    handle to tstep (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end