(57)【特許請求の範囲】
【請求項１】
眼底の視細胞の画像を処理する画像処理装置であって、
前記視細胞の画像を周波数変換した空間周波数画像を取得する周波数変換手段と、
前記空間周波数画像における複数の方向の強度情報を取得する取得手段と、
を有することを特徴とする画像処理装置。
【請求項２】
前記取得手段が、前記空間周波数画像における前記複数の方向の強度分布を前記複数の方向の強度情報として取得することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記取得手段が、前記空間周波数画像における環状構造の歪みを示す値を前記複数の方向の強度情報として取得することを特徴とする請求項１あるいは２に記載の画像処理装置。
【請求項４】
前記歪みを示す値が、前記空間周波数画像における所定の方向の積算値と該所定の方向に交差する方向の積算値とからなることを特徴とする請求項３に記載の画像処理装置。
【請求項５】
前記複数の方向の強度情報を示す表示形態を表示手段に表示させる表示制御手段を有することを特徴とする請求項１乃至４のいずれか１項に記載の画像処理装置。
【請求項６】
前記表示制御手段が、前記複数の方向の強度情報を示す表示形態と、前記視細胞の画像と前記空間周波数画像とのうち少なくとも一方の画像とを該複数の方向それぞれに対応付けて前記表示手段に表示させることを特徴とする請求項５に記載の画像処理装置。
【請求項７】
前記周波数変換手段が、前記眼底の視細胞の複数の画像を周波数変換し、複数の空間周波数画像を取得し、
前記取得手段が、前記複数の空間周波数画像それぞれにおける前記複数の方向の強度情報を取得することを特徴とする請求項１乃至６のいずれか１項に記載の画像処理装置。
【請求項８】
前記強度情報に基づいて複数の画像を分類する画像分類手段と、
前記複数の画像の分類の結果に基づいて画像を選択する画像選択手段と、
を有することを特徴とする請求項７に記載の画像処理装置。
【請求項９】
前記選択された画像を重ね合わせる重ね合わせ手段を有し、
前記画像選択手段は、基準画像と、該基準画像に重ね合わせる画像を選択することを特徴とする請求項８に記載の画像処理装置。
【請求項１０】
眼底の視細胞の画像を処理する画像処理装置であって、
前記視細胞の画像を前記視細胞の周期性を示す画像に変換する変換手段と、
前記周期性を示す画像における複数の方向の強度情報を取得する取得手段と、
を有することを特徴とする画像処理装置。
【請求項１１】
前記取得手段は、前記周期性を示す画像の第１の方向に沿った画素の強度と前記第１の方向に交差する第２の方向に沿った画素の強度とを用いて、前記周期性を示す画像から前記強度情報を取得することを特徴とする請求項１０に記載の画像処理装置。
【請求項１２】
眼底の視細胞の画像を処理する画像処理装置であって、
前記視細胞の画像を前記視細胞の周期性を示す画像に変換する変換手段と、
前記周期性を示す画像の環状構造の歪みに基づいて、前記周期性を示す画像から特徴量を抽出する抽出手段と、
を有することを特徴とする画像処理装置。
【請求項１３】
請求項１乃至１２のいずれか１項に記載の画像処理装置に通信可能に接続された補償光学ＳＬＯ装置により前記眼底の視細胞を撮影して得た前記視細胞の画像を取得する画像取得手段を更に有することを特徴とする画像処理装置。
【請求項１４】
請求項１乃至１２のいずれか１項に記載の画像処理装置にネットワークを介して接続されたデータベースから前記データベースに保存された前記視細胞の画像を取得する画像取得手段を更に有することを特徴とする画像処理装置。
【請求項１５】
眼底の視細胞の画像を処理する画像処理方法であって、
前記視細胞の画像を周波数変換し、空間周波数画像を取得する工程と、
前記空間周波数画像における複数の方向の強度情報を取得する工程と、
を有することを特徴とする画像処理方法。
【請求項１６】
前記取得する工程では、前記空間周波数画像における前記複数の方向の強度分布を前記複数の方向の強度情報として取得することを特徴とする請求項１５に記載の画像処理方法。
【請求項１７】
前記取得する工程では、前記空間周波数画像における環状構造の歪みを示す値を前記複数の方向の強度情報として取得することを特徴とする請求項１５あるいは１６に記載の画像処理方法。
【請求項１８】
眼底の視細胞の画像を処理する画像処理方法であって、
前記視細胞の画像を前記視細胞の周期性を示す画像に変換する工程と、
前記周期性を示す画像における複数の方向の強度情報を取得する工程と、
を有することを特徴とする画像処理方法。
【請求項１９】
前記取得する工程では、前記周期性を示す画像の第１の方向に沿った画素の強度と前記第１の方向に交差する第２の方向に沿った画素の強度とを用いて、前記周期性を示す画像から前記強度情報を取得することを特徴とする請求項１８に記載の画像処理方法。
【請求項２０】
眼底の視細胞の画像を処理する画像処理方法であって、
前記視細胞の画像を前記視細胞の周期性を示す画像に変換する工程と、
前記周期性を示す画像の環状構造の歪みに基づいて、前記周期性を示す画像から特徴量を抽出する工程と、
を有することを特徴とする画像処理方法。
【請求項２１】
請求項１５乃至２０のいずれか１項に記載の画像処理方法の各工程をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、眼底の画像を処理する画像処理装置及び画像処理方法に関する。
【背景技術】
【０００２】
生活習慣病や失明原因の上位を占める疾病の早期診断を目的として、被検眼の眼底の検査が広く行われている。共焦点レーザー顕微鏡の原理を利用した眼科装置である走査型レーザー検眼鏡（ＳＬＯ：Ｓｃａｎｎｉｎｇ  Ｌａｓｅｒ  Ｏｐｈｔｈａｌｍｏｓｃｏｐｅ）は、測定光であるレーザーを眼底に対してラスタースキャンを行い、その戻り光の強度から平面画像を高分解能かつ高速に得る装置である。近年、被検眼の収差を波面センサでリアルタイムに測定し、被検眼にて発生する測定光やその戻り光の収差を波面補正デバイスで補正する補償光学系を有する補償光学ＳＬＯが開発され、高横分解能な平面画像の取得を可能にしている。さらに、取得した網膜の平面画像を用いて網膜における視細胞に対応する箇所を抽出し、その密度や分布の解析から疾病の診断や薬剤応答の評価が試みられている。
【０００３】
補償光学ＳＬＯを用いて視細胞の可視化を行った例として、非特許文献１においては、網膜の平面画像を取得し、視細胞に対応する箇所の自動抽出を可能とした眼科撮影装置が提案されている。この眼科撮影装置では、高横分解能の網膜の平面画像を撮像し、画像に描出される視細胞の配置の周期性を利用して高周波成分を除去することで画像の前処理を行い、視細胞の自動検出を行っている。また視細胞の検出結果を用いて、視細胞の密度や視細胞間の距離を計測し、空間内分布のボロノイ解析による分析を行っている。
【０００４】
視細胞解析を精度よく行うためには、画質の高い画像を取得することが重要である。そのための一つの手法として、同一箇所を連続して撮影した複数の画像を重ね合わることで、ランダムノイズを除去する方法が知られている。画像の重ね合わせ処理を行うためには、複数枚の画像の中から基準となる画像（以下基準画像とする）を選択し、該基準画像に対してそれ以外の画像の位置あわせ処理を行った後に、重ね合わせ処理を行う方法が一般的である。
【０００５】
ここで、被検眼には、固視微動と呼ばれるランダムな動きが常に存在する。このため、同一箇所を連続して撮影した場合でも、撮影位置は一定ではない。そこで、基準画像を用いて、複数の断層画像の位置合わせを正確に行う手法が、特許文献１に開示されている。
【０００６】
ところで、固視微動には、数μ程度の動きで振動するトレモアや、数１０μ程度のフリック、常に流れるように移動するドリフト、非常に速い動きとなるサッケードなどが知られている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００８−２９３１８５号公報
【非特許文献】
【０００８】
【非特許文献１】Ｋａｃｃｉｅ  Ｙ．  Ｌｉ  ａｎｄ  Ａｕｓｔｉｎ  Ｒｏｏｒｄａ，  “Ａｕｔｏｍａｔｅｄ  ｉｄｅｎｔｉｆｉｃａｔｉｏｎ  ｏｆ  ｃｏｎｅ  ｐｈｏｔｏｒｅｃｅｐｔｏｒｓ  ｉｎ  ａｄａｐｔｉｖｅ  ｏｐｔｉｃｓ  ｒｅｔｉｎａｌ  ｉｍａｇｅｓ”Ｊ．  Ｏｐｔ．  Ｓｏｃ．  Ａｍ．  Ａ，  Ｍａｙ  ２００７，  Ｖｏｌ．  ２４，  Ｎｏ．  ５，  １３５８
【発明の概要】
【発明が解決しようとする課題】
【０００９】
固視微動のある眼底画像群に対して画像の重ね合わせ処理を行うためには、基準画像としてフリック等の大きな動きの最中に撮影された画像を選択しないことや、フリックとフリックの間のトレモア中に撮影された画像群を重ね合わせ画像として選択することが重要になる。このため、眼底の視細胞の画像の位置ずれの有無を簡単に判別することが求められる。
【課題を解決するための手段】
【００１０】
本発明に係る画像処理装置の一つは、
眼底の視細胞の画像を処理する画像処理装置であって、
前記視細胞の画像を周波数変換した空間周波数画像を取得する周波数変換手段と、
前記空間周波数画像における複数の方向の強度情報を取得する取得手段と、
を有することを特徴とする。
【発明の効果】
【００１１】
本発明によれば、眼底の視細胞の画像を周波数変換した周波数画像における複数の方向の強度情報を取得することができる。このとき、例えば、この強度情報を示す表示形態をモニタ等の表示部に表示させることができる。これにより、眼底の視細胞の画像の位置ずれの有無を簡単に判別することができる。
【図面の簡単な説明】
【００１２】
【図１】第１の実施形態に係る画像処理装置のブロック図である。
【図２】第１の実施形態に係る画像処理装置の画像処理のフロー図である。
【図３】第１の実施形態に係る補償光学ＳＬＯ装置により、視細胞を撮影した高精細な平面画像の模式図である。
【図４】第１の実施形態に係る画像処理装置により、平面画像を周波数変換して得られるフーリエ画像の一例を説明する図である。
【図５】第１の実施形態に係る画像処理装置により、歪みの少ない平面画像に対応するフーリエ画像から特徴量を抽出する一例を説明する図である。
【図６】第１の実施形態に係る画像処理装置により、歪みの大きな平面画像に対応するフーリエ画像から特徴量を抽出する一例を説明する図である。
【図７】第２の実施形態に係る画像処理装置の画像処理のフロー図である。
【図８】第２の実施形態に係る画像処理装置において、特徴量に基づいて平面画像群を被検眼の固視微動の種類毎に分類することを説明する図である。
【図９】第３の実施形態に係る補償光学ＳＬＯ装置の一例である複合装置の全体の構成について説明する図である。
【図１０】第３の実施形態に係る補償光学ＳＬＯ装置の一例である複合装置により、断層画像を取得する方法を説明する図である。
【図１１】第３の実施形態に係る補償光学ＳＬＯ装置の一例である複合装置により、平面画像を取得する方法を説明する図である。
【発明を実施するための形態】
【００１３】
本実施形態に係る画像処理装置は、被検眼の眼底の視細胞の画像を周波数変換した周波数画像を取得する周波数変換部を有する。また、本実施形態に係る画像処理装置は、上記周波数画像における複数の方向（例えば、ＸＹ方向）の強度情報を取得する取得部を有する。ここで、上記強度情報は、例えば、図５や図６のＸＹ方向それぞれの強度分布のことである。また、上記強度情報は、例えば、環状構造（リング構造）の歪みを示す値（画像の歪みを示す特徴量や指標Ｄ）のことである。これにより、視細胞の画像内における位置ずれの有無を簡単に判別することができる。すなわち、重ね合わせに用いる画像として、位置ずれの無い画像を選択することができる。
【００１４】
ここで、上記複数の方向の強度情報を示す表示形態をモニタ等の表示部に表示させる表示制御部を有することが好ましい。これにより、ユーザは、視細胞の画像内における位置ずれの有無を簡単に判別することができる。このとき、この位置ずれの有無を表示部に表示させても良い。また、上記複数の方向の強度情報を示す表示形態と、視細胞の画像と周波数画像とのうち少なくとも一方の画像とを該複数の方向それぞれに対応付けて表示部に表示させることが好ましい。これにより、ユーザは、例えば、ＸＹ方向のいずれに位置ずれが生じたのかを簡単に判別することができる。
【００１５】
（第１の実施形態）
本実施形態に係る補償光学ＳＬＯ装置は、被検眼の網膜の同一箇所を連続して撮影することにより、被検眼の眼底の複数の画像を取得することができる。このとき、できるだけ同一箇所が撮影された画像同士で、かつ一枚の画像内における被検眼の固視微動による歪みが少ない画像同士で重ね合わせすることが望ましい。これにより、高画質な画像を取得することができる。
【００１６】
このとき、画像同士の位置合わせのための基準画像や重ね合わせる複数の画像を選択することが好ましい。具体的には、補償光学ＳＬＯにより取得された複数の眼底平面画像（以下、平面画像）の空間周波数画像を離散フーリエ変換によりそれぞれ取得する（取得された画像を以下、フーリエ画像とする）。該取得されたフーリエ画像から、画像内での歪みを反映した特徴量を取得し、該特徴量の分布から規準画像、及び重ね合わせに用いる画像群を選択する。こうして得られた基準画像及び重ね合わせ画像を用いることにより、重ね合わせ処理によって高画質な画像を取得することが可能になる。
【００１７】
＜平面画像＞
図３に、補償光学ＳＬＯで撮影された平面画像を模式的に示す。図３に示すように、輝度が比較的大きい小領域として、各視細胞ＰＲが一つ一つ区別される形で描出される。また、視細胞の輝度に比較して輝度の低い血管領域Ｖが描出されることもある。この血管領域Ｖは視細胞より上層にある血管の影が描出されているものである。
【００１８】
＜フーリエ画像＞
図４（ａ）に、上記平面画像の空間周波数成分を離散フーリエ変換により取得したフーリエ画像の模式図を示す。図４（ａ）に示すように、視細胞が周期的に配置していることを反映して、視細胞の周期に相当するリング構造が見られる。
【００１９】
＜画像処理装置の構成＞
図１は、本実施形態に係る画像処理装置２０の機能構成を示したものである。
【００２０】
まず、１０は画像取得部であり、補償光学ＳＬＯ装置から網膜の同一箇所を撮影した複数の平面画像を取得する。１１は入力情報取得部であり、補償光学ＳＬＯによる平面画像の撮影の際の、被検眼情報を取得する。取得した画像は制御部１２を通じて記憶部１３に記憶される。１４は画像処理部であり、周波数変換部２１、特徴抽出部２２、画像選択部２３、重ね合わせ部２４を含む。画像処理部１４は、取得した平面画像よりフーリエ画像を生成し、フーリエ画像より抽出した特徴量から重ね合わせの基準画像及び重ね合わせに用いる画像群を選択する。それらの画像を用いて重ね合わせ画像を作成し、記憶部１３に記憶する。１５は出力部であり、作成された重ね合わせ画像をモニタ等に出力するほか、記憶部１３に記憶されている処理結果をデータベースに出力する。
【００２１】
ここで、画像処理装置２０は、補償光学ＳＬＯから直接平面画像を取得する構成としたが、ネットワークを介してそれらの画像を取得してもよい。その場合には、ネットワークを介して接続されたデータベースに補償光学ＳＬＯで撮影された複数の平面画像が保存され、画像処理装置２０はそのデータベースから画像を取得する構成となる。
【００２２】
＜画像処理装置の処理手順＞
次に、図２（ａ）のフローチャートを参照して、本実施形態の画像処理装置２０の処理手順を説明する。
【００２３】
＜ステップＳ２１０＞
ステップＳ２１０において、画像取得部１０は、画像処理装置２０に接続される補償光学ＳＬＯ、もしくは該装置により撮影された平面画像を保存するデータベースから、解析対象となる複数の平面画像を取得する。取得された複数の平面画像は、制御部１２を通じて記憶部１３に記憶する。
【００２４】
またこのとき、取得された複数の平面画像が撮影された際の撮影パラメータ情報を取得し、制御部１２を通じて記憶部１３に記憶する。ここで撮影パラメータ情報とは、撮影の際の固視灯の位置情報などのことである。固視灯の位置情報などの撮影パラメータ情報は、平面画像に付加している画像撮影情報ファイルに記載されている場合もあれば、画像のタグ情報として含まれている場合もある。
【００２５】
＜ステップＳ２２０＞
ステップＳ２２０において、入力情報取得部１１は、データベースもしくは操作者の入力部（不図示）による入力により被検眼の情報を取得する。ここで被検眼の情報とは、被検眼の患者ＩＤ、氏名、年齢、性別、検査対象が右目であるか左目であるか、撮影日時等のことであり、これら取得した情報は、制御部１２を通じて記憶部１３に記憶される。
【００２６】
＜ステップＳ２３０＞
ステップＳ２３０において、周波数変換部２１は、記憶部１３に記憶されている補償光学ＳＬＯにより取得された複数の平面画像のそれぞれに対し、空間周波数画像を離散フーリエ変換により取得する。図３に示したように、平面画像の多くの領域は、輝度の高い小領域として観察される規則的な配置をもつ視細胞で構成される。よって該平面画像を空間周波数変換して得られるフーリエ画像は、図４（ａ）に示したようなリング状の構造をもつ。
【００２７】
＜ステップＳ２４０＞
ステップＳ２４０において、特徴抽出部２２は、ステップＳ２３０で得られたフーリエ画像から、画像の歪みを示す特徴量を抽出する。そうして算出した特徴量を、制御部１２を通じて抽出した平面画像と対応付けて記憶部１３に記憶する。以下にステップＳ２４０によって取得する特徴量の一例を示すが、取得する特徴量は以下の例に限定されるものではない。
【００２８】
図４（ｂ）に画面内で歪みが見られる平面画像のフーリエ画像を示す。図４（ｂ）では、図４（ａ）で示されるようなリング状の構造が大きく歪んだ状態となることがわかる。
【００２９】
そこで本実施形態では、画像の歪みを示す特徴量として、フーリエ画像のｘ方向の積算値（所定の方向の積算値の一例）とｙ方向の積算値（該所定の方向に交差する方向の積算値の一例）を考える。具体的には、図５（ａ）に示すように、フーリエ画像のサイズを縦と横のピクセル数がＮとなるＮｘＮの正方形であるとした時、フーリエ画像を縦方向に積算した関数Ｉｙｓｕｍ（ｘ）と、横方向に積算した関数Ｉｘｓｕｍ（ｙ）を考える。
【００３０】
【数１】
0005847510.tif 000002
【００３１】
ここで、Ｉ（ｘ、ｙ）はフーリエ画像の（ｘ、ｙ）の位置にあるピクセルの輝度値を表わす。そして画像の歪みの指標Ｄとして、Ｉｙｓｕｍ（ｘ）とＩｘｓｕｍ（ｙ）の差を考える。
【００３２】
【数２】
0005847510.tif 000003
【００３３】
ここで、｜｜は絶対値を表わし、ｆａｖｅｒａｇｅは元の平面画像の、輝度の平均値を示す。またフーリエ画像内でｘ＝０、ｙ＝０に対応する値にはエッジの影響等のノイズが多くのるため、上記和を計算する場合に除くなどの処理をしている。
【００３４】
図５（ｂ）に示すように、図４（ａ）のような歪みのない画像の場合のフーリエ画像は円形のリング構造が現れるために、Ｉｙｓｕｍ（ｘ）と、Ｉｘｓｕｍ（ｙ）の値がほぼ等しくなるためＤは小さな値となる。これに対して図６に示すように、図４（ｂ）のような歪みのある画像に対応するフーリエ画像は、リング構造が円形ではなく歪んだ形となるために、Ｉｙｓｕｍ（ｘ）と、Ｉｘｓｕｍ（ｙ）の値が大きくずれることとなり、指標Ｄは大きな値となる。
【００３５】
＜ステップＳ２５０＞
ステップＳ２５０において、画像選択部２３は、ステップＳ２４０で取得された特徴量から、複数の同一箇所を撮影した補償光学ＳＬＯ画像の中から、重ね合わせの基準となる基準画像を選択する。そうして選択した基準画像を、制御部１２を通じて記憶部１３に記憶する。
【００３６】
ステップＳ２４０において、画像の歪みを示す特徴量として指標Ｄを取得した。基準画像としてはＤが小さなものが望ましい。よって基準画像の選択方法としては、Ｄが最も小さくなる画像を選択する。なお、画像の画質を示す他の指標、例えばコントラストや輝度の平均値などの指標を用いて、基準画像の候補画像を選択し、その中から指標Ｄが大きなものを除いてもよい。
【００３７】
具体的には、複数の同一箇所を撮影した平面画像群の枚数がｍ枚であるとする。そのとき、各平面画像に対してステップＳ２４０で示した指標Ｄを求め、平面画像群ｍ内での平均値と分散を取得する。
【００３８】
【数３】
0005847510.tif 000004
【００３９】
ここで、各平面画像に対して画像の画質を示す指標をＣｉとする。基準画像の選択は、Ｃｉが最も高いもので、かつＤｉがＤの平均値＋標準偏差以下であるものを選択するとする。
【００４０】
更に、選択した基準画像に重ね合わせる画像群を、各画像のＤに基づいて選択する。具体的には、各画像のＤがＤの平均値＋標準偏差以下である画像を重ね合わせに用いる画像として選択する。
【００４１】
＜ステップＳ２６０＞
ステップＳ２６０において、重ね合わせ部２４は、ステップＳ２５０で選択された基準画像に対して他の平面画像の位置合わせ処理を行う。そして取得された位置合わせ後の画像群から重ね合わせ画像を作成し、制御部１２を通じて記憶部１３に記憶する。
【００４２】
ここで、同一箇所を撮影した複数画像の位置合わせ手法については多くの手法が知られている。画像に歪みが生じない剛体変換の場合には、移動方向のみを取得する位相限相関法が知られている。多少の歪みを生じる場合には、アフィン変換によって歪みの補正を行う手法がある。被写体がより大きな変形を伴う場合には、様々な非剛体の位置合わせ手法が知られていると同時、医用画像処理の分野で研究がさかんである。本実施形態と同様に眼球の固視微動に伴う歪みを課題とする位置合わせとしては、例えば、Ｓｃｏｔｔ  Ｂ．  Ｓｔｅｖｅｎｓｏｎ  ａｎｄ  Ａｕｓｔｉｎ  Ｒｏｏｒｄａ，  “Ｃｏｒｒｅｃｔｉｎｇ  ｆｏｒ  ｍｉｎｉａｔｕｒｅ  ｅｙｅ  ｍｏｖｅｍｅｎｔｓ  ｉｎ  ｈｉｇｈ  ｒｅｓｏｌｕｔｉｏｎ  ｓｃａｎｎｉｎｇ  ｌａｓｅｒ  ｏｐｈｔｈａｌｍｏｓｃｏｐｙ”Ｐｒｏｃ．  ＳＰＩＥ  ５６８８，  １４５  （２００５）では、補償光学ＳＬＯで撮影した網膜画像の位置合わせをする方法が紹介されている。
【００４３】
本実施形態においては、基準画像に対して、各平面画像の画像全体の移動方向を位相限相関法により取得した後、アフィン変換により歪み補正を行う手法を用いたが、位置合わせ処理の方法はこの手法に限定されるものではなく、上記紹介したすべての手法に対して適応することが可能である。
【００４４】
このようにして取得された位置合わせ済みの画像群を用いて、重ね合わせ画像を生成する。位置合わせ済みの画像群の画像サイズがＭｘＬであるとする。位置合わせ前の平面画像群のサイズがＮｘＮであったとすると、Ｍ＞Ｎ、Ｌ＞Ｎである。基準画像は位置合わせ処理前後で変化しないため、基準画像のフレーム番号がｎであったとすると、位置合わせ済み画像群のフレームｎ番目には、画像サイズＭｘＬの中に、基準画像と同じＮｘＮの画像が含まれていて、その周囲のピクセルは値０となっている。重ね合わせ画像の作成は、このＭｘＬ中のＮｘＮの領域に対して、位置合わせ済み画像群の他のフレームの値を取得しその平均値を求めることで生成する。このとき、すべての画像を重ね合わせるのではなく、平面画像群中で例えば瞬きなどで撮影に失敗したフレームや、位置合わせ処理中、位置合わせがうまくいかなかったフレームを除去して重ね合わせを行うと、ノイズレベルの低い重ね合わせ画像を生成することができる。
【００４５】
＜ステップＳ２７０＞
ステップＳ２７０において、出力部１７０は、ステップＳ２６０において記憶部に保存された重ね合わせ画像を、出力部１７０を通じてモニタに表示する。さらにステップＳ２１０〜ステップＳ２６０で記憶部１３に保存された特徴量その他をデータベースに保存する。
【００４６】
以上の構成により、補償光学ＳＬＯ装置による取得される、網膜の同一位置を撮影した複数の平面画像に対して、重ね合わせに適した基準画像を選択することが可能になる。このような基準画像の選択により、画質のよい重ね合わせ画像を生成することが可能になる。
【００４７】
（第２の実施形態）
第１の実施形態では、被検眼の網膜の同一箇所を撮影した複数の平面画像から、重ね合わせ処理によって高い画質の視細胞像を取得するために、複数の平面画像から歪みの少ない画像を基準画像として選択する方法について説明した。しかし、眼底の固視微動はトレモア、ドリフト、フリック、サッケードといった特徴的な動きをすることが知られていて、例えばサッケードの前後では撮影箇所に大きな変化があるために、正確な位置合わせ処理を行うことが困難になる。本実施形態では、特徴量に基づいて平面画像群を固視微動の種類毎に分類することができる。これにより、トレモアのように変動が小さく、位置合わせ処理が比較的簡単な画像が多数存在するグループから基準画像を選択することができる。
【００４８】
平面画像群を取得した特徴量から、トレモアなどの小さな動きに対応する複数のグループに分類し、その分類に基づき基準画像及び重ね合わせに用いる画像を選択する。
【００４９】
図７のフローチャートを参照して、本実施形態の画像処理装置２０の処理手順を説明する。ここでステップＳ２１０、Ｓ２２０、Ｓ２３０、Ｓ２４０、Ｓ２６０、Ｓ２７０に関しては第１の実施形態で説明した処理手順と変わらないために、説明は省略する。
【００５０】
第１の実施形態では、補償光学ＳＬＯによる取得される平面画像群のそれぞれに対して画像の歪みを反映する特徴量を抽出し、それらを比較することで基準画像を選択していた。本実施形態においては、平面画像群全体を該特徴量に基づき分類することで、眼に特有の動きを考慮しその中で最良の重ね合わせ画像を生成することが可能な基準画像及び重ね合わせに用いる画像を選択する。
【００５１】
以下に各ステップについてより詳細に説明する。
【００５２】
＜ステップＳ８４０＞
ステップＳ８４０において、画像処理部１４に含まれる画像分類部（不図示）は、眼底の同一箇所を補償光学ＳＬＯにより撮影された複数の平面画像を、ステップＳ２４０で取得された画像の歪みを表す特徴量に基づき分類する。そうして分類した結果を、制御部１２を通じて記憶部１３に記憶する。
【００５３】
固視微動では、トレモアＴと呼ばれる小さな動き（数μｍ程度）とフリックＦ（数１０μｍ程度）と呼ばれる比較的大きな動きが組み合わされている。ここで、トレモア（Ｔ）中に撮影される平面画像には大きな歪みは生じないが、フリック（Ｆ）中の平面画像には大きな歪みが生じることになり、ステップＳ２４０で抽出した特徴量の変化は図８のようになる。
【００５４】
ここで、フリック（Ｆ）に対応する画像は歪みに対応する特徴量Ｄが大きくなるのに対して、トレモア（Ｔ）に対応した部分は歪みに対応する特徴量Ｄは小さくなる。よって、ステップＳ２４０で求めた特徴量に基づき、フリックにあたる平面画像を境界として、平面画像群を図８に示したように分類する。具体的には、ステップＳ２４０で特徴量Ｄの平均値と分散を求めたが、Ｄの値が３σよりも大きくなる場合をフリックに対応するとして分類の境界とした。
【００５５】
＜ステップＳ８５０＞
ステップＳ８５０において、画像処理部１４に含まれる画像選択部（不図示）は、ステップＳ８４０で求められた分類に基づき重ね合わせ画像の基準画像と、重ね合わせに用いる平面画像群を選択する。そうして選択した基準画像および重ね合わせに用いる平面画像を、制御部１２を通じて記憶部１３に記憶する。
【００５６】
高い画質の重ね合わせ結果を得るためには、なるべく近い同一箇所を撮影した画像が多くあることが望ましい。ステップＳ８４０の分類（Ａ〜Ｄ）において、同一グループに属する平面画像群は近い位置を撮影した歪みの少ない平面画像群と考えられるため、各グループに含まれるフレーム数が重要になる。さらにグループＡ〜Ｄは、分類は異なるが近い箇所を撮影している可能性があるので、各グループ間の距離を求める。具体的には、各グループで撮影時間的に中央に位置する平面画像を選択し、平面画像間の距離を位相限定相関法により求める。求められた距離が画像サイズの１５％程度以下の場合、つまり画像領域の８０％が重なる場合には、２つのグループは同一のグループであるとする。図８の場合には、グループＡとグループＤが同一グループに分類される。
【００５７】
その後、各グループで含まれる平面画像数が最も多いグループを選択する。重ね合わせ処理に用いる平面画像群は、選択されたグループに含まれる平面画像とする。さらに、それら平面画像の中で基準画像を選択するが、その方法は、ステップＳ２５０に示した場合と同様、特徴量Ｄが最も小さいものであってもよいし、別のコントラストなど別の特徴量に基づくものであってもよい。
【００５８】
（第３の実施形態）
第３の実施形態に係る補償光学ＳＬＯ装置の一例として、ＳＬＯ装置とＯＣＴ装置とを複合した複合装置について、図９を用いて説明する。この装置により、網膜の高横分解能の平面画像（ＳＬＯ像）と断層画像（ＯＣＴ像）との両方の撮像を行うことができる。ここで、ＳＬＯ装置は、空間光変調器を用いることにより、被検眼の前眼部で生じる収差を補正した平面画像を取得することができる。また、ＯＣＴ装置は、フーリエドメイン方式であり、被検眼の前眼部で生じる収差を補正した断層画像を取得することができる。
【００５９】
＜全体＞
図９を用いて、まず、本実施例における複合装置１００の概略構成について、具体的に説明する。光源１０１から出射した光は光カプラー１３１によって参照光１０５と測定光１０６とに分割される。測定光１０６は、シングルモードファイバー１３０−４、空間光変調器１５９、ＸＹスキャナ１１９、Ｘスキャナ１２１、球面ミラー１６０−１〜９等を介して、観察対象である被検眼１０７に導かれる。
【００６０】
被検眼１０７によって反射あるいは散乱された戻り光１０８となり、ディテクター１３８あるいはラインセンサ１３９に入射される。ディテクター１３８は戻り光１０８の光強度を電圧に変換し、その信号を用いて、被検眼１０７の平面画像が構成される。また、ラインセンサ１３９には参照光１０５と戻り光１０８とが合波されて入射され、被検眼１０７の断層画像が構成される。さらに、複数の取得された断層画像を用いて、血管の３次元的な走行が描出される。ここでは、波面収差を補正するデバイスとして空間光変調器を用いたが、波面収差を補正できればよく、可変形状ミラー等を用いることもできる。
【００６１】
＜光源＞
つぎに、光源１０１の周辺について説明する。光源１０１は代表的な低コヒーレント光源であるＳＬＤ（Ｓｕｐｅｒ  Ｌｕｍｉｎｅｓｃｅｎｔ  Ｄｉｏｄｅ）である。波長は８３０ｎｍ、バンド幅５０ｎｍである。ここでは、スペックルノイズの少ない平面画像を取得するために、低コヒーレント光源を選択している。また、光源の種類は、ここではＳＬＤを選択したが、低コヒーレント光が出射できればよく、ＡＳＥ（Ａｍｐｌｉｆｉｅｄ  Ｓｐｏｎｔａｎｅｏｕｓ  Ｅｍｉｓｓｉｏｎ）等も用いることができる。また、波長は眼を測定することを鑑みると、近赤外光が適する。さらに波長は、得られる平面画像の横方向の分解能に影響するため、なるべく短波長であることが望ましく、ここでは８３０ｎｍとする。観察対象の測定部位によっては、他の波長を選んでも良い。
【００６２】
光源１０１から出射された光は、シングルモードファイバー１３０−１と光カプラー１３１とを介して、参照光１０５と測定光１０６とに、９６：４の割合で分割される。１５３は偏光コントローラである。
【００６３】
＜参照光路＞
つぎに、参照光１０５の光路について説明する。光カプラー１３１にて分割された参照光１０５はシングルモードファイバー１３０−２を通して、レンズ１３５−１に導かれ、ビーム径４ｍｍの平行光になるよう、調整される。  次に、参照光１０５は、ミラー１５７−１〜４によって、参照ミラーであるミラー１１４に導かれる。参照光１０５の光路長は、測定光１０６の光路長と略同一に調整されているため、参照光１０５と測定光１０６とを干渉させることができる。次に、ミラー１１４にて反射され、再び光カプラー１３１に導かれる。ここで、参照光１０５が通過した分散補償用ガラス１１５は被検眼１０７に測定光１０６が往復した時の分散を、参照光１０５に対して補償するものである。ここでは、日本人の平均的な眼球の直径として代表的な値を想定し、Ｌ１＝２３ｍｍとする。
【００６４】
さらに、１１７−１は電動ステージであり、矢印で図示している方向に移動することができ、参照光１０５の光路長を、調整・制御することができる。また、電動ステージ１１７−１はパソコン１２５からドライバ部１８１内の電動ステージ駆動ドライバ１８３を介して制御される。
【００６５】
＜測定光路＞
つぎに、測定光１０６の光路について説明する。光カプラー１３１によって分割された測定光１０６はシングルモードファイバー１３０−４を介して、レンズ１３５−４に導かれ、ビーム径４ｍｍの平行光になるよう調整される。また、偏光コントローラ１５３−１又は２は、測定光１０６の偏光状態を調整することができる。ここでは、測定光１０６の偏光状態は紙面に平行な方向の直線偏光に調整されている。
【００６６】
測定光１０６は、ビームスプリッタ１５８、可動式ビームスプリッタ１６１（分岐手段とも呼ぶ。）を通過し、球面ミラー１６０−１、１６０−２を介し、空間光変調器１５９にて入射・変調される。ここで、空間光変調器１５９は、液晶の配向性を利用して変調を行う変調器であり、紙面に平行な方向の直線偏光（Ｐ偏光）の位相を変調する向きに配置され、測定光１０６の偏光の向きと合わせている。
【００６７】
さらに、測定光１０６は偏光板１７３を通過し、球面ミラー１６０−３、１６０−４を介し、Ｘスキャナ１２１のミラーに入射される。ここで、偏光板１７３は戻り光１０８のうち紙面に平行な方向の直線偏光のみを空間光変調器１５９に導く役割がある。また、ここで、Ｘスキャナ１２１は測定光１０６を紙面に平行な方向に走査するＸスキャナであり、ここでは共振型スキャナを用いている。駆動周波数は約７．９ｋＨｚである。
【００６８】
さらに、測定光１０６は球面ミラー１６０−５〜６を介し、ＸＹスキャナ１１９のミラーに入射される。ここで、ＸＹスキャナ１１９は一つのミラーとして記したが、実際にはＸスキャン用ミラーとＹスキャン用ミラーとの２枚のミラーが近接して配置されるものである。また、測定光１０６の中心はＸＹスキャナ１１９のミラーの回転中心と一致するように調整されている。ＸＹスキャナ１１９の駆動周波数は〜５００Ｈｚの範囲で可変できる。球面ミラー１６０−７〜９は網膜１２７を走査するための光学系であり、測定光１０６を角膜１２６の付近を支点として、網膜１２７をスキャンする役割がある。
【００６９】
ここで、測定光１０６のビーム径は４ｍｍであるが、より高分解能な断層画像を取得するために、ビーム径はより大径化してもよい。また、１１７−２は電動ステージであり、矢印で図示している方向に移動することができ、付随する球面ミラーである球面ミラー１６０−８の位置を、調整・制御することができる。電動ステージ１１７−２は電動ステージ１１７−１と同様に、電動ステージ駆動ドライバ１８３によって制御される。球面ミラー１６０−８の位置を調整することで、被検眼１０７の網膜１２７の所定の層に測定光１０６を合焦し、観察することが可能になる。初期状態では、測定光１０６は平行光の状態で、角膜１２６に入射するように、球面ミラー１６０−８の位置が調整されている。
【００７０】
また、被検眼１０７が屈折異常を有している場合にも対応できる。測定光１０６は被検眼１０７に入射すると、網膜１２７からの反射や散乱により戻り光１０８となり、再び光カプラー１３１に導かれ、ラインカメラ１３９に到達する。また、戻り光１０８の一部は可動式ビームスプリッタ１６１で反射され、レンズ１３５−５を介して、ディテクター１３８に導かれる。ここで、１７２はピンホールを有する遮光板であり、戻り光１０８の内、網膜１２７に合焦していない不要な光を遮断する役割がある。また、遮光板１７２はレンズ１３５−５の合焦位置に共役に配置される。また、遮光板１７２のピンホールの直径は例えば５０μｍである。ディテクター１３８は例えば高速・高感度な光センサであるＡＰＤ（Ａｖａｌａｎｃｈｅ  Ｐｈｏｔｏ  Ｄｉｏｄｅ）が用いられる。また、ビームスプリッタ１５８にて分割される戻り光１０８の一部は、波面センサ１５５に入射される。波面センサ１５５はシャックハルトマン方式の波面センサである。
【００７１】
ここで、ＸＹスキャナ１１９、Ｘスキャナ１２１、角膜１２６、波面センサ１５５、空間光変調器１５９は光学的に共役になるよう、球面ミラー１６０−１〜９が配置されている。そのため、波面センサ１５５は被検眼１０７の収差を測定することが可能になっている。また、空間光変調器１５９は被検眼１０７の収差を補正することが可能になっている。さらに、得られた収差に基づいて、空間光変調器１５９をリアルタイムに制御することで、被検眼１０７で発生する収差を補正し、より高横分解能な断層画像の取得を可能にしている。
【００７２】
＜測定系の構成＞
つぎに、測定系の構成について説明する。複合装置１００は、断層画像（ＯＣＴ像）及び平面画像（ＳＬＯ像）を取得することができる。まず、断層画像の測定系について説明する。戻り光１０８は光カプラー１３１によって合波される。合波された光１４２は、シングルモードファイバー１３０−３、レンズ１３５−２を介して、透過型グレーティング１４１に導かれ、波長毎に分光され、レンズ１３５−３を介してラインカメラ１３９に入射される。ラインカメラ１３９は位置（波長）毎に光強度を電圧に変換し、その電圧信号はフレームグラバー１４０にてデジタル値に変換されて、パソコン１２５にて、被検眼１０７の断層画像が構成される。ここでは、ラインカメラ１３９は１０２４画素を有し、合波された光１４２の波長毎（１０２４分割）の強度を得ることができる。
【００７３】
つぎに、平面画像の測定系について説明する。戻り光１０８の一部は、可動式ビームスプリッタ１６１で反射される。反射された光は遮光板１７２によって不要な光が遮断された後、ディテクター１３８に到達し、光の強度が電気信号に変換される。得られた電気信号に対して、パソコン１２５にてＸスキャナ１２１とＸＹスキャナ１１９との走査信号と同期したデータ処理が行われ、平面画像が形成される。ビームスプリッタ１５８にて分割される戻り光１０８の一部は、波面センサ１５５に入射され、戻り光１０８の収差が測定される。波面センサ１５５にて得られた画像信号は、パソコン１２５に取り込まれ、収差が算出される。得られた収差はツェルニケ多項式を用いて表現され、これは被検眼１０７の収差を示している。ここで、ツェルニケ多項式は、チルト（傾き）の項、デフォーカス（ｄｅｆｏｃｕｓ）の項、アスティグマ（非点収差）の項、コマの項、トリフォイルの項等からなる。
【００７４】
＜ＯＣＴ像の取得方法＞
つぎに、複合装置１００を用いた断層画像（ＯＣＴ像）の取得方法について図１０を用いて説明する。ＯＣＴ装置１００は、ＸＹスキャナ１１９を制御し、Ｘスキャナ１２１を固定ミラーとして用いて、ラインカメラ１３９で干渉縞を取得することで、網膜１２７の断層画像を取得することができる。戻り光１０８がディテクター１３８に導光されないように可動式ビームスプリッタ１６１を制御する。また、Ｘスキャナ１２１、ＸＹスキャナ１１９は、パソコン１２５からドライバ部１８１内の光スキャナ駆動ドライバ１８２を介して制御される。ここでは、網膜１２７の断層画像（光軸に平行な面）の取得方法について説明する。図１０（ａ）は被検眼１０７の模式図であり、複合装置１００によって観察されている様子を示している。図１０（ａ）に示すように、測定光１０６は角膜１２６を通して、網膜１２７に入射すると様々な位置における反射や散乱により戻り光１０８となり、それぞれの位置での時間遅延を伴って、ラインカメラ１３９に到達する。ここでは、光源１０１のバンド幅が広く、コヒーレンス長が短いために、参照光路の光路長と測定光路の光路長とが略等しい場合に、ラインカメラ１３９にて、干渉縞が検出できる。上述のように、ラインカメラ１３９で取得されるのは波長軸上のスペクトル領域の干渉縞となる。
【００７５】
次に、波長軸上の情報である干渉縞を、ラインカメラ１３９と透過型グレーティング１４１との特性を考慮して、光周波数軸の干渉縞に変換する。さらに、変換された光周波数軸の干渉縞を逆フーリエ変換することで、深さ方向の情報が得られる。さらに、図１０（ｂ）に示すように、ＸＹスキャナ１１９を駆動しながら、干渉縞を検知すれば、各Ｘ軸の位置毎に干渉縞が得られ、つまり、各Ｘ軸の位置毎の深さ方向の情報を得ることができる。結果として、ＸＺ面での戻り光１０８の強度の２次元分布が得られ、それはすなわち断層画像１３２である（図１０（ｃ））。本来は、断層画像１３２は上記説明したように、該戻り光１０８の強度をアレイ状に並べたものであり、例えば該強度をグレースケールに当てはめて、表示されるものである。Ｘ方向の長さは７００μｍであり、下に説明するＳＬＯ像と同様である。ここでは得られた断層画像の境界のみ強調して表示している。ここで、１４６は網膜色素上皮層、１４７は視神経線維層である、１７８は血管である。また、各Ｙ軸の位置毎に複数の断層画像を取得すれば、血管の３次元的な走行を描出できる。
【００７６】
＜ＳＬＯ像の取得方法＞
つぎに、複合装置１００を用いた平面画像（ＳＬＯ像）の取得方法について、図１１を用いて説明する。複合装置１００は、ＸＹスキャナ１１９のＹ軸方向のみとＸスキャナ１２１とを動作・制御し、ＸＹスキャナ１１９のＸ軸方向を固定し、ディテクター１３８で戻り光１０８の強度を取得することで、網膜１２７の平面画像を取得することができる。Ｘスキャナ１２１とＸＹスキャナ１１９は、パソコン１２５からドライバ部１８１内の光スキャナ駆動ドライバ１８２を介して制御される。また、複合装置１００は、波面センサ１５５で測定した被検眼１０７の収差を用いて空間光変調器１５９を制御し、被検眼１０７等にて生じる収差を補正しながら平面画像を取得することができる。また、空間光変調器１５９をリアルタイムに制御しながら平面画像を取得することができる。
【００７７】
図１１（ａ）は被検眼１０７の模式図であり、ＳＬＯ装置１０９によって観察されている様子を示している。図１１（ａ）に示すように、測定光１０６は角膜１２６を通して、網膜１２７に入射すると様々な位置における反射や散乱により戻り光１０８となり、ディテクター１３８に到達する。さらに、図１１（ｂ）に示すように、ＸＹスキャナ１０９をＸ方向に駆動しながら、戻り光１０８の強度を検知すれば、各Ｘ軸の位置毎の情報を得ることができる。
【００７８】
さらに、図１１（ｃ）に示すように、ＸＹスキャナ１０９のＸ軸とＹ軸とを同時に駆動し、網膜１２７のある撮像範囲１９２に対して、測定光１０６を軌跡１９３のように、ラスタースキャンしながら戻り光１０８の強度を検知すれば、戻り光１０８の強度の２次元分布が得られ、それはすなわち平面画像１７７（図１１（ｄ））である。
【００７９】
ここで、測定光１０６は右上の点Ｓから、左下の点Ｅに向かってスキャンされ、その間の戻り光１０８の強度が平面画像１７７の構成に用いられる。点Ｅから点Ｓへの軌跡１９３は次の平面画像１７７の撮像の撮像のための準備である。スキャンにかかる時間は、図中の軌跡１９３に対して、Ｓ→Ｅが８４％、Ｅ→Ｓが１６％であり、この比は上述のＹスキャナの駆動波形のデューティ比に基づいている。また、図１１（ｃ）では簡単のため、軌跡１９３におけるＸ方向のスキャン回数を少なめに記している。ここで、この平面画像１７７の大きさは７００×３５０μｍ、取得に要する時間は約１５．６ｍｓである。この時間はＹスキャナの駆動周波数に基づいている。
【００８０】
また、平面画像１７７中には、戻り光１０８の強度が比較的大きい視細胞群１７９が明るく、比較的小さい血管１７８が暗く描出される。また、血管１７８に血球（不図示）が明るく描出される。また、平面画像１７７を時間的に連続して取得すれば、血球が血管１７８内を移動する様子を描出することができる。また、連続して取得した平面画像１７７から、血球が描出されている血管１７８を抽出し、撮影順に積み重ねて時空間画像を作成してもよい。
【００８１】
（その他の実施形態）
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【図１】
0005847510.tif 000005
【図２】
0005847510.tif 000006
【図３】
0005847510.tif 000007
【図７】
0005847510.tif 000008
【図８】
0005847510.tif 000009
【図９】
0005847510.tif 000010
【図１０】
0005847510.tif 000011
【図１１】
0005847510.tif 000012
【図４】
0005847510.tif 000013
【図５】
0005847510.tif 000014
【図６】
0005847510.tif 000015