WO2022176828A1

WO2022176828A1 - 眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラム

Info

Publication number: WO2022176828A1
Application number: PCT/JP2022/005828
Authority: WO
Inventors: 聡大三野
Original assignee: 株式会社トプコン
Priority date: 2021-02-19
Filing date: 2022-02-15
Publication date: 2022-08-25
Also published as: US20220265136A1; US11974806B2; JPWO2022176828A1

Abstract

眼科情報処理装置は、生成部と、抽出部と、位相補正部とを含む。生成部は、被検眼の複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化処理を行うことにより位相差プロファイルを生成する。抽出部は、生成部により生成された位相差プロファイルから位相ドリフトを抽出する。位相補正部は、抽出部により抽出された位相ドリフトに基づいて、上記の２つのＢフレームのいずれか一方の第１複素ＯＣＴデータの位相を補正する。

Description

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

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

　ＣＡＯ（Ｃｏｍｐｕｔａｔｉｏｎａｌ　Ａｄａｐｔｉｖｅ　Ｏｐｔｉｃｓ）は、被測定物や光学系に起因した収差を演算処理で補正する技術である。ＣＡＯは、波面センサ等の公知のハードウェアを用いて取得された収差情報に基づいて測定データを補正することで測定精度を向上させることが可能になる。

　例えば、非特許文献１、及び非特許文献２には、このようなＣＡＯを光コヒーレンストモグラフィ（Ｏｐｔｉｃａｌ　Ｃｏｈｅｒｅｃｅ：Ｔｏｍｏｇｒａｐｈｙ：ＯＣＴ）システムに適用した手法が開示されている。

Ｎａｔｈａｎ　Ｄ．　Ｓｈｅｍｏｎｓｋｉ　ｅｔ　ａｌ．，　"Ｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌ　ｍｏｔｉｏｎ　ｃｏｒｒｅｃｔｉｏｎ　ｕｓｉｎｇ　ｓｐｅｃｋｌｅ　ａｎｄ　ｐｈａｓｅ　ｆｏｒ　ｉｎ　ｖｉｖｏ　ｃｏｍｐｕｔｅｄ　ｏｐｔｉｃａｌ　ｉｎｔｅｒｆｅｒｏｍｅｔｒｉｃ　ｔｏｍｏｇｒａｐｈｙ"，　ＢＩＯＭＥＤＩＣＡＬ　ＯＰＴＩＣＳ　ＥＸＰＲＥＳＳ，　４　ＮＯＶ　２０１４，　Ｖｏｌ．　５，　Ｎｏ．　１２，　ｐｐ．４１３１－４１４３Ｋｅｎｓｕｋｅ　Ｏｉｋａｗａ　ｅｔ　ａｌ．，　"Ｂｕｌｋ　ｐｈａｓｅ　ｅｒｒｏｒ　ｃｏｒｒｅｃｔｉｏｎ　ｆｏｒ　ｈｏｌｏｇｒａｐｈｉｃ　ｓｉｇｎａｌ　ｐｒｏｃｅｓｓｉｎｇ　ｏｆ　ｏｐｔｉｃａｌ　ｃｏｈｅｒｅｎｃｅ　ｔｏｍｏｇｒａｐｈｙ"，　Ｐｒｏｃ．ＳＰＩＥ１１５２１，　Ｂｉｏｍｅｄｉｃａｌ　Ｉｍａｇｉｎｇ　ａｎｄ　Ｓｅｎｓｉｎｇ　Ｃｏｎｆｅｒｅｎｃｅ　２０２０，　１１５２１０Ｐ（１５　Ｊｕｎｅ　２０２０）；ｄｏｉ：１０．１１１７／１２．２５７３２３１

　ＯＣＴにより取得された測定データに含まれる位相情報は、ＯＣＴアンギオグラフィ（ＯＣＴ　Ａｎｇｉｏｇｒａｐｈｙ：ＯＣＴＡ）や、ＯＣＴにＣＡＯが適用されたＣＡＯ－ＯＣＴ等のアプリケーションにおいて有用である。この位相情報は、観察対象の動きや測定系内の位相不安定性に起因した位相ドリフトの影響を受ける。そのため、位相ドリフトの影響を受けないように位相を安定化させることにより、ＯＣＴＡやＣＡＯ－ＯＣＴの精度を向上させることができる。

　例えば、隣接する位置の測定データ間で位相差を算出し、算出された位相差を相殺することで、位相を安定化させることが可能である。これは、被検眼の網膜のような無数の散乱体からなる被測定物に対して位相差の期待値が零になると仮定できることを前提としている。しかしながら、現実的にはこの想定からの逸脱／誤差が存在し、上記の方法で位相補正を繰り返すにつれて位相誤差が積算され、積算された誤差の影響が無視できなくなる場合がある。

　本発明は、このような事情を鑑みてなされたものであり、その目的の１つは、ＯＣＴにより取得された測定データに含まれる位相を安定化させる新たな技術を提供することにある。

　実施形態の第１態様は、被検眼の複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化処理を行うことにより位相差プロファイルを生成する生成部と、前記生成部により生成された前記位相差プロファイルから位相ドリフトを抽出する抽出部と、前記抽出部により抽出された前記位相ドリフトに基づいて、前記２つのＢフレームのいずれか一方の第１複素ＯＣＴデータの位相を補正する位相補正部と、を含む、眼科情報処理装置である。

　実施形態の第２態様は、第１態様において、前記位相差プロファイルから信頼性の低い位相差データを除去する除去部を含み、前記抽出部は、前記除去部により前記信頼性の低い位相差データが除去された前記位相差プロファイルから局所的な位相差変化が除去された前記位相ドリフトを抽出する。

　実施形態の第３態様では、第２態様において、前記信頼性の低い位相差データは、信号強度が所定の閾値レベル未満であるデータ、及び前記平均化処理において所定の深さ位置における位相差による寄与が所定の閾値を超えるデータの少なくとも１つを含む。

　実施形態の第４態様では、第１態様～第３態様のいずれかにおいて、前記抽出部は、前記位相差プロファイルに対して、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つを施すことにより前記位相ドリフトを抽出する。

　実施形態の第５態様では、第４態様において、前記スムージング処理、前記多項式フィッティング、及び前記高周波カットフィルタ処理の少なくとも１つのパラメータは、前記複素ＯＣＴデータを取得する測定系の位相変化の周波数、前記位相変化の速度、及び前記測定系におけるＯＣＴスキャンに要する時間の少なくとも１つに基づいて設定される。

　実施形態の第６態様では、第１態様～第５態様のいずれかにおいて、隣接する第１のＢフレーム及び第２のＢフレームについて、前記生成部により生成された位相差プロファイルに基づいて前記第２のＢフレームの複素ＯＣＴデータの位相を補正した後、前記第２のＢフレーム及び前記第２のＢフレームに隣接する第３のＢフレームについて、前記生成部により生成された位相差プロファイルに基づいて前記第３のＢフレームの複素ＯＣＴデータの位相を補正する。

　実施形態の第７態様は、第１態様～第６態様のいずれかにおいて、前記位相補正部により少なくとも１つのＢフレームの複素ＯＣＴデータが補正された前記被検眼の複素ＯＣＴデータに基づいて前記被検眼のＯＣＴ画像を形成する画像形成部を含む。

　実施形態の第８態様は、前記被検眼に対して光コヒーレンストモグラフィを実行することにより複素ＯＣＴデータを取得する測定系と、第１態様～第７態様のいずれかの眼科情報処理装置と、を含む、眼科装置である。

　実施形態の第９態様は、被検眼の複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化処理を行うことにより位相差プロファイルを生成する生成ステップと、前記生成ステップにおいて生成された前記位相差プロファイルから位相ドリフトを抽出する抽出ステップと、前記抽出ステップにおいて抽出された前記位相ドリフトに基づいて、前記２つのＢフレームのいずれか一方の第１複素ＯＣＴデータの位相を補正する位相補正ステップと、を含む、眼科情報処理方法である。

　実施形態の第１０態様は、第９態様において、前記位相差プロファイルから信頼性の低い位相差データを除去する除去ステップを含み、前記抽出ステップは、前記除去ステップにおいて前記信頼性の低い位相差データが除去された前記位相差プロファイルから局所的な位相差変化が除去された前記位相ドリフトを抽出する。

　実施形態の第１１態様では、第１０態様において、前記信頼性の低い位相差データは、信号強度が所定の閾値レベル未満であるデータ、及び前記平均化処理において所定の深さ位置における位相差による寄与が所定の閾値を超えるデータの少なくとも１つを含む。

　実施形態の第１２態様では、第９態様～第１１態様のいずれかにおいて、前記抽出ステップは、前記位相差プロファイルに対して、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つを施すことにより前記位相ドリフトを抽出する。

　実施形態の第１３態様では、第１２態様において、前記スムージング処理、前記多項式フィッティング、及び前記高周波カットフィルタ処理の少なくとも１つのパラメータは、前記複素ＯＣＴデータを取得する測定系の位相変化の周波数、前記位相変化の速度、及び前記測定系におけるＯＣＴスキャンに要する時間の少なくとも１つに基づいて設定される。

　実施形態の第１４態様では、第９態様～第１３態様のいずれかにおいて、隣接する第１のＢフレーム及び第２のＢフレームについて、前記生成ステップにおいて生成された位相差プロファイルに基づいて前記第２のＢフレームの複素ＯＣＴデータの位相を補正した後、前記第２のＢフレーム及び前記第２のＢフレームに隣接する第３のＢフレームについて、前記生成ステップにおいて生成された位相差プロファイルに基づいて前記第３のＢフレームの複素ＯＣＴデータの位相を補正する。

　実施形態の第１５態様は、第９態様～第１４態様のいずれかにおいて、前記位相補正ステップにおいて少なくとも１つのＢフレームの複素ＯＣＴデータが補正された前記被検眼の複素ＯＣＴデータに基づいて前記被検眼のＯＣＴ画像を形成する画像形成ステップを含む。

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

　なお、上記した複数の態様に係る構成を任意に組み合わせることが可能である。

　本発明によれば、ＯＣＴにより取得された測定データに含まれる位相を安定化させる新たな技術を提供することができる。

実施形態に係る眼科システムの構成の一例を表す概略図である。実施形態に係る眼科システムの構成の一例を表す概略図である。実施形態に係る眼科システムの構成の一例を表す概略図である。実施形態に係る眼科装置の光学系の構成の一例を表す概略図である。実施形態に係る眼科装置の光学系の構成の一例を表す概略図である。実施形態に係る眼科装置の処理系の構成の一例を表す概略図である。実施形態に係る眼科装置が実行する処理を説明するための概略図である。実施形態に係る眼科装置が実行する処理を説明するための概略図である。実施形態に係る眼科装置が実行する処理を説明するための概略図である。実施形態に係る眼科装置の処理系の構成の一例を表す概略図である。実施形態に係る眼科装置の処理系の構成の一例を表す概略図である。実施形態に係る眼科装置が実行する処理を説明するための概略図である。実施形態に係る眼科装置が実行する処理を説明するための概略図である。実施形態に係る眼科装置が実行する処理を説明するための概略図である。実施形態に係る眼科装置の処理系の構成の一例を表す概略図である。実施形態に係る眼科装置の動作例を表すフローチャートである。実施形態に係る眼科装置の動作例を表すフローチャートである。実施形態に係る眼科装置の動作例を表すフローチャートである。実施形態に係る眼科装置が実行する処理を説明するための概略図である。実施形態に係る眼科装置の動作例を表すフローチャートである。実施形態に係る眼科装置の動作例を表すフローチャートである。実施形態に係る眼科装置の動作例を表すフローチャートである。実施形態に係る眼科装置の動作例を表すフローチャートである。実施形態に係る眼科装置の動作例を表すフローチャートである。

　この発明に係る眼科情報処理装置、眼科装置、眼科情報処理方法、プログラム、及び記録媒体の実施形態の例について、図面を参照しながら詳細に説明する。なお、この明細書において引用された文献の記載内容や任意の公知技術を、以下の実施形態に援用することが可能である。

　実施形態に係る眼科情報処理装置は、外部に設けられた眼科装置（ＯＣＴ装置）からＯＣＴデータ（測定データ、複素ＯＣＴデータ）を取得する。外部に設けられた眼科装置は、被検眼に対してＯＣＴを実行することにより複素ＯＣＴデータを取得する。眼科情報処理装置は、被検眼の複素ＯＣＴデータに対して位相安定化処理を施すことが可能である。具体的には、眼科情報処理装置は、複素ＯＣＴデータからＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化された位相差プロファイルを生成し、生成された位相差プロファイルから位相ドリフトを抽出する。眼科情報処理装置は、抽出された位相ドリフトに基づいて複素ＯＣＴデータの位相を補正する。このとき、眼科情報処理装置は、位相差プロファイルから信頼性の低い位相差データを除去し、信頼性の低い位相差データが除去された位相差プロファイルから位相ドリフトを抽出することが可能である。これにより、ノイズや構造由来の位相変化の影響を受けることなく、位相安定化処理が施されたＯＣＴデータが得られる。その結果、被検眼の動きやＯＣＴ測定系内の位相不安定性に起因した位相変化や位相ノイズにロバストなＯＣＴデータが得られる。

　例えば、眼科情報処理装置は、上記のように位相が補正されたＯＣＴデータに対して被検眼や光学系に起因した収差を補正することが可能である。具体的には、眼科情報処理装置は、収差を補正するための最適なＣＡＯフィルタ（ＣＡＯフィルタのパラメータ）を探索し、探索されたＣＡＯフィルタを用いてＯＣＴデータに対して収差補正処理を実行する。例えば、外部に設けられた眼科装置（屈折力測定装置）から被検眼の屈折力を表す屈折力情報を取得可能に構成される場合、眼科情報処理装置は、取得された屈折力情報に基づいてＣＡＯフィルタのパラメータ（フィルタ情報）を算出することが可能である。眼科情報処理装置は、探索されたＣＡＯフィルタ（パラメータ）に対応した収差情報（具体的には、被検眼の瞳面における収差情報）を表示手段に表示させることが可能である。

　実施形態に係る眼科装置は、ＯＣＴ計測機能及び屈折力測定機能の少なくとも一方の機能に加えて、実施形態に係る眼科情報処理装置の機能を備える。

　実施形態に係る眼科情報処理方法は、実施形態に係る眼科情報処理装置においてプロセッサ（コンピュータ）により実行される処理を実現するための１以上のステップを含む。実施形態に係るプログラムは、プロセッサに、実施形態に係る眼科情報処理方法の各ステップを実行させる。実施形態に係る記録媒体は、実施形態に係るプログラムが記録されたコンピュータ読み取り可能な非一時的な記録媒体（記憶媒体）である。

　本明細書において「プロセッサ」は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ　Ｓｐｅｃｉｆｉｃ　Ｉｎｔｅｇｒａｔｅｄ　Ｃｉｒｃｕｉｔ）、プログラマブル論理デバイス（例えば、ＳＰＬＤ（Ｓｉｍｐｌｅ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｌｏｇｉｃ　Ｄｅｖｉｃｅ）、ＣＰＬＤ（Ｃｏｍｐｌｅｘ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｌｏｇｉｃ　Ｄｅｖｉｃｅ）、ＦＰＧＡ（Ｆｉｅｌｄ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｇａｔｅ　Ａｒｒａｙ））等の回路を意味する。プロセッサは、例えば、記憶回路や記憶装置に格納されているプログラムを読み出し実行することで、実施形態に係る機能を実現する。

　以下、実施形態では、ＯＣＴを用いた計測（又は撮影）においてスウェプトソースタイプのＯＣＴの手法を用いる場合について特に詳しく説明する。しかしながら、他のタイプ（例えば、スペクトラルドメインタイプ）のＯＣＴを用いた計測を行う眼科装置に対して、実施形態に係る構成を適用することも可能である。

　この明細書では、ＯＣＴによって取得される画像をＯＣＴ画像と総称することがある。また、ＯＣＴ画像を形成するための計測動作をＯＣＴ計測と呼ぶことがある。

　いくつかの実施形態に係る眼科装置は、眼科撮影装置と、眼科測定装置と、眼科治療装置とのうちのいずれか１つ以上を含む。いくつかの実施形態の眼科装置に含まれる眼科撮影装置は、例えば、眼底カメラ、走査型レーザー検眼鏡、スリットランプ検眼鏡、手術用顕微鏡等のうちのいずれか１つ以上である。また、いくつかの実施形態の眼科装置に含まれる眼科測定装置は、例えば、眼屈折検査装置、眼圧計、スペキュラーマイクロスコープ、ウェーブフロントアナライザ、視野計、マイクロペリメータ等のうちのいずれか１つ以上である。また、いくつかの実施形態の眼科装置に含まれる眼科治療装置は、例えば、レーザー治療装置、手術装置、手術用顕微鏡等のうちのいずれか１つ以上である。

＜眼科システム＞
　実施形態に係る眼科システムは、実施形態に係る眼科情報処理装置の機能と、実施形態に係る眼科装置の機能とを備える。

　図１に、実施形態に係る眼科システムの第１構成例のブロック図を示す。

　実施形態の第１構成例に係る眼科システム１０００は、眼科装置としてのＯＣＴ装置１１００と、眼科装置としての屈折力測定装置１１１０と、眼科情報処理装置１２００と、表示部１３００と、操作部１４００とを含む。眼科情報処理装置１２００は、表示部１３００及び操作部１４００の少なくとも一方の機能を備えていてもよい。

　眼科情報処理装置１２００とＯＣＴ装置１１００とは、有線又は無線の通信路を介して通信可能に接続されている。また、眼科情報処理装置１２００と屈折力測定装置１１１０とは、有線又は無線の通信路を介して通信可能に接続されている。

　ＯＣＴ装置１１００は、被検眼に対してＯＣＴ（ＯＣＴ計測）を実行するＯＣＴ光学系１１０１を含み、ＯＣＴを実行することにより得られた複素ＯＣＴデータを眼科情報処理装置１２００に送信する。例えば、ＯＣＴ光学系１１０１は、特開２０１９－１５４９８５号公報に開示されているような公知の構成を有していてよい。

　屈折力測定装置１１１０は、被検眼の屈折力を測定する屈折力測定光学系１１１１を含み、屈折力測定により得られた被検眼の屈折力を表す屈折力情報を眼科情報処理装置１２００に送信する。例えば、屈折力測定光学系１１１１は、特開２０１９－１５４９８５号公報に開示されているような公知の構成を有していてよい。

　眼科情報処理装置１２００は、ＯＣＴ装置１１００からの被検眼の複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化された位相差プロファイルを生成する。眼科情報処理装置１２００は、生成された位相差プロファイルから位相ドリフトを抽出し、抽出された位相ドリフトに基づいて、上記の２つのＢフレームのいずれか一方の複素ＯＣＴデータの位相を補正する。

　更に、眼科情報処理装置１２００は、位相が補正された複素ＯＣＴデータに対して収差を補正するための最適なＣＡＯフィルタ（ＣＡＯフィルタのパラメータ）を探索し、探索されたＣＡＯフィルタを用いて複素ＯＣＴデータの収差補正処理を実行する。眼科情報処理装置１２００は、屈折力測定装置１１１０からの被検眼の屈折力情報（屈折力の測定結果）を用いてＣＡＯフィルタのパラメータを算出することで、上記のＣＡＯフィルタの探索処理の簡略化又は省略化を図ることが可能である。

　眼科情報処理装置１２００は、探索されたＣＡＯフィルタを用いて収差補正が行われた複素ＯＣＴデータに基づいて形成されたＯＣＴ画像や、探索されたＣＡＯフィルタに対応した被検眼の瞳面における収差情報を表示部１３００に表示させることが可能である。

　眼科情報処理装置１２００は、操作部１４００からのユーザの操作内容に対応する操作情報に基づいて表示部１３００を制御したり、装置内部で実行される処理を制御したりすることが可能である。例えば、眼科情報処理装置１２００は、操作部１４００からの操作情報に基づいてＣＡＯフィルタのパラメータを変更し、パラメータが変更されたＣＡＯフィルタを用いて複素ＯＣＴデータに対して収差補正処理を実行する。眼科情報処理装置１２００は、収差が補正された複素ＯＣＴデータに基づいてＯＣＴ画像を形成する。例えば、眼科情報処理装置１２００は、操作部１４００からの操作情報に基づいてＣＡＯフィルタのパラメータを変更し、パラメータが変更されたＣＡＯフィルタに対応した被検眼の瞳面における収差情報を表示部１３００に表示させる。

　図１において、眼科システム１０００は、眼科情報処理装置１２００がＯＣＴ計測機能と屈折力測定機能とを有する眼科装置からＯＣＴデータと屈折力情報とを取得するように構成されていてもよい。

　いくつかの実施形態では、眼科装置は、実施形態に係る眼科情報処理装置の機能を含む。

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

　実施形態の第２構成例に係る眼科システム１０００ａは、第１構成例に係る眼科システム１０００の機能を単体で実現する眼科装置１５００ａを含む。眼科装置１５００ａは、眼科情報処理部１２００ａと、表示部１３００と、操作部１４００と、測定光学系１６００ａとを含む。測定光学系１６００ａは、ＯＣＴ光学系１１０１と、屈折力測定光学系１１１１とを含む。眼科情報処理部１２００ａは、図１の眼科情報処理装置１２００の機能を有し、眼科装置１５００ａの各部を制御する。

　図２において、ＯＣＴ計測機能と実施形態に係る眼科情報処理装置の機能とを有する眼科装置が、屈折力測定機能を有する眼科装置から屈折力情報を取得するように構成されていてもよい。

　図３に、実施形態に係る眼科システムの第３構成例のブロック図を示す。なお、図３において、図１又は図２と同様の部分には同一符号を付し、適宜説明を省略する。

　実施形態の第３構成例に係る眼科システム１０００ｂは、眼科装置１５００ｂと、屈折力測定装置１７００ｂとを含む。眼科装置１５００ｂと屈折力測定装置１７００ｂとは、有線又は無線の通信路を介して通信可能に接続されている。

　眼科装置１５００ｂは、眼科情報処理部１２００ｂと、表示部１３００と、操作部１４００と、測定光学系１６００ｂとを含む。測定光学系１６００ｂは、ＯＣＴ光学系１１０１を含む。眼科情報処理部１２００ｂは、眼科装置１５００ｂの各部を制御する。

　屈折力測定装置１７００ｂは、屈折力測定光学系１１１１を含み、屈折力測定により得られた被検眼の屈折力を表す屈折力情報を眼科装置１５００ｂに送信する。

　図３において、屈折力測定機能と実施形態に係る眼科情報処理装置の機能とを有する眼科装置が、ＯＣＴ計測機能を有する眼科装置からＯＣＴデータを取得するように構成されていてもよい。

　以下、第３構成例に係る眼科装置を例に、実施形態の構成及び動作を説明する。

　以下の実施形態に係る眼科装置は、ＯＣＴ計測が可能なＯＣＴ装置と眼底カメラとを含む。しかしながら、以下の実施形態に係る構成を、単体のＯＣＴ装置に組み込むことも可能である。

　また、以下では、主に被検眼の眼底に対するＯＣＴ計測が可能な眼科装置を例に説明するが、実施形態に係る眼科装置は、被検眼の前眼部に対してＯＣＴ計測が可能であってよい。いくつかの実施形態では、測定光の焦点位置を変更するレンズを移動することで、ＯＣＴ計測の範囲や計測部位を変更する。いくつかの実施形態では、１以上のアタッチメント（対物レンズ、前置レンズ等）を加えることで、眼底に対するＯＣＴ計測と、前眼部に対するＯＣＴ計測と、眼底及び前眼部を含む全眼球に対するＯＣＴ計測とが可能な構成を有する。いくつかの実施形態では、眼底計測用の眼科装置において、対物レンズと被検眼との間に前置レンズを配置することで平行光束にされた測定光を被検眼に入射させることで前眼部に対するＯＣＴ計測を行う。

＜眼科装置＞
［構成］
（光学系）
　図４に、実施形態に係る眼科装置１５００ｂの光学系の構成例を示す。以下、対物レンズ２２の光軸の方向（後述の測定光ＬＳの進行方向）をｚ方向とし、ｚ方向に直交する水平方向をｘ方向とし、ｚ方向に直交する垂直方向をｙ方向とする。

　眼科装置１５００ｂは、眼底カメラユニット２、ＯＣＴユニット１００及び演算制御ユニット２００を含む。眼底カメラユニット２には、被検眼Ｅの正面画像を取得するための光学系や機構が設けられている。ＯＣＴユニット１００には、ＯＣＴを実行するための光学系や機構の一部が設けられている。ＯＣＴを実行するための光学系や機構の他の一部は、眼底カメラユニット２に設けられている。演算制御ユニット２００は、各種の演算や制御を実行する１以上のプロセッサを含む。これらに加え、被検者の顔を支持するための部材（顎受け、額当て等）や、ＯＣＴの対象部位を切り替えるためのレンズユニット（例えば、前眼部ＯＣＴ用アタッチメント）等の任意の要素やユニットが眼科装置１５００ｂに設けられてもよい。いくつかの実施形態では、レンズユニットが手動で被検眼Ｅと後述の対物レンズ２２との間に挿脱されるように構成される。いくつかの実施形態では、演算制御ユニット２００（後述の制御部２１０）からの制御を受け、レンズユニットが被検眼Ｅと後述の対物レンズ２２との間に自動で挿脱されるように構成される。

（眼底カメラユニット２）
　眼底カメラユニット２には、被検眼Ｅの眼底Ｅｆを撮影するための光学系が設けられている。取得される眼底Ｅｆの画像（眼底像、眼底写真等と呼ばれる）は、観察画像、撮影画像等の正面画像である。観察画像は、近赤外光を用いた動画撮影により得られる。撮影画像は、フラッシュ光を用いた静止画像である。更に、眼底カメラユニット２は、被検眼Ｅの前眼部Ｅａを撮影して正面画像（前眼部像）を取得することができる。

　眼底カメラユニット２は、照明光学系１０と撮影光学系３０とを含む。照明光学系１０は被検眼Ｅに照明光を照射する。撮影光学系３０は、被検眼Ｅからの照明光の戻り光を検出する。ＯＣＴユニット１００からの測定光は、眼底カメラユニット２内の光路を通じて被検眼Ｅに導かれ、その戻り光は、同じ光路を通じてＯＣＴユニット１００に導かれる。

　照明光学系１０の観察光源１１から出力された光（観察照明光）は、曲面状の反射面を有する反射ミラー１２により反射され、集光レンズ１３を経由し、可視カットフィルタ１４を透過して近赤外光となる。更に、観察照明光は、撮影光源１５の近傍にて一旦集束し、ミラー１６により反射され、リレーレンズ１７、１８、絞り１９及びリレーレンズ２０を経由する。そして、観察照明光は、孔開きミラー２１の周辺部（孔部の周囲の領域）にて反射され、ダイクロイックミラー４６を透過し、対物レンズ２２により屈折されて被検眼Ｅ（眼底Ｅｆ又は前眼部Ｅａ）を照明する。被検眼Ｅからの観察照明光の戻り光は、対物レンズ２２により屈折され、ダイクロイックミラー４６を透過し、孔開きミラー２１の中心領域に形成された孔部を通過し、ダイクロイックミラー５５を透過する。ダイクロイックミラー５５を透過した戻り光は、撮影合焦レンズ３１を経由し、ミラー３２により反射される。更に、この戻り光は、ハーフミラー３３Ａを透過し、ダイクロイックミラー３３により反射され、集光レンズ３４によりイメージセンサ３５の受光面に結像される。イメージセンサ３５は、所定のフレームレートで戻り光を検出する。なお、撮影光学系３０のフォーカスは、眼底Ｅｆ又は前眼部Ｅａに合致するように調整される。

　撮影光源１５から出力された光（撮影照明光）は、観察照明光と同様の経路を通って眼底Ｅｆに照射される。被検眼Ｅからの撮影照明光の戻り光は、観察照明光の戻り光と同じ経路を通ってダイクロイックミラー３３まで導かれ、ダイクロイックミラー３３を透過し、ミラー３６により反射され、集光レンズ３７によりイメージセンサ３８の受光面に結像される。

　ＬＣＤ（Ｌｉｑｕｉｄ　Ｃｒｙｓｔａｌ　Ｄｉｓｐｌａｙ）３９は固視標や視力測定用視標を表示する。ＬＣＤ３９から出力された光束は、その一部がハーフミラー３３Ａにて反射され、ミラー３２に反射され、撮影合焦レンズ３１及びダイクロイックミラー５５を経由し、孔開きミラー２１の孔部を通過する。孔開きミラー２１の孔部を通過した光束は、ダイクロイックミラー４６を透過し、対物レンズ２２により屈折されて眼底Ｅｆに投射される。

　ＬＣＤ３９の画面上における固視標の表示位置を変更することにより、被検眼Ｅの固視位置を変更できる。固視位置の例として、黄斑を中心とする画像を取得するための固視位置や、視神経乳頭を中心とする画像を取得するための固視位置や、黄斑と視神経乳頭との間の眼底中心を中心とする画像を取得するための固視位置や、黄斑から大きく離れた部位（眼底周辺部）の画像を取得するための固視位置などがある。いくつかの実施形態に係る眼科装置１５００ｂは、このような固視位置の少なくとも１つを指定するためのＧＵＩ（Ｇｒａｐｈｉｃａｌ　Ｕｓｅｒ　Ｉｎｔｅｒｆａｃｅ）等を含む。いくつかの実施形態に係る眼科装置１５００ｂは、固視位置（固視標の表示位置）をマニュアルで移動するためのＧＵＩ等を含む。

　移動可能な固視標を被検眼Ｅに呈示するための構成はＬＣＤ等の表示装置には限定されない。例えば、光源アレイ（発光ダイオード（ＬＥＤ）アレイ等）における複数の光源を選択的に点灯させることにより、移動可能な固視標を生成することができる。また、移動可能な１以上の光源により、移動可能な固視標を生成することができる。

　また、眼科装置１５００ｂには、１以上の外部固視光源が設けられてもよい。１以上の外部固視光源の１つは、被検眼Ｅの僚眼に固視光を投射することが可能である。僚眼における固視光の投射位置は、変更可能である。僚眼に対する固視光の投射位置を変更することにより、被検眼Ｅの固視位置を変更することができる。外部固視光源による固視位置は、ＬＣＤ３９を用いた被検眼Ｅの固視位置と同様であってよい。例えば、複数の外部固視光源を選択的に点灯させることにより、移動可能な固視標を生成することができる。また、移動可能な１以上の外部固視光源により、移動可能な固視標を生成することができる。

　アライメント光学系５０は、被検眼Ｅに対する光学系のアライメントに用いられるアライメント指標を生成する。ＬＥＤ５１から出力されたアライメント光は、絞り５２及び５３並びにリレーレンズ５４を経由し、ダイクロイックミラー５５により反射され、孔開きミラー２１の孔部を通過する。孔開きミラー２１の孔部を通過した光は、ダイクロイックミラー４６を透過し、対物レンズ２２により被検眼Ｅに投射される。アライメント光の角膜反射光は、観察照明光の戻り光と同じ経路を通ってイメージセンサ３５に導かれる。その受光像（アライメント指標像）に基づいてマニュアルアライメントやオートアライメントを実行できる。

　フォーカス光学系６０は、被検眼Ｅに対するフォーカス調整に用いられるスプリット指標を生成する。フォーカス光学系６０は、撮影光学系３０の光路（撮影光路）に沿った撮影合焦レンズ３１の移動に連動して、照明光学系１０の光路（照明光路）に沿って移動される。反射棒６７は、照明光路に対して挿脱可能である。フォーカス調整を行う際には、反射棒６７の反射面が照明光路に傾斜配置される。ＬＥＤ６１から出力されたフォーカス光は、リレーレンズ６２を通過し、スプリット指標板６３により２つの光束に分離され、二孔絞り６４を通過し、ミラー６５により反射され、集光レンズ６６により反射棒６７の反射面に一旦結像されて反射される。更に、フォーカス光は、リレーレンズ２０を経由し、孔開きミラー２１に反射され、ダイクロイックミラー４６を透過し、対物レンズ２２により屈折されて眼底Ｅｆに投射される。フォーカス光の眼底反射光は、アライメント光の角膜反射光と同じ経路を通ってイメージセンサ３５に導かれる。その受光像（スプリット指標像）に基づいてマニュアルフォーカスやオートフォーカスを実行できる。

　ダイクロイックミラー４６は、眼底撮影用光路とＯＣＴ用光路とを合成する。ダイクロイックミラー４６は、ＯＣＴに用いられる波長帯の光を反射し、眼底撮影用の光を透過させる。ＯＣＴ用光路（測定光の光路）には、ＯＣＴユニット１００側からダイクロイックミラー４６側に向かって順に、コリメータレンズユニット４０、光路長変更部４１、光スキャナー４２、ＯＣＴ合焦レンズ４３、ミラー４４、及びリレーレンズ４５が設けられている。

　光路長変更部４１は、図４に示す矢印の方向に移動可能とされ、ＯＣＴ用光路の長さを変更する。この光路長の変更は、眼軸長に応じた光路長補正や、干渉状態の調整などに利用される。光路長変更部４１は、コーナーキューブと、これを移動する機構とを含む。

　光スキャナー４２は、被検眼Ｅの瞳孔と光学的に共役な位置に配置される。光スキャナー４２は、ＯＣＴ用光路を通過する測定光ＬＳを偏向する。光スキャナー４２は、測定光ＬＳを１次元的又は２次元的に偏向することが可能である。

　１次元的に偏向する場合、光スキャナー４２は、所定の偏向方向に所定の偏向角度範囲で測定光ＬＳを偏向するガルバノスキャナーを含む。２次元的に偏向する場合、光スキャナー４２は、第１ガルバノスキャナーと、第２ガルバノスキャナーとを含む。第１ガルバノスキャナーは、後述のＯＣＴユニット１００に含まれる干渉光学系から対物レンズ２２に至る経路における光学系であるＯＣＴ光学系８の光軸に直交する水平方向に撮影部位（眼底Ｅｆ又は前眼部）をスキャンするように測定光ＬＳを偏向する。第２ガルバノスキャナーは、ＯＣＴ光学系８の光軸に直交する垂直方向に撮影部位をスキャンするように、第１ガルバノスキャナーにより偏向された測定光ＬＳを偏向する。光スキャナー４２による測定光ＬＳのスキャンモードとしては、例えば、水平スキャン、垂直スキャン、十字スキャン、放射スキャン、円スキャン、同心円スキャン、螺旋スキャンなどがある。

　ＯＣＴ合焦レンズ４３は、ＯＣＴ用の光学系のフォーカス調整を行うために、測定光ＬＳの光路に沿って移動される。ＯＣＴ合焦レンズ４３は、被検眼Ｅの眼底Ｅｆ又はその近傍に測定光ＬＳの焦点位置を配置するための第１レンズ位置と、被検眼Ｅに照射される測定光ＬＳを平行光束にするための第２レンズ位置とを含む移動範囲で移動可能である。撮影合焦レンズ３１の移動、フォーカス光学系６０の移動、及びＯＣＴ合焦レンズ４３の移動を連係的に制御することができる。

（ＯＣＴユニット１００）
　ＯＣＴユニット１００の構成の一例を図５に示す。ＯＣＴユニット１００には、被検眼Ｅに対してＯＣＴ計測（又はＯＣＴ撮影）を行うための光学系が設けられている。この光学系は、波長掃引型（波長走査型）光源からの光を測定光と参照光とに分割し、被検眼Ｅからの測定光の戻り光と参照光路を経由した参照光とを干渉させて干渉光を生成し、この干渉光を検出する干渉光学系である。干渉光学系による干渉光の検出結果（検出信号）は、干渉光のスペクトルを示す干渉信号であり、演算制御ユニット２００に送られる。

　光源ユニット１０１は、一般的なスウェプトソースタイプの眼科装置と同様に、出射光の波長を掃引（走査）可能な波長掃引型（波長走査型）光源を含んで構成される。波長掃引型光源は、共振器を含むレーザー光源を含んで構成される。光源ユニット１０１は、人眼では視認できない近赤外の波長帯において、出力波長を時間的に変化させる。

　光源ユニット１０１から出力された光Ｌ０は、光ファイバ１０２により偏波コントローラ１０３に導かれてその偏波状態が調整される。偏波コントローラ１０３は、例えばループ状にされた光ファイバ１０２に対して外部から応力を与えることで、光ファイバ１０２内を導かれる光Ｌ０の偏波状態を調整する。

　偏波コントローラ１０３により偏波状態が調整された光Ｌ０は、光ファイバ１０４によりファイバカプラ１０５に導かれて測定光ＬＳと参照光ＬＲとに分割される。

　参照光ＬＲは、光ファイバ１１０によりコリメータ１１１に導かれて平行光束に変換され、光路長補正部材１１２及び分散補償部材１１３を経由し、光路長変更部１１４に導かれる。光路長補正部材１１２は、参照光ＬＲの光路長と測定光ＬＳの光路長とを合わせるよう作用する。分散補償部材１１３は、参照光ＬＲと測定光ＬＳとの間の分散特性を合わせるよう作用する。

　光路長変更部１１４は、図５に示す矢印の方向に移動可能とされ、参照光ＬＲの光路長を変更する。この移動により参照光ＬＲの光路の長さが変更される。この光路長の変更は、被検眼Ｅの眼軸長に応じた光路長の補正や、干渉状態の調整などに利用される。光路長変更部１１４は、例えばコーナーキューブと、これを移動する移動機構とを含んで構成される。この場合、光路長変更部１１４のコーナーキューブは、コリメータ１１１により平行光束とされた参照光ＬＲの進行方向を逆方向に折り返す。コーナーキューブに入射する参照光ＬＲの光路と、コーナーキューブから出射する参照光ＬＲの光路とは平行である。

　光路長変更部１１４を経由した参照光ＬＲは、分散補償部材１１３及び光路長補正部材１１２を経由し、コリメータ１１６によって平行光束から集束光束に変換され、光ファイバ１１７に入射する。光ファイバ１１７に入射した参照光ＬＲは、偏波コントローラ１１８に導かれてその偏波状態が調整され、光ファイバ１１９によりアッテネータ１２０に導かれて光量が調整され、光ファイバ１２１によりファイバカプラ１２２に導かれる。

　なお、図４及び図５に示す構成においては、測定光ＬＳの光路（測定光路、測定アーム）の長さを変更するための光路長変更部４１と、参照光ＬＲの光路（参照光路、参照アーム）の長さを変更するための光路長変更部１１４の双方が設けられている。しかしながら、光路長変更部４１及び１１４の一方だけが設けられていてもよい。また、これら以外の光学部材を用いて、参照光路長と測定光路長との差を変更することも可能である。

　一方、ファイバカプラ１０５により生成された測定光ＬＳは、光ファイバ１２７によりに導かれ、コリメータレンズユニット４０により平行光束とされる。平行光束にされた測定光ＬＳは、光路長変更部４１、光スキャナー４２、ＯＣＴ合焦レンズ４３、ミラー４４及びリレーレンズ４５を経由してダイクロイックミラー４６に導かれる。ダイクロイックミラー４６に導かれてきた測定光ＬＳは、ダイクロイックミラー４６により反射され、対物レンズ２２により屈折されて被検眼Ｅに照射される。測定光ＬＳは、被検眼Ｅの様々な深さ位置において散乱（反射を含む）される。このような後方散乱光を含む測定光ＬＳの戻り光は、往路と同じ経路を逆向きに進行してファイバカプラ１０５に導かれ、光ファイバ１２８を経由してファイバカプラ１２２に到達する。

　ファイバカプラ１２２は、光ファイバ１２８を介して入射された測定光ＬＳと、光ファイバ１２１を介して入射された参照光ＬＲとを合成して（干渉させて）干渉光を生成する。ファイバカプラ１２２は、所定の分岐比（例えば１：１）で、測定光ＬＳと参照光ＬＲとの干渉光を分岐することにより、一対の干渉光ＬＣを生成する。ファイバカプラ１２２から出射した一対の干渉光ＬＣは、それぞれ光ファイバ１２３、１２４により検出器１２５に導かれる。

　検出器１２５は、例えば一対の干渉光ＬＣをそれぞれ検出する一対のフォトディテクタを有し、これらによる検出結果の差分を出力するバランスドフォトダイオード（Ｂａｌａｎｃｅｄ　Ｐｈｏｔｏ　Ｄｉｏｄｅ）である。検出器１２５は、その検出結果（干渉信号）をＤＡＱ（Ｄａｔａ　Ａｃｑｕｉｓｉｔｉｏｎ　Ｓｙｓｔｅｍ）１３０に送る。ＤＡＱ１３０には、光源ユニット１０１からクロックＫＣが供給される。クロックＫＣは、光源ユニット１０１において、波長掃引型光源により所定の波長範囲内で掃引（走査）される各波長の出力タイミングに同期して生成される。光源ユニット１０１は、例えば、各出力波長の光Ｌ０を分岐することにより得られた２つの分岐光の一方を光学的に遅延させた後、これらの合成光を検出した結果に基づいてクロックＫＣを生成する。ＤＡＱ１３０は、クロックＫＣに基づき、検出器１２５の検出結果をサンプリングする。ＤＡＱ１３０は、サンプリングされた検出器１２５の検出結果を演算制御ユニット２００に送る。演算制御ユニット２００は、例えば一連の波長走査ごとに（Ａラインごとに）、検出器１２５により得られた検出結果に基づくスペクトル分布にフーリエ変換等を施すことにより、各Ａラインにおける複素振幅プロファイル又は反射強度プロファイルを形成する。更に、演算制御ユニット２００は、各Ａラインの複素振幅プロファイル又は反射強度プロファイルを画像化することにより画像データを形成する。

（処理系）
　眼科装置１５００ｂの処理系（制御系）は、演算制御ユニット２００を中心に構成される。

（演算制御ユニット２００）
　演算制御ユニット２００は、ＤＡＱ１３０から入力される検出信号から眼底ＥｆのＯＣＴ画像を形成する。そのための演算処理は、従来のスウェプトソースタイプのＯＣＴ装置と同様である。

　演算制御ユニット２００は、眼底カメラユニット２、表示装置３、及びＯＣＴユニット１００の各部を制御する。以下、表示装置３の機能は、表示部及び操作部を備えたユーザインターフェイスにより実現されるものとする。

　図６に、眼科装置１５００ｂの処理系の構成例の機能ブロック図を示す。図６において、眼科装置１５００ｂに含まれる構成要素の一部が省略されている。

　演算制御ユニット２００は、制御部２１０と、画像形成部２２０と、データ処理部２３０とを含む。演算制御ユニット２００の機能は、１以上のプロセッサにより実現される。いくつかの実施形態では、演算制御ユニット２００の機能は、制御部２１０の機能を実現する制御プロセッサと、画像形成部２２０の機能を実現する画像形成プロセッサと、データ処理部２３０の機能を実現するデータ処理プロセッサとにより実現される。

（制御部２１０）
　制御部２１０は、各種の制御を実行する。制御部２１０は、主制御部２１１と記憶部２１２とを含む。

（主制御部２１１）
　主制御部２１１は、プロセッサを含み、眼科装置１５００ｂの各部を制御する。

　眼底カメラユニット２の制御として、主制御部２１１は、観察光源１１、撮影光源１５、ＣＣＤイメージセンサ３５、３８、及びＬＥＤ５１、６１の動作制御、ＬＣＤ３９の動作制御、撮影合焦レンズ３１の移動制御、ＯＣＴ合焦レンズ４３の移動制御、反射棒６７の移動制御、アライメント光学系５０の制御、フォーカス光学系６０の移動制御、光路長変更部４１の移動制御、光スキャナー４２の動作制御などを行う。

　ＯＣＴユニット１００の制御として、主制御部２１１は、光源ユニット１０１の動作制御、光路長変更部１１４の移動制御、アッテネータ１２０の動作制御、偏波コントローラ１０３、１１８の動作制御、検出器１２５の動作制御、ＤＡＱ１３０の動作制御などを行う。

　ユーザインターフェイス２４０の制御として、主制御部２１１は、被検眼Ｅの観察画像、被検眼Ｅの撮影画像、ＯＣＴユニット１００を用いて取得された被検眼ＥのＯＣＴ画像、測定結果、ＧＵＩ等を表示部２４０Ａに表示させる。また、主制御部２１１は、ユーザによる操作部２４０Ｂに対する操作内容に対応した操作信号を受信し、受信された操作信号に基づいて眼科装置１５００ｂの各部を制御する。

　また、主制御部２１１は、被検眼Ｅに対して光学系全体を移動する移動機構１５０を制御することが可能である。

　例えば、主制御部２１１は、手動又は自動で設定された固視位置に対応するＬＣＤ３９の画面上の位置に固視標を表示する。また、主制御部２１１は、ＬＣＤ３９に表示されている固視標の表示位置を（連続的に又は段階的に）変更することができる。それにより、固視標を移動することができる（つまり、固視位置を変更することができる）。固視標の表示位置や移動態様は、マニュアルで又は自動的に設定される。マニュアルでの設定は、例えばＧＵＩを用いて行われる。自動的な設定は、例えば、データ処理部２３０により行われる。

　合焦駆動部３１Ａは、撮影光学系３０の光軸方向に撮影合焦レンズ３１を移動させるとともに、照明光学系１０の光軸方向にフォーカス光学系６０を移動させる。それにより、撮影光学系３０の合焦位置が変更される。合焦駆動部３１Ａは、撮影合焦レンズ３１を移動させる機構と、フォーカス光学系６０を移動させる機構とを個別に有していてよい。合焦駆動部３１Ａは、フォーカス調整を行うときなどに制御される。

　合焦駆動部４３Ａは、測定光路の光軸方向にＯＣＴ合焦レンズ４３を移動させる。それにより、測定光ＬＳの合焦位置が変更される。例えば、ＯＣＴ合焦レンズ４３を第１レンズ位置に移動させることにより、測定光ＬＳの合焦位置を眼底Ｅｆ又はその近傍に配置することができる。例えば、ＯＣＴ合焦レンズ４３を第２レンズ位置に移動させることにより、測定光ＬＳの合焦位置を遠点位置に配置して測定光ＬＳを平行光束にすることができる。測定光ＬＳの合焦位置は、測定光ＬＳのビームウェストの深さ位置（ｚ位置）に相当する。

　移動機構１５０は、例えば、少なくとも眼底カメラユニット２（光学系）を３次元的に移動する。典型的な例において、移動機構１５０は、少なくとも眼底カメラユニット２をｘ方向（左右方向）に移動するための機構と、ｙ方向（上下方向）に移動するための機構と、ｚ方向（奥行き方向、前後方向）に移動するための機構とを含む。ｘ方向に移動するための機構は、例えば、ｘ方向に移動可能なｘステージと、ｘステージを移動するｘ移動機構とを含む。ｙ方向に移動するための機構は、例えば、ｙ方向に移動可能なｙステージと、ｙステージを移動するｙ移動機構とを含む。ｚ方向に移動するための機構は、例えば、ｚ方向に移動可能なｚステージと、ｚステージを移動するｚ移動機構とを含む。各移動機構は、パルスモータ等のアクチュエータを含み、主制御部２１１からの制御を受けて動作する。

　移動機構１５０に対する制御は、アライメントやトラッキングにおいて用いられる。トラッキングとは、被検眼Ｅの眼球運動に合わせて装置光学系を移動させるものである。トラッキングを行う場合には、事前にアライメントとフォーカス調整が実行される。トラッキングは、装置光学系の位置を眼球運動に追従させることにより、アライメントとピントが合った好適な位置関係を維持する機能である。いくつかの実施形態では、参照光の光路長（よって、測定光の光路と参照光の光路との間の光路長差）を変更するために移動機構１５０の制御を行うように構成される。

　マニュアルアライメントの場合、光学系に対する被検眼Ｅの変位がキャンセルされるようにユーザが後述のユーザインターフェイス２４０に対して操作することにより光学系と被検眼Ｅとを相対移動させる。例えば、主制御部２１１は、ユーザインターフェイス２４０に対する操作内容に対応した制御信号を移動機構１５０に出力することにより移動機構１５０を制御して光学系と被検眼Ｅとを相対移動させる。

　オートアライメントの場合、光学系に対する被検眼Ｅの変位がキャンセルされるように主制御部２１１が移動機構１５０を制御することにより光学系と被検眼Ｅとを相対移動させる。いくつかの実施形態では、主制御部２１１は、光学系の光軸が被検眼Ｅの軸に略一致し、かつ、被検眼Ｅに対する光学系の距離が所定の作動距離になるように制御信号を移動機構１５０に出力することにより移動機構１５０を制御して光学系と被検眼Ｅとを相対移動させる。ここで、作動距離とは、対物レンズ２２のワーキングディスタンスとも呼ばれる既定値であり、光学系を用いた測定時（撮影時）における被検眼Ｅと光学系との間の距離に相当する。

　主制御部２１１は、眼底カメラユニット２等を制御することにより眼底撮影及び前眼部撮影を制御する。また、主制御部２１１は、眼底カメラユニット２及びＯＣＴユニット１００等を制御することによりＯＣＴ計測を制御する。

　演算制御ユニット２００は、例えば、従来のコンピュータと同様に、プロセッサ、ＲＡＭ、ＲＯＭ、ハードディスクドライブ、通信インターフェイスなどを含んで構成される。ハードディスクドライブ等の記憶装置には、眼科装置１５００ｂを制御するためのコンピュータプログラムが記憶されている。演算制御ユニット２００は、各種の回路基板、例えばＯＣＴ画像を形成するための回路基板を備えていてもよい。また、演算制御ユニット２００は、キーボードやマウス等の操作デバイス（入力デバイス）や、ＬＣＤ等の表示デバイスを備えていてもよい。

　眼底カメラユニット２、ユーザインターフェイス２４０、ＯＣＴユニット１００、及び演算制御ユニット２００は、一体的に（つまり単一の筺体内に）構成されていてもよいし、２つ以上の筐体に別れて構成されていてもよい。

（記憶部２１２）
　記憶部２１２は、各種のデータを記憶する。記憶部２１２に記憶されるデータとしては、例えば、ＯＣＴ画像の画像データ、眼底像の画像データ、前眼部像の画像データ、被検眼情報などがある。被検眼情報は、患者ＩＤや氏名などの被検者に関する情報や、左眼／右眼の識別情報などの被検眼に関する情報を含む。

　また、記憶部２１２には、眼科装置１５００ｂを動作させるための各種プログラムやデータが記憶されている。

　眼科装置１５００ｂには、ユーザからの操作を受け付けたり、ユーザに情報を提示したりするためのユーザインターフェイス２４０が設けられている。制御部２１０は、ユーザインターフェイス２４０を制御することによりユーザとのインターフェイス処理を司ることができる。

（ユーザインターフェイス２４０）
　ユーザインターフェイス２４０には、表示部２４０Ａと操作部２４０Ｂとが含まれる。表示部２４０Ａは、演算制御ユニット２００の表示デバイスや表示装置３を含んで構成される。操作部２４０Ｂは、前述した演算制御ユニット２００の操作デバイスを含んで構成される。操作部２４０Ｂには、眼科装置１５００ｂの筐体や外部に設けられた各種のボタンやキーが含まれていてもよい。例えば眼底カメラユニット２が従来の眼底カメラと同様の筺体を有する場合、操作部２４０Ｂは、この筺体に設けられたジョイスティックや操作パネル等を含んでいてもよい。また、表示部２４０Ａは、眼底カメラユニット２の筺体に設けられたタッチパネルなどの各種表示デバイスを含んでいてもよい。

　なお、表示部２４０Ａと操作部２４０Ｂは、それぞれ個別のデバイスとして構成される必要はない。例えばタッチパネルのように、表示機能と操作機能とが一体化されたデバイスを用いることも可能である。その場合、操作部２４０Ｂは、このタッチパネルとコンピュータプログラムとを含んで構成される。操作部２４０Ｂに対する操作内容は、電気信号として制御部２１０に入力される。また、表示部２４０Ａに表示されたグラフィカルユーザインターフェイス（ＧＵＩ）と、操作部２４０Ｂとを用いて、操作や情報入力を行うようにしてもよい。

（画像形成部２２０）
　画像形成部２２０は、検出器１２５からの検出信号をＤＡＱ１３０でサンプリングすることにより得られたサンプリングデータに基づいて、被検眼ＥのＯＣＴ画像（画像データ）を形成する。この処理には、従来のスウェプトソースタイプのＯＣＴと同様に、ノイズ除去（ノイズ低減）、フィルタ処理、分散補償、ＦＦＴ（Ｆａｓｔ　Ｆｏｕｒｉｅｒ　Ｔｒａｎｓｆｏｒｍ）などの処理が含まれている。他のタイプのＯＣＴ装置の場合、画像形成部２２０は、そのタイプに応じた公知の処理を実行する。

　画像形成部２２０は、例えば、前述のプロセッサ、ＲＡＭ、ＲＯＭ、ハードディスクドライブ、回路基板等を含んで構成される。ハードディスクドライブ等の記憶装置には、上記機能をプロセッサに実行させるコンピュータプログラムがあらかじめ格納されている。なお、この明細書では、「画像データ」と、それに基づく「画像」とを同一視することがある。

（データ処理部２３０）
　データ処理部２３０は、被検眼Ｅの撮影やＯＣＴ計測により取得されたデータを処理する。

　データ処理部２３０は、画像形成部２２０により形成された画像に対して各種の画像処理や解析処理を施す。例えば、データ処理部２３０は、画像の輝度補正等の各種補正処理を実行する。また、データ処理部２３０は、眼底カメラユニット２により得られた画像（眼底像、前眼部像等）に対して各種の画像処理や解析処理を施す。

　例えば、データ処理部２３０は、断層像の間の画素を補間する補間処理などの公知の画像処理を実行して、眼底Ｅｆの３次元画像の画像データを形成する。なお、３次元画像の画像データとは、３次元座標系により画素の位置が定義された画像データを意味する。３次元画像の画像データとしては、３次元的に配列されたボクセルからなる画像データがある。この画像データは、ボリュームデータ或いはボクセルデータなどと呼ばれる。ボリュームデータに基づく画像を表示させる場合、データ処理部２３０は、このボリュームデータに対してレンダリング処理（ボリュームレンダリングやＭＩＰ（Ｍａｘｉｍｕｍ　Ｉｎｔｅｎｓｉｔｙ　Ｐｒｏｊｅｃｔｉｏｎ：最大値投影）など）を施して、特定の視線方向から見たときの擬似的な３次元画像の画像データを形成する。ユーザインターフェイス２４０（表示部２４０Ａ）には、この擬似的な３次元画像が表示される。

　また、３次元画像の画像データとして、複数の断層像のスタックデータを形成することも可能である。スタックデータは、複数のスキャンラインに沿って得られた複数の断層像を、スキャンラインの位置関係に基づいて３次元的に配列させることで得られる画像データである。すなわち、スタックデータは、元々個別の２次元座標系により定義されていた複数の断層像を、１つの３次元座標系により表現する（つまり１つの３次元空間に埋め込む）ことにより得られる画像データである。

　データ処理部２３０は、取得された３次元データセット（ボリュームデータ、スタックデータ等）に各種のレンダリングを施すことで、任意断面におけるＢモード画像（縦断面像、軸方向断面像）、任意断面におけるＣモード画像（横断面像、水平断面像）、プロジェクション画像、シャドウグラムなどを形成することができる。Ｂモード画像やＣモード画像のような任意断面の画像は、指定された断面上の画素（ピクセル、ボクセル）を３次元データセットから選択することにより形成される。プロジェクション画像は、３次元データセットを所定方向（ｚ方向、深さ方向、軸方向）に投影することによって形成される。シャドウグラムは、３次元データセットの一部（たとえば特定層に相当する部分データ）を所定方向に投影することによって形成される。Ｃモード画像、プロジェクション画像、シャドウグラムのような、被検眼の正面側を視点とする画像を正面画像（ｅｎ－ｆａｃｅ画像）と呼ぶ。

　データ処理部２３０は、ＯＣＴにより時系列に収集されたデータ（例えば、Ｂスキャン画像データ）に基づいて、網膜血管や脈絡膜血管が強調されたＢモード画像や正面画像（血管強調画像、アンギオグラム）を構築することができる。例えば、被検眼Ｅの略同一部位を反復的にスキャンすることにより、時系列のＯＣＴデータを収集することができる。

　いくつかの実施形態では、データ処理部２３０は、略同一部位に対するＢスキャンにより得られた時系列のＢスキャン画像を比較し、信号強度の変化部分の画素値を変化分に対応した画素値に変換することにより当該変化部分が強調された強調画像を構築する。更に、データ処理部２３０は、構築された複数の強調画像から所望の部位における所定の厚さ分の情報を抽出してｅｎ－ｆａｃｅ画像として構築することでＯＣＴＡ像を形成する。

　データ処理部２３０により生成された画像（例えば、３次元画像、Ｂモード画像、Ｃモード画像、プロジェクション画像、シャドウグラム、ＯＣＴＡ像）もまたＯＣＴ画像に含まれる。

　更に、データ処理部２３０は、ＯＣＴ計測により得られた干渉光の検出結果、又は当該検出結果に基づいて形成されたＯＣＴ画像に対して所定の解析処理を行う。所定の解析処理には、被検眼Ｅにおける所定の部位（組織、病変部）の特定；指定された部位間の距離（層間距離）、面積、角度、比率、密度の算出；指定された計算式による演算；所定の部位の形状の特定；これらの統計値の算出；計測値、統計値の分布の算出；これら解析処理結果に基づく画像処理などがある。所定の組織には、血管、視神経乳頭、中心窩、黄斑などがある。所定の病変部には、白斑、出血などがある。

　また、データ処理部２３０は、被検眼の動きやＯＣＴ光学系８の位相不安定性に起因した位相ドリフトの影響を低減するように、ＯＣＴ計測により得られた被検眼ＥのＯＣＴデータ（例えば、３次元複素ＯＣＴデータ、複素ＯＣＴボリュームデータ）に対して位相安定化処理を実行することが可能である。

　例えば、画像形成部２２０又はデータ処理部２３０は、位相安定化処理が施された複素ＯＣＴデータの絶対値に基づいてＯＣＴ画像を形成する。複素ＯＣＴデータの実数部に加えて虚数部を用いて画像化することができるので、より一層高精細なＯＣＴ画像を取得することができる。

　例えば、データ処理部２３０は、位相安定化処理が施された複素ＯＣＴデータの絶対値に基づいて位相画像（位相情報）やＯＣＴＡ画像を形成する。複素ＯＣＴデータの虚数部から位相ドリフトの影響を受けた成分が除去されているため、より一層高精度な位相画像やＯＣＴＡ画像を取得することができる。

　更に、データ処理部２３０は、位相安定化処理が施されたＯＣＴデータに対して、被検眼Ｅや光学系に起因した収差を補正する収差補正処理を実行することが可能である。

　図７に、データ処理部２３０により実行される収差補正処理の動作説明図を示す。

　まず、データ処理部２３０は、ＯＣＴ計測により得られた被検眼Ｅの複素ＯＣＴデータＤＴ１からフィルタ情報としてのＣＡＯフィルタＣＦのパラメータを探索する。以下、フィルタ情報としての「ＣＡＯフィルタのパラメータの探索」を「ＣＡＯフィルタの探索」と表記する場合がある。

　次に、データ処理部２３０は、探索されたＣＡＯフィルタＣＦを用いて複素ＯＣＴデータＤＴ１に対して収差補正処理ＰＣを実行し、収差補正が行われた複素ＯＣＴデータＤＴ２を取得する。いくつかの実施形態では、収差補正処理ＰＣは、周波数領域において実行される。いくつかの実施形態では、収差補正処理ＰＣは、空間領域において実行される。

　それにより、被検眼Ｅや光学系に起因した収差に影響を受けることなくＯＣＴ計測結果の精度を向上させることができる。これは、収差補正が行われた複素ＯＣＴデータの実数部と虚数部の精度の向上を意味し、ＯＣＴ画像（ＯＣＴＡ画像を含む）の高画質化や画像間の位相の精度向上が可能になる。

　図８に、実施形態に係るＣＡＯフィルタの説明図を示す。

　ＣＡＯフィルタは、被検眼Ｅの瞳面における波面収差の分布に対応したゼルニケ多項式（Ｚｅｒｎｉｋｅ　ｐｏｌｙｎｏｍｉａｌ）により表すことができる。ゼルニケ多項式の各次数の項の係数（パラメータ）を調整することにより、ＣＡＯフィルタのフィルタ特性を調整することが可能である。

　ＣＡＯフィルタのフィルタ特性は、あらかじめ決められた波面収差の分布に対応した１以上のフィルタを、各フィルタに対応した係数に応じて合成することにより決定される。実施形態では、図８に示すように、ＣＡＯフィルタＣＦは、デフォーカス項Ｆ_ｄｅと、ＨＶ方向（０－９０度方向）の非点収差項Ｆ_ａｓＨＶと、対角方向（±４５度方向）の非点収差項Ｆ_ａｓＤＧとを含む多項式で表現される。ＣＡＯフィルタＣＦは、各項の係数Ｃ_ｄｅ、Ｃ_{ａｓ＿ＨＶ}、Ｃ_{ａｓ＿ＤＧ}をパラメータとして調整することで所望のフィルタ特性を実現する。なお、ＣＡＯフィルタＣＦは、コマ収差項、球面収差項などのより高次の項を含んでよい。

　図９に、実施形態に係るＣＡＯフィルタの探索処理の一例の説明図を示す。

　データ処理部２３０は、互いにフィルタ特性が異なる２以上のＣＡＯフィルタを被検眼Ｅの複素ＯＣＴデータＤＴ１に順次に適用して、収差補正が行われた２以上の複素ＯＣＴデータＤＴ１１、ＤＴ１２、・・・、ＤＴ１Ｍ（Ｍは２以上の整数）を取得する。データ処理部２３０は、取得された２以上の複素ＯＣＴデータＤＴ１１～ＤＴ１Ｍのそれぞれの品質を評価し、最も品質が高い複素ＯＣＴデータに適用されたＣＡＯフィルタを探索結果のＣＡＯフィルタとして特定する。

　複素ＯＣＴデータの品質を評価する手法の例として、当該複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価する手法がある。ＯＣＴ画像の例として、断層像、正面画像（プロジェクション画像、ｅｎ－ｆａｃｅ画像）、ＯＣＴＡ画像などがある。例えば、データ処理部２３０は、形成されたＯＣＴ画像の画質の程度を表す画質評価値に基づいて、最も品質が高い複素ＯＣＴデータを特定することが可能である。

　いくつかの実施形態では、データ処理部２３０は、被検眼Ｅの屈折力を表す屈折力情報から上記のＣＡＯフィルタのパラメータ（係数）を生成することが可能である。これにより、ＣＡＯフィルタ（フィルタ情報）の探索処理を簡略化、又は省略化を図ることができるようになる。

　以下、このような機能を実現するデータ処理部２３０の各部について説明する。

　図１０に、図６のデータ処理部２３０の構成例のブロック図を示す。

　データ処理部２３０は、位相安定化部２３１と、フィルタ情報生成部２３２と、フィルタ情報探索部３００と、収差補正部２３４と、プロジェクション画像形成部２３５とを含む。

（位相安定化部２３１）
　位相安定化部２３１は、計測対象である被検眼Ｅの動きやシステム（ＯＣＴ光学系８）内の位相不安定性に起因した位相ドリフトの影響をキャンセルするように被検眼Ｅの複素ＯＣＴデータに対して位相安定化処理を行う。

　図１１に、図１０の位相安定化部２３１の構成例のブロック図を示す。

　位相安定化部２３１は、位相差プロファイル生成部２３１Ａと、位相ドリフト抽出部２３１Ｂと、位相差データ除去部２３１Ｃと、位相補正部２３１Ｄとを含む。

（位相差プロファイル生成部２３１Ａ）
　位相差プロファイル生成部２３１Ａは、３次元複素ＯＣＴデータからＡラインごとに平均化された位相差を表す位相データがＢライン方向に配列された位相差プロファイルを生成する。具体的には、位相差プロファイル生成部２３１Ａは、複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに位相差を算出し、算出された位相差を深さ方向（ｚ方向）に平均化し、平均化して得られた位相差をＢライン方向に配列して位相差プロファイルを生成する。

　図１２Ａに、位相差プロファイル生成部２３１Ａの動作説明図を示す。図１２Ａは、横方向をｘ方向とし、縦方向をｙ方向として、複素ＯＣＴデータＯＤＴを模式的に表したものである。

　例えば、位相差プロファイル生成部２３１Ａは、式（１）に示すように、隣接する２つのＢフレームＢ_ｋ、Ｂ_ｋ＋１間のＡラインの深さ位置ごとに算出された位相差をｚ方向に積算し、積算された位相差を平均化して、Ａラインごとに位相差を求める。位相差プロファイル生成部２３１Ａは、次の２つのＢフレームＢ_ｋ＋１、Ｂ_ｋ＋２間で同様に、Ａラインごとに位相差を求める。これを順次に繰り返し、３次元複素ＯＣＴデータＯＤＴについて２つのＢフレーム間の位相差プロファイルが順次に求められる。

（位相ドリフト抽出部２３１Ｂ）
　位相ドリフト抽出部２３１Ｂは、位相差プロファイル生成部２３１Ａにより生成された位相差プロファイルから位相ドリフトを抽出する。具体的には、位相ドリフト抽出部２３１Ｂは、位相差プロファイルに対して、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つを施すことにより位相ドリフトを抽出する。

　図１２Ｂに、位相ドリフト抽出部２３１Ｂの動作説明図を示す。図１２Ｂにおいて、横軸はＡラインの位置（Ｂスキャン方向に走査される測定光の入射位置）に相当するＡライン番号を表し、縦軸は位相差を表す。

　例えば、位相ドリフト抽出部２３１Ｂは、位相差プロファイル生成部２３１Ａにより生成された位相差プロファイルＲＤ０に対してスムージング処理を行うことで、位相ドリフトＤＦ０を抽出する。

　位相ドリフト抽出部２３１Ｂにより実行されるスムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つのパラメータは、複素ＯＣＴデータＯＤＴを取得するＯＣＴ光学系８の位相変化の周波数、位相変化の速度、及びＯＣＴ光学系８におけるＯＣＴスキャンに要する時間の少なくとも１つに基づいて設定される。

（位相差データ除去部２３１Ｃ）
　位相差データ除去部２３１Ｃは、上記の位相ドリフトの抽出の前に、位相差プロファイル生成部２３１Ａにより生成された位相差プロファイルから信頼性の低い位相差データを除去する。位相ドリフト抽出部２３１Ｂは、位相差データ除去部２３１Ｃにより信頼性の低い位相差データが除去された位相差プロファイルから位相ドリフトを抽出する。それにより、局所的な位相差変化が除去された位相ドリフトが得られる。

　図１３に、位相差データ除去部２３１Ｃの動作説明図を示す。図１３において、図１２Ｂと同様に、横軸はＡラインの位置に相当するＡライン番号を表し、縦軸は位相差を表す。

　例えば、位相差データ除去部２３１Ｃは、図１２Ｂに示す位相差プロファイルＲＤ０から信頼性の低い位相差データＮ１～Ｎ４を除去して（間引いて）、図１３に示す位相差プロファイルＲＤ１を生成する。位相ドリフト抽出部２３１Ｂは、図１３に示す位相差プロファイルＲＤ１から位相ドリフトＤＦ１を抽出する。これにより、位相ドリフトＤＦ０と比較して、局所的な位相変化（ローカルエラー）による影響を抑えることができる。

　信頼性の低い位相差データとして、信号強度が弱い位相差データ、式（１）の位相差の算出処理（平均化処理）において特定の深さ位置における位相差データの寄与度が大きい位相差データなどがある。信号強度が弱い位相差データは、信号強度が強い位相差データが有する高信頼性を低下させる可能性がある。式（１）の位相差の算出処理において寄与度が大きい位相差データは、位相差の期待値が零になるという前提条件を崩す可能性が高くなる。

　例えば、位相差データ除去部２３１Ｃは、式（２）を満たす位相差データを信号強度が所定の閾値レベル未満であるデータとして特定し、特定されたデータを位相差プロファイルから除去する。

　例えば、位相差データ除去部２３１Ｃは、式（３）を満たす位相差データを式（１）の位相差を算出するときの平均化処理において所定の深さ位置における位相差による寄与が所定の閾値を超えるデータとして特定し、特定されたデータを位相差プロファイルから除去する。

　位相差データ除去部２３１Ｃは、信号強度が所定の閾値レベル未満であるデータ、及び平均化処理において所定の深さ位置における位相差による寄与が所定の閾値を超えるデータの少なくとも１つを信頼性の低いデータとして位相差プロファイルから除去することが可能である。

（位相補正部２３１Ｄ）
　位相補正部２３１Ｄは、位相ドリフト抽出部２３１Ｂにより抽出された位相ドリフトに基づいて、位相差プロファイルが生成された上記の２つのＢフレームのいずれか一方の複素ＯＣＴデータの位相を補正する。

　例えば、図１２Ａに示すＢフレームＢ_ｋ、Ｂ_ｋ＋１間で位相補正を行う場合、位相補正部２３１Ｄは、Ａラインごとに位相ドリフトにおける位相差φを用いて、周波数領域におけるＢフレームＢ_ｋ＋１の複素ＯＣＴデータに対してｅｘｐ（－ｉφ）を乗算する。いくつかの実施形態では、位相補正部２３１Ｄは、空間領域における複素ＯＣＴデータに対して、位相差φを用いた畳み込み演算を行う。

　例えば、位相安定化部２３１は、３次元複素ＯＣＴデータのＢフレームに対して上記の位相補正処理を順次に実行する。すなわち、位相安定化部２３１は、隣接する第１のＢフレーム及び第２のＢフレームについて、位相差プロファイル生成部２３１Ａにより生成された位相差プロファイルに基づいて第２のＢフレームの複素ＯＣＴデータの位相を補正した後、第２のＢフレーム及び第２のＢフレームに隣接する第３のＢフレームについて、位相差プロファイル生成部２３１Ａにより生成された位相差プロファイルに基づいて第３のＢフレームの複素ＯＣＴデータの位相を補正するように構成される。

（フィルタ情報生成部２３２）
　図１０に示すフィルタ情報生成部２３２は、フィルタ情報として、被検眼Ｅの複素ＯＣＴデータに適用するＣＡＯフィルタのパラメータを生成する。図８に示す例では、ＣＡＯフィルタのパラメータは、ゼルニケ多項式の各項の係数である。

　フィルタ情報生成部２３２は、探索処理の初期値となるＣＡＯフィルタのパラメータをフィルタ情報として生成する。例えば、フィルタ情報生成部２３２は、記憶部２１２にあらかじめ記憶されたＣＡＯフィルタのパラメータの初期値を参照してフィルタ情報を取得する。

　いくつかの実施形態では、フィルタ情報生成部２３２は、外部の眼科装置において実行された被検眼Ｅの屈折力測定により得られた被検眼Ｅの屈折力情報からＣＡＯフィルタのパラメータをフィルタ情報として生成する。すなわち、フィルタ情報生成部２３２は、被検眼Ｅの屈折力情報に基づいてフィルタ情報を基準フィルタ情報として生成することができる。ＣＡＯフィルタのパラメータは、生成された基準フィルタ情報に基づいて複素ＯＣＴデータの品質が所定のレベルになるように探索される。

　等価球面度数ＳＥ、球面度数Ｓ、乱視度数Ｃ、及び乱視軸角度Ａは、ゼルニケ多項式の各項の係数を用いて式（４）～式（７）により算出可能である。

　式（４）において、Ｓ_ｍｏｖｅは、固視移動分の球面度数を表し、ｒは瞳径を表し、ｃ_２ ^０は、デフォーカス項の係数（図８の係数Ｃ_ｄｅに相当）を表す。また、式（６）、式（７）において、ｃ_２ ^－２は、ＤＧ方向の非点収差項の係数（図８の係数Ｃ_{ａｓ＿ＤＧ}に相当）を表し、ｃ_２ ^２は、ＨＶ方向の非点収差項の係数（図８の係数Ｃ_{ａｓ＿ＨＶ}に相当）を表す。

　係数がＲＭＳ（Ｒｏｏｔ　Ｍｅａｎ　Ｓｑｕａｒｅ）値の場合、等価球面度数ＳＥ、球面度数Ｓ、乱視度数Ｃ、及び乱視軸角度Ａは、式（８）～式（１１）により算出可能である。

　式（４）～式（７）又は式（８）～式（１１）に従って、フィルタ情報生成部２３２は、被検眼Ｅの屈折力情報からゼルニケ多項式の各項の係数を算出することでＣＡＯフィルタのパラメータをフィルタ情報として生成することができる。

（フィルタ情報探索部３００）
　フィルタ情報探索部３００は、被検眼Ｅの複素ＯＣＴデータの品質が所定レベルになるように、複素ＯＣＴデータの収差を補正するためのフィルタ情報を探索する。具体的には、フィルタ情報探索部３００は、複素ＯＣＴデータの品質が所定レベルになるように、フィルタ情報に基づいて補正された複素ＯＣＴデータの品質の評価と、フィルタ情報の更新とを繰り返す。更に、フィルタ情報探索部３００は、品質の評価とフィルタ情報の更新とを繰り返すことによりフィルタ情報を粗く探索し、探索された探索範囲内でフィルタ情報の更新と品質の評価とを繰り返すことによりフィルタ情報を細密に探索することが可能である。

　図１４に、図１０のフィルタ情報探索部３００の構成例のブロック図を示す。

　フィルタ情報探索部３００は、粗探索部３１０と、細密探索部３２０とを含む。

（粗探索部３１０）
　粗探索部３１０は、複素ＯＣＴデータの品質の評価とフィルタ情報の更新とを繰り返すことによりフィルタ情報を粗く探索する。このような粗探索部３１０は、フィルタテスト部３１１と、フィルタ情報更新部３１２とを含む。

（フィルタテスト部３１１）
　フィルタテスト部３１１は、複素ＯＣＴデータに対してＣＡＯフィルタを適用し、収差が補正された複素ＯＣＴデータの品質を評価する。このようなフィルタテスト部３１１は、フィルタ処理部３１１Ａと、画質評価部３１１Ｂとを含む。

（フィルタ処理部３１１Ａ）
　フィルタ処理部３１１Ａは、粗探索処理中のＣＡＯフィルタを用いて、複素ＯＣＴデータに対して公知の手法でフィルタ処理を実行する。例えば、フィルタ処理部３１１Ａは、周波数領域における複素ＯＣＴデータに対して、ＣＡＯフィルタのパラメータに対応したｅｘｐ（－ｉφ）を乗算する。例えば、フィルタ処理部３１１Ａは、空間領域における複素ＯＣＴデータに対して、ＣＡＯフィルタのパラメータに対応したフィルタ係数を用いた畳み込み演算を行う。

（画質評価部３１１Ｂ）
　画質評価部３１１Ｂは、フィルタ処理部３１１Ａにより収差が補正された複素ＯＣＴデータに基づいて形成されるＯＣＴ画像（眼底像）の画質を評価する。例えば、ＯＣＴ画像は、画像形成部２２０により形成される。

　画質評価部３１１Ｂは、ＯＣＴ画像の画質評価値を算出することが可能である。画質評価値は、ＯＣＴ画像の画質の程度を示す数値であり、画質が良好か否かの評価に用いられる。

　画質評価値は、任意の手法で算出可能な任意の種類の値であってよい。典型的な画質評価値は、ＯＣＴ画像中の２種類の画像領域に関する信号対ノイズ比（Ｓ／Ｎ比）で表される。画質評価値の具体例は、眼の部位を表す画像領域の信号と、眼の部位ではない背景を表す画像領域のノイズとに関するＳ／Ｎ比である。眼の部位を表す画像領域を信号領域と呼び、背景を表す画像領域を背景領域と呼ぶ。

　信号領域及び背景領域に関するＳ／Ｎ比として表現される画質評価値の算出方法は任意である。その具体例を以下に説明する。

　まず、画質評価部３１１Ｂは、ＯＣＴ画像を構成する複数のＡスキャン画像のそれぞれについて、最大輝度の画素と最小輝度の画素とを特定する。次に、画質評価部３１１Ｂは、特定された最大輝度の画素を含む所定範囲の画素群（例えば、前後４０画素分）の輝度値に基づいて、輝度値のヒストグラムを作成する。同様に、画質評価部３１１Ｂは、特定された最小輝度の画素を含む所定範囲の画素群の輝度値に基づいて、輝度値のヒストグラムを作成する。

　続いて、画質評価部３１１Ｂは、最小輝度の画素を含む画素群に対応するヒストグラムにおいて、頻度値が０を超える最大の位置（輝度値）を探索する。更に、最大輝度の画素を含む画素群に対応するヒストグラムにおいて、上記で探索された輝度値以下の範囲に含まれる合計画素数（Ｎ）と、探索された輝度値より上から２５５番目の輝度値に含まれている合計画素数（Ｓ）とを求める。そして、画質評価部３１１Ｂは、当該画像において信号とみなせる部分（つまりノイズではないとみなせる部分）が全体の何％になるのかを、次の演算式で評価する：１００×Ｓ÷（Ｓ＋Ｎ）。画質評価部３１１Ｂは、このような一連の演算を各Ａスキャン画像に対して適用することで、複数のＡスキャン画像に対応する複数の数値を得る。そして、画質評価部３１１Ｂは、これら数値の平均値を求め、これを画質評価値とする。

（フィルタ情報更新部３１２）
　フィルタ情報更新部３１２は、粗探索部３１０により探索されるＣＡＯフィルタのパラメータを所定の増加分だけ増加させたり所定の減少分だけ減少させたりして、当該パラメータを更新する。増加分又は減少分は、パラメータの種別に応じて異なってもよいし、２以上の種別のパラメータで共通であってよい。いくつかの実施形態では、フィルタ情報更新部３１２は、ＣＡＯフィルタのパラメータの１つだけを更新する。いくつかの実施形態では、フィルタ情報更新部３１２は、ＣＡＯフィルタのパラメータのうち２以上のパラメータを同時に更新する。

　以上のように、粗探索部３１０では、フィルタテスト部３１１は、例えば、周波数領域において、フィルタ処理部３１１ＡによりＣＡＯフィルタを適用して収差が補正された複素ＯＣＴデータを取得する。その後、フィルタテスト部３１１は、空間領域において、画質評価部３１１Ｂにより複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価する。粗探索部３１０は、フィルタ情報更新部３１２によるＣＡＯフィルタのパラメータの更新と、フィルタテスト部３１１による画質の評価とを繰り返すことで、ＣＡＯフィルタを粗く探索する。

（細密探索部３２０）
　細密探索部３２０は、粗探索部３１０により探索された探索範囲内でフィルタ情報の更新と品質の評価とを繰り返すことによりフィルタ情報を細密に探索する。このような細密探索部３２０は、フィルタテスト部３２１と、フィルタ情報更新部３２２と、勾配算出部３２３とを含む。

（フィルタテスト部３２１）
　フィルタテスト部３２１は、フィルタテスト部３１１と同様に、複素ＯＣＴデータに対してＣＡＯフィルタを適用し、収差が補正された複素ＯＣＴデータの品質を評価する。このようなフィルタテスト部３２１は、フィルタ処理部３２１Ａと、画質評価部３２１Ｂとを含む。

（フィルタ処理部３２１Ａ）
　フィルタ処理部３２１Ａは、細密探索処理中のＣＡＯフィルタを用いて、複素ＯＣＴデータに対して公知の手法でフィルタ処理を実行する。フィルタ処理部３２１Ａによる処理は、フィルタ処理部３１１Ａと同様である。

（画質評価部３２１Ｂ）
　画質評価部３２１Ｂは、フィルタ処理部３２１Ａにより収差が補正された複素ＯＣＴデータに基づいて形成されるＯＣＴ画像（眼底像）の画質を評価する。画質評価部３２１Ｂは、画質評価部３１１Ｂと同様に、ＯＣＴ画像の画質評価値を算出することが可能である。

（フィルタ情報更新部３２２）
　フィルタ情報更新部３２２は、細密探索部３２０により探索されるＣＡＯフィルタのパラメータを増加又は減少させて、当該パラメータを更新する。フィルタ情報更新部３２２は、後述する勾配算出部３２３により算出された勾配に応じてＣＡＯフィルタのパラメータを更新することが可能である。いくつかの実施形態では、フィルタ情報更新部３２２は、ＣＡＯフィルタのパラメータの１つだけを更新する。いくつかの実施形態では、フィルタ情報更新部３２２は、ＣＡＯフィルタのパラメータのうち２以上のパラメータを同時に更新する。

（勾配算出部３２３）
　勾配算出部３２３は、画質評価部３２１Ｂにより算出された画質評価値の勾配を算出する。

　例えば、勾配算出部３２３は、フィルタ情報の更新によって変化する画質評価値の勾配を算出する。フィルタ情報更新部３２２は、勾配算出部３２３により算出された画質評価値の勾配に基づいてフィルタ情報を更新する。具体的には、フィルタ情報更新部３２２は、算出された画質評価値の勾配が緩やかなほど増加分又は減少分が小さくなるようにフィルタ情報を更新し、勾配が急峻なほど増加分又は減少分が大きくなるようにフィルタ情報を更新する。

　以上のように、細密探索部３２０では、フィルタテスト部３２１は、例えば、周波数領域において、フィルタ処理部３２１ＡによりＣＡＯフィルタを適用して収差が補正された複素ＯＣＴデータを取得する。その後、フィルタテスト部３２１は、空間領域において、画質評価部３２１Ｂにより複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価する。細密探索部３２０では、勾配算出部３２３は、フィルタ情報の更新によって変化する画質評価値の勾配を算出し、フィルタ情報更新部３２２は、算出された画質評価値の勾配に基づいてフィルタ情報を更新する。細密探索部３２０は、画質評価値の勾配に応じたフィルタ情報更新部３２２によるＣＡＯフィルタのパラメータの更新と、フィルタテスト部３２１による画質の評価とを繰り返すことで、ＣＡＯフィルタを細密に探索する。

（収差補正部２３４）
　図１０に示す収差補正部２３４は、フィルタ情報探索部３００により探索されたＣＡＯフィルタ（パラメータ）に基づいて、被検眼Ｅの複素ＯＣＴデータの収差を補正する。収差補正部２３４は、ＣＡＯフィルタを用いて複素ＯＣＴデータに対して公知の手法でフィルタ処理を実行する。例えば、収差補正部２３４は、周波数領域における複素ＯＣＴデータに対して、ＣＡＯフィルタのパラメータに対応したｅｘｐ（－ｉφ´）を乗算する。例えば、収差補正部２３４は、空間領域における複素ＯＣＴデータに対して、ＣＡＯフィルタのパラメータに対応したフィルタ係数を用いた畳み込み演算を行う。

　フィルタ処理部３１１Ａ及びフィルタ処理部３２１Ａの少なくとも一方の機能は、収差補正部２３４により実現されてよい。

（プロジェクション画像形成部２３５）
　プロジェクション画像形成部２３５は、被検眼Ｅの複素ＯＣＴデータに基づいてプロジェクション画像を形成する。例えば、プロジェクション画像形成部２３５は、収差補正部２３４により収差が補正された複素ＯＣＴデータに基づいてプロジェクション画像を形成する。

　プロジェクション画像形成部２３５は、被検眼Ｅの３次元画像をｚ方向に投影することによりプロジェクション画像を形成する。例えば、プロジェクション画像は、複素ＯＣＴデータの品質の良否を判断するために用いられる。

　以上のような構成を有する眼科装置１５００ｂにおいて、実施形態に係る主制御部２１１は、表示制御部として、探索されたＣＡＯフィルタのパラメータに対応した被検眼Ｅの瞳面における収差情報を表示部２４０Ａに表示させることが可能である。収差情報は、被検眼Ｅの屈折力を表す情報、ゼルニケ係数、及び光の位相ずれの分布情報（波面収差情報）の少なくとも１つを含む。屈折力を表す情報は、等価球面度数（ＳＥ）、球面度数（Ｓ）、乱視度数（Ｃ）、及び乱視軸角度（Ａ）の少なくとも１つを含む。

　以上のように機能するデータ処理部２３０は、例えば、前述のプロセッサ、ＲＡＭ、ＲＯＭ、ハードディスクドライブ、回路基板等を含んで構成される。ハードディスクドライブ等の記憶装置には、上記機能をプロセッサに実行させるコンピュータプログラムがあらかじめ格納されている。

　眼科情報処理装置１２００、眼科情報処理部１２００ａ、１２００ｂ、データ処理部２３０は、実施形態に係る「眼科情報処理装置」の一例である。位相差プロファイル生成部２３１Ａは、実施形態に係る「生成部」の一例である。位相ドリフト抽出部２３１Ｂは、実施形態に係る「抽出部」の一例である。位相差データ除去部２３１Ｃは、実施形態に係る「除去部」の一例である。画像形成部２２０又はデータ処理部２３０は、実施形態に係る「画像形成部」の一例である。ＯＣＴ光学系１１０１又はＯＣＴ光学系８は、実施形態に係る「測定系」の一例である。

［動作］
　実施形態に係る眼科装置１５００ｂの動作について説明する。

　図１５～図２３に、実施形態に係る眼科装置１５００ｂの動作例を示す。図１５～図１７、図１９～図２３は、実施形態に係る眼科装置１５００ｂの動作例のフローチャートを表す。図１６は、図１５のステップＳ４の動作例のフローチャートを表す。図１７は、図１６のステップＳ１１の動作例のフローチャートを表す。図１９は、図１６のステップＳ１３の動作例のフローチャートを表す。図２０は、図１９のステップＳ２２～ステップＳ２４のそれぞれの動作例のフローチャートを表す。図２１は、図２０のステップＳ３１の動作例のフローチャートを表す。図２２は、図１９のステップＳ２５の動作例のフローチャートを表す。図２３は、図２２のステップＳ５４の動作例のフローチャートを表す。記憶部２１２には、図１５～図１７、図１９～図２３に示す処理を実現するためのコンピュータプログラムが記憶されている。主制御部２１１は、このコンピュータプログラムに従って動作することにより、図１５～図１７、図１９～図２３に示す処理を実行する。図１８は、図１６のステップＳ１２～ステップＳ１５の動作説明図を表す。

（Ｓ１：スキャン条件を設定）
　まず、主制御部２１１は、ユーザによるスキャン条件の指定を受け付ける。

　ユーザは、操作部２４０Ｂに対する操作によりスキャンモード及びスキャン範囲を指定することが可能である。操作部２４０Ｂに対する操作によりユーザによってスキャンモード及びスキャン範囲が指定されたとき、主制御部２１１は、操作部２４０Ｂからの操作情報を解析して、指定されたスキャンモード及びスキャン範囲を特定する。

　いくつかの実施形態では、主制御部２１１は、ユーザにより指定された動作モードに基づいて、動作モードに対応したスキャンモード及びスキャン範囲の少なくとも１つを設定する。

（Ｓ２：アライメント）
　次に、主制御部２１１は、アライメントを実行する。

　例えば、主制御部２１１は、アライメント光学系５０を制御して、被検眼Ｅにアライメント指標を投影させる。このとき、被検眼Ｅには、ＬＣＤ３９による固視標も投影される。主制御部２１１は、例えばイメージセンサ３５により取得された受光像に基づいて特定された光学系の移動量に基づいて移動機構１５０を制御し、被検眼Ｅに対して光学系を当該移動量だけ相対的に移動させる。主制御部２１１は、この処理を繰り返し実行させる。これにより、光学系の光軸が被検眼Ｅの軸に略一致し、かつ、被検眼Ｅに対する光学系の距離が所定の作動距離になるように光学系が配置される。

（Ｓ３：ＯＣＴスキャン）
　続いて、主制御部２１１は、ＯＣＴユニット１００を制御して、ステップＳ１において設定されたスキャン条件でＯＣＴスキャンを実行させる。

　具体的には、主制御部２１１は、光スキャナー４２を制御することにより、光源ユニット１０１から出射された光Ｌ０に基づいて生成された測定光ＬＳを偏向し、偏向された測定光ＬＳで被検眼Ｅの所定部位（例えば眼底）をスキャンさせる。測定光ＬＳのスキャンにより得られた干渉光の検出結果は、クロックＫＣに同期してサンプリングされる。干渉光の検出結果から、被検眼Ｅの３次元複素ＯＣＴデータが取得される。

（Ｓ４：収差補正）
　次に、主制御部２１１は、データ処理部２３０を制御して、ステップＳ３において取得された複素ＯＣＴデータに対して収差補正処理を実行させる。ステップＳ４では、後述するように、ＣＡＯフィルタ（パラメータ）を探索し、探索により得られたＣＡＯフィルタが複素ＯＣＴデータに適用される。ステップＳ４の詳細は、後述する。

（Ｓ５：表示？）
　次に、主制御部２１１は、ステップＳ４において探索されたＣＡＯフィルタに対応した収差情報を表示部２４０Ａに表示させるか否かを判定する。例えば、主制御部２１１は、ユーザによる操作部２４０Ｂに対する操作内容に基づいて、収差情報を表示部２４０Ａに表示させるか否かを判定する。例えば、主制御部２１１は、あらかじめ指定された動作モードに応じて、収差情報を表示部２４０Ａに表示させるか否かを判定する。

　ステップＳ５において収差情報を表示部２４０Ａに表示させると判定されたとき（ステップＳ５：Ｙ）、眼科装置１５００ｂの動作はステップＳ６に移行する。ステップＳ５において収差情報を表示部２４０Ａに表示させないと判定されたとき（ステップＳ５：Ｎ）、眼科装置１５００ｂの動作はステップＳ７に移行する。

（Ｓ６：収差情報を表示）
　ステップＳ５において収差情報を表示部２４０Ａに表示させると判定されたとき（ステップＳ５：Ｙ）、主制御部２１１は、ステップＳ４において探索されたＣＡＯフィルタのパラメータに対応した被検眼Ｅの瞳面における収差情報を表示部２４０Ａに表示させる。収差情報は、被検眼Ｅの瞳面における光の位相ずれの分布マップ（位相分布マップ）を含む。このとき、主制御部２１１は、ゼルニケ係数や、探索されたＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成された断層像及びプロジェクション画像を分布マップと同一画面に表示させてもよい。

（Ｓ７：ＯＣＴ画像を形成）
　ステップＳ６に続いて、又はステップＳ５において収差情報を表示部２４０Ａに表示させないと判定されたとき（ステップＳ５：Ｎ）、主制御部２１１は、画像形成部２２０又はデータ処理部２３０を制御して、ステップＳ４において探索されたＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成されたＯＣＴ画像を形成させる。ＯＣＴ画像には、断層像、プロジェクション画像、ｅｎ－ｆａｃｅ画像、ＯＣＴＡ画像などがある。

　いくつかの実施形態では、ユーザによる操作部２４０Ｂに対する操作内容に基づいて、ステップＳ７の処理が開始される。

　以上で、眼科装置１５００ｂの動作は終了である（エンド）。

　図１５のステップＳ４の処理は、図１６に示すように実行される。

（Ｓ１１：位相安定化）
　ステップＳ４において被検眼Ｅの複素ＯＣＴデータが取得されると、主制御部２１１は、位相安定化部２３１を制御して、取得された複素ＯＣＴデータに対して位相安定化処理を実行させる。

　位相安定化部２３１は、上記のように、被検眼Ｅの３次元複素ＯＣＴデータに対して位相安定化処理を実行する。ステップＳ１１の詳細は、後述する。

（Ｓ１２：フーリエ変換）
　続いて、主制御部２１１は、ステップＳ１１において位相安定化処理が行われた３次元複素ＯＣＴデータに対して公知のフーリエ変換処理を施し、周波数領域のＯＣＴデータに変換する。

（Ｓ１３：収差評価）
　次に、主制御部２１１は、データ処理部２３０（フィルタ情報探索部３００）を制御して、ＣＡＯフィルタを探索するための収差評価を実行させる。

　フィルタ情報探索部３００は、ＣＡＯフィルタのパラメータを更新しつつ、ＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価することで、最適なＣＡＯフィルタのパラメータを特定する。ステップＳ１３の詳細は、後述する。

（Ｓ１４：収差補正）
　次に、主制御部２１１は、収差補正部２３４を制御して、ステップＳ１３において探索されたＣＡＯフィルタを図１５のステップＳ３において取得された被検眼Ｅの３次元複素ＯＣＴデータに適用することで、３次元複素ＯＣＴデータの収差を補正する。

（Ｓ１５：逆フーリエ変換）
　続いて、主制御部２１１は、ステップＳ１４において収差が補正された３次元複素ＯＣＴデータに対して公知の逆フーリエ変換処理を施し、空間領域のＯＣＴデータに変換する。

　以上で、図１５のステップＳ４の処理は終了である（エンド）。

　図１６のステップＳ１１は、図１７に示すように実行される。

（Ｓ１１－１：参照データをセット）
　まず、位相安定化部２３１は、３次元複素ＯＣＴデータに基づいて形成される所定のＢフレームのデータを参照データとしてセットする。例えば、参照データとして、３次元複素ＯＣＴデータにおける上端のＢフレームのデータＢ_１が抽出される（図１２Ａ参照）。

（Ｓ１１－２：ターゲットデータをセット）
　次に、位相安定化部２３１は、３次元複素ＯＣＴデータに基づいて形成される所定のＢフレームのデータから抽出された参照データＢ_１にｙ方向に隣接するＢフレームのデータＢ_２をターゲットデータとしてセットする。

（Ｓ１１－３：位相差プロファイルを生成）
　続いて、位相差プロファイル生成部２３１Ａは、隣接する２つのＢフレームのデータＢ_１、Ｂ_２間のＡラインの深さ位置ごとに位相差を算出する。位相差プロファイル生成部２３１Ａは、算出された位相差をｚ方向に平均化処理を行うことにより得られたＡラインごとの位相差をＢライン方向に配列した位相差プロファイルを生成する。

（Ｓ１１－４：位相差データを除去）
　続いて、位相差データ除去部２３１Ｃは、ステップＳ１１－３において生成された位相差プロファイルから信頼性の低い位相差データを除去する。

　例えば、位相差データ除去部２３１Ｃは、ステップＳ１１－３において生成された位相差プロファイルから式（２）又は式（３）を満たす位相差データを除去する。

（Ｓ１１－５：位相ドリフトを抽出）
　続いて、位相ドリフト抽出部２３１Ｂは、ステップＳ１１－４において信頼性の低い位相差データが除去された位相差プロファイルに対して、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つを実行して位相ドリフト抽出する。

（Ｓ１１－６：位相補正）
　続いて、位相補正部２３１Ｄは、ステップＳ１１－５において抽出された位相ドリフトに基づいて、上記のように、ターゲットデータの位相を補正する。

（Ｓ１１－７：終了？）
　次に、位相安定化部２３１（又は主制御部２１１）は、位相安定化処理を終了するか否かを判定する。例えば、位相安定化部２３１は、３次元複素ＯＣＴデータにおけるすべてのＢフレームのデータに対して位相安定化処理を実行したか否かを判定することにより、位相安定化処理を終了するか否かを判定する。

　ステップＳ１１－７において位相安定化処理を終了すると判定されたとき（ステップＳ１１－７：Ｙ）、ステップＳ１１の動作は終了である（エンド）。ステップＳ１１－７において位相安定化処理を終了しないと判定されたとき（ステップＳ１１－７：Ｎ）、ステップＳ１１の動作はステップＳ１１－８に移行する。

（Ｓ１１－８：ターゲットデータを参照データにセット）
　ステップＳ１１－７において位相安定化処理を終了しないと判定されたとき（ステップＳ１１－７：Ｎ）、位相安定化部２３１は、ターゲットデータを参照データにセットする。その後、位相安定化部２３１の動作は、ステップＳ１１－２に移行する。

　例えば、位相安定化部２３１は、ＢフレームのデータＢ_２を参照データにセットする。これにより、次のステップＳ１１－２では、ＢフレームのデータＢ_２にｙ方向に隣接するＢフレームのデータＢ_３がターゲットデータにセットされ、ＢフレームのデータＢ_２、Ｂ_３の間で位相補正が行われる。上記の位相補正処理が、３次元複素ＯＣＴデータにおけるすべてのＢフレームのデータに対して順次に実行される。

　図１８に、図１６のステップＳ１２～ステップＳ１５の動作説明図を示す。

　ステップＳ１２では、ステップＳ１１において位相安定化処理が行われた複素ＯＣＴデータＯＤ１は、周波数領域の複素ＯＣＴデータＦＤ１に変換される。ステップＳ１３及びステップＳ１４では、周波数領域において探索されたＣＡＯフィルタＣＦが複素ＯＣＴデータＦＤ１に適用され、複素ＯＣＴデータの収差が補正される。ステップＳ１５では、収差が補正された複素ＯＣＴデータが空間領域の複素ＯＣＴデータＯＤ２に変換される。

　図１６のステップＳ１３の処理は、図１９に示すように実行される。

（Ｓ２１：フィルタ情報を生成）
　主制御部２１１は、フィルタ情報生成部２３２を制御して、ＣＡＯフィルタのパラメータを生成させる。

　フィルタ情報生成部２３２は、あらかじめ決められた初期値を取得してＣＡＯフィルタのパラメータとして生成する。外部の眼科装置から被検眼Ｅの屈折力情報が取得される場合、フィルタ情報生成部２３２は、屈折力情報からＣＡＯフィルタのパラメータを生成する。過去のＣＡＯフィルタを流用する場合、フィルタ情報生成部２３２は、記憶部２１２に記憶された過去のＣＡＯフィルタのパラメータを読み出すことで、ＣＡＯフィルタのパラメータを生成する。

　続いて、主制御部２１１は、粗探索部３１０を制御してステップＳ２２～ステップＳ２４の粗探索処理を実行させた後、細密探索部３２０を制御してステップＳ２５の細密探索処理を実行させる。

（Ｓ２２：デフォーカス項の係数を探索）
　まず、主制御部２１１は、粗探索部３１０を制御して、図８のＣＡＯフィルタのデフォーカス項の係数Ｃ_ｄｅを探索させる。ステップＳ２２の詳細は、後述する。

（Ｓ２３：ＨＶ方向の非点収差項の係数を探索）
　続いて、主制御部２１１は、粗探索部３１０を制御して、図８のＣＡＯフィルタのＨＶ方向の非点収差項の係数Ｃ_{ａｓ＿ＨＶ}を探索させる。ステップＳ２３の処理は、ステップＳ２２と同様に実行される。

（Ｓ２４：対角方向の非点収差項の係数を探索）
　続いて、主制御部２１１は、粗探索部３１０を制御して、図８のＣＡＯフィルタの対角方向の非点収差項の係数Ｃ_{ａｓ＿ＤＧ}を探索させる。ステップＳ２４の処理は、ステップＳ２２と同様に実行される。

（Ｓ２５：細密探索）
　次に、主制御部２１１は、細密探索部３２０を制御して、図８のＣＡＯフィルタの上記の係数の細密探索処理を実行させる。ステップＳ２５の詳細は、後述する。

　以上で、図１６のステップＳ１３の処理は終了である（エンド）。

　図１９のステップＳ２２～ステップＳ２４のそれぞれは、図２０に示すように実行される。すなわち、ステップＳ２２におけるデフォーカス項の係数Ｃ_ｄｅの探索処理は、図２０に示すように実行され、ステップＳ２３におけるＨＶ方向の非点収差項の係数Ｃ_{ａｓ＿ＨＶ}の探索処理は、図２０に示すように実行される。また、ステップＳ２４における対角方向の非点収差項の係数Ｃ_{ａｓ＿ＤＧ}の探索処理もまた、図２０に示すように実行される。

（Ｓ３１：フィルタテスト）
　粗探索部３１０は、フィルタテスト部３１１において、現段階で係数が設定されているＣＡＯフィルタを複素ＯＣＴデータに適用し、ＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価するフィルタテストを実行する。フィルタテストでは、形成されたＯＣＴ画像の画質評価値が算出される。ステップＳ３１の詳細は、後述する。

（Ｓ３２：終了？）
　続いて、粗探索部３１０（又は主制御部２１１）は、順次に係数を更新しつつ実行されるフィルタテストを終了するか否かを判定する。例えば、粗探索部３１０は、あらかじめ決められた係数の探索範囲内であるときフィルタテストを継続すると判定し、当該探索範囲を超えたときフィルタテストを終了すると判定する。例えば、粗探索部３１０は、あらかじめ決められた係数の更新回数が閾値以下であるときフィルタテストを継続すると判定し、当該更新回数が閾値を超えたときフィルタテストを終了すると判定する。

　ステップＳ３２においてフィルタテストを終了すると判定されたとき（ステップＳ３２：Ｙ）、粗探索部３１０の動作はステップＳ３４に移行する。ステップＳ３２においてフィルタテストを終了しないと判定されたとき（ステップＳ３２：Ｎ）、粗探索部３１０の動作はステップＳ３３に移行する。

（Ｓ３３：係数を更新）
　ステップＳ３２においてフィルタテストを終了しないと判定されたとき（ステップＳ３２：Ｎ）、フィルタ情報更新部３１２は、係数を所定の増加分だけ増加、又は所定の減少分だけ減少させる。その後、粗探索部３１０の動作はステップＳ３１に移行する。

（Ｓ３４：係数を特定）
　ステップＳ３２においてフィルタテストを終了すると判定されたとき（ステップＳ３２：Ｙ）、粗探索部３１０は、ＣＡＯフィルタの係数を特定する。

　粗探索部３１０は、ステップＳ３１におけるフィルタテストを繰り返し実行するにより得られた複数の画質評価値に基づいて、最も画質がよいＯＣＴ画像に適用されたＣＡＯフィルタの係数を特定する。例えば、粗探索部３１０は、画質評価値が最大又は最小であるＯＣＴ画像に適用されたＣＡＯフィルタの係数を特定する。

　以上で、図１９のステップＳ２２～ステップＳ２４のそれぞれは終了である（エンド）。

　図２０のステップＳ３１の処理は、図２１に示すように実行される。

（Ｓ４１：フィルタ処理）
　フィルタテスト部３１１は、フィルタ処理部３１１Ａにおいて、ＣＡＯフィルタを複素ＯＣＴデータに適用する。

（Ｓ４２：逆フーリエ変換）
　続いて、フィルタテスト部３１１（データ処理部２３０）は、ステップＳ４１においてフィルタ処理が施された複素ＯＣＴデータに対して逆フーリエ変換を施し、空間領域の複素ＯＣＴデータを生成する。

（Ｓ４３：プロジェクション画像を形成）
　次に、プロジェクション画像形成部２３５は、ステップＳ４３において生成された空間領域の複素ＯＣＴデータに基づいてプロジェクション画像を形成する。

（Ｓ４４：画質評価値を算出）
　続いて、画質評価部３１１Ｂは、ステップＳ４３において形成されたプロジェクション画像の画質評価値を算出する。

　以上で、図２０のステップＳ３１の処理は、終了である（エンド）。

　図１９のステップＳ２５の処理は、図２２に示すように実行される。

（Ｓ５１：フィルタテスト）
　細密探索部３２０は、フィルタテスト部３２１において、ＣＡＯフィルタを複素ＯＣＴデータに適用する。その後、フィルタテスト部３２１は、ＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価するフィルタテストを実行する。ステップＳ５１は、ステップＳ３１（ステップＳ４１～ステップＳ４４）と同様に実行される（図２１参照）。すなわち、ステップＳ５１においても、画質評価値が算出される。

（Ｓ５２：画質評価値を保存）
　続いて、細密探索部３２０は、ステップＳ５１において算出された画質評価値を記憶部２１２等の記憶部に保存する。

（Ｓ５３：終了？）
　次に、細密探索部３２０（又は主制御部２１１）は、順次に係数を更新しつつ実行されるフィルタテストを終了するか否かを判定する。例えば、細密探索部３２０は、ステップＳ５２において保存された画質評価値に基づいて画質の改善の程度があらかじめ決められた閾値以下のとき、これ以上の画質の改善を図ることが困難であると判断してフィルタテストを終了すると判定し、画質の改善の程度が閾値を超えるとき、より一層の画質の改善を図ることが可能であると判断してフィルタテストを継続すると判定する。画質の改善の程度は、画質評価値の勾配に対応する。また、例えば、細密探索部３２０は、あらかじめ決められた係数の更新回数が閾値以下であるときフィルタテストを継続すると判定し、当該更新回数が閾値を超えたときフィルタテストを終了すると判定する。

　ステップＳ５３においてフィルタテストを終了すると判定されたとき（ステップＳ５３：Ｙ）、細密探索部３２０の動作はステップＳ５６に移行する。ステップＳ５３においてフィルタテストを終了しないと判定されたとき（ステップＳ５３：Ｎ）、細密探索部３２０の動作はステップＳ５４に移行する。

（Ｓ５４：勾配を解析）
　ステップＳ５３においてフィルタテストを終了しないと判定されたとき（ステップＳ５３：Ｎ）、細密探索部３２０は、勾配算出部３２３において、ＣＡＯフィルタの各係数を微小変化させたときの勾配を算出する。ステップＳ５４の詳細は、後述する。

（Ｓ５５：係数を更新）
　続いて、フィルタ情報更新部３２２は、ステップＳ５４において算出された勾配に応じた変化分だけＣＡＯフィルタの係数を更新する。その後、細密探索部３２０の動作はステップＳ５１に移行する。

（Ｓ５６：係数を特定）
　ステップＳ５３においてフィルタテストを終了すると判定されたとき（ステップＳ５３：Ｙ）、細密探索部３２０は、ＣＡＯフィルタの係数を特定する。

　例えば、ステップＳ５３において、画質の改善の程度があらかじめ決められた閾値以下のときフィルタテストを終了すると判定された場合、細密探索部３２０は、ステップＳ５２において保存された画質評価値を有するＯＣＴ画像に適用されたＣＡＯフィルタの係数を特定する。

　例えば、ステップＳ５３において、あらかじめ決められた係数の更新回数が閾値を超えたときフィルタテストを終了すると判定された場合、細密探索部３２０は、ステップＳ５１におけるフィルタテストを繰り返し実行するにより得られた複数の画質評価値に基づいて、最も画質がよいＯＣＴ画像に適用されたＣＡＯフィルタの係数を特定する。例えば、細密探索部３２０は、画質評価値が最大又は最小であるＯＣＴ画像に適用されたＣＡＯフィルタの係数を特定する。

　以上で、図１９のステップＳ２５の処理は、終了である（エンド）。

　図２２のステップＳ５４の処理は、図２３に示すように実行される。図２３は、ＣＡＯフィルタのデフォーカス項、ＨＶ方向の非点収差項、対角方向の非点収差項、及び球面収差項を探索する処理例を表す。なお、図２３は、図８に示すＣＡＯフィルタの複数の係数を順次に探索する処理例を表すが、図８に示すＣＡＯフィルタの複数の係数を並列に探索するように処理が実行されていてもよい。

（Ｓ６１：Ｃ_ｄｅ？）
　細密探索部３２０は、探索対象の係数がデフォーカス項の係数Ｃ_ｄｅであるか否かを判定する。

　探索対象の係数がデフォーカス項の係数Ｃ_ｄｅであると判定されたとき（ステップＳ６１：Ｙ）、細密探索部３２０は、ステップＳ６１－１～ステップＳ６１－５の処理を実行する。探索対象の係数がデフォーカス項の係数Ｃ_ｄｅではないと判定されたとき（ステップＳ６１：Ｎ）、細密探索部３２０の動作は、ステップＳ６２に移行する。

（Ｓ６１－１：係数を（Ｃ_ｄｅ＋δ）に更新）
　ステップＳ６１において探索対象の係数がデフォーカス項の係数Ｃ_ｄｅであると判定されたとき（ステップＳ６１：Ｙ）、フィルタ情報更新部３２２は、係数Ｃ_ｄｅを「Ｃ_ｄｅ＋δ」に更新する。

（Ｓ６１－２：フィルタテスト）
　続いて、フィルタテスト部３２１は、ステップＳ６１－１において係数が更新されたＣＡＯフィルタを複素ＯＣＴデータに適用し、ＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価するフィルタテストを実行する。ステップＳ６１－１は、ステップＳ３１（ステップＳ４１～ステップＳ４４）と同様に実行される（図２１参照）。ステップＳ６１－２において、画質評価値が算出される。

（Ｓ６１－３：係数を（Ｃ_ｄｅ－δ）に更新）
　続いて、フィルタ情報更新部３２２は、係数Ｃ_ｄｅを「Ｃ_ｄｅ－δ」に更新する。

（Ｓ６１－４：フィルタテスト）
　続いて、フィルタテスト部３２１は、ステップＳ６１－３において係数が更新されたＣＡＯフィルタを複素ＯＣＴデータに適用し、ＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価するフィルタテストを実行する。ステップＳ６１－３は、ステップＳ３１（ステップＳ４１～ステップＳ４４）と同様に実行される（図２１参照）。ステップＳ６１－３において、画質評価値が算出される。

（Ｓ６１－５：勾配を算出）
　続いて、勾配算出部３２３は、ステップＳ６１－２において算出された画質評価値とステップＳ６１－４において算出された画質評価値とに基づいて画質評価値の勾配を算出する。その後、細密探索部３２０の動作は、ステップＳ６１に移行する。

（Ｓ６２：Ｃ_{ａｓ＿ＨＶ}？）
　ステップＳ６１において探索対象の係数がデフォーカス項の係数Ｃ_ｄｅではないと判定されたとき（ステップＳ６１：Ｎ）、細密探索部３２０は、探索対象の係数がＨＶ方向の非点収差項の係数Ｃ_{ａｓ＿ＨＶ}であるか否かを判定する。

　ステップＳ６２において探索対象の係数がＨＶ方向の非点収差項の係数Ｃ_{ａｓ＿ＨＶ}であると判定されたとき（ステップＳ６２：Ｙ）、細密探索部３２０は、ステップＳ６２－１～ステップＳ６２－５の処理を実行する。探索対象の係数が非点収差項の係数Ｃ_{ａｓ＿ＨＶ}ではないと判定されたとき（ステップＳ６２：Ｎ）、細密探索部３２０の動作は、ステップＳ６３に移行する。

（Ｓ６２－１：係数を（Ｃ_{ａｓ＿ＨＶ}＋δ）に更新）
　ステップＳ６２において探索対象の係数がＨＶ方向の非点収差項の係数Ｃ_{ａｓ＿ＨＶ}であると判定されたとき（ステップＳ６２：Ｙ）、フィルタ情報更新部３２２は、ステップＳ６１－１と同様に、係数Ｃ_{ａｓ＿ＨＶ}を「Ｃ_{ａｓ＿ＨＶ}＋δ」に更新する。ステップＳ６２－１における「δ」は、ステップＳ６１－１における「δ」と異なってもよい。

（Ｓ６２－２：フィルタテスト）
　続いて、フィルタテスト部３２１は、ステップＳ６１－１と同様に、ステップＳ６２－１において係数が更新されたＣＡＯフィルタを複素ＯＣＴデータに適用し、ＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価するフィルタテストを実行する（図２１参照）。

（Ｓ６２－３：係数を（Ｃ_{ａｓ＿ＨＶ}－δ）に更新）
　続いて、フィルタ情報更新部３２２は、ステップＳ６１－３と同様に、係数Ｃ_{ａｓ＿ＨＶ}を「Ｃ_{ａｓ＿ＨＶ}－δ」に更新する。

（Ｓ６２－４：フィルタテスト）
　続いて、フィルタテスト部３２１は、ステップＳ６１－４と同様に、ステップＳ６２－３において係数が更新されたＣＡＯフィルタを複素ＯＣＴデータに適用し、ＣＡＯフィルタが適用された複素ＯＣＴデータに基づいて形成されたＯＣＴ画像の画質を評価するフィルタテストを実行する（図２１参照）。

（Ｓ６２－５：勾配を算出）
　続いて、勾配算出部３２３は、ステップＳ６１－５と同様に、ステップＳ６２－２において算出された画質評価値とステップＳ６２－４において算出された画質評価値とに基づいて画質評価値の勾配を算出する。その後、細密探索部３２０の動作は、ステップＳ６１に移行する。

（Ｓ６３：Ｃ_{ａｓ＿ＤＧ}？）
　ステップＳ６２において探索対象の係数がＨＶ方向の非点収差項の係数Ｃ_{ａｓ＿ＨＶ}ではないと判定されたとき（ステップＳ６２：Ｎ）、細密探索部３２０は、探索対象の係数が対角方向の非点収差項の係数Ｃ_{ａｓ＿ＤＧ}であるか否かを判定する。

　ステップＳ６３において探索対象の係数が対角方向の非点収差項の係数Ｃ_{ａｓ＿ＤＧ}であると判定されたとき（ステップＳ６３：Ｙ）、細密探索部３２０は、ステップＳ６３－１～ステップＳ６３－５の処理を実行する。探索対象の係数が対角方向の非点収差項の係数Ｃ_{ａｓ＿ＤＧ}ではないと判定されたとき（ステップＳ６３：Ｎ）、細密探索部３２０の動作は、ステップＳ６４に移行する。

（Ｓ６３－１：係数を（Ｃ_{ａｓ＿ＤＧ}＋δ）に更新）
（Ｓ６３－２：フィルタテスト）
（Ｓ６３－３：係数を（Ｃ_{ａｓ＿ＤＧ}－δ）に更新）
（Ｓ６３－４：フィルタテスト）
（Ｓ６３－５：勾配を算出）
　ステップＳ６３において探索対象の係数が対角方向の非点収差項の係数Ｃ_{ａｓ＿ＤＧ}であると判定されたとき（ステップＳ６３：Ｙ）、細密探索部３２０は、ステップＳ６１－１～ステップＳ６１－５と同様に、対角方向の非点収差項の係数Ｃ_{ａｓ＿ＤＧ}について画質評価値の勾配を算出する（ステップＳ６３－１～ステップＳ６３－５）。その後、細密探索部３２０の動作は、ステップＳ６１に移行する。ステップＳ６３－１～ステップＳ６３－５は、ステップＳ６１－１～ステップＳ６１－５と同様であるため、詳細な説明を省略する。

（Ｓ６４：Ｃ_ｓｐｈ？）
　ステップＳ６３において探索対象の係数が対角方向の非点収差項の係数Ｃ_{ａｓ＿ＤＧ}ではないと判定されたとき（ステップＳ６３：Ｎ）、細密探索部３２０は、探索対象の係数が球面収差項の係数Ｃ_ｓｐｈであるか否かを判定する。

　ステップＳ６４において探索対象の係数が球面収差項の係数Ｃ_ｓｐｈであると判定されたとき（ステップＳ６４：Ｙ）、細密探索部３２０は、ステップＳ６４－１～ステップＳ６４－５の処理を実行する。

（Ｓ６４－１：係数を（Ｃ_ｓｐｈ＋δ）に更新）
（Ｓ６４－２：フィルタテスト）
（Ｓ６４－３：係数を（Ｃ_ｓｐｈ－δ）に更新）
（Ｓ６４－４：フィルタテスト）
（Ｓ６４－５：勾配を算出）
　ステップＳ６４において探索対象の係数が球面収差項の係数Ｃ_ｓｐｈであると判定されたとき（ステップＳ６４：Ｙ）、細密探索部３２０は、ステップＳ６１－１～ステップＳ６１－５と同様に、球面収差項の係数Ｃ_ｓｐｈについて画質評価値の勾配を算出する（ステップＳ６４－１～ステップＳ６４－５）。その後、細密探索部３２０の動作は、ステップＳ６１に移行する。ステップＳ６４－１～ステップＳ６４－５は、ステップＳ６１－１～ステップＳ６１－５と同様であるため、詳細な説明を省略する。

　ステップＳ６４において、探索対象の係数が球面収差項の係数Ｃ_ｓｐｈではないと判定されたとき（ステップＳ６４：Ｎ）、細密探索部３２０の動作は、次の項の係数について同様の処理を実行する。なお、すべての項の係数が終了したとき、図２２のステップＳ５４の処理は、終了である。

［変形例］
　上記の実施形態において、被検眼ＥにおけるＯＣＴ計測領域ごと（画角ごと）に、実施形態に係るＣＡＯフィルタのパラメータ（フィルタ情報）を探索するように構成されていてもよい。この場合、被検眼ＥにおけるＯＣＴ計測領域ごとに探索されたＣＡＯフィルタのパラメータに対応した被検眼Ｅの瞳面における収差情報を表示部２４０Ａに表示させてもよい。

　いくつかの実施形態では、上記の眼科情報処理方法をコンピュータに実行させるためのプログラムが提供される。このようなプログラムを、コンピュータによって読み取り可能な任意の非一時的な記録媒体に記憶させることができる。この記録媒体としては、たとえば、半導体メモリ、光ディスク、光磁気ディスク（ＣＤ－ＲＯＭ／ＤＶＤ－ＲＡＭ／ＤＶＤ－ＲＯＭ／ＭＯ等）、磁気記憶媒体（ハードディスク／フロッピー（登録商標）ディスク／ＺＩＰ等）などを用いることが可能である。また、インターネットやＬＡＮ等のネットワークを通じてこのプログラムを送受信することも可能である。

［作用］
　実施形態に係る眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラムについて説明する。

　実施形態に係る眼科情報処理装置（眼科情報処理装置１２００、眼科情報処理部１２００ａ、１２００ｂ、データ処理部２３０）は、生成部（位相差プロファイル生成部２３１Ａ）と、抽出部（位相ドリフト抽出部２３１Ｂ）と、位相補正部（２３１Ｄ）とを含む。生成部は、被検眼（Ｅ）の複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化処理を行うことにより位相差プロファイルを生成する。抽出部は、生成部により生成された位相差プロファイルから位相ドリフトを抽出する。位相補正部は、抽出部により抽出された位相ドリフトに基づいて、上記の２つのＢフレームのいずれか一方の第１複素ＯＣＴデータの位相を補正する。

　このような態様によれば、被検眼の動きやＯＣＴ測定系内の位相不安定性に起因した位相変化の影響を受けないＯＣＴデータを簡素な処理で取得することができる。それにより、ＯＣＴデータから得られる位相情報の精度を向上させることが可能になる。

　いくつかの実施形態は、位相差プロファイルから信頼性の低い位相差データを除去する除去部（位相差データ除去部２３１Ｃ）を含む。抽出部は、除去部により信頼性の低い位相差データが除去された位相差プロファイルから局所的な位相差変化が除去された位相ドリフトを抽出する。

　このような態様によれば、局所的な位相ノイズにロバストなＯＣＴデータを取得することが可能になる。

　いくつかの実施形態では、信頼性の低い位相差データは、信号強度が所定の閾値レベル未満であるデータ、及び平均化処理において所定の深さ位置における位相差による寄与が所定の閾値を超えるデータの少なくとも１つを含む。

　このような態様によれば、局所的な位相ノイズにロバストなＯＣＴデータを、簡素な処理で取得することが可能になる。

　いくつかの実施形態では、抽出部は、位相差プロファイルに対して、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つを施すことにより位相ドリフトを抽出する。

　このような態様によれば、被検眼の動きやＯＣＴ測定系内の位相不安定性に起因した位相変化の影響を受けないＯＣＴデータを、より一層簡素な処理で取得することができる。

　いくつかの実施形態では、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つのパラメータは、複素ＯＣＴデータを取得する測定系の位相変化の周波数、位相変化の速度、及び測定系におけるＯＣＴスキャンに要する時間の少なくとも１つに基づいて設定される。

　このような態様によれば、複素ＯＣＴデータの取得環境に対応した最適な条件で位相ドリフトを抽出することが可能になる。

　いくつかの実施形態では、隣接する第１のＢフレーム及び第２のＢフレームについて、生成部により生成された位相差プロファイルに基づいて第２のＢフレームの複素ＯＣＴデータの位相を補正した後、第２のＢフレーム及び第２のＢフレームに隣接する第３のＢフレームについて、生成部により生成された位相差プロファイルに基づいて第３のＢフレームの複素ＯＣＴデータの位相を補正する。

　このような態様によれば、簡素な処理で、３次元複素ＯＣＴデータの位相を安定化させることが可能になる。

　いくつかの実施形態は、位相補正部により少なくとも１つのＢフレームの複素ＯＣＴデータが補正された被検眼の複素ＯＣＴデータに基づいて被検眼のＯＣＴ画像を形成する画像形成部（画像形成部２２０、データ処理部２３０）を含む。

　このような態様によれば、被検眼の動きやＯＣＴ測定系内の位相不安定性に起因した位相変化にロバストなＯＣＴ画像を形成することが可能になる。

　実施形態に係る眼科装置（１５００ａ、１５００ｂ）は、被検眼に対して光コヒーレンストモグラフィを実行することにより複素ＯＣＴデータを取得する測定系（ＯＣＴ光学系１１０１、ＯＣＴ光学系８）と、上記のいずれかに記載の眼科情報処理装置と、を含む。

　このような態様によれば、被検眼の動きやＯＣＴ測定系内の位相不安定性に起因した位相変化の影響を受けないＯＣＴデータを簡素な処理で取得することができる眼科装置を提供することができる。

　実施形態に係る眼科情報処理方法は、生成ステップと、抽出ステップと、位相補正ステップとを含む。生成ステップは、被検眼（Ｅ）の複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化処理を行うことにより位相差プロファイルを生成する。抽出ステップは、生成ステップにおいて生成された位相差プロファイルから位相ドリフトを抽出する。位相補正ステップは、抽出ステップにおいて抽出された位相ドリフトに基づいて、上記の２つのＢフレームのいずれか一方の第１複素ＯＣＴデータの位相を補正する。

　いくつかの実施形態は、位相差プロファイルから信頼性の低い位相差データを除去する除去ステップを含む。抽出ステップは、除去ステップにおいて信頼性の低い位相差データが除去された位相差プロファイルから局所的な位相差変化が除去された位相ドリフトを抽出する。

　いくつかの実施形態では、抽出ステップは、位相差プロファイルに対して、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つを施すことにより位相ドリフトを抽出する。

　いくつかの実施形態では、隣接する第１のＢフレーム及び第２のＢフレームについて、生成ステップにおいて生成された位相差プロファイルに基づいて第２のＢフレームの複素ＯＣＴデータの位相を補正した後、第２のＢフレーム及び第２のＢフレームに隣接する第３のＢフレームについて、生成ステップにおいて生成された位相差プロファイルに基づいて第３のＢフレームの複素ＯＣＴデータの位相を補正する。

　いくつかの実施形態は、位相補正ステップにおいて少なくとも１つのＢフレームの複素ＯＣＴデータが補正された被検眼の複素ＯＣＴデータに基づいて被検眼のＯＣＴ画像を形成する画像形成ステップを含む。

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

＜その他＞
　上記の実施形態では、スウェプトソースタイプのＯＣＴを用いる場合について説明したが、スペクトラルドメインタイプのＯＣＴを用いてもよい。この場合、光源ユニット１０１において波長掃引光源に代えて低コヒーレンス光源（例えば、ＳＬＤ光源など）が用いられ、干渉光学系において検出器１２５に代えて分光器及び撮像素子（例えば、ＣＣＤなど）が用いられる。

　本発明のいくつかの実施形態を説明したが、これら実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると共に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

８、１１０１　ＯＣＴ光学系
１００　ＯＣＴユニット
２１０　制御部
２１１　主制御部
２１２　記憶部
２２０　画像形成部
２３０　データ処理部
２３１　位相安定化部
２３１Ａ　位相差プロファイル生成部
２３１Ｂ　位相ドリフト抽出部
２３１Ｃ　位相差データ除去部
２３１Ｄ　位相補正部
２３４　収差補正部
３００　フィルタ情報探索部
１０００、１０００ａ、１０００ｂ　眼科システム
１１１１　屈折力測定光学系
１２００　眼科情報処理装置
１２００ａ、１２００ｂ　眼科情報処理部
１５００ａ、１５００ｂ　眼科装置
Ｅ　被検眼

Claims

　被検眼の複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化処理を行うことにより位相差プロファイルを生成する生成部と、
　前記生成部により生成された前記位相差プロファイルから位相ドリフトを抽出する抽出部と、
　前記抽出部により抽出された前記位相ドリフトに基づいて、前記２つのＢフレームのいずれか一方の第１複素ＯＣＴデータの位相を補正する位相補正部と、
　を含む、眼科情報処理装置。
　前記位相差プロファイルから信頼性の低い位相差データを除去する除去部を含み、
　前記抽出部は、前記除去部により前記信頼性の低い位相差データが除去された前記位相差プロファイルから局所的な位相差変化が除去された前記位相ドリフトを抽出する
　ことを特徴とする請求項１に記載の眼科情報処理装置。
　前記信頼性の低い位相差データは、信号強度が所定の閾値レベル未満であるデータ、及び前記平均化処理において所定の深さ位置における位相差による寄与が所定の閾値を超えるデータの少なくとも１つを含む
　ことを特徴とする請求項２に記載の眼科情報処理装置。
　前記抽出部は、前記位相差プロファイルに対して、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つを施すことにより前記位相ドリフトを抽出する
　ことを特徴とする請求項１～請求項３のいずれか一項に記載の眼科情報処理装置。
　前記スムージング処理、前記多項式フィッティング、及び前記高周波カットフィルタ処理の少なくとも１つのパラメータは、前記複素ＯＣＴデータを取得する測定系の位相変化の周波数、前記位相変化の速度、及び前記測定系におけるＯＣＴスキャンに要する時間の少なくとも１つに基づいて設定される
　ことを特徴とする請求項４に記載の眼科情報処理装置。
　隣接する第１のＢフレーム及び第２のＢフレームについて、前記生成部により生成された位相差プロファイルに基づいて前記第２のＢフレームの複素ＯＣＴデータの位相を補正した後、
　前記第２のＢフレーム及び前記第２のＢフレームに隣接する第３のＢフレームについて、前記生成部により生成された位相差プロファイルに基づいて前記第３のＢフレームの複素ＯＣＴデータの位相を補正する
　ことを特徴とする請求項１～請求項５のいずれか一項に記載の眼科情報処理装置。
　前記位相補正部により少なくとも１つのＢフレームの複素ＯＣＴデータが補正された前記被検眼の複素ＯＣＴデータに基づいて前記被検眼のＯＣＴ画像を形成する画像形成部を含む
　ことを特徴とする請求項１～請求項６のいずれか一項に記載の眼科情報処理装置。
　前記被検眼に対して光コヒーレンストモグラフィを実行することにより複素ＯＣＴデータを取得する測定系と、
　請求項１～請求項７のいずれか一項に記載の眼科情報処理装置と、
　を含む、眼科装置。
　被検眼の複素ＯＣＴデータの隣接する２つのＢフレーム間のＡラインの深さ位置ごとに算出して得られた位相差を深さ方向に平均化処理を行うことにより位相差プロファイルを生成する生成ステップと、
　前記生成ステップにおいて生成された前記位相差プロファイルから位相ドリフトを抽出する抽出ステップと、
　前記抽出ステップにおいて抽出された前記位相ドリフトに基づいて、前記２つのＢフレームのいずれか一方の第１複素ＯＣＴデータの位相を補正する位相補正ステップと、
　を含む、眼科情報処理方法。
　前記位相差プロファイルから信頼性の低い位相差データを除去する除去ステップを含み、
　前記抽出ステップは、前記除去ステップにおいて前記信頼性の低い位相差データが除去された前記位相差プロファイルから局所的な位相差変化が除去された前記位相ドリフトを抽出する
　ことを特徴とする請求項９に記載の眼科情報処理方法。
　前記信頼性の低い位相差データは、信号強度が所定の閾値レベル未満であるデータ、及び前記平均化処理において所定の深さ位置における位相差による寄与が所定の閾値を超えるデータの少なくとも１つを含む
　ことを特徴とする請求項１０に記載の眼科情報処理方法。
　前記抽出ステップは、前記位相差プロファイルに対して、スムージング処理、多項式フィッティング、及び高周波カットフィルタ処理の少なくとも１つを施すことにより前記位相ドリフトを抽出する
　ことを特徴とする請求項９～請求項１１のいずれか一項に記載の眼科情報処理方法。
　前記スムージング処理、前記多項式フィッティング、及び前記高周波カットフィルタ処理の少なくとも１つのパラメータは、前記複素ＯＣＴデータを取得する測定系の位相変化の周波数、前記位相変化の速度、及び前記測定系におけるＯＣＴスキャンに要する時間の少なくとも１つに基づいて設定される
　ことを特徴とする請求項１２に記載の眼科情報処理方法。
　隣接する第１のＢフレーム及び第２のＢフレームについて、前記生成ステップにおいて生成された位相差プロファイルに基づいて前記第２のＢフレームの複素ＯＣＴデータの位相を補正した後、
　前記第２のＢフレーム及び前記第２のＢフレームに隣接する第３のＢフレームについて、前記生成ステップにおいて生成された位相差プロファイルに基づいて前記第３のＢフレームの複素ＯＣＴデータの位相を補正する
　ことを特徴とする請求項９～請求項１３のいずれか一項に記載の眼科情報処理方法。
　前記位相補正ステップにおいて少なくとも１つのＢフレームの複素ＯＣＴデータが補正された前記被検眼の複素ＯＣＴデータに基づいて前記被検眼のＯＣＴ画像を形成する画像形成ステップを含む
　ことを特徴とする請求項９～請求項１４のいずれか一項に記載の眼科情報処理方法。
　コンピュータに、請求項９～請求項１５のいずれか一項に記載の眼科情報処理方法の各ステップを実行させることを特徴とするプログラム。