(57)【要約】ウエーブレット変換を用いて複合信号中のノイズ寄与を減少する２部法が説明されている。この方法は、ダイナミックなプロセス（決定論的ノイズまたは信号自体）により形成された再構成された状態空間内のサブ空間の特定と、種々のサブ空間の分離と、確率的ノイズからのサブ区間の分離とを含む。この方法は、決定論的構造を特定および記述するために、ウエーブレット変換を用いて再構成された状態空間内で局部的投射を行うことにより磁気心電計時系列または心電計時系列の非直線ノイズ除去（ＮＬＤ）を行うために用いられる。任意の決定論的プロセスにより発生されたサブ空間の場所が定められ、それの源とは独立に分離される。
2005503855.tif 000002
【特許請求の範囲】
【請求項１】
決定論的な構造を特定および記述するためにウエーブレット変換を用いて局部的投射を再構成された状態空間内で実行することを備える、磁気心電計または心電計の時系列信号のノイズを非線形的に除去する方法。
【請求項２】
状態空間内で磁気心電計または心電計に関連するサブ空間を特定および記述するためにウエーブレット変換を適用することを備える、磁気心電計または心電計の時系列信号のノイズを非線形的に除去する方法。
【請求項３】
磁気心電計または心電計からの多様な信号から決定論的ノイズを除去する方法であって、
信号を主信号センサと、主信号センサとは別の少なくとも１つの基準センサとに記録することと、
Ｔａｋｅｎｓ理論を用いて、測定された時系列に働きかける状態空間を再構成することと、
状態空間ベクトルをウエーブレット基準システムに変換することと、
基準センサの時系列中のノイズの状態空間ベクトルの特徴を特定し、この特徴を源（信号）センサの時系列内に再び配置することと、
状態空間内のノイズの状態空間ベクトルを差し引くことと、
挾雑物を除かれた磁気心電計信号を再構成することと、
を備える磁気心電計または心電計からの多様な信号から決定論的ノイズを除去する方法。
【請求項４】
信号の源が磁気心電計装置または心電計装置内に適切に配置されている適当なセンサからのものである請求項１に記載の方法。
【請求項５】
複合信号のウエーブレット特定ステップと複合信号の分離ステップとの少なくとも一方が個別の双直交および連続ウエーブレットに対して使用される請求項１に記載の方法。
【請求項６】
使用されるウエーブレットが、Ｈａａｒウエーブレットタイプと、Ｍｏｒｌｅｔウエーブレットタイプと、Ｍｅｘｉｃａｎ  Ｈａｔウエーブレットタイプと、双直交スプライン・ウエーブレットタイプと、Ｄａｕｂｅｃｈｉｅｓウエーブレットタイプと、Ｍａｌｖａｒウエーブレットタイプと、Ｌｅｍａ ｒｉｅウエーブレットタイプと、Ｃｏｉｆｌｅｔウエーブレットタイプと、Ｍｅｙｅｒウエーブレットタイプと、Ｓｙｍｌｅｔウエーブレットタイプとより成る群から選択される請求項１に記載の方法。
【請求項７】
使用されるウエーブレットタイプがＣｏｉｆｌｅｔウエーブレットタイプである請求項６に記載の方法。
【請求項８】
使用されるウエーブレットタイプがフィルタ次数６を用いるＣｏｉｆｌｅｔウエーブレットタイプである請求項６に記載の方法。
【請求項９】
複数の複合信号の特定と分離との少なくとも一方を、信号をいくつかのサブインターバル信号に分割し、回帰ウエーブレット変換を各サブインターバル信号に適用することにより行う請求項１に記載の方法。
【請求項１０】
確率的かつ決定論的ノイズを含んでいる磁気心電計または心電計からの標本化された複合信号をそれの信号成分とノイズ成分とに分離する方法であって、
（ａ）信号を主信号センサと、主信号センサとは別の少なくとも１つの基準センサとに記録することと、
（ｂ）サブ空間を確率的ノイズから分離することと、
（ｃ）個々のダイナミック・プロセスに属している種々のサブ空間を分離することと、
（ｃ）状態空間ベクトルをウエーブレット基準システムに変換することと、
（ｄ）源センサの時系列中のノイズの状態空間ベクトルの特徴を特定することと、
（ｅ）再構成されたノイズを信号時間系列から差し引くことと、
（ｆ）挾雑物を除かれた磁気心電計信号を再構成することと、
を備える確率的かつ決定論的ノイズを含んでいる磁気心電計または心電計からの標本化された複合信号をそれの信号成分とノイズ成分とに分離する方法。
【請求項１１】
ＳＱＵＩＤ磁気心電計装置から受けた信号のノイズを除去するための請求項１０に記載の方法。
【請求項１２】
遮蔽されていない磁気心電計からの信号のノイズを除去するための請求項１に記載の方法。
【請求項１３】
不整脈と、心筋梗塞と、虚血とで構成されている群から選択された心臓疾患によってつくられた状態空間内にサブ空間が存在することを決定するための請求項１に記載の方法の使用方法。
【請求項１４】
携帯型心電計を用いて２４時間記録するための請求項１に記載の方法の使用方法。
【請求項１５】
母と胎児との磁気心電計からの信号を分離するための請求項１に記載の方法の使用方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は全体として磁気心電計検査および心電計検査の分野に関する。更に詳しくいえば、本発明はウエーブレット変換を用いる磁気心電図の非線形ノイズ低減に関する。
【背景技術】
【０００２】
磁気心電計検査（ＭＣＧ）は心臓の筋肉の電気的に活性な細胞により、小さい電流から心臓によって発生された磁界を測定するものである。それはありふれた臨床診断にまだ導入されていない非侵襲診断法である。
【０００３】
磁気心電計検査は、心臓内の電気生理学的過程により胴（または胎児の磁気心電計検査における母の腹部）の上部に生じた時間的に変化する磁界の測定より成る。胴の上部のそれらの磁界の測定によって、心電計検査により得られる胴の上部のそれらの磁界の測定を補う情報が得られ、それは心臓の機能の異常診断に特に用いられる。
【０００４】
それらの信号の強度が極めて弱い（地磁気の１００万分１またはそれ以下）ために、現在は超伝導量子干渉素子（ＳＱＵＩＤｓ）のみがそのような作業を行える。そのようなセンサは非常に低い温度（極低温）でのみ動作するので、特殊な容器（低温槽）内に置かねばならない。それの外壁は室温にあり、内部は、通常は低温度液体（冷却剤）を充填することにより、低温にされる。ＳＱＵＩＤセンサが設けられている低温槽は人体に接触することなしに、人体に近接して配置される。原理的には、高温超伝導ＳＱＵＩＤｓ（ＨＴＳＱＵＩＤｓ）およびグラジオメーター（ｇｒａｄｉｏｍｅｔｅｒ）を利用することにより、液体ヘリウムの代わりに液体窒素を使用することが可能であり、これは取扱いを極めて簡単にするが、今日まで発表されているＨＴＳを基礎とするＭＣＧ装置は完全には実用になっていない。
【０００５】
心電計検査と比較して磁気心電計検査にはいくつかの利点があるが、実際的な臨床使用のための飛躍的な進歩は依然としてなされていない。したがって、医師を納得させる魅力的な結果をもたらし、ＳＱＵＩＤシステムのコストを低減することが必要である。両方とも改良したノイズ打ち消し法を基にして達成できる。
【０００６】
非常に強力な自然に存在する外部信号と人が発生する外部信号が存在することによって信号対ノイズ比の値が極めて低く（１００万分の１より低い）なる。それらの信号は、外部信号（ノイズ）抑制法が採用されなければ、使用できない。
【０００７】
最も効果的ではあるが、最も費用を要しかつ硬直したノイズ抑制法は高度に磁気遮蔽された部屋の内部でこの方法を実施することである。しかし、そのような部屋は心臓病学実地診療には受け容れることができないことが判明している。磁気遮蔽の外部で、高度に平衡されたＳＱＵＩＤグラジオメーター装置なしで動作することが臨床での使用には重要である。したがって、最近の努力は、磁気遮蔽された部屋なしに使用できるＭＣＧ装置の開発に向けて傾注されている。使用されている主なノイズ抑制技術は高次のグラジオメトリー（ｇｒａｄｉｏｍｅｔｒｙ）である。
【０００８】
診断応用では、心筋の病的な機能変化または器質的変化に起因する心臓の電気的活動の複雑な変化を特定するために、心臓の磁界を空間的と時間的との少なくとも一方で解析できる。それらは虚血、心筋梗塞、心室内の容積または圧力の変化、または不整脈から起きることがある。
【０００９】
ＳＱＵＩＤセンサのアレイによる磁気心電計画像形成は虚血や心筋の活性（休止している組織と壊死している組織の識別）の診断と、不整脈の危険の分析とに使用するためにますます研究されている。興奮伝導系における病巣活動を識別するために、生物磁気発生場所特定を心臓学において使用できる。特に、Ｗｏｌｆ‐Ｐａｒｋｉｎｎｓｏｎ‐Ｗｈｉｔｅ症候群 におけるように付属通路、心室性期外収縮または心室性頻脈の源の場所をミリメートルの精度で非侵襲的に定めることができる。
【００１０】
ＭＣＧの潜在的な重要性は、それが完全な非侵襲の非接触診断および機能的画像形成法であって、その方法では、何百人もの心臓動脈疾患患者を対象としたいくつかの臨床研究において非常に高い感度および特効性が示されている。
【００１１】
磁気遮蔽の外で測定された磁気心電計は心臓の信号に重畳された環境ノイズの悪影響を受ける。３つの種類のノイズ、すなわち、一様なノイズ（たとえば、地球の磁界）、確率的ノイズ（白色ノイズ、狭スペクトラムノイズ（ｃｏｌｏｒｅｄ  ｎｏｉｓｅ）、１／ｆノイズ）、および決定論的ノイズ（たとえば、電力スペクトラムにおけるピークが５０／６０Ｈｚにある電力線擾乱）、を識別できる。一様なノイズ成分および決定論的ノイズ成分は信号より何桁も大きいことがしばしばある。また、確率的ノイズおよび決定論的ノイズは時間的に変化するので適応ノイズ打ち消しが求められる。
【００１２】
決定論的ノイズ成分は低周波、中間周波または高周波のことがある。低周波の決定論的ノイズ（０．１ないし１Ｈｚ）は運転中のエレベーター、動いている金属製ドア、動いている金属製椅子、またはその他の動いている金属（磁性）物体に通常起因する。除細動器、ペースメーカー、胸骨ワイヤまたは歯科治療材料などの磁性埋め込み物は患者の呼吸数で振動することがある。呼吸は磁性部材を動かし、その結果として通常は大きな振幅の心臓時系列がずれることになる。更に、人体内の磁性部材は心臓の機械的なポンプ作用によって震動することがある。そうするとその震動の数は心拍に強く相関させられて、「相関ノイズ」と一般に呼ばれているものが生ずることになる。
【００１３】
中間周波の決定論的ノイズ（１ないし２０Ｈｚ）は回っているファン、空調機、またはその他の臨床機器によって通常発生される。建物や装置自体の震動や磁束の急激な変化もこの中間周波範囲内で擾乱をひき起こすことがある。
【００１４】
高周波ノイズ（＞２０Ｈｚ）は電源、モニタ周波数、またはその他の電子装置に起因するのがほとんどである。
【００１５】
それら種々の決定論的ノイズ源は、磁気心電計検査解析のために求められる有 用で歪みのない磁気心電計検査［ＭＣＧ］データを取出すことを困難にする。
【００１６】
そのような望ましくないノイズを減少または信号から除去することを試みて多くの技術が利用されてきた。利用される最も一般的なノイズ減少法は、低域フィルタ、高域フィルタまたはノッチフィルタを用いる旧式の濾波に組合わせたハードウェアおよびソフトウェアのグラジオメトリー技術である。
【００１７】
各種の旧式フィルタはこの領域では性能が良くなかった。フィルタは適応的ではなく、それらを使用すると、特にノッチフィルタの場合に、信号の保持が不十分になる結果となる。
【００１８】
一次または高次のハードウェア・グラジオメトリーがより低い次数の一様な磁界または傾斜磁界を抑制するために利用されてきた。この方法は、たとえば、一様な地球磁界の影響を効率的に減少し、心臓の信号にほんの僅かな影響を及ぼすだけである。しかし、近くの源から生じ、空間的に大きく傾いている決定論的ノイズ成分および確率的ノイズ成分は、高精度の高次グラジオメーターによってさえも十分には抑制されない。また、そのグラジオメーターは製造が困難であり、したがって高価である。
【００１９】
決定論的ノイズを除去するために最も成功した先行技術の方法は多数の基準センサを使用することである。相互相関技術を種々のやり方で適応的に適用することにより（Ｒｏｂｉｎｓｏｎ、「生体磁気測定のための環境ノイズ打ち消し（Ｅｎｖｉｒｏｎｍｅｎｔａｌ  Ｎｏｉｓｅ  Ｃａｎｃｅｌｌａｔｉｏｎ  Ｆｏｒ  Ｂｉｏｍａｇｎｅｔｉｃ  Ｍｅａｓｕｒｅｍｅｎｔｓ）」（１９８９）、およびＲｕｅｄｅｒｓ他、「周波数依存グラジオメトリー：磁気心電計検査に応用した新規な改良した非侵襲ノイズ打ち消し法（Ｆｒｅｑｕｅｎｃｙ  Ｄｅｐｅｎｄｅｎｔ  Ｇｒａｄｉｏｍｅｔｒｙ：Ａ  Ｎｅｗ  Ｎｏｎ‐Ｉｍｖａｓｉｖｅ  Ｍｅｔｈｏｄ  Ｏｆ  Ｉｍｐｒｏｖｅｄ  Ｎｏｉｓｅ  Ｃａｎｃｅｌｌａｔｉｏｎ  Ａｐｐｌｉｅｄ  Ｔｏ  Ｍａｇｎｅｔｏｃａｒｄｉｏｇｒａｐｈｙ）（１９８９）」）、ノイズのピークが相関させられているならば、ノイズピークを信号センサから差し引くことが可能である。この文脈では、相関させられているというのは（決定論的）ノイズが自己相関させられているが、信号とは相関させられていないことを意味する。
【００２０】
多重センサ技術での問題は、十分なノイズ傾き抑制のためには、少なくとも７個、および２５個までの基準センサを必要とすることである。
【００２１】
更に、多重基準センサは、相互相関信号処理に結び付けられたとしても、確率的ノイズの問題である、信号識別および解析における大きな問題を解決できない。確率的ノイズは、全く相関しないので、多数の基準センサ手法で取扱える。
【００２２】
確率的ノイズを信号から減少または除去するために多くの試みが行われてきた。状態空間内での局部的投射と共分散マトリックスを用いる方法（「心電図用の非線形ノイズ減少（Ｎｏｎｌｉｎｅａｒ  Ｎｏｉｓｅ  Ｒｅｄｕｃｔｉｏｎ  Ｆｏｒ  Ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｍｓ）」と題するＳｃｈｒｅｉｂｅｒ他の論文、Ｃｈａｏｓ  ６：８７、（１９９５）におけるように）が、確率的ノイズの減少に有用であることが示されている。この方法では、信号の特徴が状態空間内の場所に配置されて、共分散マトリックスの最大の固有値により示されているノイズなしサブ空間に投射される。この方法は良く機能するが、それは状態空間内の信号のサブ空間の次元が既知である時のみである。
【００２３】
一般に、ノイズレベルが高い場合には、状態空間内の信号のサブ空間の次元は既知ではなく、固有値のスペクトラムは平坦である。
【００２４】
磁気心電計検査（ＭＣＧ）では、非線形ビーム形成器などの空間濾波に対する数学的手法、特に合成開口磁力計法、が、Ｓ．Ｅ．ＲｏｂｉｎｓｏｎおよびＪ．Ｖｒｂａ、平均されていないＭＥＧのためのＳＡＭおよびＭＵＳＩＣ性能の比較（Ｃｏｍｐａｒｉｓｏｎ  ｏｆ  ＳＡＭ  ａｎｄ  ＭＵＳＩＣ  Ｐｅｒｆｏｒｍａｎｃｅ  ｆｏｒ  Ｕｎａｖｅｒａｇｅｄ  ＭＥＧ）と、フィンランドＨＵＴ、Ｅｓｐｏｏにおいて開催された、第１２回国際生体磁気会議（１２ｔｈ  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｂｉｏｍａｇｎｅｔｉｓｍ）の生体磁気２０００議事録（Ｐｒｏｃｅｅｄｉｎｇｓ  ｏｆ  Ｂｉｏｍａｇ  ２０００）、所載の、Ｊ．ＶｒｂａおよびＳ．Ｅ．Ｒｏｂｉｎｓｏｎ、合成開口磁力計法と線形ビーム形成器の違い（Ｄｉｆｆｅｒｅｎｃｅｓ  ｂｅｔｗｅｅｎ  Ｓｙｎｔｈｅｔｉｃ  Ａｐｅｒｔｕｒｅ  Ｍａｇｎｅｔｏｍｅｔｒｙ  ａｎｄ  Ｌｉ ｎｅａｒ  Ｂｅａｍｆｏｒｍｅｒｓ）に記述されているように、脳内の電気的活動源および磁気的活動源の場所を特定するために用いられてきた。
【００２５】
そのような方法はノイズ分離のためにも役立つ。しかし、著者らによれば、合成開口磁力計法または類似の方法はＭＣＧには容易に応用されそうもない。主な理由は人の心臓は、脳内の非常に局限された活動源とは異なって、少なくとも心臓サイクルのＱＲＳ期間とＳＴ期間においては、空間的に拡げられた電気源および磁気源を表すからである。
【００２６】
したがって、現在までは、磁気心電計データ中のあらゆる種類のノイズを十分に減少するために満足できる技術は利用できなかった。
【００２７】
この分野におけるいくつかの先行技術の特許および文献について以下に述べる。いくつかの特許は信号からノイズを除去するためにウエーブレット変換を利用する。
【００２８】
ゼネラル・エレクトリック・カンパニーに譲渡された、超音波画像形成における点在ノイズ濾波（Ｓｐｅｃｋｌｅ  Ｎｏｉｓｅ  Ｆｉｌｔｅｒｉｎｇ  ＩｎＵｌｔｒａｓｏｕｎｄ  Ｉｍａｇｉｎｇ）という名称の米国特許第５，４９７，７７７号で、Ａｂｄｅｌ‐Ｍａｌｅｋ他は、ウエーブレット変換を用いて対象とする信号からノイズを濾波する方法を開示している。ノイズ成分が寄与するウエーブレット変換係数のいくつかが無くされ、真の信号に属する係数のみが逆変換される。その逆変換はノイズ成分の無い真の信号に近似するものを回復する。
【００２９】
しかし、この手法は、ノイズと信号が種々の係数により表され、かつ、特に、いくつかの係数に重畳していないという仮定を基にしている。また、ノイズに属している係数のみを除去するためには、どの係数が信号情報を含んでおり、どの係数がノイズを含んでいるかについて知ることが求められる。
【００３０】
カウンシル・オブ・サイエンティフィック＆インダストリアル・リサーチ（Ｃｏｕｎｃｉｌ  ｏｆ  ｓｃｉｅｎｔｉｆｉｃ  ＆  Ｉｎｄｕｓｔｒｉａｌ  Ｒｅｓｅａｒｃｈ）に譲渡された、複素複合信号の識別と、複素複合信号の決定論的成分およびノイズ成分への分離との少なくとも一方を行う方法および装置（Ｍｅｔｈｏｄ  Ａｎｄ  Ａｎ  Ａｐｐａｒａｔｕｓ  Ｆｏｒ  Ｔｈｅ  Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ  Ａｎｄ／Ｏｒ  Ｓｅｐａｒａｔｉｏｎ  Ｉｎｔｏ  Ｉｔｓ  Ｄｅｔｅｒｍｉｎｉｓｉｔｉｃ  Ａｎｄ  Ｎｏｉｓｙ  Ｃｏｍｐｏｎｅｎｔｓ）という名称の米国特許第６，２０８，９５１号で、Ｋｕｍａｒ他が、ウエーブレット変換を用いて対象とする信号からノイズを濾波する方法をまた開示している。ウエーブレット変換の諸特性およびそれの種々の次元を用いてノイズ成分が無くされる前に、複合信号がウエーブレット変換されて真の信号とノイズ信号を分離し、所望の信号を回復する。
【００３１】
この手法の問題は、ウエーブレット変換を実行する前に信号とノイズを分離することを要することである。これは、測定されるＭＣＧ時系列ではそうではない。したがって、ウエーブレット変換を実行するために信号とノイズを予め分離することを要しない技術を要する。必要とするものは、ウエーブレット変換を適用すると所望の分離（鋭い固有スペクトラム）を行えるようにするやり方で時系列を再編成する技術である。
【００３２】
フォード・グローバル・テクノロジーズ社（Ｆｏｒｄ  Ｇｌｏｂａｌ  Ｔｅｃｈｎｏｌｏｇｉｅｓ，Ｉｎｃ．）に譲渡された、信号中のパルス状成分および非パルス状成分の分離方法および装置（Ｍｅｔｈｏｄ  Ａｎｄ  Ａｐｐａｒａｔｕｓ  Ｆｏｒ  Ｓｅｐａｒａｔｉｏｎ  Ｏｆ  Ｉｍｐｕｌｓｉｖｅ  Ａｎｄ  Ｎｏｎ‐Ｉｍｐｕｌｓｉｖｅ  Ｃｏｍｐｏｎｅｎｔｓ  Ｉｎ  Ａ  Ｓｉｇｎａｌ）という名称の米国特許第６，２４９，７４９号で、Ｔｒａｎ他が、ウエーブレット変換係数に対して統計解析を行い、種々の信号へのそれの寄与を検出することにより、複合信号内の２つの信号を分離する方法を開示している。いずれかの信号に寄与する係数は、各信号を個々に回復するために別々に逆変換される。
【００３３】
Ｋｕｍａｒ他のように、Ｔｒａｎは信号成分とノイズ成分を予め分離することを要する。それらの特許はノイズからの信号の分離のためにウエーブレット変換を使用する最も簡単なやり方を反映している。
【００３４】
非特許文献にもノイズ減少技術が開示されている。Ｌ．Ｒｅｂｏｌｌｏ‐Ｎｅｉｒａ、Ａ．Ｃｏｓｔａｎｔｉｎｉｄｅｓ、Ｔ．Ｓｔａｔｈａｋｉが「フレームをベースとするウエーブレットによる圧縮およびノイズ減少のための信号表現（Ｓｉｇｎａｌ  Ｒｅｐｒｅｓｅｎｔａｔｉｏｎ  Ｆｏｒ  Ｃｏｍｐｒｅｓｓｉｏｎ  Ａｎｄ  Ｎｏｉｓｅ  Ｒｅｄｕｃｔｉｏｎ  Ｔｈｒｏｕｇｈ  Ｆｒａｍｅ‐ Ｂａｓｅｄｅ  Ｗａｖｅｌｅｔｓ）」、ＩＥＥＥ  Ｔｒａｎｓ．Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ  ４６（３）：５８７−５９７（１９９８）において、ウエーブレット変換を用いるノイズ減少法について論じている。この論文は、信号がノイズで汚染されていることがあると仮定されているくつかのウエーブレットサブ空間の抑制について述べている。
【００３５】
Ｌｅｄｅｒ他の「遮蔽されていない環境におけるＨＴＳ‐ＳＱＵＩＤ磁気心電計検査の再生可能性（Ｒｅｐｒｏｄｕｃｔｉｖｉｔｙ  ｏｆ  ＨＴＳ‐ＳＱＵＩＤ  Ｍａｇｎｅｔｏｃａｒｄｉｏｇｒａｐｈｙ  ｉｎ  ａｎ  Ｕｎｓｈｉｅｌｄｅｄ  Ｅｎｖｉｒｏｎｍｅｎｔ）」が、高温超伝導（ＨＴＳ）センサを用いて人の心臓の磁界を測定する技術を開示している。それらのセンサは液体窒素の温度で電磁遮蔽なしで動作させられる。この論文は依然として見付かっていない万能ノイズ打ち消し技術に対する需要を強調している。
【００３６】
ＨＴＳ  ＳＱＵＩＤ技術は遮蔽の外部の磁気心電計を測定するためにはまだ適切ではない。いくつかの有望な結果は存在するが、高温超伝導体は低温超伝導体と比較して感度が低い（４〜５分の１）。これはシステムの性能を常に低下するので、拍動の磁気的特徴（ｓｉｇｎａｔｕｒｅ）における詳細はＨＴＳ装置では現れない。それは胎児のＭＣＧに対しては一層悪い。その理由は胎児のＭＣＧの磁界の強さが成人におけるそれより少なくとも１桁弱いからである。
【００３７】
Ｋｏｃｈが、ＳＱＵＩＤ磁気心電計検査：現状および見通し（ＳＱＵＩＤ  ｍａｇｎｅｔｏｃａｒｄｉｏｇｒａｐｈｙ：Ｓｔａｔｕｓ  ａｎｄ  ｐｅｒｓｐｅｃｔｉｖｅｓ）、ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ａｐｐｌｉｅｄ  Ｓｕｐｅｒｃｏｎｄｕｃｔｉｖｉｔｙ  １１：（１）４９−５１、第１部（２０００年３月）で、改良したノイズ抑制技術などのＳＱＵＩＤシステムにおける最近の進歩について詳述し、より高い磁界感度（特にＨＴＳＱＵＩＤｓに対して）、リアルタイム選択、ベクトル磁力計および新規な信号解析手法が、実用的な臨床使用への磁気心電計検査技術の実現を最近まで妨げてきた技術的制約を顕著に低減した。この論文はＳＱＵＩＤ磁気心電計検査における現在の状態を要約している。
【００３８】
Ｚｈａｎｇ他が、遮蔽されていない環境内での磁気心電計検査のための第２種 、高温超伝導グラジオメーター（Ｓｅｃｏｎｄ‐Ｏｒｄｅｒ，Ｈｉｇｈ‐Ｔｅｍｐｅｒａｔｕｒｅ  Ｓｕｐｅｒｃｏｎｄｕｃｔｉｎｇ  Ｇｒａｄｉｏｍｅｔｅｒ  Ｆｏｒ  Ｍａｇｎｅｔｏｃａｒｄｉｏｇｒａｐｈｙ  Ｉｎ  Ｕｎｓｈｉｅｌｄｅｄ  Ｅｎｖｉｒｏｎｍｅｎｔ）、Ａｐｐｌｉｅｄ  Ｐｈｙｓｉｃｓ  Ｌｅｔｔｅｒｓ  ７６：（７）９０６−９０８、２０００年２月１４日で、遮蔽されていない環境内での磁気心電計検査のための第２種、高温超伝導グラジオメーターを開示している。この高温ＳＱＵＩＤシステムは信号対ノイズ比、空間分解度、周波数帯域幅、環境擾乱の除去、および長期間安定性等についての諸考察の面で磁気心電計検査のために診断的に関連するものであることが示されている。Ｚｈａｎｇは、患者のベッドで直接使用できる、携帯可能なデュワー内の遮蔽なし単一チャネルシステムを開示している。低温超伝導体ＳＱＵＩＤ性能と比較して、それは非常に弱い。しかし、それの性能は梗塞の患者におけるＳＴ‐セグメントに対するそれの狭い意図された使用には十分であるかもしれない。
【００３９】
Ｒｏｂｉｎｓｏｎが、生体磁気測定のための環境ノイズ打ち消し（Ｅｎｖｉｒｏｎｍｅｎｔａｌ  Ｎｏｉｓｅ  Ｃａｎｃｅｌｌａｔｉｏｎ  Ｆｏｒ  Ｂｉｏｍａｇｎｅｔｉｃ  Ｍｅａｓｕｒｅｍｅｎｔｓ）、生体磁気における進歩（Ａｄｖａｎｃｅｓ  ｉｎ  Ｂｉｏｍａｇｎｅｔｉｓｂ）、Ｐｌｅｎｕｍ  Ｐｒｅｓｓ，Ｎｅｗ  Ｙｏｒｋ  １９８９で，生体磁気のノイズ除去技術の現状についての全般的な説明を行っている。この論文は基準センサの使用と、相互相関技術を基にしたノイズ打ち消しについて記述している。
【００４０】
この論文は相互相関ノイズ除去のための主要な典拠である。この文献に記述されている手法はノイズ除去のために現在商業的に利用されている。この論文では、必要な基準センサの最少数は７であることが判明している。
【００４１】
ソフトしきい値形成法によるノイズ除去（Ｄｅｎｏｉｓｅ  Ｂｙ  Ｓｏｆｔ‐Ｔｈｒｅｓｈｈｏｌｄｉｎｇ），ＩＥＥＥ  Ｔｒａｎｓ．Ｉｎｆｏｒｍ．Ｔｈｅｏｒｙ  ４１：６１３，（１９９５）、はウエーブレット係数のハードしきい値またはソフトしきい値が、Ｅｆｆｅｒｎ他、「事象に関連する可能性へ応用する過渡信号の非線形ノイズ除去（Ｎｏｎｌｉｎｅａｒ  Ｄｅｎｏｉｓｉｎｇ  Ｏｆ  Ｔｒａｎｓｉｅｎｔ  Ｓｉｇｎａｌｓ  Ｗｉｔｈ  Ａｐｐｌｉｃａｔｉｏｎ  Ｔｏ  Ｅｖｅｎｔ  Ｒｅｌａｔｅｄ  Ｐｏｔｅｎｔｉａｌｓ）」、Ｐｈｙｓｉｃｓ  Ｄ  １４０（３−４）、６月１５日（２０００）により記述されているように、状態空間内でさえも信号回復に良く適する。
【００４２】
この論文はＥＥＧ時系列から過渡事象を除去することを提案している。患者に刺激を加えることにより事象に関連する可能性（ＥＲＰｓ）がひき出される。脳内の対応する領域が特定の応答波形を示す。その波形はそれの極性および事象後の時間に従って、分類される。Ｅｆｆｅｒｎは、信号対ノイズ比（ＳＮＲ）が１よりはるかに小さい非常に弱い波であるＰ３００を解析した。Ｐ３００は通常は何ミリ秒かの間だけ生ずるのでノイズ除去は非常に困難である。
【００４３】
Ｅｆｆｅｒｎの重要な概念はいわゆる円形埋め込みである。彼は、彼が単一のＰ３００時系列の全てを１つの「大きな」時系列に連続して加え合わせることによりつくった人工時系列を埋め込むためにＴａｋｅｎの理論を用いた。埋め込まれたベクトルのウエーブレット変換は過渡事象を彼が識別することを支援し、その後で彼はそれを除去した。
【００４４】
上記先行技術の明確な長所、特徴および利点がどのようなものであっても、先行技術のいずれも磁気心電計データ中のあらゆる種類のノイズを十分に減少するための共通の技術を開示していない。
【００４５】
それの他の用途に加えて、胎児の磁気心電計検査に使用するためにそのような手法を開発することが極めて望ましい。胎児の磁気心電計検査は代わりの胎児監視法としての可能性を持つ。胎児の心臓信号は成人のそれの１０分の１弱いので、より高い磁界分解度が求められる（成人では＜１０ｆｔ／Ｈｚ1/2対＜５０ｆｔ／Ｈｚ1/2）。幸いなことに、２５Ｈｚというようなかなり限られた信号帯域幅で通常十分である。
【００４６】
これまでは、磁気遮蔽されている部屋（ＭＳＲ）の内部での胎児の磁気心電計検査のみが文献に納得できるようにして示され、報告されてきた。遮蔽なしでグラジオメーター、特にＨＴＳグラジオメーター、を使用する試みは、これまで、あまり成功しなかった。妊娠の第３の３か月では、腹腔リードを用いて電気的活動を測定することは確実には可能ではない。その理由は、この期間中は胎児上に電気絶縁層、ｖｅｒｎｉｘ  ｃａｓｅｏｓａ、が存在するからである。磁界は体 組織を通じて余り乱されることなく伝播するので、胎児の磁気心電計検査を胎児のＥＣＧよりも一層精密に記録することが可能である。
【００４７】
胎児の磁気心電計検査は信号の形、心臓時間間隔および心拍可変性を調べるために使用できる。こうすると胎児の心臓興奮伝導系、不整脈、先天性心臓疾患、成長、自律神経系の発達、アシドーシスおよび胎児のストレスを評価できる。
【００４８】
胎児の磁気心電計検査を基にした胎児の心臓診断の現在の状態の概観が、１９９９年９月２６日にドイツ、Ｊｅｎａにおいて開催された第４回Ｈａｎｓ  Ｂｅｒｇｅｒ会議（４th  Ｈａｎｓ  Ｂｅｒｇｅｒ  Ｃｏｎｆｅｒｅｎｃｅ）のサテライト・シンポジウム（ＳＡＴＥＬＬＩＴＥ  Ｓｙｍｐｏｓｉｕｍ）のＡ．Ｌ．Ｐａｓｔｕｚａｒｋ編の胎児の健康最前線における胎児生体磁気（ＦｅｔａｌＢｉｏｍａｇｎｅｔｉｓｍ  ｉｎ  Ｆｒｏｎｔｉｅｒｓ  ｉｎ  Ｆｅｔａｌ  Ｈｅａｌｔｈ）１：（５）に与えられている。
【００４９】
胎児の磁気心電計検査の重要性はそれの独特のモニタ性能および診断性能に存する。胎児の磁気心電計検査の種々の報告された可能な診断用途は２つの適用期間、すなわち、妊娠中と分娩時、に分けることができる。
【００５０】
妊娠中は胎児の磁気心電計検査は、１］特に心臓不整脈または興奮伝導系の擾乱（ＡＶブロック）が疑われる時の、心臓の律動の解析、２］胎児におけるＰＲ間隔の解析および危険状態にある（紅斑性狼瘡、自己免疫疾患、等）胎児個体の第１度ＡＶブロックの診断、３］ＱＲＳコンプレックス（ＱＲＳ  ｃｏｍｐｌｅｘ）の振幅の解析と心室が肥大している胎児（糖尿病の母親、ステロイドを投与されている等の母親の胎児）の診断および追跡診断、４］再分極段階（ｒｅｐｏｌａｒｉｚａｔｉｏｎ  ｐｈａｓｅ）（たとえば、胎児の虚血に関連するＳＴセグメント変化）の解析、５］良好な胎児（心拍数が可変である）の評価、および６］胎児の磁気心電計検査が利用できる唯一の方法であるような長いＱＴ症候群から危険である胎児の検出、に使用できる。
【００５１】
分娩期間中は胎児の磁気心電計検査を１］分娩の種々の段階中に良好な胎児の評価（ＨＲＶ研究）、２］良好である／苦しんでいる胎児についての有用な情報を得るためにＡＶ伝導（ＰＲ間隔）の直接解析、および３］胎児が苦しんでいる間の心臓虚血についての有用な情報を得るためにＳＴセグメント解析、に使用できる。
【００５２】
したがって、本発明の目的は妊娠期間中および分娩期間中に胎児の健康を診断および解析するための手法を得ることである。
【００５３】
また、本発明の目的は、測定された磁気心電計時系列または心電計時系列から決定論的ノイズおよび確率的ノイズをほとんど無くす効果的な装置および方法を得ることである。
【００５４】
本発明の目的は、磁気心電計検査または心電計検査から得られた信号のノイズ除去に特に関して、適応ノイズ打ち消し法を得ることである。
【００５５】
本発明の別の目的は確率的ノイズを除去するためにただ１つの基準センサを利用する適応ノイズ打ち消し法を得ることである。
【００５６】
本発明の目的は状態空間内の信号サブ空間が既知でないような適応ノイズ打ち消し法を得ることである。
【００５７】
本発明の別の目的は、時系列の簡単なウエーブレット変換が、ノイズ成分からの信号の分離を通常妨げるものである平坦な固有スペクトラムを持つような場合の適応ノイズ打ち消し法を得ることである。
【発明の開示】
【課題を解決するための手段】
【００５８】
本発明は、測定された磁気心電計または心電計からの決定論的ノイズまたは確率的ノイズを、既知の先行技術方法よりも効果的にほとんど無くす装置および方法を提供するものである。それは、信号がほぼ決定論的であることを要するのみである。これは、持続時間が４秒またはそれより長い磁気心電計または心電計の時間セグメントが使用されるような場合である。
【発明を実施するための最良の形態】
【００５９】
私は、ウエーブレット係数のハードしきい値またはソフトしきい値が、状態空間内での信号回復に良く適していることを発見し、この技術を磁気心電計時系列信号または心電計時系列信号のノイズ除去に適用した。
【００６０】
本発明は、決定論的構造を特定および記述するためにウエーブレット変換を用いて、再構成された状態空間内で局部的な投射を行うことにより、磁気心電計時系列信号または心電計時系列信号の非線形ノイズ除去（ＮＬＤ）方法および装置を提供するものである。したがって、目標はそれの源（心臓のノイズまたは信号 であろうとも）とは独立に任意の決定論的なプロセスにより発生されたサブ空間の場所を特定し、分離することである。この方法はサブ空間を確率的ノイズからまず分離することより成る。
【００６１】
観察される系のダイナミックな諸特性を表すためには、時間領域の代わりに再構成された状態空間内で動作することが有用である（Ｆ．Ｔａｋｅｎｓ、「擾乱している奇妙な誘因物の検出（Ｄｅｔｅｃｔｉｎｇ  Ｓｔｒａｎｇｅ  Ａｔｔｒａｃｔｏｒｓ  ｉｎ  Ｔｕｒｂｕｌｅｎｃｅ）」、Ｌｅｃｔｕｒｅ  ｎｏｔｅｓ  ｉｎ  ｍａｔｈ．，Ｎｅｗ  Ｙｏｒｋ、１９８１）（図１のａ））。少なくとも近似的な決定論的システムは密に存在している軌道になり、かつサブ空間（図１のｂ））に制約され、一方、確率的プロセスが全状態空間内で無作為に分布させる。決定論的信号に白色ノイズを重畳すると密に存在している主軌道が歪まされる（図１のｃ））。
【００６２】
状態空間内で決定論的構造を特定および記述するためには、状態空間ベクトルを適切な基準システムに変換すると有用である。「適切な」というのは基準関数を見出すことを試みることを意味する。それは存在する特定の決定論的構造に最も良く適応する。
【００６３】
向き情報を圧縮できるという事実のために、新しい基準システムの領域内でほんのわずかな係数によって決定論を記述することが可能である。対照的に、確率的ノイズは圧縮できず、したがって、完全な基準係数のセットを再現することを必要とする。
【００６４】
ウエーブレット変換はウエーブレットと呼ばれる多数の非常に適応的な基準関数を提供する。信号ｘ（ｔ）でたたみこまれた基準関数（ウエーブレット）の翻訳と拡張によりそれは定められる。付加換算係数（ほとんどは２の冪）がウエーブレット変換を数学的顕微鏡のように機能させる；それは種々の尺度にある１つの観察信号情報をそれの場所に依存させる。正確にはこの特性は有用である。その理由はウエーブレット領域内で状態空間ベクトルの多重分解度表現を得るからである（図１のｄ））。
【００６５】
ウエーブレット変換の一般的な数学的背景は周知であって、多重分解度信号分解理論：ウエーブレット表現（Ａ  Ｔｈｅｏｒｙ  Ｆｏｒ  Ｍｕｌｔｉｒｅｓｏ ｌｕｔｉｏｎ  Ｓｉｇｎａｌ  Ｄｅｃｏｍｐｏｓｉｔｉｏｎ：Ｔｈｅ  Ｗａｖｅｌｅｔ  Ｒｅｐｒｅｓｅｎｔａｔｉｏｎ），ＩＥＥＥ  Ｔｒａｎｓ．Ｐａｔ．Ｒｅｃ．Ｍａｃｈ．Ｉｎｔｅｌ．  １１：６７４，（１９８９）、という標題のＭａｌｌａｔ他の論文に概略を見出すことができる。利用できる文献およびウエーブレット応用を包含している包括的なデータベースがＡｍａｒａ  Ｇｒａｐｓウエブサイト：ｈｔｔｐ：／／ｗｗｗ．ａｍａｒａ．ｃｏｍ／ｃｕｒｒｅｎｔ／ｗａｖｅｌｅｔ．ｈｔｍｌ．に提示されている。
【００６６】
最適なウエーブレットを選択することは重要である。最適に選択されたウエーブレット（解析関数）は信号を最も良く表すものである。たとえば、高速フーリェ変換では、解析関数正弦波および余弦波である。純正弦波に適用すると、高速フーリェ変換はスペクトラム中に１つのピークを生ずる。しかし、長方形パターンに高速フーリェ変換を適用するには、このパターンを適切に記述するために膨大な数の係数を要する。このことはウエーブレットでも同じであって、ウエーブレットが対象とする関数（ここでは心臓の鼓動）に良く合致するほど良い。問題指向ウエーブレットを構成することが可能であり、それはウエーブレット変換の大きな利点である。本発明の目的のために、この場合の最良の選択はフィルタ次数６を用いる周知のＣｏｉｆｌｅｔである。Ｈａａｒウエーブレットタイプ、Ｍｏｒｌｅｔウエーブレットタイプ、Ｍｅｘｉｃａｎ  Ｈａｔウエーブレットタイプ，双直交スプライン・ウエーブレットタイプ、Ｄａｕｂｅｃｈｉｅｓウエーブレットタイプ、Ｍａｌｖａｒウエーブレットタイプ、Ｌｅｍａｒｉｅウエーブレットタイプ、Ｍｅｙｅｒウエーブレットタイプ、およびＳｙｍｌｅｔウエーブレットタイプなどの他のＣｏｉｆｌｅｔウエーブレット変換を使用することができる。
【００６７】
最適に選択されたウエーブレットは信号に関連する向きを表すウエーブレット係数に高いエントリを与え、確率的ノイズに関連する向きを表すウエーブレット係数に低いエントリを与える（図１のｅ））。これによって、信号の決定によりひき起こされる最大分散の向きへ向かって投射するための収縮条件を定めることができるようにされる。最後に、逆ウエーブレット変換が、挟雑物をなくされた時系列を再構成できるものである状態空間ベクトルを回復する。
【００６８】
決定論的ノイズが付加サブ空間を充たす。それは多様な信号から分離すべきである。ノイズに関連するサブ空間は、ノイズを付加基準センサに記録し、状態空間ベクトルをウエーブレット基準システムに変換することにより、場所を定められかつ記述される。その後で、源センサの時系列中の状態空間ベクトルの特徴が特定され、状態空間内での簡単な差し引きが実行される。この手順は一般の相互相関技術より優れている。その理由は決定論的ノイズのダイナミックな諸特性が考慮されているためである。この目的のためにはウエーブレット変換は決して用いられておらず、特に基準センサに結び付いては用いられなかったと信ぜられる。
【００６９】
説明しているノイズ減少法は磁気心電計から有用なデータを得るのに特に有用である。挟雑物をなくされた信号の特に有益な使用は妊娠している母親の、特に人の、体内の胎児の良好な状態の判定にある。妊娠のある段階中は胎児のＥＣＧを記録することは、胎児の絶縁脂肪層のために、非常に困難である。組織の透磁率は自由空間内のそれであるので、胎児のＭＣＧはその困難を免れる。しかし、これまでは、胎児の信号が非常に弱いことと、使用できないほど低い信号対ノイズ比のために、遮蔽の外部のＳＱＵＩＤ装置を使用しての胎児の心臓異常の存在を診断することは不可能であった。ここで説明している技術を用いると、母親から受けた信号を胎児のそれから分離することと、胎児の心拍異常を判定することが今や可能である。
【００７０】
開示されているＮＬＤ技術はＳＱＵＩＤ技術とともに大きな利点を提供もする。ＳＱＵＩＤ磁気心電計検査では遮蔽室は不要である。しかし遮蔽がないとノイズが増加して、ここで説明しているような強力なノイズ打ち消し技術を要する結果となる。
【００７１】
本発明の方法の重要な面の１つは適応しきい値化を使用することである。ここで使用しているように、しきい値化とはウエーブレット係数の固有スペクトラムを分割することである。
【００７２】
時系列を状態空間内に埋め込んだ後で、単一の状態空間ベクトル「ｘ」のおのおのに対して最も近いものの探索が行われ、ウエーブレット変換が適用される。その後で、「ｘ」に最も近いものである変換された全てのベクトルから平均を求 めることにより、質量中心ウエーブレットがつくられる。
【００７３】
ノイズと信号により占められているサブ空間を完全に分離することが可能であるならば、強力なしきい値化を実行できる。その場合にはノイズに属する全ての係数が零にセットされ、残りがそのまま保たれる。しかし、一般に、サブ空間は重なり合っているので、適応しきい値化を求められる。それはある係数が信号情報とノイズ情報を含んでいるという事実を説明する。
【００７４】
ソフトしきい値化では、ノイズ係数は零（ハード）にはセットされず、ある値、たとえば、平均値（ソフト）、にセットされる。これはそれら特定の係数のある情報を保持するがそれの重要性は減ずる。信号のサブ空間と重なり合うノイズが多いとそれらを分離することが一層困難であり、適応しきい値化が一層重要となる。
【００７５】
ＮＬＤの数学的方法論を支える概念は、決定論的信号構造を特定および記述するためにウエーブレット変換を用いる再構成された状態空間内の局部的な投射の性能である。目標は、それの源（それが心臓のノイズであれ、信号であれ）とは独立な任意の決定論的プロセスにより発生されたサブ空間の場所を定め、分離することである。この手順は、（１）確率的ノイズからのサブ空間の分離、（２）種々のサブ空間の分離、の２つの部分より成る。それらについては下で説明する。
【００７６】
観察されているシステムのダイナミックな諸特性を表現するためには、時間領域の代わりに再構成された状態空間内で動作すると有用である。図１のａ）は標本時系列のｘ成分の時間領域のグラフを示す。それは
【００７７】
【数１】
2005503855.tif 000003
【００７８】
として知られており次のように定義される。
【００７９】
ｘn+1＝１．４−ｘn2＋０．３ｙn
ｙn+1＝ｘn
明らかに、基礎を成す（決定論的）システムのいかなるダイナミック特性も認識することが困難である。
【００８０】
【数２】
2005503855.tif 000004
【００８１】
の時間遅延埋め込みによって次の状態空間ベクトルが導かれる。
【００８２】
【数３】
2005503855.tif 000005
【００８３】
ここにτは時間遅延を示し、ｍは埋め込みの次元を示す。
【００８４】
【数４】
2005503855.tif 000006
【００８５】
の状態空間は図１のｂ）に与えられている。ここで、τ＝１とｍ＝２を用いると状態空間ベクトルの成分が、成分ｘn+τをｘnに対して描くことにより二次元グラフで示されている。少なくとも近似的な決定論的システムが密に存在する軌道を導き、それは、確率的過程が状態空間全体にわたって無作為な分布を行わせるようなサブ空間に制約される。
【００８６】
決定論的信号に白色ノイズを重ねると密に存在する主軌道が歪まされる。図１のｃ）の左のグラフは図１のｂ）のいくつかの（集められた）軌道の抜粋したものを示す。ノイズをこの抜粋に重ねた効果が図１のｃ）の右の部分に示されている。
【００８７】
次のステップは状態空間内で決定論的構造を特定しかつ記述することである。この目的のためには、状態空間ベクトルを適切な基準システムに変換すると有用である。「適切」というのは決定論的構造に最も良く適応する基準関数を見出すことを試みることを意味する。この場合には新しい基準システムの領域内で僅かな係数によって決定論を記述することが可能である。これは、向き情報を圧縮できるという事実による。対照的に、確率論的ノイズは圧縮できず、したがって、基礎係数の完全なセットを再生する必要がある。
【００８８】
ウエーブレット変換はウエーブレットと呼ばれる多数の極めて適応的な基準関数を与える。それは信号ｘ（ｔ）でたたみこまれた基準関数（ウエーブレット）の翻訳と拡張により定められる。付加換算係数（ほとんどは２の冪）がウエーブレット変換を数学的顕微鏡のように機能させる；それは種々の尺度にある１つの観察信号情報をそれの場所に依存させることを意味する。正確にはこの特性は有用である。その理由はウエーブレット領域内で状態空間ベクトルの多重分解度表現を得るからである（図１のｄ））。
【００８９】
最適に選択されたウエーブレットで、信号に関連する向きを表すウエーブレット係数内の高いエントリと、確率的ノイズに関連する向きを表すウエーブレット係数内の低いエントリを期待できる（図１のｅ））。これによって、信号の決定論によりひき起こされる最大分散の向きへ向かって投射するための収縮条件を定めることができるようにされる。最後に、逆ウエーブレット変換が、挟雑物をなくされた時系列を再構成できるものである状態空間ベクトルを回復する。
【００９０】
ウエーブレット係数の適応（ハードまたはソフト）しきい値化が、状態空間内 においても信号回復に良く適しており、かつＭＣＧ時系列信号またはＥＣＧ時系列信号のノイズ除去において重要である。
【００９１】
決定論的ノイズが付加サブ空間を充たす。それは多様な信号から分離すべきである。ＭＣＧへの適用においては、ノイズに関連するサブ空間は、ノイズを付加基準センサに記録し、かつ状態空間ベクトルをウエーブレット基準システムに変換することにより、場所を定められかつ記述される。その後で、源センサの時系列中の状態空間ベクトルの特徴が特定され、状態空間内での簡単な差し引きが実行される。この手順は一般の相互相関技術より優れている。その理由は決定論的ノイズのダイナミックな諸特性が考慮されているためである。
【００９２】
ＮＬＤの重要性は、信号の徹底的な平均化と濾波（両方とも測定すべき信号を歪ませる）に頼ることなく、何桁も強いノイズから弱い有用な生体電気や生体磁気を分離できる潜在的なそれの性能にある。
【００９３】
この新規なノイズ除去法の効果を示すために、主センサにより記録された純粋な信号として２００Ｈｚで記録された、健康な患者の心電計データを用いて模擬された信号にそれを適用した。
例  １
健康な心臓の５秒間ＥＣＧ記録から出発して、模擬された信号にＮＬＤが適用され、２００Ｈｚの帯域幅で記録され、主センサから純粋な信号として取られた。このＥＣＧは５０Ｈｚノッチフィルタと１００Ｈｚの第２種低域フィルタにより予め濾波された（図２ａと図２ｂ）。
【００９４】
それに続いて、心電計の分散と呼ばれている３０％の振幅の分散と、決定論的ノイズが白色ノイズに加えられる。決定論的ノイズの周波数ピークは１６  ２／３Ｈｚ、５０Ｈｚ（ヨーロッパにおけるレール電源および低調波）、および６０Ｈｚ（信号解析システム）で、振幅分散は１００％であった。
【００９５】
決定論的ノイズの周波数ピークは１６  ２／３Ｈｚ、５０Ｈｚ（ヨーロッパにおける電源および低調波）、および６０Ｈｚ（信号解析システム）で、振幅分散は１００％であった（図２ｃおよび図２ｄ参照）。上記パラメータと同じパラメータを用いるが、更に、決定論的ノイズ成分に対して振幅の分散と一定の位相推移がある、基準ノイズ時系列がつくられた。
【００９６】
図２ｃは白色ノイズが重畳されている信号を示し、図２ｄは結果としてのノイズスペクトラムを示す。基準時系列は、上記パラメータと同じパラメータを用いるが、更に、決定論的ノイズ成分に対して振幅の分散と一定の位相推移があるノイズをつくることにより発生される。
【００９７】
ウエーブレット変換と状態空間内での差し引きの後で、図２ｅは挟雑物のない時系列を示す。基準時系列は、上記パラメータと同じパラメータを用いるが、更に、決定論的ノイズ成分に対して振幅の分散と一定の位相推移があるノイズをつくることにより発生される。図２ｆは、ウエーブレット変換と状態空間内での差し引きの後の心電図の周波数スペクトラムを示す。
【００９８】
心拍（ノイズ除去の質の良い指標）の間の基線にはノイズがほとんど無いことを推論できる。したがって、本発明の方法は信号保持とかなりのノイズ減少を行った。
例  ２
測定される信号データの例として、研究用ＨＴＳＱＵＩＤ装置を用いて遮蔽室の外部で記録された、健康な患者の磁気心電計から得られたデータが図３に示されている。図３ａに示されているように、［５］秒の磁気心電計信号が得られた。患者の心拍は図３ａではかろうじて見えるだけである。同時に記録されたノイズ時系列が図３ｃに示されているように記録された。
【００９９】
図３ａに示されている信号の周波数スペクトラムが図３ｂに示されている。図３ｃに示されている信号の周波数スペクトラムが図３ｄに示されている。スペクトラム中の５０Ｈｚピークの幅のためにノッチフィルタは使用しなかった。
【０１００】
この測定のために、基線が７ｃｍである第１種の２つの軸線方向グラジオメーターを７ｃｍの高さを隔てて相互に上下に配置した。この例では、上のグラジオメーターが基準信号を記録した（図３ｃと図３ｄ）。
【０１０１】
図３ｅと図３ｆは時系列を、このノイズ除去法の結果として得られたそれの対応するフーリェスペクトラムと共に示す。図３ｅの再構成された磁気心電図（ＭＣＧ）では、心拍のかなりの詳細が明らかにされている。また、心拍の間の基線にはノイズがほとんど無い。
例  ３
図４ａ〜４ｂは先行技術のノイズ除去技術の１つに対する本発明の装置および方法の優秀性を示す。
【０１０２】
この方法を基にした例２のデータセットの解析が図４ａと図４ｂに示されている。「非線形時系列解析（Ｎｏｎｌｉｎｅａｒ  Ｔｉｍｅ  Ｓｅａｒｉｅｓ  Ａｎａｌｙｓｉｓ（ＴＩＳＥＡＮ），）という標題のＨｅｇｇｅｒ他の論文、参照することによりここに含まれる、に記述されているツール「ｇｈｋｓｓ」がそのデータセットを解析するために用いられる。これは「心電図用非線形ノイズ減少（Ｎｏｎｌｉｎｅａｒ  Ｎｏｉｓｅ  Ｒｅｄｕｃｔｉｏｎ  Ｆｏｒ  Ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｍｓ）」（Ｃｈａｏｓ６：８７，１９９５）のアルゴリズム形式である。
【０１０３】
ツール「ｇｈｋｓｓ」はそのデータセットに適用され、図４ｃと図４ｄに示されている結果を得た。明らかに、ＮＬＤはこの場合にはより良好なノイズ減少質に到達し、それはそれぞれの残留物により明らかにされている（図４ｅと図４ｆ参照）。これは、「ｇｈｋｓｓ」が状態空間内で重なり合っているサブ空間を分離できないという事実による。これはＮＬＤの最も重要な特徴の１つである。
【０１０４】
これを示すために、同じデータの解析が、「心電図用非直線ノイズ減少（Ｎｏｎｌｉｎｅａｒ  Ｎｏｉｓｅ  Ｒｅｄｕｃｔｉｏｎ  Ｆｏｒ  Ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｍｓ）」（Ｃｈａｏｓ６：８７，１９９５）という標題のＳｃｈｒｅｉｂｅｒ他の論文、参照することによりここに含まれる、に記述されている技術を基にして実行される。要約した形式で、この方法は共分散マトリックスを用いて状態空間内の局部的投射を行うことにより確率的ノイズを減少する。信号の特徴の場所が状態空間内で定められて、共分散マトリックスの最大固有値により示されているノイズなしサブ空間上に投射される。この方法は良く機能するが、状態空間内の信号のサブ空間の次元が既知である場合に限られる。
【０１０５】
解析の結果が図４ｃと図４ｄに示されている。ＮＬＤはこの場合により良好な質のノイズ減少に到達する。これはそれぞれの残留ノイズによっても示されている（図４ｅと図４ｆ）。ＮＬＤの残留物は「ｇｈｋｓｓ」のそれよりはるかに少ない。その理由は、ＮＬＤとは対照的に、「ｇｈｋｓｓ」は状態空間内の重なり合っているサブ空間を分離できないが、ＮＬＤはできる、という事実による。そ の分離性能はＮＬＤの最も重要な特徴である。
【０１０６】
ＮＬＤを他の既存の技術とも比較した。周波数依存グラジオメトリー（ＦＤＧ）とＮＬＤを同じＭＣＧ資料に適用したが、ＮＬＤははるかに優れたノイズ除去を行ったことが判明した。
例  ４
この例はＭＣＧを用いた胎児の心拍の測定へのこの方法の適用可能性を示す。図５ａは、遮蔽内でのＬＴＳＱＵＩＤによる妊娠中の女性から記録された時系列の３秒間の抜粋を示す。図５ｂでは５０Ｈｚにおける典型的なノイズピークが失われている。これは遮蔽室の使用を示す。最初のＮＬＤステップでは決定論的ノイズ成分が除去される。
【０１０７】
図５ｃは第２のＮＬＤステップ後の結果を示す。母親のＭＣＧは人の低頻度（呼吸）動作で依然として汚染されていることを見ることができる。それは、観察時間を長くすることにより除去できる。図５ｄにおけるそれのパワースペクトラムにはノイズのピークはなく、白色ノイズのレベルが低下させられたことを示す。
【０１０８】
母親のＭＣＧを時系列から除去し、ＮＬＤを再び適用すると、胎児のＭＣＧのＱＲＳ複合物が図５ｅに示されているように得られる。以前に重ね合わされていた心拍でも分離されることに注目されたい。図５ｆに示されている胎児のＭＣＧのスペクトルエネルギーははるかに小さいが、図５ｄに示されている母親のそれのスペクトルエネルギーと同じ帯域幅内に含まれている。これは極めて適応性に富むノイズ除去法の重要性を更に示すものである。
【０１０９】
この発明のプログラミングはデジタル信号処理技術の当業者によって実現できる。
【０１１０】
上記例はウエーブレット変換を用いての磁気心電図の非線形ノイズ除去法の効果的な実現を示すものである。種々の好適な実施形態を示しかつ説明したが、そのような開示により本発明を限定する意図はなく、むしろ、添付した特許請求の範囲で規定されているように、本発明の要旨および範囲内に入るあらゆる変更および代替構造を含めることを意図するものである。
【図面の簡単な説明】
【０１１１】
【図１】ａ）時間的に見た観察されている系を表す。ｂ）再構成された状態空間に関して見た観察されている系を表し、少なくとも近似的に決定論的な系の密に存在する軌道を示す。ｃ）信号にノイズが混じる前および混じった後の図１ｂの状態空間の一部を示す。ｄ）ウエーブレット領域内の状態空間ベクトルの多重分解表現である。ｅ）信号に関連する向きを表すウエーブレット係数における高エントリと確率的ノイズに関連する向きを表すウエーブレット係数における低エントリを示す。
【図２ａ】主センサにより記録された純粋な信号として２００Ｈｚに記録されている５秒間の心電計データを表す。
【図２ｂ】５０Ｈｚノッチフィルタおよび１００Ｈｚ第２種低域フィルタによる事前濾波後の心電計の周波数スペクトラムを示す。
【図２ｃ】白色ノイズが重畳されている図２ｂの信号を示す。
【図２ｄ】図２ｃの結果としてのノイズスペクトラムを示す。
【図２ｅ】図２ｃの信号の状態空間内でのウエーブレット変換および差し引き後の、挟雑物のない時系列を示す。
【図２ｆ】状態空間内でのウエーブレット変換および差し引き後の心電計周波数スペクトラムを示す。
【図３ａ】遮蔽室外で記録された、心臓信号の主成分（Ｒ波）のみが見える、５秒間の磁気心電計信号を表す。
【図３ｂ】図３ａの信号の周波数スペクトラムを示す。
【図３ｃ】図３ａの同時に記録された５秒間のノイズ信号を表す。
【図３ｄ】図３ｃに示されている信号のフーリェスペクトラムを表す。
【図３ｅ】このノイズ除去法の結果得られた時系列を示す。
【図３ｆ】図３ｅの時系列に対応するフーリェスペクトラムを示す。
【図４ａ】例１のデータを用いる元の時系列を示す。
【図４ｂ】図４ａの信号の周波数スペクトラムを示す。
【図４ｃ】ｇｈｋｓｓによるノイズ減少後の時系列を示す。
【図４ｄ】図４ｃの信号のパワースペクトラムを示す。
【図４ｅ】このノイズ除去法を用いる図４ａの信号中の残留ノイズを示す。
【図４ｆ】「ｇｈｋｓｓ」でのノイズ減少後の図４ａの信号中の残留ノイズを示す。
【図５ａ】遮蔽内での低温ＳＱＵＩＤによる妊娠中の女性から記録された３秒間の時系列の抜粋を示す。
【図５ｂ】遮蔽室の使用を示す、５０Ｈｚにおける失われている典型的なノイズピークのいくらかを示す。
【図５ｃ】観察時間を長くすることにより除去できる、見えるが、人の低頻度（呼吸）動作で汚染されている母親のＭＣＧを示す、ＮＬＤ適用後の結果を示す。
【図５ｄ】白色ノイズレベルが低下していることを示す、ノイズピーク無しのパワースペクトラムを示す。
【図５ｅ】時系列からの母親のＭＣＧの除去とＮＬＤの再度の適用後の胎児のＭＣＧのＱＲＳ複合物のスペクトラムを示して、以前に重ね合わされた心拍が分離されていることを示す。
【図５ａ】母親のＭＣＧのスペクトルエネルギーを示す。
【図５ｆ】母親のＭＣＧのスペクトルエネルギー（図５ｄ）よりもはるかに低いが、母親のＭＣＧのスペクトルエネルギー（図５ｄ）と同じ帯域幅内にあり、非常に適応性に富むノイズ除去法の重要性を示す、胎児のＭＣＧのスペクトルエネルギーを示す。
【図１】
2005503855.tif 000007
【図２ａ】
2005503855.tif 000008
【図２ｂ】
2005503855.tif 000009
【図２ｃ】
2005503855.tif 000010
【図２ｄ】
2005503855.tif 000011
【図２ｅ】
2005503855.tif 000012
【図２ｆ】
2005503855.tif 000013
【図３ａ】
2005503855.tif 000014
【図３ｂ】
2005503855.tif 000015
【図３ｃ】
2005503855.tif 000016
【図３ｄ】
2005503855.tif 000017
【図３ｅ】
2005503855.tif 000018
【図３ｆ】
2005503855.tif 000019
【図４ａ】
2005503855.tif 000020
【図４ｂ】
2005503855.tif 000021
【図４ｃ】
2005503855.tif 000022
【図４ｄ】
2005503855.tif 000023
【図４ｅ】
2005503855.tif 000024
【図４ｆ】
2005503855.tif 000025
【図５ａ】
2005503855.tif 000026
【図５ｂ】
2005503855.tif 000027
【図５ｃ】
2005503855.tif 000028
【図５ｄ】
2005503855.tif 000029
【図５ｅ】
2005503855.tif 000030
【図５ｆ】
2005503855.tif 000031
【国際公開パンフレット】
2005503855.tif 000032
2005503855.tif 000033
2005503855.tif 000034
2005503855.tif 000035
2005503855.tif 000036
2005503855.tif 000037
2005503855.tif 000038
2005503855.tif 000039
2005503855.tif 000040
2005503855.tif 000041
2005503855.tif 000042
2005503855.tif 000043
2005503855.tif 000044
2005503855.tif 000045
2005503855.tif 000046
2005503855.tif 000047
2005503855.tif 000048
2005503855.tif 000049
2005503855.tif 000050
2005503855.tif 000051
2005503855.tif 000052
2005503855.tif 000053
2005503855.tif 000054
2005503855.tif 000055
2005503855.tif 000056
2005503855.tif 000057
2005503855.tif 000058
2005503855.tif 000059
2005503855.tif 000060
2005503855.tif 000061
【国際公開パンフレット（コレクトバージョン）】
2005503855.tif 000062
2005503855.tif 000063
【国際調査報告】
2005503855.tif 000064