JPH0374101B2 - - Google Patents

Info

Publication number
JPH0374101B2
JPH0374101B2 JP60196661A JP19666185A JPH0374101B2 JP H0374101 B2 JPH0374101 B2 JP H0374101B2 JP 60196661 A JP60196661 A JP 60196661A JP 19666185 A JP19666185 A JP 19666185A JP H0374101 B2 JPH0374101 B2 JP H0374101B2
Authority
JP
Japan
Prior art keywords
images
pulse
image
scan
magnetic resonance
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.)
Expired - Lifetime
Application number
JP60196661A
Other languages
English (en)
Other versions
JPS6264348A (ja
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed filed Critical
Priority to JP60196661A priority Critical patent/JPS6264348A/ja
Priority to US06/833,224 priority patent/US4694250A/en
Priority to DE19863606220 priority patent/DE3606220A1/de
Priority to GB08604693A priority patent/GB2173598B/en
Publication of JPS6264348A publication Critical patent/JPS6264348A/ja
Publication of JPH0374101B2 publication Critical patent/JPH0374101B2/ja
Granted legal-status Critical Current

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Description

【発明の詳細な説明】 (産業上の利用分野) 本発明は、核磁気共鳴撮像装置(以下核磁気共
鳴をNMRと略す)に関し、特にスキヤン方法お
よび診断の際のデータ処理の方法に関する。
(従来の技術) NMR画像の特徴の一つは、その濃度がプロト
1H密度(以後ρと記す)だけでなく、緩和時
間T1、T2や更にはパルス系列およびRFパルス間
隔等のスキヤンパラメータの関数となつていて、
測定パラメータを様々に変えて種々のコントラス
トの画像が得られ、疾患に関する多様な情報が得
られることにある。
(発明が解決しようとする問題点) しかし、逆に臨床データ間の相互比較が難しく
なり、画像に客観性がなくなる点があることも否
めない。
現在行なわれている撮像法は、第9図のフロー
チヤートに示されているように、目的とする疾患
による画像コントラストが大きいパルス系列やス
キヤンパラメータを経験的に選んでいる。また、
臨床例のない疾患に対しては、第10図のフロー
チヤートに示されているように、撮像条件を変え
て複数の画像を得、画像コントラストがつく撮像
条件と疾患を経験的に対応付けている。このよう
な臨床データの蓄積方法は非能率であるばかりで
なく、実際の臨床では、多数の撮像条件でスキヤ
ンしなければならないこと、更に過去の臨床デー
タのない疾患では最適コントラストが得られるス
キヤンができない等、実用上大きな問題がある。
一方、複数のNMR像から画像間演算を行い、
画像に客観的に意味のある純粋なT1、T2像を得
て、診断に用いるという試みがある。しかしT1
T2値からだけでは総べての疾患が診断できない
こと、また現在のT1、T2計算画像はS/N比、
T1、T2値の精度が悪く、実用に耐えないこと等
大きな問題点がある。
本発明の目的は、この様な点に鑑み、あらゆる
疾患に対して同一のスキヤン条件で撮像し、スキ
ヤン時間の短縮化を図るようにする一方、得られ
る画像の物理的意味が明確であると共に診断に有
効な画像コントラストがあり客観性のある画像を
得ることのできる核磁気共鳴撮像装置を提供する
ことにある。
(問題点を解決するための手段) この様な目的を達成するために本発明では、対
象物に高周波パルスおよび磁場を印加して核磁気
共鳴信号を発生させ、この信号を用いて対象物の
組織に関する画像を得るようにした核磁気共鳴撮
像装置において、 あらかじめ設定された最適なパルス系列とスキ
ヤンパラメータでスキヤンし、得られたデータか
らT1、T2、ρ計算画像を計算し、次いでパルス
系列の選択とスキヤンパラメータの設定を行うた
めの操作卓の入力に応じ診断に都合のよい画像を
再合成する制御・演算手段を具備したことを特徴
とする。
(実施例) 以下図面において本発明を詳しく説明する。第
1図は本発明に係るNMR撮像装置の一実施例を
示す要部構成図である。図において、1はマグネ
ツトアセンブリで、内部には対象物を挿入するた
めの空間部分(孔)が設けられ、この空間部分を
取巻くようにして、対象物に一様静磁場H0を印
加する主磁場コイル2と、勾配磁場を発生するた
めの勾配磁場コイル3(個別に勾配磁場を発生す
ることができるように構成されたx勾配磁場コイ
ル、y勾配磁場コイル、z勾配磁場コイルより構
成される)と、対象物内の原子核のスピンを励起
するための高周波パルスを与えるRF送信コイル
4と、対象物からのNMR信号を検出する受信用
コイル5等が配置されている。
主磁場コイルは静磁場制御回路15に、Gx,
Gy,Gz各勾配磁場コイルは勾配磁場制御回路1
4に、RF送信コイルは電力増幅器18に、そし
てNMR信号の受信用コイルはプリアンプ19
に、それぞれ接続されている。
13はコントローラで、勾配磁場や高周波磁場
の発生シーケンスを制御すると共に得られた
NMR信号を波形メモリ21に取込むために必要
な制御を行う。
17はゲート変調回路、16は高周波信号を発
生する高周波発振器である。ゲート変調回路17
は、コントローラ13からの制御信号により高周
波発振器16が出力した高周波信号を適宜に変調
し、所定の位相の高周波パルスを生成する。この
高周波パルスはRF電力増幅器18を通してRF送
信コイル4に加えられる。
19は検出コイル5から得られるNMR信号を
増幅するプリアンプ、20は高周波発振器の出力
信号を参照してNMR信号を位相検波する位相検
波回路、21は位相検波されたプリアンプからの
波形信号を記憶する波形メモリで、ここにはA/
D変換器を含んでいる。
11は波形メモリ21からの信号を受け、所定
の信号処理を施して断層像を得るコンピユータ、
12は得られた断層像を表示するテレビジヨンモ
ニタのような表示器である。
22は高速画像演算装置で、パイプライン方式
や並列演算により例えば数分の1秒で画像の演算
が完了するものである。なお、この高速画像演算
装置22はコンピユータ11側に含め、一体的な
構成としてもよい。
30は操作卓で、高速画像演算装置22と連結
され、本装置に必要な各種の情報、例えばパルス
系列、RFパルス間隔等のスキヤンパラメータ、
その他を入力するための入力手段である。
この様な構成における計算画像作成の手順につ
いて次に説明する。
ここでは、IR3SE法で求めた3画像と、
SR4SE法により求めた4画像の計7画像から、
T1、T2、ρ像を計算する場合を例にとつて説明
する。
なお、IR3SE法は、第2図に示すように、
IRSE法{反転回復(Inversion Recovery)法と
スピンエコー(Spin Echo)法を組合せたパルス
シーケンス}に準ずるパルスシーケンスであり、
一つのビユーにおいて、180°パルスを繰返し3回
印加して、3つのエコー信号を得るようにした方
式である。
また、SR4SE法とは、第3図に示すように、
SRSE法{飽和回復(Saturation Recovery)法
とスピンエコー(Spin Echo)法を組合せたパル
スシーケンス}に準じたものであり、180°パルス
を4回印加し1ビユーつき4個のエコー信号を得
るマルチエコー法である。
以下第4図に示すスキヤン手順を参照し、主と
して本発明の特徴とする動作について説明する。
(A) あらかじめ最適化されたパルス系列およびス
キヤンパラメータを設定しておく、本装置は設
定されたパルス系列およびスキヤンパラメータ
でスキヤンするようにプログラムされており、
いかなる疾患を撮影する目的でも同一スキヤン
で行われる。
ただし、所望の撮像部位、例えば頭部、胴、
脚等に関連して最適なパルス系列のスキヤンパ
ラメータは変り、パルス系列のスキヤン条件は
以下のようにして決定される。
パルス系列としては、撮像したい部分のT1
T2、ρの範囲で、撮像時間が短く、高S/N
比、高精度のT1、T2、ρ計算画像が得られる
IR3SE法とSR4SE法とを組合せたパルス系列
を用いるのが望ましい。そのIR3SE法と
SR4SE法との組合せの場合の最適なパラメー
タは次のようにして求められる。
(1) パルスシーケンスについて 第2図及び第3図のパルスシーケンスにお
いて、90°パルスと第1エコー信号の中心ま
での間隔Ts1、第1エコー信号以後の各エコ
ー信号の中心間隔Ts2、180°(1)パルスと90°パ
ルスまでの間隔Td、繰返し時間Trはそれぞ
れ任意に選ぶことができる。これらの時間管
理はコントローラ13で行われる。
180°(1)パルスは、スピン反転用の180°パル
ス、180°(2)はスピンエコー用180°パルスで、
パルス誤差を小さくするためにどちらも
90°-45・270°45・90°-45のコンポジツト・パル
スを使用している。
各ビユーごと180°パルス数はIR3SE法、
SR4SE法共に偶数である。
なお、コンポジツトパルスの度数に付した
サフイツクス値は励起用90°パルスとの位相
差を表わし、これらのパルスは非選択パルス
である。
励起用90°パルスは、選択パルスであり、
ガウシアン変調されたものである。
このような90°パルスないし180°パルスの
印加は次のようにして行われる。すなわち、
コントローラ13の制御のもとにゲート変調
回路17を通して得た所定の90°パルス又は
180°パルス信号を電力増幅器18を介して
RF送信コイル4に与え、対象物に印加する
RF磁場を発生させる。
他方、勾配磁場については次の通りであ
る。x方向の勾配磁場Gxは、プロジエクシ
ヨン勾配で、aは180°パルスによるスライス
面外のノイズを消去するためのスポイラであ
る。
z方向勾配磁場Gzはスライス勾配、y方
向勾配磁場Gyはワープ勾配で、bは180°パル
ス誤差によるアーテイファクトを消去するた
めのスポイラである。また、cはビユー間の
相関を取除くためのスポイラである。
各勾配磁場の印加はコントローラ13によ
り制御される。
上記のようなパルスシーケンスにより発生
した各エコー信号は受信コイル5で検出され
る。受信コイルで検出されたスピンエコー信
号は、プリアンプ19、位相検波回路20を
経て波形メモリ21に蓄えられる。
(2) 信号強度式について IR3SE法の信号強度式は I0=CIR3(T1/Tr){1−2・exp(−Td/T1
+2・exp(−Tr/T1 +Ts1/T1+3Ts2/2T1)−2・exp(−Tr/T1
+Ts1/T1+Ts2/2T1) +2・exp(−Tr/T1+Ts1/2T1)−exp(−Tr
/T1)}・ρ として、第1エコーは、 I0・exp(−Ts1/T2) 第2エコーは、 I0・exp(−Ts1/T2−Ts1/T2) 第3エコーは、 I0・exp(−Ts1/T2−2Ts1/T2) である。
SR4SE法の信号強度式は I0=CSR4(T1/Tr){1−2・exp(−Tr/T1+T
s1/T1+5Ts2/2T1) +2・exp(−Tr/T1+Ts1/T1+3Ts2/2T1) −2・exp(−Tr/T1+Ts1/T1+Ts2
2T1) +2・exp(−Tr/T1+Ts1/2T1)−exp(−Tr
T1)}・ρ として、第1エコーは、 I0・exp(−Ts1/T2) 第2エコーは、 I0・exp(−Ts1/T2−Ts2/T2) 第3エコーは、 I0・exp(−Ts1/T2−2Ts2/T2) 第4エコーは、 I0・exp(−Ts1/T2−3Ts2/T2) である。
ここで、CIR3、CSR4はスライス形状の影
響を表わす関数で、次のように求められ
る。
スライス形状の影響を含まない信号強度
式をFo(Tr、Ts、T1、T2、ρ)とする。
(イ) 90°パルスにより磁化の倒れる角度がα°のと
きの信号強度はパルスシーケンスが1個
の90°パルスと奇数個の180°パルスから
構成されたものである場合には、 sinα/1+cosα・exp(−Tr/T1
・Fo(Tr、Ts、T1、T2、ρ)……(1) となる。
(ロ) ガウシアン90°パルスを用いていれば、スラ
イス中央から距離zの点でのαは α=(π/2)・exp(−Z2) ……(2) となる。
(ハ) (1)式を(2)式によりzで積分すればスライス形
状の影響を含んだ信号強度が求まり、次
式となる。
Fo(Tr、Ts、T1、T2、ρ)∫sin{(π/2)exp
(−Z2)}/1+cos{(π/2)exp(−Z2)}exp(
−Tr/T1)dZ……(3) (3)式の積分は(T1/Tr)のみの関数
であるので、この値をCpdd((T1/Tr
と書く。
(ニ) CpddはT1/Trのみの関数なので、必
要なT1/Trの範囲で数値積分により
Cpddを求め、この値からCoddをT1/Tr
の多項式として求めることができる。
以上から、スライス形式の影響を含ん
だ信号強度式Fs(Tr、Ts、T1、T2、ρ)
は、Foと、スライス形状の影響を表わ
す係数Cpddとの積として求まる。
Fs(Tr、Ts、T1、T2、ρ)=Fo(Tr、Ts
T1、T2、ρ)Cpdd(T1/Tr) ここで、Cpddは、例えば0.2<T1/Tr
<10.0の場合には Cpdd=8.1537E−6(T1/Tr6−2.95086E−4(T1
Tr5+4.27675E−3(T1/Tr4 −3.17902E−2(T1/Tr3+1.29262E−1(T1/Tr
2−2.8554E−1(T1/Tr)+1.0557 パルスシーケンスが1つの90°パルス
と偶数個の180°パルスから構成されたも
のである場合には、スライス形状の影響
を含まない信号強度式をFoとすれば、
磁化が倒れる角度がα°のときの信号強度
は sinα/1−cosα・exp(−Tr/T1)・Fo……(4) となり、以下上述の場合と同様に計算可
能である。
例えば、ガウシアン90°パルスを用い
ていれば、0.2<T1/T2<10.0で、スラ
イス形状の影響を表わす係数Ceveoは Ceveo=−2.4203E−5(T1/Tr5+5.6861E−4(T1
/Tr4−3.6523E−3(T1/Tr3 −1.0071E−2(T1/Tr2+3.2162E−1(T1/Tr
+0.9178 である。
ここで求めたスライス形状の影響を含
んだ信号強度式を用いることにより、ス
ライス形状による系統誤差を除去するこ
とができる。
以上がガウシアン90°パルスを用いた
場合の計算であるが、他の90°パルスを
用いた場合でもスライス中央から距離z
の点での90°パルスにより磁化が倒れる
角度αが求まれば同様に計算できる。
(3) 90°パルス誤差による系統誤差について 90°パルスがRFパルス付均一等のためβ°パ
ルスとなつた場合、得られる信号強度は前記
(2)の場合と同様に(1)式又は(4)式を(5)式により
zで積分すれば求まる。
α=β・exp(−z2) ……(5) この影響を表わす係数をC〓とすれば、
Cpdd・Ceveoと同様にC〓は(T1/Tr)のみの
関数となり、またその形は180°パルス数の奇
数により決まる。
IR3SEとSR4SEの180°パルス数は共に偶数
のため、Trを等しくすれば、IR3SEと
SR4SEのC〓は一致する。したがつて、90°パ
ルス誤差の影響は総べての画像でC〓をかける
こととなり、T1、T2、ρ計算値へと影響は
ρ値にC〓がかかるだけとなる。
以上から、IR3SEとSR4SEのTrを等しく
すれば、T1、T2値への90°パルス誤差による
系統誤差を除くことができる。
Trを等しくしたことによる他の効果とし
て次のことがある。スライス形状の影響を表
わす係数CeveoもIR3SEとSR4SEで一致する
ため、スライス形状による系統誤差はやはり
ρ値にCeveoがかかるだけとなる。したがつ
て、ρ値を精度良く求める必要がなければ、
スライス形状の影響を含まない通常の信号強
度式を用いてもよい。またスライス形状が不
明の場合にも有効である。
(4) スキヤンパラメータの最適化 人体のT1、T2、ρ計算画像の評価関数が
最良となるスキヤンパラメータを誤差伝播の
法則により計算する。信号強度式には前記(2)
の信号強度式を用いる。
ここで、信号強度の理論式と、求まるT1
T2、ρ値とから計算画像の評価関数を最良
にするスキヤンパラメータを求める手法につ
いて説明する。ここでは評価関数として正規
化した標準偏差の和、すなわち、 ρT1/T1+σT2/T2+σ〓/ρ ただし、σT1、σT2、σ〓はT1、T2、ρの標準
偏差 を用いる。
7つの画像のスキヤンパラメータを〓1
2、…、〓7、信号強度式をF1、F2、…、F7
とすれば、画像から最小2乗法により計算し
たT1、T2、ρの値の共分散行列VT1T2〓は VT1T2〓=(ATV123 -1A)-1 ただし、V123は原画像の共分散行列で原画
像の分散σ2は、平均値をn、サンプリング時
間をTaとして、σ2∝n-1Ta -1で表わされ、ま
たAは、 となる。
したがつて、T1、T2、ρの値の分散はVT1T2
の対角要素として求まる。
以上から、計算画像の評価関数が〓1、〓2
…、〓7、Ta1、Ta2、…、Ta7、n1、n2、…、n7
関数として求まる。
このような原理に基づき、次のような手順によ
り適切なスキヤンパラメータが求められる。
信号強度の理論式を定める。
理論式と、測定したいT1、T2、ρの範囲と、
原画像の分散から、計算画像の評価関数を
スキヤンパラメータの関数として求める。
上記において計算画像の評価関数がスキヤ
ンパラメータの多変数関数として求まつた
ので、多変数関数の極値を求める方法(シ
ンプレツクス法等)により評価関数が最良
となるスキヤンパラメータを求める。
このようにして求めたスキヤンパラメータ
の一例を示せば次の通りである。なお、
IR3SE法とSR4SE法でのTrは等しくしてい
る。
トータルスキヤンタイムを600秒でスキヤン
する場合には IR3SE法においては Tr=2.36秒 Td=0.56秒 Ts1=0.02秒 Ts2=0.02秒 平均回数(AVE)=1 SR4SE法においては Tr=2.36秒 Ts1=0.02秒 Ts2=0.079秒 平均回数(AVE)=1 また、トータルスキヤンタイムは300秒
でスキヤンする場合には IR3SE法では Tr=1.18秒 Td=0.41秒 Ts1=0.02秒 Ts2=0.02秒 平均回数(AVE)=1 SR4SE法においては Tr=1.18秒 Ts1=0.02秒 Ts2=0.074秒 平均回数(AVE)=1 (5) 前記(4)のスキヤンパラメータで撮像する。
すなわち、IR3SE法において、ワープ勾配
(勾配磁場Gy)の互いに異なる所定のビユー
(ビユー数は例えば127)にわたつて上記のパ
ラメータでスキヤンし、エコー信号を測定す
る。測定採取したエコー信号を第1、第2及
び第3エコー信号群ごとに分け、コンピユー
タ11を使つてそれぞれ2次元画像に再構成
して3枚の原画像を得る。
次にSR4SE法において、ワープ勾配(勾
配磁場Gy)の互いに異なる所定のビユー
(ビユー数は例えば127)にわたつて上記のパ
ラメータでスキヤンし、同様にエコー信号を
測定し波形メモリ21に格納する。得られた
データを第1、第2、第3及び第4エコー信
号群ごとに分け、同様にコンピユータ11を
使つてそれぞれ2次元画像に再構成して4枚
の原画像を得る。
(B) 前記(5)で得られた7枚の原画像を用い、画像
間演算を行い非線形最小2乗法によりT1、T2
ρ計算画像を得る(コンピユータ11にて演算
により求める。)。
なお、以上の手法において、勾配磁場Gx
Gy、Gzと、スライス、プロジエクシヨン、ワ
ープの関係は任意である。
このようにして各画素に相当する部位の精密
なT1、T2、ρ値が分かつたことになる。
(C) 操作卓30は再合成画像のためのパラメータ
入力部を持つており、ここからパラメータを入
力することができる。高速画像演算装置22は
任意のパルス系列やスキヤンパラメータの画像
をT1、T2、ρ値を用いて高速に演算し、得ら
れた画像は表示装置12で表示される。
例えば、まずIRSE法のパルス系列を選択す
る。次にスキヤンパラメータであるTd、Ts
Trの設定ボタンを適当時間押圧する。押圧時
間長さに対応して各Td、Ts、Trが設定され、
その設定値の再合成画像が表示される。Tr
ボタンを更に押し続けると、Trの変化に伴い
画像コントラストが変化することが分り、診断
に最適なコントラストの画像をリアルタイムで
見出すことができる。
また、Trが無限大等のように実際には実現
し得ない条件での画像や、非常に長い撮像時間
を要する画像も容易に得ることができる。この
ときの演算において、IRSE法ではその信号強
度の理論式として IIR=ρ・exp(−Ts/T2){1−2・exp(−Td
/T1)+2・exp(−Tr/T1+Ts/2T1) −exp(−Tr/T1)} を用いる。他のパルスシーケンス例えばSR法、
SE法等についても、理論式を前もつてプログ
ラムしておけば、同様に再構成画像を得ること
ができる。
更にパルス系列とは全く関係のない架空の画
像を再合成し、疾患を識別することもできる。
例えば、フイツシヤー(Fisher)の識別関数を
最大にするT1、T2、ρの線形結合画像を用い
ることもできる。更にこの場合T1、T2、ρの
各押しボタンがあり、押圧時間に対応した重み
で線形結合することができる。
(他の実施例) 本発明は上記実施例に限定されるものではなく
次のようにすることもできる。
(1) 撮像条件(パルス系列、スキヤンパラメー
タ)は、短い撮像時間で、高S/N比、高精度
なT1、T2、ρ計算画像が得られるものであれ
ばよい。なぜなら、結果として得られる再合成
画像が従来法のスキヤン画像より高品質でなけ
れば意味がないからである。
他のパルス系列としては、IRnSE法(nは整
数)とFRmSE法(FRはFast Recovery法の略
で、FR法とSE法の組合せで1ビユーごとにm
個のエコー信号を得るようにしたもの。なお
FR4SE法のパルスシーケンスは第5図に示す
通りである。)を用いてもよい。また、IR3SE
法とSR4SE法に限らず、一般にIRnSE法と
SRmSE法(n、mは整数)であつてもよい。
第6図には各パルス系列の組合せとT1、T2
ρ係数画像の分散の和について、また第7図に
は各パルス系列におけるT1、T2係数画像の
T1、T2値の分散について、また第8図には各
パルス系列の用いる画像数とT1、T2、ρ計算
画像の分散の和についてそれぞれ示してある。
また、撮像条件はどんな疾患に対しても同一
としたが、更に効果的なスキヤンをするめに
は、いくつかのスキヤンプログラムが準備され
ており、任意に選択可能に構成されていてもよ
い。要するに、診断は総べて再合成画像を用い
たことにある。
(発明の効果) 以上説明したように、本発明によれば、次のよ
うな効果がある。
(1) 総べての疾患に対しても、同一のパルス系
列、スキヤンパラメータで撮像するため、医師
が撮像条件を選択する必要がない。したがつ
て、疾患と撮像条件の関係を知らなくても画一
的な撮像が可能である。撮像条件を誤つたた
め、画像から疾患を識別できず、再スキヤンを
行う等の無駄な操作が不要となる。また、未知
の疾患に対しても、何の配慮も要せず、同様に
撮像することができる。このような点から、本
発明は実用に供してNMR装置の広い普及に大
きな効果がある。
(2) 疾患によつては、従来長い撮像時間を要する
ものがあつた、本発明によれば撮像時間は短縮
される。
(3) 同一の撮像条件を用いるため、画像間に客観
性があり、相互比較することができる。総べて
のNMR像が同一撮像条件であれば、情報交換
やデータ蓄積に画期的な意味がある。
(4) 再合成画像で診断するため、常に最良の画像
コントラストが得られ、診断が迅速かつ正確に
行われる。
(5) 従来の装置では各種のパルス系列やスキヤン
パラメータで自由にスキヤンできるようにして
いたためハードウエアが複雑であつたが、本発
明では唯一の撮像条件(又は少数の撮像条件)
であるため、ハードウエアが簡素化できる。
(6) あらゆる疾患や複合した疾患にも対応でき
る。疾患の種類により、最大コントラストが得
られる撮像条件が異なるため、従来装置では困
難であつたが本発明ではいかなる条件の場合の
画像でも再合成画像により得ることができる。
(7) 画像コントラストを用いた形態診断だけでな
く、高精度なT1、T2、ρ値が得られるので、
これら生体情報から病理学的診断も可能であ
る。
【図面の簡単な説明】
第1図は本発明に係る核磁気共鳴撮像装置の一
実施例を示す要部構成図、第2図はIR3SE法のパ
ルスシーケンスを示す図、第3図はSR4SE法の
パルスシーケンスを示す図、第4図はスキヤン手
順を示す図、第5図はFR4SE法のパルスシーケ
ンスを示す図、第6図は各パルス系列の組合せと
T1、T2、ρ計算画像の分散の和を示す図、第7
図は各パルス系列におけるT1、T2計算画像の
T1、T2値の分散を示す図、第8図は各パルス系
列の用いる画像数とT1、T2、ρ計算画像の分散
の和を示す図、第9図ないし第10図は従来のス
キヤン手順について説明するためのフローチヤー
トである。 1……マグネツトアセンブリ、2……主磁場コ
イル、3……勾配磁場コイル、4……RF送信コ
イル、5……受信用コイル、11……コンピユー
タ、12……表示器、13……コントローラ、1
4……勾配磁場制御回路、15……静磁場制御回
路、16……高周波発振器、17……ゲート変調
回路、18……電力増幅器、19……プリアン
プ、20……位相検波回路、21……波形メモ
リ、30……操作卓。

Claims (1)

  1. 【特許請求の範囲】 1 対象物に高周波パルスおよび磁場を印加して
    核磁気共鳴信号を発生させ、この信号を用いて対
    象物の組織に関する画像を得るようにした核磁気
    共鳴撮像装置において、 撮像したい部分のT1、T2、ρの範囲で、撮像
    時間が短くしかも高S/N比、高精度なT1、
    T2、ρ計算画像が得られるようなパルス系列と
    下記ないしの方法により求めたスキヤンパラ
    メータにより、同一の撮像部位についてはいかな
    る疾患を撮像する目的に対しても同一スキヤンで
    撮像し、 この撮影により得られたデータからT1、T2、
    ρ計算画像を求め、操作卓からの入力信号に応じ
    て画像を再合成する機能を有する制御・演算手段
    を具備したことを特徴とする核磁気共鳴撮像装
    置。 記 信号強度の理論式を定める。 この理論式と、測定したいT1、T2、ρの範
    囲と、原画像の分散から、計算画像の評価関数
    をスキヤンパラメータの多変数関数として求め
    る。 多変数関数の極値を求める方法により前記評
    価関数が最良となるスキヤンパラメータを求め
    る。 2 前記スキヤンの際に用いるパルス系列とし
    て、IRnSE法とSRmSE法(n、mは整数)、ま
    たはIRnSE法とFRmSE法(n、mは整数)を使
    用するように構成したことを特徴とする請求項1
    に記載の核磁気共鳴撮像装置。 3 前記スキヤンの際に用いるパルス系列とし
    て、IR3SE法とSR4SE法、またはIR3SE法と
    FR4SE法を使用するように構成したことを特徴
    とする請求項1に記載の核磁気共鳴撮像装置。 4 前記操作卓は、画像再合成の際、パルス系列
    の選択と、スキヤンパラメータ設定の連続可変が
    可能であり、しかも再合成画像表示がリアルタイ
    ムに追従するような設定が可能に構成されたこと
    を特徴とする請求項1に記載の核磁気共鳴撮像装
    置。
JP60196661A 1985-02-27 1985-09-05 核磁気共鳴撮像装置 Granted JPS6264348A (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP60196661A JPS6264348A (ja) 1985-09-05 1985-09-05 核磁気共鳴撮像装置
US06/833,224 US4694250A (en) 1985-02-27 1986-02-24 Nuclear magnetic resonance imaging device
DE19863606220 DE3606220A1 (de) 1985-02-27 1986-02-26 Nmr-abbildungsvorrichtung
GB08604693A GB2173598B (en) 1985-02-27 1986-02-26 Nuclear magnetic resonance imaging device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP60196661A JPS6264348A (ja) 1985-09-05 1985-09-05 核磁気共鳴撮像装置

Publications (2)

Publication Number Publication Date
JPS6264348A JPS6264348A (ja) 1987-03-23
JPH0374101B2 true JPH0374101B2 (ja) 1991-11-25

Family

ID=16361491

Family Applications (1)

Application Number Title Priority Date Filing Date
JP60196661A Granted JPS6264348A (ja) 1985-02-27 1985-09-05 核磁気共鳴撮像装置

Country Status (1)

Country Link
JP (1) JPS6264348A (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4709212A (en) * 1986-01-03 1987-11-24 General Electric Company Method of enhancing image signal-to-noise ratio by combining NMR images of differing pulse sequence timing
JPH04241839A (ja) * 1991-01-08 1992-08-28 Fuji Electric Co Ltd Mri画像処理方法
JP4316078B2 (ja) * 1999-11-29 2009-08-19 株式会社東芝 Mri装置
JP4768907B2 (ja) * 2000-09-22 2011-09-07 株式会社東芝 磁気共鳴イメージング装置および記憶媒体
JP6687383B2 (ja) 2015-12-21 2020-04-22 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59116534A (ja) * 1982-12-24 1984-07-05 Hitachi Ltd Nmrイメ−ジング装置
JPS61249457A (ja) * 1985-04-26 1986-11-06 デユ−ク・ユニバ−シテイ− 核磁気共鳴画像を自動合成する装置
JPS61276073A (ja) * 1985-05-31 1986-12-06 Shimadzu Corp Nmr−ct装置の画像作成法

Also Published As

Publication number Publication date
JPS6264348A (ja) 1987-03-23

Similar Documents

Publication Publication Date Title
US6171241B1 (en) Method for measuring myocardial motion and the like
US5604435A (en) Spiral scanning method for monitoring physiological changes
US6064203A (en) Method and apparatus for determining or imaging longitudinal spin relaxation time or producing images which substantially reflect longitudinal spin relaxation time contrast
US6546275B2 (en) Determination of the arterial input function in dynamic contrast-enhanced MRI
US6806709B2 (en) Flow imaging using balanced phase contrast steady state free precession magnetic resonance imaging
JP5366484B2 (ja) 磁気共鳴イメージング装置およびこの磁気共鳴イメージング装置における脂肪抑制効果の分析方法
US7418286B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP2007117765A (ja) Epiシーケンスにおける傾斜誘起された交差項磁場の測定および補正
CN101354433A (zh) 记录和处理时间上相继的图像数据组的序列的方法及设备
US7609059B2 (en) Magnetic resonance method and apparatus with nuclear spins type-specific signal suppression
US6198960B1 (en) Flip angle modulated magnetic resonance angiography
US20110237931A1 (en) Magnetic resonance imaging apparatus and synchronous imaging method
US8217649B2 (en) System for use in MR imaging using tissue mechanical resonance
JP2011200342A (ja) 磁気共鳴イメージング装置
US6759848B2 (en) Magnetic resonance method and apparatus for spatially resolved presentation of change in the functional activities of a brain
JPH0374101B2 (ja)
JP2005111059A (ja) 磁気共鳴イメージング装置
JPH0549610A (ja) 磁気共鳴診断装置
US9915713B2 (en) Determining positions of a magnetic field probe in a magnetic resonance measurement
JP3903276B2 (ja) 磁気共鳴撮像装置
JP3907944B2 (ja) 磁気共鳴イメージング方法及び装置
JP2003325477A (ja) 磁気共鳴イメージング装置
JP5116257B2 (ja) 磁気共鳴撮影装置
CN112824922B (zh) 磁共振拍摄装置及其控制方法
JP2004148024A (ja) N―アセチルアスパレート、グルタミンおよびグルタメートの定量化方法並びに磁気共鳴撮像装置

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees