EEG Processor Users Manual
変更履歴
2007/08/10 ver 0.1A (Sako)
第一版
2007/08/23 ver 0.1B (Sako)
BDFファイルロード時のBDF_SPEC構造体にfigure_info構造体フィールドを追加した。
この文書では、VBMEGが提供するEEGデータ処理プログラムの使い方について解説します。
BDFファイルを持つEEGユーザが対象です。
以下の文書中では、このEEGデータ処理プログラムをEEGプロセッサと呼びます。
現状の実装で想定している環境は、次の通りです。
l BIOSEMI社のEEG測定器によるBDF形式ファイルを入力ファイルとする
l MATLAB6.5以上
EEGプロセッサは、次の処理をすることができます。
1) BDFファイルから、チャネル毎のデータファイルを作成する
2) 対象となるチャネル、サンプルを指定してEEGデータを抽出する
3) トリガ情報を取得する(自動または手動)
4) トリガ情報を基に、EEGデータを分割する
サンプルデータを用いて、典型的な使い方を実行してみましょう。
チュートリアルの実行に先立ち、以下の環境を確認します。
l vbmegパスを通します(vbmeg関数を使用しているため)
l
以下に置いてあるサンプルデータを使用します
/home/cbi-data/sakot/eeg/sample_data/
l
チュートリアルは次のテストプログラムを実行します
$VBMEG/script/eeg_tutorial.m
このプログラムが存在することを確認します
パラメータ設定例は、このテストプログラムの記述を参考にしてください。パラメータ設定については後述します。
チュートリアルの実行は、MATLABコンソールからコマンドを入力することで行います。
(例)「コード1で実行する」
% eeg_tutorial(1)
処理内容は次の通りです。
1) BDFファイルを読む
2) チャネルデータファイルを作成する
3) EEG-MATファイルを作成する
コード1で実行します。
% eeg_tutorial(1)
内部ではvb_make_msrmnt_data関数が呼ばれます。vb_make_msrmnt_data関数については後述します。
bdf_spec.figure_info を新規に追加しました。
これは従来の”digit_file” および “face_file” に代わるもので、それぞれのファイルから取得していたパラメータを直接指定することができます。これにより、ファイルがなくても動かすことができるようになりました。この構造体の仕様は次のとおりです。
l .SensorPosition : センサー位置座標 [Nsensor x 3] (digit_file)
l .CoordType : コードタイプ<文字列> (例) “SPM_Right_m” 等 (digit_file)
l ChannelLabel : チャネルラベル <文字列のセル配列> {NChannel x 1} (digit_file)
l .Vcenter : 中心座標 [1 x 3] (face_file)
l bdf_spec.bdf_file = BDFファイル(.bdf) : (sample) sample_mrmt.bdf
l bdf_spec.digit_file = DIGITIZERファイル(.digit.mat) : (sample) sample_mrmt.digit.mat
l bdf_spec.face_file = FACEファイル(.face.mat) : (sample) sample_mrmt.face.mat
l bdf_spec.device = [] (default)
l bdf_spec.primary_path = ‘../test_dir/test’
l bdf_spec.label_type = [] (default)
l bdf_spec.eeginfo_version = [] (default)
l ch_spec = [] (チャネルデータファイルの読み込みは行わない)
カレントディレクトリ以下の’test_dir’ディレクトリに、以下のデータが作成されます。
l test.primary.eeg.mat (初期EEG-MATファイル)
l test_${channel_label}.ch.dat
センサの位置情報は、DIGITファイルが保持しています。一方、各センサの測定データはBDFファイルが保持しています。センサのことをチャネルとも呼びます。それぞれのファイルで、それぞれのチャネルに名前が付けられています。これを「チャネルラベル」と言います。
EEGプロセッサでは、DIGITファイルに記録されているチャネルラベルとBDFファイルに記録されているチャネルラベルとを照会することにより、位置情報と測定データとを統合します。DIGITファイルのチャネルラベルとBDFファイルのチャネルラベルとが一致しない場合には、位置座標として無効値が入ります。
2007/08/10現在、このチュートリアルが提供しているデータは、チャネルラベルが不一致のものであるため、成果物に含まれる座標値が正しくありません。とりあえず正しい座標値を入れるためには、Matlabデバッガを使い、次のように操作します。
l vb_read_bdf_save_ch_bin.m(L.97)、inner_rearrange_channel_label関数の呼び出し箇所にブレイクポイントを張る
l Matlabコンソール(Command Window)で次のように入力する
l > ChannelLabel_hardware(1:64) = ChannelLabel_standard’;
処理内容は次の通りです。
1) 初期EEG-MATファイルを読む
2) チャネルデータファイルからEEGデータを抽出する
3) EEG-MATファイルを作成する
テストコード2で実行します。
% eeg_tutorial(2)
内部ではvb_make_msrmnt_data関数が呼ばれます。vb_make_msrmnt_data関数については後述します。
l 初期EEG-MATファイル(EEGinfo構造体をロードして使用する)
l
ch_spec.old_info = 初期EEG-MATファイル(.primary.eeg.mat)からロードしたEEGinfo構造体
./test_dir/test.primary.eeg.mat
1) ch_spec.output_file = 出力されるEEG-MATファイル
2) ch_spec.to_read_swt = [] (default)
3) ch_spec.ch_list = [] (default)
4) ch_spec.time_from = [] (default)
l ch_spec.time_to = [] (default)
l bdf_spec = [] (BDFファイル読み込み処理をスキップする)
ch_specの初期設定では、全チャネルの全サンプルデータを読み出します。
カレントディレクトリ以下の’test_dir’ディレクトリに、以下のデータが作成されます。
l test.eeg.mat
取得したEEGデータから、トリガ情報を抽出、取得することができます。
この処理には大きく次の三つのモードがあります。
1) (自動) 自動でトリガチャネル(トリガ情報出力チャネル)を認識し、トリガポイントを取得する
2) (半自動) ユーザが指定したトリガチャネルからトリガポイントを取得する
3) (手動) トリガチャネルに依らず、ユーザの目視によりトリガポイントを取得する
内部ではvb_mark_trigger関数が呼ばれます。vb_mark_trigger関数については後述します。
一つのトリガ信号は数ミリ秒の幅を持っています。どの時点のサンプルを以て「トリガ」とするかを決めなければいけません。
デフォルトでは「開始サンプル番号」が適用されています。
開始サンプル番号以外には、終了サンプル番号、全サンプル番号が用意されていますが、現在、切替用ユーザインタフェイスを備えていません。変更する場合は、プログラムを一部書き換える(１行だけですが)必要があります。
指定したEEGデータのeeg_data('A1','B2', 'GSR1')をプロットしたfigureが表示されます。
X軸がサンプル、Y軸がEEG値です。
[RANGE]/[LINE]ボタンを[LINE]にするとマウスクリックでトリガをポイントすることができます。トリガラインを再クリックすると削除できます。
[SAVE]ボタンでトリガリストを保存することができます。
処理内容は次の通りです。
4) 指定されたチャネルのEEGデータをプロットする
5) EEGinfoヘッダを読んで、自動でトリガチャネルを認識する
6) トリガチャネルからトリガポイント情報を取得する
7) 取得したトリガポイントにラインを引く
8) トリガリストファイルを作成する(手動)
1)
コード3で実行します。
% eeg_tutorial(3)
2) 指定したEEGデータのeeg_data('A1','B2','GSR1')をプロットしたfigureが表示されます。X軸がサンプル、Y軸がEEG値です。
3) そのグラフの上に、垂直線(X値一定)でトリガの線が引かれます。
4) [SAVE]ボタンをクリックして、トリガ情報を保存します。保存ファイルを’./test_dir/trigger_list.mat’としてください。
l measurement_file = EEG-MATファイル("./test_dir/test.eeg.mat" c.f. 2-3-4.)
l target_list = {'A1', 'B2', 'GSR1'}
l to_read_swt = ON
l output_file = [] (unused)
l auto_info.auto_detect = ON
カレントディレクトリ以下の’test_dir’ディレクトリに、以下のデータが作成されます
trigger_list.mat (トリガリストファイル)
処理内容は次の通りです。
5) 指定されたチャネルのEEGデータをプロットする
6) 指定されたトリガチャネルを読み、トリガポイント情報を取得する
7) 取得したトリガポイントにラインを引く
8) トリガリストファイルを作成する(手動)
コード4で実行します。
% eeg_tutorial(4)
トリガマーク用のGUIが表示されます。
l measurement_file = EEG-MATファイル("./test_dir/test.eeg.mat" c.f. 2-3-4.)
l target_list = {'A1', 'B2', 'GSR1'}
l to_read_swt = ON
l output_file = [] (unused)
l auto_info.auto_detect = OFF
l auto_info.trigger_channel = 'Status'
カレントディレクトリ以下の’test_dir’ディレクトリに、以下のデータが作成されます
(トリガリストファイル)
処理内容は次の通りです。
9) 指定されたチャネルのEEGデータをプロットする
10) 手作業でトリガポイントをクリックし、ラインを引く
11) トリガリストファイルを作成する(手動)
コード5で実行します。
% eeg_tutorial(5)
トリガマーク用のGUIが表示されます。
l measurement_file = EEG-MATファイル("./test_dir/test.eeg.mat" c.f. 2-3-4.)
l target_list = {'A1', 'B2', 'GSR1'}
l to_read_swt = ON
l output_file = [] (unused)
l auto_info.auto_detect = OFF
l auto_info.trigger_channel = []
カレントディレクトリ以下の’test_dir’ディレクトリに、以下のデータが作成されます
(トリガリストファイル)
処理内容は次の通りです。
1) ベースとなるEEGデータを読む
2) トリガポイント、プレトリガ、ポストトリガ情報から、トライアル単位にデータを分割する
3) 分割して取り出したデータを、新規ファイルとしてEEG-MATファイルに保存する
% output : .\test_dir\test_segmented.eeg.mat
%
% --------------------------------------------------------------------- %
コード6で実行します。
% eeg_tutorial(6)
<spec>
: 切り出す幅は trial_specの二つのフィールド
: pre_trigger(トリガ前)
: post_trigger(トリガ後)
: で指定できます。
: 切り出す時間幅を指定できる segmented_ti (time interval)フィールド
: も用意してあります。
l EEG-MATファイル (./test_dir/test.eeg.mat) EEGinfoをロードする
l トリガリストファイル (./test_dir/trigger_list.mat – c.f. 2-4-3-4) trigger_listをロードする
ch_specにtrial_spec構造体を渡します。
l bdf_spec = [] (BDFファイル読み込み処理はスキップする)
l BASE_EEG_FILE = [BASE_DIR filesep TEST_PREFIX '.eeg.mat']
l load(BASE_EEG_FILE, 'EEGinfo')
l ch_spec.old_info = EEGinfo
l ch_spec.output_file = './test_dir/test_segmented.eeg.mat'
l ch_spec.to_read_swt = OFF
l ch_spec.ch_list = []
l ch_spec.time_from = []
l ch_spec.time_to = []
l trial_spec.trigger_list = trigger_sample_list
l trial_spec.pre_trigger = 100
l trial_spec.post_trigger = 500
l trial_spec.segmented_ti = []
l trial_spec.segment_rule = 1
l ch_spec.trial_spec = trial_spec
カレントディレクトリ以下の’test_dir’ディレクトリに、以下のデータが作成されます
test_segmented.eeg.mat
このファイルに保存されている eeg_dataデータは、各チャネル毎に、サンプル数が601(トリガ前100、トリガ後500)、トライアル数186個のデータに分割されています。
vb_make_msrmnt_data関数は、大きく二つの処理を担当します。
1) BDFファイルを読んでチャネルデータファイルを作成する
2) チャネルデータファイルから、指定チャネル、指定サンプルのデータを抽出する
BDFファイルを読む処理の仕様を定義する構造体。これがemptyの場合は、BDFファイルを読む処理がスキップされる。フィールドは次の通り。
· . bdf_file : BDFファイル(.bdf)
l .digit_file : DIGITIZERファイル(.digit.mat)
l .face_file : FACEファイル(.face.mat)
l .device : デバイス文字列、初期値は’BIOSEMI’
l
.primary_path : 初期EEG-MATファイルおよびチャネルデータのパス＋プレフィクス
(例) ‘./test_dir/test’
ならば、./test_dir/test.primary.eeg.mat
ファイルと、./test_dir/test_${チャネルラベル}.ch.datファイルができる
l .label_type : 使用するラベルタイプ、現在使用されていません
l .eeginfo_version : EEGinfo構造体のバージョン、現在、有効活用されていません
チャネルデータファイルを読む処理の仕様を定義する構造体。これがemptyの場合は、チャネルデータファイルを読む処理がスキップされる。フィールドは次の通り。
l .old_info : 基礎となるEEGinfo構造体
l .output_file : 出力EEG-MATファイル
l .to_read_swt : 読取スイッチ、ch_listフィールドで渡されるチャネルリストが読取なのか無視なのかを指定するスイッチ、boolean、trueならばch_listは読取チャネルリスト
l .ch_list : チャネルラベル文字列を羅列したリスト
l .time_from : 読み取る開始サンプル番号
l .time_to : 読み取る終了サンプル番号
l .trial_spec : トライアル分割時に分割仕様を規定する構造体、フィールドは次の通り
· .pre_trigger : トリガ前サンプル数 (A)
· .post_trigger : トリガ後サンプル数 (B)
· .segmented_ti : セグメント時間長 (C)
·
.segment_rule : セグメントルール(1 or 2 or 3)
切り出すデータのサンプル数は(A)(B)(C)の三つのうち二つが与えられていることが必要。
セグメントルールでは、三つが与えられた場合に適用されるルールを指定する。
1) (A)と(B)
2) (A)と(C)
3) (B)と(C)
l MSRMNTfile : 更新されたファイル名
l MSRMNTinfo : 更新された情報構造体(EEGinfo)
トリガを取得する際に用いる関数です。測定データとトリガポイントを表示・操作できるGUIを起動します。
l msrmnt_file : 測定ファイルを指定します、通常はEEG-MAT(.eeg.mat)ファイルです
l target_list : 表示する測定データのチャネルラベルを文字列を羅列して指定します
l to_read_swt : target_listで指定されるチャネルを読み取るのか無視するのかを指定するboolean値です、trueならば読み取ります
l output_file : 現在使用されていません
l auto_info : トリガを自動的にポイントする際の仕様を規定する構造体です、フィールドは次の通りです
· .auto_detect : 自動取得スイッチです、ON or OFF
· .trigger_channel : トリガ情報が出力されたチャネルを指定します。
· .trigger_file : トリガ情報MATファイルを指定します
·
.how : トリガ取得方法を以下の文字列で指定します、指定されない場合は他のフィールド指定状況から判断します
‘FULL’: 自動
‘SEMI’: 半自動
‘FILE’: トリガ情報MATファイル
‘MANUAL’ : 手動
戻り値はありません。