JPS61196147A - 核磁気共鳴画像装置 - Google Patents
核磁気共鳴画像装置Info
- Publication number
- JPS61196147A JPS61196147A JP60038236A JP3823685A JPS61196147A JP S61196147 A JPS61196147 A JP S61196147A JP 60038236 A JP60038236 A JP 60038236A JP 3823685 A JP3823685 A JP 3823685A JP S61196147 A JPS61196147 A JP S61196147A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- image
- pulse
- gradient magnetic
- signal
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
- G06F3/011—Arrangements for interaction with the human body, e.g. for user immersion in virtual reality
- G06F3/014—Hand-worn input/output arrangements, e.g. data gloves
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Computer Interaction (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
（産業上の利用分野）
本発明は、核磁気共鳴画像装置（以下核磁気共鳴をＮＭ
Ｒと略す）に関し、特に緩和時間Ｔ＋。
Ｒと略す）に関し、特に緩和時間Ｔ＋。
Ｔ２およびプロトン密度ρの計算画像を求める手段の改
善に関するものである。
善に関するものである。
（従来の技術）
従来より、ＮＭＲ画像装置において、測定した画像から
医学上有用とされている縦緩和時間Ｔ１値に関する画１
１（Ｔｚ像）や横緩和時間Ｔ２値に関する画像（Ｔ２像
）を求める技法があった。
医学上有用とされている縦緩和時間Ｔ１値に関する画１
１（Ｔｚ像）や横緩和時間Ｔ２値に関する画像（Ｔ２像
）を求める技法があった。
（発明が解決しようとする問題点）
しかしながら、ＴＩ像とＴ２像は次のように別々の方法
により求められていた。
により求められていた。
■Ｔ１像については、例えば、次のようにして計算され
る。第７図に示すような反転回復法（Ｉｎｖｅｒｓｉｏ
ｎ Ｒｅｃｏｖｅｒｙ法：以下ＩＲ法と略す）とスピ
ンエコー法（３ｐｉｎ Ｅ ｃｈｏ法：以下ＳＥ法と
略す）とを併せて適用したＩ Ｒ２Ｈ法と、第７図に示
すような飽和回復法（Ｓ ａｔｕｒａｔｉｏｎ Ｒｅ
ｃｏｖｅｒｙ ：以下ＳＲ法と略す）とＳＥ法とを併せ
て適用した５Ｒ８Ｅ法により、各１枚ずつの原画像を得
、この２枚の画像と、信号強度の近似式を用いて計算す
る。
る。第７図に示すような反転回復法（Ｉｎｖｅｒｓｉｏ
ｎ Ｒｅｃｏｖｅｒｙ法：以下ＩＲ法と略す）とスピ
ンエコー法（３ｐｉｎ Ｅ ｃｈｏ法：以下ＳＥ法と
略す）とを併せて適用したＩ Ｒ２Ｈ法と、第７図に示
すような飽和回復法（Ｓ ａｔｕｒａｔｉｏｎ Ｒｅ
ｃｏｖｅｒｙ ：以下ＳＲ法と略す）とＳＥ法とを併せ
て適用した５Ｒ８Ｅ法により、各１枚ずつの原画像を得
、この２枚の画像と、信号強度の近似式を用いて計算す
る。
５Ｒ８Ｅ法は第８図に示すように９０”パルス印加の後
に１８０°パルスを印加してエコー信号を得るようにし
たパルスシーケンスで、９０°パルスからエコー信号の
中心までの時間をＴ、 、９Ｇ’パルス印加から次のビ
ューでの９０°パルス印加までの時開をＴｒとしている
。
に１８０°パルスを印加してエコー信号を得るようにし
たパルスシーケンスで、９０°パルスからエコー信号の
中心までの時間をＴ、 、９Ｇ’パルス印加から次のビ
ューでの９０°パルス印加までの時開をＴｒとしている
。
また、Ｉ Ｒ２Ｈ法は第７図に示すように第８図の５Ｒ
８Ｅ法の各９０°パルスの前にインバージョン・リカバ
リ用の１８００パルスを印加するようにしたパルスシー
ケンスで、インバージョン・リカバリ用の１８０°パル
スの印加から９０°パルスの印加までの時間を７ｄ、ｓ
ｏ°パルスからエコー信号の中心までの時間をＴｇ、イ
ンバージョン・リカバリ用の１８０°パルスの印加から
次のビューでの１８００パルス印加までの時間をＴｒと
している。
８Ｅ法の各９０°パルスの前にインバージョン・リカバ
リ用の１８００パルスを印加するようにしたパルスシー
ケンスで、インバージョン・リカバリ用の１８０°パル
スの印加から９０°パルスの印加までの時間を７ｄ、ｓ
ｏ°パルスからエコー信号の中心までの時間をＴｇ、イ
ンバージョン・リカバリ用の１８０°パルスの印加から
次のビューでの１８００パルス印加までの時間をＴｒと
している。
５Ｒ８Ｅ法での信号強度の理論式ＩｓｇはｈＲ＝ Ｌ・
ａｒｐ（−Ｔｓ／ＴＩ ）ｆ ｌ−２、ｅｘｐ（−Ｔｒ
／Ｔ ＋Ｔｓ ／２Ｔ＋ ）寸ｅｘｐ（−Ｔｒ／工）ｊ また、ＩＲ８Ｅ法での信号強度の理論式ｒｕはＩｍ−Ｉ
ｏ・ｅＸＰ（−Ｔ−／Ｔｚ）（１−２ａＸＰ（−７１／
Ｔ； ）＋２−ｅｇｐ（−ＴｒｌＴｒ ｆＴ、１２石）
−ｅｘｐ（−Ｔ／Ｔｉ）７この理論式に対し、ここで、
７ ｒ＞ Ｔ、としてｅＸｐ （−ＴｒｌＴｒ ）＝０
とすれば、ＩＳＲ寺Ｌ−ｆＪｐ（−Ｔｓ／Ｔ２ ）ｈａ
／ｌ５ｙｔ＝ ｌ−２’ｅＸｒ（−Ｔａ／Ｔｉ）ＴＩ
＝ Ｔ、Ｉ／ｊＡ（２シ／ （Ｉｓｔ −Ｉ工、））こ
の式からＴＩ値を求める。
ａｒｐ（−Ｔｓ／ＴＩ ）ｆ ｌ−２、ｅｘｐ（−Ｔｒ
／Ｔ ＋Ｔｓ ／２Ｔ＋ ）寸ｅｘｐ（−Ｔｒ／工）ｊ また、ＩＲ８Ｅ法での信号強度の理論式ｒｕはＩｍ−Ｉ
ｏ・ｅＸＰ（−Ｔ−／Ｔｚ）（１−２ａＸＰ（−７１／
Ｔ； ）＋２−ｅｇｐ（−ＴｒｌＴｒ ｆＴ、１２石）
−ｅｘｐ（−Ｔ／Ｔｉ）７この理論式に対し、ここで、
７ ｒ＞ Ｔ、としてｅＸｐ （−ＴｒｌＴｒ ）＝０
とすれば、ＩＳＲ寺Ｌ−ｆＪｐ（−Ｔｓ／Ｔ２ ）ｈａ
／ｌ５ｙｔ＝ ｌ−２’ｅＸｒ（−Ｔａ／Ｔｉ）ＴＩ
＝ Ｔ、Ｉ／ｊＡ（２シ／ （Ｉｓｔ −Ｉ工、））こ
の式からＴＩ値を求める。
■Ｔ２像を求める場合は、例えば、刊行物「映像情報（
Ｍ）Ｊ１９８４年６月号（Ｖｏ ｌ、 １６ Ｎｏ
。
Ｍ）Ｊ１９８４年６月号（Ｖｏ ｌ、 １６ Ｎｏ
。
１１）の第５７０頁ないし第５７６頁に記載されたｃＰ
ＭＧ法により複数個のエコーデータからＴｌ 、ρを消
去して最小２乗法によりＴ２値を求めるようにしている
。
ＭＧ法により複数個のエコーデータからＴｌ 、ρを消
去して最小２乗法によりＴ２値を求めるようにしている
。
１回のデータ収集で複数個のエコーデータを連続的に取
り出しＴ２値を求められるようにしたＣＰ法には、印加
するパルスの長さが不完全であればその誤差がエコーを
得るに従−い累積され、結果としてＴ２値に誤差を生ず
ると言う欠点があったが、ＣＰＭＧ法はこれを解決した
もので、第９図に示すように９０”パルスの後に１８ｏ
°パルスを０回繰返し印加してｎ個のエコーを発生させ
るようにしたパルスシーケンスである。
り出しＴ２値を求められるようにしたＣＰ法には、印加
するパルスの長さが不完全であればその誤差がエコーを
得るに従−い累積され、結果としてＴ２値に誤差を生ず
ると言う欠点があったが、ＣＰＭＧ法はこれを解決した
もので、第９図に示すように９０”パルスの後に１８ｏ
°パルスを０回繰返し印加してｎ個のエコーを発生させ
るようにしたパルスシーケンスである。
このような手法による従来の方法においては次のような
欠点があった。
欠点があった。
（１）ＴＩ像とＴ２像がそれぞれ別個に求められており
、Ｔ Ｉ ｓ Ｔ２ 、ρの計算画像が同時に得られな
い。
、Ｔ Ｉ ｓ Ｔ２ 、ρの計算画像が同時に得られな
い。
（２）近似式を用いているため正確な値が求まらない。
（３） Ｔ ？″＞ Ｔ、の条件のためＴｒを長くしな
ければならず１．全スキャンタイムが長い。
ければならず１．全スキャンタイムが長い。
・本発明の目的は、この様な点に鑑み、ＮＭＲ画像装置
において、３枚以上の原画像からＴＩ 、 Ｔ２、ρ計
算画像を正確かつ同時に求めるようなＮＭＲ画像装置を
提供することにある。
において、３枚以上の原画像からＴＩ 、 Ｔ２、ρ計
算画像を正確かつ同時に求めるようなＮＭＲ画像装置を
提供することにある。
本発明の他の目的は、任意のパルスシーケンスから得た
原画像を使用しても同様にＴ ＋ ＊ Ｔ２ ＋ρ計算
画像を正確かつ同時に求め得るようにしたＮＭＲ画像装
置を提供することにある。
原画像を使用しても同様にＴ ＋ ＊ Ｔ２ ＋ρ計算
画像を正確かつ同時に求め得るようにしたＮＭＲ画像装
置を提供することにある。
（問題点を解決するための手段）
この様な目的を達成するために本発明では、核磁気共鳴
画像装置において、 緩和時間およびプロトン密度に基づく計算画像を得るた
めの下記（イ）ないし（ロ）機能を有する制御手段を具
備したことを特徴とする。
画像装置において、 緩和時間およびプロトン密度に基づく計算画像を得るた
めの下記（イ）ないし（ロ）機能を有する制御手段を具
備したことを特徴とする。
記
（イ）決定されたパルスシーケンスを実行し、パルスパ
ラメータの異なる少なくとも３枚の原画像を得る。
ラメータの異なる少なくとも３枚の原画像を得る。
（ロ）前記少なくとも３枚の原画像と、各画像のスキャ
ン条件に応じた各画像の信号強度の理論式を用い、緩和
時間ＴＩ、Ｔ２およびプロトン密度ρの値を同時に求め
、これらの値をそれぞれ画像化して各計算画像を得る。
ン条件に応じた各画像の信号強度の理論式を用い、緩和
時間ＴＩ、Ｔ２およびプロトン密度ρの値を同時に求め
、これらの値をそれぞれ画像化して各計算画像を得る。
（実施例）
以下図面を用いて本発明の詳細な説明する。第１図は本
発明の方法を実施するためのＮＭＲ画像装置の一実施例
を示す要部構成図である。図において、１はマグネット
アセンブリで、内部には対象物を挿入するための空間部
分（孔）が設けられ、この空間部分を取巻くようにして
、対象物に一様静磁場Ｈｏを印加する主磁場コイル２と
、勾配磁場を発生するための勾配磁場コイル３（個別に
勾配磁場を発生することができるように構成されたＸ勾
配磁場コイル、ｙ勾配磁場コイル、−２勾配磁場コイル
より構成される）と、対象物内の原子核のスピンを励起
するための高周波パルスを与えるＲＦ送信コイル４と、
対象物からのＮＭＲ信号を検出する受信用コイル５等が
配置されている。
発明の方法を実施するためのＮＭＲ画像装置の一実施例
を示す要部構成図である。図において、１はマグネット
アセンブリで、内部には対象物を挿入するための空間部
分（孔）が設けられ、この空間部分を取巻くようにして
、対象物に一様静磁場Ｈｏを印加する主磁場コイル２と
、勾配磁場を発生するための勾配磁場コイル３（個別に
勾配磁場を発生することができるように構成されたＸ勾
配磁場コイル、ｙ勾配磁場コイル、−２勾配磁場コイル
より構成される）と、対象物内の原子核のスピンを励起
するための高周波パルスを与えるＲＦ送信コイル４と、
対象物からのＮＭＲ信号を検出する受信用コイル５等が
配置されている。
主磁場コイルは静磁場制御回路１５に、ＧＸ。
Ｇｙ、Ｇｚ各勾配磁場コイルは勾配磁場制御回路１４に
、ＲＦ送信コイルは電力増幅器１８に、そしてＮＭＲ信
号の受信用コイルはプリアンプ１９に、それぞれ接続さ
れている。
、ＲＦ送信コイルは電力増幅器１８に、そしてＮＭＲ信
号の受信用コイルはプリアンプ１９に、それぞれ接続さ
れている。
１３はコントローラで、勾配磁場や高周波磁場の発生シ
ーケンスを制御すると共に得られたＮＭＲ信号を波形メ
モリ２１に取込むために必要な制御を行う。
ーケンスを制御すると共に得られたＮＭＲ信号を波形メ
モリ２１に取込むために必要な制御を行う。
１７はゲート変調回路、１６は高周波信号を発生する高
周波発振器である。ゲート変調回路１７は、コントロー
ラ１３からの制御信号により高周波発振器１６が出力し
た高周波信号を適宜に変調し、所定の位相の高周波パル
スを生成する。この高周波パルスはＲＦ電力増幅器１８
を通してＲＦ送信コイル４に加えられる。
周波発振器である。ゲート変調回路１７は、コントロー
ラ１３からの制御信号により高周波発振器１６が出力し
た高周波信号を適宜に変調し、所定の位相の高周波パル
スを生成する。この高周波パルスはＲＦ電力増幅器１８
を通してＲＦ送信コイル４に加えられる。
１９は検出コイル５から得られるＮＭＲ信号を増幅する
プリアンプ、２０は高周波発振器の出力信号を参照して
ＮＭＲ信号を位相検波する位相検波回路、２１は位相検
波されたプリアンプからの波形信号を記憶する波形メモ
らで、ここにはＡ／Ｄ変換器を含んでいる。
プリアンプ、２０は高周波発振器の出力信号を参照して
ＮＭＲ信号を位相検波する位相検波回路、２１は位相検
波されたプリアンプからの波形信号を記憶する波形メモ
らで、ここにはＡ／Ｄ変換器を含んでいる。
１１は波形メモリ２１からの信号を受け、所定の信号処
理を施して断層像を得るコンピュータ、１２は得られた
断層像を表示するテレビジョンモニタのような表示器で
ある。
理を施して断層像を得るコンピュータ、１２は得られた
断層像を表示するテレビジョンモニタのような表示器で
ある。
この様な構成における計算画像作成の手順について次に
説明する。
説明する。
■先ず、所望のパルスシーケンスおよびスキャンパラメ
ータを決定する。ここでは原画像として５Ｒ８Ｅ法によ
り得られる３枚の画像を使用する場合を例にとる。３画
像のスキャンパラメータ等は、計算画像の分散又は標準
偏差を最小にするように選ぶ。
ータを決定する。ここでは原画像として５Ｒ８Ｅ法によ
り得られる３枚の画像を使用する場合を例にとる。３画
像のスキャンパラメータ等は、計算画像の分散又は標準
偏差を最小にするように選ぶ。
■前記■で決定されたパルス・シーケンスおよびスキャ
ンパラメータ等のこれらの条件はコントローラ１３に設
定される。コントローラ１３の制御に基づきゲート変調
回路１７をとおして第２図（イ）に示すような９０°パ
ルスを得、電力増幅器１８を介してＲＦ送信コイル４に
与え、対象物を励起する。この時同時に勾配磁場Ｇｚも
印加して（同図（ロ））、特定のスライス面内にあるス
ピンのみを選択励起する。
ンパラメータ等のこれらの条件はコントローラ１３に設
定される。コントローラ１３の制御に基づきゲート変調
回路１７をとおして第２図（イ）に示すような９０°パ
ルスを得、電力増幅器１８を介してＲＦ送信コイル４に
与え、対象物を励起する。この時同時に勾配磁場Ｇｚも
印加して（同図（ロ））、特定のスライス面内にあるス
ピンのみを選択励起する。
次に、勾配磁場Ｇｙにより位相エンコードを行い、それ
と同時に勾配磁場Ｇｘを印加して（同図（ニ））、エコ
ーを観測する準備をしておく。
と同時に勾配磁場Ｇｘを印加して（同図（ニ））、エコ
ーを観測する準備をしておく。
続いて、勾配磁場の印加を停止し、１８０°パルスを印
加しスピンを反覧させる。その後同図（ニ）に示すよう
にＧｘを印加しながら発生するエコー信号（同図（ホ）
）を受信コイル５で検出し、観測する。受信コイルで検
出されたスピンエコー信号は、プリアンプ１９１位相検
波回路２０を経て波形メモリ２１に蓄えられる。
加しスピンを反覧させる。その後同図（ニ）に示すよう
にＧｘを印加しながら発生するエコー信号（同図（ホ）
）を受信コイル５で検出し、観測する。受信コイルで検
出されたスピンエコー信号は、プリアンプ１９１位相検
波回路２０を経て波形メモリ２１に蓄えられる。
■上記■の方法で位相エンコード量の異なる多数のビュ
ーにおけるエコーを観測し、コンピュータ１１により原
画像を得る。
ーにおけるエコーを観測し、コンピュータ１１により原
画像を得る。
■上記■により得た３枚の原画像と理論式から反復補正
最小２乗法によりＴｌ ＊ Ｔ２ ｔρの計算画像を得
る。
最小２乗法によりＴｌ ＊ Ｔ２ ｔρの計算画像を得
る。
ここで、反復補正最小２乗法による計算例について説明
する。
する。
３枚の原画像のスキャンパラメータＴｒ、Ｔｇに関し、
第１の画像については、Ｔｒｌ −ＴＳＩＳ２Ｏ２像に
ついては、Ｔｒｚ ｌ ＴＳ２第３の画像については、
Ｔｙ３ 、 Ｔ５３とし、また５Ｒ８Ｅ法における信号
強度の理論式％式％ 次のように表わされる。
ついては、Ｔｒｚ ｌ ＴＳ２第３の画像については、
Ｔｙ３ 、 Ｔ５３とし、また５Ｒ８Ｅ法における信号
強度の理論式％式％ 次のように表わされる。
ここでＮ Ｃ３ｔＳＥ はスライスの影響を表わす係数
で、ＴＩ／ＴＴの関数である。例えばガウシアン９０°
パルスを用いればＣ５Ｒ５Ｉ”は次式で表わされ得る。
で、ＴＩ／ＴＴの関数である。例えばガウシアン９０°
パルスを用いればＣ５Ｒ５Ｉ”は次式で表わされ得る。
０．２＜ Ｔ／Ｔｒ＜ ＩＯ，ＯＴ”Ｃ報−８，Ｆ；
３７ Ｅ−６（王／１）５−２．’？５０８６ Ｆ−４
（工／Ｔｒ）”＋４．２７Ａ’７５ Ｅ３ （Ｔ；／Ｔ
ｒ−一−３，Ｉｒ７１？０２ Ｅ−２（Ｔ；／Ｔｒ）３
十＋、ｚｑ２６２１ニー１ （工／Ｔ）２−２．８’５
５４ Ｅ−１（工／Ｔ、ン寸１．０５５７ このスライスの影響とは、スライス中央以外では９０°
パルスにより磁化が９０°倒れないことによる影響であ
る。
３７ Ｅ−６（王／１）５−２．’？５０８６ Ｆ−４
（工／Ｔｒ）”＋４．２７Ａ’７５ Ｅ３ （Ｔ；／Ｔ
ｒ−一−３，Ｉｒ７１？０２ Ｅ−２（Ｔ；／Ｔｒ）３
十＋、ｚｑ２６２１ニー１ （工／Ｔ）２−２．８’５
５４ Ｅ−１（工／Ｔ、ン寸１．０５５７ このスライスの影響とは、スライス中央以外では９０°
パルスにより磁化が９０°倒れないことによる影響であ
る。
次に、第３図に示す計算処理のフローを参照して計算の
手順を述べる。初期近似値ｐ（ＴＩ、丁２、ρ）のまわ
りで理論式をテーラ展開し、その−次までをとれば、各
画像における信号の強度１１＋１２＋Ｉコは次のように
なる。
手順を述べる。初期近似値ｐ（ＴＩ、丁２、ρ）のまわ
りで理論式をテーラ展開し、その−次までをとれば、各
画像における信号の強度１１＋１２＋Ｉコは次のように
なる。
±（’２１Ｊ”ｎ下ｎＰ・４％ 十（此）％ＩＴｓｓＰ
’乙！ｙ 続いて、３つの画像の値と上式から最小２乗法によりΔ
Ｔ＋＋ ΔＴｅｌ Δρを求める。
’乙！ｙ 続いて、３つの画像の値と上式から最小２乗法によりΔ
Ｔ＋＋ ΔＴｅｌ Δρを求める。
ヘ ヘ へ
次に、Ｔ １ ＊ Ｔ２ ＊ρについての新しい近似値
Ｔ’＋ ＋ΔＴ１ Ｔ２ ＋Δ丁２ ρ ＋ Δ ρ とし、ＴＩ ＊ Ｔ２ ＋ρが収束するまで反復する。
Ｔ’＋ ＋ΔＴ１ Ｔ２ ＋Δ丁２ ρ ＋ Δ ρ とし、ＴＩ ＊ Ｔ２ ＋ρが収束するまで反復する。
収束不能の場合、初期値を変えて初めから計算を繰返す
。収束したと判断されたときは、その時の＾ハへ Ｔｌｏ丁２．ρが求めるＴＩ ａ Ｔ２ ＋ρである。
。収束したと判断されたときは、その時の＾ハへ Ｔｌｏ丁２．ρが求めるＴＩ ａ Ｔ２ ＋ρである。
このスライスの影響とは、スライス中央以外では９０°
パルスにより磁化が９゛０°倒れないことによる影響で
ある。
パルスにより磁化が９゛０°倒れないことによる影響で
ある。
以上の手順により、Ｔ ＋ ＋ Ｔ２ ＊ρの計算画像
を同時に得ることができる。
を同時に得ることができる。
なお、本発明では原画像の枚数は３枚に限ることはなく
、それ以上であってもよい。また、異なるパルスシーケ
ンスにおける混合３画像であってもよく、パルスシーケ
ンスの組合せは自由である。
、それ以上であってもよい。また、異なるパルスシーケ
ンスにおける混合３画像であってもよく、パルスシーケ
ンスの組合せは自由である。
各パルスシーケンスとそのときの信号強度の理論式は次
のとおりである。
のとおりである。
第４図ないし第６図は他のパルスシーケンスを示す図で
ある。ただし、勾配磁場の印加の様子は省略しである。
ある。ただし、勾配磁場の印加の様子は省略しである。
第４図は、５Ｒ８Ｅ法において各ビューごとに１８０°
パルスを２回ずつ加えるようにした５Ｒ２ＳＥ法の場合
であり、この場合の理論式は、Ｍｏを下式として、 士２・ｅｘｐ（−’ｒｌ下十Ｔｓ＋７２ ＴＩ） −ｅ
にＰ（−１’、−／Ｔ、 ） ｌ・／’ここでＣｇδε
はスライスの影響を表わす係数で、ＴｒｌＴｒの関数
である。
パルスを２回ずつ加えるようにした５Ｒ２ＳＥ法の場合
であり、この場合の理論式は、Ｍｏを下式として、 士２・ｅｘｐ（−’ｒｌ下十Ｔｓ＋７２ ＴＩ） −ｅ
にＰ（−１’、−／Ｔ、 ） ｌ・／’ここでＣｇδε
はスライスの影響を表わす係数で、ＴｒｌＴｒの関数
である。
第１のエコー信号の強度は
Ｍ ｏ −ｅＸｌ） （−Ｔ５ｔ ／ Ｔ２ ）第
２のエコー信号の強度は Ｍｏ−ｅＸｐ（ＴＳＩ ／ＴＺ Ｔ５２／Ｔ２ ）
で表わされる。
２のエコー信号の強度は Ｍｏ−ｅＸｐ（ＴＳＩ ／ＴＺ Ｔ５２／Ｔ２ ）
で表わされる。
ただし、Ｔ５１ は９０″パルス印加から第１のエコー
信号のピーク時までの時間、ＴＳ２は第１のエコー信号
のピーク時から第２のエコー信号のピーク時までの時間
である。
信号のピーク時までの時間、ＴＳ２は第１のエコー信号
のピーク時から第２のエコー信号のピーク時までの時間
である。
また、第５図はエコー信号発生後に１８０＠パルスおよ
び９０°パルスを印加して次のビュー開始までの待時間
を短縮するようにしたパルスシーケンスのもので（ＦＲ
８Ｅ法）、この場合の信号強度の理論式は である。
び９０°パルスを印加して次のビュー開始までの待時間
を短縮するようにしたパルスシーケンスのもので（ＦＲ
８Ｅ法）、この場合の信号強度の理論式は である。
また、第６図は各ビューごとに２回エコー信号を発生す
るようにした方式（ＦＲ２ＳＥ法）で、第５図のパルス
シーケンスに更に１８０°パルスを゛付加したものであ
る。
るようにした方式（ＦＲ２ＳＥ法）で、第５図のパルス
シーケンスに更に１８０°パルスを゛付加したものであ
る。
この場合の信号強度の理論式は、Ｍｏを下式として
第１のエコー信号の強度は
Ｍｏ −ｅＸＥｌ （Ｔｓ＋ ／Ｔ２ ）第２のエコ
ー信号の強度は Ｍｏ −ｅＸＤ （−Ｔｓ＋ ／Ｔ２−Ｔｓ２／Ｔ２
）で表わされる。
ー信号の強度は Ｍｏ −ｅＸＤ （−Ｔｓ＋ ／Ｔ２−Ｔｓ２／Ｔ２
）で表わされる。
（発明の効果）
以上説明したように、本発明によれば、次のような効果
がある。
がある。
■近似式を用いていないため正確な値を求めることがで
きる。
きる。
■理論式同士の演算により、Ｔ ｌ ｅ Ｔ ２ ｅρ
を消去することなく求められるので、Ｔｌ ＊ Ｔ２
＊ρの計算画像が同時に得られる。
を消去することなく求められるので、Ｔｌ ＊ Ｔ２
＊ρの計算画像が同時に得られる。
■理論式の形の特徴を使っていないので、パルスシーケ
ンス、スキャンパラメータは自由に選ぶことができる。
ンス、スキャンパラメータは自由に選ぶことができる。
従って、全スキャンタイムを短くすることができる。
■用いる画像数を多くすれば、広いＴＩ、Ｔ２の範囲で
分散を小さくすることができる。
分散を小さくすることができる。
第１図は本発明に係るＮＭＲ画像装置の一実施例を示す
要部構成図、第２図はパルスシーケンスの一実施例を示
す図、第３図は計算処理のフローを示す図、第４図ない
し第９図は各種のパルスシーケンスを示す図である。 １・・・マグネットアセンブリ、２・・・主磁場コイル
、３・・・勾配磁場コイル、４・・・ＲＦ送信コイル、
５・・・受信用コイル、１１・・・コンピュータ、１２
−・・表示器、１３・・・コントローラ、１４・・・勾
配磁場制御回路、１５・・・静磁場ＩＪ ｗＪ口路、１
６・・・高周波発振器、１７・・・ゲート変ＷＩＩ＠路
、１８・・・電力層ｍ！１，１９・・・プリアンプ、２
０・・・位相検波回路、２１・・・波形メモリ。 第３図
要部構成図、第２図はパルスシーケンスの一実施例を示
す図、第３図は計算処理のフローを示す図、第４図ない
し第９図は各種のパルスシーケンスを示す図である。 １・・・マグネットアセンブリ、２・・・主磁場コイル
、３・・・勾配磁場コイル、４・・・ＲＦ送信コイル、
５・・・受信用コイル、１１・・・コンピュータ、１２
−・・表示器、１３・・・コントローラ、１４・・・勾
配磁場制御回路、１５・・・静磁場ＩＪ ｗＪ口路、１
６・・・高周波発振器、１７・・・ゲート変ＷＩＩ＠路
、１８・・・電力層ｍ！１，１９・・・プリアンプ、２
０・・・位相検波回路、２１・・・波形メモリ。 第３図
Claims (1)
- 【特許請求の範囲】 １）対象物に高周波パルスおよび勾配磁場を印加して核
磁気共鳴信号を発生させ、この信号を検出し、検出した
信号を使って対象物の組織に関する画像を得るようにし
た核磁気共鳴画像装置において、 緩和時間およびプロトン密度に基づく計算画像を得るた
めの下記（イ）ないし（ロ）機能を有する制御手段を具
備したことを特徴とする核磁気共鳴画像装置。 記 （イ）決定されたパルスシーケンスを実行し、パルスパ
ラメータの異なる少なくとも３枚の原画像を得る。 （ロ）前記少なくとも３枚の原画像と、各画像のスキャ
ン条件に応じた各画像の信号強度の理論式を用い、緩和
時間Ｔ＿１、Ｔ＿２およびプロトン密度ρの値を同時に
求め、これらの値をそれぞれ画像化して各計算画像を得
る。 ２）前記信号強度の理論式は、画像に対してスライスに
よる影響のある場合にはスライスの影響を含んだ信号強
度の理論式を用いて計算画像を求めるようにしたことを
特徴とする特許請求の範囲第１項記載の核磁気共鳴画像
装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60038236A JPS61196147A (ja) | 1985-02-27 | 1985-02-27 | 核磁気共鳴画像装置 |
| US06/833,224 US4694250A (en) | 1985-02-27 | 1986-02-24 | Nuclear magnetic resonance imaging device |
| GB08604693A GB2173598B (en) | 1985-02-27 | 1986-02-26 | Nuclear magnetic resonance imaging device |
| DE19863606220 DE3606220A1 (de) | 1985-02-27 | 1986-02-26 | Nmr-abbildungsvorrichtung |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60038236A JPS61196147A (ja) | 1985-02-27 | 1985-02-27 | 核磁気共鳴画像装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS61196147A true JPS61196147A (ja) | 1986-08-30 |
| JPH0318453B2 JPH0318453B2 (ja) | 1991-03-12 |
Family
ID=12519662
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP60038236A Granted JPS61196147A (ja) | 1985-02-27 | 1985-02-27 | 核磁気共鳴画像装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS61196147A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104412118B (zh) * | 2012-07-02 | 2017-03-29 | 合成Mr公司 | 改进的磁共振采集的方法和系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6024463A (ja) * | 1983-07-20 | 1985-02-07 | Toshiba Corp | 核磁気共鳴映像法 |
-
1985
- 1985-02-27 JP JP60038236A patent/JPS61196147A/ja active Granted
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6024463A (ja) * | 1983-07-20 | 1985-02-07 | Toshiba Corp | 核磁気共鳴映像法 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104412118B (zh) * | 2012-07-02 | 2017-03-29 | 合成Mr公司 | 改进的磁共振采集的方法和系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0318453B2 (ja) | 1991-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5530259B2 (ja) | 磁気共鳴ｒｆ磁場計測のためのシステム、方法及び装置 | |
| US9316707B2 (en) | System and method of receive sensitivity correction in MR imaging | |
| JPH07163544A (ja) | Ｍｒｉシステムでのプレスキャン | |
| JPH0448451B2 (ja) | ||
| JPWO2004004563A1 (ja) | 磁気共鳴イメージング装置及び渦電流補償導出方法 | |
| US6288545B1 (en) | Method and apparatus for calibration of RF and gradient field time delays | |
| JP2945048B2 (ja) | 磁気共鳴イメージング装置 | |
| JPH0365971B2 (ja) | ||
| JPH0454448B2 (ja) | ||
| JPH0556140B2 (ja) | ||
| JPS61196147A (ja) | 核磁気共鳴画像装置 | |
| JPS61202286A (ja) | 核磁気共鳴画像装置 | |
| JPS6264345A (ja) | 核磁気共鳴撮像装置 | |
| JP3167038B2 (ja) | 磁気共鳴イメージング装置 | |
| JPS6267433A (ja) | Ｎｍｒイメ−ジング装置 | |
| JPH11290289A (ja) | 核磁気共鳴を用いた検査装置 | |
| JPS61266945A (ja) | 核磁気共鳴撮像装置 | |
| JPH11216124A (ja) | 核磁気共鳴検査装置 | |
| JPH0467848A (ja) | 磁気共鳴イメージング装置及びその静磁場強度測定表示方法 | |
| JPS60146140A (ja) | 核磁気共鳴による検査方法及びその装置 | |
| JPS61196146A (ja) | 核磁気共鳴画像装置 | |
| JP2928760B2 (ja) | 核磁気共鳴を用いた検査装置 | |
| JPH0421491B2 (ja) | ||
| JPS61180130A (ja) | Ｎｍｒ検査装置の静磁場強度分布測定方法 | |
| JP2004229866A (ja) | 磁気共鳴イメージング装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EXPY | Cancellation because of completion of term |