JPH03184530A - Nmr装置 - Google Patents

Nmr装置

Info

Publication number
JPH03184530A
JPH03184530A JP2276578A JP27657890A JPH03184530A JP H03184530 A JPH03184530 A JP H03184530A JP 2276578 A JP2276578 A JP 2276578A JP 27657890 A JP27657890 A JP 27657890A JP H03184530 A JPH03184530 A JP H03184530A
Authority
JP
Japan
Prior art keywords
nmr
signal
motion
data
pulse sequence
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
Application number
JP2276578A
Other languages
English (en)
Other versions
JPH0622494B2 (ja
Inventor
Gary H Glover
ガリー・ハロルド・グローバー
Stephen W Flax
ステファン・ウエイン・フラックス
Ann Shimakawa
アン・シマカワ
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 JPH03184530A publication Critical patent/JPH03184530A/ja
Publication of JPH0622494B2 publication Critical patent/JPH0622494B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5676Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/541Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Physiology (AREA)
  • Power Engineering (AREA)
  • Pulmonology (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Signal Processing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 発明の背景 この発明は核磁気共鳴作像方法に関する。更に具体的に
云えば、この発明は例えばNMR走査の進行中の被検体
の動きによる略周期的なNMR信号の変動が原因で起る
像の人為効果(アーチファクト)を制御する方法に関す
る。
NMRは、患者の解εり学的な特徴の像を求める様に開
発されている。こう云う像は、核スピンの分布(典型的
には水及び組織に関連する陽子)、スピン−格子緩和時
間T1、及び/又はスピン−スピン緩和1侍間T2を表
わし、医学的に診断価値がある。像を構成する為のNM
Rデータは、多重角度投影再生及びフーリエ変換(FT
)の杵な利用し得る多数の方式の1つを用いて、収集す
ること力咄来る。典型的には、こう云う方式は、逐次的
に求められる複数個のビュー(View)で構成すれた
パルス順序で構成される。各々のビュー(よ1回又は更
に多くのNMR実験を含んでいることがあり、各々の吏
験は少なくともRF励振)くルスと、その結果得られる
NMR信号に空間情報を符号化する為の磁界勾配パルス
とで構成される。周知の様に、NMR信号は自由誘導減
衰(FID)信号であってもよいし、或いはこの方が好
ましいが、スピン・エコー信号である。
次に、「スピン捩れ(spin−warp )形」とよ
く呼ばれる周知のFT方式の変形の場合について、この
発明の好ましい実施例を詳しく説明する。然し、この発
明の方法がFT作像方法に制限されず、米国特許第4.
471.’306号に記載されている多重角度投影再生
、並びに米国特許第4,070.611号に記載されて
いるFT方式の別の変形の様な他の方式と共に有利に実
施することが出来ることを承知されたい。スピン捩れ形
方式は、フィジックス・イン・メディスン・アンド・バ
イオロジ誌、第25巻第751頁乃至第756頁(19
80年)所載のW、 A、ニーデルシュタイン他の論文
「スピン捩れ形NMR作像と全身作像への応用」に記載
されている。
簡単に云うと、スピン捩れ形方式は、NMRスピン・エ
コー信号を収集する前に、可変銀輪の位相符号化磁界勾
配パルスを用いて、この勾配の方向に空間情報を位相符
号化する。例えば、2次元形式(2DFT)では、位相
符号化勾配(Gy)をその方向に印加することによって
、空間情報を一方の方向に符号化し、その後位相符号化
方向に対して直交する方向に磁界勾配(G1)が存在す
る状態でスピン・エコー信号を観11111する。スピ
ン・エコーの間に存在する勾配が、空間情報を直交方向
に?1月化する。典型的な2DFT/くルス順序ては、
位相符号化勾配パルスcyyの大きさを、像全体を再坐
し得る様な1組のNMRデータを発生する為に収集され
る一連のビューで、単調に増分(ΔG、)的に州側する
NMR像データの収集中に物体の変位があると、位相符
号化方向にぼけと「ゴースト」の両方が坐する。動きが
周期的であるか或いは大体周期的である侍、特に別々の
ゴーストが目につき、動きが非周期的である時、スミア
リングが起る。心臓及び呼吸の運動を含む大抵の生理学
的な動きでは、NMR信号の各々のビューは、収集窓の
間に物体が不動と見なしてもよい位短い期朋内に収集さ
れる。従って、ぼけ及びゴーストは主に物体がビュ毎に
一貫しない形で見えることによるものであり、特に動き
によるNMR信号の振幅並びに/又は位相の変化によっ
て起る。
データの収集を物体の機能的なサイクルと同期させれば
、ぼけもゴーストも少なくすることが出来る。この方法
がゲート式NMR走査の名前で知られており、その目的
は、各々のビューで物体が同じ様に「見える」様に相次
ぐ機能的なサイクルの間、同じ点のNMRデータを収集
することである。ゲート作用の欠点は、物体の機能的な
サイクルの小さな一部分の間しか、NMRデータを収集
することが出来ないことであり、受入れることの出来る
様な最も短いパルス順序を用いた11!iでも、このゲ
ート方式では、データ収集時間が大幅に長くなることが
ある。
ゴーストの人為効果を除く為に提案された1つの方法が
、米国特許第4.567.893号に記載されている。
この米国特許では、動きが周期的である時、ゴーストと
作像する物体との間の像の中での距離は、NMRパルス
順序の繰返し時間が、(米国特許第4,443,760
号に記載されている様に、1つのビュー当たりに2つの
位相交番のRF励振パルスを使う場合)周期的な信号の
変化の持続■、+i間の174の奇数倍である1j7、
最大になることが認識されている。この比を使って、呼
吸運動によるゴーストを軽減することが出来ることが判
っている。この方法は、実際に像の品質を改善するが、
これはNMRパルス順序の繰返し時間に拘束を加えるも
のであり、その結果合計の走査間開が一層長くなる場合
が多い。更にこれは、動きが周期的であると仮定してい
る。彼検体の呼吸が不規則である時、ゴーストがぼけ、
関心のある像領域と重なることがあるので、その効果が
減少する。
周期的な信号の変化による望ましくない影響を減少する
別の方法が、米国特許第4. 706. 026号に記
載されている。この方法の1例では、(例えば患者の呼
吸による)信号の変化の周期に成る仮定を下し、ビュー
の順序を普通の様に単調に増加する位相符号化勾配から
予め選ばれた順序に変更する。これは、勾配パラメータ
、即ち、(スピン捩れ形方法に於ける)位相符号化勾配
パルスの振幅又は読出し勾配パルスの方向(多重角度投
影再生方法に於ける)を用いる順序を定めることを必要
とする。所定の信号の変動の周期に対し、ビューの順序
は、位相符号化振幅(又は勾配の方向)の関数としての
NMR信号の変動が所望の周波数になる様に選ぶ。1丈
施例では、動きの受動周期が合計のNMR走査時間(低
周波数)と等しくなる様にビューの順序を選んで、ゴー
ストの人為効果が出来るだけ作像する物体に接近する様
にする。別の史施例(高周波数)では、変動周期を出来
るだけ短くして、ゴーストの人為効果を出来るだけ物体
から遠ざける様に、ビューの順序が選ばれる。
従来のこの方法は、人為効果を少なくする効果があり、
変動が割合規則的でその周波数がtl+つていれば、成
る点では理想的である。他方、動きの時間的な周期に対
して下した仮定が成立しない(例えば思考の呼吸パター
ンが変化するか或いは不規則である為に)場合、この方
法は余りよく成立しない。そう云うことが起ると、ゴー
ストを出来るだけ物体に近付けるとか或いは物体から出
来るだけ遠ざけると云う集束作用がぼける為に、その効
果の幾分かゾ失われる。この問題に対する解決策が米国
特許第4,663,591号に記載されている。この方
法では、走査を実行する時、単調でないビューの順序を
決定し、これは周期の変化に応答する様にして、信号の
変動と勾配パラメータの間に所望の関係(低周波数又は
高周波数)を生ずる様にする。勿論、この方法の有効性
は、患者の動き、特にその動きの周期の振幅変動を感知
するのに使われる手段の精度によってた右される。
N M R走査を受ける患者の呼吸運動を監視する為に
使われる普通の1つの手段は、患者の胸の周りに取付け
た圧力センサを有する空気圧ベローで構成される。患者
が呼吸すると、前脚及び胃壁が膨張並びに収縮し、その
結果ベローに起る圧力変化により、圧力変換器からの電
気信号の大きさが呼吸位相の関数として変化する。空気
圧ベローは連続的な保守を必要とする複雑な電気機械的
な装置であるばかりでなく、こう云うモニタによって発
生される信号は、NMR像に動きによる人為効果を招く
呼吸運動を必ずしも正確に表示しないことが判った。
呼吸を監視する別の方式は、米国特許第4,564.0
17号に記載されている様に、走査中に得られたNMR
信号を用いて、呼吸の速度を表わす信号を発生すること
である。この方法は、NMRエコー信号の大きさの変化
を測定することにより、呼吸の速度の「オン・ライン」
の目安を生ずるが、呼吸サイクルの位相を示す正確なオ
ン・ラインのデータとはならない。
発明の要約 この発明は、動きによる人為効果を抑圧する方式に使う
為、NMR走査の間、呼吸位相情報を求める改良された
装置に関する。具体的に云うと、この発明は、NMR走
査の間、NMR運動信号を発生する手段と、NMR運動
信号を処理して、走査される彼検体の特徴の位1dを示
す1組のデータを発生する手段と、1組のデータ中の特
定の特徴の位置を突ILめる手段と、特定の特徴の現在
の位置をその特徴の基準位置と相関させることによって
、呼吸信号を発生する手段とを含む。
この発明の全般的な目的は、NMR走査の間、解ill
学的な描込の動きの一層正確な表示を発生することであ
る。動きを検出するl′bの1 illのNMRデータ
は、像の為の1組のN MRデータを収集するのに使わ
れるのと同じ走査の間に収集され、全く同し装置を用い
て収集される。呼吸運動を監視するのに全く同じNMR
装置を使うことにより、像を再生するハに使われるNM
R信号に対する動きの影響を一層iE 6&に測定する
ことが出来る。多数の、動きによる人為効果の抑圧力l
ノ、のどれに使った時でも、この改良された呼吸位相信
号により、曹通の呼吸モニタからの位相情報を用い?:
*rr生された像よりも、−層よい像を再生することが
出来る。
この発明の別の目的は、帯紐によって患者に取イ;1け
なければならない電気機械的な装置の必要をなくすこと
である。この発明は、像の為のNMRデータを収集する
のに必要な装置以外の装置を何″、9必要としない。N
MRパルス順序は、動きを検出する為の1組のN M 
Rデータを収集する様に女史されるが、このハに余分の
装置を必要としない。
患者に対するIIIF−の影響は、走査全体を文施する
のに必要な時間が若干長くなることである。
この発明の上記並びにその他の目的及び利点は、以下説
明する所から明らかになろう。この説明では、この発明
の好ましい実施例を示した図面を参j1(くする。然し
、この実施例は必ずしもこの発明の範1J11全体を表
わすものではなく、この発明の範囲の解釈に当たっては
特許請求の範囲を参照されたい。
好ましい起径例の説[jll 第ii”4+よ1.・ハ仔凹力好ましい実施例をL12
いるNMR作像装置の簡略ブロック図である。装置はパ
ルス制御モジュール112を持ち、これがホスト・コン
ピュータ114の制御のもとに、包括的に116に示し
た磁界勾配電源に文・lする正しいタイミングのパルス
波形信号を供給する。電源116が、全体をブロック1
18で示した勾配コイル集成体の一部分をなす勾配コイ
ルを付勢する。集成体は、デカルト座標系のx、  y
、  z方向の向きを持つ磁W勾配C,,Cy、C2を
夫々発生するコイルを含む。NMR作像の用途で勾配C
,,Cy。
G2をどの様に使うかは、後で第2図について説明する
第1図の説明を続けると、パルス制御モジュール112
がRF合成器120に作動パルスを供給する。この合成
基はRFトランシーバの一部分であり、その一部が破線
のブロック122の中に囲い込まれている。パルス制御
モジュール112は、RF周波数合成器120の出力を
変調する変21器124にも信号を供給する。変、iR
F信号が、RF電力僧幅器128及び送信/受信スイッ
チ130を介してRFコイル集成体126に印加される
RF倍信号使って、作像するサンプル物体(図面に示し
てない)内の咳スピンを励振する。
励振された核スピンからのNMR信号がRFコイル集成
体126によって感知され、送信/受信スイッチ130
を介してRF前置増幅器132に印加される。増幅され
たNMR信号が直角位相検波器134に印加され、検波
された信号がA/D変換器136によってディジタル化
され、周知の様に記憶並びに処理の為に、コンピュータ
114に印加される。
第2図は2次元フーリエ変換(2DFT)と呼ばれるが
、2次元「スピン捩れ形」とも呼ばれる従来の作像パル
ス順序の2つのビューを示す。このパルス順序は、周知
の様に、検査する物体の像を再生する為の作像NMRデ
ータを求めるのに役立つ。2つのビューを“A”及びB
′で示してあり、これらは位相符号化勾配磁界Gyの他
は同一である。各々のビューは、位相交番RF励振パル
スを利用するパルス順序であり、こう云う位相交番RF
励振パルスにより、米国特許第4,443,760号に
記載されている様に、位相交番のNMR信倍信+  (
t)及びS+’  (t)が発生され、NMR装置に於
ける成るベースライン誤差を相殺する。
第2図のビューAについて説明すると、期間1(横軸に
沿って示す)に、正のG!磁界勾配パルスの7’j (
Lのちとに印加される選択性90°RF励振パルスが示
されている。パルス制御モジュール112(第1図)が
、周波数合成櫓120及び変、LJe124に必要な制
御信号を供給して、その結果iすられる励振パルスが、
作像する物体の予定の領域にある核スピンだけを励振す
るのに正しい位相及び周波数になる様にする。典型的に
は、励振パルスは(sln  x) /x関数によって
振幅変調することが出来る。合成器120の周波数は、
周知のラーモアの方程式に従って、印加された分極磁W
の強度と作像される特定のN M R種口とに関係スる
。パルス制御モジュール112は勾配電源116にも作
動信号を印加して、この場合はG:!勾配パルスを売上
する。
第2図の説明を続けると、G、、G、及びG工勾配パル
スが期間2に同時に印加される。期間2のG2勾配は位
相戻しパルスであって、典型的には期間2にわたる勾配
波形の時間積分が、期間1にわたるG2勾配波形の間開
積分の一1/2倍に大体等しくなる様に選ばれる。負の
Gエバルスの作用は、期間1に励振された核スピンの位
相戻しをすることである。Gy勾配パルスは、勾配の方
向に空間情報を符号化する為に、ビューA、B。
・・・・・・等の各々で異なる振幅を持つ様に選ばれた
位相符号化パルスである。相異なるGy勾配の振幅の数
は、典型的には、再生像が位相符号化方向(Y)に持つ
画素分解要素の数と少なくとも等しくなる様に選ばれる
。典型的には、128個、256個又は512個の異な
る勾配Gyの振幅が選ばれ、典型的なNMR装置では、
NMR走査が完了するまで、Gyの値は、あるビューと
次のビューで、一定ヱた1プ僧加する。
期間2の61勾配パルスは、スピン・エコー信号S+ 
 (t)の発生時刻を期間4に遅延させる為に、励振さ
れた核スピンを予定工だけ位相外しするのに必要な位相
外しパルスである。スピン・エコー信号が、典型的には
期間3に180’RFパルスを印加することによって発
生される。周知の様に、180″RFパルスは、スピン
・エコー信号を発生する様にスピンの位相外しの方向を
逆転するパルスである。スピン・エコー信号が、期間4
に、勾配パルスG1の存在のもとに標本化され、この勾
配の方向(X)に空間情報を符号化する。
前に述べた様に、各々のビューで追加のNMR7II1
1定を用いることにより、ベースライン誤差成分が除か
れる。この2回目の測定は、ビューAの期+::+Sの
RF励振パルスが、ビューAの期間1の励振パルスに比
べて(負の符号で示す様に)1000位…がずれる様に
選ばれることを別とすると、1回11と路間−である。
その結果、期間8のスピン・エコー信号s + ’  
(t )は、期間4のスピン・エコー信号S+  (t
)に対して1806位相がずれている。信号S+’  
(t)をSt  (t)から減算すれば、信号の内、信
号S+   (t)で符号が逆転した成分だけが残る。
この2′51ベ一スライン誤差成分が相殺される。
k:ユAについて上に述べた過程が、位相符号化勾配G
yの全ての振幅に対し、ビューB等で繰返される。この
走査の間に収集されたNMR像デタをホスト・コンピュ
ータ114に記憶し、そこで処理して、CRT表示装置
を制御するのに適した像データを発生する。
−にに述べた普通のNMR走査を実施するn;f、作像
される物体の平面又はスライス内の物理的な全ての場所
からのNMRデータが収集される。正確な像を再生しよ
うとする場合、N M R走査全体を完了するのに必要
な時間の間、物体及び測定条件の両方が安定であるか、
固定でなければならない。
この究明は、そうならない実際的な場合を取上げるもの
である。実際的には、測定条件は成る循環的な或いは循
環的に近い形で変化する。
人体の腹を通る像を発生する場合、この様な状態が起る
。この場合、作像する物質の大部分は、人体の呼吸のれ
に動いており、像全体の為のNMRデータを収集するの
に要する旧聞が、多数の呼吸サイクルを越える場合が多
い。呼吸サイクル全体にわたってNMRデータを連続的
に収集すると、ビュー毎に肢検体の配置が異なり、再坐
像は動きによる多くの人為効果を持つことになる。
この発明は、呼吸による披検体の運動を示すNMRデー
タを収集することが出来る様にすると共に、このNMR
運動データを処理して呼吸モニタ信号を発生する様なN
MR装置の改良である。前に述べた様に、呼吸モニタ信
号は、公知の多数の、動きによる人為効果を抑圧する方
式のどれにも用いることが出来る。
走査の間のNMR運動情報を収集するハ、第2図の像デ
ータ・パルス順序の各々のビューの間に、別個のパルス
順序をインターリーブ形に入れる。
この運動用のNMRパルス順序か第3図に示されており
、これは第4図に示す様に、像データ・パルス順序の各
々のビューの直前に実行される。第3図について具体的
に説明すると、運動用NMRパルス順序は像用データ・
パルス順序とは幾つかの点で異なる。例えば、スライス
中のスピンを励振する為に、RF励振パルス141を印
加する間、スライス選択勾配パルス140を印加するが
、位相符号化勾配パルスは印加しない。その代りに、N
MR信号142は読出し勾配パルス143の間に収集す
る。好ましい実施例では、パルス143は磁W勾配G、
によってY軸に沿って発生される。
X軸磁W勾配G、が印加されず、従って、この結果得ら
れるNMR信号はX軸に沿った位置情報を持っていない
。この発明の好ましい実施例では、RF励振パルス14
1がスピンに対して90°のはしき(Nip)角を加え
、スライス選択勾配パルス140が、像データを収集す
るのと同じスライスを選択する。
走査の間に収集されたデータが第6図に示されている。
収集された像データが2次元配列50に記憶される。こ
の配列は、発生しようとする像の選ばれた分解能によっ
て決定された寸法(例えば128X128.256X2
56又は512X512個の画素)を持つ。当業者であ
れば、実際にはこの様な配列を2つ収集することが理解
されよう。1つは同相収集データ(I)を持ち、他方は
直角位t0データ(Q)を持つ。配列50にあるブタの
各行が、1つのビューの間に収集されたディジタル化さ
れたデータを表わし、各々の行は、位相7:)号化磁昇
勾配Gyの異なる値を用いて収集される。位#+1?:
)9化勾配Gyの各々に対するデータが収集された哨、
走査が完了し、配列50を処理して、表示装置に於ける
各々の画素の輝度を示す像配列51を作成する。好まし
い実施例では、この処理は、配列50にある像データの
読出し次元に沿った第1のフーリエ変換を実施し、その
後こうして得られたデータ配列に位相符号化次元に沿っ
た第2のフーリエ変換を実施することを含む。
この結果得られたデータ配列中の要素の値が、像配列5
1に於ける対応する画素の輝度を示す。
配列50にある像データの各々の行を収集する直前に、
第3図の運動用NMRパルス順序を用いて、運動データ
の1次元配列52を収集する。第5図及び第6図につい
て説明すると、この運動データは作像するスライス54
から収集され、垂直又はY軸方向の所望の位置7II1
1定の分解能が得られる様にするのに十分な数のサンプ
ル(例えば、128個、256個、512個等)を持っ
ている。
好ましい実施例では、運動データは、像データと同じ垂
直及び水平の視野から収集され、同じ数のサンプルを有
する。然し、これはこの発明を実施するのに必要なこと
ではない。
運動データの各々の配列52を収集した時、それを直ち
に処理して動きを表わす信号を発生する。
この信号は、例えば、次のビューで使われる位相符号化
勾配Gyの値を決定するハに使うことか出来る。好まし
い実施例では、監蜆される動きは、前側腹壁の動きであ
る。これは、呼吸の関数として、身体の他の描込の動き
を非常に正確に示すと共に、収集された運動用NMRデ
ータを確認するのが比較的容易であるからである。腹壁
は、NMR信号にピークを発生する脂肪層を有し、これ
は殆んど信号を発生しない患者の身体の外側の空気に隣
接している。
第3図のパルス順序が、運動用NMR信号を収集するl
′hの好ましい手段であるが、この他の代りも考えられ
るし、場合によっては好ましいことがある。例えば、運
動用パルス順序によって選ばれるスライスは、像パルス
順序によって選択されるスライスと同じである必要はな
く、スライスは同じ軸に沿った向きである必要さえない
。909RF励振パルスを使うことが出来るが、はじき
角(x)を−IW小さくしてもよく、それでも作像用パ
ルス順序との相互作用が一層少なくなる。第3図の運動
用パルス順序では180°RFパルスを用いているが、
勾配パルス143の前縁に負のローブを印加することに
より、勾配エコー信号を発生することも可能である。1
80°RF励振によって発生されるエコーに較べた、勾
配によって誘起されるエコーの利点は、運動用NMRデ
ータを収集するのに、時間がかなり短縮されることであ
る。実際、像データと同じ順序で運動データを収集すれ
ば、合計走査時間は更に短縮することが出来る。例えば
、勾配によって誘起されたエコー信号は、第2図の期間
2で、位相符号化勾配パルスの直前に、又は第2図の期
間4で、エコー信号S+  (t)を収集した直後に収
集することが出来る。
この発明による運動データの処理が、記憶されているプ
ログラム命令に応答して、コンピュータ114によって
実行される。こう云うプログラム命令が第8図にフロー
チャートで示されている。
運動NMRデータを処理する肋の最初の工程は、配列5
2にあるデータの1つの次元、即ち垂直読出し次元に沿
って、プロセス・ブロック65で示す様に、フーリエ変
換を史施することである。
この変換は要素256の配列に対しては、数学的に次の
様に表わすことが出来る。
こ\でS はデータ配列52の要素、S 、は変換y されたデータ配列の対応する要素である。
その後、プロセス・ブロック66で、運動データのこう
して得られた1次元配列の大きさを計算して、前側腹壁
の物理的な場所の表示を発生する。
この計算は、変換された配列の各々の要素に対し、数学
的に次の様に表わすことが出来る。
1m 、 −(1=)  +(Q 、)    (2)
x        y          yこ\でI
 、は変換された配列の要素S 、の同相のy 値、Q 、は変換された配列の要素S 、の直角位相y の値である。
この様な変換された運動データの1例が第7図に示され
ており、グラフ60は、最初の運動データ・パルス順序
に対する変換されたデータの大きさを示し、グラフ61
は同じ走査の間にこの後収集された運動データに対して
変換されたデータの大きさを示す。最初のグラフ60を
基準として使い、この後の各々のグラフ61の位置はこ
の基準に対して測定する。例えば、前側腹壁によって発
生された信号60のピークが62に示されており、この
後の1組のNMR運動データの中にある同じピークが6
3に示されている。これらのピークの位置のずれ(ΔY
)が、前側腹壁の運動の直接的な測定値であり、この量
に倍率をかけて、呼吸モニタのディジタル出力に直接的
に代るものにする。
第1図に戻って説明すると、上に述べた運動NMRデー
タがコンピュータ114に印加され、そこで呼吸モニタ
信号ΔYを発生する様に処理される。こう云う呼吸モニ
タ信号は、内部で人為効果を抑圧する為に使ってもよい
し、或いは64から出力して外部で使ってもよい。例え
ば、64に出力されるΔY倍信号、前に引用した米国特
許第4゜663.591号の第18図に示される呼吸計
変換器、ま曽幅器及びA/D変換器の代りに使うことが
出来る。
第7図及び第8図に戻って説明すると、位置のずれΔY
の測定は2つの工程に分けて行なわれる。
即ち、窓ぎめ工程と相互相関工程とである。窓ぎめは、
配列52にある運動データから、前側腹壁によって発生
された信号のピークを含む小さな部分又は窓を隔離する
様に作用する。この窓が第7図に示されており、これは
数ある周知の方法の任意の1つによって選ぶことが出来
る。好ましい実施例では、第8図のプロセス・ブロック
70に示す様に、ポインタを定めて、変換された配列I
m、中の最後の要素を読出す。この要素をスライス54
のてっぺんに於けるスピン密度と関係づけるが、典型的
には、この場所には空気しか存在しない2′)に非常に
小さい。その後、ループに入り、プロセス・ブロック7
1で示す様に、i11定ブロック72で211定して予
め設定した閾値を越える1つの要素が見つかるまで、変
換された配列の相次ぐ要素を検品する。ブロック73で
ポインタをデクレメントして、探索がスライス54のて
っぺんから始まり、前側腹壁を示す信号に出会うまで、
下向きに作業を続ける。出会った時、変換された配列の
中にあって、選ばれた窓の外側にある全ての値は、プロ
セス・ブロック74て示す様にゼロに設定される。好ま
しい実施例では、要素60個の窓寸法を使い、予め設定
した閾値を最初に越えた要素を中心とする。
次の工程は、基準のグラフ60から得られた窓ぎめして
求めたデータを、現在のグラフ61からの窓ぎめして求
めたデータと相関させることである。一方の波形を他方
の波形と出来るだけ密接に相関させ又は整合させる為に
、一方の波形をどれだけずらさなければならないかを決
定する相関方式は、数多く知られている。第7図及び第
8図について云うと、好ましい実施例では、プロセス・
ブロック75に示す様に、窓ぎめした両方の組のデータ
を最初に微分し、その後プロセス・ブロック76に示す
様に、微分しなかった値に定数を加えたもので除して、
投影の縁を表わす、振幅が小さくて変化の激しい信号に
有利となる様に、加重する。
こ\でCは基準の運動データの組のピーク値の一定の分
数である。この基準波形の窓で得られた要素(q 、)
をに右にずらして、その要素と、現在の運動データの窓
ぎめした部分にある対応する要素との積の和が、最大に
なる位置をプロセス・ブロック77で示す様に決定する
。云い換えれば、窓ぎめした基準要素を中心の左から右
までに及ぶ一連の位置にわたってずらし、各々の位置で
、窓ぎめした2つのデータ配列の対応する要素を乗算し
て加算する。こうして計算された値が、相互相関曲線を
定め、この曲線は、窓ぎめした基準データを相関が最も
よい位置へずらした時に、ピークになる。ずらした各々
の位置が、Y軸に沿った距離の増分を表わし、中心から
のずらしの合計の回数が、基準位置からの前側腹壁の所
望の変位(ΔY)を示す。更に高い分解能を希望する場
合、補間方法を用いて、通常は2つの別々のずらし位置
の間にある相関山1線のピークを決定することが出来る
。何れにせよ、変位の値ΔYが、プロセス・ブロック7
8に示す様に、現在の腹壁のt−u 置を表示するもの
として出力され、それをコンピュータ114西で用いて
、人ハ効果を抑圧する方式に利用してもよいし、或いは
前に述べた様に使う為に出力してもよい。
【図面の簡単な説明】
第1図はこの発明を用いたNMR装置の電気的なブロッ
ク図、第2図は第1図の装置が実行する作像パルス順序
の1例を示すグラフ、第3図は第2図の作像パルス順序
とインターリーブ形に使われるこの発明の運動パルス順
序の好ましい1実施例を示すグラフ、第4図は運動パル
ス順序と像パルス順序がインターリーブになる様子を示
す11、r間線図、第5図は作像される披検体の見取図
で、腹部を通るスライスを示しており、このスライスか
らのNMRデータを収集する。第6図はデータ収集及び
処理工程の間に得られるデータ構造を示す見取図、第7
図は所望の変位出力の値を発生するのに必要な相関を例
示する2組の処理済み運動ブタを示すグラフ、第8図は
この発明を実施するのに用いられるプログラムのフロー
チャー1・である。 [主な符号の説明] 51;像配列、 52:運動データ配列、 60゜ 62゜  12 14 22 26 61、運動データのグラフ、 63、ピーク、 :パルス制御モジュール、 :ボス!・・コンピュータ1 、I・ランシーバ、 RFコイル集成体。 特1;′「出願人

Claims (1)

  1. 【特許請求の範囲】 1、被検体のNMR走査の間、呼吸位相情報を発生する
    装置に於て、NMR像データの走査の間、NMR運動パ
    ルス順序を発生する手段と、該NMR運動パルス順序に
    よって発生されたNMR運動信号を収集する手段と、該
    NMR運動信号を処理して、被検体に於ける特徴の位置
    を示す1組のデータを発生する手段と、該1組のデータ
    中の特定の特徴の位置を突止める手段と、特定の特徴の
    現在の位置を該特徴に対する基準位置と相関させること
    によって呼吸信号を発生する手段とを有する装置。 2、NMR運動信号を処理する手段が、それに沿った特
    徴の位置を測定する軸に対応する次元に沿ってフーリエ
    変換を行なう手段を含む請求項1記載の装置。 3、NMR運動パルス順序を発生する手段が、被検体中
    にあるスピンに横方向磁化を発生するRF励振パルスを
    発生する手段と、前記軸に沿って磁界勾配を発生する手
    段とを含み、NMR運動信号を収集する手段が、前記磁
    界勾配を発生している間、前記NMR運動信号を受取っ
    てディジタル化する手段を含む請求項2記載の装置。 4、被検体の特定の特徴が人体の前側腹壁である請求項
    1記載の装置。 5、NMR走査が一連のNMR像パルス順序で構成され
    、NMR運動パルス順序が各々のNMR像パルス順序の
    前に発生され、それで収集したデータを使って、次に続
    くNMR像パルス順序を発生する前に前記呼吸信号を発
    生する請求項1記載の装置。 6、NMR走査が一連のNMR像パルス順序で構成され
    、各々のNMR像パルス順序に続いてNMR運動パルス
    順序が発生され、それで収集したデータを使って、次に
    続くNMR像パルス順序を発生する前に、前記呼吸信号
    を発生する請求項1記載の装置。
JP2276578A 1989-10-27 1990-10-17 Nmr装置 Expired - Lifetime JPH0622494B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US07/427,401 US4994743A (en) 1989-10-27 1989-10-27 Method for monitoring respiration with acquired NMR data
US427401 1989-10-27

Publications (2)

Publication Number Publication Date
JPH03184530A true JPH03184530A (ja) 1991-08-12
JPH0622494B2 JPH0622494B2 (ja) 1994-03-30

Family

ID=23694712

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2276578A Expired - Lifetime JPH0622494B2 (ja) 1989-10-27 1990-10-17 Nmr装置

Country Status (2)

Country Link
US (1) US4994743A (ja)
JP (1) JPH0622494B2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008307287A (ja) * 2007-06-15 2008-12-25 Hitachi Medical Corp 磁気共鳴イメージング装置

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9025344D0 (en) * 1990-11-21 1991-01-02 Marconi Gec Ltd Nuclear magnetic resonance methods and apparatus
US5200700A (en) * 1990-11-30 1993-04-06 General Electric Reduction of NMR artifacts caused by time varying linear geometric distortion
US5245283A (en) * 1991-08-07 1993-09-14 Picker International, Inc. Technique for shifting out-of-slice artifacts to the edge of the field of view
JPH0775627A (ja) * 1993-06-11 1995-03-20 Hitachi Ltd 磁気共鳴診断装置における体動追従計測方法
JPH0723945A (ja) * 1993-07-07 1995-01-27 Toshiba Corp 医用画像撮影装置
US5800354A (en) * 1994-11-23 1998-09-01 U.S. Phillips Corporation Method of and device for magnetic resonance imaging
US5613492A (en) * 1995-03-31 1997-03-25 New York University Method and product for improved images in magnetic resonance imaging using multiple breatholding
DE19811376C1 (de) * 1998-03-16 2000-01-05 Siemens Ag Verfahren zur zeitlich hochaufgelösten Magnetresonanztomographie
US6675036B2 (en) 2001-07-18 2004-01-06 Ge Medical Systems, Inc. Diagnostic device including a method and apparatus for bio-potential noise cancellation utilizing the patient's respiratory signal
JP3639825B2 (ja) * 2002-04-03 2005-04-20 キヤノン株式会社 動画像表示方法、プログラム、コンピュータ可読記憶媒体、及び動画像表示装置
JP4133348B2 (ja) * 2003-01-07 2008-08-13 株式会社日立メディコ 核磁気共鳴を用いた検査装置
US20110275933A1 (en) * 2008-10-22 2011-11-10 University Of Massachusetts Motion compensation in nuclear imaging
WO2011034538A1 (en) * 2009-09-18 2011-03-24 Analogic Corporation Rf power transmitter

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63164943A (ja) * 1986-09-03 1988-07-08 株式会社日立製作所 Nmrイメ−ジング方式
JPS6432851A (en) * 1987-07-30 1989-02-02 Mitsubishi Electric Corp Magnetic resonance imaging method
JPS6434342A (en) * 1987-07-30 1989-02-03 Mitsubishi Electric Corp Magnetic resonance imaging method
JPH01254156A (ja) * 1988-04-01 1989-10-11 Hitachi Ltd Nmrイメージング装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
USRE32701E (en) * 1983-01-04 1988-06-21 Wisconsin Alumni Research Foundation NMR scanner with motion zeugmatography
US4564017A (en) * 1984-11-21 1986-01-14 General Electric Company Method and apparatus for respiration monitoring with an NMR scanner
US4567893A (en) * 1984-11-21 1986-02-04 General Electric Company Method of eliminating breathing artifacts in NMR imaging
US4724386A (en) * 1985-09-30 1988-02-09 Picker International, Inc. Centrally ordered phase encoding
US4710717A (en) * 1986-12-29 1987-12-01 General Electric Company Method for fast scan cine NMR imaging
US4761613A (en) * 1987-08-12 1988-08-02 Picker International, Inc. Monitored echo gating for the reduction of motion artifacts
US4857846A (en) * 1988-03-31 1989-08-15 The Regents Of The University Of California Rapid MRI using multiple receivers producing multiply phase-encoded data derived from a single NMR response
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63164943A (ja) * 1986-09-03 1988-07-08 株式会社日立製作所 Nmrイメ−ジング方式
JPS6432851A (en) * 1987-07-30 1989-02-02 Mitsubishi Electric Corp Magnetic resonance imaging method
JPS6434342A (en) * 1987-07-30 1989-02-03 Mitsubishi Electric Corp Magnetic resonance imaging method
JPH01254156A (ja) * 1988-04-01 1989-10-11 Hitachi Ltd Nmrイメージング装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008307287A (ja) * 2007-06-15 2008-12-25 Hitachi Medical Corp 磁気共鳴イメージング装置

Also Published As

Publication number Publication date
US4994743A (en) 1991-02-19
JPH0622494B2 (ja) 1994-03-30

Similar Documents

Publication Publication Date Title
EP0218838B1 (en) Method and apparatus for reducing image artifacts due to periodic signal variations in nmr imaging
CN102078196B (zh) 借助磁共振采集相应于对象的呼吸运动的信号数据的方法
US4710717A (en) Method for fast scan cine NMR imaging
EP0212526B1 (en) Method for reducing image artifacts due to periodic signal variations in nmr imaging
JP2974759B2 (ja) Nmr画像における動きアーチファクトの低減装置および方法
US4614195A (en) Method for reduction of motion artifacts in Fourier transform NMR imaging techniques
US4751462A (en) Method for acquiring NMR data which is subject to periodic variations
US5195525A (en) Noninvasive myocardial motion analysis using phase contrast mri maps of myocardial velocity
US5897496A (en) Method and apparatus for producing magnetic resonance angiogram
US7372269B2 (en) Magnetic resonance imaging method and apparatus
CN102597795B (zh) 使用导航器的mr成像
US4730620A (en) Nuclear magnetic resonance method and apparatus for reducing motion artifacts
US4994743A (en) Method for monitoring respiration with acquired NMR data
US5200700A (en) Reduction of NMR artifacts caused by time varying linear geometric distortion
CN101283288A (zh) 扩散加权mri中的高度约束背投重构过程
US10094900B2 (en) 4D velocity profile image recording with a magnetic resonance system
JPH0616766B2 (ja) Nmrイメージング装置を使用した高速流れ測定方式
US4994744A (en) Method for combining acquired NMR data to suppress motion artifacts
US4786871A (en) NMR imaging method and apparatus
JPH06327649A (ja) 核磁気共鳴イメージング装置における画像再構成方法
WO2024081779A2 (en) Imaging slow flow dynamics of cerebrospinal fluid using magnetic resonance imaging
EP0415683A2 (en) NMR system
EP0415682A2 (en) NMR system
JPS63111847A (ja) 信号を発生する方法と装置
JPH0237770B2 (ja)

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090330

Year of fee payment: 15

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100330

Year of fee payment: 16

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100330

Year of fee payment: 16

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110330

Year of fee payment: 17

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110330

Year of fee payment: 17