JPS63119745A - 迅速nmr像形成方法及び装置 - Google Patents

迅速nmr像形成方法及び装置

Info

Publication number
JPS63119745A
JPS63119745A JP62197246A JP19724687A JPS63119745A JP S63119745 A JPS63119745 A JP S63119745A JP 62197246 A JP62197246 A JP 62197246A JP 19724687 A JP19724687 A JP 19724687A JP S63119745 A JPS63119745 A JP S63119745A
Authority
JP
Japan
Prior art keywords
waveform
applying
image
spatial frequency
region
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
JP62197246A
Other languages
English (en)
Other versions
JPH0698124B2 (ja
Inventor
アルバート マコヴスキ
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.)
Leland Stanford Junior University
Original Assignee
Leland Stanford Junior University
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 Leland Stanford Junior University filed Critical Leland Stanford Junior University
Publication of JPS63119745A publication Critical patent/JPS63119745A/ja
Publication of JPH0698124B2 publication Critical patent/JPH0698124B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related 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
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • 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
    • 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/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

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

Description

【発明の詳細な説明】 (産業上の利用分野) この発明は、NMR(核磁気共鳴)像形成方法及び装置
に関する。主な用途において本発明は、NMRデータを
迅速に取得し、緩和時間と比較的独立に高信頬の像を生
成することに関する。
(従来の技術と問題点) 既存の磁気共鳴作像系は、2DFT (2次元フーリエ
変換)つまりスピンワープ方式を用いている。この方式
は特に、W、 S、 tlinsha−及び^、H0L
entによる一般的な論文「ブロッホ方程式から作像方
程式まで;NMR像形成の概論J 、Proceedi
ngsof the IEEE、、Vol、 71 、
N13.1963年3月に記されている。この方式では
、一連の励起が使われ、各励起が空間周波数領域つまり
にスペース内のラインを与える。これは、読出前に加え
られ、kスペースの走査を所望のラインに移動させる一
傾斜軸上の可変領域位相コード化傾斜によって達成され
る。次いで続出中、kスペース内の選択ラインを走査す
るのに固定振巾の傾斜が使われる。
nxn像の場合、各々約1秒であるn個の別々の走査が
必要なので、比較的長い時間で高品質の像を与える。
より高速の走査のため、多くの手法が講じられてきた。
恐らく最も早かったのは、P、 Mansfield及
び1. L、 Pykettによる論文、Journa
l of MagneticResonances V
ol、 29 、p、 355.1978、に記された
エコー平面方式であろう。この方法では、平面の励起後
に、一定の傾斜が一軸上に、また矩形波が他軸上に加え
られる。これが1回の励起で、kスペースの迅速走査を
与える。走査の時間は妥協の結果を表わす。つまり走査
時間が比較的短いと、長い緩和時間T2を持つ物質のS
N比が妥協される。逆に比較的長いと、短いT2を持つ
物質かにスペースの部分的走査だけとなり、歪められる
迅速な像形成のための別の手法は、C,B、 Ahn。
J、 It、Kim及びZ、 Il、 Choの論文「
高速ラセン走査エコー平面式NMR像形成J 、IEE
E Trans、 onMedical Imagin
g、 Vol、M I −5,1986年3月によって
提案されている。これでは、kスペースを走査するのに
一定の周波数ラセンが使われる。
この場合にも、走査時間においてSN比と歪みの間で同
様の妥協が生じる。但し、走査は円形に対称なので、歪
みまたはぼけ関数が対称である。また、kスペースの走
査は非一様な線形速度を持ち、高い方の空間周波数で最
大となる。これは高周波数の応答に損をもたらす。
高速の磁気共鳴像形成に関する一般的手法は、1980
年8月14日提出のヨーロッパ特許出願第801048
08.3号、公告番号第0−024−640号「移動す
る傾斜ツーグマトグラフィ」 (発明者Richard
S、 Likes)に見られる。迅速なデータ取得のた
め、1回の励起後に空間周波数領域を走査する多くの方
法が示されている。これらの中には、サイン波によって
発生されるリササジュパターンとロゼツト、異なる初期
位相によるラセン、及び同心円が含まれている。こうし
た方法は、空間周波数領域を適度に満たす。しかし、比
較的短い緩和時間T2を持つ材料は、空間周波数領域の
一部だけが走査されるので歪んだ像を生じる。これらの
物質を含めようと走査を短くすると、積分時間が減少す
るため、長いT!値を持つ物質のSN比が減じる。
迅速な像形成に関する一般的な手法は、A。
Macovskiによる論文「時間変化傾斜による体積
式%式% Vo1.2.1985年2月に記されている。これに4
は、サイン波及びその他の傾斜波形を用いた方式が示し
である。また、パルス応答を補正し、時間変化傾斜をス
ペクトロスコピーに加える方法も示されている。
(問題点を解決するための手段) 体積中の一領域の磁気共鳴活性度が、空間周波数領域内
に一定速度のラセンを与える時間変化傾斜を用いて像形
成される。比較的低いTt4fLに対し高い信頼性、比
較的高いTt値に対し高いSN比、及びTt感度の良い
像を与えるため、ラセンは起点へと引き返される。傾斜
の条件は、別々の励起を用いた一連の周期すらしラセン
を使うことで最小限化される。更なる時間削減のため、
非一様性の影響を減少させる反転励起を用いて、空間周
波数領域の半分だけを走査可能である。
発明のより完全な開示のため、以下添付の図面を参照し
ながら幾つかの例示実施例について詳しく説明する。
(実施例) 発明の広い観点での理解は、対象物10のNMRパラメ
ータが像形成される第1図の系を参照することで得られ
よう。図面を見易くするため、標準的なコイル構造は全
て示してない。何故なら、それらのコイル構造はどのN
MR作像系にとっても一般的なものだからである。通例
に従い、Boは2方向で示されている。これは0.1〜
2. OTe5la間の、静的な主磁場である。それぞ
れの傾斜磁場G、 、G、及びG8は約1.0ガウス/
c11で、2方向磁場に対するxSy及び2方向の変化
分を表わす。コイル11が無線周波数の励起と受信を表
わす、一部の系では励起と受信用に別々のコイルが使わ
れ、また別の系では第1図に示すように同一のコイルが
使われる。スイッチ12が、2モ一ド間の切り換えに使
われる。スイッチ12はまず位置13に接続され、そこ
で発信器15が対象領域つまり体積10全体を励起する
のに用いられる無線周波数バーストを発生する。励起バ
ーストの後、スイッチ12は位置14に切り換えられ、
そこで励起対象物10から発せられる自由誘導減衰(F
 I D)信号を受信する。この信号がプロセッサ16
で処理され、処理後の像がディスプレイ17上に表示さ
れる。
2DFTを用いた既存の市販系では、各々が空間周波数
領域つまりにスペース内のラインを与える一連の励起が
使われる。次いで、このライン配列が所望の像を形成す
るように変換される。これらのラインは例えばG、軸等
、傾斜軸中の1つに沿った静的傾斜によって発生される
。任意の時点におけるに、値は傾斜の積分で、従ってX
方向に一様の走査を形成する。2D(次元)配列中にお
ける各ラインの位置は、続出間隔前にGyパルスの面積
によって求められる。これで各ラインのk。
位置が決まる。
この手法は高品質の像を与えるが、各励起間隔は約1秒
のオーダなので、長い作像時間を与える。
そこで本発明者等は、各励起後ににスペースの全部また
は主要部を走査することによって、極めて迅速にデータ
を収集する手法を提供する。
kスペースを迅速に走査するその他の方法も発表されて
いる。しかし、これらの方法は緩和時間T2のために歪
みを蒙っていた。像内の物質が比較的短い減衰時間T2
を有すると、その放射信号かにスペース走査の一部だけ
に存在する結果、パルス応答の劣化と像の歪みをもたら
す。こうした迅速走査を記した方式には、前述したMa
nsf 1eld。
Cho、 Likes及びMacovskiの各論文が
含まれる。
Choの論文は、起点から始まって半径方向外側に移動
するラセン波形を用いている。しかし、このラセンは一
定の角周波数なので、線形走査速度は半径の増大につれ
減少する。従って、走査はより高い空間周波数で少ない
時間を費し、それらの振巾を減じる。これはパルス応答
を広げ、解像度の劣化を引き起こす。
好ましい方式を第2図に示す。こ\では、波形が次式で
表わされる一般的な形を持つ:Gx (t) =Asi
n(BfD及びGy (t) =Acos(B−J′T
″>これは、はゾ円形の−様な分布を表わす一定速度の
ラセンを与える。半径はJ「と共に増加し、角度も−r
Vにつれて変化し、はり一定の線形速度を生じる。
第2図は、x−y平面内の断面像を発生するのに使われ
る迅速な本走査方式を示している。無線周波数励起信号
である信号13は、選択された平面の共鳴周波数に対応
したバーストである。平面の選択は、傾斜G、を用い通
常の方法で行われる。
バースト20中に生じる波形21が、共鳴を与える2位
置を表わす所望の断面を表わす。平面の励起後、スイッ
チ12が位置14に切り換えられて信号を受信し、前述
した2つの傾斜波形が加えられる。つまり、Gx(t)
がAs1n(B IT’)の波形23、C,(t)がA
cos(B−/T′)の波形34である。
緩和を無視すれば、受信信号14は次式で与えられる: 但しTは磁気回転比。この式は、空間周波数領域つまり
にスペースが走査される2次元のフーリエ変換を表わす
。こ\で、瞬間的な空間周波数kX(t) とky(t
)は次式で与えられる:起点の近い区域を除き、空間周
波数領域の半径方向成分は次式で与えられ; また角度は次式で与えられる: これは、kスペースを等しく加重する一定速度のラセン
を表わす。この点は第3図に実線27で示しである。点
線を無視すれば、kスペースの走査は起点からスタート
し、所望のラセンを生じる。
定数AとBは、連続状の円が相互の接触を避けて充分に
接近し、且つスペース全体が所望の時間内で走査される
のに適切なサンプリングを与えるように調整される。
こ\で生じる主な問題の一つは、緩和時間T2に関する
ものである。kスペースが走査されるといるサンプル中
の異なる物質は、異なった割合で減衰する。従って、長
い減衰時間を持つ一部の物質はにスペース全部に遭遇し
てパルス応答が短くなる一方、比較的短いT2値を持つ
他の物質は主に低い空間周波数に遭遇しそれに応じてパ
ルス応答が広がる。
この問題の一つの解決策は、全ての物質が小さい減衰と
なるような比較的高いB値で、極めて迅速な走査を行う
ことである。しかしこうすると、長時間信号を放射し続
けるTtの高い物質のSN比が損われる。また、得られ
た像からT2感度が取り除かれてしまう。この問題を解
消するのはにスペースの冗長走査である。最初の走査後
、経路が反転され、反対方向に走査が繰り返される。前
進及び逆転走査全体は各励起後1回以上反復でき、得ら
れた各値が加算される。これは高いT2成分に対し高い
SN比を保証すると共に、得られる像の強度値が時定数
T2により変化することを保証する。
kスペースの反転は、両傾斜波形の極性を反対にし、波
形を逆戻りすることによって達成される。
従って、2回の冗長走査について次の式が得られる: G(t) =G(t)      O< t <T=−
G (2T−t)  T<t<27=G (t−2T)
   2T<t<3T=−G (4T−t)   3T
<t<47但しTはにスペース全体を走査するのに使わ
れる時間。
第2図はそれぞれ後半区分25.26でGX(t)及び
C,(t)波形の反転と逆戻りを示している。
これらの波形は各々、反転波形全体が完了し、第3図中
実線27上の矢印で示された起点に走査を戻すまで続く
。像の再構成は、プロセッサ16を用いて達成される。
再構成の一般的手法は、前記A、 Macovskiに
よる論文「時間変化傾斜による体積式NMR像形成J 
、Magnetic Re5onance inMed
icine  Vol、2.1985年2月に記されて
いる。これでは、受信信号14が次式で与えられる傾斜
変調波形の共役と乗算される: m(xo、yo) ky(t)(31−Vo))) dxdydt但し?n
 (xo 、  )’o )はXo+  yoにおける
磁性の予測値で、系のパルス応答によるm (x、y)
のたたき込み(合成積)である。つまり、再構成像内の
各点Xo*  yo毎に、上記のように適切な関数を乗
算し、プロセッサ16で積分する。次いで、処理後の磁
性m (xo、  )’o )が17に表示される。
冗長走査の方式には、数多くの変形が可能である。応答
の平坦さを改善するため、受信動作で各種の連続的な冗
長走査を加重できる。しかし、これはSN比に悪影響を
及ぼすことがある。また、後続の励起後に反転走査を行
い、より平坦な応答を与えることもできる。
ラセン中に多くのループを含めた高解像度の冗長走査を
用いれば、 〜β〒x”+cv”−に比例する走査速度
が上昇し、傾斜増巾器からより多くのパワーを必要とす
る。また、kスペースが短時間のうちに高密度でサンプ
リングされねばならないので、必要なサンプリング速度
も上昇する。これは、より精巧なA/D変換器とより多
くのバッファ段を必要とする。つまり、迅速、冗長、高
解像度の走査はハードウェア上の負担を系に加える。
上記の負担は、第3図に示すように周期ずらしを行うこ
とによって減少できる。これはSN比も高める。一連の
励起がなされ、各励起の後ににスペースの走査が続く。
励起シーケンス全体は、人工的な動きを避けるための息
止の期間内に充分大る数秒間とし得る。心臓の動きが含
まれるか、あるいは脈動を含む血液の場合、励起はEK
G (心電図)波形にゲートできる。いずれにせよ、k
スペース内に周期ずらしセットを形成するように、各々
のシーケンスがそれぞれの初期タイミングにおいて変更
される。第3図は連続する励起によってなされる周期す
らしラセン走査を示し、第2図の波形が傾斜波形GX、
GVの初期位相を変えて繰り返されている。走査の初期
勾配α、つまり起点での勾配は次式で与えられる: kX          tJX αの値を対称的に配分して、周期ずらし走査を達成でき
る。
傾斜波形の積分の共役と乗算するという前述の再構成方
式は、計算上集中させられる。計算上の負担を最小限化
するため、周知なように、nビクセルについてn2でな
く約n log nの演算を必要とする2OFFTを利
用することもできる。高速フーリエ変換つまりFFTは
、J、 W、 Cooley及びJ、 W、 Tuke
yの論文、Mathematical Co+mput
ing+vo1.19.1965年4月、297〜30
1頁によって最初に論じられた。明解な説明は、G、 
D。
Berglandの論文、5pectrua+、 vo
l、 6.1969年7月、41〜43頁に見い出せる
。これでは、第5図の“矩形”ラセンに示しであるよう
に、2D変換が−様なデカルト系グリッド上でサンプリ
ングされる。曲線37は起点からスタートし、一定の速
度で図示の矩形ラセンを発生する。これは第4図に示す
ように、GxとGy両軸間を交互に繰り返す一連の順次
増加する長さの傾斜パルスによって一般化される。傾斜
に先行する平面励起の波形は、第2図のものと等しい。
 ゾ石x” 4Q、2は一定なので、kスペースの軌跡
は定速度で移動する。
この走査もあらゆる点で、第2及び3図の場合と同様に
処理できる。kスペースを逆戻りする反転の冗長走査で
は、時間T後に傾斜が反転され、点線の波形36で示す
ように、傾斜が逆の順序で繰り返される。また、同じく
傾斜波形の初期位相を変化することで発生される点線3
8のように、周期ずらし走査も発生可能である。1コの
時間装置で最初のパルス30と33の巾を半分にし、更
にその後の各パルスの巾も減少させれば、2対1の割合
の周期ずらし走査が得られる。それより大きい割合の周
期ずらしも同様に達成される。
処理動作では、走査が終了したところで、所望なm(x
o、yo)の2次元変換を記憶してお(。次いでその変
換値が、付随の高速及び低コストで、2DFFTを用い
てプロセッサ16により再構成される。
何、れの磁気共鳴作像系でも、励起後時間の経過した時
点の信号振巾は、不均等性及び化学的シフトによって減
じられる。異なる対象地点の時間に伴う位相シフトの増
加が、コヒーレンスを失なわせる。この位相外れは周知
のように、励起の段階に応じ、Xまたはy軸を中心とし
た180°の励起によって生じるスピンエコーを用いて
補正できる。つまり、位相外れのスピンが、位相の外れ
るのに要した時間と等しい時間内に再び位相合せつまり
集束させられる。これらのスピンエコーは第1〜5図の
系に加えられ、コヒーレンスの消失を防止して、緩和時
間の全間隔中データの収集を可能とする。
スピンエコーの一使用法を第6〜7図に示ス。
こ\で、スピンエコーはにスペースの上方半分内に矩形
ラセンを発生するのに使われる。この系も、スイッチ端
子13を介した発振器15からのバースト50によって
励起される。スライス選択用の傾斜G2は示してない。
スライス選択が所望なら、第2図の波形21と22を使
える。励起の後に一連の傾斜パルス55.53及び56
が続き、kスペースの軌跡40を起点から上方、右方更
に下方に位置41へと移動させる。この時点で、180
゜バースl−51かにスペースの軌跡を共役位置45に
移動する。この後に傾斜パルス57.54及び58が続
き、再び地点45から地点42へと上方、右方更に下方
に進む走査がなされる。またこの後に180″のバース
ト52が続き、kスペースの軌跡を42から46に移す
。このシーケンスが再度繰り返され、傾斜パルス59か
ら始って軌跡を46から43に移動させる。後続の各シ
ーケンスが、G、の上方及び下方傾斜パルスに1単位の
巾を付は加え、GXの“右方”傾斜パルスに2単位の巾
を付は加える。つまり、不図示の180°バーストを次
に続け、軌跡を43から 47、更に44へ移動した後
、別の1806バーストで44から48へと移動させる
。このようにして、牢乎面全体が発生される。
傾斜Gyを全て反転して同じシーケンスを繰り返せば、
kスペースの他の半分を与え、保全体の2DFTを行う
ことができる。次いで、これがプロセッサ16内で2D
FFTに付され、ディスプレイ17に像を表示する。か
かる180’バーストの使用は、第2及び3図の円形ラ
センにも使える。つまり、波形24のピークからピーク
まで等Gyの各完了サイクル後、180°バーストを用
いてkX軸を旋回させる。この後に波形23の正の半サ
イクルが続く。波形23の負の半サイクルは使われない
。従って、Gyの完全なサイクル、180°バースト及
びGXの半サイクルというシーケンスが得られ、第2図
に示すように周波数が減少する。
一定の条件下では、像を再構成するのに第6図に示した
にスペースの牢乎面で充分である。これは、走査時間を
半分に削減するので望ましい。所望の像m (x、y)
が実数値なら、そのフーリエ変換はエルミート行列で、
実数部分が偶数、虚数部分が奇数である。このため、正
の半分中の各M(kX、ky)について、負の半分では
M (−kX、−ky)=M“(kX、ky)とし得る
。但し*は複素共役数である。しかし、多(の磁気共鳴
系は非均質性、化学的シフト、渦巻電流等による位相シ
フトを示す。これらが非エルミート行列への変換を生じ
る。つまり、kスペースの負の牢乎面を推定するときは
、エラーが発生する。これらのエラーは、例えば一定の
位相を減じて各位相を調整し、起点がゼロの位相シフト
を持つようにすれば減少できる。
第8図は、位相シフトの影響を最小限とする別のにスペ
ース走査を示す。これでは、低周波数領域が完全な走査
を有する一方、高い周波数領域は上側の牢乎面に制限さ
れる。完全な低周波数領域は、第5図のように180°
の励起を含めないでか、あるいは第6〜7図に示すよう
に180’の励起を含め、負の半分を別に行うことによ
って生成できる。高い周波数領域は第6〜7図と同様に
、66.67.68及び69で生じそれぞれ軌跡を70
.71.72及び73へ移行させる180@の励起によ
って生成される。一般に像は、高い周波数の位相エラー
よりも低い周波数の位相エラーに対して感度が高い。k
スペース全体が存在するれ、θ(x、y)は位相シフト
のエラーである。
単純に振巾を取れば、現在の各方式で行われているよう
に歪みのないl m(x、 y) lを得られる。第8
図の方式では、逆変換に続く振巾演算でl ml(X、
 y) l  (但しmtはmの低周波数領域を表わす
)を求めることによって、まず低周波数の領域を再構成
できる。スペクトルの残りは、エルミート近似を用いた
にスペースの失なわれている半分を推定することで再構
成される。次いでこれらを組合せれば、高周波数領域の
エラーは低周波数領域よりも許容できるので許容可能な
像を生成できる。
第8図のデータ取得法を用いた別の手法では、cos 
θ(x、 y)とsinθ(x、 y)が“ゆっくり変
化する”関数で、低い周波数スペクトルを表わすこ変換
は、第8図の対称的な低周波数領域とほぼ匹敵するかま
たは小さいサイズのスペクトルL(kX、に、)である
。こう仮定すれば、対称られ、さらに変換することによ
ってL(kX、ky)を得る。θ(x、  y)とL(
kX、に、)を得れば、M(kX、に、)の残りを補正
してエルミート行列にできる。
換したM(kX、ky)本L (kx 、ky )であ
る。このM*Lは、L(ラプラス変換)であるたたき込
みでき、奇数部分の極性が反転される。
M(k、、に、)値を得れば、エルミート特性を用いて
M (−に、、−に、)=M” (k、、に? )であ
る負の半分を再構成可能である。
同じ論旨を用いたにスペースの部分走査に対する別の手
法は、周期ずらし走査を組合せることを含む。最初の走
査で、第6図のようなにスペースの上半分を生成できる
。次の励起では、第5図の方式に示すような周期ずらし
走査に続くが、kXX上下方点線で表わしであるにスペ
ースの下半分だけを生成する。こうして、kスペースの
異なる半分に1に81及び1kylの各値をずらしては
さみ込む。次いで、内挿法を用い全てのにスペースを推
定する。kスペース内の各点は、kスペースの反対側半
分に基づくエルミート反射から得られる隣接点と平均化
できる。
第3及び5図の方式は、起点からの距離が順次増加する
連続ループを含むラセン形走査の一般的方式の例に過ぎ
ない。例えば、傾斜波形を変更すれば、六角形や六角形
のパターンを発生できる。
ラセン形走査の何れにおいても、緩和時間T2による歪
みを最小限化し、T2感度を持つ像を与えるように、軌
跡の方向反転逆戻りを用いることができる。また第3及
び5図に示すように、連続走査に周期ずらし走査を加え
て、傾斜の強さとサンプリング要求を最小限化可能であ
る。さらに、どのラセン形走査においても、第6図に示
すような180”の反転を使え、また第6及び8図に示
したようなにスペースの一部の走査も使える。
以上説明した各種の高速走査実施例は、励起領域におけ
る磁気共鳴活性度の断面または投影像を与える。しかし
一部のケースでは、血管だけの投影像を与えてデータを
取得するのが望ましい。この種の方式はA、 Maco
vskiの論文、1B!!E Tran−sactio
ns on Medical I+saging、 v
ol、 Ml −1,42〜48頁(1982)に記さ
れている。高速走査は、心臓や呼吸の動きを考慮せずに
、血管像を迅速に得るのを可能とする。複数の周期ずら
し走査の実施例を使う場合には、心臓の動きによるエラ
ーを避けるように、各励起をEKG信号にタイミング合
せできる。−触に血管造影法では、非選択的励起を用い
た投影作像が使われる。つまり、第2図の励起は傾斜G
、がなければ非選択となる。
前述した各論文に記された多くの手法を用いて、静止組
織を消去し、移動する物質だけを作像するのも可能であ
る。励起は、静止物質を平衡状態に戻し、移動する物質
だけを励起するような“消去”式励起ともし得る。ある
いは、血液の異なる速度に対応して心臓サイクルの2つ
の異なる部分でデータを取得し、その後差し引くという
時間的減法も使える。
【図面の簡単な説明】 第1図は磁気共鳴作像系の概略ブロッホ図;第2図は本
発明の実施例の一組の波形図;第3図は第2図の実施例
の空間周波数図;第4図は本発明の別の実施例の一組の
傾斜波形図;第5図は第4図の実施例の空間周波数図;
第6図はにスペース平面の半分を利用した実施例の空間
周波数図;第7図は第6図の実施例に対応した一組の波
形図;及び第8図は本発明の別の実施例の空間周波数図
である。 10・・・・・・体積(対象物) 11・・・・・・励起/受信手段(コイル)15・・・
・・・発振器 16・・・・・・処理手段(プロセッサ)、20.50
〜52・・・・・・バースト23〜26;31〜36;
53〜58 :・・・・・・磁気傾斜波形 27i37・・・・・・ループ 図面のj’4”:::( 白°゛1“:二゛くニガなヒ) ぐ             −) ・ 8                   ■■ 手続補正書(方式) 特許庁長官 小 川 邦 夫 殿 1、事件の表示   昭和62年特許願第197246
号2、発明の名称   迅速NMR像形成方法及び装置
3、補正をする者 事件との関係  出願人 4、代理人

Claims (1)

  1. 【特許請求の範囲】 1、体積中の一領域の磁気共鳴像を与える方法において
    : 前記領域を無線周波数バーストで励起する工程; 磁気傾斜波形を2つの直交軸に沿って加える工程で、2
    つの傾斜の大きさがほゞ一定で、像の空間周波数を表わ
    す傾斜波形の積分が起点から順次増大する距離のループ
    配列を形成する;傾斜波形の積分を前記とほゞ同じルー
    プ配列に沿って起点に戻す磁気傾斜波形を加える工程;
    傾斜波形が加えられている間に、磁気共鳴信号を受信す
    る工程;及び 前記磁気共鳴信号を処理して前記領域の像を生成する工
    程;を含む方法。 2、体積中の一領域の磁気共鳴像を与える装置において
    : 前記領域を無線周波数バーストで励起する手段; 磁気傾斜波形を2つの直交軸に沿って加える手段で、2
    つの傾斜の大きさがほゞ一定で、像の空間周波数を表わ
    す傾斜波形の積分が起点から順次増大する距離のループ
    配列を形成する;傾斜波形の積分を前記とほゞ同じルー
    プ配列に沿って起点に戻す磁気傾斜波形を加える手段;
    傾斜波形が加えられている間に、磁気共鳴信号を受信す
    る手段;及び 前記磁気共鳴信号を処理して前記領域の像を生成する手
    段;を備えた装置。 3、前記傾斜波形を加える手段が、半径方向の空間周波
    数成分と周方向の空間周波数成分が実質上時間の平方根
    につれて変化し、空間周波数スペース内にラセン状軌跡
    が生じるように成す波形を加える手段を含む特許請求の
    範囲第2項記載の装置。 4、前記傾斜波形を加える手段が、一方の傾斜がオンで
    ある間他方の傾斜がオフになり、波形の持続時間が順次
    増大して空間周波数スペース内に“矩形”ラセンが生成
    されるように、一連の定振巾矩形波を2つの直交する傾
    斜軸に加える手段を含む特許請求の範囲第2項記載の装
    置。
JP62197246A 1986-08-06 1987-08-06 迅速nmr像形成装置 Expired - Fee Related JPH0698124B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US893793 1986-08-06
US06/893,793 US4748410A (en) 1986-08-06 1986-08-06 Rapid NMR imaging system

Publications (2)

Publication Number Publication Date
JPS63119745A true JPS63119745A (ja) 1988-05-24
JPH0698124B2 JPH0698124B2 (ja) 1994-12-07

Family

ID=25402101

Family Applications (1)

Application Number Title Priority Date Filing Date
JP62197246A Expired - Fee Related JPH0698124B2 (ja) 1986-08-06 1987-08-06 迅速nmr像形成装置

Country Status (4)

Country Link
US (1) US4748410A (ja)
EP (1) EP0256779B1 (ja)
JP (1) JPH0698124B2 (ja)
DE (1) DE3751909T2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006317435A (ja) * 2005-04-14 2006-11-24 Tokyo Electric Power Co Inc:The 位相・振幅検出装置および方法

Families Citing this family (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62139641A (ja) * 1985-12-16 1987-06-23 横河メディカルシステム株式会社 Nmrイメ−ジング装置
JP2644744B2 (ja) * 1987-02-04 1997-08-25 株式会社日立製作所 Nmrスペクトロスコピックイメージング装置
US4833407A (en) * 1987-06-24 1989-05-23 Picker International, Inc. Scan time reduction using conjugate symmetry and recalled echo
US4912413A (en) * 1987-08-14 1990-03-27 Picker International, Inc. Conjugate symmetry magnetic resonance imaging
GB2228090B (en) * 1989-02-01 1993-09-22 Gen Electric Co Plc Magnetic resonance methods
GB8915090D0 (en) * 1989-06-30 1989-08-23 Cottrell Stephen P Nuclear magnetic resonance imaging methods
US5072182A (en) * 1989-07-13 1991-12-10 Diasonics, Inc. Fourier shifting of voxel boundaries in MR chemical shift imaging spectroscopy
US5150053A (en) * 1989-07-28 1992-09-22 The Board Of Trustees Of The Leland Stanford Junior University Magnetic resonance imaging of short T2 species with improved contrast
US5025216A (en) * 1989-07-28 1991-06-18 The Board Of Trustees Of The Leland Stanford Junior University Magnetic resonance imaging of short T2 species
US5233301A (en) * 1989-07-28 1993-08-03 Board Of Trustees Of The Leland Stanford Junior University High resolution/reduced slice width magnetic resonance imaging and spectroscopy by signal combination and use of minimum phase excitation pulses
GB8918105D0 (en) * 1989-08-08 1989-09-20 Nat Res Dev Echo planar imaging using 180grad pulses
DE4035410C2 (de) * 1989-11-20 2000-03-16 Siemens Ag Pulssequenz nach dem Echoplanarverfahren
US5192909A (en) * 1991-01-07 1993-03-09 General Electric Company Spectroscopic localization using pinwheel nmr excitation pulses
US5243284A (en) * 1991-07-24 1993-09-07 Board Of Trustees Of The Leland Stanford Junior University Method of magnetic resonance reconstruction imaging from projections using partial data collected in k-space
US5270653A (en) * 1992-04-10 1993-12-14 The Board Of Trustees Of The Leland Stanford Junior University Selective three dimensional excitation for magnetic resonance applications
US5420803A (en) * 1992-11-18 1995-05-30 Technology Development Group, Inc. Enhanced resolution wafer thickness measurement system
US5349294A (en) * 1993-05-12 1994-09-20 Picker International Inc. Two and three-dimensionally selective RF pulses for magnetic resonance imaging
US5402067A (en) * 1993-08-04 1995-03-28 Board Of Trustees Of The Leland Stanford Junior University Apparatus and method for rare echo imaging using k-space spiral coverage
JP2713160B2 (ja) * 1994-03-31 1998-02-16 株式会社島津製作所 Mrイメージング装置
JP3566787B2 (ja) * 1994-08-03 2004-09-15 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Mr法
US5539313A (en) * 1995-08-04 1996-07-23 Board Of Trustees Of The Leland Stanford Junior University Full echo spiral-in/spiral-out magnetic resonance imaging
US5652516A (en) * 1996-01-22 1997-07-29 The Board Of Trustees Of The Leland Stanford Junior University Spectroscopic magnetic resonance imaging using spiral trajectories
DE19814950A1 (de) * 1998-04-03 1999-10-07 Univ Ludwigs Albert Verfahren der Kernspintomographie zur geräuschlosen Durchführung von Untersuchungen insbesondere der Hirnaktivierung
US6281681B1 (en) 1999-01-28 2001-08-28 General Electric Company Magnetic resonance imaging with interleaved Fibonacci spiral scanning
DE10016234B4 (de) 2000-03-31 2006-12-14 Siemens Ag Verfahren zur Bildgebung mittels Kernspinresonanz bei einer gekrümmt verlaufenden K-Raum-Trajektorie
US6411089B1 (en) * 2000-11-22 2002-06-25 Philips Medical Systems (Cleveland), Inc. Two-dimensional phase-conjugate symmetry reconstruction for 3d spin-warp, echo-planar and echo-volume magnetic resonance imaging
US6980846B2 (en) * 2001-03-30 2005-12-27 General Electric Company Robust coronary MR angiography without respiratory navigation
US6683454B2 (en) * 2002-03-28 2004-01-27 Ge Medical Systems Global Technology Company, Llc Shifting of artifacts by reordering of k-space
WO2004080302A1 (ja) * 2003-03-14 2004-09-23 Hitachi Medical Corporation 磁気共鳴イメージング装置
WO2004095050A1 (en) * 2003-04-23 2004-11-04 Philips Intellectual Property & Standards Gmbh Mr imaging method
DE102004019394B4 (de) * 2004-04-19 2008-04-03 Forschungszentrum Jülich GmbH Bildgebungsverfahren sowie Kernspintomograph zur Erfassung der longitudinalen Spin-Gitter Relaxationszeit
DE102004061509B4 (de) * 2004-12-21 2007-02-08 Siemens Ag Verfahren und Gerät zur beschleunigten Spiral-kodierten Bildgebung in der Magnetresonanztomographie
FR2898685A1 (fr) * 2006-03-14 2007-09-21 Centre Nat Rech Scient Traitement d'image obtenue par resonance magnetique a bandes passantes variables
DE602007003081D1 (de) * 2007-05-21 2009-12-17 Emid Alain S R Schnelles, geräuscharmes Echtzeitverfahren der bildgebenden magnetischen Resonanz mit hoher Auflösung
WO2013154834A1 (en) * 2012-04-11 2013-10-17 Dignity Health System and method for magnetic resonance imaging using three-dimensional, distributed, non-cartesian sampling trajectories
US10996299B2 (en) * 2018-08-22 2021-05-04 General Electric Company Systems and methods for waveform optimization for oblique scans

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5630636A (en) * 1979-08-20 1981-03-27 Gen Electric Conjugate recording method for movement inclination
JPS59105548A (ja) * 1982-12-09 1984-06-18 Yokogawa Hokushin Electric Corp 核磁気共鳴による検査方法及び検査装置
JPS6288949A (ja) * 1985-10-16 1987-04-23 Hitachi Ltd 核磁気共鳴を用いた磁化分布の計測方法
JPS62231642A (ja) * 1986-03-31 1987-10-12 株式会社島津製作所 Mriデ−タ収集法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
US4642567A (en) * 1984-06-04 1987-02-10 Indiana University Foundation Methods for two dimensional nuclear magnetic resonance imaging
US4707660A (en) * 1984-09-10 1987-11-17 U.S. Philips Corporation Fast method and device for determining an NMR distribution in a region of a body
NL8402959A (nl) * 1984-09-28 1986-04-16 Philips Nv Snelle werkwijze en inrichting voor het bepalen van een kernmagnetisatieverdeling in een deel van een lichaam.
US4651096A (en) * 1985-02-19 1987-03-17 Resonex Imaging method for nuclear magnetic resonance utilizing balanced gradient fields
US4727325A (en) * 1985-10-16 1988-02-23 Hitachi, Ltd. NMR imaging method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5630636A (en) * 1979-08-20 1981-03-27 Gen Electric Conjugate recording method for movement inclination
JPS59105548A (ja) * 1982-12-09 1984-06-18 Yokogawa Hokushin Electric Corp 核磁気共鳴による検査方法及び検査装置
JPS6288949A (ja) * 1985-10-16 1987-04-23 Hitachi Ltd 核磁気共鳴を用いた磁化分布の計測方法
JPS62231642A (ja) * 1986-03-31 1987-10-12 株式会社島津製作所 Mriデ−タ収集法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006317435A (ja) * 2005-04-14 2006-11-24 Tokyo Electric Power Co Inc:The 位相・振幅検出装置および方法

Also Published As

Publication number Publication date
EP0256779A3 (en) 1990-01-17
JPH0698124B2 (ja) 1994-12-07
US4748410A (en) 1988-05-31
DE3751909T2 (de) 1997-02-06
EP0256779A2 (en) 1988-02-24
DE3751909D1 (de) 1996-10-24
EP0256779B1 (en) 1996-09-18

Similar Documents

Publication Publication Date Title
JPS63119745A (ja) 迅速nmr像形成方法及び装置
JP2750880B2 (ja) 磁気共鳴画像化装置
EP0678754B1 (en) Magnetic resonance imaging methods and apparatus
JP3529446B2 (ja) Epi及びgrase mriにおける読み出し傾斜磁界極性の補正方法
US4689567A (en) NMR Fourier imaging from multiple echoes
JP4395516B2 (ja) 脂肪抑制を用いる血管の等方性画像形成
JPH0763455B2 (ja) 磁気共鳴映像装置
JPH07323021A (ja) Mrイメージング装置
JPH0256093B2 (ja)
JPH10502858A (ja) 対象の磁気共鳴画像化方法及び装置
JPH0775628A (ja) Grase磁気共鳴画像化方法及び装置
US4940941A (en) Method of high-speed magnetic resonance imaging
EP0783113B1 (en) Magnetic resonance imaging system
JP3136630B2 (ja) 核磁気共鳴装置
JP2713160B2 (ja) Mrイメージング装置
JPS61142448A (ja) Nmr像作成におけるベースライン誤差成分の影響を除く方法と装置
JPH09510898A (ja) 核スピントモグラフィ装置用のパルスシーケンス
JP2724830B2 (ja) 核磁気共鳴信号を用い、対象物から高速度で影像情報を得る装置
JPH0241967B2 (ja)
US4695800A (en) Non harmonic NMR spin echo imaging
EP0182873A1 (en) Nmr fourier imaging from multiple echoes
JPH0921853A (ja) Nmr画像におけるベースライン誤差アーティファクトを除去する方法
JPH04246327A (ja) ダイナミックmrイメージング法
US6181135B1 (en) MRI system with fractional decimation of acquired data using linear phase shift
JPH0856930A (ja) 磁気共鳴画像形成装置

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees