特許 6786670 磁気共鳴イメージング装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6786670

(24)【登録日】2020年10月30日

(45)【発行日】2020年11月18日

(54)【発明の名称】磁気共鳴イメージング装置

(51)【国際特許分類】

   A61B 5/055 20060101AFI20201109BHJP

【ＦＩ】

   A61B5/055 382

   A61B5/055 376

   A61B5/055 311

【請求項の数】14

【全頁数】36

(21)【出願番号】特願2019-124132(P2019-124132)

(22)【出願日】2019年7月3日

(62)【分割の表示】特願2014-77341(P2014-77341)の分割

【原出願日】2014年4月3日

(65)【公開番号】特開2019-188205(P2019-188205A)

(43)【公開日】2019年10月31日

【審査請求日】2019年7月3日

(31)【優先権主張番号】特願2013-78873(P2013-78873)

(32)【優先日】2013年4月4日

(33)【優先権主張国】JP

(73)【特許権者】

【識別番号】594164542

【氏名又は名称】キヤノンメディカルシステムズ株式会社

(74)【代理人】

【識別番号】110001771

【氏名又は名称】特許業務法人虎ノ門知的財産事務所

(72)【発明者】

【氏名】木村 徳典

【審査官】 安田 明央

(56)【参考文献】

【文献】 特開２０００−０４１９６７（ＪＰ，Ａ）

【文献】 特開２００４−０２４６３７（ＪＰ，Ａ）

【文献】 国際公開第２０１２／０６０１９２（ＷＯ，Ａ１）

【文献】 米国特許出願公開第２０１３／０２０７６５３（ＵＳ，Ａ１）

【文献】 特開平１０−３２３３３６（ＪＰ，Ａ）

【文献】 特開２０１１−１３１０４１（ＪＰ，Ａ）

【文献】 米国特許出願公開第２０１１／０１３０６６８（ＵＳ，Ａ１）

【文献】 米国特許出願公開第２０１０／００３００６２（ＵＳ，Ａ１）

【文献】 特開２００８−０８６７４８（ＪＰ，Ａ）

【文献】 米国特許出願公開第２００８／００７１１６６（ＵＳ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ ５／０５５

Ｇ０１Ｒ ３３／２０−３３／５８

(57)【特許請求の範囲】

【請求項1】

被検体の撮像領域に流入する流体を標識化するためのＲＦ波を印加した後に、待ち時間の異なる複数のタイミングでデータ収集を行う場合に、前記待ち時間の長さに応じて、前記複数のタイミングを少なくとも２つの時間的な範囲に分割する分割部と、
前記撮像領域の全体について、第１の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が前記第１の読み出しシーケンスと異なる第２の読み出しシーケンスそれぞれを用いてデータ収集を行う収集部と、
前記第１の読み出しシーケンスによって収集されたデータから生成された画像と、前記第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成する合成部と
を備える、磁気共鳴イメージング装置。

【請求項2】

被検体の撮像領域に流入する流体を標識化するためのＲＦ波を印加して、待ち時間が経過した後にデータ収集を行うパターンを、前記待ち時間を異ならせて複数回実行する場合に、前記待ち時間の長さに応じて、複数のタイミングを少なくとも２つの時間的な範囲に分割する分割部と、
前記撮像領域の全体について、第１の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が前記第１の読み出しシーケンスと異なる第２の読み出しシーケンスそれぞれを用いてデータ収集を行う収集部と、
前記第１の読み出しシーケンスによって収集されたデータから生成された画像と、前記第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成する合成部と
を備える、磁気共鳴イメージング装置。

【請求項3】

被検体の撮像領域に流入する流体を標識化するためのＲＦ波を印加した後に、待ち時間の異なる複数のタイミングでデータ収集を行う場合に、前記待ち時間に基づいて、前記複数のタイミングを閾値により区切られた少なくとも２つの時間的な範囲に分割する分割部と、
前記撮像領域の全体について、第１の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が前記第１の読み出しシーケンスと異なる第２の読み出しシーケンスそれぞれを用いてデータ収集を行う収集部と、
前記第１の読み出しシーケンスによって収集されたデータから生成された画像と、前記第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成する合成部と
を備える、磁気共鳴イメージング装置。

【請求項4】

被検体の撮像領域に流入する流体を標識化するためのＲＦ波を印加して、待ち時間が経過した後にデータ収集を行うパターンを、前記待ち時間を異ならせて複数回実行する場合に、前記待ち時間に基づいて、前記待ち時間に応じた複数のタイミングを閾値により区切られた少なくとも２つの時間的な範囲に分割する分割部と、
前記撮像領域の全体について、第１の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が前記第１の読み出しシーケンスと異なる第２の読み出しシーケンスそれぞれを用いてデータ収集を行う収集部と、
前記第１の読み出しシーケンスによって収集されたデータから生成された画像と、前記第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成する合成部と
を備える、磁気共鳴イメージング装置。

【請求項5】

前記収集部は、前記時間的な範囲のうち、前記待ち時間が短い方の範囲については前記第１の読み出しシーケンスを用いてデータ収集を行い、前記待ち時間が長い方の範囲については前記第２の読み出しシーケンスを用いてデータ収集を行う、
請求項１〜４のいずれか一つに記載の磁気共鳴イメージング装置。

【請求項6】

前記収集部は、前記撮像領域における所定の位置に標識化された流体が到達したことを検出した場合に、データ収集用の読み出しシーケンスを前記第１の読み出しシーケンスから前記第２の読み出しシーケンスに変更する、
請求項１〜５のいずれか一つに記載の磁気共鳴イメージング装置。

【請求項7】

前記分割部は、前記流体として血液を標識化する撮像法で前記被検体の頭部を撮像する場合に、標識化された血液が主幹動脈を流れる範囲と、末梢動脈を流れる範囲とに前記撮像領域を分割する、
請求項１〜６のいずれか一つに記載の磁気共鳴イメージング装置。

【請求項8】

前記分割部は、所定の定量値を空間マッピングするための撮像法で前記被検体を撮像する場合に、対象部位の特性に応じて前記撮像領域を分割する、
請求項１〜６のいずれか一つに記載の磁気共鳴イメージング装置。

【請求項9】

前記分割部は、撮像対象の部位に関する解剖学的な情報、シミング用に測定された静磁場強度の分布、又は、撮像領域の位置決めに用いられる画像に基づいて、前記撮像領域を分割する、
請求項１〜６のいずれか一つに記載の磁気共鳴イメージング装置。

【請求項10】

前記分割部は、複数の待ち時間のうち少なくとも１つの待ち時間を閾値として設定することで、前記複数のタイミングを前記少なくとも２つの時間的な範囲に分割する、
請求項１〜６のいずれか一つに記載の磁気共鳴イメージング装置。

【請求項11】

前記分割部は、前記複数の待ち時間のうち、ＭＣＡ（Middle-Cerebral Artery）の左右に分かれる分岐を第１分枝とした場合の第２分枝に至る辺りに標識化された血液が到達するときの待ち時間を前記閾値として設定する、
請求項１０に記載の磁気共鳴イメージング装置。

【請求項12】

前記第２の読み出しシーケンスは、前記第１の読み出しシーケンスより単位収集時間あたりのＳＮ比が高い読み出しシーケンスである、
請求項１〜１１のいずれか一つに記載の磁気共鳴イメージング装置。

【請求項13】

前記第１の読み出しシーケンスは、前記第２の読み出しシーケンスより空間分解能が高い又は静磁場不均一性に強いシーケンスである、
請求項１〜１２のいずれか一つに記載の磁気共鳴イメージング装置。

【請求項14】

前記合成部は、前記第１の読み出しシーケンスによって収集されたデータから生成された画像と前記第２の読み出しシーケンスによって収集されたデータから生成された画像の少なくとも一方について、ゲイン補正、Ｔ１緩和補正及び画像補間のうち少なくとも１つの補正処理を行った後に、各画像を合成する、
請求項１〜１３のいずれか一つに記載の磁気共鳴イメージング装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明の実施形態は、磁気共鳴イメージング装置に関する。

【背景技術】

【0002】

従来、磁気共鳴イメージング装置で用いられるパルスシーケンスに含まれる読み出しシーケンスの種類には、大きく分類すると、ＳＥ（Spin Echo）系とＦＥ（Field Echo）（ＧＲＥ（Gradient Echo））系とがあり、さらに、ＦＥ系には、ＦＬＡＳＨ（Fast Low-Angle shot）系とＳＳＦＰ（Steady State Free Precession）系とがある。一般的に、これらの読み出しシーケンスは、前者ほど、単位収集時間あたりのＳＮ比（ＳＮＲ：Signal-to-Noise Ratio）は低いが磁場不均一性（歪みや信号抜け）に強く、後者ほど、単位収集時間あたりのＳＮ比は高いが磁場不均一性に弱いことが知られている。このため、磁気共鳴イメージング装置の操作者は、各読み出しシーケンスの特性を把握したうえで、撮像対象の部位や撮像の目的に応じて、読み出しシーケンスの種類や撮像条件を適宜に使い分けている。

【先行技術文献】

【非特許文献】

【0003】

【非特許文献1】Edelmann RR et al、Radiology 192:513-519(1994)

【非特許文献2】木村徳典: Modified STAR using asymmetric inversion slabs（ＡＳＴＡＲ）法による非侵襲血流イメージング、日磁医誌 2001;20(8)、374-385

【非特許文献3】Kwong KK, Chesler DA, koff RM, Donahue KM, et al、MR perfusion studies with T1-weighted echo planar imaging、Magn Reson Med 1995; 34:878-887

【非特許文献4】Mani S et al. MRM 37:898-905 (1997)

【非特許文献5】Wells JA at al. In vivo hadamard encoded continuous arterial spin labeling（Ｈ−ＣＡＳＬ）. MRM 63:1111-1118 (2010)

【発明の概要】

【発明が解決しようとする課題】

【0004】

本発明が解決しようとする課題は、撮像対象の部位や撮像の目的に応じて、より高画質な画像を得ることができる磁気共鳴イメージング装置を提供することである。

【課題を解決するための手段】

【0005】

実施形態に係る磁気共鳴イメージング（Magnetic Resonance Imaging：ＭＲＩ）装置は、分割部と、収集部と、合成部とを備える。分割部は、被検体の撮像領域に流入する流体を標識化するためのＲＦ波を印加した後に、待ち時間の異なる複数のタイミングでデータ収集を行う場合に、前記待ち時間に基づいて、前記複数のタイミングを少なくとも２つの時間的な範囲に分割する。収集部は、前記撮像領域の全体について、第１の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が前記第１の読み出しシーケンスと異なる第２の読み出しシーケンスそれぞれを用いてデータ収集を行う。合成部は、前記第１の読み出しシーケンスによって収集されたデータから生成された画像と、前記第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。

【図面の簡単な説明】

【0006】

【図1】図１は、第１の実施形態に係るＭＲＩ装置の構成例を示す図である。

【図2】図２は、従来のＡＳＬ法の一例を説明するための図（１）である。

【図3】図３は、従来のＡＳＬ法の一例を説明するための図（２）である。

【図4】図４は、従来のＡＳＬ法の一例を説明するための図（３）である。

【図5】図５は、従来の３ＤＴＯＦ−ＭＲＡを説明するための図（１）である。

【図6】図６は、従来の３ＤＴＯＦ−ＭＲＡを説明するための図（２）である。

【図7】図７は、第１の実施形態に係るＦＥ系の読み出しシーケンスの一例を示す図である。

【図8】図８は、第１の実施形態に係るＳＳＦＰ系の読み出しシーケンスの一例を示す図である。

【図9】図９は、第１の実施形態に係るＭＲＩ装置の詳細な構成例を示す図である。

【図10】図１０は、従来のＳＴ−ＭＩ−ＡＳＬによるデータ収集を示す図である。

【図11】図１１は、第１の実施形態に係る収集部によって行われるＳＴ−ＭＩ−ＡＳＬを示す図である。

【図12】図１２は、第１の実施形態に係る合成部によって行われる画像補正を説明するための図である。

【図13】図１３は、第１の実施形態に係るＭＲＩ装置によって行われる処理の流れを示すフローチャートである。

【図14】図１４は、第２の実施形態に係るＭＲＩ装置によって行われるマルチスラブ法を用いたＡＳＬ−ＭＲＡを説明するための図である。

【図15】図１５は、第２の実施形態に係るＭＲＩ装置によって行われるＳＴ−ＭＩ−ＡＳＬを用いたＡＳＬ−ＭＲＡを説明するための図である。

【図16】図１６は、第３の実施形態に係るＭＲＩ装置によって行われるＰＡＳＬ法を説明するための図（１）である。

【図17】図１７は、第３の実施形態に係るＭＲＩ装置によって行われるＰＡＳＬ法を説明するための図（２）である。

【図18】図１８は、第３の実施形態に係るＭＲＩ装置によって行われるＰＡＳＬ法を説明するための図（３）である。

【図19】図１９は、第３の実施形態に係るＭＲＩ装置によって行われるマルチスラブ法を用いたＰＡＳＬ法を説明するための図である。

【図20】図２０は、第４の実施形態に係るＭＲＩ装置によって行われるＨＥ−ＭＴ−ＰＡＳＬを説明するための図である。

【図21】図２１は、第４の実施形態に係るＭＲＩ装置によって行われるＳＴ−ＭＲ−ＡＳＬを説明するための図（１）である。

【図22】図２２は、第４の実施形態に係るＭＲＩ装置によって行われるＳＴ−ＭＲ−ＡＳＬを説明するための図（２）である。

【図23】図２３は、第５の実施形態に係るＭＲＩ装置によって行われるマルチスラブ３ＤＴＯＦ−ＭＲＡを説明するための図である。

【図24】図２４は、第６の実施形態に係るＭＲＩ装置によって行われるマルチスラブＡＳＬを説明するための図である。

【発明を実施するための形態】

【0007】

以下では添付図面を参照して、ＭＲＩ装置の実施形態を詳細に説明する。

【0008】

（第１の実施形態）
まず、第１の実施形態について説明する。図１は、第１の実施形態に係るＭＲＩ装置の構成例を示す図である。図１に示すように、ＭＲＩ装置１００は、静磁場磁石１、傾斜磁場コイル２、傾斜磁場電源３、寝台４、寝台制御部５、送信ＲＦコイル６、送信部７、受信ＲＦコイル８、受信部９、シーケンサ１０、ＥＣＧ（Electrocardiogram）センサ２１、ＥＣＧユニット２２、及び計算機システム３０を備える。

【0009】

静磁場磁石１は、中空の円筒形状に形成された磁石であり、内部の空間に一様な静磁場を発生する。この静磁場磁石１としては、例えば永久磁石、超伝導磁石等が使用される。

【0010】

傾斜磁場コイル２は、中空の円筒形状に形成されたコイルであり、静磁場磁石１の内側に配置される。この傾斜磁場コイル２は、互いに直交するＸ，Ｙ，Ｚの各軸に対応する３つのコイルが組み合わされて形成されており、これら３つのコイルは、後述する傾斜磁場電源３から個別に電流供給を受けて、Ｘ，Ｙ，Ｚの各軸に沿って磁場強度が変化する傾斜磁場を発生させる。なお、Ｚ軸方向は、静磁場と同方向とする。傾斜磁場電源３は、傾斜磁場コイル２に電流を供給する。

【0011】

ここで、傾斜磁場コイル２によって発生するＸ，Ｙ，Ｚ各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ及びリードアウト用傾斜磁場Ｇｒにそれぞれ対応する。スライス選択用傾斜磁場Ｇｓは、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場Ｇｅは、空間的位置に応じて磁気共鳴信号の位相を変化させるために利用される。リードアウト用傾斜磁場Ｇｒは、空間的位置に応じて磁気共鳴信号の周波数を変化させるために利用される。

【0012】

寝台４は、被検体Ｐが載置される天板４ａを備え、後述する寝台制御部５による制御のもと、被検体Ｐが載置された状態で天板４ａを傾斜磁場コイル２の空洞（撮像口）内へ挿入する。通常、この寝台４は、長手方向が静磁場磁石１の中心軸と平行になるように設置される。寝台制御部５は、制御部３６による制御のもと、寝台４を制御する装置であり、寝台４を駆動して、天板４ａを長手方向及び上下方向へ移動する。

【0013】

送信ＲＦコイル６は、傾斜磁場コイル２の内側に配置され、送信部７からＲＦ（Radio Frequency）パルスの供給を受けて、被検体ＰにＲＦ波を印加する。送信部７は、ラーモア周波数に対応するＲＦ波を送信ＲＦコイル６に送信する。

【0014】

受信ＲＦコイル８は、傾斜磁場コイル２の内側に配置され、ＲＦ波の影響によって被検体Ｐから放射される磁気共鳴信号を受信する。この受信ＲＦコイル８は、磁気共鳴信号を受信すると、その磁気共鳴信号を受信部９へ出力する。

【0015】

受信部９は、受信ＲＦコイル８から出力される磁気共鳴信号に基づいてＭＲ（Magnetic Resonance：磁気共鳴）データを生成する。具体的には、この受信部９は、受信ＲＦコイル８から出力される磁気共鳴信号をデジタル変換することによってＭＲデータを生成する。このＭＲデータは、前述したスライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ及びリードアウト用傾斜磁場Ｇｒによって、ＰＥ（Phase Encode）方向、ＲＯ（Read Out）方向、ＳＥ（Slice Encode）方向の空間周波数の情報が対応付けられることで、ｋ空間に対応するデータとして生成される。そして、ＭＲデータを生成すると、受信部９は、そのＭＲデータをシーケンサ１０へ送信する。なお、受信部９は、静磁場磁石１や傾斜磁場コイル２などを備える架台装置側に備えられていてもよい。

【0016】

シーケンサ１０は、計算機システム３０から送信されるシーケンス実行情報にしたがって傾斜磁場電源３、送信部７及び受信部９を駆動することによって、被検体Ｐの画像を生成するためのデータ収集を行う。ここで、シーケンス実行情報とは、被検体ＰからＭＲデータを収集するためのパルスシーケンスを定義する情報である。具体的には、シーケンス実行情報は、傾斜磁場電源３が傾斜磁場コイル２に供給する電源の強さや電源を供給するタイミング、送信部７が送信ＲＦコイル６に送信するＲＦ信号の強さやＲＦ信号を送信するタイミング、受信部９が磁気共鳴信号を検出するタイミングなど、データ収集を行うための手順を定義した情報である。

【0017】

なお、シーケンサ１０は、傾斜磁場電源３、送信部７及び受信部９を駆動して被検体Ｐをスキャンした結果、受信部９からＭＲデータが送信されると、そのＭＲデータを計算機システム３０へ転送する。

【0018】

ＥＣＧセンサ２１は、被検体Ｐの体表に付着され、被検体Ｐの心拍、脈波、呼吸などを示すＥＣＧ信号を電気信号として検出する。ＥＣＧユニット２２は、ＥＣＧセンサ２１により検出されたＥＣＧ信号にＡ／Ｄ変換処理やディレー処理を含む各種処理を施してゲート信号を生成し、生成したゲート信号をシーケンサ１０に送信する。

【0019】

計算機システム３０は、ＭＲＩ装置１００の全体制御を行う。例えば、計算機システム３０は、上述した各部を駆動することで、データ収集や画像再構成などを行う。この計算機システム３０は、インタフェース部３１、画像再構成部３２、記憶部３３、入力部３４、表示部３５、及び制御部３６を有する。

【0020】

インタフェース部３１は、計算機システム３０とシーケンサ１０との間でやり取りされる各種信号の送受信を制御する。例えば、このインタフェース部３１は、シーケンサ１０に対してシーケンス実行情報を送信し、シーケンサ１０からＭＲデータを受信する。ＭＲデータを受信すると、インタフェース部３１は、各ＭＲデータを被検体Ｐごとに記憶部３３に格納する。

【0021】

画像再構成部３２は、記憶部３３によって記憶されたＭＲデータに対して後処理すなわちフーリエ変換等の再構成処理を施すことで、被検体Ｐの体内が描出された画像データを生成する。

【0022】

記憶部３３は、インタフェース部３１により受信されたＭＲデータや、画像再構成部３２により生成された画像データなどを被検体Ｐごとに記憶する。

【0023】

入力部３４は、操作者からの各種指示や情報入力を受け付ける。この入力部３４としては、マウスやトラックボールなどのポインティングデバイス、モード切替スイッチ等の選択デバイス、あるいはキーボード等の入力デバイスを適宜に利用可能である。

【0024】

表示部３５は、制御部３６による制御のもと、スペクトラムデータあるいは画像データ等の各種の情報を表示する。この表示部３５としては、液晶表示器などの表示デバイスを利用可能である。

【0025】

制御部３６は、図示していないＣＰＵ（Central Processing Unit）やメモリ等を有し、ＭＲＩ装置１００の全体制御を行う。具体的には、この制御部３６は、入力部３４を介して操作者から受け付けられた各種指示に基づいてシーケンス実行情報を生成し、生成したシーケンス実行情報をシーケンサ１０に送信することによってスキャンを制御したり、スキャンの結果としてシーケンサ１０から送られるＭＲデータに基づいて行われる画像の再構成を制御したりする。

【0026】

以上、第１の実施形態に係るＭＲＩ装置１００の構成例について説明した。このような構成のもと、ＭＲＩ装置１００は、各種の読み出しシーケンスを用いて、様々な撮像法で被検体を撮像することが可能である。なお、ここでは、励起用のＲＦパルスを印加してから磁気共鳴信号を取得するまでのシーケンスを「読み出しシーケンス」と呼び、読み出しシーケンスの前に印加される標識化用のＲＦパルスなどの各種プリパルスと読み出しシーケンスとを含んだ全体的なシーケンスを「パルスシーケンス」と呼ぶ。

【0027】

ここで、従来、ＭＲＩ装置で用いられるパルスシーケンスに含まれる読み出しシーケンスの種類には、大きく分類すると、ＳＥ系とＦＥ（ＧＲＥ）系とがあり、さらに、ＦＥ系には、ＦＬＡＳＨ系とＳＳＦＰ系とがある。一般的に、これらの読み出しシーケンスは、前者ほど、単位収集時間あたりのＳＮ比は低いが磁場不均一性（歪みや信号抜け）に強く、後者ほど、単位収集時間あたりのＳＮ比は高いが磁場均一性に弱いことが知られている。このため、磁気共鳴イメージング装置の操作者は、各読み出しシーケンスの特性を把握したうえで、撮像対象の部位や撮像の目的に応じて、読み出しシーケンスの種類や撮像条件を適宜に使い分けている。

【0028】

例えば、ＭＲＩ装置による撮像法の一つに、ＡＳＬ（Arterial Spin Labeling）法がある。ＡＳＬ法は、撮像領域に流入する流体を標識化するためのＲＦ波を印加してから所定の待ち時間が経過した後にデータ収集を行うことによって、非造影で血流を画像化する撮像法である。

【0029】

図２〜４は、従来のＡＳＬ法を説明するための図である。一般的に、ＡＳＬ法では、背景信号を抑制するために、タグ画像と呼ばれる画像と、コントロール画像と呼ばれる画像とがそれぞれ生成され、両画像の差分画像が生成される。

【0030】

例えば、図２に示すように、ＡＳＬ法では、撮像領域（imaging slab）の上流に設定されたタグ領域（tag slab）に、撮像領域に流入する血液（arterial blood flow）を標識化するためのＲＦ波が印加される。そして、図３に示すように、ＲＦ波（tag）が印加されてから所定の待ち時間ＴＩが経過した後にデータ収集（imaging）が行われ、収集されたデータからタグ画像が生成される。さらに、ＡＳＬ法では、図２に示すように、撮像領域を挟んでタグ領域の反対側に設定されたコントロール領域（control slab）にもＲＦ波が印加される。そして、図３に示すように、ＲＦ波（control）が印加されてから、タグ画像用のデータ収集が行われたときと同じ待ち時間ＴＩが経過した後にデータ収集（imaging）が行われ、収集されたデータからコントロール画像が生成される。その後、図４に示すように、生成されたタグ画像（tag image）とコントロール画像（control image）との差分画像を生成することで、背景信号を抑制したＡＳＬ画像（ASL image）が得られる。

【0031】

また、ＡＳＬ法は、その用途から、大きく、ＡＳＬ−ＭＲＡ（MR Angio）とＡＳＬ−ＭＲＰ（MR Perfusion）とに分けられる。ＡＳＬ−ＭＲＡは、血管画像を得るための方法であり、ＡＳＬ−ＭＲＰは、灌流画像を得るための方法である。さらに、近年では、ＡＳＬ−ＭＲＡ又はＡＳＬ−ＭＲＰの時間的な変化（hemodynamics）を観察するために、標識化用のＲＦ波を印加してからデータ収集までの待ち時間を動的に変えて、複数回撮像を行うことも行われている。

【0032】

ここで、ＡＳＬ−ＭＲＡとＡＳＬ−ＭＲＰとは、どちらも撮像法の種類はＡＳＬ法であるが、技術的なポイントが異なる。例えば、ＡＳＬ−ＭＲＡでは、血管のＳＮ比が大きいため加算回数は少なくてよいが、時間分解能及び空間分解能が要求される。また、ＡＳＬ−ＭＲＰでは、ＡＳＬ−ＭＲＡに比べて血流信号のＳＮ比が背景の静止組織よりもはるかに低い（１／１０００以下）画像が得られるため、タグ画像とコントロール画像との差分を行う場合に、それぞれの画像でＭＴＣ（Magnetization Transfer Contrast）効果を同じにしたり、低マトリクスでの複数回の加算平均が行われたりする。さらに、ＡＳＬ−ＭＲＰの中でも、血流量は腎や肺などは高く、脳や前立腺、筋などは低いため、撮像対象の臓器によって読み出しシーケンスの撮像条件が異なる。

【0033】

一方、一般的に最も広く利用されている３ＤＴＯＦ（3 Dimensional Time-Of-Flight）−ＭＲＡでは、血管流のタイムオブフライト（time-of-flight）効果を用いるため、１回の撮像当たりの撮像領域（imaging slab）の厚さは薄いほどＳＮ比が高くなる。そのため、３ＤＴＯＦ−ＭＲＡでは、撮像領域を複数の分割領域に分割して当該複数の分割領域ごとにデータ収集を行うマルチスラブ（multi-slab）法が用いられる。

【0034】

図５及び６は、従来の３ＤＴＯＦ−ＭＲＡを説明するための図である。図５は、シングルスラブ（single-slab）法による撮像の一例を示しており、図６は、マルチスラブ（multi-slab）法による撮像の一例を示している。

【0035】

例えば、図５に示すように、シングルスラブ法では、撮像領域（imaging slab）に照射するＲＦ波に傾斜をつける（流入部から距離的に遠いほどＲＦ波のフリップ角（Flip Angle）を浅くする）ＴＯＮＥ（Tilted Optimized Non-saturation Excitation）法（ＩＳＣＥ（Inclined Slab for Contrast Enhancement）法）が併用される場合が多い。また、図６に示すように、マルチスラブ法でＴＯＮＥ法が併用される場合もある。このような３Ｄ−ＭＲＡで頭部を撮像する場合は、ウィリス動脈輪（Willis ring）部の近傍にある主幹動脈部の抜け（void）が少ない画像が得られるＦＥ（ＧＲＥ）系の読み出しシーケンスが用いられるのが一般的である。

【0036】

また、ＭＲＩ装置による他の撮像法の一つに、拡散イメージング（diffusion imaging）がある。この拡散イメージングでは、シングルショットのＳＥ−ＥＰＩ（Echo Planar Imaging）系の読み出しシーケンスが主に用いられるが、歪みの大きな部位についてはＦＳＥ（Fast Spin Echo）系の読み出しシーケンスも用いられている。また、ＥＰＩの中でも、ＴＥ（Echo Time）、収集マトリクス（matrix）、ＳＥＮＳＥ（Sensitivity Encoding）−Ｆａｃｔｏｒ、ショット数が使い分けられている。さらには、ｋ空間の軌跡（trajectory）の違いであるＰＲＯＰＥＬＬＥＲ（Periodically Rotated Overlapping Parallel Lines with Enhanced Reconstruction）、ｒａｄｉａｌ系などのサンプリング（sampling）法も使い分けられている。

【0037】

さらに、最近の傾向として、Ｔ１やＴ２、プロトン（proton）密度、ＡＤＣ（Apparent Diffusion Coefficient）、血流量、ケミカルシフト、磁化率、温度などのＭＲＩに関する各種パラメータを強調画像ではなく定量値画像（パラメータ画像）に変換する手法が盛んに用いられている。これらの各種パラメータの定量値画像を得るための撮像法では、様々な読み出しシーケンスが利用可能になっており、ＳＮ比や歪み、時間などの特性によって、使用する読み出しシーケンスの種類や撮像条件が決められるようになっている。

【0038】

このように、ＭＲＩ装置では、各種の撮像法で被検体を撮像することが可能であり、各撮像法において、撮像対象の部位や撮像の目的に応じて読み出しシーケンスの種類や撮像が使い分けられている。

【0039】

例えば、前述したＡＳＬ法では、エコー信号を収集するための読み出しシーケンス（imaging sequence）として、ＦＥ（ＧＲＥ）、ＳＳＦＰ、ＦＥ−ＥＰＩなどのＦＥ（ＧＲＥ）系の読み出しシーケンスが用いられ、さらには、ＦＳＥ、ＳＥ−ＥＰＩなどのＳＥ系の読み出しシーケンスも用いられる。また、ｋ空間を充填する方法も、Ｃａｒｔｅｓｉａｎ、ｎｏｎ−Ｃａｒｔｅｓｉａｎ（ｓｐｉｒａｌ、ｒａｄｉａｌ）などの様々な方法が用いられる。これらの中には、対象部位（対象臓器）や使用目的の全ての場合に最適な方法、すなわちオールマイティな方法はなく、ＭＲＡかＭＲＰかに応じて、又は、撮像対象の臓器に応じて、適宜使い分けられている。

【0040】

ここで、前述したように、ＦＥ（ＧＲＥ）系はＳＥ系と比べてＳＮ比は高いが磁場不均一性に弱く、さらにＦＥ系の中でも、ＳＳＦＰ（Fast Imaging with Steady-state Precession：ＦＩＳＰとも呼ばれる）系はＦＬＡＳＨ系と比べてＳＮ比は高いが、さらに磁場不均一性に弱い。そのため、ＡＳＬ−ＭＲＡで脳の主幹部を撮像する場合は、副鼻腔や側頭骨などがあるためＦＥ系が主に用いられ、腹部の臓器を撮像する場合は、ＳＳＦＰ系が用いられることが多い。また、脳の中でも、頭頂部の末梢動脈を撮像する場合は磁場不均一性がさほど問題にならないため、ＳＮ比が高いＳＳＦＰ系が優れていることが知られている。

【0041】

また、Ｔｉｍｅ−Ｒｅｓｏｌｖｅ ＭＲＡ（ｔｒＭＲＡ）では、標識化された血液は、ＴＩが短いときは主幹動脈にあり、ＴＩが長くなると末梢動脈に移行していく。そのため、ＴＩが短いときはＦＥ（ＧＲＥ）系が、ＴＩが長いときはＳＳＦＰ系が適していることになる。また、３ＤＴＯＦ−ＭＲＡでは、ウィリス動脈輪（Willis ring）部の近傍はＦＥ系で、頭頂部はＳＳＦＰ系が適していることになる。

【0042】

このように、ＭＲＩ装置による撮像法では、被検体を時系列に撮像する際の時間的な位置や撮像領域の空間的な位置に応じて、それぞれの位置に適した読み出しシーケンスの種類や撮像条件を使い分けるのが望ましい。しかしながら、一般的に、操作者が時間的な位置や空間的な位置に応じて読み出しシーケンスの種類や撮像条件を切り替えることは難しく、装置側で制御する仕組みが求められていた。

【0043】

上述した従来技術の課題に対し、第１の実施形態に係るＭＲＩ装置１００は、被検体を撮像する際の時間的な範囲又は空間的な範囲を少なくとも２つの範囲に分割し、第１の範囲については第１の読み出しシーケンスを用いてデータ収集を行い、第２の範囲についてはシーケンスの種類又は撮像条件が第１の読み出しシーケンスと異なる第２の読み出しシーケンスを用いてデータ収集を行う。そして、ＭＲＩ装置１００は、第１の読み出しシーケンスによって収集されたデータから生成された画像と第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。

【0044】

このような構成によれば、第１の読み出しシーケンスの種類又は撮像条件を第１の範囲を撮像するために適したものとし、第２の読み出しシーケンスの種類又は撮像条件を第２の範囲を撮像するために適したものとすることで、被検体を撮像する際の時間的な位置又は空間的な位置に応じて、適切な読み出しシーケンスを用いて撮像が行われることになる。すなわち、第１の実施形態に係るＭＲＩ装置１００によれば、被検体を撮像する際の時間的な位置又は空間的な位置に応じて、各位置に適した読み出しシーケンスの種類又は撮像条件を用いて撮像を行うことで、より高画質な画像を容易に得ることができる。

【0045】

以下では、第１の実施形態に係るＭＲＩ装置１００の詳細な構成、及び、第１の実施形態に係るＭＲＩ装置１００によって行われる処理の流れについて説明する。なお、本実施形態では、被検体の撮像領域に流入する流体を標識化するための１回のＲＦ波を印加してから複数の待ち時間が経過するごとに連続的にデータ収集を行うＳＴ−ＭＩ（Single-Tag Multiple-TI）−ＡＳＬで被検体の頭部血管を撮像する場合の例を説明する。また、本実施形態では、第１の読み出しシーケンスとして、一般的なＦＥ系（以下、単にＦＥ系と呼ぶ）の読み出しシーケンスを用い、第２の読み出しシーケンスとして、ＳＳＦＰ系の読み出しシーケンスを用いる場合の例を説明する。

【0046】

図７は、第１の実施形態に係るＦＥ系の読み出しシーケンスの一例を示す図である。また、図８は、第１の実施形態に係るＳＳＦＰ系の読み出しシーケンスの一例を示す図である。なお、図７及び８に例示する各読み出しシーケンスは、いずれも３Ｄ撮像用の読み出しシーケンスである。具体的には、図７に例示するＦＥ系の読み出しシーケンスは、リードアウト用傾斜磁場Ｇｒを反転させることでエコー信号を発生させるものであり、さらに、位相エンコード用傾斜磁場Ｇｅについて、印加した傾斜磁場をデータ収集後にキャンセルするためのリワインド用傾斜磁場を印加するものである。また、図８に例示するＳＳＦＰ系の読み出しシーケンスは、図７に示す読み出しシーケンスにおいて、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ及びリードアウト用傾斜磁場Ｇｒの全てについて、繰り返し時間ＴＲの中で傾斜磁場の積分値がゼロになるようにリワインド用傾斜磁場を印加するものである。

【0047】

ここで、図８に例示するＳＳＦＰ系の読み出しシーケンスは、図７に例示するＦＥ系の読み出しシーケンスと比べると、単位収集時間あたりのＳＮ比は高いが、３軸方向にリワインド用傾斜磁場を印加して傾斜磁場の積分値を調整するため、磁場不均一性に弱い。一方、図７に例示するＦＥ系の読み出しシーケンスは、図８に例示するＳＳＦＰ系の読み出しシーケンスと比べると、磁場不均一性に強く、かつ、空間分解能が高いが、単位収集時間あたりのＳＮ比は低い。

【0048】

そこで、本実施形態では、ＭＲＩ装置１００は、ＳＴ−ＭＩ−ＡＳＬで被検体頭部血管を撮像する場合に、標識化用のＲＦ波を印加してからデータ収集を行うまでの待ち時間に基づいて、撮像時の時間的な範囲を少なくとも２つの範囲に分割する。そして、ＭＲＩ装置１００は、分割した２つの範囲のうち、待ち時間が短い方の範囲については、磁場不均一性に強いＦＥ系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲については、ＳＮ比が高いＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行う。

【0049】

図９は、第１の実施形態に係るＭＲＩ装置１００の詳細な構成例を示す図である。図９は、図１に示したシーケンサ１０及び計算機システム３０を示している。また、図９は、計算機システム３０が有する機能部のうち、インタフェース部３１、画像再構成部３２、記憶部３３、入力部３４、表示部３５、及び制御部３６を示している。

【0050】

図９に示すように、記憶部３３は、撮像条件記憶部３３ａと、分割情報記憶部３３ｂと、ＭＲデータ記憶部３３ｃと、画像データ記憶部３３ｄとを有する。

【0051】

撮像条件記憶部３３ａは、被検体の撮像に用いられる読み出しシーケンスの種類ごとに撮像条件を記憶する。なお、ここでいう撮像条件とは、読み出しシーケンスの生成に必要な各種撮像パラメータのパラメータ値であり、例えば、収集マトリクス、ＮＡＱ、スラブ厚、１スラブ内のスライス枚数、ショット数、ＳＥＮＳＥ−Ｆａｃｔｏｒなどである。また、撮像条件には、ＰＲＯＰＥＬＬＥＲ、ｒａｄｉａｌ系などのサンプリング法の種類も含まれる。本実施形態では、撮像条件記憶部３３ａは、少なくとも、ＦＥ系の読み出しシーケンスに関する撮像条件と、ＳＳＦＰ系の読み出しシーケンスに関する撮像条件とを記憶する。

【0052】

分割情報記憶部３３ｂは、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割するための基準となる情報を記憶する。例えば、分割情報記憶部３３ｂは、撮像対象の部位に関する解剖学的な情報を記憶する。本実施形態では、分割情報記憶部３３ｂは、ＳＴ−ＭＩ−ＡＳＬで被検体の頭部血管を撮像する場合に、標識化された血液が主幹動脈を流れる範囲と末梢動脈を流れる範囲とに時間的な範囲を分割するための情報として、ＭＣＡ（Middle-Cerebral Artery）の左右に分かれる分岐を第１分枝とした場合の第２分枝に至る辺りに標識化された血液が到達するときの待ち時間を記憶する。

【0053】

なお、分割情報記憶部３３ｂが記憶する情報は、解剖学的な情報に限られない。例えば、過去に行われた撮像において基準として用いられた待ち時間のうち、操作者が最も適していると判断した待ち時間を記憶しておいてもよい。すなわち、待ち時間に関する操作者の経験値を記憶しておいてもよい。

【0054】

ＭＲデータ記憶部３３ｃは、インタフェース部３１を介してシーケンサ１０から受信されたＭＲデータを記憶する。画像データ記憶部３３ｄは、画像再構成部３２によってＭＲデータから再構成された画像を記憶する。

【0055】

また、制御部３６は、撮像条件設定部３６ａと、分割部３６ｂと、収集部３６ｃと、血流画像生成部３６ｄと、合成部３６ｅとを有する。

【0056】

撮像条件設定部３６ａは、被検体の撮像に用いられる読み出しシーケンスの種類ごとに撮像条件を設定する。具体的には、撮像条件設定部３６ａは、入力部３４を介して、操作者から各種撮像パラメータのパラメータ値の入力を受け付け、受け付けたパラメータ値を撮像条件として撮像条件記憶部３３ａに登録する。本実施形態では、撮像条件設定部３６ａは、少なくとも、ＦＥ系の読み出しシーケンスに関する撮像条件と、ＳＳＦＰ系の読み出しシーケンスに関する撮像条件とを設定する。

【0057】

分割部３６ｂは、被検体を撮像する際の時間的な範囲又は空間的な範囲を少なくとも２つの範囲に分割する。具体的には、分割部３６ｂは、分割情報記憶部３３ｂによって記憶された情報に基づいて、時間的な範囲又は空間的な範囲を分割する。本実施形態では、分割部３６ｂは、ＳＴ−ＭＩ−ＡＳＬで被検体の頭部血管を撮像する場合に、標識化用のＲＦ波を印加してからデータ収集を行うまでの待ち時間に基づいて、撮像時の時間的な範囲を少なくとも２つの範囲に分割する。

【0058】

例えば、分割部３６ｂは、複数の待ち時間のうち少なくとも１つの待ち時間を閾値として設定することで、撮像時の時間的な範囲を少なくとも２つの範囲に分割する。具体的には、分割部３６ｂは、ＳＴ−ＭＩ−ＡＳＬで用いられる複数の待ち時間のうち、分割情報記憶部３３ｂによって記憶されている待ち時間と一致する待ち時間を閾値として設定する。本実施形態では、分割部３６ｂは、ＭＣＡの左右に分かれる分岐を第１分枝とした場合の第２分枝に至る辺りに標識化された血液が到達するときの待ち時間を閾値として設定することで、標識化された血液が主幹動脈を流れる範囲と末梢動脈を流れる範囲とに時間的な範囲を分割する。

【0059】

収集部３６ｃは、分割部３６ｂによって分割された少なくとも２つの範囲のうち、第１の範囲については第１の読み出しシーケンスを用いてデータ収集を行い、第２の範囲については読み出しシーケンスの種類又は撮像条件が第１の読み出しシーケンスと異なる第２の読み出しシーケンスを用いてデータ収集を行う。

【0060】

本実施形態では、収集部３６ｃは、ＳＴ−ＭＩ−ＡＳＬで被検体の頭部血管を撮像する場合に、分割された２つの範囲のうち、待ち時間が短い方の範囲についてはＦＥ系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲についてはＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行う。具体的には、収集部３６ｃは、撮像条件記憶部３３ａによって記憶されている撮像条件に基づいて、ＦＥ系のシーケンス情報及びＳＳＦＰ系のシーケンス情報をそれぞれ生成する。そして、収集部３６ｃは、生成した各シーケンス情報をシーケンサ１０に送信し、それらの各シーケンス情報を用いてＳＴ−ＭＩ−ＡＳＬを実行するようにシーケンサ１０を制御する。

【0061】

図１０は、従来のＳＴ−ＭＩ−ＡＳＬによるデータ収集を示す図である。また、図１１は、第１の実施形態に係る収集部３６ｃによって行われるＳＴ−ＭＩ−ＡＳＬを示す図である。図１０に示すように、従来のＳＴ−ＭＩ−ＡＳＬでは、複数の待ち時間ＴＩ１〜ＴＩ６ごとに、同じ種類の読み出しシーケンス（例えば、ＦＥ系の読み出しシーケンス）を用いてデータ収集が行われていた（Ｓ１〜Ｓ６）。

【0062】

これに対し、本実施形態では、収集部３６ｃは、分割部３６ｂによって設定された待ち時間の閾値に基づいて、データ収集に用いる読み出しシーケンスを切り替える。例えば、待ち時間の閾値がＴＩ３であった場合には、収集部３６ｃは、図１１に示すように、待ち時間がＴＩ３以下であるときは、ＦＥ系の読み出しシーケンスを用いてデータ収集を行い、待ち時間がＴＩ３を超えているときは、ＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行うように、シーケンサ１０を制御する。

【0063】

ここで、本実施形態では、分割部３６ｂによって、ＭＣＡ（Middle-Cerebral Artery）の左右に分かれる分岐を第１分枝とした場合の第２分枝に至る辺りに標識化された血液が到達するときの待ち時間が閾値として設定される。これにより、本実施形態では、標識化された血液が主幹動脈を流れる範囲については、ＦＥ系の読み出しシーケンスでデータ収集が行われ、標識化された血液が末梢動脈を流れる範囲については、ＳＳＦＰ系の読み出しシーケンスでデータ収集を行われることになる。

【0064】

なお、収集部３６ｃは、ＦＥ系の読み出しシーケンス及びＳＳＦＰ系の読み出しシーケンスのどちらを用いてデータ収集を行う場合も、各待ち時間において、タグ画像用のデータ収集とコントロール画像用のデータ収集とをそれぞれ行うようにシーケンサ１０を制御する。そして、収集部３６ｃは、インタフェース部３１を介してシーケンサ１０からＭＲデータを受信すると、受信したＭＲデータをＭＲデータ記憶部３３ｃに格納する。この結果、画像再構成部３２によって、ＳＴ−ＭＩ−ＡＳＬで用いられる複数の待ち時間ごとに、シーケンサ１０によって収集されたＭＲデータからタグ画像及びコントロール画像がそれぞれ生成されることになる。

【0065】

図９に戻って、血流画像生成部３６ｄは、第１の読み出しシーケンスによって収集されたデータ及び第２の読み出しシーケンスによって収集されたデータそれぞれから血流画像を生成する。本実施形態では、血流画像生成部３６ｄは、ＳＴ−ＭＩ−ＡＳＬで用いられる複数の待ち時間ごとに、画像再構成部３２によって生成されたタグ画像及びコントロール画像の差分画像を血流画像として生成する。すなわち、血流画像生成部３６ｄは、ＦＥ系の読み出しシーケンスによって収集されたデータ及びＳＳＦＰ系の読み出しシーケンスによって収集されたデータそれぞれについて、待ち時間ごとに血流画像を生成する。

【0066】

合成部３６ｅは、第１の読み出しシーケンスによって収集されたデータから生成された画像と第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。本実施形態では、合成部３６ｅは、ＳＴ−ＭＩ−ＡＳＬが実行される場合に、血流画像生成部３６ｄによって待ち時間ごとに生成された各血流画像を合成する。すなわち、合成部３６ｅは、ＦＥ系の読み出しシーケンスによって収集されたデータから生成された血流画像と、ＳＳＦＰ系の読み出しシーケンスによって収集されたデータから生成された血流画像とを合成する。そして、合成部３６ｅは、合成した血流画像を表示部３５に表示させる。

【0067】

ここで、合成部３６ｅは、第１の読み出しシーケンスによって収集されたデータから生成された血流画像と第２の読み出しシーケンスによって収集されたデータから生成された血流画像の少なくとも一方について、ゲイン補正、Ｔ１緩和補正及び画像補間のうち少なくとも１つの補正処理を行った後に、各画像を合成してもよい。本実施形態では、ＦＥ系の読み出しシーケンスを用いて得られた血流画像及びＳＳＦＰ系の読み出しシーケンスを用いて得られた血流画像の少なくとも一方について補正処理を行う。

【0068】

図１２は、第１の実施形態に係る合成部３６ｅによって行われる画像補正を説明するための図である。図１２の上段、中段及び下段に示す図において、横軸はＳＴ−ＭＩ−ＡＳＬで用いられる待ち時間ＴＩを示しており、縦軸は血流画像における血流部分の信号強度Ｓ．Ｉ．（Signal Intensity）を示している。

【0069】

例えば、図１２の上段に示すように、０＜ＴＩ＜ＴＩ_thの範囲ではＦＥ系の読み出しシーケンスで、ＴＩ＝ＴＩ_thではＦＥ系及びＳＳＦＰ系の両方の読み出しシーケンスで、ＴＩ_th＜ＴＩ＜ＴＩ_maxの範囲ではＳＳＰＦ系の読み出しシーケンスでデータ収集を行うとする。この場合には、ＦＥ系及びＳＳＦＰ系それぞれの読み出しシーケンスでデータ収集を行う際の受信ゲイン及びノイズＳＤ（Noise S.D.）値を同じにしたとすると、血流部分の信号強度は、磁場不均一性が問題とならない範囲（ＴＩが長い範囲）では、ＦＥ系の読み出しシーケンスで得られた血流画像よりも、ＳＳＦＰ系の読み出しシーケンスで得られた血流画像のほうが大きい。なお、ここでいうノイズＳＤ値は、画像に含まれる画素ごとの信号強度の標準偏差（Standard Deviation）であり、ノイズの特性を示す指標値である。

【0070】

そこで、合成部３６ｅは、ＦＥ系の読み出しシーケンスを用いて得られた血流画像及びＳＳＦＰ系の読み出しシーケンスを用いて得られた血流画像の少なくとも一方について、ゲイン補正を行う。例えば、図１２の中段に示すように、ＴＩ＝ＴＩ_thとしたときにＦＥ系及びＳＳＦＰ系それぞれの読み出しシーケンスで収集されるデータの信号強度が一致するように、少なくとも一方の読み出しシーケンスでデータ収集を行う際の受信ゲインを調整する。

【0071】

なお、図１２の中段に示す例では、ＦＥ系の読み出しシーケンスでデータ収集を行う際の受信ゲインに、ＳＳＦＰ系の読み出しシーケンスでデータ収集を行う際の受信ゲインを合わせた場合を示している。しかし、逆に、ＳＳＦＰ系の読み出しシーケンスでデータ収集を行う際の受信ゲインに、ＦＥ系の読み出しシーケンスでデータ収集を行う際の受信ゲインを合わせてもよい。また、ＦＥ系の読み出しシーケンスでデータ収集を行う際の受信ゲイン、及び、ＳＳＦＰ系の読み出しシーケンスでデータ収集を行う際の受信ゲインの両方を調整してもよい。

【0072】

また、図１２に示す例では、ＴＩ＝ＴＩ_thとしてＦＥ系及びＳＳＦＰ系の両方の読み出しシーケンスでデータ収集を行うこととしたが、必ずしも、このように異なる同じＴＩで重複してデータ収集を行わなくてもよい。その場合には、例えば、合成部３６ｅは、ＦＥ系の読み出しシーケンスで収集される最後のＴＩのデータと、ＳＳＦＰ系の読み出しシーケンスで収集される最初のＴＩのデータとで信号強度が一致するように、少なくとも一方の読み出しシーケンスでデータ収集を行う際の受信ゲインを調整する。

【0073】

また、合成部３６ｅは、上述したゲイン補正を行った後に、さらにＴ１緩和補正を行ってもよい。その場合には、例えば、図１２の下段に示すように、合成部３６ｅは、ＴＩ＝０としてＦＥ系の読み出しシーケンスで得られた血流画像における血流部分の信号強度に合うように、他のＴＩで得られた血流画像における血流部分の信号強度を補正する。例えば、合成部３６ｅは、血液のＴ１値をＴ１ｂとした場合に、各ＴＩで得られた血流画像それぞれの信号値に｛１／ｅｘｐ（−ＴＩ／Ｔ１ｂ）｝を乗じることによって、各血流画像における信号強度を補正する。例えば、このＴ１緩和補正は、ＡＳＬ−ＭＲＰで血液の流量を定量化する際に行われることが望ましい。

【0074】

また、合成部３６ｅは、上述したゲイン補正を行った後に、さらに画像補間を行ってもよい。例えば、合成部３６ｅは、待ち時間ＴＩごとに生成された血流画像について、連続するＴＩの２つの血流画像を用いて、重み付き補間により画像を生成し、生成した画像を当該２つの血流画像の間に補間する。

【0075】

このように、ＦＥ系の読み出しシーケンスを用いて得られた血流画像及びＳＳＦＰ系の読み出しシーケンスを用いて得られた血流画像の少なくとも一方についてゲイン補正を行うことによって、一方の読み出しシーケンスで得られた血流画像の信号強度を他方の読み出しシーケンスで得られた血流画像の信号強度に合わせることができ、合成画像における血管の連続性を向上させることができる。また、ゲイン補正に加えてＴ１緩和補正及び画像補間を行うことによって、合成画像における血管の連続性をさらに向上させることができる。なお、上述したゲイン補正、Ｔ１緩和補正及び画像補間は必ずしも全ての補正処理を行わなくてもよく、いずれか１つ又は２つの補正処理だけを行ってもよい。

【0076】

図１３は、第１の実施形態に係るＭＲＩ装置１００によって行われる処理の流れを示すフローチャートである。図１３に示すように、ＭＲＩ装置１００では、撮像条件設定部３６ａが、被検体の撮像に用いられる読み出しシーケンスの種類ごとに撮像条件を設定する。本実施形態では、撮像条件設定部３６ａは、少なくとも、ＦＥ系の読み出しシーケンスに関する撮像条件と、ＳＳＦＰ系の読み出しシーケンスに関する撮像条件とを設定する（ステップＳ１０１）。

【0077】

そして、操作者から撮像開始の指示を受け付けると（ステップＳ１０２，Ｙｅｓ）、分割部３６ｂが、被検体を撮像する際の時間的な範囲又は空間的な範囲を少なくとも２つの範囲に分割する。本実施形態では、分割部３６ｂは、ＳＴ−ＭＩ−ＡＳＬで被検体の頭部血管を撮像する場合に、標識化用のＲＦ波を印加してからデータ収集を行うまでの待ち時間に基づいて、撮像時の時間的な範囲を少なくとも２つの範囲に分割する（ステップＳ１０３）。

【0078】

続いて、収集部３６ｃが、分割部３６ｂによって分割された少なくとも２つの範囲のうち、第１の範囲については第１の読み出しシーケンスを用いてデータ収集を行い、第２の範囲については第１の読み出しシーケンスと種類又は撮像条件が異なる第２の読み出しシーケンスを用いてデータ収集を行う。本実施形態では、収集部３６ｃは、分割部３６ｂによって分割された２つの範囲のうち、待ち時間が短い方の範囲についてはＦＥ系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲についてはＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行う（ステップＳ１０４）。

【0079】

血流画像生成部３６ｄは、第１の読み出しシーケンスによって収集されたデータ及び第２の読み出しシーケンスによって収集されたデータそれぞれから血流画像を生成する。本実施形態では、血流画像生成部３６ｄは、ＳＴ−ＭＩ−ＡＳＬで用いられる複数の待ち時間ごとに、画像再構成部３２によって生成されたタグ画像及びコントロール画像の差分画像を血流画像として生成する（ステップＳ１０５）。

【0080】

合成部３６ｅは、第１の読み出しシーケンスによって収集されたデータから生成された画像と第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。本実施形態では、ＦＥ系の読み出しシーケンスによって収集されたデータから生成された血流画像と、ＳＳＦＰ系の読み出しシーケンスによって収集されたデータから生成された血流画像とを合成する（ステップＳ１０６）。そして、合成部３６ｅは、合成した血流画像を表示部３５に表示させる（ステップＳ１０７）。

【0081】

上述したように、第１の実施形態に係るＭＲＩ装置１００は、ＳＴ−ＭＩ−ＡＳＬで被検体の頭部血管を撮像する場合に、撮像時の時間的な範囲を少なくとも２つの範囲に分割し、待ち時間が短い方の範囲については、磁場均一性に強く、かつ、空間分解能が高いＦＥ系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲については、ＳＮ比が高いＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行うので、主幹動脈及び末梢動脈の両方を高精度に描出した画像を得ることができる。したがって、第１の実施形態に係るＭＲＩ装置１００によれば、撮像対象の部位や撮像の目的に応じて、より高画質な画像を得ることができる。

【0082】

なお、第１の実施形態では、ＳＴ−ＭＩ−ＡＳＬで被検体を撮像する場合について説明したが、撮像法がこれに限られない。例えば、本実施形態は、ＳＴ−ＳＩ（Single-Tag Single-TI）−ＡＳＬを用いる場合にも同様に実施が可能である。ＳＴ−ＳＩ−ＡＳＬは、被検体の撮像領域に流入する流体を標識化するためのＲＦ波を１回印加してから所定の待ち時間が経過した後にデータ収集を行うパターンを、待ち時間を変えて複数回実行する撮像法である。

【0083】

このようなＳＴ−ＳＩ−ＡＳＬで被検体を撮像する場合に、例えば、分割部３６ｂは、ＳＴ−ＭＩ−ＡＳＬを用いる場合と同様に、標識化用のＲＦ波を印加してからデータ収集を行うまでの待ち時間に基づいて、撮像時の時間的な範囲を少なくとも２つの範囲に分割する。そして、収集部３６ｃは、分割部３６ｂによって分割された少なくとも２つの範囲のうち、第１の範囲については第１の読み出しシーケンスを用いてデータ収集を行い、第２の範囲については第１の読み出しシーケンスと種類又は撮像条件が異なる第２の読み出しシーケンスを用いてデータ収集を行う。例えば、収集部３６ｃは、ＳＴ−ＭＩ−ＡＳＬを用いる場合と同様に、分割された２つの範囲のうち、待ち時間が短い方の範囲についてはＦＥ系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲についてはＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行う。

【0084】

さらに、１回のＲＦ波の印加でｋ空間全体のデータをサンプリングすることが難しい場合には、ｋ空間を複数の領域に分割し、分割した領域ごとに、標識化用のＲＦ波を印加してデータ収集を行ってもよい。その場合には、分割した領域全てについてデータを収集するまでは、同じ待ち時間でデータ収集を繰り返し、分割した領域全てについてデータを収集した後に、待ち時間を変更して、次のｋ空間のデータを収集する。このような場合も、同様に、分割部３６ｂが、待ち時間に基づいて、撮像時の時間的な範囲を少なくとも２つの範囲に分割する。また、収集部３６ｃが、第１の範囲については第１の読み出しシーケンスを用いてデータ収集を行い、第２の範囲については第１の読み出しシーケンスと種類又は撮像条件が異なる第２の読み出しシーケンスを用いてデータ収集を行う。

【0085】

また、第１の実施形態では、収集部３６ｃが、分割部３６ｂによって設定された待ち時間の閾値に基づいて、データ収集に用いる読み出しシーケンスを変更する場合について説明したが、読み出しシーケンスの制御はこれに限られない。例えば、収集部３６ｃは、撮像領域における所定の位置に標識化された流体が到達したことを検出した場合に、データ収集用の読み出しシーケンスを第１の読み出しシーケンスから第２の読み出しシーケンスに変更してもよい。

【0086】

例えば、収集部３６ｃは、ＳＴ−ＭＩ−ＡＳＬで被検体の頭部血管を撮像する場合に、血流画像生成部３６ｄによって血流画像が生成されるごとに、生成された血流画像における標識化された血液の位置を検出する。そして、収集部３６ｃは、検出した血液がＭＣＡの第２分枝に至る辺りに到達したことを検出した場合に、データ収集に用いる読み出しシーケンスをＦＥ系の読み出しシーケンスからＳＳＦＰ系の読み出しシーケンスに変更する。

【0087】

また、第１の実施形態では、収集部３６ｃが、ＳＴ−ＭＩ−ＡＳＬで被検体の頭部血管を撮像する場合に、データ収集用に用いる読み出しシーケンスの種類をＦＥ系からＳＳＦＰ系に変更する場合について説明したが、実施形態はこれに限られない。例えば、収集部３６ｃは、読み出しシーケンスの種類を変更するのではなく、読み出しシーケンスの撮像条件を変更してもよい。

【0088】

例えば、収集部３６ｃは、ＳＴ−ＭＩ−ＡＳＬで被検体の頭部血管を撮像する場合に、同じ種類の読み出しシーケンスを用いて、分割部３６ｂによって分割された２つの範囲のうち、待ち時間が短い方の範囲については空間分解重視の撮像条件でデータ収集を行い、待ち時間が長い方の範囲についてはＳＮ比重視の撮像条件でデータ収集を行う。この場合には、撮像条件記憶部３３ａが、同一種類の読み出しシーケンスについて、空間分解重視の撮像条件とＳＮ比重視の撮像条件とをそれぞれ記憶する。なお、空間分解重視の撮像条件とＳＮ比重視の撮像条件とは、それぞれ、収集マトリクス、ＮＡＱ、スラブ厚、１スラブ内のスライス枚数、ショット数などを適宜に調整することによって設定される。

【0089】

このように読み出しシーケンスの撮像条件を変更することによって、被検体を時系列に撮像する際の時間的な位置や撮像領域の空間的な位置に応じて、撮像法をＡＳＬ−ＭＲＡからＡＳＬ−ＭＲＰへの切り替えることも可能になる。ＡＳＬ−ＭＲＡで空間分解能重視のまま待ち時間ＴＩを延長していくと、灌流はＳＮ比の限界のためノイズレベル以下になって描出されないが、途中でＳＮ比重視モードに変更することによって、標識化された血液が待ち時間ＴＩの延長とともに毛細血管床に移行する様子をシームレスに画像化することができる。

【0090】

なお、収集部３６ｃは、データ収集に用いるパルシーケンスの種類及び撮像条件のいずれか一方を変更するだけではなく、パルシーケンスの種類及び撮像条件の両方を変更してもよい。これにより、被検体を時系列に撮像する際の時間的な位置や撮像領域の空間的な位置に応じて、より適切な方法でデータを収集することができるようになり、さらに高画質な画像を得ることができる。

【0091】

（第２の実施形態）
次に、第２の実施形態について説明する。第２の実施形態では、被検体を撮像する際の時間的な範囲又は空間的な範囲の全体について、種類が異なる２つの読み出しシーケンスそれぞれを用いてデータ収集を行う場合の例を説明する。なお、第２の実施形態に係るＭＲＩ装置の構成は、基本的には、図１及び９に示したものと同じである。

【0092】

本実施形態では、マルチスラブ法を用いたＡＳＬ−ＭＲＡで被検体の頭部血管を撮像する場合の例と、ＳＴ−ＭＩ−ＡＳＬを用いたＡＳＬ−ＭＲＡで被検体の頭部血管を撮像する場合の例とを説明する。ここで、マルチスラブ法は、被検体の撮像領域を複数の分割領域に分割して当該複数の分割領域ごとにデータ収集を行う撮像法である。

【0093】

第２の実施形態では、収集部３６ｃは、被検体を撮像する際の時間的な範囲又は空間的な範囲の全体については第１の読み出しシーケンス及び第２の読み出しシーケンスそれぞれを用いてデータ収集を行う。また、合成部３６ｅは、第１の読み出しシーケンスによって収集されたデータから生成された画像のうち第１の範囲における画像と、第２の読み出しシーケンスによって収集されたデータから生成された画像のうち第２の範囲における画像とを合成する。

【0094】

図１４は、第２の実施形態に係るＭＲＩ装置によって行われるマルチスラブ法を用いたＡＳＬ−ＭＲＡを説明するための図である。図１４では、マルチスラブ法を用いたＡＳＬ−ＭＲＡで被検体の頭部血管を撮像する場合に、主幹動脈側の分割領域１４１及び末梢動脈側の分割領域１４２の２つの領域に撮像領域を分割してデータ収集を行う場合の例を示している。

【0095】

例えば、図１４の左側に示すように、収集部３６ｃは、待ち時間ＴＩごとに、ＦＥ系の読み出しシーケンス及びＳＳＦＰ系の読み出しシーケンスそれぞれでデータ収集を行う。そして、図１４の右側に示すように、合成部３６ｅは、ＦＥ系の読み出しシーケンスを用いて得られた血流画像のうち主幹動脈側の分割領域１４１の画像と、ＳＳＦＰ系の読み出しシーケンスを用いて得られた血流画像のうち末梢動脈側の分割領域１４２の画像とを合成する。これにより、主幹動脈から末梢動脈までが良好に描出された血流画像が得られる。

【0096】

このとき、例えば、合成部３６ｅは、ＦＥ系の読み出しシーケンスで得られた血流画像については、主幹動脈側の分割領域１４１における信号値の重みを末梢動脈側の分割領域１４２における信号値の重みより大きくし、ＳＳＦＰ系の読み出しシーケンスで得られた血流画像については、末梢動脈側の分割領域１４２における信号値の重みを主幹動脈側の分割領域１４１における信号値の重みより大きくしたうえで、それぞれの血流画像を合成する。

【0097】

具体的な例として、例えば、合成部３６ｅは、ＦＥ系の読み出しシーケンスで得られた血流画像については、分割領域１４１における信号値の重み：分割領域１４２における信号値の重み＝１：０とし、ＳＳＦＰ系の読み出しシーケンスで得られた血流画像については、分割領域１４１における信号値の重み：分割領域１４２における信号値の重み＝０：１としたうえで、それぞれの血流画像を合成する。このとき、例えば、合成部３６ｅは、重み付き加算又は最大値投影（Maximum Intensity Projection：ＭＩＰ）により、血流画像を合成する。

【0098】

図１５は、第２の実施形態に係るＭＲＩ装置によって行われるＳＴ−ＭＩ−ＡＳＬを用いたＡＳＬ−ＭＲＡを説明するための図である。図１５では、ＳＴ−ＭＩ−ＡＳＬを用いたＡＳＬ−ＭＲＡで、被検体の頭部血管を撮像する場合の例を示している。図１５において、左上の図は、ＦＥ系の読み出しシーケンスでＴＩ＝６００ｍｓとしたときに得られる画像を示しており、右上の図は、ＦＥ系の読み出しシーケンスでＴＩ＝１４００ｍｓとしたときに得られる画像を示している。また、図１５において、左下の図は、ＳＳＦＰ系の読み出しシーケンスでＴＩ＝６００ｍｓとしたときに得られる画像を示しており、右下の図は、ＳＳＦＰ系の読み出しシーケンスでＴＩ＝１４００ｍｓとしたときに得られる画像を示している。なお、図１５に示す例では、標識化された血液が、ＴＩ＝６００ｍｓとしたときに主幹動脈に達し、ＴＩ＝１４００ｍｓとしたときに末梢動脈に達することとしている。

【0099】

例えば、図１５に示すように、収集部３６ｃは、複数の待ち時間ＴＩごとに、ＦＥ系の読み出しシーケンス及びＳＳＦＰ系の読み出しシーケンスそれぞれでデータ収集を行う。そして、合成部３６ｅは、ＴＩ＝６００ｍｓについてはＦＥ系の読み出しシーケンスで得られた血流画像（図１５の左上）を用い、ＴＩ＝１４００ｍｓについてはＳＳＦＰ系の読み出しシーケンスで得られた血流画像（図１５の右下）を用いて、複数のＴＩに渡るダイナミック画像を生成する。これにより、主幹動脈から末梢動脈までが良好に描出されたダイナミック画像が得られる。

【0100】

このとき、例えば、合成部３６ｅは、ＴＩ＝６００ｍｓについては、ＦＥ系の読み出しシーケンスで得られた血流画像における信号値の重みをＳＳＦＰ系の読み出しシーケンスで得られた血流画像における信号値の重みより大きくし、ＴＩ＝１４００ｍｓについては、ＳＳＦＰ系の読み出しシーケンスで得られた血流画像における信号値の重みをＦＥ系の読み出しシーケンスで得られた血流画像における信号値の重みより大きくしたうえで、ＴＩごとに、それぞれの血流画像を合成する。

【0101】

具体的な例として、例えば、合成部３６ｅは、ＴＩ＝６００ｍｓについては、ＦＥ系の読み出しシーケンスで得られた血流画像における信号値の重み：ＳＳＦＰ系の読み出しシーケンスで得られた血流画像における信号値の重み＝１：０とし、ＴＩ＝１４００ｍｓについては、ＦＥ系の読み出しシーケンスで得られた血流画像における信号値の重み：ＳＳＦＰ系の読み出しシーケンスで得られた血流画像における信号値の重み＝０：１としたうえで、ＴＩごとに、それぞれの血流画像を合成する。このとき、例えば、合成部３６ｅは、重み付き加算又は最大値投影（Maximum Intensity Projection：ＭＩＰ）により、血流画像を合成する。

【0102】

（第３の実施形態）
次に、第３の実施形態について説明する。第３の実施形態では、標識化用のＲＦ波としてパルス波を用いるＰＡＳＬ（Pulsed ASL）法で被検体の頭部血管を撮像する場合の例を説明する。なお、第３の実施形態に係るＭＲＩ装置の構成は、基本的には、図１及び９に示したものと同じである。

【0103】

例えば、ＰＡＳＬ法で被検体を撮像する場合に、標識化用のＲＦ波を印加する領域の大きさを血管の走行方向に長く設定すると、待ち時間ＴＩを長くしたときに、標識化された血液は、末梢動脈に流入した状態となるとともに、主幹動脈にも流入した状態となる。そのため、待ち時間ＴＩが長いときにＳＳＦＰ系の読み出しシーケンスを用いると、主幹動脈の部分でフローボイド（flow-void）が出現したり、血管以外の部分にゴーストなどのアーチファクトが出現したりする場合がある。このことから、血管画像を合成する場合に、待ち時間ＴＩを長くしたときの血管画像には、主幹動脈に流入した血液を残さないようにするのが望ましい。

【0104】

そこで、本実施形態に係るＭＲＩ装置は、ＰＡＳＬ法で被検体を撮像する際に、標識化用のＲＦ波をタグ領域に印加してから所定の時間が経過した後にタグ領域にサチュレーションパルスを印加することで、ＲＦ波によって標識化された血液の後端における血液信号を低減させる。ここで用いられるサチュレーションパルスは、ＴＥＣ（Tag-End-Cut）パルスと呼ばれる。このように、標識化された血液の後端部分の血液信号を低減させることで、待ち時間ＴＩを長くしたときの血管画像において、主幹動脈に流入した血液を除去することができる。

【0105】

図１６〜１８は、第３の実施形態に係るＭＲＩ装置によって行われるＰＡＳＬ法を説明するための図である。図１６は、第３の実施形態に係るＰＡＳＬ法の読み出しシーケンスを示している。また、図１７は、第３の実施形態に係るＰＡＳＬ法によって得られる血流画像を示している。また、図１８は、第３の実施形態に係るＰＡＳＬ法によって得られる合成画像を示している。

【0106】

本実施形態に係るＰＡＳＬ法の読み出しシーケンスでは、標識化用のＲＦ波を印加してから一定の時間Ｔ_tecが経過した後に、標識化用のＲＦ波を印加する時間幅と等価なＴＥＣパルスを印加する。そして、標識化用のＲＦ波を印加してからデータ収集を行うまでの待ち時間ＴＩを、Ｔ_tecの整数倍となるように段階的に延長する。例えば、待ち時間ＴＩを、短い方から、ＴＩ（１）＝Ｔ_tec、ＴＩ（２）＝２＊Ｔ_tec、・・・、ＴＩ（ｎ−１）＝（ｎ−１）＊Ｔ_tec、ＴＩ（ｎ）＝ｎ＊Ｔ_tecと変える（ｎ：自然数）。ここで、例えば、図１６に示すように、ＴＩを変える幅を少し短めにして、（Ｔ_tec−ａ）としてもよい（ａ：定数、ａ＜Ｔ_tec）。

【0107】

これにより、例えば、図１７に示すように、所定の幅を有するボーラス（bolus）状の血流が待ち時間ＴＩごとに段階的に進んでいく血流画像が得られる。そして、これらの血流画像を合成することで、例えば、図１８に示すように、主幹動脈及び末梢動脈を含んだ頭部全体の血流画像が得られる。なお、図１６に示すようにＴＩを変える幅を（Ｔ_tec−ａ）とした場合には、血流画像を合成する際に、標識化された血液の境界をオーバーラップさせることができる。これにより、ＴＩが異なる血流画像の境界で血管の端部に無信号の部分が生じるのを防ぐことができ、合成画像における血管の連続性を向上させることができる。

【0108】

そして、本実施形態では、このようなＰＡＳＬ法で被検体を撮像する場合に、収集部３６ｃが、所定の待ち時間ＴＩを境に、データ収集用の読み出しシーケンスを、ＦＥ系の読み出しシーケンスからＳＳＦＰ系の読み出しシーケンスに変更する。例えば、図１６に示すように、収集部３６ｃは、待ち時間の閾値がＴＩ３であった場合には、収集部３６ｃは、図１１に示すように、待ち時間がＴＩ３以下であるときは、ＦＥ系の読み出しシーケンスを用いてデータ収集（Ｓ１〜Ｓ３）を行い、待ち時間がＴＩ３を超えているときは、ＳＳＦＰ系の読み出しシーケンスを用いてデータ収集（Ｓ４〜Ｓ６）を行うように、シーケンサ１０を制御する。

【0109】

なお、ここでは、ＰＡＳＬ法を用いる場合について説明したが、本実施形態は、例えば、ＣＡＳＬ（Continuous ASL）法やＰＣＡＳＬ（Pulsed Continuous ASL）法を用いる場合も同様に実施が可能である。ＣＡＳＬ法は、標識化用のＲＦ波として連続波を用いる方法であり、ＰＣＡＳＬ法は、ＣＡＳＬ法の実用化を目的としたもので、標識化用のＲＦ波として、短いパルス波を複数用いる方法である。

【0110】

ＣＡＳＬ法又はＰＣＡＳＬ法では、Ｔ_tecは、標識化用のＲＦ波の印加時間Ｔ_tagに相当し、ＴＥＣパルスを印加してからデータ収集を行うまでの待ち時間は、標識化用のＲＦ波を印加してからデータ収集を行うまでの待ち時間Ｔ_pld（post labeling delay）に相当する。そのため、例えば、ＣＡＳＬ法又はＰＣＡＳＬ法を用いる場合は、Ｔ_pldをＴ_tagの整数倍となるように段階的に延長する。そして、収集部３６ｃは、ＰＡＳＬ法を用いる場合と同様に、所定の待ち時間ＴＩを境に、データ収集用の読み出しシーケンスを、ＦＥ系の読み出しシーケンスからＳＳＦＰ系の読み出しシーケンスに変更する。

【0111】

さらに、本実施形態は、上述したＰＡＳＬ法とマルチスラブ法とを組み合わせて用いる場合も同様に実施が可能である。この場合には、ＭＲＩ装置は、被検体の撮像領域を複数の分割領域に分割し、撮像領域に流入する流体を標識化するためのＲＦ波を印加してから複数の分割領域にそれぞれに対応する複数の待ち時間が経過するごとにデータ収集を行う撮像法で被検体を撮像する撮像法を行う。そして、分割部３６ｂは、標識化用のＲＦ波を印加してからデータ収集を行うまでの待ち時間に基づいて、空間的な範囲を少なくとも２つの範囲に分割する。また、収集部３６ｃは、分割された２つの範囲のうち、待ち時間が短い方の範囲については第１の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲については第２の読み出しシーケンスを用いてデータ収集を行う。そして、合成部３６ｅは、第１の読み出しシーケンスによって収集されたデータから生成された画像と第２の読み出しシーケンスによって収集されたデータから生成された画像とを合成することで撮像領域の画像を生成する。つまり、合成部３６ｅは、複数の分割領域の画像を合成することで、分割前の撮像領域の画像を生成する。

【0112】

図１９は、第３の実施形態に係るＭＲＩ装置によって行われるマルチスラブ法を用いたＰＡＳＬ法を説明するための図である。図１９では、マルチスラブ法を用いたＰＡＳＬ法で被検体の頭部血管を撮像する場合に、主幹動脈側の分割領域１９１と末梢動脈側の分割領域１９２とに撮像領域を分割してデータ収集を行う場合の例を示している。この場合に、例えば、分割部３６ｂは、被検体の頭部に設定された撮像領域を、待ち時間ＴＩに基づいて、待ち時間Ｔ１（１）及びＴＩ（２）に対応する主幹動脈側の分割領域１９１と、待ち時間Ｔ１（３）及びＴＩ（４）に対応する末梢動脈側の分割領域１９２とに分割する。

【0113】

そして、この場合に、例えば、収集部３６ｃは、待ち時間Ｔ１（１）及びＴＩ（２）についてはＦＥ系の読み出しシーケンスを用いてデータ収集を行い、待ち時間Ｔ１（３）及びＴＩ（４）についてはＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行うように、シーケンサ１０を制御する。この場合も、例えば、図１７に示したように、所定の幅を有するボーラス状の血流が待ち時間ＴＩごとに段階的に進んでいく血流画像が得られる。そして、これらの血流画像を合成することで、例えば、図１８に示したように、主幹動脈及び末梢動脈を含んだ頭部全体の血流画像が得られる。

【0114】

このように、短い待ち時間ＴＩに対応する主幹動脈側の分割領域については、磁場均一性に強く、かつ、空間分解能が高いＦＥ系の読み出しシーケンスを用いてデータ収集を行い、長い待ち時間ＴＩに対応する末梢動脈側の分割領域については、ＳＮ比が高いＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行うことによって、主幹動脈及び末梢動脈の両方を高精度に描出した画像を得ることができる。

【0115】

なお、ここでは、収集部３６ｃが、ＦＥ系の読み出しシーケンスからＳＳＦＰ系の読み出しシーケンスに読み出しシーケンスの種類を変更する場合の例を説明したが、例えば、読み出しシーケンスの条件を変更してもよい。その場合には、例えば、収集部３６ｃは、同じ種類の読み出しシーケンスを用いて、短い待ち時間ＴＩに対応する主幹動脈側の分割領域については、空間分解重視の撮像条件でデータ収集を行い、長い待ち時間ＴＩに対応する末梢動脈側の分割領域については、ＳＮ比重視の撮像条件でデータ収集を行う。

【0116】

さらに、例えば、収集部３６ｃは、単一の読み出しシーケンスで、血管信号が多い分割領域からデータを収集し、血管信号が少ない分割領域からはデータを収集しないように、シーケンサ１０を制御してもよい。なお、この場合には、各分割領域は、データを収集する領域としない領域とを細かく区分けできるように、それぞれ血管の走行方向における幅を狭く設定するのが望ましい。

【0117】

（第４の実施形態）
次に、第４の実施形態について説明する。第４の実施形態では、被検体の撮像領域に流入する流体を標識化するためのＲＦ波を印加してから所定の待ち時間が経過した後にデータ収集を行うパターンを待ち時間を変えて複数回実行する撮像法で被検体を撮像する場合の例を説明する。なお、第４の実施形態に係るＭＲＩ装置の構成は、基本的には、図１及び９に示したものと同じである。

【0118】

例えば、このような撮像法の一例として、Ｈａｄａｍａｒｄ Ｅｎｃｏｄｉｎｇ（ＨＥ）法を用いたＭｕｌｔｉ−Ｔａｇ（ＭＴ）−ＰＡＳＬであるＨＥ−ＭＴ−ＰＡＳＬがある。ＨＥ法については、例えば、ＣＡＳＬ法と組み合わせた方法が、Wells JA et al. In vivo hadamard encoded continuous arterial spin labeling (H-CASL). MRM 63:1111-1118 (2010)により提案されている。ＨＥ−ＭＴ−ＰＡＳＬでは、タグ画像用のＲＦ波とコントロール画像用のＲＦ波とを所定の時間間隔ごとに効率よく配置してデータを収集し、収集した画像を加減算することで、待ち時間ＴＩが異なる複数のＡＳＬ画像が得られる。この方法によれば、ＳＴ−ＳＩ−ＡＳＬやＳＴ−ＭＩ−ＡＳＬと比べて、ＳＮ比をより向上させることができる。

【0119】

図２０は、第４の実施形態に係るＭＲＩ装置によって行われるＨＥ−ＭＴ−ＰＡＳＬを説明するための図である。ＨＥ−ＭＴ−ＰＡＳＬでは、例えば、異なる待ち時間ＴＩの組み合わせ数Ｎが、３、７、１５・・・など、すなわち、Ｎ＝２ⁿ−１（ｎ：自然数）に設定される。図２０では、異なる待ち時間ＴＩの組み合わせ数をＮ＝７（ｎ＝３）とした場合の例を示している。

【0120】

なお、図２０に示す複数のＲＦ波（ｔａｇ／ｃｏｎｔｒｏｌ ＲＦ）のうち、左下がりの斜線が付されたものはタグ画像用のＲＦ波であり、右下がりの斜線が付されたものはコントロール画像用のＲＦ波である。図２０に示すように、例えば、Ｎ＝７とした場合には、２種類の画像用のＲＦ波の順序が異なる８種類（（Ｎ＋１）種類）の読み出しシーケンスの組み合わせが用いられる。また、各読み出しシーケンスにおける待ち時間ＴＩ（ｉ）は、ｋ＝１〜８の各パターンで同じであれば、任意の時間間隔でよい。

【0121】

ここで、ｋ＝１〜８の各パターンで収集されるデータの複素信号をＳｋとし、各ＴＩ（ｉ）における４倍の差分信号Ｓ｛ＴＩ（ｉ）｝、ｉ＝１〜７は、以下の式で求められる。

【0122】

Ｓ｛ＴＩ（ｉ）｝＝４［Ｓ｛ＴＩ（ｉ）｝−Ｓｃｏｎｔ｛ＴＩ（ｉ）｝］

【0123】

すなわち、ｉ＝１〜７とした場合には、Ｓ｛ＴＩ（１）｝〜Ｓ｛ＴＩ（７）｝は、以下の式により求められる。

【0124】

Ｓ｛ＴＩ（１）｝＝Ｓ１−Ｓ２＋Ｓ３−Ｓ４＋Ｓ５−Ｓ６＋Ｓ７−Ｓ８
Ｓ｛ＴＩ（２）｝＝Ｓ１＋Ｓ２−Ｓ３−Ｓ４＋Ｓ５＋Ｓ６−Ｓ７−Ｓ８
Ｓ｛ＴＩ（３）｝＝Ｓ１−Ｓ２−Ｓ３＋Ｓ４＋Ｓ５−Ｓ６−Ｓ７＋Ｓ８
Ｓ｛ＴＩ（４）｝＝Ｓ１＋Ｓ２＋Ｓ３＋Ｓ４−Ｓ５−Ｓ６−Ｓ７−Ｓ８
Ｓ｛ＴＩ（５）｝＝Ｓ１−Ｓ２＋Ｓ３−Ｓ４＋Ｓ５−Ｓ６＋Ｓ７−Ｓ８
Ｓ｛ＴＩ（６）｝＝Ｓ１＋Ｓ２−Ｓ３−Ｓ４−Ｓ５−Ｓ６＋Ｓ７＋Ｓ８
Ｓ｛ＴＩ（７）｝＝Ｓ１−Ｓ２−Ｓ３＋Ｓ４−Ｓ５＋Ｓ６＋Ｓ７−Ｓ８

【0125】

そして、この場合のＳＮ比は、１回あたりの差分のＮＡＱ＝４回加算に相当するｓｑｒｔ（４）＝２倍となる。

【0126】

また、１回の繰り返し時間をＴ_repeatとすると、撮像時間は、Ｎ組のＴＩに対して、（Ｎ＋１）＊Ｔ_repeatとなる。一方、標識用のＲＦ波を１回印加するごとにデータ収集を行う場合には、撮像時間は、Ｎ組のＴＩに対して、２＊（Ｎ＋１）／２＊Ｎ＊Ｔ_repeat＝Ｎ（Ｎ＋１）＊Ｔ_repeatとなる。したがって、ＨＥ−ＭＴ−ＳＩ法では、標識用のＲＦ波を１回印加するごとにデータ収集を行う場合に比べて、撮像時間比はＮ＋１／Ｎ（Ｎ＋１）＝１／Ｎとなる。つまり、Ｎ＝７とした場合には、標識用のＲＦ波を１回印加するごとにデータ収集を行う場合と同等のＳＮ比のＴＩの異なる差分画像が、１／７の時間で得られることになる。

【0127】

なお、ＡＳＬ−ＭＲＰにおいて一般的なＳＴ−ＰＡＳＬでは、十分なＳＮ比を確保するためには数回程度の加算平均が行われる。これに対し、上述したＨＥ−ＭＴ−ＰＡＳＬによれば、同一のデータが１回の収集のみで得られ、かつ、数段階のＴＩの加算平均と等価な画像が得られるので、収集時間を有効に短縮することができる。なお、コントロール画像用のＲＦ波は、ＭＲＰ用でＭＴＣが問題となる場合のみ印加されればよい。したがって、血管のＭＲＡ用などでＭＴＣが問題とならない場合には、コントロール画像用のＲＦ波は省略されてもよい。

【0128】

また、上述したＨＥ−ＭＴ−ＰＡＳＬでは、血液の流速に応じたパルス間隔の制御など、タグ条件に制約がある。さらに、ＨＥ法では、血液の流速が時間的に変動する場合、例えば拍動流下で用いられる場合には、同一時刻でのＲＦ波のタイミングで血液の流速が同じでないと、画像を加減算する際、特に減算の際に、不要なタイミングで標識化された流体信号が意図した通りに消えない場合もある。これは、ＰＡＳＬ法と組み合わせた場合でも、ＣＡＳＬ法と組み合わせた場合でも同様である。そこで、例えば、収集部３６ｃは、同一時刻のＲＦ波が同一の心位相で行われるように、ＥＣＧセンサ２１及びＥＣＧユニット２２によって検出される心拍のゲート信号におけるＲ波をトリガ信号として、最初のＲＦ波を印加するようにしてもよい。

【0129】

そして、本実施形態では、収集部３６ｃは、上述したＨＥ−ＭＴ−ＰＡＳＬを用いて被検体を撮像する場合に、ｋ＝１〜８の８種類のパターンのうち、いずれか４種類のパターンについてはＦＥ系の読み出しシーケンスでデータ収集を行い、残りの４種類のパターンについてはＳＳＦＰ系の読み出しシーケンスでデータ収集を行う。この場合には、平均効果によって、ＳＳＦＰ系の読み出しシーケンスで得られる画像とＦＥ系の読み出しシーケンスで得られる画像とが混じった画像が得られる。

【0130】

例えば、図２０に示すように、収集部３６ｃは、８種類のパターンを順に実行する際に、ＦＥ系の読み出しシーケンスとＳＳＦＰ系の読み出しシーケンスとを交互に用いるように、シーケンサ１０を制御する。すなわち、収集部３６ｃは、ｋ＝１、３、５、７のパターンについてはＳＳＦＰ系の読み出しシーケンスでデータ収集を行い、ｋ＝２、４、６、８のパターンについてはＦＥ系の読み出しシーケンスでデータ収集を行うように制御する。

【0131】

または、例えば、収集部３６ｃは、ＴＩが短い側にタグ画像用のＲＦ波を多く配置したパターンについてはＦＥ系の読み出しシーケンスを用い、ＴＩが長い側にタグ画像用のＲＦ波を多く配置したパターンについてはＦＥ系の読み出しシーケンスを用いるようにしてもよい。

【0132】

さらに、被検体の撮像領域に流入する流体を標識化するためのＲＦ波を印加してから所定の待ち時間が経過した後にデータ収集を行うパターンを待ち時間を変えて複数回実行する撮像法の他の例として、ＳＴ−ＭＩ−ＡＳＬを複数回実行する撮像法がある。この撮像法では、標識化用のＲＦ波を印加してから最初のデータ収集を行うまでの待ち時間を変えながら、ＳＴ−ＭＩ−ＡＳＬのパターンを複数回実行し、収集された複数のデータを、ＲＦ波が印加されてからデータ収集を行うまでの待ち時間が同じものごとに加算する例について説明する。

【0133】

図２１及び２２は、第４の実施形態に係るＭＲＩ装置によって行われるＳＴ−ＭＲ−ＡＳＬを説明するための図である。図２１に示すように、例えば、収集部３６ｃは、ＳＴ−ＭＩ−ＡＳＬによるデータ収集をＭ回実行する。この場合に、例えば、収集部３６ｃは、標識用のＲＦ波を印加してから最初のデータ収集までの待ち時間を所定の時間幅Ｔ_intの整数倍となるように変えながら、ＳＴ−ＭＩ−ＡＳＬによるデータ収集をＮ回実行する。また、例えば、収集部３６ｃは、ＳＴ−ＭＩ−ＡＳＬによるデータ収集を実行するごとに、標識化用のＲＦ波を印加してから所定の時間幅Ｔ_intの整数倍の時刻ごとに、Ｎ個のＭＲデータを収集する。なお、本実施形態では、標識用のＲＦ波が印加されてから最初のデータ収集を行うまでの待ち時間と、２回目のデータ収集において標識用のＲＦ波が印加されてからデータ収集を行うまでの経過時間とを、それぞれＴＩと呼ぶ。

【0134】

ここで、例えば、収集部３６ｃは、各ＴＩをＴＩｍｎと表した場合に、以下に示すＴＩ１１、ＴＩ２１、ＴＩ３１、ＴＩ４１のように、標識用のＲＦ波が印加されてから最初のデータ収集を行うまでの待ち時間を変化させながら、ＳＴ−ＭＩ−ＰＡＳＬによるデータ収集を実行する。そして、収集部３６ｃは、１回目については、待ち時間ＴＩ１１が経過した後にデータ収集を行い、その後、ＴＩ１２、ＴＩ１３、ＴＩ１４が経過するごとにデータ収集を実行する。また、収集部３６ｃは、２回目については、待ち時間ＴＩ２１が経過した後にデータ収集を行い、その後、ＴＩ２２、ＴＩ２３、ＴＩ２４が経過するごとにデータ収集を実行する。また、収集部３６ｃは、３回目については、待ち時間ＴＩ３１が経過した後にデータ収集を行い、その後、ＴＩ３２、ＴＩ３３、ＴＩ３４が経過するごとにデータ収集を実行する。また、収集部３６ｃは、４回目については、待ち時間ＴＩ４１が経過した後にデータ収集を行い、その後、ＴＩ４２、ＴＩ４３、ＴＩ４４が経過するごとにデータ収集を実行する。

【0135】

ＴＩ１４＝７Ｔ_int
ＴＩ１３＝ＴＩ２４＝６Ｔ_int
ＴＩ１２＝ＴＩ２３＝ＴＩ３４＝５Ｔ_int
ＴＩ１１＝ＴＩ２２＝ＴＩ３３＝ＴＩ４４＝４Ｔ_int
ＴＩ２１＝ＴＩ３２＝ＴＩ４３＝３Ｔ_int
ＴＩ３１＝ＴＩ４２＝２Ｔ_int
ＴＩ４１＝１Ｔ_int

【0136】

そして、収集部３６ｃは、収集された複数のＭＲデータを補正したうえで、標識用のＲＦ波が印加されてからデータ収集を行うまでの経過時間が同じものごとに加算する。例えば、収集部３６ｃは、加算されるＭＲデータをＳ（ＴＩ）とし、ＴＩに依存したＴ１減衰を補正するための補正係数をａ（ＴＩ）とし、データを収集する順番ｎに依存したＴ１減衰を補正するための補正係数をｂｎとした場合に、以下の式によりＳ（ＴＩ）を求める。

【0137】

S（7T_int）＝a（7T_int）＊b4＊S14
S（6T_int）＝a（6T_int）＊｛b3＊S13＋b4＊S24｝／2
S（5T_int）＝a（5T_int）＊｛b2＊S12＋b3＊S23＋b4＊S34｝／3
S（4T_int）＝a（4T_int）＊｛b1＊S11＋b2＊S22＋b3＊S33＋b4＊S44｝／4
S（3T_int）＝a（5T_int）＊｛b1＊S21＋b2＊S32＋b3＊S43｝／3
S（2T_int）＝a（2T_int）＊｛b1＊S31＋b2＊S42｝／2
S（1T_int）＝a（1T_int）＊b1＊S41

【0138】

この方法によれば、図２１に示すように、Ｍ＋Ｎ−１個の異なるＴＩのＭＲデータが得られ、かつ、中間のＴＩほどＳＮ比が向上する。なお、ＴＩが中間より短い場合にはもともとＳＮ比が大きいので問題ないが、ＴＩが中間より長い場合にはＳＮ比は低下し、最も長いＴＩのＳＮ比は収集されるデータが１個なので、１組のＭＴ−ＳＩ−ＰＡＳＬと同じになる。そのため、Ｔ_intごとに標識用のＲＦ波とデータ収集とを配置するのが望ましい。さらに、Ｎ回のデータ収集の組と最初のデータ収集までの待ち時間ＴＩをＴ_intずつシフトさせたｍ組の独立したｔａｇのＭ組として、ｍ＝ｎとするのがより望ましい。なお、ｍ、ｎは異なってもよい。

【0139】

また、ＴＩが最も短くなる最後（ｍ＝Ｍ）のＲＦ波と１回目（ｎ＝１）のデータ収集までの時間はＴ_intである必要はない。また、ＴＩは、Ｔ_intの整数倍であればよく、ここでいう整数は必ずしも連続していなくてもよい。抜けたＴＩについては、合成計算時にゼロとしてないものとすればよい。

【0140】

なお、ＡＳＬでは、ＴＩが長いほどＳＮ比は低下するので加算回数を増やしたいが、ＳＴ−ＭＩ−ＰＡＳＬでは時間的に後のデータほどＳＮ比は低下するため、あまり多くのｄａｔａは収集できない。したがって、ｎは数個にしておき、独立したＲＦ波で収集したｍ組のデータから合成して中間以降の長いＴＩでＳＮ比の小さなデータを捨てたとしても、処理効率を向上できる。

【0141】

そして、本実施形態では、収集部３６ｃは、上述した撮像法を用いて被検体を撮像する場合に、標識用のＲＦ波が印加されてから最初のデータ収集を行うまでの待ち時間が短いパターンについてはＦＥ系の読み出しシーケンスでデータ収集を行い、標識用のＲＦ波が印加されてから最初のデータ収集を行うまでの待ち時間が長いパターンについてはＳＳＦＰ系の読み出しシーケンスでデータ収集を行う。

【0142】

例えば、図２１に示すように、収集部３６ｃは、標識用のＲＦ波が印加されてから最初のデータ収集を行うまでの待ち時間をＴ_intずつ短くしながらＳＴ−ＭＴ−ＡＳＬのパターンを４回実行する場合に、１回目及び２回目のパターンについてはＳＳＦＰ系の読み出しシーケンスでデータ収集を行い、３回目及び４回目のパターンについてはＳＳＦＰ系の読み出しシーケンスでデータ収集を行うように、シーケンサ１０を制御する。

【0143】

または、収集部３６ｃは、各パターンにおけるデータ収集ごとに、標識用のＲＦ波が印加されてからデータ収集を行うまでの待ち時間が短いものについてはＦＥ系の読み出しシーケンスを用い、標識用のＲＦ波が印加されてからデータ収集を行うまでの待ち時間が長いものについてはＳＳＦＰ系の読み出しシーケンスを用いるようにしてもよい。

【0144】

例えば、図２２に示すように、収集部３６ｃは、４回のパターンについて、標識用のＲＦ波が印加されてからデータ収集を行うまでの待ち時間が４Ｔ_int以下であるデータ収集Ｓ１１、Ｓ２１、Ｓ２２、Ｓ３１〜Ｓ３３、Ｓ４１〜Ｓ４４についてはＦＥ系の読み出しシーケンスを用い、標識用のＲＦ波が印加されてからデータ収集を行うまでの待ち時間が４Ｔ_intより長いデータ収集Ｓ１２〜Ｓ１４、Ｓ２３〜Ｓ２４、Ｓ３４についてはＳＳＦＰ系の読み出しシーケンスを用いるように、シーケンサ１０を制御する。

【0145】

上述した第１〜４の実施形態に係るＭＲＩ装置によれば、単一の読み出しシーケンスを用いる場合と比べて、全ての待ち時間ＴＩにおいて、磁化率効果の影響が少ない、かつ、ＳＮ比に優れた、臓器全体のＭＲＡ画像又はＭＲＰ画像を得ることができる。また、第１〜４の実施形態に係るＭＲＩ装置によれば、例えば脳の動脈血管を撮像する場合に、主幹動脈だけでなく、末梢動脈における血流の描出能を向上させることができる。また、第１〜４の実施形態に係るＭＲＩ装置によれば、標識化用のＲＦ波を印加してから標識化された血液が到達するまでの遅延時間が長い組織の部分について、ＳＮ比が高い灌流画像が得られる。

【0146】

（第５の実施形態）
次に、第５の実施形態について説明する。第５の実施形態では、被検体の撮像領域を複数の分割領域に分割して当該複数の分割領域ごとにデータ収集を行う撮像法で被検体の頭部血管を撮像する場合の例を説明する。なお、第５の実施形態に係るＭＲＩ装置の構成は、基本的には、図１及び９に示したものと同じである。

【0147】

例えば、このような撮像法の一例として、マルチスラブ法を適用した３ＤＴＯＦ−ＭＲＡであるマルチスラブ３ＤＴＯＦ−ＭＲＡがある。マルチスラブ３ＤＴＯＦ−ＭＲＡでは、標識化用のＲＦ波を印加する代わりに、縦磁化が十分に回復した血液のインフロー（in-flow）効果と、静止部における縦磁化の抑制効果とによるコントラスト差を利用するため、撮像領域（imaging slab）のみが設定される。

【0148】

このようなマルチスラブ３ＤＴＯＦ−ＭＲＡで被検体を撮像する場合に、分割部３６ｂは、分割領域に基づいて、被検体を撮像する際の空間的な範囲を少なくとも２つの範囲に分割する。また、収集部３６ｃは、分割部３６ｂによって分割された２つの範囲のうち、主幹動脈を含む範囲については第１の読み出しシーケンスを用いてデータ収集を行い、末梢動脈を含む範囲については第２の読み出しシーケンスを用いてデータ収集を行う。

【0149】

図２３は、第５の実施形態に係るＭＲＩ装置によって行われるマルチスラブ３ＤＴＯＦ−ＭＲＡを説明するための図である。図２３に示すように、例えば、マルチスラブ３ＤＴＯＦ−ＭＲＡでは、被検体の撮像領域を４つの分割領域に分割して、それぞれの分割領域についてデータ収集が行われる。

【0150】

この場合に、例えば、分割部３６ｂは、標識化された血液が主幹動脈を流れる範囲と、末梢動脈を流れる範囲とに空間的な範囲を分割する。また、例えば、収集部３６ｃは、ウィリス動脈輪部の近傍にある主幹動脈部については、磁場不均一性に強いＦＥ系の読み出しシーケンスを用いてデータ収集を行い、頭頂などの末梢動脈部については、ＳＮ比に優れたＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行うように、シーケンサ１０を制御する。

【0151】

なお、上述したマルチスラブ３ＤＴＯＦ−ＭＲＡでは、撮像領域に照射するＲＦ波のフリップ角に傾斜をつけることで末梢動脈の描出能の向上を図るＴＯＮＥ法が併用されてもよい。また、２Ｄマルチスライス法でも適用が可能である。その場合、読み出しシーケンスの撮像条件を変えてもよい。ＴＯＦ法は基本的にはＧＲＥであるが、その中の撮像条件うち、ＴＯＮＥ法での傾斜は対象スラブ内の血流速に応じて変化させてもよい。すなわち、流速が遅いほど傾斜を急にする。またＴＥとＧＭＮ（Gradient Moment Nulling）の組み合わせを、短いＴＥ＋０次ＧＭＮと、長いＴＥ＋１次ＧＭＮとで使い分ける。脳の主幹動脈部など磁化率が大きく流速が速ければ前者、抹消で逆なら後者にすることで、位相分散による信号低下を抑制する。さらに、ＢＷ（BandWidth）を変えたものを組み合わせても良い。抹消は狭いＢＷでＳＮ比を向上させ、主幹部は広めのＢＷで位相分散による信号低下を抑制する。

【0152】

（第６の実施形態）
次に、第６の実施形態について説明する。第６の実施形態では、マルチスラブ法による撮像をＴＯＦ−ＭＲＡ及びＡＳＬ−ＭＲＡ／ＭＲＰと組み合わせたマルチスラブＡＳＬで被検体の頭部血管を撮像する場合の例を説明する。なお、第２の実施形態に係るＭＲＩ装置の構成は、基本的には、図１及び９に示したものと同じである。

【0153】

マルチスラブＡＳＬでは、被検体の撮像領域を複数の分割領域に分割し、各分割領域について、当該分割領域から一定の距離だけ流体の上流側に離れた標識化領域に当該分割領域に流入する流体を標識化するためのＲＦ波を印加してから一定の待ち時間が経過した後にデータ収集を行う撮像法で被検体を撮像する。

【0154】

図２４は、第６の実施形態に係るＭＲＩ装置によって行われるマルチスラブＡＳＬを説明するための図である。図２４に示すように、例えば、マルチスラブＡＳＬでは、撮像領域（imaging slab）を４つの分割領域１２１〜１２４に分割する。そして、４つの分割領域それぞれのデータを収集する際に、各分割領域から一定の距離だけ血流の上流側に離れたタグ領域（tag slab）に標識化用のＲＦ波を印加する。すなわち、マルチスラブＡＳＬでは、データ収集を行う分割領域を移動しながら、分割領域の移動にともなって、タグ領域を分割領域から一定の距離を保った状態で追随させて移動する。なお、図２４では図示を省略しているが、コントロール領域についても、それぞれの分割領域に対して相対的な位置を同じにして追随させる。

【0155】

このようなマルチスラブＡＳＬによれば、複数の分割領域についてタグ領域の位置を撮像領域全体の上流側に固定する場合と比べて、下流側の分割領域になるほど、標識化された血液が到達するまでの時間が短くなり、ＳＮ比を向上させることができる。そのため、標識化用のＲＦ波を印加してからデータ収集を行うまでの待ち時間ＴＩを延長しなくても、タグ領域から遠い分割領域内にある血管（ＭＲＡ）や灌流（ＭＲＰ）も良好に描出することができる。

【0156】

このようなマルチスラブＡＳＬで被検体の頭部血管を撮像する場合に、例えば、分割部３６ｂは、分割領域に基づいて、撮像領域を少なくとも２つの範囲に分割する。また、収集部３６ｃは、分割部３６ｂによって分割された範囲のうち、主幹動脈を含む範囲については第１の読み出しシーケンスを用いてデータ収集を行い、末梢動脈を含む範囲については第２の読み出しシーケンスを用いてデータ収集を行う。

【0157】

例えば、図２４に示すように、収集部３６ｃは、４つの分割領域のうち、主幹動脈を含む上流側の分割領域１２１及び１２２については、ＦＥ系の読み出しシーケンスを用いてデータ収集を行う。また、収集部３６ｃは、末梢動脈を含む下流側の分割領域１２３及び１２４については、ＳＳＦＰ系の読み出しシーケンスを用いてデータ収集を行う。これにより、被検体の頭部血管が全体にわたって良好に描出することができる。

【0158】

上述した第５及び６の実施形態に係るＭＲＩ装置によれば、単一の読み出しシーケンスを用いる場合と比べて、撮像時間を延長することなく、磁化率効果の影響の少ない、かつ、ＳＮ比が優れた画像を得ることができる。また、第５及び６の実施形態に係るＭＲＩ装置によれば、例えば脳の動脈血管を撮像する場合に、主幹動脈だけでなく、末梢動脈における血流の描出能を向上させることができる。

【0159】

（第７の実施形態）
次に、第７の実施形態について説明する。第７の実施形態では、拡散イメージングで被検体を撮像する場合の例を説明する。なお、第７の実施形態に係るＭＲＩ装置の構成は、基本的には、図１及び９に示したものと同じである。

【0160】

近年、拡散イメージングは、脳神経領域だけでなく、腫瘍検出用に全身臓器でも使用されてきている。また、拡散イメージングでは、一般的に、ＳＥ系のシングルショットＳＥ−ＥＰＩ（Ｓ−ＥＰＩ）の読み出しシーケンスが用いられる。

【0161】

ここで、脳神経領域では、小脳、前髄など脳底部や頸部、胸部の脊髄などの部分、腹部領域では、肺や肺に近い肝臓部や直腸や大腸に近い部分は、磁場不均一性が大きくなることが知られている。しかし、Ｓ−ＥＰＩの読み出しシーケンスは磁場不均一性に弱く、歪みが大きくなったり、別のボクセルに描出されるべき信号が１つのボクセル内で重なって描出されてしまったりする場合もある。

【0162】

また、拡散イメージングで用いられる他の読み出しシーケンスには、ＦＳＥ系（例えば、ＨＡＳＴＥ：Half-Fourier Single shot Turbo spin Echo）や、ナビゲータエコー（navigator echo）を同時に収集して動きによるショットごとの位相差を補正するマルチショットＥＰＩ（Ｍ−ＥＰＩ）がある。さらに、ｋ空間の軌跡（trajectory）の違いによって、ｎｏｎ−Ｃａｒｔｅｓｉａｎ系のＰＲＯＰＥＬＬＥＲ（ＪＥＴ）、ｒａｄｉａｌ系などがある。しかし、これらの読み出しシーケンスは、Ｓ−ＥＰＩの読み出しシーケンスと比べて、単位時間当たりのＳＮ比が低いなどの理由で、１種類の読み出しシーケンスや撮像条件で全身などの広い領域を撮像するのは難しかった。

【0163】

このような拡散イメージングで被検体を撮像する場合に、分割部３６ｂは、静磁場強度（Ｂ０）の分布に基づいて、被検体を撮像する際の空間的な範囲を少なくとも２つの範囲に分割する。例えば、分割部３６ｂは、ｄｅｌｔａＢ０の標準偏差に関する所定の閾値に基づいて、被検体を撮像する際の撮像領域を２つに分割する。

【0164】

また、収集部３６ｃは、分割部３６ｂによって分割された２つの範囲のうち、静磁場不均一性が大きい方の範囲については第１の読み出しシーケンスを用いてデータ収集を行い、静磁場不均一性が小さい方の範囲については第２の読み出しシーケンスを用いてデータ収集を行う。例えば、収集部３６ｃは、２つの範囲のうち、ｄｅｌｔａＢ０の標準偏差が大きい方を静磁場不均一性が大きい範囲とし、ｄｅｌｔａＢ０の標準偏差が小さい方を静磁場不均一性が小さい範囲とする。なお、撮像対象が頭部であれば、静磁場不均一性が大きい方の範囲は頸部や脳底部などに相当し、静磁場不均一性が小さい方の範囲は頭頂部などに相当する。また、撮像対象が腹部であれば、静磁場不均一性が大きい方の範囲は肺や肝臓の肺に近い部分などに相当し、静磁場不均一性が小さい方の範囲は腔腹膜臓器などに相当する。

【0165】

例えば、収集部３６ｃは、静磁場不均一性が大きい方の範囲については、歪みに強いＦＳＥ系の読み出しシーケンスを用いてデータ収集を行う。また、収集部３６ｃは、静磁場不均一性が小さい方の範囲については、ＦＳＥ系と比べて歪みには強くないが、ＳＮ比が高いＳ−ＥＰＩ系の読み出しシーケンスを用いてデータ収集を行う。

【0166】

また、例えば、収集部３６ｃは、静磁場不均一性が大きい方の範囲については、ＦＳＥ系の読み出しシーケンスでＰＲＯＰＥＬＬＥＲ（ＪＥＴ）を用いてデータ収集を行い、静磁場不均一性が小さい方の範囲については、Ｓ−ＥＰＩの読み出しシーケンスでデータ収集を行ってもよい。

【0167】

なお、例えば、収集部３６ｃは、静磁場不均一性が大きい方の範囲については、同じ読み出しシーケンス又は異なる読み出しシーケンスを用いて、歪みに強い撮像条件でデータ収集を行い、静磁場不均一性が小さい方の範囲については、ＳＮ比が高い撮像条件でデータ収集を行ってもよい。このとき、例えば、収集部３６ｃは、歪みに強い撮像条件として、ＳＥＮＳＥ−Ｆａｃｔｏｒを高くしたものを用い、ＳＮ比が高い撮像条件として、ＳＥＮＳＥ−Ｆａｃｔｏｒを低くしたものを用いる。

【0168】

また、例えば、収集部３６ｃは、同じＥＰＩの読み出しシーケンスを用いて、静磁場不均一性が大きい方の範囲については、ショット数を多くした撮像条件でデータ収集を行い、静磁場不均一性が小さい方の範囲については、ショット数を少なくした撮像条件でデータ収集を行ってもよい。例えば、収集部３６ｃは、静磁場不均一性が大きい方の範囲については、Ｍ−ＥＰＩの読み出しシーケンスを用いてデータ収集を行い、静磁場不均一性が小さい方の範囲については、Ｓ−ＥＰＩの読み出しシーケンスでデータ収集を行う。

【0169】

また、拡散イメージングでは、ＭＰＧ強度やｂ−ｖａｌｕｅが大きくなると、渦電流による歪みも大きくなることが知られている。そこで、例えば、収集部３６ｃは、静磁場不均一性が大きい方の範囲については、ＭＰＧ強度やｂ−ｖａｌｕｅを小さくした撮像条件でデータ収集を行い、静磁場不均一性が小さい方の範囲については、ＭＰＧ強度やｂ−ｖａｌｕｅを大きくした撮像条件でデータ収集を行ってもよい。

【0170】

（第８の実施形態）
次に、第８の実施形態について説明する。第８の実施形態では、定量値画像（パラメータ画像）を撮像する場合の例を説明する。なお、第８の実施形態に係るＭＲＩ装置の構成は、基本的には、図１及び９に示したものと同じである。

【0171】

近年、検査間での比較が容易なことから、Ｔ１やＴ２、プロトン密度、ＡＤＣ、血流量、ケミカルシフト、磁化率、温度などのＭＲＩに関する各種パラメータを画像化することが行われている。

【0172】

このように、所定の定量値を空間マッピングするための撮像法で被検体を撮像する場合に、分割部３６ｂは、対象部位の特性に応じて、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割する。そして、収集部３６ｃは、分割部３６ｂによって分割された複数の範囲ごとに、各範囲に含まれる対象部位の特性に応じて、データ収集用の読み出しシーケンスの種類又は撮像条件を変更する。

【0173】

以上、第１〜８の実施形態に係るＭＲＩ装置について説明したが、実施形態はこれに限られない。例えば、上述した実施形態では、分割部３６ｂが、撮像対象の部位に関する解剖学的な情報に基づいて、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割する場合について説明したが、分割の方法はこれに限られない。

【0174】

例えば、分割部３６ｂは、シミング用に測定された静磁場強度の分布に基づいて、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割してもよい。例えば、分割部３６ｂは、ｄｅｌｔａＢ０に関する標準偏差や平均値などの統計値に基づいて、時間的な範囲又は空間的な範囲を少なくとも２つの範囲に分割する。この場合には、分割部３６ｂは、ｄｅｌｔａＢ０に関する統計値について所定の閾値を設定し、その閾値を境に、時間的な範囲又は空間的な範囲を分割する。

【0175】

また、例えば、分割部３６ｂは、撮像領域の位置決めに用いられる画像に基づいて、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割してもよい。例えば、分割部３６ｂは、位置決め画像として撮像されたアキシャル像、コロナル像又はサジタル像に対して操作者が指定した位置を境に、時間的な範囲又は空間的な範囲を分割する。この場合、例えば、分割部３６ｂは、３ＤＴＯＦの画像を取得済みであれば、３ＤＴＯＦの画像を位置決め用の画像として用い、未取得であれば、２ＤＰＣ（Phase Contrast）法で取得した画像を位置決め用の画像として用いてもよい。

【0176】

また、上述した実施形態では、ＦＥ系の読み出しシーケンスとＳＳＦＰ系の読み出しシーケンスとを用いた場合の例を中心に説明したが、読み出しシーケンスの種類はこれに限られない。すなわち、収集部３６ｃは、分割部３６ｂによって分割された少なくとも２つの範囲について、撮像対象の部位や撮像の目的に応じて、読み出しシーケンスの種類を適宜に使い分ける。例えば、収集部３６ｃは、２つの範囲のうち、一方の範囲でＳＮ比が求められる場合は、その範囲については、他方の範囲で用いる読み出しシーケンスよりも単位収集時間あたりのＳＮ比が高い読み出しシーケンスを用いる。また、例えば、収集部３６ｃは、２つの範囲のうち、一方の範囲で空間分解能が求められる場合は、その範囲については、他方の範囲で用いる読み出しシーケンスよりも空間分解能が高い読み出しシーケンスを用いる。また、例えば、収集部３６ｃは、２つの範囲のうち、一方の範囲が静磁場不均一性の影響を受けやすい場合は、その範囲については、他方の範囲で用いる読み出しシーケンスよりも静磁場不均一性に強い読み出しシーケンスを用いる。

【0177】

また、収集部３６ｃは、分割部３６ｂによって分割された少なくとも２つの範囲について、撮像対象の部位や撮像の目的に応じて、読み出しシーケンスの撮像条件を適宜に使い分けてもよい。例えば、収集部３６ｃは、２つの範囲のうち、一方の範囲でＳＮ比が求められる場合は、その範囲については、他方の範囲で用いる撮像条件よりも単位収集時間あたりのＳＮ比が高い撮像条件を用いる。また、例えば、収集部３６ｃは、２つの範囲のうち、一方の範囲で空間分解能が求められる場合は、その範囲については、他方の範囲で用いる撮像条件よりも空間分解能が高い撮像条件を用いる。また、例えば、収集部３６ｃは、２つの範囲のうち、一方の範囲が静磁場不均一性の影響を受けやすい場合は、その範囲については、他方の範囲で用いる撮像条件よりも静磁場不均一性に強い撮像条件を用いる。

【0178】

また、上述した実施形態では、被検体の頭部血管を撮像する場合の例を中心に説明したが、撮像対象となる部位は頭部に限られない。例えば、撮像対象の部位は、腹部などの他の部位であってもよい。また、上述した実施形態では、血流画像を生成する場合の例を中心に説明したが、撮像対象となる流体は血液に限らない。例えば、撮像対象の流体は、ＣＳＦ（cerebrospinal fluid）などの他の流体であってもよい。

【0179】

前述したように、読み出しシーケンスの撮像条件は、所定の撮像領域を撮像する際に用いられる撮像パラメータのパラメータ値である。つまり、撮像条件が異なるということは、あるスライス位置のデータを収集する際の撮像パラメータのパラメータ値が異なることを意味する。例えば、撮像領域を時間的な範囲に分割する場合には、同じスライス位置について、分割された時間的な範囲ごとに、読み出しシーケンスの撮像パラメータのパラメータ値を変えてデータ収集が行われる。また、例えば、撮像領域を空間的な範囲に分割する場合には、異なるスライス位置について、分割された空間的な範囲ごとに、読み出しシーケンスの撮像パラメータのパラメータ値を変えてデータ収集が行われる。

【0180】

上述した第１〜８の実施形態に係るＭＲＩ装置によれば、時間的な範囲又は空間的な範囲の異なる全ての部分について、単一の読み出しシーケンス種又は単一の撮像条件で撮像する場合と比べて、ほぼ同一の撮像時間で、よりアーチファクトの少ない、空間分解能やＳＮ比に優れたＭＲＩ画像を得ることができる。

【0181】

また、上述した第１〜８の実施形態では、撮像領域を時間的又空間的な範囲に分割し、分割された範囲ごとに、読み出しシーケンスの種類及び撮像条件の少なくとも一方を変える場合の例を説明したが、実施形態はこれに限られない。例えば、分割された範囲ごとに、読み出しシーケンスの前に印加されるプリパルスの条件を変えてもよい。

【0182】

例えば、分割された範囲ごとに、ＡＳＬ法で用いられる標識化用のＲＦパルスのタグ条件を変えてもよい。例えば、第６の実施形態で説明したように、収集部が、マルチスラブＡＳＬで被検体を撮像する場合に、時間的又空間的に分割された範囲ごとに、データ収集を行う分割領域の移動に追随させて、標識化用のＲＦパルスを印加する領域すなわちタグ領域の位置を移動させる。このとき、例えば、収集部は、タグ領域と分割領域との相対位置が変わらないように、タグ領域の絶対位置を変える。

【0183】

また、例えば、アディアバティックパルスを用いた撮像が行われる場合に、分割された範囲ごとに、アディアバティックパルスである反転パルスのＲＦパワー（フリップ角）を変えてもよい。アディアバティックパルスは、理論的には、縦磁化が１８０度倒れる以上のＲＦパワー（フリップ角）に対しても縦磁化は１８０度以上回転することはない。しかし、ＲＦパワーが大きくなれば、送信ＲＦコイルに送信されるＲＦ波を増幅するＲＦアンプの負荷が大きくなるので、短時間に複数のアディアバティックパルスを印加する場合は、デューティサイクルを考慮した最適な大きさにＲＦパワーを配分する必要がある。例えば、送信ＲＦコイルの感度が不均一であることを考慮して、収集部が、分割された範囲のうち、送信ＲＦパワーが低めになるＲＦコイルの中心から遠い範囲については、反転パルスのフリップ角を３００度とし、送信ＲＦパワーが十分なＲＦコイルの中心に近い範囲については、反転パルスのフリップ角を１８０度とするように制御する。

【0184】

上述した少なくとも１つの実施形態に係るＭＲＩ装置によれば、撮像対象の部位や撮像の目的に応じて、より高画質な画像を得ることができる。

【0185】

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

【符号の説明】

【0186】

１００ ＭＲＩ（磁気共鳴イメージング）装置
３０ 計算機システム
３６ 制御部
３６ｂ 分割部
３６ｃ 収集部
３６ｅ 合成部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】

【図21】

【図22】

【図23】

【図24】