(57)【要約】      （修正有）【課題】より正確な静脈出力関数を算出することが可能な血流動態装置、その血流動態解析装置を有する磁気共鳴イメージング装置、および被検体の血流動態を解析する血流動態解析方法を提供する。【解決手段】各スライスＳ１〜Ｓｎごとに得られたピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップから、ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃｐのジョイントヒストグラムを作成し、作成したジョイントヒストグラムの静脈範囲Ａβにおいて、ボーラス通過後濃度Ｃｐが負の値を有するドットの中から、ピーク濃度Ｃmaxの大きい順に３個のドットＤ２１、Ｄ２２、およびＤ２３を選択する。【選択図】図１４
2010104671.tif 000002
【特許請求の範囲】
【請求項１】
造影剤が注入された被検体に設定された各スライス内の複数の領域からデータを収集し、前記データに基づいて、前記被検体の血流動態を解析する血流動態解析装置であって、
前記各スライス内の前記領域ごとに、前記造影剤の濃度の時間変化を表すプロファイルを作成する造影剤濃度プロファイル作成手段と、
前記プロファイルから、前記造影剤のボーラスが前記領域を通過した後の前記領域における前記造影剤の濃度を算出する特徴量算出手段と、
前記造影剤のボーラスが前記領域を通過した後の前記領域における前記造影剤の濃度を含むヒストグラムを作成するヒストグラム作成手段と、
前記ヒストグラムに基づいて、前記複数の領域の中から、静脈が存在する領域を決定する領域決定手段と、
前記静脈が存在する領域における前記造影剤の濃度の時間変化に基づいて、静脈出力関数を算出する関数算出手段と、
を有する、血流動態解析装置。
【請求項２】
前記特徴量算出手段は、
前記造影剤濃度プロファイルから、前記造影剤の濃度がピークになるときのピーク濃度と、前記造影剤の濃度が前記ピーク濃度に到達するまでに掛かるピーク濃度到達時間とを算出する、請求項１に記載の血流動態解析装置。
【請求項３】
前記ヒストグラム作成手段は、
前記造影剤のボーラスが前記領域を通過した後の前記領域における前記造影剤の濃度と、前記ピーク濃度と、前記ピーク濃度到達時間とを含むヒストグラムを作成する、請求項２に記載の血流動態解析装置。
【請求項４】
前記領域決定手段は、
前記ヒストグラムに基づいて、前記複数の領域の中から、第１の範囲内のピーク濃度到達時間を有する複数の第１の領域を抽出する領域抽出部と、
前記複数の第１の領域の中から、ノイズとなる領域を除去するノイズ除去部と、
前記造影剤のボーラスが前記領域を通過した後の前記領域における前記造影剤の濃度に基づいて、前記ノイズとなる領域が除去された前記複数の第１の領域の中から、静脈が存在する領域を選択する領域選択部と、
を有する請求項３に記載の血流動態解析装置。
【請求項５】
前記領域選択部は、
前記造影剤のボーラスが前記領域を通過した後の前記領域における前記造影剤の濃度と、前記ピーク濃度とに基づいて、前記ノイズとなる領域が除去された前記複数の第１の領域の中から、静脈が存在する領域を選択する、請求項４に記載の血流動態解析装置。
【請求項６】
前記領域抽出部は、前記ヒストグラムに基づいて、前記複数の領域の中から、第２の範囲内のピーク濃度到達時間を有する複数の第２の領域を抽出し、
前記ノイズ除去部は、前記複数の第２の領域の中から、ノイズとなる領域を除去し、
前記領域選択部は、前記ピーク濃度に基づいて、前記ノイズとなる領域が除去された前記複数の第２の領域の中から、動脈が存在する領域を選択する、請求項４又は５に記載の血流動態解析装置。
【請求項７】
前記特徴量算出手段は、
前記造影剤のボーラスが前記領域を通過した後の前記領域における前記造影剤の濃度を、前記造影剤の濃度がピークになった時刻から、所定時間が経過した後の造影剤の濃度に基づいて算出する、請求項１〜６のうちのいずれか一項に記載の血流動態解析装置。
【請求項８】
前記所定時間は、５秒〜１０秒の時間である、請求項７に記載の血流動態解析装置。
【請求項９】
前記データは、磁気共鳴イメージング装置によって収集されたデータである、請求項１〜８のうちのいずれか一項に記載の血流動態解析装置。
【請求項１０】
請求項１〜９のうちのいずれか一項に記載の血流動態解析装置を有する磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体の血流動態を解析する血流動態解析装置、およびその血流動態解析装置を有する磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
脳梗塞を発症した患者に対しては、迅速に適切な処置をする必要がある。迅速な処理をするためには、被検体の脳の血流状態を定量化する重要性が高まっている。そこで、造影剤を使って被検体の脳を断層撮影し、脳の血流状態を定量化する種々の方法が提案されている（特許文献１参照）。
【特許文献１】特開2006-326078号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
脳の血流状態を定量化する場合、動脈入力関数（ＡＩＦ：Arterial Input Function）および静脈出力関数（ＶＯＦ：Venous Output Function）を算出する必要がある。動脈入力関数は、被検体の断層像の動脈における造影剤の濃度の時間変化から算出される関数であり、静脈出力関数は、被検体の断層像の静脈における造影剤の濃度の時間変化から算出される関数である。しかし、被検体の脳の動脈の内径が、病巣などが原因で細くなっている場合、動脈における造影剤の濃度の時間変化が、静脈における造影剤の濃度の時間変化に類似した挙動を示すことがある。この場合、従来の方法では、動脈における造影剤の濃度の時間変化に基づいて、静脈出力関数を算出してしまうことがあり、正確な静脈出力関数を算出することが困難となる。
【０００４】
本発明は、上記の事情に鑑み、より正確な静脈出力関数を算出することが可能な血流動態装置、その血流動態解析装置を有する磁気共鳴イメージング装置、および被検体の血流動態を解析する血流動態解析方法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
上記の問題を解決する本発明の血流動態解析装置は、
造影剤が注入された被検体に設定された各スライス内の複数の領域からデータを収集し、上記データに基づいて、上記被検体の血流動態を解析する血流動態解析装置であって、
上記各スライス内の上記領域ごとに、上記造影剤の濃度の時間変化を表すプロファイルを作成する造影剤濃度プロファイル作成手段と、
上記プロファイルから、上記造影剤のボーラスが上記領域を通過した後の上記領域における上記造影剤の濃度を算出する特徴量算出手段と、
上記造影剤のボーラスが上記領域を通過した後の上記領域における上記造影剤の濃度を含むヒストグラムを作成するヒストグラム作成手段と、
上記ヒストグラムに基づいて、上記複数の領域の中から、静脈が存在する領域を決定する領域決定手段と、
上記静脈が存在する領域における上記造影剤の濃度の時間変化に基づいて、静脈出力関数を算出する出力関数算出手段と、
を有している。
【発明の効果】
【０００６】
本発明では、各スライス内の領域ごとに作成された造影剤濃度プロファイルから、各スライス内の各領域ごとに、造影剤のボーラス（Bolus）が通過した後の造影剤の濃度を算出する。造影剤のボーラスが通過した後の造影剤の濃度は、造影剤のボーラスが動脈を流れるのか、静脈を流れるのかによって、異なる値になる傾向がある。したがって、造影剤のボーラスが通過した後の造影剤の濃度を含むヒストグラムを作成することにより、各スライス内において、静脈を含む領域と、動脈を含む領域とを区別することが可能になり、より正確な静脈出力関数を算出することができる。
【発明を実施するための最良の形態】
【０００７】
以下、図面を参照しながら、発明を実施するための最良の形態を詳細に説明する。
【０００８】
図１は、本発明の一実施形態の磁気共鳴イメージング装置１の概略図である。
【０００９】
磁気共鳴イメージング装置（以下、ＭＲＩ（Magnetic
Resonance Imaging）装置と呼ぶ）１は、ガントリ２と、クレードル３と、受信コイル４と、造影剤注入装置５と、制御装置６と、入力装置７とを有している。
【００１０】
ガントリ２は、被検体８が収容されるボア２１と、超伝導コイル２２と、勾配コイル２３と、送信コイル２４とを有している。超伝導コイル２２は静磁場Ｂ0を印加し、勾配コイル２３は勾配パルスを印加し、送信コイル２４はＲＦパルスを送信する。
【００１１】
クレードル３は、ｚ方向および−ｚ方向に移動するように構成されている。クレードル３がｚ方向に移動することによって、被検体８がボア２１に搬送される。クレードル３が−ｚ方向に移動することによって、ボア２１に搬送された被検体８は、ボア２１から搬出される。
【００１２】
造影剤注入装置５は、被検体８に造影剤を注入する。
【００１３】
受信コイル４は、被検体８の頭部８ａに取り付けられている。受信コイル４が受信したＭＲ（Magnetic Resonance）信号は、制御装置６に伝送される。
【００１４】
制御装置６は、コイル制御手段６１〜関数算出手段７０を有している。
【００１５】
コイル制御手段６１は、オペレータ９によって入力装置７から入力された被検体８の撮影命令に応答して、被検体８を撮影するためのパルスシーケンスが実行されるように、送信コイル２４および勾配コイル２３を制御する。
【００１６】
造影剤濃度プロファイル作成手段６２は、スライスＳ１〜Ｓｎの領域ごとに、被検体８の頭部８ａの造影剤の濃度の時間変化を表す造影剤濃度−時間曲線ＧａおよびＧｖを算出する（図５参照）。
【００１７】
特徴量算出手段６３は、造影剤濃度−時間曲線ＧａおよびＧｖ（図５（ｂ）および（ｃ）参照）の３つの特徴量（ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃp）を算出する。
【００１８】
マップ作成手段６４は、各スライスＳ１〜Ｓｎについて、ピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップを作成する（図６および図７参照）。
【００１９】
ジョイントヒストグラム作成手段６５は、各スライスＳ１〜Ｓｎごとに得られたピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップから（図７参照）、ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃｐのジョイントヒストグラム（図８参照）を作成する。
【００２０】
領域決定手段６６は、ジョイントヒストグラム（図８参照）に基づいて、スライスＳ１〜Ｓｎが有する全ての領域Ｒの中から、静脈が存在する領域と、動脈が存在する領域とを決定する。領域決定手段６６は、領域抽出部６７と、ノイズ除去部６８と、領域選択部６９とを有している。
【００２１】
領域抽出部６７は、ジョイントヒストグラム（図８参照）に基づいて、スライスＳ１〜Ｓｎが有する全ての領域Ｒの中から、ｔ１〜ｔ２の範囲内のピーク濃度到達時間（図１１参照）を有する複数の領域Ｒ１およびＲ２（図１２（ｃ）参照）と、ｔ３〜ｔ４の範囲内のピーク濃度到達時間（図１１参照）を有する複数の領域Ｒ３およびＲ４（図１２（ｃ）参照）とを抽出する。
【００２２】
ノイズ除去部６８は、複数の領域Ｒ１およびＲ２（図１２（ｃ）参照）の中から、ノイズとなる領域Ｒ２を除去するとともに、複数の領域Ｒ３およびＲ４（図１２（ｃ）参照）の中から、ノイズとなる領域Ｒ４を除去する。
【００２３】
領域選択部６９は、ボーラス通過後濃度Ｃｐに基づいて、複数の領域Ｒ１の中から、動脈が存在する領域Ｒ１１、Ｒ１２、およびＲ１３（図１４（ｃ）参照）を選択する。更に、領域選択部６９は、ボーラス通過後濃度Ｃｐに基づいて、複数の領域Ｒ３の中から、静脈が存在する領域Ｒ２１、Ｒ２２、およびＲ２３（図１４（ｃ）参照）を選択する。
【００２４】
関数算出手段７０は、選択された３つの領域Ｒ１１、Ｒ１２、およびＲ１３における造影剤濃度−時間曲線Ｃ１１、Ｃ１２、およびＣ１３を加算平均することにより、動脈入力関数Ｆａを算出する（図１６参照）。更に、関数算出手段７０は、選択された３つの領域Ｒ２１、Ｒ２２、およびＲ２３における造影剤濃度−時間曲線Ｃ２１、Ｃ２２、およびＣ２３を加算平均することにより、静脈出力関数Ｆｖを算出する（図１６参照）。
【００２５】
入力装置７は、オペレータ９の操作に応じて、種々の命令を制御装置６に入力する。
【００２６】
図２は、磁気共鳴イメージング装置１の処理フローを示す図である。
【００２７】
ステップＳ１では、被検体８の頭部８ａの造影撮影が行われる。オペレータは、入力装置７を操作して、被検体８にスライスを設定する。
【００２８】
図３は、被検体８に設定されたスライスの一例である。
【００２９】
被検体８には、ｎ枚のスライスＳ１〜Ｓｎが設定されている。スライスの枚数は、例えば、ｎ＝１２である。スライスの枚数は、必要に応じて、任意の枚数を設定することができる。各スライスＳ１〜Ｓｎごとに、被検体８の頭部８ａの撮影領域が決定される。
【００３０】
オペレータ９は、スライスＳ１〜Ｓｎを設定した後、造影剤注入装置５に造影剤の注入命令を伝送するとともに、ＭＲＩ装置のコイル制御手段６１（図１参照）に、被検体８を撮影する撮影命令を伝送する。コイル制御手段６１は、撮影命令に応答して、被検体８の頭部８ａを撮影するためのパルスシーケンスが実行されるように、送信コイル２４および勾配コイル２３を制御する。
【００３１】
本実施形態では、マルチスライススキャンにより、各スライスから、連続撮影されたｍ枚のフレーム画像を得るためのパルスシーケンスが実行される。したがって、１枚のスライスにつき、ｍ枚のフレーム画像が得られる。例えば、フレーム画像の枚数ｍ＝８５枚である。パルスシーケンスを実行することにより、被検体８の頭部８ａからデータが収集される。
【００３２】
図４は、スライスＳ１〜ＳｎのうちのスライスＳｘから得られるフレーム画像を示す概念図である。
【００３３】
スライスＳｘの断面からは、ｍ枚のフレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）が得られる。各フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）は、α×β個の画素Ｐを有している。図４には、各フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）の中の２つの画素ＰｉおよびＰｊの位置が図示されている。
【００３４】
尚、図４では、スライスＳｘで得られるフレーム画像のみが示されているが、他のスライスでも、スライスＳｘと同様に、ｍ枚のフレーム画像が得られる。
【００３５】
ステップＳ１を実行した後、ステップＳ２に進む。
【００３６】
ステップＳ２では、造影剤濃度プロファイル作成手段６２（図１参照）が、スライスＳ１〜Ｓｎの断面領域において、被検体８の頭部８ａの造影剤の濃度の時間変化を算出する。以下に、造影剤濃度プロファイル作成手段６２が、造影剤の濃度の時間変化をどのようにして算出するかについて、図５を参照しながら説明する。
【００３７】
図５は、被検体８の頭部８ａに設定されたスライスＳxの断面領域において、造影剤の濃度の時間変化を算出する方法を説明する図である。
【００３８】
図５（ａ）には、被検体８のスライスＳｘの断面（図４参照）と、スライスＳｘのフレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）が示されている。各フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）は、時刻ｔ１、ｔ２、・・・ｔｍ（時間間隔Δｔ）におけるスライスＳｘの画像である。
【００３９】
スライスＳｘの断面は、α×β個の領域Ｒに分割されている。図５（ａ）には、スライスＳｘの断面の中の２つの領域ＲａおよびＲｖの位置が具体的に図示されている。領域Ｒａは動脈を含む領域であり、一方、領域Ｒｖは静脈を含む領域である。フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）の各画素Ｐは、時刻ｔ１、ｔ２、・・・ｔｍにおけるスライスＳｘの各領域Ｒを撮影したものである。例えば、フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）の画素Ｐａは、時刻ｔ１、ｔ２、・・・ｔｍにおけるスライスＳｘの領域Ｒａを撮影したものである。また、フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）の画素Ｐｖは、時刻ｔ１、ｔ２、・・・ｔｍにおけるスライスＳｘの領域Ｒｖを撮影したものである。
【００４０】
図５（ｂ）には、スライスＳｘの領域Ｒａにおける造影剤濃度−時間曲線（造影剤濃度プロファイル）Ｇａの概略図が示されている。横軸は時間ｔであり、縦軸はスライスＳｘの領域Ｒａにおける造影剤の濃度値Ｃｍ（ｔ）である。図５（ｂ）では、スライスＳｘの撮影開始時刻を、ｔ＝０に定めている。造影剤濃度−時間曲線Ｇａは、各フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）の画素Ｐａにおける造影剤の濃度値Ｃｍ（ｔ）のプロットを線で結ぶことによって作成されている。
【００４１】
図５（ｃ）には、スライスＳｘの領域Ｒｖにおける造影剤濃度−時間曲線（造影剤濃度プロファイル）Ｇｖの概略図が示されている。横軸は時間ｔであり、縦軸はスライスＳｘの領域Ｒｖにおける造影剤の濃度値Ｃｍ（ｔ）である。図５（ｃ）でも、図５（ｂ）と同様に、スライスＳｘの撮影開始時刻を、ｔ＝０に定めている。造影剤濃度−時間曲線Ｇｖは、各フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）の画素Ｐｖにおける造影剤の濃度値Ｃｍ（ｔ）のプロットを線で結ぶことによって作成されている。
【００４２】
動脈を含む領域Ｒａにおける造影剤濃度−時間曲線Ｇａ（図５（ｂ）参照）では、フレーム画像Ｉ（ｘ，１）〜Ｉ（ｘ，ｍ）のうちのどのフレーム画像であっても、造影剤の濃度Ｃｍ（ｔ）は、正の値である。しかし、静脈を含む領域Ｒｖにおける造影剤濃度−時間曲線Ｇｖ（図５（ｃ）参照）では、フレーム画像Ｉ（ｘ，ｙ）〜Ｉ（ｘ，ｍ）において、造影剤の濃度Ｃｍ（ｔ）は、負の値になっている。一般的に、ＭＲＩ装置で取得されたフレーム画像（ｘ，１）〜Ｉ（ｘ，ｍ）のデータに基づいて造影剤の濃度Ｃｍ（ｔ）を算出すると、動脈を含む領域Ｒａでは、造影剤の濃度Ｃｍ（ｔ）は正の値として算出される傾向があるが、静脈を含む領域Ｒｖでは、遅い時刻に取得されたフレーム画像Ｉ（ｘ，ｙ）〜Ｉ（ｘ，ｍ）における造影剤の濃度Ｃｍ（ｔ）は、負の値として算出される傾向があることが、経験的に知られている。
【００４３】
図５には、スライスＳｘの２つの領域ＲａおよびＲｖにおける造影剤濃度−時間曲線が示されているが、スライスＳｘ内の他の領域についても、造影剤濃度−時間曲線が作成される。更に、スライスＳｘ以外の他のスライスについての各領域についても、同様に、造影剤濃度−時間曲線が作成される。
【００４４】
造影剤濃度−時間曲線を作成した後、ステップＳ３に進む。
【００４５】
ステップＳ３では、特徴量算出手段６３（図１参照）が、造影剤濃度−時間曲線ＧａおよびＧｖ（図５（ｂ）および（ｃ）参照）の特徴量を算出する。本実施形態では、造影剤濃度−時間曲線ＧａおよびＧｖの特徴量として、図５（ｂ）および（ｃ）に示すように、３つの特徴量（ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃp）が算出される。
【００４６】
ピーク濃度Ｃmaxは、造影剤の濃度Ｃｍ（ｔ）のピーク値を表している。
【００４７】
ピーク濃度到達時刻ＴＴＰは、撮影開始時刻（ｔ＝０）からピーク濃度Ｃmaxに到達するまでの時間を表している。ただし、ピーク濃度到達時間ＴＴＰの開始時刻は、他の時刻（例えば、造影剤を投与した時刻）とすることもできる。
【００４８】
ボーラス通過後濃度Ｃpは、造影剤のボーラスが各スライスの各領域Ｒを通過した後の造影剤の濃度である。本実施形態では、造影剤のボーラスが各領域Ｒを通過した時刻を、ピーク濃度Ｃmaxに到達した時刻から所定時間Ｔｓが経過した時刻に設定している。所定時間Ｔｓは、経験的に決定することができる値であり、例えば、５秒〜１０秒程度の値とすることができるが、造影剤の種類やその他の条件に応じて、適宜決定すればよい。また、本実施形態では、ボーラス通過後濃度Ｃpを、所定時間Ｔｓが経過した後の複数枚（例えば５枚）のフレーム画像Ｉｆにおける造影剤の濃度Ｃｍ（ｔ）の平均値としている。しかし、所定時間Ｔｓが経過した直後の１枚のフレーム画像における造影剤の濃度Ｃｍ（ｔ）を、ボーラス通過後濃度Ｃｐと定めてもよい。
【００４９】
静脈を含む領域Ｒｖにおけるピーク濃度Ｃmax（図５（ｃ）参照）は、動脈を含む領域Ｒａにおけるピーク濃度Ｃmax（図５（ｂ）参照）よりも時間的に遅れて現れる特性がある。また、上記のように、動脈を含む領域Ｒａでは、造影剤の濃度Ｃｍ（ｔ）は正の値として算出される傾向があるので、動脈を含む領域Ｒａにおけるボーラス通過後濃度Ｃｐも、正の値として算出される傾向がある（図５（ｂ）参照）。しかし、静脈を含む領域Ｒｖでは、遅い時刻に取得されたフレーム画像Ｉ（ｘ，ｙ）〜Ｉ（ｘ，ｍ）における造影剤の濃度Ｃｍ（ｔ）は、負の値として算出される傾向があるので、静脈を含む領域Ｒｖにおけるボーラス通過後濃度Ｃｐは、負の値になる傾向がある（図５（ｃ）参照）。
【００５０】
尚、図５には、スライスＳｘの２つの領域Ｒａおよび領域Ｒｖにおけるピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃｐが示されているが、他の領域Ｒについても、同様に、ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃｐが算出される。更に、スライスＳｘ以外のほかのスライスの各領域についても、同様に、ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃｐが算出される。
【００５１】
上記のように、ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃｐを算出した後、ステップＳ４に進む。
【００５２】
ステップＳ４では、マップ作成手段６４（図１参照）が、各スライスＳ１〜Ｓｎについて、ピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップを作成する（図６および図７参照）。
【００５３】
図６は、スライスＳｘについてのピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップを概略的に示す図である。
【００５４】
図６（ａ）は、スライスＳｘについてのピーク濃度Ｃmaxのマップ、図６（ｂ）は、スライスＳｘについてのピーク濃度到達時間ＴＴＰのマップ、および図６（ｃ）は、スライスＳｘについてのボーラス通過後濃度Ｃｐのマップである。
【００５５】
図６（ａ）には、スライスＳｘの中の領域Ｒによるピーク濃度Ｃmaxの大きさの違いが、領域Ｒの濃淡で示されている。ピーク濃度Ｃmaxが大きければ大きいほど、領域Ｒの色は白に近くなり、ピーク濃度Ｃmaxが小さければ小さいほど、領域Ｒの色は黒に近くなる。図６（ａ）では、説明の便宜上、領域Ｒの濃度の違いを３段階（白、灰色、黒）でしか示していないが、実際は、領域Ｒの濃度の違いは、もっと細かい多段階で示されている。
【００５６】
図６（ｂ）には、スライスＳｘの中の領域Ｒによるピーク濃度到達時間ＴＴＰの大きさの違いが、領域Ｒの濃淡で示されている。ピーク濃度到達時間ＴＴＰが大きければ大きいほど、領域Ｒの色は白に近くなり、ピーク濃度到達時間ＴＴＰが大きければ大きいほど、領域Ｒの色は黒に近くなる。図６（ｂ）では、説明の便宜上、図６（ａ）と同様に、領域Ｒの濃度の違いを３段階（白、灰色、黒）で示してある。
【００５７】
図６（ｃ）には、スライスＳｘの中の領域Ｒによるボーラス通過後濃度Ｃｐの大きさの違いが、領域Ｒの濃淡で示されている。ボーラス通過後濃度Ｃｐが大きければ大きいほど、領域Ｒの色は白に近くなり、ボーラス通過後濃度Ｃｐが小さければ小さいほど、領域Ｒの色は黒に近くなる。図６（ｃ）では、説明の便宜上、図６（ａ）および（ｂ）と同様に、領域Ｒの濃度の違いを３段階（白、灰色、黒）で示してある。
【００５８】
図６に示すようなマップを作成することにより、スライスＳｘの中の領域Ｒによるピーク濃度Ｃmaxの違い、ピーク濃度到達時間ＴＴＰの違い、およびボーラス通過後濃度Ｃｐの違いを理解することができる。例えば、領域Ｒｅの色は、図６（ａ）〜（ｃ）のどのマップであっても白色であるので、Ｃmax、ＴＴＰ、およびＣｐの値は大きいことがわかる。一方、領域Ｒｆの色は、図６（ａ）〜（ｃ）のどのマップであっても灰色であるので、領域ＲｅにおけるＣmax、ＴＴＰ、およびＣｐの値よりも小さいことがわかる。尚、図６のマップには、被検体の頭部の外側にも、濃度の違う領域Ｒが示されているが、これらは、マップを作成するときに現れるノイズである。
【００５９】
図６には、スライスＳｘについての３つのマップ（ピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップ）が示されているが、他のスライスについても、スライスＳｘと同様に、３つのマップ（ピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップ）が作成される（図７参照）。
【００６０】
図７は、各スライスＳ１〜Ｓｎの３つのマップ（ピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップ）を概略的に示す図である。
【００６１】
図７に示すように、ステップＳ４において、スライスＳ１〜Ｓｎごとに、ピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップが作成される。これらのマップを作成した後、ステップＳ５に進む。
【００６２】
ステップＳ５では、ジョイントヒストグラム作成手段６５（図１参照）が、スライスＳ１〜Ｓｎごとに得られたピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップから（図７参照）、ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃｐのジョイントヒストグラムを作成する。
【００６３】
図８は、ジョイントヒストグラムの一例を概略的に示す図である。
【００６４】
ジョイントヒストグラムの中には、多数のドット（白丸）が示されている。各ドットが、スライスＳ１〜Ｓｎの中の各領域Ｒにおける３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係を表している。例えば、ドットＤｅは、スライスＳｘの中の領域Ｒｅ（図６参照）における３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係を表しており、ドットＤｆは、スライスＳｘの中の領域Ｒｆ（図６参照）における３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係を表している。また、ドットＤ１は、スライスＳ８の中の領域Ｒｍにおける３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係を表しており、ドットＤ２は、スライスＳ３の中の領域Ｒｎにおける３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係を表している。
【００６５】
更に、ドットＤ３は、スライスＳ７の中の一つの領域Ｒｏにおける３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係と、スライスＳ３の中の一つの領域Ｒｐにおける３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係とを表している。異なる領域Ｒであっても、３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係が同じになる場合は、ドットＤ３のように、一つのドットが、異なる領域Ｒにおける３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係を表すことになる。
【００６６】
図８に示すように、ジョイントヒストグラムを作成した後、ステップＳ６に進む。
【００６７】
ステップＳ６では、領域抽出部６７（図１参照）が、スライスＳ１〜Ｓｎが有する全ての領域Ｒの中から、動脈を含む領域Ｒと、静脈を含む領域Ｒを抽出する。動脈を含む領域Ｒは、後述する動脈入力関数を算出するために抽出され、静脈を含む領域Ｒは、後述する静脈出力関数を算出するために抽出される。しかし、スライスＳ１〜Ｓｎには、動脈を含む領域Ｒと、静脈を含む領域Ｒが混在しており、更に、脳組織を含む領域Ｒや、被検体８の頭部８ａの外側に位置する領域Ｒなど、様々な領域Ｒが存在している。そこで、領域抽出部６７は、動脈を含む領域Ｒと、静脈を含む領域Ｒを抽出できるように、図８に示すジョイントヒストグラムを使用している。以下に、ジョイントヒストグラムを使用して、スライスＳ１〜Ｓｎの全領域Ｒの中から、どのようにして、動脈を含む領域Ｒと静脈を含む領域Ｒを抽出しているかについて説明する。
【００６８】
領域抽出部６７は、ジョイントヒストグラムをＣmax−ＴＴＰ面に投影する。
【００６９】
図９は、ジョイントヒストグラムをＣmax−ＴＴＰ面に投影することにより得られた２次元プロファイルである。図９では、プロットＤ１、Ｄ２、およびＤ３の位置が具体的に示されている。
【００７０】
図９のような２次元プロファイルを得た後、領域抽出部６７は、２次元プロファイルの輪郭線を求める。
【００７１】
図１０は、図９に示す２次元プロファイルの輪郭線の一例を示す図である。
【００７２】
輪郭線Ｌｏには、２つのピークαおよびβが現れる。ピークαのピーク濃度到達時間ｔαは、ピークβのピーク濃度到達時間ｔβよりも短くなっている。動脈を含む領域Ｒのピーク濃度到達時間ＴＴＰは、ｔαを中心として、ピーク到達時間ｔ１〜ｔ２の範囲内に現れる可能性が高いことがわかっている（ｔ１〜ｔ２の時間間隔は、例えば、５秒程度）。そこで、領域抽出部６７は、動脈を含む領域Ｒを抽出するため、ジョイントヒストグラムに対して、ピーク到達時間ｔ１〜ｔ２の範囲（以下、「動脈範囲」と呼ぶ）Ａαを決定する。また、静脈を含む領域Ｒのピーク濃度到達時間ＴＴＰは、ｔβを中心として、ピーク到達時間ｔ３〜ｔ４の範囲内に現れる可能性が高いこともわかっている（ｔ３〜ｔ４の時間間隔は、例えば、十秒程度）。そこで、領域抽出部６７は、静脈を含む領域Ｒを抽出するため、ジョイントヒストグラムに対して、ピーク到達時間ｔ３〜ｔ４の範囲（以下、「静脈範囲」と呼ぶ）Ａβを決定する。
【００７３】
図１１は、ジョイントヒストグラムに対して決定された動脈範囲Ａαおよび静脈範囲Ａβを示す図である。
【００７４】
ジョイントヒストグラムの動脈範囲Ａαに含まれている各ドットは、ｔ１〜ｔ２のピーク到達時間ＴＴＰを有する領域Ｒの３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係を表している。一方、ジョイントヒストグラムの静脈範囲Ａβに含まれている各ドットは、ｔ３〜ｔ４のピーク到達時間ＴＴＰを有する領域Ｒの３つの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）の関係を表している。
【００７５】
次に、ジョイントヒストグラムの動脈範囲Ａαおよび静脈範囲Ａβに含まれている各ドットが、スライスＳ１〜Ｓｎの中のどの領域Ｒに対応しているかについて、図１２を参照しながら説明する。
【００７６】
図１２（ａ）は、図１１に示すジョイントヒストグラムの動脈範囲ＡαのＣmax−Ｃｐ面への投影図、図１２（ｂ）は、図１１に示すジョイントヒストグラムの静脈範囲ＡβのＣmax−Ｃｐ面への投影図、図１２（ｃ）は、ジョイントヒストグラムの動脈範囲Ａαおよび静脈範囲Ａβに含まれている各ドットが、スライスＳ１〜Ｓｎの中のどの領域Ｒの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表しているかを説明する図である。
【００７７】
図１２（ｃ）には、説明の便宜上、スライスＳ１〜Ｓｎのうちの４枚のスライスＳ１、Ｓ３、Ｓn-4、Ｓn-1の断面のみが示されている。ジョイントヒストグラムの動脈範囲Ａα（図１２（ａ）参照）に含まれる各ドットは、スライスＳ１、Ｓ３、Ｓn-4、Ｓn-1の中の黒丸で示されている領域Ｒ１又はＲ２の特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表している。領域Ｒ１は、被検体の頭部８ａの内側の領域であり、領域Ｒ２は、被検体の頭部８ａの外側の領域である。一方、ジョイントヒストグラムの静脈範囲Ａβ（図１２（ｂ）参照）に含まれるドットは、スライスＳ１、Ｓ３、Ｓn-4、Ｓn-1の「×」で示されている領域Ｒ３又はＲ４の特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表している。領域Ｒ３は、被検体の頭部８ａの内側の領域であり、領域Ｒ４は、被検体の頭部８ａの外側の領域である。
【００７８】
図１２（ｃ）には、スライスＳ１〜Ｓｎのうちの４枚のスライスＳ１、Ｓ３、Ｓn-4、Ｓn-1についてのみ、領域Ｒ１〜Ｒ４が示されているが、他のスライスについても、スライスＳ１、Ｓ３、Ｓn-4、Ｓn-1と同様に、領域Ｒ１〜Ｒ４を有している。
【００７９】
領域抽出部６７は、スライスＳ１〜Ｓｎの全ての領域Ｒの中から、領域Ｒ１〜Ｒ４を抽出する。
【００８０】
しかし、図１２（ｃ）のスライスＳ１、Ｓ３、Ｓn-4、Ｓn-1を参照すると、被検体８の頭部８ａの外側にも、領域Ｒ２および領域Ｒ４が存在していることがわかる。領域Ｒ２およびＲ４は、被検体８の頭部８ａの外側の領域であるので、明らかに、動脈を含む領域Ｒ１でもなければ、静脈を含む領域Ｒ３でもない。したがって、領域Ｒ２およびＲ４は、動脈入力関数および静脈出力関数の算出には使用することができないノイズである。そこで、ノイズとなる領域Ｒ２およびＲ４を除去するために、ステップＳ７に進む。
【００８１】
ステップＳ７では、ノイズ除去部６８（図１参照）が、ジョイントヒストグラムの動脈範囲Ａαから、領域Ｒ２に対応するドットを除去するとともに、ジョイントヒストグラムの静脈範囲Ａβから、領域Ｒ４に対応するドットを除去する。ノイズ除去部６８は、領域Ｒ２およびＲ４に対応するドットを除去するために、先ず、ジョイントヒストグラムの動脈範囲ＡαおよびＡβの中から、領域Ｒ２およびＲ４の特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表しているドットを検出する。本実施形態では、ノイズ除去部６８は、ジョイントヒストグラムの動脈範囲Ａαの中の一部の範囲Ｎ１〜Ｎ４に含まれているドット（図１２（ａ）参照）を、領域Ｒ２の特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表しているドットとして検出する。更に、ノイズ除去部６８は、ジョイントヒストグラムの静脈範囲Ａβの中の一部の範囲Ｎ５〜Ｎ７に含まれているドット（図１２（ｂ）参照）を、領域Ｒ４の特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表しているドットとして検出する。
【００８２】
上記のように、ジョイントヒストグラムの範囲Ｎ１〜Ｎ７に含まれているドットを検出した後、ノイズ除去部６８は、ジョイントヒストグラムの範囲Ｎ１〜Ｎ７に含まれているドットを除去する。
【００８３】
図１３は、ジョイントヒストグラムの範囲Ｎ１〜Ｎ７からドットを除去した後の様子を説明する図である。
【００８４】
図１３（ａ）は、ドット除去後にのジョイントヒストグラムの動脈範囲Ａαに残ったドットのＣmax−Ｃｐ面への投影図、図１３（ｂ）は、ドット除去後にジョイントヒストグラムの静脈範囲Ａβに残ったドットのＣmax−Ｃｐ面への投影図、図１３（ｃ）は、ジョイントヒストグラムの動脈範囲Ａαおよび静脈範囲Ａβに残ったドットが、スライスＳ１〜Ｓｎの中のどの領域Ｒの特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表しているかを説明する図である。
【００８５】
図１３（ａ）を参照すると、ジョイントヒストグラムの範囲Ｎ１〜Ｎ４に存在していたドットが除去されていることがわかる。ジョイントヒストグラムの範囲Ｎ１〜Ｎ４に存在していたドットが除去されることにより、ステップＳ６において抽出された領域Ｒ（図１２参照）は、被検体８の頭部８ａの内側にしか残っていないことがわかる。
【００８６】
ノイズとなるドットを除去した後、ステップＳ８に進む。
【００８７】
ステップＳ８では、領域選択部６９（図１参照）が、各スライスＳ１〜Ｓｎの被検体８の頭部８ａ内に残った領域Ｒ１およびＲ３の中から、動脈入力関数および静脈出力関数を決定するために使用される領域Ｒ１およびＲ３を絞り込む。動脈入力関数を決定するために使用される領域Ｒ１は、Ｃmaxができるだけ大きいことが望ましく、静脈出力関数を決定するために使用される領域Ｒ３も、Ｃmaxができるだけ大きいことが望ましい。そこで、領域選択部６９は、各スライスＳ１〜Ｓｎの被検体８の頭部８ａ内に残った領域Ｒ１およびＲ３の中から、Ｃmaxの大きい領域を選択する。
【００８８】
図１４は、各スライスＳ１〜Ｓｎの被検体８の頭部８ａ内に残った領域Ｒの中から、どの領域Ｒが選択されるかを説明する図である。
【００８９】
本実施形態では、領域選択部６９は、図１４（ａ）に示すように、ジョイントヒストグラムの動脈範囲Ａαから、ピーク濃度Ｃmaxの大きい順に３個のドットＤ１１、Ｄ１２、およびＤ１３を選択する。ジョイントヒストグラムの動脈範囲Ａαからドットを選択する場合は、ボーラス通過後濃度Ｃｐが正の値であるか負の値であるかとは無関係に、ピーク濃度Ｃmaxの大きい順に３個のドットＤ１１、Ｄ１２、およびＤ１３が選択される。ドットＤ１１、Ｄ１２、およびＤ１３は、それぞれ、スライスＳ３の領域Ｒ１１、スライスＳn-4の領域Ｒ１２、およびスライスＳn-1の領域Ｒ１３の特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表している。図１４（ｃ）には、各領域Ｒ１１、Ｒ１２、およびＲ１３の造影剤濃度−時間曲線Ｃ１１、Ｃ１２、およびＣ１３が概略的に示されている。
【００９０】
上記のように、ジョイントヒストグラムの動脈範囲Ａαからドットを選択する場合は、ボーラス通過後濃度Ｃｐが正の値であるか負の値であるかとは無関係に、ピーク濃度Ｃmaxの大きい順に３個のドットＤ１１、Ｄ１２、およびＤ１３が選択される。しかし、ジョイントヒストグラムの静脈範囲Ａβからは、ボーラス通過後濃度Ｃｐが正の値を有するドットＤｎは選択せずに、ボーラス通過後濃度Ｃｐが負の値を有するドットの中から、Ｃmaxの大きい順に３個のドットＤ２１、Ｄ２２、およびＤ２３を選択する。以下に、この理由について説明する。
【００９１】
図１５は、ジョイントヒストグラムの静脈範囲Ａβからは、ボーラス通過後濃度Ｃｐが正の値を有するドットを選択しない理由を説明する図である。
【００９２】
図１５（ａ）は、スライスＳn-4の頭部８ａ中の一部の部位ＴＡを示す図であり、図１５（ｂ）は、一部の部位ＴＡの拡大図である。
【００９３】
部位ＴＡの中には、動脈ＡＲと静脈ＶＥが示されている。の領域Ｒa1には、動脈血Ｂｆが流入する。通常、動脈血Ｂｆが流入する領域Ｒa1は、静脈血Ｂｖが流入する領域Ｒv1よりも、ピーク濃度Ｃmaxが先に現れる特性がある（図５（ｂ）および（ｃ）参照）。
【００９４】
しかし、部位ＴＡの動脈ＡＲに病巣ＦＯがあり、動脈ＡＲの内径が狭くなっている場合、病巣ＦＯを流れる動脈血Ｂｓは、動脈血Ｂｆよりも血流が遅くなり、この血流の遅い動脈血Ｂｓが、領域Ｒa2内に流入する。したがって、領域Ｒa2におけるピーク濃度Ｃmaxの現れる時刻が、領域Ｒa1におけるピーク濃度Ｃmaxの現れる時刻よりも遅くなり、その結果、領域Ｒa2におけるピーク濃度Ｃmaxの現れる時刻が、静脈を含む領域Ｒv1のＣmaxの現れる時刻に近くなる。したがって、領域Ｒa2のピーク濃度Ｃmaxだけを調べても、領域Ｒa2が、動脈を含む領域であるのか、静脈を含む領域であるのか、区別をすることができない。
【００９５】
しかし、動脈を含む領域Ｒa1およびＲa2では、ボーラス通過後濃度Ｃｐは正の値として算出される傾向があるが、一方、静脈を含む領域Ｒv1では、ボーラス通過後濃度Ｃｐは負の値として算出される傾向があることが一般に知られている（図５（ｂ）および（ｃ）参照）。したがって、図１４（ｂ）を参照すると、ジョイントヒストグラムの静脈範囲Ａβに、ボーラス通過後濃度Ｃｐが正の値を有するドットＤｎが含まれているが、このドットＤｎは、動脈を含む画素Ｒa2（図１５参照）の可能性がある。そこで、領域選択部６９は、ジョイントヒストグラムの静脈範囲Ａβにおいて、ボーラス通過後濃度Ｃｐが負の値を有するドットの中から、ピーク濃度Ｃmaxの大きい順に３個のドットＤ２１、Ｄ２２、およびＤ２３を選択する。ドットＤ２１、Ｄ２２、およびＤ２３は、それぞれ、スライスＳ１の領域Ｒ２１、スライスＳ３の領域Ｒ２２、およびスライスＳn-4の領域Ｒ２３の特徴量（Ｃmax、ＴＴＰ、Ｃｐ）を表している。図１４（ｃ）には、各領域Ｒ２１、Ｒ２２、およびＲ２３の造影剤濃度−時間曲線Ｃ２１、Ｃ２２、およびＣ２３が概略的に示されている。
【００９６】
このようにして、各スライスＳ１〜Ｓｎの被検体８の頭部８ａ内に残った領域Ｒ１およびＲ３の中から、領域Ｒ１１、Ｒ１２、Ｒ１３、Ｒ２１、Ｒ２２、およびＲ２３が選択される。
【００９７】
領域Ｒ１１、Ｒ１２、Ｒ１３、Ｒ２１、Ｒ２２、およびＲ２３を選択した後、ステップＳ９に進む。
【００９８】
ステップＳ９では、関数算出手段７０（図１参照）が、選択された３つの領域Ｒ１１、Ｒ１２、およびＲ１３における造影剤濃度−時間曲線Ｃ１１、Ｃ１２、およびＣ１３を加算平均することにより、動脈入力関数を算出する。更に、関数算出手段７０は、選択された３つの領域Ｒ２１、Ｒ２２、およびＲ２３における造影剤濃度−時間曲線Ｃ２１、Ｃ２２、およびＣ２３を加算平均することにより、静脈出力関数を算出する。
【００９９】
図１６は、算出された動脈入力関数と、静脈出力関数を示す図である。
【０１００】
動脈入力関数Ｆａおよび静脈出力関数Ｆｖを算出したら、プログラムを終了する。
【０１０１】
本実施形態では、ピーク濃度Ｃｐが負の値を有するドットの中から、Ｃmaxの大きい順に３個のドットを選択することによって、スライスＳ１〜Ｓｎの中から、静脈出力関数Ｆｖを算出するために使用される領域Ｒ２１、Ｒ２２、およびＲ２３を選択している。したがって、領域Ｒ２１、Ｒ２２、およびＲ２３の中に、動脈を含む領域Ｒa2（図１５参照）が含まれないようにすることができ、より正確な静脈出力関数Ｆｖを算出することができる。
【０１０２】
尚、本実施形態では、３つの領域Ｒ２１、Ｒ２２、およびＲ２３における造影剤濃度−時間曲線Ｃ２１、Ｃ２２、およびＣ２３（図１４（ｃ）参照）を加算平均することにより、動脈入力関数Ｆａおよび静脈出力関数Ｆｖを算出しているが、１つの領域、２つの領域、又は４つ以上の領域における造影剤濃度−曲線を用いて、動脈入力関数Ｆａおよび静脈出力関数Ｆｖを算出してもよい。また、本実施形態では、造影剤濃度−曲線を加算平均することにより動脈入力関数Ｆａおよび静脈出力関数Ｆｖを算出しているが、造影剤濃度−曲線を加算平均する以外の方法で動脈入力関数Ｆａおよび静脈出力関数Ｆｖを算出してもよい。
【０１０３】
本実施形態では、造影剤濃度−時間曲線ＧａおよびＧｖ（図５（ｂ）および（ｃ）参照）の特徴量として、３つの特徴量（ピーク濃度Ｃmax、ピーク濃度到達時間ＴＴＰ、およびボーラス通過後濃度Ｃp）を算出している。しかし、算出される特徴量は、動脈を含む領域と、静脈を含む領域とを区別するのに有用な特徴量であれば、他の特徴量を使用することもできる。
【０１０４】
本実施形態では、ノイズ除去部６８は、領域抽出部６７でスライスＳ１〜Ｓｎから領域が抽出された後に、ノイズとなる領域を除去している。しかし、領域抽出部６７でスライスＳ１〜Ｓｎから領域を抽出する前に、除去してもよい。
【０１０５】
また、本実施形態では、ＭＲＩ装置１（図１参照）が、被検体８から収集したフレーム画像のデータに基づいて、動脈入力関数Ｆａおよび静脈出力関数Ｆｖを算出している。しかし、動脈入力関数Ｆａおよび静脈出力関数Ｆｖの算出は、ＭＲＩ装置１ではなく、動脈入力関数Ｆａおよび静脈出力関数Ｆｖを算出する別の装置で行ってもよい。
【０１０６】
更に、本実施形態では、被検体８の頭部８ａの血流動態を解析しているが、本発明は、頭部８ａ以外の他の部位の血流動態を解析するのにも適用可能である。
【図面の簡単な説明】
【０１０７】
【図１】本発明の一実施形態の磁気共鳴イメージング装置１の概略図である。
【図２】磁気共鳴イメージング装置１の処理フローを示す図である。
【図３】被検体８に設定されたスライスの一例である。
【図４】スライスＳ１〜ＳｎのうちのスライスＳｘから得られるフレーム画像を示す概念図である。
【図５】被検体８の頭部８ａに設定されたスライスＳxの断面領域において、造影剤の濃度の時間変化を算出する方法を説明する図である。
【図６】スライスＳｘについてのピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップを概略的に示す図である。
【図７】各スライスＳ１〜Ｓｎの３つのマップ（ピーク濃度Ｃmaxのマップ、ピーク濃度到達時間ＴＴＰのマップ、およびボーラス通過後濃度Ｃｐのマップ）を概略的に示す図である。
【図８】ジョイントヒストグラムの一例を概略的に示す図である。
【図９】ジョイントヒストグラムをＣmax−ＴＴＰ面に投影することにより得られた２次元プロファイルである。
【図１０】図９に示す２次元プロファイルの輪郭線の一例を示す図である。
【図１１】ジョイントヒストグラムに対して決定された動脈範囲Ａαおよび静脈範囲Ａβを示す図である。
【図１２】Ｃmax−Ｃｐ面への投影図を示す図である。
【図１３】ジョイントヒストグラムの範囲Ｎ１〜Ｎ７からドットを除去した後の様子を説明する図である。
【図１４】各スライスＳ１〜Ｓｎの被検体８の頭部８ａ内に残った領域Ｒの中から、どの領域Ｒが選択されるかを説明する図である。
【図１５】ジョイントヒストグラムの静脈範囲Ａβからは、ボーラス通過後濃度Ｃｐが正の値を有するドットを選択しない理由を説明する図である。
【図１６】算出された動脈入力関数と、静脈出力関数を示す図である。
【符号の説明】
【０１０８】
１  ＭＲＩ装置
２  ガントリ
３  クレードル
４  受信コイル
５  造影剤注入装置
６  制御装置
７  入力装置
８  被検体
８ａ  頭部
９  オペレータ
２１  ボア
２２  超伝導コイル
２３  勾配コイル
２４  送信コイル
６１  コイル制御手段
６２  造影剤濃度プロファイル作成手段
６３  特徴量算出手段
６４  マップ作成手段
６５  ジョイントヒストグラム作成手段
６６  領域決定手段
６７  領域抽出部
６８  ノイズ除去部
６９  領域選択部
７０  関数算出手段
【図１】
2010104671.tif 000003
【図２】
2010104671.tif 000004
【図３】
2010104671.tif 000005
【図４】
2010104671.tif 000006
【図５】
2010104671.tif 000007
【図６】
2010104671.tif 000008
【図７】
2010104671.tif 000009
【図８】
2010104671.tif 000010
【図９】
2010104671.tif 000011
【図１０】
2010104671.tif 000012
【図１１】
2010104671.tif 000013
【図１２】
2010104671.tif 000014
【図１３】
2010104671.tif 000015
【図１４】
2010104671.tif 000016
【図１５】
2010104671.tif 000017
【図１６】
2010104671.tif 000018