(57)【要約】【課題】  静磁場の不均一性の乱れに伴う共鳴周波数の変位により発生する不具合を解消するための磁気共鳴イメージング装置を提供すること。【解決手段】  水よりも小さな時定数と水と異なる共鳴周波数とを有するサンプルに、静磁場及び高周波を印加することにより、当該サンプルに発生したＦＩＤ信号を検出する周波数変位検出部１６と、検出されたＦＩＤ信号に基づいて、サンプルの共鳴周波数の基準時刻からの時間的推移を測定する周波数推移演算部２０５と、サンプルの共鳴周波数の基準時刻からの時間的推移に基づいて、所定の空間内の磁場不均一性による影響を補正するように、再構成部２０４、傾斜磁場コイル駆動装置１７、送信部１８、シムコイル制御部等を制御する制御部２０２と、を具備する。【選択図】      図１
2005270327.tif 000002
【特許請求の範囲】
【請求項１】
所定の空間内に配置された被検体に静磁場を印加すると共に、高周波及び傾斜磁場を所定のパルスシーケンスに従って印加し、前記被検体において発生する磁気共鳴信号を検出する撮影動作を繰り返し実行する磁気共鳴イメージング装置において、
前記被検体組織よりも小さな時定数（Ｔ１，Ｔ２）と前記被検体組織と異なる共鳴周波数とを有する物質に、前記静磁場及び高周波を印加することにより、前記物質の共鳴周波数の基準時刻からの時間的推移を測定する周波数推移測定手段と、
前記物質の共鳴周波数の基準時刻からの時間的推移に基づいて、前記所定の空間内の磁場不均一性による影響を補正する補正手段と、
を具備することを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記補正手段は、前記被検体から検出された前記磁気共鳴信号に対し、前記時間的推移に従って周波数を推移させた画像再構成を実行することを特徴とする請求項１記載の磁気共鳴イメージング装置。
【請求項３】
前記所定の空間内の磁場不均一性を制御するための電流シムコイルをさらに具備し、
前記補正手段は、
前記時間的推移に基づいて球面調和関数の展開係数を計算し、
前記展開係数に基づいて、前記高周波を印加するための電流、前記傾斜磁場を印加するための電流、前記電流シムコイルに印加される電流のうち、少なくとも一つを制御すること、
を特徴とする請求項１記載の磁気共鳴イメージング装置。
【請求項４】
前記所定の空間内の磁場不均一性を制御するための電流シムコイルをさらに具備し、
前記補正手段は、
前記被検体から検出された前記磁気共鳴信号に対し、前記時間的推移に従って周波数を推移させた画像再構成を実行し、
前記時間的推移に基づいて球面調和関数の展開係数を計算し、
前記展開係数に基づいて、前記高周波を印加するための電流、前記傾斜磁場を印加するための電流、前記電流シムコイルに印加される電流のうち、少なくとも一つを制御すること、
を特徴とする請求項１記載の磁気共鳴イメージング装置。
【請求項５】
前記周波数推移測定手段は、前記物質と、当該物質の近傍に配置された高周波コイルとからなり、前記被検体からの磁気共鳴信号を受信する受信コイル近傍の異なる少なくとも八箇所に配置されるセンサを有することを特徴とする請求項１乃至４のうちいずれか一項記載の磁気共鳴イメージング装置。
【請求項６】
前記物質はＭｅ４Ｓｉ又はＡＢＳ樹脂であることを特徴とする請求項１乃至５のうちいずれか一項記載の磁気共鳴イメージング装置。
【請求項７】
前記基準時刻は、前記電流シムコイルを利用したシミング時であることを特徴とする請求項１乃至６のうちいずれか一項記載の磁気共鳴イメージング装置。
【請求項８】
前記周波数推移測定手段は、前記撮影動作毎に前記時間的推移の測定を実行し、
前記補正手段は、前記撮影動作毎に前記補正を実行すること、
を特徴とする請求項１乃至７のうちいずれか一項記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、静磁場不均一性の乱れに伴う共鳴周波数の変化により発生する不具合を解消するための磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
磁気共鳴イメージング装置は、固有の磁気モーメントを持つ核の集団が一様な静磁場中に置かれたときに、特定の周波数で回転する高周波磁場のエネルギーを共鳴的に吸収する現象を利用して、物質の化学的及び物理的な微視的情報を映像化し、あるいは化学シフトスペクトラムを観測する装置である。
【０００３】
この磁気共鳴イメージング装置を使用した画像診断の手法として、例えば、ｆＭＲＩ（functional MRI）と呼ばれるものがある。このｆＭＲＩは、神経活動に伴う局所血行動態の変化に基づいて、脳の機能を画像化する技術であり、傾斜磁場の高入力、長時間連続スキャンを実行する。この様なスキャンにおいては、傾斜磁場コイルの中にある鉄シムの温度が上昇し、これによって透磁率が変化し、静磁場不均一性の乱れを発生させる。この静磁場不均一性の乱れは、例えばＥＰＩ（Echo Planer Imaging）画像の歪み、位置シフトとしてイメージングに影響を及ぼす。
【０００４】
この不具合を解消するための技術として、収集したエコーデータから磁場分布の０次成分（画像全体の位置シフトに対応）を推定し、位置シフトを補正する方法が提案されている。また、ハードウェアにて共鳴周波数のずれを補正する方法も別途提案されている。
【０００５】
しかしながら、上記技術は磁場分布の０次成分を補正するのみであり、一次以上の高次成分を補正するものではない。また、ｆＭＲＩ撮影後のスキャンに与える影響（１Ｈｚ／分程度）に対応する共鳴周波数の変化には対応していない。特に後者は、特に脂肪抑制（Chess・SPIR）不良となるおそれがある。
【発明の開示】
【発明が解決しようとする課題】
【０００６】
本発明は、上記事情を鑑みてなされたもので、磁場分布の一次以上の高次成分による影響を含む、静磁場の不均一性の乱れに伴う共鳴周波数の変位により発生する不具合を解消することができる磁気共鳴イメージング装置を提供することを目的としている。
【課題を解決するための手段】
【０００７】
本発明は、上記目的を達成するため、次のような手段を講じている。
【０００８】
本発明の視点は、所定の空間内に配置された被検体に静磁場を印加すると共に、高周波及び傾斜磁場を所定のパルスシーケンスに従って印加し、前記被検体において発生する磁気共鳴信号を検出する撮影動作を繰り返し実行する磁気共鳴イメージング装置において、前記被検体組織よりも小さな時定数（Ｔ１，Ｔ２）と前記被検体組織と異なる共鳴周波数とを有する物質に、前記静磁場及び高周波を印加することにより、前記物質の共鳴周波数の基準時刻からの時間的推移を測定する周波数推移測定手段と、前記物質の共鳴周波数の基準時刻からの時間的推移に基づいて、前記所定の空間内の磁場不均一性による影響を補正する補正手段と、を具備することを特徴とする磁気共鳴イメージング装置である。
【発明の効果】
【０００９】
以上本発明によれば、静磁場の不均一性の乱れに伴う共鳴周波数の変位により発生する不具合を解消するための磁気共鳴イメージング装置を実現できる。
【発明を実施するための最良の形態】
【００１０】
以下、本発明の実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００１１】
図１は、本実施形態に係る磁気共鳴イメージング装置１０の構成を示した図である。図１に示すように、本磁気共鳴イメージング装置１０は、静磁場磁石１１、冷却系制御部１２、傾斜磁場コイル１３、全身用高周波（RF）コイル１４、高周波受信コイル１５、周波数変位検出部１６、傾斜磁場コイル駆動装置１７、送信部１８、受信部１９、演算装置２０、表示部２４を具備している。
【００１２】
静磁場磁石１１は、静磁場を発生する磁石であり、一様な静磁場を発生する。この静磁場磁石１１には、例えば永久磁石、超伝導磁石等が使用され、図示していない冷却系によって冷却される。
【００１３】
冷却系制御部１２は、静磁場磁石１１を冷却するための冷却系（図示せず）を制御する。
【００１４】
傾斜磁場コイル１３は、静磁場磁石１１よりも短軸な磁場コイルであり、静磁場磁石１１の内側に設けられる。傾斜磁場コイル１３は、傾斜磁場コイル駆動装置１７から供給されるパルス電流に基づいて、互いに直交するＸ，Ｙ，Ｚの三方向に線形傾斜磁場分布を持つ傾斜磁場を形成する。この傾斜磁場コイル１３が発生する傾斜磁場によって、信号発生部位（位置）が特定される。
【００１５】
なお、Ｚ軸方向は、本実施形態では静磁場の方向と同方向（被検体の体軸方向）にとるものとする。また、本実施形態において、傾斜磁場コイル１３及び静磁場磁石１１は円筒形をしているものとする。また、傾斜磁場コイル１３は、所定の支持機構によって真空中に配置されていてもよい。これは、静音化の観点から、パルス電流の印加によって発生する傾斜磁場コイル１３の振動を、音波として外部に伝播させないためである。
【００１６】
全身用RFコイル１４は、被検体の撮像領域に対して、磁気共鳴信号を発生させるための高周波パルスを印加するコイルである。また、例えば腹部等を撮影する場合には、受信コイルとしても使用される。
【００１７】
高周波受信コイル（ＲＦ受信コイル）１５は、部位別に専用の形状を有した移動可能なサーフェスコイルである。高周波受信コイル１５は、所定の器具や撮影者の手によって撮影領域近傍に配置され、被検体に発生した磁気共鳴信号を受信する。
【００１８】
周波数変位検出部１６は、高周波受信コイル１５に複数設けられ、基準時刻からの周波数変位を検出する。周波数変位検出部１６は、周波数変位を検出するためのサンプルとしてのコア部１６０と、受信コイルとして機能する微小コイル部１６１とを有している。
【００１９】
図２は、周波数変位検出部１６の構成を、図３は、周波数変位検出部１６の高周波受信コイル１５への設置例を、それぞれ図示している。
【００２０】
図２に示すように、周波数変位検出部１６は、例えば直径５ｍｍ程度の微小な円柱形状のコア部１６０、コア部１６０を取り巻き、当該コア部１６０に発生するＮＭＲ信号を受信するためのアンテナである微小コイル部１６１を有している。コア部１６０は、Ｔ１、Ｔ２が短く、磁気共鳴によるプロトン信号を発生し、化学シフト量が水に比して大きな素材から構成される。好適な例としては、Ｍｅ４Ｓｉ（化学シフト量：４．３７ｐｐｍ）、ＡＢＳ樹脂（Ｔ２＝１．６ｍｓ）が挙げられる。また、コア部１６０の大きさは、当該コア部１６０が発生するプロトン信号を高周波受信コイル１５が検出することができない程度に、十分に小さくする必要がある。
【００２１】
また、図３に示すように、周波数変位検出部１６は、例えば高周波受信コイル１５の開口周縁部に設置される。磁場の分布状況（静磁場の不均一性の乱れ）の検出が目的であることから、周波数変位検出部１６は、対称に設置されることが好ましい。また、各座標軸方向に関する磁場不均一性の一次以上の成分を補正するため、周波数変位検出部１６は、少なくとも８個以上設置されることが好ましい。
【００２２】
傾斜磁場コイル装置電源１７は、傾斜磁場を形成するためのパルス電流を発生し、傾斜磁場コイル１３に供給する。また、傾斜磁場コイル装置電源１７は、制御部２０２の制御に従って、傾斜磁場コイル１３に供給するパルス電流の向きを切替えることにより、傾斜磁場の極性を制御する。
【００２３】
送信部１８は、発振部、位相選択部、周波数変換部、振幅変調部、高周波電力増幅部（それぞれ図示せず）を有しており、ラーモア周波数に対応する高周波パルスを全身用RFコイル１４に送信する。当該送信によって全身用RFコイル１４から発生した高周波によって、被検体の所定原子核の磁化は、励起状態となる。
【００２４】
受信部１９は、増幅部、中間周波数変換部、位相検波部、フィルタ、Ａ／Ｄ変換器（それぞれ図示せず）を有し、高周波受信コイル１５又は周波数変位検出部１６から受信した各磁気共鳴信号（高周波信号）に対して、個別に所定の信号処理を施す。すなわち、受信部１９は、核の磁化が励起状態から基底状態に緩和するとき放出する磁気共鳴信号に対して、増幅処理、発信周波数を利用した中間周波数変換処理、位相検波処理、フィルタ処理、Ａ／Ｄ変換処理を施す。
【００２５】
演算装置２０は、記憶部２０１、制御部２０２、データ収集部２０３、再構成部２０４、周波数推移演算部２０５、入力部２０７を有している。
【００２６】
記憶部２０１は、受信部１９を介して得られた再構成前の磁気共鳴信号データ再構成後の磁気共鳴画像データ等を患者毎に記憶する。
【００２７】
制御部２０２は、図示していないＣＰＵ、メモリ等を有しており、システム全体の制御中枢として、本磁気共鳴イメージング装置を静的又は動的に制御する。特に、制御部２０２は、パラレルイメージングを行う場合には、感度分布の異なる複数のRF受信コイルにより、並列的に磁気共鳴信号を受信・処理するための制御を実行する。
【００２８】
また、制御部２０２は、周波数推移演算部２０５によって推定された撮影断面における共鳴周波数のシフト量に基づいて、画像再構成部２０４、及び傾斜磁場コイル装置電源１７、送信部１８、図示していないシムコイル制御部をそれぞれ制御する。
【００２９】
データ収集部２０３は、受信部１９によってサンプリングされたディジタル信号を収集する。
【００３０】
再構成部２０４は、後処理すなわちフーリエ変換等の再構成を実行し、被検体内の所望核スピンのスペクトラムデータあるいは画像データを求める。また、パラレルイメージングを行った場合には、再構成部２０４は、各コイルからの磁気共鳴信号からそれぞれ画像を再構成し、その後、各コイルの感度分布を使って、得られた複数枚の画像の後処理として展開処理を行い１枚の画像を生成する。ここで、パラレルイメージングとは、感度分布の異なる複数のRF受信コイルを用いて、位相エンコードを間引いたシーケンスを実行し、行列演算によって折り返しアーチファクトを除去する展開処理を行うことにより撮影時間を短縮化する技術である。
【００３１】
なお、再構成部２０４が実行する後処理は、必要に応じて、周波数推移演算部２０５によって得られた周波数変位を考慮して実行される。すなわち、再構成部２０４は、周波数推移演算部２０５によって得られた周波数変位を反映したフーリエ変換等を実行する。
【００３２】
周波数推移演算部２０５は、基準時（例えば、シミング時）において、各周波数変位検出部１６によって検出される中心周波数（以下、ＣＦ）を記録し、診断画像撮影時において各周波数変位検出部１６によって検出されるＣＦとの差分値を演算する。周波数推移演算部２０５は、この差分値に基づいて、コア部１６０の基準時刻からの共鳴周波数の推移を計算する。
【００３３】
また、周波数推移演算部２０５は、得られたコア部１６０の基準時刻からの共鳴周波数の推移に基づいて、上述する手法により撮影断面における共鳴周波数のシフト量を推定する。
【００３４】
入力部２０７は、オペレータからの各種指示・命令・情報をとりこむため入力装置（マウスやトラックボール、モード切替スイッチ、キーボード等）を有している。
【００３５】
表示部２４は、制御部２０２を介して演算装置２２から入力したスペクトラムデータあるいは画像データ等を表示する出力装置である。
【００３６】
（静磁場不均一性の乱れによる影響を補正する補正処理）
次に、本磁気共鳴イメージング装置が実行する、静磁場不均一性の乱れによる影響を補正する補正処理について説明する。この補正処理は、周波数変位検出部１６を用いた共鳴周波数シフト量の測定、撮影断面における共鳴周波数シフト量の推定、共鳴周波数シフト量の補正の各処理によって実現されるものである。以下、各処理について説明する。
【００３７】
（サンプルの共鳴周波数シフト量の測定）
周波数変位検出部１６を用いたサンプル（コア部１６０）の共鳴周波数シフト量の測定は、基準時刻からの受信信号（ＦＩＤ信号）の位相変位に基づいて実行される。すなわち、例えば本実施形態において基準時刻をシミング時とすれば、当該シミング時において、被検体の組織を励起せず、コア部１６０のみを励起する空間非選択パルスが印加されると、コア部１６０は短いＴ１、Ｔ２のＦＩＤ信号（プロトン信号）を発生し、微小コイル部１６１はこれを検出する。この微小コイル部１６１による信号の検出タイミングは、図示していないＰｉｎダイオードによって制御される。
【００３８】
なお、上記の様に基準時刻をシミング時としたのは、シミング時が、最も静磁場均一性が高い状態であるからである。
【００３９】
次に、本撮影前に空間非選択パルス、より具体的には図４に示すようなタイミングチャートに従うＲＦプリパルスを印加し、これによって発生する各コア部１６０からのＦＩＤ信号を、各微小コイル部１６１によって測定する。その後、各コア部１６０からのＦＩＤ信号以外の信号を残存させないためのスポイラー(spoiler)パルスが印加される。
【００４０】
なお、プリパルス印加に伴う各コア部１６０からのＦＩＤ信号の測定は、精度向上の観点から、所定の時間間隔で複数回、好ましくは本撮影毎の直前において実行されることが好ましい。また、各コア部１６０からのＦＩＤ信号のデータ収集時間（ＴＥ：９０度パルス印加時刻からデータ収集完了までの時間）は、１ｍｓ〜２ｍｓであるとする。
【００４１】
データ収集後、得られた検波後の各ＦＩＤ信号の１〜２ｍｓ間（すなわち、１ｍｓあたり）の複素位相差Δφ［degree］を算出し、さらに、基準時刻から一定時間経過後のコア部１６０の共鳴周波数シフト量Δfｃｏｒｅ[Hz]を、基準時における複素位相差をΔφ0、一定時間経過後の複素位相差をΔφ1として、次の式（１）により算出する。
【数１】
2005270327.tif 000003
【００４２】
但し、ｔ＝０．００１ｓ、Δφ０、Δφ１の単位はdegreeとする。従って、例えばｔ＝０．００１ｓ、Δφ１−Δφ０＝２．７８［Hz/degree］である場合には、エリアシングなしに測定できる範囲は−５００Hzから＋５００Hzとなる。
【００４３】
この様な構成により、基準時刻からの、コア部１６０の共鳴周波数の経時的なシフト量（推移）Δfｃｏｒｅを把握することができる。特に、各本撮影の直前に実施される空間非選択パルス印加タイミングに合わせてコア部１６０からのＦＩＤ信号を測定することで、直後に実施される本撮影における共鳴周波数の経時的なシフト量を容易かつ迅速に把握することができる。
【００４４】
また、コア部１６０には、磁気共鳴によるプロトン信号を発生し、図５に示すように化学シフト量が水と異なり（共鳴周波数が水と異なる）、Ｔ１、Ｔ２値が水より短い素材（Ｍｅ４Ｓｉ、ＡＢＳ樹脂等）が使用される。従って、プリパルスを印加した場合であっても、撮像対象の被検体組織に影響を与えずにコア部１６０の周波数推移を測定することが可能である。
【００４５】
また、周波数変位検出部１６を用いたＦＩＤ信号の測定では、コア部１６０が直近にあることから、各コア部１６０を取り巻く微小コイル１６１は当該コア部１６０が発生するプロトン信号を十分に検出することができる。一方、コア部１６０が発生するプロトン信号は、コア部１６０のサイズが十分微小であることから、高周波受信コイル１５が検出するには非常に小さい。従って、診断画像イメージングに対して問題とはならない。
【００４６】
さらに、診断画像イメージング用のパルスでは、コア部１６０のＴ２に比してＴＥが長い。このため、周波数変位検出部１６は、診断画像イメージングにより発生する信号を検出することはなく、コア部１６０からのプロトン信号検出については問題とならない。
【００４７】
なお、本測定処理においては、ハードウェア簡素化のため、全てのタイミングにおいて全センサのデータを同時に収集する必要は無く、１つのセンサの収集系に複数のセンサを接続して、センサを切り替えて使用しても良い。
【００４８】
（撮影断面における共鳴周波数シフト量の推定）
撮影断面における共鳴周波数シフト量の推定については、コア部１６０の共鳴周波数シフト量の測定重み付け加算による手法と、静磁場分布を示す球面楕円関数の展開係数算出による手法とがある。
【００４９】
コア部１６０の共鳴周波数シフト量の測定重み付け加算による推定手法は、n箇所（本実施形態では、ｎ＝８）に設置したセンサの共鳴周波数シフト値を重み付け加算するものである。各センサの検出値の重みに対応する重み付け係数は、センサの位置と撮影断面までの距離に応じて定義される。
【００５０】
すなわち、n個のセンサにおける周波数シフトをΔfi[Hz](i=1,2,…,n)とすれば、撮影断面における周波数シフトΔfｓｌｉｃｅ[Hz]は、次の式（２）によって表現することができる。
【数２】
2005270327.tif 000004
【００５１】
こうして推定された撮影断面における周波数シフトΔfｓｌｉｃｅは、後述する共鳴周波数シフト量の補正処理において利用される。
【００５２】
一方、静磁場分布を示す球面楕円関数の展開係数算出による手法は、n個のセンサにおける周波数シフトをΔfi[Hz](i=1,2,…,n)を利用して、直交座表系における球面調和関数によりｎ次（n=1,2,…）の磁場分布成分に展開するものである。ここでは、０次（定数）、１次（Ｘ，Ｙ，Ｚ）、２次（Ｚ２，Ｘ２，Ｙ２，ＺＸ，ＺＹ，ＸＹ）の各磁場分布成分の展開係数を求めるものとするが、実際に展開係数を求める成分の次数は、周波数変位検出部１６の個数、設置場所、シムコイルの有無等によって決定される。なお、０次成分は励起周波数（中心周波数）のオフセット、１次成分は傾斜磁場出力のオフセット、２次（以降）の高次成分はシムコイルの電流変化にそれぞれ対応する。
【００５３】
各磁場分布成分は、次の式（３）によって示される行列計算の係数マトリックスＡｉｊの各成分を算出することで、求めることができる。
【００５４】
Ｂｉ＝Ａｉｊｘｊ                                （３）
ここで、Ｂｉは測定点における周波数シフト量（i=1,2,・・・、測定点数）であり、Ａｉｊは係数マトリックス（j=1,2,・・・、磁場の分布成分数）であり、ｘｊは磁場分布成分である。
【００５５】
こうして得られた磁場分布成分（今の場合、０次、１次、２次の各成分）の係数を用いて、後述する共鳴周波数シフト量の補正処理として、本撮影時における励起周波数のオフセット、傾斜磁場出力のオフセット、シムコイルの電流をそれぞれ制御する。
【００５６】
（共鳴周波数シフト量の補正処理）
周波数推移演算部２０５によって推定された、撮影断面における共鳴周波数シフト量を用いた補正処理について説明する。この補正処理は、本撮影によって得られた磁気共鳴信号に対して行う事後的なものと、本撮影前に行うものとの二つに分類することができる。
【００５７】
本撮影後に行う補正処理とは、既述の式（２）に従って計算された撮影断面における周波数シフトΔfを反映させた後処理、すなわちフーリエ変換による画像再構成を行うものである。また、本撮影前に行う補正処理とは、式（３）に従って計算された０次磁場分布成分を用いて励起周波数（中心周波数）のオフセットを、１次磁場分布成分を用いて傾斜磁場出力のオフセットを、１次磁場分布成分を用いてシムコイルへの印加電流をそれぞれ制御し、静磁場の歪みを補正するものである。
【００５８】
この二つの補正処理は、原則として両立するものではない。しかし、例えば既述の式（２）に従って計算された撮影断面における周波数シフトΔfｓｌｉｃｅを用いて０次磁場分布成分の補正を行い、残りの１次、２次の磁場分布成分については式（３）に従って計算し、これを用いて傾斜磁場出力のオフセット及びシムコイルへの印加電流を制御する等、矛盾しない範囲で一部を組み合わせて実行する構成であってもよい。
【００５９】
（撮影動作）
次に、上記静磁場不均一性の乱れを補正する補正処理を含む、一連の撮影動作について説明する。
【００６０】
図６は、本磁気共鳴イメージング装置による一連の撮影動作において実行される各処理の流れを示したフローチャートである。同図に示すように、まず、天板へ被検体を搭載し、静磁場磁石１１内の適切な位置に当該被検体を配置すると共に、各種撮影条件を入力部２０７より入力する（ステップＳ１）。
【００６１】
次に、静磁場分布を調整するためのシミングを実行し（ステップＳ２）、当該シミングと連動して図４に示す空間非選択パルスを印加して、基準時刻（シミング時）におけるコア部１６０からのＦＩＤ信号（基準ＦＩＤ信号）を収集する（ステップＳ３）。
【００６２】
次に、例えば共鳴周波数シフト量測定を目的とした空間非選択パルス印加によるプレスキャンが実行され（ステップＳ４）、その直後に検出されたコア部１６０から検出されたＦＩＤ信号及び基準ＦＩＤ信号とに基づいて、共鳴周波数シフト量の計算が実行される（ステップＳ５）。
【００６３】
次に、計算されたコア部１６０の共鳴周波数シフト量に基づいて、撮影断面に関する補正量、すなわち撮影断面における周波数シフト量及び２次までの各展開係数を計算する（ステップＳ６）。また、得られた各展開係数に基づいて、基準時からずれた磁場分布の不均一性を当該基準時の状態に戻すように、励起周波数のオフセット、傾斜磁場出力のオフセット、シムコイルの電流をそれぞれ制御し、本撮影を実行する（ステップＳ７）。このとき、撮影断面における周波数シフト量を０次成分の展開係数とする。
【００６４】
次に、本撮影によって得られた磁気共鳴信号に対し画像再構成等の後処理が実行され（ステップＳ８）、得られた画像が表示部２４に表示される（ステップＳ９）。
【００６５】
以下、本撮影を繰り返し実行する場合には、Ｓ４乃至Ｓ９までの各処理が繰り返し実行され、各本撮影において静磁場不均一性の乱れによる影響を補正する補正処理が実行されることになる。
【００６６】
以上述べた構成によれば、以下の効果を得ることができる。
【００６７】
本磁気共鳴イメージング装置によれば、サンプル（コア部）を用いて、本撮影直前に基準時からの撮影断面における共鳴周波数のシフト量を推定し、これに基づいて当該シフト量を考慮した後処理、及び撮影断面の励起周波数のオフセット、傾斜磁場出力のオフセット、シムコイルへの印加電流の制御のうち、少なくとも一方を実行する。従って、撮影中に静磁場不均一性により撮影断面の共鳴周波数が推移した場合であっても、これに影響をうけることなく、歪みの程度がシミング直後に撮影したものに相当するＭＲＩ画像を常に提供することができる。特に、ＥＰＩ以外のシーケンスにおいても、上記補正処理を実行することで、安定した脂肪抑制を実現することができる。
【００６８】
また、本磁気共鳴イメージング装置では、上記補正処理を各本撮影の直前において実行する。従って、直後に実施される本撮影に適合した、好適な補正処理を実現することができる。
【００６９】
さらに、本磁気共鳴イメージング装置では、周波数変位検出コイルを用いて、静磁場不均一性による撮影断面の共鳴周波数推移をハード的に把握している。従って、シーケンスに影響されることはなく、どのようなスキャンシーケンスであっても本補正処理を用いことができる。その結果、歪みのない好適なＭＲＩ画像を常に提供することができる。
【００７０】
また、図６では、プレスキャンにて本撮影前に１回だけ共鳴周波数シフト量を測定しているが、時間的に連続して行うダイナミックスキャンにおいて、毎時相又は、数時相毎に共鳴周波数シフト量を測定し、補正を行ってもよい。
【００７１】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【図面の簡単な説明】
【００７２】
【図１】図１は、本実施形態に係る磁気共鳴イメージング装置１０の構成を示した図である。
【図２】図２は、周波数変位検出部１６の構成を図示している。
【図３】図３は、周波数変位検出部１６の設置の様子を図示している。
【図４】図４は、コア部１６０の共鳴周波数シフト量の測定において印加されるプリパルスのシーケンス図である。
【図５】図５は、コア部１６０の共鳴周波数、及び水の共鳴周波数の分布を示した図である。
【図６】図６は、本磁気共鳴イメージング装置による一連の撮影動作において実行されえる各処理の流れを示したフローチャートである。
【符号の説明】
【００７３】
１０…磁気共鳴イメージング装置、１１…静磁場磁石、１３…傾斜磁場コイル、１４…全身用高周波（RF）コイル、１５…高周波受信コイル、１６…周波数変位検出部、１７…傾斜磁場コイル駆動装置、１８…送信部、１９…受信部、２０…演算装置、２４…表示部、１６０…コア部、１６１…微小コイル部
【図１】
2005270327.tif 000005
【図２】
2005270327.tif 000006
【図３】
2005270327.tif 000007
【図４】
2005270327.tif 000008
【図５】
2005270327.tif 000009
【図６】
2005270327.tif 000010