function h1 = imageplot(M,str, a,b,c)% imageplot - diplay an image and a title%% Example of usages:%   imageplot(M);%   imageplot(M,title);%   imageplot(M,title,1,2,1);   % to make subplot(1,2,1);%%   imageplot(M,options);%%   If you want to display several images:%       imageplot({M1 M2}, {'title1', 'title2'});%%   Copyright (c) 2007 Gabriel Peyreif nargin<2    str = [];endoptions.null = 0;if isstruct(str)    options = str;    str = '';endnbdims = nb_dims(M);if iscell(M)    q = length(M);    if nargin<5        c = 1;    end    if nargin<4        a = ceil(q/4);    end    if nargin<3        b = ceil(q/a);    end    if (c-1+q)>(a*b)        warning('a and c parameters not large enough');        a = ceil((c-1+q)/4);        b = ceil((c-1+q)/a);    end    for i=1:q        if iscell(str)            str1 = str{i};        else            str1 = str;        end        h{i} = imageplot(M{i},str1, a,b,c-1+i);    end    global axlist;    if not(isempty(axlist))        linkaxes(axlist, 'xy');    end    if nargout>0        if exist('h')            h1 = h;        else            h1 = [];        end    end    return;endif nargin==5    global axlist;    global imageplot_size;    if c==1 || isempty(imageplot_size) || imageplot_size~=size(M,1)        clear axlist;         global axlist;         axlist = [];        imageplot_size = size(M,1);    end    axlist(end+1) = subplot(a,b,c);endif nbdims==1    h = plot(M); axis tight;elseif size(M,3)<=3    % gray-scale or color image    if size(M,3)==2        M = cat(3,M, zeros(size(M,1),size(M,2)));    end    if not(isreal(M))        if size(M,3)==1            % turn into color matrix            M = cat(3, real(M), imag(M), zeros(size(M,1),size(M,2)));        else            warning('Complex data');            M = real(M);        end    end    if size(M,3)==1        colormap gray(256);    else        colormap jet(256);    end    h = imagesc(rescale(M)); axis image; axis off;else    if not(isfield(options, 'center') )        options.center = .5; % here a value in [0,1]    end    if not(isfield(options, 'sigma'))        options.sigma = .08; % control the width of the non-transparent region    end    a = compute_alpha_map('gaussian', options); % you can plot(a) to see the alphamap    % volumetric image    h = vol3d('cdata',rescale(M),'texture','2D');    view(3);    axis tight; % daspect([1 1 .4])    colormap bone(256);%    alphamap('rampup');    % alphamap(.06 .* alphamap);    vol3d(h);endif not(isempty(str))    title(str);endif nargout>0     if exist('h')        h1 = h;    else         h1 = [];    endendif nargin==5 && c==a*b    linkaxes(axlist, 'xy');endfunction d = nb_dims(x)% nb_dims - debugged version of ndims.%%   d = nb_dims(x);%%   Copyright (c) 2004 Gabriel Peyr