JPH0351176B2 - - Google Patents
Info
- Publication number
- JPH0351176B2 JPH0351176B2 JP61270030A JP27003086A JPH0351176B2 JP H0351176 B2 JPH0351176 B2 JP H0351176B2 JP 61270030 A JP61270030 A JP 61270030A JP 27003086 A JP27003086 A JP 27003086A JP H0351176 B2 JPH0351176 B2 JP H0351176B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- gradient
- region
- phase
- encoding
- 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 - Lifetime
Links
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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【発明の詳細な説明】
発明の背景
この発明は磁気共鳴MR方法に関する。更に特
定すれば、この発明は磁界勾配流れ補償パルスを
適当に選ぶことにより、勾配の不規則な位相の累
積が原因で、流れる原子核から生ずる像の人為効
果（アーテイフアクト）を減少する方法に関す
る。
定すれば、この発明は磁界勾配流れ補償パルスを
適当に選ぶことにより、勾配の不規則な位相の累
積が原因で、流れる原子核から生ずる像の人為効
果（アーテイフアクト）を減少する方法に関す
る。
従来、磁気共鳴現象は、化学的な組成の構造を
解析する為に、構造化学者によつて高解像度の磁
気共鳴分光装置で用いられて来た。最近になつ
て、MRは解剖学的な部分の作像並びに非侵入型
分光解析を生体内で行なう様な用途を持つ医学用
の診断様式が開発された。現在ではよく知られて
いるが、均質な分極磁界の中に配置された患者の
様なサンプル物体をラーモア周波数の無線周波
RFエネルギで照射することにより、この物体内
にMR現象を励振することが出来る。医学用の診
断の用途では、これは検査しようとする患者を円
筒形を持つRFコイルの磁界の中に位置ぎめし、
RFコイルをRF電力増幅器で付勢することによつ
て行なうのが典型的である。RF励振が中断した
時、同じRFコイル又は異なるRFコイルを用い
て、RFコイルの磁界の中に横たわつている患者
の容積から出て来るMR信号、大抵はスピンエコ
ーを検出する。完全なMR走査の過程で、複数個
のMR信号を観測するのが典型的である。こうい
う信号を使つて、検査する物体のMR作像情報又
はそれに関する分光学的な情報を取出す。
解析する為に、構造化学者によつて高解像度の磁
気共鳴分光装置で用いられて来た。最近になつ
て、MRは解剖学的な部分の作像並びに非侵入型
分光解析を生体内で行なう様な用途を持つ医学用
の診断様式が開発された。現在ではよく知られて
いるが、均質な分極磁界の中に配置された患者の
様なサンプル物体をラーモア周波数の無線周波
RFエネルギで照射することにより、この物体内
にMR現象を励振することが出来る。医学用の診
断の用途では、これは検査しようとする患者を円
筒形を持つRFコイルの磁界の中に位置ぎめし、
RFコイルをRF電力増幅器で付勢することによつ
て行なうのが典型的である。RF励振が中断した
時、同じRFコイル又は異なるRFコイルを用い
て、RFコイルの磁界の中に横たわつている患者
の容積から出て来るMR信号、大抵はスピンエコ
ーを検出する。完全なMR走査の過程で、複数個
のMR信号を観測するのが典型的である。こうい
う信号を使つて、検査する物体のMR作像情報又
はそれに関する分光学的な情報を取出す。
典型的な検査では、空間情報を信号に符号化す
る為に使われるパルス形の直線的な磁界勾配の存
在の下にMR信号を観測するのが普通である。検
査を受ける物体の予定の領域にある核スピンを励
振する為に選択的なRFパルス（例えば90゜及び
180゜パルス）と共にパルス形磁界勾配も用いる。
MR検査の過程で、普通のデカルト座標系の夫々
ｘ，ｙ及びｚ方向のパルス形磁界勾配Gx，Gy及
びGzを印加するのが望ましい場合が多い。然し、
実際には、磁界勾配パルスを印加する方向がどん
な形でも制限されず、必要に応じて任意の方向を
選ぶことが出来ることが理解されよう。
る為に使われるパルス形の直線的な磁界勾配の存
在の下にMR信号を観測するのが普通である。検
査を受ける物体の予定の領域にある核スピンを励
振する為に選択的なRFパルス（例えば90゜及び
180゜パルス）と共にパルス形磁界勾配も用いる。
MR検査の過程で、普通のデカルト座標系の夫々
ｘ，ｙ及びｚ方向のパルス形磁界勾配Gx，Gy及
びGzを印加するのが望ましい場合が多い。然し、
実際には、磁界勾配パルスを印加する方向がどん
な形でも制限されず、必要に応じて任意の方向を
選ぶことが出来ることが理解されよう。
血流による人為効果が、腹部及び胸部のMR作
像に於ける画像の重大な欠陥の１つであることが
判つた。こういう人為効果は（スピン捩れ形の様
なフーリエ変換作像方法を用いる時）位相符号化
方向の高い空間周波数の不規則な縞となつて現わ
れ、大動脈や、心臓や、大容積の血流が起るその
他の領域から出て来る。こういう縞が解剖学的な
表示をぼかす場合が多く、画像の一部分を使いも
のにならなくすることがある。
像に於ける画像の重大な欠陥の１つであることが
判つた。こういう人為効果は（スピン捩れ形の様
なフーリエ変換作像方法を用いる時）位相符号化
方向の高い空間周波数の不規則な縞となつて現わ
れ、大動脈や、心臓や、大容積の血流が起るその
他の領域から出て来る。こういう縞が解剖学的な
表示をぼかす場合が多く、画像の一部分を使いも
のにならなくすることがある。
この発明では、この様な縞の人為効果が、ビユ
ー毎に相関性を持たない、測定されたスピンエコ
ー信号の振幅及び位相の変動に由ることが判つ
た。こういう信号の変動の元を辿れば、流速の変
動と共に、パルス順序に応じて、１）選択的なパ
ルスの間の飛行時間の変動、２）選択的な
180゜RFパルスによる位相歪み及び３）勾配内の
移動するスピンの位相の累積の内の１つ又は更に
多くのメカニズムにたどりつく。
ー毎に相関性を持たない、測定されたスピンエコ
ー信号の振幅及び位相の変動に由ることが判つ
た。こういう信号の変動の元を辿れば、流速の変
動と共に、パルス順序に応じて、１）選択的なパ
ルスの間の飛行時間の変動、２）選択的な
180゜RFパルスによる位相歪み及び３）勾配内の
移動するスピンの位相の累積の内の１つ又は更に
多くのメカニズムにたどりつく。
この発明の主な目的は、勾配内の移動するスピ
ンの位相累積誤差を原因とする人為効果を減少す
る方法を提供することである。
ンの位相累積誤差を原因とする人為効果を減少す
る方法を提供することである。
発明の要約
不動の核スピン及び流れる核スピンを持つ物体
の領域から収集された磁気共鳴信号の不規則な位
相の累積を補償する方法を提供する。この方法
は、領域をラーモア周波数のRF励振パルスで照
射してMR信号を発生し、調節自在な振幅及び持
続時間のパラメータを持つパルス形磁界勾配を前
記領域の少なくとも１つの次元の軸線に沿つて印
加することを夫々含む様な複数個のビユーを求め
ることを含む。流れるスピンの流速のビユー毎の
変動による不規則な位相の累積が、振幅及び磁束
時間のパラメータの内の少なくとも一方を変更し
て、磁界勾配の少なくとも１つの高次モーメント
をゼロにすることによつて補償される。
の領域から収集された磁気共鳴信号の不規則な位
相の累積を補償する方法を提供する。この方法
は、領域をラーモア周波数のRF励振パルスで照
射してMR信号を発生し、調節自在な振幅及び持
続時間のパラメータを持つパルス形磁界勾配を前
記領域の少なくとも１つの次元の軸線に沿つて印
加することを夫々含む様な複数個のビユーを求め
ることを含む。流れるスピンの流速のビユー毎の
変動による不規則な位相の累積が、振幅及び磁束
時間のパラメータの内の少なくとも一方を変更し
て、磁界勾配の少なくとも１つの高次モーメント
をゼロにすることによつて補償される。
この発明に特有と考えられる特徴は特許請求の
範囲に記載してあるが、この発明自体、その構成
と作用並びにその他の目的及び利点は、以下図面
について説明する所から最もよく理解されよう。
範囲に記載してあるが、この発明自体、その構成
と作用並びにその他の目的及び利点は、以下図面
について説明する所から最もよく理解されよう。
発明の詳しい記載
勾配の中で移動するスピンの位相の累積がどの
様にして起るかを直感的に理解するには、奇数個
の陽子又は中性子を持つ原子核で核磁気共鳴現象
が起ることを最初に考えるのがよいと思われる。
陽子及び中性子のスピンの為、各々の原子核が正
味の磁気モーメントを持ち、こういう原子核で構
成されたサンプルが均質な静磁界B0の中に配置
された時、多数の核磁気モーメントが磁界と整合
して、磁界の方向に正味の巨視的な磁化Ｍを発生
する。磁界Boの作用により、磁気モーメントが、
印加磁界の強さと原子核の特性に関係する周波数
で、磁界の軸線の周りの歳差運動をする。歳差運
動の角周波数ωはラーモア周波数とも呼ばれる
が、ラーモア方程式ω＝γBで表わされる。ここ
でγは磁気回転比であつて、各々のMRアイソト
ープにとつて一定であり、Ｂは核スピンに作用す
る磁界（Boを含む）である。従つて、共鳴周波
数が、サンプルをその中においた磁界の強さに関
係することが理解されよう。
様にして起るかを直感的に理解するには、奇数個
の陽子又は中性子を持つ原子核で核磁気共鳴現象
が起ることを最初に考えるのがよいと思われる。
陽子及び中性子のスピンの為、各々の原子核が正
味の磁気モーメントを持ち、こういう原子核で構
成されたサンプルが均質な静磁界B0の中に配置
された時、多数の核磁気モーメントが磁界と整合
して、磁界の方向に正味の巨視的な磁化Ｍを発生
する。磁界Boの作用により、磁気モーメントが、
印加磁界の強さと原子核の特性に関係する周波数
で、磁界の軸線の周りの歳差運動をする。歳差運
動の角周波数ωはラーモア周波数とも呼ばれる
が、ラーモア方程式ω＝γBで表わされる。ここ
でγは磁気回転比であつて、各々のMRアイソト
ープにとつて一定であり、Ｂは核スピンに作用す
る磁界（Boを含む）である。従つて、共鳴周波
数が、サンプルをその中においた磁界の強さに関
係することが理解されよう。
第１図にはｘ及びｙ軸によつて定められた横平
面内の線１０で示す磁界勾配Gxが示されている。
Gx勾配がラーモア方程式に従つて、歳差運動の
周波数又は位相とｘ軸に沿つた位置の間の一意的
な関係を定める。つまり、核スピン１２は所定の
歳差運動周波数f1を持ち、これは（他の条件が変
わらないと仮定すれば）不動のスピンでは時間に
対して一定のまゝである。然し、（例えば血流に
より）スピンが矢印１４で示す様に、磁界勾配の
方向に沿つて移動すると、MR実験の過程で１２
から１６に移動したスピンは、その移動の間、連
続的に変化する歳差運動周波数を持ち、勾配に沿
つた新しい位置で周波数f2で終る。所定の期間に
わたつてスピンに累積した位相はこの期間にわた
るスピン歳差運動周波数の時間積分に関係する。
従つて、この明細書では「位相」及び「周波数」
と云う言葉は互換性を以て用いている。上の説明
はｘ軸の場合であるが、周波数、従つて位相に対
するこういう影響は、勾配を印加したあらゆる方
向に起ることが理解されよう。
面内の線１０で示す磁界勾配Gxが示されている。
Gx勾配がラーモア方程式に従つて、歳差運動の
周波数又は位相とｘ軸に沿つた位置の間の一意的
な関係を定める。つまり、核スピン１２は所定の
歳差運動周波数f1を持ち、これは（他の条件が変
わらないと仮定すれば）不動のスピンでは時間に
対して一定のまゝである。然し、（例えば血流に
より）スピンが矢印１４で示す様に、磁界勾配の
方向に沿つて移動すると、MR実験の過程で１２
から１６に移動したスピンは、その移動の間、連
続的に変化する歳差運動周波数を持ち、勾配に沿
つた新しい位置で周波数f2で終る。所定の期間に
わたつてスピンに累積した位相はこの期間にわた
るスピン歳差運動周波数の時間積分に関係する。
従つて、この明細書では「位相」及び「周波数」
と云う言葉は互換性を以て用いている。上の説明
はｘ軸の場合であるが、周波数、従つて位相に対
するこういう影響は、勾配を印加したあらゆる方
向に起ることが理解されよう。
上の説明から、磁界勾配の方向に沿つて動くス
ピンは、速度に正比例する位相の累積を持つと或
る間考えられていた。関心のある読者は、J.
Phys.E.：Sci.Instr誌（1978年）第281頁所載のJ.
R.シンガーの輪文「NMRによる拡散と流れの測
定及びスピン位相グラフの案内を参照されたい。
然し、この発明では、作像パルス順序内で速度が
ビユー毎に変動する場合、位相も同じ様に変動
し、画像に縞の人為効果となつて寄与することを
発見した。
ピンは、速度に正比例する位相の累積を持つと或
る間考えられていた。関心のある読者は、J.
Phys.E.：Sci.Instr誌（1978年）第281頁所載のJ.
R.シンガーの輪文「NMRによる拡散と流れの測
定及びスピン位相グラフの案内を参照されたい。
然し、この発明では、作像パルス順序内で速度が
ビユー毎に変動する場合、位相も同じ様に変動
し、画像に縞の人為効果となつて寄与することを
発見した。
スピン戻れ形とよく呼ばれるフーリエ変換作像
方法の場合について、この発明の好ましい実施例
を説明するが、この発明はこの場合に制限されな
い。一般的に、この発明は磁界勾配の中で動くス
ピンを持つ物体に関する情報を収集する任意の
MR方法に用いることが出来る。これに制約する
つもりはないが、こういう方法の例としては、２
次元及び３次元形のフーリエ変換（例えば米国特
許第4431968号参照）、多重角度投影再生作像、並
びにMR分光法（例えば係属中の米国特許出願通
し番号第626941号に記載されている）がある。何
れの場合も、以下の説明から判るが、検査される
物体の少なくとも１つの軸線に沿つた勾配パルス
の分布をこの発明に従つて変更して、磁界勾配の
少なくとも若干の高次モーメントをゼロに等しく
する。
方法の場合について、この発明の好ましい実施例
を説明するが、この発明はこの場合に制限されな
い。一般的に、この発明は磁界勾配の中で動くス
ピンを持つ物体に関する情報を収集する任意の
MR方法に用いることが出来る。これに制約する
つもりはないが、こういう方法の例としては、２
次元及び３次元形のフーリエ変換（例えば米国特
許第4431968号参照）、多重角度投影再生作像、並
びにMR分光法（例えば係属中の米国特許出願通
し番号第626941号に記載されている）がある。何
れの場合も、以下の説明から判るが、検査される
物体の少なくとも１つの軸線に沿つた勾配パルス
の分布をこの発明に従つて変更して、磁界勾配の
少なくとも若干の高次モーメントをゼロに等しく
する。
第２図は普通の２次元スピン捩れ形作像順序の
１つのビユーを示すものであるが、このビユー
で、横軸に示した期間１に、正のGz勾配パルス
が印加されることが判る。Cz勾配の方向はデカル
ト座標系のｚ軸の正の方向に任意に選ぶが、これ
は磁界B0の方向と一致する。更に期間１に、Gz
勾配の存在の下に選択性90゜RFパルスを印加し
て、検査する物体の平面状容積（図に示してな
い）内にある核スピンを励振する。選ばれるスラ
イスの厚さがGz勾配の振幅と選択性90゜RFパルス
の周波数成分とによつて決定される。この発明の
好ましい実施例では、RFパルスがsinc関数（sin
ｘ／ｘ）によつて変調され、略矩形の輪郭を持つ
作像スライス内の核スピンを選択的に励振する。
90゜RFパルス及びその他の選択性RFパルスはガ
ウス形包絡線関数で変調してもよく、この場合ス
ライス１０５の輪郭は近似的にガウス形になる。
１つのビユーを示すものであるが、このビユー
で、横軸に示した期間１に、正のGz勾配パルス
が印加されることが判る。Cz勾配の方向はデカル
ト座標系のｚ軸の正の方向に任意に選ぶが、これ
は磁界B0の方向と一致する。更に期間１に、Gz
勾配の存在の下に選択性90゜RFパルスを印加し
て、検査する物体の平面状容積（図に示してな
い）内にある核スピンを励振する。選ばれるスラ
イスの厚さがGz勾配の振幅と選択性90゜RFパルス
の周波数成分とによつて決定される。この発明の
好ましい実施例では、RFパルスがsinc関数（sin
ｘ／ｘ）によつて変調され、略矩形の輪郭を持つ
作像スライス内の核スピンを選択的に励振する。
90゜RFパルス及びその他の選択性RFパルスはガ
ウス形包絡線関数で変調してもよく、この場合ス
ライス１０５の輪郭は近似的にガウス形になる。
Gz勾配をターンオフした時、励振された核ス
ピンは同じ周波数で歳差運動をするが、勾配の位
相外し効果の為に、互いに位相がずれているφz。
励振された核スピンの位相コヒーレンスが、期間
２にこの場合は負のGz勾配パルスを印加するこ
とによつて再び設定される。典型的には、スピン
の位相戻しに必要な、期間２にわたるGz勾配の
波形の時間積分は、期間１に於けるGz勾配の波
形の時間積分の半分に負の符号をつけたものに大
体等しくなる様に選ぶ。ｚ軸方向のスピンの位相
戻しがグラフ２０によつて示されており、これか
ら、スピンの位相φzが期間２の終りにゼロに回
復することが判る。更に期間２の間、正のGx勾
配パルスを印加するのと同時に、位相符号化Gy
勾配を印加する。Gy勾配は、この順序の期間１
乃至７で構成された例としての１つのビユーの
間、１個のピーク振幅を持つている。然し、相次
ぐ各々のビユーでは、第２図の破線の曲線２１で
示す様に、異なる振幅の位相符号化勾配が選択さ
れる。Gy勾配が空間情報をｙ軸方向に位相符号
化する。この後のビユーではGy勾配の各々の相
異なる振幅が異なる程度の位相符号化を持ち込
む。プログラム可能なGy勾配の振幅の数は、再
生像がｙ軸方向に持つ画素の数（典型的には128，
256，512）に等しく選ぶ。この発明をGy勾配の
プログラム可能な振幅について説明するが、重要
なのは、パルスの期間にわたる勾配の波形の時間
積分であるから、位相符号化勾配の振幅並びに／
又は持続時間を変えることによつても位相符号化
を達成することが出来る。
ピンは同じ周波数で歳差運動をするが、勾配の位
相外し効果の為に、互いに位相がずれているφz。
励振された核スピンの位相コヒーレンスが、期間
２にこの場合は負のGz勾配パルスを印加するこ
とによつて再び設定される。典型的には、スピン
の位相戻しに必要な、期間２にわたるGz勾配の
波形の時間積分は、期間１に於けるGz勾配の波
形の時間積分の半分に負の符号をつけたものに大
体等しくなる様に選ぶ。ｚ軸方向のスピンの位相
戻しがグラフ２０によつて示されており、これか
ら、スピンの位相φzが期間２の終りにゼロに回
復することが判る。更に期間２の間、正のGx勾
配パルスを印加するのと同時に、位相符号化Gy
勾配を印加する。Gy勾配は、この順序の期間１
乃至７で構成された例としての１つのビユーの
間、１個のピーク振幅を持つている。然し、相次
ぐ各々のビユーでは、第２図の破線の曲線２１で
示す様に、異なる振幅の位相符号化勾配が選択さ
れる。Gy勾配が空間情報をｙ軸方向に位相符号
化する。この後のビユーではGy勾配の各々の相
異なる振幅が異なる程度の位相符号化を持ち込
む。プログラム可能なGy勾配の振幅の数は、再
生像がｙ軸方向に持つ画素の数（典型的には128，
256，512）に等しく選ぶ。この発明をGy勾配の
プログラム可能な振幅について説明するが、重要
なのは、パルスの期間にわたる勾配の波形の時間
積分であるから、位相符号化勾配の振幅並びに／
又は持続時間を変えることによつても位相符号化
を達成することが出来る。
期間２のGx勾配パルスの効果は、核スピンを
予定量だけ位相外しして、（普通のスピン捩れ形
順序の場合に）非選択性180゜RFパルスを期間３
に印加した時、スピンエコー信号が90゜RFパルス
から時間TE後の期間４に発生される様にするこ
とである。スピンエコー信号の発生時刻は、期間
２に印加されたGx勾配の強さ並びに期間４に印
加される直線的なGx勾配の強さによつて決定さ
れる。一般的に、あらゆる勾配の方向（この例で
はｘ，ｙ及びｚ）のスピンの位相（不動及び動く
両方のスピンに対し）がゼロに戻る時点にスピン
エコーが発生する。これを、ｘ軸方向の位相φx
について、第２図の点２２に示してある。ｘ軸方
向では、全てのスピンが不動であると仮定すれ
ば、期間２のパルスの時間積分（斜線を施した面
積２４）が期間４の勾配の時間積分（斜線を施し
た面積２６）と等しくなる様に、Gx勾配パルス
を選択することにより、こういうことが達成され
る。Gx勾配により、空間情報がｘ軸方向に符号
化され、前に第１図について説明した様に、ｘ軸
に対するその位置に特有な周波数で核スピンが共
鳴を起す。期間４に、典型的には再生像がｘ軸方
向に持つ画素の数（典型的には128，256，512）
に等しい回数だけ、スピンエコー信号を標本化す
る。周知の方法で、期間５に180゜RFパルスを印
加することにより、同じビユーの中で２番目のス
ピンエコー信号を期間６に発生することが出来
る。（２次元フーリエ変換方式の場合）２次元フ
ーリエ変換を用いて、標本化された信号から周知
の方法で画像の画素値が得られる。
予定量だけ位相外しして、（普通のスピン捩れ形
順序の場合に）非選択性180゜RFパルスを期間３
に印加した時、スピンエコー信号が90゜RFパルス
から時間TE後の期間４に発生される様にするこ
とである。スピンエコー信号の発生時刻は、期間
２に印加されたGx勾配の強さ並びに期間４に印
加される直線的なGx勾配の強さによつて決定さ
れる。一般的に、あらゆる勾配の方向（この例で
はｘ，ｙ及びｚ）のスピンの位相（不動及び動く
両方のスピンに対し）がゼロに戻る時点にスピン
エコーが発生する。これを、ｘ軸方向の位相φx
について、第２図の点２２に示してある。ｘ軸方
向では、全てのスピンが不動であると仮定すれ
ば、期間２のパルスの時間積分（斜線を施した面
積２４）が期間４の勾配の時間積分（斜線を施し
た面積２６）と等しくなる様に、Gx勾配パルス
を選択することにより、こういうことが達成され
る。Gx勾配により、空間情報がｘ軸方向に符号
化され、前に第１図について説明した様に、ｘ軸
に対するその位置に特有な周波数で核スピンが共
鳴を起す。期間４に、典型的には再生像がｘ軸方
向に持つ画素の数（典型的には128，256，512）
に等しい回数だけ、スピンエコー信号を標本化す
る。周知の方法で、期間５に180゜RFパルスを印
加することにより、同じビユーの中で２番目のス
ピンエコー信号を期間６に発生することが出来
る。（２次元フーリエ変換方式の場合）２次元フ
ーリエ変換を用いて、標本化された信号から周知
の方法で画像の画素値が得られる。
勾配内で移動するスピンに位相の累積が起るメ
カニズムを次に詳しく説明する。
カニズムを次に詳しく説明する。
このメカニズムを検討する為に、時刻ｔに座標
r〓（α＝ｘ，ｙ，ｚ）にあるスピンを考える。こ
の時、勾配G〓内でスピンが持つ位相φは時間の
級数展開で表わすことが出来る。
r〓（α＝ｘ，ｙ，ｚ）にあるスピンを考える。こ
の時、勾配G〓内でスピンが持つ位相φは時間の
級数展開で表わすことが出来る。
φ(t)＝∞
〓p=0
γ／ｐ！∫t 0G〓(t′)t′pdpr〓／dtpdt′(1)
ここでr〓のゼロでない時間微分が、スピンの運動
を特徴づけ、γはスピンの磁気回転比である。
を特徴づけ、γはスピンの磁気回転比である。
式(1)は次の様に書き直すことが出来る。
φ（ｔ）＝∞
〓p=0
φp（ｔ） (2)
ここでφp（ｔ）＝γ／ｐ！∫t 0G〓（t′）t′pdpr〓／
dtpdt′(3) ここでパルス順序の間、種々の位置微分（例えば
速度、加速度、加速度の変化率等）が一定である
と仮定する。この時 φp＝γ／ｐ！dpr〓(0)／dtp∫t 0G〓（t′）t′pdt′
(4) 式(4)中の積分項の部分が勾配分布G〓（ｔ）のｐ
次モーメントであることが理解されよう。
dtpdt′(3) ここでパルス順序の間、種々の位置微分（例えば
速度、加速度、加速度の変化率等）が一定である
と仮定する。この時 φp＝γ／ｐ！dpr〓(0)／dtp∫t 0G〓（t′）t′pdt′
(4) 式(4)中の積分項の部分が勾配分布G〓（ｔ）のｐ
次モーメントであることが理解されよう。
一般的に前にも触れたが、スピンエコーが形成
される為には、全てのスピンが同じ位相を持つて
いなければならない。第２図に示す様に便宜的に
この位相をゼロに選ぶことが出来る。この時、式
(2)から、全てのｐに対してφp＝０であることが
必要である。例えば不動のスピンに対しては、こ
れは次のことを意味する。
される為には、全てのスピンが同じ位相を持つて
いなければならない。第２図に示す様に便宜的に
この位相をゼロに選ぶことが出来る。この時、式
(2)から、全てのｐに対してφp＝０であることが
必要である。例えば不動のスピンに対しては、こ
れは次のことを意味する。
γr〓∫TE 0G〓（t′）dt′＝０、α〓ｙ
γry∫t7 t6
Gy（t′）dt′α＝ｙの場合
ここでTEはスピンエコーの発生時刻であり
（すなわち90゜RFパルスの照射時からスピンエコ
ーの発生までの時間を表わす）、t6及びt7は第３
図の期間３の位相符号化勾配の期間を定めるもの
である。スピンエコーを発生する様に、不動のス
ピンの位相戻しが行なわれる為には、普通のパル
ス順序で、α≠ｙの時、式(5)の第１の条件を充た
さなければならないことが理解されよう。これを
０次モーメントの補正と呼ぶ。不動のスピン以外
のスピンに対するこの発明の補正を高次補正と呼
ぶ。
（すなわち90゜RFパルスの照射時からスピンエコ
ーの発生までの時間を表わす）、t6及びt7は第３
図の期間３の位相符号化勾配の期間を定めるもの
である。スピンエコーを発生する様に、不動のス
ピンの位相戻しが行なわれる為には、普通のパル
ス順序で、α≠ｙの時、式(5)の第１の条件を充た
さなければならないことが理解されよう。これを
０次モーメントの補正と呼ぶ。不動のスピン以外
のスピンに対するこの発明の補正を高次補正と呼
ぶ。
或るスピンが速度（dr〓／dt）＝v〓を持つていれ
ば、動くスピンの位相戻しを行なう為には、勾配
が次の式を充たさなければならない。
ば、動くスピンの位相戻しを行なう為には、勾配
が次の式を充たさなければならない。
γv〓∫TE 0G〓（t′）t′dt′＝０ (6)
これが１番目の高次モーメントの補正である。
例えばスピンの加速度（２次）及び加速度の変化
率（３次）の様な他のモーメントに対して、追加
の補正を行なうことが出来る。例えばスピンの加
速度の影響を補償する条件は、 γa〓∫TE 0G〓（t′）t′2dt′＝０ (7) こゝでａはスピンの加速度である。勿論、一般
的には、勾配の全ての高次（０次より高い）モー
メントをゼロにすれば、ビユー毎の夫々の運動成
分の変動が、データに位相ジツタを誘起すること
はなく、従つてこのメカニズムによつて縞の人為
効果が発生しない。実際には、流れの人為効果を
減少する為には、全ての運動成分の補償は不要で
あることが判つた。例えば、速度成分に対する補
償だけでも流れの人為効果が目立つて減少するこ
とが判つた。
例えばスピンの加速度（２次）及び加速度の変化
率（３次）の様な他のモーメントに対して、追加
の補正を行なうことが出来る。例えばスピンの加
速度の影響を補償する条件は、 γa〓∫TE 0G〓（t′）t′2dt′＝０ (7) こゝでａはスピンの加速度である。勿論、一般
的には、勾配の全ての高次（０次より高い）モー
メントをゼロにすれば、ビユー毎の夫々の運動成
分の変動が、データに位相ジツタを誘起すること
はなく、従つてこのメカニズムによつて縞の人為
効果が発生しない。実際には、流れの人為効果を
減少する為には、全ての運動成分の補償は不要で
あることが判つた。例えば、速度成分に対する補
償だけでも流れの人為効果が目立つて減少するこ
とが判つた。
勾配のモーメントをゼロにするこの発明の方法
を、選択的な180゜RFパルス順序（TE1＝
25msec，TE2＝80msec）を用いたエコー２つの
多重スライス順序で実施した。その１つのビユー
が第３図に示してあり、０次及び１次モーメント
をゼロにする。この発明の方法は２つより多くの
スピンエコー信号を持つ順序に適用し得るが後掲
係属中の米国特許出願に記載されているようなゴ
ースト抑圧方法を保存する必要がある為に、実施
するのが一層困難である。一般的に、このパルス
順序は第２図について説明したパルス順序と同様
である。即ち、スライスの選択が、矩形のGz勾
配パルスの存在の下に選択性90゜RFパルスを印加
することによつて期間１に行なわれる。正弦の半
分で構成された位相符号化Gy勾配が期間２に印
加され、選択性180゜RFパルスが期間４及び１０
に矩形Gx勾配パルスの存在の下に印加され、こ
の結果期間６及び１３に生ずるスピンエコー信号
を直線的なGx勾配読出パルスの存在の下に観測
する。
を、選択的な180゜RFパルス順序（TE1＝
25msec，TE2＝80msec）を用いたエコー２つの
多重スライス順序で実施した。その１つのビユー
が第３図に示してあり、０次及び１次モーメント
をゼロにする。この発明の方法は２つより多くの
スピンエコー信号を持つ順序に適用し得るが後掲
係属中の米国特許出願に記載されているようなゴ
ースト抑圧方法を保存する必要がある為に、実施
するのが一層困難である。一般的に、このパルス
順序は第２図について説明したパルス順序と同様
である。即ち、スライスの選択が、矩形のGz勾
配パルスの存在の下に選択性90゜RFパルスを印加
することによつて期間１に行なわれる。正弦の半
分で構成された位相符号化Gy勾配が期間２に印
加され、選択性180゜RFパルスが期間４及び１０
に矩形Gx勾配パルスの存在の下に印加され、こ
の結果期間６及び１３に生ずるスピンエコー信号
を直線的なGx勾配読出パルスの存在の下に観測
する。
式(5)及び(6)によつて定められた条件を充たす様
な形で、勾配内で移動するスピンの位相の累積を
補償する為の勾配パルスも、第３図のパルス順序
内で、ｘ，ｙ及びｚ軸に沿つて印加される。スピ
ンエコー１に対して補償を行なう為の勾配パルス
が、この図の数字３０で示した線の左側で印加さ
れ、スピンエコー２に対する補償パルスがその右
側で印加される。この為、スピンエコー１に対し
ては、補償用のGx，Gy及びGz勾配パルスが期間
３及び５に印加される。２番目のスピンエコーに
対する補償用のGx及びGz勾配パルスが期間７，
９及び１１に印加される。夫々の場合、Gy及び
Gz補償パルスは正弦の半分の形であり、Gx補償
パルスが矩形の形であることが認められよう。例
として云うと、スピンエコー１に対するGz方向
の補償（期間６）では、式(5)及び(6)に相当する特
定の式は次の通りである。
な形で、勾配内で移動するスピンの位相の累積を
補償する為の勾配パルスも、第３図のパルス順序
内で、ｘ，ｙ及びｚ軸に沿つて印加される。スピ
ンエコー１に対して補償を行なう為の勾配パルス
が、この図の数字３０で示した線の左側で印加さ
れ、スピンエコー２に対する補償パルスがその右
側で印加される。この為、スピンエコー１に対し
ては、補償用のGx，Gy及びGz勾配パルスが期間
３及び５に印加される。２番目のスピンエコーに
対する補償用のGx及びGz勾配パルスが期間７，
９及び１１に印加される。夫々の場合、Gy及び
Gz補償パルスは正弦の半分の形であり、Gx補償
パルスが矩形の形であることが認められよう。例
として云うと、スピンエコー１に対するGz方向
の補償（期間６）では、式(5)及び(6)に相当する特
定の式は次の通りである。
−Gz1t1−Gz2（t2−t1）−２／πGz3（t3−t2）
＋２／πGz5（t5−t4）−０ (8)
−G1t2／1／２−G22／２（t2 2−t2 1）−２／πGz3／２
（t2 3−t2 2） −Gz4t3t4＋２／πGz5／２（t2 5−t2 4）＝０ (9) 同様に、スピンエコー１に対するGy勾配方向
の補償では、式は次の様になる。
（t2 3−t2 2） −Gz4t3t4＋２／πGz5／２（t2 5−t2 4）＝０ (9) 同様に、スピンエコー１に対するGy勾配方向
の補償では、式は次の様になる。
−Gy2（t9−t8）＋Gy3（t11−t10）＝０ (10)
−Gy1／２（t2 7−t2 6）−Gy2／２（t2 9−t2 8）＋Gy3／
２（t2 11−t2 10） ＝０ (11) Gx勾配方法の補償では、 −Gx1（t13−t12）＋Gx2（t15−t14）＋Gx3（TE−t16）
＝
０ (12) −Gx1／２（t2 13−t2 12）＋Gx2／２（t2 15−t2 14）＋G
×3／２（TE2 −t2 16）＝０ （13） 式(7)及至(12)で量G〓i（α＝ｘ，ｙ，ｚ，ｉ＝１，
２，３…）は勾配パルスの振幅を表わし、tj（ｊ
＝１，２，３…）は第３図に示す様に、勾配パル
スの初め及び終りの時を表わし、ｔ＝０から始ま
る。特定の式(8)、(10)及び(12)が式(5)に対応し、式
(9)、(11)及び（13）が(6)に対応する。
２（t2 11−t2 10） ＝０ (11) Gx勾配方法の補償では、 −Gx1（t13−t12）＋Gx2（t15−t14）＋Gx3（TE−t16）
＝
０ (12) −Gx1／２（t2 13−t2 12）＋Gx2／２（t2 15−t2 14）＋G
×3／２（TE2 −t2 16）＝０ （13） 式(7)及至(12)で量G〓i（α＝ｘ，ｙ，ｚ，ｉ＝１，
２，３…）は勾配パルスの振幅を表わし、tj（ｊ
＝１，２，３…）は第３図に示す様に、勾配パル
スの初め及び終りの時を表わし、ｔ＝０から始ま
る。特定の式(8)、(10)及び(12)が式(5)に対応し、式
(9)、(11)及び（13）が(6)に対応する。
各々の軸に対してこれらの式を別々に解く計算
機のプログラムを工夫することが出来ることが理
解されよう。これらの式を解く時、一般的には、
変更しないパルス順序で存在する勾配パルス（例
えばスライス選択用のGz、位相符号化用のGy、
読出用のGx）の他に、各々の軸に沿つて、２つ
の勾配パルスを選択し、これらの２つの勾配パル
スは作像の条件によつて拘束されず、そのパラメ
ータの値（振幅及び持続時間）が、これらの式の
解によつて決定される様にすることが望ましい。
実際には、スピンに対する勾配の影響が振幅と持
続時間の積（即ち、勾配パルスの波形の下にある
面積）に関係するから、振幅及び持続時間の特定
の値は重要ではない。或るパルス順序では、現存
の勾配パルスの振幅及び／又は持続時間の値を調
節することにより、式(5)及び(6)の条件を充たすこ
とが出来るので、こういう追加のパルスを選択す
ることは必要ではないことがある。
機のプログラムを工夫することが出来ることが理
解されよう。これらの式を解く時、一般的には、
変更しないパルス順序で存在する勾配パルス（例
えばスライス選択用のGz、位相符号化用のGy、
読出用のGx）の他に、各々の軸に沿つて、２つ
の勾配パルスを選択し、これらの２つの勾配パル
スは作像の条件によつて拘束されず、そのパラメ
ータの値（振幅及び持続時間）が、これらの式の
解によつて決定される様にすることが望ましい。
実際には、スピンに対する勾配の影響が振幅と持
続時間の積（即ち、勾配パルスの波形の下にある
面積）に関係するから、振幅及び持続時間の特定
の値は重要ではない。或るパルス順序では、現存
の勾配パルスの振幅及び／又は持続時間の値を調
節することにより、式(5)及び(6)の条件を充たすこ
とが出来るので、こういう追加のパルスを選択す
ることは必要ではないことがある。
４種類の異なる直径の管を直径に配置して、
夫々両方の流れの方向を持つ様にした流れのフア
ントムを用いて、この発明の方法による実験デー
タを収集した。直径が一番大きい管でも非層流が
発生する位に高い速度で、硫酸銅（CuSO4）でド
ープした水を圧送した。一番小さい管に於ける流
量は100cm／秒を優に越えたが、正確な流量は測
定しなかつた。
夫々両方の流れの方向を持つ様にした流れのフア
ントムを用いて、この発明の方法による実験デー
タを収集した。直径が一番大きい管でも非層流が
発生する位に高い速度で、硫酸銅（CuSO4）でド
ープした水を圧送した。一番小さい管に於ける流
量は100cm／秒を優に越えたが、正確な流量は測
定しなかつた。
この発明の補償方法を用いる第４図の順序を用
いた初期の実験は、180゜RFパルス（期間４）が
選択性でない時は成功したが、選択勾配をターン
オンした時、補償はそれ程有効でなかつた。
180゜RFパルスによる追加の位相雑音が問題であ
ることが判つた。その解決策は、モーメントの長
い勾配ロープ（第４図の期間２のGx及びGz）を
取去り、第３図の様に、流れ補償パルスを出来る
だけ180゜RFパルスに近づけることであつた。そ
の様にした時、選択性パルスを用いても、補償は
満足し得るものであることが判つた。
いた初期の実験は、180゜RFパルス（期間４）が
選択性でない時は成功したが、選択勾配をターン
オンした時、補償はそれ程有効でなかつた。
180゜RFパルスによる追加の位相雑音が問題であ
ることが判つた。その解決策は、モーメントの長
い勾配ロープ（第４図の期間２のGx及びGz）を
取去り、第３図の様に、流れ補償パルスを出来る
だけ180゜RFパルスに近づけることであつた。そ
の様にした時、選択性パルスを用いても、補償は
満足し得るものであることが判つた。
どの方向の補償が必要であるかを決定する実験
も行なつた。この実験は、特定の方向（ｘ，ｙ又
はｚ）の補償をターンオブした時、パルス順序が
その軸で普通のものに戻る様に構成した。
も行なつた。この実験は、特定の方向（ｘ，ｙ又
はｚ）の補償をターンオブした時、パルス順序が
その軸で普通のものに戻る様に構成した。
流れが公称ｚ方向（磁石の軸線）であつたが、
人為効果を実質的に少なくする為には、ｘ及びｚ
の両方の補償を利用しなければならないことが判
つた。この意味するところは、ｚ方向と同じく横
方向にも実質的な変動があるということである。
ｙ（位相符号化）方向の補償はそれ程重要である
様に思われないが、これはｘ（読出）方向よりも、
ｙ方向の動きが少ない又は流速の変動が少ない為
ではない。この様に影響が低下するのは、ｙ軸に
は90゜RFパルスに接近した丁度１つの勾配ローブ
があり、従つて、ｙ方向の速度変動による位相誤
差が累積する発展時間が殆んどない為である。
人為効果を実質的に少なくする為には、ｘ及びｚ
の両方の補償を利用しなければならないことが判
つた。この意味するところは、ｚ方向と同じく横
方向にも実質的な変動があるということである。
ｙ（位相符号化）方向の補償はそれ程重要である
様に思われないが、これはｘ（読出）方向よりも、
ｙ方向の動きが少ない又は流速の変動が少ない為
ではない。この様に影響が低下するのは、ｙ軸に
は90゜RFパルスに接近した丁度１つの勾配ローブ
があり、従つて、ｙ方向の速度変動による位相誤
差が累積する発展時間が殆んどない為である。
軸平面（ｘ−ｙ平面）及びサジツタル平面（ｘ
−ｚ平面）の両方で志願者に対する検査も行なつ
た。走査は心臓ゲート形ではなかつたが、非逐次
形位相符号化順序を用いて、呼吸による人為効果
を減少した。こういう順序が出願人の係属中の米
国特許出願通し番号第766733号及び同第766842号
に記載されている。補償を用いた時、大動脈及び
大静脈の流れによる人為効果が著しく減少し、血
管は暗く見える代りに、明るく見えた。サジツタ
ル像では、補償を用いた時、血液（普通の像では
見えないのが普通である）が可視的になつた。
80msecに於ける２番目のエコーの人為効果の抑
圧は、１番目のエコー程よくないことが判つた。
これは２つの180゜RFパルスの間の一層長い時間
の為に重要な飛行時間誤差の為である。
−ｚ平面）の両方で志願者に対する検査も行なつ
た。走査は心臓ゲート形ではなかつたが、非逐次
形位相符号化順序を用いて、呼吸による人為効果
を減少した。こういう順序が出願人の係属中の米
国特許出願通し番号第766733号及び同第766842号
に記載されている。補償を用いた時、大動脈及び
大静脈の流れによる人為効果が著しく減少し、血
管は暗く見える代りに、明るく見えた。サジツタ
ル像では、補償を用いた時、血液（普通の像では
見えないのが普通である）が可視的になつた。
80msecに於ける２番目のエコーの人為効果の抑
圧は、１番目のエコー程よくないことが判つた。
これは２つの180゜RFパルスの間の一層長い時間
の為に重要な飛行時間誤差の為である。
この発明を特定の実施例及び例について説明し
たが、以上の説明から、当業者にはこの他の変更
が考えられよう。従つて、特許請求の範囲内で、
この発明はこゝに具体的に説明しなかつた形で実
施することが出来ることを承知されたい。
たが、以上の説明から、当業者にはこの他の変更
が考えられよう。従つて、特許請求の範囲内で、
この発明はこゝに具体的に説明しなかつた形で実
施することが出来ることを承知されたい。
第１図は動くスピンの位相の累積が起る様子を
理解するのに役立つ横平面（ｘ−ｙ平面）内の勾
配を示すグラフ、第２図はｘ及びｚ軸方向のスピ
ンの位相のグラフを含めて、エコー２個のスピン
捩れ形パルス順序の１つのビユーを示すグラフ、
第３図は勾配のモーメントをゼロにすることを用
いたこの発明によるスピン流れ補償を含むエコー
２個のスピン捩れ形パルス順序の１つのビユーを
示すグラフ、第４図は180゜RFパルスから離れて
発生する大きなGx及びGz磁界勾配パルスの為に、
流れの補償効果が幾分低下するスピン捩れ形パル
ス順序の１つのビユーを示すグラフである。
理解するのに役立つ横平面（ｘ−ｙ平面）内の勾
配を示すグラフ、第２図はｘ及びｚ軸方向のスピ
ンの位相のグラフを含めて、エコー２個のスピン
捩れ形パルス順序の１つのビユーを示すグラフ、
第３図は勾配のモーメントをゼロにすることを用
いたこの発明によるスピン流れ補償を含むエコー
２個のスピン捩れ形パルス順序の１つのビユーを
示すグラフ、第４図は180゜RFパルスから離れて
発生する大きなGx及びGz磁界勾配パルスの為に、
流れの補償効果が幾分低下するスピン捩れ形パル
ス順序の１つのビユーを示すグラフである。
Claims (1)
- 【特許請求の範囲】 １ 時刻TEにスピンエコー信号を発生させるよ
うに不動の核スピン並びに流れる核スピンを持つ
物体の領域をラーモア周波数のRF励振パルスで
照射し、前記領域の少なくとも１つの次元の軸線
に沿つて調節自在の振幅及び持続時間のパラメー
タを持つパルス形磁界勾配を印加することを何れ
も含む様な複数個のビユーを求める場合におい
て、前記物体の領域から収集される磁気共鳴信号
に、前記流れるスピンの流速がビユー毎に変動す
ることによつて起る不規則な位相の累積を補償す
る方法に於て、 αをｘ、ｙ又はｚとしv〓をαの方向の核スピン
の流速、TEをスピンエコー信号の発生時刻、G〓
をα方向の勾配、γを磁気回転比として、前記流
れる核スピンについて次の式 γv〓∫TE 0G〓（t′）t′dt′＝０ を充たす様に、前記振幅及び持続時間のパラメー
タの内の少なくとも１つを変更して、前記磁気勾
配の少なくとも１つの高次モーメントをゼロに等
しくする工程を含む方法。 ２ 特許請求の範囲１に記載した方法に於て、更
に、前記領域の第１の次元の軸線に沿つて少なく
とも１つの磁界勾配を印加して、前記第１の磁界
の方向に、空間情報を前記スピンエコー信号に符
号化する工程を含む方法。 ３ 特許請求の範囲２に記載した方法に於て、情
報を符号化する為の磁界勾配が、前記ビユーを求
める過程で、前記領域内の複数個の方向の内の１
つを持つ様に選択される方法。 ４ 特許請求の範囲３に記載した方法に於て、情
報を符号化する為の勾配が、前記領域内で互いに
垂直である少なくとも２つの磁界勾配のベクトル
加算の合成磁界勾配で構成される方法。 ５ 特許請求の範囲３に記載した方法に於て、前
記情報を符号化する為の勾配が、その内の少なく
とも１つが他の磁界勾配と同一平面でない様な複
数個の磁界勾配のベクトル加算の合成磁界勾配で
構成される方法。 ６ 特許請求の範囲２に記載した方法に於て、更
に、照射してスピンエコー信号を発生する工程の
前に、前記複数個のビユーを求める過程で前記ス
ピンエコー信号に空間情報を位相符号化する為
に、複数個の振幅及び持続時間の積の内の１つを
持つ少なくとも１つの位相符号化磁界勾配を印加
する工程を含む方法。 ７ 特許請求の範囲６に記載した方法に於て、前
記スピンエコー信号が略直線的な磁界勾配の存在
の下に標本化され、該勾配が前記領域内で前記位
相符号化磁界勾配の方向に対して垂直である方
法。 ８ 特許請求の範囲７に記載した方法に於て、前
記１つの位相符号化磁界勾配に対して直交する様
に選ばれた別の位相符号化磁界勾配を印加する工
程を含み、該別の位相符号化磁界勾配は複数個の
振幅と持続時間の積を持つており、 前記１つの勾配の振幅及び持続時間の積を進め
る前に、前記別の位相符号化勾配の全ての振幅と
持続時間の積を順次進める間、１つの位相符号化
勾配の振幅と持続時間の積を一定に保つ工程を含
む方法。 ９ 特許請求の範囲１に記載した方法に於て、 (a) 前記領域の第１の次元の軸線に沿つて印加し
た第１の勾配の存在の下に、前記領域内の核ス
ピンを選択的に励振して、スピンエコー信号を
発生させるようにし、 (b) 前記領域の第２の次元の軸線に沿つて位相符
号化磁界勾配を印加して、該第２の軸線の方向
の空間情報を前記スピンエコー信号に付号化
し、 (c) 前記軸線に直交する前記領域の第３の次元の
軸線に沿つて読出磁界勾配の存在の下に前記ス
ピンエコー信号を観測して、前記第３の軸線に
沿つた空間情報を符号化する工程を含む方法。 １０ 時刻TEにスピンエコー信号を発生させる
ように不動の核スピン並びに流れる核スピンを持
つ物体の領域をラーモア周波数のRF励振パルス
で照射し、前記領域の少なくとも１つの次元の軸
線に沿つて調節自在の振幅及び持続時間のパラメ
ータを持つパルス形磁界勾配を印加することを何
れも含む様な複数個のビユーを求める場合におい
て、前記物体の領域から収集される磁気共鳴信号
に、前記流れるスピンの流速がビユー毎に変動す
ることによつて起る不規則な位相の累積を補償す
る装置に於て、 αをｘ，ｙ又はｚとし、v〓をαの方向の核スピ
ンの流速、TEをスピンエコー信号の発生時刻、
G〓をα方向の勾配、γを磁気回転比として、前
記流れる核スピンについて次の式 γv〓∫TE 0G〓（t′）t′dt′＝０ を充たす様に、前記振幅及び持続時間のパラメー
タの内の少なくとも１つを変更して、前記磁界勾
配の少なくとも１つの高次モーメントをゼロに等
しくする手段を含む装置。 １１ 特許請求の範囲１０に記載した装置に於
て、更に、前記領域の第１の次元の軸線に沿つて
少なくとも１つの磁界勾配を印加して、前記第１
の磁界の方向に、空間情報を前記スピンエコー信
号に符号化する手段を含む装置。 １２ 特許請求の範囲１１に記載した装置に於
て、情報を符号化する為の磁界勾配が、前記ビユ
ーを求める過程で、前記領域内の複数個の方向の
内の１つを持つ様に選択される装置。 １３ 特許請求の範囲１２に記載した装置に於
て、情報を符号化する為の勾配が、前記領域内で
互いに垂直である少なくとも２つの磁界勾配のベ
クトル加算の合成磁界勾配で構成される装置。 １４ 特許請求の範囲１２に記載した装置に於
て、前記情報を符号化する為の勾配が、その内の
少なくとも１つが他の磁界勾配と同一平面でない
様な複数個の磁界勾配のベクトル加算の合成磁界
勾配で構成される装置。 １５ 特許請求の範囲１１に記載した装置に於
て、更に、照射してスピンエコー信号を発生する
前に、前記複数個のビユーを求める過程で前記ス
ピンエコー信号に空間情報を位相符号化する為
に、複数個の振幅及び持続時間の積の内の１つを
持つ少なくとも１つの位相符号化磁界勾配を印加
する手段を含む装置。 １６ 特許請求の範囲１５に記載した装置に於
て、前記スピンエコー信号が略直線的な磁界勾配
の存在の下に標本化され、該勾配が前記領域内で
前記位相符号化磁界勾配の方向に対して垂直であ
る装置。 １７ 特許請求の範囲１６に記載した装置に於
て、前記１つの位相符号化磁界勾配に対して直交
する様に選ばれた別の位相符号化磁界勾配を印加
する手段を含み、該別の位相符号化磁界勾配は複
数個の振幅と持続時間の積を持つており、 前記１つの勾配の振幅及び持続時間の積を進め
る前に、前記別の位相符号化勾配の全ての振幅と
持続時間の積を順次進める間、１つの位相符号化
勾配の振幅と持続時間の積を一定に保つ手段を含
む装置。 １８ 特許請求の範囲１０に記載した装置に於
て、 (a) 前記領域の第１の次元の軸線に沿つて印加し
た第１の勾配の存在の下に、前記領域内の核ス
ピンを選択的に励振して、スピンエコー信号を
発生させるようにする手段、 (b) 前記領域の第２の次元の軸線に沿つて位相符
号化磁界勾配を印加して、該第２の軸線の方向
の空間情報を前記スピンエコー信号に符号化す
る手段、 (c) 前記軸線に直交する前記領域の第３の次元の
軸線に沿つて読出磁界勾配の存在の下に前記ス
ピンエコー信号を観測して、前記第３の軸線に
沿つた空間情報を符号化する手段を含む装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US06/798,544 US4731583A (en) | 1985-11-15 | 1985-11-15 | Method for reduction of MR image artifacts due to flowing nuclei by gradient moment nulling |
| US798544 | 1985-11-15 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS62161354A JPS62161354A (ja) | 1987-07-17 |
| JPH0351176B2 true JPH0351176B2 (ja) | 1991-08-06 |
Family
ID=25173671
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP61270030A Granted JPS62161354A (ja) | 1985-11-15 | 1986-11-14 | 不規則な位相の累積を補償する方法と装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4731583A (ja) |
| EP (1) | EP0222325A3 (ja) |
| JP (1) | JPS62161354A (ja) |
| IL (1) | IL80134A0 (ja) |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4728890A (en) * | 1985-08-16 | 1988-03-01 | Picker International, Inc. | Motion artifact suppression technique of magnetic resonance imaging |
| US4937526A (en) * | 1988-11-23 | 1990-06-26 | Mayo Foundation For Medical Education And Research | Adaptive method for reducing motion and flow artifacts in NMR images |
| JP2776844B2 (ja) * | 1988-11-30 | 1998-07-16 | 株式会社日立製作所 | 磁気共鳴動体イメージング装置 |
| JPH02237542A (ja) * | 1989-03-10 | 1990-09-20 | Yokogawa Medical Syst Ltd | Ｍｒｉシーケンスの面外ごみ除去方法 |
| EP0401908B1 (en) * | 1989-06-07 | 1996-05-15 | Koninklijke Philips Electronics N.V. | Magnetic resonance method and device for acquiring short echo time images |
| EP0401429A1 (en) * | 1989-06-07 | 1990-12-12 | Koninklijke Philips Electronics N.V. | Magnetic resonance method and device for acquiring short echotime images |
| US4959611A (en) * | 1989-06-29 | 1990-09-25 | Picker International, Inc. | Out-of-slice artifact reduction technique for magnetic resonance imagers |
| US4973906A (en) * | 1989-08-17 | 1990-11-27 | General Electric Company | Flow compensated NMR fast pulse sequence |
| US5093620A (en) * | 1990-08-09 | 1992-03-03 | General Electric | Encoding for nmr phase contrast flow measurement |
| US5280244A (en) * | 1992-03-19 | 1994-01-18 | General Electric Company | Gradient moment nulling in a fast spin echo NMR pulse sequence |
| US5729140A (en) * | 1996-10-04 | 1998-03-17 | Mayo Foundation For Medical Education And Research | Suppression of artifacts in NMR images by correlation of two NEX acquistions |
| US5926022A (en) * | 1997-04-09 | 1999-07-20 | Mayo Foundation For Medical Education And Research | Method for flow compensating a multi-shot echo-planar MRI pulse sequence |
| US6184682B1 (en) | 1998-04-15 | 2001-02-06 | Mayo Foundation For Medical Education And Research | Correction of MR images for motion artifacts using navigator echoes and autocorrection |
| US6281681B1 (en) | 1999-01-28 | 2001-08-28 | General Electric Company | Magnetic resonance imaging with interleaved Fibonacci spiral scanning |
| US6647134B1 (en) * | 2000-03-30 | 2003-11-11 | Mayo Foundation For Medical Education And Research | Autocorrection of MR projection images |
| US7558612B2 (en) * | 2002-04-16 | 2009-07-07 | The Board Of Trustees Of The Leland Stanford Junior University | Motion compensated spiral FISP MRI |
| US7457655B2 (en) * | 2004-10-08 | 2008-11-25 | Mayo Foundation For Medical Education And Research | Motion correction of magnetic resonance images using moments of spatial projections |
| US10098563B2 (en) * | 2006-11-22 | 2018-10-16 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus |
| JP2010525892A (ja) * | 2007-05-04 | 2010-07-29 | カリフォルニア インスティテュート オブ テクノロジー | 低磁場ｓｑｕｉｄ−ｍｒｉ装置、コンポーネント、及び方法 |
| WO2009023303A2 (en) * | 2007-05-04 | 2009-02-19 | California Institute Of Technology | Geometries for superconducting sensing coils for squid-based systems |
| US8981775B2 (en) * | 2009-09-17 | 2015-03-17 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and gradient magnetic field application method |
| WO2018160839A1 (en) * | 2017-03-01 | 2018-09-07 | The Regents Of The University Of California | System and method for improved diffusion-weighted imaging |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4431968A (en) * | 1982-04-05 | 1984-02-14 | General Electric Company | Method of three-dimensional NMR imaging using selective excitation |
| JPS6038644A (ja) * | 1983-08-12 | 1985-02-28 | Hitachi Ltd | 核磁気共鳴を用いた検査装置 |
| US4595879A (en) * | 1983-11-14 | 1986-06-17 | Technicare Corporation | Nuclear magnetic resonance flow imaging |
| US4616180A (en) * | 1983-11-14 | 1986-10-07 | Technicare Corporation | Nuclear magnetic resonance imaging with reduced sensitivity to motional effects |
| US4629988A (en) * | 1984-07-02 | 1986-12-16 | General Electric Company | Method of imaging by depth-resolved surface coil spectroscopy |
| DE3673107D1 (de) * | 1985-08-13 | 1990-09-06 | Gen Electric Cgr | Modulationsverfahren des geschwindigkeitseffekts von beweglichen teilchen eines koerpers in einer dichtemessung mittels der kernmagnetischen resonanz. |
| US4663591A (en) * | 1985-08-16 | 1987-05-05 | General Electric Company | Method for reducing image artifacts due to periodic signal variations in NMR imaging |
-
1985
- 1985-11-15 US US06/798,544 patent/US4731583A/en not_active Expired - Lifetime
-
1986
- 1986-09-24 IL IL80134A patent/IL80134A0/xx not_active IP Right Cessation
- 1986-11-06 EP EP86115389A patent/EP0222325A3/en not_active Withdrawn
- 1986-11-14 JP JP61270030A patent/JPS62161354A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| EP0222325A3 (en) | 1989-06-07 |
| IL80134A0 (en) | 1986-12-31 |
| EP0222325A2 (en) | 1987-05-20 |
| US4731583A (en) | 1988-03-15 |
| JPS62161354A (ja) | 1987-07-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH0351176B2 (ja) | ||
| US4532473A (en) | NMR method for measuring and imaging fluid flow | |
| CN1138156C (zh) | 用于扩散加权的成像的快速自旋回波脉冲序列 | |
| US7372269B2 (en) | Magnetic resonance imaging method and apparatus | |
| JP5325795B2 (ja) | 位相差合成画像を用いた高度に制約された逆投影再構成法を利用した造影ｍｒａ | |
| JP3878176B2 (ja) | インターリービングされた投影−再構成データを使用する三次元位相コントラスト磁気共鳴イメージング | |
| US7706857B2 (en) | Methods and apparatus for mapping internal and bulk motion of an object with phase labeling in magnetic resonance imaging | |
| JPH0614914B2 (ja) | 多次元再構成技術を使用する高速ｎｍｒ映像化方法及びその装置 | |
| US20070167729A1 (en) | Highly constrained magnetic resonance spectroscopy image reconstruction method | |
| EP1145028B1 (en) | Phase contrast mr flow imaging using angularly interleaved projection data | |
| JPH0644904B2 (ja) | 磁界強度の空間的な変動を決定する装置，非均質性マップからシム勾配を計算する装置，非均質性マップで表されている磁界を補正する装置，及び中心外れの２次成分を有する磁界を補正する装置 | |
| JP2769473B2 (ja) | Ｎｍｒにより身体内の体積素の運動を表示する装置 | |
| JPH0616766B2 (ja) | Ｎｍｒイメージング装置を使用した高速流れ測定方式 | |
| US5798642A (en) | Magnetic resonance imaging apparatus | |
| JPH0373300B2 (ja) | ||
| EP1693680B1 (en) | Continuous table motion MRI involving phase correction | |
| JP5291852B2 (ja) | パルス列、核磁気共鳴断層撮影装置及びイメージング法 | |
| JP3588690B2 (ja) | 磁気共鳴装置 | |
| JP2002000579A (ja) | Ｍｒ拡散画像撮影装置 | |
| EP0153703A2 (en) | NMR imaging apparatus | |
| JP3322695B2 (ja) | 磁気共鳴イメージング装置 | |
| Müller-Bierl | MRI Signal Formation: From Spins to Measured Signal | |
| JP3152690B2 (ja) | 磁気共鳴映像装置 | |
| JPH0470013B2 (ja) | ||
| JP3249114B2 (ja) | Ｍｒｉ装置及びｍｒｉ装置における傾斜磁場印加方法 |