trackplot.m 25.07 KiB
function [tout, Nout] = trackplot(varargin)
%TRACKPLOT Estimates the line of synchronization in a cross recurrence plot.
% TRACKPLOT(X [,DX, DY, FP, PARAM]) estimates the line of synchronization in a
% cross recurrence plot X. The resulted path can be saved to the
% workspace variable T_OUT. This command allows interactive
% changing of estimation parameters.
%
% [A B]=TRACKPLOT(X,DX,DY) estimates the LOS in a recurrence
% plot X and stores it in A. The number of recurrence points
% met by the LOS is stored in B(1) and the number of lacks in
% the LOS is stored in B(2).
%
% TRACKPLOT(X,DX,DY,FP) estimates a LOS whith some fixed points
% given by the two-column vector FP (eg. [194 0; 201 10]).
%
% Parameter: The search of the LOS can be forced with the
% parameters DX and DY. PARAM can be used to suppress the GUI
% (useful in order to use this programme by other programmes).
%
% Suppressing the GUI.
% gui - Creates the GUI and the output plot.
% nogui - Suppresses the GUI and the output plot.
% silent - Suppresses all output.
%
% Examples: y = sin([1:900]*2*pi/67)';
% y2 = sin(.01*([1:900]*2*pi/67).^2)';
% x = crp_big(y,y2,3,12,.1,'fan','nogui');
% trackplot(x,2,2)
%
% See also CRP2, CRP and CRP_BIG.
%
% References:
% Marwan, N., Thiel, M., Nowaczyk, N.: Cross Recurrence Plot Based
% Synchronization of Time Series, Nonlin. Proc. Geophys., 2001.
% Copyright (c) 1998-2005 by AMRON
% Norbert Marwan, Falko Zetsche, Potsdam University, Germany
% http://www.agnld.uni-potsdam.de
%
% $Date$
% $Revision$
%
% $Log$
% Revision 4.9 2005/09/02 08:02:57 marwan
% line fitting algorithm improved (linear interpolation between set points)
%
% Revision 4.8 2005/03/16 11:19:02 marwan
% help text modified
%
% Revision 4.7 2004/11/10 07:09:31 marwan
% initial import
%
%
% This program is part of the new generation XXII series.
%
% This program is free software; you can redistribute it and/or
% modify it under the terms of the GNU General Public License
% as published by the Free Software Foundation; either version 2
% of the License, or any later version.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% programme properties
global errcode props
init_properties
errcode=0; mflag=0;fixp=[]; X=[];Dmax1=1;Dmax2=1;
set(0,'ShowHidden','On')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% check and read the input
warning off
error(nargchk(1,5,nargin)) % error, if < 1 oder > 4 input variables
if nargout>2, error('Too many output arguments'), end
if isnumeric(varargin{1}) % read commandline input
vargin{4}=[];
i_double=find(cellfun('isclass',varargin,'double'));
i_char=find(cellfun('isclass',varargin,'char'));
% check the text input parameter for gui
check_gui={'gui','nog','sil'};
temp_gui=0;
if ~isempty(i_char)
for i=1:length(i_char),
varargin{i_char(i)}(4)='0';
temp_gui=temp_gui+strcmpi(varargin{i_char(i)}(1:3),check_gui');
end
nogui=min(find(temp_gui))-1;
if isempty(nogui), nogui=0; end
if nogui>2,nogui=2; end
else
nogui=0;
end
% get the parameters for creating RP
if max(size(varargin{1}))<2
error('Too less values in data X.')
end
X=double(varargin{1});
Nx = size(X,2); Ny = size(X,1);
if nargin<3 ,Dmax1=1; Dmax2=1;
else
if isnumeric(varargin{2}), Dmax1=double(varargin{2});
else Dmax1=1;
end
if isnumeric(varargin{3}), Dmax2=double(varargin{3});
else Dmax2=1;
end
end
if nargin>=4
if isnumeric(varargin{4}), fixp=fliplr(varargin{4});
fixp = sortrows(fixp,1);
else fixp=[];
end
end
action='start_gui';
else
%%%%%%%%%%%%%%%%%%%%%%%%%% read input from the GUI
action=varargin{1};
nogui=0;
h=get(gcf,'Name');
h=h(findstr(h,'(')+1:findstr(h,')')-1);
hTP=findobj('Name',['TrackPlot (' h ')']);
hTPCtrl=findobj('Name',['TPControl (' h ')']);
h=str2num(h);
temp=get(findobj('Tag','fixp','Parent',findobj('Parent',hTP,'Tag','TPPlot')),'XData');
if ~isempty(temp)
for i=1:length(temp), if iscell(temp), fixp(i,1)=temp{i}; else, fixp(i,1)=temp(i); end, end
temp=get(findobj('Tag','fixp','Parent',findobj('Parent',hTP,'Tag','TPPlot')),'YData');
for i=1:length(temp), if iscell(temp), fixp(i,2)=temp{i}; else, fixp(i,2)=temp(i); end, end
fixp = sortrows(fixp,1);
else fixp=[];
end
if ( nogui == 0 & ~isempty(hTPCtrl))
Dmax1=str2num(get(findobj('Tag','LOSwidthX','Parent',hTPCtrl),'string'));
Dmax2=str2num(get(findobj('Tag','LOSwidthY','Parent',hTPCtrl),'string'));
X=double(get(findobj('Tag','RP','Parent',findobj('Parent',hTP,'Tag','TPPlot')),'UserData'));
Nx=size(X,2); Ny=size(X,1);
end
if isempty(findobj('Tag','TPFig')) & nogui==0
action='start_gui';
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% splash the GPL
splash_gpl('trackplot');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% nogui
if nogui>0
hTP=9999;
action='LOSsearch';
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% switch routines
try
switch(action)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% start_gui
case 'start_gui'
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% create GUI
errcode=2;
if isempty(X)
lasterr(['??? Error using ==> trackplot',10,...
'No properly input. Check the class of input.',10,...
'Only INT and DOUBLE are allowed, but not LOGICAL.'])
disp(lasterr)
return
end
scr=get(0,'ScreenSize');
root_ud=get(0,'UserData');
if isstruct(root_ud)
if isfield(root_ud,'tp')
if ~isempty(root_ud.tp)
root_ud.tp=[root_ud.tp max(root_ud.tp)+1];
else
root_ud.tp=1;
end
h=num2str(root_ud.tp(end));
else
root_ud.tp=1;
h=num2str(1);
end
else
root_ud.old=root_ud;
root_ud.tp=1;
h=num2str(1);
end
set(0,'UserData',root_ud)
%%%%%%%%%%%%%%%%% Trackplot Figure
h8=figure('Tag','TPFig',... %Plot Figure
'Position',[scr(3)/4 scr(4)/4 3*scr(3)/8 3*scr(3)/8 ],...
'NumberTitle','off',...
'Name',['TrackPlot (' h ')'],...
'Color',props.window.Color,...
'DeleteFcn','trackplot smartclose',...
'PaperType','a4',...
'PaperOrientation','portrait');
set(h8,props.window,'Units','Norm')
h1=axes(props.axes,'Parent',h8, ... % Initialize CRP Plot
'Units','normalized', ...
'Position',[.1 .1 .8 .8], ...
'Color',[1 1 1], ...
'Tag','TPPlot', ...
'XColor',[0 0 0], ...
'YColor',[0 0 0], ...
'ZColor',[0 0 0]);
h2=imagesc('Tag','RP','Parent',h1); % Plot Data
X=(X-min(min(X)))/(max(max(X))-min(min(X)));
% set(h2,'cdata',64*(-(X-max(max(X)))),'UserData',X);
axis tight
minvalue=min(min(X));maxvalue=max(max(X));
set(h1,'tickdir','out','box','on','layer','top')
h4=line('Parent',h1,'visible','off','Tag','Diagonal',...
'color',[1 0 0],'LineWidth',1);
cm={'hsv';'hot';'gray';'french';'bone';'copper';... % Colormap
'pink';'flag';'lines';'colorcube';...
'jet';'prism';'cool';'autumn';...
'spring';'winter';'summer'};
h0=uimenu('Label','Colormap','Tag','cm');
if (length(find(X==minvalue))+length(find(X==maxvalue))==length(X(:)))
set(h2,'cdata',(-X),'UserData',X);
colormap((gray(2))); cmflag=0;
else
set(h2,'cdata',(X),'UserData',X);
colormap(french(256)); cmflag=1;
end
set(h1,'XLim',[0 size(X,2)], 'YLim',[0 size(X,1)])
for i=1:length(cm);
h1=uimenu(h0,'Label',cm{i},'Checked','Off',...
'Tag',num2str(i),...
'Callback','trackplot colormap');
if i==18 & cmflag==1, set(h1,'Checked','On'), end
end
h1=uimenu(h0,'Label','b/w','Tag','18','Callback','trackplot colormap');
if cmflag==0, set(h1,'Checked','On'), end
h1=uimenu(h0,'Label','inverse','Tag','19','Separator','On','Callback','trackplot colormap');
h1=uimenu('Label','SmartClose',... % SmartClose
'Callback','trackplot smartclose');
%%%%%%%%%%%%%%%%% Control Figure
errcode=3;
h3=figure('Tag','TPFig',... % Control Figure
'Position',[5*scr(3)/8+10 scr(4)/4 1*scr(3)/6 3*scr(4)/9.8],...
'NumberTitle','off',...
'Color',props.window.Color,...
'Name',['TPControl (' h ')'],...
'DeleteFcn','trackplot handlevisON',...
'MenuBar','None',...
'Resize','Off');
set(h3,props.window,'Units','Norm')
uicontrol(props.frame, ... % Frame LOSsearch
'Units','Normalized',...
'Position',[.1 .51 .8 .42]);
h0=uicontrol(props.text,... % Text LOSsearch
'Units','Normalized',...
'FontAngle','italic', ...
'String','LOS Search',...
'Tag','TextLOSsearch',...
'Enable','on',...
'Position',[.13 .846 .35 .065]);
h1=get(h0,'Extent'); h2=get(h0,'Position'); set(h0,'Position',[h2(1) h2(2) h1(3) h1(4)])
uicontrol(props.button,... % Set Point
'Units','Normalized',...
'Tag','SetPoint',...
'String','Set Point',...
'Position',[.16 .765 .314 .08],...
'Enable','on',...
'ToolTip','Set a Point on LOS.',...
'CallBack','trackplot LOSset');
uicontrol(props.button,... % Clear Point
'Units','Normalized',...
'Tag','ClearPoint',...
'String','Clear P',...
'Position',[.52 .765 .314 .08],...
'Enable','on',...
'ToolTip','Clear a Point on LOS.',...
'CallBack','trackplot LOSclear');
uicontrol(props.button,... % Clear All Point
'Units','Normalized',...
'Tag','ClearAllPoint',...
'String','Clear All',...
'Position',[.52 .65 .314 .08],...
'Enable','on',...
'ToolTip','Clear All Points on LOS.',...
'CallBack','trackplot LOSallclear');
h0=uicontrol(props.text,... % Text LOSwidthX
'Units','Normalized',...
'Tag','LOSwidthXtext',...
'String','dx:',...
'Position',[.16 .542 .35 .07],...
'Enable','on');
h1=get(h0,'Extent'); h2=get(h0,'Position'); set(h0,'Position',[h2(1) h2(2) h1(3) h1(4)])
h0=uicontrol(props.edit,... % Input LOSwidthX
'Units','Normalized',...
'Tag','LOSwidthX',...
'Position',[.31 .55 .14 .06],...
'String',num2str(Dmax1),...
'Enable','on',...
'ToolTip','Insert the LOS search width in X-direction.' );
h1=get(h0,'Extent'); h2=get(h0,'Position'); set(h0,'Position',[h2(1) h2(2) h2(3) h1(4)])
h0=uicontrol(props.text,... % Text LOSwidthY
'Units','Normalized',...
'Tag','LOSwidthYtext',...
'String','dy:',...
'Position',[.52 .542 .35 .07],...
'Enable','on');
h1=get(h0,'Extent'); h2=get(h0,'Position'); set(h0,'Position',[h2(1) h2(2) h1(3) h1(4)])
h0=uicontrol(props.edit,... % Input LOSwidthY
'Units','Normalized',...
'Tag','LOSwidthY',... 'Position',[.68 .55 .14 .06],...
'String',num2str(Dmax2),...
'Enable','on',...
'ToolTip','Insert the LOS search width in Y-direction.' );
h1=get(h0,'Extent'); h2=get(h0,'Position'); set(h0,'Position',[h2(1) h2(2) h2(3) h1(4)])
if ~isunix
h0=uicontrol(props.frame, ... % Frame Embedding
'Units','Normalized',...
'Position',[.1 .365 .38 .11]);
end
uicontrol(props.checkbox, ... % Button Stretch Plot
'Units','Normalized',...
'String','Stretch', ...
'Position',[.1 .365 .38 .11], ...
'Tag','Stretch',...
'CallBack','trackplot stretch',...
'Value',1,...
'ToolTip','Streches the plotted CRP to a squared plot.' );
uicontrol(props.button,... % Button Stores LOS
'Units','Normalized',...
'String','Store',...
'Position',[.52 .365 .38 .11],...
'Tag','Store2',...
'Callback','trackplot store2',...
'Enable','off',...
'ToolTip','Stores the LOS.');
uicontrol(props.button,... % Button ApplyLOSsearch
'Units','Normalized',...
'String','Find LOS',...
'Position',[.1 .22 .8 .11],...
'Tag','Apply2',...
'Callback','trackplot LOSsearch',...
'Enable','on',...
'ToolTip','Searches the LOS.');
uicontrol(props.button,... % Button Help
'Units','Normalized',...
'String','Help',...
'Position',[.1 .075 .38 .11],...
'Tag','Help',...
'Callback','helpwin trackplot',...
'ToolTip','Opens the helpwindow.');
uicontrol('Parent',h3, ...
props.button,... % Button Close
'Units','Normalized',...
'String','Close',...
'Position',[.52 .075 .38 .11],...
'Tag','Close',...
'Callback','trackplot close',...
'ToolTip','Closes Trackplot windows.');
set(h8, 'HandleVis','CallBack')
set(h3, 'HandleVis','CallBack')
if nargout>=1
tout=h0;
end
if nargout>=2
Nout=h2;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% change colormap
case 'colormap'
errcode=81;
c=str2num(get(gcbo,'Tag'));
if c~=19,
set(get(get(gcbo,'Parent'),'Children'),'Checked','Off')
set(gcbo,'Checked','On')
end
cm_old=get(hTP,'Colormap');
cm=[{hsv(256)}; {hot(256)}; {gray(256)};...
{french(256)}; {bone(256)}; {copper(256)}; {pink(256)};...
{flag(256)}; {lines(256)}; {colorcube(256)};...
{jet(256)}; {prism(256)}; {cool(256)};...
{autumn(256)}; {spring(256)}; {winter(256)};...
{summer(256)}; {(gray(2))}; {flipud(cm_old)}];
set(hTP,'Colormap',(cm{c}))
clear cm cm_old c
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% handlevisON
case 'handlevisON'
set(hTP, 'HandleVis','on')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% close
case 'close'
errcode=101;
if ~isempty(findobj('Tag','TPFig')), delete(findobj('Tag','TPFig')), end
root_ud=get(0,'UserData');
if isstruct(root_ud)
if isfield(root_ud,'tp')
root_ud=rmfield(root_ud,'tp');
if length(fieldnames(root_ud))==1
if isfield(root_ud,'old'); root_ud=root_ud.old; end
end
end
end
try, set(0,'UserData',root_ud,props.root), end
clear all
disp('Thank you for using CRP toolbox.')
case 'smartclose'
errcode=102;
[h h1]=strtok(get(hTP,'Name'),'(');
h1([1, end])=[];
if ishandle(hTPCtrl), delete(hTPCtrl), end
if ishandle(hTP), delete(hTP), end
root_ud=get(0,'UserData');
if isstruct(root_ud)
if isfield(root_ud,'tp')
root_ud.tp(root_ud.tp==str2num(h1))=[];
if isempty(root_ud.tp), root_ud=rmfield(root_ud,'tp'); end
end
if length(fieldnames(root_ud))==1
if isfield(root_ud,'old'); root_ud=root_ud.old; end
end
end
try, set(0,'UserData',root_ud,props.root), end
clear all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% store2
case 'store2'
errcode=15;
t=get(findobj('Tag','Apply2','Parent',hTPCtrl),'UserData');
if isempty(t)
warndlg('The LOS vector is still empty. Please start the computation of the LOS before storing.','No LOS')
waitforbuttonpress
else
assignin('base','t_out', t)
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% stretch
case 'stretch'
errcode=7;
if get(findobj('Tag','Stretch','Parent',gcf),'value')==1
set(findobj('Tag','TPPlot','Parent',hTP),'PlotBoxAspectRatio',[1 1 1])
elseif get(findobj('Tag','Stretch','Parent',gcf),'value')==0
set(findobj('Tag','TPPlot','Parent',hTP),'PlotBoxAspectRatio',[Nx Ny 1])
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% LOSmove
case 'LOSmove'
errcode=16;
if isempty(get(gco,'UserData'))
set(gco,'UserData',1,'ButtonDownFcn','trackplot LOSmove_end')
set(gcf,'WindowButtonMotionFcn','trackplot LOSmove')
end
h1 = round(get(gca,'CurrentPoint'));
set(gco, 'XData', h1(1,1), 'YData', h1(1,2))
clear h1
case 'LOSmove_end'
if ~isempty(get(gco,'UserData'))
set(gco,'UserData',[],'ButtonDownFcn','')
set(gcf,'WindowButtonMotionFcn','')
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% LOSclear
case 'LOSclear'
errcode=17;
if gcf==hTPCtrl
figure(hTP)
k=waitforbuttonpress;
h1 = round(get(gca,'CurrentPoint')); h1(2,:)=[]; h1(:,3)=[];
finalRect = rbbox;
h2 = round(get(gca,'CurrentPoint')); h2(2,:)=[]; h2(:,3)=[];
h(1,1)=min(h1(:,1),h2(:,1));h(2,1)=max(h1(:,1),h2(:,1));
h(1,2)=min(h1(:,2),h2(:,2));h(2,2)=max(h1(:,2),h2(:,2));
i=find(fixp(:,1)>=h(1,1) & fixp(:,2)>=h(1,2) & fixp(:,1)<=h(2,1) & fixp(:,2)<=h(2,2));
for j=1:length(i); delete(findobj('tag','fixp','Parent',findobj('Parent',hTP,'Tag','TPPlot'),'xdata',fixp(i(j),1))), end
else
h(1,1)=get(gco,'XData');h(1,2)=get(gco,'YData');
delete(gco)
end
clear h
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% LOSallclear
case 'LOSallclear'
errcode=171;
delete(findobj('tag','fixp','Parent',findobj('Parent',hTP,'Tag','TPPlot')))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% LOSset
case 'LOSset'
errcode=18;
if gcf==hTPCtrl
figure(hTP)
ginput(1);
end
h=round(get(gca,'currentp'));
fixp(end+1,1)=h(1,1); fixp(end,2)=h(1,2);
[i j]=sort(fixp(:,1));
fixp=fixp(j,:);
h0=uicontextmenu;
line(h(1,1),h(1,2),1000,...
'MarkerSize',12,...
'Marker','.',...
'Color',[1 0 0],...
'Tag','fixp',...
'UIContextMenu',h0)
uimenu(h0, 'Label', 'Set Point', 'Callback', 'trackplot LOSset')
uimenu(h0, 'Label', 'Move Point', 'Callback', 'trackplot LOSmove')
uimenu(h0, 'Label', 'Clear Point', 'Callback', 'trackplot LOSclear')
clear h h0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% LOSsearch
case 'LOSsearch'
errcode=19;
x0=0;
y0=0;
N=size(X);
Nleer=0; Nvoll=0;
minvalue=min(min(X));maxvalue=max(max(X));
if (length(find(X==minvalue))+length(find(X==maxvalue))==length(X(:)))
flagpoint=0;
% supress output for 'silence'
if nogui~=2
h_wait=waitbar(0,'Compute LOS ...');
set(h_wait,'HandleVisibility','on',props.window);
end
if nogui==0
setptr([hTP,hTPCtrl],'watch')
h=findobj('Tag','Apply2','Parent',hTPCtrl);
obj.children=get(hTPCtrl,'Children');
obj.enable=get(obj.children,'Enable');
set(obj.children,'Enable','off')
set(h(1),'String','Stop',...
'ToolTip','Stops the computation.',...
'Enable','on',...
'Callback','set(gcbo,''String'',''Stopped'')')
end
% looks for the beginning of the diagonal LOS
errcode=191;
for i=1:N(2);
if nogui==0
if check_stop_LOS(hTP,hTPCtrl,nogui,obj), break, end
end
if nogui~=2
waitbar(i/N(2)) end
if i<=N(1)
if ~isempty(fixp), if i>=fixp(1,1), x0=fixp(1,1); y0=fixp(1,2); end, end
if X(i+y0,1+x0)~=0
v=i+y0;
h=1+x0;
break
end
end
if X(1+y0,i+x0)~=0
h=i+x0;
v=1+y0;
break
end
end
if ~isempty(fixp) & nogui~=0
h=fixp(1,1);v=fixp(1,2);x0=0;y0=0;
end
Nleer=i-1;
% t(1:v)=h;
t(1+x0:h)=v;
% start estimation of the LOS
errcode=192;
mflag=0;
i=2;
while h<N(2)-1 & v<N(1)-1 & mflag~=1,
if nogui==0
if check_stop_LOS(hTP,hTPCtrl,nogui,obj), break, end
end
if nogui~=2
waitbar(h/N(2))
end
dw=1;
if v < 1, v = 1; end
if h < 1, h = 1; end
W=X(v:v+dw,h:h+dw);
W(1,1)=0;
if ~isempty(fixp)
if find(h==fixp(:,1))
dv=fixp(find(h==fixp(:,1)),2)-v;
dh=1;
if dv <= 0
h0 = min(find(t > v+dv));
if isempty(h0), h0 = h; end
t0 = t(end) + dv;
fixold = find(h==fixp(:,1));
if ~isempty(fixold) & fixold > 1; fixold = fixp(fixold - 1,1); else fixold = 0; end
h1 = max([h0-min([Dmax1,Dmax2]);1;fixold]);
h2 = max([h0;1;fixold]);
% h1 = max([h-2*(h-h0);1])
if length(t) > h1 & length(t) > h2
t1 = t(h1);
t2 = t(h2);
p = polyfit([h1 h2 h],[t1 t2 t0],1);
t(h1:h) = polyval(p,h1:h);
% t(h1:h) = spline([h1 h2 h],[t1 t2 t0],h1:h);
end
v = v + dv;
dv = 1; dw = 0;
end
flagpoint=1;
end end
if flagpoint==0;
% looks for the existence of the next recurrence point in diagonal direction
errcode=193;
while sum(sum(W))==0 & mflag==0 & flagpoint==0,
Nleer=Nleer+1;
if ~isempty(fixp)
if find(h+dw==fixp(:,1))
W=1; dh=dw;
flagpoint=1;
break
end
end
dw=dw+1;
if v+dw < N(1) & h+dw < N(2)
W=X(v:v+dw,h:h+dw);
W(1,1)=0;
else
mflag=1;
end
end
if mflag==1
break
end
% determines the coordinates of the next recurrence point
errcode=194;
dh0=min(find(sum(W)));
dv0=min(find(sum(W')));
dh=min(find(W(dv0,:)));
dv=min(find(W(:,dh0)));
if dh>dh0, dh=dh0;end
if dv>dv0, dv=dv0;end
% determines the local width of the diagonal LOS
errcode=195;
dh1=dh;
dv1=dv;
% neues Fenster
WL1=Dmax1;
WL2=Dmax2;
if v+dv1-2+WL2>=N(1), WL2=N(1)-(v+dv1-2); end
if h+dh1-2+WL1>=N(2), WL1=N(2)-(h+dh1-2); end
Wn=X(v+dv1-1:v+dv1-2+WL2,h+dh-1:h+dh-2+WL1);
% Schwerpunkt davon ausrechnen
if sum(sum(Wn))~=0, Sh=sum(sum(Wn).*[1:WL1])/sum(sum(Wn));else Sh=WL1/2;end
if sum(sum(Wn'))~=0,Sv=sum(sum(Wn').*[1:WL2])/sum(sum(Wn'));else Sv=WL2/2;end
% neue Dmax berechnen
if Sh>=Sv, Dmax2n=WL1*Sv/Sh; Dmax1n=WL1; end
if Sv>Sh, Dmax1n=WL2*Sh/Sv; Dmax2n=WL2; end
% while X(v+dv1-1,h:dh-1)==1 & v+dv1-1<N(1) &dv1<Dmax2
while X(v+dv1-1,h+dh-1)==1 & v+dv1-1<N(1) &dv1<Dmax2n
dv1=dv1+1;
end
% while X(v:dv-1,h+dh1-1)==1 & h+dh1-1<N(2) & dh1<Dmax1
while X(v+dv-1,h+dh1-1)==1 & h+dh1-1<N(2) & dh1<Dmax1n
if ~isempty(fixp)
if find(h+dh1-1==fixp(:,1)) & find(h+fix((dh1+dh)/2)==fixp(:,1))
dv=fixp(find((h+dh1-1)==fixp(:,1)),2)-v;
if isempty(dv), dv=fixp(find((h+fix((dh1+dh)/2))==fixp(:,1)),2)-v; end
flagpoint=1; break
end
end
dh1=dh1+1;
end
% compute the mean of the diagonal LOS
errcode=196;
if flagpoint==0
dh=fix((dh1+dh)/2);
dv=fix((dv1+dv)/2);
if dh>0
dh=dh-1;
end
if dv>0
dv=dv-1;
end
% dh=dh+dh1;
end
end % flagpoint end
flagpoint=0;
% output
errcode=197;
if dh~=0 & dv~=0
t(h:h+dh)=v:dv/dh:v+dv;
elseif dv~=0
t(h)=v+dv;
elseif dv==0
t(h:h+dh)=v;
end
Nvoll=Nvoll+sqrt(dv^2+dh^2);
if dh==0 & dv==0
dh=1;
end
% moves the startpoint for further looking
h=h+dh;
v=v+dv;
end
else
% DTW algorithm
if nogui~=2
h_wait=waitbar(0,'Compute LOS ...');
set(h_wait,'HandleVisibility','on',props.window);
end
t=1;
i=y0+1; j=x0+1;
while i<N(1)-Dmax1-1 & j<N(2)-Dmax2-1
errcode=197;
if nogui~=2
waitbar(i/N(1)), j0=j;
end
% [temp pos]=min([sum(sum(X(i:i+Dmax1,j+1:j+1+Dmax2))),sum(sum(X(i+1:i+1+Dmax1,j+1:j+1+Dmax2))),sum(sum(X(i+1:i+1+Dmax1,j:j+Dmax2)))]);
[temp pos]=min([(mean(X(i,j+1:j+1+Dmax2))),(mean(diag(X(i+1:i+1+Dmax1,j+1:j+1+Dmax2)))),(mean(X(i+1:i+1+Dmax1,j)))]);
switch(pos)
case 1
j=j+1;
flag1=1;
case 2 i=i+1; j=j+1;
flag1=1;
flag2=1;
case 3
i=i+1;
flag2=1;
end
if ~isempty(fixp) & flag1==1 & flag2==1
errcode=1981;
h=find(fixp(:,1)==j);
if ~isempty(h) & i<max(fixp(h,2));
h2=find(i<=fixp(h,2));
i=fixp(h(h2(1)),2); flag1=0; end
end
if ~isempty(fixp) & flag2==1 & flag1==1
errcode=1982;
h=find(fixp(:,2)==i);
if ~isempty(h) & j<max(fixp(h,1));
h2=find(j<=fixp(h,1));
j=fixp(h(h2(1)),1)
t(j0:j-1)=t(j0);
flag2=0; end
end
t(j)=i;
end
t(find(~t))=1;
Nvoll=NaN; Nleer=NaN;
end
errcode=199;
if nogui~=2
delete(h_wait)
end
if nogui==0
h1=findobj('Parent',hTP,'Tag','TPPlot');
h2=findobj('Tag','Diagonal','Parent',h1);
if isempty(h2)
h2=line('Parent',h1,'visible','on','Tag','Diagonal','color',[1 0 0],...
'LineWidth',1);
end
set(h2,'xdata',[1:length(t)],'ydata',t,'visible','on')
h=findobj('Tag','Apply2','Parent',hTPCtrl);
set(h(1),'String','Apply',...
'ToolTip','Searches the LOS.',...
'Callback','trackplot LOSsearch')
for i=1:length(obj.enable), set(obj.children(i),'Enable',obj.enable{i}); end
set(findobj('Tag','Store2','Parent',hTPCtrl),'Enable','On')
set(findobj('Tag','Apply2','Parent',hTPCtrl),'UserData',t)
end
if nargout~=0, tout=t; end
if nargout==2
tout=t;
Nout(1)=Nvoll;
Nout(2)=Nleer;
end
if nogui==0
setptr([hTP,hTPCtrl],'arrow')
end
end
warning on
try, set(0,props.root), end
set(0,'ShowHidden','Off')
%%%%%%% error handling
%if 0
catch
z=whos;x=lasterr;y=lastwarn;in=varargin{1};
print_error('trackplot',z,x,y,in,[],action)
try, set(0,props.root), end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function out=check_stop_LOS(hCRP,hCtrl,nogui,obj)
global errcode
errcode=errcode+.02;
out=0;
if nogui==0
h=findobj('Tag','Apply2','Parent',hCtrl);
if strcmpi(get(h(1),'String'),'stopped')
set(h(1),'String','Apply',...
'ToolTip','Searches the LOS.',...
'Callback','trackplot LOSsearch')
for i=1:length(obj.enable), set(obj.children(i),'Enable',obj.enable{i}); end
set(findobj('Tag','Status','Parent',findobj('Parent',hCRP,'Tag','CRPPlot')'),'String','Stopped'),drawnow
setptr([hCRP,hCtrl],'arrow')
out=1;
end
end