vb_check_sensor_direction

 plot pick and Qpick on the same figure and plot Qpick and pick vector
 
-- Syntax
 vb_check_sensor_direction(megfile)

-- Input 
 megfile : meg.mat file name  

-- Output
 figure1 : pick and Qpick
 figure2 : pick-Qpick vector

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

0001 function    vb_check_sensor_direction(megfile)
0002 % plot pick and Qpick on the same figure and plot Qpick and pick vector
0003 %
0004 %-- Syntax
0005 % vb_check_sensor_direction(megfile)
0006 %
0007 %-- Input
0008 % megfile : meg.mat file name
0009 %
0010 %-- Output
0011 % figure1 : pick and Qpick
0012 % figure2 : pick-Qpick vector
0013 %
0014 % Copyright (C) 2011, ATR All Rights Reserved.
0015 % License : New BSD License(see VBMEG_LICENSE.txt)
0016 
0017 [pick,Qpick,Wsensor,Vcenter] = vb_load_sensor(megfile);
0018 
0019 figure
0020 pick_num = size(pick,1)/2;
0021 plot3(pick(1:pick_num,1),pick(1:pick_num,2),pick(1:pick_num,3),'ob');
0022 hold on; 
0023 plot3(pick(pick_num+1:end,1),pick(pick_num+1:end,2),pick(pick_num+1:end,3),'or');
0024 title('pick(blue point) and Qpick(red point)')
0025 axis equal
0026 
0027 figure
0028 Qpick_pick=0.005*Qpick(1:pick_num,:)+pick(1:pick_num,:);
0029 plot3(pick(1:pick_num,1),pick(1:pick_num,2),pick(1:pick_num,3),'ob');
0030 hold on
0031 for n=1:pick_num,
0032     plot3([pick(n,1);Qpick_pick(n,1)],[pick(n,2);Qpick_pick(n,2)],[pick(n,3);Qpick_pick(n,3)],'-');
0033 end
0034 title('pick-Qpick vector')
0035 axis equal