JPS62194845A - 不所望の応答信号を抑圧する方法と装置 - Google Patents

不所望の応答信号を抑圧する方法と装置

Info

Publication number
JPS62194845A
JPS62194845A JP62015361A JP1536187A JPS62194845A JP S62194845 A JPS62194845 A JP S62194845A JP 62015361 A JP62015361 A JP 62015361A JP 1536187 A JP1536187 A JP 1536187A JP S62194845 A JPS62194845 A JP S62194845A
Authority
JP
Japan
Prior art keywords
signal
pulse
pulse signal
response
midpoint
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
JP62015361A
Other languages
English (en)
Other versions
JPH044892B2 (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.)
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 JPS62194845A publication Critical patent/JPS62194845A/ja
Publication of JPH044892B2 publication Critical patent/JPH044892B2/ja
Granted 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/4828Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • 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/483NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • 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

Landscapes

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

Abstract

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

Description

【発明の詳細な説明】 発明の背景 この発明は核磁気共鳴(NMR)作像及び分光法、更に
具体的に云えば、検査するサンプル中に存在する不所望
の少なくとも1種類の原子核種目の共鳴による応答信号
の発生を抑圧する為に、1次以外の量子コヒーレンスを
発生する新規なパルス順序に関する。
天然に比較的多皿にある為、並びに磁気回転比が高い為
(約42. 58Ml1z /T) 、原子核の内で最
も観測し易い種目が1Hであることはよく知られている
。生物学的な材料が多量の1H原子核を含んでいて、そ
の大部分が水(H2O)分子中に入っていることもよく
知られている。従って、1Hを含むこの他のどの原子核
種目の量並びに/又は定数T I 、T 2を測定する
よりも、水分子中の1H原子核の量並びにそれに関連す
る定数T1、T2を測定する方がずっと容易である。他
の種目は、関心が持たれる区域/器官にある多くの代謝
物質の存在並びに天然に発生するその相対濃度に関する
情報が得られるのが典型的であり、こういう種目を観測
することが望ましい場合も多いが、実際には、こういう
他の原子核種目を観測することは困難であるのが普通で
ある。今日のNMR作像装置及び並びに/又は分光計は
、他の種目から放射されるNMR応答信号に同調させる
ことが出来るが、同じ生物学的なサンプル中にあるH2
0分子によって放射されるNMR応答信号の強度の為に
、他の種目のN M R応答信号は観測出来なくなる。
この為、水分子からの’HNMR応答信号の様な、不所
望の原子核種目からのNMR応答信号を選択的に除去し
て、他の種目のNMR応答信号、特に同じ種類の原子核
の他の種目のNMR応答信号を更に容易に観測出来る様
にする様な手順を提供することが非常に望ましい。
従来技術 NMRスペクトル中の水応答信号を抑圧することは、過
去数年間活発な活動がなされた1つの分野であった。こ
れまでの所、この為に利用し得る水抑圧方法の大多数は
、大別して2種類になる。
1番目の種類では、選択的な励振又は交差相関方法の内
の1つにより、程度によって差はあるが、水分子の共鳴
を励振することを避けることである。
2番目の種類では、共鳴を励振すると共に、所望の種目
からのNMR応答信号の収集が行なわれると想定される
時刻に、水分子の応答信号を除去しようとすることであ
る。この除去する試みは、普通は、(1)Hz O共鳴
の選択的な飽和、又は(2)H20共鳴を反転し、その
後不所望の共鳴のNMR応答信号が発生するのに十分な
時間(例えば0.693XT+)が経過する様にして、
データ収集の際に略振幅Oになる様にすることによって
行なわれている。高分解能の分光法に対する現存の水抑
圧方法の概観が、ジャーナル・オブ・マグネティック・
レゾナンス誌第55巻第283頁(1983年)所載の
P、  J、ホアーの論文「フーリエ変換核磁気共鳴に
於ける溶媒の抑圧」に記載されている。実際には、こう
いう方法は完全に抑圧するのではなく、進行中の実験を
阻害して、水の共鳴の領域に於けるNMR応答信号のス
ペクトルを歪める場合が多い。更に、(T +ゼロ方式
を除いて)こういう全ての方法は、水の共鳴の化学シフ
トに対して弁別するものであるから、こういう方法はN
MR像中のHz0のNMR応答信号を抑圧するのに特に
役立つものではない。
発明の要約 この発明では、互いに結合された同様な原子核を含むサ
ンプルから所望のNMR信号を求めながら、不所望のN
MR応答信号を抑圧する方法が、各々のNMR応答信号
を引出す励振信号順序の前に、無線周波(RF)パルス
を持つパルス副順序を加え、その各々のパルスの位相は
、NMRスピン系内のスピン・ポピュレーションのゼロ
量子及び多重量子コヒーレンスの内の少なくとも選ばれ
た一方を発生する様に予定の形で選ばれる工程を含む。
水の1H共鳴の様に、別のスピン共鳴と結合されていな
いスピン共鳴は、ゼロ/多重量子挙動を示さず、従って
ゼロ/多重量子側順序により、異なる強度及び位相にな
る。多重/ゼロ量子コヒーレンスは、直接的に観測する
ことは出来ないが、ゼロ/多重量子コヒーレンス側順序
の後にサンプルに印加される検出パルスに応答して、観
測することが出来る。検出パルス信号及び受信信号の位
相は、単一量子コヒーレンス信号が相殺され、且つゼロ
量子コヒーレンス又は多重量子コヒーレンスの内の選ば
れた一方を持つ信号が積極的に加算される様に操作され
る。この為、単一量子以外の挙動を持つNMR応答信号
が選択的に検出されると共に、H20から出て来る様な
単一量子信号が抑圧される。この所望の抑圧を強化する
為に、ゼロ量子/多重量子励振副順序は、(1)不所望
の、例えば水の共鳴に使われる化学選択性飽和パルス信
号と共に、又は(2)不所望の原子核種目のT1応答を
ゼロにする様に予定の形で選ばれたその後の遅延期間(
即ち、遅延時間が不所望の水の中の水素原子核に対する
0、693XT、に等しい)を伴なう反転パルス信号と
共に用いることが出来る。同様に、この他の予備飽和及
び/又は選択性励振方式をこの発明の新規な単一量子以
外のコヒーレンス発生側順序と組合せて、像及びスペク
トルの両方を収集する際のその効果を更に高めることが
出来る。
従って、この発明の目的は、単一量子コヒーレンスを持
つ原子核種目からの不所望のNMR応答信号を抑圧する
と共に、ゼロ量子又は多重童子コヒーレンスの内の一方
を持つ原子核種目からのNMR応答信号を強化する方法
を提供することである。
この発明の上記並びにその他の目的は、以下図面につい
てこの発明の詳細な説明する所から、当業者に明らかに
なろう。
不所望のNMR共鳴応答信号を抑圧するこの発明の方法
が十分に理解される様に、最初に、NMR応答信号のス
ペクトル特性を左右する幾つかの要因について幾分簡単
に説明しておく。
最初に第1a図及び第1b図について説明すると、一番
簡単な場合、分子が1個の水素原子を持ち、その原子核
の1個の陽子が1/2のスピンを持っている。このスピ
ン1/2の原子核は、他の全ての原子核から孤立してい
る。即ち、結合されていない。磁界の中に置いた時に、
この原子核が取得るエネルギ状態は2つしかない。即ち
、例えば、スピンが下向き又は(−)の状態(状態1)
であって、状態図10で第1のエネルギ・レベル10−
1として示すものになるか、或いはスピンが上向き又は
(+)の状態(状態2)で、第2のエネルギ・レベル1
0−2に示すものになるかである。レベル10−1及び
10−2の間の量子遷移11の1つしか起り得ない。原
子核をエネルギの低いレベル10−2からエネルギの高
いレベル10−1に高める為に必要な刺激エネルギEは
、この原子核がエネルギの高いレベル10−1からエネ
ルギの低いレベル10−2へ下降する時に放出されるエ
ネルギEと同じ大きさであり、一定量E−h・νである
。こ\でhはブランクの定数、νは1個のエネルギm子
の周波数π1である。この遷移の正味の磁化は±1であ
るから、単一量子遷移と呼ぶ。単一量子遷移は、NMR
像で現在直接的に観測し得る唯一の遷移である。H2O
にある2つの陽子は同一であり、その為、その2対のエ
ネルギ・レベルが整合して、許される唯一の遷移11に
より、遷移に関係する電磁エネルギの周波数スペクトル
12では、1本のスペクトル線12aになる。これらの
原子核の出カスベクトル12に周波数偏移を招く効果が
あれば、スペクトル12に他の線が存在しないので、そ
れは他のスペクトル線によってぼけることがないことが
理解されよう。
次に第2a図及び第2b図について説明すると、更に複
雑な系では、(脂質又は同様な材料の様な)分子は、1
対の同様な原子を持ち、その原子核が互いに結合されて
いる。この為、各々の原子核が他の全ての原子核から孤
立していない。説明を簡単にする為、各々の原子核が1
/2のスピンを持つと考える。(N−2)原子核の各々
が取得るエネルギ状態(E)は2つしかなく、この為、
1対の結合された原子核の組合せは、どんな時でも、取
得るNXE−2X2−4個の状態の内の1つにしかなり
得ない。両方の原子核が、状態図10′の第1のエネル
ギ・レベル10’−1で示すスピン下向き状態(状態1
)にあるか、第2のエネルギ・レベル10’−2で示す
様に、第1の原子核がスピン上向き(+)状態にあり、
第2の原子核がスピン下向き(−)状態にあるか(状態
2)、第3のエネルギ・レベル10’−3で示す様に、
第2の原子核がスピン上向き(+)状態にあり、 ゛第
1の原子核がスピン下向き(−)状態にある(状態3)
か、又は第4のエネルギ・レベル10′−4で示す様に
、両方の原子核がスピン上向き(+)状態にある(状態
4)かである。この時、合計6個の遷移が起り得る。最
初に、考えられる4個の単一遷移を考える一第1の単一
遷移11′−1(例えばレベル10’−1と10’−3
の間)は、第2の原子核だけが関係し、第1の原子核は
スピン下向き(−)状態にあり、エネルギ差Aを持つこ
とを特徴とし、これに対して第1のスペクトル線14a
が共鳴スペクトル14内で対応している。今度は第1の
原子核がスピン上向き(+)状態にある時の同じ第2の
原子核の単一の遷移11’ −2(例えばレベル10’
−2と10’−4の間)は、エネルギ差Aを特徴として
おり、これに対して第2のスペクトル線14bが対応す
る。
第1の原子核が状態を変える唯一の原子核である場合、
異なる1対の遷移が起る。即ち、第2の原子核がスピン
下向き(−)状態にある様な、第1の原子核だけの単一
の遷移11’−3(例えばレベル10’−1及び10’
−2の間)は、エネルギ差Bを特徴としており、周波数
スペクトル14で第3のスペクトル線14cがこれに対
応している。今度は第2の原子核がスピン上向き(+)
状態にある時の同じ第1の原子核の単一の遷移11’ 
−4(例えばレベル10’−3及び10′−4の間)は
、エネルギ差Cを特徴としており、これに対して第4ス
ペクトル線14dが対応する。
この為、遷移11’−1及び11’−2の二重子が同じ
(第2)のスピンから生じ、これは結合された(第1の
)スピンの分割作用の為に、中心周波数π2 (スピン
結合による分割作用がない時の、第2の原子核から得ら
れる単一の遷移による1種類の共鳴周波数)を中心とし
て変位していて、スピン−スピン結合定数Jの大きさに
よって決定される周波数の隔たりΔFAを持ち、周波数
π2に対して略対称的な1対の周波数スペクトル応答線
14a、14bを発生する。同様に、遷移11′−3及
び11’−4の二重子が同じ(第1の)スピンから生じ
、これは結合された(第2の)スピンの分割作用の為、
中心周波数π1 (スピン結合の分割作用がない状態で
第1の原子核の単一の遷移によって得られる1種類の共
鳴周波数)を中心として変位していて、やはり結合定数
Jの大きさによって決定される周波数の隔たりΔFBを
持ち、周波数π1に対して略対称的な1対の周波数スペ
クトル応答線14c、14dを発生する。これが、NM
R像で現在直接的に観測し得る唯一の遷移である。
この他に2つの遷移が起り得る。これらの遷移は、両方
の原子核のスピンが同時に反転することを必要とする。
スピンの磁化の2つの異なる状態の間の相関性がコヒー
レンスとして知られている。
1個の原子核の1回の反転によって発生される単i子コ
ヒーレンスと異なり、複数個の原子核のスピンが同時に
変化することにより、単一量子以外のコヒーレンス、即
ち、Oの正味の磁化変化(これをゼロ量子コヒーレンス
と呼ぶ)又は±M単位の正味の磁化変化(こ\でMは1
より大きい整数であり、これを多重量子コヒーレンスと
呼ぶ)を生ずる。1対の原子核の各々が、二重反転事象
より前に他方の原子核のスピンの方向と反対のスピンの
方向を持つ場合(例えば、状態10’−3で示す様に、
第1の原子核がスピン下向き(−)状態、第2の原子核
がスピン上向き(+)状態にある場合)、そして夫々が
二重反転事象の後も、他方の原子核と依然として反対の
異なるスピンの方向を持つ場合(例えば、状態10’−
2に示す様に、第1の原子核が今度はスピン上向き(+
)状態にあり、第2の原子核がスピン下向き(−)状態
にある場合)、考えられる5番目の遷移11′−5によ
り、0次量子コヒーレンスが発生される。ゼロ量子遷移
エネルギは、関係する原子核のラーモア周波数の差に対
応するのではなく、化学シフトの差に対応する(即ち、
励振信号の周波数に対する応答信号のピーク・オフセッ
ト周波数)。ゼロ量子遷移は、線拡がり効果(磁界の非
均質性等による)が相殺される為に、特に有用である。
このため、ゼロ量子遷移の線幅は、T2fではなく、1
2時定数によって決定される。1対の原子核の各々が、
二重反転事象の前の、他方の原子核のスピン方向と同一
のスピン方向を持つ場合(例えば、状態10’−1に示
すように、第1及び第2の原子核が共にスピン下向き状
態にある場合)、そして各々が二重反転事象の後にも他
方の原子核と依然として同じ別のスピン方向を持つ場合
(例えば状態10’−4に示すように、第1及び第2の
原子核が共にスピン上向き状態にある場合)、考えられ
る6番目の遷移11’−6により、2次量子コヒーレン
スが発生される。二重量子遷移エネルギは、関係する1
対の原子核のラーモア周波数の和ではなく、正確な励振
送信周波数に対する関連する化学シフトの和に対応する
。この為、こういう0次及び多重次量子コヒーレンスは
比較的低い周波数の応答信号を発生し、これは(スピン
系の横平面内の磁化を発生する単一量子遷移とは異なり
)、現在では直接的に観/l1llすることが出来ない
。N個の原子核を持つ系は最大N次のコヒーレンスを持
ち得る。原子核2個の系では、2次より高次のコヒーレ
ンスは存在し得ない。
次に第3a図及び第3b図について説明すると、(単純
な有機物質、例えば乳酸の様な)更に複雑な分子は、各
々の分子内にかなり多数の原子を持つことがある。例と
して云うに過ぎないが、結合された4個の原子を持ち、
AB3形、即ち1個の原子Aと3つの同質原子核原子B
を持つ分子に対する単純化したスピン密度マトリクスを
考える。
磁化状態  :S:T:U:V:W  :X  :Y:
Z:S:Ml=+372.M2−+l/2:−−:le
: 2: 3:  la : 2  : 3: 4:T
二重1−+172.M2−+l/2:    :−−:
le:   2:   0   :   lb  : 
  2:   3:U:Ml−−172,M2−+l/
2:  :  :−−:le:  l” : O:1c
: 2:V:Ml−3/2.M2−+1/2:    
:    :    ニー:   2    :   
1”   :   旧 1d:W:Ml−+372.M
2−1/2:  :  :  :  ニー:  if 
: 2: 3:X:Ml−+1/2.M2−1/2: 
 :  :  :  :   ニー:l(’:  2:
Y:Ml−1/2.M2−1/2:  :  :  :
  :   :   ニー:1t:Z:Ml−3/2.
M2−1/2:  :  :  :  、   、  
 、  、  。
こ\でスピン密度状態Oがゼロ量子コヒーレンスを表わ
し、状態1a乃至1fは、相異なる6つの単純な単一量
子コヒーレンスを表わし、状態1ゞは組合せ単一量子コ
ヒーレンスを表わし、状態2゜3及び4は夫々二重、三
重及び四重量子コヒーレンス表わす。第3a図に別の形
で示されている様に、4つの原子核が、エネルギ・レベ
ル線図1O″に合計NXE=4X2=8個のエネルギ・
レベル10’−1乃至10’−8を発生するが、合計2
8個の遷移を持ち得る。こういう遷移全体の中には、原
子・核Aに対する単純な単一量子遷移の四重子(即ち、
遷移11’−1乃至11’−4)、各々の原子核Bに対
する単−皿子遷移11’−5乃至11’−7の3個1組
(原子核Aの第1のスピン状態に対し)、原子核Aの他
方のスピン状態に対する、各々の原子核Bに対する単一
量子遷移11’−8乃至11’−10の3個1組、原子
核Aと3つの原子核Bの内の相異なる1つとの同時の反
対向きの反転に関係するゼロ量子コヒーレンスを発生す
る遷移11’−11乃至11’−13の3個1組、何れ
も正味の単一量子遷移になる様な複数個のスピン反転を
持つ1対の単一量子の組合せ15−1.15−2)(原
子核Aの第1のスピン状態に対して起る)第1の対16
−1a、16−1b及び(原子核Aの他方のスピン状態
に対して起る)第2の対16−1c、16−1dを含む
2対の二重量子コヒーレンスを発生する遷移、4つの組
合せ二重量子遷移16−1乃至16−1as4つの三重
量子コヒーレンスを発生する遷移16−2a乃至16−
2b、及び(関係する原子核の数N−4に等しい最高次
である)四重量子コヒーレンスを発生する遷移16−3
が含まれる。
通常の単一量子コヒーレンスを発生する遷移だけが、直
接的に求めることが出来、且っNMR作像及び/又は分
光法に利用することの出来る応答スペクトル信号部分を
発生する。この為、4つの単一遷移11’−1乃至11
’−4が、原子核Aの非結合時の共鳴周波数π3を中心
とし、結合定数Jの倍数の周波数の隔たり(典型的には
陽子では1.000Hz未満)を持つ共鳴スペクトル線
18a乃至18dの4個1組を発生する。原子核Bの夫
々3個1組の遷移11’−5乃至11’−7又は11’
−8乃至11’−10は、原子核Bの共鳴周波数π4を
中心とし、異なる原子核の間のスピン−スピン結合定数
Jに関係する周波数の隔たりを持つ二重共鳴スペクトル
線19a、19b(これらは実際には三重縮退共鳴であ
る)の1つを設定する。種々の原子核の間の空間的な結
合関係が、単一量子以外のコヒーレンスを設定するエネ
ルギ差の大きさを決定するから、残りのコヒーレンスの
内の少なくとも1つから応答データを求め、検査される
生物学的なサンプル中に大量の水等の物質があることに
よる例外的に強い単一量子コヒーレンスによる前に述べ
た問題を避けることが望ましい。
第4図には、(公知の様に、デカルト座標系を持つNM
R作像/分光装置の夫々xSy、z軸に沿って)3個1
組の互いに直交する磁界勾配G、X。
Gy、G:!を発生する様な核磁気共鳴装置に対する、
単一量子以外のコヒーレンス応答を求める為の現在好ま
しいと考えられる1組の信号波形が示されている。更に
、限定されたサンプル・スライス内にある物質の原子核
を励振する為の無線周波(RF)信号の振幅特性が磁界
勾配と共に示されている。これらの信号は、単一量子以
外のコヒーレンスを刺激して応答信号を発生するだけで
なく、その応答信号をNMR装置のRF分光計受信機で
受信する為の準単−量子形式に変換する(真の単一量子
信号自体が相殺される形に変換されるので)。この発明
の方法は、従来公知の大きなグループを成す応答信号を
引出すNMR順序22(例えば図示の様な標準形のスピ
ン捩れ形励振順序)のどれより前に、ゼロ量子/多重量
子信号副順序20を用いる。
この発明の1つの考えでは、副順序20は、単一量子以
外の遷移周波数によって決定される速度で、スピンの磁
化の位相を正弦状に変化させるのに必要な状態を作る。
その後、所望の単一量子以外のコヒーレンスの瞬時位相
が検出パルスによって変換される。図示の実施例の副順
序は、時刻toに第1のRF信号25の前縁25aから
始まる。この信号はπ/2又は90″の信号パルスであ
る。サンプルがパルス状RF信号25によって照射され
、3つの磁界勾配は全て振幅が略ゼロである。RFパル
ス信号は略四角の包絡線を持つと共に、基準位相に対し
て、複数個の位相の内の1つに選択的に設定された位相
φaを持っている。
RF倍信号位相を選択することが出来る様な1つの装置
が、1985年9月23日に出願された係属中の米国特
許出願通し番号第779,338号に記載されている。
パルス信号25は時刻t1に後縁25bで略ゼロの振幅
に戻る。2番目のRF信号パルス27がそれより後の時
刻t4に前縁27aを持ち、この信号パルス27はπ又
は180°信号パルスである。このスピン反転RFパル
ス信号も略四角の包路線を持ち、位相φ5を持ってε)
る。位相φbは、利用し得る複数個の位相の内の独立に
選択された1つである。2番目のRF信号パルスは時刻
t5に後縁27bで終る。2番目のRFパルス信号27
のタイミングは、その時間的な中点27cが、1番目の
RF信号パルスの時間的な中点25cから大体始まって
持続時間Tを持つ第1の期間の終りに来る様になってい
る。
この副順序の期間の持続時間Tはスペクトル線の周波数
の差Jの逆数に比例する様に、即ちT−N/(4J)、
(N−1,3,5・・・)に設定される。
N−1では、例として結合定数Jの範囲が約1゜乃至1
.000Hzとすると、期間Tは約0.25乃至25ミ
リ秒の範囲になる。期間Tの実際の長さは特に重要では
ない。これは、前に述べた様に、(同等スピン緩和がな
い場合)多重量子コヒーレンスの大きさが、Tの関数と
して正弦状に変化するからである。(例えば、単一量子
以外のコヒーレンスの大きさMQC(t)がMQC(t
)−ACOS (ωt)で表わされる場合(tが単一量
子以外のコヒーレンスが設定されてからの経過時間、ω
が共鳴オフセット周波数(ラーモア周波数以外の)、多
重量子遷移では和、又はゼロ量子遷移では差である)L
を正の整数(0を含む)として、T−(2L+1)/4
 Jである時に最大になる。)3M目のRF信号29の
前縁29aが時刻t8に発生する。この信号はRF信号
パルスであって、パルス25の位相に対しである余分の
移相θを持つことが出来るが、そうする様にする必要は
なく、その持続時間はπ/4の倍数Kに選び、K・45
6の原子核回転効果を持たせる。このRF信号パルスは
、略矩形の包路線を持っていて、それより後の時刻t9
に後縁29bを有する。図示の実施例の副順序では、K
−2であり、パルス信号29は90’RFパルスである
。この3番目のRF信号パルスの時間的な中点29cは
、2番目のパルスの中心27cから、パルスの時間的な
中心25c及び27cの間の期間Tと同じ期間Tの大体
終りに発生する様に選ばれる。単一量子以外のコヒーレ
ンスは、共鳴オフセットの程度に対する依存性がある為
に、効率は低下するが、T−1/4Jだけ離れた90″
パルス25及び29だけで構成された副順序によっても
励振することが出来る。
π信号パルス27を追加することにより、この依存性が
なくなる。この為、パルス25.27.29は単一量子
以外のコヒーレンスを励振するのに役立つが、この後の
スピン捩れ形順序22は単一量子コヒーレンス形の信号
にだけ作用し得る。単一量子以外のコヒーレンス形の信
号を単一量子コヒーレンス形の信号に変換する方法がな
ければならない。
このπ/2−T−π−T−π/2パルス順序にスピン・
ベクトル・モデルを適用すると、スピンエコ一応答信号
が発生され、全てのスピンの磁化は2番目のπ/2RF
パルス信号によって反転状態になることが予測される。
この結果は、H2O中・1H原子核等の様に、結合され
ていないスピンについてだけ成立する。結合されたスピ
ンでは、挙動が異なる。例えば、この列順序のみに応答
して、エコーが形成されない。更に、普通のスピン舎ベ
クトル・モデルは、スピンが結合された原子核の挙動を
適切に予測するものではない。もとのスピン密度マトリ
クスを適当な回転演算子及びスカラー乗算の作用にかけ
、第2a図及び第2b図について説明したのと同様な、
比較的単純なABスピン1/2系に対して第4a図乃至
第4f図に示す様な種々のスピン密度マトリクスσ[t
lを求めることにより、時間的な解析を行なわなければ
ならない。スピン密度マトリクスσ[1]は方形マトリ
クスであり、1辺に(m+1)X (’n+1)個の要
素を持っている。ニーでm及びnはAmBo系の次数の
添字であり、各々のマトリクス要素は系の初期エネルギ
状態及び最終的なエネルギ状態に考えられる1対の間の
相関性を持っている。時刻t。の列順序20の始めに、
第1のスピン密度マトリクスσ[tolに示す様に、核
スピンは平衡状態にあり、コヒーレンスは存在しない。
スピン密度マトリクスの対角線から外れた全ての要素は
ゼロであり、ゼロでない対角線上の要素は、位置[1,
11及び[4,4]にあるものだけであり、反対向きの
スピン、例えば水平方向のスピンを持っている。例とし
て位相Ooの第1のπ/2パルス25に応答して、マト
リクスσ[to ]に適当な回転演算子を適用すると、
マトリクスσ[t1]になり、これは単一量子マトリク
ス要素の全部、即ち、要素[1,2]、[1゜3] 、
 [2,11、[2,4コ 、 [3,11、[3,4
]、[4,2]及び[4,3]にコヒーレンスが発生さ
れたことを示す。考えられる二重量子要素[1,4]及
び[4,11は依然として作動されないでいる。パルス
に応答して、全てのコヒーレンスはY軸、即ち垂直方向
に沿った磁化に対応していなければならない。この為、
マトリクスの対角線より右側のスピンは上向きであり、
対角線より左側にある像となるスピンは全部下向きであ
る。応答信号は第4g図に示す1対の二重子で構成され
る(信号34a、34bの第1の二重子Aは単一量子コ
ヒーレンスa、bに夫々対応し、信号34c及び34d
の第2の二重子Bは夫々単一量子コヒーレンスc、dに
対応する)。期間T−1/4Jが経過した後、時刻t4
のスピン密度マトリクスσ[t4]が得られる。種々の
コヒーレンスの夫々の位相が変化しており、その絶対的
な位相はその周波数オフセットの関数である。
共鳴の間の間隔Jの為、各々の二重子A及びB(第4g
図)の成分の間の相対的な位相差は900のま\である
。igo’パルス信号27がスピンの位相を第4d図の
スピン・マトリクスσ[t5]に示す位相に回転させ、
この後の期間Tに於けるその発展により、コヒーレンス
が第4e図のマトリクスσ[t8]のX軸スピン磁化に
なる。然し、二重子の1つのピークに対応する各々のコ
ヒーレンスは、他方の二重子のコヒーレンスの極性(こ
れはスピン中ベクトル・モデルの反位相磁化と同様なも
のと考えることが出来る)に対して反対の極性を持って
いる。2番目の90°パルス信号31が回転演算子とな
り、その結果第4f図のスピン密度マトリクスσ[t9
]になる。
この時、各々の単一量子コヒーレンスの二重子の1対の
成分(例えば、スピン要素36aと36C1又はスピン
要素36bと36d)は、反対向きの位相の回転の為、
互いに実質的に相殺する様に操作することが出来る。然
し、今度は二重量子コヒーレンス38が発生され、π信
号パルス29が存在する為、その振幅は共鳴オフセット
の量に略無関係である。
この発明の別の考えとして、3番目のこの列順序のパル
ス29が中断した後に存在する単一量子以外のコヒーレ
ンス信号が、検査されるサンプルに検出RF信号パルス
を印加することによって、単一量子コヒーレンス状の検
出可能な信号に変換される。このパルスは、単−量子コ
ヒーレンスによって発生された信号の位相を互いに反対
に設定しながら、単一量子以外の信号を受信平面(例え
ばX−Y横平面)内に整合させると共に、互いに同相に
整合させる様に作用する。サンプルの一連のこの副順序
/順序の照会により、この副順序のパルス信号の位相が
予定の形(pをゼロより大きい整数として、典型的には
2p個の繰返しを持つ)で操作することにより、特定次
数の単一量子以外の信号を積極的に相加わる様にすると
共に、単一量子信号が互いに相殺する様にすることが出
来る。
一般的に、単−量子コヒーレンスによって発生される信
号の相殺の程度は、位相の組合せの数が増加した位相プ
ログラムで増加するが、繰返しプログラムを増加するに
は、完了までの時間を長くすることが必要である。この
様な多重工程位相プログラムに役立つ1つの検出パルス
31が図示のパルスであり、これは略矩形の包絡線を持
つことが出来、前縁31aから始まり、その後、パルス
3Iヲ180’又はπパルスにするのに十分な時間後に
後縁31bを持つことが出来る。検出パルスの時間的な
中心点31cは、この副順序の3番目のパルスの時間的
な中心29cから選ばれた検出期間Td後に発生する。
検出パルス31は標準的なスピン捩れ形作像順序の始め
イこ使われるスライス選択性RFパルス信号であってよ
いが、更に一般的な実験には、検出パルス31が終了し
た後の開始時刻t6に始まるスピン捩れ形作像順序22
を使うことが出来る。
多数の振幅の内の所要の1つを持つY軸磁界勾配Gy信
号40と所要の振幅を持つX軸磁界勾配G1信号42が
、作像しようとするサンプルの位相符号化された「縞」
を限定する。この後の180°RFパルス信号44がZ
軸磁界勾配G2信号部分のパルス46aと共に現れ、作
像しようとするサンプルの横方向スライスを選択する。
5inCで変調されたRFπパルス、スライスの上下の
縁を比較的明確に限定す゛るのに役立つ。公知の形で、
スピン捩れ形順序には02位相戻し部分のパルス46b
を使うことが出来る。この信号は、検出パルス31によ
って単一量子状の信号に変換された単一量子以外のコヒ
ーレンス信号が、位相外しを含めた同じ効果を全て持つ
様にすることが出来る点で、非常に役に立つことがある
。パルス31とこの後のスライス選択パルス44の間の
タイミングは、(n−1)Tであるこの期間が、整数個
のサイクル、即ちqを整数として、nT=q+(1/2
J)サイクルより大きな(1/2 J)サイクルの全磁
化回転を持つのに十分になる様にし、応答信号を収集す
る前に、各々の周波数スペクトルの多重子に適正な位相
が復元するようにする。
この為、交番位相であって、その為に反対向きを向いて
いた多重子のスピンが、同じ位相を達成するのに十分な
遅延1/2Jを経由する。従って、この結果前られるス
ピンエコ一応答信号50は、選択パルスの中心時刻t1
より同じ期間nT後に発生する時刻tyを中心とする。
NMR分光計の受信機は、応答期間の中心時刻tyを中
心とする時間的な窓52の間だけ、サンプル応答信号5
0を受入れる様にゲートするのが有利である。受信した
応答信号の位相を、前に引用した係属中の米国特許出願
に記載されている様な装置を用いて、単一量子応答信号
を相殺して、単−量子以外の応答信号を受信する為にそ
の特使われている多重工程位相サイクル・プログラムの
内の1つによって決定された量だけ、受信機内で移相す
る。単一量子以外の副順序20及び作像順序22の2p
個の位相サイクルの繰返しの内の現在の1つの初めと、
こういう繰返しの次のものとの間の繰返し期間T「があ
る。
2p工程の位相サイクル・プログラムの一番簡単な場合
が、p−1とした2工程プログラムであり、1例は次の
通りである。
工程の番号:φ25:φ27:φ29:φ31::φr
:1   :   0:+90:   Q:  −::
  Q:2       :+90:/ど180:  
 +90:    −: ・ 180 :こ\で全ての
角度φは度数であり、φ、は受信した応答信号に対して
受信機で導入した位相変化である。検出パルスの位相(
φ3、)はこの順序で移相する必要がないことに注意さ
れたい。更に高度に単一量子コヒーレンス信号を相殺す
る数多くのプログラムが考えられる。役に立つI)−2
の1つの4工程位相サイクル・プログラム(全ての角度
φはやはり度数)は次の通りである。
工程の番号:φ25:φ27:φ29:φ31::φ「
:1   :   o:+qo:   O:  −::
  Q:2      :  +90:/180:  
 +90:    −::  180:3   :  
180:  −9o:  180:  −::  Q:
4   :  −90:   Q:  −9o:  −
:: 180:この他にも多くの位相サイクル・プログ
ラムがある。次数のパラメータpが増加するにつれて、
その数が増加する。この様なプログラムは多数あるが、
少なくとも1つの位相に1組の初期並びに/又は最終条
件を課すことにより、p−tとした、2p“4−16エ
程の1例の位相サイクル・プログラムが得られる。この
様な1つのプログラムは、全ての角度φを度数を表わし
て、次の通りである。
工程の番号:φ25:φ27:φ29:φ31::φr
:1   :   O:+90:   O:   0:
:  O:2   :   0:  +90:   0
:  +90:: −9o:3      :    
0:   +90:     0:   180:: 
 180:4   :   0;  +90:   0
:  −90:: +90:5二+90:180:+9
0:0::Q:6   :  +9o:  180: 
 +90:  +90:: −90ニア   :  +
90:  180:  +90:  180:: 18
0:8   :  +90:  180:  +90:
  −90:: +90:9   :  180:  
−9o:  180:   0::  0:10   
:  180:  −90:  180:  +90:
: −90:11:180ニー90:180:180:
二180:12   :  180:  −90:  
180:  −90:: +90:13   :  −
9o:   Q:  −9o:   O::  Q:1
 4      :  −9o:    Q:  −9
o:  +90::  +90:15   :  −9
o:   0:  −90:  180:: 180:
16   :  −90:   Q:  −90:  
−90:: +90:こういう各々の位相サイクル・プ
ログラムで、後の励振によって得られた応答信号を(受
信して処理した後)前の励振から得られた応答信号(既
に受信して処理しである)から減算する。p−1の2工
程プログラムでは、2番目の応答を収集し、最初に収集
された応答から減算する。単一量子コヒーレンスには移
相は何の影響もないから、それらの同様な応答が減算に
よって実質的に相殺される。偶数次の多重量子コヒーレ
ンスが励振順序の90″の移相によって反転されるから
、多重量子コヒーレンスを経由したスピンからの応答信
号は、前の応答信号及び後の応答信号を減算する時に、
組合さる。位相が異なる2つより多くの副順序を持つ位
相サイクル・プログラムを使うことが、パルスの不完全
さの影響を最小限に抑える為に通常要求される。この発
明では、p−6まで(即ち64工程まで)の上に述べた
位相サイクル・プログラムを及び同様な位相サイクル・
プログラムを使って、多重量子信号フィルタを求めた。
p−5(即ち32工程)フィルタを用いて、隣接した水
をベースとする1H応答信号の典型的な排除率が約2,
500倍、時+、:ハ13. 000倍と云つ様に高く
することが出来た。p−9又はそれ以上にした更に長い
プログラムも、必要に応じて、計算して使うことが出来
ることを承知されたい。
NMR装置の中心軸線(今の場合はZ軸)に沿った磁界
勾配(今の場合は勾配G:りに、信号パルス56.58
で表わされる様な1対のパルス信号を追加して、引出そ
うとするコヒーレンスの次数Nを選択することが出来る
。1番目の磁界勾配パルス56は、時間的には、引出し
期間の終りに、但し検出パルス31の前に発生する様に
定める。
この為、パルス56は時間的にはRF信号パルス29及
び31の間にある。パルス56の前縁は、順序の1番目
のRF信号パルスの後縁29bが発生する時刻t9より
後の時刻t。に発生する。1番目のG:!パルスの後縁
56bは後の時刻tdに終了する。時刻1dは、検出パ
ルスの前縁31aが発生する時刻teより前である。2
番目の磁界勾配パルス58は、時間的には、検出パルス
31より後、但し作像順序22の初めより前に発生する
様に定める。パルス58は前縁58aがら始まるが、こ
れは、検出パルスの後縁31bが発生ずる時刻tfより
後の時刻t9に発生する。2番目の02パルスの後縁5
8bはそれより後の時刻thに終了するが、これは主作
像順序22が開始する時刻tsより前である。2番目の
パルス58の振幅と持続時間の積は、1番目のパルス5
6の振幅と持続時間の積のN倍である。N次量子コヒー
レンスを持つスピンだけが検出される。これは他の次数
の全てのスピンが位相外れになるがらである。
理論的には不必要であるが、現実のNMR作像装置では
、中心軸の磁界勾配(例えばZ軸に沿った勾配G:りに
、RF信号πパルス27に対して対称的に配置して、信
号パルス60.62によって表わされる様な釣合った1
対のパルス信号を追加するのが役立つことが判った。こ
ういうパルスを用いて、横方向(例えばX−Y)平面内
の擬似的な磁化等の影響を除くことが出来る。こういう
影響は、再集束パルス27が完全に矩形でないこと等に
よって起ることがある。即ち、1番目の07Lパルス6
0は、1番目のπ/2パルス25が終る時刻t1より後
の時刻t2に前縁60aを持ち、πパルス27が開始す
る前の時刻t3に後縁60bを持っている。この為、パ
ルス60の持続時間は(t3   tz)である。これ
は2番目のパルス62と略同じ持続時間である。2番目
のパルスの前縁は、πパルス27が終了する時刻t5よ
り後の時刻t6であり、後縁62bは、K・π/4パル
ス29が始まる前の時刻t7である。両方ノハルス60
.62は略同じ振幅を持っている。
単一量子以外のコヒーレンスの変調により、そのコヒー
レンスからの信号が減衰したり、或いは失われることが
ある。減衰の程度は、前に述べた様に、パルス25,2
7.29の時間関係の関数として変化する。この順序の
単一量子以外のコヒーレンスを発生する部分の終り(例
えばK・π/4パルスの終り)とその後の作像順序の間
の期間内に、単一量子以外のエコー信号を形成する為の
再集束効果を用いることが出来る。検出パルスより前に
再集束パルスを置くことによって、多重量子エコー中の
時刻に検出パルスを印加することが出来る様にし、こう
して最大限のコヒーレンスを検出することが出来る様に
するのがを利である。
この発明では、理論的には必要ではないが、或いは上に
述べた方法を実際に利用するのに不可欠ではないが、引
出し期間Td中の再集束用のπRF信号パルス64が単
一量子以外の遷移を位相戻しして、時刻teに於けるそ
の位相が時刻t9に於ける位相と略同じである様にし、
こうして励振共鳴周波数のオフセットに依存する位相誤
差を除くことが出来ることが判った。即ち、パルス64
の前縁64aはこの副順序の2番目のπ/2パルス29
が終る時刻t9より後の時刻taにあり、後R64bは
検出パルスの前縁31aより前の時刻t、にある。再集
束パルスの時間的な中点64cと2番目のパルスの時間
的な中点29cの間の期間Tcは、期間Tdの半分にす
るのが理想的である。今説明している実施例では、信号
パルス64は最初の位相戻しパルス27の位相と略同じ
位相を持っているが、考えられるある位相サイクル・プ
ログラムでは、パルス信号64の位相を段階的に変える
ことが出来る。
次に第5a図及び第5b図について説明すると、P+ 
 (π/2+α)−T−π−T−P2  (π/2+β
)の副順序によって発生された二重量子コヒーレンス(
第5a図)及び単一量子コヒーレンス(第5b図)の正
規化された量が、パルスP1(縦軸70に示す)のはじ
き角度誤差(α)及びパルスP2  (横軸71に示す
)のはじき角度誤差(β)の回転効果の関数として、グ
ラフに示されている。こ\でT−1/4Jである。即ち
、二重量子コヒーレンスの二は、P+ −90’及びP
2−90° (即ち、α−β−θ″の点72)で最大で
あり(即ち1.0の正規化された振幅を持つ)が、振幅
は中位のはじき角度誤差にあまり影響されない。−谷内
側の輪郭線73a上又はその中でのはじき角度P1及び
P2のあらゆる組合せに対し、依然として5/6 (0
,833・・・)未満の二重量子コヒーレンスの振幅が
得られる。他の輪郭線73d乃至73eに関連した二重
量子コヒーレンスの振幅を記入しである。振幅が0. 
5に減少するまでには、比較的大きなはじき角度誤差が
必要であることが判る。第5b図では、単一量子コヒー
レンス応答信号の振幅が、縦軸74に示したはじき角度
誤差α及び横軸75に示したはじき角度誤差βの組合せ
に対してグラフとして示されているが、単一量子コヒー
レンス信号の最小(略ゼロの)振幅が、中心点76a 
(P+ −π/2−90″、P2−π/2−90”)で
得られ、この時α−β−0@であり、二重量子コヒーレ
ンス信号が最大であることが判る。最小の単一量子コヒ
ーレンス応答は他の点76b乃至76eでも得られるが
、これらの4つの点は何れも、はじき角度の両方ではな
いとしても、少なくとも一方がO″又は180’である
ことを必要とし、こういう角度では、二重量子コヒーレ
ンス応答信号の関連する振幅が実質的にゼロである。単
一量子コヒーレンス応答信号の振幅輪郭線78a乃至7
8dは、こういう応答が、一方のはじき角度が90’で
あり、他方のはじき角度が00又は180°である点7
9a乃至79dでだけ、0.5の正規化された最大値に
達することを示している。この副順序によっては、目に
つく様な二のゼロ量子コヒーレンス応答信号が発生され
ない。
第6a図及び第6b図には、P2をパルスP1の位相に
対して90°移相したP+  (π/2+α)−丁−π
−TP2(π/2+β)の副順序によって発生されるゼ
ロ量子及び二重量子コヒーレンス(第6a図)及び単一
量子コヒーレンス(第6b図)の正規化した量が、パル
スP1のはじき角度誤差(α)(縦軸80に示す)及び
パルスP2のはじき角度誤差(β)(横軸81に示す)
の回転効果の関数としてグラフに示しである。点82a
−1(又は82a−2)に於けるPi−yr/2−90
01及びに−1として、P2−π/4−45° (又は
に−3として、P2−3π/4−135@)の所望のは
じき角度により、ゼロ及び二重量子コヒーレンス応答信
号の正規化された最大値が0.5になる。第5a図及び
第5b図のπ/2−T−π−T−π/2副順序(この場
合、2番目のπ/2信号の位相は2番目のπ/2信号パ
ルスの位相と同じであった)と異なり、このπ/2−T
−π−T−π/4副順序により、略同量のゼロ量子及び
二重量子コヒーレンスが発生される。この為、ゼロ量子
コヒーレンス応答の信号対雑音比は、対応する二重量子
コヒーレンス応答の信号対雑音比の半分である。他のあ
らゆる種類のコヒーレンスは拡がるが、ゼロ量子遷移は
均質でない磁界によって拡がることがないから、ホモス
ポイル・パルス56(第4図参照)を(例えばに−1と
して、π/4パルス29の後に)使って、ゼロ量子コヒ
ーレンス応答信号の大きさに実質的に影響せずに、望ま
しくない単−量子及び多重量子コヒーレンス応答信号を
破壊することが出来る。線82bに沿って発生される単
一量子以外のコヒーレンス応答信号の大きさが実質的に
ゼロになり、P1パルスのはじき角が略π/2−90”
に等しいことに注意されたい。輪郭線84a乃至84e
の大きな卵形パターンは、PIのはじき角度の比較的大
きな誤差α及びはじき角度P2のそれより幾分小さい誤
差βがあっても、単一量子以外の所望のコヒーレンス応
答信号の振幅に目立った減少がないことを示している。
第6b図では、はじき角度誤差αを縦軸85、はじき角
度誤差βを横軸°86に示しであるが、この図に示す様
に、この副順序によって発生される単一量子コヒーレン
ス応答信号は予想点87a−1及び87a−2で(P+
−π/2及びP2−π/4又は3π/4で、単一量子以
外のコヒーレンス応答信号の振幅が最大になる点で)、
並びにゼロ振幅線87aの全体に沿って、及び点87a
−2及び87a−3で、振幅が略ゼロである。単一量子
コヒーレンス応答信号は点88aで振幅0.707・・
・の正規化された絶対最大値を持ち、1対の点88b、
88cで振幅0. 5を持つ相対的な極大値を持ってい
る。
第7a図及び第7b図を見れば、主たる静磁界が1.5
テスラである研究用NMR装置で、前に引用した係属中
の米国特許出願の位相操作装置を用いて収集された像に
対する、これまで説明した方法の実用的な利点が理解さ
れよう。この順序の2p個の相異なる位相の組合せの内
の1つで励振する為に、同相又は反対位相である多重量
子コヒーレンスが発生される。各々の同相応答で収集さ
れたデータを第1のデータ・バッファに加算し、これに
対して各々の反対位相応答で収集されたデータを第2の
バッファに加算する。2つのバッファの内容を加算すれ
ば、単一量子コヒーレンス信号が残るが、多重量子コヒ
ーレンスの位相特性を持つ信号が相殺される。2つのバ
ッファの内容の減算によって、単一量子応答信号が抑圧
され、多重量子成分が残る。
1対のファントムを作った。第1のファントムは着色し
た(CuSO4をドープした)水を入れたバイアルであ
り、第2のファントムは、脂質擬似体(機械用軽油)を
大体半分溝たしたポリエチレンびんである。1対のMQ
Cパルス25.29、再集束パルス27及び検出パルス
31を持つ多重量子フィルタ副順序を先行させたスピン
捩れ形パルス順序を用いて、これらのファントムの作像
をした。10Hzの結合定数Jを仮定したが、この結果
パルス間遅延Tは約25ミリ秒になる。両方のスピン捩
れ形作像順序で約35ミリ秒の遅延期間Tdを用いた。
順序を繰返し間隔T2)は600ミリ°秒にした。一番
簡単な2工程(p−1)の位相サイクルを用いた。
第7a図の写真に示す加算像は、単一量子コヒーレンス
を示すスピンから得られた情報を示している。像91は
、水のファントム中の結合されていないスピンによるも
のと思われる。油は結合されたスピン及び結合されてい
ないスピンの両方を持つ原子核で構成されているから、
この写真には油像92が現れている。油サンプル中にあ
る結合されていないスピンの数は比較的少ないから、氷
像がずっと強いことが判る。第7b図の差像は、結合さ
れているスピンだけを表示の為にフィルタにかけた時、
氷像93が実質的に消え(この場合は1/100以下に
抑圧され)、これに対して油像94は、それ以上ではな
いとしても同じ位に明白である。
第7c図には、健康な成人のボランチアの前脚の横断面
の多重量子フィルタ像の写真が示されている。実験の条
件は第7a図及び第7bに使った条件と同様であるが、
Trは1.OOOミリ秒に長くした。和像95は、筋肉
部分と脂肪成分にある両方の結合されていない水からの
相加的な応答をはっきりと示している。差像96では、
単一量子の水応答信号が排除され、実質的に多重量子の
脂質及び骨髄の応答信号だけが目に見える情報に寄与す
る。
例として、量子コヒーレンスを発生することによって、
NMR作像及び/又は分光法で望ましくない共鳴を抑圧
する為のこの発明の方法並びにパルス順序の幾つかの実
施例を説明したが、当業者には、これまでの説明から種
々の変更が考えられよう。従って、この発明の範囲が、
こ\に具体的に説明した細部及び道具によって制限され
ず特許請求の範囲のみによって限定されることを承知さ
れだい。
【図面の簡単な説明】
’1518図及び第1b図は、結合されていないスピン
1/2の原子核のスピン状態線図とこのスピンのNMR
応答特性を周波数領域で示す振幅のグラフ、 第2a図及び第2b図は1対の結合されたスピン1/2
の原子核のスピン状態線図と、この1対のスピンのNM
R応答特性を周波数領域で示す振幅のグラフ、 第3a図及び第3bはAB3スピン系内の4個のスピン
1/2の原子核のスピン状態線図とこのAB3スピン系
のNMR応答特性を周波数領域で示す振幅のグラフ、 第4図はこの発明の現在好ましいと考えられる実施例の
単一量子以外のコヒーレンスを発生する副順序及びそれ
と共に使われるスピン捩れ影信号作像順序中の磁界勾配
、無線周波励振信号、データ収集信号及びゲート信号の
時間軸を合せた1組の信号波形図、 第4a図乃至第4f図は1組のスピン密度マトリクスを
示す図で、第4図の順序の選ばれた時刻に於けるABス
ピン系の状態を示す。 第4g図は第4b図のスピン密度マトリクスから予想さ
れる応答の周波数領域と振幅のグラフ、第5a図及び第
5bはABスピン系に対し、第1及び第2の90RFパ
ルス信号のはじき角度に誤差の種々の大きさを用いた時
の、第4図の励振側順序から得られる二重量子及び単一
量子コヒーレンス信号の夫々の量を示す線図、 第6a図及び第6b図は同じ励振側順序であるが、第2
のRFパルス信号の位相を90”変更した場合のゼロ量
子/二重量子及び単一量子コヒーレンス信号の夫々の量
を示す線図、 第7a図及び第7b図は同じ原子核の異なる種目の百分
率が高いファントム物体から発生されたNMR像の写真
(X線写真に相当するものであり、図面代用写真を提出
する。)で、多重量子フィルタ作用を行なうこの発明の
方法によって得られる可視効果が相対的に異なる特性を
例示している。 第7C図はボランチアである被検体の前脚の1対の横方
向の図を示す写真で、多重量子コヒーレンスを持つ所望
の原子核種目の応答を強調する為に、単一量子コヒーレ
ンスを持つ望ましくない核物質種目の含有回が高い組織
からの応答を除くこの発明の効果を例示している。

Claims (1)

  1. 【特許請求の範囲】 1)互いにスピン結合された原子核を含むサンプルから
    受取ることが出来る所望のNMR信号を求めながら、核
    磁気共鳴(NMR)作像で不所望の応答信号を抑圧する
    方法に於て、 (a)励振信号順序に応答してサンプルからNMR応答
    信号を引出し、 (b)各々のNMR励振信号順序の前に、サンプル内の
    選ばれたスピン結合系のスピン・ポピュレーションのゼ
    ロ量子及び多重量子コヒーレンスの内の少なくとも選ば
    れた一方を発生する様に予定の形で選ばれた信号副順序
    を加え、 (c)該信号副順序中に少なくとも1つの検出信号を用
    いて、前記少なくとも一方の選ばれたコヒーレンスをサ
    ンプルから受信することが出来る所望のNMR応答信号
    に変換すると同時に、サンプルからの全応答信号中に存
    在する少なくとも1つの単一量子コヒーレンス信号の最
    終的な応答振幅を実質的に減少させ、 (d)前記少なくとも一方の選ばれたコヒーレンスに応
    答してサンプルの像を発生する工程を含む方法。 2)特許請求の範囲1)に記載した方法に於て、各々の
    NMR励振信号順序の前に信号副順序を加える工程(b
    )が、略90°の持続時間を持つ第1のRF磁界パルス
    信号を発生し、Kを0より大きい正の整数として、K・
    45°の持続時間を持つ第2のRF磁界パルス信号を発
    生し、前記90°パルス信号及び前記K・45°パルス
    信号の時間的な中点の間に、TをN/(4J)に略等し
    くし、Nを正の奇数、Jを所望のスピン系のスピン−ス
    ピン結合定数として、2Tに略等しい期間を設けること
    を含み、前記信号副順序中に少なくとも1つの検出信号
    を用いる工程が、前記K・45°パルス信号の時間的な
    中点から1/(2J)に略等しい検出期間T_d後に発
    生する時間的な中点を持つ、持続時間が略90°のRF
    磁界パルス信号として前記検出信号を発生することを含
    む方法。 3)特許請求の範囲2)に記載した方法に於て、工程(
    a)の前記NMR応答信号を引出す励振信号順序が、ス
    ライス選択性RF磁界パルス信号を用い、更に、N及び
    qを夫々正の整数として、前記検出信号の時間的な中点
    からnT=q+(1/2J)後の期間に、前記スライス
    選択性RFパルス信号の時間的な中点を位置ぎめする工
    程を含む方法。 4)特許請求の範囲3)に記載した方法に於て、スライ
    ス選択性RF信号をsinc変調形RF搬送波信号とし
    て発生する工程を含む方法。 5)特許請求の範囲2)に記載した方法に於て、更に、
    前記第1及び第2のRFパルス信号の時間的な中点の間
    の略中点にある時間的な中点を持ち、持続時間が略18
    0°の再集束RFパルス信号を発生する工程を含む方法
    。 6)特許請求の範囲5)に記載した方法に於て、ゼロ量
    子コヒーレンスによる応答を求め、更に、整数K=1に
    設定して、第2のRFパルス信号が略45°の持続時間
    を持つ様にする工程を含む方法。 7)特許請求の範囲6)に記載した方法に於て、更に、
    NMR装置の静磁界中に磁界勾配が実質的に存在しない
    状態で、前記第1、第2、再集束及び検出RFパルス信
    号の全てを発生する工程を含む方法。 8)特許請求の範囲5)に記載した方法に於て、多重量
    子コヒーレンスによる応答を求め、更に、整数K=2に
    設定して第2のRFパルス信号が略90°の持続時間を
    持つ様にする工程を含む方法。 9)特許請求の範囲8)に記載した方法に於て、更に、
    NMR装置の静磁界中に磁界勾配が実質的に存在しない
    状態で、前記第1、第2、再集束及び検出RFパルス信
    号の全てを発生する工程を含む方法。 10)特許請求の範囲8)に記載した方法に於て、更に
    、NMR装置の静磁界と平行な方向に、前記勾配磁界中
    に第1及び第2のパルス信号を用い、第1の磁界勾配パ
    ルス信号は、振幅及び持続時間の第1の積を持つと共に
    、発生する時間的な位置が前記第2のRFパルス信号よ
    り後で、前記検出RFパルス信号より前になる様にし、
    Nを1より大きい多重量子次数を表わす整数として、第
    2の磁界勾配パルス信号は、前記第1の磁界勾配パルス
    信号の振幅と持続時間の積のN倍の振幅と持続時間の積
    を持つと共に、発生する時間的な位置が検出RF信号パ
    ルスより後でNMR順序の開始より前である方法。 11)特許請求の範囲10)に記載した方法に於て、更
    に、NMR装置の静磁界と平行な方向に第1及び第2の
    追加の磁界勾配パルス信号を用い、各々の追加の勾配パ
    ルス信号は略同じ振幅Mを持ち、第1の追加のパルス信
    号は第1のRFパルス信号より後且つ再集束RFパルス
    信号より前に発生させると共に、再集束パルスの時間的
    な中点より時間T_gだけ前の時間的な中点を持つ様に
    し、第2の追加のパルス信号が再集束RFパルス信号よ
    り後且つ第2のRFパルス信号より前に発生すると共に
    、時間的な中点が、再集束パルスの時間的な中点から同
    じ時間T_9後になる様にする工程を含む方法。 12)特許請求の範囲11)に記載した方法に於て、更
    に、略180°の持続時間を持つと共に、時間的な中点
    が、第2のRFパルス信号及び検出RFパルス信号の時
    間的な中点の間の大体中点の時刻T_cになる第2の再
    集束パルスRF信号を用いる工程を含む方法。 13)特許請求の範囲10)に記載した方法に於て、更
    に、略180°の持続時間を持つと共に、時間的な中点
    が第2及び検出RFパルス信号の時間的な中点の間の大
    体中点の時刻T_cに発生する第2の再集束180°パ
    ルス信号を用いる工程を含む方法。 14)特許請求の範囲8)に記載した方法に於て、更に
    、NMR装置の静磁界と平行な方向に第1及び第2の追
    加の磁界勾配パルス信号を用い、各々の追加の勾配パル
    ス信号は略同じ振幅Mを持つ様にし、第1の追加のパル
    ス信号は第1のRFパルス信号より後且つ再集束RFパ
    ルス信号より前に発生して時間的な中点が再集束パルス
    の時間的な中点より時間T_gだけ前になる様にし、第
    2の追加のパルス信号は再集束RFパルス信号より後且
    つ第2のRFパルス信号より前に発生して、時間的な中
    点が再集束パルスの時間的な中点より同じ時間T_g後
    になる様にする工程を含む方法。 15)特許請求の範囲14)に記載した方法に於て、更
    に、略180°の持続時間を持ち、時間的な中点が第2
    及び検出RFパルス信号の時間的な中点の間の大体中点
    である時刻T_cに発生する第2の再集束パルスRF信
    号を用いる工程を含む方法。 16)特許請求の範囲8)に記載した方法に於て、略1
    80°の持続時間を持つと共に、時間的な中点が第2及
    び検出RFパルス信号の時間的な中点の間の大体中点の
    時刻T_cに発生する第2の再集束RFパルス信号を用
    いる工程を含む方法。 17)特許請求の範囲2)に記載した方法に於て、基準
    位相に対する第1及び第2のRFパルス信号の位相、及
    びサンプルから引出されたものであるが、検出する前の
    応答信号の位相の複数個(pを0より大きい正の整数と
    して、2^p)の異なる組合せに対し、前記信号副順序
    と励振信号順序を用いる工程を繰返し、各々の位相は、
    検出されたサンプルの応答中の少なくとも1つの単一量
    子応答信号を実質的に相殺する様に選ばれた予定のパタ
    ーンで、夫々量φ_a、φ_b、φ_rだけ移相する工
    程を含む方法。 18)特許請求の範囲17)に記載した方法に於て、p
    が10より小さい方法。 19)特許請求の範囲2)に記載した方法に於て、更に
    、検出信号及び受信した応答信号を含めて、励振信号順
    序の少なくとも1つの信号の位相を変え、変えた位相を
    操作して、単一量子コヒーレンス信号が実質的に相殺さ
    れる様にすると同時に、ゼロ量子コヒーレンス及び多重
    量子コヒーレンスの選ばれた一方に関連する受信した応
    答信号の成分が積極的に加算される様にする工程を含む
    方法。 20)特許請求の範囲1)に記載した方法に於て、更に
    検出信号及び受信した応答信号を含めて、励振信号順序
    の少なくとも1つの信号の位相を変え、変えた位相を操
    作して、単一量子コヒーレンス信号が実質的に相殺され
    る様にすると同時に、ゼロ量子コヒーレンス及び多重量
    子コヒーレンスの選ばれた一方に関連する受信した応答
    信号の成分が積極的に加算される様にする工程を含む方
    法。
JP62015361A 1986-01-27 1987-01-27 不所望の応答信号を抑圧する方法と装置 Granted JPS62194845A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US06/822,730 US4680546A (en) 1986-01-27 1986-01-27 Methods of, and pulse sequences for, the supression of undesired resonances by generation of quantum coherence in NMR imaging and spectroscopy
US822730 1986-01-27

Publications (2)

Publication Number Publication Date
JPS62194845A true JPS62194845A (ja) 1987-08-27
JPH044892B2 JPH044892B2 (ja) 1992-01-29

Family

ID=25236809

Family Applications (1)

Application Number Title Priority Date Filing Date
JP62015361A Granted JPS62194845A (ja) 1986-01-27 1987-01-27 不所望の応答信号を抑圧する方法と装置

Country Status (6)

Country Link
US (1) US4680546A (ja)
EP (1) EP0239724B1 (ja)
JP (1) JPS62194845A (ja)
DE (1) DE3766735D1 (ja)
FI (1) FI870295A (ja)
IL (1) IL80281A0 (ja)

Families Citing this family (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4701708A (en) * 1986-08-01 1987-10-20 General Electric Company Polarization transfer by selective homonuclear technique for suppression of uncoupled spins in NMR spectroscopy
GB8718723D0 (en) * 1987-08-07 1987-09-16 Hall L D Investigating sample using nmr
NL8701889A (nl) * 1987-08-12 1989-03-01 Philips Nv Volume selektieve spektroscopie door middel van gerefokusseerde echo's.
NL8702700A (nl) * 1987-11-12 1989-06-01 Philips Nv Werkwijze en inrichting voor automatische fasecorrectie van complexe nmr spectra.
DE3739856A1 (de) * 1987-11-25 1989-06-08 Philips Patentverwaltung Kernresonanz-spektroskopieverfahren
WO1989010553A1 (en) * 1988-04-20 1989-11-02 Laurence David Hall Investigating a sample using nmr
NL8900506A (nl) * 1989-03-02 1990-10-01 Philips Nv Magnetische resonantiewerkwijze voor het verkrijgen van geselekteerde spektra.
DE3914302C2 (de) * 1989-04-29 1994-01-13 Bruker Analytische Messtechnik Verfahren zur Aufnahme von Spinresonanzspektren
DE3914301A1 (de) * 1989-04-29 1990-10-31 Bruker Medizintech Verfahren zur aufnahme von spinresoneanzspektren und zur spinresonanz-bildgebung
EP0422170B1 (de) * 1989-04-29 1997-01-15 University Of Queensland Verfahren zur aufnahme von spinresonanzspektren
DE3914351A1 (de) * 1989-04-29 1990-10-31 Bruker Medizintech Verfahren zur aufnahme von spinresonanzspektren
DE3920433A1 (de) * 1989-06-22 1991-01-03 Philips Patentverwaltung Kernresonanzabbildungsverfahren
DE4018683A1 (de) * 1989-06-23 1991-01-10 Siemens Ag Schichtprofiloptimierung fuer ein mit einer gradientenechosequenz betriebenes kernspin-tomographiegeraet
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
NL8902003A (nl) * 1989-08-04 1991-03-01 Philips Nv 2-quantum selectieve nmr sequentie voor het selectief bepalen van een kernmagnetisatieverdeling van een metaboliet.
DE3925972A1 (de) * 1989-08-05 1991-02-07 Philips Patentverwaltung Kernresonanz-spektroskopieverfahren und anordnung zur durchfuehrung des verfahrens
FR2653901B1 (fr) * 1989-10-27 1991-12-27 Commissariat Energie Atomique Flacon a double paroi pour sonde a rmn.
US5162734A (en) * 1991-03-08 1992-11-10 Elbit-Ati, Ltd. Method and apparatus for obtaining an nmr signal having a preselected frequency domain
US5168229A (en) * 1991-04-02 1992-12-01 General Electric Company Multidimensional nmr spectroscopy using switched acquisition time gradients for multiple coherence transfer pathway detection
US5221899A (en) * 1991-04-29 1993-06-22 The Trustees Of Columbia University In The City Of New York Signal acquisition in magnetic resonance analysis
US5262723A (en) * 1991-07-09 1993-11-16 General Electric Company Method and apparatus for obtaining pure-absorption two-dimensional lineshape data for multidimensional NMR spectroscopy using switched acquisition time gradients
US5214380A (en) * 1991-09-30 1993-05-25 General Electric Company Velocity-resolved NMR spectroscopy
US5321359A (en) * 1993-03-29 1994-06-14 General Electric Company Selective imaging among three or more chemical species
US5627588A (en) * 1995-03-30 1997-05-06 International Business Machines Corporation Video signal processing stripe color demodulator system and method
US5510713A (en) * 1995-06-30 1996-04-23 General Electric Company Technique for selective fat saturation in MR imaging
US6144199A (en) * 1997-08-18 2000-11-07 Ramot University Authority For Applied Research And Industrial Development Ltd. Method for imaging vessel wall anatomy and strain
US6605942B1 (en) * 1998-03-20 2003-08-12 The Trustees Of Princeton University Method and apparatus for improvement of magnetic resonance imaging contrast
US7474095B2 (en) 2003-07-28 2009-01-06 University Of Southampton NMR spectroscopy using spin states with long lifetimes
US7777485B2 (en) * 2006-08-15 2010-08-17 General Electric Company Method for multiplexed MR tracking
US8583213B2 (en) * 2006-09-12 2013-11-12 General Electric Company Combined MR imaging and tracking
US9733326B2 (en) * 2007-11-07 2017-08-15 University of Pittsburgh—of the Commonwealth System of Higher Education Selective zero-quantum coherence transfer (Sel-ZQC) method for metabolite imaging in a poorly shimmed magnet field without susceptibility artifact
WO2011130633A2 (en) * 2010-04-16 2011-10-20 University Of Pittsburgh - Of The Commonwealth System Of Higher Education Simultaneous mapping of multiple chemicals with suppression of unwanted signals via molecular specific coherence (msc)-selmqc (selective multiple quantum coherence)
US9075122B2 (en) * 2009-11-06 2015-07-07 New York University Method, system and computer-accessible medium for providing multiple-quantum-filtered imaging
US8378679B2 (en) * 2009-12-04 2013-02-19 Siemens Medical Solutions Usa, Inc. Coherent signal acquisition system for MR imaging and spectroscopy
RU2477907C2 (ru) * 2010-12-23 2013-03-20 Государственное образовательное учреждение высшего профессионального образования "Российский университет дружбы народов" (РУДН) Способ формирования спиновых волн
CN110325268A (zh) * 2016-12-29 2019-10-11 通用电气公司 混合系统
CN110516811B (zh) * 2019-08-30 2023-04-18 南方科技大学 一种量子态确定方法、装置、设备和存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60168041A (ja) * 1983-12-14 1985-08-31 ゼネラル・エレクトリツク・カンパニイ Nmr作像方法と装置
JPS62142541A (ja) * 1985-12-12 1987-06-25 エヌ・ベ−・フイリツプス・フル−イランペンフアブリケン 磁気共鳴断層撮影方法および装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4134058A (en) * 1977-11-28 1979-01-09 Varian Associates, Inc. Selective detection of multiple quantum transitions in nuclear magnetic resonance
DE2833853C2 (de) * 1978-08-02 1985-11-21 Spectrospin AG, Fällanden, Zürich Verfahren zur Aufnahme von Kernspinresonanzspektren und Vorrichtung zur Durchführung des Verfahrens
US4238735A (en) * 1979-02-21 1980-12-09 Varian Associates, Inc. Indirect detection of nuclear spins of low gyromagentic ratio coupled to spins of high gyromagnetic ratio
DE3143625C2 (de) * 1981-11-04 1985-11-28 Spectrospin AG, Fällanden, Zürich Verfahren zur Aufnahme von Kernresonanzspektren von Molekülen in einem zweidimensionalen Frequenzbereich und Kernresonanzspektrometer zur Durchführung des Verfahrens

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60168041A (ja) * 1983-12-14 1985-08-31 ゼネラル・エレクトリツク・カンパニイ Nmr作像方法と装置
JPS62142541A (ja) * 1985-12-12 1987-06-25 エヌ・ベ−・フイリツプス・フル−イランペンフアブリケン 磁気共鳴断層撮影方法および装置

Also Published As

Publication number Publication date
JPH044892B2 (ja) 1992-01-29
US4680546A (en) 1987-07-14
EP0239724A3 (en) 1988-01-27
IL80281A0 (en) 1987-01-30
FI870295A0 (fi) 1987-01-23
DE3766735D1 (de) 1991-01-31
EP0239724B1 (en) 1990-12-19
FI870295A (fi) 1987-07-28
EP0239724A2 (en) 1987-10-07

Similar Documents

Publication Publication Date Title
JPS62194845A (ja) 不所望の応答信号を抑圧する方法と装置
KR880001529B1 (ko) 화학 이동성 nmr 주파수를 갖는 핵을 포함한 샘플의 영상화 방법
JPS6382351A (ja) Nmr分光方法
Clayton et al. 1H spectroscopy without solvent suppression: characterization of signal modulations at short echo times
JP2007322443A (ja) サンプルの核四極子共鳴テスト方法および装置、および核四極子共鳴を励起するためのパルス・シーケンス
De Graaf et al. B1-insensitive, single-shot localization and water suppression
US4238735A (en) Indirect detection of nuclear spins of low gyromagentic ratio coupled to spins of high gyromagnetic ratio
JPS6046447A (ja) 核磁気共鳴画像形成方法およびそれに用いる画像形成装置
JPS6047946A (ja) Nmr作像方法
JPS635734A (ja) 望ましくない応答信号を抑圧する方法
JPH0287050A (ja) 結合スピンを含む新陳代謝物質によるnmr信号の二次元スペクトルを作る方法と装置
Bottomley Basics of NMR
JPS635696B2 (ja)
Dumoulin et al. Water suppression in 1H magnetic resonance images by the generation of multiple‐quantum coherence
JPH0213430A (ja) 結合ピストンを含む新陳代謝物質による局在化スペクトルのnmr信号を発生する方式
Mackin et al. Phase-sensitive two-dimensional HMQC and HMQC-TOCSY spectra obtained using double pulsed-field-gradient spin echoes
US5111143A (en) Magnetic resonance spectroscopy method and apparatus for performing the method
JPH01156647A (ja) 磁気共鳴分光法とそのデバイス
JPH0146824B2 (ja)
EP0213614A2 (en) Method for obtaining nuclear magnetic resonance information data
JP3212119B2 (ja) 磁気共鳴分光装置
EP0398440B1 (en) Method of suppressing water resonance in a magnetic proton resonance spectrum, and device for performing such a method
JP3169640B2 (ja) 2次元核磁気共鳴分光用方法
US4712066A (en) Method for the selective excitation of a volume in an object
US4706028A (en) NMR imaging method