JP5965880B2 - 磁気共鳴装置およびプログラム - Google Patents

Description

本発明は、血液のデータを収集する磁気共鳴装置、およびこの磁気共鳴装置に適用することができるプログラムに関する。

血流を撮影する方法として、流入（inflow）効果を利用した方法が知られている（特許文献１参照）。

特開２００８−０８６７４８号公報

また、流入効果を利用して血管画像を取得する手法として、ＩＦＩＲ（Inhance Inflow IR）シーケンスを用いた方法が知られている。ＩＦＩＲシーケンスでは、血管を撮影したい領域の磁化を反転させるための反転パルスを印加する。そして、反転パルスを印加してから反転時間が経過するのを待ち、データ収集を実行する。この方法では、反転パルスを印加してからデータを収集するまでに間に、縦磁化の十分に大きい血液が血管を撮影したい領域に流入するので、血管が描出された画像を取得することができる。

この方法では、血液の流入効果を利用して血管を撮影するので、反転パルスを印加してからデータを収集するまでの間に、縦磁化の十分に大きい血液を、血管を撮影したい領域の全体に流入させる必要がある。しかし、反転時間の設定値によっては、縦磁化の十分に大きい血液を、血管を撮影したい領域の全体に十分に流入させることができず、所望の血管画像を得ることができないことがある。また、反転時間を長くしすぎた場合、データ収集前に背景組織が回復してしまい、画像のコントラストが低下するという問題がある。

したがって、縦磁化の十分に大きい血液を、血管を撮影したい領域の全体に十分に流入させることができ、更に、画像のコントラストを十分に高くすることができる撮影方法が要求されている。

本発明の第１の観点は、血液を含む領域内のスピンをフリップさせるための第１のＲＦパルスと前記領域から血液のデータを収集するデータ収集シーケンスとを有するイメージングシーケンスを繰り返し実行する磁気共鳴装置であって、
前記血液と背景組織とのコントラストと、前記第１のＲＦパルスから前記データ収集シーケンスまでの第１の時間と、イメージングシーケンスが終了してから次のイメージングシーケンスが開始されるまでの第２の時間との対応関係を記憶する記憶部と、
前記血液の流速に基づいて、前記イメージングシーケンスを繰り返し実行するときの前記第１の時間を決定する第１の決定手段と、
前記第１の決定手段により決定された前記第１の時間と、前記対応関係とに基づいて、前記イメージングシーケンスを繰り返し実行するときの前記第２の時間を決定する第２の決定手段とを有する磁気共鳴装置である。

本発明の第２の観点は、血液を含む領域内のスピンをフリップさせるための第１のＲＦパルスと前記領域から血液のデータを収集するデータ収集シーケンスとを有するイメージングシーケンスを繰り返し実行する磁気共鳴装置であって、前記血液と背景組織とのコントラストと、前記第１のＲＦパルスから前記データ収集シーケンスまでの第１の時間と、イメージングシーケンスが終了してから次のイメージングシーケンスが開始されるまでの第２の時間との対応関係を記憶する磁気共鳴装置に適用されるプログラムであって、
前記血液の流速に基づいて、前記イメージングシーケンスを繰り返し実行するときの前記第１の時間を決定する第１の決定処理と、
前記第１の決定手段により決定された前記第１の時間と、前記対応関係とに基づいて、前記イメージングシーケンスを繰り返し実行するときの前記第２の時間を決定する第２の決定処理とを計算機に実行させるためのプログラムである。

血液の流速に基づいて第１の時間を決定するので、撮影したい領域の全体に渡って血管を描出することができる。また、上記対応関係に基づいて第２の時間を決定するので、コントラストの高い血管画像を取得することができる。

本発明の一形態の磁気共鳴装置の概略図である。 被検体１２をスキャンするときの説明図である。 スキャンを実行するときのフローを示す図である。 反転時間ＴＩの決定方法の説明図である。 待ち時間Ｔｗの許容範囲を求めるときに使用されるコントラストマップＭ１を示す図である。 待ち時間Ｔｗの許容範囲を求めるときに使用されるコントラストマップＭ２を示す図である。 ステップＳＴ４で実行されるフローを示す図である。 ステップＳＴ４０の具体的なフローを示す図である。 抽出されたマップデータＤ１を示す図である。 マップデータＤ１の中に含まれるｒ＝０．５の範囲ｖ１を示す図である。 ｒ＝０．５の範囲ｖ１に対応する待ち時間Ｔｗの範囲Ｈ１を示す図である。 抽出されたマップデータＤ２を示す図である。 マップデータＤ２の中に含まれるｓ＝０．５〜０．５２の範囲ｖ２を示す図である。 ｓ＝０．５〜０．５２の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を示す図である。 コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２とを比較して示す図である。 ＴＩ＝１５００ｍｓおよびＴｗ＝２０００ｍｓの場合のシーケンスＰＳを示す図である。 抽出されたマップデータＤ１を示す図である。 マップデータＤ１の中に含まれるｒ＝０．５の範囲ｖ１を示す図である。 ｒ＝０．５の範囲ｖ１に対応する待ち時間Ｔｗの範囲Ｈ１を示す図である。 抽出されたマップデータＤ２を示す図である。 マップデータＤ２の中に含まれるｓ＝０．５〜０．５４の範囲ｖ２を示す図である。 ｓ＝０．５〜０．５４の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を示す図である。 コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２とを比較して示す図である。 ＴＩ＝１４００ｍｓおよびＴｗ＝２１００ｍｓの場合のシーケンスＰＳを示す図である。 抽出されたマップデータＤ１を示す図である。 ｒ＝０．４５に対応する待ち時間Ｔｗの範囲Ｈ１を示す図である。 抽出されたマップデータＤ２を示す図である。 マップデータＤ２の中に含まれるｓ＝０．５〜０．５４の範囲ｖ２を示す図である。 ｓ＝０．５〜０．５４の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を示す図である。 コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２とを比較して示す図である。 ＴＩ＝１３００ｍｓおよびＴｗ＝２２００ｍｓの場合のシーケンスＰＳを示す図である。 α°パルスを用いた場合のイメージングシーケンスＰＳの一例を示す図である。 時間Ｔαを、α°パルスとデータ収集シーケンスＤＡＱの終了時点ｔｅとの間の時間に設定した例を示す図である。 時間Ｔαを、α°パルスと、データ収集シーケンスＤＡＱが実行されている途中の時点ｔｍとの間の時間に設定した例を示す図である。

以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。

図１は、本発明の一形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。

マグネット２は、被検体１２が収容されるボア２１を有している。また、マグネット２には、超伝導コイル、ＲＦコイル、勾配コイルなどが内蔵されている。

テーブル３は、被検体１２を支持するクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１２はボア２１に搬送される。

受信コイル４は、被検体１２に取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。

ＭＲ装置１００は、更に、生体信号処理部５、送信器６、勾配磁場電源７、受信器８、記憶部８０、制御部９、操作部１０、および表示部１１などを有している。

心拍信号処理部５は、被検体１２に取り付けられたセンサ５ａから信号を受け取り、被検体の心拍数やＲＲ間隔を求める。

送信器６はＲＦコイルに電流を供給し、勾配磁場電源７は勾配コイルに電流を供給する。受信器８は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。

記憶部８０は、コントラストマップＭ１およびＭ２（図５および図６参照）を記憶する。コントラストマップＭ１およびＭ２については後述する。

制御部９は、表示部１１に必要な情報を伝送したり、受信器８から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部９は、流速算出手段９１〜待ち時間決定手段９３などを有している。

流速算出手段９１は、スキャンにより取得された流速情報に基づいて動脈血の流速ｖを求める。

反転時間決定手段９２は、動脈血の流速ｖに基づいて、イメージングシーケンスを実行するときの反転時間ＴＩを決定する。

待ち時間決定手段９３は、反転時間決定手段９２により決定された反転時間ＴＩと、コントラストマップＭ１およびＭ２（図５および図６参照）とに基づいて、イメージングシーケンスＰＳを実行するときの待ち時間Ｔｗを決定する。

尚、制御部９は、流速算出手段９１〜待ち時間決定手段９３を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。

操作部１０は、オペレータにより操作され、種々の情報を制御部９に入力する。表示部１１は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
上記のように構成されたＭＲＩ装置１００を用いて、被検体１２を撮影する。

図２は、被検体１２をスキャンするときの説明図である。
図２の上側には、被検体の心拍信号ＣＳと、被検体をスキャンするときに使用されるイメージングシーケンスＰＳが示されている。図２の下側には、被検体の撮影領域Ｒと、ｋ空間（ｋｙ−ｋｚ面）が示されている。
本形態では、頭部および頸部を流れる動脈血を描出するためのイメージングシーケンスＰＳが繰返し実行される。

各イメージングシーケンスＰＳは、選択反転パルスＳＩＲ（Selective Inversion Recovery）、脂肪抑制パルスＦ、データ収集シーケンスＤＡＱ、非選択反転パルスＮＩＲを有している。

選択反転パルスＳＩＲは、被検体１２の頭部および頸部を含む領域Ｒ_ｉｎｖの組織（動脈血、静脈血、脂肪、筋肉など）の縦磁化を反転させるためのＲＦパルスである。選択反転パルスＳＩＲは、心拍信号ＣＳのＲ波から遅延時間ＴＤが経過した時点で印加される。

選択反転パルスＳＩＲから反転時間ＴＩが経過した時点で、撮影領域Ｒのデータを収集するためのデータ収集シーケンスＤＡＱが実行される。データ収集シーケンスＤＡＱは、例えば３Ｄ ＦＳＥ（Fast Spin Echo）や、ＦＩＥＳＴＡ（Fast Imaging Employing Steady state Acquisition）である。選択反転パルスＳＩＲによって領域Ｒ_ｉｎｖの組織の縦磁化は反転するので、反転時間ＴＩの間に、領域Ｒ_ｉｎｖ内の各組織の縦磁化はヌルポイントに近づく。一方、心臓は領域Ｒ_ｉｎｖの外側に位置しているので、選択反転パルスＳＩＲを印加しても、心臓の中の動脈血は、縦磁化Ｍ＝１のままである。したがって、反転時間ＴＩの間に、心臓から、縦磁化Ｍ＝１の動脈血が頸部および頭部に流入する。このように、反転時間ＴＩの間に縦磁化Ｍ＝１の動脈血を頸部および頭部に流入させてから、データ収集シーケンスＤＡＱを実行するので、動脈血が強調して描出されるとともに背景組織（静脈血など）が抑制されたＭＲ画像を得ることができる。

また、データ収集シーケンスＤＡＱの直前には、脂肪抑制パルスＦが印加されている。したがって、撮影領域Ｒの脂肪信号を効果的に抑制することができる。尚、脂肪抑制パルスＦは、例えばＳＰＥＣＩＲ（Spectrally Selected IR）や、ＳＴＩＲ（Short-TI IR）である。

更に、データ収集シーケンスＤＡＱの直後には、非選択反転パルスＮＩＲが印加される。非選択反転パルスＮＩＲは、被検体内の各組織の磁化を反転させるために印加されるパルスである。

非選択反転パルスＮＩＲを印加した後、待ち時間Ｔｗが経過した時点で、次のイメージングシーケンスＰＳが実行される。

次のイメージングシーケンスＰＳでも、選択反転パルスＳＩＲおよび脂肪抑制パルスＦを印加した後、データ収集シーケンスＤＡＱを実行し、非選択反転パルスＮＩＲを印加する。以下同様に、イメージングシーケンスＰＳを繰り返し実行する。本形態では、１回のイメージングシーケンスＰＳで、ｋｙ−ｋｚ面内の１本のｋｚビューのデータを収集する。したがって、ｍ回のイメージングシーケンスＰＳを実行することにより、ｋｙ−ｋｚ面内の全てのｋｚビュー１〜ｍを収集することができる。

次に、繰り返し時間ＴＲについて説明する。
繰り返し時間ＴＲは、以下の式で表すことができる。
ＴＲ＝ＴＩ＋Ｔａ＋Ｔｗ ・・・（１）
ここで、ＴＩ：反転時間
Ｔａ：データ収集期間ＤＡＱが開始してから非選択反転パルスＮＩＲが印加されるまでの時間
Ｔｗ：待ち時間

また、ＴＲは、以下の条件を満たすように設定される。
ＴＲ＝ＲＲ×ｎ ・・・（２）
ここで、ＲＲ：ＲＲ間隔
ｎ：整数

本形態ではｎの値は、反転時間ＴＩおよび待ち時間Ｔｗに依存して決まる値である。ｎの決定方法については後述する。図２では、ｎの一例として、ｎ＝４の場合が示されている。

上記のイメージングシーケンスＰＳを実行してｋ空間のデータを収集した後、フーリエ変換を実行することにより、動脈血の画像を得ることができる。ただし、反転時間ＴＩが短すぎると、血液の流入効果が十分に得ることができず、動脈血を高信号で描出することができないという問題がある。一方、反転時間ＴＩが長すぎると、背景組織が回復し、画像のコントラストが低下する問題がある。更に、血流速や撮影部位によって、反転時間ＴＩの最適値は異なる。したがって、高品質な血流画像を得るためには、各スキャンに適した反転時間ＴＩを設定する必要がある。また、背景組織と動脈血とのコントラストは、待ち時間Ｔｗにも依存するので、良好なコントラストが得られるように待ち時間Ｔｗを設定する必要もある。そこで、本形態では、背景組織と動脈血とのコントラストを大きくすることができるように反転時間ＴＩおよび待ち時間Ｔｗの値を決定し、イメージングシーケンスＰＳを実行している。以下に、反転時間ＴＩおよび待ち時間Ｔｗの値を決定し、イメージングシーケンスＰＳを実行するときのフローについて説明する。

図３は、反転時間ＴＩおよび待ち時間Ｔｗの値を決定し、イメージングシーケンスＰＳを実行するときのフローを示す図である。
本形態のフローは、大きく２つのステップＳＴ１０およびＳＴ２０に分けられる。ステップＳＴ１０は、反転時間ＴＩおよび待ち時間Ｔｗの値を決定するためのステップである。ステップＳＴ２０では、ステップＳＴ１０で決定された反転時間ＴＩおよび待ち時間Ｔｗに基づいて、イメージングシーケンスを実行するステップである。

ステップＳＴ１０はステップＳＴ１〜ＳＴ４を有している。したがって、ステップＳＴ１０の説明に当たっては、ステップＳＴ１〜ＳＴ４について順に説明する。

ステップＳＴ１では、被検体の動脈血の流速情報を取得するためのスキャンを実行する。動脈血の流速情報は、反転時間ＴＩを決定するときに用いられる情報である。動脈血の流速情報を用いて反転時間ＴＩを決定する具体的な手順についてはステップＳＴ２およびＳＴ３で詳しく説明する。動脈血の流速情報を取得するためのスキャンとしては、位相コントラストＭＲＡなどを用いることができる。このスキャンを実行した後、ステップＳＴ２に進む。

ステップＳＴ２では、流速算出手段９１（図１参照）が、ステップＳＴ１のスキャンにより取得された流速情報に基づいて、動脈血の流速ｖを求める。動脈血の流速ｖを求めた後、ステップＳＴ３に進む。

ステップＳＴ３では、反転時間決定手段９２（図１参照）が、動脈血の流速ｖに基づいて、反転時間ＴＩを決定する。

図４は、反転時間ＴＩの決定方法の説明図である。
図４では、心臓から血管Ｖを通って頭部に流れる動脈血Ｂが矢印で表されている。

先ず、撮影領域Ｒの中で、特に高品質でイメージングしたい血液の領域ＲＱを考える。本形態では、領域ＲＱは、頸動脈を含む領域であるとする。

次に、動脈血Ｂが、位置Ｐｅから位置Ｐｃまで移動するのに必要な時間Ｔｍについて考える。位置Ｐｅは領域Ｒ_ｉｎｖの心臓側のエッジｅの位置であり、位置Ｐｃは領域ＲＱの右端の位置である。

位置Ｐｅと位置Ｐｃとの間では、血管Ｖの形状が、ＳＩ方向に延びる直線形状に近似できるとする。この場合、動脈血Ｂが位置Ｐｅから位置Ｐｃまで移動するのに必要な時間Ｔｍは、以下の式で表すことができる。
Ｔｍ＝Ｌ／ｖ ・・・（３）
ここで、Ｌ：位置Ｐｅから位置Ｐｃまでの距離
ｖ：動脈血の流速

領域ＲＱ内の血液を高品質でイメージングするためには、動脈血Ｂは、時間Ｔｍの間に、位置Ｐｅから位置Ｐｃに到達していることが望ましい。そこで、本形態では、時間Ｔｍ（＝Ｌ／ｖ）を、反転時間ＴＩとする。したがって、反転時間ＴＩは、式（４）で表することができる。
ＴＩ＝Ｌ／ｖ ・・・（４）

式（４）の動脈血の流速ｖはステップＳＴ２で求められているので、距離Ｌの値が決まれば、反転時間ＴＩを求めることができる。距離Ｌは、例えば、撮影領域ＲのＳＩ方向の長さに応じて決定してもよいし、固定値としてもよい。本形態では、距離Ｌは、撮影領域ＲのＳＩ方向の長さによって決まる値とする。したがって、流速ｖおよび距離Ｌは既知であるので、式（４）から反転時間ＴＩを決定することができる。尚、上記のように、ＴＩは流速ｖおよび距離Ｌによって決まるので（式（４）参照）、ＴＩはｖおよびＬの値に応じて種々の値を取り得る。本形態では、説明の便宜上、ＴＩの値として、以下の３つの値を考える。
ＴＩ＝１５００ｍｓ
ＴＩ＝１４００ｍｓ
ＴＩ＝１３００ｍｓ

そして、ＴＩ＝１５００ｍｓ、１４００ｍｓ、１３００ｍｓの３つの場合に分けて、フローの説明を続ける。

（１）ステップＳＴ３においてＴＩ＝１５００ｍｓと決定された場合
ステップＳＴ３においてＴＩ＝１５００ｍｓと決定されると、ステップＳＴ４に進む。

ステップＳＴ４では、待ち時間決定手段９３（図１参照）が、ステップＳＴ３で決定した反転時間ＴＩと、コントラストマップとに基づいて、待ち時間Ｔｗを決定する。以下に、待ち時間Ｔｗを決定するときに使用されるコントラストマップについて説明する（図５および図６参照）。

図５および図６は、待ち時間Ｔｗを決定するときに使用されるコントラストマップＭ１およびＭ２を示す図である。
コントラストマップＭ１（図５参照）は、動脈血とＣＳＦ（脳脊髄液）とのコントラストｒを表すマップであり、被検体をスキャンする前に予め用意されている。コントラストマップＭ１の横軸は反転時間ＴＩ、縦軸は待ち時間Ｔｗを表しており、反転時間ＴＩの値と待ち時間Ｔｗの値との組合せに対してコントラストｒの値が対応付けられている。コントラストｒの値は、例えば、ブロッホ（Bloch）方程式を用いて動脈血の信号値とＣＳＦの信号値とを計算し、これらの信号値の比を取ることにより求めることができる。

コントラストｒは０．２５≦ｒ≦０．５５の範囲内の値を有しており、コントラストマップＭ１ではコントラストｒの値の違いがグレーの濃淡で表現されている。ｒ＝０．２５は黒に対応し、ｒ＝０．５５は白に対応している。コントラストｒが０．２５から大きくなるにつれて、コントラストマップＭ１の色は黒から次第に白に近づく。

一方、コントラストマップＭ２（図６参照）は、動脈血と静脈血とのコントラストｓを表すマップであり、被検体をスキャンする前に予め用意されている。コントラストマップＭ２の横軸は反転時間ＴＩ、縦軸は待ち時間Ｔｗを表しており、反転時間ＴＩの値と待ち時間Ｔｗの値との組合せに対してコントラストｓの値が対応付けられている。コントラストｓの値は、例えば、ブロッホ（Bloch）方程式を用いて動脈血の信号値と静脈血の信号値とを計算し、これらの信号値の比を取ることにより求めることができる。

コントラストｓは０．２５≦ｓ≦０．５５の範囲内の値を有しており、コントラストマップＭ２ではコントラストｓの値の違いがグレーの濃淡で表現されている。ｓ＝０．２５は黒に対応し、ｓ＝０．５５は白に対応している。コントラストｓが０．２５から大きくなるにつれて、コントラストマップＭ２の色は黒から次第に白に近づく。

ステップＳＴ４では、図５および図６に示すコントラストマップＭ１およびＭ２を用いて、待ち時間Ｔｗが決定される。以下に、ステップＳＴ４の動作について説明する。

図７は、ステップＳＴ４で実行されるフローを示す図である。
ステップＳＴ４は、２つのステップＳＴ４０およびＳＴ４１に分けられる。ステップＳＴ４０では、待ち時間Ｔｗとして使用可能な時間範囲を決定するための処理が実行される。ステップＳＴ４０の具体的なフローは図８に示されている。図８については後述する。

ステップＳＴ４１では、ステップＳＴ４０で決定された時間範囲の中から、シーケンスを実行するときの待ち時間Ｔｗを決定する処理が実行される。

以下に、ステップＳＴ４の各ステップＳＴ４０およびＳＴ４１について説明する。
ステップＳＴ４では、先ず、ステップＳＴ４０が実行される（図８参照）。

図８はステップＳＴ４０の詳細なフローを示す図である。
ステップ４０ａでは、待ち時間決定手段９３が、コントラストマップＭ１（図５参照）の中から、ステップＳＴ３で決定された反転時間ＴＩ＝１５００（ｍｓ）におけるマップデータを抽出する。図９に抽出されたマップデータＤ１を示す。マップデータＤ１を抽出した後、ステップ４０ｂに進む。

ステップ４０ｂでは、待ち時間決定手段９３が、マップデータＤ１の中に、以下の条件を満たすコントラストｒが含まれているか否かを判断する。
ｒ≧ｒ_０ ・・・（５）
ここで、ｒ_０：動脈血とＣＳＦとのコントラストｒの値として許容される下限値

本形態では、ｒ_０＝０．５に設定されているとする。この場合、式（５）は以下の式で表される。
ｒ≧０．５ ・・・（６）

したがって、待ち時間決定手段９３は、待ち時間ＴＩ＝１５００（ｍｓ）におけるマップデータＤ１の中に、式（６）を満たすコントラストｒが含まれているか否かを判断する。マップデータＤ１の中に式（６）を満たすコントラストｒが含まれている場合、ステップ４０ｄに進む。一方、マップデータＤ１の中に式（６）を満たすコントラストｒが含まれていない場合、ステップ４０ｃに進む。図９を参照すると、マップデータＤ１はコントラストｒ＝０．５の領域を横切っていることが分かる。したがって、マップデータＤ１の中に式（６）を満たすコントラストｒが含まれている。図１０に、マップデータＤ１の中に含まれるｒ＝０．５の範囲ｖ１を示す。マップデータＤ１の中に式（６）を満たすコントラストｒが含まれているので、ステップ４０ｄに進む。

ステップ４０ｄでは、待ち時間決定手段９３が、ｒ＝０．５の範囲ｖ１に対応する待ち時間Ｔｗの範囲を特定する。図１１に、ｒ＝０．５の範囲ｖ１に対応する待ち時間Ｔｗの範囲Ｈ１を示す。ここでは、待ち時間Ｔｗの範囲Ｈ１は、以下の範囲とする。
１１００（ｍｓ）≦Ｔｗ≦２２００（ｍｓ） ・・・（７）

コントラストマップＭ１における待ち時間Ｔｗの範囲Ｈ１を求めた後、ステップ４０ｅに進む。

ステップ４０ｅでは、待ち時間決定手段９３が、コントラストマップＭ２（図６参照）の中から、ステップＳＴ３で決定された反転時間ＴＩ＝１５００（ｍｓ）におけるマップデータを抽出する。図１２に抽出されたマップデータＤ２を示す。マップデータＤ２を抽出した後、ステップ４０ｆに進む。

ステップ４０ｆでは、待ち時間決定手段９３が、マップデータＤ２の中に、以下の条件を満たすコントラストｓが含まれているか否かを判断する。
ｓ≧ｓ_０ ・・・（８）
ここで、ｓ_０：動脈血と静脈血とのコントラストｓの値として許容される下限値

本形態では、ｓ_０＝０．５に設定されているとする。この場合、式（８）は以下の式で表される。
ｓ≧０．５ ・・・（９）

したがって、待ち時間決定手段９３は、待ち時間ＴＩ＝１５００（ｍｓ）におけるマップデータＤ２の中に、式（９）を満たすコントラストｓが含まれているか否かを判断する。マップデータＤ２の中に式（９）を満たすコントラストｓが含まれている場合、ステップ４０ｈに進む。一方、マップデータＤ２の中に式（９）を満たすコントラストｓが含まれていない場合、ステップ４０ｇに進む。図１２を参照すると、マップデータＤ２はコントラスト０．５〜０．５２の領域を横切っていることが分かる。したがって、マップデータＤ２の中に式（９）を満たすコントラストｓが含まれている。図１３に、マップデータＤ２の中に含まれるｓ＝０．５〜０．５２の範囲ｖ２を示す。マップデータＤ２の中に式（９）を満たすコントラストｓが含まれているので、ステップ４０ｈに進む。

ステップ４０ｈでは、待ち時間決定手段９３が、ｓ＝０．５〜０．５２の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を特定する。図１４に、ｓ＝０．５〜０．５２の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を示す。ここでは、待ち時間Ｔｗの範囲Ｈ２は、以下の範囲とする。
１８００（ｍｓ）≦Ｔｗ≦４０００（ｍｓ） ・・・（１０）

コントラストマップＭ２における待ち時間Ｔｗの範囲Ｈ２を求めた後、ステップ４０ｉに進む。

ステップ４０ｉでは、待ち時間決定手段９３が、コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１（図１１参照）と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２（図１４参照）に、重なる部分があるか否かを判断する。図１５に、コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２とを比較して示す。

待ち時間Ｔｗの範囲Ｈ１とＨ２との間に待ち時間Ｔｗの重なる範囲Ｈがある場合はステップ４０ｋに進み、一方、待ち時間Ｔｗの重なる範囲Ｈが無い場合はステップ４０ｊに進む。図１５では、１８００ｍｓ≦Ｔｗ≦２２００ｍｓの範囲で待ち時間Ｔｗが重なっている。したがって、ステップ４０ｋに進む。

ステップ４０ｋでは、待ち時間決定手段９３は、待ち時間Ｔｗの重なっている範囲Ｈを、待ち時間Ｔｗとして使用可能な時間範囲と決定する。ここでは、１８００ｍｓ≦Ｔｗ≦２２００ｍｓの範囲で待ち時間Ｔｗが重なっている。したがって、待ち時間Ｔｗとして使用可能な時間範囲Ｈは、以下の範囲となる。
１８００ｍｓ≦Ｔｗ≦２２００ｍｓ ・・・（１１）

待ち時間Ｔｗとして使用可能な時間範囲が決定されたら、ステップＳＴ４０のフローを終了する。ステップＳＴ４０が終了したら、ステップＳＴ４１（図７参照）に進む。

ステップＳＴ４１では、待ち時間決定手段９３が、待ち時間Ｔｗとして使用可能な時間範囲Ｈの中から、イメージングシーケンスＰＳを実行するときの待ち時間Ｔｗを決定する。以下に待ち時間Ｔｗの決定方法について説明する。

式（１）より、待ち時間Ｔｗは以下の式で表すことができる。
Ｔｗ＝ＴＲ−（ＴＩ＋Ｔａ） ・・・（１２）

式（２）を式（１２）に代入すると、以下の式が得られる。
Ｔｗ＝ＲＲ×ｎ×−（ＴＩ＋Ｔａ） ・・・（１３）

式（１３）において、ＲＲは心拍信号から算出することができる値である。ここでは、ＲＲ＝１０００ｍｓとする。また、時間Ｔａは、データ収集シーケンスＤＡＱの開始時点から非選択反転パルスＮＩＲの印加時点までの時間である。時間ＴａはイメージングシーケンスＰＳによって決まる固定値であるので、時間Ｔａは既知の値である。ここでは、Ｔａ＝５００ｍｓとする。したがって、式（１３）は、以下の式で表される。
Ｔｗ＝１０００ｎ−（ＴＩ＋５００） ・・・（１４）

また、反転時間ＴＩは、ステップＳＴ３でＴＩ＝１５００ｍｓと算出されている。したがって、ＴＩ＝１５００ｍｓを式（１４）に代入すると、以下の式が得られる。
Ｔｗ＝１０００ｎ−（ＴＩ＋５００）
＝１０００ｎ−（１５００＋５００）
＝１０００ｎ−２０００ ・・・（１５）

また、ステップＳＴ４０において、待ち時間Ｔｗとして使用可能な時間範囲Ｈは１８００ｍｓ≦Ｔｗ≦２２００ｍｓとして求められている（式（１１）参照）。したがって、式（１５）で表される待ち時間Ｔｗは、以下の式を満たす必要がある。
１８００ｍｓ≦Ｔｗ＝１０００ｎ−２０００≦２２００ｍｓ ・・・（１６）

待ち時間決定手段９３は、式（１６）を満たすｎを求める。式（１６）を満たすｎは、ｎ＝４である。ｎ＝４を求めた後、ｎ＝４を式（１５）に代入する。ｎ＝４を式（１５）に代入することによって、以下のように、待ち時間Ｔｗを求めることができる。
Ｔｗ＝１０００ｎ−２０００
＝１０００×４−２０００
＝２０００ｍｓ

したがって、待ち時間Ｔｗ＝２０００ｍｓと決定することができる。
以上説明したように、ステップＳＴ３およびＳＴ４の処理を実行することによって、ＴＩ＝１５００ｍｓおよびＴｗ＝２０００ｍｓが求められる。図１６に、ＴＩ＝１５００ｍｓおよびＴｗ＝２０００ｍｓの場合のイメージングシーケンスＰＳを示す。尚、ｎ＝４であるので、ｎ＝４を式（２）に代入するとＴＲ＝４ＲＲとなる。
待ち時間Ｔｗを決定した後、ステップＳＴ２０（図３参照）に進む。

ステップＳＴ２０では、図１６に示すイメージングシーケンスＰＳに従ってｋ空間のデータを収集し、フローを終了する。

本形態では、動脈血Ｂが位置Ｐｅから位置Ｐｃまで移動するのに必要な時間Ｔｍを求め（図４および式（３）参照）、時間Ｔｍを反転時間ＴＩとして設定している。したがって、縦磁化が十分に大きい動脈血Ｂが領域ＲＱ（図４参照）の全体に流入したときにデータ収集シーケンスＤＡＱが実行されるので、領域ＲＱから十分に大きな動脈血信号を得ることができる。

また、本形態では、コントラストマップＭ１およびＭ２に基づいて、動脈血とＣＳＦとのコントラストｒがｒ≧０．５になるときの待ち時間Ｔｗの範囲Ｈ１（図１１参照）と、動脈血と静脈血とのコントラストｓがｓ≧０．５になるときの待ち時間Ｔｗの範囲Ｈ２（図１４参照）とを求める。そして、両方の範囲Ｈ１およびＨ２に重なる待ち時間Ｔｗの範囲Ｈを特定し、この範囲Ｈの中から、イメージングシーケンスＰＳを実行するときの待ち時間Ｔｗを決定している。したがって、動脈血とＣＳＦとのコントラストｒが十分に大きく、且つ動脈血と静脈血とのコントラストｓも十分に大きくすることができるので、良好なコントラストの血管画像を得ることができる。

尚、式（１３）では、ＲＲ間隔を用いて待ち時間Ｔｗを規定しているが、ＲＲ間隔の代わりに心拍数を用いて待ち時間Ｔｗを規定してもよい。心拍数をＢＰＭとすると、ＲＲ間隔と心拍数との間にはＲＲ＝（６０／ＢＰＭ）×１０^３（ｍｓ）の関係があるので、ＲＲ間隔の代わりに心拍数ＢＰＭを算出しても、待ち時間Ｔｗを決定することができる。

上記の説明では、ステップＳＴ３においてＴＩ＝１５００ｍｓと決定された例が述べられている。次に、ＴＩ＝１４００ｍｓの場合について説明する。

（２）ステップＳＴ３においてＴＩ＝１４００ｍｓと決定された場合
ステップＳＴ３においてＴＩ＝１４００ｍｓと決定されると、ＴＩ＝１５００ｍｓの場合と同様に、ステップ４０ａ（図８参照）に進む。

ステップ４０ａでは、待ち時間決定手段９３が、コントラストマップＭ１（図５参照）の中から、ステップＳＴ３で決定された反転時間ＴＩ＝１４００（ｍｓ）におけるマップデータを抽出する。図１７に抽出されたマップデータＤ１を示す。マップデータＤ１を抽出した後、ステップ４０ｂに進む。

ステップ４０ｂでは、待ち時間決定手段９３が、マップデータＤ１の中に、ｒ≧０．５（式（６）参照）を満たすコントラストｒが含まれているか否かを判断する。図１７を参照すると、マップデータＤ１はコントラスト０．５の領域を僅かに横切っていることが分かる。したがって、マップデータＤ１の中にｒ≧０．５を満たすコントラストｒが含まれている。図１８に、マップデータＤ１の中に含まれるｒ＝０．５の範囲ｖ１を示す。マップデータＤ１の中にｒ≧０．５を満たすコントラストｒが含まれているので、ステップ４０ｄに進む。

ステップ４０ｄでは、待ち時間決定手段９３が、ｒ＝０．５の範囲ｖ１に対応する待ち時間Ｔｗの範囲を特定する。図１９に、ｒ＝０．５の範囲ｖ１に対応する待ち時間Ｔｗの範囲Ｈ１を示す。ここでは、待ち時間Ｔｗの範囲Ｈ１は、以下の範囲とする。
１０００（ｍｓ）≦Ｔｗ≦１５００（ｍｓ） ・・・（１７）

コントラストマップＭ１における待ち時間Ｔｗの範囲Ｈ１を求めた後、ステップ４０ｅに進む。

ステップ４０ｅでは、待ち時間決定手段９３が、コントラストマップＭ２（図６参照）の中から、ステップＳＴ３で決定された反転時間ＴＩ＝１４００（ｍｓ）におけるマップデータを抽出する。図２０に抽出されたマップデータＤ２を示す。マップデータＤ２を抽出した後、ステップ４０ｆに進む。

ステップ４０ｆでは、待ち時間決定手段９３が、マップデータＤ１の中に、ｓ≧０．５（式（９）参照）を満たすコントラストｓが含まれているか否かを判断する。図２０を参照すると、マップデータＤ２はコントラスト０．５〜０．５４の領域を横切っていることが分かる。したがって、マップデータＤ２の中にｓ≧０．５を満たすコントラストｓが含まれている。図２１に、マップデータＤ２の中に含まれるｓ＝０．５〜０．５４の範囲ｖ２を示す。マップデータＤ２の中にｓ≧０．５を満たすコントラストｓが含まれているので、ステップ４０ｈに進む。

ステップ４０ｈでは、待ち時間決定手段９３が、ｓ＝０．５〜０．５４の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を特定する。図２２に、ｓ＝０．５〜０．５４の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を示す。ここでは、待ち時間Ｔｗの範囲Ｈ２は、以下の範囲とする。
１７００（ｍｓ）≦Ｔｗ≦４０００（ｍｓ） ・・・（１８）

コントラストマップＭ２における待ち時間Ｔｗの範囲Ｈ２を求めた後、ステップ４０ｉに進む。

ステップ４０ｉでは、待ち時間決定手段９３が、コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１（図１９参照）と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２（図２２参照）に、重なる部分があるか否かを判断する。図２３に、コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２とを比較して示す。図２３では、待ち時間Ｔｗが重なる範囲は存在していない。したがって、ステップ４０ｊに進む。

ステップ４０ｊでは、待ち時間決定手段９３は、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２を、待ち時間Ｔｗとして使用可能な時間範囲Ｈと決定する。ここでは、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２は、１７００ｍｓ≦Ｔｗ≦４０００ｍｓであるので、待ち時間Ｔｗとして使用可能な時間範囲Ｈは、以下の範囲となる。
１７００ｍｓ≦Ｔｗ≦４０００ｍｓ ・・・（１９）

待ち時間Ｔｗとして使用可能な時間範囲Ｈを式（１９）の範囲に設定することにより、動脈血と静脈血とのコントラストｓをｓ≧０．５以上にすることができる。尚、待ち時間Ｔｗとして使用可能な時間範囲Ｈを式（１９）の範囲に設定した場合、動脈血とＣＳＦとのコントラストｒはｒ＜０．５になるので（図１９参照）、動脈血とＣＳＦとのコントラストｒが小さくなってしまうことが考えられる。しかし、ＣＳＦの位置は、診断に重要な頸動脈の位置に対してＲＬ方向にずれているので、動脈血とＣＳＦとのコントラストｒが小さくても、頸動脈の診断にはそれほど妨げにならないと考えられる。

待ち時間Ｔｗとして使用可能な時間範囲が決定されたら、ステップＳＴ４０のフローを終了する。ステップＳＴ４０が終了したら、ステップＳＴ４１（図７参照）に進む。

ステップＳＴ４１では、待ち時間決定手段９３が、待ち時間Ｔｗとして使用可能な時間範囲の中から、イメージングシーケンスＰＳを実行するときの待ち時間Ｔｗを決定する。以下に待ち時間Ｔｗの決定方法について説明する。

先に説明したように、ＲＲ＝１０００ｍｓ、Ｔａ＝５００ｍｓとすると、待ち時間Ｔｗは式（１４）で表すことができる。以下に式（１４）を再び示しておく。
Ｔｗ＝１０００ｎ−（ＴＩ＋５００） ・・・（１４）

ここで、反転時間ＴＩは、ステップＳＴ３でＴＩ＝１４００ｍｓと算出されている。したがって、ＴＩ＝１４００ｍｓを式（１４）に代入すると、以下の式が得られる。
Ｔｗ＝１０００ｎ−（ＴＩ＋５００）
＝１０００ｎ−（１４００＋５００）
＝１０００ｎ−１９００ ・・・（２０）

また、ステップＳＴ４０において、待ち時間Ｔｗとして使用可能な時間範囲Ｈは１７００ｍｓ≦Ｔｗ≦４０００ｍｓとして求められている（式（１９）参照）。したがって、式（２０）で表される待ち時間Ｔｗは、以下の式を満たす必要がある。
１７００ｍｓ≦Ｔｗ＝１０００ｎ−１９００≦４０００ｍｓ ・・・（２１）

待ち時間決定手段９３は、式（２１）を満たすｎを求める。式（２１）を満たすｎは、ｎ＝４および５である。ｎ＝４を式（２０）に代入すると、待ち時間Ｔｗは以下の値となる。
Ｔｗ＝１０００ｎ−１９００
＝１０００×４−１９００
＝２１００ｍｓ

一方、ｎ＝５を式（２０）に代入すると、待ち時間Ｔｗは以下の値となる。
Ｔｗ＝１０００ｎ−１９００
＝１０００×５−１９００
＝３１００ｍｓ

したがって、２通りの待ち時間Ｔｗ＝２１００ｍｓ、Ｔｗ＝３１００ｍｓが得られる。このように、式（２１）を満たす待ち時間Ｔｗが複数存在する場合は、スキャン時間が短くなるように、最小の待ち時間Ｔｗ（＝２１００ｍｓ）を、イメージングシーケンスＰＳを実行するときの待ち時間Ｔｗとして決定する。

したがって、ステップＳＴ３およびＳＴ４の処理を実行することによって、ＴＩ＝１４００ｍｓおよびＴｗ＝２１００ｍｓが求められる。図２４に、ＴＩ＝１４００ｍｓおよびＴｗ＝２１００ｍｓの場合のイメージングシーケンスＰＳを示す。尚、ｎ＝４であるので、ｎ＝４を式（２）に代入すると、ＴＲ＝４ＲＲとなる。

待ち時間Ｔｗを決定した後、ステップＳＴ２０（図３参照）に進み、図２４に示すイメージングシーケンスＰＳに従ってｋ空間のデータを収集し、フローを終了する。
最後に、ＴＩ＝１３００ｍｓの場合について説明する。

（３）ステップＳＴ３においてＴＩ＝１３００ｍｓと決定された場合
ステップＳＴ３においてＴＩ＝１３００ｍｓと決定されると、ＴＩ＝１５００ｍｓおよび１４００ｍｓの場合と同様に、ステップ４０ａ（図８参照）に進む。

ステップ４０ａでは、待ち時間決定手段９３が、コントラストマップＭ１（図５参照）の中から、ステップＳＴ３で決定された反転時間ＴＩ＝１３００（ｍｓ）におけるマップデータを抽出する。図２５に抽出されたマップデータＤ１を示す。マップデータＤ１を抽出した後、ステップ４０ｂに進む。

ステップ４０ｂでは、待ち時間決定手段９３が、マップデータＤ１の中に、ｒ≧０．５（式（６）参照）を満たすコントラストｒが含まれているか否かを判断する。図２５を参照すると、マップデータＤ１はコントラスト０．４５の領域を横切っているが、コントラスト０．５の領域は横切っていない。したがって、マップデータＤ１の中にｒ≧０．５を満たすコントラストｒが含まれていない。この場合、ステップ４０ｃに進む。

ステップ４０ｃでは、待ち時間決定手段９３は、先ず、マップデータＤ１の中からコントラストｒの最大値を求める。図２５を参照すると、マップデータＤ１におけるコントラストｒの最大値は、０．４５である。コントラストｒの最大値０．４５を求めた後、待ち時間決定手段９３は、コントラストｒの最大値０．４５に対応する待ち時間Ｔｗの範囲Ｈ１を特定する。図２６に、ｒ＝０．４５に対応する待ち時間Ｔｗの範囲Ｈ１を示す。ここでは、待ち時間Ｔｗの範囲Ｈ１は、以下の範囲とする。
６００（ｍｓ）≦Ｔｗ≦２３００（ｍｓ） ・・・（２２）

コントラストマップＭ１における待ち時間Ｔｗの範囲Ｈ１を求めた後、ステップ４０ｅに進む。

ステップ４０ｅでは、待ち時間決定手段９３が、コントラストマップＭ２（図６参照）の中から、ステップＳＴ３で決定された反転時間ＴＩ＝１３００（ｍｓ）におけるマップデータを抽出する。図２７に抽出されたマップデータＤ２を示す。マップデータＤ２を抽出した後、ステップ４０ｆに進む。

ステップ４０ｆでは、待ち時間決定手段９３が、マップデータＤ２の中に、ｓ≧０．５（式（９）参照）を満たすコントラストｓが含まれているか否かを判断する。図２７を参照すると、マップデータＤ２はコントラスト０．５〜０．５４の領域を横切っていることが分かる。したがって、マップデータＤ２の中にｓ≧０．５を満たすコントラストｒが含まれている。図２８に、マップデータＤ２の中に含まれるｓ＝０．５〜０．５４の範囲ｖ２を示す。マップデータＤ２の中にｓ≧０．５を満たすコントラストｒが含まれているので、ステップ４０ｈに進む。

ステップ４０ｈでは、待ち時間決定手段９３が、ｓ＝０．５〜０．５４の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を特定する。図２９に、ｓ＝０．５〜０．５４の範囲ｖ２に対応する待ち時間Ｔｗの範囲Ｈ２を示す。ここでは、待ち時間Ｔｗの範囲Ｈ２は、以下の範囲とする。
１６００（ｍｓ）≦Ｔｗ≦４０００（ｍｓ） ・・・（２３）

コントラストマップＭ２における待ち時間Ｔｗの範囲Ｈ２を求めた後、ステップ４０ｉに進む。

ステップ４０ｉでは、待ち時間決定手段９３が、コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１（図２６参照）と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２（図２９参照）に、重なる部分があるか否かを判断する。図３０に、コントラストマップＭ１の待ち時間Ｔｗの範囲Ｈ１と、コントラストマップＭ２の待ち時間Ｔｗの範囲Ｈ２とを比較して示す。図３０では、１６００ｍｓ≦Ｔｗ≦２３００ｍｓの範囲で待ち時間Ｔｗが重なっている。したがって、ステップ４０ｋに進む。

ステップ４０ｋでは、待ち時間決定手段９３は、待ち時間Ｔｗの重なっている範囲Ｈを、待ち時間Ｔｗとして使用可能な時間範囲と決定する。ここでは、１６００ｍｓ≦Ｔｗ≦２３００ｍｓの範囲で待ち時間Ｔｗが重なっている。したがって、待ち時間Ｔｗとして使用可能な時間範囲Ｈは、以下の範囲となる。
１６００ｍｓ≦Ｔｗ≦２３００ｍｓ ・・・（２４）

待ち時間Ｔｗとして使用可能な時間範囲が決定されたら、ステップＳＴ４０のフローを終了する。ステップＳＴ４０が終了したら、ステップＳＴ４１（図７参照）に進む。

ステップＳＴ４１では、待ち時間決定手段９３が、待ち時間Ｔｗとして使用可能な時間範囲の中から、イメージングシーケンスＰＳを実行するときの待ち時間Ｔｗを決定する。以下に待ち時間Ｔｗの決定方法について説明する。

先に説明したように、ＲＲ＝１０００ｍｓ、Ｔａ＝５００ｍｓとすると、待ち時間Ｔｗは式（１４）で表すことができる。以下に式（１４）を再び示しておく。
Ｔｗ＝１０００ｎ−（ＴＩ＋５００） ・・・（１４）

ここで、反転時間ＴＩは、ステップＳＴ３でＴＩ＝１３００ｍｓと算出されている。したがって、ＴＩ＝１３００ｍｓを式（１４）に代入すると、以下の式が得られる。
Ｔｗ＝１０００ｎ−（ＴＩ＋５００）
＝１０００ｎ−（１３００＋５００）
＝１０００ｎ−１８００ ・・・（２５）

また、ステップＳＴ４０において、待ち時間Ｔｗとして使用可能な時間範囲Ｈは１６００ｍｓ≦Ｔｗ≦２３００ｍｓとして求められている（式（２４）参照）。したがって、式（２５）で表される待ち時間Ｔｗは、以下の式を満たす必要がある。
１６００ｍｓ≦Ｔｗ＝１０００ｎ−１８００≦２３００ｍｓ ・・・（２６）

待ち時間決定手段９３は、式（２６）を満たすｎを求める。式（２６）を満たすｎは、ｎ＝４および５である。ｎ＝４を式（２５）に代入すると、待ち時間Ｔｗは以下の値となる。
Ｔｗ＝１０００ｎ−１８００
＝１０００×４−１８００
＝２２００ｍｓ

一方、ｎ＝５を式（２５）に代入すると、待ち時間Ｔｗは以下の値となる。
Ｔｗ＝１０００ｎ−１８００
＝１０００×５−１８００
＝３２００ｍｓ

したがって、２通りの待ち時間Ｔｗ＝２２００ｍｓ、Ｔｗ＝３２００ｍｓが得られる。このように、式（２６）を満たす待ち時間Ｔｗが複数存在する場合は、スキャン時間が短くなるように、最小の待ち時間Ｔｗ（＝２２００ｍｓ）を、イメージングシーケンスＰＳを実行するときの待ち時間Ｔｗとして決定する。

したがって、ステップＳＴ３およびＳＴ４の処理を実行することによって、ＴＩ＝１３００ｍｓおよびＴｗ＝２２００ｍｓが求められる。図３１に、ＴＩ＝１３００ｍｓおよびＴｗ＝２２００ｍｓの場合のイメージングシーケンスＰＳを示す。尚、ｎ＝４であるので、ｎ＝４を式（２）に代入すると、ＴＲ＝４ＲＲとなる。

待ち時間Ｔｗを決定した後、ステップＳＴ２０（図３参照）に進み、図３１に示すイメージングシーケンスＰＳに従ってｋ空間のデータを収集し、フローを終了する。

ＴＩ＝１３００ｍｓの場合、コントラストマップＭ１から抽出されるマップデータＤ１には、ｒ≧０．５のコントラストは含まれていない。そこで、マップデータＤ１の中からコントラストｒの最大値０．４５を求め、コントラストの最大値０．４５に対応する待ち時間Ｔｗの範囲Ｈ１（６００ｍｓ≦Ｔｗ≦２３００ｍｓ）を特定している（図２６参照）。その後、待ち時間Ｔｗの重なる範囲Ｈを求め（図３０参照）、この範囲Ｈの中から、待ち時間Ｔｗを決定している。したがって、コントラストｒの値を、ＴＩ＝１３００ｍｓのときに取り得る最大値に設定することができる。

尚、本形態では、イメージングシーケンスＰＳは反転パルスＳＩＲ（つまり、フリップ角が１８０°のＲＦパルス）を有している。しかし、本発明は反転パルスＳＩＲに限定されることはなく、α°パルス（αは任意の角度）を用いることができる。図３２に、α°パルスを用いた場合のイメージングシーケンスＰＳの一例を示す。図３２では、α°パルスとデータ収集シーケンスＤＡＱとの間の時間は「Ｔα」で示されている。この場合、時間Ｔαは、反転時間ＴＩと同様に、動脈血の流速ｖに基づいて算出すればよい。また、待ち時間Ｔｗは、時間Ｔαと待ち時間Ｔｗとコントラストｒ（又はコントラストｓ）との対応関係を規定するコントラストマップを用いて決定することができる。

また、時間Ｔαは、図３３に示すように、α°パルスとデータ収集シーケンスＤＡＱの終了時点ｔｅとの間の時間としてもよいし、図３４に示すように、α°パルスと、データ収集シーケンスＤＡＱが実行されている途中の時点ｔｍとの間の時間としてもよい。

尚、本形態では、２つのコントラストマップＭ１およびＭ２を用いて待ち時間Ｔｗを決定したが、いずれか一方のコントラストマップのみを用いて待ち時間Ｔｗを決定してもよい。例えば、動脈血と静脈血とのコントラストを大きくすることが重要な場合は、コントラストマップＭ２のみを用いて待ち時間Ｔｗを決定してもよいし、一方、動脈血とＣＳＦとのコントラストを大きくすることが重要な場合は、コントラストマップＭ１のみを用いて待ち時間Ｔｗを決定してもよい。

また、本形態では、イメージングシーケンスＰＳは、選択反転パルスＳＩＲと、脂肪抑制パルスＦと、データ収集シーケンスＤＡＱと、非選択反転パルスＮＩＲとを備えている。しかし、イメージングシーケンスＰＳは、選択反転パルスＳＩＲ、脂肪抑制パルスＦ、および非選択反転パルスＮＩＲとは別のパルスを有していてもよいし、データ収集シーケンスＤＡＱとは別のシーケンスを有していてもよい。また、イメージングシーケンスＰＳは、脂肪抑制パルスＦおよび非選択反転パルスＮＩＲを備えているが、これらのパルスＦおよびＮＩＲを備えないようにしてもよい。

本形態では、動脈血を撮影する例について説明したが、本発明は、静脈血を撮影する場合にも適用することができる。

２ マグネット
３ テーブル
３ａ クレードル
４ 受信コイル
５ 送信器
６ 勾配磁場電源
７ 受信器
８ 制御部
８ 操作部
１０ 表示部
１１ 被検体
２１ ボア
８０ 記憶部
９１ 流速算出手段
９２ 反転時間決定手段
９３ 待ち時間決定手段
１００ ＭＲ装置

Claims (16)

1. 血液を含む領域内のスピンをフリップさせるための第１のＲＦパルスと前記領域から血液のデータを収集するデータ収集シーケンスとを有するイメージングシーケンスを繰り返し実行する磁気共鳴装置であって、
   前記血液と背景組織とのコントラストと、前記第１のＲＦパルスから前記データ収集シーケンスまでの第１の時間と、イメージングシーケンスが終了してから次のイメージングシーケンスが開始されるまでの第２の時間との対応関係を記憶する記憶部と、
   前記血液の流速に基づいて、前記イメージングシーケンスを繰り返し実行するときの前記第１の時間を決定する第１の決定手段と、
   前記第１の決定手段により決定された前記第１の時間と、前記対応関係とに基づいて、前記イメージングシーケンスを繰り返し実行するときの前記第２の時間を決定する第２の決定手段と、
   を有する、磁気共鳴装置。
2. 前記記憶部は、
   前記血液と第１の背景組織とのコントラストと、前記第１のＲＦパルスから前記データ収集シーケンスまでの第１の時間と、イメージングシーケンスが終了してから次のイメージングシーケンスが開始されるまでの第２の時間との対応関係を規定する第１のマップを記憶し、
   前記第２の決定手段は、
   前記第１の決定手段により決定された前記第１の時間と前記第１のマップと基づいて、前記第２の時間として使用可能な時間範囲を決定し、前記時間範囲の中から、前記イメージングシーケンスを繰り返し実行するときの前記第２の時間を決定する、請求項１に記載の磁気共鳴装置。
3. 前記血液は動脈血であり、前記第１の背景組織は静脈血である、請求項２に記載の磁気共鳴装置。
4. 前記第２の決定手段は、
   前記第１のマップの中から、前記第１の決定手段により決定された前記第１の時間における第１のデータを特定し、前記第１のデータに基づいて、前記第２の時間として使用可能な時間範囲を決定する、請求項２又は３に記載の磁気共鳴装置。
5. 前記記憶部は、
   前記血液と第２の背景組織とのコントラストと、前記第１のＲＦパルスから前記データ収集シーケンスまでの第１の時間と、イメージングシーケンスが終了してから次のイメージングシーケンスが開始されるまでの第２の時間との第２の対応関係を規定する第２のマップを記憶し、
   前記第２の決定手段は、
   前記第２のマップの中から、前記第１の決定手段により決定された前記第１の時間における第２のデータを特定し、前記第１のデータと前記第２のデータとに基づいて、前記第２の時間として使用可能な時間範囲を決定する、請求項４に記載の磁気共鳴装置。
6. 前記第２の決定手段は、
   前記第１のデータに基づいて、第１の所定値以上のコントラストに対応する前記第２の時間の第１の範囲を特定し、
   前記第２のデータに基づいて、第２の所定値以上のコントラストに対応する前記第２の時間の第２の範囲を特定し、
   前記第２の時間の第１の範囲と前記第２の時間の第２の範囲とに基づいて、前記第２の時間として使用可能な時間範囲を決定する、請求項５に記載の磁気共鳴装置。
7. 前記第２の決定手段は、
   前記第２の時間の第１の範囲と前記第２の時間の第２の範囲との間に、前記第２の時間の重なる範囲が存在する場合、前記第２の時間の重なる範囲を、前記第２の時間として使用可能な時間範囲と決定する、請求項６に記載の磁気共鳴装置。
8. 前記第２の決定手段は、
   前記第２の時間の第１の範囲と前記第２の時間の第２の範囲との間に、前記第２の時間の重なる範囲が存在しない場合、前記第１の範囲を、前記第２の時間として使用可能な時間範囲と決定する、請求項６に記載の磁気共鳴装置。
9. 前記第２の背景組織は脳脊髄液である、請求項５〜８のうちのいずれか一項に記載の磁気共鳴装置。
10. 前記第２の決定手段は、
    心拍数又はＲＲ間隔と、前記第１の決定手段により決定された前記第１の時間とに基づいて、前記第２の時間として使用可能な時間範囲の中から、前記イメージングシーケンスを繰り返し実行するときの前記第２の時間を決定する、請求項２〜９のうちのいずれか一項に記載の磁気共鳴装置。
11. 前記第１の時間は、前記第１のＲＦパルスから前記データ収集シーケンスの開始時点までの時間である、請求項１〜１０のうちのいずれか一項に記載の磁気共鳴装置。
12. 前記第１のＲＦパルスは反転パルスであり、前記第１の時間は反転時間である、請求項１１に記載の磁気共鳴装置。
13. 前記イメージングシーケンスは、前記第１のＲＦパルスと前記データ収集シーケンスとの間に、第２のＲＦパルスを有する、請求項１〜１２のうちのいずれか一項に記載の磁気共鳴装置。
14. 前記第２のＲＦパルスは脂肪抑制パルスである、請求項１３に記載の磁気共鳴装置。
15. 前記血液の流速を求めるためのシーケンスを実行する、請求項１〜１４のうちのいずれか一項に記載の磁気共鳴装置。
16. 血液を含む領域内のスピンをフリップさせるための第１のＲＦパルスと前記領域から血液のデータを収集するデータ収集シーケンスとを有するイメージングシーケンスを繰り返し実行する磁気共鳴装置であって、前記血液と背景組織とのコントラストと、前記第１のＲＦパルスから前記データ収集シーケンスまでの第１の時間と、イメージングシーケンスが終了してから次のイメージングシーケンスが開始されるまでの第２の時間との対応関係を記憶する磁気共鳴装置に適用されるプログラムであって、
    前記血液の流速に基づいて、前記イメージングシーケンスを繰り返し実行するときの前記第１の時間を決定する第１の決定処理と、
    前記第１の決定処理により決定された前記第１の時間と、前記対応関係とに基づいて、前記イメージングシーケンスを繰り返し実行するときの前記第２の時間を決定する第２の決定処理と、
    を計算機に実行させるためのプログラム。