sensor_pos_neuromag

 sensor position of neuromag
 [usage]

 [input]
   MEGid : <required> MEG channel type. A type of this varible must be a
            charactor string or integer as followed.
            char)
              'VV_PLANAR_T1' or 'VV_PLANAR_T2', 'VV_PLANER_T3'
              'VV_MAG_T1' or 'VV_MAG_T2', 'VV_MAG_T3'
              'Neuromag_122'
            integer)
               3012 or 3013 or 3014
               3022 or 3023 or 3024
               2
       U : <required> transformation matrix [4 x 4]
 [output]
    pick : coil position [n_coil x 3]
   Qpick : coil direction [n_coil x 3]
       w : 
 [note]
 [history]
　　Masa-aki Sato 2009-4-13
    M.Fukushima   2012/11/15
      * modified to import Henson dataset
    rhayashi      2013-9-9   Added support VV_MAG_T3, VV_PLANER_T3

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

0001 function    [pick, Qpick, w] = sensor_pos_neuromag(MEGid,U)
0002 % sensor position of neuromag
0003 % [usage]
0004 %
0005 % [input]
0006 %   MEGid : <required> MEG channel type. A type of this varible must be a
0007 %            charactor string or integer as followed.
0008 %            char)
0009 %              'VV_PLANAR_T1' or 'VV_PLANAR_T2', 'VV_PLANER_T3'
0010 %              'VV_MAG_T1' or 'VV_MAG_T2', 'VV_MAG_T3'
0011 %              'Neuromag_122'
0012 %            integer)
0013 %               3012 or 3013 or 3014
0014 %               3022 or 3023 or 3024
0015 %               2
0016 %       U : <required> transformation matrix [4 x 4]
0017 % [output]
0018 %    pick : coil position [n_coil x 3]
0019 %   Qpick : coil direction [n_coil x 3]
0020 %       w :
0021 % [note]
0022 % [history]
0023 %　　Masa-aki Sato 2009-4-13
0024 %    M.Fukushima   2012/11/15
0025 %      * modified to import Henson dataset
0026 %    rhayashi      2013-9-9   Added support VV_MAG_T3, VV_PLANER_T3
0027 %
0028 % Copyright (C) 2011, ATR All Rights Reserved.
0029 % License : New BSD License(see VBMEG_LICENSE.txt)
0030 
0031 %fprintf('MEGid=')
0032 if ischar(MEGid),
0033     switch    MEGid
0034     case    {'VV_PLANAR_T1','VV_PLANAR_T2','VV_PLANER_T3'}
0035         % 3012 : Vectorview type 1 planar gradiometer
0036         % 3013 : Vectorview type 2 planar gradiometer
0037         % ( +/-8.4 ,0.0 , 0.3 ) mm
0038         % +/- 1/16.8mm
0039         
0040         % sensor position
0041         x1 = [  8.4 , 0.0 , 0.3]*0.001;
0042         x2 = [ -8.4 , 0.0 , 0.3]*0.001;
0043         z  = [  0.0 , 0.0 , 1.0];
0044         
0045         w     = [ 1/16.8 , -1/16.8]*1000;
0046         pick  = [x1, 1; x2, 1] * U(:,1:3);
0047         Qpick = [z; z] * U(1:3,1:3);
0048     %    fprintf('%d\n',MEGid)
0049     case    {'VV_MAG_T1','VV_MAG_T2'}
0050         % 3022 : Vectorview type 1 magnetometer
0051         % 3023 : Vectorview type 2 magnetometer
0052         % (±6.45,±6.45,0.3) mm
0053         % 1/4
0054         
0055         % sensor position
0056         x1 = [  6.45 , 6.45 , 0.3]*0.001;
0057         x2 = [  6.45 ,-6.45 , 0.3]*0.001;
0058         x3 = [ -6.45 , 6.45 , 0.3]*0.001;
0059         x4 = [ -6.45 ,-6.45 , 0.3]*0.001;
0060         z  = [  0.0 , 0.0 , 1.0];
0061         
0062         w     = [ 1/4 ,1/4 ,1/4 ,1/4 ];
0063         pick  = [x1, 1; x2, 1; x3, 1; x4, 1] * U(:,1:3);
0064         Qpick = [z; z; z; z] * U(1:3,1:3);
0065     %    fprintf('%d\n',MEGid)
0066     case {'VV_MAG_T3'}
0067         % 3024 : Vectorview type 3 magnetometer
0068         % (±5.25,±5.25,0.3) mm
0069         % 1/4
0070         
0071         % sensor position
0072         x1 = [  5.25 , 5.25 , 0.3]*0.001;
0073         x2 = [  5.25 ,-5.25 , 0.3]*0.001;
0074         x3 = [ -5.25 , 5.25 , 0.3]*0.001;
0075         x4 = [ -5.25 ,-5.25 , 0.3]*0.001;
0076         z  = [  0.0 , 0.0 , 1.0];
0077         
0078         w     = [ 1/4 ,1/4 ,1/4 ,1/4 ];
0079         pick  = [x1, 1; x2, 1; x3, 1; x4, 1] * U(:,1:3);
0080         Qpick = [z; z; z; z] * U(1:3,1:3);
0081     %    fprintf('%d\n',MEGid)
0082     case    {'Neuromag_122'}
0083         % 2 : Neuromag-122 : planar gradiometer
0084         % ±(8.1, 0, 0)mm
0085         % ±1/16.2mm
0086         %
0087         
0088         % sensor position
0089         x1 = [  8.1 , 0.0 , 0.0]*0.001;
0090         x2 = [ -8.1 , 0.0 , 0.0]*0.001;
0091         z  = [  0.0 , 0.0 , 1.0];
0092         
0093         w     = [ 1/16.2 , -1/16.2]*1000;
0094         pick  = [x1, 1; x2, 1] * U(:,1:3);
0095         Qpick = [z; z] * U(1:3,1:3);
0096     %    fprintf('%d\n',MEGid)
0097     otherwise
0098         w      = [];
0099         pick  = [];
0100         Qpick = [];
0101     %    fprintf('\n')
0102     end
0103 else
0104     switch    MEGid
0105     case    {3012,3013,3014}
0106         % 3012 : Vectorview type 1 planar gradiometer
0107         % 3013 : Vectorview type 2 planar gradiometer
0108         % 3014 : Vectorview type 3 planar gradiometer
0109         % ( +/-8.4 ,0.0 , 0.3 ) mm
0110         % +/- 1/16.8mm
0111         
0112         % sensor position
0113         x1 = [  8.4 , 0.0 , 0.3]*0.001;
0114         x2 = [ -8.4 , 0.0 , 0.3]*0.001;
0115         z  = [  0.0 , 0.0 , 1.0];
0116         
0117         w     = [ 1/16.8 , -1/16.8]*1000;
0118         pick  = [x1, 1; x2, 1] * U(:,1:3);
0119         Qpick = [z; z] * U(1:3,1:3);
0120     %    fprintf('%d\n',MEGid)
0121     case    {3022,3023}
0122         % 3022 : Vectorview type 1 magnetometer
0123         % 3023 : Vectorview type 2 magnetometer
0124         % (±6.45,±6.45,0.3) mm
0125         % 1/4
0126         
0127         % sensor position
0128         x1 = [  6.45 , 6.45 , 0.3]*0.001;
0129         x2 = [  6.45 ,-6.45 , 0.3]*0.001;
0130         x3 = [ -6.45 , 6.45 , 0.3]*0.001;
0131         x4 = [ -6.45 ,-6.45 , 0.3]*0.001;
0132         z  = [  0.0 , 0.0 , 1.0];
0133         
0134         w     = [ 1/4 ,1/4 ,1/4 ,1/4 ]; % fT
0135         pick  = [x1, 1; x2, 1; x3, 1; x4, 1] * U(:,1:3);
0136         Qpick = [z; z; z; z] * U(1:3,1:3);
0137     %    fprintf('%d\n',MEGid)
0138     case {3024}
0139         % 3024 : Vectorview type 3 magnetometer
0140         % (±5.25,±5.25,0.3) mm
0141         % 1/4
0142         
0143         % sensor position
0144         x1 = [  5.25 , 5.25 , 0.3]*0.001;
0145         x2 = [  5.25 ,-5.25 , 0.3]*0.001;
0146         x3 = [ -5.25 , 5.25 , 0.3]*0.001;
0147         x4 = [ -5.25 ,-5.25 , 0.3]*0.001;
0148         z  = [  0.0 , 0.0 , 1.0];
0149         
0150         w     = [ 1/4 ,1/4 ,1/4 ,1/4 ]; % fT
0151         pick  = [x1, 1; x2, 1; x3, 1; x4, 1] * U(:,1:3);
0152         Qpick = [z; z; z; z] * U(1:3,1:3);
0153     %    fprintf('%d\n',MEGid)
0154     case    {2}
0155         % 2 : Neuromag-122 : planar gradiometer
0156         % ±(8.1, 0, 0)mm
0157         % ±1/16.2mm
0158         %
0159         
0160         % sensor position
0161         x1 = [  8.1 , 0.0 , 0.0]*0.001;
0162         x2 = [ -8.1 , 0.0 , 0.0]*0.001;
0163         z  = [  0.0 , 0.0 , 1.0];
0164         
0165         w     = [ 1/16.2 , -1/16.2]*1000;
0166         pick  = [x1, 1; x2, 1] * U(:,1:3);
0167         Qpick = [z; z] * U(1:3,1:3);
0168     otherwise
0169         w      = [];
0170         pick  = [];
0171         Qpick = [];
0172     %    fprintf('\n')
0173     end
0174 end