JP4236352B2 - 生体信号計測装置 - Google Patents
生体信号計測装置 Download PDFInfo
- Publication number
- JP4236352B2 JP4236352B2 JP30413299A JP30413299A JP4236352B2 JP 4236352 B2 JP4236352 B2 JP 4236352B2 JP 30413299 A JP30413299 A JP 30413299A JP 30413299 A JP30413299 A JP 30413299A JP 4236352 B2 JP4236352 B2 JP 4236352B2
- Authority
- JP
- Japan
- Prior art keywords
- component
- signal
- independent
- noise
- detection signal
- 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 - Fee Related
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Description
【０００１】
【発明の属する技術分野】
この発明は、脳磁計や脳波計などの生体信号計測装置に係り、特に生体信号測定用のセンサにより得られた検出信号からノイズ成分を除去し、脳内活動を推定するための技術に関する。
【０００２】
【従来の技術】
生体内に流れる生体活動電流により微小な生体磁気（生体磁界）が生体から発生する。例えば、脳から発生する生体磁気は脳磁と呼ばれ、生体に刺激を与えることにより発生する誘発脳磁や、α波やてんかんのスパイク波のように自然に発生する自発脳磁などがある。
【０００３】
近年、生体から出る微小な生体磁気を測定できる磁束計として、ＳＱＵＩＤ（Superconducing Quantum Interference Device：超電導量子干渉計）を用いたマルチチャンネルＳＱＵＩＤセンサが開発されている。マルチチャンネルＳＱＵＩＤセンサは、デュアーと呼ばれる容器内に多数個のＳＱＵＩＤセンサを液体窒素などの冷媒に浸漬・収納した構成となっている。
【０００４】
このマルチチャンネルＳＱＵＩＤセンサ（以下、適宜「磁束計」と略記）を備えた生体信号計測装置、つまり生体磁気計測装置の場合、磁束計を被検体の関心部位である例えば頭部の傍らに置くと、頭部内に生じた生体活動電流源から発生する微小な生体磁気が磁束計内の各ＳＱＵＩＤセンサで無侵襲で測定されて磁気検出信号として出力されるとともに、ＳＱＵＩＤセンサからの磁気検出信号に基づいて生体磁気解析が行われて、生体活動電流源の状態（例えば、位置や向き或いは大きさ等）を把握することができるという構成になっている（例えば特開平７−３２７９４３号公報参照）。
【０００５】
【発明が解決しようとする課題】
しかしながら、従来の生体磁気計測装置には、各ＳＱＵＩＤセンサによって検出された磁気検出信号に含まれるノイズ成分を十分に除去することが難しいという問題がある。測定対象である生体磁気は非常に微弱であることから生体活動電流源以外の別の磁気発生源から出る（環境ノイズとも称するような）ノイズ磁気の混入が避けられない。したがって、各ＳＱＵＩＤセンサからの磁気検出信号には混入するノイズ磁気によるノイズ成分が含まれており、磁気検出信号からノイズ成分を十分に除去しなければ、生体磁気の正確な解析はおぼつかないことになる。
【０００６】
そこで、生体磁気測定用のＳＱＵＩＤセンサとは別のノイズ磁気検出専用の磁気センサでノイズ磁気だけを同時に測定することにより得られたノイズ磁気検出信号を利用し、生体磁気測定用のＳＱＵＩＤセンサの磁気検出信号に含まれているノイズ成分を除去する補正処理を行うようなことも提案されてはいる。
【０００７】
この場合、ノイズ磁気検出専用の磁気センサと生体磁気測定用のＳＱＵＩＤセンサとは設置位置が異なっており、ノイズ磁気検出専用の磁気センサで得るノイズ磁気検出信号は、生体磁気測定用のＳＱＵＩＤセンサの磁気検出信号に含まれているノイズ成分と正確に対応しているわけではないので、ノイズ磁気検出専用の磁気センサで得たノイズ磁気検出信号に基づき、生体磁気測定用のＳＱＵＩＤセンサの磁気検出信号に含まれているノイズ成分を推定することになる。しかし、ノイズ磁気が複雑な様相を呈するものであることなどから、空間的に異なる位置のノイズ成分を正確に求めることは非常に難しく、その結果、生体磁気測定用のＳＱＵＩＤセンサの磁気検出信号からノイズ成分を十分に除去することは、やはり望めない。
【０００８】
そこで、本出願人は、磁気検出信号のノイズ成分に対処すべく検討を重ねた結果、いわゆるＩＣＡ（独立成分分析）手法に従って、複数個のセンサにより検出された原検出信号（観測信号）を、生体活動の各電流源とその他の各信号源ごとの独立成分に分解した上で、各独立成分について各独立成分自体の状態のみに基づきノイズ成分であるか否かを判定し、ノイズ成分と判定された独立成分は除いた後、残りの非ノイズ独立成分によって検出信号を各々の情報源毎に分離して復元する構成を備えた点を特徴とする生体信号計測装置を、先に特願平１１−１７３８３９号として提案している。この出願人の先願に係る生体信号計測装置は、ノイズ検出専用のセンサを別途に設けずとも、ノイズ成分が十分に除去された生体信号解析用の復元検出信号に基づいて正確な信号解析が可能である。
【０００９】
しかしながら、先願の生体信号計測装置でも、原検出信号の数（即ちセンサの個数）が、各生体活動電流源（信号源）と、その他のノイズ源などを含む各発生源（信号源）とを合計した発生源全個数より大幅に多い場合や、ノイズレベルが大きい（Ｓ／Ｎ比が低い）場合は、独立成分への分解が適正に行われないという問題がある。センサの個数が発生源全個数より大幅に多い場合は、分解後は一つの独立成分だけに纏まっているはずの成分が複数の独立成分に分割されたり、またノイズレベルが大きい場合は、分解後は単独の成分のみからなるはずの独立成分が複数の独立成分を含んでいて、独立成分への分解が適正には行われず、信号解析の信頼度が悪くなるという不都合がある。
【００１０】
この発明は、上記の事情に鑑み、センサの個数が発生源全個数より多い場合や、ノイズレベルが大きい場合であっても、独立成分への分解が適正に行われて正確な信号解析結果を得ることができる生体信号計測装置を提供することを課題とする。
【００１１】
【課題を解決するための手段】
前記課題を解決するために、この発明に係る生体信号計測装置は、被検体の診断対象領域において生体活動電流源により生じる微小な生体信号を測定する複数個のセンサを備えた生体信号計測装置において、前記複数個のセンサにより検出された検出信号における独立成分の個数を求出する独立成分個数求出手段と、検出信号に乗算処理を行って、検出信号を（ａ）複数個のセンサにより検出された検出信号において生体活動電流源による生体信号成分と（ｂ）ノイズ成分との和とすることにより、生体信号成分とノイズ成分との相互依存性を解消する（生体信号成分とノイズ成分との間を無相関化する）成分間相互依存解消手段と、生体信号成分とノイズ成分との相互依存性が解消された検出信号を前記独立成分個数求出手段により求出した個数の数の独立成分に分解する信号分解手段と、この信号分解手段により分解された独立成分の中のノイズ成分を独立成分自体の状態のみに基づいて判定して取り除くノイズ成分除去手段と、各非ノイズ独立成分に基づき検出信号を復元する信号復元手段と、復元された各検出信号の各々の独立成分に対応する脳内の活動（位置、方向、強さ）を得る信号解析手段とを備えている。
【００１２】
〔作用〕
次に、この発明の生体信号計測装置により生体信号の計測を行う時のノイズ成分の除去作用について説明する。
この発明の装置により生体信号計測を実行する際は、先ず複数個のセンサが被検体の診断対象領域の直ぐ傍にセットされて、各センサによって例えば生体活動電流源により発生する微小な生体磁気が複数個のセンサにより各々測定されて原検出信号として出力される。
【００１３】
そして、原検出信号における独立成分の個数が独立成分個数求出手段により求出されるとともに、成分間相互依存解消手段により、原検出信号において生体活動電流源による生体信号成分とノイズ成分との相互依存性が解消される（つまり生体信号成分とノイズ成分との間が無相関化される）。ついで、信号分解手段により、生体信号成分とノイズ成分との相互依存性が解消された検出信号は、独立成分個数求出手段により予め求出した個数の数の独立成分に分解される。その後、ノイズ成分除去手段により、各独立成分の中のノイズ成分が各独立成分自体の状態のみに基づいて判定されて取り除かれてから、信号復元手段によってノイズ成分が除去された残りの非ノイズ独立成分に基づき検出信号（復元検出信号）が復元されて信号解析手段へ送出される。そして、信号解析手段によって、復元検出信号に基づいて生体解析が行われ、個々の生体信号の場所と活動波形が的確に把握される。
【００１４】
この発明の装置では、出願人の先願に係る発明の装置と同様、統計的に独立性の高い複数の信号に分解する独立成分分析（ＩＣＡ：Independence Component Analysis)手法に従って、センサにより検出された原検出信号（観測信号）が生体活動の各電流源とその他の各発生源ごとの独立成分に分解された上で、各独立成分について各独立成分自体の状態のみに基づきノイズ成分であるか否かが判定される。ノイズ成分と判定された独立成分が除かれた後、残りの各非ノイズ独立成分によって検出信号が復元される。その結果、ノイズ検出専用のセンサを別途に設けずとも、生体信号測定用のセンサで得られた信号だけで、ノイズ成分が十分に除去された生体信号解析用の復元検出信号が各成分に分解した形で得られる。
【００１５】
そして、この発明の装置の場合、加えて、独立成分に分解される前に、原検出信号における生体活動電流源による生体信号成分とノイズ成分との間が無相関化されるので、ハイレベルの（大きい）ノイズ成分であっても、独立成分の分解の際にノイズ成分が生体信号成分の側に紛れ込み難くなっていて、分解後は単独の成分のみからなるはずの独立成分が複数の独立成分を含む未独立状態となる心配がない。また、独立成分の分解は予め求出しておいた原検出信号の独立成分の個数の独立成分に検出信号が分解されるので、センサの個数が発生源全個数より大幅に多い場合であっても、検出信号の分解個数が事前に十分に絞られていて、分解後は一つの独立成分だけに纏まっているはずの成分が複数の独立成分に分割されてしまう心配もない。
すなわち、この発明では、センサの個数が発生源全個数より大幅に多い場合や、ノイズレベルが大きい場合であっても、検出信号の独立成分への分解は適正に行われるのである。
【００１６】
【発明の実施の形態】
続いて、この発明の一実施例を図面を参照しながら説明する。図１は本発明に係る生体信号計測装置の一例である生体磁気計測装置の全体構成を示すブロック図である。
【００１７】
実施例の生体磁気計測装置は、図１に示すように、被検体（患者）Ｍの診断対象領域において生体活動電流源により生じる微小な生体磁気を測定するマルチチャンネルＳＱＵＩＤセンサ１と、マルチチャンネルＳＱＵＩＤセンサ１で得られる出力データ（出力信号）を適宜に変換・収集して原磁気検出信号として出力するデータ変換収集部２と、データ変換収集部２から送られてくる原磁気検出信号のノイズ成分除去処理や生体磁気解析処理などを行う信号処理部３とを備えている。
【００１８】
マルチチャンネルＳＱＵＩＤセンサ１は、デュアーと呼ばれる容器内に微小な生体磁気の測定に適した多数のＳＱＵＩＤセンサ１ａが縦横にアレイ状に配列された形で液体窒素などの冷媒に浸漬・収納されている構成になっている。実施例装置のマルチチャンネルＳＱＵＩＤセンサ１には、１２８チャンネル分のＳＱＵＩＤセンサ１ａが配備されており、各チャンネルの１回当たりの生体磁気の測定時間は例えば５１２ｍｓｅｃ（ミリ秒）である。なお、多数のＳＱＵＩＤセンサ１ａがこの発明における複数個のセンサに相当する。
【００１９】
また、マルチチャンネルＳＱＵＩＤセンサ１の後段のデータ変換収集部２は、各ＳＱＵＩＤセンサ１ａの出力信号をディジタル信号に変換して収集し、これを原磁気検出信号として信号処理部３へ送り込む。実施例装置の場合、原磁気検出信号は１２８行５１２列の行列形態の信号として扱われる構成となっている。つまり、各ＳＱＵＩＤセンサ１ａの出力信号を１ｋＨのサンプリング周期（１ミリ秒間隔）で５１２回続けて取り込み、データ数５１２個の原磁気検出信号を、１２８チャンネル分得ており、原磁気検出信号は１２８行５１２列の行列の形態で扱われることになるのである。したがって、原磁気検出信号の行列における行数や列数は、ＳＱＵＩＤセンサ１ａのチャンネル数や出力信号のサンプリング回数に応じて変わることになる。
【００２０】
信号処理部３は、この発明の生体磁気計測装置における特徴的な構成部分であって、原磁気検出信号に基づいて原磁気検出信号における独立成分の個数ｍ0 、および後で信号処理に用いる後述の二つの変換行列Ｗp ，Ｗq を予め先に求出する予備処理部４と、原磁気検出信号において生体活動電流源による生体信号成分とノイズ成分との相互依存性を解消する成分間相互依存解消部５と、生体信号成分とノイズ成分との相互依存性が解消された磁気検出信号を複数の独立成分に分解する磁気信号分解部６と、磁気信号分解部６により分解された独立成分の中のノイズ成分を独立成分自体の状態のみに基づいて判定して取り除くノイズ成分除去部７と、ノイズ成分が除去された残りの非ノイズ独立成分に基づき元の磁気検出信号の形へ戻して各成分毎の復元磁気検出信号として個別に出力する磁気信号復元部８とを備えているとともに、磁気信号復元部８からの復元磁気検出信号に基づいて生体磁気解析を行う磁気解析部９を備えている他、生体内の分極（ダイポール）を破壊して生体活動電流を流すための刺激を被検体Ｍ（生体）に与えるための刺激付与部ＤＴなどを備えている。以下、この信号処理部３の各部構成を、より詳しく説明する。
【００２１】
すなわち、予備処理部４は、計測実行に伴って得られる１２８行５１２列の行列形態の原磁気検出信号Ｂにおける独立成分の個数ｍ0 を求出する演算を行うのに加え、成分間相互依存解消部５で用いるFactor Loading Matrix の疑似逆行列としてのｍ0 行１２８列の第１変換行列Ｗp 、および、磁気信号分解部６で用いるRotation Matrix（回転行列）としてのｍ0 行１２８列の第２変換行列Ｗq とを求める演算を行うよう構成されている。独立成分の個数ｍ0 や各変換行列Ｗp ，Ｗq を求める演算プロセスは後で具体的に説明する。
【００２２】
一方、成分間相互依存解消部５は、生体信号成分とノイズ成分との相互依存性を解消する処理として、第１変換行列Ｗp と原磁気検出信号Ｂとの積（Ｗp Ｂ）を求める行列演算を行うよう構成されている。この演算は原磁気検出信号Ｂに対してＦＡ手法によるSphering処理を施すことに相当する。（Ｗp Ｂ）はｍ0 行５１２列の行列となる。
他方、磁気信号分解部６は、成分間相互依存解消後の磁気検出信号である（Ｗp Ｂ）をｍ0 個の独立成分に分解する処理として、第２変換行列Ｗq と（Ｗp Ｂ）との積〔Ｗq （Ｗp Ｂ）〕を求める行列演算を行うよう構成されている。この演算はＩＣＡ手法によるRotation処理を施すことに相当している。磁気信号分解部６による行列演算の結果、〔Ｗq （Ｗp Ｂ）〕はｍ0 行５１２列のＩＣＡ（独立成分分析）信号Ｘとなる。なお、原磁気検出信号Ｂの独立成分の個数ｍ0 は、発生源全個数であって、また変換行列Ｗp ，Ｗq の行数ともなっている。
【００２３】
続いて、第１変換行列Ｗp および独立成分の個数ｍ0 を求める演算プロセスを説明する。実施例装置の予備処理部４の場合、第１変換行列Ｗp と個数ｍ0 は、ＦＡ（Factor Analysis ＝因子分析）手法に従って求められる。ＦＡ手法の場合、先ず次の（１）式のように、原磁気検出信号Ｂと、Factor Loading Matrix （因子負荷量行列）Ａ（＝Ｗp -1）とが、原磁気検出信号Ｂの行列と同一行列構成の因子ベクトルｆおよび外乱ε（独自因子分散行列）を介して関係づけられるよう設定される。因子ベクトルｆは生体信号成分に相当し、εはノイズ成分に相当する。
Ｂ＝Ａｆ＋ε ・・・（１）
因子負荷量行列Ａは、センサの数と同じ行数（１２８行）と因子の数と同じ列数（ｍ列）の行列である。そして、ＦＡ手法では統計処理で因子負荷量行列Ａを求出することになるが、例えば次の（２）式で示すＭＬＥ(Maxｉmum Likelihood Estimate)評価関数Ｌ（Ａ，σ）を用いる方法が考えられる。
【００２４】
Ｌ（Ａ，σ）＝−｛tr〔Ｃ（σ＋ＡＡT ）-1〕＋Log det(σ＋ＡＡT ）＋128Log2 π｝・・・（２）
但し、σはεに対応する対角行列，trはトレース（対角要素の和），Ｃは原磁気検出信号Ｂの相関行列の平均値（Σi=1 N ＢＢT ／Ｎ：（但し、Ｎはサンプリング数，ＢT はＢの転置行列），ＡT はＡの転置行列，det は行列式の値である。
そして、ＭＬＥ評価関数Ｌ（Ａ，σ）が最大となる時のＡM ，σM が求める行列として決定される。なお、ＭＬＥ評価関数Ｌ（Ａ，σ）の値を求出するのには、例えばニューラルネットワーク法に属するＥＭアルゴリズムを用いる方法が考えられる。
【００２５】
一方、第１変換行列Ｗp の行数でもある独立成分の個数ｍ0 は、次の（３）式で示すＭＤＬ評価関数が用いられる。
ＭＤＬ＝−Ｌ（ＡM ，σM ）＋（ LogＮ）÷Ｎ×｛128(m+1)-m(m-1)/2 ｝・・・（３）
この（３）式において、ｍ＝１，ｍ＝２，・・・, ｍ＝ (1/2)｛2 ×128 −√(8×128 ＋1)｝の時の各ＭＤＬ評価関数の値を個々に求めて、ＭＤＬ値が最小値となる時のｍを、独立成分の個数ｍ0 と決定して因子負荷量行列Ａを完成させた後、完成した１２８行ｍ0 列の因子負荷量行列Ａを逆行列へ変換することにより、ｍ0 行１２８列の第１変換行列Ｗp を求出する。
【００２６】
つまり、上の（１）式に示すように、第１変換行列Ｗp の元である因子負荷量行列Ａは、原磁気検出信号Ｂを生体信号成分に相当する因子ベクトルｆと、ノイズ成分に相当する外乱εとに別かれていることを前提として設定されているので、逆に第１変換行列Ｗp で原磁気検出信号Ｂを乗算処理（Ｗp Ｂ）することにより、原磁気検出信号Ｂにおいて生体信号成分とノイズ成分の間の相互依存性を解消（無相関化）することができるのである。
このように、原磁気検出信号が独立成分に分解される前に、原磁気検出信号における生体活動電流源による生体信号成分とノイズ成分との間が無相関化されていれば、大きなノイズ成分であっても独立成分の分解の際に生体信号成分の側にノイズ成分が紛れ込み難くなり、また、分解後は単独の成分のみからなるはずの独立成分が複数の独立成分を含む未独立状態となる心配がなく、Ｓ／Ｎ比が低い場合でも、磁気検出信号の独立成分への分解は適正になされる。
【００２７】
次に、第２変換行列Ｗq を求める演算プロセスを説明する。実施例装置の予備処理部４の場合、第２変換行列Ｗq はＩＣＡ手法に従って求められる。第２変換行列Ｗq は、無相関化された磁気検出信号をｍ0 個の独立成分に分解するものであって、先に求出した第１変換行列Ｗp と原磁気検出信号Ｂ並びに次の（４）式に従って求出される。
Σk=1 r Σi NO j｜（Ｗq Ｍk Ｗq T ）ij｜2 ・・（４）
即ち、この（４）式の値（対角要素の和）を最小にする時のＷq が求めるｍ0 行１２８列の第２変換行列Ｗq として求出される。
但し、Ｗq T はＷq の転置行列，（Ｗq Ｍk Ｗq T ）ijは行列（Ｗq Ｍk Ｗq T ）のｉｊ成分，「ｉ ＮＯ ｊ」はｉ≠ｊ，Ｍk は次の（５）式で示す通りのものである。
Ｍk ＝〈Ｗp Ｂ（ｔ）Ｗp Ｂ（ｔ＋τk ）〉 ・・（５）
但し、τk はデータサンプリングの時間差であって、ｋ＝１，・・・，ｒであり、具体的には，１ｍｓｅｃ，２ｍｓｅｃ，・・・，ｒｍｓｅｃである。また、〈 〉はアンサンブル平均であることを示す。
【００２８】
実施例装置の磁気信号分解部６は、独立成分の分解の際、予め求出しておいた個数ｍ0 の独立成分に原磁気検出信号を分解するので、センサの個数が発生源全個数より大幅に多い場合であっても、検出信号の分解個数が事前に十分に絞られていて、分解後は一つの独立成分だけに纏まって含まれているはずの成分が複数の独立成分に分割されてしまう心配がなく、磁気検出信号の独立成分への分解は適正になされる。
【００２９】
また、実施例装置の場合、磁気信号分解部６は、同一事象についてＩＣＡ信号Ｘを複数回繰り返し求めて加算平均する構成にもなっている。なお、実施例の装置の場合、ＩＣＡ信号Ｘを加算平均する代わりに、予備処理部４で原磁気検出信号Ｂを加算平均することもできるよう構成されている。
ここでの繰り返し回数は、数回〜数百回までの間の適宜の回数が選ばれる。この構成により、例えば、音を聞いた時の脳の反応を検査する場合、ＩＣＡ信号Ｘまたは原磁気検出信号Ｂの加算平均処理によって、目の筋肉から発生する（スパイク波的な）磁気や脳から定常的に発生するα波による磁気の他、量子ノイズなどの不要成分を除去できる。
なお、逆に、目の筋肉から発生する磁気や脳から定常的に発生するα波による磁気を残したい場合には、加算平均の繰り返し回数を少なくするか、或いは加算平均しないようにすればよい。
【００３０】
続いて、ノイズ成分除去部７より後段の構成を具体的に説明する。ノイズ成分除去部７は、磁気信号分解部６によって求められたＩＣＡ信号Ｘにおけるｍ0 個の各独立成分についてノイズ成分であるか否かを先ず判別する。実施例装置の場合、ＩＣＡ信号Ｘの各行ベクトルについて、全測定時間５１２ｍｓｅｃのうち測定開始から被検体に刺激付与部ＤＴにより刺激が与えられる時点までの非検査対象区間（例えば０〜１００ｍｓｅｃ）の標準偏差値Ｍａと、刺激付与部ＤＴにより刺激が与えられた時点以降の検査対象区間（例えば１００ｍｓｅｃ〜５１２ｍｓｅｃ）の標準偏差値Ｍｂの比Ｍａ／Ｍｂを求め、これが一定値以上の場合、その行ベクトルに対応する独立成分はノイズ成分であると判別し、一定値未満の場合、その行ベクトルに対応する独立成分は真の信号成分（生体信号成分）であると判別するように構成されている。
【００３１】
非ノイズ独立成分（生体信号成分）を決定づける生体磁気は、刺激付与部ＤＴにより刺激が与えられた時点以降に発生するので、行ベクトルの要素は刺激が与えられた時点以降に大きくなる。逆に、ノイズ成分に対応するノイズ磁気は、刺激付与部ＤＴにより刺激が与えられた時点以降とは直接関係がなく、行ベクトルの要素は刺激が与えられた時点の前後で変化が少ない。したがって、非検査対象区間の標準偏差値Ｍａと検査対象区間の標準偏差値Ｍｂの比Ｍａ／Ｍｂについては、真の信号成分である独立成分の方の比Ｍａ／Ｍｂは小さく、ノイズ独立成分の方の比Ｍａ／Ｍｂは大きくなり、比Ｍａ／Ｍｂの大小を監視することで独立成分がノイズ成分であるか否かの判別が可能となる。
【００３２】
すなわち、ノイズ成分除去部７は、ノイズ成分（の中でも特に環境ノイズと称するようなノイズ成分）であると判別されたものを除去する。ノイズ成分の除去は、ノイズ成分と判別された独立成分に対応する行ベクトルを０に置換することにより行われる。ノイズ成分除去部７によるノイズ成分除去に伴ってＩＣＡ信号ＸはＩＣＡ信号Ｘａとなる。
【００３３】
磁気信号復元部８は、ノイズ成分除去部７で求められたＩＣＡ信号Ｘａと、ＩＣＡ（独立成分分析）行列としての第２変換行列Ｗq の逆行例である１２８行ｍ0 列の逆ＩＣＡ行列Ｗq -1とを用いて、１２８行５１２列の復元磁気検出信号Ｂａを求める演算を行う。つまり、磁気信号復元部８において、Ｗq -1Ｘａなる行列演算が行われて、原磁気検出信号Ｂからノイズ成分が十分に除去された復元磁気検出信号Ｂａ（＝Ｗq -1Ｘａ）が求められるのである。
【００３４】
なお、磁気信号復元部８においては、独立信号源ごとに磁気検出信号を復元することも可能である。つまり、この場合には、ノイズ成分除去部７で求められたＩＣＡ信号Ｘａの中で、復元独立成分に対応する行ベクトルのみを残し、その他の要素を全て『０』に置換したＩＣＡ信号Ｘａ’を用いてＷq -1Ｘａ’なる演算を各独立成分毎に行うことにより、復元磁気検出信号Ｂａを求めることもできるのである。
【００３５】
そして、磁気解析部９は、復元磁気検出信号Ｂａに基づいて生体磁気解析を行う。具体的には、復元磁気検出信号Ｂａの１ダイポール解析の結果、生体活動電流源の重心位置分布が求められたり、復元磁気検出信号ＢａのSpatial Filterの結果、生体活動電流源の空間分布が求められたりして、生体活動電流源の状態が把握できる。
【００３６】
実施例装置によれば、ＩＣＡ（独立成分分析）手法に従って、原磁気検出信号Ｂが、各磁気発生源ごとの独立成分に分解された上で、各独立成分について各独立成分自体の状態のみに基づきノイズ成分であるか否かが判定され、ノイズ成分と判定された独立成分が除かれた後、非ノイズ独立成分に従って磁気検出信号が復元されており、ノイズ磁気検出専用の磁気センサを別途に設けずとも、生体磁気測定用の磁気センサで得られた信号だけで、ノイズ成分が十分に除去された生体磁気解析用の復元磁気検出信号が容易に得られている。
【００３７】
なお、実施例装置の信号処理部３は、コンピュータおよびその制御プログラム等を中心に構成されているものである。
【００３８】
さらに、実施例装置は、解析結果を画面に映し出す表示モニタ１０および解析結果をシートに印刷して出力するプリンター１１といった出力機器類を備えており、必要に応じて磁気解析部９で得られた生体活動電流源の重心位置や空間分布を表示モニタ１０に表示させたり、プリンター１１で印刷させたりできる構成にもなっている。
【００３９】
続いて、以上に詳述した構成を有する実施例の生体磁気計測装置により、生体磁気の計測を行う時の装置動作を、図面を参照しながら具体的に説明する。図２は実施例装置による生体磁気の計測実行の様子を経時的に示すフローチャートであり、図３はＩＣＡ信号Ｘの一部を示すグラフである。
なお、以下の測定の場合、図３に示すように、測定開始から１００ｍｓｅｃ経過した時点ＴＭで被検体Ｍに「絵」を見せて被検体Ｍの後頭部側大脳内下方寄りに生体活動電流源が集中して生じるような刺激を付与した。また、計測現場は、ノイズレベルが相当に高くてＳ／Ｎ比は低い状況にある。
【００４０】
〔ステップＳ１〕被検体Ｍの傍らに１２８チャンネルＳＱＵＩＤセンサ１をセットして、測定を所定回数繰り返し実行することにより、１２８行５１２列の原磁気検出信号Ｂが信号処理部３へ所定回数繰り返し入力される。
【００４１】
〔ステップＳ２〕予備処理部４により、原磁気検出信号Ｂが加算平均されるとともに、加算平均された原磁気検出信号Ｂに基づき原磁気検出信号Ｂの独立成分の個数ｍ0 と第１，第２の両変換行列Ｗp ，Ｗq が求出される。ここでは独立成分の個数ｍ0 は１９であり、変換行列Ｗp ，Ｗq は１９行１２８列となった。
原磁気検出信号Ｂの加算平均処理によって、不規則的に発生するノイズ成分や量子ノイズなどの不要成分除去が可能となる。
【００４２】
〔ステップＳ３〕成分間相互依存解消部５により、Ｗp Ｂなる行列演算が行われることで、原磁気検出信号Ｂにおいて生体信号成分とノイズ成分との間が無相関化される。
【００４３】
〔ステップＳ４〕磁気信号分解部６により、Ｗq （Ｗp Ｂ）なる行列演算が行われることにより、磁気検出信号が１９個の独立成分に分解されてＩＣＡ信号Ｘが求められる。図３には１９個の独立成分のうち９個分の独立成分ＩＣＡ１〜ＩＣＡ９を図示する。
【００４４】
〔ステップＳ５〕ノイズ成分除去部７によって、ＩＣＡ信号Ｘにおける各独立成分がノイズ成分であるか否かが判定されて、ノイズ成分と判定された独立成分は除かれ、ＩＣＡ信号Ｘａが求められる。例えば、図３の場合、独立成分ＩＣＡ８，ＩＣＡ９は電源などによるノイズ成分として除かれる。
【００４５】
〔ステップＳ６〕磁気信号復元部８によってＷq -1Ｘａなる行列演算が行われて復元磁気検出信号Ｂａが求められる。
【００４６】
〔ステップＳ７〕得られた復元磁気検出信号Ｂａに従って磁気解析部９によりSpatial Filter法に準拠した生体磁気解析処理が行われ、被検体Ｍの大脳内における生体活動電流源の空間分布が求められるとともに、表示モニタ１０やプリンター１１により解析結果が出力されて、計測は終了となる。
【００４７】
図４は解析処理で得られた結果を示す図であって、被検体Ｍの後頭部側大脳内における生体活動電流源の発生状況を示す模式図である。図４の場合、各矢印がが各生体活動電流源に相当し、矢印の長さは電流源の強度に比例し、矢印の方向が電流源の向きを示す。長い矢印が集中している箇所に真の電流源が存在していると推定される。
また、図５は、図４の中に一点鎖線で示す９個の各ポイント（１０，１０）〜（１２，１２）に対応する解析電流源の時間変化を示すグラフである。例えば，図４と図５は、図４のポイント（１２，１２）の解析電流源は図５の（１２，１２）の直ぐ右のグラフが示すという対応付けになっている。図４の各グラフは横軸が時間，縦軸が電流源強度を示し、略中央の縦直線が刺激付与後１００ｍｓｅｃ経過した時点を示す。
【００４８】
なお、図６は従来装置により得られた解析結果を示す図であって、図４と同様、被検体Ｍの後頭部側大脳内における生体活動電流源の発生状況を示す模式図である。また、図７は、図５と同様の９個のポイントに対応する従来装置の各解析電流源の時間変化を示すグラフである。
【００４９】
実施例装置による解析結果は、図４に示すように、生体活動電流源が後頭部側大脳内下方寄りによく集中して発生しており、被検体Ｍに加えられた刺激によく対応した正確な解析結果が得られている。また、復元磁気検出信号も、正確な結果に呼応して、図５に示すように、刺激付与から検出信号強度が最大となる予定時点である１００ｍｓｅｃ経過後の位置にピークが単一で出現するという適正な波形となっている。
【００５０】
これに対して、従来装置の場合、図６に示すように、生体活動電流源が後頭部側大脳内略中央近傍に二つに別れる感じで非集中的に発生しており、被検体Ｍに加えられた刺激と対応しておらず、解析結果は不正確である。
また、従来装置の場合、解析電流源の時間変化も、図７に示すように、刺激付与から検出信号強度が最大となる予定時点である１００ｍｓｅｃ経過後の位置以外のところにも、顕著なピークが出現するという不適正な波形となっている。従来装置の場合、センサの個数が１２８個と電流源１９個（独立成分の個数に相当）より大幅に多く、またＳ／Ｎ比が低いために、独立成分への分解が適正に行われず、正確な解析結果が得られないと推察される。
【００５１】
しかし、図５および図６の結果から分かるように、実施例装置によれば、センサの個数が１２８個と電流源１９個より大幅に多くても、またＳ／Ｎ比が低くても、磁気検出信号において生体信号成分とノイズ成分の間を無相関化するとともに、予め求めた独立成分の個数ｍ0 の独立成分に磁気検出信号を分解することにより、磁気検出信号の独立成分への分解が適正に行われるので、正確な解析結果を得ることができるのである。
【００５２】
この発明は、上記実施の形態に限られることはなく、下記のように変形実施することができる。
【００５３】
（１）実施例装置では、ＩＣＡ信号Ｘの独立成分における非検査対象区間の標準偏差値Ｍａと検査対象区間の標準偏差値Ｍｂの比Ｍａ／Ｍｂに基づき独立成分がノイズ成分であるか否かを判別する構成であったが、ＩＣＡ信号Ｘの独立成分の信号強度に基づき独立成分がノイズ成分であるか否かを判別する構成の装置が、変形例として挙げられる。例えば１ｐＴ（ピコテスラ）以上のものはノイズと判別する。但し、この変形例の場合、判定対象の独立成分の信号強度は正規化された値ではなく、生の磁場強度に対応する値に変換される必要がある。
【００５４】
（２）また、ＩＣＡ信号Ｘの独立成分の周波数に基づき独立成分がノイズ成分であるか否かを判別する構成の装置も、変形例として挙げられる。例えば、ＩＣＡ信号ＸにＦＦＴをかけて一番支配的な周波数を求め、例えば１００Ｈｚ以上であればノイズと判別する。
【００５５】
（３）さらに、実施例を含めて前述した三つのノイズ成分判別方式の二つあるいは三つの方式を併用して独立成分がノイズ成分であるか否かを判別する構成の装置も、変形例として挙げられる。
【００５６】
また、上述した実施例では、脳磁計による生体磁気計測装置を例に採って説明したが、脳波計などを用いた生体信号計測装置であっても同様の効果を得ることができる。
【００５７】
【発明の効果】
以上に詳述したように、この発明の生体信号計測装置によれば、所謂ＩＣＡ（独立成分分析）手法に従って、複数個のセンサにより検出された原検出信号（観測信号）は、生体活動の各電流源とその他の各信号源ごとの独立成分に分解された上で、各独立成分について各独立成分自体の状態のみに基づきノイズ成分であるか否かが判定され、ノイズ成分と判定された独立成分が除かれた後、残りの非ノイズ独立成分によって検出信号が各々の情報源毎に分離して復元される。したがって、ノイズ検出専用のセンサを別途に設けずとも、生体信号測定用のセンサで得られた信号だけで、ノイズ成分が十分に除去された生体信号解析用の復元検出信号が容易に得られる。その結果、復元生体検出信号に基づいて行われる信号解析も正確なものとなる。
【００５８】
加えて、この発明の生体信号計測装置によれば、独立成分に分解される前に、原検出信号における生体活動電流源による生体信号成分とノイズ成分との間が無相関化される構成を備えているので、大きなノイズ成分であっても、独立成分の分解の際にノイズ成分が生体信号成分の側に紛れ込み難くなる。また、分解後は単独の成分のみからなるはずの独立成分が複数の独立成分を含む未独立状態となる心配がなくなる。さらに、予め求出しておいた独立成分の個数に等しい個数の独立成分に原検出信号が分解されるので、センサの個数が発生源全個数より大幅に多い場合であっても、検出信号の分解個数が事前に十分に絞られていて、分解後は一つの独立成分だけに纏まって含まれているはずの成分が複数の独立成分に分割されてしまう心配がなくなる。したがって、センサの個数が発生源全個数より大幅に多い場合や、ノイズレベルが大きい場合でも、検出信号の独立成分への分解は適正になされることから、正確な信号解析が行われることとなる。
【図面の簡単な説明】
【図１】実施例に係る生体磁気計測装置の全体構成を示すブロック図である。
【図２】実施例装置による生体磁気の計測実行の様子を経時的に示すフローチャートである。
【図３】実施例装置により得られたＩＣＡ信号Ｘの一部を示すグラフである。
【図４】実施例装置により解析された後頭部側大脳内の生体活動電流源の発生状況を示す模式図である。
【図５】実施例装置により計測された後頭部側大脳内の特定の９ポイントに対応する解析電流源の時間変化を示すグラフである。
【図６】従来装置により解析された後頭部側大脳内の生体活動電流源の発生状況を示す模式図である。
【図７】従来装置により計測された後頭部側大脳内の特定の９ポイントに対応する解析電流源の時間変化を示すグラフである。
【符号の説明】
１ …マルチチャンネルＳＱＵＩＤセンサ
１ａ …ＳＱＵＩＤセンサ
２ …データ変換収集部
３ …信号処理部
４ …予備処理部
５ …成分間相互依存解消部
６ …磁気信号分解部
７ …ノイズ成分除去部
８ …磁気信号復元部
９ …磁気解析部
Ｍ …被検体
【発明の属する技術分野】
この発明は、脳磁計や脳波計などの生体信号計測装置に係り、特に生体信号測定用のセンサにより得られた検出信号からノイズ成分を除去し、脳内活動を推定するための技術に関する。
【０００２】
【従来の技術】
生体内に流れる生体活動電流により微小な生体磁気（生体磁界）が生体から発生する。例えば、脳から発生する生体磁気は脳磁と呼ばれ、生体に刺激を与えることにより発生する誘発脳磁や、α波やてんかんのスパイク波のように自然に発生する自発脳磁などがある。
【０００３】
近年、生体から出る微小な生体磁気を測定できる磁束計として、ＳＱＵＩＤ（Superconducing Quantum Interference Device：超電導量子干渉計）を用いたマルチチャンネルＳＱＵＩＤセンサが開発されている。マルチチャンネルＳＱＵＩＤセンサは、デュアーと呼ばれる容器内に多数個のＳＱＵＩＤセンサを液体窒素などの冷媒に浸漬・収納した構成となっている。
【０００４】
このマルチチャンネルＳＱＵＩＤセンサ（以下、適宜「磁束計」と略記）を備えた生体信号計測装置、つまり生体磁気計測装置の場合、磁束計を被検体の関心部位である例えば頭部の傍らに置くと、頭部内に生じた生体活動電流源から発生する微小な生体磁気が磁束計内の各ＳＱＵＩＤセンサで無侵襲で測定されて磁気検出信号として出力されるとともに、ＳＱＵＩＤセンサからの磁気検出信号に基づいて生体磁気解析が行われて、生体活動電流源の状態（例えば、位置や向き或いは大きさ等）を把握することができるという構成になっている（例えば特開平７−３２７９４３号公報参照）。
【０００５】
【発明が解決しようとする課題】
しかしながら、従来の生体磁気計測装置には、各ＳＱＵＩＤセンサによって検出された磁気検出信号に含まれるノイズ成分を十分に除去することが難しいという問題がある。測定対象である生体磁気は非常に微弱であることから生体活動電流源以外の別の磁気発生源から出る（環境ノイズとも称するような）ノイズ磁気の混入が避けられない。したがって、各ＳＱＵＩＤセンサからの磁気検出信号には混入するノイズ磁気によるノイズ成分が含まれており、磁気検出信号からノイズ成分を十分に除去しなければ、生体磁気の正確な解析はおぼつかないことになる。
【０００６】
そこで、生体磁気測定用のＳＱＵＩＤセンサとは別のノイズ磁気検出専用の磁気センサでノイズ磁気だけを同時に測定することにより得られたノイズ磁気検出信号を利用し、生体磁気測定用のＳＱＵＩＤセンサの磁気検出信号に含まれているノイズ成分を除去する補正処理を行うようなことも提案されてはいる。
【０００７】
この場合、ノイズ磁気検出専用の磁気センサと生体磁気測定用のＳＱＵＩＤセンサとは設置位置が異なっており、ノイズ磁気検出専用の磁気センサで得るノイズ磁気検出信号は、生体磁気測定用のＳＱＵＩＤセンサの磁気検出信号に含まれているノイズ成分と正確に対応しているわけではないので、ノイズ磁気検出専用の磁気センサで得たノイズ磁気検出信号に基づき、生体磁気測定用のＳＱＵＩＤセンサの磁気検出信号に含まれているノイズ成分を推定することになる。しかし、ノイズ磁気が複雑な様相を呈するものであることなどから、空間的に異なる位置のノイズ成分を正確に求めることは非常に難しく、その結果、生体磁気測定用のＳＱＵＩＤセンサの磁気検出信号からノイズ成分を十分に除去することは、やはり望めない。
【０００８】
そこで、本出願人は、磁気検出信号のノイズ成分に対処すべく検討を重ねた結果、いわゆるＩＣＡ（独立成分分析）手法に従って、複数個のセンサにより検出された原検出信号（観測信号）を、生体活動の各電流源とその他の各信号源ごとの独立成分に分解した上で、各独立成分について各独立成分自体の状態のみに基づきノイズ成分であるか否かを判定し、ノイズ成分と判定された独立成分は除いた後、残りの非ノイズ独立成分によって検出信号を各々の情報源毎に分離して復元する構成を備えた点を特徴とする生体信号計測装置を、先に特願平１１−１７３８３９号として提案している。この出願人の先願に係る生体信号計測装置は、ノイズ検出専用のセンサを別途に設けずとも、ノイズ成分が十分に除去された生体信号解析用の復元検出信号に基づいて正確な信号解析が可能である。
【０００９】
しかしながら、先願の生体信号計測装置でも、原検出信号の数（即ちセンサの個数）が、各生体活動電流源（信号源）と、その他のノイズ源などを含む各発生源（信号源）とを合計した発生源全個数より大幅に多い場合や、ノイズレベルが大きい（Ｓ／Ｎ比が低い）場合は、独立成分への分解が適正に行われないという問題がある。センサの個数が発生源全個数より大幅に多い場合は、分解後は一つの独立成分だけに纏まっているはずの成分が複数の独立成分に分割されたり、またノイズレベルが大きい場合は、分解後は単独の成分のみからなるはずの独立成分が複数の独立成分を含んでいて、独立成分への分解が適正には行われず、信号解析の信頼度が悪くなるという不都合がある。
【００１０】
この発明は、上記の事情に鑑み、センサの個数が発生源全個数より多い場合や、ノイズレベルが大きい場合であっても、独立成分への分解が適正に行われて正確な信号解析結果を得ることができる生体信号計測装置を提供することを課題とする。
【００１１】
【課題を解決するための手段】
前記課題を解決するために、この発明に係る生体信号計測装置は、被検体の診断対象領域において生体活動電流源により生じる微小な生体信号を測定する複数個のセンサを備えた生体信号計測装置において、前記複数個のセンサにより検出された検出信号における独立成分の個数を求出する独立成分個数求出手段と、検出信号に乗算処理を行って、検出信号を（ａ）複数個のセンサにより検出された検出信号において生体活動電流源による生体信号成分と（ｂ）ノイズ成分との和とすることにより、生体信号成分とノイズ成分との相互依存性を解消する（生体信号成分とノイズ成分との間を無相関化する）成分間相互依存解消手段と、生体信号成分とノイズ成分との相互依存性が解消された検出信号を前記独立成分個数求出手段により求出した個数の数の独立成分に分解する信号分解手段と、この信号分解手段により分解された独立成分の中のノイズ成分を独立成分自体の状態のみに基づいて判定して取り除くノイズ成分除去手段と、各非ノイズ独立成分に基づき検出信号を復元する信号復元手段と、復元された各検出信号の各々の独立成分に対応する脳内の活動（位置、方向、強さ）を得る信号解析手段とを備えている。
【００１２】
〔作用〕
次に、この発明の生体信号計測装置により生体信号の計測を行う時のノイズ成分の除去作用について説明する。
この発明の装置により生体信号計測を実行する際は、先ず複数個のセンサが被検体の診断対象領域の直ぐ傍にセットされて、各センサによって例えば生体活動電流源により発生する微小な生体磁気が複数個のセンサにより各々測定されて原検出信号として出力される。
【００１３】
そして、原検出信号における独立成分の個数が独立成分個数求出手段により求出されるとともに、成分間相互依存解消手段により、原検出信号において生体活動電流源による生体信号成分とノイズ成分との相互依存性が解消される（つまり生体信号成分とノイズ成分との間が無相関化される）。ついで、信号分解手段により、生体信号成分とノイズ成分との相互依存性が解消された検出信号は、独立成分個数求出手段により予め求出した個数の数の独立成分に分解される。その後、ノイズ成分除去手段により、各独立成分の中のノイズ成分が各独立成分自体の状態のみに基づいて判定されて取り除かれてから、信号復元手段によってノイズ成分が除去された残りの非ノイズ独立成分に基づき検出信号（復元検出信号）が復元されて信号解析手段へ送出される。そして、信号解析手段によって、復元検出信号に基づいて生体解析が行われ、個々の生体信号の場所と活動波形が的確に把握される。
【００１４】
この発明の装置では、出願人の先願に係る発明の装置と同様、統計的に独立性の高い複数の信号に分解する独立成分分析（ＩＣＡ：Independence Component Analysis)手法に従って、センサにより検出された原検出信号（観測信号）が生体活動の各電流源とその他の各発生源ごとの独立成分に分解された上で、各独立成分について各独立成分自体の状態のみに基づきノイズ成分であるか否かが判定される。ノイズ成分と判定された独立成分が除かれた後、残りの各非ノイズ独立成分によって検出信号が復元される。その結果、ノイズ検出専用のセンサを別途に設けずとも、生体信号測定用のセンサで得られた信号だけで、ノイズ成分が十分に除去された生体信号解析用の復元検出信号が各成分に分解した形で得られる。
【００１５】
そして、この発明の装置の場合、加えて、独立成分に分解される前に、原検出信号における生体活動電流源による生体信号成分とノイズ成分との間が無相関化されるので、ハイレベルの（大きい）ノイズ成分であっても、独立成分の分解の際にノイズ成分が生体信号成分の側に紛れ込み難くなっていて、分解後は単独の成分のみからなるはずの独立成分が複数の独立成分を含む未独立状態となる心配がない。また、独立成分の分解は予め求出しておいた原検出信号の独立成分の個数の独立成分に検出信号が分解されるので、センサの個数が発生源全個数より大幅に多い場合であっても、検出信号の分解個数が事前に十分に絞られていて、分解後は一つの独立成分だけに纏まっているはずの成分が複数の独立成分に分割されてしまう心配もない。
すなわち、この発明では、センサの個数が発生源全個数より大幅に多い場合や、ノイズレベルが大きい場合であっても、検出信号の独立成分への分解は適正に行われるのである。
【００１６】
【発明の実施の形態】
続いて、この発明の一実施例を図面を参照しながら説明する。図１は本発明に係る生体信号計測装置の一例である生体磁気計測装置の全体構成を示すブロック図である。
【００１７】
実施例の生体磁気計測装置は、図１に示すように、被検体（患者）Ｍの診断対象領域において生体活動電流源により生じる微小な生体磁気を測定するマルチチャンネルＳＱＵＩＤセンサ１と、マルチチャンネルＳＱＵＩＤセンサ１で得られる出力データ（出力信号）を適宜に変換・収集して原磁気検出信号として出力するデータ変換収集部２と、データ変換収集部２から送られてくる原磁気検出信号のノイズ成分除去処理や生体磁気解析処理などを行う信号処理部３とを備えている。
【００１８】
マルチチャンネルＳＱＵＩＤセンサ１は、デュアーと呼ばれる容器内に微小な生体磁気の測定に適した多数のＳＱＵＩＤセンサ１ａが縦横にアレイ状に配列された形で液体窒素などの冷媒に浸漬・収納されている構成になっている。実施例装置のマルチチャンネルＳＱＵＩＤセンサ１には、１２８チャンネル分のＳＱＵＩＤセンサ１ａが配備されており、各チャンネルの１回当たりの生体磁気の測定時間は例えば５１２ｍｓｅｃ（ミリ秒）である。なお、多数のＳＱＵＩＤセンサ１ａがこの発明における複数個のセンサに相当する。
【００１９】
また、マルチチャンネルＳＱＵＩＤセンサ１の後段のデータ変換収集部２は、各ＳＱＵＩＤセンサ１ａの出力信号をディジタル信号に変換して収集し、これを原磁気検出信号として信号処理部３へ送り込む。実施例装置の場合、原磁気検出信号は１２８行５１２列の行列形態の信号として扱われる構成となっている。つまり、各ＳＱＵＩＤセンサ１ａの出力信号を１ｋＨのサンプリング周期（１ミリ秒間隔）で５１２回続けて取り込み、データ数５１２個の原磁気検出信号を、１２８チャンネル分得ており、原磁気検出信号は１２８行５１２列の行列の形態で扱われることになるのである。したがって、原磁気検出信号の行列における行数や列数は、ＳＱＵＩＤセンサ１ａのチャンネル数や出力信号のサンプリング回数に応じて変わることになる。
【００２０】
信号処理部３は、この発明の生体磁気計測装置における特徴的な構成部分であって、原磁気検出信号に基づいて原磁気検出信号における独立成分の個数ｍ0 、および後で信号処理に用いる後述の二つの変換行列Ｗp ，Ｗq を予め先に求出する予備処理部４と、原磁気検出信号において生体活動電流源による生体信号成分とノイズ成分との相互依存性を解消する成分間相互依存解消部５と、生体信号成分とノイズ成分との相互依存性が解消された磁気検出信号を複数の独立成分に分解する磁気信号分解部６と、磁気信号分解部６により分解された独立成分の中のノイズ成分を独立成分自体の状態のみに基づいて判定して取り除くノイズ成分除去部７と、ノイズ成分が除去された残りの非ノイズ独立成分に基づき元の磁気検出信号の形へ戻して各成分毎の復元磁気検出信号として個別に出力する磁気信号復元部８とを備えているとともに、磁気信号復元部８からの復元磁気検出信号に基づいて生体磁気解析を行う磁気解析部９を備えている他、生体内の分極（ダイポール）を破壊して生体活動電流を流すための刺激を被検体Ｍ（生体）に与えるための刺激付与部ＤＴなどを備えている。以下、この信号処理部３の各部構成を、より詳しく説明する。
【００２１】
すなわち、予備処理部４は、計測実行に伴って得られる１２８行５１２列の行列形態の原磁気検出信号Ｂにおける独立成分の個数ｍ0 を求出する演算を行うのに加え、成分間相互依存解消部５で用いるFactor Loading Matrix の疑似逆行列としてのｍ0 行１２８列の第１変換行列Ｗp 、および、磁気信号分解部６で用いるRotation Matrix（回転行列）としてのｍ0 行１２８列の第２変換行列Ｗq とを求める演算を行うよう構成されている。独立成分の個数ｍ0 や各変換行列Ｗp ，Ｗq を求める演算プロセスは後で具体的に説明する。
【００２２】
一方、成分間相互依存解消部５は、生体信号成分とノイズ成分との相互依存性を解消する処理として、第１変換行列Ｗp と原磁気検出信号Ｂとの積（Ｗp Ｂ）を求める行列演算を行うよう構成されている。この演算は原磁気検出信号Ｂに対してＦＡ手法によるSphering処理を施すことに相当する。（Ｗp Ｂ）はｍ0 行５１２列の行列となる。
他方、磁気信号分解部６は、成分間相互依存解消後の磁気検出信号である（Ｗp Ｂ）をｍ0 個の独立成分に分解する処理として、第２変換行列Ｗq と（Ｗp Ｂ）との積〔Ｗq （Ｗp Ｂ）〕を求める行列演算を行うよう構成されている。この演算はＩＣＡ手法によるRotation処理を施すことに相当している。磁気信号分解部６による行列演算の結果、〔Ｗq （Ｗp Ｂ）〕はｍ0 行５１２列のＩＣＡ（独立成分分析）信号Ｘとなる。なお、原磁気検出信号Ｂの独立成分の個数ｍ0 は、発生源全個数であって、また変換行列Ｗp ，Ｗq の行数ともなっている。
【００２３】
続いて、第１変換行列Ｗp および独立成分の個数ｍ0 を求める演算プロセスを説明する。実施例装置の予備処理部４の場合、第１変換行列Ｗp と個数ｍ0 は、ＦＡ（Factor Analysis ＝因子分析）手法に従って求められる。ＦＡ手法の場合、先ず次の（１）式のように、原磁気検出信号Ｂと、Factor Loading Matrix （因子負荷量行列）Ａ（＝Ｗp -1）とが、原磁気検出信号Ｂの行列と同一行列構成の因子ベクトルｆおよび外乱ε（独自因子分散行列）を介して関係づけられるよう設定される。因子ベクトルｆは生体信号成分に相当し、εはノイズ成分に相当する。
Ｂ＝Ａｆ＋ε ・・・（１）
因子負荷量行列Ａは、センサの数と同じ行数（１２８行）と因子の数と同じ列数（ｍ列）の行列である。そして、ＦＡ手法では統計処理で因子負荷量行列Ａを求出することになるが、例えば次の（２）式で示すＭＬＥ(Maxｉmum Likelihood Estimate)評価関数Ｌ（Ａ，σ）を用いる方法が考えられる。
【００２４】
Ｌ（Ａ，σ）＝−｛tr〔Ｃ（σ＋ＡＡT ）-1〕＋Log det(σ＋ＡＡT ）＋128Log2 π｝・・・（２）
但し、σはεに対応する対角行列，trはトレース（対角要素の和），Ｃは原磁気検出信号Ｂの相関行列の平均値（Σi=1 N ＢＢT ／Ｎ：（但し、Ｎはサンプリング数，ＢT はＢの転置行列），ＡT はＡの転置行列，det は行列式の値である。
そして、ＭＬＥ評価関数Ｌ（Ａ，σ）が最大となる時のＡM ，σM が求める行列として決定される。なお、ＭＬＥ評価関数Ｌ（Ａ，σ）の値を求出するのには、例えばニューラルネットワーク法に属するＥＭアルゴリズムを用いる方法が考えられる。
【００２５】
一方、第１変換行列Ｗp の行数でもある独立成分の個数ｍ0 は、次の（３）式で示すＭＤＬ評価関数が用いられる。
ＭＤＬ＝−Ｌ（ＡM ，σM ）＋（ LogＮ）÷Ｎ×｛128(m+1)-m(m-1)/2 ｝・・・（３）
この（３）式において、ｍ＝１，ｍ＝２，・・・, ｍ＝ (1/2)｛2 ×128 −√(8×128 ＋1)｝の時の各ＭＤＬ評価関数の値を個々に求めて、ＭＤＬ値が最小値となる時のｍを、独立成分の個数ｍ0 と決定して因子負荷量行列Ａを完成させた後、完成した１２８行ｍ0 列の因子負荷量行列Ａを逆行列へ変換することにより、ｍ0 行１２８列の第１変換行列Ｗp を求出する。
【００２６】
つまり、上の（１）式に示すように、第１変換行列Ｗp の元である因子負荷量行列Ａは、原磁気検出信号Ｂを生体信号成分に相当する因子ベクトルｆと、ノイズ成分に相当する外乱εとに別かれていることを前提として設定されているので、逆に第１変換行列Ｗp で原磁気検出信号Ｂを乗算処理（Ｗp Ｂ）することにより、原磁気検出信号Ｂにおいて生体信号成分とノイズ成分の間の相互依存性を解消（無相関化）することができるのである。
このように、原磁気検出信号が独立成分に分解される前に、原磁気検出信号における生体活動電流源による生体信号成分とノイズ成分との間が無相関化されていれば、大きなノイズ成分であっても独立成分の分解の際に生体信号成分の側にノイズ成分が紛れ込み難くなり、また、分解後は単独の成分のみからなるはずの独立成分が複数の独立成分を含む未独立状態となる心配がなく、Ｓ／Ｎ比が低い場合でも、磁気検出信号の独立成分への分解は適正になされる。
【００２７】
次に、第２変換行列Ｗq を求める演算プロセスを説明する。実施例装置の予備処理部４の場合、第２変換行列Ｗq はＩＣＡ手法に従って求められる。第２変換行列Ｗq は、無相関化された磁気検出信号をｍ0 個の独立成分に分解するものであって、先に求出した第１変換行列Ｗp と原磁気検出信号Ｂ並びに次の（４）式に従って求出される。
Σk=1 r Σi NO j｜（Ｗq Ｍk Ｗq T ）ij｜2 ・・（４）
即ち、この（４）式の値（対角要素の和）を最小にする時のＷq が求めるｍ0 行１２８列の第２変換行列Ｗq として求出される。
但し、Ｗq T はＷq の転置行列，（Ｗq Ｍk Ｗq T ）ijは行列（Ｗq Ｍk Ｗq T ）のｉｊ成分，「ｉ ＮＯ ｊ」はｉ≠ｊ，Ｍk は次の（５）式で示す通りのものである。
Ｍk ＝〈Ｗp Ｂ（ｔ）Ｗp Ｂ（ｔ＋τk ）〉 ・・（５）
但し、τk はデータサンプリングの時間差であって、ｋ＝１，・・・，ｒであり、具体的には，１ｍｓｅｃ，２ｍｓｅｃ，・・・，ｒｍｓｅｃである。また、〈 〉はアンサンブル平均であることを示す。
【００２８】
実施例装置の磁気信号分解部６は、独立成分の分解の際、予め求出しておいた個数ｍ0 の独立成分に原磁気検出信号を分解するので、センサの個数が発生源全個数より大幅に多い場合であっても、検出信号の分解個数が事前に十分に絞られていて、分解後は一つの独立成分だけに纏まって含まれているはずの成分が複数の独立成分に分割されてしまう心配がなく、磁気検出信号の独立成分への分解は適正になされる。
【００２９】
また、実施例装置の場合、磁気信号分解部６は、同一事象についてＩＣＡ信号Ｘを複数回繰り返し求めて加算平均する構成にもなっている。なお、実施例の装置の場合、ＩＣＡ信号Ｘを加算平均する代わりに、予備処理部４で原磁気検出信号Ｂを加算平均することもできるよう構成されている。
ここでの繰り返し回数は、数回〜数百回までの間の適宜の回数が選ばれる。この構成により、例えば、音を聞いた時の脳の反応を検査する場合、ＩＣＡ信号Ｘまたは原磁気検出信号Ｂの加算平均処理によって、目の筋肉から発生する（スパイク波的な）磁気や脳から定常的に発生するα波による磁気の他、量子ノイズなどの不要成分を除去できる。
なお、逆に、目の筋肉から発生する磁気や脳から定常的に発生するα波による磁気を残したい場合には、加算平均の繰り返し回数を少なくするか、或いは加算平均しないようにすればよい。
【００３０】
続いて、ノイズ成分除去部７より後段の構成を具体的に説明する。ノイズ成分除去部７は、磁気信号分解部６によって求められたＩＣＡ信号Ｘにおけるｍ0 個の各独立成分についてノイズ成分であるか否かを先ず判別する。実施例装置の場合、ＩＣＡ信号Ｘの各行ベクトルについて、全測定時間５１２ｍｓｅｃのうち測定開始から被検体に刺激付与部ＤＴにより刺激が与えられる時点までの非検査対象区間（例えば０〜１００ｍｓｅｃ）の標準偏差値Ｍａと、刺激付与部ＤＴにより刺激が与えられた時点以降の検査対象区間（例えば１００ｍｓｅｃ〜５１２ｍｓｅｃ）の標準偏差値Ｍｂの比Ｍａ／Ｍｂを求め、これが一定値以上の場合、その行ベクトルに対応する独立成分はノイズ成分であると判別し、一定値未満の場合、その行ベクトルに対応する独立成分は真の信号成分（生体信号成分）であると判別するように構成されている。
【００３１】
非ノイズ独立成分（生体信号成分）を決定づける生体磁気は、刺激付与部ＤＴにより刺激が与えられた時点以降に発生するので、行ベクトルの要素は刺激が与えられた時点以降に大きくなる。逆に、ノイズ成分に対応するノイズ磁気は、刺激付与部ＤＴにより刺激が与えられた時点以降とは直接関係がなく、行ベクトルの要素は刺激が与えられた時点の前後で変化が少ない。したがって、非検査対象区間の標準偏差値Ｍａと検査対象区間の標準偏差値Ｍｂの比Ｍａ／Ｍｂについては、真の信号成分である独立成分の方の比Ｍａ／Ｍｂは小さく、ノイズ独立成分の方の比Ｍａ／Ｍｂは大きくなり、比Ｍａ／Ｍｂの大小を監視することで独立成分がノイズ成分であるか否かの判別が可能となる。
【００３２】
すなわち、ノイズ成分除去部７は、ノイズ成分（の中でも特に環境ノイズと称するようなノイズ成分）であると判別されたものを除去する。ノイズ成分の除去は、ノイズ成分と判別された独立成分に対応する行ベクトルを０に置換することにより行われる。ノイズ成分除去部７によるノイズ成分除去に伴ってＩＣＡ信号ＸはＩＣＡ信号Ｘａとなる。
【００３３】
磁気信号復元部８は、ノイズ成分除去部７で求められたＩＣＡ信号Ｘａと、ＩＣＡ（独立成分分析）行列としての第２変換行列Ｗq の逆行例である１２８行ｍ0 列の逆ＩＣＡ行列Ｗq -1とを用いて、１２８行５１２列の復元磁気検出信号Ｂａを求める演算を行う。つまり、磁気信号復元部８において、Ｗq -1Ｘａなる行列演算が行われて、原磁気検出信号Ｂからノイズ成分が十分に除去された復元磁気検出信号Ｂａ（＝Ｗq -1Ｘａ）が求められるのである。
【００３４】
なお、磁気信号復元部８においては、独立信号源ごとに磁気検出信号を復元することも可能である。つまり、この場合には、ノイズ成分除去部７で求められたＩＣＡ信号Ｘａの中で、復元独立成分に対応する行ベクトルのみを残し、その他の要素を全て『０』に置換したＩＣＡ信号Ｘａ’を用いてＷq -1Ｘａ’なる演算を各独立成分毎に行うことにより、復元磁気検出信号Ｂａを求めることもできるのである。
【００３５】
そして、磁気解析部９は、復元磁気検出信号Ｂａに基づいて生体磁気解析を行う。具体的には、復元磁気検出信号Ｂａの１ダイポール解析の結果、生体活動電流源の重心位置分布が求められたり、復元磁気検出信号ＢａのSpatial Filterの結果、生体活動電流源の空間分布が求められたりして、生体活動電流源の状態が把握できる。
【００３６】
実施例装置によれば、ＩＣＡ（独立成分分析）手法に従って、原磁気検出信号Ｂが、各磁気発生源ごとの独立成分に分解された上で、各独立成分について各独立成分自体の状態のみに基づきノイズ成分であるか否かが判定され、ノイズ成分と判定された独立成分が除かれた後、非ノイズ独立成分に従って磁気検出信号が復元されており、ノイズ磁気検出専用の磁気センサを別途に設けずとも、生体磁気測定用の磁気センサで得られた信号だけで、ノイズ成分が十分に除去された生体磁気解析用の復元磁気検出信号が容易に得られている。
【００３７】
なお、実施例装置の信号処理部３は、コンピュータおよびその制御プログラム等を中心に構成されているものである。
【００３８】
さらに、実施例装置は、解析結果を画面に映し出す表示モニタ１０および解析結果をシートに印刷して出力するプリンター１１といった出力機器類を備えており、必要に応じて磁気解析部９で得られた生体活動電流源の重心位置や空間分布を表示モニタ１０に表示させたり、プリンター１１で印刷させたりできる構成にもなっている。
【００３９】
続いて、以上に詳述した構成を有する実施例の生体磁気計測装置により、生体磁気の計測を行う時の装置動作を、図面を参照しながら具体的に説明する。図２は実施例装置による生体磁気の計測実行の様子を経時的に示すフローチャートであり、図３はＩＣＡ信号Ｘの一部を示すグラフである。
なお、以下の測定の場合、図３に示すように、測定開始から１００ｍｓｅｃ経過した時点ＴＭで被検体Ｍに「絵」を見せて被検体Ｍの後頭部側大脳内下方寄りに生体活動電流源が集中して生じるような刺激を付与した。また、計測現場は、ノイズレベルが相当に高くてＳ／Ｎ比は低い状況にある。
【００４０】
〔ステップＳ１〕被検体Ｍの傍らに１２８チャンネルＳＱＵＩＤセンサ１をセットして、測定を所定回数繰り返し実行することにより、１２８行５１２列の原磁気検出信号Ｂが信号処理部３へ所定回数繰り返し入力される。
【００４１】
〔ステップＳ２〕予備処理部４により、原磁気検出信号Ｂが加算平均されるとともに、加算平均された原磁気検出信号Ｂに基づき原磁気検出信号Ｂの独立成分の個数ｍ0 と第１，第２の両変換行列Ｗp ，Ｗq が求出される。ここでは独立成分の個数ｍ0 は１９であり、変換行列Ｗp ，Ｗq は１９行１２８列となった。
原磁気検出信号Ｂの加算平均処理によって、不規則的に発生するノイズ成分や量子ノイズなどの不要成分除去が可能となる。
【００４２】
〔ステップＳ３〕成分間相互依存解消部５により、Ｗp Ｂなる行列演算が行われることで、原磁気検出信号Ｂにおいて生体信号成分とノイズ成分との間が無相関化される。
【００４３】
〔ステップＳ４〕磁気信号分解部６により、Ｗq （Ｗp Ｂ）なる行列演算が行われることにより、磁気検出信号が１９個の独立成分に分解されてＩＣＡ信号Ｘが求められる。図３には１９個の独立成分のうち９個分の独立成分ＩＣＡ１〜ＩＣＡ９を図示する。
【００４４】
〔ステップＳ５〕ノイズ成分除去部７によって、ＩＣＡ信号Ｘにおける各独立成分がノイズ成分であるか否かが判定されて、ノイズ成分と判定された独立成分は除かれ、ＩＣＡ信号Ｘａが求められる。例えば、図３の場合、独立成分ＩＣＡ８，ＩＣＡ９は電源などによるノイズ成分として除かれる。
【００４５】
〔ステップＳ６〕磁気信号復元部８によってＷq -1Ｘａなる行列演算が行われて復元磁気検出信号Ｂａが求められる。
【００４６】
〔ステップＳ７〕得られた復元磁気検出信号Ｂａに従って磁気解析部９によりSpatial Filter法に準拠した生体磁気解析処理が行われ、被検体Ｍの大脳内における生体活動電流源の空間分布が求められるとともに、表示モニタ１０やプリンター１１により解析結果が出力されて、計測は終了となる。
【００４７】
図４は解析処理で得られた結果を示す図であって、被検体Ｍの後頭部側大脳内における生体活動電流源の発生状況を示す模式図である。図４の場合、各矢印がが各生体活動電流源に相当し、矢印の長さは電流源の強度に比例し、矢印の方向が電流源の向きを示す。長い矢印が集中している箇所に真の電流源が存在していると推定される。
また、図５は、図４の中に一点鎖線で示す９個の各ポイント（１０，１０）〜（１２，１２）に対応する解析電流源の時間変化を示すグラフである。例えば，図４と図５は、図４のポイント（１２，１２）の解析電流源は図５の（１２，１２）の直ぐ右のグラフが示すという対応付けになっている。図４の各グラフは横軸が時間，縦軸が電流源強度を示し、略中央の縦直線が刺激付与後１００ｍｓｅｃ経過した時点を示す。
【００４８】
なお、図６は従来装置により得られた解析結果を示す図であって、図４と同様、被検体Ｍの後頭部側大脳内における生体活動電流源の発生状況を示す模式図である。また、図７は、図５と同様の９個のポイントに対応する従来装置の各解析電流源の時間変化を示すグラフである。
【００４９】
実施例装置による解析結果は、図４に示すように、生体活動電流源が後頭部側大脳内下方寄りによく集中して発生しており、被検体Ｍに加えられた刺激によく対応した正確な解析結果が得られている。また、復元磁気検出信号も、正確な結果に呼応して、図５に示すように、刺激付与から検出信号強度が最大となる予定時点である１００ｍｓｅｃ経過後の位置にピークが単一で出現するという適正な波形となっている。
【００５０】
これに対して、従来装置の場合、図６に示すように、生体活動電流源が後頭部側大脳内略中央近傍に二つに別れる感じで非集中的に発生しており、被検体Ｍに加えられた刺激と対応しておらず、解析結果は不正確である。
また、従来装置の場合、解析電流源の時間変化も、図７に示すように、刺激付与から検出信号強度が最大となる予定時点である１００ｍｓｅｃ経過後の位置以外のところにも、顕著なピークが出現するという不適正な波形となっている。従来装置の場合、センサの個数が１２８個と電流源１９個（独立成分の個数に相当）より大幅に多く、またＳ／Ｎ比が低いために、独立成分への分解が適正に行われず、正確な解析結果が得られないと推察される。
【００５１】
しかし、図５および図６の結果から分かるように、実施例装置によれば、センサの個数が１２８個と電流源１９個より大幅に多くても、またＳ／Ｎ比が低くても、磁気検出信号において生体信号成分とノイズ成分の間を無相関化するとともに、予め求めた独立成分の個数ｍ0 の独立成分に磁気検出信号を分解することにより、磁気検出信号の独立成分への分解が適正に行われるので、正確な解析結果を得ることができるのである。
【００５２】
この発明は、上記実施の形態に限られることはなく、下記のように変形実施することができる。
【００５３】
（１）実施例装置では、ＩＣＡ信号Ｘの独立成分における非検査対象区間の標準偏差値Ｍａと検査対象区間の標準偏差値Ｍｂの比Ｍａ／Ｍｂに基づき独立成分がノイズ成分であるか否かを判別する構成であったが、ＩＣＡ信号Ｘの独立成分の信号強度に基づき独立成分がノイズ成分であるか否かを判別する構成の装置が、変形例として挙げられる。例えば１ｐＴ（ピコテスラ）以上のものはノイズと判別する。但し、この変形例の場合、判定対象の独立成分の信号強度は正規化された値ではなく、生の磁場強度に対応する値に変換される必要がある。
【００５４】
（２）また、ＩＣＡ信号Ｘの独立成分の周波数に基づき独立成分がノイズ成分であるか否かを判別する構成の装置も、変形例として挙げられる。例えば、ＩＣＡ信号ＸにＦＦＴをかけて一番支配的な周波数を求め、例えば１００Ｈｚ以上であればノイズと判別する。
【００５５】
（３）さらに、実施例を含めて前述した三つのノイズ成分判別方式の二つあるいは三つの方式を併用して独立成分がノイズ成分であるか否かを判別する構成の装置も、変形例として挙げられる。
【００５６】
また、上述した実施例では、脳磁計による生体磁気計測装置を例に採って説明したが、脳波計などを用いた生体信号計測装置であっても同様の効果を得ることができる。
【００５７】
【発明の効果】
以上に詳述したように、この発明の生体信号計測装置によれば、所謂ＩＣＡ（独立成分分析）手法に従って、複数個のセンサにより検出された原検出信号（観測信号）は、生体活動の各電流源とその他の各信号源ごとの独立成分に分解された上で、各独立成分について各独立成分自体の状態のみに基づきノイズ成分であるか否かが判定され、ノイズ成分と判定された独立成分が除かれた後、残りの非ノイズ独立成分によって検出信号が各々の情報源毎に分離して復元される。したがって、ノイズ検出専用のセンサを別途に設けずとも、生体信号測定用のセンサで得られた信号だけで、ノイズ成分が十分に除去された生体信号解析用の復元検出信号が容易に得られる。その結果、復元生体検出信号に基づいて行われる信号解析も正確なものとなる。
【００５８】
加えて、この発明の生体信号計測装置によれば、独立成分に分解される前に、原検出信号における生体活動電流源による生体信号成分とノイズ成分との間が無相関化される構成を備えているので、大きなノイズ成分であっても、独立成分の分解の際にノイズ成分が生体信号成分の側に紛れ込み難くなる。また、分解後は単独の成分のみからなるはずの独立成分が複数の独立成分を含む未独立状態となる心配がなくなる。さらに、予め求出しておいた独立成分の個数に等しい個数の独立成分に原検出信号が分解されるので、センサの個数が発生源全個数より大幅に多い場合であっても、検出信号の分解個数が事前に十分に絞られていて、分解後は一つの独立成分だけに纏まって含まれているはずの成分が複数の独立成分に分割されてしまう心配がなくなる。したがって、センサの個数が発生源全個数より大幅に多い場合や、ノイズレベルが大きい場合でも、検出信号の独立成分への分解は適正になされることから、正確な信号解析が行われることとなる。
【図面の簡単な説明】
【図１】実施例に係る生体磁気計測装置の全体構成を示すブロック図である。
【図２】実施例装置による生体磁気の計測実行の様子を経時的に示すフローチャートである。
【図３】実施例装置により得られたＩＣＡ信号Ｘの一部を示すグラフである。
【図４】実施例装置により解析された後頭部側大脳内の生体活動電流源の発生状況を示す模式図である。
【図５】実施例装置により計測された後頭部側大脳内の特定の９ポイントに対応する解析電流源の時間変化を示すグラフである。
【図６】従来装置により解析された後頭部側大脳内の生体活動電流源の発生状況を示す模式図である。
【図７】従来装置により計測された後頭部側大脳内の特定の９ポイントに対応する解析電流源の時間変化を示すグラフである。
【符号の説明】
１ …マルチチャンネルＳＱＵＩＤセンサ
１ａ …ＳＱＵＩＤセンサ
２ …データ変換収集部
３ …信号処理部
４ …予備処理部
５ …成分間相互依存解消部
６ …磁気信号分解部
７ …ノイズ成分除去部
８ …磁気信号復元部
９ …磁気解析部
Ｍ …被検体
Claims (1)
- 被検体の診断対象領域において生体活動電流源により生じる微小な生体信号を測定する複数個のセンサを備えた生体信号計測装置において、前記複数個のセンサにより検出された検出信号における独立成分の個数を求出する独立成分個数求出手段と、前記検出信号に乗算処理を行って、前記検出信号を
（ａ）前記複数個のセンサにより検出された検出信号において生体活動電流源による生体信号成分と
（ｂ）ノイズ成分
との和とすることにより、前記生体信号成分と前記ノイズ成分との相互依存性を解消する（前記生体信号成分と前記ノイズ成分との間を無相関化する）成分間相互依存解消手段と、前記生体信号成分と前記ノイズ成分との相互依存性が解消された検出信号を前記独立成分個数求出手段により求出した個数の数の独立成分に分解する信号分解手段と、この信号分解手段により分解された独立成分の中のノイズ成分を独立成分自体の状態のみに基づいて判定して取り除くノイズ成分除去手段と、各非ノイズ独立成分に基づき検出信号を復元する信号復元手段と、復元された各検出信号の各々の独立成分に対応する脳内の活動（位置、方向、強さ）を得る信号解析手段とを備えていることを特徴とする生体信号計測装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP30413299A JP4236352B2 (ja) | 1999-10-26 | 1999-10-26 | 生体信号計測装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP30413299A JP4236352B2 (ja) | 1999-10-26 | 1999-10-26 | 生体信号計測装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001120511A JP2001120511A (ja) | 2001-05-08 |
| JP4236352B2 true JP4236352B2 (ja) | 2009-03-11 |
Family
ID=17929435
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP30413299A Expired - Fee Related JP4236352B2 (ja) | 1999-10-26 | 1999-10-26 | 生体信号計測装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4236352B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102540130A (zh) * | 2012-01-10 | 2012-07-04 | 中国科学院半导体研究所 | 应用于电流传感器角差在线监测***的信号分离方法 |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4352315B2 (ja) * | 2002-10-31 | 2009-10-28 | 日本光電工業株式会社 | 信号処理方法／装置及びそれを用いたパルスフォトメータ |
| JP4807598B2 (ja) * | 2002-10-31 | 2011-11-02 | 日本光電工業株式会社 | パルスフォトメータ |
| 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 |
| US7025728B2 (en) * | 2003-06-30 | 2006-04-11 | Nihon Kohden Corporation | Method for reducing noise, and pulse photometer using the method |
| JP4876988B2 (ja) * | 2007-03-13 | 2012-02-15 | 日本電気株式会社 | ブレイン・コンピュータ・インタフェース装置 |
| JP2012146116A (ja) * | 2011-01-12 | 2012-08-02 | Kyushu Institute Of Technology | スピーチ内容を識別する装置及び方法 |
| WO2017143319A1 (en) * | 2016-02-19 | 2017-08-24 | The Regents Of The University Of California | Systems and methods for artifact reduction in recordings of neural activity |
| JP6785446B2 (ja) | 2016-04-14 | 2020-11-18 | パナソニックＩｐマネジメント株式会社 | 生体信号計測システム |
| CA3055631A1 (en) | 2017-03-15 | 2018-09-20 | The Regents Of The University Of California | Removal of stimulation artifact in multi-channel neural recordings |
| JP6996095B2 (ja) * | 2017-03-17 | 2022-01-17 | 株式会社リコー | 情報表示装置、生体信号計測システムおよびプログラム |
-
1999
- 1999-10-26 JP JP30413299A patent/JP4236352B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102540130A (zh) * | 2012-01-10 | 2012-07-04 | 中国科学院半导体研究所 | 应用于电流传感器角差在线监测***的信号分离方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2001120511A (ja) | 2001-05-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4068763B2 (ja) | 生体信号計測装置 | |
| US7729740B2 (en) | Noise cancellation in magnetoencephalography and electroencephalography with isolated reference sensors | |
| Ghuman et al. | A wavelet-based method for measuring the oscillatory dynamics of resting-state functional connectivity in MEG | |
| Van Dijk et al. | Prestimulus oscillatory activity in the alpha band predicts visual discrimination ability | |
| JP4236352B2 (ja) | 生体信号計測装置 | |
| US6195576B1 (en) | Quantitative magnetoencephalogram system and method | |
| Larson et al. | Reducing sensor noise in MEG and EEG recordings using oversampled temporal projection | |
| Kruggel et al. | Recording of the event‐related potentials during functional MRI at 3.0 Tesla field strength | |
| EP2624748B1 (en) | Device for use in electro-biological signal measurement in the presence of a magnetic field | |
| Yamasaki et al. | Effects of distraction on pain-related somatosensory evoked magnetic fields and potentials following painful electrical stimulation | |
| JP3473210B2 (ja) | 生体磁気計測装置 | |
| Jaušovec et al. | Differences in EEG current density related to intelligence | |
| Tang et al. | Independent components of magnetoencephalography: single-trial response onset times | |
| JP4804647B2 (ja) | 磁気的核スピン共鳴により生物の機能的脳活動を時間および位置的に分解して表示する方法および装置 | |
| Ichikawa et al. | Subjective and model-estimated reward prediction: Association with the feedback-related negativity (FRN) and reward prediction error in a reinforcement learning task | |
| US20120043969A1 (en) | Impedance Tomography Apparatus | |
| Yaple et al. | Feedback-related negativity reflects omission of monetary gains: evidence from ERP gambling study | |
| Wang et al. | Detection of P300 brain waves using a Magneto-Impedance sensor | |
| Westin et al. | Detection of interictal epileptiform discharges: A comparison of on-scalp MEG and conventional MEG measurements | |
| Takemi et al. | Is event-related desynchronization a biomarker representing corticospinal excitability? | |
| Tuunanen et al. | Comparison of BOLD fMRI and MEG characteristics to vibrotactile stimulation | |
| Tervaniemi et al. | Test–retest stability of the magnetic mismatch response (MMNm) | |
| US20210186400A1 (en) | Apparatuses, systems, and methods for suppression of artifacts in non-invasive electromagnetic recordings | |
| Hamada et al. | Effect of interstimulus interval on attentional modulation of cortical activities in human somatosensory areas | |
| Anurova et al. | Cortical generators of slow evoked responses elicited by spatial and nonspatial auditory working memory tasks |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20040303 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060621 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080623 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080729 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080926 |
|
| 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: 20081209 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20081216 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111226 Year of fee payment: 3 |
|
| LAPS | Cancellation because of no payment of annual fees |