特許 6234214 磁気共鳴イメージング装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6234214

(24)【登録日】2017年11月2日

(45)【発行日】2017年11月22日

(54)【発明の名称】磁気共鳴イメージング装置

(51)【国際特許分類】

   A61B 5/055 20060101AFI20171113BHJP

【ＦＩ】

   A61B5/05 370

   A61B5/05 311

【請求項の数】4

【全頁数】13

(21)【出願番号】特願2013-265444(P2013-265444)

(22)【出願日】2013年12月24日

(65)【公開番号】特開2015-119845(P2015-119845A)

(43)【公開日】2015年7月2日

【審査請求日】2016年10月26日

(73)【特許権者】

【識別番号】000005108

【氏名又は名称】株式会社日立製作所

(74)【代理人】

【識別番号】110000350

【氏名又は名称】ポレール特許業務法人

(72)【発明者】

【氏名】平井 甲亮

【審査官】 安田 明央

(56)【参考文献】

【文献】 特開２００９−２２６１９９（ＪＰ，Ａ）

【文献】 特開２０１１−２４５２７３（ＪＰ，Ａ）

【文献】 特開昭６３−２５２１４４（ＪＰ，Ａ）

【文献】 特開２００３−１０２７０２（ＪＰ，Ａ）

【文献】 特開平０５−０５６９４９（ＪＰ，Ａ）

【文献】 米国特許第０５３３４９３６（ＵＳ，Ａ）

【文献】 特開２００４−３２９６１３（ＪＰ，Ａ）

【文献】 特開平０６−０２２９２６（ＪＰ，Ａ）

【文献】 米国特許第０５４９７７７３（ＵＳ，Ａ）

【文献】 米国特許第０４６１６１８２（ＵＳ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ ５／０５５

(57)【特許請求の範囲】

【請求項1】

静磁場発生手段と、傾斜磁場発生手段と、高周波磁場発生手段と、データ取得シーケンスに従って被検体からのエコー信号を検出する検出手段と、前記エコー信号検出手段により検出されたエコー信号に基づき、前記被検体の画像を再構成する信号処理手段とを備え、
前記データ取得シーケンスは、励起ＲＦパルスと再収束ＲＦパルスとを含み、
前記信号処理手段は、事前計測された前記静磁場発生手段と前記傾斜磁場生手段に起因した磁場歪みによって生じる画像歪み位置と撮影ＦＯＶを含む撮像条件とに基づいて、前記撮影ＦＯＶの範囲内に折り返しアーチファクトとして生じる位置を予測する判定基準値を算出し、
前記判定基準値に基づき、前記再収束ＲＦパルスのＲＦ照射位相を前記データ取得シーケンスが有するパルス条件に応じて変化させることにより、前記アーチファクトの位置を前記撮影ＦＯＶの範囲外に移動させることを特徴とする磁気共鳴イメージング装置。

【請求項2】

前記撮影ＦＯＶの範囲外への移動は、前記ＲＦ照射位相の調整により、前記ＦＯＶサイズの１/２移動させることを特徴とする請求項１記載の磁気共鳴イメージング装置。

【請求項3】

前記折り返しアーチファクトの位置が前記画像の関心領域上に出現する場合、前記被検体の画像撮影を偶数回行い、
前記エコー信号を取得してなるｋ空間におけるエコー配列の全てに対して、偶数加算目の撮影ではＲＦ照射位相が、奇数加算目のＲＦ照射位相に対して１８０度反転させることで前記折り返しアーチファクトを消失させることを特徴とする請求項１記載の磁気共鳴イメージング装置。

【請求項4】

前記信号処理手段で再構成された画像を表示する表示部を備え、
ユーザーが前記画像の関心領域上に前記折り返しアーチファクトを視認した場合、前記表示部上のユーザーパラメータからパラメータ選択によって、前記ＲＦ照射位相を変化させることにより、前記アーチファクトの位置を関心領域外に移動させることを特徴とする請求項１記載の磁気共鳴イメージング装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、磁気共鳴イメージング装置（以下、ＭＲＩ装置という）において、静磁場不均一および傾斜磁場非線形性によって生じるアーチファクトを抑制する手段を備えたＭＲＩ装置に関する。

【背景技術】

【0002】

水平磁場ＭＲＩにおける被検者の開放感を重視するために、ガントリのＺ軸長を短縮したショートガントリ型ＭＲＩが主流になっている。しかしながら、ショートガントリ型ＭＲＩは、静磁場均一空間、及び、傾斜磁場線形領域が狭まることで、磁場歪みによるアーチファクトが生じやすくなる問題がある。特に、Ｚ軸を位相エンコード方向に設定したＳＡＧ及びＣＯＲ断面のFast Spin Echo撮影において、撮影ＦＯＶ内の関心領域に輝点アーチファクトの折り返しが生じて診断の妨げとなる事がある。

【0003】

本アーチファクトを回避するために、ハードウェア、撮影シーケンス、再構成処理、及び、撮影方法の見地からいくつかの対策が検討されている。

【0004】

・ハードウェアの見地からは、フェーズドアレイコイル自身のデザイン最適化や、受信エレメント選択の最適化によって、磁場歪み領域からのＮＭＲ信号の受信を制限する方法がある。しかしながら、関心領域の信号低下にも繋がり、最適化の限界がある。

【0005】

・撮影シーケンスの見地からは、励起ＲＦパルスと再収束ＲＦパルスの２つの励起断面が、磁場歪み領域で重ならないようにする事で、Spin Echo信号の発生を抑制する方法がある（特許文献１、非特許文献１）。しかしながら、励起ＲＦパルスと再収束ＲＦパルスの励起帯域の違いによって、再収束ＲＦパルス単一の影響で生じるＦＩＤ(Free Induction Decay)信号に対しては効果が見込めない事に加え、Spin Echo信号を抑制する事で、逆にＦＩＤ信号によるアーチファクトが増強する事が考えられる。

【0006】

・再構成処理の見地からは、事前計測で取得したフェーズドアレイコイルの感度分布データを利用してアーチファクトを抑制する方法がある（特許文献２、３、非特許文献２、３）。しかしながら、アーチファクト強度は画像種や計測条件の依存性が高い場合には、アーチファクトを除去しきれない。また、パラレルイメージングの手法を使う場合、関節撮影等、高分解能撮影において、細かな折り返しを展開するには展開精度の限界があり、現実的な手法ではない。

【0007】

・撮影方法の見地からは、ＲＦブランケットで磁場歪み領域に覆う事でＲＦ照射を遮蔽する方法も考えられるが、ＳＡＲ(Specific Absorption Rate)等、安全面の問題が残る。

【先行技術文献】

【特許文献】

【0008】

【特許文献1】米国特許出願公開第２０１２／００２５８２６号明細書

【特許文献2】米国特許第７、２５０、７６２号明細書

【特許文献3】米国特許第６、１３４、４６５号明細書

【非特許文献】

【0009】

【非特許文献1】Novena Rangwala et al. “Reduction of Fast Spin Echo Cusp Artifact Using a Slice-Tilting Gradient”, Magn Reson Med 2010; 64: 220-228

【非特許文献2】David Atkinson et al.” Coil-Based Artifact Reduction” Magn Reson Med 2004; 52: 825-830

【非特許文献3】Larkman et al. “Elimination of Magnetic Field Foldover Artifacts in MR Images”, J Magn Reson Imaging 2000; 12: 795-797.

【発明の概要】

【発明が解決しようとする課題】

【0010】

Fast Spin Echoシーケンスを用いた場合、磁場歪み領域における再収束ＲＦパルスの影響で折り返しアーチファクト生じることがある。関心領域上に重畳されて診断の妨げとなるアーチファクトを拡張ＦＯＶに移動することができれば、見かけ上はアーチファクトを抑制する事ができる。

【0011】

上述した従来技術では、関心領域の信号低下や、ＦＩＤ信号によるアーチファクトの増大や、画像種や計測条件の依存性が高いことや、ＳＡＲ(Specific Absorption Rate)等の安全面の課題が解消されない。

【0012】

そこで、本発明の目的は、関心領域の信号低下がなく、ＦＩＤ信号によるアーチファクトの増大を抑制し、画像種や計測条件の依存性が少なく、ＳＡＲ等の問題が生じない撮像が可能な磁気共鳴イメージング装置を提供することである。

【課題を解決するための手段】

【0013】

上記課題を解決するために、本発明になる磁気共鳴イメージング装置の主な特徴は以下の通りである。

【0014】

静磁場発生手段と、傾斜磁場発生手段と、高周波磁場発生手段と、データ取得シーケンスに従って被検体からのエコー信号を検出する検出手段と、エコー信号検出手段により検出されたエコー信号に基づき、被検体の画像を再構成する信号処理手段とを備え、データ取得シーケンスは、励起ＲＦパルスと再収束ＲＦパルスとを含み、信号処理手段は、事前計測された静磁場発生手段と傾斜磁場生手段に起因した磁場歪みによって生じる画像歪み位置と撮影ＦＯＶを含む撮像条件とに基づいて、撮影ＦＯＶ内に折り返しアーチファクトとして生じる位置を予測する判定基準値を算出し、判定基準値に基づき、再収束ＲＦパルスの照射位相をデータ取得シーケンスが有するパルス条件に応じて変化させることにより、アーチファクトを前記撮影ＦＯＶの範囲外に移動させることを特徴とする。

【発明の効果】

【0015】

本発明によれば、関心領域の信号低下がなく、ＦＩＤ信号によるアーチファクトの増大を抑制し、画像種や計測条件の依存性が少なく、ＳＡＲ等の問題が生じない撮像が可能な磁気共鳴イメージング装置を提供できる。

【図面の簡単な説明】

【0016】

【図1】本発明に係るＭＲＩ装置の全体概要を示すブロック図。

【図2】Fast Spin Echoのシーケンス図。

【図3】（ａ）は撮影FOV、及び、アーチファクト発生位置を示す図、（ｂ）は撮影FOV、拡張FOV、アーチファクト発生位置、及び、折り返し位置を示す図、（ｃ）は撮影FOV、拡張FOV、アーチファクト発生位置、折り返し位置、及び、折り返し移動位置を示す図。

【図4】FOV/2だけFID信号の折りかえしを移動させるためのｋ空間上の位相情報を示す図。

【図5】Sequential計測時のエコー配列と位相反転（網掛）を示す図。

【図6】Centric計測時のエコー配列と位相反転（網掛）を示す図。

【図7】偶数回加算時のFID信号を消失させるための、ｋ空間上の位相情報を示す図。

【図8】アーチファクト回避のためのＲＦ照射位相計算のフローチャート。

【図9】(ａ)は奇数のShot数、（ｂ）は偶数のShot数の場合のSequential計測時の照射位相を示す図。

【図10】(ａ)はShot数/2が奇数、（ｂ）はShot数/2が偶数の場合のCentric計測時の照射位相を示す図。

【発明を実施するための形態】

【実施例1】

【0017】

＜ＭＲＩ装置のブロック図＞
以下、添付図面に従って本発明のＭＲＩ装置の好ましい実施例について詳説する。なお、発明の実施例を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。

【0018】

最初に、本発明に係るＭＲＩ装置を図１に基づいて説明する。図１は、本発明に係るＭＲＩ装置の一実施例の全体構成を示すブロック図である。

【0019】

このＭＲＩ装置は、ＮＭＲ現象を利用して被検体101の断層画像を得るもので、図1に示すように、静磁場発生磁石102と、傾斜磁場コイル103及び傾斜磁場電源109と、ＲＦ送信コイル104及びＲＦ送信部110と、ＲＦ受信コイル105及び信号処理部107と、計測制御部111と、全体制御部112と、表示・操作部118と、被検体101を搭載する天板を静磁場発生磁石102の内部に出し入れするベッド106と、を備えて構成される。

【0020】

静磁場発生磁石102は、垂直磁場方式であれば被検体101の体軸と直交する方向に、水平磁場方式であれば体軸方向に、それぞれ均一な静磁場を発生させるもので、被検体101の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。

【0021】

傾斜磁場コイル103は、ＭＲＩ装置の実空間座標系(静止座標系)であるＸ、Ｙ、Ｚの３軸方向に巻かれたコイルであり、それぞれの傾斜磁場コイルは、それを駆動する傾斜磁場電源109に接続され電流が供給される。具体的には、各傾斜磁場コイルの傾斜磁場電源109は、それぞれ後述の計測制御部111からの命令に従って駆動されて、それぞれの傾斜磁場コイルに電流を供給する。これにより、Ｘ、Ｙ、Ｚの３軸方向に傾斜磁場Ｇｘ、Ｇｙ、Ｇｚが発生する。この傾斜磁場コイル103と傾斜磁場電源109とを含めて傾斜磁場発生部となる。

【0022】

２次元スライス面の撮像時には、スライス面（撮像断面）に直交する方向にスライス傾斜磁場パルス（Ｇｓ）が印加されて被検体101に対するスライス面が設定され、そのスライス面に直交して且つ互いに直交する残りの２つの方向に位相エンコード傾斜磁場パルス（Ｇｐ）と周波数エンコード（読み出し）傾斜磁場パルス（Ｇｆ）が印加されて、核磁気共鳴信号（エコー信号）にそれぞれの方向の位置情報がエンコードされる。

【0023】

ＲＦ送信コイル104は、被検体101にＲＦパルスを照射するコイルであり、ＲＦ送信部110に接続され高周波パルス電流が供給される。これにより、被検体101の生体組織を構成する原子のスピンにＮＭＲ現象が誘起される。具体的には、ＲＦ送信部110が、後述の計測制御部111からの命令に従って駆動されて、高周波パルスを振幅変調し、増幅した後に被検体101に近接して配置されたＲＦ送信コイル104に供給することにより、ＲＦパルスが被検体101に照射される。このＲＦ送信コイル104とＲＦ送信部110とを含めてＲＦパルス発生部となる。

【0024】

ＲＦ受信コイル105は、被検体101の生体組織を構成するスピンのＮＭＲ現象により放出されるエコー信号を受信するコイルであり、信号処理部107に接続されて受信したエコー信号が信号処理部107に送られる。

【0025】

信号処理部107は、ＲＦ受信コイル105で受信されたエコー信号の検出処理を行う。具体的には、後述の計測制御部111からの命令に従って、信号処理部107が、受信されたエコー信号を増幅し、直交位相検波により直交する二系統の信号に分割し、それぞれを所定数(例えば、１２８、２５６、５１２等)サンプリングし、各サンプリング信号をＡ／Ｄ変換してディジタル量に変換する。従って、エコー信号は所定数のサンプリングデータからなる時系列のデジタルデータ(以下、エコーデータという)として得られる。そして、信号処理部107は、エコーデータに対して各種処理を行い、処理したエコーデータを計測制御部111に送る。

【0026】

計測制御部111は、被検体101の断層画像の再構成に必要なエコーデータ収集のための種々の命令を、主に、傾斜磁場電源109と、ＲＦ送信部110と、信号処理部107に送信してこれらを制御する制御部である。具体的には、計測制御部111は、後述する全体制御部112の制御で動作し、ある所定のパルスシーケンスの制御データに基づいて、傾斜磁場電源109、ＲＦ送信部110及び信号処理部107を制御して、被検体101へのＲＦパルスの照射及び傾斜磁場パルスの印加と、被検体101からのエコー信号の検出と、を繰り返し実行し、被検体101の撮像領域についての画像の再構成に必要なエコーデータの収集を制御する。繰り返しの際には、２次元撮像の場合には位相エンコード傾斜磁場の印加量を、３次元撮像の場合には更にスライスエンコード傾斜磁場の印加量も、変えて行なう。位相エンコードの数は通常1枚の画像あたり１２８、２５６、５１２等の値が選ばれ、スライスエンコードの数は、通常１６、３２、６４等の値が選ばれる。これらの制御により信号処理部107からのエコーデータを全体制御部112に出力する。

【0027】

全体制御部112は、計測制御部111の制御、及び、各種データ処理と処理結果の表示及び保存等の制御を行うものであって、演算処理部（ＣＰＵ）114と、メモリー113と、磁気ディスク等の内部記憶部115と、外部ネットワークとのインターフェースを行うネットワークＩＦ116と、を有して成る。また、全体制御部112には、光ディスク等の外部記憶部117が接続されていても良い。

【0028】

具体的には、計測制御部111に撮像シーケンスの実行によりエコーデータを収集させ、計測制御部111からのエコーデータが入力されると、演算処理部114がそのエコーデータに印加されたエンコード情報に基づいて、メモリー113内のｋ空間に相当する領域に記憶させる。以下、エコーデータをｋ空間に配置する旨の記載は、エコーデータをメモリー113内のｋ空間に相当する領域に記憶させることを意味する。また、メモリー113内のｋ空間に相当する領域に記憶されたエコーデータ群をｋ空間データともいう。

【0029】

演算処理部114は、このｋ空間データに対して信号処理やフーリエ変換による画像再構成等の処理を実行し、その結果である被検体101の画像を、後述の表示・操作部118に表示させ、内部記憶部115や外部記憶部117に記録させたり、ネットワークＩＦ116を介して外部装置に転送したりする。

【0030】

表示・操作部118は、再構成された被検体101の画像を表示する表示部と、ＭＲＩ装置の各種制御情報や上記全体制御部112で行う処理の制御情報を入力するトラックボール又はマウス及びキーボード等の操作部と、から成る。この操作部は表示部に近接して配置され、操作者が表示部を見ながら操作部を介してインタラクティブにＭＲＩ装置の各種処理を制御する。

【0031】

現在ＭＲＩ装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を２次元もしくは３次元的に撮像する。

【0032】

図２に、本発明の対象とする複数のエコー信号を取得するＦＳＥ(Fast Spin Echo)法のパルスシーケンスを示す。
本図において、ＲＦ、Ｇｓ、Ｇｐ、Ｇｒ、Ｓｉｇｎａｌはそれぞれ、ＲＦパルスの印加、スライス選択傾斜磁場パルスの印加、位相エンコード傾斜磁場パルスの印加、周波数エンコード傾斜磁場パルスの印加、エコー信号の取得を示す軸を表す。また、ここでは、一例として、１回の９０°励起ＲＦパルス308について、エコー信号322群を取得するＦＳＥシーケンスを示す。

【0033】

ＦＳＥシーケンスでは、まず、撮影面内のスピンに高周波磁場を与える９０°励起ＲＦパルス308とともに、スライス選択傾斜磁場パルス309を印加する。その後、スピンをスライス面内で反転するための１８０°反転ＲＦパルス314を繰り返し印加する。上述のように、1回の９０°励起ＲＦパルス308についてエコー信号322群を取得するため、１８０°反転ＲＦパルス314を複数回印加する。そのため、Shot毎の再収束ＲＦパルスの照射位相を反転させることで、それぞれの再収束ＲＦパルスから生じるＦＩＤ信号に対し、１８０度反転した位相情報を付与できる。

【0034】

ＦＳＥシーケンスは、１回のＲＦパルスによる励起後に複数のエコーデータを計測し、１つのk空間中にそのエコーデータを配置するシーケンスである。計測されたエコーデータはエコー番号ごとに、計測されたエコー順と同じ順序でk空間に配置され充填される。

【0035】

＜磁場歪みによるアーチファクトの発生位置と撮影ＦＯＶから折り返し位置を特定＞
ショートガントリ型ＭＲＩは、静磁場均一空間、及び、傾斜磁場線形領域が狭まることで、磁場歪みによるアーチファクトが生じやすくなる問題がある。

【0036】

今、仮にＺ軸方向（被検体の体軸方向）にＦＯＶ中心から２５０ｍｍの位置に磁場歪みによる画像歪みが生じる領域がある事を想定すると、４００ｍｍ程度のＦＯＶサイズを使用するSpine撮影などでは１５０％程度ＦＯＶを内部的に拡張する折り返し除去機能で、磁場歪み領域を十分に拡張ＦＯＶ内に収める事ができるため、課題となるアーチファクトは問題にならない。

【0037】

しかし、KneeやAnkleなどの関節部撮影では、ＦＯＶサイズは２００ｍｍ以下となる事が多く、磁場歪み領域を含むように拡張ＦＯＶを広げるようとすると３００〜４００倍と大幅に広げる必要があるため、現実的な方法ではない。

【0038】

図３（ａ）は、ＦＯＶ１５０ｍｍのKnee SAG撮影において、Body側Ｚ＝２５０ｍｍの位置に画像歪みによって高輝度化して輝点となったアーチファクトが発生する例を示している。

【0039】

内部的に２００％に拡張させる折り返し除去を使用することを想定すると、ＦＯＶは３００ｍｍとなり、画像歪みはＺ＝２５０ｍｍの位置にあることから、この輝点アーチファクトは、図３（ｂ）に示したようにＦＯＶ中心から−５０ｍｍの位置に折り返しが生じる事になる。

【0040】

折り返し位置が表示ＦＯＶ内に表示されるかが否かを判定する必要があるため、判定方法について次に説明する。
撮影ＦＯＶ内のアーチファクトの発生位置［ｍｍ］をＺ_alias、内部的に拡張した後のＦＯＶサイズ［ｍｍ］をFOV、Ｚ_aliasをFOVで割った小数点以下の剰余をR_alias、ＦＯＶ中心のＺ座標［ｍｍ］をZ_FOVとして、(Z_FOV＋ FOV/2)をFOVで割った小数点以下の剰余をR_FOVとし、判定値ｄを（式１）と定義する。
d ＝ R_alias − R_FOV ………（式１）
拡張率は２００％固定とするが、撮影におけるアーチファクトの発生位置が、拡張ＦＯＶ内か、表示ＦＯＶ内か以下の判定式から判断できる。
拡張ＦＯＶ内: 0 ＜ ｄ ＜ 0.25、 もしくは 0.75 ＜ ｄ ＜ 1.0 ………（式２）
表示ＦＯＶ内: 0.25 ≦ ｄ ≦ 0.75 ………（式３）

【0041】

図３（ａ）の例では、Z_alias＝250、Z_FOV＝0、FOV＝300であるため、R_alias＝0.833、R_FOV＝0.5、d＝0.333となるため、表示FOV内に発生する事になる。FOVをZ方向75mmの位置にオフセンターさせた撮影の場合には、Z_FOV ＝75、R_alias＝0.833、R_FOV＝0.75、d＝0.083となるため、拡張FOV内にアーチファクトが発生する事になる。

【0042】

表示ＦＯＶ内にアーチファクトが発生する場合に、関心領域にアーチファクトが重畳されて、診断の妨げとなる場合があるため、この状態を回避するためにアーチファクトの位置をＦＯＶ／２だけ移動して、拡張ＦＯＶ内に来る状態の説明図を図３（ｃ）に示す。
この図から、図３（ｂ）で折り返し位置に現われるアーチファクトが、図３（ｃ）では、矢印で示すように、上方に移動して拡張ＦＯＶ内に来ることが分かる。

【0043】

尚、磁場歪み（画像歪み）の位置は、装置据付時に十分に大きなファントムと折り返しの生じない十分に大きな５００ｍｍ以上のＦＯＶサイズで事前計測を行い、磁場歪みの座標情報をシステム情報として記憶しておく。事前計測はFast Spin Echoシーケンスを用いるのが望ましいが、その限りではない。

【0044】

＜折り返しアーチファクトの移動＞
ＦＩＤによる折り返しアーチファクトの位置をＦＯＶ／２だけ移動するためには、ｋ空間上でky(位相エンコード方向) １step毎に、１８０度の反転位相を付与する必要がある（図４）。図４は、ｋ空間上のＦＩＤ信号の位相を示し、step毎に交互に位相が反転した状態を示している。

【0045】

ｋ空間の充填方法は２つに分類され、上側から順に（もしくは下側から順に）充填していくSequential計測と、中心から上側に向かって充填後、中心から下側に充填（もしくは、中心から下側に充填後、中心から上側に充填）していくCentric計測がある。各充填方法では、Segment分割数（Shot数）に応じて、再収束ＲＦパルスの照射位相を最適化させる必要がある。

【0046】

図５にSequential計測時のk空間上のエコー配列を示す。本図（ａ）は、Shot数が奇数（例、Shot数3、Echo数4）の場合を示し、（ｂ）は、Shot数が偶数（例、Shot数2、Echo数6）の場合を示す。本図では、エコー配列は、下側から順に充填される場合を図示している。このエコー配列において、ky １step毎にＦＩＤ信号に反転位相を付与するためのＲＦ照射位相情報の詳細は、図９に示す。

【0047】

図９に示す表において、励起ＲＦパルスの照射位相は、図２で示すパルス３０８の条件を示し、再収束ＲＦパルスの照射位相は、図２で示すパルス314の条件を示す。また、表中のx,yはそれぞれx方向へのパルス印加、y方向へのパルス印加を示す。図１０に示す表も同様である。

【0048】

総Shot数が偶数個の場合は、偶数番目のShotの再収束ＲＦパルスは常に照射位相を反転させる。一方で、総Shot数が奇数個の場合、再収束ＲＦパルスの照射位相はパルス毎に反転させ、かつ、偶数番目のShotは、奇数番目のShotに対して１８０度反転するように設定する。

【0049】

次に、Centric計測時のk空間上のエコー配列を図６に示す。本図（ａ）は、Shot数/2が奇数（例、Shot数6、Echo数2）の場合を示し、（ｂ）は、Shot数/2が偶数（例、Shot数4、Echo数3）の場合を示す。本図では、k空間を半分にして、計測前半に上側半分を充填した後、下側半分を充填する配列を示している。このエコー配列において、ky １step毎にＦＩＤ信号の位相を反転させるためのＲＦ照射位相情報は、図１０に示す。

【0050】

総Shot数を２で割った商が偶数の場合、k空間の上側を充填する前半Shot群は、偶数番目のShotにおける再収束ＲＦパルスの照射位相を反転させ、後半の下側半分を充填するShot群では、奇数番目のShotにおける再収束ＲＦパルスの照射位相を反転させる。

【0051】

一方で、総Shot数を2で割った商が奇数の場合は、k空間の上側を充填する前半のShot群は、再収束ＲＦパルスの照射位相をパルス毎に反転させ、かつ、奇数番目のShotと偶数番目のShotでは照射位相が反転するようにする。更に、後半の下側半分を充填するShot群では、前半のShot群の照射位相を全反転したものになる。

【0052】

このように再収束ＲＦパルスの照射位相を変化させる事で、ＦＩＤ信号によるアーチファクトの位置をＦＯＶ／２だけ移動させる事ができる。

【0053】

＜偶数回加算撮影によるアーチファクトの信号抑制＞
また、偶数回の加算撮影を行う場合、前記全てのエコー配列に共通して、偶数加算目の撮影ではＲＦ照射位相が、奇数加算目のＲＦ照射位相に対して１８０度反転させる事でアーチファクト自体を消失させることができる。図７に示したようにk空間の同一kyラインにおいて、ＦＩＤの位相が奇数加算目と偶数加算目の撮影で位相反転した状態になるため、複素加算することで、アーチファクトを発生させる信号を消失させる事ができる。

【0054】

本実施例によれば、従来技術では判断できなかった磁場歪により発生するアーチファクトとなる輝点の位置の判断が可能となる。特に、マルチエコーを発生させるFast Spin Echoなどでは、輝点の輝度がより高くなるので、精度のよい画像が取得できる。

【実施例2】

【0055】

図８に、アーチファクト回避のためのＲＦ照射位相計算のフローチャートを示す。本フローチャートは、アーチファクトの発生位置を自動判定する処理の場合を示す。

【0056】

まず、磁場歪みの位置、ＦＯＶの位置、およびＦＯＶのサイズから式１を用いて判定値ｄを算出する（801）。
ここで、磁場歪みの位置は、事前に計測を行い、その計測結果は、例えば、メモリー113に保管しておく。ＦＯＶの位置、およびＦＯＶのサイズはユーザーが随時設定する。

【0057】

次に、判定値ｄが、0 ＜ d ＜ 0.25 、もしくは0.75 ＜ d ＜ 1.0のいずれかの範囲に入っているか判定する（802）。この範囲にあれば、ステップ803へ進み、範囲内になければ、ＲＦ照射位相計算は終了する。

【0058】

次に、Sequential計測か否かを判定(803)し、YESならば、Shot数が奇数か否か判定し（804）、図９（ａ）ＲＦ照射位相(805a)、又は図９（ｂ）ＲＦ照射位相(805b)の設定処理を行う。
NOならば、Shot数/2が奇数か否か判定し（808）、図１０（ａ）ＲＦ照射位相(805c)、又は図１０（ｂ）ＲＦ照射位相805d)の設定処理を行う。

【0059】

次に、ステップ805a〜805dのいずれかの処理において、偶数番目の加算か否か判定(806)し、YESならば、再収束ＲＦのＲＦ照射位相を全反転（807）し、NOならばＲＦ照射位相計算は終了する。
なお、上記の自動判定の処理は、演算処理装置（ＣＰＵ）１１４で実行される。

【0060】

上記ＲＦ照射位相計算では、アーチファクトの発生位置を自動判定で処理する場合を示したが、ユーザー判断ですることも可能である。例えば、アーチファクトが関心領域に重なって診断の妨げとなる場合に、撮影条件上でアーチファクト移動処理を発動するパラメータを選択する事も可能である。

【0061】

すなわち、本ＲＦ照射位相計算のフローチャートは、ステップ802の判定値dによる判定を自動で設定する手順を示しているが、手動設定することも可能である。その場合、上記フローチャートにおいてステップ802の処理後に判定値ｄを表示するステップ（表示せず）を設け、表示部（例えば、図１の118など）に判定値を表示する。ユーザーは表示部に表示された判定値と撮像画像上の関心領域とを参照しながら、フローチャートの流れをYESのフローにするかNOのフローにするかを入力部に設けた設定ボタンを操作することで決定することが可能となる。

【0062】

本実施例によれば、アーチファクトの発生位置を自動で判定でき、効率良く画像を取得可能である。また、手動モードを用いることにより、アーチファクトが関心領域にあるか判断しながら、現場の状況に応じた処理が可能となる。

【符号の説明】

【0063】

101：被検体、102：静磁場発生磁石、103：傾斜磁場コイル、104：ＲＦ送信コイル、105：ＲＦ受信コイル、106：ベッド、107：信号処理部、109：傾斜磁場電源、110：ＲＦ送信部、111：計測制御部、112：全体制御部、113：メモリー、114：演算処理部（ＣＰＵ）、115：内部記憶部、116：ネットワークＩＦ、117：外部記憶部、118：表示・操作部。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】