JP3507586B2 - 磁気共鳴イメージング装置 - Google Patents
磁気共鳴イメージング装置Info
- Publication number
- JP3507586B2 JP3507586B2 JP14352295A JP14352295A JP3507586B2 JP 3507586 B2 JP3507586 B2 JP 3507586B2 JP 14352295 A JP14352295 A JP 14352295A JP 14352295 A JP14352295 A JP 14352295A JP 3507586 B2 JP3507586 B2 JP 3507586B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- phase
- magnetic resonance
- waveform
- resonance imaging
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、スピンエコー法等によ
りエコーを収集し、このエコーに基づいて画像を生成す
る磁気共鳴イメージング装置に関する。
りエコーを収集し、このエコーに基づいて画像を生成す
る磁気共鳴イメージング装置に関する。
【０００２】
【従来の技術】磁気共鳴イメージングにおいては、パル
スシーケンスの傾斜磁場波形はハードウエアの特定、特
に渦電流の影響で歪む。このためエコー中心（ｔ＝Ｔ
Ｅ）では０次の位相がゼロにならず、このため位相画像
では位相値エラーとして現象し、また通常の絶対値画像
では形状の歪みや濃度不均一性として現象する。従来こ
の渦電流の影響を低減するために、傾斜磁場コイルに印
加する電流波形をオーバシュートやアンダーシュートさ
せて磁場波形の一次遅れ波形を理想的な矩形波に近付け
る第１の方法、Active Shielded Gradient Coil により
渦電流そのものを抑制する第２の方法、収集した磁気共
鳴信号を信号処理上で位相補正する第３の方法が用いら
れている。
スシーケンスの傾斜磁場波形はハードウエアの特定、特
に渦電流の影響で歪む。このためエコー中心（ｔ＝Ｔ
Ｅ）では０次の位相がゼロにならず、このため位相画像
では位相値エラーとして現象し、また通常の絶対値画像
では形状の歪みや濃度不均一性として現象する。従来こ
の渦電流の影響を低減するために、傾斜磁場コイルに印
加する電流波形をオーバシュートやアンダーシュートさ
せて磁場波形の一次遅れ波形を理想的な矩形波に近付け
る第１の方法、Active Shielded Gradient Coil により
渦電流そのものを抑制する第２の方法、収集した磁気共
鳴信号を信号処理上で位相補正する第３の方法が用いら
れている。
【０００３】しかし、このような方法では位相誤差を十
分抑制できず、例えばフェーズコントラスト法（Phase
Contrast法）や高速スピンエコー法（fast Spin Echo）
法等のような位相誤差に敏感なパルスシーケンスには対
応できないのが現状である。例えば、フェーズコントラ
スト法では依然、定量性の低下が見られ、また高速スピ
ンエコー法ではマルチスライスの画像間の濃度むらが見
られる。
分抑制できず、例えばフェーズコントラスト法（Phase
Contrast法）や高速スピンエコー法（fast Spin Echo）
法等のような位相誤差に敏感なパルスシーケンスには対
応できないのが現状である。例えば、フェーズコントラ
スト法では依然、定量性の低下が見られ、また高速スピ
ンエコー法ではマルチスライスの画像間の濃度むらが見
られる。
【０００４】また、第２の方法はActive Shielded Grad
ient Coil （ASGC）を随所に設置することが必要とされ
るためコストがかかるし、またASGCを設置しても十分で
ない。また、第３の方法では２次以上の成分については
単純に補正することができない。
ient Coil （ASGC）を随所に設置することが必要とされ
るためコストがかかるし、またASGCを設置しても十分で
ない。また、第３の方法では２次以上の成分については
単純に補正することができない。
【０００５】
【発明が解決しようとする課題】本発明の目的は、渦電
流によるエコー中心での位相誤差を低減することのでき
る磁気共鳴診断装置を提供することである。
流によるエコー中心での位相誤差を低減することのでき
る磁気共鳴診断装置を提供することである。
【０００６】
【課題を解決するための手段】本発明は、傾斜磁場コイ
ルに所定のパルスシーケンスに従った波形の電流を供給
し、時間的に変化する磁場分布を形成し、エコーを収集
し、このエコーに基づいて画像を生成する磁気共鳴イメ
ージング装置において、渦電流の時定数と前記電流の波
形とに基づいて空間的な磁場分布の時間変化を求め、こ
の時間変化に基づいてエコー中心における位相の絶対値
が極小化するように前記電流の波形を変更することを特
徴とする。
ルに所定のパルスシーケンスに従った波形の電流を供給
し、時間的に変化する磁場分布を形成し、エコーを収集
し、このエコーに基づいて画像を生成する磁気共鳴イメ
ージング装置において、渦電流の時定数と前記電流の波
形とに基づいて空間的な磁場分布の時間変化を求め、こ
の時間変化に基づいてエコー中心における位相の絶対値
が極小化するように前記電流の波形を変更することを特
徴とする。
【０００７】
【作用】本発明によれば、電流波形の変更によりエコー
中心における位相の絶対値が極小化され、渦電流による
エコー中心での位相誤差を低減することができる。ASGC
を入れた上でさらに位相誤差を小さくできる。
中心における位相の絶対値が極小化され、渦電流による
エコー中心での位相誤差を低減することができる。ASGC
を入れた上でさらに位相誤差を小さくできる。
【０００８】
【実施例】以下、図面を参照して本発明による磁気共鳴
イメージング装置の一実施例を説明する。図１に本実施
例による磁気共鳴イメージング装置の構成を示す。被検
体Ｐを収容できるように円筒状の内部空間を有するコイ
ルアセンブリ２０には、静磁場磁石１、Ｘ・Ｙ・Ｚ軸傾
斜磁場コイル２、ＲＦコイル３が装備される。常電導磁
石又は超電導磁石である静磁場磁石１は、静磁場制御装
置４から電流供給を受けて円筒内部に通常、Ｚ軸に沿っ
て静磁場を形成可能に構成されている。Ｘ・Ｙ・Ｚ軸傾
斜磁場コイル２は、Ｘ、Ｙ、Ｚ軸各々に対応する３組の
コイルから構成されていて、それぞれ傾斜磁場電源７，
８，９から電流の供給を受けて、任意に撮影断面を決め
たり、磁気共鳴信号（ここではエコー）に空間的位置情
報を与えるためのＸ、Ｙ、Ｚ各軸の傾斜磁場を形成す
る。多くの場合、、Ｘ、Ｙ、Ｚ各軸の傾斜磁場はそれぞ
れ、リード方向（周波数エンコード方向）に沿って磁場
強度が変化するリードアウト用傾斜磁場Ｇr 、位相エン
コード方向に沿って磁場強度が変化する位相エンコード
用傾斜磁場Ｇe 、スライス方向に沿って磁場強度が変化
するスライス用傾斜磁場Ｇs として用いられる。これら
３方向の磁場強度が全て線形に変化する領域（撮影領
域）内で磁気共鳴信号の収集がなされる。磁気共鳴信号
の収集時には、被検体Ｐは寝台１３の天板に載置された
状態で、天板のスライドに伴って撮像領域に挿入され
る。
イメージング装置の一実施例を説明する。図１に本実施
例による磁気共鳴イメージング装置の構成を示す。被検
体Ｐを収容できるように円筒状の内部空間を有するコイ
ルアセンブリ２０には、静磁場磁石１、Ｘ・Ｙ・Ｚ軸傾
斜磁場コイル２、ＲＦコイル３が装備される。常電導磁
石又は超電導磁石である静磁場磁石１は、静磁場制御装
置４から電流供給を受けて円筒内部に通常、Ｚ軸に沿っ
て静磁場を形成可能に構成されている。Ｘ・Ｙ・Ｚ軸傾
斜磁場コイル２は、Ｘ、Ｙ、Ｚ軸各々に対応する３組の
コイルから構成されていて、それぞれ傾斜磁場電源７，
８，９から電流の供給を受けて、任意に撮影断面を決め
たり、磁気共鳴信号（ここではエコー）に空間的位置情
報を与えるためのＸ、Ｙ、Ｚ各軸の傾斜磁場を形成す
る。多くの場合、、Ｘ、Ｙ、Ｚ各軸の傾斜磁場はそれぞ
れ、リード方向（周波数エンコード方向）に沿って磁場
強度が変化するリードアウト用傾斜磁場Ｇr 、位相エン
コード方向に沿って磁場強度が変化する位相エンコード
用傾斜磁場Ｇe 、スライス方向に沿って磁場強度が変化
するスライス用傾斜磁場Ｇs として用いられる。これら
３方向の磁場強度が全て線形に変化する領域（撮影領
域）内で磁気共鳴信号の収集がなされる。磁気共鳴信号
の収集時には、被検体Ｐは寝台１３の天板に載置された
状態で、天板のスライドに伴って撮像領域に挿入され
る。
【０００９】ＲＦコイル３は、ＲＦパルス（高周波磁場
または回転磁場ともいう）を被検体に送信し、被検体か
らの磁気共鳴信号を受信するためのコイルである。この
ように送受信にＲＦコイル３を兼用するのではなく、送
信用コイルと受信用コイルとを別体で設けてもよい。送
信器５は、対象原子核に固有のラーモア周波数に応じた
高周波パルスをＲＦコイル３に供給して、対象原子核の
スピンを励起状態にするためのものであり、ここではス
ライス選択励起法に対応して、さらに高周波パルスを周
波数調整、及びシンク関数等に基づいて振幅調整する機
能を有している。受信器６は、励起されたスピンが緩和
する過程で放出される高周波の磁気共鳴信号をＲＦコイ
ル３を介して受信し、これを増幅検波し、さらにアナロ
グ／ディジタル変換する機能を有している、シーケンサ
１０は、送信器５、受信器６、ＸＹＺ各軸の傾斜磁場電
源７，８，９を制御して所定のパルスシーケンスを実行
するものである。シーケンサ１０は、傾斜磁場電源７，
８，９それぞれに対して出力電流の波形（時間振幅変
化）を整形するための制御信号を供給する。
または回転磁場ともいう）を被検体に送信し、被検体か
らの磁気共鳴信号を受信するためのコイルである。この
ように送受信にＲＦコイル３を兼用するのではなく、送
信用コイルと受信用コイルとを別体で設けてもよい。送
信器５は、対象原子核に固有のラーモア周波数に応じた
高周波パルスをＲＦコイル３に供給して、対象原子核の
スピンを励起状態にするためのものであり、ここではス
ライス選択励起法に対応して、さらに高周波パルスを周
波数調整、及びシンク関数等に基づいて振幅調整する機
能を有している。受信器６は、励起されたスピンが緩和
する過程で放出される高周波の磁気共鳴信号をＲＦコイ
ル３を介して受信し、これを増幅検波し、さらにアナロ
グ／ディジタル変換する機能を有している、シーケンサ
１０は、送信器５、受信器６、ＸＹＺ各軸の傾斜磁場電
源７，８，９を制御して所定のパルスシーケンスを実行
するものである。シーケンサ１０は、傾斜磁場電源７，
８，９それぞれに対して出力電流の波形（時間振幅変
化）を整形するための制御信号を供給する。
【００１０】コンピュータシステム１１は、受信器６で
ディジタル化された磁気共鳴信号（磁気共鳴データ）を
取り込み、これを例えば２次元フーリエ変換することに
より磁気共鳴画像を再構成する。この画像は表示部１２
に表示される。また、コンピュータシステム１１は、詳
細は後述するが、予め計算により又は実測により求めら
れた装置固有の渦電流の時定数を保持していて、この渦
電流の時定数と、所定のパルスシーケンス上の傾斜磁場
波形に対応する電流の波形とに基づいて空間的な磁場分
布の時間変化を求め、この時間変化に基づいてエコー中
心における位相の絶対値が極小化するように上記電流波
形を変更し、この変更した電流波形を示す情報をシーケ
ンサ１０に供給する。
ディジタル化された磁気共鳴信号（磁気共鳴データ）を
取り込み、これを例えば２次元フーリエ変換することに
より磁気共鳴画像を再構成する。この画像は表示部１２
に表示される。また、コンピュータシステム１１は、詳
細は後述するが、予め計算により又は実測により求めら
れた装置固有の渦電流の時定数を保持していて、この渦
電流の時定数と、所定のパルスシーケンス上の傾斜磁場
波形に対応する電流の波形とに基づいて空間的な磁場分
布の時間変化を求め、この時間変化に基づいてエコー中
心における位相の絶対値が極小化するように上記電流波
形を変更し、この変更した電流波形を示す情報をシーケ
ンサ１０に供給する。
【００１１】次に電流波形の求める方法について説明す
る。本実施例による渦電流による位相誤差低減について
説明する。図２にここで取り扱う変数ｒ，Ａを示してい
る。エコー中心（ｔ＝ＴＥ；ＴＥはエコー時間）での空
間座標（Ｚ，ｒ）における渦電流のみによる位相成分△
φ(Z,r,TE)は（１）式で与えられる。
る。本実施例による渦電流による位相誤差低減について
説明する。図２にここで取り扱う変数ｒ，Ａを示してい
る。エコー中心（ｔ＝ＴＥ；ＴＥはエコー時間）での空
間座標（Ｚ，ｒ）における渦電流のみによる位相成分△
φ(Z,r,TE)は（１）式で与えられる。
【００１２】
【数１】
【００１３】渦電流の作る磁場（渦磁場）の強度の空間
的・時間的変化関数Ｈe(Z,r,t)は、ASGCの装備有無で相
違し、ASGCを装備した場合を（２）式で、ASGCを装備し
ない場合を（３）式でそれぞれ示す。
的・時間的変化関数Ｈe(Z,r,t)は、ASGCの装備有無で相
違し、ASGCを装備した場合を（２）式で、ASGCを装備し
ない場合を（３）式でそれぞれ示す。
【００１４】
【数２】
【００１５】ASGCの装備有無に応じて、（２）式又は
（３）式を、（１）式に代入し、△φ(Z,r,TE)を極小化
（最小化）するように、傾斜磁場の強度の空間的・時間
的変化関数を決定し、この関数にしたがって傾斜磁場の
波形を定める。これにより、渦電流が作る磁場に起因す
るエコー中心での位相誤差を低減することができる。
（３）式を、（１）式に代入し、△φ(Z,r,TE)を極小化
（最小化）するように、傾斜磁場の強度の空間的・時間
的変化関数を決定し、この関数にしたがって傾斜磁場の
波形を定める。これにより、渦電流が作る磁場に起因す
るエコー中心での位相誤差を低減することができる。
【００１６】例えば、スピンエコー法において、スライ
ス方向の傾斜磁場（スライス選択用傾斜磁場）の変化に
よる渦電流に関して見ると、図３に示すように、位相補
償用の傾斜磁場Ｇsa、Ｇsbの面積比（面積＝振幅×印加
時間）を、△φ(Z,r,TE)を極小化（最小化）するように
調整すればよい。勿論この調整はコンピュータシステム
１１でなされる。
ス方向の傾斜磁場（スライス選択用傾斜磁場）の変化に
よる渦電流に関して見ると、図３に示すように、位相補
償用の傾斜磁場Ｇsa、Ｇsbの面積比（面積＝振幅×印加
時間）を、△φ(Z,r,TE)を極小化（最小化）するように
調整すればよい。勿論この調整はコンピュータシステム
１１でなされる。
【００１７】次に渦磁場の強度の空間的・時間的変化関
数Ｈe(Z,r,t)について説明する。ここではスライス方向
の傾斜磁場の変化による渦電流が作る磁場（渦磁場）に
起因する位相誤差を補償する場合を考える。スライス方
向の傾斜磁場の変化による渦磁場は、スライス平面内で
２次分布を示すと実測された。
数Ｈe(Z,r,t)について説明する。ここではスライス方向
の傾斜磁場の変化による渦電流が作る磁場（渦磁場）に
起因する位相誤差を補償する場合を考える。スライス方
向の傾斜磁場の変化による渦磁場は、スライス平面内で
２次分布を示すと実測された。
【００１８】ここで、傾斜磁場波形、渦電流の時定数、
位相の空間的分布を関数式にしたがってシュミレートす
る。説明の便宜上、空間的位置をｒのみで扱うものとす
る。渦磁場Ｈe(r,t)は（４）式で与えられる。
位相の空間的分布を関数式にしたがってシュミレートす
る。説明の便宜上、空間的位置をｒのみで扱うものとす
る。渦磁場Ｈe(r,t)は（４）式で与えられる。
【００１９】
【数３】
【００２０】この場合、渦磁場Ｈe(r,t)、渦磁場Ｈe(r,
t)と渦磁場ゼロのときの原磁場Ｈo(r,t)とを合わせた観
測磁場Ｈm(r,t)を求め、エコー中心（ｔ＝ＴＥ）におけ
る位相φ（r,TE）を（５）式にしたがって算出する。
t)と渦磁場ゼロのときの原磁場Ｈo(r,t)とを合わせた観
測磁場Ｈm(r,t)を求め、エコー中心（ｔ＝ＴＥ）におけ
る位相φ（r,TE）を（５）式にしたがって算出する。
【００２１】
【数４】
【００２２】（５）式において、渦電流の時定数τを変
えながら、位相φ（r,TE）の空間的変化を、バイポーラ
型(bipolar-type)とユニポーラ型(unipolar-type) それ
ぞれについて計算しグラフ化したものが図５（ａ）乃至
図５（ｆ）である。
えながら、位相φ（r,TE）の空間的変化を、バイポーラ
型(bipolar-type)とユニポーラ型(unipolar-type) それ
ぞれについて計算しグラフ化したものが図５（ａ）乃至
図５（ｆ）である。
【００２３】傾斜磁場波形が矩形波である場合の渦磁場
及び観測磁場の時間変化を考える。図６にこの波形の一
例を示す。傾斜磁場の変化時刻を順にＴ0 ，Ｔ1 ，Ｔ2
，…，Ｔi ，Ｔi+1 とする。各時刻の渦磁場ゼロのと
きの原磁場強度を順にＨ0 ，Ｈ1 ，Ｈ2 ，…，Ｈi とす
る。ｉは傾斜磁場の変化回数を示すパラメータである。
各時刻の渦磁場強度は（６）式、（７）式、（８）式で
与えられる。
及び観測磁場の時間変化を考える。図６にこの波形の一
例を示す。傾斜磁場の変化時刻を順にＴ0 ，Ｔ1 ，Ｔ2
，…，Ｔi ，Ｔi+1 とする。各時刻の渦磁場ゼロのと
きの原磁場強度を順にＨ0 ，Ｈ1 ，Ｈ2 ，…，Ｈi とす
る。ｉは傾斜磁場の変化回数を示すパラメータである。
各時刻の渦磁場強度は（６）式、（７）式、（８）式で
与えられる。
【００２４】
【数５】
これら（６）乃至（８）式を一般化すると、（９）式が
与えられる。
与えられる。
【００２５】
【数６】
一方観測磁場の時間変化は（１０）式で与えられる。
【００２６】
【数７】
【００２７】次に上記時刻Ｔ0 ，Ｔ1 ，Ｔ2 ，…，Ｔi
，Ｔi+1 における各期間における渦磁場の積分値を考
える。渦磁場Ｈe(t)は、ゲインＡ、時定数τの一次遅れ
系とみなせるから、（１１）式で表現できる。
，Ｔi+1 における各期間における渦磁場の積分値を考
える。渦磁場Ｈe(t)は、ゲインＡ、時定数τの一次遅れ
系とみなせるから、（１１）式で表現できる。
【００２８】
【数８】
次に原磁場Ｈo(t)が矩形波とした場合の各期間の渦磁場
の累積値は（１２）式、（１３）式、（１４）式で与え
られる。
の累積値は（１２）式、（１３）式、（１４）式で与え
られる。
【００２９】
【数９】
これら（１２）乃至（１４）式を一般化すると、（１
５）式が与えられる。
５）式が与えられる。
【００３０】
【数１０】
（１５）式の（Ｂ）項は（１６）式で与えられる。
【００３１】
【数１１】
したがって、各期間の渦磁場の累積値は（１７）式で与
えられる。
えられる。
【００３２】
【数１２】
以上をまとめると、傾斜磁場の変化回数に応じて（１
９）式、（２０）式が得られる。
９）式、（２０）式が得られる。
【００３３】
【数１３】
【００３４】次に、以上のような関数式にしたがって位
相変化をシュミレートする。なお、ここではＡ(r) ＝0.
001 ×ｒ2 、τ＝２として与えた。ただし、０≦Ａ(r)
≦１０である。図７（ａ）にユニポーラ型の位相の時間
変化を示し、図７（ｂ）にバイポーラ型の位相の時間変
化を示し、図７（ｃ）に両型の位相の空間的変化を示
す。
相変化をシュミレートする。なお、ここではＡ(r) ＝0.
001 ×ｒ2 、τ＝２として与えた。ただし、０≦Ａ(r)
≦１０である。図７（ａ）にユニポーラ型の位相の時間
変化を示し、図７（ｂ）にバイポーラ型の位相の時間変
化を示し、図７（ｃ）に両型の位相の空間的変化を示
す。
【００３５】また、τを変化させて位相の空間的変化を
関数志木にしたがって求めたものを図５（ａ）乃至図５
（ｆ）に示した。以上のシュミレートから、（１）エコ
ー中心での位相の空間的変化（ｒ vs φ(r,TE)）は、傾
斜磁場波形、渦電流の時定数に依存する、（２）ユニポ
ーラ型のスピンエコーシーケンスではτ＝２のとき、位
相φ(r,TE)は渦磁場の影響は少ない、（３）バイポーラ
型のスピンエコーシーケンスではτ＝２のとき、位相φ
(r,TE)は渦磁場の影響を受けて、ｒに依存して発散す
る、（４）バイポーラ型のスピンエコーシーケンスでは
τ＝４のとき、位相φ(r,TE)は渦磁場の影響を受けず、
ｒにほとんど依存しない、ことが検証された。
関数志木にしたがって求めたものを図５（ａ）乃至図５
（ｆ）に示した。以上のシュミレートから、（１）エコ
ー中心での位相の空間的変化（ｒ vs φ(r,TE)）は、傾
斜磁場波形、渦電流の時定数に依存する、（２）ユニポ
ーラ型のスピンエコーシーケンスではτ＝２のとき、位
相φ(r,TE)は渦磁場の影響は少ない、（３）バイポーラ
型のスピンエコーシーケンスではτ＝２のとき、位相φ
(r,TE)は渦磁場の影響を受けて、ｒに依存して発散す
る、（４）バイポーラ型のスピンエコーシーケンスでは
τ＝４のとき、位相φ(r,TE)は渦磁場の影響を受けず、
ｒにほとんど依存しない、ことが検証された。
【００３６】このように、パルスシーケンスの傾斜磁場
波形、渦電流の時定数に依存して位相は変化するが、図
８に示すように、位相がエコー中心でほぼゼロとなる例
えば検証（２）、（４）のような条件が存在する。この
条件が成立するときはスライス面内での渦磁場の空間的
分布の形状には依存しない。
波形、渦電流の時定数に依存して位相は変化するが、図
８に示すように、位相がエコー中心でほぼゼロとなる例
えば検証（２）、（４）のような条件が存在する。この
条件が成立するときはスライス面内での渦磁場の空間的
分布の形状には依存しない。
【００３７】したがって、渦電流の時定数、傾斜磁場の
強度分布のもとで、傾斜磁場の型を選択し、そして位相
補償用の傾斜磁場Ｇsa、Ｇsbの面積比（面積＝振幅×印
加時間）を、△φ(Z,r,TE)を極小化（最小化）するよう
に調整することにより電流による位相誤差を最小化する
ことができる。
強度分布のもとで、傾斜磁場の型を選択し、そして位相
補償用の傾斜磁場Ｇsa、Ｇsbの面積比（面積＝振幅×印
加時間）を、△φ(Z,r,TE)を極小化（最小化）するよう
に調整することにより電流による位相誤差を最小化する
ことができる。
【００３８】このように本実施例では傾斜磁場波形の最
適条件を計算により求めることができるのでパルスシー
ケンス調整の時間が短縮されオペレータの負担が軽減さ
れ、また渦電流による位相エコーが小さくなるので画像
の分解能が向上し濃度むらが低減され、また渦電流の空
間分布が変化しても位相エラーによる影響を抑えること
ができ、またASGCが装備されていなくても或る程度位相
エラーを低減できコスト低減に寄与し、またASGCを装備
した場合にはASGCだけの場合より位相誤差を小さくで
き、フローリフェーズのシーケンスで動きに対して高次
の位相エラーまで解消してモーションアーチファクトを
低減させことができる。
適条件を計算により求めることができるのでパルスシー
ケンス調整の時間が短縮されオペレータの負担が軽減さ
れ、また渦電流による位相エコーが小さくなるので画像
の分解能が向上し濃度むらが低減され、また渦電流の空
間分布が変化しても位相エラーによる影響を抑えること
ができ、またASGCが装備されていなくても或る程度位相
エラーを低減できコスト低減に寄与し、またASGCを装備
した場合にはASGCだけの場合より位相誤差を小さくで
き、フローリフェーズのシーケンスで動きに対して高次
の位相エラーまで解消してモーションアーチファクトを
低減させことができる。
【００３９】なお上記説明ではSE法を例に説明したが、
Fast SE 法、Fast GSE(GRASE) 法、EPI （エコープレナ
ー）法、PC（フェーズコントラスト）法、水脂肪分離法
等にも適用できる。本発明は上記実施例に限定されるこ
となく種々変形して実施可能である。
Fast SE 法、Fast GSE(GRASE) 法、EPI （エコープレナ
ー）法、PC（フェーズコントラスト）法、水脂肪分離法
等にも適用できる。本発明は上記実施例に限定されるこ
となく種々変形して実施可能である。
【００４０】
【発明の効果】本発明によれば、電流波形の変更により
エコー中心における位相の絶対値が極小化され、渦電流
によるエコー中心での位相誤差を低減することができ
る。ASGCを入れた上で、さらに位相誤差を小さくでき
る。
エコー中心における位相の絶対値が極小化され、渦電流
によるエコー中心での位相誤差を低減することができ
る。ASGCを入れた上で、さらに位相誤差を小さくでき
る。
【図１】本発明に係る磁気共鳴イメージング装置の一実
施例の構成図。
施例の構成図。
【図２】変数ｒ、Ａの説明図。
【図３】位相補償用の傾斜磁場波形を示す図。
【図４】傾斜磁場波形の型を示す図。
【図５】エコー中心における位相の型別の空間的変化を
示す図。
示す図。
【図６】矩形波の傾斜磁場波形を示す図。
【図７】位相の時間的、空間的変化を示す図。
【図８】位相の時間変化を示す図。
１…静磁場磁石、 ２…傾斜磁場コイ
ル、３…ＲＦコイル、 ４…静磁場制
御装置、５…送信器、 ６…受信
器、７，８，９…傾斜磁場電源、 １０…シーケン
サ、１１…コンピュータシステム、 １２…表示部、
１３…寝台、 ２０…コイルアセン
ブリ。
ル、３…ＲＦコイル、 ４…静磁場制
御装置、５…送信器、 ６…受信
器、７，８，９…傾斜磁場電源、 １０…シーケン
サ、１１…コンピュータシステム、 １２…表示部、
１３…寝台、 ２０…コイルアセン
ブリ。
─────────────────────────────────────────────────────
フロントページの続き
(58)調査した分野(Int.Cl.7，ＤＢ名)
A61B 5/055
Claims (4)
- 【請求項１】 傾斜磁場コイルに所定のパルスシーケン
スに従った波形の電流を供給し、時間的に変化する磁場
分布を形成し、エコーを収集し、このエコーに基づいて
画像を生成する磁気共鳴イメージング装置において、 渦電流の時定数と前記電流の波形とに基づいて空間的な
磁場分布の時間変化を求め、この時間変化に基づいてエ
コー中心における位相の絶対値が極小化するように前記
電流の波形を変更することを特徴とする磁気共鳴イメー
ジング装置。 - 【請求項２】 前記空間的な磁場分布の時間変化は２次
成分において一次遅れ応答であることを特徴とする請求
項１記載の磁気共鳴イメージング装置。 - 【請求項３】 前記傾斜磁場コイルはスライス方向軸に
関して磁場分布を形成することを特徴とする請求項１記
載の磁気共鳴イメージング装置。 - 【請求項４】 前記パルスシーケンスはスピンエコー法
と高速スピンエコー法のいずれかであることを特徴とす
る請求項１記載の磁気共鳴イメージング装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP14352295A JP3507586B2 (ja) | 1995-06-09 | 1995-06-09 | 磁気共鳴イメージング装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP14352295A JP3507586B2 (ja) | 1995-06-09 | 1995-06-09 | 磁気共鳴イメージング装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08336506A JPH08336506A (ja) | 1996-12-24 |
| JP3507586B2 true JP3507586B2 (ja) | 2004-03-15 |
Family
ID=15340703
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP14352295A Expired - Fee Related JP3507586B2 (ja) | 1995-06-09 | 1995-06-09 | 磁気共鳴イメージング装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3507586B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5864233A (en) * | 1997-02-22 | 1999-01-26 | General Electric Company | Method to reduce eddy current effects in diffusion-weighted echo planar imaging |
| US8008915B2 (en) | 2007-12-28 | 2011-08-30 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
-
1995
- 1995-06-09 JP JP14352295A patent/JP3507586B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08336506A (ja) | 1996-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3676853B2 (ja) | 表面コイル配設により核磁化分布を決定するｍｒ法 | |
| JP4229487B2 (ja) | マクスウェル項誤差を補正する方法 | |
| US8111070B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| EP1585991B1 (en) | Magnetic resonance method and device | |
| US8093894B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method for improving uniformity in sensitivity map | |
| EP1371327A1 (en) | Magnetic resonance imaging apparatus and method | |
| US8890530B2 (en) | Magnetic resonance imaging method, magnetic resonance imaging apparatus, and control device of magnetic resonance imaging apparatus | |
| JPWO2002056767A1 (ja) | 高精度コイル感度マップを用いたパラレルｍｒイメージング | |
| US7696753B2 (en) | MRI apparatus and method with B0 corrections using reference and object average units | |
| US10114097B2 (en) | Medical imaging apparatus having multiple subsystems, and operating method therefor | |
| EP3191862B1 (en) | Zero echo time mr imaging | |
| JP4343726B2 (ja) | 磁気共鳴イメージング装置及び不整磁場補正方法 | |
| WO2004095048A1 (en) | Mr imaging method | |
| JP3507586B2 (ja) | 磁気共鳴イメージング装置 | |
| US6919722B2 (en) | Image quality improvement for SENSE with low signal regions | |
| EP1367407A2 (en) | K-space correction of gradient non-linearities in MRI | |
| US20220057467A1 (en) | Epi mr imaging with distortion correction | |
| WO1998047015A1 (en) | Diffusion-weighted mri method | |
| JP6377378B2 (ja) | 磁気共鳴イメージング装置 | |
| US20040165758A1 (en) | MRI apparatus and method for adjusting MR image display parameters | |
| EP4343356A1 (en) | Mr imaging with water/fat/b0 mapping | |
| US20230400545A1 (en) | Spin echo mr imaging with spiral acquisition | |
| JPH09509A (ja) | Ｍｒ画像生成方法及びｍｒｉ装置 | |
| JP3171872B2 (ja) | 磁気共鳴診断装置 | |
| WO2022171608A1 (en) | Dixon-type water/fat separation mr imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20031216 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20031219 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071226 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081226 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091226 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091226 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101226 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101226 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111226 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121226 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131226 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |