0001 function vb_check_sphere(brainfile,megfile,Vradius)
0002
0003
0004
0005
0006
0007
0008 if ~exist('Vradius','var') || isempty(Vradius)
0009 Vradius = 0.08;
0010 else
0011 Vradius = Vradius;
0012 end;
0013
0014
0015
0016
0017
0018
0019 [pick,Qpick,Wsensor,Vcenter] = vb_load_sensor(megfile);
0020
0021
0022 [V,F] = vb_load_cortex(brainfile);
0023 F = F.F3;
0024
0025
0026 NY = 2;
0027 NX = 2;
0028
0029
0030 fclr = [0.9 0.7 0.6];
0031 eclr = 'none';
0032
0033 light_mode = 1;
0034 max_mode = 0;
0035
0036
0037 vangle = [ 0 90; 90 0 ; 180 0 ; 0 0 ];
0038 NV = size(vangle,1);
0039
0040 [sx,sy,sz]=sphere;
0041 sdim = size(sx);
0042 sx = sx*Vradius + repmat(Vcenter(1),sdim);
0043 sy = sy*Vradius + repmat(Vcenter(2),sdim);
0044 sz = sz*Vradius + repmat(Vcenter(3),sdim);
0045
0046 figure;
0047
0048 for n=1:NV,
0049 subplot(NY,NX,n);
0050
0051
0052 vb_plot_surf(V, F, fclr,eclr,light_mode,max_mode)
0053 hold on;
0054 surf(sx,sy,sz);
0055 view( vangle(n,:) );
0056 xlabel('X')
0057 ylabel('Y')
0058 zlabel('Z')
0059 alpha(0.2);
0060 end
0061
0062 subplot(NY,NX,1);
0063 title('Sphere model overlapped on brain model')
0064