JPH0265843A - 核磁化分布像を得る方法 - Google Patents
核磁化分布像を得る方法Info
- Publication number
- JPH0265843A JPH0265843A JP1174335A JP17433589A JPH0265843A JP H0265843 A JPH0265843 A JP H0265843A JP 1174335 A JP1174335 A JP 1174335A JP 17433589 A JP17433589 A JP 17433589A JP H0265843 A JPH0265843 A JP H0265843A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- distribution
- spin resonance
- spin
- time constant
- 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
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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo 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
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
【発明の詳細な説明】
（産業上の利用分野）
本発明はスピン共鳴信号を発生させるために静磁場の存
在の下で身体中で核スピンを励起し、その後静磁場に重
畳する傾斜磁場の存在の下でスピン共鳴信号をサンプル
化し、かくしてサンプル化したスピン共鳴信号を信号変
換により核磁化分布像に変換するようにした身体の一部
分の核磁化分布像を得る方法に関するものである。
在の下で身体中で核スピンを励起し、その後静磁場に重
畳する傾斜磁場の存在の下でスピン共鳴信号をサンプル
化し、かくしてサンプル化したスピン共鳴信号を信号変
換により核磁化分布像に変換するようにした身体の一部
分の核磁化分布像を得る方法に関するものである。
サンプル化されたスピン共鳴信号は上記方法に従つパル
スエコー信号とすることができる。
スエコー信号とすることができる。
また、本発明は均一静磁場を発生すると共に傾斜磁場を
発生する磁石系を具え、かつスピン共鳴信号を発生しか
つ検出するＲＦ送信器および受信器と、発生したスピン
共鳴信号をサンプリングするサンプリング装置と、サン
プル化したスピン共鳴信号を変換する信号変換装置と、
変換した信号を表示する表示装置と、前記磁石系、ＲＦ
送信器および受信器と、サンプリング装置と、信号変換
装置と、表示装置とを制御する制御装置とを具えるＭＲ
Ｉ装置に関するものである。
発生する磁石系を具え、かつスピン共鳴信号を発生しか
つ検出するＲＦ送信器および受信器と、発生したスピン
共鳴信号をサンプリングするサンプリング装置と、サン
プル化したスピン共鳴信号を変換する信号変換装置と、
変換した信号を表示する表示装置と、前記磁石系、ＲＦ
送信器および受信器と、サンプリング装置と、信号変換
装置と、表示装置とを制御する制御装置とを具えるＭＲ
Ｉ装置に関するものである。
（従来の技術及び発明が解決しようとする課題）この種
方法および装置は米国特許明細書第４，０７０．６１１
号から既知であり、検査すべき物体の核磁化分布像の密
度その他の特性を基にビデオ情報を得る場合にしばしば
用いる。更に、ＰＣＴ出願第１０８５１０２２６４号に
は物体、または、その一部分のかかる像をシンセサイズ
する方法が記載されている。この場合、像の一つの点の
信号強度は局部磁化ＭＯ１局部縦方向および横方向弛緩
時定数Ｔ１およびＴ２並びに測定方法に用いるパラメー
タ、例えば、測定サイクルＴＲの繰返し時間並びに測定
サイクルの始端時の励起およびスピン共鳴信号のサンプ
リング瞬時の間に経過する待ち時間ＴＥによって決まる
。
方法および装置は米国特許明細書第４，０７０．６１１
号から既知であり、検査すべき物体の核磁化分布像の密
度その他の特性を基にビデオ情報を得る場合にしばしば
用いる。更に、ＰＣＴ出願第１０８５１０２２６４号に
は物体、または、その一部分のかかる像をシンセサイズ
する方法が記載されている。この場合、像の一つの点の
信号強度は局部磁化ＭＯ１局部縦方向および横方向弛緩
時定数Ｔ１およびＴ２並びに測定方法に用いるパラメー
タ、例えば、測定サイクルＴＲの繰返し時間並びに測定
サイクルの始端時の励起およびスピン共鳴信号のサンプ
リング瞬時の間に経過する待ち時間ＴＥによって決まる
。
測定パラメータＴＲおよび測定パラメータＴＥが変化す
る箇所で多数のイメージを決めることにより、磁化ＭＯ
並びに縦方向および横方向弛緩時定数ＴＩおよびＴ２を
を各画素に対して計算することができる。このデータが
任意の測定パラメータＴＲおよびＴＥに対し既知である
ため、合成像は、これら像が測定パラメータＴＲおよび
ＴＥにより決まる測定中に得られたかのように再構成す
ることができる。これがため、像の種々の部分間のコン
トラストを最適化し得る像をシンセサイズすることがで
きる。
る箇所で多数のイメージを決めることにより、磁化ＭＯ
並びに縦方向および横方向弛緩時定数ＴＩおよびＴ２を
を各画素に対して計算することができる。このデータが
任意の測定パラメータＴＲおよびＴＥに対し既知である
ため、合成像は、これら像が測定パラメータＴＲおよび
ＴＥにより決まる測定中に得られたかのように再構成す
ることができる。これがため、像の種々の部分間のコン
トラストを最適化し得る像をシンセサイズすることがで
きる。
従って、かかる方法は、合成像が再構成される前に、多
数の完全な像を常時必要とする。これは必要な測定時間
に対し、即ち、物体がＭＲＩ装置内に留どまる必要のあ
る間の時間に対する直接の結果である。
数の完全な像を常時必要とする。これは必要な測定時間
に対し、即ち、物体がＭＲＩ装置内に留どまる必要のあ
る間の時間に対する直接の結果である。
本発明の目的は、物体の画素の種々のスピン共鳴パラメ
ータ、即ち、スピン密度（または平衡磁化ＭＯ）および
横方向弛緩時定数Ｔ２を迅速に決めることができ、スピ
ン密度分布に対してはただ１つのエコーマトリックスの
みを必要とし、即ち、ただ１つの像層波数マトリックス
を必要とするＭＲ１方法および装置を提供せんとするに
ある。
ータ、即ち、スピン密度（または平衡磁化ＭＯ）および
横方向弛緩時定数Ｔ２を迅速に決めることができ、スピ
ン密度分布に対してはただ１つのエコーマトリックスの
みを必要とし、即ち、ただ１つの像層波数マトリックス
を必要とするＭＲ１方法および装置を提供せんとするに
ある。
（課題を解決するための手段）
本発明核磁化分布像を得る方法はスピン共鳴信号を発生
させるために静磁場の存在の下で身体中で核スピンを励
起し、その後静磁場に重畳する傾斜磁場の存在の下でス
ピン共鳴信号をサンプル化し、かくしてサンプル化した
スピン共鳴信号を信号変換により核磁化分布像に変換す
るようにした身体の一部分の核磁化分布像を得るに当た
り、この核磁化分布を横方向弛緩時定数分布とし、スピ
ンエコー信号を発生する励起スピンの横方向弛緩により
非対称性を生ぜしめるサンプル化スピン共鳴信号の相対
的な非対称性から、部分横方向弛緩時定数の測定値を取
出すようにしたことを特徴とする。
させるために静磁場の存在の下で身体中で核スピンを励
起し、その後静磁場に重畳する傾斜磁場の存在の下でス
ピン共鳴信号をサンプル化し、かくしてサンプル化した
スピン共鳴信号を信号変換により核磁化分布像に変換す
るようにした身体の一部分の核磁化分布像を得るに当た
り、この核磁化分布を横方向弛緩時定数分布とし、スピ
ンエコー信号を発生する励起スピンの横方向弛緩により
非対称性を生ぜしめるサンプル化スピン共鳴信号の相対
的な非対称性から、部分横方向弛緩時定数の測定値を取
出すようにしたことを特徴とする。
本発明方法では横方向弛緩が時間の関数としてスピンエ
コー信号に影響を与えると云う事実を用いる。この影響
は、Ｔ２を横方向弛緩時定数とする場合における指数関
数（ｅｘｐ−ｔ／Ｔ２）である。この時定数７２（ｘ、
ｙ、ｚ）は位置に依存し、再構成すべき像の一方の画素
から他方の画素に亘って相違させることができる。
コー信号に影響を与えると云う事実を用いる。この影響
は、Ｔ２を横方向弛緩時定数とする場合における指数関
数（ｅｘｐ−ｔ／Ｔ２）である。この時定数７２（ｘ、
ｙ、ｚ）は位置に依存し、再構成すべき像の一方の画素
から他方の画素に亘って相違させることができる。
本発明スピン密度分布および横方向弛緩時定数分布をき
める方法の変形例では、スピン共鳴信号を発生させるた
めに静磁場の存在の下で身体中で核スピンを励起し、そ
の後静磁場に重畳すると共に第１の方向を有する傾斜磁
場の存在の下でスピン共鳴信号をサンプル化し、次いで
、がくしでサンプル化したスピン共鳴信号からの信号変
換により核磁化分布像を測定するようにした身体の一部
分のスピン密度分布および横方向弛緩時定数Ｔ２の分布
を測定するに当たり、スピン共鳴信号の信号サンプルの
フーリエ変換後、前記第１方向における変換信号の虚数
部分を積分し、その後実数部分と変換信号の積分虚数部
分との商を決めるようにしたことを特徴とする。
める方法の変形例では、スピン共鳴信号を発生させるた
めに静磁場の存在の下で身体中で核スピンを励起し、そ
の後静磁場に重畳すると共に第１の方向を有する傾斜磁
場の存在の下でスピン共鳴信号をサンプル化し、次いで
、がくしでサンプル化したスピン共鳴信号からの信号変
換により核磁化分布像を測定するようにした身体の一部
分のスピン密度分布および横方向弛緩時定数Ｔ２の分布
を測定するに当たり、スピン共鳴信号の信号サンプルの
フーリエ変換後、前記第１方向における変換信号の虚数
部分を積分し、その後実数部分と変換信号の積分虚数部
分との商を決めるようにしたことを特徴とする。
横方向弛緩のない場合に、即ち、初期の大きな横方向弛
緩時定数Ｔ２に対し、信号サンプルのフーリエ変換が純
粋に実数となる位相の信号を得ることができる。
緩時定数Ｔ２に対し、信号サンプルのフーリエ変換が純
粋に実数となる位相の信号を得ることができる。
他の好適な例ではスピン共鳴信号の位相を補正して、ツ
ー１フエ変換後に信号が横方向弛緩のない場合のフーリ
エ変換後に純粋に実数となるような条件を満足し、その
後これら位相補正信号サンプルのフーリエ変換に続き、
変換信号の虚数部分を積分し、次いで変換信号の実数部
分および積分した虚数部分の商を決め得るようにする。
ー１フエ変換後に信号が横方向弛緩のない場合のフーリ
エ変換後に純粋に実数となるような条件を満足し、その
後これら位相補正信号サンプルのフーリエ変換に続き、
変換信号の虚数部分を積分し、次いで変換信号の実数部
分および積分した虚数部分の商を決め得るようにする。
本発明方法の好適な例では、変換した信号の積分虚数部
分が前記間を決めるには小さすぎる画素に対しては、商
の決定を、周囲の画素に基づく推定、または結果を零に
等置する等の他の演算で行い得るようにする。
分が前記間を決めるには小さすぎる画素に対しては、商
の決定を、周囲の画素に基づく推定、または結果を零に
等置する等の他の演算で行い得るようにする。
他の好適な例では、フーリエ変換後の信号の積分は信号
処理視点からのフーリエ変換に等価の演算、即ち、重み
付き関数（１／ｔ）に等しいか又はこれから取出したフ
ィルタ関数を有するフィルタリングと置換し得るように
する。
処理視点からのフーリエ変換に等価の演算、即ち、重み
付き関数（１／ｔ）に等しいか又はこれから取出したフ
ィルタ関数を有するフィルタリングと置換し得るように
する。
本発明はスピン共鳴信号が横方向弛緩時定数Ｔ２により
影響を受けると云う事実を用いる。まず最初、絶対信号
値はｅｘｐ（−ＴＥ／Ｔ２）に比例する。第２に、この
横方向弛緩時定数Ｔ２のため、サンプル化したスピン共
鳴信号は非対称となる。この非対称性からスピン密度ρ
（ｘ、 Ｙ＋ ｚ）の像を決めるに充分な情報を抽出し
得ると共に同一のＭＲＩ測定から横方向弛緩時定数Ｔ２
（ＸＩ Ｌ ｚ）の像を決めるに充分な情報を抽出し得
るようにする。サンプリング周期中に得られたスピン共
鳴信号のフーリエ変換後、実数および虚数成分の周波数
分布を得るようにする。
影響を受けると云う事実を用いる。まず最初、絶対信号
値はｅｘｐ（−ＴＥ／Ｔ２）に比例する。第２に、この
横方向弛緩時定数Ｔ２のため、サンプル化したスピン共
鳴信号は非対称となる。この非対称性からスピン密度ρ
（ｘ、 Ｙ＋ ｚ）の像を決めるに充分な情報を抽出し
得ると共に同一のＭＲＩ測定から横方向弛緩時定数Ｔ２
（ＸＩ Ｌ ｚ）の像を決めるに充分な情報を抽出し得
るようにする。サンプリング周期中に得られたスピン共
鳴信号のフーリエ変換後、実数および虚数成分の周波数
分布を得るようにする。
スピン共鳴信号を正しい位相でサンプル化するものとす
ると、この周波数分布の実数部分は次式で表わすことが
できる。
ると、この周波数分布の実数部分は次式で表わすことが
できる。
ｐ （ｘ、ｙ、ｚ）＊ｅｘｐ（−ＴＥ／Ｔ２（ｘＸｙ、
ｚ））通常、虚数部分は無視するか、またはモジュラス
像を計算する場合には、この項の影響はその説明のため
無視し得るようにする。しかし、これは、横方向弛緩時
定数Ｔ２が同等影響を与えない場合にのみ完全に正しい
ものである。
ｚ））通常、虚数部分は無視するか、またはモジュラス
像を計算する場合には、この項の影響はその説明のため
無視し得るようにする。しかし、これは、横方向弛緩時
定数Ｔ２が同等影響を与えない場合にのみ完全に正しい
ものである。
その理由はかかる変換された信号の虚数部分を次式で表
わし得るからである。
わし得るからである。
ｄ／ｄｘ （Ｃ＊ｏ （ｘ、ｙ、ｚ）＊１／Ｔ２（ｘ、
ｙ、ｚ）＊ｅｘｐ（−ＴＥ／Ｔ２（ｘ、ｙ、ｚ）））即
ち、Ｘの導関数として、読出し傾度の強度に依存する定
数Ｃと、検出された核のジャイロマグネチック比、スピ
ン密度ρ（ＸＪＹ％Ｚ）、重み係数ｅＸｐ（−ＴＥ／Ｔ
２（ｘ、ｙ、ｚ））、横方向弛緩時定数Ｔ２の（Ｘ％Ｙ
％Ｚ）の逆数の／積の読出し傾度の方向に関連する座標
である。この虚数部分が座標Ｘの関数として積分される
場合には次式で示す信号が発生する。
ｙ、ｚ）＊ｅｘｐ（−ＴＥ／Ｔ２（ｘ、ｙ、ｚ）））即
ち、Ｘの導関数として、読出し傾度の強度に依存する定
数Ｃと、検出された核のジャイロマグネチック比、スピ
ン密度ρ（ＸＪＹ％Ｚ）、重み係数ｅＸｐ（−ＴＥ／Ｔ
２（ｘ、ｙ、ｚ））、横方向弛緩時定数Ｔ２の（Ｘ％Ｙ
％Ｚ）の逆数の／積の読出し傾度の方向に関連する座標
である。この虚数部分が座標Ｘの関数として積分される
場合には次式で示す信号が発生する。
０ （ｘ、ｙ、ｚ）＊１／Ｔ２（ｘＳｙ、ｚ）＊ｅｘｐ
（−ＴＥ／Ｔ２（ｘ、ｙ、ｚ））従ってフーリエ変換の
結果の実数部分はこれにより除算され、この商は横方向
弛緩時定数Ｔ２（ｘ、ｙ、ｚ）の極めて正確な測定値と
なる。この定数が既知になると、これから、フーリエ変
換されたスピン共鳴信号の実数部と相俟ってスピン密度
ρ（Ｘ％Ｙ％Ｚ）を極めて正確に決めることができる。
（−ＴＥ／Ｔ２（ｘ、ｙ、ｚ））従ってフーリエ変換の
結果の実数部分はこれにより除算され、この商は横方向
弛緩時定数Ｔ２（ｘ、ｙ、ｚ）の極めて正確な測定値と
なる。この定数が既知になると、これから、フーリエ変
換されたスピン共鳴信号の実数部と相俟ってスピン密度
ρ（Ｘ％Ｙ％Ｚ）を極めて正確に決めることができる。
本発明の好適な変形例では、スピン共鳴信号の信号サン
プルのフーリエ変換後、前記変換信号の虚数部分を積分
し、その後実数部分と変換信号の積分虚数部分との商を
決めるようにし、かつ、横方向弛緩のない場合に、即ち
、初期の長さＴ２に対し信号のフーリエ変換が純粋に実
数となるような位相のスピン共鳴信号を得るようにする
。
プルのフーリエ変換後、前記変換信号の虚数部分を積分
し、その後実数部分と変換信号の積分虚数部分との商を
決めるようにし、かつ、横方向弛緩のない場合に、即ち
、初期の長さＴ２に対し信号のフーリエ変換が純粋に実
数となるような位相のスピン共鳴信号を得るようにする
。
本発明方法の更に他の例では、フーリエ変換後、横方向
弛緩のない場合に純粋に実数となる要求をフーリエ変換
信号が満足するまで、スピン共鳴信号の位相の補正を継
続し、位相の補正された信号サンプルのフーリエ変換に
ついで変換信号の虚数部の積分を行い、その後この目的
のために変換された信号の積分虚数部が小さすぎないか
ぎり、変換信号の実数部と積分虚数部との商を決め得る
ようにする。前記変換信号の積分虚数部分が前記間を決
めるには小さすぎる画素に対し、前記間の決定を周囲の
画素に基づく推定または結果を零に等置するような他の
演算で置換し得るようにする。
弛緩のない場合に純粋に実数となる要求をフーリエ変換
信号が満足するまで、スピン共鳴信号の位相の補正を継
続し、位相の補正された信号サンプルのフーリエ変換に
ついで変換信号の虚数部の積分を行い、その後この目的
のために変換された信号の積分虚数部が小さすぎないか
ぎり、変換信号の実数部と積分虚数部との商を決め得る
ようにする。前記変換信号の積分虚数部分が前記間を決
めるには小さすぎる画素に対し、前記間の決定を周囲の
画素に基づく推定または結果を零に等置するような他の
演算で置換し得るようにする。
本発明方法の更に他の例では、フーリエ変換後の信号の
積分は信号処理視点からのフーリエ変換に等価の演算、
即ち、重み付き関数に等しいか又はこれから取出したフ
ィルタ関数を有するフィルタリングと置換し得るように
する。
積分は信号処理視点からのフーリエ変換に等価の演算、
即ち、重み付き関数に等しいか又はこれから取出したフ
ィルタ関数を有するフィルタリングと置換し得るように
する。
本発明の好適な変形例では、前記スピン共鳴信号を傾斜
が第１方向にある磁場の存在下でサンプル化し、この際
、フーリエ変換後横方向弛緩の存在しない場合に純粋な
実数結果が得られるように前記スピン共鳴信号をサンプ
ル化するかまたは補正し、横方向弛緩時定数Ｔ２の分布
の決定はかくして得た信号のフーリエ変換の虚数成分を
上記方向に積分して行い、各画素に対し当該画素に関連
するフーリエ変換の実数成分を電流積分によって除算し
得るようにする。
が第１方向にある磁場の存在下でサンプル化し、この際
、フーリエ変換後横方向弛緩の存在しない場合に純粋な
実数結果が得られるように前記スピン共鳴信号をサンプ
ル化するかまたは補正し、横方向弛緩時定数Ｔ２の分布
の決定はかくして得た信号のフーリエ変換の虚数成分を
上記方向に積分して行い、各画素に対し当該画素に関連
するフーリエ変換の実数成分を電流積分によって除算し
得るようにする。
この方法によって決められた横方向弛緩時定数Ｔ２の分
布を用いることによりフーリエ変換の実数部分に関連す
る像を補正して同一の解に対し、横方向弛緩時定数Ｔ”
の分布が決まるものと同一の身体の部分に対して純粋の
スピン密度分布ρを見いだし得るようにする。
布を用いることによりフーリエ変換の実数部分に関連す
る像を補正して同一の解に対し、横方向弛緩時定数Ｔ”
の分布が決まるものと同一の身体の部分に対して純粋の
スピン密度分布ρを見いだし得るようにする。
また、本発明は均一静磁場を発生すると共に傾斜磁場を
発生する磁石系を具え、かつスピン共鳴信号を発生しか
つ検出するＲＦ送信器および受信器と、発生したスピン
共鳴信号をサンプリングするサンプリング装置と、サン
プル化したスピン共鳴信号を変換する信号変換装置と、
変換した信号を表示する表示装置と、前記磁石系、ＲＦ
送信器および受信器と、サンプリング装置と、信号変換
装置と、表示装置とを制御する制御装置とを具えるＭＲ
Ｉ装置において、前記信号変換装置はサンプル化された
スピン共鳴信号の信号変換の虚数部分を積分すると共に
信号変換および積分部分の実数成分の商を決める算術演
算手段を具える弧とを特徴とする。
発生する磁石系を具え、かつスピン共鳴信号を発生しか
つ検出するＲＦ送信器および受信器と、発生したスピン
共鳴信号をサンプリングするサンプリング装置と、サン
プル化したスピン共鳴信号を変換する信号変換装置と、
変換した信号を表示する表示装置と、前記磁石系、ＲＦ
送信器および受信器と、サンプリング装置と、信号変換
装置と、表示装置とを制御する制御装置とを具えるＭＲ
Ｉ装置において、前記信号変換装置はサンプル化された
スピン共鳴信号の信号変換の虚数部分を積分すると共に
信号変換および積分部分の実数成分の商を決める算術演
算手段を具える弧とを特徴とする。
本発明装置では、スピン密度分布および横方向弛緩時定
数の分布をソフトウェアおよび／またはハードウェアの
最小量のみを加えることによって、信号測定から、即ち
、単一に得られる像層波数マトリックスから決めること
ができる。
数の分布をソフトウェアおよび／またはハードウェアの
最小量のみを加えることによって、信号測定から、即ち
、単一に得られる像層波数マトリックスから決めること
ができる。
本発明ＭＲＩ装置の好適な例では、前記信号変換装置は
これによってフーリエ変換を行うと共に前記スピン共鳴
信号の実数成分および記憶する記憶手段を具え、他に虚
数部分の現在の積分値を記憶し、および／またはフーリ
エ変換の実数成分およびこの現在の積分の関連する値の
商を記憶する記憶手段を具えるようにする。このＭＲＩ
装置の好適な例によれば、１つの画像に関連する記憶手
段および画素当たり、原理的には、横方向弛緩時定数Ｔ
２の分布を決める２つの計算のみを追加するだけで充分
である。
これによってフーリエ変換を行うと共に前記スピン共鳴
信号の実数成分および記憶する記憶手段を具え、他に虚
数部分の現在の積分値を記憶し、および／またはフーリ
エ変換の実数成分およびこの現在の積分の関連する値の
商を記憶する記憶手段を具えるようにする。このＭＲＩ
装置の好適な例によれば、１つの画像に関連する記憶手
段および画素当たり、原理的には、横方向弛緩時定数Ｔ
２の分布を決める２つの計算のみを追加するだけで充分
である。
（実施例）
図面につき本発明を説明する。
第１図に示す本発明ＭＲＩ装置ｌは均一静磁場および傾
斜磁場を発生する磁石系ｌＯを具える。かかる磁場は検
査すべき物体、例えば、患者を位置させる検査スペース
１１に発生させるようにする。磁石系１０は発生器１２
および１３により駆動し、この際発生器１２によって静
磁場を発生させるが、発生器１３によって任意の方向を
有する傾斜磁場を発生させるようにする。均一の静磁場
を発生させるためには、磁石系１０には発生器１２によ
り駆動される抵抗性コイルまたは超伝導コイルを設ける
ようにする。
斜磁場を発生する磁石系ｌＯを具える。かかる磁場は検
査すべき物体、例えば、患者を位置させる検査スペース
１１に発生させるようにする。磁石系１０は発生器１２
および１３により駆動し、この際発生器１２によって静
磁場を発生させるが、発生器１３によって任意の方向を
有する傾斜磁場を発生させるようにする。均一の静磁場
を発生させるためには、磁石系１０には発生器１２によ
り駆動される抵抗性コイルまたは超伝導コイルを設ける
ようにする。
また、ＭＲＩ装置１には送／受信器１４を設け、これに
より送／受信コイル１５を駆動し得るようにする。
より送／受信コイル１５を駆動し得るようにする。
送信コイル１５の作動後検査スペース１１に位置する物
体にスピン共鳴信号を発生し、この信号を受信コイル１
５によって検出すると共に送／受信器１４の受信器に供
給する。実行すべき測定技術に応じて同一のコイルを送
信コイルおよび受信コイルの双方に用いることができ、
また、送受信に別個のコイルを用いることもできる。検
出したスピン共鳴信号はサンプル化し、ディジタル化し
て中央処理ユニット１６に供給する。受信器１４による
スピン共鳴信号の発生、検出およびディジタル化につい
ては後に詳細に説明する。中央処理ユニット１６は信号
変換装置１６Ａおよび制御兼調整ユニッｌ−１６Ｂを具
え、これにより発生器１２および１３並びに送／受信器
１４を制御すると共にこの作動を同期化し得るようにす
る。信号変換装置１６Ａによって測定したスピン共鳴信
号を処理して物体の励起部分の像を得るようにする。ま
た、再構成した像を表示するために、例えばモニタのよ
うな表示装置１７を設ける。
体にスピン共鳴信号を発生し、この信号を受信コイル１
５によって検出すると共に送／受信器１４の受信器に供
給する。実行すべき測定技術に応じて同一のコイルを送
信コイルおよび受信コイルの双方に用いることができ、
また、送受信に別個のコイルを用いることもできる。検
出したスピン共鳴信号はサンプル化し、ディジタル化し
て中央処理ユニット１６に供給する。受信器１４による
スピン共鳴信号の発生、検出およびディジタル化につい
ては後に詳細に説明する。中央処理ユニット１６は信号
変換装置１６Ａおよび制御兼調整ユニッｌ−１６Ｂを具
え、これにより発生器１２および１３並びに送／受信器
１４を制御すると共にこの作動を同期化し得るようにす
る。信号変換装置１６Ａによって測定したスピン共鳴信
号を処理して物体の励起部分の像を得るようにする。ま
た、再構成した像を表示するために、例えばモニタのよ
うな表示装置１７を設ける。
第２図はスピンエコ一方法と称される既知の測定方法を
示す。このスピンエコ一方法は後述するように、本発明
による本発明の一部分を示す。瞬時１＝ａ前では、制御
手段１６ｂにより発生器１３を作動させて傾斜磁場Ｇｚ
を発生し得るようにし、従って制御手段１６ｂが送信器
１４を瞬時１＝０で作動させて、通常９０°パルスであ
るＲＦパルスＰａを発生させる場合にはスピンが任意の
方向に励起されて検査スペース１１に配置された物体を
選択スライスし得るようにする。ＲＦパルスＰ１の終了
後で、次のＲＦパルスＰ２の発生前に傾斜Ｇｙを有する
第２磁場を供給すると共に傾斜Ｇｘを有する第３磁場を
も供給し、傾斜磁場ＧＹの時間積分が各測定サイクルに
対し相違し得るようにする。既知のように、傾斜Ｇｙを
有する磁場の時間積分を用いて像層波数ｋＹを選択する
と共にこの周波数ｋｙの各値に対しスピン共鳴信号を像
層波数ｋｙの関数として測定し得るようにする。
示す。このスピンエコ一方法は後述するように、本発明
による本発明の一部分を示す。瞬時１＝ａ前では、制御
手段１６ｂにより発生器１３を作動させて傾斜磁場Ｇｚ
を発生し得るようにし、従って制御手段１６ｂが送信器
１４を瞬時１＝０で作動させて、通常９０°パルスであ
るＲＦパルスＰａを発生させる場合にはスピンが任意の
方向に励起されて検査スペース１１に配置された物体を
選択スライスし得るようにする。ＲＦパルスＰ１の終了
後で、次のＲＦパルスＰ２の発生前に傾斜Ｇｙを有する
第２磁場を供給すると共に傾斜Ｇｘを有する第３磁場を
も供給し、傾斜磁場ＧＹの時間積分が各測定サイクルに
対し相違し得るようにする。既知のように、傾斜Ｇｙを
有する磁場の時間積分を用いて像層波数ｋＹを選択する
と共にこの周波数ｋｙの各値に対しスピン共鳴信号を像
層波数ｋｙの関数として測定し得るようにする。
ここに、ｋｙは傾斜Ｇｘを有する磁場の時間積分によっ
て決めるものとする。瞬時ｔｌではエコーパルスと称さ
れ２通常１８０°パルスであるＲＦパルスＰ２を発生し
、従って、エコー信号Ｅがほぼ瞬時ｔ２に発生し、瞬時
ｔ＝ａおよびＬ２間に経過する時間周期および瞬時ｔ１
およびｔ２間に経過する時間周期は等しくなり、このエ
コー信号を傾斜磁場Ｇｘの存在下でサンプル化し得るよ
うにする。この信号Ｅの振幅はスピン密度ρおよび弛緩
時定数Ｔ１およびＴ２の関数である。
て決めるものとする。瞬時ｔｌではエコーパルスと称さ
れ２通常１８０°パルスであるＲＦパルスＰ２を発生し
、従って、エコー信号Ｅがほぼ瞬時ｔ２に発生し、瞬時
ｔ＝ａおよびＬ２間に経過する時間周期および瞬時ｔ１
およびｔ２間に経過する時間周期は等しくなり、このエ
コー信号を傾斜磁場Ｇｘの存在下でサンプル化し得るよ
うにする。この信号Ｅの振幅はスピン密度ρおよび弛緩
時定数Ｔ１およびＴ２の関数である。
エコー信号Ｅの信号サンプルを取出すサンプリング周期
Ｔｍ後、新たなエコー信号を次の１８０°パルスによっ
て発生し、その後サンプリングし得るようにする。更に
、傾斜磁場ａｙを上述したように進行方向に一時的に作
動させる場合には、充分な数の信号サンプルを取出すに
要する時間周期を有利に減少させるようにする。この方
法も既知である。
Ｔｍ後、新たなエコー信号を次の１８０°パルスによっ
て発生し、その後サンプリングし得るようにする。更に
、傾斜磁場ａｙを上述したように進行方向に一時的に作
動させる場合には、充分な数の信号サンプルを取出すに
要する時間周期を有利に減少させるようにする。この方
法も既知である。
実際上、上述した方法はパルスエコーまたはフィールド
エコーが発生する場合の測定方法の全部に対し用いるこ
とができる。かかる測定方法により発生したエコー信号
は関数ｅｘｐ（−ｔ／Ｔ２）の時間の関数として重み付
けされるようになる。第２図はこの重み係数を関数Ｗと
して示す。
エコーが発生する場合の測定方法の全部に対し用いるこ
とができる。かかる測定方法により発生したエコー信号
は関数ｅｘｐ（−ｔ／Ｔ２）の時間の関数として重み付
けされるようになる。第２図はこの重み係数を関数Ｗと
して示す。
第３図のフローチャートは第２図に示す測定方法によっ
て得られた信号サンプルの処理を示す。
て得られた信号サンプルの処理を示す。
これら信号サンプルは、個別の傾斜磁場ＧｘおよびＧＶ
の強度の時間積分によって決まる像層波数ｋｘおよびｋ
Ｙの関数として得ることができ、この時間積分は信号サ
ンプルを取出す瞬時まで励起パルスから延在し得るよう
にする。
の強度の時間積分によって決まる像層波数ｋｘおよびｋ
Ｙの関数として得ることができ、この時間積分は信号サ
ンプルを取出す瞬時まで励起パルスから延在し得るよう
にする。
フローチャートの第１ステツプ１００では、信号サンプ
ル５（ｋｘ、ｋｙ）の位相を、所望に応じ既知の式ｅｘ
ｐ（ｉφ）で表すように補正する。
ル５（ｋｘ、ｋｙ）の位相を、所望に応じ既知の式ｅｘ
ｐ（ｉφ）で表すように補正する。
フローチャートの第２ステツプ２００に示すように、分
布Ｓ（Ｘ％Ｙ）はかくして得た像層波数マド１ルックス
５２（ｋｘ、ｋｙ）からのフーリエ変換によって決まる
。横方向弛緩時定数Ｔ２が有効でない場合にはスピン密
度ρ（ｘ、 ｙ）は測定されるようになる。その理由は
５（ｋｘ、 ｋｙ）および密度ρ（ｘＳｙ）がフーリエ
対を形成するからである。実際上、時定数Ｔ２は常時有
効である場合には得られた像層波数マトリックス５（ｋ
ｘ、ｋｙ）は密度ｐ （ｘ、ｙ）及び重み関数Ｗ＝ｅｘ
ｐ（−ｔ／Ｔ２（ｘ、ｙ））のフーリエ変換の積となる
。
布Ｓ（Ｘ％Ｙ）はかくして得た像層波数マド１ルックス
５２（ｋｘ、ｋｙ）からのフーリエ変換によって決まる
。横方向弛緩時定数Ｔ２が有効でない場合にはスピン密
度ρ（ｘ、 ｙ）は測定されるようになる。その理由は
５（ｋｘ、 ｋｙ）および密度ρ（ｘＳｙ）がフーリエ
対を形成するからである。実際上、時定数Ｔ２は常時有
効である場合には得られた像層波数マトリックス５（ｋ
ｘ、ｋｙ）は密度ｐ （ｘ、ｙ）及び重み関数Ｗ＝ｅｘ
ｐ（−ｔ／Ｔ２（ｘ、ｙ））のフーリエ変換の積となる
。
この重み関数Ｗはｅｘｐ（（ｒ ＋ＴＥ／Ｔ２））とし
て表わすことができる。ここにＴＥは瞬時１＝０におけ
るパルスＰ１とエコー信号Ｅの中心（−ｔ２）との間の
時間経過であり、その双方を第２図に示す。時間τは瞬
時ｔ＝ＴＥからとるため、これは正および負の値になる
ものとする。かくして得た重み関数Ｗはτの関数として
、対称部分ＳＤおよび非対称部分ＡＤに分解することが
できる：即ち、Ｗ＝ＳＤ＋＾Ｄ、ここに５Ｄ＝ｅｘｐ（
−ＴＥ／Ｔ２）＊ｃｏｓ（ｒ ／Ｔ２）およびＡＤ＝ｅ
ｘｐ（−ＴＥ／Ｔ２）＊５ｉｎ（ｒ／Ｔ２）である。重
み関数の対称部分ＳＤ重み非対称部分ＡＤに対するこれ
らの式は簡単化することができる。その理由はこれによ
り重み付けされた信号がτ−０の外側でＯとなり、５Ｄ
−ｅｘｐ（−ＴＥ／Ｔ２）およびＡＤ＝（−ｒ ／Ｔ２
）＊ｅｘｐ（−ＴＥ／Ｔ２）が得られるからである。
て表わすことができる。ここにＴＥは瞬時１＝０におけ
るパルスＰ１とエコー信号Ｅの中心（−ｔ２）との間の
時間経過であり、その双方を第２図に示す。時間τは瞬
時ｔ＝ＴＥからとるため、これは正および負の値になる
ものとする。かくして得た重み関数Ｗはτの関数として
、対称部分ＳＤおよび非対称部分ＡＤに分解することが
できる：即ち、Ｗ＝ＳＤ＋＾Ｄ、ここに５Ｄ＝ｅｘｐ（
−ＴＥ／Ｔ２）＊ｃｏｓ（ｒ ／Ｔ２）およびＡＤ＝ｅ
ｘｐ（−ＴＥ／Ｔ２）＊５ｉｎ（ｒ／Ｔ２）である。重
み関数の対称部分ＳＤ重み非対称部分ＡＤに対するこれ
らの式は簡単化することができる。その理由はこれによ
り重み付けされた信号がτ−０の外側でＯとなり、５Ｄ
−ｅｘｐ（−ＴＥ／Ｔ２）およびＡＤ＝（−ｒ ／Ｔ２
）＊ｅｘｐ（−ＴＥ／Ｔ２）が得られるからである。
密度ρの所定の分布に対しては、ステップ２００で得ら
れた信号Ｓ（ｘ、ｙ）の実数部Ａは重み関数Ｗの対称部
分ＳＤによってのみ決まるようになる；即ち、この実数
部Ａは密度ρ＊ｅｘｐ（−ＴＥ／Ｔ２）に等しくなる。
れた信号Ｓ（ｘ、ｙ）の実数部Ａは重み関数Ｗの対称部
分ＳＤによってのみ決まるようになる；即ち、この実数
部Ａは密度ρ＊ｅｘｐ（−ＴＥ／Ｔ２）に等しくなる。
信号Ｓ（ｘ、 ｙ）の虚数部Ｂは密度ρの所定の分布に
対し重み関数Ｗの非対称部分へ〇によってのみ決まる。
対し重み関数Ｗの非対称部分へ〇によってのみ決まる。
ＡＤに対し最後に与えられた式からこの虚数部Ｂがｄ／
ｄｘ （Ｃ＊ｊ／Ｔ２ｔｅｘｐ（−ＴＥ／Ｔ２）に等し
くなることを導くことができる。ここにＣは読出し傾斜
の大きさおよび検出された核のジャイロマグネチック比
に依存する定数、Ｘは読出し傾斜の方向に関連する座標
である。
ｄｘ （Ｃ＊ｊ／Ｔ２ｔｅｘｐ（−ＴＥ／Ｔ２）に等し
くなることを導くことができる。ここにＣは読出し傾斜
の大きさおよび検出された核のジャイロマグネチック比
に依存する定数、Ｘは読出し傾斜の方向に関連する座標
である。
フローチャートの第３ステツプ３００に従って、信号Ｓ
（ｘ、 ｙ）の虚数部Ｂを座標Ｘに対して積分し、かつ
係数−Ｃで除算する。信号Ｂの現在の積分値によって、
式ｐ ／Ｔ２＊ｅｘｐ（−ＴＥ／Ｔ２）により示される
新たなイ言号１ （ｘ、ｙ）を形成する。
（ｘ、 ｙ）の虚数部Ｂを座標Ｘに対して積分し、かつ
係数−Ｃで除算する。信号Ｂの現在の積分値によって、
式ｐ ／Ｔ２＊ｅｘｐ（−ＴＥ／Ｔ２）により示される
新たなイ言号１ （ｘ、ｙ）を形成する。
フローチャートの第４ステ・ツブ４００（こ示すように
、画素当たりの信号＾（ｘ、ｙ）および積分信号１（ｘ
、ｙ）の商を決めるようにする。これを画素の全部Ｉこ
対して行う。これらの信号に対し与えられた式ｈＡら明
らかなように、この商は横方向弛緩時定数Ｔ２に等しく
なる。信号１（ｘ、ｙ）に対し用０られた武力（正しく
読出されないと云う事実のため、この商とＴ２の実際の
値との間の差は小さくなる。
、画素当たりの信号＾（ｘ、ｙ）および積分信号１（ｘ
、ｙ）の商を決めるようにする。これを画素の全部Ｉこ
対して行う。これらの信号に対し与えられた式ｈＡら明
らかなように、この商は横方向弛緩時定数Ｔ２に等しく
なる。信号１（ｘ、ｙ）に対し用０られた武力（正しく
読出されないと云う事実のため、この商とＴ２の実際の
値との間の差は小さくなる。
かくして得た時定数Ｔ２の像は、次のステ・ノブ５００
中表示装置１７（第１図参照）に表示すること力（でき
る。
中表示装置１７（第１図参照）に表示すること力（でき
る。
第１図は本発明装置の１例の構成を示すブロック図、
第２図は本発明測定サイクルの１例を示す説明図、
第３図は第２図の測定サイクルの信号の処理例を示すフ
ローチャートである。 ■ ・・・ １１ ・・・ １４ ・・・ １５ ・・・ １６ ・・・ １６Ａ ・・・ １６Ｂ ・・・ １７ ・・・ ＭＲＩ装置、 １０ ・・・ 磁石系検査スペース、
１２、Ｉ３ ・・・ 発生器送／受信器 送／受信コイル 中央処理装置 信号変換装置 制御兼調整ユニット 表示装置
ローチャートである。 ■ ・・・ １１ ・・・ １４ ・・・ １５ ・・・ １６ ・・・ １６Ａ ・・・ １６Ｂ ・・・ １７ ・・・ ＭＲＩ装置、 １０ ・・・ 磁石系検査スペース、
１２、Ｉ３ ・・・ 発生器送／受信器 送／受信コイル 中央処理装置 信号変換装置 制御兼調整ユニット 表示装置
Claims (1)
- 【特許請求の範囲】 １、スピン共鳴信号を発生させるために静磁場の存在の
下で身体中で核スピンを励起し、その後静磁場に重畳す
る傾斜磁場の存在の下でスピン共鳴信号をサンプル化し
、かくしてサンプル化したスピン共鳴信号を信号変換に
より核磁化分布像に変換するようにした身体の一部分の
核磁化分布像を得るに当たり、この核磁化分布を横方向
弛緩時定数分布とし、スピンエコー信号を発生する励起
スピンの横方向弛緩により非対称性を生ぜしめるサンプ
ル化スピン共鳴信号の相対的な非対称性から、部分横方
向弛緩時定数の測定値を取出すようにしたことを特徴と
する核磁化分布像を得る方法。 ２、スピン共鳴信号を発生させるために静磁場の存在の
下で身体中で核スピンを励起し、その後静磁場に重畳す
ると共に第１の方向を有する傾斜磁場の存在の下でスピ
ン共鳴信号をサンプル化し、次いで、かくしてサンプル
化したスピン共鳴信号からの信号変換により核磁化分布
像を測定するようにした身体の一部分のスピン密度分布
および横方向弛緩時定数Ｔ２の分布を測定するに当たり
、スピン共鳴信号の信号サンプルのフーリエ変換後、前
記第１方向における変換信号の虚数部分を積分し、その
後実数部分と変換信号の積分虚数部分との商を決めるよ
うにしたことを特徴とするスピン密度分布および横方向
弛緩時定数分布をきめる方法。 ３、フーリエ変換後、信号が横方向弛緩のない場合に純
粋な実数となるような位相を有する信号を得るようにし
たことを特徴とする請求項２に記載のスピン密度分布お
よび横方向弛緩時定数分布をきめる方法。 ４、前記変換信号の積分虚数部分が前記商を決めるには
小さすぎる画素に対し、前記商の決定を周囲の画素から
の補間のような他の演算で置換することを特徴とする請
求項２に記載のスピン密度分布および横方向弛緩時定数
分布をきめる方法。 ５、フーリエ変換後の信号の積分は信号処理視点からの
フーリエ変換に等価の演算、即ち、重み付き関数（１／
ｔ）に等しいか又はこれから取出したフィルタ関数を有
するフィルタリングと置換するようにしたことを特徴と
する請求項２に記載のスピン密度分布および横方向弛緩
時定数分布をきめる方法。 ６、スピン共鳴信号の位相をフーリエ変換前又は後に補
正して、横方向弛緩のない場合に前記スピン共鳴信号の
フーリエ変換が純粋に実数となるようにしたことを特徴
とする請求項２に記載のスピン密度分布および横方向弛
緩時定数分布をきめる方法。 ７、各画素に対し決められた商を横方向弛緩時定数Ｔ２
の推測値として用いるようにしたことを特徴とする請求
項２、３、４、５、または６に記載のスピン密度分布お
よび横方向弛緩時定数分布をきめる方法。 ８、前記スピン共鳴信号を傾斜が第１方向にある磁場の
存在下でサンプル化し、この際、フーリエ変換後横方向
弛緩の存在しない場合に純粋な実数結果が得られるよう
に前記スピン共鳴信号をサンプル化するかまたは補正し
、横方向弛緩時定数Ｔ２の分布の決定はかくして得た信
号のフーリエ変換の虚数成分を上記方向に積分して行い
、各画素に対し当該画素に関連するフーリエ変換の実数
成分を電流積分によって除算するようにしたことを特徴
とする請求項１〜７の何れかの項に記載のスピン密度分
布および横方向弛緩時定数分布をきめる方法。 ９、均一静磁場を発生すると共に傾斜磁場を発生する磁
石系を具え、かつスピン共鳴信号を発生しかつ検出する
ＲＦ送信器および受信器と、発生したスピン共鳴信号を
サンプリングするサンプリング装置と、サンプル化した
スピン共鳴信号を変換する信号変換装置と、変換した信
号を表示する表示装置と、前記磁石系、ＲＦ送信器およ
び受信器と、サンプリング装置と、信号変換装置と、表
示装置とを制御する制御装置とを具えるＭＲＩ装置にお
いて、前記信号変換装置はサンプル化されたスピン共鳴
信号の信号変換の虚数部分を積分すると共に信号変換お
よび積分部分の実数成分の商を決める算術演算手段を具
えることを特徴とするＭＲＩ装置。 １０、前記信号変換装置はこれによってフーリエ変換を
行うと共に前記スピン共鳴信号の実数成分および記憶す
る記憶手段を具え、他に虚数部分の現在の積分値を記憶
し、および／またはフーリエ変換の実数成分およびこの
現在の積分の関連する値の商を記憶する記憶手段を具え
ることを特徴とする請求項９に記載のＭＲＩ装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| NL8801731A NL8801731A (nl) | 1988-07-08 | 1988-07-08 | Werkwijze en kernspinresonantie-inrichting om de dwarsrelaxatietijdkonstante t2 snel te bepalen. |
| NL8801731 | 1988-07-08 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0265843A true JPH0265843A (ja) | 1990-03-06 |
Family
ID=19852593
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1174335A Pending JPH0265843A (ja) | 1988-07-08 | 1989-07-07 | 核磁化分布像を得る方法 |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US5068610A (ja) |
| EP (1) | EP0350120B1 (ja) |
| JP (1) | JPH0265843A (ja) |
| CN (1) | CN1039484A (ja) |
| DE (1) | DE58907744D1 (ja) |
| FI (1) | FI893275A7 (ja) |
| IL (1) | IL90877A0 (ja) |
| NL (1) | NL8801731A (ja) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0525858A2 (en) * | 1991-07-17 | 1993-02-03 | Koninklijke Philips Electronics N.V. | Method and device for determining variations in the position of a movable object by means of ultrasound |
| DE19619572A1 (de) * | 1996-05-15 | 1997-12-18 | Klaus Lehmann | Verfahren zur Detektion von Unstetigkeiten in Meßsignalen |
| DE19930707C2 (de) * | 1999-07-02 | 2003-04-10 | Forschungszentrum Juelich Gmbh | Meßverfahren, Meßvorrichtung sowie Auswerteelektronik |
| US6898455B2 (en) * | 1999-10-29 | 2005-05-24 | The Mclean Hospital Corporation | Method for providing optimal drug dosage |
| US6400978B1 (en) | 1999-10-29 | 2002-06-04 | The Mclean Hospital Corporation | Method and apparatus for detecting mental disorders |
| US6804384B2 (en) | 2001-06-12 | 2004-10-12 | Mclean Hospital Corporation | Color magnetic resonance imaging |
| WO2007051224A1 (en) * | 2005-11-01 | 2007-05-10 | Resonance Health Analysis Services Pty Ltd | A method of characterising an area of interest within a body |
| JP5113761B2 (ja) * | 2005-11-27 | 2013-01-09 | アキュイタス・メディカル・リミテッド | 空間周波数解析を使用する構造評価 |
| DE102006012405A1 (de) * | 2006-03-17 | 2007-09-27 | Siemens Ag | Verfahren und Referenzsignalisierungsanordnung zur Versorgung einer Mehrzahl von Systemkomponenten mit einem gemeinsamen Referenzsignal |
| WO2011114264A1 (en) * | 2010-03-18 | 2011-09-22 | Koninklijke Philips Electronics N.V. | Simultaneous and dynamic determination of longitudinal and transversal relaxation times of a nuclear spin system |
| US11122988B2 (en) * | 2015-08-06 | 2021-09-21 | Hitachi, Ltd. | Magnetic resonance imaging apparatus |
| US11287501B1 (en) * | 2021-02-25 | 2022-03-29 | Wisconsin Alumni Research Foundation | System and method for multiple contrast encoding with magnetic resonance imaging |
| CN115993562B (zh) * | 2022-12-21 | 2026-02-24 | 上海电气控股集团有限公司智惠医疗装备分公司 | 确定弛豫时间的方法、装置、设备及存储介质 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4521733A (en) * | 1983-05-23 | 1985-06-04 | General Electric Company | NMR Imaging of the transverse relaxation time using multiple spin echo sequences |
| US4573014A (en) * | 1983-11-09 | 1986-02-25 | Duke University | NMR Imaging method and apparatus |
| DE3514530A1 (de) * | 1985-04-22 | 1986-10-23 | Siemens AG, 1000 Berlin und 8000 München | Kernspinresonanzgeraet |
| US4720679A (en) * | 1985-12-31 | 1988-01-19 | Picker International, Inc. | Magnetic resonance imaging with phase encoded chemical shift correction |
| IL80727A (en) * | 1986-11-21 | 1990-06-10 | Elscint Ltd | Noise reduction in magnetic resonance images |
| US4866386A (en) * | 1988-04-26 | 1989-09-12 | Picker International, Inc. | High signal-to-noise, rapid calibration of RF pulses |
-
1988
- 1988-07-08 NL NL8801731A patent/NL8801731A/nl not_active Application Discontinuation
-
1989
- 1989-05-26 US US07/358,327 patent/US5068610A/en not_active Expired - Fee Related
- 1989-07-03 DE DE58907744T patent/DE58907744D1/de not_active Expired - Fee Related
- 1989-07-03 EP EP89201752A patent/EP0350120B1/de not_active Expired - Lifetime
- 1989-07-05 FI FI893275A patent/FI893275A7/fi not_active Application Discontinuation
- 1989-07-05 CN CN89104674.7A patent/CN1039484A/zh active Pending
- 1989-07-05 IL IL90877A patent/IL90877A0/xx unknown
- 1989-07-07 JP JP1174335A patent/JPH0265843A/ja active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| EP0350120A1 (de) | 1990-01-10 |
| CN1039484A (zh) | 1990-02-07 |
| FI893275L (fi) | 1990-01-09 |
| US5068610A (en) | 1991-11-26 |
| FI893275A0 (fi) | 1989-07-05 |
| EP0350120B1 (de) | 1994-06-01 |
| DE58907744D1 (de) | 1994-07-07 |
| IL90877A0 (en) | 1990-02-09 |
| FI893275A7 (fi) | 1990-01-09 |
| NL8801731A (nl) | 1990-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3402916B2 (ja) | 核スピントモグラフ装置のマグネット装置のシム調整方法及び方法を実施する装置 | |
| JP4201993B2 (ja) | Ｎｍｒイメージング操作中に誘導される渦電流の測定および補償のための装置及びシステム | |
| JP4672826B2 (ja) | 核スピン共鳴装置における渦電流を検出するための方法 | |
| JP2716889B2 (ja) | 磁石の高速補正法 | |
| EP3060116B1 (en) | Mr imaging with temperature mapping | |
| EP0132358B1 (en) | Imaging method and apparatus for obtaining images using nmr | |
| US4520828A (en) | Nuclear magnetic resonance method and apparatus | |
| EP3295199B1 (en) | Magnetic resonance examination system with field probes | |
| US4553096A (en) | Nuclear magnetic resonance method and apparatus | |
| US5860921A (en) | Method for measuring the reversible contribution to the transverse relaxation rate in magnetic resonance imaging | |
| US5068610A (en) | Mri method and device for fast determination of the transverse relaxation time constant t2 | |
| US5451877A (en) | Method for the compensation of eddy currents caused by gradients in a nuclear magnetic resonance apparatus | |
| US6127826A (en) | EPI image based long term eddy current pre-emphasis calibration | |
| US4684892A (en) | Nuclear magnetic resonance apparatus | |
| US20020050816A1 (en) | MR imaging method, phase error measuring method, and MRI apparatus | |
| EP0106472B1 (en) | Nuclear magnetic resonance method and apparatus | |
| EP4139699B1 (en) | Mr imaging with t1 compensated b1 mapping | |
| WO2006040866A1 (ja) | 磁気共鳴撮影装置及び磁気共鳴撮影方法 | |
| JP2528864B2 (ja) | 核磁気共鳴を用いた検査装置 | |
| GB2127155A (en) | Flow determination by nuclear magnetic resonance | |
| JPH08215167A (ja) | 空間磁場分布測定方法及び空間磁場分布測定装置 | |
| JPS6151581A (ja) | 空間磁場分布測定方法 | |
| JPS61196146A (ja) | 核磁気共鳴画像装置 | |
| JPH038215B2 (ja) | ||
| Zur et al. | MR angiography without subtraction |