image_levels

PURPOSE ^

IMAGE_LEVELS(img, levels, clim ) show slices at levels of an image

SYNOPSIS ^

function image_levels(img, levels, clim )

DESCRIPTION ^

 IMAGE_LEVELS(img, levels, clim  ) show slices at levels of an image
 img    = EIDORS image struct
 levels = array of vertical levels
  OR
 levels = [ [z_lev1 ,xpos,ypos], ...

 clim   = colourmap limit (or default if not specified)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function image_levels(img, levels, clim )
0002 % IMAGE_LEVELS(img, levels, clim  ) show slices at levels of an image
0003 % img    = EIDORS image struct
0004 % levels = array of vertical levels
0005 %  OR
0006 % levels = [ [z_lev1 ,xpos,ypos], ...
0007 %
0008 % clim   = colourmap limit (or default if not specified)
0009 
0010 % (C) 2005 Andy Adler. License: GPL version 2 or version 3
0011 % $Id: image_levels.m 3174 2012-06-27 15:14:07Z aadler $
0012 
0013 
0014 set(gcf,'NumberTitle','off');
0015 set(gcf,'Name', img.name);
0016 fwd_mdl= img.fwd_model;
0017 vtx=  fwd_mdl.nodes;
0018 simp= fwd_mdl.elems;
0019 img_data = img.elem_data;
0020 
0021 
0022 if size(levels,1)==3
0023    ll = size(levels,2);
0024    img_cols = max(levels(2,:));
0025    img_rows = max(levels(3,:));
0026 else
0027    ll = length( levels );
0028    img_cols = ceil( sqrt( ll ));
0029    img_rows = ceil( ll/ img_cols );
0030 end
0031    subplot(img_rows,img_cols,1);
0032 
0033 % Get geometry Fc
0034 fc = eidors_obj('get-cache', fwd_mdl, 'slicer_plot_fc');
0035 if ~isempty( fc )
0036     eidors_msg('image_levels: using cached value', 3);
0037 else
0038    [fc] = slicer_plot_n(levels(1),img_data,vtx,simp);
0039    eidors_obj('set-cache', fwd_mdl, 'slicer_plot_fc', fc);
0040    eidors_msg('image_levels: setting cached value', 3);
0041 end
0042 
0043 % Set mapped colours
0044 global eidors_colours;
0045 mpc= eidors_colours.mapped_colour;
0046 eidors_colours.mapped_colour = 128;
0047 
0048 for idx= 1:length(levels);
0049     subplot(img_rows,img_cols,idx);
0050     lev= levels(idx);
0051     slicer_plot_n(lev,img_data,vtx,simp, fc);
0052     view(2);
0053     axis('off');
0054 if ~exist('OCTAVE_VERSION'); % octave gives ft_render warnings
0055     axis equal
0056 end
0057     title(sprintf('z=%4.2f',lev));
0058 end
0059 
0060 % Reset mapped colours
0061 eidors_colours.mapped_colour = mpc;

Generated on Wed 29-May-2013 17:11:47 by m2html © 2005