JP3983989B2 - 脳運動機能解析診断装置 - Google Patents
脳運動機能解析診断装置 Download PDFInfo
- Publication number
- JP3983989B2 JP3983989B2 JP2001079251A JP2001079251A JP3983989B2 JP 3983989 B2 JP3983989 B2 JP 3983989B2 JP 2001079251 A JP2001079251 A JP 2001079251A JP 2001079251 A JP2001079251 A JP 2001079251A JP 3983989 B2 JP3983989 B2 JP 3983989B2
- Authority
- JP
- Japan
- Prior art keywords
- brain
- time
- potential
- motor function
- dynamic spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Eye Examination Apparatus (AREA)
Description
【発明の属する技術分野】
本発明は、運動に関連して変化する脳電位の信号成分を脳波データから取り出して、脳運動機能の解析診断を行なう脳運動機能解析診断装置に関する。
【0002】
【従来の技術】
ヒトの随意運動に伴い発生する脳電位の変化は、運動関連脳電位と呼ばれており、この電位は脳波の信号レベルに比べて非常に小さいことから、従来は、被験者に随意運動を反復させ、その際の筋電図による運動開始時点を基準に多数の脳波データを加算平均し、定常的に変動する大きなレベルの周期信号成分をキャンセルすることで求めている。しかし、その平均化のためには、かなり多数のデータを取得しなければならず、かつ、その処理には相当の時間がかかるという問題があった。また、それは脳波の波形の形状から運動に起因する変化の成分を同定しているために、定量的な判断が難しいという問題があった。
【0003】
一方、公開特許公報「特開2000−126148号(発明の名称:脳波データ処理装置及び記録媒体)」には、脳波デジタルデータに対してのみウェーブレット変換を行い、その結果から特徴的な変化を捉えるようにした手法が開示されている。しかし、実際に筋肉が動いたかどうかを確認する手段がないので、運動開始時間を筋電位変動に基づいて決定することができないところに問題がある。
【0004】
【発明が解決しようとする課題】
本発明の課題は、運動関連脳電位の同定に際して、筋電図などから運動開始点を判断してそこから多数の脳波データを収集し、平均化するといった人為的で繁雑な処理の必要をなくすことにある。
【0005】
【課題を解決するための手段】
本発明は、脳波データ中の定常的な周期信号成分は、フィルタリングにより比較的容易に除去することができ、他方、運動関連脳電位のレベルは低くても、変化があれば周波数の変化として検出しやすくなることに着目して、脳波データ中に埋没している運動関連脳電位の微小な信号成分を、ウェーブレット変換により、周波数分布の時間変化の形で容易に観測できるようにするものである。
【0006】
本発明は、被験者であるヒトの頭皮上に付着させた電極により得られた脳電位デジタルデータに対して、ウェーブレット変換を用いた時間周波数成分同定法を適用し、それにより周波数の時間変化を表わすパラメータを得て、運動関連脳電位の信号成分を定量的に同定することを可能にする。本発明による運動関連脳電位の同定では、必要とされる加算の平均回数が従来例と比較して極端に少なくて済むことから、脳運動機能解析診断をより短時間で行うことができる。
【0007】
さらに、本発明では、被験者であるヒトの皮膚上に付着させた電極により得られる単一試行の随意運動に伴って生じる筋電位デジタルデータに対してウェーブレット変換を適用し、周波数の時間変化(ダイナミックスペクトラム)を求め、その広帯域で生じる変化の生起時間を検出することで、運動開始時間を自動的に検出することができる。そのため、本発明では、この筋電位デジタルデータによるダイナミックスペクトラムを利用して検出した運動開始時間を、脳電位時間周波数成分を求める上での基準時間とすることができる。
【0008】
また本発明では、筋電位、脳電位、眼球電位のデジタルデータに対してそれぞれダイナミックスペクトラムを求めることにより、筋電位と脳電位の関係、及び、眼球運動に基づくノイズを同時に検出できるので、より精度の高い解析と、運動関連脳電位の時間周波数成分同定とが可能となる。
【0009】
図1により、本発明の概要を説明する。図1は、本発明による、ヒト運動関連脳電位のウェーブレット変換を用いた時間周波数成分同定機能を内蔵した簡易型ヒト脳運動機能解析診断装置のブロック構成図である。図において、Ag/AgCl(銀/塩化銀)皿電極を用いて、ヒトの頭皮1上の脳電位2、耳朶電位3、皮膚上の筋電位4、顔面上の眼球電位5を取得し、それぞれの信号を脳・筋・眼球電位信号伝送ケーブル6により、ヒト脳運動機能解析診断装置7へ入力する。
【0010】
ヒト脳運動機能解析診断装置7は、増幅・A/D変換器8に送り、デジタルデータに変換する。このデジタルデータをデジタル信号伝送ケーブル9により、デジタルフィルタリング処理部(IIRデジタルフィルタリング)10に送り、電源周波数成分(50Hz/60Hz任意に選択決定)、低周波成分(1.6Hz以下、任意に決定)、高周波成分(30Hz以上、任意に決定)を遮断する。そのフィルタリング処理されたデータをデジタルフィルタリング処理信号伝送ケーブル11により、本発明のウェーブレット変換を用いた時間周波数成分同定システム部12に送り、任意のデータ数(可変)にて脳電位デジタルデータの加算平均、時間周波数成分同定の各処理を自動的に行なう。なお、加算平均、時間周波数同定を行う際に基準となる時間は、単一試行時の筋電位デジタルデータに対してウェーブレット変換を適用することで得られるダイナミックスペクトラムにより求める。なお、時間周波数成分同定システム部12には、統計処理を可能とするために記憶装置が設けられており、多人数の、多試行によって求め得たデータを保存し、そのデータと個人データとの比較を行なうことで、脳運動機能の統計的データに基づく診断を可能とする。
【0011】
【発明の実施の形態】
図2は、図1に示す時間周波数成分同定システム部の1実施の形態における処理フローである。
【0012】
図2において、被験者であるヒトの頭皮から取得された脳電位のデジタルデータと、眼球電位のデジタルデータ、耳朶電位のデジタルデータ、皮膚から取得された筋電位のデジタルデータは、それぞれタイムマーカとともにいったんファイル記憶装置21に格納される。次に加算平均処理部22は、ファイル記憶装置21から脳電位デジタルデータの複数波形分を読み出して、加算平均処理する。加算平均された脳電位デジタルデータは、次にウェーブレット変換処理部23においてウェーブレット変換される。ダイナミックスペクトラム生成処理部24は、ウェーブレット変換結果からダイナミックスペクトラムを生成し、ファイル記憶装置25に格納する。
【0013】
一方、脳電位と同時に皮膚から取得された眼球電位のデジタルデータ、耳朶電位のデジタルデータ、筋電位のデジタルデータは、それぞれウェーブレット変換処理部26でウェーブレット変換される。ダイナミックスペクトラム生成処理部27は、それぞれのウェーブレット変換結果からダイナミックスペクトラムを生成し、ファイル記憶装置25に格納する。変化生起時間検出処理部28は、運動の開始に伴って生起される筋電位の変化はパルス性であり、広帯域の周波数が含まれることから、ダイナミックスペクトラム上で広帯域の周波数成分の同時的変化を検出して、その時間位置を基準時間として出力する。脳運動機能解析処理部29は、その基準時間に基づいて脳電位のダイナミックスペクトラムを解析し、運動関連脳電位を同定する。
【0014】
図1により取得した単一試行による筋電位デジタル波形に対してウェーブレット変換を行なうことにより周波数の時間変化(ダイナミックスペクトラム)が求まり、その筋電位デジタルデータのダイナミックスペクトラムにおいて、広帯域で生じる変化の生起時間を検出することで、運動開始時点を自動的に検出する。そして、その運動開始時点を基準に、運動前2秒、運動後1秒の計3秒を切り出し、その脳電位デジタルデータ20波形分を加算平均した波形に対してウェーブレット変換を行うことにより、ダイナミックスペクトラムが図3のように求めることができる。すなわち運動関連脳電位を時間周波数成分として同定可能となる。この時間周波数成分はデジタルデータとして記憶装置に保存される。
【0015】
なお、同図は国際10−20電極配置法に基づいて配置した電極Czでの脳電位の時間周波数成分同定結果(上図:脳波の20波形分の加算平均結果、下図:ダイナミックスペクトラム)を示しており、ダイナミックスペクトラムは、エネルギーレベルが高い部分を黒色、低い部分を白色と階調表現している。また、他の箇所(C3、C4、Fz、Pz等)でも同様なものが取得できる。データを取得するための測定点数、測定箇所も上記以外にて任意に決定できる。
【0016】
ここで、ウェーブレット変換は、
【0017】
【数1】
【0018】
a:スケール・パラメータ
b:トランスレート・パラメータ
と定義されており、但し、信号f(t)が、脳電位デジタルデータ、筋電位デジタルデータ、眼球電位デジタルデータに相当し、ψ(t)が、マザーウェーブレットである。本発明のウェーブレット変換では、マザーウェーブレットとして、時間と周波数に関する不確定性が最小の関数(時間周波数空間において局在性が良い関数)であることから、〔数2〕に示す信号の周波数を探り出すのに適しているガボール、又は、ガボアー(Gabor )のマザー・ウェーブレットを用いている。
【0019】
【数2】
【0020】
なお、
a=1/2πf
f:周波数〔Hz〕
とし、この周波数fを可変することにより、様々な時間周波数成分を同定することができる。また、σは任意に決定できる。
【0021】
同図から、運動開始時0sから10Hz前後(α波帯域)の周波数成分が強く表れていることを明らかにできることが分かる。他の電極位置においても求め得た。なお、これまでに報告されている運動関連脳電位において、Czでは運動開始時に約0.1sの波形(P−50)が報告されており、すなわち、約10Hz(=1/0.1s)の周波数と概算できることから、本手法において求められたものと一致すると考えられる。また、Fzでは運動開始後約0.5sにおいて、同様に約10Hzの周波数成分がより強く表れており、この成分はN+50を検出したものと考えられる。
【0022】
運動開始前1.0s〜0.5sでは、約5Hz前後の周波数成分が生じ、0.2s前にその成分が減少することがわかる。これは、BP(Beritschaftspotential)、IS(Intermediate Slope)成分に含まれるものと考えられる。
【0023】
また、本発明の脳運動機能解析診断装置により得られたダイナミックスペクトラムの頭皮上分布により、空間−時間周波数マッピングをデジタルデータとして表示することもできる。
【0024】
【発明の効果】
本発明によるウェーブレット変換を用いた時間周波数成分同定機能を内蔵した脳運動機能解析診断装置は、ヒトの随意運動による脳電位成分を周波数の時間変化として同定することができ、従来の報告されている方法では、例えば、1被験者に対して200回、14名の試行を総加算し、平均化するものと比較して、本発明の図2では1被験者のみ20回試行(さらに少ない数にすることも可能。但し、検出精度が低下する。)と、その加算回数が少ないということから、解析、診断を短時間で容易に行うことができる。
【0025】
従来の報告では運動開始時間を筋電位波形の変化により検出するものや、公開特許公報第2000−126148号では、脳波データのみのウェーブレット変換による結果から特徴的な変化を捉えていたが、本発明では、単一試行の随意運動時の筋電位デジタルデータに対してダイナミックスペクトラムを求め、広帯域で生じる変化の生起時間を検出することで、運動開始時間を自動的に検出し、その時間を基準として脳電位の時間周波数成分同定を行うことから、時間基準という点から高い精度を保証することができる。本発明では筋電位、脳電位、眼球電位の各々のダイナミックスペクトラムを個別に求めていることから、筋電位と脳電位の時間周波数成分の関係、及び、眼球運動によるアーチファクトの時間周波数成分をデジタルデータとして同時に検出でき、より精度の高い解析、運動関連脳電位の時間周波数成分同定が可能となる。
【0026】
また本発明装置では、運動関連脳電位の成分を、波形の形状ではなく、周波数と時間にて数値化できることからより客観的な指標を求めることができ、その値をデジタルデータとして表示することができる。なお、加算平均回数は任意に決定できることから、多人数による統計的処理等も可能である。個人の評価を行うにあたり、統計的データが必要であることから、多人数の多数回の試行データを標準データとして本発明の装置内に保存しておき、そのデータと個人との比較が可能となる。
【0027】
さらに本発明装置により得られたダイナミックスペクトラムの頭皮上分布により、空間−時間周波数マッピングをつくることができることから、脳運動機能異常箇所を特定できる。すなわち、空間−時間周波数成分マッピングのパターンの異常による脳運動機能解析診断が可能となる。
【図面の簡単な説明】
【図1】 この発明のヒト運動関連脳電位のウェーブレット変換を用いた時間周波数成分同定法を内蔵した簡易型ヒト脳運動機能解析診断装置の概要構成図である。
【図2】 時間周波数成分同定システム部の1実施の形態における処理フローを示す説明図である。
【図3】 本発明の簡易型ヒト脳運動機能解析診断装置により得られた結果の一部の写真である。
【符号の説明】
1: ヒトの頭表面(頭頂)
2: Ag/AgCl(銀/塩化銀)皿電極(脳電位測定)
3: Ag/AgCl(銀/塩化銀)皿電極(耳朶電位測定、基準用)
4: Ag/AgCl(銀/塩化銀)皿電極(筋電位測定)
5: Ag/AgCl(銀/塩化銀)皿電極(眼球電位測定)
6: 脳・筋電、眼球電位信号伝送ケーブル
7: 簡易型ヒト脳運動機能解析診断装置
8: 増幅・A/D変換器
9: デジタル信号伝送ケーブル
10: デジタルフィルタリング処理部
11: デジタルフィルタリング処理信号伝送ケーブル
12: ウェーブレット変換を用いた時間周波数成分同定システム部
Claims (6)
- 被験者から取得される脳電位データに対してウェーブレット変換を行なうことにより、周波数成分の時間変化を表わすダイナミックスペクトラムを生成して運動関連脳電位の成分同定を行なう脳運動機能解析診断装置において、
上記被験者から同時に取得される筋電位データにウェーブレット変換を行なうことによりダイナミックスペクトラムを生成し、そのダイナミックスペクトラム上での広帯域の周波数成分の変化に基づいて運動開始時間を検出し、検出された運動開始時間を基準にして運動関連脳電位の成分同定を行なう手段を備えていることを特徴とする脳運動機能解析診断装置。 - 請求項1において、ウェーブレット変換のマザーウェーブレットとして、時間周波数空間において局在性が良い関数を用いることを特徴とする脳運動機能解析診断装置。
- 請求項2において、ウェーブレット変換のマザーウェーブレットとして、時間周波数空間において局在性が良い関数は、ガボールまたはガボアーの関数であることを特徴とする脳運動機能解析診断装置。
- 請求項1ないし請求項3のいずれかにおいて、被験者から取得される筋電位データは、皮膚上に付着させた電極により得られた単一試行の随意運動による筋電位データであることを特徴とする脳運動機能解析診断装置。
- 請求項1において、頭皮上、顔面上、皮膚上に付着された電極により得られた脳電位、眼球電位、筋電位のデータに対してそれぞれウェーブレット変換を行なって、ダイナミックスペクトラムを個別に、かつ、同時に求めて、運動に関連して変化する脳電位信号の時間周波数成分を、筋電位のダイナミックスペクトラムにより得られた運動開始時間を基準として脳波データから取り出すとともに、眼球運動によるアーチファクトの検出と削除を行なうことを特徴とする脳運動機能解析診断装置。
- 被験者から取得される脳電位データに対してウェーブレット変換を行なうことにより、周波数成分の時間変化を表わすダイナミックスペクトラムを生成して運動関連脳電位の成分同定を行なう脳運動機能解析診断装置において、
上記被験者から同時に取得される筋電位データにウェーブレット変換を行なうことによりダイナミックスペクトラムを生成し、そのダイナミックスペクトラム上での広帯域の周波数成分の変化に基づいて運動開始時間を検出する手段を備えていることを特徴とする脳運動機能解析診断装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001079251A JP3983989B2 (ja) | 2001-03-19 | 2001-03-19 | 脳運動機能解析診断装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001079251A JP3983989B2 (ja) | 2001-03-19 | 2001-03-19 | 脳運動機能解析診断装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002272692A JP2002272692A (ja) | 2002-09-24 |
JP3983989B2 true JP3983989B2 (ja) | 2007-09-26 |
Family
ID=18935735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001079251A Expired - Lifetime JP3983989B2 (ja) | 2001-03-19 | 2001-03-19 | 脳運動機能解析診断装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3983989B2 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102783946A (zh) * | 2012-08-20 | 2012-11-21 | 湖州康普医疗器械科技有限公司 | 一种自动脑源定位方法及装置 |
CN102824172A (zh) * | 2012-08-28 | 2012-12-19 | 浙江大学 | 一种脑电特征提取方法 |
US8639320B2 (en) | 2010-06-16 | 2014-01-28 | Sony Corporation | Muscle-activity diagnosis apparatus, method, and program |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10312519A1 (de) * | 2003-03-20 | 2004-10-28 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Verfahren zum Auslösen insassenunterstützter Maßnahmen in einem Fahrzeug |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3156777B2 (ja) * | 1998-10-28 | 2001-04-16 | 日本電気株式会社 | 脳波データ処理装置及び記録媒体 |
-
2001
- 2001-03-19 JP JP2001079251A patent/JP3983989B2/ja not_active Expired - Lifetime
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8639320B2 (en) | 2010-06-16 | 2014-01-28 | Sony Corporation | Muscle-activity diagnosis apparatus, method, and program |
CN102783946A (zh) * | 2012-08-20 | 2012-11-21 | 湖州康普医疗器械科技有限公司 | 一种自动脑源定位方法及装置 |
CN102783946B (zh) * | 2012-08-20 | 2014-06-25 | 湖州康普医疗器械科技有限公司 | 一种自动脑源定位方法及装置 |
CN102824172A (zh) * | 2012-08-28 | 2012-12-19 | 浙江大学 | 一种脑电特征提取方法 |
CN102824172B (zh) * | 2012-08-28 | 2014-03-19 | 浙江大学 | 一种脑电特征提取方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2002272692A (ja) | 2002-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7672717B1 (en) | Method and system for the denoising of large-amplitude artifacts in electrograms using time-frequency transforms | |
Fatourechi et al. | EMG and EOG artifacts in brain computer interface systems: A survey | |
Mak et al. | An automated ECG-artifact removal method for trunk muscle surface EMG recordings | |
US20070179396A1 (en) | Method and System for Detecting and Classifying Facial Muscle Movements | |
Sawangjai et al. | EEGANet: Removal of ocular artifacts from the EEG signal using generative adversarial networks | |
Sakai et al. | Data augmentation methods for machine-learning-based classification of bio-signals | |
WO2014150684A1 (en) | Artifact as a feature in neuro diagnostics | |
CA3101371A1 (en) | Quantifying motor function using eeg signals | |
Taelman et al. | Wavelet-independent component analysis to remove electrocardiography contamination in surface electromyography | |
Talsma et al. | Methods for the estimation and removal of artifacts and overlap in ERP waveforms | |
Zou et al. | Automatic EEG artifact removal based on ICA and Hierarchical Clustering | |
CN108451494B (zh) | 使用瞳孔反应来检测时域心脏参数的方法及系统 | |
JP3983989B2 (ja) | 脳運動機能解析診断装置 | |
Haddix et al. | A comparison of eeg alpha rhythm detection by tripolar concentric ring electrodes and conventional disk electro des | |
Heravi et al. | Recurrence quantification analysis of electrooculography signal to a control question test: A new approach for the detection of deception | |
JP3661512B2 (ja) | 脳波信号処理装置及び方法 | |
Duhan et al. | Study of signal processing techniques for EMG analysis | |
JP2006110234A (ja) | 脳活動情報データ処理方法およびその装置 | |
Kiamini et al. | A wavelet based algorithm for ocular artifact detection in the EEG signals | |
Kukreti | EEG-based Stress Detection using DEAP Dataset and Support Vector Machines | |
Egambaram et al. | Automated and online eye blink artifact removal from electroencephalogram | |
JP3724524B2 (ja) | 眼球制御系情報検出装置および眼球制御系の解析方法 | |
JP2004350797A (ja) | 脳波の評価方法及びその評価装置 | |
JP3799889B2 (ja) | 瞬目動作解析装置及び方法 | |
Arjunan et al. | Changes in decibel scale wavelength properties of EEG with alertness levels while performing sustained attention tasks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20040210 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050826 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070109 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070312 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070410 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070611 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070703 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070705 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100713 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |