1(57)【特許請求の範囲】【請求項１】 所定間隔毎の断層像に基づいて補間演算を行って立体像データを得、処理対象となる空間を所定の小立体領域に区分して各小立体領域に含まれる立体像データを計数し、計数値と所定の閾値との大小に基づいて小立体領域を２値化することを特徴とする断層像処理方法。
【請求項２】 所定間隔毎の断層像に基づいて補間演算を行って立体像データを得る断層像補間手段（３）と、処理対象となる空間を所定の小立体領域に区分する空間 10
区分手段（７）と、区分された各小立体領域に含まれる
立体像データを計数する立体像データ計数手段（８）と、計数値と所定の閾値との大小に基づいて小立体領域
を２値化する２値化手段（９）（１０）とを含むことを特徴とする断層像処理装置。

2【請求項３】 請求項１の断層像処理方法により得られた小立体領域の頂点の位置情報に基づいて異なる磁場源
情報に基づく複数のビオ・サバールの法則の演算を行
い、各演算結果を累積加算して得られる磁場演算値と測定した磁場計測値との差を算出し、算出した差に基づいて各演算式に含まれる磁場情報源を補正し、差が十分に小さくなるまで上記一連の処理を反復してから各演算式に含まれる磁場源情報を生体磁場測定結果として出力することを特徴とする生体磁場測定方法。
【請求項４】 算出した差に基づく補正に代えて、各小立体領域の各頂点の各座標成分の修正量を算出し、各座標成分の修正量の符号に基づいて修正対象頂点を得るとともに、各座標成分の修正量に基づいて、得られた修正対象頂点に対して磁場源情報測定値の再配分を行う処理を採用する請求項３に記載の生体磁場測定方法。

(2)3【請求項５】 請求項２の断層像処理装置により得られた小立体領域の頂点の位置情報に基づいて異なる磁場源
情報に基づく複数のビオ・サバールの法則の演算を行う
演算手段（１１１）（１１２）・・・（１１ｍ）（１０１）（１０２）・・・（１０ｍ）と、各演算手段（１１１）（１１２）・・・（１１ｍ）（１０１）（１０２）・・・（１０ｍ）から出力される演算結果を累積加算する累積加算手段（２０）（２００）と、累積加算手段（２０）（２００）から出力される累積加算結果と磁場計測値とを入力として誤差を算出する誤差算出手段（３０）（３００）と、算出誤差に基づいて各演算手段（１１１）（１１２）・・・（１１ｍ）（１０１）（１０
２）・・・（１０ｍ）における磁場源情報の補正を行う
補正手段（１１１ａ）（１１２ａ）・・・（１１ｍａ）（１２１）（１２２）・・・（１２ｍ）と、補正手段（１１１ａ）（１１２ａ）・・・（１１ｍａ）（１２１）（１２２）・・・（１２ｍ）による補正が施された結果を収集して生体磁場源測定結果として出力する補正
結果出力手段（４０）（５００）とを含むことを特徴とする生体磁場測定装置。
【請求項６】 補正手段として、各小立体領域の各頂点の各座標成分の修正量を算出する修正量算出手段（４１ｉ）と、各座標成分の修正量の符号に基づいて修正対象頂点を抽出する修正対象頂点抽出手段（４２ｉ）と、各座標成分の修正量に基づいて、得られた修正対象頂点に対して磁場源情報推定値の再配分を行う再配分手段（４３ｉ）（４４ｉ）（４５ｉ）（４６ｉ）とを含むものを用いる請求項５に記載の生体磁場測定装置。
【請求項７】 請求項３または請求項４の生体磁場測定
方法によって得られた生体磁場測定結果を請求項１の断層像処理方法によって得られた生体構造参照像と共に表示する方法であって、生体内の複数の点のそれぞれに対応する生体磁場測定値の絶対値に基づいて単位図形のサイズを設定し、生体構造参照像を表示するとともに、生体構造参照像の該当位置に設定されたサイズの単位図形を表示することを特徴とする生体磁場測定結果表示方法。
【請求項８】 単位図形が円である請求項７に記載の生
体磁場測定結果表示方法。
【請求項９】 請求項５または請求項６の生体磁場測定
装置によって得られた生体磁場測定結果を請求項２の断層像処理装置によって得られた生体構造参照像と共に表示する装置であって、生体内の複数の点のそれぞれに対応する生体磁場測定値の絶対値に基づいて単位図形のサイズを設定する単位図形サイズ設定手段（７００）と、生体構造参照像を表示する生体像表示手段（８００）と、生体構造参照像の該当位置に設定されたサイズの単位図形を表示する単位図形表示手段（９００）とを含むことを特徴とする生体磁場測定結果表示装置。
【請求項１０】単位図形が円である請求項９に記載の生

10

20

30

40

50

特許３０８４９５１
4
体磁場測定結果表示装置【発明の詳細な説明】【０００１】【産業上の利用分野】この発明は断層像処理方法、生体磁場測定方法、生体磁場測定結果表示方法およびこれらの装置に関し、さらに詳細にいえば、所定間隔毎の断層像に基づいてきめ細かい像データを得るための新規な断層像処理方法およびその装置、断層像処理方法により得られたきめ細かい像データに基づいて生体磁場源の解析を行なうための新規な生体磁場測定方法およびその装置、ならびに生体磁場源の解析結果を把握し易い状態で表示するための新規な生体磁場測定結果表示方法およびその装置に関する。
【０００２】【従来の技術】従来から生体磁場解析結果を把握し易い
状態で表示するために、ＭＲＩ画像等の断層像を可視的に表示した状態で磁場源を可視的に重畳して表示する方法が採用されている。また、磁場源を可視的に重畳表示する方法として、磁場源に最も近い断層像と重畳した状態で磁場源を２次元的に表示する方法および複数の断層像を３次元的に表示し、磁場源をも３次元的に表示する方法が提案されている。
【０００３】また、生体磁場測定を測定する方法として、磁場源がシングル・ダイポールであると仮定した場合に対処するための測定方法および磁場源が複数電流素片であると仮定した場合に対処するための測定方法が提案されている。前者の方法は、体表近傍における複数点の磁場測定結果に基づいてシングル・ダイポールの位置
および方向性等を算出する方法である。また、後者の方
法は、複数の電流素片（通常は３０〜１００程度の電流素片）の存在を仮定しておき、電流素片の数よりも著しく多い測定点における磁場測定結果に基づいて演算を行ない、上記複数の電流素片の位置および方向性を算出する方法である。
【０００４】また、生体磁気計測結果に基づく磁場源解
析を行なうことにより得られた電流素片を電流ベクトルとして表示するアローマップ表示、または、等磁線に基づくコンターマップ表示が知られている。そして、医療解析に用いる場合には、ＭＲＩ画像上に電流ベクトルを示す矢印、または等磁線をオーバーラップさせて表示することにより、臓器中における電流素片または等磁線を可視的に表示することが提案されている。
【０００５】【発明が解決しようとする課題】上記磁場源表示方法のうち、２次元的に表示する方法を採用した場合には、磁場源が何れかの断層像上に位置する可能性が著しく低いとともに、選択された断層像に対する磁場源の相対位置（特に奥行き方向の相対位置）が全く把握できないのであるから、医師等が２次元的な重畳表示を見ても磁場源の位置を正確に把握することが非常に困難になってしま(3)5うという不都合がある。また、断層像では生体内部の臓器の有無を示しているのであるが、臓器の奥行き方向の形状が断層像の形状からずれてしまう可能性が高いのであるから、実際には断層像上において臓器が存在しない
位置に磁場源が表示されてしまう可能性が高くなり、この点からも磁場源の位置を正確に把握することが非常に困難になってしまう。
【０００６】また、３次元的に表示する方法を採用した
場合には、基準になる画像として複数の断層像が存在するだけであるから、臓器との位置関係が非常に把握しにくいのみならず、断層像以外に目印がないのであるか
ら、磁場源の深さ、奥行きが非常に把握しにくいという不都合がある。上記生体磁場測定を測定する方法のう
ち、磁場源がシングル・ダイポールであると仮定した場合に対処するための測定方法を採用した場合には、測定結果としてシングル・ダイポールが得られるだけであるから、臓器の運動時に流れる体積電流に代表される分布性磁場源の測定には対処できないという不都合があるの
みならず、深すぎる解に収束してしまい勝ちであるという不都合もある。
【０００７】また、磁場源が複数電流素片であると仮定した場合に対処するための方法を採用した場合には、分布性磁場源の測定が可能になるのである。しかし、一般的に、３０〜１００程度の電流素片を想定することが必要になるのであるから、全ての電流素片の測定を行なおうとすれば、２００〜１０００チャネル程度のアレイセンサを設けて多数の磁場計測値を得、得られた多数の磁場計測値に基づいて各電流素片の推定処理のための演算を行なう必要がある。そして、上述の場合に対する推定処理のための演算負荷はテラ・フロップスが要求されるのであるから、システム全体として著しく大規模になるとともに、高価になってしまうという不都合がある。また、推定処理においてローカルミニマに収束してしまう
確率が高く、正確な電流素片の推定を達成できる確率が著しく減少してしまうという不都合もある。
【０００８】また、上記アローマップ表示、等磁線表示の何れも特定のタイミングにおける生体磁場測定値、即ち分布を可視的に表示するものであり、動的な生体磁場
測定値の表示には不適切である。さらに詳細に説明すると、アローマップ表示により動的に電流素片を表示することは可能であるが、一般的に臓器における電流素片の数は、体積電流を表示するのに十分な数に設定しなければならない。即ち、著しく多数の電流素片がベクトル表
示されなければならない。そして、これら全ての電流素
片の向き、大きさが時々刻々と変化するのであるから、ＭＲＩ画像とオーバーラップ表示された矢印を見ても、例えば、臓器の興奮電流の伝播状態を簡単には把握できないという不都合がある。
【０００９】コンターマップ表示を動的な生体磁場測定
値の表示に適用しても、等磁線はもともと個々の電流素

10

20

30

40

50

特許３０８４９５１
6片の表示には不向きであるから、到底臓器の興奮電流の伝播状態の把握のために採用することはできない。
【００１０】【発明の目的】この発明は上記の問題点に鑑みてなされたものであり、所定間隔毎の断層像に基づいてきめ細かい像データを得ることができる新規な断層像処理方法およびその装置を提供することを第１の目的とし、断層像処理方法により得られたきめ細かい像データに基づいて簡単にかつ正確に生体磁場源の解析を行なうことができる生体磁場測定方法およびその装置を提供することを第２の目的とし、動的な生体物理量測定値を把握し易い状
態で表示できる生体物理量計測結果表示方法およびその装置を提供することを第３の目的としている。
【００１１】【課題を解決するための手段】上記第１の目的を達成するための、請求項１の断層像処理方法は、所定間隔毎の断層像に基づいて補間演算を行って立体像データを得、処理対象となる空間を所定の小立体領域に区分して各小立体領域に含まれる立体像データを計数し、計数値と所定の閾値との大小に基づいて小立体領域を２値化する方法である。
【００１２】請求項２の断層像処理装置は、所定間隔毎の断層像に基づいて補間演算を行って立体像データを得
る断層像補間手段と、処理対象となる空間を所定の小立体領域に区分する空間区分手段と、区分された各小立体領域に含まれる立体像データを計数する立体像データ計
数手段と、計数値と所定の閾値との大小に基づいて小立
体領域を２値化する２値化手段とを含んでいる。
【００１３】上記第２の目的を達成するための、請求項３の生体磁場測定方法は、請求項１の断層像処理方法により得られた小立体領域の頂点の位置情報に基づいて異なる磁場源情報に基づく複数のビオ・サバールの法則の演算を行ない、各演算結果を累積加算して得られる磁場演算値と測定した磁場計測値との差を算出し、算出した差に基づいて各演算式に含まれる磁場源情報を補正し、差が十分に小さくなるまで上記一連の処理を反復してから各演算式に含まれる磁場源情報を生体磁場測定結果として出力する方法である。
【００１４】請求項４の生体磁場測定方法は、算出した差に基づく補正に代えて、各小立体領域の各頂点の各座
標成分の修正量を算出し、各座標成分の修正量の符号に基づいて修正対象頂点を得るとともに、各座標成分の修正量に基づいて、得られた修正対象頂点に対して磁場源
情報推定値の再配分を行なう処理を採用する方法である。
【００１５】請求項５の生体磁場測定装置は、請求項２の断層像処理装置により得られた小立体領域の頂点の位
置情報に基づいて異なる磁場源情報に基づく複数のビオ・サバールの法則の演算を行う演算手段と、各演算手段
から出力される演算結果を累積加算する累積加算手段

(4)7と、累積加算手段から出力される累積加算結果と磁場計
測値とを入力として誤差を算出する誤差算出手段と、算出誤差に基づいて各演算手段における磁場源情報の補正を行う補正手段と、補正手段による補正が施された結果を収集して生体磁場源測定結果として出力する補正結果
出力手段とを含んでいる。
【００１６】請求項６の生体磁場測定装置は、補正手段
として、各小立体領域の各頂点の各座標成分の修正量を算出する修正量算出手段と、各座標成分の修正量の符号に基づいて修正対象頂点を抽出する修正対象頂点抽出手
段と、各座標成分の修正量に基づいて、得られた修正対
象頂点に対して磁場源情報推定値の再配分を行なう再配分手段とを含むものを用いている。
【００１７】上記第３の目的を達成するための、請求項７の生体磁場測定結果表示方法は、請求項３または請求
項４の生体磁場測定方法によって得られた生体磁場測定
結果を請求項１の断層像処理方法によって得られた生体構造参照像と共に表示する方法であって、生体内の複数の点のそれぞれに対応する生体磁場測定値の絶対値に基づいて単位図形のサイズを設定し、生体構造参照像を表示するとともに、生体構造参照像の該当位置に設定されたサイズの単位図形を表示する方法である。請求項８の
生体磁場測定結果表示方法は、円からなる単位図形を表示する方法である。
【００１８】請求項９の生体磁場測定結果表示装置は、
請求項５または請求項６の生体磁場測定装置によって得られた生体磁場測定結果を請求項２の断層像処理装置によって得られた生体構造参照像と共に表示する装置であって、生体内の複数の点のそれぞれに対応する生体磁場
測定値の絶対値に基づいて単位図形のサイズを設定する
単位図形サイズ設定手段と、生体構造参照像を表示する
生体像表示手段と、生体構造参照像の該当位置に設定されたサイズの単位図形を表示する単位図形表示手段とを含んでいる。
【００１９】【作用】請求項１の断層像処理方法であれば、所定間隔毎の断層像に基づいて補間演算を行って立体像データを得、処理対象となる空間を所定の小立体領域に区分して各小立体領域に含まれる立体像データを計数し、計数値と所定の閾値との大小に基づいて小立体領域を２値化するのであるから、所定間隔毎の断層像しか得られていな
くても、補間演算を行うことにより全空間に対応して立体像データを得、しかも区分された小立体領域に含まれる立体像データの数と所定の閾値との大小に基づいて小立体領域を２値化することにより単純化された立体像を得ることができる。したがって、可視的に表示した場合に、臓器等に相当する箇所、空間に相当する箇所を簡単に認識できる。
【００２０】尚、上記補間演算としては、直線補間演
算、スプライン補間演算等から適宜選択された補間演算

10

20

30

40

50

特許３０８４９５１8を採用すればよい。また、小立体領域を２値化する場合
における立体像データの数についても断層像の種類、対象臓器等に適合させるべく適宜閾値を設定すればよい。
さらに、小立体領域としては、空間を密に充填し得る立体領域であればよいが、後処理の簡素化等を考慮すれば、立方体領域であることが最も好ましい。
【００２１】請求項２の断層像処理装置であれば、所定間隔毎の断層像に基づいて断層像補間手段により補間演算を行って立体像データを得る。また、空間区分手段に
より、処理対象となる空間を所定の小立体領域に区分しておく。そして、立体像データ計数手段により、区分された各小立体領域に含まれる立体像データを計数し、２値化手段により、計数値と所定の閾値との大小に基づいて小立体領域を２値化する。したがって、可視的に表示した場合に、臓器等に相当する箇所、空間に相当する箇所を簡単に認識できる。請求項３の生体磁場測定方法であれば、請求項１の断層像処理方法により得られた小立
体領域の頂点の位置情報に基づいて異なる磁場源情報に基づく複数のビオ・サバールの法則の演算を行い、各演算結果を累積加算して得られる磁場演算値と測定した磁場計測値との差を算出し、算出した差に基づいて各演算式にに含まれる磁場源情報を補正し、差が十分に小さくなるまで上記一連の処理を反復してから、複数のビオ・サバールの法則の演算結果の累積加算値が磁場計測値と高精度に近似されるまで補正を行うことにより、各小立
体領域の頂点における磁場源情報を得ることができる。
また、小立体領域の頂点の位置情報は既知であるから、推定すべき未知数の数を大幅に低減できる。
【００２２】請求項４の生体磁場測定方法であれば、各小立体領域の各頂点の各座標成分の修正量を算出し、各座標成分の修正量の符号に基づいて修正対象頂点を得る
とともに、各座標成分の修正量に基づいて、得られた修正対象頂点に対して磁場源情報推定値の再配分を行なう
ことにより、算出した差に基づく補正を行なうのであるから、互に近接する状態で割り当てられた、大きさが互に等しくかつ向きが互に逆の磁場源ベクトルが存在していても、これらを確実に補正でき、一層高精度に推定された磁場源情報を得ることができる。また、隣接する頂点間で再配分を行なうことにより推定を行なうのであるから、推定対象となる磁場源情報の数に見合って従来必要とされていた測定点数よりも少ない測定点に基づいて高精度の磁場源推定を達成できる。
【００２３】請求項５の生体磁場測定装置であれば、請求項２の断層像処理装置により得られた小立体領域の頂点の位置情報に基づいて、演算手段により異なる磁場源
情報に基づく複数のビオ・サバールの法則の演算を行って、累積加算手段により、各演算手段から出力される演
算結果を累積加算する。そして、累積加算手段から出力される累積加算結果と磁場計測値とに基づいて誤差算出
手段により誤差を算出し、算出誤差に基づいて補正手段

(5)9により各演算手段における磁場源情報の補正を行う。補正手段による補正処理を必要回数だけ行った後に、補正結果収集手段により、補正手段による補正が施された結果を収集して生体磁場源測定結果として出力することにより、各小立体領域の頂点における磁場源情報を得ることができる。また、小立体領域の頂点の位置情報は既知であるから、推定すべき未知数の数を大幅に低減でき、装置全体として構成を簡素化できる。
【００２４】請求項６の磁場源測定装置であれば、補正
量算出手段により各小立体領域の各頂点の各座標成分の修正量を算出し、各座標成分の修正量の符号に基づいて修正対象頂点抽出手段により修正対象頂点を抽出し、再配分手段により、各座標成分の修正量に基づいて、得られた修正対象頂点に対して磁場源情報推定値の再配分を行なうことにより、算出した差に基づく補正を行なうのであるから、互に近接する状態で割り当てられた、大きさが互に等しくかつ向きが互に逆の磁場源ベクトルが存在していても、これらを確実に補正でき、一層高精度に推定された磁場源情報を得ることができる。また、隣接する頂点間で再配分を行なうことにより推定を行なうのであるから、推定対象となる磁場源情報の数に見合って従来必要とされていた測定点数よりも少ない測定点に基づいて高精度の磁場源推定を達成できる。
【００２５】生体磁場測定方法についてさらに詳細に説明すると、本件発明者が生体磁場の測定に関して鋭意研究を重ねた結果、次のことを見出した。各小立体領域の各頂点は断層像処理方法により得られているのであるか
ら、座標値は既知である。したがって、各頂点における磁場源情報を得るために必要な情報は電流素片のベクトル成分のみである。
【００２６】また、１つの小立体領域の１つの頂点は、隣合う他の小立体領域の頂点を兼ねているのであるから、小立体領域の数が増加しても、頂点数が小立体領域の増加に比例して増加することはなく、小立体領域の数が十分に大きければ、頂点数と小立体領域の数はほぼ等しくなる。さらに、１つの頂点の電流素片は近接する他の電流素片と相関があるので、１つの頂点の電流素片当りの情報量は一般的な未知数量の推定を行なう場合と比較して著しく少なくできる。
【００２７】さらにまた、磁場源でない頂点における電流素片は０ベクトルに収束するので、最終的に推定を行
なわなければならない電流素片の数は著しく少なくな
る。請求項３から請求項６の発明は上記の知見に基づいてなされたものであり、測定開始時点においては未知の電流素片（数も未知）の推定を行なうに当って、未知の電流素片数を考慮することなく、計測システムの能力、容量等に基づいて定まる所定数の磁場計測値を得、所定数の磁場計測値に基づいて磁場源情報の補正処理を反復することにより未知の電流素片の高精度の推定を達成でき、推定結果を生体磁場測定結果として採用できる。

10

20

30

40

50

特許３０８４９５１
10【００２８】請求項７の生体磁場測定結果表示方法であれば、請求項３または請求項４の生体磁場測定方法によって得られた生体磁場測定結果を請求項１の断層像処理
方法によって得られた生体構造参照像と共に表示するに当たって、生体内の複数の点のそれぞれに対応する生体
磁場測定値の絶対値に基づいて単位図形のサイズを設定
し、生体構造参照像を表示するとともに、生体構造参照
像の該当位置に設定されたサイズの単位図形を表示するのであるから、矢印の方向と長さで磁場測定結果を表示する場合と比較して、例えば、最大の生体磁場測定結果の分布状態を簡単に認識できる。但し、生体磁場の方向成分に関しては把握できないことになるが、生体磁場測定結果を動的に表示すれば、例えば、最大の生体磁場測定結果の分布状態の時間的変動を簡単に認識できるのであるから、特に不都合はない。生体磁場測定結果の表示が最も意味を持つ医療解析においては、一般に最大の生体磁場測定結果の分布および時間的変動を簡単に把握できることが要望されるのであるから、この生体磁場測定
結果表示方法はこの要望に最適の方法である。
【００２９】請求項８の生体磁場測定結果表示方法であれば、請求項７の作用に加え、単位図形が円であるか
ら、単位図形同士の隙間をかなり大きくでき、背景となる生体構造参照像の把握が容易になる。請求項９の生体磁場測定結果表示装置であれば、請求項５または請求項６の生体磁場測定装置によって得られた生体磁場測定結果を請求項２の断層像処理装置によって得られた生体構
造参照像と共に表示するに当たって、生体内の複数の点のそれぞれに対応する生体磁場測定値の絶対値に基づいて単位図形サイズ設定手段により単位図形のサイズを設定し、生体像表示手段により生体構造参照像を表示するとともに、単位図形表示手段により生体構造参照像の該
当位置に設定されたサイズの単位図形を表示するのであるから、矢印の方向と長さで磁場測定結果を表示する場合と比較して、例えば、最大の生体磁場測定結果の分布
状態を簡単に認識できる。但し、生体磁場の方向成分に関しては把握できないことになるが、生体磁場測定結果
を動的に表示すれば、例えば、最大の生体磁場測定結果の分布状態の時間的変動を簡単に認識できるのであるか
ら、特に不都合はない。生体磁場測定結果の表示が最も意味を持つ医療解析においては、一般に最大の生体磁場測定結果の分布および時間的変動を簡単に把握できることが要望されるのであるから、この生体磁場測定結果表
示装置はこの要望に最適の装置である。
【００３０】請求項１０の生体磁場測定結果表示装置であれば、請求項９の作用に加え、単位図形が円であるか
ら、単位図形同士の隙間をかなり大きくでき、背景となる生体構造参照像の把握が容易になる。
【００３１】【実施例】以下、実施例を示す添付図面によって詳細に説明する。図１はこの発明の断層像処理方法の一実施例

(6)11を説明するフローチャートであり、ステップＳＰ１において断層像が与えられるまで待ち、ステップＳＰ２において、与えられた断層像に基づく補間処理を行なって補間断層像を得、ステップＳＰ３において、与えられた断層像および補間断層像を予め設定されたサイズの正方形に量子化し、ステップＳＰ４において、量子化された正方形の中心を中心とし、かつ１辺が正方形の１辺と等しい立方体を得、ステップＳＰ５において、何れかの立方体について、与えられた断層像を構成する画素および補間断層像を構成する画素の数を計数し、ステップＳＰ６において計数値が予め設定されている閾値（例えば、全平均値の平均値）よりも大きいか否かを判別する。そして、計数値が閾値よりも大きいと判別された立方体についてはステップＳＰ７において表示されるべき立方体（以下、表示立方体と称する）に割り当て、逆にステップＳＰ６において計数値が閾値よりも大きくないと判別された立方体についてはステップＳＰ８において表示されない立方体（以下、非表示立方体と称する）に割り当てる。ステップＳＰ７またはステップＳＰ８の処理が行なわれた後は、ステップＳＰ９において全ての立方体について処理が完了したか否かを判別し、処理が行なわれていない立方体が存在すると判別された場合には、ステップＳＰ１０において他の立方体を選択し、再びステッ
プＳＰ５の処理を行なう。また、ステップＳＰ９において全ての立方体について処理が完了したと判別された場
合には、ステップＳＰ１１において、表示立方体に割り当てられた立方体のみを可視的に表示して一連の処理を終了する。
【００３２】但し、立方体に含まれる画素数に代えて、立方体中の全画素数に対する断層像画素数の割合および対応する閾値を採用してもよいことはもちろんである。
図２から図４を参照しながら断層像処理方法をさらに説明する。図２は断層像の一例を示す概略図であり、ぬり
つぶし表示された部分が人体の臓器を示している。図２に示す断層像について量子化された正方形の中心が小さいドットで示されるように与えられるので、与えられた
ドットを中心とする立方体が定まることになる。そして、定められた立方体に含まれる画素数が閾値よりも大きい場合には図３に平面視して示すように表示立方体と非表示立方体とに割り当てられる。図４は図３に対応する斜視図であり、表示立方体に割り当てられた立方体のみが表示されている。したがって、表示された立方体を見ることにより、臓器の立体形状を簡単に把握できる。
【００３３】尚、この一連の処理を全ての断層像および補間断層像について行なうことにより、生体臓器をソリッド・モデルとして可視的に表示できる。但し、磁場源
解析結果を重畳表示する場合には、磁場源解析結果を確認し易くするために、表示立方体に割り当てられた立方体を半透明表示することが好ましい。また、断層像中における各臓器を指定することにより、該当する臓器のみ

10

20

30

40

50

特許３０８４９５１
12の表示、複数の臓器の色分け表示をも簡単に達成できる。さらに、非表示立方体は臓器の空洞等に対応するのであるから、後述する生体磁場源解析等に当っては表示
立方体のみについて処理を行なえばよくなり、後処理に必要な演算負荷を大幅に低減できる。
【００３４】【実施例２】図５はこの発明の断層像処理装置の一実施
例を示すブロック図であり、複数の断層像を保持する断層像保持部１と、断層像保持部１から処理対象となる断層像を選択する第１断層像選択部２と、選択された断層像に基づく補間処理を行なって補間断層像を得る断層像
補間部３と、補間断層像を保持する補間断層像保持部４と、断層像選択部２および断層像補間部３を必要回数だ
け反復動作させる第１反復制御部５と、第１反復制御部
５による反復動作制御が終了したことに応答して、断層像保持部１または補間断層像保持部４から順次断層像を選択する第２断層像選択部６と、選択された断層像に対応して立方体を割り当てる立方体割り当て部７と、割り当てられた立方体中に存在する画素数を計数する画素数
計数部８と、計数された画素数が所定の閾値よりも大きいか否かを判別する判別部９と、画素数が所定の閾値よりも大きいことを示す判別部９の判別結果に応答して該当する立方体を表示立方体に割り当てる表示立方体割当部１０と、画素数が所定の閾値よりも大きくないことを示す判別部９の判別結果に応答して該当する立方体を非表示立方体に割り当てる非表示立方体割当部１１と、断層像、補間断層像に対応して立方体の割当結果を保持する割当結果保持部１２と、第２断層像選択部６、立方体
割当部７、画素計数部８、判別部９、表示立方体割当部
１０および非表示立方体割当部１１を必要回数だけ反復動作させる第２反復制御部１３と、第２反復制御部１３により反復動作制御が終了したことに応答して表示立方体に割り当てられた立方体のみを可視的に表示する立方体表示部１４とを有している。
【００３５】したがって、図５の構成の断層像処理装置を採用することにより、図２から図４に示す断層像の処理および表示立方体に割り当てられた立方体のみに基づ
く可視的表示を達成できる。
【００３６】【実施例３】図６はこの発明の生体磁場測定装置の一実施例を示すブロック図であり、複数個の、ビオ・サバー
ルの法則の演算を行なうビオ・サバール演算ユニット１
１１，１１２，・・・，１１ｍと、ビオ・サバール演算
ユニット１１１，１１２，・・・，１１ｍから出力される演算結果を累積加算するシグマ・ユニット２０と、シグマ・ユニット２０から出力される累積加算結果Ｏｊ（ｔ）と教師パターンとしての磁場計測値Ｓｊ（ｔ）とを入力として両者の差を算出する誤差演算器３０と、算出された差に基づいて、ビオ・サバール演算ユニット１
１１，１１２，・・・，１１ｍにおいて推定されている

(7)

特許３０８４９５１1314
電流素片のベクトル成分を補正する補正部１１１ａ，１＊ たは図５の断層像処理装置により得られているのである
１２ａ，・・・，１１ｍａと、ビオ・サバール演算ユニから、３次元座標値は全て既知の値である。
ット１１１，１１２，・・・，１１ｍにおいて推定され【００３８】したがって、時刻ｔおよび各小立方体の各ている電流素片のベクトル成分を収集して解析結果とし頂点の３次元座標値をｍ個のビオ・サバール演算ユニッて出力する情報収集ユニット４０とを有している。尚、
ト１１１，１１２，・・・，１１ｍに供給してそれぞれ
上記ビオ・サバール演算ユニット１１１，１１２，・・
関数ｇ１，ｇ２，・・・，ｇｍの演算を行なって関数値
・，１１ｍは、時刻ｔ、図１の断層像処理方法もしくはを算出し、算出された全ての関数値をシグマ・ユニット
図５の断層像処理装置により得られた、表示立方体に割２０に供給することにより累積加算値Ｏｊ（ｔ）を得る
り当てられた立方体の頂点の位置情報が既知情報としてことができる。但し、当初は各電流素片のベクトル成分供給されたことに応答して、既知情報に基づいて各ビオ 10 が適当に設定されているのであるから、得られる累積加・サバール演算ユニットに設定されている、ビオ・サバ
算値Ｏｊ（ｔ）は実際の計測値Ｓｊ（ｔ）とは異なる。
ールの法則の演算を行なうとともに、誤差演算器３０かしたがって、誤差演算器３０において実際の計測値Ｓｊ
ら出力される推定誤差ｄｊ（ｔ）｛＝Ｓｊ（ｔ）−Ｏｊ（ｔ）と累積加算値Ｏｊ（ｔ）との差を算出し、算出さ（ｔ）｝が供給されたことに応答して演算式に含まれるれた差を推定誤差ｄｊ（ｔ）としてビオ・サバール演算
暖流素片のベクトル成分を推定誤差が少なくなるように
ユニット１１１，１１２，・・・，１１ｍの補正部１１補正する。また、ビオ・サバール演算ユニット１１１，
１ａ，１１２ａ，・・・，１１ｍａにフィードバック
１１２，・・・，１１ｍはそれぞれ同期的に動作するよし、推定誤差ｄｊ（ｔ）が小さくなるように各ビオ・サ
うに制御してもよく、また、非同期的に動作するように
バール演算ユニットのベクトル成分を変化させる。
制御してもよい。
【００３９】上記一連の処理を反復すれば推定誤差ｄｊ【００３７】上記の構成の生体磁場測定装置の作用は次 20 （ｔ）が小さくなり、ついには推定誤差ｄｊ（ｔ）がほのとおりである。解析対象となる磁場Ｏｊ（ｔ）は、時ぼ０になるので、この時点においてビオ・サバール演算
刻ｔと各小立方体の頂点の３次元座標値ｘｉ，ｙｉ，ｚ
ユニットのベクトル成分の値を情報収集ユニット４０に
ｉ（ｉは正の整数、以下同じ）および電流素片のベクトより収集して出力することにより各小立方体の各頂点にル成分Ｐｘｉ，Ｐｙｉ，Ｐｚｉを持つ関数ｇｉ（ビオ・おける電流素片のベクトル成分を得ることができる。まサバールの法則に基づいて定まる関数）の線形和であた、推定誤差評価関数Ｅｊ（ｔ）を次式で定義すれば、る。但し、磁場センサとして平面内のベクトル成分のみ
数１が得られる。
2を計測可能なものを使用する場合には、座標系をセンサ
Ｅｊ（ｔ）＝（１／２）｛Ｓｊ（ｔ）−Ｏｊ（ｔ）｝に適合させることにより何れかのベクトル成分をオミッ【００４０】
トできる。また、各小立方体は図１の断層像処理方法ま＊【数１】【００４１】そして、各ビオ・サバール演算ユニットに※ 但し、εｋはベクトル成分ａｉｋの学習ゲイン（補正ゲ
おけるベクトル成分の補正を最急降下法に基づいて行なイン）である。
うこととすれば、、推定誤差評価関数が最小になるベク【００４２】
トル成分の推定は数２に基づいて行なうことができる。※【数２】【００４３】また、累積加算値Ｏｊ（ｔ）のａｉｋによ 40★ 【００４４】
る偏微分値は数３で与えられるので、数２は数４と表現
【数３】
できる。
★

【００４５】

【数４】

(8)
15

【００４６】したがって、数４の処理を行なうことによりベクトル成分の推定精度を高め、より正確なベクトル成分を得ることができる。尚、推定誤差ｄｊ（ｔ）の傾きが正の場合には補正値を負に、推定誤差ｄｊ（ｔ）の傾きが負の場合には補正値を正に設定すればよい。以上のようにして各頂点の電流素片のベクトル成分が得られれば、８頂点における電流素片に基づいてベクトルの空間補間処理を行なって各微小点が発生する磁束密度を得、各微小点が発生する磁束密度を積分することにより、各小立方体が発生する磁界を得ることができる。したがって、得られた結果に基づいて分布電流を霧状に表示でき、または矢印によるアロー図で表示できる。
【００４７】また、図６の装置に基づく処理は図１または図５の実施例により得られた小立方体（表示立方体に割り当てられた小立方体）のみに基づいて行なえばよいのであるから、臓器等を対象とする場合に、臓器の内部
空洞部分等は当初から処理対象に含まれないことにな
り、処理負荷を低減できる。また、各小立方体の各頂点はかなり高い確率で他の小立方体の頂点を兼ねているのであるから、小立方体の数が増加割合に比して処理対象となる頂点数の増加割合が小さくなり、この点からも処
理負荷を低減できる。
【００４８】さらに、磁場計測値を得るための測定点数については、一般的な逆問題の解析を行なう場合に未知数の数よりも多い数の測定点が要求されると思われていたのであるが、比較的少ない測定点数で十分なベクトル成分の推定を達成できることを見出した。したがって、
測定点数を比較的少なくできることに伴なって全体としての処理負荷を低減できることになる。
【００４９】この点についてさらに詳細に説明する。従来は臓器等を小立方体の集合として定義することは全く行なわれておらず、したがって、例えば、推定すべき電流素片の数よりも少なくない測定点が必要であった。この結果、電流素片の数が定まっていなければ、十分な安全性を見込んで測定点数を設定していた。具体的には、臓器に流れる電流は時間の経過に伴なって位置および向き、大きさが変化するのであるから、電流素片の数は膨大にならざるを得ない。しかし、この実施例においては、各小立方体の位置が予め固定されているのであるか
ら、測定点に影響を及ぼさない小立方体における電流素
片は０ベクトルに収束してしまうことになる。そして、０ベクトルに収束した電流素片については以後の推定処理が不要になるのであるから、比較的少ない測定点で得られる磁場計測値のみに基づいて高精度のベクトル成分の推定を達成できる。
【００５０】【実施例３】図７はこの発明の生体磁場測定装置の他の特許３０８４９５１
16

10

20

30

40

50

実施例を示すブロック図、図８は１つの階層型パーセプ
トロンに対応する部分の構成を詳細に示す概略図であり、多入力１出力の複数個の階層型パーセプトロン１０ｐ（ｐ＝１，２，・・・，ｍ）と、階層型パーセプトロ
ン１０ｐからの関数値出力ｇｉｊを累積加算するシグマ・ユニット２００と、シグマ・ユニット２００から出力される累積加算結果Ｏｊ（ｔ）と教師パターンとしての計測値Ｓｊ（ｔ）とを入力として両者の差を算出する関
数値誤差演算器３００と、各階層型パーセプトロン１０ｐからの関数値出力ｇｉｊに対応する偏微分値を算出す
る偏微分値算出部１０ｐｒ（ｒ＝１，２，・・・ｎ）と、各階層型パーセプトロン１０ｐの偏微分値算出部１０ｐｒからの出力と数値微分法等により予め算出されている偏微分値（偏微分値教師パターン）とを入力として両者の差を算出する偏微分値誤差演算器３０ｐｒと、関数値誤差演算器３００により算出される差および偏微分値誤差演算器３０ｐｒにより算出される差に基づいて該当する階層型パーセプトロン１０ｐに偏微分関数の学習を行なわせる偏微分関数学習部６０ｐと、関数値誤差演算器３００により算出された累積加算結果Ｏｊ（ｔ）と教師パターンとしての計測値Ｓｊ（ｔ）との差および偏微分値算出部１０ｐｒにより算出された偏微分値に基づいて階層型パーセプトロン１０ｐの入力層における入力を補正する補正部１２ｐと、補正部１２ｐと偏微分関数
学習部６０ｐとを選択するとともに、選択された補正部
１２ｐまたは偏微分関数学習部６０ｐを所定回数（差を十分に小さくできる回数）だけ反復動作させる制御部４００と、補正部１２ｐにより補正が所定回数だけ反復されたことを条件として各階層型パーセプトロン１０ｐの未知のベクトル成分に対応する入力を未知数量推定結果
として出力する情報収集ユニット５００とを有している。 尚、上記各階層型パーセプトロン１０ｐは、既知の入力パターンおよび対応する教師パターンを与えて十分にビオ・サバールの法則に基づく関数の学習を行なったものである。また、階層型パーセプトロン１０ｐはそれぞれ同期的に動作するように制御してもよく、また、
非同期的に動作するように制御してもよい。
【００５１】また、偏微分関数学習部６０ｐは、例えば、関数値誤差演算器３００により算出される差に基づく階層型パーセプトロン１０ｐの学習（例えば、バック
プロパゲーション学習）および偏微分値誤差演算器３０ｐｒにより算出される差に基づく階層型パーセプトロン
１０ｐの学習を交互に反復させるものであり、関数の学
習結果を余り損なうことなく偏微分関数の学習を行なわせることができる。
【００５２】さらに、補正部１２ｐはそれぞれ複数の補正部を有しており、関数値誤差演算器３００により算出

(9)17された差および偏微分値算出部１０ｐｋにより算出された差を入力として該当する未知数量を補正する。上記の構成の未知数量推定装置の作用は次のとおりである。各
階層型パーセプトロン１０ｐにおいてビオ・サバールの法則に基づく学習が完了しているのであるから、未知のベクトル成分に対応する既知の入力の何れか１つを微小量だけ変化させて出力パターンを得るとともに、偏微分
値算出部１０ｐｒにより対応する偏微分値を得る。尚、
未知のベクトル成分に対応する既知の入力を微小量だけ変化させた場合における偏微分値を数値微分法等により予め算出して偏微分値教師パターンとして与えておく。
この状態において制御部４０により偏微分関数学習部６
０ｐを選択し、選択した偏微分関数学習部６０ｐを所定回数だけ反復動作させることにより再びバックプロパゲーション学習を行ない、ビオ・サバールの法則に基づく関数およびこの関数の偏微分関数の演算と等価な処理を達成できるように各階層型パーセプトロン１０ｐを構成するニューロン素子の荷重、閾値を決定する。
【００５３】尚、何れかの階層型パーセプトロン１０ｐ
により得られた学習結果としての荷重、閾値を他の階層型パーセプトロンの荷重、閾値としてそのまま採用すれば、他の階層型パーセプトロンについては関数の学習および偏微分関数の学習を共に不要にでき、学習所要時間を大幅に短縮できる。以上のようにして必要な学習が完了した後は、偏微分関数学習部６０Ｐに代えて制御部４００により補正部１２ｐを選択するとともに情報収集ユニット５００を動作させればよく、図６の実施例と同様にベクトル成分を精度よく推定できる。
【００５４】【実施例４】図９はこの発明の生体磁場測定方法のさらに他の実施例を示すブロック図であり、図７および図８の実施例と異なる点は、図７の実施例において階層型パーセプトロン毎に設けられている偏微分値算出部１０ｐｒに代えて、１入力のみを微小量だけ変化させる入力変化部１５ｐと、１入力の変化前後における出力を保持す
る出力保持部１６ｐと、両出力の差を算出する差算出部
１７ｐを設けた点のみである。
【００５５】したがって、この実施例の場合には、偏微分演算を行なわなくても、１入力のみの微小変化に対応
する出力変化を算出でき、この出力変化が偏微分値に相当するのであるから、上記実施例と同様の作用を達成できる。また、この実施例では偏微分関数の学習が不要になるので、学習所要時間を大幅に短縮できる。
【００５６】【実施例５】図１０はこの発明の磁場源測定装置のさらに他の実施例を示すブロック図であり、１ユニット分のみを示している。複数個の、ビオ・サバールの法則の演算を行なうビオ・サバール演算ユニット１１ｉと、ビオ・サバール演算ユニット１１ｉにおいて推定されている
電流素片のベクトル成分を収集して解析結果として出力

10

20

30

40

50

特許３０８４９５１18する電流素片収集ユニット４０ａと、図１の方法または
図５の装置により得られた３次元座標値を基準として修
正量Δｘ，Δｙ，Δｚを算出する修正量算出部４１ｉと、算出された修正量Δｘ，Δｙ，Δｚの符号に基づいて、推定されている電流素片のベクトル成分の再配分対象となる隣接頂点を抽出する頂点抽出部４２ｉと、修正
量Δｘ，Δｙ，Δｚに基づいて定まる修正ベクトルを正規化する正規化部４３ｉと、正規化された修正ベクトルの各成分の絶対値に基づいて再配分比率を算出する再配分比率算出部４４ｉと、推定されている電流素片のベク
トル成分および算出された再配分比率に基づいて抽出された頂点に関して推定されている電流素片のベクトル成分に推定値の再配分を行なう再配分部４５ｉと、再配分
部４５ｉにより再配分された再配分値に基づいて該当する電流素片のベクトル成分の推定値を更新する推定値更
新部４６ｉと、全ての頂点について修正量Δｘ，Δｙ，Δｚが十分に小さくなるまでビオ・サバール演算ユニッ
ト１１ｉ、修正量算出部４１ｉ、頂点抽出部４２ｉ、正規化部４３ｉ、再配分比率算出部４４ｉ、再配分部４５ｉおよび推定値更新部４６ｉの処理を反復させる反復制
御部４７とを有している。
【００５７】尚、上記推定値更新部４６ｉにおいては、
複数の隣接頂点から配分を受ける場合に、各隣接頂点からの再配分比率の加算を行なうのではなく、全ての再配分比率の平均値を算出して現在の値（他の隣接頂点に再配分した残りの割合）に加算することにより新たな推定値を得るようにしている。上記修正量算出部４１ｉは、
例えば、ビオ・サバール演算ユニットにより算出された関数値の偏微分値を算出することにより修正量Δｘ，Δｙ，Δｚを得るものである。
【００５８】上記の構成の生体磁場測定装置の作用は次のとおりである。各ビオ・サバール演算ユニットにおいては、
ｇｉｊ＝（μ０／４π）｛（Ｙｊｉ×ｐｘｉ−Ｘｊｉ×3ｐｙｉ）／Ｒｊｉ ｝の演算を行なう。ここで、μ０は透磁率を示し、ｉ（ｉ＝１，２，・・・，Ｎ）は各小立方体の各頂点における電流素片を、ｊ（ｊ＝１，２，・・・，ｍ）は測定点をそれぞれ示し、測定点の計測条件をＰｊ＝（ｘｊ，ｙｊ，ｚｊ）、電流素片の未知の物理量をＵｉ＝（ｘｉ，ｙｉ，ｚｉ，ｐｘｉ，ｐｙｉ）とした場合にＸｊｉ＝ｘ
ｊ−ｘｉ、Ｙｊｉ＝ｙｊ−ｙｉ、Ｚｊｉ＝ｚｊ−ｚｉで
2
2
2
1/2
あり、Ｒｊｉ＝（Ｘｊｉ ＋Ｙｊｉ ＋Ｚｊｉ ） である。
【００５９】したがって、ビオ・サバール演算ユニット
により得られる結果を対応する修正量算出部４１に供給すれば、上記演算式を各未知数で偏微分した値が得られる。これらの偏微分値を得るための演算式は数５から数９で与えられる。
【００６０】

(10)

特許３０８４９５１

19

20

【数５】

【００６１】【００６２】【数７】＊ ＊ 【数６】

10【００６３】【数８】【００６４】【数９】20【００６５】即ち、数５から数７の偏微分演算とＳｊ−
Ｏｊに基づいて修正量Δｘ，Δｙ，Δｚを得ることができる。修正量Δｘ，Δｙ，Δｚが得られれば、頂点抽出部４２ｉにより、算出された修正量Δｘ，Δｙ，Δｚの符号に基づいて、推定されている電流素片のベクトル成
分の再配分対象となる隣接頂点を抽出する。具体的には、図１１に破線で示すように修正ベクトルが得られるのであるから、二重丸で示される頂点が再配分対象となる隣接頂点として抽出される。そして、正規化部４３ｉ 30
により修正ベクトルを正規化し、正規化された修正ベク
トルの各成分の絶対値に基づいて再配分比率算出部４４ｉにより再配分比率を算出し、再配分部４５ｉにより、推定されている電流素片のベクトル成分および頂点抽出
部４２ｉにより抽出された頂点に関して推定されている
電流素片のベクトル成分に推定値の再配分を行なう。この場合において、各電流素片のベクトル成分は複数の修正ベクトルに基づく再配分が行なわれる可能性があるが、これらを単純に加算すると再配分前の値と再配分後の値とが大幅に異なり、推定の安定性を損なう危険性が 40あるので、推定値更新部４６ｉにおいて、該当する現在の推定値（他の推定値に配分して残りの割合）に対して再配分値の平均（複数の修正ベクトルに基づく再配分比率の平均）を加算することにより新たなベクトル成分推定値を得るようにしている。
【００６６】以下、全ての頂点について修正量Δｘ，Δｙ，Δｚが十分に小さくなるまで反復制御部４７によりビオ・サバール演算ユニット１１ｉ、座標値収集ユニッ
ト４０ａ、修正量算出部４１ｉ、頂点抽出部４２ｉ、正規化部４３ｉ、再配分比率算出部４４ｉ、再配分部４５ 50

ｉおよび推定値更新部４６ｉの処理を反復させることにより高精度の電流素片の推定を達成できる。尚、反復制
御部４７により推定処理を反復させる場合において、各小立方体の各頂点の３次元座標値は固定されているのであるから、上記修正量Δｘ，Δｙ，Δｚが直接数５、数６、数７に反映されることはないが、電流素片のベクト
ル成分は再配分処理により変更されるのであるから、ベクトル成分を通して間接的に反映されることになる。即ち、数５、数６、数７、数８および数９におけるＸｊｉ，Ｙｊｉ，Ｚｊｉ，Ｒｊｉは推定開始当初から推定終了まで変化せず、ｐｘｉ，ｐｙｉのみが変化し、最終的に高精度に推定されたｐｘｉ，ｐｙｉを得ることができる。この結果、例えば図１２（Ａ）に示すように、互にキャンセルし合う電流素片が互に近接する位置に割り当
てられていても、再配分処理を行なうことにより図１２（Ｂ）に示すように推定が行なわれ、電流が存在していない位置に互にキャンセルし合う電流が残留するという不都合を確実に解消できる。尚、図１２中、黒丸が処理
対象頂点を、白丸が隣接する頂点を示している。
【００６７】図１４から図１８は１チャネルのマグネト
メータを用いて３cm間隔、６×６点で計測された心磁図
を心電図を用いて同期加算して得た磁界分布図に基づいてＱ波近傍時刻の磁場源推定結果を示す図であり、医学的に確認されている心臓の体積電流に高精度に近似できる推定結果が得られていることが分る。尚、図１３はＭ
ＲＩ画像に基づいて得られた深さが６cmの頂点情報であり、図１４から図１８はそれぞれＲ波出現前３５msec、
３３msec、３０msec、２５msec、２０msecの時刻に対応する推定結果であり、大きい黒丸が心臓の壁を示しているとともに、黒丸を起点とする線分で電流素片を示している。
【００６８】尚、この実施例において、各ビオ・サバー
ル演算ユニットに代えて階層型パーセプトロンを採用することが可能である。
【００６９】【実施例７】図１９はこの発明の生体磁場測定結果表示
方法の一実施例を説明するフローチャートであり、ステップＳＰ１において、生体断層像に対応する平面内における、所定の時点における電流素片の位置および大きさ（興奮電流の絶対値）が得られるまで待ち、ステップＳ
Ｐ２において複数の電流素片のそれぞれの大きさに対応する直径の円を割り当て、ステップＳＰ３において生体

(11)
21
断層像を表示し、ステップＳＰ４において、割り当てられた円を対応箇所に生体断層像とオーバーラップした状
態で表示し、ステップＳＰ５において興奮電流表示の終了が指示されたか否かを判別し、終了が指示されていな
い場合には、ステップＳＰ６において異なる時点における電流素片の位置および大きさが得られるまで待ってから再びステップＳＰ２の処理を行なう。逆に、ステップ
ＳＰ５において終了が指示された場合にはそのまま一連の処理を終了する。但し、ステップＳＰ２における円の割り当てについては、電流素片の大きさに正確に対応する直径の円を割り当てるようにしてもよいが、電流素片の大きさを複数段階に区分し、各区分に対応して段階的に変化するように円の直径を割り当ててもよい。
【００７０】【実施例８】図２０はこの発明の生体磁場測定測結果表
示方法の他の実施例を説明するフローチャートであり、ステップＳＰ１において、生体断層像に対応する平面内における臓器等の存在箇所が検出されるまで待ち、ステップＳＰ２において臓器等の存在箇所に対応して所定の直径の通常の円を割り当て、、ステップＳＰ３において
上記平面内における、所定の時点における電流素片の位置および大きさが得られるまで待ち、ステップＳＰ４において、該当する円に代えて複数の電流素片のそれぞれの大きさに対応する直径の強調表示円を割り当て、ステップＳＰ５において生体断層像を表示し、ステップＳＰ
６において、割り当てられた円を対応箇所に生体断層像とオーバーラップした状態で表示し、ステップＳＰ７において興奮電流表示の終了が指示されたか否かを判別し、終了が指示されていない場合には、ステップＳＰ８において異なる時点における電流素片の位置および大きさが得られるまで待ってから再びステップＳＰ４の処理を行なう。逆に、ステップＳＰ７において終了が指示された場合にはそのまま一連の処理を終了する。
【００７１】図２１はこの実施例による興奮電流の表示
例を示す図であり、著しく高い視認性を発揮できていることが分る。また、興奮電流が存在していない箇所においても臓器等に対応して通常の円が表示されているとともに、円同士の間にはかなりの間隙が存在し、この間隙を通して断層像を簡単に確認できる。特に、時系列的に興奮電流を表示する場合には、例えば最大の興奮電流が変化する状態を簡単に把握できることになる。
【００７２】【実施例９】図２２はこの発明の生体磁場測定結果表示
方法のさらに他の実施例を説明するフローチャートであり、ステップＳＰ１において所定の時点における電流素
片の位置および大きさが得られるまで待ち、ステップＳ
Ｐ２において複数の電流素片のそれぞれの大きさに対応する直径の球を割り当て、ステップＳＰ３において生体構造参照像を３次元的に半透明表示し、ステップＳＰ４において、割り当てられた球を対応箇所に生体構造参照

10

20

30

40

50

特許３０８４９５１
22
像とオーバーラップした状態で表示し、ステップＳＰ５において興奮電流表示の終了が指示されたか否かを判別し、終了が指示されていない場合には、ステップＳＰ６において異なる時点における電流素片の位置および大きさが得られるまで待ってから再びステップＳＰ２の処理を行なう。逆に、ステップＳＰ５において終了が指示された場合にはそのまま一連の処理を終了する。
【００７３】したがって、この実施例の場合には、３次元的な生体構造参照像に対応付けて電流素片の大きさに対応する直径の球を表示でき、時系列的に球の表示を行なうことにより興奮電流の３次元的な伝播方向を簡単に把握できることになる。
【００７４】【実施例１０】図２３はこの発明の生体磁場測定結果表
示装置の一実施例を示すブロック図であり、電流素片解
析結果を保持する解析結果保持部６００と、電流素片の大きさに対応する直径の円を割り当てる円割り当て部７
００と、生体断層像を表示する断層像表示部８００と、断層像とオーバーラップさせた状態で割り当てられた直径の円を表示する円表示部９００とを有している。
【００７５】したがって、生体断層像の表示とオーバーラップさせて電流素片の大きさに対応する円を表示でき、興奮電流の視認性を著しく高めることができる。そして、時系列的に解析結果に基づいて円のオーバーラッ
プ表示を行なうことにより興奮電流の伝播方向を簡単に把握できる。尚、上記実施例７〜１０において、例えば、生体像の視認性の多少の低下を許容できる場合には、円に代えて正方形、球に代えて立方体を表示することが可能であるほか、この発明の要旨を変更しない範囲内において種々の設計変更を施すことが可能である。
【００７６】【発明の効果】以上のように請求項１の発明は、所定間隔毎の断層像しか得られていなくても、補間演算を行なうことにより全空間に対応して立体像データを得、しかも区分された小立体領域に含まれる立体像データの数に基づいて小立体領域を２値化することにより単純化された立体像を得ることができ、ひいては、可視的に表示した場合に、臓器等に相当する箇所、空間に相当する箇所
を簡単に認識できるという特有の効果を奏する。
【００７７】請求項２の発明も、所定間隔毎の断層像しか得られていなくても、補間演算を行なうことにより全空間に対応して立体像データを得、しかも区分された小立体領域に含まれる立体像データの数に基づいて小立体領域を２値化することにより単純化された立体像を得ることができ、ひいては、可視的に表示した場合に、臓器等に相当する箇所、空間に相当する箇所を簡単に認識できるという特有の効果を奏する。
【００７８】請求項３の発明は、複数のビオ・サバールの法則の演算結果の累積加算値が磁場計測値と高精度に近似されるまで補正を行なうことにより、各小立体領域

(12)
23
の頂点における磁場源情報を得ることができ、全体として処理負荷を大幅に低減できるという特有の効果を奏する。請求項４の発明は、互に近接する状態で割り当てら
れた、大きさが互に等しくかつ向きが互に逆の磁場源ベクトルが存在していても、これらを確実に補正でき、一層高精度に推定された磁場源情報を得ることができるとともに、隣接する頂点間で再配分を行なうことにより推定を行なうのであるから、推定対象となる磁場源情報の数に見合って従来必要とされていた測定点数よりも少ない測定点に基づいて高精度の磁場源推定を達成できるという特有の効果を奏する。
【００７９】請求項５の発明も、複数のビオ・サバールの法則の演算結果の累積加算値が磁場計測値と高精度に近似されるまで補正を行なうことにより、各小立体領域の頂点における磁場源情報を得ることができ、全体として処理負荷を大幅に低減できるという特有の効果を奏する。請求項６の発明も、互に近接する状態で割り当てら
れた、大きさが互に等しくかつ向きが互に逆の磁場源ベクトルが存在していても、これらを確実に補正でき、一層高精度に推定された磁場源情報を得ることができるとともに、隣接する頂点間で再配分を行なうことにより推定を行なうのであるから、推定対象となる磁場源情報の数に見合って従来必要とされていた測定点数よりも少ない測定点に基づいて高精度の磁場源推定を達成できるという特有の効果を奏する。
【００８０】請求項７の発明は、生体構造参照像を表示
するとともに、生体構造参照像の該当位置に設定されたサイズの単位図形を表示するのであるから、矢印の方向と長さで物理量計測結果を表示する場合と比較して、例
えば、最大の生体物理量計測結果の分布状態を簡単に認識できるとともに、生体物理量計測結果を動的に表示することにより、例えば、最大の生体物理量計測結果の分
布状態の時間的変動を簡単に認識できるという特有の効果を奏する。
【００８１】請求項８の発明は、請求項７の効果に加え、単位図形が円であるから、単位図形同士の隙間をかなり大きくでき、背景となる生体構造参照像の把握が容易になるという特有の効果を奏する。請求項９の発明
は、生体構造参照像を表示するとともに、生体構造参照
像の該当位置に設定されたサイズの単位図形を表示するのであるから、矢印の方向と長さで物理量計測結果を表示する場合と比較して、例えば、最大の生体物理量計測
結果の分布状態を簡単に認識できるとともに、生体物理量計測結果を動的に表示することにより、例えば、最大の生体物理量計測結果の分布状態の時間的変動を簡単に認識できるという特有の効果を奏する。
【００８２】請求項１０の発明は、請求項９の効果に加え、単位図形が円であるから、単位図形同士の隙間をかなり大きくでき、背景となる生体構造参照像の把握が容易になるという特有の効果を奏する。

10

20

30

40

50

特許３０８４９５１
24【図面の簡単な説明】【００８３】【図１】この発明の断層像処理方法の一実施例を説明す
るフローチャートである。
【００８４】【図２】断層像の一例を示す概略図である。
【００８５】【図３】表示立方体と非表示立方体との割り当て状態を概略的に示す平面図である。
【００８６】【図４】図３に対応する斜視図である。
【００８７】【図５】この発明の断層像処理装置の一実施例を示すブロック図である。
【００８８】【図６】この発明の生体磁場測定装置の一実施例を示す
ブロック図である。
【００８９】【図７】この発明の生体磁場測定装置の他の実施例を示
すブロック図である。
【００９０】【図８】１つの階層型パーセプトロンに対応する部分の構成を詳細に示す概略図である。
【００９１】【図９】この発明の生体磁場測定方法のさらに他の実施例を示すブロック図である。
【００９２】【図１０】この発明の磁場源測定装置のさらに他の実施例を示すブロック図である。
【００９３】【図１１】推定されている電流素片のベクトル成分の再配分対象となる隣接頂点を抽出する処理を説明する概略図である。
【００９４】【図１２】再配分による電流素片の推定処理を説明する
概略図である。
【００９５】【図１３】ＭＲＩ画像に基づいて得られた深さが６cmの頂点情報を示す図である。
【００９６】【図１４】Ｒ波出現前３５msecに対応する推定結果を示す図である。
【００９７】【図１５】Ｒ波出現前３３msecに対応する推定結果を示す図である。
【００９８】【図１６】Ｒ波出現前３０msecに対応する推定結果を示す図である。
【００９９】【図１７】Ｒ波出現前２５msecに対応する推定結果を示

(13)

特許３０８４９５１2526す図である。
＊ 【０１０６】【０１００】【符号の説明】【図１８】Ｒ波出現前２０msecに対応する推定結果を示
３ 断層像補間部
７ 立方体割り当て部
８ 画す図である。
素数計数部【０１０１】
９ 判別部
１０ 表示立方体割当部
２０，２０【図１９】この発明の生体磁場測定結果表示方法の一実
０ シグマ・ユニット
施例を説明するフローチャートである。
３０，３００ 誤差演算器
４１ｉ 修正量算出部【０１０２】
４２ｉ 頂点抽出部【図２０】この発明の生体磁場測定結果表示方法の他の４３ｉ 正規化部
４４ｉ 再配分比率算出部
４
実施例を説明するフローチャートである。
10 ５ｉ 再配分部
【０１０３】
４６ｉ 推定値更新部
４０，５００ 情報収集ユニ【図２１】図２０の実施例による興奮電流の表示例を示
ット
す図である。
１０１，１０２，・・・，１０ｍ 階層型パーセプトロ【０１０４】
ン【図２２】この発明の生体磁場測定結果表示方法のさら
１１１，１１２，・・・，１１ｍ ビオ・サバール演算に他の実施例を説明するフローチャートである。
ユニット【０１０５】
１１１ａ，１１２ａ，・・・，１１ｍａ，１２１，１２【図２３】この発明の生体磁場測定結果表示装置の一実
２，・・・，１２ｍ 補正部
施例としての興奮電流表示装置を示すブロック図であ７００ 円割り当て部
８００ 断層像表示部る。
＊20 ９００ 円表示部
【図２】

【図１１】

【図１４】

【図１２】

【図３】

【図１３】

【図４】

【図２１】

(14)

【図１】

特許３０８４９５１

(15)

【図５】

特許３０８４９５１

(16)特許３０８４９５１

【図６】

【図１５】

【図１６】

【図１７】

(17)

【図７】

特許３０８４９５１

(18)

【図８】

特許３０８４９５１

(19)

【図９】

特許３０８４９５１

(20)

【図１０】

特許３０８４９５１

(21)

【図１８】

【図１９】

特許３０８４９５１

(22)

【図２０】

特許３０８４９５１

(23)

【図２２】

【図２３】

特許３０８４９５１

(24)フロントページの続き

7(58)調査した分野(Int.Cl. ，ＤＢ名)
G06T 1/00
A61B 5/055
A61B 6/03

特許３０８４９５１