plot_subplot_sequence.m

Go to the documentation of this file.

function p = plot_subplot_sequence(varargin)
%function p = plot_subplot_sequence(data,grid,params[,callbackfigure, cbhandle])
% plotting a sequence of data slices on polygonal 2d grid (constructed from
% params if empty) and providing a slider for changing between the data slices.
%
% A new figure is opened and the handle returned in p.  If further parameters
% are set, the call is assumed to stem from a callback-function
%
% Every column of data is interpreted as one discrete function
% dof vector forwarded to the params.plot() function.
%
% parameters:
%  varargin: usually called with 3 arguments:
%            @code plot_sequence(data, grid, params) @endcode
%   - data - cell of data vectors to be plotted
%   - grid - the underlying grid
%   - params - plotting parameters
%   .
%            Alternatively, there can be 5 arguments:
%            @code plot_sequence(data, grid, params, callbackfigure, cbhandle) @endcode
%            This function is set as a callback function for the time slider
%            uicontrol. The additional arguments are:
%   - callbackfigure - handle to the figure
%   - cbhandle       - handle to the object calling the callback function. This
%                      is usually the time slider.
%
% return values:
%  p : figure handle of plot
%
% required fields in params:
%  - 'plot' -- pointer to the plot-function performing the plotting of a
%              single slice, e.g.  plot_element_data(), plot_vertex_data(),
%              fv_plot(), ldg_plot().
%
% optional field of params:
%  - 'title' -- string indicating the title of the newly opened figure
%  - 'clim'  -- 2-vector giving the range of the colormap. If this is set,
%               identical range is used for all slices. Default is the min and
%               max of all slices.
%  - 'clim_tight' -- if this flag is set, the colorbar is set tightly to the
%                    range of every single data slice. Clearly only one of
%                    'clim' or 'clim_tight' should be set.
%
% see also the chosen 'params.plot_function' for its further params-options

% Bernard Haasdonk 9.5.2007

% open figure
if nargin < 4 % no callback, open new figure

  %  keyboard;
  f = figure;
  di = 0.05;
  textwidth= 0.15;
  textheight = 0.03;
  panel = uipanel(f, 'Title','Control Panel','FontSize', 12, ...
                  'Units', 'normalized', ...
                  'Position', [0.5*di, di, 1-di, 2.0*textheight]);
  txt = uicontrol('Parent', panel,'Style','text','Units','normalized',...
                  'Position',[0.5*di,di,textwidth,0.9],...
                  'String','Slice 1','Tag','text1');
  slider = uicontrol('Parent', panel,'Style','slider','Units','normalized',...
                     'Position',[di+textwidth, di ,1-3*di-textwidth, 0.9],...
                     'Tag','slider1',...
                     'Callback','plot_subplot_sequence([],[],[],[],gcbf,gcbo)');
  ax = uipanel(f, 'Title', 'Plots', 'Units', 'normalized', 'FontSize', 12, ...
               'Position', [0.5*di, 2*di+textheight, 1-di, 1-3*di+textheight], ...
               'Tag', 'axesarea');

  %  cb = colorbar;
  ud = [];
  ud.data = varargin{1};
  ud.grid = varargin{2};
  ud.params = varargin{3};

  if isempty(ud.grid)
    ud.grid = construct_grid(ud.params);
  end;

  if isfield(ud.params,'title')
    set(f,'Name',ud.params.title);
  end;

  if ~isfield(ud.params,'clim_tight')
    if isfield(ud.params, 'clim')
      ud.params.clim_tight = 0;
    else
      ud.params.clim_tight = 1;
    end
  end

  set(f,'Userdata',ud);

  sl = findobj(f,'Tag','slider1');
  %set(sl,'min',0);
  %set(sl,'max',ud.params.nt);
  %set(sl,'sliderstep',[1/ud.params.nt,1]);
  set(sl,'min',1);
  if size(ud.data{1},2)>1
    set(sl,'max',size(ud.data{1},2));
    set(sl,'sliderstep',[1/(size(ud.data{1},2)-1),1]);
  else
    set(sl,'max',size(ud.data{1},2)+eps);
    set(sl,'visible','off');
  end;
  set(sl,'value',1);
  th = findobj(f,'Tag','text1');
  set(th,'String',('data slice 1') );
  replot(f,sl,1);

  p = f;
end

if nargin >=4 % callback-function
  cbf = varargin{4};
  cbo = varargin{5};
  %  disp('performing callback call of plot_fv_data')
  th = findobj(gcbf,'Tag','text1');
  v = round(get(gcbo,'value'));
  set(gcbo,'value',v)
  set(th,'String',['data slice ',num2str(v)]);
  replot(gcbf,gcbo,v);
  p = gcbf;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function replot(fh,oh,v)
ud = get(fh,'Userdata');

if ~iscell(ud.params.clim_tight)
  clt = cell(1, length(ud.data));
  for i = 1:length(clt)
    clt{i} = ud.params.clim_tight;
  end
  ud.params.clim_tight = clt;
end

ax = findobj(fh,'Tag','axesarea');
for i = 1:length(ud.data);
  subplot(ud.params.sp_x, ud.params.sp_y, ud.params.index{i}, 'Parent', ax);
  d = ud.data{i};
  if isstruct(d)
    d.t = v;
  else
    d = d(:,v);
  end

  % if clim_tight, adapt color-range:
  if ud.params.clim_tight{i} && isnumeric(d)
    ud.params.clim{i} = [min(min(d)), max(max(d))];
    if ud.params.clim{i}(1) == ud.params.clim{i}(2)
      ud.params.clim{i} = ud.params.clim{i} + [-eps, +eps];
    end
  elseif ~isfield(ud.params, 'clim') || ~iscell(ud.params.clim) || length(ud.params.clim) < i
    ud.params.clim{i} = [];
  end

  plot_params = structfun(@(x) extract_plot_param(x, i), ud.params, 'UniformOutput', false);

  g = ud.params.plot{i}(ud.grid, d, plot_params);

%  if ~isempty(ud.params.clim)
%    set(get(g, 'Parent'),'Clim',ud.params.clim{i});
%  end
end

function param = extract_plot_param(x,index)

if iscell(x)
  param = x{index};
else
  param = x;
end