(57)【要約】      （修正有）【課題】リグリッディング処理を従来よりも高精度に実行することができる磁気共鳴イメージング装置（ＭＲＩ）を提供する。【解決手段】ＭＲＩ装置１０は、第１収集部と、第２収集部と、画像再構成部６２とを有する。第１収集部は、読み出し方向傾斜磁場の印加の下で撮像領域の少なくとも一部を含む領域から収集されたＭＲ信号をサンプリングし、第１ｋ空間データを生成する。第２収集部は、第１収集部による収集時とは収集領域が読み出し方向にずれるように、読み出し方向傾斜磁場の印加の下で撮像領域の少なくとも一部を含む領域から収集されたＭＲ信号をサンプリングし、第２ｋ空間データを生成する。画像再構成部は、第１ｋ空間データと第２ｋ空間データとの位相差を示す位相差データ、及び、本スキャンで収集されたＭＲ信号に基づいて生成されるｋ空間データから画像データを再構成する。【選択図】図１
2015000161.tif 000002
【特許請求の範囲】
【請求項１】
本スキャンの撮像領域から収集された複数の核磁気共鳴信号をサンプリングすることで複数のマトリクス要素で構成される本スキャンｋ空間データを生成し、前記本スキャンｋ空間データに基づいて前記撮像領域の画像データを再構成する磁気共鳴イメージング装置であって、
読み出し方向傾斜磁場の印加が含まれる第１パルスシーケンスの実行により、前記撮像領域の少なくとも一部を含む領域から収集された核磁気共鳴信号をサンプリングすることで、第１ｋ空間データを生成する第１収集部と、
前記読み出し方向傾斜磁場の印加が含まれる第２パルスシーケンスの実行により、前記第１収集部による核磁気共鳴信号の収集時とは収集領域が読み出し方向にずれるように、前記撮像領域の少なくとも一部を含む領域から収集された核磁気共鳴信号をサンプリングすることで、第２ｋ空間データを生成する第２収集部と、
前記第１ｋ空間データと前記第２ｋ空間データとの間の前記読み出し方向の位相差を示す位相差データを算出する位相差データ算出部と、
前記読み出し方向傾斜磁場及び位相エンコード方向傾斜磁場の印加が含まれる本スキャンパルスシーケンスを実行することで、前記撮像領域から複数の核磁気共鳴信号を収集する本スキャン実行部と、
前記本スキャン実行部により収集された複数の核磁気共鳴信号と、前記位相差データとに基づいて前記本スキャンｋ空間データを生成することで、前記画像データを再構成する画像再構成部と
を備えることを特徴とする磁気共鳴イメージング装置。
【請求項２】
請求項１記載の磁気共鳴イメージング装置において、
前記位相差データ算出部は、前記読み出し方向傾斜磁場が印加される際の経過時間の関数として前記位相差データを算出し、
前記画像再構成部は、
前記読み出し方向傾斜磁場の波形を前記位相差データに基づいて算出後、前記本スキャン実行部により収集された複数の核磁気共鳴信号と、前記読み出し方向傾斜磁場の波形とに基づいて、各々の前記マトリクス要素に対応するサンプリング期間までの前記読み出し方向傾斜磁場の強度の時間積分値が等間隔になるように、前記本スキャンｋ空間データを生成するリグリッディング処理部と、
前記本スキャンｋ空間データに対して、フーリエ変換が含まれる画像再構成処理を施すことで前記画像データを再構成するフーリエ変換部とを有する
ことを特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項２記載の磁気共鳴イメージング装置において、
前記位相差データ算出部は、前記第１ｋ空間データにおける実効エコー時間のタイミングで収集された核磁気共鳴信号に対応するデータと、前記第２ｋ空間データにおける実効エコー時間のタイミングで収集された核磁気共鳴信号に対応するデータとに基づいて、前記位相差データを算出する
ことを特徴とする磁気共鳴イメージング装置。
【請求項４】
請求項２又は請求項３記載の磁気共鳴イメージング装置において、
前記第１収集部は、前記読み出し方向傾斜磁場の極性反転の繰り返しによって収集された核磁気共鳴信号をサンプリングすることで前記第１ｋ空間データを生成し、
前記第２収集部は、前記読み出し方向傾斜磁場の極性反転の繰り返しによって収集された核磁気共鳴信号をサンプリングすることで前記第２ｋ空間データを生成し、
前記本スキャン実行部は、前記本スキャンパルスシーケンスとして、前記読み出し方向傾斜磁場の極性反転を繰り返すエコープラナーイメージングのパルスシーケンスを実行する
ことを特徴とする磁気共鳴イメージング装置。
【請求項５】
請求項４記載の磁気共鳴イメージング装置において、
前記リグリッディング処理部は、複数の核磁気共鳴信号に対する各サンプリング期間の各代表時刻を積分期間の終期とする前記時間積分値が互いに等間隔になるように、不等間隔なサンプリング期間を定めてサンプリングすることで、前記本スキャンｋ空間データを生成する
ことを特徴とする磁気共鳴イメージング装置。
【請求項６】
請求項４記載の磁気共鳴イメージング装置において、
前記リグリッディング処理部は、複数の核磁気共鳴信号を時間的に等間隔でサンプリングすることで前記本スキャンｋ空間データを生成後、前記本スキャンｋ空間データの各マトリクス要素に対応するサンプリング期間の各代表時刻を積分期間の終期とする前記時間積分値が互いに等間隔となるように、前記本スキャンｋ空間データを再配列し、
前記フーリエ変換部は、再配列後の前記本スキャンｋ空間データに対して、フーリエ変換が含まれる画像再構成処理を施すことで前記画像データを再構成する
ことを特徴とする磁気共鳴イメージング装置。
【請求項７】
請求項２乃至請求項６のいずれか１項に記載の磁気共鳴イメージング装置において、
複数の核磁気共鳴信号を検出する受信ＲＦコイルと、前記受信ＲＦコイルにより検出された複数の核磁気共鳴信号を収集して前記画像再構成部側に送信するＲＦ受信器とをさらに備え、
前記第２収集部は、前記第１収集部による核磁気共鳴信号の収集時とは前記ＲＦ受信器の検波周波数をずらして核磁気共鳴信号を収集する
ことを特徴とする磁気共鳴イメージング装置。
【請求項８】
請求項２乃至請求項７のいずれか１項に記載の磁気共鳴イメージング装置において、
前記第１収集部及び前記第２収集部は、前記本スキャン実行部による複数の核磁気共鳴信号の収集領域よりも読み出し方向に拡張された領域から核磁気共鳴信号をそれぞれ収集する
ことを特徴とする磁気共鳴イメージング装置。
【請求項９】
請求項２乃至請求項８のいずれか１項に記載の磁気共鳴イメージング装置において、
前記読み出し方向傾斜磁場の印加が含まれる第３パルスシーケンスの実行により前記撮像領域の少なくとも一部を含む領域から収集された核磁気共鳴信号をサンプリングすることで、第３ｋ空間データを生成する第３収集部をさらに備え、
前記第１収集部及び前記第２収集部は、前記第３収集部による核磁気共鳴信号の収集時とは、実効エコー時間での前記読み出し方向傾斜磁場の極性が逆になるように、前記第１パルスシーケンス及び前記第２パルスシーケンスをそれぞれ実行し、
前記画像再構成部は、前記第１ｋ空間データ及び前記第３ｋ空間データから得られる位相補正データに基づいて、前記本スキャン実行部により収集された複数の核磁気共鳴信号に含まれる位相誤差を補正しつつ、前記画像データを再構成する
ことを特徴とする磁気共鳴イメージング装置。
【請求項１０】
請求項９記載の磁気共鳴イメージング装置において、
前記第３収集部は、励起パルスの印加タイミングを基準とした前記読み出し方向傾斜磁場の印加開始タイミング、前記実効エコー時間、及び、前記実効エコー時間での前記読み出し方向傾斜磁場の極性が前記本スキャンパルスシーケンスと同じになるように、前記第３パルスシーケンスを実行し、
前記第３パルスシーケンスにおいて前記実効エコー時間での核磁気共鳴信号が奇数番目のエコーとなる場合、前記第１収集部及び前記第２収集部は、前記実効エコー時間での核磁気共鳴信号が偶数番目のエコーとなるように、前記第１パルスシーケンス及び前記第２パルスシーケンスをそれぞれ実行し、
前記第３パルスシーケンスにおいて前記実効エコー時間での核磁気共鳴信号が偶数番目のエコーとなる場合、前記第１収集部及び前記第２収集部は、前記実効エコー時間での核磁気共鳴信号が奇数番目のエコーとなるように、前記第１パルスシーケンス及び前記第２パルスシーケンスをそれぞれ実行する
ことを特徴とする磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、磁気共鳴イメージングに関する。
【背景技術】
【０００２】
ＭＲＩは、静磁場中に置かれた被検体の原子核スピンをラーモア周波数のＲＦパルスで磁気的に励起し、この励起に伴って発生するＭＲ信号から画像を再構成する撮像法である。なお、上記ＭＲＩは磁気共鳴イメージング(Magnetic Resonance Imaging)の意味であり、ＲＦパルスは高周波パルス(radio frequency pulse)の意味であり、ＭＲ信号は核磁気共鳴信号(nuclear magnetic resonance signal)の意味である。
【０００３】
ＭＲＩには、エコープラナーイメージング(EPI: Echo Planar Imaging)と呼ばれる高速撮像法がある（例えば、特許文献１及び特許文献２参照）。ＥＰＩは、１回の核磁気励起に対して読み出し方向傾斜磁場を高速で連続的に反転させ、連続的にエコー（ＭＲ信号）を生じさせてスキャンを行うものである。具体的には、ＥＰＩでは、励起パルスを印加後、ｘｙ平面内の磁化が横緩和により減衰して消滅する前に、位相エンコード量を順次変えつつ連続的なグラジエントエコーを発生させて、画像再構成に必要な全てのデータを収集する。
【０００４】
ＥＰＩには、励起パルス及び再収束パルスの後に発生するＭＲ信号を収集するスピンエコー法(SE: spin echo)を用いたスピンエコー系ＥＰＩや、励起パルスの印加後に発生するＭＲ信号を収集するフィールドエコー法（FE: field echo）を用いたフィールドエコー系ＥＰＩや、ＦＦＥ(Fast FE)法を用いたＦＦＥ系ＥＰＩがある。
【０００５】
また、複数回に亘る励起パルスを印加して得られるエコートレインのデータを合わせて１枚分の画像データを生成するＥＰＩは、マルチショットＥＰＩと呼ばれる。これに対して、１回の励起パルスの印加のみで画像を再構成するＥＰＩは、シングルショット(SS: single shot)ＥＰＩと呼ばれる。
【０００６】
ＥＰＩの読み出し方向傾斜磁場のパルス波形は、他の撮像法に比べると、パルス幅が短く、パルス周期が短いので、周波数成分が高いと言える。
【０００７】
他方、傾斜磁場パルスは、傾斜磁場コイルにパルス状の電流を供給することで生成される。傾斜磁場コイルに供給されるパルス電流の波形は、理想的には矩形波であるが、実際には、立ち上がり領域と立ち下り領域とを有する台形波となる。この結果、傾斜磁場のパルス波形も、理想的な矩形波とはならず、立ち上がり領域と立ち下り領域を有する台形波となる。
【０００８】
一般に、ＥＰＩ等の高速撮像法では、傾斜磁場パルスのパルス幅は短く、パルス両端における立ち上り領域及び立ち下り領域のパルス幅全体に対する比率は高くなる。このため、パルスの平坦な領域だけでなく、立ち上り及び立ち下り領域においてもデータをサンプリングし、画像再構成用のデータとして使用する手法が提案されている。
【０００９】
立ち上り領域及び立ち下り領域においてデータをサンプリングする手法は、ランプサンプリング(Ramp Sampling)と呼ばれる。傾斜磁場強度が平坦な領域のみをサンプリングする手法に比べると、ランプサンプリングではデータ収集期間が短縮される。
【００１０】
しかし、立ち上り領域及び立ち下り領域において時間的に等間隔でサンプリングされるＭＲ信号の生データは、傾斜磁場が変化しているときにサンプリングされるため、ｋ空間上では等間隔にならない。従って、サンプリングされたＭＲ信号の生データを画像再構成前にｋ空間上で等間隔となるように再配列することが望ましい。この再配列の処理は、リグリッディング（regridding）と呼ばれる。
【００１１】
特許文献３の従来技術では、表皮効果や渦電流などが考慮された実際の傾斜磁場発生システムに近い等価回路モデルに基づいて傾斜磁場波形を算出している。そして、このように算出された傾斜磁場波形に基づいてリグリッディング処理を実行することで、リグリッディング処理の精度の向上を図っている。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２００１−１４９３４２号公報
【特許文献２】特開２０１２−５５６８４号公報
【特許文献３】特開２０１３−１７８１１号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
特許文献３の発明は、上述のように優れた作用効果を有するものの、画質を向上させるためには、リグリッディング処理をできる限り高精度に実行することが望ましい。
【００１４】
このため、ＭＲＩにおいて、リグリッディング処理を従来よりも高精度に実行する新技術が要望されていた。
【課題を解決するための手段】
【００１５】
一実施形態のＭＲＩ装置は、本スキャンの撮像領域から収集された複数のＭＲ信号をサンプリングすることで複数のマトリクス要素で構成される本スキャンｋ空間データを生成し、本スキャンｋ空間データに基づいて撮像領域の画像データを再構成する。このＭＲＩ装置は、第１収集部と、第２収集部と、位相差データ算出部と、本スキャン実行部と、画像再構成部とを有する。
第１収集部は、読み出し方向傾斜磁場の印加が含まれる第１パルスシーケンスの実行により、撮像領域の少なくとも一部を含む領域から収集されたＭＲ信号をサンプリングすることで、第１ｋ空間データを生成する。
第２収集部は、読み出し方向傾斜磁場の印加が含まれる第２パルスシーケンスの実行により、第１収集部によるＭＲ信号の収集時とは収集領域が読み出し方向にずれるように、撮像領域の少なくとも一部を含む領域から収集されたＭＲ信号をサンプリングすることで、第２ｋ空間データを生成する。
位相差データ算出部は、第１ｋ空間データと第２ｋ空間データとの間の読み出し方向の位相差を示す位相差データを算出する。
本スキャン実行部は、読み出し方向傾斜磁場及び位相エンコード方向傾斜磁場の印加が含まれる本スキャンパルスシーケンスを実行することで、撮像領域から複数のＭＲ信号を収集する。
画像再構成部は、本スキャン実行部により収集された複数のＭＲ信号と、位相差データとに基づいて本スキャンｋ空間データを生成することで、画像データを再構成する。
【図面の簡単な説明】
【００１６】
【図１】本実施形態におけるＭＲＩ装置の全体構成を示すブロック図。
【図２】位相誤差補正用のテンプレートスキャンＡ、Ｂ、及び、本スキャンの各パルスシーケンスの一例を示すタイミング図。
【図３】位相エンコード方向及び周波数エンコード方向のマトリクス要素数が２５６×２５６の場合のｋ空間データの生成方法の一例を示す模式図。
【図４】リグリッディング処理用のテンプレートスキャンＢ、Ｃ、及び、本スキャンの各パルスシーケンスの一例を示すタイミング図。
【図５】読み出し方向傾斜磁場が非線形な領域において等時間間隔でサンプリングされたＭＲ信号が、ｋ空間上では不等間隔になることを示す概念図。
【図６】本実施形態におけるリグリッディング処理の第１方法の概念を示す模式図。
【図７】本実施形態におけるリグリッディング処理の第２方法の概念を示す模式図。
【図８】本実施形態において、リグリッディング処理の第２方法が採用される場合のＭＲＩ装置の動作の流れの一例を示すフローチャート。
【図９】本実施形態において、リグリッディング処理の第１方法が採用される場合のＭＲＩ装置の動作の流れの一例を示すフローチャート。
【発明を実施するための形態】
【００１７】
以下、ＭＲＩ装置及びＭＲＩ方法の実施形態について、添付図面に基づいて説明する。なお、各図において同一要素には同一符号を付し、重複する説明を省略する。
【００１８】
＜本実施形態の構成＞
図１は、本実施形態におけるＭＲＩ装置１０の全体構成を示すブロック図である。ここでは一例として、ＭＲＩ装置１０の構成要素を寝台ユニット２０、ガントリ３０、制御装置４０の３つに分けて説明する。
【００１９】
第１に、寝台ユニット２０は、寝台２１と、天板２２と、寝台２１内に配置される天板移動機構２３とを有する。天板２２の上面には、被検体Ｐが載置される。また、天板２２内には、被検体ＰからのＭＲ信号を検出する受信ＲＦコイル２４が配置される。さらに、天板２２の上面には、装着型のＲＦコイル装置１００が接続される接続ポート２５が複数配置される。
【００２０】
寝台２１は、天板２２を水平方向（装置座標系のＺ軸方向）に移動可能に支持する。天板移動機構２３は、天板２２がガントリ３０外に位置する場合に、寝台２１の高さを調整することで、天板２２の鉛直方向の位置を調整する。また、天板移動機構２３は、天板２２を水平方向に移動させることで天板２２をガントリ３０内に入れ、撮像後には天板２２をガントリ３０外に出す。
【００２１】
第２に、ガントリ３０は、例えば円筒状に構成され、撮像室に設置される。ガントリ３０は、静磁場磁石３１と、シムコイルユニット３２と、傾斜磁場コイルユニット３３と、ＲＦコイルユニット３４とを有する。
【００２２】
静磁場磁石３１は、例えば超伝導コイルであり、円筒状に構成される。静磁場磁石３１は、後述の制御装置４０の静磁場電源４２から供給される電流により、撮像空間に静磁場を形成する。撮像空間とは例えば、被検体Ｐが置かれて、静磁場が印加されるガントリ３０内の空間を意味する。なお、静磁場電源４２を設けずに、静磁場磁石３１を永久磁石で構成してもよい。
【００２３】
シムコイルユニット３２は、例えば円筒状に構成され、静磁場磁石３１の内側において、静磁場磁石３１と軸を同じにして配置される。シムコイルユニット３２は、後述の制御装置４０のシムコイル電源４４から供給される電流により、静磁場を均一化するオフセット磁場を形成する。
【００２４】
傾斜磁場コイルユニット３３は、例えば円筒状に構成され、シムコイルユニット３２の内側に配置される。傾斜磁場コイルユニット３３は、Ｘ軸傾斜磁場コイル３３ｘと、Ｙ軸傾斜磁場コイル３３ｙと、Ｚ軸傾斜磁場コイル３３ｚとを有する。
【００２５】
本明細書では、特に断りのない限り、Ｘ軸、Ｙ軸、Ｚ軸は装置座標系であるものとする。ここでは一例として、装置座標系のＸ軸、Ｙ軸、Ｚ軸を以下のように定義する。
まず、鉛直方向をＹ軸方向とし、天板２２は、その上面の法線方向がＹ軸方向となるように配置される。天板２２の水平移動方向をＺ軸方向とし、ガントリ３０は、その軸方向がＺ軸方向となるように配置される。Ｘ軸方向は、これらＹ軸方向、Ｚ軸方向に直交する方向であり、図１の例では天板２２の幅方向である。
【００２６】
Ｘ軸傾斜磁場コイル３３ｘは、後述のＸ軸傾斜磁場電源４６ｘから供給される電流に応じたＸ軸方向の傾斜磁場Ｇｘを撮像領域に形成する。同様に、Ｙ軸傾斜磁場コイル３３ｙは、後述のＹ軸傾斜磁場電源４６ｙから供給される電流に応じたＹ軸方向の傾斜磁場Ｇｙを撮像領域に形成する。同様に、Ｚ軸傾斜磁場コイル３３ｚは、後述のＺ軸傾斜磁場電源４６ｚから供給される電流に応じたＺ軸方向の傾斜磁場Ｇｚを撮像領域に形成する。
【００２７】
そして、スライス選択方向傾斜磁場Ｇｓｓ、位相エンコード方向傾斜磁場Ｇｐｅ、及び、読み出し方向（周波数エンコード方向）傾斜磁場Ｇｒｏは、装置座標系の３軸方向の傾斜磁場Ｇｘ、Ｇｙ、Ｇｚの合成により、任意の方向に設定可能である。
【００２８】
上記撮像領域は、例えば、１画像又は１セットの画像の生成に用いられるＭＲ信号の収集範囲の少なくとも一部であって、画像となる領域である。撮像領域は例えば、撮像空間の一部として装置座標系で３次元的に規定される。例えば折り返しアーチファクトを防止するために、画像化される領域よりも広範囲でＭＲ信号が収集される場合、撮像領域はＭＲ信号の収集範囲の一部である。一方、ＭＲ信号の収集範囲の全てが画像となり、ＭＲ信号の収集範囲と撮像領域とが合致する場合もある。また、上記「１セットの画像」は、例えばマルチスライス撮像などのように、１のパルスシーケンスで複数画像のＭＲ信号が一括的に収集される場合の複数画像である。
【００２９】
ＲＦコイルユニット３４は、例えば円筒状に構成され、傾斜磁場コイルユニット３３の内側に配置される。ＲＦコイルユニット３４は、例えば、ＲＦパルスの送信及びＭＲ信号の受信を兼用する全身用コイルや、ＲＦパルスの送信のみを行う送信ＲＦコイルを含む。
【００３０】
第３に、制御装置４０は、静磁場電源４２と、シムコイル電源４４と、傾斜磁場電源４６と、ＲＦ送信器(Radio Frequency Transmitter)４８と、ＲＦ受信器(Radio Frequency Receiver)５０と、シーケンスコントローラ５８と、演算装置６０と、入力装置７２と、表示装置７４と、記憶装置７６とを有する。
【００３１】
傾斜磁場電源４６は、Ｘ軸傾斜磁場電源４６ｘと、Ｙ軸傾斜磁場電源４６ｙと、Ｚ軸傾斜磁場電源４６ｚとを有する。Ｘ軸傾斜磁場電源４６ｘ、Ｙ軸傾斜磁場電源４６ｙ、Ｚ軸傾斜磁場電源４６ｚは、傾斜磁場Ｇｘ、Ｇｙ、Ｇｚを形成するための各電流を、Ｘ軸傾斜磁場コイル３３ｘ、Ｙ軸傾斜磁場コイル３３ｙ、Ｚ軸傾斜磁場コイル３３ｚにそれぞれ供給する。
【００３２】
ＲＦ送信器４８は、シーケンスコントローラ５８から入力される制御情報に基づいて、核磁気共鳴を起こすラーモア周波数のＲＦ電流パルスを生成し、これをＲＦコイルユニット３４に送信する。このＲＦ電流パルスに応じたＲＦパルスが、ＲＦコイルユニット３４から撮像領域に送信される。
【００３３】
ＲＦコイルユニット３４の全身用コイル、受信ＲＦコイル２４は、被検体Ｐ内の原子核スピンがＲＦパルスによって励起されることで発生したＭＲ信号を検出し、検出されたＭＲ信号は、ＲＦ受信器５０に入力される。
【００３４】
ＲＦ受信器５０は、受信したＭＲ信号に所定の信号処理を施した後、Ａ／Ｄ(analog to digital)変換を施すことで、デジタル化されたＭＲ信号の複素データである生データを生成する。ＲＦ受信器５０は、ＭＲ信号の生データを演算装置６０（の画像再構成部６２及び位相差データ算出部６５）に入力する。
【００３５】
シーケンスコントローラ５８は、演算装置６０の指令に従って、傾斜磁場電源４６、ＲＦ送信器４８及びＲＦ受信器５０の駆動に必要な制御情報を記憶する。ここでの制御情報とは、例えば、傾斜磁場電源４６に印加すべきパルス電流の強度や印加時間、印加タイミング等の動作制御情報を記述したシーケンス情報である。シーケンスコントローラ５８は、記憶した所定のシーケンスに従って傾斜磁場電源４６、ＲＦ送信器４８及びＲＦ受信器５０を駆動させることで、傾斜磁場Ｇｘ、Ｇｙ、Ｇｚ及びＲＦパルスを発生させる。
【００３６】
演算装置６０は、システム制御部６１と、システムバスＳＢと、画像再構成部６２と、画像データベース６３と、画像処理部６４と、位相差データ算出部６５とを有する。
【００３７】
システム制御部６１は、本スキャンの撮像条件の設定、撮像動作及び撮像後の画像表示において、システムバスＳＢ等の配線を介してＭＲＩ装置１０全体のシステム制御を行う。
【００３８】
上記撮像条件とは例えば、どの種類のパルスシーケンスにより、どのような条件でＲＦパルス等を送信し、どのような条件で被検体ＰからＭＲ信号を収集するかを意味する。撮像条件の例としては、撮像空間内の位置的情報としての撮像領域、フリップ角、繰り返し時間ＴＲ(Repetition Time)、スライス数、撮像部位、スピンエコー法やパラレルイメージング等のパルスシーケンスの種類などが挙げられる。
【００３９】
上記撮像部位とは、例えば、頭部、胸部、腹部などの被検体Ｐのどの部分を撮像領域として画像化するかを意味する。
上記「本スキャン」は、Ｔ１強調画像などの、目的とする診断画像の撮像のためのスキャンであって、位置決め画像用のＭＲ信号収集のスキャンや、較正スキャンを含まないものとする。
スキャンとは、ＭＲ信号の収集動作を指し、画像再構成を含まないものとする。
【００４０】
較正スキャンとは例えば、本スキャンの撮像条件の内の未確定のものや、画像再構成処理や画像再構成後の補正処理に用いられる条件やデータを決定するために、本スキャンとは別に行われるスキャンを指す。較正スキャンの例としては、ＥＰＩの位相補正データを得る「テンプレートスキャン」などが挙げられる（特許文献２参照）。なお、特許文献２では、ＥＰＩにおける較正スキャンを「テンプレートショット」と記載している。
【００４１】
また、システム制御部６１は、撮像条件の設定画面情報を表示装置７４に表示させ、入力装置７２からの指示情報に基づいて撮像条件を設定し、設定した撮像条件をシーケンスコントローラ５８に入力する。また、システム制御部６１は、撮像後には、生成された表示用画像データが示す画像を表示装置７４に表示させる。
入力装置７２は、撮像条件や画像処理条件を設定する機能をユーザに提供する。
【００４２】
画像再構成部６２は、リグリッディング処理部６２ａと、フーリエ変換部６２ｂとを有する。リグリッディング処理部６２ａは、位相エンコードステップ数及び周波数エンコードステップ数に応じて、ＲＦ受信器５０から入力されるＭＲ信号の生データをｋ空間データとして配置及び保存する。ｋ空間とは、周波数空間の意味である。
【００４３】
リグリッディング処理部６２ａは、リグリッディング処理が施されたｋ空間データを生成する。リグリッディング処理の詳細については、後述する。
【００４４】
フーリエ変換部６２ｂは、後述の位相誤差の補正やフーリエ変換が含まれる画像再構成処理をｋ空間データに施すことで、被検体Ｐの画像データを生成する。画像再構成部６２は、生成した画像データを画像データベース６３に保存する。
【００４５】
なお、ＭＲＩの画像データは、例えば各画素が画素値を有することで構成される。画素値は、例えば、その画素が表示される際の輝度レベル（その画素に対応する被検体領域から検出されたＭＲ信号の強度）を示す。スライスの場合、ＭＲＩの画像データは、縦横の画素数が例えば位相エンコードステップ数×周波数エンコードステップ数となる。
【００４６】
画像処理部６４は、画像データベース６３から画像データを取り込み、これに所定の画像処理を施し、画像処理後の画像データを表示用画像データとして記憶装置７６に保存する。
記憶装置７６は、上記の表示用画像データに対し、その表示用画像データの生成に用いた撮像条件や被検体Ｐの情報（患者情報）等を付帯情報として付属させて記憶する。
【００４７】
なお、演算装置６０、入力装置７２、表示装置７４、記憶装置７６の４つを１つのコンピュータとして構成し、例えば制御室に設置してもよい。
また、上記説明では、ＭＲＩ装置１０の構成要素をガントリ３０、寝台ユニット２０、制御装置４０の３つに分類したが、これは一解釈例にすぎない。例えば、天板移動機構２３は、制御装置４０の一部として捉えてもよい。
【００４８】
或いは、ＲＦ受信器５０は、ガントリ３０外ではなく、ガントリ３０内に配置されてもよい。この場合、例えばＲＦ受信器５０に相当する電子回路基盤がガントリ３０内に配設される。そして、ＲＦコイル装置１００や受信ＲＦコイル２４等によって電磁波からアナログの電気信号に変換されたＭＲ信号は、当該電子回路基盤内のプリアンプで増幅され、デジタル信号としてガントリ３０外に出力され、画像再構成部６２及び位相差データ算出部６５に入力される。ガントリ３０外への出力に際しては、例えば光通信ケーブルを用いて光デジタル信号として送信すれば、外部ノイズの影響が軽減されるので望ましい。
【００４９】
＜本実施形態の原理説明＞
本実施形態では一例として、各スライス又は各スラブに対して較正スキャンとしての３つのテンプレートスキャンＡ、Ｂ、Ｃが実行される。テンプレートスキャンＡは、位相補正データの取得のみに用いられ、テンプレートスキャンＣは、リグリッディング処理のみに用いられ、テンプレートスキャンＢは、位相補正データの取得及びリグリッディング処理に用いられる。
【００５０】
図２は、位相誤差補正に用いられるテンプレートスキャンＡ、Ｂ、及び、本スキャンの各パルスシーケンスの一例を示すタイミング図である。図２の上段、中段、下段はそれぞれ、本スキャン(MAIN SCAN)、テンプレートスキャンＡ(TEMPLATE SCAN A)、テンプレートスキャンＢ(TEMPLATE SCAN B)の各パルスシーケンスを示し、各横軸は経過時間ｔを示す。
【００５１】
図２の上段、中段、下段において、ＲＦはＲＦパルス、Ｇｓｓはスライス選択方向傾斜磁場、Ｇｐｅは位相エンコード方向傾斜磁場、Ｇｒｏは読み出し方向傾斜磁場、MR SIGNALはエコー（ＭＲ信号）をそれぞれ示す。
【００５２】
図２の上段の本スキャンでは一例として、スピンエコー系のシングルショットＥＰＩを示す。この例では、フリップ角９０°の励起パルスの印加後、フリップ角１８０°の再収束パルスが印加される。その後、読み出し方向傾斜磁場パルスの極性の反転が繰り返されることで、ＭＲ信号が収集される。このとき、位相エンコード方向傾斜磁場Ｇｐｅとして、１のプレパルスの印加後、プレパルスとは極性が反対のブリップパルス(BLIP PULSE)が読み出し方向傾斜磁場パルスの極性の反転時に繰り返し印加されることで、位相エンコード量が順次加算される。
【００５３】
位相エンコード傾斜磁場Ｇｐｅのプレパルスの面積（強度の絶対値の時間積分値）は、始めの４つのブリップパルスの合算面積に等しい。即ち、４つ目のブリップパルスの印加後のタイミングにおいて、位相エンコード量がゼロとなる。この例では、本スキャンで５番目に収集されるＭＲ信号は、位相エンコードステップがゼロとなる上に、励起パルス及び再収束パルスの印加タイミングで決まる実効エコー時間(EFFECTIVE ECHO TIME)に該当するので、このＭＲ信号の強度が最大になる。
【００５４】
テンプレートスキャンＡは、ここでは一例として、以下の２点を除いて本スキャンと同じパルスシーケンスである。
【００５５】
第１の相違点として、テンプレートスキャンＡでは、以下の理由で位相エンコード傾斜磁場Ｇｐｅが印加されない。位相誤差の補正において、ｋ空間の中心ライン近辺の複数ラインが使用される場合を考える。この場合、位相エンコード方向傾斜磁場Ｇｐｅが印加されると、位相エンコード方向傾斜磁場Ｇｐｅの位相ずれも各ＭＲ信号に含まれ、それがｋ空間において位相エンコード軸方向の位置毎に配列されるＭＲ信号毎に異なるから、読み出し方向の位相ずれを正確に算出できない。
【００５６】
但し、本実施形態では一例として、ｋ空間の中心ラインのデータのみが位相誤差の補正に用いられるので、各ラインのＭＲ信号同士の位相エンコード方向傾斜磁場Ｇｐｅの位相ずれの違いは、あまり影響しない。従って、ｋ空間の１ラインのＭＲ信号のみが位相誤差の補正に用いられる場合、位相エンコード方向傾斜磁場Ｇｐｅを印加してもよいが、本実施形態のように位相エンコード傾斜磁場Ｇｐｅが印加されない方が望ましい。
【００５７】
第２の相違点として、テンプレートスキャンＡのＭＲ信号の収集領域（ＦＯＶ:Field Of View: 撮像視野)は、本スキャンのＦＯＶを例えば中央に含むように、読み出し方向にのみ本スキャンよりも拡大される。拡大の程度については、後述する。ここで、位相エンコード傾斜磁場Ｇｐｅが印加されないテンプレートスキャンＡで得られたｋ空間データから画像を再構成することはできない。従って、厳密には、テンプレートスキャンＡについて、ＭＲ信号の収集領域を「ＦＯＶ」と称することは正確ではない。
【００５８】
相違点は上記２点のみであるから、テンプレートスキャンＡは、位相エンコード方向及びスライス選択方向におけるＭＲ信号の収集領域が本スキャンと同じである（テンプレートスキャンＡでは位相エンコード方向傾斜磁場Ｇｐｅが印加されないので、テンプレートスキャンＡの位相エンコード方向は、厳密には読み出し方向及びスライス選択方向に直交する方向である）。
【００５９】
また、テンプレートスキャンＡは、奇数番目である５番目に収集されるＭＲ信号ＡＣ（図２の中段参照）が最大強度であり、ｋ空間データの中心ラインに配置される。この点は、励起パルス又は再収束パルスの印加タイミングを基準とした実効エコー時間も含めて、本スキャンと同じである。
【００６０】
テンプレートスキャンＢは、以下の２点を除いてテンプレートスキャンＡと同じパルスシーケンスである。
【００６１】
第１の相違点として、テンプレートスキャンＢでは、偶数番目である４番目に収集されるＭＲ信号ＢＣが最大強度となり、ｋ空間データの中心ラインに配置される（図２の下段参照）。
即ち、テンプレートスキャンＢでは、読み出し方向傾斜磁場Ｇｒｏの印加開始タイミングをテンプレートスキャンＡよりも「Ｇｒｏ反転間隔」の時間幅で遅らす。
これにより、テンプレートスキャンＢにおける９０°励起パルスの印加時を基準としたＭＲ信号の発生開始タイミングは、テンプレートスキャンＡよりも「Ｇｒｏ反転間隔」の時間幅で遅れる。
【００６２】
上記「Ｇｒｏ反転間隔」とは、読み出し方向傾斜磁場Ｇｒｏの極性が反転してから再反転するまでの時間幅である。また、ここでの「印加開始タイミング」とは、例えば、励起パルスの印加開始時刻を基準とするものである。ここでは一例として、本スキャン及び各テンプレートスキャンＡ、Ｂ、Ｃについて、励起パルス、再収束パルスの印加開始タイミングは共通である。
【００６３】
第２の相違点として、テンプレートスキャンＢでは、実効エコー時間のタイミングにおける読み出し方向傾斜磁場Ｇｒｏの極性がテンプレートスキャンＡとは逆になるように、読み出し方向傾斜磁場Ｇｒｏの印加開始時の極性が決定される。この条件の場合には静磁場不均一性による位相誤差成分を消去し易く、位相補正データの取得する上で望ましいからである。
【００６４】
より詳細には、ＥＰＩにおける位相誤差の主な原因は、静磁場不均一性と、傾斜磁場のスイッチングにより発生する渦磁場の２つと考えられている。実効エコー時間のタイミングでの読み出し方向傾斜磁場Ｇｒｏの極性を反転させた２つのテンプレートスキャンＡ、Ｂのデータの差分により、静磁場不均一性による位相誤差成分が消去されるので、それ以外の理由による位相誤差成分を抽出できる（特許文献２参照）。
【００６５】
上記条件を満たすため、図２の縦の破線で示す実効エコー時間のタイミングでの読み出し方向傾斜磁場Ｇｒｏの極性は、テンプレートスキャンＡではマイナスであり（図２の中段の５番目のＭＲ信号ＡＣのタイミング）、テンプレートスキャンＢではプラスである（図２の下段の４番目のＭＲ信号ＡＣの発生タイミング）。
【００６６】
次に、テンプレートスキャンＡ、ＢにおけるＭＲ信号の収集領域の拡大の程度について説明する。例えば、本スキャンのＦＯＶにおける位相エンコード方向及び読み出し方向の長さが互いに等しい場合、テンプレートスキャンＡ、Ｂでは、ＭＲ信号の収集領域が読み出し方向において本スキャンの２倍に拡大される。
【００６７】
読み出し方向のＦＯＶの拡大に際しては、後述するテンプレートスキャンＣのＭＲ信号の収集領域を読み出し方向にシフトさせることによって折り返し(wraparound)が生じない程度、の拡大であることが望ましい。従って、拡大の程度の目安としては、例えば、本スキャンの２倍が望ましい。
【００６８】
なお、ＦＯＶ内における被検体の領域がおよそ中央にあり、ＦＯＶに占める被検体の領域が少なければ、読み出し方向にＦＯＶを拡大する割合は、２倍よりも小さくてよい。
また、読み出し方向のＦＯＶの拡大に際しては、読み出し方向にＦＯＶを拡大しても、本スキャンと同じ分解能が維持されることが望ましい。理由は以下である。
【００６９】
テンプレートスキャンＢは、前述のように、テンプレートスキャンＣと共にリグリッディング処理にも用いられる。テンプレートスキャンＢ、Ｃの目的は、本スキャンにおける読み出し方向傾斜磁場Ｇｒｏの波形のなまり方を正確に算出することである。そうすると、読み出し方向傾斜磁場Ｇｒｏの反転の時間間隔や波高（最大強度）などのテンプレートスキャンＢ、Ｃの各条件が本スキャン実行時と同じである方が、本スキャンの読み出し方向傾斜磁場Ｇｒｏの波形をテンプレートスキャンＢ、Ｃの実行結果から正確に算出し易い。
【００７０】
また、テンプレートスキャンＢは、テンプレートスキャンＡと共に位相誤差の補正に用いられるので、テンプレートスキャンＡ、ＢはＭＲ信号の収集領域を互いに揃えることが望まれる。
【００７１】
以上の理由から、テンプレートスキャンＡ、ＢにおいてＭＲ信号の収集領域が本スキャンよりも読み出し方向に拡大される場合、テンプレートスキャンＡ、Ｂでは、周波数エンコードステップ数を本スキャンの２倍にすることが望ましい。
【００７２】
本実施形態では一例として、本スキャンでは位相エンコードステップ数及び周波数エンコードステップ数が共に２５６であり、テンプレートスキャンＡ、Ｂでは周波数エンコードステップ数が５１２であるものとする。従って、テンプレートスキャンＡ、Ｂでは、読み出し方向傾斜磁場Ｇｒｏの反転によるＭＲ信号の収集が、どちらも本スキャンの約２倍となる。
【００７３】
次に、位相誤差の補正用の位相補正データの算出方法の説明の前に、その算出の概念に関わるｋ空間データの生成方法の一例について説明する。
【００７４】
図３は、位相エンコード方向及び周波数エンコード方向のマトリクス要素数が２５６×２５６の場合のｋ空間データの生成方法の一例を示す模式図である。図３の上段は、ｋ空間データの実数部分又は虚数部分を示し、図３の下段は、読み出し方向傾斜磁場パルスの一例を示す。
【００７５】
図３の上段において、ＴＲは繰り返し時間(Repetition Time)であり、縦軸は、ｋ空間における位相エンコード軸(PHASE ENCODE AXIS)である。縦軸の右側に付した各数字は、本スキャンのように位相エンコード方向傾斜磁場Ｇｐｅが印加される場合の位相エンコードステップ番号（Phase Encode Step）である。図３の上段における横軸は、ｋ空間における周波数エンコード軸(FREQUENCY ENCODE AXIS)である。
【００７６】
また、図３において、横方向のＴｓはサンプリング時間(Sampling Time)である。図３の下段において、横軸は経過時間ｔ、即ち、サンプリング時刻ｔであり、縦軸は読み出し方向傾斜磁場パルスの磁場強度を示す。
【００７７】
通常のスピンエコー法のようなパルスシーケンスでは、位相エンコードを例えば２５６回変えて収集した２５６ラインのＭＲ信号からそれぞれ、搬送周波数の余弦関数を差し引く。この処理後の２５６のＭＲ信号を、図３の上段のように、位相エンコードステップ順に下から−１２７，−１２６，・・・−１，０，１，・・・１２７，１２８のように並べる。
【００７８】
さらに、図３の上段の周波数エンコード軸の方向（サンプリング時刻ｔの方向）に、各ＭＲ信号のライン（この例では２５６ライン）を等間隔に２５６分割したΔＴｓ毎に、ＭＲ信号の強度をマトリクス値にする。これにより、２５６×２５６のマトリクス要素からなるマトリクスデータ、即ち、ｋ空間データの実数部分が得られる。
【００７９】
また、搬送周波数の余弦関数の代わりに搬送周波数の正弦関数を引く点を除いて、上記同様に処理することで、２５６×２５６のマトリクス要素からなるｋ空間データの虚数部分が得られる。
【００８０】
一方、図２の上段のようなシングルショットのＥＰＩでは、位相エンコードステップがゼロのタイミングの前に４ライン分のＭＲ信号しか収集できないので、収集数は（２５６／２）＋４＝１３２ラインとなる。この場合、収集されなかった１２４ラインは、ｋ空間上では例えばデータとしてゼロが入る。
【００８１】
なお、テンプレートスキャンＡ、Ｂでは位相エンコード方向傾斜磁場Ｇｐｅが印加されないので、ｋ空間に埋められる全２５６のＭＲ信号のライン同士で、位相エンコード量は変わらない。
【００８２】
次に、位相誤差補正用の位相補正データの算出方法について説明する。以下の説明では、テンプレートスキャンＡで収集されたＭＲ信号から得られるｋ空間データをテンプレートデータＡと表記する。同様に、テンプレートスキャンＢ、Ｃでそれぞれ収集されたＭＲ信号から得られる各ｋ空間データをテンプレートデータＢ、テンプレートデータＣと表記する。
【００８３】
まず、テンプレートデータＡの実数部分、即ち、テンプレートスキャンＡの各ＭＲ信号から搬送周波数の正弦波を差し引くことで得られるｋ空間データの実数部分を考える。ここでは一例として、ｋ空間データの中心ラインのみを用いる。テンプレートスキャンＡでは、周波数エンコードステップ数が５１２である。
【００８４】
従って、テンプレートデータＡの実数部分の中心ライン（図２の実効エコー時間のＭＲ信号ＡＣに対応）を、ｋ空間の周波数エンコード軸方向に５１２分割し、ＭＲ信号の強度をサンプリングする。
【００８５】
これは、読み出し方向傾斜磁場Ｇｒｏの印加開始時刻ｔ０（ＥＰＩの場合は反転時刻）から、印加終了時刻ｔ５１１（ＥＰＩの場合は再反転時刻）までの時間幅を５１２分割することと捉えてもよい。なお、これと同様の概念を示す図３の下段では、周波数エンコードステップ数が２５６なので、サンプリング時刻ｔ０〜ｔ２５５に２５６分割されている。
【００８６】
このようにして、各サンプリング時刻ｔ０〜ｔ５１１における、実効エコー時間のＭＲ信号の各強度を、テンプレートデータＡの中心ラインの実数データKAreal(kr)と定義する。ここで、krは、読み出し方向の空間周波数[radian / meter]であり、(kr)は、空間周波数krの関数という意味である（以下の説明で出てくる他のパラメータも同様である）。
【００８７】
従って、KAreal(kr)は、正確には読み出し方向の空間周波数krの関数であるが、上記のようにサンプリング時刻ｔの関数として捉えてもよい。サンプリング時刻がｋ空間の中心に近いｔ２５５付近では、ＭＲ信号の強度が高く、情報量が大きく、読み出し方向の空間周波数ｋｒが低周波領域のデータとなる。一方、サンプリング時刻がｋ空間の端側であるｔ０付近やｔ５１１付近では、ＭＲ信号の強度が低く、高周波領域のデータとなる。
【００８８】
テンプレートデータＡの虚数部分の中心ラインのＭＲ信号を、上記同様にｋ空間の周波数エンコード軸方向に５１２分割してサンプリングすることで、テンプレートデータＡの中心ラインの虚数データKAimag(kr)が得られる。従って、テンプレートデータＡの中心ラインの複素データKA(kr)は、jを虚数単位とすれば次式で与えられる。
【００８９】
KA(kr) ＝ KAreal(kr) ＋ ｛j × KAimag(kr)｝  …（１）
【００９０】
上記同様にして、テンプレートデータＢの中心ラインの実数データKBreal(kr)と、テンプレートデータＢの中心ラインの虚数データKBimag(kr)とを算出後、次式によりテンプレートデータＢの中心ラインの複素データKB(kr)が得られる。
【００９１】
KB(kr) ＝ KBreal(kr) ＋ ｛j × KBimag(kr)｝  …（２）
【００９２】
次に、テンプレートデータＡの中心ラインの複素データKA(kr)をｋ空間の周波数エンコード軸方向に１次元フーリエ逆変換することで、その実空間データRA(xr)が得られる。
同様に、テンプレートデータＢの中心ラインの複素データKB(kr)をｋ空間の周波数エンコード軸方向に１次元フーリエ逆変換することで、その実空間データRB(xr)が得られる。
【００９３】
実空間データRA(xr)，RB(xr)における(xr)は、読み出し方向の位置xr[meter]の関数という意味である（以下の説明で出てくる他のパラメータも同様である）。実空間データRA(xr)，RB(xr)は、どちらも複素データである。
ここで、実空間データRA(xr)，RB(xr)の生成元のテンプレートデータＡ、Ｂは、実効エコー時間が等しく、且つ、ＭＲ信号の収集時の読み出し方向傾斜磁場Ｇｒｏの極性が逆である。従って、実空間データRA(xr)，RB(xr)の位相成分にそれぞれ含まれる静磁場不均一性による位相誤差成分は、互いに等しくなる。
【００９４】
次に、｛QQ｝＊をQQの共役複素数と定義する。そうすると、位相誤差補正用の位相補正データdV(xr)は、以下の（３）式で与えられる。
【００９５】
【数１】
2015000161.tif 000003
【００９６】
上記（３）式は、２つのベクトルRA(xr)，RB(xr)の位相差にあたる位相を持ち、絶対値が１の単位ベクトルを求めることに相当する。従って、（３）式で与えられる位相補正データdV(xr)により、静磁場不均一性による位相誤差成分は消去され、その他の位相誤差成分のみを抽出できる。
【００９７】
次に、例えば図３で示す手法により得られる本スキャンのｋ空間データをKM(shot, echo, kr）と定義する。shotは、そのｋ空間データを収集したショット番号（スキャンの順番を示す番号）である。echoは、励起パルスの印加後、何番目に収集されたＭＲ信号のデータかを意味する。例えば図２の上段の実効エコー時間で収集されるＭＲ信号の場合、echo＝５となる。
【００９８】
次に、本スキャンのｋ空間データKM(shot, echo, kr）に対して、ｋ空間の周波数エンコード軸方向に１次元フーリエ逆変換を施すことで得られる実空間データをRM(shot, echo, xr）と定義する。
【００９９】
本スキャンの実空間データRM(shot, echo, xr）に対して、位相誤差補正を施した実空間データをRM'(shot, echo, xr）と定義する。
【０１００】
ここでは一例として、本スキャンの実空間データRM(shot, echo, xr）において、echoが奇数のデータについては位相誤差の補正が行われず、echoが偶数のデータについてのみ位相誤差の補正が行われる。echoが偶数の場合、位相補正データdV(xr)の共役複素数データを用いた次式により、位相誤差が補正される。
【０１０１】
RM'(shot, echo, xr） ＝ RM(shot, echo, xr）× ｛dV(xr)｝*  …（４）
【０１０２】
即ち、本スキャンの実空間データRM(shot, echo, xr）において、偶数番目のＭＲ信号（エコー）に対してのみ、求められた誤差成分を（４）式により戻すことで、各ＭＲ信号の位相の不連続が減少する。静磁場不均一性以外の原因による位相誤差は、主に実空間の読み出し方向に１次の傾斜を持ち、ＭＲ信号の収集中の読み出し方向傾斜磁場Ｇｒｏの極性により、その位相傾斜の方向が反転する。上記補正方法は、この規則性を前提とする（例えば特開平９−２７６２４３号公報等を参照）。
【０１０３】
位相補正データdV(xr)の算出方法は新規であるが、位相補正データdV(xr)の算出後の位相誤差の補正については従来技術と同様でよいため、さらなる説明を省略する。従って、上記の位相誤差の補正方法は、一例にすぎず、他の手法でもよい。
【０１０４】
図４は、リグリッディング処理用のテンプレートスキャンＢ、Ｃ、及び、本スキャンの各パルスシーケンスの一例を示すタイミング図である。図４の上段、中段、下段はそれぞれ、本スキャン(MAIN SCAN)、テンプレートスキャンＢ(TEMPLATE SCAN B)、テンプレートスキャンＣ(TEMPLATE SCAN C)の各パルスシーケンスを示す。図４の表記方法は、図２と同様であり、図４の本スキャン及びテンプレートスキャンＢのパルスシーケンスは図２と同じである。
【０１０５】
テンプレートスキャンＣは、以下の１点を除いてテンプレートスキャンＢと同じである。テンプレートデータＣのＭＲ信号の収集領域は、テンプレートデータＢのＭＲ信号の収集領域とは、読み出し方向に所定間隔ずれる（シフトする）。
読み出し方向におけるＭＲ信号の収集領域のシフト量は、以下の２つの理由により、例えば５画素程度から１０画素程度が望ましい。
【０１０６】
第１に、１画素程度のシフト量だと、テンプレートデータＢ、Ｃの間での読み出し方向の位相差が十分に大きく表れないために、両者の位相差データΔθ（ｔ）を十分正確に算出できないおそれがある。
第２に、例えば画像の半分程度など、読み出し方向のシフト量が大きすぎると、位相データの折り返しが多く発生し、両者の位相差データΔθ（ｔ）を十分正確に算出できない。
【０１０７】
ＭＲ信号の収集領域を読み出し方向にシフトさせるためには、読み出し方向へのシフト量に比例した周波数だけ、ＭＲ信号の位相検波時に用いる搬送波の周波数をシフトさせればよい。例えば、ＭＲ信号を収集中の読み出し方向傾斜磁場パルスの強度をIGro［Tesla／meter］、ＭＲ信号の収集領域の読み出し方向でのシフト量をr0［meter］、水素原子の磁気回転比をγ［radian／(Tesla・second）］とすると、搬送波の周波数シフト量Δｆ［Ｈｚ］は次式で表される（例えば、特許文献１の段落［０００８］等を参照）。
【０１０８】
Δｆ＝γ／｛2π × IGro × r0｝  …（５）
【０１０９】
即ち、テンプレートスキャンＢでＭＲ信号を収集中のＲＦ受信器５０における検波周波数（搬送波の周波数）が例えば磁場中心のラーモア周波数である場合、テンプレートスキャンＣでの検波周波数は、磁場中心のラーモア周波数からΔｆだけシフトすればよい。
【０１１０】
本実施形態では、読み出し方向傾斜磁場パルスの印加開始時刻からの経過時間ｔの関数として、リグリッディング処理用の読み出し方向の位相差データΔθ（ｔ）を算出し、位相差データΔθ（ｔ）から読み出し方向傾斜磁場パルスの波形を正確に算出する。本実施形態では、このように正確に算出された読み出し方向傾斜磁場パルスの波形に基づいてリグリッディング処理を実行することで、リグリッディング処理の精度を向上する。
【０１１１】
なお、（ｔ）は、サンプリング時刻ｔの関数、という意味である（以下の説明の他のパラメータも同様）。また、リグリッディング処理用の位相差データΔθ（ｔ）は、ｋ空間データを用いるので、厳密には読み出し方向の空間周波数krの関数である。しかし、ＭＲ信号の収集期間内で空間周波数を変化させながらＭＲ信号を収集しているので、位相差データΔθ（ｔ）は、後述の（１０）式からサンプリング時刻ｔの関数の関数とも言える。ここでは、磁場波形の観点から、位相差データΔθ（ｔ）をサンプリング時刻ｔの関数として扱う。
【０１１２】
以下、読み出し方向の位相差データΔθ（ｔ）の算出方法、読み出し方向傾斜磁場パルスの波形の算出方法、リグリッディング処理の方法、の順に説明する。
【０１１３】
まず、ｋ空間データの実数部分における、位相エンコードステップＰＳのラインの時刻ｔでのマトリクス値を、ＫＲ（ＰＳ，ｔ）と定義する。
同様に、ｋ空間データの虚数部分における、位相エンコードステップＰＳのラインの時刻ｔでのマトリクス値を、ＫＩ（ＰＳ，ｔ）と定義する。
【０１１４】
上記時刻ｔは、図３の下段におけるサンプリング時刻ｔである。テンプレートスキャンＢ、Ｃでは、本スキャンよりも読み出し方向にＭＲ信号の収集領域が例えば２倍に拡張されるので、前述の分解能の観点から周波数エンコードステップ数が本スキャンの２倍にされる。
【０１１５】
サンプリング時刻ｔの早い順に時刻ｔ０，ｔ１，ｔ２，…，ｔ５１１の合計５１２のサンプリング時刻がそれぞれ、周波数エンコードステップの−２５５〜２５６に対応する。時刻ｔ０は、読み出し方向傾斜磁場パルスの印加開始時刻であり、時刻ｔ５１１は、当該読み出し方向傾斜磁場パルスの印加終了時刻である。なお、図３は本スキャンに対応するので、図３の下段において、周波数エンコードステップ数は５１２ではなく、２５６で記載されている。
【０１１６】
従って、例えばｋ空間データの実数部分の中心ラインにおける、ＭＲ信号の検出タイミングが最先のマトリクス値は、ＫＲ（０，ｔ０）と表される。なお、ｋ空間データの中心ラインは、実効エコー時間のタイミングで収集されたＭＲ信号に対応し、本スキャンの場合には位相エンコードステップがゼロのラインでもある。
【０１１７】
同様に、ｋ空間データの実数部分の中心ラインにおける、ＭＲ信号の検出タイミングが最後のマトリクス値は、ＫＲ（０、ｔ５１１）と表される。
同様に、ｋ空間データの虚数部分の中心ラインにおける、ＭＲ信号の検出タイミングが最先のマトリクス値は、ＫＩ（０，ｔ０）と表される。
【０１１８】
ｋ空間データでは、中心ラインにおいて、検出時のＭＲ信号の強度が最も高く、ＳＮ比(Signal to Noise ratio:信号／ノイズ比)が最も高い。そこで、ここでは一例として、ｋ空間データの実数部及び虚数部における中心ラインを用いて読み出し方向の位相データθ（ｔ）が算出される。
【０１１９】
例えば時刻ｔ０、時刻ｔ１における各位相データθ（ｔ０），θ（ｔ１）は、逆正接arctanを用いた以下の（６）式、（７）式で算出される。
【０１２０】
θ（ｔ０）＝arctan｛ＫＩ（０，ｔ０）／ＫＲ（０，ｔ０）｝  …（６）
θ（ｔ１）＝arctan｛ＫＩ（０，ｔ１）／ＫＲ（０，ｔ１）｝  …（７）
【０１２１】
上記同様にして、サンプリング時刻ｔ０〜ｔ５１１まで、５１２の位相データθ（ｔ）をサンプリング時刻ｔの関数として算出できる。
そこで、テンプレートデータＢの中心ラインから上記のように得られる読み出し方向の位相データをθｂ（ｔ）と定義し、テンプレートデータＣの中心ラインから上記のように得られる読み出し方向の位相データをθｃ（ｔ）とする。
【０１２２】
ここで、読み出し方向の位相データθｂ（ｔ）、θｃ（ｔ）は、共に中心ラインであり、共に偶数番目（図４では４番目）のエコーであるから、対応するＭＲ信号の検出期間における読み出し方向傾斜磁場Ｇｒｏの極性が同じである。従って、両者の単純な差分によって、磁場不均一性による位相への影響を相殺できるので、次式で読み出し方向の位相差データΔθ（ｔ）を算出できる。
【０１２３】
Δθ（ｔ）＝θｂ（ｔ）−θｃ（ｔ）  …（８）
【０１２４】
次に、読み出し方向の空間周波数kr（ｔ）は次式で表される。
【０１２５】
kr（ｔ）＝Δθ（ｔ）／r0  …（９）
【０１２６】
また、読み出し方向傾斜磁場パルスの印加開始時刻をサンプリング時刻ｔの起点としつつ、読み出し方向傾斜磁場パルスの磁場強度[Tesla / meter]を、サンプリング時刻ｔの関数としてＧｒｏ（ｔ）と定義する。このとき、任意のサンプリング時刻ｔにおける読み出し方向の空間周波数kr（ｔ）は、次式で表される（特許文献１の段落［００１６］等を参照）。
【０１２７】
【数２】
2015000161.tif 000004
【０１２８】
上記（１０）式において、γは前記水素原子の磁気回転比である。（１０）式をサンプリング時刻ｔで時間微分すれば、次式が成り立つ。
【０１２９】
Ｇｒｏ（ｔ）＝｛ｄ kr（ｔ）/ ｄｔ｝／γ  …（１１）
【０１３０】
（１１）式に（９）式を代入すれば、次式が成り立つ。
【０１３１】
Ｇｒｏ（ｔ）＝｛ｄ Δθ（ｔ）/ ｄｔ｝／｛γ × r0｝  …（１２）
【０１３２】
テンプレートスキャンＢ、Ｃから位相差データΔθ（ｔ）を算出すれば、（１２）式により、読み出し方向傾斜磁場パルスＧｒｏ（ｔ）の波形を正確に算出できる。
【０１３３】
なお、上記の例では、信号強度が最大となるｋ空間データの中心ラインのデータを用いるが、本発明の実施形態は、かかる態様に限定されるものではない。
例えば、中心ライン近辺の複数ライン（信号強度が高い低周波領域の複数ライン）のデータを用いて読み出し方向の位相差データΔθ（ｔ）をそれぞれ算出後、それらを平均することで最終的な位相差データΔθ（ｔ）を決定してもよい。この変形例の場合、平均により、位相差データΔθ（ｔ）の精度をさらに向上できる。
【０１３４】
また、テンプレートスキャンＢ、Ｃでは位相エンコード方向傾斜磁場Ｇｐｅが印加されない例を述べたが、本発明の実施形態は、かかる態様に限定されるものではない。前述のように、ｋ空間データの中心ラインのみからリグリッディング処理用の読み出し方向の位相差データΔθ（ｔ）を算出する場合、テンプレートスキャンＢ、Ｃにおいて位相エンコード方向傾斜磁場Ｇｐｅを印加してもよい。
【０１３５】
但し、ｋ空間データの複数ラインが位相差データΔθ（ｔ）の算出に用いられる場合、位相エンコード方向傾斜磁場Ｇｐｅを印加しないことが望ましい。位相エンコード方向傾斜磁場Ｇｐｅが印加されると、位相エンコード方向の位相ずれもＭＲ信号に含まれ、位相エンコード方向の位相ずれの程度がエコー毎に異なるから、読み出し方向の位相ずれを正確に算出しにくくなる。
【０１３６】
次に、リグリッディング処理の方法について、２つの例を説明する。
読み出し方向傾斜磁場Ｇｒｏの印加の下で等間隔又は不等間隔にＭＲ信号をサンプリングすることで得られたｋ空間データの各マトリクス値は、ｋ空間上では、読み出し方向傾斜磁場Ｇｒｏの強度（の絶対値）を時間軸方向に積分した量、即ち、読み出し方向傾斜磁場Ｇｒｏの０次モーメントに対応する。
【０１３７】
傾斜磁場が平坦な領域では、０次モーメントは直線状に変化するが、傾斜磁場が平坦ではない領域の０次モーメントは、非線形となる。画像再構成は、サンプリングされたデータがｋ空間上で線形な領域にあることを前提とするため、非線形なサンプリングされたデータをｋ空間上で線形となるように補正することが望まれる。
【０１３８】
図５は、読み出し方向傾斜磁場Ｇｒｏが非線形な領域において等時間間隔でサンプリングされたＭＲ信号が、ｋ空間上では不等間隔になることを示す概念図である。
【０１３９】
図５の上段は、前述のリグリッディング処理用の読み出し方向の位相差データΔθ（ｔ）により、正確に算出された読み出し方向傾斜磁場Ｇｒｏの波形の一例を示す。
即ち、図５の上段において、横軸は、読み出し方向傾斜磁場Ｇｒｏのパルスの印加開始時刻からの経過時間ｔ（上記サンプリング時刻ｔと同じ）を示し、縦軸は読み出し方向傾斜磁場Ｇｒｏの磁場強度を示す。
【０１４０】
図５の中段は、図５の上段の読み出し方向傾斜磁場Ｇｒｏの磁場強度の絶対値を時間積分したものである。積分期間の始期は、共通して、読み出し方向傾斜磁場Ｇｒｏのパルスの印加開始時刻である。従って、図５の中段において、横軸は積分期間の終期時刻を示し、縦軸は、読み出し方向傾斜磁場Ｇｒｏの磁場強度の絶対値の時間積分値、即ち、０次モーメント(0-ORDER MOMENT)を示す。
【０１４１】
図５の下段は、時間的に等間隔でサンプリングされる場合における、１つの位相エンコードステップ分のＭＲ信号（１ラインのＭＲ信号）に対する各サンプリング期間の模式図である。
【０１４２】
図５の下段において、横軸は、上段と同様に、読み出し方向傾斜磁場Ｇｒｏのパルスの印加開始時刻からの経過時間ｔを示し、縦軸は、ＭＲ信号の強度を示す。この例では、周波数エンコードステップ数が２５６の例を示し、２５６のサンプリング期間ＳＰ１，ＳＰ２，ＳＰ３，ＳＰ４，…，ＳＰ２５６が設定される。即ち、１ラインのＭＲ信号は、図５の縦方向の一点鎖線で示すように、２５６のサンプリング期間ＳＰ１〜ＳＰ２５６に等間隔に分割される。
【０１４３】
「傾斜磁場が平坦ではない領域」は、傾斜磁場の０次モーメントが非線形な領域であり、「傾斜磁場が平坦な領域」は、傾斜磁場の０次モーメントが線形な領域である。従って、図５の上段、中段、下段から分かるように、読み出し方向傾斜磁場Ｇｒｏが平坦な領域及び平坦ではない領域を含めて、ＭＲ信号が時間的に等間隔でサンプリングされる場合、生成されるｋ空間データは、ｋ空間上では不等間隔になる。
【０１４４】
読み出し方向傾斜磁場Ｇｒｏの印加の下でサンプリングされたＭＲ信号は、ｋ空間上では、読み出し方向傾斜磁場Ｇｒｏの０次モーメントに対応するところ、この０次モーメントは、図５の中段の横方向の複数の一点鎖線で示すように、互いに不等間隔となるからである。
【０１４５】
なお、以下の説明では、単に「０次モーメント」と言った場合、読み出し方向傾斜磁場Ｇｒｏの０次モーメントを指すものとする。
【０１４６】
図６は、本実施形態におけるリグリッディング処理の第１方法の概念を示す模式図である。図６の上段は図５の上段と同じである。図６の下段は、不等間隔でサンプリングされる場合における、１つの位相エンコードステップ分のＭＲ信号に対する各サンプリング期間の模式図である。
【０１４７】
第１方法では、１ラインのＭＲ信号は、図６の中段及び下段の縦方向の一点鎖線で示すように、２５６のサンプリング期間ＳＰ１’，ＳＰ２’，ＳＰ３’，…，ＳＰ２５６’に不等間隔で分割される。
【０１４８】
上記の不等間隔なサンプリング期間ＳＰ１’〜ＳＰ２５６’の定め方を示すのが、図６の中段である。図６の中段は、図５と同様の０次モーメントを示すが、図中の横方向の複数の一点鎖線のみが図５とは異なる。即ち、０次モーメントが等間隔で上昇するように、横方向の各一点鎖線が引かれている。横方向の各一点鎖線と、０次モーメントを示す太線との各交点を通るように、縦方向の各一点鎖線が引かれている。
【０１４９】
第１方法では、各サンプリング期間ＳＰ１’〜ＳＰ２５６’の各代表時刻を積分期間の終期とする読み出し方向傾斜磁場Ｇｒｏの磁場強度の絶対値の各時間積分値が互いに等間隔となるように、各サンプリング期間ＳＰ１’〜ＳＰ２５６’は定められる。上記の「時間積分値」における積分期間の始期は例えば、共通に、読み出し方向傾斜磁場Ｇｒｏのパルスの印加開始時刻である。また、上記「代表時刻」とは、例えば、各サンプリング期間ＳＰ１’〜ＳＰ２５６’の終了時刻でもよいし、中心時刻でもよい。
【０１５０】
第１方法では、上記のように定められたサンプリング期間に従って、ＭＲ信号が時間的に不等間隔でサンプリングされることで、ｋ空間データが生成される。このようにして生成されるｋ空間データの各マトリクス要素は、ｋ空間で等間隔に配置される。上記「ｋ空間で等間隔に配置される」とは、各サンプリング期間に対応する各０次モーメントの値が、図６の中段の横方向の複数の一点鎖線で示すように、互いに等間隔になることを意味する。
【０１５１】
換言すれば、第１方法では、ＭＲ信号における、各サンプリング期間に対応する部分の収集時刻（受信時刻）での各０次モーメントが等間隔になるように、ＭＲ信号が時間的に不等間隔でサンプリングされる。
【０１５２】
図７は、本実施形態におけるリグリッディング処理の第２方法の概念を示す模式図である。図７の上段は、図５の下段のように時間的に等間隔でサンプリングされることで生成されるｋ空間データの１ライン分の各マトリクス要素のマトリクス値ＭＥ１，ＭＥ２，ＭＥ３，ＭＥ４，…，ＭＥ２５６を示す。
【０１５３】
ここでは一例として、周波数エンコードステップ数を２５６で考えるので、１ライン分のマトリクス要素数も２５６である。各マトリクス値ＭＥ１，ＭＥ２，ＭＥ３，…，ＭＥ２５６は、その上に示すＭＲ信号における各サンプリング期間ＳＰ１，ＳＰ２，ＳＰ３，…，ＳＰ２５６にそれぞれ対応する。
【０１５４】
第２方法では、図７の上段のように、時間的に等間隔でＭＲ信号をサンプリングすることで、一旦ｋ空間データが生成される。この後、ｋ空間データは、各マトリクス要素に対応するサンプリング期間の各代表時刻までの読み出し方向傾斜磁場Ｇｒｏのパルス強度の各時間積分値が等間隔となるように再配列されて、新たなｋ空間データとなる。再配列に際しては補間等の処理を用いればよく、上記「代表時刻」は、第１方法と同様である。
【０１５５】
図７の中段は、再配列後のｋ空間データの各マトリクス値ＭＥ１’，ＭＥ２’，ＭＥ３’，…，ＭＥ２５６’を上側に示し、下側に元のＭＲ信号（図７の上段と同じ）を示す。これらマトリクス値ＭＥ１’〜ＭＥ２５６’は、元のＭＲ信号における各サンプリング期間ＳＰ１’〜ＳＰ２５６’（図６の下段と同じ）の各信号強度からそれぞれ生成されたはずの値である。
【０１５６】
図７の下段は、図６の中段と同じ０次モーメントを示す。図７の下段の縦軸である０次モーメントを等間隔で分割する横方向の複数の一点鎖線に示すように、再配列後のｋ空間データの各マトリクス要素に対応する０次モーメントは、等間隔となる。
【０１５７】
即ち、第２方法では、ＭＲ信号における、各マトリクス要素に対応する部分の各収集時刻（受信時刻）での各０次モーメントが等間隔になるように、ｋ空間データは再配列される。
【０１５８】
＜本実施形態の動作＞
図８は、本実施形態においてリグリッディング処理の第２方法が採用される場合のＭＲＩ装置１０の動作の流れの一例を示すフローチャートである。以下、前述した各図を適宜参照しながら、図８に示すステップ番号に従って、ＭＲＩ装置１０の動作を説明する。
【０１５９】
［ステップＳ１］システム制御部６１（図１参照）は、入力装置７２を介してＭＲＩ装置１０に対して入力された撮像条件に基づいて、ＭＲＩ装置１０の本スキャンの撮像条件の一部を設定する。ここでは一例として、ＥＰＩが本スキャンとして設定されるものとする。また、プレスキャンなどによってＲＦパルスの中心周波数等が設定される。また、例えば、天板２２上の被検体ＰにＲＦコイル装置１００がセットされる。
この後、ステップＳ２に進む。
【０１６０】
［ステップＳ２］天板移動機構２３は、シーケンスコントローラ５８の制御に従って、ガントリ３０内の磁場中心に被検体Ｐの撮像部位が位置するように、天板２２を移動させる。
【０１６１】
次に、システム制御部６１は、位置決め画像のデータ収集をＭＲＩ装置１０に実行させる。具体的には、静磁場電源４２により励磁された静磁場磁石３１によって撮像空間に静磁場が形成される。また、シムコイル電源４４からシムコイルユニット３２に電流が供給されて、撮像空間に形成された静磁場が均一化される。
【０１６２】
そして、システム制御部６１は、パルスシーケンスをシーケンスコントローラ５８に入力する。シーケンスコントローラ５８は、入力されたパルスシーケンスに従って傾斜磁場電源４６、ＲＦ送信器４８及びＲＦ受信器５０を駆動させることで、被検体Ｐの撮像部位が含まれる撮像領域に傾斜磁場を形成させると共に、ＲＦコイルユニット３４からＲＦパルスを発生させる。
【０１６３】
このため、被検体Ｐ内の核磁気共鳴により生じたＭＲ信号は、受信ＲＦコイル２４やＲＦコイル装置１００により検出されて、ＲＦ受信器５０に入力される。ＲＦ受信器５０は、ＭＲ信号に前述の処理を施すことでＭＲ信号の生データを生成し、これら生データを画像再構成部６２に入力する。
【０１６４】
画像再構成部６２は、ＭＲ信号の生データをｋ空間データとして配置及び保存する。画像再構成部６２のフーリエ変換部６２ｂは、ｋ空間データにフーリエ変換を含む画像再構成処理を施すことで画像データを再構成し、得られた画像データを画像データベース６３に保存する。
【０１６５】
画像処理部６４は、画像データベース６３から画像データを取り込み、これに所定の画像処理を施すことで２次元の表示用画像データを生成し、この表示用画像データを記憶装置７６に保存する。この後、システム制御部６１は、表示用画像データが示す位置決め画像を表示装置７４に表示させる。
この後、ステップＳ３に進む。
【０１６６】
[ステップＳ３]システム制御部６１は、表示される位置決め画像に基づいてユーザが選択したＦＯＶなどの撮像条件に応じて、パルスシーケンスを含む本スキャンの残りの撮像条件を設定する。
【０１６７】
次に、システム制御部６１は、前述の２点のみ本スキャンのパルスシーケンスから変更されたテンプレートスキャンＡのパルスシーケンスを設定する（図２の中段参照）。
また、システム制御部６１は、前述の２点のみテンプレートスキャンＡから変更されたテンプレートスキャンＢのパルスシーケンスを設定する（図２の下段参照）。
【０１６８】
また、システム制御部６１は、ＭＲ信号の収集領域がテンプレートスキャンＢの場合よりも読み出し方向に所定間隔ずれたテンプレートスキャンＣのパルスシーケンスを設定する（図４の下段参照）。
そして、システム制御部６１は、テンプレートスキャンＡ〜Ｃ、及び、本スキャンの各パルスシーケンスをシーケンスコントローラ５８（図１参照）に入力する。
この後、ステップＳ４に進む。
【０１６９】
［ステップＳ４］ＭＲＩ装置１０は、以下の＜１＞〜＜４＞のフローに従って、図２及び図４で説明したテンプレートスキャンＡ〜Ｃ、及び、本スキャンを実行する。
【０１７０】
＜１＞システム制御部６１は、ＭＲＩ装置１０の各部を制御することで、位置決め画像のＭＲ信号の収集時と同様にして、テンプレートスキャンＡのパルスシーケンスを実行させる。これにより、テンプレートスキャンＡで収集されたＭＲ信号がサンプリングされて、位相差データ算出部６５においてｋ空間データとして保存される。
【０１７１】
＜２＞上記同様にテンプレートスキャンＢのパルスシーケンスが実行され、これにより収集されたＭＲ信号がサンプリングされ、位相差データ算出部６５においてｋ空間データとして保存される。
【０１７２】
＜３＞上記同様にテンプレートスキャンＣのパルスシーケンスが実行され、これにより収集されたＭＲ信号の生データがサンプリングされ、位相差データ算出部６５においてｋ空間データとして保存される。
【０１７３】
ここで、上記テンプレートスキャンＡ〜Ｃの実行に際して、以下の２点を補足する。
第１に、テンプレートスキャンＣの実行時において、ＲＦ受信器５０における読み出し方向の検波周波数は、テンプレートスキャンＢの実行時よりもΔｆだけシフトされる。これにより、テンプレートスキャンＣにおけるＭＲ信号の収取領域は、テンプレートスキャンＢとは所定間隔（例えば５画素〜１０画素）だけ読み出し方向にずれる。
【０１７４】
第２に、位相誤差補正用の位相補正データdV(xr)、及び、リグリッディング処理用の位相差データΔθ（ｔ）の算出では、前述の例ではｋ空間データの中心ライン（実効エコー時間でのＭＲ信号）のみ用いられる。従って、ｋ空間データの中心ラインに相当するＭＲ信号まで収集すれば十分であり、その後に発生するＭＲ信号は、収集しないでもよい。
【０１７５】
＜４＞上記同様に本スキャンのパルスシーケンスが実行され、これにより収集されたＭＲ信号がサンプリングされ、画像再構成部６２においてｋ空間データとして保存される。
【０１７６】
ＭＲＩ装置１０は、以上の＜１＞〜＜４＞の処理を１スライス（又は１スラブ）のデータ収集として、全スライス（全スラブ）のデータ収集を行う。
【０１７７】
なお、上記＜１＞〜＜４＞の各ＭＲ信号の収集では、一例として、ＭＲ信号を時間的に等間隔でサンプリングすることでｋ空間データが一旦生成される。
【０１７８】
また、上記＜１＞〜＜４＞の順番は一例にすぎず、変更してもよい。
この後、ステップＳ５に進む。
【０１７９】
［ステップＳ５］位相差データ算出部６５は、画像再構成部６２に保存されたテンプレートデータＡ、Ｂと、前述の（３）式等に基づいて、全てのスライス（又はスラブ）に対してそれぞれ、位相誤差補正用の位相補正データdV(xr)を算出する。位相差データ算出部６５は、位相補正データdV(xr)を画像再構成部６２のフーリエ変換部６２ｂに入力する。
【０１８０】
また、位相差データ算出部６５は、画像再構成部６２に保存されたテンプレートデータＢ、Ｃと、前述の（５）〜（１２）式に基づいて、全てのスライス（又はスラブ）に対してそれぞれ、リグリッディング処理用の読み出し方向の位相差データΔθ（ｔ）を算出する。
【０１８１】
画像再構成部６２のリグリッディング処理部６２ａは、位相差データΔθ（ｔ）を位相差データ算出部６５から取得し、位相差データΔθ（ｔ）に基づいて本スキャンの実行時における読み出し方向傾斜磁場Ｇｒｏの波形を算出（再現）する。
この後、ステップＳ６に進む。
【０１８２】
［ステップＳ６］リグリッディング処理部６２ａは、ステップＳ５においてＭＲ信号を等時間間隔でサンプリングすることで生成された本スキャンの各スライス（又はスラブ）のｋ空間データに対してそれぞれ、前述の第２方法（図７参照）によりリグリッディング処理を施す。
【０１８３】
即ち、一旦生成されたｋ空間データの各マトリクス要素に対応するサンプリング期間の各代表時刻を積分期間の終期とする読み出し方向傾斜磁場Ｇｒｏの磁場強度の時間積分値が互いに等間隔となるように、ｋ空間データが再配列され、リグリッディング処理後のｋ空間データとなる。
【０１８４】
フーリエ変換部６２ｂは、リグリッディング処理後のｋ空間データをリグリッディング処理部６２ａから取得し、取得したｋ空間データに対して、前述の位相補正データdV(xr)に基づく位相誤差の補正処理が含まれた画像再構成処理を施す（この処理は、スライス毎又はスラブ毎に実行され、画像再構成処理にはフーリエ変換が含まれる）。
【０１８５】
これにより、位相誤差が補正された本スキャンの画像データがスライス（又はスラブ）毎に生成される。フーリエ変換部６２ｂは、生成した本スキャンの画像データを画像データベース６３に保存する。
【０１８６】
画像処理部６４は、画像データベース６３から画像データを取り込み、これに所定の画像処理を施すことで２次元の表示用画像データを生成し、この表示用画像データを記憶装置７６に保存する。
この後、ステップＳ７に進む。
【０１８７】
［ステップＳ７］システム制御部６１は、表示用画像データが示す画像を表示装置７４に表示させる。以上が図８のフローの説明である。
【０１８８】
図９は、本実施形態においてリグリッディング処理の第１方法が採用される場合のＭＲＩ装置１０の動作の流れの一例を示すフローチャートである。以下、前述した各図を適宜参照しながら、図９に示すステップ番号に従って、ＭＲＩ装置１０の動作を説明する。
【０１８９】
［ステップＳ１１〜Ｓ１３］図８のステップＳ１〜Ｓ３と同じである。この後、ステップＳ１４に進む。
【０１９０】
［ステップＳ１４］ＭＲＩ装置１０は、各スライス（又は１スラブ）に対して、テンプレートスキャンＡ、Ｂ、Ｃを実行し、テンプレートデータＡ、Ｂ、Ｃが位相差データ算出部６５においてｋ空間データとして保存される。本スキャンのパルスシーケンスが実行されない点を除いて、図８のステップＳ４と同様である。
この後、ステップＳ１５に進む。
【０１９１】
［ステップＳ１５］図８のステップＳ５と同様に、位相差データ算出部６５は、リグリッディング処理用の読み出し方向の位相差データΔθ（ｔ）を算出し、これをリグリッディング処理部６２ａに入力する。リグリッディング処理部６２ａは、前述同様に、位相差データΔθ（ｔ）に基づいて、本スキャン実行時の読み出し方向傾斜磁場Ｇｒｏの波形を算出する。
この後、ステップＳ１６に進む。
【０１９２】
［ステップＳ１６］ＭＲＩ装置１０は、各スライス（又は１スラブ）に対して、本スキャンのパルスシーケンスを実行し、収集されたＭＲ信号サンプリングしてｋ空間データを生成し、画像再構成部６２においてｋ空間データを保存する。
このとき、画像再構成部６２のリグリッディング処理部６２ａは、ステップＳ１５で算出された読み出し方向傾斜磁場Ｇｒｏの波形に基づいて、前述の第１方法（図６参照）によりＭＲ信号を時間的に不等間隔でサンプリングすることでｋ空間データを生成する。
【０１９３】
即ち、収集されたＭＲ信号の各ラインは、読み出し方向傾斜磁場の０次モーメントが等間隔となるように、時間的に不等間隔でサンプリングされる。このようにして、リグリッディング処理済のｋ空間データが生成される。
この後、ステップＳ１７に進む。
【０１９４】
［ステップＳ１７］図８のステップＳ５と同様に、位相差データ算出部６５は、位相誤差補正用の位相補正データdV(xr)を算出し、フーリエ変換部６２ｂに入力する。フーリエ変換部６２ｂは、ステップＳ１６で生成されたリグリッディング処理済のｋ空間データに対して、前述のステップＳ６と同様に、位相補正データdV(xr)に基づく位相誤差の補正が含まれる画像再構成処理を施す。フーリエ変換部６２ｂは、これにより生成された本スキャンの画像データを画像データベース６３に保存する。
【０１９５】
画像処理部６４は、上記画像データに所定の画像処理を施すことで２次元の表示用画像データを生成し、この表示用画像データを記憶装置７６に保存する。
この後、ステップＳ１８に進む。
【０１９６】
［ステップＳ１８］システム制御部６１は、表示用画像データが示す画像を表示装置７４に表示させる。
以上が本実施形態のＭＲＩ装置１０の動作説明である。
【０１９７】
＜本実施形態の効果＞
本実施形態のリグリッディング処理の原理を概念的に説明すると、以下になる。
読み出し方向傾斜磁場Ｇｒｏの強度が一定のサンプリング時刻ｔｘ〜ｔｙにおいて読み出し方向の検波周波数をΔｆずらす場合を考える。この場合、各サンプリング時刻において読み出し方向の位相が等間隔にずれ、各サンプリング時刻のＭＲ信号に付与される読み出し方向の位置情報は、等間隔でずれる。このような理想的な場合、リグリッディング処理は不要である。
【０１９８】
しかし、実際には読み出し方向傾斜磁場Ｇｒｏの波形は矩形波ではなく、立ち上がり及び立ち下がりの期間でランプサンプリングされる場合もある上、オーバーシュートやアンダーシュートなどがある。このように読み出し方向傾斜磁場Ｇｒｏの強度が変化している期間に読み出し方向の検波周波数をΔｆずらすと、各時刻での読み出し方向の位相は、「各時刻の読み出し方向傾斜磁場Ｇｒｏの強度」と、Δｆとの積に等価となり、各時刻において読み出し方向の位相が不等間隔でずれる。
【０１９９】
ここで、テンプレートデータＢのみでは、収集されたＭＲ信号が受信ＲＦコイル２４やＲＦコイル装置１００で検出された期間における読み出し方向傾斜磁場Ｇｒｏの波形を再現できない。しかし、読み出し方向の検波周波数がΔｆだけずれたテンプレートスキャンＢ、Ｃの各ｋ空間データを適切に用いれば、両者の差分で磁場不均一性による位相への影響を相殺できるので、読み出し方向の位相差データΔθ（ｔ）を算出できる。
【０２００】
本実施形態では、読み出し方向の位相差データΔθ（ｔ）によって読み出し方向傾斜磁場Ｇｒｏの波形が正確に再現され、このように正確に再現された読み出し方向傾斜磁場Ｇｒｏに基づいてリグリッディング処理が実行される。この結果、リグリッディング処理の精度を向上できる。
【０２０１】
また、本実施形態は、サーチコイル及び積分器による計測、又は、シミュレーション、又は、ＭＲ信号による計測によって読み出し方向傾斜磁場Ｇｒｏを計算するものではない。本実施形態では、ＭＲ信号の収集領域を読み出し方向に空間的シフトさせた２つのテンプレートスキャンＢ、Ｃの各ｋ空間データから読み出し方向傾斜磁場Ｇｒｏを計算するので、計測に時間を要することはなく、容易且つ正確にリグリッディング処理を実行できる。
【０２０２】
また、位相誤差補正用と、リグリッディング処理用とに２つずつのテンプレートスキャンを実行する場合、テンプレートスキャンの数が４となり、撮像時間が長くなる。本実施形態では、テンプレートスキャンＢは、テンプレートスキャンＡと共に位相誤差補正にも用いられる。テンプレートスキャンの数が３で済むため、撮像時間をそれほど長くすることなく、良好な画像が得られる。
【０２０３】
以上説明した実施形態によれば、ＭＲＩにおいて、リグリッディング処理を従来よりも高精度に実行できる。
【０２０４】
＜本実施形態の補足事項＞
［１］上記実施形態では、テンプレートスキャンＢ、Ｃの各ｋ空間データの中心ラインから読み出し方向の位相差データΔθ（ｔ）を算出し、この位相差データΔθ（ｔ）に基づいてリグリッディング処理を実行する例を述べた。
【０２０５】
これは、ｋ空間の中心ラインのＭＲ信号を収集時の読み出し方向傾斜磁場Ｇｒｏの波形の歪みと、ｋ空間の端側の各ラインのＭＲ信号を収集時の読み出し方向傾斜磁場Ｇｒｏの波形の歪みとが同じであるという前提に基づく。本発明の実施形態は、かかる態様に限定されるものではない。
【０２０６】
例えば、傾斜磁場電源４６内の傾斜磁場アンプ（図示せず）などのデバイスの性能が十分ではない場合を考える。この場合、本スキャンの読み出し方向傾斜磁場Ｇｒｏの波形は、時間的に始めの方に収集されるエコー（ＭＲ信号）と、最後の方で収集されるＭＲ信号とで異なることがある。
【０２０７】
このように、ｋ空間のライン毎に、そのラインに該当するエコー（ＭＲ信号）を収集時の読み出し方向傾斜磁場Ｇｒｏの波形が異なる場合、以下のようにリグリッディング処理を実行してもよい。
【０２０８】
具体的には、１の画像用に収集されたｋ空間データに対するリグリッディング処理は、図６、図７で説明したように、このｋ空間データのライン毎に実行できる。即ち、位相エンコードステップ数及び読出しエンコードステップ数が共に２５６のｋ空間データであれば、２５６のライン毎に実行できる。
【０２０９】
そこで、例えば、図４のテンプレートスキャンＢ、Ｃで収集される全てのＭＲ信号を用いて、ｋ空間のライン毎に読み出し方向の位相差データΔθ（ｔ）を算出する。そして、ｋ空間のライン毎に、当該ライン用に算出された位相差データΔθ（ｔ）に基づいてリグリッディング処理（図６、図７参照）を実行する。
【０２１０】
より詳細には、例えば、各テンプレートスキャンＢ、Ｃの始めに収集されるＭＲ信号同士で、ｋ＝−３のライン（仮に位相エンコード方向傾斜磁場Ｇｐｅが印加されていれば位相エンコードステップ−３に相当）の位相差データΔθｋ＝−３（ｔ）を算出する。この位相差データΔθｋ＝−３（ｔ）に基づいて、ｋ空間データにおけるｋ＝−３のラインのリグリッディング処理を実行できる。
【０２１１】
同様に、各テンプレートスキャンＢ、Ｃの２番目に収集されるＭＲ信号同士で、ｋ＝−２のラインの位相差データΔθｋ＝−２（ｔ）を算出し、これに基づいてｋ空間データにおけるｋ＝−２のラインのリグリッディング処理を実行できる。
【０２１２】
以下、ｋ＝−１のラインから、最後に収集されるラインまで、同様の処理を繰り返せばよい。
【０２１３】
［２］本実施形態では、テンプレートスキャンＡ、Ｂ、Ｃ、及び、本スキャンの実行後に、位相誤差補正用の位相補正データ及びリグリッディング処理用の位相差データが算出されてから、位相誤差補正及びリグリッディング処理が含まれる画像再構成が実行される例を述べた。これは、説明を分かり易くするための便宜的な順序にすぎない。
【０２１４】
位相誤差補正用の位相補正データの算出、リグリッディング処理用の位相差データの算出、及び、位相誤差補正及びリグリッディング処理が含まれる画像再構成は、１スライスに対するテンプレートスキャンＡ〜Ｃ及び本スキャンの終了後、他のスライスに対するテンプレートスキャンＡ〜Ｃ及び本スキャンの実行と並行して実行してもよい。
【０２１５】
［３］本実施形態では、２つのテンプレートスキャンＢ、Ｃにより、リグリッディング処理用の位相差データが算出される例を述べた。本発明の実施形態は、かかる態様に限定されるものではない。読み出し方向にＭＲ信号の収集領域を互いにずらした３つのテンプレートスキャン、或いは、４つ以上のテンプレートスキャンにより、リグリッディング処理用の位相差データを算出してもよい。
【０２１６】
例えば、ＭＲ信号の収集領域をテンプレートスキャンＢから読み出し方向にシフトさせる量を、テンプレートスキャンＣ、Ｄ、Ｅ、Ｆ、Ｇ、Ｈにおいてそれぞれ、５画素、６画素、７画素、８画素、９画素、１０画素としてもよい。
【０２１７】
この場合、テンプレートデータＢ、Ｃから算出される第１の位相差データ、テンプレートデータＢ、Ｄから算出される第２の位相差データ、テンプレートデータＢ、Ｅから算出される第３の位相差データ、テンプレートデータＢ、Ｆから算出される第４の位相差データ、テンプレートデータＢ、Ｇから算出される第５の位相差データ、及び、テンプレートデータＢ、Ｈから算出される第６の位相差データの６つの平均により、リグリッディング処理用の最終的な位相差データが算出される。
【０２１８】
［４］上記実施形態では、テンプレートスキャンＡを実行することで、位相誤差の補正を実行する例を述べた。本発明の実施形態は、かかる態様に限定されるものではない。
位相誤差があまり大きくない場合、テンプレートスキャンＡを省略し、位相補正データの算出（位相誤差の補正）を実行しないシーケンスとしてもよい。この場合、テンプレートスキャンの回数を２回に減らせるので、スキャン時間をさらに短縮できると共に、消費電力をさらに低減できる。
【０２１９】
［５］上記実施形態では、本スキャン及びテンプレートスキャンＡにおいて、実効エコー時間でのＭＲ信号が奇数番目のエコーであり、実効エコー時間での読み出し方向傾斜磁場Ｇｒｏの極性がマイナスである例を述べた。
【０２２０】
この場合、前述の理由で、テンプレートスキャンＢ、Ｃでは、実効エコー時間でのＭＲ信号が偶数番目のエコーであり、実効エコー時間での読み出し方向傾斜磁場Ｇｒｏの極性がプラスであることが望ましい。本発明の実施形態は、かかる態様に限定されるものではない。
【０２２１】
本スキャン及びテンプレートスキャンＡにおいて、実効エコー時間でのＭＲ信号が奇数番目のエコーであり、実効エコー時間での読み出し方向傾斜磁場Ｇｒｏの極性がプラスであってもよい。この場合、テンプレートスキャンＢ、Ｃでは、実効エコー時間でのＭＲ信号が偶数番目のエコーであり、実効エコー時間での読み出し方向傾斜磁場Ｇｒｏの極性がマイナスであることが望ましい。
【０２２２】
或いは、本スキャン及びテンプレートスキャンＡにおいて、実効エコー時間でのＭＲ信号が偶数番目のエコーであり、実効エコー時間での読み出し方向傾斜磁場Ｇｒｏの極性がマイナスであってもよい。この場合、テンプレートスキャンＢ、Ｃでは、実効エコー時間でのＭＲ信号が奇数番目のエコーであり、実効エコー時間での読み出し方向傾斜磁場Ｇｒｏの極性がプラスであることが望ましい。
【０２２３】
或いは、本スキャン及びテンプレートスキャンＡにおいて、実効エコー時間でのＭＲ信号が偶数番目のエコーであり、実効エコー時間での読み出し方向傾斜磁場Ｇｒｏの極性がプラスであってもよい。この場合、テンプレートスキャンＢ、Ｃでは、実効エコー時間でのＭＲ信号が奇数番目のエコーであり、実効エコー時間での読み出し方向傾斜磁場Ｇｒｏの極性がマイナスであることが望ましい。
［６］ＭＲ信号の収集時間を短縮すると、ｋ空間データ生成時のサンプリング間隔と対比して、読み出し方向傾斜磁場Ｇｒｏの波形が一定とはみなせない場合が出てくる。その場合、リグリッディング処理により、サンプリングされたＭＲ信号の生データを画像再構成前にｋ空間上で等間隔となるように再配列することが望ましい。
【０２２４】
従って、ＥＰＩのパルスシーケンスに限らず、アンダーシュートやオーバーシュートの部分がある場合などにも、上記実施形態の手法で本スキャンの読み出し方向傾斜磁場の波形を正確に算出し、この傾斜磁場波形に基づいてリグリッディング処理を実行することで画質を向上できる。
【０２２５】
即ち、ＭＲ信号の検出期間中における読み出し方向傾斜磁場Ｇｒｏの波形において、その磁場強度が一定ではない期間が少しでも含まれるパルスシーケンスであれば、上記実施形態の手法を適用可能であり、リグリッディング処理の精度を向上できる。
【０２２６】
［７］請求項の用語と実施形態との対応関係を説明する。なお、以下に示す対応関係は、参考のために示した一解釈であり、本発明を限定するものではない。
【０２２７】
ガントリ３０内の各構成要素、及び、制御装置４０の全体（図１参照）が本スキャンのパルスシーケンスの実行により被検体ＰからＭＲ信号を収集する機能は、請求項記載の本スキャン実行部の一例である。
【０２２８】
ガントリ３０内の各構成要素、及び、制御装置４０の全体がテンプレートスキャンＢの実行によりＭＲ信号を収集し、これをサンプリングしてｋ空間データを生成する機能は、請求項記載の第１収集部の一例である。即ち、テンプレートスキャンＢは、請求項記載の第１パルスシーケンスの一例である。
【０２２９】
ガントリ３０内の各構成要素、及び、制御装置４０の全体がテンプレートスキャンＣの実行によりＭＲ信号を収集し、これをサンプリングしてｋ空間データを生成する機能は、請求項記載の第２収集部の一例である。即ち、テンプレートスキャンＣは、請求項記載の第２パルスシーケンスの一例である。
【０２３０】
ガントリ３０内の各構成要素、及び、制御装置４０の全体がテンプレートスキャンＡの実行によりＭＲ信号を収集し、これをサンプリングしてｋ空間データを生成する機能は、請求項記載の第３収集部の一例である。即ち、テンプレートスキャンＡは、請求項記載の第３パルスシーケンスの一例である。
ＲＦコイル装置１００、及び、受信ＲＦコイル２４は、請求項記載の受信ＲＦコイルの一例である。
【０２３１】
［８］本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０２３２】
１０：ＭＲＩ装置，
２０：寝台ユニット，２２：天板，
３１：静磁場磁石，３２：シムコイルユニット，３３：傾斜磁場コイルユニット，
３４：ＲＦコイルユニット，４０：制御装置，６０：演算装置
【図１】
2015000161.tif 000005
【図２】
2015000161.tif 000006
【図３】
2015000161.tif 000007
【図４】
2015000161.tif 000008
【図５】
2015000161.tif 000009
【図６】
2015000161.tif 000010
【図７】
2015000161.tif 000011
【図８】
2015000161.tif 000012
【図９】
2015000161.tif 000013