WO2023189262A1

WO2023189262A1 - プログラム、情報処理方法および情報処理装置

Info

Publication number: WO2023189262A1
Application number: PCT/JP2023/008441
Authority: WO
Inventors: 聖清水; 耕一井上; 祐次郎松下
Original assignee: テルモ株式会社
Priority date: 2022-03-31
Filing date: 2023-03-07
Publication date: 2023-10-05

Abstract

撮像対象の大きさに関する情報を含む画像を生成できるプログラム等を提供すること。　プログラムは、走査面を軸方向に移動させながら画像を取得する画像取得用カテーテルを用いた三次元走査により生成された複数の断層像（４１）を取得し、それぞれの前記断層像（４１）について代表線（４２）を決定し、それぞれの前記断層像（４１）を構成するそれぞれの走査線と前記代表線（４２）との交点に基づいて決定した計算範囲における、前記走査線の指標値を算出し、前記断層像（４１）の位置と、前記交点の位置とにより定めた座標に、前記指標値に基づいて定めた色を表示する処理をコンピュータに実行させる。

Description

プログラム、情報処理方法および情報処理装置

　本発明は、プログラム、情報処理方法および情報処理装置に関する。

　三次元走査型ＩＶＵＳ（Intravascular Ultrasound）用カテーテルを使用して取得した一組の断層像に基づいて、それぞれの走査線の輝度パターンと、あらかじめ定義された分類パターンとの類似度を、走査線の角度とカテーテルの軸方向との２軸で配列した二次元画像を生成する診断支援装置が提案されている（特許文献１）。

特開２０２１－４１０２９号公報

　ユーザは、特許文献１の診断支援装置で生成された画像を観察することにより、たとえば血管に留置したステントの素線がどのように配置されているかを、血管を切り開いたような態様の画像で確認できる。

　しかしながら、特許文献１の装置で生成された画像は、走査線に含まれる撮像対象物の大きさの絶対量を評価できないという課題がある。たとえば、プラークの面積、ステントの拡張量、血管径の変化等は、特許文献１の装置で生成された画像においては観察できない。

　一つの側面では、撮像対象の大きさに関する情報を含む画像を生成できるプログラム等を提供することを目的とする。

　プログラムは、走査面を軸方向に移動させながら画像を取得する画像取得用カテーテルを用いた三次元走査により生成された複数の断層像を取得し、それぞれの前記断層像について代表線を決定し、それぞれの前記断層像を構成するそれぞれの走査線と前記代表線との交点に基づいて決定した計算範囲における、前記走査線の指標値を算出し、前記断層像の位置と、前記交点の位置とにより定めた座標に、前記指標値に基づいて定めた色を表示する処理をコンピュータに実行させる。

　一つの側面では、撮像対象の大きさに関する情報を含む画像を生成できるプログラム等を提供できる。

画像生成方法の概要を説明する説明図である。情報処理装置の構成を説明する説明図である。指標値ストライプの生成方法を説明する説明図である。指標値の例－１の算出方法を説明する説明図である。指標値の例－１の算出方法を説明する説明図である。指標値の例－１の算出方法を説明する説明図である。指標値の例－１の算出方法の変形例を説明する説明図である。指標値の例－２の算出方法を説明する説明図である。指標値の例－２の算出方法を説明する説明図である。指標値の例－３の算出方法を説明する説明図である。指標値の例－３の算出方法を説明する説明図である。指標値の例－３の算出方法を説明する説明図である。角度と長さとの変換方法を説明する説明図である。角度と長さとの変換方法を説明する説明図である。角度と長さとの変換方法を説明する説明図である。角度と長さとの変換方法を説明する説明図である。代表線の決定方法を説明する説明図である。代表線の決定方法を説明する説明図である。断層像ＤＢのレコードレイアウトを説明する説明図である。最大勾配点ＤＢのレコードレイアウトを説明する説明図である。代表線ＤＢのレコードレイアウトを説明する説明図である。指標値ＤＢのレコードレイアウトを説明する説明図である。プログラムの処理の流れを説明するフローチャートである。仮代表線算出のサブルーチンの処理の流れを説明するフローチャートである。代表線算出のサブルーチンの処理の流れを説明するフローチャートである。指標値算出のサブルーチンの処理の流れを説明するフローチャートである。第１指標値算出のサブルーチンの処理の流れを説明するフローチャートである。第２指標値算出のサブルーチンの処理の流れを説明するフローチャートである。第３指標値算出のサブルーチンの処理の流れを説明するフローチャートである。表示のサブルーチンの処理の流れを説明するフローチャートである。変形例の角度と面積との関係を説明する説明図である。変形例の角度と面積との関係を説明する説明図である。変形例の角度と面積との関係を説明する説明図である。実施の形態２の画面例である。実施の形態３のカテーテルシステムの構成を説明する説明図である。実施の形態４の情報処理装置の構成を説明する説明図である。

［実施の形態１］
　図１は、画像生成方法の概要を説明する説明図である。医師等のユーザは、画像取得用カテーテル２８（図２５参照）を患部付近に挿入する。画像取得用カテーテル２８は、三次元走査用であり、走査面を軸方向に少しずつ変更した複数の断層像４１を連続して撮影できる。

　断層像４１同士の間隔を、間隔Ｔで示す。以後の説明では、一回の三次元走査で作成された断層像４１を、一組の断層像４１と記載する場合がある。

　それぞれの断層像４１について、代表線４２が決定される。代表線４２は、たとえば断層像４１においてステントが留置された深さ、画像取得用カテーテル２８が挿入された血管の外弾性板、または、輝度が大きく変化する深さ等に基づいて定められた閉曲線である。代表線４２は、ユーザにより指定された閉曲線であってもよい。

　以下の説明では、代表線４２が画像取得用カテーテル２８を囲む円形である場合を例にして説明する。それぞれの断層像４１に対して、代表線４２の中心位置および半径が決定されている。代表線４２の決定方法の具体例については、後述する。

　断層像４１を形成するそれぞれの走査線４５（図３参照）、について、指標値が算出される。指標値は、走査線４５の特性および走査線４５上の交点４８（図３参照）の位置に基づいて算出される定数である。指標値の算出方法の具体例については後述する。

　走査線４５と代表線４２との交点４８に、それぞれの走査線４５について算出された指標値を割り当てて、指標値ストライプ４３を作成できる。指標値ストライプ４３は、指標値に基づいて着色した代表線４２を図１における下側の一か所、すなわち走査角が１８０度の位置で切断して作成して伸ばして形成される、仮想的な帯である。なお、交点４８の密度が低い場所に関しては、線形補間等の公知の補間手法を用いた補間処理が行なわれる。

　指標値ストライプ４３は、断層像４１同士の間隔Ｔと同じ幅を有し、代表線４２の長さＬと同じ長さを有する。図１において代表線４２および指標値ストライプ４３に接する短い矢印は、断層像４１において画像取得用カテーテル２８の真上、すなわち走査角が０度の位置を示す。指標値ストライプ４３を生成する方法の詳細については、後述する。

　それぞれの断層像４１に対応する指標値ストライプ４３を、たとえば走査角が０度の位置が一直線状に並ぶように配置することにより、指標値画像４４が生成される。指標値画像４４は、それぞれの断層像４１に対応する代表線４２により形成された立体の表面に指標値をマッピングし、走査角１８０度の位置で切り開いて平面に構成した画像である。

　ユーザは、指標値画像４４により代表線４２の大きさ、および代表線４２の存在する深さ周辺の指標値を直観的に把握し、診断および治療を行なえる。

　なお、代表線４２を切り開いて一本の指標値ストライプ４３を形成する位置は、走査角が１８０度の位置に限定しない。たとえば、患部等の注目領域が断層像４１の下側に描出されている場合には、１８０度の位置で代表線４２を切り開くことにより患部が指標値画像４４の上下に分かれてしまう。しかしながら、走査角が０度の位置で代表線４２を切り開いた状態の指標値ストライプ４３が、走査角が１８０度の位置が一直線状に並ぶように配置されることにより、注目領域を観察しやすい指標値画像４４が生成される。

　同様に指標値ストライプ４３を配置する基準の走査角度は、０度または１８０度の位置に限定しない。注目領域が存在する領域に対応する走査角度が一直線状に並ぶように代表線４２を配置することにより、注目する角度の領域が観察しやすい指標値画像４４が生成される。

　代表線４２を切り開く走査角度、および、指標値ストライプ４３を配置する際の基準とする走査角度は、ユーザが適宜選択できてもよい。

　図２は、情報処理装置２００の構成を説明する説明図である。情報処理装置２００は、制御部２０１、主記憶装置２０２、補助記憶装置２０３、通信部２０４、表示部２０５、入力部２０６およびバスを備える。制御部２０１は、本実施の形態のプログラムを実行する演算制御装置である。制御部２０１には、一または複数のＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）、またはマルチコアＣＰＵ等が使用される。制御部２０１は、バスを介して情報処理装置２００を構成するハードウェア各部と接続されている。

　主記憶装置２０２は、ＳＲＡＭ（Static Random Access Memory）、ＤＲＡＭ（Dynamic Random Access Memory）、フラッシュメモリ等の記憶装置である。主記憶装置２０２には、制御部２０１が行なう処理の途中で必要な情報および制御部２０１で実行中のプログラムが一時的に保存される。

　補助記憶装置２０３は、ＳＲＡＭ、フラッシュメモリ、ハードディスクまたは磁気テープ等の記憶装置である。補助記憶装置２０３には、断層像ＤＢ（Database）３６、最大勾配点ＤＢ３７、代表線ＤＢ３８、指標値ＤＢ３９、制御部２０１に実行させるプログラム、およびプログラムの実行に必要な各種データが保存される。断層像ＤＢ３６、最大勾配点ＤＢ３７、代表線ＤＢ３８および指標値ＤＢ３９は、情報処理装置２００に接続された外部の大容量記憶装置に保存されていてもよい。通信部２０４は、情報処理装置２００とネットワークとの間の通信を行なうインターフェースである。

　表示部２０５は、たとえば液晶表示パネルまたは有機ＥＬ（electro-luminescence）パネル等である。入力部２０６は、たとえばキーボードまたはマウス等である。表示部２０５と入力部２０６とは、積層されてタッチパネルを構成していてもよい。

　情報処理装置２００は、汎用のパソコン、タブレット、大型計算機、大型計算機上で動作する仮想マシン、または、量子コンピュータである。情報処理装置２００は、分散処理を行なう複数のパソコン、または大型計算機等のハードウェアにより構成されても良い。情報処理装置２００は、クラウドコンピューティングシステムにより構成されても良い。情報処理装置２００は、連携して動作する複数のパソコン、または大型計算機等のハードウェアにより構成されてもよい。

　図３は、指標値ストライプ４３の生成方法を説明する説明図である。図３を使用して、一枚の断層像４１に基づいて一本の指標値ストライプ４３を生成する方法を説明する。

　断層像４１は、ＸＹ形式断層像４１１とＲＴ形式断層像４１２の両方の形式で表現可能である。ＸＹ形式断層像４１１は、実際の形状に合わせて構築された画像である。ＸＹ形式断層像４１１の中心は、画像取得用カテーテル２８の中心に対応している。ＲＴ形式断層像４１２は、走査線４５を走査角度順に並行に並べて構築された画像である。ＲＴ形式断層像４１２の左端は、画像取得用カテーテル２８の中心に対応している。

　それぞれの走査線４５は、画像取得用カテーテル２８の中心からの距離と、当該距離における断層像４１の輝度とを関連づけた輝度データにより表される。画像取得用カテーテル２８を用いて、走査線４５ごとの輝度データを算出する方法、輝度データを用いて断層像４１を作成する方法、および、ＲＴ形式とＸＹ形式との間の変換方法は、いずれも公知であるため、説明を省略する。

　以後の説明においては、走査角度が－１８０度から１８０度までの範囲の走査線４５により一枚の断層像４１が形成されている場合を例にして説明する。ＸＹ形式断層像４１１においては、上方向が走査角度０度であり、時計回りが正の走査角度、反時計回りが負の走査角度であると定義して、説明を続ける。

　制御部２０１は、一枚の断層像４１を取得する。制御部２０１は断層像４１に基づいて、代表線４２を決定する。代表線４２は、ＸＹ形式断層像４１１においては画像取得用カテーテル２８を囲む円であり、ＲＴ形式断層像４１２においては画像の上端から下端まで伸びる曲線である。

　走査線４５は、ＸＹ形式断層像４１１においては画像の中心から放射状に延びる直線であり、ＲＴ形式断層像４１２においては水平線である。図３においては、走査角度９０度の走査線４５を例示する。ＸＹ形式断層像４１１とＲＴ形式断層像４１２のいずれにおいても一本の走査線４５は一か所で代表線４２と交差する。したがって、一本の走査線４５に対して一個の交点４８が存在する。

　制御部２０１は、それぞれの走査線４５について、後述する手順に基づいて指標値を算出する。指標値の例を、角度－指標値グラフ５４に図示する。角度－指標値グラフ５４の横軸はそれぞれの走査線４５の走査角度θであり、最小値は－１８０度、最大値は１８０度である。走査線４５の縦軸は指標値である。後述するように指標値には複数の定義が存在する。指標値の単位は、その定義により異なる。

　制御部２０１は、角度－指標値グラフ５４の横軸を代表線４２の弧の長さに変換する。変換した状態を模式的に説明する位置－指標値グラフ５５を図示する。位置－指標値グラフ５５の横軸は、走査角度０度の位置から走査角度θの位置までの代表線４２の弧の長さである。位置－指標値グラフ５５の縦軸は指標値である。変換方法の詳細については、後述する。

　制御部２０１は、位置－指標値グラフ５５の横軸の各位置における指標値の大小を色に割り当てて、指標値ストライプ４３を生成する。色は、白から黒までのモノトーンでも、カラーでもよい。色割り当ての具体例については後述する。

　以上の手順により、制御部２０１は、一組の断層像４１を構成するそれぞれの断層像４１に基づいて、一本の指標値ストライプ４３を生成する。制御部２０１は、複数の指標値ストライプ４３に基づいて、図１を使用して説明したように指標値画像４４を生成する。

［指標値の例－１］
　ガウス関数を使用して算出する指標値の例を説明する。図４Ａおよび図４Ｂは、指標値の例－１の算出方法を説明する説明図である。図４Ａは、一本の走査線４５に含まれる輝度データを示す輝度グラフ５９を示す。輝度グラフ５９の横軸は、画像取得用カテーテル２８の中心からの距離を示す。なお、以下の説明においては、画像取得用カテーテル２８の中心からの距離を、輝度グラフ５９上のデータ数で示す。すなわち、輝度グラフ５９の横軸の単位はデータの個数である。

　輝度グラフ５９の縦軸は、輝度を示す。縦軸の単位は０から２５５までの整数に正規化した輝度値である。以下の説明においては、画像取得用カテーテル２８の中心からｘ番目の輝度データを、Ｂ（ｘ）と記載する。

　横軸のＲ１は、輝度グラフ５９に表示した走査線４５が代表線４２と交差する位置を示す。ｂは定数である。定数ｂは、たとえば（１）式で定義される。

　Ｍag は、固定倍率である。
　Ｐeffは、有効画素数である。
　Ｄepは、画像取得深度である。

　制御部２０１は、ユーザによるＭagの変更指示を受け付けてもよい。なお、（１）式は定数ｂの定義の例示である。定数ｂは、ユーザが適宜指定した値であってもよい。なお、図４Ａにおいては、横軸が（Ｒ１－ｂ／２）未満の領域については、データ点の図示を省略している。

　図４Ｂは、平均値がＲ１、標準偏差がｂであるガウス関数グラフ５１を示す。ガウス関数グラフ５１の横軸は、図４Ａに示す輝度グラフ５９と同様に、画像取得用カテーテル２８の中心からの距離を示す。ガウス関数グラフ５１の縦軸は、ガウス関数である。画像取得用カテーテル２８の中心からｘ番目のデータに対応するガウス関数の値Ｎ（ｘ）は、（２）式により算出される。

　ガウス関数グラフ５１のピーク位置は、走査線４５と代表線４２との交点４８の位置に対応している。平均値Ｒ１および標準偏差ｂは、第１統計量の例示である。

　図５Ａは、指標値の例－１の算出方法を説明する説明図である。図５Ａは、図４Ａに示した輝度を、図４Ｂに示したガウス関数により重み付けした状態を示すグラフである。図５Ａの横軸は、図４Ａに示す輝度グラフ５９と同様に、画像取得用カテーテル２８の中心からの距離を示す。図５Ａの縦軸は、重み付け輝度である。画像取得用カテーテル２８の中心からｘ番目のデータに対応する重み付け輝度Ｂｗ（ｘ）は（３）式により算出される。
　　　Ｂｗ（ｘ）＝Ｎ（ｘ）・Ｂ（ｘ）　‥‥‥　（３）

　画像取得用カテーテル２８の中心からの距離が、（Ｒ１－ｂ／２）以上、画像取得範囲の最大値以下である計算範囲４７における重み付け輝度Ｂｗ（ｘ）の平均値が、指標値に使用される。平均値は、たとえば相加平均値である。平均値は、相乗平均値または調和平均値等であってもよい。制御部２０１は、それぞれの走査線４５に関する指標値を算出する。以後の説明においては、計算範囲４７における重み付け輝度Ｂｗ（ｘ）の平均値を第１指標値と記載する場合がある。

　計算範囲４７における重み付け輝度Ｂｗ（ｘ）の標準偏差または分散等の任意の統計量が指標値に使用されてもよい。統計量の倍数またはべき乗等の値が、第１指標値に使用されてもよい。

［変形例－１］
　図５Ｂは、指標値の例－１の算出方法の変形例を説明する説明図である。本変形例における計算範囲４７は、画像取得用カテーテル２８の中心からの距離が（Ｒ１－ｂ／２）以上、（Ｒ１＋ｂ）以下の範囲である。制御部２０１は、計算範囲４７における重み付け輝度Ｂｗ（ｘ）の平均値を算出して、第１指標値に使用する。

　本変形例によると、画像取得用カテーテル２８から離れた領域の影響を受けにくい第１指標値を算出して、使用できる。なお、画像取得用カテーテル２８がＩＶＵＳ用カテーテルである場合には、走査対象領域に石灰化病変やステントなどの音響インピーダンスの強い反射体が存在する場合に、外側に多重反射などのアーティファクトが生じる可能性がある。画像取得用カテーテル２８から離れた領域の影響を受けにくくすることにより、多重反射等のアーティファクトの影響を避けることができる。なお、計算範囲４７の範囲は、図５Ａおよび図５Ｂに限定しない。

　たとえば計算範囲４７は、画像取得用カテーテル２８の中心からの距離が（Ｒ１±ｂ／２）、（Ｒ１±ｂ／４）、または、（Ｒ１±ｂ）の範囲であってもよい。計算範囲４７は、画像取得用カテーテル２８の中心からの距離がＲ１以上（Ｒ１＋ｂ／２）以下、または、（Ｒ１－ｂ／２）以上Ｒ１以下の範囲であってもよい。

　以上に例示したように、走査線４５と代表線４２との交点４８付近を含むように計算範囲４７を定めることにより、断層像４１のうち代表線４２付近の情報を強く反映した指標値を算出できる。

［指標値の例－２］
　レイリー分布を使用して算出する指標値の例を説明する。指標値の例－１と共通する部分については、説明を省略する。

　図６Ａおよび図６Ｂは、指標値の例－２の算出方法を説明する説明図である。図６Ａは、一本の走査線４５の輝度を示す輝度グラフ５９を示す。輝度グラフ５９の横軸は、画像取得用カテーテル２８の中心からの距離を示す。輝度グラフ５９の縦軸は、輝度を示す。計算範囲４７は、画像取得用カテーテル２８の中心からの距離が（Ｒ１－ｂ／２）以上、（Ｒ１＋ｂ）以下の範囲である。

　図６Ｂは、計算範囲４７内の輝度データから推定したレイリー分布グラフ５２を示す。レイリー分布グラフ５２の横軸は、輝度グラフ５９と同様に画像取得用カテーテル２８の中心からの距離を示す。レイリー分布グラフ５２の縦軸は、レイリー分布の確率密度関数を示す。画像取得用カテーテル２８の中心からｘ番目のデータに対応するレイリー分布の確率密度関数ＰＬ（ｘ）は、（４）式により算出される。

　σは、レイリー分布のパラメータである。

　（４）式で示すレイリー分布の期待値は（５）式で、分散は（６）式である。

　パラメータσの最尤推定値σｇは、計算範囲４７内の輝度データを使用して（７）式により算出される。

　ｎは、計算範囲内の輝度データの数である。
　ｎ１は、計算範囲におけるｘの最小値である。
　ｎ２は、計算範囲におけるｘの最大値である。

　制御部２０１は、（８）式によりそれぞれの走査線４５に関する指標値を算出する。以後の説明においては、（８）式により算出される指標値を第２指標値と記載する場合がある。

　（６）式と（８）式とから明らかであるように、（８）式で示す第２指標値は、パラメータσの最尤推定値に基づいて算出した、レイリー分布の分散である。なお、第２指標値はレイリー分布の分散に限定しない。

　たとえば（８）式で算出される分散の平方根である標準偏差が第２指標に使用されてもよい。パラメータσの最尤推定値に基づいて算出した、（５）式で示すレイリー分布の期待値が第２指標値に使用されてもよい。（７）式に示すパラメータσの最尤推定値σｇが第２指標値に使用されてもよい。これらの値に基づいて算出された値が第２指標値に使用されてもよい。

［指標値の例－３］
　輝度ヒストグラム５３（図７Ｂ参照）を使用して算出する指標値の例を説明する。指標値の例－１と共通する部分については、説明を省略する。

　図７Ａおよび図７Ｂは、指標値の例－３の算出方法を説明する説明図である。図７Ａは、一本の走査線４５の輝度を示す輝度グラフ５９を示す。輝度グラフ５９の横軸は、画像取得用カテーテル２８の中心からの距離を示す。輝度グラフ５９の縦軸は、輝度を示す。計算範囲４７は、画像取得用カテーテル２８の中心からの距離が（Ｒ１－ｂ／２）以上、画像取得範囲の最大値以下の範囲である。

　図７Ｂは、計算範囲４７内の輝度データに基づいて算出した輝度ヒストグラム５３を示す。輝度ヒストグラム５３の横軸は輝度である。輝度ヒストグラム５３の縦軸は度数である。以下の説明においては、輝度ヒストグラム５３のうち、輝度が所定の第１閾値未満である部分を低輝度領域５３１と記載する。第１閾値は、たとえば６４である。

　制御部２０１は、計算範囲４７の輝度データから、輝度が第１閾値未満であるデータを抽出する。制御部２０１は、抽出したデータの平均値μＬおよび標準偏差σＬを算出する。制御部２０１は、（９）式に基づいて高輝度側の閾値である第２閾値を算出する。
　　　第２閾値＝μＬ＋３σＬ‥‥‥　（９）

　図８は、指標値の例－３の算出方法を説明する説明図である。図８は、図７Ｂと同様に、計算範囲４７内の輝度データに基づいて算出した輝度ヒストグラム５３を示す。図８においては、縦軸を０から１２の範囲に拡大して示す。図８の横軸に、（９）式により算出した第２閾値を示す。以下の説明においては、輝度ヒストグラム５３のうち、輝度が所定の第２閾値以上である部分を高輝度領域５３２と記載する。

　制御部２０１は、計算範囲４７の輝度データから高輝度領域５３２に含まれるデータを抽出する。制御部２０１は、抽出したデータの標準偏差σＨを算出する。標準偏差σＨは、第２統計量の例示である。以後の説明においては、標準偏差σＨを第３指標値と記載する場合がある。なお、第３指標値はσＨに限定しない。

　たとえば、高輝度領域５３２に含まれるデータの分散が、第３指標値に使用されてもよい。高輝度領域５３２に含まれるデータの、相加平均値、相乗平均値、調和平均値、中央値または最頻値等の任意の統計量が第３指標値に使用されてもよい。

　第２閾値の算出式は（９）式に限定しない。第２閾値は、所定の定数であってもよい。計算範囲４７は、画像取得用カテーテル２８の中心からの距離が（Ｒ１－ｂ／２）以上、画像取得範囲の最大値以下の範囲に限定しない。たとえば計算範囲４７は、画像取得用カテーテル２８の中心からの距離が（Ｒ１－ｂ／２）以上、（Ｒ１＋ｂ）、（Ｒ１±ｂ／２）、（Ｒ１±ｂ／４）、または、（Ｒ１±ｂ）の範囲であってもよい。計算範囲４７は、画像取得用カテーテル２８の中心からの距離がＲ１以上（Ｒ１＋ｂ／２）以下、または、（Ｒ１－ｂ／２）以上Ｒ１以下の範囲であってもよい。

　第１指標値、第２指標値、および第３指標値を算出する際の計算範囲４７は、共通の範囲であっても、互いに異なる範囲であっても良い。

　指標値は、センサ２８２（図２５参照）により取得した信号を輝度データに変換する途中段階のデータに基づいて算出されてもよい。指標値は、センサ２８２により取得した信号自体に基づいて算出されてもよい。

　それぞれの走査線４５の指標値には、たとえば、第１指標値、第２使用値または第３指標値のいずれか一つを選択して使用できる。このようにする場合、指標値ストライプ４３は、指標値を黒から白までのグレースケールに割り当てた白黒濃淡の帯であり、指標値画像４４は白黒の画像であるように定めることができる。

　第１指標値、第２使用値または第３指標値の三つの指標値が同時に使用されてもよい。たとえば制御部２０１は、第１指標値、第２使用値および第３指標値を、それぞれ０から２５５までの整数値に規格化する。制御部２０１は、たとえば規格化したそれぞれの指標値を、Ｒ（赤）、Ｇ（緑）、Ｂ（青）の輝度に重複しない任意の順番で割り当てる。指標値の番号と、色との組み合わせは任意である。

　このようにしてＲＧＢそれぞれの輝度を割り当てる場合、指標値ストライプ４３は指標値を２４ビットのいわゆるフルカラーで着色した帯であり、指標値画像４４はフルカラーの画像であるように定めることができる。

　なお、ＲＧＢは、表示部２０５に表示されるカラー画像を構成するカラーチャネルの例示である。制御部２０１は、ＲＧＢの代わりに、たとえばＨＬＳ色空間、または、ＨＳＶ（Hue, Saturation, Value／色相、明度、彩度）色空間等の、任意の色空間に対応するチャネルにそれぞれの指標値を割り当ててもよい。

　制御部２０１は、四個の指標値を算出し、ＲＧＢに加えてアルファチャンネルを加えたＲＧＢＡに割り当ててもよい。制御部２０１は、ＣＹＭＫ（Cyan Magenta Yellow Black）の４チャネルに指標値を割り当ててもよい。

　図９Ａから図９Ｃは、角度と長さとの変換方法を説明する説明図である。図９Ａは、ＸＹ形式断層像４１１に代表線４２を重畳した画像を模式的に示す。ＸＹ形式断層像４１１を構成する走査線４５のうち、第１走査線４５１と第２走査線４５２との二本の走査線４５を図９Ａに示す。ＸＹ形式断層像４１１において、画像取得用カテーテル２８の中心から上方向に延びる基準線と、代表線４２との交点４８を基準交点４８９で示す。基準交点４８９は、代表線４２上に定めた基準点の例示である。

　第１走査線４５１は、基準線から時計回りにθ１の位置にある。第１走査線４５１と代表線４２との交点４８を第１交点４８１で示す。代表線４２に沿って基準交点４８９から第１交点４８１までの長さをＬ１で示す。

　同様に第２走査線４５２は、基準線から反時計回りにθ２の位置にある。第２走査線４５２と代表線４２との交点４８を第２交点４８２で示す。代表線４２に沿って基準交点４８９から第２交点４８２までの長さをＬ２で示す。

　図９Ｂは、角度－指標値グラフ５４を示す。第１走査線４５１に対応する指標値をＡ１で、第２走査線４５２に対応する指標値をＡ２でそれぞれ示す。図９Ｃは、位置－指標値グラフ５５を示す。図９Ａにおいて、代表線４２の中心は画像取得用カテーテル２８の中心に対して右上に存在するため、位置－指標値グラフ５５において基準交点４８９に対応する横軸状の位置は、横軸の中央よりも左側に寄った位置に存在する。位置－指標値グラフ５５において、横軸の最小値と最大値との間の距離は、指標値ストライプ４３の長さＬである。

　図１０は、角度と長さとの変換方法を説明する説明図である。図１０を使用して、図９Ｂに示す角度－指標値グラフ５４を、図９Ｃに示す位置－指標値グラフ５５に変換する方法を説明する。

　図１０は、図９Ａの右上の部分を拡大した状態を示す。画像取得用カテーテル２８の中心と、第１交点４８１との間の距離をＲ１１で示す。破線で示す第３走査線４５３は、第１走査線４５１の隣の走査線４５である。第３走査線４５３と代表線４２との交点４８を第３交点４８３で示す。

　第１走査線４５１と第３走査線４５３との成す角度をΔθで示す。ＸＹ形式断層像４１１を構成するすべての走査線４５において、Δθは一定である。たとえば、一枚の断層像４１が５１２本の走査線４５で構成される場合、Δθは、約０．７度である。図１０においては、Δθを大きな角度で模式的に示す。代表線４２に沿って第１交点４８１から第３交点４８３までの長さをΔＬで示す。

　Δθは小さい角度であるため、ΔＬは、（１０）式で近似できる。
　　　ΔＬ＝Ｒ１１×ｓｉｎΔθ　‥‥‥　（１０）
　同様に、ｎ番目の走査線４５とｎ＋１番目の走査線４５との間の弧の長さΔＬｎは、（１１）式で近似できる。
　　　ΔＬｎ＝Ｒ１ｎ×ｓｉｎΔθ　‥‥‥　（１１）
　Ｒ１ｎは、ｎ番目の走査線４５と代表線４２との交点４８と、画像取得用カテーテル２８の中心との間の距離である。

　走査線４５を基準線から第１走査線４５１まで順次移動させながら、（１０）式で示すΔＬを加算することにより、代表線４２に沿った長さＬ１を算出できる。以上の処理により、制御部２０１は角度－指標値グラフ５４に基づいて位置－指標値グラフ５５を生成できる。位置－指標値グラフ５５で縦軸に示されている指標値を、前述のように色に変換することで、制御部２０１は一枚の断層像４１に対応する一本の指標値ストライプ４３を生成できる。

［代表線の決定－１］
　制御部２０１は、たとえばＸＹ形式断層像４１１上でユーザにより指定された三点を通る円を生成して、代表線４２に使用する。制御部２０１は、円を示すテンプレートをＸＹ形式断層像４１１に重畳表示し、ユーザによる位置および半径の変更を受け付け、ユーザが確定した円を代表線４２に使用してもよい。そのほか、任意のユーザインターフェイスを介して、制御部２０１はユーザによる代表線４２の指定を受け付けてもよい。

［代表線の決定－２］
　図１１および図１２は、代表線４２の決定方法を説明する説明図である。図１１および図１２を使用して、制御部２０１が自動的に代表線４２を決定する方法の例を説明する。

　制御部２０１は、ＲＴ形式断層像４１２をＲ軸と平行な線により、たとえば２５個のブロック４６に分割する。一枚のＲＴ形式断層像４１２が５１２本の走査線４５により構成されている場合、一つのブロック４６は２０本または２１本の走査線４５を含む。

　図１２に移動して、三次元走査により取得した複数の断層像４１の関係を説明する。ハッチングは、処理中の断層像４１を示す。制御部２０１は、処理中の断層像４１を、その前後Ｋ枚の断層像４１と組み合わせた断層像４１のグループを処理する。以下の説明では、Ｋが３であり、制御部２０１はｎ－３枚目の断層像４１からｎ＋３枚目の断層像４１までの合計７枚の断層像４１を組み合わせて処理する。

　なお、図１１においては、ブロック４６が八個である場合を模式的に示す。走査角度が－１８０度の側から順に第１ブロック４６１から第８ブロック４６８の符号をつけて説明する。制御部２０１は、それぞれのブロック４６ごとに、最大勾配点Ｇを算出する。最大勾配点Ｇは、Ｒ方向に輝度が大きくなる方向の変化、すなわち勾配が最大である位置である。

　座標ｚに対応するＲＴ形式断層像４１２上の各点における輝度の勾配は、（１２）式により算出できる。なお、（１２）式は、三次元走査により取得した複数のＲＴ形式断層像４１２を、走査角度θの方向に切断した平面に、Ｓｏｂｅｌ演算子Ｆを適用する畳み込み演算を意味する式である。

（ｒ，θ，ｚ）は、座標ｚに対応するＲＴ形式断層像における座標である。
　ｒは、画像取得用カテーテル２８の中心からの距離方向の位置である。
　θは、走査角方向の位置である。
　ｚは、三次元走査時の画像取得用カテーテル２８の軸方向における座標である。
　ｇ（ｒ，θ，ｚ）は、座標（ｒ，θ，ｚ）における輝度の勾配である。
　Ｉ（ｒ，θ，ｚ）は、（ｒ，θ，ｚ）座標ｚに対応するＲＴ形式断層像の輝度である。
　Ｆ（ｋ，ｊ）は、Ｓｏｂｅｌ演算子Ｆの（ｋ，ｊ）要素である。
　ｊ、ｋ、ｗは、畳み込み演算のパラメータである。

　（１２）式で使用するＳｏｂｅｌ演算子Ｆは、水平方向の輪郭検出に用いる演算子である。前述のように、７枚の断層像４１を組み合わせて処理する場合に使用するＳｏｂｅｌ演算子の一例を、（１３）式に示す。（１３）式のｊ行ｋ列の要素が、（１２）式のＦ（ｋ，ｊ）である。

　（１３）式に示す７×７のＳｏｂｅｌ演算子Ｆを使用する場合、（１２）式のｗは当該Ｓｏｂｅｌ演算子Ｆのカーネル半幅である３である。なお、Ｓｏｂｅｌ演算子Ｆは（１３）式に限定しない。水平方向の輪郭検出が可能な任意のＳｏｂｅｌ演算子Ｆを使用できる。前後Ｋ枚の断層像４１を組み合わせて処理する場合のＳｏｂｅｌ演算子Ｆは、縦（２Ｋ＋１）、横（２Ｋ＋１）の正方行列である。

　制御部２０１は、それぞれの座標ｚに対応するＲＴ形式断層像４１２について、ブロック４６ごとに（１２）式で算出したｇ（ｒ，θ，ｚ）が最大である座標（ｒ，θ，ｚ）を抽出する。抽出された座標（ｒ，θ，ｚ）が、ブロック４６の最大勾配点Ｇである。図１１に、第１ブロック４６１から第８ブロック４６８までのそれぞれのブロック４６で抽出されたＧ１からＧ８の最大勾配点Ｇを、ＲＴ形式断層像４１２上に模式的に示す。

　制御部２０１は、公知の座標変換により、ＲＴ形式断層像４１２をＸＹ形式断層像４１１に変換できる。図１１においては、ＸＹ形式断層像４１１上で最大勾配点Ｇ１から最大勾配点Ｇ８までが略円周形に配置されている。

　制御部２０１は、ＸＹ形式断層像４１１上で最大勾配点Ｇ１から最大勾配点Ｇ８の重心Ｃの座標を算出する。平面状の複数の点の座標から重心Ｃの座標を算出する方法は公知であるため、詳細については説明を省略する。制御部２０１は、算出した重心Ｃと、それぞれの最大勾配点Ｇとの間の距離Ｄを算出する。制御部２０１は、算出した距離Ｄの代表値である中央値Ｄｃを算出する。

　以上の処理により、中心が最大勾配点Ｇの重心Ｃであり、半径が重心Ｃとそれぞれの最大勾配点Ｇとの間の距離Ｄの中央値Ｄｃである仮代表線４２２が決定される。距離Ｄの代表値には、平均値が使用されてもよい。

　なお、断層像４１内のノイズおよびアーティファクト等により、抽出された最大勾配点Ｇの位置が不適切になる場合がある。不適切な最大勾配点Ｇが混在することによる悪影響を避けるため、制御部２０１は最大勾配点Ｇの外れ値を除去する処理を行なうことが望ましい。

　たとえば二枚目以降の断層像４１においては、制御部２０１は重心Ｃとそれぞれの最大勾配点Ｇとの間の距離Ｄを算出した後に、一枚前の断層像４１における距離Ｄの代表値Ｄｃとの差が所定の閾値を超える距離Ｄを除去する。その後、制御部２０１は残った距離Ｄの代表値Ｄｃを算出する。以上の処理により、制御部２０１は不適切な最大勾配点Ｇによる影響を受けにくい仮代表線４２２を決定できる。

　図１２に移動して、それぞれの断層像４１について仮代表線４２２を決定した後の処理を説明する。以後の説明においては、ｎ枚目の断層像４１における仮代表線４２２の中心をＣｎ、仮代表線４２２の半径をＤｎと記載する。

　制御部２０１は、ｎ枚目の断層像４１について、前後Ｋ枚の断層像４１における中心Ｃの座標および半径Ｄの移動平均値を算出する。たとえばＫが３である場合、制御部２０１はｎ－３枚目の断層像４１からｎ＋３枚目の断層像４１までの合計７枚の断層像４１の仮代表線４２２について、中心Ｃの座標の平均値Ｃｎavgと、半径Ｄの平均値Ｄｎavgとを算出する。制御部２０１は、算出したＣｎavgおよびＤｎavgを、ｎ枚目の断層像４１における代表線４２の中心および半径に使用する。

　以上の説明では、（１２）式および（１３）式を使用して説明した勾配の算出時と、ＣｎavgおよびＤｎavgの算出時とで、同一枚数の断層像４１を使用する場合を例にして説明したが、両者で異なる枚数の断層像４１を使用してもよい。たとえば制御部２０１は、Ｋ＝３で勾配を算出し、Ｋ＝７でＣｎavgおよびＤｎavgを算出してもよい。

　なお、一組の断層像４１の両端近傍に位置する断層像４１については、たとえば仮代表線４２２の中心および半径を代表線４２の中心および半径に使用する。一組の断層像４１の両端近傍に位置する断層像４１については、一番端の断層像４１までの断層像４１の枚数をＫの値に使用してもよい。

　以上の処理により、断層像４１内のノイズおよびアーティファクト等の影響を避け、一組の断層像４１により形成される三次元空間内で滑らかに繋がる立体形状を構成可能な代表線４２が作成される。

　代表線の決定－１、および、代表線の決定－２で説明した代表線４２の決定方法は例示であり、これらに限定しない。たとえば、パターンマッチング等の手法によりＸＹ形式断層像４１１上に代表線４２を決定してもよい。

　図１３Ａは、断層像ＤＢ３６のレコードレイアウトを説明する説明図である。断層像ＤＢ３６は、三次元走査により作成された断層像４１が記録されたデータベースである。断層像ＤＢ３６は、３Ｄ走査ＩＤフィールド、断層番号フィールドおよび断層像フィールドを有する。断層像フィールドは、ＸＹ形式フィールドおよびＲＴ形式フィールドを有する。

　３Ｄ走査ＩＤフィールドには、三次元走査ごとに付与される３Ｄ走査ＩＤが記録されている。断層番号フィールドには、一回の三次元走査で作成した断層像４１の順番を示す番号が記録されている。ＸＹ形式フィールドには、ＸＹ形式断層像４１１が記録されている。ＲＴ形式フィールドには、ＲＴ形式断層像４１２が記録されている。断層像ＤＢ３６は、一枚の断層像４１について一つのレコードを有する。

　なお、断層像ＤＢ３６にはＲＴ形式断層像４１２のみが記録されており、必要に応じて制御部２０１が座標変換によりＸＹ形式断層像４１１を作成してもよい。断層像４１が記録された断層像ＤＢ３６の代わりに、断層像４１が作成される前段階のたとえば音線データ等が記録されたデータベースに基づいて、制御部２０１が必要な断層像４１を随時作成してもよい。

　断層像ＤＢ３６は、断層番号フィールドを使用する代わりに、三次元走査の起点から断層像４１を取得した位置までの距離を記録するフィールドを有してもよい。制御部２１１は、たとえばＭＤＵ２８９（図２５参照）から三次元走査の起点からそれぞれの断層像４１を取得した位置までの距離を取得できる。

　図１３Ｂは、最大勾配点ＤＢ３７のレコードレイアウトを説明する説明図である。最大勾配点ＤＢ３７は、図１１を使用して説明した最大勾配点Ｇに関する情報が記録されたデータベースである。最大勾配点ＤＢ３７は、３Ｄ走査ＩＤフィールド、断層番号フィールド、距離代表値フィールド、ブロック番号フィールド、最大勾配位置フィールド、距離フィールドおよびフラグフィールドを有する。最大勾配位置フィールドは、ＲフィールドおよびＴフィールドを有する。

　３Ｄ走査ＩＤフィールドには、三次元走査ごとに付与される３Ｄ走査ＩＤが記録されている。断層番号フィールドには、一回の三次元走査で作成した断層像４１の順番を示す番号が記録されている。距離代表値フィールドには、図１１を使用して説明した重心Ｃとそれぞれの最大勾配点Ｇとの間の距離Ｄの代表値が記録されている。

　ブロック番号フィールドには、ブロック４６の番号が記録されている。Ｒフィールドにはブロック４６内における最大勾配点ＧのＲ座標、すなわち画像取得用カテーテル２８の中心からの距離が記録されている。Ｔフィールドにはブロック４６内における最大勾配点ＧのＴ座標、すなわち走査角度が記録されている。

　距離フィールドには、重心Ｃと最大勾配位置Ｇとの距離Ｄが記録されている。フラグフィールドには、最大勾配位置Ｇが外れ値であるか否かを示すフラグが記録されている。「１」は、最大勾配位置Ｇが外れ値ではなく、距離代表値フィールドに記録される代表値の算出に使用されることを意味する。「０」は、最大勾配位置Ｇが外れ値であり、距離代表値フィールドに記録される代表値の算出に使用されないことを意味する。最大勾配点ＤＢ３７は、一つのブロック４６について、一つのレコードを有する。

　図１３Ｃは、代表線ＤＢ３８のレコードレイアウトを説明する説明図である。代表線ＤＢ３８は、代表線４２に関する情報を記録するデータベースである。代表線ＤＢ３８は、３Ｄ走査ＩＤフィールド、断層番号フィールド、仮代表線フィールドおよび代表線フィールドを有する。仮代表線フィールドおよび代表線フィールドは、それぞれ重心フィールドおよび距離代表値フィールドを有する。仮代表線フィールドおよび代表線フィールドの重心フィールドは、それぞれＸフィールドおよびＹフィールドを有する。

　３Ｄ走査ＩＤフィールドには、三次元走査ごとに付与される３Ｄ走査ＩＤが記録されている。断層番号フィールドには、一回の三次元走査で作成した断層像４１の順番を示す番号が記録されている。仮代表線フィールドの重心フィールドには、図１１を使用して説明した最大勾配点Ｇの重心ＣのＸ座標およびＹ座標が記録されている。仮代表線フィールドの距離代表値フィールドには、距離Ｄの代表値Ｄｃが記録されている。

　代表線フィールドの重心フィールドには、図１２を使用して説明した重心Ｃの移動平均値である重心ＣｎavgのＸ座標およびＹ座標の移動平均値が記録されている。代表線フィールドの距離代表値フィールドには、図１２を使用して説明した距離代表値Ｄの移動平均である距離代表値Ｄｎavgが記録されている。代表線ＤＢ３８は、一枚の断層像４１について一つのレコードを有する。

　図１４は、指標値ＤＢ３９のレコードレイアウトを説明する説明図である。指標値ＤＢ３９は、それぞれの走査線４５に対応する指標値を記録したデータベースである。指標値ＤＢ３９は、３Ｄ走査ＩＤフィールド、断層番号フィールド、走査線番号フィールド、交点フィールドおよび指標値フィールドを有する。

　交点フィールドは、Ｒ１フィールドおよびＬフィールドを有する。指標値フィールドは、第１指標値フィールド、第２指標値フィールドおよび第３指標値フィールドを有する。なお、指標値フィールドはサブフィールドを有さなくてもよい。指標値フィールドは四個以上のサブフィールドを有してもよい。

　３Ｄ走査ＩＤフィールドには、三次元走査ごとに付与される３Ｄ走査ＩＤが記録されている。断層番号フィールドには、一回の三次元走査で作成した断層像４１の順番を示す番号が記録されている。走査線番号フィールドには、走査線の番号が記録されている。Ｒ１フィールドには、画像取得用カテーテル２８の中心と、交点４８との間の距離の距離Ｒ１が記録されている。Ｌフィールドには、代表線４２に沿って基準交点４８９から交点４８までの長さＬが記録されている。第１指標値フィールド、第２指標値フィールドおよび第３指標値フィールドには、それぞれ第１指標値、第２指標値および第３指標値が記録されている。

　図１５は、プログラムの処理の流れを説明するフローチャートである。制御部２０１は、断層像ＤＢ３６から一枚の断層像４１を取得する（ステップＳ５０１）。制御部２０１は、仮代表線算出のサブルーチンを起動する（ステップＳ５０２）。仮代表線算出のサブルーチンは、図１１を使用して説明した仮代表線４２２の中心および半径を算出するサブルーチンである。仮代表線算出のサブルーチンの処理の流れは後述する。

　制御部２０１は一組の断層像４１に含まれる断層像４１の処理を終了したか否かを判定する（ステップＳ５０４）。終了していないと判定した場合（ステップＳ５０４でＮＯ）、制御部２０１はステップＳ５０１に戻る。終了したと判定した場合（ステップＳ５０４でＹＥＳ）、制御部２０１は代表線算出のサブルーチンを起動する（ステップＳ５０５）。代表線算出のサブルーチンは、図１２を使用して説明したように仮代表線４２２の移動平均値に基づいて代表線４２の中心および半径を算出するサブルーチンである。代表線算出のサブルーチンの処理の流れは後述する。

　制御部２０１は、指標値を算出する断層像４１を選択する（ステップＳ５０６）。ステップＳ５０６で選択される断層像４１はステップＳ５０１からステップＳ５０４までのループで取得した断層像４１の中の一枚である。

　制御部２０１はステップＳ５０６で選択した断層像４１を構成する走査線４５のうちの一本を選択する（ステップＳ５０７）。なお、制御部２０１はステップＳ５０７を実行するたびに、図９Ａおよび図１０を使用して説明した基準交点４８９から時計回りに１本ずつ走査線４５を選択する。

　制御部２０１は、指標値算出のサブルーチンを起動する（ステップＳ５０８）。指標値算出のサブルーチンは、一本の走査線４５に基づいて指標値を算出するサブルーチンである。指標値算出のサブルーチンの処理の流れは後述する。

　制御部２０１は、指標値ＤＢ３９の、ステップＳ５０７で選択した走査線４５に対応するレコードを抽出する。制御部２０１は、Ｒ１フィールド、Ｌフィールドおよび指標値フィールドに、指標値算出のサブルーチンで算出した距離Ｒ１、長さＬおよび指標値をそれぞれ記録する（ステップＳ５１０）。

　制御部２０１は、ステップＳ５０６で取得した断層像４１を構成する走査線４５の処理を終了したか否かを判定する（ステップＳ５１１）。終了していないと判定した場合（ステップＳ５１１でＮＯ）、制御部２０１はステップＳ５０７に戻る。終了したと判定した場合（ステップＳ５１１でＹＥＳ）、制御部２０１は一組の断層像４１の処理を終了したか否かを判定する（ステップＳ５１２）。

　終了していないと判定した場合（ステップＳ５１２でＮＯ）、制御部２０１はステップＳ５０６に戻る。終了したと判定した場合（ステップＳ５１２でＹＥＳ）、制御部２０１は表示のサブルーチンを起動する（ステップＳ５１３）。表示のサブルーチンは、指標値ＤＢ３９に基づいて指標値画像４４を表示するサブルーチンである。表示のサブルーチンの処理の流れは後述する。その後、制御部２０１は処理を終了する。

　図１６は、仮代表線算出のサブルーチンの処理の流れを説明するフローチャートである。仮代表線算出のサブルーチンは、一枚の断層像４１に基づいて図１１を使用して説明した仮代表線４２２の中心および半径を算出するサブルーチンである。

　制御部２０１は、図１１を使用して説明したように、処理中の断層像４１を所定の数のブロック４６に分割する（ステップＳ５２１）。制御部２０１は、一つのブロック４６を選択する（ステップＳ５２２）。制御部２０１は、ブロック４６内の最大勾配点Ｇの座標を算出する（ステップＳ５２３）。具体的には、制御部２０１は（１２）式に基づいてブロック４６の各位置における輝度の勾配を算出する。制御部２０１は、算出した輝度の勾配が最大である位置の座標を抽出する。

　制御部２０１は、ステップＳ５２１で分割したブロック４６の処理を終了したか否かを判定する（ステップＳ５２４）。処理を終了していないと判定した場合（ステップＳ５２４でＮＯ）、制御部２０１はステップＳ５２２に戻る。

　処理を終了したと判定した場合（ステップＳ５２４でＹＥＳ）、制御部２０１は複数の最大勾配点Ｇの重心Ｃの座標を算出する（ステップＳ５２５）。制御部２０１はそれぞれの最大勾配点と重心Ｃとの間の距離Ｄを算出する（ステップＳ５２６）。

　制御部２０１は、処理中の断層像４１が一組の断層像４１のうちの一枚目であるか否かを判定する（ステップＳ５２７）。一枚目であると判定した場合（ステップＳ５２７でＹＥＳ）、制御部２０１はステップＳ５２６で算出した距離Ｄの代表値Ｄｃを算出する（ステップＳ５２８）。代表値Ｄｃは、たとえば中央値である。

　一枚目ではないと判定した場合（ステップＳ５２７でＮＯ）、制御部２０１はステップＳ５２６で算出した距離Ｄから外れ値を除外する（ステップＳ５３１）。外れ値は、たとえば一つ前の断層像４１で算出した代表値Ｄｃとの差が所定の閾値を超える距離Ｄである。制御部２０１は、距離Ｄの代表値Ｄｃを算出する（ステップＳ５３２）。

　ステップＳ５２８またはステップＳ５３２の終了後、制御部２０１はステップＳ５２５で算出した重心ＣおよびステップＳ５２８またはステップＳ５３２で算出した代表値Ｄｃを代表線ＤＢ３８に記録する（ステップＳ５３３）。

　具体的には制御部２０１は、代表線ＤＢ３８の、処理中の断層像４１に対応するレコードを抽出する。制御部２０１は、抽出したレコードの仮代表線フィールドの重心フィールドに、ステップＳ５２５で算出した重心ＣのＸ座標およびＹ座標をそれぞれ記録する。制御部２０１は、抽出したレコードの仮代表線フィールドの代表値フィールドに、ステップＳ５２８またはステップＳ５３２で算出した代表値Ｄｃを記録する。その後、制御部２０１は処理を終了する。

　図１７は、代表線算出のサブルーチンの処理の流れを説明するフローチャートである。制御部２０１は、代表線ＤＢ３８から一組の断層像４１分の仮代表線フィールドに記録された情報を取得する（ステップＳ５４１）。制御部２０１は、処理を行なう断層像４１の断層番号を選択する（ステップＳ５４２）。たとえば制御部２０１は、ステップＳ５４２を実行するたびに断層番号を１から順に１ずつ増加させる。

　制御部２０１は、ステップＳ５４２で選択した断層番号が一組の断層像４１の両端近傍を示すか否かを判定する（ステップＳ５４３）。具体的には、図１２を使用して説明したように、断層像４１の前後Ｋフレームでの移動平均値を使用する場合、ステップＳ５４２で選択した断層像番号がＫ以下である場合、および、（断層像４１の総数－Ｋ）以上である場合、制御部２０１は選択した断層像番号が両端近傍を示すと判定する。

　両端近傍を示すと判定した場合（ステップＳ５４３でＹＥＳ）、制御部２０１は仮代表線フィールドに記録された情報と同一の情報を、代表線フィールドに記録する（ステップＳ５４４）。両端近傍を示さないと判定した場合（ステップＳ５４３でＮＯ）、制御部２０１は処理中の断層番号の前後所定数の断層番号のデータを使用して、重心Ｃと、代表値Ｄｃの移動平均をそれぞれ算出する（ステップＳ５４５）。

　制御部２０１は、代表線ＤＢ３８の代表線フィールドに、算出した移動平均を記録する（ステップＳ５４６）。具体的には制御部２０１は、代表線ＤＢ３８の、処理中の断層像４１に対応するレコードを抽出する。制御部２０１は、抽出したレコードの代表線フィールドの重心フィールドに、ステップＳ５４５で算出した重心Ｃの移動平均を記録する。制御部２０１は、抽出したレコードの代表線フィールドの代表値フィールドに、ステップＳ５４５で算出した代表値Ｄｃの移動平均を記録する。

　制御部２０１は、ステップＳ５４１で取得した仮代表線４２２に関する処理を終了したか否かを判定する（ステップＳ５４７）。終了していないと判定した場合（ステップＳ５４７でＮＯ）、制御部２０１はステップＳ５４２に戻る。終了したと判定した場合（ステップＳ５４７でＹＥＳ）、制御部２０１は処理を終了する。

　図１８は、指標値算出のサブルーチンの処理の流れを説明するフローチャートである。指標値算出のサブルーチンは、一本の走査線４５に基づいて指標値を算出するサブルーチンである。

　制御部２０１は、処理中の走査線４５と、代表線４２との交点４８の座標を算出する（ステップＳ５５１）。ここで制御部２０１はＸＹ座標系と、ＲＴ座標系との両方で交点４８の座標を算出する。制御部２０１は、（１０）式により一つ前の走査線４５上の交点４８からの弧の長さであるΔＬを算出する（ステップＳ５５２）。

　制御部２０１は、過去に算出したΔＬの総和である基準交点４８９からの弧の長さＬを算出する（ステップＳ５５３）。なお、走査角度が１８０度以上の走査線４５については、制御部２０１は代表線４２の周長からΔＬの総和を減算して弧の長さＬを算出する。ステップＳ５５３により、図３を使用して説明した角度－指標値グラフ５４の横軸を変換して位置－指標値グラフ５５を算出する処理が実現される。

　制御部２０１は、処理中の走査線４５に関する走査線データを取得する（ステップＳ５５４）。制御部２０１は第１指標値算出のサブルーチンを起動する（ステップＳ５５５）。第１指標値算出のサブルーチンは、図４Ａ、図４Ｂおよび図５Ａを使用して説明した第１指標値を算出するサブルーチンである。第１指標値算出のサブルーチンの処理の流れは後述する。

　制御部２０１は第２指標値算出のサブルーチンを起動する（ステップＳ５５６）。第２指標値算出のサブルーチンは、図６Ａおよび図６Ｂを使用して説明した第２指標値を算出するサブルーチンである。第２指標値算出のサブルーチンの処理の流れは後述する。

　制御部２０１は第３指標値算出のサブルーチンを起動する（ステップＳ５５７）。第３指標値算出のサブルーチンは、図７Ａ、図７Ｂおよび図８を使用して説明した第３指標値を算出するサブルーチンである。第３指標値算出のサブルーチンの処理の流れは後述する。その後、制御部２０１は処理を終了する。

　図１９は、第１指標値算出のサブルーチンの処理の流れを説明するフローチャートである。第１指標値算出のサブルーチンは、図４Ａ、図４Ｂおよび図５Ａを使用して説明した第１指標値を算出するサブルーチンである。

　制御部２０１は、図４Ｂおよび（２）式を使用して説明したガウス関数を生成する（ステップＳ５６１）。制御部２０１は、（３）式を使用して説明した、ガウス関数により重み付けした重み付け輝度Ｂｗを算出する（ステップＳ５６２）。制御部２０１は、重み付け輝度Ｂｗから、図５Ａを使用して説明した計算範囲４７に含まれるデータを抽出する（ステップＳ５６３）。

　制御部２０１は、計算範囲４７に含まれるデータに基づいて第１指標値を算出する（ステップＳ５６４）。前述の通り、第１指標値はたとえば重み付け輝度Ｂｗの平均値である。その後、制御部２０１は処理を終了する。

　図２０は、第２指標値算出のサブルーチンの処理の流れを説明するフローチャートである。第２指標値算出のサブルーチンは、図６Ａおよび図６Ｂを使用して説明した第２指標値を算出するサブルーチンである。

　制御部２０１は、輝度データから図６Ａを使用して説明した計算範囲４７に含まれるデータを抽出する（ステップＳ５７１）。制御部２０１は、（８）式に基づいて第２指標値を算出する（ステップＳ５７２）。その後、制御部２０１は処理を終了する。

　図２１は、第３指標値算出のサブルーチンの処理の流れを説明するフローチャートである。第３指標値算出のサブルーチンは、図７Ａ、図７Ｂおよび図８を使用して説明した第３指標値を算出するサブルーチンである。

　制御部２０１は、輝度データから図７Ａを使用して説明した計算範囲４７に含まれるデータを抽出する（ステップＳ５８１）。制御部２０１は抽出したデータの度数分布を算出する（ステップＳ５８２）。多数のデータに基づいて度数分布を算出する処理は公知であるため、詳細については説明を省略する。

　制御部２０１は、図７Ｂを使用して説明した低輝度領域５３１のデータを抽出する。制御部２０１は、抽出したデータの平均値μＬおよび標準偏差σＬの２つの統計量を算出する（ステップＳ５８３）。制御部２０１は、（９）式に基づいて高輝度側の閾値である第２閾値を算出する（ステップＳ５８４）。

　制御部２０１は、輝度データから図８を使用して説明した高輝度領域５３２のデータを抽出する。制御部２０１は、抽出したデータの標準偏差σＨである第３指標値を算出する（ステップＳ５８５）。その後、制御部２０１は処理を終了する。

　図２２は、表示のサブルーチンの処理の流れを説明するフローチャートである。表示のサブルーチンは、指標値ＤＢ３９に基づいて指標値画像４４を表示するサブルーチンである。

　制御部２０１は、指標値ＤＢ３９の第１指標値フィールドから、一組の断層像４１を構成する走査線４５に関する第１指標値を取得する。制御部２０１は、取得した第１指標値を０から２５５の整数に正規化する（ステップＳ５９１）。

　制御部２０１は、指標値ＤＢ３９の第２指標値フィールドから、一組の断層像４１を構成する走査線４５に関する第２指標値を取得する。制御部２０１は、取得した第２指標値を０から２５５の整数に正規化する（ステップＳ５９２）。

　制御部２０１は、指標値ＤＢ３９の第３指標値フィールドから、一組の断層像４１を構成する走査線４５に関する第３指標値を取得する。制御部２０１は、取得した第３指標値を０から２５５の整数に正規化する（ステップＳ５９３）。

　制御部２０１は、それぞれの走査線４５について、正規化した第１指標値から第３指標値までをカラーコードに変換する（ステップＳ５９４）。具体的には、第１指標値、第２指標値、第３指標値のうち、Ｒ（赤）に割り当てる指標値をＡＲ、Ｇ（緑）に割り当てる指標値ＡＧ、Ｂ（青）に割り当てる指標値をＡＢで表す場合、カラーコードはＲＧＢ（ＡＲ，ＡＧ、ＡＢ）と表現できる。

　制御部２０１は、それぞれの走査線４５について、指標値ＤＢ３９のＬフィールドに記録された長さＬと、断層番号フィールドに記録された断層番号に間隔Ｔを積算した断層位置Ｚとで定まる座標位置に、カラーコードにより表現される色を設定する。制御部２０１は、走査線４５同士の隙間に相当する座標位置のカラーコードを補間により決定する（ステップＳ５９５）。

　制御部２０１は、それぞれの座標位置のカラーコードにより生成される画像を表示部２０５に表示する（ステップＳ５９６）。制御部２０１はユーザによる指示に応じて、画像の回転、拡大、縮小等を行なう。制御部２０１は、ユーザが画像を操作する際に使用可能なインジケータおよびハンドル等を、表示部２０５に表示してもよい。その後、制御部２０１は処理を終了する。

　本実施の形態によると、走査線４５の大きさに関する情報を含む画像を生成できる。走査線４５のうち、代表線４２付近の部分の情報に重点を置いた指標を展開図のような態様で表示する指標値画像４４により、医師等のユーザが注目したい部分を直観的に把握できる。

　たとえば代表線４２を血管内腔に対応する位置に設定した場合、それぞれの代表線４２の長さが、図１に示す指標値画像４４の縦方向の長さで表現されているため、ユーザは、注目したい領域が狭窄しているのか、拡張しているのかを容易に把握できる。

　算出方法が異なる３つの指標を、Ｒ，Ｇ、Ｂに割り当てて１枚の指標値画像４４に表示することにより、ユーザは多くの情報を短時間に把握できる。なお、使用する指標の種類、および指標と色との組み合わせは、ユーザが適宜変更できても良い。

　制御部２０１は、指標値画像４４と断層像４１とを並べて表示部２０５に表示してもよい。ユーザが指標値画像４４上の一点をクリックした場合に、制御部２０１はクリックされた位置を含む断層像４１を表示部２０５に表示してもよい。

［変形例－２］
　制御部２０１は、図１０を使用して説明した、断層像４１に沿った長さＬおよびΔＬを、代表線４２の中心点を基準にして定めてもよい。具体的には、制御部２０１は図１０を使用して説明したように第１交点４８１、第２交点４８２等、それぞれの走査線４５と代表線４２との交点について、指標値を算出する。

　制御部２０１は、代表線４２の中心点を算出する。制御部２０１は、中心点から一定の角度ごとにΔＬを定めて、指標値ストライプ４３を生成する。すなわち、代表線４２全体にわたって、ΔＬは一定の値である。

　本変形例によると、狭窄等により血管のサイズ変化を判りやすく表示するカテーテルシステム１０を提供できる。

［変形例－３］
　本変形例においては、代表線４２の弧に沿った長さＬの代わりに面積を使用する。図２３Ａから図２３Ｃは、変形例の角度と面積との関係を説明する説明図である。

　図２３Ａは、ＸＹ形式断層像４１１に代表線４２を重畳した画像を模式的に示す。本変形例においては、基準交点４８９から第１交点４８１までの長さＬ１の代わりに、代表線４２と、第１走査線４５１とＹ軸とにより囲まれた略扇形の面積Ｓ１を使用する。同様に、基準交点４８９から第２交点４８２までの長さＬ２の代わりに、代表線４２と、第２走査線４５２と、Ｙ軸とにより囲まれた略扇形の面積Ｓ２を使用する。

　図２３Ｂは、図９Ｂと同様の角度－指標値グラフ５４を示す。図２３Ｃは、面積－指標値グラフ５６を示す。図２３の横軸は、図２３Ａを使用して説明した面積である。図２３の縦軸は指標値である。

　図２３Ａから明らかであるように走査角度と面積とは一対一対応している。面積Ｓ１は、たとえば図２３Ａにおいて左下がりのハッチングを行なった部分に含まれる画素数と、１画素の面積との積により算出できる。したがって制御部２０１は、角度－指標値グラフ５４の横軸を面積に変換して、面積－指標値グラフ５６を作成できる。

　面積－指標値グラフ５６に基づいて指標値ストライプ４３を生成することにより、面積Ｓと断層位置Ｚとで定まる座標位置に、それぞれの走査線４５に対応する指標値を設定した指標値画像４４を生成できる。

［変形例－４］
　本変形例においては、それぞれの断層像４１について同一の長さの指標値ストライプ４３を生成した後に、代表線４２に関する値に基づいて全体的に伸縮させる。

　具体的には、図３を使用して説明した生成方法において、制御部２０１は位置－指標値グラフ５５を作成せず、角度－指標値グラフ５４に基づいて指標値ストライプ４３を生成する。ここで生成する指標値ストライプ４３の長さは、３６０度の角度に対応しており、すべての断層像４１に対して同一の長さの指標値ストライプ４３が生成される。

　制御部２０１は、それぞれの断層像４１について、代表線４２に関する値を算出する。代表線４２に関する値は、たとえば代表線４２の長さ、代表線４２の内側の面積、代表線４２の平均径等である。制御部２０１は、指標値ストライプ４３を代表線４２に関する値に基づいて伸縮させる。以上の手順により、制御部２０１は、図１を使用して説明したように断層像４１ごとに異なる長さの指標値ストライプ４３を生成する。

［変形例－５］
　制御部２０１は、代表線４２に関する値に基づいて、指標値画像４４に指標を表示してもよい。たとえば制御部２０１は、代表線４２の平均径が最小である断層像４１の位置に対応する指標を、指標値画像４４の縁に表示する。制御部２０１は、代表線４２の平均径が採用である断層像４１に対応する指標値ストライプ４３を、他の指標値ストライプ４３とは異なる色で着色して表示してもよい。

［変形例－６］
　制御部２０１は、複数の代表線４２を作成して、指標値を算出してもよい。たとえば制御部２０１は、血管内腔に対応する代表線４２と、外弾性板に対応する代表線４２との２本の代表線４２を作成する。制御部２０１は、外弾性板に対応する代表線４２の面積から、血管内腔に対応する代表線４２の面積を減算する。以上により制御部２０１は、それぞれの断層像４１についてプラーク面積に対応する指標値を算出できる。

［実施の形態２］
　図２４は、実施の形態２の画面例である。制御部２０１は、長さＬと、断層位置Ｚに第３の軸Ｕを加えた三次元空間により定まる座標位置に、指標値に対応する色を設定する。指標値画像４４は、三次元空間内の曲面の態様で表示される。

　第３の軸Ｕは、たとえば図９ＡにおけるＸＹ形式断層像４１１の中心、すなわち画像取得用カテーテル２８の中心から、交点４８までの距離である。そのほか任意のパラメータを第３の軸Ｕに設定してもよい。

　本変形例によると、さらに多数の情報を同時に表現した指標値画像４４を表示する情報処理装置２００を提供できる。

［実施の形態３］
　本実施の形態はリアルタイムで断層像４１を取得して、指標値画像４４を表示するカテーテルシステム１０に関する。実施の形態１と共通する部分については、説明を省略する。

　図２５は、実施の形態３のカテーテルシステム１０の構成を説明する説明図である。カテーテルシステム１０は、画像処理装置２１０と、カテーテル制御装置２７とＭＤＵ（Motor Driving Unit）２８９と、画像取得用カテーテル２８とを備える。画像取得用カテーテル２８は、ＭＤＵ２８９およびカテーテル制御装置２７を介して画像処理装置２１０に接続されている。

　画像処理装置２１０は、制御部２１１、主記憶装置２１２、補助記憶装置２１３、通信部２１４、表示部２１５、入力部２１６およびバスを備える。制御部２１１は、本実施の形態のプログラムを実行する演算制御装置である。制御部２１１には、一または複数のＣＰＵ、ＧＰＵ、またはマルチコアＣＰＵ等が使用される。制御部２１１は、バスを介して画像処理装置２１０を構成するハードウェア各部と接続されている。

　主記憶装置２１２は、ＳＲＡＭ、ＤＲＡＭ、フラッシュメモリ等の記憶装置である。主記憶装置２１２には、制御部２１１が行なう処理の途中で必要な情報、および、制御部２１１で実行中のプログラムが一時的に保存される。

　補助記憶装置２１３は、ＳＲＡＭ、フラッシュメモリ、ハードディスクまたは磁気テープ等の記憶装置である。補助記憶装置２１３には、断層像ＤＢ３６、最大勾配点ＤＢ３７、代表線ＤＢ３８、指標値ＤＢ３９、制御部２１１に実行させるプログラム、およびプログラムの実行に必要な各種データが保存される。断層像ＤＢ３６、最大勾配点ＤＢ３７、代表線ＤＢ３８および指標値ＤＢ３９、画像処理装置２１０に接続された外部の大容量記憶装置等に記憶されていてもよい。外部の大容量記憶装置は、いわゆるクラウドストレージであってもよい。

　通信部２１４は、画像処理装置２１０とネットワークとの間の通信を行なうインターフェースである。表示部２１５は、たとえば液晶表示パネルまたは有機ＥＬパネル等である。入力部２１６は、たとえばキーボードおよびマウス等である。表示部２１５に入力部２１６が積層されてタッチパネルを構成していてもよい。表示部２１５は、画像処理装置２１０に接続された表示装置であってもよい。

　画像処理装置２１０は、汎用のパソコン、タブレット、大型計算機、または、大型計算機上で動作する仮想マシンである。画像処理装置２１０は、分散処理を行なう複数のパソコン、または大型計算機等のハードウェアにより構成されても良い。画像処理装置２１０は、クラウドコンピューティングシステムにより構成されても良い。画像処理装置２１０とカテーテル制御装置２７とは、一体のハードウェアを構成していてもよい。

　画像取得用カテーテル２８は、シース２８１と、シース２８１の内部に挿通されたシャフト２８３と、シャフト２８３の先端に配置されたセンサ２８２とを有する。ＭＤＵ２８９は、シース２８１の内部でシャフト２８３およびセンサ２８２を回転および進退させる。

　センサ２８２は、たとえば超音波の送受信を行なう超音波トランスデューサ、または、近赤外光の照射と反射光の受光とを行なうＯＣＴ（Optical Coherence Tomography）用の送受信部である。以下の説明では、画像取得用カテーテル２８は循環器の内側から超音波断層像を撮影する際に用いられるＩＶＵＳ（Intravascular Ultrasound）用カテーテルである場合を例にして説明する。

　カテーテル制御装置２７は、センサ２８２の一回転ごとに１枚の断層像４１を作成する。ＭＤＵ２８９がセンサ２８２を軸方向に引っ張りながら、または軸方向に押し込みながら回転させる操作により、カテーテル制御装置２７はシース２８１に略垂直な複数枚の断層像４１を連続的に作成する。

　制御部２１１は、カテーテル制御装置２７から断層像４１を逐次取得して、断層像ＤＢ３６に記録する。以上により、いわゆる三次元走査が行なわれ、一組の断層像４１が断層像ＤＢ３６に記録される。

　センサ２８２の進退操作には、画像取得用カテーテル２８全体を進退させる操作と、シース２８１の内部でセンサ２８２を進退させる操作との両方を含む。進退操作は、ＭＤＵ２８９により所定の速度で自動的に行なわれても、ユーザにより手動で行なわれても良い。

　なお、画像取得用カテーテル２８は機械的に回転および進退を行なう機械走査方式に限定しない。たとえば、複数の超音波トランスデューサを環状に配置したセンサ２８２を用いた、電子ラジアル走査型の画像取得用カテーテル２８であってもよい。

　制御部２１１は、以上の動作により記録した断層像ＤＢ３６を用いて、実施の形態１で説明した処理を行ない、表示部２１５に指標値画像４４を表示する。

［実施の形態４］
　図２６は、実施の形態４の情報処理装置２００の構成を説明する説明図である。本実施の形態は、汎用のコンピュータ９０と、プログラム９７とを組み合わせて動作させることにより、本実施の形態の情報処理装置２００を実現する形態に関する。実施の形態１と共通する部分については、説明を省略する。

　コンピュータ９０は、前述の制御部２０１、主記憶装置２０２、補助記憶装置２０３、通信部２０４、表示部２０５、入力部２０６およびバスに加えて読取部２０９を備える。

　プログラム９７は、可搬型記録媒体９６に記録されている。制御部２０１は、読取部２０９を介してプログラム９７を読み込み、補助記憶装置２０３に保存する。また制御部２０１は、コンピュータ９０内に実装されたフラッシュメモリ等の半導体メモリ９８に記憶されたプログラム９７を読出してもよい。さらに、制御部２０１は、通信部２０４および図示しないネットワークを介して接続される図示しない他のサーバコンピュータからプログラム９７をダウンロードして補助記憶装置２０３に保存してもよい。

　プログラム９７は、コンピュータ９０の制御プログラムとしてインストールされ、主記憶装置２０２にロードして実行される。以上により、実施の形態１で説明した情報処理装置２００が実現される。本実施の形態のプログラム９７は、プログラム製品の例示である。

　各実施例および変形例で記載されている技術的特徴（構成要件）はお互いに組合せ可能であり、組み合わせすることにより、新しい技術的特徴を形成することができる。
　今回開示された実施の形態はすべての点で例示であって、制限的なものでは無いと考えられるべきである。本発明の範囲は、上記した意味では無く、請求の範囲によって示され、請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。

　１０　　カテーテルシステム
　２００　情報処理装置
　２０１　制御部
　２０２　主記憶装置
　２０３　補助記憶装置
　２０４　通信部
　２０５　表示部
　２０６　入力部
　２０９　読取部
　２１０　画像処理装置
　２１１　制御部
　２１２　主記憶装置
　２１３　補助記憶装置
　２１４　通信部
　２１５　表示部
　２１６　入力部
　２７　　カテーテル制御装置
　２８　　画像取得用カテーテル
　２８１　シース
　２８２　センサ
　２８３　シャフト
　２８９　ＭＤＵ
　３６　　断層像ＤＢ
　３７　　最大勾配点ＤＢ
　３８　　代表線ＤＢ
　３９　　指標値ＤＢ
　４１　　断層像
　４１１　ＸＹ形式断層像
　４１２　ＲＴ形式断層像
　４２　　代表線
　４２２　仮代表線
　４３　　指標値ストライプ
　４４　　指標値画像
　４５　　走査線
　４５１　第１走査線
　４５２　第２走査線
　４５３　第３走査線
　４６　　ブロック
　４６１　第１ブロック
　４６８　第８ブロック
　４７　　計算範囲
　４８　　交点
　４８１　第１交点
　４８２　第２交点
　４８９　基準交点
　５１　　ガウス関数グラフ
　５２　　レイリー分布グラフ
　５３　　輝度ヒストグラム
　５３１　低輝度領域
　５３２　高輝度領域
　５４　　角度－指標値グラフ
　５５　　位置－指標値グラフ
　５６　　面積－指標値グラフ
　５９　　輝度グラフ
　９０　　コンピュータ
　９６　　可搬型記録媒体
　９７　　プログラム
　９８　　半導体メモリ

Claims

　走査面を軸方向に移動させながら画像を取得する画像取得用カテーテルを用いた三次元走査により生成された複数の断層像を取得し、
　それぞれの前記断層像について代表線を決定し、
　それぞれの前記断層像を構成するそれぞれの走査線と前記代表線との交点に基づいて決定した計算範囲における、前記走査線の指標値を算出し、
　前記断層像の位置と、前記交点の位置とにより定めた座標に、前記指標値に基づいて定めた色を表示する
　処理をコンピュータに実行させるプログラム。
　前記代表線は、前記画像取得用カテーテルを囲む閉曲線である
　請求項１に記載のプログラム。
　前記代表線は、前記画像取得用カテーテルを囲む円形である
　請求項１または請求項２に記載のプログラム。
　前記座標は、前記断層像の位置と、前記代表線に関する値に基づいて定められる
　請求項１から請求項３のいずれか一つに記載のプログラム。
　前記代表線に関する値は、前記代表線上に定めた基準点から前記交点までの前記代表線に沿った長さである
　請求項４か一つに記載のプログラム。
　前記指標値は、互いに異なる算出方法により算出される第１指標値、第２指標値および第３指標値のうち１つ以上に基づいて定められる
　請求項１から請求項５のいずれか一つに記載のプログラム。
　前記第１指標値は、
　　前記走査線における輝度値を、前記交点をピーク位置とするガウス関数により重み付けし、
　　重み付けした輝度値を、前記計算範囲において平均して算出された値である
　請求項６に記載のプログラム。
　前記第２指標値は、前記走査線のうち前記計算範囲における輝度値に基づいて推定されたレイリー分布のパラメータに関する最尤推定値の関数である
　請求項６または請求項７に記載のプログラム。
　前記第３指標値は、
　　前記走査線における輝度値の度数分布を算出し、
　　前記度数分布のうち、輝度値が所定の第１閾値より低い範囲について、輝度値の第１統計量を算出し、
　　算出した前記第１統計量に基づいて、第２閾値を算出し、
　　前記度数分布のうち、輝度値が前記第２閾値よりも高い範囲について算出した第２統計量である
　請求項６から請求項８のいずか一つに記載のプログラム。
　前記第１統計量は、平均値および標準偏差であり、
　前記第２統計量は、標準偏差である
　請求項９に記載のプログラム。
　前記指標値に基づいて定めた色は、
　　前記第１指標値、前記第２指標値および前記第３指標値をそれぞれ０から２５５までの整数に正規化し、
　　正規化した前記第１指標値、前記第２指標値および前記第３指標値を、カラー画像を構成する複数のカラーチャネルそれぞれの輝度に、重複しない任意の順番で割り当てて生成される
　請求項６から請求項１０のいずか一つに記載のプログラム。
　前記計算範囲は、前記交点を含む範囲である
　請求項１から請求項１１のいずれか一つに記載のプログラム。
　前記画像取得用カテーテルの中心と前記交点との間の距離をＲ１、所定の定数をｂで示す場合に、前記計算範囲は、前記画像取得用カテーテルの中心からの距離が、（Ｒ１－ｂ／２）以上、画像取得範囲の最大値以下である
　請求項１から請求項１２のいずれか一つに記載のプログラム。
　前記画像取得用カテーテルの中心と前記交点との間の距離をＲ１、所定の定数をｂで示す場合に、前記計算範囲は、前記画像取得用カテーテルの中心からの距離が、（Ｒ１－ｂ／２）以上、（Ｒ１＋ｂ）以下である
　請求項１から請求項１２のいずれか一つに記載のプログラム。
　前記断層像の位置および前記交点の位置に加えて、第３の軸により定めた三次元空間における座標に、前記指標値に基づいて定めた色を表示する
　請求項１から請求項１４のいずれか一つに記載のプログラム。
　前記第３の軸は、前記画像取得用カテーテルと前記交点との間の距離である
　請求項１５に記載のプログラム。
　走査面を軸方向に移動させながら画像を取得する画像取得用カテーテルを用いた三次元走査により生成された複数の断層像を取得し、
　それぞれの前記断層像について代表線を決定し、
　それぞれの前記断層像を構成するそれぞれの走査線と前記代表線との交点に基づいて決定した計算範囲における、前記走査線の指標値を算出し、
　前記断層像の位置と、前記交点の位置とにより定めた座標に、前記指標値に基づいて定めた色を表示する
　処理をコンピュータが実行する情報処理方法。
　制御部を備える情報処理装置であって、
　前記制御部は、
　　走査面を軸方向に移動させながら画像を取得する画像取得用カテーテルを用いた三次元走査により生成された複数の断層像を取得し、
　　それぞれの前記断層像について代表線を決定し、
　　それぞれの前記断層像を構成するそれぞれの走査線と前記代表線との交点に基づいて決定した計算範囲における、前記走査線の指標値を算出し、
　　前記断層像の位置と、前記交点の位置とにより定めた座標に、前記指標値に基づいて定めた色を表示する
　情報処理装置。