JP2012016617A

JP2012016617A - 眼底画像表示装置及びその制御方法、コンピュータプログラム

Info

Publication number: JP2012016617A
Application number: JP2011231091A
Authority: JP
Inventors: 康文 ▲高▼間; Yasufumi Takama
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 2011-10-20
Filing date: 2011-10-20
Publication date: 2012-01-26
Anticipated expiration: 2028-12-19
Also published as: JP5484425B2

Abstract

【課題】眼底血管の３次元の走行を医師が観察できるように可視化する技術を提供する。
【解決手段】網膜の断層像を示す眼底画像において網膜層のいずれかの境界位置を特定する特定手段と、不透明度のピーク位置が前記網膜内の所定位置となるように、特定された前記境界位置からの距離を不透明度を表すパラメータに変換するための距離伝達関数を設定する距離伝達関数設定手段と、前記眼底画像の輝度値を、不透明度を表すパラメータに変換する輝度伝達関数を設定する輝度伝達関数設定手段と、前記距離伝達関数と前記輝度伝達関数とを用いて、前記断層像の各位置の不透明度を算出し、ボリュームレンダリングにより半透明表示画像を生成する可視化手段とを備える。
【選択図】図５

Description

本発明は、眼底画像表示装置及びその制御方法、コンピュータプログラムに関する。

眼科医療分野においては眼底写真を撮影する眼底カメラが一般的に使われてきた。近年では、さらに、光干渉断層計(Optical Coherence Tomography、以下、OCTとする)が登場したことによって、網膜の断層像を撮影することが可能になった。この結果、網膜が層構造を形成していることや、その層構造が病気の進行によって崩れることなど、眼科医学として新しい知見が発見されつつある。現在では、特許文献１に開示されているように、撮影された網膜の複数の断層像から再構成された網膜ボリュームデータから層構造を抽出し、その情報を眼科の病気を診断するために利用したりしている。以下、眼底写真と網膜ボリュームデータを眼底画像と呼ぶ。

従来、眼科医は網膜の断層像やボリュームデータを用いて網膜の層構造を読影したり、眼底写真を用いて眼底血管や白斑などの状態を読影したりしている。特に、眼底血管は、体外から観察できる唯一の血管であり、眼底血管からは様々な疾病の兆候や現象が確認できる。中でも、眼底血管の動脈硬化によって塞がった血管が膨張し、近傍の血管に接触してしまう交叉現象は、血管破裂などを引き起こし、最悪、失明につながる恐れがある現象として知られている。したがって、眼科医にとって眼底血管の走行がわかることは診断をつけるために有意義であると言える。さらに、患者にとっては、病気の早期発見につながり、早期回復や失明などの回避が期待できる。このような理由から、特許文献２では、眼底写真や輝度値を奥行き方向に積算して生成した積算画像から２次元の血管領域を抽出する方法が開示されている。さらに、眼底写真から血管領域や白斑領域を抽出する方法は、それぞれ非特許文献１や非特許文献２のような方法が開示されている。

しかしながら、眼底写真では眼底血管の２次元的な走行しか観察することができなかった。それゆえ、交叉現象にみられる３次元的な眼底血管の重なりを直接観察することはできなかった。OCTが登場し高解像度の断層像から網膜ボリュームデータが再構成できる現在では、眼底血管の３次元的な走行が観察でき、交叉現象が直接観察できる可能性があると考えられる。

網膜ボリュームデータを観察するためには、ボクセルが持つ値を伝達関数によって不透明度や色に変換することによって、ボリュームデータを半透明に表示できるボリュームレンダリングと呼ばれる方法が有効である。透明度を左右する伝達関数は、例えば、図１に示すように横軸に輝度値を、縦軸に不透明度を定義した関数として表現される。通常、ユーザは、伝達関数の形やピークの位置、幅をユーザインターフェースによって手動で設定することができる。また特許文献３のように、ＣＴに写る臓器や血管などのＣＴ値ヒストグラムが器官毎に山を見せるのを利用し、ＣＴ値ヒストグラムにガウス関数をフィッティングし、その平均値と分散から計算したＣＴ値の範囲を不透明表示するように自動設計できる。
特開２００８−０７３０９９号公報特開２００７−３２５８３１公報特開２００８−６２７４公報エリサ・リッチ、レンゾ・パルフェッティ、"ラインオペレータ及びサポートベクトル分類を利用した網膜血管セグメンテーション"、ＩＥＥＥトランザクション、メディカルイメージング、第２６巻、第１０号、１３５７〜１３６５頁、２００７年（Elisa Ricci, Renzo Perfetti, "Retinal Blood Vessel Segmentation Using Line Operators and Support Vector Classification, " IEEE Transactions on Medical Imaging, Vol.26, No.10, pp.1357-1365, 2007.）トーマス・ウォルター、ジャン・クロード・クライン、パスカル・マッシン、アリ・エルジニー、"糖尿病性網膜症の診断に対する画像処理補助"ＩＥＥＥトランザクション、メディカルイメージング、第２１巻、第１０号、１２３６〜１２４３頁、２００２年１０月（Thomas Walter, Jean-Claude Klein, Pascale Massin and Ali Erginay: "A Contribution of Image Processing to the Diagnosis of Diabetic Retinopathy - Detection of Exudates in Color Fundus Images of the Human Retina, " IEEE Transactions on Medical Imaging, Vol.21, No.10, pp.1236-pp.1243, Oct. 2002.）

網膜ボリュームデータを用いて眼底血管の３次元走行を可視化するためには、次のような課題がある。

特許文献２や非特許文献１による方法は、２次元の眼底写真や積算画像の眼底血管を抽出する方法であって、眼底血管の３次元位置は特定できない。

特許文献３による方法では、ＯＣＴの輝度値ヒストグラムにおいて血管領域は山を見せないので、自動で眼底血管の３次元走行を可視化することは困難である。また、図２に示す網膜断層像２０１において、眼底血管２０２が走行する神経線維層下端付近は高輝度領域になる。さらに、眼底血管２０２では光がよく反射し、輝度値は高くなる。つまり、眼底血管２０２が走行する領域はコントラストが低く、輝度値を基に手動で伝達関数を設定しても、眼底血管だけを可視化することは困難である。

本発明は、上記の課題に鑑みてなされたものであり、網膜ボリュームデータから血管の３次元走行や白斑の３次元分布等を可視化することを可能とすることを目的とする。

上記の目的を達成するための本発明は、眼底画像表示装置であって、
網膜の断層像を示す眼底画像において網膜層のいずれかの境界位置を特定する特定手段と、
不透明度のピーク位置が前記網膜内の所定位置となるように、特定された前記境界位置からの距離を不透明度を表すパラメータに変換するための距離伝達関数を設定する距離伝達関数設定手段と、
前記眼底画像の輝度値を、不透明度を表すパラメータに変換する輝度伝達関数を設定する輝度伝達関数設定手段と、
前記距離伝達関数と前記輝度伝達関数とを用いて、前記断層像の各位置の不透明度を算出し、ボリュームレンダリングにより半透明表示画像を生成する可視化手段と
を備える。

本発明の構成により、網膜ボリュームデータから血管の３次元走行や白斑の３次元分布等を可視化することが可能となる。

ボリュームレンダリングのための伝達関数の一例を示す図である。網膜断層像を説明するための図である。発明の実施形態に係る眼底画像表示システムの構成の一例を示す図である。発明の実施形態に係る眼底写真、網膜ボリュームデータ及び積算画像の関係を説明するための図である。発明の実施形態に対応する眼底画像表示装置１の機能構成の一例を示す図である。発明の実施形態に対応する血管領域抽出部５０３の機能構成の一例を示す図である。発明の実施形態に対応する距離伝達関数設定部５０５の機能構成の一例を示す図である。発明の第１実施形態の眼底画像表示装置１における眼底血管の３次元走行をボリュームレンダリングによって可視化する処理手順を示すフローチャートである。発明の第１実施形態における、血管領域抽出部５０３による血管領域の抽出処理の詳細を示すフローチャートである。発明の第２実施形態に係る眼底画像表示装置１０００の機能構成の一例を示す図である。発明の第２実施形態に係る眼底画像表示装置１０００での網膜ボリュームデータの白斑を可視化する処理手順を示すフローチャートである。発明の実施形態に係る眼底画像表示装置のハードウェア構成の一例を示す図である。

以下、添付図面に従って本発明に係る眼底画像表示装置の好ましい実施形態について詳説する。ただし、発明の範囲は図示例に限定されるものではない。

［第１実施形態：血管の可視化］
図３は、本実施形態に係る眼底画像表示システムの構成の一例を示す図である。本実施形態において、眼底画像表示装置１は、ＬＡＮ３を介してデータベース２から眼底画像を読み出すことができる。或いは、眼底画像表示装置１に記憶装置、例えばＦＤＤ、ＣＤ−ＲＷドライブ、ＭＯドライブ、ＺＩＰドライブ等を接続し、それらのドライブから眼底画像を読み込むようにしても良い。また、ＬＡＮ３を経由して眼底画像撮影装置４から直接に医用画像等を取得してもよい。データベース２には、眼底画像の他にも、患者の名前や所見情報、図４に示す網膜４０３、神経線維層４０５、外網状層４０６、網膜色素上皮４０７の厚み情報が格納されている。図４は、眼底写真、網膜ボリュームデータ及び積算画像の関係を説明するための図である。

眼底画像撮影装置４としては、眼底の断層像を撮影できるＯＣＴや眼底の写真が撮影できる眼底カメラが挙げられる。また、ＯＣＴの種類としては、タイムドメイン方式のＴＤ−ＯＣＴやフーリエドメイン方式のＦＤ−ＯＣＴが挙げられる。一般的に、ＴＤ−ＯＣＴでは断層像が一枚撮像される。一方、高速撮像が可能なＦＤ−ＯＣＴでは、１回の撮影で複数の断層像が得られ、これらの断層像を順番に並べることによって、網膜ボリュームデータが再構成できる。眼底画像撮影装置４は、ユーザ（技師や医師）による操作に応じて被検者（患者）の眼底画像を撮像し、得られた眼底画像を眼底画像表示装置１へ出力する。また、眼底画像表示装置１は、眼底画像撮影装置４により得られた眼底画像等を格納するデータベース２に接続され、そこから必要な眼底画像等を取得するように構成されてもよい。なお、これらの機器との接続は、ＵＳＢやＩＥＥＥ１３９４等のインターフェイスを介して行ってもよい。また、インターネット等の外部ネットワークを介して接続される構成であってもよい。

次に、図５を用いて、本実施形態に眼底画像表示装置１の機能構成を説明する。図５は、発明の実施形態に対応する眼底画像表示装置１の機能構成の一例を示す図である。図５において、眼底画像表示装置１は、眼底画像入力部５０１、層抽出部５０２、血管領域抽出部５０３、距離算出部５０４、距離伝達関数設定部５０５、輝度伝達関数設定部５０６、可視化部５０７を備える。次に各部の動作を説明する。

眼底画像入力部５０１は、眼底画像撮影装置４から出力された眼底画像や、データベース２に蓄えられている眼底画像を選択し、入力する。図４に眼底画像として入力される網膜ボリュームデータ４０１と眼底写真４２１の模式図を示す。網膜内は層構造を成していることが医学的に知られている。網膜ボリュームデータ４０１の座標系は、網膜断層像４０２（模式図）の横方向をＸ軸、縦方向をＺ軸、網膜断層像４０２が並べられる方向にＹ軸とする。また、眼底写真４２１の座標系は、横方向をＸ軸、縦方向をＹ軸とする。

層抽出部５０２は、網膜ボリュームデータ４０１から、神経線維層４０５、外網状層４０６、網膜色素上皮４０７を抽出する。そして、内境界膜４０４、神経線維層４０５の下端、外網状層４０６の上端、網膜色素上皮４０７の上端、内境界膜４０４から網膜色素上皮４０７の下端までで定義される網膜４０３を求める。抽出する層は全ての層であっても良いし、１つもしくは、複数の層を選択しても良い。抽出された層の境界位置を示す層境界情報は、距離算出部５０４へ出力される。

血管領域抽出部５０３は、入力された眼底画像から血管領域を抽出する。図６に示すように、血管領域抽出部５０３は、積算画像血管領域抽出部６０１と、眼底写真血管領域抽出部６０２、血管位置合わせ部６０３から構成される。積算画像血管領域抽出部６０１は、網膜ボリュームデータ４０１のボクセルの輝度値をｚ軸方向へ積算した積算画像４１１から眼底の血管領域４１２を抽出する。さらに、眼底画像入力部５０１において眼底写真が入力された場合、眼底写真血管領域抽出部６０２は、眼底写真４２１から血管領域４２２を抽出する。血管位置合わせ部６０３は、眼底画像入力部５０１において眼底写真が入力されている場合、血管領域４１２と血管領域４２２を用いて眼底写真の座標系(X,Y)と積算画像の座標系(x,y)の位置合わせを行う。

これらの処理の結果、血管領域抽出部５０３は、血管領域を座標系(x,y)で表現されたボクセルの集合として、距離算出部５０４、距離伝達関数設定部５０５、輝度伝達関数設定部５０６へ出力する。

距離算出部５０４は、まず、層境界情報から１つの層境界を特定し、その境界位置を距離＝０とする。例えば、内境界膜が選択により特定された場合、図４に示す内境界膜４０４として抽出された層境界上の全ボクセルの距離を０とする。次に、距離算出部５０４は、血管領域抽出部５０３で得られた座標系(x,y)で表現されたボクセルの集合である血管領域４１２を用いて、網膜ボリュームデータの座標(x,y)を持つボクセルの集合である血管投影領域４０８を定義する。さらに、距離算出部５０４は、血管投影領域４０８の各ボクセルと、血管投影領域４０８と層境界との交点に存在する基準４０９との距離を算出し、算出した距離を各ボクセルに与える。このとき、基準４０９よりもｚが小さい座標を持つボクセルの距離の符号を「−」、ｚが大きい座標を持つボクセルの距離の符号を「＋」とする。

距離伝達関数設定部５０５は、距離算出部５０４で算出した各ボクセルの距離を不透明度に変換する距離伝達関数を設定する。距離伝達関数は、例えば図１において横軸に距離、縦軸に不透明度を定義した関数であり、血管や白斑など表示対象によって自動で設定する。距離伝達関数設定部５０５は、血管領域抽出部５０３で抽出した血管投影領域４０８において、基準４０９から眼底血管が走行すると考えられる眼底血管２０２までの距離を距離伝達関数のピーク位置に設定する。

ここで、眼底血管２０２が通っていると考えられる位置は、光を強く反射し、図２に示す血管の影領域２０３の輝度値が周囲に比べて低くなる。そこで、血管投影領域４０８内で最大輝度値を持つボクセル、もしくは、任意の大きさの矩形を設定しその矩形内の分散が大きくなるボクセルのｚ座標とする。距離伝達関数の分散は、血管領域から算出できる血管径から設定することも可能である。このようにして自動で設定された距離伝達関数は可視化部５０７へ出力される。

輝度伝達関数設定部５０６は、網膜ボリュームデータの輝度値を不透明度に変換する輝度伝達関数を設定する。輝度伝達関数は、例えば図１において横軸に輝度値、縦軸に不透明度を定義した関数である。この関数の設定は、ユーザがユーザインターフェースを用いて手動で行っても良いし、特許文献３のように輝度値のヒストグラムから自動で設定しても良い。設定された輝度伝達関数は、可視化部５０７へ出力される。

可視化部５０７は、輝度伝達関数によって算出した不透明度と距離伝達関数によって算出した不透明度を下記の式１により合成し、ボリュームレンダリング時にボクセルの不透明度を設定するために使用する。ここで、座標xにおけるボクセルの輝度値をｖ(x)、距離をｄ(x)、輝度伝達関数をｆ_v(v(x))、距離伝達関数をｆ_d(d(x))とする。合成された不透明度α(ｖ(x), ｄ(x))は、式１に示すように、輝度伝達関数と距離伝達関数の合成比βを用いて、輝度伝達関数と距離伝達関数の線形和によって得られる。

α(ｖ(x), ｄ(x)) =β* ｆ_v(ｖ(x)) + (1-β)*ｆ_d(ｄ(x)) （式１）
なお、βは、表示対象毎にあらかじめ設定しておいても良い。血管を可視化したい場合、血管付近ではコントラストが低いことから、βを低く設定し、距離伝達関数の重みを重くすることが考えられる。最後に、可視化部５０７は、ボリュームレンダリングの結果画像として、半透明表示画像を生成する。

次に、以上で説明した本実施形態の眼底画像表示装置１における眼底血管の３次元走行をボリュームレンダリングによって可視化する処理手順を図８に示すフローチャートを用いて説明する。なお、図８のフローチャートによって示される処理は、図１２に示すＣＰＵ５が主メモリ６に格納されているプログラムを実行することにより実現される。

ステップＳ８０１において、眼底画像入力部５０１は、網膜ボリュームデータや眼底写真といった眼底画像を眼底画像表示装置１に入力する。入力された眼底画像は、層抽出部５０２、血管領域抽出部５０３、距離伝達関数設定部５０５、輝度伝達関数設定部５０６、可視化部５０７へ出力される。

ステップＳ８０２において、血管領域抽出部５０３は、網膜ボリュームデータ及び眼底画像から血管領域４１２及び４２２を抽出する。本ステップにおける処理は、図９を参照して後段でより詳細に説明する。

ステップＳ８０３において、層抽出部５０２は、ステップＳ８０１で受け取った網膜ボリュームデータから神経線維層４０５、外網状層４０６、網膜色素上皮４０７を抽出する。層抽出部５０２は、これらの層の情報から内境界膜４０４、神経線維層４０５の下端、外網状層４０６の上端、網膜４０３を求める。不図示の指示部により、ユーザは基準としたい層境界を選択し、層境界の情報はステップＳ８０４へ出力される。

ステップＳ８０４において、距離算出部５０４は、選択された層境界を基準として、ステップ８０２で抽出した血管領域内で、層境界からｚ軸方向にある各ボクセルまでの距離を算出する。算出された距離情報は、可視化部５０７へ出力される。なお、血管領域以外のボクセルまでの距離を計算しない場合、つまり、距離＝０とする場合、血管領域以外のボクセルの不透明度は０になり、ボリュームレンダリングを高速に計算することが可能である。また、血管領域以外のボクセルまでの距離も計算する場合、層からの距離に応じた不透明度を計算しながら、ボリュームレンダリングすることにより半透明表示画像を生成することが可能である。

ステップＳ８０５において、輝度伝達関数設定部５０６と距離伝達関数設定部５０５とは、輝度伝達関数と距離伝達関数とをそれぞれ設定する。距離伝達関数設定部５０５は、血管の影領域を抽出し、内境界膜から血管が走行すると考えられる所定位置の座標までの距離をピーク位置とする。設定された輝度伝達関数と距離伝達関数は、可視化部５０７へ出力される。

ステップＳ８０６において、可視化部５０７は、輝度伝達関数と距離伝達関数とから得られる不透明度を式１により掛け合わせながらボリュームレンダリングする。

図９は、ステップＳ８０２における処理を詳しく説明したフローチャートである。まず、ステップＳ９０１において、積算画像血管領域抽出部６０１は、入力された網膜ボリュームデータ４０１の積算画像４１１を生成する。さらに、積算画像血管領域抽出部６０１は、ステップＳ９０２において、積算画像４１１から血管領域４１２を抽出し、血管位置合わせ部６０３へ出力する。

ここで、ステップＳ９０２における血管領域４１２の抽出方法としては、公知の任意の手法を用いて行うことが可能である。たとえば、積算画像の画素値を解析して、隣接する画素の画素値の差を演算し、この差が所定値よりも大きい隣接する画素を探索することにより、血管領域とその他の領域との境界領域を検出する。それにより、積算画像から血管領域を抽出することができる。この抽出処理では、積算画像における、血管領域とその他の領域との画素値（輝度値等）の差を利用する。

ステップＳ９０３において、眼底画像入力部５０１は、眼底写真４２１がステップＳ８０１において入力されているかどうかを判断する。もし、眼底写真４２１が入力されていれば（Ｓ９０３で「ＹＥＳ」）、ステップＳ９０４に移行する。もし、眼底写真４２１が入力されていなければ（Ｓ９０３で「ＮＯ」）、ステップＳ９０６に移行する。

ステップＳ９０４では、眼底写真血管領域抽出部６０２は眼底写真４２１から血管領域４２２を抽出し、血管位置合わせ部６０３へ出力する。ステップＳ９０３における血管領域４２２の抽出方法も、ステップＳ９０２におけるものと同様である。即ち、眼底写真４２１の画素値を解析して、隣接する画素の画素値の差を演算し、この差が所定値よりも大きい隣接する画素を探索することにより、血管領域とその他の領域との境界領域を検出する。それにより、眼底写真４２１から血管領域を抽出することができる。その後ステップＳ９０５において、血管位置合わせ部６０３は、血管領域４１２と血管領域４２２を位置合わせする。

続くステップＳ９０６において、血管位置合わせ部６０３は、位置合わせ結果を用いて、眼底写真４２１の血管領域４１２と積算画像４１１の血管領域４２２を統合し、統合された血管領域をステップＳ８０３へ渡す。血管領域４１２と血管領域４２２の統合は、例えば、二つの領域の論理積もしくは論理和を計算することによって実現できる。もし、眼底写真４２１が入力されていなければ、積算画像４１１から抽出した血管領域４１２をステップＳ８０３に渡す。ここでステップＳ８０２が終了する。

以上のように本実施形態では、眼底の血管が神経線維層の下端付近を走行するという知見を距離伝達関数で表現し、網膜断層像では眼底血管の輝度値が高くなるという特徴を輝度伝達関数で表現し、ボリュームレンダリングする。これにより、網膜ボリュームデータの眼底血管の３次元走行が可視化できる。

［第２実施形態：白斑の可視化］
上述の第１実施形態では、眼底画像から血管領域と層境界を抽出し、層境界からの距離と網膜ボリュームデータの輝度値から算出した不透明度を用いて血管をボリュームレンダリングすることによって、眼底血管の３次元走行を可視化した。これに対して、第２実施形態では、血管ではなく、眼底の特徴的な病変である白斑を可視化することを目的とする。

本実施形態に係る眼底画像表示装置及びそれに接続される機器の構成は、図２に示す第１実施形態の構成と同様である。但し、本実施形態の眼底画像表示装置の機能ブロックは、図１０に示すように第１実施形態の機能構成に白斑領域抽出部１００１を追加して構成される。

以下では、白斑領域抽出部１００１と、距離伝達関数設定部５０５の機能構成の一つである白斑中心探索部７０４について説明し、その他の機能構成は第１実施形態と同様の処理のため説明を省略する。

白斑領域抽出部１００１は、入力された眼底画像から白斑領域を抽出する。入力された眼底画像が網膜ボリュームデータ４０１の場合、網膜層の情報と輝度値をもとに白斑領域を抽出する。ここで、白斑が出現する位置は外網状層４０６の付近であるという医学的知見をもとに予測することができる。この知見を利用するために、まず、神経線維層４０５と網膜色素上皮４０７とを抽出する。さらに、網膜断層像４０２における白斑は周囲に比べて高い輝度値を持つことから、神経線維層４０５と網膜色素上皮４０７の間の高輝度領域を二値化などの画像処理により抽出し、網膜ボリュームデータ４０１の白斑領域とする。

また、網膜ボリュームデータ４０１の他に眼底写真４２１が入力されている場合、眼底写真４２１の方が白斑を良く捉えられる場合があることから、眼底写真４２１において、２次元の白斑領域を先に抽出する。ここでの抽出方法も、白斑は周囲に比べて高い輝度値を持つことを利用して、高輝度領域を二値化などの画像処理を行うことで実現できる。

次に、血管領域抽出部５０３で算出できる位置合わせ情報をもとに、眼底写真４２１の白斑領域を網膜ボリュームデータ４０１に逆投影することによって血管投影領域と同じような白斑投影領域を求める。さらに、白斑投影領域内のボクセルから前述したように、白斑領域を抽出する。この方法によって、網膜ボリュームデータ４０１に写らないが、眼底写真４２１に写る白斑を白斑領域に含めることが可能である。

距離伝達関数設定部５０５の白斑中心探索部７０４は、白斑領域抽出部１００１によって抽出した白斑領域において、白斑領域内のボクセルをｚ軸方向に探索する。その際、層境界情報から１つの層境界を選択し、選択した層境界（例えば内境界膜４０４）から白斑領域の中点に存在するボクセルまでの距離を算出し、距離伝達関数のピーク位置とする。

次に、図１１に示すフローチャートを用いて、網膜ボリュームデータの白斑を可視化する手順について説明する。なお、図１１のフローチャートによって示される処理は、ＣＰＵ５が主メモリ６に格納されているプログラムを実行することにより実現される。ステップＳ１１０１とステップＳ１１０２、ステップＳ１１０３、ステップＳ１１０５、ステップＳ１１０７は、夫々、図８のステップＳ８０１とステップＳ８０２、ステップＳ８０３、ステップＳ８０４、ステップＳ８０６と同じである。よって、ここでの説明は省略する。

ステップＳ１１０４において、入力された眼底画像が網膜ボリュームデータ４０１だけの場合、網膜ボリュームデータ４０１から３次元白斑領域を抽出する。網膜ボリュームデータ４０１に加えて、眼底写真４２１も入力された場合、まず、眼底写真４２１から白斑領域を抽出する。次に、血管位置合わせ部６０３において算出した血管領域位置合わせ情報をもとに、眼底写真４２１の白斑領域を網膜ボリュームデータ４０１に逆投影する。さらに、ステップＳ１１０３で抽出した層の情報を用いて、白斑の３次元領域を抽出する。抽出した白斑領域の情報は、距離算出部５０４へ出力される。

ステップＳ１１０６において、輝度伝達関数設定部５０６と距離伝達関数設定部５０５は、輝度伝達関数と距離伝達関数を設定する。距離伝達関数設定部５０５は、抽出した白斑領域の白斑中心に相当する所定位置を求める。そして、内境界膜４０４から白斑領域の各(x,y)座標上におけるｚ座標の該所定位置を距離伝達関数のピーク位置とする。設定された輝度伝達関数と距離伝達関数は、可視化部５０７へ出力される。

以上のように第２実施形態では、白斑が外網状層付近に出現するという知見を距離伝達関数で表現し、網膜断層像では白斑の輝度値が高くなるという特徴を輝度伝達関数で表現し、ボリュームレンダリングする。これにより、網膜ボリュームデータの白斑の３次元分布が可視化できる。

［第３実施形態］
本実施形態に係る眼底画像表示装置及びそれに接続される機器の構成は、図２に示す第１実施形態の構成と同様である。また、本実施形態に係る眼底画像表示装置の各部の機能をソフトウェアで実現するためのコンピュータの基本構成も、図１２に示す第１実施形態の構成と同様である。

なお、第１実施形態、第２実施形態の距離伝達関数設定部５０５は、図７に示すように、平均網膜厚取得部７０１と、関数パラメータ入力部７０２、影領域抽出部７０３、白斑中心探索部７０４、層境界取得部７０５を備えている。第１実施形態では、影領域抽出部７０３を、第２実施形態では、白斑中心探索部７０４を用いて、距離伝達関数を設定していた。本実施形態では、平均網膜厚取得部７０１、関数パラメータ入力部７０２、層境界取得部７０５のいずれか１つを用いた距離伝達関数設定方法について説明する。

平均網膜厚取得部７０１は、データベース２から各層の平均の厚みを取得し、該厚みに基づいて距離伝達関数のピーク位置に相当する所定位置を決定する。眼底血管を可視化する場合、眼底血管が神経線維層４０５の下端付近に出現することから、神経線維層４０５の平均の厚みを取得する。白斑を可視化する場合、外網状層４０６付近に白斑が出現することから、内境界膜４０４から外網状層４０６の上端までの平均の厚みを取得する。これら取得した厚みを有する所定位置を不透明度のピーク位置とし、不透明度が１．０になるように距離伝達関数を設定する。

関数パラメータ入力部７０２は、ユーザインターフェースを用いて手動でピーク位置を設定する。層境界取得部７０５は、層抽出部５０２から、各層の境界情報を取得し、距離伝達関数のピーク位置とする。眼底血管を可視化する場合、神経線維層の境界情報を取得し、神経線維層の厚みが不透明度＝１．０となるようにピーク位置を設定する。白斑を可視化する場合、外網状層の境界情報を取得し、内境界膜から外網状層上端までの厚みが不透明度＝１．０となるようにピーク位置を設定する。

このように距離伝達関数設定部５０５を構成する、平均網膜厚取得部７０１と、関数パラメータ入力部７０２、影領域抽出部７０３、白斑中心探索部７０４、層境界取得部７０５のいずれか１つを用いて距離伝達関数を設定する。これにより、眼底血管の３次元走行や白斑の３次元分布が可視化できる。

［その他の実施形態］
次に、第１実施形態、第２実施形態及び第３実施形態の眼底画像表示装置の各ブロックの機能をコンピュータプログラムで実現するためのコンピュータの基本構成を、図１２を用いて説明する。

ＣＰＵ５は主メモリ６に格納されたデータやプログラムを用いてコンピュータ全体を制御する。また、眼底画像表示装置１の各部に対応するソフトウェアの実行を制御して、各部の機能を実現する。

主メモリ６は、ＣＰＵ５が実行する制御プログラムを格納したり、ＣＰＵ５によるプログラム実行時の作業領域を提供したりする。

磁気ディスク７は、オペレーティングシステム（ＯＳ）、周辺機器のデバイスドライブ、眼底画像を表示するためのプログラムを含む各種アプリケーションソフト等を格納する。表示メモリ８は、モニタ９のための表示用データを一時記憶する。

モニタ９は、例えばＣＲＴモニタや液晶モニタ等であり、表示メモリ８からのデータに基づいて画像を表示する。マウス１０及びキーボード１１はユーザによるポインティング入力及び文字等の入力をそれぞれ行う。操作者はこれらを用いて、各種の指示を眼底画像表示装置１に与えることができる。上記各構成要素は共通バス１２により互いに通信可能に接続されている。

以上、実施形態を詳述したが、本発明は、例えば、システム、装置、方法、プログラムもしくは記憶媒体等としての実施態様をとることが可能である。具体的には、複数の機器から構成されるシステムに適用しても良いし、また、一つの機器からなる装置に適用しても良い。

尚、本発明は、ソフトウェアプログラムをシステム或いは装置に直接或いは遠隔から供給し、システム或いは装置のコンピュータが該供給されたプログラムコードを読み出して実行することによって前述した実施形態の機能が達成される場合を含む。この場合、供給されるプログラムは実施形態で図に示したフローチャートに対応したコンピュータプログラムである。

従って、本発明の機能処理をコンピュータで実現するために、該コンピュータにインストールされるプログラムコード自体も本発明を実現するものである。つまり、本発明は、本発明の機能処理を実現するためのコンピュータプログラム自体も含まれる。その場合、プログラムの機能を有していれば、オブジェクトコード、インタプリタにより実行されるプログラム、ＯＳに供給するスクリプトデータ等の形態であっても良い。

コンピュータプログラムを供給するためのコンピュータ読み取り可能な記憶媒体としては以下が挙げられる。例えば、フロッピー（登録商標）ディスク、ハードディスク、光ディスク、光磁気ディスク、ＭＯ、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＣＤ−ＲＷ、磁気テープ、不揮発性のメモリカード、ＲＯＭ、ＤＶＤ（ＤＶＤ−ＲＯＭ、ＤＶＤ−Ｒ）などである。

その他、プログラム供給方法には、クライアントコンピュータのブラウザを用いてインターネットのホームページに接続し、該ホームページから本発明のコンピュータプログラムをハードディスク等の記録媒体にダウンロードすることが含まれる。この場合、ダウンロードされるプログラムは、圧縮され自動インストール機能を含むファイルであってもよい。また、本発明のプログラムを構成するプログラムコードを複数のファイルに分割し、それぞれのファイルを異なるホームページからダウンロードすることによっても実現可能である。つまり、本発明の機能処理をコンピュータで実現するためのプログラムファイルを複数のユーザに対してダウンロードさせるＷＷＷサーバも、本発明に含まれるものである。

また、本発明のプログラムを暗号化してＣＤ−ＲＯＭ等の記憶媒体に格納してユーザに配布するという形態をとることもできる。この場合、所定の条件をクリアしたユーザに、インターネットを介してホームページから暗号を解く鍵情報をダウンロードさせ、その鍵情報を使用して暗号化されたプログラムを実行し、プログラムをコンピュータにインストールさせるようにもできる。

また、コンピュータが、読み出したプログラムを実行することによって、前述した実施形態の機能が実現される他、そのプログラムの指示に基づき、コンピュータ上で稼動しているＯＳなどとの協働で実施形態の機能が実現されてもよい。この場合、ＯＳなどが、実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される。

さらに、記録媒体から読み出されたプログラムが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書き込まれて前述の実施形態の機能の一部或いは全てが実現されてもよい。この場合、機能拡張ボードや機能拡張ユニットにプログラムが書き込まれた後、そのプログラムの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部または全部を行う。

１眼底画像表示装置
２データベース
３ LAN
４眼底画像表示装置
５ CPU
６主メモリ
７磁気ディスク
８表示メモリ
９モニタ
１０マウス
１１キーボード
１２共通バス

Claims

網膜の断層像を示す眼底画像において網膜層のいずれかの境界位置を特定する特定手段と、
不透明度のピーク位置が前記網膜内の所定位置となるように、特定された前記境界位置からの距離を不透明度を表すパラメータに変換するための距離伝達関数を設定する距離伝達関数設定手段と、
前記眼底画像の輝度値を、不透明度を表すパラメータに変換する輝度伝達関数を設定する輝度伝達関数設定手段と、
前記距離伝達関数と前記輝度伝達関数とを用いて、前記断層像の各位置の不透明度を算出し、ボリュームレンダリングにより半透明表示画像を生成する可視化手段と
を備えることを特徴とする眼底画像表示装置。
前記眼底画像において、前記網膜内の眼底血管が存在する領域を抽出する血管領域抽出手段をさらに備え、
前記所定位置は、前記網膜内の眼底血管の存在する領域の位置であることを特徴とする請求項１に記載の眼底画像表示装置。
前記眼底画像には、眼底写真が更に含まれ、
前記血管領域抽出手段は前記眼底写真において眼底血管の位置を特定し、該特定した眼底血管の位置を利用して、前記眼底画像において前記網膜内の眼底血管が存在する領域を抽出することを特徴とする請求項２に記載の眼底画像表示装置。
前記眼底画像において、前記網膜内の白斑が存在する領域を抽出する白斑領域抽出手段をさらに備え、
前記所定位置は、前記網膜内の白斑の存在する領域の位置であることを特徴とする請求項１に記載の眼底画像表示装置。
前記眼底画像には、眼底写真が更に含まれ、
前記白斑領域抽出手段は前記眼底写真において白斑の位置を特定し、該特定した白斑の位置を利用して、前記眼底画像において前記網膜内の白斑が存在する領域を抽出することを特徴とする請求項４に記載の眼底画像表示装置。
前記網膜層を構成する各層の厚みを登録するデータベースをさらに備え、
前記距離伝達関数設定手段は、前記網膜内の所定位置を前記各層の厚みにより特定し、該所定位置が不透明度のピークとなるように、前記特定された境界位置からの距離を不透明度を表すパラメータに変換することを特徴とする請求項１に記載の眼底画像表示装置。
前記所定位置は、前記網膜内の他の境界の位置であることを特徴とする請求項１に記載の眼底画像表示装置。
網膜の断層像を示す眼底画像において網膜層のいずれかの境界位置を特定する特定工程と、
不透明度のピーク位置が前記網膜内の所定位置となるように、特定された前記境界位置からの距離を不透明度を表すパラメータに変換するための距離伝達関数を設定する距離伝達関数設定工程と、
前記眼底画像の輝度値を、不透明度を表すパラメータに変換する輝度伝達関数を設定する輝度伝達関数設定工程と、
前記距離伝達関数と前記輝度伝達関数とを用いて、前記断層像の各位置の不透明度を算出し、ボリュームレンダリングにより半透明表示画像を生成する可視化工程と
を備えることを特徴とする眼底画像表示装置の制御方法。
コンピュータを請求項１乃至請求項７のいずれか１項に記載の眼底画像表示装置を機能させるためのコンピュータプログラム。