vb_plot_mask

 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function    vb_plot_mask(v,zlist,mode,X0,r0)
0002 %
0003 % Copyright (C) 2011, ATR All Rights Reserved.
0004 % License : New BSD License(see VBMEG_LICENSE.txt)
0005 
0006 if nargin<3, mode=0; end;
0007 
0008 NZ=length(zlist);
0009 
0010 NX=ceil(sqrt(NZ*4/3));
0011 NY=ceil(NZ/NX);
0012 
0013 if nargin>=4,
0014     x0=X0(1);
0015     y0=X0(2);
0016     z0=X0(3);
0017 end;
0018 
0019 if nargin<5
0020     r0=10;
0021 end;
0022 
0023 N=100;
0024 t=1:N;
0025 t=(2*pi/N)*t;
0026 
0027 
0028 for n=1:NZ
0029     zt=zlist(n);
0030     subplot(NY,NX,n);
0031     
0032     xy=squeeze(v(:,:,zt));
0033     
0034     switch    mode,
0035     case    0,
0036         h=image(xy,'CDataMapping','scaled');
0037         hold on
0038         if nargin >= 4,
0039             r1= r0^2 - (zt-z0)^2;
0040             if r1 > 0,
0041                 r=sqrt(r1);
0042                 xx=r*sin(t) + x0;
0043                 yy=r*cos(t) + y0;
0044                 plot(xx,yy,'-r');
0045             end;
0046         end;
0047         colormap('gray'); 
0048     case    1,
0049         mesh(xy);
0050         colormap('jet'); 
0051         zlim([0 1]);
0052     end;
0053     
0054     tlabel=sprintf('Z = %d',zt);
0055     title(tlabel);
0056 
0057 end;