function [Xpoint,Ypoint,XYhist] = crosshair_subplots(action);
% CROSSHAIR: A gui interface for reading (x,y) values from a plot.
%
% [Xpoint,Ypoint,XYhist] = crosshair_subplots(action);
%
% A set of mouse driven crosshairs is placed on the current axes,
% and displays the current (x,y) values of the line plot. There is an
% option to provide data specific values or interpolated values. The
% resolution of the data specific values depends on both the data
% resolution and the GUI interface (mainly mouse movement resolution).
% The interpolated values appear to provide a more continuous function,
% however they too depend on the GUI interface resolution. There are
% no options for extrapolation.
%
% For multiple traces, plots with the same length(xdata) are
% tracked. Each mouse click returns Xpoint,Ypoint values and selecting
% 'done' will remove the GUI and restore the mouse buttons to previous
% values. Selecting 'exit' will remove the GUI and close the figure.
%
% In this version (Dec 2002), the crosshairs can track multiple subplots
% and there are new options to define functions of X/Y and monitor
% their results. There is also a new STORE button that creates and
% updates an XYhist struct in the base workspace, which contains value
% labels and values. This version has better controls of X/Y movements,
% including better interpolation movement options. This version attempts
% to respond correctly to keyboard entries, including TAB between subplots,
% RETURN to 'save', ESC for 'done' and all the arrow keys to move the
% crosshairs.
%
% Some further help is given in the tool tips of the GUI.
%
% Note: crosshair should always update the Xpoint,Ypoint in the base
% workspace. Here is an example of how to get return values within
% a script/function after pressing the exit button of crosshair:
% function [X,Y] = crosshair_returnXY
% x = [1:10]; y(1,:) = sin(x); y(2,:) = cos(x);
% figure; plot(x,y); crosshair;
% uiwait
% X = evalin('base','Xpoint');
% Y = evalin('base','Ypoint');
% return
% Copy this text to a function .m file and then call it from the
% base workspace with [X,Y] = crosshair_returnXY
%
% Useage: wt=0:0.01:2.5*pi;
% t=wt;
% subplot(2,1,1),plot(t,1.0*sin(wt),t,1.0*sin(wt-110/180*pi),t,1.0*sin(wt-250/180*pi));
% subplot(2,1,2), plot(t,sin(2*wt));
% crosshair_subplots;
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Licence: GNU GPL, no express or implied warranties
% History: 03/96, Richard G. Cobb <cobbr@plk.af.mil>
% 08/01, Darren.Weber@flinders.edu.au
% replaced obsolete 'table1' with 'interp1'; fixed bug
% with number of 'traces'; rationalized calculations into
% a common subfunction for x,y point calc in 'down','up',
% & 'move' button functions; added option to turn on/off
% interpolation and the exit button; simplified updates
% to graphics using global GUI handle structure.
% 11/01, Darren.Weber@flinders.edu.au
% added tooltips for several GUI handles
% added multiple interpolation methods
% added GUI for data matrix indices (given no interpolation)
% added option to select trace nearest to mouse click point
% reversed order of lines in data matrix to be consistent
% with the value returned from the nearest trace subfunction
% create crosshair lines after finding all plot lines to
% avoid confusing them with the plot lines
% 01/02, Darren.Weber@flinders.edu.au
% should now work across multiple plot figures, given
% that all gui handles and data are now stored in the
% plot figure 'userdata' handle.
% added functionality to move smoothly from interpolation
% back to the delimited data via the "next/previous"
% buttons.
% 06/02, Darren.Weber@flinders.edu.au
% learned how to get values back to a script/function with
% evalin command and updated help above.
% 12/02, Darren.Weber@flinders.edu.au
% added Store uicontrol and associated XYhist variable
% updates in base workspace to provide storage of
% consecutive XY values (thanks to C.A.Swenne@lumc.nl &
% H.van_de_Vooren@lumc.nl for their suggestion/assistance).
% added keyboard controls to left/right arrow and return.
% added prev/next X interpolation interval.
% added handling of subplots in crosshair_subplots version.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if ~exist('action','var')
action = 'init';
elseif isempty(action)
action = 'init';
end
XHR_HANDLES = get(gcf,'userdata');
if gca & ~isempty(XHR_HANDLES),
if ~isequal(gca,XHR_HANDLES.axis),
% Ensure we have the right axis data
XHR_HANDLES.axis = gca;
XHR_HANDLES.data = get_data;
end
end
% Check for specific keys and assign reasonable actions
if strcmp(action, 'keypress'),
CC = get(XHR_HANDLES.gui,'CurrentCharacter');
cc = double(CC);
if cc,
switch cc,
case 9, action = 'axes'; % TAB, switch axes, if possible
case 27, action = 'done'; % ESC
case 28, action = 'prevx'; % left
case 29, action = 'nextx'; % right
case 30, action = 'ygt'; % up
case 31, action = 'ylt'; % down
case 13, action = 'store'; % return/enter
otherwise, action = 'up'; % all other keys
end
end
end
action = lower(action);
switch action,
case 'init',
% Paint GUI
XHR_HANDLES = INIT;
% Update and return values
XHR_HANDLES = updateDATA(XHR_HANDLES);
Xpoint = get(XHR_HANDLES.xvalue,'Value');
Ypoint = get(XHR_HANDLES.yvalue,'Value');
%updateXYhistory(Xpoint,Ypoint);
case 'axes',
% TAB between axes in subplots of a figure
ax = findobj(gcf,'type','axes');
if length(ax) > 1,
nax = find(ax == gca);
if nax < length(ax),
axes(ax(nax+1));
else
axes(ax(1));
end
end
XHR_HANDLES.axis = gca;
XHR_HANDLES.data = get_data;
XHR_HANDLES = updateGUI(XHR_HANDLES);
Xpoint = get(XHR_HANDLES.xvalue,'Value');
Ypoint = get(XHR_HANDLES.yvalue,'Value');
% Mouse Click Down
case 'down',
set(XHR_HANDLES.gui,'WindowButtonMotionFcn','crosshair_subplots(''move'');');
set(XHR_HANDLES.gui,'WindowButtonUpFcn','[Xpoint,Ypoint] = crosshair_subplots(''up'');');
XHR_HANDLES = updateDATA(XHR_HANDLES);
Xpoint = get(XHR_HANDLES.xvalue,'Value');
Ypoint = get(XHR_HANDLES.yvalue,'Value');
%updateXYhistory(Xpoint,Ypoint);
% Mouse Drag Motion
case 'move',
XHR_HANDLES = updateDATA(XHR_HANDLES);
Xpoint = get(XHR_HANDLES.xvalue,'Value');
Ypoint = get(XHR_HANDLES.yvalue,'Value');
%updateXYhistory(Xpoint,Ypoint);
% Mouse Click Up
case 'up',
set(XHR_HANDLES.gui,'WindowButtonMotionFcn',' ');
set(XHR_HANDLES.gui,'WindowButtonUpFcn',' ');
XHR_HANDLES = updateDATA(XHR_HANDLES);
Xpoint = get(XHR_HANDLES.xvalue,'Value');
Ypoint = get(XHR_HANDLES.yvalue,'Value');
%updateXYhistory(Xpoint,Ypoint);
% Next or Previous X point
case {'nextx','prevx','changex','nexty','prevy','changey','ylt','ygt'}, % Change X/Y
XHR_HANDLES = moveXY(XHR_HANDLES,action);
Xpoint = get(XHR_HANDLES.xvalue,'Value');
Ypoint = get(XHR_HANDLES.yvalue,'Value');
%updateXYhistory(Xpoint,Ypoint);
% Store XY values into a history array
case 'store',
Xpoint = get(XHR_HANDLES.xvalue,'Value');
Ypoint = get(XHR_HANDLES.yvalue,'Value');
updateXYhistory(XHR_HANDLES);
% Exit crosshairs GUI
case {'done','exit'},
Xpoint = get(XHR_HANDLES.xvalue,'Value');
Ypoint = get(XHR_HANDLES.yvalue,'Value');
%updateXYhistory(Xpoint,Ypoint);
handles = fieldnames(XHR_HANDLES);
for i=1:length(handles),
switch handles{i},
case {'axis','datalines','gui'},
otherwise,
h = getfield(XHR_HANDLES,handles{i});
if ishandle(h), delete(h); end
end
end
if strcmp(action,'exit');
if ishandle(XHR_HANDLES.gui),
close(XHR_HANDLES.gui);
end
else
lines = findobj(gcf,'tag','XHR_XLINE');
for l = 1:length(lines),
line = lines(l);
if ishandle(line), delete(line); end
end
lines = findobj(gcf,'tag','XHR_YLINE');
for l = 1:length(lines),
line = lines(l);
if ishandle(line), delete(line); end
end
set(XHR_HANDLES.gui,'WindowButtonUpFcn','');
set(XHR_HANDLES.gui,'WindowButtonMotionFcn','');
set(XHR_HANDLES.gui,'WindowButtonDownFcn',' ');
refresh(XHR_HANDLES.gui);
end
clear XHR_HANDLES;
return;
end
if ishandle(XHR_HANDLES.axis),
set(XHR_HANDLES.axis,'userdata',XHR_HANDLES.data);
end
if ishandle(XHR_HANDLES.gui),
set(XHR_HANDLES.gui,'userdata',XHR_HANDLES);
figure(XHR_HANDLES.gui);
end
return;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function updateXYhistory(H),
Ch = get(H.yindex,'Value');
X = get(H.xvalue,'Value');
Y = get(H.yvalue,'Value');
fX = get(H.fxvalue,'Value');
fY = get(H.fyvalue,'Value');
fXeq = get(H.fxeq,'String');
fYeq = get(H.fyeq,'String');
XY.labels = {'Channel','X','Y',['f(x) = ',fXeq],['f(y) = ',fYeq]};
XY.data = [Ch,X,Y,fX,fY];
if evalin('base','exist(''XYhist'',''var'');'),
XYhist = evalin('base','XYhist');
% get the current history set
set = getfield(XYhist,['set',num2str(XYhist.set)]);
if isequal(set.labels,XY.labels),
set.data(end+1,:) = XY.data;
XYhist = setfield(XYhist,['set',num2str(XYhist.set)],set);
assignin('base','XYhist',XYhist);
else
fprintf('\nWarning: creating new set of XYhist in base workspace.\n\n');
XYhist.set = XYhist.set + 1;
XYhist = setfield(XYhist,['set',num2str(XYhist.set)],XY);
assignin('base','XYhist',XYhist);
end
else
XYhist.set = 1;
XYhist.set1 = XY;
assignin('base','XYhist',XYhist);
end
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [H] = moveXY(H,move)
interp = get(H.interp,'Value');
xinterp = get(H.xinterp,'Value');
if (xinterp > 0) & (interp > 1),
% Use incremental interpolation of x values
switch move,
case 'nexty', H.data.yindex = H.data.yindex + 1;
case 'prevy', H.data.yindex = H.data.yindex - 1;
case 'ygt',
ydata = interpYall(H);
[ysort,yi] = sort(ydata);
currYI = find(ysort > H.data.ypoint);
if min(currYI),
H.data.yindex = yi(min(currYI));
end
case 'ylt',
ydata = interpYall(H);
[ysort,yi] = sort(ydata);
currYI = find(ysort < H.data.ypoint);
if max(currYI),
H.data.yindex = yi(max(currYI));
end
case 'nextx',
H.data.xpoint = H.data.xpoint + xinterp;
case 'prevx',
H.data.xpoint = H.data.xpoint - xinterp;
end
H = checkdatarange(H);
H = interpY(H);
updateGUI(H);
return
end
% No interpolation of x values...
if (interp > 1)
xdata = H.data.xdata(:,H.data.yindex);
[H.data.xindex] = NearestXYArrayPoint( xdata, H.data.xpoint, move );
end
switch move,
case 'nextx',
% Increase current xindex by one
if(interp == 1), H.data.xindex = H.data.xindex + 1; end
case 'prevx',
% Decrease current xindex by one
if(interp == 1), H.data.xindex = H.data.xindex - 1; end
case 'nexty', H.data.yindex = H.data.yindex + 1;
case 'prevy', H.data.yindex = H.data.yindex - 1;
case 'ygt',
ydata = H.data.ydata(H.data.xindex,:);
[ysort,yi] = sort(ydata);
currYI = find(ysort == H.data.ypoint);
if currYI < length(yi),
H.data.yindex = yi(currYI+1);
end
case 'ylt',
ydata = H.data.ydata(H.data.xindex,:);
[ysort,yi] = sort(ydata);
currYI = find(ysort == H.data.ypoint);
if currYI > 1,
H.data.yindex = yi(currYI-1);
end
otherwise
end
H = checkdatarange(H);
% Get x/y value at new x/y index
H.data.xpoint = H.data.xdata(H.data.xindex,H.data.yindex);
H.data.ypoint = H.data.ydata(H.data.xindex,H.data.yindex);
set(H.interp,'Value',1);
H = updateGUI(H);
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ H ] = checkdatarange(H),
% Ensure that x/y index is within data range
s = size(H.data.xdata,1);
if( H.data.xindex < 1 ),
H.data.xindex = 1;
elseif( H.data.xindex >= s ),
H.data.xindex = s;
end
s = size(H.data.ydata,2);
if( H.data.yindex < 1 ),
H.data.yindex = 1;
elseif( H.data.yindex >= s ),
H.data.yindex = s;
end
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ H ] = updateDATA( H )
% Only update if mouse pointer is in the
% axis limits
set(H.gui,'units','normalized');
axpos = get(H.axis,'position');
figcp = get(H.gui,'Currentpoint');
axlim = axpos(1) + axpos(3);
aylim = axpos(2) + axpos(4);
if or(figcp(1) > (axlim+.01), figcp(1) < (axpos(1)-.01)),
return;
elseif or(figcp(2) > (aylim+.01), figcp(2) < (axpos(2)-.01)),
return;
end
CurrentPoint = get(H.axis,'Currentpoint');
H.data.xpoint = CurrentPoint(1,1);
H.data.ypoint = CurrentPoint(1,2);
if get(H.traceNearest,'Value'),
% Get new yindex for nearest trace
[ H.data.xpoint, ...
H.data.xindex, ...
H.data.ypoint, ...
H.data.yindex ] = NearestXYMatrixPoint( H.data.xdata,...
H.data.ydata,...
H.data.xpoint,...
H.data.ypoint);
end
CurrentPoint = get(H.axis,'Currentpoint');
H.data.xpoint = CurrentPoint(1,1);
H.data.ypoint = CurrentPoint(1,2);
H = interpY(H);
H = updateGUI(H);
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ H ] = updateGUI( H )
InterpMethod = get(H.interp,'Value');
if (InterpMethod > 1)
% There is no specific matrix x-index for
% an interpolated point, but the nearest xindex
% is always returned from the interp function below
% so that next/prev move function works correctly
set(H.xindex,'String','interp');
else
set(H.xindex,'String',num2str(H.data.xindex));
end
set(H.xindex,'Value',H.data.xindex);
tracestr = sprintf('%d',H.data.yindex);
set(H.yindex,'String',tracestr,'Value',uint16(H.data.yindex));
%set(H.trace,'Value',uint16(H.data.yindex));
% Create the crosshair lines on the figure, crossing at the x,y point
x_rng = get(H.axis,'Xlim');
y_rng = get(H.axis,'Ylim');
set(H.data.xline,'Xdata',[H.data.xpoint H.data.xpoint],'Ydata',y_rng);
set(H.data.yline,'Ydata',[H.data.ypoint H.data.ypoint],'Xdata',x_rng);
% Update the x,y values displayed for the x,y point
xstring = sprintf('%14.6f',H.data.xpoint);
ystring = sprintf('%14.6f',H.data.ypoint);
set(H.xvalue,'String',xstring,'Value',H.data.xpoint);
set(H.yvalue,'String',ystring,'Value',H.data.ypoint);
% Calculate the f(x) function
fyeq = get(H.fyeq,'string');
fyval = eval(strrep(fyeq,'y',ystring));
fystr = sprintf('%14.6f',fyval);
set(H.fyvalue,'String',fystr,'Value',fyval);
% Calculate the f(y) function
fxeq = get(H.fxeq,'string');
fxval = eval(strrep(fxeq,'x',xstring));
fxstr = sprintf('%14.6f',fxval);
set(H.fxvalue,'String',fxstr,'Value',fxval);
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ H ] = interpY( H )
% In this function, xdata & ydata are arrays, not matrices
xdata = H.data.xdata(:,H.data.yindex);
ydata = H.data.ydata(:,H.data.yindex);
if H.data.xpoint >= max(xdata)
H.data.xpoint = max(xdata);
H.data.xindex = find(xdata == max(xdata));
H.data.ypoint = ydata(H.data.xindex);
return;
elseif H.data.xpoint <= min(xdata)
H.data.xpoint = min(xdata);
H.data.xindex = find(xdata == min(xdata));
H.data.ypoint = ydata(H.data.xindex);
return;
end
% 'none|nearest|linear|spline|cubic'
interp = get(H.interp,'Value');
switch interp
case 1
% Given that xdata & ydata are the same length arrays,
% we can find the ypoint given the nearest xpoint.
[H.data.xindex, H.data.xpoint] = NearestXYArrayPoint( xdata, H.data.xpoint );
H.data.ypoint = ydata(H.data.xindex);
case 2
H.data.ypoint = interp1( xdata, ydata, H.data.xpoint, 'nearest' );
case 3
H.data.ypoint = interp1( xdata, ydata, H.data.xpoint, 'linear' );
case 4
H.data.ypoint = interp1( xdata, ydata, H.data.xpoint, 'spline' );
case 5
H.data.ypoint = interp1( xdata, ydata, H.data.xpoint, 'cubic' );
otherwise
%use default (linear in matlabR12)
H.data.ypoint = interp1( xdata, ydata, H.data.xpoint );
end
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ Yall ] = interpYall( H )
xdata = H.data.xdata(:,H.data.yindex);
Yall = H.data.ydata;
if H.data.xpoint >= max(xdata),
H.data.xpoint = max(xdata);
H.data.xindex = find(xdata == max(xdata));
Yall = ydata(:,H.data.xindex);
return;
elseif H.data.xpoint <= min(xdata),
H.data.xpoint = min(xdata);
H.data.xindex = find(xdata == min(xdata));
Yall = ydata(:,H.data.xindex);
return;
end
% 'none|nearest|linear|spline|cubic'
interp = get(H.interp,'Value');
switch interp,
case 1
% do nothing in this case
case 2
Yall = interp1( xdata, Yall, H.data.xpoint, 'nearest' );
case 3
Yall = interp1( xdata, Yall, H.data.xpoint, 'linear' );
case 4
Yall = interp1( xdata, Yall, H.data.xpoint, 'spline' );
case 5
Yall = interp1( xdata, Yall, H.data.xpoint, 'cubic' );
otherwise
%use default (linear in matlabR12)
Yall = interp1( xdata, Yall, H.data.xpoint );
end
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ index, point ] = NearestXYArrayPoint( data_array, point, type )
if ~exist('type','var') type = ''; end
% In this function, input data_array is an array, not a matrix.
% This function returns the data point in the array
% that has the closest value to the value given (point). In
% the context of 'crosshair' the point is a mouse position.
if point >= max(data_array)
point = max(data_array);
index = find(data_array == point);
return;
elseif point <= min(data_array)
point = min(data_array);
index = find(data_array == point);
return;
end
data_sorted = sort(data_array);
greater = find(data_sorted > point);
greater_index = greater(1);
lesser = find(data_sorted < point);
lesser_index = lesser(length(lesser));
greater_dif = data_sorted(greater_index) - point;
lesser_dif = point - data_sorted(lesser_index);
if strcmp(type,'nextx'),
index = greater_index;
elseif strcmp(type,'prevx'),
index = lesser_index;
else
if (greater_dif < lesser_dif)
index = find(data_array == data_sorted(greater_index));
else
index = find(data_array == data_sorted(lesser_index));
end
end
point = data_array(index);
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ xpoint, xindex, ypoint, yindex ] = NearestXYMatrixPoint( Xdata, Ydata, xpoint, ypoint )
% In this function, Xdata & Ydata are matrices of the same dimensions.
% This function attempts to find the nearest values in Xdata & Ydata
% to the mouse position (xpoint, ypoint).
% It is assumed that Xdata has identical columns, so we only really
% need the first column to find the nearest value to xpoint.
[ xindex, xpoint ] = NearestXYArrayPoint( Xdata(:,1), xpoint );
% Now, given the xpoint, we can select just that row of the
% Ydata matrix corresponding to the xpoint.
ydata = Ydata(xindex,:);
% The ydata array is searched in same manner as the xdata
% array for the nearest value.
[ yindex, ypoint ] = NearestXYArrayPoint( ydata, ypoint );
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [ data ] = get_data,
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Get Line Data from Current Axes
if gca,
% see if we have stored any data in userdata already
data = get(gca,'userdata');
if isfield(data,'lineobj'),
if data.lineobj,
return;
end
end
data.ylim = get(gca,'ylim');
% Get all the line objects from the current axis
data.lineobj = findobj(gca,'Type','line');
if data.lineobj,
data.lines = length(data.lineobj);
data.xdata = [];
data.ydata = [];
data.xpoint = [];
data.ypoint = [];
data.xindex = 1;
data.yindex = 1;
for i = data.lines:-1:1,
% put line data into a columns of data
data.xdata(:,i) = get(data.lineobj(i),'XData').';
data.ydata(:,i) = get(data.lineobj(i),'YData').';
end
% Set X,Y cross hair lines
% Do this after finding all the line axis children
% to avoid confusing these lines with those of the
% plot itself (counted above).
x_rng = get(gca,'Xlim');
y_rng = get(gca,'Ylim');
data.xline = line(x_rng,[y_rng(1) y_rng(1)]);
data.yline = line(x_rng,[y_rng(1) y_rng(1)]);
set(data.xline,'Color','r','EraseMode','xor','tag','XHR_XLINE');
set(data.yline,'Color','r','EraseMode','xor','tag','XHR_YLINE');
set(gca,'userdata',data);
else
fprintf('No lines in current axes\n');
end
end
return
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [H] = INIT
H.gui = gcf; % Get current figure handles
H.axis = gca; % Get current axis handles
H.data = get_data;
% store current button fcn
H.button = get(H.gui,'WindowButtonDownFcn');
% set XHR button down fcn
set(H.gui,'WindowButtonDownFcn','crosshair_subplots(''down'');');
set(H.gui,'KeyPressFcn','crosshair_subplots(''keypress'');');
% Match background figure colour
bgcolor = get(H.gui,'Color');
% Try to adapt the foreground colour a little
black = find(bgcolor <= .6);
fgcolor = [0 0 0]; %black text
if length(black)>2, fgcolor = [1 1 1]; end
Font.FontName = 'Helvetica';
Font.FontUnits = 'Pixels';
Font.FontSize = 10;
Font.FontWeight = 'normal';
Font.FontAngle = 'normal';
if findobj(gcf,'tag','XHR_YFLIP'),
H.yflip = findobj(gcf,'tag','XHR_YFLIP');
else
H.yflip = uicontrol(H.gui,'tag','XHR_YFLIP','Style','pushbutton',...
'Units','Normalized',Font,...
'Position',[.00 .95 .08 .05],...
'TooltipString','Flip Y Axis', ...
'String','Flip',...
'Callback',strcat('XHR = get(gcf,''userdata'');',...
'ydir = get(XHR.axis,''YDir'');',...
'if isequal(ydir,''normal''),',...
' set(XHR.axis,''YDir'',''reverse'');',...
'else,',...
' set(XHR.axis,''YDir'',''normal'');',...
'end;',...
'set(XHR.gui,''userdata'',XHR); figure(XHR.gui); clear XHR;'));
end
if findobj(gcf,'tag','XHR_YRESET'),
H.yreset = findobj(gcf,'tag','XHR_YRESET');
else
H.yreset = uicontrol(H.gui,'tag','XHR_YRESET','Style','pushbutton',...
'Units','Normalized',Font,...
'Position',[.00 .90 .08 .05],...
'TooltipString','Reset Axis Limits', ...
'String','Reset',...
'Callback',strcat('XHR = get(gcf,''userdata'');',...
'XHR.data.ylim = get(XHR.axis,''ylim'');',...
'set(XHR.ymin,''string'',sprintf(''%7.1f'',XHR.data.ylim(1)));',...
'set(XHR.ymax,''string'',sprintf(''%7.1f'',XHR.data.ylim(2)));',...
'set(XHR.gui,''userdata'',XHR); figure(XHR.gui); clear XHR;'));
end
if findobj(gcf,'tag','XHR_YMIN'),
H.ymin = findobj(gcf,'tag','XHR_YMIN');
else
H.ymin = uicontrol(H.gui,'tag','XHR_YMIN','Style','edit',...
'Units','Normalized',Font,...
'Position',[.00 .85 .08 .05],...
'HorizontalAlign','left',...
'TooltipString','Set Y min', ...
'String',sprintf('%7.1f',H.data.ylim(1)),...
'Callback',strcat('XHR = get(gcf,''userdata'');',...
'ymin = str2num(get(XHR.ymin,''string''));',...
'XHR.data.ylim(1) = ymin;',...
'set(XHR.axis,''ylim'',XHR.data.ylim);',...
'set(XHR.ymin,''string'',sprintf(''%7.1f'',XHR.data.ylim(1)));',...
'set(XHR.gui,''userdata'',XHR); figure(XHR.gui); clear XHR ymin;'));
end
if findobj(gcf,'tag','XHR_YMAX'),
H.ymax = findobj(gcf,'tag','XHR_YMAX');
else
H.ymax = uicontrol(H.gui,'tag','XHR_YMAX','Style','edit',...
'Units','Normalized',Font,...
'Position',[.00 .80 .08 .05],...
'HorizontalAlign','left',...
'TooltipString','Set Y max', ...
'String',sprintf('%7.1f',H.data.ylim(2)),...
'Callback',strcat('XHR = get(gcf,''userdata'');',...
'ymax = str2num(get(XHR.ymax,''string''));',...
'XHR.data.ylim(2) = ymax;',...
'set(XHR.axis,''ylim'',XHR.data.ylim);',...
'set(XHR.ymax,''string'',sprintf(''%7.1f'',XHR.data.ylim(2)));',...
'set(XHR.gui,''userdata'',XHR); figure(XHR.gui); clear XHR ymax;'));
end
if findobj(gcf,'tag','XHR_GRID'),
H.grid = findobj(gcf,'tag','XHR_GRID');
else
H.grid = uicontrol(H.gui,'tag','XHR_GRID','Style','checkbox',...
'Units','Normalized',Font,...
'Position',[.00 .75 .08 .05],...
'BackgroundColor',bgcolor,'ForegroundColor',fgcolor,...
'TooltipString','Toggle plot grid on/off.', ...
'String','grid',...
'Callback',strcat('XHR = get(gcf,''userdata'');',...
'grid(XHR.axis);',...
'set(XHR.gui,''userdata'',XHR); figure(XHR.gui); clear XHR;'));
end
if findobj(gcf,'tag','XHR_XVALUE'),
H.xvalue = findobj(gcf,'tag','XHR_XVALUE');
else
H.xvalue = uicontrol(H.gui,'tag','XHR_XVALUE','Style','edit',...
'Units','Normalized',Font,...
'Position',[.13 .95 .15 .05],...
'BackGroundColor',[ 0 .9 0],'ForeGroundColor',[ 0 0 0],...
'TooltipString','X value (Read Only)',...
'String',' ');
end
if findobj(gcf,'tag','XHR_YVALUE'),
H.yvalue = findobj(gcf,'tag','XHR_YVALUE');
else
H.yvalue = uicontrol(H.gui,'tag','XHR_YVALUE','Style','edit',...
'Units','Normalized',Font,...
'Position',[.28 .95 .15 .05],...
'BackGroundColor',[ 0 0 .9],'ForeGroundColor',[ 1 1 1],...
'TooltipString','Y value (Read Only)',...
'String',' ');
end
if findobj(gcf,'tag','XHR_FXEQ'),
H.fxeq = findobj(gcf,'tag','XHR_FXEQ');
else
H.fxeq = uicontrol(H.gui,'tag','XHR_FXEQ','Style','edit',...
'Units','Normalized',...
'Position',[.45 .95 .10 .05],...
'TooltipString','Define f(x) equation here',...
'BackGroundColor',[ 0 .9 0],...
'ForeGroundColor',[ 0 0 0],'String','x');
end
if findobj(gcf,'tag','XHR_FXVALUE'),
H.fxvalue = findobj(gcf,'tag','XHR_FXVALUE');
else
H.fxvalue = uicontrol(H.gui,'tag','XHR_FXVALUE','Style','edit',...
'Units','Normalized',...
'Position',[.55 .95 .15 .05],...
'TooltipString','f(x) result',...
'BackGroundColor',[ 0 .9 0],...
'ForeGroundColor',[ 0 0 0],'String',' ');
end
if findobj(gcf,'tag','XHR_FYEQ'),
H.fyeq = findobj(gcf,'tag','XHR_FYEQ');
else
H.fyeq = uicontrol(H.gui,'tag','XHR_FYEQ','Style','edit',...
'Units','Normalized',...
'Position',[.70 .95 .10 .05],...
'TooltipString','Define f(y) equation here',...
'BackGroundColor',[ 0 0 .9],...
'ForeGroundColor',[ 1 1 1],'String','y');
end
if findobj(gcf,'tag','XHR_FYVALUE'),
H.fyvalue = findobj(gcf,'tag','XHR_FYVALUE');
else
H.fyvalue = uicontrol(H.gui,'tag','XHR_FYVALUE','Style','edit',...
'Units','Normalized',...
'Position',[.80 .95 .15 .05],...
'TooltipString','f(y) result',...
'BackGroundColor',[ 0 0 .9],...
'ForeGroundColor',[ 1 1 1],'String',' ');
end
if findobj(gcf,'tag','XHR_YINDEX'),
H.yindex = findobj(gcf,'tag','XHR_YINDEX');
else
H.yindex = uicontrol(H.gui,'tag','XHR_YINDEX','Style','edit',...
'Units','Normalized',Font,...
'Position',[.92 .87 .08 .05],...
'BackGroundColor',[ 0 0 .9],'ForeGroundColor',[ 1 1 1],...
'TooltipString','Enter Y index into plot data matrix. Same as trace number.',...
'String','1','Value',1,...
'Callback',strcat('XHR = get(gcf,''userdata'');',...
'yi = str2num(get(XHR.yindex,''String''));',...
'XHR.data.yindex = yi;',...
'set(XHR.gui,''userdata'',XHR); clear XHR yi; ',...
'[Xpoint,Ypoint] = crosshair_subplots(''changey'');'));
end
if findobj(gcf,'tag','XHR_YPREV'),
H.yprev = findobj(gcf,'tag','XHR_YPREV');
else
H.yprev = uicontrol(H.gui,'tag','XHR_YPREV','Style','Push',...
'Units','Normalized',Font,...
'Position',[.92 .82 .04 .05],...
'String','<',...
'TooltipString','Goto Previous Y Index (channel).',...
'CallBack','[Xpoint,Ypoint] = crosshair_subplots(''prevy'');');
end
if findobj(gcf,'tag','XHR_YNEXT'),
H.ynext = findobj(gcf,'tag','XHR_YNEXT');
else
H.ynext = uicontrol(H.gui,'tag','XHR_YNEXT','Style','Push',...
'Units','Normalized',Font,...
'Position',[.96 .82 .04 .05],...
'String','>',...
'TooltipString','Goto Next Y Index (channel).',...
'CallBack','[Xpoint,Ypoint] = crosshair_subplots(''nexty'');');
end
if findobj(gcf,'tag','XHR_YLT'),
H.yLT = findobj(gcf,'tag','XHR_YLT');
else
H.yLT = uicontrol(H.gui,'tag','XHR_YLT','Style','Push',...
'Units','Normalized',Font,...
'Position',[.92 .77 .04 .05],...
'String','LT',...
'TooltipString','Goto next Y Less Than current Y.',...
'CallBack','[Xpoint,Ypoint] = crosshair_subplots(''yLT'');');
end
if findobj(gcf,'tag','XHR_YGT'),
H.yGT = findobj(gcf,'tag','XHR_YGT');
else
H.yGT = uicontrol(H.gui,'tag','XHR_YGT','Style','Push',...
'Units','Normalized',Font,...
'Position',[.96 .77 .04 .05],...
'String','GT',...
'TooltipString','Goto next Y Greater Than current Y.',...
'CallBack','[Xpoint,Ypoint] = crosshair_subplots(''yGT'');');
end
if findobj(gcf,'tag','XHR_XINDEX'),
H.xindex = findobj(gcf,'tag','XHR_XINDEX');
else
H.xindex = uicontrol(H.gui,'tag','XHR_XINDEX','Style','edit',...
'Units','Normalized',Font,...
'Position',[.92 .70 .08 .05],...
'BackGroundColor',[ 0 .9 0],'ForeGroundColor',[ 0 0 0],...
'TooltipString','Enter X index into plot data matrix. Only available for interpolation = ''none''.',...
'String','1','Value',1,...
'Callback',strcat('XHR = get(gcf,''userdata'');',...
'xi = str2num(get(XHR.xindex,''String''));',...
'XHR.data.xindex = xi;',...
'set(XHR.xinterp,''value'',0); ',...
'set(XHR.xinterp,''string'',''0''); ',...
'set(XHR.interp, ''value'',1); ',...
'set(XHR.gui,''userdata'',XHR); clear XHR; ',...
'[Xpoint,Ypoint] = crosshair_subplots(''changex'');'));
end
if findobj(gcf,'tag','XHR_XPREV'),
H.xprev = findobj(gcf,'tag','XHR_XPREV');
else
H.xprev = uicontrol(H.gui,'tag','XHR_XPREV','Style','Push',...
'Units','Normalized',Font,...
'Position',[.92 .65 .04 .05],...
'String','<',...
'TooltipString','Goto Previous X Index (no interpolation).',...
'CallBack','[Xpoint,Ypoint] = crosshair_subplots(''prevx'');');
end
if findobj(gcf,'tag','XHR_XNEXT'),
H.xnext = findobj(gcf,'tag','XHR_XNEXT');
else
H.xnext = uicontrol(H.gui,'tag','XHR_XNEXT','Style','Push',...
'Units','Normalized',Font,...
'Position',[.96 .65 .04 .05],...
'String','>',...
'TooltipString','Goto Next X Index (no interpolation).',...
'CallBack','[Xpoint,Ypoint] = crosshair_subplots(''nextx'');');
end
if findobj(gcf,'tag','XHR_XINTERP'),
H.xinterp = findobj(gcf,'tag','XHR_XINTERP');
else
H.xinterp = uicontrol(H.gui,'tag','XHR_XINTERP','Style','Edit',...
'Units','Normalized',...
'Position',[.92 .60 .08 .05],...
'String','0','Value',0,...
'TooltipString','Interpolation X increment (zero = nearest X).',...
'Callback',strcat('XHR = get(gcf,''userdata''); ',...
'xint = str2num(get(XHR.xinterp,''string'')); ',...
'set(XHR.xinterp,''value'',xint); ',...
'set(XHR.gui,''userdata'',XHR); figure(XHR.gui); clear XHR; '));
end
if findobj(gcf,'tag','XHR_INTERP'),
H.interp = findobj(gcf,'tag','XHR_INTERP');
else
interpstr = 'none|nearest|linear|spline|cubic';
H.interp = uicontrol(H.gui,'tag','XHR_INTERP','Style','popup',...
'Units','Normalized',Font,...
'Position',[.92 .55 .08 .05],...
'TooltipString','INTERP1 methods (none = raw values).', ...
'String',interpstr,...
'Callback',strcat('XHR = get(gcf,''userdata'');',...
'xint = get(XHR.xinterp,''Value''); ',...
'if xint == 0, ',...
' xint = 0.5; ',...
' set(XHR.xinterp,''value'',xint); ',...
' set(XHR.xinterp,''string'',num2str(xint)); ',...
' set(XHR.xindex, ''string'',''interp''); ',...
'end; ',...
'set(XHR.gui,''userdata'',XHR); figure(XHR.gui); clear XHR; '));
end
if findobj(gcf,'tag','XHR_STORE'),
H.store = findobj(gcf,'tag','XHR_STORE');
else
H.store = uicontrol(H.gui,'tag','XHR_STORE','Style','Push',...
'Units','Normalized',...
'Position',[.92 .40 .08 .05],...
'String','Store',...
'TooltipString','Store current Channel & XY values into base XYhist array.', ...
'CallBack','crosshair_subplots(''store'');');
end
if findobj(gcf,'tag','XHR_DONE'),
H.done = findobj(gcf,'tag','XHR_DONE');
else
H.done = uicontrol(H.gui,'tag','XHR_DONE','Style','Push',...
'Units','Normalized',...
'Position',[.80 .00 .10 .05],...
'BackgroundColor',[.8 .5 0],...
'ForegroundColor',[1 1 1],...
'FontWeight','bold',...
'String','Done',...
'TooltipString','Close crosshair', ...
'CallBack','crosshair_subplots(''done'');');
end
if findobj(gcf,'tag','XHR_EXIT'),
H.exit = findobj(gcf,'tag','XHR_EXIT');
else
H.exit = uicontrol(H.gui,'tag','XHR_EXIT','Style','Push',...
'Units','Normalized',...
'Position',[.90 .00 .10 .05],...
'BackgroundColor',[.8 0 0],...
'ForegroundColor',[1 1 1],...
'FontWeight','bold',...
'String','Exit',...
'TooltipString','Close crosshair and Figure', ...
'CallBack','crosshair_subplots(''exit'');');
end
if findobj(gcf,'tag','XHR_TRACENEAREST'),
H.traceNearest = findobj(gcf,'tag','XHR_TRACENEAREST');
else
H.traceNearest = uicontrol(H.gui,'tag','XHR_TRACENEAREST','Style','checkbox',...
'Units','Normalized',...
'Position',[.60 .00 .19 .05],...
'BackgroundColor',bgcolor,'ForegroundColor',fgcolor,...
'String','Nearest Trace','Value',1,...
'TooltipString','Trace nearest to mouse click; switch off to keep trace constant.');
end
% % Switch off||on Trace Selection GUI
% Vis = 'Off';
% if H.data.lines > 1,
% Vis = 'On';
% elseif H.data.lines == 0
% error('No lines found in the current axes (gca)\n');
% end
%
% % 'traces' string variable must be in ascending order
% traces = '';
% i = 1;
% while i <= lines;
% if i < lines
% tracelabel = sprintf('Column %4d|',i);
% else
% tracelabel = sprintf('Column %4d',i);
% end
% traces = strcat(traces,tracelabel);
% i = i + 1;
% end
%
% % If more than one line, provide GUI for line selection
%
%
% if findobj(gcf,'tag','XHR_TRACESWITCH'),
% H.trace = findobj(gcf,'tag','XHR_TRACESWITCH');
% set(H.trace,'String',traces);
% else
% H.trace = uicontrol(H.gui,'tag','XHR_TRACESWITCH','Style','Popup',...
% 'Units','Normalized',...
% 'Position',[.15 .00 .20 .05],...
% 'BackGroundColor','w','String',traces,...
% 'TooltipString','Select trace to follow with crosshairs.',...
% 'Visible',Vis,...
% 'CallBack','crosshair_subplots(''up'');');
% end
set(H.axis,'userdata',H.data);
set(H.gui,'userdata',H);
return
