JPS63214246A - 核磁気共鳴信号を用い、対象物から高速度で影像情報を得る装置 - Google Patents
核磁気共鳴信号を用い、対象物から高速度で影像情報を得る装置Info
- Publication number
- JPS63214246A JPS63214246A JP62306827A JP30682787A JPS63214246A JP S63214246 A JPS63214246 A JP S63214246A JP 62306827 A JP62306827 A JP 62306827A JP 30682787 A JP30682787 A JP 30682787A JP S63214246 A JPS63214246 A JP S63214246A
- Authority
- JP
- Japan
- Prior art keywords
- pulse
- gradient
- high speed
- image information
- space
- 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
Links
- 238000000034 method Methods 0.000 title claims description 49
- 238000005481 NMR spectroscopy Methods 0.000 title claims description 20
- 238000003384 imaging method Methods 0.000 title description 12
- 238000005070 sampling Methods 0.000 claims description 22
- 230000000694 effects Effects 0.000 claims description 14
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 14
- 239000000126 substance Substances 0.000 claims description 11
- 230000001629 suppression Effects 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000005686 electrostatic field Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 230000005311 nuclear magnetism Effects 0.000 claims 1
- 230000036961 partial effect Effects 0.000 description 11
- 230000005284 excitation Effects 0.000 description 8
- 230000002829 reductive effect Effects 0.000 description 8
- 238000002592 echocardiography Methods 0.000 description 7
- 230000003068 static effect Effects 0.000 description 7
- 230000010363 phase shift Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 230000009467 reduction Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000002238 attenuated effect Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000013016 damping Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 150000002632 lipids Chemical group 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 208000032005 Spinocerebellar ataxia with axonal neuropathy type 2 Diseases 0.000 description 1
- 238000004026 adhesive bonding Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 208000033361 autosomal recessive with axonal neuropathy 2 spinocerebellar ataxia Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000001225 nuclear magnetic resonance method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
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/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
-
- 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
-
- 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/561—Image 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
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
-
- 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/561—Image 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
-
- 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/561—Image 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
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
本発明は、核磁気共鳴信号を用い、対象物から高速度で
影像情報を得る方法に関し、より詳しくいうと、約10
0ms以下の速度で影像をつくる方法に関する。
影像情報を得る方法に関し、より詳しくいうと、約10
0ms以下の速度で影像をつくる方法に関する。
いろいろな理由から、非常に短い時間内に影像を得るこ
とができる磁気共鳴影像(MRI)技術の確立が望まれ
ている。現在、影像診断分野において、MHIの発展を
はばんでいる要因のうち主なものは、原価効率に係る問
題のほかに、不随意的生理運動、特に心臓の動きがある
ため、臓器系の診断への適用が制約されてしまう問題で
ある。
とができる磁気共鳴影像(MRI)技術の確立が望まれ
ている。現在、影像診断分野において、MHIの発展を
はばんでいる要因のうち主なものは、原価効率に係る問
題のほかに、不随意的生理運動、特に心臓の動きがある
ため、臓器系の診断への適用が制約されてしまう問題で
ある。
長いデータの受取り時間(通常、数分)を小さくできれ
ば、前述の問題は、いずれも(特に、後者)改善される
。つまり、それにより、患者の全診察時間が短縮され、
患者の処理量が増大するため、原価効率の適性化がはか
れる。もし、全走査時間を、心臓周期の極一部程度に縮
めることができれば、心臓の運動を含むすべての生理学
的運動による本来の動きは、無視される。
ば、前述の問題は、いずれも(特に、後者)改善される
。つまり、それにより、患者の全診察時間が短縮され、
患者の処理量が増大するため、原価効率の適性化がはか
れる。もし、全走査時間を、心臓周期の極一部程度に縮
めることができれば、心臓の運動を含むすべての生理学
的運動による本来の動きは、無視される。
従来得られている高速の全生体影像は、立体的分解能が
よくないし、またSハ比が低すぎる。(本明細書におい
て使用する「高速」とは、完全な断層撮影陽子像の受取
りに対し、約100msを意味する。)この低品質は、
−面、使用される相当に低い磁場の強さの反対にある。
よくないし、またSハ比が低すぎる。(本明細書におい
て使用する「高速」とは、完全な断層撮影陽子像の受取
りに対し、約100msを意味する。)この低品質は、
−面、使用される相当に低い磁場の強さの反対にある。
しかし、公知のように、S/N比は、磁場の強さととも
にほぼ直線的に増大するため、良質の像は、原則として
、より高い強度の磁場で得られるようになる。
にほぼ直線的に増大するため、良質の像は、原則として
、より高い強度の磁場で得られるようになる。
高強度磁場において、高速の影像形成を行なう際に、い
くつかの問題が派生する。本発明は、これらの問題を解
決することである。
くつかの問題が派生する。本発明は、これらの問題を解
決することである。
磁気共鳴影像に対し高磁場を用いる際の第1の大きな問
題は、磁場における不均質性が、実際に不可避的である
点と、観察対象の磁気共鳴影像に幾何学的位相歪みがも
たらされる点である。このような磁気の不均質性は、2
つの異なる供給源からもたらされる。
題は、磁場における不均質性が、実際に不可避的である
点と、観察対象の磁気共鳴影像に幾何学的位相歪みがも
たらされる点である。このような磁気の不均質性は、2
つの異なる供給源からもたらされる。
第1に、静磁場は、特に、磁場強度に比例して変化する
固有の不均質性を持っている。実際問題として、このよ
うな不均質性は、それを除去することが極めて困難であ
り、かつ膨大な費用を必要とするため、静磁場において
は不可避的であると見做されている。
固有の不均質性を持っている。実際問題として、このよ
うな不均質性は、それを除去することが極めて困難であ
り、かつ膨大な費用を必要とするため、静磁場において
は不可避的であると見做されている。
第2に、投影される対象物の磁気感受性が、システムに
おける磁気の不均質性に対して一層の拍車をかけること
である。
おける磁気の不均質性に対して一層の拍車をかけること
である。
更に、高い周波数では、人体のような弱R@体を貫通す
る際、RF磁場は、減衰するとともに、位相シフトを受
ける。この影響は、85MHzにおいて明白に現われる
。この値は、約2テスラでの陽子共鳴周波数である。
る際、RF磁場は、減衰するとともに、位相シフトを受
ける。この影響は、85MHzにおいて明白に現われる
。この値は、約2テスラでの陽子共鳴周波数である。
高強度磁場において影像がつくられる際の第2の大きな
問題点は、陽子の「化学シフト(chemical 5
hift)J効果である。この種の問題は、像における
あらゆる点からの信号が、一般に、2つの近接した周波
数、一つは水の陽子からのもの、もう一つは、脂質の陽
子からのものを含むために発生する。
問題点は、陽子の「化学シフト(chemical 5
hift)J効果である。この種の問題は、像における
あらゆる点からの信号が、一般に、2つの近接した周波
数、一つは水の陽子からのもの、もう一つは、脂質の陽
子からのものを含むために発生する。
MHIにおいて使用される影像手順は、特定の容積要素
から得られる信号を、固有の周波数に基づく最終像にお
ける対応の画素へ入れるようになっているため、組織の
同じ容積要素における脂質および水から得られる2つの
異なる周波数は、異なる画素へ誤って入れられてしまう
。それが、2重像を生じさせる。
から得られる信号を、固有の周波数に基づく最終像にお
ける対応の画素へ入れるようになっているため、組織の
同じ容積要素における脂質および水から得られる2つの
異なる周波数は、異なる画素へ誤って入れられてしまう
。それが、2重像を生じさせる。
この問題は、通常、低強度磁場にするか、遅速影像法を
用い、像における空間識別を与える磁界勾配強度を、勾
配が画素ごとに重ねら九るバンド幅が、脂肪および水の
陽子共鳴の周波数分離幅より実質的に大きくなる程度に
簡単に増大させることによって、解決される。
用い、像における空間識別を与える磁界勾配強度を、勾
配が画素ごとに重ねら九るバンド幅が、脂肪および水の
陽子共鳴の周波数分離幅より実質的に大きくなる程度に
簡単に増大させることによって、解決される。
高速影像の場合、画素ごとに、このような高いバンド幅
を使用することは、実行しにくく、避けるべきである。
を使用することは、実行しにくく、避けるべきである。
また、バンド幅は、磁場強度とともに必然的に増大され
るようになり、それによって、あらゆるSハ比を減らす
ことが、磁場強度の平方根に依存して増大して行く。
るようになり、それによって、あらゆるSハ比を減らす
ことが、磁場強度の平方根に依存して増大して行く。
強い磁場を用いる高速磁気共鳴影像法における上記の問
題点を解消するため、本発明者は、対象物から磁気影像
を迅速に取得しうる方法および装置の発明をした。
題点を解消するため、本発明者は、対象物から磁気影像
を迅速に取得しうる方法および装置の発明をした。
対象物は、第1の軸(本明細書において、Z軸と呼ぶ。
)に沿う静磁場に設定される。このZ軸に沿う磁場は、
約0.5〜5テスラの範囲の強度を有する。
約0.5〜5テスラの範囲の強度を有する。
核スピンは、対象物の選択面における像領域において、
Z方向における静磁場に第1の磁界勾配(スライス選択
勾配と呼ぶ。)を重畳させ、同時に分光的に変えられた
RF励起パルスを加えることによって励起され、そこで
、自由誘導崩壊信号が、励起核スピンによってつくられ
る。
Z方向における静磁場に第1の磁界勾配(スライス選択
勾配と呼ぶ。)を重畳させ、同時に分光的に変えられた
RF励起パルスを加えることによって励起され、そこで
、自由誘導崩壊信号が、励起核スピンによってつくられ
る。
RF励起が持続している間、核スピンの位相ずれが、選
択面の幅にわたって生ずるため、核スピンの位相ずれの
戻しが必要になる。位相ずれの戻しは、第1の磁界勾配
を逆転させるか、180°パルス(後段で説明する)を
加えた後で、同じ第1の勾配波形を再度加えることによ
って達成される。
択面の幅にわたって生ずるため、核スピンの位相ずれの
戻しが必要になる。位相ずれの戻しは、第1の磁界勾配
を逆転させるか、180°パルス(後段で説明する)を
加えた後で、同じ第1の勾配波形を再度加えることによ
って達成される。
第2および第3の磁気勾配(これらを符号化勾配および
読出し勾配と呼ぶ。)は、それぞれ、X方向およびY方
向に加えられる。通常使用される如く、x、Yおよび2
の方向は、互いに直交している。選択領域は、X−Y面
内にあり、この面は、通常、横断面と呼ばれる。
読出し勾配と呼ぶ。)は、それぞれ、X方向およびY方
向に加えられる。通常使用される如く、x、Yおよび2
の方向は、互いに直交している。選択領域は、X−Y面
内にあり、この面は、通常、横断面と呼ばれる。
時間変動型の位相符号化勾配および読出し勾配は、空間
変動式に一時的な核信号の位相を変化させる働き、つま
り、データ受取り時、特定の時間にサンプルされる空間
周波数定義域における点を決定する働きをする。結果と
して得られるに空間におけるサンプル点の軌線は、k軌
線として知られている。
変動式に一時的な核信号の位相を変化させる働き、つま
り、データ受取り時、特定の時間にサンプルされる空間
周波数定義域における点を決定する働きをする。結果と
して得られるに空間におけるサンプル点の軌線は、k軌
線として知られている。
デコード・プロセスは、いくつかの初期処理段階を含み
、それに引き続き、2次元フーリエ変換が行なわれて、
最初の空間分布の概略である不連続像を生成する。
、それに引き続き、2次元フーリエ変換が行なわれて、
最初の空間分布の概略である不連続像を生成する。
サンプリングが開始される点をに空間に設定するため、
サンプリングの前に、エンコード・シーケンスが生ずる
。それから、RF倍信号、静電場における不均質性によ
り位相がずれているあらゆる従属的核スピンの位相ずれ
を戻すべく、対象物に加えられる。次に、データは、位
相符号化勾配および読出し勾配が交互に加えられるシー
ケンスによりサンプルされる。
サンプリングの前に、エンコード・シーケンスが生ずる
。それから、RF倍信号、静電場における不均質性によ
り位相がずれているあらゆる従属的核スピンの位相ずれ
を戻すべく、対象物に加えられる。次に、データは、位
相符号化勾配および読出し勾配が交互に加えられるシー
ケンスによりサンプルされる。
本発明の好適実施例においては、読出し勾配は、連続正
弦波として加えられ、位相符号化勾配と読出し勾配との
間に若干の重畳が生ずるようになっている。符号化時間
、およびRFパルスのアプリケーションを含む全サンプ
リング・シーケンスは、約100ms以下影像時間と共
に生じる。
弦波として加えられ、位相符号化勾配と読出し勾配との
間に若干の重畳が生ずるようになっている。符号化時間
、およびRFパルスのアプリケーションを含む全サンプ
リング・シーケンスは、約100ms以下影像時間と共
に生じる。
既に述べた技術を使用しての走査時間については、S/
N比または空間分解能におけるコスト削減の点だけで、
相当の改良がなされている。身体の動きを問題にしてい
ない情況にあっては、このような削減は、容認できない
。かかる状態においては、データ取得に多くの時間をか
け、かつ空間分解能を取り戻すことが好ましい。
N比または空間分解能におけるコスト削減の点だけで、
相当の改良がなされている。身体の動きを問題にしてい
ない情況にあっては、このような削減は、容認できない
。かかる状態においては、データ取得に多くの時間をか
け、かつ空間分解能を取り戻すことが好ましい。
そのため、既に述べた高速影像技術に関し、次の変形を
開発し、医者が、選択でき、かつ臨床使用時に調整しう
る空間分解能、全観察時間、およびS/N比の中から、
1組の選択ができるようにした。
開発し、医者が、選択でき、かつ臨床使用時に調整しう
る空間分解能、全観察時間、およびS/N比の中から、
1組の選択ができるようにした。
もし、磁場の不均質性の結果として、または他の原因に
よりもたらされる位相歪みが全くないとしたら、k軌線
が、単に、k空間の正か負の半分におけるラインを覆え
ば十分である。k空間の反対側半分においてサンプルさ
れた対応データラインは、同じように共役データを含む
ため、データの冗長度を生ずるだけである。しかし、位
相歪みが存在する場合、k空間の両側半分における対応
データラインは、同じではない。
よりもたらされる位相歪みが全くないとしたら、k軌線
が、単に、k空間の正か負の半分におけるラインを覆え
ば十分である。k空間の反対側半分においてサンプルさ
れた対応データラインは、同じように共役データを含む
ため、データの冗長度を生ずるだけである。しかし、位
相歪みが存在する場合、k空間の両側半分における対応
データラインは、同じではない。
サンプリング軌線を、k空間の半分以上に拡張すること
から得られる情報が、位相のエラーを補償するのに使用
できることを発見した。必要なサンプリングの拡張は、
導入される位相のエラーの程度に左右される。例えば、
位相エラーが、ノ像の上でゆっくり変化する場合(つま
り、エラーが、低い空間周波数をとっている場合)、軌
線を更に拡張してサンプルを行なう必要がない。
から得られる情報が、位相のエラーを補償するのに使用
できることを発見した。必要なサンプリングの拡張は、
導入される位相のエラーの程度に左右される。例えば、
位相エラーが、ノ像の上でゆっくり変化する場合(つま
り、エラーが、低い空間周波数をとっている場合)、軌
線を更に拡張してサンプルを行なう必要がない。
適切な位相補正を行なうのに必要な付加的ラインのみを
サンプルすることによって、全受取り時間は、最低に抑
えられ、最初の励起パルスに続くデータサンプリングの
開始時間は、最小になる。
サンプルすることによって、全受取り時間は、最低に抑
えられ、最初の励起パルスに続くデータサンプリングの
開始時間は、最小になる。
後者の結果は、T2の結果的な信号崩壊が、それにより
最小になる、という点で望ましい。
最小になる、という点で望ましい。
更に、もし、T2崩壊によるSハ比損失が、より少ない
データラインの受取りから結果として生ず=16− るSハ比損失よりも大きければ、出力のS/N比は、最
適になる。
データラインの受取りから結果として生ず=16− るSハ比損失よりも大きければ、出力のS/N比は、最
適になる。
部分に空間影像の一変形例において、k空間の部分受取
りは、X面かY面のいずれかの面でなされる。部分走査
におけるサンプリング軌線は、100%以下であるが、
k空間の半分以上に拡張するのが好ましい。
りは、X面かY面のいずれかの面でなされる。部分走査
におけるサンプリング軌線は、100%以下であるが、
k空間の半分以上に拡張するのが好ましい。
サンプリング軌線を元のものの」二に拡張することによ
り取得される情報は、導入されるあらゆる位相エラーを
補償するのに使用される。k空間の各ラインは、約40
0msを必要とするので、部分的にχ受取りにおけるに
χラインを省きうろことは、影像時間を短縮する。
り取得される情報は、導入されるあらゆる位相エラーを
補償するのに使用される。k空間の各ラインは、約40
0msを必要とするので、部分的にχ受取りにおけるに
χラインを省きうろことは、影像時間を短縮する。
部分的に空間受取りは、同数の点が、以前の全に空間受
取りにおけるようにして集められるよう行なわれ、従っ
て、信号のバンド幅を維持させながら、空間周波数応答
を増大させることができる。
取りにおけるようにして集められるよう行なわれ、従っ
て、信号のバンド幅を維持させながら、空間周波数応答
を増大させることができる。
そのため、この技術的方法は、「ワン・ショット」式に
より高い分解能の像を生成する。しかし、若干低いSハ
比を伴う。
より高い分解能の像を生成する。しかし、若干低いSハ
比を伴う。
本発明による別の変形例において、2つ、4つ、あるは
それ以上の部分的受取りは、過渡応答(TR)の相互受
取り遅延により、k空間の異なる部分」二で行なわれる
。取得された部分は、フーリエ変換の前に、k空間のパ
ッチワーク、即ち「モザイク」的に合わせられる。
それ以上の部分的受取りは、過渡応答(TR)の相互受
取り遅延により、k空間の異なる部分」二で行なわれる
。取得された部分は、フーリエ変換の前に、k空間のパ
ッチワーク、即ち「モザイク」的に合わせられる。
そのため、この「モザイク走査」技術におけるに空間の
有効範囲は、100%となりうるが、多重受取りが必要
である。S/N比損失は、(a)k空間のすべてが充た
されるため、(b)低空間周波数が効果的に信号平均化
されるため、かつ(c)いくつかの受取りが一つの像に
好影響を与えているため、部分的に取り戻される。
有効範囲は、100%となりうるが、多重受取りが必要
である。S/N比損失は、(a)k空間のすべてが充た
されるため、(b)低空間周波数が効果的に信号平均化
されるため、かつ(c)いくつかの受取りが一つの像に
好影響を与えているため、部分的に取り戻される。
既に言及した化学シフトの問題を排除するため、2つの
好適実施例について説明する。
好適実施例について説明する。
第1の実施例において、スライス選択を行なう前に、R
F磁場の狹帯域抑制パルスが加えられ、引き続き、スピ
ン位相ずれにおいて、影像シーケンスにおけるあらゆる
後続勾配より相当大きい磁界勾配のパルスが続く。これ
は、一つの化学成分、脂質、または水から、信号を有効
に排除する。
F磁場の狹帯域抑制パルスが加えられ、引き続き、スピ
ン位相ずれにおいて、影像シーケンスにおけるあらゆる
後続勾配より相当大きい磁界勾配のパルスが続く。これ
は、一つの化学成分、脂質、または水から、信号を有効
に排除する。
第2の実施例において、]80°再焦点パルスは、その
周波数スペクトルおよび「オフ共鳴」効果が、成分の一
方若しくは他方の位相ずれの戻りを排除するべく、結合
しうるよう選択される。
周波数スペクトルおよび「オフ共鳴」効果が、成分の一
方若しくは他方の位相ずれの戻りを排除するべく、結合
しうるよう選択される。
次に、添付図面を参照して、本発明の実施例を詳細に説
明する。
明する。
第1図は、磁石(14)内に配置された人体等の被写体
(12)内の選択された部位(10)の磁気共鳴(MR
)影像を得るための図である。第1図に示すよう1、二
選択された部位(10)は、コイル(16)によって発
生されるx、y、X方向の磁界の勾配の向きによっ゛τ
画定さ九る平面内1.二位置する。
(12)内の選択された部位(10)の磁気共鳴(MR
)影像を得るための図である。第1図に示すよう1、二
選択された部位(10)は、コイル(16)によって発
生されるx、y、X方向の磁界の勾配の向きによっ゛τ
画定さ九る平面内1.二位置する。
以下の説明のために、静磁界(BO)の方向は、通常の
慣習のように、Z軸を決定する。また、第1の磁界勾配
(スライス−セレクション勾配)は、2方向に適用され
、第2の磁界勾配(位相符号化勾配)は、X方向に適用
さ九、第3の磁界勾配(請出し勾配)は、Y方向に適用
される。、x、y、z方向は、互いに直交するように決
められている。
慣習のように、Z軸を決定する。また、第1の磁界勾配
(スライス−セレクション勾配)は、2方向に適用され
、第2の磁界勾配(位相符号化勾配)は、X方向に適用
さ九、第3の磁界勾配(請出し勾配)は、Y方向に適用
される。、x、y、z方向は、互いに直交するように決
められている。
従って、選択された部位(10)は、X−Y平面に入る
。人体等の被写体の縦軸が、Z軸と平行ならば、このX
−Y平面は、通常、交軸面と呼ばれる。
。人体等の被写体の縦軸が、Z軸と平行ならば、このX
−Y平面は、通常、交軸面と呼ばれる。
前記の第1、第2及び第3磁界勾配X、Y−Z方向への
座標指定は、絶対ではない。例えば、第1磁界勾配(ス
ライス−セレクション)が、それぞれ又もしくはY方向
に適用され、第2(位相符号化)及び第3(読出し)勾
配が、それぞれXとYもしくはYと2方向に適用される
と仮定することによって、X−7もしくはY−Z平面に
おける画定部位を選択しうる。選択部位は、冠状面もし
くは矢状面となる。
座標指定は、絶対ではない。例えば、第1磁界勾配(ス
ライス−セレクション)が、それぞれ又もしくはY方向
に適用され、第2(位相符号化)及び第3(読出し)勾
配が、それぞれXとYもしくはYと2方向に適用される
と仮定することによって、X−7もしくはY−Z平面に
おける画定部位を選択しうる。選択部位は、冠状面もし
くは矢状面となる。
更に、第1、第2及び第3勾配は、通常、必らずしも、
前記のX、YもしくはZ方向に沿って位置する必要はな
い。第1、第2及び第3の勾配を常に直交させるために
、これら3つの勾配を、xl、Yl及びZ□の変換座標
系を向くように回転させることができ、このXIY、Z
□系の中の任意に指向される「斜め」の平面が、選択さ
れ、矢状面と垂直な視点=20= から映し出される。
前記のX、YもしくはZ方向に沿って位置する必要はな
い。第1、第2及び第3の勾配を常に直交させるために
、これら3つの勾配を、xl、Yl及びZ□の変換座標
系を向くように回転させることができ、このXIY、Z
□系の中の任意に指向される「斜め」の平面が、選択さ
れ、矢状面と垂直な視点=20= から映し出される。
第17、第2及び第3の勾配は、必ずしも互いに直角で
ある必要はない。その場合、例えば、任意に指向される
斜めの平面が選択され、ついで斜めに映し出される。す
なわち、矢状面に対して斜めの視点から映し出される。
ある必要はない。その場合、例えば、任意に指向される
斜めの平面が選択され、ついで斜めに映し出される。す
なわち、矢状面に対して斜めの視点から映し出される。
これらの座標割当の変化は、当業者には公知である。
前記の磁界勾配の特定の座標割当のために、勾配Gx、
Gy及びGzは、次のように決定される。すな本発明
によって研究されている対象物は、常にZ軸に沿った連
続磁界(B、)に左右される。選択された磁界の強度は
、十分な信号対雑音比(SNR)が、十分な空間解像度
を有して最終影像の中で得られるのに十分な強さを持た
なければならない。更に、十分な空間解像度を有する場
所は、位相符号化と読出しの勾配の強さを必要とし、読
出し勾配は、信号を取得する無線周波数アンテナの所望
の信号帯域幅を必要とし、この信号は、最適のSNRに
対して、できるだけ高いのを保持しなければならない。
Gy及びGzは、次のように決定される。すな本発明
によって研究されている対象物は、常にZ軸に沿った連
続磁界(B、)に左右される。選択された磁界の強度は
、十分な信号対雑音比(SNR)が、十分な空間解像度
を有して最終影像の中で得られるのに十分な強さを持た
なければならない。更に、十分な空間解像度を有する場
所は、位相符号化と読出しの勾配の強さを必要とし、読
出し勾配は、信号を取得する無線周波数アンテナの所望
の信号帯域幅を必要とし、この信号は、最適のSNRに
対して、できるだけ高いのを保持しなければならない。
妥当なSNRを約10=1〜1.00:1の範囲に、成
人にとっての妥当な空間解像度を、約5mm〜]−rr
mの範囲に収めるには、本出願人の方法を使用すると、
約0.5〜5.0テスラの範囲の強さの磁界が必要とな
る。
人にとっての妥当な空間解像度を、約5mm〜]−rr
mの範囲に収めるには、本出願人の方法を使用すると、
約0.5〜5.0テスラの範囲の強さの磁界が必要とな
る。
本発明の好適実施例の磁界の強さは、2.0テスラであ
る。この磁界の強さは、磁界の不均一性が対象物の導入
によって起こり、この不均一性は、後で説明される技術
によって補われるのに十分な強さである。
る。この磁界の強さは、磁界の不均一性が対象物の導入
によって起こり、この不均一性は、後で説明される技術
によって補われるのに十分な強さである。
強力な磁界は、特定の分子化合物から発生する信号のス
ペクトル分解の化学シフト問題を解決する。スペクトル
分解を行なうための技術も、後で説明する。
ペクトル分解の化学シフト問題を解決する。スペクトル
分解を行なうための技術も、後で説明する。
第2図示のパルスシーケンスにおける符号(X□)〜(
X4)は、時間間隔を示す。時間間隔(xl)において
は、核スピンは、研究されている対象物の平面内の選択
された部位(10)の中で励起される。これは、静磁界
の上に、例えば、望ましくは、(114)のような波形
を有する第2磁界勾配Gzを重ねることによって行なわ
れる。
X4)は、時間間隔を示す。時間間隔(xl)において
は、核スピンは、研究されている対象物の平面内の選択
された部位(10)の中で励起される。これは、静磁界
の上に、例えば、望ましくは、(114)のような波形
を有する第2磁界勾配Gzを重ねることによって行なわ
れる。
分光的に適合される90°高周波数(RF)パルス(1
,16)は、勾配Gzが定量的に適用されている間に適
用される。90°以下もしくは90°以上のRFパルス
を適用することもできるが、90°のパルスが望ましい
。このRF励起の際には、核スピンの位相のずれが、選
択された平面の幅で発生するおそれがあり、ついで、最
大信号振幅および最小位相歪を得るために位相のずれを
戻す必要がある。
,16)は、勾配Gzが定量的に適用されている間に適
用される。90°以下もしくは90°以上のRFパルス
を適用することもできるが、90°のパルスが望ましい
。このRF励起の際には、核スピンの位相のずれが、選
択された平面の幅で発生するおそれがあり、ついで、最
大信号振幅および最小位相歪を得るために位相のずれを
戻す必要がある。
位相のす九を戻すには、波形(114)に示された第1
の磁界勾配を(114a)の部分で反転させて行なう。
の磁界勾配を(114a)の部分で反転させて行なう。
あるいは、位相のすわの戻しは、例えば時間間隔(X4
)(180°パルス適用後、後で説明する)まで遅延さ
せることができ、その場合は、(114b)で示された
波形を使用することができる。Gz用の別の適切な波形
、及び90°パルス用の分光的に適合される各種の略図
は、当業者には容易に理解できるであろう。
)(180°パルス適用後、後で説明する)まで遅延さ
せることができ、その場合は、(114b)で示された
波形を使用することができる。Gz用の別の適切な波形
、及び90°パルス用の分光的に適合される各種の略図
は、当業者には容易に理解できるであろう。
本発明の符号化シーケンスは、時間間隔(x2)の間に
発生する。信号の符号化は、データ取得が始まるに一ス
ペースの点(すなわち、第2図のパルスシーケンスにお
いては共役点)の測定に影響を及ぼす。
発生する。信号の符号化は、データ取得が始まるに一ス
ペースの点(すなわち、第2図のパルスシーケンスにお
いては共役点)の測定に影響を及ぼす。
図示の特定の波形(120a)(120b)及び(13
0)は、絶対的ではない。第13勾配〔;χ及び第3勾
配Gyに適した多数の波形があることは、当業者には明
白であろう。例えば、第2勾配の波形は、(120b)
で示すように台形でもよい。あるいは、類似の部分(A
2)は、パルスシーケンスの時間間隔(x4)で所望さ
れるのと同一の、複数の短い勾配パルスによって掃引さ
れうる。
0)は、絶対的ではない。第13勾配〔;χ及び第3勾
配Gyに適した多数の波形があることは、当業者には明
白であろう。例えば、第2勾配の波形は、(120b)
で示すように台形でもよい。あるいは、類似の部分(A
2)は、パルスシーケンスの時間間隔(x4)で所望さ
れるのと同一の、複数の短い勾配パルスによって掃引さ
れうる。
基本パラメーターは、第2の磁界勾配に対しては、(A
、)で示された波形の下エリア、第3の磁界勾配に対し
ては、(A3)で示された波形の下のエリアである。と
いうのは、これらのエリアは、k空間軌線と数学的に関
連があるからである。従って、実際的な要件は、時間間
隔(x4)で特定の波形で示され、符号化期間(X2)
で使用される波形は、多数の適当な形状のいずれでもよ
い。
、)で示された波形の下エリア、第3の磁界勾配に対し
ては、(A3)で示された波形の下のエリアである。と
いうのは、これらのエリアは、k空間軌線と数学的に関
連があるからである。従って、実際的な要件は、時間間
隔(x4)で特定の波形で示され、符号化期間(X2)
で使用される波形は、多数の適当な形状のいずれでもよ
い。
第3図は、l(空間、すなわち空間周波数定義域を示し
ており、k軌線が見本となる線が、離散の空間周波数座
標に相当する。垂直軸に2は、第2磁界勾配(位相符号
化)に相当する空間周波数座標を表わす。水平軸に、は
、第3磁界勾配(読出し)の空間周波数座標を表わす。
ており、k軌線が見本となる線が、離散の空間周波数座
標に相当する。垂直軸に2は、第2磁界勾配(位相符号
化)に相当する空間周波数座標を表わす。水平軸に、は
、第3磁界勾配(読出し)の空間周波数座標を表わす。
RF励起間隔(X□)の直後で、か符号化の前の時間1
=16は、k空間内のに軌線の原点に2=に3=0を定
義し、第3図にAで示す。時間間隔(x2)の符号化シ
ーケンスは、軌線を描き、軌線の時間依存性は、実際の
波形(120a) (120b) (130)に左右さ
れる。
=16は、k空間内のに軌線の原点に2=に3=0を定
義し、第3図にAで示す。時間間隔(x2)の符号化シ
ーケンスは、軌線を描き、軌線の時間依存性は、実際の
波形(120a) (120b) (130)に左右さ
れる。
しかし、符号化期間(X2)の終端における軌線の位置
は、エリア(A2) (A3)にのみ左右され、かつこ
れらのエリアは、任意に設定できる。エリア(A2)(
A3)の好適設定例は、前符号化期間(X、)の終端に
おける軌線の位置が、第3図にBで示すように、k、=
(A2)i= −Mかつに3(A3 )J = N
/ 2となる場合である。
は、エリア(A2) (A3)にのみ左右され、かつこ
れらのエリアは、任意に設定できる。エリア(A2)(
A3)の好適設定例は、前符号化期間(X、)の終端に
おける軌線の位置が、第3図にBで示すように、k、=
(A2)i= −Mかつに3(A3 )J = N
/ 2となる場合である。
第2図の時間間隔(x3)で、RFパルス(117)が
対象物に適用される。RFパルス(117)の適用の主
目的は、静磁界(Bo’)の不均一性によって位相がず
れた核スピンの位相のずれを元に戻すことである。
対象物に適用される。RFパルス(117)の適用の主
目的は、静磁界(Bo’)の不均一性によって位相がず
れた核スピンの位相のずれを元に戻すことである。
主磁界強度を測定すると、この不均一性は、0.5〜5
.0テスラの磁界を発生させる標準的な磁石と類似して
いる。すなわち、通常、不均一性は、映し出される対象
物の区域に適用された磁界の強度(SO)のIgloo
万一10/100万である。従って、絶対項における不
均一性は、強磁界においては、もっと厳しい。
.0テスラの磁界を発生させる標準的な磁石と類似して
いる。すなわち、通常、不均一性は、映し出される対象
物の区域に適用された磁界の強度(SO)のIgloo
万一10/100万である。従って、絶対項における不
均一性は、強磁界においては、もっと厳しい。
RFパルス(117)によって位相ずれを戻すには、パ
ルスを180°パルスとして形成するのが望ましい。
ルスを180°パルスとして形成するのが望ましい。
180″パルス(1,17)を使用することは、最初の
90°RFパルス(116)から生じるオリジナル肝信
号の鏡像、もしくはスピンエコーとなるので、周知であ
る。このエコー効果は、第2図に鎖線(200)で示し
ており、時定数T2 と共に増大しかつ衰退する。T2
は、不均一磁界における核スピンの相違干渉性の損失
用の時定数であり、先行技術で定義されている。
90°RFパルス(116)から生じるオリジナル肝信
号の鏡像、もしくはスピンエコーとなるので、周知であ
る。このエコー効果は、第2図に鎖線(200)で示し
ており、時定数T2 と共に増大しかつ衰退する。T2
は、不均一磁界における核スピンの相違干渉性の損失
用の時定数であり、先行技術で定義されている。
対象物に関連のある平面内で選択された核スピンの位相
のずれを直すために使用する180°RFパルス(1]
、7)に対する入力として、(B)から共役位置(c)
へ移る軌線は、k空間にもたらされる。このに空間では
、(c)は、空間周波数座標に2=(kz)に十Mと(
k3 )j = 十N/2によって、画定される。
のずれを直すために使用する180°RFパルス(1]
、7)に対する入力として、(B)から共役位置(c)
へ移る軌線は、k空間にもたらされる。このに空間では
、(c)は、空間周波数座標に2=(kz)に十Mと(
k3 )j = 十N/2によって、画定される。
標準的なRF送信コイルによって発生されるRF磁界は
、均一ではないので、RFパルス(117)の実際のパ
ルス角度は、対象物によって占有される体積に応じて変
化するおそれがある。特に、(最初の90°励起パルス
(116)を受信しなかった)選択された平面の外側の
パルス角度は、平面内に比べて不正確であってもよい。
、均一ではないので、RFパルス(117)の実際のパ
ルス角度は、対象物によって占有される体積に応じて変
化するおそれがある。特に、(最初の90°励起パルス
(116)を受信しなかった)選択された平面の外側の
パルス角度は、平面内に比べて不正確であってもよい。
選択された平面の外側の核スピンは、以前の90°パル
ス(1,16)は受信しなかったが、RFパルス(11
,7)のパルス角度が180°でない場合は、不必要な
信号を発生させてしまう。この不必要な信号を排除する
ために、勾配の等しいパルスが180°RFパルス(1
17)の両側に適用される。例えば、第2図に(132
a) (132b)で示す第3勾配によって発生した勾
配パルスを適用してもよい。
ス(1,16)は受信しなかったが、RFパルス(11
,7)のパルス角度が180°でない場合は、不必要な
信号を発生させてしまう。この不必要な信号を排除する
ために、勾配の等しいパルスが180°RFパルス(1
17)の両側に適用される。例えば、第2図に(132
a) (132b)で示す第3勾配によって発生した勾
配パルスを適用してもよい。
もし、勾配波形(132a) (132b)の下部エリ
ア(Q)(Ql)がほぼ同一であり、そのことが所望の
場合は、有効位相符号化は、180°パルスが波形の間
に反転して適用されるので、波形(132b)が完成し
た後では、関連の平面内ではスピンには発生しない。
ア(Q)(Ql)がほぼ同一であり、そのことが所望の
場合は、有効位相符号化は、180°パルスが波形の間
に反転して適用されるので、波形(132b)が完成し
た後では、関連の平面内ではスピンには発生しない。
しかし、RF磁界の不均一性の結果、RFパルス(11
7)によって選択された平面の外側に発生した信号は、
パルス(132b)によって位相がずれ、パルスシーケ
ンスの残部は、これら特別のスピンの位相のずれの戻し
が起らないように設計されている。
7)によって選択された平面の外側に発生した信号は、
パルス(132b)によって位相がずれ、パルスシーケ
ンスの残部は、これら特別のスピンの位相のずれの戻し
が起らないように設計されている。
この方法によって、十分な同一性を有し、適度に広がっ
たサンプル内の「平面外」の信号がすべて効果的に除去
される。
たサンプル内の「平面外」の信号がすべて効果的に除去
される。
代わりの方法、もしくは前記の方法以外に(114a)
の代わりに、(114b)の勾配パルスを使用すること
も可能である。この場合、スライス内の位相のずれたス
ピンは、勾配(114b)によって180’ RFパル
スの後でのみ位相のずれが戻され、パルスシーケンスも
、RFパルス(↑17)によって選択ぎわだ平面の外側
に発生する信号のすべての位相のずれを戻す。
の代わりに、(114b)の勾配パルスを使用すること
も可能である。この場合、スライス内の位相のずれたス
ピンは、勾配(114b)によって180’ RFパル
スの後でのみ位相のずれが戻され、パルスシーケンスも
、RFパルス(↑17)によって選択ぎわだ平面の外側
に発生する信号のすべての位相のずれを戻す。
本発明の別の実施例では、180°RFパルスは省略さ
れ、すべてのエコーは、勾配の応用のみによって形成さ
れている。この場合、第3図に示すように、B点から0
点へ行く軌線はないが、k空間の横移動部がB点から始
まり、正方向(K2)へ伸びている。k空間を部分的に
満たす場合(後で説明する)には、横移動部は、B′点
から始まる。k3に沿っての第1の横断は、正方向に3
を向いており、これは、180°パルスが適用される場
合と反対である。これは、第3(読出し)勾配は、反対
の極性に適用しなければならないということを意味する
。
れ、すべてのエコーは、勾配の応用のみによって形成さ
れている。この場合、第3図に示すように、B点から0
点へ行く軌線はないが、k空間の横移動部がB点から始
まり、正方向(K2)へ伸びている。k空間を部分的に
満たす場合(後で説明する)には、横移動部は、B′点
から始まる。k3に沿っての第1の横断は、正方向に3
を向いており、これは、180°パルスが適用される場
合と反対である。これは、第3(読出し)勾配は、反対
の極性に適用しなければならないということを意味する
。
また、k空間の移動は、正方向に2へ続いているので、
第2(位相符号化)勾配を、18(1°パルスを含む場
合の極性と反対の極性の中へ適用しなければならない。
第2(位相符号化)勾配を、18(1°パルスを含む場
合の極性と反対の極性の中へ適用しなければならない。
本発明の重要な特徴によると、MR影像用の完全なデー
タは、約100msec以内で得られ、関連の平面内で
選択された核の特性減衰時間T2より短い時間であるこ
とが望ましい。
タは、約100msec以内で得られ、関連の平面内で
選択された核の特性減衰時間T2より短い時間であるこ
とが望ましい。
このような短時間でのデータ取得は、心臓等の人体器官
の動きによって起る影像への悪影響をほぼ0にするため
に必須であった。このような迅速なデータ収集は、第2
(位相符号化)及び第3(読出し)勾配を映し出される
対象物に交互に適用することによって、第2図のパルス
シーケンスの時間間隔(X4)の中で行なわれる。
の動きによって起る影像への悪影響をほぼ0にするため
に必須であった。このような迅速なデータ収集は、第2
(位相符号化)及び第3(読出し)勾配を映し出される
対象物に交互に適用することによって、第2図のパルス
シーケンスの時間間隔(X4)の中で行なわれる。
本発明によると、読出し勾配波形(140)は、先行技
術文献で示されたように、台形とはなり得ない。台形の
読出し勾配を使用すると、重大な問題が発生する。とい
うのは、勾配を生じさせるコイルの電気インダクタンス
が、所望の磁界勾配の強度を達成するために、必然的に
大きくなるからである。その結果、勾配波形の立上がり
時間は、前述の時間枠の中で必要なデータを取得するの
に十分速くなければならず、実際上は、空論の動力を使
用して達成されうる。
術文献で示されたように、台形とはなり得ない。台形の
読出し勾配を使用すると、重大な問題が発生する。とい
うのは、勾配を生じさせるコイルの電気インダクタンス
が、所望の磁界勾配の強度を達成するために、必然的に
大きくなるからである。その結果、勾配波形の立上がり
時間は、前述の時間枠の中で必要なデータを取得するの
に十分速くなければならず、実際上は、空論の動力を使
用して達成されうる。
この高速変化磁界勾配を発生させるための好ましい方法
は、1986年12月9日に特許され、本発明の出願人
に譲渡された米国特許第4,628,264号明細書に
記載されており、ほぼ正弦曲線である勾配波形線分(1
32a) (132b)を生じさせる。
は、1986年12月9日に特許され、本発明の出願人
に譲渡された米国特許第4,628,264号明細書に
記載されており、ほぼ正弦曲線である勾配波形線分(1
32a) (132b)を生じさせる。
第2図において、もし、勾配波形(130)の下部エリ
ア(A、が、波形(140)の下部エリア(A4)の約
半分であれば、第3磁界勾配の方向において、核信号の
位相のずれの戻しが起こり、この方向における信号の最
大の位相のずれの戻しが、勾配波形(140)のほぼピ
ーク時に起こる。
ア(A、が、波形(140)の下部エリア(A4)の約
半分であれば、第3磁界勾配の方向において、核信号の
位相のずれの戻しが起こり、この方向における信号の最
大の位相のずれの戻しが、勾配波形(140)のほぼピ
ーク時に起こる。
信号は、同じエリアが勾配波形(140)によって掃引
された後で、各標本点が取得されるような方法で適切な
時期にサンプリングされうる。よって。
された後で、各標本点が取得されるような方法で適切な
時期にサンプリングされうる。よって。
核スピンによって、有効位相回転が等しく増加させられ
る。
る。
サンプリング用「コーム」の時間依存性が、勾配波形(
140)と完全に一致させられるこの非線形サンプリン
グ法は、先行技術に記載されている。例えば、オーディ
ッジ及びマンスフイールドによる米国特許第4,509
.0+、5号明細書「核磁気共鳴方法」を参照されたい
。
140)と完全に一致させられるこの非線形サンプリン
グ法は、先行技術に記載されている。例えば、オーディ
ッジ及びマンスフイールドによる米国特許第4,509
.0+、5号明細書「核磁気共鳴方法」を参照されたい
。
別のデータサンプリングおよびデータ処理方法も、もち
ろん可能である。その一つの方法では、データは、標本
点の間で等時間増分と共に得られ、信号の最大周波数成
分を表示するのに十分速い。
ろん可能である。その一つの方法では、データは、標本
点の間で等時間増分と共に得られ、信号の最大周波数成
分を表示するのに十分速い。
ついで、データは、k空間のステップかに3における増
加とほぼ同じとなるように、例えば補間及び再サンプリ
ングによって処理される。k3におけるデータの過サン
プリングは、すべての機器的影響を補うために、等時間
間隔及び非等時間間隔のどちらの場合も必要となるおそ
れがある。
加とほぼ同じとなるように、例えば補間及び再サンプリ
ングによって処理される。k3におけるデータの過サン
プリングは、すべての機器的影響を補うために、等時間
間隔及び非等時間間隔のどちらの場合も必要となるおそ
れがある。
更に、勾配の形状は渦電流効果のために、空間内で変化
しうるので、どの非等サンプリング間隔も対象物の回折
パターンを正確に測定するのに十分とはなり得ない。こ
の場合、空間的依存勾配形状は、渦電流効果を補うため
の情報を提供するために、前もって測定されうる。
しうるので、どの非等サンプリング間隔も対象物の回折
パターンを正確に測定するのに十分とはなり得ない。こ
の場合、空間的依存勾配形状は、渦電流効果を補うため
の情報を提供するために、前もって測定されうる。
勾配波形(140)の間に生じるスピンエコーのサンプ
リングの終わりに、k軌線は、k2、k3=(M、N/
2−1)からに2、K3=(M、N/2)、すなわち、
第3図のD点へ移動した。特に、データ取得中の=32
− 非線形サンプリング中には、k軌線は、K、線をCから
Dへ等増分の中でサンプリングされる。(第3勾配方向
における最大位相戻しは、軌線かに、=Oの線を通過す
る際に起こる)。
リングの終わりに、k軌線は、k2、k3=(M、N/
2−1)からに2、K3=(M、N/2)、すなわち、
第3図のD点へ移動した。特に、データ取得中の=32
− 非線形サンプリング中には、k軌線は、K、線をCから
Dへ等増分の中でサンプリングされる。(第3勾配方向
における最大位相戻しは、軌線かに、=Oの線を通過す
る際に起こる)。
従って、サンプリングは、ステップが八に3=トに、に
沿った軌線のポイントで行なわれる。
沿った軌線のポイントで行なわれる。
ここで、X3は、第3方向の視界の空間磁界であり、G
3は、第3勾配の強度であり、tj−ti+、は、i番
目とi十□番Hの標本点の間の間隔である。
3は、第3勾配の強度であり、tj−ti+、は、i番
目とi十□番Hの標本点の間の間隔である。
ついで、位相符号化波形(150)が第2の方向に適用
され、k軌線を第3図のE点に移動させる。
され、k軌線を第3図のE点に移動させる。
られる。ここで、X2は、第2方向の視界の空間磁界で
あり、G2は、第2勾配の強度であり、t2−t□は、
適用の間隔である。この期間は、データサンプリングは
行なわれない。
あり、G2は、第2勾配の強度であり、t2−t□は、
適用の間隔である。この期間は、データサンプリングは
行なわれない。
次の読出し勾配パルス(142)が、前の勾配波形(1
40)と同一強度だが反対の極性の中に適用されると、
第3図のに3=Oに相当する位置における第3勾配の方
向で核信号の位相戻しが行なわれる。
40)と同一強度だが反対の極性の中に適用されると、
第3図のに3=Oに相当する位置における第3勾配の方
向で核信号の位相戻しが行なわれる。
前述の方法によるデータサンプリングによって。
第3図のに空間にEからFへ平行な軌線が示される。こ
の軌線は、第2勾配方向のテップにの位相符号化によっ
てCからDへ移動する第1の軌線から変換される。
の軌線は、第2勾配方向のテップにの位相符号化によっ
てCからDへ移動する第1の軌線から変換される。
次に、第2及び第3磁界勾配の波形を繰り返し適用する
と、k空間を完全におおうに軌線ラスターが示される。
と、k空間を完全におおうに軌線ラスターが示される。
第2図示の読出し勾配(140)は、不連続である。
すなおち、位相符号化勾配の短かいパルス(150)が
適用されるときは、勾配(140)の振幅はOとなる。
適用されるときは、勾配(140)の振幅はOとなる。
第5図示の好適実施例では、読出し勾配は、連続正弦波
(160)として適用されている。ついで、位相符号化
勾配(150)が、前述のように適用されるが、読出し
勾配(160)がOを通過して交差する点の回りに集め
られる。
(160)として適用されている。ついで、位相符号化
勾配(150)が、前述のように適用されるが、読出し
勾配(160)がOを通過して交差する点の回りに集め
られる。
位相符号化勾配(1,50)の短かいパルスのために、
位相符号化及び読出し勾配が存在するときに収集される
いくつかのデータポイントは、影像の質の低下を伴わず
、フーリエ変換の市に削除されうる。
位相符号化及び読出し勾配が存在するときに収集される
いくつかのデータポイントは、影像の質の低下を伴わず
、フーリエ変換の市に削除されうる。
この実施例の利点は、読出し勾配を適用する際にギャッ
プがないので、取得時間がより短かいことである。従っ
て、高速移動体から伝わる人工的な動きは更に減少され
、取得の際の減衰度(T2)(T2*)は低下し、エコ
ーエンベロープ(200)は、最初の90°励起パルス
(11,6)の近くに集められる。
プがないので、取得時間がより短かいことである。従っ
て、高速移動体から伝わる人工的な動きは更に減少され
、取得の際の減衰度(T2)(T2*)は低下し、エコ
ーエンベロープ(200)は、最初の90°励起パルス
(11,6)の近くに集められる。
その結果、減衰T2より生じる信号損失が減少する。
あるいは、より短かい取得は、オリジナルの不連続読出
し波形の際に得られるのと同じ時間でもっと多くの波形
のに、ラインを得ることができるということを意味する
。更に、連続勾配波形に実施は、勾配波長の観点からは
、技術的にはより簡単と考えられ、渦電流補償は、より
簡単になる。すなわち、簡単に成就できる。
し波形の際に得られるのと同じ時間でもっと多くの波形
のに、ラインを得ることができるということを意味する
。更に、連続勾配波形に実施は、勾配波長の観点からは
、技術的にはより簡単と考えられ、渦電流補償は、より
簡単になる。すなわち、簡単に成就できる。
完全な(ほぼ1.00msec)取得の後で、k軌線は
、第3図のB点で終わる。軌線ラスターかに2=0を通
過すると、第2勾配方向におけるスピンの最大再集束が
起こり、軌線かに空間の原点、すなわち第3図のA点を
通過する場所では、第2及び第3勾配方向における最大
再某来が起こり、最大振幅に近づくスピンエコー信号が
発生する。
、第3図のB点で終わる。軌線ラスターかに2=0を通
過すると、第2勾配方向におけるスピンの最大再集束が
起こり、軌線かに空間の原点、すなわち第3図のA点を
通過する場所では、第2及び第3勾配方向における最大
再某来が起こり、最大振幅に近づくスピンエコー信号が
発生する。
第3図から、b3(すなわち読出し)方向のに軌線がす
べて交互に反転していることがわかる。これは、k空間
を適切におおうために、第3の勾配を反転させる必要が
あるからである。
べて交互に反転していることがわかる。これは、k空間
を適切におおうために、第3の勾配を反転させる必要が
あるからである。
従って、k空間の交互に交差するラインは、反転しなけ
ればならない。この操作は、適切なデータの時間反転に
相当する。従って、第3図示の軌線E、E+L、E+2
、・・・F−2、F−1、Fが、F、F−1、F−2・
E+2、E+]、Eとなる。この手続は、k空間の交互
の線に対して繰り返される。
ればならない。この操作は、適切なデータの時間反転に
相当する。従って、第3図示の軌線E、E+L、E+2
、・・・F−2、F−1、Fが、F、F−1、F−2・
E+2、E+]、Eとなる。この手続は、k空間の交互
の線に対して繰り返される。
更に、好ましくない器械効果によって発生するおそれの
あるデータのベースラインオフセットを補正するのが望
ましい。これは、影像シーケンスの最初のRFパルスを
適用する直前に、受信機の直角位相出力からくるデータ
をサンプリングすることによって達成しうる。
あるデータのベースラインオフセットを補正するのが望
ましい。これは、影像シーケンスの最初のRFパルスを
適用する直前に、受信機の直角位相出力からくるデータ
をサンプリングすることによって達成しうる。
データをフーリエ交換する前に、位相補正及び濾過も行
なう必要がある。これらのステップは、先行技術では公
知であり、説明の必要はない。
なう必要がある。これらのステップは、先行技術では公
知であり、説明の必要はない。
次に、k空間データの2次元フーリエ変換を行なうこと
によって、影像が発生し、フーリエ変換で対象物の回折
パターンを離散方法で測定する。
によって、影像が発生し、フーリエ変換で対象物の回折
パターンを離散方法で測定する。
影像は、奇数エコー(k空間で得られる1、3.5、・
・等)と偶数エコー(k空間で得られる2、4.6、・
・・等)のフーリエ変換の重ね合わせによりなる。
・等)と偶数エコー(k空間で得られる2、4.6、・
・・等)のフーリエ変換の重ね合わせによりなる。
従って、奇数及び偶数エコーを分離して部分影像に変換
したり、これらのエコーを組み合わせる前に変形させる
ことができる。事実、これらの処理ステップは、一定の
状況では器械的理由に対して望ましい。
したり、これらのエコーを組み合わせる前に変形させる
ことができる。事実、これらの処理ステップは、一定の
状況では器械的理由に対して望ましい。
部分影像の変形は、すべての器械的誤差を除去するため
に位相及び振幅を変化させることを伴う。
に位相及び振幅を変化させることを伴う。
一定の状況では、肝影像は、k空間有効範囲に関して不
完全な空間周波数データセットから発生されうる。通常
、(主磁石の不完全さが原因の)磁界の不均一性、突接
は効果及び渦電流が原因の研究中の対象物内の磁石の感
受率の変化及び無線周波数の磁界変化)によって発生す
るすべての位相誤差は、第2図及び第3図に関連して説
明したように、完全な空間周波数データセット(完全な
に空間の取得)を取得することによって補償されうる。
完全な空間周波数データセットから発生されうる。通常
、(主磁石の不完全さが原因の)磁界の不均一性、突接
は効果及び渦電流が原因の研究中の対象物内の磁石の感
受率の変化及び無線周波数の磁界変化)によって発生す
るすべての位相誤差は、第2図及び第3図に関連して説
明したように、完全な空間周波数データセット(完全な
に空間の取得)を取得することによって補償されうる。
しかし、もし、位相誤差が低空間周波数のみによって影
像上を変化するならば、完全取得が行なわれない場合で
も、このような位相誤差を補正することができる。
像上を変化するならば、完全取得が行なわれない場合で
も、このような位相誤差を補正することができる。
このようなに空間の部分取得のために、第2図のパルス
シーケンスは、次のように変形される。
シーケンスは、次のように変形される。
符号化期間(X2)中は、波形(]、20a)又は(1
20b)のエリアは、k軌線が、AからB□、すなわち
に2=(k2)j、<Mへ移動するうに減衰する。RF
パルス(117)を適用すると、軌線が共役点C1へ再
移動する。ついで、データ取得が以前のように行なわれ
るが、この場合は、k空間のに2次元におけるラインは
、最終エンドポイントBへ到達するめに、少ないライン
が得られる。
20b)のエリアは、k軌線が、AからB□、すなわち
に2=(k2)j、<Mへ移動するうに減衰する。RF
パルス(117)を適用すると、軌線が共役点C1へ再
移動する。ついで、データ取得が以前のように行なわれ
るが、この場合は、k空間のに2次元におけるラインは
、最終エンドポイントBへ到達するめに、少ないライン
が得られる。
k空間の完全な取得の場合のように、第2及び第3次元
における位相戻しは、k軌線が、k空間の原点、すなわ
ち第3図のA点と交差する際に最大となる。
における位相戻しは、k軌線が、k空間の原点、すなわ
ち第3図のA点と交差する際に最大となる。
位相誤差を最小にし、かつ共役対称を最大にするには、
第3図のA点に相当する第2及び第3勾配方向の位相戻
しポイントが、極めて適切な時期に鎖線(120b)で
示された180°RFパルスを適用して得られる「磁界
エコー」に相当するように、k空間軌線を配置するのが
望ましい。
第3図のA点に相当する第2及び第3勾配方向の位相戻
しポイントが、極めて適切な時期に鎖線(120b)で
示された180°RFパルスを適用して得られる「磁界
エコー」に相当するように、k空間軌線を配置するのが
望ましい。
もし、第1RFパルス(116)と第2RFパルス(1
17)を適用する間の時間をt□とすると、磁界エコー
が、180’ RFパルス(117)の後でt1時に発
生する。従って、縮小に空間を取得するには、前述の最
大信号位相戻しと磁界エコーのピークとを一致させるた
めに、180″RFパルス(117)の前の時間t□は
、適度に短縮しなければならない。
17)を適用する間の時間をt□とすると、磁界エコー
が、180’ RFパルス(117)の後でt1時に発
生する。従って、縮小に空間を取得するには、前述の最
大信号位相戻しと磁界エコーのピークとを一致させるた
めに、180″RFパルス(117)の前の時間t□は
、適度に短縮しなければならない。
従って、k空間の部分取得は、影像の位相変化を測定す
るのに十分なに空間部分のサンプリングを行ないつつ、
k空間の正もしくは負の部分の符号化ステップをいくつ
か省略することによって、達成される。あるいは、省略
されるに空間部分は、直角方向(k3)でもよいが、影
像の位相変化を測定するのに十分な対称的に空間のサン
プリングを行なう。
るのに十分なに空間部分のサンプリングを行ないつつ、
k空間の正もしくは負の部分の符号化ステップをいくつ
か省略することによって、達成される。あるいは、省略
されるに空間部分は、直角方向(k3)でもよいが、影
像の位相変化を測定するのに十分な対称的に空間のサン
プリングを行なう。
k空間の全データより少ないデータを取得する際に節約
される時間は、 (1)影像の解像度を向上させるために、k空間のデー
タの追加の非対称ラインを取得するため、もしくは、 (2)動きに関する誤差もしくは減衰時間t2に起因す
る信号損失を減少させるために、全取得時間をできるだ
け短かくするために使用しうる。
される時間は、 (1)影像の解像度を向上させるために、k空間のデー
タの追加の非対称ラインを取得するため、もしくは、 (2)動きに関する誤差もしくは減衰時間t2に起因す
る信号損失を減少させるために、全取得時間をできるだ
け短かくするために使用しうる。
つまり、後述の理由のために、X方向のに空間の完全な
データセットより少ないデータセットを収集するのが望
ましい。
データセットより少ないデータセットを収集するのが望
ましい。
(1)収集されるにχラインの数が減少する。よって、
データ取得時間が短縮される。その結果、動きのある対
象物の表示がもっと信頼できるようになる。しかし、よ
り少ないデータが得られるということは、最終影像内の
信号対雑音弁(SNR)が減少するということを意味す
る。しかし、これは、以下のことで補償される。
データ取得時間が短縮される。その結果、動きのある対
象物の表示がもっと信頼できるようになる。しかし、よ
り少ないデータが得られるということは、最終影像内の
信号対雑音弁(SNR)が減少するということを意味す
る。しかし、これは、以下のことで補償される。
4O−
(2)本発明の高速影像技法におけるデータは、180
°再集束パルスの後で、スピンエコーのエンベロープ内
で取得されるので、スピンエコーのエンベロープの中心
と最初の90°励起RFパルスの間の時間は、2t□で
あり、tlは、最初の90°パルスと180″再集束パ
ルスの間の時間である。スピンエコーのエンベロープ内
の信号の最大振幅は、T2緩和の結果、減衰し、その減
衰は、先□の指数関数である。従って、最良のSNRは
、tlの最小値より得られる。
°再集束パルスの後で、スピンエコーのエンベロープ内
で取得されるので、スピンエコーのエンベロープの中心
と最初の90°励起RFパルスの間の時間は、2t□で
あり、tlは、最初の90°パルスと180″再集束パ
ルスの間の時間である。スピンエコーのエンベロープ内
の信号の最大振幅は、T2緩和の結果、減衰し、その減
衰は、先□の指数関数である。従って、最良のSNRは
、tlの最小値より得られる。
完全なに空間取得の際に、エコーエンベロープのピーク
が、k空間の中心線(kx=0)と適切な時間で一致す
ると、エコーエンベロープの中心の前方にMラインが得
られ、M−1ラインはその後方に得られる。もし、Mラ
インが、取得するのにtMmsecかかると、tlはt
Mより短かいはずはない。しかし、Mラインより少ない
ラインが、エコーエンベロープの前方に得られるなら、
tM及びt□は、短縮されうる。
が、k空間の中心線(kx=0)と適切な時間で一致す
ると、エコーエンベロープの中心の前方にMラインが得
られ、M−1ラインはその後方に得られる。もし、Mラ
インが、取得するのにtMmsecかかると、tlはt
Mより短かいはずはない。しかし、Mラインより少ない
ラインが、エコーエンベロープの前方に得られるなら、
tM及びt□は、短縮されうる。
最初の90°パルスとスピンエコーエンベ口−プの間の
時間は、toの2倍であるので、部分的に空間取得の結
果得られる時間のかからないゲインは、完全な高速実験
によって効果的に2倍にされる。
時間は、toの2倍であるので、部分的に空間取得の結
果得られる時間のかからないゲインは、完全な高速実験
によって効果的に2倍にされる。
また、スピンエコーエンベロープのピークは、早めに達
成するので、t2より小さい信号減衰が起こり、ついで
SNRが高くなる。
成するので、t2より小さい信号減衰が起こり、ついで
SNRが高くなる。
部分に空間取得を行なう方法は、前述のように、90°
−180’ RFパルスの間の前符号化の期間の中の符
号化勾配(120) (130)の変形を含む。
−180’ RFパルスの間の前符号化の期間の中の符
号化勾配(120) (130)の変形を含む。
本発明の部分に空間影像の再現は、後述の理論と計算で
行なわれる。
行なわれる。
□Lゴ如遡の緩変似く本ゑ影像の再現1−影像の空間分
布は、複関数F (X)で説明されうる。−次元分布に
ついて考慮されたい(2次元の場合にも一般性があるこ
とは明白である)。密度は、F (X)の大きさで与え
られる。
布は、複関数F (X)で説明されうる。−次元分布に
ついて考慮されたい(2次元の場合にも一般性があるこ
とは明白である)。密度は、F (X)の大きさで与え
られる。
空間分布関数のフーリエ変換
Fn=ΣF (X) e−”′y″ (
1)は、2つのDFT式NMRM像実験で得られた時間
データを表わす。(1)式において、X−、(m ”
Oy・・・、iリ−1)であり、Nは、絵素の数であり
、Kn=2ル ーである。Lは、結像エリアの長さである。
1)は、2つのDFT式NMRM像実験で得られた時間
データを表わす。(1)式において、X−、(m ”
Oy・・・、iリ−1)であり、Nは、絵素の数であり
、Kn=2ル ーである。Lは、結像エリアの長さである。
(T)式から、関数Fnは、Nの期間を有するnの周期
関数であることがわかる。よって、実験で得られた時間
データの範囲外でFnを定義することができる。
関数であることがわかる。よって、実験で得られた時間
データの範囲外でFnを定義することができる。
Fnの逆フーリエ変換によって、空間分布関数Fもし、
F (X)が実数ならば、F(X)は、k空間(フーリ
エ領域)の半分の領域内の時間データを使用もし、F
(X)は、実数ではないが、その位相は空間における緩
変化関数ならば、f (X)は、k空間の半分のデータ
以外に、フーリエ調和関数を使用して、Fnから再現す
ることができる。
F (X)が実数ならば、F(X)は、k空間(フーリ
エ領域)の半分の領域内の時間データを使用もし、F
(X)は、実数ではないが、その位相は空間における緩
変化関数ならば、f (X)は、k空間の半分のデータ
以外に、フーリエ調和関数を使用して、Fnから再現す
ることができる。
F (X)の大きさG (X)は、次のように表わされ
る。
る。
G(X) = ’? (x)F (X、)
(4/フーリエ領域では、この方程式は、
次のようになる。
(4/フーリエ領域では、この方程式は、
次のようになる。
Gn= Σ♀ F
(ら)n/ rl’ h −
n’ (x)変動の長さ1が、絵素サイズL/Nよりはるかに
大きいと仮定すると、lnl≦Tである(Pn#の唯一
の低空間周波数の調和関数は、(5)式においては、不
可欠であり、そこでこの合計は、Gn= Σ
’f’F、
(6)In’ l <hn n’
h−hまで減少されうる。
(ら)n/ rl’ h −
n’ (x)変動の長さ1が、絵素サイズL/Nよりはるかに
大きいと仮定すると、lnl≦Tである(Pn#の唯一
の低空間周波数の調和関数は、(5)式においては、不
可欠であり、そこでこの合計は、Gn= Σ
’f’F、
(6)In’ l <hn n’
h−hまで減少されうる。
ここで、丁ζ乙n << Nである。
G(X)は、解像度による実関数である。G(X)を再
生するにはO’n≦トシタの範囲のGnを知るだけでよ
い。これらのGnは、Fnが−nun≦NA + 4n
の範囲にあり、千〇が−Δn≦n≦4nの範囲にあるこ
とがわかれば、(6)式から得られる。所望のFnは、
実験から得なければならない。次に、どのようにして、
取得したFnから?nを取得するかを示す。
生するにはO’n≦トシタの範囲のGnを知るだけでよ
い。これらのGnは、Fnが−nun≦NA + 4n
の範囲にあり、千〇が−Δn≦n≦4nの範囲にあるこ
とがわかれば、(6)式から得られる。所望のFnは、
実験から得なければならない。次に、どのようにして、
取得したFnから?nを取得するかを示す。
最初に、関数を導入する。
14X)−一Σ F ei lc 、XN−、。っ
、54□ ・ (I)
ここで、 5(X) =Sin[”(6n+’/2 )xコ/s’
+n X (g)し であり、関数F (X)は、X=0において、ピークS
o= 2.an十lを有する。ピークS (X)の外側
は、はるかに小さい装置能力を有する。もし、長さ□2
乙rN−1 の重要な周期構造が、影像源に存在しなれければ、(7
)式のX′を超える合計に対する寄与は、×−△×から
X+△Xの範囲から発生する。
、54□ ・ (I)
ここで、 5(X) =Sin[”(6n+’/2 )xコ/s’
+n X (g)し であり、関数F (X)は、X=0において、ピークS
o= 2.an十lを有する。ピークS (X)の外側
は、はるかに小さい装置能力を有する。もし、長さ□2
乙rN−1 の重要な周期構造が、影像源に存在しなれければ、(7
)式のX′を超える合計に対する寄与は、×−△×から
X+△Xの範囲から発生する。
ここで、T−’(X’)の変数X′をXと変換する。
’P (X)は、(9)式から簡単に得られる。
宇(X)=に(×ゾ1白×)1 (り方
程式(6)(7)(10)によって、緩変化位相を有す
る影像用の影像再現の問題が解決される。この技法を使
用すると、光干渉効果が観測されることがある。これは
、影像が長さC/、、、、オ、((8)図参照)の重要
な周期構造を有するときに、発生しうる。
程式(6)(7)(10)によって、緩変化位相を有す
る影像用の影像再現の問題が解決される。この技法を使
用すると、光干渉効果が観測されることがある。これは
、影像が長さC/、、、、オ、((8)図参照)の重要
な周期構造を有するときに、発生しうる。
2、再現計算
次に、(6)(7)(10)式に基づく影像再現の計算
を2つ説明する。取得データはに空間の半分と△n低周
波数ラインをふさぎ、k空間の残りの部分は、Oになる
とみなす6(6)式を使用するには、△n高周波数ライ
ンを知ることも必要である。しかし、はとんどの場合、
高周波数ラインは、通常、非常に小さいので、無視する
ことができる。エコー内の符号化は、水平方向に相当す
るとみなされる。
を2つ説明する。取得データはに空間の半分と△n低周
波数ラインをふさぎ、k空間の残りの部分は、Oになる
とみなす6(6)式を使用するには、△n高周波数ライ
ンを知ることも必要である。しかし、はとんどの場合、
高周波数ラインは、通常、非常に小さいので、無視する
ことができる。エコー内の符号化は、水平方向に相当す
るとみなされる。
術−算□L
8.1水平方向のRFT(逆フーリエ変換)S、22△
n+1低周波数ラインを使用して影像を作り出す。
n+1低周波数ラインを使用して影像を作り出す。
S、38.2用の位相緩変化数甲を見出す。
8.4垂直方向におけるS、3のDFT (垂直フーリ
エ変換) 5.5(6)式を使用してGn(n =O,N/2)を
計算する。
エ変換) 5.5(6)式を使用してGn(n =O,N/2)を
計算する。
におけるRFT。
S、73.6の実部を取り出す。
L」■非
3.1水平方向のRFT
S、22Δn+1低周波数ラインを使用して影像を作り
出す。
出す。
S、33.2用の位相関数を見出す。
3.4時間データのすべてのラインを使用して影像を作
り出す。
り出す。
S、53.4に3.3を乗じる。
S、68.5のDFT
5.7 ラインn=1、n/2用の3.6の値を2倍に
し、ラインn =N/2+1、N−1を0にする。
し、ラインn =N/2+1、N−1を0にする。
S、8 S、7のRFT。
S、95.8の実部を取り出す。
計算■と計算■は、関数的に完全に同等である。
主に、X方向の位相符号化ラインの100%未満の取得
に関して部分に空間取得を説明したが、符号化前のGy
勾配のエリアを減少させることによって、k空間のy次
元の100%未満の取得も、同様に可能である。k空間
のX方向における走査は、約300μsec L/かか
からないので、明らかなT2減衰は起こらない。
に関して部分に空間取得を説明したが、符号化前のGy
勾配のエリアを減少させることによって、k空間のy次
元の100%未満の取得も、同様に可能である。k空間
のX方向における走査は、約300μsec L/かか
からないので、明らかなT2減衰は起こらない。
第6図示の原点の回りのky定走査移動させ、か47一
つ部分ky定走査おけるポイントと同数のポイントを収
集し、それによって信号帯域幅を同一に維持することに
より、高空間周波数を描くことができる。従って、この
変量によって、「ワン・ショット」で低SNRの高解像
度影像が産み出される。
集し、それによって信号帯域幅を同一に維持することに
より、高空間周波数を描くことができる。従って、この
変量によって、「ワン・ショット」で低SNRの高解像
度影像が産み出される。
モザイク走渣−訳
身体を静止させることができ、かつ「ワン・ショット」
以上取り出すのに十分な時間が得られるような状況では
、高周波数は、k空間を2.4、もしくはそれ以上の部
分に分割することによって描くことができ、各部分は、
部分に空間取得でおおわれる。
以上取り出すのに十分な時間が得られるような状況では
、高周波数は、k空間を2.4、もしくはそれ以上の部
分に分割することによって描くことができ、各部分は、
部分に空間取得でおおわれる。
第7図に、この考えをに空間の「2−タイル」モザイク
走査で示す。各単一取得部でおおわれた全エリアは、ほ
ぼ一定である(Kxのやや大きい有効範囲が第7図に示
されている)が、第2図示のパルスシーケンスの906
RFパルス(1,16)と180’ RFパルス(11
7)の間の適切な前符号化によって置換される。
走査で示す。各単一取得部でおおわれた全エリアは、ほ
ぼ一定である(Kxのやや大きい有効範囲が第7図に示
されている)が、第2図示のパルスシーケンスの906
RFパルス(1,16)と180’ RFパルス(11
7)の間の適切な前符号化によって置換される。
第8A図乃至第8C図に示すように、2つの取得部は、
所望のエリアをおおうのに十分であり、解像度を向上さ
せる。取得されたに空間部は、「モザイク」を形成する
。このモザイクは、フーリエ変換の前に継ぎ合わされる
。従って、合成影像は、本発明の単一の高速取得で得ら
れるものより解像度が高い(「瞬間走査」と称される)
。第8C図に示すように「タイル」の間のオーバーラツ
プは、k空間の間の部分を正確に結合するために備える
場合には有効であるが、必らずしも必要ではない。
所望のエリアをおおうのに十分であり、解像度を向上さ
せる。取得されたに空間部は、「モザイク」を形成する
。このモザイクは、フーリエ変換の前に継ぎ合わされる
。従って、合成影像は、本発明の単一の高速取得で得ら
れるものより解像度が高い(「瞬間走査」と称される)
。第8C図に示すように「タイル」の間のオーバーラツ
プは、k空間の間の部分を正確に結合するために備える
場合には有効であるが、必らずしも必要ではない。
k空間有効範囲を真っすぐ延長するのが可能なのは、k
y定走査312μsee L/かかからないからである
。このような短時間では、T2*減衰は無視できるので
、伝達関数は、本質的に一様である。
y定走査312μsee L/かかからないからである
。このような短時間では、T2*減衰は無視できるので
、伝達関数は、本質的に一様である。
第9図は、2つの連続走査用のパルスシーケンス及び信
号を示しており、Gxがない場合に示される第]−の走
査、すなわち走査1は、適切な前符号化によって得られ
る変換エコーを示す。
号を示しており、Gxがない場合に示される第]−の走
査、すなわち走査1は、適切な前符号化によって得られ
る変換エコーを示す。
走査2には、2つの選択がある。第1は、走査1のよう
に、各kxラインに関しては同一方向でkyラインを走
査することであり、第2は、2−タイルの場合のために
、反対方向で走査することであり、各走査の開始に関し
て同時間に得られるky=Oとなる。何れの場合も、k
yラインの各線分は、常に勾配パルスの中心から±15
6μsec以内で変換され、多数の「タイル」を接着さ
せ、各線分は、影像内の高空間周波数を示す。
に、各kxラインに関しては同一方向でkyラインを走
査することであり、第2は、2−タイルの場合のために
、反対方向で走査することであり、各走査の開始に関し
て同時間に得られるky=Oとなる。何れの場合も、k
yラインの各線分は、常に勾配パルスの中心から±15
6μsec以内で変換され、多数の「タイル」を接着さ
せ、各線分は、影像内の高空間周波数を示す。
前述のように、1回の取得でkx及びky内で得られる
ポイント数は、それぞれ54 X 128である。X及
びY方向において、2に近づく係数(絵素エリアでの4
倍の減少を意味する)は、4−タイルモザイク走査によ
って達成される。この走査では、k空間の各象限からの
に空間取得は、第10図に示すにように接着される。
ポイント数は、それぞれ54 X 128である。X及
びY方向において、2に近づく係数(絵素エリアでの4
倍の減少を意味する)は、4−タイルモザイク走査によ
って達成される。この走査では、k空間の各象限からの
に空間取得は、第10図に示すにように接着される。
M E S H走査法
モザイク走査の特徴は、各成分取得は、それに関連する
非対称のT2減衰を有し、取得部がモザイクで作られる
場合は、このT、減衰は、不連続関数を形成する。この
減衰は、kx次元において最も顕著であり、このkx次
元では、取得は実時間で起る。
非対称のT2減衰を有し、取得部がモザイクで作られる
場合は、このT、減衰は、不連続関数を形成する。この
減衰は、kx次元において最も顕著であり、このkx次
元では、取得は実時間で起る。
この問題を解決するには、もし、T2減衰が比較的短か
く、かつ磁界が、かなり不均一ならば、第10図に示す
ように4−タイルモザイク走査を使用し、X方向におけ
るに空間の反対部分のTacqの間に適用されたGx勾
配の方向を反対にすることである。
く、かつ磁界が、かなり不均一ならば、第10図に示す
ように4−タイルモザイク走査を使用し、X方向におけ
るに空間の反対部分のTacqの間に適用されたGx勾
配の方向を反対にすることである。
反対の勾配、すなわち最初に+kx、次に−kxに沿っ
て走査すると、X方向におけるT2減衰は、kx=oの
回りでは対称的であることがわかる。
て走査すると、X方向におけるT2減衰は、kx=oの
回りでは対称的であることがわかる。
最初の軌線A−B−Cは、kx=oの近くにある。
180°パルスによって発生する磁界エコーは、+kx
/2付近に現われる。従って、4−タイルモザイク走査
は、高均一性の状況に対して望ましい。
/2付近に現われる。従って、4−タイルモザイク走査
は、高均一性の状況に対して望ましい。
あるいは、T2が比較的長く、磁界がかなり不均一な場
合は、kx次元を、第10図に示すように、ky次元が
MESII(Mosaicked−Echo 5can
1lybrjd、モザイク式エコー走査ハイブリッド
)と呼ばれるモザイク化によっておおわれる際に、次の
取得部からのラインを間に差し込む、すなわち、「かみ
合わせる」ことによっておおうことができる。kx次元
の間に差し込むと、T2減衰は、確実にkxの連続関数
になる。
合は、kx次元を、第10図に示すように、ky次元が
MESII(Mosaicked−Echo 5can
1lybrjd、モザイク式エコー走査ハイブリッド
)と呼ばれるモザイク化によっておおわれる際に、次の
取得部からのラインを間に差し込む、すなわち、「かみ
合わせる」ことによっておおうことができる。kx次元
の間に差し込むと、T2減衰は、確実にkxの連続関数
になる。
第1図示のように、有効範囲が2倍になる場合は、最初
の軌線A−B□−G□は、k空間の端まで進む。パルス
勾配Gxを倍加することによって、ステップDEが倍加
され、所望のkxエリア全体が、半分の強度で1回の走
査でおおわれる。
の軌線A−B□−G□は、k空間の端まで進む。パルス
勾配Gxを倍加することによって、ステップDEが倍加
され、所望のkxエリア全体が、半分の強度で1回の走
査でおおわれる。
kxに沿った次の走査A−82−C2はわずかなオフセ
ットと共に始まり、前の有効範囲を差し込む、すなわち
、「かみ合わせる」。この場合には、T2減衰は、kx
Xoのまわりに集められ、かつ対称的な磁界エコーを有
する全に空間取得におけるのと同等に示されている。
ットと共に始まり、前の有効範囲を差し込む、すなわち
、「かみ合わせる」。この場合には、T2減衰は、kx
Xoのまわりに集められ、かつ対称的な磁界エコーを有
する全に空間取得におけるのと同等に示されている。
わずかな非対称性が短かいT2の間に生じるおそれがあ
るが、実際には、これは磁界エコーであり、特に高電界
においては、重大関心事であり、対称的にする必要があ
る。MESHを行なうにはより強いGx勾配を使用する
必要があるが、実際にはこれは問題ではない。というの
は、この勾配は、常にy勾配より弱いからである。また
、MESHを使用することによって、kxの有効範囲を
増加させるために、2より大きい係数を得ることができ
る。
るが、実際には、これは磁界エコーであり、特に高電界
においては、重大関心事であり、対称的にする必要があ
る。MESHを行なうにはより強いGx勾配を使用する
必要があるが、実際にはこれは問題ではない。というの
は、この勾配は、常にy勾配より弱いからである。また
、MESHを使用することによって、kxの有効範囲を
増加させるために、2より大きい係数を得ることができ
る。
モザイク及びMESH法においては、信号の帯域幅は一
定であるので、走査の平均化効果によって、SNRは、
絵素エリアの減少分の平方根に比例して確実に減少する
。従って、シングルショット影像と同等のSNRを達成
するには、平均化の適切な量が必要なだけである。例え
ば、4−タイルモザイクもしくはMESHの場合には、
単一の低解像度取得のSNRへ再び到達するために、全
部で4つの平均値もしくは、16の走査を行なわなけれ
ばならない。
定であるので、走査の平均化効果によって、SNRは、
絵素エリアの減少分の平方根に比例して確実に減少する
。従って、シングルショット影像と同等のSNRを達成
するには、平均化の適切な量が必要なだけである。例え
ば、4−タイルモザイクもしくはMESHの場合には、
単一の低解像度取得のSNRへ再び到達するために、全
部で4つの平均値もしくは、16の走査を行なわなけれ
ばならない。
モザイク走査、及びMESHにおけるマルチ−スライス
影像法を使用すると、複数の平面が続々と得られ、3次
元影像に備えうる。
影像法を使用すると、複数の平面が続々と得られ、3次
元影像に備えうる。
モザイク走査及びMESHを使用する直接斜面走査は、
前述のように、非直角勾配によって磁気共鳴研究を行な
うことによって達成される。
前述のように、非直角勾配によって磁気共鳴研究を行な
うことによって達成される。
化学シフト問題を解決するには、次の2つの方法のいず
れかが行なわれる。
れかが行なわれる。
第1の方法では、第4図における化学成分(水分もしく
は、脂肪)の共鳴周波数と同等の搬送周波数及びT2*
に起因する広がりをおおうのに十分な帯域幅を有する9
0°RF制御パルス(180)が、第1(スライス−セ
レクション)勾配を適用する前に、適用される。
は、脂肪)の共鳴周波数と同等の搬送周波数及びT2*
に起因する広がりをおおうのに十分な帯域幅を有する9
0°RF制御パルス(180)が、第1(スライス−セ
レクション)勾配を適用する前に、適用される。
十分な強度の勾配パルスによって、選択された成分が、
k空間の関連部分から完全に除去され、第2図のパルス
シーケンスが適用されると読出さなくなるように、選択
された成分の位相応答かに空間の中でシフトされる。こ
の種の抑制シーケンスは、この技術分野においては公知
であり、これ以上の説明は必要としないであろう。
k空間の関連部分から完全に除去され、第2図のパルス
シーケンスが適用されると読出さなくなるように、選択
された成分の位相応答かに空間の中でシフトされる。こ
の種の抑制シーケンスは、この技術分野においては公知
であり、これ以上の説明は必要としないであろう。
化学シフトを解決するための第2の方法は、180°位
相戻しパルスの周波数スペクトルを、所望の化学成分の
共鳴周波数をカバーするのに丁度良くなるように減少さ
せることである。従って、望まれていない成分は、位相
戻しは行なわれない(すなわち、第3図のに空間におい
ては、対角線的にはシフトされない)し、サンプリング
中は、k空間の関連部分から移動され、(読出しは行な
わない)。
相戻しパルスの周波数スペクトルを、所望の化学成分の
共鳴周波数をカバーするのに丁度良くなるように減少さ
せることである。従って、望まれていない成分は、位相
戻しは行なわれない(すなわち、第3図のに空間におい
ては、対角線的にはシフトされない)し、サンプリング
中は、k空間の関連部分から移動され、(読出しは行な
わない)。
このような方法は、コンピュータ断層撮影法認(Jou
rnal of Computer As5isted
Tomography)。
rnal of Computer As5isted
Tomography)。
9 (4):651−658ページ、7月/8月、ロー
エン・プレス、ニューヨーク(1985)のピータ−・
エム・ジョセフ(Peter M、Joseph)によ
る[スピンエコー化学シフI−MR影像技法](”A
5pin Echo ChemicalShift M
RImaging Technique”)という論文
の中に開示されている。
エン・プレス、ニューヨーク(1985)のピータ−・
エム・ジョセフ(Peter M、Joseph)によ
る[スピンエコー化学シフI−MR影像技法](”A
5pin Echo ChemicalShift M
RImaging Technique”)という論文
の中に開示されている。
本発明を、複数の好適実施例に基づいて説明したが、別
の変更及び変形も、当業者には明白となるであろう。従
って、本発明は、特定の開示ではなく、特許請求の範囲
のみに限定されるものである。
の変更及び変形も、当業者には明白となるであろう。従
って、本発明は、特定の開示ではなく、特許請求の範囲
のみに限定されるものである。
第1図は、磁界勾配を通常の空間座標系x、y、Zで示
し、影像となる対象物が、磁石の内部に設定される様子
を示す概略図である。 第2図は、本発明による好適なパルスシーケンスを示す
ものである。 第3図は、観察されている第1図示の対象物の像の空間
的もしくはフーリエ定義域周波数を示すグラフである。 第4図は、化学シフトを排除するため、スライス選択前
に与えられる抑制パルスシーケンスを示す図である。 第5図は、第2図示の取得パルスシーケンスを拡大した
ものであるが、読出し勾配は、連続正弦波として与えら
れる。 第6図は、部分ky定走査おける空間周波数定義域の有
効範囲を示すものである。 第7図は、最初の軌線、および2−タイルモザイク走査
における空間周波数の有効範囲を示すものである。 第8A図乃至第8C図は、2−タイルモザイク走査にお
ける空間周波数定義域のシーケンスの有効範囲を示すも
のである。 第9図は、第7図および第8A図乃至第8c図示の2−
タイルモザイク走査に対するパルスシーケンスおよび信
号を示すものである。 第10図は、最初の軌線、および4−タイルモザイク走
査 イク走査における空間周波数定義域の有効範囲を示すも
のである。 第11図は、最初の軌線、および4−パルスMESH走
査における空間周波数定義域の有効範囲を示すものであ
る。 (10)部位 (12)被写体(1
4)磁石 (16)コイル(B、)
磁界 (x□) (X2 ) (X3
) (X、 )時間間隔(114) (114a)
(114b)波形 (116)高周波(120a)
(120b) (130) (140)波形 (A2
)(A3)エリア(117)RFパルス
(132a) (132b)波形(Q)(Q□)エリア
(150)波形(1,42)パルス
(160)正弦波(180)抑制パルス
(200)エコーエンベ口−プ手続補正帯動
式) 昭和63年3月7日 特許庁長官 小 川 邦 夫 殿 1、事件の表示 昭和62年特許願第306827号 2、発明の名称 核磁気共鳴信号を用い、対象物から高速度で影像情報を
得る方法 3、補正をする者 事件との関係 特許出願人 名 称 アドバンスト エヌエムアールシステム
ズ インコーホレーテッド 4、代理人 (発送日 昭和63年2月23日) 6、補正の対象 (1)願書の特許出願人の欄(2
)委任状及びその訳文 (3)明細書 (4)図面 7、補正の内容 (1) (2)別紙の通り。 (3)願書に最初に添付した明細書の浄書・別紙の通り
(内容に変更なし)
し、影像となる対象物が、磁石の内部に設定される様子
を示す概略図である。 第2図は、本発明による好適なパルスシーケンスを示す
ものである。 第3図は、観察されている第1図示の対象物の像の空間
的もしくはフーリエ定義域周波数を示すグラフである。 第4図は、化学シフトを排除するため、スライス選択前
に与えられる抑制パルスシーケンスを示す図である。 第5図は、第2図示の取得パルスシーケンスを拡大した
ものであるが、読出し勾配は、連続正弦波として与えら
れる。 第6図は、部分ky定走査おける空間周波数定義域の有
効範囲を示すものである。 第7図は、最初の軌線、および2−タイルモザイク走査
における空間周波数の有効範囲を示すものである。 第8A図乃至第8C図は、2−タイルモザイク走査にお
ける空間周波数定義域のシーケンスの有効範囲を示すも
のである。 第9図は、第7図および第8A図乃至第8c図示の2−
タイルモザイク走査に対するパルスシーケンスおよび信
号を示すものである。 第10図は、最初の軌線、および4−タイルモザイク走
査 イク走査における空間周波数定義域の有効範囲を示すも
のである。 第11図は、最初の軌線、および4−パルスMESH走
査における空間周波数定義域の有効範囲を示すものであ
る。 (10)部位 (12)被写体(1
4)磁石 (16)コイル(B、)
磁界 (x□) (X2 ) (X3
) (X、 )時間間隔(114) (114a)
(114b)波形 (116)高周波(120a)
(120b) (130) (140)波形 (A2
)(A3)エリア(117)RFパルス
(132a) (132b)波形(Q)(Q□)エリア
(150)波形(1,42)パルス
(160)正弦波(180)抑制パルス
(200)エコーエンベ口−プ手続補正帯動
式) 昭和63年3月7日 特許庁長官 小 川 邦 夫 殿 1、事件の表示 昭和62年特許願第306827号 2、発明の名称 核磁気共鳴信号を用い、対象物から高速度で影像情報を
得る方法 3、補正をする者 事件との関係 特許出願人 名 称 アドバンスト エヌエムアールシステム
ズ インコーホレーテッド 4、代理人 (発送日 昭和63年2月23日) 6、補正の対象 (1)願書の特許出願人の欄(2
)委任状及びその訳文 (3)明細書 (4)図面 7、補正の内容 (1) (2)別紙の通り。 (3)願書に最初に添付した明細書の浄書・別紙の通り
(内容に変更なし)
Claims (19)
- (1)核磁気共鳴信号を用い、対象物から高速度で影像
情報を得る方法において、 (a)ある軸に沿い約0.5〜5テスラの強さを有する
連続静電場に、対象物を設定する段階と、(b)対象物
に対し、スライス選択勾配を含む選択面に垂直な第1磁
界勾配とともに第1のRFパルスを与えることにより、
対象物の選択面における核スピンを励起し、そのスピン
によって、自由誘導崩壊信号を生成させる段階と、 (c)対象物に対し、前記面と平行な方向を有する第2
の磁界勾配と、同じく前記面に平行であるとともに前記
第2磁界勾配に垂直である方向を有する第3の磁界勾配
とを含む符号化シーケンスを与える段階と、 (d)前記段階(b)において与えられるスライス選択
勾配の全デフェージング効果が、可能な限りゼロに近く
なるよう、対象物に対し、更にスライス選択勾配を与え
る段階と、 (e)位相符号化勾配を決める第2磁界勾配と、読出し
勾配を決める第3磁界勾配とを、対象物に対して時間ご
とに交互に加えるようになっているシーケンスを与え、
かつ第2および第3の磁界勾配の前記シーケンスを、約
100ms以下の時間内に完了させる段階と、 (f)磁気共鳴影像をつくるため、時間定義域における
一連のデータ値をつくるべく読出し勾配が与えられてい
る間、対象物からの核磁気共鳴信号出力をサンプリング
することによってデータを取得し、前記時間定義域デー
タを、対象物の空間周波数を推定判断しうる変更データ
にフォーマット化し、かつ前記変更データを、対象物を
像として表示させるための空間定義域データに変換させ
る段階 とからなることを特徴とする核磁気共鳴信号を用い、対
象物から高速度で影像情報を得る方法。 - (2)前記段階(c)及び(d)の間に、対象物に対し
、第2のRFパルスを与える段階が入ることを特徴とす
る特許請求の範囲第(1)項に記載の核磁気共鳴信号を
用い、対象物かから高速度で影像情報を得る方法。 - (3)周波数成分の狹帯域のみを含むように、第2パル
スの長さを選択し、前記第2RFパルスを与えることに
より、化学的にシフトされた2つの成分のうち一つだけ
をリフェーズさせることを特徴とする特許請求の範囲第
(2)項に記載の核磁気共鳴信号を用い、対象物から高
速度で影像情報を得る方法。 - (4)前記段階(d)において加えられるスライス選択
勾配を、むしろ、反対の極性を有する第2RFパルスの
前に加えることを特徴とする特許請求の範囲第(2)項
に記載の核磁気共鳴信号を用い、対象物から高速度で影
像情報を得る方法。 - (5)前記段階(b)の前に、化学シフト抑制パルスシ
ーケンスを与える段階からなる抑制時間を更に含み、前
記シーケンス付与段階は、 RFパルスのバンド幅がシフトされた成分を実質的に含
み、かつ前記パルスの中心周波数が抑制されるべきライ
ンと同じ周波数になっているあらゆる磁界勾配を欠いた
状態で、核スピンにRFパルスを与える段階と、 第1、第2および第3の勾配のどれよりも大きいデフェ
ージング効果を有する磁気勾配を与える段階と、 スライス選択勾配とともに、別のRFパルスを与えるこ
とにより、前記段階(b)を繰り返し行なう段階と、 デフェージング効果を実質的に排除するべく、更にスラ
イス選択勾配を与えることにより、前記段階(b)を繰
り返し行なう段階 とからなることを特徴とする特許請求の範囲第(1)項
に記載の核磁気共鳴信号を用い、対象物から高速度で影
像情報を得る方法。 - (6)第2のRFパルスにおける不完全さにより間違っ
て励起されたあらゆるスピンの後続信号を破壊するため
、第2のRFパルスが、その前後に与えられる等しい大
きさの磁界勾配を有することを特徴とする特許請求の範
囲第(2)項に記載の核磁気共鳴信号を用い、対象物か
ら高速で影像情報を得る方法。 - (7)第2のRFパルスが、空間的に選択されるパルス
であり、かつあらゆる方向の磁界勾配とともに加えられ
ることを特徴とする特許請求の範囲第(2)項に記載の
核磁気共鳴信号を用い、対象物から高速度で影像情報を
得る方法。 - (8)第2のRFパルスが、全影像面の未変更周波数ス
ペクトルよりも広い周波数スペクトルを有する広帯域パ
ルスであることを特徴とする特許請求の範囲第(2)項
に記載の核磁気共鳴信号を用い、対象物から高速度で影
像情報を得る方法。 - (9)位相符号化勾配および読出し勾配が、適切なタイ
ミンで部分的に重畳して与えられ、前記読出し勾配は、
ほぼ正弦状の連続波として与えられるようになっており
、また、前記読出し勾配および位相符号化勾配が重畳さ
れている間に取得されるデータが、変換を行なう前に消
去されることを特徴とする特許請求の範囲第(1)項に
記載の核磁気共鳴信号を用い、対象物から高速度で影像
情報を得る方法。 - (10)第1のRFパルスを与える前に、レシーバベー
スラインレベルの複雑なサンプルをとり、それを、順次
得られるデータのベースラインオフセットを補正するた
めに使用する段階を含むことを特徴とする特許請求の範
囲第(1)項に記載の核磁気共鳴信号を用い、対象物か
ら高速度で影像情報を得る方法。 - (11)第1のRFパルスが、ほぼ90°のRFパルス
からなることを特徴とする特許請求の範囲第(1)項に
記載の核磁気共鳴信号を用い、対象物から高速度で影像
情報を得る方法。 - (12)第2のRFパルスが、ほぼ180°のRFパル
スからなることを特徴とする特許請求の範囲第(2)項
に記載の核磁気共鳴信号を用い、対象物から高速度で影
像情報を得る方法。 - (13)予め決められた大きさの第2および第3の磁界
勾配が、符号化の間に対象物に与えられ、また、位相符
号化勾配および読出し勾配を交互に与えることにより、
前記の予め決められた大きさの第2および第3の磁界勾
配に基づくk空間の一部のみを貫通する軌線が生ずるこ
とを特徴とする特許請求の範囲第(1)項に記載の核磁
気共鳴信号を用い、対象物から高速度で影像情報を得る
方法。 - (14)与えられる符号化シーケンスおよび位相符号化
勾配は、k空間の50%以上が前記位相符号化勾配の方
向でサンプルされるようになっており、かつ、サンプル
されたk空間の対称部分は、像の位相を適切に表示させ
るのに十分であることを特徴とする特許請求の範囲第(
13)項に記載の核磁気共鳴信号を用い、対象物から高
速度で影像情報を得る方法。 - (15)与えられる符号化シーケンスおよび読出し勾配
は、k空間の50%以上が前記読出し勾配の方向でサン
プルされるようになっており、かつ、サンプルされたに
空間の対称部分は、像の位相を適切に表示させるのに十
分であることを特徴とする特許請求の範囲第(13)項
に記載の核磁気共鳴信号を用い、対象物から高速度で影
像情報を得る方法。 - (16)与えられる符号化シーケンスおよび勾配は、k
空間の一部のみがサンプルされる際、同数の点が、k空
間の完全なサンプリングにおける如く集められるように
なっているため、空間周波数応答は、信号のバンド幅が
維持されている間、増加するようになっていることを特
徴とする特許請求の範囲第(13)項に記載の核磁気共
鳴信号を用い、対象物から高速度で影像情報を得る方法
。 - (17)少なくとも2つの取得を、k空間の異なる部分
に関して行い、取得された部分は、変換前にk空間モザ
イクにおいて、k空間の連続有効範囲をつくるため、合
わせられることを特徴とする特許請求の範囲第(16)
項に記載の核磁気共鳴信号を用い、対象物から高速度で
影像情報を得る方法。 - (18)k空間の各象限において一つずつ、少なくとも
4つの重畳若しくは隣接している取得がなされ、その取
得された部分は、変換前にk空間モザイクに合わせられ
ることを特徴とする特許請求の範囲第(16)項に記載
の核磁気共鳴信号を用い、対象物から高速度で影像情報
を得る方法。 - (19)第1セットの2つの重畳若しくは隣接する取得
が、読出し勾配の方向でk空間の両側に関してなされ、
次に、第2セットの重畳若しくは隣接する取得が、読出
し勾配の方向でなされ、前記第2セットの取得は、位相
符号化勾配の方向で、前記第1セットの取得とインター
リーブされるようになっていることを特徴とする特許請
求の範囲第(16)項に記載の核磁気共鳴信号を用い、
対象物から高速度で影像情報を得る方法。
Applications Claiming Priority (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US937,529 | 1986-12-03 | ||
| US06/937,529 US4740748A (en) | 1986-12-03 | 1986-12-03 | Method of high-speed magnetic resonance imaging |
| US85,568 | 1987-08-17 | ||
| US07/085,568 US4767991A (en) | 1986-12-03 | 1987-08-17 | Method of high speed imaging with improved spatial resolution using partial k-space acquisitions |
| US109,091 | 1987-10-16 | ||
| US07/109,091 US4818942A (en) | 1986-12-03 | 1987-10-16 | Method of high-speed magnetic resonance imaging employing continuous wave readout gradient |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS63214246A true JPS63214246A (ja) | 1988-09-06 |
| JP2724830B2 JP2724830B2 (ja) | 1998-03-09 |
Family
ID=27375095
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP62306827A Expired - Lifetime JP2724830B2 (ja) | 1986-12-03 | 1987-12-03 | 核磁気共鳴信号を用い、対象物から高速度で影像情報を得る装置 |
Country Status (4)
| Country | Link |
|---|---|
| EP (1) | EP0270320B1 (ja) |
| JP (1) | JP2724830B2 (ja) |
| CA (1) | CA1256492A (ja) |
| DE (1) | DE3781869T2 (ja) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4833407A (en) * | 1987-06-24 | 1989-05-23 | Picker International, Inc. | Scan time reduction using conjugate symmetry and recalled echo |
| US4851779A (en) * | 1987-08-14 | 1989-07-25 | Picker International, Inc. | Three dimensional volume imaging with conjugate symmetrization |
| GB8909270D0 (en) * | 1989-04-24 | 1989-06-07 | Hafslund Nycomed Innovation | Method |
| DE4107515C2 (de) * | 1991-03-08 | 1995-03-09 | Siemens Ag | Verfahren zum Betrieb eines Kernspintomographiegeräts zur Erzeugung von Schnittbildern in beliebiger Lage und Schaltungsanordnung zur Druchführung des Verfahrens |
| JP3403751B2 (ja) * | 1992-03-31 | 2003-05-06 | 株式会社東芝 | 磁気共鳴映像装置 |
| US5229717A (en) * | 1992-05-22 | 1993-07-20 | General Electric Company | Simultaneous two-contrast fast spin echo NMR imaging |
| DE10108341C1 (de) * | 2001-02-21 | 2002-09-12 | Siemens Ag | Magnetresonanz-Spektroskopieverfahren mit einem Variieren von Phasen von HF-Pulsen |
| CN103961097B (zh) * | 2013-02-04 | 2016-08-31 | 上海联影医疗科技有限公司 | 一种磁共振扫描短te成像方法及磁共振扫描系统 |
| CN107076818B (zh) * | 2014-09-12 | 2019-11-05 | 皇家飞利浦有限公司 | 零回波时间mr成像 |
| WO2016081677A1 (en) * | 2014-11-19 | 2016-05-26 | The Johns Hopkins University | A steady pulsed imaging and labeling scheme for non-invasive perfusion |
| CN111329479B (zh) * | 2018-12-19 | 2023-09-29 | 西门子保健有限责任公司 | 用于提供mr图像的方法、磁共振成像设备和介质 |
| CN116299108A (zh) * | 2023-03-27 | 2023-06-23 | 广东省人民医院 | 基于3d-buda采集和联合低秩约束重建的mr成像方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5968656A (ja) * | 1982-09-10 | 1984-04-18 | エヌ・ベ−・フイリツプス・フル−イランペンフアブリケン | 核スピン分布決定装置 |
| JPS59105548A (ja) * | 1982-12-09 | 1984-06-18 | Yokogawa Hokushin Electric Corp | 核磁気共鳴による検査方法及び検査装置 |
| JPS6148753A (ja) * | 1984-08-17 | 1986-03-10 | Mitsubishi Electric Corp | 核磁気共鳴映像装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2079463B (en) * | 1980-03-14 | 1984-10-10 | Nat Res Dev | Nmr imaging method |
| 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 |
| JPH0763455B2 (ja) * | 1986-03-31 | 1995-07-12 | 株式会社東芝 | 磁気共鳴映像装置 |
-
1987
- 1987-11-27 DE DE19873781869 patent/DE3781869T2/de not_active Expired - Lifetime
- 1987-11-27 EP EP87310510A patent/EP0270320B1/en not_active Expired - Lifetime
- 1987-12-02 CA CA000553393A patent/CA1256492A/en not_active Expired
- 1987-12-03 JP JP62306827A patent/JP2724830B2/ja not_active Expired - Lifetime
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5968656A (ja) * | 1982-09-10 | 1984-04-18 | エヌ・ベ−・フイリツプス・フル−イランペンフアブリケン | 核スピン分布決定装置 |
| JPS59105548A (ja) * | 1982-12-09 | 1984-06-18 | Yokogawa Hokushin Electric Corp | 核磁気共鳴による検査方法及び検査装置 |
| JPS6148753A (ja) * | 1984-08-17 | 1986-03-10 | Mitsubishi Electric Corp | 核磁気共鳴映像装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2724830B2 (ja) | 1998-03-09 |
| CA1256492A (en) | 1989-06-27 |
| DE3781869D1 (de) | 1992-10-29 |
| EP0270320B1 (en) | 1992-09-23 |
| EP0270320A3 (en) | 1989-05-24 |
| DE3781869T2 (de) | 1993-01-07 |
| EP0270320A2 (en) | 1988-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4599565A (en) | Method and apparatus for rapid NMR imaging using multi-dimensional reconstruction techniques | |
| US4740748A (en) | Method of high-speed magnetic resonance imaging | |
| US4748410A (en) | Rapid NMR imaging system | |
| Posse et al. | High speed 1H spectroscopic imaging in human brain by echo planar spatial‐spectral encoding | |
| US5786692A (en) | Line scan diffusion imaging | |
| JP3529446B2 (ja) | Epi及びgrase mriにおける読み出し傾斜磁界極性の補正方法 | |
| US4767991A (en) | Method of high speed imaging with improved spatial resolution using partial k-space acquisitions | |
| US4521733A (en) | NMR Imaging of the transverse relaxation time using multiple spin echo sequences | |
| EP0212526B1 (en) | Method for reducing image artifacts due to periodic signal variations in nmr imaging | |
| JP4427152B2 (ja) | Mriシステムによる画像形成方法及びmriシステム | |
| US4940941A (en) | Method of high-speed magnetic resonance imaging | |
| US4731583A (en) | Method for reduction of MR image artifacts due to flowing nuclei by gradient moment nulling | |
| US5055789A (en) | Magnetic resonance imaging system | |
| US4818942A (en) | Method of high-speed magnetic resonance imaging employing continuous wave readout gradient | |
| EP1053484B1 (en) | Ghost artifact reduction | |
| US4612504A (en) | Method for removing the effects of baseline error components in NMR imaging applications | |
| US4728893A (en) | Increased signal-to-noise ratio in magnetic resonance images using synthesized conjugate symmetric data | |
| JPS63214246A (ja) | 核磁気共鳴信号を用い、対象物から高速度で影像情報を得る装置 | |
| JPH02149250A (ja) | 磁気共鳴映像装置 | |
| Block et al. | RARE spiral T2‐weighted imaging | |
| JPH0263435A (ja) | 磁気共鳴イメージング装置 | |
| JP3432506B2 (ja) | 磁気共鳴映像装置 | |
| JP3499927B2 (ja) | 磁気共鳴イメージング装置 | |
| JP2000175882A (ja) | Mrイメージング装置 | |
| JP2591405B2 (ja) | Mrイメージング装置のバイポーラグラジェントの調整法 |
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 |
|
| EXPY | Cancellation because of completion of term |