JP5947497B2

JP5947497B2 - 眼科画像処理装置、眼科画像処理方法、及び該方法を実行させるためのプログラム

Info

Publication number: JP5947497B2
Application number: JP2011152860A
Authority: JP
Inventors: 康文高間
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 2011-07-11
Filing date: 2011-07-11
Publication date: 2016-07-06
Anticipated expiration: 2031-07-11
Also published as: JP2013017619A

Description

本発明は、眼科診療などに用いられる眼科画像処理装置の画像生成技術に関し、より詳細には、眼科画像処理装置、眼科画像処理方法、及び該方法をコンピュータにより実行させるためのコンピュータプログラムに関する。

光干渉断層計（ＯｐｔｉｃａｌＣｏｈｅｒｅｎｃｅＴｏｍｏｇｏｒａｐｈｙ、以下、ＯＣＴとする）の登場により、網膜の複数の２次元断層画像（以下、断層像とする）から成るボリューム画像を得ることが可能となった。眼科の医療現場においてユーザは、ボリューム画像から層構造や、病変部の状態や分布、眼底血管の走行、網膜層厚の経時変化などを観察することができる。

ＯＣＴに限らず、ボリューム画像から病変の状態や分布、血管の走行などを観察するためには、視線方向に垂直な投影面に、ボクセルの輝度値を投影した投影像が有効である。特許文献１では、時系列心臓ボリューム画像の表面形状の経時変化を観察するために、異なる時相の心臓の関心領域の境界面から投影面までの距離を視線方向に沿って計算し、その差分に応じて色を付けて表示している。特許文献２では、対象組織の異常運動を診断するために、対象組織の基準点（重心）と対象組織表面上の各部位との距離を計算し、その変化を変位量として表示する方法が開示されている。

特開２００７−２８９５６９医用画像処理装置及び医用画像処理方法特開２００４−１９５０８２超音波診断装置

網膜の白斑や血管などの高反射な領域は、ＯＣＴボリューム画像において高輝度領域として観察される。これら高輝度な組織を観察するためには、３次元空間に置かれた投影面に視線上のボクセルの中から、最大輝度値を持つボクセルを投影する最大値投影法が有効である。

しかしながら、硝子体側に位置する神経線維層や脈絡膜側の網膜色素上皮の輝度値は血管や白斑と同程度に高輝度である。そのため、最大値投影法を用いると、白斑や血管は神経線維層や網膜色素上皮と同じような輝度値で表示され、観察することが困難である。

本発明は上記課題に鑑みてなされたものであり、ボリューム空間中に設定した基点から所定の値を持つボクセルまでの距離を用いて、白斑の分布や血管の走行が観察できる画像を生成することを目的とする。

上記の目的を達成するための、本発明の一態様による眼科画像処理装置は、以下の構成を備える。

すなわち、複数のボクセルからなる眼部のボリューム画像を取得するボリューム画像取得手段と、前記複数のボクセル内から閾値以上の輝度値を持つボクセルを所定の方向に探索することにより検出するボクセル検出手段と、前記ボリューム画像において前記眼部の層境界のうち少なくとも内境界膜と視細胞内節外節接合部との境界を検出する層境界検出手段と、前記視細胞内節外節接合部の連続性に応じて、前記内境界膜と前記視細胞内節外節接合部とのいずれかの各ボクセルに基点を設定する基点設定手段と、前記所定の方向に対して交差する平面に対して、前記設定された基点から前記検出されたボクセルまでの前記所定の方向における距離に応じて変化する表示形態を反映して得た平面画像を生成する画像生成手段と、を備えることを特徴とする。

本発明によれば、基点から所定の値を持つボクセルまでの距離に応じた色を画素値とする画像を生成する。それにより、ユーザが網膜の血管走行や白斑分布を容易に観察できる画像を生成することが可能になる。さらに、色情報から観察対象が網膜の奥行き方向に分布する状況を観察することも可能になる。

実施例１に係る眼科画像処理装置１の機能構成を示す図実施例に係る眼科画像処理装置１と接続される機器の構成を示す図実施例１に係る眼科画像処理装置１の処理手順を示すフローチャート実施例１に係る基点と検出したボクセルとの距離を示す図実施例１に係る血管走行と白斑分布を示す画像の模式図実施例２に係る眼科画像処理装置１の機能構成を示す図実施例２に係る眼科画像処理装置１の処理手順を示すフローチャート網膜の層構造を示す模式図実施例２に係るＩＳ／ＯＳの連続性が失われた場合の基点と検出したボクセルとの距離を示す図実施例２に係るＩＳ／ＯＳの連続性を保持する場合の基点と検出したボクセルとの距離を示す図

［実施例１］
本実施例では、断層像の上端に設定した基点から最大輝度値を持つボクセルまでの距離に応じた表示形態の画像を生成する方法について説明する。この画像により、ユーザは血管走行や白斑分布のような高輝度な値を持つ領域を観察することが可能になる。以下に、本実施例の一態様を説明する。

図２は、本実施例に係る眼科画像処理装置１に接続される機器の構成図である。図２に示すように、眼科画像処理装置１は、断層像撮像装置２及びデータサーバ３と、イーサネット（登録商標）等によるローカル・エリア・ネットワーク（ＬＡＮ）４を介して接続されている。ここで断層像撮像装置２は、眼部の断層像を撮像する装置であり、例えばＯＣＴが挙げられる。ＯＣＴは、１回の撮像で複数の断層像を取得するため、これらの断層像を順番に並べることによって、網膜のボリューム画像を取得することができる。断層像撮像装置２は、ユーザ（技師や医師）による操作に応じて被検者（患者）の断層像を撮像し、取得したボリューム画像を眼科画像処理装置１へ出力する。また、眼科画像処理装置１は、断層像撮像装置２により得られたボリューム画像を格納するデータサーバ３に接続され、そこから必要なボリューム画像を取得するように構成されてもよい。なお、これらの機器との接続は、ＵＳＢやＩＥＥＥ１３９４等のインターフェイスを介して行ってもよい。また、ＬＡＮ４によって、インターネット等の外部ネットワークを介して接続される構成であってもよい。

図１に本実施例における眼科画像処理装置１の構成を示す。眼科画像処理装置１は、ボリューム画像取得部１１、基点設定部１２、ボクセル検出部１３、距離計算部１４、画像生成部１５から構成される。

次に、図３のフローチャートを利用して、本実施例における眼科画像処理装置１の処理手順を説明する。

ステップＳ３０１において、ボリューム画像取得部１１は、断層像撮像装置２により撮像されたボリューム画像や、データサーバ３に蓄えられているボリューム画像を取得し、ボクセル検出部１３に出力する。該ボリューム画像は複数のボクセルによって構成されている。

ステップＳ３０２において、基点設定部１２は、距離の基点を設定し、距離計算部１４へ出力する。本実施例では、基点４０１は図４に示す断層像４の上端、つまりＺ＝０におけるＸＹ平面上に設定する。なお、後述する他の実施形態にも示すように、起点の設定位置としての該ＸＹ平面は本実施形態の場合であり、当該基点は取得されたボリューム画像中の任意の条件に従った点として設定可能である。

ステップＳ３０３において、ボクセル検出部１３は、図４に示す断層像４のＺ軸に沿って基点から網膜の奥行き方向にボクセルを走査しながら所定の値を持つボクセル４０２を特定のボクセルとして探索し、検出する。ここで、所定の値を持つボクセルとは、例えば、奥行き方向に走査した複数のボクセルの内、最大輝度値を持つボクセルとする。これは、ＯＣＴで撮像した断層像において、観察対象の血管や白斑が高輝度の組織として描出されるからである。一方、図４に示すＩＳ／ＯＳ（視細胞内節外節接合部）４０５とＲＰＥ（網膜色素上皮）４０６間の組織の輝度値はさらに高い。しかしながら、血管や白斑が存在する領域は光がＩＳ／ＯＳ４０５まで届かないために、ＩＳ／ＯＳ４０５とＲＰＥ４０６間の組織は輝度値が低く描出される。よって、ボクセル検出部１３は、血管や白斑が存在しない領域ではＩＳ／ＯＳとＲＰＥ間のボクセルを、それらが存在する領域では血管や白斑のボクセルを最大輝度値ボクセルとして検出する。検出されたボクセルの座標は距離計算部１４へ出力される。

ステップＳ３０４において、距離計算部１４は、ステップＳ３０２において設定した基点４０１の座標と、ステップＳ３０３において検出した特定のボクセル４０２の座標との差を計算し、図４に示す断層像４における距離４０３を求める。求めた距離は画像生成部１５へ出力される。

ステップＳ３０５において、ステップＳ３０４において求めた距離に応じた表示形態の画素値を得る。画像生成部１５は、距離に応じて実際の画像表示の際に用いる表示形態を設定する。例えば、表示形態に色を用いた場合、距離に応じた色を設定するために色テーブルを用いる。距離に応じて色が線形に変化するように色テーブルを設定すれば、画像の色情報から最大輝度値ボクセルが存在している奥行き情報を読みとることが可能となる。また、表示形態に濃淡を用いても良い。この場合、距離が長いほど濃く短いほど淡く色を設定すれば、画像の濃淡情報から奥行き情報を読み取ることが可能となる。

ステップＳ３０６において、最初の基点を設定したＸＹ平面上の全てのボクセルに対して基点を設定し終えたか否かを確認する。設定し終えていない場合にステップＳ３０２に戻る。ステップＳ３０２からステップＳ３０５までの処理を繰り返し、Ｚ＝０のＸＹ平面上のすべてのボクセルを基点に設定しながら、Ｚ軸に沿って最大輝度値ボクセルを検出し、２つのボクセルの距離に応じた色を画素値に設定する。従って、該最大輝度値ボクセルは、基点の数に応じて複数求められる。また、前述して設定された表示形態に応じた色調、濃淡等を設定する。

ＸＹ平面上の全てのボクセルに対しての基点設定と該ボクセルに基づいた表示形態の設定が終了したと認められるとステップＳ３０７に移行する。ステップＳ３０７では、設定された全ての表示形態を各々整合させ、新たなボリューム画像を生成する。このようにして、図４の断層像４のＸＹ平面に平行な面にその距離情報を投影した画像を生成する。

以上で述べた構成によれば、設定した基点から最大輝度値ボクセルまでの距離に応じた色を画素値とする画像を生成する。血管や白斑が存在する領域と存在しない領域において、最大輝度値ボクセルが検出されるＺ座標が異なる。したがって、図５に示すように、距離に応じて諧調値が異なる画像が生成され、白斑５０１や血管５０２の位置を識別することが可能である。

（変形例１−１）
上記の実施例では、ステップＳ３０３において、所定の値を最大輝度値とした。本変形例では、設定した閾値以上の値とした場合について説明する。

ステップＳ３０３において、ボクセル検出部１３は、ボクセルの輝度に関しての閾値を設定し、閾値以上の輝度値を持つボクセルを１つ以上検出し、例えば、それらのＺ座標の平均値を計算する。計算した値は距離計算部１４へ出力される。従って、この場合検出される特定のボクセルは複数となり、距離計算は基点とこれら特性のボクセル各々の座標との間で複数回行われる。

ステップＳ３０４において、距離計算部１４は、ステップＳ３０３において計算した値と基点の座標との差を距離として計算し、画像生成部１５へ出力する。

このように実施することによって、ＯＣＴの断層像に含まれる高輝度なノイズの影響を抑えたり、厚みを持った高輝度領域のＺ座標の平均値までの距離を示したりすることが可能である。

［実施例２］
実施例１では、基点をＺ＝０のＸＹ平面上のボクセルとした場合について、血管走行や白斑分布を観察できる画像を生成する方法について述べた。本実施例では、基点を網膜の層境界上のボクセルとした場合について、その方法を述べる。

図６に本実施例における眼科画像処理装置１の構成を示す。眼科画像処理装置１は、ボリューム画像取得部６１、基点設定部６２、ボクセル検出部６３、距離計算部６４、画像生成部６５、層境界検出部６６から構成される。

次に、図７のフローチャートを利用して、本実施例における眼科画像処理装置１の処理手順を説明する。

ステップＳ７０１において、ボリューム画像取得部６１は、断層像撮像装置２により撮像されたボリューム画像や、データサーバ３に蓄えられているボリューム画像を取得し、ボクセル検出部６３と層境界検出部６６に出力する。

ステップＳ７０２において、層境界検出部６６は、ボリューム画像から、網膜層内部の各層を抽出する。網膜層内部の層抽出について図８を用いて説明する。図８は網膜層内の各層８０１〜８１０と、ｚ軸方向に並行な１列の画素列Ａ（Ａ−ｓｃａｎ）を示している。８０１はＩＬＭ（内境界膜）、８０２はＮＦＬ（神経線維層）、８０３はＧＣＬ（神経節細胞層）、８０４はＩＰＬ（内網状層）、８０５はＩＮＬ（内顆粒層）、８０６はＯＰＬ（外網状層）、８０７はＯＮＬ（外顆粒層）、８０８はＥＬＭ（外境界膜）、８０９はＩＳ/ＯＳ（視細胞内節外節接合部）、８１０はＲＰＥ（網膜色素上皮層）を表す。また、層と層の境界は、上側層／下側層、のように表す。本実施例において、層境界検出部６６はＩＬＭ８０１とＮＦＬ８０２／ＧＣＬ８０３、ＩＰＬ８０４／ＩＮＬ８０５、ＯＰＬ８０６／ＯＮＬ８０７、ＩＳ／ＯＳ８０９、ＲＰＥ８１０を少なくとも抽出する。なお、層境界検出部６６が層を抽出するということと、層と層の境界を抽出することは同じである。例えば、硝子体と網膜との境界のＩＬＭ８０１と、ＮＦＬ８０２／ＧＣＬ８０３との境界を抽出することで、ＮＦＬ８０２を抽出することになる。

網膜層の抽出について、まず、断層像に対して、メディアンフィルタとＳｏｂｅｌフィルタをそれぞれ適用して画像を作成する（以下、メディアン画像、Ｓｏｂｅｌ画像とする）。次に、変換したメディアン画像とＳｏｂｅｌ画像から、Ａ−ｓｃａｎ毎にプロファイルを作成する。メディアン画像では輝度値のプロファイル、Ｓｏｂｅｌ画像では勾配のプロファイルとなる。そして、Ｓｏｂｅｌ画像から作成したプロファイル内のピークを検出する。検出したピークの前後やピーク間に対応するメディアン画像のプロファイルを参照することで、網膜層を検出する。検出された層境界情報を基点設定部６２とボクセル検出部６３へ出力する。

ステップＳ７０３において、基点設定部６２はＩＳ／ＯＳの連続性を確認する。図９に示すように、ＯＣＴの断層像９において、高輝度領域９０７の存在により、それよりも奥のボクセルの輝度値は暗くなる。また、視神経乳頭部ではＩＳ／ＯＳが描出されない。よって、高輝度領域が存在するボリューム画像や視神経乳頭部のボリューム画像においては、図９に示すように、ＩＳ／ＯＳ９０５が分断される。そこで、基点設定部６２は入力された層境界の情報を利用して、ＩＳ／ＯＳの連続性を確認する。例えば、層境界が３次元座標の点列で表現されている場合、隣接する点が一定間隔以上離れている場合、連続性が失われていると判断する。

ＩＳ／ＯＳの連続性が失われている場合、ステップＳ７０４において、多くのＯＣＴの断層像に存在するＩＬＭを基点として設定する。例えば、図９に示すように、ＩＬＭ９０４上の基点９０１を利用する。ＩＬＭをＸＹ平面に射影し、ＩＬＭを平坦にすることにより、ＩＬＭの凹凸の影響を軽減した画像を生成することが可能である。

ステップＳ７０６において、ボクセル検出部６３は、図９に示す断層像９のＺ軸に沿って奥行き方向にボクセルを走査しながら最大輝度値を持つボクセル９０２を検出する。

ステップＳ７０７において、距離計算部６４は、ステップＳ７０４において設定した基点の座標と、ステップＳ７０６において検出したボクセルの座標の差を計算し、図９に示す断層像９における距離９０３を求める。求めた距離は画像生成部６５へ出力される。

ステップＳ７０８において、ステップＳ７０７において求めた距離９０３に応じた色を画素値とする。また、求めた距離９０３に応じて表示形態を設定する。ＯＣＴの断層像において血管はＩＬＭよりに、白斑は網膜の内部に多く分布する。このことから、色テーブルは、例えば、距離が短いほど赤、距離が長いほど黒、その中間を白に設計する。このように設計することで、ＩＬＭよりに存在する血管を赤く、網膜の内部に存在する白斑を白く、観察対象ではないＩＳ／ＯＳを黒く表現することが可能である。

ステップＳ７０９において、基点の設定を行った層境界上における基点設定の可能な領域の全域のボクセルに対して基点を設定し終えたか否かを確認する。設定し終えていない場合にステップＳ７０３に戻る。

ステップＳ７０３からステップＳ７０９までの処理を繰り返し、ＩＬＭ上のすべてのボクセルを基点に設定しながら、Ｚ軸に沿って最大輝度値ボクセルを検出し、２つのボクセルの距離に応じた色を画素値に設定する。

層境界上の全てのボクセルに対しての基点設定と該ボクセルに基づいた表示形態の設定が終了したと認められるとステップＳ７１０に移行する。ステップＳ７１０では、設定された全ての表示形態を各々整合させ、新たなボリューム画像を生成する。このようにして、図９の断層像９のＸＹ平面に平行な面にその距離情報を投影した画像を生成する。

ＩＳ／ＯＳの連続性が保持されている場合、ステップＳ７０５において、ＩＳ／ＯＳはＩＬＭよりも起伏が少ないことから、ＩＳ／ＯＳ上のボクセルを基点として設定する。例えば、図１０に示すように、ＩＳ／ＯＳ１００５上の基点１００１を利用する。網膜全体が湾曲している場合、ＩＳ／ＯＳをＸＹ平面に射影し、ＩＳ／ＯＳを平坦にすることにより、網膜の湾曲の影響を軽減した画像を生成することが可能である。

ステップＳ７０６において、ボクセル検出部６３は、図１０に示す断層像１０のＺ軸に沿って、ＩＬＭ１００４からＩＳ／ＯＳ１００５までの間に存在するボクセルを走査しながら最大輝度値を持つボクセル１００２を検出する。

ステップＳ７０７において、距離計算部６４は、ステップＳ７０５において設定した基点の座標と、ステップＳ７０６において検出したボクセルの座標の差を計算し、図１０に示す断層像１０における距離１００３を求める。求めた距離は画像生成部６５へ出力される。

ステップＳ７０８において、ステップＳ７０７において求めた距離に応じた色を画素値とする。また、求めた距離９０３に応じて表示形態を設定する。距離に応じた色を設定するための色テーブルは、例えば、距離が長いほど赤、距離が短いほど黒、その中間を白に設計する。このように設計することにより、血管を赤く、白斑を白く、ＩＳ／ＯＳを黒く表現することが可能である。

ステップＳ７０３からステップＳ７０９までの処理を繰り返し、ＩＳ／ＯＳ上のすべてのボクセルを基点に設定しながら、Ｚ軸に沿って最大輝度値ボクセルを検出し、２つのボクセルの距離に応じた色を設定する。

層境界上の全てのボクセルに対しての基点設定と該ボクセルに基づいた表示形態の設定が終了したと認められるとステップＳ７１０に移行する。ステップＳ７１０では、設定された全ての表示形態を各々整合させ、新たなボリューム画像を生成する。このようにして、図１０の断層像１０のＸＹ平面に平行な面にその距離情報を投影した画像を生成する。

層境界情報を用いることにより、層境界に色を設定した色テーブルを自動で設計することが可能になる。また、網膜外の輝度値に影響されることもない。そのため、網膜内における血管走行や白斑分布を示す、より見易い画像を生成することが可能となる。

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

１眼科画像処理装置
１１ボリューム画像取得部
１２基点設定部
１３ボクセル検出部
１４距離計算部
１５画像生成部

Claims

複数のボクセルからなる眼部のボリューム画像を取得するボリューム画像取得手段と、
前記複数のボクセル内から閾値以上の輝度値を持つボクセルを所定の方向に探索することにより検出するボクセル検出手段と、
前記ボリューム画像において前記眼部の層境界のうち少なくとも内境界膜と視細胞内節外節接合部との境界を検出する層境界検出手段と、
前記視細胞内節外節接合部の連続性に応じて、前記内境界膜と前記視細胞内節外節接合部とのいずれかの各ボクセルに基点を設定する基点設定手段と、
前記所定の方向に対して交差する平面に対して、前記設定された基点から前記検出されたボクセルまでの前記所定の方向における距離に応じて変化する表示形態を反映して得た平面画像を生成する画像生成手段と、
を備えることを特徴とする眼科画像処理装置。
前記ボクセル検出手段は、前記基点から前記眼部の奥行き方向に前記ボクセルを探索して前記閾値以上の輝度値を持つボクセルを検出することを特徴とする請求項１に記載の眼科画像処理装置。
前記閾値を設定する閾値設定手段を更に備え、
前記ボクセル検出手段は、前記設定された閾値以上の輝度値を持つボクセルを検出することを特徴とする請求項１又は２に記載の眼科画像処理装置。
前記設定された基点から前記検出されたボクセルまでの距離を計算する距離計算手段を更に備え、
前記画像生成手段は、前記平面に対して、前記計算された距離に応じて変化する表示形態を反映して得た平面画像を生成することを特徴とする請求項１乃至３の何れか一項に記載の眼科画像処理装置。
前記検出されたボクセルは複数あり、
前記距離計算手段は前記複数のボクセル各々の座標に基づいて得られる特定の座標と前記設定された基点とから前記距離を計算することを特徴とする請求項４に記載の眼科画像処理装置。
前記基点設定手段は、前記視細胞内節外節接合部の連続性が失われていないと判断された場合には前記視細胞内節外節接合部の各ボクセルに前記基点を設定し、前記視細胞内節外節接合部の連続性が失われていると判断された場合には前記内境界膜の各ボクセルに前記基点を設定することを特徴とする請求項１乃至５の何れか一項に記載の眼科画像処理装置。
前記表示形態は、色および濃淡の少なくとも何れかであることを特徴とする請求項１乃至６の何れか一項に記載の眼科画像処理装置。
複数のボクセルからなる眼部のボリューム画像を取得するボリューム画像取得工程と、
前記複数のボクセル内から閾値以上の輝度値を持つボクセルを所定の方向に探索することにより検出するボクセル検出工程と、
前記ボリューム画像において前記眼部の層境界のうち少なくとも内境界膜と視細胞内節外節接合部との境界を検出する層境界検出工程と、
前記視細胞内節外節接合部の連続性に応じて、前記内境界膜と前記視細胞内節外節接合部とのいずれかの各ボクセルに基点を設定する基点設定工程と、
前記所定の方向に対して交差する平面に対して、前記設定された基点から前記検出されたボクセルまでの前記所定の方向における距離に応じて変化する表示形態を反映して得た平面画像を生成する画像生成工程と、
を有することを特徴とする眼科画像処理方法。
前記ボクセル検出工程では、前記基点から眼部の奥行き方向に前記ボクセルを探索して前記閾値以上の輝度値を持つボクセルを検出することを特徴とする請求項８に記載の眼科画像処理方法。
前記閾値を設定する閾値設定工程を更に有し、
前記ボクセル検出工程では、前記設定された閾値以上の輝度値を持つボクセルを検出することを特徴とする請求項８又は９に記載の眼科画像処理方法。
前記設定された基点から前記検出されたボクセルまでの距離を計算する距離計算工程を更に有し、
前記画像生成工程では、前記平面に対して、前記計算された距離に応じて変化する表示形態を反映して得た平面画像を生成することを特徴とする請求項８乃至１０の何れか一項に記載の眼科画像処理方法。
前記検出されたボクセルは複数あり、
前記距離計算工程では、前記複数のボクセル各々の座標に基づいて得られる特定の座標と前記設定された基点とから前記距離を計算することを特徴とする請求項１１に記載の眼科画像処理方法。
前記基点設定工程において、前記視細胞内節外節接合部の連続性が失われていないと判断された場合には前記視細胞内節外節接合部の各ボクセルに前記基点を設定し、前記視細胞内節外節接合部の連続性が失われていると判断された場合には前記内境界膜の各ボクセルに前記基点を設定することを特徴とする請求項８乃至１２の何れか一項に記載の眼科画像処理方法。
前記表示形態は、色および濃淡の少なくとも何れかであることを特徴とする請求項８乃至１３の何れか一項に記載の眼科画像処理方法。
請求項８乃至１４の何れか一項に記載の眼科画像処理方法をコンピュータに実行させるためのコンピュータプログラム。