0001 function image_levels(img, levels, clim )
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
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
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
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');
0055 axis equal
0056 end
0057 title(sprintf('z=%4.2f',lev));
0058 end
0059
0060
0061 eidors_colours.mapped_colour = mpc;