JPS6264345A - 核磁気共鳴撮像装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 （産業上の利用分野） 本発明は、核磁気共鳴撮像装置（以下核磁気共鳴をＮＭ
Ｒと略す）に関し、特に緩和時間Ｔ＋＋Ｔ２およびプロ
トン密度ρの計算画像を求める手段の改善に関するもの
である。

（従来の技術） 従来より、ＮＭＲｆｉ像装置において、測定した画像か
ら医学上有用とされている！緩和時間Ｔｔ値に関する画
＠（Ｔ＋鍮）や横緩和時間Ｔ２値に関する画像（Ｔ２像
）を求める技法があった。

（発明が解決しようとする問題点） しかしながら、このＴ１像とＴ２像は次のように別々の
方法により求められていた。

■Ｔ＋（ｉについては、例えば、次のようにして計算さ
れる。第５図に示すような反転回復法（ｌ　ｎｖｅｒｓ
ｉｏｎ　　Ｒｅｃｏｖｃｒｙ法：以下ＩＲ法と略す）と
スピンエコー法（Ｓｐｉｎ　　Ｅｃｈｏ法二以下ＳＥ法
と略す）とを併せて適用したＩＲ８Ｅ法と、第６図に示
すような飽和回復法（Ｓ　ａｔｕｒａｔｉｏｎ　　Ｒｅ
ｃｏｖｃｒｙ　：以下ＳＲ法と略す）とＳＥ法とを併せ
て適用した５Ｒ８Ｅ法により、各１枚ずつの原画像を得
、この２枚の画像と、信号強度の近似式を用いて計算す
る。

５Ｒ８Ｅ法は第６図に示すように９０’パルス印加の後
に１８０°パルスを印加してエコー信号を得るようにし
たパルスシーケンスで、　９０’パルスからエコー信号
の中心までの時間をＴｓ、９０”パルス印加から次のビ
ューでの９０′パルス印加までの時間をＴｒとしている
。

また、ＩＲ８Ｅ法は第５図に示すように第６図の５Ｒ８
Ｅ法の各９０°パルスの前にインバージョン・リカバリ
用の１８０８パルスを印加するようにしたパルスシーケ
ンスで、インバージョン・リカバリ用の１８０°パルス
の印加から９０°パルスの印加までの時間をＴａ、９０
’パルスからエコー信号の中心までの時間をＴｓＮイン
パージ３ン・リカバリ用の１８０′パルスの印加から次
のビューでの１８０６パルスの印加までの時間をＴｒと
している。

寸−ｅＸｐ（−Ｔヒ／丁、）） である。

このＩ！Ｉｌ！論式に対し、ここで、Ｔｒ＞Ｔ＋として
ｅｘｐ　　（−Ｔｒ／Ｔ＋　）−０とすれば、工銘十工
。・４）ｔＰ（−τｓ／Ｔλ）ＩＨ：］−０・−ｅｔＰ
（−Ｔｓ／　ＴＬ）ｌ　ｌ−２・ｅ−Ｘｉ’（−Ｔ、７
．／Ｔ＋）１ゆえに、 エエＦＬ／Ｉ吠＝１−１−θＰ（−７４／ＴＯＴ、　＝
　Ｔシム１２工鯨／（丁５に一エエＲ）１この式からＴ
１値を求める。

■Ｔ２像を求める場合は、例えば、刊行物「映像情報（
Ｍ）４１９８４年６月号（ＶＯｌ、１６　　Ｎｏ。

１１）の第５７０頁ないし第５７６頁に記載されたＣＰ
ＭＧ法により複数個のエコーデータからＴＩ＋ρを消去
して最小２乗法によりＴ２値を求めるようにしている。

なお、１回のデータ収集で複数個のエコーデータを連続
的に取り出しＴ２（ａを求め得るようにしたＣＰ法では
、印加するパルスの長さが不完全であればその誤差がエ
コーを得るに従い累積され、結果としてＴ２ＩＩｔＩに
誤差を生ずると言う欠点がありたが、ＣＰＭＧ法と呼ば
れるパルスシーケンスはこれを解決したもので、第７図
に示すように９００パルスの侵に１８０”パルスをｎ回
繰返し印加してｎ個のエコーを発生させるようにしたパ
ルスシーケンスである。

このような手法による従来の方法においては次のような
欠点があった。

（１）Ｔ＋象とＴ２像がそれぞれ別個に求められており
、Ｔ＋　Ｎ　Ｔ２　Ｎρの計算画像が同時に得られない
。

（２）近似式を用いているため正確な値が求まらない。

（３）Ｔｒ＞Ｔ、の条件のためＴｒを艮くしなければな
らず、全スキャンタイムが長い。

（４）原画−のスキャンパラメータが最適化されておら
ず、与えられた条件下で最良の計算画像が求まらない。

（５）スライス形状の影響のため系統誤差がある。

本発明の目的は、この様な点に鑑み、ＮＭＲ県像装置に
おいて、複数枚の原画像から短時間に高精度なＴ　、＊
　Ｔ　２　ｅρ計算画像を同時に求め得るようなＮＭＲ
撮像装置を提供することにある。

（問題点を解決するための手段） この様な目的を達成するために本Ｒ明では、核磁気共鳴
ｌ１ｉ）ｅ装置において、計算画像用の新しいパルスシ
ーケンスを考案し、ＩＲｎＳＥ法とＦＲ撮像Ｅ法（ＦＲ
はＦ　ａｓｔ　Ｒｅｃｏｖｅｒｙ法（７）１８）ｋ：Ｊ
：り撮像し、この撮像から（ｎ十ｍ）枚の画像を求め、
この画像からＴ　Ｉ　＊　Ｔ　２　ｇρ像を計算するも
のとし、かつこの場合のＩＲｎＳＥ法とＦＲｍＳＥ法に
おける各スキャンパラメータは、注目する範囲のＴ　１
＋　Ｔ　２１ρ計算画像の評価関数が最良となるように
選定されていることを特徴とする。

〈実施例） 以下図面を用いて本発明の詳細な説明する。第１図は本
発明に係るＮＭＲ１ｉ像装置の一実施例を示す要部構成
図である。図において、１はマグネットアセンブリで、
内部には対象物を挿入するための空間部分（孔）が設け
られ、この空間部分を取巻くようにして、対客物に一様
静磁場Ｈｏを印加する主磁場コイル２と、勾配磁場を発
生するための勾配磁場コイル３（個別に勾配磁場を発生
することができるように構成されたＸ勾配磁場コイル、
ｙ勾配！１ｓコイル、２勾配磁場コイルより構成される
）と、対家物内の原子核のスピンを励起するための高周
波パルスを与えるＲＦ送信コイル４と、対象物からのＮ
ＭＲ信号を検出する受信用コイル５Ｗ−が配置されてい
る。

主磁場コイルは静磁場１１．Ｉｊ　′ｗＪ回路１５に、
Ｇｘ。

Ｇｙ、Ｇｚ各勾配磁場コイルは勾配磁場ちり御回路１４
に、ＲＦ送信コイルは電力増幅器１８に、モしてＮＭＲ
信号の受信用コイルはプリアンプ１９に、それぞれ接続
されている。

１３はコントローラで、勾配磁場や高周波磁場の発生シ
ーケンスを制御すると共に得られたＮＭＲ信号を波形メ
モリ２１に取込むために必要な制御を行う。

１７はゲート変調回路、１６は高周波信号を発生する高
周波発振器である。ゲート変調回路１７は、コントロー
ラ１３からの制御信号により高周波発振器１６が出力し
た高周波信号を適宜に変調し、所定の位相の高周波パル
スを生成する。この高周波パルスはＲＦｉｔ力増幅＠１
８を通してＲＦ送信コイル４に加えられる。

１つは検出コイル５から得られるＮＭＲ信号を増幅する
プリアンプ、２０は高周波発振器の出力信号を参照して
ＮＭＲ信号を位相検波する位相検波回路、２１は位相検
波されたプリアンプからの波形信号を記憶する波形メモ
リで、ここにはＡ／Ｄ変換器を含んでいる。

１１は波形メモリ２１からの信号を受け、所定の信号処
理を施してｌ１ｌｉＩ像を得るコンピュータ、１２は得
られた断層像を表示するテレビジョンモニタのような表
示器である。

３０は操作卓で、コンピュータ１１と連結され、本装置
に必要な各種の情報を入力するための入力手段である。

この様な構成における計算画像作成の手順について次に
説明する。

ここでは、本発明による新規なＩＲｎＳＥ法とＦＲｍＳ
Ｅ法のパルスシーケンスを用い、ＩＲ３ＳＥ法で求めた
３画像と、ＦＲ４ＳＥ法により求めた４画像の計７画像
から、ＴＩＮＴ２、ρ像を計算する場合を例にとって説
明する。

なお、ｌＲａ５Ｅ法は、第２図に示すように、ＩＲ８Ｅ
法に準するパルスシーケンスであるが、一つのビューに
おいて、　　１８０’″（２）パルスを繰返し３回印加
して、３つのエコー信号を得るようにした方式である。

また、ＦＲ４ＳＥ法とは、第３図に示すようなパルスシ
ーケンスを採用した方式である。すなわち、５Ｒ８Ｅ法
において最後（この場合は、第４番目）のエコー信号を
得た後、１８００パルスの印加と、更にそれに続り９０
°パルスと１８０°パルスの印加により磁化を強制的に
主磁場方向へ向け、磁化の平衡状態への回復時間を短縮
するようにしたことを特徴とするパルスシーケンスによ
りエコー信号を得るマルチエコー法である。

（１）パルスシーケンス 第２図及び第３図のパルスシーケンスについて更に詳し
く説明すれば次の通りである。

９０°パルスと第１エコー信号の中心までの間隔Ｔ　Ｂ
　Ｉ　、第１エコー信号以後の各エコー信号の中心間隔
下ｓ２、第４エコー信号の中心とその後に与えられる９
０°パルスまでの間隔Ｔ　Ｂ　３　、第１の１８０°パ
ルス（インバージョン・リカバリ用の１８０°パルス）
から９０１パルスの印加までの間隔Ｔｄ、繰返し時ＩＴ
丁はそれぞれ任意に選ぶことができる。これらの時間管
理はコントローラ１３で行われ、その時間設定は操作車
３０を使用して行うことができる。

１８０°〔１ンパルスはスピン反転用の１８０°パルス
、１′８０°（コ）はスピンエコー用１８０＠パルスで
、パルス誤差を小さくするためにどちらも９０　”−４
５・　２７０°愕　・９０″−付のコンポジット・パル
スを使用している。

ｌＲａ５Ｅ法では各ビューごとの１８０°パルス数は偶
数である。

なお、各パルスの度数に付したサフィックス値は励起用
９０°パルスとの位相差を表わし、これらのパルスは非
選択パルスである。

励起用９０′パルスは、選択パルスであり、ガウシアン
変調されたものである。

このような９０°パルスないし１８００パルスの印加は
次のようにして行われる。すなわら、コントローラ１３
の制御のもとにゲート変調回路１７を通して得た所定の
９０°パルス又は１８０°パルス信号を電力増幅器１８
を介してＲＦ送信コイル４に与え、対象物に印加するＲ
Ｆ磁場を発生させる。

他方、勾配磁場については次の通りである。Ｘ方向の勾
配磁場Ｇｘは、プロジェクション勾配で、ａは１８０°
パルスによるスライス面外のノイズを消去するためのス
ポイラである。

２方向勾配磁場Ｇｚはスライス勾配、ｙ方向勾配磁場Ｇ
ｙはワーブ勾配で、ｂは１８０°パルス誤差によるアー
ティファクトを消去するためのスポイラである。また、
Ｃはビュー間の相関を取除くためのスポイラである。

各勾配ｍ場の印加はコントローラ１３により制御される
。

上記のようなパルスシーケンスにより発生する各エコー
信号は受信コイル５で検出される。受信コイルで検出さ
れたスピンエコー信号は、プリアンプ１９１位相検波回
路２０を経て波形メモリ２１に蓄えられる。

（２）信号強度式について 十丁’−＋／Ｔ＋　十３　Ｔｚ／２Ｔ＋　）　−２−ｅ
ｘ　Ｐ　（−Ｔｉ−７丁を十Ｔ５＋／丁＋　＋Ｔｓｚ　
／　Ｚ　Ｔ　＋）　＋２６ｐ’　（−７１，／７　（と
して、第１エコーは、 １ｏ　　・ｅＸｌ］　　（Ｔｓ　＋　／Ｔ２　）第２エ
コーは、 １ｏ　・ｅＸｒ）　（−Ｔｓ　＋　／Ｔ２−ＴＩ２　／
Ｔ２　）第３エコーは、 １ｏ　−ｅＸＩ）　（Ｔｓ　＋　／Ｔ２−２７！Ｂ　２
　／Ｔ２　）である。

ここでＣＴＦＺ３　　はスライス形状の影響を表わす関
数で、次のように求められる。

スライス形状の影響を含まない信号強度式を二）Ｃ−は
Ｔ　＋　／　Ｔ　ｒのみの関数なので、必要なＴ　＋　
／　Ｔ　ｒの範囲で数値積分によりＣｏｄ誌を求め、こ
の値からＣ０Ｕ　をＴ　＋　／　Ｔ　ｒの多項式として
求めることができる。

以上から、スライス形状の影響を含んだ信号強度式Ｆ　
ｓ　（ＴｒｚＴｓ、Ｔ＋　２ＴＺ／　Ｐ　）は、Ｆｎと
、スライス形状の影響を表わす係数００走と　との積と
して求まる。

Ｆｓ（Ｔ？ｚ丁’ｈｌ　’Ｔｉ　／　Ｔｚ　／　Ｐ　）
−Ｆｎ　　（”’５−Ｔｓ、ＴｕＴｇ−）ＣｏｔＬｄ　
　　（Ｔ＋　　／　Ｔｒ　　）、ここで、Ｃｏａｄ　　
は、例えば０．２＜’Ｔ　＋　／　Ｔ　ｒ　＜１０．０
の場合には Ｃ，ム　−　８．１５３７　Ｅ　−６（Ｔ＋　／Ｔｒ　
）　’−２，９５０８８Ｅ　−４（Ｔ＋　／Ｔｒ　）　
’＋　４．２７６７５Ｅ　−３（Ｔ＋　／Ｔｒ　）　’
−３，１７９０２Ｅ　−２（Ｔ　＋　／　Ｔｒ　）　’
＋　１．２９２６２Ｅ　−ｊ　（Ｔ＋　／Ｔｒ　）　２
−２．８５５４２−１　（Ｔ　Ｉ／　Ｔ、？　）＋　１
．０５５７ パルスシーケンスが１つの９０”パルスと偶数個Ｆ　ｎ
　（Ｔｒ、Ｔｓ＞Ｔ１、ＴＬ）＃）とする。

イ〉　磁化の倒れる角度がＣ１のときの信号強度はパル
スシーケンスが１個の９０°パルスと奇数個の１８０°
パルスから構成されたものである場合には、 となる。

口）ガウシアン９０°パルスを用いていれば、スライス
中央から距離２の点でのαは 医＝（π／ｚ　）−ｅｘ　Ｐ　（−Ｚリ　　　　　　　
　　　・口（乙）となる。

ハ）　（１）式を（２）式により２で積分すればスライ
ス形状の影響を含んだ信号強度が求まり、次式となる。

・・・（３） （３）式の積分は（Ｔ＋／Ｔｒ）のみの関数であるので
、この値をＣ０ｕ（Ｔ＋／Ｔｒ）と書く。

の１８０°パルスから構成されたものである場合には、
スライス形状の彰菅を含まない信号強度式をＦｎとすれ
ば、磁化が倒れる角度がＣ０のときの信号強度は となり、以下上述の場合と同様に計算可能である。

例えば、ガウシアン９０’″パルスを用いていれば、０
．２＜　Ｔ　Ｉ／　Ｔ　２　＜　１０．０で、スライス
形状の影響を表わす係数Ｃ４Ｖ４Ｎは Ｃ４ｖｅａ−−２，４２０３Ｅ　−５（Ｔ＋　／Ｔｒ　
）　’＋　５．６８６１　Ｅ　−４（Ｔ＋　／Ｔｒ　）
　’−３，６５２３Ｅ　−３（Ｔ＋　／Ｔτ）８１−１
．００７１　Ｅ　−２＜Ｔ＋　／Ｔｒ　）　’＋　３．
２１６２　Ｅ　−１（Ｔ＋　／　Ｔｒ　）＋　０．９１
７８ である。

以上がガウシアン９０°パルスを用いた場合の計棹であ
るが、他の９０”パルスを用いた場合でもスライス中央
から距離２の点での９０°パルスにより磁化が倒れる角
度αが求まれば同様に計算できる。

■　ＦＲ４ＳＥ法の信号強度式は 工ｏ　＝　［１−ｅｔｙ　（−（Ｔｒ　−Ｔｓｒ　迅ｚ
−Ｔｓｉ）／Ｔ＋ロリ／［１−σｉ’（−（Ｔｈ−Ｔぎ
１−３丁≦２−Ｔｓｔ）／Ｔ＋−（Ｔｓけ？Ｔｓｚ＋”
Ｔｈ３’）／Ｔ４として、第１エコーは、 Ｎｏ　−ｅＸｐ　　（−Ｔｓ　＋　／Ｔ２　）第２エコ
ーは、 Ｉｏ　−ｅＸｌｌ　　（−Ｔｇ　＋　／Ｔ２　　Ｔｇ　
２　／Ｔ２　）第３エコーは、 ｉｏ　−ｆ３ＸＩ）　　（Ｔｓ　＋　／Ｔ２−２Ｔ！！
　２　／Ｔ２　）第４エコーは、 Ｉｏ−ｅＸｐ（−Ｔｓ　＋　／Ｔ２　２Ｔ！＄　２　／
Ｔ２　）である。

（３）　　スキャンパラメータの最適化人体のＴＩＩＴ
２１ρ計算画像の評価関数が最良となるスキャンパラメ
ータをＷ４差伝播の法則により計算する。信号強汝式に
は前記（２）の信号強度式を用いる。

ここで、信号強度の理論式と、求めるＴ１、Ｔｇ、ρ値
とから計算ｉｉｉ像の評価関数を最良にするスキャンパ
ラメータを求める手法について説明する。ここでは評価
関数として正規化した標準偏差の和、すなわち、 〕　ση／　Ｔ　＋＋σ７ｚ　／　Ｔ　２＋σ４／ρた
だし、六１．σＴ２　、σ、はＴ　Ｉ　＊　Ｔｇ　＊　
Ｉ）の標準偏差 を用いる。

７つの画像のスキャンパラメータをｆｌ　／　Ｐｚｚ・
・°／ＰＴｓ信号強度式をＦ　Ｉ　＊　Ｆ　２　＊　−
−、ｒ　Ｆ　Ｔとすれば、画像から最小２乗法によりｆ
ｆｆ算したＴ＋。

Ｔ２ｅρの値の共分散行列’Ｔ＋　Ｔｚ　Ｐは■工＋Ｘ
　Ｐ　−（Ａ”　■Ｉ□３　Ａ＞−’ただし、ｖ１２３
　　は原画像の共分散行列で、原画像の分散σ２は、平
均値をｎ、サンプリング時間をＴａとして、σ２ｏｃｎ
−’　Ｔａ−’　ｔ’表わされ、まとなる。

したがってＮ　Ｔ　Ｉ　＋　７２　＊ρの値の分散はＶ
ＴＩＴＺｐの対角要素として求まる。

以上から、計算画像の評価関数がＰＩ　／Ｐｚ　ｔ・・
・ｆ’７　／Ｔ（Ｌｌ／ＴｌＬＸ度−／丁’ｃｂｑｒ％
＋ｉ’ｒしｌ　ｔ　＋　Ｌ　　Ｌ　／　りＬＹの関数と
して求まる。

このような原理に基づき、次のような手順により適切な
スキャンパラメータが求められる。

■信号強度の理論式を定める。

■理論式と、測定したいＴ＋　、Ｔｇ−ρの範囲と、原
画像の分散から、計算画像の評価関数をスキャンパラメ
ータの関数として求める。

■上記■において計算画像の評価関数がスキャンパラメ
ータの多変数関数として求まったので、多変数関数の極
値を求める方法（シンプレックス法等〉によりｉ！ｌＦ
価関数が最良となるスキャンパラメータを求める。

このようにして求めたスキャンパラメータの一例を示せ
ば次の通りである。

■トータルスキャンタイムを６００秒でスキャンする場
合には ｌＲａ５Ｅ法においては ＴＴ−２，１秒 Ｔｄ＝０．４５７秒 Ｔ、　、　＝０．０２秒 Ｔｓ　２　＝　０．０２秒 平均回数（ＡＶＥ）−１ ＦＲ４ＳＥ法においては Ｔｒ＝１．３秒 Ｔｓ　Ｉ＝０．０２秒 Ｔ、　、　＝０．０７３秒 Ｔｓ　３−０．０２秒 平均回数（ＡＶＥ）−２ ■トータルスキャンタイムは３００秒の場合ＩＲ３ＳＥ
法では Ｔｒｘ　１，２８秒 ７ｄ　＝０．３６６秒 Ｔｓ　＋　＝０．０２秒 Ｔｓ　２−０．０２秒 平均回数（ＡＶＥ）＝　１ ＦＲ４ＳＥ法では Ｔｒ＝１．１秒 Ｔｓ　＋　＝　０．０２秒 Ｔ、　２−　０．０７３秒 Ｔｓ　３　””　０．０２秒 平均回数（ＡＶＥ）−２ （４）　　前記（３）のスキャンパラメータで撮像する
。

すなわち、ＩＲ３ＳＥ法において、ワープ勾配（勾配計
１１ｉＧｙ）の互いに異なる所定のビュー（ビュー数は
例えば１２７）にわたって上記のパラメータでスキャン
し、エコー信号を測定する。測定採取したエコー信号を
第１、第２及び第３エコー信号群ごとに分け、コンピュ
ータ１１を使ってそれぞれ２次元画像に再構成して３枚
の原画像を１ｑる。

次にＦＲ４ＳＥ法において、ワープ勾配（勾配磁ＪｉＧ
ｙ）の互いに異なる所定のビュー（ビュー数は例えば１
２１）にわたって上記のパラメータでスキャンし、同様
にエコー信号を測定し波形メモリ２１に格納する。得ら
れたデータを第１、第２、第３及び第４エコー信号群ご
とに分け、同様にコンピュータ１１を使ってそれぞれ２
次元画像に再構成して４枚の原画像を得る。

（５）　　前記（４）で得られた７枚の原画像を用い、
非線形最小２乗法によりＴ　Ｉ　Ｔ　Ｔ　２１ρ鍮（計
算画像）を求める（コンピュータ１１にて演算により求
める。）。

以上の手法において、勾配１１１１Ｇｘ　、　Ｇｙ　、
　Ｇ２と、スライス、プロジェクション、ワープの関係
は任意である。

以上の手順により、ＴＩ　＊　７２　＊ρの計算画像を
正確かつ同時に得ることができる。

（他の実施例） なお、本発明は上記実施例に限らず次のようにすること
ができる。

（１）パルスシーケンスとして、実施例では、９０゜パ
ルスから第１エコー信号までの間隔をＴ！ｉ１、第１エ
コー信号以後の各エコー信号間隔をＴｓ２としたが、こ
れを、 ９０°パルスから第１エコー信号までの間隔と、第１エ
コー信号から第２エコー信号までの間隔をそれぞれＴｓ
　Ｉ　ｓ第２エコー信号から第３エコー信号までの間隔
及び第３エコー信号から第４エコー信号までの間隔をＴ
ｓ２、それ以後の第（２ｎ−２）エコー信号から第（２
ｎ−１＞エコー信号までの間隔及び第（２ｎ−１）エコ
ー信号から第２ｎエコー信号までの間隔（ただしｎ≧３
）をＴｓｎとする。この場合１８０°パルス誤差を偶数
エコーで除くことができる。

更に、各エコー信号間隔を全く任意にしてもよい。

（２）ガウシアン変調の９０°パルス、コンポジット１
８０°パルスも実施例に限定されるものではなく、他の
９０°パルス、１８０”パルスを用いてもよく、同様に
計算画像を得ることができる。

（３）　　実施例では、ｌＲ３８０法とＦＲ４ＳＥ法に
よる！ｆｉ像から７つの画像を得、これをもとにＴ１゜
Ｔ２＋ρの各計算画像を求める場合を示したが、これに
限定することなく、一般にＩＲｎＳＥ法とＦＲｍＳＥ法
によるｍ像から（ｎ＋ｍ）枚の画像を得、これをもとに
Ｔ　Ｉ　＋　７２　＊ρの各計算画像を求めるようにし
てもよい。

（発明の効果） 以上説明したように、本発明によれば、次のような効果
がある。

■近似式を用いていないため正確な値を求めることがで
きる。

■計算画像の分散を最小とするスキャンパラメータを用
いることにより、所定の撮像時間で最良の計算画像が求
まる。

■ＩＲｎＳＥ法の１８０°パルス数を偶数とすることに
より計算画像の標準偏差を小さくすることができる。第
４図は、種々のパルスシーケンスの組合せについて、人
体のＴＩ、Ｔ２値の代表点での正規化した標準１２差の
総和を示したものである。

このグラフから明らかなように、ＩＲｎＳＥ法（ｎは整
数）の１８０°パルス数が偶数の場合、標準偏差が小さ
いことが分る。すなわち、ｌＲ４５Ｅ法よりもＩＲ３Ｓ
Ｅ法を使用する方が良いことが分る。

■ＦＲ法を用いることにより、他の方法より更に標準偏
差を小さくすることができる。第４図において、ｌＲａ
５Ｅ法、ＦＲ４ＳＥ法が最良である。

■また、ＦＲ法を用いることにより、Ｔ１．Ｔ２の広い
範囲で分散を小さくすることができる。

■１８０°（２）パルスにコンポジット・パルス９０’
−４５２７０°確　・９Ｇ”−４５を使用し、各エコー
信号間の開隔を、隣接する２つずつが等しくなるように
設定することにより、１８０°パルス誤差を偶数エコー
信号で除去することができる。

【図面の簡単な説明】

第１図は本発明に係るＮＭＲＩｌ像装置の一実施例を示
す要部構成図、第２図はＩＲ３ＳＥ法のパルスシーケン
スを示す図、第３図はＦＲ４ＳＥ法のパルスシーケンス
を示す図、第４図は種々のパルスシーケンスの組合せに
ついて、人体のＴＩ＋７２Ｍの代表点での正規化した標
準偏差の総和を示す図、第５図ないし第７図は従来のパ
ルスシーケンスの一例を示す図である。 １・・・マグネットアセンブリ、２・・・主磁場コイル
、３・・・勾配磁場コイル、４・・・ＲＦ送信コイル、
５・・・受信用コイル、１１・・・５ンビユータ、１２
・・・表示器、１３・・・コントローラ、１４・・・勾
配磁場＄り御回路、１５・・・静磁場制御回路、１６・
・・高周波発振器、１７・・・ゲート変調回路、１８・
・・電力増幅器、１９・・・プリアンプ、２０・・・位
相検波回路、２１・・・波形メモリ、３０・・・操作卓
。

Claims (1)

1. 【特許請求の範囲】 １）対象物に高周波パルスおよび磁場を印加して核磁気
   共鳴信号を発生させ、この信号を用いて対象物の組織に
   関する画像を得るようにした核磁気共鳴撮像装置におい
   て、 緩和時間（Ｔ＿１、Ｔ＿２）又はプロトン密度（ρ）の
   少なくともいずれか一つに関する計算画像を得るための
   下記（イ）ないし（ハ）の機能を有する制御・演算手段
   を具備したことを特徴とする核磁気共鳴撮像装置。 記 （イ）注目する範囲のＴ＿１、Ｔ＿２、ρ計算画像の評
   価関数が最良となるＩＲｎＳＥ法（ｎは整数）とＦＲｍ
   ＳＥ法（ｍは整数）のスキャンパラメータを求める。 （ロ）ＩＲｎＳＥ法とＦＲｍＳＥ法とによるパルスシー
   ケンスにより前記（イ）で求めたスキャンパラメータで
   撮像する。 （ハ）ＩＲｎＳＥ法とＦＲｍＳＥ法による撮像から求め
   た（ｎ＋ｍ）枚の画像からＴ＿１、Ｔ＿２、ρ像を計算
   する。 ２）前記制御演算手段は、ＩＲｎＳＥ法においてそのパ
   ルスシーケンス中で１８０°パルスが偶数回印加される
   ように構成されたことを特徴とする特許請求の範囲第１
   項記載の核磁気共鳴撮像装置。 ３）前記制御・演算手段は、ＩＲ３ＳＥ法とＦＲ４ＳＥ
   法により撮像を行い７つの画像を得、この７つの画像か
   らＴ＿１、Ｔ＿２、ρ像を計算するように構成したこと
   を特徴とする特許請求の範囲第１項記載の核磁気共鳴撮
   像装置。 ４）前記制御・演算手段は、ＩＲ３ＳＥ法とＦＲ４ＳＥ
   法により撮像を行い７つの画像を得、この７つの画像か
   らＴ＿１、Ｔ＿２、ρ像を計算し、スキャンパラメータ
   としては下記の値を用いるように構成したことを特徴と
   する特許請求の範囲第１項記載の核磁気共鳴撮像装置。 記 （１）ＩＲ３ＳＥ法では Ｔ＿ｒ＝２．１秒 Ｔ＿ｄ＝０．４６秒 Ｔ＿ｓ＿１＝０．０２秒 Ｔ＿ｓ＿２＝０．０２秒 として所定のビュー数スキャンする。 （２）ＦＲ４ＳＥ法では Ｔ＿ｒ＝１．３秒 Ｔ＿ｓ＿１＝０．０２秒 Ｔ＿ｓ＿２＝０．０７３秒 Ｔ＿ｓ＿３＝０．０２秒 として所定のビュー数スキャンする。 （３）求める７つの画像は、ＩＲ３ＳＥ法とＦＲ４ＳＥ
   法を１対２の割合で平均して求める。 ５）前記制御・演算手段は、ＩＲ３ＳＥ法とＦＲ４ＳＥ
   法により撮像を行い７つの画像を得、この７つの画像か
   らＴ＿１、Ｔ＿２、ρ像を計算し、スキャンパラメータ
   としては下記の値を用いるように構成したことを特徴と
   する特許請求の範囲第１項記載の核磁気共鳴撮像装置。 記 （１）ＩＲ３ＳＥ法では Ｔ＿ｒ＝１．２８秒 Ｔ＿ｄ＝０．３７秒 Ｔ＿ｓ＿１＝０．０２秒 Ｔ＿ｓ＿２＝０．０２秒 として所定のビュー数スキャンする。 （２）ＦＲ４ＳＥ法では Ｔ＿ｒ＝１．１秒 Ｔ＿ｓ＿１＝０．０２秒 Ｔ＿ｓ＿２＝０．０７３秒 Ｔ＿ｓ＿３＝０．０２秒 として所定のビュー数スキャンする。 （３）求める７つの画像は、ＩＲ３ＳＥ法とＦＲ４ＳＥ
   法を１対１の割合で平均して求める。