Reject the selected channels. This is interface function.
[MEGinfo,bexp,pick,Qpick,CoilWeight] = ...
vb_channel_rejection(MEGinfo,bexp,pick,Qpick,CoilWeight,reject_channel)
[IN]
MEGinfo : struct with MEG acquisition information
.MEGch_id Active MEG channel index (absolute channel inndex)
.sensor_weight % weight to calculate lead field
.Nchannel % number of active gardiometers
bexp(ch,t,trial) : MEG data (Channel, Time, Trial)
pick(ch, 1:3) : position of detector coils.
(ch+Nchannel, 1:3) : position of conpensation coils.
Qpick(ch, 1:3) : direction of detector coils.
(ch+Nchannel, 1:3) : direction of conpensation coils.
CoilWeight : for SBI.
reject_channel : vector of rejected channel number .
absolute channel inndex
[OUT]
MEGinfo : struct with MEG acquisition information after rejection
.MEGch_id Active MEG channel index (absolute channel inndex)
.sensor_weight % weight to calculate lead field
.Nchannel % number of active gardiometers after rejection
bexp(ch,t,trial) : MEG data (Channel, Time, Trial) after rejection
pick(ch, 1:3) : position of detector coils.
(ch+Nchannel, 1:3) : position of conpensation coils.
Qpick(ch, 1:3) : direction of detector coils.
(ch+Nchannel, 1:3) : direction of conpensation coils.
CoilWeight : for SBI.
MEG channel index
MEGinfo.MEGch_id : Active MEG channel index (absolute channel inndex)
bexp(ch,:) <-> MEGinfo.MEGch_id(n)
pick(ch, 1:3) <-> MEGinfo.MEGch_id(n)
(ch+Nchannel, 1:3) <-> MEGinfo.MEGch_id(n)
Qpick(ch, 1:3) <-> MEGinfo.MEGch_id(n)
(ch+Nchannel, 1:3) <-> MEGinfo.MEGch_id(n)
CoilWeight
2006-05-08 Masanori Osako
2006-05-14 M. Sato
2006-06-19 D. Kawawaki modified based on "yokogawa_channel_rejection".
Copyright (C) 2011, ATR All Rights Reserved.
License : New BSD License(see VBMEG_LICENSE.txt)0001 function [MEGinfo,bexp,pick,Qpick,CoilWeight] = ... 0002 vb_channel_rejection(MEGinfo,bexp,pick,Qpick,reject_channel,CoilWeight) 0003 % 0004 % Reject the selected channels. This is interface function. 0005 % 0006 % [MEGinfo,bexp,pick,Qpick,CoilWeight] = ... 0007 % vb_channel_rejection(MEGinfo,bexp,pick,Qpick,CoilWeight,reject_channel) 0008 % [IN] 0009 % MEGinfo : struct with MEG acquisition information 0010 % .MEGch_id Active MEG channel index (absolute channel inndex) 0011 % .sensor_weight % weight to calculate lead field 0012 % .Nchannel % number of active gardiometers 0013 % bexp(ch,t,trial) : MEG data (Channel, Time, Trial) 0014 % pick(ch, 1:3) : position of detector coils. 0015 % (ch+Nchannel, 1:3) : position of conpensation coils. 0016 % Qpick(ch, 1:3) : direction of detector coils. 0017 % (ch+Nchannel, 1:3) : direction of conpensation coils. 0018 % CoilWeight : for SBI. 0019 % reject_channel : vector of rejected channel number . 0020 % absolute channel inndex 0021 % [OUT] 0022 % MEGinfo : struct with MEG acquisition information after rejection 0023 % .MEGch_id Active MEG channel index (absolute channel inndex) 0024 % .sensor_weight % weight to calculate lead field 0025 % .Nchannel % number of active gardiometers after rejection 0026 % bexp(ch,t,trial) : MEG data (Channel, Time, Trial) after rejection 0027 % pick(ch, 1:3) : position of detector coils. 0028 % (ch+Nchannel, 1:3) : position of conpensation coils. 0029 % Qpick(ch, 1:3) : direction of detector coils. 0030 % (ch+Nchannel, 1:3) : direction of conpensation coils. 0031 % CoilWeight : for SBI. 0032 % 0033 % MEG channel index 0034 % MEGinfo.MEGch_id : Active MEG channel index (absolute channel inndex) 0035 % bexp(ch,:) <-> MEGinfo.MEGch_id(n) 0036 % pick(ch, 1:3) <-> MEGinfo.MEGch_id(n) 0037 % (ch+Nchannel, 1:3) <-> MEGinfo.MEGch_id(n) 0038 % Qpick(ch, 1:3) <-> MEGinfo.MEGch_id(n) 0039 % (ch+Nchannel, 1:3) <-> MEGinfo.MEGch_id(n) 0040 % CoilWeight 0041 % 0042 % 2006-05-08 Masanori Osako 0043 % 2006-05-14 M. Sato 0044 % 2006-06-19 D. Kawawaki modified based on "yokogawa_channel_rejection". 0045 % 0046 % Copyright (C) 2011, ATR All Rights Reserved. 0047 % License : New BSD License(see VBMEG_LICENSE.txt) 0048 0049 if ~isfield(MEGinfo, 'device') 0050 device = 'SBI'; 0051 else 0052 device = MEGinfo.device; 0053 end 0054 0055 switch (lower(MEGinfo.device)) 0056 case 'sbi' 0057 [MEGinfo,bexp,pick,Qpick,CoilWeight] = ... 0058 sbi_channel_rejection(MEGinfo,bexp,pick,Qpick,CoilWeight,reject_channel); 0059 case 'yokogawa' 0060 [MEGinfo,bexp,pick,Qpick] = ... 0061 vb_yokogawa_channel_rejection(MEGinfo,bexp,pick,Qpick,reject_channel); 0062 CoilWeight = MEGinfo.sensor_weight; 0063 end 0064 0065 return