(57)【要約】【課題】電流素片数が不自然に制限されず、あり得ない局所解に陥ったりせずに、脳磁気検出装置により測定した磁場データから電流素片分布を推定する。【解決手段】空間分布した電流素片を発生するステップと、電流素片分布から計算される計算値と測定値の誤差が十分小さいように電流素片のパラメータを変更するステップと、パラメータを変更した電流素片の各々がいずれかのカバリング・ボールに含まれるようなカバリング・ボールを発生するステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくするステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更するステップの一部または全部を繰り返して電流素片を局在化させる。【選択図】図８
2005058558.tif 000002
【特許請求の範囲】
【請求項１】
３次元空間に分布したＡ個の電流素片を発生する電流素片発生ステップと、電流素片の分布から計算される計算値と測定値の誤差が十分小さいように電流素片のパラメータを変更する第１探索ステップと、前記第１探索ステップでパラメータを変更したＡ個の電流素片の各々がいずれかのカバリング・ボールに含まれるようなＢ（＜Ａ）個の半径εのカバリング・ボールを発生するカバリング・ボール発生ステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくする半径縮小ステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更する第２探索ステップとを有し、前記ステップの一部または全部を繰り返して電流素片をＢ箇所に局在化させることを特徴とする電流素片局在化処理方法。
【請求項２】
３次元空間に分布したＡ個の電流素片の各々がいずれかのカバリング・ボールに含まれるようなＢ（＜Ａ）個の半径εのカバリング・ボールを発生するカバリング・ボール発生ステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくする半径縮小ステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更する探索ステップとを有し、前記ステップの一部または全部を繰り返して電流素片をＢ箇所に局在化させることを特徴とする電流素片局在化処理方法。
【請求項３】
請求項１または請求項２に記載の電流素片局在化処理方法において、前記測定値が、脳磁気検出装置により測定した磁場データであることを特徴とする電流素片局在化処理方法。
【請求項４】
３次元空間に分布したＡ個の電流素片を発生する電流素片発生手段と、電流素片の分布から計算される計算値と測定値の誤差が十分小さいように電流素片のパラメータを変更する第１探索手段と、前記第１探索手段でパラメータを変更したＡ個の電流素片の各々がいずれかのカバリング・ボールに含まれるようなＢ（＜Ａ）個の半径εのカバリング・ボールを発生するカバリング・ボール発生手段と、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくする半径縮小手段と、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更する第２探索手段とを具備し、前記手段の一部または全部を繰り返し作動させて電流素片をＢ箇所に局在化させることを特徴とする電流素片局在化処理装置。
【請求項５】
３次元空間に分布したＡ個の電流素片の各々がいずれかのカバリング・ボールに含まれるようなＢ（＜Ａ）個の半径εのカバリング・ボールを発生するカバリング・ボール発生手段と、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくする半径縮小手段と、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更する探索手段とを具備し、前記手段の一部または全部を繰り返し作動させて電流素片をＢ箇所に局在化させることを特徴とする電流素片局在化処理装置。
【請求項６】
請求項４または請求項５に記載の電流素片局在化処理装置において、前記測定値が、脳磁気検出装置により測定した磁場データであることを特徴とする電流素片局在化処理装置。
【請求項７】
空間に分布したＡ個の信号源を発生する信号源発生ステップと、信号源の分布から計算される計算値と目標値の誤差が十分小さいように信号源の位置を動かす第１探索ステップと、前記第１探索ステップで動かしたＡ個の信号源の各々がいずれかのカバリング・セルに含まれるようなＢ（＜Ａ）個の所定サイズのカバリング・セルを発生するカバリング・セル発生ステップと、各カバリング・セルが含む信号源が該カバリング・セルから飛び出さない範囲で各カバリング・セルのサイズを小さくするサイズ縮小ステップと、各カバリング・セルが含む信号源が該カバリング・セルから出ないという制約下で信号源の位置を動かす第２探索ステップとを有し、前記ステップの一部または全部を繰り返して信号源をＢ箇所に局在化させることを特徴とする信号源局在化処理方法。
【請求項８】
空間に分布したＡ個の信号源の各々がいずれかのカバリング・セルに含まれるようなＢ（＜Ａ）個の所定サイズのカバリング・セルを発生するカバリング・セル発生ステップと、各カバリング・セルが含む信号源が該カバリング・セルから飛び出さない範囲で各カバリング・セルのサイズを小さくするサイズ縮小ステップと、各カバリング・セルが含む信号源が該カバリング・セルから出ないという制約下で信号源の位置を動かす探索ステップとを有し、前記ステップの一部または全部を繰り返して信号源をＢ箇所に局在化させることを特徴とする信号源局在化処理方法。
【請求項９】
空間に分布したＡ個の信号源を発生する信号源発生手段と、信号源の分布から計算される計算値と目標値の誤差が十分小さいように信号源の位置を動かす第１探索手段と、前記第１探索手段で動かしたＡ個の信号源の各々がいずれかのカバリング・セルに含まれるようなＢ（＜Ａ）個の所定サイズのカバリング・セルを発生するカバリング・セル発生手段と、各カバリング・セルが含む信号源が該カバリング・セルから飛び出さない範囲で各カバリング・セルのサイズを小さくするサイズ縮小手段と、各カバリング・セルが含む信号源が該カバリング・セルから出ないという制約下で信号源の位置を動かす第２探索手段とを具備し、前記手段の一部または全部を繰り返し作動させて信号源をＢ箇所に局在化させることを特徴とする信号源局在化処理装置。
【請求項１０】
空間に分布したＡ個の信号源の各々がいずれかのカバリング・セルに含まれるようなＢ（＜Ａ）個の所定サイズのカバリング・セルを発生するカバリング・セル発生手段と、各カバリング・セルが含む信号源が該カバリング・セルから飛び出さない範囲で各カバリング・セルのサイズを小さくするサイズ縮小手段と、各カバリング・セルが含む信号源が該カバリング・セルから出ないという制約下で信号源の位置を動かす探索手段とを有し、前記手段の一部または全部を繰り返し作動させて信号源をＢ箇所に局在化させることを特徴とする信号源局在化処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電流素片局在化処理方法及び装置並びに信号源局在化処理方法及び装置に関し、さらに詳しくは、空間に分布したＡ個の電流素片または信号源をＣ（＜Ａ）箇所に局在化させるための方法および装置に関する。
【背景技術】
【０００２】
従来、脳磁気検出装置により測定した磁場データから、脳に存在すると仮定した電流素片（等価電流双極子）の分布を推定する方法および装置が種々提案されている（例えば、非特許文献１，非特許文献２参照。）。
【非特許文献１】原宏・栗城真也共著「脳磁気科学、−ＳＱＵＩＤ計測と医学応用−」オーム社出版局、平成９年１月２５日、ｐ．１０１−１２０
【非特許文献２】電子情報通信学会編、武田常広著「電子情報通信レクチャーシリーズＤ−２４、脳工学」コロナ社、２００３年４月１１日、ｐ．１２８−１４０
【発明の開示】
【発明が解決しようとする課題】
【０００３】
従来の方法および装置では、電流素片の数が計算の都合上から不自然に制限されたり、電流素片の数を増やせても実際にはありそうもない局所解（例えば存在しそうもない場所に電流素片が存在しているような解）に陥ったりする等、満足な推定結果が得られない問題点があった。
そこで、本発明の目的は、電流素片または信号源の数が不自然に制限されず、また、実際にはありそうもない局所解に陥ったりせずに、電流素片または信号源の分布を推定することが出来る電流素片局在化処理方法及び装置並びに信号源局在化処理方法及び装置を提供することにある。
【課題を解決するための手段】
【０００４】
第１の観点では、本発明は、３次元空間に分布したＡ個の電流素片を発生する電流素片発生ステップと、電流素片の分布から計算される計算値と測定値の誤差が十分小さいように電流素片のパラメータを変更する第１探索ステップと、前記第１探索ステップでパラメータを変更したＡ個の電流素片の各々がいずれかのカバリング・ボールに含まれるようなＢ（＜Ａ）個の半径εのカバリング・ボールを発生するカバリング・ボール発生ステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくする半径縮小ステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更する第２探索ステップとを有し、前記ステップの一部または全部を繰り返して電流素片をＢ箇所に局在化させることを特徴とする電流素片局在化処理方法を提供する。
上記第１の観点による電流素片局在化処理方法では、発生させる電流素片の数Ａが不自然に制限されることはない。また、カバリング・ボールを用いることで、Ａ個の電流素片をＢ個にグループ化すると共に各グループ内の電流素片を１カ所に集めることが出来るようなグループ化パターンを探索する。この結果、実際にありそうなＢ箇所に電流素片を局在化させることが出来る。なお、１カ所に集中した電流素片を１つの電流素片に合成してもよい。
【０００５】
第２の観点では、本発明は、３次元空間に分布したＡ個の電流素片の各々がいずれかのカバリング・ボールに含まれるようなＢ（＜Ａ）個の半径εのカバリング・ボールを発生するカバリング・ボール発生ステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくする半径縮小ステップと、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更する探索ステップとを有し、前記ステップの一部または全部を繰り返して電流素片をＢ箇所に局在化させることを特徴とする電流素片局在化処理方法を提供する。
上記第２の観点による電流素片局在化処理方法では、電流素片の数Ａが不自然に制限されることはない。また、カバリング・ボールを用いることで、Ａ個の電流素片をＢ個にグループ化すると共に一つのグループ内の電流素片を１カ所に集めることが出来るようなグループ化パターンを探索する。この結果、実際にありそうなＢ箇所に電流素片を局在化させることが可能になる。なお、１カ所に集中した電流素片を１つの電流素片に合成してもよい。
【０００６】
第３の観点では、本発明は、上記構成の電流素片局在化処理方法において、前記測定値が、脳磁気検出装置により測定した磁場データであることを特徴とする電流素片局在化処理方法を提供する。
上記第３の観点による電流素片局在化処理方法では、脳磁気検出装置により測定した磁場データから、脳に存在すると仮定した電流素片の分布を推定することが出来る。
【０００７】
第４の観点では、本発明は、３次元空間に分布したＡ個の電流素片を発生する電流素片発生手段と、電流素片の分布から計算される計算値と測定値の誤差が十分小さいように電流素片のパラメータを変更する第１探索手段と、前記第１探索手段でパラメータを変更したＡ個の電流素片の各々がいずれかのカバリング・ボールに含まれるようなＢ（＜Ａ）個の半径εのカバリング・ボールを発生するカバリング・ボール発生手段と、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくする半径縮小手段と、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更する第２探索手段とを具備し、前記手段の一部または全部を繰り返し作動させて電流素片をＢ箇所に局在化させることを特徴とする電流素片局在化処理装置を提供する。
上記第４の観点による電流素片局在化処理装置では、上記第１の観点による電流素片局在化処理方法を好適に実施できる。
【０００８】
第５の観点では、本発明は、３次元空間に分布したＡ個の電流素片の各々がいずれかのカバリング・ボールに含まれるようなＢ（＜Ａ）個の半径εのカバリング・ボールを発生するカバリング・ボール発生手段と、各カバリング・ボールが含む電流素片が該カバリング・ボールから飛び出さない範囲で各カバリング・ボールの半径を小さくする半径縮小手段と、各カバリング・ボールが含む電流素片が該カバリング・ボールから出ないという制約下で電流素片のパラメータを変更する探索手段とを具備し、前記手段の一部または全部を繰り返し作動させて電流素片をＢ箇所に局在化させることを特徴とする電流素片局在化処理装置を提供する。
上記第５の観点による電流素片局在化処理装置では、上記第２の観点による電流素片局在化処理方法を好適に実施できる。
【０００９】
第６の観点では、本発明は、上記第４または第５の観点による電流素片局在化処理装置において、前記測定値が、脳磁気検出装置により測定した磁場データであることを特徴とする電流素片局在化処理装置を提供する。
上記第６の観点による電流素片局在化処理装置では、上記第３の観点による電流素片局在化処理方法を好適に実施できる。
【００１０】
第７の観点では、本発明は、空間に分布したＡ個の信号源を発生する信号源発生ステップと、信号源の分布から計算される計算値と目標値の誤差が十分小さいように信号源の位置を動かす第１探索ステップと、前記第１探索ステップで動かしたＡ個の信号源の各々がいずれかのカバリング・セルに含まれるようなＢ（＜Ａ）個の所定サイズのカバリング・セルを発生するカバリング・セル発生ステップと、各カバリング・セルが含む信号源が該カバリング・セルから飛び出さない範囲で各カバリング・セルのサイズを小さくするサイズ縮小ステップと、各カバリング・セルが含む信号源が該カバリング・セルから出ないという制約下で信号源の位置を動かす第２探索ステップとを有し、前記ステップの一部または全部を繰り返して信号源をＢ箇所に局在化させることを特徴とする信号源局在化処理方法を提供する。
上記構成において、信号源としては、電流素片や磁極のような磁気発生源，音発生源，エコー発生源，熱発生源，振動発生源が挙げられる。また、カバリング・セルとしては、２次元空間ではカバリング・サークル，３次元空間ではカバリング・ボールが挙げられる。
上記第７の観点による信号源局在化処理方法では、発生させる信号源の数Ａが不自然に制限されることはない。また、カバリング・セルを用いることで、Ａ個の信号源をＢ個にグループ化すると共に各グループ内の信号源を１カ所に集めることが出来るようなグループ化パターンを探索する。この結果、実際にありそうなＢ箇所に信号源を局在化させることが出来る。なお、１カ所に集中した信号源を１つの信号源に合成してもよい。
【００１１】
第８の観点では、本発明は、空間に分布したＡ個の信号源の各々がいずれかのカバリング・セルに含まれるようなＢ（＜Ａ）個の所定サイズのカバリング・セルを発生するカバリング・セル発生ステップと、各カバリング・セルが含む信号源が該カバリング・セルから飛び出さない範囲で各カバリング・セルのサイズを小さくするサイズ縮小ステップと、各カバリング・セルが含む信号源が該カバリング・セルから出ないという制約下で信号源の位置を動かす探索ステップとを具備し、前記ステップの一部または全部を繰り返して信号源をＢ箇所に局在化させることを特徴とする信号源局在化処理方法を提供する。
上記第８の観点による信号源局在化処理方法では、信号源の数Ａが不自然に制限されることはない。また、カバリング・セルを用いることで、Ａ個の信号源をＢ個にグループ化すると共に一つのグループ内の信号源を１カ所に集めることが出来るようなグループ化パターンを探索する。この結果、実際にありそうなＢ箇所に信号源を局在化させることが可能になる。なお、１カ所に集中した信号源を１つの信号源に合成してもよい。
【００１２】
第９の観点では、本発明は、空間に分布したＡ個の信号源を発生する信号源発生手段と、信号源の分布から計算される計算値と目標値の誤差が十分小さいように信号源の位置を動かす第１探索手段と、前記第１探索手段で動かしたＡ個の信号源の各々がいずれかのカバリング・セルに含まれるようなＢ（＜Ａ）個の所定サイズのカバリング・セルを発生するカバリング・セル発生手段と、各カバリング・セルが含む信号源が該カバリング・セルから飛び出さない範囲で各カバリング・セルのサイズを小さくするサイズ縮小手段と、各カバリング・セルが含む信号源が該カバリング・セルから出ないという制約下で信号源の位置を動かす第２探索手段とを具備し、前記手段の一部または全部を繰り返し作動させて信号源をＢ箇所に局在化させることを特徴とする信号源局在化処理装置を提供する。
上記第９の観点による信号源局在化処理装置では、上記第７の観点による信号源局在化処理方法を好適に実施できる。
【００１３】
第１０の観点では、本発明は、空間に分布したＡ個の信号源の各々がいずれかのカバリング・セルに含まれるようなＢ（＜Ａ）個の所定サイズのカバリング・セルを発生するカバリング・セル発生手段と、各カバリング・セルが含む信号源が該カバリング・セルから飛び出さない範囲で各カバリング・セルのサイズを小さくするサイズ縮小手段と、各カバリング・セルが含む信号源が該カバリング・セルから出ないという制約下で信号源の位置を動かす探索手段とを具備し、前記手段の一部または全部を繰り返し作動させて信号源をＢ箇所に局在化させることを特徴とする信号源局在化処理装置を提供する。
上記第１０の観点による信号源局在化処理装置では、上記第８の観点による信号源局在化処理方法を好適に実施できる。
【発明の効果】
【００１４】
本発明の電流素片局在化処理方法及び装置並びに信号源局在化処理方法及び装置によれば、電流素片または信号源の数が不自然に制限されず、また、実際にはありそうもない局所解に陥ったりせずに、電流素片または信号源の分布を推定することが出来る。
【発明を実施するための最良の形態】
【００１５】
以下、図に示す実施の形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
【実施例１】
【００１６】
図１は、本発明の実施例１に係る電流素片推定装置１００を示す構成図である。
この電流素片推定装置１００は、被検体Ｈに刺激を与えた時に脳から発生する磁場データを脳磁気検出装置Ｎを介して収集するデータ収集部１と、前記磁場データから脳に存在すると仮定した電流素片の分布を推定する電流素片推定処理部２と、推定した電流素片を視覚化して表示する表示部３とを具備している。
【００１７】
図２は、電流素片推定装置１００による電流素片推定処理の手順を示すフロー図である。
ステップＳ１では、データ収集部１は、図４に示すように、多数の磁気センサの３次元空間位置（図４の黒丸）とそれら磁気センサで測定した磁場ベクトル（図４の矢印）とを対応させた磁場データを得る。
【００１８】
ステップＳ２では、電流素片推定処理部２は、図５に示すように、脳が存在する３次元空間内にランダムに分布させてＡ個（例えば、Ａ＝３０〜４０）の電流素片Ｊを発生する。各電流素片Ｊのモーメントについては、電流方向はランダムに、電流値はありえそうな範囲内でランダムに定める。
【００１９】
ステップＳ３では、電流素片推定処理部２は、図６に示すように、多数の磁気センサの３次元空間位置（図６の黒丸）にＡ個の電流素片が発生させる磁気ベクトル（図６の矢印）を算出し、３次元空間位置（図６の黒丸）と磁気ベクトル（図６の矢印）とを対応させた計算磁場を得る。
【００２０】
ステップＳ４では、電流素片推定処理部２は、公知の評価関数を用いて磁場データ（図４）と計算磁場（図６）の誤差を計算し、誤差が十分小さくないならステップＳ５へ進み、誤差が十分小さいなら図３のステップＳ１１へ進む。
【００２１】
ステップＳ５では、電流素片推定処理部２は、誤差が小さくなる見込みがあるか否かを判定し（例えば、ステップＳ５を通った回数が１万回未満なら誤差が小さくなる見込みがあると判定し、１万回以上なら誤差が小さくなる見込みがないと判定する）、見込みがあるならステップＳ６へ進み、見込みがないならステップＳ２に戻って電流素片の発生からやり直す。
【００２２】
ステップＳ６では、電流素片推定処理部２は、ランダムに一つの電流素片Ｊを選択する。
ステップＳ７では、電流素片推定処理部２は、ランダムに一つのパラメータすなわち電流素片の位置またはモーメントを選択する。
ステップＳ８では、電流素片推定処理部２は、選択した電流素片Ｊの、選択したパラメータを、公知の最適化方法を用いて、誤差を小さくする方向に、少し変更する。そして、ステップＳ３に戻る。
【００２３】
上記ステップＳ２〜Ｓ８を十分な回数繰り返すことにより、誤差が十分小さくなる電流素片Ｊの分布と各電流素片のモーメントが定まるが、これは未だ実際にはありそうもない局所解（例えば存在しそうもない場所に電流素片Ｊが存在しているような解）である可能性が高い。そこで、図３のステップＳ１１以下を実行し、実際にありそうな解へと導く。
【００２４】
図３のステップＳ１１では、電流素片推定処理部２は、Ｂ個（例えば、Ｂ＝３〜１０）の半径εのカバリング・ボールをランダムに発生することを、各電流素片Ｊがいずれかのカバリング・ボールに入っている状態になるまで繰り返し、図７に示すように、各電流素片Ｊがいずれかのカバリング・ボールＶに含まれるような、Ｂ個の、半径εのカバリング・ボールＶを発生する。
【００２５】
ステップＳ１２では、電流素片推定処理部２は、図８に示すように、各カバリング・ボールＶが含む電流素片Ｊが該カバリング・ボールＶから飛び出さない範囲で各カバリング・ボールＶの半径を小さくする（カバリング・ボールＶの中心位置は変えない）。
【００２６】
ステップＳ１３では、電流素片推定処理部２は、全てのカバリング・ボールＶの半径が十分小さくないならステップＳ１４へ進み、十分小さいならステップＳ２１へ進む。
【００２７】
ステップＳ１４では、電流素片推定処理部２は、カバリング・ボールＶの半径が小さくなる見込みがあるか否かを判定し（例えば、ステップＳ１４を通った回数が１万回未満なら半径が小さくなる見込みがあると判定し、１万回以上なら半径が小さくなる見込みがないと判定する）、見込みがあるならステップＳ１５へ進み、見込みがないなら図２のステップＳ２に戻って電流素片の発生からやり直す。
【００２８】
ステップＳ１５では、電流素片推定処理部２は、ランダムに一つのカバリング・ボールＶを選択する。
ステップＳ１６では、電流素片推定処理部２は、選択したカバリング・ボールＶ中の、一つの電流素片Ｊをランダムに選択する。
ステップＳ１７では、電流素片推定処理部２は、ランダムに一つのパラメータすなわち電流素片Ｊの位置またはモーメントを選択する。
ステップＳ１８では、電流素片推定処理部２は、選択した電流素片Ｊがカバリング・ボールＶから出ないという制約下で、選択した電流素片Ｊの、選択したパラメータを、少し変更する。
【００２９】
ステップＳ１９では、電流素片推定処理部２は、ステップＳ３と同様に、多数の磁気センサの３次元空間位置にＡ個の電流素片が発生させる磁気ベクトルを算出し、３次元空間位置と磁気ベクトルとを対応させた計算磁場を得る。
【００３０】
ステップＳ２０では、電流素片推定処理部２は、公知の評価関数を用いて磁場データ（図４）と計算磁場の誤差を計算し、誤差が十分小さいならステップＳ１２に戻り、誤差が十分小さくないなら図２のステップＳ６に戻って電流素片Ｊのパラメータ変更からやり直す。
【００３１】
上記ステップＳ１２〜Ｓ２０を十分な回数繰り返すことにより、半径が十分小さくなるカバリング・ボールＶが定まるが、これは１つのカバリング・ボールＶに含まれる複数の電流素片Ｊが該カバリング・ボールＶの中心に集まった状態であり、図１０に示すように、それら複数の電流素片Ｊを合成して１つの電流素片Ｊ’と見なすことが出来る。つまり、元のＡ個の電流素片Ｊの分布を、Ｂ個の電流素片Ｊ’の分布に変換できる。このＢ個の電流素片Ｊ’の分布は、実際にありそうな解である可能性が高い。
なお、本発明は、３個以上（Ｂ≧３）の電流素片を仮定した場合に有用であり、２個以下の電流素片を仮定する場合は従来公知の方法で足る。
【００３２】
ステップＳ２１では、電流素片推定処理部２は、図１１に示すように、磁場データとＢ個の電流素片Ｊ’の分布とを視覚化して表示する。
【実施例２】
【００３３】
実施例１のステップＳ１４で、見込みがないならステップＳ１１に戻ってカバリング・ボールＶの発生からやり直してもよい。
【実施例３】
【００３４】
実施例１では、脳の磁場データから電流素片の分布を推定するのに本発明を適用したが、これに限定されず、例えば超音波エコーのデータから金属内部やコンクリート内部の傷（エコー発生源）の分布を推定したり、地震波のデータから地中や海中の地震源の分布を推定するのにも本発明を適用できる（特に、エコー発生源や地震源が３箇所以上である場合に有用である）。
【実施例４】
【００３５】
信号源の分布が３次元空間である場合には球体のカバリング・ボールを用いるが、信号源の分布が２次元空間である場合には円領域のカバリング・サークルを用いればよい。
【産業上の利用可能性】
【００３６】
例えば、脳磁気検出装置により測定した磁場データから、脳に存在すると仮定した電流素片の分布を推定するのに利用できる。
【図面の簡単な説明】
【００３７】
【図１】実施例１に係る電流素片推定装置を示す構成図である。
【図２】電流素片推定処理の手順を示すフロー図である。
【図３】図２の続きのフロー図である。
【図４】磁場データの概念図である。
【図５】発生した電流素片を示す概念図である。
【図６】計算磁場を示す概念図である。
【図７】発生したカバリング・ボールを示す概念図である。
【図８】カバリング・ボールの半径縮小を示す概念図である。
【図９】半径を十分小さくしたカバリング・ボールと電流素片を示す概念図である。
【図１０】半径を十分小さくしたカバリング・ボール中の複数の電流素片の合成を示す概念図である。
【図１１】視覚化した磁場データと電流素片の分布を示す概念図である。
【符号の説明】
【００３８】
１          データ収集部
２          電流素片推定処理部
３          表示部
１００      電流素片推定装置
Ｎ          脳磁気検出装置
【図１】
2005058558.tif 000003
【図２】
2005058558.tif 000004
【図３】
2005058558.tif 000005
【図４】
2005058558.tif 000006
【図５】
2005058558.tif 000007
【図６】
2005058558.tif 000008
【図７】
2005058558.tif 000009
【図８】
2005058558.tif 000010
【図９】
2005058558.tif 000011
【図１０】
2005058558.tif 000012
【図１１】
2005058558.tif 000013