(57)【要約】【課題】精度の高い経時比較を行うことを可能とする。【解決手段】  被検眼による収差が補正された前記被検眼からの戻り光に基づいて前記被検眼の第１眼底画像を取得する第１取得手段と、前記被検眼の第２眼底画像を取得する第２取得手段と、前記第１眼底画像と、第２眼底画像の位置合わせを行う位置合わせ手段と、前記位置合わせに基づき解析対象候補領域を設定する領域設定手段とを備えたことを特徴とする。【選択図】  図１
2014128367.tif 000002
【特許請求の範囲】
【請求項１】
被検眼による収差が補正された前記被検眼からの戻り光に基づいて前記被検眼の第１眼底画像を取得する第１取得手段と、
前記被検眼の第２眼底画像を取得する第２取得手段と、
前記第１眼底画像と、前記第２眼底画像との位置合わせを行う位置合わせ手段と、
前記位置合わせに基づき解析対象候補領域を設定する領域設定手段と
を備えたことを特徴とする眼科装置。
【請求項２】
前記解析対象候補領域が、前記位置合わせの結果少なくとも２つの眼底画像が重なる領域であることを特徴とする請求項１に記載の眼科装置。
【請求項３】
前記解析対象候補領域が、前記位置合わせの結果すべての眼底画像の重なる領域であることを特徴とする、請求項１に記載の眼科装置。
【請求項４】
前記解析対象候補領域が、前記位置合わせの結果、重なる眼底画像の枚数に応じて識別されて提示されることを特徴とする、請求項１に記載の眼科装置。
【請求項５】
前記解析対象候補領域から、血管領域を含む領域を除外して解析対象領域を選択することを特徴とする、請求項１に記載の眼科装置。
【請求項６】
前記解析対象候補領域から、選択された解析対象領域を含み、比較対象となる眼底画像が識別されて提示されることを特徴とする、請求項１に記載の眼科装置。
【請求項７】
請求項１〜６のいずれか１項に記載の機能をコンピュータに実行させるためのプログラム。
【請求項８】
請求項７に記載のプログラムを記憶した記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は眼科装置に関する。
【背景技術】
【０００２】
生活習慣病や失明原因の上位を占める疾病の早期診断を目的として、眼底部の検査が広く行われている。共焦点レーザー顕微鏡の原理を利用した眼科装置である走査型レーザー検眼鏡（ＳＬＯ：Ｓｃａｎｎｉｎｇ  Ｌａｓｅｒ  Ｏｐｈｔｈａｌｍｏｓｃｏｐｅ）は、測定光であるレーザーを眼底に対してラスタースキャンを行い、その戻り光の強度から平面画像を高分解能かつ高速に得る装置である。近年、被検眼の収差を波面センサーでリアルタイムに測定し、被検眼にて発生する測定光やその戻り光の収差を波面補正デバイスで補正する補償光学系を有する補償光学ＳＬＯ（ＡＯ−ＳＬＯ：ａｄａｐｔｉｖｅ  Ｏｐｔｉｃｓ  ＳＬＯ）が開発され、高横分解能な平面画像の取得を可能にしている。さらに、取得した網膜の平面画像を用いて網膜における視細胞を抽出し、その密度や分布の解析から疾病の診断や薬剤応答の評価が試みられている。
【０００３】
疾病の進行や薬剤応答による視細胞の変化を評価する場合には、経時的に取得したＡＯ−ＳＬＯ像を観察する必要がある。特許文献１では、経時的に取得された眼底断層像の比較をするために、深度方向における位置合わせを行う手段について記載されている。また特許文献２では、経時的に取得された断層像及び眼底像を並置して提示し、対応関係を示す表示法について記載されている。さらに特許文献３では経時的に取得された眼底像の、先に取得された画像上で検出された病変位置を、後に取得された画像上に提示する表示法について記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１２−１６６２０  眼底観察装置
【特許文献２】特開２０１０−２２０７７１  眼科観察装置
【特許文献３】特開２００８−２９７３２  眼科装置
【非特許文献】
【０００５】
【非特許文献１】“Ｒｅｔｉｎａｌ  ｍｏｔｉｏｎ  ｅｓｔｉｍａｔｉｏｎ  ｉｎ  ａｄａｐｔｉｖｅ  ｏｐｔｉｃｓ  ｓｃａｎｎｉｎｇ  ｌａｓｅｒ  ｏｐｈｔｈａｌｍｏｓｃｏｐｙ”Ｃｕｒｔｉｓ  Ｒ．  Ｖｏｇｅｌ  ｅｔ  ａｌ．  Ｏｐｔｉｃｓ  Ｅｘｐｒｅｓｓ，  ｖｏｌ．１４，  ｉｓｓｕｅ  ２，  ４８７−４９７  （２００６）
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ＡＯ−ＳＬＯ像における視細胞の経時比較は、集合として捉えた視細胞群の密度の変化や配置構造の乱れを比較する必要がある。このため比較する解析対象領域は、経時的に取得されたＡＯ−ＳＬＯ像間で、同一位置に置かれた同一形状の領域である必要がある。経時的に複数枚のＡＯ−ＳＬＯ像が取得された場合には、すべての撮影に含まれる領域から解析対象領域を選択するか、選択領域が一部含まれないＡＯ−ＳＬＯ像は比較対象からはずす処理をする必要がある。このとき、選択した解析対象領域に対して、比較対象となるＡＯ−ＳＬＯ像が把握できないという課題があった。このような経時比較の課題は先行技術では指摘されていなかった．
【課題を解決するための手段】
【０００７】
被検眼による収差が補正された前記被検眼からの戻り光に基づいて前記被検眼の第１眼底画像を取得する第１取得手段と、前記被検眼の第２眼底画像を取得する第２取得手段と、前記第１眼底画像と、第２眼底画像の位置合わせを行う位置合わせ手段と、前記位置合わせに基づき解析対象候補領域を設定する領域設定手段とを備えたことを特徴とする。
【発明の効果】
【０００８】
本発明によれば、経時的に取得されたＡＯ−ＳＬＯ像に関して、視細胞の密度や配置構造の乱れを比較する際に、比較対象となるＡＯ−ＳＬＯ像を把握しながら解析対象領域を選択することが可能となり、より精度の高い経時比較を行うことが可能となる。
【図面の簡単な説明】
【０００９】
【図１】実施例１に係る画像処理装置１０の機能構成の一例を示す図。
【図２】実施例１に係る画像処理装置１０の処理手順の一例を示すフローチャート。
【図３】固視灯の提示位置を操作する固視灯マップの一例を示す模式図。
【図４】撮影されたＷＦ−ＳＬＯ像、ＡＯ−ＳＬＯ像群の一例を示す図。
【図５】経時比較を行うために提示されるＡＯ−ＳＬＯ像群の一例を示す図。
【図６】図２の視細胞解析の一例の詳細を示すフローチャート。
【図７】経時評価ＡＯ−ＳＬＯ像群の重畳表示の一例。
【図８】重畳表示上に、重なり枚数に応じた表示を行う一例。
【図９】領域選択の一例を示す図。
【図１０】領域選択後に比較対象となる経時評価ＡＯ−ＳＬＯ像群が明示される一例。
【図１１】視細胞密度の経時変化を示すグラフの一例。
【発明を実施するための形態】
【００１０】
［実施例］
（実施例１）
本実施例では、補償光学ＳＬＯにより経時的に網膜を撮影した画像を取得し、視細胞密度等の指標の変化を提示する処理について説明する。
【００１１】
具体的には、黄斑周辺の複数箇所を撮影する検査を、数か月おきに複数回行った画像群を取得する。ここで、一か所の撮影は、撮影時間、フレームレートで指定される複数枚の画像を取得するので、これら一か所の撮影によって得られる画像群をＡＯ−ＳＬＯ像と呼ぶ場合がある。ここでＡＯ−ＳＬＯ像は撮影範囲を変更することが可能であるため、異なる解像度のＡＯ−ＳＬＯ像を取得することが可能である。さらに、一つの被検眼に対して取得される複数箇所、異なる解像度のＡＯ−ＳＬＯ像をまとめて、被検眼に対するＳｔｕｄｙと呼ぶ場合がある。
【００１２】
経時観察をするために、一つのＳｔｕｄｙに含まれる複数のＡＯ−ＳＬＯ像から解析対象とするＡＯ−ＳＬＯ像を選択すると、このＡＯ−ＳＬＯ像と同一固視位置で異なる時期に撮影されたものが並べて提示される。同一固視位置で撮影されていても、固視微動等の影響で撮影されたＡＯ−ＳＬＯ像には位置ずれが生じているため、ずれを補正するための位置合わせを行い、その結果重なる部分から解析対象領域を選択する。各画像で視細胞解析を行い、選択領域中の解析結果を比較することで経時的な評価を行うことができる。
【００１３】
＜画像処理装置の構成＞
図１は、本実施例に係る画像処理装置１０の機能構成を示したものである。
【００１４】
図１において、１００は画像取得部であり、補償光学ＳＬＯ装置により取得されたＡＯ−ＳＬＯ像を取得する。すなわち、画像取得部１００は、被検眼による収差が補正された被検眼からの戻り光に基づいて被検眼の第１眼底画像を取得する第１取得手段の一例に相当する。また、被検眼の第２眼底画像を取得する第２取得手段の一例に相当する。
【００１５】
取得したＡＯ−ＳＬＯ像は制御部１２０を通じて記憶部１３０に記憶される。１１０は入力情報取得部であり、ユーザーによる入力を取得する。１４０は画像処理部であり、視細胞解析部１４１、位置合わせ部１４２、領域選択部１４３、比較部１４４を含む。画像処理部１４０は、異なる時期に撮影された各ＡＯ−ＳＬＯ像の視細胞解析を行い、相互に位置合わせして重なり領域を求める。前記重なり領域中から選択した解析対象領域で、前記視細胞解析結果を比較することで、経時変化の評価を行う。１５０は出力部であり、算出した指標や比較結果をモニタ等に出力するほか、記憶部１３０に記憶されている処理結果を、不図示のデータベース等外部記憶に出力する。
【００１６】
＜画像処理装置の処理手順＞
次に、図２のフローチャートを参照して、本実施形態の画像処理装置１０の処理手順を説明する。
【００１７】
＜ステップＳ２１０＞
ステップＳ２１０において、画像取得部１００は、画像処理装置１０に接続される補償光学ＳＬＯにより、被検眼網膜の複数の高画角（Ｗｉｄｅ  Ｆｉｅｌｄ）ＳＬＯ像（以下、ＷＦ−ＳＬＯ像）及びＡＯ−ＳＬＯ像を取得する。
【００１８】
ここで、視細胞解析を行うために、一つの被検眼に対して行う撮影方法の一例を示す。本実施例で用いた収差補正ＳＬＯ装置は、固視灯の位置を変えることで、被検眼が異なる位置を凝視した状態で撮影することにより、網膜の異なる位置の撮影を行う。図３に固視灯の提示位置を操作する固視灯マップを示す。
【００１９】
まずは図３の固視灯マップで中心を選択した状態で固視灯を提示する。以下ではこの位置を基準位置と呼ぶ。このとき、提示される固視灯を凝視した被検眼を撮影すると、黄斑付近の撮影を行うことができる。この状態で、ＷＦ−ＳＬＯ像と、解像度の異なる複数のＡＯ−ＳＬＯ像を撮影する。ここでＷＦ−ＳＬＯ像とは画像サイズが８ｍｍｘ６ｍｍ、ピクセルサイズは５３３ｘ４００で、網膜の広い範囲を撮影することで網膜の全体イメージを取得する。なお、画像サイズおよびピクセルサイズは上記の値に限定されるものではなく、他の値としてもよい。以下についても同様である。
【００２０】
ＡＯ−ＳＬＯ像とＷＦ−ＳＬＯ画像と対応づけることにより、画角の狭いＡＯ−ＳＬＯ像が網膜全体のどの位置にあるかが提示される。ＡＯ−ＳＬＯ像としては、撮影領域のサイズが１．７ｍｍｘ１．７ｍｍ、０．８２ｍｍｘ０．８２ｍｍ、０．３４ｍｍｘ０．３４ｍｍ、ピクセルサイズはすべて共通で４００ｘ４００となる、３種類の解像度で撮影をする。なお、解像度のバリエーションは３種類以上であっても良いし、２種類であっても良い。ここで撮影領域が１．７ｍｍｘ１．７ｍｍであるＡＯ−ＳＬＯ像をＬ像、０．８２ｍｍｘ０．８２ｍｍであるＡＯ−ＳＬＯ像をＭ像、０．３４ｍｍｘ０．３４ｍｍであるＡＯ−ＳＬＯ像をＳ像とする。ＡＯ−ＳＬＯ像の撮影時間及びフレームレートは変更することが可能であり、ここではフレームレートが毎秒３２フレーム、撮影時間が２秒で、６４枚の画像より構成される。なお、フレームレート等は他の値であってもよい。
【００２１】
次に固視位置を０．５ｍｍ  Ｓｕｐｅｒｉｏｒの位置に移動し，移動した固視灯を凝視してもらう。その状態でＭ像及びＳ像を撮影する。同様に、０．５ｍｍ  Ｎａｓａｌ、０．５ｍｍ  Ｉｎｆｅｒｉｏｒ、０．５ｍｍ  Ｔｅｍｐｏｒａｌの順に固視位置を移動し、Ｍ像、Ｓ像の撮影を行う。
【００２２】
さらに、固視位置を１．０ｍｍ  Ｓｕｐｅｒｉｏｒの位置に移動し，移動した固視灯を凝視してもらう。その状態でＭ像及びＳ像を撮影する。同様に、１．０ｍｍ  Ｎａｓａｌ、１．０ｍｍ  Ｉｎｆｅｒｉｏｒ、１．０ｍｍ  Ｔｅｍｐｏｒａｌの順に固視位置を移動し、Ｍ像、Ｓ像の撮影を行う。
【００２３】
このようにして撮影された、１被検眼に対するＷＦ−ＳＬＯ像、ＡＯ−ＳＬＯ像群の模式図を図４に示す。図４では、ＷＦ−ＳＬＯ像上に、撮影されたＬ，Ｍ，Ｓ像が、固視位置の情報のみに基づいて重ねて表示されている。図４に示すように、一つのＷＦ−ＳＬＯ像と、１９のＡＯ−ＳＬＯ像（Ｌ像１、Ｍ像９、Ｓ像９）を一つのＳｔｕｄｙの基準セットとする。なお、基準セットは上記のものに限定されるものではない。
【００２４】
さらに、上記Ｓｔｕｄｙで、同一被検眼に対して異なる時期に撮影されたＳｔｕｄｙも取得する。
【００２５】
取得された被検眼のＷＦ−ＳＬＯ像、ＡＯ−ＳＬＯ像群は、制御部１２０を通じて記憶部１３０に保存される。
【００２６】
＜ステップＳ２２０＞
ステップＳ２２０において、入力情報取得部１１０は、記憶部１３０に保存されている複数の同一被検眼のＳｔｕｄｙの中から、ユーザーが選択したＡＯ−ＳＬＯ像の情報を取得する。ここでＡＯ−ＳＬＯ像の情報とは、固視位置や解像度などである。ここでは解像度がＳ像の場合とする。そして、選択されたＡＯ−ＳＬＯ像と同一被検眼、同一固視位置、同一解像度で、他のＳｔｕｄｙに含まれるＡＯ−ＳＬＯ像を記憶部１３０より取得し、制御部１２０を通じて記憶部１３０に保存するとともに、出力部１５０を通じてモニタ等に提示する。図５にモニタに提示された複数のＡＯ−ＳＬＯ像の模式図を示す。
【００２７】
ＡＯ−ＳＬＯ像は最近に撮影されたＡＯ−ＳＬＯ像（ｓｅｌｅｃｔｅｄ像）が左上に配置され、そこから右方向、さらに下段の左から右方向に新しい順に配置されている。なお、配置は上記の配置に限定さるものではなく、任意の配置に変更可能である。各ＡＯ−ＳＬＯ像には、撮影日時とｓｅｌｅｃｔｅｄ  Ｓ像の撮影日時からの差が提示される。図５に示された同一被検眼、同一固視位置、同一解像度の複数のＡＯ−ＳＬＯ像群を経時評価ＡＯ−ＳＬＯ像群と呼ぶ場合がある。
【００２８】
＜ステップＳ２３０＞
ステップＳ２３０において、視細胞解析部１４１は、ステップＳ２２０で記憶部１３０に保存された経時評価ＡＯ−ＳＬＯ像群の視細胞解析を行う。
【００２９】
図６に一つのＡＯ−ＳＬＯ像に対する視細胞解析の詳細を説明するためのフローチャートを示す。すべての経時評価ＡＯ−ＳＬＯ像群に対して以下のステップを繰り返して、視細胞解析を行う。
【００３０】
＜ステップＳ６１０＞
ステップＳ６１０において、入力情報取得部１１０は、ＡＯ−ＳＬＯ像を構成するフレームの中から、基準フレームを選択する。
【００３１】
ＡＯ−ＳＬＯ像は同一箇所を２秒間撮影した６４枚のフレームより構成される。しかし被検眼の固視微動のため、６４枚の撮影位置にはズレがあり、フレーム中にも歪みが生じている。この６４枚のフレームの中から、歪みが少なく撮影状態が良好なフレームを、ユーザーが基準フレームとして選択する。
【００３２】
ここでは基準フレームをユーザーが選択する場合を示したが、ソフトウェア的に選択する方法も考えられる。例えば輝度の平均値や分散を算出し、平均値や分散が大きなフレームを選択することが考えられる。さらには、周波数解析を行うことで、視細胞の存在を示すリング構造が明確となるフレームを選択することも考えられる。
【００３３】
こうしてＡＯ−ＳＬＯ像に対して選択された基準フレームが、制御部１２０を通じて記憶部１３０に保存される。
【００３４】
＜ステップＳ６２０＞
ステップＳ６２０において、視細胞解析部１４１は、ステップＳ６１０で取得されたＡＯ−ＳＬＯ像の基準フレームに基づき、視細胞解析の前処理を行う。前処理には複数の方法が考えられるが、ここでは周波数解析によるノイズ除去について示す。具体的には、基準フレームに周波数変換を行い、高周波を除去するフィルタを適応した後に逆変換する。ここでノイズとして除去する高周波のカットオフ値は、視細胞の大きさが最も小さい（視細胞の密度が最も高い）黄斑部付近で２μｍ程度であることが知られているため、２μｍより短い周期の振動を除去する。さらに、視細胞の密度は黄斑から離れるに従って低くなることが知られているため、撮影時の固視位置から得られる黄斑からの距離に基づきカットオフ周波数を変化させる。例えば、黄斑中心部では２μｍ、黄斑中心から１．０ｍｍの位置までは黄斑からの距離０．５ｍｍに対して１μｍずつ増加、黄斑から１．０ｍｍより離れた位置では４μｍ（固定）となるように決め、黄斑からの距離に応じたカットオフ周波数とする。
【００３５】
さらに、ノイズ除去の方法として、ＡＯ−ＳＬＯ像として取得された複数のフレームを用いて、その重ね合わせを行う方法がある。具体的には、各ＡＯ−ＳＬＯ像の６４枚のフレームを、アフィン変換などでレジストレーション処理を行った後、基準フレームに相当する領域に対して平均値処理を行う。この手法の精度は、レジストレーション処理の精度に依存する。そして上記レジストレーション処理を行った後に、上記周波数変換による高周波成分の除去処理を行う。
【００３６】
こうして取得された画像を前処理後画像と呼ぶ場合がある。
【００３７】
＜ステップＳ６３０＞
ステップＳ６３０において、視細胞解析部１４１は、ステップＳ６２０で取得された前処理後画像に対して、視細胞検出を行う。
【００３８】
具体的には、前処理後画像の輝度の極大値検出を行う方法がある。このとき、極大値として検出された点同士の距離が、知識として知られる視細胞の大きさよりも小さい場合には、ノイズの影響であるとして検出点同士を統合する処理を行うことで、検出のロバスト性を向上させることができる。このときに用いる視細胞の大きさとして、ステップＳ６２０と同様、解析対象となるＡＯ−ＳＬＯ像の黄斑からの距離に基づき算出することで、より精度の高い検出を行うことができる。
【００３９】
さらに、こうして取得された検出点のうち、指定されたある閾値以上の値をもつ点を視細胞とする。閾値としては、画像の最低輝度値を用いる方法（すべての検出点を視細胞とする）や、画像の輝度平均値を用いる方法などが考えられる。
【００４０】
ここでは視細胞検出の一例を示したが、検出法は上記に限定されるものではなく、様々な方法が考えられる。例えば、輝度値がある閾値以上のピクセルのみ選択し、該選択ピクセルが連結している領域の重心を求める方法がある。他には、小領域の特徴量を算出して、パターン認識の手法を用いて検出することも考えられる。具体的には、例えば１１ｘ１１ピクセルの小領域に対して、ガボール特徴量を算出する。視細胞とされる検出点を中心に含む複数の小領域と、含まない複数の小領域から得られるガボール特徴量ベクトルを用いて、サポートベクターマシンによる学習を行う。新たに対象とする小領域に対してガボール特徴量を算出し、上記学習結果を用いることで、該小領域の中心が視細胞であるか否かを判断する。
【００４１】
さらに、上記ソフトウェアによる検出に加えて、ユーザーによりマニュアルで修正することも可能である。このとき入力情報取得部１１１は、ステップＳ６３０で検出した検出点に対してユーザーが修正した検出点の位置を取得する。
【００４２】
＜ステップＳ６４０＞
ステップＳ６４０において、視細胞解析部１４１は、ステップＳ６３０で検出された検出点に対して、ボロノイ解析を行う。
【００４３】
具体的には、ステップＳ６３０で検出された全検出点に対して、近傍にある検出点同士の垂直二等分線によって画像内の領域を分割することで、各検出点に属するボロノイ領域を算出する。
【００４４】
こうして視細胞解析を行った結果を、制御部１２０を通じて記憶部１３０に保存した後、ステップＳ２３０に戻る。
【００４５】
＜ステップＳ２４０＞
ステップＳ２４０において、位置合わせ部１４２は、ステップＳ２２０で記憶部１３０に保存された経時評価ＡＯ−ＳＬＯ像群間の位置合わせを行う。すなわち、位置合わせ部１４２は第１眼底画像と、第２眼底画像との位置合わせを行う位置合わせ手段の一例に相当する。ここで位置合わせは、ステップＳ６１０で各ＡＯ−ＳＬＯ像に対して選択された基準フレームを用いて行う。
【００４６】
具体的には、経時評価ＡＯ−ＳＬＯ像群の選択の基準となったＳｅｌｅｃｔｅｄ像と、その他の各経時評価ＡＯ−ＳＬＯ像群間の移動量を、位相限定相関法により算出する。なお、移動量の算出は位相限定相関法以外の方法を用いることとしてもよい。
【００４７】
比較したいＡＯ−ＳＬＯ像の基準フレームを、それぞれＦＦＴにより周波数変換する。ここで周波数変換された画像をそれぞれＦ（ｕ，ｖ），Ｇ（ｕ，ｖ）と表すと、以下の式で示される位相限定相関係数Ｃ（ｕ，ｖ）を計算する。
【００４８】
【数１】
2014128367.tif 000003
【００４９】
Ｃ（ｕ，ｖ）の逆ＦＦＴを求め、ピーク値を検出することにより移動量を算出することができる。そして算出した移動量を、制御部１２０を通じて記憶部１３０に保存する。
図７に、算出した移動量に基づき経時評価ＡＯ−ＳＬＯ像群を重ねて表示した例を示す。図７で示した表示を、経時評価ＡＯ−ＳＬＯ像群の重畳表示と呼ぶ。図７ではＳｅｌｅｃｔｅｄ像の枠のみ太線で示している。
【００５０】
ここではソフトウェアにより位置合わせをする方法の一例について述べたが、マニュアルにより位置を合わせることも可能である。具体的には、血管等の特徴を参考にマウスのクリック＆ドラッグにより各ＡＯ−ＳＬＯ像の位置を移動する方法である。マニュアルのみで位置合わせをすることも可能であり、また上記ソフトウェアでの位置合わせの後、マニュアルで修正することも可能である。さらに、上記ではステップＳ６１０で選択された基準フレームを用いて位置合わせを行ったが、ＡＯ−ＳＬＯ像を構成するフレームより生成した合成フレームを用いてもよい。合成フレームの生成法は、非特許文献１に一例が示されている。さらに、ステップＳ６２０で取得された前処理後画像を用いて位置合わせを行ってもよい。
【００５１】
＜ステップＳ２５０＞
ステップＳ２５０において、領域選択部１４３は、ステップＳ２４０で求めた経時評価ＡＯ−ＳＬＯ像群間の位置合わせ結果に基づき、解析対象となる領域の選択を行う。すなわち、領域選択部１４３は位置合わせに基づき解析対象候補領域を設定する領域設定手段の一例に相当する。
【００５２】
ここで解析対象領域とは、ステップＳ２３０で取得された視細胞解析結果を経時的に比較し、疾病の進行による変化を評価するための指標を算出する領域となる。このため，すべての経時評価ＡＯ−ＳＬＯ像群間で比較するためには、すべての経時評価ＡＯ−ＳＬＯ像群に含まれる領域であることが求められる。ここで指標の例としては、ステップＳ６３０で求められた検出点の数、検出点数を領域の面積で割った密度、ステップＳ６４０で求められた最近傍検出点までの距離、ボロノイ領域が六角形である割合などがあげられる。
【００５３】
図８（ａ）に、図７で求めた位置合わせ結果に基づき、すべての経時評価ＡＯ−ＳＬＯ像群が重なりをもつ部分を太枠で示した図を示す。すなわち、太枠は位置合わせの結果すべての眼底画像の重なる領域を示している。なお、太枠は、すべての経時評価ＡＯ−ＳＬＯ像群を用いて比較をする場合の、選択可能な最大領域であるため、解析対象候補領域の最小領域となる。図８（ｂ）には、比較のために少なくとも２つの経時評価ＡＯ−ＳＬＯ像が重なりをもつ領域を太枠で示した図を示す。すなわち、図８（ｂ）は位置合わせの結果少なくとも２つの眼底画像が重なる領域を解析対象候補領域としたものであるこの領域は、解析対象候補領域の最大領域となる。さらに、図８（ｂ）上で重なりをもつ経時評価ＡＯ−ＳＬＯ像の数に応じて、領域を異なる色で示すなど明示することも可能である。こうすることで、選択したい特徴物が存在する場合に、その特徴物を含む解析領域を選択する場合に、いくつの経時評価ＡＯ−ＳＬＯ像群間での比較が可能かを知ることができる。
【００５４】
図８（ａ），（ｂ）に示したように、位置合わせの結果、重なる眼底画像の枚数に応じて解析対象候補領域が識別されて提示されるようになっている。
【００５５】
上記太枠で示した選択領域内から、ユーザーが解析対象領域を選択することができる。図９にユーザーによる選択領域の例を示す。図９（ａ）は、経時評価ＡＯ−ＳＬＯ像群の重畳表示上に選択領域を示した場合であり、図９（ｂ）は個別のＡＯ−ＳＬＯ像上に示した場合である。なお、図９（ａ）と図９（ｂ）表示は表示制御部により表示部に同時に表示させることとしてもよいし、ユーザーからの指示に基づいて切り替えて表示させることとしてもよい。図９（ａ）では、すべての経時評価ＡＯ−ＳＬＯ像群が重なりを持つ部分に領域を選択している。領域は矩形や多角形、楕円等を選択することが可能であり、ユーザーによるマウスクリックで設定できる。血管等が含まれる場合、血管の下の輝度は低くなるため、視細胞の検出精度が悪くなる。よって血管は避けて領域を選択する。すなわち、血管領域を含む領域を除外して解析対象領域を選択する。なお、血管領域を自動認識してユーザーにより血管領域を選択出来ないようにしてもよい。図９（ｂ）では、位置合わせの精度が足りない場合に備えて、各画像に選択領域の位置をユーザーが修正することができる。このとき選択領域の大きさや形は変化しない状態で、マウスによるクリック＆ドラッグにより移動することができる。なお、図９（ｂ）を表示部に表示させるとともに、図４に示す基準セットを並べて表示することとしてもよい。また、図４のうちＷＦ−ＳＬＯ像を抜いたＡＯ−ＳＬＯ像を図９（ｂ）および図４と並べて表示することとしてもよい。
【００５６】
図８（ｂ）で、重なりをもつ経時評価ＡＯ−ＳＬＯ像の数が明示される表示について説明したが、図１０ではさらに、図１０（ａ）上でマウスがポイントする位置に重なりをもつ経時評価ＡＯ−ＳＬＯ像が、図１０（ｂ）上に示される。そして選択した解析対象領域を含むＡＯ−ＳＬＯ像が図１０（ｂ）上に明示される。すなわち、解析対象候補領域から、選択された解析対象領域を含み、比較対象となる眼底画像が識別されて提示される。これにより、比較対象となるＡＯ−ＳＬＯ像の数や撮影時期を認識しながら、解析対象領域を選択することが可能になる。
【００５７】
こうして領域選択を行った結果を、制御部１２０を通じて記憶部１３０に保存する。
【００５８】
＜ステップＳ２６０＞
ステップＳ２６０において、比較部１４４は、ステップＳ２２０で取得された経時評価ＡＯ−ＳＬＯ像群間の比較を、ステップＳ２３０〜ステップＳ２５０で行った解析結果に基づき行う。
【００５９】
図１１に、図９の選択領域で計測した視細胞密度の経時変化を示すグラフの例を示す。ここで比較対象となるので、ステップＳ２３０で取得された視細胞解析の結果から得られる指標であり、具体的には、視細胞数、視細胞数を領域の面積で割った密度、最近傍視細胞までの距離、ボロノイ領域の面積や六角形率などがあげられるが、これらに限定されるものではない。
【００６０】
＜ステップＳ２７０＞
ステップＳ２７０において、制御部１２０は、算出した指標や比較結果など記憶部１３０に記憶されている処理結果を、データベースに保存する。
【００６１】
以上の構成により、補償光学ＳＬＯ装置により経時的に取得した複数のＡＯ−ＳＬＯ像を解析して経時変化を評価する場合に、眼底の同一領域を対象に評価を行うことが可能になる。
【００６２】
以上の構成により、補償光学ＳＬＯ装置により経時的に取得した複数の経時評価ＡＯ−ＳＬＯ像群を解析して比較を行う場合に、比較したいＡＯ−ＳＬＯ像の数や撮影時期を考慮しながら解析対象領域を選択することが可能になる。
【００６３】
（その他の実施形態）
本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記憶した記憶媒体を、システムあるいは装置に供給する。そして、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読出し実行することによっても、達成されることは言うまでもない。
【符号の説明】
【００６４】
１０  画像処理装置
１００  画像取得部
１１０  入力情報取得部
１２０  制御部
１３０  記憶部
１４０  画像処理部
１４１  視細胞解析部
１４２  位置合わせ部
１４３  領域選択部
１４４  比較部
１５０  出力部
【図１】
2014128367.tif 000004
【図２】
2014128367.tif 000005
【図３】
2014128367.tif 000006
【図４】
2014128367.tif 000007
【図５】
2014128367.tif 000008
【図６】
2014128367.tif 000009
【図７】
2014128367.tif 000010
【図８】
2014128367.tif 000011
【図９】
2014128367.tif 000012
【図１０】
2014128367.tif 000013
【図１１】
2014128367.tif 000014