JPH1071135A - Ｍｒ方法及び装置 - Google Patents
Ｍｒ方法及び装置Info
- Publication number
- JPH1071135A JPH1071135A JP9209338A JP20933897A JPH1071135A JP H1071135 A JPH1071135 A JP H1071135A JP 9209338 A JP9209338 A JP 9209338A JP 20933897 A JP20933897 A JP 20933897A JP H1071135 A JPH1071135 A JP H1071135A
- Authority
- JP
- Japan
- Prior art keywords
- gradient
- magnetic field
- signal
- gradient magnetic
- field
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 26
- 230000005415 magnetization Effects 0.000 claims abstract description 9
- 238000009826 distribution Methods 0.000 claims abstract description 4
- 230000004044 response Effects 0.000 claims description 13
- 230000007423 decrease Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 230000036962 time dependent Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 13
- 230000010363 phase shift Effects 0.000 abstract description 3
- 239000011159 matrix material Substances 0.000 description 9
- 230000002123 temporal effect Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 8
- 230000009466 transformation Effects 0.000 description 8
- 238000012937 correction Methods 0.000 description 5
- 230000005284 excitation Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56518—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56554—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】
【課題】 アイソセンターの外側に位置する傾けられた
スライスのＭＲ画像内でアーティファクトを除去するＭ
Ｒ方法を提供する。 【解決手段】 これらのアーティファクトの除去のため
に位相誤差の時間変動は課された（所望の）傾斜磁界の
変動及びこれらの磁界を発生する傾斜磁界コイルシステ
ムのパルス応答から計算され、測定されたＭＲ信号はそ
れにより補正される。
スライスのＭＲ画像内でアーティファクトを除去するＭ
Ｒ方法を提供する。 【解決手段】 これらのアーティファクトの除去のため
に位相誤差の時間変動は課された（所望の）傾斜磁界の
変動及びこれらの磁界を発生する傾斜磁界コイルシステ
ムのパルス応答から計算され、測定されたＭＲ信号はそ
れにより補正される。
Description
【０００１】
【発明の属する技術分野】本発明は検査領域のスライス
での核磁化が均一で定常な磁界の存在下で励起され、そ
の後にスライスに接する第一の方向に延在する傾斜を有
する少なくとも一つの第一の傾斜磁界が発生され、スラ
イスに接する第二の方向に延在する傾斜を有する可能な
第二の傾斜磁界が発生され、ＭＲ信号は検査領域から受
信され、少なくとも２つの傾斜磁界コイルシステムの傾
斜磁界は第一と第二の傾斜磁界を発生するために組み合
わされる少なくとも一つのシーケンスを含むＭＲ方法に
関する。
での核磁化が均一で定常な磁界の存在下で励起され、そ
の後にスライスに接する第一の方向に延在する傾斜を有
する少なくとも一つの第一の傾斜磁界が発生され、スラ
イスに接する第二の方向に延在する傾斜を有する可能な
第二の傾斜磁界が発生され、ＭＲ信号は検査領域から受
信され、少なくとも２つの傾斜磁界コイルシステムの傾
斜磁界は第一と第二の傾斜磁界を発生するために組み合
わされる少なくとも一つのシーケンスを含むＭＲ方法に
関する。
【０００２】
【従来の技術】この種のＭＲ方法はＺｈｏｕ等及びＧａ
ｔｅｈｏｕｓｅ等によるによる文献ＩＳＭＲＭ Ａｂｓ
ｔｒａｃｔ， ３８６，１４７７又は１４８１頁から知
られている。これらの文献はＥＰＩ方法を記載し、それ
は位相エンコーディング傾斜磁界及び読み取り傾斜磁界
の両方がｘ方向とｙ方向に関してそれぞれオブリークに
延在し、それによりｘ方向とｙ方向に対する傾斜磁界コ
イルシステムの傾斜磁界がこれらの傾斜磁界を発生する
よう組み合わされなければならない。上記の文献からＥ
ＰＩのようなシーケンスの位相エンコーディング傾斜磁
界は読み取り傾斜磁界により影響され、それによりいわ
ゆるゴースト画像がそのような影響が補正されない場合
に発生する。そのような影響はそれぞれの傾斜磁界コイ
ルシステムを有するチャンネルの時間的な振る舞いが異
なる場合にのみ発生する。
ｔｅｈｏｕｓｅ等によるによる文献ＩＳＭＲＭ Ａｂｓ
ｔｒａｃｔ， ３８６，１４７７又は１４８１頁から知
られている。これらの文献はＥＰＩ方法を記載し、それ
は位相エンコーディング傾斜磁界及び読み取り傾斜磁界
の両方がｘ方向とｙ方向に関してそれぞれオブリークに
延在し、それによりｘ方向とｙ方向に対する傾斜磁界コ
イルシステムの傾斜磁界がこれらの傾斜磁界を発生する
よう組み合わされなければならない。上記の文献からＥ
ＰＩのようなシーケンスの位相エンコーディング傾斜磁
界は読み取り傾斜磁界により影響され、それによりいわ
ゆるゴースト画像がそのような影響が補正されない場合
に発生する。そのような影響はそれぞれの傾斜磁界コイ
ルシステムを有するチャンネルの時間的な振る舞いが異
なる場合にのみ発生する。
【０００３】本発明の元となった研究により、提案され
た補正段階は限定された効果のみを有するものであるこ
とが判明した。垂直又は第一又は第二の傾斜磁界の方向
のいずれもがそれぞれの傾斜磁界コイルシステムの一つ
により発生された傾斜磁界の傾斜が向けられる一の方向
と一致しないスライスの検査中にゴースト画像はＭＲ信
号の位相が傾斜磁界により影響される故になお発生す
る。これらの位相誤差はＥＰＩのような位相感応性ＭＲ
方法の場合には画像アーティファクト妨害を発生させ
る。
た補正段階は限定された効果のみを有するものであるこ
とが判明した。垂直又は第一又は第二の傾斜磁界の方向
のいずれもがそれぞれの傾斜磁界コイルシステムの一つ
により発生された傾斜磁界の傾斜が向けられる一の方向
と一致しないスライスの検査中にゴースト画像はＭＲ信
号の位相が傾斜磁界により影響される故になお発生す
る。これらの位相誤差はＥＰＩのような位相感応性ＭＲ
方法の場合には画像アーティファクト妨害を発生させ
る。
【０００４】
【発明が解決しようとする課題】本発明の目的は該位相
誤差が補正されるような方法で上記の種類の方法を提供
することにある。
誤差が補正されるような方法で上記の種類の方法を提供
することにある。
【０００５】
【課題を解決するための手段】この目的はＭＲ信号に含
まれる位相誤差の時間変動は傾斜磁界コイルシステムの
パルス応答からと同様に第一及び／又は第二の傾斜磁界
の課された時間変動から計算され、ＭＲ信号はそれによ
り補正される本発明により達成される。本発明は以下の
思考に基づいている。
まれる位相誤差の時間変動は傾斜磁界コイルシステムの
パルス応答からと同様に第一及び／又は第二の傾斜磁界
の課された時間変動から計算され、ＭＲ信号はそれによ
り補正される本発明により達成される。本発明は以下の
思考に基づいている。
【０００６】時間的な挙動が同一ではない２つ（又は３
つ）の傾斜磁界コイルの傾斜磁界の組み合わせの場合に
は、例えば所望の又は課された（第一又は第二の）傾斜
磁界の代わりに渦電流効果により方向及び可能なその強
度が課された傾斜磁界から逸脱した傾斜磁界が遷移相中
に発生される。所望の傾斜磁界からの実際の傾斜磁界の
この時間的偏差はスライスに接する成分とスライスの法
線（垂直な）方向の成分とに分解されうる傾斜磁界誤差
であると考えられる。
つ）の傾斜磁界コイルの傾斜磁界の組み合わせの場合に
は、例えば所望の又は課された（第一又は第二の）傾斜
磁界の代わりに渦電流効果により方向及び可能なその強
度が課された傾斜磁界から逸脱した傾斜磁界が遷移相中
に発生される。所望の傾斜磁界からの実際の傾斜磁界の
この時間的偏差はスライスに接する成分とスライスの法
線（垂直な）方向の成分とに分解されうる傾斜磁界誤差
であると考えられる。
【０００７】スライスに接する成分及びそれらの補正の
効果は上記文献に記載されている（ＥＰＩシーケンスに
関して）。本発明はスライスに垂直な成分及びその影響
の補正に関する。スライスに垂直な時間的に変動する傾
斜磁界は検査されているスライスがＭＲシステムのアイ
ソセンターの外側に存在する場合に（アイソセンターは
傾斜磁界コイルシステムにより発生された傾斜磁界がゼ
ロ値を有する空間の点である）ＭＲ信号での位相誤差を
引き起こす。これらの位相誤差はＭＲ方法（ＭＲ信号か
らのＭＲ画像の再構成を含む）が位相感応性である場合
にはＭＲ画像にアーティファクトを引き起こす（例えば
ＥＰＩシーケンスの場合にはゴースト画像）。本発明は
ＭＲ信号の位相誤差が補正されることによりこれらのア
ーティファクトを回避する。位相誤差の補正に対してそ
の時間的変動は傾斜磁界コイルシステムのパルス応答か
ら良く知られているように第一及び／又は第二の傾斜磁
界の所定の時間的変動から計算される。
効果は上記文献に記載されている（ＥＰＩシーケンスに
関して）。本発明はスライスに垂直な成分及びその影響
の補正に関する。スライスに垂直な時間的に変動する傾
斜磁界は検査されているスライスがＭＲシステムのアイ
ソセンターの外側に存在する場合に（アイソセンターは
傾斜磁界コイルシステムにより発生された傾斜磁界がゼ
ロ値を有する空間の点である）ＭＲ信号での位相誤差を
引き起こす。これらの位相誤差はＭＲ方法（ＭＲ信号か
らのＭＲ画像の再構成を含む）が位相感応性である場合
にはＭＲ画像にアーティファクトを引き起こす（例えば
ＥＰＩシーケンスの場合にはゴースト画像）。本発明は
ＭＲ信号の位相誤差が補正されることによりこれらのア
ーティファクトを回避する。位相誤差の補正に対してそ
の時間的変動は傾斜磁界コイルシステムのパルス応答か
ら良く知られているように第一及び／又は第二の傾斜磁
界の所定の時間的変動から計算される。
【０００８】本発明はＭＲ信号に含まれる位相情報がＭ
Ｒ画像の再構成に対して必要とされる全てのＭＲ方法に
対して原理的には用いられる。斯くして本発明の更なる
バージョンは第一の傾斜磁界は交番する極性を有する読
み取り傾斜磁界であり、第二の傾斜磁界は読み取り磁界
の一の極性から他の極性への遷移で専ら作用される位相
エンコード傾斜磁界であり、極性の逆転の後に発生する
ＭＲ信号はＭＲ画像を形成するためにその位相位置の補
正後に取得され、用いられる。これはいわゆるＥＰＩ方
法であり、位相誤差が本発明により回避されるゴースト
画像を引き起こす。
Ｒ画像の再構成に対して必要とされる全てのＭＲ方法に
対して原理的には用いられる。斯くして本発明の更なる
バージョンは第一の傾斜磁界は交番する極性を有する読
み取り傾斜磁界であり、第二の傾斜磁界は読み取り磁界
の一の極性から他の極性への遷移で専ら作用される位相
エンコード傾斜磁界であり、極性の逆転の後に発生する
ＭＲ信号はＭＲ画像を形成するためにその位相位置の補
正後に取得され、用いられる。これはいわゆるＥＰＩ方
法であり、位相誤差が本発明により回避されるゴースト
画像を引き起こす。
【０００９】本発明の更なるバージョンでは少なくとも
２つの傾斜磁界コイルシステムの傾斜磁界は傾斜磁界の
重畳が振幅が時間の関数として増加又は減少する回転読
み取り磁界をもたらすような時間依存的な方法で組み合
わされる。このバージョンでは、ｋ空間は螺旋状にスキ
ャンされ、螺旋形が不均一な角速度でトラックされる場
合には本発明の利点は特に顕著である。
２つの傾斜磁界コイルシステムの傾斜磁界は傾斜磁界の
重畳が振幅が時間の関数として増加又は減少する回転読
み取り磁界をもたらすような時間依存的な方法で組み合
わされる。このバージョンでは、ｋ空間は螺旋状にスキ
ャンされ、螺旋形が不均一な角速度でトラックされる場
合には本発明の利点は特に顕著である。
【００１０】本発明の更なるバージョンでは検査領域に
作用する複数のシーケンス及び少なくとも２つの傾斜磁
界コイルシステムの傾斜磁界はそれから得られる読み取
り傾斜磁界の方向が一のシーケンスから他のシーケンス
で変動する一方でその絶対値は一定のままであるような
方法で一のシーケンスから他のシーケンスへと組み合わ
され、変動され、この読み取り磁界と関連して発生する
ＭＲ信号（Ｓ（ｔ））はＭＲ画像を形成するためにその
位相位置の補正の後に取得され、用いられる。このＭＲ
方法は投影再構成法として知られており、また「絶対
値」画像の形成を可能にし、ここでＭＲ信号に含まれる
位相情報は問題とはならないが、本発明が好ましく用い
られる位相感応性の再構成をまた可能にする。
作用する複数のシーケンス及び少なくとも２つの傾斜磁
界コイルシステムの傾斜磁界はそれから得られる読み取
り傾斜磁界の方向が一のシーケンスから他のシーケンス
で変動する一方でその絶対値は一定のままであるような
方法で一のシーケンスから他のシーケンスへと組み合わ
され、変動され、この読み取り磁界と関連して発生する
ＭＲ信号（Ｓ（ｔ））はＭＲ画像を形成するためにその
位相位置の補正の後に取得され、用いられる。このＭＲ
方法は投影再構成法として知られており、また「絶対
値」画像の形成を可能にし、ここでＭＲ信号に含まれる
位相情報は問題とはならないが、本発明が好ましく用い
られる位相感応性の再構成をまた可能にする。
【００１１】上記の方法を実施するＭＲ装置は ａ） 均一で、一定の磁界を発生する主磁界磁石と、 ｂ） 異なる方向に延在する傾斜を有する傾斜磁界を発
生する傾斜磁界コイルシステムと、 ｃ） 少なくとも一つのＲＦパルスにより検査領域での
スライスでの核磁化を励起するＲＦコイルシステムと、 ｄ） 検査領域に発生されたＭＲ信号を取得する受信機
と、 ｅ） 補正されたＭＲ測定データから検査領域での核磁
化分布を再構成する再構成ユニットと、 ｆ） 傾斜磁界コイルシステム、ＲＦコイルシステム、
受信機を制御する制御ユニットとを含み、制御ユニット
は、少なくとも２つの傾斜磁界コイルシステムが制御ユ
ニットにより課された時間変動で傾斜磁界を同時に発生
するようプログラムされており、ＭＲ信号に含まれる位
相誤差の時間変動を計算する手段が設けられており、そ
の手段は傾斜磁界コイルシステムのパルス応答からと同
様に傾斜磁界が課された時間変動から位相誤差を計算
し、それによりＭＲ信号を補正することを特徴とする。
生する傾斜磁界コイルシステムと、 ｃ） 少なくとも一つのＲＦパルスにより検査領域での
スライスでの核磁化を励起するＲＦコイルシステムと、 ｄ） 検査領域に発生されたＭＲ信号を取得する受信機
と、 ｅ） 補正されたＭＲ測定データから検査領域での核磁
化分布を再構成する再構成ユニットと、 ｆ） 傾斜磁界コイルシステム、ＲＦコイルシステム、
受信機を制御する制御ユニットとを含み、制御ユニット
は、少なくとも２つの傾斜磁界コイルシステムが制御ユ
ニットにより課された時間変動で傾斜磁界を同時に発生
するようプログラムされており、ＭＲ信号に含まれる位
相誤差の時間変動を計算する手段が設けられており、そ
の手段は傾斜磁界コイルシステムのパルス応答からと同
様に傾斜磁界が課された時間変動から位相誤差を計算
し、それによりＭＲ信号を補正することを特徴とする。
【００１２】
【発明の実施の形態】本発明を以下に図を参照して詳細
に説明する。図１は検査領域（図示せず）内のｚ方向に
例えば約１．５テスラの強度の定常で本質的に均一な磁
界を発生する主磁界磁石１の概略を示す。ｚ方向はＭＲ
検査中に患者が収容される検査テーブル（図示せず）の
長手方向に向けられている。
に説明する。図１は検査領域（図示せず）内のｚ方向に
例えば約１．５テスラの強度の定常で本質的に均一な磁
界を発生する主磁界磁石１の概略を示す。ｚ方向はＭＲ
検査中に患者が収容される検査テーブル（図示せず）の
長手方向に向けられている。
【００１３】傾斜磁界コイルシステム２がまた設けら
れ、これは傾斜磁界Ｇx ，Ｇy ，Ｇzを発生させること
を可能にする３つのコイルシステムを含み、これらはｚ
方向に延在し、ｘ，ｙ，ｚ方向にそれぞれ傾斜磁界を有
する。この傾斜コイルシステムがｘ，ｙ，ｚ方向に一致
しない傾斜を有する傾斜磁界を発生するために用いられ
る場合には２つまたは３つの傾斜磁界コイルシステムは
同時に作用されなければならない。傾斜磁界の時間的変
動は波形発生器により換言すれば各ｘ，ｙ，ｚ方向に対
して別々に課される。波形発生器４は演算及び制御ユニ
ット５により制御され、これは検査方法（例えばＥＰＩ
シーケンス）及びその幾何的パラメータ（例えば画像化
されるスライスの位置）が入力された後に波形発生器に
ロードするために傾斜磁界Ｇx ，Ｇy ，Ｇz の必要な時
間変動が計算される。ＭＲ検査中にこれらの信号は波形
発生器から読み取られ、傾斜磁界増幅装置３に印加さ
れ、これはそれに基づいて傾斜磁界コイルシステム２に
対して必要とされる電流を発生する。
れ、これは傾斜磁界Ｇx ，Ｇy ，Ｇzを発生させること
を可能にする３つのコイルシステムを含み、これらはｚ
方向に延在し、ｘ，ｙ，ｚ方向にそれぞれ傾斜磁界を有
する。この傾斜コイルシステムがｘ，ｙ，ｚ方向に一致
しない傾斜を有する傾斜磁界を発生するために用いられ
る場合には２つまたは３つの傾斜磁界コイルシステムは
同時に作用されなければならない。傾斜磁界の時間的変
動は波形発生器により換言すれば各ｘ，ｙ，ｚ方向に対
して別々に課される。波形発生器４は演算及び制御ユニ
ット５により制御され、これは検査方法（例えばＥＰＩ
シーケンス）及びその幾何的パラメータ（例えば画像化
されるスライスの位置）が入力された後に波形発生器に
ロードするために傾斜磁界Ｇx ，Ｇy ，Ｇz の必要な時
間変動が計算される。ＭＲ検査中にこれらの信号は波形
発生器から読み取られ、傾斜磁界増幅装置３に印加さ
れ、これはそれに基づいて傾斜磁界コイルシステム２に
対して必要とされる電流を発生する。
【００１４】制御ユニット５はまた波形発生器６と協働
する。ワークステーションはＭＲ画像の表示用のモニタ
７を含む。入力はキーボード８又は例えばライトペンに
より対話的入力ユニット９を介してなされる；そのよう
な入力は例えばＭＲ検査中に用いられるＭＲ方法では、
画像化されるスライスの位置等々である。検査領域での
核磁化はＲＦコイル１０からのＲＦパルスにより励起さ
れえ、これはＲＦ送信機１２の出力信号を増幅するＲＦ
増幅器１１に接続される。ＲＦ送信機１２は発振器１３
により供給される搬送波発振でＲＦパルスの（複素）包
絡線を変調し、その周波数はラーマー周波数に対応する
（１．５テスラの主磁界の場合には約６３ＭＨｚ）。計
算及び制御ユニット５は（複素）包絡線を発振器１４に
ロードし、それは送信機１２に組み合わされている。検
査領域で発生されたＭＲ信号は受信コイル２０によりピ
ックアップされ、増幅器２１により増幅される。増幅さ
れたＭＲ信号は発振器の２つの９０゜オフセットされた
搬送波発振により直交復調器２２で復調され、これらは
複素ＭＲ信号の実数部分と虚数部分と見なしうる。これ
らの信号はアナログ／デジタル変換器２３に印加され、
これはＭＲデータを形成し、制御ユニット５により抑制
されずに供給される。ＭＲデータは再構成ユニット２４
で記憶され、補正され、これは複数のＭＲ信号から得ら
れたＭＲデータから検査領域の核磁化を表すＭＲ画像を
再構成するためにワークステーション６と協働する。
する。ワークステーションはＭＲ画像の表示用のモニタ
７を含む。入力はキーボード８又は例えばライトペンに
より対話的入力ユニット９を介してなされる；そのよう
な入力は例えばＭＲ検査中に用いられるＭＲ方法では、
画像化されるスライスの位置等々である。検査領域での
核磁化はＲＦコイル１０からのＲＦパルスにより励起さ
れえ、これはＲＦ送信機１２の出力信号を増幅するＲＦ
増幅器１１に接続される。ＲＦ送信機１２は発振器１３
により供給される搬送波発振でＲＦパルスの（複素）包
絡線を変調し、その周波数はラーマー周波数に対応する
（１．５テスラの主磁界の場合には約６３ＭＨｚ）。計
算及び制御ユニット５は（複素）包絡線を発振器１４に
ロードし、それは送信機１２に組み合わされている。検
査領域で発生されたＭＲ信号は受信コイル２０によりピ
ックアップされ、増幅器２１により増幅される。増幅さ
れたＭＲ信号は発振器の２つの９０゜オフセットされた
搬送波発振により直交復調器２２で復調され、これらは
複素ＭＲ信号の実数部分と虚数部分と見なしうる。これ
らの信号はアナログ／デジタル変換器２３に印加され、
これはＭＲデータを形成し、制御ユニット５により抑制
されずに供給される。ＭＲデータは再構成ユニット２４
で記憶され、補正され、これは複数のＭＲ信号から得ら
れたＭＲデータから検査領域の核磁化を表すＭＲ画像を
再構成するためにワークステーション６と協働する。
【００１５】図２は本発明が好ましく用いられるＥＰＩ
タイプのＭＲシーケンスを示す。このシーケンスは２つ
のＲＦ励起パルスＲＦ１，ＲＦ２（第一のライン）から
なり、それはｗ方向での傾斜磁界Ｇw （第二のライン）
と共にｗ方向に垂直なスライスの核磁化をその定常状態
からそれぞれ９０゜と１８０゜回転する。続いてｕ方向
に延在する交番読み取り傾斜磁界Ｇu （第三のライン）
がｖ方向に延在し、読み取り傾斜磁界Ｇu の極性逆転中
に専ら作用する位相エンコード傾斜磁界Ｇv （第四のラ
イン）と同様に印加される。読み取り傾斜磁界Ｇu が正
又は負の極性を有するときはいつもＭＲ信号が取得され
（第五のライン）；斯くして計算及び制御ユニット５は
受信されたＭＲ信号がＭＲデータに変換されるようアナ
ログ／デジタル変換器２３を作動させる。
タイプのＭＲシーケンスを示す。このシーケンスは２つ
のＲＦ励起パルスＲＦ１，ＲＦ２（第一のライン）から
なり、それはｗ方向での傾斜磁界Ｇw （第二のライン）
と共にｗ方向に垂直なスライスの核磁化をその定常状態
からそれぞれ９０゜と１８０゜回転する。続いてｕ方向
に延在する交番読み取り傾斜磁界Ｇu （第三のライン）
がｖ方向に延在し、読み取り傾斜磁界Ｇu の極性逆転中
に専ら作用する位相エンコード傾斜磁界Ｇv （第四のラ
イン）と同様に印加される。読み取り傾斜磁界Ｇu が正
又は負の極性を有するときはいつもＭＲ信号が取得され
（第五のライン）；斯くして計算及び制御ユニット５は
受信されたＭＲ信号がＭＲデータに変換されるようアナ
ログ／デジタル変換器２３を作動させる。
【００１６】図４はｘ，ｙ，ｚ座標系での相互に垂直な
３つの方向ｕ，ｖ，ｗを示す。読み取り傾斜磁界の方向
ｕ及び位相エンコーディング傾斜磁界のｖ方向はパルス
ＲＦ１，ＲＦ２により励起されたスライスへの垂線が延
在するｗ方向に関してと同様にｚ方向に関して９０゜以
外の角度をなす。これは傾斜磁界コイルシステム２の
ｘ，ｙ，ｚ方向に対してコイルシステムを傾斜磁界
Ｇu ，Ｇv ，Ｇw の一つを発生させるために同時に作用
させなければならない。図５の変換マトリックスはそれ
に要求される寄与の比を示す。
３つの方向ｕ，ｖ，ｗを示す。読み取り傾斜磁界の方向
ｕ及び位相エンコーディング傾斜磁界のｖ方向はパルス
ＲＦ１，ＲＦ２により励起されたスライスへの垂線が延
在するｗ方向に関してと同様にｚ方向に関して９０゜以
外の角度をなす。これは傾斜磁界コイルシステム２の
ｘ，ｙ，ｚ方向に対してコイルシステムを傾斜磁界
Ｇu ，Ｇv ，Ｇw の一つを発生させるために同時に作用
させなければならない。図５の変換マトリックスはそれ
に要求される寄与の比を示す。
【００１７】一方でこの変換マトリックスは傾斜磁界Ｇ
u ，Ｇv ，Ｇw を形成する成分を示す。このマトリック
スによれば例えば傾斜磁界Ｇu は傾斜磁界Ｇx ，Ｇy ，
Ｇzがこの傾斜磁界の強度の成分ａ1 ，ｂ1 ，ｃ1 でそ
れぞれ作用する場合には所定の傾斜磁界強度で得られ
る。変換マトリックスはまた傾斜磁界Ｇx ，Ｇy ，又は
Ｇz の一つのみが作用される場合に傾斜磁界Ｇu ，
Ｇv ，又はＧw がどのくらい大きいかを示す（例えば強
度１の傾斜磁界がｚ方向で作用された場合には傾斜磁界
の成分はｕ，ｖ，ｗ方向で相対的強度ｃ1 ，ｃ2 ，ｃ3
で発生する）。要素ａ 1 ，ａ2 ，ａ3 は軸ｘに関して軸
ｕ，ｖ，ｗがなす角のそれぞれのコサインにより決定さ
れ、同様に要素ｂ1 ，ｂ2 ，ｂ3 は軸ｙに関する軸ｕ，
ｖ，ｗの方向の方向余弦を表し、ｃ1 ，ｃ2 ，ｃ3 は一
方でｕ，ｖ，ｗと他方でｚとの間の方向余弦を表す。知
られているように変換マトリックスの一つの行又は列に
含まれる３つの要素の二乗の和は１である。
u ，Ｇv ，Ｇw を形成する成分を示す。このマトリック
スによれば例えば傾斜磁界Ｇu は傾斜磁界Ｇx ，Ｇy ，
Ｇzがこの傾斜磁界の強度の成分ａ1 ，ｂ1 ，ｃ1 でそ
れぞれ作用する場合には所定の傾斜磁界強度で得られ
る。変換マトリックスはまた傾斜磁界Ｇx ，Ｇy ，又は
Ｇz の一つのみが作用される場合に傾斜磁界Ｇu ，
Ｇv ，又はＧw がどのくらい大きいかを示す（例えば強
度１の傾斜磁界がｚ方向で作用された場合には傾斜磁界
の成分はｕ，ｖ，ｗ方向で相対的強度ｃ1 ，ｃ2 ，ｃ3
で発生する）。要素ａ 1 ，ａ2 ，ａ3 は軸ｘに関して軸
ｕ，ｖ，ｗがなす角のそれぞれのコサインにより決定さ
れ、同様に要素ｂ1 ，ｂ2 ，ｂ3 は軸ｙに関する軸ｕ，
ｖ，ｗの方向の方向余弦を表し、ｃ1 ，ｃ2 ，ｃ3 は一
方でｕ，ｖ，ｗと他方でｚとの間の方向余弦を表す。知
られているように変換マトリックスの一つの行又は列に
含まれる３つの要素の二乗の和は１である。
【００１８】Ｇu ，Ｇv ，Ｇw の時間変動は図２のライ
ン２から４に示されるように制御ユニット５又は波形発
生器４を介して課されるようなこれらの量の時間変動を
表す。傾斜磁界の実際の変化は部分３、２（図１を参
照）の周波数応答の故にそれから逸脱する。この影響及
びその結果は図３の（ａ），（ｂ），（ｃ）を参照して
以下に詳細に説明する。図３の（ａ）の実線はＧu の時
間変動を示し、即ち所望の又は予め設定（制御ユニット
５又は波形発生器４により）されたｕ方向の傾斜磁界の
変動を表す。図３の（ａ）の破線は傾斜磁界の実際の変
化Ｇurを示し、これは種々の傾斜磁界システムが遅延な
しに課された時間変動に追従できないためである。異な
る遅延がｘ方向、ｙ方向、ｚ方向に対するそれぞれのチ
ャンネルで発生する場合に問題が発生する。
ン２から４に示されるように制御ユニット５又は波形発
生器４を介して課されるようなこれらの量の時間変動を
表す。傾斜磁界の実際の変化は部分３、２（図１を参
照）の周波数応答の故にそれから逸脱する。この影響及
びその結果は図３の（ａ），（ｂ），（ｃ）を参照して
以下に詳細に説明する。図３の（ａ）の実線はＧu の時
間変動を示し、即ち所望の又は予め設定（制御ユニット
５又は波形発生器４により）されたｕ方向の傾斜磁界の
変動を表す。図３の（ａ）の破線は傾斜磁界の実際の変
化Ｇurを示し、これは種々の傾斜磁界システムが遅延な
しに課された時間変動に追従できないためである。異な
る遅延がｘ方向、ｙ方向、ｚ方向に対するそれぞれのチ
ャンネルで発生する場合に問題が発生する。
【００１９】図５の変換マトリックスからわかるように
ｕ方向に延在する傾斜を有する傾斜磁界は傾斜磁界
Ｇx ，Ｇy ，Ｇz の強度が変換マトリックスの要素
ａ1 ，ｂ1 ，ｃ1 に関係する場合にのみ発生する。この
条件は傾斜磁界コイルシステム２により発生された傾斜
磁界Ｇx ，Ｇy ，Ｇz が波形発生器４により課された値
に到達せず、異なる時間的挙動でこの値に接近する場合
に限り満たされない。故にこれらの繰り返しの間にｖ方
向又はｗ方向に傾斜を有する傾斜磁界の成分がまた存在
する。これは例えば問題の瞬間にｗ方向に傾斜磁界が課
されない場合にすら傾斜磁界Ｇu は選択されたスライス
に垂直に延在するような傾斜磁界Ｇw により加えられる
ことを意味する。
ｕ方向に延在する傾斜を有する傾斜磁界は傾斜磁界
Ｇx ，Ｇy ，Ｇz の強度が変換マトリックスの要素
ａ1 ，ｂ1 ，ｃ1 に関係する場合にのみ発生する。この
条件は傾斜磁界コイルシステム２により発生された傾斜
磁界Ｇx ，Ｇy ，Ｇz が波形発生器４により課された値
に到達せず、異なる時間的挙動でこの値に接近する場合
に限り満たされない。故にこれらの繰り返しの間にｖ方
向又はｗ方向に傾斜を有する傾斜磁界の成分がまた存在
する。これは例えば問題の瞬間にｗ方向に傾斜磁界が課
されない場合にすら傾斜磁界Ｇu は選択されたスライス
に垂直に延在するような傾斜磁界Ｇw により加えられる
ことを意味する。
【００２０】この傾斜磁界Ｇw （ｔ）の時間変動は以下
の関係に従う： Ｇw （ｔ）＝ｃ3 （ｃ1 Ｇu （ｔ）＋ｃ2 Ｇv （ｔ））＊（Ｉz −Ｉx ）＋ ｂ3 （ｂ1 Ｇu （ｔ）＋ｂ2 Ｇv （ｔ））＊（Ｉy −Ｉx ） （１） ここで、ｃ1 ．．．ｃ3 及びｂ1 ．．．ｂ3 は図５の変
換マトリックスの要素であり、一方でＧu （ｔ）及びＧ
v （ｔ）はそれぞれｕ及びｖ方向の傾斜磁界の課された
（所望の）時間変化を表す。演算子＊はコンボリューシ
ョンを表す。Ｉx，Ｉy ，Ｉz はユニット２、３のｘ，
ｙ，ｚ方向に対するチャンネルのパルス応答である。
パルス応答は傾斜磁界システム２及び増幅器システム３
の構成により課される。それはいわゆるディラックのパ
ルス、即ち急激に高い値に増加し再び急速に消失するパ
ルスが入力に与えられたときに傾斜磁界の時間変動を表
す。そのようなディラックのパルスは実際には発生でき
ない。何故ならばそれは無限の振幅とゼロの持続時間を
有さなければならないからである。パルス応答は傾斜磁
界増幅器３の一つの入力上の適切な励起変数（入力変
数）の時間変化及び傾斜磁界増幅器に接続された傾斜磁
界コイル２により順次発生された磁界傾斜（出力変数）
の時間変化を測定することにより決定されうる。傾斜磁
界増幅器３及び傾斜磁界コイル２からなるシステムが線
形であること（この条件は実際に満たされる）を仮定す
ると線形システムの出力変数は入力変数とパルス応答と
のコンボリューションにより与えられるという事実を用
いることができる。パルス応答は例えばコンボリューシ
ョン操作を打ち消す「デコンボリューション」によりそ
れから計算可能である。ｘ，ｙ，ｚ方向に対するチャン
ネルに対してこのように決定されたパルス応答Ｉx ，Ｉ
y ，Ｉz は記憶される。
の関係に従う： Ｇw （ｔ）＝ｃ3 （ｃ1 Ｇu （ｔ）＋ｃ2 Ｇv （ｔ））＊（Ｉz −Ｉx ）＋ ｂ3 （ｂ1 Ｇu （ｔ）＋ｂ2 Ｇv （ｔ））＊（Ｉy −Ｉx ） （１） ここで、ｃ1 ．．．ｃ3 及びｂ1 ．．．ｂ3 は図５の変
換マトリックスの要素であり、一方でＧu （ｔ）及びＧ
v （ｔ）はそれぞれｕ及びｖ方向の傾斜磁界の課された
（所望の）時間変化を表す。演算子＊はコンボリューシ
ョンを表す。Ｉx，Ｉy ，Ｉz はユニット２、３のｘ，
ｙ，ｚ方向に対するチャンネルのパルス応答である。
パルス応答は傾斜磁界システム２及び増幅器システム３
の構成により課される。それはいわゆるディラックのパ
ルス、即ち急激に高い値に増加し再び急速に消失するパ
ルスが入力に与えられたときに傾斜磁界の時間変動を表
す。そのようなディラックのパルスは実際には発生でき
ない。何故ならばそれは無限の振幅とゼロの持続時間を
有さなければならないからである。パルス応答は傾斜磁
界増幅器３の一つの入力上の適切な励起変数（入力変
数）の時間変化及び傾斜磁界増幅器に接続された傾斜磁
界コイル２により順次発生された磁界傾斜（出力変数）
の時間変化を測定することにより決定されうる。傾斜磁
界増幅器３及び傾斜磁界コイル２からなるシステムが線
形であること（この条件は実際に満たされる）を仮定す
ると線形システムの出力変数は入力変数とパルス応答と
のコンボリューションにより与えられるという事実を用
いることができる。パルス応答は例えばコンボリューシ
ョン操作を打ち消す「デコンボリューション」によりそ
れから計算可能である。ｘ，ｙ，ｚ方向に対するチャン
ネルに対してこのように決定されたパルス応答Ｉx ，Ｉ
y ，Ｉz は記憶される。
【００２１】スライス面に延在する傾斜磁界Ｇ
u （ｔ），Ｇv （ｔ）への印加おいてスライスに垂直に
延在するＧw （ｔ）はｘ，ｙ，ｚ方向に対する３つのチ
ャンネルが同一のパルス応答Ｉx ，Ｉy ，Ｉz を有する
場合には発生しないことは式（１）から明らかである。
この条件は通常は満足されない。同一の時間変化が通常
ｘとｙ方向に対して発生する（即ち式（１）の２つの加
数の二番目が消える）場合にすらｘ，ｙ方向に対する傾
斜磁界システム２、３の時間的挙動はｚ方向に対する時
間的挙動から逸脱する。式（１）から、後者の場合には
スライス面に垂直に延在する傾斜磁界はｃ3 がゼロ、即
ちスライスの法線ｗがｚ方向に垂直に延在するいずれか
の場合に発生しないことがまた明らかである。
u （ｔ），Ｇv （ｔ）への印加おいてスライスに垂直に
延在するＧw （ｔ）はｘ，ｙ，ｚ方向に対する３つのチ
ャンネルが同一のパルス応答Ｉx ，Ｉy ，Ｉz を有する
場合には発生しないことは式（１）から明らかである。
この条件は通常は満足されない。同一の時間変化が通常
ｘとｙ方向に対して発生する（即ち式（１）の２つの加
数の二番目が消える）場合にすらｘ，ｙ方向に対する傾
斜磁界システム２、３の時間的挙動はｚ方向に対する時
間的挙動から逸脱する。式（１）から、後者の場合には
スライス面に垂直に延在する傾斜磁界はｃ3 がゼロ、即
ちスライスの法線ｗがｚ方向に垂直に延在するいずれか
の場合に発生しないことがまた明らかである。
【００２２】図３の（ｂ）は図３の（ａ）に従って傾斜
磁界Ｇu により引き起こされた傾斜磁界Ｇw （ｔ）の時
間変化を示す。攪乱傾斜磁界Ｇx （ｔ）はｕ方向での実
際の傾斜磁界Ｇurと同時に発生し、即ち後者は課された
時間変化Ｇu から逸脱することがわかる。攪乱成分Ｇw
（ｔ）の発生はいわゆるｋ空間がｕ，ｖ平面に沿ってス
キャンされず、平坦でない面に沿ってスキャンされるこ
とを意味する。ｋ空間での成分ｋ w （ｔ）に対して以下
のようになる。
磁界Ｇu により引き起こされた傾斜磁界Ｇw （ｔ）の時
間変化を示す。攪乱傾斜磁界Ｇx （ｔ）はｕ方向での実
際の傾斜磁界Ｇurと同時に発生し、即ち後者は課された
時間変化Ｇu から逸脱することがわかる。攪乱成分Ｇw
（ｔ）の発生はいわゆるｋ空間がｕ，ｖ平面に沿ってス
キャンされず、平坦でない面に沿ってスキャンされるこ
とを意味する。ｋ空間での成分ｋ w （ｔ）に対して以下
のようになる。
【００２３】
【数１】
【００２４】ここでγは核磁気回転比（約４２．５ＭＨ
ｚ／Ｔ）である。励起されたスライスが傾斜磁界コイル
システムのアイソセンターを通って延在しない場合には
位相シフトφ（ｔ）が以下の式に従って発生する。 φ（ｔ）＝ｋw （ｔ）Ｌ （５） ここでＬは励起されたスライスとアイソセンターとの間
の距離である。
ｚ／Ｔ）である。励起されたスライスが傾斜磁界コイル
システムのアイソセンターを通って延在しない場合には
位相シフトφ（ｔ）が以下の式に従って発生する。 φ（ｔ）＝ｋw （ｔ）Ｌ （５） ここでＬは励起されたスライスとアイソセンターとの間
の距離である。
【００２５】斯くしてこれはアイソセンターから離れて
位置する励起されたオブリークスライスに現れるＭＲ信
号Ｓ（ｔ）は以下の関係に従うことを意味する：
位置する励起されたオブリークスライスに現れるＭＲ信
号Ｓ（ｔ）は以下の関係に従うことを意味する：
【００２６】
【数２】
【００２７】ここでＳ0 （ｔ）は位相シフトが攪乱傾斜
磁界Ｇw （ｔ）により引き起こされない場合に発生する
ＭＲ信号である。図３の（ｃ）は式（２）、（３）から
得られた位相誤差φの時間変化を示す。Ｇurが増加中に
位相誤差φは一定の値に増加し、Ｇurが減少すると共に
再び値ゼロに減少する。原理的に同じことがＧu と同様
にＧv にも言える。しかしながらこれらのブリップのリ
ーディング及びトレーリングエッジが相互に直接続く故
に位相誤差はＭＲ信号の取得（図２の第五ライン）です
でに消失している。故にＧv に基づく成分はＥＰＩシー
ケンスの場合には無視しうる。
磁界Ｇw （ｔ）により引き起こされない場合に発生する
ＭＲ信号である。図３の（ｃ）は式（２）、（３）から
得られた位相誤差φの時間変化を示す。Ｇurが増加中に
位相誤差φは一定の値に増加し、Ｇurが減少すると共に
再び値ゼロに減少する。原理的に同じことがＧu と同様
にＧv にも言える。しかしながらこれらのブリップのリ
ーディング及びトレーリングエッジが相互に直接続く故
に位相誤差はＭＲ信号の取得（図２の第五ライン）です
でに消失している。故にＧv に基づく成分はＥＰＩシー
ケンスの場合には無視しうる。
【００２８】上記の考察により図６を参照して記載され
るように受信されたＭＲデータＳ（ｔ）に含まれる位相
誤差の補正に対するフローチャートが得られる。ブロッ
ク１０１は初期化及び就中ＲＦパルスＲＦ１，ＲＦ２の
発生からなる。続いてそれぞれのＭＲ信号は読み取り傾
斜磁界Ｇu の負の極性と同様に正の極性で取得され、こ
の信号は一連のＭＲデータに変換される（ブロック１０
２）。スライス面に垂直に延在する傾斜磁界Ｇw （ｔ）
の時間変化はブロック１０３で式（１）に従って計算さ
れる。斯くして計算された値から得られた値ｋw （ｔ）
が式（２）により計算され、この値はブロック１０５で
式（３）に従って位相φに対する値を計算するために用
いられる。ブロック１０６では信号Ｓ（ｔ）の測定され
た時間変化は計算された位相誤差と測定されたＭＲ信号
Ｓ（ｔ）の時間変化から得られた信号位相との加算によ
り計算された時間変化に従って補正される。式（４）に
従ってこの結果はＭＲ信号Ｓ0 （ｔ）で得られ、これか
ら位相誤差が除去される。
るように受信されたＭＲデータＳ（ｔ）に含まれる位相
誤差の補正に対するフローチャートが得られる。ブロッ
ク１０１は初期化及び就中ＲＦパルスＲＦ１，ＲＦ２の
発生からなる。続いてそれぞれのＭＲ信号は読み取り傾
斜磁界Ｇu の負の極性と同様に正の極性で取得され、こ
の信号は一連のＭＲデータに変換される（ブロック１０
２）。スライス面に垂直に延在する傾斜磁界Ｇw （ｔ）
の時間変化はブロック１０３で式（１）に従って計算さ
れる。斯くして計算された値から得られた値ｋw （ｔ）
が式（２）により計算され、この値はブロック１０５で
式（３）に従って位相φに対する値を計算するために用
いられる。ブロック１０６では信号Ｓ（ｔ）の測定され
た時間変化は計算された位相誤差と測定されたＭＲ信号
Ｓ（ｔ）の時間変化から得られた信号位相との加算によ
り計算された時間変化に従って補正される。式（４）に
従ってこの結果はＭＲ信号Ｓ0 （ｔ）で得られ、これか
ら位相誤差が除去される。
【００２９】計算段階１０３及び１０５は値
（ａ1 ．．．ｃ3 ；Ｇu （ｔ）．．．；Ｉx ．．．
Ｉz ）に基づき、これはＭＲ検査の始めに既に知られて
いる。故にこの計算はＭＲ信号の取得（ブロック１０
２）の前に既になされていることが可能であり、それに
よりＭＲ信号は測定された後にすぐに補正されうる。ブ
ロック１０７で励起されたスライス内の核磁化分布Ｍ
（ｕ，ｖ）は補正されたＭＲ信号Ｓ0 （ｔ）から再構成
される。攪乱傾斜磁界Ｇw （ｔ）を読み取り磁界Ｇ
u （ｔ）に関して時間的な関係にするのと同じ効果がま
たｖ方向に延在する望ましくない傾斜磁界が発生され、
位相エンコードを誤らせることを考慮に入れなければな
らない。この誤差の導入は上記の引用文献に記載される
ように、従って傾斜磁界Ｇv のブリップを拡大又は減少
することにより、又は例えば知られているグリッディン
グ方法のようなｋ空間内のＭＲデータの均一な分布を仮
定しない再構成方法を用いることにより対策が可能であ
る。この方法の実行はブロック１０８で完了する。
（ａ1 ．．．ｃ3 ；Ｇu （ｔ）．．．；Ｉx ．．．
Ｉz ）に基づき、これはＭＲ検査の始めに既に知られて
いる。故にこの計算はＭＲ信号の取得（ブロック１０
２）の前に既になされていることが可能であり、それに
よりＭＲ信号は測定された後にすぐに補正されうる。ブ
ロック１０７で励起されたスライス内の核磁化分布Ｍ
（ｕ，ｖ）は補正されたＭＲ信号Ｓ0 （ｔ）から再構成
される。攪乱傾斜磁界Ｇw （ｔ）を読み取り磁界Ｇ
u （ｔ）に関して時間的な関係にするのと同じ効果がま
たｖ方向に延在する望ましくない傾斜磁界が発生され、
位相エンコードを誤らせることを考慮に入れなければな
らない。この誤差の導入は上記の引用文献に記載される
ように、従って傾斜磁界Ｇv のブリップを拡大又は減少
することにより、又は例えば知られているグリッディン
グ方法のようなｋ空間内のＭＲデータの均一な分布を仮
定しない再構成方法を用いることにより対策が可能であ
る。この方法の実行はブロック１０８で完了する。
【００３０】本発明はＥＰＩシーケンスを参照して記載
されてきた。しかしながら本発明は画像情報が位相感応
的な方法により得られる他のＭＲ方法に対してもまた用
いられる。望ましくない傾斜磁界Ｇw （ｔ）により導入
されたアーティファクト（ゴースト画像、不鮮明、等
々）は問題のＭＲ方法に依存する。
されてきた。しかしながら本発明は画像情報が位相感応
的な方法により得られる他のＭＲ方法に対してもまた用
いられる。望ましくない傾斜磁界Ｇw （ｔ）により導入
されたアーティファクト（ゴースト画像、不鮮明、等
々）は問題のＭＲ方法に依存する。
【図１】本発明を実施するために適切なＭＲ装置のブロ
ック図である。
ック図である。
【図２】本発明が好ましく用いられるシーケンスでの種
々の信号の時間変化を示す。
々の信号の時間変化を示す。
【図３】拡大した尺度でそのようなシーケンスの更なる
信号の時間変化を示す。
信号の時間変化を示す。
【図４】ｘ，ｙ，ｚ座標系に関して傾斜磁界の方向ｕ，
ｖ，ｗを示す。
ｖ，ｗを示す。
【図５】２つの座標系相互の変換に対する変換マトリッ
クスを示す。
クスを示す。
【図６】本発明による信号処理を示すフローチャートで
ある。
ある。
１ 磁石 ２ 傾斜磁界コイルシステム ４ 波形発生器 ５ 制御ユニット ３ 傾斜磁界増幅装置 ７ モニタ ８ キーボード ９ 対話的入力ユニット １０ ＲＦコイル １２ ＲＦ送信機 １１ ＲＦ増幅器 １３、１４ 発振器 ２０ 受信コイル ２１ 増幅器 ２２ 直交復調器 ２３ アナログ／デジタル変換器 Ｇu ，Ｇv ，Ｇw 傾斜磁界 ＲＦ１，ＲＦ２ パルス
Claims (5)
- 【請求項１】 検査領域のスライスでの核磁化が均一で
定常な磁界の存在下で励起され、その後にスライスに接
する第一の方向（ｕ）に延在する傾斜を有する少なくと
も一つの第一の傾斜磁界（Ｇu ）が発生され、スライス
に接する第二の方向（ｖ）に延在する傾斜を有する可能
な第二の傾斜磁界（Ｇv ）が発生され、ＭＲ信号（Ｓ
（ｔ））は検査領域から受信され、少なくとも２つの傾
斜磁界コイルシステムの傾斜磁界（Ｇx ，Ｇy ，Ｇz ）
は第一と第二の傾斜磁界を発生するために組み合わされ
る少なくとも一つのシーケンスを含むＭＲ方法であっ
て、ＭＲ信号に含まれる位相誤差（φ）の時間変動は、
傾斜磁界コイルシステムのパルス応答（Ｉx ，Ｉy ，Ｉ
z ）からと同様に第一及び／又は第二の傾斜磁界
（Ｇ u 、Ｇv ）の課された時間変動から計算され、また
ＭＲ信号（Ｓ（ｔ））はそれにより補正されることを特
徴とするＭＲ方法。 - 【請求項２】 第一の傾斜磁界（Ｇu ）は交番する極性
を有する読み取り傾斜磁界であり、第二の傾斜磁界（Ｇ
v ）は読み取り磁界の一の極性から他の極性への遷移で
専ら作用される位相エンコード傾斜磁界であり、極性反
転の後に発生するＭＲ信号はＭＲ画像を形成するために
その位相位置の補正後に取得され、用いられることを特
徴とする請求項１記載のＭＲ方法。 - 【請求項３】 少なくとも２つの傾斜磁界コイルシステ
ムの傾斜磁界（Ｇx，Ｇy ，Ｇz ）は傾斜磁界の重畳が
振幅が時間の関数として増加又は減少する回転読み取り
磁界（Ｇu ）をもたらすような時間依存的な方法で組み
合わされることを特徴とする請求項１記載のＭＲ方法。 - 【請求項４】 複数のシーケンスが検査領域に作用し、
少なくとも２つの傾斜磁界コイルシステムの傾斜磁界
（Ｇx ，Ｇy ，Ｇz ）はそれから得られる読み取り傾斜
磁界の方向が一のシーケンスから他のシーケンスに変わ
るように一のシーケンスから他のシーケンスへと組み合
わされて変えられ、一方、その絶対値は一定のままであ
り、この読み取り磁界と関連して発生するＭＲ信号（Ｓ
（ｔ））はＭＲ画像を形成するためにその位相位置の補
正の後に取得され、用いられることを特徴とする請求項
１記載のＭＲ方法。 - 【請求項５】ａ） 均一で、一定の磁界を発生する主磁
界磁石（１）と、 ｂ） 異なる方向に延在する傾斜を有する傾斜磁界を発
生する傾斜磁界コイルシステム（２．．．４）と、 ｃ） 少なくとも一つのＲＦパルス（ＲＦ）により検査
領域のスライスでの核磁化を励起するＲＦコイルシステ
ム（１０）と、 ｄ） 検査領域に発生されたＭＲ信号を取得する受信機
（９）と、 ｅ） 補正されたＭＲ測定データから検査領域での核磁
化分布を再構成する再構成ユニット（２４）と、 ｆ） 傾斜磁界コイルシステム（２．．．４）、ＲＦコ
イルシステム（１０）、受信機（２１、２２）を制御す
る制御ユニット（５）とを含む請求項１記載の方法を実
施するＭＲ装置であって、制御ユニット（５）は、少な
くとも２つの傾斜磁界コイルシステムが制御ユニットに
より課された時間変動で傾斜磁界（Ｇu 、Ｇv ）を同時
に発生するようプログラムされており、ＭＲ信号に含ま
れる位相誤差（φ）の時間変動を計算する手段が設けら
れており、その手段は傾斜磁界コイルシステムのパルス
応答（Ｉx ，Ｉy ，Ｉz ）からと同様に傾斜磁界から課
された時間変動から位相誤差（φ）を計算し、それによ
りＭＲ信号（Ｓ（ｔ））を補正することを特徴とするＭ
Ｒ装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE19631845:9 | 1996-08-07 | ||
| DE19631845 | 1996-08-07 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JPH1071135A true JPH1071135A (ja) | 1998-03-17 |
| JPH1071135A5 JPH1071135A5 (ja) | 2005-05-26 |
| JP3898294B2 JP3898294B2 (ja) | 2007-03-28 |
Family
ID=7802006
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP20933897A Expired - Fee Related JP3898294B2 (ja) | 1996-08-07 | 1997-08-04 | Ｍｒ方法及び装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5929638A (ja) |
| EP (1) | EP0823642B1 (ja) |
| JP (1) | JP3898294B2 (ja) |
| DE (1) | DE59709890D1 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6201987B1 (en) * | 1998-05-26 | 2001-03-13 | General Electric Company | Error compensation for device tracking systems employing electromagnetic fields |
| DE19826864A1 (de) * | 1998-06-17 | 1999-12-23 | Philips Patentverwaltung | MR-Verfahren |
| CA2648310C (en) * | 2006-04-12 | 2016-10-18 | Thomson Licensing | Virtual dvd on demand and electronic dvd rental/buy/burn |
| EP3447520A1 (en) * | 2017-08-22 | 2019-02-27 | Koninklijke Philips N.V. | Data-driven correction of phase depending artefacts in a magnetic resonance imaging system |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4585995A (en) * | 1984-04-19 | 1986-04-29 | Technicare Corporation | Nuclear magnetic resonance eddy field suppression apparatus |
| EP0216523A3 (en) * | 1985-08-27 | 1989-04-05 | Resonex, Inc. | Process for non-orthogonal nmr imaging |
| US4698591A (en) * | 1986-01-03 | 1987-10-06 | General Electric Company | Method for magnetic field gradient eddy current compensation |
-
1997
- 1997-07-18 US US08/896,759 patent/US5929638A/en not_active Expired - Fee Related
- 1997-08-04 JP JP20933897A patent/JP3898294B2/ja not_active Expired - Fee Related
- 1997-08-04 EP EP97202430A patent/EP0823642B1/de not_active Expired - Lifetime
- 1997-08-04 DE DE59709890T patent/DE59709890D1/de not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| EP0823642B1 (de) | 2003-04-23 |
| EP0823642A1 (de) | 1998-02-11 |
| JP3898294B2 (ja) | 2007-03-28 |
| DE59709890D1 (de) | 2003-05-28 |
| US5929638A (en) | 1999-07-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4229487B2 (ja) | マクスウェル項誤差を補正する方法 | |
| CN104345289B (zh) | 用于确定水的t1时间和脂肪的t1时间的方法和磁共振设备 | |
| JP5198859B2 (ja) | 位相マッピングと、位相基準として用いる基準媒体が関係するｍｒｉ温度測定 | |
| US4672320A (en) | Imaging apparatus and method using nuclear magnetic resonance | |
| US6259250B1 (en) | Method and apparatus for reducing artifacts in echo planar imaging | |
| CN100415161C (zh) | 磁共振摄像装置以及磁共振摄像方法 | |
| US6285187B1 (en) | Method and apparatus for reducing artifacts in echo planar imaging | |
| JP2000070239A (ja) | 磁気共鳴イメ―ジング方法および磁気共鳴イメ―ジング装置 | |
| JPWO2004004563A1 (ja) | 磁気共鳴イメージング装置及び渦電流補償導出方法 | |
| US6906515B2 (en) | Magnetic resonance imaging device and method | |
| JP4052721B2 (ja) | 核磁気共鳴システムにより発生されるマクスウェル項位相誤差を補正するシステム及び核磁気共鳴システムにより位相コントラスト・アンジオグラムを形成するシステム | |
| JP2773840B2 (ja) | 渦電流起因位相歪修正システム | |
| JPH0747023B2 (ja) | 核磁気共鳴を用いた検査装置 | |
| JPS62254744A (ja) | 物体の１部分の核磁化分布を測定する方法及び装置 | |
| JP4229476B2 (ja) | Ｍｒ方法 | |
| US8493067B2 (en) | Magnetic resonance system and method to create a magnetic resonance image data set by radial scanning of a magnetic resonance system | |
| JP2003052668A (ja) | 核スピントモグラフィ装置及びその作動方法 | |
| JP3898294B2 (ja) | Ｍｒ方法及び装置 | |
| JP3588690B2 (ja) | 磁気共鳴装置 | |
| JPH0687845B2 (ja) | 磁気共鳴反転回復像の位相修正方法とその装置 | |
| US7187171B2 (en) | Magnetic resonance imaging apparatus and its control method | |
| JP3112930B2 (ja) | 磁気共鳴装置 | |
| JPH11290289A (ja) | 核磁気共鳴を用いた検査装置 | |
| JPH0568672A (ja) | 磁気共鳴イメージング装置における傾斜磁場印加方法 | |
| JP4146735B2 (ja) | 磁気共鳴イメージング装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040803 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040803 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060721 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060801 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061018 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20061121 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20061221 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |