(57)【特許請求の範囲】
【請求項１】
被検体を励起パルスで励起し、発生する磁気共鳴信号を収集し、前記磁気共鳴信号に基づいて画像を生成する磁気共鳴イメージング装置において、
オペレータがショット数と位相エンコードマトリクスサイズとを指定するための入力手段と、
データの対称性を利用して略半分のエンコード数で撮影を行うハーフイメージング法において収集する０エンコード付近を占めるエコーデータの数、前記指定されたショット数及び前記指定された位相エンコードマトリクスサイズに基づいてトータルのエコー数を決定し、前記励起パルスによる励起後、前記決定したトータルのエコー数のエコーを収集するように前記励起パルスから始まる磁気共鳴信号を収集するための一連のシーケンスを構築し、当該シーケンスを前記ショット数の回数繰り返し実行する実行手段と、
を具備することを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記０エンコード付近を占めるエコーデータの数をContrast ETL、前記ショット数をShot of Image 、前記位相エンコードマトリクスサイズをPE Matrix Sizeとしたとき、前記トータルのエコー数Total ETL は、
Total ETL ＝Contrast ETL＋PE Matrix Size／(Shot of Image×２) 
により決定されることを特徴とする請求項１記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、磁気共鳴イメージング装置に関する。
【０００２】
【従来の技術】
磁気共鳴イメージングではそれが登場した当初から、撮影時間の短縮が重要な課題の１つとされており、現在までに、磁化励起後、リフォーカスＲＦパルス( 180 °パルス) を繰り返し印加しながら複数の磁気共鳴信号（エコー）を収集する一連のシーケンスを一定の繰り返し周期（繰り返し時間ＴＲ）で繰り返す高速スピンエコー法(fastSE)や、磁化励起後、リフォーカスＲＦパルスの代わりに傾斜磁場反転により複数のエコーを収集する一連のシーケンスを一定の繰り返し時間ＴＲで繰り返すエコープレナーイメージング法(EPI) 、さらにはデータの対称性を利用して半分のエンコード数で撮影を行うことで撮影時間の短縮を図るハーフイメージング法をfastSEに適用したfastASE 等様々な種類の高速イメージング法が開発されている。
【０００３】
ところでfastASE やEPI において、体動の影響を受ける腹部等の部位を息止めの状態で撮影したり、心臓の動きを捕らえるようなダイナミック撮影を行う場合には、撮影時間や時間分解能が重要である。
【０００４】
撮影時間Ｔは、ショット数（励起パルスから始まる一連のシーケンスの繰り返し回数）をＮとし、繰り返し時間をＴＲとして、
Ｔ＝Ｎ×ＴＲ                                              …（１）
で与えられる。オペレータは、ショット数と繰り返し時間ＴＲとから（１）式により撮影時間Ｔを容易に理解できる。
【０００５】
従来の磁気共鳴イメージング装置では、ユーザインタフェースとして、１回の励起で収集されるエコーの数（エコー数）ETL と、位相エンコードマトリクスサイズ（エンコード数）Matrixをオペレータが指定できるようになっており、指定されたETL とMatrixとに基づいて、以下の（２）式にしたがって装置側でショット数Ｎを自動的に計算し、計算したショット数の回数だけ一連のシーケンスを繰り返し実行するようになっている。
【０００６】
Ｎ＝int(Matrix／ETL )                                     …（２）
例えば、Matrix＝２５６、ETL ＝１０とすると、Ｎは２５で与えられる。
したがって、オペレータは（１）式で撮影時間を理解する前に、ETL とMatrixとから（２）式からショット数を求める必要があり、非常に面倒であり、直感的な理解が困難であるという問題があった。撮影時間を変更する場合も同様であり、所望の撮影時間とエコー数とから（１）式にしたがって位相エンコードマトリクスサイズMatrixを求め、これを指定する必要があった。このように従来の磁気共鳴イメージング装置では、ユーザインタフェースの制限上、オペレータが撮影時間を理解しにくいという問題があった。
【０００７】
また、従来の磁気共鳴イメージング装置では次のような問題もあった。図８（ａ）はfastSEにおける空打ちパルスシーケンスを示し、図８（ｂ）は実際に画像生成に用いるエコーを収集するためのパルスシーケンスを示している。通常、少なくとも１回目のシーケンスはいわゆる空打ちとして、エコーは収集しない。つまり、一連のシーケンスの繰り返しのうち、ｎ≧２として、ｎ回目以降のシーケンスで収集したエコーに基づいて画像を生成している。これは次の理由による。ｎ−１回目以前のシーケンスでは、磁化が完全に戻りきった、つまり繰り返し時間ＴＲが無限大の状態でエコーが収集される。一方、ｎ回目以降のシーケンスではエコーは磁化が完全に戻りきっていない状態で収集される。したがって、ｎ−１回目以前のシーケンスで収集したエコーの信号値はｎ回目以降のシーケンスで収集されるエコーの信号値と明かに相違することになる。この信号値の相違によるアーチファクトを抑制する目的で、ｎ回目以降のシーケンスで収集したエコーに基づいて画像を生成する。したがって、ｎ−１回目以前のシーケンスは、エコーを収集する必要がないので、図８（ａ）に示すようにリード傾斜磁場及び位相エンコード傾斜磁場の印加、エコーの収集を削除して簡略化されている。
【０００８】
また、磁場不均一による位相誤差に起因するアーチファクトを是正することを目的に、エコーの位相補正が不可欠とされている。磁場分布は被検体がコイルアセンブリに挿入されると変化するし、また被検体によっても個々に異なる。したがって、被検体がコイルアセンブリに挿入された状態で位相補正情報を取得する必要がある。このため、上述した一連のシーケンスを開始する前に、この位相補正情報を収集するためのシーケンスをプレスキャンとして実行することが要求されている。
これら空打ち及び位相補正情報を収集するためのシーケンスの実行に要する時間が撮影時間の短縮を疎外している。
【０００９】
【発明が解決しようとする課題】
本発明の第１の目的は、オペレータが撮影時間を理解しやすいような改良されたユーザインタフェースを有する磁気共鳴イメージング装置を提供することである。第２の目的は、撮影時間の短縮化を実現する磁気共鳴イメージング装置を提供することである。
【００１０】
【課題を解決するための手段】
本発明は、被検体を励起パルスで励起し、発生する磁気共鳴信号を収集し、前記磁気共鳴信号に基づいて画像を生成する磁気共鳴イメージング装置において、オペレータがショット数と位相エンコードマトリクスサイズとを指定するための入力手段と、データの対称性を利用して略半分のエンコード数で撮影を行うハーフイメージング法において収集する０エンコード付近を占めるエコーデータの数、前記指定されたショット数及び前記指定された位相エンコードマトリクスサイズに基づいてトータルのエコー数を決定し、前記励起パルスによる励起後、前記決定したトータルのエコー数のエコーを収集するように前記励起パルスから始まる磁気共鳴信号を収集するための一連のシーケンスを構築し、当該シーケンスを前記ショット数の回数繰り返し実行する実行手段とを具備する。
【００１２】
【作用】
請求項１に係る発明によれば、オペレータがショット数を指定すると、一連のシーケンスがこの指定されたショット数の回数だけ繰り返し実行される。したがって、オペレータは、ショット数Ｎと一連のシーケンスの繰り返し時間ＴＲとから、Ｎ×ＴＲにより撮影時間を容易に認識することができる。
【００１４】
【実施例】
以下、図面を参照して本発明による磁気共鳴イメージング装置の一実施例を説明する。
図１に本実施例による磁気共鳴イメージング装置の構成を示す。被検体Ｐを収容できるように円筒状の内部空間が形成されたコイルアセンブリ２０には、静磁場磁石１、傾斜磁場コイル２、ＲＦコイル３が装備される。常電導磁石又は超電導磁石である静磁場磁石１は、静磁場制御装置４から電流供給を受けて円筒内部に通常、Ｚ軸に沿って静磁場を形成可能に構成されている。傾斜磁場コイル２は、Ｘ、Ｙ、Ｚ各軸の傾斜磁場を作る３組のコイルから構成される。Ｘ軸傾斜磁場電源７，Ｙ軸傾斜磁場電源８，Ｚ軸傾斜磁場電源９それぞれは、傾斜磁場コイル２の対応する軸のコイルに電流を供給する。例えば、Ｘ、Ｙ、Ｚ各軸の傾斜磁場はそれぞれ、リード傾斜磁場read、位相エンコード傾斜磁場phase encode、スライス選択傾斜磁場slice selectとして用いられる。これら３方向の磁場強度が全て線形に変化する領域（撮影領域）内でエコーの収集（撮影）が可能である。被検体Ｐは寝台１３の天板に載置された状態で、天板のスライドに伴って撮影領域に挿入される。
【００１５】
ＲＦコイル３は、ＲＦパルスを被検体に送信し、被検体からの磁気共鳴信号を受信するためのコイルである。ここではＲＦコイル３は送受信で兼用されるとして説明するが、送信用と受信用とで別々に設けてもよい。送信器５は、高周波電流をパルス信号としてＲＦコイル３に供給する。このパルス長の変化で励起パルス（90°パルス）と再結像パルス（180 °パルス）とが区別される。受信器６は、プリアンプ、位相検波部、アナログ／ディジタル変換器とを有し、ＲＦコイル３を介して受信したエコー信号を増幅し、位相検波し、さらにアナログ／ディジタル変換器を介してディジタル信号として出力する。
【００１６】
コンピュータシステム１１は、受信器６からのディジタル信号を取り込み、これを位相補正処理及びＴ2 補正処理に供し、その後、２次元フーリエ変換処理（２ＤＦＴ）に供することにより磁気共鳴画像を生成する。この画像は表示部１２に表示される。また、コンピュータシステム１１は、エコーデータについて後述するtemplateデータに基づいて位相補正、Ｔ2 補正を行う。
【００１７】
コンピュータシステム１１には、キーボード等の入力部１３が接続される。オペレータによる入力部１３の操作により、ショット数、位相エンコードマトリクスサイズ（エンコード数）及びその他撮影に必要なデータが装置側に入力され、コンピュータシステム１１を介してシーケンサ１０に送られる。
【００１８】
シーケンサ１０は、ショット数及び位相エンコードマトリクスサイズに基づいてエコー数を決定し、ショット数、位相エンコードマトリクスサイズ及びエコー数に基づいてパルスシーケンス（パルス系列）を組み立てる。シーケンサ１０は、送信器５、受信器６、ＸＹＺ各軸の傾斜磁場電源７，８，９の各動作タイミングを制御して、組み立てたパルスシーケンスを実行する。
【００１９】
なお、ショット数とは、、fastSE、fastASE 、EPI における励起パルスから始まる一連のパルスシーケンスの繰り返し数のことをいう。また、エコー数とは、この一連のパルスシーケンスで収集されるエコーの数のことをいう。
【００２０】
次に本実施例の動作について説明する。図２にシーケンサ１０を主体としたパルスシーケンスの組み立てまでのフローチャートが示されている。オペレータによる入力部１３の操作により、少なくともショット数及び位相エンコードマトリクスサイズ（エンコード数）が装置側に入力され、コンピュータシステム１１を介してシーケンサ１０に送られる（ステップＳ１）。シーケンサ１０により、ショット数及び位相エンコードマトリクスサイズに基づいてエコー数が決定され（ステップＳ２）、ショット数、位相エンコードマトリクスサイズ及びエコー数に基づいてパルスシーケンスが組み立てられる（ステップＳ３）。そして、シーケンサ１０による送信器５、受信器６、ＸＹＺ各軸の傾斜磁場電源７，８，９の各動作タイミングの制御により、組み立てられたパルスシーケンスが所定の繰り返し時間ＴＲでショット数の回数だけ繰り返し実行される。
【００２１】
エコー数の決定は次のように行われる。ショット数をShot of Image 、位相エンコードマトリクスサイズをPE Matrix Sizeとすると、エコー数ETL は次の（３）式により決定される。
ETL ＝PE Matrix Size／Shot of Image                       …（３）
このようにショット数をオペレータが入力し、このショット数に応じて装置側でエコー数を決定することにより、オペレータは上述した（１）式により撮影時間を容易に把握することができ、また時間分解能を容易に調整できるようになる。また、ショット数と位相エンコードマトリクスサイズをオペレータが自由に設定できるので、撮影時間と分解能（コントラスト分解能、空間分解能）を直観的なパラメータとして調整できるようになる。
【００２２】
また、データの対称性を利用して半分のエンコード数で撮影を行うことで撮影時間の短縮を図るハーフイメージング法をfastSEに適用したfastASE であれば、次の（４）式により合計エコー数Total ETL が決定される。なお、画像コントラストを決めるにあたって支配的なのが０エンコード付近を占めるエコーデータであり、Contrast ETLとは、０エンコード付近に入るエコーまでのエコー数、つまり０エンコード付近を占めるエコーデータを何番目のエコーで収集するかをいう。
【００２３】
Total ETL ＝Contrast ETL＋PE Matrix Size／(Shot of Image×２) …（４）
この（４）式によりエコー数を決定することで、位相エンコードマトリクスサイズの総数に対する実際に収集するエンコード数の比率を変化させて、Contrast ETLを変化させることなく一定のコントラストを得ることができるようになる。
【００２４】
従来のようにエコー数を固定してパルスシーケンスを組み立てると、ショット数や位相エンコードマトリクスサイズの変化に伴って、大幅に無駄なエコーを収集することがあり、シーケンス長が延長したり、マルチスライス枚数及びＴ1 緩和時間が減少する等の不具合を引き起こす。また、位相エンコードマトリクスサイズに対して過少のエコーを収集することがあり、０データが詰められるエンコードの増加に伴って、Ｓ／Ｎ及び分解能の低下が引き起こされることがある。本実施例では、ショット数や位相エンコードマトリクスサイズの変化に伴って、エコー数を可変としてパルスシーケンスを組み立てるので、このような不具合は生じない。
【００２５】
次にショット数可変に対応した位相制御について説明する。
この位相制御により、ＤＣアーチファクト及びＦＩＤアーチファクトが低減される。ＤＣアーチファクトとは、受信器６内における直流アーチファクトをいう。また、ＦＩＤアーチファクトとは、 180°パルスを印加したとき、スライスの端の方では完全に磁化を 180°倒すことができないので、縦磁化から横磁化になってＦＩＤ信号が発生することがあり、この信号は画像化に寄与するエコー（本エコー）とは異なる位相シフトをするため画像上に本当のエコー成分と同様の形の虚像としてのアーチファクトとなって表れるが、このアーチファクトのことをいう。これらＤＣアーチファクト及びＦＩＤアーチファクトを低減するために、ＲＦパルスの位相及び受信器６内のアナログディジタル変換器（ＡＤ）の極性を制御することが一般的に行われている。
【００２６】
図３にＤＣアーチファクト低減の原理説明図であり、同図（ａ）には１回目の収集におけるＡＤ後の或るエンコードのエコー成分及びＤＣ成分と、２回目の収集におけるＡＤ後の同じエンコードのエコー成分及びＤＣ成分を示している。ＤＣ成分はエコー成分の位相に関わらず、一定の極性でＡＤに入力されるので、１回目と２回目とで、ＡＤの極性を反転させ、そして同図（ｂ）に示すように、両者を加算することにより、ＤＣ成分を打ち消すことができる。
【００２７】
図４（ａ）にCPMG(Carr-Purcell Meiboom-Gill) パルスシーケンス、図５（ａ）に位相反転CP(Carr-Purcell)パルスシーケンスのパターンＡ、図６（ａ）に位相反転CPパルスシーケンスのパターンＢそれぞれのＲＦパルスの位相、各エコーに対するＡＤの極性、ＡＤ前のエコー成分とＦＩＤ成分とＤＣ成分の極性、ＡＤ後のエコー成分とＦＩＤ成分とＤＣ成分の極性を示し、また各図（ｂ）に１NAQ 時の画像上における各成分の現象を示し、各図（ｃ）に偶数NAQ 時の画像上における各成分の現象を示している。
【００２８】
位相反転CPパルスシーケンスとは、全てのＲＦパルスを同じｘ´軸の正方向に沿ってかけるCPパルスシーケンスに対して、 180°パルスの位相を交互に 180°かえるパルスシーケンスをいい、これにより 180°パルスの精度誤差の累積的効果を抑える効果がある。ここではそのパターンＡ、Ｂ共に、１回目の収集時には最初の90°パルスをｘ´軸の正方向に沿ってかけ、２回目の収集時には最初の90°パルスをｘ´軸の負方向に沿ってかける。さらに、そのパターンＡでは１回目と２回目の収集間で各エコーに対して 180°パルスを同じ位相でかけ、一方、パターンＢでは１回目と２回目の収集間で各エコーに対して 180°パルスの位相を 180°かえて印加する。
【００２９】
また、CPMGパルスシーケンスとは、基本的にはCPパルスシーケンスと同じであるが、 180°パルスはｙ´軸の正方向、つまり最初の90°パルスと位相が90°ずれた方向にかけられるパルスシーケンスをいい、これによっても 180°パルスの精度誤差の累積的効果を抑える効果がある。ここでは、１回目の収集時には最初の90°パルスをｘ´軸の正方向に沿ってかけ、２回目の収集時には最初の90°パルスをｘ´軸の負方向に沿ってかける。
【００３０】
次に、ＡＤ極性を各パルスシーケンスで次のように制御する。CPMGパルスシーケンスでは１回目の収集では全てにエコーに対してＡＤ極性を正極性に設定し、２回目の収集では全てにエコーに対してＡＤ極性を負極性に設定する。位相反転CPパルスシーケンスではパターンＡ，Ｂとで同じように、ＡＤ極性を正／負交互に切り替え、さらに１回目と２回目とで同じ番目のエコーについて逆極性になるようにＡＤ極性を設定する。
【００３１】
シーケンサ１０はCPMGパルスシーケンス、位相反転CPパルスシーケンスのパターンＡ、位相反転CPパルスシーケンスのパターンＢをNAQ 数（アベレージング数）とショット数に応じて適宜選択的に使用する。
【００３２】
CPMGパルスシーケンスでは、奇数NAQ の場合、ＦＩＤ成分及びＤＣ成分共に、とぶが、消えない。偶数NAQ の場合、ＦＩＤ成分及びＤＣ成分共に相殺され、消える。
【００３３】
位相反転CPパルスシーケンスのパターンＡでは、奇数NAQ の場合、ＦＩＤ成分及びＤＣ成分共に、とぶが、消えない。偶数NAQ の場合、K-space の１ライン毎にＦＩＤ成分及びＤＣ成分の位相が反転するため最高周波数にエンコードされ、ＦＩＤ成分及びＤＣ成分共に相殺され消える。
【００３４】
位相反転CPパルスシーケンスのパターンＢでは、奇数NAQ の場合、ＦＩＤ成分はエコー成分に重なり、消えない。ＤＣ成分は、奇数ショットのとき、K-space の最高周波数にエンコードされ、とぶが消えない。偶数NAQ の場合、ＤＣ成分は変則的なエンコード順序であればとぶ。
【００３５】
これら各シーケンスのＦＩＤ成分、ＤＣ成分の現象を考察するに、偶数NAQ の場合、ショット数に関わらず、図４（ａ）のCPMGパルスシーケンス又は図５（ａ）の位相反転CPパルスシーケンスのパターンＡが選択される。また、奇数NAQ で、ショット数が偶数の場合、図４（ａ）のCPMGパルスシーケンスが選択される。奇数NAQ で、ショット数が奇数の場合、図６（ａ）の位相反転CPパルスシーケンスのパターンＢが選択される。
【００３６】
このようなパルスシーケンス及びパターンの選択により、ＤＣアーチファクト及びＦＩＤアーチファクトを好適に低減することができる。
次に撮影時間の短縮に関してfastSEを一例に説明する。一連のパルスシーケンスは一定の繰り返し時間ＴＲでショット数だけ繰り返されるが、従来説明にもあたように、ｎ≧２（通常はｎ＝２に設定される）として、ｎ−１回目までのパルスシーケンスはいわゆる空打ちとして、エコーは収集しないで、ｎ回目以降のパルスシーケンスで画像を生成するためのエコーを収集する。
【００３７】
本実施例では、この空打ちにおいて、位相補正及びＴ2 補正の情報を収集する、換言すると、空打ちと、位相補正及びＴ2 補正の情報を収集するためのパルスシーケンスとを共有することにより、空打ちと、位相補正及びＴ2 補正の情報を収集するためのパルスシーケンスとを別々に行う従来よりも撮影時間の短縮化を図る。
【００３８】
図７（ａ）に本実施例による空打ちのパルスシーケンスを示し、同図（ｂ）に実際に画像生成に使用されるエコーを収集するパルスシーケンスを示す。本実施例による空打ちのパルスシーケンスは、位相エンコード方向の傾斜磁場を全て０にして位相エンコードをかけずにエコーデータ（このデータは一般にtemplateデータと呼ばれる）、つまり位相補正及びＴ2 補正の情報を収集するためのパルスシーケンスに改良されている。この空打ちのパルスシーケンスは、位相エンコード方向の傾斜磁場を全て０にする点のみ相違し、他のＲＦパルス、スライス選択傾斜磁場、リード傾斜磁場及びエコー収集については同図（ｂ）の実際のパルスシーケンスと同じである。
【００３９】
これにより、空打ちとは別に、位相補正及びＴ2 補正の情報を収集するためのパルスシーケンスを実行する必要がなく、撮影時間の短縮化が図れる。
ここで、フーリエ変換前のエコーデータをｆ（kx,ky ）、フーリエ変換前のtemplateデータをｆt(kx) 、フーリエ変換後のエコーデータをＦ（ξ，η）、フーリエ変換後のtemplateデータをＦt(ξ，η）とする。
【００４０】
Ｔ2 補正後のエコーデータｆ´（kx,ky ）は、コンピュータシステム１１において１次元高速フーリエ変換処理（1st FFT ）に供される前のエコーデータｆ（kx,ky ）とtemplateデータｆt(kx) とにより、（５）式にしたがって計算される。なお、ｆt(kx) 0 は、０エンコードを占めるエコー番号のtemplateデータを示している。
【００４１】
ｆ´（kx,ky ）＝ｆ（kx,ky ）×｛Σｆt(kx) 0 ／Σｆt(kx) ｝…（５）
また、位相補正は、コンピュータシステム１１において１次元高速フーリエ変換処理（1st FFT ）に供された後のエコーデータとtemplateデータにより行われる。まず、フーリエ変換後のエコーデータとtemplateデータはそれぞれ（６）式、（７）式で与えられる。
【００４２】
Ｆ（ξ，η）＝Ｒe(Ｆ) ＋ｉ・Ｉm(Ｆ)                       …（６）
Ｆt(ξ，η）＝Ｒe(Ｆt)＋ｉ・Ｉm(Ｆt)                      …（７）
そして、位相補正は次の（８）式により行われる。
0003940442.tif 000002
このように空打ちにおいて、位相補正及びＴ2 補正の情報を収集する、換言すると、空打ちと、位相補正及びＴ2 補正の情報を収集するためのパルスシーケンスとを共有することにより、空打ちと、位相補正及びＴ2 補正の情報を収集するためのパルスシーケンスとを別々に行う従来よりも撮影時間の短縮化を図ることができる。
本発明は、上述した実施例に限定されることなく種々変形して実施可能であるのは勿論である。
【００４３】
【発明の効果】
請求項１に係る発明によれば、オペレータがショット数を指定すると、一連のシーケンスがこの指定されたショット数の回数だけ繰り返し実行される。したがって、オペレータは、ショット数Ｎと一連のシーケンスの繰り返し時間ＴＲとから、Ｎ×ＴＲにより撮影時間を容易に認識することができる。
【図面の簡単な説明】
【図１】本発明の一実施例による磁気共鳴イメージング装置の構成図。
【図２】シーケンサを主体としたパルスシーケンスの組み立てまでのフローチャートを示す図。
【図３】ＤＣアーチファクト低減の原理説明図。
【図４】  CPMGパルスシーケンスの場合のＲＦパルスの位相と各エコーに対するＡＤの極性との制御を示す図。
【図５】位相反転CPパルスシーケンスのパターンＡの場合のＲＦパルスの位相と各エコーに対するＡＤの極性との制御を示す図。
【図６】位相反転CPパルスシーケンスのパターンＢの場合のＲＦパルスの位相と各エコーに対するＡＤの極性との制御を示す図。
【図７】本実施例による空打ちのパルスシーケンスと実際に画像生成に用いるエコーを収集するパルスシーケンスを示す図。
【図８】従来の空打ちのパルスシーケンスと実際に画像生成に用いるエコーを収集するパルスシーケンスを示す図。
【符号の説明】
１…静磁場磁石、                ２…Ｘ・Ｙ・Ｚ軸傾斜磁場コイル、
３…ＲＦコイル、                ４…静磁場制御装置、
５…送信器、                    ６…受信器、
７…Ｘ軸傾斜磁場電源、          ８…Ｙ軸傾斜磁場電源、
９…Ｚ軸傾斜磁場電源、        １０…シーケンサ
１１…コンピュータシステム、    １２…表示部、
１３…入力部、                  ２０…コイルアセンブリ。
【図１】
0003940442.tif 000003
【図２】
0003940442.tif 000004
【図３】
0003940442.tif 000005
【図４】
0003940442.tif 000006
【図５】
0003940442.tif 000007
【図６】
0003940442.tif 000008
【図７】
0003940442.tif 000009
【図８】
0003940442.tif 000010