JPH06225863A - 磁気共鳴スペクトル画像撮影方法 - Google Patents
磁気共鳴スペクトル画像撮影方法Info
- Publication number
- JPH06225863A JPH06225863A JP5016146A JP1614693A JPH06225863A JP H06225863 A JPH06225863 A JP H06225863A JP 5016146 A JP5016146 A JP 5016146A JP 1614693 A JP1614693 A JP 1614693A JP H06225863 A JPH06225863 A JP H06225863A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic resonance
- measurement
- spectral
- dimensional
- image
- 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.)
- Pending
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】
【構成】計測１および２は独立したＭＲＳＩパルスシー
ケンスを表す。時間Ｔｗは計測１の終了と計測２の開始
を隔てる時間であり、スペクトル中の最長のＴ１よりも
十分長くする。ＭＲＳＩの計測を、繰返し時間ＴＲを
長，短に変えて２回以上実行し（１１および１２）、ス
ペクトルを規格化し（１３）、計測間のピーク強度の比
を画素値にとってコントラストを作成する（１４および
１５）。また、ＭＲＳＩとして、ケミカルシフトと空間
座標の一つを同時にエンコードする計測法を用いて短時
間で測定する。 【効果】生体内物質の縦緩和時間の相対値を化学種毎に
分布像として高速で得る。
ケンスを表す。時間Ｔｗは計測１の終了と計測２の開始
を隔てる時間であり、スペクトル中の最長のＴ１よりも
十分長くする。ＭＲＳＩの計測を、繰返し時間ＴＲを
長，短に変えて２回以上実行し（１１および１２）、ス
ペクトルを規格化し（１３）、計測間のピーク強度の比
を画素値にとってコントラストを作成する（１４および
１５）。また、ＭＲＳＩとして、ケミカルシフトと空間
座標の一つを同時にエンコードする計測法を用いて短時
間で測定する。 【効果】生体内物質の縦緩和時間の相対値を化学種毎に
分布像として高速で得る。
Description
【０００１】
【産業上の利用分野】本発明は、医療用磁気共鳴診断装
置を用いて、生体内の特定の化学種の分布を計測する方
法に関する。
置を用いて、生体内の特定の化学種の分布を計測する方
法に関する。
【０００２】
【従来の技術】従来、磁気共鳴イメージング装置（ＭＲ
Ｉ）を用いて、生体組織の縦緩和時間（Ｔ１）を測定す
る手法として、反転回復法や飽和回復法が用いられてい
る。これらの手法では、磁化反転用πパルスまたは飽和
用π／２パルスと計測用π／２パルスとの間隔を変えて
測定した複数の画像から、演算によりＴ１を画素毎に得
ることができる。しかし、これらの手法で得られるＴ１
はいわゆる密度画像に基づくものであり、画素に含まれ
る個々の化学種毎のＴ１ではない。これに対し、近年、
生体をin vivo 計測できるＭＲＳ（磁気共鳴スペクトロ
スコピー）の臨床応用が進展した結果、31Ｐを用いたＭ
ＲＳにおいて対象組織の特定の化学種のＴ１が、正常組
織と病変を生じた部位とでは異なることが見出された。
Ｉ）を用いて、生体組織の縦緩和時間（Ｔ１）を測定す
る手法として、反転回復法や飽和回復法が用いられてい
る。これらの手法では、磁化反転用πパルスまたは飽和
用π／２パルスと計測用π／２パルスとの間隔を変えて
測定した複数の画像から、演算によりＴ１を画素毎に得
ることができる。しかし、これらの手法で得られるＴ１
はいわゆる密度画像に基づくものであり、画素に含まれ
る個々の化学種毎のＴ１ではない。これに対し、近年、
生体をin vivo 計測できるＭＲＳ（磁気共鳴スペクトロ
スコピー）の臨床応用が進展した結果、31Ｐを用いたＭ
ＲＳにおいて対象組織の特定の化学種のＴ１が、正常組
織と病変を生じた部位とでは異なることが見出された。
【０００３】例えば、間島らは、表面コイルを用いて領
域限定した31ＰーＭＲＳを脳腫瘍，アルツハイマー病等
の患者に実施し、繰返し時間（ＴＲ）を１秒および１５
秒の２段階に変えて得たスペクトル中のＰＭＥ（リン酸
モノエステル）とβーＡＴＰ（アデノシン３リン酸）の
比が、ＴＲに応じてどのように変化するかを測定した。
その結果、腫瘍，アルツハイマー病等の患者ではＰＭＥ
のＴ１が長くなっていることを見出した（日本磁気共鳴
医学会誌第１１巻，補足１，１９４頁および同誌第１０
巻，補足２，１８７頁等）。また、成瀬らは、ＩＳＩＳ
法を用いて領域限定した31ＰーＭＲＳを脳神経疾患の患
者に実施し、得られたスペクトルの各ピークの値が、Ｔ
Ｒを２秒および１５秒の２段階に変えた場合にどのよう
に変化するかを調べた。その結果、やはり脳腫瘍，変性
疾患ではＰＭＥのＴ１が延長していることを見出した
(日本磁気共鳴医学会誌第１０巻，補足２，１８４頁
等)。
域限定した31ＰーＭＲＳを脳腫瘍，アルツハイマー病等
の患者に実施し、繰返し時間（ＴＲ）を１秒および１５
秒の２段階に変えて得たスペクトル中のＰＭＥ（リン酸
モノエステル）とβーＡＴＰ（アデノシン３リン酸）の
比が、ＴＲに応じてどのように変化するかを測定した。
その結果、腫瘍，アルツハイマー病等の患者ではＰＭＥ
のＴ１が長くなっていることを見出した（日本磁気共鳴
医学会誌第１１巻，補足１，１９４頁および同誌第１０
巻，補足２，１８７頁等）。また、成瀬らは、ＩＳＩＳ
法を用いて領域限定した31ＰーＭＲＳを脳神経疾患の患
者に実施し、得られたスペクトルの各ピークの値が、Ｔ
Ｒを２秒および１５秒の２段階に変えた場合にどのよう
に変化するかを調べた。その結果、やはり脳腫瘍，変性
疾患ではＰＭＥのＴ１が延長していることを見出した
(日本磁気共鳴医学会誌第１０巻，補足２，１８４頁
等)。
【０００４】このように、化学種毎、言い変えれば、ス
ペクトルの成分毎にＴ１を計測することは、疾病の診断
に際して有意義な知見を提供するものといえる。
ペクトルの成分毎にＴ１を計測することは、疾病の診断
に際して有意義な知見を提供するものといえる。
【０００５】
【発明が解決しようとする課題】上記従来技術では、表
面コイル法や、ＩＳＩＳ（イメージ セレクテッド イ
ンビボ スペクトロスコピー（Image selected invivo
Spectroscopy))法を用いることにより生体の限定された
単一領域からのスペクトルを測定しているため、測定し
た微小な領域の外の組織ではＴ１がどのように変化して
いるかを知ることができないという問題があった。
面コイル法や、ＩＳＩＳ（イメージ セレクテッド イ
ンビボ スペクトロスコピー（Image selected invivo
Spectroscopy))法を用いることにより生体の限定された
単一領域からのスペクトルを測定しているため、測定し
た微小な領域の外の組織ではＴ１がどのように変化して
いるかを知ることができないという問題があった。
【０００６】本発明の目的は、広い領域にわたるＴ１の
値を分布情報として得ることにより、上記問題を解決す
るものである。また、本発明の目的は、Ｔ１の分布情報
を得る測定では、一般に測定時間が長くなるという問題
を解決する。
値を分布情報として得ることにより、上記問題を解決す
るものである。また、本発明の目的は、Ｔ１の分布情報
を得る測定では、一般に測定時間が長くなるという問題
を解決する。
【０００７】
【課題を解決するための手段】上記課題を解決するため
に本発明では、一回の測定で複数の領域のスペクトルが
得られるスペクトロスコピック・イメージング法（以
下、ＭＲＳＩ法と略記する）を採用し、これと複数の繰
返し時間（ＴＲ）を組み合わせる。例えば代表的なＭＲ
ＳＩ法である３Ｄ−ＣＳＩ法（３次元ケミカルシフト・
イメージング）法あるいは４Ｄ−ＣＳＩ（４次元ケミカ
ルシフト・イメージング）法を、ＴＲを１秒および１５
秒として２回実施して、ＰＭＥ等の着目するピークの強
度の比を画素値に採って画像を作成する。これによれ
ば、ＰＭＥのＴ１の相対的な長さが組織毎にどのように
分布しているかが容易に把握できる。
に本発明では、一回の測定で複数の領域のスペクトルが
得られるスペクトロスコピック・イメージング法（以
下、ＭＲＳＩ法と略記する）を採用し、これと複数の繰
返し時間（ＴＲ）を組み合わせる。例えば代表的なＭＲ
ＳＩ法である３Ｄ−ＣＳＩ法（３次元ケミカルシフト・
イメージング）法あるいは４Ｄ−ＣＳＩ（４次元ケミカ
ルシフト・イメージング）法を、ＴＲを１秒および１５
秒として２回実施して、ＰＭＥ等の着目するピークの強
度の比を画素値に採って画像を作成する。これによれ
ば、ＰＭＥのＴ１の相対的な長さが組織毎にどのように
分布しているかが容易に把握できる。
【０００８】更に、高速な計測が要求される場合には、
空間の１次元もしくは２次元の座標を位相エンコード
し、残りの１次元の空間座標とスペクトル座標を傾斜磁
場の矩形波駆動や正弦波駆動により同時に周波数エンコ
ードする、いわゆる、高速スペクトロスコピック・イメ
ージング法と複数のＴＲを組み合わせる。
空間の１次元もしくは２次元の座標を位相エンコード
し、残りの１次元の空間座標とスペクトル座標を傾斜磁
場の矩形波駆動や正弦波駆動により同時に周波数エンコ
ードする、いわゆる、高速スペクトロスコピック・イメ
ージング法と複数のＴＲを組み合わせる。
【０００９】
【作用】本発明ではシングルボクセル法の代わりにマル
チボクセル法を採用しているため、組織の各部分におけ
る相対的なＴ１が分布画像として得られる。このため、
単一部位のＴ１の値を数値で評価する場合に比べて、病
変の有無や分布をはるかに容易に把握できる。また、シ
ングルボクセル法とは異なり、測定領域を選択する手続
きが不要であり、病変の初期段階でのスクリーニングで
も異常な部位を見逃す恐れが少ない。また、高速スペク
トロスコピック・イメージング法を用いた場合には、短
時間で相対的なＴ１の分布画像が得られるため、病態の
経時変化を追跡することが可能になる。
チボクセル法を採用しているため、組織の各部分におけ
る相対的なＴ１が分布画像として得られる。このため、
単一部位のＴ１の値を数値で評価する場合に比べて、病
変の有無や分布をはるかに容易に把握できる。また、シ
ングルボクセル法とは異なり、測定領域を選択する手続
きが不要であり、病変の初期段階でのスクリーニングで
も異常な部位を見逃す恐れが少ない。また、高速スペク
トロスコピック・イメージング法を用いた場合には、短
時間で相対的なＴ１の分布画像が得られるため、病態の
経時変化を追跡することが可能になる。
【００１０】
【実施例】図１は、本発明の基本的な手順を示す。同図
において、計測１および２は独立したＭＲＳＩパルスシ
ーケンスを表す。時間Ｔｗは計測１の終了と計測２の開
始を隔てる時間であり、スペクトル中の最長のＴ１より
も十分長いことが望ましい。具体的には 1Ｈでは５秒、
31Ｐでは１５秒程度とする。
において、計測１および２は独立したＭＲＳＩパルスシ
ーケンスを表す。時間Ｔｗは計測１の終了と計測２の開
始を隔てる時間であり、スペクトル中の最長のＴ１より
も十分長いことが望ましい。具体的には 1Ｈでは５秒、
31Ｐでは１５秒程度とする。
【００１１】ＭＲＳＩのシーケンスとして、例えば、図
３に示す３Ｄ−ＣＳＩ法を用いることができる。図１の
処理１１において、ＴＲ＝１秒として３Ｄ−ＣＳＩの測
定を行い第一のスペクトル群（ボクセル毎のスペクト
ル）を得る。次に、処理１２においてＴＲ＝１５秒とし
て３Ｄ−ＣＳＩを行い、第二のスペクトル群を得る。次
いで、処理１３においてスペクトル線の内の一つ、例え
ば、31Ｐの場合であれば、ＰＭＥの強度を基準とし、他
のスペクトルを規格化する。規格化されたスペクトルを
μ1(ｉ）とする。ここにｉ＝１，２，…はスペクトル線
の番号を、上付き添字は計測時のＴＲを表す。処理１４
において第一のスペクトルで第二のスペクトルを除し、
比のスペクトルＫ(ｉ）＝μ15(ｉ)／μ1(ｉ)を得る。最
後に処理１５においてＫ(ｉ)を画素の濃度に対応させた
分布画像を作成する。
３に示す３Ｄ−ＣＳＩ法を用いることができる。図１の
処理１１において、ＴＲ＝１秒として３Ｄ−ＣＳＩの測
定を行い第一のスペクトル群（ボクセル毎のスペクト
ル）を得る。次に、処理１２においてＴＲ＝１５秒とし
て３Ｄ−ＣＳＩを行い、第二のスペクトル群を得る。次
いで、処理１３においてスペクトル線の内の一つ、例え
ば、31Ｐの場合であれば、ＰＭＥの強度を基準とし、他
のスペクトルを規格化する。規格化されたスペクトルを
μ1(ｉ）とする。ここにｉ＝１，２，…はスペクトル線
の番号を、上付き添字は計測時のＴＲを表す。処理１４
において第一のスペクトルで第二のスペクトルを除し、
比のスペクトルＫ(ｉ）＝μ15(ｉ)／μ1(ｉ)を得る。最
後に処理１５においてＫ(ｉ)を画素の濃度に対応させた
分布画像を作成する。
【００１２】μ15(ｉ)はＴ１に比べて十分長い繰返し時
間（ＴＲ）を設定して得たものであるため、各スペクト
ル線は濃度に比例した強度を示す。これに対して、μ
1(ｉ)はＴＲが短く、縦磁化が熱平衡の値に復帰する余
裕がないため部分的な飽和が生じ、Ｔ１の長いスペクト
ル線ほど強度が減衰する。一例として、31Ｐスペクトル
の変化を模式的に図２に示す。ヒト脳の２.０ テスラに
おけるＴ１はボスカら（マグネティクレゾナンス・イン
・メディスン誌，第１３巻，２２８−２３８頁、１９９
０年；M.D.Boska et.al.,Magn.Reson.Med.,13,228(199
0))によると、３種類のＡＴＰのＴ１が０.６〜１秒 で
あるのに対し、ＰＭＥとクレアチンリン酸（ＰＣｒ）で
はそれぞれ１.７ｓ，２.７ｓである。従って後者はＴＲ
を１秒とした場合の減衰量が大きい（図２のピーク２１
と２２）。
間（ＴＲ）を設定して得たものであるため、各スペクト
ル線は濃度に比例した強度を示す。これに対して、μ
1(ｉ)はＴＲが短く、縦磁化が熱平衡の値に復帰する余
裕がないため部分的な飽和が生じ、Ｔ１の長いスペクト
ル線ほど強度が減衰する。一例として、31Ｐスペクトル
の変化を模式的に図２に示す。ヒト脳の２.０ テスラに
おけるＴ１はボスカら（マグネティクレゾナンス・イン
・メディスン誌，第１３巻，２２８−２３８頁、１９９
０年；M.D.Boska et.al.,Magn.Reson.Med.,13,228(199
0))によると、３種類のＡＴＰのＴ１が０.６〜１秒 で
あるのに対し、ＰＭＥとクレアチンリン酸（ＰＣｒ）で
はそれぞれ１.７ｓ，２.７ｓである。従って後者はＴＲ
を１秒とした場合の減衰量が大きい（図２のピーク２１
と２２）。
【００１３】なお、ＴＲを１秒としたスペクトルではピ
ーク間の比が変化するばかりでなく、全体的に強度が低
下する。この効果は図２の表現では省いてあるが、図１
ではスペクトルを計測毎に規格化することにより排除し
ている。
ーク間の比が変化するばかりでなく、全体的に強度が低
下する。この効果は図２の表現では省いてあるが、図１
ではスペクトルを計測毎に規格化することにより排除し
ている。
【００１４】また、スペクトルのピーク強度の減衰率は
ＴＲとＴ１の比に依存し、１−exp（−ＴＲ／Ｔ１）で
与えられる。
ＴＲとＴ１の比に依存し、１−exp（−ＴＲ／Ｔ１）で
与えられる。
【００１５】なお、31ＰではＴ１の比較的長い化学種が
あるため、完全な縦緩和を達成するには約１５秒という
長いＴＲが必要になるが、 1Ｈでは生体内の主要な化学
種のＴ１は１〜２秒であり、より短いＴＲを用いること
ができる。
あるため、完全な縦緩和を達成するには約１５秒という
長いＴＲが必要になるが、 1Ｈでは生体内の主要な化学
種のＴ１は１〜２秒であり、より短いＴＲを用いること
ができる。
【００１６】３Ｄ−ＣＳＩ法を用いた場合の短所は、繰
返し時間（ＴＲ）を挟んで空間座標に対する位相エンコ
ードのステップ数だけ計測を繰り返すため、全体の計測
時間が長くなりがちな点である。即ち、空間のｘ，ｙ２
方向の座標を位相エンコードするため、ｘ，ｙ方向のマ
トリクス数Ｎｘ，Ｎｙを各々１６と比較的少数に抑えた
場合でも、繰返しは１６×１６＝２５６回となる。ＴＲ
＝１５秒とすると、１５×１６×１６＝３８４０秒（６
４分）かかることになり、臨床応用には不向きとなる。
図４に示すような４Ｄ−ＣＳＩ法では、この問題は更に
深刻になる。
返し時間（ＴＲ）を挟んで空間座標に対する位相エンコ
ードのステップ数だけ計測を繰り返すため、全体の計測
時間が長くなりがちな点である。即ち、空間のｘ，ｙ２
方向の座標を位相エンコードするため、ｘ，ｙ方向のマ
トリクス数Ｎｘ，Ｎｙを各々１６と比較的少数に抑えた
場合でも、繰返しは１６×１６＝２５６回となる。ＴＲ
＝１５秒とすると、１５×１６×１６＝３８４０秒（６
４分）かかることになり、臨床応用には不向きとなる。
図４に示すような４Ｄ−ＣＳＩ法では、この問題は更に
深刻になる。
【００１７】測定時間を短縮する方法として、第二のス
ペクトル群の測定に於いて完全な縦緩和を行わず、より
短いＴＲを用いる方法が考えられる。例えば、計測２で
ＴＲを１５秒の代わりに７秒とすれば測定時間を短縮で
きるが、Ｔ１の精度は低くなる。計測１のＴＲをＴＲ
１、計測２のＴＲをＴＲ２とすれば、２回の計測間での
スペクトルのピーク強度低下率Ｒ（＜１）は次式で与え
られる。但し、Ｔ２による減衰は無視する。
ペクトル群の測定に於いて完全な縦緩和を行わず、より
短いＴＲを用いる方法が考えられる。例えば、計測２で
ＴＲを１５秒の代わりに７秒とすれば測定時間を短縮で
きるが、Ｔ１の精度は低くなる。計測１のＴＲをＴＲ
１、計測２のＴＲをＴＲ２とすれば、２回の計測間での
スペクトルのピーク強度低下率Ｒ（＜１）は次式で与え
られる。但し、Ｔ２による減衰は無視する。
【００１８】
【数１】
【００１９】ここにＴ１(ｉ)はｉ番目のピークのＴ１を
表す。理想的な測定条件はＴＲ２＞Ｔ１(ｉ)となるよう
にすることであるが、この条件が満たされない場合には
Ｒは相対的に大きくなり、Ｔ１のコントラストが小さく
なる。言い替えると測定の精度が低下する。
表す。理想的な測定条件はＴＲ２＞Ｔ１(ｉ)となるよう
にすることであるが、この条件が満たされない場合には
Ｒは相対的に大きくなり、Ｔ１のコントラストが小さく
なる。言い替えると測定の精度が低下する。
【００２０】第２の実施例では、ＴＲをＴ１よりも小さ
くするとともにＲＦパルスのフリップ角を９０度よりも
小さくすることにより、信号強度の低下を最小限に留め
る。フリップ角をαとした時、信号強度Ｉは次式、
くするとともにＲＦパルスのフリップ角を９０度よりも
小さくすることにより、信号強度の低下を最小限に留め
る。フリップ角をαとした時、信号強度Ｉは次式、
【００２１】
【数２】
【００２２】に従い、また、Ｉを最大にするα（エルン
スト角）は、
スト角）は、
【００２３】
【数３】
【００２４】で与えられることが知られている。ただ
し、Ｅ＝exp(−ＴＲ／Ｔ１）である。
し、Ｅ＝exp(−ＴＲ／Ｔ１）である。
【００２５】ここでは、スペクトル中の最も注目するピ
ーク（第ｉ番目のピークとする）を最大にするフリップ
角を２回の計測毎に別個に設定する。即ち、計測１での
フリップ角α１を
ーク（第ｉ番目のピークとする）を最大にするフリップ
角を２回の計測毎に別個に設定する。即ち、計測１での
フリップ角α１を
【００２６】
【数４】
【００２７】により、また計測２でのフリップ角を、
【００２８】
【数５】
【００２９】により設定する。このようにすれば、部分
的な飽和による信号低下を最小限に留めつつ、ＴＲを短
縮することができる。但し、コントラスト（２種類のＩ
の比）はＴ１の複雑な関数になる。
的な飽和による信号低下を最小限に留めつつ、ＴＲを短
縮することができる。但し、コントラスト（２種類のＩ
の比）はＴ１の複雑な関数になる。
【００３０】さて、計測時間を大きく短縮するには、位
相エンコードのための反復回数を低減することが最も有
効である。
相エンコードのための反復回数を低減することが最も有
効である。
【００３１】第３の実施例では、高速ＭＲＳＩ法として
提案されているマンスフィールドらの方法、松井らの方
法やウェブらの方法（ジャーナル．オブ．マグネティッ
クレゾナンス誌第６７巻，４７６−４９０頁、１９８６
年：S.Matsui,K.Sekihara，and H.Khono,Journal of Ma
gnetic Resonance，67,476,(1986)）を利用する。これ
らの方法に共通する特徴は、傾斜磁場の内の一つ（以
下、便宜的にこれをｘ方向の傾斜磁場、Ｇｘとする。）
を連続的に正負に切り替え、勾配エコー列を計測するこ
とによりｘ座標に関する情報とケミカルシフトに関する
情報を同時に得、スペクトロスコピック・イメージング
を高速に実行している点である。傾斜磁場の駆動方法
は、正弦波駆動を用いる方法も知られている。以下、こ
れらを一括して多重エンコード法という。
提案されているマンスフィールドらの方法、松井らの方
法やウェブらの方法（ジャーナル．オブ．マグネティッ
クレゾナンス誌第６７巻，４７６−４９０頁、１９８６
年：S.Matsui,K.Sekihara，and H.Khono,Journal of Ma
gnetic Resonance，67,476,(1986)）を利用する。これ
らの方法に共通する特徴は、傾斜磁場の内の一つ（以
下、便宜的にこれをｘ方向の傾斜磁場、Ｇｘとする。）
を連続的に正負に切り替え、勾配エコー列を計測するこ
とによりｘ座標に関する情報とケミカルシフトに関する
情報を同時に得、スペクトロスコピック・イメージング
を高速に実行している点である。傾斜磁場の駆動方法
は、正弦波駆動を用いる方法も知られている。以下、こ
れらを一括して多重エンコード法という。
【００３２】図５に一例として、松井らによるシーケン
スを示す。ここでは、計測１，２において位相エンコー
ドを変えてそれぞれＮｙ回測定を繰り返す。マトリクス
数を１６×１６，ＴＲ１を１秒、ＴＲ２を１５秒に取っ
た場合でも計測の所要時間は５分程度で済む。但し、多
重エンコード法では計測スペクトル帯域幅は空間座標の
標本点数と関連しており、強い制限を受ける。従って、
31Ｐ，13Ｃなどのスペクトル帯域の広い核種へ多重エン
コード法を適用する場合には、選択励起，選択飽和等の
スペクトル軸上の折り返しアーチファクトを防止する手
段や、ハイブリッド・スキャン等の計測スペクトル帯域
を拡大する手段を併用することが望ましい。
スを示す。ここでは、計測１，２において位相エンコー
ドを変えてそれぞれＮｙ回測定を繰り返す。マトリクス
数を１６×１６，ＴＲ１を１秒、ＴＲ２を１５秒に取っ
た場合でも計測の所要時間は５分程度で済む。但し、多
重エンコード法では計測スペクトル帯域幅は空間座標の
標本点数と関連しており、強い制限を受ける。従って、
31Ｐ，13Ｃなどのスペクトル帯域の広い核種へ多重エン
コード法を適用する場合には、選択励起，選択飽和等の
スペクトル軸上の折り返しアーチファクトを防止する手
段や、ハイブリッド・スキャン等の計測スペクトル帯域
を拡大する手段を併用することが望ましい。
【００３３】ここまではＴＲを２通りに変える場合につ
いて述べてきたが、３通り以上に変えてもよい。その場
合には複数の数１あるいは、数２に適合するＴ１を求め
ることにより、Ｔ１の相対値だけではなく絶対値を精度
良く求めることが可能になる。
いて述べてきたが、３通り以上に変えてもよい。その場
合には複数の数１あるいは、数２に適合するＴ１を求め
ることにより、Ｔ１の相対値だけではなく絶対値を精度
良く求めることが可能になる。
【００３４】
【発明の効果】本発明では化学種毎の相対的なＴ１の分
布像を得ることができるため、疾患の生じている部位を
容易に発見できる。また、高速で分布像を得ることがで
きる。
布像を得ることができるため、疾患の生じている部位を
容易に発見できる。また、高速で分布像を得ることがで
きる。
【図１】本発明による相対Ｔ１分布測定の手順を示す説
明図。
明図。
【図２】31ＰスペクトルのＴＲによる変化を模式的に示
す波形図。
す波形図。
【図３】３Ｄ−ＣＳＩ法のパルスシーケンスを示す図。
【図４】４Ｄ−ＣＳＩ法のパルスシーケンスを示す図。
【図５】高速ＭＲＳＩ法のシーケンスの一例を示す図。
１１…第一回目のＭＲＳＩシーケンス、１２…第二回目
のＭＲＳＩシーケンス。
のＭＲＳＩシーケンス。
Claims (4)
- 【請求項１】磁気共鳴画像撮影装置において、パルス系
列の繰返し時間を変えて複数のスペクトロスコピック画
像計測を行い、各画像計測間で各スペクトル成分の比を
計算しこの値を画素値とすることにより、各々の前記ス
ペクトル成分の縦緩和時間を反映した分布画像を得るこ
とを特徴とする磁気共鳴スペクトル画像撮影方法。 - 【請求項２】請求項１において、前記磁気共鳴スペクト
ロスコピック画像の撮影方法が、３次元あるいは４次元
ケミカルシフト・イメージング法である磁気共鳴スペク
トル画像撮影方法。 - 【請求項３】請求項１において、前記磁気共鳴スペクト
ロスコピック画像の撮影方法が、１次元もしくは２次元
の空間座標の位相エンコードと、傾斜磁場の矩形波もし
くは正弦波駆動による残りの１次元の空間座標およびス
ペクトル座標の同時エンコードとを用いる磁気共鳴スペ
クトル画像撮影方法。 - 【請求項４】請求項２において、前記磁気共鳴スペクト
ロスコピック画像の撮影方法が、９０度未満のフリップ
角を用いる３次元あるいは４次元ケミカルシフト・イメ
ージング法である磁気共鳴スペクトル画像撮影方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5016146A JPH06225863A (ja) | 1993-02-03 | 1993-02-03 | 磁気共鳴スペクトル画像撮影方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5016146A JPH06225863A (ja) | 1993-02-03 | 1993-02-03 | 磁気共鳴スペクトル画像撮影方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH06225863A true JPH06225863A (ja) | 1994-08-16 |
Family
ID=11908365
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5016146A Pending JPH06225863A (ja) | 1993-02-03 | 1993-02-03 | 磁気共鳴スペクトル画像撮影方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH06225863A (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009532163A (ja) * | 2006-04-06 | 2009-09-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁気共鳴装置及び方法 |
| JP2009534101A (ja) * | 2006-04-19 | 2009-09-24 | ウイスコンシン アラムナイ リサーチ フオンデーシヨン | 多重ピークを備えた種の磁気共鳴スペクトロスコピー |
| JP2019144051A (ja) * | 2018-02-19 | 2019-08-29 | 日本ポリプロ株式会社 | 熱可塑性樹脂の定量測定方法 |
-
1993
- 1993-02-03 JP JP5016146A patent/JPH06225863A/ja active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009532163A (ja) * | 2006-04-06 | 2009-09-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁気共鳴装置及び方法 |
| JP2009534101A (ja) * | 2006-04-19 | 2009-09-24 | ウイスコンシン アラムナイ リサーチ フオンデーシヨン | 多重ピークを備えた種の磁気共鳴スペクトロスコピー |
| JP2019144051A (ja) * | 2018-02-19 | 2019-08-29 | 日本ポリプロ株式会社 | 熱可塑性樹脂の定量測定方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2021201998B2 (en) | Selective sampling magnetic resonance-based method for assessing structural spatial frequencies | |
| JP3512482B2 (ja) | 磁気共鳴映像装置 | |
| US4585992A (en) | NMR imaging methods | |
| JP4936864B2 (ja) | 磁気共鳴イメージング装置 | |
| US6094591A (en) | Measurement of coronary flow reserve with MR oximetry | |
| US6704593B2 (en) | Realtime MR scan prescription using physiological information | |
| US6377834B1 (en) | Real time in vivo measurement of temperature changes with contrast enhanced NMR imaging | |
| US5604435A (en) | Spiral scanning method for monitoring physiological changes | |
| US8498688B2 (en) | Magnetic resonance device and method | |
| US9116218B2 (en) | System and method for tissue specific MR imaging of metabolites using spectral-spatially formed stimulated echo | |
| US10191132B2 (en) | Method and apparatus for multi-slice imaging of T2-relaxation time | |
| US8315450B2 (en) | Method and system for display of medical image data | |
| US8112145B2 (en) | MRI method for assessing myocardial viability | |
| JP4101300B2 (ja) | スペクトロスコピー磁気共鳴画像化方法及び装置 | |
| US7705593B1 (en) | Apparatus and method for detecting and classifying atherosclerotic plaque hemorrhage | |
| JP5472945B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3699963B2 (ja) | 磁気共鳴映像装置 | |
| JP2005007181A (ja) | Ｍｒ分光法において代謝産物の信号分離を改良させた方法及び装置 | |
| JPH06225863A (ja) | 磁気共鳴スペクトル画像撮影方法 | |
| JP3153573B2 (ja) | 磁気共鳴装置 | |
| Lai et al. | FAIR exempting separate T 1 measurement (FAIREST): a novel technique for online quantitative perfusion imaging and multi‐contrast fMRI | |
| JP4718698B2 (ja) | 磁気共鳴映像装置 | |
| JP3690874B2 (ja) | 磁気共鳴イメージング装置 | |
| WO2008132686A1 (en) | Quantification for mr parameters such as t1 or t2 in a sub-region of a subject | |
| JP2005161102A (ja) | 磁気共鳴診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20060427 |