JPH0424051B2 - - Google Patents

Info

Publication number
JPH0424051B2
JPH0424051B2 JP63027838A JP2783888A JPH0424051B2 JP H0424051 B2 JPH0424051 B2 JP H0424051B2 JP 63027838 A JP63027838 A JP 63027838A JP 2783888 A JP2783888 A JP 2783888A JP H0424051 B2 JPH0424051 B2 JP H0424051B2
Authority
JP
Japan
Prior art keywords
sample
response
pair
data
nmr
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
Application number
JP63027838A
Other languages
English (en)
Other versions
JPS63226343A (ja
Inventor
Rushian Deyumorin Chaaruzu
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Electric Co
Original Assignee
General Electric Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by General Electric Co filed Critical General Electric Co
Publication of JPS63226343A publication Critical patent/JPS63226343A/ja
Publication of JPH0424051B2 publication Critical patent/JPH0424051B2/ja
Granted legal-status Critical Current

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/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

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)
  • Measuring Volume Flow (AREA)

Description

【発明の詳細な説明】 発明の背景 この発明はサンプル中の流体の流れを作像する
核磁気共鳴（NMR）造影法、更に具体的に云え
ば、医学的に重要な流体の流れの解剖像を非侵入
形で求めるため等のために、サンプルの励振毎に
多数の応答が発生されるような新規のNMR造影
法に関する。
医学的な診断のために、種々の身体通路を通る
体液の流れを示すNMR造影法のデータ像を提供
することが知られている。このような像を発生す
る方法が、1986年３月３日に出願された係属中の
米国特許出願通し番号第835683号（米国特許第
4714081号）に記載されている。この方法は、真
の投影像を作り、選ばれた投影軸線に沿つた静脈
及び動脈構造の品質の高いNMR造影図が得ら
れ、サンプル中の選ばれた流れの方向に対する感
度を持つが、NMR造影図の情報を更に増やすこ
とが依然として非常に望ましい。例えば、何れも
直交する流れの成分に敏感な複数個のNMR造影
図を組合せることができれば、流れ全体のNMR
造影図を作ることができるので、多数の投影軸線
に沿つた投影像が見られることが望ましい。更
に、（１つの投影軸線に沿つてでも）NMR造影
図の信号対雑音比を改善して、動きの人為効果
（motion artifact）の抑圧を良くすると共に、異
なる範囲の流速に何れも敏感な一連のNMR造影
図を求めることが望ましい。しかし、心臓サイク
ル当たり１つのパルス順序しか使わない心臓ゲー
ト形NMR造影図では、１つの流れの方向だけで
スピン捩れ形作像順序に128個の位相符号化工程
を必要とするが、最低の作像時間は、心室収縮率
に関係するが、現在、約４分を必要とする。流れ
の方向が２つの場合、２倍の励振順序を必要とす
るが、患者（サンプル）からのデータを収集する
には、約16分を必要とすることがある。こういう
期間にわたつて患者を動かない状態に保つことは
極めて困難であり、余分のビユー（view）に対
して必要な更に長い期間に対しては、なおさらで
ある。従つて、実質的な余分の収集時間を必要と
せず、且つサンプルの動きによる整合外れの影響
を受けずに、改善された特徴を持つNMR造影図
を提供する方法が非常に望ましい。
発明の簡単な要約 この発明では、サンプル内の流体の流れに伴う
動くスピンだけを実質的に表わす核磁気共鳴造影
像を得る方法と装置が提供される。この発明で
は、サンプルを種静磁界の中に浸漬し、サンプル
の選ばれた部分の多数（Ｓ個）の領域の各々に対
する作像順序の逐次的な対の内の第１の順序及び
第２の順序の各々の初期部分で、選ばれた種目の
全ての原子核のスピンを章動させ、動く原子核の
スピンから得られるNMR応答エコー信号が、略
不動の原子核のスピンから得られるNMR応答エ
コー信号と異なるように選ばれた第１の方向に、
サンプルに印加される第１の磁界勾配中に極性が
交互に変わる一対の流れ符号化信号パルスを印加
し、各対の第１の順序中の一対の流れ符号化パル
スは、各対の第２の順序中の一対の流れ符号化パ
ルスとは極性が逆になつており、前記第１の方向
に対して略直交する第２の方向に前記サンプルに
印加された読取磁界勾配に応答して、第１及び第
２の順序の各々の複数個（Ｎ個）の応答データ収
集期間の夫々相異なる１つの期間中に、サンプル
の少なくとも前記部分から誘起されたNMR応答
エコー信号から一組のデータを夫々収集し、１≦
ｊ≦Ｎとして、前記第１の順序で収集されたＮ組
のNMR応答信号データの内の夫々ｊ番目の組の
データを、前記第２の順序収集された同じ番号
（ｊ番目）の組のデータから減算して、不動の原
子核によつて得られる応答データが実質的に除去
されたＮ組の差データの内の対応するｊ番目の組
の差データを夫々発生し、Ｎ組の差データの全て
に応答して、何れも前記第１及び第２の方向に対
して予め選ばれた関係を持つ平面内にある少なく
とも１つの造影投影像の各々を発生する工程を含
む。
現在好ましいと考えている実施例の多重エコー
NMR造影法では、流れ符号化工程は、１個の
NMR造影像を求めるのに必要なのと略同じ時間
内に、共通の原点に対して何れも異なる角度であ
つて、必要によつては、その組の内の少なくとも
１つの像の信号対雑音比を高めた複数個の像の
組、信号対雑音比を高めた１個の像、流れ全体の
造影像、動きの人為効果を一層抑圧した少なくと
も１つの像、又は何れも流速が異なる複数個の像
がこの順序によつて得られるような形で実施され
る。
従つて、この発明の目的は、生体の解剖学的な
部分のようなサンプルの所望の部分を通る血管の
ような管を流れる動く流体の多重エコー核磁気共
鳴造影作像法を提供することである。
この発明の上記並びにその他の目的は、以下図
面についてこの発明を詳しく説明するところから
明らかになろう。
発明の詳しい説明 最初に第１図について説明すると、この発明の
方法は、サンプル、例えば患者が搬送手段１１に
載せられて、作像領域内に運びこまれるような核
磁気共鳴（NMR）装置１０で実施される。この
領域は、磁石手段（図に示していない）の中孔１
２の中にあるように示されているが、この領域で
は高度に均質で比較的強度の強い静磁界B0が発
生される。図示の静磁界を形成する手段の中孔１
２は円筒形であり、静磁界B0が選ばれた容積の
軸線、例えばその中心が磁石手段の中孔の中にあ
るデカルト座標系のＺ軸に沿つている。主静磁界
を形成する手段の中孔１２に付設された磁界勾配
を形成する手段１２ａを使つて、静磁界に重畳さ
れる１組の略直交する磁界勾配を形成する。図示
のデカルト座標系では、勾配Gx，Gy，Gzは Gz＝∂B0／∂Z， Gy＝∂B0／∂Y， Gx＝∂B0／∂x， ということができる。静磁界B0及び勾配コイル
手段１２ａ（典型的には磁石手段の中孔１２の中
に配置される）によつて発生される磁界勾配Gx，
Gy及びGzの他に、少なくとも１つのコイル手段
１２ｂ（これも典型的には中孔１２の中に配置さ
れ、勾配コイル１２ｂよりも中孔の中で更に内側
に配置される）が（図面に示していない手段によ
つて）付勢されて、共鳴周波数又はラーモア周波
数ω＝γB0で回転する無線周波（RF）磁界B1を
発生する。ここでγは作像しようとする特定の原
子核種目に対する磁気回転比である。静磁界B0、
磁界勾配Gx，Gy，Gz及びRF磁界B1は、何れも周
知の手段により、周知の方法で作られる。NMR
作像装置の他の部分（図面に示してない）の他
に、NMR作像装置心臓（EKG）信号を供給する
モニタ導線手段１４を作像しようとする患者の胸
区域に設けることができる。前記の米国特許出願
に説明されているように、NMR造影法の順序の
心臓ゲート動作を使つて、周期的な動きによる人
為効果、即ち、ゴーストが最終的な像に現われな
いようにし、並びに／又は心臓サイクル内、造影
図を撮影する関心の持たれる特定の部分を選択す
ることができる。
第２図について説明すると、前に引用した米国
特許出願に記載されるNMR造影法は、各々の収
集期間に１個の造影投影を収集することができる
ようにしている。磁界勾配及びRF磁界パルス信
号の逐次的な各対の印加により、１つの投影軸線
Ａ（１つだけ）に沿つて見た、投影平面Ｐに投影
される像の多数（Ｓ個、典型的には128又は256
個）の平行な縞の内の１つに対し、選ばれた原子
核種目のスピンの動きの濃度に関する情報を持つ
相異なる１組の差のNMR応答信号が得られる。
このため、何れも相異なる投影軸線Anに沿つて、
相異なる投影平面Pnに対する複数個（Ｎ個）の
造影投影には、（Ｎ×Ｓ）個の逐次的な対のパル
ス順序を印加することが必要であり、その各々が
相異なる逐次的な、重ならない収集期間にある。
例えば、夫々の軸線A1，A2，A3及びA4に沿つ
て、（＋Ｙ軸に対して）θ1＝0°，θ2＝30°，θ3＝60
°
及びθ4＝90°の夫々の角度θで、Ｎ＝４の造影投
影情報を収集するためには、相次ぐ４回の収集期
間が必要である。典型的には各々の期間が２乃至
20分の範囲であるから、この手順は８乃至80分を
必要とし、動きによる人為効果を少なくするため
には、患者を１つの姿勢に比較的拘束しなければ
ならなくなる。
この発明の一面では、各々の励振順序に対する
複数個の、情報を持つエコー応答を収集する。
RFπ（180°）再集束パルス又は勾配磁界再集束パ
ルスの何れかを使つて、追加のエコーを発生して
検出する。追加の各々のエコーは、各々の順序で
最初のエコー応答で得られた情報を補う余分の情
報を提供するために利用される。偶数個の順序
（これは順序の対で極性が交互に変わる流れ符号
化勾配パルスを使うために必要である）の各々
が、心臓サイクル中の同じ点で開始されれば、順
序の合計の持続時間が、それ自体は患者の心臓サ
イクルによつて設定される、即ち、患者の脈率の
逆数に設定される順序の繰返し期間TR未満であ
る限り、任意の数（Ｎ個）の多重エコー応答信号
を形成して収集することができる。
この発明の多重エコーNMR造影法の種々の好
ましい実施例を説明する前に、スピン磁化の位相
を監視することにより、巨視的なスピンの運動を
作像するために用いられる理論を簡単に説明す
る。横方向のスピン磁化の位相変化に関する必要
な情報は、２モードの流れ符号化勾配パルスの
組、即ち、各々の励振順序で反対の極性を持つ一
対のパルスから容易に得られる。交互の順序では
極性が交互に変わる。磁界勾配が存在する時のラ
ーモア周波数ωは、応答信号を発生したスピンの
位置に関係する。即ち、ω（ｚ）＝γ（B0＋ｚ・
Gz）である。ここでGZは所望の方向、例えばＺ
方向の磁界勾配の強さである。時刻ｔ＝０ニ開始
して、別の時刻ｔ＝Tgに終る印加された勾配パ
ルスに応答して、横方向のスピン磁化の位相変化
φは次の式のようになる。
φ＝γ∫Tg 0ｚ（ｔ）Gz（ｔ）dt ……(1) ここでｚ（ｔ）及びGz（ｔ）は夫々時間の関数
としてのスピンの位置及び勾配磁界の強さであ
る。動く及び動かない両方の励振されたスピンが
勾配磁界の作用を受ける。動く励振されたスピン
だけが一定でないスピン位置時間関数ｚ（ｔ）を
持つ。励振されたスピンが、印加された磁界勾配
の方向に絶えず動く場合、このスピン位置時間関
数はｚ（ｔ）＝Z0−Vtである。ここでZ0は時刻ｔ
＝０に於けるスピン位置であり、Ｖはスピンの速
度である。このスピン位置時間関数を式(1)に代入
すると、動くスピンの位相変化は次のようになる
ことが分かる。
φ＝γ∫Tg 0Z0Gz（ｔ）dt −γTg 0VtGz（ｔ）dt ……(2) 最初の勾配パルスの始めから時間Ｔ後に２番目
の勾配パルスG′z（ｔ）をサンプルに印加すると、
この一対の勾配磁界パルスによつて導入された移
相φ全体は次の式で表わされる。
φ＝γ∫Tg 0Z0Gz（ｔ）dt −γ∫Tg 0VtGz（ｔ）dt ＋γ∫T+Tg TZ1G′z（ｔ）dt −γ∫T+Tg TVtG′z（ｔ）dt ……(3) ここでZ1は２番目の勾配パルスG′z（ｔ）が開始
する時のスピンの位置である。２番目の勾配パル
スの形及び振幅を１番目の勾配パルスの形及び振
幅と同一であるが、反対の極性にすると、即ち、
G′z（ｔ）＝−Gz及びT′g＝Tg′とすると、式(3)は次
のようになる。
φ＝γAg（Z0−Z1） ……(4) ここでAgは各々の勾配パルスの下にある合計
面積である。スピンの速度を一定と仮定したか
ら、Z0−Z1＝VT、従つてφ＝γAgVTである。
これが動くスピン磁化を選択的に検出するための
根拠になる。２モードの一対の勾配パルスによ
り、速度Ｖ＝０を持つ（定義により）不動のスピ
ンのスピン磁化に誘起される移相は、それ自身が
０であることが分かる。逆に、動くスピン磁化に
よつて誘起される移相は、スピン速度Ｖ、パルス
間の遅延期間Ｔ及び勾配パルスの面積（Ag）に
直線的な関係を持つことが分かる。従つて、これ
らの３つの変数（Ｖ，Ｔ又はAg）の内のどれか
の２つの異なる状態に対して収集されたデータを
互いに減算すれば、その結果は、動くスピンに対
するスピン密度情報だけを持つ。前に引用した出
願に記載されるように、NMR造影法は、一対の
内の一方の順序に２モードの勾配パルスが無い場
合でも作用するが、交互の励振で流れ符号化勾配
パルスの極性を反転すれば、更に良い結果が得ら
れる。これから説明する例のパルス順序ではこの
形を使う。
以上の簡単な理論的な説明から（上に述べたよ
うな）２つのエコーの相異なる１つから夫々収集
された２組のデータの複素数の差は、スピン速度
の正弦関数である係数を持つことが分かる。これ
は、あるスピン速度では観測される信号が得られ
ず、或いはそれより遅いスピン速度と識別するこ
とができないというようなエイリアシング状態の
可能性が生ずる。エイリアシングの問題は、誘起
される移相をπラジアン未満に制限すれば、避け
ることができる。実際、誘起された移相φが１ラ
ジアン未満に制限されれば、像の強度はスピン速
度に対してだいたい直線的である。測定される信
号の強度も、スピンの数の直線的な関数であるか
ら、像の画素の強度は容積（従つて、質量）流量
に比例し、このため、血管にわたつて信号の強度
を積分することにより、血管内の血液の合計流量
を測定することができる。このパルス順序によつ
て、流れ符号化勾配を印加した方向に於ける流れ
の成分だけが作像されるが、何れも互いに直交す
ると共に投影軸線に対しても直交する流れの方向
を感知するような２つの別々のNMR造影図を収
集し、式It＝（IA 2＋IB 2）1/2（ここでItは合計の強
度、IA及びIBは直交する２つのNMR造影図に於
ける対応する画像の強度である）を使つて、２つ
の個別のNMR造影図を組合せることにより、全
体的な流れのNMR造影図を求めることができ
る。最後に、任意の標準的な作像手順に流れ符号
化勾配パルスを用いることにより、NMR造影図
を発生することができることを承知されたい。こ
こで利用する作像手順は勾配再集束形スピン捩れ
形順序であるが、スピンエコー再集束形スピン捩
れ形等の手順も用いることができる。
多数の視角（angle of view） 投影造影図の投影軸線Ａは、作像装置の座標系
の主軸の内の選ばれた第１の軸（例えば、＋Ｙ）
と軸線Ａの間に形成される角度θと、この第１の
軸線に対して直交する座標系の別の軸（例えば＋
Ｘ）と軸線Ａの間に形成される角度θとによつて
限定される。角度θ及びφは、読取勾配磁界及び
位相符号化勾配磁界の向きによつて決定される。
読取用及び位相符号化用勾配磁界は直交してい
る。作像の多くの用途では、直交する読取用及び
位相符号化用の勾配は、(1)パルス状であり、(1)デ
カルト座標軸の相異なる軸に沿つて印加される。
座標系の軸に対して斜めのビユー（view）は、
何れも座標系の２つ又は３つの軸に沿つた成分を
持つと共に、三角関数から容易に導かれる振幅を
持つ読取用勾配パルス及び／又は位相符号化用勾
配パルスを印加することによつて求めることがで
きる。例えば、＋Ｙ軸に対して角度θk（並びにＺ軸
に対してφk＝90°）に設定された投影軸線Akに対
して直交するある平面Pkに投影された造影図で
は、Ｘ軸読取用勾配成分はGx＝cos（θk）Gx,naxで
あり、Ｙ軸の読取成分はGx＝sin（θk）Gy,naxでる。
これによつて、作像する被検体を動かすことを必
要としないで、投影角度θを変える時、一定の分
解能でNMR造影像を求めることができる。
この発明の一面では、各NMR造影順序対の
夫々の作像順序の間、同じ複数個（Ｎ個）の読取
期間Tkの関連した、何れも異なる視角θk（１≦ｋ
≦Ｎ）を持つ複数個（Ｎ個）のNMR造影図の各
対に対するデータをインターリーブした形で受取
る。多数のビユーの多数の読取データの組を求め
るための一般的な形のRF／勾配パルス順序が第
２ａ図に示されている。この一般的な順序は、そ
の順序の１個のRFパルスの直後に一対の流れ符
号化勾配パルスを印加することを必要とし、この
後の複数個の多重エコーの各々は、読取勾配の異
なる向きを用いて得られる。具体的に云うと、多
重視角を求める方法は、対の励振応答順序を使
い、第２ａ図の順序を使い、その後は、１番目の
順序の破線部分に示すように、２番目の順序で流
れ符号化勾配の極性を交互に変えるが、これを後
で説明する。必要な場合、各々の読取勾配磁界パ
ルスの前に位相外しパルス部分がくるようにし
て、関連するエコーが関連する読取パルスの中央
で発生するようにすることができる。その後、
各々の読取勾配磁界パルスの後に位相戻しパルス
部分を印加して、関連するエコーの発生中に累積
した移相を除くことができる。このため、最初に
サンプルに対し位相不変のRF信号パルス２０を
印加して、スピンを90°にわたつて章動させる。
一般的にπ／2RFパルス２０は容積選択性にはし
ない。しかし、少なくともパルス２０と同じ期間
（時刻t0から時刻t1まで）にわたつてスライス選
択性勾配Ｇ（スライス）部分２１を用いることに
よつて、パルス２０をスライス選択性にし、サン
プルの励振される容積を制限することができる。
周知のように、スライス選択パルスを使う場合、
RFパルス２０は、例えばsin（ｘ）／ｘの形にし、
励振される容積の限界を実質的に鋭くするため
に、特定の奇数個（例えば５個）のローブを持つ
ように切取ることができる。やはり良く知られて
いるように、スライス選択性勾配部分２１を使う
場合、時刻t1から時刻t2までの期間内の関連する
位相戻し部分２２を使うのが有利である。その
後、時刻t3から開始して、流れ符号化勾配Ｇ（流
れ）には、持続時間Tg（時刻t3時刻t4まで）、第１
の極性及び第１の振幅Ａを持つ第１のパルス２４
と、第２の持続時間T′g（時刻T5から時刻t6ま
で）、第２の極性（第１の極性に対して反対）及
び第２の振幅A′を持つ第２のパルス２５を用い
る。パルス２４及び２５の中点はパルス間間隔Ｔ
だけ離れている。具体的に云うと、選ばれた順序
の対は、第１の順序に、正の極性の第１のパルス
２４および負の極性の第２のパルス２５を用い
る。（この順序の対の内、第２の順序の破線で示
したパルスは、負の極性の第１のパルス２４′及
び正の極性の第２のパルス２５′を持ち、第２の
順序のパルスが第１の順序のパルスに対して極性
が反転するようにする。）両方のパルス２４及び
２５（並びにパルス２４′及び２５′）が同じ振幅
｜Ａ｜及び同じ持続時間Tg＝T′gを持つのが有
利である。時刻t6の後、何れも相異なる投影軸線
の１組のデータを供給する多重エコーの読取りを
行なう。１つの投影角度θk並びにそれに関連する
投影補角φk（ここで（θk−φk）＝90°）を持つ各々
の読取りの図に対し、位相外し読取勾配Ｇ（読取）
部分３１と略同時に、位相符号化勾配Ｇ（位相符
号化）パルス３０を印加する。スピンの位相に対
する各々のパルス３１の影響は、関連する、極性
を反転した位相戻し読取勾配部分３５の影響によ
り、この後で実質的に相殺される。投影位相外し
勾配パルス３２を利用して、検出された信号のダ
イナミツク・レンジを制限する。スピンの位相に
対する各々パルス３２の影響は、この後、極性を
反転した関連する位相戻し投影勾配パルス３６の
影響によつて、実質的に相殺される。更に、偶数
番号のエコーの各々の位相符号化勾配パルス３０
ｂ，３０ｄ……は、奇数番号のエコーの位相符号
化勾配パルス３０ａ，３０ｃ……の極性に対して
反転して、一対の位相符号化勾配によつて起こ
る、速度によつて誘起される移相を相殺する。こ
のため、（第１の角度θ1で第１の軸線A1に沿つた
１番目の投影に対する情報を持つ）第１のエコー
に対し、（周知のように、複数個の位相符号化の
値の内の１つの値を持つ）位相符号化勾配部分３
０ａを、（角度θ1の方向の）第１の読取勾配部分
の位相外し部分３１ａ及び（角度θ1に直交する補
角φ1の方向の）第１の投影位相外し勾配パルス
３２ａと略同時に用いる。勾配部分３０ａ，３１
ａ及び３２ａが落着くのを待つための短い期間の
後、第１の読取期間のNMR応答データを収集す
るため、第１の読取符号化勾配部分３４ａの存在
のもとにデータ・ゲートを開く（部分３３ａにあ
る）。その後、データ・ゲートを閉じる。読取勾
配３４ａが落着かせるための別の短い期間の後、
第１の読取勾配位相戻しパルス部分３５ａが発生
して、スピン磁化の位相戻しをする。位相戻し読
取パルス部分３５ａは、最初の投影位相外しパル
ス３２ａによつて起こるスピンの位相外し効果を
相殺するために、第１の補角の軸線φ1に沿つて
印加される、極性を反転した第１の投影位相外し
パルス部分３６ａと略同時に発生する。パルス３
５ａ及び３６ａと同時に要求される、反転した位
相符号化勾配部分はこの図に示していないが、こ
れは、この時負の極性を持つこのパルスが、偶数
番号の第２の読取期間の始めにある反転した極性
（例えば、この時負の極性）の位相符号化勾配部
分に相加つて、２倍振幅の位相符号化パルス３０
ｂを形成するからである。即ち、１つの読取期間
の終りにある反転した位相符号化パルス部分は次
の読取期間の始めにある最初のパルス部分と同じ
極性であり（これは交互のエコーに対し、位相符
号化勾配の極性を交互に変えるから）、このため、
位相符号化パルス３０ｂは２倍の負の振幅であ
る。同時に、第２の投影角θ2の勾配読取位相外し
パルス部分３１ｂが第２の投影位相外しパルス部
分３２ｂと共に発生する。従つて、１回の読取り
の終りにある最終パルス部分を、次のエコーの読
取りの開始時にある準備パルス部分と重なるよう
にする。この時間的な重なりが、次の心臓サイク
ルが次の順序をトリガする前に残つている時間内
に、各々の順序の終りに更に多くの数の異なるビ
ユーの読出しを逐次的に行なうことができるよう
にする。第２のNMRエコー信号のデータは、(1)
読取勾配部分３４ｂによつて符号化され、(2)パル
ス部分３０ｂ，３１ｂ及び３２ｂが終了してから
短い時間後に受信され、(3)その後、データ・ゲー
トによつて付能された部分３３ｂの間に収集され
る。その後、第２のビユーの読取りは、略同時の
位相戻し勾配パルス部分３５ｂ及び３６ｂと、今
度は正の極性の位相符号化勾配部分３０ｃで終了
する。同時に、第３の多重視角エコー読取部分
が、位相符号化パルス部分３０ｃの正の極性を持
つ半分、読取勾配位相外し部分３１ｃ及び投影位
相外しパルス３２ｃと共に開始する。読取後勾配
部分３４ｃによつて角度θ3に符号化された第３の
ビユーの読取りは、ゲート部分３３ｃの間に収集
され、その後位相戻しパルス部分３５ｃ及び３６
ｃ及び位相符号化パルス３０ｄの半分を用いて終
了する。この手順が、何れも関連した異なるビユ
ーのデータを得るために、必要なだけ多くの追加
のエコーに対して続けられる。従つて、（図面の
示すように）４つのビユーしか必要としない場
合、４番目の、即ち最後のビユーのエコーの読取
りは、位相符号化パルス部分３０ｄ、読取位相外
しパルス３１ｄ及び投影位相外しパルス部分３０
ｄの関連した半分によつて準備される。第４のビ
ユーのデータは、データ・ゲートの付能部分３３
ｄの間、第４の（θ4）読取勾配部分３４ｄが存在
する間に収集される。最後のエコーの読取勾配部
分３４ｄに応答して、最後のデータ・ゲートの付
能部分３３ｄの間の最後のエコー・データを収集
した後、この順序が終了する。追加のビユーから
のデータを求めないから、余分の読取位相戻しパ
ルス部分３５又は極性を反転した投影位相外しパ
ルス部分３６を必要としない。同様に、最後のエ
コーに対する最初のパルス符号化勾配部分３０ｄ
の影響は問題ではないから、最後のエコーの読取
りが終了する時、位相符号化勾配の極性を反転し
た部分は必要ではない。一対の内の２番目の順序
も略同様であるが、流れ符号化パルス２４′及び
２５′の反転だけが１番目の順序の特性とは異な
る。
読取勾配の向きの角度θi、投影位相外し勾配の
向きφi及び位相符号化方向は全て互いに直交して
いる。従つて、原則として、各々のエコーに対す
る位相符号化勾配の向きは、そのエコーに対する
読取勾配の向きと同じように変えることができ
る。位相符号化勾配パルスを主座標系の１つの
軸、例えばデカルト座標系のＺ軸に沿つて印加す
る場合、各々のエコーに対する位相外しローブは
読取（Ｒ）位相外し部分及び投影（Ｐ）位相外し
部分を含む。投影角度θnを用いて収集されるＮ
番目のエコーでは、Ｘ軸勾配コイル及びＹ軸勾配
コイルに同時に夫々印加されるパルスの振幅Ix及
びIyは次の通りである。
Ix＝IRsin（θN）＋IPcos（θN） ＋IRsin（θN-1）−IPcos（θN-1） ……（5a） Iy＝IRcos（θN）−IPsin（θN） ＋IRsin（θN-1）−IPcos（θN-1） ……（5b） 第２ａ図の種々の勾配が、デカルト座標系の主
軸に沿つた勾配に対して正しく割当てられた時、
第２ｂ図に示すような１組の波形が得られる。
第２図に示す例の場合、θ＝0°，θ2が大体30°、
θ3が約60°及びθ4を90°にして、４つの投影造影図
を求める。順序の対の内の各々の順序は、スライ
ス選択性部分Gs（第２ａ図に示したＧ（スライス）
部分に対応する。）で始まり、これは(a)時刻t0か
ら時刻t1までの、90°RFパルス２０′及びそのスラ
イス選択性Ｚ軸勾配Gz部分２１′及び(b)時刻t1か
ら時刻t2までの期間内のスライス選択性位相戻し
勾配部分２２′で構成される。その後、第２ａ図
のＧ（流れ）部分に対応する流れ符号化部分Gfが
Ｚ軸勾配Gz中に発生し、これは振幅＋Ａ並びに
持続時間Tg（時刻t3と時刻t4の間）を持つ第１の
正の極性を持つパルス２４ａと、振幅−Ａを持ち
時刻t5からt6までの期間Tg内に発生する負の極性
のパルス２５ａとを用いる。パルス２４ａ及び２
５ａの時間的な中点はパルス間期間Ｔだけ離れて
いる。その後、４つの多重のビユーの多重エコー
受信期間T1乃至T4が発生する。その各々は、開
始時刻ta1〜ta4に開始し、終了時刻te1（第２の期間
の開始時刻ta2と同一）乃至te4（この後の続くエコ
ーの応答データを必要とする場合の５番目のエコ
ー開始時刻ta5と同一になる）で終る。ある一対
の順序に対して利用される、複数個の値の内の１
つの値で、各々のエコー収集期間T1乃至T4に位
相符号化勾配パルス３７ａ乃至３７ｄの内の１つ
を利用する。関連する応答期間が開始してから短
い時間後に、時刻tb1から時刻tc1までの期間内に、
夫々Gy及びGk勾配中に読取位相外しパルス部分
３８ａ及び３９ａ（後者は今の場合は振幅ゼロで
ある）が発生する。各々のパルス３８ｎ及び４０
ｎは振幅Gy＝Gy,nax（sinθo）を持ち、各々のパル
ス３９ｎ及び４１ｎは振幅Gx＝Gx,nax（cosθo）を
持つ。その後、時刻td1から時刻te1まで、Ｙ軸及
びＸ軸読取勾配部分４０ａ及び４１ａ（後者は今
の場合振幅ゼロでる）が発生して、第１の順序の
第１の投影読取データを符号化する。この期間の
間、データ・ゲートの付能信号４２ａが存在し、
第１の投影のNMR応答エコー信号４３ａを収集
して処理する。次の（第２の）ビユー読取期間
T2が開始されてから間もなく、時刻tb2から時刻
tc2までの期間内に、（振幅が２倍で極性を反転し
た）次の位相符号化パルス３７ｂと共に、読取位
相外しパルス部分３８ｂより３９ｂが夫々Gy及
びGx勾配中に発生する。その後、時刻td2から時
刻te2まで、Ｙ軸及びＸ軸読取勾配部分４０ａ及
び４１ａが発生して、第１の順序の第２の投影の
読取データを符号化する。この期間の間、デー
タ・ゲート付能信号４２ｂが存在し、第２の投影
NMR応答エコー信号４３ｂを収集して処理す
る。次の（第３の）ビユーの読取期間T3が開始
されてからまもなく、時刻tb3から時刻tc3までの
期間内に、位相符号化パルス３７ｃ及び読取位相
外しパルス部分３８ｃ及び３９ｃが夫々Gy及び
Gx勾配中に発生する。時刻td3から時刻te3まで、
Ｙ軸及びＸ軸読取勾配部分４０ｃ及び４１ｃが発
生して、第１の順序の第３の投影の読取データを
符号化する。この期間の間、データ・ゲート付能
信号４２ｃが存在し、第３の投影のNMR応答エ
コー信号４３ｃを収集して処理する。最後の（第
４の）ビユーの読取期間T4が開始されてからま
もなく、時刻tb4から時刻tc4までの期間内に、位
相符号化パルス３７ｂ及読取位相外しパルス部分
３８ｄ及び３９ｄ（前者の読取位相外しパルス部
分は今の場合振幅ゼロである）が夫々Gy及びGx
勾配中に発生する。その後、時刻td4からte4まで、
Ｙ軸及びＸ軸読取勾配部分４０ｄ及び４１ｄ（前
者は今の場合振幅ゼロである）が発生して、第１
の順序の第４の投影の読取データを符号化する。
この期間の間、データ・ゲート付能信号４２ｄが
存在し、第４の投影のNMR応答エコー信号４３
ｄを収集して処理する。第１の順序の開始時刻t0
から繰返し期間TRが経過した時刻t′0に、第２の
順序が開始される。第２の順序のGf部分だけが
第１の順序の同様な部分とは異なる。流れ符号化
パルス２４′ａ（負の極性、振幅Ａ）は、正の極性
の振幅Ａのパルス２５′ａより前に発生する。パ
ルス３７の同じ「縞形」の符号化の値が存在し、
勾配Gy中の同じパルス３８及び４０と勾配Gx中
のパルス３９及び４１の全てが発生する。典型的
なNMR造影検査では、患者の脈搏数が毎分60乃
至120であると、繰返期間TRは約0.5秒乃至約1.0
秒になる。各々の測定の間の一貫した速度を保証
する為に心臓との同期を用いた場合でも、その期
間全体に亘る、流れる血液のT* 2によつて、後の
方のビユーの品質が幾分低下するけれども、５mm
秒程度のビユー標本化期間Tnにより、典型的に
は10個までの相異なるビユーに対するデータを求
めることができる。
このモードの変形では、各々に対する視角θn
を変更せず、読取位相外し部分４６ａ乃至４６ｄ
（Gyに対する）は何れも同じ振幅であり、部分４
７ａ乃至４７ｄ（Gxに対する）も同様である。各
組の読取勾配４８ａ乃至４８ｄ及び４９ａ乃至４
９ｄは一つの振幅で同じ極性である。ｊ番目のエ
コーに対する投影位相外しパルスは、（ｊ−１）
番目のエコーに対する投影位相戻しパルスによつ
て実質的に相殺することができる。しかし、式(5)
に従つて、ｊ番目のエコーの読取位相外しパルス
と（ｊ−１）番目のエコーに対する読取位相戻し
パルスは相加わる。第２ｃ図に示すように、交互
の読取勾配の極性を逆にすれば、こう云う読取位
相戻し／位相外しパルスを全く省略することがで
きる。即ち、RFパルス２０（並びに使う場合に
はスライス選択性勾配パルス２１及び２２）とそ
の対の流れ符号パルス２４ａ／２５ａ（各対の第
１の順序に対する）又は２４′ａ／２５′ａ（各対
の第２の順序に対する）の後、読取部分が、何れ
も第１の位相符号化パルス３７ａと略同時の１個
のＧ（読取）位相外しパルス５０及び１個のＧ（投
影−位相外し）パルス５１と共に開始される。読
取勾配にも投影勾配にも、この他の位相外し又は
位相戻しパルスは現われない。夫々データ・ゲー
トによつて付能された期間４２ａ，４２ｂ，４２
ｃ及び４２ｄと一致する極性が交互に変わる読取
勾配パルス部分５２ａ，５２ｂ，５２ｃ及び５２
ｄで十分である。この様な一般的なパルス順序を
使つて、一連のＭ個の同一のNMR造影図を収集
し、それを加算して信号対雑音比を高めることが
できる。第２ｃ図では、Ｍ＝４のデータ収集の全
部が同じ投影角度であるが、一つの角度θ1で第１
の数M1個の読取を行ない、別の角度θ2で別の数
M2個（これはM1と同じである必要はない）を行
なうということも可能である。これらの方式を混
合させることにより、各々の順序の対にある各々
の順序で、インターリーブ形のMi個の逐次形及
び非逐次形の読取を投影角度θiで収集し、Mj個
の逐次的及び非逐次的な読取りを角度θjで行な
い、Mk個の逐次又は非逐次形読取りを角度θk等
へ行なうことができる。従つて、Ｍが１より大き
ければ、多重角度の１組の投影でも、或る投影角
度に対して、改善された信号雑音比を利用するこ
とができる。
直交する流れ方向の像の同時の収集 前に述べた様に、個別の各々のNMR造影図に
対する復元された組のデータは、一方向だけの１
組の流れ成分を持つている。従つて、互いに直交
する方向Ａ及びＢで撮影した一対のNMR造影図
を組合せて、各々の画素が次の式で表わされる、
合計強度Itを持つ様な全体的な流れのNMR造影
図を作られなけれならない。
It＝（IA 2＋IB 2）1/2 ……(6) ここでIA及びIBは夫々の方向Ａ及びＢにおける
流れの強度である。一対のNMR造影図の夫々に
対し、何れも多数（128個又は256個）の順序の次
いで構成された従来の単一エコー形のNMR造影
パルス順序は、逐次形（すなわち、方向Ａの
NMR造影図に対する128個又は256個の必要な順
序の対全部を収集し、その後方向ＢのNMR造影
図に対する必要な、例えば128個又は256個の順序
の対の全部を収集する形式）で収集しても良い
し、又はインターリーブ形（すなわち、１番目の
順序の対では、第１のＡの方向で収集し、その後
１番目の順序の対でＢ方向の収集をし、その後必
要な数、例えば128個又は256個の順序の対の内の
２番目、３番目……等の順序の対に対し、逐次的
に方向Ａ及びＢで順序の対の収集をする形式）で
収集しても良い。例えば１番目の像は、読取方向
（例えばＸ軸）に沿つて、流れの感度を持つ様に
収集し、２番目の像はパルス符号化方向（Ｚ軸に
沿つて）流れの感度を持つ様に収集することがで
きる。どの対の軸を選んでも、逐次的な２つの単
一方向の流れのNMR造影図で構成しても、或い
はインターリーブ形の同時の２つの単一方向の流
れのNMR造影図で構成しても、全体的な流れの
NMR造影図は、単一の流れのNMR造影図の２
倍の収集時間を必要とする。この欠点を解決する
為、第３ａ図に示す一般的なインターリーブ形多
重エコー方式を利用することにより、２つの流れ
の成分を略同時に収集することにより、２つの流
れ成分を持つ全体的なNMR造影図を作ることが
できる。この順序は、各々の縞に対し、２対の順
序を必要とする代りに、全体的な像に対して必要
な（例えば128又は256個の）全部の「縞」の夫々
１つに対する一対の順序で、方向Ａ及びＢの両方
の流れを感知するエコーを収集する。
この発明の別の一面として、全体的な流れの
NMR造影図を表示するのに必要な、互いに直交
する一対のNMR造影図の各々に対するデータ
は、各NMR造影順序対の内の各々の作像順序の
間、一対の読取期間の夫々１つでインターリーブ
形で収集する。第１のエコーに対して印加される
流れ符号化勾配の後、第１のエコー応答データの
読取りの後に、第１の流れ符号化勾配に対して極
性を反転した一対の流れ符号化勾配が続く。一対
の第２の勾配の内の一方は、第１の流れ符号化勾
配パルスの対によつて生じた速度によつて誘起さ
れる移相を除くために、第１の勾配方向に、極性
を反転した第２の流れ符号化勾配パルスの対を持
つている。他方の対の流れ符号化勾配パルスは、
第１の対の方向に対して略直交する第２の流れ符
号化方向に印加する。この後に第２のエコー応答
データを、第１のエコー応答データの流れ感知方
向に対して直交する流れ感知方向から受信する。
直交する流れ感知方向を、夫々読取及び位相符号
化勾配の方向と一致する様に割当てて、IA＝I〓及
びIB＝I〓となる様にすることが、必要ではないが
便利であることが理解されよう。即ち、読取勾配
をデカルト座標系のＸ軸に印加し、位相符号化勾
配をそのＺ軸に印加した場合、第１の流れ勾配は
Ｚ軸方向に印加することができ、第２の流れ勾配
はＸ軸方向に印加することができる。インターリ
ーブ形で収集された２つのエコーを独立に処理し
て、流れの成分に対するNMR造影図の強度の組
を求め、それらをこの後式(6)を用いて組合せる。
共鳴のオフセツト等による第１及び第２のエコー
NMR造影図の間の移相は、加算されるのが各々
のNMR造影図の大きさであるから、重要ではな
いことが理解されよう。更に、エコーが時間的に
比較的接近して（典型的には各々の順序のデータ
に対して、互いに10mm秒程度）収集されるから、
エコーの対のT2 *信号の強度の差は比較的小さい
ことが理解されよう。更に、この順序は第３のエ
コー又は場合によつては第４のエコーも含む様に
拡張することもでき、これらのエコーを処理し、
第１及び／又は第２のエコーの処理済みエコーと
比較して、T2 *強度の減衰を補正するのに役立つ
係数を導き出すことができることも理解されよ
う。
直交する流れ方向の像を同時に収集する為の一
般的な形のRF／勾配パルスの順序が第３ａ図に
示されている。スライス選択性勾配パルス２１及
び位相戻しパルス２２を伴なう、又はスライス選
択性勾配を（希望に応じて）伴なわない章動用の
90°RF信号パルス２０が、時刻t0及びt2の間の予
想期間内に発生する。第１の対の流れ符号化勾配
パルス２４ａ及び２５ａ（これは勾配Ｇ（流れ１）
にある。）が時刻t3及び時刻t6の間に発生する。
その後、第１の読取期間が、投影位相戻し勾配パ
ルス５６及び第１の位相符号化勾配パルス３７′
ａと略同時に、前の位相戻しパルス部分５５ａか
ら始まる。関連する読取用勾配パルス部分５７ａ
が存在する第１の読取期間の中央で発生する様に
した関連するエコーを、データ・ゲートの付能部
分５９ａの間に受信する。その後、読取勾配に位
相戻しパルス６０を用いて、第１のエコーを発生
する間に累積した移相を取除く。読取期間はパル
ス６０が終了した時に終る。極性を反転した位相
符号化勾配パルス３７′ｂが位相戻しパルス６０
と略同時に現われて、パルス３７′ａの為に位相
符号化方向に生じた移相を補正する。パルス６０
及び３７′ｂが終つてから短い時間の後、パルス
２４ａ，２５ａの第１の流れ符号化勾配の対（又
は一対の内の第２の順序に於ける極性を反転した
パルス２４′ａ／２５′ａ）の効果は、第１の流れ
勾配の極性を反転した第２の流れ符号パルスの対
によつて位相が補償される。この対のパルス２４
ｂ，２５ｂ（又は一対の内の第２の順序に於ける
極性を反転したパルス２４′ｂ，２５′ｂ）は、
夫々同じ期間Tg（夫々時間t′3とt′4の間又は時刻t′5
とt′6の間）に、第１のパルス対の期間Ｔと略等
しいパルス間期間T′で発生する。略同時に、第
２の直交する勾配Ｇ（流れ２）に一対の流れ符号
化パルスを用いる。この対のパルスの極性は、第
１の流れ符号勾配方向に於ける第１の流れ符号化
パルス２４ａ及び２５ａ（又は２４′ａ及び２５′
ａ）の極性と同じである。即ち、第１の順序で第
１の正の極性のパルス６２の後に負の極性のパル
ス６４が続く（又は一対の第２の順序では負の極
性のパルス６２′の後に正の極性のパルス６４′が
続く）。パルス６２及び６４によつて第２の流れ
方向を符号化した後、先行する位相外しパルス５
５ｂの発生により、第２のエコーの読取が準備さ
れる。読取勾配部分５８ｂによつて読取方向に符
号化された直交する流れ方向のデータが、デー
タ・ゲートの付能部分５９ｂの間に収集される。
流れ符号化パルス（例えば、パルス２４′ａ／２
５′ａ，２４′ｂ／２５′ｂ及び６２′／６４′）を
反転してこの順序全体を繰返した後、第１の順序
の第１のエコーからのデータを、第２の順序の第
１のエコーからのデータから減算して、第１の方
向の像データを取出す。第１の順序の第２のエコ
ーからのデータを第２の順序の第２のエコーから
のデータから減算して、直交する方向の像に対す
る１組のデータを発生する。その後、２組のデー
タを式 Cij｛（Aij）2＋（Bij）2｝1/2 を用いて、２次元フーリエ交換して加算する。こ
の式でCijは２次元像マトリクス内の点ｉ，ｊの
最終データであり、Aijは各々の順序の第１のエ
コー応答に応答して発生された像内の同じ点ｉ，
ｊに対して得られた１組のデータであり、Bijは、
各々の順序の第２のエコー応答から発生された像
内の同じ点ｉ，ｊに対する１組のデータである。
デカルト座標系に於ける特定の場合が第３ｂ図
に示されており、この図は、流れ符号化パルスの
対は、一対の内の両方の順序で反転しても正しく
動作すること、並びに「正常」な順序の対（即
ち、第１の順序の対が第１の正の極性のパルスを
持つ）の第２の順序の流れ符号化パルスの対をも
示している。スライス選択性Gz勾配パルス２
１′並びにそれに関連する位相戻しパルス２２′と
共に又は伴なわずに、RF信号パルス２０′の後、
第１の方向に対して選択性の（例えばＺ方向の）
流れ符号化用の１組の勾配パルスGfzが、時刻t3
及びt4の間の第１の負の極性のパルス６６ａ、及
び時刻t5及びt6の間の第２の正の極性のパルス６
８ａを含む。それより若干後の時刻t7に、位相符
号化勾配Gpe部分が順序の対のデータの組にある
多数の縞の内の１つに対する勾配パルス７０ａ
の、多重の値の内の１つの値と共に開始する。そ
れと同時に、時刻t7及びt8の間の期間に、Gx信号
中のダイナミツク・レンジを制限する投影位相外
しパルス７４と共に、読取位相外しパルス７２ａ
が発生される。位相符号化、投影位相外し及び読
取位相外しパルスの後、Ｙ方向読取勾配部分７６
ａが時刻t9からtaまで存在し、この期間中に、デ
ータ・ゲート７８ａを付能し、受信応答信号８０
ａを収集し、デイジタル化し、特定の装置によつ
て要求される通りに処理する。その後、読取勾配
Gy中の位相戻しパルス８２が時刻tbに開始し、時
刻tcに終るが、極性を判定した位相符号化相殺用
パルス７０ｂも同じである。その少し後、一対の
第１の方向（Ｚ方向）の相殺用パルス６６ｂ，６
８ｂが、時刻td及びtgの間の期間に、相殺勾配
G-fzに対して発生し、これと略同時に（第２の流
れ勾配Gfyに対する）第２の方向（Ｙ方向）の流
れ符号化パルス８６及び８８が発生する。Ｚ方向
の位相変化パルス６６ｂ／６８ｂによつて補正さ
れ、流れ符号化がＹ方向に付能されているので、
時刻thから時刻tiまでの期間に、Ｙ軸の読取予備
位相外しパルス７２ｂと略同時に新しい位相符号
化パルス７０ｃが発生する。その後、時刻tjから
時刻tkまでの期間に、読取勾配部分７６ｂが存在
し、データ・ゲート７８ｂが開いて、第１の応答
信号８０ａのＺ方向に対して直交する方向の、第
２のＹ方向の応答信号８０ｂを受入れる。Ｚ及び
Ｙ投影データの組をフーリエ交換し、式(6)の操作
にかけて、必要に応じて表示する為に、全体的な
NMR造影図のデータの組として記憶する。
動きによる人為効果の強い抑圧 この発明の別の特徴として、振幅ゼロの部分８
４ａ及び８４ｂで示す様に、第１の一対の流れ符
号化パルス６６ａ／６８ａの振幅をゼロにした場
合、第３ｂ図に示す信号順序が各々の順序の対の
内の第１の順序（前に説明した様に第２の順序で
はなく）であり、この時、動きによる人為効果が
比較的強く抑圧された、単一の流れ方向の１つの
像が得られる変形になる。即ち、図示の抑圧の強
い順序は、(a)各対の第１及び第２の順序の何れに
於いても、第１の読取データ・ゲート７８ａより
前に現われる振幅ゼロの部分８４ａ，８４ｂ、及
び(b)の第１の順序の第１及び第２のデータ・ゲー
トの間の正の極性の第１のパルス６６ｂ及び負の
極性６８ｂ、又は各対の第２の順序に於ける第１
及び第２のデータ・ゲートの間の負の極性のパル
ス６６ｃ（破線）及び正の極性のパルス６８ｃ（破
線）を持つていて、両方のエコーを同じ視角で求
める勾配再集束形スピン捩れ形順序になる。第１
のエコー応答データ信号８０ａは流れによつて左
右される位相変調を受けないが、第２のエコー信
号８０ｂ（これは一対の流れ符号化パルス６６ｃ，
６８ｃより後に収集される）は流れに影響されて
変調されているから、第１の順序の第２の応答信
号８０ｂのデータの組S12を第１の順序の第１の
応答信号８０ａのデータ組S11から減算すれば、
各々の順序に対する不動のスピンを相殺したデー
タの組S1＝（S11−S12）が得られる。順序の対の
内、第２の順序の第２のデータの組（S22）及び
第１のデータの組（S21）の差で不動のピンが相
殺される結果、別のデータの組S2＝（S21−S22）
が得られ、これを第１の順序のデータの組S1から
減算して、不動のスピンを更に抑圧した最終的な
データの組Ｓ＝S1−S2が得られる。表示される像
は、収集されたエコーの組の間の差の差を求める
ことにより、位置ｉ，ｊにある画素に対する特定
のデータの値が形成される。即ち、Eij＝（Aij−
Bij）−（Cij−Dij）である。ここでEijは２次元デー
タ・マトリクスに於ける点ｉ，ｊのデータであ
り、Ａ，Ｂ，Ｃ及びＤは、第１の順序の第１のエ
コー、第１の第２のエコー、第２の順序の第１の
エコー及び第２の順序の第２のエコーの夫々のデ
ータ点であり、その何れもフーリエ交換データ処
理が完了した後に求められる。第２及び第３の位
相外し勾配パルス８６，８８は、第２の読取勾配
パルス７６ｂが極性を反転したパルス７６ｃとし
て使われれば、省略することができることが理解
されよう。更に、第１の順序が第１の流れ符号化
パルス６６ｄ／６８ｄ及び第２の流れ符号化パル
ス６６ｃ／６８ｄを持ち、第２の順序が第１の流
れ符号化パルス６８ａ／６８ａ及び第２の流れ符
号化パルス６６ｂ／６６ｂを持つ順序の対に対し
ても、同様な効用があることが理解されよう。
多重速度の像 この発明の更に別の一面として、第３ｂ図のデ
カルト座標の順序に示してあるが、別の変形は、
流れ及び流れ符号化勾配パルスによつて応答信号
中に誘起される位相変化の相加性を利用して、
別々の１組の順序の対で、別々に各範囲を収集す
ることを必要とせず、速度の異なる複数個の範囲
の各々に対する１組の像データを発生することが
できる。一般的に言うと、Ｎ個の範囲の異なるエ
コーの各々より前に、Ｎ個の同一の対の流れ符号
化勾配パルスの内の１つが発生する。全体的な位
相変化は相加的であり、Ｎ番目のパルスの対の後
のＮ番目のエコー中の合計の位相変化は、同じ流
速に対する第１のエコーに於ける位相変化のＮ倍
の大きさである。従つて、Ｎ番目のエコーに対し
ては、検出された速度を係数Ｎで割る。即ち、第
１のエコーに対する合計の移相が、流速V1に対
応してφ1であれば、速度V1に対するＮ番目のエ
コーの移相φNはＮ・φ1であり、基本的な位相量
φ1に対する流速はV1／Ｎである。一対の内の第
１の順序の第１のエコーを収集（時刻t9から時刻
taまで）する前に、正の極性を持つ第１のパルス
６６ｄ（鎖線）、及び負の極性を持つ第２のパルス
６８ｄ（鎖線）が（Ｎ＝１に対して）発生する。
その後、これらより後の、最初は正の極性のパル
ス６６ｂ及び極性を反転した第２のパルス６８ｂ
の同じ極性の対に応答して、時刻tjから時刻tkま
での間に、第２のエコー（Ｎ＝２）が収集され
る。夫々３番目……Ｎ番目の対の最初は正の極
性、その後の負の極性のパルスを使つて、３番目
……Ｎ番目のエコーを収集することができる。
各々の順序の対（各対は異なる値の位相符号化パ
ルス７０を持つている）の第２の順序では、第１
の対のパルスは負の極性のパルス６６ａ及び正の
極性のパルス６８ａで構成され、第２の対の負の
パルス６６ｃ及び正のパルス６８ｃで構成される
という様になる。こうして、Ｎ個の像を収集し、
後の各々の像は異なる一層低い流速を感知する。
像データの組を「３次元」で（即ち、応答エコー
信号の数Ｎを関数として）フーリエ交換すれば、
真の速度選択性NMR造影図が得られる。
この発明を現在好ましいと考えられる幾つかの
実施例について説明したが、当業者には種々の変
更が考えられよう。例えば、不動のスピンによつ
て発生される信号の大きな振幅を少なくする為
（こうして、最終的なNMR造影図に於ける不動
スピンの残留成分を少なくする為）に、投影方向
のスピン磁化を位相外しする為に、勾配信号を用
いているが、関連する血管の物理的な場所の為
に、所望の応答信号の位相相殺が起ると思われる
場合、投影位相外し勾配を省くことができる。こ
の為に起る人為効果は、認識するのが割合い容易
である。同様に、位相情報を保存する為に、励振
パルスより後の同一の時刻に収集すること、幾何
学的な歪みを最少限に抑える為に、位相符号化パ
ルスと読取勾配の中心の間に同一の期間を保つこ
と、速度によつて誘起される移相の累積を少なく
する為に、交互のエコーで勾配の極性を反転する
こと等を含めて、上で述べた順序を変更すること
が考えられる。従つて、この発明は特許請求の範
囲によつて限定されるものであつて、ここに例と
して示した細部によつて制限されないことを承知
されたい。
【図面の簡単な説明】
第１図はNMR作像装置内に静磁界を発生する
為の磁石の中孔の斜視図で、種々の磁界勾配及び
RF磁界を発生する手段の中にサンプルを位置決
めする様子を示すと共に、種々の符号の説明にも
なつている。第２図は人間の解剖学的な一部分
と、それに対して設定された複数個の投影方向の
軸線並びに関連する投影平面図を示す略図、第２
ａ図は、一般的なNMR装置で、多重エコー
NMR造影法により、夫々相異なる投影軸線に沿
つた複数個の相異なるインターリーブ形のNMR
造影像を求める現在好ましいと考えられる一実施
例の方法の完全な１回の繰返しの間に使われる無
線周波、データ・ゲート及び磁界勾配の各信号を
示す時間座標を合せた１組のグラフ、第２ｂ図は
デカルト座標の磁界勾配を持つ装置内で現在好ま
しいと考えられる多重視覚方法、又は単一の投影
軸線に沿つて求められたNMR造影像の信号対雑
音比を高める現在好ましいと考えられる方法の完
全な１回の繰返しの間に使われる別の１組のRF、
データ・ゲート及び磁界勾配の各信号を示す時間
軸を合せた１組のグラフ、第２ｃ図は一般的な
NMR装置で、単一の投影軸線に沿つて信号対雑
音比を改善する為の現在好ましいと考えられる方
法の完全な１回の繰返しで利用されるRF、デー
タ・ゲート及び磁界勾配の各信号を示す時間軸を
合せた３番目の１組のグラフ、第３ａ図は直交す
る一対の流れ方向に沿つた情報を同時に収集する
為の現在好ましいと考えられる実施例の多重エコ
ーNMR造影順序の完全な１回の繰返しの間に使
われるRF、データ・ゲート及び磁界勾配の各信
号を示す時間軸を合せた１組のグラフ、第３ｂ図
はデカルト座標系の軸に沿つて勾配磁界を配置し
た装置で、第３図の一般的な方法の原理を利用し
て行なわれる、この発明の現在好ましいと考えら
れる一実施例の造影像を得る方法で使われるRF、
データ・ゲート及び磁界勾配の各信号を示す時間
軸を合せた１組のグラフである。

Claims (1)

  1. 【特許請求の範囲】 １ サンプルの少なくとも選ばれた部分にある流
    れる材料の複数個（Ｎ個）の核磁気共鳴
    （NMR）造影像を得る方法に於いて、 (a) サンプルを主静磁界の中に浸漬し、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させ、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するめに選ばれた第１の方向に、前記
    サンプルに印加された第１の磁界勾配中に交互
    に極性が変わる一対の流れ符号化信号パルスを
    印加し、各順序対の第１の順序での一対の流れ
    符号化パルスは、各順序対の第２の順序での一
    対の流れ符号化パルスとは極性が逆になつてお
    り、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプルの少なくとも前
    記部分によつて誘起されたNMR応答エコー信
    号から１組のデータを収集し、 (e) 前記第１及び第２の順序の内の一方で収集さ
    れたNMR応答信号データのｊ番目の組（１≦
    ｊ≦Ｎ）のデータを、前記第１及び第２の順序
    の内の他方で収集された同じｊ番目の組のデー
    タから減算して、不動の原子核からの応答デー
    タが実質的に取除かれた、複数組（Ｎ組）の差
    データの内の対応するｊ番目の組の差データを
    夫々発生し、 (f) Ｎ組全ての差データに応答して、何れも前記
    第１及び第２の方向に対して予め選ばれた関係
    を持つ平面内にある少くとも１つの造影投影像
    の夫々を発生する工程を含む方法。 ２ 前記工程(d)が、ある順序の複数個（Ｎ個）の
    収集期間の夫々ｊ番目で、第２の方向を相異なる
    ように変えることを含み、第２の方向は全ての対
    の両方の順序の同じ番号（ｊ番目）の各々の収集
    期間で同じであり、工程(e)が一方の順序のｊ番目
    の収集期間に収集されたデータを、各対の他方の
    順序の同じｊ番目の収集期間に収集されたデータ
    から減算して、ｊ番目の対の差のデータの組を発
    生する工程を含み、工程(f)がＮ個の相異なる造影
    像を発生する工程を含み、ｊ番目の像が、全ての
    ｊ番目の対の差のデータの組から発生される請求
    項１記載の方法。 ３ 工程(d)が更に、各々の順序のＮ回の収集期間
    の夫々より前に、読取勾配に位相外しパルスを発
    生し、各々の順序の最初の（Ｎ−１）個の収集期
    間の各々の後、読取勾配中に位相戻しパルスを発
    生し、各々の位相戻しパルスには、関連する位相
    外しパルスと反対の極性であるが、略同じ振幅を
    持たせる工程を含む請求項２記載の方法。 ４ 工程(d)が、更に、読取収集期間より後の位相
    戻しパルスを、次の読取収集期間より前の位相外
    しパルスと略一致するようにタイミングを定める
    工程を含む請求項３記載の方法。 ５ 工程(d)が更に、前記第１及び第２の両方の方
    向に対して略直交する第３の方向に、略不動のス
    ピンを持つ原子核から誘起される応答信号のダイ
    ナミツク・レンジを制限すると共に該応答信号の
    抑圧を増大するための投影位相外し勾配信号を印
    加する工程を含む請求項４記載の方法。 ６ 更に、読取勾配の位相外しパルスと略一致し
    た投影位相外し勾配信号を位相外しパルスとして
    発生すると共に、該投影位相外し勾配の関連する
    位相戻しパルスとは反対の極性であるが、略同じ
    振幅を持ち、読取勾配の位相外しパルスと時間的
    に略一致する位置にある位相戻しパルスを発生す
    る工程を含む請求項５記載の方法。 ７ 更に工程(d)が、前記第１及び第２の方向の両
    方と略直交する第３の方向に、略不動のスピンを
    持つ原子核から誘起される応答信号のダイナミツ
    ク・レンジを制限すると共に該応答信号の抑圧を
    増大するための投影位相外し勾配信号を印加する
    工程を含む請求項３記載の方法。 ８ 更に、読取勾配の位相外しパルスと略一致す
    る投影位相外し勾配信号を位相外しパルスとして
    発生すると共に、投影位相外し勾配の関連する位
    相戻しパルスと反対の極性であるが略同じ振幅を
    持ち且つ読取勾配の位相外しパルスと略一致する
    位相戻しパルスを発生する工程を含む請求項７記
    載の方法。 ９ 工程(d)が、読取勾配の方向に対して直交する
    方向の磁界勾配中に位相符号化パルスを発生し、
    各対の順序で、それに対して多数（Ｓ個）の相異
    なる値の内の１つを割当てる工程を含む請求項３
    記載の方法。 １０ 第１の方向が主静磁界の方向であり、相異
    なるＮ個の第２の方向が何れもそれと直交する平
    面内にある請求項５記載の方法。 １１ 第１の方向がデカルト座標系のＺ軸であ
    り、Ｎ個の第２の方向の夫々ｊ番目はＸ−Ｙ平面
    内にあつて、座標系のＹ軸に対して角度θを持
    ち、工程(d)が、座標系のＸ軸と平行な方向の磁界
    勾配Gxを主静磁界中に発生し、座標系のＹ軸と
    平行な方向の磁界勾配Gyを主静磁界中に発生し、
    工程(d)の各々の収集期間の少なくとも読取部分の
    間、夫々θ＝0°，θ＝90°に対する勾配の夫々の
    振幅をGx,nax及びGy,naxとして、勾配の振幅をGx
    ＝Gx,naxcos（θj），Gy＝Gy,naxsin（θj）に設定する
    工程を含む請求項１０記載の方法。 １２ 更に工程(d)が、各々の順序のＮ個の収集期
    間の夫々より前に、Gx及びGy勾配の夫々の中に
    位相外しパルスを発生し、各々の位相外しパルス
    の振幅を、同じ勾配の読取部分の振幅と反対の極
    性であるが、その一定の倍数になるように位取り
    する工程を含む請求項１１記載の方法。 １３ 工程(d)が、Ｚ軸方向の磁界勾配Gz中に位
    相符号化パルスを発生し、各対の順序で、この位
    相符号化パルスに多数（Ｓ個）の相異なる値の内
    の１つを割当てる工程を含む請求項１２記載の方
    法。 １４ 工程(a)が章動動作をサンプルの部分の予定
    の容積に制限する工程を含む請求項１記載の方
    法。 １５ サンプルが周期的な動作サイクルを持つも
    のであり、更に、該サイクル中の選定された点に
    対して、各順序の開始をゲートする工程を含む請
    求項１記載の方法。 １６ サンプルの少なくとも選ばれた部分にある
    流れる材料の少なくとも１つの核磁気共鳴
    （NMR）造影像を得る方法に於いて、 (a) サンプルを主静磁界の中に浸漬し、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させ、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するために選ばれた第１の方向に、前
    記サンプルに印加された第１の磁界勾配中に交
    互に極性が変わる一対の流れ符号化信号パルス
    を印加し、各順序対の第１の順序での一対の流
    れ符号化パルスは、各順序対の第２の順序での
    一対の流れ符号化パルスとは極性が逆になつて
    おり、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプル部分によつて誘
    起されたNMR応答エコー信号から１組のデー
    タを収集し、 (e) 前記第１及び第２の順序の内の一方で収集さ
    れたNMR応答信号データのｊ番目の組（１≦
    ｊ≦Ｎ）のデータを、前記第１及び第２の順序
    の内の他方で収集された同じｊ番目の組のデー
    タから減算して、不動の原子核からの応答デー
    タが実質的に取除かれた、複数組（Ｎ組）の差
    データの内の対応するｊ番目の組の差データを
    夫々発生し、 (f) Ｎ組全ての差データに応答して、何れも前記
    第１及び第２の方向に対して予め選ばれた関係
    を持つ平面内にある少なくとも１つの造影投影
    像の夫々を発生する工程を含み、 前記工程(d)が、各々の順序のＮ個の収集期間の
    １番目の期間より前に、読取勾配中に位相外しパ
    ルスを発生し、該１個の位相外しパルスの振幅
    を、第１の読取勾配と反対の極性であるが、その
    振幅の一定倍数になるように位取りし、その後、
    Ｎ個の読取勾配の夫々ｊ番目の振幅に係数（−
    １）(j+1)を乗じて、この順序中の全ての読取勾配
    の極性を交互に変える工程を含む方法。 １７ 工程(d)が、Ｎ個の読取勾配の内の（ｊ−
    １）番目より後並びに夫々ｊ番目より前に、読取
    勾配の方向に対して略直交する選ばれた方向に磁
    界勾配中に位相符号化パルスを発生し、最初の
    （Ｎ−１）個の読取勾配の内の夫々ｊ番目より後
    且つ夫々（ｊ＋１）番目より前に、前記選ばれた
    方向の磁界勾配中に極性が反転した位相符号化パ
    ルスを発生し、各対の順序で、位相符号化パルス
    の振幅に、多数（Ｓ個）の相異なる値の内の１つ
    を割当てる工程を含む請求項１６記載の方法。 １８ 更に工程(d)が、前記第１及び第２の方向の
    両方と略直交する第３の方向に、略不動のスピン
    を持つ原子核から誘起される応答信号の抑圧を増
    大するための投影位相外し勾配信号を印加する工
    程を含む請求項１７記載の方法。 １９ 工程(a)が章動動作をサンプルの部分の予定
    の容積に制限する工程を含む請求項１６記載の方
    法。 ２０ サンプルが周期的な動作サイクルを持つも
    のであり、更に、該サイクル中の選定された点に
    対して、各順序の開始をゲートする工程を含む請
    求項１６記載の方法。 ２１ サンプルの少なくとも選ばれた部分にある
    流れる材料の少なくとも１つの核磁気共鳴
    （NMR）造影像を得る方法に於いて、 (a) サンプルを主静磁界の中に浸漬し、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させ、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するために選ばれた第１の方向に、前
    記サンプルに印加された第１の磁界勾配中に交
    互に極性が変わる一対の流れ符号化信号パルス
    を印加し、各順序対の第１の順序での一対の流
    れ符号化パルスは、各順序対の第２の順序での
    一対の流れ符号化パルスとは極性が逆になつて
    おり、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプル部分によつて誘
    起されたNMR応答エコー信号から１組のデー
    タを収集し、１≦ｊ≦Ｎとして、前記第２の方
    向はＮ個の収集期間のｊ番目の期間において、
    夫々他の期間におけるものとは相異なるように
    変えられ、 該工程(d)は、各々の順序で、１より大きく且
    つＮより小さいＭの少なくとも１つの値に対
    し、少なくとも１つの相異なる予め選ばれた第
    ２の方向の各々に於いて、該順序内の同じ方向
    に於けるＭ個のエコーの相異なる１つから夫々
    求められた複数個（Ｍ個）の応答データの組を
    収集し、該データの組は各順序に於けるＭ個の
    エコーの同じ組から収集され、各々順序に於け
    るＭ個の応答データの組全部のデータの値の和
    をとつて前記予め選ばれた第２の方向及び該順
    序に対する１個の和データの組を求める工程を
    含み、 (e) 前記第１及び第２の順序の内の一方の和デー
    タの組を、前記対の他方の順序から得られた和
    データの組から減算して、該対に対する差デー
    タの組を発生し、 (f) 各組の差データから、対応する相異なる第２
    の方向ならびに第１の方向に対して予め選ばれ
    た関係を持つ平面内にある相異なる１つの造影
    像を夫々発生する工程を含む方法。 ２２ 工程(a)が章動動作をサンプルの部分の予定
    の容積に制限する工程を含む請求項２１記載の方
    法。 ２３ サンプルが周期的な動作サイクルを持つも
    のであり、更に、該サイクル中の選定された点に
    対して、各順序の開始をゲートする工程を含む請
    求項２１記載の方法。 ２４ サンプルの少なくとも選ばれた部分にある
    流れる材料の少なくとも１つの核磁気共鳴
    （NMR）造影像を得る方法に於いて、 (a) サンプルを主静磁界の中に浸漬し、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させ、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するために選ばれた第１の方向に、前
    記サンプルに印加された第１の磁界勾配中に交
    互に極性が変わる一対の流れ符号化信号パルス
    を印加し、各順序対の第１の順序での一対の流
    れ符号化パルスは、各順序対の第２の順序での
    一対の流れ符号化パルスとは極性が逆になつて
    おり、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｍ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプルの少なくとも前
    記部分によつて誘起されたNMR応答エコー信
    号から１組のデータを収集し、 該工程(d)は、１個の共通の読取用の第２の方
    向に対し、各々の順序で応答データの組のＭ個
    全部を求め、各順序に於けるＭ個の応答データ
    の組全部に対するデータの値の和をとつて、該
    共通の第２の方向及び該順序に対する１個の和
    データの組を求める工程を含み、 (e) 前記第１及び第２の順序の内の一方の和デー
    タの組を、該対の他方の順序の和データの組か
    ら減算して、該対に対する差データの組を発生
    し、 (f) 前記の差データの組に応答して、前記第１及
    び第２の方向に対して予め選ばれた関係を持つ
    平面内にある少なくとも１つの造影像を夫々発
    生する工程を含む方法。 ２５ 工程(a)が章動動作をサンプルの部分の予定
    の容積に制限する工程を含む請求項２４記載の方
    法。 ２６ サンプルが周期的な動作サイクルを持つも
    のであり、更に、該サイクル中の選定された点に
    対して、各順序の開始をゲートする工程を含む請
    求項２４記載の方法。 ２７ サンプルの少なくとも選ばれた部分にある
    流れ材料の少なくとも１つの核磁気共鳴
    （NMR）造影像を得る方法に於いて、 (a) サンプルを主静磁界の中に浸漬し、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させ、 （c1） 各々の順序に於ける複数個（Ｎ個）の応
    答データ収集期間の最初の期間より前に、一対
    の流れ符号化方向の第１の方向で、サンプルに
    印加される第１の磁界勾配中に極性が交互に変
    わる一対の流れ符号化信号パルスを印加し、 （c2） 各々の順序の第１の応答期間の後且つ第
    ２の応答データ収集期間より前に、前記第１の
    流れ符号化方向と同じ方向の第１の磁界勾配中
    に、工程（c1）の一対のパルスの反転である極
    性が交互に変わる別の一対の流れ符号化信号パ
    ルスを印加し、 （c3） 各順序に於ける第２の応答期間より前
    に、前記第１の磁界勾配と直交する第２の磁界
    勾配中に他方の流れ符号化方向に極性が交互に
    変わる一対の流れ符号化信号パルスを印加し、
    各順序対の第１の順序における一対の流れ符号
    化パルスは各順序対の第２の順序における一対
    の流れ符号化パルスとは極性が逆になつてお
    り、各流れ符号化方向は、動く原子核のスピン
    から得られるNMR応答エコー信号が、実質的
    に不動の原子核のスピンから得られるNMR応
    答エコー信号と異なるように、対応するNMR
    造影像の１つの軸を設定するために選ばれてお
    り、 (d) 前記流れ符号化方向とは独立の１つの方向に
    前記サンプルに印加された読取磁界勾配に応答
    して、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプルの少なくとも前
    記部分によつて誘起されたNMR応答エコー信
    号から１組のデータを収集し、 (e) 前記第１の順序で収集されたNMR応答エコ
    ー信号による第１及び第２組のデータを、前記
    第２の順序で収集された第１及び第２の組のデ
    ータから夫々に減算して、不動の原子核からの
    応答データが実質的に取除かれた、一対の第１
    組IA及び第２組IBの差データを夫々発生し、 (f) 前記第１の方向の差のデータIAの組及び前記
    第２の方向の差のデータIBの組を処理して、 It＝｛（IA）2＋（IB）2｝1/2 によつて表わされる値を持つ最終的な差のデー
    タの組Itを求めて、流れ全体のNMR管造影図
    を表示する工程を含む方法。 ２８ 工程（c2）及び（c3）が略同時に行なわれ
    る請求項２７記載の方法。 ２９ 更に工程(d)が各順序に於ける２つの応答収
    集期間の各々より前に、読取勾配中に位相外しパ
    ルスを発生し、各順序に於ける最初の応答収集期
    間より後に該読取勾配中に位相戻しパルスを発生
    し、該位相戻しパルスには、前記順序の２つの位
    相外しパルスの略相等しい振幅と略同じ振幅を持
    たせる工程を含む請求項２７記載の方法。 ３０ 更に工程(d)が、前記第１の読取勾配の位相
    外しパルスと略同時に、読取勾配の方向に対して
    略直交する方向に、略不動のスピンを持つ原子核
    から誘起される応答信号のダイナミツク・レンジ
    を制限することにより抑圧を増大するための投影
    位相外し勾配信号パルスを印加する工程を含む請
    求項２９記載の方法。 ３１ 第１の方向が主静磁界の方向であり、２つ
    の異なる第２の方向が何れもそれに対して直交す
    る平面内にある請求項３０記載の方法。 ３２ 工程(d)が、読取勾配の方向に対して直交す
    る方向の磁界勾配中に位相符号化パルスを発生
    し、各対の順序で、それに対して多数（Ｓ個）の
    相異なる値の内の１つを割当てる工程を含む請求
    項２７記載の方法。 ３３ 工程(a)が章動動作サンプルの部分の予定の
    容積に制限する工程を含む請求項２７記載の方
    法。 ３４ サンプルが周期的な動作サイクルを持つも
    のであり、更に、該サイクル中の選定された点に
    対して、各順序の開始をゲートする工程を含む請
    求項２７記載の方法。 ３５ サンプルの少なくとも選ばれた部分にある
    流れる材料の核磁気共鳴（NMR）造影像を得る
    方法に於いて、 (a) サンプルを主静磁界の中に浸漬し、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させ、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するために選ばれた第１の方向に、前
    記サンプルに印加された第１の磁界勾配中に交
    互に極性が変わる一対の流れ符号化信号パルス
    を印加し、各順序の第１の順序での一対の流れ
    符号化パルスは、各順序対の第２の順序での一
    対の流れ符号化パルスとは極性が逆になつてお
    り、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の２つの応
    答データ収集期間の夫々相異なる１つの期間内
    に、前記サンプルの少なくとも前記部分によつ
    て誘起されたNMR応答エコー信号から１組の
    データを収集し、 (e) 前記第１及び第２の順序の内の一方で収集さ
    れたNMR応答信号データのｊ番目の組（１≦
    ｊ≦２）のデータを、前記第１及び第２の順序
    の他方で収集された同じｊ番目の組のデータか
    ら減算して、不動の原子核からの応答データが
    実質的に取除かれた、２組の差データ内の対応
    するｊ番目の組の差データを夫々発生し、 (f) 前記全ての組の差データに応答して、前記第
    １及び第２の方向に対して予め選ばれた関係を
    持つ平面内にある造影投影像を発生する工程を
    含み、 前記工程(c)が（c1）各順序に於ける最初の応答
    データ収集期間より前に、流れ符号化信号を印加
    せず、（c2）各順序に於ける最初の応答期間より
    後且つ第２の応答データ収集期間より前に、第１
    の流れ符号化勾配方向に、第１の磁界勾配中に極
    性が交互に変わる一対の流れ符号化信号パルスを
    印加し、（c3）各対の第２の順序で前記工程（c1）
    及び（c2）を繰返す工程を含み、 前記工程(e)が、（e1）第１の順序の２番目の応
    答収集期間から得られたデータの組S12を、第１
    の順序の最初の応答収集期間から得られたデータ
    の組S11から減算してデータの組S1を形成し、
    （e2）第２の順序中の２番目の応答収集期間から
    えられたデータの組S22を第２の順序中の最初の
    応答収集期間から得られたデータの組S21から減
    算してデータの組S2を形成し、（e3）データの組
    S2をデータの組S1から減算して差データの組を求
    める工程を含み、 前記工程(f)が動きの人為効果の抑圧を強めた
    NMR造影図の像を表示する工程を含む方法。 ３６ 工程(a)が章動動作をサンプルの部分の予定
    の容積に制限する工程を含む請求項３５記載の方
    法。 ３７ サンプルが周期的な動作サイクルを持つも
    のであり、更に、該サイクル中の選定された点に
    対して、各順序の開始をゲートする工程を含む請
    求項３５記載の方法。 ３８ サンプルの少なくとも選ばれた部分にある
    流れる材料の複数の核磁気共鳴（NMR）造影像
    を得る方法に於いて、 (a) サンプルを主静磁界の中に浸漬し、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させ、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核スピンから得られるNMR応答エコー
    信号と異なるように、NMR造影像の１つの軸
    を設定するために選ばれた第１の方向に、前記
    サンプルに印加された第１の磁界勾配中に交互
    に極性が変わる一対の流れ符号化信号パルスを
    複数個（Ｎ個）の応答収集期間の各々より前に
    印加し、各順序対の第１の順序での一対の流れ
    符号化パルスは、各順序対の第２の順序での一
    対の流れ符号化パルスとは極性が逆になつてお
    り、流れ符号化信号の全てのパルス対は第１の
    パルスが同じ極性であり、前記第１の方向は流
    れ符号化方向であり、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプルの少なくとも前
    記部分によつて誘起されたNMR応答エコー信
    号から１組のデータを収集し、 (e) 前記第１及び第２の順序の内の一方で収集さ
    れたNMR応答信号データのｊ番目の組（１≦
    ｊ≦Ｎ）のデータを、前記第１及び第２の順序
    の内の他方で収集された同じｊ番目の組のデー
    タから減算して、不動の原子核からの応答デー
    タが実質的に取除かれた、複数組（Ｎ組）の差
    データの内の対応するｊ番目の組の差データを
    夫々発生し、 (f) Ｎ組全ての差データに応答して、何れも前記
    第１及び第２の方向に対して予め選ばれた関係
    を持つ平面内にある複数の造影投影像の夫々を
    発生し、この結果得られたＮ個の差データの組
    の各々から、他のどの１つの差データの組から
    表示される速度情報とも異なる、サンプルの部
    分の流体の流速に関する情報を持つNMR造影
    図の像を表示する工程を含む方法。 ３９ 工程(a)が章動動作をサンプルの部分の予定
    の容積に制限する工程を含む請求項３８記載の方
    法。 ４０ サンプルが周期的な動作サイクルを持つも
    のであり、更に、該サイクル中の選定された点に
    対して、各順序の開始をゲートする工程を含む請
    求項３８記載の方法。 ４１ サンプルの少なくとも選ばれた部分にある
    流れる材料の複数個（Ｎ個）の核磁気共鳴
    （NMR）造影像を得る装置に於いて、 (a) サンプルを主静磁界の中に浸漬する手段と、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させる手段と、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するために選ばれた第１の方向に、前
    記サンプルに印加された第１の磁界勾配中に交
    互に極性が変わる一対の流れ符号化信号パルス
    を印加する手段であつて、各順序対の第１の順
    序での一対の流れ符号化パルスが、各順序対の
    第２の順序での一対の流れ符号化パルスとは極
    性が逆になつている手段と、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプルの少なくとも前
    記部分によつて誘起されたNMR応答エコー信
    号から１組のデータを収集する手段と、 (e) 前記第１及び第２の順序の内の一方で収集さ
    れたNMR応答信号データのｊ番目の組（１≦
    ｊ≦Ｎ）のデータを、前記第１及び第２の順序
    の内の他方で収集された同じｊ番目の組のデー
    タから減算して、不動の原子核からの応答デー
    タが実質的に取除かれた、複数組（Ｎ組）の差
    データの内の対応するｊ番目の組の差データを
    夫々発生する手段と、 (f) Ｎ組全ての差データに応答して、何れも前記
    第１及び第２の方向に対して予め選ばれた関係
    を持つ平面内にある少なくとも１つの造影投影
    像の夫々を発生する手段とを含む装置。 ４２ 前記手段(d)が、ある順序の複数個（Ｎ個）
    の収集期間の夫々ｊ番目で、第２の方向を相異な
    るように変えることを含み、第２の方向は全ての
    対の両方の順序の同じ番号（ｊ番目）の各々の収
    集期間で同じであり、手段(e)が一方の順序のｊ番
    目の収集期間に収集されたデータを、各対の他方
    の順序の同じｊ番目の収集期間に収集されたデー
    タから減算して、ｊ番目の対の差のデータの組を
    発生する手段を含み、手段(f)がＮ個の相異なる造
    影像を発生する手段を含み、ｊ番目の像が、全て
    のｊ番目の対の差のデータの組から発生される請
    求項４１記載の装置。 ４３ 手段(d)が更に、各々の順序のＮ回の収集期
    間の夫々より前に、読取勾配に位相外しパルスを
    発生し、各々の順序の最初の（Ｎ−１）個の収集
    期間の各々の後、読取勾配中に位相戻しパルスを
    発生し、各々の位相戻しパルスには、関連する位
    相外しパルスと反対の極性であるが、略同じ振幅
    を持たせる手段を含む請求項４２記載の装置。 ４４ 手段(d)が、更に、読取収集期間より後の位
    相戻しパルスを、次の読取収集期間より前の位相
    外しパルスと略一致するようにタイミングを定め
    る手段を含む請求項４３記載の装置。 ４５ 手段(d)が更に、前記第１及び第２の両方の
    方向に対して略直交する第３の方向に、略不動の
    スピンを持つ原子核から誘起される応答信号のダ
    イナミツク・レンジを制限すると共に該応答信号
    の抑圧を増大するための投影位相外し勾配信号を
    印加する手段を含む請求項４４記載の装置。 ４６ 更に、読取勾配の位相外しパルスと略一致
    した投影位相外し勾配信号を位相外しパルスとし
    て発生すると共に、該投影位相外し勾配の関連す
    る位相戻しパルスとは反対の極性であるが、略同
    じ振幅を持ち、読取勾配の位相外しパルスと時間
    的に略一致する位置にある位相戻しパルスを発生
    する手段を含む請求項４５記載の装置。 ４７ 更に手段(d)が、前記第１及び第２の方向の
    両方と略直交する第３の方向に、略不動のスピン
    を持つ原子核から誘起される応答信号のダイナミ
    ツク・レンジを制限すると共に該応答信号の抑圧
    を増大するめの投影位相外し勾配信号を印加する
    手段を含む請求項４３記載の装置。 ４８ 更に、読取勾配の位相外しパルスと略一致
    する投影位相外し勾配信号を位相外しパルスとし
    て発生すると共に、投影位相外し勾配の関連する
    位相戻しパルスと反対の極性であるが略同じ振幅
    を持ち且つ読取勾配の位相外しパルスと略一致す
    る移動戻しパルスを発生する手段を含む請求項４
    ７記載の装置。 ４９ 手段(d)が、読取勾配の方向に対して直交す
    る方向の磁界勾配中に位相符号化パルスを発生
    し、各対の順序で、それに対して多数（Ｓ個）の
    相異なる値の内の１つを割当てる手段を含む請求
    項４３記載の装置。 ５０ 第１の方向が主静磁界の方向であり、相異
    なるＮ個の第２の方向が何れもそれと直交する平
    面内にある請求項４５記載の装置。 ５１ 第１の方向がデカルト座標系のＺ軸であ
    り、Ｎ個の第２の方向の夫々ｊ番目はＸ−Ｙ平面
    内にあつて、座標系のＹ軸に対して角度θを持
    ち、手段(d)が、座標系のＸ軸と平行な方向の磁界
    勾配Gxを主静磁界中に発生し、座標系のＹ軸と
    平行な方向の磁界勾配Gyを主静磁界中に発生し、
    手段(d)の各々の収集期間の少なくとも読取部分の
    間、夫々θ＝0°，θ＝90°に対する勾配の夫々の
    振幅をGx,nax及びGy,naxとして、勾配の振幅をGx
    ＝Gx,naxcos（θj），Gy＝Gy,naxsin（θj）に設定する
    手段を含む請求項５０記載の装置。 ５２ 更に手段(d)が、各々の順序のＮ個の収集期
    間の夫々より前に、Gx及びGy勾配の夫々の中に
    位相外しパルスを発生し、各々の位相外しパルス
    の振幅を、同じ勾配の読取部分の振幅と反対の極
    性であるが、その一定の倍数になるように位取り
    する手段を含む請求項５１記載の装置。 ５３ 手段(d)が、Ｚ軸方向の磁界勾配Gz中に位
    相符号化パルスを発生し、各対の順序で、この位
    相符号化パルスに多数（Ｓ個）の相異なる値の内
    の１つを割当てる手段を含む請求項５２記載の装
    置。 ５４ 手段(a)が章動動作をサンプルの部分の予定
    の容積に制限する手段を含む請求項４１記載の装
    置。 ５５ サンプルが心臓サイクルを持つ生体の一部
    分であり、更に、該生体への心臓サイクル中の選
    定された点に対して、各順序の開始をゲートする
    手段を含む請求項４１記載の装置。 ５６ サンプルの少なくとも選ばれた部分にある
    流れる材料の少なくとも１つの核磁気共鳴
    （NMR）造影像を得る装置において、 (a) サンプルを主静磁界の中に浸漬する手段と、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させる手段と、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するめに選ばれた第１の方向に、前記
    サンプルに印加された第１の磁界勾配中に交互
    に極性が変わる一対の流れ符号化信号パルスを
    印加する手段であつて、各順序対の第１の順序
    での一対の流れ符号化パルスが、各順序対の第
    ２の順序での一対の流れ符号化パルスとは極性
    が逆になつている手段と、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプル部分によつて誘
    起されたNMR応答エコー信号から１組のデー
    タを収集する手段と、 (e) 前記第１及び第２の順序の内の一方で収集さ
    れたNMR応答信号データのｊ番目の組（１≦
    ｊ≦Ｎ）のデータを、前記第１及び第２の順序
    の内の他方で収集された同じｊ番目の組のデー
    タから減算して、不動の原子核からの応答デー
    タが実質的に取除かれた、複数組（Ｎ組）の差
    データの内の対応するｊ番目の組の差データを
    夫々発生する手段と、 (f) Ｎ組全ての差データに応答して、何れも前記
    第１及び第２の方向に対して予め選ばれた関係
    を持つ平面内にある少なくとも１つの造影投影
    像の夫々を発生する手段とを含み、 前記手段(d)が、各々の順序のＮ個の収集期間
    の１番目だけより前に、読取勾配中に位相外し
    パルスを発生し、該１個の位相外しパルスの振
    幅を、第１の読取勾配と反対の極性であるが、
    その振幅の一定倍数になるように位取りし、そ
    の後、Ｎ個の読取勾配の夫々ｊ番目の振幅に係
    数（−１）(j+1)を乗じて、この順序中の全ての
    読取勾配の極性を交互に変える手段を含む装
    置。 ５７ 手段(d)が、Ｎ個の読取勾配の内の（ｊ−
    １）番目より後並びに夫々ｊ番目より前に、読取
    勾配の方向に対して略直交する選ばれた方向の磁
    界勾配中に位相符号化パルスを発生し、最初の
    （Ｎ−１）個の読取勾配の内の夫々ｊ番目より後
    且つ夫々（ｊ＋１）番目より前に、前記選ばれた
    方向の磁界勾配中に極性が反転した位相符号化パ
    ルスを発生し、各対の順序で、位相符号化パルス
    の振幅に、多数（Ｓ個）の相異なる値の内の１つ
    を割当てる手段を含む請求項５６記載の装置。 ５８ 更に手段(d)が、前記第１及び第２の方向の
    両方と略直交する第３の方向に、略不動のスピン
    を持つ原子核から誘起される応答信号の抑圧を増
    大するための投影位相外し勾配信号を印加する手
    段を含む請求項５６記載の装置。 ５９ 手段(a)が章動動作をサンプルの部分の予定
    の容積に制限する手段を含む請求項５６記載の装
    置。 ６０ サンプルが心臓サイクルを持つ生体の一部
    分であり、更に、該生体の心臓サイクル中の選定
    された点に対して、各順序の開始をゲートする手
    段を含む請求項５６記載の装置。 ６１ サンプルの少なくとも選ばれた部分にある
    流れる材料の少なくとも１つの核磁気共鳴
    （NMR）造影像を得る装置に於いて、 (a) サンプルを主静磁界の中に浸漬する手段と、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させる手段と、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するために選ばれた第１の方向に、前
    記サンプルに印加された第１の磁界勾配中に交
    互に極性が変わる一対の流れ符号化信号パルス
    を印加する手段であつて、各順序対の第１の順
    序での一対の流れ符号化パルスが、各順序対の
    第２の順序での一対の流れ符号化パルスとは極
    性が逆になつている手段と、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプル部分によつて誘
    起されたNMR応答エコー信号から１組のデー
    タを収集する手段であつて、１≦ｊ≦Ｎとし
    て、前記第２の方向はＮ個の収集期間のｊ番目
    の期間において夫々他の期間におけるものとは
    異なるように変えられ、 更に、各々の順序で、１より大きく且つＮよ
    り小さいＭの少なくとも１つの値に対し、少な
    くとも１つの相異なる予め選ばれた第２の方向
    の各々に於いて、該順序内の同じ方向に於ける
    Ｍ個のエコーの相異なる１つから夫々求められ
    た複数個（Ｍ個）の応答データの組を収集し、
    該データの組は各順序に於けるＭ個のエコーの
    同じ組から収集され、各々順序に於けるＭ個の
    応答データの組全部のデータの値の和をとつて
    前記予め選ばれた第２の方向及び該順序に対す
    る１個の和データの組を求める手段を含む手段
    と、 (e) 前記第１及び第２の順序の内の一方の和デー
    タの組を、前記対の他方の順序から得られた和
    データの組から減算して、該対に対する差デー
    タの組を発生する手段と、 (f) 各組の差データから、対応する相異なる第２
    の方向ならびに第１の方向に対して予め選ばれ
    た関係を持つ平面内にある相異なる１つの造影
    像を夫々発生する手段とを含む装置。 ６２ 手段(a)が章動動作をサンプルの部分の予定
    の容積に制限する手段を含む請求項６１記載の装
    置。 ６３ サンプルが心臓サイクルを持つ生体の一部
    分であり、更に、該生体の心臓サイクル中の選定
    された点に対して、各順序の開始をゲートする手
    段を含む請求項６１記載の装置。 ６４ サンプルの少なくとも選ばれた部分にある
    流れる材料の少なくとも１つの核磁気共鳴
    （NMR）造影像を得る装置に於いて、 (a) サンプルを主静磁界の中に浸漬する手段と、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させる手段と、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するために選ばれた第１の方向に、前
    記サンプルに印加された第１の磁界勾配中に交
    互に極性が変わる一対の流れ符号化信号パルス
    を印加する手段であつて、各順序対の第１の順
    序での一対の流れ符号化パルスが、各順序対の
    第２の順序での一対の流れ符号化パルスとは極
    性が逆になつている手段と、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｍ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプル部分によつて誘
    起されたNMR応答エコー信号から１組のデー
    タを収集する手段であつて、１個の共通の読取
    用の第２の方向に対し、各々の順序で応答デー
    タの組のＭ個全部を求め、各順序に於けるＭ個
    の応答データの組全部に対するデータの値の和
    をとつて、該共通の第２の方向及び該順序に対
    する１個の和データの組を求める手段を含む手
    段と、 (e) 前記第１及び第２の順序の内の一方の和デー
    タの組を、該対の他方の順序の和データの組か
    ら減算して、該対に対する差データの組を発生
    する手段を含む手段と、 (f) 各組の差データの組に応答して、前記第１及
    び第２の方向に対して予め選ばれた関係を持つ
    平面内にある少なくとも１つの造影像を夫々発
    生する手段とを含む装置。 ６５ 手段(a)が章動動作をサンプルの部分の予定
    の容積に制限する手段を含む請求項６４記載の装
    置。 ６６ サンプルが心臓サイクルを持つ生体の一部
    分であり、更に、該生体の心臓サイクル中の選定
    された点に対して、各順序の開始をゲートする手
    段を含む請求項６４記載の装置。 ６７ サンプルの少なくとも選ばれた部分にある
    流れる材料の少なくとも１つの核磁気共鳴
    （NMR）造影像を得る装置に於いて、 (a) サンプルを主静磁界の中に浸漬する手段と、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させる手段と、 （c1） 各々の順序に於ける最初の応答データ収
    集期間より前に、一対の流れ符号化方向の第１
    の方向で、サンプルに印加される第１の磁界勾
    配中に極性が交互に変わる一対の流れ符号化信
    号パルスを印加する手段と、 （c2） 各々の順序の第１の応答期間の後且つ第
    ２の応答データ収集期間より前に、前記第１の
    流れ符号化方向と同じ方向の第１の磁界勾配中
    に、手段（c1）の一対のパルスの反転である極
    性が交互に変わる別の一対の流れ符号化信号パ
    ルスを印加する手段と、 （c3） 各順序に於ける第２の応答期間より前
    に、前記第１の磁界勾配中に他方の流れ符号化
    方向に極性が交互に変わる一対の流れ符号化信
    号パルスを印加する手段であつて、各対の第１
    の順序における一対の流れ符号化パルスは各対
    の第２の順序における一対の流れ符号化パルス
    とは極性が逆になつており、各流れ符号化方向
    は、動く原子核のスピンから得られるNMR応
    答エコー信号が、実質的に不動の原子核のスピ
    ンから得られるNMR応答エコー信号と異なる
    ように、対応するNMR造影像の１つの軸を設
    定するために選ばれており、 (d) 前記流れ符号化方向とは独立の１つの方向に
    前記サンプルに印加された読取磁界勾配に応答
    して、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプルの少なくとも前
    記部分によつて誘起されたNMR応答エコー信
    号から１組のデータを収集する手段と、 (e) 前記第１の順序で収集されたNMR応答エコ
    ー信号による第１及び第２の組のデータを、前
    記第２の順序で収集された第１及び第２の組の
    データから減算して、不動の原子核からの応答
    データが実質的に取除かれた、一対の第１組IA
    及び第２組IBの差データを夫々発生する手段
    と、 (f) 前記第１の方向の差のデータIAの組及び第２
    の方向の差のデータIBの組を処理して、 It＝｛（IA）2＋（IB）2｝1/2 によつて表わされる値を持つ最終的な差のデー
    タの組Itを求めて、流れ全体のNMR造影図を
    表示する手段とを含む装置。 ６８ 手段（c2）及び（c3）が略同時に行なわれ
    る請求項６７記載の装置。 ６９ 更に手段(d)が各順序に於ける２つの応答収
    集期間の各々より前に、読取勾配中に位相外しパ
    ルスを発生し、各順序に於ける最初の応答収集期
    間より後に該読取勾配中に位相戻しパルスを発生
    し、該位相戻しパルスには、前記順序の２つの位
    相外しパルスの略相等しい振幅と略同じ振幅を持
    たせる手段を含む請求項６７記載の装置。 ７０ 更に手段(d)が、前記第１の読取勾配の位相
    外しパルスと略同時に、読取勾配方向に対して略
    直交する方向に、略不動のスピンを持つ原子核か
    ら誘起される応答信号のダイナミツク・レンジを
    制限することにより抑圧を高めるための投影位相
    外し勾配信号パルスを印加する手段を含む請求項
    ６９記載の装置。 ７１ 第１の方向が主静磁界の方向であり、２つ
    の相異なる第２の方向が何れもそれに対して直交
    する平面内にある請求項７０記載の装置。 ７２ 手段(d)が、読取勾配の方向に対して直交す
    る方向の磁界勾配中に位相符号化パルスを発生
    し、各対の順序で、それに対して多数（Ｓ個）の
    相異なる値の内の１つを割当てる手段を含む請求
    項６７記載の装置。 ７３ 手段(a)が章動動作をサンプルの部分の予定
    の容積に制限する手段を含む請求項６７記載の装
    置。 ７４ サンプルが心臓サイクルを持つ生体の一部
    分であり、更に、該生体の心臓サイクル中の選定
    された点に対して、各順序の開始をゲートする手
    段を含む請求項６７記載の装置。 ７５ サンプルの少なくとも選ばれた部分にある
    流れる材料の核磁気共鳴（NMR）造影像を得る
    装置に於いて、 (a) サンプルを主静磁界の中に浸漬する手段と、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させる手段と、 (c) その結果として、動く原子核のスピンから得
    れるNMR応答エコー信号が、実質的に不動の
    原子核のスピンから得れるNMR応答エコー信
    号と異なるように、NMR造影像の１つの軸を
    設定するために選ばれた第１の方向に、前記サ
    ンプルに印加された第１の磁界勾配中に交互に
    極性が変わる一対の流れ符号化信号パルスを印
    加する手段であつて、各順序対の第１の順序で
    の一対の流れ符号化パルスが、各順序対の第２
    の順序での一対の流れ符号化パルスとは極性が
    逆になつている手段と、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の２つの応
    答データ収集期間の夫々相異なる１つの期間内
    に、前記サンプルの少なくとも前記部分によつ
    て誘起されたNMR応答エコー信号から１組の
    データを収集する手段と、 (e) 前記第１及び第２の順序の内の一方で収集さ
    れたNMR応答信号データのｊ番目の組（１≦
    ｊ≦２）のデータを、前記第１及び第２の順序
    の内の他方で収集された同じｊ番目の組のデー
    タから減算して、不動の原子核からの応答デー
    タが実質的に取除かれた、２組の差データの内
    の対応するｊ番目の組の差データを夫々発生す
    る手段と、 (f) 全ての組の差データに応答して、何れも前記
    第１及び第２の方向に対して予め選ばれた関係
    を持つ平面内にある造影投影像を発生する手段
    とを含み、 前記手段(c)が（c1）各順序に於ける最初の応答
    データ収集期間より前に、流れ符号化信号を印加
    させない手段、（c2）各順序に於ける最初の応答
    期間より後且つ第２の応答データ収集期間より前
    に、第１の流れ符号化勾配方向に、第１の磁界勾
    配中に極性が交互に変わる一対の流れ符号化信号
    パルスを印加する手段、（c3）各対の第２の順序
    で前記手段（c1）及び（c2）を繰返し作動する手
    段を含み、 前記手段(e)が、（e1）第１の順序の２番目の応
    答収集期間から得られたデータの組S12を、第１
    の順序の最初の応答収集期間から得られたデータ
    の組S11から減算してデータの組S1を形成する手
    段、（e2）第２の順序中の２番目の応答収集期間
    から得られたデータの組S22を第２の順序中の最
    初の応答収集期間から得られたデータの組S21か
    ら減算してデータの組S2を形成する手段、（e3）
    データの組S2をデータの組S1から減算して差デー
    タの組を求める手段を含み、 前記手段(f)が動きの人為効果の抑圧を強めた
    NMR造影図の像を表示する手段を含む装置。 ７６ 手段(a)が章動動作をサンプルの部分の予定
    の容積に制限する手段を含む請求項７５記載の装
    置。 ７７ サンプルが心臓サイクルを持つ生体の一部
    分であり、更に、該生体の心臓サイクル中の選定
    された点に対して、各順序の開始をゲートする手
    段を含む請求項７５記載の装置。 ７８ サンプルの少なくとも選ばれた部分にある
    流れる材料の複数の核磁気共鳴（NMR）造影像
    を得る装置に於いて、 (a) サンプルを主静磁界の中に浸漬する手段と、 (b) 前記サンプルの選ばれた部分の多数（Ｓ個）
    の領域の各々に対する逐次的な一対の作像順序
    の内の第１の順序及び第２の順序の各々の初期
    部分の間、選ばれた種目の全ての原子核のスピ
    ンを章動させる手段と、 (c) その結果として、動く原子核のスピンから得
    られるNMR応答エコー信号が、実質的に不動
    の原子核のスピンから得られるNMR応答エコ
    ー信号と異なるように、NMR造影像の１つの
    軸を設定するために選ばれた第１の方向に、前
    記サンプルに印加された第１の磁界勾配中に交
    互に極性が変わる一対の流れ符号化信号パルス
    を、複数個（Ｎ個）の応答収集期間の各々より
    前に印加する手段であつて、各順序対の第１の
    順序での一対の流れ符号化パルスが、各順序対
    の第２の順序での一対の流れ符号化パルスとは
    極性が逆になつており、流れ符号化信号の全て
    のパルス対は第１のパルスが同じ極性であり、
    前記第１の方向は流れ符号化方向である手段
    と、 (d) 前記第１の方向と略直交する第２の方向に前
    記サンプルに印加された読取磁界勾配に応答し
    て、前記第１及び第２の順序の夫々の複数個
    （Ｎ個）の応答データ収集期間の夫々相異なる
    １つの期間内に、前記サンプルの少なくとも前
    記部分によつて誘起されたNMR応答エコー信
    号から１組のデータを収集する手段と、 (e) 前記第１及び第２の順序の内の一方で収集さ
    れたNMR応答信号データのｊ番目の組（１≦
    ｊ≦Ｎ）のデータを、前記第１及び第２の順序
    の内の他方で収集された同じｊ番目の組のデー
    タから減算して、不動の原子核からの応答デー
    タが実質的に取除かれた、複数組（Ｎ組）の差
    データの内の対応するｊ番目の組の差データを
    夫々発生する手段と、 (f) Ｎ組全ての差データに応答して、何れも前記
    第１及び第２の方向に対して予め選ばれた関係
    を持つ平面内にある複数の造影投影像の夫々を
    発生し、この結果得られたＮ個の差データの組
    の各々から、他のどの１つの差データの組から
    表示される速度情報とも異なる、サンプルの部
    分の流体の流速に関する情報を持つNMR造影
    図の像を表示する手段とを含む装置。 ７９ 手段(a)が章動動作をサンプルの部分の予定
    の容積に制限する手段を含む請求項７８記載の装
    置。 ８０ サンプルが心臓サイクルを持つ生体の一部
    分であり、更に、該生体の心臓サイクル中の選定
    された点に対して、各順序の開始をゲートする手
    段を含む請求項７８記載の装置。
JP63027838A 1987-02-11 1988-02-10 核磁気共鳴造影像を得る方法と装置 Granted JPS63226343A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US07/013,592 US4796635A (en) 1987-02-11 1987-02-11 Multiple-echo, multiple-view NMR angiography
US13592 1987-02-11

Publications (2)

Publication Number Publication Date
JPS63226343A JPS63226343A (ja) 1988-09-21
JPH0424051B2 true JPH0424051B2 (ja) 1992-04-24

Family

ID=21760727

Family Applications (1)

Application Number Title Priority Date Filing Date
JP63027838A Granted JPS63226343A (ja) 1987-02-11 1988-02-10 核磁気共鳴造影像を得る方法と装置

Country Status (5)

Country Link
US (1) US4796635A (ja)
EP (1) EP0278271A3 (ja)
JP (1) JPS63226343A (ja)
FI (1) FI880585A (ja)
IL (1) IL85126A (ja)

Families Citing this family (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5022397A (en) * 1987-02-11 1991-06-11 General Electric Company Multiple-echo NMR angiography with multiple-velocity images
US5025788A (en) * 1987-02-11 1991-06-25 General Electric Company Method of acquiring NMR angiograms in selected flow component directions
US5022398A (en) * 1987-02-11 1991-06-11 General Electric Company Multiple-echo NMR angiography for enhancement of signal-to noise ratio
DE3726932A1 (de) * 1987-08-13 1989-02-23 Spectrospin Ag Verfahren zum kodieren von n parametern bei der mehrdimensionalen fourier-nmr-spektroskopie
US4978918A (en) * 1988-04-24 1990-12-18 Mitsubishi Denki Kabushiki Kaisha Magnetic resonance imaging method
DE3823961A1 (de) * 1988-07-15 1990-01-18 Philips Patentverwaltung Kernspintomographieverfahren und kernspintomograph zur durchfuehrung des verfahrens
JPH0685769B2 (ja) * 1988-09-09 1994-11-02 富士電機株式会社 磁気共鳴イメージング装置
US4952876A (en) * 1988-11-23 1990-08-28 General Electric Company Variable bandwidth multiecho NMR imaging
US5115812A (en) * 1988-11-30 1992-05-26 Hitachi, Ltd. Magnetic resonance imaging method for moving object
US4918386A (en) * 1988-12-23 1990-04-17 General Electric Company Method for simultaneously obtaining three-dimensional NMR angiograms and stationary tissue NMR images
US4888552A (en) * 1989-01-05 1989-12-19 Elscint Ltd. Magnetic resonance imaging
DE4024161A1 (de) * 1989-08-11 1991-02-14 Siemens Ag Pulssequenz zur schnellen ermittlung von bildern der fett- und wasserverteilung in einem untersuchungsobjekt mittels der kernmagnetischen resonanz
US4973906A (en) * 1989-08-17 1990-11-27 General Electric Company Flow compensated NMR fast pulse sequence
AT399646B (de) * 1989-12-01 1995-06-26 Wach Paul Dipl Ing Dr Techn Einrichtung zur erfassung des durchflusses des blutgefässsystems
JPH03228745A (ja) * 1990-02-05 1991-10-09 Toshiba Corp Ｍｒｉ装置
EP0447970B1 (en) * 1990-03-20 1997-09-10 Kabushiki Kaisha Toshiba Method and apparatus for imaging blood vessels employing magnetic resonance
US5221898A (en) * 1990-11-30 1993-06-22 Hitachi, Ltd. Flow imaging method using an MRI apparatus
US5133357A (en) * 1991-02-07 1992-07-28 General Electric Company Quantitative measurement of blood flow using cylindrically localized fourier velocity encoding
JP3161750B2 (ja) * 1991-06-05 2001-04-25 株式会社日立製作所 磁気共鳴診断装置および画像データ処理方法
US5233302A (en) * 1991-06-21 1993-08-03 General Electric Company Masking motion ghost artifacts in NMR images
US5270654A (en) * 1991-07-05 1993-12-14 Feinberg David A Ultra-fast multi-section MRI using gradient and spin echo (grase) imaging
US5281914A (en) * 1991-08-09 1994-01-25 The Johns Hopkins University Method of vector magnetic resonance measurement and imaging and associated apparatus
US5233298A (en) * 1992-02-20 1993-08-03 General Electric Company Quantitative measurement of blood flow at multiple positions using comb excitation and fourier velocity encoding
US5277192A (en) * 1992-09-18 1994-01-11 General Electric Company Imaging of turbulence with magnetic resonance
DE4432575C2 (de) * 1993-09-14 2003-04-10 Toshiba Kawasaki Kk Verfahren zur Bildgebung der Gehirnfunktion mittels einer Kernspinresonanz-Vorrichtung und hieran angepasste Kernspinresonanz-Vorrichtung
US6122540A (en) * 1994-06-08 2000-09-19 The Regents Of The University Of California Noninvasive measurement of renal hemodynamic functions using magnetic resonance imaging
US5602891A (en) * 1995-11-13 1997-02-11 Beth Israel Imaging apparatus and method with compensation for object motion
US6188922B1 (en) * 1999-01-08 2001-02-13 Wisconsin Alumni Research Foundation Phase contrast imaging using interleaved projection data
WO2002067202A1 (en) * 2001-02-16 2002-08-29 The Government Of The United States Of America, Represented By The Secretary, Department Of Health And Human Services Real-time, interactive volumetric mri
JP4309632B2 (ja) * 2002-10-08 2009-08-05 株式会社東芝 磁気共鳴イメージング装置
US7646198B2 (en) * 2007-03-09 2010-01-12 Case Western Reserve University Methods for fat signal suppression in magnetic resonance imaging
WO2009094304A2 (en) * 2008-01-23 2009-07-30 The Regents Of The University Of Colorado Susceptibility weighted magnetic resonance imaging of venous vasculature
KR101253024B1 (ko) * 2011-05-06 2013-04-16 연세대학교 산학협력단 자기공명영상을 이용한 삼차원 영상 복원장치 및 방법
US11191503B2 (en) 2018-07-17 2021-12-07 International Business Machines Corporation Fluid-injector for a simultaneous anatomical and fluid dynamic analysis in coronary angiography

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58165046A (ja) * 1981-12-21 1983-09-30 アルバ−ト・マコフスキ 核磁気共鳴を用いた血管の投影像形成方法及び装置
JPS6188132A (ja) * 1984-08-10 1986-05-06 ザ・ゼネラル・ホスピタル・コ−ポレ−シヨン 核磁気共鳴の手段による流動的画像形成
JPS61119253A (ja) * 1984-03-10 1986-06-06 ゼネラル・エレクトリツク・カンパニイ 磁気共鳴作像方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4458203A (en) * 1980-12-11 1984-07-03 Picker International Limited Nuclear magnetic resonance imaging
US4443760A (en) * 1982-07-01 1984-04-17 General Electric Company Use of phase alternated RF pulses to eliminate effects of spurious free induction decay caused by imperfect 180 degree RF pulses in NMR imaging
JPS5940843A (ja) * 1982-08-31 1984-03-06 株式会社東芝 診断用核磁気共鳴装置
US4602641A (en) * 1983-08-15 1986-07-29 The Regents Of The University Of California Method and apparatus for NMR detection and imaging of flowing fluid nuclei
US4595879A (en) * 1983-11-14 1986-06-17 Technicare Corporation Nuclear magnetic resonance flow imaging
DE3504734C2 (de) 1985-02-12 1998-12-10 Max Planck Gesellschaft Verfahren und Vorrichtung zur Aufnahme von Spinresonanzdaten
US4654591A (en) * 1985-07-29 1987-03-31 Wisconsin Alumni Research Foundation NMR flow imaging using bi-phasic excitation field gradients
US4714081A (en) * 1986-03-03 1987-12-22 General Electric Company Methods for NMR angiography

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS58165046A (ja) * 1981-12-21 1983-09-30 アルバ−ト・マコフスキ 核磁気共鳴を用いた血管の投影像形成方法及び装置
JPS61119253A (ja) * 1984-03-10 1986-06-06 ゼネラル・エレクトリツク・カンパニイ 磁気共鳴作像方法
JPS6188132A (ja) * 1984-08-10 1986-05-06 ザ・ゼネラル・ホスピタル・コ−ポレ−シヨン 核磁気共鳴の手段による流動的画像形成

Also Published As

Publication number Publication date
FI880585A0 (fi) 1988-02-09
EP0278271A2 (en) 1988-08-17
IL85126A (en) 1996-06-18
FI880585A (fi) 1988-08-12
IL85126A0 (en) 1988-06-30
US4796635A (en) 1989-01-10
JPS63226343A (ja) 1988-09-21
EP0278271A3 (en) 1990-06-13

Similar Documents

Publication Publication Date Title
JPH0424051B2 (ja)
US4714081A (en) Methods for NMR angiography
US7706857B2 (en) Methods and apparatus for mapping internal and bulk motion of an object with phase labeling in magnetic resonance imaging
EP0498570B1 (en) Quantitative measurement of blood flow using cylindrically localized fourier velocity encoding
EP2699926B1 (en) Spatially encoded phase-contrast mri
US5034692A (en) Magnetic resonance imaging method for acquiring flux-compensated, T2 -w
US8148984B2 (en) Method for magnitude constrained phase contrast magnetic resonance imaging
US5588431A (en) Inversion recovery MRI
US5038783A (en) Multiple gradient echo pulse sequence for acquisition of NMR angiograms
JPH0420617B2 (ja)
US20080048657A1 (en) Method for aligning multiple MR images acquired with alternating readout gradient
US5038784A (en) Multiple-echo angiography with enhanced signal-to-noise ratio
JPH0263437A (ja) ３次元核磁気共鳴造影図を作る方法と装置
US4918386A (en) Method for simultaneously obtaining three-dimensional NMR angiograms and stationary tissue NMR images
US5022397A (en) Multiple-echo NMR angiography with multiple-velocity images
US7233818B1 (en) Methods and apparatus for mapping internal and bulk motion of an object with phase labeling in magnetic resonance imaging
US5092335A (en) Method for the imaging of intra-voxel movements by nmr in a body
JPH0616766B2 (ja) Ｎｍｒイメージング装置を使用した高速流れ測定方式
US5025788A (en) Method of acquiring NMR angiograms in selected flow component directions
JPH0373300B2 (ja)
US4855679A (en) Magnetic resonance studies of restricted volumes
US5022398A (en) Multiple-echo NMR angiography for enhancement of signal-to noise ratio
KR930000103B1 (ko) 핵 자기 공명 혈관 조영 영상 획득 장치
JPH0622934A (ja) 磁気共鳴イメージング装置及び方法
EP4097498B1 (en) Mr imaging using dixon-type water/fat separation with suppression of flow-induced leakage and/or swapping artifacts