JPH1075940A - 位相シフト測定方法、位相シフト補正方法およびｍｒｉ装置 - Google Patents

Abstract

(57)【要約】 【課題】 エンコード勾配等に起因する渦電流や残留磁
化の影響による位相シフトを正確に測定する。 【解決手段】 励起パルスＲを送信し、反転パルスＰ１
を送信し、エンコード勾配ｇｙｎを位相軸に印加し、通
常のリード勾配の前半分のリード勾配ｇｘｒをリード軸
に印加し、補償用パルスpypnを組み込んだリワインド勾
配gyrnを位相軸に印加し、反転パルスＰ２を送信し、デ
フェーザ勾配ｇywdnを位相軸に印加し、リードアウト勾
配gywnを位相軸に印加しながらエコーecho２からデータ
を収集し、そのデータを１次元フーリエ変換して得られ
た位相データを基にエンコード勾配ｇｙｎやリワインド
勾配gyrnの影響によるエコーecho２の位相シフト量を求
める。求めた位相シフト量から補正値を決定し、それに
より補償用パルスpypnを更新する。これをエンコード勾
配ｇｙｎを変えて、必要な回数だけ繰り返す。 【効果】 エンコード勾配等に起因する渦電流や残留磁
化の影響による位相シフトを高精度に補正でき、イメー
ジの画質を改善できる。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、位相シフト測定方
法、位相シフト補正方法およびＭＲＩ（Ｍagnetic Ｒes
onance Imaging）装置に関する。さらに詳しくは、先行
するエンコード勾配等に起因する渦電流や残留磁化の影
響による後続のエコーの位相シフトを測定する位相シフ
ト測定方法、その位相シフトを補正する位相シフト補正
方法およびそれらの方法を実施するＭＲＩ装置に関す
る。

【０００２】

【従来の技術】スプリットエコートレイン（split echo
train）法は、一つのエコートレイン中の複数のエコー
を前半エコーと後半エコーとに分け、前半エコーから第
１のイメージ（proton weighted image）を生成し、後
半エコーから第２のイメージ（T2 weighted image）を
生成するものである。図３７は、スプリットエコートレ
イン法のパルスシーケンス図である。このパルスシーケ
ンスでは、まず、励起パルスＲとスライス勾配ｓｓを印
加する。次に、前半の第１の反転パルスＰ１とスライス
勾配ｓｓを印加する。次に、エンコード勾配gy1iを位相
軸に印加する。次に、リード勾配ｇｘｗを印加しなが
ら、エコーecho１からＮＭＲ信号を受信する。その後、
前記エンコード勾配gy1iと時間積分値が等しく逆極性の
リワインド勾配ｇyr1iを位相軸に印加する。なお、ｉは
繰り返し番号であり、ｉ＝１，２，…，Ｉ（たとえばＩ
＝２５６）である。

【０００３】次に、前半の第２の反転パルスＰ２とスラ
イス勾配ｓｓを印加し、エンコード勾配gy2iを位相軸に
印加し、リード勾配ｇｘｗを印加しながらエコーecho２
からＮＭＲ信号を受信し、その後、前記エンコード勾配
gy2iと時間積分値が等しく逆極性のリワインド勾配ｇyr
2iを位相軸に印加する。

【０００４】次に、後半の第１の反転パルスＰ’１とス
ライス勾配ｓｓを印加し、エンコード勾配ｇy'1iを位相
軸に印加し、リード勾配ｇｘｗを印加しながらエコーec
ho’１からＮＭＲ信号を受信し、その後、前記エンコー
ド勾配ｇy'1iと時間積分値が等しく逆極性のリワインド
勾配gyr'1iを位相軸に印加する。

【０００５】次いで、後半の第２の反転パルスＰ’２と
スライス勾配ｓｓを印加し、エンコード勾配ｇy'2iを位
相軸に印加し、リード勾配ｇｘｗを印加しながらエコー
echo’２からＮＭＲ信号を受信し、その後、前記エンコ
ード勾配ｇy'2iと時間積分値が等しく逆極性のリワイン
ド勾配gyr'2iを位相軸に印加する。

【０００６】このように前半の反転パルスＰの送信から
リワインド勾配ｇyrの印加までをＮ／２（Ｎは、エコー
トレイン数。図３７ではＮ＝４）回反復し、次いで後半
の反転パルスＰ’の送信からリワインド勾配gyr'の印加
までをＮ／２（図３７ではＮ＝４）回反復した後、横磁
化をスポイル（ｓpoil）するためのキラー勾配Ｋilを印
加する。なお、図３７では、図示の都合上、Ｎ＝４とし
たが、一般には、Ｎ＞４とする。さらに、上記パルスシ
ーケンスを繰返し時間ＴＲでｉ＝１，２，…，Ｉについ
て実行し、前半のエコーechoから図３８に示すＫ−空間
ＫSPのデータを収集し、後半のエコーecho’から図３８
に示すＫ−空間ＫSP’のデータを収集する。エンコード
勾配ｇｙによってＫ−空間ＫSPにおけるデータの位相軸
方向の位置が決まり、エンコード勾配ｇｙ’によってＫ
−空間ＫSP’におけるデータの位相軸方向の位置が決ま
る。位相軸のエンコード数は、Ｉ（Ｎ／２）となる。そ
して、Ｋ−空間ＫSPのデータから第１のイメージを作成
し、Ｋ−空間ＫSP’のデータから第２のイメージを作成
する。

【０００７】図３９は、従来の高速ＳＥ（Ｓpin Ｅch
o）法のパルスシーケンスである。このパルスシーケン
スでは、まず、励起パルスＲとスライス勾配ｓｓを印加
する。次に、第１の反転パルスＰ１とスライス勾配ｓｓ
を印加する。次に、エンコード勾配gy1iを位相軸に印加
する。次に、リード勾配ｇｘｗを印加しながら、第１エ
コーecho１からＮＭＲ信号を受信する。その後、前記エ
ンコード勾配gy1iと時間積分値が等しく逆極性のリワイ
ンド勾配ｇyr1iを位相軸に印加する。なお、ｉは繰返し
番号であり、ｉ＝１，２，…，Ｉ（例えばＩ＝１２８）
である。

【０００８】次に、第２の反転パルスＰ２とスライス勾
配ｓｓを印加し、エンコード勾配gy2iを位相軸に印加
し、リード勾配ｇｘｗを印加しながら第２エコーecho２
からＮＭＲ信号を受信し、その後、前記エンコード勾配
gy2iと時間積分値が等しく逆極性のリワインド勾配ｇyr
2iを位相軸に印加する。

【０００９】次に、第３の反転パルスＰ３とスライス勾
配ｓｓを印加し、エンコード勾配gy3iを位相軸に印加
し、リード勾配ｇｘｗを印加しながら第３エコーecho３
からＮＭＲ信号を受信し、その後、前記エンコード勾配
gy3iと時間積分値が等しく逆極性のリワインド勾配ｇyr
3iを位相軸に印加する。

【００１０】次に、第４の反転パルスＰ４とスライス勾
配ｓｓを印加し、エンコード勾配gy4iを位相軸に印加
し、リード勾配ｇｘｗを印加しながら第４エコーecho４
からＮＭＲ信号を受信し、その後、前記エンコード勾配
gy4iと時間積分値が等しく逆極性のリワインド勾配ｇyr
4iを位相軸に印加する。

【００１１】このように反転パルスＰの送信からリワイ
ンド勾配ｇyrの受信までをＮ（図３９ではＮ＝４）回反
復した後、キラー勾配Ｋilを印加する。なお、図３９で
は、図示の都合上、Ｎ＝４としたが、実際は、Ｎ＞４と
するのが普通である。さらに、上記パルスシーケンスを
繰返し時間ＴＲでｉ＝Ｉまで繰り返し、エコーechoによ
り図４０に示すＫ−空間ＫSPのデータを収集する。エン
コード勾配ｇｙによってＫ−空間ＫSPにおけるデータの
位相軸方向の位置が決まる。位相軸のエンコード数はＩ
・Ｎとなる。そして、Ｋ−空間ＫSPのデータからイメー
ジを作成する。

【００１２】図３７に示すスプリットエコートレイン法
のパルスシーケンスでは、繰返し時間ＴＲを短縮するた
めに、エンコード勾配ｇｙ，ｇｙ’の波形の振幅を大き
くし、それに合せて時間幅ｔgyを短くしている。また、
図３９に示す高速ＳＥ法のパルスシーケンスでも、繰返
し時間ＴＲを短縮するために、エンコード勾配ｇｙの波
形の振幅を大きくし、それに合せて時間幅ｔgyを短くし
ている。ところが、勾配パルスの振幅を大きくし時間幅
を短くするほど、勾配コイルの周辺の導体に渦電流を生
じる。そして、渦電流を生じると、それらの影響により
エコー毎の位相シフトが大きくなる。そして、位相シフ
トは、イメージ上に位相軸方向のゴーストとなって現わ
れ、アーチファクトを発生させる。このようなゴースト
を低減させる従来技術として、プリスキャンでエンコー
ド勾配を印加せずにプリスキャンデータを収集し、その
プリスキャンデータをＫ−空間の周波数軸方向に１次元
フーリエ変換して位相データを得て、その位相データを
基にリード勾配や反転パルスＰのオフセット位相を補正
してイメージ作成用の本スキャンを行なう方法（「R.Sc
ott Hinks, et al., Proc.ＳＭＲ, 1995, p634」参照）
が提案されている。また、イメージ作成用の本スキャン
を行った後に、エンコード勾配を印加せずに補正用デー
タを収集し、画像再構成計算時に前記補正用データを基
に位相補正処理を行う方法（「Xin Wan, Dennis L.Park
er, etal., Mag. Reso. in Med., 34, pp632-638, 199
5」参照）も提案されている。

【００１３】

【発明が解決しようとする課題】永久磁石型ＭＲＩ装置
においても、上記渦電流による位相シフトの問題がある
が、渦電流の外に残留磁化による位相シフトの問題も生
じる。すなわち、パルスの振幅を大きくするほど、例え
ば整磁板に残留磁化を生じ、この残留磁化が位相シフト
に大きな影響を与える。しかし、上述の従来技術では、
エンコード勾配を印加せずに補正用データを収集するた
め、その補正用データは残留磁化の影響を含んでおら
ず、残留磁化の影響による位相シフトを補正することが
できない問題点がある。そこで、本発明の第１の目的
は、先行するエンコード勾配等に起因する渦電流や残留
磁化の影響による後続のエコーの位相シフトを測定する
位相シフト測定方法を提供することにある。また、本発
明の第２の目的は、その位相シフトを補正して、エンコ
ード勾配等に起因する渦電流や残留磁化の影響によるイ
メージの画質の劣化を防止することが出来る位相シフト
補正方法を提供することにある。さらに、本発明の第３
の目的は、上記位相シフト測定方法および位相シフト補
正方法を実施しうるＭＲＩ装置を提供することにある。

【００１４】

【課題を解決するための手段】第１の観点では、本発明
は、励起パルスを送信し、反転パルスを送信し、エンコ
ード勾配を位相軸に印加し、リード勾配をリード軸に印
加し、リワインド勾配を位相軸に印加することを１回以
上行い、続いて、反転パルスを送信し、デフェーザ（de
phaser）勾配を位相軸に印加し、リードアウト（read o
ut）勾配を位相軸に印加しながらエコーからデータを収
集し、そのデータを１次元フーリエ変換して得られた位
相データを基に前記エンコード勾配の影響による位相シ
フト量を求めることを特徴とする位相シフト測定方法を
提供する。

【００１５】上記第１の観点による位相シフト測定方法
では、リード勾配を印加せずに、デフェーザ勾配および
リードアウト勾配を位相軸に印加しながらエコーからデ
ータを収集する。このデータは、Ｋ−空間における位相
軸に沿った軌跡（trajectory）上のデータになってい
る。そこで、このデータを１次元フーリエ変換して得ら
れた位相データから１次の位相シフト（Ｋ−空間の位相
軸に沿って位相シフト量が１次の変化を示す位相シフト
成分）を正確に知ることが出来る。従って、エンコード
勾配，リード勾配およびリワインド勾配を印加した後、
これを行えば、エンコード勾配等に起因する渦電流や残
留磁化の影響による位相シフトを正確に測定できること
になる。

【００１６】第１の変形として、本発明は、上記第１の
観点による位相シフト測定方法において、前記デフェー
ザ勾配と前記リードアウト勾配の極性を反転してデータ
を収集し、そのデータを１次元フーリエ変換して得られ
た位相データと極性を反転する前に得られた位相データ
とを基に前記エンコード勾配の影響による位相シフト量
を求めることを特徴とする位相シフト測定方法を提供す
る。デフェーザ勾配によっても渦電流が生じるので、そ
の影響を除去しないと本スキャンの場合と異なる測定結
果が出かねない。ところで、デフェーザ勾配の極性を逆
にしても、それによる渦電流の作用は１次の位相シフト
については同方向に作用する。一方、先行するエンコー
ド勾配等による位相シフトは、デフェーザ勾配の極性を
逆にすると、１次の位相シフトについては逆方向に作用
する。そこで、デフェーザ勾配の極性を反転したときの
データと反転していないときのデータの差をとれば、デ
フェーザ勾配による渦電流の影響を除去でき、先行する
エンコード勾配等による位相シフトを正確に測定でき
る。

【００１７】第２の観点では、本発明は、励起パルスを
送信し、反転パルスを送信し、エンコード勾配を位相軸
に印加し、リード勾配をリード軸に印加し、リワインド
勾配を位相軸に印加することを２回以上繰り返し、続い
て、反転パルスを送信し、デフェーザ勾配を位相軸に印
加し、リードアウト勾配を位相軸に印加しながらエコー
から第１のデータを収集し、リフェーザ（rephaser）勾
配を位相軸に印加し、続いて、反転パルスを送信し、デ
フェーザ勾配を位相軸に印加し、リードアウト勾配を位
相軸に印加しながらエコーから第２のデータを収集し、
前記第１のデータを１次元フーリエ変換して得られた位
相データと前記第２のデータを１次元フーリエ変換して
得られた位相データとを基に前記エンコード勾配の影響
による位相シフト量を求めることを特徴とする位相シフ
ト測定方法を提供する。エンコード勾配が“０”でも位
相シフトは“０”でなく、オフセットの位相シフト分が
存在する。このオフセットの位相シフト分は、第１のデ
ータでも第２のデータでも同方向である。一方、第１の
データと第２のデータを収集する間には反転パルスがあ
るので、先行するエンコード勾配等による位相シフト
は、第１のデータと第２のデータとで逆方向になる。そ
こで、第１のデータと第２のデータの差をとれば、オフ
セットの位相シフト分を除去でき、先行するエンコード
勾配等による位相シフトを正確に測定できる。

【００１８】第２の変形として、本発明は、上記第２の
観点による位相シフト測定方法において、前記デフェー
ザ勾配と前記リードアウト勾配と前記リフェーザ勾配の
極性を反転して第１のデータおよび第２のデータを収集
し、それらのデータを１次元フーリエ変換して得られた
各位相データと極性を反転する前に得られた各位相デー
タとを基に前記エンコード勾配の影響による位相シフト
量を求めることを特徴とする位相シフト測定方法を提供
する。デフェーザ勾配とリフェーザ勾配によっても渦電
流が生じるので、その影響を除去しないと本スキャンの
場合と異なる測定結果が出かねない。ところで、デフェ
ーザ勾配とリフェーザ勾配の極性を逆にしても、それに
よる渦電流の作用は１次の位相シフトについては同方向
に作用する。一方、先行するエンコード勾配等による位
相シフトは、デフェーザ勾配の極性を逆にすると、１次
の位相シフトについては逆方向に作用する。そこで、デ
フェーザ勾配とリフェーザ勾配の極性を反転したときの
データと反転していないときのデータの差をとれば、デ
フェーザ勾配とリフェーザ勾配による渦電流の影響を除
去でき、先行するエンコード勾配等による位相シフトを
正確に測定できる。

【００１９】第３の変形として、本発明は、上記第１の
観点から第２の観点の変形までの構成の位相シフト測定
方法において、スティミュレーティッドエコー（stimul
atedecho）を除去するためのクラッシャ（crusher）勾
配をスライス軸に印加することを特徴とする位相シフト
測定方法を提供する。位相シフトを測定したいエコーに
対して、スティミュレーティッドが重なると、位相シフ
トを正確に測定する妨げとなる。そこで、クラッシャ勾
配をスライス軸に印加してスティミュレーティッドエコ
ーを除去すれば、位相シフトを正確に測定できるように
なる。

【００２０】第３の観点では、本発明は、励起パルスを
送信し、反転パルスを送信し、エンコード勾配を位相軸
に印加し、リード勾配をリード軸に印加し、リワインド
勾配を位相軸に印加することを１回以上行い、続いて、
反転パルスを送信し、位相軸には勾配を印加せず、リー
ド勾配をリード軸に印加しながらエコーからデータを収
集し、そのデータを１次元フーリエ変換して得られた位
相データを基に前記エンコード勾配の影響による０次の
位相シフト量を求めることを特徴とする位相シフト測定
方法を提供する。前記第１の観点から第３の変形までの
位相シフト測定方法でも０次の位相シフトを測定しうる
が、デフェーザ勾配やリードアウト勾配が０次の位相シ
フトに影響を与えてしまう可能性がある。上記第３の観
点による位相シフト測定方法では、位相軸に勾配を印加
せずにデータを収集するため、０次の位相シフト量を正
確に求めることが出来る。

【００２１】第４の観点では、本発明は、励起パルスを
送信し、反転パルスを送信し、エンコード勾配を位相軸
に印加し、リード勾配をリード軸に印加し、リワインド
勾配を位相軸に印加することを２回以上繰り返し、続い
て、反転パルスを送信し、位相軸には勾配を印加せず、
リード勾配をリード軸に印加しながらエコーから第１の
データを収集し、続いて、反転パルスを送信し、位相軸
には勾配を印加せず、リード勾配をリード軸に印加しな
がらエコーから第２のデータを収集し、前記第１のデー
タを１次元フーリエ変換して得られた位相データと前記
第２のデータを１次元フーリエ変換して得られた位相デ
ータとを基に前記エンコード勾配の影響による０次の位
相シフト量を求めることを特徴とする位相シフト測定方
法を提供する。エンコード勾配が“０”でも０次の位相
シフトは“０”でなく、オフセットの位相シフト分が存
在する。このオフセットの位相シフト分は、第１のデー
タでも第２のデータでも同方向である。一方、第１のデ
ータと第２のデータを収集する間には反転パルスがある
ので、先行するエンコード勾配等による位相シフトは、
第１のデータと第２のデータとで逆方向になる。そこ
で、第１のデータと第２のデータの差をとれば、オフセ
ットの位相シフト分を除去でき、先行するエンコード勾
配等による位相シフトを正確に測定できる。また、前記
第１の変形による位相シフト測定方法でも０次の位相シ
フトを測定しうるが、デフェーザ勾配やリードアウト勾
配が０次の位相シフトに影響を与えてしまう可能性があ
る。上記第４の観点による位相シフト測定方法では、位
相軸に勾配を印加せずにデータを収集するため、０次の
位相シフト量を正確に求めることが出来る。

【００２２】第５の観点では、本発明は、励起パルスを
送信し、第ｊの反転パルスを送信し、第ｊのエンコード
勾配を位相軸に印加し、リード勾配をリード軸に印加し
ながら第ｊのエコーからデータを収集し、第ｊのリワイ
ンド勾配を位相軸に印加することを、ｊ＝１〜Ｎについ
て繰り返し、前半のエコーから収集したデータを基に第
１のイメージを作成し、後半のエコーから収集したデー
タを基に第２のイメージを作成するスプリットエコート
レイン法において、前半のエンコード勾配およびリワイ
ンド勾配の影響による後半のエコーの１次の位相シフト
を補正するための補償用パルスを、前半の最後のエコー
の前のエンコード勾配に組み込むか、そのエンコード勾
配の直前または直後の一方または両方に付加するか、前
半の最後のエコーの後のリワインド勾配に組み込むか、
そのリワインド勾配の直前または直後の一方または両方
に付加することを特徴とする位相シフト補正方法を提供
する。上記第５の観点による位相シフト補正方法では、
前半の最後のエンコード勾配やリワインド勾配のところ
に補償用パルスを入れるため、前半のエンコード勾配等
の作るデフェーズ分が後半のエコーの間に大きな位相シ
フトを生じることを防止できる。従って、第２のイメー
ジの画質の劣化を防止することが出来る。

【００２３】第４の変形として、本発明は、上記第５の
観点の位相シフト補正方法において、上記第１の観点か
ら第３の変形による位相シフト測定方法により測定した
位相シフト量を基に、前記補償用パルスによる補正値を
決定することを特徴とする位相シフト補正方法を提供す
る。上記第４の変形による位相シフト補正方法では、上
記第１の観点から第３の変形による位相シフト測定方法
により正確に測定した位相シフト量を基に補正値を決め
るから、前半のエンコード勾配の作るデフェーズ分が後
半のエコーの間に大きな位相シフトを生じることを高精
度に防止できる。従って、第２のイメージの画質の劣化
を好適に防止することが出来る。

【００２４】第５の変形として、本発明は、上記第４の
変形の位相シフト補正方法において、前半のエンコード
勾配の極性のパターンの種類が異なるものは少なくとも
別のグループとし、各グループを代表する１以上のエン
コード勾配についての補正値は、前記第１の観点から第
３の変形による位相シフト測定方法により位相シフト量
を測定して決定し、各グループを代表しないエンコード
勾配についての補正値は、同じグループに属するエンコ
ード勾配についての前記決定した補正値を基に計算で求
めるか又は前記決定した補正値をそのまま流用すること
を特徴とする位相シフト補正方法を提供する。前半のエ
ンコード勾配の極性のパターンの種類が異なると、それ
に起因する渦電流や残留磁化の発生具合も異なり、その
影響による位相シフトも異なってくる。そこで、上記第
５の変形による位相シフト補正方法では、前半のエンコ
ード勾配の極性のパターンの種類が異なるものについて
は、それぞれ別個に上記第１の観点から第３の変形によ
る位相シフト測定方法により位相シフト量を測定する。
一方、前半のエンコード勾配の極性のパターンの種類が
同じものについては、実際の測定を行った結果を基に計
算（例えば比例計算）するか又はそのまま流用する。こ
れにより、プリスキャン時間の短縮と補正の精度を調和
させることが出来る。

【００２５】第６の観点では、本発明は、励起パルスを
送信し、第ｊの反転パルスを送信し、第ｊのエンコード
勾配を位相軸に印加し、リード勾配をリード軸に印加し
ながら第ｊのエコーからデータを収集し、第ｊのリワイ
ンド勾配を位相軸に印加することを、ｊ＝１〜Ｎについ
て繰り返し、前半のエコーから収集したデータを基に第
１のイメージを作成し、後半のエコーから収集したデー
タを基に第２のイメージを作成するスプリットエコート
レイン法において、前半のエンコード勾配およびリワイ
ンド勾配の影響による後半のエコーの０次の位相シフト
を補正するためのオフセット位相を、後半の反転パルス
に与えるか、エコーの検波位相に与えるか、両方に与え
ることを特徴とする位相シフト補正方法を提供する。上
記第６の観点による位相シフト補正方法では、後半の反
転パルスにオフセット位相を与えて送信するか、受信後
のエコーの位相検波時にオフセット位相を与えてか、両
方を行うことにより、後半のエコーに０次の位相シフト
を生じることを防止できる。従って、第２のイメージの
画質の劣化を防止することが出来る。

【００２６】第６の変形として、本発明は、前記第６の
観点の位相シフト補正方法において、前記第３の観点ま
たは第４の観点による位相シフト測定方法により測定し
た０次の位相シフト量を基に、オフセット位相を決定す
ることを特徴とする位相シフト補正方法を提供する。上
記第６の変形による位相シフト補正方法では、上記第３
の観点または第４の観点の位相シフト測定方法により正
確に測定した位相シフト量を基にオフセット位相を決め
るから、後半のエコーに０次の位相シフトを生じること
を高精度に防止できる。従って、第２のイメージの画質
の劣化を好適に防止することが出来る。

【００２７】第７の変形として、本発明は、上記第６の
変形の位相シフト補正方法において、前半のエンコード
勾配の極性のパターンの種類が異なるものは少なくとも
別のグループとし、各グループを代表する１以上のエン
コード勾配についてのオフセット位相は、前記第３の観
点または第４の観点の位相シフト測定方法により０次の
位相シフト量を測定して決定し、各グループを代表しな
いエンコード勾配についてのオフセット位相は、同じグ
ループに属するエンコード勾配についての前記決定した
オフセット位相を基に計算で求めるか又は前記決定した
オフセット位相をそのまま流用することを特徴とする位
相シフト補正方法を提供する。前半のエンコード勾配の
極性のパターンの種類が異なると、それに起因する渦電
流や残留磁化の発生具合も異なり、その影響による０次
の位相シフトも異なってくる。そこで、上記第７の変形
による位相シフト補正方法では、前半のエンコード勾配
の極性のパターンの種類が異なるものについては、それ
ぞれ別個に上記第３の観点または第４の観点の位相シフ
ト測定方法により０次の位相シフト量を測定する。一
方、前半のエンコード勾配の極性のパターンの種類が同
じものについては、実際の測定を行った結果を基に計算
（例えば比例計算）するか又はそのまま流用する。これ
により、プリスキャン時間の短縮と補正の精度を調和さ
せることが出来る。

【００２８】第７の観点では、本発明は、励起パルスを
送信した後、反転パルスを送信し、エンコード勾配を位
相軸に印加し、リード勾配をリード軸に印加しながらエ
コーからデータを収集し、リワインド勾配を位相軸に印
加することを、複数回繰り返し、収集したデータをイメ
ージを作成する高速スピンエコー法において、先行する
エンコード勾配およびリワインド勾配の影響による後続
のエコーの１次の位相シフトを補正するための補償用パ
ルスを、各エンコード勾配に組み込むか、各エンコード
勾配の直前または直後の一方または両方に付加するか、
各リワインド勾配に組み込むか、各リワインド勾配の直
前または直後の一方または両方に付加することを特徴と
する位相シフト補正方法を提供する。上記第７の観点に
よる位相シフト補正方法では、高速スピンエコー法にお
ける複数のエンコード勾配やリワインド勾配のそれぞれ
のところに補償用パルスを入れるため、先行するエンコ
ード勾配等の作るデフェーズ分が後続のエコーに位相シ
フトを生じることを防止できる。従って、イメージの画
質の劣化を防止することが出来る。

【００２９】第８の変形として、本発明は、前記第７の
観点の位相シフト補正方法において、前記第１の観点ま
たは第１の変形の位相シフト測定方法により測定した位
相シフト量を基に、前記補償用パルスによる補正値を決
定することを特徴とする位相シフト補正方法を提供す
る。上記第８の変形による位相シフト補正方法では、上
記第１の観点または第１の変形の位相シフト測定方法に
より正確に測定した位相シフト量を基に補正値を決める
から、位相シフトを高精度に防止できる。従って、イメ
ージの画質の劣化を好適に防止することが出来る。

【００３０】第９の変形として、本発明は、上記第８の
変形の位相シフト補正方法において、エンコード勾配の
極性のパターンの種類が異なるものは少なくとも別のグ
ループとし、各グループを代表する１以上のエンコード
勾配についての補正値は、前記第１の観点または第１の
変形のいずれかに記載の位相シフト測定方法により位相
シフト量を測定して決定し、各グループを代表しないエ
ンコード勾配についての補正値は、同じグループに属す
るエンコード勾配についての前記決定した補正値を基に
計算で求めるか又は前記決定した補正値をそのまま流用
することを特徴とする位相シフト補正方法を提供する。
エンコード勾配の極性のパターンの種類が異なると、そ
れに起因する渦電流や残留磁化の発生具合も異なり、そ
の影響による位相シフトも異なってくる。そこで、上記
第９の変形による位相シフト補正方法では、エンコード
勾配の極性のパターンの種類が異なるものについては、
それぞれ別個に上記第１の観点または第１の変形による
位相シフト測定方法により位相シフト量を測定する。一
方、エンコード勾配の極性のパターンの種類が同じもの
については、実際の測定を行った結果を基に計算（例え
ば比例計算）するか又はそのまま流用する。これによ
り、プリスキャン時間の短縮と補正の精度を調和させる
ことが出来る。

【００３１】第１０の変形として、本発明は、上記第８
の観点または第９の変形の位相シフト補正方法におい
て、イメージを作成する高速スピンエコー法のパルスシ
ーケンスの２番目のリード勾配を印加せず、その代り
に、前記デフェーザ勾配を位相軸に印加し、前記リード
アウト勾配を位相軸に印加しながらエコーからデータを
収集し、そのデータにより２番目のエコーのための補償
用パルスを決定し、これを３番目以降のリード勾配につ
いても１つ１つ行って、各エコーのための補償用パルス
を決定することを特徴とする位相シフト補正方法を提供
する。上記第１０の変形による位相シフト補正方法で
は、イメージ作成用の本スキャンのパルスシーケンスと
最小限異なるプリスキャンのパルスシーケンスを用い
て、各エコーのための補償用パルスを１つ１つ決定する
ため、プリスキャン時間は長くなるが、補正の精度を高
くすることが出来る。

【００３２】第８の観点では、本発明は、ＲＦパルス送
信手段と、勾配磁場印加手段と、ＮＭＲ信号受信手段と
を備えたＭＲＩ装置において、励起パルスを送信し、反
転パルスを送信し、エンコード勾配を位相軸に印加し、
リード勾配をリード軸に印加し、リワインド勾配を位相
軸に印加することを１回以上行い、続いて、反転パルス
を送信し、デフェーザ勾配を位相軸に印加し、リードア
ウト勾配を位相軸に印加しながらＮＭＲ信号を受信して
データを収集し、そのデータを１次元フーリエ変換し、
得られた位相データを基に前記エンコード勾配の影響に
よる位相シフトを補正するための補正値を決定する位相
シフト補正値決定手段を具備し、前記決定した補正値を
基に、１次の位相シフトを補正するための補償用パルス
を、イメージ作成用のパルスシーケンスのエンコード勾
配に組み込むか、エンコード勾配の直前または直後の一
方または両方に付加するか、リワインド勾配に組み込む
か、リワインド勾配の直前または直後の一方または両方
に付加することを特徴とするＭＲＩ装置を提供する。上
記第８の観点のＭＲＩ装置では、前記第１の観点から第
３の変形による位相シフト測定方法を好適に実施でき
る。また、前記第５の観点から第５の変形または第７の
観点から第１０の変形による位相シフト補正方法を好適
に実施できる。従って、位相シフトによるアーチファク
トのないイメージを得ることが出来る。

【００３３】第９の観点では、本発明は、ＲＦパルス送
信手段と、勾配磁場印加手段と、ＮＭＲ信号受信手段と
を備えたＭＲＩ装置において、励起パルスを送信し、反
転パルスを送信し、エンコード勾配を位相軸に印加し、
リード勾配をリード軸に印加し、リワインド勾配を位相
軸に印加することを１回以上行い、続いて、反転パルス
を送信し、位相軸には勾配を印加せず、リード勾配をリ
ード軸に印加しながらエコーからデータを収集し、その
データを１次元フーリエ変換して得られた位相データを
基に前記エンコード勾配の影響による０次の位相シフト
を補正するための補正値を決定する位相シフト補正値決
定手段を具備し、前記決定した補正値を基に、０次の位
相シフトを補正するためのオフセット位相を、反転パル
スに与えるか、エコーの検波位相に与えるか、両方に与
えることを特徴とするＭＲＩ装置を提供する。上記第９
の観点のＭＲＩ装置では、前記第３の観点または第４の
観点による位相シフト測定方法を好適に実施できる。ま
た、前記第６の観点から第７の変形による位相シフト補
正方法を好適に実施できる。従って、位相シフトによる
アーチファクトのないイメージを得ることが出来る。

【００３４】第１０の観点では、本発明は、ＲＦパルス
送信手段と、勾配磁場印加手段と、ＮＭＲ信号受信手段
とを備えたＭＲＩ装置において、励起パルスを送信し、
第１の反転パルスを送信し、位相軸には勾配を印加せ
ず、リード勾配をリード軸に印加しながら第１のエコー
からデータを収集し、そのデータを１次元フーリエ変換
して得られた位相データを基に第１の０次項φ０を求め
ると共に、励起パルスを送信し、第１の反転パルスを送
信し、エンコード勾配を位相軸に印加し、リード勾配を
リード軸に印加し、リワインド勾配を位相軸に印加し、
続いて、第２の反転パルスを送信し、位相軸には勾配を
印加せず、リード勾配をリード軸に印加しながら第２の
エコーからデータを収集し、そのデータを１次元フーリ
エ変換して得られた位相データを基に第２の０次項φ１
を求め、前記第１の０次項φ０と前記第２の０次項φ１
を基に前記エンコード勾配の影響による０次の位相シフ
トを補正するための補正値を決定する位相シフト補正値
決定手段を具備し、前記決定した補正値を基に、０次の
位相シフトを補正するためのオフセット位相を、第２の
反転パルス以降の反転パルスに与えるか、エコーの検波
位相に与えるか、両方に与えることを特徴とするＭＲＩ
装置を提供する。なお、第１のエコーからデータを収集
する代りに第２のエコーからデータを収集し、そのデー
タを１次元フーリエ変換して得られた位相データを基に
前記第１の０次項φ０を求めても等価であり、本発明の
範囲に含まれる。デフェーザ勾配とリフェーザ勾配を加
えると、それらによっても渦電流が生じるので、その影
響を除去しないと本スキャンの場合と異なる測定結果が
出かねない。そこで、上記第１０の観点のＭＲＩ装置で
は、デフェーザ勾配とリフェーザ勾配を全く用いないよ
うにした。これにより、デフェーザ勾配とリフェーザ勾
配による渦電流の影響を除去でき、先行するエンコード
勾配等による位相シフトを正確に測定できる。従って、
位相シフト補正を好適に実施でき、位相シフトによるア
ーチファクトのないイメージを得ることが出来る。

【００３５】

【発明の実施の形態】以下、図に示す実施形態により本
発明をさらに詳しく説明する。なお、これにより本発明
が限定されるものではない。

【００３６】図１は、本発明の第１の実施形態〜第９の
実施形態に共通するＭＲＩ装置のブロック図である。こ
のＭＲＩ装置１００において、マグネットアセンブリ１
は、内部に被検体を挿入するための空間部分（孔）を有
し、この空間部分を取りまくようにして、被検体に一定
の主磁場を印加する永久磁石１ｐと、スライス軸，位相
軸，リード軸の勾配磁場を発生するための勾配磁場コイ
ル１ｇと、被検体内の原子核のスピンを励起するための
ＲＦパルスを与える送信コイル１ｔと、被検体からのＮ
ＭＲ信号を検出する受信コイル１ｒとが配置されてい
る。前記勾配磁場コイル１ｇ，送信コイル１ｔおよび受
信コイル１ｒは、それぞれ勾配磁場駆動回路３，ＲＦ電
力増幅器４および前置増幅器５に接続されている。

【００３７】シーケンス記憶回路８は、計算機７からの
指令に従い、記憶しているパルスシーケンスに基づいて
勾配磁場駆動回路３を操作し、前記マグネットアセンブ
リ１の勾配磁場コイル１ｇから勾配磁場を発生させると
共に、ゲート変調回路９を操作し、ＲＦ発振回路１０の
搬送波出力信号を所定タイミング・所定包絡線形状のパ
ルス状信号に変調し、それをＲＦパルスとしてＲＦ電力
増幅器４に加え、ＲＦ電力増幅器４でパワー増幅した
後、前記マグネットアセンブリ１の送信コイル１ｔに印
加し、目的のスライス領域を選択励起する。

【００３８】前置増幅器５は、マグネットアセンブリ１
の受信コイル１ｒで検出された被検体からのＮＭＲ信号
を増幅し、位相検波器１２に入力する。位相検波器１２
は、ＲＦ発振回路１０の搬送波出力信号を参照信号と
し、前置増幅器５からのＮＭＲ信号を位相検波して、Ａ
／Ｄ変換器１１に与える。Ａ／Ｄ変換器１１は、位相検
波後のアナログ信号をディジタル信号に変換して、計算
機７に入力する。計算機７は、Ａ／Ｄ変換器１１からデ
ータを読み込み、画像再構成演算を行い、目的のスライ
ス領域のイメージを生成する。このイメージは、表示装
置６にて表示される。また、計算機７は、操作卓１３か
ら入力された情報を受け取るなどの全体的な制御を受け
持つ。

【００３９】−第１の実施形態− 第１の実施形態では、スプリットエコートレイン法にお
ける１次および０次の位相シフトを補正する。スプリッ
トエコートレイン法では、第１のイメージを生成するた
めに印加する前半のエンコード勾配ｇｙの作るデフェー
ズ分が、第２のイメージを生成するための後半のエコー
収集時間（リード勾配ｇｘｗを印加する時間）内に加わ
り、後半のエコーecho’の間に大きな位相シフトを生じ
る。この位相シフトは、Ｋ−空間の位相軸に対して１次
の変化をする位相シフト量を持つので、１次の位相シフ
トと呼ぶ。また、ＲＦシールドやコイル配置の不整合な
どにより勾配磁場が０次項（Ｂ０成分）を含む。その勾
配磁場の０次項は、場所依存性のない均一成分ではある
が、結果的に渦電流のような指数関数で落ちる時間特性
を有し、位相シフトの原因となる。この位相シフトを０
次の位相シフトと呼ぶ。

【００４０】図２および図３は、スプリットエコートレ
イン法における１次および０次の位相シフトを補正する
ために加える補正値を決定するための位相シフト補正値
決定処理を示すフロー図である。

【００４１】ステップＳ１では、図４，図６，図７，図
８に示すプリスキャンのパルスシーケンスにおける補償
用パルスｇyp1iの振幅ａ_gyp1i の初期値を適当な値にセ
ットする。ここで、ｉは繰返し番号である。どのような
ｉを選ぶかについては後述する。

【００４２】ステップＳ２では、エンコード勾配ｇｙの
極性が正の場合は図４に示すプリスキャンのパルスシー
ケンスにより、エンコード勾配ｇｙの極性が負の場合は
図７に示すプリスキャンのパルスシーケンスにより、後
半のエコーecho'1，echo'2からそれぞれデータを収集す
る。図４，図７に示すプリスキャンのパルスシーケンス
では、まず、励起パルスＲとスライス勾配ｓｓを印加す
る。次に、前半の第１の反転パルスＰ１とスライス勾配
ｓｓを印加する。次に、エンコード勾配gy1iを位相軸に
印加する。次に、リード勾配ｇｘｗを印加する。その
後、前記エンコード勾配gy1iと時間積分値が等しく逆極
性のリワインド勾配ｇyr1iを位相軸に印加する。図４は
エコートレイン数Ｎ＝４の場合であるが、Ｎ＞４の場合
には、同様に、反転パルスＰの送信からリワインド勾配
ｇyrの印加までを（Ｎ／２−１）回反復する。次に、前
半の最後の反転パルスＰ２とスライス勾配ｓｓを印加
し、エンコード勾配gy2iを位相軸に印加する。次に、前
記リード勾配ｇｘｗの前半分のリード勾配ｇｘrを印加
し、リード勾配をリフェーズ（rephase）する。この後
は、リード勾配は“０”とする。次いで、前記エンコー
ド勾配gy2iと時間積分値が等しく逆極性のリワインド勾
配ｇyr2iを位相軸に印加する。なお、破線で示すよう
に、リード勾配ｇｘｒの代りに通常と同じリード勾配ｇ
ｘｗを印加し、その後で、リード勾配ｇｘｗの半分の面
積のリード軸リワインド勾配ｇｘｒ’を印加しても等価
である。

【００４３】次に、エンコード勾配ｇｙと逆極性の補償
用パルスｇyp1iを位相軸に印加する。続いて、スティミ
ュレーテッドエコーを取り除くためのクラッシャ勾配ｇ
zi，ｇzｒiを前後に付けて、後半の第１の反転パルス
Ｐ’１とスライス勾配ｓｓを印加する。次に、前記リワ
インド勾配ｇyr2iと等しいデフェーザ勾配gywd1iを位相
軸に印加する。次に、リードアウト勾配ｇyw1iを位相軸
に印加しながらエコーecho’１からＮＭＲ信号を受信
し、その後、前記デフェーザ勾配gywd1iと等しいリフェ
ーザ勾配gywr1iを位相軸に印加する。図５に、図４のパ
ルスシーケンスの前半の最後のエンコード勾配gy2iから
リードアウト勾配ｇyw1iまでのＫ−空間ＫSPでの軌跡を
示す。エコーecho’１からデータを収集する時には、Ｋ
ｙ軸上のＮＭＲ信号成分のみを受信することになる。こ
の意味から、前記デフェーザ勾配gywd1i，リードアウト
勾配ｇyw1iおよびリフェーザ勾配gywr1iを位相軸プロジ
ェクション（Y projection）パルスＹpg1iという。

【００４４】次に、後半の第２の反転パルスＰ’２とス
ライス勾配ｓｓを印加する。次に、前記位相軸プロジェ
クションパルスＹpg1iと等しい位相軸プロジェクション
パルスＹpg2iを位相軸に印加しながらエコーecho'2から
ＮＭＲ信号を受信する。そして、最後に、横磁化をスポ
イルするためのキラー勾配Ｋilを印加する。

【００４５】図４，図７におけるτi1，τi2は、リード
アウト勾配ｇywの中心とエコーecho’の中心の時間ずれ
である。この時間ずれτi1，τi2は、エンコード勾配gy
1i，gy2iに起因する渦電流および残留磁化の影響を表し
ている。フーリエ変換の性質から、収集したデータを１
次元フーリエ変換した結果の位相の１次項として、時間
ずれτi1，τi2を正確に知ることが出来る。すなわち、
エンコード勾配gy1i，gy2iに起因する渦電流および残留
磁化の影響を正確に知ることが出来る。

【００４６】図３に戻り、ステップＳ３では、エコーec
ho’１，echo’２から得られた各データを１次元フーリ
エ変換し、その結果の位相データを最小二乗法等を用い
てフィッティング（curve fitting） し、その結果の各
１次項ｄi1，ｄi2と０次項Ｂ０i1，Ｂ０i2を求める。ス
テップＳ４では、次式により１次の位相シフト量φi+を
求める。 φi+＝｛（ｄi1−ｄi2）／２｝Ｘ_res・１０⁶／（γ・ｆov） …(１) ここで、（ｄi1−ｄi2）／２は、エンコード勾配ｇｙが
“０”でも生じるオフセットの位相シフト分を打ち消す
操作である。また、Ｘ_res は、エコーecho’のサンプリ
ングポイント数である。また、γは、磁気回転比であ
る。また、ｆovは、撮像視野（field of view）の大き
さ（ｃｍ）である。

【００４７】ステップＳ５では、次式により０次の位相
シフト量ΔＢ０i+を求める。 ΔＢ０i+＝（Ｂ０i1−Ｂ０i2）／２ これは、エンコード勾配ｇｙが“０”でも生じるオフセ
ットの位相シフト分を打ち消す操作である。

【００４８】図３に進み、ステップＳ６では、エンコー
ド勾配ｇｙの極性が正の場合は図６に示すプリスキャン
のパルスシーケンスにより、エンコード勾配ｇｙの極性
が負の場合は図８に示すプリスキャンのパルスシーケン
スにより、後半のエコーecho'1，echo'2からそれぞれデ
ータを収集する。図６，図８に示すプリスキャンのパル
スシーケンスは、図４，図７に示すプリスキャンのパル
スシーケンスの位相軸プロジェクションパルスＹpgの極
性を逆にしたものである。

【００４９】図３に戻り、ステップＳ７では、エコーec
ho’１，echo’２から得られた各データを１次元フーリ
エ変換し、その結果の位相データを最小二乗法等を用い
てフィッティングし、その結果の各１次項ｄi1，ｄi2と
０次項Ｂ０i1，Ｂ０i2を求める。ステップＳ８では、次
式により１次の位相シフト量φi-を求める。 φi-＝｛（ｄi1−ｄi2）／２｝Ｘ_res・１０⁶／（γ・ｆ
ov） ここで、（ｄi1−ｄi2）／２は、エンコード勾配ｇｙが
“０”でも生じるオフセットの位相シフト分を打ち消す
操作である。また、Ｘ_res は、エコーecho’のサンプリ
ングポイント数である。また、γは、磁気回転比であ
る。また、ｆovは、撮像視野の大きさである。

【００５０】ステップＳ９では、次式により０次の位相
シフト量ΔＢ０i-を求める。 ΔＢ０i-＝（Ｂ０i1−Ｂ０i2）／２ これは、エンコード勾配ｇｙが“０”でも生じるオフセ
ットの位相シフト分を打ち消す操作である。

【００５１】ステップＳ１０では、次式により１次の位
相シフト量φｉを求める。 φｉ＝｛（φi+）−（φi-）｝／２ これは、位相軸プロジェクションパルスＹpgにより発生
する渦電流の影響を除去する操作である。極性が逆の２
つの位相軸プロジェクションパルスＹpgによる各渦電流
は、１次の位相シフト量については逆方向に作用するの
で、１次の位相シフト量φi+，φi-の差を１／２するこ
とにより影響を除去できる。

【００５２】ステップＳ１１では、次式により０次の位
相シフト量ΔＢ０ｉを求める。 ΔＢ０ｉ＝｛（ΔＢ０i+）−（ΔＢ０i-）｝／２ これは、位相軸プロジェクションパルスＹpgにより発生
する渦電流の影響を除去する操作である。

【００５３】ステップＳ１２では、補償用パルスｇyp1i
の振幅ａ_gyp1i を次式により更新する。 _{new_}ａ_gyp1i ＝（１＋φｉ／gyp1iarea）_{old_}ａ_gyp1i …(２) ここで、_{new_}ａ_gyp1iは更新後の振幅であり、_{old_}ａ
_gyp1iは更新前の振幅であり、gyp1iarea は更新前の補
償用パルスｇyp1iの面積（デフェーズ量に相当する）で
ある。

【００５４】ステップＳ１３では、上記ステップＳ２〜
Ｓ１２を設定回数繰り返す。その後、ステップＳ１４へ
進む。ステップＳ１４では、イメージ作成用の本スキャ
ンにおける全ての繰り返し番号ｉまたは適当なｉについ
て上記ステップＳ１〜Ｓ１３を繰り返す。その後、処理
を終了する。

【００５５】イメージ作成用の本スキャンにおける全て
の繰り返し番号ｉについて上記ステップＳ１〜Ｓ１３を
繰り返した場合は、多数のエンコード勾配について最適
の補償用パルスｇyp1iを決定できる。但し、プリスキャ
ン時間が長くなる。

【００５６】一方、適当なｉについて上記ステップＳ１
〜Ｓ１３を繰り返した場合は、プリスキャン時間を短く
出来る。但し、上記ステップＳ１〜Ｓ１３を実行してい
ないｉに対応するエンコード勾配についての補償用パル
スｇyp1iは、決定した補償用パルスｇyp1iから補間計算
で求めるか、決定した補償用パルスｇyp1iの最近傍のも
のを流用する必要がある。適当なｉの選択方法は、パル
スシーケンスの前半のエンコード勾配ｇｙの極性のパタ
ーンの種類が同じものは１つのグループにしても，２以
上のグループに分けてもよいが、パルスシーケンスの前
半のエンコード勾配ｇｙの極性のパターンの種類が異な
るものは少なくとも別のグループとし、各グループのう
ちの平均的なエンコード勾配に対応するｉを選ぶように
する。

【００５７】図９，図１０は、Ｋ−空間ＫSPのセグメン
ト分割方法の説明図である。図９は、セントリックビュ
ーオーダリング（Ｃentric Ｖiew Ｏrdering）と呼ばれ
る方法で、エンコード勾配ｇy1，ｇy2の極性を正にし
て、エコーecho１，echo’１によりセグメントＳgp１の
データを収集し、エコーecho２，echo’２によりセグメ
ントＳgp２のデータを収集する。また、エンコード勾配
ｇy1，ｇy2の極性を負にして、エコーecho１，echo’１
によりセグメントＳgn１のデータを収集し、エコーecho
２，echo’２によりセグメントＳgn２のデータを収集す
る。この場合、エンコード勾配の極性のパターンは２種
類である。図１０は、シーケンシャルビューオーダリン
グ（Ｓequential Ｖiew Ｏrdering）と呼ばれる方法
で、エンコード勾配ｇy1の極性を正にし，ｇy２の極性
を負にして、エコーecho１，echo’１によりセグメント
Ｓgp１のデータを収集し、エコーecho２，echo’２によ
りセグメントＳgn２のデータを収集する。また、エンコ
ード勾配ｇy1の極性を負，ｇy2の極性を正にして、エコ
ーecho１，echo’１によりセグメントＳgn１のデータを
収集し、エコーecho２，echo’２によりセグメントＳgp
２のデータを収集する。この場合、エンコード勾配の極
性のパターンは２種類である。

【００５８】図９の場合（Ｎ＝４でセントラルビューオ
ーダリング方法の場合）、エンコード勾配ｇｙの極性が
正のパターンと、エンコード勾配ｇｙの極性が負のパタ
ーンとを少なくとも別のグループにし、少なくとも各グ
ループを代表するｉについて補正値を決定する必要があ
る。実際には、同一パターンでも２以上のグループに分
けてグループ数を増やし、各グループを代表するｉにつ
いて補正値を決定するのが好ましい。図１０の場合（Ｎ
＝４でシーケンシャルビューオーダリング方法の場
合）、エンコード勾配ｇy1の極性が正，ｇy2の極性が負
のパターンと、エンコード勾配ｇy1の極性が負，ｇy2の
極性が正のパターンとを少なくとも別のグループにし、
各グループを代表するｉについて補正値を決定する必要
がある。実際には、同一パターンでも２以上のグループ
に分けてグループ数を増やし、各グループを代表するｉ
について補正値を決定するのが好ましい。

【００５９】図１１は、第１のイメージおよび第２のイ
メージを得るための本スキャンのパルスシーケンスの第
１例である。このパルスシーケンスでは、前記位相シフ
ト補正値決定処理により決定した補償用パルスｇyp1i
を、スプリットエコートレイン法のパルスシーケンスの
前半の最後のリワインド勾配ｇyr2iの後に付加して、１
次の位相シフトを補正する。また、前記位相シフト補正
値決定処理で求めた０次の位相シフト量ΔＢ０ｉに相当
する量だけオフセット位相を与えた反転パルスＰ”を送
信して、０次の位相シフトを補正する。図１２は、第１
のイメージおよび第２のイメージを得るための本スキャ
ンのパルスシーケンスの第２例である。このパルスシー
ケンスでは、前記位相シフト補正値決定処理により決定
した補償用パルスｇyp1iを、スプリットエコートレイン
法のパルスシーケンスの前半の最後のリワインド勾配ｇ
yr2iの前に、リード勾配ｇｘｗとオーバラップしないよ
うに、付加して、１次の位相シフトを補正する。また、
前記位相シフト補正値決定処理で求めた０次の位相シフ
ト量ΔＢ０ｉに相当する量だけオフセット位相を与えた
反転パルスＰ”を送信して、０次の位相シフトを補正す
る。なお、図１１と図１２の組み合わせとして、補償用
パルスｇyp1iを、リワインド勾配ｇyr2iの前後に分けて
付加してもよい。

【００６０】図１３は、第１のイメージおよび第２のイ
メージを得るための本スキャンのパルスシーケンスの第
３例である。このパルスシーケンスでは、前記位相シフ
ト補正値決定処理により決定した補償用パルスｇyp1i
を、スプリットエコートレイン法のパルスシーケンスの
前半の最後のリワインド勾配ｇyr2iに組み込んで、１次
の位相シフトを補正する。また、前記位相シフト補正値
決定処理で求めた０次の位相シフト量ΔＢ０ｉに相当す
る量だけオフセット位相を与えた反転パルスＰ”を送信
して、０次の位相シフトを補正する。図１４は、第１の
イメージおよび第２のイメージを得るための本スキャン
のパルスシーケンスの第４例である。このパルスシーケ
ンスでは、前記位相シフト補正値決定処理により決定し
た補償用パルスｇyp1iを、スプリットエコートレイン法
のパルスシーケンスの前半の最後のエンコード勾配gy2i
に組み込んで、１次の位相シフトを補正する。また、前
記位相シフト補正値決定処理で求めた０次の位相シフト
量ΔＢ０ｉに相当する量だけオフセット位相を与えた反
転パルスＰ”を送信して、０次の位相シフトを補正す
る。

【００６１】なお、ΔＢ０ｉに相当する量だけオフセッ
ト位相を与えた反転パルスＰ”を送信する代りに、位相
検波器１２での検波位相にΔＢ０ｉに相当する量だけオ
フセット位相を与えて、０次の位相シフトを補正しても
よい。あるいは、オフセット位相を与えた反転パルスを
送信するのに加えて、位相検波器１２での検波位相にオ
フセット位相を与えて、ΔＢ０ｉに相当する量だけ０次
の位相シフトを補正してもよい。

【００６２】図１５は、１次の位相シフトの補正の効果
の説明図である。図９のセグメントＳgp１，Ｓgp２，Ｓ
gn１，Ｓgn２の中央に相当するｉを選んで決定した各報
償用パルスｇpy1iを各セグメントの１次の位相シフトの
補正に用いたため、各セグメントの中央では１次の位相
シフト量は“０”であるが、各セグメントの中央から離
れると位相シフトが残っている。

【００６３】図１６，図１７は、０次の位相シフトの補
正の効果の説明図である。図１６は、オフセット位相を
与えた反転パルスを送信することのみにより補正した場
合である。図１７は、オフセット位相を与えた反転パル
スを送信するのに加えて、位相検波器１２での検波位相
にオフセット位相を与えた場合である。０次の位相シフ
トは、ビューによる差が小さければよく、“０”でなく
ても構わない。

【００６４】以上により、前半のエンコード勾配ｇｙに
起因する渦電流または残留磁化の影響による第２のイメ
ージの画質の劣化を防止することが出来る。

【００６５】−第２の実施形態− 上記第１の実施形態では、エンコード勾配ｇｙが“０”
でも生じるオフセットの位相シフト分を打ち消すため、
２つの位相軸プロジェクションパルスＹpg1i，Ｙpg2iを
用いてそれぞれ得た位相シフト量の差をとったが、１つ
の位相軸プロジェクションパルスＹpg1iのみを用いても
よい。この場合でも、位相軸プロジェクションパルスＹ
pg1iの極性を逆転してそれぞれ得た位相シフト量の差を
とる時にオフセットの位相シフト分が打ち消される。

【００６６】なお、エンコード勾配ｇｙを“０”にした
時の１次項ｄi0を測定すれば、それがオフセットの位相
シフト分であるから、それを引くようにしても、オフセ
ットの位相シフト分を打ち消すことが出来る。

【００６７】−第３の実施形態− 第３の実施形態では、位相軸プロジェクションパルスＹ
pgを用いないで０次の位相シフト補正値を決定する。図
４，図６〜図８のプリスキャンのパルスシーケンスにお
ける位相軸プロジェクションパルスＹpg（特に、デフェ
ーザ勾配gywdおよびリワインド勾配gywr）は、それ自身
が０次の位相シフトの要因となりうる。そこで、０次の
位相シフト補正値の決定に際しては、位相軸プロジェク
ションパルスＹpgを用いないことが好ましい。

【００６８】図１８，図１９は、０次位相シフト補正値
決定処理のフロー図である。ステップＱ１では、図２０
に示すプリスキャンのパルスシーケンスにより後半のエ
コーecho'1，echo'2から、それぞれデータを収集する。
図２０に示すプリスキャンのパルスシーケンスは、従来
のスプリットエコートレイン法のパルスシーケンスで、
エンコード勾配ｇｙの極性が正のものにおいて、パルス
シーケンスの後半は、エンコード勾配ｇｙ’を“０”に
し、オフセット位相を与えた反転パルスＰ”を送信する
ようにしたものである。

【００６９】ステップＱ２では、エコーecho’１，ech
o’２から得られた各データを１次元フーリエ変換し、
その結果の位相データを最小二乗法等を用いてフィッテ
ィングし、その結果の各０次項Ｂ０i1，Ｂ０i2を求め
る。ステップＱ３では、次式により０次の位相シフト量
ΔＢ０i+を求める。 ΔＢ０i+＝（Ｂ０i1−Ｂ０i2）／２ これは、エンコード勾配ｇｙが“０”でも生じるオフセ
ットの位相シフト分を打ち消す操作である。ステップＱ
４では、後半の反転パルスＰ”のオフセット位相をΔＢ
０i+にセットする。ステップＱ５では、上記ステップＱ
１〜Ｑ４を設定回数繰り返す。その後、図１９のステッ
プＱ６へ進む。図１９のステップＱ６では、図２１に示
すプリスキャンのパルスシーケンスにより後半のエコー
echo'1，echo'2から、それぞれデータを収集する。図２
１に示すプリスキャンのパルスシーケンスは、従来のス
プリットエコートレイン法のパルスシーケンスで、エン
コード勾配ｇｙの極性が負のものにおいて、パルスシー
ケンスの後半は、エンコード勾配ｇｙ’を“０”にし、
オフセット位相を与えた反転パルスＰ”を送信するよう
にしたものである。

【００７０】ステップＱ７では、エコーecho’１，ech
o’２から得られた各データを１次元フーリエ変換し、
その結果の位相データを最小二乗法等を用いてフィッテ
ィングし、その結果の各０次項Ｂ０i1，Ｂ０i2を求め
る。ステップＱ８では、次式により０次の位相シフト量
ΔＢ０i-を求める。 ΔＢ０i-＝（Ｂ０i1−Ｂ０i2）／２ これは、エンコード勾配ｇｙが“０”でも生じるオフセ
ットの位相シフト分を打ち消す操作である。ステップＱ
９では、後半の反転パルスＰ”のオフセット位相をΔＢ
０i-にセットする。ステップＱ１０では、上記ステップ
Ｑ６〜Ｑ９を設定回数繰り返す。その後、ステップＱ１
１へ進む。

【００７１】ステップＱ１１では、全てのｉまたは適当
なｉについて上記ステップＱ１〜ＱＳ１０を繰り返す。
その後、ステップＱ１２へ進む。

【００７２】ステップＱ１２では、センタービューの前
後のビューにおける０次の位相シフト量ΔＢ０i+とΔＢ
０i-の差ΔＢ００を求める。 ΔＢ００＝（ΔＢ０i+）−（ΔＢ０i-） ステップＱ１３では、位相検波器１２の検波位相にオフ
セット位相ΔＢ００を与える。そして、処理を終了す
る。

【００７３】以上により決定した０次の位相シフトの補
正値により補正を行えば、０次の位相シフトに起因する
第２のイメージの画質の劣化を防止することが出来る。

【００７４】−第４の実施形態− 第４の実施形態では、高速ＳＥ法における１次および０
次の位相シフトを補正する。高速ＳＥ法では、先に結像
したエコーに対して印加したエンコード勾配ｇｙの作る
デフェーズ分が、次に観測するエコーに加わり、次に観
測するエコーに位相シフトを生じる。この位相シフト
は、場所の１次関数になるので、１次の位相シフトと呼
ぶ。また、ＲＦシールドやコイル配置の不整合などによ
り勾配磁場が０次項（Ｂ０成分）を含む。その勾配磁場
の０次項は、場所依存性のない均一成分ではあるが、結
果的に渦電流のような指数関数で落ちる時間特性を有
し、位相シフトの原因となる。この位相シフトを０次の
位相シフトと呼ぶ。

【００７５】図２２および図２３は、高速ＳＥ法におけ
る１次および０次の位相シフトを補正するために加える
補正値を決定するための位相シフト補正値決定処理を示
すフロー図である。

【００７６】ステップＦ１では、図２４〜図２７に示す
プリスキャンのパルスシーケンスにおける補償用パルス
gyp(K-1)iの振幅ａ_gyp(K-1)iの初期値を適当な値にセッ
トする。ここで、Ｋは、位相シフトを補正するエコーの
番号であり、Ｋ＝２，…，Ｎ（エコートレイン数）であ
る。ｉは繰返し番号である。どのようなｉを選ぶかにつ
いては後述する。

【００７７】ステップＦ２では、高速ＳＥ法のパルスシ
ーケンスの第Ｋエコーのエンコード勾配gyKiおよびリワ
インド勾配ｇyrKiの代りに位相軸プロジェクションパル
スＹgpKiを挿入したプリスキャンのパルスシーケンスに
よりエコーechoＫからデータを収集する。図２４に、Ｋ
＝２，Ｎ＝４のときのプリスキャンのパルスシーケンス
を例示する。このプリスキャンのパルスシーケンスで
は、まず、励起パルスＲとスライス勾配ｓｓを印加す
る。次に、第１の反転パルスＰ１とスライス勾配ｓｓを
印加する。次に、エンコード勾配gy1iを位相軸に印加す
る。次に、通常のリード勾配ｇｘｗの前半分のリード勾
配ｇｘrを印加し、リード勾配をリフェーズ（rephase）
する。この後は、エコーecho３からＮＭＲ信号を受信す
る時までリード勾配を“０”とする。次に、前記エンコ
ード勾配gy1iと時間積分値が等しく逆極性のリワインド
勾配ｇyr1iを位相軸に印加する。なお、破線で示すよう
に、リード勾配ｇｘｒの代りに通常と同じリード勾配ｇ
ｘｗを印加し、その後で、リード勾配ｇｘｗの半分の面
積のリード軸リワインド勾配ｇｘｒ’を印加しても等価
である。

【００７８】次に、エンコード勾配gy1iと逆極性の補償
用パルスｇyp1iを位相軸に印加する。次に、スティミュ
レーテッドエコーを取り除くためのクラッシャ勾配ｇz
i，ｇzｒiを前後に付けて、第２の反転パルスＰ２とス
ライス勾配ｓｓを印加する。

【００７９】次に、前記リワインド勾配ｇyr1iと等しい
デフェーザ勾配gywd2iを位相軸に印加する。次に、リー
ドアウト勾配ｇyw2iを位相軸に印加しながらエコーecho
２からＮＭＲ信号を受信し、その後、前記デフェーザ勾
配gywd2iと等しいリフェーザ勾配gywr2iを位相軸に印加
する。前記デフェーザ勾配gywd2i，リードアウト勾配ｇ
yw2iおよびリフェーザ勾配gywr2iを位相軸プロジェクシ
ョンパルスＹpg2iという。

【００８０】次に、第３の反転パルスＰ３とスライス勾
配ｓｓを印加し、エンコード勾配gy3iを位相軸に印加
し、通常のリード勾配ｇｘｗの後半分のリード勾配ｇｘ
ｄを印加する。これは、次のエコーecho４を結像させる
ためである。なお、破線で示すように、リード勾配ｇｘ
ｗの半分の面積のリード軸デフェーズ勾配ｇｘｄ’を印
加しておき、通常と同じリード勾配ｇｘｗを印加しても
等価である。その後、前記エンコード勾配gy3iと時間積
分値が等しく逆極性のリワインド勾配ｇyr3iを位相軸に
印加する。次いで、第４の反転パルスＰ４とスライス勾
配ｓｓを印加し、エンコード勾配gy4iを位相軸に印加
し、リード勾配ｇｘｗを印加し、その後、前記エンコー
ド勾配gy4iと時間積分値が等しく逆極性のリワインド勾
配ｇyr4iを位相軸に印加する。そして、最後に、横磁化
をスポイルするためのキラー勾配Ｋilを印加する。

【００８１】図２４におけるτｉは、リードアウト勾配
ｇyw2iの中心とエコーecho２の中心の時間ずれである。
この時間ずれτｉは、エンコード勾配gy1iに起因する渦
電流および残留磁化の影響を表している。フーリエ変換
の性質から、収集したデータを１次元フーリエ変換した
結果の位相の１次項として、時間ずれτｉを正確に知る
ことが出来る。すなわち、エンコード勾配gy1iに起因す
る渦電流および残留磁化の影響を正確に知ることが出来
る。

【００８２】図２２に戻り、ステップＦ３では、エコー
echoＫから得られたデータを１次元フーリエ変換し、そ
の結果の位相データを最小二乗法等を用いてフィッティ
ングし、その結果の１次項ｄi1と０次項Ｂ０i+を求め
る。ステップＦ４では、次式により１次の位相シフト量
φi+を求める。 φi+＝ｄi1・Ｘ_res・１０⁶／（γ・ｆov） ここで、Ｘ_res は、エコーechoＫのサンプリングポイン
ト数である。また、γは、磁気回転比である。また、ｆ
ovは、撮像視野の大きさ（ｃｍ）である。

【００８３】ステップＦ５では、前記ステップＦ２で用
いたプリスキャンのパルスシーケンスの位相軸プロジェ
クションパルスの極性を反転したプリスキャンのパルス
シーケンスによりエコーechoＫからデータを収集する。
図２５に、図２４のパルスシーケンスの位相軸プロジェ
クションパルスＹgp2iの極性を反転したパルスシーケン
スを示す。

【００８４】図２２に戻り、ステップＦ６では、エコー
echoＫから得られたデータを１次元フーリエ変換し、そ
の結果の位相データを最小二乗法等を用いてフィッティ
ングし、その結果の１次項ｄi1と０次項Ｂ０i-を求め
る。ステップＦ７では、次式により１次の位相シフト量
φi-を求める。 φi-＝ｄi1・Ｘ_res・１０⁶／（γ・ｆov） ここで、Ｘ_res は、エコーecho’のサンプリングポイン
ト数である。また、γは、磁気回転比である。また、ｆ
ovは、撮像視野の大きさである。

【００８５】図２３へ進み、ステップＦ８では、次式に
より１次の位相シフト量φｉを求める。 φｉ＝｛（φi+）−（φi-）｝／２ これは、位相軸プロジェクションパルスＹpgにより発生
する渦電流の影響を除去する操作である。極性が逆の２
つの位相軸プロジェクションパルスＹpgによる各渦電流
は、１次の位相シフト量については逆方向に作用するの
で、１次の位相シフト量φi+，φi-の差を１／２するこ
とにより影響を除去できる。

【００８６】ステップＦ９では、次式により０次の位相
シフト量ΔＢ０ｉを求める。 ΔＢ０ｉ＝｛（Ｂ０i+）−（Ｂ０i-）｝／２ これは、位相軸プロジェクションパルスＹpgにより発生
する渦電流の影響を除去する操作である。

【００８７】ステップＦ１０では、補償用パルスｇyp(K
-1)iの振幅ａ_gyp(K-1)i を次式により更新する。_{new_} ａ_gyp(K-1)i ＝（１＋φｉ／gyp(K-1)iarea）_{old_}
ａ_gyp(K-1)i ここで、_{new_}ａ_gyp(K-1)iは更新後の振幅であり、_{old_}
ａ_gyp(K-1)iは更新前の振幅であり、gyp(K-1)iarea は
更新前の補償用パルスｇyp(K-1)iの面積（デフェーズ量
に相当する）である。

【００８８】ステップＦ１１では、上記ステップＦ２〜
Ｆ１０を設定回数繰り返す。その後、ステップＦ１２へ
進む。

【００８９】ステップＦ１２では、全てのＫについて上
記ステップＦ２〜Ｆ１１を繰り返す。図２６に、Ｋ＝
３，Ｎ＝４のときのプリスキャンのパルスシーケンスを
例示する。また、図２７に、Ｋ＝４，Ｎ＝４のときのプ
リスキャンのパルスシーケンスを例示する。Ｋ＝Ｎまで
上記ステップＦ２〜Ｆ１１を繰り返したら、ステップＦ
１３へ進む。

【００９０】ステップＦ１３では、イメージ作成用の本
スキャンにおける全ての繰り返し番号ｉまたは適当なｉ
について上記ステップＦ１〜Ｆ１２を繰り返す。その
後、処理を終了する。

【００９１】イメージ作成用の本スキャンにおける全て
の繰り返し番号ｉについて上記ステップＦ１〜Ｆ１２を
繰り返した場合は、多数のエンコード勾配について最適
の補償用パルスｇyp1iを決定できる。但し、プリスキャ
ン時間が長くなる。

【００９２】一方、適当なｉについて上記ステップＦ１
〜Ｆ１２を繰り返した場合は、プリスキャン時間を短く
出来る。但し、上記ステップＦ１〜Ｆ１２を実行してい
ないｉに対応するエンコード勾配についての補償用パル
スｇyp1iは、決定した補償用パルスｇyp1iから補間計算
で求めるか、決定した補償用パルスｇyp1iの最近傍のも
のを流用する必要がある。適当なｉの選択方法は、パル
スシーケンスのエンコード勾配ｇy1〜ｇyKの極性のパタ
ーンの種類が同じものは１つのグループにしても，２以
上のグループに分けてもよいが、パルスシーケンスのエ
ンコード勾配ｇy1〜ｇyKの極性のパターンの種類が異な
るものは少なくとも別のグループとし、各グループのう
ちの平均的なエンコード勾配に対応するｉを選ぶように
する。

【００９３】図２８，図２９は、Ｋ−空間ＫSPのセグメ
ント分割方法の説明図である。図２８のセグメント分割
方法では、エンコード勾配ｇy1の極性を正，エンコード
勾配ｇy2の極性を正，エンコード勾配ｇy3の極性を負，
エンコード勾配ｇy4の極性を負にして、エコーecho１，
echo２，echo３，echo４により、セグメントＳｇ１，Ｓ
ｇ２，Ｓｇ３，Ｓｇ４のデータをそれぞれ収集する。こ
の場合、エンコード勾配の極性のパターンは１種類であ
る。図２９のセグメント分割方法では、エンコード勾配
ｇy1の極性を正，エンコード勾配ｇy2の極性を負，エン
コード勾配ｇy3の極性を負，エンコード勾配ｇy4の極性
を正にして、エコーecho１，echo２，echo３，echo４に
より、セグメントＳｇ11，Ｓｇ21，Ｓｇ31，Ｓｇ41のデ
ータをそれぞれ収集する。また、エンコード勾配ｇy1の
極性を正，エンコード勾配ｇy2の極性を正，エンコード
勾配ｇy3の極性を負，エンコード勾配ｇy4の極性を負に
して、エコーecho１，echo２，echo３，echo４により、
セグメントＳｇ12，Ｓｇ22，Ｓｇ32，Ｓｇ42のデータを
それぞれ収集する。この場合、エンコード勾配の極性の
パターンは２種類である。

【００９４】図２８の場合、エンコード勾配ｇy1の極性
が正，エンコード勾配ｇy2の極性が正，エンコード勾配
ｇy3の極性が負，エンコード勾配ｇy4の極性が負のパタ
ーンについて補正値を決定すればよい。但し、実際に
は、２以上のグループに分けてグループ数を増やし、各
グループを代表するｉについて補正値を決定するのが好
ましい。図２９の場合、エンコード勾配ｇy1の極性が
正，エンコード勾配ｇy2の極性が負，エンコード勾配ｇ
y3の極性が負，エンコード勾配ｇy4の極性が正のパター
ンと、エンコード勾配ｇy1の極性が正，エンコード勾配
ｇy2の極性が正，エンコード勾配ｇy3の極性が負，エン
コード勾配ｇy4の極性が負のパターンとを少なくとも別
のグループにし、各グループを代表するｉについて補正
値を決定する必要がある。実際には、同一パターンでも
２以上のグループに分けてグループ数を増やし、各グル
ープを代表するｉについて補正値を決定するのが好まし
い。

【００９５】図３０は、イメージ作成用の本スキャンの
パルスシーケンスの例である。このパルスシーケンスで
は、前記位相シフト補正値決定処理により決定した補償
用パルスｇyp1i，ｇyp2i，ｇyp3iを、高速ＳＥ法のパル
スシーケンスのリワインド勾配ｇyr1i，ｇyr2i，ｇyr3i
の後に付加して、１次の位相シフトを補正する。また、
前記位相シフト補正値決定処理で求めた０次の位相シフ
ト量ΔＢ０ｉに相当する量だけオフセット位相を与えた
反転パルスＰ”を送信して、０次の位相シフトを補正す
る。

【００９６】なお、前記第１の実施形態で説明したよう
に、補償用パルスｇyp1i，ｇyp2i，ｇyp3iを、高速ＳＥ
法のパルスシーケンスのリワインド勾配ｇyr1i，ｇyr2
i，ｇyr3iの前に付加してもよい。また、補償用パルス
ｇyp1i，ｇyp2i，ｇyp3iを、高速ＳＥ法のパルスシーケ
ンスのリワインド勾配ｇyr1i，ｇyr2i，ｇyr3iの前後に
分けて付加してもよい。さらに、補償用パルスｇyp1i，
ｇyp2i，ｇyp3iを、リワインド勾配ｇyr1i，ｇyr2i，ｇ
yr3iに組み込んでもよい。また、補償用パルスｇyp1i，
ｇyp2i，ｇyp3iを、エンコード勾配gy1i，gy2i，gy3iに
組み込んでもよい。

【００９７】なお、ΔＢ０ｉに相当する量だけオフセッ
ト位相を与えた反転パルスＰ”を送信する代りに、位相
検波器１２での検波位相にΔＢ０ｉに相当する量だけオ
フセット位相を与えて、０次の位相シフトを補正しても
よい。あるいは、オフセット位相を与えた反転パルスを
送信するのに加えて、位相検波器１２での検波位相にオ
フセット位相を与えて、ΔＢ０ｉに相当する量だけ０次
の位相シフトを補正してもよい。

【００９８】以上により、エンコード勾配ｇｙに起因す
る渦電流または残留磁化の影響によるイメージの画質の
劣化を防止することが出来る。

【００９９】−第５の実施形態− 前記第４の実施形態では、位相軸プロジェクションパル
スＹpg(K-1)iの極性を逆転してそれぞれ得た位相シフト
量の差をとることで、エンコード勾配ｇｙが“０”でも
生じるオフセットの位相シフト分を打ち消しているが、
エンコード勾配ｇｙを“０”にした時の１次項ｄi0を測
定すれば、それがオフセットの位相シフト分であるか
ら、それを引くようにしても、オフセットの位相シフト
分を打ち消すことが出来る。

【０１００】−第６の実施形態− 前記実施形態では、イメージ作成用の本スキャンのパル
スシーケンスに準じたプリスキャンのパルスシーケンス
を用いたため、補正量を高精度に決定できたが、スキャ
ン時間が長くかかった。第６の実施形態では、精度を下
げる代りに、スキャン時間を短縮できるようにする。

【０１０１】図３１は、第６の実施形態にかかる位相シ
フト補正値決定処理を示すフロー図である。ステップＥ
１では、図３２に示すプリスキャンのパルスシーケンス
における補償用パルスgypnの振幅ａ_gypnの初期値を適当
な値にセットする。ここで、ｎは、補正値を決定したい
エンコード勾配の番号である。

【０１０２】ステップＥ２では、図３２のパルスシーケ
ンスによりエコーecho２からデータを収集する。図３２
のプリスキャンのパルスシーケンスでは、まず、励起パ
ルスＲとスライス勾配ｓｓを印加する。次に、第１の反
転パルスＰ１とスライス勾配ｓｓを印加する。次に、補
正値を決定したいエンコード勾配ｇｙｎを位相軸に印加
する。次に、通常のリード勾配ｇｘｗの前半分のリード
勾配ｇｘｒを印加する。この後は、リード勾配を“０”
とする。次に、前記エンコード勾配ｇynと時間積分値が
等しく逆極性のリワインド勾配gyrnに補償用パルスgypn
を組み込んだパルスを位相軸に印加する。次に、スティ
ミュレーテッドエコーを取り除くためのクラッシャ勾配
ｇｚ，ｇｚｒを前後に付けて第２の反転パルスＰ２とス
ライス勾配ｓｓを印加する。次に、前記リワインド勾配
gyrnと等しいデフェーザ勾配ｇywdnを位相軸に印加す
る。次に、リードアウト勾配gywnを位相軸に印加しなが
らエコーecho２からＮＭＲ信号を受信し、その後、前記
デフェーザ勾配ｇywdnと等しいリフェーザ勾配ｇywrnを
位相軸に印加する。前記デフェーザ勾配ｇywdn，リード
アウト勾配gywnおよびリフェーザ勾配ｇywrnを位相軸プ
ロジェクションパルスYpgnという。図３２におけるτｉ
は、リードアウト勾配gywnの中心とエコーecho２の中心
の時間ずれである。この時間ずれτｉは、エンコード勾
配ｇynに起因する渦電流および残留磁化の影響を表して
いる。フーリエ変換の性質から、収集したデータを１次
元フーリエ変換した結果の位相の１次項として、時間ず
れτｉを正確に知ることが出来る。すなわち、エンコー
ド勾配ｇynに起因する渦電流および残留磁化の影響を正
確に知ることが出来る。

【０１０３】図３１に戻り、ステップＥ３では、エコー
echo２から得られたデータを１次元フーリエ変換し、そ
の結果の位相データを最小二乗法等を用いてフィッティ
ングし、その結果の１次項ｄｎと０次項Ｂ０ｎを求め
る。ステップＥ４では、次式により１次の位相シフト量
φｎを求める。 φｎ＝ｄｎ・Ｘ_res・１０⁶／（γ・ｆov） ここで、Ｘ_res は、エコーecho２のサンプリングポイン
ト数である。また、γは、磁気回転比である。また、ｆ
ovは、撮像視野の大きさ（ｃｍ）である。

【０１０４】ステップＥ５では、補償用パルスgypnの振
幅ａ_gypnを更新する。_{new_} ａ_gypn＝（１＋φｉ／gypnarea）_{old_}ａ_gypn ここで、_{new_}ａ_gypnは更新後の振幅であり、_{old_}ａ_gypn
は更新前の振幅であり、gypnareaは更新前の補償用パル
スgypnの面積（デフェーズ量に相当する）である。

【０１０５】ステップＥ６では、上記ステップＥ２〜Ｅ
５を設定回数繰り返す。その後、ステップＥ７へ進む。
ステップＥ７では、必要なｎについて上記ステップＥ１
〜Ｅ７を繰り返す。その後、処理を終了する。以上によ
り、必要なエンコード勾配ｇynについて補正値を決定で
きる。

【０１０６】−第７の実施形態− 前記実施形態では、デフェーザ勾配ｇywdnおよびリフェ
ーザ勾配ｇywrnを用いたため、これらによる渦電流の影
響があり、本スキャンの場合と異なる測定結果が出かね
ない。そこで、第７の実施形態では、デフェーザ勾配ｇ
ywdnおよびリフェーザ勾配ｇywrnを全く用いないで補正
値を決定できるようにする。

【０１０７】図３３は、第７の実施形態にかかる位相シ
フト補正値決定処理を示すフロー図である。ステップＱ
５１では、図３４に示すプリスキャンのパルスシーケン
スにより、第１のエコーecho１からデータを収集する。
図３４のプリスキャンのパルスシーケンスでは、まず、
励起パルスＲとスライス勾配ｓｓを印加する。次に、第
１の反転パルスＰ１とスライス勾配ｓｓを印加する。次
に、リード勾配ｇｘｗを印加しながら第１のエコーecho
１からデータを収集する。位相軸には、勾配を印加しな
い。図３３に戻り、ステップＱ５２では、前記ステップ
Ｑ５１で収集したデータを一次元フーリエ変換し、得ら
れた位相データをフィッティングし、勾配磁場が“０”
の場所の位相を第１の０次項φ０とする。

【０１０８】ステップＱ５３では、図３５に示すプリス
キャンのパルスシーケンスにより、第２のエコーecho２
からデータを収集する。図３５のプリスキャンのパルス
シーケンスでは、まず、励起パルスＲとスライス勾配ｓ
ｓを印加する。次に、第１の反転パルスＰ１とスライス
勾配ｓｓを印加する。次に、補正値を決定したいエンコ
ード勾配ｇynを位相軸に印加する。次に、リード勾配ｇ
ｘｗを印加する。次に、前記エンコード勾配ｇynと時間
積分値が等しく逆極性のリワインド勾配gyrnを位相軸に
印加する。次に、第２の反転パルスＰ２とスライス勾配
ｓｓを印加する。次に、リード勾配ｇｘｗを印加しなが
ら第２のエコーecho２からデータを収集する。図３３に
戻り、ステップＱ５４では、前記ステップＱ５３で収集
したデータを一次元フーリエ変換し、得られた位相デー
タをフィッティングし、勾配磁場が“０”の場所の位相
を第２の０次項φ１とする。

【０１０９】ステップＱ５５では、第１の０次項φ０と
第２の０次項φ１から補正値Φを求める。すなわち、Φ
＝（φ０−φ１）／２とする。図３６に、上記演算の原
理を示す。図３４のプリスキャンのパルスシーケンスで
第１の反転パルスＰ１を印加した後の磁化Ｍxy＿P1が、
何らかの原因で回転軸ｘに一致せず、回転軸ｘから反時
計回りにφ０だけ回転した位置になったと仮定する。さ
らに、図３５のプリスキャンのパルスシーケンスで位相
軸に印加したエンコード勾配ｇynとリワインド勾配gyrn
の渦電流の影響により、回転軸ｘが回転軸ｘ’へと移動
したとする。このとき、前記φ０は、前記第１の０次項
として観測される。また、図３５のプリスキャンのパル
スシーケンスで第１の反転パルスＰ１を印加した後の磁
化もＭxy＿P1であるから、第２の反転パルスＰ２を印加
した後の磁化Ｍxy＿P2は、回転軸ｘ’の回りに前記磁化
Ｍxy＿P1を回転させた位置にある。そして、前記第２の
０次項φ１は、回転軸ｘから時計回りに磁化Ｍxy＿P2を
見た位相を意味することになる。そこで、（φ０−φ
１）は、磁化Ｍxy＿P2から磁化Ｍxy＿P1を時計回りに見
た位相になる。そして、図３６から判るように、（φ０
−φ１）／２は、回転軸ｘ’から磁化Ｍxy＿P1を時計回
りに見た位相Φになる。よって、この位相Φを補正値と
すれば、エンコード勾配ｇynとリワインド勾配gyrnの渦
電流の影響により移動した回転軸ｘ’が磁化Ｍxy＿P1の
位置まで移動したことと等価となる。すると、回転軸が
磁化Ｍxy＿P1に一致するから、反転パルスＰ２を加えた
後の磁化Ｍxy＿P2は磁化Ｍxy＿P1に一致することとな
り、エンコード勾配ｇynとリワインド勾配gyrnの渦電流
の影響が打ち消される。

【０１１０】図３３に戻り、ステップＱ５６では、前記
補正量Φを図３４，図３５のプリスキャンパルスシーケ
ンスの第２の反転パルスＰ２以降の反転パルスにオフセ
ット位相として加える。ステップＱ５７では、上記ステ
ップＱ５１〜Ｑ５６を設定回数繰り返し、最適の補正値
Φに集束させる。その後、処理を終了する。以上によ
り、必要なエンコード勾配ｇynについて補正値を決定で
きる。

【０１１１】

【発明の効果】本発明の位相シフト測定方法によれば、
エンコード勾配等に起因する渦電流や残留磁化の影響に
よる１次および０次の位相シフトを正確に測定できる。
また、本発明の位相シフト補正方法によれば、エンコー
ド勾配等に起因する渦電流や残留磁化の影響による１次
および０次の位相シフトを高精度に補正でき、画質の劣
化のないイメージを得ることが出来る。また、本発明の
ＭＲＩ装置によれば、上記位相シフト測定方法および位
相シフト補正方法を好適に実施できる。また、エンコー
ド勾配等に起因する渦電流または残留磁化の影響による
画質の劣化のないイメージを得ることが出来る。

【図面の簡単な説明】

【図１】本発明の実施形態にかかるＭＲＩ装置を示すブ
ロック図である。

【図２】第１の実施形態にかかる位相シフト補正値決定
処理を示すフロー図である。

【図３】図２の続きのフロー図である。

【図４】第１の実施形態にかかるプリスキャンのパルス
シーケンスの例示図である。

【図５】図４のパルスシーケンスにおけるＫ−空間上の
軌跡の説明図である。

【図６】図４のパルスシーケンスの位相軸プロジェクシ
ョンパルスを反転したパルスシーケンスの例示図であ
る。

【図７】第１の実施形態にかかるプリスキャンのパルス
シーケンスの別の例示図である。

【図８】図７のパルスシーケンスの位相軸プロジェクシ
ョンパルスを反転したパルスシーケンスの例示図であ
る。

【図９】セントラルビューオーダリングのセグメントと
エコーの関係説明図である。

【図１０】シーケンシャルビューオーダリングのセグメ
ントとエコーの関係説明図である。

【図１１】補償用パルスを後半の最後のリワインド勾配
の後に付加したスプリットエコートレイン法のパルスシ
ーケンスの例示図である。

【図１２】補償用パルスを後半の最後のリワインド勾配
の前に付加したスプリットエコートレイン法のパルスシ
ーケンスの例示図である。

【図１３】補償用パルスを後半の最後のリワインド勾配
に組み込んだスプリットエコートレイン法のパルスシー
ケンスの例示図である。

【図１４】補償用パルスを後半の最後のエンコード勾配
に組み込んだスプリットエコートレイン法のパルスシー
ケンスの例示図である。

【図１５】補償用パルスによる補正効果の説明図であ
る。

【図１６】反転パルスにオフセット位相を与えることに
よる補正効果の説明図である。

【図１７】反転パルスおよび検波位相にオフセット位相
を与えることによる補正効果の説明図である。

【図１８】第３の実施形態にかかる０次位相シフト補正
値決定処理を示すフロー図である。

【図１９】図１８の続きのフロー図である。

【図２０】第３の実施形態にかかるプリスキャンのパル
スシーケンスの例示図である。

【図２１】第３の実施形態にかかるプリスキャンのパル
スシーケンスの別の例示図である。

【図２２】第４の実施形態にかかる位相シフト補正値決
定処理を示すフロー図である。

【図２３】図２２の続きのフロー図である。

【図２４】第４の実施形態にかかるプリスキャンのパル
スシーケンスの例示図である。

【図２５】図２４のパルスシーケンスの位相軸プロジェ
クションパルスを反転したパルスシーケンスの例示図で
ある。

【図２６】第４の実施形態にかかるプリスキャンのパル
スシーケンスの別の例示図である。

【図２７】第４の実施形態にかかるプリスキャンのパル
スシーケンスのさらに別の例示図である。

【図２８】セグメントとエコーの関係説明図である。

【図２９】セグメントとエコーの別の関係説明図であ
る。

【図３０】補償用パルスを各リワインド勾配の後に付加
した高速ＳＥ法のパルスシーケンスの例示図である。

【図３１】第６の実施形態にかかる位相シフト補正値決
定処理を示すフロー図である。

【図３２】第６の実施形態にかかるプリスキャンのパル
スシーケンスの例示図である。

【図３３】第７の実施形態にかかる位相シフト補正値決
定処理を示すフロー図である。

【図３４】第６の実施形態にかかるプリスキャンのパル
スシーケンスの例示図である。

【図３５】第６の実施形態にかかる別のプリスキャンの
パルスシーケンスの例示図である。

【図３６】第６の実施形態にかかる補正値の演算の原理
を示す説明図である。

【図３７】従来のスプリットエコートレイン法のパルス
シーケンスの例示図である。

【図３８】スプリットエコートレイン法におけるＫ−空
間でのデータ収集軌跡の説明図である。

【図３９】従来の高速ＳＥ法のパルスシーケンスの例示
図である。

【図４０】高速ＳＥ法におけるＫ−空間でのデータ収集
軌跡の説明図である。

【符号の説明】

１００ ＭＲＩ装置 １ マグネットアセンブリ １ｇ 勾配磁場コイル １ｔ 送信コイル １ｐ 永久磁石 ７ 計算機 ８ シーケンス記憶回路 １２ 位相検波器

Claims (10)

【特許請求の範囲】
1. 【請求項１】 励起パルスを送信し、反転パルスを送信
   し、エンコード勾配を位相軸に印加し、リード勾配をリ
   ード軸に印加し、リワインド勾配を位相軸に印加するこ
   とを１回以上行い、続いて、反転パルスを送信し、デフ
   ェーザ勾配を位相軸に印加し、リードアウト勾配を位相
   軸に印加しながらエコーからデータを収集し、そのデー
   タを１次元フーリエ変換して得られた位相データを基に
   前記エンコード勾配の影響による位相シフト量を求める
   ことを特徴とする位相シフト測定方法。
2. 【請求項２】 励起パルスを送信し、反転パルスを送信
   し、エンコード勾配を位相軸に印加し、リード勾配をリ
   ード軸に印加し、リワインド勾配を位相軸に印加するこ
   とを２回以上繰り返し、続いて、反転パルスを送信し、
   デフェーザ勾配を位相軸に印加し、リードアウト勾配を
   位相軸に印加しながらエコーから第１のデータを収集
   し、リフェーザ勾配を位相軸に印加し、続いて、反転パ
   ルスを送信し、デフェーザ勾配を位相軸に印加し、リー
   ドアウト勾配を位相軸に印加しながらエコーから第２の
   データを収集し、前記第１のデータを１次元フーリエ変
   換して得られた位相データと前記第２のデータを１次元
   フーリエ変換して得られた位相データとを基に前記エン
   コード勾配の影響による位相シフト量を求めることを特
   徴とする位相シフト測定方法。
3. 【請求項３】 励起パルスを送信し、反転パルスを送信
   し、エンコード勾配を位相軸に印加し、リード勾配をリ
   ード軸に印加し、リワインド勾配を位相軸に印加するこ
   とを１回以上行い、続いて、反転パルスを送信し、位相
   軸には勾配を印加せず、リード勾配をリード軸に印加し
   ながらエコーからデータを収集し、そのデータを１次元
   フーリエ変換して得られた位相データを基に前記エンコ
   ード勾配の影響による０次の位相シフト量を求めること
   を特徴とする位相シフト測定方法。
4. 【請求項４】 励起パルスを送信し、反転パルスを送信
   し、エンコード勾配を位相軸に印加し、リード勾配をリ
   ード軸に印加し、リワインド勾配を位相軸に印加するこ
   とを２回以上繰り返し、続いて、反転パルスを送信し、
   位相軸には勾配を印加せず、リード勾配をリード軸に印
   加しながらエコーから第１のデータを収集し、続いて、
   反転パルスを送信し、位相軸には勾配を印加せず、リー
   ド勾配をリード軸に印加しながらエコーから第２のデー
   タを収集し、前記第１のデータを１次元フーリエ変換し
   て得られた位相データと前記第２のデータを１次元フー
   リエ変換して得られた位相データとを基に前記エンコー
   ド勾配の影響による０次の位相シフト量を求めることを
   特徴とする位相シフト測定方法。
5. 【請求項５】 励起パルスを送信し、第ｊの反転パルス
   を送信し、第ｊのエンコード勾配を位相軸に印加し、リ
   ード勾配をリード軸に印加しながら第ｊのエコーからデ
   ータを収集し、第ｊのリワインド勾配を位相軸に印加す
   ることを、ｊ＝１〜Ｎについて繰り返し、前半のエコー
   から収集したデータを基に第１のイメージを作成し、後
   半のエコーから収集したデータを基に第２のイメージを
   作成するスプリットエコートレイン法において、前半の
   エンコード勾配およびリワインド勾配の影響による後半
   のエコーの１次の位相シフトを補正するための補償用パ
   ルスを、前半の最後のエコーの前のエンコード勾配に組
   み込むか、そのエンコード勾配の直前または直後の一方
   または両方に付加するか、前半の最後のエコーの後のリ
   ワインド勾配に組み込むか、そのリワインド勾配の直前
   または直後の一方または両方に付加することを特徴とす
   る位相シフト補正方法。
6. 【請求項６】 励起パルスを送信し、第ｊの反転パルス
   を送信し、第ｊのエンコード勾配を位相軸に印加し、リ
   ード勾配をリード軸に印加しながら第ｊのエコーからデ
   ータを収集し、第ｊのリワインド勾配を位相軸に印加す
   ることを、ｊ＝１〜Ｎについて繰り返し、前半のエコー
   から収集したデータを基に第１のイメージを作成し、後
   半のエコーから収集したデータを基に第２のイメージを
   作成するスプリットエコートレイン法において、前半の
   エンコード勾配およびリワインド勾配の影響による後半
   のエコーの０次の位相シフトを補正するためのオフセッ
   ト位相を、後半の反転パルスに与えるか、エコーの検波
   位相に与えるか、両方に与えることを特徴とする位相シ
   フト補正方法。
7. 【請求項７】 励起パルスを送信した後、反転パルスを
   送信し、エンコード勾配を位相軸に印加し、リード勾配
   をリード軸に印加しながらエコーからデータを収集し、
   リワインド勾配を位相軸に印加することを、複数回繰り
   返し、収集したデータをイメージを作成する高速スピン
   エコー法において、先行するエンコード勾配およびリワ
   インド勾配の影響による後続のエコーの１次の位相シフ
   トを補正するための補償用パルスを、各エンコード勾配
   に組み込むか、各エンコード勾配の直前または直後の一
   方または両方に付加するか、各リワインド勾配に組み込
   むか、各リワインド勾配の直前または直後の一方または
   両方に付加することを特徴とする位相シフト補正方法。
8. 【請求項８】 ＲＦパルス送信手段と、勾配磁場印加手
   段と、ＮＭＲ信号受信手段とを備えたＭＲＩ装置におい
   て、 励起パルスを送信し、反転パルスを送信し、エンコード
   勾配を位相軸に印加し、リード勾配をリード軸に印加
   し、リワインド勾配を位相軸に印加することを１回以上
   行い、続いて、反転パルスを送信し、デフェーザ勾配を
   位相軸に印加し、リードアウト勾配を位相軸に印加しな
   がらＮＭＲ信号を受信してデータを収集し、そのデータ
   を１次元フーリエ変換し、得られた位相データを基に前
   記エンコード勾配の影響による位相シフトを補正するた
   めの補正値を決定する位相シフト補正値決定手段を具備
   し、 前記決定した補正値を基に、１次の位相シフトを補正す
   るための補償用パルスを、イメージ作成用のパルスシー
   ケンスのエンコード勾配に組み込むか、エンコード勾配
   の直前または直後の一方または両方に付加するか、リワ
   インド勾配に組み込むか、リワインド勾配の直前または
   直後の一方または両方に付加することを特徴とするＭＲ
   Ｉ装置。
9. 【請求項９】 ＲＦパルス送信手段と、勾配磁場印加手
   段と、ＮＭＲ信号受信手段とを備えたＭＲＩ装置におい
   て、 励起パルスを送信し、反転パルスを送信し、エンコード
   勾配を位相軸に印加し、リード勾配をリード軸に印加
   し、リワインド勾配を位相軸に印加することを１回以上
   行い、続いて、反転パルスを送信し、位相軸には勾配を
   印加せず、リード勾配をリード軸に印加しながらエコー
   からデータを収集し、そのデータを１次元フーリエ変換
   して得られた位相データを基に前記エンコード勾配の影
   響による０次の位相シフトを補正するための補正値を決
   定する位相シフト補正値決定手段を具備し、 前記決定した補正値を基に、０次の位相シフトを補正す
   るためのオフセット位相を、反転パルスに与えるか、エ
   コーの検波位相に与えるか、両方に与えることを特徴と
   するＭＲＩ装置。
10. 【請求項１０】 ＲＦパルス送信手段と、勾配磁場印加
    手段と、ＮＭＲ信号受信手段とを備えたＭＲＩ装置にお
    いて、 励起パルスを送信し、第１の反転パルスを送信し、位相
    軸には勾配を印加せず、リード勾配をリード軸に印加し
    ながら第１のエコーからデータを収集し、そのデータを
    １次元フーリエ変換して得られた位相データを基に第１
    の０次項φ０を求めると共に、励起パルスを送信し、第
    １の反転パルスを送信し、エンコード勾配を位相軸に印
    加し、リード勾配をリード軸に印加し、リワインド勾配
    を位相軸に印加し、続いて、第２の反転パルスを送信
    し、位相軸には勾配を印加せず、リード勾配をリード軸
    に印加しながら第２のエコーからデータを収集し、その
    データを１次元フーリエ変換して得られた位相データを
    基に第２の０次項φ１を求め、前記第１の０次項φ０と
    前記第２の０次項φ１を基に前記エンコード勾配の影響
    による０次の位相シフトを補正するための補正値を決定
    する位相シフト補正値決定手段を具備し、 前記決定した補正値を基に、０次の位相シフトを補正す
    るためのオフセット位相を、第２の反転パルス以降の反
    転パルスに与えるか、エコーの検波位相に与えるか、両
    方に与えることを特徴とするＭＲＩ装置。