JPH10332645A - 渦電流補正方法及び装置 - Google Patents
渦電流補正方法及び装置Info
- Publication number
- JPH10332645A JPH10332645A JP10040783A JP4078398A JPH10332645A JP H10332645 A JPH10332645 A JP H10332645A JP 10040783 A JP10040783 A JP 10040783A JP 4078398 A JP4078398 A JP 4078398A JP H10332645 A JPH10332645 A JP H10332645A
- Authority
- JP
- Japan
- Prior art keywords
- gradient
- eddy current
- planar imaging
- echo
- magnetic resonance
- 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 abstract description 36
- 238000009792 diffusion process Methods 0.000 claims abstract description 58
- 238000012937 correction Methods 0.000 claims abstract description 35
- 238000003384 imaging method Methods 0.000 claims abstract description 32
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 25
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 238000013480 data collection Methods 0.000 claims description 2
- 239000000470 constituent Substances 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 8
- 230000036962 time dependent Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 230000002829 reductive effect Effects 0.000 description 4
- 208000006011 Stroke Diseases 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 201000006474 Brain Ischemia Diseases 0.000 description 1
- 206010008120 Cerebral ischaemia Diseases 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 206010008118 cerebral infarction Diseases 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000002597 diffusion-weighted imaging Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
- 229940124597 therapeutic agent Drugs 0.000 description 1
- 230000001225 therapeutic effect Effects 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/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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56341—Diffusion imaging
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Signal Processing (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Vascular Medicine (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
(57)【要約】 （修正有）
【課題】 シフト、ずれ、収縮又は膨張及び信号強度損
失等の望ましくない渦電流の影響を大幅に減少させるこ
とできる渦電流補正方法及び装置を提供する。 【解決手段】 ＭＲシステム１０に用いられる渦電流補
正方法が提供される。この補正方法は、拡散加重用勾配
によって生じる渦電流に誘起される磁場の各々の成分と
関連した渦電流パラメータを導出する工程と、その各々
が上述のパラメータの関数であるような一組の補正項を
発生する工程とを含んでいる。この方法は、一組の理想
的なエコー・プラナ・イメージング勾配、並びに受信器
１２の位相及び周波数を修正して、渦電流に誘起される
磁場勾配及び渦電流に誘起される空間依存的でないＢ0
磁場をそれぞれ相殺する更なる工程を含んでいる。この
修正は２段階で行われる。
失等の望ましくない渦電流の影響を大幅に減少させるこ
とできる渦電流補正方法及び装置を提供する。 【解決手段】 ＭＲシステム１０に用いられる渦電流補
正方法が提供される。この補正方法は、拡散加重用勾配
によって生じる渦電流に誘起される磁場の各々の成分と
関連した渦電流パラメータを導出する工程と、その各々
が上述のパラメータの関数であるような一組の補正項を
発生する工程とを含んでいる。この方法は、一組の理想
的なエコー・プラナ・イメージング勾配、並びに受信器
１２の位相及び周波数を修正して、渦電流に誘起される
磁場勾配及び渦電流に誘起される空間依存的でないＢ0
磁場をそれぞれ相殺する更なる工程を含んでいる。この
修正は２段階で行われる。
Description
【０００１】
【産業上の利用分野】ここに開示されると共に請求され
る方法は一般的には、拡散加重（diffusion-weighted）
エコー・プラナ・イメージング（ＥＰＩ）として知られ
ている磁気共鳴（ＭＲ）イメージング手法に関するもの
である。より具体的には、本発明は、拡散加重ＥＰＩ手
法に従って収集されるＭＲ画像において、拡散加重用勾
配によって誘起される渦電流に起因する誤差を大幅に減
少させる方法に関する。
る方法は一般的には、拡散加重（diffusion-weighted）
エコー・プラナ・イメージング（ＥＰＩ）として知られ
ている磁気共鳴（ＭＲ）イメージング手法に関するもの
である。より具体的には、本発明は、拡散加重ＥＰＩ手
法に従って収集されるＭＲ画像において、拡散加重用勾
配によって誘起される渦電流に起因する誤差を大幅に減
少させる方法に関する。
【０００２】
【従来の技術】ＭＲイメージングの分野の業者には周知
のように、拡散加重画像は、例えば、１８０°リフォー
カシング（再収束）無線周波数（ＲＦ）パルスの前後に
配置された一対の拡散加重用勾配を用いて取得すること
ができる。拡散加重用勾配を用いると、相異なる拡散係
数を有しているスピンは、式Ｓ＝Ｓ0 ｅ-bD に従って相
異なる信号損失度を示す。この式において、Ｓ及びＳ0
はそれぞれ、拡散加重用勾配を伴う信号及び伴わない信
号であり、Ｄは、所与の組織についての拡散係数であ
り、ｂは、拡散加重用勾配の振幅の平方に比例する「ｂ
因子」として知られている。明らかに異なる拡散係数を
有している組織がＭＲ画像内で適切なコントラストを示
すようにするためには、拡散加重用勾配の振幅は、十分
に大きなｂ因子を確実に得るように十分強くなければな
らない。
のように、拡散加重画像は、例えば、１８０°リフォー
カシング（再収束）無線周波数（ＲＦ）パルスの前後に
配置された一対の拡散加重用勾配を用いて取得すること
ができる。拡散加重用勾配を用いると、相異なる拡散係
数を有しているスピンは、式Ｓ＝Ｓ0 ｅ-bD に従って相
異なる信号損失度を示す。この式において、Ｓ及びＳ0
はそれぞれ、拡散加重用勾配を伴う信号及び伴わない信
号であり、Ｄは、所与の組織についての拡散係数であ
り、ｂは、拡散加重用勾配の振幅の平方に比例する「ｂ
因子」として知られている。明らかに異なる拡散係数を
有している組織がＭＲ画像内で適切なコントラストを示
すようにするためには、拡散加重用勾配の振幅は、十分
に大きなｂ因子を確実に得るように十分強くなければな
らない。
【０００３】
【発明が解決しようとする課題】拡散加重ＭＲイメージ
ング手法は、脳の虚血の早期検出に特に有用である。脳
卒中の発症から僅か数時間以内に、拡散加重イメージン
グは優れたコントラストで虚血領域を強調表示すること
ができるが、他の手法では、遥かに遅れた時間に虚血を
検出するか、又は全く検出できないかのいずれかであ
る。早期の脳卒中の検出は、ＴＰＡ等の治療薬が、典型
的には数時間持続する比較的狭い「治療ウィンドウ」
（therapeutic window）の範囲内でしか有効でないの
で、特に重要である。有効な脳卒中検出を行うために、
６００ｓ／（ｍｍ）2 〜１０００ｓ／（ｍｍ）2 の間の
ｂ因子が頻繁に用いられるので、結果的に、大きな拡散
加重用勾配（例えば、２Ｇ／ｃｍ）が長時間（例えば、
３０ｍｓ）にわたって作用する。このような拡散加重用
勾配のため、パルス・シーケンスは患者の動き及び渦電
流の影響に対して極端に敏感になる。動きの問題は、Ｅ
ＰＩ手法を用いて効果的に除去され得るが、渦電流の問
題は、ＥＰＩでは一層深刻になり、画像シフト、形状の
歪み及び強度の低下を含めた劣化効果を招く。
ング手法は、脳の虚血の早期検出に特に有用である。脳
卒中の発症から僅か数時間以内に、拡散加重イメージン
グは優れたコントラストで虚血領域を強調表示すること
ができるが、他の手法では、遥かに遅れた時間に虚血を
検出するか、又は全く検出できないかのいずれかであ
る。早期の脳卒中の検出は、ＴＰＡ等の治療薬が、典型
的には数時間持続する比較的狭い「治療ウィンドウ」
（therapeutic window）の範囲内でしか有効でないの
で、特に重要である。有効な脳卒中検出を行うために、
６００ｓ／（ｍｍ）2 〜１０００ｓ／（ｍｍ）2 の間の
ｂ因子が頻繁に用いられるので、結果的に、大きな拡散
加重用勾配（例えば、２Ｇ／ｃｍ）が長時間（例えば、
３０ｍｓ）にわたって作用する。このような拡散加重用
勾配のため、パルス・シーケンスは患者の動き及び渦電
流の影響に対して極端に敏感になる。動きの問題は、Ｅ
ＰＩ手法を用いて効果的に除去され得るが、渦電流の問
題は、ＥＰＩでは一層深刻になり、画像シフト、形状の
歪み及び強度の低下を含めた劣化効果を招く。
【０００４】
【課題を解決するための手段】本発明は、拡散加重ＥＰ
Ｉパルス・シーケンスを発生してＭＲ画像を収集するよ
うに設けられているＭＲシステムに関連して用いられる
渦電流補正方法を提供する。このパルス・シーケンス
は、従来のＥＰＩパルス・シーケンス内に、１８０°リ
フォーカシングＲＦパルスの前後に配置された一対の拡
散加重用勾配と、他の成分とを含んでいる。この方法
は、拡散加重用勾配によって生ずる渦電流磁場成分にそ
れぞれ関係している一組のパラメータを導き出す工程を
含んでいる。これらのパラメータから一組の補正項が発
生され、これらの補正項に従って、一組の理想的なＥＰ
Ｉ勾配、並びに受信器の周波数及び位相がそれぞれ修正
される。修正されたＥＰＩ勾配、並びに修正された受信
器の周波数及び位相を上述のパルス・シーケンス内に含
めてＭＲデータを収集するので、渦電流磁場によって生
ずる収集されたデータ・セットの不完全性が実質的に減
少する。
Ｉパルス・シーケンスを発生してＭＲ画像を収集するよ
うに設けられているＭＲシステムに関連して用いられる
渦電流補正方法を提供する。このパルス・シーケンス
は、従来のＥＰＩパルス・シーケンス内に、１８０°リ
フォーカシングＲＦパルスの前後に配置された一対の拡
散加重用勾配と、他の成分とを含んでいる。この方法
は、拡散加重用勾配によって生ずる渦電流磁場成分にそ
れぞれ関係している一組のパラメータを導き出す工程を
含んでいる。これらのパラメータから一組の補正項が発
生され、これらの補正項に従って、一組の理想的なＥＰ
Ｉ勾配、並びに受信器の周波数及び位相がそれぞれ修正
される。修正されたＥＰＩ勾配、並びに修正された受信
器の周波数及び位相を上述のパルス・シーケンス内に含
めてＭＲデータを収集するので、渦電流磁場によって生
ずる収集されたデータ・セットの不完全性が実質的に減
少する。
【０００５】
【本発明の目的】本発明の目的は、拡散加重エコー・プ
ラナ・イメージングの改良された方法を提供することに
ある。もう１つの目的は、シフト、ずれ（shear）、収
縮（又は膨張）及び信号強度損失等の望ましくない渦電
流の影響を大幅に減少させる拡散加重ＥＰＩ走査のため
の方法を提供することにある。
ラナ・イメージングの改良された方法を提供することに
ある。もう１つの目的は、シフト、ずれ（shear）、収
縮（又は膨張）及び信号強度損失等の望ましくない渦電
流の影響を大幅に減少させる拡散加重ＥＰＩ走査のため
の方法を提供することにある。
【０００６】もう１つの目的は、ＥＰＩ走査に要求され
る読み出し勾配、位相エンコーディング勾配及びスライ
ス選択勾配、並びに受信器の周波数及び位相が、それぞ
れの理想的な波形から修正され、これにより、渦電流の
悪影響を減少させるような上述の形式の方法を提供する
ことにある。本発明のこれらの目的及びその他の目的
は、以下の記載及び図面からより容易に明らかになろ
う。
る読み出し勾配、位相エンコーディング勾配及びスライ
ス選択勾配、並びに受信器の周波数及び位相が、それぞ
れの理想的な波形から修正され、これにより、渦電流の
悪影響を減少させるような上述の形式の方法を提供する
ことにある。本発明のこれらの目的及びその他の目的
は、以下の記載及び図面からより容易に明らかになろ
う。
【０００７】
【実施例】図１を参照すると、同図には、本明細書に記
載されている拡散加重ＥＰＩ走査を実行するように動作
し得るＭＲシステム１０の基本的な構成要素が示されて
いる。システム１０は、ＲＦコイル１２と、イメージン
グされる被検体１６を収容している円筒マグネットのボ
ア（中孔）内に主磁場又は静磁場Ｂ0 を発生するマグネ
ット１４とを含んでいる。システム１０は更に、直交し
ているｘ基準軸、ｙ基準軸及びｚ基準軸に対してＧx 磁
場勾配、Ｇy 磁場勾配及びＧz 磁場勾配をそれぞれ発生
する勾配コイル１８、２０及び２２を含んでいる。図１
は、増幅器２４、２６及び２８によってそれぞれ駆動さ
れる勾配コイル１８、２０及び２２の各々を示してお
り、ＲＦコイル１２は、送受信増幅器３０に結合されて
いる。図１を更に詳細に参照すると、システム１０はパ
ルス・シーケンス制御器３２が設けられているものとし
て示されており、パルス・シーケンス制御器３２は、Ｒ
Ｆ増幅器及び勾配増幅器を制御するように動作し、これ
により、拡散加重ＥＰＩパルス・シーケンスを発生し
て、一組のＭＲデータを発生すると共に収集する。シス
テム１０は又、イメージング平面３６において取得され
る被検体１６の画像を構成するように、収集されたＭＲ
データを処理する画像処理電子回路３４を含んでいる。
載されている拡散加重ＥＰＩ走査を実行するように動作
し得るＭＲシステム１０の基本的な構成要素が示されて
いる。システム１０は、ＲＦコイル１２と、イメージン
グされる被検体１６を収容している円筒マグネットのボ
ア（中孔）内に主磁場又は静磁場Ｂ0 を発生するマグネ
ット１４とを含んでいる。システム１０は更に、直交し
ているｘ基準軸、ｙ基準軸及びｚ基準軸に対してＧx 磁
場勾配、Ｇy 磁場勾配及びＧz 磁場勾配をそれぞれ発生
する勾配コイル１８、２０及び２２を含んでいる。図１
は、増幅器２４、２６及び２８によってそれぞれ駆動さ
れる勾配コイル１８、２０及び２２の各々を示してお
り、ＲＦコイル１２は、送受信増幅器３０に結合されて
いる。図１を更に詳細に参照すると、システム１０はパ
ルス・シーケンス制御器３２が設けられているものとし
て示されており、パルス・シーケンス制御器３２は、Ｒ
Ｆ増幅器及び勾配増幅器を制御するように動作し、これ
により、拡散加重ＥＰＩパルス・シーケンスを発生し
て、一組のＭＲデータを発生すると共に収集する。シス
テム１０は又、イメージング平面３６において取得され
る被検体１６の画像を構成するように、収集されたＭＲ
データを処理する画像処理電子回路３４を含んでいる。
【０００８】上述のＭＲシステム１０のそれぞれの構成
要素の構成、作用及び相互関係は周知であり、１９９２
年９月２９日にMaier等に付与された米国特許第５，１
６１，６５６号等の従来技術に十分詳細に記載されてい
る。このような特許は、ＥＰＩ走査に特に指向してお
り、その教示はここに参照されるべきものである。図２
を参照すると、同図には、９０°ＲＦパルス３８と、１
８０°ＲＦパルス４０とを含んでいる拡散加重ＥＰＩパ
ルス・シーケンスが示されている。これらのＲＦパルス
は、ＲＦコイル１２によって送信することができ、空間
情報を含めてエンコードされ得るエコー信号４２を発生
する。エコー信号も又、コイル１２によって受信するこ
とができ、画像を再構成するのに用いられる。
要素の構成、作用及び相互関係は周知であり、１９９２
年９月２９日にMaier等に付与された米国特許第５，１
６１，６５６号等の従来技術に十分詳細に記載されてい
る。このような特許は、ＥＰＩ走査に特に指向してお
り、その教示はここに参照されるべきものである。図２
を参照すると、同図には、９０°ＲＦパルス３８と、１
８０°ＲＦパルス４０とを含んでいる拡散加重ＥＰＩパ
ルス・シーケンスが示されている。これらのＲＦパルス
は、ＲＦコイル１２によって送信することができ、空間
情報を含めてエンコードされ得るエコー信号４２を発生
する。エコー信号も又、コイル１２によって受信するこ
とができ、画像を再構成するのに用いられる。
【０００９】エコー・プラナ・イメージングに従ってエ
コー信号を空間的にエンコードするために、図２に示す
シーケンスは更に、読み出し勾配Ｇroと、位相エンコー
ディング勾配Ｇpeと、スライス選択勾配Ｇslとを含んで
いる。読み出し勾配Ｇroは、プリ・フェイジング（予備
位相合わせ）・パルス４４と、読み出しパルス４６とを
含んでいる。同様に、位相エンコーディング勾配Ｇ
peは、プリ・フェイジング・パルス４８と、位相エンコ
ーディング・パルス５０とを含んでいる。スライス選択
勾配Ｇslは、９０°ＲＦパルスに対するスライス選択パ
ルス５２と、１８０°ＲＦパルスに対するスライス選択
パルス５４と、更にスライス・リフォーカシングのため
のパルス５２ａとを含んでいる。
コー信号を空間的にエンコードするために、図２に示す
シーケンスは更に、読み出し勾配Ｇroと、位相エンコー
ディング勾配Ｇpeと、スライス選択勾配Ｇslとを含んで
いる。読み出し勾配Ｇroは、プリ・フェイジング（予備
位相合わせ）・パルス４４と、読み出しパルス４６とを
含んでいる。同様に、位相エンコーディング勾配Ｇ
peは、プリ・フェイジング・パルス４８と、位相エンコ
ーディング・パルス５０とを含んでいる。スライス選択
勾配Ｇslは、９０°ＲＦパルスに対するスライス選択パ
ルス５２と、１８０°ＲＦパルスに対するスライス選択
パルス５４と、更にスライス・リフォーカシングのため
のパルス５２ａとを含んでいる。
【００１０】図２を更に詳細に参照すると、同図には、
拡散加重ＥＰＩ走査に要求される拡散加重用勾配Ｇd が
示されている。拡散加重用勾配Ｇd は、１８０°ＲＦパ
ルスの両側に配置された２つの等価な台形パルス５６及
び５８を含んでいる。図１に示すイメージング平面３６
が基準軸のうちの１つに直交しており且つ他の２つの基
準軸に平行であるならば、読み出し勾配、位相エンコー
ディング勾配及びスライス選択勾配は、例えば、ｘ勾配
軸、ｙ勾配軸及びｚ勾配軸に沿ってそれぞれ位置すると
いうように、物理的基準軸のうちの１つに沿ってそれぞ
れ位置することになる。更に、このような取り決めを行
った場合には、Ｇx 勾配は読み出し勾配Ｇroを排他的に
含み、Ｇy 勾配は位相エンコーディング勾配Ｇpeを排他
的に含み、Ｇz 勾配はスライス選択勾配Ｇslを排他的に
含むことになる。勾配コイル１８、２０及び２２によっ
てそれぞれ発生されるＧx 勾配、Ｇy 勾配及びＧ z 勾配
を、以下では物理的勾配と呼ぶ。Ｇro勾配、Ｇpe勾配及
びＧsl勾配を、以下ではそれぞれの論理的勾配と呼ぶ。
拡散加重ＥＰＩ走査に要求される拡散加重用勾配Ｇd が
示されている。拡散加重用勾配Ｇd は、１８０°ＲＦパ
ルスの両側に配置された２つの等価な台形パルス５６及
び５８を含んでいる。図１に示すイメージング平面３６
が基準軸のうちの１つに直交しており且つ他の２つの基
準軸に平行であるならば、読み出し勾配、位相エンコー
ディング勾配及びスライス選択勾配は、例えば、ｘ勾配
軸、ｙ勾配軸及びｚ勾配軸に沿ってそれぞれ位置すると
いうように、物理的基準軸のうちの１つに沿ってそれぞ
れ位置することになる。更に、このような取り決めを行
った場合には、Ｇx 勾配は読み出し勾配Ｇroを排他的に
含み、Ｇy 勾配は位相エンコーディング勾配Ｇpeを排他
的に含み、Ｇz 勾配はスライス選択勾配Ｇslを排他的に
含むことになる。勾配コイル１８、２０及び２２によっ
てそれぞれ発生されるＧx 勾配、Ｇy 勾配及びＧ z 勾配
を、以下では物理的勾配と呼ぶ。Ｇro勾配、Ｇpe勾配及
びＧsl勾配を、以下ではそれぞれの論理的勾配と呼ぶ。
【００１１】当業界で周知のように、イメージング平面
３６が基準軸のうちの２つ又は３つすべてに関して非平
行で且つ非直交の関係にあるならば、そこで収集される
画像を斜方（oblique）画像と呼ぶ。斜方走査において
は、論理的勾配のうちの２つ又は３つすべての勾配が、
物理的勾配のうちの２つ又は３つすべての成分から形成
されている。拡散加重用勾配Ｇd も又、基準軸のうちの
１つの軸に沿って配向していてもよいし、又は代替的
に、より一般的な場合では、２つ若しくは３つの物理的
勾配成分を含んでいる。
３６が基準軸のうちの２つ又は３つすべてに関して非平
行で且つ非直交の関係にあるならば、そこで収集される
画像を斜方（oblique）画像と呼ぶ。斜方走査において
は、論理的勾配のうちの２つ又は３つすべての勾配が、
物理的勾配のうちの２つ又は３つすべての成分から形成
されている。拡散加重用勾配Ｇd も又、基準軸のうちの
１つの軸に沿って配向していてもよいし、又は代替的
に、より一般的な場合では、２つ若しくは３つの物理的
勾配成分を含んでいる。
【００１２】
【数１】
【００１３】渦電流によって誘起される上述の磁場は、
Ｂ0 磁場ばかりでなくＥＰＩ勾配磁場とも相互作用する
可能性があり、多数の方式で画質に悪影響を及ぼす。第
１に、ｂm0（ｔ）の存在下では、拡散加重ＥＰＩ画像は
位相エンコーディング方向に沿ってシフトする。ｂ
m0（ｔ）渦電流磁場は又、９０°スライス・プロファイ
ルと１８０°スライス・プロファイルとの間に相対的な
シフトを発生し、その結果、画像強度が低下する。第２
に、拡散加重用勾配軸から読み出し軸に向けての線形渦
電流磁場ｇmn（ｔ）は、位相エンコーディング方向に沿
った画像の歪み、即ち、ずれを発生する。渦電流磁場
は、渦電流を感受する軸が読み出し軸である限りでは、
軸上又は交差項のいずれかであり得る。第３に、拡散加
重用勾配から位相エンコーディング勾配に向けての線形
渦電流は、渦電流勾配と位相エンコーディング勾配との
間の相対的な極性に応じて、位相エンコーディング方向
に沿った画像の収縮又は膨張をもたらす。第４に、拡散
加重用勾配によって誘起され、スライス選択勾配上で作
用する線形渦電流磁場は、不完全なスライス・リフォー
カシングを導入し、その結果、画像の強度が低下する。
Ｂ0 磁場ばかりでなくＥＰＩ勾配磁場とも相互作用する
可能性があり、多数の方式で画質に悪影響を及ぼす。第
１に、ｂm0（ｔ）の存在下では、拡散加重ＥＰＩ画像は
位相エンコーディング方向に沿ってシフトする。ｂ
m0（ｔ）渦電流磁場は又、９０°スライス・プロファイ
ルと１８０°スライス・プロファイルとの間に相対的な
シフトを発生し、その結果、画像強度が低下する。第２
に、拡散加重用勾配軸から読み出し軸に向けての線形渦
電流磁場ｇmn（ｔ）は、位相エンコーディング方向に沿
った画像の歪み、即ち、ずれを発生する。渦電流磁場
は、渦電流を感受する軸が読み出し軸である限りでは、
軸上又は交差項のいずれかであり得る。第３に、拡散加
重用勾配から位相エンコーディング勾配に向けての線形
渦電流は、渦電流勾配と位相エンコーディング勾配との
間の相対的な極性に応じて、位相エンコーディング方向
に沿った画像の収縮又は膨張をもたらす。第４に、拡散
加重用勾配によって誘起され、スライス選択勾配上で作
用する線形渦電流磁場は、不完全なスライス・リフォー
カシングを導入し、その結果、画像の強度が低下する。
【００１４】拡散加重用勾配は比較的単純な形態を有し
ているので、拡散加重用勾配によって発生される空間的
に一定であると共に線形である渦電流磁場（ｂm0（ｔ）
及びｇmn（ｔ））は、解析的に算出することができる。
次いで、ｂm0（ｔ）及びｇmn（ｔ）が既知になれば、ア
ーティファクトが大幅に減少した画像を取得し得るよう
に、パルス・シーケンスを設計する際に渦電流の影響を
補償することができる。これが、本発明の渦電流補正方
法の背後にある基本的な考え方である。以下の議論で
は、Ｇro、Ｇpe及びＧslを、理想的な読み出し勾配、位
相エンコーディング勾配及びスライス選択勾配をそれぞ
れ表すために、即ち、ＥＰＩ走査に要求される論理的勾
配を表すために用いる。上述のように、項ｂm0（ｔ）
は、拡散勾配Ｇd によって誘起される時間依存的な渦電
流Ｂ0 磁場を表すために用いられ、項ｇmn（ｔ）は、や
はりＧd によって誘起される時間依存的な渦電流勾配磁
場のそれぞれの成分を表すために用いられる。最初の下
付き文字ｍは、渦電流ドナー（供与側）の軸（即ち、拡
散加重用勾配軸）を表し、第２の下付き文字ｎは、渦電
流アクセプタ（受容側）の軸（即ち、読み出し軸、位相
エンコーディング軸又はスライス選択軸）を指す。
ているので、拡散加重用勾配によって発生される空間的
に一定であると共に線形である渦電流磁場（ｂm0（ｔ）
及びｇmn（ｔ））は、解析的に算出することができる。
次いで、ｂm0（ｔ）及びｇmn（ｔ）が既知になれば、ア
ーティファクトが大幅に減少した画像を取得し得るよう
に、パルス・シーケンスを設計する際に渦電流の影響を
補償することができる。これが、本発明の渦電流補正方
法の背後にある基本的な考え方である。以下の議論で
は、Ｇro、Ｇpe及びＧslを、理想的な読み出し勾配、位
相エンコーディング勾配及びスライス選択勾配をそれぞ
れ表すために、即ち、ＥＰＩ走査に要求される論理的勾
配を表すために用いる。上述のように、項ｂm0（ｔ）
は、拡散勾配Ｇd によって誘起される時間依存的な渦電
流Ｂ0 磁場を表すために用いられ、項ｇmn（ｔ）は、や
はりＧd によって誘起される時間依存的な渦電流勾配磁
場のそれぞれの成分を表すために用いられる。最初の下
付き文字ｍは、渦電流ドナー（供与側）の軸（即ち、拡
散加重用勾配軸）を表し、第２の下付き文字ｎは、渦電
流アクセプタ（受容側）の軸（即ち、読み出し軸、位相
エンコーディング軸又はスライス選択軸）を指す。
【００１５】従来技術の教示によれば、例えば、Journa
l of Magnetic Resonance誌、第９０巻、第２６４頁〜
第２７８頁（１９９０年）のJehenson等の「ＮＭＲシス
テムにおけるパルス型磁場勾配によって誘起される渦電
流の影響を補償するための解析的方法」と題された刊行
物に述べられているように、渦電流磁場を以下のように
表すことができる。
l of Magnetic Resonance誌、第９０巻、第２６４頁〜
第２７８頁（１９９０年）のJehenson等の「ＮＭＲシス
テムにおけるパルス型磁場勾配によって誘起される渦電
流の影響を補償するための解析的方法」と題された刊行
物に述べられているように、渦電流磁場を以下のように
表すことができる。
【００１６】
【数２】
【００１７】式（１）において、ｍは印加される拡散加
重用勾配の物理的軸であり、時間変化する拡散勾配Ｇd
が、ｘ軸、ｙ軸又はｚ軸に沿って配向した成分を有して
いるか否かに応じて、ｘ、ｙ又はｚのいずれかであり得
る。ｎは、作用を受ける空間エンコーディング勾配の軸
であり、やはりｘ、ｙ又はｚのいずれかであり得る。量
αmnj 及びτmnj は、ＭＲシステムの隣接した構造体内
で拡散加重用勾配Ｇm（ｔ）によって誘起されるｊ番目
の線形渦電流成分のそれぞれ振幅及び時定数である。同
様に、αm0j 及びτm0j は、主磁場に作用するｊ番目の
Ｂ0 渦電流成分のそれぞれ振幅及び時定数である。過去
に、このような渦電流成分の各々が任意の所与のＭＲシ
ステムについて経験的に測定可能であることが発見され
ている。印加される所与の勾配Ｇm （ｔ）についての係
数αmnj 、αm0j 、τmnj 及びτ m0j のそれぞれの値を
決定する手法は従来技術において公知であり、例えば、
１９８６年１月３日にGlover等に付与された米国特許第
４，６９８，５９１号に記載されている。
重用勾配の物理的軸であり、時間変化する拡散勾配Ｇd
が、ｘ軸、ｙ軸又はｚ軸に沿って配向した成分を有して
いるか否かに応じて、ｘ、ｙ又はｚのいずれかであり得
る。ｎは、作用を受ける空間エンコーディング勾配の軸
であり、やはりｘ、ｙ又はｚのいずれかであり得る。量
αmnj 及びτmnj は、ＭＲシステムの隣接した構造体内
で拡散加重用勾配Ｇm（ｔ）によって誘起されるｊ番目
の線形渦電流成分のそれぞれ振幅及び時定数である。同
様に、αm0j 及びτm0j は、主磁場に作用するｊ番目の
Ｂ0 渦電流成分のそれぞれ振幅及び時定数である。過去
に、このような渦電流成分の各々が任意の所与のＭＲシ
ステムについて経験的に測定可能であることが発見され
ている。印加される所与の勾配Ｇm （ｔ）についての係
数αmnj 、αm0j 、τmnj 及びτ m0j のそれぞれの値を
決定する手法は従来技術において公知であり、例えば、
１９８６年１月３日にGlover等に付与された米国特許第
４，６９８，５９１号に記載されている。
【００１８】図３を参照すると、同図には、印加される
勾配Ｇm （ｔ）を構成している拡散加重用勾配パルス５
６及び５８が極く詳細に図示されている。これらのパル
スは、時刻０に開始して、時刻ｔ8 に終了している。図
３によれば、勾配Ｇm （ｔ）は、時刻ｔ1 と時刻ｔ2 と
の間及び時刻ｔ6 と時刻ｔ7 との間の時間周期δ2 には
最大振幅にある。印加される勾配は、−∞と時刻０との
間及び時刻ｔ8 と＋∞との間、並びにパルス５６とパル
ス５８とを分離している時刻ｔ3 と時刻ｔ5 との間の時
間周期δ3 にはゼロ振幅にある。勾配Ｇm （ｔ）は、時
刻０と時刻ｔ1との間及び時刻ｔ5 と時刻ｔ6 との間の
時間δ1 には、スルー・レートＲ（Ｒ＝Ｇd ／ｔ1 ）で
増大している。勾配Ｇm （ｔ）は、時刻ｔ2 と時刻ｔ3
との間及び時刻ｔ7 と時刻ｔ8 との間の時間周期δ1 に
は、スルー・レート−Ｒで減少している。
勾配Ｇm （ｔ）を構成している拡散加重用勾配パルス５
６及び５８が極く詳細に図示されている。これらのパル
スは、時刻０に開始して、時刻ｔ8 に終了している。図
３によれば、勾配Ｇm （ｔ）は、時刻ｔ1 と時刻ｔ2 と
の間及び時刻ｔ6 と時刻ｔ7 との間の時間周期δ2 には
最大振幅にある。印加される勾配は、−∞と時刻０との
間及び時刻ｔ8 と＋∞との間、並びにパルス５６とパル
ス５８とを分離している時刻ｔ3 と時刻ｔ5 との間の時
間周期δ3 にはゼロ振幅にある。勾配Ｇm （ｔ）は、時
刻０と時刻ｔ1との間及び時刻ｔ5 と時刻ｔ6 との間の
時間δ1 には、スルー・レートＲ（Ｒ＝Ｇd ／ｔ1 ）で
増大している。勾配Ｇm （ｔ）は、時刻ｔ2 と時刻ｔ3
との間及び時刻ｔ7 と時刻ｔ8 との間の時間周期δ1 に
は、スルー・レート−Ｒで減少している。
【００１９】図３に示す拡散加重用勾配Ｇm （ｔ）の時
間関係から、その時間導関数が以下のようになることが
直ちにわかる。 ｔ′∈（−∞，０］又は（ｔ1 ，ｔ2 ］又は（ｔ3 ，ｔ5 ］又は（ｔ6 ，ｔ 7 ］又は（ｔ8 ，＋∞］のときに、 ｄＧm （ｔ′）／ｄｔ′＝０ ｔ′∈（０，ｔ1 ］又は（ｔ5 ，ｔ6 ］のときに、 ｄＧm （ｔ′）／ｄｔ′＝Ｒ ｔ′∈（ｔ2 ，ｔ3 ］又は（ｔ7 ，ｔ8 ］のときに、 ｄＧm （ｔ′）／ｄｔ′＝−Ｒ 式（２） 式（１ａ）と式（２）とを組み合わせることにより、渦
電流の各成分の累積効果によって生ずる渦電流勾配ｇmn
（ｔ）は、以下のように表される。
間関係から、その時間導関数が以下のようになることが
直ちにわかる。 ｔ′∈（−∞，０］又は（ｔ1 ，ｔ2 ］又は（ｔ3 ，ｔ5 ］又は（ｔ6 ，ｔ 7 ］又は（ｔ8 ，＋∞］のときに、 ｄＧm （ｔ′）／ｄｔ′＝０ ｔ′∈（０，ｔ1 ］又は（ｔ5 ，ｔ6 ］のときに、 ｄＧm （ｔ′）／ｄｔ′＝Ｒ ｔ′∈（ｔ2 ，ｔ3 ］又は（ｔ7 ，ｔ8 ］のときに、 ｄＧm （ｔ′）／ｄｔ′＝−Ｒ 式（２） 式（１ａ）と式（２）とを組み合わせることにより、渦
電流の各成分の累積効果によって生ずる渦電流勾配ｇmn
（ｔ）は、以下のように表される。
【００２０】
【数３】
【００２１】式（３）に示す式の右辺では、単純化のた
めに下付き文字ｍｎが削除されていることを特記してお
く。
めに下付き文字ｍｎが削除されていることを特記してお
く。
【００２２】
【数４】
【００２３】Δｋmni ＝Δｋmn0 ＋Δｋ′mni である。図２及び図３を併せて参照すると、ＭＲデータ
の収集は、時刻ｔ9 の後になるまでは開始しないことが
わかる。従って、Δｋmn0 を収集前誤差と呼び、Δｋ′
mni を収集時誤差と呼ぶ。収集時誤差Δｋ′mni は時間
依存的であるが、収集前誤差Δｋmn0 は時間依存的では
なく、ｋ空間データの収集中は一定に留まっている。
の収集は、時刻ｔ9 の後になるまでは開始しないことが
わかる。従って、Δｋmn0 を収集前誤差と呼び、Δｋ′
mni を収集時誤差と呼ぶ。収集時誤差Δｋ′mni は時間
依存的であるが、収集前誤差Δｋmn0 は時間依存的では
なく、ｋ空間データの収集中は一定に留まっている。
【００２４】１８０°ＲＦパルスの位相反転効果のた
め、１８０°パルスの前の（即ち、０≦ｔ≦ｔ4 ）第１
の勾配ローブによって発生されるｋ空間誤差は、ｔ5 と
ｔ5 ＋ｔ4 との間の時間間隔中の第２の勾配ローブによ
って打ち消されることになる。従って、正味の収集前ｋ
空間誤差Δｋmn0 を決定するためには、ｔ4 からｔ9 ま
での第１の拡散加重用勾配ローブと、ｔ5 ＋ｔ4 からｔ
9 までの第２の勾配ローブとを考えるだけでよい。即
ち、
め、１８０°パルスの前の（即ち、０≦ｔ≦ｔ4 ）第１
の勾配ローブによって発生されるｋ空間誤差は、ｔ5 と
ｔ5 ＋ｔ4 との間の時間間隔中の第２の勾配ローブによ
って打ち消されることになる。従って、正味の収集前ｋ
空間誤差Δｋmn0 を決定するためには、ｔ4 からｔ9 ま
での第１の拡散加重用勾配ローブと、ｔ5 ＋ｔ4 からｔ
9 までの第２の勾配ローブとを考えるだけでよい。即
ち、
【００２５】
【数５】
【００２６】ここで、式（３）から、
【００２７】
【数６】
【００２８】である。式（５）の積分を実行後、以下の
式が得られる。
式が得られる。
【００２９】
【数７】
【００３０】上の各式でも、単純化のために右辺から下
付き文字ｍｎが削除されている。拡散加重用勾配がｘ軸
に沿って印加される成分を有しており、且つ読み出し勾
配、位相エンコーディング勾配及びスライス選択勾配が
ｘ軸、ｙ軸及びｚ軸に沿ってそれぞれ印加されているな
らば、これらの勾配にそれぞれ対応するｋ空間収集前誤
差を式（６）から容易に決定することができ、このｋ空
間収集前誤差は、以下の各式で構成される。
付き文字ｍｎが削除されている。拡散加重用勾配がｘ軸
に沿って印加される成分を有しており、且つ読み出し勾
配、位相エンコーディング勾配及びスライス選択勾配が
ｘ軸、ｙ軸及びｚ軸に沿ってそれぞれ印加されているな
らば、これらの勾配にそれぞれ対応するｋ空間収集前誤
差を式（６）から容易に決定することができ、このｋ空
間収集前誤差は、以下の各式で構成される。
【００３１】
【数８】
【００３２】式（７）〜式（９）において、Ｒ、α及び
τの各項は、式（１）及び式（２）の場合と同じ量をそ
れぞれ表している。式（７）〜式（９）は又、図３に示
されているタイミング・パラメータを含んでいる。ｍ及
びｎの他の値についての他のｋ空間収集前誤差が導出さ
れ得ることは容易に認められよう。これは、例えば、拡
散加重用勾配がｙ軸若しくはｚ軸に沿って印加されてい
るときに、又は拡散加重用勾配がこれらの物理的軸に沿
った成分を有しているときに起こる。
τの各項は、式（１）及び式（２）の場合と同じ量をそ
れぞれ表している。式（７）〜式（９）は又、図３に示
されているタイミング・パラメータを含んでいる。ｍ及
びｎの他の値についての他のｋ空間収集前誤差が導出さ
れ得ることは容易に認められよう。これは、例えば、拡
散加重用勾配がｙ軸若しくはｚ軸に沿って印加されてい
るときに、又は拡散加重用勾配がこれらの物理的軸に沿
った成分を有しているときに起こる。
【００３３】前述のように、静的なＢ0 磁場は、渦電流
によって誘起される磁場ｂm0（ｔ）による影響を受け
る。このような成分は、式（１ｂ）から、Ｇm （ｔ）並
びに係数α及びτの関数として容易に決定され得る。式
（１ｂ）を式（２）と組み合わせると、Ｂ0 磁場誤差に
ついての式を導出することができる。勾配誤差がｋ空間
誤差を発生する場合と同様に、ｂm0（ｔ）は、ＭＲ信号
に位相誤差φm0（ｔ）を発生する。この位相誤差も又、
時間依存的な収集前誤差φm00 と、時間依存的な収集時
誤差φm0′（ｔ）とに分割することができる。即ち、φ
m0（ｔ）＝φm00＋φm0′（ｔ）である。ｘ軸に沿って
印加される拡散加重用勾配については、収集前位相誤差
は以下の式によって与えられる。
によって誘起される磁場ｂm0（ｔ）による影響を受け
る。このような成分は、式（１ｂ）から、Ｇm （ｔ）並
びに係数α及びτの関数として容易に決定され得る。式
（１ｂ）を式（２）と組み合わせると、Ｂ0 磁場誤差に
ついての式を導出することができる。勾配誤差がｋ空間
誤差を発生する場合と同様に、ｂm0（ｔ）は、ＭＲ信号
に位相誤差φm0（ｔ）を発生する。この位相誤差も又、
時間依存的な収集前誤差φm00 と、時間依存的な収集時
誤差φm0′（ｔ）とに分割することができる。即ち、φ
m0（ｔ）＝φm00＋φm0′（ｔ）である。ｘ軸に沿って
印加される拡散加重用勾配については、収集前位相誤差
は以下の式によって与えられる。
【００３４】
【数９】
【００３５】同様に、物理的なｙ軸又はｚ軸に沿った拡
散加重用勾配に対応する収集前位相誤差も求めることが
できる。式（７）〜式（９）にそれぞれ記述されている
一定の収集前ｋ空間オフセットは、それぞれ読み出し方
向及び位相エンコーディング方向にあるプリ・フェイジ
ング・パルス４４及び４８、並びにスライス選択勾配方
向にあるリフォーカシング・パルス５２ａを修正するこ
とにより、厳密に除去することができる。例えば、プリ
・フェイジング・パルス又はリフォーカシング勾配の元
の振幅をＧppと置き、そのパルス幅をＴとし、その波形
は無次元の関数ψ（ｔ）によって定義されるものとする
と、このようなプリ・フェイジング・パルス又はリフォ
ーカシング勾配の振幅を式（１１）によって与えられる
新たな値に修正することができる。
散加重用勾配に対応する収集前位相誤差も求めることが
できる。式（７）〜式（９）にそれぞれ記述されている
一定の収集前ｋ空間オフセットは、それぞれ読み出し方
向及び位相エンコーディング方向にあるプリ・フェイジ
ング・パルス４４及び４８、並びにスライス選択勾配方
向にあるリフォーカシング・パルス５２ａを修正するこ
とにより、厳密に除去することができる。例えば、プリ
・フェイジング・パルス又はリフォーカシング勾配の元
の振幅をＧppと置き、そのパルス幅をＴとし、その波形
は無次元の関数ψ（ｔ）によって定義されるものとする
と、このようなプリ・フェイジング・パルス又はリフォ
ーカシング勾配の振幅を式（１１）によって与えられる
新たな値に修正することができる。
【００３６】
【数１０】
【００３７】ここでもやはり、下付き文字ｍ及びｎは、
拡散加重用勾配はどの軸に対して印加されているのか、
又、収集補正はどの軸を目標としているのかに応じて、
ｘ、ｙ又はｚのいずれかであり得る。代替的には、プリ
・フェイジング勾配パルス又はリフォーカシング勾配パ
ルスを、パルス幅Ｔを変更することにより修正すること
もできる。但しこのとき、修正後のパルスの勾配−時間
面積は、元のパルスの面積からｋmn0 ／γを減じたもの
に等しくする。
拡散加重用勾配はどの軸に対して印加されているのか、
又、収集補正はどの軸を目標としているのかに応じて、
ｘ、ｙ又はｚのいずれかであり得る。代替的には、プリ
・フェイジング勾配パルス又はリフォーカシング勾配パ
ルスを、パルス幅Ｔを変更することにより修正すること
もできる。但しこのとき、修正後のパルスの勾配−時間
面積は、元のパルスの面積からｋmn0 ／γを減じたもの
に等しくする。
【００３８】一定の位相オフセットφm00 は、ＭＲシス
テムの受信器の位相を調節することにより除去すること
もできるし、又は代替的には、ＭＲ画像の強度に影響を
与えるわけではないので放置しておくこともできる。ｘ
軸に沿って印加される拡散加重用勾配について、又、ｘ
軸、ｙ軸及びｚ軸に沿ってそれぞれ配向している読み出
し勾配、位相エンコーディング勾配及びスライス選択勾
配について、収集時誤差に関する式を式（１ａ）、式
（１ｂ）及び式（３）から以下のように導出することが
できる。
テムの受信器の位相を調節することにより除去すること
もできるし、又は代替的には、ＭＲ画像の強度に影響を
与えるわけではないので放置しておくこともできる。ｘ
軸に沿って印加される拡散加重用勾配について、又、ｘ
軸、ｙ軸及びｚ軸に沿ってそれぞれ配向している読み出
し勾配、位相エンコーディング勾配及びスライス選択勾
配について、収集時誤差に関する式を式（１ａ）、式
（１ｂ）及び式（３）から以下のように導出することが
できる。
【００３９】
【数１１】
【００４０】上述の収集時誤差の補正は、時間依存性の
ため、収集前誤差の補正よりも困難である。厳密な補正
は、論理的勾配及び受信器の周波数を、ＥＰＩ読み出し
トレインの全体にわたって時点ごとの（point-by-poin
t）形式で動的に調節することにより初めて可能にな
る。このような解決法は、パルス・シーケンスのプログ
ラミングに甚だしい複雑さを導入する。
ため、収集前誤差の補正よりも困難である。厳密な補正
は、論理的勾配及び受信器の周波数を、ＥＰＩ読み出し
トレインの全体にわたって時点ごとの（point-by-poin
t）形式で動的に調節することにより初めて可能にな
る。このような解決法は、パルス・シーケンスのプログ
ラミングに甚だしい複雑さを導入する。
【００４１】比較的単純でありながら有効であることが
判明した渦電流補正方法では、図４に示すように、一般
化された渦電流ｇmn（ｔ）についての多元指数関数的減
衰曲線６０を区分的（piece-wise）定数勾配６２によっ
て近似する。例えば、減衰曲線が、式（１３）に記述さ
れているような読み出し勾配の収集時誤差ｇxx（ｔ）を
表すとして、勾配近似６２は、セグメント６４で構成さ
れる。各々のセグメント６４は、読み出し勾配パルス４
６のうちに１つに対応しており、ｅｓｐに等しい時間長
を有している。ここで、ｅｓｐは、読み出し勾配のエコ
ー間間隔、即ち、図４に示されているす隣接した読み出
しパルス４６ａと読み出しパルス４６ｂとの間の間隔で
ある。各々のセグメント６４の振幅は、式（１３）に従
って、セグメント６４に対応している読み出しパルスの
中点ｔp において、ｇxx（ｔ）から決定される。従っ
て、パルス４６ａ、４６ｂ、４６ｃ及び４６ｄにそれぞ
れ対応しているセグメント６４の振幅は、図４に示すよ
うに、ｇ1 、ｇ2 、ｇ3 及びｇ4 を含んでいる。曲線６
０によって表される渦電流の影響を補償するために、各
々のセグメント６４の振幅を、その対応している読み出
し勾配パルスから減じることにより、理想的な読み出し
勾配の波形を修正する。修正後の読み出し勾配は破線６
６として図示されており、Ｇ−ｇ1 及び−Ｇ−ｇ2 等の
パルス振幅を有している。Ｇ及び−Ｇは、理想的な勾配
パルス４６ａ及び４６ｂのそれぞれの振幅である。
判明した渦電流補正方法では、図４に示すように、一般
化された渦電流ｇmn（ｔ）についての多元指数関数的減
衰曲線６０を区分的（piece-wise）定数勾配６２によっ
て近似する。例えば、減衰曲線が、式（１３）に記述さ
れているような読み出し勾配の収集時誤差ｇxx（ｔ）を
表すとして、勾配近似６２は、セグメント６４で構成さ
れる。各々のセグメント６４は、読み出し勾配パルス４
６のうちに１つに対応しており、ｅｓｐに等しい時間長
を有している。ここで、ｅｓｐは、読み出し勾配のエコ
ー間間隔、即ち、図４に示されているす隣接した読み出
しパルス４６ａと読み出しパルス４６ｂとの間の間隔で
ある。各々のセグメント６４の振幅は、式（１３）に従
って、セグメント６４に対応している読み出しパルスの
中点ｔp において、ｇxx（ｔ）から決定される。従っ
て、パルス４６ａ、４６ｂ、４６ｃ及び４６ｄにそれぞ
れ対応しているセグメント６４の振幅は、図４に示すよ
うに、ｇ1 、ｇ2 、ｇ3 及びｇ4 を含んでいる。曲線６
０によって表される渦電流の影響を補償するために、各
々のセグメント６４の振幅を、その対応している読み出
し勾配パルスから減じることにより、理想的な読み出し
勾配の波形を修正する。修正後の読み出し勾配は破線６
６として図示されており、Ｇ−ｇ1 及び−Ｇ−ｇ2 等の
パルス振幅を有している。Ｇ及び−Ｇは、理想的な勾配
パルス４６ａ及び４６ｂのそれぞれの振幅である。
【００４２】上述の区分的定数化の手順を、位相エンコ
ーディング勾配及びスライス選択勾配について実行する
ことは遥かに容易である。なぜならば、元の勾配が単極
であるか、又は存在していないかのいずれかであるから
である。図４を更に参照すると、同図には、上述の補償
方法及び式（１４）に従って理想的な位相エンコーディ
ング波形から修正された位相エンコーディング波形を含
んでいる曲線６８が示されている。同様に、曲線７０
は、上述の補償方法及び式（１５）に従って修正された
スライス選択波形を示している。
ーディング勾配及びスライス選択勾配について実行する
ことは遥かに容易である。なぜならば、元の勾配が単極
であるか、又は存在していないかのいずれかであるから
である。図４を更に参照すると、同図には、上述の補償
方法及び式（１４）に従って理想的な位相エンコーディ
ング波形から修正された位相エンコーディング波形を含
んでいる曲線６８が示されている。同様に、曲線７０
は、上述の補償方法及び式（１５）に従って修正された
スライス選択波形を示している。
【００４３】Ｂ0 収集時誤差を補償するためには、受信
器の位相は、異なるエコーについて調節されなければな
らず、又、受信器の周波数は、各々のエコー収集につい
て変更されなければならない。受信器の位相オフセット
は、以下の式を用いて各々の収集ウィンドウの中央にお
いて算出することができる。
器の位相は、異なるエコーについて調節されなければな
らず、又、受信器の周波数は、各々のエコー収集につい
て変更されなければならない。受信器の位相オフセット
は、以下の式を用いて各々の収集ウィンドウの中央にお
いて算出することができる。
【００４４】
【数１２】
【００４５】ここで、ｌ（エル）はｌ番目のエコーであ
り、ｔl はｌ番目のエコーの中央での時刻である。周波
数オフセットΔｆx0（ｔ）は、式（１２）及び式（１
７）からｔ＝ｔl として直接的に得ることができる。 Δｆx0（ｔ）＝（γ／２π）ｂx0（ｔ） 式（１７） 周波数の修正７２は、図４に示されている。
り、ｔl はｌ番目のエコーの中央での時刻である。周波
数オフセットΔｆx0（ｔ）は、式（１２）及び式（１
７）からｔ＝ｔl として直接的に得ることができる。 Δｆx0（ｔ）＝（γ／２π）ｂx0（ｔ） 式（１７） 周波数の修正７２は、図４に示されている。
【００４６】
【数１３】
【００４７】更に図５を参照すると、同図には、上述の
方法に従って改変された勾配補償値を含んでいる修正後
の位相エンコーディング波形７６、スライス選択波形７
８、及び周波数波形８０が示されている。これらの値
は、時刻ｔc において、式（１４）、式（１５）及び式
（１７）からそれぞれ算出された。上述の議論では、拡
散加重用勾配を物理的なｘ軸に沿って印加するものと仮
定していたが、区分的定数近似を用いるものにせよ又は
定数近似を用いるものにせよ、この渦電流補正方法は、
いかなる物理的軸における拡散加重用勾配にも対応する
ように容易に一般化することができる。
方法に従って改変された勾配補償値を含んでいる修正後
の位相エンコーディング波形７６、スライス選択波形７
８、及び周波数波形８０が示されている。これらの値
は、時刻ｔc において、式（１４）、式（１５）及び式
（１７）からそれぞれ算出された。上述の議論では、拡
散加重用勾配を物理的なｘ軸に沿って印加するものと仮
定していたが、区分的定数近似を用いるものにせよ又は
定数近似を用いるものにせよ、この渦電流補正方法は、
いかなる物理的軸における拡散加重用勾配にも対応する
ように容易に一般化することができる。
【００４８】用途によっては、拡散加重用勾配Ｇd を任
意の方向に印加して、組織拡散の異方性を研究すること
ができる。本発明の補正方法を拡張して、この一般的な
例に応用することができる。３つの論理的勾配軸ｒｏ、
ｐｅ及びｓｌに関して任意の方向に拡散勾配Ｇd が印加
されるとする。この基準フレームでは、Ｇd を３つの直
交成分Ｇd,ro、Ｇd,pe及びＧd,slによって表すことがで
きる。
意の方向に印加して、組織拡散の異方性を研究すること
ができる。本発明の補正方法を拡張して、この一般的な
例に応用することができる。３つの論理的勾配軸ｒｏ、
ｐｅ及びｓｌに関して任意の方向に拡散勾配Ｇd が印加
されるとする。この基準フレームでは、Ｇd を３つの直
交成分Ｇd,ro、Ｇd,pe及びＧd,slによって表すことがで
きる。
【００４９】
【数１４】
【００５０】論理的勾配成分Ｇd,ro、Ｇd,pe及びＧd,sl
は、回転行列Ｒによって物理的勾配に関係付けられる。
は、回転行列Ｒによって物理的勾配に関係付けられる。
【００５１】
【数１５】
【００５２】各々の独立した物理的勾配Ｇd,x 、Ｇd,y
及びＧd,z は、Ｂ0 磁場及び３つの線形勾配磁場、並び
に前述のより高次の空間依存的な磁場を形成するような
渦電流を発生する可能性がある。式（１）によれば、こ
れらの渦電流磁場は、誘導性勾配に比例している。この
ように、Ｇd,x 勾配、Ｇd,y 勾配及びＧd,z 勾配によっ
て発生される総合的なＢ0 磁場は、
及びＧd,z は、Ｂ0 磁場及び３つの線形勾配磁場、並び
に前述のより高次の空間依存的な磁場を形成するような
渦電流を発生する可能性がある。式（１）によれば、こ
れらの渦電流磁場は、誘導性勾配に比例している。この
ように、Ｇd,x 勾配、Ｇd,y 勾配及びＧd,z 勾配によっ
て発生される総合的なＢ0 磁場は、
【００５３】
【数１６】
【００５４】として得られる。ここで、ξx0、ξy0及び
ξz0は、以下のように定義される。
ξz0は、以下のように定義される。
【００５５】
【数１７】
【００５６】式（２１）において、ｍはｘ、ｙ又はｚで
あり得る。又、σm0j は、
あり得る。又、σm0j は、
【００５７】
【数１８】
【００５８】として与えられる。同様に、渦電流に誘起
される勾配磁場は、 ｇx ＝Ｇd,x ξxx＋Ｇd,y ξyx＋Ｇd,z ξzx ｇy ＝Ｇd,x ξxy＋Ｇd,y ξyy＋Ｇd,z ξzy ｇz ＝Ｇd,x ξxz＋Ｇd,y ξyz＋Ｇd,z ξzz として、又はマトリクス形式では、
される勾配磁場は、 ｇx ＝Ｇd,x ξxx＋Ｇd,y ξyx＋Ｇd,z ξzx ｇy ＝Ｇd,x ξxy＋Ｇd,y ξyy＋Ｇd,z ξzy ｇz ＝Ｇd,x ξxz＋Ｇd,y ξyz＋Ｇd,z ξzz として、又はマトリクス形式では、
【００５９】
【数１９】
【００６０】として得られる。ここで、式（２２）のξ
mnは以下のように定義される。
mnは以下のように定義される。
【００６１】
【数２０】
【００６２】式（２３）において、ｍ及びｎはｘ、ｙ又
はｚであり得る。又、σmnj は、
はｚであり得る。又、σmnj は、
【００６３】
【数２１】
【００６４】によって与えられる。式（１９）及び式
（２０）から、データ収集「中に」用いられるべき周波
数補償値を容易に導き出すことができる。
（２０）から、データ収集「中に」用いられるべき周波
数補償値を容易に導き出すことができる。
【００６５】
【数２２】
【００６６】ここで、γはｒａｄ／（ｇａｕｓｓ ｓｅ
ｃ）単位での磁気回転比である。データ収集「中の」勾
配補償パラメータの導出は、もう少し複雑なだけであ
る。この補正方法は論理的勾配レベルでの渦電流勾配を
補正するものであるので、式（２２）の物理的渦電流勾
配ｇx 、ｇy 及びｇz を論理的勾配に変換しなければな
らない。これは、回転行列の逆行列Ｒ-1を用いて達成さ
れ得る。ユニタリ回転行列についてはＲ-1はＲ+ （Ｒの
転置行列）と等しいことがわかっているので、以下の式
が得られる。
ｃ）単位での磁気回転比である。データ収集「中の」勾
配補償パラメータの導出は、もう少し複雑なだけであ
る。この補正方法は論理的勾配レベルでの渦電流勾配を
補正するものであるので、式（２２）の物理的渦電流勾
配ｇx 、ｇy 及びｇz を論理的勾配に変換しなければな
らない。これは、回転行列の逆行列Ｒ-1を用いて達成さ
れ得る。ユニタリ回転行列についてはＲ-1はＲ+ （Ｒの
転置行列）と等しいことがわかっているので、以下の式
が得られる。
【００６７】
【数２３】
【００６８】式（１９）、式（２２）及び式（２５）を
組み合わせると、エコー・トレイン収集中の勾配補償パ
ラメータについての変換式が、以下のように導き出され
る。
組み合わせると、エコー・トレイン収集中の勾配補償パ
ラメータについての変換式が、以下のように導き出され
る。
【００６９】
【数２４】
【００７０】全く同じ手順を用いて、ｋ空間内の「収集
前の」補償パラメータも導き出すことができる。
前の」補償パラメータも導き出すことができる。
【００７１】
【数２５】
【００７２】ここで、ηmnは以下のように定義される。
【００７３】
【数２６】
【００７４】同様に、位相補償パラメータは、
【００７５】
【数２７】
【００７６】として与えられる。ここで、ηm0はｎを０
で置き換えて式（２８）によって定義される。この時点
で、任意の拡散加重用勾配（即ち、任意の配向及び任意
の振幅）を収集前の補償パラメータ（式（２７）及び式
（２９））、並びに収集時の補償パラメータ（式（２
４）及び式（２６））に変換するための方程式が導き出
された。これらの誤差を、前述の方法を用いてパルス・
シーケンス内で補償することができる。
で置き換えて式（２８）によって定義される。この時点
で、任意の拡散加重用勾配（即ち、任意の配向及び任意
の振幅）を収集前の補償パラメータ（式（２７）及び式
（２９））、並びに収集時の補償パラメータ（式（２
４）及び式（２６））に変換するための方程式が導き出
された。これらの誤差を、前述の方法を用いてパルス・
シーケンス内で補償することができる。
【００７７】明らかに、上述の教示に照らして、本発明
の他の多くの改変及び変形が可能である。従って、開示
された発明の概念の範囲内で、本発明を特定的に記載さ
れたものとは他の方式で実行してもよいことを理解すべ
きである。
の他の多くの改変及び変形が可能である。従って、開示
された発明の概念の範囲内で、本発明を特定的に記載さ
れたものとは他の方式で実行してもよいことを理解すべ
きである。
【図１】本発明に従って拡散加重ＥＰＩ走査を実行する
ＭＲシステムの基本的な構成要素を示す概略図である。
ＭＲシステムの基本的な構成要素を示す概略図である。
【図２】拡散加重ＥＰＩ走査のための理想的な勾配波形
及びＲＦ波形を示すパルス・シーケンス・ダイアグラム
である。
及びＲＦ波形を示すパルス・シーケンス・ダイアグラム
である。
【図３】図２に示す拡散加重用勾配を、実施例の議論で
用いられているタイミング・パラメータと共に示す詳細
図である。
用いられているタイミング・パラメータと共に示す詳細
図である。
【図４】本発明の第１の実施例に従って修正された図２
に示す波形の部分を示す図である。
に示す波形の部分を示す図である。
【図５】本発明の第２の実施例に従って修正された図２
に示す波形の部分を示す図である。
に示す波形の部分を示す図である。
１０ ＭＲシステム １２ ＲＦコイル １４ マグネット １６ 被検体 １８、２０、２２ 勾配コイル ２４、２６、２８ 増幅器 ３０ 送受信増幅器 ３２ パルス・シーケンス制御器 ３４ 画像処理電子回路 ３６ イメージング平面 ３８ ９０°ＲＦパルス ４０ １８０°ＲＦパルス ４２ エコー信号 ４４、４８ プリ・フェイジング・パルス ４６、４６ａ、４６ｂ、４６ｃ、４６ｄ 読み出しパル
ス ５０ 位相エンコーディング・パルス ５２、５４ スライス選択パルス ５２ａ スライス・リフォーカシング・パルス ５６、５８ 拡散加重用勾配パルス ６０ 減衰曲線 ６２ 区分的定数勾配 ６４ セグメント ６６、７４ 修正された読み出し勾配 ６８、７６ 修正された位相エンコーディング波形 ７０、７８ 修正されたスライス選択波形 ７２、８０ 修正された周波数波形
ス ５０ 位相エンコーディング・パルス ５２、５４ スライス選択パルス ５２ａ スライス・リフォーカシング・パルス ５６、５８ 拡散加重用勾配パルス ６０ 減衰曲線 ６２ 区分的定数勾配 ６４ セグメント ６６、７４ 修正された読み出し勾配 ６８、７６ 修正された位相エンコーディング波形 ７０、７８ 修正されたスライス選択波形 ７２、８０ 修正された周波数波形
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ジョセフ・ケネス・メイア アメリカ合衆国、ウィスコンシン州、ミル ウォーキー、ノース・ウェイバリイ・プレ イス・ナンバー705、1121番 (72)発明者 グレン・レイノルズ・ハモンド アメリカ合衆国、ウィスコンシン州、ワウ ケシャ、ウェスト・ウッドローン・コー ト、4751番
Claims (11)
- 【請求項１】 拡散加重エコー・プラナ・イメージング
・パルス・シーケンスを発生して磁気共鳴画像を収集す
るように設けられている磁気共鳴システム（１０）にお
いて、前記パルス・シーケンスは、拡散加重用勾配を含
んでいると共に、多数のエコー・プラナ・イメージング
空間エンコーディング勾配を要求しており、 前記拡散加重用勾配により生ずる渦電流磁場にそれぞれ
関連している渦電流パラメータを導出する工程と、 当該補正項の各々が前記パラメータの関数であるような
一組の補正項を発生する工程と、 該補正項に従って一組の理想的なエコー・プラナ・イメ
ージング勾配を選択的に修正する工程と、 前記渦電流磁場により生ずる収集される前記磁気共鳴デ
ータのアーティファクトを実質的に減少させるために、
修正された前記エコー・プラナ・イメージング勾配（７
４、７６、７８）を、前記パルス・シーケンスに対して
要求される前記エコー・プラナ・イメージング論理的勾
配として採用する工程とを備えた渦電流補正方法。 - 【請求項２】 前記方法は、前記磁気共鳴システム（１
０）のデータ収集システムの受信器（１２）の位相及び
周波数を修正する工程と、該修正された位相及び周波数
（８０）を前記磁気共鳴データの収集時に採用する工程
を含んでいる請求項１に記載の渦電流補正方法。 - 【請求項３】 前記一組の補正項は、第１及び第２の補
正項のサブセットを含んでおり、 前記理想的なエコー・プラナ・イメージング勾配の各々
は、前記データの収集の前に発生する第１の勾配成分
と、前記データの収集中に発生する第２の勾配成分とを
含んでおり、 前記理想的な勾配のうちの所与の１つの勾配についての
前記第１（４４）及び第２（４６）の勾配成分は、前記
第１及び第２のサブセットの補正項によりそれぞれ修正
される請求項１に記載の渦電流補正方法。 - 【請求項４】 前記第１のサブセットの各々の補正項
は、定数のオフセット値を含んでおり、 前記第２のサブセットの各々の補正項は、前記渦電流パ
ラメータの時間変化する関数を含んでいる請求項３に記
載の渦電流補正方法。 - 【請求項５】 前記理想的なエコー・プラナ・イメージ
ング勾配は、相互に直交している読み出し勾配と、位相
エンコーディング勾配と、スライス選択勾配とを含んで
いる請求項４に記載の渦電流補正方法。 - 【請求項６】 前記読み出しパルスは、特定の時間間隔
を有しており、 前記磁気共鳴データの収集中の前記渦電流成分の累積的
な影響は、指数関数的な減衰曲線（６０）により表さ
れ、 前記所与の理想的な勾配の前記第２の勾配成分（４６）
は、一連のセグメント（６４）を含んでいる前記減衰曲
線の区分的近似により修正され、各々のセグメント（６
４）は、前記特定の間隔に等しい時間長と、前記第２の
サブセットの対応する補正項から算出される振幅とを有
している請求項５に記載の渦電流補正方法。 - 【請求項７】 前記所与の理想的な勾配は、前記データ
収集時間のｋ空間の中央において前記サブセットの対応
する補正項から算出される平均勾配値により修正される
請求項５に記載の渦電流補正方法。 - 【請求項８】 前記渦電流パラメータは、前記渦電流成
分のそれぞれの振幅及び時定数を含んでいる請求項１に
記載の渦電流補正方法。 - 【請求項９】 エコー・プラナ・イメージング・パルス
・シーケンスを発生して磁気共鳴画像を収集するように
設けられている磁気共鳴システム（１０）において、前
記パルス・シーケンスは、双極の勾配を含んでいると共
に、多数のエコー・プラナ・イメージング空間エンコー
ディング勾配を要求しており、 前記双極の勾配により生ずる渦電流磁場にそれぞれ関連
している渦電流パラメータを導出する工程と、 前記補正項に従って一組の理想的なエコー・プラナ・イ
メージング勾配を選択的に修正する工程と、 前記渦電流磁場により生ずる収集される前記磁気共鳴デ
ータのアーティファクトを実質的に減少させるために、
修正された前記エコー・プラナ・イメージング勾配（７
４、７６、７８）を、前記パルス・シーケンスに対して
要求される前記エコー・プラナ・イメージング論理的勾
配として採用する工程とを備えた渦電流補正方法。 - 【請求項１０】 拡散加重エコー・プラナ・イメージン
グ・パルス・シーケンスを発生して磁気共鳴画像を収集
するように設けられている磁気共鳴システム（１０）に
おいて、前記パルス・シーケンスは、拡散加重用勾配を
含んでいると共に、多数のエコー・プラナ・イメージン
グ空間エンコーディング勾配を要求しており、 前記拡散加重用勾配により生ずる渦電流磁場にそれぞれ
関連している渦電流パラメータを導出する手段と、 当該補正項の各々が前記パラメータの関数であるような
一組の補正項を発生する手段と、 該補正項に従って一組の理想的なエコー・プラナ・イメ
ージング勾配を選択的に修正する手段と、 前記渦電流磁場により生ずる収集される前記磁気共鳴デ
ータのアーティファクトを実質的に減少させるために、
修正された前記エコー・プラナ・イメージング勾配（７
４、７６、７８）を、前記パルス・シーケンスに対して
要求される前記エコー・プラナ・イメージング論理的勾
配として採用する手段とを備えた渦電流補正装置。 - 【請求項１１】 エコー・プラナ・イメージング・パル
ス・シーケンスを発生して磁気共鳴画像を収集するよう
に設けられている磁気共鳴システム（１０）において、
前記パルス・シーケンスは、双極の勾配を含んでいると
共に、多数のエコー・プラナ・イメージング空間エンコ
ーディング勾配を要求しており、 前記双極の勾配により生ずる渦電流磁場にそれぞれ関連
している渦電流パラメータを導出する手段と、 前記補正項に従って一組の理想的なエコー・プラナ・イ
メージング勾配を選択的に修正する手段と、 前記渦電流磁場により生ずる収集される前記磁気共鳴デ
ータのアーティファクトを実質的に減少させるために、
修正された前記エコー・プラナ・イメージング勾配（７
４、７６、７８）を、前記パルス・シーケンスに対して
要求される前記エコー・プラナ・イメージング論理的勾
配として採用する手段とを備えた渦電流補正装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US08/803,808 US5864233A (en) | 1997-02-22 | 1997-02-22 | Method to reduce eddy current effects in diffusion-weighted echo planar imaging |
| US08/803808 | 1997-02-22 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10332645A true JPH10332645A (ja) | 1998-12-18 |
| JP4121175B2 JP4121175B2 (ja) | 2008-07-23 |
Family
ID=25187477
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP04078398A Expired - Fee Related JP4121175B2 (ja) | 1997-02-22 | 1998-02-23 | 渦電流補正方法及び装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5864233A (ja) |
| JP (1) | JP4121175B2 (ja) |
| DE (1) | DE19807306B4 (ja) |
| IL (1) | IL123254A (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005270303A (ja) * | 2004-03-24 | 2005-10-06 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置および勾配磁場形成装置 |
| JP2007068996A (ja) * | 2005-09-02 | 2007-03-22 | Toshiba America Mri Inc | 拡散磁気共鳴イメージングの方法と装置 |
| JP2007252653A (ja) * | 2006-03-23 | 2007-10-04 | Ziosoft Inc | 拡散強調画像処理装置及び画像処理プログラム |
| JP2009172360A (ja) * | 2007-12-28 | 2009-08-06 | Toshiba Corp | 磁気共鳴イメージング装置および磁気共鳴イメージング装置の制御プログラム |
Families Citing this family (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6078176A (en) * | 1996-11-08 | 2000-06-20 | General Electric Company | Fast spin echo pulse sequence for diffusion weighted imaging |
| US6175236B1 (en) * | 1998-03-25 | 2001-01-16 | General Electric Company | Method for acquiring spatially and spectrally selective MR images |
| US6192263B1 (en) * | 1998-09-01 | 2001-02-20 | General Electric Company | Phase-sensitive inversion recovery method of MR imaging |
| WO2000031559A1 (en) * | 1998-11-25 | 2000-06-02 | Koninklijke Philips Electronics N.V. | Magnetic resonance method and device |
| US6323646B1 (en) | 1999-05-21 | 2001-11-27 | General Electric Company | Method and apparatus for producing diffusion weighted MR images |
| US6271666B1 (en) * | 1999-05-21 | 2001-08-07 | General Electric Company | Filtering of oscillatory B0 behavior in NMR devices |
| US6239599B1 (en) | 1999-05-21 | 2001-05-29 | General Electric Company | Method and apparatus for identifying errors in magnetic resonance imaging examinations |
| US6313630B1 (en) * | 1999-08-25 | 2001-11-06 | Ge Medical Systems Global Technology Company Llc | Modular gradient system for MRI system |
| US6329821B1 (en) | 1999-11-15 | 2001-12-11 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus to compensate for image artifacts caused by magnet vibration in an MR imaging system |
| US6380738B1 (en) | 1999-11-15 | 2002-04-30 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for reducing image artifacts caused by magnet vibration in an MR imaging system |
| WO2001038895A1 (en) | 1999-11-24 | 2001-05-31 | Board Of Regents, The University Of Texas System | Methods and systems for generating tractograms |
| US6448773B1 (en) | 2000-02-24 | 2002-09-10 | Toshiba America Mri, Inc. | Method and system for measuring and compensating for eddy currents induced during NMR imaging operations |
| US6933720B2 (en) * | 2001-12-11 | 2005-08-23 | Toshiba America Mri, Inc. | Sequence preconditioning for ultra-fast magnetic resonance imaging |
| US6670812B1 (en) * | 2002-09-13 | 2003-12-30 | Ge Medical Systems Global Technology, Llc | B-value calculation and correction using a linear segment gradient waveform model |
| DE10250922B4 (de) * | 2002-10-31 | 2008-01-31 | Siemens Ag | Verfahren zur Ermittlung des ADC-Koeffizienten in der Diffusionsgewichteten Magnetresonanz-Bildgebung bei Verwendung von Steady-State-Sequenzen |
| JP3930439B2 (ja) * | 2003-02-06 | 2007-06-13 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 渦電流補正方法および磁気共鳴撮影装置 |
| US7112964B2 (en) * | 2004-08-02 | 2006-09-26 | General Electric Company | Eddy current measurement and correction in magnetic resonance imaging systems with a static phantom |
| US7239140B1 (en) * | 2005-12-21 | 2007-07-03 | General Electric Company | Method and apparatus for vibration-related artifact reduction |
| DE102009003889B3 (de) * | 2009-01-02 | 2010-09-02 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Reduktion von Verzerrungen in der Diffusionsbildgebung |
| DE102010035539B4 (de) | 2010-08-26 | 2012-04-05 | Siemens Aktiengesellschaft | Verfahren zur Kompensation von Wirbelstromfeldern in Magnetresonanzaufnahmen und Magnetresonanzeinrichtung |
| JP5971683B2 (ja) | 2011-03-15 | 2016-08-17 | 東芝メディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
| US8942945B2 (en) | 2011-04-19 | 2015-01-27 | General Electric Company | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
| JP6139119B2 (ja) * | 2012-01-13 | 2017-05-31 | 東芝メディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
| CN104854471B (zh) * | 2012-12-12 | 2018-10-26 | 皇家飞利浦有限公司 | 针对磁共振弥散加权成像(dwi)的运动检测与校正方法 |
| US9709653B2 (en) | 2013-02-19 | 2017-07-18 | Toshiba Medical Systems Corporation | Mapping eddy current fields in MRI system |
| JP2022515825A (ja) | 2018-12-28 | 2022-02-22 | ハイパーファイン，インコーポレイテッド | 磁気共鳴撮像におけるヒステリシスの補正 |
| US20220123947A1 (en) | 2019-01-18 | 2022-04-21 | Zeu Technologies, Inc. | A Method for Generating Random Numbers in Blockchain Smart Contracts |
| CN109765512B (zh) * | 2019-01-18 | 2021-09-17 | 上海联影医疗科技股份有限公司 | 磁共振梯度系统及其涡流补偿方法及装置 |
| CN114556128A (zh) * | 2019-08-15 | 2022-05-27 | 海珀菲纳运营有限公司 | 涡流缓解系统和方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62189057A (ja) * | 1986-01-03 | 1987-08-18 | ゼネラル・エレクトリツク・カンパニイ | 磁界勾配の歪みを補償する方法 |
| JPH0531099A (ja) * | 1991-07-31 | 1993-02-09 | Toshiba Corp | 磁気共鳴イメージング装置 |
| JPH0647014A (ja) * | 1992-07-31 | 1994-02-22 | Yokogawa Medical Syst Ltd | Ｍｒによるイメージング方法 |
| JPH06261877A (ja) * | 1992-10-26 | 1994-09-20 | Philips Electron Nv | 磁気共鳴画像化における渦電流補償 |
| JPH08336506A (ja) * | 1995-06-09 | 1996-12-24 | Toshiba Corp | 磁気共鳴イメージング装置 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2580747B2 (ja) * | 1988-11-28 | 1997-02-12 | 株式会社島津製作所 | Ｍｒｉ装置の傾斜磁場印加方法 |
| JPH02264636A (ja) * | 1988-12-23 | 1990-10-29 | Picker Internatl Ltd | 核磁気共鳴方法 |
| US4965521A (en) * | 1989-08-11 | 1990-10-23 | Spectroscopy Imaging Systems | Method and apparatus for compensating eddy current effects in a magnetic resonance device having pulsed magnetic field gradients |
| US5270654A (en) * | 1991-07-05 | 1993-12-14 | Feinberg David A | Ultra-fast multi-section MRI using gradient and spin echo (grase) imaging |
| EP0538668B1 (en) * | 1991-10-25 | 1996-06-26 | The University Of Queensland | Correction of signal distortion in an NMR apparatus |
| DE69326202T2 (de) * | 1992-06-29 | 2000-04-13 | Koninklijke Philips Electronics N.V., Eindhoven | Verfahren und Gerät zur Kernresonanzabbildung |
| EP0595391B1 (en) * | 1992-10-26 | 2001-01-17 | Koninklijke Philips Electronics N.V. | Eddy current compensation in magnetic resonance imaging |
| US5378985A (en) * | 1993-07-15 | 1995-01-03 | General Electric Company | Fast spin echo prescan for MRI system |
| DE4325031C1 (de) * | 1993-07-26 | 1994-11-03 | Siemens Ag | Verfahren zur Erfassung von durch Gradienten verursachten Wirbelstrom-Magnetfeldern in einem Kernspinresonanzgerät |
| EP0644437B1 (en) * | 1993-09-16 | 1999-05-12 | Koninklijke Philips Electronics N.V. | Correction of read-gradient polarity in EPI and GRASE MRI |
| GB9513544D0 (en) * | 1995-07-04 | 1995-09-06 | Marconi Gec Ltd | Magnetic resonance methods and apparatus |
-
1997
- 1997-02-22 US US08/803,808 patent/US5864233A/en not_active Expired - Fee Related
-
1998
- 1998-02-10 IL IL12325498A patent/IL123254A/xx not_active IP Right Cessation
- 1998-02-20 DE DE19807306A patent/DE19807306B4/de not_active Expired - Fee Related
- 1998-02-23 JP JP04078398A patent/JP4121175B2/ja not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62189057A (ja) * | 1986-01-03 | 1987-08-18 | ゼネラル・エレクトリツク・カンパニイ | 磁界勾配の歪みを補償する方法 |
| JPH0531099A (ja) * | 1991-07-31 | 1993-02-09 | Toshiba Corp | 磁気共鳴イメージング装置 |
| JPH0647014A (ja) * | 1992-07-31 | 1994-02-22 | Yokogawa Medical Syst Ltd | Ｍｒによるイメージング方法 |
| JPH06261877A (ja) * | 1992-10-26 | 1994-09-20 | Philips Electron Nv | 磁気共鳴画像化における渦電流補償 |
| JPH08336506A (ja) * | 1995-06-09 | 1996-12-24 | Toshiba Corp | 磁気共鳴イメージング装置 |
Non-Patent Citations (1)
| Title |
|---|
| P.JEHENSON ET AL.: "Analytical Method for the Compensation of Eddy-Current Effects Induced by Pulsed Magnetic Field Grad", JOURNAL OF MAGNETIC RESONANCE, vol. 90, JPN6007008858, 1990, pages 264 - 278, XP023961036, ISSN: 0000931323, DOI: 10.1016/0022-2364(90)90133-T * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005270303A (ja) * | 2004-03-24 | 2005-10-06 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置および勾配磁場形成装置 |
| JP2007068996A (ja) * | 2005-09-02 | 2007-03-22 | Toshiba America Mri Inc | 拡散磁気共鳴イメージングの方法と装置 |
| JP2007252653A (ja) * | 2006-03-23 | 2007-10-04 | Ziosoft Inc | 拡散強調画像処理装置及び画像処理プログラム |
| US7974451B2 (en) | 2006-03-23 | 2011-07-05 | Ziosoft, Inc. | Diffusion weighted image processing apparatus and recording medium for recording image analysis program |
| JP2009172360A (ja) * | 2007-12-28 | 2009-08-06 | Toshiba Corp | 磁気共鳴イメージング装置および磁気共鳴イメージング装置の制御プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4121175B2 (ja) | 2008-07-23 |
| IL123254A (en) | 2000-10-31 |
| DE19807306A1 (de) | 1998-08-27 |
| IL123254A0 (en) | 1998-09-24 |
| US5864233A (en) | 1999-01-26 |
| DE19807306B4 (de) | 2010-10-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4121175B2 (ja) | 渦電流補正方法及び装置 | |
| EP0634664B1 (en) | Fast spin echo prescan for MRI system | |
| US5923168A (en) | Correction of artifacts caused by Maxwell terms in slice offset echo planar imaging | |
| US5942897A (en) | Magnetic resonance imaging apparatus | |
| JP3529446B2 (ja) | Ｅｐｉ及びｇｒａｓｅ ｍｒｉにおける読み出し傾斜磁界極性の補正方法 | |
| US5825185A (en) | Method for magnetic resonance spin echo scan calibration and reconstruction | |
| US6160397A (en) | Fast spin echo prescan for magnetic resonance imaging systems | |
| EP0735379B1 (en) | MRI apparatus | |
| EP1089087B1 (en) | Reduction of ghost artifacts due to a Maxwell term phase error caused by a data acquisition read gradient in MRI | |
| US6664787B2 (en) | Method, apparatus and recording medium for measuring and correcting phase error in spins in a phase axis direction | |
| US6472872B1 (en) | Real-time shimming of polarizing field in magnetic resonance system | |
| US6815952B1 (en) | Magnetic resonance diffusion imaging with eddy-current compensation | |
| JPH0856928A (ja) | 表面コイル配設により核磁化分布を決定するｍｒ法 | |
| JP4519967B2 (ja) | 磁気共鳴イメージング・システム内の残留磁気を補償する方法及び磁気共鳴イメージング・システム | |
| US6836113B2 (en) | Measurement and correction of gradient-induced cross-term magnetic fields in an EPI sequence | |
| EP1083438A2 (en) | MR imaging method and MRI apparatus | |
| US6265872B1 (en) | Diffusion detection by means of magnetic resonance | |
| Caprihan et al. | Flow-velocity imaging from linear regression of phase images with techniques for reducing eddy-current effects | |
| US6528996B1 (en) | Diffusion-weighted imaging method and apparatus for fast pulse sequence with MPG pulses | |
| JP2000157509A (ja) | Ｍｒｉシステムの較正方法、ｍｒｉシステム及びリセット勾配発生方法 | |
| US6586935B1 (en) | Magnetic resonance image artifact correction using navigator echo information | |
| JPH03264046A (ja) | 核磁気共鳴映像法及び装置 | |
| US5157330A (en) | Method and apparatus for compensating magnetic field inhomogeneity artifact in MRI | |
| US6229309B1 (en) | MR method | |
| US4714884A (en) | Method of eliminating effects of spurious NMR signals caused by imperfect 180 degree RF pulses |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050221 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071127 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080125 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080219 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080304 |
|
| 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: 20080401 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080428 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110509 Year of fee payment: 3 |
|
| 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 |