(57)【要約】  テーブル連続移動撮影においてエコープレナー法などの高速撮影を、画質の劣化を招くことなく実現し、テーブル移動の高速化、すなわち撮影の高速化を図る。  磁気共鳴撮影装置の磁場発生手段、移動手段および信号処理手段を制御する撮影制御手段は、移動手段を移動させながら、１回の励起用高周波磁場パルス印加後、複数の読み出し傾斜磁場を印加して複数の核磁気共鳴信号の計測を行なう撮影シーケンスを実行する。この際、移動手段の移動に起因して前記複数の核磁気共鳴信号に与えられる読み出し傾斜磁場の位置ずれを補正用データとして算出しておき、計測した核磁気共鳴信号を補正用データで補正する。
2007094174.tif 000002
【特許請求の範囲】
【請求項１】
静磁場、高周波磁場および傾斜磁場の各磁場を発生させる磁場発生手段と、静磁場中で被検体を移動させる移動手段と、前記被検体から発生する核磁気共鳴信号を計測し、被検体の画像を作成する信号処理手段と、前記磁場発生手段、移動手段および信号処理手段とを制御し、前記移動手段を移動させながら、１回の励起用高周波磁場パルス印加後、複数の読み出し傾斜磁場を印加して複数の核磁気共鳴信号の計測を行なう撮影シーケンスを実行する撮影制御手段とを備えた磁気共鳴撮影装置において、
前記撮影制御手段は、前記移動手段の移動に起因して前記複数の核磁気共鳴信号に与えられる読み出し傾斜磁場の位置ずれを算出する算出手段と、前記位置ずれを補正する補正手段を備えたことを特徴とする磁気共鳴撮影装置。
【請求項２】
請求項１記載の磁気共鳴撮影装置において、
前記補正手段は、前記撮影シーケンスを実行するためのパラメータを変更して、前記位置ずれを補正することを特徴とする磁気共鳴撮影装置。
【請求項３】
請求項１記載の磁気共鳴撮影装置において、
前記補正手段は、前記複数の核磁気共鳴信号のそれぞれについて、前記位置ずれを補正することを特徴とする磁気共鳴撮影装置。
【請求項４】
請求項１に記載の磁気共鳴撮影装置において、
前記算出手段は、前記移動手段の移動速度および前記複数の核磁気共鳴信号の計測間隔をもとに、読み出し傾斜磁場の位置ずれを補正する補正値を算出することを特徴とする磁気共鳴撮影装置。
【請求項５】
請求項１に記載の磁気共鳴撮影装置であって、
前記撮影制御手段は、前記撮影シーケンスと同じシーケンス形状であって位相エンコード傾斜磁場パルスの印加を含まない補正用データ取得用シーケンスを実行し、
前記算出手段は、前記補正用データ取得用シーケンスで得られた核磁気共鳴信号から、読み出し傾斜磁場の位置ずれを補正するための補正値を算出することを特徴とする磁気共鳴撮影装置。
【請求項６】
請求項２に記載の磁気共鳴撮影装置であって、
前記補正手段は、前記読み出し傾斜磁場の変化に応じて核磁気共鳴信号を計測する際の参照周波数を補正する周波数補正手段である磁気共鳴撮影装置。
【請求項７】
請求項２に記載の磁気共鳴撮影装置であって、
前記補正する手段は、前記読み出し傾斜磁場の変化に応じて静磁場のオフセット値を補正する静磁場オフセット補正手段である磁気共鳴撮影装置。
【請求項８】
請求項３に記載の磁気共鳴撮影装置であって、
前記補正手段は、前記読み出し傾斜磁場の変化により生じた核磁気共鳴信号の位相変化を補正する位相補正手段である磁気共鳴撮影装置。
【請求項９】
請求項１ないし８いずれか１項に記載の磁気共鳴撮影装置であって、
前記撮影シーケンスは、エコープレナー法による撮影シーケンスであることを特徴とする磁気共鳴撮影装置。
【請求項１０】
請求項９に記載の磁気共鳴撮影装置であって、
前記撮影制御手段は、前記撮影シーケンスと同じシーケンス形状であって位相エンコード傾斜磁場パルスの印加を含まない補正用データ取得用シーケンスを、前記テーブルを移動させることなく実行し、
前記補正手段は、前記補正用データ取得用シーケンスで得られた核磁気共鳴信号から位相補正値を作成し、当該位相補正値を用いて、複数の核磁気共鳴信号のうち奇数番目の信号に生じる位相誤差と偶数番目の信号に生じる位相誤差を補正することを特徴とする磁気共鳴撮影装置。
【請求項１１】
請求項１ないし８いずれか１項に記載の磁気共鳴撮影装置であって、
前記撮影シーケンスは、ファーストスピンエコー法による撮影シーケンスであることを特徴とする磁気共鳴撮影装置。
【請求項１２】
請求項９又は１１に記載の磁気共鳴撮影装置であって、
前記撮影シーケンスは、複数回の励起で全位相エンコードのデータを取得するマルチショット撮影シーケンスであって、前記撮影制御手段は、信号の計測順序を位相エンコード方向に連続するように制御することを特徴とする磁気共鳴撮影装置。
【請求項１３】
請求項９又は１０に記載の磁気共鳴撮影装置であって、
前記撮影シーケンスは、第１の軸および第２の軸の位相エンコードを含む３Ｄエコープレナー法による撮影シーケンスであって、前記撮影制御手段は、前記第１の軸および第２の軸の位相エンコードのうち、移動手段の移動方向に直交する方向の一つの軸の位相エンコードについて、他の軸の全位相エンコードの信号を１回の励起で取得するように計測順序を制御することを特徴とする磁気共鳴撮影装置。
【請求項１４】
静磁場、高周波磁場および傾斜磁場の各磁場を発生させる磁場発生手段と、静磁場中で被検体を移動させる移動手段と、前記被検体から発生する核磁気共鳴信号を計測し、被検体の画像を作成する信号処理手段と、前記磁場発生手段、移動手段および信号処理手段とを制御する撮影制御手段とを備え、
前記撮影制御手段は、前記移動手段を移動させながら、１回の励起用高周波磁場パルス印加後、第１の軸の位相エンコード傾斜磁場を印加するとともに、複数の読み出し傾斜磁場と第２の軸の位相エンコード傾斜磁場を印加して複数の核磁気共鳴信号の計測を行なう３Ｄ撮影シーケンスを実行し、前記第２の軸の位相エンコード傾斜磁場によりエンコードされる全位相エンコードの信号を１回の励起で取得するように計測順序を制御することを特徴とする磁気共鳴撮影装置。
【請求項１５】
請求項１４に記載の磁気共鳴撮影装置であって、
前記第２の軸の位相エンコードは、スライスエンコードであることを特徴とする磁気共鳴撮影装置。
【請求項１６】
請求項１４に記載の磁気共鳴撮影装置であって、
前記撮影制御手段は、前記第１の軸および第２の軸の位相エンコードのうち、位相エンコード数の小さい位相エンコードについて、全位相エンコードの信号を１回の励起で取得するように計測順序を制御することを特徴とする磁気共鳴撮影装置。
【請求項１７】
被検体を搭載した移動手段を移動させながら、１回の励起用高周波磁場パルス印加後、複数の読み出し傾斜磁場を印加して、前記被検体から複数の核磁気共鳴信号の計測を行なうテーブル移動磁気共鳴撮影方法において、
前記移動手段の移動に起因して前記複数の核磁気共鳴信号に与えられる読み出し傾斜磁場の位置ずれを算出するステップ（１）、
前記ステップ（１）で算出した位置ずれから、核磁気共鳴信号毎の位相補正値を算出するステップ（２）、および
計測した複数の核磁気共鳴信号に対し、それぞれ前記位相補正値を用いて位相補正するステップ（３）を含むことを特徴とする磁気共鳴撮影方法。
【請求項１８】
被検体を搭載した移動手段を移動させながら、１回の励起用高周波磁場パルス印加後、複数の読み出し傾斜磁場を印加して、前記被検体から複数の核磁気共鳴信号の計測を行なうテーブル移動磁気共鳴撮影方法において、
前記移動手段の移動に起因して前記複数の核磁気共鳴信号に与えられる読み出し傾斜磁場の位置ずれを算出するステップ（１）、
前記ステップ（１）で算出した位置ずれから、核磁気共鳴信号毎の周波数差を算出するステップ（２’）、および
前記ステップ（２’）で算出された周波数差により、前記複数の核磁気共鳴信号を計測する際の参照周波数を補正するステップ（３’）を含むことを特徴とする磁気共鳴撮影方法。
【請求項１９】
被検体を搭載した移動手段を移動させながら、１回の励起用高周波磁場パルス印加後、複数の読み出し傾斜磁場を印加して、前記被検体から複数の核磁気共鳴信号の計測を行なうテーブル移動磁気共鳴撮影方法において、
前記移動手段の移動に起因して前記複数の核磁気共鳴信号に与えられる読み出し傾斜磁場の位置ずれを算出するステップ（１）、
前記ステップ（１）で算出した位置ずれから、核磁気共鳴信号毎に読み出し傾斜磁場強度の差を算出するステップ（２”）、および
前記複数の核磁気共鳴信号を計測する際に、前記ステップ（２”）で算出された読み出し傾斜磁場強度の差を静磁場のオフセットとして印加するステップ（３”）を含むことを特徴とする磁気共鳴撮影方法。
【請求項２０】
被検体を搭載した移動手段を移動させながら、１回の励起用高周波磁場パルス印加後、複数の読み出し傾斜磁場を印加して、前記被検体から複数の核磁気共鳴信号の計測を行なうテーブル移動磁気共鳴撮影方法であって、
前記計測は、第１の軸および第２の軸の位相エンコード傾斜磁場の印加を含み、
第１の軸および第２の軸の位相エンコードのうちエンコード数の少ない方の全位相エンコードの信号を１回の励起で取得するように計測順序を制御することを特徴とする磁気共鳴撮影方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、核磁気共鳴現象を利用して被検体の検査部位の断層画像を得る磁気共鳴イメージング（以下、ＭＲＩという）装置に関し、特にテーブルを連続移動しながら、限定された装置の視野よりも広い視野の画像を高速で取得可能なＭＲＩ装置に関する。
【背景技術】
【０００２】
ＭＲＩ装置では良好な画像を得ることができる静磁場均一度の高い空間は比較的狭い領域に限定されており、従来、撮像視野もその領域内に限定されていた。近年、このような制約下でも、被検体が載置されたテーブルを撮影中に連続移動させることによって、均一な静磁場空間以上の大きな領域、もしくは全身の画像を取得する方法（以下、テーブル移動撮影と呼ぶ）が提案されている（例えば、特許文献１）。この特許文献１には、テーブルをステップ状或いは連続的に移動させながら、周波数エンコード方向をテーブルの移動方向と同一の方向として信号を取得する動作を繰り返し、最終的にテーブル移動全工程で得られた信号から画像を再構成する手法が開示されている。データ収集法としては、エコープレナー型（ＥＰＩ）での収集も可能であることが記載されている。
【特許文献１】特開２００３−１３５４２９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
一般にＭＲＩにおいてはエコー信号を計測するために読み出し傾斜磁場が利用される。高周波磁場の印加によって励起されたスラブの核磁化は、その収集以前に印加された読み出し傾斜磁場（ディフェイズ傾斜磁場）と、極性が逆で強度が同じ読み出し傾斜磁場（リフェイズ傾斜磁場）が印加されたときに位相が揃い、エコー信号として計測される。テーブル連続移動撮影においてテーブルの移動方向を周波数エンコード方向（読み出し方向）とした場合、ディフェイズ傾斜磁場からリフェイズ傾斜磁場までの間に、テーブルの移動に伴いリフェイズしようとする磁化の位置がずれることとなり、完全に位相が揃わない、すなわち計測される信号は位置ずれに対応する位相回転を含むという問題を生じる。この問題は、１回の励起で１つの信号を計測するグラディエントエコー法等のパルスシーケンスを採用した場合には、計測される全ての信号について同じ条件となるので再構成した画像に与える影響は無視することができる。【０００４】
しかし例えばエコープレナー法のように１回の励起後、読み出し傾斜磁場の反転を利用して複数の信号を取得する高速撮影方法をテーブル移動撮影に適用した場合には、複数の信号毎に印加される読み出し傾斜磁場の条件（信号が収集されるまでの読み出し傾斜磁場の履歴）が異なることになる。このようなＭＲ信号毎の読み出し傾斜磁場履歴のばらつきは、アーチファクトの発生や画質の劣化につながる。【０００５】
本発明は、テーブル移動撮影において１回の励起後に複数の信号を計測し、これら複数の信号を用いて画像再構成する際に、テーブルの移動により、読み出し傾斜磁場の被検体に対する相対的な位置にずれを生じ、そのことにより発生するアーチファクトや画質劣化を抑制することを目的とする。
【課題を解決するための手段】
【０００６】
上記目的を達成するため、本発明の磁気共鳴撮影装置は、テーブルを連続移動させながら、１回の励起後に複数の核磁気共鳴信号を計測する撮影シーケンスを実行するに際し、移動手段の移動に起因して複数の核磁気共鳴信号に与えられる読み出し傾斜磁場の位置ずれを算出する算出手段と、当該位置ずれを補正する補正手段とを備える。
読み出し傾斜磁場の位置ずれを算出する算出手段は、例えば、移動手段の移動速度および複数の核磁気共鳴信号の計測間隔をもとに、読み出し傾斜磁場の位置ずれ或いは当該位置ずれから読み出し傾斜磁場強度の差、周波数差或いは位相差を算出する。或いは、算出手段は、本撮影計測と同じ撮影シーケンスからなり位相エンコードを含まない予備計測シーケンスによって得られた計測データから、読み出し傾斜磁場の位置ずれ或いは当該位置ずれから読み出し傾斜磁場強度の差、周波数差或いは位相差を算出する。【０００７】
補正手段は、前記撮影シーケンスを実行するためのパラメータを変更して、前記位置ずれを補正する。例えば、核磁気共鳴信号を計測する際に印加される読み出し傾斜磁場に、算出手段で算出された読み出し傾斜磁場強度の差と同じ強度の静磁場オフセットを与えて読み出し傾斜磁場を補正する。或いは、算出手段で算出された周波数差を用いて、核磁気共鳴信号を計測するための参照周波数を補正する。
また補正手段は、複数の核磁気共鳴信号のそれぞれについて、前記位置ずれを補正する。例えば、算出手段で算出された位相差を用いて、計測された各核磁気共鳴信号を補正する。【０００８】
これら補正する手段によって、画像再構成前に読み出し傾斜磁場の印加位置のずれによる影響が解消されるので、アーチファクトの発生や画質の劣化を防止できる。【０００９】
また本発明の磁気共鳴撮影装置は、テーブルを連続移動させながら、１回の励起後に第１の軸および第２の軸の位相エンコードを与えて複数の核磁気共鳴信号を計測する３Ｄ撮影シーケンスを実行するに際し、第１の軸および第２の軸の位相エンコードのうち一方の全位相エンコードの信号を１回の励起で計測するように第１の軸および第２の軸の位相エンコードの順序を制御する。【００１０】
撮影制御手段が実行する撮影シーケンスは、例えば、エコープレナー法による撮影シーケンス、ファーストスピンエコー法による撮影シーケンス、複数回の励起で全位相エンコードのデータを取得するマルチショット撮影シーケンス、或いは第１の軸および第２の軸の位相エンコードを含む３Ｄエコープレナー法による撮影シーケンスである。マルチショット撮影シーケンスの場合、撮影制御手段は、信号の計測順序を位相エンコード方向に連続するように制御する。３Ｄエコープレナー法による撮影シーケンスの場合、第１の軸および第２の軸の位相エンコードのうち、移動手段の移動方向に直交する方向の位相エンコードについて、全位相エンコードの信号を１回の励起で取得するように計測順序を制御する。
このように計測順序を制御することによりＫ空間において連続したデータを取得することができので、広い領域のデータをつなげて広視野画像を再構成することができる。【００１１】
また本発明の磁気共鳴撮影装置において、撮影制御手段は、前記撮影シーケンスと同じシーケンス形状であって位相エンコード傾斜磁場パルスの印加を含まない補正用データ取得用シーケンスを実行し、前記補正する手段は、前記補正用データ取得用シーケンスで得られた核磁気共鳴信号から補正用データを作成し、当該補正用データを用いて読み出し傾斜磁場の変化を補正する。この場合、補正用データは関数フィッティングした後、補正に用いてもよい。【００１２】
補正用データ取得用シーケンスを実行することにより得た補正用データを用いることにより、読み出し傾斜磁場の印加位置のずれの補正と同時に、エコープレナー法の撮影パルスシーケンスに生じるN/2アーチファクトも補正することができる。関数フィッティングしたデータを用いることにより、補正用データ取得時のテーブル位置（すなわち被検体の形状）によらず精度のよい補正を行うことができる。
【発明の効果】
【００１３】
本発明によれば、テーブル移動撮影において１回の励起後に複数の信号を計測し、これら複数の信号を用いて画像再構成する際に、テーブルの移動により、読み出し傾斜磁場の被検体に対する相対的な位置にずれを生じ、そのことにより発生するアーチファクトや画質劣化を抑制することができる。
【発明を実施するための最良の形態】
【００１４】
以下、本発明のＭＲＩ装置の実施の形態を説明する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。【００１５】
最初に、本発明のＭＲＩ装置の概略を図１に基づいて説明する。図１(ａ)は本発明のＭＲＩ装置の全体構成を示すブロック図、図１(ｂ)は各要素の詳細を示すブロック図である。この磁気共鳴イメージング装置は、主たる構成として、静磁場発生系２と、傾斜磁場発生系３と、シーケンサ４と、送信系５と、受信系６と、信号処理系７と、テーブル駆動系８と、中央処理装置（ＣＰＵ）９を備えている。ＣＰＵ９は、信号処理系７の処理に必要な演算を行なう演算機能と、シーケンサ４等の装置全体を制御する制御機能を有している。【００１６】
静磁場発生系２は、被検体１の周りの空間にその体軸方向または体軸と直交する方向に均一な静磁場を発生させるもので、被検体１の周りに永久磁石方式または常電導方式あるいは超電導方式の磁場発生手段が配置されている。静磁場発生系２には、静磁場の均一度を向上するために必要に応じてシムコイル(図示せず)が備えられる。【００１７】
傾斜磁場発生系３は、Ｘ，Ｙ，Ｚの３軸方向に巻かれた傾斜磁場コイル３１と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源３２とから成り、シ−ケンサ４からの命令に従ってそれぞれのコイルの傾斜磁場電源３２を駆動することにより、Ｘ，Ｙ，Ｚの３軸方向の傾斜磁場を被検体１に印加する。具体的には、Ｘ，Ｙ，Ｚの１ないし複数の傾斜磁場を組み合わせて印加することにより、所望の方向にスライス方向傾斜磁場パルス（Ｇｓ）を印加して被検体1に対するスライス面を設定し、スライス方向に直交する２つの方向に位相エンコード方向傾斜磁場パルス（Ｇｐ）と周波数エンコード方向傾斜磁場パルス（Ｇｆ）を印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。【００１８】
本発明においても、これら傾斜磁場パルスの軸は装置座標系（Ｘ，Ｙ，Ｚ）に対し任意であるが、周波数エンコード方向がテーブルの移動方向と同一である場合に本発明は好適である。【００１９】
シーケンサ４は、高周波磁場パルス（以下、ＲＦパルスという）と傾斜磁場パルスを所定のパルスシーケンスで繰り返し印加する撮影制御手段で、ＣＰＵ９の制御で動作し、被検体1の断層画像のデータ収集に必要な種々の命令を送信系５、傾斜磁場発生系３、受信系６およびテーブル駆動系８に送る。シーケンサ４の動作については後に詳述する。【００２０】
送信系５は、被検体１の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるためにＲＦパルスを照射するもので、高周波発振器５１と変調器５２と高周波増幅器５３と送信側の高周波コイル５４とから成る。高周波発振器５１から出力された高周波パルスをシーケンサ４からの指令によるタイミングで変調器５２により振幅変調し、この振幅変調された高周波パルスを高周波増幅器５３で増幅した後に被検体１に近接して配置された高周波コイル５４に供給することにより、電磁波（ＲＦパルス）が被検体１に照射される。【００２１】
受信系６は、被検体１の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号（ＮＭＲ信号）を検出するもので、受信側の高周波コイル６１と、増幅器６２と、直交位相検波器６３と、Ａ／Ｄ変換器６４とから成る。送信側の高周波コイル５４から照射された電磁波によって誘起される被検体１の応答の電磁波（ＮＭＲ信号）が被検体１に近接して配置された高周波コイル６１で検出され、増幅器６２で増幅された後、シーケンサ４からの指令によるタイミングで直交位相検波器６３により直交する二系統の信号に分割され、それぞれがＡ／Ｄ変換器６４でディジタル量に変換されて、信号処理系７に送られる。【００２２】
信号処理系７は、ＣＰＵ９の一部、光ディスク７１、磁気ディスク７２等の外部記憶装置と、ＣＲＴ等からなるディスプレイ７３と、操作部７４とを有し、受信系６からのデータがＣＰＵ９に入力されると、ＣＰＵ９が信号処理、画像再構成等の処理を実行し、その結果である被検体１の断層画像をディスプレイ７３に表示すると共に、外部記憶装置の磁気ディスク７２等に記録する。操作部７４は、ＭＲＩ装置の各種制御情報や信号処理系７で行う処理の制御情報を入力するもので、トラックボール又はマウス７４１、キーボード７４２などから成る。【００２３】
テーブル駆動系８は、被検体１を載せるテーブル８１と、テーブル８１を垂直方向および水平方向に移動させるテーブル移動機構８２と、シーケンサ４或いはＣＰＵ９からの制御によりパルスシーケンスと関連してテーブルの移動を制御するテーブル制御部８３とを備えている。またテーブル移動機構８２には、テーブル８１の位置を検出する位置検出部（図示せず）が備えられている。テーブル制御部８３は、この位置検出部からの位置情報を利用してテーブルの移動制御を行なうとともに、ＭＲ信号に必要な補正量などの情報を信号処理系７に与える。【００２４】
本実施の形態では、静磁場発生系２が発生する均一度の高い静磁場空間（撮影空間）２１に対し、テーブル８１を移動させて被検体１の広い範囲（上記撮影空間より広い範囲）を撮影する場合を説明するが、静磁場発生系２とテーブル８１との関係は相対的であり、静磁場発生系２自体が移動機構を有し、停止するテーブルに対し移動する装置であってもよい。【００２５】
シーケンサ４は、予め設定されたパルスシーケンスに従い、傾斜磁場発生系３、送信系５、受信系６、信号処理系７およびテーブル駆動系８を制御し、テーブルを連続的に移動させながら、撮影空間に位置する被検体１の特定領域の励起とその領域からの信号の収集を繰り返す。パルスシーケンスは、撮影方法によって異なる種々のパルスシーケンスが、ＣＰＵ９の記憶部に格納されており、操作部７４を介して所望のパルスシーケンスを選択するとともに撮影に必要なパラメータを設定するようになっている。本発明で実行されるパルスシーケンスは、高周波磁場パルスの印加による１回の励起後、複数の信号を計測するマルチエコーシーケンスであり、複数の信号を計測する間隔（エコー間隔）とその間に移動するテーブルの移動量に応じた信号の補正が組み込まれる。【００２６】
次に上記構成のＭＲＩ装置の動作を説明する。
まず本発明の理解をしやすくするために、図１０および図１１を参照して、従来のテーブル連続移動撮影の概略を説明する。図１０（a）は従来のテーブル連続移動撮影で採用されるパルスシーケンスの一例を示す図である。このパルスシーケンスは公知のグラディエントエコー法のパルスシーケンスで、1回のＲＦパルス照射201につき1個のエコー信号206を取得する。ＲＦパルス照射201から１個のエコー信号取得までを繰り返し時間間隔ＴＲ(208)で繰り返して画像再構成に必要な全エコー信号を取得する。図中、Ｇｓ、Ｇｐ、Ｇｒはそれぞれスライス傾斜磁場、位相エンコード傾斜磁場、読み出し傾斜磁場（周波数エンコード傾斜磁場）の印加を示すラインであり、202はスライス傾斜磁場パルス、203は位相エンコード傾斜磁場パルス、204は読み出し傾斜磁場パルス、205は、ＡＤ変換を行なうサンプリングウィンドウである。一般のＭＲＩではこれら傾斜磁場は任意の軸に設定することができるが、図１１に示すテーブル連続移動撮影では、読み出し傾斜磁場をテーブル移動方向（ｘ方向）とし、それと直交する２方向をスライス方向、位相エンコードとしている。スライス方向、位相エンコード方向は、ｘ方向と直交する方向であれば任意であるが、図１１（a）、（b）では、ｚ方向がスライス方向である場合、すなわちテーブル８１に対して水平に撮影断面が設定された場合を示している。後に説明する本発明の実施の形態でも、スライス方向および位相エンコード方向を図１１と同じ方向にする場合を説明するが、それに限定されない。【００２７】
図１０（a）のパルスシーケンスを、全位相エンコードのエコー信号206が得られるまで繰り返すことにより、図１０(b)に示すようにＫ空間（Ky-Kx空間）を埋めるデータ206-1〜206-11が収集計測される。ここでテーブル８１を移動しながら撮影した場合には、静磁場中心に対する被検体１の撮影位置は、撮影開始の位置301からテーブル移動方向（図１１における矢印Ａ方向）と逆方向（図１１における矢印Ｂ方向）に移動する。これにより読み出し傾斜磁場は、被検体１が停止しているときに印加される部位よりずれて印加され、その結果、全位相エンコードの信号を読み出し方向にフーリエ変換したデータは、図１１(c)に示すように、テーブル移動方向と逆方向にずれたデータとなる。【００２８】
ここでテーブルの移動速度ＶをＶ≦FOVx÷（TR×N）（ただしFOVxは読み出し方向の視野サイズ、Ｎは位相エンコード数）としておくことにより、図１１(d)に示すように、読み出し方向に連続したハイブリッドデータ（Ky-x空間データ）を得ることができる。このハイブリッドデータを位相エンコード方向にフーリエ変換することにより、本来の視野FOVよりも読み出し方向に広い視野の画像を作成することができる。【００２９】
このように従来のテーブル連続移動撮影では、テーブル移動方向を読み出し方向とすることにより、比較的速い速度で広い視野を撮影できるという利点があるが、読み出し方向に連続したデータを得るために、テーブル移動速度はパルスシーケンスの繰り返し時間ＴＲと繰り返し回数Ｎによる制約がある。一方、高速パルスシーケンスとしてＲＦパルスによる１回の励起（１ショット）で複数の信号を収集するＥＰＩなどのパルスシーケンスが知られている。このようなＥＰＩ等の高速パルスシーケンスをテーブル連続移動撮影に適用することによりテーブル移動速度の高速度化、撮影の高速化を図ることが可能となる。しかし、単純にＥＰＩ等のパルスシーケンスを組み合わせた場合には、画質が劣化する。これは１回のRF照射後に複数のエコー信号を計測する際にエコー間でもテーブルは移動し、本来読み出し傾斜磁場によってリフェイズすべき核磁化がリフェイズしないままエコーを計測することによる。そこで本発明では、テーブル連続移動撮影にＥＰＩ等のシングルショットマルチエコーの撮影方法を適用すると共に、読み出し傾斜磁場がずれた位置に印加されることによる影響を実質的に排除する手段を提供する。【００３０】
＜第１の実施の形態＞
以下、本発明の第１の実施の形態として、パルスシーケンスがマルチショットＥＰＩ（エコープレナーイメージング）の場合を例にしてＭＲＩ装置の動作を説明する。図２に動作の手順を示す。本実施の形態では、位相補正値を算出するステップ101、本計測データを計測するステップ102、本計測データの位相を補正するステップ103、補正後の本計測データを読み出し方向にフーリエ変換するステップ104、読み出し方向のデータを結合するステップ105、予め設定したテーブル移動範囲の最終位置までデータを収集した後に画像を作成するステップ106、107とからなる。【００３１】
これら動作を実現するための信号処理系（ＣＰＵ９）の構成を図１２に示す。図示するように、ＣＰＵ９は、撮影制御部９０と、受信系６で計測したエコー信号に対しフーリエ変換や補正計算などを施し画像を再構成する画像演算部９１と、各エコーの位相補正値を算出する位相算出部９２と、位相算出部９２が算出した位相補正値を保存する補正値メモリ部９３と、主制御部９４とを備えている。【００３２】
撮影制御部９０は、シーケンサ４を介して撮影動作を制御する。位相算出部９２は、補正用位相算出ステップ101において、シーケンサ４に設定されたパルスシーケンスのパラメータ及びテーブル制御部８３に設定されたテーブル移動速度などの情報を用いて、ＥＰＩ計測を行なった際にテーブル移動に起因して各エコーに生じる位相ずれ（位相補正値）を算出する。画像演算部９１は、テーブル移動撮影により得られたエコーを用いて画像再構成する（ステップ104及びステップ105）。その際、位相補正ステップ103を実行し、位相算出部９２が算出した位相補正値を用いて各エコーを補正する。主制御部９４は、撮影制御部９０、画像演算部９１、位相算出部９２および補正値メモリ部９３の動作を制御する。【００３３】
これら補正用位相算出ステップ101および位相補正ステップ103以外のステップは従来のテーブル移動撮影とほぼ同じである。まず本計測データ計測ステップ102について説明する。【００３４】
図３は、本計測で行なうテーブル移動撮影の概要を説明する図であり、被検体１と撮影空間301の関係を示している。この実施の形態では、被検体１はその体軸がテーブル８１の移動方向（x方向とする）に対し平行になるようにテーブル８１上に載置されている。テーブル８１を移動させることにより撮影空間301と被検体1との位置関係が変わるので、被検体１の異なる部位の画像を取得できる。【００３５】
本計測で採用するパルスシーケンスの一例を図４に示す。このパルスシーケンスは、テーブルの定速移動下で実行されることを除き、公知のＥＰＩシーケンスと同じである。図中、ＲＦ、Ｇｓ、Ｇｐ、Ｇｒ、ＡＤ/echoはそれぞれ、ＲＦパルス、スライス傾斜磁場、位相エンコード傾斜磁場、周波数エンコード（読み出し）傾斜磁場、Ａ/Ｄ変換（サンプリングウィンドウ）/エコー信号の軸を表している。【００３６】
簡単に説明すると、まず検知する磁化を含む被検体にＲＦパルス401を照射すると同時にスライス（スラブ）選択のためのスライス傾斜磁場402を印加し、画像化するスライスを選択する。次いで、位相エンコードのオフセットを与える位相エンコード傾斜磁場403と読み出し傾斜磁場のオフセットを与える読み出し傾斜磁場405を印加し、位相エンコード傾斜磁場パルス404を離散的に印加しながら反転する読み出し傾斜磁場406の各周期内でエコー信号群407を取得する（図示する例では６つのエコー信号を計測している）。図示するシーケンスを一定の繰り返し時間(ＴＲ)408で、位相エンコード傾斜磁場403の大きさを変化させながら所定の回数（全位相エンコード数÷１回の励起で収集される信号の数）繰り返し、全位相エンコードのエコー信号を収集する。ここで重要なことは、Ｋ空間データの取得順序を位相エンコード方向のデータ配列順序と一致させることである。つまり、Ｋ空間上の端から順々に、間引きを行うことなくデータを取得していくことが望ましいと考えられる。即ち、位相エンコード傾斜磁場404の大きさは位相エンコードの１ステップ分に調整され、位相エンコード傾斜磁場403の変化（オフセット）は、１回の励起で収集される信号の数×１ステップ分に調整される。これにより、位相エンコード方向に連続するＫ空間データを取得することができる。【００３７】
またこのパルスシーケンスの実行中におけるテーブルの移動速度Ｖは、全位相エンコードのデータを収集する時間Ｔ（＝TR×N/M：ただしＮは全位相エンコード数、Ｍは１ショットで計測される信号数）にテーブルが移動する距離が、読み出し方向（ｘ方向）のＦＯＶの長さFOVxと等しいかそれ以下であるように設定されている。すなわちＶ≦FOVx÷Ｔである。この関係は、テーブル移動速度をパルスシーケンスの実行時間（例えば、繰り返し時間）に対応して制御、あるいは逆に、パルスシーケンスの実行時間（例えば、繰り返し時間）をテーブル移動速度に対応して制御することによって実現する。具体的には、操作部７４を介してパルスシーケンスのパラメータ（図４の繰り返し時間408やエコー間隔412など）がシーケンサ４に設定され、パルスシーケンスが決定されると、テーブル制御部８３は、設定されたパラメータからテーブル移動速度を算出し、算出されたテーブル移動速度でテーブル８１を移動するようにテーブル移動機構８２を制御する。或いは操作部７４を介してテーブル移動速度がテーブル制御部８３に設定されると、ＣＰＵ９は設定されたテーブル移動速度に基づきパルスシーケンスの繰り返し時間を算出しシーケンサ４に設定する。【００３８】
図５は、このようにして取得したエコー信号をＫ空間に配置した結果を示した図である。図５(a)はKy-Kx空間データ、図５(b)は(a)のデータを読み出し方向にフーリエ変化した後のデータ（ハイブリッドデータ）である。図では、単純化し、１回の繰り返しで計測される信号の数が４つで繰り返し回数が３回の場合を示している。図５(a)に示すように、本実施の形態では、矢印の方向に計測が進み、1回の繰り返し(=ショット)時間408内で取得した４つのエコー信号群は、Ｋ空間の下側から上側に向って配置される。そして、ショット毎に与える位相エンコード傾斜磁場のオフセット403の大きさが異なることから、各ショットでＫ空間の開始位置が異なり、複数のショット408（408-1〜408-3）によりＫ空間を埋めることができる。ここで、テーブル移動方向に加えられる読み出し傾斜磁場（周波数エンコード傾斜磁場）は、テーブル移動方向と逆方向にシフトして加えられるのと同じ結果となり、図５（b）に示すようなハイブリッドデータ407’-1〜407’-3となる。図５中、408-1、408-2、408-3はそれぞれ図４のショット408に対応し、エコートレイン407-1、407-2、407-3はそれぞれ図４の信号407に対応している。なお、図５（ｂ）では、各ハイブリッドデータのｘ方向のシフト量が、同一ショット内のエコー毎のｘ方向シフト量より多い場合を示しているが、図４のパルスシーケンスにおいて、ＴＲあるいはショット毎の励起位置を制御することにより、図５（ｂ）における各ハイブリッドデータのｘ方向のシフト量を、同一ショット内のエコー毎のシフト量と同じにすることも可能である。この場合には、図１１（ｃ）と同様のハイブリッド空間配置となる。【００３９】
次にこうして計測した本計測データの位相を補正するステップ103および補正ステップ103に用いるための位相補正値を算出するステップ101について説明する。【００４０】
まずエコープレナー法にテーブル連続移動撮影を適用した場合の影響を、図６を参照して説明する。図６(a)は図４の読み出し傾斜磁場軸Grとエコー信号部Echoを拡大したものである。エコープレナー法では、読み出し傾斜磁場パルス406の極性を反転することから、テーブルの移動方向が一定でも、読み出し傾斜磁場の極性によってその影響が異なる。すなわち、図６(a)に示す期間412-1で正の読み出し傾斜磁場406-1-1が印加されたとする。このときの被検体601と傾斜磁場との位置関係を、図６（b）に示す。次の期間412-2では負の読み出し傾斜磁場406-1-2が印加される。ここでテーブルの移動方向が読み出し傾斜磁場の軸方向(x)が一致し、テーブルが図の左から右へ移動したとすると、テーブルとともに被検体601は移動し、被検体601と傾斜磁場との位置関係は、図６（c）に示すように、図６（b）に示す位置よりもｘ方向にずれることになる。この位置のずれ602は、期間412の長さに比例する。【００４１】
一般に、エコープレナー法で連続的にエコー信号を取得できるのは、読み出し傾斜磁場の反転により、読み出し傾斜磁場方向の横磁化を繰り返しリフェイズできるからである。従って、被検体の各位置(x)で、傾斜磁場パルス406-1-1と406-1-2で加えられた磁化のモーメントが０となる必要がある。すなわち傾斜磁場パルス406-1-2は図６(c)の点線で示す位置に印加されなければならない。しかし、撮影中に被検体が移動した場合、被検体の位置(x)が変わることから、傾斜磁場パルス406-1-2は図６(c)に点線で示す位置よりずれて印加されることとなり、傾斜磁場パルス406-1-1と406-1-2を印加した後の磁化のモーメントが０とならない。この結果、取得したエコー信号407-1-1と407-1-2との間には位相の差が生じる。同様にエコー信号407-1-1と407-1-3との間には、期間412の２倍の長さに比例する位置ずれ（602+603）に対応した位相の差が生じる。このような位相の差があるまま画像再構成すると、位相エンコード方向のフーリエ変換時にアーチファクトを生じる。【００４２】
本実施の形態では、位相算出部９２において、テーブルの移動によりエコー信号に生じる位相の差を予め補正値として求めておき（図２、ステップ101）、この補正値を用いて本計測時に取得したエコー信号407から位相の差を除去する（ステップ103）。補正値は、以下のようにして算出される。テーブルの移動速度をＶ、エコー信号取得間隔412をIETとしたとき、IETの間に移動する移動量ΔＸは、式（１）で表される。
ΔＸ=Ｖ×IET      （１）
一方、読み出し傾斜磁場の出力Gr(mT/m)は、読み出し傾斜磁場方向の撮影視野FOVxと、帯域BW、磁気回転比γを用いて、式（２）で表され、
Gr=BW/(γ×FOVx)    （２）
位置のシフトΔＸによって生じる磁場出力の差は、式（３）で表される。
Gr×ΔＸ=BW/(γ×FOVx)×ΔＸ  （３）
これを角周波数の差Δωに換算すると式（４）となる。
Δω=γ×Gr×ΔＸ=BW×ΔＸ/FOVx    （４）
この位置ずれによって生じた角周波数の差Δωは、式（５）で表されるエコー信号の位相θとなる。
θ=Δω×Δｔ MOD 2π        （５）
（式(５)中、Δｔは読み出し傾斜磁場パルスの印加時間、a MOD bはaをbで除した余りを表す。）【００４３】
このようにテーブル移動速度Ｖ、エコー信号取得間隔412の長さIETおよび読み出し傾斜磁場の条件（出力、印加時間）がわかっていれば、式（１）〜（５）から、位置ずれによって生じる位相差を算出することができる。これら数値Ｖ、IET、読み出し傾斜磁場条件は、パルスシーケンスが決まれば決まる値であるので、ステップ101では設定されたパルスシーケンスを基に補正値を算出する。算出された補正値は、例えばＣＰＵ９内の補正値メモリ部９３に記憶される。【００４４】
画像演算部９１は、このように補正用位相算出ステップ101で求めた位相（補正値）を用いて、取得した本計測データの位相を補正する（図２、本計測データ位相補正ステップ103）。次いで、位相補正後の本計測データを読み出し方向にフーリエ変換し（ステップ104）、得られるデータ（x,ky）を位置xに対応して画像用の空間（ハイブリッド空間：横軸が位置x、縦軸が位相エンコードＫｙ）に配置する（ステップ105）。エコープレナー法では、図４に示すように、各ショットでエコー信号取得開始までの時間411と、エコー信号取得後から次のＲＦパルス照射401までの時間413があるが、この間もテーブルは速度Vで移動し続けるので、その時間(即ち411と413の合計)に移動したテーブルの位置を求め、正しい位置にデータを配置する。【００４５】
上述の本計測ステップ102から読み出し方向データ結合ステップ105までを、最終データを取得するまで繰り返し行い、画像再構成に必要な全てのデータを得る（ステップ106）。既に述べたように、テーブルの移動速度Ｖは読み出し方向（ｘ方向）のＦＯＶの長さFOVxを全位相エンコードのデータを収集する時間Ｔで除した値以下に設定されているので、ステップ106で得られるハイブリッドデータは、図１１（d）に示したのと同様に、テーブル移動方向に連続したデータとなる。【００４６】
最終的に結合されたハイブリッド空間データから、必要に応じて、位相エンコード方向のデータが揃わないｘ方向両端のデータを除去した後、位相エンコード方向にフーリエ変換することにより、装置で規定されるＦＯＶより広い範囲の画像を作成する（ステップ106、107）。【００４７】
本実施の形態によれば、テーブル移動撮影にエコープレナー法を適用することにより、テーブルの移動速度を、従来法に比べ１/M（Ｍは１ショットで計測される信号の数）だけ高速化でき、しかも、その際、生じるエコー信号間の位相差を予め算出された補正値で補正することにより、アーチファクトの低減された広視野画像を得ることができる。【００４８】
なお本実施の形態では、エコー信号取得中（サンプリングウィンドウ内）の位置変化は考慮されていない。即ち、エコー信号を取得している間にもテーブル移動により位置は変化するので、サンプリングされた時系列データのデータ間にもわずかであるが位相ずれは生じている。しかしエコープレナー法では、画像の歪みが生じないよう、エコー間隔IETは通常数msに設定されるので、テーブル移動速度が100mm/s(=0.1mm/ms)と非常に速い場合であっても、IETの間に移動する距離は1mm以下である。MRIの一般的な空間分解能は、0.5mm程度であるので、撮影中に位置が移動しても、エコー内のデータは劣化しないと考えられる。【００４９】
しかしながら、エコー信号取得中の位置変化も除去する場合には、厳密解として、式（３）に時間の変数の考えを取り入れた式（３’）を用い、式（３’）〜式（５’）に基づいて計算すればよい。
Gr×ΔＸ(t)=BW/(γ×FOVx)×ΔＸ(t)          （３’）
Δω(t)=γ×Gr×ΔＸ(t)=BW×ΔＸ(t)/FOVx    （４’）
θ(t)=∫Δω(t)dt MOD 2π                    （５’）【００５０】
また以上の実施の形態では、１回の励起後に収集される複数のエコー信号のうち１番目のエコーを基準として１番目のエコーからの位相差を算出し、位相補正する場合を説明したが、より好適には、１番目のエコーについても位相補正する。すなわち、読み出し傾斜磁場のオフセットを与えるパルス405から１番目の読み出し傾斜磁場パルス406までの経過時間を求め、その間に移動したテーブルの距離より、1番目のエコー信号計測時の位相も算出し、これを補正値として１番目のエコーを補正する。これにより、より精度の高い補正を行なうことができる。【００５１】
また以上の実施の形態では、補正用位相算出ステップ101で求めた位相を用いて、計測後の信号を補正する場合を説明したが、このように事後的に補正するのではなく、本計測のエコー計測時の受信条件あるいは読み出し傾斜磁場条件を変更することにより、各エコーに、テーブル移動に起因して生じる誤差を低減することが可能である。この場合には、図２の補正用位相算出ステップ101は、周波数差算出ステップ或いは静磁場オフセット算出ステップとなる。【００５２】
具体的には、位相算出部９２は、テーブルの移動速度Ｖ、エコー信号取得間隔IET、読み出し傾斜磁場の出力Gr、読み出し傾斜磁場方向の撮影視野FOVx、帯域BWおよび磁気回転比γを用いて、前述した式（１）〜（４）から周波数差Δωを算出する（周波数差算出ステップ）。次に本計測データ計測ステップ102でエコー信号を受信する際に、エコー信号を検出時の参照周波数を、上記式（４）で算出したΔωだけずらして受信する。これによっても、テーブルの移動に起因する画質劣化を低減するも可能である。【００５３】
あるいは位相算出部９２は、前述した式（３）により磁場出力の差GrΔＸを算出する。次に本計測データ計測ステップ102でエコー信号を検出する際（読み出し傾斜磁場406の印加時）に上記式（３）で求めた磁場出力の差に相当する静磁場オフセットを印加する。このような静磁場オフセットは、例えば静磁場発生系がシムコイルを備える場合には、シムコイルを駆動することにより実現できる。この場合にも、テーブルの移動に起因する画質劣化を低減するが可能になる。
ただし、上述した３つの補正方法、すなわち、事後的に位相を補正する方法、エコー信号の検出時の参照周波数をずらして受信する方法、静磁場オフセットを印加する方法は、単独で用いてもよいし、２つ或いは３つを組み合わせて実行してもよい。【００５４】
さらに上記実施の形態では、本発明をマルチショットＥＰＩに適用した場合を説明したが、本発明は１回の励起（シングルショット）で全位相エンコードの信号を計測するシングルショットＥＰＩにも同様に適用できることは言うまでもない。また上記実施の形態は２Ｄ撮影を例に説明したが、３Ｄ撮影についても適用することができる。【００５５】
＜第２の実施の形態＞
以下、本発明の第２の実施の形態として、３Ｄ撮影の場合を３Ｄ−ＥＰＩを例に取り説明する。
３Ｄ−ＥＰＩのパルスシーケンスを図７に示す。３Ｄ−ＥＰＩのパルスシーケンスでは、スライス傾斜磁場方向(Kz)と位相エンコード傾斜磁場方向(Ky)に、それぞれエンコードパルスを用いる。従来の３Ｄ−ＥＰＩでは、スライスエンコードパルスと位相エンコードパルスの印加順序は、任意に設定することができ、図４の２Ｄ−ＥＰＩのパルスシーケンスに対し、スライス軸（Ｇｓ）に403と同様のＴＲごとに印加量が変化するスライスエンコードループを追加するのが一般的である。
これに対し、本発明のテーブル連続移動撮影では、同一のスライスエンコードで計測される、位相エンコードの異なる複数の信号が、同じ読み出し傾斜磁場の条件となるように、スライスエンコードパルスと位相エンコードパルスの印加順序を設定する。【００５６】
即ち、本実施の形態では、図７（a）に示すように、スライス選択傾斜磁場701に対して、エコートレイン407毎に傾斜磁場パルス702を用いることで、1回のショットでスライス(Kz)方向に必要な全てのエコー信号を取得する。１回のショットで得られたエコー信号のKz-Ky空間データを図７(b)に示す。ショット毎にスライス方向(Kz)のデータ列703を取得し、各ショットで位相エンコード方向(Ky)のオフセットを与えるパルス403の量を変えながらシーケンスを繰り返し実行し、Ky方向にデータ703-1,703-2,703-3,…を取得する。繰り返しによって得られたエコー信号のKz-Ky空間データを図７(c)に示す。図７(c)の703-1、703-2・・・は、それぞれ１ショットで得られたエコー信号を示し、Kz方向に縦に並ぶ複数の点は、それぞれ一つのエコー信号に対応する。ただし、図７(c)は図１１におけるz-Ky平面に垂直な方向から見た平面図である。
このパルスシーケンス実行中、テーブル移動速度Ｖは、Ｖ≦FOVx/（Ｎ×TR）（ただし、FOVxはテーブル移動方向（ｘ方向）の視野サイズ、Ｎは位相エンコード数、TRは繰り返し時間408である）とする。【００５７】
このようにパルスシーケンスにおけるスライスエンコードパルスと位相エンコードパルスの印加順序を制御することにより、同一のスライスエンコードの信号についてみると、全位相エンコードの信号が同じ読み出し傾斜磁場の条件のときに計測されることになる。即ち、第１スライスエンコードでは、全位相エンコードの信号がエコートレインの最初に計測され、第２スライスエンコードでは、全位相エンコードの信号がエコートレイン407の２番目の信号として計測され、第３スライスエンコードでは、全位相エンコードの信号がエコートレインの３番目の信号として計測される（以下、同様）。取得したデータをKz方向にフーリエ変換した後のKy-x空間の配置は図１１(c)と同様になる。【００５８】
このように３Ｄ撮影では、データの計測順序を制御すると共に、テーブル移動速度Ｖを上式のように制御することにより、テーブル移動方向に連続したデータを得ることができ、また同じスライスエンコードにおいてはエコー信号毎に読み出し傾斜磁場の条件が同じとなり位置ずれの問題が生じないので、位相補正を行なわなくても良好な画像を得ることができる。【００５９】
ここでスライスエンコード方向については、１ショットで計測された信号にはテーブル移動による位置シフトにより傾斜磁場強度差が生じているが、３Ｄ撮影では、同一スライスエンコード内で傾斜磁場強度差をなくすことによりテーブル移動撮影時の画質を改善することができる。また各スライスエンコード同士の信号についても読み出し傾斜磁場の条件を揃えることにより、更に画質を向上させることができる。この場合の撮影の手順は図２に示すフローにおいて本計測データ計測ステップ102が３Ｄ−ＥＰＩであることを除き、第１の実施の形態と同様であり、１ショットで得られる複数のエコー信号について第１の実施の形態と同様の補正を行なう。【００６０】
すなわち、エコー信号取得間隔412の間に移動したテーブル（被検体）の移動量ΔＸから、式（３）〜（５）または式（３’）〜（５’）に基づき位相差を求めておき、取得した信号を補正する。この場合にも、第１の実施の形態と同様に、読み出し傾斜磁場のオフセットを与えるパルス405から１番目の読み出し傾斜磁場パルス406までの経過時間を求め、その間に移動したテーブルの距離より、1番目のエコー信号計測時の位相も算出し、これを補正値として１番目のエコーを補正してもよい。【００６１】
また位相補正値を用いて取得後の信号を補正する代わりに、位相差に対応して信号受信時の参照周波数をシフトするか静磁場オフセットを与えることにより、エコー信号取得間隔412の被検体の移動によりエコー信号にずれて与えられる読み出し傾斜磁場406の影響を補正する。その後、取得したデータを読み出し方向にフーリエ変換し、最終データまで結合し、スライス方向および位相エンコード方向にフーリエ変換して３Ｄ−画像データを作成する。【００６２】
本実施の形態によれば、テーブル連続移動撮影においても、データの連続性を損なうことなく、高速に３次元のデータを取得できる。
なお、３次元撮影では、任意の方向をスライス傾斜磁場方向、位相エンコード傾斜磁場方向、読み出し傾斜磁場方向とすることができるが、本実施の形態では、テーブル移動方向（＝読み出し傾斜磁場方向）に直交する２つの軸のうち、データ取得点数の少ない方向をスライス方向として選択するのが良い。【００６３】
＜第３の実施の形態＞
次に、ＦＳＥ（高速スピンエコー）法によるパルスシーケンスを用いた場合の実施の形態を説明する。本実施の形態においても、撮影の手順は図２に示すフローにおいて本計測データ計測ステップ102がＦＳＥシーケンスであることを除き、第１の実施の形態と同様である。【００６４】
ＦＳＥシーケンスでは、図８(a)に示すようにＲＦパルス801で磁化を励起後、時間間隔802で反転ＲＦパルス803-1〜803-3を連続的に印加する間に、読み出し傾斜磁場パルス804-1〜804-3と位相エンコード傾斜磁場パルス805-1〜805-6を所定の間隔で印加することにより、連続的にエコー信号806-1〜806-3を発生させる。【００６５】
このパルスシーケンスでも、各エコー信号の計測時点における読み出し傾斜磁場方向(X)に与えた傾斜磁場パルス強度の積分量は、信号毎に異なる。撮影中のテーブル移動がない場合およびテーブル移動がある場合について、読み出し傾斜磁場方向(X)に与えた傾斜磁場パルス強度の積分量を図８(b)、(c)に示す。図中、808（808-1〜808-3）は反転RFパルス803の印加タイミング、809（809-1〜809-3）はエコー信号のピーク位置となる時刻を示している。撮影中のテーブル移動がない場合には、(b)に示すように、エコー信号のピーク位置の時間809では、読み出し傾斜磁場方向(X)の傾斜磁場の積分量は０となる。これに対し、テーブル移動がある場合は、(c)に示すように、テーブル位置の変化に伴って傾斜磁場の印加位置が変化するため、本来のエコー信号ピーク時809で傾斜磁場の積分量が０とならない。このため、エコー信号のピーク位置がシフトする。このことは、図８(b)に示すタイミング809において積分量を０とするために与えられるべき読み出し傾斜磁場より大きい又は小さい強度の傾斜磁場が与えられたことを意味する。【００６６】
この傾斜磁場強度の差（積分量を０とするために与えられるべき傾斜磁場強度と実際に印加された傾斜磁場強度との差）は、最初に計測されるエコー806-1を基準にした場合には、エコー信号取得間隔とテーブル移動速度から求められる位置ずれΔＸを用いて前述した式（３）から算出することができる。読み出し傾斜磁場のオフセットを与えるパルス804-0から最初に計測されるエコー806-1までの経過時間を用いた場合には、その間のテーブル移動距離から１番目のエコー806-1についても傾斜磁場強度の差が算出できる。こうして算出した傾斜磁場強度の差から、さらに前述した式（３）〜（５）又は式（３’）〜式（５’）を用いて位相差を算出し、この位相差を用いて計測したエコー信号を信号毎に補正する。補正後の信号を用いて、読み出し方向にフーリエ変換し、データ結合した後、位相エンコード方向にフーリエ変換し画像を作成することは第１および第２の実施の形態と同じである。
なお本実施の形態でも、計測した信号を補正する代わりに、式（４）の角周波数差を用いて信号受信時の参照周波数を補正するか式（３）の静磁場オフセットを与えて、傾斜磁場強度の差によるピークずれを解消した信号を収集してもよい。【００６７】
本実施の形態によれば、テーブル連続移動撮影にＦＳＥシーケンスを適用し、高速に且つアーチファクトのない広領域画像を得ることができる。【００６８】
＜第４の実施の形態＞
次に第４の実施の形態として、補正値の取得方法が異なる実施の形態を説明する。以上説明した実施の形態では、予め補正値を撮影条件（パルスシーケンスのパラメータやテーブル移動速度）から算出する場合を説明したが、本実施の形態では補正値を予備計測によって取得し、それを計測した信号の補正或いは信号受信時の周波数補正或いは静磁場オフセット補正に用いる。
本実施の形態においても、信号処理系の構成は、図１２に示す構成と同じであるが、シーケンサ４が予備計測のシーケンスを実行すること、位相算出部９２がシーケンサ４に設定されたパラメータではなく、予備計測によって得られた信号を用いて補正値を算出することが異なる。以下、予備計測により補正値を求める手法を説明する。【００６９】
図９(a)に手順を示す。本実施の形態は、本計測データ計測102の前に、補正用データ計測ステップ108および補正用位相算出ステップ109があることが、図２に示す手順と異なっている。補正用データ計測ステップ108では、本計測と同様にテーブルを連続に移動しながら、補正用データ取得シーケンスを実行して補正用データを取得する。補正用データ取得シーケンスには、本計測シーケンスと同じシーケンス形状で、位相エンコード（３Ｄシーケンスの場合には位相エンコードとスライスエンコード）傾斜磁場軸のパルスのみを出力しないシーケンスを用いる。本計測シーケンスが、図４のパルスシーケンスである場合には、位相エンコード傾斜磁場403を用いずエコー信号を計測する。この場合、テーブルを連続移動しながら１組のデータ（図４の例では６個の信号）を取得し、この１組のデータを補正用データとしてもよいし、テーブルを連続移動しながら所定の間隔で複数回の補正用データ計測を行ない、複数のテーブル位置の補正用データを取得してもよい。【００７０】
補正用位相算出ステップ109では、補正用データ取得シーケンスで計測したエコー信号から、テーブル移動に起因する読み出し傾斜磁場の位置ずれに相当する位相回転を算出し、補正値とする。ここで所定の間隔で補正用データを取得した場合には、本計測データを計測した位置と最も近い位置で取得した補正用データを選択して使用する。これにより精度のよい補正を行なうことができる。しかし１組のデータのみ（例えば頭部撮影位置によるデータ）を取得した場合には、そのデータから算出した位相を異なるテーブル位置（例えば腹部撮影位置）で取得した本計測データに用いて補正すると、正しく補正できない。そこで、算出した位相を関数でフィッティングし、それを補正値とする。これにより１組のデータであっても被検体の存在する範囲の影響を受けず位相補正を行なうことができる。【００７１】
補正用位相算出ステップ109で算出した補正値を用いて、本計測後のステップ103で本計測データの位相補正を行なうこと、或いは本計測データ受信時に参照周波数或いは静磁場オフセットを位相補正量に応じて変化させることは第１の実施の形態と同様である。【００７２】
本実施の形態では、本計測および補正用データ計測のパルスシーケンスが、読み出し傾斜磁場を反転しながら時系列的に複数のエコー信号を取得するエコープレナー法のシーケンスの場合には、本発明の目的であるテーブル移動に起因する読み出し傾斜磁場の位置ずれの補正と併せて、N/2アーチファクトの補正を行なうことができる。【００７３】
N/2アーチファクトは、傾斜磁場パルスの不完全性や出力応答等の要因によって、パルスシーケンスの計算と実際の出力との間に誤差が生じ、サンプリングの中心とエコー信号の発生位置がずれた場合に発生する。Ｋ空間上でエコー信号のピーク位置が中央からずれると、エコープレナー法のシーケンスでは、読み出し方向のフーリエ変換後の偶数番目のエコー、奇数番目のエコーで位相の変化が逆になるため、位相エンコード方向のフーリエ変換時に画像にアーチファクトが生じる。このアーチファクトがN/2アーチファクトである。本実施の形態において、本計測シーケンスと同じ形状の補正用データ取得シーケンスで得られた補正用データには、N/2アーチファクトの原因である位相誤差成分も含まれる。従ってこの補正用データを用いることにより、テーブル移動に起因する読み出し傾斜磁場の位置ずれの補正と同時にN/2アーチファクトの補正を行なうことが可能となる。【００７４】
以上説明したように本実施の形態によれば、本計測とは別にテーブルを連続移動しながら信号を取得し、その位相値を算出して補正用データとすることにより、精度のよい補正を行なうことができる。またテーブル移動に起因する読み出し傾斜磁場の位置ずれの補正と同時に、N/2アーチファクトの補正を行なうことができる。特に、所定の間隔で補正用データを取得した場合には、本計測データを計測した位置と最も近い位置で取得した補正用データを選択して使用することにより、被検体形状や位置に応じてより精度の高い補正を行なうことができる。また１組の補正用データを取得した場合でも、関数フィッティングすることにより、テーブル位置の相違（被検体の形状の違い）の影響を受けることなく補正を行なうことができる。【００７５】
＜第５の実施の形態＞
次に第５の実施の形態として、予備計測によりN/2アーチファクト補正用データのみを取得する実施の形態を説明する。
上記第４の実施の形態では、テーブルを連続移動しながら補正用データ計測を行なっているので、テーブル移動に起因する読み出し傾斜磁場の位置ずれ(位相回転)の補正用データには必然的にN/2アーチファクト補正用データが含まれることになる。しかし、N/2アーチファクト補正用データは別途取得することも可能である。本実施の形態では、N/2アーチファクト補正用データのみを予備計測により取得する。その動作手順を、図９(b)に示す。このフローはステップ118、119の内容が第４の実施の形態（図９(a)のフロー）のステップ108,109と異なる。他のステップは同様である。
本実施の形態では、補正用データ計測ステップ118で、テーブルを停止した状態でN/2アーチファクト補正用データを得るためのエコープレナー法のパルスシーケンスを実行し、N/2アーチファクト補正用データを取得する。このパルスシーケンスは、本計測のパルスシーケンスと形状が同じで、位相エンコード傾斜磁場軸のパルスのみを出力しないシーケンスである。本計測のパルスシーケンスは、シングルショットＥＰＩ、マルチショットＥＰＩ、３Ｄ−ＥＰＩのいずれでもよい。【００７６】
補正用位相算出ステップ119では、第１或いは第２の実施の形態のステップ101と同様に、例えば前述の式（３）〜（５）又は式（３’）〜式（５’）を用いてテーブル移動に起因する読み出し傾斜磁場の位置ずれ(位相回転)の補正用データ算出する。次いで本計測を行ない、本計測データを補正用位相算出ステップ119で算出された補正用データおよび補正用データ計測ステップ118で得られたN/2アーチファクト補正用データを用いて補正する。【００７７】
本実施の形態においても、N/2アーチファクト用の補正用データは関数でフィッティングしたものを用いることが好ましい。一般にN/2アーチファクトは、傾斜磁場パルスの不完全性や出力応答等が主な要因であり、テーブルの位置（被検体位置）に依存しないと考えられるが、関数フィッティングしたデータを用いることにより、N/2アーチファクト補正用データを取得したテーブル位置と本計測データを取得したテーブル位置とが異なることによる補正精度の低下を防止できる。
以上、説明したように、本発明の各実施の形態によれば、テーブル連続移動撮影においてエコープレナー法などの高速撮影を画質の劣化を招くことなく実現することができ、テーブル移動の高速化、すなわち撮影の高速化を図ることができる。【００７８】
本発明は上記実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々に変形して実施できる。本発明は、グラディエントエコーとファーストスピンエコーを組み合わせたGRASE（Gradient echo And Spin Echo）といったシーケンスにも適用できることは言うまでもない。
【図面の簡単な説明】
【００７９】
【図１】本発明が適用される磁気共鳴撮影装置の全体構成を示す図
【図２】本発明の第1の実施の形態の動作を示すフロー
【図３】本発明におけるテーブル移動撮影を説明する図
【図４】第1の実施の形態で採用するパルスシーケンスを示す図
【図５】図４のパルスシーケンスにより取得されたデータのＫ空間配置を示す図
【図６】（a）は図４のパルスシーケンスの一部を示す図、（b）〜（d）は読み出し傾斜磁場の位置ずれを説明する図
【図７】（a）は第２の実施の形態で採用するパルスシーケンスを示す図、（b）および(c)は計測順序を示す図
【図８】（a）は第３の実施の形態で採用するパルスシーケンスを示す図、（b）および(c)はテーブル停止撮影時およびテーブル移動撮影時の読み出し傾斜磁場の印加積分量を示す図
【図９】（a）は本発明の第４の実施の形態の動作を示すフロー、（b）は本発明の第５の実施の形態の動作を示すフロー
【図１０】従来のテーブル移動撮影で採用するパルスシーケンスの一例を示す図
【図１１】従来のテーブル移動撮影を説明する図
【図１２】本発明の信号処理手段（ＣＰＵ）の構成例を示す図
【符号の説明】
【００８０】
１・・・被検体、２・・・静磁場発生系、３・・・傾斜磁場発生系、４・・・シーケンサ、５・・・送信系、６・・・受信系、７・・・信号処理系、８・・・テーブル駆動系、９・・・中央処理装置。
【図１】
2007094174.tif 000003
【図２】
2007094174.tif 000004
【図３】
2007094174.tif 000005
【図４】
2007094174.tif 000006
【図５】
2007094174.tif 000007
【図６】
2007094174.tif 000008
【図７】
2007094174.tif 000009
【図８】
2007094174.tif 000010
【図９】
2007094174.tif 000011
【図１０】
2007094174.tif 000012
【図１１】
2007094174.tif 000013
【図１２】
2007094174.tif 000014
【国際調査報告】
2007094174.tif 000015
2007094174.tif 000016
2007094174.tif 000017
2007094174.tif 000018
2007094174.tif 000019
2007094174.tif 000020
2007094174.tif 000021