特許 7416083 画像処理方法、画像処理装置、およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-01-09

(45)【発行日】2024-01-17

(54)【発明の名称】画像処理方法、画像処理装置、およびプログラム

(51)【国際特許分類】

   A61B 3/12 20060101AFI20240110BHJP

   G06T 7/00 20170101ALI20240110BHJP

【ＦＩ】

A61B3/12

G06T7/00 612

【請求項の数】 10

(21)【出願番号】P 2021552011

(86)(22)【出願日】2019-10-15

(86)【国際出願番号】 JP2019040483

(87)【国際公開番号】W WO2021074962

(87)【国際公開日】2021-04-22

【審査請求日】2022-04-15

(73)【特許権者】

【識別番号】000004112

【氏名又は名称】株式会社ニコン

(74)【代理人】

【識別番号】110001519

【氏名又は名称】弁理士法人太陽国際特許事務所

(72)【発明者】

【氏名】廣川 真梨子

(72)【発明者】

【氏名】田邉 泰士

【審査官】山口 裕之

(56)【参考文献】

【文献】特開２００８－０２９７３２（ＪＰ，Ａ）

【文献】特開２０１４－０８３２８５（ＪＰ，Ａ）

【文献】国際公開第２０１９／１８１９８１（ＷＯ，Ａ１）

【文献】特開２０１８－１６１２３０（ＪＰ，Ａ）

【文献】OHNO-MATSUI Kyoko et al.，"The long-term follow-up of a highly myopic patient with a macular vortex vein"，Acta Ophthalmologica Scandinavica，1997年06月30日，Vol.75, No.3，pp.329-332

【文献】MORIYAMA Muka et al.，"Detection of posterior vortex veins in eyes with pathologic myopia by ultra-widefield indocyanine g，British Journal of Ophthalmology，英国，2017年09月18日，Volume 101, Issue 9，pp.1179-1184

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ３／１２

(57)【特許請求の範囲】

【請求項1】

プロセッサが、第１眼底画像から渦静脈の第１位置を特定することと、
前記プロセッサが、第２眼底画像から前記渦静脈の第２位置を特定することと、
前記プロセッサが、前記第１位置および前記第２位置を表示する渦静脈マップのデータを生成することと、
を含む、画像処理方法であって、
前記第１位置を特定することは、前記第１眼底画像において脈絡膜血管の走行方向を解析し、前記脈絡膜血管の集中する度合いに基づいて前記第１位置を特定することであり、
前記第２位置を特定することは、前記第２眼底画像において脈絡膜血管の走行方向を解析し、前記脈絡膜血管の集中する度合いに基づいて前記第２位置を特定することである、
画像処理方法。

【請求項2】

前記第１眼底画像および前記第２眼底画像は、眼底を異なるタイミングで撮影した画像である、請求項１に記載の画像処理方法。

【請求項3】

前記渦静脈の前記第１位置および前記第２位置は、
眼底の第１構造物と前記第１構造物と異なる第２構造物とを結ぶ第１の線と前記第１構造物と前記渦静脈とを結ぶ第２の線との第１のなす角度、および、前記渦静脈と前記第１構造物との第１の距離の第１の組と、
前記渦静脈の前記第１位置および前記第２位置は、前記第１の線と前記第２構造物と前記渦静脈とを結ぶ第３の線との第２のなす角度、および、前記渦静脈と前記第２構造物との第２の距離の第２の組と、
の少なくとも一方により、特定される、
請求項１また請求項２に記載の画像処理方法。

【請求項4】

前記渦静脈マップは、前記渦静脈の前記第１位置および前記第２位置を、前記第１の組および前記第２の組の少なくとも一方を用いた極座標系、２次元直交座標系、および３次元直交座標系の少なくとも１つの座標系で表示したマップである、
請求項３に記載の画像処理方法。

【請求項5】

前記極座標系で表示したマップは、レーダチャートであり、
前記２次元直交座標系で表示したマップは、２次元直交軸で規定される散布図であり、
前記３次元直交座標系で表示したマップは、３次元直交軸で規定される散布図である、
請求項４に記載の画像処理方法。

【請求項6】

前記第１眼底画像および前記第２眼底画像の各々には、複数の前記渦静脈が存在し、
前記プロセッサが、複数の前記渦静脈の中心位置を求め、
前記渦静脈マップに前記中心位置が表示される、
請求項１から請求項５の何れか１項に記載の画像処理方法。

【請求項7】

前記第１眼底画像および前記第２眼底画像の各々には、前記渦静脈が複数存在し、
前記渦静脈マップのデータは、前記複数の渦静脈から選択された渦静脈を対象として生成される、
請求項１から請求項６の何れか１項に記載の画像処理方法。

【請求項8】

前記プロセッサが、前記渦静脈マップと、前記第１眼底画像および前記第２眼底画像と、を表示する表示ビューを生成することを更に含む請求項１から請求項７の何れか１項に記載の画像処理方法。

【請求項9】

メモリと、前記メモリに接続されたプロセッサとを備え、
前記プロセッサは、
第１眼底画像から渦静脈の第１位置を特定し、
第２眼底画像から前記渦静脈の第２位置を特定し、
前記第１位置および前記第２位置を表示する渦静脈マップのデータを生成する、
画像処理装置であって、
前記第１位置を特定することは、前記第１眼底画像において脈絡膜血管の走行方向を解析し、前記脈絡膜血管の集中する度合いに基づいて前記第１位置を特定することであり、
前記第２位置を特定することは、前記第２眼底画像において脈絡膜血管の走行方向を解析し、前記脈絡膜血管の集中する度合いに基づいて前記第２位置を特定することである、
画像処理装置。

【請求項10】

コンピュータに、
第１眼底画像から渦静脈の第１位置を特定し、
第２眼底画像から前記渦静脈の第２位置を特定し、
前記第１位置および前記第２位置を表示する渦静脈マップのデータを生成する、
ことを実行させるプログラムであって、
前記第１位置を特定することは、前記第１眼底画像において脈絡膜血管の走行方向を解析し、前記脈絡膜血管の集中する度合いに基づいて前記第１位置を特定し、
前記第２位置を特定することは、前記第２眼底画像において脈絡膜血管の走行方向を解析し、前記脈絡膜血管の集中する度合いに基づいて前記第２位置を特定することである、
プログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、画像処理方法、画像処理装置、およびプログラムに関する。

【背景技術】

【0002】

米国特許８６３６３６４号明細書には、眼底画像から渦静脈の位置を特定することが開示されている。

【0003】

渦静脈位置を経過観察することが望まれている。

【発明の概要】

【0004】

本開示の技術の第１の態様の画像処理方法は、プロセッサが、第１眼底画像から渦静脈の第１位置を特定することと、前記プロセッサが、第２眼底画像から前記渦静脈の第２位置を特定することと、前記プロセッサが、前記第１位置および前記第２位置を表示する渦静脈マップのデータを生成することと、を含む。

【0005】

本開示の技術の第２の態様の画像処理装置は、メモリと、前記メモリに接続されたプロセッサとを備え、前記プロセッサは、第１眼底画像から渦静脈の第１位置を特定し、第２眼底画像から前記渦静脈の第２位置を特定し、前記第１位置および前記第２位置を表示する渦静脈マップのデータを生成する。

【0006】

本開示の技術の第３の態様のプログラムは、コンピュータに、第１眼底画像から渦静脈の第１位置を特定し、第２眼底画像から前記渦静脈の第２位置を特定し、前記第１位置および前記第２位置を表示する渦静脈マップのデータを生成する、ことを実行させる。

【図面の簡単な説明】

【0007】

【図1】眼科システム１００のブロック図である。

【図2】眼科装置１１０の全体構成を示す概略構成図である。

【図3】サーバ１４０の電気系の構成のブロック図である。

【図4】サーバ１４０のＣＰＵ２６２の機能のブロック図である。

【図5】サーバ１４０による画像処理のフローチャートである。

【図6】図５のステップ５０６の渦静脈のフローチャートである。

【図7A】図５のステップ５０６の渦静脈の相対位置算出処理における渦静脈と視神経乳頭ＯＮＨとの距離ｒの計算処理のフローチャートである。

【図7B】図５のステップ５０６の渦静脈の相対位置算出処理における角度θの算出処理のフローチャートである。

【図8】ＲＧＢカラー眼底画像Ｇ１を示す図である。

【図9】脈絡膜血管画像Ｇ２を示す図である。

【図10】４個の渦静脈ＶＶ１からＶＶ４が検出された様子を示す図である。

【図11】黄斑Ｍの位置－視神経乳頭ＯＮＨの位置－ＶＶの位置を示す図である。

【図12】４象限の各々のｘの符号を示す図である。

【図13】仮想球面を示す図である。

【図14A】渦静脈の位置を、視神経乳頭ＯＮＨを基準とした角度θ及び距離ｒにより特定した図である。

【図14B】渦静脈の位置を、角度θ及び距離ｒを用いた極座標系で表示したマップ（レーダチャート）を示す図である。

【図14C】渦静脈の位置を、角度θ及び距離ｒを用いた２次元直交座標系で表示したマップを示す図である。

【図14D】渦静脈の位置を、角度θ及び距離ｒを用いた３次元直交座標系で表示したマップを示す図である。

【図15A】渦静脈の位置を、黄斑Ｍを基準とした角度φ及び距離ｓにより特定した図である。

【図15B】渦静脈の位置を、角度φ及び距離ｓを用いた極座標系で表示したマップ（レーダチャート）を示す図である。

【図15C】渦静脈の位置を、角度φ及び距離ｓを用いた２次元直交座標系で表示したマップを示す図である。

【図15D】渦静脈の位置を、角度φ及び距離ｓを用いた３次元直交座標系で表示したマップを示す図である。

【図16】第１眼底画像表示画面１０００Ａを示す図である。

【図17】第２眼底画像表示画面１０００Ｂを示す図である。

【図18】第３眼底画像表示画面１０００Ｃを示す図である。

【図19】各撮影日時の渦静脈の４つの位置の中心Ｃ１からＣ４、中心Ｃ１からＣ４の中心ＣＣを表示した図である。

【発明を実施するための形態】

【0008】

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

【0009】

図１を参照して、眼科システム１００の構成を説明する。図１に示すように、眼科システム１００は、眼科装置１１０と、眼軸長測定器１２０と、管理サーバ装置（以下、「サーバ」という）１４０と、画像表示装置（以下、「ビューワ」という）１５０と、を備えている。眼科装置１１０は、眼底画像を取得する。眼軸長測定器１２０は、患者の眼軸長を測定する。サーバ１４０は、眼科装置１１０によって患者の眼底が撮影されることにより得られた眼底画像を、患者のＩＤに対応して記憶する。ビューワ１５０は、サーバ１４０から取得した眼底画像などの医療情報を表示する。

【0010】

眼科装置１１０、眼軸長測定器１２０、サーバ１４０、およびビューワ１５０は、ネットワーク１３０を介して、相互に接続されている。

【0011】

次に、図２を参照して、眼科装置１１０の構成を説明する。

【0012】

説明の便宜上、走査型レーザ検眼鏡（Ｓｃａｎｎｉｎｇ Ｌａｓｅｒ Ｏｐｈｔｈａｌｍｏｓｃｏｐｅ）を「ＳＬＯ」と称する。また、光干渉断層計（Ｏｐｔｉｃａｌ Ｃｏｈｅｒｅｎｃｅ Ｔｏｍｏｇｒａｐｈｙ）を「ＯＣＴ」と称する。

【0013】

なお、眼科装置１１０が水平面に設置された場合の水平方向を「Ｘ方向」、水平面に対する垂直方向を「Ｙ方向」とし、被検眼１２の前眼部の瞳孔の中心と眼球の中心とを結ぶ方向を「Ｚ方向」とする。従って、Ｘ方向、Ｙ方向、およびＺ方向は互いに垂直である。

【0014】

眼科装置１１０は、撮影装置１４および制御装置１６を含む。撮影装置１４は、ＳＬＯユニット１８、ＯＣＴユニット２０、および撮影光学系１９を備えており、被検眼１２の眼底の眼底画像を取得する。以下、ＳＬＯユニット１８により取得された二次元眼底画像をＳＬＯ画像と称する。また、ＯＣＴユニット２０により取得されたＯＣＴデータに基づいて作成された網膜の断層画像や正面画像（ｅｎ－ｆａｃｅ画像）などをＯＣＴ画像と称する。

【0015】

制御装置１６は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ（中央処理装置））１６Ａ、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）１６Ｂ、ＲＯＭ（Ｒｅａｄ－Ｏｎｌｙ ｍｅｍｏｒｙ）１６Ｃ、および入出力（Ｉ／Ｏ）ポート１６Ｄを有するコンピュータを備えている。

【0016】

制御装置１６は、Ｉ／Ｏポート１６Ｄを介してＣＰＵ１６Ａに接続された入力／表示装置１６Ｅを備えている。入力／表示装置１６Ｅは、被検眼１２の画像を表示したり、ユーザから各種指示を受け付けたりするグラフィックユーザインターフェースを有する。グラフィックユーザインターフェースとしては、タッチパネル・ディスプレイが挙げられる。

【0017】

また、制御装置１６は、Ｉ／Ｏポート１６Ｄに接続された画像処理装置１６Ｇを備えている。画像処理装置１６Ｇは、撮影装置１４によって得られたデータに基づき被検眼１２の画像を生成する。制御装置１６はＩ／Ｏポート１６Ｄに接続された通信インターフェース（Ｉ／Ｆ）１６Ｆを備えている。眼科装置１１０は、通信インターフェース（Ｉ／Ｆ）１６Ｆおよびネットワーク１３０を介して眼軸長測定器１２０、サーバ１４０、およびビューワ１５０に接続される。

【0018】

上記のように、図２では、眼科装置１１０の制御装置１６が入力／表示装置１６Ｅを備えているが、本開示の技術はこれに限定されない。例えば、眼科装置１１０の制御装置１６は入力／表示装置１６Ｅを備えず、眼科装置１１０とは物理的に独立した別個の入力／表示装置を備えるようにしてもよい。この場合、当該表示装置は、制御装置１６のＣＰＵ１６Ａの制御下で動作する画像処理プロセッサユニットを備える。画像処理プロセッサユニットが、ＣＰＵ１６Ａが出力指示した画像信号に基づいて、ＳＬＯ画像等を表示するようにしてもよい。

【0019】

撮影装置１４は、制御装置１６のＣＰＵ１６Ａの制御下で作動する。撮影装置１４は、ＳＬＯユニット１８、撮影光学系１９、およびＯＣＴユニット２０を含む。撮影光学系１９は、第１光学スキャナ２２、第２光学スキャナ２４、および広角光学系３０を含む。

【0020】

第１光学スキャナ２２は、ＳＬＯユニット１８から射出された光をＸ方向、およびＹ方向に２次元走査する。第２光学スキャナ２４は、ＯＣＴユニット２０から射出された光をＸ方向、およびＹ方向に２次元走査する。第１光学スキャナ２２および第２光学スキャナ２４は、光束を偏向できる光学素子であればよく、例えば、ポリゴンミラーや、ガルバノミラー等を用いることができる。また、それらの組み合わせであってもよい。

【0021】

広角光学系３０は、共通光学系２８を有する対物光学系（図２では不図示）、およびＳＬＯユニット１８からの光とＯＣＴユニット２０からの光を合成する合成部２６を含む。

【0022】

なお、共通光学系２８の対物光学系は、楕円鏡などの凹面ミラーを用いた反射光学系や、広角レンズなどを用いた屈折光学系、あるいは、凹面ミラーやレンズを組み合わせた反射屈折光学系でもよい。楕円鏡や広角レンズなどを用いた広角光学系を用いることにより、視神経乳頭や黄斑が存在する眼底中心部だけでなく眼球の赤道部や渦静脈が存在する眼底周辺部の網膜を撮影することが可能となる。

【0023】

楕円鏡を含むシステムを用いる場合には、国際公開ＷＯ２０１６／１０３４８４あるいは国際公開ＷＯ２０１６／１０３４８９に記載された楕円鏡を用いたシステムを用いる構成でもよい。国際公開ＷＯ２０１６／１０３４８４の開示および国際公開ＷＯ２０１６／１０３４８９の開示の各々は、その全体が参照により本明細書に取り込まれる。

【0024】

広角光学系３０によって、眼底において広い視野（ＦＯＶ：Ｆｉｅｌｄ ｏｆ Ｖｉｅｗ）１２Ａでの観察が実現される。ＦＯＶ１２Ａは、撮影装置１４によって撮影可能な範囲を示している。ＦＯＶ１２Ａは、視野角として表現され得る。視野角は、本実施の形態において、内部照射角と外部照射角とで規定され得る。外部照射角とは、眼科装置１１０から被検眼１２へ照射される光束の照射角を、瞳孔２７を基準として規定した照射角である。また、内部照射角とは、眼底へ照射される光束の照射角を、眼球中心Ｏを基準として規定した照射角である。外部照射角と内部照射角とは、対応関係にある。例えば、外部照射角が１２０度の場合、内部照射角は約１６０度に相当する。本実施の形態では、内部照射角は２００度としている。

【0025】

内部照射角の２００度は、本開示の技術の「所定値」の一例である。

【0026】

ここで、内部照射角で１６０度以上の撮影画角で撮影されて得られたＳＬＯ眼底画像をＵＷＦ－ＳＬＯ眼底画像と称する。なお、ＵＷＦとは、ＵｌｔｒａＷｉｄｅ Ｆｉｅｌｄ（超広角）の略称を指す。もちろん、内部照射角で１６０度未満の撮影画角で眼底を撮影することにより、ＵＷＦではないＳＬＯ画像を取得することができる。

【0027】

ＳＬＯシステムは、図２に示す制御装置１６、ＳＬＯユニット１８、および撮影光学系１９によって実現される。ＳＬＯシステムは、広角光学系３０を備えるため、広いＦＯＶ１２Ａでの眼底撮影を可能とする。

【0028】

ＳＬＯユニット１８は、複数の光源、例えば、Ｂ光（青色光）の光源４０、Ｇ光（緑色光）の光源４２、Ｒ光（赤色光）の光源４４、およびＩＲ光（赤外線（例えば、近赤外光））の光源４６と、光源４０、４２、４４、４６からの光を、反射または透過して１つの光路に導く光学系４８、５０、５２、５４、５６とを備えている。光学系４８、５０、５６は、ミラーであり、光学系５２、５４は、ビームスプリッタ―である。Ｂ光は、光学系４８で反射し、光学系５０を透過し、光学系５４で反射し、Ｇ光は、光学系５０、５４で反射し、Ｒ光は、光学系５２、５４を透過し、ＩＲ光は、光学系５６、５２で反射して、それぞれ１つの光路に導かれる。

【0029】

ＳＬＯユニット１８は、Ｇ光、Ｒ光、およびＢ光を発するモードと、赤外線を発するモードなど、波長の異なるレーザ光を発する光源あるいは発光させる光源の組合せを切り替え可能に構成されている。図２に示す例では、Ｂ光（青色光）の光源４０、Ｇ光の光源４２、Ｒ光の光源４４、およびＩＲ光の光源４６の４つの光源を備えるが、本開示の技術は、これに限定されない。例えば、ＳＬＯユニット１８は、さらに、白色光の光源をさらに備え、白色光のみを発するモード等の種々のモードで光を発するようにしてもよい。

【0030】

ＳＬＯユニット１８から撮影光学系１９に入射された光は、第１光学スキャナ２２によってＸ方向およびＹ方向に走査される。走査光は広角光学系３０および瞳孔２７を経由して、被検眼１２の後眼部に照射される。眼底により反射された反射光は、広角光学系３０および第１光学スキャナ２２を経由してＳＬＯユニット１８へ入射される。

【0031】

ＳＬＯユニット１８は、被検眼１２の後眼部（例えば、眼底）からの光の内、Ｂ光を反射し且つＢ光以外を透過するビームスプリッタ６４、ビームスプリッタ６４を透過した光の内、Ｇ光を反射し且つＧ光以外を透過するビームスプリッタ５８を備えている。ＳＬＯユニット１８は、ビームスプリッタ５８を透過した光の内、Ｒ光を反射し且つＲ光以外を透過するビームスプリッタ６０を備えている。ＳＬＯユニット１８は、ビームスプリッタ６０を透過した光の内、ＩＲ光を反射するビームスプリッタ６２を備えている。

【0032】

ＳＬＯユニット１８は、複数の光源に対応して複数の光検出素子を備えている。ＳＬＯユニット１８は、ビームスプリッタ６４により反射したＢ光を検出するＢ光検出素子７０、およびビームスプリッタ５８により反射したＧ光を検出するＧ光検出素子７２を備えている。ＳＬＯユニット１８は、ビームスプリッタ６０により反射したＲ光を検出するＲ光検出素子７４、およびビームスプリッタ６２により反射したＩＲ光を検出するＩＲ光検出素子７６を備えている。

【0033】

広角光学系３０および第１光学スキャナ２２を経由してＳＬＯユニット１８へ入射された光（眼底により反射された反射光）は、Ｂ光の場合、ビームスプリッタ６４で反射してＢ光検出素子７０により受光され、Ｇ光の場合、ビームスプリッタ６４を透過し、ビームスプリッタ５８で反射してＧ光検出素子７２により受光される。上記入射された光は、Ｒ光の場合、ビームスプリッタ６４、５８を透過し、ビームスプリッタ６０で反射してＲ光検出素子７４により受光される。上記入射された光は、ＩＲ光の場合、ビームスプリッタ６４、５８、６０を透過し、ビームスプリッタ６２で反射してＩＲ光検出素子７６により受光される。ＣＰＵ１６Ａの制御下で動作する画像処理装置１６Ｇは、Ｂ光検出素子７０、Ｇ光検出素子７２、Ｒ光検出素子７４、およびＩＲ光検出素子７６で検出された信号を用いてＵＷＦ－ＳＬＯ画像を生成する。

【0034】

ＵＷＦ－ＳＬＯ画像（後述するようにＵＷＦ眼底画像、オリジナル眼底画像ともいう）には、眼底がＧ色で撮影されて得られたＵＷＦ－ＳＬＯ画像（Ｇ色眼底画像）と、眼底がＲ色で撮影されて得られたＵＷＦ－ＳＬＯ画像（Ｒ色眼底画像）とがある。ＵＷＦ－ＳＬＯ画像には、眼底がＢ色で撮影されて得られたＵＷＦ－ＳＬＯ画像（Ｂ色眼底画像）と、眼底がＩＲで撮影されて得られたＵＷＦ－ＳＬＯ画像（ＩＲ眼底画像）とがある。

【0035】

また、制御装置１６が、同時に発光するように光源４０、４２、４４を制御する。Ｂ光、Ｇ光およびＲ光で同時に被検眼１２の眼底が撮影されることにより、各位置が互いに対応するＧ色眼底画像、Ｒ色眼底画像、およびＢ色眼底画像が得られる。Ｇ色眼底画像、Ｒ色眼底画像、およびＢ色眼底画像からＲＧＢカラー眼底画像が得られる。制御装置１６が、同時に発光するように光源４２、４４を制御し、Ｇ光およびＲ光で同時に被検眼１２の眼底が撮影されることにより、各位置が互いに対応するＧ色眼底画像およびＲ色眼底画像が得られる。Ｇ色眼底画像およびＲ色眼底画像からＲＧカラー眼底画像が得られる。

【0036】

このようにＵＷＦ－ＳＬＯ画像として、具体的には、Ｂ色眼底画像、Ｇ色眼底画像、Ｒ色眼底画像、ＩＲ眼底画像、ＲＧＢカラー眼底画像、ＲＧカラー眼底画像がある。ＵＷＦ－ＳＬＯ画像の各画像データは、入力／表示装置１６Ｅを介して入力された患者の情報と共に、通信インターフェース（Ｉ／Ｆ）１６Ｆを介して眼科装置１１０からサーバ１４０へ送信される。ＵＷＦ－ＳＬＯ画像の各画像データと患者の情報とは、記憶装置２５４に、対応して記憶される。なお、患者の情報には、例えば、患者名ＩＤ、氏名、年齢、視力、右眼／左眼の区別等がある。患者の情報はオペレータが入力／表示装置１６Ｅを介して入力する。

【0037】

ＯＣＴシステムは、図２に示す制御装置１６、ＯＣＴユニット２０、および撮影光学系１９によって実現される。ＯＣＴシステムは、広角光学系３０を備えるため、上述したＳＬＯ眼底画像の撮影と同様に、広いＦＯＶ１２Ａでの眼底撮影を可能とする。ＯＣＴユニット２０は、光源２０Ａ、センサ（検出素子）２０Ｂ、第１の光カプラ２０Ｃ、参照光学系２０Ｄ、コリメートレンズ２０Ｅ、および第２の光カプラ２０Ｆを含む。

【0038】

光源２０Ａから射出された光は、第１の光カプラ２０Ｃで分岐される。分岐された一方の光は、測定光として、コリメートレンズ２０Ｅで平行光にされた後、撮影光学系１９に入射される。測定光は、第２光学スキャナ２４によってＸ方向およびＹ方向に走査される。走査光は広角光学系３０および瞳孔２７を経由して、眼底に照射される。眼底により反射された測定光は、広角光学系３０および第２光学スキャナ２４を経由してＯＣＴユニット２０へ入射され、コリメートレンズ２０Ｅおよび第１の光カプラ２０Ｃを介して、第２の光カプラ２０Ｆに入射する。

【0039】

光源２０Ａから射出され、第１の光カプラ２０Ｃで分岐された他方の光は、参照光として、参照光学系２０Ｄへ入射され、参照光学系２０Ｄを経由して、第２の光カプラ２０Ｆに入射する。

【0040】

第２の光カプラ２０Ｆに入射されたこれらの光、即ち、眼底で反射された測定光と、参照光とは、第２の光カプラ２０Ｆで干渉されて干渉光を生成する。干渉光はセンサ２０Ｂで受光される。ＣＰＵ１６Ａの制御下で動作する画像処理装置１６Ｇは、センサ２０Ｂで検出されたＯＣＴデータに基づいて断層画像やｅｎ－ｆａｃｅ画像などのＯＣＴ画像を生成する。

【0041】

ここで、内部照射角で１６０度以上の撮影画角で撮影されて得られたＯＣＴ眼底画像をＵＷＦ－ＯＣＴ画像と称する。もちろん、内部照射角で１６０度未満の撮影画角でＯＣＴ眼底画像データを取得することができる。

【0042】

ＵＷＦ－ＯＣＴ画像の画像データは、患者の情報と共に、通信インターフェース（Ｉ／Ｆ）１６Ｆを介して眼科装置１１０からサーバ１４０へ送信される。ＵＷＦ－ＯＣＴ画像の画像データと患者の情報とは、記憶装置２５４に、対応して記憶される。

【0043】

なお、本実施の形態では、光源２０Ａが波長掃引タイプのＳＳ－ＯＣＴ（Ｓｗｅｐｔ－Ｓｏｕｒｃｅ ＯＣＴ）を例示するが、ＳＤ－ＯＣＴ（Ｓｐｅｃｔｒａｌ－Ｄｏｍａｉｎ ＯＣＴ）、ＴＤ－ＯＣＴ（Ｔｉｍｅ－Ｄｏｍａｉｎ ＯＣＴ）など、様々な方式のＯＣＴシステムであってもよい。

【0044】

次に、眼軸長測定器１２０を説明する。眼軸長測定器１２０は、被検眼１２の眼軸方向の長さである眼軸長を測定する第１のモードと第２のモードとの２つのモードを有する。第１のモードは、図示しない光源からの光を被検眼１２に導光した後、眼底からの反射光と角膜からの反射光との干渉光を受光し、受光した干渉光を示す干渉信号に基づいて眼軸長を測定する。第２のモードは、図示しない超音波を用いて眼軸長を測定するモードである。

【0045】

眼軸長測定器１２０は、第１のモードまたは第２のモードにより測定された眼軸長をサーバ１４０に送信する。第１のモードおよび第２のモードにより眼軸長を測定してもよく、この場合には、双方のモードで測定された眼軸長の平均を眼軸長としてサーバ１４０に送信する。サーバ１４０は、患者の眼軸長を患者名ＩＤに対応して記憶する。

【0046】

次に、図３を参照して、サーバ１４０の電気系の構成を説明する。図３に示すように、サーバ１４０は、コンピュータ本体２５２を備えている。コンピュータ本体２５２は、バス２７０により相互に接続されたＣＰＵ２６２、ＲＡＭ２６６、ＲＯＭ２６４、および入出力（Ｉ／Ｏ）ポート２６８を有する。入出力（Ｉ／Ｏ）ポート２６８には、記憶装置２５４、ディスプレイ２５６、マウス２５５Ｍ、キーボード２５５Ｋ、および通信インターフェース（Ｉ／Ｆ）２５８が接続されている。記憶装置２５４は、例えば、不揮発メモリで構成される。入出力（Ｉ／Ｏ）ポート２６８は、通信インターフェース（Ｉ／Ｆ）２５８を介して、ネットワーク１３０に接続されている。従って、サーバ１４０は、眼科装置１１０、およびビューワ１５０と通信することができる。記憶装置２５４には、後述する画像処理プログラムが記憶されている。なお、画像処理プログラムを、ＲＯＭ２６４に記憶してもよい。

【0047】

画像処理プログラムは、本開示の技術の「プログラム」の一例である。記憶装置２５４、ＲＯＭ２６４は、本開示の技術の「メモリ」、「コンピュータ可読記憶媒体」の一例である。ＣＰＵ２６２は、本開示の技術の「プロセッサ」の一例である。

【0048】

サーバ１４０の後述する処理部２０８（図５も参照）は、眼科装置１１０から受信した各データを、記憶装置２５４に記憶する。具体的には、処理部２０８は記憶装置２５４に、ＵＷＦ－ＳＬＯ画像の各画像データおよびＵＷＦ－ＯＣＴ画像の画像データと患者の情報（上記のように患者名ＩＤ等）とを対応して記憶する。また、患者の被検眼に病変がある場合や病変部分に手術がされた場合には、眼科装置１１０の入力/表示装置１６Ｅを介して病変の情報が入力され、サーバ１４０に送信される。病変の情報は患者の情報と対応付けられて記憶装置２５４に記憶される。病変の情報には、病変部分の位置の情報、病変の名称、病変部分に手術がされている場合には手術名や手術日時等がある。

【0049】

ビューワ１５０は、ＣＰＵ、ＲＡＭ、ＲＯＭ等を備えたコンピュータとディスプレイとを備え、ＲＯＭには、画像処理プログラムがインストールされており、ユーザの指示に基づき、コンピュータは、サーバ１４０から取得した眼底画像などの医療情報が表示されるようにディスプレイを制御する。

【0050】

次に、図４を参照して、サーバ１４０のＣＰＵ２６２が画像処理プログラムを実行することで実現される各種機能について説明する。画像処理プログラムは、表示制御機能、画像処理機能（渦静脈解析機能、比較画像生成機能）、および処理機能を備えている。ＣＰＵ２６２がこの各機能を有する画像処理プログラムを実行することで、ＣＰＵ２６２は、図４に示すように、表示制御部２０４、画像処理部２０６（渦静脈解析部２０６０、比較画像生成部２０６２）、および処理部２０８として機能する。

【0051】

次に、図５を用いて、サーバ１４０による画像処理を詳細に説明する。サーバ１４０のＣＰＵ２６２が画像処理プログラムを実行することで、図５のフローチャートに示された画像処理が実現される。当該画像処理は、眼科装置１１０によりＵＷＦ眼底画像（ＵＷＦ－ＳＬＯ画像）が、取得され、患者名ＩＤと共にサーバ１４０に送信され、サーバ１４０が、患者名ＩＤおよびＵＷＦ眼底画像を受信した場合に、スタートする。

【0052】

ステップ５０２で、処理部２０８は、記憶装置２５４から、ＵＷＦ眼底画像Ｇ１（なお、図８には、ＲＧＢカラー眼底画像が示されている。）を取得する。

【0053】

ステップ５０４で、処理部２０８は、脈絡膜血管画像を、以下のようにして生成する。

【0054】

まず、ＵＷＦ眼底画像の中のＲ色眼底画像とＧ色眼底画像とに含まれる情報を説明する。

【0055】

眼の構造は、硝子体を、構造が異なる複数の層が覆うようになっている。複数の層には、硝子体側の最も内側から外側に、網膜、脈絡膜、強膜が含まれる。Ｒ光は、網膜を通過して脈絡膜まで到達する。よって、Ｒ色眼底画像には、網膜に存在する血管（網膜血管）の情報と脈絡膜に存在する血管（脈絡膜血管）の情報とが含まれる。これに対し、Ｇ光は、網膜までしか到達しない。よって、Ｇ色眼底画像には、網膜に存在する血管（網膜血管）の情報のみが含まれる。

【0056】

処理部２０８は、ブラックハットフィルタ処理をＧ色眼底画像に施すことにより、Ｇ色眼底画像から網膜血管を抽出する。次に、処理部２０８は、Ｒ色眼底画像から、Ｇ色眼底画像から抽出した網膜血管を用いてインペインティング処理により、網膜血管を除去する。つまり、Ｇ色眼底画像から抽出された網膜血管の位置情報を用いてＲ色眼底画像の網膜血管構造を周囲の画素と同じ値に塗りつぶす処理を行う。そして、処理部２０８は、網膜血管が除去されたＲ色眼底画像の画像データに対し、適応ヒストグラム均等化処理（ＣＬＡＨＥ、Ｃｏｎｔｒａｓｔ Ｌｉｍｉｔｅｄ Ａｄａｐｔｉｖｅ Ｈｉｓｔｏｇｒａｍ Ｅｑｕａｌｉｚａｔｉｏｎ）を施すことにより、Ｒ色眼底画像において、脈絡膜血管を強調する。これにより、図９に示すように、脈絡膜血管画像Ｇ２が得られる。生成された脈絡膜血管画像は記憶装置２５４に記憶される。

【0057】

また、Ｒ色眼底画像とＧ色眼底画像から脈絡膜血管画像を生成しているが、処理部２０８は、Ｒ色眼底画像Ｒ色眼底画像あるはＩＲ光で撮影されたＩＲ眼底画像を用いて脈絡膜血管画像を生成してもよい。

【0058】

脈絡膜眼底画像を生成する方法について、２０１８年３月２０日に出願された特願２０１８－０５２２４６の開示は、その全体が参照により、本明細書に取り込まれる。

【0059】

ステップ５０６で、渦静脈解析部２０６０は、脈絡膜血管画像を解析し、渦静脈を位置を検出する。そして、渦静脈位置と黄斑や視神経乳頭などの眼底構造物との位置関係を解析する。

【0060】

ここで、図６を参照して、ステップ５０６の渦静脈（Ｖｏｒｔｅｘ Ｖｅｉｎ）の検出処理を説明する。

【0061】

ステップ６０２で、渦静脈解析部２０６０は、記憶装置２５４から脈絡膜血管画像を読み込む。

【0062】

ステップ６０４で、渦静脈解析部２０６０は、渦静脈の位置を、以下のようにして検出する。

【0063】

渦静脈とは、脈絡膜に流れ込んだ血流の流出路であり、眼球の赤道部の後極寄りに４～６個存在する。渦静脈の位置は、脈絡膜血管の走行方向に基づいて算出される。これは、脈絡膜血管をたどると必ず流出路である渦静脈につながることにもとづく算出方法である。脈絡膜血管画像から血管の走行方向をたどり、複数の血管が合流する位置を渦静脈の位置とする。

【0064】

渦静脈解析部２０６０は、脈絡膜血管画像Ｇ２における各脈絡膜血管の移動方向（血管走行方向）を設定する。具体的には、第１に、渦静脈解析部２０６０は、脈絡膜血管画像の各画素について、下記の処理を実行する。即ち、渦静脈解析部２０６０は、画素に対して、当該画素を中心とした領域（セル）を設定し、セル内の各画素における輝度の勾配方向のヒストグラムを作成する。次に、渦静脈解析部２０６０は、各セルにおけるヒストグラムにおいて、最もカウントが少なかった勾配方向を各セルの内の画素における移動方向とする。この勾配方向が、血管走行方向に対応する。なお、最もカウントが少なかった勾配方向が血管走行方向であるとなるのは、次の理由からである。血管走行方向には輝度勾配が小さく、一方、それ以外の方向には輝度勾配が大きい（例えば、血管と血管以外のものでは輝度の差が大きい）。したがって、各画素の輝度勾配のヒストグラムを作成すると、血管走行方向に対するカウントは少なくなる。以上の処理により、脈絡膜血管画像の各画素における血管走行方向が設定される。

【0065】

渦静脈解析部２０６０は、Ｍ（自然数）×Ｎ（自然数）（＝Ｌ）個の仮想粒子の初期位置を脈絡膜血管画像上に設定する。具体的には、渦静脈解析部２０６０は、脈絡膜血管画像上に等間隔に、縦方向にＭ個、横方向にＮ個、合計Ｌ個の初期位置を設定する。

【0066】

渦静脈解析部２０６０は、渦静脈の位置を推定する。具体的には、渦静脈解析部２０６０は、Ｌ個の各々の位置について以下の処理を行う。即ち、渦静脈解析部２０６０は、最初の位置（Ｌ個の何れか）の血管走行方向を取得し、取得した血管走行方向に沿って所定距離だけ、仮想粒子を移動させ、移動した位置において、再度、血管走行方向を取得し、取得した血管走行方向に沿って所定距離だけ、仮想粒子を移動させる。このように血管走行方向に沿って所定距離移動させることを予め設定した移動回数、繰り返す。以上の処理を、Ｌ個の全ての位置において実行する。その時点で仮想粒子が一定個数以上集まっている点を渦静脈の位置とする。図１０には、４個の渦静脈ＶＶ１からＶＶ４が検出された様子が示されている。

【0067】

渦静脈の位置情報（渦静脈の個数や、脈絡膜血管画像上での座標など）は、記憶装置２５４に記憶される。

【0068】

ステップ６０６で、渦静脈解析部２０６０は、記憶装置２５４から、脈絡膜血管画像Ｇ２とＧ色眼底画像を読み込む。

【0069】

ステップ６０８で、渦静脈解析部２０６０は、黄斑と視神経乳頭との各々の位置（座標）を検出する。

【0070】

黄斑はＧ色眼底画像において暗い領域である。そこで、渦静脈解析部２０６０は、上記読み出したＧ色眼底画像において画素値が最も小さい所定数の画素の領域を黄斑の位置として検出する。

【0071】

渦静脈解析部２０６０は、Ｇ色眼底画像から視神経乳頭の位置を検出する。具体的には、渦静脈解析部２０６０は、上記読み出したＧ色眼底画像に対して、予め定まる視神経乳頭の画像のパターンマッチングをすることにより、Ｇ色眼底画像において視神経乳頭を検出する。また、視神経乳頭はＧ色眼底画像において最も明るい領域であるので、上記読み出したＧ色眼底画像において画素値が最も大きい所定数の画素の領域を視神経乳頭の位置として検出するようにしてもよい。

【0072】

ところで、脈絡膜血管画像は、上記のようにＲ色眼底画像およびＧ色眼底画像を処理することにより、作られる。従って、Ｇ眼底画像の座標系を脈絡膜血管画像の座標系に重ねると、Ｇ眼底画像の座標系の各位置は、脈絡膜血管画像の座標系の各位置と同じである。よって、Ｇ色眼底画像から検出された黄斑および視神経乳頭の各々の位置に相当する脈絡膜血管画像上の各位置は、黄斑および視神経乳頭の各々の位置である。

【0073】

よって、ステップ６０８の処理では、Ｇ色眼底画像に代えて脈絡膜血管画像から黄斑の位置を検出するようにしてもよい。同様にステップ６０８の処理では、Ｇ色眼底画像に代えて脈絡膜眼底画像から視神経乳頭の位置を検出するようにしてもよい。

【0074】

ステップ６１０で、渦静脈解析部２０６０は、上記検出された各渦静脈の相対位置（渦静脈と黄斑や視神経乳頭などの眼底構造物との位置関係）を算出する。

【0075】

渦静脈ＶＶ３は、図１４Ａに示すように、視神経乳頭ＯＮＨと黄斑Ｍとを結ぶ第１の線と視神経乳頭ＯＮＨと渦静脈（例えば、ＶＶ３）とを結ぶ第２の線とのなす角度θ、および、渦静脈ＶＶ３と視神経乳頭ＯＮＨとの距離ｒという位置関係になっている。ステップ６１０で、渦静脈解析部２０６０は、当該位置関係を算出する。また、渦静脈ＶＶ３は、図１５Ａに示すように、視神経乳頭ＯＮＨと黄斑Ｍとを結ぶ第１の線と第１の線と黄斑Ｍと渦静脈ＶＶ３とを結ぶ第３の線とのなす角度φ、および、渦静脈と黄斑との距離ｓという位置関係になっている。このような位置関係の算出を検出され各渦静脈に対して行う。

【0076】

角度θおよび距離ｒの組、角度φおよび距離ｓの組は、本開示の技術の「第１の組」および「第２の組」の一例である。

【0077】

渦静脈解析部２０６０は、各渦静脈について、角度θ、距離ｒ、角度φ、および、距離ｓを計算する。なお、角度φおよび距離ｓの計算方法は、角度θおよび距離ｒの計算方法と、視神経乳頭ＯＮＨに代えて黄斑Ｍを基準にする点が異なり、略同様である。そこで、以下、角度θおよび距離ｒの計算方法を説明し、角度φおよび距離ｓの計算方法の説明を省略する。

【0078】

まず、図７Ａを参照して、渦静脈（以下、ＶＶという）と視神経乳頭ＯＮＨとの距離（ＶＶ距離という）ｒの計算方法を説明する。

【0079】

ステップ２２１で、渦静脈解析部２０６０は、脈絡膜眼底画像において視神経乳頭ＯＮＨ、黄斑Ｍ、ＶＶ位置の各座標を取得する。

【0080】

次に、ステップ２２３で、渦静脈解析部２０６０は、視神経乳頭ＯＮＨ、黄斑Ｍ、ＶＶ位置の各座標を図１３に示す仮想球面に投影する。図１３に示す仮想球面は、眼球の中心をＣ、半径をＬ（眼軸長を２Ｌとする）とした球面である。この球面に、ＶＶの位置をＶ、視神経乳頭ＯＮＨの位置をＯ、黄斑の位置をＭとして投影する。

【0081】

この仮想球面を眼球モデルと考え、ステップ２２５で、渦静脈解析部２０６０は、球面上の２点間の大円距離をＶＶ距離ｒとして計算する。つまり、球の中心Ｏを通るように球を切ったときの切り口を大円と定義し、大円距離とは、球面上の距離計測対象の２地点（ＶＶ位置：Ｖと、視神経乳頭位置：Ｏ）を結ぶ大円の弧の長さと定義する。ＶＶ位置：Ｖの仮想球面上での緯度経度を（緯度θ１，経度φ１），視神経乳頭位置：Ｏの緯度経度（緯度θ２，経度φ２）とすると、本ステップ２２３で、渦静脈解析部２０６０は、ＶＶ位置と視神経乳頭位置との間のＶＶ距離ｒ、すなわち大円距離ＯＶを、球面三角法の公式から、計算する。

【0082】

【数1】

【0083】

ステップ２２７で、渦静脈解析部２０６０は、算出されたＶＶ距離ｒ、即ち、ＶＶ位置と視神経乳頭位置との間のＶＶ距離ｒ（大円距離ＯＶ）を、記憶装置２５４に記憶する。

【0084】

図７Ａに示す処理が終了すると、相対位置算出処理は、図７Ｂに示す角度θ（以下、ＶＶ角度という）算出処理に進む。

【0085】

ＶＶ角度は、図１１および図１４Ａに示すように、黄斑Ｍの位置－視神経乳頭ＯＮＨの位置－ＶＶの位置がなす角度θである。

【0086】

黄斑Ｍの位置－視神経乳頭ＯＮＨの位置－ＶＶの位置がなす角度θの算出方法には、次の正角図法から算出する方法と、球面三角法から算出する方法とがある。

【0087】

最初に、図７Ｂを参照して、角度θを、正角図法から算出する方法を説明する。

【0088】

通常の内積から算出する方法では、算出された角度の正負の区別がつかず、上半球にあるＶＶ（図１１の（ｘ_３、ｙ_３））なのか、下半球にあるＶＶ‘（（ｘ_３、－ｙ_３））なのかを区別できない。また逆正接関数を用いる方法では、正負の区別はつくが、θの算出方向を常に一定（例えば反時計回り）にしてしまうので、解剖学的特徴（鼻側・耳側）を基準にした場合、左右眼で上下半球の値が反転してしまう。そこで、本実施の形態では、左右眼調整符号ｆ_ｓｉｇｎを用いて、算出した角度の正負を調整するようにしている。

【0089】

図７Ｂのステップ２２２で、渦静脈解析部２０６０は、左右眼調整符号ｆ_ｓｉｇｎを算出する。図１１に示すように、黄斑Ｍの位置を（ｘ_１、ｙ_１）、視神経乳頭ＯＮＨの位置を（ｘ_２、ｙ_２）、ＶＶの位置を（ｘ_３、ｙ_３）とする。

【0090】

左右眼調整符号ｆ_ｓｉｇｎは、

【0091】

ｆ_ｓｉｇｎ＝＋１ （ｘ_１＞ｘ_２の場合）

【0092】

ｆ_ｓｉｇｎ＝－１ （ｘ_１＜ｘ_２の場合）

【0093】

と設定さる。これは、黄斑と視神経乳頭の位置から解剖学的に、ｘ_１＞ｘ_２となるのは左眼であり、ｘ_１＜ｘ_２となるのは右眼と判断できるからである。また、仮にｘ_１＝ｘ_２となる場合は、ｆ_ｓｉｇｎ＝＋１とする。

【0094】

ステップ２２４で、渦静脈解析部２０６０は、ｃｏｓθおよびｓｉｎθを、ベクトルの内積と外積の定義に基づいて、数式２と数式３を用いて算出する。ここで、斑Ｍの位置－視神経乳頭ＯＮＨの位置－ＶＶの位置がなす角度θはベクトルＯＭ（視神経乳頭位置Ｏと黄斑位置Ｍを結ぶベクトル）とベクトルＯＶ（視神経乳頭位置Ｏと渦静脈位置Ｖを結ぶベクトル）とがなす角である。

【0095】

【数2】

【0096】

【数3】

【0097】

ステップ２２６で、渦静脈解析部２０６０は、４象限逆正接関数によるθを、次のように算出する。

【0098】

【数4】

【0099】

４象限逆正接関数によって求められたθは、図１２に示すように、ｙ／ｘの値だけでなく、４象限の各々のｘの符号も考慮されている。

【0100】

ステップ２２８で、渦静脈解析部２０６０は、算出されたθの符号を、次のようにｆ_ｓｉｇｎを用いて左右眼補正をする。ステップ２３０で、渦静脈解析部２０６０は、このように求められたθの値をＶＶ角度として記憶装置２５４に記憶する。

【0101】

【数5】

【0102】

次に、角度θを球面三角法から算出する方法を説明する。

【0103】

図１３に示すように、黄斑の位置Ｍ、視神経乳頭ＯＮＨの位置Ｏ、ＶＶの位置Ｖは、眼球中心Ｃを中心とし、眼軸長を２Ｌとしたときの半径Ｌの球の表面に位置する。黄斑の位置Ｍ、視神経乳頭ＯＮＨの位置Ｏ、ＶＶの位置Ｖを頂点とする三角形ＯＭＶの頂角Ｏの角度をαとすると、

【0104】

【数6】

【0105】

【数7】

【0106】

となり、α（ただし、αは［０，π］の範囲内）を算出できる。

【0107】

αの値を、θ（開区間［－π，π］の範囲内）に変更するには、ｇ_ｓｉｇｎ（黄斑Ｍと視神経乳頭ＯＮＨの位置関係，ＶＶの位置Ｖと黄斑Ｍの位置関係）＝｛１，－１｝を考慮し、θ＝α・ｇ_ｓｉｇｎで算出する。このように求められたθの値をＶＶ角度としてメモリ１６４に記憶する。

【0108】

画像処理プログラムには、図７Ｂに示した正角図法からＶＶ角度を算出するプログラムか、上記の球面三角法から算出するプログラムか、の何れか一方を備えている。

【0109】

図７Ｂに示した算出処理が終了すると、図６の渦静脈の解析処理（図５のステップ５０６）が終了し、画像処理は、ステップ５０８に進む。

【0110】

ステップ５０８で、処理部２０８は、患者名ＩＤに対応する過去画像に関するデータを読み出す。患者名ＩＤに対応する過去画像に関するデータは、患者名ＩＤに対応する過去のＵＷＦ眼底画像に基づく渦静脈の相対位置である。

【0111】

過去画像に関するデータは、今回とは異なるタイミングにおいて図５の画像処理が過去に実行されて記憶されているデータである。

【0112】

図５の画像処理が今回実行されて得られた今回のＵＷＦ眼底画像に基づく渦静脈の相対位置と、図５の画像処理が過去に実行されて得られた過去画像に基づく渦静脈の相対位置とは、本開示の技術の「第１位置」及び「第２位置」の一例である。
なお、本開示の技術の「第１位置」及び「第２位置」は、上記相対位置に限定されず、渦静脈の位置する座標で特定されてもよい。

【0113】

ステップ５１０で、比較画像生成部２０６２は、渦静脈マップのデータを生成する。

【0114】

ここで、渦静脈マップは、異なるタイミング（日時）で撮影された複数の眼底画像から抽出した渦静脈の位置をプロットすることにより得られるマップである。このように異なるタイミングで撮影された複数の眼底画像から抽出した渦静脈の位置がマップにプロットされているので、マップを見ることにより、眼科医等は、新たに撮影されたＵＷＦ眼底画像に基づく渦静脈の相対位置と過去のＵＷＦ眼底画像に基づく渦静脈の相対位置とを比較することができる。

【0115】

今回のＵＷＦ眼底画像に基づく渦静脈の相対位置は、画像処理が今回実行され、今回取得されたＵＷＦ眼底画像に基づく渦静脈の相対位置である。過去のＵＷＦ眼底画像に基づく渦静脈の相対位置は、ステップ５０８で読み出された、患者名ＩＤに対応する過去のＵＷＦ眼底画像に基づく渦静脈の相対位置である。

【0116】

上記のように渦静脈の相対位置は、図１４Ａに示すように、視神経乳頭ＯＮＨを基準とした角度θおよび距離ｒと、図１５Ａに示すように、黄斑Ｍを基準とした角度φおよび距離ｓにより、特定される。

【0117】

角度θおよび距離ｒを用いた渦静脈マップおよび角度φおよび距離ｓを用いた渦静脈マップの各々には、複数のマップがある。

【0118】

例えば、角度θおよび距離ｒを用いた渦静脈マップには、第１に、今回および過去の渦静脈の位置を、図１４Ｂに示すように、角度θおよび距離ｒを用いた極座標系で表示したマップ（レーダチャート）がある。渦静脈マップには、第２に、今回および過去の渦静脈の位置を、図１４Ｃに示すように、角度θおよび距離ｒを用いた２次元直交座標系で表示したマップがある。渦静脈マップには、第３に、今回および過去の渦静脈の位置を、図１４Ｄに示すように、角度θおよび距離ｒを用いた３次元直交座標系で表示したマップがある。
図１４Ｃに示す２次元直交座標系で表示したマップは、本開示の技術の「２次元直交軸で規定される散布図」の一例であり、図１４Ｄに示す３次元直交座標系で表示したマップは、本開示の技術の「３次元直交軸で規定される散布図」の一例である。
比較画像生成部２０６２は、図１４Ｂから図１４Ｄの各マップのデータを、次のように、生成する。

【0119】

図１４Ｂから図１４Ｄの各マップは、今回のＵＷＦ眼底画像に基づく渦静脈の相対位置が、「・」でプロットして、生成されており、過去３回の眼底撮影で得られたＵＷＦ眼底画像に基づく渦静脈の相対位置が、「〇」、「□」、および「◎」でプロットして、作成されている。

【0120】

図１４Ｄに示す３次元直交座標系のマップでは、眼球中心は、３次元（Ｘ，Ｙ，Ｚ）の中心（０．０．０）であり、頂点（０．０．１）は視神経乳頭ＯＮＨの位置である。３次元直交座標系のマップでは、角度θおよび距離ｒで特定された渦静脈の位置が、眼球モデルの、角度θおよび距離ｒに対応する位置にプロットされている。眼球モデルは、眼軸長により、大きさが補正されている。

【0121】

図１４Ｂから図１４Ｄの各マップに示す例では、患者名ＩＤで識別される患者の渦静脈の位置が変化していることが示されている。

【0122】

一方、角度φおよび距離ｓを用いた渦静脈マップは、図１５Ｂから図１５Ｄに示されている。図１５Ｂから図１５Ｄの各マップそれぞれは、図１４Ｂから図１４Ｄの各マップに対応するので、その説明を省略する。

【0123】

ステップ５１２で、処理部２０８は、ＵＷＦ眼底画像、脈絡膜眼底画像、各渦静脈の位置（角度θ、距離ｒ、角度φ、および距離ｓ）、および各マップのデータを、患者名ＩＤに対応して記憶装置２５４に保存する。

【0124】

ところで、図５の画像処理は、サーバ１４０が、眼科装置１１０から、患者名ＩＤおよびＵＷＦ眼底画像を受信する毎に実行している。本開示の技術は、これに限定されず、例えば、画像処理は、オペレータが患者名ＩＤを入力し、図示しないスタートボタンが操作された場合に実行するようにしてもよい。サーバ１４０が、眼科装置１１０から、ＵＷＦ眼底画像を受信しても、渦静脈の相対位置が算出されていないＵＷＦ眼底画像が存在する態様がある。

【0125】

この場合、第１に、ステップ５０２では、処理部２０８は、記憶装置２５４から、渦静脈の相対位置が算出されていないＵＷＦ眼底画像Ｇ１を取得する。渦静脈の相対位置が算出されていない複数のＵＷＦ眼底画像Ｇ１が取得された場合、取得された各ＵＷＦ眼底画像Ｇ１について、ステップ５０４、５０６を実行する。

【0126】

第２に、渦静脈の相対位置が算出されているか否かに拘わらず、患者名ＩＤに対応する全てのＵＷＦ眼底画像Ｇ１について、ステップ５０４、５０６を実行する。

【0127】

眼科医は、患者の被検眼を診察する際、患者名ＩＤをビューワ１５０に入力する。ビューワ１５０はサーバ１４０に、患者名ＩＤに対応する各データを送信するように指示する。サーバ１４０は、患者名ＩＤに対応する、患者名、患者の年齢、患者の視力、左眼か右眼かの情報、眼軸長、撮影日、および各データを、患者名ＩＤと共に、ビューワ１５０に送信する。各データは、上記ＵＷＦ眼底画像、脈絡膜眼底画像、各渦静脈の位置（角度θ、距離ｒ、角度φ、および距離ｓ）、および各マップの画面のデータが含まれる。

【0128】

患者名ＩＤ、患者名、患者の年齢、患者の視力、左眼か右眼かの情報、眼軸長、撮影日、およびの各データを受信したビューワ１５０は、図１６に示す第１眼底画像表示画面１０００Ａをディスプレイに表示する。

【0129】

図１６に示すように、第１眼底画像表示画面１０００Ａは、患者情報表示欄１００２と、第１眼底画像情報表示欄１００４Ａとを備えている。
第１眼底画像情報表示欄１００４Ａは、本開示の技術の「表示ビュー」の一例である。

【0130】

患者情報表示欄１００２は、患者名ＩＤ、患者名、患者の年齢、患者の視力、左眼か右眼かの情報、眼軸長を表示するための、各々の表示欄１０１２から１０２２と、画面切替ボタン１０２４とを有する。表示欄１０１２から１０２２に、受信した患者名ＩＤ、患者名、患者の年齢、患者の視力、左眼か右眼かの情報、眼軸長が表示される。

【0131】

第１眼底画像情報表示欄１００４Ａは、ＵＷＦ眼底画像表示欄１０３２Ａ、脈絡膜血管画像表示欄１０３２Ｂ、および現在および過去の脈絡膜血管画像表示欄１０３５ＢＴ４から１０３５ＢＴ１、スライダー１０３７、およびスライダーバー１０３５を備えている。第１眼底画像情報表示欄１００４Ａは、渦静脈マップ表示欄１０３２Ｃ、ＶＶマップ切替ボタン１０３６、および情報表示欄１０３４を備えている。

【0132】

脈絡膜血管画像表示欄１０３５ＢＴ４から１０３５ＢＴ１には、眼底が撮影された撮影日時（ＹＹＹ/ＭＭ/ＤＤ）のデータが表示される。

【0133】

ＵＷＦ眼底画像表示欄１０３２Ａおよび脈絡膜血管画像表示欄１０３２Ｂには、最初は、現在に一番近い時に眼底が撮影されて得られたＵＷＦ眼底画像および脈絡膜血管画像が表示される。ＵＷＦ眼底画像、即ち、オリジナル眼底画像としては、例えば、ＲＧＢカラー眼底画像である。

【0134】

眼科医がスライダー１０３７を移動させると、ＵＷＦ眼底画像表示欄１０３２Ａおよび脈絡膜血管画像表示欄１０３２Ｂには、スライダー１０３７が位置する撮影日時に眼底が撮影されて得られたＵＷＦ眼底画像および脈絡膜血管画像が表示される。

【0135】

渦静脈マップ表示欄１０３２Ｃには、渦静脈マップが表示される。具体的には、ＶＶマップ切替ボタン１０３６が操作されると、渦静脈マップ（図１４Ｂから図１４Ｄおよび図１５Ｂから図１５Ｄ）を選択するためのプルダウンメニュー１０３６ＰＤＭが表示される。

【0136】

プルダウンメニュー１０３６ＰＤＭには、「レーダチャート 視神経乳頭中心」、「レーダチャート 黄斑中心」、「プロット図 視神経乳頭中心」、「プロット図 黄斑中心」、「３Ｄ 視神経乳頭中心」、「３Ｄ 黄斑中心」の６個の項目が含まれる。

【0137】

「レーダチャート 視神経乳頭中心」が選択されると、図１４Ｂに示すように、視神経乳頭ＯＮＨを基準とした角度θおよび距離ｒを用いた極座標系で表示したマップ（レーダチャート）が表示されることが選択される。

【0138】

「レーダチャート 黄斑中心」が選択されると、図１５Ｂに示すように、黄斑Ｍを基準とした角度θおよび距離ｒを用いた極座標系で表示したマップ（レーダチャート）が表示されることが選択される。

【0139】

「プロット図 視神経乳頭中心」が選択されると、図１４Ｃに示すように、角度θおよび距離ｒを用いた２次元直交座標系で表示したマップが表示されることが選択される。

【0140】

「プロット図 黄斑中心」が選択されると、図１５Ｃに示すように、角度θおよび距離ｒを用いた２次元直交座標系で表示したマップが表示されることが選択される。

【0141】

「３Ｄ 視神経乳頭中心」が選択されると、図１４Ｄに示すように、角度θおよび距離ｒを用いた３次元直交座標系で表示したマップが表示されることが選択される。

【0142】

「３Ｄ 黄斑中心」が選択されると、図１５Ｄに示すように、角度θおよび距離ｒを用いた３次元直交座標系で表示したマップが表示されることが選択される。

【0143】

例えば、プルダウンメニュー１０３６ＰＤＭの中で、視神経乳頭を中心としたレーダチャート（「レーダチャート 視神経乳頭中心」の項目）が選択された場合、渦静脈マップ表示欄１０３２Ｃには、スライダー１０３７が位置する撮影日時に対応する、視神経乳頭を中心としたレーダチャート（図１４Ｂ参照）が表示される。

【0144】

プルダウンメニュー１０３６ＰＤＭの項目が選択されると、対応する渦静脈マップが表示され、表示された渦静脈マップには、複数（上記例では４個）の全ての渦静脈の相対位置が表示される。本開示の技術はこれに限定されない。渦静脈マップに、複数の渦静脈の中から選択された１つの渦静脈の相対位置が表示される。例えば、図１６に示す例では、例えば、ＵＷＦ眼底画像表示欄１０３２Ａに表示されている４個の渦静脈の中から、眼科医が、１つの渦静脈を選択する。例えば、左上の渦静脈が選択されると、渦静脈マップ表示欄１０３２Ｃには、選択された渦静脈のみの相対位置が表示される。なお、本開示の技術は、４個の渦静脈の中から、眼科医が、単一の渦静脈を選択することに限定されない。例えば、ビューワ１５０が、４個の渦静脈について予め定められた順、例えば、左上、右上、右下、左下の順で、１つずつ渦静脈を選択し、渦静脈マップ表示欄１０３２Ｃに、選択された渦静脈のみの相対位置が表示される。あるいは、２つの渦静脈を選択し、選択した２つの渦静脈の相対位置が表示されるようにしてもよい。例えば、左上と左下の２つの渦静脈をユーザが選択することにより、黄斑と視神経乳頭を結ぶ線分に対してと左上の渦静脈の位置と左下の渦静脈の位置が、経過観察期間において対称性を維持しているか、ある時点から対称性が崩れ非対称となったかを簡単に判別できる。

【0145】

情報表示欄１０３４には、ユーザ（眼科医）の診察時のコメントやメモがテキストとして表示される。

【0146】

現在および過去の脈絡膜血管画像表示欄１０３５ＢＴ４から１０３５ＢＴ１の大きさは、ＵＷＦ眼底画像表示欄１０３２Ａおよび脈絡膜血管画像表示欄１０３２Ｂの大きさより小さい。

【0147】

渦静脈マップ表示欄１０３２Ｃの大きさは、ＵＷＦ眼底画像表示欄１０３２Ａおよび脈絡膜血管画像表示欄１０３２Ｂの大きさより小さい。

【0148】

図１６の画面切替ボタン１０２４が操作された場合、ビューワ１５０は、図１７に示す第２眼底画像表示画面１０００Ｂをディスプレイに表示する。

【0149】

第１眼底画像表示画面１０００Ａと第２眼底画像表示画面１０００Ｂとは略同様の内容であるので、同一の部分には同一の符号を付してその説明を省略し、異なる部分のみ説明する。

【0150】

第２眼底画像表示画面１０００Ｂの眼底画像情報表示欄１００４Ｂの渦静脈マップ表示欄１０３２Ｃの大きさは、図１６の渦静脈マップ表示欄１０３２Ｃの大きさより大きい。第２眼底画像表示画面１０００Ｂの渦静脈マップ表示欄１０３２Ｃの位置は、図１６の渦静脈マップ表示欄１０３２Ｃの位置と異なる。

【0151】

図１７の渦静脈マップ表示欄１０３２Ｃにおいて、何れかのプロットがカーソル（黒→）で選択されると、ビューワ１５０は、当該プロットの色を他の色に変化させ、同じ脈絡膜血管画像から検出された他の渦静脈の位置のプロットの色も当該他の色に変化させる。また、ビューワ１５０は、当該プロットに対応する撮影日時に対応するように、スライダー１０３７を移動させる。ビューワ１５０は、ＵＷＦ眼底画像表示欄１０３２Ａおよび脈絡膜血管画像表示欄１０３２Ｂには、スライダー１０３７が位置する撮影日時に眼底が撮影されて得られたＵＷＦ眼底画像および脈絡膜血管画像を表示する。

【0152】

図１７の画面切替ボタン１０２４が操作された場合、ビューワ１５０は、図１８に示す第２眼底画像表示画面１０００Ｃをディスプレイに表示する。

【0153】

第２眼底画像表示画面１０００Ｂと第３眼底画像表示画面１０００Ｃとは略同様の内容であるので、同一の部分には同一の符号を付してその説明を省略し、異なる部分のみ説明する。

【0154】

第３眼底画像表示画面１０００Ｃの第３眼底画像情報表示欄１００４Ｃは、渦静脈マップ表示欄１０３８、カラーバー１０３５Ｃ、および情報表示欄１０３４のみを備えている。

【0155】

渦静脈マップ表示欄１０３８には、例えば、今回および過去の渦静脈の位置を、黄斑を中心とした２次元直交座標系で表示したマップ（図１５Ｃ）が表示される。何れかのプロット（例えば、「・」）がカーソル（黒→）で選択されると、ビューワ１５０は、当該プロットの色を他の色に変化させ、同じ脈絡膜血管画像から検出された他の渦静脈の位置のプロットの色も当該他の色に変化させる。ビューワ１５０は、渦静脈マップ表示欄１０３８の左上に、選択されたプロット（例えば、「・」）に対応する撮影日時を表示する。ビューワ１５０は、カラーバー１０３５Ｃの、現在（左端）から「・」の撮影日時まで色を変化させる。

【0156】

以上説明したように本実施の形態では、今回から過去の各眼底画像から得られた渦静脈の位置を比較する比較マップを表示するので、渦静脈の位置の変化をわかりやすく表示することができる。よって、渦静脈位置の経過観察を行うことができ、渦静脈位置の変化をユーザにわかりやすく表示できる。

【0157】

ところで、図１６および図１７の渦静脈マップ表示欄１０３２Ｃおよび図１８の渦静脈マップ表示欄１０３８に表示される渦静脈マップには、渦静脈の位置のみがプロットされる。本開示の技術は、これに限定されず、ビューワ１５０は、更に次の情報も表示するようにしてもよい。例えば、図１９に示すように、ビューワ１５０は、第１に、各撮影日時の渦静脈の４つの位置の中心Ｃ１からＣ４、第２に、中心Ｃ１からＣ４の中心ＣＣを表示する。

【0158】

更に、本開示の技術ではさらに、各撮影日時の各渦静脈までの距離ｒ、ｓの時間変化を示すグラフを表示してもよい。
また、上記実施の形態では、渦静脈の位置を、角度θおよび距離ｒの組、角度φおよび距離ｓにより特定しているが、本開示の技術は、これに限定されない。例えば、脈絡膜眼底画像の各画素の位置と、３次元空間の位置（Ｘ，Ｙ，Ｚ）との対応関係を予め求めておく。脈絡膜血管画像から特定された渦静脈の位置から、当該対応関係から、渦静脈の３次元空間の位置を特定し、特定された渦静脈の３次元空間の位置から、上記角度及び距離を計算する。

【0159】

以上説明した例では、図５の画像処理は、サーバ１４０が実行しているが、本開示の技術の技術はこれに限定されない。例えば、眼科装置１１０またはビューワ１５０が実行したり、ネットワーク１３０に更に別の画像処理装置を接続させ、当該画像処理装置が実行したり、してもよい。

【0160】

以上説明した各例では、コンピュータを利用したソフトウェア構成により画像処理が実現される場合を例示したが、本開示の技術はこれに限定されるものではない。例えば、コンピュータを利用したソフトウェア構成に代えて、ＦＰＧＡ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）またはＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）等のハードウェア構成のみによって、画像処理が実行されるようにしてもよい。画像処理のうちの一部の処理がソフトウェア構成により実行され、残りの処理がハードウェア構成によって実行されるようにしてもよい。

【0161】

このように本開示の技術は、コンピュータを利用したソフトウェア構成により画像処理が実現される場合と、コンピュータを利用したソフトウェア構成でない構成で画像処理が実現される場合とを含むので、以下の第１技術および第２技術を含む。

【0162】

（第１技術）
第１眼底画像から渦静脈の第１位置を特定し且つ第２眼底画像から前記渦静脈の第２位置を特定する特定部と、
前記第１位置および前記第２位置が表示される画面のデータを生成する生成部と、
を含む画像処理装置。

【0163】

なお、上記実施の形態の渦静脈解析部２０６０は、上記第１技術の「特定部」の一例であり、上記実施の形態の比較画像生成部２０６２は、上記第１技術の「生成部」の一例である。

【0164】

上記のように以上の開示内容から以下の第２技術が提案される

【0165】

（第２技術）
特定部が、第１眼底画像から渦静脈の第１位置を特定し且つ第２眼底画像から前記渦静脈の第２位置を特定し、
生成部が、前記第１位置および前記第２位置が表示される画面のデータを生成する、
画像処理方法。

【0166】

以上の開示内容から以下の第３技術が提案される。

【0167】

（第３技術）
画像処理するためのコンピュータープログラム製品であって、
前記コンピュータープログラム製品は、それ自体が一時的な信号ではないコンピュータ可読記憶媒体を備え、
前記コンピュータ可読記憶媒体には、プログラムが格納されており、
前記プログラムは、
コンピュータに、
第１眼底画像から渦静脈の第１位置を特定し、
第２眼底画像から前記渦静脈の第２位置を特定し、
前記第１位置および前記第２位置が表示される画面のデータを生成する、
ことを実行させる、
コンピュータープログラム製品。

【0168】

以上説明した画像処理はあくまでも一例である。従って、主旨を逸脱しない範囲内において不要なステップを削除したり、新たなステップを追加したり、処理順序を入れ替えたりしてもよいことは言うまでもない。

【0169】

本明細書に記載された全ての文献、特許出願、および技術規格は、個々の文献、特許出願、および技術規格が参照により取り込まれることが具体的にかつ個々に記載された場合と同様に、本明細書中に参照により取り込まれる。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7A】

【図7B】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14A】

【図14B】

【図14C】

【図14D】

【図15A】

【図15B】

【図15C】

【図15D】

【図16】

【図17】

【図18】

【図19】