(57)【要約】【課題】体液の流動による信号強度の変化を容易に把握することができる。【解決手段】流速画像作成部２２ｃが、被検体内を流れる体液の流速成分が得られるＥＰＩを複数回繰り返すことによって得られた複数の画像それぞれについて流速成分の分布を表す流速画像を作成する。また、流速分散画像作成部２２ｄが、作成された複数の流速画像を用いて、各流速画像の同一位置ごとに時系列に沿った流速成分の分散を算出し、算出した流速成分の分散の分布を表す流速分散画像を作成する。そして、重畳画像処理部２２ｅが、流速分散画像における流速成分の分散の分布を平均絶対値画像に重畳させ、画像表示制御部２６ａが、重畳画像を表示する。【選択図】        図２
2011173020.tif 000002
【特許請求の範囲】
【請求項１】
トリガ信号に同期して被検体の撮像領域に励起パルスを印加して該撮像領域に含まれるスピンに関するエコー信号を収集し、収集したエコー信号に基づいて磁気共鳴画像を生成する磁気共鳴イメージング装置であって、
前記被検体内を流れる脳脊髄液を得るための撮像を複数回繰り返すことによって脳脊髄液に関する複数の再構成画像を作成する画像作成手段と、
前記画像作成手段によって作成された複数の再構成画像を用いて、各再構成画像における同一位置ごとに前記脳脊髄液の流れに関する情報を算出し、当該算出された情報に基づく計算画像を作成する計算画像作成手段と、
を備えたことを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記撮像は、前記脳脊髄液の流速成分を得るためのフローエンコード傾斜磁場を含むことを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記撮像は、前記脳脊髄液にラベリングパルスを印加することによって得られたエコー信号を収集することを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【請求項４】
前記撮像は、一回の励起パルスで一つの画像を再構成するためのエコー信号を収集することを特徴とする請求項１〜３のいずれか一つに記載の磁気共鳴イメージング装置。
【請求項５】
前記撮像は、前記ラベリングパルスを印加してから前記励起パルスを印加するまでの時間を固定して一つの画像を再構成するために必要なエコー信号を収集する動作を複数回繰り返すことによって、脳脊髄液に関する複数の再構成画像を作成することを特徴とする請求項３に記載の磁気共鳴イメージング装置。
【請求項６】
前記計算画像作成手段は、前記複数の再構成画像における同一位置ごとに前記脳脊髄液の流速の最大値または最小値を算出することを特徴とする請求項１〜３のいずれか一つに記載の磁気共鳴イメージング装置。
【請求項７】
前記計算画像作成手段は、前記複数の再構成画像における同一位置ごとに前記脳脊髄液の流速を大きい順に並べた場合に全体に対して上位から所定の割合に入る流速成分の最大値または最小値を算出することを特徴とする請求項１〜３のいずれか一つに記載の磁気共鳴イメージング装置。
【請求項８】
前記計算画像作成手段は、前記複数の再構成画像における同一位置ごとの信号強度の分散を算出することを特徴とする請求項１〜３のいずれか一つに記載の磁気共鳴イメージング装置。
【請求項９】
被検体内を流れる脳脊髄液を得るための撮像を複数回繰り返すことによって作成された脳脊髄液に関する複数の再構成画像を取得する画像取得手段と、
前記画像取得手段によって取得された複数の再構成画像を用いて、各再構成画像における同一位置ごとに前記脳脊髄液の流れに関する情報を算出し、当該算出された情報に基づく計算画像を作成する計算画像作成手段と、
を備えたことを特徴とする画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、被検体内を撮像した画像を処理する磁気共鳴イメージング装置および画像処理装置に関する。
【背景技術】
【０００２】
従来、磁気共鳴イメージング装置を用いてＣＳＦや血液等の体液の動態を観測する方法として、フローエンコード傾斜磁場を利用した「フェイズシフト(phase shift)法」または「フェイズコントラスト法」と呼ばれる方法がある（例えば、非特許文献１参照）。
【０００３】
図１４は、従来のフェイズシフト法のパルスシーケンスを示す図である。図１４に示すように、通常、フェイズシフト法では、グラジエントエコー法が用いられ、励起パルス（図１４に示すＰ１）とエコー信号（図１４に示すＰ２）との間でフローエンコード傾斜磁場（Ｇｆｅ）（図１４に示すＰ３）を加えることによって、流速に比例した位相変化が画像に与えられる。
【０００４】
このフェイズシフト法は、基本的には、対象の流速変化が心電波形（図１４に示すＥＣＧ）に強く相関していることを前提としている。すなわち、図１４に示すように、心電同期波形のＲ波から一定の待ち時間（図１４に示すＴｄｅｌａｙ）が経過した後に高周波励起が行われて、エコー信号が収集される。ここで、通常、１つの画像を再構成するためには、１２８個から２５６個程度のエコー信号が必要である。そのため、一般的には、フェイズシフト法では、１回のＲ波ごとに１つのエコー信号を収集する手順が、位相エンコード傾斜磁場パルスを段階的に変化させながら１２８回から２５６回程度繰り返される。
【０００５】
こうして収集されたエコー信号に離散フーリエ変換などの再構成処理を施すことによって得られた画像Ｓ（ｒａ）は、流速に比例した位相変化を有する。この画像Ｓ（ｒａ）に対して、装置の不完全性や静磁場不均一性による補正処理を施した後に、例えばＶ（ｒａ）＝ｋ（ｖｅｎｃ）・ａｒｇ｛Ｓ（ｒａ）｝で表される処理（ここで、ｋ（ｖｅｎｃ）は、フローエンコード傾斜磁場のパルスの形状によって変化する比例係数）を施すことによって、流速の分布を表す流速画像Ｖ（ｒａ）が得られる。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Matt A.Bernstein、Kevin F.King、Xiaohong Joe Zhou 著、「Handbook of MRI Pulse Sequences」、Elsevier Academic Press、ｐ．６５９−６７７、２００４年
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、従来技術では、体液の流動による信号強度の変化を容易に把握することができない場合があった。
【０００８】
この発明は、上述した従来技術による課題を解決するためになされたものであり、体液の流動による信号強度の変化を容易に把握することができる磁気共鳴イメージング装置および画像処理装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上述した課題を解決し、目的を達成するため、本発明は、トリガ信号に同期して被検体の撮像領域に励起パルスを印加して該撮像領域に含まれるスピンに関するエコー信号を収集し、収集したエコー信号に基づいて磁気共鳴画像を生成する磁気共鳴イメージング装置であって、前記被検体内を流れる脳脊髄液を得るための撮像を複数回繰り返すことによって脳脊髄液に関する複数の再構成画像を作成する画像作成手段と、前記画像作成手段によって作成された複数の再構成画像を用いて、各再構成画像における同一位置ごとに前記脳脊髄液の流れに関する情報を算出し、当該算出された情報に基づく計算画像を作成する計算画像作成手段と、を備えたことを特徴とする。
【００１０】
また、本発明は、被検体内を流れる脳脊髄液を得るための撮像を複数回繰り返すことによって作成された脳脊髄液に関する複数の再構成画像を取得する画像取得手段と、前記画像取得手段によって取得された複数の再構成画像を用いて、各再構成画像における同一位置ごとに前記脳脊髄液の流れに関する情報を算出し、当該算出された情報に基づく計算画像を作成する計算画像作成手段と、を備えたことを特徴とする。
【発明の効果】
【００１１】
本発明によれば、体液の流動による信号強度の変化を容易に把握することができるという効果を奏する。
【図面の簡単な説明】
【００１２】
【図１】図１は、本実施例１に係るＭＲＩ装置の全体構成を示す図である。
【図２】図２は、図１に示した画像再構成部、記憶部および制御部の詳細を示す機能ブロック図である。
【図３】図３は、本実施例１で用いられるＥＰＩ法のパルスシーケンスを示す図である。
【図４】図４は、図３に示したパルスシーケンスによる撮像の流れを示す図である。
【図５】図５は、本実施例１に係るＭＲＩ装置による画像処理の処理手順を示すフローチャートである。
【図６】図６は、本実施例１に係るＭＲＩ装置による画像作成の流れを示す図である。
【図７】図７は、本実施例１で用いられるパルスシーケンスの変形例を示す図である。
【図８】図８は、本実施例２における画像再構成部、記憶部および制御部の詳細を示す機能ブロック図である。
【図９】図９は、本実施例２で用いられるＴｉｍｅ−ＳＬＩＰ法のパルスシーケンスを示す図である。
【図１０】図１０は、本実施例２におけるＴｉｍｅ−ＳＬＩＰ法によって行われるラベリングの位置の一例を示す図である。
【図１１】図１１は、本実施例２における絶対値画像作成部によって作成される絶対値画像の一例を示す図である。
【図１２】図１２は、本実施例２に係るＭＲＩ装置による画像処理の処理手順を示すフローチャートである。
【図１３】図１３は、本実施例２に係るＭＲＩ装置による画像作成の流れを示す図である。
【図１４】図１４は、従来のフェイズシフト法のパルスシーケンスを示す図である。
【発明を実施するための形態】
【００１３】
以下に添付図面を参照して、この発明に係る磁気共鳴イメージング装置および画像処理装置の好適な実施例を詳細に説明する。なお、以下に示す実施例では、本発明を磁気共鳴イメージング装置（以下、「ＭＲＩ（Magnetic Resonance Imaging）装置」と呼ぶ）に適用した場合について説明する。
【実施例１】
【００１４】
まず、本実施例１に係るＭＲＩ装置の全体構成について説明する。図１は、本実施例１に係るＭＲＩ装置の全体構成を示す図である。図１に示すように、本実施例１に係るＭＲＩ装置１００は、静磁場磁石１、傾斜磁場コイル２、傾斜磁場電源３、寝台４、寝台制御部５、送信ＲＦコイル６、送信部７、受信ＲＦコイル８、受信部９、シーケンス制御部１０、架台部１５および計算機システム２０を備える。
【００１５】
静磁場磁石１は、中空の円筒形状に形成された磁石であり、内部の空間に一様な静磁場を発生する。この静磁場磁石１としては、例えば永久磁石、超伝導磁石等が使用される。
【００１６】
傾斜磁場コイル２は、中空の円筒形状に形成されたコイルであり、静磁場磁石１の内側に配置される。この傾斜磁場コイル２は、互いに直交するＸ，Ｙ，Ｚの各軸に対応する３つのコイルが組み合わされて形成されており、これら３つのコイルは、後述する傾斜磁場電源３から個別に電流供給を受けて、Ｘ，Ｙ，Ｚの各軸に沿って磁場強度が変化する傾斜磁場を発生させる。なお、Ｚ軸方向は、静磁場と同方向とする。傾斜磁場電源３は、傾斜磁場コイル２に電流を供給する装置である。
【００１７】
ここで、傾斜磁場コイル２によって発生するＸ，Ｙ，Ｚ各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅおよびリードアウト用傾斜磁場Ｇｒにそれぞれ対応している。スライス選択用傾斜磁場Ｇｓは、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場Ｇｅは、空間的位置に応じてエコー信号（磁気共鳴信号）の位相を変化させるために利用される。リードアウト用傾斜磁場Ｇｒは、空間的位置に応じてエコー信号の周波数を変化させるために利用される。
【００１８】
寝台４は、被検体Ｐが載置される天板４ａを備えた装置であり、後述する寝台制御部５による制御のもと、天板４ａを、被検体Ｐが載置された状態で傾斜磁場コイル２の空洞（撮像口）内へ挿入する。通常、この寝台４は、長手方向が静磁場磁石１の中心軸と平行になるように設置される。寝台制御部５は、寝台４を制御する装置であり、計算機システム２０による制御のもと、寝台４を駆動して天板４ａを長手方向および上下方向へ移動する。
【００１９】
送信ＲＦコイル６は、傾斜磁場コイル２の内側に配置されたコイルであり、送信部７から高周波パルスの供給を受けて、高周波磁場を発生する。送信部７は、ラーモア周波数に対応する高周波パルスを送信ＲＦコイル６に送信する装置である。
【００２０】
受信ＲＦコイル８は、傾斜磁場コイル２の内側に配置されたコイルであり、上記の高周波磁場の影響によって被検体Ｐから発せられるエコー信号を受信する。この受信ＲＦコイル８は、エコー信号を受信すると、そのエコー信号を受信部９へ出力する。
【００２１】
受信部９は、受信ＲＦコイル８から出力されるエコー信号に基づいてｋ空間データを生成する装置である。具体的には、この受信部９は、受信ＲＦコイル８から出力されるエコー信号をデジタル変換することによってｋ空間データを生成する。このｋ空間データには、前述したスライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅおよびリードアウト用傾斜磁場Ｇｒによって、ＰＥ方向、ＲＯ方向、ＳＥ方向の空間周波数の情報が対応付けられている。そして、ｋ空間データを生成すると、受信部９は、そのｋ空間データをシーケンス制御部１０へ送信する。
【００２２】
シーケンス制御部１０は、計算機システム２０から送信されるシーケンス情報に基づいて、傾斜磁場電源３、送信部７および受信部９を駆動することによって、被検体Ｐのスキャンを行う装置である。ここで、シーケンス情報とは、傾斜磁場電源３が傾斜磁場コイル２に供給する電源の強さや電源を供給するタイミング、送信部７が送信ＲＦコイル６に送信するＲＦ信号の強さやＲＦ信号を送信するタイミング、受信部９がエコー信号を検出するタイミングなど、スキャンを行うための手順を定義した情報である。
【００２３】
なお、シーケンス制御部１０は、傾斜磁場電源３、送信部７および受信部９を駆動して被検体Ｐをスキャンした結果、受信部９からｋ空間データが送信されると、そのｋ空間データを計算機システム２０へ転送する。
【００２４】
架台部１５は、前述した静磁場磁石１、傾斜磁場コイル２、送信ＲＦコイル６および受信ＲＦコイル８を有し、静磁場中に置かれた被検体に高周波磁場を照射し、当該被検体から発せられるＮＭＲ信号を検出する。
【００２５】
計算機システム２０は、ＭＲＩ装置１００の全体制御や、データ収集、画像再構成などを行う装置であり、インタフェース部２１、画像再構成部２２、記憶部２３、入力部２４、表示部２５および制御部２６を有している。
【００２６】
インタフェース部２１は、シーケンス制御部１０との間で授受される各種信号の入出力を制御する。例えば、このインタフェース部２１は、シーケンス制御部１０に対してシーケンス情報を送信し、シーケンス制御部１０からｋ空間データを受信する。ｋ空間データを受信すると、インタフェース部２１は、各ｋ空間データを被検体Ｐごとに記憶部２３に格納する。
【００２７】
画像再構成部２２は、記憶部２３によって記憶されたｋ空間データに対して、後処理、すなわちフーリエ変換処理等の再構成処理を施すことによって、被検体Ｐ内における所望核スピンのスペクトラムデータあるいは画像データを生成する。なお、かかる画像再構成部２２の構成については、後に詳細に説明する。
【００２８】
記憶部２３は、インタフェース部２１によって受信されたｋ空間データと、画像再構成部２２によって生成された画像データなどを、被検体Ｐごとに記憶する。なお、かかる記憶部２３の構成については、後に詳細に説明する。
【００２９】
入力部２４は、操作者からの各種指示や情報入力を受け付ける装置である。この入力部２４としては、マウスやトラックボールなどのポインティングデバイス、モード切替スイッチ等の選択デバイス、あるいはキーボード等の入力デバイスを適宜に利用可能である。
【００３０】
表示部２５は、制御部２６による制御のもと、スペクトラムデータあるいは画像データ等の各種の情報を表示する装置である。この表示部２５としては、液晶表示器などの表示デバイスを利用可能である。
【００３１】
制御部２６は、図示していないＣＰＵやメモリ等を有し、ＭＲＩ装置１００の全体制御を行う。具体的には、この制御部２６は、入力部２４を介して操作者から入力される撮像条件に基づいてシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御部１０に送信することによってスキャンを制御したり、スキャンの結果としてシーケンス制御部１０から送られるｋ空間データに基づいて行われる画像の再構成を制御したりする。なお、かかる制御部２６の構成については、後に詳細に説明する。
【００３２】
次に、図１に示した画像再構成部２２、記憶部２３および制御部２６の構成について説明する。図２は、図１に示した画像再構成部２２、記憶部２３および制御部２６の構成を示す機能ブロック図である。
【００３３】
図２に示すように、記憶部２３は、特に、ｋ空間データ記憶部２３ａおよび画像データ記憶部２３ｂを有する。ｋ空間データ記憶部２３ａは、インタフェース部２１によって受信されたｋ空間データを記憶する。画像データ記憶部２３ｂは、画像再構成部２２によって生成された画像データを記憶する。
【００３４】
画像再構成部２２は、特に、フーリエ変換処理部２２ａ、絶対値画像作成部２２ｂ、流速画像作成部２２ｃ、流速分散画像作成部２２ｄおよび重畳画像処理部２２ｅを有する。
【００３５】
フーリエ変換処理部２２ａは、ｋ空間データ記憶部２３ａに記憶されたｋ空間データに離散２次元フーリエ変換などの再構成処理を施すことによって画像を再構成する。
【００３６】
ここで、本実施例１では、被検体内を流れる体液の流速成分が得られる撮像法による撮像を複数回繰り返すことによって得られた複数の画像が用いられる。具体的には、本実施例１では、１回の励起パルスで１画像を再構成し、かつ、再構成した画像に対して流速に比例した位相変化を与えることが可能なＥＰＩ法による撮像を被検体の心電波形に同期して複数回繰り返すことによって得られた複数の画像が用いられる。
【００３７】
図３は、本実施例１で用いられるＥＰＩ法のパルスシーケンスを示す図である。図３に示すように、本実施例１で用いられるＥＰＩ法では、心電波形（図３に示すＥＣＧ）のＲ波から一定の待ち時間（図３に示すＴｄｅｌａｙ）が経過した後に、９０度励起パルス（図３に示すＰ１１）を印加することによって撮像断面内が励起される。この９０度励起パルスとエコー信号収集部分の間（図３に示すＰ１３）でフローエンコード傾斜磁場（図３に示すＰ１２）を加えることによって、得られる画像の位相が撮像対象の流速に比例して変化する。
【００３８】
なお、ここでは、ＳＥ（Spin Echo）タイプのＥＰＩを一例として示したが、例えば、ＦＥ（Field Echo）タイプのＥＰＩや、その他、位相画像が得られる撮像法であればいずれの撮像法が用いられてもよい。望ましくは、なるべく短い時間で画像が得られる撮像方法を用いるのがよい。例えば、「Proc. of Annual Meeting, Society of Magnetic Resonance in Medicine, No.2138, 1998.」に記載されている撮像法もその一例である。撮像方法や装置の状態によっては、静磁場不均一性や渦磁場による位相変化を低減するため、位相補正や他のデータ補正処理を施すのが望ましい。
【００３９】
図４は、図３に示したパルスシーケンスによる撮像の流れを示す図である。図４に示すように、図３に示したパルスシーケンスによれば、心電波形（図４に示すＥＣＧ）のＲ波から一定の待ち時間（図４に示すＴｄｅｌａｙ）が経過した後に、パルスシーケンスの１ショット分のスキャンが行われる。
【００４０】
本実施例１では、フーリエ変換処理部２２ａは、上記で説明したように１ショットで得られたエコー信号のｋ空間データそれぞれに離散２次元フーリエ変換などの再構成処理を施すことによって、図４に示すように複数の画像Ｓ（ｒａ，ｉ）を再構成する。ここで、ｒａは、画像上の位置を示すベクトルであり、ｉは、収集の順番を示す添え字である。各画像Ｓ（ｒａ，ｉ）では、流速に比例して画素の位相が変化している。
【００４１】
絶対値画像作成部２２ｂは、フーリエ変換処理部２２ａによって再構成された画像から絶対値画像を作成する。具体的には、絶対値画像作成部２２ｂは、フーリエ変換処理部２２ａによって再構成された複数の画像Ｓ（ｒａ，ｉ）それぞれについて、以下に示す式（１）で与えられる絶対値画像Ｉ（ｒａ，ｉ）を作成する。
【００４２】
Ｉ（ｒａ，ｉ）＝ａｂｓ｛Ｓ（ｒａ，ｉ）｝        ・・・（１）
【００４３】
さらに、絶対値画像作成部２２ｂは、作成した複数の絶対値画像Ｉ（ｒａ，ｉ）をもとに、以下に示す式（２）で与えられる平均絶対値画像Ｉａｖｇ（ｒａ）を作成する。
【００４４】
【数１】
2011173020.tif 000003
【００４５】
流速画像作成部２２ｃは、フーリエ変換処理部２２ａによって再構成された画像から流速成分の分布を含む計算画像を作成する。本実施例では、流速画像作成部２２ｃは、計算画像として、流速成分の分布を表す流速画像を作成する。具体的には、流速画像作成部２２ｃは、フーリエ変換処理部２２ａによって再構成された複数の画像Ｓ（ｒａ，ｉ）それぞれについて、以下に示す式（３）で与えられる流速画像Ｖ（ｒａ，ｉ）を作成する。なお、式（３）において、ｋ（ｖｅｎｃ）は、フローエンコード傾斜磁場のパルスの形状によって変化する比例係数である。
【００４６】
Ｖ（ｒａ，ｉ）＝ｋ（ｖｅｎｃ）・ａｒｇ｛Ｓ（ｒａ，ｉ）｝        ・・・（３）
【００４７】
さらに、流速画像作成部２２ｃは、作成した複数の流速画像Ｖ（ｒａ，ｉ）をもとに、以下に示す式（４）で与えられる平均流速画像Ｖａｖｇ（ｒａ）を作成する。
【００４８】
【数２】
2011173020.tif 000004
【００４９】
流速分散画像作成部２２ｄは、流速画像作成部２２ｃによって作成された複数の流速画像を用いて、各流速画像の同一位置ごとに体液の流速変化を示す統計値を算出し、算出した統計値の分布を表す統計値画像を作成する。本実施例では、流速分散画像作成部２２ｄは、統計値として、位相変化を利用して時系列に沿った流速成分の分散を算出し、統計値画像として、流速成分の分散の分布を表す流速分散画像を作成する。
【００５０】
具体的には、流速分散画像作成部２２ｄは、流速画像作成部２２ｃによって作成された複数の流速画像Ｖ（ｒａ，ｉ）と平均流速画像Ｖａｖｇ（ｒａ）とから、以下に示す式（５）で与えられる流速分散画像Ｖａｒｉ，Ｖ（ｒａ）を作成する。
【００５１】
【数３】
2011173020.tif 000005
【００５２】
重畳画像処理部２２ｅは、流速分散画像作成部２２ｄによって作成された流速分散画像における流速成分の分散の分布を表す重畳画像を作成する。具体的には、重畳画像処理部２２ｅは、流速分散画像Ｖａｒｉ，Ｖ（ｒａ）における流速成分の分散の分布を、絶対値画像作成部２２ｂによって作成された平均絶対値画像Ｉａｖｇ（ｒａ）に重畳させる。
【００５３】
このとき、重畳画像処理部２２ｅは、流速成分の分散値の値に応じて、流速成分の分散の分布の表示形態を変化させる。例えば、重畳画像処理部２２ｅは、平均絶対値画像Ｉａｖｇ（ｒａ）をグレースケールの濃淡画像とし、流速分散画像Ｖａｒｉ，Ｖ（ｒａ）における流速成分の分散の分布を以下に示す式（６）〜（８）で与えられるＲＧＢ（Red-Green-Blue）値で表す。なお、式（６）〜（８）において、Ｐ，ｒｅｄは赤（Ｒｅｄ）の輝度を、Ｐ，ｇｒｅｅｎは緑（Ｇｒｅｅｎ）の輝度を、Ｐ，ｂｌｕｅは青（Ｂｌｕｅ）の輝度を示している。また、ｋＩおよびｋＶは、それぞれ所定の比例係数である。
【００５４】
Ｐ，ｒｅｄ    ＝ｋＩ・Ｉａｖｇ（ｒａ）＋ｋＶ・Ｖａｒｉ，Ｖ（ｒａ）    ・・・（６）
Ｐ，ｇｒｅｅｎ＝ｋＩ・Ｉａｖｇ（ｒａ）    ・・・（７）
Ｐ，ｂｌｕｅ  ＝ｋＩ・Ｉａｖｇ（ｒａ）    ・・・（８）
【００５５】
これにより、重畳画像において流速成分の分散の分布が赤色で表されるとともに、分散値の値に応じて画素ごとに輝度が変化するので、操作者は体液の流速変化を容易に把握することができるようになる。なお、ここでは、重畳画像処理部２２ｅが、流速成分の分散の分布を赤色で表す場合について説明したが、例えば、緑色または青色で表すようにしてもよい。
【００５６】
図２にもどって、制御部２６は、特に、画像表示制御部２６ａを有する。画像表示制御部２６ａは、重畳画像処理部２２ｅによって作成された重畳画像を表示部２５に表示する。
【００５７】
次に、本実施例１に係るＭＲＩ装置１００による画像処理の処理手順について説明する。図５は、本実施例１に係るＭＲＩ装置１００による画像処理の処理手順を示すフローチャートであり、図６は、本実施例１に係るＭＲＩ装置１００による画像作成の流れを示す図である。
【００５８】
図５に示すように、このＭＲＩ装置１００では、制御部２６が、入力部２４を介して操作者からＥＰＩ法による撮像の撮像条件を受け付ける（ステップＳ１０１）。その後、操作者によって撮像開始が指示されると（ステップＳ１０２，Ｙｅｓ）、制御部２６が、ＥＰＩ法によるパルスシーケンスのシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御部１０に送信することによって、エコー信号に基づくｋ空間データを収集する（ステップＳ１０３）。
【００５９】
続いて、フーリエ変換処理部２２ａが、収集されたｋ空間データに離散２次元フーリエ変換などの再構成処理を施すことによって、図６の（１）に示すように、複数の画像Ｓ（ｒａ，ｉ）を再構成する（ステップＳ１０４）。
【００６０】
これに続いて、絶対値画像作成部２２ｂが、フーリエ変換処理部２２ａによって再構成された複数の画像Ｓ（ｒａ，ｉ）それぞれについて、図６の（２）に示すように、絶対値画像Ｉ（ｒａ，ｉ）を作成する（ステップＳ１０５）。さらに、絶対値画像作成部２２ｂは、作成した複数の絶対値画像Ｉ（ｒａ，ｉ）をもとに、図６の（３）に示すように、平均絶対値画像Ｉａｖｇ（ｒａ）を作成する（ステップＳ１０６）。
【００６１】
この一方で、流速画像作成部２２ｃが、フーリエ変換処理部２２ａによって再構成された複数の画像Ｓ（ｒａ，ｉ）それぞれについて、図６の（４）に示すように、流速画像Ｖ（ｒａ，ｉ）を作成する（ステップＳ１０７）。さらに、流速画像作成部２２ｃは、作成した複数の流速画像Ｖ（ｒａ，ｉ）をもとに、図６の（５）に示すように、平均流速画像Ｖａｖｇ（ｒａ）を作成する（ステップＳ１０８）。
【００６２】
これに続いて、流速分散画像作成部２２ｄが、流速画像作成部２２ｃによって作成された複数の流速画像Ｖ（ｒａ，ｉ）と平均流速画像Ｖａｖｇ（ｒａ）とから、図６の（６）に示すように、流速分散画像Ｖａｒｉ，Ｖ（ｒａ）を作成する（ステップＳ１０９）。
【００６３】
そして、重畳画像処理部２２ｅが、流速分散画像Ｖａｒｉ，Ｖ（ｒａ）の分散成分を色づけし（ステップＳ１１０）、図６の（７）に示すように、色づけした流速分散画像Ｖａｒｉ，Ｖ（ｒａ）を平均絶対値画像Ｉａｖｇ（ｒａ）に重畳させ、画像表示制御部２６ａが、重畳画像処理部２２ｅによって作成された重畳画像を表示部２５に表示する（ステップＳ１１１）。
【００６４】
上述してきたように、本実施例１では、流速画像作成部２２ｃが、被検体内を流れる体液の流速成分が得られるＥＰＩを複数回繰り返すことによって得られた複数の画像それぞれについて流速成分の分布を表す流速画像を作成する。また、流速分散画像作成部２２ｄが、作成された複数の流速画像を用いて、各流速画像の同一位置ごとに時系列に沿った流速成分の分散を算出し、算出した流速成分の分散の分布を表す流速分散画像を作成する。そして、重畳画像処理部２２ｅが、流速分散画像における流速成分の分散の分布を平均絶対値画像に重畳させ、画像表示制御部２６ａが、重畳画像を表示部２５に表示する。したがって、本実施例１によれば、ＣＳＦの流速変化のように流速変化に周期性がない体液であっても、その流速変化の分布を忠実に画像化することができる。
【００６５】
なお、本実施例１では、重畳画像処理部２２ｅが、平均絶対値画像Ｉａｖｇ（ｒａ）をグレースケールの濃淡画像として用いる場合について説明したが、例えば、平均絶対値画像Ｉａｖｇ（ｒａ）の替わりに他の撮像で得られた画像を用いてもよい。その場合、重畳画像処理部２２ｅは、例えば、Ｔ１強調画像やＴ２強調画像などを用いる。
【００６６】
また、本実施例１では、重畳画像処理部２２ｅが、流速成分の分散の分布を１つの色で表す場合について説明したが、さらに、体液が流れる方向に応じて、統計値の分布の表示態様を変化させるようにしてもよい。その場合、重畳画像処理部２２ｅは、例えば、超音波診断装置におけるパルスドップラーモードのように、血流の流れ方向を赤青の色彩情報で表し、血流の流速を明度で表し、血流の流速成分の分散を緑色で表すことによって、二次元血流画像を作成する。
【００６７】
また、本実施例１では、心電波形に同期させて撮像を繰り返すことによって得られた複数の画像を用いる場合について説明したが（図３を参照）、本発明はこれに限られるわけではない。図７は、本実施例１で用いられるパルスシーケンスの変形例を示す図である。例えば、図７に示すように、心電波形に同期させずに、一定の周期で撮像を複数回繰り返すことによって得られた複数の画像を用いてもよい。ＣＳＦの流速変化は心電同期との相関が低いため、心電波形に同期させずにエコー信号を収集しても支障は少ない。また、一定の周期で撮像を繰り返す場合には、同期撮像で必要な待ち時間（例えば、図３に示したＴｄｅｌａｙ）が不要となるので、より短時間で撮像を行うことができる。
【実施例２】
【００６８】
ところで、上記実施例１では、フローエンコード傾斜磁場を用いたＥＰＩ撮像を繰り返すことによって得られた複数の画像を用いて流速成分の分散の分布を画像化する場合について説明したが、本発明はこれに限られるわけではない。例えば、反転励起パルスを印加することによって撮像断面内にある一定の領域をラベリングし、当該領域から流れ出る体液の信号を検出して画像を再構成する撮像方法を用いてもよい（例えば、特開２００１−２５２２６３号公報を参照）。
【００６９】
そこで、以下では、かかる撮像方法を用いる場合を実施例２として説明する。なお、本実施例２では、この撮像方法を「Ｔｉｍｅ−ＳＬＩＰ法」と呼ぶ。また、本実施例２に係るＭＲＩ装置の全体構成は図１に示したものと同様であるので、ここでは説明を省略し、本実施例２における画像再構成部、記憶部および制御部の詳細について説明する。
【００７０】
図８は、本実施例２における画像再構成部、記憶部および制御部の構成を示す機能ブロック図である。なお、ここでは説明の便宜上、図２に示した各部と同様の役割を果たす機能部については、同一符号を付すこととしてその詳細な説明を省略する。
【００７１】
図８に示すように、本実施例２における画像再構成部１２２は、特に、フーリエ変換処理部１２２ａ、絶対値画像作成部１２２ｂ、信号強度分散画像作成部１２２ｄ、重畳画像処理部１２２ｅを有する。
【００７２】
フーリエ変換処理部１２２ａは、ｋ空間データ記憶部２３ａに記憶されたｋ空間データに離散２次元フーリエ変換などの再構成処理を施すことによって画像を再構成する。
【００７３】
ここで、本実施例２では、被検体内を流れる体液の流速成分が得られる撮像法として、Ｔｉｍｅ−ＳＬＩＰ法が用いられる。図９は、本実施例２で用いられるＴｉｍｅ−ＳＬＩＰ法のパルスシーケンスを示す図であり、図１０は、本実施例２におけるＴｉｍｅ−ＳＬＩＰ法によって行われるラベリングの位置の一例を示す図である。
【００７４】
図９に示すように、本実施例２で用いられるＴｉｍｅ−ＳＬＩＰ法では、心電波形（図９に示すＥＣＧ）のＲ波から一定の待ち時間（図９に示すＴＤｔａｇ１）が経過した後に、例えば、図１０に示すように、撮像断面内にある一定の領域を反転励起させるＩＲパルス（図９に示すＰ２１）が印加される。このとき、励起位置に合わせて励起周波数（図９に示すΔｆ）、選択励起の傾斜磁化（図９に示すＧｓｓ、Ｇｒｏ、Ｇｐｅ）を変化させる。
【００７５】
その後、一定時間（図９に示すＴＩ１）が経過した後に、１画像を再構成するのに必要なエコー信号を収集する（図９に示すＰ２２）。なお、図９は、高速スピンエコー法を用いてエコー信号を収集する場合を示している。本実施例２で用いられるＴｉｍｅ−ＳＬＩＰ法では、この収集を複数回繰り返すことによって複数画像分のエコー信号を収集する。
【００７６】
本実施例２では、フーリエ変換処理部１２２ａは、上記で説明したように収集されたエコー信号のｋ空間データそれぞれに離散２次元フーリエ変換などの再構成処理を施すことによって、複数の画像Ｓ（ｒａ，ｉ）を再構成する。
【００７７】
絶対値画像作成部１２２ｂは、フーリエ変換処理部１２２ａによって再構成された画像から信号強度の分布を表す絶対値画像を作成する。具体的には、絶対値画像作成部１２２ｂは、フーリエ変換処理部１２２ａによって再構成された複数の画像Ｓ（ｒａ，ｉ）それぞれについて、以下に示す式（９）で与えられる絶対値画像Ｉ（ｒａ，ｉ）を作成する。
【００７８】
Ｉ（ｒａ，ｉ）＝ａｂｓ｛Ｓ（ｒａ，ｉ）｝        ・・・（９）
【００７９】
図１１は、本実施例２における絶対値画像作成部１２２ｂによって作成される絶対値画像Ｉ（ｒａ，ｉ）の一例を示す図である。図１１に示すように、ラベリングされた領域から流れ出たＣＳＦが画像化されるが、各絶対値画像のもとになるエコー信号が収集された際のＣＳＦの流速がそれぞれ異なるため、各絶対値画像において、ラベリングされた領域から流れ出たＣＳＦによって信号強度が変化している領域は異なる。
【００８０】
さらに、絶対値画像作成部１２２ｂは、作成した複数の絶対値画像Ｉ（ｒａ，ｉ）をもとに、以下に示す式（１０）で与えられる平均絶対値画像Ｉａｖｇ（ｒａ）を作成する。
【００８１】
【数４】
2011173020.tif 000006
【００８２】
図８にもどって、信号強度分散画像作成部１２２ｄは、絶対値画像作成部１２２ｂによって作成された複数の絶対値画像を用いて、各絶対値画像の同一位置ごとに体液の流速変化を示す統計値を算出し、算出した統計値の分布を表す統計値画像を作成する。本実施例では、信号強度分散画像作成部１２２ｄは、統計値として、体液の信号値に基づいて時系列に沿った信号強度の分散を算出し、統計値画像として、信号強度の分散の分布を表す信号強度分散画像を作成する。
【００８３】
具体的には、信号強度分散画像作成部１２２ｄは、絶対値画像作成部１２２ｂによって作成された複数の絶対値画像Ｉ（ｒａ，ｉ）と平均絶対値画像Ｉａｖｇ（ｒａ）とから、以下に示す式（１１）で与えられる信号強度分散画像Ｖａｒｉ，Ｉ（ｒａ）を作成する。
【００８４】
【数５】
2011173020.tif 000007
【００８５】
重畳画像処理部１２２ｅは、信号強度分散画像作成部１２２ｄによって作成された信号強度分散画像における信号強度の分散の分布を表す重畳画像を作成する。具体的には、重畳画像処理部１２２ｅは、信号強度分散画像Ｖａｒｉ，Ｉ（ｒａ）における信号強度の分散の分布を、絶対値画像作成部１２２ｂによって作成された平均絶対値画像Ｉａｖｇ（ｒａ）に重畳させる。
【００８６】
このとき、重畳画像処理部１２２ｅは、信号強度の分散値の値に応じて、信号強度の分散の分布の表示形態を変化させる。例えば、重畳画像処理部１２２ｅは、平均絶対値画像Ｉａｖｇ（ｒａ）をグレースケールの濃淡画像とし、信号強度分散画像Ｖａｒｉ，Ｉ（ｒａ）における信号強度の分散の分布を以下に示す式（１２）〜（１４）で与えられるＲＧＢ（Red-Green-Blue）値で表す。なお、式（１２）〜（１４）において、Ｐ，ｒｅｄは赤（Ｒｅｄ）の輝度を、Ｐ，ｇｒｅｅｎは緑（Ｇｒｅｅｎ）の輝度を、Ｐ，ｂｌｕｅは青（Ｂｌｕｅ）の輝度を示している。また、ｋＩおよびｋＶは、それぞれ所定の比例係数である。
【００８７】
Ｐ，ｒｅｄ    ＝ｋＩ・Ｉａｖｇ（ｒａ）＋ｋＶ・Ｖａｒｉ，Ｉ（ｒａ）  ・・・（１２）
Ｐ，ｇｒｅｅｎ＝ｋＩ・Ｉａｖｇ（ｒａ）    ・・・（１３）
Ｐ，ｂｌｕｅ  ＝ｋＩ・Ｉａｖｇ（ｒａ）    ・・・（１４）
【００８８】
これにより、重畳画像において信号強度の分散の分布が赤色で表されるとともに、分散値の値に応じて画素ごとに輝度が変化するので、操作者は体液の流動による信号強度の変化を容易に把握することができるようになる。なお、ここでは、重畳画像処理部１２２ｅが、信号強度の分散の分布を赤色で表す場合について説明したが、例えば、緑色または青色で表すようにしてもよい。
【００８９】
次に、本実施例２に係るＭＲＩ装置１００による画像処理の処理手順について説明する。図１２は、本実施例２に係るＭＲＩ装置１００による画像処理の処理手順を示すフローチャートであり、図１３は、本実施例２に係るＭＲＩ装置１００による画像作成の流れを示す図である。
【００９０】
図１２に示すように、このＭＲＩ装置１００では、制御部２６が、入力部２４を介して操作者からＴｉｍｅ−ＳＬＩＰ法による撮像の撮像条件を受け付ける（ステップＳ２０１）。その後、操作者によって撮像開始が指示されると（ステップＳ２０２，Ｙｅｓ）、制御部２６が、Ｔｉｍｅ−ＳＬＩＰ法によるパルスシーケンスのシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御部１０に送信することによって、エコー信号に基づくｋ空間データを収集する（ステップＳ２０３）。
【００９１】
続いて、フーリエ変換処理部１２２ａが、収集されたｋ空間データに離散２次元フーリエ変換などの再構成処理を施すことによって、図１３の（１）に示すように、複数の画像Ｓ（ｒａ，ｉ）を再構成する（ステップＳ２０４）。
【００９２】
これに続いて、絶対値画像作成部１２２ｂが、フーリエ変換処理部１２２ａによって再構成された複数の画像Ｓ（ｒａ，ｉ）それぞれについて、図１３の（２）に示すように、絶対値画像Ｉ（ｒａ，ｉ）を作成する（ステップＳ２０５）。さらに、絶対値画像作成部１２２ｂは、作成した複数の絶対値画像Ｉ（ｒａ，ｉ）をもとに、図１３の（３）に示すように、平均絶対値画像Ｉａｖｇ（ｒａ）を作成する（ステップＳ２０６）。
【００９３】
その後、信号強度分散画像作成部１２２ｄが、絶対値画像作成部１２２ｂによって作成された複数の絶対値画像Ｉ（ｒａ，ｉ）と平均絶対値画像Ｉａｖｇ（ｒａ）とから、図１３の（４）に示すように、信号強度分散画像Ｖａｒｉ，Ｉ（ｒａ）を作成する（ステップＳ２０７）。
【００９４】
そして、重畳画像処理部１２２ｅが、信号強度分散画像Ｖａｒｉ，Ｉ（ｒａ）の分散成分を色づけし（ステップＳ２０８）、図１３の（５）に示すように、色づけした信号強度分散画像Ｖａｒｉ，Ｉ（ｒａ）を平均絶対値画像Ｉａｖｇ（ｒａ）に重畳させ、画像表示制御部２６ａが、重畳画像処理部１２２ｅによって作成された重畳画像を表示部２５に表示する（ステップＳ２０９）。
【００９５】
上述してきたように、本実施例２では、絶対値画像作成部１２２ｂが、ラベリングされた領域から流れ出る体液の信号を検出して画像を再構成するＴｉｍｅ―ＳＬＩＰ法によって得られた複数の画像それぞれについて信号強度の分布を表す絶対値画像を作成する。また、信号強度分散画像作成部１２２ｄが、作成された複数の絶対値画像を用いて、各絶対値画像の同一位置ごとに時系列に沿った信号強度の分散を算出し、算出した信号強度の分散の分布を表す信号強度分散画像を作成する。そして、重畳画像処理部１２２ｅが、信号強度分散画像における信号強度の分散の分布を平均絶対値画像に重畳させ、画像表示制御部２６ａが、重畳画像を表示部２５に表示する。したがって、本実施例２によれば、実施例１と同様、ＣＳＦの流速変化のように流速変化に周期性がない体液であっても、その流速変化の分布を忠実に画像化することができる。
【００９６】
なお、上記実施例１および２では、統計値として、流速成分あるいは信号強度の分散を算出する場合について説明したが、本発明はこれに限られるわけではない。例えば、統計値として、複数の計算画像上の同一位置における流速成分あるいは信号強度の最大値を算出し、つまり最大値画像Ｖｍａｘ（ｒａ）を作成し、算出した最大値に応じて流速成分を色づけするようにしてもよい。
【００９７】
または、統計値として、各計算画像における同一位置ごとに大きい順に流速成分あるいは信号強度を並べた場合に、全体に対して上位から所定の割合（例えば、上位１０％など）に入る流速成分あるいは信号強度の平均値を算出し、算出した平均値に応じて流速成分を色づけするようにしてもよい。この場合、具体的には、計算画像から同一ピクセルの画素値のヒストグラムを求め、全体に対して上位から所定の割合に入る画素値の平均を算出する。
【００９８】
なお、上記実施例１および２では、ＭＲＩ装置によって撮像された画像を用いる場合について説明したが、本発明はこれに限定されるものではなく、例えばＸ線ＣＴ装置やＸ線診断装置など、他の診断装置によって撮像された画像を用いる場合にも同様に適用することができる。すなわち、被検体内を流れる体液の流速成分が得られる撮像を複数回繰り返すことによって得られた複数の画像を用いる場合には、本発明を同様に適用することが可能である。
【００９９】
また、上記実施例１および２に係る画像処理方法は、十分な繰り返し回数の撮像が前提となるが、流速の周期的な変化を前提としないため、ＣＳＦのように心電同期とはあまり相関のない動きの画像上での分布を簡便に、かつ安定して得ることができる。また、１ショットの撮像を使用しているため、通常のフェイズシフト法のように平均加算効果によって流速変化が過小評価されるのを防ぐことができる。また、同期撮像を行わない場合には、不規則な動きに関する俯瞰的な情報を比較的簡便に得ることができる。
【符号の説明】
【０１００】
１    静磁場磁石
２    傾斜磁場コイル
３    傾斜磁場電源
４    寝台
４ａ    天板
５    寝台制御部
６    送信ＲＦコイル
７    送信部
８    受信ＲＦコイル
９    受信部
１０    シーケンス制御部
１５    架台部
２０    計算機システム
２１    インタフェース部
２２，１２２    画像再構成部
２２ａ，１２２ａ    フーリエ変換処理部
２２ｂ，１２２ｂ    絶対値画像作成部
２２ｃ    流速画像作成部
２２ｄ    流速分散画像作成部
１２２ｄ    信号強度分散画像作成部
２２ｅ，１２２ｅ    重畳画像処理部
２３    記憶部
２３ａ    ｋ空間データ記憶部
２３ｂ    画像データ記憶部
２４    入力部
２５    表示部
２６    制御部
２６ａ    画像表示制御部
１００    ＭＲＩ装置（磁気共鳴イメージング装置）
【図１】
2011173020.tif 000008
【図２】
2011173020.tif 000009
【図３】
2011173020.tif 000010
【図４】
2011173020.tif 000011
【図５】
2011173020.tif 000012
【図６】
2011173020.tif 000013
【図７】
2011173020.tif 000014
【図８】
2011173020.tif 000015
【図９】
2011173020.tif 000016
【図１０】
2011173020.tif 000017
【図１１】
2011173020.tif 000018
【図１２】
2011173020.tif 000019
【図１３】
2011173020.tif 000020
【図１４】
2011173020.tif 000021