(57)【特許請求の範囲】
【請求項１】
方法であって、
第１の位置に置かれた第１の磁界放射器及び第２の位置に置かれた第２の磁界放射器から、ある領域において磁界を発生させることと、
前記領域内で、多数の頂点を有するある体積を設定することと、
前記多数の頂点において、前記磁界のそれぞれの値を測定することと、
前記磁界のそれぞれの値に応じて、前記多数の頂点のそれぞれに対して、前記第１の磁界放射器に対応する有効双極子モーメントとして第１の双極子モーメントを算出し、前記第２の磁界放射器に対応する有効双極子モーメントとして第２の双極子モーメントを算出し、
前記多数の頂点のそれぞれにおいて算出された前記第１及び第２の双極子モーメント全てに基づき、前記体積内の任意の点における前記第１及び第２の磁界放射器に対応した双極子モーメントを算出し、
前記体積内の任意の点における前記第１及び第２の磁界放射器に対応した双極子モーメントに基づき、前記任意の点における磁界の値を計算することと、
を含む、方法。
【請求項２】
前記任意の点における磁界の値を計算することに続いて、前記磁界の前記値を測定するように構成されたプローブを前記領域に挿入することと、前記測定された値に応じて、前記領域内の前記プローブの場所を決定することとを含む、請求項１に記載の方法。
【請求項３】
前記体積が、８つの頂点を有する立方体である、請求項１に記載の方法。
【請求項４】
前記第１の磁界放射器及び第２の磁界放射器が逆二乗則に従う極を有する単純な双極子として動作する、請求項１に記載の方法。
【請求項５】
前記第１の位置及び前記第２の位置によって定められる座標系の原点からの前記多数の頂点までの距離に応じて、前記多数の頂点のそれぞれに対して、前記第１及び第２の双極子モーメントを算出する、請求項１に記載の方法。
【請求項６】
前記体積内のある点についての前記磁界の前記値が、前記多数の頂点それぞれに対して算出された前記第１の双極子モーメントの第１の平均値と、前記多数の頂点それぞれに対して算出された前記第２の双極子モーメントの第２の平均値と、に基づいて計算される、請求項１に記載の方法。
【請求項７】
前記第１の平均値及び第２の平均値が、前記体積内の任意の点の場所に基づいて計算された、それぞれの線形加重平均値である、請求項６に記載の方法。
【請求項８】
前記第１の磁界放射器及び前記第２の磁界放射器が、それぞれ、第１の交番磁界を第１の周波数で送信し、第２の交番磁界を前記第１の周波数とは異なる第２の周波数で送信する、請求項１に記載の方法。
【請求項９】
装置であって、
第１の位置に置かれた第１の磁界放射器及び第２の位置に置かれた第２の磁界放射器であって、ある領域において磁界を発生させるように構成された、第１の磁界放射器及び第２の磁界放射器と、
プロセッサであって、
前記領域内で、多数の頂点を有するある体積を設定し、
前記多数の頂点において、前記磁界のそれぞれの値を測定し、
前記磁界のそれぞれの値に応じて、前記多数の頂点のそれぞれに対して、前記第１の磁界放射器に対応する有効双極子モーメントとして第１の双極子モーメントを算出し、前記第２の磁界放射器に対応する有効双極子モーメントとして第２の双極子モーメントを算出し、
前記多数の頂点のそれぞれにおいて算出された前記第１及び第２の双極子モーメント全てに基づき、前記体積内の任意の点における前記第１及び第２の磁界放射器に対応した双極子モーメントを算出し、
前記体積内の任意の点における前記第１及び第２の磁界放射器に対応した双極子モーメントに基づき、前記任意の点における磁界の値を計算するように構成された、プロセッサと、
を備える、装置。
【請求項１０】
前記任意の点における磁界の値を計算することに続いて、前記領域に挿入されるプローブを備え、前記プローブが、前記磁界の前記値を測定するように構成され、前記プロセッサが、前記測定された値に応じて、前記領域内の前記プローブの場所を決定するように構成される、請求項９に記載の装置。
【請求項１１】
前記体積が、８つの頂点を有する立方体である、請求項９に記載の装置。
【請求項１２】
前記第１の磁界放射器及び第２の磁界放射器が、逆二乗則に従う極を有する単純な双極子として動作する、請求項９に記載の装置。
【請求項１３】
前記プロセッサが、前記第１の位置及び前記第２の位置によって定められる座標系の原点からの前記多数の頂点までの距離に応じて、前記多数の頂点のそれぞれに対して、前記第１及び第２の双極子モーメントを算出する、請求項９に記載の装置。
【請求項１４】
前記体積内のある点についての前記磁界の前記値が、前記多数の頂点それぞれに対して算出された前記第１の双極子モーメントの第１の平均値と、前記多数の頂点それぞれに対して算出された前記第２の双極子モーメントの第２の平均値と、に基づいて計算される、請求項９に記載の装置。
【請求項１５】
前記第１の平均値及び第２の平均値が、前記体積内の任意の点の場所に基づいて計算された、それぞれの線形加重平均値である、請求項１４に記載の装置。
【請求項１６】
前記第１の磁界放射器及び前記第２の磁界放射器が、それぞれ、第１の交番磁界を第１の周波数で送信し、第２の交番磁界を前記第１の周波数とは異なる第２の周波数で送信する、請求項９に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、全般的には磁界発生に関し、詳細には、発生した磁界のモデリングに関する。
【背景技術】
【０００２】
広範な医療処置は、身体内に、センサ、チューブ、カテーテル、分配装置、及び移植片のような対象物を設置することを伴う。これらの医療処置において、医師が対象物及びその周囲を可視化することを補助するために、リアルタイムイメージング法がしばしば用いられる。いくつかの方法では、磁界を利用して対象物を追跡する。しかし、磁界での擾乱が追跡において誤差を発生させる可能性がある。
【０００３】
参照により本特許出願に組み込まれる文書は、いずれかの用語が、それらの組み込まれた文書内で、本明細書で明示的又は暗示的に行われる定義と相反するように定められる場合を除き、本出願の一体部分と見なされるべきであり、本明細書における定義のみが考慮されるべきである。
【発明の概要】
【課題を解決するための手段】
【０００４】
本発明のある実施形態は、方法を提供するものであり、その方法は、
第１の位置に置かれた第１の磁界放射器及び第２の位置に置かれた第２の磁界放射器から、ある領域において磁界を発生させることと、
この領域内で、多数の頂点を有するある体積の輪郭を描くことと、
多数の頂点において、磁界のそれぞれの値を測定することと、
それぞれの値に応じて、それぞれの第１の双極子モーメントを第１の磁界放射器に割り当て、それぞれの第２の双極子モーメントを第２の磁界放射器に割り当てることと、
第１の双極子モーメント及び第２の双極子モーメントに基づいて、その体積内の磁界の値を計算することと、を含む。
【０００５】
本方法は、典型的には、値を計算することに続いて、前記磁界の前記値を測定するよう構成されたプローブをこの領域に挿入することと、前記測定された値に応じて、この領域内のプローブの場所を決定することとを含む。
【０００６】
開示された実施形態では、この体積は、８つの頂点を有する立方体である。
【０００７】
更なる開示された実施形態では、磁界の値を計算することは、逆二乗則に従う極を有する単純な双極子として動作する第１及び第２の磁界放射器を含む。
【０００８】
本方法は、第１の位置及び第２の位置によって定められる座標系の原点からの多数の頂点の変位に応じて、それぞれの第１の双極子モーメントを第１の磁界放射器に割り当てることと、それぞれの第２の双極子モーメントを第２の磁界放射器に割り当てることと、を更に含んでもよい。
【０００９】
代替実施形態では、この体積内のある点についての磁界の値は、多数の頂点に対して割り当てられた第１の双極子モーメントの第１の平均値と、多数の頂点に対して割り当てられた第２の双極子モーメントの第２の平均値と、に基づいて計算される。典型的には、第１及び第２の平均値は、この体積内の点の場所に基づいて計算された、それぞれの線形加重平均値である。
【００１０】
更なる代替実施形態では、第１の磁界放射器及び第２の磁界放射器は、それぞれ、第１の交番磁界を第１の周波数で送信し、第２の交番磁界を第１の周波数とは異なる第２の周波数で送信する。
【００１１】
本発明のある実施形態により、装置が更に提供され、この装置は、
第１の位置に置かれた第１の磁界放射器及び第２の位置に置かれた第２の磁界放射器であって、ある領域において磁界を発生させるように構成された、第１の磁界放射器及び第２の磁界放射器と、
プロセッサであって、
この領域内で、多数の頂点を有するある体積の輪郭を描き、
この多数の頂点において、磁界のそれぞれの値を測定し、
それぞれの値に応じて、それぞれの第１の双極子モーメントを第１の磁界放射器に割り当て、かつそれぞれの第２の双極子モーメントを第２の磁界放射器に割り当て、
第１の双極子モーメント及び第２の双極子モーメントに基づいて、この体積内の磁界の値を計算するように構成された、プロセッサと、を含む。
【００１２】
以下の本開示の発明を実施するための形態を、図面と併せて読むことによって、本開示のより深い理解が得られるであろう。
【図面の簡単な説明】
【００１３】
【図１】本発明のある実施形態による、磁界モデリングシステムの概略描写図である。
【図２】本発明のある実施形態による、図１のシステムによって形成されるモデルの製作を説明する工程の流れ図である。
【図３】本発明のある実施形態による、工程を示す概略図である。
【発明を実施するための形態】
【００１４】
概説
磁界放射器によって送信された、ある領域における磁界をモデリングするための従来技術のシステムは、計算的に重い球面調和関数などの数学的手順を用いる。この数学的手順は、単純な双極子放射器としては挙動しない実践上の放射器を説明するために必要であり、この手順は、放射器についての単純な双極子モデルからの偏差を補償する。
【００１５】
本発明の実施形態は、異なるアプローチを採用し、放射器が単純な双極子放射器として挙動するモデルを想定する。しかしながら、所定の放射器に対して想定される双極子モーメントは、この領域内の場所の関数であると想定される。このモデルを用いることは、モデリングの従来技術のシステムと比べて、このモデルによって予測される磁界の精度の低下を伴うことなく、簡略化された計算をもたらす。
【００１６】
したがって、多数の頂点を有する領域内のある体積について、複数の放射器からの磁界が、それぞれの頂点で測定される。所定の放射器について、放射器に関するそれぞれの双極子モーメントが、頂点のそれぞれに割り当てられる。この体積内のある点における磁界の値が、この体積のそれぞれの頂点のそれぞれの双極子モーメントに基づいて計算される。計算は、多数の頂点に割り当てられた双極子モーメントの平均値、典型的には線形加重平均値を見出すことによって、行われる。
【００１７】
本発明者は、この領域を満たす連続的体積に対して計算を適用することが、領域内の磁界をモデリングするための効率的、迅速、かつ正確な方法を提供することを見出した。更に、この計算は、必ずしも体積境界において可微分ではなく、全領域にわたって連続するモデル化された磁界をもたらす。
【００１８】
システムの説明
ここで、本発明のある実施形態による、磁界モデリングシステム２０の概略描写図である、図１を参照する。システム２０は、領域２４に送信する、複数の概ね類似した交番磁界放射器２２によって発生した磁界のモデルを製作する。領域２４は、図において楕円形により概略的に表現されている。
【００１９】
典型的には、一旦磁界がシステム２０によってモデル化されると、患者が磁界内に配置される。カテーテルプローブ２６が、医療専門家によって患者の器官内に挿入され、プローブにおける磁界が測定される。（明確化のために、プローブ２６は、点線で図１に示されている。）測定された値が、モデルを導出した磁界の値と比較され、この比較が患者の器官内のプローブの位置及び向きを決定する。患者は、典型的には、医療専門家によって実施される処置の間、（患者の器官が領域２４内にある状態で）手術台３０上に横たわる。明確化及び簡略化のために、図１には患者及び医療専門家の双方共示されていない。
【００２０】
本発明の実施形態においては、ｐ個の放射器２２があり、このｐは、１を超える整数である。放射器２２は、本明細書で、Ｒ１、Ｒ２、．．．Ｒｎ、．．．Ｒｐとしても称され、ここで１＜ｎ≦ｐであり、ｎは整数である。放射器は、手術台３０に対して固定され、固定された放射器は、直交するｘｙｚ軸のセットを有する放射器の座標系を定める。上述したように、放射器２２は、領域に交番磁界を送信し、放射器Ｒ１、Ｒ２、．．．Ｒｎ、．．．Ｒｐは、それぞれの周波数ｆ１、ｆ２、．．．ｆｎ、．．．ｆｐにおいて放射すると想定される。
【００２１】
一実施形態では、この放射器は、位置パッド４０（これが今度は手術台に固定される）に固定された３つの放射器セット３４、３６、及び３８に分類されると想定される。各セットは、３つの３軸コイルを備え、３軸コイルは互いに直交し、３つの磁界放射器として動作する。したがって、セット３４は、放射器Ｒ１、Ｒ２、Ｒ３を備え、セット３６は、放射器Ｒ４、Ｒ５、Ｒ６を備え、セット３８は、放射器Ｒ７、Ｒ８、Ｒ９を備える。その開示が参照により本明細書に組み込まれる米国特許第６，４８４，１１８号（Ｇｏｖａｒｉ）は、磁界放射器として用いられる、このようなコイルの配置を記載している。例として、放射器座標系のｘｙｚ軸は、その原点を位置パッド４０内に有すると想定され、ｚ軸は、位置パッドに対して垂直であると想定される。（明確化のために、図１は、位置パッドから離隔されたｘｙｚ軸を示す。）
【００２２】
しかしながら、上記で想定された９つの放射器の配置は一例であり、それゆえ、本発明の他の実施形態が、放射器の他の個数及び／又は配置を有することが理解される。上記で参照した米国特許第６，４８４，１１８号は、放射器の他の配置を記載している。磁気共鳴画像診断（ＭＲＩ）施設においてなどの放射器からの磁界が歪む動作環境においては、開示された実施形態は、１５個の放射器を用いる。
【００２３】
制御ユニット５０は、システム２０を操作する。制御ユニット５０は、プロセッサ５２、典型的には、適切な信号処理回路を備えるコンピュータを含む。ユニット５０は、プロセッサ５２が放射器２２を駆動するために用いる、磁界発生器モジュール６０を備える制御ユニットはまた、放射器２２によって発生した磁界内に配置された１つ又は２つ以上のプローブから信号を受信するように、及び典型的にはプロセッサ５２と一緒に、１つ又は２つ以上のプローブにおいて磁界を評価するように構成される、磁界測定モジュール６４を備える。プロセッサは、システム２０を操作するためのデータが記憶されている、揮発性データ記憶装置及び非揮発性データ記憶装置の両方を典型的に含む、メモリ５４を使用する。プロセッサは、典型的には、画像表示５６をシステム２０の操作者に提供するように連結される。
【００２４】
典型的には、プロセッサ５２は、本明細書に記載される機能を実行するようにソフトウェアでプログラムされる。そのソフトウェアは、例えば、ネットワークを介して電子形式でプロセッサにダウンロードされてもよく、あるいは、それに代わって、若しくはそれに加えて、磁気的、光学的、又は電子的メモリなどの非一時的な有形媒体上に提供及び／又は記憶されてもよい。
【００２５】
システム２０は、領域２４内に放射器２２によって送信された磁界を測定することが可能な装置も備える。本発明の一実施形態では、マッパー７０が、磁界を測定するために用いられ、マッパーは、プラスチックのシートなどの固体ベース上の既知の位置に固定するよう据え付けられたｑ個の磁界検出器７２の配列を備える。マッパー７０は、位置パッド４０に対して既知の所定の位置及び向きで手術台３０上に配置できるように構成されている。一実施形態では、検出器７２は、これらが領域２４を満たす立方体体積の頂点にあるように配置される。一実施形態では、立方体体積の縁部は２ｃｍの長さであり、それにより、頂点は２ｃｍで離隔される。しかしながら、縁部の任意の他の好都合な長さが用いられてもよい。
【００２６】
検出器７２は、Ｈａｌｌプローブ若しくは１つ又は２つ以上のコイルなどのような、交番磁界の大きさ及び方向を測定するための任意の好都合なセンサを備えてもよい。検出器からの信号は、典型的にはケーブル７４によってモジュール６４に送信されるが、無線送信などの任意の他の好都合な送信法が用いられてもよい。モジュール６４及び／又はプロセッサ５２は、この信号を用いて検出器における磁界を測定する。
【００２７】
マッパー７０に代わる磁界測定システムは、当業者には明らかであろう。例えば、検出器７２の配列の代わりに、検出器７２と実質的に同様な１つ又は２つ以上の検出器が、三次元で既知の距離まで平行移動するように構成されるジグ上に取り付けられてもよく、磁界の測定は、１つ又は２つ以上の検出器が領域２４を横断すると同時に行われる。当業者には明白である、他のこのようなシステムは、本発明の範囲内に含まれることが想定される。
【００２８】
以下により詳細に説明されるように、マッパー７０は、システム２０の校正段階で用いられてもよい。システム２０がその操作段階にあるときには、マッパー及びそのケーブルは取り外される。
【００２９】
図２は、本発明のある実施形態による、システム２０によって形成されるモデルの製作を説明する工程の流れ図であり、図３は、本発明のある実施形態による、工程を示す概略図である。最初の工程２００において、放射器２２は、正しい位置に固定され、放射器Ｒ１、Ｒ２、．．．Ｒｎ，．．．Ｒｐのそれぞれの場所が測定される。図３は、それぞれの位置ベクトルｒ１、ｒ２を有する、２つの放射器Ｒ１、Ｒ２の場所を示す。
【００３０】
組立工程２０２では、領域２４内の既知の場所において、放射器２２によって発生した磁界を測定するためのシステムが、設定される。明確化のために、本明細書の説明において、システムの設定は、マッパー７０を、立方体の頂点で配置された検出器７２と共に備えると想定される。図３は、８つの立方体頂点Ｖ１、Ｖ２、．．．Ｖ８（これらは、それぞれの位置ベクトルｃ１、ｃ２、．．．ｃ８を有すると想定される）を有する１つの立方体を示す。マッパー７０のｑ個の検出器について、それぞれが位置ベクトルｃ１、ｃ２、．．．ｃｑを有するｑ個の頂点Ｖ１、Ｖ２、．．．Ｖｑが存在する。以下の説明では、領域２４内の一般的頂点は、位置ベクトルｃｍ（このｍは添え字である）を有することが想定され、頂点はまた、それらの位置ベクトルによって本明細書で称されてもよい。
【００３１】
測定工程２０４において、放射器２２が起動され、各放射器２２からのそれぞれの磁界の大きさ及び方向が、各検出器７２によって測定される。したがって、９個の放射器がある場合、各検出器は、放射器からの磁界に対して９個の大きさ及び９個の方向を測定する。一般に、ｑ個の検出器のそれぞれは、ｐ個の放射器からの磁界のｐ個の大きさ及びｐ個の方向を測定する。
【００３２】
計算工程２０６において、各頂点に対して、プロセッサ５２は、各放射器Ｒｎについてそれぞれの有効双極子モーメントを計算する。この計算は、各放射器が、双極子の磁極が従う逆二乗則から導出された双極子場方程式を有する単純な双極子であることを想定する。単純な双極子場方程式は、等式（１）によってもたらされ、
【数１】
0006621656.tif 000002
式中、ｍは、放射器の双極子モーメントであり、
ｒは、放射器からの点の変位であり、
並びにＨは、この点において放射器によって発生される磁界である。Ｈ、ｍ、及びｒは、ベクトルである。
【００３３】
周波数ｆｎで頂点ｃｍに対して放射する放射器Ｒｎについて、及び頂点が放射器からの変位ｒｎｍを有する場合、等式（１）は、次のように書き換えてもよい。
Ｈｎｍ＝ｆ（ｍｎｍ，ｒｎｍ）                          （２）
式中、ｍｎｍは、頂点ｃｍにおける放射器Ｒｎの有効双極子モーメントであり、
ｒｎｍは、放射器からの頂点ｃｍの変位であり、
Ｈｎｍは、この頂点において放射によって発生される磁界であり、
ｆは、等式（１）により得られた関数である。
【００３４】
等式（２）は、等式（３）として書き換えられてもよい。
ｍｎｍ＝ｇ（Ｈｎｍ，ｒｎｍ）                          （３）
式中、ｍｎｍ、Ｈｎｍ、ｒｎｍは、等式（２）に関して上記で定められた通りであり、
ｇは、関数ｆから導出された関数である。
【００３５】
したがって、工程２０６において、プロセッサは、各頂点ｃｍ、各放射器Ｒｎについてのそれぞれの有効双極子モーメントを計算するために、等式（３）を用いる。プロセッサは、流れ図の残りの工程で使用するために、双極子モーメントの値を記憶する。
【００３６】
工程２０６は、モデリングシステム２０の校正段階の完結工程であり、ここでは放射器Ｒｎからの磁界が、領域２４内に位置する校正頂点について決定されている。流れ図の次の工程は、工程２０６で見出された結果に基づき、領域２４内の任意の点における磁界についての表現を見出す補間プロセスを説明する。
【００３７】
特定工程２０８において、任意の点Ａ（ｘ，ｙ，ｚ）が位置する立方体の頂点が特定される。図３を参照すると、点Ａは、頂点Ｖ１、Ｖ２、．．．Ｖ８を有する立方体内にあると想定される。
【００３８】
以下の説明において、ｘｙｚ軸の座標の原点は、その初期原点から頂点Ｖ１における原点まで平行移動されているものと想定され、立方体縁部は、軸の方向を定めると想定される。平行移動は、既知の量によるものであることが理解されるであろう。加えて、立方体縁部は、１の公称長さを有すると想定される。この場合、点Ａの座標は、特性０≦ｘ，ｙ，ｚ≦１を有する。上記の想定は、以下の説明を簡略化するが、本発明の実施形態は、立方体の頂点を定めるために使用される軸に対して、又は立方体縁部の長さに対して制限を課すものではないことが理解されるであろう。（明確化のために、軸の頂点Ｖ１への平行移動を、図３には示していない。）
【００３９】
第１の双極子モーメント工程２１０において、プロセッサは、立方体の頂点のそれぞれに対して、選択された放射器Ｒｎの有効双極子モーメントを呼び戻す。この場合、双極子モーメントは、ｍｎ１、ｍｎ２、ｍｎ３、．．．ｍｎ８である。
【００４０】
次いでプロセッサは、点Ａを囲む頂点のそれぞれの有効双極子モーメントを、点Ａについての有効双極子モーメントとして計算する。本明細書に記載の開示された実施形態では、この平均値は、線形加重平均値であるが、更に単純な平均値を含む他の平均値が、本発明の範囲内に含まれることが理解されるであろう。
【００４１】
等式（４）は、線形加重平均値に応じる。
【数２】
0006621656.tif 000003
式中、
【数３】
0006621656.tif 000004
は、頂点Ｖ１、Ｖ２、．．．Ｖ８を有する立方体内の点Ａ（ｘ，ｙ，ｚ）についての放射器Ｒｎの有効双極子モーメントである。
【００４２】
プロセッサは、等式（４）を使用して、ｐ個の放射器Ｒ１、．．．、Ｒｐのそれぞれに対して点Ａ（ｘ，ｙ，ｚ）において使用されることになる有効双極子モーメントを計算する。
【００４３】
第２の双極子モーメント工程２１２において、プロセッサは、工程２１０において評価されたｐ個の有効双極子モーメントを用いて、Ａ点（ｘ，ｙ，ｚ）における放射器Ｒ１、．．．、Ｒｐから発生したｐ個のそれぞれの磁界を計算する。磁界のそれぞれは、等式（５）を用いて計算され、等式（５）は、等式（１）の応用版である。
【数４】
0006621656.tif 000005
式中、ｒＡは、座標の初期原点から点Ａに至るベクトルであり、
ＨＡｎ（ｒＡ）は、点Ａにおける放射器Ｒｎからのベクトル場であり、
【数５】
0006621656.tif 000006
は、等式（４）を参照して定められる。
【００４４】
最終の合計工程２１４において、プロセッサは、等式（６）に従って、等式（５）の適用から計算された全てのｐ個の磁界ＨＡ１（ｒＡ）、．．．ＨＡｎ（ｒＡ）、．．．ＨＡｐ（ｒＡ）をベクトル的に合計する。
Ｈ（ｒＡ）＝ＨＡ１（ｒＡ）＋．．＋ＨＡｎ（ｒＡ）＋．．＋ＨＡｐ（ｒＡ）  （６）
式中、Ｈ（ｒＡ）は、点Ａ（ｘ，ｙ，ｚ）における全ての放射器Ｒ１、．．．、Ｒｐからの全体の磁界である。
【００４５】
図１に戻ると、一旦図２の流れ図の手順が実施されると、磁界は領域２４の全てに対してモデル化される。患者は、検査される患者の部位が領域２４内にあるように配置され得、プローブ２６はこの領域において患者に挿入され得る。
【００４６】
プローブは、磁界測定モジュール６４と一緒に、ベクトル成分のそれぞれが異なる周波数で送信されるという特性を用いて、Ｈ（ｒＡ）（等式（６））中に含まれる磁界の全ての異なるベクトルを測定する。全ての異なる成分の値を認識することから、及び等式（５）から導出された成分の値のモデリングから、プロセッサ５２は、位置、すなわち、プローブ２６の場所及び向きを導出する。
【００４７】
上記説明は、検出器７２が立方体の頂点にあることを想定しているが、この検出器は、四面体などの他の実体積の頂点において配置されてもよいことを理解されるであろう。この実体積は、計算の容易さのために、典型的には同一タイプであるが、この制限は必ずしも必要ではなく、それにより、例えば、検出器は、立方体及び四面体の混合体の頂点にあると想定されてもよい。したがって、検出器７２にそれぞれする頂点は、領域２４を満たす実体積の任意の組み合わせの頂点であってもよい。このような組み合わせは、異なる体積についての加重平均値を導出するために、典型的には、１つ又は２つ以上の等式（４）とは異なる等式を必然的に伴うが、等式に必要な変更は、当業者には明らかであろう。
【００４８】
上述の実施形態は、例として引用したものであり、本発明は、上に具体的に示し、記載した内容に限定されるものではないということが理解されるであろう。むしろ、本発明の範囲は、上で説明した様々な特徴の組み合わせと部分的組み合わせの両方、並びにそれらの変形形態及び修正形態を含み、これらは、上述の説明を読めば当業者には想到されるであろうものであり、従来技術では開示されていないものである。
【００４９】
〔実施の態様〕
（１）  方法であって、
第１の位置に置かれた第１の磁界放射器及び第２の位置に置かれた第２の磁界放射器から、ある領域において磁界を発生させることと、
前記領域内で、多数の頂点を有するある体積の輪郭を描くことと、
前記多数の頂点において、前記磁界のそれぞれの値を測定することと、
前記それぞれの値に応じて、それぞれの第１の双極子モーメントを前記第１の磁界放射器に割り当て、それぞれの第２の双極子モーメントを前記第２の磁界放射器に割り当てることと、
前記第１の双極子モーメント及び前記第２の双極子モーメントに基づいて、前記体積内の前記磁界の値を計算することと、
を含む、方法。
（２）  前記値を計算することに続いて、前記磁界の前記値を測定するように構成されたプローブを前記領域に挿入することと、前記測定された値に応じて、前記領域内の前記プローブの場所を決定することとを含む、実施態様１に記載の方法。
（３）  前記体積が、８つの頂点を有する立方体である、実施態様１に記載の方法。
（４）  前記磁界の前記値を計算することが、逆二乗則に従う極を有する単純な双極子として動作する前記第１の磁界放射器及び第２の磁界放射器を含む、実施態様１に記載の方法。
（５）  前記第１の位置及び前記第２の位置によって定められる座標系の原点からの前記多数の頂点の変位に応じて、前記それぞれの第１の双極子モーメントを前記第１の磁界放射器に割り当て、前記それぞれの第２の双極子モーメントを前記第２の磁界放射器に割り当てることと、を含む、実施態様１に記載の方法。
【００５０】
（６）  前記体積内のある点についての前記磁界の前記値が、前記多数の頂点に対して割り当てられた前記第１の双極子モーメントの第１の平均値と、前記多数の頂点に対して割り当てられた前記第２の双極子モーメントの第２の平均値と、に基づいて計算される、実施態様１に記載の方法。
（７）  前記第１の平均値及び第２の平均値が、前記体積内の前記点の場所に基づいて計算された、それぞれの線形加重平均値である、実施態様６に記載の方法。
（８）  前記第１の磁界放射器及び前記第２の磁界放射器が、それぞれ、第１の交番磁界を第１の周波数で送信し、第２の交番磁界を前記第１の周波数とは異なる第２の周波数で送信する、実施態様１に記載の方法。
（９）  装置であって、
第１の位置に置かれた第１の磁界放射器及び第２の位置に置かれた第２の磁界放射器であって、ある領域において磁界を発生させるように構成された、第１の磁界放射器及び第２の磁界放射器と、
プロセッサであって、
前記領域内で、多数の頂点を有するある体積の輪郭を描き、
前記多数の頂点において、前記磁界のそれぞれの値を測定し、
前記それぞれの値に応じて、それぞれの第１の双極子モーメントを前記第１の磁界放射器に割り当て、それぞれの第２の双極子モーメントを前記第２の磁界放射器に割り当て、
前記第１の双極子モーメント及び前記第２の双極子モーメントに基づいて、前記体積内の前記磁界の値を計算するように構成された、プロセッサと、
を備える、装置。
（１０）  前記値を計算することに続いて、前記領域に挿入されるプローブを備え、前記プローブが、前記磁界の前記値を測定するように構成され、前記プロセッサが、前記測定された値に応じて、前記領域内の前記プローブの場所を決定するように構成される、実施態様９に記載の装置。
【００５１】
（１１）  前記体積が、８つの頂点を有する立方体である、実施態様９に記載の装置。
（１２）  前記第１の磁界放射器及び第２の磁界放射器が、逆二乗則に従う極を有する単純な双極子として動作する、実施態様９に記載の装置。
（１３）  前記プロセッサが、前記第１の位置及び前記第２の位置によって定められる座標系の原点からの前記多数の頂点の変位に応じて、前記それぞれの第１の双極子モーメントを前記第１の磁界放射器に割り当て、前記それぞれの第２の双極子モーメントを前記第２の磁界放射器に割り当てるように構成される、実施態様９に記載の装置。
（１４）  前記体積内のある点についての前記磁界の前記値が、前記多数の頂点に対して割り当てられた前記第１の双極子モーメントの第１の平均値と、前記多数の頂点に対して割り当てられた前記第２の双極子モーメントの第２の平均値と、に基づいて計算される、実施態様９に記載の装置。
（１５）  前記第１の平均値及び第２の平均値が、前記体積内の前記点の場所に基づいて計算された、それぞれの線形加重平均値である、実施態様１４に記載の装置。
【００５２】
（１６）  前記第１の磁界放射器及び前記第２の磁界放射器が、それぞれ、第１の交番磁界を第１の周波数で送信し、第２の交番磁界を前記第１の周波数とは異なる第２の周波数で送信する、実施態様９に記載の装置。
【図１】
0006621656.tif 000007
【図２】
0006621656.tif 000008
【図３】
0006621656.tif 000009