function lvmClassVisualisePath(call)

% LVMCLASSVISUALISEPATH Latent variable model path drawing in latent space.
%
%	Description:
%	lvmClassVisualisePath(call)
%% 	lvmClassVisualisePath.m SVN version 24
% 	last update 2007-11-08T09:36:53.000000Z

global visualiseInfo


switch call
 case 'click'
  visualiseInfo.clicked = ~visualiseInfo.clicked;
  if visualiseInfo.clicked
    [x, y]  = localCheckPointPosition(visualiseInfo);  
    if ~isempty(x) 
      visualiseInfo.latentPos = [x, y];
      visualiseInfo.posTrail = [x, y];
      visualiseInfo.lastTime = clock;
      visualiseInfo.timeTrail = 0;
    end
  end
 case 'move'
  if visualiseInfo.clicked 
    [x, y]  = localCheckPointPosition(visualiseInfo);  
    if ~isempty(x) 
      visualiseInfo.latentPos=[x, y];
      % Draw connecting lines
      iplot(visualiseInfo.plotAxes, x, y);
      visualiseInfo.posTrail = [visualiseInfo.posTrail; [x, y]];
      timeNow = clock;
      visualiseInfo.timeTrail = [visualiseInfo.timeTrail; visualiseInfo.timeTrail(end)+etime(timeNow, ...
                                                        visualiseInfo.lastTime)];
      visualiseInfo.lastTime = timeNow;
    end
  end
end




function point = localGetNormCursorPoint(figHandle)

point = get(figHandle, 'currentPoint');
figPos = get(figHandle, 'Position');
% Normalise the point of the curstor
point(1) = point(1)/figPos(3);
point(2) = point(2)/figPos(4);

function [x, y] = localGetNormAxesPoint(point, axesHandle)

position = get(axesHandle, 'Position');
x = (point(1) - position(1))/position(3);
y = (point(2) - position(2))/position(4);
lim = get(axesHandle, 'XLim');
x = x*(lim(2) - lim(1));
x = x + lim(1);
lim = get(axesHandle, 'YLim');
y = y*(lim(2) - lim(1));
y = y + lim(1);


function [x, y] = localCheckPointPosition(visualiseInfo)

% Get the point of the cursor
point = localGetNormCursorPoint(gcf);

% get the position of the axes
position = get(visualiseInfo.plotAxes, 'Position');


% Check if the pointer is in the axes
if point(1) > position(1) ...
      & point(1) < position(1) + position(3) ...
      & point(2) > position(2) ...
      & point(2) < position(2) + position(4);
  
  % Rescale the point according to the axes
  [x y] = localGetNormAxesPoint(point, visualiseInfo.plotAxes);

  % Find the nearest point
else
  % Return nothing
  x = [];
  y = [];
end
