WO2015136853A1

WO2015136853A1 - 画像処理装置、画像処理方法及びプログラム

Info

Publication number: WO2015136853A1
Application number: PCT/JP2015/000874
Authority: WO
Inventors: 真透坂本
Original assignee: テルモ株式会社
Priority date: 2014-03-14
Filing date: 2015-02-23
Publication date: 2015-09-17
Also published as: JPWO2015136853A1

Abstract

　管状物体を断層撮影して得られた情報に基づいて、管状物体から分岐する分岐管部についての情報を自動的に提供する。　プローブを用いて第１の管状物体内を走査することにより得られた管状物体画像のうち、第１の管状物体に相当する部位と、第１の管状物体から分岐する第２の管状物体に相当する部位と、を区別する情報を取得する。情報を使用して、第１の管状物体からの分岐部分における第２の管状物体の形態を示す定量的な情報を生成する。

Description

画像処理装置、画像処理方法及びプログラム

　本発明は、画像処理装置、画像処理方法及びプログラムに関するものである。

　従来より、血管等の管状物体の画像診断のために、管状物体内にプローブを挿入して断層撮影を行う装置が知られている。例えば、動脈硬化の診断や、バルーンカテーテル又はステント等の高機能カテーテルによる血管内治療時の術前診断、あるいは、術後の結果確認のために、このような画像診断装置が広く利用されている。代表的な画像診断装置としては、例えば血管内超音波診断装置（ＩＶＵＳ）、光干渉断層画像診断装置（ＯＣＴ／ＯＦＤＩ）等が開発されている。

　これらの装置により撮影された血管画像を、医師が容易に確認できるようにする技術もまた開発されている。引用文献１には、ＩＶＵＳ又はＯＣＴにより得られた断層像から再構成された３次元画像を、ＣＴ又はＭＲＩにより得られた３次元画像と対応づける技術が開示されている。引用文献１に記載の技術においては、対応付けのために、ＩＶＵＳ又はＯＣＴにより得られた画像から、血管の分岐部が自動的に抽出される。

特開２０１２－０７５７０２号公報

　プローブを挿入した管状物体についての情報だけでなく、管状物体から分岐する分岐管部についての情報を確認することが望まれる場合がある。例えば、血管の分岐部にステントを留置する場合には、主血管の情報のみならず、主血管から分岐している分岐血管の情報を確認することが望まれる。しかしながら、分岐管部の情報を容易に取得することを可能とする技術は開発されていなかった。

　本発明は、管状物体を断層撮影して得られた情報に基づいて、管状物体から分岐する分岐管部についての情報を自動的に提供することを目的とする。

　本発明の課題を解決するために、例えば、本発明の一実施形態に係る画像処理装置は以下の構成を備える。すなわち、
　プローブを用いて第１の管状物体内を走査することにより得られた管状物体画像のうち、前記第１の管状物体に相当する部位と、前記第１の管状物体から分岐する第２の管状物体に相当する部位と、を区別する情報を取得する情報取得手段と、
　前記情報を使用して、前記第１の管状物体からの分岐部分における前記第２の管状物体の形態を示す定量的な情報を生成する生成手段と、
　を備える。

　管状物体を断層撮影して得られた情報に基づいて、管状物体から分岐する分岐管部についての情報を自動的に提供することができる。

　本発明のその他の特徴及び利点は、添付図面を参照とした以下の説明により明らかになるであろう。なお、添付図面においては、同じ若しくは同様の構成には、同じ参照番号を付す。

　添付図面は明細書に含まれ、その一部を構成し、本発明の実施の形態を示し、その記述と共に本発明の原理を説明するために用いられる。
図１は本発明の一実施形態に係る画像処理装置１００の概略構成を示す図である。図２は一実施形態における境界抽出部２００の概略構成を示す図である。図３は一実施形態において境界抽出部２００が行う処理のフローチャートである。図４Ａは一実施形態において境界抽出部２００が行う処理を説明する図である。図４Ｂは一実施形態において境界抽出部２００が行う処理を説明する図である。図４Ｃは一実施形態において境界抽出部２００が行う処理を説明する図である。図４Ｄは一実施形態において境界抽出部２００が行う処理を説明する図である。図４Ｅは一実施形態において境界抽出部２００が行う処理を説明する図である。図４Ｆは一実施形態において境界抽出部２００が行う処理を説明する図である。図５Ａは一実施形態において判定更新部２５０が行う処理を説明する図である。図５Ｂは一実施形態において判定更新部２５０が行う処理を説明する図である。図５Ｃは一実施形態において判定更新部２５０が行う処理を説明する図である。図６Ａは一実施形態において判定更新部２５０が行う処理を説明する図である。図６Ｂは一実施形態において判定更新部２５０が行う処理を説明する図である。図６Ｃは一実施形態において判定更新部２５０が行う処理を説明する図である。図７Ａは一実施形態において判定更新部２５０が行う処理を説明する図である。図７Ｂは一実施形態において判定更新部２５０が行う処理を説明する図である。図７Ｃは一実施形態において判定更新部２５０が行う処理を説明する図である。図８Ａは一実施形態において判定更新部２５０が行う処理を説明する図である。図８Ｂは一実施形態において判定更新部２５０が行う処理を説明する図である。図９は境界抽出部２００による判定結果を表示する画面の一例を示す図である。図１０は一実施形態における情報算出部１０００の概略構成を示す図である。図１１は一実施形態において情報算出部１０００が行う処理のフローチャートである。図１２Ａは主血管と分岐血管との境界点群の近似方法の一例を説明する図である。図１２Ｂは主血管と分岐血管との境界点群の近似方法の一例を説明する図である。図１２Ｃは主血管と分岐血管との境界点群の近似方法の一例を説明する図である。図１２Ｄは主血管と分岐血管との境界点群の近似方法の一例を説明する図である。図１２Ｅは主血管と分岐血管との境界点群の近似方法の一例を説明する図である。図１３は分岐血管の分岐角度の定義の一例を説明する図である。図１４Ａは分岐血管の断面へと投影される点の選択方法を説明する図である。図１４Ｂは分岐血管の断面へと投影される点の選択方法を説明する図である。図１４Ｃは分岐血管の断面へと投影される点の選択方法を説明する図である。図１４Ｄは分岐血管の断面へと投影される点の選択方法を説明する図である。図１５は情報算出部１０００により生成された情報の表示画面の一例を示す図である。図１６は一実施形態におけるコンピュータの基本構成を示す図である。

　以下、本発明の実施形態を図面に基づいて説明する。ただし、本発明の範囲は以下の実施形態に限定されるものではない。

　図１は、本発明の一実施形態に係る画像処理装置１００を示す。画像処理装置１００は、表示制御部１１０と、境界抽出部２００と、情報算出部１０００とを備える。画像処理装置１００は、血管画像１９０を取得する撮影装置（不図示）を備える撮影システムの一部であってもよい。

　画像処理装置１００には、プローブを用いて管状物体内を走査することにより得られた管状物体画像が入力される。管状物体の種類は特に限定されないが、以下では管状物体は血管であり、画像処理装置１００には血管画像１９０が入力されるものとして説明する。

　血管画像１９０は、血管の形態を示す画像情報であって、その形態は特に限定されない。血管画像の取得方法は特に限定されず、例えば既存の超音波診断装置（ＩＶＵＳ）、光干渉断層画像診断装置（ＯＣＴ）、又は光干渉断層画像診断装置（ＯＦＤＩ）等を用いることができる。以下では、プローブが挿入された血管のことを主血管と呼ぶ。

　血管画像１９０は、以下のようにして取得することができる。すなわち、血管にプローブを挿入して血管内を走査することにより、複数のラインデータが取得される。それぞれのラインは１走査に対応し、１つのラインデータは、プローブから血管の深さ方向への距離と得られた信号強度との関係を示す。プローブの向き及びプローブの血管長さ方向の位置を変えながら走査を行うことにより、複数のラインデータを得ることができる。一実施形態において、血管画像１９０はこの複数のラインデータにより構成されている。例えば、血管画像１９０は、ラインデータにより示される１次元の画像を横に並べたものであってもよい。また、血管画像１９０は、ラインデータにより示される１次元の画像を円状に並べることにより得られる複数の横断断面画像（血管軸を横断する方向の断面像、例えば血管軸に垂直な断層像）により構成されていてもよい。また、血管画像１９０に対しては、ゲイン補正、コントラスト補正、又はγ補正等が行われていてもよい。本実施形態においては、血管画像１９０は、ゲイン及びコントラストが調整された、複数の横断断面画像で構成されている。

　血管画像１９０は複数の横断断面画像の組み合わせには限られず、複数の血管軸方向断面画像（血管軸と平行な断面での断層像）の組み合わせであってもよいし、血管の３Ｄ画像であってもよい。既に知られているように、これらの横断断面画像、血管軸方向断面画像、及び３Ｄ画像は、互いに変換することが可能である。

　境界抽出部２００は、血管画像１９０から、主血管に相当する部分を検出する。また、境界抽出部２００は、血管画像１９０から、主血管から分岐している分岐血管に相当する部分を検出する。このようにして、境界抽出部２００は、主血管と分岐血管との境界部分を検出する。検出結果に従って、境界抽出部２００は、血管画像１９０のうち主血管に相当する部位を示す情報と分岐血管に相当する部位とを区別する情報とを、情報算出部１０００へと出力する。もっとも、後述するように、主血管と分岐血管との境界部分を検出することに先立って、主血管に相当する部分と分岐血管に相当する部分とを検出することは必須ではない。境界抽出部２００の詳しい構成については後述する。ここで、分岐血管とは、プローブが挿入された主血管に接続している構造体であって、管状のものだけでなく、瘤状のものをも含む。

　情報算出部１０００は、境界抽出部２００から得た情報を参照して、主血管からの分岐部分における分岐血管の形態を示す定量的な情報を生成する。そして、情報算出部１０００は、生成された情報を表示制御部１１０を介して表示部１２０に表示させる。

　表示制御部１１０は、表示部１２０を制御して所望の情報を表示させる。表示部１２０は、画像を表示可能な装置である。表示部１２０の種類は特に限定されず、例えば液晶ディスプレイ等でありうる。

［境界抽出部］
　まず、境界抽出部２００の概略構成について、図２を参照して説明する。境界抽出部２００は、画像取得部２１０、位置取得部２２０、検出部２３０、及び判定部２４０を備える。

　画像取得部２１０は、血管画像１９０を取得する。本実施形態において、画像取得部２１０は複数の断層像を取得する。それぞれの断層像に対して、後述する各部による処理が行われ、主血管に相当する部位を示す情報と分岐血管に相当する部位とを区別する情報が生成される。また、画像取得部２１０は、取得した断層像に対する前処理も行う。

　位置取得部２２０は、断層像上における主血管の中心点の推定位置を取得する。本実施形態において、位置取得部２２０は、断層像に対する画像処理により、主血管の中心点の推定位置を算出する。しかしながら、入力部（不図示）を介したユーザ入力により、主血管の中心点が指定されてもよい。

　一実施形態において、位置取得部２２０は、断層像上におけるガイドワイヤ及びカテーテルシースの像の推定位置を取得する。位置取得部２２０は、断層像に対する画像処理により、これらの推定位置を算出する。しかしながら、入力部（不図示）を介したユーザ入力により、これらの位置が指定されてもよい。

　検出部２３０は、断層像上における血管の内表面を示す複数の点を検出する。検出方法としては、テンプレートマッチング等の様々な方法が考えられるが、本実施形態において、検出部２３０は、断層像を走査することにより断層像上における血管内壁の位置を抽出する。具体的には、検出部２３０は、主血管の中心点から外側に進んだ際における血管との交差点を、血管内壁の位置として抽出する。

　判定部２４０は、抽出部２３０によって検出された複数の点のそれぞれが、主血管を示しているか否かを判定する。また、判定部２４０は、抽出部２３０によって検出された交差点が、分岐血管を示しているか否かを判定する。また、判定部２４０は、抽出部２３０によって検出された複数の点のそれぞれが、主血管と分岐血管との境界を示すか、を判定することもできる。これらの判定は、抽出部２３０によって検出された複数の点の位置に基づいて行われる。一実施形態において、これらの判定は、抽出部２３０によって検出された、連続した複数の点同士の位置関係に基づいて行われる。例えば、後述するように、連続した複数の点についての主血管の中心からの距離の変化に基づいて、抽出部２３０によって検出された点が主血管を示すか及び分岐血管を示すかを判定することができる。後述するように、抽出部２３０によって検出された連続した複数の点についての主血管の中心からの距離と閾値との比較に基づいて、それぞれの点が主血管を示すか及び分岐血管を示すかを判定することも、連続した複数の点同士の位置関係に基づく判定に含まれる。また、連続した複数の点が途切れている箇所には分岐血管が存在すると判定することができ、連続した複数の点の末端は主血管と分岐血管との境界を示すと判定することができる。

　また、これらの判定は、独立に行うことができる。すなわち、判定部２４０は、抽出部２３０によって検出された点が主血管を示すか及び分岐血管を示すか、又は抽出部２３０によって検出された点が主血管と分岐血管との境界を示すか、のうちの少なくとも一方を判定することができる。これらの判定により、主血管と分岐血管との境界が検出される。

　判定部２４０による判定には、主血管の中心点と交差点との距離が主に用いられる。しかしながら、判定部２４０は、他の様々な基準を用いてこの判定を行う。例えば、一実施形態において、判定部２４０は断層像上におけるガイドワイヤ及びカテーテルシースの像の位置を考慮して判定を行う。

　一実施形態において、境界抽出部２００は判定更新部２５０をさらに備える。判定更新部２５０は、判定部２４０による判定結果を更新する。例えば、判定更新部２５０は、判定部２４０による判定結果を修正することができる。また、判定更新部２５０は、新たな交差点を検出し、検出された交差点が主血管又は分岐血管を示すと判定することができる。また、判定更新部２５０は、血管画像１９０を構成する１つの断層像についての判定部２４０による判定結果に基づいて、別の断層像についての判定部２４０による判定結果を更新することができる。

　次に、図３のフローチャートを参照して、境界抽出部２００が行う処理について詳しく説明する。

　ステップＳ３０５において、画像取得部２１０は、血管画像１９０を取得する。本実施例において、画像取得部２１０は、複数の断層像で構成される血管画像１９０を一度に取得するものとする。しかしながら、画像取得部２１０は、プローブを用いた血管の走査を行いながら順次生成された断層像を、順次取得してもよい。以下のステップＳ３１０以降の処理は、それぞれの断層像に対して順次行われる。以下の説明において処理対象となる断層像を図４Ａに示す。

　ステップＳ３１０において、画像取得部２１０は、取得した断層像に対する前処理を行う。もっとも、画像取得部２１０は既に前処理が行われている断層像を取得してもよい。前処理の種類は特に限定されないが、一実施形態においては、フィルタ処理及び２値化処理が行われる。用いられるフィルタの種類は特に限定されないが、例えば平滑化フィルタ処理を用いることができ、断層像に対して平滑化フィルタ処理を行うことにより、ノイズとして含まれる孤立点を低減することができる。フィルタ処理後の断層像を図４Ｂに示す。

　また、所定の閾値を用いて２値化処理を行うことにより、後述する処理を簡略化することができる。２値化処理においては、血管内壁が白画素となり、血管内腔が黒画素となるような閾値が用いられる。このような閾値は、予め設定されていてもよいし、断層像のヒストグラムを参照して画像取得部２１０により自動的に決定されてもよい。また、閾値は入力部（不図示）を介してユーザにより入力されてもよい。以下では、ステップＳ３１０において２値化処理を行ったことを前提として説明する。しかしながら、２値化処理を行うことは必須ではない。２値化処理を行う場合は、閾値により規定される所定範囲内の画素値を有する画素が白画素であると判定され、所定範囲外の画素値を有する画素は黒画素であると判定される。また、２値化処理を行わない場合であっても、例えば後述するステップＳ３２５において、所定範囲内の画素値を有する画素が白画素であると判定し、所定範囲外の画素値を有する画素は黒画素であると判定することができる。２値化後の断層像を図４Ｃに示す。

　ステップＳ３１５において、位置取得部２２０は、断層像上における主血管の中心点の推定位置を取得する。前述したように、中心点の推定位置はユーザにより入力されてもよいが、本実施形態において、位置取得部２２０は、断層像に対する画像処理により中心点の推定位置を算出する。

　中心点の推定位置の算出方法は特に限定されない。一実施形態においては、断層像上の主血管が検出され、その重心が主血管の中心点の推定位置として用いられる。本実施形態においては、主血管の検出方法としてはハフ変換が用いられる。すなわち、位置取得部２２０は、２値化後の断層像に対してハフ変換処理を行うことにより、主血管の内壁形状を近似する円を検出する。そして、位置取得部２２０は、検出された円の中心位置を中心点の推定位置として取得する。ハフ変換により検出された円の例を図４Ｄに示す。

　ハフ変換を用いた中心点の推定方法について、以下に詳しく説明する。まず、位置取得部２２０は、黒画素の円（輪郭の内側が黒画素であり、輪郭の外側（背景）が白画素である円）をハフ変換により抽出する。そして、位置取得部２２０は、断層像の中心を内包する円のうち最大のものを、主血管の内壁を示す円（以下、血管円と呼ぶ）として特定する。この血管円の中心が、主血管の中心点の推定位置として扱われる。また、位置取得部２２０は、血管円の半径も算出する。

　もっとも、主血管の内腔が略円形状ではない場合には、ハフ変換により血管円を抽出できない可能性がある。この場合、位置取得部２２０は、主血管の他の位置における断層像から抽出された血管円の中心を、中心点の推定位置として用いることができる。血管円の半径についても同様である。一実施形態においては、ハフ変換により血管円を抽出できた断層像であって、最も近い位置の断層像から抽出された血管円の中心が、中心点の推定位置として用いられる。別の実施形態においては、隣接する位置の断層像における中心点の推定位置が、中心点の推定位置として用いられる。

　さらなる方法として、ハフ変換により血管円を抽出できない場合、位置取得部２２０は血管の内接円を検出し、この内接円の中心を主血管の中心点の推定位置として扱ってもよい。例えば、２値化後の断層像上において、３点以上の白画素に内接し、白画素を通らない円を、血管の内接円として検出することができる。また例えば、ユークリッド距離変換を用いて、血管に内接する円のうち最大の円を検出することができる。距離変換はユークリッドに限らず、また重み付けを行ってもよい。さらに別の方法として、抽出部２３０によって検出された複数の点の重心を、主血管の中心点の推定位置として扱ってもよい。

　後述するように断層像における主血管を示す交差点が決定された後に、主血管を示す交差点の重心位置を算出することもできる。この場合、位置取得部２２０は、主血管の他の位置における断層像について算出された主血管を示す交差点の重心位置を、中心点の推定位置として用いてもよい。さらなる実施形態において、位置取得部２２０は、ハフ変換により血管円を抽出できなかった場合に、主血管の他の位置における断層像について主血管を示す交差点の重心位置を算出し、中心点の推定位置として用いることができる。

　このように、主血管の他の位置における断層像から抽出された血管円の中心を、中心点の推定位置として用いる場合、交差点の誤判定が生じる可能性がある。そこで、このような断層像については、判定部２４０は、隣接する断層像について検出されたワイヤ影境界ペアを参照して、ワイヤ影境界ペアとなる交差点を決定してもよい。具体的な一例としては、判定部２４０は、隣接する断層像について検出されたワイヤ影境界ペアと同じ角度方向に、ワイヤ影境界ペアが存在すると判定することができる。また、判定部２４０は、隣接する２つの断層像について検出されたワイヤ影境界のそれぞれの角度方向の中央に、ワイヤ影境界が存在すると判定することができる。

　ステップＳ３２０において、位置取得部２２０は、断層像からガイドワイヤ及びカテーテルシースの像の推定位置をさらに検出する。通常、血管画像を撮影する際には、ガイドワイヤを介してガイディングカテーテル（不図示）が主血管内へと挿入され、カテーテルシースで覆われているプローブはガイディングカテーテル内へと挿入される。このため、断層像には、主血管に挿入されたガイドワイヤ及びカテーテルシースの像が含まれている。本実施形態においては、判定部２４０は断層像上におけるガイドワイヤ及びカテーテルシースの像の位置を考慮して判定を行う。

　具体的には、位置取得部２２０は、白画素の円（輪郭の内側が白画素であり、輪郭の外側が黒画素である円）をハフ変換により抽出する。そして、位置取得部２２０は、断層像の中心に最も近い円を、カテーテルシースを示す円（以下、シース円と呼ぶ）として特定する。位置取得部２２０は、抽出された円のうち、主血管の内壁を示す円の内側にあるもののみを、カテーテルシースを示す円の候補として扱ってもよい。このシース円の中心が、カテーテルシース像の推定位置として扱われる。

　また、位置取得部２２０は、２値化前の断層像において、シース円と血管円との間にある画素のうち、最高輝度を有する画素を検出する。検出された画素は、ガイドワイヤの一部を示す画素（以下、ガイドワイヤ画素と呼ぶ）として扱われ、ガイドワイヤ画素の位置が、ガイドワイヤ像の推定位置として扱われる。

　ステップＳ３２５において、検出部２３０は、主血管の中心点から外側に進んだ際における血管との交差点を、血管内壁の位置として抽出する。以下に、その具体的な処理について説明する。

　まず、検出部２３０は、２値化された断層像に対する前処理を行う。具体的には、検出部２３０は、血管円の中心から所定距離内にある画素を全て黒画素とする。この所定距離は、血管円の半径以下となるように定められる。この処理により、血管円の中心に近く、血管内壁ではない画素が、誤って血管内壁として検出される可能性を低くすることができる。一方で、この所定距離を血管円の半径よりも小さくすることにより、血管円付近に位置する血管内壁を示す画素の検出に失敗する可能性を低くすることができる。この所定距離は特に限定されず、例えば血管円の半径の０．９５倍以下であってもよいし、血管円の半径の０．５０倍以上であってもよい。

　この処理により、例えば、ガイドワイヤ像、カテーテルシース像、及び血管内腔のノイズ等を、断層像から消すことができる。このため、主血管の中心点からの距離が所定距離以下である領域からは交差点が検出されなくなる。ガイドワイヤ像、カテーテルシース像、及び血管内腔のノイズを消すためには、別の方法を用いることもできる。例えば、テンプレートマッチングを用いてガイドワイヤ像又はカテーテルシース像等を認識し、認識した像を消すことができる。さらに、形態学的処理を用いることにより、ノイズ又はカテーテルシース像等を認識して消してもよい。交差点検出の直前の断層像を図４Ｅに示す。

　その後、検出部２３０は、主血管の中心点から各角度方向へと走査を行い、血管との交差点、すなわち白画素との交差点を検出する。言い換えれば、検出部２３０は、主血管の中心点から各角度方向に伸びる複数の半直線のそれぞれについて、半直線上の白画素のうち中心点に最も近い画素を、この角度方向に対応する交差点として検出する。３６０°の走査を行うことにより、血管内壁を示す画素が抽出される。検出部２３０は、検出されたそれぞれの交差点に、角度方向と、主血管の中心点からの距離と、を対応づけて保存する。保存された情報を用いて、横軸が走査角度を示し、縦軸が主血管の中心点からの距離を示すグラフを作成することができる。

　検出部２３０による検出の例を図４Ｆに示す。図４Ｆにおいては、血管円の中心４００から、各角度方向４１１，４１２，４１３のそれぞれについて走査が行われている。その結果、交差点４２１，４２２，４２３が検出されている。

　ステップＳ３３０において、判定部２４０は、主血管の中心点から交差点までの距離に基づいて、交差点が主血管を示すかを判定する。また、判定部２４０は、交差点が分岐血管を示すかを判定する。一般的に、断層像において分岐血管は主血管から離れていくので、分岐血管の内壁は主血管の内壁と比べて主血管の中心点からより離れていると考えられる。そこで、判定部２４０は、主血管の中心点から交差点までの距離が閾値以下である場合に、この交差点が主血管の内壁を示すと判定する。この閾値は、血管円の半径に基づいて決定される。具体的な閾値の決定方法は特に限定されず、例えば血管円の半径の０．９０倍以上であってもよいし、血管円の半径の１．１０倍以下であってもよい。

　また、一実施形態においては、判定部２４０は、主血管の中心点から交差点までの距離が閾値を超える場合に、この交差点が分岐血管の内壁を示すと判定する。しかしながら、プローブから見てガイドワイヤの後ろに位置する主血管の内壁は断層像に写らない。そして、図４Ｆに示す写らない領域４３０（以下、ガイドワイヤの影と呼ぶ）が分岐血管として誤認識される可能性がある。そこで、本実施形態において、判定部２４０は、さらにカテーテルシースの像及びガイドワイヤの像の位置を考慮して、主血管の中心点からの距離が閾値を超える交差点が、分岐血管を示すかガイドワイヤの影を示すのかを判定する。

　具体的には、判定部２４０は、シース円の中心からガイドワイヤ画素の位置への延長上付近にある交差点はガイドワイヤの影を示し、この延長線上付近にはない交差点は分岐血管を示すと判定する。例えば、判定部２４０は、主血管の中心点と交差点との間の距離が閾値を超え、かつ断層像上でカテーテルシースの像からガイドワイヤの像へと向かう方向から所定の角度範囲内に交差点が存在する場合に、この交差点はガイドワイヤの影を示すと判定する。また、判定部２４０は、主血管の中心点と交差点との間の距離が閾値を超え、かつ断層像上でカテーテルシースの像からガイドワイヤの像へと向かう方向から所定の角度範囲内にない交差点は分岐血管を示すと判定する。この角度範囲は特に限定されず、適宜設定することができる。

　判定部２４０による判定方法はこの方法には限定されない。例えば、判定部２４０は、隣接する角度方向に位置する２つの交差点について、主血管の中心点から交差点までの距離が大きく変化している場合に、この２つの交差点はガイドワイヤの影又は分岐血管を示すと判定することができる。そして、判定部２４０は、ガイドワイヤの影又は分岐血管を示すと判定されなかった交差点は、主血管を示すと判定することができる。この方法は、例えば、ハフ変換により血管円が検出できなかった断層像に対して適用することができる。

　ステップＳ３３５において、判定部２４０は、判定結果に従って、主血管と分岐血管との境界を示す交差点を検出する。また、判定部２４０は、主血管とガイドワイヤの影との境界を示す交差点を検出する。本実施形態においては、隣接する角度方向にある交差点について、一方が主血管を示し、他方が分岐血管を示すと判定されている場合、判定部２４０は、主血管を示す交差点が、主血管と分岐血管との境界を示すと判定する。もっとも、この場合、分岐血管を示す交差点が、主血管と分岐血管との境界を示すと判定してもよい。検出部２３０による角度方向への走査の解像度が十分に高ければ、どちらの方法を採用しても得られる情報に大差はない。また、隣接する角度方向にある交差点について、一方が主血管を示し、他方がガイドワイヤの影を示すと判定されている場合、判定部２４０は、主血管を示す交差点が、主血管とガイドワイヤの影との境界を示すと判定する。もっとも、この場合、ガイドワイヤの影を示す交差点が、主血管とガイドワイヤの影との境界を示すと判定してもよい。

　また、主血管の内壁が途切れている部分に、分岐血管又はガイドワイヤの影を示す交差点が検出されない可能性もある。例えば、図４Ｆに示す場合、ガイドワイヤの影が位置する角度方向において、ガイドワイヤの影を示す交差点は検出されないかもしれない。また、主血管の内壁と分岐血管の内壁とがなす角度によっては、分岐位置の周辺で分岐血管を示す交差点が検出されない可能性がある。そこで、判定部２４０は、断層像上において主血管の内表面が途切れている箇所の位置に基づいて、内表面を示す点が主血管と分岐血管（又はガイドワイヤの影）との境界に位置するかを判定することができる。

　具体的な一例としては、判定部２４０は、隣接する角度方向について、一方の方向に主血管を示す交差点が存在するが、他方の方向には主血管を示す交差点が存在しない場合に、この交差点は主血管と分岐血管又はガイドワイヤの影との境界を示すと判定することができる（これを境界候補点とする）。より具体的には、判定部２４０は、カテーテルシースの像からガイドワイヤの像へと向かう方向から所定の角度範囲内に境界候補点が存在する場合、この境界候補点は主血管とガイドワイヤの影との境界を示すと判定する。また、判定部２４０は、カテーテルシースの像からガイドワイヤの像へと向かう方向から所定の角度範囲内に境界候補点が存在しない場合、この境界候補点は主血管と分岐血管との境界を示すと判定する。このような方法により、判定部２４０は、抽出部２３０によって検出された複数の点のそれぞれが、主血管と分岐血管との境界を示すかを判定することができる。本実施形態においては、抽出部２３０によって検出された点が主血管と分岐血管との境界を示すかの判定は、抽出部２３０によって検出された点が主血管を示すか及び分岐血管を示すかを判定した後に行われる。しかしながら、抽出部２３０によって検出された点が主血管と分岐血管との境界を示すかの判定は、抽出部２３０によって検出された点が主血管を示すか及び分岐血管を示すかの判定とは独立に行うことができる。

　ステップＳ３４０において、判定部２４０は、分岐血管部分を挟む、主血管と分岐血管との境界を示す交差点ペアを検出する。具体的には、判定部２４０は、主血管と分岐血管との境界を示し、異なる２つの角度方向に存在する２つの交差点の組であって、この２つの角度方向の間に存在する交差点が全て分岐血管を示す（又は交差点が存在しない）交差点の組を検出する。以下では、こうして検出された交差点のペアのことを、分岐部境界ペアと呼ぶ。また、判定部２４０は、主血管とガイドワイヤの影との境界を示し、異なる２つの角度方向に存在する２つの交差点の組であって、この２つの角度方向の間に存在する交差点が全てガイドワイヤの影を示す（又は交差点が存在しない）交差点の組を検出する。以下では、こうして検出された交差点ペアのことを、ワイヤ影境界ペアと呼ぶ。このような検出された分岐部境界ペア及びワイヤ影境界ペアは、判定更新部２５０による処理のために用いられる。

　ステップＳ３４５において、判定更新部２５０は、ステップＳ３３０における判定部２４０による判定結果を更新する。ここで、判定更新部２５０は、主血管と分岐血管との境界を示す交差点の位置と、主血管とガイドワイヤの影との境界を示す交差点の位置と、の少なくとも一方を参照する。例えば、判定更新部２５０は、これらの境界を示す交差点の位置に基づいて定まるパラメータが所定の条件に合致するか否かを判定し、合致する場合に、判定結果を修正することができる。また、判定更新部２５０は、これらの境界を示す交差点の位置に基づいて新たな交差点を検出し、検出された交差点に対して判定を行うことができる。以下では判定結果の更新方法について５つの例を説明するが、判定結果の更新方法には様々な形態が考えられ、以下の例には限定されない。

（１）三日月状主血管についての誤検出
　主血管が図５Ａに示す形状（三日月形状）を有している場合、主血管の血管内壁の一部について交差点が検出されないか、又は分岐血管を示す交差点として検出される可能性がある。例えば図５Ａの場合、交差点５０２は主血管の中心点からの距離が大きい、もしくはその距離の変化が大きいため、分岐血管を示す交差点として検出される可能性がある。このような場合には、主血管の中心点から、分岐部境界ペア５０１，５０３のそれぞれへの距離の差が大きくなることが考えられる。一方で、図５Ｂに示すように、主血管が略円形状であり、分岐血管を示す交差点５１２が検出された場合には、主血管の中心点から、分岐部境界ペア５１１，５１３のそれぞれへの距離は同程度であることが多い。

　そこで、本実施形態において、判定更新部２５０は、主血管と分岐血管との境界を示す交差点の位置を参照して、所定の条件に適合する場合、判定部２４０によって分岐血管を示すと判定された交差点は、主血管を示すと判定する。そして、判定更新部２５０は、判定部２４０による判定を更新する。具体的には、判定更新部２５０は、主血管の中心点から分岐部境界ペアのそれぞれまでの距離の差が閾値以上である場合に、判定部２４０によって分岐血管を示すと判定された交差点は、主血管を示すと判定する。図５Ａの場合であれば、判定更新部２５０は、判定部２４０によって分岐血管を示すと判定された交差点５０２は、主血管を示すと判定する。また、判定更新部２５０は、交差点５０１，５０３は主血管と分岐血管との境界を示さず、分岐部境界ペアではないと判定する。

　主血管が図５Ａに示す形状を有している場合、主血管の中心点から各角度方向へと走査を行っても、交差点５０２付近で検出される、主血管を示す交差点の数は少なくなる。そこで、本実施形態において、判定更新部２５０は、分岐部境界ペアではないと判定された交差点５０１，５０３の間に存在する、主血管を示す交差点を追加検出する。追加検出方法の一例を、図５Ｃを参照して説明する。例えば、判定更新部２５０は、分岐部境界ペアではないと判定された交差点５０１，５０３を通る直線５０２を、主血管の中心点から遠ざかる方向に所定距離だけ平行移動することにより得られる直線５０５を設定する。そして、判定更新部２５０は、平行移動された直線５０５と血管との交差点５０４，５０６を主血管を示す交差点として検出する。具体的には、判定更新部２５０は、交差点５０１，５０３の中点５０２ｍを、交差点５０１，５０３を通る直線５０２と垂直に、主血管の中心点から遠ざかる方向に所定距離だけ移動した点５０５ｓを基準点として設定する。そして、判定更新部２５０は、設定された基準点５０５ｓから両方向へと直線５０５上で走査を行い、血管との交差点のペア、すなわち白画素との交差点のペア５０４，５０６を検出する。

　こうして検出された交差点５０４，５０６は、主血管を示す交差点である可能性がある。そこで、判定更新部２５０は、上流側又は下流側の横断断面画像における血管の大きさを考慮して、検出された交差点５０４，５０６が主血管を示すとの判定が妥当であるか否かを判定する。一実施形態においては、略円形である上流側又は下流側の横断断面画像が参照される。例えば、図５Ａに示す形状を有するために主血管を示す交差点が追加検出された横断断面画像は、この判定においては考慮されない。判定方法の一例として、交差点５０４，５０６が主血管を示す交差点であると判定を更新することにより、更新された血管の大きさが上流側又は下流側の横断断面画像における血管のおおよその大きさに近づく場合には、判定が妥当であると判定することができる。判定更新部２５０は、妥当であると判定した場合、交差点５０４，５０６は主血管を示すと、判定部２４０による判定を更新する。

　判定更新部２５０は、このように新たに検出された主血管を示す交差点のペアを通る直線の設定、平行移動、及び走査を行うことにより、順次主血管を示す交差点を検出することができる。すなわち、判定更新部２５０は、交差点５０４，５０６を通る直線５０５を、主血管の中心点から遠ざかる方向に所定距離だけ平行移動することにより得られる直線５０８を設定する。また、判定更新部２５０は、交差点５０４，５０６の中点５０５ｍを、交差点５０４，５０６を通る直線５０５と垂直に、主血管の中心点から遠ざかる方向に所定距離だけ移動した点５０８ｓを基準点として設定する。そして、判定更新部２５０は、設定された基準点５０８ｓから両方向へと直線５０８上で走査を行い、血管との交差点のペア、すなわち白画素との交差点のペア５０７，５０９を検出する。判定更新部２５０は、検出された交差点５０７，５０９を、主血管を示す交差点と判定する。この処理を繰り返した際に基準点が白画素となった場合には、繰り返しを終了してもよいし、平行移動させる距離をより小さくしてもよい。

（２）ひょうたん状主血管についての誤検出
　主血管が狭窄しており、図６Ｂに示す形状（ひょうたん形状）を有している場合にも、主血管の血管内壁の一部について交差点が検出されないか、又は分岐血管を示す交差点として検出される可能性がある。例えば図６Ｂに示す断層像６２０の場合、交差点６２２，６２３は主血管の中心点からの距離が大きいため、分岐血管を示す交差点として検出される可能性がある。この場合、交差点６２１，６２４は、分岐部境界ペアであると判定される。

　このように、狭窄部位付近では複数の断層像について分岐血管を示す交差点が誤検出される可能性がある。狭窄部位を前後の断層像と比較してみると、一般的には、上流側から狭窄部位に近づくにつれて主血管が狭くなり、狭窄部位から下流側に進むにつれて主血管は広くなる。図６Ｂは、狭窄が最も高度な部位における断層像６２０を示し、６２５は、断層像６２０における分岐部境界ペア６２１，６２４間の距離を示す。一方で、図６Ａは、断層像６２０よりも上流における断層像６１０を示し、６１２は、断層像６１０における分岐部境界ペア６１１，６１３間の距離を示す。また、図６Ｃは、断層像６２０よりも下流における断層像６３０を示し、６３２は、断層像６３０における分岐部境界ペア６３１，６３３間の距離を示す。このように、分岐部境界ペア間の距離は、上流側から狭窄部位に近づくにつれて短くなり、狭窄部位から下流側に進むにつれて長くなる。一方で、主血管が略円形状であって分岐血管が存在する場合、分岐部境界ペア間の距離は、上流側から分岐血管の中心に近づくにつれて長くなり、分岐血管の中心から下流側に近づくにつれて短くなる。

　そこで、本実施形態では、判定更新部２５０は、主血管と分岐血管との境界を示す交差点の位置を参照して、所定の条件に適合する場合、判定部２４０によって分岐血管を示すと判定された交差点は、主血管を示すと判定する。具体的には、判定更新部２５０は、断層像６１０から検出された分岐部境界ペア間の距離６１２が、断層像６２０から検出された分岐部境界ペア間の距離６２５よりも大きいか否かを判定する。また、判定更新部２５０は、断層像６３０から検出された分岐部境界ペア間の距離６３２が、断層像６２０から検出された分岐部境界ペア間の距離６２５よりも大きいか否かを判定する。このどちらの条件も満たされる場合、断層像６１０，６２０，６３０において、判定部２４０によって分岐血管を示すと判定された交差点は、主血管を示す可能性がある。また、断層像６１０と断層像６３０との間の断層像において、判定部２４０によって分岐血管を示すと判定された交差点についても、主血管を示す可能性がある。

　そして、判定更新部２５０は、上流側又は下流側の横断断面画像における血管の大きさを考慮して、検出された交差点６２２，６２３が主血管を示すとの判定が妥当であるか否かを判定する。この判定は、図５Ａ～Ｃを参照して説明した交差点５０４，５０６についての判定と同様に行うことができる。判定更新部２５０は、妥当であると判定した場合、交差点６２２，６２３は主血管を示すと、判定部２４０による判定を更新する。また、判定更新部２５０は、交差点６２１，６２４は主血管と分岐血管との境界を示さず、分岐部境界ペアではないと判定する。

　主血管が図６Ｂに示す形状を有している場合も、主血管の中心点から各角度方向へと走査を行っても、交差点６２２，６２３付近で検出される、主血管を示す交差点の数は少なくなる。そこで、本実施形態において、判定更新部２５０は、分岐部境界ペアではないと判定された交差点６２１，６２４の間に存在する、主血管を示す交差点を追加検出する。追加検出は、主血管が図５Ａに示す形状を有している場合と同様に行うことができる。

（３）分岐血管とワイヤ影とが重なっている場合の誤検出
　分岐血管とガイドワイヤの影とが重なっている場合、ガイドワイヤの位置によっては、本来は分岐血管を示す交差点が、ガイドワイヤの影を示すと誤判定される可能性がある。例えば図７Ｂに示す断層像７２０において、交差点７２３は分岐血管を示すが、ガイドワイヤの影を示すと誤判定されている。図７Ｂにおいては、交差点７２２もガイドワイヤの影を示すと判定されており、交差点７２１，７２４がワイヤ影境界ペアとして判定されている。このような誤判定が行われた場合、ワイヤ影境界ペア７２１，７２４により示されるガイドワイヤの影の大きさは大きくなる。

　一方で、互いに隣接する位置の断層像において、通常ガイドワイヤの影の大きさはあまり異ならない。図７Ａは、図７Ｂに示す断層像７２０の近隣の位置における断層像７１０である。また、図７Ｃも、図７Ｂに示す断層像７２０の近隣の位置における断層像７３０である。断層像７１０及び断層像７３０では、分岐血管とガイドワイヤの影とは重なっていない。この場合、ワイヤ影境界ペア７１１，７１２により示されるガイドワイヤの影の大きさは、ワイヤ影境界ペア７３１，７３２によるガイドワイヤの影の大きさとは大きく異ならない。

　そこで、判定更新部２５０は、主血管とガイドワイヤの影との境界を示す交差点の位置を参照して、判定部２４０によってガイドワイヤの影を示すと判定された交差点は、分岐血管を示すと判定する。具体的には、判定更新部２５０は、近隣の位置の断層像と比較してガイドワイヤの影の大きさが大きい場合に、ガイドワイヤの影と分岐血管とが重なっていると判定し、ガイドワイヤの影を示すと判定された交差点のうち一部について、分岐血管を示すと判定する。

　以下に、判定更新部２５０が行う具体的な処理の一例を説明する。まず判定更新部２５０は、検出されたワイヤ影境界ペアについて、角度幅及び中心角度を算出する。角度幅とは、主血管の中心点からワイヤ影境界ペアの一方への角度方向と、主血管の中心点からワイヤ影境界ペアの他方への角度方向との差分を表す。また、中心角度とは、主血管の中心点からワイヤ影境界ペアの一方への角度方向と、主血管の中心点からワイヤ影境界ペアの他方への角度方向と、の中央に位置する角度方向を表す。図７Ａ，Ｃにおいて、７１３，７３３は角度幅を示し、７１４，７３４は中心角度を示す。

　そして、判定更新部２５０は、明らかにワイヤ影境界ペアの角度幅が広い断層像を検出する。具体的には、判定更新部２５０は、ワイヤ影境界ペアの角度幅が閾値を超えている断層像を検出する。閾値は、予め設定されていてもよいし、判定更新部２５０によって算出されてもよい。例えば、この閾値は、検出対象となっている断層像に隣り合っている所定数の断層像について算出された角度幅の平均値に基づいて算出することができる。一例においては、閾値は、こうして算出された平均値に所定値を加算したものでありうる。以下では、図７Ｂに示す断層像７２０について、判定更新部２５０は、角度幅７２５が閾値を超えていると判定したものとする。

　ワイヤ影境界ペアの角度幅が閾値を超えている断層像７２０が検出されると、判定更新部２５０は、この断層像７２０においてガイドワイヤの影が位置する部分、すなわちワイヤ影境界ペアの位置を推定する。判定更新部２５０は、この推定を、処理対象となる断層像の近隣に位置する断層像７１０，７３０について算出されたワイヤ影境界ペアの位置に基づいて行う。具体的には、判定更新部２５０は、断層像７１０から検出されたワイヤ影境界ペア７１１，７１２の位置と、断層像７３０から検出されたワイヤ影境界ペア７３１，７３２の位置とを参照して、第２の断層像におけるワイヤ影境界ペアの位置を推定する。ここで、断層像７１０、断層像７２０、及び断層像７３０はそれぞれ血管の第１、第２、及び第３の位置における断層像であり、第２の位置は第１の位置と第３の位置との間に存在する。

　より具体的には、判定更新部２５０は、断層像７１０から検出されたワイヤ影境界ペアの中心角度７１４と、断層像７３０から検出されたワイヤ影境界ペアの中心角度７３４との線形補間によって値を算出することができる。この場合、判定更新部２５０は、算出された値を、断層像７２０についてのワイヤ影境界ペアの中心角度７２７と推定することができる。また、判定更新部２５０は、断層像７１０から検出されたワイヤ影境界ペアの角度幅７１３と、断層像７３０から検出されたワイヤ影境界ペアの角度幅７３３との線形補間によって値を算出することができる。この場合、判定更新部２５０は、算出された値を、断層像７２０についてのワイヤ影境界ペアの角度幅７２８と推定することができる。ここで、線形補間は距離に応じて重み付けされていてもよいし、また３つ以上のフレームからの非線形補間であってもよい。

　こうして推定されたワイヤ影境界ペアの中心角度７２７及び角度幅７２６を用いて、判定更新部２５０は、断層像７２０においてガイドワイヤの影を示すと判定された交差点のうち一部を選択する。具体的には、判定更新部２５０は、断層像７２０においてガイドワイヤの影を示すと判定された交差点のうち、推定されたワイヤ影境界ペアの中心角度７２７及び角度幅７２６により特定される領域７２８内に存在しない交差点７２３を選択することができる。もっとも、判定更新部２５０は、推定されたワイヤ影境界ペアの角度幅７２６の算出誤差を考慮して、この角度幅７２６に対して所定の角度値を加算してから、交差点７２３の選択を行ってもよい。そして、判定更新部２５０は、選択された交差点７２３は分岐血管を示すと判定する。また、判定更新部２５０は、新たな判定に従って、主血管と分岐血管との境界点を示す交差点７２４を、判定部２４０と同様に検出することができる。

　判定更新部２５０は、ワイヤ影境界ペアの位置を正確に推定するために、分岐血管とガイドワイヤの影とが重なっていない断層像を、推定に用いる断層像として選択することができる。すなわち、一実施形態において、検出部２３０が検出した断層像７１０におけるワイヤ影境界ペアの角度幅７１３は閾値以下であり、かつ断層像７３０におけるワイヤ影境界ペアの角度幅７３３は閾値以下である。

（４）分岐血管の分岐起始部付近における追加検出
　検出部２３０は所定の角度方向解像度に従って主血管の中心点から各角度方向に進んだ際における血管との交差点を検出する。しかしながら、解像度が粗い場合、分岐血管の分岐起始部付近、すなわち分岐血管の中心から離れた位置においては、分岐血管を抽出できない可能性がある。そこで、本実施形態において、判定更新部２５０は、主血管と分岐血管との境界を示す交差点の位置を参照して、検出されていない分岐血管が存在する可能性がある部分を検出する。そして、判定更新部２５０は、検出された部分についてより高い角度方向解像度に従って交差点の検出を行い、検出された交差点が分岐血管を示すか否かを判定する。

　具体的には、まず、判定更新部２５０は、分岐部境界ペアが検出された断層像８１０と、分岐部境界ペアが検出されなかった断層像８２０と、の組み合わせを検出する。ここで、断層像８１０は、血管の第１の位置における断層像であり、断層像８２０は、第１の位置に隣接する第２の位置における断層像である。分岐血管は互いに隣接する複数の断層像に写るため、断層像８２０からは分岐血管は検出されなかったものの分岐血管が写っている可能性がある。

　次に、判定更新部２５０は、図８Ａに示すように断層像８１０から算出された分岐部境界ペア８１１，８１２の角度方向８１４を参照して、図８Ｂに示すように断層像８２０における検出範囲８２５を設定する。この検出範囲８２５は、検出開始角度及び検出終了角度により規定される。具体的には、検出範囲８２５は、角度方向８１４を含むように設定される。設定される検出範囲８２５の広さは特に制限されない。一実施形態においては、検出範囲８２５の広さは、分岐部境界ペア８１１，８１２の角度幅８１３と一致するように設定される。すなわち、検出開始角度は主血管の中心点から交差点８１１への角度方向に一致し、検出終了角度は主血管の中心点から交差点８１２への角度方向に一致する。別の実施形態においては、主血管の中心点から交差点までの距離の変化のトレンドを確認するために、検出範囲８２５の広さは、分岐部境界ペア８１１，８１２の角度幅８１３よりも所定値だけ大きくなるように設定される。分岐部境界ペアの角度幅及び角度方向は、ワイヤ影境界ペアの角度幅及び角度方向と同様に算出できる。

　そして、判定更新部２５０は、設定された検出範囲８２５において、検出部２３０と同様の方法を用いながらも、より感度の高い検出方法を使用して、交差点を検出する。以下に、交差点の検出感度を高める例について説明するが、感度の高め方は以下の例に限られるわけではない。

　一例として、判定更新部２５０は、検出部２３０が走査を行った際の角度方向解像度よりも高い角度方向解像度に従って、交差点を検出することができる。一例において、検出部２３０が１°ごとに走査を行うのであれば、判定更新部２５０は０．２５°ごとに走査を行うことができる。そして、判定更新部２５０は、判定部２４０と同様に、検出された交差点が主血管を示すのか分岐血管を示すのかを判定する。

　また、上述のように、判定部２４０は、隣接する角度方向に位置する２つの交差点について、主血管の中心点から交差点までの距離が大きく変化している場合に、この２つの交差点は分岐血管（又はガイドワイヤの影）を示すと判定することができる。判定更新部２５０も、同様の方法を用いて交差点が分岐血管を示すと判定することができるが、この場合に用いる判定の閾値を、判定部２４０が判定を行う場合の閾値よりも小さくすることができる。例えば、判定部２４０は、主血管の中心点から２つの交差点までの距離差が第１の閾値よりも大きい場合にこの２つの交差点は分岐血管（又はガイドワイヤの影）を示すと判定することができる。また、判定更新部２５０は、主血管の中心点から２つの交差点までの距離差が第２の閾値よりも大きい場合にこの２つの交差点は分岐血管（又はガイドワイヤの影）を示すと判定することができる。ここで、第２の閾値は第１の閾値よりも小さい。

　別の方法においては、隣接する角度方向に位置する交差点について、主血管の中心点からの距離の変化のトレンドが検出され、距離の変化が閾値よりも大きくなる交差点が分岐血管を示すと判定される。例えば、隣接する角度方向に位置する第１、第２、及び第３交差点について、主血管の中心点から第１及び第２の交差点までの距離差と、主血管の中心点から第２及び第３の交差点までの距離差と、の差分が計算される。判定部２４０は、この差分が第１の閾値よりも大きい場合に、第１の交差点は分岐血管（又はガイドワイヤの影）を示すと判定することができる。また、判定更新部２５０は、この差分が第２の閾値よりも大きい場合に、第１の交差点は分岐血管を示すと判定することができる。ここでも、第２の閾値は第１の閾値よりも小さい。

（５）分岐血管の追加抽出
　検出部２３０が用いる、血管の中心点から各角度方向へと走査することにより交差点を検出する方法では、分岐血管部分においては検出される交差点の数が少なくなることが考えられる。そこで、判定更新部２５０は、主血管と分岐血管との境界を示す交差点の位置、具体的には分岐部境界ペアの位置を参照して、分岐血管を示すさらなる交差点を検出し、検出された交差点について分岐血管を示すと判定する。

　具体的な方法は、図５Ｃを参照して説明した、主血管を示す交差点を追加検出する方法と同様である。すなわち、判定更新部２５０は、分岐部境界ペアであると判定された交差点を通る直線を、主血管の中心点から遠ざかる方向に所定距離だけ平行移動することにより得られる直線を設定する。そして、判定更新部２５０は、平行移動された直線と血管との交差点を、分岐血管を示す交差点として検出する。

　さらなる実施形態において、判定更新部２５０は、分岐部境界ペアを形成していない、主血管と分岐血管との境界を示す交差点の位置を参照して、分岐血管を示すさらなる交差点を検出し、検出された交差点について分岐血管を示すと判定する。例えば、図７Ｂに示されるように分岐血管とガイドワイヤの影とが重なっている場合、分岐部境界ペアを形成していない、主血管と分岐血管との境界を示す交差点が生じうる。この場合、判定更新部２５０は、主血管と分岐血管との境界を示す交差点７２４を通り、推定されたワイヤ影境界ペアの中心角度７２７により示される方向に直交する直線を、主血管の中心点から遠ざかる方向に所定距離だけ平行移動することにより得られる直線を設定する。そして、判定更新部２５０は、平行移動された直線と血管との交差点を、分岐血管を示す交差点として検出することができる。

　この場合、判定更新部２５０は、交差点７２４から、中心角度７２７により示される方向に直交する直線に沿って、主血管の中心点へと近づく方向に所定距離離れた点を決定する。そして、判定更新部２５０は、決定された点を、中心角度７２７により示される方向に、主血管の中心点から遠ざかる方向に所定距離だけ移動した点を基準点として設定する。さらに、判定更新部２５０は、設定された基準点から交差点７２４に近づく方向へと設定された直線上で走査を行い、血管との交差点、すなわち白画素との交差点を、分岐血管を示す交差点として検出することができる。

　別の実施形態においては、判定更新部２５０は、分岐部境界ペアであると判定された交差点の代わりに、主血管とガイドワイヤの影との境界を示す交差点７２１と、主血管と分岐血管との境界を示す交差点７２４と、のペアを用いることもできる。すなわち、判定更新部２５０は、このペアを通る直線を、主血管の中心点から遠ざかる方向に所定距離だけ平行移動することにより得られる直線を設定する。そして、判定更新部２５０は、平行移動された直線と血管との交差点のペアのうち、交差点７２４に近い方の交差点を、分岐血管を示す交差点として検出する。

　判定更新部２５０は、同様の方法により、ワイヤ影境界ペアの位置を参照して、ガイドワイヤの影を示すさらなる交差点を検出し、検出された交差点についてガイドワイヤの影を示すと判定してもよい。

　判定更新部２５０は、以上の（１）～（５）の方法のうち１つのみを用いてもよいし、複数の方法を組み合わせて用いてもよい。また、以上の方法を適用する順序は特に制限されない。一実施形態においては、分岐血管を示す交差点を多く抽出するために、判定更新部２５０は、（１）～（４）の処理を行った後に、（５）の処理を行う。方法（２）～（４）では、処理対象となる断層像とは異なる位置での断層像における、主血管と分岐血管との境界を示す交差点の位置と、主血管とガイドワイヤの影との境界を示す交差点の位置と、の少なくとも一方を参照して処理が行われる。このように、複数の断層像を参照することにより、より高い精度で主血管と分岐血管との境界を検出することができる。

　境界抽出部２００は、こうして得られた判定結果を、分岐血管についての定量的な情報の算出のために、情報算出部１０００へと送る。具体的には、それぞれの交差点についての位置情報及び判定部２４０又は判定更新部２５０による判定結果が情報算出部１０００へと送られる。一実施形態において、この位置情報には、主血管の中心点からの角度方向及び中心点からの距離が含まれる。この場合、境界抽出部２００は、断層像における主血管の中心点の推定位置を、情報算出部１０００に送ることができる。別の実施形態において、この位置情報には、断層像上における交差点のＸＹ座標が含まれていてもよい。上述したように、境界抽出部２００は断層像ごとに複数の交差点についての判定結果を生成し、情報算出部１０００へと送る。

　もっとも、全ての判定結果を送ることは必須ではなく、境界抽出部２００は、少なくとも、主血管に相当する部位と、主血管から分岐する分岐血管に相当する部位と、を区別する情報を、情報算出部１０００へと送る。例えば、境界抽出部２００は、主血管と分岐血管との境界を示す交差点についての位置情報のみを、情報算出部１０００へと送ってもよい。また、境界抽出部２００は、分岐血管を示す交差点についての位置情報のみを、情報算出部１０００へと送ってもよい。

　また、境界抽出部２００は、表示制御部１１０を介して、判定結果を表示部１２０に表示させることができる。表示画面の例を図９に示す。図９に示す表示画面９００には、横断断面画像９１０及び血管軸方向断面画像９２０が表示されている。横断断面画像９１０には、検出された交差点を示し、互いに異なる色で表示されたドット９１１～９１５が重畳表示されている。それぞれのドット９１１～９１５は、主血管を示す交差点、分岐血管を示す交差点、ガイドワイヤの影を示す交差点、主血管と分岐血管との境界点を示す交差点、及び主血管とガイドワイヤの影との境界点を示す交差点を表している。ユーザは、例えばチェックボックス９１６を操作することにより、これらのドット９１１～９１５の表示を停止することができる。血管軸方向断面画像にも、同様にドット９１１～９１５が重畳表示されている。このように、判定部２４０又は判定更新部２５０により自動的に生成された判定結果を表示することにより、ユーザは血管画像のうち分岐血管又はガイドワイヤの影に相当する部分を容易に識別することが可能となる。

［情報算出部］
　次に、情報算出部１０００の概略構成について、図１０を参照して説明する。情報算出部１０００は、情報取得部１０１０及び生成部１０２０を備える。

　情報取得部１０１０は、血管画像１９０のうち、主血管に相当する部位と、主血管から分岐する分岐血管に相当する部位と、を区別する情報を取得する。上述したように、情報取得部１０１０は、境界抽出部２００により得られた判定結果を取得する。本実施形態において、情報取得部１０１０は、少なくとも、主血管と分岐血管との境界を示す交差点についての位置情報、又は分岐血管を示す交差点についての位置情報を、境界抽出部２００から取得する。情報取得部１０１０は、さらに、血管画像１９０を取得することもできる。

　生成部１０２０は、情報取得部１０１０が取得した情報を使用して、主血管からの分岐部分における分岐血管の形態を示す定量的な情報を生成する。この定量的な情報としては、例えば、主血管からの分岐部分における分岐血管の大きさを示す情報、及び主血管からの分岐血管の分岐方向を示す情報等が挙げられる。本実施形態において、生成部１０２０は、分岐血管の大きさを示す情報と分岐方向を示す情報との双方を生成するが、生成部１０２０は、このうち一方の情報を生成してもよい。

　次に、図１１のフローチャートを参照して、情報算出部１０００が行う処理について詳しく説明する。

　ステップＳ１１１０において、情報取得部１０１０は、上述したように、血管画像１９０のうち、主血管に相当する部位と、主血管から分岐する分岐血管に相当する部位と、を区別する情報を取得する。

　ステップＳ１１２０において、生成部１０２０は、血管画像１９０に含まれる分岐血管の分岐部分を検出する。生成部１０２０は、分岐血管が存在する連続した断層像について、１つの分岐部分が存在するものと判定する。断層像に分岐血管が存在するか否かは、主血管と分岐血管との境界を示す交差点、又は分岐血管を示す交差点の存在を確認することにより、容易に行うことができる。複数の分岐部分が検出された場合、以下の処理はそれぞれの分岐部分について行われる。

　ステップＳ１１３０～Ｓ１１４０において、生成部１０２０は、ステップＳ１１２０で検出された分岐部分についての分岐面情報を算出する。分岐面とは、主血管と分岐血管とを接続する面であり、分岐面情報の例としては、分岐面の位置、向き及び大きさ等が挙げられる。分岐面の大きさは、主血管からの分岐部分における分岐血管の大きさに相当する。分岐面の大きさには、分岐面の面積だけではなく、分岐面の特徴を示す長さ（円であれば半径や直径、楕円であれば短径や長径など）及び分岐面の周囲の長さも含まれる。

　ステップＳ１１３０において、生成部１０２０は、主血管と分岐血管との境界点群を近似する近似多角形又は近似楕円を導出する。例えば、図１２Ａには、主血管と分岐血管との境界を示す交差点群である境界点群１２０１が表されている。本実施形態において、生成部１０２０は、図１２Ｂに示すように、この境界点群１２０１を近似する近似平面１２０２を導出する。さらに、生成部１０２０は、図１２Ｃに示すように、それぞれの境界点群１２０１を近似平面１２０２へと垂直に投影することにより、投影点群１２０３を得る。こうして、図１２Ｄに示すように、投影点群１２０３により構成される近似多角形１２０４が導出される。近似多角形１２０４は隣り合う投影点群１２０３を線形に結んで得られる。しかしながら、近似多角形１２０４の代わりに、投影点群１２０３を例えばスプライン曲線等を用いて非線形に結んで得られる近似曲線を導出してもよい。そして、生成部１０２０は、図１２Ｅに示すように、近似多角形１２０４を近似する近似楕円１２０５を導出する。近似に用いるアルゴリズムは特に限定されず、最小二乗近似等の従来知られている方法を用いることができる。

　ステップＳ１１４０において、生成部１０２０は、導出された近似多角形（又は近似曲線）の大きさ、導出された近似多角形（又は近似曲線）の最大内接円の大きさ、又は導出された近似楕円の大きさを、定量的な情報として算出する。例えば、生成部１０２０は、導出された近似多角形の面積を算出することができる。また、生成部１０２０は、導出された近似多角形の最大内接円の半径、直径、又は面積を算出することができる。さらに、生成部１０２０は、導出された近似楕円の短軸長、長軸長、又は面積等を算出することができる。また、生成部１０２０は、導出された近似多角形又は近似楕円の周長を算出することもできる。生成部１０２０は、導出された近似多角形の最大内接円の中心、導出された近似楕円の長軸と短軸との交点、又は導出された近似多角形若しくは近似楕円の重心をさらに算出してもよい。ここで、近似楕円の短軸長は、分岐部分における分岐血管の直径に近似していると考えられる。また、近似多角形又は近似楕円の面積は、分岐面の面積に近似していると考えられる。

　こうして求められた分岐面情報は、分岐血管へと挿入するバルーンやステント等のデバイスを選択するための参考情報として用いられる。ステップＳ１１３０及びステップＳ１１４０において、生成部１０２０は、主血管と分岐血管との境界点群に内接する最大内接球を導出し、導出された最大内接球の大きさ及び位置を算出してもよい。こうして算出された値も分岐情報であり、分岐部分の大きさ及び位置を表している。

　ステップＳ１１５０～Ｓ１１７０において、生成部１０２０は、ステップＳ１１２０で検出された分岐部分についての、主血管からの分岐血管の分岐方向を算出する。

　ステップＳ１１５０において、生成部１０２０は、分岐部分における主血管の向きを算出する。例えば、生成部１０２０は、分岐部分の近傍での２以上の断層像における主血管の内腔の重心を算出し、重心の位置を近似するベクトルを導出することができる。この計算には、境界抽出部２００から取得した、主血管を示す交差点についての位置情報を使用することができる。このようにして導出されたベクトルは、主血管の向きを表している。重心の算出は、主血管を示す交差点についての位置情報を用いて行うことができる。主血管の向きを算出するために用いる断層像は、分岐部分が写っている断層像であってもよいし、分岐部分から所定の範囲内にある断層像であってもよい。もっとも、主血管の内腔の重心の代わりに、主血管の内腔の近似楕円の中心、又は主血管の内腔の最大内接円の中心等を用いることもできる。

　分岐部分の近傍に病変が存在する場合、病変部においては主血管の内腔が狭くなり、重心の位置も移動する。一実施形態においては、分岐部分の近傍に病変が存在する場合、病変よりも分岐部分に近い断層像が用いられる。この場合、主血管の向きを算出する際における病変の影響を抑えることができる。病変の検出方法は特に限定されないが、例えば主血管の内腔が閾値よりも小さい場合に、断層像には病変が写っていると判定することができる。別の実施形態において、生成部１０２０は、病変部が写っている断層像のうち、主血管の内腔の最大内接円が最も小さい断層像を特定する。そして、特定された断層像における内接円の中心を、主血管の内腔の重心の代わりに用いることができる。このような構成によれば、算出される分岐方向は、血管におけるカテーテル等の進行方向に対する分岐血管の進行方向に近似するため、治療計画を立てる上で有用となる。

　ステップＳ１１６０において、生成部１０２０は、分岐部分における分岐血管の分岐方向を算出する。本実施形態において、生成部１０２０は、分岐血管を示す交差点についての位置情報を参照して、分岐血管の分岐方向を算出する。例えば、生成部１０２０は、分岐血管に内接する内接球を導出することができる。そして、生成部１０２０は、ステップＳ１１４０で導出された近似多角形（又は近似曲線）又は近似楕円の重心から、導出された内接球の中心へと向かう方向に基づいて、分岐血管の分岐方向を算出することができる。

　具体的には、生成部１０２０は、分岐血管を示す交差点群に内接する最大内接球を導出する。ここで、最大内接球は、分岐血管を示す交差点群により囲まれる空間から、分岐部分側又は分岐部分から離れた開放部側（下流側）においてはみ出ないように導出される。また、分岐血管を示す交差点群に、分岐血管の始点を表す点群と、分岐部分から離れた開放部側（下流側）の点群と、を加えてから最大内接球を導出してもよい。そして、生成部１０２０は、ステップＳ１１４０で算出された近似多角形（又は近似曲線）又は近似楕円の重心から最大内接球の重心へと向かうベクトルを、分岐血管の分岐方向を示すベクトルとして算出する。もっとも、生成部１０２０は、近似多角形（又は近似曲線）又は近似楕円の重心の代わりに、ステップＳ１１４０で導出された近似多角形（又は近似曲線）の最大内接円の中心を用いてもよい。一実施形態においては、算出精度を高めるために、分岐血管を示す交差点群の中から、分岐部分からより離れている所定数の点群が選択され、選択された点群に内接する最大内接球が導出される。

　別の実施形態において、生成部１０２０は、分岐血管を示す交差点群に内接する２以上の内接球を導出し、一方の内接球の中心から、他方の内接球の中心へと向かう方向を、分岐血管の分岐方向として算出してもよい。この場合、分岐血管を示す交差点についての位置情報のみに基づいて分岐血管の分岐方向を算出することができる。一実施形態においては、分岐血管の分岐方向の算出精度を高めるために、一方の内接球が最大内接球であり、２つの内接球の中心間の距離が閾値以上であり、２つの内接球の半径の差が閾値以下となるように、２つの内接球が導出される。

　分岐血管に病変部が存在する場合には、主血管に病変が存在する場合と同様の処理を行うことができる。生成部１０２０は、こうして得られた分岐部分における分岐血管の分岐方向を、角度を用いて定量的に表すことができる。一実施形態においては、この分岐方向は、血管画像１９０を撮影した際の走査方向に対する角度で表すことができる。このように求められた分岐方向は、ガイドワイヤ等を分岐血管へと挿入する際の参考情報として用いることができる。

　一方で、本実施形態のステップＳ１１７０において、生成部１０２０は、より有用な情報を提供するために、主血管からの分岐血管の分岐角度を算出する。こうして算出された分岐角度は、主血管からの分岐部分における分岐血管の形態を示す定量的な情報に相当する。分岐角度の算出方法は特に限定されない。例えば、生成部１０２０は、主血管の向きを表すベクトルと、分岐血管の分岐方向を表すベクトルと、の間の角度を算出してもよい。

　別の実施形態においては、ユーザが分岐角度をより感覚的に理解しやすいように、図１３に示す角度Φ及びΘが算出される。図１３においては、主血管の向きを表すベクトル１３０１がｘｙ平面上にあり、ｘ軸と平行となるように、ｘｙｚ座標系（右手系のデカルト座標系）が設定されている。また、分岐血管の分岐方向を表すベクトル１３０２は、始点が原点に位置し、終点がｘｙ平面の第一象限に位置するように、このｘｙｚ座標系は設定されている。図１３において、投影点１３０３は、ベクトル１３０２の終点をｘｙ平面上に投影して得られる点である。

　こうして求められた分岐角度は、ガイドワイヤ等を分岐血管へと挿入する際の参考情報として用いることができる。

　ステップＳ１１８０において、生成部１０２０は、ステップＳ１１６０で算出された分岐方向と直交する断面における、分岐血管の大きさを示す情報を算出する。こうして算出された情報も、主血管からの分岐部分における分岐血管の形態を示す定量的な情報に相当する。分岐血管の大きさを示す情報としては、例えば、断面積、周長、又は最大内接円の半径若しくは直径等が挙げられる。生成部１０２０は、分岐血管を示す交差点群に従って分岐血管の大きさを示す情報を生成することができる。ここで、交差点群の輝度を投影することにより、分岐血管の断面画像を生成することもできる。投影される交差点群の輝度とは、血管画像１９０（例えば横断断面画像）における対応する画素の輝度値である。

　例えば、生成部１０２０は、分岐血管を示す交差点群を分岐方向と直交する断面へと投影してもよく、投影された点の位置に従って分岐血管の大きさを示す情報を算出することができる。この場合には、分岐血管を示す交差点群のうち、断面からの距離が所定範囲内にある点群を、断面へと投影することができる。

　別の実施形態においては、以下のように投影される点群が選択される。この選択方法を、図１４Ａ～Ｄを参照して説明する。分岐血管の分岐方向を表すベクトル１４０１が面上に来て、血管画像１９０の撮影時の撮影断面１４０２と直交する面をＸＹ面とする。Ｘ軸は、血管画像１９０の撮影時の走査方向を表す。また、ＸＹ面上において、走査方向Ｘと分岐方向を表すベクトルとがなす角をΩとする。ここで、血管画像１９０の撮影時のｉ番目の撮影断面を断面ｉと呼び、分岐方向を表すベクトル１４０１と直交する投影先の断面１４０３を断面ｊと呼ぶ。これらの関係は図１４Ａに示されている。

　まず、各断面ｉについて、分岐血管を示す交差点群のうち断面ｊに最も近い点を選択する。選択された点をＰ_ｉ，ｊと呼ぶ。次に、点Ｐ_ｉ，ｊとＰ_{ｉ＋１，ｊ}とをＸＹ面へと垂直に投影した際に、２つの投影点間を区切る垂直二等分面１４０４が算出される。図１４Ｂは、垂直二等分面１４０４の位置を示している。そして、垂直二等分面１４０４と、垂直二等分面１４０４に平行でＰ_ｉ，ｊを通る面と、の間にある断面ｉ上の分岐血管を示す交差点群が、断面ｊへと投影される点として選択される。また、垂直二等分面１４０４と、垂直二等分面１４０４に平行でＰ_{ｉ＋１，ｊ}を通る面と、の間にある断面ｉ＋１上の分岐血管を示す交差点群が、断面ｊへと投影される点として選択される。図１４Ｃには、選択された点１４０５が表されている。

　ただし、断面画像を複数作成する場合であって、断面ｊの間隔Ｄｓと、断面ｉの間隔Ｄｍとが、Ｄｓ＜Ｄｍ・ｃｏｓΩの関係を満たす場合には、分岐血管を示す交差点群のうち一部は、断面ｉと断面ｉ＋１との双方に投影される。このことを避けるために、以下の処理を行ってもよい。すなわち、Ｐ_ｉ，ｊとＰ_{ｉ，ｊ＋１}とをＸＹ面へと垂直に投影した際に、２つの投影点間を区切る垂直二等分面１４０６が算出される。そして、図１４Ｄに示すように、垂直二等分面１４０６と、垂直二等分面１４０６に平行でＰ_ｉ，ｊを通る面と、の間にある断面ｉ上の分岐血管を示す交差点群は、断面ｊへと投影されるが、断面ｊ＋１へは投影されない点として選択される。また、垂直二等分面１４０６と、垂直二等分面１４０６に平行でＰ_{ｉ，ｊ＋１}を通る面と、の間にある断面ｉ上の分岐血管を示す交差点群は、断面ｊ＋１へと投影されるが、断面ｊへは投影されない点として選択される。

　もっとも、血管画像１９０からの任意の断面における断層像の作成方法は多断面再構成法として公知である。したがって、生成部１０２０が、血管画像１９０を再構成することにより、分岐血管の分岐方向と直交する断面における断層像を作成する方法は、以上の方法に限られない。

　分岐方向と直交する断面における分岐血管の大きさを示す情報の算出方法の具体例としては、投影された点を近似する近似楕円、近似曲線及び近似多角形を導出してその面積や周囲の長さを求める方法、投影された点の最大内接円を導出してその半径、直径、若しくは面積を求める方法、等が挙げられる。

　上述したように、生成部１０２０は、分岐血管の分岐方向と直交する断面における断層像を作成し、出力することもできる。例えば、生成部１０２０は、ステップＳ１１４０で算出された近似多角形（若しくは近似曲線）又は近似楕円の重心を通る断層像を作成することができる。また、生成部１０２０は、分岐面のうち最も上流側又は最も下流側の位置を通る断層像を作成することもできる。生成部１０２０は、分岐血管の分岐方向に沿った視線方向、又は分岐方向に直交する視線方向における、分岐血管のサーフェスレンダリング画像又はボリュームレンダリング画像等を生成することもできる。

　このような構成によれば、バルーンやステント等を用いて主血管に対して治療を行った後に、プラークが分岐血管方向へと移動することによる分岐血管の狭窄の度合いを容易に確認することができる。また、以上の方法によれば、例えば分岐血管とガイドワイヤの影とが重なっている等の理由で、一部の断層像について分岐血管を示す交差点が十分に検出できない場合であっても、一部の交差点を用いて分岐血管の形態を示す定量的な情報が得られる可能性がある。

　分岐部分にステントが留置されている場合、ステップＳ１１４０において、生成部１０２０は、ステントストラット（ステントの支柱）の大きさ、例えば幅又は厚みを算出してもよい。ステントストラットは、断層像上には輝度の高い点として現れる。したがって、ガイドワイヤの像の推定位置と同様に、検出部２３０はステントストラットを形態学的にオブジェクトとして抽出することができ、その位置を検出することができる。この場合、生成部１０２０は、ステントストラットを表すオブジェクトを、分岐面へと垂直に投影することができる。そして、生成部１０２０は、分岐面をステントストラットにより分割することにより生じた分割領域の面積、又は最大内接円の大きさを算出することができる。こうして得られた情報により、主血管にステントを配置した後におけるステントストラットの幅又は厚みを確認することができる。

　また、ステップＳ１１８０において、生成部１０２０は、ステントストラットを表すオブジェクトを、分岐血管の分岐方向と直交する断面へと投影することができる。そして、生成部１０２０は、この断面上で、ステントストラットにより分割された分割領域の面積、周長、又は最大内接円の大きさを算出することができる。こうして得られた情報は、ステントストラットの間隙を通るようにカテーテルを操作する際における、カテーテルに垂直な方向のステントストラットの幅に対応する。したがって、この情報は、ステントストラットの間隙を通るようにガイドワイヤ等を操作する際、例えばキッシングバルーンテクニックを行うために用いるバルーンを選択するために、有用である。

　さらに、分岐部分にステントが配置されている場合、生成部１０２０は、ステントの大きさを示す情報を用いて、主血管からの分岐部分における分岐血管の形態を示す定量的な情報を補正することができる。ステントの大きさを示す情報は、情報取得部１０１０によって取得される。ステントの大きさを示す情報は、ユーザによって入力部（不図示）を介して入力されてもよいし、予め情報算出部１０００に記録されていてもよい。

　例えば、情報取得部１０１０は、ステントストラットの幅又は厚み（Ｉ１）を示す情報を取得することができる。この場合、生成部１０２０は、算出したステントストラットの幅又は厚み（Ｉ２）を用いて、補正パラメータ（Ｉ１／Ｉ２）を算出できる。算出されるステントストラットの幅又は厚みは、複数のステントストラットのそれぞれについて計測された値の平均値であってもよい。生成部１０２０は、こうして算出した補正パラメータ（Ｉ１／Ｉ２）を、分岐血管の形態を示す定量的な情報、例えばステップＳ１１４０で算出された近似楕円の短軸長等に乗算することにより、より正確な値を求めることができる。例えば、ステントストラットの幅を用いる場合、血管軸方向の長さを補正することができる。また、ステントストラットの厚みを用いる場合、走査ライン方向の長さを補正することができる。

　情報算出部１０００は、表示制御部１１０を介して、得られた分岐血管の形態を示す定量的な情報を表示部１２０に表示させることができる。例えば、表示制御部１１０は、定量的な情報を含む表示画面を表示部１２０に表示させることができる。また、表示制御部１１０は、血管の横断断面画像、軸方向断面画像、及び３次元画像のうちの少なくとも１つが含む表示画面を表示部１２０に表示させることもできる。定量的な情報と、血管の画像との双方を含む表示画面の例を図１５に示す。もっとも、定量的な情報と、血管の画像との双方を、表示画面が含むことは必須ではない。

　表示画面１５００には、血管軸方向断面画像１５１０が表示されており、血管軸方向断面画像１５１０には、検出された主血管かから分岐血管への分岐部分ＳＢ１～ＳＢ５が示されている。ユーザは、入力部（不図示）を介して、情報を知りたい分岐部分を指定することができる。表示制御部１１０は、複数の分岐部分のうちの１つを指定する入力に応じて、指定された分岐部分における分岐血管の形態を示す定量的な情報を含むように、表示画面を更新することができる。また、表示制御部１１０は、複数の分岐部分のうちの１つを指定する入力に応じて、指定された分岐部分を含む血管の横断断面画像、軸方向断面画像、及び３次元画像のうちの少なくとも１つを含むように、表示画面を更新することができる。

　例えば、表示画面１５００においては、分岐部分ＳＢ５が選択されている。図１５に示す表示画面１５００の領域１５２０には、分岐血管の分岐角度（Θ５，Φ５）、分岐血管の直径（Ｄ５）、及び分岐面の面積（Ｓ５）が表示されている。また、領域１５２０には、ステントストラットの幅又は厚みによって補正された後の、分岐血管の直径（Ｄ５）、及び分岐面の面積（Ｓ５）がさらに表示されている。

　また、表示画面１５００には、分岐部分ＳＢ５を含む血管軸方向断面画像の拡大画像１５３０，１５４０，１５５０と、分岐部分ＳＢ５を含む横断断面画像１５６０が表示されている。表示される拡大画像１５３０の断面位置は、画面奥行き方向へと変えることができる。また、拡大画像１５３０の断面の向きは、表示画面１５００の右下のユーザインタフェース１５８０を用いてコントロールできる。拡大画像１５４０は、図１３に示すｘｚ平面に平行で、中心（プローブ）を通る血管軸方向断面である。拡大画像１５５０には、検出されたステントストラットの位置が強調表示されている。表示画面１５００には、主血管についての血管軸方向断面画像と横断断面画像とが表示されているが、分岐部分ＳＢ５における分岐血管についての血管軸方向断面画像又は横断断面画像を表示してもよい。

　３Ｄ画像１５７０は、分岐部分ＳＢ５を含む血管の３次元画像が表示されている。一実施形態においては、表示される３次元画像の視線方向は、指定された分岐部分における分岐血管の分岐方向に応じて決定される。例えば、表示制御部１１０は、選択された分岐部分ＳＢ５における分岐血管の分岐方向のベクトルと一致するように、３Ｄ画像１５７０の視線方向を決定することができる。視線方向の決定方法はこの方法に限られず、例えば、表示制御部１１０は、選択された分岐部分ＳＢ５における分岐血管の分岐方向のベクトルと直交するように、３Ｄ画像１５７０の視線方向を決定してもよい。

［その他の実施形態］
　以上、境界抽出部２００及び情報算出部１０００を備える画像処理装置１００について説明した。しかしながら、境界抽出部２００を備える画像処理装置と、情報算出部１０００を備える画像処理装置とが、別個の装置であってもよい。この場合、境界抽出部２００を備える画像処理装置は、血管画像１９０を取得し、判定部２４０又は判定更新部２５０による判定結果を出力することができる。出力された情報に基づき、ユーザは血管画像１９０のうち分岐血管又はガイドワイヤの影に相当する部分を容易に識別することが可能となる。また、情報算出部１０００を備える画像処理装置は、境界抽出部２００が用いる方法とは別の方法で生成された、主血管に相当する部位と、分岐血管に相当する部位と、を区別する情報を取得してもよい。このような情報を用いても、情報算出部１０００を備える画像処理装置は、得られた分岐血管の形態を示す定量的な情報を生成及び出力することができる。

　図１に示した画像処理装置１００が備える各部の機能は、汎用のコンピュータを用いて実現することができる。上述のように、境界抽出部２００を備える画像処理装置と、情報算出部１０００を備える画像処理装置とが、別個の装置である場合も、それぞれの装置の機能は汎用のコンピュータを用いて実現することができる。

　図１６はコンピュータの基本構成を示す図である。図１６においてプロセッサ１６１０は、例えばＣＰＵであり、コンピュータ全体の動作をコントロールする。メモリ１６２０は、例えばＲＡＭであり、プログラム及びデータ等を一時的に記憶する。コンピュータが読み取り可能な記憶媒体１６３０は、例えばハードディスク又はＣＤ－ＲＯＭ等であり、プログラム及びデータ等を長期的に記憶する。本実施形態においては、記憶媒体１６３０が格納している、各部の機能を実現するプログラムが、メモリ１６２０へと読み出される。そして、プロセッサ１６１０が、メモリ１６２０上のプログラムを実行することにより、上述の各ステップの処理が行われ、各部の機能が実現される。

　図１６において、入力インタフェース１６４０は外部の装置から情報を取得するためのインタフェースである。また、出力インタフェース１６５０は外部の装置へと情報を出力するためのインタフェースであり、例えば表示部１２０等に接続されている。バス１６６０は、上述の各部を接続し、データのやりとりを可能とする。

　本発明は上記実施の形態に制限されるものではなく、本発明の精神及び範囲から離脱することなく、様々な変更及び変形が可能である。従って、本発明の範囲を公にするために、以下の請求項を添付する。

　本願は、２０１４年３月１４日提出の日本国特許出願特願２０１４－０５２３７１を基礎として優先権を主張するものであり、その記載内容の全てを、ここに援用する。

Claims

　プローブを用いて第１の管状物体内を走査することにより得られた管状物体画像のうち、前記第１の管状物体に相当する部位と、前記第１の管状物体から分岐する第２の管状物体に相当する部位と、を区別する情報を取得する情報取得手段と、
　前記情報を使用して、前記第１の管状物体からの分岐部分における前記第２の管状物体の形態を示す定量的な情報を生成する生成手段と、
　を備えることを特徴とする、画像処理装置。
　前記生成手段は、前記第１の管状物体からの分岐部分における前記第２の管状物体の大きさを、前記定量的な情報として算出することを特徴とする、請求項１に記載の画像処理装置。
　前記生成手段は、前記第１の管状物体と前記第２の管状物体との境界点群を近似する近似多角形又は近似楕円を導出し、前記近似多角形の大きさ、前記近似多角形の最大内接円の大きさ、又は前記近似楕円の大きさを、前記定量的な情報として算出することを特徴とする、請求項２に記載の画像処理装置。
　前記情報取得手段はさらに、前記管状物体画像のうち、前記第１の管状物体から前記第２の管状物体への分岐部分に留置されたステントの支柱に相当する部位を示す情報を取得し、
　前記生成手段は、前記ステントの支柱の間隙の大きさを算出することを特徴とする、請求項１乃至３の何れか１項に記載の画像処理装置。
　前記情報取得手段はさらに、前記ステントの支柱の間隙の大きさを示す情報を取得し、
　前記生成手段は、前記情報取得手段が取得した前記ステントの支柱の間隙の大きさを示す情報と、前記算出された前記ステントの支柱の間隙の大きさとに基づいて、前記定量的な情報を補正することを特徴とする、請求項４に記載の画像処理装置。
　前記生成手段は、前記第１の管状物体からの前記第２の管状物体の分岐方向を算出することを特徴とする、請求項１に記載の画像処理装置。
　前記生成手段は、
　　前記第１の管状物体と前記第２の管状物体との境界点群を近似する近似多角形又は近似楕円を導出し、
　　前記第２の管状物体に内接する内接球を導出し、
　　前記近似多角形又は近似楕円の重心から前記内接球の中心へと向かう方向に基づいて、前記分岐方向を算出する
　ことを特徴とする、請求項６に記載の画像処理装置。
　前記生成手段は、前記定量的な情報として、前記第１の管状物体からの前記第２の管状物体の分岐角度を前記分岐方向に基づいて算出することを特徴とする、請求項６又は７に記載の画像処理装置。
　前記生成手段は、前記定量的な情報として、前記分岐方向と直交する断面における前記第２の管状物体の大きさを示す情報を算出することを特徴とする、請求項６乃至８の何れか１項に記載の画像処理装置。
　前記生成手段は、前記分岐方向と直交する断面における前記第２の管状物体の断面画像を生成することを特徴とする、請求項６乃至９の何れか１項に記載の画像処理装置。
　前記情報取得手段はさらに、前記管状物体画像のうち、前記第１の管状物体から前記第２の管状物体への分岐部分に留置されたステントの支柱に相当する部位を示す情報を取得し、
　前記生成手段は、前記分岐方向と直交する断面へと前記ステントの支柱を垂直に投影した場合における、当該断面上での前記ステントの支柱の間隙の大きさを算出することを特徴とする、請求項６乃至１０の何れか１項に記載の画像処理装置。
　前記第１の管状物体及び前記第２の管状物体は血管であることを特徴とする、請求項１乃至１１の何れか１項に記載の画像処理装置。
　画像処理装置が行う画像処理方法であって、
　プローブを用いて第１の管状物体内を走査することにより得られた管状物体画像のうち、前記第１の管状物体に相当する部位と、前記第１の管状物体から分岐する第２の管状物体に相当する部位と、を区別する情報を取得する情報取得工程と、
　前記情報を使用して、前記第１の管状物体からの分岐部分における前記第２の管状物体の形態を示す定量的な情報を生成する生成工程と、
　を備えることを特徴とする、画像処理方法。
　請求項１３に記載の画像処理方法の各工程をコンピュータに実行させるためのプログラム。