(57)【特許請求の範囲】
【請求項１】
血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部と、
前記血管領域内における血流速度および血流方向を表す流速ベクトルを取得する流速ベクトル取得部と、
前記血管領域の延伸方向に延びる中心軸を設定し、前記流速ベクトルの前記中心軸の方向成分を時間軸方向に積分した積分値を算出し、前記積分値の変化に基づいて、前記血管領域内の血流の逆流の有無を特定する逆流特定部とを備えた血流解析装置。
【請求項２】
前記逆流特定部が、前記積分値が単調増加または単調減少でない場合に、前記逆流が有ると特定する請求項１記載の血流解析装置。
【請求項３】
前記逆流の箇所を前記血管領域の座標軸上にマッピングしたマップ画像を表示部に表示させる表示制御部を備えた請求項１または２記載の血流解析装置。
【請求項４】
前記表示制御部が、前記流速ベクトルを予め設定された期間だけ積分した前記積分値を用いて前記マップ画像を生成する請求項３記載の血流解析装置。
【請求項５】
前記予め設定された期間が、１心拍の期間である請求項４記載の血流解析装置。
【請求項６】
前記表示制御部が、前記流速ベクトルを時間軸方向に積分する毎に前記マップ画像を時系列に生成して動画として前記表示部に表示させる請求項３記載の血流解析装置。
【請求項７】
前記表示制御部が、前記マップ画像を２次元または３次元の血管画像に重畳表示させる請求項３から６いずれか１項記載の血流解析装置。
【請求項８】
前記表示制御部が、前記積分値に基づいて、前記逆流の量の情報を前記表示部にさせる請求項３から７いずれか１項記載の血流解析装置。
【請求項９】
前記表示制御部が、前記マップ画像以外の血流を示す指標を前記表示部に表示させる請求項３から８いずれか１項記載の血流解析装置。
【請求項１０】
前記血流を示す指標が、前記流速ベクトル、流線または流跡線である請求項９記載の血流解析装置。
【請求項１１】
前記逆流特定部が、前記血管領域のボクセル毎または複数のボクセル毎に前記積分値を算出する請求項１から１０いずれか１項記載の血流解析装置。
【請求項１２】
前記逆流特定部が、前記血管領域の中心軸方向に直交する断面毎に前記積分値を算出する請求項１から１０いずれか１項記載の血流解析装置。
【請求項１３】
前記流速ベクトル取得部が、時系列に撮影された３次元の血管画像または任意の時点において撮影された３次元の血管画像に基づいて、前記流速ベクトルを取得する請求項１から１２いずれか１項記載の血流解析装置。
【請求項１４】
前記流速ベクトル取得部が、前記血管領域の数値流体力学(Computational  Fluid  Dynamics)モデル、３次元シネ位相コントラスト磁気共鳴法または３次元の超音波画像を用いて前記流速ベクトルを取得する請求項１から１２いずれか１項記載の血流解析装置。
【請求項１５】
血管を含む被写体を撮影した血管画像から血管領域を抽出し、
前記血管領域内における血流速度および血流方向を表す流速ベクトルを取得し、
前記血管領域の延伸方向に延びる中心軸を設定し、前記流速ベクトルの前記中心軸の方向成分を時間軸方向に積分した積分値を算出し、前記積分値の変化に基づいて、前記血管領域内の血流の逆流の有無を特定する血流解析方法。
【請求項１６】
コンピュータを、
血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部と、
前記血管領域内における血流速度および血流方向を表す流速ベクトルを取得する血流ベクトル取得部と、
前記血管領域の延伸方向に延びる中心軸を設定し、前記流速ベクトルの前記中心軸の方向成分を時間軸方向に積分した積分値を算出し、前記積分値の変化に基づいて、前記血管領域内の血流の逆流の有無を特定する逆流特定部として機能させる血流解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用画像から抽出された血管領域を用いて、その血管領域内を流れる血流の解析を行う血流解析装置および方法並びにプログラムに関するものである。
【背景技術】
【０００２】
近年、心臓などを撮影した医用画像を用いて、心臓などの血管内の血流を解析することが行われている。【０００３】
このような医用画像を用いた血流解析方法としては、たとえば３次元シネ位相コントラスト磁気共鳴法によって撮影されたＭＲＩ（Magnetic Resonance Imaging）画像を用いて、ボクセル毎、ピクセル毎または領域毎に流速ベクトルを算出し、これを表示する方法が提案されている。【０００４】
また、特開２０１４−１８８３２３号公報においては、医用画像から血管モデルを生成し、その血管モデルに基づいて、ＣＦＤ（Computational Fluid Dynamics）を行って流速ベクトルを求めることが提案されている。【０００５】
ここで、血流解析においては、血流が逆流している箇所を特定することが、診断上重要である。上述したように流速ベクトルを表示した場合、ユーザが、その向きと大きさを観察することによって逆流箇所を特定することができる。【０００６】
また、医用画像から抽出された血管領域に関心領域を設定し、その関心領域に流れる血流を表す流線を表示させて観察することにより、逆流を確認することもできる。【０００７】
また、超音波診断の分野でも逆流を確認することが行われている。たとえば特開２００９−１５３９１９号公報および特開平８−１５０１４２号公報においては、超音波診断装置を用いてカラードップラ計測を行って、逆流を可視化することが提案されている。
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、上述した３次元シネ位相コントラスト磁気共鳴法およびＣＦＤを用いて逆流箇所を確認する場合、ユーザが関心のある箇所付近の流速ベクトルの向きを確認したり、関心領域を設定した際の流線の形成具合を確認したりすることによって、逆流箇所を確認することができるが、この場合、ユーザが目的を持って観察しようとした場所以外で逆流があった場合、逆流を見落とす可能性がある。【０００９】
また、超音波診断によって逆流を確認する場合、ユーザが関心のある領域にプローブを当ててカラードップラ計測が行われるが、超音波画像は２次元的であり、２次元画像上での流れの正負を表示させることはできるが、血管形状に沿った流れを抽出することは困難であり、逆流として画像化することは難しい。【００１０】
本発明は、上記事情に鑑み、想定外の箇所に逆流があった場合でも、その逆流箇所をとらえることができ、逆流箇所の見落としリスクを減らすことができる血流解析装置および方法並びにプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【００１１】
本発明の血流解析装置は、血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部と、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する流速ベクトル取得部と、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する逆流特定部とを備える。【００１２】
また、上記本発明の血流解析装置において、逆流特定部は、積分値が単調増加または単調減少でない場合に、逆流が有ると特定することができる。【００１３】
また、上記本発明の血流解析装置においては、逆流の箇所を血管領域の座標軸上にマッピングしたマップ画像を表示部に表示させる表示制御部を備えることができる。【００１４】
また、上記本発明の血流解析装置において、表示制御部は、流速ベクトルを予め設定された期間だけ積分した積分値を用いてマップ画像を生成することができる。【００１５】
また、上記本発明の血流解析装置において、上記予め設定された期間は、１心拍の期間とすることができる。【００１６】
また、上記本発明の血流解析装置において、表示制御部は、流速ベクトルを時間軸方向に積分する毎にマップ画像を時系列に生成して動画として表示部に表示させることができる。【００１７】
また、上記本発明の血流解析装置において、表示制御部は、マップ画像を２次元または３次元の血管画像上に重畳表示させることができる。【００１８】
また、上記本発明の血流解析装置において、表示制御部は、積分値に基づいて、逆流の量の情報を表示部にさせることができる。【００１９】
また、上記本発明の血流解析装置において、表示制御部は、マップ画像以外の血流を示す指標を表示部に表示させることができる。【００２０】
また、上記本発明の血流解析装置においては、血流を示す指標として、流速ベクトル、流線または流跡線を表示させることができる。【００２１】
また、上記本発明の血流解析装置において、逆流特定部は、血管領域のボクセル毎または複数のボクセル毎に積分値を算出することができる。【００２２】
また、上記本発明の血流解析装置において、逆流特定部は、血管領域の中心軸方向に直交する断面毎に積分値を算出することができる。【００２３】
また、上記本発明の血流解析装置において、流速ベクトル取得部は、時系列に撮影された３次元の血管画像または任意の時点において撮影された３次元の血管画像に基づいて、流速ベクトルを取得することができる。【００２４】
また、上記本発明の血流解析装置において、流速ベクトル取得部は、血管領域の数値流体力学(Computational  Fluid  Dynamics)モデル、３次元シネ位相コントラスト磁気共鳴法または３次元の超音波画像を用いて流速ベクトルを取得することができる。【００２５】
本発明の血流解析方法は、血管を含む被写体を撮影した血管画像から血管領域を抽出し、血管領域内における血流速度および血流方向を表す流速ベクトルを取得し、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する。【００２６】
本発明の血流解析プログラムは、コンピュータを、血管を含む被写体を撮影した血管画像から血管領域を抽出する血管領域抽出部と、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する血流ベクトル取得部と、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する逆流特定部として機能させる。【００２７】
本発明の他の血流解析装置は、コンピュータに実行させるための命令を記憶するメモリと、記憶された命令を実行するよう構成されたプロセッサを備え、プロセッサが、血管を含む被写体を撮影した血管画像から血管領域を抽出する処理と、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する処理と、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する処理とを実行するように構成されている。
【発明の効果】
【００２８】
本発明の血流解析装置および方法並びにプログラムによれば、血管を含む被写体を撮影した血管画像から血管領域を抽出し、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する。そして、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する。このように流速ベクトルの中心軸方向成分の積分値に基づいて逆流の有無を特定することによって、想定外の箇所に逆流があった場合でも、その逆流箇所をとらえることができ、逆流箇所の見落としリスクを減らすことができる。
【図面の簡単な説明】
【００２９】
【図１】本発明の血流解析装置の一実施形態を用いた血流解析システムの概略構成を示すブロック図
【図２】血管領域内における流速ベクトルを模式的に示す図
【図３】逆流箇所の特定方法を説明するための図
【図４】積分値が単調増加または単調減少しない例を示す図
【図５】断面毎に積分値を算出した場合の逆流箇所の表示例を示す図
【図６】本発明の血流解析装置の一実施形態を用いた血流解析システムの作用を説明するためのフローチャート
【図７】複数のボクセル毎に積分値を算出した場合の逆流箇所の表示例を示す図
【発明を実施するための形態】
【００３０】
以下、図面を参照して本発明の血流解析装置の一実施形態を用いた血流解析システムについて説明する。図１は、本実施形態の血流解析システムの概略構成を示すブロック図である。
【００３１】
本実施形態の血流解析システムは、具体的には、図１に示すように、血流解析装置１と、医用画像保管サーバ２と、表示部３と、入力部４とを備えている。
【００３２】
血流解析装置１は、コンピュータに本実施形態の血流解析プログラムをインストールしたものである。 
【００３３】
血流解析装置１は、中央処理装置（ＣＰＵ（central processing unit））、半導体メモリ、およびハードディスクやＳＳＤ(Solid State Drive)等のストレージデバイスを備えている。ストレージデバイスには、本実施形態の血流解析プログラムがインストールされており、この血流解析プログラムが中央処理装置によって実行されることによって、図１に示す画像取得部１０、血管領域抽出部１１、流速ベクトル取得部１２、逆流特定部１３および表示制御部１４が動作する。
【００３４】
血流解析プログラムは、ＤＶＤ（Digital Versatile Disc）およびＣＤ−ＲＯＭ（Compact Disc Read Only Memory）などの記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。または、血流解析プログラムは、ネットワークに接続されたサーバコンピュータの記憶装置もしくはネットワークストレージに対して、外部からアクセス可能な状態で記憶される。そして、外部からの要求に応じてコンピュータにダウンロードされ、インストールされる。
【００３５】
画像取得部１０は、予め撮影された患者の医用画像を取得するものである。３次元医用画像６は、本実施形態においては、たとえばＣＴ（Computed Tomography）装置、ＭＲＩ
（Magnetic Resonance Imaging）装置および超音波撮影装置などによって撮影された３次元医用画像６を取得する。本実施形態では、患者の心臓（本発明の被写体に相当する）の３次元医用画像６（本発明の血管画像に相当する）を取得する場合について説明するが、これに限らず、肺、肝臓および頭部などその他の臓器でもよい。
【００３６】
３次元医用画像６は、医用画像保管サーバ２に患者の識別情報とともに予め保管されており、画像取得部１０は、入力部４などを用いてユーザによって入力された患者の識別情報に基づいて、その識別情報を有する３次元医用画像６を医用画像保管サーバ２から読み出して一時記憶するものである。
【００３７】
血管領域抽出部１１は、３次元医用画像６から血管領域を抽出するものである。具体的には、本実施形態の血管領域抽出部１１は、心臓の３次元医用画像６から冠動脈の領域を血管領域として抽出する。血管領域抽出部１１は、たとえば心臓の３次元医用画像６に対して多重解像度変換を行い、各解像度の画像に対してヘッセ行列の固有値解析を行い、各解像度の画像における解析結果を統合することによって、心臓領域中の様々なサイズの線構造（血管）の集合体として、冠動脈の領域を抽出する（たとえばY Sato, et al.、「Three-dimensional multi-scale line filter for segmentation and visualization of curvilinear structures in medical images.」、Medical Image Analysis、1998年6月、Vol.2、No.2、p.p.143-168など参照）。また、さらに最小全域木アルゴリズム等を用いて、抽出された各線構造の中心点を連結することにより、冠動脈を表す木構造のデータを生成し、抽出された冠動脈の中心点を結ぶ芯線上の各点（木構造データの各ノード）において、芯線に直交する断面を求め、各断面において、グラフカット法などの公知のセグメンテーション手法を用いて冠動脈の輪郭を認識し、その輪郭を表す情報を木構造データの各ノードに関連づけることによって、冠動脈の領域を抽出するようにしてもよい。
【００３８】
なお、冠動脈の領域の抽出方法としては上記の方法に限らず、領域拡張法などのその他の公知な手法を用いるようにしてもよい。
【００３９】
流速ベクトル取得部１２は、血管領域内における血流速度および血流方向を表す流速ベクトルを取得するものである。流速ベクトルの取得方法としては、種々の方法を採用することができる。具体的には、たとえば血管領域抽出部１１によって抽出された血管領域を用いてＣＦＤ（Computational Fluid Dynamics）による血流解析を行うことによって流速ベクトルを取得することができる。なお、この際、任意の時点において撮影された３次元医用画像６から抽出された血管領域を用いて血流解析を行うようにしてもよいし、時系列に撮影された３次元医用画像６のそれぞれから抽出された血管領域を用いて血流解析を行うようにしてもよい。
【００４０】
また、画像取得部１０において３次元シネ位相コントラスト磁気共鳴法（3D cine PC MRI）によって撮影されたＭＲＩ画像を取得し、そのＭＲＩ画像に基づいて取得された血管領域内の速度情報を用いて流速ベクトルを取得するようにしてもよい。
【００４１】
また、画像取得部１０においてドップラー計測によって時系列に撮影された３次元の超音波画像を取得し、その超音波画像に基づいて取得された血管領域内の速度情報を用いて流速ベクトルを取得するようにしてもよい。図２は、血管領域Ｂ内における流速ベクトルＶを模式的に示した図である。
【００４２】
逆流特定部１３は、流速ベクトル取得部１２によって取得された血管領域内の流速ベクトルに基づいて、血管領域内の血流の逆流の有無を特定するものである。
【００４３】
逆流特定部１３は、具体的には、図２に示すように、血管領域Ｂの延伸方向に延びる中心軸Ｃを設定し、流速ベクトルＶの中心軸Ｃの方向成分を時間軸方向に積分した積分値を算出し、その積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する。
【００４４】
流速ベクトルＶの中心軸Ｃの方向成分については、たとえば図２および図３に示すように、中心軸Ｃ上の各点について、中心軸Ｃに直交する断面Ｓを設定し、その断面Ｓ内の流速ベクトルＶの中心軸方向成分Ｖｃを算出する。そして、中心軸方向成分Ｖｃを時間軸方向に積分した積分値を算出し、その積分値が、単調増加するかまたは単調減少するかを判定する。すなわち、たとえば血管領域内において逆流が生じていない場合には、上述した積分値は単調増加するか、または単調減少する。これに対し、血管領域内において逆流が生じている場合には、流速ベクトルＶの中心軸方向成分Ｖｃが逆方向に変化するので、中心軸方向成分Ｖｃの積分値は、図４ＩまたはＩＩに示すように、単調増加または単調減少とはならない。
【００４５】
したがって、逆流特定部１３は、中心軸方向成分Ｖｃの積分値が単調増加するかまたは単調減少する場合には、逆流は無いと判定し、中心軸方向成分Ｖｃの積分値が単調増加しないまたは単調減少しない場合には、逆流が有ると判定する。なお、中心軸方向成分Ｖｃの積分値としては、断面Ｓ内の各流速ベクトルＶの中心軸方向成分Ｖｃの積分値の平均値を用いるようにしてもよいし、中央値を用いるようにしてもよい。
【００４６】
なお、血管領域Ｂの中心軸Ｃに直交する断面Ｓの中心軸方向の設定間隔については、ユーザによって任意の間隔に設定可能である。また、積分値を算出する際の積分期間についても、ユーザによって任意の間隔に設定可能である。本実施形態のように、心臓の血管内を流れる血流の逆流の有無を特定する場合には、積分期間を１心拍の期間とすることが好ましい。
【００４７】
図１に戻り、表示制御部１４は、逆流特定部１３によって特定された逆流箇所の情報を表示部３に表示させるものである。具体的には、表示制御部１４は、逆流箇所を血管領域の座標軸上にマッピングしたマップ画像を生成し、表示部３に表示させる。図５は、血管領域Ｂの座標軸上に逆流箇所Ｒをマッピングしたマッピング画像の一例を示す図である。本実施形態においては、上述したように血管領域の中心軸に直交する断面を設定し、断面毎に積分値を算出するので、図５に示すように、逆流箇所Ｒは断面単位で表示される。
【００４８】
また、断面毎に算出される積分値は逆流量を表す。したがって、積分値に応じて表示を変更することによって、逆流量を表す画像とすることができる。この際、全積分期間の積分値に基づいて逆流量を表す画像を生成して表示することによって、静止画として表示するようにしてもよいし、流速ベクトルを時間軸方向に積分する毎に、すなわち積分期間内における各積分過程の積分値毎に、逆流量を表す画像を順次生成して表示することによって、動画として表示するようにしてもよい。
【００４９】
このように逆流量を表す画像を表示させることによって、直観的に逆流量をとらえることができ、診断や治療の際に役立てることができる。
【００５０】
図５に示すマップ画像は、積分値が大きいほど濃度を大きくしたマップ画像の一例である。マップ画像としては、カラー表示することが望ましく、その場合には、積分値に応じて色相を変化させたり、彩度を変化させたりすることが望ましい。
【００５１】
また、マップ画像を３次元医用画像６上に重畳表示させるようにしてもよい。なお、マップ画像を重ね合わせる画像としては、３次元医用画像６に限らず、２次元の断面画像上にマップ画像を重畳表示するようにしてもよい。
【００５２】
図１に戻り、医用画像保管サーバ２は、各種データを保存して管理するコンピュータであり、大容量記憶装置およびデータベース管理用プログラムを備えている。医用画像保管サーバ２は、予め撮影された３次元医用画像６などをネットワーク経由でＣＴ装置またはＭＲＩ装置などの撮影装置から取得し、大容量記憶装置に保存して管理する。
【００５３】
入力部４は、ユーザによる種々の設定入力を受け付けるものであり、キーボードやマウスなどの入力デバイスを備えたものである。入力部４は、たとえば患者の識別情報の設定入力などを受け付けるものである。
【００５４】
次に、本実施形態の血流解析システムの作用について、図６に示すフローチャートを参照しながら説明する。
【００５５】
まず、ユーザによる患者の識別情報などの設定入力に応じて、その患者の医用画像が画像取得部１０によって取得される（Ｓ１０）。本実施形態においては、上述したように心臓の３次元医用画像６が取得される。
【００５６】
そして、画像取得部１０によって取得された心臓の３次元医用画像６は、血管領域抽出部１１に入力され、血管領域抽出部１１は、入力された３次元医用画像６から血管領域を抽出する（Ｓ１２）。
【００５７】
そして、流速ベクトル取得部１２において、たとえば血管領域を用いてＣＦＤにより血流解析が行われ、血管領域内の流速ベクトルが取得される（Ｓ１４）。そして、逆流特定部１３において、血管領域内の流速ベクトルに基づいて、血流の逆流の有無が特定される。具体的には、逆流特定部１３は、上述したように血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、その積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する（Ｓ１６）。
【００５８】
そして、逆流特定部１３によって特定された逆流箇所の情報が表示制御部１４に出力され、表示制御部１４は、入力された情報に基づいてマップ画像を生成して表示部３に表示させる（Ｓ１８）。
【００５９】
上記実施形態の血流解析システムによれば、血管を含む被写体を撮影した血管画像から血管領域を抽出し、血管領域内における血流速度および血流方向を表す流速ベクトルを取得する。そして、血管領域の延伸方向に延びる中心軸を設定し、流速ベクトルの中心軸の方向成分を時間軸方向に積分した積分値を算出し、積分値の変化に基づいて、血管領域内の血流の逆流の有無を特定する。このように流速ベクトルの中心軸方向成分の積分値に基づいて逆流の有無を特定することによって、想定外の箇所に逆流があった場合でも、その逆流箇所をとらえることができ、逆流箇所の見落としリスクを減らすことができる。
【００６０】
また、上記実施形態の血流解析システムにおいては、血管領域の中心軸方向に直交する断面毎に流速ベクトルの積分値を算出するようにしたが、これに限らず、血管領域のボクセル毎または複数のボクセル毎に積分値を算出するようにしてもよい。図７は、複数のボクセル毎に積分値を算出し、逆流箇所Ｒ１，Ｒ２を特定してマップ画像を生成した例を示す図である。複数のボクセル毎の積分値としては、たとえば複数のボクセルのそれぞれについて算出された積分値の平均値を用いるようにしてもよいし、中央値などを用いるようにしてもよい。このようにボクセル毎または複数のボクセル毎に積分値を算出することによって、逆流箇所をより詳細に特定することができる。
【００６１】
また、上記実施形態の血流解析システムにおいては、表示制御部１４が、マップ画像以外の血流を示す指標を表示部３に表示させるようにしてもよい。具体的には、血流を示す指標として、流速ベクトルを表す矢印をマップ画像に重畳表示したり、血流を示す流線または流跡線などをマップ画像に重畳表示したりしてもよい。
【００６２】
また、上記実施形態の血流解析システムにおいては、３次元医用画像６を用いて逆流箇所を特定するようにしたが、必ずしも３次元的な画像でなくてもよく、２次元の断面画像を用いて、上記実施形態と同様にして流速ベクトルを取得し、その血管領域の中心軸方向成分の積分値の変化に基づいて逆流箇所を特定するようにしてもよい。
【符号の説明】
【００６３】
１    血流解析装置
２    医用画像保管サーバ
３    表示部
４    入力部
６    ３次元医用画像
１０  画像取得部
１１  血管領域抽出部
１２  流速ベクトル取得部
１３  逆流特定部
１４  表示制御部
Ｂ    血管領域
Ｃ    中心軸
Ｒ    逆流箇所
Ｒ１，Ｒ２  逆流箇所
Ｓ    断面
Ｖ    流速ベクトル
Ｖｃ  流速ベクトルの中心軸方向成分
【図１】
0006690053.tif 000002
【図２】
0006690053.tif 000003
【図３】
0006690053.tif 000004
【図４】
0006690053.tif 000005
【図５】
0006690053.tif 000006
【図６】
0006690053.tif 000007
【図７】
0006690053.tif 000008