磁気共鳴イメージング装置
【課題】ディフェーズまたはリフェーズにより収集された磁気共鳴信号からでは得られなかった有益な情報を得ることを可能とする。
【解決手段】収集手段、再構成手段、定量化手段および生成手段を備える。収集手段は、被検体から放射される磁気共鳴信号を収集する。再構成手段は、収集手段により収集された磁気共鳴信号に基づいてディフェーズ画像およびリフェーズ画像を少なくとも1枚ずつ再構成する。定量化手段は、再構成手段により再構成されたディフェーズ画像およびリフェーズ画像の双方に基づいて前記被検体に関する特性を定量化する。生成手段は、定量化手段により定量化された特性を表す複数種の定量化画像を生成する。
【解決手段】収集手段、再構成手段、定量化手段および生成手段を備える。収集手段は、被検体から放射される磁気共鳴信号を収集する。再構成手段は、収集手段により収集された磁気共鳴信号に基づいてディフェーズ画像およびリフェーズ画像を少なくとも1枚ずつ再構成する。定量化手段は、再構成手段により再構成されたディフェーズ画像およびリフェーズ画像の双方に基づいて前記被検体に関する特性を定量化する。生成手段は、定量化手段により定量化された特性を表す複数種の定量化画像を生成する。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、被検体から放射される磁気共鳴信号に基づいて医用診断のための情報を得る磁気共鳴イメージング装置に関する。
【背景技術】
【0002】
磁気共鳴信号の収集には、スピンをディフェーズさせて行う方法と、スピンをリフェーズさせて行う方法とがある。これら2つの収集方法は、用途に応じて選択的に使用される。そして、いずれかの方法で収集された磁気共鳴信号に基づいて、医用診断のための情報が得られる。
【0003】
また、いずれかの方法で収集された磁気共鳴信号に基づいて得られた情報に対して、別の処理を施すことによってより有用な情報を得る技術が知られている。例えば、T2*強調よりさらに磁化率変化に鋭敏な撮像法として、リフェーズで収集された磁気共鳴信号に基づいて再構成された絶対値画像に対して位相強調処理を施す手法が非特許文献1で提案されている。
【先行技術文献】
【非特許文献】
【0004】
【非特許文献1】Magn Reson Med 52:612-618, 2004.
【発明の概要】
【発明が解決しようとする課題】
【0005】
しかしながら、ディフェーズまたはリフェーズにより収集された磁気共鳴信号からでは得られる情報に限りがあり、例えば血管などを良好に描出するのに必要な情報を得ることができなかった。
【0006】
本発明はこのような事情を考慮してなされたものであり、その目的とするところは、ディフェーズまたはリフェーズにより各々単独で収集された磁気共鳴信号からでは得られなかった有益な情報を得ることが可能な磁気共鳴イメージング装置を提供することにある。
【課題を解決するための手段】
【0007】
以上の目的を達成するために本発明は、被検体から放射される磁気共鳴信号を収集する収集手段と、収集された前記磁気共鳴信号に基づいてディフェーズ画像およびリフェーズ画像を少なくとも1枚ずつ再構成する再構成手段と、再構成された前記ディフェーズ画像および前記リフェーズ画像の双方に基づいて前記被検体に関する特性を定量化する定量化手段と、前記定量化された特性を表す複数種の定量化画像を生成する生成手段とを磁気共鳴イメージング装置に備えた。
【図面の簡単な説明】
【0008】
【図1】本発明の一実施形態に係る磁気共鳴イメージング装置の構成を示す図。
【図2】図1に示す磁気共鳴イメージング装置における撮像の際の処理手順を示す図。
【図3】リフェーズTEとディフェーズTEとを1点ずつ収集する場合におけるサンプル点の関係を示す図。
【図4】ディフェーズTEを2点収集する場合におけるサンプル点の関係を示す図。
【図5】リフェーズTEを2点収集する場合におけるサンプル点の関係を示す図。
【図6】GREマルチエコーシーケンスの一例を示す図。
【図7】4点法におけるサンプル点の関係を示す図。
【図8】解析パラメータ画像の算出処理の概念を示す図。
【図9】リフェーズTEを2点とディフェーズTEを1点収集する場合における解析パラメータ画像の算出処理の流れを示す図。
【図10】図9中のステップSb2において作成する振幅画像に相当するサンプル点と実際の収集に係るサンプル点との関係を示す図。
【図11】リフェーズTEを1点とディフェーズTEを2点収集する場合における解析パラメータ画像の算出処理の流れを示す図。
【図12】図11中のステップSc2において作成する振幅画像に相当するサンプル点と実際の収集に係るサンプル点との関係を示す図。
【図13】非対称型のマスクを表す図。
【図14】対称型のマスクを表す図。
【図15】カラー合成処理の具体例を示す図。
【図16】kに対するTR,TEの変化の例を示す図。
【発明を実施するための形態】
【0009】
以下、図面を参照して本発明の実施形態について説明する。
【0010】
図1は本実施形態に係る磁気共鳴イメージング装置(以下、MRI装置と称する)の構成を示す図である。
【0011】
このMRI装置は、被検体としての患者Pを載せる寝台部と、静磁場を発生させる静磁場発生部と、静磁場に位置情報を付加するための傾斜磁場発生部と、高周波信号を送受信する送受信部と、システム全体のコントロールおよび画像再構成を担う制御・演算部と、被検体Pの心時相を表す信号としてのECG信号を計測する心電計測部と、患者Pに息止めを指令するための息止め指令部とを備えている。
【0012】
静磁場発生部は、例えば超電導方式の磁石1と、この磁石1に電流を供給する静磁場電源2とを含み、被検体Pが遊挿される円筒状の開口部(診断用空間)の軸方向(Z軸方向)に静磁場H0を発生させる。なお、この磁石部にはシムコイル14が設けられている。このシムコイル14には、後述するホスト計算機の制御下で、シムコイル電源15から静磁場均一化のための電流が供給される。寝台部は、被検体Pを載せた天板を磁石1の開口部に退避可能に挿入できる。
【0013】
傾斜磁場発生部は、傾斜磁場コイルユニット3を含む。この傾斜磁場コイルユニット3は、互いに直交するX、YおよびZ軸方向の傾斜磁場を発生させるための3組のコイル3x、コイル3y,コイル3zを備える。傾斜磁場発生部はまた、コイル3x〜3zに電流を供給する傾斜磁場電源4を含む。傾斜磁場電源4は、後述するシーケンサ5の制御のもと、傾斜磁場を発生させるためのパルス電流をコイル3x〜3zに供給する。
【0014】
傾斜磁場電源4からコイル3x〜3zに供給されるパルス電流を調整することにより、物理軸であるX,Y,Z方向の各軸の傾斜磁場を合成して、互いに直交するスライス方向傾斜磁場Gs、位相エンコード方向傾斜磁場Ge、および読出し方向(周波数エンコード方向)傾斜磁場Grの各論理軸方向を任意に設定することができる。スライス方向、位相エンコード方向および読出し方向の各傾斜磁場は、静磁場H0に重畳される。
【0015】
送受信部は、磁石1内の撮影空間にて被検体Pの近傍に配設されるRFコイル7と、このコイル7に接続された送信器8Tおよび受信器8Rとを含む。送信器8Tおよび受信器8Rは、後述するシーケンサ5の制御の下で動作する。送信器8Tは、核磁気共鳴(NMR)を起こさせるためのラーモア周波数のRFパルスをRFコイル7に供給する。受信器8Rは、RFコイル7が受信したエコー信号(高周波信号)を取り込み、これに前置増幅、中間周波変換、位相検波、低周波増幅、フィルタリングなどの各種の信号処理を施した後、A/D変換してエコー信号に応じたデジタル量のエコーデータ(原データ)を生成する。
【0016】
制御・演算部は、シーケンサ(シーケンスコントローラとも呼ばれる)5、ホスト計算機6、演算ユニット10、記憶ユニット11、表示器12、入力器13および音声発生器16を含む。このうち、ホスト計算機6は、記憶したソフトウエア手順により、シーケンサ5にパルスシーケンス情報を指令するとともに、装置全体の動作を統括する機能を有する。
【0017】
ホスト計算機6は、位置決め用スキャンなどの準備作業に引き続いて、イメージングスキャンを実施する。イメージングスキャンは、画像再構成に必要なエコーデータの組を収集するスキャンであり、ここでは2次元スキャンに設定されている。イメージングスキャンは、ECG信号に依るECGゲート法を併用して行うことができる。なお、このECGゲート法は場合によっては併用しなくても良い。
【0018】
シーケンサ5は、CPUおよびメモリを備えており、ホスト計算機6から送られてきたパルスシーケンス情報を記憶し、この情報にしたがって傾斜磁場電源4、送信器8Tおよび受信器8Rのそれぞれの動作を制御するとともに、受信器8Rが出力したエコーデータを一旦入力し、これを演算ユニット10に転送する。パルスシーケンス情報は、一連のパルスシーケンスにしたがって傾斜磁場電源4、送信器8Tおよび受信器8Rを動作させるために必要な全ての情報であり、例えばコイル3x〜3zに印加するパルス電流の強度、印加時間、印加タイミングなどに関する情報を含む。
【0019】
演算ユニット10は、受信器8Rが出力したエコーデータをシーケンサ5を介して入力する。演算ユニット10は、その内部メモリ上のフーリエ空間(k空間または周波数空間とも呼ばれる)にエコーデータを配置し、このエコーデータを各組毎に2次元または3次元のフーリエ変換に付して実空間の画像データに再構成する。また演算ユニットは、再構成された画像に基づく解析パラメータ画像の作成や合成処理などの画像処理を行う機能を有する。
【0020】
記憶ユニット11は、再構成された画像データのみならず、上述の合成処理や差分処理が施された画像データを保管することができる。表示器12は、ホスト計算機6の制御の下に画像を表示する。入力器13を介して、術者が希望する撮影条件、パルスシーケンス、画像合成や差分演算に関する情報をホスト計算機6に入力できる。
【0021】
息止め指令部の一要素として音声発生器16を備える。この音声発生器16は、ホスト計算機6から指令の下に、息止め開始および息止め終了のメッセージを音声として発することができる。
【0022】
心電計測部は、被検体の体表に付着させてECG信号を電気信号として検出するECGセンサ17と、このセンサ信号にデジタル化処理を含む各種の処理を施してホスト計算機6およびシーケンサ5に出力するECGユニット18とを含む。心電計測部による計測信号は、イメージングスキャンを実行するときにシーケンサ5により用いられる。これにより、ECGゲート法(心電同期法)による同期タイミングを適切に設定でき、この同期タイミングに基づくECGゲート法のイメージングスキャンを行ってデータ収集できるようになっている。
【0023】
次に以上のように構成されたMRI装置の動作について詳細に説明する。
【0024】
(A)理論的準備
具体的な動作について説明する前の理論的準備として、不均一磁化率と流れのあるボクセルのMR信号モデルと、リフェーズ(Rephase)/ディフェーズ(Dephase)の信号モデルとについて定義する。
【0025】
(A-1)不均一磁化率と流れのあるボクセルのMR信号モデル
まず、不均一磁化率と流れのあるボクセルのMR信号モデルに関するパラメータを次のように定義する。
M0:プロトン密度
AT1:T1に依存した減衰(T1 dependent attenuation)
AT1=1−exp(−TR/T1)
AT2:T2に依存した減衰(T2 dependent attenuation)
AT2=exp(−TE/T2)
AD:拡散に依存した減衰(Diffusion dependent attenuation)
AD=exp[−bD]
Asus:磁化率に依存した減衰(Susceptibility dependent attenuation)
Asus=exp[−TE(γΔB0σ)]*
φsus:磁化率に依存した位相(Susceptibility dependent phase)
φsus=−TE(γΔB0m)
但しLorenzian modelの場合。
T2*:T2と磁化率効果による成分を含む緩和時間
1/T2*=1/T2+γΔB0 (T2*<T2)
なお、T2*を用いると、AT2 Asus=exp(−TE/T2*)となる。
Aflow:フローに依存した減衰(Flow dependent attenuation)
Aflow=exp[−bDflow]
b:拡散で定義される傾斜磁場パターンから算出される係数(gradient factor)
Dflow:フローによる位相分散係数(流れのない部分で定義される拡散係数に等価な係数に相当する)
Φflow:フローに依存した位相シフト(Flow dependent phase shift)
【数1】
【0026】
V:不均一磁化率かつ流れを含む組織。ΔB0<>0であり、かつF<>0である組織であり、例えば静脈などがこれに相当する。
【0027】
この場合、組織VからのMR信号Sは、次の式(1)のように一般化される。
S=(M0 AT1 AT2 AD Asus Aflow) exp[i(Φ0+Φsus+Φflow)] …(1)
ここで磁化率やフロー(flow)の効果のうち、コヒーレント(coherent)な成分は位相に、インコヒーレント(incoherent)な成分は振幅減衰に生じる。すなわち、コントラストには、ボクセル内でのΔB0分布が大きいほど、あるいはIVIM(intravoxel incoherent motion)成分が多い(例えば、静脈→細静脈→毛細血管の順に大きい)ほど、振幅項が位相項より支配的に作用することになる。
【0028】
(A-2)リフェーズ(Rephase)/ディフェーズ(Dephase)の信号モデル
次に、リフェーズ/ディフェーズの信号モデルについて説明する。
【0029】
理想的なモデルでは、リフェーズでは振幅・位相ともフローによる成分はキャンセルされ、磁化率による成分のみとなる。ただし、実際はリフェーズでもGMN(gradient moment nulling)オーダーや乱流に依存した変化が生じる。一方、ディフェーズでは、振幅・位相とも磁化率成分にフロー成分が加わるので、リフェーズ/ディフェーズの各々の理想的なモデルでの信号Sre,Sdeは下記の式(2)および式(3)のように表される。
Sre=(M0 AT1 AT2 AD Asus)exp[i(Φ0+Φsus)] …(2)
Sde=(M0 AT1 AT2 AD Asus Aflow)exp[i(Φ0+Φsus+Φflow)] …(3)
ここでもし、リフェーズ/ディフェーズのシーケンス条件のうちの繰り返し時間TRおよびエコー時間TEが同一で、ディフェーズのb値は十分小さくてADの効果は無視できるとすれば、Sdeは次の式(4)となる。
Sde=Sre Aflow exp[i Φflow] …(4)
すなわち、SdeはSreにフロー効果が加わったものになる。
【0030】
MR信号(振幅・位相)の変動には次のような性質がある。
・ボクセル内のフロー成分比Mflow/(Mflow+Mst)が多いほど大きい。
・b値が大きいほど大きい。
・IVIM成分が大きければ、dephase grad.による振幅減衰効果が大きい。
・IVCM成分が大きければ、dephase grad.による位相変化効果が大きい。
(B)ディフェーズ/リフェーズのデータ収集および処理法
従来のリフェーズのみやディフェーズのみの場合に対し、ディフェーズとリフェーズを組み合わせて用いることで新たな応用が可能になる。上述のように、リフェーズはT1、T2緩和と磁化率効果のみによる信号変化であるが、ディフェーズはリフェーズに加えフローによる効果が加わっている。従って、リフェーズとディフェーズとのデータ間の演算により、静止部分と磁化率効果などフロー効果以外の成分とフロー効果が分離された画像、さらには磁化率やフロー効果のより定量的な普遍的なパラメータ画像を得ることができる。臨床的には、ディフェーズおよびリフェーズの物理的なフロー効果および磁化率効果の大小を用いた、静脈と動脈の分離などである。ここでは収集の工夫やその臨床応用例をあげる。
【0031】
(B-1)理論
具体的な処理手順を説明するのに先立って、理論について説明する。
【0032】
リフェーズおよびディフェーズのMR信号は、振幅および位相を有する複素信号であり、MR信号のゲインをK、A0=M0、AT1、AT2として、式(2),(3)を任意のゲインKとシーケンスの可変パラメータであるTR,TE,b値の関数として、以下の式(2')および式(3')のように表せる。
Sre(TR,TE)=Are(TR,TE) exp[iΦre(TE)]=K A0(TR,TE) Asus(TE) exp[i{Φ0(TE)+Φsus(TE)}] …(2')
Sde(TE,TE,b)=Ade(TE,TE) exp[iΦde(TE)]=K A0(TE,TE) Asus(TE) Aflow(b) exp[i{Φ0(TE)+Φsus(TE)+Φflow(b)}] …(3')
なお、ここではb値はフローによる信号減衰効果をもたらせば良いため十分小さく、A0に占める分子拡散の効果ADは無視できる。すなわちA0は、リフェーズ,ディフェーズで同一とした。また縦磁化Mzは、M0 AT1として求まり、AT1はGRE(gradient echo)の場合にはインフロー効果も含めTRで決まる。
【0033】
(B-1-1)生体組織と磁化率、フロー
磁化率とフローの性質の違いを生体内の組織別にみると、動脈はオキシヘモグロビン(oxyHb)を多く含むため磁化率効果が無視できるので、リフェーズでは動脈は描出されない。一方、静脈は動脈よりはフロー効果は小さく、さらにデオキシヘモグロビン(deoxyHb)による磁化率効果が加わる。従って、血管部については、リフェーズではフローによる位相変化が完全にリフォーカスされれば(集まれば)、静脈信号(磁化率効果成分)が支配的に描出され、ディフェーズでは動脈および静脈(磁化率効果+フロー成分)のいずれもが描出される。静止部については、リフェーズおよびディフェーズのいずれでも同等に描出される。
【0034】
静脈の選択的描出にはフローが遅いため、磁化率効果以外に良い方法がない。また、リフェーズのみでは動脈と静止部との分離はコントラスト差が十分でないため困難である。
【0035】
(B-1-2)最適シーケンスパラメータ:TE
TEについては、振幅、位相ともにCNR(contrast to noise ratio)を最大にするTEの最適な条件は、TE=T2*である。すなわち、対象とする組織のT2*に等しい場合に、その近傍のT2*を有する組織との間で最大のCNRを与える。収集条件や人体組織の種類に依存してT2*には分布があるが、ボクセルサイズを一定にし、血管描出に重要な血液も磁化率の範囲があるし、白質や灰白質はほぼ一定とみなせるのでほぼ決めることは可能である。
【0036】
(B-1-3)フロー効果と磁化率効果の分離
受信ゲイン、TR、TEおよびbなどの撮像パラメータを固定すれば、そのままでもリフェーズ・ディフェーズ間の相互の比較が可能であるが、より定量的な指標にするには、フロー成分はDflow、磁化率成分はT2*やΔχなどから計算して表示したほうが良い。その場合は、b=0はリフェーズで代用し、TEは2つ以上と最低3画像必要となる。
【0037】
(B-1-3-1)フロー効果の定量化
ディフェーズ/リフェーズ間の振幅比と位相差をとると、(4)式に基づいて次の式(5)および式(6)に示すようにフロー効果を分離できる。
ディフェーズ/リフェーズ振幅比:A(de)/A(re)=Aflow …(5)
ディフェーズ/リフェーズ位相差:Φ(de)−Φ(re)=Φflow …(6)
ここで、位相差ΦflowはPS-MRA(phase contrast MR angiography)法のように速度エンコードパルスVENCの3方向用いて測定し求めても良いが、ここではあまり現実的ではないので省略し、次の式(7)のように表される関係にあり、多様な流速や多様な方向への流れによるランダムな位相分散の程度を表すDflowを求める。
Aflow=exp[−b* Dflow] …(7)
ここで、ディフェーズのbファクタ(b-factor)をb(de)、ディフェーズおよびリフェーズの信号強度(振幅)をそれぞれA(de)、A(re)とすると、b=0をリフェーズで代用できるので、Dflowは次の式(8)により算出される。
Dflow[mm2/sec]=−ln[A(de)/A(re)]/b(de) …(8)
(B-1-3-2)磁化率効果の定量化
まず磁化率効果の定量化指標を振幅を用いて算出する場合について述べる。
【0038】
磁化率効果による減衰項Asusは、TE依存であり、単一TEではT1緩和およびT2緩和などの寄与が入ったA0が消去できないので、γΔB0σまたはT2を含んだT2*を求めたほうがより普遍化できる。
【0039】
リフェーズでの2つのTE(以下、それぞれをTE1,TE2と称し、TE2>TE1であることとする)を用いる場合、TE1およびTE2のそれぞれでの信号強度(振幅)A1,A2は、次の式(9)および式(10)により定まる。
A1=K A0 exp[−TE1/T2*] …(9)
A2=K A0 exp[−TE2/T2*] …(10)
これら式(9)および式(10)より、磁化率効果による振幅減衰効果の指標は次の式(11)で表される。
T2*=(TE2−TE1)/ln[A1(TE1)/A2(TE2)] …(11)
またT2*の算出には、3エコー以上のマルチエコーを取得して最小二乗近似を用いても良い。
【0040】
一方、位相を用いる場合について述べる。
【0041】
磁化率効果の別の表現としてリフェーズの単一TEまたは2エコーの位相差から低周波の位相成分を減じた位相からΦsusを求める。単一エコーの場合および2エコーの場合のΦsusは、それぞれ次の式(12)および式(12')により表される。
Φsus=−2πγΔχB0(cos2θ−1/3)TE …(12)
Φsus=−2πγΔχB0(cos2θ−1/3)(TE2−TE1) …(12')
これに基づいて、単一エコーの場合および2エコーの場合のΔχは、それぞれ次の式(13)および式(13')により表される。
Δχ[ppm]=−Φsus/{2πγ B0(cos2θ−1/3)TE} …(13)
Δχ[ppm]=−Φsus/{2πγ B0(cos2θ−1/3)(TE2−TE1)} …(13')
これはボクセル内の平均の磁化率をあらわすとみなせる。
【0042】
ここで、γは磁気回転比、B0は静磁場強度を表し、装置により一意に決まる。しかし、θはB0方向と血管の走行角度であるので、ボクセル間の相関を調べるなどの方法で測定する必要があり、振幅から求めるよりは困難になる。
【0043】
なお、磁化率分布のみならば、リフェーズのみでTEが2段階以上あれば実現可能であり、この技術に関しては公知である。
【0044】
次に、具体的な処理手順について図2に従って説明する。
【0045】
(B-2)ディフェーズデータおよびリフェーズデータの収集
ステップSa1においてはシーケンサ5が、ディフェーズデータおよびリフェーズデータを収集する。
【0046】
同一TEのディフェーズおよびリフェーズのシーケンスを独立させて連続的に撮像しても良い。あるいは動きの影響を最小にするために、1データセットを複数のセグメントに分割して交互に収集して合成しても良い。分割するセグメントは、k−spaceで1ライン(TR)単位としたり、2D面単位とするなど任意である。
【0047】
b値の段階は複数でも良い。TEも複数取得すればT2*の正確な算出にも役立つ。複数のTEは、Gradient Echo法ならTR内でマルチエコーを収集すれば一回で取得できる。また複数のb値もLook−Locker法を用いれば一回で取得できる。なおLook−Locker法は、「"Measurement of Gd-DTPA dialysis clearance rates by using a look-locker imaging technique.",Magn Reson,Med. 1996 Oct;36(4):571-8.」により知られている。
【0048】
(B-2-1)ディフェーズ・リフェーズ交互分割収集(ディフェーズ,リフェーズの一組は同一TE)
フロー効果Dflowのみの算出のためには、シーケンサ5はリフェーズおよびディフェーズを同一TEで2画像収集する。Dflowに加えてT2*を算出する場合には、シーケンサ5はリフェーズ側を2エコーでTE1、TE2とし、ディフェーズはTE1,TE2のいずれかにする。その場合、ディフェーズおよびリフェーズのシーケンスを分割させて別々に撮像しても良い。あるいは、画像間の動きの影響を最小にするために、1データセットを複数のセグメントに分割して交互に収集して合成しても良い。分割するセグメントは、k−spaceで1ライン(TR)単位としたり、2D面単位とするなど任意である。
【0049】
(B-2-2)GREマルチエコーによるリフェーズ,ディフェーズ混合連続収集の場合(TEがすべて違う)
GREでリフェーズおよびディフェーズ組み合わせた2点以上のマルチエコーで処理する実施例を示す。
【0050】
上記の(B-2-1)のようにリフェーズおよびディフェーズで同一TEになるように別々のデータとして取得する場合に比べ、マルチエコーで収集すれば、リフェーズおよびディフェーズのTEは同一にできないが、複数のTEで収集し計算により同一にでき、定量的なパラメータが得られる。マルチエコーは1回のRF励起で、すなわち同一のTR内で取得できるため、スキャン時間が1エコーの収集時間と同程度であるという最大のメリットがあるのに加え、独立に時間をおいて収集する場合に比べてデータ間の体動も無視できるため、異種データ間の演算時の誤差が軽減できるメリットもある。もちろん別々のデータとして取得しても、後述する解析パラメータ画像の算出処理は共通にできる。
【0051】
(B-2-2-1)2点法
(B-2-2-1-1)リフェーズTEとディフェーズTEとを1点ずつ収集する場合
図3に示すように、TE=TE1のAde(TE1)と、TE=TE2のAre(TE2)とをそれぞれ収集する。
【0052】
この場合は、同一モードのデータが2点以上ないため、2つの画像をそのまま観察するのが主体で、定量パラメータのうちT2*,Dflowなどは算出できない。しかし、条件をTE2=2*TE1に設定すれば、流れの位相項Φflowのみは定量化可能である。
【0053】
TE(de)<TE(re)とすることで、ディフェーズは磁化率効果を抑制したフロー強調になり、リフェーズはフロー効果を抑制した磁化率強調になる。TE(re)<TE(de)といった具合にTE(re)を短く(<10ms)設定すれば、リフェーズは動脈描出用として通常ルーチンで収集しているTOF−MRA(time of flight-magnetic resonance angiography)の代用となり、通常のTOF−MRAが不要になる可能性がある。さらに、ディフェーズは、フローと磁化率効果も強調した動静脈描出用にできるし、TEの設定しだいで静脈の描出は制御できる。また位相情報も用いれば背景組織に対する静脈のCNRを強調できる。さらに、後述するようにディフェーズとリフェーズとの間で相互に演算処理を施すことにより、動脈と静脈の分離表示も可能である。なお、2エコーの場合でもリフェーズおよびディフェーズのTEは厳密に合わせる必要はなく、十分近くT2*の差が大きくなければ半定量的な表示は可能である。
【0054】
(B-2-2-1-2)ディフェーズTEを2点収集する場合
図4に示すように、TE=TE1に関するAde(TE1)およびTE=TE2に関するAde(TE2)をそれぞれ収集する。
【0055】
De TE=TE1,TE2 (TE1<TE2)の2つのエコーのb値を同じにすれば、Aflow=exp[−bDflow]、またGradient momentも同じにすれば位相のフロー効果分Φflowもほぼ同じであるから、次の式(14)からフロー効果がキャンセルされることが明らかである。
【0056】
Sde(TE2,b)/Sde(TE1,b)=exp[−(TE2−TE1)/T2*] exp[−i(TE2−TE1)γΔB0m] …(14)
そしてTE1,TE2は既知なので、式(13)の振幅からT2*、さらには位相項からΔB0をそれぞれ算出可能であり、Δχが算出可能である。
【0057】
(B-2-2-1-3)リフェーズTEを2点収集する場合
図5に示すように、TE=TE1に関するAre(TE1)およびTE=TE2に関するAre(TE2)をそれぞれ収集する。
【0058】
A2でb=0とみなせるので、Aflow=1 Φflow=0となる以外は(B-2-2-1-2)における添え字deをreに置き換えたものになる。
【0059】
(B-2-2-2)3点法
リフェーズおよびディフェーズをまぜて3点収集すれば、T2*またはΔB0に加えて、すなわちΔχに加えて、フロー効果分も算出可能であるし、またGREのマルチエコーを用いるにもかかわらずに同一TEのリフェーズ画像およびディフェーズ画像が作成可能である。3点のTEは、どう組み合わせても良い。
【0060】
3点法では各エコーがリフェーズおよびディフェーズのいずれであるかの組み合わせで2×2×2=8通り、順番を問題にしなければ4通りがある。このうちリフェーズ,リフェーズ,ディフェーズと、ディフェーズ,ディフェーズ,リフェーズとの2通りを取り上げる。
【0061】
(B-2-2-2-1)リフェーズTEを2点とディフェーズTEを1点収集する場合
最初のリフェーズの2点よりT2*を算出し、これを3点目のディフェーズに代入してDflowを算出する。またリフェーズよりT2*が分かれば、リフェーズの任意のTEの信号を生成可能なので、ディフェーズと同一TEのリフェーズが求まる。すなわち同一のT2緩和効果および磁化率効果でフロー効果のみ異なる画像が得られる。
【0062】
このときのGREマルチエコーシーケンスの一例を図6に示す。
【0063】
(B-2-2-2-2)ディフェーズTEを2点とリフェーズTEを1点収集する場合
上記の(B-2-2-1-2)と同様にして、最初のb値が同じでTEが異なるディフェーズの2点よりT2*を算出し、3点目のリフェーズとからDflowを算出する。またディフェーズよりT2*が分かればディフェーズの任意のTEの信号を生成可能なので、リフェーズと同一TEのディフェーズが求まる。すなわち同一のT2緩和効果および磁化率効果でフロー効果のみ異なる画像が得られる。
【0064】
(B-2-2-3)4点以上法
ディフェーズ、リフェーズとも各々2点以上ずつ収集し、T2*,同一TEの画像を求める。未知数2個で4点以上では最小二乗近似となる。もちろん必要ならば未知パラメータのM0、T2やDなども算出可能である。
【0065】
図7は4点法におけるデータ収集の一例を示す図であり、TE=TE1に関するAre(TE1)、TE=TE2に関するAde(TE2)、TE=TE3に関するAre(TE3)、TE=TE4に関するAde(TE4)をそれぞれ収集する。
【0066】
なお、以上に説明した各種のシーケンスモードのいずれにおいても、マルチエコーのTEを、T2*算出に用いる最適TE(=T2*)を含むように設定する。対象のT2*が長すぎたり(時間延長、SNR低下)、短かすぎたり(RFが入らない、傾斜磁場が出ない、十分なb値が出ないなど)で最適TEの設定困難な場合は、計算して作る任意のTEの画像からTE=T2*の画像を作成すれば良い。
【0067】
また、位相算出においては、すべての場合で折り返しがないか、折り返しを補正することが好ましい。GREのマルチエコーでもリフェーズおよびディフェーズのTEは厳密にあわせる必要はなく、十分近いTEでT2*の差が大きくなければ半定量的な表示は可能である。シーケンスもGREのみならず、1回のRF励起後に実効TEを数段階変えたmulti−shot EPI(echo planar imaging)でk−spaceをセグメント分割するなどして短時間化と高分解能化を適当にコントロールしても良い。
【0068】
またTEに関しては、水と脂肪とが混在する組織では、水脂肪が同位相になるように設定することが重要である。静磁場強度のもとでの水(プロトン)とケミカルシフトδppmの物質との位相差は、ΔΦ=2πγδB0TEと表されるが、それが同位相となる条件はnを整数としてδΦ=n2πなので、TE=n/(γδB0)の倍数に設定すれば良い。脳実質では脂肪があまりないためにあまり問題とはならないが、ボクセル内で水と脂肪とが混在する骨髄や腹部臓器では問題となる場合があるのでその条件が必要になる。γ=42.6MHz/T、脂肪はγ=3.6ppm、B0=1.5TとするとTE=n 4.3 msとなる。さらに脂肪以外でも磁化率の異なる物質がボクセル内に混在すると位相差が生じるが、酸素濃度によるケミカルシフトはδ=0.1ppm程度とされ、TEを適当に選んでも定量化パラメータが振幅のT2*ではあまり問題とならないし、位相では厳密に定量化する場合は問題となり得るので、必要ならδを既知としてその寄与による位相を求め補正すれば良い。
【0069】
(B-3)画像再構成
ステップSa2においては演算ユニット10が、以上のような各種の手法を適宜に採用して収集されたデータのそれぞれを使用して周知の再構成処理を行うことによって、m枚のディフェーズ画像とn枚のリフェーズ画像とを再構成する。なお、mおよびnの値は、いずれも0を含む整数で、採用するデータ収集の手法により定まる。
【0070】
(B-4)解析パラメータ画像の算出
ステップSa3においては演算ユニット10が、再構成された1乃至m枚のディフェーズ画像および1乃至n枚のリフェーズ画像を使用して解析パラメータ画像を算出する。
【0071】
図8は解析パラメータ画像の算出処理の概念を示す図である。
【0072】
演算ユニット10は、リフェーズ画像およびディフェーズ画像を使用してフローパラメータ算出処理P1を行うことによってよってDflow画像を算出する。演算ユニット10は、リフェーズ画像およびディフェーズ画像を使用して磁化率パラメータ算出処理P2を行うことによってT2*画像およびΔχ画像を算出する。演算ユニット10は、リフェーズ画像、ディフェーズ画像およびT2*画像を使用して任意TE画像作成処理P3を行うことによって、任意TEに関するリフェーズ画像およびディフェーズ画像を算出する。
【0073】
以下に、上記した各種のデータ収集手法のそれぞれに応じた解析パラメータ画像の算出処理の具体例をそれぞれ説明する。
【0074】
(B-4-1)2点法
(B-4-1-1)ディフェーズ2点の場合
(a) 演算ユニット10は、同一b値のディフェーズに関するTE=TE1の振幅画像A1(TE1)およびTE=TE2の振幅画像A2(TE2)より、次の式(15)および式(16)によりT2*およびKdeをそれぞれ算出する。
T2*=(TE2−TE1)/ln[Ade(TE1)/Ade(TE2)] …(15)
Kde=Ade(TE1)/exp[−TE1/T2*] …(16)
(b) 演算ユニット10は、任意TEに関して、次の式(17)および式(18)によりディフェーズ振幅および位相をそれぞれ算出する。
Ade(TE)=Kde exp[−TE/T2*] …(17)
Φde(TE)=(TE/TE1)Φde(TE1) …(18)
(B-4-1-2)リフェーズ2点の場合
(a) 演算ユニット10は、同一b値のディフェーズに関するTE=TE1の振幅画像A1(TE1)およびTE=TE2の振幅画像A2(TE2)より、次の式(19)および式(20)によりT2*およびKreをそれぞれ算出する。
T2*=(TE2−TE1)/ln[Are(TE1)/Are(TE2)] …(19)
Kre=Are(TE1)/exp[−TE1/T2*] …(20)
(b) 演算ユニット10は、任意TEに関して、次の式(21)および式(22)によりディフェーズ振幅および位相をそれぞれ算出する。
Are(TE)=Kre exp[−TE/T2*] …(21)
Φre(TE)=(TE/TE1)Φre(TE1) …(22)
すなわち、(B-4-1)における添え字deをreに置き換えた式を適用する。
【0075】
(B-4-2)3点法
(B-4-2-1)リフェーズTEを2点とディフェーズTEを1点収集する場合
この場合の処理の流れを図9に示す。
【0076】
(a) ステップSb1において演算ユニット10は、リフェーズに関するTE=TE1の振幅画像Are(TE1)およびTE=TE2の振幅画像Are(TE2)より、次の式(23)および式(24)によりT2*およびKreをそれぞれ算出する。
T2*=(TE2−TE1)/ln[Are(TE1)/Are(TE2)] …(23)
Kre=Are(TE1)/exp[−TE1/T2*] …(24)
ここで求めたT2*により、T2*画像が得られる。
【0077】
(b) ステップSb2において演算ユニット10は、次の式(25)により、例えば図10に示すようなリフェーズに関するTE=TE3の振幅画像Are(TE3)を作成する。また、ステップSb3において演算ユニット10は、TEに依存しないバックグラウンド位相がない場合には式(26)により、TEに依存しないバックグラウンド位相がある場合には式(27)により、位相画像Φre(TE3)を作成する。
Are(TE3)=Kre exp[−TE3/T2*] …(25)
Φre(TE3)=(TE3/TE1)Φre(TE1) …(26)
Φre(TE3)={TE3/(TE2−TE1)} {Φre(TE2)−Φre(TE1)} …(27)
なお、ステップSb2およびステップSb3において演算ユニット10は、TEがTE1,TE2,TE3のいずれでもない任意の値であるリフェーズの振幅画像および位相画像を作成することもできる。これは、上記の式(25)乃至(27)におけるTE3に任意の欲しいTEを代入することにより計算することにより実現できる。
【0078】
(c) ステップSb4において演算ユニット10は、TE=TE3の振幅画像Are(TE3),Ade(TE3)より、次の式(28)および式(29)によりフロー分散係数Dflow画像を算出する。
Aflow=Ade(TE3)/Are(TE3) …(28)
Dflow[mm2/sec]=−ln[Aflow]/bde …(29)
(d) ステップSb5において演算ユニット10は、血管とB0との方向角度θを算出する。
【0079】
(e) さらにステップSb5において演算ユニット10は、θとTE=TE3のリフェーズの位相画像より、次の式(30)によりΔχ[ppm]を算出する。
Δχ[ppm]=−Φre(TE3)/{2πγB0(cos2θ−1/3)TE3} …(30)
ここで求めたΔχにより、Δχ画像が得られる。
【0080】
(f) ステップSb6において演算ユニット10は、任意TEのディフェーズの振幅画像を作成する。
ディフェーズの任意TEの振幅画像は式(28)により求めたAflowを用いて、次の式(31)により計算する。
Ade(TE)=Are(TE)Aflow …(31)
ステップSb7において演算ユニット10は、Φde(TE3)を用いて、次の式(32)によりディフェーズの任意TEの位相画像を算出する。これは、TEに依存しないバックグラウンド位相がない場合のみ算出可能である。
Φde(TE)=(TE/TE3)Φde(TE3) …(32)
なお、上記の(c)〜(f)で求める各種パラメータは、必要な場合にのみ算出すれば良い。
【0081】
(B-4-2-2)ディフェーズTEを2点とリフェーズTEを1点収集する場合
この場合の処理の流れを図11に示す。
【0082】
(a) ステップSc1において演算ユニット10は、ディフェーズに関するTE=TE1の振幅画像Ade(TE1)およびTE=TE2の振幅画像Ade(TE2)より、次の式(33)および式(34)によりT2*およびKdeを算出する。
T2*=(TE2−TE1)/ln[Ade(TE1)/Ade(TE2)] …(33)
Kde=Ade(TE1)/exp[−TE1/T2*] …(34)
ここで求めたT2*により、T2*画像が得られる。
【0083】
(b) ステップSc2において演算ユニット10は、次の式(35)により、図12に示すようにディフェーズに関するTE=TE3の振幅画像Ade(TE3)を作成する。また、ステップSc3において演算ユニット10は、次の式(36)により、位相画像Φde(TE3)を作成する。
Ade(TE3)=Kdeexp[−TE3/T2*] …(35)
Φde(TE3)=(TE3/TE1)Φde(TE1) …(36)
なお、ステップSc2およびステップSc3において演算ユニット10は、TEがTE1,TE2,TE3のいずれでもない任意の値であるディフェーズの振幅画像および位相画像を作成することもできる。これは、上記の式(35)および(36)におけるTE3に任意の欲しいTEを代入することにより計算することにより実現できる。
【0084】
(c) ステップSc4において演算ユニット10は、TE=TE3の振幅画像Are(TE3),Ade(TE3)より、次の式(38)および式(39)によりフロー分散係数Dflow画像を算出する。
Aflow=Ade(TE3)/Are(TE3) …(37)
Dflow[mm2/sec]=−ln[Aflow]/bde …(38)
(d) ステップSc5において演算ユニット10は、血管とB0との方向角度θを算出する。
【0085】
(e) さらにステップSc5において演算ユニット10は、θとTE=TE3の位相画像より、次の式(39)によりΔχ[ppm]を算出する。
Δχ[ppm]=Φre(TE3)/{2πγB0(cos2θ−1/3)TE3} …(39)
ここで求めたΔχにより、Δχ画像が得られる。
【0086】
(f) ステップSc6において演算ユニット10は、任意TEのリフェーズの振幅画像を作成する。
リフェーズの任意TEの振幅画像は式(37)により求めたAflowを用いて、次の式(40)により計算する。
Are(TE)=Ade(TE)/Aflow …(40)
ステップSc7において演算ユニット10は、Φre(TE3)を用いて、次の式(41)によりリフェーズの任意TEの位相画像を算出する。これは、TEに依存しないバックグラウンド位相がない場合のみ算出可能である。
Φre(TE)=(TE/TE3)Φre(TE3) …(41)
なお、(c)〜(f)で求めるパラメータは、必要な場合にのみ算出すれば良い。
【0087】
さらに、リフェーズおよびディフェーズの各々が3点以上の場合は、未知パラメータをモデルに即して最小二乗法で算出する。Kre,Kde,T2*は1次の指数減衰関数のモデルを用い、Δχは一次関数のモデルを用いれば良い。
【0088】
以上のようにして、少なくとも1つずつのリフェーズ画像およびディフェーズ画像(それぞれ振幅画像および位相画像を含む)、T2*画像、Δχ画像およびDflow画像などがそれぞれ算出される。なお以下においては、これらの画像を解析パラメータ画像と総称する。
【0089】
なお、マルチエコーではなく、1エコーずつ別々に収集して画像化したものを使用しても良い。あるいは、2Dマルチスライス収集や、3DFT法などによるボリューム収集でも良い。さらに、パルスシーケンスとしては、GREに代えて非対称スピンエコー(ASE)法を利用しても良い。
【0090】
(B-5)解析パラメータ画像の合成処理および表示
各種の解析パラメータ画像を算出した後、これらの解析パラメータ画像をそのまま表示して観察に供することができる。また、Aflow、Φflow、Δflowなどのパラメータ自体を表示しても良い。
【0091】
さらに診断の補助とするためにステップSa4において演算ユニット10は、各種の解析パラメータ画像のカラーで分離した2Dの合成画像を作成して、2D画像のまま表示することができる。このステップSa4の合成処理は、必要時だけに行うようにしても良い。
【0092】
あるいは、特に血管画像を目的とする場合にはステップSa5において演算ユニット10は、血管を連続した管として表現するために3D処理を行う。その代表的なものは、最大値投影(MIP)または最小値投影(minIP)である。演算ユニット10は、ステップSa5にて複数の3D画像を生成した上で、ステップSa6にてこれら複数の3D画像を用いてフュージョン処理を行うこともできる。
【0093】
(B-5-1)位相合成処理および表示
任意TEのリフェーズ画像およびディフェーズ画像より位相合成画像を作成する機能について説明する。
【0094】
(a) 通常のMR画像:S=A exp[iΘ]と、ローパスフィルタをかけたMR画像:S1=A1 exp[iΘ1]とをリフェーズおよびディフェーズのおのおのについて作成する。
【0095】
(b) 位相アーチファクト補正
次の式(42)によりリフェーズのMR画像の位相マップよりローパスフィルタをかけた画像の位相マップを減算する。
Φ=Θ-Θl=arg[S]−arg[S1] …(42)
この式(42)の計算後、位相飛び補正により−π<Φ<=πに収める処理を行う。すなわち、Φが−π以下であるならばΦ+πに、またΦがπよりも大きい場合にはΦ−πにΦを置き換える。
【0096】
なお、Φは、式(42)に代えて次の式(43)により算出しても良い。
Φ=arg[S/S1] …(43)
(c) 位相マスク(phase masking)
図13に示す非対称型(asymmetric type)または図14に示す対称型(symmetric type)のどちらかのマスクMを作成し、振幅画像とマスクMとのn重の積Iを次の式(44)により求める。
I=A×Mn …(44)
なお、Aはディフェーズの位相から、Mはリフェーズの位相から作成したものを用いる。
【0097】
これにより位相飛びやフローと磁化率由来の位相キャンセルのない位相マスクが可能になる。
【0098】
(B-5-2)カラー合成処理および表示
上記のように得られた解析パラメータ画像をカラーで色分けしてフュージョン表示する。例えば、フロー成分のAflow、ΦflowおよびDflowは動脈および静脈に支配的に大きくなり、磁化率成分のT2*およびΔχは静脈や、静止部なら出血などの部分で大きくなる。そこで図15に示す画像21のようなフロー成分の画像を赤、画像22のような磁化率成分の画像を青とするなどの色分けをしてカラー合成処理を行うことにより、例えば画像21としてフュージョン表示する。このときに表示されるのは、物理的にはフロー成分(赤)と磁化率成分(青)との合成画像である。血管部のみを表す画像の場合には、動脈は赤により、静脈は紫により表される。磁化率成分部は、出血や静止部の磁化率がゼロでない部分のアーチファクトも含まれるので、必要なら閾値処理などによる静脈抽出を行った後にフュージョン表示すれば良い。
【0099】
また、画像21および画像22から次の式(45)または式(46)によって動脈および静脈を抽出した画像24を作成して、この画像24を画像22とカラー合成することによって画像23を得ることもできる。
Aflow=Ade/Are …(45)
Φflow=Φde−Φre …(46)
色の混じり具合は、フロー成分と磁化率成分の比で決まる。上記の色の割り当て例では、赤→紫→青の変化を示すが、赤に近いほどフロー成分が、青に近いほど磁化率成分が多いことを意味する。脳梗塞などの疾患時は必ずしも各々動脈および静脈に対応せず、酸素代謝状態を反映する指標となる。色の割り当てはこれに限らず、対比できる組み合わせならなんでも良い。TOF−MRAで動脈を描出する方法もあるが、原理的に遅い動脈は描出されない。本法ではディフェーズによる信号減衰または位相変化を用いるため、動脈も細いのや、上方から回り込んでくる側副血行路も描出することができ、脳梗塞の診断などに向けた臨床的に重要な情報を提供できる。またある程度の長いTEに設定すれば、磁化率効果により血栓や出血などの情報も同時に得られるので、治療計画を立てる上でも重要な情報を提供できる。
【0100】
なお、画像21としては、ディフェーズ振幅画像、ディフェーズ位相画像、ディフェーズの振幅位相合成画像およびDflow画像などが利用できる。画像22としては、リフェーズ振幅画像、リフェーズ位相画像、リフェーズの振幅位相合成画像およびT2*画像などが利用できる。
【0101】
(B-5-3)3D処理および表示
ディフェーズの振幅画像では、動脈および静脈のいずれも周囲組織に対して低い画像値になるのでminIPが妥当である。リフェーズの振幅画像では、静脈は低信号なのでminIPが妥当だが、動脈はTOF(Time of Flight)効果で高信号になるためにMIPを併用しても良い。表面を抽出するなどしてボリュームレンダリングやサーフェスレンダリングも用いることができる。また原画像信号を見たい場合は、単純に断面変換(MPR)も用途により使用することが有効である。フュージョン表示は、前記のように2Dのままでも良いが、3D画像を作成してから行っても良い。
【0102】
なお、各種画像の表示は、表示器12により行うことが可能であるが、外部のビューワ装置などにおいて表示させても良い。
【0103】
以上のように本実施形態によれば、血流などのフローの効果を分離して定量化することができる。そしてこの定量化結果に基づいて、医用診断に有用な情報や画像を提供することが可能である。
【0104】
本実施形態によれば、磁化率やフローに関して、機種やシーケンスに依存しないパラメータが作成可能である。このため、共通のデータベースになり、重要な知見がエビデンスとして蓄積され易くなる。
【0105】
本実施形態によれば、リフェーズとディフェーズとを組み合わせることで動脈と静脈との分離が可能である。
【0106】
本実施形態によれば、リフェーズおよびディフェーズを別々に収集する方法に比べ、マルチエコーを用いることにより1回の収集で双方の画像が取得可能であるので、収集時間が短くなり、かつ複数画像間の時間差が無視できるのでより動きの影響が生じにくい。
【0107】
本実施形態によれば、2画像を用いるため、速い流速の血液のtime of fright (TOF)効果が抑制可能である。血液は、磁化率またはフローのみに依存したコントラストが得られる。
【0108】
本実施形態によれば、側副血行路の描出が可能である。
【0109】
本実施形態によれば、磁化率とフローとが分離可能である。
【0110】
本実施形態によれば、ディフェーズで得られた振幅画像では、血管内孔が選択的に描出されるという手術などでは際めて重要な特長を有する。
【0111】
この実施形態は、次のような種々の変形実施が可能である。
【0112】
位相マスク画像を作成する場合には、ディフェーズの振幅とリフェーズの位相とを、あるいはディフェーズの位相とリフェーズの振幅とを組み合わせてもよい。このようにすることで、位相の磁化率とフローとの相殺が生じないため、振幅の位相マスク画像における静脈信号がより低下し、動脈と静脈の分離能が向上する。
【0113】
撮像時間とSNR向上を目的としてk-spaceにおける周波数毎にTRまたはTEを可変しても良い。例えば、シーケンスによる収集時、低周波は短いTEで、中間から高周波は長いTEで収集する。その間は、TEを滑らかに変える。これにより、低周波成分が支配的な静磁場不均一による位相は小さくなるのでアーチファクトが低減する。位相の場合には縦磁化は無関係なので、TRも最短にして良い。なおシーケンス種はGRE(FE)タイプで、マルチエコーおよび1エコーのいずれでも適用可能で、k-space trajectoryは、spin warp、spiral、あるいはEPIなどのいずれでも良い。図16にkに対するTR,TEの変化の例を示す。
【0114】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【符号の説明】
【0115】
1…磁石、2…静磁場電源、3…傾斜磁場コイルユニット、3x,3y,3z…コイル、4…傾斜磁場電源、5…シーケンサ、6…ホスト計算機、7…RFコイル、8R…受信器、8T…送信器、10…演算ユニット、11…記憶ユニット、12…表示器、13…入力器、14…シムコイル、15…シムコイル電源、16…音声発生器、17…ECGセンサ、18…ECGユニット。
【技術分野】
【0001】
本発明は、被検体から放射される磁気共鳴信号に基づいて医用診断のための情報を得る磁気共鳴イメージング装置に関する。
【背景技術】
【0002】
磁気共鳴信号の収集には、スピンをディフェーズさせて行う方法と、スピンをリフェーズさせて行う方法とがある。これら2つの収集方法は、用途に応じて選択的に使用される。そして、いずれかの方法で収集された磁気共鳴信号に基づいて、医用診断のための情報が得られる。
【0003】
また、いずれかの方法で収集された磁気共鳴信号に基づいて得られた情報に対して、別の処理を施すことによってより有用な情報を得る技術が知られている。例えば、T2*強調よりさらに磁化率変化に鋭敏な撮像法として、リフェーズで収集された磁気共鳴信号に基づいて再構成された絶対値画像に対して位相強調処理を施す手法が非特許文献1で提案されている。
【先行技術文献】
【非特許文献】
【0004】
【非特許文献1】Magn Reson Med 52:612-618, 2004.
【発明の概要】
【発明が解決しようとする課題】
【0005】
しかしながら、ディフェーズまたはリフェーズにより収集された磁気共鳴信号からでは得られる情報に限りがあり、例えば血管などを良好に描出するのに必要な情報を得ることができなかった。
【0006】
本発明はこのような事情を考慮してなされたものであり、その目的とするところは、ディフェーズまたはリフェーズにより各々単独で収集された磁気共鳴信号からでは得られなかった有益な情報を得ることが可能な磁気共鳴イメージング装置を提供することにある。
【課題を解決するための手段】
【0007】
以上の目的を達成するために本発明は、被検体から放射される磁気共鳴信号を収集する収集手段と、収集された前記磁気共鳴信号に基づいてディフェーズ画像およびリフェーズ画像を少なくとも1枚ずつ再構成する再構成手段と、再構成された前記ディフェーズ画像および前記リフェーズ画像の双方に基づいて前記被検体に関する特性を定量化する定量化手段と、前記定量化された特性を表す複数種の定量化画像を生成する生成手段とを磁気共鳴イメージング装置に備えた。
【図面の簡単な説明】
【0008】
【図1】本発明の一実施形態に係る磁気共鳴イメージング装置の構成を示す図。
【図2】図1に示す磁気共鳴イメージング装置における撮像の際の処理手順を示す図。
【図3】リフェーズTEとディフェーズTEとを1点ずつ収集する場合におけるサンプル点の関係を示す図。
【図4】ディフェーズTEを2点収集する場合におけるサンプル点の関係を示す図。
【図5】リフェーズTEを2点収集する場合におけるサンプル点の関係を示す図。
【図6】GREマルチエコーシーケンスの一例を示す図。
【図7】4点法におけるサンプル点の関係を示す図。
【図8】解析パラメータ画像の算出処理の概念を示す図。
【図9】リフェーズTEを2点とディフェーズTEを1点収集する場合における解析パラメータ画像の算出処理の流れを示す図。
【図10】図9中のステップSb2において作成する振幅画像に相当するサンプル点と実際の収集に係るサンプル点との関係を示す図。
【図11】リフェーズTEを1点とディフェーズTEを2点収集する場合における解析パラメータ画像の算出処理の流れを示す図。
【図12】図11中のステップSc2において作成する振幅画像に相当するサンプル点と実際の収集に係るサンプル点との関係を示す図。
【図13】非対称型のマスクを表す図。
【図14】対称型のマスクを表す図。
【図15】カラー合成処理の具体例を示す図。
【図16】kに対するTR,TEの変化の例を示す図。
【発明を実施するための形態】
【0009】
以下、図面を参照して本発明の実施形態について説明する。
【0010】
図1は本実施形態に係る磁気共鳴イメージング装置(以下、MRI装置と称する)の構成を示す図である。
【0011】
このMRI装置は、被検体としての患者Pを載せる寝台部と、静磁場を発生させる静磁場発生部と、静磁場に位置情報を付加するための傾斜磁場発生部と、高周波信号を送受信する送受信部と、システム全体のコントロールおよび画像再構成を担う制御・演算部と、被検体Pの心時相を表す信号としてのECG信号を計測する心電計測部と、患者Pに息止めを指令するための息止め指令部とを備えている。
【0012】
静磁場発生部は、例えば超電導方式の磁石1と、この磁石1に電流を供給する静磁場電源2とを含み、被検体Pが遊挿される円筒状の開口部(診断用空間)の軸方向(Z軸方向)に静磁場H0を発生させる。なお、この磁石部にはシムコイル14が設けられている。このシムコイル14には、後述するホスト計算機の制御下で、シムコイル電源15から静磁場均一化のための電流が供給される。寝台部は、被検体Pを載せた天板を磁石1の開口部に退避可能に挿入できる。
【0013】
傾斜磁場発生部は、傾斜磁場コイルユニット3を含む。この傾斜磁場コイルユニット3は、互いに直交するX、YおよびZ軸方向の傾斜磁場を発生させるための3組のコイル3x、コイル3y,コイル3zを備える。傾斜磁場発生部はまた、コイル3x〜3zに電流を供給する傾斜磁場電源4を含む。傾斜磁場電源4は、後述するシーケンサ5の制御のもと、傾斜磁場を発生させるためのパルス電流をコイル3x〜3zに供給する。
【0014】
傾斜磁場電源4からコイル3x〜3zに供給されるパルス電流を調整することにより、物理軸であるX,Y,Z方向の各軸の傾斜磁場を合成して、互いに直交するスライス方向傾斜磁場Gs、位相エンコード方向傾斜磁場Ge、および読出し方向(周波数エンコード方向)傾斜磁場Grの各論理軸方向を任意に設定することができる。スライス方向、位相エンコード方向および読出し方向の各傾斜磁場は、静磁場H0に重畳される。
【0015】
送受信部は、磁石1内の撮影空間にて被検体Pの近傍に配設されるRFコイル7と、このコイル7に接続された送信器8Tおよび受信器8Rとを含む。送信器8Tおよび受信器8Rは、後述するシーケンサ5の制御の下で動作する。送信器8Tは、核磁気共鳴(NMR)を起こさせるためのラーモア周波数のRFパルスをRFコイル7に供給する。受信器8Rは、RFコイル7が受信したエコー信号(高周波信号)を取り込み、これに前置増幅、中間周波変換、位相検波、低周波増幅、フィルタリングなどの各種の信号処理を施した後、A/D変換してエコー信号に応じたデジタル量のエコーデータ(原データ)を生成する。
【0016】
制御・演算部は、シーケンサ(シーケンスコントローラとも呼ばれる)5、ホスト計算機6、演算ユニット10、記憶ユニット11、表示器12、入力器13および音声発生器16を含む。このうち、ホスト計算機6は、記憶したソフトウエア手順により、シーケンサ5にパルスシーケンス情報を指令するとともに、装置全体の動作を統括する機能を有する。
【0017】
ホスト計算機6は、位置決め用スキャンなどの準備作業に引き続いて、イメージングスキャンを実施する。イメージングスキャンは、画像再構成に必要なエコーデータの組を収集するスキャンであり、ここでは2次元スキャンに設定されている。イメージングスキャンは、ECG信号に依るECGゲート法を併用して行うことができる。なお、このECGゲート法は場合によっては併用しなくても良い。
【0018】
シーケンサ5は、CPUおよびメモリを備えており、ホスト計算機6から送られてきたパルスシーケンス情報を記憶し、この情報にしたがって傾斜磁場電源4、送信器8Tおよび受信器8Rのそれぞれの動作を制御するとともに、受信器8Rが出力したエコーデータを一旦入力し、これを演算ユニット10に転送する。パルスシーケンス情報は、一連のパルスシーケンスにしたがって傾斜磁場電源4、送信器8Tおよび受信器8Rを動作させるために必要な全ての情報であり、例えばコイル3x〜3zに印加するパルス電流の強度、印加時間、印加タイミングなどに関する情報を含む。
【0019】
演算ユニット10は、受信器8Rが出力したエコーデータをシーケンサ5を介して入力する。演算ユニット10は、その内部メモリ上のフーリエ空間(k空間または周波数空間とも呼ばれる)にエコーデータを配置し、このエコーデータを各組毎に2次元または3次元のフーリエ変換に付して実空間の画像データに再構成する。また演算ユニットは、再構成された画像に基づく解析パラメータ画像の作成や合成処理などの画像処理を行う機能を有する。
【0020】
記憶ユニット11は、再構成された画像データのみならず、上述の合成処理や差分処理が施された画像データを保管することができる。表示器12は、ホスト計算機6の制御の下に画像を表示する。入力器13を介して、術者が希望する撮影条件、パルスシーケンス、画像合成や差分演算に関する情報をホスト計算機6に入力できる。
【0021】
息止め指令部の一要素として音声発生器16を備える。この音声発生器16は、ホスト計算機6から指令の下に、息止め開始および息止め終了のメッセージを音声として発することができる。
【0022】
心電計測部は、被検体の体表に付着させてECG信号を電気信号として検出するECGセンサ17と、このセンサ信号にデジタル化処理を含む各種の処理を施してホスト計算機6およびシーケンサ5に出力するECGユニット18とを含む。心電計測部による計測信号は、イメージングスキャンを実行するときにシーケンサ5により用いられる。これにより、ECGゲート法(心電同期法)による同期タイミングを適切に設定でき、この同期タイミングに基づくECGゲート法のイメージングスキャンを行ってデータ収集できるようになっている。
【0023】
次に以上のように構成されたMRI装置の動作について詳細に説明する。
【0024】
(A)理論的準備
具体的な動作について説明する前の理論的準備として、不均一磁化率と流れのあるボクセルのMR信号モデルと、リフェーズ(Rephase)/ディフェーズ(Dephase)の信号モデルとについて定義する。
【0025】
(A-1)不均一磁化率と流れのあるボクセルのMR信号モデル
まず、不均一磁化率と流れのあるボクセルのMR信号モデルに関するパラメータを次のように定義する。
M0:プロトン密度
AT1:T1に依存した減衰(T1 dependent attenuation)
AT1=1−exp(−TR/T1)
AT2:T2に依存した減衰(T2 dependent attenuation)
AT2=exp(−TE/T2)
AD:拡散に依存した減衰(Diffusion dependent attenuation)
AD=exp[−bD]
Asus:磁化率に依存した減衰(Susceptibility dependent attenuation)
Asus=exp[−TE(γΔB0σ)]*
φsus:磁化率に依存した位相(Susceptibility dependent phase)
φsus=−TE(γΔB0m)
但しLorenzian modelの場合。
T2*:T2と磁化率効果による成分を含む緩和時間
1/T2*=1/T2+γΔB0 (T2*<T2)
なお、T2*を用いると、AT2 Asus=exp(−TE/T2*)となる。
Aflow:フローに依存した減衰(Flow dependent attenuation)
Aflow=exp[−bDflow]
b:拡散で定義される傾斜磁場パターンから算出される係数(gradient factor)
Dflow:フローによる位相分散係数(流れのない部分で定義される拡散係数に等価な係数に相当する)
Φflow:フローに依存した位相シフト(Flow dependent phase shift)
【数1】
【0026】
V:不均一磁化率かつ流れを含む組織。ΔB0<>0であり、かつF<>0である組織であり、例えば静脈などがこれに相当する。
【0027】
この場合、組織VからのMR信号Sは、次の式(1)のように一般化される。
S=(M0 AT1 AT2 AD Asus Aflow) exp[i(Φ0+Φsus+Φflow)] …(1)
ここで磁化率やフロー(flow)の効果のうち、コヒーレント(coherent)な成分は位相に、インコヒーレント(incoherent)な成分は振幅減衰に生じる。すなわち、コントラストには、ボクセル内でのΔB0分布が大きいほど、あるいはIVIM(intravoxel incoherent motion)成分が多い(例えば、静脈→細静脈→毛細血管の順に大きい)ほど、振幅項が位相項より支配的に作用することになる。
【0028】
(A-2)リフェーズ(Rephase)/ディフェーズ(Dephase)の信号モデル
次に、リフェーズ/ディフェーズの信号モデルについて説明する。
【0029】
理想的なモデルでは、リフェーズでは振幅・位相ともフローによる成分はキャンセルされ、磁化率による成分のみとなる。ただし、実際はリフェーズでもGMN(gradient moment nulling)オーダーや乱流に依存した変化が生じる。一方、ディフェーズでは、振幅・位相とも磁化率成分にフロー成分が加わるので、リフェーズ/ディフェーズの各々の理想的なモデルでの信号Sre,Sdeは下記の式(2)および式(3)のように表される。
Sre=(M0 AT1 AT2 AD Asus)exp[i(Φ0+Φsus)] …(2)
Sde=(M0 AT1 AT2 AD Asus Aflow)exp[i(Φ0+Φsus+Φflow)] …(3)
ここでもし、リフェーズ/ディフェーズのシーケンス条件のうちの繰り返し時間TRおよびエコー時間TEが同一で、ディフェーズのb値は十分小さくてADの効果は無視できるとすれば、Sdeは次の式(4)となる。
Sde=Sre Aflow exp[i Φflow] …(4)
すなわち、SdeはSreにフロー効果が加わったものになる。
【0030】
MR信号(振幅・位相)の変動には次のような性質がある。
・ボクセル内のフロー成分比Mflow/(Mflow+Mst)が多いほど大きい。
・b値が大きいほど大きい。
・IVIM成分が大きければ、dephase grad.による振幅減衰効果が大きい。
・IVCM成分が大きければ、dephase grad.による位相変化効果が大きい。
(B)ディフェーズ/リフェーズのデータ収集および処理法
従来のリフェーズのみやディフェーズのみの場合に対し、ディフェーズとリフェーズを組み合わせて用いることで新たな応用が可能になる。上述のように、リフェーズはT1、T2緩和と磁化率効果のみによる信号変化であるが、ディフェーズはリフェーズに加えフローによる効果が加わっている。従って、リフェーズとディフェーズとのデータ間の演算により、静止部分と磁化率効果などフロー効果以外の成分とフロー効果が分離された画像、さらには磁化率やフロー効果のより定量的な普遍的なパラメータ画像を得ることができる。臨床的には、ディフェーズおよびリフェーズの物理的なフロー効果および磁化率効果の大小を用いた、静脈と動脈の分離などである。ここでは収集の工夫やその臨床応用例をあげる。
【0031】
(B-1)理論
具体的な処理手順を説明するのに先立って、理論について説明する。
【0032】
リフェーズおよびディフェーズのMR信号は、振幅および位相を有する複素信号であり、MR信号のゲインをK、A0=M0、AT1、AT2として、式(2),(3)を任意のゲインKとシーケンスの可変パラメータであるTR,TE,b値の関数として、以下の式(2')および式(3')のように表せる。
Sre(TR,TE)=Are(TR,TE) exp[iΦre(TE)]=K A0(TR,TE) Asus(TE) exp[i{Φ0(TE)+Φsus(TE)}] …(2')
Sde(TE,TE,b)=Ade(TE,TE) exp[iΦde(TE)]=K A0(TE,TE) Asus(TE) Aflow(b) exp[i{Φ0(TE)+Φsus(TE)+Φflow(b)}] …(3')
なお、ここではb値はフローによる信号減衰効果をもたらせば良いため十分小さく、A0に占める分子拡散の効果ADは無視できる。すなわちA0は、リフェーズ,ディフェーズで同一とした。また縦磁化Mzは、M0 AT1として求まり、AT1はGRE(gradient echo)の場合にはインフロー効果も含めTRで決まる。
【0033】
(B-1-1)生体組織と磁化率、フロー
磁化率とフローの性質の違いを生体内の組織別にみると、動脈はオキシヘモグロビン(oxyHb)を多く含むため磁化率効果が無視できるので、リフェーズでは動脈は描出されない。一方、静脈は動脈よりはフロー効果は小さく、さらにデオキシヘモグロビン(deoxyHb)による磁化率効果が加わる。従って、血管部については、リフェーズではフローによる位相変化が完全にリフォーカスされれば(集まれば)、静脈信号(磁化率効果成分)が支配的に描出され、ディフェーズでは動脈および静脈(磁化率効果+フロー成分)のいずれもが描出される。静止部については、リフェーズおよびディフェーズのいずれでも同等に描出される。
【0034】
静脈の選択的描出にはフローが遅いため、磁化率効果以外に良い方法がない。また、リフェーズのみでは動脈と静止部との分離はコントラスト差が十分でないため困難である。
【0035】
(B-1-2)最適シーケンスパラメータ:TE
TEについては、振幅、位相ともにCNR(contrast to noise ratio)を最大にするTEの最適な条件は、TE=T2*である。すなわち、対象とする組織のT2*に等しい場合に、その近傍のT2*を有する組織との間で最大のCNRを与える。収集条件や人体組織の種類に依存してT2*には分布があるが、ボクセルサイズを一定にし、血管描出に重要な血液も磁化率の範囲があるし、白質や灰白質はほぼ一定とみなせるのでほぼ決めることは可能である。
【0036】
(B-1-3)フロー効果と磁化率効果の分離
受信ゲイン、TR、TEおよびbなどの撮像パラメータを固定すれば、そのままでもリフェーズ・ディフェーズ間の相互の比較が可能であるが、より定量的な指標にするには、フロー成分はDflow、磁化率成分はT2*やΔχなどから計算して表示したほうが良い。その場合は、b=0はリフェーズで代用し、TEは2つ以上と最低3画像必要となる。
【0037】
(B-1-3-1)フロー効果の定量化
ディフェーズ/リフェーズ間の振幅比と位相差をとると、(4)式に基づいて次の式(5)および式(6)に示すようにフロー効果を分離できる。
ディフェーズ/リフェーズ振幅比:A(de)/A(re)=Aflow …(5)
ディフェーズ/リフェーズ位相差:Φ(de)−Φ(re)=Φflow …(6)
ここで、位相差ΦflowはPS-MRA(phase contrast MR angiography)法のように速度エンコードパルスVENCの3方向用いて測定し求めても良いが、ここではあまり現実的ではないので省略し、次の式(7)のように表される関係にあり、多様な流速や多様な方向への流れによるランダムな位相分散の程度を表すDflowを求める。
Aflow=exp[−b* Dflow] …(7)
ここで、ディフェーズのbファクタ(b-factor)をb(de)、ディフェーズおよびリフェーズの信号強度(振幅)をそれぞれA(de)、A(re)とすると、b=0をリフェーズで代用できるので、Dflowは次の式(8)により算出される。
Dflow[mm2/sec]=−ln[A(de)/A(re)]/b(de) …(8)
(B-1-3-2)磁化率効果の定量化
まず磁化率効果の定量化指標を振幅を用いて算出する場合について述べる。
【0038】
磁化率効果による減衰項Asusは、TE依存であり、単一TEではT1緩和およびT2緩和などの寄与が入ったA0が消去できないので、γΔB0σまたはT2を含んだT2*を求めたほうがより普遍化できる。
【0039】
リフェーズでの2つのTE(以下、それぞれをTE1,TE2と称し、TE2>TE1であることとする)を用いる場合、TE1およびTE2のそれぞれでの信号強度(振幅)A1,A2は、次の式(9)および式(10)により定まる。
A1=K A0 exp[−TE1/T2*] …(9)
A2=K A0 exp[−TE2/T2*] …(10)
これら式(9)および式(10)より、磁化率効果による振幅減衰効果の指標は次の式(11)で表される。
T2*=(TE2−TE1)/ln[A1(TE1)/A2(TE2)] …(11)
またT2*の算出には、3エコー以上のマルチエコーを取得して最小二乗近似を用いても良い。
【0040】
一方、位相を用いる場合について述べる。
【0041】
磁化率効果の別の表現としてリフェーズの単一TEまたは2エコーの位相差から低周波の位相成分を減じた位相からΦsusを求める。単一エコーの場合および2エコーの場合のΦsusは、それぞれ次の式(12)および式(12')により表される。
Φsus=−2πγΔχB0(cos2θ−1/3)TE …(12)
Φsus=−2πγΔχB0(cos2θ−1/3)(TE2−TE1) …(12')
これに基づいて、単一エコーの場合および2エコーの場合のΔχは、それぞれ次の式(13)および式(13')により表される。
Δχ[ppm]=−Φsus/{2πγ B0(cos2θ−1/3)TE} …(13)
Δχ[ppm]=−Φsus/{2πγ B0(cos2θ−1/3)(TE2−TE1)} …(13')
これはボクセル内の平均の磁化率をあらわすとみなせる。
【0042】
ここで、γは磁気回転比、B0は静磁場強度を表し、装置により一意に決まる。しかし、θはB0方向と血管の走行角度であるので、ボクセル間の相関を調べるなどの方法で測定する必要があり、振幅から求めるよりは困難になる。
【0043】
なお、磁化率分布のみならば、リフェーズのみでTEが2段階以上あれば実現可能であり、この技術に関しては公知である。
【0044】
次に、具体的な処理手順について図2に従って説明する。
【0045】
(B-2)ディフェーズデータおよびリフェーズデータの収集
ステップSa1においてはシーケンサ5が、ディフェーズデータおよびリフェーズデータを収集する。
【0046】
同一TEのディフェーズおよびリフェーズのシーケンスを独立させて連続的に撮像しても良い。あるいは動きの影響を最小にするために、1データセットを複数のセグメントに分割して交互に収集して合成しても良い。分割するセグメントは、k−spaceで1ライン(TR)単位としたり、2D面単位とするなど任意である。
【0047】
b値の段階は複数でも良い。TEも複数取得すればT2*の正確な算出にも役立つ。複数のTEは、Gradient Echo法ならTR内でマルチエコーを収集すれば一回で取得できる。また複数のb値もLook−Locker法を用いれば一回で取得できる。なおLook−Locker法は、「"Measurement of Gd-DTPA dialysis clearance rates by using a look-locker imaging technique.",Magn Reson,Med. 1996 Oct;36(4):571-8.」により知られている。
【0048】
(B-2-1)ディフェーズ・リフェーズ交互分割収集(ディフェーズ,リフェーズの一組は同一TE)
フロー効果Dflowのみの算出のためには、シーケンサ5はリフェーズおよびディフェーズを同一TEで2画像収集する。Dflowに加えてT2*を算出する場合には、シーケンサ5はリフェーズ側を2エコーでTE1、TE2とし、ディフェーズはTE1,TE2のいずれかにする。その場合、ディフェーズおよびリフェーズのシーケンスを分割させて別々に撮像しても良い。あるいは、画像間の動きの影響を最小にするために、1データセットを複数のセグメントに分割して交互に収集して合成しても良い。分割するセグメントは、k−spaceで1ライン(TR)単位としたり、2D面単位とするなど任意である。
【0049】
(B-2-2)GREマルチエコーによるリフェーズ,ディフェーズ混合連続収集の場合(TEがすべて違う)
GREでリフェーズおよびディフェーズ組み合わせた2点以上のマルチエコーで処理する実施例を示す。
【0050】
上記の(B-2-1)のようにリフェーズおよびディフェーズで同一TEになるように別々のデータとして取得する場合に比べ、マルチエコーで収集すれば、リフェーズおよびディフェーズのTEは同一にできないが、複数のTEで収集し計算により同一にでき、定量的なパラメータが得られる。マルチエコーは1回のRF励起で、すなわち同一のTR内で取得できるため、スキャン時間が1エコーの収集時間と同程度であるという最大のメリットがあるのに加え、独立に時間をおいて収集する場合に比べてデータ間の体動も無視できるため、異種データ間の演算時の誤差が軽減できるメリットもある。もちろん別々のデータとして取得しても、後述する解析パラメータ画像の算出処理は共通にできる。
【0051】
(B-2-2-1)2点法
(B-2-2-1-1)リフェーズTEとディフェーズTEとを1点ずつ収集する場合
図3に示すように、TE=TE1のAde(TE1)と、TE=TE2のAre(TE2)とをそれぞれ収集する。
【0052】
この場合は、同一モードのデータが2点以上ないため、2つの画像をそのまま観察するのが主体で、定量パラメータのうちT2*,Dflowなどは算出できない。しかし、条件をTE2=2*TE1に設定すれば、流れの位相項Φflowのみは定量化可能である。
【0053】
TE(de)<TE(re)とすることで、ディフェーズは磁化率効果を抑制したフロー強調になり、リフェーズはフロー効果を抑制した磁化率強調になる。TE(re)<TE(de)といった具合にTE(re)を短く(<10ms)設定すれば、リフェーズは動脈描出用として通常ルーチンで収集しているTOF−MRA(time of flight-magnetic resonance angiography)の代用となり、通常のTOF−MRAが不要になる可能性がある。さらに、ディフェーズは、フローと磁化率効果も強調した動静脈描出用にできるし、TEの設定しだいで静脈の描出は制御できる。また位相情報も用いれば背景組織に対する静脈のCNRを強調できる。さらに、後述するようにディフェーズとリフェーズとの間で相互に演算処理を施すことにより、動脈と静脈の分離表示も可能である。なお、2エコーの場合でもリフェーズおよびディフェーズのTEは厳密に合わせる必要はなく、十分近くT2*の差が大きくなければ半定量的な表示は可能である。
【0054】
(B-2-2-1-2)ディフェーズTEを2点収集する場合
図4に示すように、TE=TE1に関するAde(TE1)およびTE=TE2に関するAde(TE2)をそれぞれ収集する。
【0055】
De TE=TE1,TE2 (TE1<TE2)の2つのエコーのb値を同じにすれば、Aflow=exp[−bDflow]、またGradient momentも同じにすれば位相のフロー効果分Φflowもほぼ同じであるから、次の式(14)からフロー効果がキャンセルされることが明らかである。
【0056】
Sde(TE2,b)/Sde(TE1,b)=exp[−(TE2−TE1)/T2*] exp[−i(TE2−TE1)γΔB0m] …(14)
そしてTE1,TE2は既知なので、式(13)の振幅からT2*、さらには位相項からΔB0をそれぞれ算出可能であり、Δχが算出可能である。
【0057】
(B-2-2-1-3)リフェーズTEを2点収集する場合
図5に示すように、TE=TE1に関するAre(TE1)およびTE=TE2に関するAre(TE2)をそれぞれ収集する。
【0058】
A2でb=0とみなせるので、Aflow=1 Φflow=0となる以外は(B-2-2-1-2)における添え字deをreに置き換えたものになる。
【0059】
(B-2-2-2)3点法
リフェーズおよびディフェーズをまぜて3点収集すれば、T2*またはΔB0に加えて、すなわちΔχに加えて、フロー効果分も算出可能であるし、またGREのマルチエコーを用いるにもかかわらずに同一TEのリフェーズ画像およびディフェーズ画像が作成可能である。3点のTEは、どう組み合わせても良い。
【0060】
3点法では各エコーがリフェーズおよびディフェーズのいずれであるかの組み合わせで2×2×2=8通り、順番を問題にしなければ4通りがある。このうちリフェーズ,リフェーズ,ディフェーズと、ディフェーズ,ディフェーズ,リフェーズとの2通りを取り上げる。
【0061】
(B-2-2-2-1)リフェーズTEを2点とディフェーズTEを1点収集する場合
最初のリフェーズの2点よりT2*を算出し、これを3点目のディフェーズに代入してDflowを算出する。またリフェーズよりT2*が分かれば、リフェーズの任意のTEの信号を生成可能なので、ディフェーズと同一TEのリフェーズが求まる。すなわち同一のT2緩和効果および磁化率効果でフロー効果のみ異なる画像が得られる。
【0062】
このときのGREマルチエコーシーケンスの一例を図6に示す。
【0063】
(B-2-2-2-2)ディフェーズTEを2点とリフェーズTEを1点収集する場合
上記の(B-2-2-1-2)と同様にして、最初のb値が同じでTEが異なるディフェーズの2点よりT2*を算出し、3点目のリフェーズとからDflowを算出する。またディフェーズよりT2*が分かればディフェーズの任意のTEの信号を生成可能なので、リフェーズと同一TEのディフェーズが求まる。すなわち同一のT2緩和効果および磁化率効果でフロー効果のみ異なる画像が得られる。
【0064】
(B-2-2-3)4点以上法
ディフェーズ、リフェーズとも各々2点以上ずつ収集し、T2*,同一TEの画像を求める。未知数2個で4点以上では最小二乗近似となる。もちろん必要ならば未知パラメータのM0、T2やDなども算出可能である。
【0065】
図7は4点法におけるデータ収集の一例を示す図であり、TE=TE1に関するAre(TE1)、TE=TE2に関するAde(TE2)、TE=TE3に関するAre(TE3)、TE=TE4に関するAde(TE4)をそれぞれ収集する。
【0066】
なお、以上に説明した各種のシーケンスモードのいずれにおいても、マルチエコーのTEを、T2*算出に用いる最適TE(=T2*)を含むように設定する。対象のT2*が長すぎたり(時間延長、SNR低下)、短かすぎたり(RFが入らない、傾斜磁場が出ない、十分なb値が出ないなど)で最適TEの設定困難な場合は、計算して作る任意のTEの画像からTE=T2*の画像を作成すれば良い。
【0067】
また、位相算出においては、すべての場合で折り返しがないか、折り返しを補正することが好ましい。GREのマルチエコーでもリフェーズおよびディフェーズのTEは厳密にあわせる必要はなく、十分近いTEでT2*の差が大きくなければ半定量的な表示は可能である。シーケンスもGREのみならず、1回のRF励起後に実効TEを数段階変えたmulti−shot EPI(echo planar imaging)でk−spaceをセグメント分割するなどして短時間化と高分解能化を適当にコントロールしても良い。
【0068】
またTEに関しては、水と脂肪とが混在する組織では、水脂肪が同位相になるように設定することが重要である。静磁場強度のもとでの水(プロトン)とケミカルシフトδppmの物質との位相差は、ΔΦ=2πγδB0TEと表されるが、それが同位相となる条件はnを整数としてδΦ=n2πなので、TE=n/(γδB0)の倍数に設定すれば良い。脳実質では脂肪があまりないためにあまり問題とはならないが、ボクセル内で水と脂肪とが混在する骨髄や腹部臓器では問題となる場合があるのでその条件が必要になる。γ=42.6MHz/T、脂肪はγ=3.6ppm、B0=1.5TとするとTE=n 4.3 msとなる。さらに脂肪以外でも磁化率の異なる物質がボクセル内に混在すると位相差が生じるが、酸素濃度によるケミカルシフトはδ=0.1ppm程度とされ、TEを適当に選んでも定量化パラメータが振幅のT2*ではあまり問題とならないし、位相では厳密に定量化する場合は問題となり得るので、必要ならδを既知としてその寄与による位相を求め補正すれば良い。
【0069】
(B-3)画像再構成
ステップSa2においては演算ユニット10が、以上のような各種の手法を適宜に採用して収集されたデータのそれぞれを使用して周知の再構成処理を行うことによって、m枚のディフェーズ画像とn枚のリフェーズ画像とを再構成する。なお、mおよびnの値は、いずれも0を含む整数で、採用するデータ収集の手法により定まる。
【0070】
(B-4)解析パラメータ画像の算出
ステップSa3においては演算ユニット10が、再構成された1乃至m枚のディフェーズ画像および1乃至n枚のリフェーズ画像を使用して解析パラメータ画像を算出する。
【0071】
図8は解析パラメータ画像の算出処理の概念を示す図である。
【0072】
演算ユニット10は、リフェーズ画像およびディフェーズ画像を使用してフローパラメータ算出処理P1を行うことによってよってDflow画像を算出する。演算ユニット10は、リフェーズ画像およびディフェーズ画像を使用して磁化率パラメータ算出処理P2を行うことによってT2*画像およびΔχ画像を算出する。演算ユニット10は、リフェーズ画像、ディフェーズ画像およびT2*画像を使用して任意TE画像作成処理P3を行うことによって、任意TEに関するリフェーズ画像およびディフェーズ画像を算出する。
【0073】
以下に、上記した各種のデータ収集手法のそれぞれに応じた解析パラメータ画像の算出処理の具体例をそれぞれ説明する。
【0074】
(B-4-1)2点法
(B-4-1-1)ディフェーズ2点の場合
(a) 演算ユニット10は、同一b値のディフェーズに関するTE=TE1の振幅画像A1(TE1)およびTE=TE2の振幅画像A2(TE2)より、次の式(15)および式(16)によりT2*およびKdeをそれぞれ算出する。
T2*=(TE2−TE1)/ln[Ade(TE1)/Ade(TE2)] …(15)
Kde=Ade(TE1)/exp[−TE1/T2*] …(16)
(b) 演算ユニット10は、任意TEに関して、次の式(17)および式(18)によりディフェーズ振幅および位相をそれぞれ算出する。
Ade(TE)=Kde exp[−TE/T2*] …(17)
Φde(TE)=(TE/TE1)Φde(TE1) …(18)
(B-4-1-2)リフェーズ2点の場合
(a) 演算ユニット10は、同一b値のディフェーズに関するTE=TE1の振幅画像A1(TE1)およびTE=TE2の振幅画像A2(TE2)より、次の式(19)および式(20)によりT2*およびKreをそれぞれ算出する。
T2*=(TE2−TE1)/ln[Are(TE1)/Are(TE2)] …(19)
Kre=Are(TE1)/exp[−TE1/T2*] …(20)
(b) 演算ユニット10は、任意TEに関して、次の式(21)および式(22)によりディフェーズ振幅および位相をそれぞれ算出する。
Are(TE)=Kre exp[−TE/T2*] …(21)
Φre(TE)=(TE/TE1)Φre(TE1) …(22)
すなわち、(B-4-1)における添え字deをreに置き換えた式を適用する。
【0075】
(B-4-2)3点法
(B-4-2-1)リフェーズTEを2点とディフェーズTEを1点収集する場合
この場合の処理の流れを図9に示す。
【0076】
(a) ステップSb1において演算ユニット10は、リフェーズに関するTE=TE1の振幅画像Are(TE1)およびTE=TE2の振幅画像Are(TE2)より、次の式(23)および式(24)によりT2*およびKreをそれぞれ算出する。
T2*=(TE2−TE1)/ln[Are(TE1)/Are(TE2)] …(23)
Kre=Are(TE1)/exp[−TE1/T2*] …(24)
ここで求めたT2*により、T2*画像が得られる。
【0077】
(b) ステップSb2において演算ユニット10は、次の式(25)により、例えば図10に示すようなリフェーズに関するTE=TE3の振幅画像Are(TE3)を作成する。また、ステップSb3において演算ユニット10は、TEに依存しないバックグラウンド位相がない場合には式(26)により、TEに依存しないバックグラウンド位相がある場合には式(27)により、位相画像Φre(TE3)を作成する。
Are(TE3)=Kre exp[−TE3/T2*] …(25)
Φre(TE3)=(TE3/TE1)Φre(TE1) …(26)
Φre(TE3)={TE3/(TE2−TE1)} {Φre(TE2)−Φre(TE1)} …(27)
なお、ステップSb2およびステップSb3において演算ユニット10は、TEがTE1,TE2,TE3のいずれでもない任意の値であるリフェーズの振幅画像および位相画像を作成することもできる。これは、上記の式(25)乃至(27)におけるTE3に任意の欲しいTEを代入することにより計算することにより実現できる。
【0078】
(c) ステップSb4において演算ユニット10は、TE=TE3の振幅画像Are(TE3),Ade(TE3)より、次の式(28)および式(29)によりフロー分散係数Dflow画像を算出する。
Aflow=Ade(TE3)/Are(TE3) …(28)
Dflow[mm2/sec]=−ln[Aflow]/bde …(29)
(d) ステップSb5において演算ユニット10は、血管とB0との方向角度θを算出する。
【0079】
(e) さらにステップSb5において演算ユニット10は、θとTE=TE3のリフェーズの位相画像より、次の式(30)によりΔχ[ppm]を算出する。
Δχ[ppm]=−Φre(TE3)/{2πγB0(cos2θ−1/3)TE3} …(30)
ここで求めたΔχにより、Δχ画像が得られる。
【0080】
(f) ステップSb6において演算ユニット10は、任意TEのディフェーズの振幅画像を作成する。
ディフェーズの任意TEの振幅画像は式(28)により求めたAflowを用いて、次の式(31)により計算する。
Ade(TE)=Are(TE)Aflow …(31)
ステップSb7において演算ユニット10は、Φde(TE3)を用いて、次の式(32)によりディフェーズの任意TEの位相画像を算出する。これは、TEに依存しないバックグラウンド位相がない場合のみ算出可能である。
Φde(TE)=(TE/TE3)Φde(TE3) …(32)
なお、上記の(c)〜(f)で求める各種パラメータは、必要な場合にのみ算出すれば良い。
【0081】
(B-4-2-2)ディフェーズTEを2点とリフェーズTEを1点収集する場合
この場合の処理の流れを図11に示す。
【0082】
(a) ステップSc1において演算ユニット10は、ディフェーズに関するTE=TE1の振幅画像Ade(TE1)およびTE=TE2の振幅画像Ade(TE2)より、次の式(33)および式(34)によりT2*およびKdeを算出する。
T2*=(TE2−TE1)/ln[Ade(TE1)/Ade(TE2)] …(33)
Kde=Ade(TE1)/exp[−TE1/T2*] …(34)
ここで求めたT2*により、T2*画像が得られる。
【0083】
(b) ステップSc2において演算ユニット10は、次の式(35)により、図12に示すようにディフェーズに関するTE=TE3の振幅画像Ade(TE3)を作成する。また、ステップSc3において演算ユニット10は、次の式(36)により、位相画像Φde(TE3)を作成する。
Ade(TE3)=Kdeexp[−TE3/T2*] …(35)
Φde(TE3)=(TE3/TE1)Φde(TE1) …(36)
なお、ステップSc2およびステップSc3において演算ユニット10は、TEがTE1,TE2,TE3のいずれでもない任意の値であるディフェーズの振幅画像および位相画像を作成することもできる。これは、上記の式(35)および(36)におけるTE3に任意の欲しいTEを代入することにより計算することにより実現できる。
【0084】
(c) ステップSc4において演算ユニット10は、TE=TE3の振幅画像Are(TE3),Ade(TE3)より、次の式(38)および式(39)によりフロー分散係数Dflow画像を算出する。
Aflow=Ade(TE3)/Are(TE3) …(37)
Dflow[mm2/sec]=−ln[Aflow]/bde …(38)
(d) ステップSc5において演算ユニット10は、血管とB0との方向角度θを算出する。
【0085】
(e) さらにステップSc5において演算ユニット10は、θとTE=TE3の位相画像より、次の式(39)によりΔχ[ppm]を算出する。
Δχ[ppm]=Φre(TE3)/{2πγB0(cos2θ−1/3)TE3} …(39)
ここで求めたΔχにより、Δχ画像が得られる。
【0086】
(f) ステップSc6において演算ユニット10は、任意TEのリフェーズの振幅画像を作成する。
リフェーズの任意TEの振幅画像は式(37)により求めたAflowを用いて、次の式(40)により計算する。
Are(TE)=Ade(TE)/Aflow …(40)
ステップSc7において演算ユニット10は、Φre(TE3)を用いて、次の式(41)によりリフェーズの任意TEの位相画像を算出する。これは、TEに依存しないバックグラウンド位相がない場合のみ算出可能である。
Φre(TE)=(TE/TE3)Φre(TE3) …(41)
なお、(c)〜(f)で求めるパラメータは、必要な場合にのみ算出すれば良い。
【0087】
さらに、リフェーズおよびディフェーズの各々が3点以上の場合は、未知パラメータをモデルに即して最小二乗法で算出する。Kre,Kde,T2*は1次の指数減衰関数のモデルを用い、Δχは一次関数のモデルを用いれば良い。
【0088】
以上のようにして、少なくとも1つずつのリフェーズ画像およびディフェーズ画像(それぞれ振幅画像および位相画像を含む)、T2*画像、Δχ画像およびDflow画像などがそれぞれ算出される。なお以下においては、これらの画像を解析パラメータ画像と総称する。
【0089】
なお、マルチエコーではなく、1エコーずつ別々に収集して画像化したものを使用しても良い。あるいは、2Dマルチスライス収集や、3DFT法などによるボリューム収集でも良い。さらに、パルスシーケンスとしては、GREに代えて非対称スピンエコー(ASE)法を利用しても良い。
【0090】
(B-5)解析パラメータ画像の合成処理および表示
各種の解析パラメータ画像を算出した後、これらの解析パラメータ画像をそのまま表示して観察に供することができる。また、Aflow、Φflow、Δflowなどのパラメータ自体を表示しても良い。
【0091】
さらに診断の補助とするためにステップSa4において演算ユニット10は、各種の解析パラメータ画像のカラーで分離した2Dの合成画像を作成して、2D画像のまま表示することができる。このステップSa4の合成処理は、必要時だけに行うようにしても良い。
【0092】
あるいは、特に血管画像を目的とする場合にはステップSa5において演算ユニット10は、血管を連続した管として表現するために3D処理を行う。その代表的なものは、最大値投影(MIP)または最小値投影(minIP)である。演算ユニット10は、ステップSa5にて複数の3D画像を生成した上で、ステップSa6にてこれら複数の3D画像を用いてフュージョン処理を行うこともできる。
【0093】
(B-5-1)位相合成処理および表示
任意TEのリフェーズ画像およびディフェーズ画像より位相合成画像を作成する機能について説明する。
【0094】
(a) 通常のMR画像:S=A exp[iΘ]と、ローパスフィルタをかけたMR画像:S1=A1 exp[iΘ1]とをリフェーズおよびディフェーズのおのおのについて作成する。
【0095】
(b) 位相アーチファクト補正
次の式(42)によりリフェーズのMR画像の位相マップよりローパスフィルタをかけた画像の位相マップを減算する。
Φ=Θ-Θl=arg[S]−arg[S1] …(42)
この式(42)の計算後、位相飛び補正により−π<Φ<=πに収める処理を行う。すなわち、Φが−π以下であるならばΦ+πに、またΦがπよりも大きい場合にはΦ−πにΦを置き換える。
【0096】
なお、Φは、式(42)に代えて次の式(43)により算出しても良い。
Φ=arg[S/S1] …(43)
(c) 位相マスク(phase masking)
図13に示す非対称型(asymmetric type)または図14に示す対称型(symmetric type)のどちらかのマスクMを作成し、振幅画像とマスクMとのn重の積Iを次の式(44)により求める。
I=A×Mn …(44)
なお、Aはディフェーズの位相から、Mはリフェーズの位相から作成したものを用いる。
【0097】
これにより位相飛びやフローと磁化率由来の位相キャンセルのない位相マスクが可能になる。
【0098】
(B-5-2)カラー合成処理および表示
上記のように得られた解析パラメータ画像をカラーで色分けしてフュージョン表示する。例えば、フロー成分のAflow、ΦflowおよびDflowは動脈および静脈に支配的に大きくなり、磁化率成分のT2*およびΔχは静脈や、静止部なら出血などの部分で大きくなる。そこで図15に示す画像21のようなフロー成分の画像を赤、画像22のような磁化率成分の画像を青とするなどの色分けをしてカラー合成処理を行うことにより、例えば画像21としてフュージョン表示する。このときに表示されるのは、物理的にはフロー成分(赤)と磁化率成分(青)との合成画像である。血管部のみを表す画像の場合には、動脈は赤により、静脈は紫により表される。磁化率成分部は、出血や静止部の磁化率がゼロでない部分のアーチファクトも含まれるので、必要なら閾値処理などによる静脈抽出を行った後にフュージョン表示すれば良い。
【0099】
また、画像21および画像22から次の式(45)または式(46)によって動脈および静脈を抽出した画像24を作成して、この画像24を画像22とカラー合成することによって画像23を得ることもできる。
Aflow=Ade/Are …(45)
Φflow=Φde−Φre …(46)
色の混じり具合は、フロー成分と磁化率成分の比で決まる。上記の色の割り当て例では、赤→紫→青の変化を示すが、赤に近いほどフロー成分が、青に近いほど磁化率成分が多いことを意味する。脳梗塞などの疾患時は必ずしも各々動脈および静脈に対応せず、酸素代謝状態を反映する指標となる。色の割り当てはこれに限らず、対比できる組み合わせならなんでも良い。TOF−MRAで動脈を描出する方法もあるが、原理的に遅い動脈は描出されない。本法ではディフェーズによる信号減衰または位相変化を用いるため、動脈も細いのや、上方から回り込んでくる側副血行路も描出することができ、脳梗塞の診断などに向けた臨床的に重要な情報を提供できる。またある程度の長いTEに設定すれば、磁化率効果により血栓や出血などの情報も同時に得られるので、治療計画を立てる上でも重要な情報を提供できる。
【0100】
なお、画像21としては、ディフェーズ振幅画像、ディフェーズ位相画像、ディフェーズの振幅位相合成画像およびDflow画像などが利用できる。画像22としては、リフェーズ振幅画像、リフェーズ位相画像、リフェーズの振幅位相合成画像およびT2*画像などが利用できる。
【0101】
(B-5-3)3D処理および表示
ディフェーズの振幅画像では、動脈および静脈のいずれも周囲組織に対して低い画像値になるのでminIPが妥当である。リフェーズの振幅画像では、静脈は低信号なのでminIPが妥当だが、動脈はTOF(Time of Flight)効果で高信号になるためにMIPを併用しても良い。表面を抽出するなどしてボリュームレンダリングやサーフェスレンダリングも用いることができる。また原画像信号を見たい場合は、単純に断面変換(MPR)も用途により使用することが有効である。フュージョン表示は、前記のように2Dのままでも良いが、3D画像を作成してから行っても良い。
【0102】
なお、各種画像の表示は、表示器12により行うことが可能であるが、外部のビューワ装置などにおいて表示させても良い。
【0103】
以上のように本実施形態によれば、血流などのフローの効果を分離して定量化することができる。そしてこの定量化結果に基づいて、医用診断に有用な情報や画像を提供することが可能である。
【0104】
本実施形態によれば、磁化率やフローに関して、機種やシーケンスに依存しないパラメータが作成可能である。このため、共通のデータベースになり、重要な知見がエビデンスとして蓄積され易くなる。
【0105】
本実施形態によれば、リフェーズとディフェーズとを組み合わせることで動脈と静脈との分離が可能である。
【0106】
本実施形態によれば、リフェーズおよびディフェーズを別々に収集する方法に比べ、マルチエコーを用いることにより1回の収集で双方の画像が取得可能であるので、収集時間が短くなり、かつ複数画像間の時間差が無視できるのでより動きの影響が生じにくい。
【0107】
本実施形態によれば、2画像を用いるため、速い流速の血液のtime of fright (TOF)効果が抑制可能である。血液は、磁化率またはフローのみに依存したコントラストが得られる。
【0108】
本実施形態によれば、側副血行路の描出が可能である。
【0109】
本実施形態によれば、磁化率とフローとが分離可能である。
【0110】
本実施形態によれば、ディフェーズで得られた振幅画像では、血管内孔が選択的に描出されるという手術などでは際めて重要な特長を有する。
【0111】
この実施形態は、次のような種々の変形実施が可能である。
【0112】
位相マスク画像を作成する場合には、ディフェーズの振幅とリフェーズの位相とを、あるいはディフェーズの位相とリフェーズの振幅とを組み合わせてもよい。このようにすることで、位相の磁化率とフローとの相殺が生じないため、振幅の位相マスク画像における静脈信号がより低下し、動脈と静脈の分離能が向上する。
【0113】
撮像時間とSNR向上を目的としてk-spaceにおける周波数毎にTRまたはTEを可変しても良い。例えば、シーケンスによる収集時、低周波は短いTEで、中間から高周波は長いTEで収集する。その間は、TEを滑らかに変える。これにより、低周波成分が支配的な静磁場不均一による位相は小さくなるのでアーチファクトが低減する。位相の場合には縦磁化は無関係なので、TRも最短にして良い。なおシーケンス種はGRE(FE)タイプで、マルチエコーおよび1エコーのいずれでも適用可能で、k-space trajectoryは、spin warp、spiral、あるいはEPIなどのいずれでも良い。図16にkに対するTR,TEの変化の例を示す。
【0114】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【符号の説明】
【0115】
1…磁石、2…静磁場電源、3…傾斜磁場コイルユニット、3x,3y,3z…コイル、4…傾斜磁場電源、5…シーケンサ、6…ホスト計算機、7…RFコイル、8R…受信器、8T…送信器、10…演算ユニット、11…記憶ユニット、12…表示器、13…入力器、14…シムコイル、15…シムコイル電源、16…音声発生器、17…ECGセンサ、18…ECGユニット。
【特許請求の範囲】
【請求項1】
被検体から放射される磁気共鳴信号を収集する収集手段と、
収集された前記磁気共鳴信号に基づいてディフェーズ画像およびリフェーズ画像を少なくとも1枚ずつ再構成する再構成手段と、
再構成された前記ディフェーズ画像および前記リフェーズ画像の双方に基づいて前記被検体に関する特性を定量化する定量化手段と、
前記定量化された特性を表す複数種の定量化画像を生成する生成手段とを具備したことを特徴とする磁気共鳴イメージング装置。
【請求項2】
前記生成手段は、前記複数種の定量化画像の1つとしてフロー成分を表す画像を生成するとともに、前記複数種の定量化画像の1つとして磁化率成分を表す画像を生成することを特徴とする請求項1に記載の磁気共鳴イメージング装置。
【請求項3】
前記定量化画像を表示する手段をさらに具備することを特徴とする請求項1または請求項2に記載の磁気共鳴イメージング装置。
【請求項4】
前記収集手段は、同一のエコー時間でディフェーズおよびディフェーズの磁気共鳴信号を、3D k-spaceにおける1ライン単位または面単位で交互に収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項5】
前記収集手段は、グラディエントエコー(GRE)のマルチエコー法によって前記磁気共鳴信号を収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項6】
前記収集手段は、エコー時間を変化させながら前記磁気共鳴信号を収集することを特徴とする請求項5に記載の磁気共鳴イメージング装置。
【請求項7】
前記収集手段は、1回のRF励起後にエコープラナーイメージング(EPI)法で連続収集するLook-Locker法により前記ディフェーズの磁気共鳴信号を収集することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項8】
前記収集手段は、複数のエコー時間のそれぞれに関するディフェーズの磁気共鳴信号をそれぞれ収集することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項9】
前記収集手段は、複数のエコー時間のそれぞれに関するリフェーズの磁気共鳴信号をそれぞれ収集することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項10】
前記収集手段は、少なくとも1つのエコー時間に関するディフェーズの磁気共鳴信号と、別の少なくとも1つのエコー時間に関するリフェーズの磁気共鳴信号とを収集することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項11】
前記収集手段は、1つのエコー時間に関するディフェーズの磁気共鳴信号および2つのエコー時間に関するリフェーズの磁気共鳴信号、あるいは1つのエコー時間に関するリフェーズの磁気共鳴信号および2つのエコー時間に関するディフェーズの磁気共鳴信号を収集することを特徴とする請求項10に記載の磁気共鳴イメージング装置。
【請求項12】
前記収集手段は、2つ以上のエコー時間に関するディフェーズの磁気共鳴信号と、別の少なくとも2つ以上のエコー時間に関するリフェーズの磁気共鳴信号を収集することを特徴とする請求項10に記載の磁気共鳴イメージング装置。
【請求項13】
前記収集手段は、複数のエコー時間のうちの1つをT2*に設定することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項14】
前記収集手段は、水脂肪が同位相となるようにエコー時間を設定することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項15】
前記再構成手段は、複数のエコー時間でそれぞれ収集された磁気共鳴信号にそれぞれ基づいて再構成した複数のディフェーズ画像または複数のリフェーズ画像に基づいて、前記複数のエコー時間とは異なるエコー時間に関するディフェーズ画像またはリフェーズ画像を再構成することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項16】
前記定量化手段は、前記被検体におけるフロー効果を定量化することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項17】
前記定量化手段は、前記フロー効果を定量化するために、多様な流速や多様な方向への流れによるランダムな位相分散の程度を表す位相分散係数を算出することを特徴とする請求項16に記載の磁気共鳴イメージング装置。
【請求項18】
再構成された前記ディフェーズ画像および前記リフェーズ画像の少なくともいずれか一方に基づいて前記被検体に関する特性を定量化する第2の定量化手段をさらに備え、
前記生成手段は、前記定量化手段による定量化の結果を表す画像と前記第2の定量化手段による定量化の結果を表す画像とを合成して前記定量化画像を生成することを特徴とする請求項2に記載の磁気共鳴イメージング装置。
【請求項19】
前記収集手段は、2次元または3次元で磁気共鳴信号を収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項20】
前記収集手段は、グラディエントエコー法を用いて、リフェーズおよびディフェーズのそれぞれに関して同一のエコー時間の磁気共鳴信号を収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項21】
前記収集手段は、非対称スピンエコー法を用いて、リフェーズおよびディフェーズのそれぞれに関して同一のエコー時間の磁気共鳴信号を収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項22】
前記定量化手段は、前記ディフェーズ画像と前記リフェーズ画像との振幅比または位相差を算出することによって前記被検体に関する特性を定量化することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項23】
前記生成手段は、振幅単独、位相単独、あるいは振幅位相合成したスライス画像を複数生成するとともに、この複数のスライス画像に基づく3D処理によって3Dの定量化画像を生成することを特徴とする請求項2または請求項3に記載の磁気共鳴イメージング装置。
【請求項24】
前記生成手段は、振幅と位相とを別シーケンス種とすることを特徴とする請求項23に記載の磁気共鳴イメージング装置。
【請求項25】
前記生成手段は、振幅はディフェーズのものを使用し、位相はリフェーズのものを使用することを特徴とする請求項24に記載の磁気共鳴イメージング装置。
【請求項26】
前記生成手段は、前記リフェーズ画像と前記ディフェーズ画像との位相合成画像を生成することを特徴とする請求項2または請求項3に記載の磁気共鳴イメージング装置。
【請求項27】
前記生成手段は、それぞれが定量化された特性を表す複数の定量化画像どうし、あるいは前記リフェーズ画像または前記ディフェーズ画像に基づく画像と前記定量化画像とをカラー合成した画像を生成することを特徴とする請求項2または請求項3に記載の磁気共鳴イメージング装置。
【請求項28】
前記生成手段は、定量化されたフロー成分を表す定量化画像と、前記リフェーズ画像または前記ディフェーズ画像に基づく磁化率成分の画像とをカラー合成した画像を生成することを特徴とする請求項27に記載の磁気共鳴イメージング装置。
【請求項29】
前記収集手段は、k空間での周波数毎に繰り返し時間またはエコー時間を変化させることを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項30】
前記収集手段は、k空間での周波数が低い場合には高い場合に比べて繰り返し時間またはエコー時間を短くすることを特徴とする請求項29に記載の磁気共鳴イメージング装置。
【請求項1】
被検体から放射される磁気共鳴信号を収集する収集手段と、
収集された前記磁気共鳴信号に基づいてディフェーズ画像およびリフェーズ画像を少なくとも1枚ずつ再構成する再構成手段と、
再構成された前記ディフェーズ画像および前記リフェーズ画像の双方に基づいて前記被検体に関する特性を定量化する定量化手段と、
前記定量化された特性を表す複数種の定量化画像を生成する生成手段とを具備したことを特徴とする磁気共鳴イメージング装置。
【請求項2】
前記生成手段は、前記複数種の定量化画像の1つとしてフロー成分を表す画像を生成するとともに、前記複数種の定量化画像の1つとして磁化率成分を表す画像を生成することを特徴とする請求項1に記載の磁気共鳴イメージング装置。
【請求項3】
前記定量化画像を表示する手段をさらに具備することを特徴とする請求項1または請求項2に記載の磁気共鳴イメージング装置。
【請求項4】
前記収集手段は、同一のエコー時間でディフェーズおよびディフェーズの磁気共鳴信号を、3D k-spaceにおける1ライン単位または面単位で交互に収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項5】
前記収集手段は、グラディエントエコー(GRE)のマルチエコー法によって前記磁気共鳴信号を収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項6】
前記収集手段は、エコー時間を変化させながら前記磁気共鳴信号を収集することを特徴とする請求項5に記載の磁気共鳴イメージング装置。
【請求項7】
前記収集手段は、1回のRF励起後にエコープラナーイメージング(EPI)法で連続収集するLook-Locker法により前記ディフェーズの磁気共鳴信号を収集することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項8】
前記収集手段は、複数のエコー時間のそれぞれに関するディフェーズの磁気共鳴信号をそれぞれ収集することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項9】
前記収集手段は、複数のエコー時間のそれぞれに関するリフェーズの磁気共鳴信号をそれぞれ収集することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項10】
前記収集手段は、少なくとも1つのエコー時間に関するディフェーズの磁気共鳴信号と、別の少なくとも1つのエコー時間に関するリフェーズの磁気共鳴信号とを収集することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項11】
前記収集手段は、1つのエコー時間に関するディフェーズの磁気共鳴信号および2つのエコー時間に関するリフェーズの磁気共鳴信号、あるいは1つのエコー時間に関するリフェーズの磁気共鳴信号および2つのエコー時間に関するディフェーズの磁気共鳴信号を収集することを特徴とする請求項10に記載の磁気共鳴イメージング装置。
【請求項12】
前記収集手段は、2つ以上のエコー時間に関するディフェーズの磁気共鳴信号と、別の少なくとも2つ以上のエコー時間に関するリフェーズの磁気共鳴信号を収集することを特徴とする請求項10に記載の磁気共鳴イメージング装置。
【請求項13】
前記収集手段は、複数のエコー時間のうちの1つをT2*に設定することを特徴とする請求項6に記載の磁気共鳴イメージング装置。
【請求項14】
前記収集手段は、水脂肪が同位相となるようにエコー時間を設定することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項15】
前記再構成手段は、複数のエコー時間でそれぞれ収集された磁気共鳴信号にそれぞれ基づいて再構成した複数のディフェーズ画像または複数のリフェーズ画像に基づいて、前記複数のエコー時間とは異なるエコー時間に関するディフェーズ画像またはリフェーズ画像を再構成することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項16】
前記定量化手段は、前記被検体におけるフロー効果を定量化することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項17】
前記定量化手段は、前記フロー効果を定量化するために、多様な流速や多様な方向への流れによるランダムな位相分散の程度を表す位相分散係数を算出することを特徴とする請求項16に記載の磁気共鳴イメージング装置。
【請求項18】
再構成された前記ディフェーズ画像および前記リフェーズ画像の少なくともいずれか一方に基づいて前記被検体に関する特性を定量化する第2の定量化手段をさらに備え、
前記生成手段は、前記定量化手段による定量化の結果を表す画像と前記第2の定量化手段による定量化の結果を表す画像とを合成して前記定量化画像を生成することを特徴とする請求項2に記載の磁気共鳴イメージング装置。
【請求項19】
前記収集手段は、2次元または3次元で磁気共鳴信号を収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項20】
前記収集手段は、グラディエントエコー法を用いて、リフェーズおよびディフェーズのそれぞれに関して同一のエコー時間の磁気共鳴信号を収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項21】
前記収集手段は、非対称スピンエコー法を用いて、リフェーズおよびディフェーズのそれぞれに関して同一のエコー時間の磁気共鳴信号を収集することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項22】
前記定量化手段は、前記ディフェーズ画像と前記リフェーズ画像との振幅比または位相差を算出することによって前記被検体に関する特性を定量化することを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項23】
前記生成手段は、振幅単独、位相単独、あるいは振幅位相合成したスライス画像を複数生成するとともに、この複数のスライス画像に基づく3D処理によって3Dの定量化画像を生成することを特徴とする請求項2または請求項3に記載の磁気共鳴イメージング装置。
【請求項24】
前記生成手段は、振幅と位相とを別シーケンス種とすることを特徴とする請求項23に記載の磁気共鳴イメージング装置。
【請求項25】
前記生成手段は、振幅はディフェーズのものを使用し、位相はリフェーズのものを使用することを特徴とする請求項24に記載の磁気共鳴イメージング装置。
【請求項26】
前記生成手段は、前記リフェーズ画像と前記ディフェーズ画像との位相合成画像を生成することを特徴とする請求項2または請求項3に記載の磁気共鳴イメージング装置。
【請求項27】
前記生成手段は、それぞれが定量化された特性を表す複数の定量化画像どうし、あるいは前記リフェーズ画像または前記ディフェーズ画像に基づく画像と前記定量化画像とをカラー合成した画像を生成することを特徴とする請求項2または請求項3に記載の磁気共鳴イメージング装置。
【請求項28】
前記生成手段は、定量化されたフロー成分を表す定量化画像と、前記リフェーズ画像または前記ディフェーズ画像に基づく磁化率成分の画像とをカラー合成した画像を生成することを特徴とする請求項27に記載の磁気共鳴イメージング装置。
【請求項29】
前記収集手段は、k空間での周波数毎に繰り返し時間またはエコー時間を変化させることを特徴とする請求項1乃至請求項3のいずれか1項に記載の磁気共鳴イメージング装置。
【請求項30】
前記収集手段は、k空間での周波数が低い場合には高い場合に比べて繰り返し時間またはエコー時間を短くすることを特徴とする請求項29に記載の磁気共鳴イメージング装置。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【公開番号】特開2012−86056(P2012−86056A)
【公開日】平成24年5月10日(2012.5.10)
【国際特許分類】
【出願番号】特願2012−173(P2012−173)
【出願日】平成24年1月4日(2012.1.4)
【分割の表示】特願2006−315823(P2006−315823)の分割
【原出願日】平成18年11月22日(2006.11.22)
【出願人】(000003078)株式会社東芝 (54,554)
【出願人】(594164542)東芝メディカルシステムズ株式会社 (4,066)
【Fターム(参考)】
【公開日】平成24年5月10日(2012.5.10)
【国際特許分類】
【出願日】平成24年1月4日(2012.1.4)
【分割の表示】特願2006−315823(P2006−315823)の分割
【原出願日】平成18年11月22日(2006.11.22)
【出願人】(000003078)株式会社東芝 (54,554)
【出願人】(594164542)東芝メディカルシステムズ株式会社 (4,066)
【Fターム(参考)】
[ Back to top ]