JPH0654831A - 磁気共鳴機能イメージング装置 - Google Patents
磁気共鳴機能イメージング装置Info
- Publication number
- JPH0654831A JPH0654831A JP4212588A JP21258892A JPH0654831A JP H0654831 A JPH0654831 A JP H0654831A JP 4212588 A JP4212588 A JP 4212588A JP 21258892 A JP21258892 A JP 21258892A JP H0654831 A JPH0654831 A JP H0654831A
- Authority
- JP
- Japan
- Prior art keywords
- image
- magnetic resonance
- resonance imaging
- phase
- imaging apparatus
- 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
【０００１】
【産業上の利用分野】本発明は核磁気共鳴現象を利用し
た体内断層撮像装置ＭＲＩに関し、特に組織の磁化率の
時間的変化を映像化する手法に関するもので、脳におけ
る活動部位を観察する機能イメージング等に利用する。
た体内断層撮像装置ＭＲＩに関し、特に組織の磁化率の
時間的変化を映像化する手法に関するもので、脳におけ
る活動部位を観察する機能イメージング等に利用する。
【０００２】
【従来の技術】ＭＲＩを利用した脳の機能イメージング
手法の代表的なものに次のようなものがある。
手法の代表的なものに次のようなものがある。
【０００３】（１）造影剤を使用せず、ディフュージョ
ン強調パルスの強度を変化させて、ディフュージョン／
パフュージョン画像を分離して得る方法。
ン強調パルスの強度を変化させて、ディフュージョン／
パフュージョン画像を分離して得る方法。
【０００４】（２）造影剤を使用せず、ディフュージョ
ン強調画像を得る方法。
ン強調画像を得る方法。
【０００５】（３）造影剤を使用して、Ｔ２ 強調画像
を高速撮影シーケンスで得る方法。
を高速撮影シーケンスで得る方法。
【０００６】それぞれの代表的な手法の文献を下記に示
す。
す。
【０００７】（１）D.Le Bihan,E.Breton,D.Lallemand,
M.L.Aubin,J.Vignaud,and M Laval-Jeantet,'Separatio
n of Diffusion and Perfusion in Intravoxel Incoher
ent Motion MR Imaging', Radiolgy 168,pp.497-505(19
88) （２）M.E.Moseley,Y.Cohen,J.Mintorovitch,J.Kucharc
rzyk,P.R.Weinstein,andD.Norman,'Diffusion-weighted
MR Imaging in Experimental Brain Research', Sylla
bus Special Course:MR 1990,pp.63-67(1990) （３）B.R.Rosen,J.W.Belliveau,D.Chien,M.S.Cohen,an
d R.M.Weisskopf,'MR Perfusion Imaging',Syllabus Sp
ecial Course:MR 1990,pp.69-84(1990) 従来の手法には、ディフュージョン強調パルス（大きな
フローエンコードパルス）を印加することで局所的な流
れのあるディフュージョン／パフュージョン領域の信号
をディフェーズさせ、黒くコントラストを付ける（１）
（２）の方法と、造影剤を注入し、ディフュージョン／
パフュージョン領域で磁場の不均一性を発生させ、その
効果を強調するＴ２強調撮影法で、黒くコントラストを
発生させる（３）の方法があった。いずれの方法もター
ボＦｌａｓｈや超高速イメージングと組合せ、時間的に
連続した画像をえる工夫がなされている。
M.L.Aubin,J.Vignaud,and M Laval-Jeantet,'Separatio
n of Diffusion and Perfusion in Intravoxel Incoher
ent Motion MR Imaging', Radiolgy 168,pp.497-505(19
88) （２）M.E.Moseley,Y.Cohen,J.Mintorovitch,J.Kucharc
rzyk,P.R.Weinstein,andD.Norman,'Diffusion-weighted
MR Imaging in Experimental Brain Research', Sylla
bus Special Course:MR 1990,pp.63-67(1990) （３）B.R.Rosen,J.W.Belliveau,D.Chien,M.S.Cohen,an
d R.M.Weisskopf,'MR Perfusion Imaging',Syllabus Sp
ecial Course:MR 1990,pp.69-84(1990) 従来の手法には、ディフュージョン強調パルス（大きな
フローエンコードパルス）を印加することで局所的な流
れのあるディフュージョン／パフュージョン領域の信号
をディフェーズさせ、黒くコントラストを付ける（１）
（２）の方法と、造影剤を注入し、ディフュージョン／
パフュージョン領域で磁場の不均一性を発生させ、その
効果を強調するＴ２強調撮影法で、黒くコントラストを
発生させる（３）の方法があった。いずれの方法もター
ボＦｌａｓｈや超高速イメージングと組合せ、時間的に
連続した画像をえる工夫がなされている。
【０００８】しかしながら、（１）（２）の方法は、画
像上で黒く欠落させるために、組織上のコントラストと
の区別が必ずしも明確にできなかった。また、高速撮影
法と組み合わせない場合には、計測時に毎回ディフュー
ジョン強調パルスを印加するために動きに敏感になり、
モーションアーチファクトが発生しやすくなる問題があ
った。さらに（１）の方法は、渦電流の影響を与えやす
いディフュージョン強調パルスの撮影パラメータを変化
させて撮影した後、画像間演算をするために渦電流の影
響を直接受けてしまう問題があった。一方、（３）の方
法は、造影剤を使用するために、患者に苦痛を与えてし
まい、ＭＲＩの無侵襲性の特徴が活かせない問題があっ
た。
像上で黒く欠落させるために、組織上のコントラストと
の区別が必ずしも明確にできなかった。また、高速撮影
法と組み合わせない場合には、計測時に毎回ディフュー
ジョン強調パルスを印加するために動きに敏感になり、
モーションアーチファクトが発生しやすくなる問題があ
った。さらに（１）の方法は、渦電流の影響を与えやす
いディフュージョン強調パルスの撮影パラメータを変化
させて撮影した後、画像間演算をするために渦電流の影
響を直接受けてしまう問題があった。一方、（３）の方
法は、造影剤を使用するために、患者に苦痛を与えてし
まい、ＭＲＩの無侵襲性の特徴が活かせない問題があっ
た。
【０００９】
【発明が解決しようとする課題】本発明の目的は、「デ
ィフュージョン強調パルス」及び、「造影剤」を用いる
こと無くディフュージョン／パフュージョン領域を検出
し、秒オーダで時間的に変化する様子を映像化する機能
イメージングを実現することにある。
ィフュージョン強調パルス」及び、「造影剤」を用いる
こと無くディフュージョン／パフュージョン領域を検出
し、秒オーダで時間的に変化する様子を映像化する機能
イメージングを実現することにある。
【００１０】
【課題を解決するための手段】血中のヘモグロビンが酸
素と結合すると磁化率が変化する事を利用する。酸素が
常磁性体であるために、酸素含有量によって磁化率が変
化する。これは磁場の不均一性の変化として現われる。
ＭＲＩでは磁場の不均一性が存在すると、画像上に濃度
歪や位相歪、空間歪を発生する。従来の手法は画素内の
不均一によって生じるディフェーズ効果で濃度が低下す
ることを利用していたが、本発明では、位相歪に着目し
て機能的な変化を検出するものである。従来から磁化率
の変化を位相歪でマッピングする手法は知られている
が、時間的に変化する過程を映像化することはされてい
なかった。本発明は、位相マッピング技術と高速イメー
ジング技術を組み合わせて、従来不可能であった時間的
に変化する過程を、連続して撮影した画像の位相の変化
から検出しようとするものである。位相の変化部分のみ
取り出し、一般の形態学的情報を表わすＭＲＩ画像に重
ねて表示することで、脳のどの部分で活動が行われてい
るかを視覚的に見ることが可能になる。
素と結合すると磁化率が変化する事を利用する。酸素が
常磁性体であるために、酸素含有量によって磁化率が変
化する。これは磁場の不均一性の変化として現われる。
ＭＲＩでは磁場の不均一性が存在すると、画像上に濃度
歪や位相歪、空間歪を発生する。従来の手法は画素内の
不均一によって生じるディフェーズ効果で濃度が低下す
ることを利用していたが、本発明では、位相歪に着目し
て機能的な変化を検出するものである。従来から磁化率
の変化を位相歪でマッピングする手法は知られている
が、時間的に変化する過程を映像化することはされてい
なかった。本発明は、位相マッピング技術と高速イメー
ジング技術を組み合わせて、従来不可能であった時間的
に変化する過程を、連続して撮影した画像の位相の変化
から検出しようとするものである。位相の変化部分のみ
取り出し、一般の形態学的情報を表わすＭＲＩ画像に重
ねて表示することで、脳のどの部分で活動が行われてい
るかを視覚的に見ることが可能になる。
【００１１】位相変化は磁化率のみならず、装置の不完
全性（渦電流等）や動きそのものでも生じるために磁化
率の変化分と区別する必要がある。このため時間的に連
続して撮影し、ある時点の位相をレファレンスとして他
の画像の位相の補正処理をリアルタイムで実施する必要
がある。この補正処理は、単純に複素数であるレファレ
ンス画像の複素共役をとり、補正したい画像に乗算する
ことで実施できる。画像に変化が無ければ、補正後の画
像の位相は全ての位置でゼロになるが、変化があるとそ
の位置の位相のみが変化して現われる。その結果をその
まま表示して医師が結果を判断しても良いが、位相の変
化した部分のみを抽出し、画像の絶対値情報に色付けし
て重ねあわせてることで、脳のどの部分が活発に活動し
ているかを直接観察することが可能になる。
全性（渦電流等）や動きそのものでも生じるために磁化
率の変化分と区別する必要がある。このため時間的に連
続して撮影し、ある時点の位相をレファレンスとして他
の画像の位相の補正処理をリアルタイムで実施する必要
がある。この補正処理は、単純に複素数であるレファレ
ンス画像の複素共役をとり、補正したい画像に乗算する
ことで実施できる。画像に変化が無ければ、補正後の画
像の位相は全ての位置でゼロになるが、変化があるとそ
の位置の位相のみが変化して現われる。その結果をその
まま表示して医師が結果を判断しても良いが、位相の変
化した部分のみを抽出し、画像の絶対値情報に色付けし
て重ねあわせてることで、脳のどの部分が活発に活動し
ているかを直接観察することが可能になる。
【００１２】
【作用】磁化率χとすると、磁束密度Ｂと磁場の強さＨ
とは次式の関係にある。
とは次式の関係にある。
【００１３】 Ｂ＝（１＋χ）Ｈ （１） 水分の磁化率は、およそ１０ｐｐｍである。酸素の結合
により磁化率が１０％変化すると言われているために、
約１ｐｐｍの磁場不均一の変化が検出されることにな
る。
により磁化率が１０％変化すると言われているために、
約１ｐｐｍの磁場不均一の変化が検出されることにな
る。
【００１４】一方、位相変化は、磁場不均一をΔＢ、励
起から計測までの時間をｔとすると Δθ＝γ・ΔＢ・ｔ （２） ここで、γ：核磁気回転比 であらわされる。たとえば１．５Ｔ場合、１ｐｐｍあた
り Δθ＝２３゜・ｔ （３） となる。グラディエントエコーのエコータイムＴＥが４
ｍｓｅｃの場合の位相変化は、９２゜になる。位相変化
は１８０゜を超えると位相ラップが発生するので、ＴＥ
を４ｍｓｅｃ前後の適当な値に設定することでは、適切
な感度で酸素含有量の変化を計測できることになり、機
能イメージングが達成されることになる。仮に、磁化率
の変化がもっと小さい場合でも、ＴＥを適切に伸ばすこ
とで対応可能である。上記計算は、１．５Ｔ場合である
が、磁場不均一は磁場強度に比例するため、０．５Ｔの
場合にも位相感度が低くなるため、ＴＥはそれに比例し
て伸ばす事で対応させることができる。
起から計測までの時間をｔとすると Δθ＝γ・ΔＢ・ｔ （２） ここで、γ：核磁気回転比 であらわされる。たとえば１．５Ｔ場合、１ｐｐｍあた
り Δθ＝２３゜・ｔ （３） となる。グラディエントエコーのエコータイムＴＥが４
ｍｓｅｃの場合の位相変化は、９２゜になる。位相変化
は１８０゜を超えると位相ラップが発生するので、ＴＥ
を４ｍｓｅｃ前後の適当な値に設定することでは、適切
な感度で酸素含有量の変化を計測できることになり、機
能イメージングが達成されることになる。仮に、磁化率
の変化がもっと小さい場合でも、ＴＥを適切に伸ばすこ
とで対応可能である。上記計算は、１．５Ｔ場合である
が、磁場不均一は磁場強度に比例するため、０．５Ｔの
場合にも位相感度が低くなるため、ＴＥはそれに比例し
て伸ばす事で対応させることができる。
【００１５】また、グラディエントエコーの代わりにい
わゆる超高速イメージングであるエコープラナーを用い
ても良い。この場合、１８０度パルスを併用する場合
は、９０度−１８０度間隔と１８０度−エコー中心間隔
の差を、先ほどの議論のように４ｍｓｅｃ程度ずらすこ
とで対応させる。
わゆる超高速イメージングであるエコープラナーを用い
ても良い。この場合、１８０度パルスを併用する場合
は、９０度−１８０度間隔と１８０度−エコー中心間隔
の差を、先ほどの議論のように４ｍｓｅｃ程度ずらすこ
とで対応させる。
【００１６】
【実施例】以下、本発明の実施例を図面にて詳細に説明
する。
する。
【００１７】図３は、本発明の一実施例であるＭＲＩ装
置の構成の概要を示すブロック図である。図において、
３０１は均一な静磁場を発生させるための静磁場発生
系、３０２はスピンを励起するための高周波磁場を発生
させる励起系、３０３は磁場の強さをｘ，ｙ，ｚ方向
に、それぞれ独立に線形に変化させることが可能な傾斜
磁場発生系、３０４は被検体から発生する電磁波を受信
し検波した後、Ａ／Ｄ変換する計測系、３０５は該計測
系からの計測データとキーボードから入力された撮影条
件に基づいて画像再生に必要な各種演算を行う処理装
置、３０６は再生結果及び撮影条件を表示するＣＲＴ、
３０７は上記構成における各系の動作のコントロール手
順を格納しておくパルスシーケンスファイル、３０８は
キーボードから入力された撮影条件に従って上記パルス
シーケンスファイルからローディングしたパルスシーケ
ンスに基づいて、各装置の動作タイミングをコントロー
ルするシーケンス制御部である。
置の構成の概要を示すブロック図である。図において、
３０１は均一な静磁場を発生させるための静磁場発生
系、３０２はスピンを励起するための高周波磁場を発生
させる励起系、３０３は磁場の強さをｘ，ｙ，ｚ方向
に、それぞれ独立に線形に変化させることが可能な傾斜
磁場発生系、３０４は被検体から発生する電磁波を受信
し検波した後、Ａ／Ｄ変換する計測系、３０５は該計測
系からの計測データとキーボードから入力された撮影条
件に基づいて画像再生に必要な各種演算を行う処理装
置、３０６は再生結果及び撮影条件を表示するＣＲＴ、
３０７は上記構成における各系の動作のコントロール手
順を格納しておくパルスシーケンスファイル、３０８は
キーボードから入力された撮影条件に従って上記パルス
シーケンスファイルからローディングしたパルスシーケ
ンスに基づいて、各装置の動作タイミングをコントロー
ルするシーケンス制御部である。
【００１８】本発明を実施する処理フローの実施例を図
１に、パルスシーケンスの一実施例を図２に示す。
１に、パルスシーケンスの一実施例を図２に示す。
【００１９】図２のパルスシーケンスはよく知られたＧ
ＦＥ（gradient field echo)シーケンスであるが、簡単
に説明しておく。まず、励起系３０２からα度（＜９０
度）の励起用パルス２０２を発生させ、水素スピンを励
起する。この時同時に、特定のスライスを選択するため
にスライス選択傾斜磁場２０２を傾斜磁場発生系３０３
にて発生させる（以下、全ての傾斜磁場は、３０３にて
発生させる）。これに引き続き、この傾斜磁場によりス
ライス方向に位相が乱れる影響をキャンセルするため、
反転傾斜磁場２０３を印加する。位相エンコード磁場２
０４は、スライス内の１軸方向（この例ではｙ方向）の
位置を分離するためのもので、シーケンスを繰り返す度
に強度を順次変化させる。読みだし磁場の反転磁場２０
５は、読みだし傾斜磁場２０６での計測信号のピークが
中央に来るようにシフトさせる傾斜磁場である。これに
引き続き、読みだし傾斜磁場２０６を印加し、この時同
時に計測系３０４を用いてＭＲＩ信号２０７を計測す
る。
ＦＥ（gradient field echo)シーケンスであるが、簡単
に説明しておく。まず、励起系３０２からα度（＜９０
度）の励起用パルス２０２を発生させ、水素スピンを励
起する。この時同時に、特定のスライスを選択するため
にスライス選択傾斜磁場２０２を傾斜磁場発生系３０３
にて発生させる（以下、全ての傾斜磁場は、３０３にて
発生させる）。これに引き続き、この傾斜磁場によりス
ライス方向に位相が乱れる影響をキャンセルするため、
反転傾斜磁場２０３を印加する。位相エンコード磁場２
０４は、スライス内の１軸方向（この例ではｙ方向）の
位置を分離するためのもので、シーケンスを繰り返す度
に強度を順次変化させる。読みだし磁場の反転磁場２０
５は、読みだし傾斜磁場２０６での計測信号のピークが
中央に来るようにシフトさせる傾斜磁場である。これに
引き続き、読みだし傾斜磁場２０６を印加し、この時同
時に計測系３０４を用いてＭＲＩ信号２０７を計測す
る。
【００２０】図１、図２を用いて本発明の実施方法を以
下に説明する。
下に説明する。
【００２１】〈ステップ１０１〉処理装置３０５に接続
されたキーボードから撮影シーケンス、撮影パラメー
タ、撮影位置などの撮影条件を入力し、それに基づきパ
ルスシーケンスファイル３０７から該当のパルスシーケ
ンスをローディングする。ここでは、そのシーケンス
が、図２のＧＦＥシーケンスとして説明する。
されたキーボードから撮影シーケンス、撮影パラメー
タ、撮影位置などの撮影条件を入力し、それに基づきパ
ルスシーケンスファイル３０７から該当のパルスシーケ
ンスをローディングする。ここでは、そのシーケンス
が、図２のＧＦＥシーケンスとして説明する。
【００２２】〈ステップ１０２〉ローディングした図２
のパルスシーケンスに基づき、励起系３０２、傾斜磁場
発生系３０３、計測系３０４をコントロールし、位相エ
ンコード数だけシーケンスを繰返してレファレンス画像
用のデータを計測し、それに基づき処理装置３０５でレ
ファレンス画像を再生する。通常、画像再生は信号計測
と並行して行う。計測時間は位相エンコード数と繰返し
時間ＴＲの積で決まる。例えば、１２８エンコードでＴ
Ｒが１０ｍｓｅｃならば、計測時間は約１秒になる。
のパルスシーケンスに基づき、励起系３０２、傾斜磁場
発生系３０３、計測系３０４をコントロールし、位相エ
ンコード数だけシーケンスを繰返してレファレンス画像
用のデータを計測し、それに基づき処理装置３０５でレ
ファレンス画像を再生する。通常、画像再生は信号計測
と並行して行う。計測時間は位相エンコード数と繰返し
時間ＴＲの積で決まる。例えば、１２８エンコードでＴ
Ｒが１０ｍｓｅｃならば、計測時間は約１秒になる。
【００２３】〈ステップ１０３〉レファレンス画像を用
いて位相補正用のデータを計算する。具体的には、複素
数であるレファレンス画像の絶対値を１に正規化した
後、複素共役演算（虚数部の符号反転）したデータを計
算する。
いて位相補正用のデータを計算する。具体的には、複素
数であるレファレンス画像の絶対値を１に正規化した
後、複素共役演算（虚数部の符号反転）したデータを計
算する。
【００２４】〈ステップ１０４〉機能計測を行うための
画像を〈ステップ１０２〉と同じく、図２のＧＦＥパル
スシーケンスに基づき計測し、再生して得る。計測処理
と再生処理も同じく、同時に行う。
画像を〈ステップ１０２〉と同じく、図２のＧＦＥパル
スシーケンスに基づき計測し、再生して得る。計測処理
と再生処理も同じく、同時に行う。
【００２５】〈ステップ１０５〉〈ステップ１０２〉の
補正用のデータに、〈ステップ１０４〉の再生画像を複
素乗算し、画像の位相補正を行う。これにより、変化の
ない部分の位相成分をゼロにし、変化のある部分のみの
位相成分を得ることができる。
補正用のデータに、〈ステップ１０４〉の再生画像を複
素乗算し、画像の位相補正を行う。これにより、変化の
ない部分の位相成分をゼロにし、変化のある部分のみの
位相成分を得ることができる。
【００２６】〈ステップ１０６〉得られた位相を濃度に
変換し、表示する。この時、位相成分のあらかじめ定め
た所定値を超えた部分のみを抽出し、再生画像の振幅値
の絶対値像に重ねあわせて表示しても良い。さらに重ね
る際、カラー合成することでより変化部分を強調して表
示する事が可能になる。
変換し、表示する。この時、位相成分のあらかじめ定め
た所定値を超えた部分のみを抽出し、再生画像の振幅値
の絶対値像に重ねあわせて表示しても良い。さらに重ね
る際、カラー合成することでより変化部分を強調して表
示する事が可能になる。
【００２７】なお、上記ステップにおいて、下記の代替
処理が考えられる。
処理が考えられる。
【００２８】〈ステップ１０２〉において、直ちに撮影
に入る前にカラ打ちをしばらく行い定常状態になってか
らレファレンス画像を撮影した方が、安定した画像が得
られる。
に入る前にカラ打ちをしばらく行い定常状態になってか
らレファレンス画像を撮影した方が、安定した画像が得
られる。
【００２９】〈ステップ１０２〉と〈ステップ１０４〉
の画像計測に用いるＧＦＥシーケンスは、図２のＴＲ短
縮シーケンスの代わりに、連続的なグラディエントエコ
ーを用いて計測信号を得る超高速シーケンス（例えば、
エコープレーナ）を用いても良い。この場合、１８０度
パルスを併用する場合は、９０度−１８０度間隔と１８
０度−エコー中心間隔の差を、前に述べたように４ｍｓ
ｅｃ程度ずらすことで対応させる。
の画像計測に用いるＧＦＥシーケンスは、図２のＴＲ短
縮シーケンスの代わりに、連続的なグラディエントエコ
ーを用いて計測信号を得る超高速シーケンス（例えば、
エコープレーナ）を用いても良い。この場合、１８０度
パルスを併用する場合は、９０度−１８０度間隔と１８
０度−エコー中心間隔の差を、前に述べたように４ｍｓ
ｅｃ程度ずらすことで対応させる。
【００３０】〈ステップ１０３〉において、レファレン
ス画像の位相を計算しておき、〈ステップ１０５〉で逆
位相をかけても良い。その場合は、フェーズアンラップ
処理が新たに必要になることに注意しなければならな
い。
ス画像の位相を計算しておき、〈ステップ１０５〉で逆
位相をかけても良い。その場合は、フェーズアンラップ
処理が新たに必要になることに注意しなければならな
い。
【００３１】〈ステップ１０５〉において、位相変化量
から、(２）式に基づき、γ、ｔが既知なので、ΔＢが
算出でき、ヘモグロビンと結合した酸素の量から、脳で
あれば機能の活性度を推定することも可能である。
から、(２）式に基づき、γ、ｔが既知なので、ΔＢが
算出でき、ヘモグロビンと結合した酸素の量から、脳で
あれば機能の活性度を推定することも可能である。
【００３２】上記例では、最初の画像をレファレンスと
して選び、その差を検出しているが、次に、医師や放射
線技師等のユーザが画像をＣＲＴ上で見ながら、適切な
部位あるいはタイミングを見計らい、指示した時点から
実際の機能観察を行う別の処理手順の例について説明す
る。
して選び、その差を検出しているが、次に、医師や放射
線技師等のユーザが画像をＣＲＴ上で見ながら、適切な
部位あるいはタイミングを見計らい、指示した時点から
実際の機能観察を行う別の処理手順の例について説明す
る。
【００３３】〈ステップ４０１〉処理装置３０５に接続
されたキーボードから撮影シーケンス、撮影パラメー
タ、撮影位置などの撮影条件を入力し、それに基づきパ
ルスシーケンスファイル３０７から該当のパルスシーケ
ンスをローディングする。この時、実際に画像を観察し
ながら指定したい場合は、撮影位置の自動移動モードを
指定しておく。撮影位置を先程と同様に、そのシーケン
スが、図２のＧＦＥシーケンスとして説明する。
されたキーボードから撮影シーケンス、撮影パラメー
タ、撮影位置などの撮影条件を入力し、それに基づきパ
ルスシーケンスファイル３０７から該当のパルスシーケ
ンスをローディングする。この時、実際に画像を観察し
ながら指定したい場合は、撮影位置の自動移動モードを
指定しておく。撮影位置を先程と同様に、そのシーケン
スが、図２のＧＦＥシーケンスとして説明する。
【００３４】〈ステップ４０２〉ローディングした図２
のパルスシーケンスに基づき、励起系３０２、傾斜磁場
発生系３０３、計測系３０４をコントロールし、位相エ
ンコード数だけシーケンスを繰返してレファレンス画像
用のデータを計測し、それに基づき処理装置３０５でレ
ファレンス画像を再生する。
のパルスシーケンスに基づき、励起系３０２、傾斜磁場
発生系３０３、計測系３０４をコントロールし、位相エ
ンコード数だけシーケンスを繰返してレファレンス画像
用のデータを計測し、それに基づき処理装置３０５でレ
ファレンス画像を再生する。
【００３５】〈ステップ４０３〉前ステップで得られた
画像（複素数）の絶対値と位相を計算する。
画像（複素数）の絶対値と位相を計算する。
【００３６】〈ステップ４０４〉絶対値及び、位相を濃
度に変換し、画像を例えば並列に表示する。その後、
〈ステップ４０２〉に戻り、処理を繰り返す。この時
〈ステップ４０１〉の撮影条件にて位置の自動移動モー
ドが設定してあれば、その時のパラメータに応じてスラ
イス位置を移動させて、〈ステップ４０２〉に戻る。移
動は励起パルスの変調周波数を変化させることで容易に
実施することができる。
度に変換し、画像を例えば並列に表示する。その後、
〈ステップ４０２〉に戻り、処理を繰り返す。この時
〈ステップ４０１〉の撮影条件にて位置の自動移動モー
ドが設定してあれば、その時のパラメータに応じてスラ
イス位置を移動させて、〈ステップ４０２〉に戻る。移
動は励起パルスの変調周波数を変化させることで容易に
実施することができる。
【００３７】〈ステップ４０５〉ユーザは、〈ステップ
４０４〉で表示された画像を見て適当な画像（タイミン
グ及び位置）と判断したときに、割込みをシステムにか
ける。
４０４〉で表示された画像を見て適当な画像（タイミン
グ及び位置）と判断したときに、割込みをシステムにか
ける。
【００３８】〈ステップ４０６〉割込みがかかった時点
に表示されているレファレンス画像を用いて位相補正用
のデータを計算し、メモリに保存する。具体的には、レ
ファレンス画像の絶対値を１に正規化した後、複素共役
演算（虚数部の符号反転）したデータを計算し、メモリ
に保存する。この時点で、スライス位置は固定する。
に表示されているレファレンス画像を用いて位相補正用
のデータを計算し、メモリに保存する。具体的には、レ
ファレンス画像の絶対値を１に正規化した後、複素共役
演算（虚数部の符号反転）したデータを計算し、メモリ
に保存する。この時点で、スライス位置は固定する。
【００３９】〈ステップ４０７〉機能計測を行うための
画像を〈ステップ４０２〉と同じく、図２のＧＦＥパル
スシーケンスに基づき計測し、再生して得る。
画像を〈ステップ４０２〉と同じく、図２のＧＦＥパル
スシーケンスに基づき計測し、再生して得る。
【００４０】〈ステップ４０８〉〈ステップ４０６〉の
保存してある補正用のデータに、〈ステップ４０７〉の
再生画像を複素乗算し、画像の位相補正を行う。これに
より、変化のない部分の位相成分をゼロにし、変化のあ
る部分のみの位相成分を得ることができる。
保存してある補正用のデータに、〈ステップ４０７〉の
再生画像を複素乗算し、画像の位相補正を行う。これに
より、変化のない部分の位相成分をゼロにし、変化のあ
る部分のみの位相成分を得ることができる。
【００４１】〈ステップ４０９〉得られた位相を濃度に
変換し、表示する。この時、位相成分のあらかじめ定め
た所定値を超えた部分のみを抽出し、再生画像の振幅値
の絶対値像に重ね合わせて表示しても良い。さらに重ね
る際、カラー合成することでより変化部分を強調して表
示する事が可能になる。
変換し、表示する。この時、位相成分のあらかじめ定め
た所定値を超えた部分のみを抽出し、再生画像の振幅値
の絶対値像に重ね合わせて表示しても良い。さらに重ね
る際、カラー合成することでより変化部分を強調して表
示する事が可能になる。
【００４２】上記実施例によれば、任意の時点の状態を
レファレンスとして、その時点からの変化を検出できる
ので、外部から光等で脳に刺激を与えた結果を容易に観
察することが可能になる。
レファレンスとして、その時点からの変化を検出できる
ので、外部から光等で脳に刺激を与えた結果を容易に観
察することが可能になる。
【００４３】また、呼吸や心臓のゲーテイングをかけて
撮影することで周期的な動きのある部位での撮影を実施
してすることも可能である。
撮影することで周期的な動きのある部位での撮影を実施
してすることも可能である。
【００４４】
【発明の効果】以上、詳細に説明したごとく、本発明に
よれば、「ディフュージョン強調パルス」及び、「造影
剤」を用いること無くディフュージョン／パフュージョ
ン領域を磁化率の変化として検出し、秒オーダで時間的
に変化する様子を映像化する機能イメージングを実現す
ることが可能になる。
よれば、「ディフュージョン強調パルス」及び、「造影
剤」を用いること無くディフュージョン／パフュージョ
ン領域を磁化率の変化として検出し、秒オーダで時間的
に変化する様子を映像化する機能イメージングを実現す
ることが可能になる。
【図１】本発明の第１の実施例の動作を示すフローチャ
ートである。
ートである。
【図２】本発明を実施するパルスシーケンスの一例であ
る。
る。
【図３】本発明を実施するブロック構成図の一例であ
る。
る。
【図４】本発明の第２の実施例の動作を示すフローチャ
ートである。
ートである。
３０１：静磁場発生系、３０２：励起系、３０３：傾斜
磁場発生系、３０４：計測系、３０５：処理装置、３０
６：ＣＲＴ、３０７：パルスシーケンスファイル、３０
８：シーケンス制御部。
磁場発生系、３０４：計測系、３０５：処理装置、３０
６：ＣＲＴ、３０７：パルスシーケンスファイル、３０
８：シーケンス制御部。
Claims (7)
- 【請求項１】静磁場、傾斜磁場、高周波磁場の発生装置
と検査対象物からの核磁気共鳴信号を取り出す検出装置
と、画像再生を含む各種演算を行う処理装置を有する磁
気共鳴イメージング装置において、連続的に画像を撮影
する手段と、その画像の位相を算出する手段と、特定の
１枚の画像の位相をレファレンスとして他の画像との位
相差を算出する手段と、その位相差を表示する手段を持
つことを特徴とする磁気共鳴イメージング装置。 - 【請求項２】連続的に画像を撮影する手段は、同じ部位
の画像を撮影する手段と連続的に位置を移動させる手段
とを含む特許請求項１の磁気共鳴イメージング装置。 - 【請求項３】連続的に画像を撮影する手段が、励起間隔
（ＴＲ）が短いＧＦＥ(gradient field echo)シーケン
スである特許請求項１の磁気共鳴イメージング装置。 - 【請求項４】画像の位相を算出する手段は、画像を撮影
する手段と同時に並行して行う特許請求項１の磁気共鳴
イメージング装置。 - 【請求項５】レファフレンスに用いる特定の１枚の画像
を、撮影中に外部から指定できる手段をもつ特許請求項
１の磁気共鳴イメージング装置。 - 【請求項６】表示する位相差を、時間的に変化する磁化
率の変化率に換算して表示する手段を持つ特許請求項１
の磁気共鳴イメージング装置。 - 【請求項７】表示する位相差を、通常の形態学的情報を
表わすＭＲＩ断層面上に重ねて表示する手段を持つ特許
請求項１の磁気共鳴イメージング装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4212588A JPH0654831A (ja) | 1992-08-10 | 1992-08-10 | 磁気共鳴機能イメージング装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4212588A JPH0654831A (ja) | 1992-08-10 | 1992-08-10 | 磁気共鳴機能イメージング装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0654831A true JPH0654831A (ja) | 1994-03-01 |
Family
ID=16625191
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP4212588A Pending JPH0654831A (ja) | 1992-08-10 | 1992-08-10 | 磁気共鳴機能イメージング装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0654831A (ja) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006524088A (ja) * | 2003-04-24 | 2006-10-26 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 繊維トラッキング磁気共鳴イメージング方法 |
| WO2010061736A1 (ja) | 2008-11-25 | 2010-06-03 | 日鉱金属株式会社 | 印刷回路用銅箔 |
| WO2011138876A1 (ja) | 2010-05-07 | 2011-11-10 | Jx日鉱日石金属株式会社 | 印刷回路用銅箔 |
| WO2012169249A1 (ja) | 2011-06-07 | 2012-12-13 | Jx日鉱日石金属株式会社 | 液晶ポリマー銅張積層板及び当該積層板に用いる銅箔 |
| KR20150119489A (ko) | 2011-03-30 | 2015-10-23 | 제이엑스 닛코 닛세키 킨조쿠 가부시키가이샤 | 인쇄 회로용 동박 |
| KR20170077280A (ko) | 2012-09-28 | 2017-07-05 | 제이엑스금속주식회사 | 캐리어 부착 동박 및 캐리어 부착 동박을 사용한 구리 피복 적층판 |
| CN106999111A (zh) * | 2014-10-01 | 2017-08-01 | 纽洛斯公司 | 用于检测不可见人类情感的系统和方法 |
| US11857323B2 (en) | 2017-10-24 | 2024-01-02 | Nuralogix Corporation | System and method for camera-based stress determination |
-
1992
- 1992-08-10 JP JP4212588A patent/JPH0654831A/ja active Pending
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006524088A (ja) * | 2003-04-24 | 2006-10-26 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 繊維トラッキング磁気共鳴イメージング方法 |
| WO2010061736A1 (ja) | 2008-11-25 | 2010-06-03 | 日鉱金属株式会社 | 印刷回路用銅箔 |
| WO2011138876A1 (ja) | 2010-05-07 | 2011-11-10 | Jx日鉱日石金属株式会社 | 印刷回路用銅箔 |
| KR20150119489A (ko) | 2011-03-30 | 2015-10-23 | 제이엑스 닛코 닛세키 킨조쿠 가부시키가이샤 | 인쇄 회로용 동박 |
| WO2012169249A1 (ja) | 2011-06-07 | 2012-12-13 | Jx日鉱日石金属株式会社 | 液晶ポリマー銅張積層板及び当該積層板に用いる銅箔 |
| US9060431B2 (en) | 2011-06-07 | 2015-06-16 | Jx Nippon Mining & Metals Corporation | Liquid crystal polymer copper-clad laminate and copper foil used for said laminate |
| KR20170077280A (ko) | 2012-09-28 | 2017-07-05 | 제이엑스금속주식회사 | 캐리어 부착 동박 및 캐리어 부착 동박을 사용한 구리 피복 적층판 |
| CN106999111A (zh) * | 2014-10-01 | 2017-08-01 | 纽洛斯公司 | 用于检测不可见人类情感的系统和方法 |
| US11857323B2 (en) | 2017-10-24 | 2024-01-02 | Nuralogix Corporation | System and method for camera-based stress determination |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9535147B2 (en) | Magnetic resonance imaging apparatus | |
| JP4444413B2 (ja) | ４重フィールドエコーシーケンスを用いて水と脂肪を定量的にｍｒ撮影する装置 | |
| EP1310211B1 (en) | Magnetic resonance imaging device and method | |
| US6195579B1 (en) | Contrast detection and guided reconstruction in contrast-enhanced magnetic resonance angiography | |
| JP5979327B2 (ja) | 磁気共鳴イメージング装置、その作動方法及び時系列画像作成プログラム | |
| JP3512482B2 (ja) | 磁気共鳴映像装置 | |
| JP3920140B2 (ja) | Ｍｒｉ装置及びフロー定量化装置 | |
| JP5942271B2 (ja) | 磁気共鳴イメージング装置及び流体強調画像取得法 | |
| EP0798566A1 (en) | Three-dimensional digital subtraction magnetic resonance angiography | |
| US20110064294A1 (en) | Magnetic resonance imaging apparatus and blood vessel image acquiring method | |
| WO2000065995A1 (en) | Magnetic resonance imaging method and device therefor | |
| US20100090696A1 (en) | Method, processor, and magnetic resonance apparatus for selective presentation of lung movement | |
| JP2008272248A (ja) | 磁気共鳴イメージング装置 | |
| JPH11322A (ja) | 磁気共鳴映像法（ｍｒｉ）における横方向緩和速度への可逆性寄与分測定方法 | |
| US6198960B1 (en) | Flip angle modulated magnetic resonance angiography | |
| US7715900B2 (en) | Quadruple inversion recovery for quantitative contrast-enhanced black blood imaging | |
| JPH0654831A (ja) | 磁気共鳴機能イメージング装置 | |
| JP3998814B2 (ja) | 磁気共鳴イメージング装置 | |
| JPS63186639A (ja) | 血流イメ−ジング方式 | |
| JP2021520951A (ja) | 磁化の反転状態の評価を伴う動脈スピンラベリング法 | |
| Hagio et al. | T2 mapping of the heart with a double-inversion radial fast spin-echo method with indirect echo compensation | |
| US20070286808A1 (en) | Method for display presentation of lymph nodes | |
| JP3907944B2 (ja) | 磁気共鳴イメージング方法及び装置 | |
| JPS61106140A (ja) | 高分解能核磁気共鳴イメ−ジング方式 | |
| JPH07108288B2 (ja) | Ｎｍｒイメ−ジング方法 |