(57)【要約】【課題】  共振スキャナによる測定光の網膜上スキャンより画像を得る構成において、画像の歪み等の画質劣化を抑制する。【解決手段】  眼底を往復走査された測定光による戻り光を電気信号に変換する光電変換手段と、一往復走査する間に走査位置を特定するための基準信号を抽出する抽出部と、測定光の往復走査による画像データを生成する画像生成手段と、基準信号を往復走査時での戻り光をサンプリングする際のサンプリング基準位置とし、且つ画像データ生成時のスキャナの走査振幅を基準振幅とする手段と、サンプリング基準位置及び基準振幅の少なくとも何れかの補正値を入力する補正値入力手段と、補正値入力手段により入力された補正値に基づいてサンプリング基準位置を補正する第一補正手段及び基準振幅を補正する第二補正手段と、を配し、画像生成手段は補正されたサンプリング基準位置及び基準振幅に基づいて画像データを再生成する構成とする。【選択図】  図１
2015139528.tif 000002
【特許請求の範囲】
【請求項１】
スキャナにより被検体を往復走査された測定光による戻り光を受光して前記戻り光を電気信号に変換する光電変換手段と、
スキャナが一往復走査する間に、前記スキャナの走査位置を特定するための基準信号を抽出する抽出部と、
前記電気信号から前記測定光の往復走査による画像データを生成する画像生成手段と、
前記基準信号を、前記往復走査時での前記戻り光をサンプリングする際のサンプリング基準位置とし、且つ前記画像データ生成時の前記往復走査の前記スキャナの走査振幅を基準振幅とする手段と、
前記サンプリング基準位置及び前記基準振幅の少なくとも何れかの補正値を入力する補正値入力手段と、
前記補正値入力手段により入力された補正値に基づいて前記サンプリング基準位置を補正する第一補正手段及び前記基準振幅を補正する第二補正手段と、を有し、
前記画像生成手段は補正された前記サンプリング基準位置及び前記基準振幅に基づいて前記画像データを再生成することを特徴とする画像生成装置。
【請求項２】
前記補正値入力手段は、
前記スキャナの往路走査による画像データと復路走査による画像データとを比較することにより、前記画像データによる画像の画質を評価する第一評価手段と、
前記第一評価手段の評価結果に応じた補正値を前記第一補正手段に入力する手段と、
前記第一補正手段により補正されたサンプリング基準位置に基づいて前記画像生成手段が再生成した画像データから、連続するスキャナ一往復走査による画像同士を比較することにより前記画像データによる画像の画質を評価する第二評価手段と、
前記第二評価手段の評価結果に基づき、前記戻り光をサンプリングする際の前記スキャナの走査振幅についての基準振幅を補正する補正値を第二補正手段に入力する手段と、を有することを特徴とする請求項１に記載の画像生成装置。
【請求項３】
前記補正値入力手段は、ユーザーからの入力により前記サンプリング基準位置及び前記基準振幅の補正パラメータを取得するスキャナ補正パラメータ取得部を有し、
前記第一補正手段及び第二補正手段は前記補正パラメータに基づき、前記サンプリング基準位置及び前記基準振幅を補正することを特徴とする請求項１に記載の画像生成装置。
【請求項４】
前記第一評価手段は、前記往復走査における往路走査における前記画像データと復路走査における前記画像データについて、隣接するそれぞれの画素同士の前記画像データを用いて前記画質の評価を行うことを特徴とする請求項２に記載の画像生成装置。
【請求項５】
前記第二評価手段は、前記往復走査における二往復の走査で得られた前記画像データのうち、一往復目の復路走査における前記画像データと二往復目の往路走査における前記画像データについて、隣接するそれぞれの画素同士の前記画像データを用いて前記画質の評価を行うことを特徴とする請求項２又は４に記載の画像生成装置。
【請求項６】
前記第一評価手段及び第二評価手段は、評価する二つの画素列の前記画像データについて、隣接するそれぞれの画素のうち輝度の高い画素と近傍の画素同士を用いて前記画質の評価を行うことを特徴とする請求項２、４及び５の何れか一項に記載の画像生成装置。
【請求項７】
前記第一評価手段及び第二評価手段は、評価する二つの画素列の前記画像データについて、各々の画素毎に比較することにより前記画質の評価を行うことを特徴とする請求項２、４及び５の何れか一項に記載の画像生成装置。
【請求項８】
前記第一評価手段及び第二評価手段は、評価する二つの画素列の前記画像データにおける輝度の比較による前記画質の評価、評価する二つの画素列の前記画像データにおける輝度差の２乗と輝度和の２乗の商の相和による前記画質の評価、評価する二つの画素列の前記画像データにおける輝度差の絶対値による前記画質の評価、評価する二つの画素列の前記画像データにおける輝度差の２乗の相和による前記画質の評価、評価する二つの画素列の正規化相互相関による前記画質の評価、の少なくとも何れかを行うことを特徴とする請求項２、４及び５の何れか一項に記載の画像生成装置。
【請求項９】
前記第一補正手段は、前記画像データより得られた画像に対する歪みの影響が最も小さい前記サンプリング基準位置を推定し、
前記推定に際して、前記第一補正手段は、前記基準信号の間の電気信号に対して前記サンプリング基準位置を決定するように前記サンプリング基準値を補正することを特徴とする請求項２及び４乃至８の何れか一項に記載の画像生成装置。
【請求項１０】
前記第二補正手段は、前記画像データより得られた画像に対する歪みの影響が最も小さい前記基準振幅を推定し、
前記推定に際して、前記第二補正手段は、前記基準信号間の電気信号に対して前記基準振幅を決定するように前記基準振幅を補正することを特徴とする請求項２及び４乃至９の何れか一項に記載の画像生成装置。
【請求項１１】
前記スキャナは共振スキャナを有し、
前記画像生成手段は、前記共振スキャナの走査に対してsin関数を用いた近似補正を行うことを特徴とする請求項１乃至１０の何れか一項に記載の画像生成装置。
【請求項１２】
前記被検体が被検眼であり、
前記被検眼の収差を測定し、前記被検眼にて発生する前記測定光及び前記戻り光における前記収差を補正する補償光学系を有する請求項１乃至１１の何れか一項に記載の画像生成装置。
【請求項１３】
前記被検体が被検眼であり、
前記被検眼の動きを追尾して撮像位置を補正する手段を有する請求項１乃至１２の何れか一項に記載の画像生成装置。
【請求項１４】
請求項１乃至１３の何れか一項に記載の画像生成装置における各構成の操作をコンピュータで実現させるためのプログラム。
【請求項１５】
請求項１４に記載のプログラムを記憶した記憶媒体。
【請求項１６】
スキャナにより被検体を往復走査された測定光による戻り光を受光し電気信号に変換する光電変換工程と、
スキャナが一往復走査する間に、前記スキャナの走査位置を特定するための基準信号を抽出する抽出工程と、
前記電気信号から往復走査による画像データを生成する画像生成工程と、
前記基準信号を、前記往復走査時での前記戻り光をサンプリングする際のサンプリング基準位置とし、且つ前記画像データ生成時の前記往復走査の前記スキャナの走査振幅を基準振幅とする工程と、
前記サンプリング基準位置及び前記基準振幅の少なくとも何れかの補正値を入力する補正値入力工程と、
前記補正値入力手段により入力された補正値に基づいて前記サンプリング基準位置及び前記基準振幅を補正する工程と、
補正された前記サンプリング基準位置及び前記基準振幅に基づいて前記画像データを再生成することを特徴とする画像生成方法。
【請求項１７】
前記補正値入力工程は、
前記スキャナの往路走査による画像データと復路走査による画像データとを比較することにより、前記画像データによる画像の画質を評価する第一評価工程と、
前記第一評価工程により得られた評価結果に応じた補正値を得る工程と、
前記補正値により補正されたサンプリング基準位置に基づいて再生成した画像データから、連続するスキャナ一往復走査による画像同士を比較することにより前記画像データによる画像の画質を評価する第二評価工程と、
前記第二評価工程により得られた評価結果に基づき、前記戻り光をサンプリングする際の前記スキャナの走査振幅についての基準振幅を補正する補正値を得る工程と、を有することを特徴とする請求項１６に記載の画像生成方法。
【請求項１８】
前記補正値入力工程は、ユーザーからの入力により前記サンプリング基準位置及び前記基準振幅の補正パラメータを取得するスキャナ補正パラメータ取得工程を有し、
前記サンプリング基準位置を補正する工程及び前記基準振幅を補正する工程は前記補正パラメータに基づき、行われることを特徴とする請求項１６に記載の画像生成方法。
【請求項１９】
請求項１６乃至１８の何れか一項に記載の画像生成方法の各工程をコンピュータに実行させることを特徴とするプログラム。
【請求項２０】
請求項１９に記載のプログラムを記憶した記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像生成装置及び画像生成方法に関し、特に眼科診療等に用いられる画像生成装置及び画像生成方法に関する。
【背景技術】
【０００２】
生活習慣病や失明原因の上位を占める疾病の早期診断を目的として、眼底部の検査が広く行われている。眼底部の検査には、例えば共焦点レーザー顕微鏡の原理を利用した眼科装置である走査型レーザー検眼鏡（ＳＬＯ：Ｓｃａｎｎｉｎｇ  Ｌａｓｅｒ  Ｏｐｈｔｈａｌｍｏｓｃｏｐｅ）が用いられる。ＳＬＯは、測定光であるレーザーを眼底に対してラスタースキャンを行い、その戻り光の強度から平面画像を高分解能かつ高速に得る装置である。
【０００３】
近年、被検眼の収差を波面センサでリアルタイムに測定し、被検眼にて発生する測定光やその戻り光の収差を波面補正デバイスで補正する補償光学系を有する補償光学ＳＬＯ（ＡＯ−ＳＬＯ：ａｄａｐｔｉｖｅ  Ｏｐｔｉｃｓ  ＳＬＯ）が開発されている。これにより、高横分解能な平面画像（以下ＡＯ−ＳＬＯ像という場合がある）の取得を可能にしている。
【０００４】
眼底画像を得る装置は、共振スキャナやガルバノスキャナを用いて測定光を２次元走査し、撮影の高速化のためスキャナの往路方向と復路方向の双方の測定光の走査時で測定を行っている。またそのスキャナの制御回路は、光検出器による検出信号を画像に変換して表示する。このために、該制御回路は、スキャナの走査に同期して画像生成用の同期信号を発生させるようにしている。このスキャナ走査位置を示す同期信号と光検出による電気的サンプリング位置をあわせることで、眼底画像を形成している。
【０００５】
ここで、スキャナの駆動波形と画像生成用の同期信号とは同じタイミングで出力される場合が多い。しかし、実際には共振スキャナの周波数、振幅等のばらつきや電気回路系の時間遅れなどが存在するので、同期信号から得られる光学的走査位置と電気的サンプリング位置が必ずしも一致しないことがある。さらに補償光学ＳＬＯは高分解能であるため、そのスキャナ位置の誤差の影響を強く受け、同期信号に基づいてサンプリングされた画像間に位置ずれが目立ってしまい、画像の歪みなどの画質劣化が発生する可能性があった。
【０００６】
このような課題を解決するために、特許文献１に開示される発明では、スキャナがスキャンしている正確な位置を専用のハードウェアを用いて取得する、或いは該ハードウェアを用いてスキャナを制御している。また特許文献２では、ピクセルクロックと呼ばれる手法が開示されている。当該手法として、画像のピクセルに対応する位置のシグナルのみを取得するように制御しているが、この手法でもスキャナ位置を正確に検出するためのハードウェアの構成が必要である。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０００−０３９５６０号公報
【特許文献２】特開２０１２−２５５９７８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
上記の特別なハードウェアの構成を追加する解決法では、装置を作成する際のコスト高となる可能性がある。また、装置の大きさに影響を与えたり、設計に制限を生じたり、部品点数の増加により組み立て時の工数にも影響する。
【０００９】
本発明は上記課題に鑑み、特別なハードウェア構成を用いずに、スキャナの特性による画像歪みが補正された画像を取得することを可能とする画像生成装置及び画像生成方法を提供する。
【課題を解決するための手段】
【００１０】
前記課題を解決するため、本発明による画像生成装置は、
スキャナにより被検眼の眼底を往復走査された測定光による戻り光を受光して前記戻り光を電気信号に変換する光電変換手段と、
スキャナが一往復走査する間に、前記スキャナの走査位置を特定するための基準信号を抽出する抽出部と、
前記電気信号から前記測定光の往復走査による画像データを生成する画像生成手段と、
前記基準信号を、前記往復走査時での前記戻り光をサンプリングする際のサンプリング基準位置とし、且つ前記画像データ生成時の前記往復走査の前記スキャナの走査振幅を基準振幅とする手段と、
前記サンプリング基準位置及び前記基準振幅の少なくとも何れかの補正値を入力する補正値入力手段と、
前記補正値入力手段により入力された補正値に基づいて前記サンプリング基準位置を補正する第一補正手段及び前記基準振幅を補正する第二補正手段と、を有し、
前記画像生成手段は補正された前記サンプリング基準位置及び前記基準振幅に基づいて前記画像データを再生成する。
【発明の効果】
【００１１】
本発明によれば、共振スキャナの特性による画像歪みが補正された画像を取得することができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の実施例１に係る画像生成装置の機能構成を示す図である。
【図２】実施例１に係る補償光学ＳＬＯの構成を示す模式図である。
【図３】実施例１に係る画像生成装置の処理手順を示すフローチャートである。
【図４】図３に示す処理手順中のスキャナ位置推定時の操作を詳細に説明するフローチャートである。
【図５】図３に示す処理手順中のスキャナ振幅推定時の操作を詳細に説明するフローチャートである。
【図６】ガルバノスキャナ及び共振スキャナのトリガ信号と反射信号との関係を示す模式図である。
【図７】sine補正を説明する図である。
【図８】画像開始点のずれの画像への影響を示す模式図である。
【図９】図５の共振スキャナ基準振幅の探索方法を説明する図である。
【図１０】本発明の実施例２に係る画像生成装置の機能構成を示す図である。
【図１１】実施例２に係る画像生成装置の処理手順を示すフローチャートである。
【図１２】実施例２に係るユーティリティアプリケーションのＵＩ画面の例を示した図である。
【発明を実施するための形態】
【００１３】
（実施例１)
本実施例では、補償光学ＳＬＯにより被検体である網膜を撮影した画像を取得する際に行われる、共振スキャナとガルバノスキャナの組み合わせによって取得したシグナル値から画像を再構成する処理について説明する。具体的には、画像横方向の一枚分に相当するガルバノスキャナのトリガ信号群と、画像縦方向の往復２列分に相当する共振スキャナのトリガ信号群を取得する。このトリガ信号発信のタイミングから推定されるスキャナ位置をサンプリング基準位置として用いて、画像を生成する。そして、画像中の各一往復画像間のずれ量が最も小さくなるサンプリング基準位置をそれぞれ算出し、共振スキャナの周波数のばらつきや電気回路系の時間遅れによるサンプリング基準位置のずれを補正する。
【００１４】
次に、求めたサンプリング基準位置に基づいて画像をさらに再構成し、画像中の隣り合う２列の画像間のずれ量が最も小さくなるサンプリング基準振幅をそれぞれ算出し、共振スキャナの振幅のばらつきを補正する。そして、算出したサンプリング基準位置と振幅から共振スキャナのトリガ信号発信時の実際のスキャナ位置とスキャナ各走査時の振幅を求める。こうして補正したパラメータを用いて画像を再構成することにより、特別なハードウェア的機能の追加を行わなくても、共振スキャナによる歪みが補正された網膜の画像を高速で取得することが可能となる。
【００１５】
＜画像生成装置の構成＞
図１は、本実施例に係る画像生成装置１０の機能構成を示したものであり、図２にその概略構成が例示される補償光学ＳＬＯ装置と接続されている。
【００１６】
図１において、信号取得部１００は、補償光学ＳＬＯ装置から共振スキャナ及びガルバノスキャナのそれぞれのトリガ信号と、網膜からの反射信号とを取得する。ここで、トリガ信号とは、共振スキャナ及びガルバノスキャナが特定の位置（基準位置）にあることを示す信号である。取得した情報は制御部１２０を通じて記憶部１３０に記憶される。画像生成部１４０は、トリガ抽出部１４１、スキャナ位置推定部１４２、スキャナ振幅推定部１４３、及び画像再構成部１４４を含む。本実施例における画像生成部１４０は、画像生成手段に相当し、光センサ２１１より得られた電気信号から測定光の往復走査による画像データを生成する。
【００１７】
トリガ抽出部１４１は、取得した信号よりガルバノスキャナ及び共振スキャナのトリガ信号を抽出し、画像生成部１４０は、抽出したトリガ信号を用いて反射信号を画像化する。本実施例におけるトリガ抽出部１４１は、スキャナが一往復走査する間に該スキャナの走査位置を特定するための基準信号を抽出する抽出部に相当する。スキャナ位置推定部１４２は、生成した画像を用いて、スキャン時の往路にて得られる画像と復路にて得られる画像間での画像ずれ量を算出する。スキャナ振幅推定部１４３は、スキャナ基準振幅を算出し、各往復の走査ごとのトリガ信号の発信時の共振スキャナの実際の位置とのずれと実際に走査した振幅とを求める。そして、得られたスキャナ位置と振幅に基づき、画像再構成部１４４が網膜の画像を再構成する。出力部１５０は、以上の処理を経て生成した画像を不図示のモニタ等に出力する。
【００１８】
次に、図２を用いて本実施例において用いた補償光学ＳＬＯの構成を説明する。補償光学ＳＬＯは、ＳＬＤ２０１、シャックハルトマン波面センサ２０６、補償光学系２０４、ビームスプリッタ（２０２、２０３）、Ｘ−Ｙ走査ミラー２０５、フォーカスレンズ２０９、絞り２１０、及び光センサ２１１を有しており、図１の画像形成装置１０と接続されている。
【００１９】
光源であるＳＬＤ（Ｓｕｐｅｒ  Ｌｕｍｉｎｅｓｃｅｎｔ  Ｄｉｏｄｅ）２０１から照射された光は眼底で反射される。反射された戻り光は、その一部が第二のビームスプリッタ２０３経由でシャックハルトマン波面センサ２０６へ、それ以外は第一のビームスプリッタ２０２経由で光センサ２１１へ入力される。シャックハルトマン波面センサ２０６は眼の収差を測定するためのデバイスであり、レンズアレイ２０７にＣＣＤ２０８が接続されている。
【００２０】
入射光がレンズアレイ２０７を透過するとＣＣＤ２０８に輝点群が現れ、該投影された輝点の位置ずれに基づき波面収差が測定される。補償光学系２０４はシャックハルトマン波面センサ２０６で測定された波面収差に基づき、収差補正デバイス(可変形状ミラーもしくは空間光位相変調器)を駆動して収差を補正する。該収差補正された光はフォーカスレンズ２０９、絞り２１０を経由し光センサ２１１にて受光される。該光センサ２１１は、光電変換手段として、被検眼の眼底上を走査された測定光による戻り光を受光してこれを電気信号に変換する。得られた電気信号は、信号取得部１００に伝達される。以上述べたように、この補償光学系は、該被検眼にて発生する測定光及び戻り光における収差を補正する。
【００２１】
スキャナであるＸ／Ｙ走査ミラー２０５を動かすことで、測定光の眼底上の走査位置を制御できる。このスキャナの制御により、操作者が予め指定した被検体の撮影対象領域、時間（フレーム数／フレームレート）のデータを取得する。このＸ／Ｙ走査ミラー２０５は、主走査であるＸ方向は共振スキャナが用いられ、副走査であるＹ方向にはガルバノスキャナが用いられる。該データを画像生成装置１０へ伝送し、画像データ（動画像もしくは静止画像）を形成する際に用いられる。
【００２２】
眼底上の特定の深さ位置にフォーカスを合わせるためには、補償光学系２０４内の収差補正デバイスを用いた調整か、光学系内に不図示のフォーカス調整用レンズを設置し、該レンズを移動することによる調整かの少なくとも何れかを用いることができる。また、眼底等の画像を得る際には、固視微動と呼ばれる不随意的な眼球運動や固視不良による眼球運動、あるいは顔の動きに伴う眼の動きの影響を軽減することを要する。このための対処として、眼の動きを追尾（トラッキング）してこれを画像に反映する場合がある。眼底の動きの計測には、画像生成装置１０によって生成した眼底画像を用いて、不図示の計算部によって特徴を持った小領域の画像であるテンプレート画像を抽出し、テンプレート画像と最も似た部位を探すパターン・マッチングによって眼底の動きを計測する方法等がある。この場合、光学系内に不図示の追尾用のスキャナを設置し、不図示の制御部を追尾用のスキャナと計算部を接続する位置に設置することで、照射光の位置を追尾用のスキャナによって計測した眼底の動きに追従させることができる。以上の構成は、本実施例において、被検眼の動きを追尾して撮像位置を補正する手段に対応する。
【００２３】
＜画像処理装置の処理手順＞
次に、図３のフローチャートを参照して、本実施例の画像生成装置１０において行われる各処理についての手順を説明する。
【００２４】
＜ステップＳ３１０＞
ステップＳ３１０において、信号取得部１００は、画像生成装置１０に接続される補償光学ＳＬＯ装置から取得された少なくとも画像１枚分の信号情報を取得する。信号情報は、網膜の撮影に用いられるガルバノスキャナのトリガ信号、共振スキャナのトリガ信号、及び撮影によって取得された網膜からの反射信号の３種類である。取得された信号情報は、制御部１２０を通じて記憶部１３０に記憶される。
【００２５】
またこのとき、取得された信号情報に付随するハードウェアの制御情報を取得し、制御部１２０を通じて記憶部１３０に記憶する。ここで制御情報とは、網膜の反射信号を取得する際にサンプリング周波数やガルバノスキャナの周波数に相当するフレームレートなどである。これらの制御情報は、信号情報に付加している撮影情報ファイルに記載されている場合もあれば、信号情報のタグ情報として含まれている場合もある。
【００２６】
＜ステップＳ３２０＞
ステップＳ３２０において、トリガ抽出部１４１は、記憶部１３０に記憶されている補償光学ＳＬＯ装置により取得された信号情報から、ガルバノスキャナ及び共振スキャナのトリガ信号の位置（トリガ位置）を取得する。そして、抽出したトリガ位置を、制御部１２０を通じて記憶部１３０に記憶する。
【００２７】
ここで、トリガ信号とは、スキャナが走査を行っている際に所定の姿勢、即ち所定角度での測定光の照射に応じて該スキャナより得られる信号、或いは所定の時間間隔にてスキャナより得られて該スキャナの姿勢を推定する際に用いられる信号にあたる。トリガ位置は、後述するように時間軸上でトリガ信号を抽出した位置に対応し、反射信号をサンプリングしてピクセル値を得るための基準となるスキャナの姿勢或いは位置に対応する。
【００２８】
図６の（ａ）にガルバノスキャナのトリガ信号の模式図、（ｂ）に共振スキャナのトリガ信号の模式図、（ｃ）に反射信号の模式図を示す。ガルバノ、共振スキャナの両トリガ信号は、それぞれ図６（ａ）、（ｂ）に示すような信号であるため、トリガ信号の抽出は、閾値処理によって行う。
【００２９】
具体的には、信号強度の変化がある閾値以上となる場合に、トリガ信号として検出する。ガルバノスキャナは画像の水平方向にレーザビームを走査していて、１つのトリガ信号から次のトリガ信号が検出される間に、水平方向の走査が1回行われ、この間に画像１枚が撮影される。共振スキャナは1つのトリガ信号が検出されている間に、画像の垂直方向２ライン分（往路と復路）の走査が行われ、この間に画像１枚分の中の画像縦方向の往復２列分が撮影される。そこで、トリガ信号検出開始のタイミングを、トリガ位置、つまりサンプリングの基準位置として抽出する。
【００３０】
また図６（ｃ）の反射信号は、眼底からの反射光の強度を図２の光センサ２１１で検出した値である。図６（ｂ）のように共振スキャナからのトリガ信号が検出されている間のみこの反射光のサンプリングを実行して、得られた画像データを画像生成に用いる。
なお、ここで述べたトリガ位置の取得法には様々な方法が考えられ、本実施例で用いた手法に限定されるものではない。
【００３１】
＜ステップＳ３３０＞
ステップＳ３３０において、スキャナ位置推定部１４２は、ステップＳ３２０で取得されたトリガ位置に基づき、画像生成を行い、生成画像を用いてトリガ信号発信時のスキャナ位置を推定する。推定したスキャナ位置は、制御部１２０を通じて記憶部１３０に記憶される。
【００３２】
図４にスキャナ位置推定の詳細を説明するためのフローチャートを示す。ここで以下のスキャナ位置推定は、連続的に撮影された複数画像の一枚一枚の共振スキャナの一往復の走査毎に対して行う。また、この処理を一画像の中の往復走査分について繰り返し、複数画像に対しても同様に行う。このようにして推定されたそれぞれのスキャナ推定位置は、制御部１２０を通じて記憶部１３０に記憶される。
【００３３】
なお、本実施例では共振スキャナ往復走査の画像２ライン分ごとに、トリガ信号発信時の共振スキャナ位置を推定している。これは、補償光学ＳＬＯ装置は高横分解能であるため、そのスキャナ周波数や電気回路系の時間遅れがスキャナ位置に及ぼすわずかな誤差が、各一往復走査ごとに不均一な画像ゆがみとして影響を及ぼすためである。また、共振スキャナの走査に関わるばらつきが、装置本体内の温度や電源の不安定性などの影響も受けるため、スキャナ位置の誤差自体のばらつきも大きいことも原因として挙げられる。
【００３４】
なお、補償光学ＳＬＯ装置のスキャナ設定や取得画像の撮影条件によっては、この画像歪みがある一定の周期ごとに現れる場合もある。この場合は、画像歪みの現れる周期に相当する数の往復走査ごとにスキャナ位置の推定を行うことで推定ステップの回数を減らすことも可能である。
【００３５】
＜ステップＳ４１０＞
ステップＳ４１０において、スキャナ位置推定部１４２は、ステップＳ３１０で信号情報と共に取得された制御情報に基づき、トリガ発信時の共振スキャナ位置を、画像開始位置の設計値の初期位置ｘ0として設定する。ここで初期値の設定法にはいくつかの方法があり、例えば本装置の設計値としては、トリガ位置よりトリガーディレイ分を考慮した点からといった形で取得される。また他にも、連続して数枚の画像を撮影する場合、その撮影画像より前に撮影した画像の往復走査ごとの平均的なスキャナ推定位置を初期位置として設定することで、推定位置の探索ステップの回数を少なくするといった方法等も用いられる。
【００３６】
＜ステップＳ４２０＞
ステップＳ４２０において、スキャナ位置推定部１４２は、ステップＳ４１０で設定された画像開始位置ｘ0に基づき、共振スキャナの一往復分の画像２ライン分のsine補正を行い画像化する。共振スキャナは走査範囲内を同一速度で走査していないため、縦軸方向の位置に均一にデータを分割すると、共振スキャナの走査位置に応じた速度によって画像にゆがみが生じてしまう。このため、共振スキャナの走査軌道をsine波に近似して、各画素に対しサンプリング数に重みづけを行うことで画像歪みを解消している。
【００３７】
図７にsine補正の模式図を示す。同グラフにおいて、横軸は時間ｔとして、縦軸は位置ｘとして示し、共振スキャナ走査の軌道はステップＳ３１０で取得された周波数のsine波であると仮定し、且つその振幅を１．０としている。このsine波の縦軸を仮に１０等分して、分割された区間をそれぞれ時間０．０からｎ＝０、１、２、・・・９と名付ける。また、縦軸を分割した直線とsine波との交点の時間をそれぞれｔ（０）、ｔ（１）、・・・ｔ（９）と名付ける。この縦軸を分割した区間を共振スキャナが走査している間にサンプリングしたデータを、この位置の画素に割り当てる。すると、サンプリングは１５ＭＨｚで、時間的に等間隔に行われるが、対応する共振スキャナの１区間を走査している、つまりサンプリングしている時間は、画素により異なっていることがわかる。よって取得される反射信号の輝度値を、対応する位置のピクセル値に対応づけ、一つの画素に複数の信号が対応する場合には、画素値を複数の信号の平均値として取得する。
【００３８】
＜ステップＳ４３０＞
ステップＳ４３０において、スキャナ位置推定部１４２は、ステップＳ４２０で生成されたライン画像における各々の画素のピクセル値に基づき、指標Ｄを算出する。この指標Ｄを算出する目的は、画像開始位置のずれによって生じる画質の劣化を、定量的に評価することである。即ち、画像データによる画像の画質を評価する。また、指標Ｄの算出方法は、たとえば、共振スキャナ一往復分の隣接する画像２ラインの左右に隣接する画素同士の輝度値差を以下の式で計算したもの等がある。
【００３９】
【数１】
2015139528.tif 000003
ここでは、ライン２i−１、２i列目（ｉ＝１，２，３，，，）の隣接するラインのｊ行目の画素の輝度をIとあらわしており、それを画像２ラインすべての画素値を計算している。
【００４０】
ここで、この演算に用いるｊ行は上記２ラインのうちの隣り合う値を持つピクセルや、ある程度の範囲を指定して行っても良い。他にも輝度の高い領域の画素に対して重みづけを行う場合は以下のような式で計算する方法もある。この方法は、輝度の高い領域のみを計算してもよい。
【数２】
2015139528.tif 000004
【００４１】
また他には、輝度差の絶対値の相和や、正規化相互相関などを用いて評価する方法もある。これらの評価指標は、どのような画像のずれ方を検討するかによって最適な方法が異なると考えられるが、焦点の合った眼底画像の視細胞や血管といった輝度差の大きい構造物の多い画像の場合は、上記二つの式によるものでも良いと考えられる。
【００４２】
図８に眼底の視細胞の撮影画像の模式図を示す。図８（ａ）のように画像のサンプリング基準位置が正しく設定されている場合には横方向に隣接する画素間のずれ量はほとんどない。一方、図８（ｂ）に示すように画像サンプリング基準位置にずれが生じると、横方向に隣接する画素間のずれ量は大きくなる。この双方の指標Ｄを比較すると、図８（ａ）のほうが隣接する画素間の輝度差も小さいので、指標Ｄは小さな値になる。つまり、指標Ｄが小さいほど画素間のずれ量は小さいことがわかる。
【００４３】
＜ステップＳ４４０＞
ステップＳ４４０において、スキャナ位置推定部１４２は、ステップＳ４１０で設定されたサンプリング基準位置から、±Δｘずれたサンプリング基準位置に対する指標Ｄを算出する。具体的には、ｘ0＋Δｘ、ｘ0−Δｘのサンプリング基準位置を用いてライン画像を作成し、それぞれ指標Ｄを算出する。Δｘの大きさは画像の解像度やスキャンレートに依存するが、ここでは１サンプリング単位とした。
【００４４】
＜ステップＳ４５０＞
ステップＳ４５０において、スキャナ位置推定部１４２は、ステップＳ４４０で求められたｘ0＋Δｘ、ｘ0、ｘ0−Δｘのサンプリング基準位置に対する指標Ｄが、以下の関係を満たすか否かを判定する。ここで、サンプリング基準位置ｘ0に対応する指標の値をＤ（ｘ0）としている。
Ｄ（ｘ0）＜Ｄ（ｘ0−Δｘ）
Ｄ（ｘ0）＜Ｄ（ｘ0＋Δｘ）
【００４５】
上記の関係を満たす場合にはステップＳ４７０に、満たさない場合にはステップＳ４６０に進む。ここで上記の関係を満たす場合とは、指標Ｄが最も小さくなるサンプリング基準位置が、ｘ0−Δｘからｘ0＋Δｘの間に含まれている場合である。このような場合には、ステップＳ４７０でさらに精度の高いサンプリング基準位置を求める処理を行う。一方、指標Ｄが最も小さくさくなるサンプリング基準位置がｘ0−Δｘからｘ0＋Δｘの間に含まれていない可能性がある場合には、ステップＳ４６０でサンプリング基準位置を探索する範囲の変更を行う。
【００４６】
＜ステップＳ４６０＞
ステップＳ４６０において、スキャナ位置推定部１４２は、ライン画像の指標Ｄが最も小さくなるサンプリング基準位置を探索する範囲をΔｘだけずらす。具体的には、Ｄ（ｘ0−Δｘ）＞Ｄ（ｘ0＋Δｘ）の場合には、ｘ0＋Δｘをサンプリング基準位置の初期値としてステップＳ８２０に戻る。Ｄ（ｘ0＋Δｘ）＞Ｄ（ｘ0−Δｘ）の場合にはｘ0−Δｘをサンプリング基準位置の初期値としてステップＳ４２０に戻る。
【００４７】
＜ステップＳ４７０＞
ステップＳ４７０において、スキャナ位置推定部１４２は、ライン画像の指標Ｄが最も小さくなるサンプリング基準位置を探索する範囲を、ｘ0もしくはｘ0−Δｘ/２、ｘ0＋Δｘ/２を中心となるサンプリング基準位置とした±Δｘ/２の範囲に限定する。つまり、Δｘの値を半分にすることに相当する。ここで、中心となるサンプリング基準位置の選択は、指標Ｄが最も小さくなるサンプリング基準位置を選択する。
【００４８】
そして、Ｓ４２０〜Ｓ４６０の処理を再度行った後、Ｓ４７０にて再度Δｘの値をさらに半分にして、Ｓ４２０〜Ｓ４６０の処理を行うという処理を繰り返すことで、指標Ｄの最も小さな画像を生成するサンプリング基準位置の精度を向上させることが可能になる。具体的には、Ｍ回の繰り返し処理により、サンプリング基準位置の推定精度は±Δｘ/２Mとなる。求めたい精度に応じて繰り返し回数は決定されるが、ここではＭ＝３とする。
【００４９】
＜ステップＳ４８０＞
ステップＳ４８０において、スキャナ位置推定部１４２は、ステップＳ４７０の処理が何回繰り返されたかを判断する。繰り返し回数がＭ回未満の場合には、Δｘ＝Δｘ/２としてステップＳ４２０に戻る。繰り返し回数がＭ回に達した場合には、その時の画像開始位置を取得する。そうして抽出したサンプリング基準位置を、トリガ信号発信時の実際の共振スキャナ位置として制御部１２０を通じて記憶部１３０に記憶する。
【００５０】
以上の本実施例におけるスキャナ位置推定部１４２にて行われた操作は、第一評価手段、第一補正手段、及び該第一補正手段に補正値を入力する手段により行われ、当該スキャナ位置推定部１４２はこれら各手段における操作を実行するモジュール領域を有する。第一評価手段は、前述したように、スキャナの往路走査による画像データと復路走査による画像データとを比較することにより、画像データによる画像の画質を評価する。また、以上の繰り返し操作は、第一評価手段による評価結果に応じた補正値を第一補正手段に入力する機能を呈するスキャナ位置推定部１４２内のモジュール領域により実行される。また、本実施例では、第一評価手段は、スキャナの往復走査における往路走査における画像データと復路走査における画像データについて、隣接するそれぞれの画素同士の画像データを用いて画質の評価を行っている。
【００５１】
また、以上述べた第一補正手段による操作は、画像データより得られた画像に対する歪みの影響が最も小さいサンプリング基準位置を推定し、推定に際して、第一補正手段は、基準信号の間の電気信号に対してサンプリング基準位置を決定するようにサンプリング基準値を補正するものであるとの把握も可能である。
【００５２】
＜ステップＳ３４０＞
ステップＳ３４０において、スキャナ振幅推定部１４３は、ステップＳ３３０で取得された、補正されたサンプリング基準位置に基づき、画像生成を行い、生成画像を用いて共振スキャナ走査時のスキャナ基準振幅を推定する。そうして推定したスキャナ振幅を、制御部１２０を通じて記憶部１３０に記憶する。
【００５３】
図５にスキャナ基準振幅推定の詳細を説明するためのフローチャートを示す。ここで以下のスキャナ振幅推定は、スキャナ位置推定同様に、連続的に撮影された複数画像各々について共振スキャナの二往復の走査毎で得られた画像データに対して行う。また、当該処理を一画像の中のそれぞれの二往復の走査分を繰り返し、複数画像に対しても同様に行う。以上の処理により推定されたそれぞれのスキャナ振幅は、制御部１２０を通じて記憶部１３０に記憶される。
【００５４】
本実施例では共振スキャナの二往復の走査毎ごとに、共振スキャナ基準振幅を推定している。これは、補償光学ＳＬＯ装置は高横分解能であるため、スキャナ振幅のわずかな誤差も、各一往復走査ごとに不均一な画像ゆがみとして影響を及ぼすためである。また同様に、得られる画像の再構成には、装置本体内の温度や電源の不安定性などの影響も受けるため、スキャナ基準振幅の誤差自体のばらつきも大きいことが原因である。
【００５５】
なお、補償光学ＳＬＯ装置のスキャナ設定や取得画像の撮影条件によっては、この画像歪みが往復走査数回ごとに現れる場合もある。この場合は、画像歪みの現れる往復走査ごとにスキャナ基準振幅の推定を行うことで推定ステップの回数を減らすことも可能である。
【００５６】
＜ステップＳ５１０＞
ステップＳ５１０において、スキャナ振幅推定部１４３は、ステップＳ３１０で取得された制御情報に基づき、共振スキャナ走査時のスキャナ基準振幅を、共振スキャナの設計値の初期値y0として設定する。また、すでに対象となる共振スキャナ二往復分の画像よりも一つ前にある一往復画像についてスキャナ基準振幅の探索が既に終わっており、値が決定されている場合は、一つ前にある一往復画像の算出した振幅を初期値として設定する。
【００５７】
なお、以上基準振幅の設定は、前述したサンプリング基準位置の設定とあわせ、基準信号を往復走査時での戻り光をサンプリングする際のサンプリング基準位置とし、且つ画像データ生成時の往復走査のスキャナの走査振幅を基準振幅とする手段として機能する、画像生成部１４０内のモジュール領域により実行される。
【００５８】
＜ステップＳ５２０＞
ステップＳ５２０において、スキャナ振幅推定部１４３は、ステップＳ５１０で設定されたスキャナ基準振幅y0とステップＳ３３０で取得されたスキャナ一往復画像ごとのサンプリング基準位置に基づき、共振スキャナの二往復分の画像４ライン分をステップＳ４２０と同様にsine補正を行い画像化する。
【００５９】
＜ステップＳ５３０＞
ステップＳ５３０において、スキャナ振幅推定部１４３は、ステップＳ５２０で生成されたライン画像のピクセル値に基づき、指標Ｅを算出する。この指標Ｅを算出する目的は、振幅のずれによって生じる画質の劣化を、定量的に評価することである。また、指標Ｅの算出方法としては、たとえば、共振スキャナ二往復分のうち、一往復目の復路画像列と二往復目の往路画像列に対して、隣接する画像２ラインの左右に隣接する画素同士の輝度値差を以下の式で計算したもの等がある。
【００６０】
【数３】
2015139528.tif 000005
ここでは、ライン２i、２i＋１列目（ｉ＝１，２，３，，，）の隣接するラインのｊ行目の画素の輝度をIとあらわしており、それを画像２ラインすべての画素値を計算している。この演算に用いるｊ行は上記２ラインのうちの隣り合う値を持つピクセルや、ある程度の範囲を指定して行っても良い。
【００６１】
また他には、輝度差の絶対値の相和や、正規化相互相関などを用いて評価する方法もある。これらの評価指標は、どのような画像のずれ方を検討するかによって最適な方法が異なると考えられるが、焦点の合った眼底画像の視細胞や血管といった輝度差の大きい構造物の多い画像の場合は、上記の式によるものでも良いと考えられる。
【００６２】
＜ステップＳ５４０＞
ステップＳ５４０において、スキャナ振幅推定部１４３は、ステップＳ５１０で設定されたスキャナ基準振幅から、二往復目の往復画像のみ±Δｙずれたスキャナ基準振幅に対する指標Ｅを算出する。具体的には、ｙ0＋Δｙ、ｙ0−Δｙのスキャナ基準振幅を用いて二往復目のライン画像を作成し、それぞれ指標Ｅを算出する。図９にスキャナ基準振幅探索方法を説明する図を示す。図９には共振スキャナの二往復分の駆動走査を示す曲線が示されており、初期の振幅y0は１往復目のスキャナ走査の振幅に相当する。また、それぞれ一、二往復目の往路、復路の画像列に相当するエリアは４つの四角いエリアで示した。
【００６３】
ここで、ステップＳ５３０で説明したように、指標Ｅの算出に用いるのは、一往復目の復路画像列と二往復目の往路画像列の２ラインである。この２ラインは図の上方でつながっている、つまりそれぞれの振幅の上端は必ず等しい位置にあることがわかる。そこで２往復目の振幅の下端をそれぞれ±Δｙずれた位置にあると想定して、２往復目をそれぞれの振幅で再画像化する。そして、再画像化した２往復目の往路画像を用いて、指標Ｅを算出するのである。Δｙの大きさは画像の解像度やスキャンレートに依存するが、ここでは表示する１画素あたりの長さ単位とした。
【００６４】
＜ステップＳ５５０＞
ステップＳ５５０において、スキャナ振幅推定部１４３は、ステップＳ５４０で求められたｙ0＋Δｙ、ｙ0、ｙ0−Δｙの二往復目のスキャナ基準振幅に対する指標Ｅが、以下の関係を満たすか否かを判定する。ここで、スキャナ基準振幅ｙ0に対応する指標の値をＥ（ｙ0）としている。
Ｅ（ｙ0）＜Ｅ（ｙ0−Δｙ）
Ｅ（ｙ0）＜Ｅ（ｙ0＋Δｙ）
【００６５】
上記の関係を満たす場合にはステップＳ５７０に、満たさない場合にはステップＳ５６０に進む。ここで上記の関係を満たす場合とは、指標Ｅが最も小さくなる二往復目のスキャナ基準振幅が、ｙ0−Δｙからｙ0＋Δｙの間に含まれている場合である。このような場合には、ステップＳ５７０でさらに精度の高いスキャナ基準振幅を求める処理を行う。一方、指標Ｅが最も小さくさくなる二往復目のスキャナ基準振幅がｙ0−Δｙからｙ0＋Δｙの間に含まれていない可能性がある場合には、ステップＳ５６０で二往復目のスキャナ基準振幅を探索する範囲の変更を行う。
【００６６】
＜ステップＳ５６０＞
ステップＳ５６０において、スキャナ振幅推定部１４３は、ライン画像の指標Ｅが最も小さくなる二往復目のスキャナ基準振幅を探索する範囲をΔｙだけずらす。具体的には、Ｅ（ｙ0−Δｙ）＞Ｅ（ｙ0＋Δｙ）の場合には、ｙ0＋Δｙを二往復目のスキャナ基準振幅の初期値としてステップＳ５２０に戻る。Ｅ（ｙ0＋Δｙ）＞Ｅ（ｙ0−Δｙ）の場合にはｙ0−Δｙを二往復目のスキャナ基準振幅の初期値としてステップＳ５２０に戻る。
【００６７】
＜ステップＳ５７０＞
ステップＳ５７０において、スキャナ振幅推定部１４３は、ライン画像の指標Ｅが最も小さくなるスキャナ基準振幅を探索する範囲を、y0もしくはy0−Δy/２、y0＋Δy/２を中心となるスキャナ基準振幅とした±Δy/２の範囲に限定する。つまり、Δyの値を半分にすることに相当する。ここで、中心となるスキャナ基準振幅の選択は、指標Eが最も小さくなるスキャナ基準振幅を選択する。
【００６８】
そして、Ｓ５２０〜Ｓ５６０の処理を再度行った後、Ｓ５７０にて再度Δｙの値をさらに半分にして、Ｓ５２０〜Ｓ５６０の処理を行うという処理を繰り返すことで、指標Eの最も小さな画像を生成するスキャナ基準振幅の精度を向上させることが可能になる。具体的には、Ｍ回の繰り返し処理により、スキャナ基準振幅の推定精度は±Δｙ/２Mとなる。求めたい精度に応じて繰り返し回数は決定されるが、ここではＭ＝３とする。
【００６９】
＜ステップＳ５８０＞
ステップＳ５８０において、スキャナ振幅推定部１４３は、ステップＳ５７０の処理が何回繰り返されたかを判断する。繰り返し回数がＭ回未満の場合には、Δｙ＝Δｙ/２としてステップＳ５２０に戻る。繰り返し回数がＭ回に達した場合には、その時のスキャナ基準振幅を取得する。そうして抽出したスキャナ基準振幅を、共振スキャナ走査時の実際の共振スキャナ振幅として制御部１２０を通じて記憶部１３０に記憶する。
【００７０】
以上の本実施例におけるスキャナ振幅推定部１４３にて行われた操作は、第二評価手段、第二補正手段、及び該第二補正手段に補正値を入力する手段により行われ、当該スキャナ振幅推定部１４３はこれら各手段における処理を実行するモジュール領域を有する。第二評価手段は、前述したように、第一補正手段により補正されたサンプリング基準位置に基づいて画像生成手段が再生成した画像データから、連続するスキャナ一往復走査による画像同士を比較することにより画像データによる画像の画質を評価する。また、以上の繰り返し操作は、第二評価手段の評価結果に基づき、戻り光をサンプリングする際のスキャナの走査振幅についての基準振幅を補正する補正値を第二補正手段に入力する手段として機能するスキャナ振幅推定部１４３内のモジュール領域により実行される。
【００７１】
第二評価手段による評価結果に応じた補正値を第二補正手段に入力することにより実行される。また、本実施例では、第二評価手段は、スキャナの往復走査における二往復の走査で得られた画像データのうち、一往復目の復路走査における画像データと二往復目の往路走査における画像データについて、隣接するそれぞれの画素同士の画像データを用いて画質の評価を行っている。
【００７２】
また、本実施例において述べた以上の第一評価手段、補正値を第一補正手段に入力する手段、第二評価手段、及び補正値を第二補正手段に入力する手段は、本実施例における補正値入力手段を構成する。即ち、該補正値入力手段は、サンプリング基準位置及び基準振幅の少なくとも何れかの補正値を対応する補正手段に入力する。また、第一補正手段及び第二補正手段は、補正値入力手段により入力された補正値に基づいて、サンプリング基準位置及び基準振幅を各々補正する。
【００７３】
また、以上述べた第二補正手段の操作は、画像データより得られた画像に対する歪みの影響が最も小さい基準振幅を推定し、該推定に際して、第二補正手段は、基準信号間の電気信号に対して基準振幅を決定するように基準振幅を補正するものであるとの把握も可能である。
【００７４】
＜ステップＳ３５０＞
ステップＳ３５０において、画像再構成部１４４は、ステップＳ３３０で取得されたスキャナ一往復画像ごとのサンプリング基準位置とステップＳ３４０で取得されたスキャナ一往復画像ごとのスキャナ基準振幅に基づきsin補正を行い、画像の再構成を行う。即ち、画像生成手段は補正されたサンプリング基準位置及び基準振幅に基づいて画像データを再生成し、且つ画像の再構成を実行する。その際、画像再構成部１４４は、共振スキャナの走査に対してsin関数を用いた近似補正を行うこととなる。そして再構成した補償光学ＳＬＯ装置による網膜の画像を、制御部１２０を通じて記憶部１３０に記憶する。
【００７５】
ここで画像再構成は、連続的に撮影された複数画像の一枚一枚に対して行う。つまり、ステップＳ３２０で取得されたガルバノスキャナのトリガ信号一つ一つに相当する画像一枚分に対して上記の処理を行い、それを枚数分繰り返す。そうしてすべての画像に対して画像再構成を行った後、それらを統合した画像群を構成して、それを記憶部１３０に記憶する。
【００７６】
＜ステップＳ３６０＞
ステップＳ３６０において、出力部１５０は、ステップＳ３５０において記憶部１３０に記憶された網膜の再構成画像を、出力部１５０を通じてモニタ等に表示する。さらに、ステップＳ３１０〜ステップＳ３５０で記憶部１３０に保存されたスキャナ往復走査毎のサンプリング基準位置やスキャナ基準振幅などをデータベースに保存する。
【００７７】
以上の構成により、補償光学ＳＬＯ装置により網膜の画像を取得する際に、スキャナ位置を特定する特別なハードウェア構成を持たなくても、共振スキャナ由来の歪みを補正した画像を生成することが可能になる。
【００７８】
また上記実施例１では、サンプリング基準位置を算出後、スキャナ基準振幅を算出して、画像を生成したが、振幅の変化によって最適なサンプリング基準位置が異なる場合は、再度振幅算出後、サンプリング基準位置を算出しても良い。または、計算は複雑になるが、サンプリング基準位置とスキャナ基準振幅の二つのパラメータを同時に振って最適なそれぞれの値を算出することも可能である。また画像によってはどちらかの補正による画像歪みが小さい場合、計算の量を減らすため、１枚目のみ双方のパラメータを算出して、その後の画像は、どちらかの算出処理のみ行うことも可能である。
【００７９】
また、先に述べた第一評価手段及び第二評価手段による画質の評価は、評価する二つの画素列の画像データについて、隣接するそれぞれの画素のうち輝度の高い画素と近傍の画素同士を用いている。しかし、これを評価する二つの画素列の画像データについて、各々の画素毎に比較することにより画質の評価を行うこととしても良い。或いは、評価する二つの画素列の画像データにおける輝度の比較による画質の評価、評価する二つの画素列の画像データにおける輝度差の２乗と輝度和の２乗の商の相和による画質の評価、評価する二つの画素列の画像データにおける輝度差の絶対値による画質の評価、評価する二つの画素列の画像データにおける輝度差の２乗の相和による画質の評価、評価する二つの画素列の正規化相互相関による画質の評価、の少なくとも何れかを行うこととしても良い。
【００８０】
（実施例２）
実施例１では、サンプリング基準位置とスキャナ基準振幅をそれぞれ探索して算出することで、自動的に画像の再構成を行う処理について示した。
実施例２では、サンプリング基準位置とスキャナ基準位置を（以後、スキャナ補正パラメータと呼ぶ）ユーザーが画像をみながら変更する構成とする。
【００８１】
本実施例に関わる画像生成装置１０の機能構成を、図１０に示す。但し１００、１２０、１３０、１５０の機能構成に関しては図１と同様なので説明は省略する。本実施例では画像生成部１４０が、トリガ抽出部１４１、画像再構成部１４４のみを含み、実施例１でトリガ抽出部１４１が取得したサンプリング基準位置と記憶部１３０に記憶された制御情報のみを用いて生成した画像を用いて、スキャナ補正パラメータ取得部１０１０でユーザーが入力したスキャナ補正パラメータを用いて画像の再構成を行う。
【００８２】
本実施例では、実施例１における補正値入力手段はユーザーからの入力によりサンプリング基準位置及び基準振幅の補正パラメータを取得するスキャナ補正パラメータ取得部を有し、且つ前述した第一補正手段及び第二補正手段は補正パラメータに基づいてサンプリング基準位置及び基準振幅を補正することとなる。
【００８３】
以下、図１１のフローチャートを参照して、本実施形態の画生成理装置１０の処理手順を説明する。ここでステップＳ３１０、Ｓ３２０、Ｓ３４０に関しては実施例１で説明した処理手順と変わらないために、説明は省略する。
【００８４】
＜ステップＳ１１５０＞
ステップＳ１１５０において、出力部１５０は、ステップＳ３４０において記憶部に記憶された網膜の再構成画像と、記憶部１３０に記憶されているスキャナ補正パラメータとを、出力部１５０を通じてモニタ等に表示する。この表示の一例として図１２に、モニタ画面に表示するスキャナ補正パラメータを調整するためのユーティリティアプリケーションのＵI画面構成の例を示す。このアプリケーションは、メイン画面１２１０、スキャナ補正パラメータ情報表示１２２０、網膜の再構成された画像表示１２３０、及びスキャナ補正パラメータ変更用のコマンドボタン１２４０などを有している。
【００８５】
＜ステップＳ１１６０＞
ステップＳ１１６０において、スキャナ補正パラメータ取得部１０１０は、ステップＳ１１５０においてアプリケーションのメイン画面１２１０に表示されたスキャナ補正パラメータ１２２０及び生成された網膜の再構成画像１２３０に対して、ユーザーが行うスキャナ補正パラメータの変更を取得する。変更方法は、すべての各一往復２ラインの画像に対して同じ変化量を与える場合と、それぞれの一往復２ライン画像ごとにスキャナ補正パラメータを変更する場合などがある。図１２の例では、スキャナ補正パラメータ変更用のコマンドボタン１２４０を用いて、ユーザーの入力に従って変更を実行する。画像変更がなければ処理は終了する。変更がある場合には、ステップＳ１１７０に進む。
【００８６】
＜ステップＳ１１７０＞
ステップＳ１１７０において、スキャナ補正パラメータ取得部１０１０は、ユーザーが入力したスキャナ補正パラメータを取得し、制御部１２０を通じて記憶部１３０に記憶する。その後ステップＳ３４０に戻り、新たに取得されたスキャナ補正パラメータに基づき再度画像の再構成を行う。
【００８７】
以上の構成により、補償光学ＳＬＯ装置により取得される網膜の画像を生成する際に、共振スキャナのスキャナ補正パラメータをユーザーが画像をみながら調節し、最適なスキャナ補正パラメータを選択することが可能になる。
【００８８】
（その他の実施形態）
実施例１で紹介した方法を用いて、画像歪みを生成画像から算出したスキャナ補正パラメータを用いて画像を再生成した後に、実施例２で紹介した方法を用いて、ユーザーの入力によりスキャナ補正パラメータを選択して補正を行うことで、実施例１の評価方法で画像歪みの補正が不十分な画像や誤検出により不適切に生成された画像などに対しても、適切な歪み補正を行うことができる。
【００８９】
また、本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記憶した記憶媒体を、システムあるいは装置に供給する。そして、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読出し実行することによっても、達成されることは言うまでもない。
【符号の説明】
【００９０】
１０      画像生成装置
１００    信号取得部
１２０    制御部
１３０    記憶部
１４０    画像生成部
１４１    トリガ抽出部
１４２    スキャナ位置推定部
１４３    スキャナ振幅推定部
１４４    画像再構成部
１５０    出力部
【図１】
2015139528.tif 000006
【図２】
2015139528.tif 000007
【図３】
2015139528.tif 000008
【図４】
2015139528.tif 000009
【図５】
2015139528.tif 000010
【図６】
2015139528.tif 000011
【図７】
2015139528.tif 000012
【図９】
2015139528.tif 000013
【図１０】
2015139528.tif 000014
【図１１】
2015139528.tif 000015
【図１２】
2015139528.tif 000016
【図８】
2015139528.tif 000017