(57)【特許請求の範囲】
【請求項１】
眼底の層構造をそれぞれ描写する複数の断層像を記憶する記憶手段と、
前記記憶された各断層像の画素の画素値に基づいて、網膜色素上皮層に相当する当該断層像中の層領域を特定する層領域特定手段と、
前記特定された層領域の形状に基づいて、前記眼底の奥行方向に凸な曲線を求める曲線演算手段と、
前記特定された層領域と前記求められた曲線とに基づいて、前記奥行方向の逆方向に向かって前記層領域が突出する突出領域を特定する突出領域特定手段と、
該特定された突出領域の形態を表す形態情報を生成する形態情報生成手段と、
を備えることを特徴とする眼底解析装置。
【請求項２】
前記曲線演算手段は、前記層領域特定手段により特定された層領域の形状に基づいて、前記層領域中の複数の特徴部位を特定する特徴部位特定手段を含み、該特定された複数の特徴部位に基づいて前記曲線を求める、
ことを特徴とする請求項１に記載の眼底解析装置。
【請求項３】
前記特徴部位特定手段は、前記層領域の形状に基づいて前記奥行方向における前記層領域中の最深部位を特定して前記特徴部位とし、前記最深部位を通過しかつ前記層領域に接する直線を求め、前記層領域と前記直線との接点を前記特徴部位とする、
ことを特徴とする請求項２に記載の眼底解析装置。
【請求項４】
前記特徴部位特定手段は、前記最深部位を通過する直線を前記最深部位を中心として回転させていくことにより前記接点を順次に特定する、
ことを特徴とする請求項３に記載の眼底解析装置。
【請求項５】
前記特徴部位特定手段は、前記最深部位を通過する直線を前記最深部位を中心として回転させて接点を特定し、該特定された接点を通過する直線を当該接点を中心として回転させて更なる接点を特定する、
ことを特徴とする請求項３に記載の眼底解析装置。
【請求項６】
前記曲線演算手段は、前記特定された複数の特徴部位に基づく二次曲線を前記曲線として求める、
ことを特徴とする請求項２〜請求項５のいずれか一項に記載の眼底解析装置。
【請求項７】
前記曲線演算手段は、前記複数の特徴部位を結んだ折れ線との差が最小になる二次曲線を最小二乗法によって求める、
ことを特徴とする請求項６に記載の眼底解析装置。
【請求項８】
前記曲線演算手段は、予め前記断層像に設定された座標系における前記複数の特徴部位のそれぞれの座標値を二次曲線の方程式に代入して演算することで前記方程式の係数を求める、
ことを特徴とする請求項６に記載の眼底解析装置。
【請求項９】
前記突出領域特定手段は、前記層領域と前記曲線との間の前記奥行方向における距離が所定閾値以上となる画像領域を前記突出領域として特定する、
ことを特徴とする請求項１〜請求項８のいずれか一項に記載の眼底解析装置。
【請求項１０】
前記突出領域特定手段は、前記曲線上の各点と前記層領域との間の前記奥行方向における距離を算出し、該算出された距離が前記所定閾値以上であるか判断し、前記所定閾値以上であると判断された前記曲線上の点の集合と前記層領域とに挟まれた画像領域を特定して前記突出領域とする、
ことを特徴とする請求項９に記載の眼底解析装置。
【請求項１１】
前記形態情報生成手段は、前記突出領域特定手段により前記各断層像について特定された突出領域に基づいて、前記奥行方向に直交する平面における突出領域の分布状態を表す分布画像を形成する分布画像形成手段と、該形成された分布画像の画素の画素値に基づいて当該分布画像における突出領域の連結成分を特定する連結成分特定手段とを含み、
該特定された連結成分に基づいて前記形態情報を生成する、
ことを特徴とする請求項１〜請求項１０のいずれか一項に記載の眼底解析装置。
【請求項１２】
前記形態情報生成手段は、前記連結成分特定手段により特定された連結成分の個数をカウントするカウント手段を含み、前記カウントにより得られた個数に基づいて前記形態情報を生成する、
ことを特徴とする請求項１１に記載の眼底解析装置。
【請求項１３】
前記形態情報生成手段は、前記連結成分特定手段により特定された各連結成分のサイズを算出するサイズ算出手段を含み、該算出されたサイズの分布を表すサイズ分布情報を生成して前記形態情報とする、
ことを特徴とする請求項１１に記載の眼底解析装置。
【請求項１４】
前記サイズ算出手段は、前記サイズとして前記各連結成分の面積を算出する、
ことを特徴とする請求項１３に記載の眼底解析装置。
【請求項１５】
前記サイズ算出手段は、前記サイズとして前記各連結成分の径を算出する、
ことを特徴とする請求項１３に記載の眼底解析装置。
【請求項１６】
前記サイズ算出手段は、前記各連結成分の面積を算出し、該算出された面積を有する円の径を求めて当該連結成分の径とする、
ことを特徴とする請求項１５に記載の眼底解析装置。
【請求項１７】
前記サイズ算出手段は、前記サイズとして前記各連結成分の体積を算出する、
ことを特徴とする請求項１３に記載の眼底解析装置。
【請求項１８】
前記サイズ算出手段は、前記層領域と前記曲線との間の前記奥行方向における距離を、前記各連結成分にわたって積分することにより当該連結成分の体積を算出する、
ことを特徴とする請求項１７に記載の眼底解析装置。
【請求項１９】
前記記憶手段は、前記眼底の撮影画像を更に記憶し、
前記形態情報生成手段は、前記撮影画像と前記分布画像との合成画像を形成して前記形態情報とする、
ことを特徴とする請求項１１に記載の眼底解析装置。
【請求項２０】
光源からの光を信号光と参照光とに分割し、被検眼の眼底を経由した前記信号光と参照光路を経由した前記参照光とを重畳させて干渉光を生成して検出する光学系と、
前記干渉光の検出結果に基づいて前記眼底の複数の断層像を形成する画像形成手段と、
を更に備え、
前記記憶手段は、前記画像形成手段により形成された複数の断層像を記憶する、
ことを特徴とする請求項１〜請求項１９のいずれか一項に記載の眼底解析装置。
【請求項２１】
眼底の層構造を描写する３次元画像を記憶する記憶手段と、
前記記憶された３次元画像の画素の画素値に基づいて、網膜色素上皮層に相当する前記３次元画像中の層領域を特定する層領域特定手段と、
前記特定された層領域の形状に基づいて、前記眼底の奥行方向に凸な曲面を求める曲面演算手段と、
前記特定された層領域と前記求められた曲面とに基づいて、前記奥行方向の逆方向に向かって前記層領域が突出する突出領域を特定する突出領域特定手段と、
該特定された突出領域の形態を表す形態情報を生成する形態情報生成手段と、
を備えることを特徴とする眼底解析装置。
【請求項２２】
光源からの光を信号光と参照光とに分割し、被検眼の眼底を経由した前記信号光と参照光路を経由した前記参照光とを重畳させて干渉光を生成して検出する光学系と、
前記干渉光の検出結果に基づいて前記眼底の複数の断層像を形成する画像形成手段と、
該形成された前記複数の断層像に基づいて３次元画像を形成する３次元画像形成手段と、
を更に備え、
前記記憶手段は、前記３次元画像形成手段により形成された３次元画像を記憶する、
ことを特徴とする請求項２１に記載の眼底解析装置。
【請求項２３】
前記光学系は、前記眼底に対する前記信号光の照射位置を複数の走査線に沿って順次に走査する走査手段を含み、
前記画像形成手段は、前記光学系による前記干渉光の検出結果に基づいて、前記複数の走査線のそれぞれに沿った断層像を形成する、
ことを特徴とする請求項２０又は請求項２２に記載の眼底解析装置。
【請求項２４】
眼底の層構造をそれぞれ描写する複数の断層像を記憶する記憶手段を有するコンピュータに、
前記記憶された各断層像の画素の画素値に基づいて、網膜色素上皮層に相当する当該断層像中の層領域を特定させ、
前記特定された層領域の形状に基づいて、前記眼底の奥行方向に凸な曲線を求めさせ、
前記特定された層領域と前記求められた曲線とに基づいて、前記奥行方向の逆方向に向かって前記層領域が突出する突出領域を特定させ、
該特定された突出領域の形態を表す形態情報を生成させる、
ことを特徴とする眼底解析プログラム。
【請求項２５】
眼底の層構造を描写する３次元画像を記憶する記憶手段を有するコンピュータに、
前記記憶された３次元画像の画素の画素値に基づいて、網膜色素上皮層に相当する前記３次元画像中の層領域を特定させ、
前記特定された層領域の形状に基づいて、前記眼底の奥行方向に凸な曲面を求めさせ、
前記特定された層領域と前記求められた曲面とに基づいて、前記奥行方向の逆方向に向かって前記層領域が突出する突出領域を特定させ、
該特定された突出領域の形態を表す形態情報を生成させる、
ことを特徴とする眼底解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、光コヒーレンストモグラフィ（Ｏｐｔｉｃａｌ  Ｃｏｈｅｒｅｎｃｅ  Ｔｏｍｏｇｒａｐｈｙ：ＯＣＴ）を用いて形成された眼底の画像を解析する技術に関する。
【背景技術】
【０００２】
近年、レーザ光源等からの光ビームを用いて被測定物体の表面形態や内部形態を表す画像を形成するＯＣＴが注目を集めている。ＯＣＴは、Ｘ線ＣＴのような人体に対する侵襲性を持たないことから、特に医療分野や生物学分野における応用の展開が期待されている。たとえば眼科分野においては、眼底や角膜等の画像を形成する装置が実用化段階に入っている。
【０００３】
特許文献１にはＯＣＴを適用した装置が開示されている。この装置は、測定腕が回転式転向鏡（ガルバノミラー）により物体を走査し、参照腕に参照ミラーが設置されており、その出口に計測腕及び参照腕からの光束の干渉光の強度を分光器で分析する干渉器が設けられている。更に、参照腕は、参照光光束位相を不連続な値で段階的に変えるように構成されている。
【０００４】
特許文献１の装置は、いわゆる「フーリエドメインＯＣＴ（Ｆｏｕｒｉｅｒ  Ｄｏｍａｉｎ  ＯＣＴ）」の手法を用いるものである。すなわち、被測定物体に対して低コヒーレンス光のビームを照射し、その反射光と参照光とを重ね合わせて干渉光を生成し、この干渉光のスペクトル強度分布を取得してフーリエ変換を施すことにより被測定物体の深度方向（ｚ方向）の形態を画像化するものである。なお、このタイプの手法は、特にスペクトラルドメイン（Ｓｐｅｃｔｒａｌ  Ｄｏｍａｉｎ）とも呼ばれる。
【０００５】
更に、特許文献１に記載の装置は、光ビーム（信号光）を走査するガルバノミラーを備え、それにより被測定物体の所望の測定対象領域の画像を形成するようになっている。この装置においては、ｚ方向に直交する１方向（ｘ方向）にのみ光ビームを走査するように構成されているので、この装置により形成される画像は、光ビームの走査方向（ｘ方向）に沿った深度方向（ｚ方向）の２次元断層像となる。
【０００６】
特許文献２には、信号光を水平方向（ｘ方向）及び垂直方向（ｙ方向）に走査（スキャン）することにより水平方向の２次元断層像を複数形成し、これら複数の断層像に基づいて測定範囲の３次元の断層情報を取得して画像化する技術が開示されている。この３次元画像化としては、たとえば、複数の断層像を垂直方向に並べて表示させる方法や（スタックデータなどと呼ばれる）、複数の断層像にレンダリング処理を施して３次元画像を形成する方法などが考えられる。
【０００７】
特許文献３、４には、他のタイプのＯＣＴ装置が開示されている。特許文献３には、被測定物体に照射される光の波長を走査し、各波長の光の反射光と参照光とを重ね合わせて得られる干渉光に基づいてスペクトル強度分布を取得し、それに対してフーリエ変換を施すことにより被測定物体の形態を画像化するＯＣＴ装置が記載されている。このようなＯＣＴ装置は、スウェプトソース（Ｓｗｅｐｔ  Ｓｏｕｒｃｅ）タイプなどと呼ばれる。スウェプトソースタイプはフーリエドメインタイプの一種である。
【０００８】
また、特許文献４には、所定のビーム径を有する光を被測定物体に照射し、その反射光と参照光とを重ね合わせて得られる干渉光の成分を解析することにより、光の進行方向に直交する断面における被測定物体の画像を形成するＯＣＴ装置が記載されている。このようなＯＣＴ装置は、フルフィールド（ｆｕｌｌ−ｆｉｅｌｄ）タイプ、或いはインファス（ｅｎ−ｆａｃｅ）タイプなどと呼ばれる。
【０００９】
特許文献５には、ＯＣＴを眼科分野に適用した構成が開示されている。この文献に記載された装置は、眼底を撮影して眼底像を形成する機能と、ＯＣＴを用いて眼底を計測して断層像や３次元画像（まとめてＯＣＴ画像と呼ぶことがある）を形成する機能とを備える。更に、この装置は、断層像を解析して、眼底を構成する層組織に相当する画像領域を特定する。特定対象の層組織としては、内境界膜、神経繊維層、神経節細胞層、内網状層、内顆粒層、外網状層、外顆粒層、外境界膜、視細胞層、網膜色素上皮層などがある。なお、眼底は複数の層組織が重なって構成されるので、層組織に相当する画像領域を求めることと、隣接する層組織の境界位置に相当する画像領域を求めることは同義である。なお、ＯＣＴが応用される以前から、眼底を観察するための装置としては眼底カメラが広く使用されている（たとえば特許文献６を参照）。
【００１０】
ＯＣＴを用いた装置は、高精細の画像を取得できる点、更には断層像や３次元画像を取得できる点などにおいて、眼底カメラ等に対して優位性を持つ。
【００１１】
ところで近年、加齢黄斑変性症と呼ばれる眼疾患が注目を集めている。加齢黄斑変性症は、網膜の黄斑部の機能が老化によって低下するために起こる疾患であり、視界がゆがむ、視力が低下する、視野が部分的に見えにくくなる、周囲は正常に見えるが見ようとするものが見えない、などといった症状を引き起こす。
【００１２】
加齢黄斑変性症（滲出型）は次のようなメカニズムで発症すると考えられている。正常な網膜細胞は新陳代謝を繰り返している。新陳代謝で生じる老廃物は正常であれば網膜色素上皮内で消化されて消えてしまう。しかし加齢により網膜色素上皮の働きが低下すると、未消化の老廃物がブルッフ膜と網膜色素上皮層との間に溜まる。この状態の眼底を撮影すると、老廃物はドルーゼンと呼ばれる白い塊として認識される。老廃物が蓄積すると弱い炎症反応が起こる。そうすると特定の化学物質（ケミカルメディエータ）が産出されて炎症の治癒を促す。しかし、ケミカルメディエータには血管の発生を促す因子が存在し、それにより脈絡膜から新たな血管（新生血管）が生えてくる。新生血管がブルッフ膜を突き破って網膜色素上皮層の下または上まで侵入して増殖すると、血液や血液成分の滲出が激しくなって黄斑の機能低下が顕著になる。
【００１３】
加齢黄斑変性症の診断ではドルーゼンの存在や分布が重要となる。従来、ドルーゼンの状態を把握するためには主として眼底像（眼底カメラによる撮影画像）が用いられていた（たとえば特許文献７を参照）。
【先行技術文献】
【特許文献】
【００１４】
【特許文献１】特開平１１−３２５８４９号公報
【特許文献２】特開２００２−１３９４２１号公報
【特許文献３】特開２００７−２４６７７号公報
【特許文献４】特開２００６−１５３８３８号公報
【特許文献５】特開２００８−７３０９９公報
【特許文献６】特開平９−２７６２３２号公報
【特許文献７】特開２００８−２９５８０４号公報
【発明の概要】
【発明が解決しようとする課題】
【００１５】
前述のようにドルーゼンは白い塊として眼底像に描写される。比較的大きなドルーゼンについては目視観察や画像解析により眼底像中から見つけることは可能であるが、比較的小さなドルーゼンについては周囲の正常部位との色や明るさの差が小さいために、これを発見することは困難であった。
【００１６】
加齢黄斑変性症では早期治療が重要であるが、小さなドルーゼンを発見できないことは発症の早期発見の障害となる。また、小さなドルーゼンを発見できないと、実際にドルーゼンがどのように分布しているか正確に把握することができない。
【００１７】
また、眼底にレーザ治療の跡や網膜剥離などが存在すると、眼底像ではこれらとドルーゼンとの区別が難しい場合があった。
【００１８】
この発明は、以上のような問題を解決するためになされたもので、その目的は、ドルーゼンを効果的に検出することが可能な眼底解析装置及び眼底解析プログラムを提供することにある。
【課題を解決するための手段】
【００１９】
上記目的を達成するために、請求項１に記載の発明は、眼底の層構造をそれぞれ描写する複数の断層像を記憶する記憶手段と、前記記憶された各断層像の画素の画素値に基づいて、網膜色素上皮層に相当する当該断層像中の層領域を特定する層領域特定手段と、前記特定された層領域の形状に基づいて、前記眼底の奥行方向に凸な曲線を求める曲線演算手段と、前記特定された層領域と前記求められた曲線とに基づいて、前記奥行方向の逆方向に向かって前記層領域が突出する突出領域を特定する突出領域特定手段と、該特定された突出領域の形態を表す形態情報を生成する形態情報生成手段と、を備えることを特徴とする眼底解析装置である。
また、請求項２に記載の発明は、請求項１に記載の眼底解析装置であって、前記曲線演算手段は、前記層領域特定手段により特定された層領域の形状に基づいて、前記層領域中の複数の特徴部位を特定する特徴部位特定手段を含み、該特定された複数の特徴部位に基づいて前記曲線を求める、ことを特徴とする。
また、請求項３に記載の発明は、請求項２に記載の眼底解析装置であって、前記特徴部位特定手段は、前記層領域の形状に基づいて前記奥行方向における前記層領域中の最深部位を特定して前記特徴部位とし、前記最深部位を通過しかつ前記層領域に接する直線を求め、前記層領域と前記直線との接点を前記特徴部位とする、ことを特徴とする。
また、請求項４に記載の発明は、請求項３に記載の眼底解析装置であって、前記特徴部位特定手段は、前記最深部位を通過する直線を前記最深部位を中心として回転させていくことにより前記接点を順次に特定する、ことを特徴とする。
また、請求項５に記載の発明は、請求項３に記載の眼底解析装置であって、前記特徴部位特定手段は、前記最深部位を通過する直線を前記最深部位を中心として回転させて接点を特定し、該特定された接点を通過する直線を当該接点を中心として回転させて更なる接点を特定する、ことを特徴とする。
また、請求項６に記載の発明は、請求項２〜請求項５のいずれか一項に記載の眼底解析装置であって、前記曲線演算手段は、前記特定された複数の特徴部位に基づく二次曲線を前記曲線として求める、ことを特徴とする。
また、請求項７に記載の発明は、請求項６に記載の眼底解析装置であって、前記曲線演算手段は、前記複数の特徴部位を結んだ折れ線との差が最小になる二次曲線を最小二乗法によって求める、ことを特徴とする。
また、請求項８に記載の発明は、請求項６に記載の眼底解析装置であって、前記曲線演算手段は、予め前記断層像に設定された座標系における前記複数の特徴部位のそれぞれの座標値を二次曲線の方程式に代入して演算することで前記方程式の係数を求める、ことを特徴とする。
また、請求項９に記載の発明は、請求項１〜請求項８のいずれか一項に記載の眼底解析装置であって、前記突出領域特定手段は、前記層領域と前記曲線との間の前記奥行方向における距離が所定閾値以上となる画像領域を前記突出領域として特定する、ことを特徴とする。
また、請求項１０に記載の発明は、請求項９に記載の眼底解析装置であって、前記突出領域特定手段は、前記曲線上の各点と前記層領域との間の前記奥行方向における距離を算出し、該算出された距離が前記所定閾値以上であるか判断し、前記所定閾値以上であると判断された前記曲線上の点の集合と前記層領域とに挟まれた画像領域を特定して前記突出領域とする、ことを特徴とする。
また、請求項１１に記載の発明は、請求項１〜請求項１０のいずれか一項に記載の眼底解析装置であって、前記形態情報生成手段は、前記突出領域特定手段により前記各断層像について特定された突出領域に基づいて、前記奥行方向に直交する平面における突出領域の分布状態を表す分布画像を形成する分布画像形成手段と、該形成された分布画像の画素の画素値に基づいて当該分布画像における突出領域の連結成分を特定する連結成分特定手段とを含み、該特定された連結成分に基づいて前記形態情報を生成する、ことを特徴とする。
また、請求項１２に記載の発明は、請求項１１に記載の眼底解析装置であって、前記形態情報生成手段は、前記連結成分特定手段により特定された連結成分の個数をカウントするカウント手段を含み、前記カウントにより得られた個数に基づいて前記形態情報を生成する、ことを特徴とする。
また、請求項１３に記載の発明は、請求項１１に記載の眼底解析装置であって、前記形態情報生成手段は、前記連結成分特定手段により特定された各連結成分のサイズを算出するサイズ算出手段を含み、該算出されたサイズの分布を表すサイズ分布情報を生成して前記形態情報とする、ことを特徴とする。
また、請求項１４に記載の発明は、請求項１３に記載の眼底解析装置であって、前記サイズ算出手段は、前記サイズとして前記各連結成分の面積を算出する、ことを特徴とする。
また、請求項１５に記載の発明は、請求項１３に記載の眼底解析装置であって、前記サイズ算出手段は、前記サイズとして前記各連結成分の径を算出する、ことを特徴とする。
また、請求項１６に記載の発明は、請求項１５に記載の眼底解析装置であって、前記サイズ算出手段は、前記各連結成分の面積を算出し、該算出された面積を有する円の径を求めて当該連結成分の径とする、ことを特徴とする。
また、請求項１７に記載の発明は、請求項１３に記載の眼底解析装置であって、前記サイズ算出手段は、前記サイズとして前記各連結成分の体積を算出する、ことを特徴とする。
また、請求項１８に記載の発明は、請求項１７に記載の眼底解析装置であって、前記サイズ算出手段は、前記層領域と前記曲線との間の前記奥行方向における距離を、前記各連結成分にわたって積分することにより当該連結成分の体積を算出する、ことを特徴とする。
また、請求項１９に記載の発明は、請求項１１に記載の眼底解析装置であって、前記記憶手段は、前記眼底の撮影画像を更に記憶し、前記形態情報生成手段は、前記撮影画像と前記分布画像との合成画像を形成して前記形態情報とする、ことを特徴とする。
また、請求項２０に記載の発明は、請求項１〜請求項１９のいずれか一項に記載の眼底解析装置であって、光源からの光を信号光と参照光とに分割し、被検眼の眼底を経由した前記信号光と参照光路を経由した前記参照光とを重畳させて干渉光を生成して検出する光学系と、前記干渉光の検出結果に基づいて前記眼底の複数の断層像を形成する画像形成手段と、を更に備え、前記記憶手段は、前記画像形成手段により形成された複数の断層像を記憶する、ことを特徴とする。
また、請求項２１に記載の発明は、眼底の層構造を描写する３次元画像を記憶する記憶手段と、前記記憶された３次元画像の画素の画素値に基づいて、網膜色素上皮層に相当する前記３次元画像中の層領域を特定する層領域特定手段と、前記特定された層領域の形状に基づいて、前記眼底の奥行方向に凸な曲面を求める曲面演算手段と、前記特定された層領域と前記求められた曲面とに基づいて、前記奥行方向の逆方向に向かって前記層領域が突出する突出領域を特定する突出領域特定手段と、該特定された突出領域の形態を表す形態情報を生成する形態情報生成手段と、を備えることを特徴とする眼底解析装置である。
また、請求項２２に記載の発明は、請求項２１に記載の眼底解析装置であって、光源からの光を信号光と参照光とに分割し、被検眼の眼底を経由した前記信号光と参照光路を経由した前記参照光とを重畳させて干渉光を生成して検出する光学系と、前記干渉光の検出結果に基づいて前記眼底の複数の断層像を形成する画像形成手段と、該形成された前記複数の断層像に基づいて３次元画像を形成する３次元画像形成手段と、を更に備え、前記記憶手段は、前記３次元画像形成手段により形成された３次元画像を記憶する、ことを特徴とする。
また、請求項２３に記載の発明は、請求項２０又は請求項２２に記載の眼底解析装置であって、前記光学系は、前記眼底に対する前記信号光の照射位置を複数の走査線に沿って順次に走査する走査手段を含み、前記画像形成手段は、前記光学系による前記干渉光の検出結果に基づいて、前記複数の走査線のそれぞれに沿った断層像を形成する、ことを特徴とする。
また、請求項２４に記載の発明は、眼底の層構造をそれぞれ描写する複数の断層像を記憶する記憶手段を有するコンピュータに、前記記憶された各断層像の画素の画素値に基づいて、網膜色素上皮層に相当する当該断層像中の層領域を特定させ、前記特定された層領域の形状に基づいて、前記眼底の奥行方向に凸な曲線を求めさせ、前記特定された層領域と前記求められた曲線とに基づいて、前記奥行方向の逆方向に向かって前記層領域が突出する突出領域を特定させ、該特定された突出領域の形態を表す形態情報を生成させる、ことを特徴とする眼底解析プログラムである。
また、請求項２５に記載の発明は、眼底の層構造を描写する３次元画像を記憶する記憶手段を有するコンピュータに、前記記憶された３次元画像の画素の画素値に基づいて、網膜色素上皮層に相当する前記３次元画像中の層領域を特定させ、前記特定された層領域の形状に基づいて、前記眼底の奥行方向に凸な曲面を求めさせ、前記特定された層領域と前記求められた曲面とに基づいて、前記奥行方向の逆方向に向かって前記層領域が突出する突出領域を特定させ、該特定された突出領域の形態を表す形態情報を生成させる、ことを特徴とする眼底解析プログラムである。
【発明の効果】
【００２０】
この発明に係る眼底解析装置及び眼底解析プログラムの第１の態様によれば、眼底の層構造をそれぞれ描写する複数の断層像のそれぞれについて、その断層像の画素の画素値に基づいて網膜色素上皮層に相当する層領域を特定し、この層領域の形状に基づいて眼底の奥行方向に凸な曲線を求め、この層領域とこの曲線とに基づいてこの層領域の突出領域を特定し、この突出領域の形態を表す形態情報を生成するように構成されているので、ドルーゼンの可能性がある網膜色素上皮層の突出部分をそれぞれ個別に把握することができ、それにより小さなドルーゼンであっても検出することが可能となる。また、断層像に基づいてドルーゼンに特徴的な突出領域を特定するように構成されているので、ドルーゼンを効果的に検出することが可能である。
【００２１】
この発明に係る眼底解析装置及び眼底解析プログラムの第２の態様によれば、眼底の層構造を描写する３次元画像の画素の画素値に基づいて網膜色素上皮層に相当する層領域を特定し、この層領域の形状に基づいて眼底の奥行方向に凸な曲面を求め、この層領域とこの曲面とに基づいてこの層領域の突出領域を特定し、この突出領域の形態を表す形態情報を生成するように構成されているので、ドルーゼンの可能性がある網膜色素上皮層の突出部分をそれぞれ個別に把握することができ、それにより小さなドルーゼンであっても検出することが可能となる。また、３次元画像に基づいてドルーゼンに特徴的な突出領域を特定するように構成されているので、ドルーゼンを効果的に検出することが可能である。
【図面の簡単な説明】
【００２２】
【図１】この発明に係る眼底解析装置の実施形態の構成の一例を表す概略図である。
【図２】この発明に係る眼底解析装置の実施形態の構成の一例を表す概略図である。
【図３】この発明に係る眼底解析装置の実施形態の構成の一例を表す概略ブロック図である。
【図４】この発明に係る眼底解析装置の実施形態の構成の一例を表す概略ブロック図である。
【図５】この発明に係る眼底解析装置の実施形態が実行する処理の一例を説明するための概略図である。
【図６】この発明に係る眼底解析装置の実施形態が実行する処理の一例を説明するための概略図である。
【図７】この発明に係る眼底解析装置の実施形態が実行する処理の一例を説明するための概略図である。
【図８】この発明に係る眼底解析装置の実施形態が実行する処理の一例を説明するための概略図である。
【図９】この発明に係る眼底解析装置の実施形態が実行する処理の一例を説明するための概略図である。
【図１０】この発明に係る眼底解析装置の実施形態の動作の一例を表すフローチャートである。
【図１１】この発明に係る眼底解析装置の実施形態が表示する解析結果表示画面の一例を表す概略図である。
【発明を実施するための形態】
【００２３】
この発明に係る眼底解析装置及び眼底解析プログラムの実施形態の一例について、図面を参照しながら詳細に説明する。
【００２４】
この発明に係る眼底解析装置は、眼底のＯＣＴ画像（断層像、３次元画像）を解析するコンピュータであってもよいし、光コヒーレンストモグラフィを用いて眼底の断層像を形成することが可能なＯＣＴ装置であってもよい。後者のＯＣＴ装置には前者のコンピュータが含まれる。よって、以下においては後者のＯＣＴ装置について特に詳しく説明する。
【００２５】
ＯＣＴ装置としての眼底解析装置は、被検眼の眼底のＯＣＴ画像を形成できるものであればどのようなタイプであってもよい。以下の実施形態では、フーリエドメインタイプについて特に詳しく説明する。なお、この発明の特徴の中心は眼底のＯＣＴ画像を解析する処理にあるので、スウェプトソースタイプやインファスタイプなどの他のタイプのＯＣＴ装置であっても同様に構成することが可能である。ＯＣＴ画像を形成するための計測動作をＯＣＴ計測と呼ぶことがある。
【００２６】
以下に説明する眼底解析装置は、特許文献５に開示された装置と同様に、眼底のＯＣＴ画像及び眼底像の双方を取得可能とされる。
【００２７】
［構成］
図１及び図２に示すように、眼底解析装置１は、眼底カメラユニット２、ＯＣＴユニット１００及び演算制御ユニット２００を含んで構成される。眼底カメラユニット２は、従来の眼底カメラとほぼ同様の光学系を有する。ＯＣＴユニット１００には、眼底のＯＣＴ画像を取得するための光学系が設けられている。演算制御ユニット２００は、各種の演算処理や制御処理等を実行するコンピュータを具備している。
【００２８】
〔眼底カメラユニット〕
図１に示す眼底カメラユニット２には、被検眼Ｅの眼底Ｅｆの表面形態を表す２次元画像（眼底像）を形成するための光学系が設けられている。眼底像には、観察画像や撮影画像などが含まれる。観察画像は、たとえば、近赤外光を用いて所定のフレームレートで形成されるモノクロの動画像である。撮影画像は、たとえば、可視光をフラッシュ発光して得られるカラー画像である。なお、眼底カメラユニット２は、これら以外の画像、たとえばフルオレセイン蛍光画像やインドシアニングリーン蛍光画像や自発蛍光画像などを取得可能に構成されていてもよい。
【００２９】
眼底カメラユニット２には、従来の眼底カメラと同様に、被検者の顔が動かないように支えるための顎受けや額当てが設けられている。更に、眼底カメラユニット２には、従来の眼底カメラと同様に照明光学系１０と撮影光学系３０が設けられている。照明光学系１０は眼底Ｅｆに照明光を照射する。撮影光学系３０は、この照明光の眼底反射光を撮像装置（ＣＣＤイメージセンサ３５、３８）に導く。また、撮影光学系３０は、ＯＣＴユニット１００からの信号光を眼底Ｅｆに導くとともに、眼底Ｅｆを経由した信号光をＯＣＴユニット１００に導く。
【００３０】
照明光学系１０の観察光源１１は、たとえばハロゲンランプにより構成される。観察光源１１から出力された光（観察照明光）は、曲面状の反射面を有する反射ミラー１２により反射され、集光レンズ１３を経由し、可視カットフィルタ１４を透過して近赤外光となる。更に、観察照明光は、撮影光源１５の近傍にて一旦集束し、ミラー１６により反射され、リレーレンズ１７、１８、絞り１９及びリレーレンズ２０を経由する。そして、観察照明光は、孔開きミラー２１の周辺部（孔部の周囲の領域）にて反射され、対物レンズ２２を経由して眼底Ｅｆを照明する。
【００３１】
観察照明光の眼底反射光は、対物レンズ２２により屈折され、孔開きミラー２１の中心領域に形成された孔部を通過し、ダイクロイックミラー５５を透過し、合焦レンズ３１を経由し、ダイクロイックミラー３２により反射される。更に、この眼底反射光は、ハーフミラー４０を透過し、ダイクロイックミラー３３により反射され、集光レンズ３４によりＣＣＤイメージセンサ３５の受光面に結像される。ＣＣＤイメージセンサ３５は、たとえば所定のフレームレートで眼底反射光を検出する。表示装置３には、ＣＣＤイメージセンサ３５により検出された眼底反射光に基づく画像（観察画像）Ｋが表示される。
【００３２】
撮影光源１５は、たとえばキセノンランプにより構成される。撮影光源１５から出力された光（撮影照明光）は、観察照明光と同様の経路を通って眼底Ｅｆに照射される。撮影照明光の眼底反射光は、観察照明光のそれと同様の経路を通ってダイクロイックミラー３３まで導かれ、ダイクロイックミラー３３を透過し、ミラー３６により反射され、集光レンズ３７によりＣＣＤイメージセンサ３８の受光面に結像される。表示装置３には、ＣＣＤイメージセンサ３８により検出された眼底反射光に基づく画像（撮影画像）Ｈが表示される。なお、観察画像Ｋを表示する表示装置３と撮影画像Ｈを表示する表示装置３は、同一のものであってもよいし、異なるものであってもよい。
【００３３】
ＬＣＤ（Ｌｉｑｕｉｄ  Ｃｒｙｓｔａｌ  Ｄｉｓｐｌａｙ）３９は、固視標や視力測定用視標を表示する。固視標は被検眼Ｅを固視させるための視標であり、眼底撮影時やＯＣＴ計測時などに使用される。
【００３４】
ＬＣＤ３９から出力された光は、その一部がハーフミラー４０にて反射され、ダイクロイックミラー３２に反射され、合焦レンズ３１及びダイクロイックミラー５５を経由し、孔開きミラー２１の孔部を通過し、対物レンズ２２により屈折されて眼底Ｅｆに投影される。
【００３５】
ＬＣＤ３９の画面上における固視標の表示位置を変更することにより、被検眼Ｅの固視位置を変更できる。被検眼Ｅの固視位置としては、たとえば従来の眼底カメラと同様に、眼底Ｅｆの黄斑部を中心とする画像を取得するための位置や、視神経乳頭を中心とする画像を取得するための位置や、黄斑部と視神経乳頭との間の眼底中心を中心とする画像を取得するための位置などがある。
【００３６】
更に、眼底カメラユニット２には、従来の眼底カメラと同様にアライメント光学系５０とフォーカス光学系６０が設けられている。アライメント光学系５０は、被検眼Ｅに対する装置光学系の位置合わせ（アライメント）を行うための視標（アライメント視標）を生成する。フォーカス光学系６０は、眼底Ｅｆに対してフォーカス（ピント）を合わせるための視標（スプリット視標）を生成する。
【００３７】
アライメント光学系５０のＬＥＤ（Ｌｉｇｈｔ  Ｅｍｉｔｔｉｎｇ  Ｄｉｏｄｅ）５１から出力された光（アライメント光）は、絞り５２、５３及びリレーレンズ５４を経由してダイクロイックミラー５５により反射され、孔開きミラー２１の孔部を通過し、対物レンズ２２により被検眼Ｅの角膜に投影される。
【００３８】
アライメント光の角膜反射光は、対物レンズ２２及び上記孔部を経由し、その一部がダイクロイックミラー５５を透過し、合焦レンズ３１を通過し、ダイクロイックミラー３２により反射され、ハーフミラー４０を透過し、ダイクロイックミラー３３に反射され、集光レンズ３４によりＣＣＤイメージセンサ３５の受光面に投影される。ＣＣＤイメージセンサ３５による受光像（アライメント視標）は、観察画像Ｋとともに表示装置３に表示される。ユーザは、従来の眼底カメラと同様の操作を行ってアライメントを実施する。また、演算制御ユニット２００がアライメント視標の位置を解析して光学系を移動させることによりアライメントを行ってもよい。
【００３９】
フォーカス調整を行う際には、照明光学系１０の光路上に反射棒６７の反射面が斜設される。フォーカス光学系６０のＬＥＤ６１から出力された光（フォーカス光）は、リレーレンズ６２を通過し、スプリット視標板６３により二つの光束に分離され、二孔絞り６４を通過し、ミラー６５に反射され、集光レンズ６６により反射棒６７の反射面に一旦結像されて反射される。更に、フォーカス光は、リレーレンズ２０を経由し、孔開きミラー２１に反射され、対物レンズ２２により眼底Ｅｆに結像される。
【００４０】
フォーカス光の眼底反射光は、アライメント光の角膜反射光と同様の経路を通ってＣＣＤイメージセンサ３５により検出される。ＣＣＤイメージセンサ３５による受光像（スプリット視標）は、観察画像とともに表示装置３に表示される。演算制御ユニット２００は、従来と同様に、スプリット視標の位置を解析して合焦レンズ３１及びフォーカス光学系６０を移動させてピント合わせを行う。また、スプリット視標を視認しつつ手動でピント合わせを行ってもよい。
【００４１】
ダイクロイックミラー３２の後方には、ミラー４１、コリメータレンズ４２、及びガルバノミラー４３、４４を含む光路が設けられている。この光路はＯＣＴユニット１００に導かれている。
【００４２】
ガルバノミラー４４は、ＯＣＴユニット１００からの信号光ＬＳをｘ方向に走査する。ガルバノミラー４３は、信号光ＬＳをｙ方向に走査する。これら二つのガルバノミラー４３、４４により、信号光ＬＳをｘｙ平面上の任意の方向に走査することができる。
【００４３】
〔ＯＣＴユニット〕
ＯＣＴユニット１００には、眼底Ｅｆの断層像を取得するための光学系が設けられている（図２を参照）。この光学系は、従来のフーリエドメインタイプのＯＣＴ装置と同様の構成を有する。すなわち、この光学系は、低コヒーレンス光を参照光と信号光に分割し、眼底Ｅｆを経由した信号光と参照光路を経由した参照光とを干渉させて干渉光を生成し、この干渉光のスペクトル成分を検出するように構成されている。この検出結果（検出信号）は演算制御ユニット２００に送られる。
【００４４】
光源ユニット１０１は広帯域の低コヒーレンス光Ｌ０を出力する。低コヒーレンス光Ｌ０は、たとえば、近赤外領域の波長帯（約８００ｎｍ〜９００ｎｍ程度）を含み、数十マイクロメートル程度の時間的コヒーレンス長を有する。なお、人眼では視認できない波長帯、たとえば１０５０〜１０６０ｎｍ程度の中心波長を有する近赤外光を低コヒーレンス光Ｌ０として用いてもよい。
【００４５】
光源ユニット１０１は、スーパールミネセントダイオード（Ｓｕｐｅｒ  Ｌｕｍｉｎｅｓｃｅｎｔ  Ｄｉｏｄｅ：ＳＬＤ）や、ＬＥＤや、ＳＯＡ（Ｓｅｍｉｃｏｎｄｕｃｔｏｒ  Ｏｐｔｉｃａｌ  Ａｍｐｌｉｆｉｅｒ）等の光出力デバイスを含んで構成される。
【００４６】
光源ユニット１０１から出力された低コヒーレンス光Ｌ０は、光ファイバ１０２によりファイバカプラ１０３に導かれて信号光ＬＳと参照光ＬＲに分割される。なお、ファイバカプラ１０３は、光を分割する手段（スプリッタ；ｓｐｌｉｔｔｅｒ）、及び、光を合成する手段（カプラ；ｃｏｕｐｌｅｒ）の双方の作用を有するが、ここでは慣用的に「ファイバカプラ」と称する。
【００４７】
信号光ＬＳは、光ファイバ１０４により導光され、コリメータレンズユニット１０５により平行光束となる。更に、信号光ＬＳは、各ガルバノミラー４４、４３により反射され、コリメータレンズ４２により集光され、ミラー４１により反射され、ダイクロイックミラー３２を透過し、ＬＣＤ３９からの光と同じ経路を通って眼底Ｅｆに照射される。信号光ＬＳは、眼底Ｅｆにおいて散乱、反射される。この散乱光及び反射光をまとめて信号光ＬＳの眼底反射光と称することがある。信号光ＬＳの眼底反射光は、同じ経路を逆向きに進行してファイバカプラ１０３に導かれる。眼底Ｅｆに対する信号光ＬＳの照射方向（＋ｚ方向）は、眼底Ｅｆの表面（網膜表面）から深部に向かう方向であるから、これを奥行方向と呼ぶことにする。また、ｚ軸に沿う方向を深度方向と呼ぶことにする。
【００４８】
参照光ＬＲは、光ファイバ１０６により導光され、コリメータレンズユニット１０７により平行光束となる。更に、参照光ＬＲは、ミラー１０８、１０９、１１０により反射され、ＮＤ（Ｎｅｕｔｒａｌ  Ｄｅｎｓｉｔｙ）フィルタ１１１により減光され、ミラー１１２に反射され、コリメータレンズ１１３により参照ミラー１１４の反射面に結像される。参照ミラー１１４に反射された参照光ＬＲは、同じ経路を逆向きに進行してファイバカプラ１０３に導かれる。なお、分散補償用の光学素子（ペアプリズム等）や、偏光補正用の光学素子（波長板等）を参照光ＬＲの光路（参照光路）に設けてもよい。
【００４９】
ファイバカプラ１０３は、信号光ＬＳの眼底反射光と、参照ミラー１１４に反射された参照光ＬＲとを合波する。これにより生成された干渉光ＬＣは、光ファイバ１１５により導光されて出射端１１６から出射される。更に、干渉光ＬＣは、コリメータレンズ１１７により平行光束とされ、回折格子１１８により分光（スペクトル分解）され、集光レンズ１１９により集光されてＣＣＤイメージセンサ１２０の受光面に投影される。図２に示す回折格子１１８は透過型であるが、反射型の回折格子を用いてもよい。
【００５０】
ＣＣＤイメージセンサ１２０は、たとえばラインセンサであり、分光された干渉光ＬＣの各スペクトル成分を検出して電荷に変換する。ＣＣＤイメージセンサ１２０は、この電荷を蓄積して検出信号を生成する。更に、ＣＣＤイメージセンサ１２０は、この検出信号を演算制御ユニット２００に送る。
【００５１】
この実施形態ではマイケルソン型の干渉計を採用しているが、たとえばマッハツェンダー型など任意のタイプの干渉計を適宜に採用することが可能である。また、ＣＣＤイメージセンサに代えて、他の形態のイメージセンサ、たとえばＣＭＯＳ（Ｃｏｍｐｌｅｍｅｎｔａｒｙ  Ｍｅｔａｌ  Ｏｘｉｄｅ  Ｓｅｍｉｃｏｎｄｕｃｔｏｒ）イメージセンサなどを用いることが可能である。
【００５２】
〔演算制御ユニット〕
演算制御ユニット２００の構成について説明する。演算制御ユニット２００は、ＣＣＤイメージセンサ１２０から入力される検出信号を解析して眼底ＥｆのＯＣＴ画像を形成する。そのための演算処理は、従来のフーリエドメインタイプのＯＣＴ装置と同様である。更に、演算制御ユニット２００は、眼底ＥｆのＯＣＴ画像に対して後述の解析処理を施す。この解析処理がこの実施形態の特徴である。
【００５３】
また、演算制御ユニット２００は、眼底カメラユニット２、表示装置３及びＯＣＴユニット１００の各部を制御する。たとえば演算制御ユニット２００は、眼底Ｅｆの断層像Ｇ（図２を参照）等のＯＣＴ画像を表示装置３に表示させる。
【００５４】
また、眼底カメラユニット２の制御として、演算制御ユニット２００は、観察光源１１、撮影光源１５及びＬＥＤ５１、６１の動作制御、ＬＣＤ３９の動作制御、合焦レンズ３１の移動制御、反射棒６７の移動制御、フォーカス光学系６０の移動制御、各ガルバノミラー４３、４４の動作制御などを行う。
【００５５】
また、ＯＣＴユニット１００の制御として、演算制御ユニット２００は、光源ユニット１０１の動作制御、参照ミラー１１４及びコリメータレンズ１１３の移動制御、ＣＣＤイメージセンサ１２０の動作制御などを行う。
【００５６】
演算制御ユニット２００は、たとえば、従来のコンピュータと同様に、マイクロプロセッサ、ＲＡＭ、ＲＯＭ、ハードディスクドライブ、通信インターフェイスなどを含んで構成される。ハードディスクドライブ等の記憶装置には、眼底解析装置１を制御するためのコンピュータプログラムが記憶されている。演算制御ユニット２００は、ＣＣＤイメージセンサ１２０からの検出信号に基づいてＯＣＴ画像を形成する専用の回路基板を備えていてもよい。また、演算制御ユニット２００は、キーボードやマウス等の操作デバイス（入力デバイス）や、ＬＣＤ等の表示デバイスを備えていてもよい。
【００５７】
眼底カメラユニット２、表示装置３、ＯＣＴユニット１００及び演算制御ユニット２００は、一体的に（つまり単一の筺体内に）構成されていてもよいし、それぞれ別体として構成されていてもよい。
【００５８】
〔制御系〕
眼底解析装置１の制御系の構成について図３及び図４を参照しつつ説明する。
【００５９】
（制御部）
眼底解析装置１の制御系は、演算制御ユニット２００の制御部２１０を中心に構成される。制御部２１０は、たとえば、前述のマイクロプロセッサ、ＲＡＭ、ＲＯＭ、ハードディスクドライブ、通信インターフェイス等を含んで構成される。
【００６０】
制御部２１０には、主制御部２１１と記憶部２１２が設けられている。主制御部２１１は、前述の各種制御を行う。特に、主制御部２１１は、眼底カメラユニット２の走査駆動部７０及び合焦駆動部８０、更にＯＣＴユニット１００の光源ユニット１０１及び参照駆動部１３０を制御する。
【００６１】
走査駆動部７０は、たとえばサーボモータを含んで構成され、ガルバノミラー４３、４４の向きを各々独立に変更する。走査駆動部７０は、ガルバノミラー４３、４４とともに、この発明の「走査手段」の一例を構成する。
【００６２】
合焦駆動部８０は、たとえばパルスモータを含んで構成され、合焦レンズ３１を光軸方向に移動させる。それにより、眼底Ｅｆに向かう光の合焦位置が変更される。
【００６３】
参照駆動部１３０は、たとえばパルスモータを含んで構成され、参照光ＬＲの進行方向に沿って、コリメータレンズ１１３及び参照ミラー１１４を一体的に移動させる。
【００６４】
また、主制御部２１１は、記憶部２１２にデータを書き込む処理や、記憶部２１２からデータを読み出す処理を行う。
【００６５】
記憶部２１２は、各種のデータを記憶する。記憶部２１２に記憶されるデータとしては、たとえば、ＯＣＴ画像の画像データ、眼底像（観察画像Ｋ、撮影画像Ｈ）の画像データ、被検眼情報などがある。被検眼情報は、患者ＩＤや氏名などの被検者に関する情報や、左眼／右眼の識別情報などの被検眼に関する情報を含む。記憶部２１２は、この発明の「記憶手段」の一例である。
【００６６】
（画像形成部）
画像形成部２２０は、ＣＣＤイメージセンサ１２０からの検出信号に基づいて、眼底Ｅｆの断層像の画像データを形成する。この処理には、従来のフーリエドメインタイプの光コヒーレンストモグラフィと同様に、ノイズ除去（ノイズ低減）、フィルタ処理、ＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）などの処理が含まれている。
【００６７】
画像形成部２２０は、たとえば、前述の回路基板や通信インターフェイス等を含んで構成される。なお、この明細書では、「画像データ」と、それに基づいて呈示される「画像」とを同一視することがある。画像形成部２２０は、この発明の「画像形成手段」の一例である。
【００６８】
（画像処理部）
画像処理部２３０は、画像形成部２２０により形成された画像に対して各種の画像処理や解析処理を施す。たとえば、画像処理部２３０は、画像の輝度補正や分散補正等の各種補正処理を実行する。画像処理部２３０は、たとえば、前述のマイクロプロセッサ、ＲＡＭ、ＲＯＭ、ハードディスクドライブ、回路基板等を含んで構成される。画像処理部２３０には３次元画像形成部２３１と解析処理部２３２が設けられている。
【００６９】
（３次元画像形成部）
３次元画像形成部２３１は、画像形成部２２０により形成された断層像の間の画素を補間する補間処理等の公知の画像処理を実行して、眼底Ｅｆの３次元画像の画像データを形成する。３次元画像形成部２３１は、この発明の「３次元画像形成手段」の一例である。
【００７０】
なお、３次元画像の画像データとは、３次元座標系により画素の位置が定義された画像データを意味する。３次元画像の画像データとしては、３次元的に配列されたボクセルからなる画像データがある。この画像データは、ボリュームデータ或いはボクセルデータなどと呼ばれる。ボリュームデータに基づく画像を表示させる場合、画像処理部２３０は、このボリュームデータに対してレンダリング処理（ボリュームレンダリングやＭＩＰ（Ｍａｘｉｍｕｍ  Ｉｎｔｅｎｓｉｔｙ  Ｐｒｏｊｅｃｔｉｏｎ：最大値投影）など）を施して、特定の視線方向から見たときの擬似的な３次元画像の画像データを形成する。表示部２４０等の表示デバイスには、この擬似的な３次元画像が表示される。
【００７１】
また、３次元画像の画像データとして、複数の断層像のスタックデータを形成することも可能である。スタックデータは、複数の走査線に沿って得られた複数の断層像を、走査線の位置関係に基づいて３次元的に配列させることで得られる画像データである。すなわち、スタックデータは、元々個別の２次元座標系により定義されていた複数の断層像を、一つの３次元座標系により表現する（つまり一つの３次元空間に埋め込む）ことにより得られる画像データである。
【００７２】
画像処理部２３０は、３次元画像の画像データに基づいて、任意の断面における断層像を形成することができる。この処理は、たとえば、手動又は自動で指定された断面に対し、この断面上に位置する画素（ボクセル等）を特定し、特定された画素を２次元的に配列させて当該断面における眼底Ｅｆの形態を表す画像データを形成することにより実行される。このような処理により、元の断層像の断面（信号光ＬＳの走査線の位置）だけでなく、所望の断面における断層像を取得することが可能となる。
【００７３】
（解析処理部）
解析処理部２３２は、眼底ＥｆのＯＣＴ画像に基づいて、眼底Ｅｆの形態や機能状態を把握するための各種の解析処理を実行する。解析対象となるＯＣＴ画像としては、眼底Ｅｆの断層像や３次元画像がある。この実施形態では、特に眼底Ｅｆの黄斑及びその周辺部位のＯＣＴ画像が解析対象となる。解析処理部２３２の構成の一例を図４に示す。解析処理部２３２には、層領域特定部２３３、曲線演算部２３４、突出領域特定部２３５及び形態情報生成部２３６が設けられている。
【００７４】
（層領域特定部）
層領域特定部２３３は、解析対象のＯＣＴ画像の画素の画素値に基づいて、眼底Ｅｆの網膜色素上皮層に相当する当該ＯＣＴ画像中の画像領域（層領域）を特定する。
【００７５】
層領域特定部２３３は、たとえば特許文献５と同様の処理を実行して層領域を特定する。この処理について簡単に説明する。層領域特定部２３３は、まず、階調変換処理、画像強調処理、しきい値処理、コントラスト変換処理、二値化処理、エッジ検出処理、画像平均化処理、画像平滑化処理、フィルタ処理などの前処理をＯＣＴ画像に施して、当該ＯＣＴ画像における層領域を明瞭化する。
【００７６】
次に、層領域特定部２３３は、前処理が施されたＯＣＴ画像を構成する画素の画素値（たとえば輝度値）を、眼底Ｅｆの深度方向（ｚ軸方向）に沿って一列ずつ解析し、隣接する層の境界位置に相当する画素を特定する。このとき、深度方向にのみ広がりを有するフィルタ（たとえば微分フィルタ）を用いて層の境界位置に相当する画素を特定することができる。なお、深度方向及びそれに直交する方向の二方向に広がるエリアフィルタを用いて画素のエッジ検出を行うようにしてもよい。
【００７７】
層領域特定部２３３は、このような処理により眼底Ｅｆの幾つかの層に相当する画像領域を特定する。更に、層領域特定部２３３は、特定された幾つかの画像領域のうちから網膜色素上皮層に相当するものを特定する。この処理について説明する。ＯＣＴ画像中において、網膜表面側から数えて幾つ目の明るい層が網膜色素上皮層に相当しているかは、臨床的に取得された多数の眼底のＯＣＴ画像から既知である。よって、解析対象のＯＣＴ画像について、網膜表面をまず特定し、この網膜表面の側から明るい層の数をカウントし、所定のカウント値に相当する層が目的の層領域となる。
【００７８】
層領域の別の特定方法として、網膜表面から奥行方向に向かって網膜色素上皮層までの標準的な距離に基づいて、解析対象のＯＣＴ画像中の層領域を特定してもよい。また、ＯＣＴ画像中において、眼底Ｅｆの各層の明るさには差異があるので、当該差異を考慮して層領域を特定することができる。たとえば、明るく描写される層のうち網膜色素上皮層がＮ番目に明るく描写される場合、解析対象のＯＣＴ画像中に特定された層に相当する画像領域のうちＮ番目に明るいものを特定して層領域とすることができる。なお、層領域の特定方法はここに例示したものには限定されず、目的の層領域を特定可能であればその手法は問わない。
【００７９】
なお、解析対象のＯＣＴ画像は断層像や３次元画像である。断層像を解析する場合、層領域は、（層領域の厚みを無視すると）略曲線状の画像領域として特定される。一方、３次元画像を解析する場合、層領域は、（層領域の厚みを無視すると）略曲面状の画像料記として特定される。ここでは、断層像を解析する場合について詳しく説明する。３次元画像を解析する場合については変形例として後述する。
【００８０】
層領域特定部２３３は、特定された層領域の情報、たとえばＯＣＴ画像中における層領域の位置情報（座標値）を生成する。なお、ＯＣＴ画像から層領域を抽出するようにしてもよい。また、特定された層領域の形状を表す画像情報（たとえばワイヤモデル等）を生成してもよい。いずれにしても、層領域特定部２３３は、少なくとも、ＯＣＴ画像中における網膜色素上皮層に相当する層領域を特定するものであればよい。
【００８１】
（曲線演算部）
眼底Ｅｆの断層像における層領域（略曲線状の画像領域）の特定結果は曲線演算部２３４に入力される。曲線演算部２３４は、この層領域の形状に基づいて、眼底Ｅｆの奥行方向に凸な曲線（基準曲線と呼ぶ）を求める。この基準曲線は、層領域の形状に基づいて推定される、当該断面にドルーゼンが存在しないと仮定した場合の網膜色素上皮層の形状を表す。曲線演算部２３４は、この発明の「曲線演算手段」の一例である。
【００８２】
当該断面にドルーゼンが存在しない場合、層領域は奥行方向に凸な曲線として断層像宙に描写される。もし当該断面にドルーゼンが存在する場合、特定される層領域にはドルーゼンに相当する凹凸が描写される。基準曲線は、このような凹凸を無視した層領域の大局的な形状を表すものである。このような基準曲線を求めるために、曲線演算部２３４には特徴部位特定部２３４ａと曲線決定部２３４ｂが設けられている。特徴部位特定部２３４ａは、この発明の「特徴部位特定手段」の一例である。
【００８３】
（特徴部位特定部）
特徴部位特定部２３４ａは、特定された層領域中の画素の画素値に基づいて、この層領域の形状に基づく複数の特徴部位を特定する。この処理の例を図５に示す。なお、図５は断層像から特定された層領域の一部を表している。
【００８４】
まず、図５（Ａ）に示すように、特徴部位特定部２３４ａは、層領域３００の形状に基づいて、奥行方向（＋ｚ方向）における層領域３００中の最深部位Ｐ０を特定する。この処理は、たとえば、層領域３００中の各画素の座標値を参照し、ｚ座標値が最大の画素を特定して最深部位Ｐ０に設定することで実行できる。異なる手法として、解析対象の断層像において、奥行方向に直交する直線を＋ｚ方向から−ｚ方向に移動させていき、この直線に最初に接する層領域中の位置を最深部位に設定するようにしてもよい。このようにして特定された最深部位Ｐ０は層領域３００の特徴部位とされる。
【００８５】
次に、特徴部位特定部２３４ａは、最深部位Ｐ０を通過しかつ層領域３００に接する直線を求め、この直線と層領域３００との接点を特徴部位とする。この処理の具体例を説明する。図５（Ｂ）に示すように、最深部位Ｐ０を通過する直線Ｌを最深部位Ｐ０を中心として回転させていく。図５（Ｂ）は、直線Ｌを時計回りに回転させる場合を示している。
【００８６】
このように直線Ｌを回転させていくと、図５（Ｃ）に示すように、或る段階で直線Ｌが層領域３００に接する。このときの直線Ｌが、上記「最深部位Ｐ０を通過しかつ層領域３００に接する直線」に相当する。この接点Ｐ１は層領域３００の特徴部位とされる。なお、他の特徴部位は全て最深部位Ｐ０よりも−ｚ側に位置するので、たとえば最深部位Ｐ０を通過しｚ座標軸に直交する位置から直線Ｌを回転させれば十分である。また、直線Ｌを逆方向（反時計回り）に回転させていくことにより、最深部位Ｐ０に対して特徴部位Ｐ１と反対側に位置する特徴部位を特定することができる。
【００８７】
このような処理を繰り返すことにより層領域３００の複数の特徴部位Ｐｊ（ｊ＝０、１、２、・・・、Ｊ）が特定される。なお、上記処理の反復方法としては、たとえば次の２つが考えられる。なお、他の方法で複数の特徴部位Ｐｊを特定することも勿論可能である。
【００８８】
第１の反復方法として、常に最深部位Ｐ０を通過する直線Ｌを考慮し、直線Ｌと層領域３００との接点を順次に特定していくことができる。その場合、直線Ｌの回転に伴う第１番目の接点、第２番目の接点、第３番目の接点、・・・というように順次に接点（特徴部位）が特定されていく。
【００８９】
第２の反復方法として、直線Ｌの回転中心を順次に変更していくことができる。つまり、最初は最深部位Ｐ０を中心として直線Ｌを回転させて第１の接点を特定する。次に、この第１の接点を中心として同様に直線Ｌを回転させて第２の接点を特定する。続いて、この第２の接点を中心として同様に直線Ｌを回転させて第３の接点を特定する。このようにして複数の接点（特徴部位）が順次に特定されていく。
【００９０】
上記のようにして特定される特徴部位の個数は任意であるが、その個数が多いほど下記の処理の精度は向上する。一方、特徴部位の特定個数が増加するほど、処理に必要なリソースが増大する。
【００９１】
（曲線決定部）
曲線決定部２３４ｂは、特徴部位特定部２３４ａにより特定された複数の特徴部位Ｐｊに基づいて基準曲線を決定する。曲線決定部２３４ｂが実行する処理の一例を図６に示す。
【００９２】
まず、曲線決定部２３４ｂは、複数の特徴部位Ｐｊの隣接するもの同士を連結して、図６（Ａ）に示すような折れ線３１０を形成する。次に、曲線決定部２３４ｂは、この折れ線３１０との差が最小になる二次曲線を最小二乗法を用いて求める。この「差」は、たとえば、折れ線３１０と二次曲線とにより挟まれる領域の面積である。それにより、図６（Ｂ）に示すような、折れ線３１０を近似する二次曲線３２０が得られる。図６（Ｃ）は、この二次曲線３２０と層領域３００との位置関係を表している。この二次曲線３２０が基準曲線である。前述のように解析対象のＯＣＴ画像（断層像）は黄斑及びその周辺部位を計測して得られたものであるので、基準曲線は奥行方向に凸になる。
【００９３】
基準曲線を決定する方法の他の例として、二次曲線の係数パラメータを変更して曲線族を形成し、この曲線族に含まれる二次曲線のうち折れ線３１０の形状に最も適合するものをパターンマッチング等によって求めることもできる。
【００９４】
更なる例として、折れ線３１０を考慮せずに、特徴部位Ｐｊの座標値に基づく方法を説明する。まず、複数の特徴部位Ｐｊの座標値を二次曲線の方程式に代入して連立方程式を形成し、この連立方程式を解いて未定係数を算出することによって目的の二次曲線を求める。この例を適用する場合、少なくとも、二次曲線の未定係数の個数に応じた個数（たとえば６個）の特徴部位が必要となる。なお、基準曲線の特徴（たとえば下に凸であることやパリティなど）を考慮して未定係数を減らして特徴部位の必要数を減少させてもよい。
【００９５】
（突出領域特定部）
突出領域特定部２３５は、層領域特定部２３３により特定された層領域と、曲線演算部２３４により求められた基準曲線とに基づいて、奥行方向の逆方向（−ｚ方向）に向かって層領域が突出する画像領域（突出領域）を特定する。突出領域特定部２３５は、この発明の「突出領域特定手段」の一例である。
【００９６】
なお、基準曲線に対して−ｚ方向に突出する画像領域全体を突出領域として特定することもできるが、この実施形態では、網膜色素上皮層の自然な凹凸やノイズ等を回避するために、基準曲線から所定距離以上突出している部分のみを検出することにより、突出領域の特定精度の向上を図る。そのために、距離算出部２３５ａ、距離判断部２３５ｂ及び画像領域特定部２３５ｃが突出領域特定部２３５に設けられる。
【００９７】
（距離算出部）
距離算出部２３５ａは、基準曲線上の各点と層領域との間の奥行方向における距離を算出する。この処理は、基準曲線上の全ての点（画素）について実行する必要はなく、所定の画素間隔毎（たとえば５ピクセル毎）に行うようにしてもよい。距離の算出処理は、たとえば、基準曲線上の点（画素）と層領域上の対応点（画素）との間の画素数をカウントし、隣接する画素の間隔を単位距離と当該カウント結果とに基づいて行うことができる。また、画像の計測倍率と距離計測対象の画素間の画像中における距離とに基づいて求めてもよい。なお、距離算出部２３５ａにより演算される距離は、画像中の距離（ｘｙｚ座標系で定義される距離や、画素の間隔）を実空間の距離に換算したものであってもよいし、画像中の距離をそのまま使用してもよい。
【００９８】
（距離判断部）
距離判断部２３５ｂは、距離算出部２３５ａにより算出された各距離が所定閾値以上であるか判断する。この閾値は、たとえば多数の臨床例に基づいて事前に設定される。また、装置の計測精度などを考慮して閾値を設定することもできる。
【００９９】
距離判断部２３５ｂは、距離が所定閾値以上と判断された基準曲線上の点（画素）に識別情報（たとえばフラグやタグ）を付与する。
【０１００】
（画像領域特定部）
画像領域特定部２３５ｃは、距離判断部２３５ｂにより識別情報が付与された基準曲線上の画素の集合、つまり距離が所定閾値以上と判断された基準曲線上の点の集合と、層領域とに挟まれた画像領域を特定して目的の突出領域とする。
【０１０１】
このようにして特定される突出領域の例を図７に示す。突出領域３３０、３４０は、それぞれ、層領域３００と基準曲線３２０とに挟まれた画像領域であって、層領域３００と基準曲線３２０との距離が所定閾値以上であるその部分領域である。図７に示すように、基準曲線３２０に対して−ｚ方向に突出する層領域３００のうち裾野の部分（距離が所定閾値未満である部分）は突出領域３３０、３４０には含まれない。
【０１０２】
（形態情報生成部）
形態情報生成部２３６は、特定された突出領域の形態を表す情報（形態情報）を生成する。突出領域の形態としては、個数、サイズ、分布状態などがある。形態情報生成部２３６は、この発明の「形態情報生成手段」の一例である。形態情報生成部２３６には、分布画像形成部２３６ａ、連結成分特定部２３６ｂ、カウント部２３６ｃ及びサイズ算出部２３６ｄが設けられている。
【０１０３】
（分布画像形成部）
分布画像形成部２３６ａは、突出領域特定部２３５により複数の断層像のそれぞれについて突出領域が特定された場合に、奥行方向に直交するｘｙ平面における突出領域の分布状態を表す分布画像を形成する。分布画像形成部２３６ａは、この発明の「分布画像形成手段」の一例である。以下、分布画像形成部２３６ａが実行する処理内容を説明する。
【０１０４】
複数の断層像は、たとえば後述の３次元スキャンを実行して得られる。３次元スキャンは、たとえばｘ方向に沿った、かつｙ方向に配列された複数の直線状の走査線に沿って信号光ＬＳの照射位置を走査するスキャン形態である。３次元スキャンにより、各走査線に沿った断面における断層像が複数枚得られる。
【０１０５】
分布画像形成部２３６ａは、これら断層像のそれぞれについて特定された突出領域に基づいて、ｘｙ平面における突出領域の分布画像を形成する。各断層像において、各突出領域はｘ方向（走査線方向）に延びる画像領域である。また、複数の断層像はｙ方向に配列されている。複数の断層像をｙ方向に配列させると、各断層像中の突出領域が組み合わされて突出領域の２次元分布（ｘｙ平面における分布）が得られる。
【０１０６】
このとき、隣接する断層像における突出領域がｙ方向に隣接している場合には、これら突出領域間の画素を突出領域に設定してもよい。この処理は、特に、隣接する断層像の間隔（走査線の間隔）が十分狭い場合に有効である。
【０１０７】
分布画像形成部２３６ａは、たとえば、突出領域に相当する画素とそれ以外の画素との画素値を相違させて表現することにより分布画像を形成する。一例として、突出領域とそれ以外の領域とを二値で区別して表現して二値画像を形成し、これを分布画像とする。
【０１０８】
このようにして形成される分布画像の例を図８に示す。この分布画像４００は、信号光ＬＳの入射方向（−ｚ方向）から眼底Ｅｆをみたときの突出領域Ｔの分布状態を表している。分布画像４００は、複数の走査線Ｒｉ（ｉ＝１〜ｍ）を断面とする複数の断層像に基づいて形成される。
【０１０９】
（連結成分特定部）
連結成分特定部２３６ｂは、分布画像の画素の画素値に基づいて、この分布画像における突出領域の各連結成分を特定する。各連結成分は一つひとつの突出領域に相当する。連結成分特定部２３６ｂは、この発明の「連結成分特定手段」の一例である。
【０１１０】
連結成分特定部２３６ｂが実行する処理の例を図９を参照しつつ説明する。解析対象の分布画像は二値画像とし、突出領域以外の領域（背景領域）は画素値（輝度値）０で表現されているものとする。図９に示す各正方形は一つの画素を表している。
【０１１１】
画素「１」（丸文字の１）を基準にすると、図９（Ａ）に示すように、画素「１」の上下左右の位置（つまり画素「１」の＋ｘ側、−ｘ側、＋ｙ側、−ｙ側）に隣接し、かつ画素「１」と同じ画素値を有する画素「２」、「３」、「４」、「５」は、それぞれ画素「１」に対して連結であると判断される。
【０１１２】
これ以外の画素、つまり画素「１」の上下左右の位置に隣接していない画素や画素「１」と異なる画素値を有する画素は、画素「１」に対して連結ではないと判断される。たとえば図９（Ｂ）に示すように、画素「１」に対して斜めの位置に隣接する画素「６」、「７」は、たとえ画素「１」と同じ画素値を有するとしても非連結であると判断される。
【０１１３】
この要領で画素の連結状態を判断すると、図９（Ｃ）に示す９個の画素は連結となる。これら９個の画素のそれぞれに対して連結な画素が存在しない場合、この９個の画素の集合は一つの連結成分を形成する。
【０１１４】
連結成分特定部２３６ｂは、このような連結性の判断基準に基づいて分布画像における突出領域の各連結成分を求める。たとえば、連結成分特定部２３６ｂは、図８に示す分布画像４００に基づいて、突出領域Ｔの各連結成分Ｔｋ（ｋ＝１〜Ｋ）が特定される。
【０１１５】
なお、連結成分を特定する処理は二値画像以外の多値画像に対しても実行できる。たとえば、所定の閾値を設け、その閾値に対する大小比較によって突出領域の画素であるか背景領域の画素であるかを判別することにより各連結成分を特定できる。また、上記のように隣接画素に基づいて連結成分を特定する代わりに、たとえばエッジ検出処理を分布画像に施して突出領域のエッジを検出することにより各連結成分を特定することも可能である。
【０１１６】
（カウント部）
カウント部２３６ｃは、連結成分特定部２３６ｂにより特定された連結成分の個数をカウントする。この処理は、複数の連結成分（画素の集合）に対して所定の順序（たとえば分布画像の左上から右下に向かう順序）で番号１、２、・・・を順次に付与し、付与された最大の番号を連結成分の個数とすることにより実行される。カウント部２３６ｃは、この発明の「カウント手段」の一例である。なお、分布画像に対してラベリング処理を施すなどして、連結成分の特定処理とカウント処理とを同時に行うようにしてもよい。
【０１１７】
（サイズ算出部）
サイズ算出部２３６ｄは、連結成分特定部２３６ｂにより特定された各連結成分のサイズを算出する。連結成分のサイズを表す指標としては、面積、径（直径、半径等）、体積などがある。なお、体積は、当該連結成分に含まれる突出領域の体積を意味する。サイズ算出部２３６ｄは、この発明の「サイズ算出手段」の一例である。以下、連結成分のサイズ算出処理の例を説明する。
【０１１８】
まず、連結成分の面積を求める処理例を説明する。各連結成分は、連結と判断された複数の画素の集合である。各画素には予め面積（単位面積）が設定されている。この単位面積は、分布画像やＯＣＴ画像に対して任意に設定されたものであってよい。たとえば、計測倍率を考慮して、一つの画素に相当する実空間における面積を単位面積として設定することができる。サイズ算出部２３６ｄは、各連結成分に含まれる画素数と単位面積との積を演算して当該連結成分の面積とする。
【０１１９】
次に、連結成分の径を求める処理例を説明する。サイズ算出部２３６ｄは、まず上記のように面積を算出する。そして、サイズ算出部２３６ｄは、この面積を有する円の直径（又は半径）を当該連結成分の径とする。なお、連結成分に含まれる最長の線分を探索し、この線分の長さを当該連結成分の径として採用することも可能である。これら以外にも、連結成分を特徴づけることが可能な任意の距離を径として採用することができる。
【０１２０】
次に、連結成分の体積を求める処理例を説明する。前述のように、基準曲線上の各点と層領域との間の奥行方向における距離は、距離算出部２３５ａにより既に算出されている。サイズ算出部２３６ｄは、各連結成分にわたって当該距離を積分することにより、この連結成分（に含まれる突出領域）の体積を算出する。
【０１２１】
なお、サイズの算出方法は上記のものに限定されるものではない。また、サイズを表す指標（次元）も上記のものに限定されるものではない。
【０１２２】
（表示部、操作部）
表示部２４０は、前述した演算制御ユニット２００の表示デバイスを含んで構成される。操作部２５０は、前述した演算制御ユニット２００の操作デバイスを含んで構成される。また、操作部２５０には、眼底解析装置１の筐体や外部に設けられた各種のボタンやキーが含まれていてもよい。たとえば眼底カメラユニット２が従来の眼底カメラと同様の筺体を有する場合、操作部２５０は、この筺体に設けられたジョイスティックや操作パネル等を含んでいてもよい。また、表示部２４０は、眼底カメラユニット２の筺体に設けられたタッチパネルモニタなどの各種表示デバイスを含んでいてもよい。
【０１２３】
なお、表示部２４０と操作部２５０は、それぞれ個別のデバイスとして構成される必要はない。たとえばタッチパネルモニタのように、表示機能と操作機能とが一体化されたデバイスを用いることも可能である。
【０１２４】
〔信号光の走査及びＯＣＴ画像〕
ここで、信号光ＬＳの走査及びＯＣＴ画像について説明しておく。
【０１２５】
眼底解析装置１による信号光ＬＳの走査態様としては、たとえば、水平スキャン、垂直スキャン、十字スキャン、放射スキャン、円スキャン、同心円スキャン、螺旋（渦巻）スキャンなどがある。これらの走査態様は、眼底の観察部位、解析対象（網膜厚など）、走査に要する時間、走査の精密さなどを考慮して適宜に選択的に使用される。
【０１２６】
水平スキャンは、信号光ＬＳを水平方向（ｘ方向）に走査させるものである。水平スキャンには、垂直方向（ｙ方向）に配列された複数の水平方向に延びる走査線に沿って信号光ＬＳを走査させる態様も含まれる。この態様においては、走査線の間隔を任意に設定することが可能である。また、隣接する走査線の間隔を十分に狭くすることにより、前述の３次元画像を形成することができる（３次元スキャン）。垂直スキャンについても同様である。
【０１２７】
十字スキャンは、互いに直交する２本の直線状の軌跡（直線軌跡）からなる十字型の軌跡に沿って信号光ＬＳを走査するものである。放射スキャンは、所定の角度を介して配列された複数の直線軌跡からなる放射状の軌跡に沿って信号光ＬＳを走査するものである。なお、十字スキャンは放射スキャンの一例である。
【０１２８】
円スキャンは、円形状の軌跡に沿って信号光ＬＳを走査させるものである。同心円スキャンは、所定の中心位置の周りに同心円状に配列された複数の円形状の軌跡に沿って信号光ＬＳを走査させるものである。円スキャンは同心円スキャンの一例である。螺旋スキャンは、回転半径を次第に小さく（又は大きく）させながら螺旋状（渦巻状）の軌跡に沿って信号光ＬＳを走査するものである。
【０１２９】
ガルバノミラー４３、４４は互いに直交する方向に信号光ＬＳを走査するように構成されているので、信号光ＬＳをｘ方向及びｙ方向にそれぞれ独立に走査できる。更に、ガルバノミラー４３、４４の向きを同時に制御することにより、ｘｙ面上の任意の軌跡に沿って信号光ＬＳを走査することが可能である。それにより、上記のような各種の走査態様を実現できる。
【０１３０】
上記のような態様で信号光ＬＳを走査することにより、走査線（走査軌跡）に沿った深度方向の断面（ｘｚ断面）の断層像を形成することができる。また、特に走査線の間隔が狭く設定される場合などには、前述の３次元画像を形成することができる。
【０１３１】
上記のような信号光ＬＳの走査対象となる眼底Ｅｆ上の領域、つまりＯＣＴ計測の対象となる眼底Ｅｆ上の領域を走査領域と呼ぶ。３次元スキャンにおける走査領域は、複数の水平スキャンが配列された矩形の領域である。この走査領域は、図８に示す分布画像４００の画像領域に相当する。また、同心円スキャンにおける走査領域は、最大径の円スキャンの軌跡により囲まれる円盤状の領域である。また、放射スキャンにおける走査領域は、各スキャンラインの両端位置を結んだ円盤状（或いは多角形状）の領域である。
【０１３２】
［動作］
眼底解析装置１の動作について説明する。眼底解析装置１の動作例を図１０に示す。
【０１３３】
まず、従来と同様に被検眼Ｅに対するアライメント及びピント合わせを行う（Ｓ１）。次に、主制御部２１１は、ＬＣＤ３９を制御して黄斑計測用の固視標を被検眼Ｅに呈示させる（Ｓ２）。この固視状態で、たとえば３次元スキャンを用いたＯＣＴ計測を実行し、眼底Ｅｆの黄斑及びその周辺の複数の断層像を形成する（Ｓ３）。主制御部２１１は、これら断層像を記憶部２１２に記憶させる。
【０１３４】
たとえば解析処理の開始要求がなされたことに対応し、主制御部２１１は、記憶部２１２から複数の断層像を読み出して解析処理部２３２に送る。
【０１３５】
層領域特定部２３３は、各断層像について、その画素の画素値に基づいて眼底Ｅｆの網膜色素上皮層に相当する層領域を特定する（Ｓ４）。この特定結果は曲線演算部２３４に送られる。
【０１３６】
曲線演算部２３４の特徴部位特定部２３４ａは、各断層像について、ステップ４で特定された層領域中の画素の画素値に基づいて、この層領域の形状に基づく複数の特徴部位を特定する（Ｓ５）。この特定結果は曲線決定部２３４ｂに送られる。
【０１３７】
曲線決定部２３４ｂは、各断層像について、ステップ５で特定された複数の特徴部位に基づいて、その層領域に対する基準曲線を決定する（Ｓ６）。決定された基準曲線の情報は突出領域特定部２３５に送られる。
【０１３８】
突出領域特定部２３５の距離算出部２３５ａは、各断層像について、その基準曲線上の各点とその層領域との間の奥行方向における距離を算出する（Ｓ７）。この算出結果は距離判断部２３５ｂに送られる。
【０１３９】
距離判断部２３５ｂは、各断層像について、ステップ７で算出された当該距離の各演算結果が所定閾値以上であるか判断する（Ｓ８）。この判断結果は画像領域特定部２３５ｃに送られる。
【０１４０】
画像領域特定部２３５ｃは、各断層像について、ステップ８で距離が所定閾値以上と判断された基準曲線上の点の集合と層領域とに挟まれた画像領域を特定して突出領域を求める（Ｓ９）。求められた突出領域の情報は形態情報生成部２３６に送られる。
【０１４１】
形態情報生成部２３６の分布画像形成部２３６ａは、ステップ９で特定された複数の断層像の突出領域に基づいて、眼底Ｅｆの奥行方向に直交するｘｙ平面における突出領域の分布状態を表す分布画像を形成する（Ｓ１０）。形成された分布画像は連結成分特定部２３６ｂに送られる。
【０１４２】
連結成分特定部２３６ｂは、ステップ１０で形成された分布画像の画素の画素値に基づいて、この分布画像における突出領域の連結成分を特定する（Ｓ１１）。この特定結果はカウント部２３６ｃとサイズ算出部２３６ｄにそれぞれ送られる。
【０１４３】
カウント部２３６ｃは、ステップ１１で特定された連結成分の個数をカウントする（Ｓ１２）。
【０１４４】
一方、サイズ算出部２３６ｄは、ステップ１１で特定された各連結成分のサイズを算出する（Ｓ１３）。この動作例では、連結成分のサイズとして、径、面積及び体積を算出するものとする。
【０１４５】
解析処理部２３２は、ステップ１２、１３で得られた形態情報を含む、上記動作例で得られた各種情報に基づいて、解析結果を呈示するための情報（解析結果呈示情報）を生成する（Ｓ１４）。
【０１４６】
主制御部２１１は、この解析結果呈示情報に基づいて表示部２４０に解析結果を表示させる（Ｓ１５）。この解析結果は、眼底Ｅｆにおけるドルーゼンの存在の有無、更には眼底Ｅｆに存在するドルーゼンのサイズや分布等を表すものである。なお、解析結果呈示情報に基づいて解析レポートを印刷出力することも可能である。また、解析結果呈示情報や上記処理により得られた情報、更には患者や被検眼に関する情報などを外部装置に送信したり、記録媒体に記録させたりすることもできる。以上で、この動作例に係る処理は終了となる。
【０１４７】
ステップ１５で表示される画面の例を図１１に示す。解析結果表示画面５００の断層像表示部５１０には、ステップ３で形成された複数の断層像のうちの一枚（代表断層像）Ｇが表示される。代表断層像Ｇは、たとえば３次元スキャンの走査領域の中央位置を通過する走査線に基づく断層像である。なお、断層像表示部５１０に表示される断層像はこれに限定されるものではない。
【０１４８】
代表断層像Ｇ上に基準曲線（図１１では破線で示してある）を表示することもできる。また、眼底表面に相当する画像領域や、網膜色素上皮層に相当する層領域など、代表断層像Ｇ中の所定の画像領域を明示するようにしてもよい。また、代表断層像Ｇ中の突出領域を明示することも可能である。このような明示化は、たとえば、輝度を変更したり、所定の表示色で表示させたり、所定の塗りつぶしパターンで表示させたりするものである。
【０１４９】
分布画像表示部５２０には、ステップ１０で形成された分布画像４００が表示される。分布画像４００は、眼底Ｅｆの黄斑近傍におけるドルーゼンの分布状態を視覚的に描写している。
【０１５０】
なお、眼底Ｅｆの眼底像（撮影画像Ｈ）を取得した場合、分布画像４００と撮影画像Ｈとの合成画像を表示させることができる。この合成画像は、たとえば、撮影画像Ｈに分布画像４００を重ねて表示させたものである。また、合成画像は、撮影画像Ｈと分布画像４００を合成して得られる単一の画像でもよい。このときの画像位置合わせ処理は、相互の画像中の特徴部位（黄斑、血管の分岐部、病変部など）を一致させるようにして行うことができる。また、各断層像の画素の画素値をｚ方向に積算して２次元画像を形成し、この２次元画像と撮影画像Ｈとの位置合わせを行い、その位置合わせ結果に基づいて分布画像と撮影画像Ｈとの位置合わせを行うことも可能である。
【０１５１】
解析結果表示部５３０には、ドルーゼンの可能性がある突出領域の解析結果が表示される。この表示例では、個数分布グラフ５３１、面積分布グラフ５３２及び体積分布グラフ５３３が表示される。これら各分布は、特定された連結成分をサイズ（直径が６４μｍ以下、６４μｍ超かつ１２５μｍ以下、１２５μｍ超かつ２５０μｍ以下、２５０μｍ超の４つの範囲）で分類して得られたものである。この分類処理は解析処理部２３２が実行する。
【０１５２】
個数分布グラフ５３１は、ステップ１２でカウントされた突出領域の連結成分の個数の分布を表す。個数分布グラフ５３１の下方には、連結成分の個数（カウント数）が表示される。連結成分の個数や個数分布は、眼底Ｅｆの黄斑近傍に存在するドルーゼンの個数や個数分布に相当する。
【０１５３】
面積分布グラフ５３２は、ステップ１３で算出された連結成分の面積の分布を表す。面積分布グラフ５３２の下方には、連結成分の面積の平均値が表示される。連結成分の面積や面積分布は、眼底Ｅｆの黄斑近傍に存在するドルーゼンの面積や面積分布に相当する。なお、この面積は、−ｚ方向から＋ｚ方向を見たときのｘｙ平面上における面積を表している。
【０１５４】
体積分布グラフ５３３は、ステップ１３で算出された連結成分の体積の分布を表す。体積分布グラフ５３３の下方には、連結成分の体積の平均値が表示される。連結成分の体積や体積分布は、眼底Ｅｆの黄斑近傍に存在するドルーゼンの体積や体積分布に相当する。
【０１５５】
［作用・効果］
以上のような眼底解析装置１の作用及び効果について説明する。
【０１５６】
眼底解析装置１は、眼底Ｅｆの層構造をそれぞれ描写する複数の断層像のそれぞれについて、その断層像の画素の画素値に基づいて網膜色素上皮層に相当する層領域を特定し、この層領域の形状に基づいて基準曲線を求め、この層領域と基準曲線とに基づいて層領域の突出領域を特定し、この突出領域の形態を表す形態情報を生成するように作用する。
【０１５７】
形態情報は、個々のドルーゼンに対応する突出領域の連結成分を特定してその個数やサイズを求めることにより形成される。また、形態情報は、連結成分（ドルーゼン）の個数やサイズの分布を表す情報（たとえば図１１の円グラフ５３１〜５３３）を含んでいてもよい。更に、形態情報は、層領域の突出形態や突出領域を描写した断層像や、突出領域の分布を表す分布画像や、眼底の撮影画像と分布画像との合成画像などの画像情報を含んでいてもよい。
【０１５８】
このような解析処理を断層像に施すことで、ドルーゼンの可能性がある網膜色素上皮層の突出部分をそれぞれ個別に把握することができる。それにより、小さなドルーゼンであっても検出することが可能となる。また、断層像に基づいてドルーゼンに特徴的な突出領域を特定するように構成されているので、ドルーゼンを効果的に検出することが可能である。
【０１５９】
なお、この発明に係る眼底解析装置は、ＯＣＴ計測機能を有する眼底解析装置１全体であってもよいし、上記解析処理を実行する演算制御ユニット２００のみであってもよい。すなわち、この発明に係る眼底解析装置は、上記解析処理を実行可能なＯＣＴ装置であってもよいし、上記解析処理を実行可能なコンピュータであってもよい。このＯＣＴ装置は、ＯＣＴ計測を行って眼底の断層像を形成し、この断層像に対して上記解析処理を施す。なお、外部装置から取得した断層像に対して上記解析処理を施すことも可能である。一方、このコンピュータは、ＯＣＴ装置、データベース等の外部装置や記憶媒体から断層像を取得して記憶し、更に、この断層像に対して上記解析処理を施す。
【０１６０】
［変形例］
以上に説明した構成は、この発明を好適に実施するための一例に過ぎない。よって、この発明の要旨の範囲内における任意の変形を適宜に施すことが可能である。
【０１６１】
上記実施形態では、眼底の断層像（２次元画像）に対して解析処理を施す場合について説明したが、３次元画像に対して同様の解析処理を施して同様の形態情報を生成することも可能である。このような眼底解析装置は、上記実施形態と同様の構成を有する。以下、この変形例に係る眼底解析装置について上記実施形態（特に図３、図４）を参照しつつ説明する。
【０１６２】
この変形例に係る眼底解析装置の記憶部２１２（記憶手段）は、眼底Ｅｆの層構造を描写する３次元画像を記憶する。この３次元画像は、３次元画像形成部２３１（３次元画像形成手段）により形成されたものでもよいし、外部装置から取得されたものでもよい。
【０１６３】
解析処理部２３２の層領域特定部２３３（層領域特定手段）は、３次元画像の画素（ボクセル等）の画素値に基づいて、網膜色素上皮層に相当する３次元画像中の層領域を特定する。特定される層領域は曲面（いくつか特異点を含む可能性はある）である。
【０１６４】
この変形例では、上記実施形態の曲線演算部２３４の代わりに曲面演算部（曲面演算手段）が設けられる。曲面演算部は、層領域特定部２３３により特定された層領域の形状に基づいて、眼底Ｅｆの奥行方向（＋ｚ方向）に凸な曲面（基準曲面）を求める。この処理は、たとえば上記実施形態と同様に最小二乗法を用いたり、２次曲面の方程式に座標値を代入して係数を求めたり、２次曲面の係数パラメータを変更して適当な２次曲面を求めたりすることによって実行できる。
【０１６５】
次に、突出領域特定部２３５（突出領域特定手段）は、層領域特定部２３３が特定した層領域と、曲面演算部が求めた曲面とに基づいて、眼底Ｅｆの奥行方向の逆方向（−ｚ方向）に向かって層領域が突出する突出領域を特定する。
【０１６６】
続いて、形態情報生成部２３６（形態情報生成手段）は、突出領域特定部２３５により特定された突出領域の形態を表す形態情報を生成する。形態情報としては、上記実施形態と同様にドルーゼン（突出領域の連結成分）の個数、サイズ、分布、代表断層像、合成画像等を呈示できるとともに、３次元画像をレンダリングして得られる擬似的な３次元画像を呈示することも可能である。
【０１６７】
このような変形例によれば、上記実施形態と同様に、ドルーゼンの可能性がある網膜色素上皮層の突出部分をそれぞれ個別に把握することができるので、小さなドルーゼンであっても検出することが可能となる。また、３次元画像に基づいてドルーゼンに特徴的な突出領域を特定するように構成されているので、ドルーゼンを効果的に検出することが可能である。
【０１６８】
なお、この変形例では、３次元画像処理を実行するので、２次元画像処理を実行する上記実施形態と比較して処理に使用されるリソースが増大する。一方、この変形例ではドルーゼンの３次元的な形態をより高精度、高確度に把握できるというメリットがある。ただし、上記実施形態において走査線の間隔を小さくすることにより、ドルーゼン検出の精度や確度の向上を図ったり、ドルーゼンの形態を把握する精度や確度の向上を図ったりすることができる。
【０１６９】
他の変形例を説明する。上記の実施形態や変形例では、網膜色素上皮層に相当する層領域と基準曲線や基準曲線との間の距離に基づいてドルーゼン（と考えられる突出領域）を検出しているが、ＯＣＴ計測の感度を向上させるなどしてブルッフ膜に相当する画像領域（膜領域と呼ぶ）を検出し、層領域と膜領域との間の距離に基づいてドルーゼンを検出するように構成してもよい。ドルーゼンはブルッフ膜と網膜色素上皮層との間に発生するものであるから、この変形例に係る処理を実行することで、ドルーゼンの形態をより高精度、高確度で把握することが可能となる。
【０１７０】
なお、上記実施形態は、ＯＣＴ画像からの検出が困難なブルッフ膜の代わりに、ドルーゼン（突出）が無いと仮定した状態の網膜色素上皮層の形態を近似する基準曲線や基準曲面を用いるものである。
【０１７１】
更なる変形例として、ＯＣＴ画像から網膜色素上皮層に相当する層領域を特定するとともに、ブルッフ膜に相当する膜領域の特定を試みることができる。膜領域が特定された場合には、膜領域と層領域との間の距離に基づいて突出領域を特定して形態情報を生成する。一方、膜領域の特定に失敗した場合には、上記実施形態のように基準曲線（又は基準曲面）を求め、層領域と基準曲線（又は基準曲面）とに基づいて突出領域を特定して形態情報を生成する。
【０１７２】
この変形例によれば、膜領域が特定された場合には高精度、高確度でドルーゼンの形態を把握できるとともに、膜領域が特定されなかった場合には基準曲線（又は基準曲面）を利用してドルーゼンの形態を把握できる。
【０１７３】
上記の実施形態においては、参照ミラー１１４の位置を変更して信号光ＬＳの光路と参照光ＬＲの光路との光路長差を変更しているが、光路長差を変更する手法はこれに限定されるものではない。たとえば、被検眼Ｅに対して眼底カメラユニット２やＯＣＴユニット１００を移動させて信号光ＬＳの光路長を変更することにより光路長差を変更することができる。また、特に被測定物体が生体部位でない場合などには、被測定物体を深度方向（ｚ軸方向）に移動させることにより光路長差を変更することも有効である。
【０１７４】
［眼底解析プログラム］
この発明に係る眼底解析プログラムについて説明する。この眼底解析プログラムは、眼底の層構造をそれぞれ描写する複数の断層像を記憶する記憶手段を有するコンピュータに後述の動作を実行させる。このコンピュータの例として、上記実施形態の演算制御ユニットがある。この眼底解析プログラムは、このコンピュータ自体に記憶されていてもよいし、このコンピュータと通信可能に接続されたサーバ等に記憶されていてもよい。
【０１７５】
この眼底解析プログラムに基づいてコンピュータが実行する動作について説明する。まず、コンピュータは、記憶手段に記憶された各断層像を読み出し、上記実施形態の層領域特定部２３３と同様に、その断層像の画素の画素値に基づいて、網膜色素上皮層に相当する当該断層像中の層領域を特定する。次に、コンピュータは、上記実施形態の曲線演算部２３４と同様に、特定された層領域の形状に基づいて、眼底の奥行方向に凸な曲線（基準曲線）を求める。続いて、コンピュータは、上記実施形態の突出領域特定部２３５と同様に、層領域と基準曲線とに基づいて、眼底の奥行方向の逆方向に向かって層領域が突出する突出領域を特定する。そして、コンピュータは、上記実施形態の形態情報生成部２３６と同様に、特定された突出領域の形態を表す形態情報を生成する。コンピュータは、生成された形態情報を表示出力或いは印刷出力する。
【０１７６】
このような眼底解析プログラムによれば、上記実施形態と同様に、ドルーゼンの可能性がある網膜色素上皮層の突出部分をそれぞれ個別に把握することができるので、小さなドルーゼンであっても検出することが可能となる。また、断層像に基づいてドルーゼンに特徴的な突出領域を特定するように構成されているので、ドルーゼンを効果的に検出することが可能である。
【０１７７】
この発明に係る他の眼底解析プログラムについて説明する。この眼底解析プログラムは、眼底の層構造を描写する３次元画像を記憶する記憶手段を有するコンピュータに後述の動作を実行させる。このコンピュータの例として、上記変形例の演算制御ユニットがある。
【０１７８】
この眼底解析プログラムに基づいてコンピュータが実行する動作について説明する。まず、コンピュータは、記憶手段に記憶された３次元画像を読み出し、上記変形例の層領域特定部２３３と同様に、この３次元画像の画素の画素値に基づいて、網膜色素上皮層に相当する３次元画像中の層領域を特定する。次に、コンピュータは、上記変形例の曲面演算部と同様に、特定された層領域の形状に基づいて、眼底の奥行方向に凸な曲面（基準曲面）を求める。続いて、コンピュータは、上記変形例の突出領域特定部２３５と同様に、層領域と基準曲面とに基づいて、眼底の奥行方向の逆方向に向かって層領域が突出する突出領域を特定する。そして、コンピュータは、上記変形例の形態情報生成部２３６と同様に、特定された突出領域の形態を表す形態情報を生成する。コンピュータは、生成された形態情報を表示出力或いは印刷出力する。
【０１７９】
このような眼底解析プログラムによれば、上記変形例と同様に、ドルーゼンの可能性がある網膜色素上皮層の突出部分をそれぞれ個別に把握することができるので、小さなドルーゼンであっても検出することが可能となる。また、３次元画像に基づいてドルーゼンに特徴的な突出領域を特定するように構成されているので、ドルーゼンを効果的に検出することが可能である。
【０１８０】
この発明に係る眼底解析プログラムは、上記実施形態や変形例で説明した各種の処理を更にコンピュータに実行させるように構成されていてもよい。
【０１８１】
この発明に係る眼底解析プログラムを、コンピュータによって読み取り可能な任意の記録媒体に記憶させることができる。この記録媒体としては、たとえば、光ディスク、光磁気ディスク（ＣＤ−ＲＯＭ／ＤＶＤ−ＲＡＭ／ＤＶＤ−ＲＯＭ／ＭＯ等）、磁気記憶媒体（ハードディスク／フロッピー（登録商標）ディスク／ＺＩＰ等）などを用いることが可能である。また、ハードディスクドライブやメモリ等の記憶装置に記憶させることも可能である。
【０１８２】
また、インターネットやＬＡＮ等のネットワークを通じてこのプログラムを送受信することも可能である。
【０１８３】
［眼底解析方法］
この発明に係る眼底解析方法について説明する。この眼底解析方法は、眼底の層構造をそれぞれ描写する複数の断層像を解析するものであり、次のようなステップを含んで構成される。
【０１８４】
第１のステップとして、複数の断層像のそれぞれについて、その断層像の画素の画素値に基づいて、網膜色素上皮層に相当する当該断層像中の層領域を特定する。次に、第２のステップとして、特定された層領域の形状に基づいて、眼底の奥行方向に凸な曲線（基準曲線）を求める。続いて、第３のステップとして、層領域と基準曲線とに基づいて、眼底の奥行方向の逆方向に向かって層領域が突出する突出領域を特定する。そして、第４のステップとして、特定された突出領域の形態を表す形態情報を生成する。
【０１８５】
このような眼底解析方法によれば、上記実施形態と同様に、ドルーゼンの可能性がある網膜色素上皮層の突出部分をそれぞれ個別に把握することができるので、小さなドルーゼンであっても検出することが可能となる。また、断層像に基づいてドルーゼンに特徴的な突出領域を特定するように構成されているので、ドルーゼンを効果的に検出することが可能である。
【０１８６】
この発明に係る他の眼底解析方法について説明する。この眼底解析方法は、眼底の層構造を描写する３次元画像を解析するものであり、次のようなステップを含んで構成される。
【０１８７】
第１のステップとして、３次元画像の画素の画素値に基づいて、網膜色素上皮層に相当する３次元画像中の層領域を特定する。次に、第２のステップとして、特定された層領域の形状に基づいて、眼底の奥行方向に凸な曲面（基準曲面）を求める。続いて、第３のステップとして、層領域と基準曲面とに基づいて、眼底の奥行方向の逆方向に向かって層領域が突出する突出領域を特定する。そして、第４のステップとして、特定された突出領域の形態を表す形態情報を生成する。
【０１８８】
このような眼底解析方法によれば、上記変形例と同様に、ドルーゼンの可能性がある網膜色素上皮層の突出部分をそれぞれ個別に把握することができるので、小さなドルーゼンであっても検出することが可能となる。また、３次元画像に基づいてドルーゼンに特徴的な突出領域を特定するように構成されているので、ドルーゼンを効果的に検出することが可能である。
【０１８９】
この発明に係る眼底解析方法は、上記実施形態や変形例で説明した各種の処理を実行するステップを更に含んでいてもよい。
【符号の説明】
【０１９０】
１  眼底解析装置
２  眼底カメラユニット
１０  照明光学系
３０  撮影光学系
４３、４４  ガルバノミラー
７０  走査駆動部
１００  ＯＣＴユニット
１０１  光源ユニット
１１４  参照ミラー
１１８  回折格子
１２０  ＣＣＤイメージセンサ
１３０  参照駆動部
２００  演算制御ユニット
２１０  制御部
２１１  主制御部
２１２  記憶部
２２０  画像形成部
２３０  画像処理部
２３１  ３次元画像形成部
２３２  解析処理部
２３３  層領域特定部
２３４  曲線演算部
２３５  突出領域特定部
２３６  形態情報生成部
２４０  表示部
２５０  操作部
３００  層領域
３２０  基準曲線（二次曲線）
３３０、３４０、Ｔ  突出領域
４００  分布画像
５００  解析結果表示画面
Ｅ  被検眼
Ｅｆ  眼底
Ｇ  断層像、代表断層像
Ｔｋ  連結成分
【図３】
0005474435.tif 000002
【図４】
0005474435.tif 000003
【図５】
0005474435.tif 000004
【図６】
0005474435.tif 000005
【図８】
0005474435.tif 000006
【図９】
0005474435.tif 000007
【図１０】
0005474435.tif 000008
【図１１】
0005474435.tif 000009
【図１】
0005474435.tif 000010
【図２】
0005474435.tif 000011
【図７】
0005474435.tif 000012