JP2013111358A - 画像処理装置、方法およびプログラム - Google Patents

Abstract

【課題】ループ状の部分を含む血管において、分枝元の枝から分枝した部分血管における血流の方向を特定する。
【解決手段】 画像処理装置１に、グラフ構造を基準点Ｓから下流側に辿っていく経路上に位置する分枝ノードＢにおいて、分枝元の枝Ｅ１上に位置する複数のノードＰ１１〜Ｐ１５の位置に基づく方向を有し、かつ、分枝ノードＢに向かう側の向きを有するように規定された第１のベクトルと、枝Ｅ１から分枝した枝Ｅ３上に位置する複数のノードＰ３１〜Ｐ３５の位置に基づく方向を有し、かつ、分枝ノードＢから離れる側の向きを有するように規定された第２のベクトルとがなす角θを用いて、枝Ｅ３における血流の方向が分枝ノードＢから離れる方向である尤度を算出し、算出された尤度を用いてその枝Ｅ３における血流の方向が分枝ノードＢから離れる方向であるか分枝ノードＢに向かう方向であるかを判定する判定手段３０を備える。
【選択図】図２

Description

本発明は、ループ状の部分を含む血管において、分枝元の枝から分枝した部分血管における血流の方向を特定する画像処理装置、方法およびプログラムに関する。

肝臓・肺などの人体の様々な臓器において手術を行う際には、手術をよりスムーズに進められるようにするため、術前に手術をどのように行うかを綿密に計画することが重要であり、手術の計画を立てる際には対象の臓器における血管の走行および各血管における血流の方向を正確に把握しておくことが必要とされている。

特許文献１には、特定の臓器において腫瘍を含む患部を切除する手術を行う際にその切除すべき領域を適切に決定できるようにするため、血管における血流の方向はその血管を表す木構造を起始部（根）から離れる側に辿っていく方向であるとの認識のもと、画像データから対象の血管が撮影された木構造の領域を抽出し、抽出された領域を画像表示装置上に表示させ、さらに、表示された領域上の任意の位置を指定する操作者の入力に応じて、その指定された位置から末梢部までの支配領域を特定し、他の領域と識別可能に表示することが提案されている。

特許４６８８３６１号公報

しかし、対象の血管が心臓、肝臓、肺、腎臓、胃等の木構造を有する血管であれば、その構造を起始部から離れる側に辿っていく経路は常に一意であることから、その辿っていく方向をその血管における血流の方向とする上記特許文献１に記載の方法により血流の方向を特定することが可能であるが、対象の血管がたとえば図７に示すようなループ状の部分を含む血管である場合、その構造を起始部から離れる側に辿っていく経路は複数通り存在し、単純にその辿っていく方向を血流の方向として特定する上記特許文献１に記載の方法では、血流の方向を正しく特定することができないという問題がある。なお、図７は大腸の血管を示すものであって、図中点線で囲んでいる部分Ａ等にループ状の部分が存在する。

本発明は、上記事情に鑑み、ループ状の部分を含む血管において、分枝元の枝から分枝した部分血管における血流の方向を特定することができる画像処理装置、方法およびプログラムを提供することを目的とするものである。

本発明の画像処理装置は、入力装置による操作者の入力に応じて、予め取得しておいたループ状の部分を含む血管の構造を表すグラフ構造上に基準点を設定する基準点設定手段と、グラフ構造の前記設定された基準点が位置する枝において、基準点より血流の下流側を決定する下流側決定手段と、グラフ構造の前記基準点から前記決定された血流の下流側に辿っていく経路上に位置する任意の分枝ノードにおいて、分枝元の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、分枝ノードに向かう側の向きを有するように規定された第１のベクトルと、分枝元の枝から分枝したいずれか一本の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、分枝ノードから離れる側の向きを有するように規定された第２のベクトルとがなす角θ（0＜θ＜180°）を用いて、その一本の枝における血流の方向が分枝ノードから離れる方向である尤度を角θが小さいほど大きい値となるように算出し、算出された尤度を用いてその一本の枝における血流の方向が分枝ノードから離れる方向であるか分枝ノードに向かう方向であるかを判定する判定手段とを備えたことを特徴とする。

ここで、グラフ構造とは、ノード群とノード間の連結関係を表すエッジ群で構成されるものをいう。

上記本発明の画像処理装置において、判定手段は、前記尤度を、分枝元の枝上に位置する複数のノードに対応する血管上の複数の位置それぞれにおける血管の径に基づいて規定された第１の代表値が、前記一本の枝上に位置する複数のノードに対応する血管上の複数の位置それぞれにおける血管の径に基づいて規定された第２の代表値より大きいほど大きい値となるように算出するものであってもよい。ここで、血管上の複数の位置それぞれにおける血管の径に基づいて規定された代表値（第１の代表値、第２の代表値）というのは、血管上の複数の位置それぞれにおける血管の径の全体的な大きさが大きくなればなるほど大きくなるような値を意味する。その典型的な例として平均値、最大値、最小値等がある。

この場合、尤度Lの算出は、たとえば下記式（１）により行うことができる。下記式（１）において、Ｖ１は第１のベクトルを表し、Ｖ２は第２のベクトルを表し、F（Ｖ１，Ｖ２）は第１のベクトルＶ１と第２のベクトルＶ２がなす角θが小さいほど大きい値となる関数を表す。Ｒ１は第１の代表値を表し、Ｒ２は第２の代表値を表し、Ｇ（Ｒ１，Ｒ２）は第１の代表値Ｒ１から第２の代表値Ｒ２を減算した差（あるいは、第２の代表値Ｒ２に対する第１の代表値Ｒ１の比）が大きいほど大きい値となる関数を表す。また、Ｗ１，Ｗ２は上記各関数を重み付けする係数である。

本発明の画像処理装置において、判定手段は、第１の代表値や第２の代表値は用いることなく、第１のベクトルＶ１と第２のベクトルＶ２がなす角θの情報のみに基づいて前記判定を行うものである場合、角θそのものを尤度として使用して前記判定を行うことができる。たとえば、判定手段は、角θが所定の閾値より大きい場合、前記一本の枝における血流の方向が分枝ノードに向かう方向であると判定し、角θが閾値を超えない場合、前記一本の枝における血流の方向が分枝ノードから離れる方向であると判定することができる。

また、本発明の画像処理装置において、判定手段は、分枝元の枝から分枝した複数の枝の各々について前記尤度を用いた判定を行った結果、複数の枝における血流の方向が全て分枝ノードに向かう方向であると判定された場合、複数の枝の中でその枝について算出された前記尤度が最も大きい一本の枝における血流の方向は分枝ノードから離れる方向であると判定するものであってもよい。

また、本発明の画像処理装置において、下流側決定手段は、基準点の前後に位置する複数のノードに対応する血管上の複数の位置それぞれにおける血管の径に基づいて、基準点の両側のうち血管の径の平均値が他方より小さい方を血流の下流側として決定するものであってもよいし、グラフ構造の基準点が位置する枝をディスプレイ等の表示装置上に表示させ、基準点の両側のうちいずれか一方を選択するマウスやキーボード等の入力装置による操作者の入力を受け付けて、その選択された方を血流の下流側として決定するものであってもよい。

また、本発明の画像処理装置は、血管を表す画像の表示態様として上流側表示を選択する操作者による入力に応じて、血管を表す画像中の血管を表す領域全体のうちグラフ構造の前記設定された基準点より上流側に位置する部分に対応する部分領域を、他の領域よりも操作者に認識しやすい態様で画像表示装置上に表示させるとともに、表示態様として下流側表示を選択する操作者による入力に応じて、血管を表す画像中の血管を表す領域全体のうちグラフ構造の前記設定された基準点が位置する枝の基準点より下流側の部分と判定手段によりその血流の方向が分枝ノードから離れる方向であると判定された枝の部分に対応する部分領域を、他の領域よりも操作者に認識しやすい態様で画像表示装置上に表示させる表示制御手段を備えたものであってもよい。

ここで、血管を表す領域全体のうち特定の部分領域を他の領域よりも操作者に認識しやすい態様で画像表示装置上に表示させる具体例としては、たとえば、血管を表す領域全体のうち特定の部分領域のみを画像表示装置上に表示させて他の領域は表示対象から除外すること、特定の部分領域を目立つ色で着色した血管を表す領域全体を画像表示装置上に表示させること等がある。

また、本発明の画像処理装置は、血管を表す画像中の血管を表す領域全体のうち、グラフ構造の前記設定された基準点より上流側に位置する部分に対応する第１の部分領域と、グラフ構造の前記設定された基準点が位置する枝の下流側の部分と判定手段によりその血流の方向が分枝ノードから離れる方向であると判定された枝の部分に対応する第２の部分領域とを互いに区別可能に異なる色で画像表示装置上に表示させる表示制御手段を備えたものであってもよい。

この表示制御手段は、血管を表す画像中の血管を表す領域全体のうち、判定手段によりその血流の方向が分枝ノードに向かう方向であると判定された枝の部分に対応する第３の部分領域を、第１の部分領域および第２の部分領域と区別可能に異なる色で画像表示装置上に表示させるものであってもよい。

また、本発明の画像処理方法は、上記画像処理装置の各手段が行う処理を、少なくとも１台のコンピュータにより実行する方法である。

本発明の画像処理プログラムは、上記画像処理方法を少なくとも１台のコンピュータに実行させるプログラムである。このプログラムは、ＣＤ−ＲＯＭ，ＤＶＤなどの記録メディアに記録され、またはサーバコンピュータに付属するストレージやネットワークストレージにダウンロード可能な状態で記録されて、ユーザに提供される。

本発明の画像処理装置、方法およびプログラムによれば、入力装置による操作者の入力に応じて、予め取得しておいたループ状の部分を含む血管の構造を表すグラフ構造上に基準点を設定し、グラフ構造の前記設定された基準点が位置する枝において、基準点より血流の下流側を決定する下流側決定手段と、グラフ構造の前記基準点から前記決定された血流の下流側に辿っていく経路上に位置する任意の分枝ノードにおいて、分枝元の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、分枝ノードに向かう側の向きを有するように規定された第１のベクトルと、分枝元の枝から分枝したいずれか一本の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、分枝ノードから離れる側の向きを有するように規定された第２のベクトルとがなす角θ（0＜θ＜180°）を用いて、その一本の枝における血流の方向が分枝ノードから離れる方向である尤度を角θが小さいほど大きい値となるように算出し、算出された尤度を用いてその一本の枝における血流の方向が分枝ノードから離れる方向であるか分枝ノードに向かう方向であるかを判定することにより、ループ状の部分を含む血管において、分枝元の枝から分枝した部分血管における血流の方向を特定することができる。

上記本発明の画像処理装置、方法およびプログラムにおいて、判定処理が、前記尤度を、分枝元の枝上に位置する複数のノードに対応する血管上の複数の位置それぞれにおける血管の径に基づいて規定された第１の代表値が、前記一本の枝上に位置する複数のノードに対応する血管上の複数の位置それぞれにおける血管の径に基づいて規定された第２の代表値より大きいほど大きい値となるように算出するものである場合には、ループ状の部分を含む血管において、分枝元の枝から分枝した部分血管における血流の方向をより正確に特定することができる。

本発明の実施の形態に係る画像処理装置の概略構成を示す図 判定手段による判定処理を説明するための図 判定手段による第１および第２のベクトルの取得方法を説明するための図 判定手段による第１および第２の代表値の取得方法を説明するための図 判定手段による判定結果の一例を示す図 表示制御手段により表示される下流側表示モードによる画像の一例を示す図 ループ状の部分を含む血管の一例を示す図するための図

以下、本発明の画像処理装置および方法ならびにプログラムの実施形態について説明する。図１は本発明の実施形態となる画像処理装置１の概略構成図である。なお、図１のような画像処理装置１の構成は、補助記憶装置に読み込まれた画像処理プログラムをコンピュータ上で実行することにより実現される。このとき、この画像処理プログラムは、ＣＤ−ＲＯＭ等の記憶媒体に記憶され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされる。図１の画像処理装置１は、たとえば医者等が手術の計画を立てるため、大腸の血管などループ状の部分を含む血管を表す画像を表示させて、血管の走行や血流の方向を観察する場合に、その観察に適した表示を提供するものであって、グラフ抽出手段１０、基準点設定手段２０、下流側決定手段３０、判定手段４０、表示制御手段５０を備えている。以下、本実施の形態では、本発明におけるループ状の部分を含む血管が大腸の血管である場合について説明する。

グラフ抽出手段１０は、大腸の血管が表された３次元画像Ｖから血管を表すグラフ構造を抽出する。具体的には、まず、３次元画像Ｖから大腸の血管の画像的特徴および／または構造的特徴を有する領域を抽出し、その領域を細線化して得られた細線を分岐点および所定の距離などで分割することにより仮のグラフ構造を生成する。ここで、３次元画像Ｖは、データ記憶手段に記憶されたたとえばＣＴ装置、ＭＲＩ装置、超音波診断装置等により取得された多数の２次元画像からなるものである。

次に、生成された仮のグラフ構造を構成する複数の候補点（ノード）の０中から、形状モデル記憶手段ＤＢに記憶された大腸の血管の形状モデルに一致もしくは最も類似する構造を形成する複数の候補点を、形状モデルを構成する複数の教師ラベル（ノード）の各々に対応する対応点として選択し、選択された複数の対応点を形状モデルに略一致するように接続することによりグラフ構造を生成し、抽出する。ここで、形状モデルは、大腸の血管の一般的な形状を表すグラフ構造であって、たとえば多数のサンプル画像を学習することによって取得することができる。

なお、上記対応点の選択は、たとえば下記式（２）で表すコスト関数Ｆの最小化を達成する最適解(対応関係)を求めることにより行う。ここで、Ｒは許容解ｘにおける候補点と教師ラベルの対応づけの集合を表し、θａは集合Ｒに属する任意の対応づけａに対するコストを表し、θabは、集合Ｒに属する２つの対応づけａとｂの組み合わせに対するコストを表す。

基準点設定手段２０は、ディスプレイ等の画像表示装置３上に表示されている３次元画像Ｖに対して、その画像表示装置上の任意の位置を指定するマウスやキーボード等の入力装置２による操作者の入力に応じて、その指定された位置に対応する３次元画像Ｖ上の位置を特定し、特定された３次元画像Ｖ上の位置が前記抽出されたグラフ構造上に有る場合にはその位置に、特定された位置がグラフ構造上にない場合にはその特定された位置から最も近いグラフ構造上の位置に、基準点Ｓを設定する。

下流側決定手段３０は、グラフ構造の前記設定された基準点Ｓの前後に位置する複数のノードに対応する血管上の複数の位置それぞれにおける血管の径に基づいて、基準点の両側のうち血管の径の平均値が他方より小さい方を血流の下流側として決定するものである。

判定手段４０は、グラフ構造を基準点Ｓから前記決定された血流の下流側に辿っていく経路上に位置する任意の分枝ノードＢにおいて、分枝元の枝から分枝した枝のそれぞれにおける血流の方向を判定するものであって、下記式（１）により、その各枝における血流の方向が分枝ノードＢから離れる方向である尤度Ｌを算出し、算出された尤度Ｌが所定の閾値より大きい場合、その枝における血流の方向が分枝ノードＢから離れる方向であると判定し、尤度Ｌが閾値を超えない場合、その枝における血流の方向が分枝ノードＢに向かう方向であると判定する。

上記式（１）において、Ｖ１は分枝元の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、分枝ノードＢに向かう側の向きを有する第１のベクトルを表し、Ｖ２は分枝元の枝から分枝した枝上に位置する複数のノードの位置に基づく方向を有し、かつ、分枝ノードＢから離れる側の向きを有する第２のベクトルを表し、F（Ｖ１，Ｖ２）は第１のベクトルＶ１と第２のベクトルＶ２がなす角θが小さいほど大きい値となる関数を表す。

また、Ｒ１は分枝元の枝上に位置する複数のノードに対応する血管上の複数の位置それぞれにおける血管の径に基づいて規定された代表値（たとえば、平均値）表し、Ｒ２は分枝元の枝から分枝した枝上に位置する複数のノードに対応する血管上の複数の位置それぞれにおける血管の径に基づいて規定された代表値（たとえば、平均値）表し、Ｇ（Ｒ１，Ｒ２）は第１の代表値Ｒ１から第２の代表値Ｒ２を減算した差が大きいほど大きい値となる関数を表す。また、Ｗ１，Ｗ２は上記各関数を重み付けする係数である。

なお、前記閾値としては、第１のベクトルＶ１と第２のベクトルＶ２がなす角θが９０°以下で、かつ、第１の代表値が第２の代表値よりも大きい場合に算出される尤度Ｌの値よりは小さく、第１のベクトルＶ１と第２のベクトルＶ２がなす角θが９０°より大きく、かつ、第１の代表値が第２の代表値よりも小さい場合に算出される尤度Ｌの値よりは大きい値を設定する。

たとえば、図２に示す分枝元の枝Ｅ１から分枝した枝Ｅ２における血流の方向を判定する際、図３に示すように、分枝元の枝Ｅ１上に位置する５つのノードＰ１１〜Ｐ１５の位置に基づく方向を有し、かつ、分枝ノードＢに向かう側の向きを有する第１のベクトルＶ１を取得するとともに、枝Ｅ２上に位置する５つのノードＰ２１〜Ｐ２５の位置に基づく方向を有し、かつ、分枝ノードＢから離れる側の向きを有する第２のベクトルＶ２を取得し、それらの第１のベクトルＶ１および第２のベクトルＶ２を関数Ｆに渡して上記式（１）の第１項の値を求める。このとき、第１項の値は第１のベクトルＶ１および第２のベクトルＶ２がなす角βが小さいほど大きい値となる。

また、図４に示すように、分枝元の枝Ｅ１上に位置する５つのノードＰ１１〜Ｐ１５に対応する血管上の複数の位置それぞれにおける血管の径Ｒ１１〜Ｒ１５の平均値を第１の代表値Ｒ１として取得するとともに、枝Ｅ２上に位置する５つのノードＰ２１〜Ｐ２５に対応する血管上の複数の位置それぞれにおける血管の径Ｒ２１〜Ｒ２５の平均値を第２の代表値Ｒ２として取得し、それらの第１の代表値Ｒ１および第２の代表値Ｒ２を関数Ｇに渡して上記式（１）の第２項の値を求める。このとき、血管の径は、３次元画像Ｖから大腸の血管の画像的特徴および／または構造的特徴を有する血管領域ＲＶを求め、その走行方向に直行する断面の径を求めることにより取得することができる（特開２０１０−２２０７４２公報参照）。そして、それらの各求められた第１項および第２項の値を足し合わせて尤度Ｌを算出し、算出された尤度Ｌを前記閾値と比較して、その枝Ｅ２における血流の方向を判定する。

また、図２に示す分枝元の枝Ｅ１から分枝した枝Ｅ３における血流の方向を判定する際には、図３に示すように、分枝元の枝Ｅ１上に位置する５つのノードＰ１１〜Ｐ１５の位置に基づく方向を有し、かつ、分枝ノードＢに向かう側の向きを有する第１のベクトルＶ１を取得するとともに、枝Ｅ３上に位置する５つのノードＰ３１〜Ｐ３５の位置に基づく方向を有し、かつ、分枝ノードＢから離れる側の向きを有する第２のベクトルＶ３を取得し、それらの第１のベクトルＶ１および第２のベクトルＶ３を（式中の変数Ｖ１に第１のベクトルＶ１が代入され、変数Ｖ２に第２のベクトルＶ３が代入されるように）関数Ｆに渡して上記式（１）の第１項の値を求める。このとき、第１項の値は第１のベクトルＶ１および第２のベクトルＶ３がなす角αが小さいほど大きい値となる。

また、図４に示すように、分枝元の枝Ｅ１上に位置する５つのノードＰ１１〜Ｐ１５に対応する血管上の複数の位置それぞれにおける血管の径Ｒ１１〜Ｒ１５の平均値を第１の代表値Ｒ１として取得するとともに、枝Ｅ３上に位置する５つのノードＰ３１〜Ｐ３５に対応する血管上の複数の位置それぞれにおける血管の径Ｒ３１〜Ｒ３５の平均値を第２の代表値Ｒ３として取得し、それらの第１の代表値Ｒ１および第２の代表値Ｒ３を（式中の変数Ｒ１に第１の代表値Ｒ１が代入され、変数Ｒ２に第２の代表値Ｒ３が代入されるように）関数Ｇに渡して上記式（１）の第２項の値を求める。そして、それらの各求められた第１項および第２項の値を足し合わせて尤度Ｌを算出し、算出された尤度Ｌを前記閾値と比較して、その枝Ｅ２における血流の方向を判定する。

また、判定手段４０は、分枝元の枝から分枝した複数の枝の各々について前記尤度Ｌを用いた判定を行った結果、複数の枝における血流の方向が全て分枝ノードＢに向かう方向であると判定された場合、複数の枝の中でその枝について算出された前記尤度Ｌが最も大きい一本の枝における血流の方向は分枝ノードＢから離れる方向であると判定する機能を備えている。

たとえば、図２に示す分枝元の枝Ｅ１から分枝した２本の枝Ｅ２，Ｅ３の各々について前記尤度Ｌを用いた判定を行った結果、両方の枝における血流の方向が全て分枝ノードＢに向かう方向であると判定された場合、それらの枝Ｅ２，Ｅ３のうちその枝について算出された尤度Ｌが大きい方の枝における血流の方向は分枝ノードＢから離れる方向であると判定する。

表示制御手段５０は、ディスプレイ等の画像表示装置３における３次元画像Ｖ等の血管を表す画像の表示を制御するものであって、たとえば図５に示すような判定手段４０による前記判定の結果に基づいて、画像中の大腸の血管を表す領域を、前記グラフ構造の基準点Ｓより上流側に位置する部分に対応する第１の部分領域ＲＵと、グラフ構造の基準点Ｓが位置する枝の基準点Ｓより下流側の部分と判定手段４０によりその血流の方向が分枝ノードＢから離れる方向であると判定された枝の部分に対応する第２の部分領域ＲＤと、判定手段４０によりその血流の方向が分枝ノードＢに向かう方向であると判定された枝の部分に対応する第３の部分領域ＲＢの３つの領域に区分し（図６参照）、第１の部分領域ＲＵと第２の部分領域ＲＤの２つの領域を、あるいは、第１の部分領域ＲＵと第２の部分領域ＲＤ、および第３の部分領域ＲＢの３つの領域を互いに区別可能に異なる色で表した画像を画像表示装置３上に表示させる。

また、表示制御手段５０は、画像中の大腸の血管を表す領域全体のうち、第１の部分領域ＲＵを第２の部分領域ＲＤや第３の部分領域ＲＢよりも操作者に認識しやすい態様で表示させる上流側表示モードと、第２の部分領域ＲＤを第１の部分領域ＲＵや第３の部分領域ＲＢよりも操作者に認識しやすい態様で表示させる下流側表示モードとを有し、マウスやキーボード等の入力装置２によりそのいずれか一方の表示モードを選択する操作者の入力に応じて、その選択された表示モードによる画像を画像表示装置３上に表示させる機能を備えている。たとえば、下流側表示モードを選択する操作者による入力を受けて、図６に示すような、第１の部分領域ＲＵや第３の部分領域ＲＢは表示対象から除外し、第２の部分領域ＲＤのみを表示した画像を画像表示装置３上に表示させる。

さらに、表示制御手段５０は、第２の部分領域ＲＤと、第３の部分領域ＲＢの分枝ノードＢから所定の範囲内に位置する部分領域との集合である第４の部分領域を、他の領域よりも操作者に認識しやすい態様で表示させる機能も備えている。この第４の部分領域は、基準点Ｓから下流側の血管部分を切除する切除手術における切除すべき領域に対応しており、第４の部分領域が他の血管を表す領域と接続している位置は、切除手術を行う際に血流をクリップで止めておくべき位置に対応しているため、このような表示を提供することにより、医者等の操作者による切除手術を計画作業を支援し、その作業効率を向上させることができる。

上記構成により、画像処理装置１においては、まず、操作者による要求に応じてデータ記憶手段に記憶されている３次元画像Ｖが読み出され、グラフ抽出手段１０が３次元画像Ｖから血管を表すグラフ構造を抽出し、表示制御手段５０が３次元画像Ｖを画像表示装置３に表示させる。その後、基準点設定手段２０が３次元画像Ｖ上の任意の位置を指定する操作者の入力を受け付け、その入力により指定された位置に対応するグラフ構造上の位置に基準点Ｓを設定する。

その後、後述する方法により基準点Ｓが位置する枝における血流の方向が特定されると、判定手段４０が、グラフ構造を基準点Ｓから下流側に辿りながらその経路上に位置する分枝ノードにおいて順次、分枝元の枝から分枝した枝のそれぞれにおける血流の方向を判定する。判定手段４０は、その判定により血流の方向が分枝ノードから離れる方向である枝を新たに特定し、その枝を通るように前記経路を延長させ、延長された経路上に位置する分枝ノードにおいてさらに判定を行う、という処理を繰り返し行う。これにより、基準点Ｓからグラフ構造の末梢部までの経路上に位置する全ての分枝ノードにおいて、分枝元の枝から分枝した枝のそれぞれにおける血流の方向を判定する。

その後、表示制御手段５０が、判定手段４０による判定の結果に基づいて、画像中の血管を表す領域を第１の部分領域ＲＵ、第２の部分領域ＲＤおよび第３の部分領域ＲＢの３つの領域に区分し、たとえばそれらの領域を互いに区別可能に異なる色で表した画像など、初期表示として予め設定しておいた表示態様による画像を画像表示装置３上に表示させる。表示制御手段５０は、表示中の画像とは異なる表示態様による画面表示を選択する操作者による入力がなされたかをさらに検出し、当該入力がなされた場合には、その選択された表示モードによる画像に表示を切り替える。

ここで、基準点Ｓが位置する枝における血流の方向は、基準点Ｓから前後にそれぞれ５点〜１０点のノードにおける血管の径を調べて、血管の径が細くなっていく方向をその枝における血流の方向として特定するようにしてもよいし、操作者にその枝のどちらか上流かの情報を入力させ、その情報に基づいて特定するようにしてもよい。

なお、上記実施の形態では、判定手段４０が分枝の角度と血管の径の両方の情報に基づいて前記判定を行うものである場合について説明したが、分枝の角度の情報のみに基づいて前記判定を行うようにしてもよい。たとえば図２に示す分枝元の枝Ｅ１から分枝した枝Ｅ２における血流の方向を判定する際、第１のベクトルＶ１と第２のベクトルＶ２とがなす角βが９０度より大きい場合、その枝Ｅ２における血流の方向が分枝ノードＢに向かう方向であると判定し、角θが９０度を超えない場合、その枝Ｅ２における血流の方向が分枝ノードから離れる方向であると判定することができる。

また、上記実施の形態では、グラフ抽出手段１０が３次元画像Ｖを画像解析することにより血管を表すグラフ構造を直接抽出するものである場合について説明したが、そのグラフ構造が既に抽出されて記録媒体に記録されている場合には、別途抽出するまでもなく、その記録媒体からグラフ構造のデータを読み出して以降の処理に提供すればよい。

１ 画像処理装置
１０ グラフ抽出手段
２０ 基準点設定手段
３０ 下流側決定手段
４０ 判定手段
５０ 表示制御手段
Ｓ 基準点
Ｂ 分枝ノード
Ｅ１ 分枝元の枝
Ｅ２，Ｅ３ 分枝した枝
Ｖ１ 第１のベクトル
Ｖ２，Ｖ３ 第２のベクトル
ＲＵ 第１の部分領域
ＲＤ 第２の部分領域
ＲＭ 第３の部分領域

Claims (12)

1. 入力装置による操作者の入力に応じて、予め取得しておいたループ状の部分を含む血管の構造を表すグラフ構造上に基準点を設定する基準点設定手段と、
   前記グラフ構造の前記設定された基準点が位置する枝において、前記基準点より血流の下流側を決定する下流側決定手段と、
   前記グラフ構造の前記基準点から前記決定された血流の下流側に辿っていく経路上に位置する任意の分枝ノードにおいて、分枝元の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、前記分枝ノードに向かう側の向きを有するように規定された第１のベクトルと、前記分枝元の枝から分枝したいずれか一本の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、前記分枝ノードから離れる側の向きを有するように規定された第２のベクトルとがなす角θを用いて、前記一本の枝における血流の方向が前記分枝ノードから離れる方向である尤度を前記なす角θが小さいほど大きい値となるように算出し、該算出された尤度を用いて前記一本の枝における血流の方向が前記分枝ノードから離れる方向であるか前記分枝ノードに向かう方向であるかを判定する判定手段と
   を備えたことを特徴とする画像処理装置。
2. 前記判定手段が、前記尤度を、前記分枝元の枝上に位置する複数のノードに対応する前記血管上の複数の位置それぞれにおける前記血管の径に基づいて規定された第１の代表値が、前記一本の枝上に位置する複数のノードに対応する前記血管上の複数の位置それぞれにおける前記血管の径に基づいて規定された第２の代表値より大きいほど大きい値となるように算出するものであることを特徴とする請求項１記載の画像処理装置。
3. 前記判定手段が、前記分枝元の枝から分枝した複数の枝の各々について前記尤度を用いた判定を行った結果、前記複数の枝における血流の方向が全て前記分枝ノードに向かう方向であると判定された場合、前記複数の枝の中でその枝について算出された前記尤度が最も大きい一本の枝における血流の方向は前記分枝ノードから離れる方向であると判定するものであることを特徴とする請求項１または２記載の画像処理装置。
4. 前記判定手段が、前記角θが所定の閾値より大きい場合、前記一本の枝における血流の方向が前記分枝ノードに向かう方向であると判定し、前記角θが前記閾値を超えない場合、前記一本の枝における血流の方向が前記分枝ノードから離れる方向であると判定するものであることを特徴とする請求項１記載の画像処理装置。
5. 前記下流側決定手段が、前記基準点の前後に位置する複数のノードに対応する前記血管上の複数の位置それぞれにおける前記血管の径に基づいて、前記基準点の両側のうち前記血管の径の平均値が他方より小さい方を前記血流の下流側として決定するものであることを特徴とする請求項１から４のいずれか１項記載の画像処理装置。
6. 前記血管を表す画像の表示態様として上流側表示を選択する操作者による入力に応じて、前記血管を表す画像中の前記血管を表す領域全体のうち前記グラフ構造の前記設定された基準点より上流側に位置する部分に対応する部分領域を、他の領域よりも前記操作者に認識しやすい態様で前記画像表示装置上に表示させるとともに、前記表示態様として下流側表示を選択する操作者による入力に応じて、前記血管を表す画像中の前記血管を表す領域全体のうち前記グラフ構造の前記設定された基準点が位置する枝の前記基準点より下流側の部分と前記判定手段によりその血流の方向が前記分枝ノードから離れる方向であると判定された枝の部分に対応する部分領域を、他の領域よりも前記操作者に認識しやすい態様で前記画像表示装置上に表示させる表示制御手段を備えたことを特徴とする請求項１から５のいずれか１項記載の画像処理装置。
7. 前記血管を表す画像中の前記血管を表す領域全体のうち、前記グラフ構造の前記設定された基準点より上流側に位置する部分に対応する第１の部分領域と、前記グラフ構造の前記設定された基準点が位置する枝の下流側の部分と前記判定手段によりその血流の方向が前記分枝ノードから離れる方向であると判定された枝の部分に対応する第２の部分領域とを互いに区別可能に異なる色で前記画像表示装置上に表示させる表示制御手段を備えたことを特徴とする請求項１から５のいずれか１項記載の画像処理装置。
8. 前記表示制御手段が、前記血管を表す画像中の前記血管を表す領域全体のうち、前記判定手段によりその血流の方向が前記分枝ノードに向かう方向であると判定された枝の部分に対応する第３の部分領域を、前記第１の部分領域および前記第２の部分領域と区別可能に異なる色で前記画像表示装置上に表示させるものであることを特徴とする請求項７記載の画像処理装置。
9. 入力装置による操作者の入力に応じて、予め取得しておいたループ状の部分を含む血管の構造を表すグラフ構造上に基準点を設定する処理と、
   前記グラフ構造の前記設定された基準点が位置する枝において、前記基準点より血流の下流側を決定する処理と、
   前記グラフ構造の前記基準点から前記決定された血流の下流側に辿っていく経路上に位置する任意の分枝ノードにおいて、分枝元の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、前記分枝ノードに向かう側の向きを有するように規定された第１のベクトルと、前記分枝元の枝から分枝したいずれか一本の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、前記分枝ノードから離れる側の向きを有するように規定された第２のベクトルとがなす角θを用いて、前記一本の枝における血流の方向が前記分枝ノードから離れる方向である尤度を前記なす角θが小さいほど大きい値となるように算出し、該算出された尤度を用いて前記一本の枝における血流の方向が前記分枝ノードから離れる方向であるか前記分枝ノードに向かう方向であるかを判定する処理と
   を少なくとも１台のコンピュータにより実行することを特徴とする画像処理方法。
10. 前記判定処理が、前記尤度を、前記分枝元の枝上に位置する複数のノードに対応する前記血管上の複数の位置それぞれにおける前記血管の径に基づいて規定された第１の代表値が、前記一本の枝上に位置する複数のノードに対応する前記血管上の複数の位置それぞれにおける前記血管の径に基づいて規定された第２の代表値より大きいほど大きい値となるように算出するものであることを特徴とする請求項９記載の画像処理方法。
11. コンピュータを、
    入力装置による操作者の入力に応じて、予め取得しておいたループ状の部分を含む血管の構造を表すグラフ構造上に基準点を設定する基準点設定手段と、
    前記グラフ構造の前記設定された基準点が位置する枝において、前記基準点より血流の下流側を決定する下流側決定手段と、
    前記グラフ構造の前記基準点から前記決定された血流の下流側に辿っていく経路上に位置する任意の分枝ノードにおいて、分枝元の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、前記分枝ノードに向かう側の向きを有するように規定された第１のベクトルと、前記分枝元の枝から分枝したいずれか一本の枝上に位置する複数のノードの位置に基づく方向を有し、かつ、前記分枝ノードから離れる側の向きを有するように規定された第２のベクトルとがなす角θを用いて、前記一本の枝における血流の方向が前記分枝ノードから離れる方向である尤度を前記なす角θが小さいほど大きい値となるように算出し、該算出された尤度を用いて前記一本の枝における血流の方向が前記分枝ノードから離れる方向であるか前記分枝ノードに向かう方向であるかを判定する判定手段として機能させるための画像処理プログラム。
12. 前記判定手段が、前記尤度を、前記分枝元の枝上に位置する複数のノードに対応する前記血管上の複数の位置それぞれにおける前記血管の径に基づいて規定された第１の代表値が、前記一本の枝上に位置する複数のノードに対応する前記血管上の複数の位置それぞれにおける前記血管の径に基づいて規定された第２の代表値より大きいほど大きい値となるように算出するものであることを特徴とする請求項１１記載の画像処理プログラム。