(57)【特許請求の範囲】
【請求項１】
血管の３次元形状を表した３次元血管データに基づいて、血管の中心線を特定する中心線特定手段と、
前記中心線特定手段によって特定された前記中心線の中から、血管の突出部の中心線を特定する突出部中心線特定手段と、
前記突出部中心線特定手段によって特定された突出部中心線から分岐する中心線を特定する分岐中心線特定手段と、
前記分岐中心線特定手段によって特定された分岐中心線に対応する血管の前記３次元血管データに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義し、該円筒座標系で表された血管データを、直交座標系の血管データに変換する変換手段と、
前記変換手段によって変換された前記直交座標系の血管データにおいて、前記中心線方向の座標軸の第一の座標値における前記径方向の座標軸の座標値と、前記中心線方向の座標軸の第二の座標値における前記径方向の座標軸の座標値とを用いて、前記第一の座標値から前記第二の座標値までの間の前記径方向の座標軸の座標値を補間する補間手段と、
前記補間手段によって補間された後の前記直交座標系の血管データを、前記円筒座標系の血管データに変換する逆変換手段とを備えることを特徴とする血管画像処理装置。
【請求項２】
請求項１に記載の血管画像処理装置において、
前記３次元血管データと、前記逆変換手段によって逆変換が行われた後の前記円筒座標系の血管データとの差分を算出することにより、前記突出部の３次元データを生成する突出部データ生成手段をさらに備えることを特徴とする血管画像処理装置。
【請求項３】
請求項１または２に記載の血管画像処理装置において、
前記中心線特定手段は、前記３次元血管データに基づいて、血管を外周側から中心側に向けて同心円状に細めていき、所定の細さの線に達したときに、該線を前記中心線として特定することを特徴とする血管画像処理装置。
【請求項４】
請求項１〜３のいずれか一項に記載の血管画像処理装置において、
前記突出部中心線特定手段は、前記中心線特定手段によって特定された前記中心線に対して、前記中心線の分岐点を特定し、少なくとも一方の端に前記分岐点がない中心線を、前記血管の突出部の中心線として特定することを特徴とする血管画像処理装置。
【請求項５】
請求項４に記載の血管画像処理装置において、
前記突出部中心線特定手段は、前記中心線特定手段によって特定された前記中心線のうち、前記３次元血管データの画像領域の外周に接している中心線を除外して、前記血管の突出部を特定することを特徴とする血管画像処理装置。
【請求項６】
請求項４または５に記載の血管画像処理装置において、
前記突出部中心線特定手段は、前記分岐点から前記中心線が途切れる点までの長さを中心線の長さとして特定し、特定した中心線の長さを、前記中心線から血管の外周までの半径で割った値が、所定の閾値以下である中心線を除外して、前記血管の突出部を特定することを特徴とする血管画像処理装置。
【請求項７】
血管の３次元形状を表した３次元血管データに基づいて、血管の中心線を特定する中心線特定手順と、
前記中心線特定手順で特定した前記中心線の中から、血管の突出部の中心線を特定する突出部中心線特定手順と、
前記突出部中心線特定手順で特定した突出部中心線から分岐する中心線を特定する分岐中心線特定手順と、
前記分岐中心線特定手順で特定した分岐中心線に対応する血管の前記３次元血管データに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義し、該円筒座標系で表された血管データを、直交座標系の血管データに変換する変換手順と、
前記変換手順で変換した前記直交座標系の血管データにおいて、前記中心線方向の座標軸の第一の座標値における前記径方向の座標軸の座標値と、前記中心線方向の座標軸の第二の座標値における前記径方向の座標軸の座標値とを用いて、前記第一の座標値から前記第二の座標値までの間の前記径方向の座標軸の座標値を補間する補間手順と、
前記補間手順で補間した後の前記直交座標系の血管データを、前記円筒座標系の血管データに変換する逆変換手順とをコンピュータに実行させるための血管画像処理プログラム。
【請求項８】
請求項７に記載の血管画像処理プログラムにおいて、
前記３次元血管データと、前記逆変換手順で逆変換を行った後の前記円筒座標系の血管データとの差分を算出することにより、前記突出部の３次元データを生成する突出部データ生成手順をさらに有することを特徴とする血管画像処理プログラム。
【請求項９】
請求項７または８に記載の血管画像処理プログラムにおいて、
前記中心線特定手順は、前記３次元血管データに基づいて、血管を外周側から中心側に向けて同心円状に細めていき、所定の細さの線に達したときに、該線を前記中心線として特定することを特徴とする血管画像処理プログラム。
【請求項１０】
請求項７〜９のいずれか一項に記載の血管画像処理プログラムにおいて、
前記突出部中心線特定手順は、前記中心線特定手順で特定した前記中心線に対して、前記中心線の分岐点を特定し、少なくとも一方の端に前記分岐点がない中心線を、前記血管の突出部の中心線として特定することを特徴とする血管画像処理プログラム。
【請求項１１】
請求項１０に記載の血管画像処理プログラムにおいて、
前記突出部中心線特定手順は、前記中心線特定手順で特定した前記中心線のうち、前記３次元血管データの画像領域の外周に接している中心線を除外して、前記血管の突出部を特定することを特徴とする血管画像処理プログラム。
【請求項１２】
請求項１０または１１に記載の血管画像処理プログラムにおいて、
前記突出部中心線特定手段は、前記分岐点から前記中心線が途切れる点までの長さを中心線の長さとして特定し、特定した中心線の長さを、前記中心線から血管の外周までの半径で割った値が、所定の閾値以下である中心線を除外して、前記血管の突出部を特定することを特徴とする血管画像処理プログラム。
【請求項１３】
中心線特定手段が、血管の３次元形状を表した３次元血管データに基づいて、血管の中心線を特定するステップと、
突出部中心線特定手段が、前記中心線特定手段が特定した前記中心線の中から、血管の突出部の中心線を特定するステップと、
分岐中心線特定手段が、前記突出部中心線特定手段が特定した突出部中心線から分岐する中心線を特定するステップと、
変換手段が、前記分岐中心線特定手段が特定した分岐中心線に対応する血管の前記３次元血管データに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義し、該円筒座標系で表された血管データを、直交座標系の血管データに変換するステップと、
補間手段が、前記変換手段が変換した前記直交座標系の血管データにおいて、前記中心線方向の座標軸の第一の座標値における前記径方向の座標軸の座標値と、前記中心線方向の座標軸の第二の座標値における前記径方向の座標軸の座標値とを用いて、前記第一の座標値から前記第二の座標値までの間の前記径方向の座標軸の座標値を補間するステップと、
逆変換手段が、前記補間手段で補間した後の前記直交座標系の血管データを、前記円筒座標系の血管データに変換するステップとで構成される血管画像処理方法。
【請求項１４】
請求項１３に記載の血管画像処理方法において、
突出部データ生成手段が、前記３次元血管データと、前記逆変換手段で逆変換を行った後の前記円筒座標系の血管データとの差分を算出することにより、前記突出部の３次元データを生成するステップをさらに有することを特徴とする血管画像処理方法。
【請求項１５】
請求項１３または１４に記載の血管画像処理方法において、
前記中心線特定手段は、前記３次元血管データに基づいて、血管を外周側から中心側に向けて同心円状に細めていき、所定の細さの線に達したときに、該線を前記中心線として特定することを特徴とする血管画像処理方法。
【請求項１６】
請求項１３〜１５のいずれか一項に記載の血管画像処理方法において、
前記突出部中心線特定手段は、前記中心線特定手段が特定した前記中心線に対して、前記中心線の分岐点を特定し、少なくとも一方の端に前記分岐点がない中心線を、前記血管の突出部の中心線として特定することを特徴とする血管画像処理方法。
【請求項１７】
請求項１６に記載の血管画像処理方法において、
前記突出部中心線特定手段は、前記中心線特定手段が特定した前記中心線のうち、前記３次元血管データの画像領域の外周に接している中心線を除外して、前記血管の突出部を特定することを特徴とする血管画像処理方法。
【請求項１８】
請求項１６または１７に記載の血管画像処理方法において、
前記突出部中心線特定手段は、前記分岐点から前記中心線が途切れる点までの長さを中心線の長さとして特定し、特定した中心線の長さを、前記中心線から血管の外周までの半径で割った値が、所定の閾値以下である中心線を除外して、前記血管の突出部を特定することを特徴とする血管画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血管画像処理装置、血管画像処理プログラム、および血管画像処理方法に関する。
【背景技術】
【０００２】
次のような動脈瘤診断支援装置が知られている。この動脈瘤診断支援装置では、３次元画像データ上で脳動脈領域の内部に設定した基準点を始点として芯線データを生成し、生成した芯線データの連続性に基づいて正常脳動脈領域と脳動脈瘤領域を分離することにより脳動脈瘤を抽出する（例えば、特許文献１）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−１１０４４４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
脳動脈を切除するに当たっては、事前に正常な主血管の形状を正確に把握しておくことが好ましい。しかしながら、従来の動脈瘤診断支援装置では、脳動脈領域における脳動脈瘤を抽出して特定することはできるが、正常な主血管の形状を特定するための技術については、何ら検討されていなかった。
【課題を解決するための手段】
【０００５】
本発明による血管画像処理装置は、血管の３次元形状を表した３次元血管データに基づいて、血管の中心線を特定する中心線特定手段と、中心線特定手段によって特定された中心線の中から、血管の突出部の中心線を特定する突出部中心線特定手段と、突出部中心線特定手段によって特定された突出部中心線から分岐する中心線を特定する分岐中心線特定手段と、分岐中心線特定手段によって特定された分岐中心線に対応する血管の前記３次元血管データに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義し、該円筒座標系で表された血管データを、直交座標系の血管データに変換する変換手段と、変換手段によって変換された直交座標系の血管データにおいて、中心線方向の座標軸の第一の座標値における径方向の座標軸の座標値と、中心線方向の座標軸の第二の座標値における径方向の座標軸の座標値とを用いて、第一の座標値から第二の座標値までの間の径方向の座標軸の座標値を補間する補間手段と、補間手段によって補間された後の直交座標系の血管データを、円筒座標系の血管データに変換する逆変換手段とを備えることを特徴とする。
本発明による血管画像処理プログラムは、血管の３次元形状を表した３次元血管データに基づいて、血管の中心線を特定する中心線特定手順と、中心線特定手順で特定した中心線の中から、血管の突出部の中心線を特定する突出部中心線特定手順と、突出部中心線特定手順で特定した突出部中心線から分岐する中心線を特定する分岐中心線特定手順と、分岐中心線特定手順で特定した分岐中心線に対応する血管の３次元血管データに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義し、該円筒座標系で表された血管データを、直交座標系の血管データに変換する変換手順と、変換手順で変換した直交座標系の血管データにおいて、中心線方向の座標軸の第一の座標値における径方向の座標軸の座標値と、中心線方向の座標軸の第二の座標値における径方向の座標軸の座標値とを用いて、第一の座標値から第二の座標値までの間の径方向の座標軸の座標値を補間する補間手順と、補間手順で補間した後の直交座標系の血管データを、円筒座標系の血管データに変換する逆変換手順とをコンピュータに実行させるためのプログラムである。
本発明による血管画像処理方法は、中心線特定手段が、血管の３次元形状を表した３次元血管データに基づいて、血管の中心線を特定するステップと、突出部中心線特定手段が、中心線特定手段が特定した中心線の中から、血管の突出部の中心線を特定するステップと、分岐中心線特定手段が、突出部中心線特定手段が特定した突出部中心線から分岐する中心線を特定するステップと、変換手段が、分岐中心線特定手段が特定した分岐中心線に対応する血管の３次元血管データに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義し、該円筒座標系で表された血管データを、直交座標系の血管データに変換するステップと、補間手段が、変換手段が変換した直交座標系の血管データにおいて、中心線方向の座標軸の第一の座標値における径方向の座標軸の座標値と、中心線方向の座標軸の第二の座標値における径方向の座標軸の座標値とを用いて、第一の座標値から第二の座標値までの間の径方向の座標軸の座標値を補間するステップと、逆変換手段が、補間手段で補間した後の直交座標系の血管データを、円筒座標系の血管データに変換するステップとで構成される。
【発明の効果】
【０００６】
本発明によれば、突出部中心線から分岐する中心線に対応する血管の直交座標系の血管データにおいて、中心線方向の座標軸の第一の座標値における径方向の座標軸の座標値と、中心線方向の座標軸の第二の座標値における径方向の座標軸の座標値とを用いて、第一の座標値から第二の座標値までの間の径方向の座標軸の座標値を補間し、補間した後の直交座標系の血管データを円筒座標系の血管データに変換するようにしたので、血管から突出部を除いた正常な主血管の３次元データを生成して、正常な主血管の形状を特定することができる。
【図面の簡単な説明】
【０００７】
【図１】血管画像処理装置１００の一実施の形態の構成を示すブロック図である。
【図２】３次元血管データの一例を模式的に示した図である。
【図３】３次元血管データに対して作成した血管の中心線を模式的に示す図である。
【図４】個別の中心線を分岐点に基づいて特定した場合の具体例を示す図である。
【図５】３次元血管データに基づいて中心線から血管の外周までの半径を算出する方法を模式的に示した図である。
【図６】中心線の一方の端が３次元血管データの画像領域の外周に接している血管と、中心線が途中で途切れている細い血管の具体例を示す図である。
【図７】突出部の中心線の特定例と、分岐中心線の特定例を模式的に示す図である。
【図８】処理対象ボクセルデータの具体例を示す図である。
【図９】処理対象ボクセルデータに対する円筒座標系の定義例を模式的に示した図である。
【図１０】処理対象ボクセルデータを模式的に示した図である。
【図１１】処理対象ボクセルデータを直交座標系のボクセルデータに展開した結果を模式的に示す図である。
【図１２】展開ボクセルデータの断面図を示す第１の図である。
【図１３】展開ボクセルデータの断面図を示す第２の図である。
【図１４】処理対象ボクセルデータの寸法例を示した図である。
【図１５】突出部を削除した展開ボクセルデータの一例を示す図である。
【図１６】処理対象ボクセルデータから正常な主血管のボクセルデータを生成するための画像処理の流れを示す図である。
【図１７】血管画像処理装置１００で実行される処理の流れを示すフローチャート図である。
【発明を実施するための形態】
【０００８】
図１は、本実施の形態における血管画像処理装置１００の一実施の形態の構成を示すブロック図である。血管画像処理装置１００としては、例えば、血管画像を画像処理するためのプログラムがインストールされたサーバ装置やパソコン等の画像処理装置が用いられ、図１は、血管画像処理装置１００としてサーバ装置を用いた場合の一実施の形態の構成を示している。血管画像処理装置１００は、接続インターフェース１０１と、制御装置１０２と、記録装置１０３とを備えている。
【０００９】
接続インターフェース１０１は、血管画像処理装置１００を他の装置や端末等の外部機器と接続するためのインターフェースである。例えば、接続インターフェース１０１には、血管画像処理装置１００に液晶ディスプレイ等の表示装置を接続するためのインターフェースや、血管画像処理装置１００をＬＡＮやインターネット等の通信回線に接続するためのインターフェースが含まれる。接続インターフェース１０１が血管画像処理装置１００に液晶ディスプレイ等の表示装置を接続するためのインターフェースである場合は、接続インターフェース１０１をされたモニタに情報を出力して表示することができる。また、接続インターフェース１０１が血管画像処理装置１００をＬＡＮやインターネット等の通信回線に接続するためのインターフェースである場合は、接続インターフェース１０１を介して、通信回線に接続された他の装置や機器と通信を行うことができる。
【００１０】
制御装置１０２は、ＣＰＵ、メモリ、およびその他の周辺回路によって構成され、血管画像処理装置１００の全体を制御する。なお、制御装置１０２を構成するメモリは、例えばＳＤＲＡＭ等の揮発性のメモリである。このメモリは、ＣＰＵがプログラム実行時にプログラムを展開するためのワークメモリや、データを一時的に記録するためのバッファメモリとして使用される。
【００１１】
記録装置１０３は、血管画像処理装置１００が蓄える種々のデータや、制御装置１０２が実行するためのプログラムのデータ等を記録するための記録装置であり、例えばＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）やＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）等が用いられる。なお、記録装置１０３に記録されるプログラムのデータは、ＣＤ−ＲＯＭやＤＶＤ−ＲＯＭなどの記録媒体に記録されて提供されたり、ネットワークを介して提供され、操作者が取得したプログラムのデータを記録装置１０３にインストールすることによって、制御装置１０２がプログラムを実行できるようになる。
【００１２】
本実施の形態における血管画像処理装置１００は、あらかじめ記録装置１０３に記録されている血管の３次元画像データに基づいて、血管の異常突出領域を特定し、特定した異常突出領域を除去した正常な主血管の形状を特定するための処理を行う。本実施の形態では、記録装置１０３には、あらかじめ脳の血管を撮影して得た３次元血管データが記録されており、脳血管疾患の一つである脳動脈瘤などの血管の異常突出領域を検出する場合の例について説明する。
【００１３】
記録装置１０３に記録される３次元血管データは、患者の診断時に撮影された画像データに基づいて生成されるものとする。一般的に、脳動脈瘤などの血管疾患の診察時には、ＣＴＡ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ  Ａｎｇｉｏｇｒａｐｈｙ）、３Ｄ−ＤＳＡ（Ｔｈｒｅｅ−Ｄｉｍｅｎｓｉｏｎａｌ  Ｄｉｇｉｔａｌ  Ｓｕｂｔｒａｃｔｉｏｎ  Ａｎｇｉｏｇｒａｐｈｙ）、ＭＲＡ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ａｎｇｉｏｇｒａｐｈｙ）による撮影診断が行われると，その撮影された画像は、ＤＩＣＯＭ（Ｄｉｇｉｔａｌ  Ｉｍａｇｉｎｇ  ａｎｄ  Ｃｏｍｍｕｎｉｃａｔｉｏｎ  ｉｎ  Ｍｅｄｉｃｉｎｅ)規格の画像データとして保存される。撮影診断に当たっては、撮影対象（例えば頭部）の断面が何層にも渡って撮影されるため、これらの断層画像を重ねることにより、撮影対象の３次元画像データを作成することができる。
【００１４】
本実施の形態では、このように撮影された撮影対象の３次元画像データの中から、血管形状のみを抽出した３次元血管データが記録装置１０３に記録されている。撮影対象の３次元画像データから血管形状のみを抽出した３次元血管データを作成する方法としては、例えば、以下のような方法が考えられる。
【００１５】
上記の断層画像には撮影対象が物質ごとに異なる濃さで映し出されているため、撮影対象の３次元画像データから、血管の濃さに対応した画素で構成される領域の画像データを抽出すれば、血管形状のみを抽出した３次元血管データを作成することができる。撮影対象の３次元画像データにおける色の濃淡は，各モダリティの撮影原理に基づき、撮影時に得られた数値の情報で決定される。例えば、ＣＴＡの場合は、Ｘ線の吸収されやすさ、すなわちＸ線透過率を検出して撮影を行う。Ｘ線透過率は水を０、空気を−１０００としたＣＴ値（単位：ＨＵ）で表され、この値の大きさによって色の濃淡が決まる。すなわち、ＣＴＡで撮影された撮影対象の３次元画像データでは、各画素には、Ｘ線透過率に応じて概ね−２０４８から＋２０００程度の値が割り振られている。
【００１６】
例えば、あらかじめ血管の濃さに対応する数値の範囲を設定しておき、この撮影対象の３次元画像データに対し、画素値が該範囲の下限値以下である画素は白、画素値が該範囲の上限値以上である画素は黒に変換するように画像処理を行う。これによって、画素値が血管の濃さに対応する数値の画素は、撮影時の色の濃さが保持され、それ以外の画素は白と黒に変換された３次元画像データを得ることができ、撮影対象の３次元画像データから血管領域を抽出することができる。ここで、血管領域を抽出するための画素値の閾値は、患者によって個人差があるため、処理実行時に患者に適した値が設定されるのが好ましい。
【００１７】
本実施の形態では、抽出した血管領域の内部、すなわち血管内腔をボクセルで詰めたボクセルデータを３次元血管データとして、記録装置１０３に記録しておくものとする。図２は、記録装置１０３に記録されている３次元血管データの一例を模式的に示した図である。
【００１８】
制御装置１０２は、記録装置１０３から３次元血管データを読み出して、血管の中心を通る中心線を作成する。３次元血管データにおいては、血管は略円筒状の形状として撮影されているため、本実施の形態では、ボクセルデータを外側から削っていき，残った中心部分のボクセルを結ぶことによって中心線を作成する。すなわち、制御装置１０２は、３次元血管データに基づいて、血管を外周側から中心側に向けて同心円状に細めていき、線の細さがあらかじめ設定された所定の細さに達したときに、その線を中心線として特定する。
【００１９】
図３は、３次元血管データに対して作成した血管の中心線を模式的に示す図である。図３に示す例では、血管３ａに対して、中心線３ｂが作成されている。
【００２０】
本実施の形態では、制御装置１０２は、作成した中心線が分岐する点を特定し、分岐点から分岐点までの間の中心線、及び分岐点から中心線が途切れるまでの間の中心線を個別の中心線として特定する。これにより、図３に示した中心線３ｂは、図４に示すように、中心線４ａ〜４ｅの５本の中心線として区別される。以下の処理では、この中心線４ａ〜４ｅを中心線と呼んで処理の対象とする。
【００２１】
制御装置１０２は、作成した中心線に基づいて、動脈瘤などの異常突出領域である可能性が高い血管の突出部を特定する。このために、制御装置１０２は、一方の端に分岐点を有さない中心線を突出部の中心線として特定する。なお、制御装置１０２は、一方の端に分岐点を有さない中心線であっても、その一方の端が３次元血管データの画像領域の外周に接している中心線は、特定候補から除外する。また、制御装置１０２は、一方の端に分岐点を有さない中心線であっても、その中心線が、途中で途切れている細い血管の中心線である場合は、特定候補から除外する。途中で途切れている細い血管の中心線であるか否かの判定は次のようにして行う。
【００２２】
制御装置１０２は、分岐点から中心線が途切れる点までの長さを中心線の長さＬ１として特定し、特定した中心線の長さＬ１と、該中心線から血管の外周までの半径Ｌ２との比が小さいものは、途中で途切れている細い血管の中心線であると判定する。例えば、中心線の長さＬ１を中心線から血管の外周までの半径Ｌ２で割った値が、所定の閾値以下、例えば１０以下である中心線を途中で途切れている細い血管の中心線であると判定して除外する。
【００２３】
なお、本実施の形態では、中心線上の各位置において、中心線から血管壁までの距離を血管の全円周方向に対して算出し、その平均値を中心線から血管の外周までの半径Ｌ２として特定する。このとき、平均値を算出するためのサンプル数をいくつにするか、すなわち、血管の全円周方向を何分割して、中心線から血管壁までの距離を算出するかは、あらかじめ設定されているものとする。一例としては、３０度間隔で１２分割して１２方向について中心線から血管壁までの距離を算出し、これらの平均を算出することで中心線から血管の外周までの半径Ｌ２とする方法が考えられる。
【００２４】
図５は、３次元血管データに基づいて中心線から血管の外周までの半径Ｌ２を算出する方法を模式的に示した図である。本実施の形態では、図５（ａ）に示すように、１、２、３の数字で示されている中心線上の各位置において、中心線から血管の外周までの半径を算出する。中心線上の各位置をどのような間隔で設定するかは、あらかじめ決められているものとする。中心線上の各位置における中心線から血管の外周までの半径は、図５（ｂ）に示すように、中心線から、Ａ、Ｂ、・・・Ｚ、ＡＡで表される血管の外周上の点までの距離を算出し、これらの平均値を算出することによって特定される。制御装置１０２は、中心線上の各位置ごとに特定した半径を平均することにより、血管の中心線から血管の外周までの半径Ｌ２を算出する。
【００２５】
上記の処理により、例えば、図６に示すように、中心線の一方の端が３次元血管データの画像領域の外周に接している血管６ａ、６ｂ、６ｃ、６ｄと、中心線が途中で途切れている細い血管６ｅは、突出部の特定候補から除外される。
【００２６】
制御装置１０２は、上記の処理によって動脈瘤などの異常突出領域である可能性が高い血管の中心線、すなわち突出部の中心線を特定することができる。本実施の形態では、上記の処理によって、図７（ａ）に示す中心線７ａが動脈瘤などの異常突出領域である可能性が高い血管の中心線、すなわち突出部の中心線であると特定された場合を例に、以下の処理について説明する。なお、本実施の形態では、中心線７ａに対する処理について説明するが、上記の処理で複数の中心線が突出部の中心線として特定された場合には、以下に説明する処理は、それぞれの突出部の中心線に対して実行される。
【００２７】
制御装置１０２は、突出部の中心線７ａから分岐している中心線を特定する。図７（ａ）に示す例では、中心線７ｂと中心線７ｃの２本が、突出部の中心線７ａから分岐している中心線として特定される。制御装置１０２は、突出部の中心線７ａから分岐している中心線を１本の中心線として定義し直し，以降の処理に使用する。これによって、中心線７ｂと中心線７ｃは、図７（ｂ）に示すように、１本の中心線７ｄとして再定義される。以下の説明においては、再定義した中心線７ｄを分岐中心線と呼んで、上記の中心線と区別する。
【００２８】
制御装置１０２は、再定義した分岐中心線７ｄを含む領域の３次元血管データを記録装置１０３から読み出す。これにより、分岐中心線７ｄを含む領域の３次元血管データとして、例えば、図８に示す血管のボクセルデータが読み出される。
【００２９】
制御装置１０２は、分岐中心線７ｄを含む領域の３次元血管データを処理対象ボクセルデータとし、この処理対象ボクセルデータに対して、分岐中心線７ｄを参照したボクセルデータの展開処理を行う。ここで、ボクセルデータの展開処理について説明する。制御装置１０２は、処理対象ボクセルデータに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義する。図９は、処理対象ボクセルデータに対して、分岐中心線７ｄから血管の周方向をｘ軸、分岐中心線７ｄから血管の径方向をｙ軸、分岐中心線７ｄ方向をｚ軸とする円筒座標系を定義した場合の具体例を模式的に示した図である。なお、図９（ａ）は、処理対象ボクセルデータを斜め方向から見た状態を示した図であり、図９（ｂ）は、処理対象ボクセルデータを分岐中心線７ｄ方向、すなわちｚ軸方向から見た状態を示した図である。
【００３０】
制御装置１０２は、図９に示すように円筒座標系で表された処理対象ボクセルデータを、直交座標系の血管データに変換する。本実施の形態では、制御装置１０２は、分岐中心線７ｄから血管の周方向の角度が０°である点をｘ軸の原点、分岐中心線７ｄ上をｙ軸の原点、及び分岐中心線７ｄの一方の端点をｚ軸の原点とした直交座標系に、処理対象ボクセルデータ上の各点をプロットしていくことにより、処理対象ボクセルデータを、直交座標系の血管データに変換する。これによって、円筒座標系で表された処理対象ボクセルデータを、直交座標系に展開したボクセルデータを得ることができる。なお、直交座標系の血管データにおけるｘ軸の原点は、元画像から引き継いだ座標系（以下、「グローバル座標系」と呼ぶ）において，中心線軸方向のベクトルと，３次元血管データのグローバル座標系のｙ軸方向ベクトルとの外積方向と定義する。このとき、中心線軸方向のベクトルとグローバル座標系のｙ軸方向ベクトルが一致する場合は、ｙ軸方向の代わりにｚ軸方向ベクトルを用いる。
【００３１】
図１０〜図１４を用いて、処理対象ボクセルデータを直交座標系のボクセルデータに展開する例について説明する。図１０は、処理対象ボクセルデータを説明のために模式的に示した図である。図１０（ａ）は、処理対象ボクセルデータの斜視図であり、図１０（ｂ）は、処理対象ボクセルデータを分岐中心線７ｄから血管の周方向の角度が１８０°の方向から見た図であり、図１０（ｃ）は、処理対象ボクセルデータを分岐中心線７ｄから血管の周方向の角度が２７０°の方向から見た図である。図１０に示す処理対象ボクセルデータでは、円筒形状の血管の上部に扇状の突出部が存在している血管の形状が模式的に表されている。
【００３２】
図１１は、図１０に示した処理対象ボクセルデータを直交座標系のボクセルデータに展開した結果を模式的に示す図である。図１１に示すように、直交座標系に展開した後の展開ボクセルデータでは、図１０に示した処理対象ボクセルデータの扇状の突出部に対応する領域のｙ値が大きくなるため、展開ボクセルデータにおいてもこの部分がｙ軸方向に突出している。
【００３３】
図１１に示した展開ボクセルデータを、図１２（ａ）において破線で示した切断面１２ａで切断をした場合の断面図を図１２（ｂ）に示す。切断面１２ａは、展開ボクセルデータの突出部を含むｚ位置に設定されているため、図１２（ｂ）の断面図においては、ｙ軸方向に突出部が表れている。
【００３４】
図１１に示した展開ボクセルデータを、図１３（ａ）において破線で示した切断面１３ａで切断をした場合の断面図を図１３（ｂ）に示す。切断面１３ａは、展開ボクセルデータの突出部を含まないｚ位置に設定されているため、図１３（ｂ）の断面図においては、ｙ軸方向に突出部は表れていない。
【００３５】
図１４は、図１０に示した処理対象ボクセルデータに対して寸法を付した図である。なお、ここに付した寸法は、説明のために例として付したものであり、実際の寸法を正確に表しているものではない。図１４（ａ）は、図１０に示した処理対象ボクセルデータを、図１０（ｂ）に示す切断線１０ａ、図１０（ｃ）に示す切断線１０ｂで切断した場合の断面を、切断線１０ａ、切断線１０ｂに付した矢印方向から見た状態を示している。図１４（ｂ）は、図１０（ｂ）に示した処理対象ボクセルデータに寸法を付した図である。
【００３６】
図１４（ａ）に示した寸法を加味すると、図１２（ｂ）に示す断面図においては、ｙ軸方向の値が低い部分の高さは２ｍｍに相当し、ｙ軸方向の値が高い部分の高さは６ｍｍに相当する。また、図１２（ｂ）に示す断面図において、ｘ軸方向の値は、扇状の突出部の開き角度に対応するため、ｘ軸方向の突出部の開始点のｘ座標値から終了点のｘ座標値までの間隔は、９０°に相当する。また、図１３（ｂ）に示す断面図においては、ｙ軸方向の高さは２ｍｍに相当する。図１４（ｂ）に示した寸法を加味すると、図１１に示した展開ボクセルデータにおいて、展開ボクセルデータのｚ軸方向の幅は１０ｍｍに相当し、突出部のｚ軸方向の幅は５ｍｍに相当する。
【００３７】
制御装置１０２は、展開ボクセルデータに対して、展開ボクセルデータ上に現われた突出部を削除するための処理を行う。具体的には、制御装置１０２は、展開ボクセルデータにおいて、分岐中心線７ｄ方向、すなわちｚ軸方向の始点と終点の２点をとり、この２点間のｙ座標値を線形補間などの補間方法を用いて補間する。これによって、図１１に示した展開ボクセルデータは、図１５に示すような突出部を削除した展開ボクセルデータに変換される。
【００３８】
制御装置１０２は、突出部を削除した展開ボクセルデータを、直交座標系から円筒座標系に変換することにより、処理対象ボクセルデータから突出部を削除した正常な主血管のボクセルデータを得ることができる。
【００３９】
上記で説明した処理対象ボクセルデータから正常な主血管のボクセルデータを生成するための画像処理の流れについて、図１６に具体例を示して説明する。制御装置１０２は、図１６（ａ）に示すように、突出部を含んだ処理対象ボクセルデータを、図１６（ｂ）に示すような直交座標系のボクセルデータに展開する。その後、制御装置１０２は、展開ボクセルデータに対して、展開ボクセルデータ上に現われた突出部を削除するために、ｚ軸方向の始点と終点の２点をとり、この２点間のｙ座標値を線形補間により補間する。これによって、展開ボクセルデータは、図１６（ｃ）に示す突出部を削除した展開ボクセルデータに変換される。制御装置１０２は、この突出部を削除した展開ボクセルデータを、直交座標系から円筒座標系に変換することにより、図１６（ｄ）に示すような、処理対象ボクセルデータから突出部を削除した正常な主血管のボクセルデータを生成する。
【００４０】
制御装置１０２は、さらに、処理対象ボクセルデータと正常な主血管のボクセルデータとの差分を算出して差分画像を得ることにより、動脈瘤などの異常突出部を抽出したボクセルデータを作成することもできる。
【００４１】
制御装置１０２は、作成した正常な主血管のボクセルデータや異常突出部を抽出したボクセルデータを記録装置１０３に記録する。また、制御装置１０２は、作成した正常な主血管のボクセルデータや異常突出部を抽出したボクセルデータを自動的に、または操作者による操作に基づいて、接続インターフェース１０１に接続された表示装置に出力する。これによって、医師等は、患者の正常な主血管の形状や異常突出部の形状を把握して、診断や治療に役立てることができる。
【００４２】
図１７は、本実施の形態における血管画像処理装置１００で実行される処理の流れを示すフローチャートである。図１７に示す処理は、血管画像処理装置１００の操作者によってプログラムの実行が指示されると起動するプログラムとして、制御装置１０２によって実行される。なお、図１７に示す処理においては、あらかじめ上述した３次元血管データが記録装置１０３に記録されているものとする。
【００４３】
ステップＳ１０において、制御装置１０２は、記録装置１０３から３次元血管データを読み出して、上述したように、血管の中心を通る中心線を作成する。その後、ステップＳ２０へ進む。
【００４４】
ステップＳ２０では、制御装置１０２は、上述した方法により、中心線に基づいて、動脈瘤などの異常突出領域である可能性が高い血管の突出部を特定する。その後、ステップＳ３０へ進む。
【００４５】
ステップＳ３０では、制御装置１０２は、上述したように、突出部の中心線７ａから分岐している中心線を特定し、これらを１本の中心線として、分岐中心線を定義する。
その後、ステップＳ４０へ進む。
【００４６】
ステップＳ４０では、制御装置１０２は、分岐中心線を含む領域の３次元血管データを処理対象ボクセルデータとし、この処理対象ボクセルデータに対して、上述したように、ボクセルデータの展開処理を行って、直交座標系のボクセルデータに変換する。その後、ステップＳ５０へ進む。
【００４７】
ステップＳ５０では、制御装置１０２は、上述したように、展開ボクセルデータにおいて、ｚ軸上に２点をとり、この２点間のｙ座標値を線形補間などの補間方法を用いて補間することによって、展開ボクセルデータから突出部を削除する。その後、ステップＳ６０へ進む。
【００４８】
ステップＳ６０では、制御装置１０２は、突出部を削除した展開ボクセルデータを、直交座標系から円筒座標系に変換する。その後、処理を終了する。
【００４９】
以上説明した本実施の形態によれば、以下のような作用効果を得ることができる。
（１）制御装置１０２は、血管の３次元形状を表した３次元血管データに基づいて、血管の中心線を特定し、中心線の中から血管の突出部の中心線を特定し、突出部中心線から分岐する中心線を特定するようにした。そして、制御装置１０２は、分岐中心線に対応する血管の３次元血管データに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義し、該円筒座標系で表された血管データを、直交座標系の血管データに変換し、直交座標系の血管データにおいて、中心線方向の座標軸の第一の座標値における径方向の座標軸の座標値と、中心線方向の座標軸の第二の座標値における径方向の座標軸の座標値とを用いて、第一の座標値から第二の座標値までの間の径方向の座標軸の座標値を補間して、補間した後の直交座標系の血管データを、円筒座標系の血管データに変換するようにした。これによって、血管データの中心線に基づいて簡易な処理で動脈瘤などの血管の突出部を検出することができる。また、直交座標系の血管データにおいて、血管の径方向の座標値を補間することにより突出部を削除することができるため、円筒座標系の血管データ上で突出部を削除するための画像処理を行う場合と比較して、簡易な処理で精度高く血管の突出部を削除することができる。さらに、補間した後の直交座標系の血管データを、円筒座標系の血管データに変換するようにしたので、突出部を削除した後の正常な主血管の形状を特定することができる。
【００５０】
（２）制御装置１０２は、３次元血管データと、突出部を削除した後の円筒座標系の血管データとの差分を算出することにより、突出部の３次元データを生成するようにした。これによって、動脈瘤などの血管の異常突出領域の形状を特定することができる。
【００５１】
（３）制御装置１０２は、３次元血管データに基づいて、血管を外周側から中心側に向けて同心円状に細めていき、所定の細さの線に達したときに、該線を中心線として特定するようにした。これによって、３次元血管データにおける血管の中心線を簡易な処理で精度高く特定することができる。
【００５２】
（４）制御装置１０２は、中心線の分岐点を特定し、少なくとも一方の端に分岐点がない中心線を、血管の突出部の中心線として特定するようにした。これにより、分岐点を有さずに途中で途切れている中心線は、突出部の中心線であることが高い可能性を加味して、精度高く突出部の中心線を特定することができる。
【００５３】
（５）制御装置１０２は、中心線のうち、３次元血管データの画像領域の外周に接している中心線を除外して、血管の突出部を特定するようにした。これによって、一方の端に分岐点を有さない中心線であっても、その端が３次元血管データの画像領域の外周に接している場合には、正常な血管である可能性が高いことを加味して、突出部の検出対象から除外することができため、正常な血管が突出部であると誤検出されることを防ぐことができる。
【００５４】
（６）制御装置１０２は、分岐点から中心線が途切れる点までの長さを中心線の長さとして特定し、特定した中心線の長さを、中心線から血管の外周までの半径で割った値が、所定の閾値以下である中心線を除外して、血管の突出部を特定するようにした。これによって、血管の長さに対する太さの比が小さい場合は、途中で途切れている細い血管である可能性が高いため、このような細い血管が突出部であると誤検出されることを防ぐことができる。
【００５５】
―変形例―
なお、上述した実施の形態の血管画像処理装置１００は、以下のように変形することもできる。
（１）上述した実施の形態では、制御装置１０２は、展開ボクセルデータにおいて、ｚ軸方向の始点と終点の２点をとり、この２点間のｙ座標値を線形補間などの補間方法を用いて補間することによって、展開ボクセルデータから突出部を削除する例について説明した。しかしながら、補間のために特定するｚ軸方向の２点は、展開ボクセルデータ上に現われた突出部の開始位置よりもｚ値が小さい点と、展開ボクセルデータ上に現われた突出部の終了位置よりもｚ値が大きい点の２点をとれば、展開ボクセルデータから突出部を削除することができる。このため、補間に用いる２点の取り方は、ｚ軸方向の始点と終点を取る方法に限定されるものではない。
【００５６】
（２）上述した実施の形態では、血管画像処理装置１００として図１に示した構成のサーバ装置を用いる例について説明した。しかしながら、血管画像処理装置１００の構成は、図１に示した構成に限定されず、上述した実施の形態における処理を実行することができる装置であれば、他の構成であってもよい。例えば、血管画像処理装置１００は、記録装置１０３を備えず、通信回線や有線ケーブルなどを介して接続される外部ストレージ機器に、種々のデータを記録するようにしてもよい。また、上述した実施の形態では、表示装置は、接続インターフェース１０１を介して接続される例について説明したが、血管画像処理装置１００が表示装置を内蔵するようにしてもよい。
【００５７】
なお、本発明の特徴的な機能を損なわない限り、本発明は、上述した実施の形態における構成に何ら限定されない。また、上述の実施の形態と複数の変形例を組み合わせた構成としてもよい。
【符号の説明】
【００５８】
１００  血管画像処理装置
１０１  接続インターフェース
１０２  制御装置
１０３  記録装置
【要約】      （修正有）【課題】血管の３次元画像から異常突出領域を削除すること。【解決手段】制御装置１０２は、血管の３次元形状を表した３次元血管データに基づいて、血管の中心線を特定し、中心線の中から血管の突出部の中心線を特定し、突出部中心線から分岐する中心線を特定する。そして、制御装置１０２は、分岐中心線に対応する血管の３次元血管データに対して、中心線から血管の周方向、中心線から血管の径方向、中心線方向を座標軸とした円筒座標系を定義し、該円筒座標系で表された血管データを、直交座標系の血管データに変換し、直交座標系の血管データにおいて、中心線方向の座標軸の第一の座標値における径方向の座標軸の座標値と、中心線方向の座標軸の第二の座標値における径方向の座標軸の座標値とを用いて、第一の座標値から第二の座標値までの間の径方向の座標軸の座標値を補間して、補間した後の直交座標系の血管データを、円筒座標系の血管データに変換する。【選択図】図１【図１】
0005890055.tif 000002
【図３】
0005890055.tif 000003
【図４】
0005890055.tif 000004
【図５】
0005890055.tif 000005
【図６】
0005890055.tif 000006
【図７】
0005890055.tif 000007
【図９】
0005890055.tif 000008
【図１０】
0005890055.tif 000009
【図１４】
0005890055.tif 000010
【図１７】
0005890055.tif 000011
【図２】
0005890055.tif 000012
【図８】
0005890055.tif 000013
【図１１】
0005890055.tif 000014
【図１２】
0005890055.tif 000015
【図１３】
0005890055.tif 000016
【図１５】
0005890055.tif 000017
【図１６】
0005890055.tif 000018