発表概要 背景 目的 研究テーマ 1 脳電図 (EEG) および脳磁図 (MEG) を用いた脳活動計測 システムの開発 2 経頭蓋磁気刺激 (TMS) を用いた脳標的刺激システムの開発 3 神経計算モデルを用いた計算機シミュレーション 研究成果

Similar documents
ヒト脳機能データ推定システムの研究開発

今日勉強すること 1. 反射弓と伸張反射 2. 屈曲反射 3. 膝蓋腱反射の調節機構 4. 大脳皮質運動野の機能

3T MRI導入に伴う 安全規程の変更について

2. ICA ICA () (Blind Source Separation BBS) 2) Fig. 1 Model of Optical Topography. ( ) ICA 2.2 ICA ICA 3) n 1 1 x 1 (t) 2 x 2 (t) n x(t) 1 x(t

末梢神経障害

M波H波解説

生物 第39講~第47講 テキスト

インターリーブADCでのタイミングスキュー影響のデジタル補正技術

CMOS リニアイメージセンサ用駆動回路 C CMOS リニアイメージセンサ S 等用 C は当社製 CMOSリニアイメージセンサ S 等用に開発された駆動回路です USB 2.0インターフェースを用いて C と PCを接続

周波数特性解析

バイバルコロナリーステント 2015 年 1 月作成第 1 版本ステントは 非臨床試験において 条件付きで MRI 検査の危険性がない MR Conditional に該当することが立証されている 下記条件にて留置直後から MRI 検査を安全に施行することができる 静磁場強度 3 テスラ以下 空間勾

フォト IC ダイオード S SB S CT 視感度に近い分光感度特性 視感度特性に近い分光感度特性をもったフォトICダイオードです チップ上には2つの受光部があり 1つは信号検出用受光部 もう1つは近赤外域にのみ感度をもつ補正用受光部になっています 電流アンプ回路中で2

脳組織傷害時におけるミクログリア形態変化および機能 Title変化に関する培養脳組織切片を用いた研究 ( Abstract_ 要旨 ) Author(s) 岡村, 敏行 Citation Kyoto University ( 京都大学 ) Issue Date URL http

Presentation Title Arial 28pt Bold Agilent Blue

Ø Ø Ø

16K14278 研究成果報告書

fMRIについて


CCD リニアイメージセンサ用駆動回路 C CCD リニアイメージセンサ (S11155/S ) 用 C は 当社製 CCDリニアイメージセンサ S11155/S 用に開発された駆動回路です S11155/S11156-

CMOS リニアイメージセンサ用駆動回路 C10808 シリーズ 蓄積時間の可変機能付き 高精度駆動回路 C10808 シリーズは 電流出力タイプ CMOS リニアイメージセンサ S10111~S10114 シリーズ S10121~S10124 シリーズ (-01) 用に設計された駆動回路です セン

情報処理学会研究報告 IPSJ SIG Technical Report Vol.2016-CE-136 No /10/15 Arduino 1,a) 1,b) 1,c) 1,d) Arduino Arduino Kengo Shimizu 1,a) Motomasa Omura 1,b

Microsoft PowerPoint - machida0206

電気的特性 (Ta=25 C) 項目 記号 条件 Min. Typ. Max. 単位 読み出し周波数 * 3 fop khz ラインレート * Hz 変換ゲイン Gc ゲイン =2-5 - e-/adu トリガ出力電圧 Highレベル Vdd V -

Microsoft PowerPoint - 集積回路工学(5)_ pptm


W 1983 W ± Z cm 10 cm 50 MeV TAC - ADC ADC [ (µs)] = [] (2.08 ± 0.36) 10 6 s 3 χ µ + µ 8 = (1.20 ± 0.1) 10 5 (Ge

Microsoft PowerPoint - 9.Analog.ppt

報道発表資料 2008 年 7 月 17 日 独立行政法人理化学研究所 単語やメロディーの切れ目に対応する脳活動の記録に成功 - 分節化進行過程の神経活動を 世界で初めて生理学的手法で観察 - ポイント 連続音声に含まれる単語やメロディーの切れ目だけに出現する脳波を発見 脳波の強さは音声分節化と統計

DVIOUT

<4D F736F F F696E74202D2082B182C782E082CC8AC888D58C7691AA964082C982E682E B092698A4B94BB92E FC92E E707074>

スライド 1

形式 :WYPD 絶縁 2 出力計装用変換器 W UNIT シリーズ パルスアイソレータ ( センサ用電源付 2 出力形 ) 主な機能と特長 パルス入力信号を絶縁して各種のパルス出力信号に変換 オープンコレクタ 電圧パルス リレー接点パルス出力を用意 センサ用電源内蔵 耐電圧 2000V AC 密着

理学療法学43_supplement 1

INTRODUCTION 各人が被験者となり 脳波の測定を実際に経験することで 測定方法や 得られたデータの評価方 法などを学ぶ METHODS 頭部の皮膚の所定の部位をアルコールでよく拭き ペーストをつけて電極を接着する 電極の位置は国際 法に従う 増幅器の時定数は 0.3 秒にする

. ) ) ) 4) ON DC 6 µm DC [4]. 8 NaPiOn 4

推奨条件 / 絶対最大定格 ( 指定のない場合は Ta=25 C) 消費電流絶対最大定格電源電圧 Icc 容量性負荷出力抵抗型名 Vcc Max. CL 電源電圧動作温度保存温度 Zo (V) 暗状態 Min. Vcc max Topr* 2 Tstg* 2 Min. Max. (ma) (pf)

あなたは磁界の中にいる あなたに磁力が襲いかかる Z BA 自由部門 :20059

Microsoft PowerPoint - 14.ppt [互換モード]

PEA_24回実装学会a.ppt

形式 :PDU 計装用プラグイン形変換器 M UNIT シリーズ パルス分周変換器 ( レンジ可変形 ) 主な機能と特長 パルス入力信号を分周 絶縁して単位パルス出力信号に変換 センサ用電源内蔵 パルス分周比は前面のスイッチで可変 出力は均等パルス オープンコレクタ 電圧パルス リレー接点パルス出力

Microsoft PowerPoint - 01_内田 先生.pptx

アクティブフィルタ テスト容易化設計

スライド 1

< B837B B835E82C982A882AF82E991CF905593AE90AB8CFC8FE382C98AD682B782E988EA8D6C8E40>

<4D F736F F F696E74202D2088E B691CC8C7691AA F C82512E B8CDD8AB B83685D>

レベルシフト回路の作成

第 5 章復調回路 古橋武 5.1 組み立て 5.2 理論 ダイオードの特性と復調波形 バイアス回路と復調波形 復調回路 (II) 5.3 倍電圧検波回路 倍電圧検波回路 (I) バイアス回路付き倍電圧検波回路 本稿の Web ページ ht

スライド 1

モジュール式アナログアンプ 形式 VT-MSPA1-1 VT-MSPA1-10 VT-MSPA1-11 RJ 形式 : 改訂 : シリーズ 1X H6833_d 特長 内容 電磁比例圧力弁の制御に適しています : DBET-6X DBEM...-7X (Z)D

特長 01 裏面入射型 S12362/S12363 シリーズは 裏面入射型構造を採用したフォトダイオードアレイです 構造上デリケートなボンディングワイヤを使用せず フォトダイオードアレイの出力端子と基板電極をバンプボンディングによって直接接続しています これによって 基板の配線は基板内部に納められて

フォトダイオードモジュール C10439 シリーズ 精密測光用フォトダイオードと低ノイズアンプを一体化 フォトダイオードモジュール C10439 シリーズは フォトダイオードと I/V アンプを一体化した高精度な光検出器です アナログ電圧出力のため 電圧計などで簡単に信号を観測することができます ま

形式 :RPPD 計装用プラグイン形変換器 M UNIT シリーズ パルスアイソレータ ( センサ用電源付 ロータリエンコーダ用 ) 主な機能と特長 ロータリエンコーダの 2 相パルス入力信号を絶縁して各種の 2 相パルス出力信号に変換 オープンコレクタ 電圧パルス パワーフォト MOS リレー R

NCB564個別00版

Microsoft Word - NJJ-105の平均波処理について_改_OK.doc

Microsoft PowerPoint _kawashima_DL

研究の背景これまで, アルペンスキー競技の競技者にかかる空気抵抗 ( 抗力 ) に関する研究では, 実際のレーサーを対象に実験風洞 (Wind tunnel) を用いて, 滑走フォームと空気抵抗の関係や, スーツを含むスキー用具のデザインが検討されてきました. しかし, 風洞を用いた実験では, レー

parameter SAR SAR mode parameter SAR mode Low Moderate High 3 2W/kg 3W/kg 4W/ kg Low SAR SAR parameter parameter SAR parameter parameter SAR W/k

スライド 1

サカナに逃げろ!と指令する神経細胞の分子メカニズムを解明 -個性的な神経細胞のでき方の理解につながり,難聴治療の創薬標的への応用に期待-

科目名授業方法単位 / 時間数必修 選択担当教員 人体の構造と機能 Ⅱ 演習 2 単位 /60 時間必修 江連和久 北村邦男 村田栄子 科目の目標 人体の構造と機能 はヒトの体が正常ではどうできていてどう働くのかを理解することを目的とする この学問は将来 看護師として 病む ということに向き合う際の

PowerPoint プレゼンテーション

NJU72501 チャージポンプ内蔵 圧電用スイッチングドライバ 概要 NJU72501はチャージポンプ回路を内蔵し 最大で3V 入力から 18Vppで圧電サウンダを駆動することができます このチャージポンプ回路には1 倍 2 倍 3 倍昇圧切り替え機能を備えており 圧電サウンダの音量を変更すること

- (20 ) 400 () 3DCG No.51 No.61 No.62 No.11 P 2 16

パルス波高値計測回路の製作

Microsoft PowerPoint - semi_ppt07.ppt [互換モード]

Transcription:

SCOPE 若手 ICT 研究者育成型研究開発 (082108003) 脳の感覚運動変換メカニズムの解明に関する研究開発 研究代表者 福田浩士 広島市立大学大学院情報科学研究科

発表概要 背景 目的 研究テーマ 1 脳電図 (EEG) および脳磁図 (MEG) を用いた脳活動計測 システムの開発 2 経頭蓋磁気刺激 (TMS) を用いた脳標的刺激システムの開発 3 神経計算モデルを用いた計算機シミュレーション 研究成果

脳の感覚運動変換メカニズムの解明に関する研究開発 背景 目的

ブレイン コンピュータ インタフェース (BCI) 脳 検出 処理 脳内の情報表現を解読 制御 環境 脳活動の発現機構 BCI 脳内情報の意味 検出 処理 脳内の情報表現に変換 制御 コンピュータ BCIの創出には脳の情報処理メカニズムの解明が必須新しい情報通信技術の創出が期待できる

BCI への応用を見据えた 脳計測 解析システム 脳電図 (EEG) 脳内の情報表現を解読 脳活動操作 解析 脳活動計測 解析 脳磁図 (MEG) 脳活動の発現機構脳内情報の意味 脳内の情報表現に変換 経頭蓋磁気刺激 (TMS)

感覚運動変換に関連する領野 補足運動野 : 運動の開始運動の順序 運動前野 : 運動の計画 運動野 : 運動の実行 体性感覚野 : 運動感覚の処理 頭頂連合野 : 感覚と運動の統合 視覚野 : 視覚情報の処理 これらの領野がどのようなタイミングでどのように情報をやり取りすることで感覚運動変換が行われるのか?

目的 ヒトの視覚誘導性到達運動における感覚運動変換に焦点を絞り 1 高密度脳電図 (EEG), 脳磁図 (MEG) を用いた脳活動計測 解析システムの開発と運動計測実験 2 経頭蓋磁気刺激 (TMS) を用いた脳標的刺激システムの開発と運動計測実験 3 神経計算モデルによる計算機シミュレーションを通して, そのメカニズムの解明にアプローチ

脳の感覚運動変換メカニズムの解明に関する研究開発 研究テーマ 1 脳電図 (EEG) および脳磁図 (MEG) を用いた脳活動計測 システムの開発と運動計測実験 1. MEG を用いた脳の活動部位推定法の開発 2. EEG を用いた運動タスクの集中度に関連した脳活動の計測 解析

脳の計測技術 脳の構造を計測 X 線 CT 核磁気共鳴画像法 (MRI) など 脳の活動を計測 脳電図 (EEG) 脳磁図 (MEG) 機能的 MRI(fMRI) 近赤外分光法 (NIRS) など

脳活動計測技術の比較 空間分解能時間分解能 EEG MEG fmri NIRS? 数十 mm < 1 ms? 2~10 mm < 1 ms 1~3 mm > 数十 ms 携帯性 安全性 8 mm 以上 > 数十 ms 脳電図 (EEG) と脳磁図 (MEG) の特長 時間分解能に優れている 空間分解能は解析法に依存する

1-1. 脳磁図 (MEG) を用いた 脳の活動部位推定法の開発 概要 大脳皮質の錐体細胞の方向を考慮した脳の活動部位推定法を開発 体性感覚誘発脳磁界 (SEF) と聴覚誘発脳磁界 (AEF) を用いてそれぞれの発生源を推定

全頭型 MEG の計測 磁気センサ 神経活動電流 脳磁界が発生 磁気センサを頭部全体に配置して 磁界の時間変化 磁界の空間分布を計測

大脳皮質の構造 頭部の前額断面 脳表 錐体細胞 電流 主な活動源は大脳皮質の錐体細胞 錐体細胞内の電流の方向は軸索の方向

提案法における処理の流れ 頭部の MR 画像 脳表の切り出し 脳の 3D 再構成 法線方向の計算 活動部位の表示 電流源推定

両手刺激 SEF の発生源の推定結果 L R A P 計測波形 活動部位の推定結果 ( 左半球 ) L R A P N20m 成分の空間分布 活動部位の推定結果 ( 右半球 )

両耳刺激 AEF の発生源の推定結果 L R A P 計測した磁界 活動部位の推定結果 ( 左半球 ) L R A P N100m 成分の空間分布 活動部位の推定結果 ( 右半球 )

テーマ 1 1 のまとめ 体性感覚誘発磁界 (SEF) の発生源推定結果 体性感覚野の手首領域に局限 聴覚誘発磁界 (AEF) の発生源推定結果 聴覚野に局限 同一被験者での推定結果のばらつき 目標であった10 mm 以内を達成

1-2.EEG を用いた運動タスクの集 中度に関連した脳活動の計測 解析 概要 EEG を用いて運動タスクの集中度に関連した脳活動を計測 集中度が高い時に運動野, 補足運動野, 運動前野, 頭頂連合野がうまく連携している可能性を示唆する脳活動を検出

集中度と運動関連領野の情報処理の関係についての研究 両手指交互動作における相転移 運動の周波数が高くなる 無意識の同時動作が発生 集中度高 集中度低 脳活動の違いを EEG で捉えられるか?

実験方法 256 チャネル EEG センサ 運動検出器 交互 (A) モード T T Right Left T R T R T L T L 同時 (S) モード 150 ms 150 ms トリガ処理 脳電計 (EEG) 右利きの被験者 両手の示指をできるだけ速く交互にタッピング (4 分間 ) トリガ信号 : 各人差し指の伸展開始 交互モード ( 集中度高 ) と同時モード ( 集中度低 ) で脳活動を比較

EEGの計測 計測対象の大脳皮質の領野 補足運動野 運動開始 運動順序 運動前野 運動計画 FC3 FC4 FCz 運動野 運動実行 C3 体性感覚野 運動感覚 CP3 C4 CP4 CPz 頭頂連合野 感覚と運動の統合 P3 大脳皮質左半球 右指の制御 P4 大脳皮質右半球 左指の制御

電極間コヒーレンスの解析結果

EEGの解析結果 交互モードで脳活動の関連が強い領野 補足運動野 運動開始 運動順序 運動前野 運動計画 FC3 FC4 FCz 運動野 運動実行 C3 体性感覚野 運動感覚 CP3 C4 CP4 CPz 頭頂連合野 感覚と運動の統合 P3 P4 領野間の関連を定量化 大脳皮質左半球 大脳皮質右半球 集中度を判定可能 右指の制御 左指の制御

EEG に関する今後の研究開発 EEGを用いた集中度のリアルタイム判定システムの開発 応用例 教育学習時の集中度の判定 運転時の集中度の判定 模擬運転時の脳活動計測

研究テーマ 1 のまとめ EEGおよびMEGを用いた脳活動計測 解析システムを開発した 推定精度 10 mmを概ね達成 運動の集中度に関する脳活動を検出した 集中度判定を行うBCIの創出につながる

脳の感覚運動変換メカニズムの解明に関する研究開発 研究テーマ 2 経頭蓋磁気刺激 (TMS) を用いた脳標的刺激システムの開発と運動計測実験

運動野の活動と筋電位 (EMG) M1 が活動 第 1 次運動野 (M1) 脊髄 運動指令が下降 M1 の体部位局在性 EMG 筋活動に応じた筋電位が発生

経頭蓋磁気刺激 (Transcranial Magnetic Stimulation: TMS) 8 の字コイル パルス磁界 パルス電流 誘導電界 非侵襲で脳への刺激が可能

運動誘発電位 (Motor Evoked Potential: MEP) 8 の字コイル M1 を刺激 第 1 次運動野 (M1) 脊髄 遠心性信号 1 mv 20 ms M1 の体部位局在性 MEP

TMS の刺激部位推定法 8 の字コイル パルス磁界 パルス電流 誘導電界 脳は複雑な構造をしているため中心直下で刺激が最大となるとは限らない 脳の構造を考慮した刺激部位推定法が必要

大脳皮質の構造 頭部の前額断面 脳表 錐体細胞 電流 主な活動源は大脳皮質の錐体細胞 錐体細胞内の電流の方向は軸索の方向

神経磁気刺激における神経興奮特性 脳表 誘導電界 神経の等価回路 脳表 細胞体 神経細胞 軸索 dv Gm dt n Cm Vr 1 G m V r Ga C G m m V r Ga L Vn Vn+ 1 V n+ 2 C m Ga En En+ 1 E n+ 2 [ G ( V + V LE ) G V ] = a n n+1 Cm n m n 細胞体の末端 軸索 実質的な刺激強度 : 誘導電界の神経長軸方向成分

TMS の刺激部位推定法 刺激コイル 刺激部位 脳表 刺激コイルに流れるパルス電流による誘導電界 錐体細胞 錐体細胞 錐体細胞 誘導電界の神経長軸方向成分が最大となる部位を刺激部位とする 誘導電界の神経長軸方向成分 ( 脳表の法線方向 )

TMS の刺激部位推定法 頭部の MR 画像 脳表の切り出し 脳の 3D 再構成 法線方向の計算 刺激部位推定 刺激部位の表示

運動誘発電位 (Motor Evoked Potential: MEP) 8 の字コイル M1 を刺激 第 1 次運動野 (M1) MEP が観測される筋と TMS で刺激された M1 の部位が一致するかを検証 脊髄 遠心性信号 1 mv 20 ms M1 の体部位局在性 MEP

刺激部位と MEP の対応 上腕ターゲット 短拇指外転筋 上腕二頭筋 掌ターゲット 短拇指外転筋 上腕二頭筋

研究テーマ 2 のまとめ TMSを用いた脳標的刺激システムを開発した 目標としていた精度 5 mm 以内を概ね達成

脳の感覚運動変換メカニズムの解明に関する研究開発 研究テーマ 3 神経計算モデルを用いた計算機シミュレーション 神経回路モデルを用いた運動誘発電位 (MEP) からの脳刺激部位の推定法の開発

運動野の活動と筋電位 (EMG) M1 が活動 第 1 次運動野 (M1) 外部から脳を操作することで 筋活動を誘発できないか脊髄? 運動指令が下降 M1 の体部位局在性 EMG 筋活動に応じた筋電位が発生

運動誘発電位 (Motor Evoked Potential: MEP) 8 の字 M1の各部位が表現している筋活動中心溝コイル (M1における運動表現) の解明にアプローチできる このためには少なくとも 1 刺激部位を正確に推定する 2MEPを刺激部位に対応づける脊髄第 1 次運動野ことが必要になる (M1) M1 を刺激 遠心性信号 1 mv 20 ms MEP

MEP のばらつき 8 の字コイル M1を刺激 1 mv 20 ms MEP 刺激コイルと頭部を固定しても MEP の振幅がトライアルごとにばらつく MEP のばらつきは不可避である (Kiers et al., 1993) ばらつきは脊髄のレベルではなく大脳皮質のレベルの運動ニューロンの興奮特性によるもの (Ellaway et al., 1998) ばらつきの成因は解明されていない 観測された MEP から刺激部位を特定するのは困難

目的 トライアルごとにばらつく MEP を TMS の刺激部位に対応づける MEP 刺激部位 1 mv 20 ms 神経回路モデルを用いて,MEP から刺激部位を判別することを試みる

実験装置 :TMS Delayed Trigger Pulse Ag AgCl Electrodes Chin Rest Stimulating Coil Stimulator 1 sec Delay 刺激コイル ヘッドレスト Head Rest Trigger Pulse Trigger Generator EMG AMP. FDI LUM FPI FDS EIP EDC Force Sensor AMP. f z A/D Convertor USB PC f xf y 多関節アーム 磁気刺激装置 Magstim 200 (Magstim 社 ) 最大出力 : 2.2 T 磁気パルス : duration 1 ms rise time 100 μs TMS の強度 : 最大強度の 70% 刺激コイル 8 の字コイル (70 mm) 7 自由度多関節アーム刺激コイルを任意の位置と向きに固定することができる 顎台とヘッドレスト被験者の頭部を固定 顎台

実験装置 : 表面筋電図 (EMG) の計測 Delayed Trigger Pulse Chin Rest Stimulating Coil Stimulator Ag AgCl Electrodes 1 sec Delay Head Rest Trigger Pulse Trigger Generator EMG AMP. FDI LUM FPI FDS EIP EDC Force Sensor AMP. f z A/D Convertor USB PC f xf y 筋電アンプ Personal EMG ( 追坂電子 ) 倍率 : 1000 電極 Ag AgCl 電極 A/D 変換器 USB6216 (National Instruments 社 ) サンプリング周波数 : 5 khz 計測対象の筋 ( 右手 右前腕 ) 1. 背側骨間筋 (FDI) 2. 虫様筋 (LUM) ( 右手背側 ) ( 右手掌側 ) 3. 掌側骨間 4. 浅指屈筋筋 (FPI) (FDS) (右腕橈側)( 右手掌側 ) 6. 総指伸筋 5. 示指伸筋 (右腕尺側)(右腕尺側)(EIP) (EDC)

実験装置 : 指力の計測 Delayed Trigger Pulse Stimulator Ag AgCl Electrodes Chin Rest Stimulating Coil 1 sec Delay Head Rest Trigger Pulse Trigger Generator EMG AMP. FDI LUM FPI FDS EIP EDC Force Sensor AMP. f z A/D Convertor USB PC f xf y 3 軸力覚センサ TFS 12 25( ニッタ ) 示指先の 3 軸方向の力を計測 A/D 変換器 USB6216 (National Instruments 社 ) サンプリング周波数 : 5 khz ディジタルフィルタ 100 Hz ローパスフィルタ 3 軸力覚センサ F y F z F x

実験結果 : 刺激部位推定 すべてのコイル位置で左側の M1 が刺激されていたことを確認

実験結果 :MEP 被験者 1 Superior Inferior すべての刺激部位で MEP がトライアルごとにばらつくことを確認

実験結果 : 指力 被験者 1 Superior Inferior すべての刺激部位で指力がトライアルごとにばらつくことを確認 MEP と比較してばらつきが大きい

神経回路モデルの構造

判別結果 Force NN EMG NN Force EMG NN 各被験者 EMG NN: 80 % 以上 Force EMG NN : 判別能力が向上 Force NN : 判別能力が低い 全被験者 EMG NN: 89.3 % Force EMG NN : 93.3 % Force NN : 73.7 % ばらつく成因が解明されていない運動誘発電位 (MEP) を, 神経回路モデルを用 MEPを用いた場合に高い判別能力を示した いることで刺激部位に対応づけることが正判別率 可能

研究テーマ 3 のまとめ 運動誘発電位 (MEP) から脳の刺激部位を推定する神経回路モデルを開発した 推定精度 90% 以上を達成した 脳操作型 BCIの基礎技術として期待できる

研究成果 H20 H21 H22 合計 査読付き論文数 5 件 (4 件 ) 4 件 (4 件 ) 2 件 (0 件 ) 11 件 (8 件 ) 口頭発表数 9 件 (6 件 ) 10 件 (2 件 ) 8 件 (1 件 ) 27 件 (9 件 ) 受賞数 1 件 (0 件 ) 1 件 (0 件 ) 0 件 (0 件 ) 2 件 (0 件 ) 査読付き論文 1. Hiroshi Fukuda et al., "Estimation Method on Multiple Sources of MEG Based on the Columnar Structure of the Cerebral Cortex", IEEE Transactions on Magnetics, Vol. 44, Issue 11, pp. 4425-4432 (Nov. 2008) 2. Hiroshi Fukuda et al., Identification of Stimulated Sites Using Artificial Neural Networks Based on Transcranial Magnetic Stimulation-Elicited Motor Evoked Potentials", 電気学会論文誌 C Vol.120 No.12 pp2204-2209(2010 年 12 月 1 日 ) など 受賞 1. 生体医工学シンポジウム2008ベストリサーチアワード (2008 年 9 月 19 日 ) 2. Young Investigator's Award of IEEE West Japan Chapter of Engineering in Medicine and Biology (Nov. 28, 2009)