(57)【要約】【課題】  精度良く視細胞を検出する。【解決手段】  被検眼の眼底画像を取得する取得手段と、眼底画像の位置に対する輝度値の変化を示す特性において、隣接する輝度値に対して凸となる輝度値を有する部分を視細胞として検出する検出手段とを備える。【選択図】  図１５
2014198223.tif 000002
【特許請求の範囲】
【請求項１】
被検眼の眼底画像を取得する取得手段と、
前記眼底画像の位置に対する輝度値の変化を示す特性において、隣接する輝度値に対して凸となる輝度値を有する部分を視細胞として検出する検出手段と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記検出手段は、隣接する輝度値の平均輝度値よりも高い輝度値を有する部分を前記凸となる輝度値を有する部分とすることを特徴とする請求項１記載の画像処理装置。
【請求項３】
前記検出手段は、所定の範囲内において凸となる輝度値を有する部分が複数存在する場合には、前記複数の凸となる輝度値を有する部分の重心を視細胞として検出することを特徴とする請求項１または請求項２に記載の画像処理装置。
【請求項４】
前記眼底画像は前記被検眼による収差が補正された画像であることを特徴とする請求項１乃至３のいずれか１項に記載の画像処理装置。
【請求項５】
被検眼の眼底画像を取得する取得工程と、
前記眼底画像の位置に対する輝度値の変化を示す特性において、隣接する輝度値に対して凸となる輝度値を有する部分を視細胞として検出する検出工程と、
を有することを特徴とする画像処理方法。
【請求項６】
請求項５記載の画像処理方法の各工程をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本件は、画像処理装置及び画像処理方法に関する。
【背景技術】
【０００２】
生活習慣病や失明原因の上位を占める疾病の早期診断を目的として、眼底部の検査が広く行われている。共焦点レーザー顕微鏡の原理を利用した眼科装置である走査型レーザー検眼鏡（ＳＬＯ：Ｓｃａｎｎｉｎｇ  Ｌａｓｅｒ  Ｏｐｈｔｈａｌｍｏｓｃｏｐｅ）は、測定光であるレーザーを眼底に対してラスタースキャンを行い、その戻り光の強度から平面画像を高分解能かつ高速に得る装置である。近年、被検眼の収差を波面センサーでリアルタイムに測定し、被検眼にて発生する測定光やその戻り光の収差を波面補正デバイスで補正する補償光学系を有する補償光学ＳＬＯ（ＡＯ−ＳＬＯ：ａｄａｐｔｉｖｅ  Ｏｐｔｉｃｓ  ＳＬＯ）が開発され、高横分解能な平面画像（以下ＡＯ−ＳＬＯ像という場合がある）の取得を可能にしている。さらに、取得した網膜の平面画像を用いて網膜における視細胞を抽出し、その密度や分布の解析から疾病の診断や薬剤応答の評価が試みられている。
【０００３】
疾病の進行や薬剤応答による視細胞の変化を評価する場合には、高解像度で取得されたＡＯ−ＳＬＯ像から視細胞を検出する必要がある。視細胞を検出する方法として、輝度の極大値部分を視細胞として検出する方法がしられている（非特許文献１参照）。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】“Ａｕｔｏｍａｔｅｄ  ｉｄｅｎｔｉｆｉｃａｔｉｏｎ  ｏｆ  ｃｏｎｅ  ｐｈｏｔｏｒｅｃｅｐｔｏｒｓ  ｉｎ  ａｄａｐｔｉｖｅ  ｏｐｔｉｃｓ  ｒｅｔｉｎａｌ  ｉｍａｇｅｓ，  Ｋａｃｃｉｅ  Ｙ．Ｌｉ  ａｎｄ  Ａｕｓｔｉｎ  Ｒｏｏｒｄａ，Ｊ．Ｏｐｔ．Ｓｏｃ．Ａｍ．Ａ，  Ｍａｙ  ２００７，  Ｖｏｌ．２４，Ｎｏ．５，１３５８”
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、本来視細胞として検出されるべき位置に極大値が現れない場合があり、極大値を用いた視細胞検出では視細胞の検出漏れが生じてしまう。。
本件の目的の一つは精度良く視細胞検出を行うことを可能とすることである。
なお、前記目的に限らず、後述する発明を実施するための形態に示す各構成により導かれる作用効果であって、従来の技術によっては得られない作用効果を奏することも本件の他の目的の１つとして位置付けることができる。
【課題を解決するための手段】
【０００６】
本画像処理装置は、被検眼の眼底画像を取得する取得手段と、前記眼底画像の位置に対する輝度値の変化を示す特性において、隣接する輝度値に対して凸となる輝度値を有する部分を視細胞として検出する検出手段と、を備える。
【発明の効果】
【０００７】
本件によれば、精度良く視細胞検出を行うことが可能となる。
【図面の簡単な説明】
【０００８】
【図１】実施例１に係る画像処理装置１０の機能構成の一例を示す図である。
【図２】実施例１に係る画像処理装置１０の処理手順の一例を示すフローチャートである。
【図３】固視灯の提示位置を操作する固視灯マップの一例を示す図である。
【図４】撮影されたＷＦ−ＳＬＯ像、ＡＯ−ＳＬＯ像群の一例を示す図である。
【図５】図２の位置合わせの一例を示すフローチャートである。
【図６】位置合わせ後のＷＦ−ＳＬＯ像、ＡＯ−ＳＬＯ像群の重畳画像に一例を示す図である。
【図７】図２の視細胞解析の一例を示すフローチャートである。
【図８】領域選択の一例を示す図である。
【図９】正常眼の分布と検出結果との関係の一例を示す図である。
【図１０】正常眼分布からのずれの一例を示す図である。
【図１１】正常眼分布からのずれの一例を示す画像を重畳画像上に表示する例である。
【図１２】実施例２に係る画像処理装置１０の処理手順の一例を示すフローチャートである。
【図１３】重畳画像に示される病変部の一例を示す図である。
【図１４】正常眼分布からのずれを示す画像を重畳画像上に表示する一例を示す図である。
【図１５】凸面検出と極大値検出との比較の一例を示す図である。
【図１６】黄斑からの距離に対する感度の特性の一例を示す図である。
【発明を実施するための形態】
【０００９】
以下に、本件の実施形態の一例を図面を参照しながら詳細に説明する。
【００１０】
［実施例］
（実施例１）
本実施例では、補償光学ＳＬＯ装置により網膜を撮影した画像を取得する際に、視細胞の密度を評価したい領域と黄斑中心（中心窩）が含まれる領域とを撮影し、密度を評価したい領域の黄斑中心からの距離を取得する。さらに、視細胞の検出を黄斑中心からの距離を反映したアルゴリズムで行い、視細胞の密度を黄斑中心からの距離の関数として取得する処理について説明する。なお、本発明において黄斑中心とは黄斑の完全な中心のみを示す概念ではなく、略中心を含む概念である。
【００１１】
具体的には、固視の位置を変化させながら黄斑中心位置、黄斑中心からの距離が０．５ｍｍ，  １．０ｍｍとなる位置の撮影を、異なる解像度の画像を組み合わせて行う。なお、黄斑中心からの距離については例示であり、これらの値に限定されるものではない。
【００１２】
収差補正ＳＬＯ画像は、高解像度であるが撮影範囲が狭いために、１つの被検眼に対して複数箇所の撮影を行う。一か所の撮影は、撮影時間、フレームレートで指定される複数枚の画像を取得する。以下、一か所の撮影によって得られる画像群をＡＯ−ＳＬＯ像と呼ぶ場合がある。ここでＡＯ−ＳＬＯ像は撮影範囲を変更することが可能であるため、異なる解像度のＡＯ−ＳＬＯ像を取得することが可能である。さらに、一つの被検眼に対して取得される複数箇所、異なる解像度のＡＯ−ＳＬＯ像をまとめて、被検眼に対するＡＯ−ＳＬＯ像群と呼ぶ場合がある。
【００１３】
ここでは被検眼に対するＡＯ−ＳＬＯ像同士の相互の位置合わせをすることで、相互の位置関係を確定させる。その後、黄斑中心が撮影されているＡＯ−ＳＬＯ像より黄斑中心を検出し、眼軸長を考慮した上で、各ＡＯ−ＳＬＯ像の黄斑中心からの距離を求める。十分な解像度をもつＡＯ−ＳＬＯ像で視細胞の検出を行い、視細胞密度等の指標を取得することで、視細胞密度を黄斑中心からの距離の関数として取得することができる。
なお、眼軸長を考慮する方法としては、例えば”  Ｂｅｎｎｅｔｔ  ＡＧ，  Ｒｕｄｎｉｃｋａ  ＡＲ，  Ｅｄｇａｒ  ＤＦ．  Ｉｍｐｒｏｖｅｍｅｎｔｓ  ｏｎ  Ｌｉｔｔｍａｎｎ’ｓ  ｍｅｔｈｏｄ  ｏｆ  ｄｅｔｅｒｍｉｎｉｎｇ  ｔｈｅ  ｓｉｚｅ  ｏｆ  ｒｅｔｉｎａｌ  ｆｅａｔｕｒｅｓ  ｂｙ  ｆｕｎｄｕｓ  ｐｈｏｔｏｇｒａｐｈｙ．  Ｇｒａｅｆｅｓ  Ａｒｃｈ  Ｃｌｉｎ  Ｅｘｐ  Ｏｐｈｔｈａｌｍｏｌ．  １９９４  Ｊｕｎ；２３２（６）：３６１−７．”
に記載されている方法を用いることが可能である。具体的には黄斑中心から眼底上の所定の位置までの距離をＬ、スキャナの振り角をθ、眼軸長をｘとして、Ｌ  ＝  ｑθ、ｑ  ＝０．０１３０６（ｘ  −  １．８２）により求めることが可能である。また、眼軸長を考慮する方法は上記の方法に限定されるものではなく他の方法を用いることとしてもよい。例えば、”  Ｌｉ  ＫＹ，  Ｔｉｒｕｖｅｅｄｈｕｌａ  Ｐ，  Ｒｏｏｒｄａ  Ａ．  Ｉｎｔｅｒｓｕｂｊｅｃｔ  Ｖａｒｉａｂｉｌｉｔｙ  ｏｆ  Ｆｏｖｅａｌ  ＣｏｎｅＰｈｏｔｏｒｅｃｅｐｔｏｒ  Ｄｅｎｓｉｔｙ  ｉｎ  Ｒｅｌａｔｉｏｎ  ｔｏ  Ｅｙｅ  Ｌｅｎｇｔｈ  Ｉｎｖｅｓｔ  Ｏｐｈｔｈａｌｍｏｌ  Ｖｉｓ  Ｓｃｉ．  ２０１０Ｄｅｃ；５１（１２）：６８５８−６７．”に記載されている方法を用いることとしてもよい。
こうして黄斑中心からの距離を算出して用いることにより、より精度の高い視細胞検出が行えると同時、黄斑中心からの距離によって変化することが知られている視細胞密度を、黄斑中心からの距離と関連づけて取得することが可能となる。
【００１４】
＜画像生成装置の構成＞
図１は、本実施例に係る画像処理装置１０の機能構成の一例を示したものである。
【００１５】
図１において、１００は画像取得部であり、補償光学ＳＬＯ装置により取得されたＡＯ−ＳＬＯ像を取得する。すなわち、画像取得部１００は、被検眼の眼底画像を取得する取得手段の一例に相当する。取得したＡＯ−ＳＬＯ像は制御部１２０を通じて記憶部１３０に記憶される。なお、画像取得部１００は補償光学ＳＬＯ装置から送信されたＡＯ−ＳＬＯ像を単に受信することでＡＯ−ＳＬＯ像を取得してもよいし、画像取得部１００が画像処理装置１０の外部に備えられたメモリ等に保存されたＡＯ−ＳＬＯ像を自ら取りにいくことで取得してもよい。すなわち、本件における「取得」とは少なくとも単に受信する場合および補償光学ＳＬＯ装置に積極的に取りにいく場合を含む概念である。
【００１６】
１１０は入力情報取得部であり、ユーザーによる入力を取得する。１４０は画像処理部であり、位置合わせ部１４１、距離算出部１４２、視細胞解析部１４３、比較部１４４を含む。
【００１７】
画像処理部１４０は、取得したＡＯ−ＳＬＯ像を相互に位置合わせして各ＡＯ−ＳＬＯ像の相対位置を求め，黄斑中心を検出するＡＯ−ＳＬＯ像と、解析対象となるＡＯ−ＳＬＯ像の位置関係から、解析対象領域の黄斑中心からの距離を取得する。画像処理部１４０は、取得した黄斑中心からの距離を反映して視細胞解析を行い、密度等の指標を算出する。
【００１８】
また、画像処理部１４０は算出した指標と記憶部１３０に保存された正常眼データとの比較を行い、比較結果が明示されるようなグラフや画像を生成する。
【００１９】
１５０は出力部であり、算出した指標や比較結果をモニタ等の表示部に出力するほか、記憶部１３０に記憶されている処理結果をデータベースに出力する。
【００２０】
＜画像処理装置の処理手順＞
次に、図２のフローチャートを参照して、本実施形態の画像処理装置１０の処理手順の一例を説明する。
【００２１】
＜ステップＳ２１０＞
ステップＳ２１０において、画像取得部１００は、画像処理装置１０に接続される補償光学ＳＬＯ装置により、例えば被検眼網膜の複数のＡＯ−ＳＬＯ像を取得する。
【００２２】
ここで、視細胞解析を行うために一つの被検眼に対して行う撮影方法の一例を示す。本実施例で用いた収差補正ＳＬＯ装置は、該装置が有する固視灯の位置を変えることで、被検眼が異なる位置を凝視した状態で撮影することにより、網膜の異なる位置の撮影を行う。図３に固視灯の提示位置を操作する固視灯マップの一例を示す。なお、この固視灯マップは例えば表示部に表示され、検査者は表示された固視灯マップを参照しながらマウス等のポインティングデバイス（指示手段）を用いて固視灯の点灯位置を所望位置に変更することが可能である。
【００２３】
動作の一例として、まず、図３の固視灯マップで中心を選択した状態で固視灯を提示する。以下ではこの位置を基準位置と呼ぶ場合がある。このとき、提示される固視灯を凝視した被検眼を撮影すると、黄斑中心付近の撮影を行うことができる。
【００２４】
この状態で、ＷＦ−ＳＬＯ像と、解像度の異なる複数のＡＯ−ＳＬＯ像を補償光学ＳＬＯ装置により撮影する。なお、撮影されたＷＦ−ＳＬＯ像は画像取得部１００により取得される。
【００２５】
ここでＷＦ−ＳＬＯ像は例えば画像サイズが８ｍｍｘ６ｍｍ、ピクセルサイズは５３３ｘ４００である。網膜の広い範囲を撮影することで網膜の全体イメージを取得する。ＷＦ−ＳＬＯ像と画角の狭いＡＯ−ＳＬＯ像とを対応づけることにより、ＡＯ−ＳＬＯ像が網膜全体のどの位置にあるかが提示される。
【００２６】
ＡＯ−ＳＬＯ像としては、撮影領域のサイズが１．７ｍｍｘ１．７ｍｍ、０．８２ｍｍｘ０．８２ｍｍ、０．３４ｍｍｘ０．３４ｍｍ、ピクセルサイズはすべて共通で４００ｘ４００の、３種類の解像度で撮影される。これらの撮影は補償光学ＳＬＯ装置により行われ撮影された画像は画像取得部１００により取得される。ここで撮影領域が１．７ｍｍｘ１．７ｍｍであるＡＯ−ＳＬＯ像をＬ像、０．８２ｍｍｘ０．８２ｍｍであるＡＯ−ＳＬＯ像をＭ像、０．３４ｍｍｘ０．３４ｍｍであるＡＯ−ＳＬＯ像をＳ像という場合がある。なお、上記の撮影領域のサイズおよびピクセルサイズは例示であり、これらの値に限定されるものではない。すなわち、他の値としてもよい。なお、Ｍ像またはＳ像は第１眼底画像の一例に相当する。さらに、ＷＦ−ＳＬＯ像またはＬ像は第２眼底画像の一例に相当し、Ｌ像またはＭ像は第３眼底画像の一例に相当する。
【００２７】
ＡＯ−ＳＬＯ像の撮影時間及びフレームレートは変更することが可能であるが、ここでは一例としてフレームレートが毎秒３２フレーム、撮影時間が２秒で、６４枚の画像より構成される。
【００２８】
次に固視位置を０．５ｍｍ  Ｓｕｐｅｒｉｏｒの位置に移動し，移動した固視灯を凝視してもらう。その状態でＭ像及びＳ像を撮影する。同様に、０．５ｍｍ  Ｎａｓａｌ、０．５ｍｍ  Ｉｎｆｅｒｉｏｒ、０．５ｍｍ  Ｔｅｍｐｏｒａｌの順に固視位置を移動し、Ｍ像、Ｓ像の撮影を行う。なお、固視位置の移動順は上記の例に限定されるものではなく、他の順序としてもよい。以下についても同様である。
【００２９】
さらに、固視位置を１．０ｍｍ  Ｓｕｐｅｒｉｏｒの位置に移動し，移動した固視灯を凝視してもらう。その状態でＭ像及びＳ像を撮影する。同様に、１．０ｍｍ  Ｎａｓａｌ、１．０ｍｍ  Ｉｎｆｅｒｉｏｒ、１．０ｍｍ  Ｔｅｍｐｏｒａｌの順に固視位置を移動し、Ｍ像、Ｓ像の撮影を行う。
【００３０】
このようにして撮影された、１の被検眼に対するＷＦ−ＳＬＯ像、ＡＯ−ＳＬＯ像群の模式図を図４に示す。図４では、ＷＦ−ＳＬＯ像上に、撮影されたＬ，Ｍ，Ｓ像が、固視位置の情報のみに基づいて重ねて表示されている。
【００３１】
図４に示されるように、０．５ｍｍの位置と１．０ｍｍの位置で撮影されたＳ像同士は重ならないが、Ｍ像同士は重なっていることがわかる。視細胞をより良く解像するためにはＳ像の解像度が好ましいが、Ｓ像のみで１．０ｍｍの領域まで重なりをもって撮影するためには、一例として少なくとも０．３ｍｍおきに撮影を行う必要があり、被験者の負荷増大につながる。そのためここでは、解析する場所を０．５ｍｍと１．０ｍｍとしその位置でＳ像の撮影を行うが、固視位置０．５ｍｍと１．０ｍｍでＭ像を撮影することで、Ｍ像同士の少なくとも一部が重なるようにする。この重なりによりＭ像相互の位置が正確に把握できるようなるためＳ像の正確な位置を把握することが可能となる。
【００３２】
ここで、非常に太い血管などの特徴がＳ像内に存在する場合には、Ｍ像および／またはＬ像の撮影を行わなくても、Ｓ像を直接ＷＦ−ＳＬＯ像に対応づけることやＭ像またはＬ像の一方を用いてＳ像をＷＦ−ＳＬＯ像に対応づけることも可能である。しかしＷＦ−ＳＬＯはＡＯ−ＳＬＯ像に比べて解像度が粗いために、正確な対応付けができないことが多い。Ｍ像やＬ像のように、収差補正ＳＬＯ画像だがＳ像よりは解像度の低い画像を組み合わせて撮影することで、被験者の負荷を減らしつつも、精度の高い位置合わせをすることが可能になる。
【００３３】
なお、取得された被検眼のＡＯ−ＳＬＯ像群は、制御部１２０を通じて記憶部１３０に保存される。
【００３４】
＜ステップＳ２２０＞
ステップＳ２２０において、入力情報取得部１１０は、記憶部１３０に保存されているＡＯ−ＳＬＯ像を構成するフレームの中から、基準フレームを選択する。
【００３５】
ステップＳ２１０で説明したように、一例としてＡＯ−ＳＬＯ像は同一箇所を２秒間撮影した６４枚のフレームより構成される。しかし被検眼の固視微動のため、６４枚の撮影位置にはズレがあり、フレーム中にも歪みが生じている。この６４枚のフレームの中から、歪みが少なく撮影状態が良好なフレームを、ユーザーが基準フレームとして選択する。
【００３６】
ここでは基準フレームをユーザーが選択する場合を示したが、ソフトウェア的に選択する方法も考えられる。例えば輝度の平均値や分散を算出し、平均値や分散が大きなフレームを選択することが考えられる。さらには、周波数解析を行うことで、視細胞の存在を示すリング構造が明確となるフレームを選択することも考えられる。
【００３７】
こうして各ＡＯ−ＳＬＯ像に対して選択された基準フレームが、制御部１２０を通じて記憶部１３０に保存される。
【００３８】
なお、ステップＳ２１０において、各固視位置で一枚のＡＯ−ＳＬＯ像を撮影するようにした場合にはステップＳ２２０は実行されない。
【００３９】
＜ステップＳ２３０＞
位置合わせ工程の一例であるステップＳ２３０において、位置合わせ部１４１は、記憶部１３０に保存されている補償光学ＳＬＯにより取得された複数のＡＯ−ＳＬＯ像の相互の位置合わせを行う。ここで位置合わせは、各ＡＯ−ＳＬＯ像に対して選択された基準フレームを用いて行う。
【００４０】
位置合わせは解像度が近い画像同士の方が高精度となるため、例えばまずはＷＦ−ＳＬＯ像とＬ像との位置合わせを行い、次にＬ像とＭ像との位置合わせ、その後にＭ像とＳ像との位置合わせを行う。位置合わせにはさまざまな方法が考えられるが、ここでは一例として位相限定相関法を用いた場合について、同一固視位置に対して撮影されたＭ像とＳ像との位置合わせの詳細を、図５のフローチャートに基づき説明する。
ここでは解像度の異なる画像間の位置合わせを行っているので、解像度が低く画角が大きい画像を大画像、解像度が高く画角が小さい画像を小画像と呼ぶ。すなわち、Ｍ像とＳ像との位置合わせの場合にはＭ像が大画像でＳ像が小画像になるが、Ｌ像とＭ像との位置合わせの場合にはＬ像が大画像でＭ像が小画像になる。
【００４１】
＜ステップＳ５１０＞
ステップＳ５１０において、位置合わせ部１４１は、ステップＳ２１０で取得された複数のＡＯ−ＳＬＯ像の中から，位置合わせを行う２種類のＡＯ−ＳＬＯ像のペアを設定する。設定の方法は様々な方法が考えられるが、例えば、まずはＷＦ−ＳＬＯ像と基準位置のＬ像とを選択し第一のペアとする。次に、上記Ｌ像と基準位置のＭ像とを選択する。さらに、固視位置０．５ｍｍの４か所（Ｓｕｐｅｒｉｏｒ，Ｎａｓａｌ，Ｉｎｆｅｒｉｏｒ，Ｔｅｍｐｏｒａｌ）の位置のＭ像のそれぞれと上記Ｌ像とを選択する。
【００４２】
その次に、Ｓｕｐｅｒｉｏｒの０．５ｍｍのＭ像と１．０ｍｍのＭ像とを選択する。同様に、Ｎａｓａｌ，Ｉｎｆｅｒｉｏｒ，Ｔｅｍｐｏｒａｌの０．５ｍｍ、１．０ｍｍのＭ画像を選択する。最後に、同一固視位置のＭ像とＳ像とを選択する。
【００４３】
このようにして設定したＡＯ−ＳＬＯ像のペアを、制御部１２０を通じて記憶部１３０に保存する。なお、上記の手順は一例であり、他の順序で画像の選択を行うこととしてもよい。
【００４４】
＜ステップＳ５２０＞
ステップＳ５２０において、位置合わせ部１４１は、ステップＳ５１０で選択された画像のペアを、制御部１２０を通じて記憶部１３０から取得する。以下では、取得した２種類の画像の解像度が異なる場合、例えば大画像がＭ像、小画像がＳ像であり、かつ同一固視位置の場合を例として詳細に説明する。
【００４５】
Ｍ像は画像サイズが約８２０μｍであり、Ｓ像は３４０μｍであり、それぞれ４００ｘ４００ピクセルサイズであることを考えると、Ｍ像内でＳ像と重なるのは中心部分の約１６６ピクセルの領域となる。位相限定相関法では各画像のＦＦＴを計算することから、高速化を考えてＭ像の中心１２８ｘ１２８ピクセルの領域を切り出す。なお、切り出す領域は上記の値に限定されるものではなく、大画像と小画像とが重なり合うサイズであればよい。また、Ｌ像とＭ像とをペアとして選択している際には、例えば固視位置に基づいて位置合わせされたＭ像の中心をＬ像から切り出す領域の中心とする。
【００４６】
＜ステップＳ５３０＞
ステップＳ５３０において、位置合わせ部１４１は、ステップＳ５２０で取得された画像ペアのうち小画像に対して、ステップＳ５２０で切り出された大画像の領域に対応する領域の解像度変換を行う。具体的には、Ｓ画像の中心から８２０ｘ１２８／４００＝２６２．４μｍの領域を用いて、１２８ｘ１２８ピクセルの画像を生成する。すなわち、位置合わせ部１４１は、位置合わせ手段による位置合わせの際に、解像度が高い画像を解像度が低い画像と同じ解像度に変換する解像度変換手段の一例に相当する。なお、ステップＳ５２０、５３０の実行順序は逆であってもよいし、同時に実行されることとしてもよい。      
【００４７】
＜ステップＳ５４０＞
ステップＳ５４０において、位置合わせ部１４１は、ステップＳ５２０で切り出された大画像（クロップト大画像）と、ステップＳ５３０で変換された小画像（低解像小画像）に対して、例えば位相限定相関法により移動量を算出する。具体的には、クロップト大画像および低解像小画像をそれぞれＦＦＴにより周波数変換する。ここで周波数変換された画像をそれぞれＦ（ｕ，ｖ），Ｇ（ｕ，ｖ）と表すと、以下の式で示される位相限定相関係数Ｃ（ｕ，ｖ）を計算する。すなわち、位置合わせ手段の一例である位置合わせ部１４１は解像度変換後の画像を用いて位置合わせを行う。
【００４８】
【数１】
2014198223.tif 000003
【００４９】
Ｃ（ｕ，ｖ）の逆ＦＦＴを求め、ピーク値を検出することにより移動量を算出することができる。
【００５０】
算出した移動量を、制御部１２０を通じて記憶部１３０に保存する。
【００５１】
＜ステップＳ５５０＞
ステップＳ５５０において、位置合わせ部１４１は、ステップＳ５１０で設定された画像ペア群のすべてに対して、移動量の算出がされたかを判断する。まだ移動量の算出が行われていないペアがあれば、該ペアに対してステップＳ５２０に戻って移動量の算出を行う。すべてのペアに関して処理が終了していれば、ステップＳ５６０に進む。
【００５２】
＜ステップＳ５６０＞
ステップＳ５６０において、位置合わせ部１４１は、ステップＳ５１０で設定されたすべての画像ペア群に対して、ステップＳ５４０で算出された移動量に基づき、画像相互の位置合わせを行う。
【００５３】
具体的には、基準位置で撮影されたＷＦ−ＳＬＯ像とＬ像との移動量（ずれ量）、Ｌ像とＭ像との移動量、Ｍ像とＳ像との移動量にもとづき、Ｌ、Ｍ、Ｓ像のＷＦ−ＳＬＯ像に対する移動量を求める。同様に、Ｌ像と固視位置０．５ｍｍのＭ像との移動量、固視位置０．５ｍｍのＭ像とＳ像との移動量に基づき、固視位置０．５ｍｍにおけるＭ、Ｓ像のＷＦ−ＳＬＯ像に対する移動量を求める。
【００５４】
同様にしてすべてのＡＯ−ＳＬＯ像の相互の位置を求め、対応する位置に例えば撮影領域の大きいものが下になるように重畳した重畳画像を生成する。すなわち、算出した移動量に基づいて位置合わせ部１４１は画像の位置合わせを行う。図６に、図４のＡＯ−ＳＬＯ像群を相互に位置合わせした、重畳画像の例を示す。なお図６に示す重畳画像は出力部１５０によりモニタ等の表示部に表示される。すなわち、出力部１５０は、位置合わせ手段により位置合わせが行われた第１画像と第２画像とを表示部に重畳表示させる。図４では固視位置のみに基づいてＡＯ−ＳＬＯ像群を配置していたために、血管等の位置がずれていたものが、図６では解像度の異なる画像を用いて段階的に位置合わせを行っているため正しい位置にＡＯ−ＳＬＯ像が配置されている。すなわち、位置合わせ部１４１は、Ｍ像およびＬ像を介してＷＦ−ＳＬＯ像とＳ像との位置合わせを行う。この位置合わせ部１４１は被検眼による収差が補正された画像である第１眼底画像と第１眼底画像よりも画角が広くかつ解像度が低い第２眼底画像とを第２眼底画像よりも画角が狭く解像度が高い第３眼底画像を用いて位置合わせを行う位置合わせ手段の一例に相当する。より具体的には、位置合わせ手段の一例である位置合わせ部１４１は、被検眼の眼底上の異なる位置を撮像した複数の第１眼底画像と第２眼底画像とを眼底上の異なる位置を撮像した複数の第３眼底画像を用いて位置合わせを行う。また、異なる観点から見れば、位置合わせ手段の一例である位置合わせ部１４１は、第２眼底画像と第３眼底画像との位置合わせ結果と、第３眼底画像と第１眼底画像との位置合わせ結果と、に基づいて第１眼底画像と前記第２眼底画像との位置合わせを行う。
【００５５】
ここではソフトウェアにより重畳画像を生成する方法の一例について述べたが、マニュアルにより位置を合わせることも可能である。具体的には、固視位置に基づき配置された図４の画像上で、血管等の特徴を参考にマウスのクリック＆ドラッグにより各ＡＯ−ＳＬＯ像の位置を移動する方法である。マニュアルのみで位置合わせをすることも可能であり、
また上記ソフトウェアでの位置合わせの後、マニュアルで修正することも可能である。このとき、例えば同一固視位置を撮影したＭ像とＳ像のペアを先に位置合わせし、そのセットをＷＦ−ＳＬＯ像及びＬ像と位置合わせするなどの方法も可能である。
【００５６】
こうして取得されたＡＯ−ＳＬＯ像群の重畳画像は、制御部１２０を通じて記憶部１３０に保存される。さらに出力部１５０を通じてモニタ等に表示した後、ステップＳ２３０に戻る。
【００５７】
＜ステップＳ２４０＞
ステップＳ２４０において、入力情報取得部１１１は、ステップＳ２２０で生成された重畳画像から、黄斑中心の位置を取得する。入力情報取得部１１１は、検出した黄斑中心位置を、制御部１２０を通じて記憶部１３０に保存する。ここで、入力情報取得部１１１は黄班を検出する検出手段の一例に相当する。
【００５８】
黄斑中心を検出する方法は複数の方法が考えられる。例えば、被験者に中心固視を凝視してもらって撮影した固視中心の画像中心を黄斑中心とするという方法がある。または、固視中心に対して撮影された画像の輝度の低い部分の中心を、ユーザーが目視により検出する方法が考えられる。これは黄斑中心部の輝度が低くなる傾向があるという知識を利用したものであるため、それをソフトウェア的に実装することも可能である。具体的には、各ピクセルの輝度の、画像内での最高値からの差を求め、その重心を求めて黄斑中心とするなどの方法がある。
【００５９】
すなわち、検出手段の一例である入力情報取得部１１１は、位置合わせ手段により位置合わせされた複数の第１眼底画像のうち黄班を含む第１眼底画像から黄班を検出する。
【００６０】
さらに、ここではＡＯ−ＳＬＯ像、ＷＦ−ＳＬＯ像のみを用いているが、同一被検眼を撮影したＯＣＴ等の他モダリティ画像が取得可能な場合には、それらの画像と比較することで、ユーザーが黄斑中心と考える位置を選択する方法などもある。
【００６１】
また、Ｓ画像は他の画像に比べて最も解像度が高いため黄斑中心を正確に検出することが可能となるが、Ｍ像、Ｌ像またはＷＦ−ＳＬＯ像からでも黄斑中心は検出可能である。
【００６２】
＜ステップＳ２５０＞
距離取得工程の一例であるステップＳ２５０において、距離算出部１４２は、ステップＳ２３０で算出された各ＡＯ−ＳＬＯ像のＷＦ−ＳＬＯ像からの移動量と、ステップＳ２４０で取得された黄斑中心位置から、各ＡＯ−ＳＬＯ像の黄斑中心からの距離を算出する。すなわち、距離算出部１４２は被検眼の黄斑（例えば黄斑中心）から位置合わせ手段により位置合わせされた第１眼底画像における所定の位置までの距離を取得する距離取得手段の一例に相当する。より具体的には距離算出部１４２は検出手段により検出された黄斑（例えば黄斑中心）から位置合わせ手段により位置合わせされた複数の第１眼底画像のうち黄班（例えば黄斑中心）を含まない第１眼底画像における所定の位置までの距離を取得する。
【００６３】
各ＡＯ−ＳＬＯ像の黄斑中心からの距離は、黄斑中心を原点として重畳画像の上側をＹ軸方向、右側をＸ軸方向として、各ＡＯ−ＳＬＯ像の例えば中心の座標を求めるにより求められる。このようにして取得される各ＡＯ−ＳＬＯ像の座標を黄斑中心座標と呼ぶ。これはステップＳ２２０で求めた各ＡＯ−ＳＬＯ像ｉの移動量をベクトルＤｉとし、ステップＳ２３０で取得された黄斑中心位置を、ＷＦ−ＳＬＯ像の中心からのベクトルＭとした場合に、Ｄｉ−Ｍで表される。ＡＯ−ＳＬＯ像ｉの中心の座標を（Ｘｉ，Ｙｉ）とすれば、その黄斑中心からの距離Ｒｉは以下のようにして求めることができる。
【００６４】
【数２】
2014198223.tif 000004
【００６５】
そして求められた各ＡＯ−ＳＬＯ像の黄斑中心座標と黄斑中心からの距離を、制御部１２０を通じて記憶部１３０に保存する。
【００６６】
＜ステップＳ２６０＞
ステップＳ２６０において、視細胞解析部１４３は、ステップＳ２１０で取得されたＡＯ−ＳＬＯ像のＳ像に対して、視細胞解析を行う。
【００６７】
図７に視細胞解析の詳細を説明するためのフローチャートを示す。
【００６８】
＜ステップＳ７１０＞
ステップＳ７１０において、視細胞解析部１４３は、ステップＳ２２０で取得されたＡＯ−ＳＬＯ像の基準フレームに基づき、視細胞解析の前処理を行う。前処理には複数の方法が考えられるが、ここでは周波数解析によるノイズ除去について示す。具体的には、基準フレームに周波数変換を行い、高周波を除去するフィルタを適応した後に逆変換する。
【００６９】
視細胞の大きさが最も小さい（視細胞の密度が最も高い）黄斑中心部付近で２μｍ程度であることが知られているため、ここでノイズとして除去する高周波のカットオフ値（カットオフ周波数）は、例えば２μｍより短い周期の振動を除去するように定められる。さらに、視細胞の密度は黄斑中心から離れるに従って低くなることが知られているため、ステップＳ２５０で取得した距離に基づきカットオフ周波数を変化させる。例えば、黄斑中心部では２μｍ、黄斑中心から１．０ｍｍの位置までは黄斑中心からの距離０．５ｍｍに対して１μｍずつ増加、黄斑中心から１．０ｍｍより離れた位置では４μｍ（固定）となるように決め、黄斑中心からの距離に応じたカットオフ周波数とする。すなわち、検出手段の一例である視細胞解析部１４３は距離取得手段により取得された距離に応じて、視細胞を検出するために用いるパラメタを変更する。
【００７０】
さらに、ノイズ除去の方法として、ＡＯ−ＳＬＯ像として取得された複数のフレームを用いて、その重ね合わせを行う方法がある。具体的には、各ＡＯ−ＳＬＯ像の６４枚のフレームを、アフィン変換などでレジストレーション処理を行った後、基準フレームに相当する領域に対して平均値処理を行う。この手法の精度は、レジストレーション処理の精度に依存する。そして上記レジストレーション処理を行った後に、上記周波数変換による高周波成分の除去処理を行う。
【００７１】
こうして取得された画像を前処理後画像と呼ぶ場合がある。
【００７２】
＜ステップＳ７２０＞
ステップＳ７２０において、視細胞解析部１４３は、ステップＳ７１０で取得された前処理後画像に対して、視細胞検出を行う。視細胞解析部１４３は第１眼底画像から視細胞を検出する検出手段の一例に相当する。
具体的には、前処理後画像の輝度の凸面検出を行う。凸検出とは、連続変数の場合には、２次微分係数が負になる領域の検出となるが、画像のように離散値の場合には、以下の式に従うピクセルを検出することに相当する。
【００７３】
【数３】
2014198223.tif 000005
【００７４】
すなわち、視細胞解析部１４３は、眼底画像の位置に対する輝度値の変化を示す特性において、隣接する輝度値に対して凸となる輝度値を有する部分を視細胞として検出する検出手段の一例に相当する。また、上記の式は、隣接する輝度値の平均輝度値よりも高い輝度値を有する部分を凸となる輝度値を有する部分とすることを示している。
【００７５】
ここで、画像は２次元なので上記凸検出が，縦横斜めすべてで成り立つ場合に凸面と判断して検出する。具体的には，下記を満たすピクセル（ｘ，  ｙ）を検出する。
【００７６】
【数４】
2014198223.tif 000006
【００７７】
このようにして凸点を検出すると，一つの視細胞（一つの視細胞し得る所定の範囲内）に対して複数のピクセルが凸点として検出される。そこで凸面として検出されたピクセルで隣接するもの同士を統合し、例えばその重心を求めることで検出点とする。すなわち、所定の範囲内において凸となる輝度値を有する部分が複数存在する場合には、複数の凸となる輝度値を有する部分の重心を視細胞として検出する。
【００７８】
このとき統合領域のピクセル数がある閾値以下の場合には検出点とみなさないという制限を設けることで，検出精度を向上することができる。具体的には，統合領域の大きさが２ピクセル以上の場合に検出点とするという処理を入れることで、ノイズにより単一ピクセルのみが周辺より大きな値をもつような場合を、検出点から除去することが可能となる．さらに，検出点の輝度値がある閾値以下の場合には検出点とみなさないという閾値処理を加えることもできる。このような処理を加えることで，血管の下部など極端に輝度が低い領域でのノイズ検出を抑制することができる。具体的には，画像全体の輝度平均値と輝度最低値の中間の値を閾値とするなどの処理が可能である。なお、上記の式では隣接するピクセルの値を用いているが精度が落ちても良い場合等には２つ隣りのピクセルの値を用いることとしても良い。すなわち、本発明は隣接する画素を用いた凸面検出に限定されるものではない。
【００７９】
つぎに，検出された点同士の距離が、知識として知られる視細胞の大きさよりも小さい場合には、ノイズの影響であるとして検出点同士を統合する処理を行う。具体的には、検出点のピクセルとその周囲のピクセル間の距離が、視細胞の大きさの半分（視細胞が円形であると考えその半径）よりも小さい場合には、その周囲のピクセルも検出点とする処理を行う。その後、検出点で隣接するもの同士を統合し、その重心を求めることで新たな検出点とする。このような処理を行うことにより、凸検出のみでは視細胞の大きさよりも小さい距離に検出された２つの検出点を統合することができ、ノイズにより一つの視細胞に複数のピークが現れるような場合の検出精度を向上させることができる。
【００８０】
ここでは視細胞検出の一例を示したが、検出法は上記に限定されるものではなく、様々な方法が考えられる。例えば、前処理後画像の輝度の極大値検出を行う。このとき極大値として検出された点同士の距離が、知識として知られる視細胞の視細胞の大きさよりも小さい場合には、ノイズの影響であるとして検出点同士を統合する処理を行う。ここで用いる視細胞の大きさとして、ステップＳ７１０と同様、ステップＳ２５０で求められた解析対象となるＡＯ−ＳＬＯ像の黄斑中心からの距離に基づき算出することで、より精度の高い検出を行うことができる。
【００８１】
さらに、こうして取得された検出点のうち、指定されたある閾値以上の値をもつ点を視細胞とする。閾値としては、画像の最低輝度値を用いる方法（すべての検出点を視細胞とする）や、画像の輝度平均値を用いる方法などが考えられる。
【００８２】
なお、上記の極大値検出と言われる方法の詳細は“Ａｕｔｏｍａｔｅｄ  ｉｄｅｎｔｉｆｉｃａｔｉｏｎ  ｏｆ  ｃｏｎｅ  ｐｈｏｔｏｒｅｃｅｐｔｏｒｓ  ｉｎ  ａｄａｐｔｉｖｅ  ｏｐｔｉｃｓ  ｒｅｔｉｎａｌ  ｉｍａｇｅｓ，  Ｋａｃｃｉｅ  Ｙ．Ｌｉ  ａｎｄ  Ａｕｓｔｉｎ  Ｒｏｏｒｄａ，Ｊ．Ｏｐｔ．Ｓｏｃ．Ａｍ．Ａ，  Ｍａｙ  ２００７，  Ｖｏｌ．２４，Ｎｏ．５，１３５８”に記載されている。
【００８３】
この方法によって、解像されて個別の輝点として存在する視細胞は大部分を検出することができる。しかし被検眼の状態や収差補正のかかり具合，前処理の精度等により、隣接する視細胞同士が明確に分離した個別の輝点とならない場合がある。特に視細胞の輝度の大きさは時間的にも空間的にもばらつきがあるため、隣接する視細胞間に大きな輝度の差があるような場合、輝度が小さい方の視細胞は極大値をとらず、上記極大値検出に基づく方法では検出できなくなる。そこで、上記の手法で検出できなかった場合、更に視細胞を検出するために上記凸面検出を行うようにしてもよい。
【００８４】
また、例えば、輝度値がある閾値以上のピクセルのみ選択し、該選択ピクセルが連結している領域の重心を求める方法がある。他には、小領域の特徴量を算出して、パターン認識の手法を用いて検出することも考えられる。具体的には、例えば１１ｘ１１ピクセルの小領域に対して、ガボール特徴量を算出する。視細胞とされる検出点を中心に含む複数の小領域と、含まない複数の小領域から得られるガボール特徴量ベクトルを用いて、サポートベクターマシンによる学習を行う。新たに対象とする小領域に対してガボール特徴量を算出し、上記学習結果を用いることで、該小領域の中心が視細胞であるか否かを判断する。
【００８５】
さらに、上記ソフトウェアによる検出に加えて、ユーザーによりマニュアルで修正することも可能である。このとき入力情報取得部１１１は、ステップＳ７２０で検出した検出点に対してユーザーが修正した検出点の位置を取得する。
【００８６】
図１５に極大値検出を用いた視細胞検出と上述した凸面検出を用いた視細胞検出との比較の一例を示す。図１５（a）、（ｂ）は極大値検出による視細胞検出結果であり、図１５（ｃ）、（ｄ）は凸面検出を用いた視細胞検出結果である。なお、各図の縦軸は輝度、横軸は眼底における位置である。図１５（a）、（ｃ）を参照すると、上述したように凸面検出を用いた視細胞検出では極大値検出を用いた視細胞検出に比べ多くの点を抽出している（図中、丸印）。また、図１５（ｂ）、（ｄ）を参照すると、極大値検出を用いた視細胞検出では検出ができていない点が凸面検出を用いた視細胞検出では検出できていることが分かる。これは視細胞の特性として隣接する視細胞であっても輝度の差があり、視細胞間の距離が近づくと図１５（ｂ）に示すように極大値が存在しなくなってしまうためである。その結果、図１５（ｂ）に示すように視細胞が存在するのにも関わらず、極大値検出を用いた視細胞検出では視細胞が検出できなくなってしまう。  図１６は黄斑からの距離に対する感度の特性の一例を示す図である。図１６のグラフは、横軸を黄斑（例えば黄斑中心）からの距離、縦軸を感度としている。図１６における感度とは、画像から人が抽出した視細胞に対する視細胞検出により検出された視細胞の割合を示すものである。例えば、人が抽出した視細胞と視細胞検出により検出された視細胞が一致する場合には１となる。なお、縦軸は視細胞検出による精度を測ることができるものではれば他の指標であっても構わない。また図１６において菱形のプロットは凸面検出を用いた視細胞検出の結果を示しバツ印のプロットは極大値検出を用いた視細胞検出の結果を示している。
【００８７】
図１６に示す実験結果からは、凸面検出を用いた視細胞検出の方が極大値検出を用いた視細胞検出よりも精度良く視細胞を検出出来ていることが分かる。また、黄斑からの距離が短くなり視細胞密度が高くなると、極大値検出を用いた視細胞検出の精度は凸面検出を用いた視細胞検出の精度に比べて大きく下がる。これは上述したように、輝度差のある隣接する視細胞の距離が近づくことで、本来視細胞として検出されるべき位置の輝度が極大とならないためである。
【００８８】
一方、凸面検出を用いた視細胞検出では、黄斑からの距離が短くなっても上述のごとく視細胞を検出できていることが分かる。
【００８９】
従って、凸面検出を用いた視細胞検出により精度良く視細胞を検出することが可能となる。より具体的には、凸面検出を用いた視細胞検出により、視細胞検出の精度の黄斑からの位置の依存性を減少させることが可能となる。
【００９０】
なお、黄斑からの距離に応じて凸面検出を用いた視細胞検出と極大値検出を用いた視細胞検出とを使い分けることとしても良い。予め決められた黄斑からの距離以下の場合には凸面検出を用いた視細胞検出を行い、予め決められた黄斑からの距離より大きい場合には極大値検出を用いた視細胞検出を行うこととしても良い。
＜ステップＳ７３０＞
ステップＳ７３０において、視細胞解析部１４３は、ステップＳ７２０で検出された検出点に対して、ボロノイ解析を行う。
【００９１】
具体的には、ステップＳ７２０で検出された全検出点に対して、近傍にある検出点同士の垂直二等分線によって画像内の領域を分割することで、各検出点に属するボロノイ領域を算出する。
【００９２】
＜ステップＳ７４０＞
ステップＳ７４０において、視細胞解析部１４３は、ステップＳ７２０、ステップＳ７３０の解析結果に対し、指標を算出する領域を選択する。
【００９３】
図８（ａ）に図６で０．５ｍｍ  Ｓｕｐｅｒｉｏｒの領域選択の例を、図８（ｂ）に図６で０．５ｍｍ  Ｔｅｍｐｏｒａｌの領域選択の例を示した。領域は矩形や多角形、楕円等を選択することが可能であり、ユーザーによるマウスクリック等の操作に応じて視細胞解析部１４３により設定される。すなわち、視細胞解析部１４３は第１眼底画像に対して領域を設定する設定手段の一例に相当する。なおユーザーの操作に応じるのではなく、視細胞解析部１４３が自動的に領域を設定することとしてもよい。また自動で領域を設定した後にユーザーの操作に応じて領域を変更可能にしてもよい。
【００９４】
なお、血管等が含まれる場合、血管の下の輝度は低くなるため、視細胞の検出精度が悪くなる。よって血管は避けて領域を選択する。すなわち、設定手段の一例である細胞解析部１４３は第１眼底画像に含まれる血管領域を含まないように領域を設定する。
【００９５】
また、視細胞密度が黄斑中心からの距離に依存して変化することから、Ｓｕｐｅｒｉｏｒでは横長に、Ｔｅｍｐｏｒａｌでは縦長に領域を選択している。これはＩｎｆｅｒｉｏｒ、Ｎａｓａｌでも同様であり、Ｉｎｆｅｒｉｏｒは横長、Ｎａｓａｌは縦長とすることで、同一領域内に含まれる視細胞の黄斑中心からの距離のばらつきを抑えることができる。つまり、細胞解析部１４３は、Ｓ像と黄斑中心とを結ぶ直線に垂直な方向における幅が前記直線に平行な方向における幅よりも大きくなるように領域を設定する。すなわち、設定手段の一例である視細胞解析部１４３は第１眼底画像と黄斑中心との位置関係に応じて領域の設定方法を変更する。より具体的には視細胞解析部１４３の一例である設定手段は第１眼底画像と黄斑中心とを結ぶ直線に垂直な方向における幅が直線に平行な方向における幅よりも大きくなるように領域を設定する。
【００９６】
領域の選択法としては、上記のように前処理後画像の中からある大きさのエリアを選択することもできるが、検出点ごとに選択領域を設定することも可能である。例えば、各検出点を中心とした大きさ２０μｍの円を考え、その円内を各検出点に対応した領域と考えることができる。さらに、検出点に限らず、適当な間隔をあけたピクセルごとに設定することも可能である。この場合も、例えば縦横１０ピクセルおきに選択したピクセルに対して、大きさ２０μの円を考え、その円内に含まれる領域を選択領域とすることができる。
【００９７】
さらに、上記選択領域の黄斑中心からの距離を算出する。ステップＳ２５０により、各ＡＯ−ＳＬＯ像の黄斑中心からの位置が算出されている。ステップＳ２５０で算出されたＡＯ−ＳＬＯ像の位置および上記選択領域のＡＯ−ＳＬＯ像内での位置から、選択領域の中心の黄斑中心からの距離を算出する。すなわち、距離取得手段の一例である距離算出部１４２は黄斑中心の位置から領域内の所定の位置（例えば領域の中心）までの距離を取得する。
【００９８】
＜ステップＳ７５０＞
ステップＳ７５０において、視細胞解析部１４３は、ステップＳ７４０で選択された領域に対して、ステップＳ７２０、ステップＳ７３０で取得された結果に基づき、指標を算出する。
【００９９】
ここで指標の例として、ステップＳ７２０で求められた検出点の数（視細胞の数）、検出点数を領域の面積で割った密度（視細胞の密度）、選択した領域における検出点当たりの平均面積、ステップＳ７３０で求められた最近傍検出点までの距離、ボロノイ領域が六角形である割合などがあげられる。また、上記指標とともに、領域の横幅、縦幅および面積を表示することとしてもよい。さらに、理想的な最近傍検出点までの距離や最近傍検出点までの距離に対する実際の最近傍検出点までの距離の比等を表示することとしてもよい。また、黄斑中心からの距離を表示することとしてもよい。
【０１００】
こうして視細胞解析を行った結果及び取得した指標を、制御部１２０を通じて記憶部１３０に保存した後、ステップＳ２６０に戻る。すなわち、ステップＳ２６０において、決定手段の一例である視細胞解析部１４３は各領域までの距離と第１眼底画像に含まれる視細胞に関する情報とから複数の領域毎に評価値を決定する。
【０１０１】
＜ステップＳ２７０＞
評価工程の一例であるステップＳ２７０において、比較部１４４は、ステップＳ２６０で取得された視細胞解析結果の指標を、記憶部１３０に保存されている基準情報の一例である正常眼データと比較することで。
【０１０２】
図９（ａ）に正常眼における、視細胞に関する情報の一例である視細胞の密度と黄斑中心からの距離の関係を示した例を示す。グラフで塗りつぶされた範囲が正常眼のばらつきの範囲とする。ここに、ステップＳ２６０で取得した、図８（ｂ）の選択領域に対する視細胞密度の指標をプロットしたものの一例を図９（ｂ）に示す。図９（ｂ）に示されているように、黄斑中心からの距離に応じた正常眼データと検出した視細胞密度とが比較される。選択領域のうち６か所は正常範囲に含まれているが、一か所は正常範囲から逸脱していることがわかる。すなわち、比較部１４４は、黄斑（例えば黄班中心）から第１眼底画像における所定の位置までの距離と第１眼底画像に含まれる視細胞に関する情報とから被検眼の状態を評価する評価手段の一例に相当する。より具体的には、評価手段の一例である比較部１４４は、領域内の所定の位置までの距離と領域における視細胞に関する情報とから被検眼の状態を評価する。また、上述したように評価手段の一例である比較部１４４は、視細胞に関する情報と距離に応じた基準情報とを比較する。
【０１０３】
正常範囲から逸脱している場合に、どの領域が正常範囲から逸脱しているかが明示されるように、図８（ｂ）上で正常範囲から逸脱した領域が図１０（ａ）のように明示されるように出力部１５０は表示部であるモニタ等に表示させる。すなわち、出力部１５０は比較の結果に基づく表示形態で第１眼底画像を表示部に表示させる表示制御手段の一例に相当する。例えば、図９（ｂ）と図１０（ａ）とを並べて表示部に表示することとすれば、正常範囲から逸脱した値と選択領域とが容易に把握可能となる。
【０１０４】
なお、図９（ｂ）のグラフと図８（ｂ）の領域は連動させることができ、図８と図９とを表示部に並べて表示し、例えば図９（ｂ）のグラフ上のポイントとクリックすると、図８（ｂ）上の対応する領域が輪郭強調されるなどすることで、対応関係をより明確に示すことができる。
【０１０５】
さらに、ステップＳ７４０の領域選択において、検出点ごと、もしくは適当な間隔のピクセルごとに選択領域を設定した場合について考える。例えば縦横１０ピクセルおきに領域選択がされている場合、その領域ごとに正常眼の分布と比較することがでるため、図８（ｂ）に対してそのような処理を行った例を図１０（ｂ）に示す。ここでは正常分布の平均値となる場合を基準に、そこからの逸脱量に応じて濃い色で表示するようにしている。なお、正常範囲から逸脱するほど薄い色としてもよいし、正常範囲から逸脱するほど濃い色としてもよい。すなわち、表示制御手段の一例である出力部１５０は、視細胞に関する情報と距離に応じた基準情報との差に応じて前記第１眼底画像の明度を変える。
【０１０６】
また、出力部１５０は、図９（ａ）で示された正常分布の範囲外となる領域の面積を計測し、提示することも可能である。すなわち、表示制御手段の一例である出力部１５０は、視細胞に関する情報と距離に応じた基準情報との差が所定値以上となる第１眼底画像上の領域の面積を提示する。
【０１０７】
図１１は、各ＡＯ−ＳＬＯ像に対して図１０（ｂ）で示された画像を生成し、それを図６の重畳画像に表示した例である。このとき、疾病の進行度を示す指標として、例えば黄斑中心からの距離の関数として、異常領域の面積の割合を示すことができる。なお、図１０（ａ）を図６の重畳画像に表示することとしてもよい。
【０１０８】
＜ステップＳ２８０＞
ステップＳ２８０において、入力情報取得部１１１は、ステップＳ２７０で提示された解析結果及びその比較結果に対して、ユーザーが修正を行うか否かの判断を取得する。ここでユーザーの修正とは、ステップＳ２３０で行った各ＡＯ−ＳＬＯ像の位置合わせや、ステップＳ２４０で行った黄斑中心の検出に対してである。
【０１０９】
具体的には、ステップＳ２３０のＡＯ−ＳＬＯ像の位置合わせが間違っていると判断した場合には、ユーザーが該ＡＯ−ＳＬＯ像の位置を変更することができる。さらに、ステップＳ２４０で決められた黄斑中心の位置を変更することができる。
【０１１０】
ユーザーによる修正がない場合には、ステップＳ２９０に進む。修正が行われた場合には、ステップＳ２５０に戻り、修正された黄斑中心位置及び各ＡＯ−ＳＬＯ像の位置に基づき再度距離計算を行い、それ以降の処理を行う。
【０１１１】
＜ステップＳ２９０＞
ステップＳ２９０において、制御部１２０は、記憶部１１０に保存されている、算出した指標や比較結果など記憶部１３０に記憶されている処理結果を、データベースに保存する。
【０１１２】
以上の構成により、補償光学ＳＬＯ装置により取得した複数のＡＯ−ＳＬＯ像を解析して指標を算出する場合に、同一ＡＯ−ＳＬＯ像内に含まれない黄斑中心との位置関係を考慮しながら、解析及を行うことが可能になる。
【０１１３】
（実施例２）
実施例１では、解像度の異なる複数のＡＯ−ＳＬＯ像から作成した重畳画像を用いて、解析対象となるＡＯ−ＳＬＯ像の黄斑中心からの距離を取得し、それを反映させた解析及び評価を行う例を示した。
【０１１４】
実施例２では、重畳画像より検出できる病変等の情報を利用して、解析及び評価を行う場合について説明する。
【０１１５】
視細胞に関連した遺伝性の疾病では、疾病の進行に伴い特有のリング状の構造物が観察されることが知られている。このようなリング構造は、高解像度かつ低画角の画像で把握することは困難であり、低解像度であっても高画角の画像により把握することが可能である。このような構造を高画角画像より抽出し、それに基づき高解像度画像における処理を変える方法について説明する。
【０１１６】
本実施例に関わる画像処理装置１０の機能構成は、図１と同様なので説明は省略する。
【０１１７】
図１２のフローチャートを参照して、本実施形態の画像処理装置１０の処理手順を説明する。ここでステップＳ２１０、Ｓ２２０、Ｓ２３０、Ｓ２８０、Ｓ２９０に関しては実施例１で説明した処理手順と変わらないために、説明は省略する。
【０１１８】
＜ステップＳ１２４０＞
ステップＳ１２４０において、入力情報取得部１１１は、ステップＳ２２０で生成された重畳画像に対して、ユーザーが検出した黄斑中心の位置及び病変位置を取得する。そうして検出した黄斑中心位置及び病変位置を、制御部１２０を通じて記憶部１３０に保存する。
【０１１９】
黄斑中心を検出する方法については実施例１で説明した方法と同一である。リング構造を検出する方法としては、重畳画像からユーザーが目視により検出する方法が考えられる。図１３にリング構造が見られる重畳画像の例を示した。但し図１３では、リング構造が見られる領域が小さいため、重畳画像のうち該リング構造を含む４枚のＬ像より作成した領域を拡大して示している。このようにして示されたリング構造の概形にそって、ユーザーが重畳画像上をマウスでクリックしながら概形を取得することができる。
【０１２０】
さらに、眼底自発蛍光等の画像よりリング構造を抽出し、血管等の位置を参考にしながら、重畳画像上の対応する位置に置くことで、リング構造を取得する方法などもある。
【０１２１】
＜ステップＳ１２５０＞
ステップＳ１２５０において、距離算出部１４２は、ステップＳ２３０で算出された各ＡＯ−ＳＬＯ像のＷＦ−ＳＬＯ像からの移動量と、ステップＳ１２４０で取得された黄斑中心位置から、各ＡＯ−ＳＬＯ像の、黄斑中心を原点とした座標を算出する。さらに、ステップＳ１２４０で取得した病変位置の、黄斑中心を原点とした座標を算出する。
【０１２２】
各ＡＯ−ＳＬＯ像の座標の求め方は実施例１で説明した手法と同じである。病変部に関しては、ステップＳ１２４０でマウスクリック等により取得した座標値を、重畳画像上の長さに変換し、かつ原点と黄斑中心位置とすることで取得できる。
【０１２３】
そして求められた各ＡＯ−ＳＬＯ像、及び病変の黄斑中心からの位置を、制御部１２０を通じて記憶部１３０に保存する。
【０１２４】
＜ステップＳ１２６０＞
ステップＳ１２６０において、視細胞解析部１４３は、ステップＳ２１０で取得された収差補正ＳＬＯ画像のＳ画像に対して、視細胞解析を行う。
【０１２５】
視細胞解析のステップは実施例１で説明した方法に従うが、ここではステップＳ１２４０で取得した病変情報に基づき、処理を変化させる箇所についてのみ記載する。
【０１２６】
一般に、リング構造の外側と内側では疾病の進行状態が異なるため、視細胞の描出のされ方に違いが見られる。正常に近い領域がリングの内側になるか外側になるかは疾病によって異なるが、疾病が進行している領域では一般に視細胞が変性しているため、視細胞以外の構造物が描出されていたり、輝度の低い領域が広く描出されていたりする場合がある。このため実施例１で示した極大値検出アルゴリズムを適応すると、多くのノイズを拾うためにアルゴリズムのロバスト性が下がってしまう。
【０１２７】
このような状況を防ぐために、視細胞解析を行う領域が、リング構造の疾病が進行している側となる場合には、ステップＳ７１０の前処理の段階で追加の平均値フィルタ処理を入れる。ここでは平均値フィルタを用いているが、フィルタの種類は平均値フィルタに限定されるものではなく、ノイズ低減につながるものであれば、ガウシアンフィルタなどさまざまな種類のフィルタを用いることが可能である。
【０１２８】
また、ステップＳ７２０の検出の際の閾値のパラメタを、リングの内側と外側で変更する。具体的には、リングの内側では実施例１の場合と同様、画像の輝度最低値とするが、リングの外側では画像の輝度平均値を用いるなどである。このように変更する理由は上記と同様、リングの内側もしくは外側で明確は視細胞が描出されていない場合の、検出のロバスト性を向上させるためである。
【０１２９】
さらに、ステップＳ７２０でパターン認識の手法による検出の例を示したが、この場合に、疾病ごとにリングの外側と内側で別に学習をした結果を用いることが望ましい。
【０１３０】
＜ステップＳ１２７０＞
ステップＳ１２７０において、比較部１４４は、ステップＳ１２６０で取得された視細胞解析結果の指標を、記憶部１３０に保存されている正常眼データと比較する。
【０１３１】
実施例１におけるステップＳ２７０と同様の処理により、図１３上で示されている各Ｓ像に対して、図１１と同様の処理を行った結果を図１４に示す。図１４では、重畳画像上に示されている病変部と、視細胞解析結果を正常眼分布と比較して作成した画像との間に矛盾が見られる。具体的には、図１４の右側のＳ像は、リング構造の内側であるにもかかわらず、正常分布からおおきく外れているのに対し、右側はリング構造付近であっても比較的正常に近い分布となっている。このような場合には、各Ｓ像の重畳画像内での位置が間違っているかもしれないし、もしくはリング構造の取得が正しく行われていないかもしれない。または、実際に重畳画像上でリング構造にように見えているものが、実際の病変を正しく反映していないかもしれない。前者の位置合わせやリング構造取得の問題であるならば、ステップＳ１２８０で修正することができる。しかし後者である場合には、他モダリティ画像との比較を通じた解析が必要であり、図１４の表示を行うことによりユーザーが病態解明を進めるサポートとなる。
【０１３２】
＜ステップＳ１２８０＞
ステップＳ１２８０において、入力情報取得部１１１は、ステップＳ１２７０で提示された解析結果及びその比較結果に対して、ユーザーが修正を行うか否かの判断を取得する。ここでユーザーの修正とは、ステップＳ２３０で行った各ＡＯ−ＳＬＯ像の位置合わせや、ステップＳ１２４０で行った黄斑中心及び病変部の検出に対してである。
【０１３３】
ユーザーが修正を行うと判断した場合には、ユーザーは再度ＡＯ−ＳＬＯ像の位置や黄斑中心、及び病変の位置を設定してステップＳ１２５０に戻り、修正された黄斑中心や病変の位置及び各ＡＯ−ＳＬＯ像の位置に基づき再度距離計算を行い、それ以降の処理を行う。ユーザーによる修正がない場合には、ステップＳ２９０に進む。
【０１３４】
（その他の実施形態）
本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記憶した記憶媒体を、システムあるいは装置に供給する。そして、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読出し実行することによっても、達成されることは言うまでもない。
【符号の説明】
【０１３５】
１０  画像処理装置
１００  画像取得部
１１０  入力情報取得部
１２０  制御部
１３０  記憶部
１４０  画像処理部
１４１  位置合わせ部
１４２  距離算出部
１４３  視細胞解析部
１４４  比較部
１５０  出力部
【図１】
2014198223.tif 000007
【図２】
2014198223.tif 000008
【図３】
2014198223.tif 000009
【図４】
2014198223.tif 000010
【図５】
2014198223.tif 000011
【図６】
2014198223.tif 000012
【図７】
2014198223.tif 000013
【図８】
2014198223.tif 000014
【図９】
2014198223.tif 000015
【図１０】
2014198223.tif 000016
【図１１】
2014198223.tif 000017
【図１２】
2014198223.tif 000018
【図１３】
2014198223.tif 000019
【図１４】
2014198223.tif 000020
【図１５】
2014198223.tif 000021
【図１６】
2014198223.tif 000022