(57)【要約】【課題】ＭＰＧパルス印加を含む撮像を行うＭＲＩ装置において、ＴＥの延長を抑えてリアルタイムで高画質な拡散強調画像を得る。【解決手段】励起ＲＦパルス印加後にＭＰＧパルスを印加し、励起ＲＦパルスから所定時間後にエコー信号を収集する磁気共鳴イメージング装置において、ＭＰＧパルス印加後にエコー信号収集までの間に、少なくともスライス方向のナビゲーションエコーを取得し、取得したナビゲーションエコーから位相補正量を算出し、算出した位相補正量に基き、ＭＰＧパルスに含まれる位相誤差を補正する補正傾斜磁場を印加する。補正傾斜磁場印加後にエコー時間を計測する。ＳＮＲに影響のあるスライス方向の補正を行うことで、目的とする拡散運動を高ＳＮＲで画像化できる。【選択図】図３
2018015196.tif 000002
【特許請求の範囲】
【請求項１】
所定のパルスシーケンスに従って、エコー信号を収集する撮像部と、
前記撮像部を制御するシーケンス制御部と、
前記撮像部が収集したエコー信号を用いて演算を行う演算部と、を備え、
前記パルスシーケンスは、励起ＲＦパルス印加とエコー信号収集との間に、ＭＰＧパルス印加を含み、
前記シーケンス制御部は、前記ＭＰＧ印加後エコー信号収集までの間に、少なくともスライス方向について、ナビゲーションエコーの読み出し及び補正傾斜磁場印加を含む補正シーケンスを追加する制御を行い、
前記演算部は、前記補正シーケンスにおける前記ナビゲーションエコーの読み出しと前記補正傾斜磁場の印加との間で、前記ナビゲーションエコーを用いて前記補正傾斜磁場の印加量を算出し、算出した補正傾斜磁場印加量を前記シーケンス制御部に渡す補正パルス算出部を備え、
前記シーケンス制御部は、前記補正パルス算出部から受け取った前記補正傾斜磁場の印加量で前記補正傾斜磁場印加を行うことを特徴とする磁気共鳴イメージング装置。
【請求項２】
請求項１に記載の磁気共鳴イメージング装置であって、
前記補正シーケンスは、２軸方向又は３軸方向のナビゲーションエコーの読み出しと各軸方向の補正傾斜磁場の印加とを含み、
前記補正パルス算出部は、前記２軸方向又は３軸方向のそれぞれについて、補正傾斜磁場の印加量を算出することを特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項２に記載の磁気共鳴イメージング装置であって、
前記２軸方向又は３軸方向のナビゲーションエコーは、スライス方向のナビゲーションエコー、及び、リードアウト方向及び位相エンコード方向の少なくとも一方のナビゲーションエコーであることを特徴とする磁気共鳴イメージング装置。
【請求項４】
請求項３に記載の磁気共鳴イメージング装置であって、
前記補正パルス算出部は、前記リードアウト方向及び位相エンコード方向の少なくとも一方について、前記ナビゲーションエコーと基準となるナビゲーションエコーとの位相差分を用いて前記補正傾斜磁場の印加量を算出することを特徴とする磁気共鳴イメージング装置。
【請求項５】
請求項４に記載の磁気共鳴イメージング装置であって、
前記パルスシーケンスは、ＲＦ印加から信号収集までを少なくとも２回繰り返すパルスシーケンスであって、
前記補正パルス算出部は、１回目の繰り返しで収集したナビゲーションエコーを前記基準となるナビゲーションエコーとして用いることを特徴とする磁気共鳴イメージング装置。
【請求項６】
請求項１に記載の磁気共鳴イメージング装置であって、
前記補正パルス算出部は、前記ナビゲーションエコーの位相の傾きを算出し、当該傾きをゼロとする位相補正量を算出し、前記位相補正量から前記補正傾斜磁場の印加量を算出することを特徴とする磁気共鳴イメージング装置。
【請求項７】
請求項１に記載の磁気共鳴イメージング装置であって、
前記補正シーケンスは、２軸方向又は３軸方向のナビゲーションエコーの収集を含み、
前記演算部は、前記２軸方向又は３軸方向のうちスライス方向を除く方向について、それぞれ、前記ナビゲーションエコーと基準エコーとの位相差分を用いて位相補正量を算出し、当該位相補正量を用いて、前記エコー信号のｋ空間座標への配置位置を補正することを特徴とする磁気共鳴イメージング装置。
【請求項８】
請求項７に記載の磁気共鳴イメージング装置であって、
前記演算部は、前記ナビゲーションエコーと基準エコーとの位相差分を用いて、位相のオフセット量及び傾きを算出し、当該オフセット量を用いて、前記エコー信号の位相を補正することを特徴とする磁気共鳴イメージング装置。
【請求項９】
請求項１に記載の磁気共鳴イメージング装置であって、
前記補正シーケンスは、ナビゲーションエコーを発生させる傾斜磁場印加を含み、当該傾斜磁場の波形の一部は、前記ＭＰＧパルスの波形の一部と時間的に重複していることを特徴とする磁気共鳴イメージング装置。
【請求項１０】
請求項１に記載の磁気共鳴イメージング装置であって、
前記補正シーケンスは、ナビゲーションエコーを発生させる２以上傾斜磁場印加を含み、当該傾斜磁場の波形の一部は、隣接する傾斜磁場の波形の一部と時間的に重複していることを特徴とする磁気共鳴イメージング装置。
【請求項１１】
請求項１に記載の磁気共鳴イメージング装置であって、
前記パルスシーケンスにおける前記補正傾斜磁場の印加時間は固定されており、前記シーケンス制御部は、前記補正パルス算出部から受け取った印加量に対応して前記補正傾斜磁場の強度を変化させることを特徴とする磁気共鳴イメージング装置。
【請求項１２】
請求項１に記載の磁気共鳴イメージング装置であって、
前記シーケンス制御部は、前記励起ＲＦパルス印加と前記ＭＰＧパルス印加との間に、基準ナビゲーションエコーの読み出しを追加する制御を行い、
前記補正パルス算出部は、前記スライス方向を含む１軸以上の方向について、前記ナビゲーションエコーと前記基準ナビゲーションエコーとの位相差分を用いて、前記補正傾斜磁場の印加量を算出することを特徴とする磁気共鳴イメージング装置。
【請求項１３】
所定のパルスシーケンスに従って、エコー信号を収集する撮像部と、
前記撮像部を制御するシーケンス制御部と、
前記撮像部が収集したエコー信号を用いて演算を行う演算部と、を備え、
前記パルスシーケンスは、第一のＲＦパルス及び第二のＲＦパルスの印加と、前記第二のＲＦパルスに対し時間的に前後に位置する第一及び第二のＭＰＧパルスの印加と、前記第一及び第二のＭＰＧパルスの位相誤差を補正する補正シーケンスを含み、
前記補正シーケンスは、前記第一のＲＦパルスと前記第一のＭＰＧパルスの間に実行される第一のナビゲーションエコーの読み出しと、前記第二のＭＰＧパルスに続いて実行される第二のナビゲーションエコーの読み出しと、補正傾斜磁場の印加とを含み、
前記演算部は、前記第一のナビゲーションエコーと前記第二のナビゲーションエコーを用いて、前記補正傾斜磁場の印加量を算出し、前記シーケンス制御部に送り、
前記シーケンス制御部は、前記演算部から受け取った前記補正傾斜磁場の印加量で前記補正傾斜磁場印加を行った後、前記信号収集を行うことを特徴とする磁気共鳴イメージング装置。
【請求項１４】
励起ＲＦパルス印加後にＭＰＧパルスを印加し、励起ＲＦパルスから所定時間後にエコー信号を収集する磁気共鳴イメージング装置の動作を制御する方法であって、
ＭＰＧパルス印加後にエコー信号収集までの間に、少なくともスライス方向のナビゲーションエコーを取得し、取得したナビゲーションエコーから位相補正量を算出し、算出した位相補正量に基き、ＭＰＧパルスに含まれる位相誤差を補正する補正傾斜磁場を印加する動作を行わせる磁気共鳴イメージング装置の制御方法。
【請求項１５】
請求項１４記載の磁気共鳴イメージング装置の制御方法であって、
前記ナビゲーションエコーの取得は、同じ方向の複数のナビゲーションエコーの取得を含み、
前記位相補正量の算出処理は、同じ方向の複数のナビゲーションエコーのうち、最初に取得したナビゲーションエコーを基準ナビゲーションとし、当該基準ナビゲーションエコーとそれ以外のナビゲーションエコーとの位相差分を算出する処理と、当該位相差分について、傾きを算出する処理とを含むことを特徴とする磁気共鳴イメージング装置の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体中の水素や燐等からの核磁気共鳴信号（以下、ＮＭＲ信号と呼ぶ）を測定し、核の密度分布や緩和時間分布等を画像化する磁気共鳴イメージング（以下、ＭＲＩと呼ぶ）装置に関し、特に、拡散強調画像の測定技術に関する。
【背景技術】
【０００２】
ＭＲＩ装置で得られる画像の一つに、拡散強調画像（ＤＷＩ：Ｄｉｆｆｕｓｉｏｎ  Ｗｅｉｇｈｔｅｄ  Ｉｍａｇｅ）がある。拡散強調画像は、ＭＰＧ：Ｍｏｔｉｏｎ  Ｐｒｏｂｉｎｇ  Ｇｒａｄｉｅｎｔと呼ばれる強度の高い傾斜磁場パルスを印加することで、水分子の拡散運動を画像のコントラストに反映させた画像である。
【０００３】
ＭＰＧパルスは、励起ＲＦパルスとして１８０°ＲＦパルスを用いる場合には、その前後に２つ印加される。２つのＭＰＧパルスは、その極性は異なるが、印加量（面積）は等しい。このため、空間位置が移動していないスピン（静止スピン）については、最初に印加されるＭＰＧパルスにより変化した位相は、後に印加されるＭＰＧパルスにより戻される。一方、空間位置がＭＰＧパルスの印加方向に移動したスピンは、各ＭＰＧパルスにより作られる位相の変化量が一致しないことから、後のＭＰＧ印加時点で位相がゼロにならず、周囲の静止スピンとの位相の不一致によりエコー信号が減衰し、移動スピンの部分、例えば血流などは画像上に低信号として表現されることになる。拡散強調イメージングは、この静止スピンと動きのあるスピンとのコントラスト差を利用して拡散運動を画像化する手法である。
【０００４】
このように拡張強調イメージングは、二つのＭＰＧパルスの間に生じた拡散運動を可視化することを目的としているが、拡散運動以外の運動もまた画像のコントラストに反映されてしまうという問題がある。例えば測定対象が人体の場合は、拍動や呼吸といった生理的な動きが問題となる。このような生理的な動きは一般的に観察対象では無いため、拡散強調画質の劣化として認識される。
【０００５】
この問題を解決する１つの方法として、リアルタイムにスピンに生じた位相誤差のマップを測定し、この位相誤差マップを用いて計測されるエコーを補正する方法がある（特許文献１）。この方法では、ＭＰＧパルス印加後に２Ｄもしくは３Ｄの空間を走査する傾斜磁場を印加し、２Ｄもしくは３Ｄの位相誤差マップを測定する。測定した位相誤差マップを用いて、位相誤差をキャンセルするＲＦパルスと傾斜磁場を計算し、エコー信号を計測する前に照射及び印加する。この他に類似の方法として、ＭＰＧパルス印加後に画像面内の方向、すなわち周波数エンコード方向及び位相エンコード方向に読み出し傾斜磁場パルスを印加し、得られたエコー信号から算出される０次及び１次の位相誤差に応じて傾斜磁場パルスと補償Ｂ０パルスを印加することで０次及び１次の位相誤差を補償する手法が提案されている（非特許文献１）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許２０１３／０２２９１７７明細書
【非特許文献】
【０００７】
【非特許文献１】Ｄａｖｉｄ  Ｇ．  Ｎｏｒｒｉｓ  ａｎｄ  Ｗｏｌｆｇａｎｇ  Ｄｒｉｅｓｅｌ，  Ｏｎｌｉｎｅ  ｍｏｔｉｏｎ  ｃｏｒｒｅｃｔｉｏｎ  ｆｏｒ  ｄｉｆｆｕｓｉｏｎ−ｗｅｉｇｈｔｅｄ  ｉｍａｇｉｎｇ  ｕｓｉｎｇ  ｎａｖｉｇａｔｏｒ  ｅｃｈｏｅｓ：Ａｐｐｌｉｃａｔｉｏｎ  ｔｏ  ＲＡＲＥ  ｉｍａｇｉｎｇ  ｗｉｔｈｏｕｔ  ｓｅｎｓｉｔｉｖｉｔｙ  ｌｏｗｗ，  Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  ｉｎ  Ｍｅｄｉｃｉｎｅ  ４５：  ７２９−７３３  （２００１）
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、上述した方法には以下の問題がある。まず、特許文献１の方法では、位相誤差マップを２Ｄもしくは３Ｄで測定するため、測定に要する時間が長い。特に３Ｄでの測定には長い時間を要する。得られた位相誤差マップを用いた補正は、エコー信号計測までに行う必要があるため、位相誤差マップの計測時間の延長は、エコー信号のＴＥ：Ｅｃｈｏ  Ｔｉｍｅの延長につながり、結果として画像のＳＮＲが低下する。位相誤差マップを２Ｄで測定する場合は３Ｄに比べると計測時間は短くなるが、測定していない軸に関しては動きの影響を観察できないため補正できない。具体的にはＸ−Ｙ平面の位相誤差マップを測定する場合、Ｚ軸方向の動きの影響を補正することができない。また、特許文献１の方法では、位相誤差マップから得られる位相誤差を補正するためにＲＦパルスを用いている。位相誤差をキャンセルするＲＦパルス及びそれに対応する傾斜磁場の計算にはフーリエ変換、もしくはフリップ角が大きい場合にはブロッホ方程式を解く必要があり処理コストが大きい。加えて、計算して得られた位相誤差をキャンセルするＲＦパルスを印加する時間も数十ｍｓのオーダーで要する。ＳＮＲの低下を考慮すると、計測条件やＭＲＩ装置にも依るがＴＥの延長は１０ｍｓ程度以下が望ましい。
【０００９】
次に、非特許文献１の方法では、傾斜磁場パルスのみでＭＰＧ印加中の生理的な動きの影響を補正するため位相誤差マップを作成する場合に比べて短い処理時間で実行することができる。しかしながら、非特許文献１の方法では、ＭＰＧパルスの印加軸をスライス方向に平行な方向にして、それと直交する二方向についてナビゲーターエコーによる位相誤差補正を行っているので、スライス方向に対する生理的な動きの影響は補正されずに残る。一般に、生理的な動きに起因する位相誤差は断面内の二方向に限らずスライス方向にも生じるが、スライス方向の位相誤差は画像のＳＮＲに与える影響が大きく、エコー信号取得後に補正することはできない。また拡散強調イメージングではＭＰＧパルスの印加軸を種々に異ならせて撮像が行われ、位相誤差を生じやすい方向とＭＰＧパルスの印加軸とに規則的な関係を見出すことはできないため、位相誤差補正を断面内の二方向においてのみ行うことを前提とした非特許文献１の方法には限界がある。
【００１０】
本発明は、上記問題点に鑑みてなされたものであり、ＭＰＧパルス印加を含む撮像を行うＭＲＩ装置において、ＴＥの延長を抑えてリアルタイムで高画質な拡散強調画像を得ることが可能なＭＲＩ装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上述した問題点を解決するため、本発明のＭＲＩ装置は、ＭＰＧパルス印加後エコー信号計測までの間で、ＭＰＧパルス印加後のエコー信号のｋ空間におけるシフト量を、少なくともスライス選択傾斜磁場方向（以下、スライス方向という）を含む軸方向に計測する手段と、前記シフト量を補正するために必要な傾斜磁場パルスを計算する手段と、前記補正用傾斜磁場パルスを印加する手段とを備える。
【００１２】
具体的には、本発明のＭＲＩ装置は、所定のパルスシーケンスに従って、エコー信号を収集する撮像部と、前記撮像部を制御するシーケンス制御部と、前記撮像部が収集したエコー信号を用いて演算を行う演算部と、を備え、前記パルスシーケンスは、励起ＲＦパルス印加とエコー信号収集との間に、ＭＰＧパルス印加を含む。前記シーケンス制御部は、前記ＭＰＧ印加後エコー信号収集までの間に、少なくともスライス方向について、ナビゲーションエコーの読み出し及び補正傾斜磁場印加を含む補正シーケンスを追加する制御を行い、前記演算部は、前記補正シーケンスにおける前記ナビゲーションエコーの読み出しと前記補正傾斜磁場の印加との間で、前記ナビゲーションエコーを用いて前記補正傾斜磁場の印加量を算出し、算出した補正傾斜磁場印加量を前記シーケンス制御部に渡す補正パルス算出部を備え、前記シーケンス制御部は、前記補正パルス算出部から受け取った前記補正傾斜磁場の印加量で前記補正傾斜磁場印加を行う。
【発明の効果】
【００１３】
本発明によれば、スライス方向を含むＭＰＧ印加後のエコー信号のｋ空間におけるシフト量を算出し、その値に基づいて前記シフト量を補正するための傾斜磁場を印加することにより、生理的な運動に起因して撮像対象に生じる大域的な運動による画像コントラストの変化を抑制することができ、かつそれら一連の処理を従来技術よりも短時間に実施できる。これによりＴＥの延長すなわちＳＮＲの低下が抑えられ、本来の観察対象である局所的な水分子の拡散運動を観察しやすい画像を提供できる。
【図面の簡単な説明】
【００１４】
【図１】本発明に係るＭＲＩ装置の一実施形態の全体構成を示すブロック図。
【図２】一般的なＤＷＩパルスシーケンスの一例を示す図。
【図３】第一実施形態のシーケンサ及び演算部の構成と処理の概要を示す図。
【図４】補正シーケンスを含むＤＷＩパルスシーケンスの一例を示す図。
【図５】第一実施形態の演算部の処理の一部を示すフロー図。
【図６】第一実施形態の演算部の処理の一部を示すフロー図。
【図７】図４の補正シーケンスの詳細を示す図。
【図８】図４の補正シーケンスの詳細を示す図。
【図９】第二実施形態のシーケンサ及び演算部の構成と処理の概要を示す図。
【図１０】第三実施形態のシーケンサ及び演算部の構成と処理の概要を示す図。
【図１１】第三実施形態で用いる補正シーケンスを含むＤＷＩパルスシーケンスの一例を示す図。
【図１２】第三実施形態の演算部の処理を示すフロー図。
【発明を実施するための形態】
【００１５】
以下、添付図面に従って本発明のＭＲＩ装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
【００１６】
最初に、本発明が適用されるＭＲＩ装置の一例の全体概要を図１に基づいて説明する。ＭＲＩ装置は、ＮＭＲ現象を利用して被検体の断層画像を得るものであり、図１に示すように、静磁場発生部２と、傾斜磁場発生部３と、送信部５と、受信部６と、信号処理部７と、シーケンサ４とを備えている。本明細書では、静磁場発生部２、傾斜磁場発生部３、送信部５及び受信部６をまとめて撮像部という。
【００１７】
静磁場発生部２は、被検体1の周りに配置された、永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源からなる。垂直磁場方式の静磁場発生源は、被検体1の周りの空間にその体軸と直交する方向に、水平磁場方式の静磁場発生源は、体軸方向に均一な静磁場を発生させる。
【００１８】
傾斜磁場発生部３は、ＭＲＩ装置の座標系（静止座標系）であるＸ、Ｙ、Ｚの３軸方向に巻かれた傾斜磁場コイル９と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源１０とから成り、後述のシ−ケンサ４からの命令に従ってそれぞれのコイルの傾斜磁場電源１０を駆動することにより、Ｘ、Ｙ、Ｚの３軸方向に傾斜磁場を印加する。３軸方向の傾斜磁場の組み合わせにより、任意の方向の傾斜磁場を形成することができる。撮影時には、スライス面（撮影断面）に直交する方向にスライス方向傾斜磁場パルスを印加して被検体１に対するスライス面を設定し、そのスライス面に直交して且つ互いに直交する残りの２つの方向に位相エンコード方向傾斜磁場パルスと周波数エンコード方向傾斜磁場パルスを印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。また所定方向の傾斜磁場パルスを印加することで、被検体１の組織を構成するスピンにその印加方向に沿って１次の位相変化を与えることができる。
【００１９】
送信部５は、被検体１の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体１に高周波磁場パルス（ＲＦパルスという）を照射するもので、高周波発振器１１と変調器１２と高周波増幅器１３と送信側の高周波コイル（送信コイル）１４ａとから成る。高周波発振器１１から出力された高周波磁場パルスをシーケンサ４からの指令によるタイミングで変調器１２により振幅変調し、この振幅変調された高周波パルスを高周波増幅器１３で増幅した後に被検体１に近接して配置された高周波コイル１４ａに供給することにより、ＲＦパルスが被検体１に照射される。
【００２０】
受信部６は、被検体１の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号（ＮＭＲ信号）を検出するもので、受信側の高周波コイル（受信コイル）１４ｂと、信号増幅器１５と、直交位相検波器１６と、Ａ／Ｄ変換器１７とから成る。送信側の高周波コイル１４ａから照射された電磁波によって誘起された被検体１の応答のＮＭＲ信号が被検体１に近接して配置された高周波コイル１４ｂで検出され、信号増幅器１５で増幅された後、シーケンサ４からの指令によるタイミングで直交位相検波器１６により直交する二系統の信号に分割され、それぞれがＡ／Ｄ変換器１７でディジタル量に変換されて、信号処理部７に送られる。
【００２１】
信号処理部７は、各種データ処理と処理結果の表示及び保存等を行うもので、種々の演算や制御を行うディジタル信号処理装置８と、光ディスク１９、磁気ディスク１８等の外部記憶装置と、ＲＯＭ２１、ＲＡＭ２２等の内部記憶媒体と、ディスプレイ２０とを有する。受信部６からのデータがディジタル信号処理装置８に入力されると、ディジタル信号処理装置８が信号処理、画像再構成等の処理を実行し、その結果である被検体１の断層画像をディスプレイ２０に表示すると共に、外部記憶装置の磁気ディスク１８等に記録する。ディジタル信号処理装置８が行う演算や制御などの処理は、ＣＰＵとそれに実装されるソフトウェアで実現してもよいし、その一部はＡＳＩＣやＦＰＧＡなどのハードウェアで実現してもよい。
【００２２】
また信号処理部７には、ＭＲＩ装置の各種制御情報や信号処理部７で行う処理の制御情報を入力するための操作部２５が備えられる。操作部２５は、トラックボール又はマウス２３、キーボード２４などを備える。操作部２５は、ディスプレイ２０に近接して配置され、操作者がディスプレイ２０を見ながら操作部２５を通してインタラクティブにＭＲＩ装置の各種処理を制御する。
【００２３】
シーケンサ４は、ＲＦパルスと傾斜磁場パルスを所定のパルスシーケンスで繰り返し印加する制御手段（シーケンス制御部）で、ディジタル信号処理装置８の制御で動作し、被検体１の断層画像のデータ収集に必要な種々の命令を撮像部（送信部５、傾斜磁場発生部３、および受信部６）に送る。パルスシーケンスは、ＲＦパルスや傾斜磁場パルスの印加の強度やタイミング及びＮＭＲ信号（エコー信号）を収集するタイミングなどを定めたタイミングチャートであり、撮像方法によって種々のパルスシーケンスがある。これらパルスシーケンスはプログラムとして信号処理部７の記憶装置に予め格納されており、シーケンサ４が所定のパルスシーケンスと撮像パラメータとを読み込むことで実行される。本実施形態のＭＲＩ装置は、パルスシーケンスとしてＭＰＧパルスを用いた拡散強調イメージング（ＤＷＩ）パルスシーケンスを実行する。
【００２４】
図２に一般的なＤＷＩパルスシーケンスの一例を示す。このＤＷＩパルスシーケンス２００では、まずＲＦパルス２０１を照射し、同時に印加されるスライス選択傾斜磁場（Ｇｓ）で決まる特定の位置のスピンを励起する。その後、第一のＭＰＧパルス２０３を印加する。図２ではＭＰＧパルスをスライス選択傾斜磁場パルスの軸（Ｇｓ）に印加しているが、他軸もしくは複数軸に印加することもある。第一のＭＰＧパルス２０３を印加した後、スピンの位相を反転させるＲＦパルス２０２を照射する。その後、第二のＭＰＧパルス２０４を印加する。第二のＭＰＧパルス２０４は、第一のＭＰＧパルス２０３と面積（印加量）が等しく、また反転ＲＦパルス２０２の後に印加しているので第一のＭＰＧパルス２０３とは逆の極性の位相回転を与える。これにより、静止スピンについては第一のＭＰＧパルス２０３によって与えられる位相回転が、第二のＭＰＧパルス２０４によって巻き戻される。その後、位相エンコード傾斜磁場パルス２０５を印加し、読出し傾斜磁場パルス２０６を印加しながらエコー２０７を計測する。このようなパルスシーケンスにより、拡散運動する水分子と静止スピンとにコントラスト差を生じさせて拡散運動を画像化する。なお図では、一つのエコーのみを計測するシングルショット撮像の場合を示しているが、位相エンコード傾斜磁場パルスと読出し傾斜磁場パルスの印加を繰り返しながら複数のエコーを計測する場合（マルチショット撮像）もある。
【００２５】
本実施形態のＭＲＩ装置は、このようなＭＰＧパルスを含むパルスシーケンスを用いて撮像したときに、水分子の拡散運動とは別の被検体の拍動や呼吸などの動き（これら拡散運動と区別して大域的運動という）に起因する画像のコントラストの変化を抑制するために、大域的運動補正処理機能を備えることが特徴である。具体的には、第一のＭＰＧパルス２０３印加時と第二のＭＰＧパルス２０４印加時における大域的運動に起因して生じる傾斜磁場印加量の過不足量を検出してリアルタイムに補正する。
【００２６】
以下、大域的運動補正処理機能の具体的な実施形態を説明する。
＜第一実施形態＞
本実施形態のＭＲＩ装置は、撮像部が、ＤＷＩパルスシーケンスのＭＰＧパルス印加からエコー取得までの間に、少なくともスライス方向のナビゲーションエコー（以下、ナビエコーという）を発生させ、このナビエコーから計算された位相誤差に相当する印加量の傾斜磁場（補正傾斜磁場パルス）を印加するシーケンス（補正シーケンス）を追加して撮像すること、及び、ディジタル信号処理装置８（演算部）が、シーケンスの実行中に計測したナビエコーを用いて位相誤差或いはそれに相当する傾斜磁場印加量を算出し、シーケンス制御部（シーケンサ）に渡すこと、が特徴である。
【００２７】
本実施形態のＭＲＩ装置は、上述した一連の処理をリアルタイムに実施する。このため、シーケンサ４、送信部５、受信部６、信号処理部７にはリアルタイムシステムを採用する。例えば、シーケンサ４と信号処理部７（ディジタル信号処理装置８）にはＲＴＯＳ（Ｒｅａｌ−ｔｉｍｅ  ｏｐｅｒａｔｉｎｇ  ｓｙｓｔｅｍ）を搭載し、送信部５と受信部６はＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｂｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）で構成することができる。またディジタル信号処理装置８は、図３に示すように、位相誤差を算出するための補正パルス算出部８１を備える。補正パルス算出部８１は、ＲＴＯＳ上で動作するソフトウェアや専用のハードウェアで実現することができる。なおディジタル信号処理装置８には、エコー信号を用いて画像を再構成する画像再構成部や、拡散強調画像のための演算を行う演算部が備えられるがここでは図示を省略している。
【００２８】
本実施形態における処理の概要を図３のブロック図を参照して説明する。
シーケンサ４は、操作部２５を介してＤＷＩパルスシーケンスとその撮像パラメータが設定されると、パルスシーケンスを計算し、実行する。ここで撮像パラメータのうちＴＥ（エコー時間）については、補正シーケンスに要する時間を考慮した時間（最短の時間）が設定される。ＤＷＩパルスシーケンスは、例えば、図２に示したように二つのＭＰＧパルス２０３，２０４を含むパルスシーケンス２００である。
【００２９】
シーケンサ４は、さらにパルスシーケンス２００に挿入して実行される補正シーケンス４００を含む。補正シーケンス４００は、ナビエコー読み出し傾斜磁場印加３０１と補正パルス印加３０３とを含み、最初にナビエコー読み出し傾斜磁場印加処理３０１が実行される。ナビエコー読み出し傾斜磁場は、ＭＰＧパルスを印加後にＭＰＧパルス印加中の被検体１の動きを測るために印加される読み出し傾斜磁場であり、ＭＰＧパルス印加後に、少なくともスライス方向について印加される。ナビエコー読み出し傾斜磁場印加処理３０１により発生するエコー信号（ナビエコー）は受信部６を経てディジタル信号処理装置８の補正パルス算出部８１に渡される。
【００３０】
補正パルス算出部８１は、ナビエコーを用いて、被検体1の動きにより生じるＭＰＧパルスの過不足量と、それを補正するために必要な傾斜磁場印加量（補正パルス印加量）を算出し、その算出結果をシーケンサ４に通知する（補正パルス算出処理３０２）。
【００３１】
シーケンサ４は、補正パルス算出部８１が算出した補正パルス印加量に従い、補正シーケンス４００の後半である補正パルス印加処理３０３を実行する。即ち、算出された印加量でナビエコーと同軸に補正用傾斜磁場パルスを印加する。その後、ＤＷＩパルスシーケンス２００の残り部分を実行し、ＤＷＩのためのエコー２０７を収集する。
【００３２】
図４に、シーケンサ４が実行する、補正シーケンスが追加されたＤＷＩパルスシーケンス２００Ａの一例を示す。図４において、図２と同じ要素は同じ符号で示し、重複する説明は省略する。このＤＷＩパルスシーケンス２００Ａは、第二のＭＰＧパルス２０４印加の後、エコー２０７の収集までに、ＭＰＧ の過不足量を補正するための補正シーケンス４００が挿入される。図４に示す例では、ナビエコー読み出し傾斜磁場をスライス傾斜磁場（Ｇｓ）方向、位相エンコード傾斜磁場（Ｇｐ）方向、読み出し傾斜磁場（Ｇｒ）方向のそれぞれに順番に印加し、各軸のナビエコー４０１を取得する（ナビエコー読み出し傾斜磁場印加処理３０１）。最後のナビエコー４０１を取得から所定時間後に位相誤差補正のための傾斜磁場パルス４０３を印加する（補正パルス印加処理３０３）。ナビエコー４０１取得から補正用傾斜磁場パルス４０３までの時間４０２は、補正パルス算出処理３０２に要する時間であり、補正用傾斜磁場パルス４０３は、補正パルス算出処理３０２で算出された印加量で印加される。
【００３３】
補正用傾斜磁場パルス４０３印加後に、位相エンコードパルス２０５及び読出し傾斜磁場パルス２０６を印加し、エコー２０７を計測することは図２のパルスシーケンス２００と同様である。
【００３４】
以下、補正パルス算出部８１が行う処理３０２の詳細について説明する。
補正パルス算出部８１は、処理３０１で読み出したナビエコーを用いて、Ｇｓ、Ｇｐ、Ｇｒの３方向におけるＭＰＧパルスの過不足量と補正パルス印加量を算出する。ＭＰＧパルスの過不足量の算出方法はＧｓ方向とＧｐ、Ｇｒ方向で異なる。Ｇｓ方向におけるＭＰＧパルスの過不足量の算出の手順を図５のフローチャートを参照して説明する。
【００３５】
まず処理５０１にて、Ｇｓ方向にナビエコー読み出し傾斜磁場を印加して得られたナビエコーNavi_Echo(t)を読み込む。処理５０２において読み込んだナビエコーNavi_Echo(t)をフーリエ変換し、複素データであるNavi (x)を得る。処理５０３では複素信号Navi (x)の実部と虚部とを用いて、位相のデータ列Phase_scalar(x)を算出する。このように得られる位相は主値回りを生じているため、処理５０４で位相アンラップ処理を行う。具体的には、次式（１）に従い、Navi (x)の位相Phase_unwrapped(x)を算出する。
【００３６】
【数１】
2018015196.tif 000003
式（１）において、Phase_unwrapped(x)はスカラーの位相データ、Phase[ ]は複素データの位相値を返却する関数、Conjugate[ ]は複素共役処理を表す関数である。
【００３７】
処理５０４において、Phase_unwrapped(x)に対して最小二乗法（式（２））により１次直線を当てはめ、その傾きを算出する。
【数２】
2018015196.tif 000004
式（２）中、ＮはナビエコーNavi_Echo(t)のデータ数である（他の方向について式（２）を用いる場合は、当該方向のデータ数である。以下、同じ）。こうして算出された位相の傾きがＭＰＧパルスの過不足量を表す。
【００３８】
処理５０５では、この位相の傾きを用いて、Ｇｓ方向における補正パルス印加量[s・T/m]を次式（３）に従い算出する。
【数３】
2018015196.tif 000005
式（３）中、Durationはナビエコーのサンプリング時間［ｓ］であり、GcAmpはナビエコー読み出し傾斜磁場[T/m]である。
【００３９】
なお高周波コイル（受信コイル）１４ｂが複数のコイルで構成されている場合、ナビエコーはコイル毎に計測される。この場合、各コイルの形状及び位置により得られるナビエコーの位相の分布が異なるため、各コイルのナビエコー毎に式（２）に従って位相の傾きを算出し、その平均値を用いて式（３）により補正パルス印加量を決定する。
【００４０】
次に、Ｇｐ、Ｇｒ方向におけるＭＰＧパルスの過不足量の算出の手順を図６のフローチャートに従い説明する。Ｇｐ方向とＧｒ方向における補正パルス印加量の算出手順は等しいため、ここではＧｐ方向についてのみ説明する。また、Ｇｐ方向とＧｒ方向の補正は、マルチショット計測もしくは加算数が２以上の場合、すなわち１枚の画像に必要なエコー信号を、図４で示すパルスシーケンスを単位として複数回の繰り返して得る場合にのみ行う。
【００４１】
まず、処理６０１にて、Ｇｐ方向にナビエコー読み出し傾斜磁場を印加して得られたナビエコーNavi_Echo(t)を読み込む。読み込んだナビエコーNavi_Echo(t)を処理６０２においてフーリエ変換し、Navi (x)を得る。分岐処理６０３にて、処理中のナビエコーNavi(x)が繰り返しの1番目のエコーであるかを判断する。Navi (x)が最初（1番目）の繰り返しで得たエコーである場合、処理６０４においてNavi (x)をＲＡＭ２２に保存する。保存したNavi (x)は、以後の処理において基準ナビエコーNavi _std(x)として用いられる。一方、分岐処理６０３にて、処理中のナビエコーNavi(x)が１番目の繰り返しのエコーではないと判断した場合、処理６０５で当該ナビエコーNavi(x)と基準ナビエコーNavi _std(x)の位相差分を次式に従い計算する。
【００４２】
【数４】
2018015196.tif 000006
式（４）中、Navi_subtracted(x)は正極データと負極データを位相差分した複素データである。
【００４３】
その後は、Ｇｓ方向の処理（５０３〜５０５）と同様に、処理６０６で、Navi_subtracted(x)の位相を上述の式（１）を用いてアンラップ処理し、位相Phase_unwrapped(x)を算出する。また処理６０７で上述の式（２）から位相分布の傾きを求め、処理６０８にて上述の式（３）に従いＧｐ方向における補正パルス印加量を算出する。
【００４４】
Ｇｐ、Ｇｒ方向における処理とＧｓ方向における処理の違いは、基準ナビエコーを用いるか否かである。Ｇｓ方向の処理では、ナビエコーの位相の傾きがゼロになるように補正パルス印加量を定める。これは、Ｇｓ方向の位相の傾きがゼロとなるのが最もＳＮＲの高くなる条件であり、ＭＰＧパルス印加中の動きによって生じた誤差を含んだエコーを基準にするとＳＮＲが低下するためである。一方、Ｇｐ、Ｇｒ方向の処理では基準ナビエコーの位相分布と一致するよう補正パルス印加量を定めている。Ｇｐ、Ｇｒ方向は画像面内の方向であるため、画像を構成する複数のエコー信号間でＧｐ、Ｇｒ方向の位相の傾きを一致させれば、絶対量として位相の傾きを持っていたとしても、画像の絶対値には影響しないからである。
【００４５】
補正パルス算出部８１は、上述したように算出したＧｓ、Ｇｐ、Ｇｒの３方向におけるＭＰＧパルスの過不足量、即ち補正パルス印加量をシーケンサ４に通知する。
【００４６】
シーケンサ４は、補正パルス算出部８１が算出したＧｓ、Ｇｐ、Ｇｒの各方向における補正パルス印加量に従い傾斜磁場パルス（図４：４０３）を印加する（補正パルス印加処理３０３）。なお傾斜磁場パルスの印加量は、印加時間と磁場強度との積で決まる。式（３）のようにパルスの面積として補正パルス印加量を算出した場合、シーケンサ４は補正パルスの印加時間を固定しておき、補正パルス印加量に応じてパルスの強度を変化させればよい。補正パルスを印加することにより第一のＭＰＧパルス２０３印加時と第二のＭＰＧパルス２０４印加時における被検体１の大域的な位置の変化（大域的運動）に起因して生じる傾斜磁場印加量の過不足を補うことができる。従ってその後に計測するエコー信号は、過不足なくＭＰＧパルスが印加されたときに計測されるエコー信号と同等となり、ＳＮＲのよい拡散強調画像を得ることができる。
【００４７】
本実施形態のＭＲＩ装置によれば、画像のＳＮＲに最も影響があり且つエコー信号取得後には補正できないスライス方向の位相誤差をゼロとする補正を行うことでＳＮＲの低下を抑制し、本来の観察対象である拡散運動を高コントラストで描出することができる。またスライス方向を含む３方向について補正を行うことで、非特許文献１に記載される技術や２Ｄ位相マップを作成するだけでは測定できない全て軸について大域的運動に起因する位相誤差を補正することができる。
【００４８】
また本実施形態によれば位相誤差補正のための処理（補正シーケンス）がＲＦパルスの印加を含まず、読出し傾斜磁場パルスを用いたナビエコーの発生と補正傾斜磁場パルスの印加だけで構成され、各軸1つのエコー信号をフーリエ変換し、その位相の傾きを傾斜磁場の面積に換算するだけであるため計算時間を含む補正処理時間が短い。これによりＴＥの延長を大幅に抑制することができる。またＲＦパルスの計算コストが不要であり、処理コストが低い。
【００４９】
以下、本実施形態のＭＲＩ装置が、ＴＥの延長を抑制できること、即ち補正シーケンスがリアルタイム撮像に適用可能であることを具体的な例を用いて説明する。
【００５０】
ナビエコー読み出し傾斜磁場の詳細シーケンスを図7に示す。エコー信号のサンプリング間隔を１．５６μｓ、サンプリング点数を６４とすると１軸当たりのサンプリング時間は約１００μｓである。読み出し傾斜磁場の強度を９ｍＴ／ｍとし、傾斜磁場の最大Ｓｌｅｗ  Ｒａｔｅが１００Ｔ／ｍ／ｓのＭＲＩ装置であれば、１軸のナビエコーの計測に要する時間は、読み出し傾斜磁場のディフェイズ及びリフェイズパルスを含めて７００μｓ以下に収めることができる。図４のパルスシーケンス２００Ａように、３軸方向にナビエコー読み出し傾斜磁場を順番に印加する場合、図７に示すように傾斜磁場のプラトー部以外は他軸の傾斜磁場パルスを同時に印加してもよいため、３軸全ての読み出し傾斜磁場の印加時間は７００＋７００＋１００＝１５００μｓとなる。さらに図８に示すように、エコー信号をサンプリングしていない部分を直前のＭＰＧパルス２０４と重ねることで、更にナビエコー読み出し傾斜磁場印加処理３０１の時間を短縮することもできる。また、補正パルス算出処理３０２はエコー信号をサンプリングし終わった後、読み出し傾斜磁場の印加終了を待たずに始めることができるので、ここでも時間を短縮することができ、ナビエコー読み出し傾斜磁場印加処理３０１に掛かる実質的な時間は９００μｓとなる。
【００５１】
一方、補正パルス算出処理３０２については、ディジタル信号処理装置８の性能にもよるが、ＭＰＧパルスの過不足量の算出とシーケンサ４への通知に要する時間は３ｍｓあれば十分可能である。従って、一般的な性能を持つＭＲＩ装置ではナビエコーの読み出しから補正パルスの印加までを５ｍｓ程度で行うことができ、図４のパルスシーケンスではＴＥの延長を１０ｍｓに抑えることができる。
【００５２】
以上、第一実施形態について図面を参照して説明したが、本実施形態は、少なくともスライス傾斜磁場方向について大域的運動補正処理をリアルタイムで行う手段（シーケンス制御部及び補正パルス算出部）を備えることが特徴であり、それ以外の要素や手順については適宜変更が可能である。
例えば、パルスシーケンスはＭＰＧパルスを含むＤＷＩパルスシーケンスであれば、図２に示すパルスシーケンスに限らず適用することができる。また、大域的運動補正処理におけるナビエコー及び補正パルスの印加方向は、Ｇｓ、Ｇｐ、Ｇｒ方向でなくとも、装置のＸ、Ｙ、Ｚ方向としてもよいし、３軸以上の方向に印加して補正を行ってもよい。またＧｓ方向のみでナビエコーを取得し補正パルスを印加する場合も本実施形態に含まれる。
【００５３】
また、Ｇｐ、Ｇｒ方向の位相補正に用いた基準ナビエコーを計測内の1番目のエコーではなく、あらかじめ用意したデータ列を用いてもよい。このようなデータ列は、例えば、図４に示すパルスシーケンスに先行して、ＲＦパルスとＧｐ、Ｇｒ方向の読出し傾斜磁場を印加してエコーを収集するプリスキャンを行い、Ｇｐ、Ｇｒ方向の各エコーを基準ナビエコーのデータ列とする。このようなプリスキャンでは、画像用の信号を取得する本スキャンに用いたパルスシーケンスと同じＲＦパルスを用い、ＲＦパルスから各エコーを取得するタイミングを同じにすることが好ましい。ＭＰＧパルスについては用いてもよいし、用いなくてもよい。ＭＰＧパルスを用いた場合には、計測内の１番目のエコーを用いた場合と同様に、撮像面内で位相の傾きを揃えることができる。またＭＰＧパルスを用いないで取得したエコーを基準にした場合には、ＭＰＧパルスに起因する位相の傾き自体をゼロにする補正を行うことができる。
【００５４】
＜第二実施形態＞
第一実施形態では、Ｇｓ、Ｇｐ、Ｇｒ方向の３方向について、リアルタイムの位相補正を行ったが、本実施形態では、リアルタイムの位相補正はＧｓ方向のみとし、Ｇｐ、Ｇｒ方向については、補正パルスを印加するのではなく、算出した補正量を用いて事後的にエコー信号を補正する。具体的には位相補正量を画像再構成におけるエコー信号のｋ空間配置座標に反映する。
【００５５】
本実施形態のＭＲＩ装置も基本的な構成は第一実施形態と同様である。ただし本実施形態では、図９に示すように、ディジタル信号処理装置８が、補正パルス算出部８１に加え、データ補正部８２を備える。以下、適宜第一実施形態で用いた図面を援用して、第一実施形態と異なる点を中心に本実施形態の処理を説明する。
【００５６】
本実施形態でも、図４に示すようなパルスシーケンス２００Ａを実行し、Ｇｓ方向について、図５の処理５０１〜５０５を行う。Ｇｐ、Ｇｒ方向については、図６に示す処理６０５〜６０７までを行い、基準ナビエコーに対する位相の傾きを算出する。データ補正部８２は、この位相の傾きを読み込み、エコー信号２０７のｋ空間配置座標の補正値を算出する。具体的には、処理６０７において式（２）で算出した１次直線の傾きFirstOrderPhaseから次式に従い補正値Δｋを算出する。
【００５７】
【数５】
2018015196.tif 000007
【００５８】
ディジタル信号処理装置８は、パルスシーケンス２００Ａにより収集したエコー信号２０７を読み込み、ｋ空間に配置する。その際、データ補正部８２で算出した補正値Δｋをエコー信号のk空間配置座標に加算する。これにより、収集前にＧｐ、Ｇｒ方向の位相補正がなされたエコー信号２０７のｋ空間配置と実質的に等しいｋ空間配置となる。補正後のｋ空間データを用いて画像再構成することにより、スライス面と平行な断面においても、ＭＰＧパルスの位相誤差の影響を抑制することができる。
【００５９】
＜第二実施形態の変形例＞
第二実施形態では、Ｇｐ、Ｇｒ方向の位相アンラップ処理後の位相傾きを求める処理（図６の処理６０７）において、Phase_unwrapped(x)に対して最小二乗法により１次直線を当てはめ、その傾きを算出したが、本変形例では、この処理６０７、及び、Ｇｓ方向の位相アンラップ後の位相傾きを求める処理５０４において、０次のオフセットも求め、その値を画像再構成時にエコー信号から差し引く補正を行う。０次のオフセットは次式（６）から算出する。
【００６０】
【数６】
2018015196.tif 000008
【００６１】
データ補正部８２は、式（６）で算出した０次のオフセットを、次式（７）に従い画像再構成に用いるエコー信号Echo(t)から位相差分する。
【数７】
2018015196.tif 000009
式（７）中、|Echo(t) |は、複素データ列であるエコー信号Echo(t)の振幅を表す。またiは虚数単位である。
【００６２】
本変形例によれば、Ｇｐ、Ｇｒ方向についても、ＭＰＧパルスによる位相誤差をゼロとする補正がなされるので、さらに高画質（ＳＮＲ／コントラスト）の拡散強調画像を得ることができる。
【００６３】
＜第三実施形態＞
第一実施形態では、ＭＰＧパルスの過不足量の算出において、Ｇｓ方向の処理では基準エコーを用いずに処理を行い、Ｇｐ、Ｇｒ方向の処理では最初の繰り返しで計測したエコーを基準とし、それとの差分をなくす補正処理を行ったが、本実施形態は、パルスシーケンス内で基準エコー群を取得する処理を加えること、及び、これら基準エコー群を用いたＭＰＧパルス過不足量の算出処理をＧｓ方向についても適用することが特徴である。
【００６４】
以下、第一実施形態と異なる部分について説明する。
図１０は本実施形態の処理の概要を示すブロックである。本実施形態のＭＲＩ装置は、図３のブロック図との比較からわかるように、シーケンサ４の処理に、基準ナビエコー読み出し傾斜磁場印加処理９０１が追加されている。またディジタル信号処理装置８における補正パルス算出処理９０２は、基準ナビエコーとナビエコーを用いた処理となる。
【００６５】
まず図１１に示すパルスシーケンス２００Ｂを参照して、シーケンサ４の処理を詳述する。本実施形態においても、基本的なパルスシーケンスとして図２に示すＤＷＩパルスシーケンスを用いる場合を説明する。図１０中、図２及び図４と同じ要素は同じ符号で示す。
【００６６】
エコー時間ＴＥを反転ＲＦパルス２０２までの前半と反転ＲＦ２０２からエコー２０７までの後半に分けると、ＴＥの後半に補正シーケンス４００を追加することにより、ＴＥの前半には補正シーケンスの時間に相当する空き時間が発生する。本実施形態では、その空き時間を利用して基準ナビエコーを収集する。具体的には、ＲＦパルス２０１と第一のＭＰＧパルス２０３との間で、Ｇｓ、Ｇｐ、Ｇｒの３方向について順番に読出し傾斜磁場印加しナビエコーを発生させる処理９０１（基準ナビエコー読み出し傾斜磁場印加処理）を実行する。残りのシーケンスは図４と同じであり、ナビエコー読出し傾斜磁場印加処理３０１を実行する。
【００６７】
ディジタル信号処理装置８（補正パルス算出部８１）は、基準ナビエコー読み出し傾斜磁場印加処理９０１で発生し、受信部６で受信したエコー信号と、パルスシーケンスの後半で実行したナビエコー読出し傾斜磁場印加処理３０１で得られたエコー信号とを用いてＭＰＧパルスの過不足量の算出（処理９０２）を行う。以下、図１２のフローを参照して、補正パルス算出部８１の処理９０２を説明する。図１２中、図６と同じ処理は同じ符号で示す。
【００６８】
まず、処理１２０１において処理９０１で取得した基準ナビエコーBase_Navi_Echo(t)を読み込む。続いて処理６０１にてナビエコーNavi_Echo(t)を読み込んだ後、処理１００２と処理６０２にて、基準ナビエコーBase_Navi_Echo(t)とナビエコーNavi_Echo(t)をそれぞれフーリエ変換し、Base_Navi (x)とNavi (x)を得る。次いで処理６０５で基準ナビエコーBase_Navi_Echo(t)とナビエコーNavi(x)の位相差分を計算する。
【００６９】
【数８】
2018015196.tif 000010
式（８）で得られるNavi_subtracted(x)は、第一実施形態の式（４）の右辺と同じであり、続く処理は第一実施形態の処理６０６〜６０８と同様である。すなわち、式（１）による位相アンラップ処理（６０６）、式（２）による位相の傾きの算出（６０７）、及び、式（３）による補正パルス印加量の算出（６０８）を行う。
【００７０】
本実施形態においても、受信コイル１４ｂが複数ある場合には、第一実施形態と同様に、各コイルのナビエコー毎に式（２）に従って位相の傾きを算出し、その平均値を用いて式（３）により補正パルス印加量を決定する。
【００７１】
本実施形態によれば、第一実施形態と同様に、全方向について大域的運動に起因する位相誤差を補正することができる。また本実施形態では、ＭＰＧパルス印加前に取得する基準ナビエコーを用いるため、高周波コイル（受信コイル）１４bの形状や位置によってナビエコーに与える位相分布も基準エコーとの位相差分の計算により取り除かれ、ＭＰＧパルス印加中の動きによってのみ生じた誤差を算出することができる。第一実施形態において用いた基準エコーは、ＭＰＧパルス印加中の動きによって生じた誤差を含んだエコーであるため、Ｇｓ方向ではＳＮＲの低下につながるため用いていない。本実施形態では、ＭＰＧパルス印加前に取得することでＭＰＧパルス印加中の動きによって生じる誤差を含まないエコーを基準とするため、Ｇｓ方向の処理にも用いることができる。従って本実施形態では、Ｇｓ方向とＧｒ、Ｇｐ方向の区別をする必要がない。
【００７２】
また、本実施形態によれば、基準ナビエコーを取得するために用いる読み出し傾斜磁場は、補正シーケンス４００の実行区間で印加するナビエコー取得のための読み出し傾斜磁場と同じ形状のため、ＴＥの延長なくパルスシーケンスを実施することができる。
【００７３】
なお本実施形態においても、第二実施形態と同様に、Ｇｒ、Ｇｐ方向については、処理６０７で算出した位相の傾きを用いて、事後的にエコー信号２０７のｋ空間配置座標を補正してもよいし、さらに第二実施形態の変形例のように、処理６０７で位相の傾きと共に０次のオフセットを算出し、それを用いてエコー信号の位相補正を行ってもよい。その他、矛盾がない限り第一実施形態で説明した変更例は本実施形態にも採用することができる。
【００７４】
以上、本発明の実施形態とその変更例を説明したが、本発明は、上述した実施形態や参考にした図面に限定されず、種々の変更を行うことが可能である。また、各実施形態の説明のために示したフローチャートの処理手順は一例であり、省略できる処理もあるし、また必要に応じて別の処理を追加することも可能である。
【符号の説明】
【００７５】
２・・・静磁場発生部、３・・・傾斜磁場発生部（撮像部）、４・・・シーケンサ（シーケンス制御部）、５・・・送信部（撮像部）、６・・・受信部（撮像部）、７・・・信号処理部、８・・・ディジタル信号処理装置（演算部）、８１・・・補正パルス算出部、８２・・・データ補正部、２００（２００Ａ、２００Ｂ）・・・ＤＷＩパルスシーケンス、４００・・・補正シーケンス。
【図１】
2018015196.tif 000011
【図２】
2018015196.tif 000012
【図３】
2018015196.tif 000013
【図４】
2018015196.tif 000014
【図５】
2018015196.tif 000015
【図６】
2018015196.tif 000016
【図７】
2018015196.tif 000017
【図８】
2018015196.tif 000018
【図９】
2018015196.tif 000019
【図１０】
2018015196.tif 000020
【図１１】
2018015196.tif 000021
【図１２】
2018015196.tif 000022