ハードウェア構成や画像補正用のチャート画像等を用いずに、スキャナの特性による画像歪みが補正された画像を取得することが可能な仕組みを提供する。【解決手段】パラメータ演算処理部141は、補償光学SLO装置200内のスキャナが被検眼の所定領域を1往復走査する間に出力された当該スキャナの基準信号を抽出し、補償光学SLO装置200内の各撮影手段ごとに、当該基準信号をサンプリング基準位置として往復走査のサンプリングデータ配列を生成し、当該往復走査のサンプリングデータ配列のうちの往路走査と復路走査とのサンプリングデータ配列を比較してこれらの相関を評価するとともに、各撮影手段ごとに得られた評価結果を比較して確からしさを評価し、これらの評価結果に基づいてサンプリング基準位置を補正する。画像構成部142は、補正されたサンプリング基準位置に基づいて、各撮影手段ごとに、画像データを構成する。請求項1乃至22のいずれか1項に記載の画像生成装置の各手段としてコンピュータを機能させるためのプログラム。各手段における処理を、複数の前記往復走査について並列に行うことを特徴とするプログラム。請求項24または25に記載のプログラムを記憶したことを特徴とするコンピュータ読み取り可能な記憶媒体。生活習慣病や失明原因の上位を占める疾病の早期診断を目的として、被検眼の眼底部の検査が広く行われている。Ophthalmoscope)がある。このSLO装置は、測定光であるレーザーを被検眼の眼底部に対してラスタースキャンを行い、その戻り光の強度から平面画像である眼底画像を高分解能かつ高速に得る装置である。SLO)装置が開発されている。この補償光学SLO装置では、高横分解能な平面画像(以下、「AO−SLO画像」という場合がある)の取得を可能にしている。【0004】さらに近年では、受光部の手前にあるピンホールの径や形状、位置を変えることにより、散乱光を取得して得られた非共焦点画像を観察する方法が用いられるようになってきている(非特許文献1参照)。この非共焦点画像では、これまで補償光学SLO装置では解像できていない領域で、被検眼の眼底部の構造物を表すような撮影画像が報告されている(非特許文献2参照)。このように、近年では、共焦点画像及び非共焦点画像の同時撮影を行う装置が非常に注目されてきている。共焦点及び非共焦点に係る眼底画像をともに得る装置は、共振スキャナやガルバノスキャナを用いて測定光を2次元走査している。そのスキャナの制御回路は、光検出器による検出信号を画像表示するためにスキャナの走査に同期して画像生成用の同期信号を発生させるようにしている。このスキャナの走査位置を示す同期信号と光検出による電気的サンプリング位置を合わせることで、眼底画像を形成している。【0006】しかしながら、スキャナの駆動信号と画像生成用の同期信号を同じタイミングで出力すると、実際には共振スキャナの周波数や電気回路系の時間遅れ等があるので、同期信号から得られる光学的走査位置と電気的サンプリング位置とが必ずしも一致しないことがある。補償光学SLO装置は高分解能であるため、そのスキャナ位置の誤差の影響を強く受け、同期信号に基づいてサンプリングされた画像間に位置ずれが目立ってしまい、画像の歪み等の画質劣化が発生する可能性がある。【0007】このような課題を解決するために、特許文献1では、専用のハードウェア構成を用いてスキャナがスキャンしている正確な位置を取得し、または制御する手法を開示している。【0010】AO−SLO画像を、高精細、かつ、広画角、高フレームレートにすればするほど、光学的走査位置と電気的サンプリング位置とのずれが非常に小さくても、画像間の位置ずれが目立ってしまう。このような場合、上述した特許文献1のように、スキャナ位置を検出するハードウェア構成を用いても、非常に小さい光学的走査位置と電気的サンプリング位置とのずれを正確かつ高い応答性で検出することは非常に困難である。【0011】また、上述した光学的走査位置と電気的サンプリング位置とのずれは、装置本体内の温度や電源の不安定性等の影響も受けるため、撮影ごとにずれ方が変化する可能性がある。そのため、上述した特許文献2のように、撮影前に画像補正用のチャート画像を撮影することで、画像歪みを予め取得して画像処理によって歪み補正をしたとしても、再度撮影画像内に歪みが生じてしまい、画質劣化が発生してしまう。眼科装置の場合、同一の共振スキャナからの測定光によって画像補正用のチャート画像と眼底画像とを同時に撮影することは装置の構成上難しいため、上述した特許文献2の手法では画像内の歪みをすべて補正することは困難である。従来の技術では、特別なハードウェア構成や画像補正用のチャート画像等を用いずに、スキャナの特性による画像歪みが補正された画像を取得することが困難であった。本発明は、このような問題点に鑑みてなされたものであり、特別なハードウェア構成や画像補正用のチャート画像等を用いずに、スキャナの特性による画像歪みが補正された画像を取得することが可能な仕組みを提供することを目的とする。より補正されたサンプリング基準位置に基づいて、前記各撮影手段ごとに、前記往復走査のサンプリングデータ配列に基づく画像データを構成する画像構成手段とを有する。した画像生成装置の各手段としてコンピュータを機能させるためのプログラム、並びに、当該プログラムを記憶したコンピュータ読み取り可能な記憶媒体を含む。本発明によれば、特別なハードウェア構成や画像補正用のチャート画像等を用いずに、スキャナの特性による画像歪みが補正された画像を取得することが可能となる。【図1】本発明の実施形態に係る画像生成装置を含む画像生成システムの概略構成の一例を示すブロック図である。【図2】図1に示す補償光学SLO装置の概略構成の一例を示すブロック図である。【図3】本発明の実施形態に係る画像生成装置による画像生成方法における処理手順の一例を示すフローチャートである。【図4】図3のステップS330におけるスキャナ位置候補算出処理の詳細な処理手順の一例を示すフローチャートである。【図5】本発明の実施形態を示し、ガルバノスキャナのトリガ信号のタイムチャート、共振スキャナのトリガ信号のタイムチャート、及び、反射信号のタイムチャートを示す図である。【図6】図4のステップS450における逆FFT処理で得られた逆FFTデータ配列の一例を示すグラフ図である。【図7】図3のステップS350における画像構成処理の際に行われるsine補正を説明するための図である。以下に、図面を参照しながら、本発明を実施するための形態(実施形態)について説明する。本発明の実施形態では、補償光学SLO装置により被検眼の所定領域として網膜(眼底部)を撮影した画像を取得する際に、共振スキャナとガルバノスキャナとの組み合わせによって取得したシグナル値に基づき画像を構成する場合について説明する。具体的に、本実施形態では、複数の撮影手段を用いて共焦点画像及び非共焦点画像の双方を同時に撮影し、撮影により得られたサンプリングデータ配列を用いて画像処理により画像歪みを補正する。共焦点画像及び非共焦点画像は、同じスキャナ走査位置ではサンプリング位置ずれ量も同じはずであるが、それぞれの画像の特徴が異なるため、異なる補正結果になる場合がある。本実施形態では、補正による確からしさを比較することにより、歪み補正の誤検出を低減することが可能である。本実施形態では、補正したサンプリング基準位置を用いて画像データを生成することにより、特別なハードウェア構成や画像補正用のチャート画像等を用いずに、スキャナ(本例では共振スキャナ)の特性による画像歪みが補正された網膜の画像(眼底画像)を取得することが可能となる。【0020】(画像生成装置の概略構成)図1は、本発明の実施形態に係る画像生成装置100を含む画像生成システム10の概略構成の一例を示すブロック図である。画像生成システム10は、図1に示すように、画像生成装置100、補償光学SLO装置200、及び、外部装置300を有して構成されている。画像生成装置100は、スキャナにより被検眼の所定領域を往復走査した測定光による戻り光を受光して電気信号に変換する光電変換手段を含む撮影手段が複数設けられた補償光学SLO装置200(眼科装置)と通信可能に接続された装置である。この画像生成装置100は、図1に示すように、信号情報取得部110、制御部120、記憶部130、画像生成部140、及び、出力部150を有して構成されている。信号情報取得部110は、補償光学SLO装置200から、共振スキャナ及びガルバノスキャナのそれぞれのトリガ信号と、被検眼の網膜からの反射信号等の信号情報を取得する。ここで、トリガ信号とは、共振スキャナ及びガルバノスキャナが特定の位置(基準位置)にあることを示す基準信号である。信号情報取得部110で取得された信号情報は、制御部120を通じて記憶部130に記憶される。制御部120は、画像生成装置100の動作を統括的に制御する。記憶部130は、信号情報取得部110で取得された信号情報や、画像生成装置100の処理に必要な各種の情報、及び、画像生成装置100の処理で得られた各種の情報(各種の画像データを含む)等を記憶する。画像生成部140は、各種の画像データの生成を行うものであり、図1に示すように、パラメータ演算処理部141、及び、画像構成部142を含み構成されている。【0026】パラメータ演算処理部141は、信号情報取得部110で取得した信号情報に基づいてガルバノスキャナ及び共振スキャナのトリガ信号(基準信号)を抽出する。具体的に、本実施形態においては、上述のスキャナが1往復走査する間に出力された当該スキャナの基準信号を抽出する。この基準信号の抽出を行うパラメータ演算処理部141は、抽出手段を構成する。パラメータ演算処理部141は、抽出したトリガ信号(基準信号)をサンプリング基準位置として、補償光学SLO装置200内の複数の撮影手段における各撮影手段ごとに、反射信号(補償光学SLO装置200の光電変換手段で得られた電気信号)に基づき往復走査におけるサンプリング配列を生成する。このサンプリングデータ配列を生成するパラメータ演算処理部141は、データ配列生成手段を構成する。パラメータ演算処理部141は、各撮影手段ごとに、生成した往復走査のサンプリングデータ配列のうちの往路走査によるサンプリングデータ配列と復路走査によるサンプリングデータ配列とを比較してこれらの相関を評価する。この相関を評価するパラメータ演算処理部141は、第1の評価手段を構成する。パラメータ演算処理部141は、上述した相関の評価の結果得られた各撮影手段ごとの第1の評価結果を比較して確からしさを評価する。この確からしさを評価するパラメータ演算処理部141は、第2の評価手段を構成する。パラメータ演算処理部141は、上述した第1の評価結果及び上述した確からしさの評価の結果得られた第2の評価結果に基づいて、上述したサンプリング基準位置を補正する。画像構成部142は、パラメータ演算処理部141により補正されたサンプリング基準位置に基づいて、各撮影手段ごとに、上述した往復走査のサンプリングデータ配列に基づく被検眼の網膜における画像データを構成する。出力部150は、画像生成部140で生成された画像データや、記憶部130に記憶されている各種の情報等を、外部装置300に出力する。補償光学SLO装置200は、上述のスキャナにより被検眼の所定領域を往復走査した測定光による戻り光を受光して電気信号に変換する光電変換手段を含む撮影手段が複数設けられた眼科装置である。補償光学SLO装置200は、被検眼において発生する測定光またはその戻り光の収差を補正するための補償光学系を更に含み構成されている。外部装置300は、画像生成装置100と通信可能に接続された装置であり、例えば、表示装置等の出力装置である。【0031】(補償光学SLO装置の概略構成)次に、図1に示す補償光学SLO装置200の概略構成について説明する。図2は、図1に示す補償光学SLO装置200の概略構成の一例を示すブロック図である。この図2には、図1に示す補償光学SLO装置200の概略構成に加えて、図1に示す画像生成装置100及び外部装置300も図示している。Diode)201、ビームスプリッタ202及び203、補償光学系204、X−Y走査ミラー205、シャックハルトマン波面センサ206、フォーカスレンズ209及び211、Annularミラー210、並びに、光センサ212及び213を有して構成されている。光源であるSLD201から照射された光は被検眼Eの眼底部Erで反射され、その一部がビームスプリッタ203を経由してシャックハルトマン波面センサ206へ導かれ、それ以外はビームスプリッタ202を経由して光センサ212,213へ導かれる。【0034】シャックハルトマン波面センサ206は、被検眼Eの収差を測定するためのデバイスであり、レンズアレイ207にCCD208が接続されて構成されている。このシャックハルトマン波面センサ206において、入射光がレンズアレイ207を透過すると、CCD208に輝点群が投影され、当該投影された輝点の位置ずれに基づき波面収差が測定される。補償光学系204は、シャックハルトマン波面センサ206で測定された波面収差に基づき、収差補正デバイス(可変形状ミラーもしくは空間光位相変調器)を駆動して収差を補正する。当該収差補正された光は、ビームスプリッタ203及び202、フォーカスレンズ209及びAnnularミラー210を経由し、共焦点側の光は光電変換手段である光センサ213に導かれ、非共焦点側の光はフォーカスレンズ211を経由して光電変換手段である光センサ212に導かれる。各光センサ212及び光センサ213に導かれた光は、それぞれ、各光センサ212及び光センサ213において電気信号に変換される。【0037】ここで、Annularミラー210には、共焦点の直系の光を反射し、非共焦点の光を透過するというようなマスクが施されている。ここで、例えば、非特許文献1に記載のように、さらに非共焦点の光を任意の方向に分割し、追加して光センサ(212)を設けることで、分割された方向ごとの非共焦点の光を受光することも可能である。本実施形態では、説明の簡略化のために、非共焦点の光を受光する1つの光センサ212を示しているが、本発明においてはこれに限定されるものではなく、複数の光センサ212を設ける形態も適用可能である。(フレーム数/フレームレート)の情報を取得することができる。向には共振スキャナが設けられ、副走査であるY方向にはガルバノスキャナが設けられている。X−Y走査ミラー205に係る信号情報は、画像生成装置100へ送信され、画像生成装置100において画像データ(動画像データもしくは静止画像データ)を生成する際に用いられる。本実施形態では、フレームレートの高速化のために、主走査である共振スキャナの往路走査及び復路走査の双方で信号情報を取得するものとするが、これに限定されない。被検眼Eの眼底部Er上の特定の深さ位置にフォーカスを合わせるためには、補償光学系204内の収差補正デバイスを用いた調整か、或いは、光学系内に不図示のフォーカス調整用レンズを設置し、当該フォーカス調整用レンズを移動することによる調整か、の少なくともいずれかを用いることができる。【0040】また、被検眼Eの固視微動と呼ばれる不随意的な眼球運動や固視不良による眼球運動或いは顔の動きに伴う被検眼Eの動きの影響を軽減するために、被検眼Eの動きを追尾(トラッキング)する場合がある。この際、眼底部Erの動きの計測は、例えば、画像生成装置100において、生成した眼底画像を用いて、特徴を持った小領域の画像であるテンプレート画像を抽出し、テンプレート画像と最も似た部位を探すパターンマッチングによって行う方法等がある。この場合、例えば、画像生成装置100の制御部120は、被検眼Eの動きを追尾して撮像位置を補正する制御を行う。この撮像位置を補正する制御を行う制御部120は、補正制御手段を構成する。補償光学SLO装置200は、画像生成装置100の制御部120による当該補正制御に基づいて、被検眼Eの動きを追尾した撮像位置の補正を行う。具体的に、例えば、補償光学SLO装置200の光学系内に不図示のスキャナを設置して、照射光の位置を当該スキャナによって計測した眼底部Erの動きに追従させる。【0041】(画像生成装置の処理手順)次に、本発明の実施形態に係る画像生成装置100による画像生成方法における処理手順について説明する。図3は、本発明の実施形態に係る画像生成装置100による画像生成方法における処理手順の一例を示すフローチャートである。【0042】<ステップS310>まず、ステップS310において、信号情報取得部110は、補償光学SLO装置200で取得された少なくとも画像1枚分の信号情報を取得する。ここで取得される信号情報は、被検眼Eの網膜の撮影に用いられるガルバノスキャナのトリガ信号、共振スキャナのトリガ信号、並びに、撮影によって取得された被検眼Eの網膜からの共焦点及び非共焦点それぞれの反射信号の4種類である。ステップS310で取得された信号情報は、制御部120を通じて記憶部130に記憶されて保存される。この際、取得された信号情報に付随するハードウェアの制御情報を取得し、制御部120を通じて記憶部130に記憶する。ここで、制御情報とは、被検眼Eの網膜の反射信号を取得する際にサンプリング周波数やガルバノスキャナの周波数に相当するフレームレート等の情報である。これらの制御情報は、信号情報に付加している撮影情報ファイルに記載されている場合もあれば、信号情報のタグ情報として含まれている場合もある。【0043】<ステップS320>続いて、ステップS320において、パラメータ演算処理部141は、記憶部130に保存されている信号情報から、ガルバノスキャナ及び共振スキャナのトリガ信号(基準信号)を抽出して、トリガ位置を取得する。ここで、トリガ位置とは、反射信号をサンプリングしてピクセル値を得るための基準となる位置(サンプリング基準位置)である。パラメータ演算処理部141は、取得したトリガ位置情報を、制御部120を通じて記憶部130に保存する。図5は、本発明の実施形態を示し、ガルバノスキャナのトリガ信号のタイムチャート、共振スキャナのトリガ信号のタイムチャート、及び、反射信号のタイムチャートを示す図である。具体的に、図5(a)にガルバノスキャナのトリガ信号のタイムチャートを示し、図5(b)に共振スキャナのトリガ信号のタイムチャートを示し、図5(c)に反射信号のタイムチャートを示す。【0045】ガルバノスキャナのトリガ信号及び共振スキャナのトリガ信号は、それぞれ、図5(a)及び図5(b)に示すような信号であるため、パラメータ演算処理部141は、トリガ信号の抽出を閾値処理によって行う。具体的には、パラメータ演算処理部141は、信号強度の変化が或る閾値以上となる場合に、トリガ信号として検出し抽出する。【0046】ガルバノスキャナは、画像の水平方向にレーザビームを走査していて、1つのトリガ信号が検出される間に、水平方向の走査が1回行われ、この間に画像1枚が撮影される。共振スキャナは、1つのトリガ信号が検出される間に、画像の垂直方向2ライン分の走査(往路走査と復路走査)が行われ、この間に画像1枚分の中の画像縦方向の往復2列分が撮影される。そこで、パラメータ演算処理部141は、トリガ信号検出開始のタイミングを、トリガ位置、つまりサンプリング基準位置として抽出する。【0047】また図5(c)に示す反射信号は、被検眼Eの眼底部Erからの反射光の強度を、図2に示す光センサ212または光センサ213で検出した値である。図5(b)に示すように、共振スキャナのトリガ信号が検出されている間のみサンプリングして、この間の信号情報を画像生成に用いる。本実施形態においては、共振スキャナが走査されているときに連続的にサンプリングされていてもよく、また、これに限定されるものではない。トリガ位置の取得方法としては様々な方法が考えられ、ここで用いた方法に限定されるものではない。)における各撮影手段ごとに、反射信号に基づき往復走査におけるサンプリングデータ配列を生成する。パラメータ演算処理部141は、各撮影手段ごとに、生成したサンプリングデータ配列を用いてトリガ信号発信時のスキャナ位置(サンプリング基準位置)の候補を算出する。パラメータ演算処理部141は、算出したスキャナ位置候補情報を、制御部120を通じて記憶部130に保存する。図4は、図3のステップS330におけるスキャナ位置候補算出処理の詳細な処理手順の一例を示すフローチャートである。それぞれにおいて連続的に撮影された複数枚画像の1枚1枚における共振スキャナの1往復走査ごとに行われ、それを1枚の画像の中の往復走査分を繰り返し、複数枚画像に対しても同様に行う。本実施形態におけるスキャナ位置候補算出処理では、各撮影手段ごとに、往復走査のサンプリングデータ配列のうちの往路走査によるサンプリングデータ配列と復路走査によるサンプリングデータ配列とのそれぞれの空間周波数領域を用いて、往路走査によるサンプリングデータ配列と復路走査によるサンプリングデータ配列とを比較してこれらの相関を評価する。Correlation法という画素の空間周波数領域を用いたパターンマッチングにより、往路走査によるサンプリングデータ配列と復路走査によるサンプリングデータ配列とを比較してこれらの相関を評価する。具体的な評価としては、往路走査によるサンプリングデータ配列と復路走査によるサンプリングデータ配列との間の相関が最も高い位置を計算する処理を行う。Correlation法を用いて評価を行う態様も適用できる。以下の図4に示す処理ステップでは、共焦点及び非共焦点の双方ともに同様の計算を行う。≪ステップS410≫まず、ステップS410において、パラメータ演算処理部141は、ステップS310で取得された制御情報に基づいて、トリガ信号の発信時の共振スキャナの位置を、サンプリング基準開始位置の初期位置x0として設定する。ここで、初期位置x0の設定方法にはいくつかの方法があり、例えば、本実施形態においては、トリガ位置よりトリガーディレイ分を考慮した位置といった形で設定される。この他にも、連続して数枚の画像を撮影する場合、その撮影画像より前に撮影した画像の各往復走査ごとの平均的なスキャナ推定位置を初期位置として設定することで、推定位置の探索範囲を狭めるといった方法等も用いることができる。続いて、ステップS420において、パラメータ演算処理部141は、ステップS410で設定したサンプリング基準開始位置の初期位置x0に基づいて、共振スキャナの1往復走査分の1次元サンプリングデータ配列2ライン分を、共焦点及び非共焦点の双方ともに生成する。続いて、ステップS430において、パラメータ演算処理部141は、ステップS420で生成した各往復2ラインごとのサンプリングデータ配列に対して、それぞれ、1次元フーリエ変換を行うことで、1次元の空間周波数領域データ配列を、共焦点及び非共焦点の双方ともに算出して生成する。ここで、空間周波数領域データ配列の算出方法としては、例えば、離散フーリエ変換を用いて算出する方法が適用できるが、より好適には、演算高速化のため高速フーリエ変換(以下、「FFT」と称する)を用いて算出する方法がある。以下の説明においては、空間周波数領域データ配列の算出方法として、FFTを用いて算出する方法を適用し、一般的な窓関数であるハン窓を使用する。【0055】また、FFTを行う領域としては、何らかの構造物等の輝度勾配のある領域が含まれていることが好適であるため、できる限り広い領域で、かつ、スキャン中心の共振スキャナの速度が速い領域を含めるものとする。スキャン中心の共振スキャナの速度が速い領域を含めることは、共振スキャナの走査による画像歪みの影響が他の範囲よりも小さい範囲の往復走査のサンプリングデータ配列を用いて、上述した評価を行うことに相当する。更なる計算の高速化のために、サンプリングデータ配列から輝度勾配の大きな領域を含むより狭い領域のみを用いることや、固定小数点によるFFTやビットシフト等により演算量の削減を行うことも可能である。このように、空間周波数領域の算出に関しては上述した方法に限定されるものではない。本実施形態においては、隣接する、往路走査によるサンプリングデータ配列と復路走査によるサンプリングデータ配列のうち、少なくとも輝度勾配の変化の大きい範囲について、上述した評価を行う態様も適用可能である。続いて、ステップS440において、パラメータ演算処理部141は、ステップS430で生成された往路走査によるサンプリングデータ配列の空間周波数領域データ配列Aと復路走査によるサンプリングデータ配列の空間周波数領域データ配列Bとを、共焦点及び非共焦点の双方ともに合成する。Correlation法での演算に用いる一般的な方法を用いて、Aに対してBの共役を掛け算することで得られる。続いて、ステップS450において、パラメータ演算処理部141は、ステップS440で得られた往復走査のサンプリングデータ配列の空間周波数領域データ配列の合成画像に対して、共焦点及び非共焦点の双方ともに逆フーリエ変換を行う。ここでも、逆フーリエ変換は、逆FFTを用いる。続いて、ステップS460において、パラメータ演算処理部141は、ステップS450で得られた共焦点及び非共焦点の双方の逆FFTデータ配列のうち、最大値を持つデータポイントの位置をそれぞれ探索する。図6は、図4のステップS450における逆FFT処理で得られた逆FFTデータ配列の一例を示すグラフ図である。Peak(Correlationピーク)の大きさ、横軸は往復走査のサンプリング配列の走査方向のOffset値である。図6のグラフを見ると、Offset値が0の時が最大値を示している。ここで、最大値を持つデータポイントの探索方法としては、共振スキャナの周波数や電気回路の遅れ等の物理的に起こりうるズレの範囲内のみを探索する等してもよく、その方法は限定されない。Correlation法を用いるため、この場合、例えば、Correlationピークの最大値の探索を行う態様を採ることができる。Correlation法を用いる形態の場合も、例えば、Correlationピークの最大値の探索を行う態様を採ることができる。続いて、ステップS470において、パラメータ演算処理部141は、ステップS460で得られた共焦点及び非共焦点の双方のデータポイントの位置近傍のデータを用いてそれぞれフィッティングを行い、逆FFTデータにおける最大ピーク値とその時のOffset値をそれぞれ計算する。Correlation法を用いるため、この場合、Correlationピークの最大値近傍での任意のフィッティングにより、Correlationピークの位置を決定する態様を採ることができる。Correlation法を用いる形態の場合も、例えば、Correlationピークの最大値近傍での任意のフィッティングにより、Correlationピークの位置を決定する態様を採ることができる。【0064】このようにフィッティングを行うことで、サブサンプリング単位で基準位置を求めることができる。このフィッティングには、任意の最小二乗法等の方法を用いることができるが、本実施形態では計算が容易なパラボラフィッティングを用いる。図6において、例えば、Offset値が−1,0,1の3点のデータを用いてパラボラフィッティングを行うことで、往復走査のサンプリングデータ配列同士の最も相関の高くなる時のOffset値がサブサンプリング単位で算出できる。このように計算したOffset値に対して共振スキャナの軌道を考慮することで、サンプリング基準開始位置からのシフト量Δxを算出し、x0+Δxをトリガ信号の発信時の実際の共振スキャナ位置である補正されたサンプリング基準位置とみなすことができる。この補正されたサンプリング基準位置情報は、パラメータ演算処理部141から、制御部120を通じて記憶部130に保存される。以上の図4のステップS410〜ステップS470の処理を経ることにより、図3のステップS330におけるスキャナ位置候補算出処理が行われる。ここで、再び、図3の説明に戻る。【0066】<ステップS340>図3のステップS330の処理が終了すると、ステップS340に進む。ステップS340に進むと、パラメータ演算処理部141は、ステップS330で算出された各往復走査ごとのスキャナ位置(サンプリング基準位置)の候補の中から、補正するスキャナ位置をそれぞれ1つ決定する。この際、本実施形態では、パラメータ演算処理部141は、ステップS330における評価の結果得られた各撮影手段ごとの第1の評価結果を比較して確からしさを評価する。具体的に、例えば、スキャナ位置(サンプリング基準位置)の候補の中から1つのスキャナ位置を決定する方法としては、単純に、ステップS330で算出した各撮影手段ごとのCorrelationピークの最大値を、共焦点及び非共焦点の双方で比較して、ピーク値の高さを位置ずれ補正時の相関の高さ(確からしさ)とみなして決定してもよい。他の例として、ステップS330で算出した各撮影手段ごとの、Correlationピークの最大値と2番目に大きい値との比を比較して、その大きさから確からしさを決定してもよい。その他の例として、Correlationピークの位置のうち、共振スキャナを含むハードウェア構成により推定されるサンプリング基準位置のずれ量の範囲(例えば最大ずれ量を考慮した範囲)よりも大きいCorrelationピークの位置を除外するエラー処理を行って、確からしさを決定してもよい。この相関の確からしさを比較する方法としては、様々な方法が考えられ、ここで用いた手法に限定されるものではない。本ステップで決定されたスキャナ位置(サンプリング基準位置)に係る情報は、パラメータ演算処理部141から、制御部120を通じて記憶部130に保存される。【0068】<ステップS350>続いて、ステップS350において、画像構成部142は、ステップS340で決定されたサンプリング基準位置に基づいて、各撮影手段ごとに、上述した往復走査のサンプリングデータ配列に基づく被検眼の網膜における画像データを構成する。具体的に、共振スキャナの1往復走査分の画像2ライン分のsine補正を行って、共焦点及び非共焦点の双方について画像データを構成する。これは、共振スキャナが同一速度で走査していないため、縦軸方向の位置に均一にデータを分割すると、共振スキャナの速度によって画像に歪みが生じてしまうので、共振スキャナの走査軌道をsine波に近似して、各画素に対しサンプリング数に重みづけを行うことで画像歪みを解消することを目的としている。画像構成部142は、構成した画像データを、制御部120を通じて記憶部130に記憶して保存する。図7は、図3のステップS350における画像構成処理の際に行われるsine補正を説明するための図である。図7において、横軸は時間tを示し、また、縦軸は、位置xとし、共振スキャナの走査軌道をステップS310で取得された周波数のsine波であると仮定して示し、振幅を1.0としたものである。この図7に示すsine波の縦軸を仮に10等分して、分割された区間をそれぞれ時間0.0からn=0,1,2,...,9と名付ける。縦軸を分割した直線とsine波との交点の時間をそれぞれt(0),t(1),...と名付ける。この縦軸を分割した区間を共振スキャナが走査している間にサンプリングしたデータを、この位置の画素に割り当てる。すると、サンプリングは数十MHzで時間的に等間隔に行われるが、対応する共振スキャナの1区間を走査している、つまりサンプリングしている時間は、画素により異なっていることがわかる。よって、取得される反射信号の輝度値を、対応する位置のピクセル値に対応付け、1つの画素に複数の信号が対応する場合には、画素値を複数の信号の平均値として取得する。本実施形態では、共振スキャナの走査軌道をsine波に近似しているが、近似方法はこれに限定されない。【0070】ここで、画像構成は、連続的に撮影された複数画像の1枚1枚に対して行う。ステップS320で取得されたガルバノスキャナのトリガ信号の1つ1つに相当する画像1枚分に対して上述した処理を行い、それを枚数分繰り返す。して、すべての画像に対して画像構成を行った後、それらを統合した画像群を構成し、これを記憶部130に記憶して保存する。【0071】<ステップS360>続いて、ステップS360において、出力部150は、ステップS350において得られた網膜の画像データを外部装置300に出力する。これにより、外部装置300において、例えば、出力部150から出力された画像データに基づく画像が表示され、例えば検者等に提示される。出力部150は、ステップS310〜ステップS340において記憶部130に記憶されたスキャナの往復走査ごとのサンプリング基準位置等の画像構成に係る情報を外部装置300(例えばデータベース)に出力する。【0072】ステップS360の処理が終了すると、図3に示すフローチャートの処理が終了する。【0073】なお、本実施形態においては、上述した図3のステップS330における評価を、隣接する往路走査と復路走査に対応する1次元サンプリングデータ配列ごとに行う態様や、複数の往復走査に対応する2次元サンプリングデータ配列ごとに行う態様を適用することが可能である。本実施形態によれば、特別なハードウェア構成や画像補正用のチャート画像等を用いずに、スキャナの特性による画像歪みが補正された画像を取得することが可能となる。【0075】[その他の実施形態]上述した本発明の実施形態では、共振スキャナの往復走査分の画像2ライン分ごとに、トリガ信号の発信時の共振スキャナ位置を推定している。これは、補償光学SLO装置200は、高横分解能であるため、そのスキャナ周波数や電気回路系の時間遅れがスキャナ位置に及ぼすわずかな誤差が、1往復走査ごとに不均一な画像歪みとして影響を及ぼすためである。これは、共振スキャナの走査に係るばらつきが、補償光学SLO装置200本体内の温度や電源の不安定性等の影響も受けるため、スキャナ位置の誤差自体のばらつきも大きいことが原因である。補償光学SLO装置200のスキャナ設定や取得画像の撮影条件、もしくは、共振スキャナ以外のスキャナの使用(ポリゴンミラー等)によっては、この画像歪みがある一定の周期ごとに現れる場合もある。この場合には、画像歪みの現れる周期に相当する数の往復走査ごとにスキャナ位置の推定を行うことで推定ステップの回数を減らすことも可能である。他には、スキャナ位置の推定を往復走査1ラインで行うのではなく、複数の往復走査ラインで同時に計算することも可能である。具体的には、複数の往路走査ライン及び復路走査ラインを、それぞれ二次元の画像として扱い、2次元FFTによる空間周波数領域の比較により、スキャナ走査方向の位置ずれを計算するという方法である。この方法を採れば、一括で、均一のサンプリング基準位置を計算することができる。【0077】また、上述した本発明の実施形態では、図3のステップS340においてサンプリング基準位置(スキャナ位置)の候補の中から自動でサンプリング基準位置を決定する例を示したが、本発明はこれに限定されるものではない。それぞれのサンプリング基準位置の候補のサンプリング基準位置を用いた画像データを画像構成部142で構成し、出力部150を通して外部装置300である表示装置に表示する。ステップS340において、パラメータ演算処理部141は、外部装置300である表示装置に表示された画像のうち、ユーザによって選択された画像に基づいて、上述した評価を行ってサンプリング基準位置(スキャナ位置)を決定する形態も適用可能である。この形態の場合には、例えばCorrelationピークの最大値の大きさ等の参考となるパラメータを画像と同時に表示してもよく、確からしさがある閾値以下のスキャナ位置のみをユーザに選択させるようにしてもよい。【0078】さらには、歪みが補正された画像を用いることで、表示する画像の画質向上や画像演算の精度向上が期待されるため、上述した実施形態の画像歪み補正を撮影後すぐに実行するための高速化をハードウェア構成と演算の並列化により実施する方法もある。任意のFPGAやGPU等のマルチコア演算器等を追加して、画像生成部140で行う演算を複数の往復走査について並列化によるパイプライン処理で行うことも可能である。【0079】また、上述の本発明の実施形態では、スキャナで往復走査する被検眼Eの所定領域として、被検眼Eの網膜(眼底部Er)を適用した例について説明を行ったが、本発明においてはこれに限定されるものではなく、被検眼Eの他の部分を適用することも可能である。本発明は、上述の実施形態の1以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける1つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。1以上の機能を実現する回路(例えば、ASIC)によっても実現可能である。このプログラム及び当該プログラムを記憶したコンピュータ読み取り可能な記憶媒体は、本発明に含まれる。【0081】なお、上述した本発明の実施形態は、いずれも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。本発明はその技術思想、又はその主要な特徴から逸脱することなく、様々な形で実施することができる。Reserved.7。