WO2019102740A1

WO2019102740A1 - 眼科情報処理装置、眼科システム、眼科情報処理方法、及びプログラム

Info

Publication number: WO2019102740A1
Application number: PCT/JP2018/038121
Authority: WO
Inventors: 寛一 ▲徳▼田; 太樹相見
Original assignee: 株式会社トプコン
Priority date: 2017-11-24
Filing date: 2018-10-12
Publication date: 2019-05-31
Also published as: JP2019092920A; US11832884B2; EP3714768A4; JP7126818B2; EP3714768A1; US20200275834A1

Abstract

眼科情報処理装置は、解析部と、表示制御部とを含む。解析部は、光コヒーレンストモグラフィを用いて収集された被検眼の眼底のデータを解析することにより眼底における１以上の萎縮領域を特定する。表示制御部は、収集タイミングが異なる複数のデータに基づいて、解析部により特定された複数の１以上の萎縮領域を表す画像を時系列に表示手段に表示させる。

Description

眼科情報処理装置、眼科システム、眼科情報処理方法、及びプログラム

　本発明は、眼科情報処理装置、眼科システム、眼科情報処理方法、及びプログラムに関する。

　視覚障害の原因疾患の１つに、加齢黄斑変性（Ａｇｅ－ｒｅｌａｔｅｄ　Ｍａｃｕｌａｒ　Ｄｅｇｅｎｅｒａｔｉｏｎ：ＡＭＤ）がある。加齢黄斑変性は、加齢により直接的に又は間接的に黄斑部が障害される疾患である。加齢黄斑変性は、滲出型加齢黄斑変性（ｅｘｕｄａｔｉｖｅ　ＡＭＤ）と萎縮型加齢黄斑変性（ａｔｒｏｐｈｉｃ　ＡＭＤ）とに分類される。滲出型加齢黄斑変性は、脈絡膜から脈絡膜新生血管が網膜色素上皮層（Ｒｅｔｉｎａｌ　Ｐｉｇｍｅｎｔ　Ｅｐｉｔｈｅｌｉｕｍ：以下、ＲＰＥ）の下層に侵入したり網膜とＲＰＥとの間に侵入したりすることにより網膜が障害される疾患である。萎縮型加齢黄斑変性は、ＲＰＥが徐々に萎縮することで網膜が障害され、視力が徐々に低下する疾患である。

　滲出型加齢黄斑変性に対する有効な治療方法として、光線力学的療法（Ｐｈｏｔｏ　Ｄｙｎａｍｉｃ　Ｔｈｅｒａｐｙ：ＰＤＴ）、薬物療法、レーザー凝固などが知られている。これに対して、萎縮型加齢黄斑変性に対する有効な治療方法は、現在のところ十分に確立されていない。従って、萎縮型加齢黄斑変性の状態を把握することは極めて重要である。

　萎縮型加齢黄斑変性には、中心窩を中心とする所定の領域内にいわゆる地図状萎縮（Ｇｅｏｇｒａｐｈｉｃ　Ａｔｒｏｐｈｙ：ＧＡ）が見られる。地図状萎縮は、眼底画像、フルオレセイン蛍光眼底造影画像、眼底自発蛍光検査画像などから特定されたり、光コヒーレンストモグラフィを用いて取得された網膜の断層像から特定されたりする（例えば、特許文献１～特許文献３）。特定された地図状萎縮を観察することで、萎縮型加齢黄斑変性の状態を把握することができる（例えば、特許文献４）。また、例えば、特許文献５及び特許文献６には、光コヒーレンストモグラフィを用いて取得されたデータを解析して得られた層厚データのトレンド解析結果を表示する手法が開示されている。

米国特許出願公開第２０１５／０２０１８２９号明細書特開２０１５－１３６６２６号公報特開２０１６－１０７１４８号公報特表２０１４－５０５５５２号公報特開２０１４－０８３２６６号公報特開２０１４－０８３２６８号公報

　萎縮型加齢黄斑変性の状態の把握には、地図状萎縮を伴う領域（地図状萎縮領域）の形態（形状、サイズ）や分布を経過観察することが有効である。しかしながら、従来の技術では、地図状萎縮領域の形態や分布を詳細に観察することができず、地図状萎縮領域の形態や分布の時間的変化を容易に把握することができなかった。それにより、萎縮型加齢黄斑変性に対して適確な診断を行うことが困難であった。

　本発明は、このような事情を鑑みてなされたものであり、その目的は、眼底における萎縮領域の形態や分布を容易に経過観察するための新たな技術を提供することにある。

　いくつかの実施形態の第１態様は、光コヒーレンストモグラフィを用いて収集された被検眼の眼底のデータを解析することにより前記眼底における１以上の萎縮領域を特定する解析部と、収集タイミングが異なる複数の前記データに基づいて、前記解析部により特定された複数の前記１以上の萎縮領域を表す画像を時系列に表示手段に表示させる表示制御部と、を含む眼科情報処理装置である。

　いくつかの実施形態の第２態様は、第１態様において、前記被検眼の眼底画像と前記１以上の萎縮領域を表す画像との位置合わせを行う位置合わせ処理部を含み、前記表示制御部は、前記位置合わせ処理部により位置合わせが行われた前記１以上の萎縮領域を表す画像が前記眼底画像にオーバーレイされた複数の萎縮領域画像を時系列に前記表示手段に表示させる。

　いくつかの実施形態の第３態様では、第２態様において、前記解析部は、前記複数の萎縮領域画像の少なくとも１つの萎縮領域画像に対して、前記複数の萎縮領域画像の１つである基準画像に対する差分を特定する差分特定処理を行い、前記表示制御部は、前記少なくとも１つの萎縮領域画像における前記特定された差分に対応する領域を識別可能に前記表示手段に表示させる。

　いくつかの実施形態の第４態様では、第３態様において、前記解析部は、直前に取得された萎縮領域画像を前記基準画像として当該萎縮領域画像の差分を特定することにより前記複数の萎縮領域画像に対して順次に前記差分特定処理を行う。

　いくつかの実施形態の第５態様では、第２態様～第４態様のいずれかにおいて、前記解析部は、前記データに基づいて前記眼底における黄斑部の位置を特定し、前記表示制御部は、前記解析部により特定された前記黄斑部の位置を表す画像を前記萎縮領域画像にオーバーレイ表示させる。

　いくつかの実施形態の第６態様では、第１態様～第５態様のいずれかにおいて、前記解析部は、前記データに基づいてＡスキャン方向の複数の層領域を特定するセグメンテーション処理部と、前記セグメンテーション処理部により特定されたブルッフ膜より強膜側の層領域のＡスキャン方向の画素値の積算値と、前記ブルッフ膜より角膜側の層領域のＡスキャン方向の画素値の積算値との比の分布情報を生成する分布情報生成部と、を含み、前記分布情報に基づいて前記萎縮領域を特定する。

　いくつかの実施形態の第７態様では、第１態様～第６態様のいずれかにおいて、前記解析部は、前記データを解析することにより前記萎縮領域の形態を表す形態情報を生成する形態情報生成部を含み、前記表示制御部は、複数の前記形態情報を時系列に前記表示手段に表示させる。

　いくつかの実施形態の第８態様では、第７態様において、前記形態情報は、前記萎縮領域の面積、外周囲長の少なくとも１つを含む。

　いくつかの実施形態の第９の態様では、第８態様において、前記表示制御部は、前記解析部により特定された複数の前記萎縮領域のそれぞれに対して、前記萎縮領域の面積、又は前記萎縮領域の外周囲長の少なくとも１つを時系列に前記表示手段に表示させる。

　いくつかの実施形態の第１０態様では、第８態様又は第９態様において、前記形態情報は、前記解析部により特定された前記１以上の萎縮領域の面積の合計値、又は前記１以上の萎縮領域の外周囲長の合計値を含む。

　いくつかの実施形態の第１１態様では、第７態様～第１０態様のいずれかにおいて、前記形態情報は、前記萎縮領域の個数を含む。

　いくつかの実施形態の第１２態様では、第７態様～第１１態様のいずれかにおいて、前記表示制御部は、複数の前記形態情報に基づいて生成された前記形態情報の時間的変化を表す画像を前記表示手段に表示させる。

　いくつかの実施形態の第１３態様は、第７態様～第１２態様のいずれかにおいて、前記解析部により特定された萎縮領域を指定するための指定部を含み、前記表示制御部は、前記指定部により指定された萎縮領域における前記形態情報の時間的変化を表す画像を前記表示手段に表示させる。

　いくつかの実施形態の第１４態様では、第１態様～第１３態様のいずれかにおいて、前記解析部は、前記データを解析することにより前記眼底における黄斑部に対する前記萎縮領域の位置又は距離を表す位置情報を生成する位置情報生成部を含み、前記表示制御部は、複数の前記位置情報に基づいて生成された前記位置情報の時間的変化を表す画像を前記表示手段に表示させる。

　いくつかの実施形態の第１５態様は、第１態様～第１２態様のいずれかにおいて、前記解析部により特定された萎縮領域を指定するための指定部を含み、前記解析部は、前記複数の１以上の萎縮領域を表す画像の少なくとも１つから前記指定部により指定された指定萎縮領域に対応する領域を特定し、前記表示制御部は、前記複数の１以上の萎縮領域を表す画像の少なくとも１つにおいて前記指定萎縮領域を識別可能に前記表示手段に表示させる。

　いくつかの実施形態の第１６態様では、第１態様～第１５態様のいずれかにおいて、前記解析部は、前記データに基づいてＡスキャン方向の１以上の層組織の層厚を特定し、特定された層厚の分布情報を生成する層厚分布情報生成部を含み、前記表示制御部は、前記層厚分布情報生成部により生成された前記層厚の分布情報を時系列に前記表示手段に表示させる。

　いくつかの実施形態の第１７態様は、光コヒーレンストモグラフィを用いて前記被検眼をスキャンすることにより前記データを収集するデータ収集部と、前記表示手段と、上記のいずれかに記載の眼科情報処理装置と、を含む眼科システムである。

　いくつかの実施形態の第１８態様は、光コヒーレンストモグラフィを用いて収集された被検眼の眼底のデータを解析することにより前記眼底における１以上の萎縮領域を特定する解析ステップと、収集タイミングが異なる複数の前記データに基づいて、前記解析ステップにおいて特定された複数の前記１以上の萎縮領域を表す画像を時系列に表示手段に表示させる表示ステップと、を含む眼科情報処理方法である。

　いくつかの実施形態の第１９態様は、第１８態様において、前記被検眼の眼底画像と前記１以上の萎縮領域を表す画像との位置合わせを行う位置合わせステップを含み、前記表示ステップは、前記位置合わせステップにおいて位置合わせが行われた前記１以上の萎縮領域を表す画像が前記眼底画像にオーバーレイされた複数の萎縮領域画像を時系列に前記表示手段に表示させる。

　いくつかの実施形態の第２０態様では、第１８態様又は第１９態様において、前記解析ステップは、前記データに基づいてＡスキャン方向の複数の層領域を特定するセグメンテーション処理ステップと、前記セグメンテーション処理ステップにおいて特定されたブルッフ膜より強膜側の層領域のＡスキャン方向の画素値の積算値と、前記ブルッフ膜より角膜側の層領域のＡスキャン方向の画素値の積算値との比の分布情報を生成する分布情報生成ステップと、を含み、前記分布情報に基づいて前記萎縮領域を特定する。

　いくつかの実施形態の第２１態様では、第１８態様～第２０態様のいずれかにおいて、前記解析ステップは、前記データを解析することにより前記萎縮領域の形態を表す形態情報を生成する形態情報生成ステップを含み、前記表示ステップは、複数の前記形態情報に基づいて生成された前記形態情報の時間的変化を表す画像を前記表示手段に表示させる。

　いくつかの実施形態の第２２態様は、コンピュータに、上記のいずれかに記載の眼科情報処理方法の各ステップを実行させるプログラムである。

　本発明に係るいくつかの実施形態によれば、眼底における萎縮領域の形態や分布を容易に経過観察するための新たな技術を提供することができるようになる。

実施形態に係る眼科システムの構成の一例を表す概略図である。実施形態に係る眼科装置の構成の一例を表す概略図である。実施形態に係る眼科情報処理装置の構成の一例を表す概略図である。実施形態に係る眼科情報処理装置の構成の一例を表す概略図である。実施形態に係る眼科情報処理装置の動作フローの一例を表す概略図である。実施形態に係る眼科情報処理装置の動作フローの一例を表す概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態に係る眼科情報処理装置の動作を説明するための概略図である。実施形態の変形例に係る眼科情報処理装置の構成の一例を表す概略図である。実施形態の変形例に係る眼科情報処理装置の動作を説明するための概略図である。実施形態の変形例に係る眼科情報処理装置の動作を説明するための概略図である。実施形態の変形例に係る眼科情報処理装置の動作フローの一例を表す概略図である。実施形態の変形例に係る眼科情報処理装置の動作を説明するための概略図である。実施形態の変形例に係る眼科装置の構成の一例を表す概略図である。

　この発明のいくつかの実施形態に係る眼科情報処理装置、眼科システム、眼科情報処理方法、及びプログラムの例について、図面を参照しながら詳細に説明する。なお、この明細書で引用する文献に記載された事項や任意の公知技術を実施形態に援用することができる。

　実施形態に係る眼科システムは、眼科情報処理装置を含む。実施形態に係る眼科情報処理方法は、眼科情報処理装置により実行される。実施形態に係る眼科情報処理方法は、プログラムにしたがってコンピュータにより実行可能である。

　実施形態に係る眼科情報処理装置は、眼科装置により光学的に取得された被検眼の眼底のデータに対して所定の解析処理や所定の表示処理を施すことが可能である。いくつかの実施形態に係る眼科装置は、被検眼の眼底の正面画像を取得する機能を備える。被検眼の眼底の正面画像を取得する機能を備えた眼科装置には、光コヒーレンストモグラフィ（Ｏｐｔｉｃａｌ　Ｃｏｈｅｒｅｎｃｅ　Ｔｏｍｏｇｒａｐｈｙ：ＯＣＴ）装置、眼底カメラ、走査型レーザー検眼鏡、スリットランプ顕微鏡、手術用顕微鏡などがある。いくつかの実施形態に係る眼科装置は、被検眼の光学的な特性を測定する機能を備える。被検眼の光学的な特性を測定する機能を備えた眼科装置には、レフラクトメーター、ケラトメーター、眼圧計、ウェーブフロントアナライザー、スペキュラーマイクロスコープ、視野計などがある。いくつかの実施形態に係る眼科装置は、レーザー治療に用いられるレーザー治療装置の機能を備える。

［眼科システム］
　図１に、実施形態に係る眼科システムの構成例のブロック図を示す。実施形態に係る眼科システム１は、眼科装置１０と、眼科情報処理装置（眼科画像処理装置、眼科解析装置）１００と、操作装置１８０と、表示装置１９０とを含む。

　眼科装置１０は、被検眼のデータを光学的に収集する。眼科装置１０は、被検眼の眼底をスキャンすることにより被検眼の眼底のデータを光学的に収集する。例えば、眼科装置１０は、ＯＣＴを用いて被検眼の眼底の３次元のＯＣＴデータを取得する。眼科装置１０は、取得された被検眼のデータから被検眼の眼底の画像を取得することが可能である。眼底の画像には、眼底の断層像及び正面画像が含まれる。眼底の断層像には、Ｂスキャン画像などがある。眼底の正面画像には、Ｃスキャン画像、シャドウグラム、又はプロジェクション画像などがある。眼科装置１０は、取得された被検眼のデータを眼科情報処理装置１００に送信する。

　いくつかの実施形態では、眼科装置１０と眼科情報処理装置１００とは、データ通信ネットワークを介して接続される。いくつかの実施形態に係る眼科情報処理装置１００は、データ通信ネットワークを介して選択的に接続された複数の眼科装置１０の１つからデータを受信する。

　眼科情報処理装置１００は、眼科装置１０によって異なる日（時間又はタイミング）に収集された被検眼の眼底の時系列の複数のデータそれぞれに対して解析処理を施す。眼科情報処理装置１００は、得られた複数の解析結果のそれぞれを被検者又は被検眼と収集日（収集日時）とに関連付けて図示しない記憶部に保存する。いくつかの実施形態に係る眼科情報処理装置１００は、得られた複数の解析結果を時系列に表示させる。いくつかの実施形態に係る眼科情報処理装置１００は、複数の解析結果から得られる統計情報を生成する。解析結果には、解析処理により得られた情報だけではなく、当該情報を加工することにより得られた情報も含まれる。

　眼科情報処理装置１００は、被検眼のデータを解析することにより地図状萎縮領域（萎縮領域）を特定し、特定された地図状萎縮領域を識別可能に表示するための画像を形成する。いくつかの実施形態に係る眼科情報処理装置１００は、眼底の正面画像や断層像における地図状萎縮領域を識別可能に表示装置１９０に表示させる。

　眼科情報処理装置１００は、取得された被検眼のデータから形成された眼底の正面画像における地図状萎縮領域に相当する領域をハイライト表示させる。いくつかの実施形態に係る眼科情報処理装置１００は、取得された被検眼のデータから眼底の正面画像を形成し、形成された正面画像と特定された地図状萎縮領域との位置合わせを行い、位置合わせが行われた地図状萎縮領域を表す画像を眼底の正面画像にオーバーレイ表示させる。

　眼科情報処理装置１００は、取得された被検眼のデータから形成された眼底の断層像における地図状萎縮領域に相当する領域をハイライト表示させる。いくつかの実施形態に係る眼科情報処理装置１００は、取得された被検眼のデータから眼底の断層像を形成し、形成された断層像と特定された地図状萎縮領域との位置合わせを行い、位置合わせが行われた地図状萎縮領域を表す画像を眼底の断層像にオーバーレイ表示させる。

　眼科情報処理装置１００は、眼底の複数のデータから特定された複数の地図状萎縮領域のそれぞれについて地図状萎縮領域の形態や分布を表すパラメータを求める。パラメータには、各地図状萎縮領域の面積や外周囲長、地図状萎縮領域の面積の合計値や外周囲長の合計値、地図状萎縮領域の個数、所定の領域内を地図状萎縮領域が占有する割合（密度）、基準位置に対する地図状萎縮領域の位置関係を表す情報などがある。いくつかの実施形態に係る眼科情報処理装置１００は、求められたパラメータを時系列に表示装置１９０に表示させる。いくつかの実施形態に係る眼科情報処理装置１００は、求められたパラメータの統計情報を表示装置１９０に表示させる。

　操作装置１８０及び表示装置１９０は、ユーザインターフェイス部として情報の表示、情報の入力、操作指示の入力など、眼科情報処理装置１００とそのユーザとの間で情報をやりとりするための機能を提供する。操作装置１８０は、レバー、ボタン、キー、ポインティングデバイス等の操作デバイスを含む。いくつかの実施形態に係る操作装置１８０は、音で情報を入力するためのマイクロフォンを含む。表示装置１９０は、フラットパネルディスプレイ等の表示デバイスを含む。いくつかの実施形態では、操作装置１８０及び表示装置１９０の機能は、タッチパネルディスプレイのような入力機能を有するデバイスと表示機能を有するデバイスとが一体化されたデバイスにより実現される。いくつかの実施形態では、操作装置１８０及び表示装置１９０は、情報の入出力を行うためのグラフィカルユーザインターフェイス（ＧＵＩ）を含む。

　眼科情報処理装置１００は、特定された複数の地図状萎縮領域のうち操作装置１８０を用いて指定された１以上の地図状萎縮領域の形態や分布を表すパラメータを時系列に表示装置１９０に表示させる。いくつかの実施形態に係る眼科情報処理装置１００は、特定された複数の地図状萎縮領域のうち操作装置１８０を用いて指定された１以上の地図状萎縮領域の形態や分布を表すパラメータの統計情報を表示装置１９０に表示させる。

［眼科装置］
　図２に、実施形態に係る眼科装置１０の構成例のブロック図を示す。

　眼科装置１０には、被検眼のＯＣＴデータを取得するための光学系が設けられている。眼科装置１０は、スウェプトソースＯＣＴを実行する機能を備えているが、実施形態はこれに限定されない。例えば、ＯＣＴの種別はスウェプトソースＯＣＴには限定されず、スペクトラルドメインＯＣＴ等であってもよい。スウェプトソースＯＣＴは、波長掃引型（波長走査型）光源からの光を測定光と参照光とに分割し、被測定物体を経由した測定光の戻り光を参照光と干渉させて干渉光を生成し、この干渉光をバランスドフォトダイオード等で検出し、波長の掃引及び測定光のスキャンに応じて収集された検出データにフーリエ変換等を施して画像を形成する手法である。スペクトラルドメインＯＣＴは、低コヒーレンス光源からの光を測定光と参照光とに分割し、被測定物体を経由した測定光の戻り光を参照光と干渉させて干渉光を生成し、この干渉光のスペクトル分布を分光器で検出し、検出されたスペクトル分布にフーリエ変換等を施して画像を形成する手法である。

　眼科装置１０は、制御部１１と、データ収集部１２と、画像形成部１３と、通信部１４とを含む。

　制御部１１は、眼科装置１０の各部を制御する。特に、制御部１１は、データ収集部１２、画像形成部１３、及び通信部１４を制御する。

　データ収集部１２は、ＯＣＴを用いて被検眼をスキャンすることにより被検眼のデータ（３次元のＯＣＴデータ）を収集する。データ収集部１２は、干渉光学系１２Ａと、スキャン光学系１２Ｂとを含む。

　干渉光学系１２Ａは、光源（波長掃引型光源）からの光を測定光と参照光とに分割し、被検眼を経由した測定光の戻り光を参照光路を経由した参照光と干渉させて干渉光を生成し、生成された干渉光を検出する。干渉光学系１２Ａは、ファイバーカプラと、バランスドフォトダイオード等の受光器とを少なくとも含む。ファイバーカプラは、光源からの光を測定光と参照光とに分割し、被検眼を経由した測定光の戻り光と参照光路を経由した参照光と干渉させて干渉光を生成する。受光器は、ファイバーカプラにより生成された干渉光を検出する。干渉光学系１２Ａは、光源を含んでよい。

　スキャン光学系１２Ｂは、制御部１１からの制御を受け、干渉光学系１２Ａにより生成された測定光を偏向することにより被検眼の眼底における測定光の入射位置を変更する。スキャン光学系１２Ｂは、例えば、被検眼の瞳孔と光学的に略共役な位置に配置された光スキャナを含む。光スキャナは、例えば、測定光を水平方向にスキャンするガルバノミラーと、垂直方向にスキャンするガルバノミラーと、これらを独立に駆動する機構とを含んで構成される。それにより、測定光を眼底平面上の任意の方向にスキャンすることができる。

　干渉光学系１２Ａによる干渉光の検出結果（検出信号）は、干渉光のスペクトルを示す干渉信号である。

　画像形成部１３は、制御部１１からの制御を受け、データ収集部１２により収集された被検眼のデータに基づいて被検眼の眼底の断層像の画像データを形成する。この処理には、ノイズ除去（ノイズ低減）、フィルタ処理、ＦＦＴ（Ｆａｓｔ　Ｆｏｕｒｉｅｒ　Ｔｒａｎｓｆｏｒｍ）などの処理が含まれている。このようにして取得される画像データは、複数のＡライン（被検眼内における各測定光の経路）における反射強度プロファイルを画像化することにより形成された一群の画像データを含むデータセットである。画質を向上させるために、同じパターンでのスキャンを複数回繰り返して収集された複数のデータセットを重ね合わせる（加算平均する）ことができる。

　画像形成部１３は、３次元のＯＣＴデータに各種のレンダリングを施すことで、Ｂスキャン画像、Ｃスキャン画像、プロジェクション画像、シャドウグラムなどを形成することができる。Ｂスキャン画像やＣスキャン画像のような任意断面の画像は、指定された断面上の画素（ピクセル、ボクセル）を３次元のＯＣＴデータから選択することにより形成される。プロジェクション画像は、３次元のＯＣＴデータを所定方向（Ｚ方向、深さ方向、Ａスキャン方向）に投影することによって形成される。シャドウグラムは、３次元のＯＣＴデータの一部（例えば特定層に相当する部分データ）を所定方向に投影することによって形成される。

　いくつかの実施形態に係る眼科装置１０には、画像形成部１３により形成された画像に対して各種のデータ処理（画像処理）や解析処理を施すデータ処理部が設けられる。例えば、データ処理部は、画像の輝度補正や分散補正等の補正処理を実行する。データ処理部は、断層像の間の画素を補間する補間処理などの公知の画像処理を実行することにより、被検眼のボリュームデータ（ボクセルデータ）を形成することができる。ボリュームデータに基づく画像を表示させる場合、データ処理部は、このボリュームデータに対してレンダリング処理を施して、特定の視線方向から見たときの擬似的な３次元画像を形成する。

　制御部１１及び画像形成部１３のそれぞれは、プロセッサを含む。プロセッサは、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ　Ｓｐｅｃｉｆｉｃ　Ｉｎｔｅｇｒａｔｅｄ　Ｃｉｒｃｕｉｔ）、プログラマブル論理デバイス（例えば、ＳＰＬＤ（Ｓｉｍｐｌｅ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｌｏｇｉｃ　Ｄｅｖｉｃｅ）、ＣＰＬＤ（Ｃｏｍｐｌｅｘ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｌｏｇｉｃ　Ｄｅｖｉｃｅ）、ＦＰＧＡ（Ｆｉｅｌｄ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｇａｔｅ　Ａｒｒａｙ））等の回路を含む。画像形成部１３の機能は、画像形成プロセッサにより実現される。いくつかの実施形態では、制御部１１及び画像形成部１３の双方の機能が１つのプロセッサにより実現される。いくつかの実施形態では、眼科装置１０にデータ処理部が設けられている場合、データ処理部の機能もまたプロセッサにより実現される。

　プロセッサは、例えば、記憶回路や記憶装置に格納されているプログラムを読み出し実行することで、実施形態に係る機能を実現する。記憶回路や記憶装置の少なくとも一部がプロセッサに含まれていてよい。また、記憶回路や記憶装置の少なくとも一部がプロセッサの外部に設けられていてよい。

　記憶装置等は、各種のデータを記憶する。記憶装置等に記憶されるデータとしては、データ収集部１２により取得されたデータ（測定データ、撮影データ等）や、被検者及び被検眼に関する情報などがある。記憶装置等には、眼科装置１０の各部を動作させるための各種のコンピュータプログラムやデータが記憶されていてよい。

　通信部１４は、制御部１１からの制御を受け、眼科情報処理装置１００との間で情報の送信又は受信を行うための通信インターフェース処理を実行する。

　いくつかの実施形態に係る眼科装置１０は、画像形成部１３により形成された被検眼の画像データを眼科情報処理装置１００に送信する。

　いくつかの実施形態に係る眼科装置１０には、被検眼の眼底の画像を取得するための眼底カメラ、走査型レーザー検眼鏡、スリットランプ顕微鏡が設けられる。いくつかの実施形態では、眼底カメラにより取得される眼底画像は、フルオレセイン蛍光眼底造影画像又は眼底自発蛍光検査画像である。

［眼科情報処理装置］
　図３及び図４に、実施形態に係る眼科情報処理装置１００の構成例のブロック図を示す。図３は、眼科情報処理装置１００の機能ブロック図を表す。図４は、図３の解析部２００の機能ブロック図を表す。

　眼科情報処理装置１００は、眼科装置１０によって異なる日に収集された被検眼の眼底の時系列の複数のデータそれぞれに対して解析処理を施す。眼科情報処理装置１００は、解析処理により得られた複数の解析結果のそれぞれを被検者又は被検眼と検査日（検査日時）とに関連付けて後述の記憶部１１２に保存する。いくつかの実施形態に係る眼科情報処理装置１００は、得られた複数の解析結果のそれぞれを被検眼と検査日とに関連付けて、所定のネットワークを介して接続された記憶装置に保存する。眼科情報処理装置１００は、解析処理により得られた複数の解析結果を時系列に表示させる。

　眼科情報処理装置１００は、眼科装置１０により取得された被検眼の眼底のデータを解析して眼底における地図状萎縮領域を特定する。眼科情報処理装置１００は、眼底の正面画像又は断層像において、特定された地図状萎縮領域を識別可能に表示装置１９０に表示させる。眼科情報処理装置１００は、検査日が異なる複数の眼底のデータに基づいて特定された複数の地図状萎縮領域を表す画像を時系列に表示装置１９０に表示させる。眼科情報処理装置１００は、特定された複数の地図状萎縮領域のそれぞれについて地図状萎縮領域の形態や分布を表すパラメータを求め、求められたパラメータを時系列に表示装置１９０に表示させる。いくつかの実施形態に係る眼科情報処理装置１００は、求められたパラメータの統計情報を表示装置１９０に表示させる。

　眼科情報処理装置１００は、制御部１１０と、画像形成部１２０と、データ処理部１３０と、通信部１４０とを含む。

　画像形成部１２０は、制御部１１０からの制御を受け、眼科装置１０により取得された３次元のＯＣＴデータからＢスキャン画像、Ｃスキャン画像、プロジェクション画像、シャドウグラムなどを形成する。画像形成部１２０は、画像形成部１３と同様に上記の画像を形成することが可能である。

　データ処理部１３０は、画像形成部１２０により形成された画像に対して各種のデータ処理（画像処理）や解析処理を施す。例えば、データ処理部１３０は、画像の輝度補正や分散補正等の補正処理を実行する。データ処理部１３０は、断層像の間の画素を補間する補間処理などの公知の画像処理を実行することにより、被検眼のボリュームデータ（ボクセルデータ）を形成することができる。ボリュームデータに基づく画像を表示させる場合、データ処理部１３０は、このボリュームデータに対してレンダリング処理を施して、特定の視線方向から見たときの擬似的な３次元画像を形成する。

　データ処理部１３０は、形成された被検眼の画像に対して所定のデータ処理を施す。データ処理部１３０は、解析部２００と、位置合わせ処理部２１０とを含む。

　解析部２００は、画像形成部１２０により形成された被検眼の眼底の画像データ（又は眼科装置１０により取得された被検眼の眼底の画像データ）に対して所定の解析処理を施す。いくつかの実施形態に係る解析処理には、眼底における地図状萎縮領域の特定処理や、地図状萎縮領域の分布情報の生成処理や、地図状萎縮領域の形態情報の生成処理や、分布情報や形態情報の統計処理や、眼底における層厚の分布情報の生成処理などが含まれる。

　解析部２００は、セグメンテーション処理部２０１と、領域特定部２０２と、分布情報生成部２０３と、形態情報生成部２０４と、層厚分布情報生成部２０５とを含む。解析部２００の各部は、被検眼の眼底の複数のデータのそれぞれに対して以下の処理を行う。

　セグメンテーション処理部２０１は、眼科装置１０により取得された被検眼のデータに基づいてＡスキャン方向の複数の層領域を特定する。いくつかの実施形態に係るセグメンテーション処理部２０１は、３次元のＯＣＴデータを解析することにより、被検眼の複数の組織に相当する複数の部分データセットを特定する。セグメンテーション処理は、特定の組織や組織境界を特定するための画像処理である。例えば、セグメンテーション処理部２０１は、ＯＣＴデータに含まれる各Ａスキャン画像における画素値（輝度値）の勾配を求め、勾配が大きい位置を組織境界として特定する。なお、Ａスキャン画像は、眼底の深さ方向にのびる１次元画像データである。なお、眼底の深さ方向は、例えば、Ｚ方向、ＯＣＴ測定光の入射方向、軸方向、干渉光学系の光軸方向などとして定義される。

　典型的な例において、セグメンテーション処理部２０１は、眼底（網膜、脈絡膜等）及び硝子体を表す３次元のＯＣＴデータを解析することにより、眼底の複数の層組織に相当する複数の部分データセットを特定する。各部分データセットは、層組織の境界によって画成される。部分データセットとして特定される層組織の例として、網膜を構成する層組織がある。網膜を構成する層組織には、内境界膜、神経繊維層、神経節細胞層、内網状層、内顆粒層、外網状層、外顆粒層、外境界膜、視細胞層、ＲＰＥがある。セグメンテーション処理部２０１は、ブルッフ膜、脈絡膜、強膜、硝子体等に相当する部分データセットを特定することができる。いくつかの実施形態に係るセグメンテーション処理部２０１は、病変部に相当する部分データセットを特定する。病変部の例として、剥離部、浮腫、出血、腫瘍、ドルーゼンなどがある。

　いくつかの実施形態に係るセグメンテーション処理部２０１は、ＲＰＥに対して強膜側の所定のピクセル数分の層組織をブルッフ膜として特定し、当該層組織に相当する部分データセットをブルッフ膜の部分データセットとして取得する。

　領域特定部２０２は、セグメンテーション処理部２０１により特定された複数の層組織の部分データセットを解析することにより、地図状萎縮領域を特定するための２つの層組織に相当する領域を特定する。いくつかの実施形態に係る領域特定部２０２は、ブルッフ膜に相当する領域より強膜側の層組織に対応する第１領域と、内境界膜に相当する領域からＲＰＥに相当する領域までの層組織に対応する第２領域とを特定する。いくつかの実施形態では、第２領域は、ブルッフ膜に相当する領域より角膜側の層組織に対応する領域である。

　分布情報生成部２０３は、各Ａスキャンについて、領域特定部２０２により特定された第１領域及び第２領域の画素値に基づいてコントラスト比を求め、眼底平面（Ａスキャン方向に直交する平面）におけるコントラスト比の２次元の分布情報を生成する。いくつかの実施形態では、分布情報生成部２０３は、各Ａスキャンについて、領域特定部２０２により特定された第１領域の画素値の積算値と第２領域の画素値の積算値との比の分布情報を生成する。いくつかの実施形態に係る分布情報生成部２０３は、各Ａスキャンについて、第１領域のＡスキャン方向の画素値の積算値に対する第２領域のＡスキャン方向の画素値の積算値の比をコントラスト比として求め、求められたコントラスト比の２次元の分布情報を生成する。以下、コントラスト比の２次元の分布情報をコントラストマップと表記する。

　解析部２００は、分布情報生成部２０３により生成されたコントラストマップにおいて、コントラスト比が大きい位置を地図状萎縮に起因して信号成分が減衰した位置として特定する。解析部２００は、特定された位置に基づいて地図状萎縮領域を特定する。例えば、解析部２００は、分布情報生成部２０３により生成されたコントラストマップにおいて、コントラスト比が所定の閾値以上である位置を含む領域を地図状萎縮領域として特定する。このような地図状萎縮領域を特定する手法に関連する技術が、米国特許出願公開第２０１５／０２０１８２９号明細書、特開２０１５－１３６６２６号公報、又は特開２０１６－１０７１４８号公報に開示されている。

　形態情報生成部２０４は、特定された地図状萎縮領域の形態を表す形態情報を生成する。形態情報の例として、地図状萎縮領域の面積や外周囲長などがある。形態情報生成部２０４は、地図状萎縮領域が描出された画像に対して公知の手法を適用することにより地図状萎縮領域の面積や外周囲長を求めることが可能である。いくつかの実施形態に係る形態情報生成部２０４は、特定された地図状萎縮領域ごとに形態情報を生成する。いくつかの実施形態に係る形態情報生成部２０４は、特定された地図状萎縮領域ごとの形態パラメータ（面積、外周囲長）の合計値を形態情報として生成する。いくつかの実施形態では、形態情報は、特定された地図状萎縮領域の個数を含む。

　層厚分布情報生成部２０５は、セグメンテーション処理部２０１により特定された複数の層組織の部分データセットを解析することにより、各層組織のＡスキャン方向の厚さを特定し、眼底平面における各層の層厚の２次元の分布情報を生成する。いくつかの実施形態に係る層厚分布情報生成部２０５は、操作装置１８０を用いて指定された１以上の層組織の層厚の２次元の分布情報（Ａスキャン方向に直交する平面の分布情報）を生成する。いくつかの実施形態に係る層厚分布情報生成部２０５は、内境界膜、神経線維層（ＮＦＬ）、神経節細胞層（ＧＣＬ）、内網状層（ＩＰＬ）、内顆粒層（ＩＮＬ）、外網状層（ＯＰＬ）、外顆粒層（ＯＮＬ）、外境界膜（ＥＬＭ）、ＲＰＥ、脈絡膜、強膜、脈絡膜－強膜境界面（ＣＳＩ）の少なくとも１つ、又は隣接する２以上の層の層厚の２次元の分布情報を生成する。

　解析部２００は、被検眼の眼底の複数のデータのそれぞれに対する上記の解析処理により得られた複数の解析処理結果に基づいてトレンド解析処理を行う。トレンド解析処理には、複数の解析処理結果の時間的変化を示す情報を生成する処理が含まれる。時間的変化を示す情報には、各解析処理結果を時系列にプロットしたグラフや、各解析処理結果を時間軸上にプロットしたグラフや、眼底平面の各位置における解析処理結果の分布を表す解析マップなどがある。いくつかの実施形態では、時間的変化を示す情報には、回帰分析により得られた回帰直線又は回帰曲線、ｐ値、未来の所定の時期おける解析結果の予測値、イベント情報、正常眼データに基づく時間的変化を示す情報などが含まれる。いくつかの実施形態では、解析処理の対象は、操作装置１８０を用いて指定される。

　位置合わせ処理部２１０は、画像形成部１２０によって形成された眼底の正面画像と解析部２００によって特定された地図状萎縮領域を表す画像との位置合わせを行う。位置合わせ処理部２１０は、画像形成部１２０によって形成された眼底の断層像と解析部２００によって特定された地図状萎縮領域を表す画像との位置合わせを行う。位置合わせ処理部２１０は、画像の拡大、縮小、回転等を行うためのアフィン変換等の公知の処理により画像の位置合わせを行うことが可能である。

　位置合わせ処理は、例えば、双方の画像から特徴部位を検出する処理と、双方の特徴部位を基準として双方の画像を位置合わせする処理とを含む。いくつかの実施形態では、位置合わせ処理は、眼底の正面画像又は断層像における地図状萎縮領域の位置情報を用いて正面画像又は断層像における地図状萎縮領域を表す画像の位置を特定する処理と、特定された地図状萎縮領域を表す画像を正面画像又は断層像に対して位置合わせする処理とを含む。

　通信部１４０は、制御部１１０からの制御を受け、眼科情報処理装置１００の通信部１４との間で情報の送信又は受信を行うための通信インターフェース処理を実行する。

　制御部１１０は、眼科情報処理装置１００の各部を制御する。特に、制御部１１０は、画像形成部１２０と、データ処理部１３０と、通信部１４０とを制御する。制御部１１０は、主制御部１１１と、記憶部１１２とを含む。主制御部１１１は、表示制御部１１１Ａを含む。

　表示制御部１１１Ａは、各種の情報を表示装置１９０に表示させる。例えば、表示制御部１１１Ａは、画像形成部１２０により形成された被検眼の眼底画像（正面画像、断層像）や、データ処理部１３０によるデータ処理結果（解析処理結果を含む）の画像を表示装置１９０に表示させる。特に、表示制御部１１１Ａは、被検眼の眼底画像を表示装置１９０に表示させ、且つ眼底画像における地図状萎縮領域に対応する領域を識別可能に表示させる。いくつかの実施形態に係る表示制御部１１１Ａは、被検眼の眼底画像を表示装置１９０に表示させ、且つ眼底画像における地図状萎縮領域に対応する領域をハイライト表示させる。例えば、表示制御部１１１Ａは、地図状萎縮領域又はその背景領域内の画素の輝度がそれ以外の領域の画素の輝度より高くなるように地図状萎縮領域又はその背景領域を表示させる。いくつかの実施形態に係る表示制御部１１１Ａは、位置合わせ処理部２１０により位置合わせが行われた地図状萎縮領域を表す画像を眼底画像にオーバーレイ表示させる。

　また、表示制御部１１１Ａは、検査日が異なる眼底の複数のデータから特定された地図状萎縮領域を表す画像を時系列に表示装置１９０に表示させることができる。いくつかの実施形態に係る表示制御部１１１Ａは、位置合わせ処理部２１０により位置合わせが行われた地図状萎縮領域を表す画像が眼底画像又は断層像にオーバーレイされた複数の萎縮領域画像を時系列に表示装置１９０に表示させる。

　また、表示制御部１１１Ａは、形態情報生成部２０４により生成された形態情報を表示装置１９０に表示させる。いくつかの実施形態に係る表示制御部１１１Ａは、形態情報生成部２０４により生成された形態情報を当該形態情報に対応する地図状萎縮領域に関連付けて表示装置１９０に表示させる。いくつかの実施形態に係る表示制御部１１１Ａは、眼底の時系列の複数のデータに基づいて形態情報生成部２０４により生成された複数の形態情報を時系列に表示装置１９０に表示させる。いくつかの実施形態に係る表示制御部１１１Ａは、操作装置１８０を用いて指定された１以上の地図状萎縮領域ごとに当該１以上の地図状萎縮領域に対応した形態情報（各地図状萎縮領域の形態情報、１以上の地図状萎縮領域の形態情報の合計値）を時系列に表示装置１９０に表示させる。

　制御部１１０は、操作装置１８０に対するユーザの操作内容に対応した操作指示信号に基づいて、眼科システム１の各部を制御する。

　制御部１１０、画像形成部１２０、及びデータ処理部１３０のそれぞれは、プロセッサを含む。画像形成部１２０の機能は、画像形成プロセッサにより実現される。データ処理部１３０の機能は、データ処理プロセッサにより実現される。いくつかの実施形態では、制御部１１０、画像形成部１２０、及びデータ処理部１３０のうち少なくとも２つの機能が１つのプロセッサにより実現される。

　記憶部１１２は、各種のデータを記憶する。記憶部１１２に記憶されるデータとしては、眼科装置１０により取得されたデータ（測定データ、撮影データ等）、画像形成部１２０により形成された画像データ、データ処理部１３０によるデータ処理結果、被検者及び被検眼に関する情報などがある。記憶部１１２には、眼科情報処理装置１００の各部を動作させるための各種のコンピュータプログラムやデータが記憶されていてよい。

　操作装置１８０は、「指定部」の一例である。表示装置１９０は、「表示手段」の一例である。地図状萎縮領域は、「萎縮領域」の一例である。

［動作例］
　いくつかの実施形態に係る眼科情報処理装置１００の動作例について説明する。

　図５に、実施形態に係る眼科情報処理装置１００の動作の一例を示す。図５は、眼科情報処理装置１００の動作例のフロー図を表す。図５では、眼科情報処理装置１００（記憶部１１２）において、眼科装置１０により取得された被検眼の３次元のＯＣＴデータが既に保存されているものとする。

（Ｓ１：被検者を選択）
　ユーザが操作装置１８０を用いて被検者ＩＤを入力することにより、被検者を選択する。

（Ｓ２：検査データを表示）
　記憶部１１２には、被検者ＩＤに対応して被検者の検査データがあらかじめ関連付けられたデータベースが記憶されている。制御部１１０は、ステップＳ１において入力された被検者ＩＤを検索キーとしてデータベースを検索し、当該被検者ＩＤに対応した検査データを取得する。表示制御部１１１Ａは、データベースを検索することにより得られた当該被検者ＩＤに対応した検査データを表示装置１９０に表示させる。検査データには、過去の検査において取得された１以上の被検眼画像が含まれる。

（Ｓ３：被検眼画像を選択）
　眼科情報処理装置１００は、ステップＳ２において表示装置１９０に表示された被検者の検査データ中の１以上の被検眼画像から解析対象の被検眼画像を選択させる。被検者は、操作装置１８０に対して操作を行うことにより、解析対象の被検眼画像を選択する。制御部１１０は、ユーザによる操作装置１８０に対する操作内容に対応した操作指示信号の入力を受け付ける。

（Ｓ４：表示）
　表示制御部１１１Ａは、ステップＳ３において入力された操作指示信号に基づいて指定された被検眼画像を選択し、選択された被検眼画像を表示装置１９０に表示させる。

（Ｓ５：領域解析実行？）
　次に、制御部１１０は、ステップＳ４において表示された被検眼画像に対して地図状萎縮領域の解析を実行するか否かを判定する。制御部１１０は、操作装置１８０に対する解析実行を指示する操作内容に対応した操作指示信号に基づいて地図状萎縮領域の解析を実行するか否かを判定することが可能である。

　地図状萎縮領域の解析を実行すると判定されたとき（Ｓ５：Ｙ）、眼科情報処理装置１００の動作はステップＳ６に移行する。地図状萎縮領域の解析を実行しないと判定されたとき（Ｓ５：Ｎ）、眼科情報処理装置１００の動作はステップＳ９に移行する。

（Ｓ６：萎縮領域を特定）
　ステップＳ５において地図状萎縮領域の解析を実行すると判定されたとき（Ｓ５：Ｙ）、制御部１１０は、解析部２００に地図状萎縮領域の解析を実行させることにより地図状萎縮領域を特定させる。ステップＳ６の詳細については後述する。制御部１１０は、眼底における地図状萎縮領域の位置や形状を特定するための領域特定情報を被検者又は被検眼と検査日とに関連付けて記憶部１１２に保存する。

（Ｓ７：形態解析）
　続いて、制御部１１０は、ステップＳ６において特定された地図状萎縮領域それぞれの面積及び外周囲長を形態情報生成部２０４に算出させる。形態情報生成部２０４は、地図状萎縮領域の面積の合計値と、外周囲長の合計値と、特定された地図状萎縮領域の個数とを含む形態情報を生成する。制御部１１０は、ステップＳ７において生成された形態情報を上記の領域特定情報と共に被検者又は被検眼と検査日とに関連付けて記憶部１１２に保存する。

　いくつかの実施形態に係る制御部１１０は、眼底における各層の層厚の２次元の分布情報を層厚分布情報生成部２０５に生成させる。制御部１１０は、ステップＳ７において生成された分布情報を上記の領域特定情報と共に被検者又は被検眼と検査日とに関連付けて記憶部１１２に保存する。

（Ｓ８：表示）
　次に、制御部１１０は、画像形成部１２０によって事前に形成された眼底の正面画像とステップＳ６において特定された地図状萎縮領域を表す画像との位置合わせを位置合わせ処理部２１０に行わせる。表示制御部１１１Ａは、当該地図状萎縮領域を表す画像を眼底の正面画像に重畳させて表示装置１９０に表示させる。ここで、眼底の正面画像は、ＲＰＥからブルッフ膜までの範囲のシャドウグラムであってよい。また、表示制御部１１１Ａは、ステップＳ７において生成された形態情報を当該形態情報に対応する地図状萎縮領域に関連付けて表示装置１９０に表示させる。

　同様に、制御部１１０は、画像形成部１２０によって事前に形成された眼底の断層像とステップＳ６において特定された地図状萎縮領域を表す画像との位置合わせを位置合わせ処理部２１０に行わせる。表示制御部１１１Ａは、当該地図状萎縮領域を表す画像を眼底の断層像に重畳させて表示装置１９０に表示させる。また、表示制御部１１１Ａは、ステップＳ７において生成された形態情報を当該形態情報に対応する地図状萎縮領域に関連付けて表示装置１９０に表示させる。以上で眼科情報処理装置１００の動作は終了である（エンド）。

（Ｓ９：トレンド解析？）
　ステップＳ５において地図状萎縮領域の解析を実行すると判定されなかったとき（Ｓ５：Ｎ）、制御部１１０は、ステップＳ４において表示された被検眼画像に対してトレンド解析を実行するか否かを判定する。制御部１１０は、操作装置１８０に対する解析実行を指示する操作内容に対応した操作指示信号に基づいてトレンド解析を実行するか否かを判定することが可能である。

　トレンド解析を実行すると判定されたとき（Ｓ９：Ｙ）、眼科情報処理装置１００の動作はステップＳ１０に移行する。トレンド解析を実行しないと判定されたとき（Ｓ９：Ｎ）、眼科情報処理装置１００の動作は終了である（エンド）。

（Ｓ１０：トレンド解析処理）
　ステップＳ９においてトレンド解析を実行すると判定されたとき（Ｓ９：Ｙ）、制御部１１０は、トレンド解析処理を解析部２００に実行させる。ステップＳ１０では、所定の被検眼のデータの解析処理結果が得られていない場合には、必要に応じてステップＳ６における地図状萎縮領域の特定及びステップＳ７における形態解析が実行されてよい。例えば、解析部２００は、地図状萎縮領域の面積や外周囲長のトレンドグラフを生成する。

（Ｓ１１：時系列表示）
　表示制御部１１１Ａは、地図状萎縮領域を表す画像が眼底画像又は断層像にオーバーレイされた複数の萎縮領域画像を時系列に表示装置１９０の所定の表示領域に表示させる。また、表示制御部１１１Ａは、層厚分布情報生成部２０５により生成された眼底平面における層厚の複数の２次元の分布情報を時系列に表示装置１９０の所定の表示領域に表示させる。

（Ｓ１２：トレンドグラフ表示）
　表示制御部１１１Ａは、ステップＳ１０において生成されたトレンドグラフを表示装置１９０の所定の表示領域に表示させる。以上で、眼科情報処理装置１００の動作は終了である（エンド）。

　次に、図６～図１２を参照しつつ図５のステップＳ６の動作例について説明する。

　図６に、図５のステップＳ６の動作例のフロー図を示す。図７は、ステップＳ２２の動作説明図を表す。図８は、ステップＳ２３の動作説明図を表す。図９Ａは、ステップＳ２４の動作説明図を表す。図９Ｂは、ステップＳ２５の動作説明図を表す。図１０は、ステップＳ２６の動作説明図を表す。図１１及び図１２は、ステップＳ２７の動作説明図を表す。

（Ｓ２１：Ｂスキャン画像を取得）
　上記のようにステップＳ５において地図状萎縮領域の解析を実行すると判定されたとき（Ｓ５：Ｙ）、制御部１１０は、記憶部１１２に保存されている被検眼の眼底のデータを読み出し、読み出されたデータに基づいてＢスキャン画像を画像形成部１２０に形成させる。いくつかの実施形態において、ステップＳ２１では、眼科装置１０からＢスキャン画像が取得される。

（Ｓ２２：セグメンテーション処理）
　制御部１１０は、ステップＳ２１において取得されたＢスキャン画像に対するセグメンテーション処理をセグメンテーション処理部２０１に実行させる。セグメンテーション処理部２０１は、ステップＳ２１において取得されたＢスキャン画像に対してＡスキャン方向の複数の層領域を特定する。セグメンテーション処理部２０１は、図７に示すように、Ｂスキャン画像ＩＭＧ１において、網膜を構成する、内境界膜３００、神経繊維層、神経節細胞層、内網状層、内顆粒層、外網状層、外顆粒層、外境界膜、視細胞層、ＲＰＥ３０１を特定する。また、セグメンテーション処理部２０１は、特定されたＲＰＥ３０１に対して強膜側の所定のピクセル数分の層組織をブルッフ膜３０２として特定する。

（Ｓ２３：コントラストマップを作成）
　続いて、制御部１１０は、ステップＳ２２におけるセグメンテーション処理の結果を用いてコントラストマップをデータ処理部１３０に作成させる。すなわち、領域特定部２０２は、セグメンテーション処理部２０１により特定された複数の層組織の部分データセットを解析することにより、ブルッフ膜３０２に相当する領域より強膜側の層組織に対応する第１領域と、内境界膜３００に相当する領域からＲＰＥ３０１に相当する領域までの層組織に対応する第２領域とを特定する。

　分布情報生成部２０３は、各Ａスキャンについて、第１領域のＡスキャン方向の画素値の積算値に対する第２領域のＡスキャン方向の画素値の積算値の比をコントラスト比として求め、求められたコントラスト比の２次元の分布情報をコントラストマップとして生成する（図８）。

（Ｓ２４：スムージング処理）
　次に、制御部１１０は、ステップＳ２３において生成されたコントラストマップに対してスムージング処理をデータ処理部１３０に実行させる。隣接する画素同士の画素値の変化が一般的に小さい傾向があり、画素値に重畳されるノイズ成分も同様であることに着目し、スムージング処理によりノイズ成分が除去されたコントラストマップが得られる（図９Ａ）。

（Ｓ２５：２値化処理）
　続いて、制御部１１０は、ステップＳ２４におけるスムージング処理後コントラストマップに対して２値化処理をデータ処理部１３０に実行させる。それにより、図９Ｂに示すような２値化マップが得られる。

（Ｓ２６：領域を探索）
　制御部１１０は、ステップＳ２５において得られた２値化マップに対して公知の領域拡張法を適用することにより領域を解析部２００に探索させる（図１０）。

（Ｓ２７：輪郭を抽出）
　制御部１１０は、ステップＳ２６における探索により得られた領域に対して公知の輪郭抽出処理を施すことにより、領域の輪郭を解析部２００に抽出させる（図１１）。解析部２００は、抽出された輪郭に基づいて地図状萎縮領域を特定する（図１２）。以上で、図５のステップＳ６の処理は終了である（エンド）。

　図１３に、いくつかの実施形態におけるステップＳ８において表示装置１９０に表示される解析情報の一例を示す。

　例えば、ステップＳ８では、表示制御部１１１Ａは、地図状萎縮領域を表す画像ＩＭＧＸをシャドウグラム（眼底の正面画像）ＩＭＧ２に重畳させて表示装置１９０に表示させる。

　また、表示制御部１１１Ａは、地図状萎縮領域の面積の合計値、外周囲長の合計値、地図状萎縮領域の数を含む形態情報３５０を表示装置１９０に表示させることができる。いくつかの実施形態に係る表示制御部１１１Ａは、各地図状萎縮領域の形態情報を対応する地図状萎縮領域に関連付けて表示装置１９０に表示させる。それにより、各地図状萎縮領域の形態を詳細に観察することができるようになる。

　図１４に、いくつかの実施形態におけるステップＳ８において表示装置１９０に表示される解析情報の他の例を示す。

　例えば、ステップＳ８では、表示制御部１１１Ａは、地図状萎縮領域を表す画像ＩＭＧＹ（Ｂスキャン断面の画像）を眼底のＢスキャン画像ＩＭＧ３に重畳させて表示装置１９０に表示させる。それにより、Ｂスキャン画像において地図状萎縮領域の形態を詳細に観察することができるようになる。

　図１５に、いくつかの実施形態におけるステップＳ１１及びステップＳ１２において表示装置１９０に表示されるトレンド解析情報の一例を示す。

　例えば、ステップＳ１１及びステップＳ１２では、表示制御部１１１Ａは、上記のように特定された地図状萎縮領域を表す画像が眼底画像（又は断層像）にオーバーレイされた複数の萎縮領域画像を時系列に表示装置１９０の第１時系列画像表示領域ＡＲ１に表示させる。また、表示制御部１１１Ａは、層厚分布情報生成部２０５により生成された眼底平面における複数の層厚マップ（層厚の２次元の分布情報）を時系列に表示装置１９０の第２時系列画像表示領域ＡＲ２に表示させる。いくつかの実施形態では、第２時系列画像表示領域ＡＲ２に表示された層厚マップは、第１時系列画像表示領域ＡＲ１に表示された萎縮領域画像に対応する表示される。例えば、同一検査日の層厚マップ及び萎縮領域画像が上下方向に並ぶように表示される。

　また、表示制御部１１１Ａは、特定された地図状萎縮領域の面積の合計値の時間的変化を表すトレンドグラフＧＰ１と、特定された地図状萎縮領域の外周囲長の合計値の時間的変化を表すトレンドグラフＧＰ２とを表示装置１９０の解析結果表示領域ＡＲ３に表示させる。

　トレンドグラフＧＰ１は、横軸に時間（例えば、被検者の年齢）を表し、縦軸に特定されたすべての地図状萎縮領域の面積の合計値を表す。トレンドグラフＧＰ１には、当該合計値の回帰直線が表示される。いくつかの実施形態では、トレンドグラフＧＰ１には、所望の時期における合計値の予測値が表示される。いくつかの実施形態では、トレンドグラフＧＰ１の横軸のスケール又はグラフの表示開始タイミング（トレンド表示される最初の検査日）は、操作装置１８０を用いてユーザにより指定される。いくつかの実施形態では、トレンドグラフＧＰ１上のプロット値が操作装置１８０を用いて指定されたとき、指定されたプロット値に対応する萎縮領域画像又は層厚マップが強調表示される。いくつかの実施形態では、トレンドグラフＧＰ２上のプロット値が操作装置１８０を用いて指定されたとき、指定されたプロット値に対応するトレンドグラフＧＰ１上のプロット位置が強調表示される。

　トレンドグラフＧＰ２は、横軸に時間（例えば、被検者の年齢）を表し、縦軸に特定されたすべての地図状萎縮領域の外周囲長の合計値を表す。トレンドグラフＧＰ２には、当該合計値の回帰直線が表示される。いくつかの実施形態では、トレンドグラフＧＰ２には、所望の時期における合計値の予測値が表示される。いくつかの実施形態では、トレンドグラフＧＰ２の横軸のスケール又はグラフの表示開始タイミング（トレンド表示される最初の検査日）は、操作装置１８０を用いてユーザにより指定される。いくつかの実施形態では、トレンドグラフＧＰ２上のプロット値が操作装置１８０を用いて指定されたとき、指定されたプロット値に対応する萎縮領域画像又は層厚マップが強調表示される。いくつかの実施形態では、トレンドグラフＧＰ１上のプロット値が操作装置１８０を用いて指定されたとき、指定されたプロット値に対応するトレンドグラフＧＰ２上のプロット位置が強調表示される。

　いくつかの実施形態では、常に同じ表示条件で表示されるように、トレンドグラフＧＰ１に対する表示条件（横軸のスケール、表示開始日）の変更内容は、トレンドグラフＧＰ２に対する表示条件に反映される。

　いくつかの実施形態に係る解析部２００は、操作装置１８０を用いてユーザにより指定された地図状萎縮領域（指定萎縮領域）を特定する。表示制御部１１１Ａは、時系列に表示された複数の萎縮領域画像の少なくとも１つにおいて、当該指定された地図状萎縮領域を識別可能に表示装置１９０に表示させる。

　いくつかの実施形態に係る表示制御部１１１Ａは、操作装置１８０を用いてユーザにより指定された地図状萎縮領域における形態情報の時間的変化を表す画像（トレンドグラフなど）を表示装置１９０に表示させる。

＜変形例＞
　いくつかの実施形態に係る構成は、上記の構成に限定されるものではない。

〔第１変形例〕
　萎縮型加齢黄斑変性の診断には、黄斑部（中心窩）に対する地図状萎縮の相対位置が有効である。いくつかの実施形態の変形例に係る眼科情報処理装置は、眼科装置により取得された被検眼の眼底のデータを解析することにより眼底における黄斑部（中心窩）に対する地図状萎縮領域の位置又は距離を表す位置情報を生成する。

　以下、いくつかの実施形態の変形例に係る眼科情報処理装置について、上記の実施形態に係る眼科情報処理装置との相違点を中心に説明する。本変形例に係る眼科情報処理装置の構成が上記の眼科情報処理装置１００の構成と異なる点は、解析部である。

　図１６に、実施形態の変形例に係る解析部２００ａの構成例のブロック図を示す。図１６において、図４と同様の部分には同一符号を付し、適宜説明を省略する。

　本変形例において、図３に示すデータ処理部１３０における解析部２００に代えて、図１６に示す変形例に係る解析部２００ａが設けられる。解析部２００ａが解析部２００と異なる点は、解析部２００に対して位置情報生成部２０６ａが追加されている点である。

　解析部２００ａは、被検眼の３次元のＯＣＴデータを公知の手法により解析することにより中心窩に相当する領域を特定し、中心窩を中心とする所定の半径の領域を黄斑部として特定する。

　位置情報生成部２０６ａは、解析部２００ａにより特定された黄斑部の代表位置に対する地図状萎縮領域の代表位置の相対位置、双方の代表位置の距離を表す位置情報、所定期間における地図状萎縮領域の代表位置の移動方向や移動距離を示すベクトル情報、又は所定期間における黄斑部の代表位置に対する地図状萎縮領域の代表位置の移動方向や移動距離を示すベクトル情報を表す位置情報を生成する。黄斑部の代表位置の例として、中心窩の位置、黄斑部の重心位置、黄斑部の輪郭部において地図状萎縮領域に最も近い（或いは最も遠い）位置などがある。地図状萎縮領域の代表位置の例として、地図状萎縮領域の中心位置、重心位置、地図状萎縮領域の輪郭部において黄斑部（又は中心窩）に最も近い（或いは最も遠い）位置などがある。いくつかの実施形態の変形例に係る表示制御部１１１Ａは、位置情報生成部２０６ａにより生成された位置情報を対応する地図状萎縮領域に関連付けて表示装置１９０に表示させる。

　いくつかの実施形態に係る表示制御部１１１Ａは、位置情報の時間的変化を表す画像（位置情報のトレンドグラフ）を表示装置１９０に表示させる。

　図１７に、実施形態の変形例において表示装置１９０に表示される解析情報の一例を示す。

　例えば、表示制御部１１１Ａは、地図状萎縮領域を表す画像ＩＭＧＸと解析部２００ａにより特定された黄斑部の位置（範囲）を表す画像ＩＭＧＺとをシャドウグラム（眼底の正面画像）ＩＭＧ２に重畳させて表示装置１９０に表示させる。画像ＩＭＧＺは、中心窩の位置を表す画像であってよい。

　また、表示制御部１１１Ａは、地図状萎縮領域の面積の合計値、外周囲長の合計値、地図状萎縮領域の数を含む形態情報に加えて、黄斑部に対する地図状萎縮領域の相対位置を表す位置情報を表示装置１９０に表示させることができる。いくつかの実施形態に係る表示制御部１１１Ａは、各地図状萎縮領域の位置情報を対応する地図状萎縮領域に関連付けて表示装置１９０に表示させる。それにより、各地図状萎縮領域の位置を詳細に観察することができるようになる。

　図１８に、実施形態の変形例において表示装置１９０に表示される解析情報の他の例を示す。

　例えば、表示制御部１１１Ａは、地図状萎縮領域を表す画像ＩＭＧＹ（Ｂスキャン断面の画像）と解析部２００ａにより特定された黄斑部の位置（範囲）を表す画像ＩＭＧＺ１を眼底のＢスキャン画像ＩＭＧ３に重畳させて表示装置１９０に表示させる。それにより、Ｂスキャン画像において黄斑部に対する地図状萎縮領域の位置を詳細に観察することができるようになる。

〔第２変形例〕
　いくつかの実施形態に係る眼科情報処理装置１００は、１以上の萎縮領域画像について基準となる萎縮領域画像に対する差分を識別可能に表示装置１９０に表示させる。

　本変形例に係る解析部は、複数の萎縮領域画像の少なくとも１つの萎縮領域画像に対して、複数の萎縮領域画像の１つである基準画像に対する差分を特定する差分特定処理を行う。本変形例に係る表示制御部は、上記の少なくとも１つの萎縮領域画像における上記の差分に対応する領域を識別可能に表示装置１９０に表示させる。いくつかの実施形態に係る変形例では、基準画像は、最先の検査日において取得された被検眼の眼底のデータを用いて生成された萎縮領域画像である。

　図１９に、本変形例に係る眼科情報処理装置の動作例のフロー図を示す。図１９において、図５と同様の部分には同一符号を付し、適宜説明を省略する。図１９においても、図５と同様に、眼科情報処理装置（記憶部１１２）において、眼科装置により取得された被検眼の３次元のＯＣＴデータが既に保存されているものとする。

　ステップＳ１～ステップＳ９は、図５と同様である。

（Ｓ３１：解析データの取得）
　ステップＳ９においてトレンド解析を実行すると判定されたとき（Ｓ９：Ｙ）、制御部１１０は、トレンド解析を行うための１以上の解析データ（萎縮領域画像を含む）を記憶部１１２から読み出す。ステップＳ３１では、ステップＳ１０と同様に、制御部１１０は、トレンド解析処理を解析部２００に実行させる。また、ステップＳ１０と同様に、所定の被検眼のデータの解析処理結果が得られていない場合には、必要に応じてステップＳ６における地図状萎縮領域の特定及びステップＳ７における形態解析が実行されてよい。例えば、解析部２００は、地図状萎縮領域の特定や地図状萎縮領域の面積や外周囲長のトレンドグラフを生成する。

（Ｓ３２：基準画像との差分を特定）
　次に、制御部１１０は、最先の検査日に取得された被検眼の眼底のデータを用いて特定された地図状萎縮領域画像を表す画像が眼底画像又は断層像にオーバーレイされた萎縮領域画像を基準画像として特定する。制御部１１０は、ステップＳ３１において取得された解析データに含まれる萎縮領域画像に対して基準画像に対する差分特定処理を解析部に実行させる。

（Ｓ３３：次の解析データ？）
　次の解析データがあるとき（Ｓ３３：Ｙ）、眼科情報処理装置１００動作はステップＳ３１に移行する。次の解析データがないとき（Ｓ３３：Ｎ）、眼科情報処理装置１００の動作はステップＳ３４に移行する。

　いくつかの実施形態に係る変形例では、ステップＳ３２において、ステップＳ３１において差分特定処理対象の萎縮領域画像の直前に取得された萎縮領域画像を基準画像とする。それにより、解析部は、直前に取得された萎縮領域画像を基準画像として当該萎縮領域画像の差分を特定することにより複数の萎縮領域画像に対して順次に上記の差分特定処理を行うことができる。

（Ｓ３４：差分領域の識別処理）
　ステップＳ３３において次の解析データがないとき（Ｓ３３：Ｎ）、制御部１１０は、ステップＳ３２において特定された差分領域の識別処理をデータ処理部等に行わせる。例えば、データ処理部は、当該差分領域が識別表示されるように当該差分領域の画素値を変更するように差分領域に対応する画像データを変更する。

（Ｓ３５：時系列表示）
　表示制御部１１１Ａは、ステップＳ３４において差分領域に対して識別処理が施された地図状萎縮領域を表す画像が眼底画像又は断層像にオーバーレイされた複数の萎縮領域画像を時系列に表示装置１９０の所定の表示領域に表示させる。また、表示制御部１１１Ａは、層厚分布情報生成部２０５により生成された眼底平面における層厚の複数の２次元の分布情報を時系列に表示装置１９０の所定の表示領域に表示させる。

（Ｓ３６：トレンドグラフ表示）
　表示制御部１１１Ａは、ステップＳ１０において生成されたトレンドグラフを表示装置１９０の所定の表示領域に表示させる。以上で、眼科情報処理装置１００の動作は終了である（エンド）。

　図２０に、いくつかの実施形態の変形例における萎縮領域画像の一例を示す。図２０において、図１３と同様の部分には同一符号を付し、適宜説明を省略する。

　例えば、表示制御部１１１Ａは、地図状萎縮領域を表す画像ＩＭＧＸ１をシャドウグラム（眼底の正面画像）ＩＭＧ２に重畳させて表示装置１９０に表示させる。地図状萎縮領域を表す画像ＩＭＧＸ１には、基準画像に対する差分領域ＩＭＧＷが識別可能に描出される。

　また、表示制御部１１１Ａは、地図状萎縮領域の面積の合計値、外周囲長の合計値、地図状萎縮領域の数を含む形態情報３５０を表示装置１９０に表示させることができる。いくつかの実施形態に係る表示制御部１１１Ａは、各地図状萎縮領域の形態情報を対応する地図状萎縮領域に関連付けて表示装置１９０に表示させる。それにより、各地図状萎縮領域の形態を詳細に観察することができるようになる。また、基準画像に対する差分領域の形状を詳細に把握することができるようになる。

〔第３変形例〕
　いくつかの実施形態に係る眼科装置は、眼科装置１０の機能に加えて、眼科情報処理装置１００の機能、操作装置１８０の機能、及び表示装置１９０の機能の少なくとも１つを備える。

　以下、いくつかの実施形態の変形例に係る眼科装置について、上記の実施形態に係る眼科装置との相違点を中心に説明する。

　図２１に、実施形態の変形例に係る眼科装置１０ｂの構成例のブロック図を示す。図２１において、図２と同様の部分には同一符号を付し、適宜説明を省略する。

　本変形例に係る眼科装置１０ｂの構成が上記の実施形態に係る眼科装置１０の構成と異なる点は、眼科装置１０ｂが、眼科情報処理装置１００の機能と、操作装置１８０の機能と、表示装置１９０の機能とを備えている点である。眼科装置１０ｂは、制御部１１ｂと、データ収集部１２と、画像形成部１３と、眼科情報処理部１５ｂと、操作部１６ｂと、表示部１７ｂとを含む。

　制御部１１ｂは、眼科装置１０ｂの各部を制御する。特に、制御部１１ｂは、データ収集部１２と、画像形成部１３と、眼科情報処理部１５ｂと、操作部１６ｂと、表示部１７ｂとを制御する。

　眼科情報処理部１５ｂは、眼科情報処理装置１００と同様の構成を有し、眼科情報処理装置１００と同様の機能を備えている。操作部１６ｂは、操作装置１８０と同様の構成を有し、操作装置１８０と同様の機能を備えている。表示部１７ｂは、表示装置１９０と同様の構成を有し、表示装置１９０と同様の機能を備えている。

　本変形例によれば、コンパクトな構成で、地図状萎縮領域の形態や分布を詳細に観察することが可能な眼科装置を提供することができる。

＜効果＞
　以下、いくつかの実施形態に係る眼科情報処理装置、眼科システム、眼科情報処理方法、及びプログラムの効果について説明する。

　いくつかの実施形態に係る眼科情報処理装置（１００）は、解析部（２００、２００ａ）と、表示制御部（１１１Ａ）とを含む。解析部は、光コヒーレンストモグラフィを用いて収集された被検眼の眼底のデータを解析することにより眼底における１以上の萎縮領域（地図状萎縮領域）を特定する。表示制御部は、収集タイミングが異なる複数のデータに基づいて、解析部により特定された複数の１以上の萎縮領域を表す画像を時系列に表示手段（表示装置１９０）に表示させる。

　このような構成によれば、収集タイミングが異なる複数の眼底のデータに基づいて特定された複数の１以上の萎縮領域を表す画像を時系列に表示手段に表示させるようにしたので、地図状萎縮領域などの萎縮領域の形態や分布を容易に経過観察することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置は、位置合わせ処理部（２１０）を含む。位置合わせ処理部は、被検眼の眼底画像と１以上の萎縮領域を表す画像との位置合わせを行う。表示制御部は、位置合わせ処理部により位置合わせが行われた１以上の萎縮領域を表す画像が眼底画像にオーバーレイされた複数の萎縮領域画像を時系列に表示手段に表示させる。

　このような構成によれば、地図状萎縮領域などの萎縮領域を表す画像が眼底画像にオーバーレイされた複数の萎縮領域画像を時系列に表示手段に表示させるようにしたので、萎縮領域の形態や分布の経過観察が容易になり、萎縮型加齢黄斑変性に対する適確な診断を補助することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置では、解析部は、複数の萎縮領域画像の少なくとも１つの萎縮領域画像に対して、複数の萎縮領域画像の１つである基準画像に対する差分を特定する差分特定処理を行う。表示制御部は、少なくとも１つの萎縮領域画像における上記の特定された差分に対応する領域を識別可能に表示手段に表示させる。

　このような構成によれば、萎縮領域画像において基準画像に対する差分を識別可能に表示手段に表示させるようにしたので、萎縮領域の形態や分布の変化を容易に把握することができるようになる。それにより、萎縮型加齢黄斑変性に対する適確な診断を補助することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置では、解析部は、直前に取得された萎縮領域画像を基準画像として当該萎縮領域画像の差分を特定することにより複数の萎縮領域画像に対して順次に差分特定処理を行う。

　このような構成によれば、複数の萎縮領域画像において、直前の萎縮領域画像に対する差分を識別可能に表示手段に表示させるようにしたので、萎縮領域の形態や分布の変化を容易に把握することができるようになる。それにより、萎縮型加齢黄斑変性に対する適確な診断を補助することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置では、解析部は、上記の眼底のデータに基づいて眼底における黄斑部の位置を特定する。表示制御部は、解析部により特定された黄斑部の位置を表す画像を萎縮領域画像にオーバーレイ表示させる。

　このような構成によれば、萎縮型加齢黄斑変性の進行度を容易に把握することができるため、進行する萎縮型加齢黄斑変性に対する適確な診断を補助することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置では、解析部は、セグメンテーション処理部（２０１）と、分布情報生成部（２０３）とを含む。セグメンテーション処理部は、上記の眼底のデータに基づいてＡスキャン方向の複数の層領域を特定する。分布情報生成部は、セグメンテーション処理部により特定されたブルッフ膜より強膜側の層領域のＡスキャン方向の画素値の積算値と、ブルッフ膜より角膜側の層領域のＡスキャン方向の画素値の積算値との比の分布情報（コントラストマップ）を生成する。解析部は、分布情報に基づいて萎縮領域を特定する。

　このような構成によれば、光コヒーレンストモグラフィを用いて取得されたデータから地図状萎縮領域などの萎縮領域を特定し、特定された萎縮領域の画像を時系列に表示手段に表示させるようにしたので、被検者の負担を軽減しつつ、地図状萎縮領域などの萎縮領域の形態や分布の経過観察が容易になる。

　いくつかの実施形態に係る眼科情報処理装置では、解析部は、形態情報生成部（２０４）を含む。形態情報生成部は、上記の眼底のデータを解析することにより萎縮領域の形態を表す形態情報を生成する。表示制御部は、複数の形態情報を時系列に表示手段に表示させる。

　このような構成によれば、地図状萎縮領域などの萎縮領域の形態や分布の経過観察が容易になり、萎縮領域の形態を詳細に観察することができるようになる。表示制御部は、形態情報を対応する萎縮領域に関連付けて表示手段に表示させてもよい。

　いくつかの実施形態に係る眼科情報処理装置では、形態情報は、萎縮領域の面積、外周囲長の少なくとも１つを含む。

　このような構成によれば、萎縮領域の面積や外周囲長の時間的変化を表示させるようにしたので、萎縮領域の形態や分布の詳細な経過観察が容易になる。

　いくつかの実施形態に係る眼科情報処理装置では、表示制御部は、解析部により特定された複数の萎縮領域のそれぞれに対して、萎縮領域の面積、又は萎縮領域の外周囲長の少なくとも１つを時系列に表示手段に表示させる。

　このような構成によれば、個々の地図状萎縮領域などの萎縮領域の形態や分布の経過観察が容易になり、萎縮領域の形態を詳細に観察することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置では、形態情報は、解析部により特定された１以上の萎縮領域の面積の合計値、又は１以上の萎縮領域の外周囲長の合計値を含む。

　このような構成によれば、萎縮領域の形態や分布の経過観察を定量的に、且つ詳細に観察することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置では、形態情報は、萎縮領域の個数を含む。

　いくつかの実施形態に係る眼科情報処理装置では、表示制御部は、複数の形態情報に基づいて生成された形態情報の時間的変化を表す画像を表示手段に表示させる。

　このような構成によれば、地図状萎縮領域などの萎縮領域の形態や分布の経過観察が容易になり、萎縮領域の形態を詳細に観察することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置は、指定部（操作装置１８０）を含む。指定部は、解析部により特定された萎縮領域を指定するために用いられる。表示制御部は、指定部により指定された萎縮領域における形態情報の時間的変化を表す画像を表示手段に表示させる。

　このような構成によれば、指定部により指定された萎縮領域の形態情報の時間的変化を容易に把握することができるので、着目する萎縮領域に対する適確な診断を補助することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置では、解析部は、位置情報生成部（２０６ａ）を含む。位置情報生成部は、上記の眼底のデータを解析することにより眼底における黄斑部に対する萎縮領域の位置又は距離を表す位置情報を生成する。表示制御部は、複数の位置情報に基づいて生成された位置情報の時間的変化を表す画像を表示手段に表示させる。

　このような構成によれば、黄斑部と萎縮領域との位置関係の変化を容易に把握することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置は、指定部（操作装置１８０）を含む。指定部は、解析部により特定された萎縮領域を指定するために用いられる。解析部は、複数の１以上の萎縮領域を表す画像の少なくとも１つから指定部により指定された指定萎縮領域に対応する領域を特定する。表示制御部は、複数の１以上の萎縮領域を表す画像の少なくとも１つにおいて指定萎縮領域を識別可能に表示手段に表示させる。

　このような構成によれば、複数の萎縮領域画像において、指定部を用いて指定された萎縮領域だけを識別表示させることができるようになる。それにより、所望の萎縮領域に着目して形態や分布の変化を把握することができるようになる。

　いくつかの実施形態に係る眼科情報処理装置は、層厚分布情報生成部（２０５）を含む。層厚分布情報生成部は、上記の眼底のデータに基づいてＡスキャン方向の１以上の層組織の層厚を特定し、特定された層厚の分布情報を生成する。表示制御部は、層厚分布情報生成部により生成された層厚の分布情報を時系列に表示手段に表示させる。

　このような構成によれば、眼底における層厚の分布の経過観察が容易になり、層厚の分布の変化を詳細に観察することができるようになる。

　いくつかの実施形態に係る眼科システムは、光コヒーレンストモグラフィを用いて被検眼をスキャンすることにより上記の眼底のデータを収集するデータ収集部（１２）と、表示手段と、上記のいずれかに記載の眼科情報処理装置と、を含む。

　このような構成によれば、収集タイミングが異なる複数の眼底のデータに基づいて特定された複数の１以上の萎縮領域を表す画像を時系列に表示手段に表示させることで、地図状萎縮領域などの萎縮領域の形態や分布を容易に経過観察することが可能な眼科システムを提供することができるようになる。

　いくつかの実施形態に係る眼科情報処理方法は、解析ステップと、表示ステップとを含む。解析ステップは、光コヒーレンストモグラフィを用いて収集された被検眼の眼底のデータを解析することにより眼底における１以上の萎縮領域を特定する。表示ステップは、収集タイミングが異なる複数の眼底のデータに基づいて、解析ステップにおいて特定された複数の１以上の萎縮領域を表す画像を時系列に表示手段（表示装置１９０）に表示させる。

　いくつかの実施形態に係る眼科情報処理方法は、位置合わせステップを含む。位置合わせステップは、被検眼の眼底画像と１以上の萎縮領域を表す画像との位置合わせを行う。表示ステップは、位置合わせステップにおいて位置合わせが行われた１以上の萎縮領域を表す画像が眼底画像にオーバーレイされた複数の萎縮領域画像を時系列に表示手段に表示させる。

　いくつかの実施形態に係る眼科情報処理方法では、解析ステップは、セグメンテーション処理ステップと、分布情報生成ステップとを含む。セグメンテーション処理ステップは、上記の眼底のデータに基づいてＡスキャン方向の複数の層領域を特定する。分布情報生成ステップは、セグメンテーション処理ステップにおいて特定されたブルッフ膜より強膜側の層領域のＡスキャン方向の画素値の積算値と、ブルッフ膜より角膜側の層領域のＡスキャン方向の画素値の積算値との比の分布情報（コントラストマップ）を生成する。解析ステップは、分布情報に基づいて萎縮領域を特定する。

　いくつかの実施形態に係る眼科情報処理方法では、解析ステップは、上記の眼底のデータを解析することにより萎縮領域の形態を表す形態情報を生成する。表示ステップは、複数の形態情報に基づいて生成された形態情報の時間的変化を表す画像を表示手段に表示させる。

　いくつかの実施形態に係るプログラムは、コンピュータに、上記のいずれかに記載の眼科情報処理方法の各ステップを実行させる。

　このような構成によれば、コンピュータに収集タイミングが異なる複数の眼底のデータに基づいて特定された複数の１以上の萎縮領域を表す画像を時系列に表示手段に表示させるようにしたので、地図状萎縮領域などの萎縮領域の形態や分布を容易に経過観察するためのプログラムを提供することができるようになる。

　いくつかの実施形態に係る眼科情報処理方法を実現するためのプログラムを、コンピュータによって読み取り可能な任意の記録媒体に記憶させることができる。記録媒体は、磁気、光、光磁気、半導体などを利用した電子媒体であってよい。典型的には、記録媒体は、磁気テープ、磁気ディスク、光ディスク、光磁気ディスク、フラッシュメモリ、ソリッドステートドライブなどである。

　また、インターネットやＬＡＮ等のネットワークを通じてコンピュータプログラムを送受信することも可能である。

　以上に説明した態様は、この発明を実施するための例に過ぎない。この発明を実施しようとする者は、この発明の要旨の範囲内における任意の変形（省略、置換、付加等）を施すことが可能である。

１　眼科システム
１０、１０ｂ　眼科装置
１１、１１ｂ、１１０　制御部
１２　データ収集部
１２Ａ　干渉光学系
１２Ｂ　スキャン光学系
１３、１２０　画像形成部
１４、１４０　通信部
１５ｂ　眼科情報処理部
１６ｂ　操作部
１７ｂ　表示部
１００　眼科情報処理装置
１１１　主制御部
１１１Ａ　表示制御部
１１２　記憶部
１３０　データ処理部
１８０　操作装置
１９０　表示装置
２００、２００ａ　解析部
２０１　セグメンテーション処理部
２０２　領域特定部
２０３　分布情報生成部
２０４　形態情報生成部
２０５　層厚分布情報生成部
２０６ａ　位置情報生成部
２１０　位置合わせ処理部

Claims

　光コヒーレンストモグラフィを用いて収集された被検眼の眼底のデータを解析することにより前記眼底における１以上の萎縮領域を特定する解析部と、
　収集タイミングが異なる複数の前記データに基づいて、前記解析部により特定された複数の前記１以上の萎縮領域を表す画像を時系列に表示手段に表示させる表示制御部と、
　を含む眼科情報処理装置。
　前記被検眼の眼底画像と前記１以上の萎縮領域を表す画像との位置合わせを行う位置合わせ処理部を含み、
　前記表示制御部は、前記位置合わせ処理部により位置合わせが行われた前記１以上の萎縮領域を表す画像が前記眼底画像にオーバーレイされた複数の萎縮領域画像を時系列に前記表示手段に表示させる
　ことを特徴とする請求項１に記載の眼科情報処理装置。
　前記解析部は、前記複数の萎縮領域画像の少なくとも１つの萎縮領域画像に対して、前記複数の萎縮領域画像の１つである基準画像に対する差分を特定する差分特定処理を行い、
　前記表示制御部は、前記少なくとも１つの萎縮領域画像における前記特定された差分に対応する領域を識別可能に前記表示手段に表示させる
　ことを特徴とする請求項２に記載の眼科情報処理装置。
　前記解析部は、直前に取得された萎縮領域画像を前記基準画像として当該萎縮領域画像の差分を特定することにより前記複数の萎縮領域画像に対して順次に前記差分特定処理を行う
　ことを特徴とする請求項３に記載の眼科情報処理装置。
　前記解析部は、前記データに基づいて前記眼底における黄斑部の位置を特定し、
　前記表示制御部は、前記解析部により特定された前記黄斑部の位置を表す画像を前記萎縮領域画像にオーバーレイ表示させる
　ことを特徴とする請求項２～請求項４のいずれか一項に記載の眼科情報処理装置。
　前記解析部は、
　前記データに基づいてＡスキャン方向の複数の層領域を特定するセグメンテーション処理部と、
　前記セグメンテーション処理部により特定されたブルッフ膜より強膜側の層領域のＡスキャン方向の画素値の積算値と、前記ブルッフ膜より角膜側の層領域のＡスキャン方向の画素値の積算値との比の分布情報を生成する分布情報生成部と、
　を含み、
　前記分布情報に基づいて前記萎縮領域を特定する
　ことを特徴とする請求項１～請求項５のいずれか一項に記載の眼科情報処理装置。
　前記解析部は、前記データを解析することにより前記萎縮領域の形態を表す形態情報を生成する形態情報生成部を含み、
　前記表示制御部は、複数の前記形態情報を時系列に前記表示手段に表示させる
　ことを特徴とする請求項１～請求項６のいずれか一項に記載の眼科情報処理装置。
　前記形態情報は、前記萎縮領域の面積、外周囲長の少なくとも１つを含む
　ことを特徴とする請求項７に記載の眼科情報処理装置。
　前記表示制御部は、前記解析部により特定された複数の前記萎縮領域のそれぞれに対して、前記萎縮領域の面積、又は前記萎縮領域の外周囲長の少なくとも１つを時系列に前記表示手段に表示させる
　ことを特徴とする請求項８に記載の眼科情報処理装置。
　前記形態情報は、前記解析部により特定された前記１以上の萎縮領域の面積の合計値、又は前記１以上の萎縮領域の外周囲長の合計値を含む
　ことを特徴とする請求項８又は請求項９に記載の眼科情報処理装置。
　前記形態情報は、前記萎縮領域の個数を含む
　ことを特徴とする請求項７～請求項１０のいずれか一項に記載の眼科情報処理装置。
　前記表示制御部は、複数の前記形態情報に基づいて生成された前記形態情報の時間的変化を表す画像を前記表示手段に表示させる
　ことを特徴とする請求項７～請求項１１のいずれか一項に記載の眼科情報処理装置。
　前記解析部により特定された萎縮領域を指定するための指定部を含み、
　前記表示制御部は、前記指定部により指定された萎縮領域における前記形態情報の時間的変化を表す画像を前記表示手段に表示させる
　ことを特徴とする請求項７～請求項１２のいずれか一項に記載の眼科情報処理装置。
　前記解析部は、前記データを解析することにより前記眼底における黄斑部に対する前記萎縮領域の位置又は距離を表す位置情報を生成する位置情報生成部を含み、
　前記表示制御部は、複数の前記位置情報に基づいて生成された前記位置情報の時間的変化を表す画像を前記表示手段に表示させる
　ことを特徴とする請求項１～請求項１３のいずれか一項に記載の眼科情報処理装置。
　前記解析部により特定された萎縮領域を指定するための指定部を含み、
　前記解析部は、前記複数の１以上の萎縮領域を表す画像の少なくとも１つから前記指定部により指定された指定萎縮領域に対応する領域を特定し、
　前記表示制御部は、前記複数の１以上の萎縮領域を表す画像の少なくとも１つにおいて前記指定萎縮領域を識別可能に前記表示手段に表示させる
　ことを特徴とする請求項１～請求項１２のいずれか一項に記載の眼科情報処理装置。
　前記解析部は、前記データに基づいてＡスキャン方向の１以上の層組織の層厚を特定し、特定された層厚の分布情報を生成する層厚分布情報生成部を含み、
　前記表示制御部は、前記層厚分布情報生成部により生成された前記層厚の分布情報を時系列に前記表示手段に表示させる
　ことを特徴とする請求項１～請求項１５のいずれか一項に記載の眼科情報処理装置。
　光コヒーレンストモグラフィを用いて前記被検眼をスキャンすることにより前記データを収集するデータ収集部と、
　前記表示手段と、
　請求項１～請求項１６のいずれか一項に記載の眼科情報処理装置と、
　を含む眼科システム。
　光コヒーレンストモグラフィを用いて収集された被検眼の眼底のデータを解析することにより前記眼底における１以上の萎縮領域を特定する解析ステップと、
　収集タイミングが異なる複数の前記データに基づいて、前記解析ステップにおいて特定された複数の前記１以上の萎縮領域を表す画像を時系列に表示手段に表示させる表示ステップと、
　を含む眼科情報処理方法。
　前記被検眼の眼底画像と前記１以上の萎縮領域を表す画像との位置合わせを行う位置合わせステップを含み、
　前記表示ステップは、前記位置合わせステップにおいて位置合わせが行われた前記１以上の萎縮領域を表す画像が前記眼底画像にオーバーレイされた複数の萎縮領域画像を時系列に前記表示手段に表示させる
　ことを特徴とする請求項１８に記載の眼科情報処理方法。
　前記解析ステップは、
　前記データに基づいてＡスキャン方向の複数の層領域を特定するセグメンテーション処理ステップと、
　前記セグメンテーション処理ステップにおいて特定されたブルッフ膜より強膜側の層領域のＡスキャン方向の画素値の積算値と、前記ブルッフ膜より角膜側の層領域のＡスキャン方向の画素値の積算値との比の分布情報を生成する分布情報生成ステップと、
　を含み、
　前記分布情報に基づいて前記萎縮領域を特定する
　ことを特徴とする請求項１８又は請求項１９に記載の眼科情報処理方法。
　前記解析ステップは、前記データを解析することにより前記萎縮領域の形態を表す形態情報を生成する形態情報生成ステップを含み、
　前記表示ステップは、複数の前記形態情報に基づいて生成された前記形態情報の時間的変化を表す画像を前記表示手段に表示させる
　ことを特徴とする請求項１８～請求項２０のいずれか一項に記載の眼科情報処理方法。
　コンピュータに、請求項１８～請求項２１のいずれか一項に記載の眼科情報処理方法の各ステップを実行させることを特徴とするプログラム。