JP2018020192A - 画像処理システム、処理方法及びプログラム - Google Patents

Abstract

【課題】被検眼が疾病であるか否かの判定基準となる指標とともに、診断対象となる被検眼の形状特徴量を表示できるようにした画像処理システム、画像処理方法及び画像処理プログラムを提供する。
【解決手段】画像処理システム１０は、眼底の断層画像を得るための断層画像撮影装置２０により取得された前記断層画像を処理して、前記眼底の網膜層の形状情報を取得する解析手段と、眼軸長で分類された複数の眼に関する網膜層の形状情報が記憶されたデータベースを参照して、前記網膜層の形状情報に関する解析情報を取得する取得手段とを具備する。
【選択図】図１

Description

本発明は、画像処理システム、処理方法及びプログラムに関する。

低コヒーレント光による干渉を利用した光干渉断層法（ＯＣＴ：Optical Coherence Tomography）を用いる断層画像撮影装置が知られている。このような断層画像撮影装置により眼底を撮影することによって、網膜層内部の状態を３次元的に観察することができる。

断層画像撮影装置による撮影は、疾病の診断をより的確に行なうのに有用な技術であるため、近年注目を集めている。このようなＯＣＴの形態としては、例えば、広帯域な光源とマイケルソン干渉計とを組み合わせたＴＤ−ＯＣＴ（Time domain OCT）が知られている。ＴＤ−ＯＣＴでは、参照アームの遅延を走査することで信号アームの後方散乱光との干渉光を計測し、深さ分解の情報を得ることができる。

しかし、ＴＤ−ＯＣＴでは、画像を高速に取得することが難しい。そのため、画像をより高速に取得する方法として、広帯域光源を用いて、分光器でインターフェログラムを取得するＳＤ−ＯＣＴ（Spectral domain OCT）が知られている。また、光源として、高速波長掃引光源を用いることで、単一チャネル光検出器でスペクトル干渉を計測するＳＳ−ＯＣＴ（Swept Source OCT）も知られている。

ここで、これらのＯＣＴで撮影された断層画像において、網膜の形態変化を計測できれば、緑内障などの疾病の進行度や治療後の回復具合を定量的に診断することができる。このような技術に関連して、特許文献１では、網膜の形態変化を定量的に計測するために、コンピュータを用いて断層画像から網膜の各層の境界を検出し、当該検出結果に基づいて層の厚みを計測する技術が開示されている。

特開２００８−０７３０９９号公報

上述した特許文献１の技術では、２次元画像上で指定された範囲の断層画像を取得し、層境界を検出して層厚を求めている。しかし、網膜層の３次元的な形状解析を行なっておらず、また、形状解析の結果を効果的に表示するものではない。

本発明は、上記課題に鑑みてなされたものであり、被検眼が疾病であるか否かの判定基準となる指標とともに、診断対象となる被検眼の形状特徴量を表示できるようにした技術を提供することを目的とする。

上記課題を解決するため、本発明の一態様による画像処理システムは、眼底の断層画像を得るための断層画像撮影装置により取得された前記断層画像を処理して、前記眼底の網膜層の形状情報を取得する画像処理システムであって、
眼軸長で分類された複数の眼に関する網膜層の形状情報が記憶されたデータベースを参照して、前記網膜層の形状情報に関する解析情報を取得する解析処理手段を備えることを特徴する。

本発明によれば、被検眼が疾病であるか否かの判定基準となる指標とともに、診断対象となる被検眼の形状特徴量を表示できる。

本発明の一実施の形態に係わる画像処理システム１０の構成の一例を示す図。 断層画像撮影画面６０の一例を示す図。 図１に示す画像処理装置３０の処理の流れの一例を示すフローチャート。 ３次元の形状解析処理の概要を説明するための図。 ３次元の形状解析処理の概要を説明するための図。 断層画像観察画面８０の一例を示す図。 断層画像観察画面８０の各部の一例を示す図。 断層画像観察画面８０の各部の一例を示す図。 断層画像観察画面８０の各部の一例を示す図。 断層画像観察画面８０の一例を示す図。 断層画像観察画面８０の各部の一例を示す図。 断層画像観察画面８０の一例を示す図。 断層画像観察画面４００の一例を示す図。 断層画像観察画面４００の各部の一例を示す図。 断層画像観察画面４００の一例を示す図。 断層画像観察画面４００の各部の一例を示す図。 実施形態３に係わる画像処理システム１０の構成の一例を示す図。 実施形態３に係わる画像処理装置３０の処理の流れの一例を示すフローチャート。 網膜層の形状特徴を示す特徴量の取得方法の概要を示す図。 網膜層の形状特徴を示す特徴量の取得方法の概要を示す図。 断層画像観察画面８０の一例を示す図。 実施形態４に係わる画像処理システム１０の構成の一例を示す図。 断層画像観察画面９００の一例を示す図。

以下、本発明に係わる実施の形態について添付図面を参照して詳細に説明する。

（実施形態１）
図１は、本発明の一実施の形態に係わる画像処理システム１０の構成の一例を示す図である。

画像処理システム１０には、画像処理装置３０と、断層画像撮影装置２０と、眼底画像撮影装置５１と、外部記憶装置５２と、表示装置５３と、入力装置５４とが具備される。

断層画像撮影装置２０は、例えば、ＳＤ−ＯＣＴやＳＳ−ＯＣＴ等で実現され、低コヒーレント光による干渉を利用した光干渉断層法（ＯＣＴ）を用いて、眼底の３次元形状を示す断層画像を撮影する。断層画像撮影装置２０には、ガルバノミラー２１と、駆動制御部２２と、パラメータ設定部２３と、内部固視灯２４と、コヒーレンスゲートステージ２５とが具備される。

ガルバノミラー２１は、測定光（照射光）を眼底に２次元的に走査させる機能を有し、断層画像撮影装置２０による眼底の撮影範囲を規定する。ガルバノミラー２１は、例えば、Ｘスキャン用のミラーとＹスキャン用のミラーとの２枚のミラーで構成され、被検眼の眼底に対して測定光を光軸に対して直交する面上に走査させる。

駆動制御部２２は、ガルバノミラー２１の駆動（走査）範囲及び速度を制御する。これにより、眼底における平面方向（測定光の光軸方向に直交する方向）の撮影範囲及び走査線数（平面方向の走査速度）が規定される。

パラメータ設定部２３は、駆動制御部２２によるガルバノミラー２１の駆動制御に用いられる各種パラメータを設定する。このパラメータによって、断層画像撮影装置２０による断層画像の撮影条件が決定される。例えば、走査線の走査位置、走査線数、撮影枚数などが決定される。その他、内部固視灯の位置、走査範囲（スキャン範囲）、走査パターン（スキャンパターン）、コヒーレンスゲート位置の設定等も行なわれる。なお、パラメータの設定は、画像処理装置３０からの指示に基づいて行なわれる。

内部固視灯２４は、断層画像の撮影中、眼球運動が起きないように輝点を視野内に置き、視点が動くことを抑制する。内部固視灯２４は、表示部２４ａとレンズ２４ｂとを具備して構成される。表示部２４ａは、例えば、複数の発光ダイオード（ＬＤ）がマトリックス状に配置されて実現される。発光ダイオードの点灯位置は、駆動制御部２２の制御により、撮影対象となる部位に合わせて変更される。表示部２４ａからの光は、レンズ２４ｂを介して被検眼に導かれる。表示部２４ａから出射される光は、例えば、５２０ｎｍで、駆動制御部２２の制御により、所望のパターンが表示（点灯）される。

コヒーレンスゲートステージ２５は、被検眼の眼軸長の相違等に対応するために設けられる。より具体的には、（測定光と干渉させるための）参照光の光路長を制御し、眼底の深さ方向（光軸方向）に沿った撮影位置を調整する。これにより、眼軸長が異なる被検眼に対しても、参照光と測定光との光路長を一致させることができる。なお、コヒーレンスゲートステージ２５は、駆動制御部２２により制御される。

ここで、コヒーレンスゲートとは、断層画像撮影装置２０における測定光と参照光との光学距離が等しい位置を示す。コヒーレンスゲートの位置の制御により、網膜層側の撮影か、或いは網膜層よりも深部側とするＥＤＩ（Enhanced Depth Imaging）法による撮影かが切り替えられる。ＥＤＩ法で撮影を行なった場合、コヒーレンスゲートの位置が網膜層の深部側に設定されるため、断層画像の上部側を越えて網膜層が撮影された場合に、網膜層が断層画像の中に折り返して写ることがなくなることになる。

眼底画像撮影装置５１は、例えば、眼底カメラやＳＬＯ（Scanning Laser Ophthalmoscope）等で実現され、眼底の（２次元）眼底画像を撮影する。

外部記憶装置５２は、例えば、ＨＤＤ（Hard Disk Drive）等で実現され、各種データを記憶する。外部記憶装置５２は、例えば、被検眼に関する情報（患者の氏名、年齢、性別など）に対して、撮影した画像データ、撮影パラメータ、画像解析パラメータ、操作者によって設定されたパラメータをそれぞれ関連付けて保持する。

入力装置５４は、例えば、マウス、キーボード、タッチ操作画面等で実現され、操作者からの各種指示を入力する。例えば、操作者は、入力装置５４を介して、画像処理装置３０、断層画像撮影装置２０及び眼底画像撮影装置５１に対して各種指示や設定等を行なう。表示装置５３は、例えば、液晶ディスプレイ等で実現され、各種情報を操作者に向けて表示（提示）する。

画像処理装置３０は、例えば、パーソナルコンピュータ等で実現され、各種画像を処理する。すなわち、画像処理装置３０には、コンピュータが内蔵されている。コンピュータには、ＣＰＵ（Central Processing Unit）等の主制御手段、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）等の記憶手段などが具備される。

ここで、画像処理装置３０には、その機能的な構成として、画像取得部３１と、記憶部３２と、画像処理部３３と、指示部３４と、表示制御部３５とが具備される。なお、記憶部３２以外の構成は、例えば、ＣＰＵが、ＲＯＭ等に記憶されたプログラムを読み出し実行することにより実現される。

画像処理部３３には、検出部４１と、解析部４２と、判定部４３と、位置合わせ部４７とが具備される。

画像取得部３１は、断層画像撮影装置２０により撮影された断層画像と、眼底画像撮影装置５１により撮影された眼底画像とを取得し、記憶部３２に格納する。なお、記憶部３２は、例えば、ＲＯＭやＲＡＭ等で実現される。

検出部４１は、記憶部３２に記憶されている断層画像から網膜の層（網膜層）を検出する。

解析部４２は、解析対象となる網膜層の解析を行なう。解析部４２には、判定部４３と、網膜層解析部４４と、解析結果作成部４５と、形状データ作成部４６とが設けられる。

判定部４３は、撮影モード（近視眼解析用の撮影モード、非近視眼解析用の撮影モード）に応じて、網膜層の３次元の形状解析処理を行なうか否かを判定する。なお、３次元の形状解析とは、３次元の形状データを作成し、当該形状データを用いて網膜層の形状解析を行なうことを指す。

網膜層解析部４４は、判定部４３の判定結果に基づいて、解析対象となる網膜層に対して解析処理を行なう。なお、本実施形態においては、３次元の形状解析処理として、近視眼の黄斑の解析を行なう場合について説明する。解析結果作成部４５では、解析結果（網膜層の状態を示す情報）を提示するための各種データを作成する。形状データ作成部４６では、撮影により得られた複数の断層画像同士を位置合わせし、３次元の形状データを作成する。すなわち、３次元の形状データは、網膜層の層情報に基づいて作成される。

位置合わせ部４７では、解析結果と眼底画像との位置合わせや、眼底画像同士の位置合わせ等を行なう。指示部３４は、断層画像撮影装置２０で設定されている撮影モードに応じて、撮影パラメータ等の情報の指示を行なう。

以上が、画像処理システム１０の構成の一例についての説明である。なお、上述した各装置に設けられた機能的な構成は、必ずしも図示した通りに実現される必要はなく、システム内におけるいずれかの装置にその全部若しくは一部が実現されていれば良い。例えば、図１においては、外部記憶装置５２、表示装置５３や入力装置５４が、画像処理装置３０の外部に設けられていたが、これら構成を画像処理装置３０の内部に設けるようにしても良い。また、例えば、画像処理装置と断層画像撮影装置２０とが一体的に構成されていても良い。

次に、図２を用いて、図１に示す表示装置５３に表示される断層画像撮影画面６０の一例について説明する。なお、この画面は、断層画像の撮影時に表示される。

断層画像撮影画面６０は、断層画像の表示領域６１と、眼底画像の表示領域６２と、撮影モード設定用のコンボボックス６３と、撮影を指示するキャプチャーボタン６４とを具備して構成される。なお、眼底画像の表示領域６２内の符号６５は、撮影領域を示すマークであり、眼底画像上に重畳表示される。また、Ｍは黄斑部を示し、Ｄは視神経乳頭部を示し、Ｖは血管を示している。

コンボボックス６３では、例えば、（黄斑部の）近視眼解析用の撮影モードか、（黄斑部の）非近視眼解析用の撮影モードかの設定が行なわれる。この場合、近視眼解析用の撮影モードが設定されている。

ここで、断層画像の表示領域６１には、眼底の断層画像が表示されている。Ｌ１は内境界膜（ＩＬＭ）を示し、Ｌ２は神経線維層（ＮＦＬ）と神経節細胞層（ＧＣＬ）との境界を示し、Ｌ３は視細胞内節外節接合部（ＩＳＯＳ）を示している。また、Ｌ４は網膜色素上皮層（ＲＰＥ）を示し、Ｌ５はブルッフ膜（ＢＭ）を示している。上述した検出部４１においては、Ｌ１〜Ｌ５の各境界のいずれかを検出する。

次に、図３（ａ）及び図３（ｂ）を用いて、図１に示す画像処理装置３０の処理の流れの一例について説明する。まず、図３（ａ）を用いて、断層画像の撮影時の全体的な処理の流れについて説明する。

［Ｓ１０１］
画像処理装置３０は、まず、被検眼を同定する情報として外部から被検者識別番号を取得する。そして、当該被検者識別番号に基づいて、外部記憶装置５２が保持している当該被検眼に関する情報を取得して記憶部３２に記憶する。

［Ｓ１０２］
画像処理装置３０は、画像取得部３１において、撮影時の撮影位置確認用のプレスキャン画像として眼底画像撮影装置５１から眼底画像を取得するとともに、断層画像撮影装置２０から断層画像を取得する。

［Ｓ１０３］
画像処理装置３０は、撮影モードの設定を行なう。撮影モードの設定は、図２で説明した撮影モード設定用のコンボボックス６３からの操作者の選択に基づいて行なわれる。ここでは、近視眼解析用の撮影モードで撮影を行なう場合について説明する。

［Ｓ１０４］
画像処理装置３０は、指示部３４において、当該コンボボックス６３から設定された撮影モードに応じた撮影パラメータを断層画像撮影装置２０に指示する。これにより、断層画像撮影装置２０は、パラメータ設定部２３において、当該指示に従って撮影パラメータを設定する。具体的には、内部固視灯の位置、走査範囲（スキャン範囲）、走査パターン（スキャンパターン）、コヒーレンスゲート位置の少なくともいずれかの設定を指示する。

ここで、近視眼解析用の撮影モード時のパラメータ設定について説明する。内部固視灯の位置のパラメータ設定では、黄斑部中心を撮影できるように内部固視灯２４の位置が設定される。画像処理装置３０は、駆動制御部２２において、当該撮影パラメータに従って表示部２４ａの発光ダイオードを制御する。なお、撮影範囲が充分に広く取れるような構成の装置であれば、黄斑部と視神経乳頭部との中心が撮影の中心となるように制御しても良い。このような制御を行なうのは、近視眼における網膜層の形状解析を行なうために、黄斑部が写る領域を撮影するためである。

また、走査範囲のパラメータ設定では、例えば、装置の撮影範囲の限界値として９〜１５ｍｍの範囲が設定される。この数値はあくまで一例であり、装置のスペックに応じて適宜変更すれば良い。なお、撮影範囲は、形状変化が起こっている箇所を漏れなく検出できるようにするため、より広い領域が好ましい。

走査パターンのパラメータ設定では、網膜層の３次元形状を撮影できるように、例えば、ラスタスキャンや放射状スキャンが設定される。

コヒーレンスゲートの位置のパラメータ設定では、ＥＤＩ法で撮影を行なえるようにゲート位置が設定される。高度近視眼においては、網膜層の湾曲の程度がきつくなり、断層画像の上部側を越えて網膜層が撮影されてしまう。その場合、断層画像の上部を越えた網膜層が断層画像内に折り返って写ってしまうので、それを防ぐためこのようなパラメータの設定が必要となる。なお、断層画像撮影装置２０として高深達のＳＳ−ＯＣＴを用いる場合には、網膜層の位置をゲート位置から遠ざけても良好な断層画像が得られるため、必ずしもＥＤＩ法で撮影を行なわなくても良い。

［Ｓ１０５］
画像処理装置３０は、指示部３４において、断層画像撮影装置２０に対して被検眼の撮影を指示する。この指示は、例えば、操作者が、入力装置５４を介して、断層画像撮影画面６０のキャプチャーボタン６４を押下すると行なわれる。この指示が行なわれると、断層画像撮影装置２０は、パラメータ設定部２３により設定された撮影パラメータに基づき、駆動制御部２２を制御させる。これにより、ガルバノミラー２１が動作し、断層画像の撮影が行なわれる。

上述した通り、ガルバノミラー２１は、水平方向用のＸスキャナと垂直方向用のＹスキャナとで構成されている。そのため、これらのスキャナの向きをそれぞれ変更することで、装置座標系における水平方向（Ｘ）及び垂直方向（Ｙ）それぞれの方向に沿って断層画像を撮影することができる。これらのスキャナの向きを同時に変更させることで、水平方向と垂直方向とを合成した方向に走査することができるため、眼底平面上の任意の方向に沿った撮影が可能となる。このとき、画像処理装置３０は、表示制御部３５において、当該撮影した断層画像を表示装置５３上に表示させる。これにより、操作者は、撮影結果を確認することができる。

［Ｓ１０６］
画像処理装置３０は、画像処理部３３において、記憶部３２に記憶された断層画像から網膜層の検出／解析を行なう。すなわち、Ｓ１０５の処理で撮影された断層画像に対して網膜層の検出／解析処理を行なう。

［Ｓ１０７］
画像処理装置３０は、断層画像の撮影を終了するか否かを判定する。この判定は、入力装置５４を介した操作者からの指示に基づいて行なわれる。すなわち、操作者が、終了を指示したか否かに基づいて断層画像の撮影を終了するか否かが判定される。

撮影終了が指示された場合、画像処理装置３０は、この処理を終了する。一方、処理を終了せずに撮影を続ける場合、画像処理装置３０は、再度、Ｓ１０２以降の処理を行なう。

なお、後述するＳ２０６やＳ２０９の処理によって、操作者が、網膜層の検出結果や、眼底画像とマップとの位置をマニュアル修正した場合等には、撮影の終了に際して、画像処理装置３０は、それに伴って変更された撮影パラメータを外部記憶装置５２に保存する。このとき、当該変更されたパラメータを保存するか否かの確認ダイアログを表示して、撮影パラメータの変更を行なうか否かを操作者に問い合わせるようにしても良い。

次に、図３（ｂ）を用いて、図３（ａ）のＳ１０６の示す網膜層の検出／解析処理について説明する。

［Ｓ２０１］
網膜層の検出／解析処理が開始すると、画像処理装置３０は、検出部４１において、断層画像から網膜層を検出する。この処理について、図２に示す断層画像（表示領域６１）を用いて具体的に説明する。黄斑部の場合、検出部４１は、まず、断層画像に対して、メディアンフィルタとＳｏｂｅｌフィルタとをそれぞれ適用して画像を作成する（以下、メディアン画像、Ｓｏｂｅｌ画像とする）。続いて、検出部４１は、当該作成したメディアン画像とＳｏｂｅｌ画像とから、Ａスキャン毎にプロファイルを作成する。メディアン画像では輝度値のプロファイルとなり、Ｓｏｂｅｌ画像では勾配のプロファイルとなる。そして、検出部４１は、Ｓｏｂｅｌ画像から作成したプロファイル内のピークを検出する。最後に、当該検出されたピークの前後やピーク間に対応するメディアン画像のプロファイルを参照し、網膜層の各領域の境界を検出する。すなわち、Ｌ１（ＩＬＭ）、Ｌ２（ＮＦＬとＧＣＬとの境界）、Ｌ３（ＩＳＯＳ）、Ｌ４（ＲＰＥ）、Ｌ５（ＢＭ）等の検出が行なわれる。なお、本実施形態においては、解析の対象となる層をＲＰＥとして説明する。

［Ｓ２０２］
画像処理装置３０は、判定部４３において、網膜層の３次元の形状解析を行なうか否かの判定を行なう。より具体的には、近視眼解析用の撮影モードで撮影が行なわれた場合、網膜層の３次元の形状解析を行なうと判定する。また、近視眼解析モードを用いずに撮影が行なわれた場合（非近視眼解析用の撮影モード）は、当該３次元の形状解析を行なわないと判定する。この場合、検出部４１の検出結果に基づく解析が行なわれる（３次元の形状データを用いない解析）。なお、近視眼解析用の撮影モードであっても、断層画像に基づいて当該断層画像内に黄斑部が写っているか否かを判定し、写っていなければ（例えば、視神経乳頭部のみが写ってる）、網膜層の（３次元の）形状解析を行なわないと判定するようにしても良い。

ここで、本実施形態においては、断層画像撮影装置２０と画像処理装置３０とが一体となった画像処理システム１０において、撮影から解析までの一連の処理を行なうものとして説明を行なうが、これに限られない。すなわち、画像処理装置３０は、３次元形状を取得するための走査パターンを用いて黄斑部が撮影された断層画像が入力された場合に、網膜層の形状解析を行なえれば良く、必ずしも、このような一体となったシステムである必要はない。そのため、断層画像撮影装置２０以外で撮影された断層画像についても、撮影時の情報に基づいて、網膜層の形状解析を行なうことができる。但し、形状解析を必要としない場合は、当該形状解析処理を行なわないようにすれば良い。

［Ｓ２０３］
画像処理装置３０は、形状データ作成部４６において、３次元の形状データを作成する。３次元の形状データは、Ｓ２０１の処理の網膜層の検出結果に基づいて形状解析を行なうために作成される。

ここで、例えば、撮影時の走査パターンがラスタスキャンであれば、複数の隣接する断層画像同士の位置合わせが行なわれる。断層画像の位置合わせでは、例えば、２つの断層画像の類似度を表す評価関数を事前に定義しておき、この評価関数の値が最も良くなるように断層画像を変形する。

評価関数としては、例えば、画素値で評価する方法（例えば、相関係数を用いて評価を行なう方法）が挙げられる。また、画像の変形処理としては、アフィン変換を用いて並進や回転を行なう処理が挙げられる。形状データ作成部４６においては、複数の断層画像同士の位置合わせ処理が終わった後、形状解析の対象となる層の３次元の形状データを作成する。３次元の形状データは、例えば、５１２×５１２×５００のボクセルデータを用意し、網膜層を検出した層データの座標値に対応する箇所にラベルを割り当てることで作成すれば良い。

また、撮影時の走査パターンが放射状スキャンであれば、この場合も同様に、形状データ作成部４６において、断層画像同士の位置合わせを行なった後、３次元の形状データを作成する。但し、この場合、隣接する断層画像の中心付近の領域情報のみを用いて深度方向（図２の断層画像（表示領域６１）でのＺ方向）の位置合わせを行なう。これは、放射状スキャンの場合、隣接する断層画像であっても、断層画像の中心付近と比べてその両端は、情報が疎であり、形状の変化が大きいので位置合わせの情報として使用しない方が良いためである。位置合わせの方法としては、上述した方法を用いれば良い。形状データ作成部４６においては、位置合わせ処理が終わった後、形状解析の対象となる層の３次元の形状データを作成する。

放射状スキャンの場合、例えば、５１２×５１２×５００のボクセルデータを用意し、まず、各断層画像の形状解析の対象となる層データを均等に円形に回転して展開する。その後、展開した層データにおいて、円周方向に隣接する形状データ間で補間処理を行なう。補間処理により、撮影していない箇所の形状データを作成する。補間処理の方法としては、線形補間、非線形補間などの処理を適用すれば良い。３次元の形状データは、網膜層を検出した層データの間を補間した座標値に対応する箇所にラベルを割り当てることで作成すれば良い。

なお、ここで示したボクセルデータの数値は、あくまで一例であり、撮影時のＡスキャン本数や実行する装置のメモリ状況に応じて適宜変更すれば良い。大きなボクセルデータの場合、解像度が高くなるため、形状データを正確に表現することができるが、実行速度が遅くなりメモリも大量に消費する。一方、小さなボクセルデータの場合、解像度は低くなるが、実行速度が速くなりメモリの消費量も少なくなる。

［Ｓ２０４］
画像処理装置３０は、網膜層解析部４４において、網膜層の３次元の形状解析を行なう。ここでは、形状解析の方法として、網膜層の面積と体積とを計測する手法について例を挙げて説明する。

この処理について、図４を用いて説明する。図４には、３次元の形状データ（ＲＰＥ）と、計測面（ＭＳ）と、面積（Ａｒｅａ）と、体積（Ｖｏｌｕｍｅ）とが模式的に示されている。

まず、面積計測について説明する。Ｓ２０３の処理で作成したＲＰＥの３次元の形状データにおいて、Ｚ方向（光軸方向）で最も深部に位置する層データの場所にフラット（平面）な計測面（ＭＳ）を用意する。そして、そこから一定の間隔ずつ浅い方向（Ｚ軸の原点方向）に計測面（ＭＳ）を移動させる。層の深部から浅い方向に計測面（ＭＳ）を移動させると、ＲＰＥの境界線を計測面（ＭＳ）が横断する。

面積（Ａｒｅａ）は、その計測面（ＭＳ）とＲＰＥとの境界線で囲まれる内部の平面領域を計測することにより得られる。より具体的には、面積（Ａｒｅａ）は、計測面（ＭＳ）とＲＰＥとの境界線との交差領域の面積を計測することにより得られる。このように面積（Ａｒｅａ）とは、３次元の網膜層形状データの断面積である。基準となる部位の位置の断面積を計測する場合、網膜の湾曲がきつい形状であれば、面積は小さくなり、網膜の湾曲が緩やかであれば、面積は大きくなる。

また、体積（Ｖｏｌｕｍｅ）は、面積（Ａｒｅａ）を計測する際に用いた計測面（ＭＳ）を利用して、計測面（ＭＳ）とＲＰＥとの境界線で囲まれる内部の全領域を計測することにより得られる。基準となる部位の位置からＺ方向に沿って下の体積を計測する場合、網膜の湾曲がきつい形状であれば、体積は大きくなり、網膜の湾曲が緩やかであれば、体積は小さくなる。ここで、基準位置は、部位を基準として、例えば、ブルッフ膜オープニングの位置とすれば良い。或いは、ＲＰＥの最深部の位置から１００μｍ、５００μｍなどと一定の高さを基準としても良い。なお、面積や体積を計測する際に、内部に含まれるボクセル数をカウントする場合、１ボクセルあたりの物理サイズを乗じて計算する。

［Ｓ２０５］
形状解析後、画像処理装置３０は、解析結果作成部４５において、当該３次元の形状解析結果に基づいて解析結果（例えば、マップ、グラフ、数値情報）を作成する。

ここで、３次元の形状解析結果として等高線マップを作成する場合について説明する。等高線マップは、面積及び体積を計測した結果を表示する際に使用される。

図５（ａ）は、等高線マップ７１の一例を示す。ここで、符号７１は、等高線マップ全体を示し、符号７２は、一定の間隔で引いた等高線を示し、符号７３は、３次元の網膜層形状データのＺ方向で最も深部に位置する箇所を示している。

等高線マップ７１を作成する場合、等高線マップ用のルックアップテーブルを用意し、当該テーブルを参照して体積に応じて色分けする。ここで、等高線マップ用のルックアップテーブルは、例えば、体積に応じて用意すれば良い。それにより、マップ上においては、形状と体積の変化を一度に把握することができる。

より具体的には、図５（ａ）の符号７３を底として等高線７２が一定の間隔で描画されており、マップの色が体積に応じて設定されている。そのため、Ｚ方向で最も深部に位置する箇所から計測面の高さ（深さ）を変更した場合、体積がどのように増加するかを把握することができる。具体的には、１ｍｍ^３の場合は青とし、２ｍｍ^３の場合は黄色として設定しておくと、Ｚ方向で最も深部に位置する箇所から計測面の高さ（深さ）が１００μｍで青なのか、３００μｍで青なのかにより、操作者は、形状と体積との関係を把握することができる。従って、操作者は、マップの最外側の輪郭の色を確認することで、計測対象の網膜層の形状を計測した時の全体の体積の値を把握することができる。

また、操作者は、内部の等高線近傍の色を確認することで高さ（深さ）に対応した体積の値を把握することができる。なお、等高線マップの色を設定する際に用いられるルックアップテーブルは、体積に限らず、面積に応じて設けられても良いし、また、計測面（ＭＳ）から網膜層データのＺ方向で最も深部に位置する箇所７３までの高さ（深さ）に応じて設けられていても良い。図示はしないが、等高線の高さ（深さ）が分かるようにそれぞれの等高線上に数値を一緒に表示しても良い。等高線を表現する距離の間隔としては、例えば、Ｚ方向に沿って１００μｍ間隔としておく。なお、等高線マップは、カラーであっても良いし、グレースケールであっても良いが、カラーの方が視認性は良くなる。

また、等高線マップは、Ｚ方向で最も深部に位置する箇所から計測面（ＭＳ）までの高さ（深さ）に応じて外形のサイズが変わる。図５（ｂ）は、計測面（ＭＳ）の高さ（深さ）を変更した場合の等高線マップ７４の一例を示している。

次に、３次元の形状解析として、網膜層の曲率を計測する場合について図２に示す断層画像（表示領域６１）を用いて説明する。ここでは、横軸をｘ座標とするとともに縦軸をｚ座標とし、解析の対象となる層（ＲＰＥ）の境界線の曲率を計算する場合について説明する。曲率κは境界線の各点において、「式１」を計算することで求めることができる。曲率κの符号で上に凸か下に凸かが分かり、数値の大きさで形状の曲がり具合が分かる。そのため、上に凸を＋、下に凸を−とした場合、各断層画像において、曲率の符号が−領域、＋領域、−領域となる場合はＷ形状となる。

［式１］

なお、ここでは、断層画像の境界線で曲率を計算する場合について説明したが、曲率計算はこれに限らず、３次元の形状データから３次元の曲率を計算するようにしても良い。この場合、形状解析後、画像処理装置３０は、解析結果作成部４５において、当該解析結果に基づいて曲率マップを作成する。

図５（ｃ）に曲率マップの一例を示す。この場合、曲率がきつい箇所の色を濃く、曲率が滑らかな箇所の色を薄く表現している。すなわち、曲率に応じて色の濃さが変更されている。なお、曲率マップに設定する色は、曲率値０を基準に正負で色を変えれば良い。それにより、操作者は、網膜形状が滑らかなのか、上に凸形状か下に凸形状なのかを当該マップを見ることで把握することができる。

［Ｓ２０６］
画像処理装置３０は、表示制御部３５において、断層画像と、検出部４１により検出された層（ＲＰＥ）の検出結果と、解析結果作成部４５により作成された各種形状解析の結果（マップ、グラフ、数値情報）とを表示装置５３に表示する。

図６は、図１に示す表示装置５３に表示される断層画像観察画面８０の一例を示す図である。この画面は、断層画像の解析が終わった後に表示される（すなわち、Ｓ２０６の処理によって表示される）。

断層画像観察画面８０は、断層画像の表示領域８１を含む断層画像表示部９１と、眼底画像の表示領域８２を含む眼底画像表示部９４とを具備して構成される。また、断層画像観察画面８０には、第１の解析結果８４を含む第１の解析結果表示部９６、第２の解析結果８５及び８６を含む第２の解析結果表示部９８も設けられる。

ここで、まず、断層画像の表示領域８１を含む断層画像表示部９１の詳細について説明する。断層画像の表示領域８１では、撮影された断層画像に対して、網膜層の各層を検出したセグメンテーション結果（Ｌ１〜Ｌ５）と計測面（ＭＳ）とが重畳して表示されている。断層画像の表示領域８１では、解析対象の網膜層（本実施形態においては、ＲＰＥ（Ｌ４））のセグメンテーション結果が強調して表示されている。

また、断層画像の表示領域８１において、計測面（ＭＳ）と解析対象の網膜層（ＲＰＥ（Ｌ４））とで囲まれる斜線領域８１ａが面積及び体積の計測対象領域となる。このとき、斜線領域８１ａでは、体積計測結果に応じた色が所定の透明度αによって表示される。設定する色は、等高線マップ用のルックアップテーブルと同じものを使用すれば良い。透明度αは、例えば、０．５とする。

コンボボックス９２は、断層画像を表示する際に、ＯＣＴ比で表示するか、１：１比で表示するかを操作者が選択するために設けられる。ここで、ＯＣＴ比とは、撮影時のＡスキャンの本数から得られる横方向（Ｘ方向）の解像度と、縦方向（Ｙ方向）の解像度とで表わした比のことである。１：１比とは、何ｍｍの範囲を何本のＡスキャンで撮影したかで得られる、横方向の１ピクセルあたりの物理サイズと、縦方向の１ピクセルあたりの物理サイズとを揃えた比のことである。

コンボボックス９３は、断層画像を表示する際に、２次元表示（２Ｄ）／３次元表示（３Ｄ）を切り替えるために設けられる。ここで、２Ｄ表示の場合、断層画像の１スライスが表示され、３Ｄ表示の場合、網膜層の境界線データから作成された網膜層の３次元形状が表示される。

具体的には、コンボボックス９３において３Ｄ表示が選択されると、断層画像の表示領域８１には、図７（ａ）〜図７（ｃ）に示す断層画像が表示される。

図７（ａ）は、ＯＣＴ比でＲＰＥを３Ｄ表示した場合の態様を示している。この場合、計測面（ＭＳ）も同時に３Ｄ表示されている。

ここで、断層画像の表示領域８１の下側には、網膜の各層に対応したチェックボックス１０１〜１０４が表示されている。具体的には、ＩＬＭ、ＰＲＥ、ＢＭ、ＭＳに対応したチェックボックスが表示されており、操作者は、これらチェックボックスを用いることで、各層の表示／非表示の切り替えを指示できる。

計測面（ＭＳ）を面で表現する場合、その透明度αは、０より大きく、且つ１よりも小さい値とする。透明度を１にした状態で、網膜層形状と計測面とを重ねると、蓋をしたような状態となり、網膜層の３次元形状が上側からでは分からなくなってしまうためである。或いは、計測面（ＭＳ）を面ではなく格子状にして表現しても良い。格子状の場合、計測面（ＭＳ）の透明度αは１としても良い。計測面（ＭＳ）の色は、等高線マップのルックアップテーブルを参照し、計測面（ＭＳ）がある位置の計測値（面積又は体積）に応じた色を選べば良い。

ここで、操作者は、入力装置５４を介して、計測面（ＭＳ）の位置を移動させることができる。そのため、計測面（ＭＳ）の位置が変更されると、画像処理装置３０は、当該変更に連動して、図５（ａ）及び図５（ｂ）で説明したように、等高線マップ形状を変化させる。

テキストボックス１０５は、数値を指定するための項目である。操作者は、入力装置５４を介して、Ｚ方向で最も深部に位置する箇所からの計測面（ＭＳ）の高さ（深さ）を示す数値をこのテキストボックス１０５に入力する。例えば、１００μｍや３００μｍといった数値が入力された場合、その場所に計測面（ＭＳ）が移動し、等高線マップもそれに合わせて変化する。これにより、操作者は、３次元形状での位置とその時の等高線マップとを同時に把握することができるとともに、その時の体積値、面積値についても把握することができる。

その他の例として、テキストボックス１０５に体積の値を入力できるようにしても良い。この場合、１ｍｍ^３や２ｍｍ^３といった数値が入力された場合、操作者は、当該体積に対応する３次元形状での計測面（ＭＳ）の位置とその時の等高線マップとを同時に把握することができる。また更に、操作者は、その時のＺ方向で最も深部に位置する箇所からの計測面（ＭＳ）の高さ（深さ）も把握できる。

続いて、図７（ｂ）は、計測面（ＭＳ）を非表示にした場合の表示態様を示している。それ以外の表示項目については、図７（ａ）と同様となる。図７（ｂ）の場合、チェックボックス１０２が選択されており、ＲＰＥのみが表示されている。このように眼底の３次元形状のみを表示することもできる。

ここで、図７（ｂ）では、ＯＣＴ比でＲＰＥの３Ｄ表示を行なっているのに対して、図７（ｃ）は、１：１比でＲＰＥの３Ｄ表示を行なっている場合の態様を示している。すなわち、コンボボックス９２において、１：１比が選択されている。

なお、本実施形態では、ＲＰＥを解析対象として説明しているため、コンボボックス９３で２Ｄ／３Ｄ表示を切り替えた場合、ＲＰＥ形状を２Ｄ／３Ｄ表示することになるが、これに限られない。例えば、解析対象をブルッフ膜（ＢＭ）とした場合、ブルッフ膜（ＢＭ）形状が２Ｄ／３Ｄ表示される。

ここで、断層画像の表示領域８１において、計測面がどの位置を測定しているかを模式的に示す表示を行なうようにしても良い。この場合の表示態様について、図８（ａ）〜図８（ｄ）を用いて説明する。

図８（ａ）は、表示中の断層画像において、計測面がどの位置を測定しているかを示すオブジェクト１１０を断層画像上に重畳して表示した態様を示している。符号ＭＳ’は、オブジェクト１１０上での計測面（模式計測面）を表している。網膜層の３次元形状は、操作者の入力装置５４を介した指示により上下左右の回転が行なわれる。そのため、計測面（ＭＳ）と網膜層との位置関係を操作者が把握できるように、オブジェクト１１０及び模式計測面ＭＳ’によって位置関係の指標を提示する。

なお、操作者が入力装置５４を介して、オブジェクト１１０及び模式計測面ＭＳ’の位置関係を変更した場合には、それに連動して網膜層の３次元形状の表示も変更される。この場合、解析対象となる網膜層の領域も変更されるため、第１の解析結果８４や第２の解析結果８５及び８６も連動して変更される。

ここで、オブジェクト１１０の表示態様をいくつか例示して説明する。例えば、図８（ｂ）及び図８（ｃ）に示すように、各断面位置に対応する断層画像を表示するようにしても良い。この場合、オブジェクト１１０においては、３次元で考えた場合の中心位置を横切る垂直位置及び水平位置の断層画像が表示される。或いは、図８（ｄ）に示すように、上方（superior）、下方（inferior）を示すＳ、Ｉなどの英略語を表示するようにしても良い。

次に、図６に示す眼底画像の表示領域８２を含む眼底画像表示部９４の詳細について説明する。眼底画像の表示領域８２では、撮影位置とその走査パターンマーク８３とが眼底画像に対して重畳されている。また、眼底画像表示部９４には、眼底画像の表示形式を切り替えるためのコンボボックス９５が設けられている。この場合、眼底画像としてＳＬＯ画像が表示されている。

ここで、図９（ａ）及び図９（ｂ）を用いて、コンボボックス９５から眼底画像の表示形式の切り替えを行なった場合について説明する。ここでは、眼底画像の表示形式として、“ＳＬＯ画像＋マップ”を同時に表示する場合と、“眼底写真（第２の眼底画像）＋ＳＬＯ画像（第１の眼底画像）＋マップ”を同時に表示する場合とを例に挙げて説明する。なお、ＳＬＯ画像（第１の眼底画像）は、断層画像と同時に撮影が行なわれた２次元の眼底画像であれば良く、例えば、断層画像を深度方向に積算することにより生成された積算画像であっても良い。また、眼底写真（第２の眼底画像）は、断層画像と異なるタイミングで撮影された２次元の眼底画像であれば良く、例えば、造影検査画像等が挙げられる。

図９（ａ）は、“ＳＬＯ画像＋マップ”がコンボボックス９５から選択された場合の表示態様を示している。符号２０１はＳＬＯ画像を示し、符号２００はマップを示している。ＳＬＯ画像２０１とマップ２００との位置合わせは、図１で説明した位置合わせ部４７により行なわれる。これらの位置合わせは、撮影時の固視灯の位置と走査範囲とに基づいてマップの位置及びサイズを設定することにより行なわれる。

ＳＬＯ画像２０１に対してマップ２００を重畳した結果を眼底画像として表示する場合、ＳＬＯ画像２０１の透明度αを１として表示し、マップ２００の透明度αを１より小さくし、例えば、０．５として表示する。これら透明度αのパラメータは、操作者が対象データを初めて選択した時に設定するパラメータである。これらのパラメータは、操作者が、入力装置５４を介して適宜変更することができる。操作者が変更したパラメータは、例えば、外部記憶装置５２に記憶される。次回以降同じ対象データを開いた時には、前回の操作者が設定したパラメータに従って表示が行なわれる。

続いて、図９（ｂ）は、“眼底写真＋ＳＬＯ画像＋マップ”がコンボボックス９５から選択された場合の表示態様を示している。符号２０２は、眼底写真（第２の眼底画像）を示している。

この場合、眼底写真２０２に対してマップ２００を重畳表示するために、ＳＬＯ画像２０１が利用されている。その理由として、眼底写真２０２は、断層画像と異なるタイミングで撮影されたものであり、マップ２００だけでは撮影位置や撮影範囲が分からないためである。そこで、ＳＬＯ画像２０１を用いることで眼底写真２０２とマップ２００との位置合わせを行なうことができる。眼底写真２０２とＳＬＯ画像２０１との位置合わせは、図１で説明した位置合わせ部４７により行なわれる。

位置合わせの方法は、例えば、血管特徴を用いて行なえば良い。血管の検出方法としては、血管は細い線状構造を有しているため、まず、線状構造を強調するフィルタを用いて血管を抽出する。線状構造を強調するフィルタとしては、線分を構造要素としたときに構造要素内での画像濃度値の平均値と構造要素を囲む局所領域内での平均値との差を計算するフィルタを利用すれば良い。勿論、これに限られず、Ｓｏｂｅｌフィルタのような差分フィルタを用いても良い。また、濃度値画像の画素毎にヘッセ行列の固有値を計算し、その結果として得られる２つの固有値の組み合わせから線分状の領域を抽出しても良い。更には、線分を単純に構造要素とするトップハット演算を用いても良い。位置合わせ部４７においては、これらの方法によって検出した血管位置情報を用いて、眼底写真２０２とＳＬＯ画像２０１との位置合わせを行なう。

ＳＬＯ画像２０１とマップ２００とは、上述した方法で位置合わせを行なうことができるので、その結果として眼底写真２０２とマップ２００との位置合わせも行なうことができる。眼底写真２０２に対してマップ２００を重畳した結果を眼底画像表示部９４の表示領域８２に表示する場合には、眼底写真２０２の透明度αを１として表示する。また、ＳＬＯ画像２０１の透明度α及びマップ２００の透明度αを１より小さくし、例えば、０．５として表示する。勿論、ＳＬＯ画像２０１とマップ２００との透明度αの値は同じである必要はなく、例えば、ＳＬＯ画像２０１の透明度αの値を０としても良い。

なお、被検眼が疾患眼である場合、位置合わせ部４７による位置合わせ処理が失敗する可能性がある。位置合わせの失敗の判定は、画像間類似度を計算する場合には、最大の類似度が閾値以上とならない。或いは、閾値以上となったとしても、解剖学的に不自然な箇所で位置合わせ処理が終了する場合、失敗したと判定できる。その場合、ＳＬＯ画像２０１及びマップ２００は、事前に設定された初期パラメータの位置（例えば、画像の中心）及び初期透過度で表示すれば良い。そして、位置合わせ処理に失敗したメッセージを表示して、操作者による入力装置５４を介した位置修正処理を促す。

ここで、位置修正や透過度αのパラメータ変更を操作者が行なう場合、操作者が、ＳＬＯ画像２０１に対して移動、拡大・縮小、回転操作を行なうと、ＳＬＯ画像２０１上のマップ２００も同時に、移動、拡大・縮小、回転する。すなわち、ＳＬＯ画像２０１及びマップ２００は、１枚の同じ画像として動作する。但し、透明度αに関しては、ＳＬＯ画像２０１及びマップ２００は、それぞれ独立に設定が行なわれる。操作者が、入力装置５４を介して変更したパラメータは、外部記憶装置５２に記憶され、次回以降は当該設定されたパラメータに従って表示が行なわれる。

このように眼底画像表示部９４における眼底画像の表示領域８２には、２次元の眼底画像の他、眼底画像上にマップ等が重畳された表示が行なわれる。なお、図９（ａ）や図９（ｂ）では、等高線マップを眼底画像上の対応する位置に関連付けて重畳表示する場合について説明したが、これに限られない。すなわち、曲率マップや層厚マップなどを表示しても良い。

また、解析対象の層を変更する場合には、操作者が、入力装置５４を介して、例えば、図２に示す表示領域６１の断層画像上からセグメンテーション結果（Ｌ１〜Ｌ５）を選択することで行なえば良い。解析対象の層を切り替えた場合、画像処理装置３０は、表示制御部３５において、それまでに強調表示されていたセグメンテーション結果の層を通常表示に戻し、新たな解析対象の層を強調表示させる。それにより、任意の層の解析結果を表示することができる。

次に、図６に示す第１の解析結果８４を含む第１の解析結果表示部９６の詳細について説明する。

第１の解析結果８４には、解析結果作成部４５により作成された形状解析マップが表示される。また、コンボボックス９７では、第１の解析結果８４のマップの種類を選択することができる。この場合、第１の解析結果８４として示されている形状解析マップは、等高線マップとなっている。

第１の解析結果８４として示されている形状解析マップの種類や、図９（ａ）や図９（ｂ）で説明した眼底画像上に重畳表示された形状解析マップの種類は、コンボボックス９７から種類を指定することにより連動して変更できる。更に、以下で説明を行なう第２の解析結果表示部９８の表示も連動して変更される。

ここで、図１０に解析結果として曲率の結果を表示させた場合の断層画像観察画面８０の一例について説明する。

解析結果の表示を、面積及び体積の解析結果から曲率解析結果に切り替えた場合、断層画像表示部９１、眼底画像表示部９４、第１の解析結果表示部９６、第２の解析結果表示部９８に表示される結果は、図１０に示す内容となる。

具体的には、断層画像の表示領域８１には、撮影された断層画像に対して、当該断層画像上に網膜層の各層を検出したセグメンテーション結果（Ｌ１〜Ｌ５）が重畳表示される。また、第１の解析結果８４としては曲率マップが表示され、第２の解析結果８８としては曲率グラフが表示される。眼底画像の表示領域８２には、ＳＬＯ画像（第１の眼底画像）２０１と、眼底写真（第２の眼底画像）２０２と、曲率マップ２０３とを重畳表示した画像が表示される。

次に、図６に示す第２の解析結果８５及び８６を含む第２の解析結果表示部９８の詳細について説明する。

第２の解析結果８５では、解析結果作成部４５により作成された形状解析グラフが表示される。この場合、面積及び体積を計測したグラフが示されており、横軸が高さ（深さ）を示し、縦軸が体積を示している。符号８７に示す実線は、体積を示している。

第２の解析結果８６では、形状解析結果が表として表示されている。表には、ある一定の基準値の高さ（例えば、１００μｍ、５００μｍ等）の時の面積及び体積が表示されるとともに、ある部位を基準とした時の高さ（ブルッフ膜オープニングの位置）に対応する面積及び体積が表示される。

ここで、図１１を用いて、面積及び体積の結果を一つのグラフ（第２の解析結果８５）に表示した場合について説明する。この場合、横軸が高さ（深さ）を示し、グラフ左側の縦軸が体積を示し、グラフ右側の縦軸が面積を示している。破線８９は面積を示すグラフであり、実線８７は体積を示すグラフである。

［Ｓ２０７］
画像処理装置３０は、網膜層解析部４４において、網膜層の検出結果に基づいて形状解析を行なう。この解析処理では、３次元の形状データ等を作成せずに、網膜層の検出結果を用いた解析が行なわれる。例えば、層厚の解析等が行なわれる。

［Ｓ２０８］
画像処理装置３０は、解析結果作成部４５において、当該解析結果に基づいて解析結果（例えば、マップ、グラフ、数値情報）を作成する。

［Ｓ２０９］
画像処理装置３０は、断層画像と、検出部４１により検出された層（ＲＰＥ、ＩＬＭ等）の検出結果と、解析結果作成部４５により作成された解析結果（マップ、グラフ、数値情報）とを表示装置５３に表示する。

図１２は、Ｓ２０９の処理によって、表示装置５３に表示される断層画像観察画面８０の一例を示す図である。

この場合、第１の解析結果３０２として、層厚マップが示されており、第２の解析結果３０１として、層厚グラフが示されている。すなわち、第１の解析結果３０２及び第２の解析結果３０１として、網膜層の検出結果を用いた解析結果が表示されている。

この場合にも、上記Ｓ２０６の処理と同様に、操作者が、入力装置５４を介して変更したパラメータは、外部記憶装置５２に記憶され、次回以降は当該設定されたパラメータに従って表示が行なわれる。

以上説明したように実施形態１によれば、近視眼解析用の撮影モードを含む複数の撮影モードを設け、断層画像の撮影時の撮影モードに応じて解析処理を切り替えて実行することができる。

より具体的には、複数の撮影モードとして、近視眼解析用の撮影モードと、非近視眼解析用の撮影モードとが少なくとも設けられる。近視眼解析用の撮影モードでは、３次元の形状解析（近視眼における黄斑部の解析）を行なうのに適した断層画像の撮影が行なわれる。

そして、近視眼解析用の撮影モードで撮影された断層画像に対しては、３次元の形状データを作成し、当該形状データに基づいて３次元の形状解析を実行する。それ以外の撮影モードで撮影された断層画像に対しては、３次元の形状データを作成せずに、網膜層の検出結果に基づく解析処理を実行する。

（実施形態２）
次に、実施形態２について説明する。実施形態２では、断層画像の２Ｄ表示と３Ｄ表示とを同時に表示する場合について説明する。より具体的には、実施形態２では、断層画像と、３次元の形状データと、解析結果とを並べて表示する場合について説明する。なお、実施形態２では、断層画像として、放射状スキャンで撮影された断層画像を例に挙げて説明する。

図１３を用いて、実施形態２に係わる断層画像観察画面４００の一例について説明する。

断層画像観察画面４００には、２次元の断層画像の表示領域４１１を含む第１の断層画像表示部４１０と、３次元の断層画像の表示領域４３１を含む第２の断層画像表示部４３０とが設けられる。また更に、断層画像観察画面４００には、第１の解析結果４４２を含む第１の解析結果表示部４４０、第２の解析結果４２１及び４２２を含む第２の解析結果表示部４２０も設けられる。

ここでは、第１の解析結果表示部４４０及び第２の解析結果表示部４２０は、上述した実施形態１を説明した図６や図９（ａ）と同様であるため、その説明については省略する。また、第２の断層画像表示部４３０も、実施形態１を説明した図７（ａ）〜図７（ｃ）と同様であるため、その説明については省略する。ここでは、第１の断層画像表示部４１０について重点的に説明する。

第１の断層画像表示部４１０には、２次元の断層画像の表示領域４１１の他、断層画像の視点位置（スライス位置）を変更するためのスライダーバー４１２と、断層画像のスライス番号を表示するエリア４１３とが設けられる。

ここで、図１４を用いて、３次元形状での視点方向５０３と、図１３に示す断層画像（表示領域４１１）との関係について説明する。符号５００は、断層画像の３次元形状を上から深部方向（Ｚ方向）に見た眼底の概要を示しており、符号５０１に示す放射状の直線は、断層画像の撮影スライス位置を示している。符号５０２に示す破線は、撮影スライス位置５０１の中で表示中の断層画像（表示領域４１１）と対応しているスライス位置を示しており、視点方向５０３に対して垂直になる。すなわち、破線５０２に示す位置の断層画像が、断層画像の表示領域４１１に表示される断層画像となる。

次に、操作者が、入力装置５４を介して、図１３で説明した断層画像観察画面４００を操作した場合について説明する。

操作者が、入力装置５４を介して、スライダーバー４１２の位置を中心から端に移動させたとする。そうすると、その操作に連動して、画像処理装置３０は、表示領域４３１に表示中の３次元の断層画像（３次元形状データ）の視点位置を変更する。また、表示領域４１１に表示中の２次元の断層画像のスライス位置をスライダーバー４１２で指定された断層画像と対応する視点位置に変更する。

なお、スライダーバー４１２が操作された場合、当該バーは、離散的に変化するのではなく、連続的にその位置が変更される。そのため、表示中の断層画像の３次元形状データの位置も連続的に変化する。具体的には、画像処理装置３０は、表示中の断層画像の３次元形状データが垂直方向の中心を基準として回転するように表示する。

図１５は、スライダーバー４１２の位置を端に移動させた場合の表示領域４１１における２次元の断層画像と、表示領域４３１の３次元の断層画像（３次元形状データ）との画面表示の一例を示している。この場合、３次元形状での視点方向５０７と、２次元の断層画像（表示領域４１１）との関係は図１６に示す関係となる。

スライダーバー４１２の位置が操作されると、画像処理装置３０は、表示領域４１１に表示中の２次元の断層画像のスライス位置を変更するとともに、そのスライス位置に対応する３次元の断層画像（３次元形状データ）の視点位置を変更する。なお、２次元の断層画像のスライス位置を変更するためのスライダーバー４１２の操作に伴って、３次元形状データの視点位置を変更させるだけでなく、その逆も可能である。つまり、操作者が、入力装置５４を介して、３次元形状データの視点位置を変更させた場合、それに対応して２次元の断層画像のスライス位置と、スライダーバー４１２の位置とを変更させるようにしても良い。

以上説明したように実施形態２によれば、２次元の断層画像と３次元の断層画像とを同時に表示し、操作者からの操作に応じて、それらの表示態様を連動させて変更させることができる。

（実施形態３）
次に、実施形態３について説明する。実施形態３では、網膜層の形状の解析結果と、データベースに保持されている統計データベースとの両方を比較可能に表示する場合について説明する。

図１７は、実施形態３に係わる画像処理システム１０の構成の一例を示す図である。ここでは、実施形態１を説明した図１の構成と相違する点について重点的に説明する。

ここで、実施形態３に係わる画像処理装置３０の画像処理部３３には、実施形態１の構成に加えて、特徴量取得部６０１と、比較演算部６０２とが新たに設けられる。

特徴量取得部６０１は、３次元の形状解析に基づいて、網膜層の形状特徴を示す特徴量（形状特徴量）の算出を行なう。

比較演算部６０２は、形状特徴量（第１の特徴量）と、外部記憶装置５２に記憶される統計データベース５５に保持された形状特徴量（第２の特徴量）とを比較し、その比較結果を出力する。なお、統計データベース５５とは、多数眼のデータから作成されるものであり、人種別、年齢別のデータを統合して作成される。すなわち、統計データベース５５には、複数の被検眼から得られた統計的な特徴量が保持される。なお、眼科の分野においては、左右眼別、眼軸長別など、眼特有のパラメータで分類されていても良い。

次に、図１８を用いて、実施形態３に係わる画像処理装置３０の処理の流れの一例について説明する。なお、断層画像の撮影時の全体的な処理は、実施形態１の図３（ａ）で説明した内容と同様となるので、ここでは、その説明は省略し、図３（ｂ）の処理（網膜層の検出／解析処理）との相違点について重点的に説明する。

［Ｓ３０１〜Ｓ３０４］
網膜層の検出／解析処理が開始すると、画像処理装置３０は、実施形態１を説明した図３（ｂ）のＳ２０１〜Ｓ２０４の処理を実施する。これにより、断層画像から網膜層が検出されるとともに、当該検出された網膜層の３次元の形状解析が行なわれる。

［Ｓ３０５］
画像処理装置３０は、特徴量取得部６０１において、Ｓ３０４の処理で行なわれた３次元の形状解析に基づいて、形状特徴量の算出を行なう。

ここで、網膜層の形状特徴としては、例えば、網膜層の面積及び体積の特徴量、網膜層の円形の（凹凸の）度合いを示す特徴量、網膜層の形状の曲率の特徴量、等が挙げられる。

まず、網膜層の面積及び体積の特徴量の求め方について説明する。網膜層の面積及び体積の特徴量は、例えば、網膜層の面積及び体積の増減率（典型的には、増加率）を算出することにより求められる。網膜層の面積及び体積は、上述した実施形態１と同様の方法で求めれば良い。

面積及び体積の増減率を算出する場合、特徴量取得部６０１は、網膜層解析部４４による３次元の形状解析の解析結果に基づいて、ＲＰＥの最深部の位置から１００μｍ、２００μｍ、３００μｍなどと一定の高さ（深さ）にあるときの面積及び体積をそれぞれ求める。そして、その結果に基づいて面積及び体積の増減率を網膜層の形状特徴として求める。

また、網膜層の円形の（凹凸の）度合いを示す特徴量は、例えば、（最終的に）面積計測の対象となる網膜層の断面の形状から算出する。ここで、図１９（ａ）及び図１９（ｂ）を用いて、網膜層の円形の（凹凸の）度合いを示す特徴量の算出方法について具体例を挙げて説明する。

図１９（ａ）には、水平方向最大弦長（ＣＨＯＲＤ Ｈ）、垂直方向最大弦長（ＣＨＯＲＤ Ｖ）、絶対最大長（ＭＡＸＩＭＵＭ ＬＥＮＧＴＨ）が示されている。図１９（ｂ）には、面積（ＡＲＥＡ）と周囲長（ＰＥＲＩＭＥＴＥＲ）とが示されており、これらの値から円形度Ｃ_Ｌを計算することができる。円形度の式を「式２」に示す。

［式２］

なお、Ａは面積（ＡＲＥＡ）を示し、Ｌは周囲長（ＰＥＲＩＭＥＴＥＲ）を示す。円形度の逆数を計算すると、凹凸の度合いを示す特徴量となる。勿論、円形の凹凸の度合いを示す特徴量の求め方は、これに限られず、例えば、重心回りのモーメント特徴などを用いても良い。

続いて、網膜層の形状の曲率の特徴量の求め方について説明する。曲率の特徴量を算出する場合、特徴量取得部６０１は、網膜層のそれぞれの箇所において曲率を計算し、全個所における平均、分散、標準偏差の計算を行なう。

また、曲率値を絶対値ではなく、符号付き曲率値で求める場合、正負で形状の凹凸が分かる。そのため、符号付き曲率値で求める場合には、曲率値として正の曲率値、負の曲率値の数をそれぞれカウントしておく。正負の数のカウントにおいて、完全に直線の場合、曲率値は０となるが、自動又は手動による網膜層の検出では多少の誤差が含まれる。そのため、曲率値０付近の値は、カウントから除き、曲率値が正負に大きい値、すなわち、特徴が大きく出ている個所だけをカウントするようにしても良い。

［Ｓ３０６］
画像処理装置３０は、比較演算部６０２において、Ｓ３０５の処理で求められた特徴量と、外部記憶装置５２に記憶されている統計データベース５５の特徴量とを比較する。なお、統計データベース５５では、正常データの９５％範囲を正常、４％範囲をボーダーライン、残り１％範囲を異常などした情報が保持される。

比較演算部６０２では、Ｓ３０４及びＳ３０５の処理で求められた解析結果や特徴量の値が、統計データベース５５に保持された特徴量に対して、どの範囲に位置するかを比較する。例えば、解析結果の面積や体積についてのみ比較を行なっても良いし、面積や体積の値と特徴量の値とを２軸として、その中でどこに位置するのかを比較するようにしても良い。また更に、特徴量が曲率の場合には、曲率値の平均値とその分散値とで２軸を作成し、その軸の中でどの範囲に含まれるかを比較するようにしても良い。

この点について更に説明する。図２０（ａ）に示すグラフでは、横軸が曲率の平均値を示し、縦軸が曲率の分散値を示す。符号７０１に示す菱形は、近視眼において後部ぶどう腫の症状がない場合の網膜層の曲率値を示し、符号７０２に示す三角は、近視眼において軽度な後部ぶどう腫の症状がある場合の網膜層の曲率値を示している。また、符号７０３に示す四角は、近視眼において後部ぶどう腫の症状がある場合の網膜層の曲率値を示す。

なお、後部ぶどう腫とは、眼球が後方に突出したように形状が変形する状態のことをいう。プロット点が複数存在するのは、症例の数であり、この場合、数十例分をプロットした例を示している。また、ここでは、Ｓ３０５の処理により得られた形状特徴量についての説明は省略するが、この形状特徴量は、統計データベース５５に保持された特徴量と識別可能な形式（例えば、表示形式（丸形状等）が変更される）で当該グラフ上にプロットされる。すなわち、Ｓ３０５の処理により得られた形状特徴量は、統計データベース５５に保持された特徴量と比較可能に表示される。

ここで、領域７１１は、後部ぶどう腫の症状がない曲率値領域を示し、領域７１２は、軽度な後部ぶどう腫の症状がある曲率値領域を示し、領域７１３は、後部ぶどう腫の症状がある曲率値領域を示している。従って、例えば、対象としている被検眼の曲率値を計測した際に、当該被検眼の曲率値が領域７１１に位置すれば後部ぶどう腫がなく、領域７１３に位置すれば後部ぶどう腫があると判定できる。

このように本実施形態においては、被検眼が疾病である可能性を段階的に示す複数の領域（すなわち、疾病の段階的に示す被検眼のカテゴリを示す表示形態）で囲って、曲率値（特徴量）をグラフ上にプロットして表示する。また、どの領域に属するかによって、曲率値を示すオブジェクトの表示形式（菱形、三角、四角）を変更して表示する。なお、ここでは説明を簡単にするため、３つの領域に分類する場合について説明したが、これに限らない。例えば、それぞれの領域を複数に分割して、合計で１０個程度に分類し、後部ぶどう腫の程度の大きさを判定できるようにしても良い。

［Ｓ３０７］
画像処理装置３０は、解析結果作成部４５において、Ｓ３０５の処理の結果に基づいて解析結果を作成する。解析結果は、例えば、対象としている被検眼の解析値が統計データベース５５のどこに位置するかが把握できる形式で作成される。

解析結果の一例としては、例えば、図２０（ａ）で説明したように、統計データベース５５を２次元で表現した場合に、当該被検眼の曲率値がどこに位置するかをプロットしたグラフが挙げられる。また、図２０（ｂ）に示すように、統計データベース５５の特徴を１次元で表現したカラーバー（棒グラフ）を作成し、そのカラーバー（棒グラフ）のどこに位置するのかを示す解析結果を作成しても良い。

また更に、曲率値のように、各箇所において一つの値が対応する場合、その箇所における統計データベース５５に保持されたいずれかの曲率値と比較を行ない、各箇所が統計データベース５５のどこに位置するかが分かる２次元のマップを作成しても良い。なお、２次元のマップにおいては、各箇所での比較ではなく、一定の範囲を平均化した値を統計データベース５５と比較するようにしても良い。

［Ｓ３０８］
画像処理装置３０は、表示制御部３５において、断層画像と、検出部４１により検出された層の検出結果と、解析結果作成部４５により作成された解析結果（マップ、グラフ、数値情報）とを表示装置５３に表示する。なお、解析結果としては、Ｓ３０７の処理で作成されたマップやグラフ等を表示すれば良い。具体的には、図２１に示す断層画像観察画面８０を表示する。

この場合、第１の解析結果８４として、統計データベース５５の特徴量（曲率）との比較結果として、網膜層の曲率マップが表示されている。この曲率マップでは、一定の範囲を平均化した値を統計データベース５５と比較した結果が示されている。色の濃淡が濃い場所ほど、曲率値が統計データベース５５の正常値から外れており、濃淡が淡いほど、曲率値が統計データベース５５の正常値の範囲内であることを示している。

第２の解析結果表示部９８には、第２の解析結果として、曲率の平均値と分散値とが統計データベース５５のどこに位置するかを表すグラフ７３０が示されている。グラフ７３０内の丸形状の点７２２は、対象としている被検眼の解析値をこのグラフの中にプロットしたものである。

なお、図２１に示す画面例では、マップやグラフをそれぞれ１つずつしか表示していないが、これに限らない。例えば、面積・体積マップと曲率マップとを第１の解析結果表示部９６内に並べて表示しても良いし、また、それと同時に、面積・体積グラフと曲率グラフとを第２の解析結果表示部９８内に並べて表示しても良い。

以上説明したように実施形態３によれば、網膜層の形状特徴を示す特徴量（形状特徴量：第１の特徴量）と、統計データベースに予め保持された複数の被検眼に対応する形状特徴量（第２の特徴量）とを操作者に表示（提示）することができる。これにより、被検眼が疾病であるか否かの判定基準となる指標とともに、診断対象となる被検眼の形状特徴量を提示することができる。

（実施形態４）
次に、実施形態４について説明する。実施形態４では、網膜層の形状特徴（形状特徴量）と断層画像又は眼底画像の画像特徴（画像特徴量）とを用いて、統計データベースに保持されている特徴量を検索する。そして、現在の被検眼の形状特徴に近い特徴量を持つ被検眼の断層画像、眼底画像、診断情報を表示する場合について説明する。

図２２は、実施形態４に係わる画像処理システム１０の構成の一例を示す図である。ここでは、実施形態３を説明した図１７の構成と相違する点について重点的に説明する。実施形態４に係わる画像処理装置３０には、実施形態３の構成に加えて、検索部３６が新たに設けられる。

また、実施形態４に係わる外部記憶装置５２には、統計データベース５５として、網膜層の特徴量と、断層画像と、眼底画像と、診断情報とがそれぞれ関連付けて記憶されている。検索部３６は、網膜層の形状特徴と、断層画像又は眼底画像内の画像特徴とを用いて、外部記憶装置５２（すなわち、統計データベース５５）を検索する。

次に、実施形態４に係わる画像処理装置３０の処理の流れの一例について説明する。実施形態４に係わる画像処理装置３０の処理は、上述した実施形態３と同様の流れで行なわれるため、ここでは、上述した図１８を用いて、相違する処理について重点的に説明する。相違する処理としては、Ｓ３０５、Ｓ３０６、Ｓ３０８の処理が挙げられる。

［Ｓ３０６］
画像処理装置３０は、特徴量取得部６０１において、上述した実施形態３で説明した形状特徴量の他、断層画像又は眼底画像から画像特徴量を抽出する。画像特徴量としては、例えば、エッジ特徴、色特徴、ヒストグラム特徴、ＳＩＦＴ（Scale-invariant feature transform）特徴量などが挙げられる。

［Ｓ３０７］
画像処理装置３０は、検索部３６において、Ｓ３０５の処理で求められた形状特徴量と画像特徴量とを多次元ベクトルとして表現し、当該特徴量と、外部記憶装置５２内の統計データベース５５に保持されたデータとを比較する。それにより、特徴量の類似する画像を検索する。探索方法としては、最近傍探索を用いれば良い。なお、処理の高速化のため、近似最近傍探索などを用いても良い。

［Ｓ３０８］
画像処理装置３０は、表示制御部３５において、撮影された断層画像と、検索部３６の検索により得られた類似する断層画像とを含む画面を表示装置５３に表示する。

具体的には、図２３に示す断層画像観察画面が表示される。断層画像観察画面９００には、現在の画像解析結果を表示する第１の解析結果表示部９３０と、過去の画像解析結果を表示する第２の解析結果表示部９４０とが設けられる。

第１の解析結果表示部９３０には、断層画像９０１、ＳＬＯ画像９０２、眼底写真９０３、（３次元の形状）解析結果マップ９０４、及び解析グラフ９０５が表示される。解析グラフ９０５内の丸形状の点９０６は、対象としている被検眼の解析値をグラフの中にプロットした結果を示している。また、星形状の点９０７は、第２の解析結果表示部９４０に表示されている過去の被検眼の解析値をグラフの中にプロットした結果を示している。この場合、解析グラフ９０５では、横軸が形状特徴量を示し、縦軸が画像特徴量を示す。勿論、これに限られず、解析結果として、実施形態３で説明した曲率値を示すグラフなどで表示しても良い。

第２の解析結果表示部９４０には、断層画像９１１及び９２１、ＳＬＯ画像９１２及び９２２、眼底写真９１３及び９２３、（３次元の形状）解析結果マップ９１４及び９２４、被検眼の診断情報９１５及び９２５が表示される。診断情報は、診断経過を示す情報であり、例えば、年齢、診断名、治療歴、視力、レフ値、ケラト値、眼軸長などが含まれる。図２３の場合、縦の列に沿って同一の被検眼の画像及び診断情報が表示されているが、勿論、横に沿って、同一の被検眼の画像及び診断情報を並べて表示するようにしても良い。

ここで、現在の被検眼と類似する画像が複数存在する場合、操作者は、入力装置５４を介してスライダーバー９５０を操作することにより、それらの結果を切り替えて表示することができる。なお、類似する過去の画像の表示の切り替えは、スライダーバー９５０以外を用いて行なわれても良い。例えば、操作者が、入力装置（例えば、マウス）５４を介して、解析グラフ９０５の星形状の点９０７をクリックした場合に、その特徴量を持つ画像を、現在の画像の横に並べて表示するようにしても良い。

以上説明したように実施形態４によれば、網膜層の形状特徴を示す特徴量（形状特徴量）と、断層画像又は眼底画像の画像特徴量とに基づいて、撮影された被検眼の断層画像に類似する断層画像やそれに関する情報を統計データベースから検索する。そして、撮影された断層画像とともに、当該類似する断層画像やそれに関する情報を並べて表示する。これにより、被検眼が疾病であるか否かの判定基準となる指標とともに、診断対象となる被検眼の形状特徴量を提示することができる。

以上が本発明の代表的な実施形態の例であるが、本発明は、上記及び図面に示す実施形態に限定することなく、その要旨を変更しない範囲内で適宜変形して実施できるものである。

例えば、上述した実施形態１〜４において説明したいくつかの画面においては、チェックボックスやコンボボックス等が配されていたが、ラジオボタンやリストボックス、ボタンなどは、用途に応じて適宜変更すれば良い。例えば、上述した説明でコンボボックスとして説明していた構成を、リストボックスとして実現しても良い。

（その他の実施形態）
本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（又はＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。

１０：画像処理システム、２０：断層画像撮影装置、３０：画像処理装置、５１：眼底画像撮影装置

Claims (12)

1. 眼底の断層画像を得るための断層画像撮影装置により取得された前記断層画像を処理して、前記眼底の網膜層の形状情報を取得する画像処理システムであって、
   眼軸長で分類された複数の眼に関する網膜層の形状情報が記憶されたデータベースを参照して、前記網膜層の形状情報に関する解析情報を取得する解析処理手段を備えることを特徴する画像処理システム。
2. 解析対象の網膜層を選択する選択手段を更に備え、
   前記解析処理手段は、前記選択手段により選択された網膜層の形状情報に関する解析情報を取得する
   ことを特徴とする請求項１に記載の画像処理システム。
3. 断層画像と、前記選択手段で選択された網膜層を識別可能に表示する表示制御手段を更に備えることを特徴とする請求項２に記載の画像処理システム。
4. 前記形状情報は、網膜層の曲率又は網膜層の層厚を含み、
   前記表示制御手段は、解析結果である曲率マップ又は層厚マップを表示することを特徴とする請求項３に記載の画像処理システム。
5. 近視眼解析用の撮影モードを含む複数の撮影モードを有し、
   前記解析処理手段は、断層画像の撮影時の撮影モードに応じて解析処理を切り替えて実行することを特徴とする請求項１または２に記載の画像処理システム。
6. 眼底の断層画像を取得する取得手段と、
   前記断層画像を処理して、前記眼底の網膜層の形状情報を取得する解析手段と、を備え、
   前記網膜層の形状情報が、上に凸であるか下に凸であるかを示す符号付きの曲率値であることを特徴する画像処理システム。
7. 前記取得された曲率値に基づき曲率マップを生成する生成手段と、
   前記曲率マップを表示する表示制御手段と、を更に備えることを特徴とする請求項６に記載の画像処理システム。
8. 前記眼底の眼底画像を取得する眼底画像取得手段と、
   前記眼底画像と前記曲率マップとを位置合わせする位置合わせ手段とを更に備え、
   前記表示制御手段は、前記位置合わせの結果に基づいて前記眼底画像に重畳された前記曲率マップを表示することを特徴とする請求項７に記載の画像処理システム。
9. 複数の被検眼から得られたデータが保持されたデータベースと、
   前記解析手段により取得された曲率値と、前記データベースに保持されたデータとを比較する比較手段とを更に備えることを特徴とする請求項６乃至８の何れか１項に記載の画像処理システム。
10. 画像処理システムの処理方法であって、
    取得手段が、眼底の断層画像を得るための断層画像撮影装置により取得された前記断層画像を処理して、前記眼底の網膜層の形状情報を取得する工程と、
    解析手段が、眼軸長で分類された複数の眼に関する網膜層の形状情報が記憶されたデータベースを参照して、前記網膜層の形状情報に関する解析情報を取得する工程と、を有することを特徴する画像処理システムの処理方法。
11. 画像処理システムの処理方法であって、
    取得手段が、眼底の断層画像を取得する工程と、
    解析手段が、前記断層画像を処理して、前記眼底の網膜層の形状情報を取得する工程と、を有し、
    前記網膜層の形状情報が、上に凸であるか下に凸であるかを示す符号付きの曲率値であることを特徴する画像処理システムの処理方法。
12. 請求項１から９の何れか１項に記載の画像処理システムによる処理をコンピュータに実行させるためのプログラム。