WO2023162586A1

WO2023162586A1 - 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体

Info

Publication number: WO2023162586A1
Application number: PCT/JP2023/002961
Authority: WO
Inventors: 嘉晃安野; 修一巻田; 達夫山口; 神之介東
Original assignee: 国立大学法人筑波大学; 株式会社トプコン
Priority date: 2022-02-24
Filing date: 2023-01-31
Publication date: 2023-08-31
Also published as: JP2023123245A

Abstract

実施形態の眼科検査装置１は、光走査を用いてサンプルのイメージングを行う走査型イメージング装置である。眼科検査装置１は、光スキャナ４４、ＯＣＴユニット１００、及び制御部２１０を用いた光走査によりサンプルのデータを冗長的に収集してデータセットを取得する。画像データ構築部２２０は、取得されたデータセットに基づいて画像セットを生成する。画像選択部２２５は、生成された画像セットから複数の画像を選択する。画像データ処理部２３１は、選択された複数の画像をそれぞれ基準とした複数のレジストレーションをこの画像セットに適用することによって、この画像セットの相対位置調整を実行する。

Description

走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体

　本開示は、走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体に関する。

　イメージング技術の１つに走査型イメージングがある。走査型イメージングとは、サンプルの複数の箇所に順次にビームを照射してデータを収集し、収集されたデータからサンプルの像を構築する技術である。

　走査型イメージングの手法には、例えば、サンプル上における投影像がスポット状のビームを用いるスポット型の走査手法、投影像がライン状のビームを用いるライン型の走査手法、及び、投影像がエリア状のビームを用いるエリア型の走査手法などがある。スポット型の走査手法の例としてフーリエドメイン型の光コヒーレンストモグラフィ（ＯＣＴ；光干渉断層撮影法）があり、ライン型の走査手法の例としてラインスキャンカメラがあり、エリア型の走査手法の例としてフルフィールド型のＯＣＴがある。

　光を利用した走査型イメージングモダリティの例としてＯＣＴが知られている。ＯＣＴは、光散乱媒質をマイクロメートルレベル又はそれ以下の分解能で画像化することが可能な技術であり、医用イメージングや非破壊検査などに応用されている。

　ＯＣＴは、低コヒーレンス干渉法に基づく技術である。医用イメージングでは、深達性や侵襲性を考慮し、近赤外領域のプローブ光（測定光）が利用されている。なお、近赤外光以外の波長帯の光（電磁波）や超音波や放射線を利用した走査型イメージングも知られている。

　ＯＣＴの応用が最も進んでいる分野として眼科分野があり、フーリエドメインＯＣＴ機能を搭載した製品がクリニックにまで普及している。眼科画像診断では、２次元的なイメージングだけでなく、３次元的なイメージング・構造解析・機能解析なども実用化されており、診断の強力なツールとして広く利用されるに至っている。また、眼科分野では、走査型レーザー検眼鏡（ＳＬＯ）、ラインスキャンカメラ、フルフィールドＯＣＴなど、フーリエドメインＯＣＴ以外の走査型イメージングも利用されている。

　ＯＣＴやＳＬＯで利用される走査モードには様々なものがあるが、モーションアーティファクト補正などを目的とした所謂「リサジュー（Lissajous）スキャン」が近年注目を集めている（例えば、特許文献１～４、非特許文献１及び２を参照）。

　リサジュースキャンは、互いに直交する２つの単振動を合成して得られるリサジュー曲線として生成される２次元パターンにしたがって実行されるスキャンである。リサジュースキャンでは、或る程度の大きさを持つ複数のループ（互いに交差する複数のサイクル）を描くように測定光を高速走査するため、１つのサイクルからのデータ取得時間の差を実質的に無視することができる。また、異なるサイクルの交差領域を参照してサイクル間の位置合わせを行えるため、サンプルの動きに起因するモーションアーティファクトを補正することが可能である。このようなリサジュースキャンの特徴に着目し、眼科分野では眼球運動に起因するモーションアーティファクトへの対処を図っている。

　非特許文献１に記載された技術においては、リサジュースキャンで収集されたデータセットを比較的大きな動きが介在しない複数のサブボリュームに分割し、各サブボリュームの正面プロジェクション画像（en face projection）を構築する。このような正面プロジェクション画像はストリップと呼ばれる。これらストリップの間のレジストレーションを行うことでモーションアーティファクトが補正された画像が得られる。

特開２０１６－１７９１５号公報特開２０１８－６８５７８号公報特開２０１８－１４０００４号公報特開２０１８－１４００４９号公報

Yiwei Chen, Young-Joo Hong, Shuichi Makita, and Yoshiaki Yasuno, "Three-dimensional eye motion correction by Lissajous scan optical coherence tomography", Biomedical Optics EXPRESS, Vol. 8, No. 3, 1 Mar 2017, PP. 1783-1802 Yiwei Chen, Young-Joo Hong, Shuichi Makita, and Yoshiaki Yasuno, "Eye-motion-corrected optical coherence tomography angiography using Lissajous scanning", Biomedical Optics EXPRESS, Vol. 9, No. 3, 1 Mar 2018, PP. 1111-1129

　この発明の１つの目的は、モーションアーティファクト補正の誤差の低減を図ることにある。

　例示的な実施形態の走査型イメージング装置は、光走査を用いてサンプルのイメージングを行う走査型イメージング装置であって、前記光走査により前記サンプルのデータを冗長的に収集してデータセットを取得するデータセット取得部と、前記データセットに基づいて画像セットを生成する画像セット生成部と、前記画像セットから複数の画像を選択する画像選択部と、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用することによって前記画像セットの相対位置調整を実行する画像位置調整部とを含んでいる。

　前記画像位置調整部は、前記複数のレジストレーションを前記画像セットに適用することによって前記画像セットの複数の相対位置情報を生成し、前記複数の相対位置情報に基づいて前記相対位置調整を実行するように構成されていてよい。前記画像位置調整部は、前記複数のレジストレーションのそれぞれにおいて、前記複数の画像の１つを初期基準画像として前記画像セットに再帰的なレジストレーションを適用することによって前記画像セットの相対位置情報を生成するように構成されていてよい。前記画像位置調整部は、前記複数の相対位置情報に基づいて前記相対位置調整のための調整量を決定するように構成されていてよい。前記画像選択部は、前記光走査の時間軸に基づいて前記画像セットから前記複数の画像を選択するように構成されていてよい。前記画像選択部は、前記画像セットを前記時間軸の複数の区間に対応する複数のサブセットに分割し、前記複数のサブセットのそれぞれから画像を選択することによって、前記画像セットから前記複数の画像を選択するように構成されていてよい。前記画像選択部は、前記光走査の適用エリアに基づいて前記画像セットから前記複数の画像を選択するように構成されていてよい。前記画像選択部は、前記画像セットを前記適用エリアの複数のサブエリアに対応する複数のサブセットに分割し、前記複数のサブセットのそれぞれから画像を選択することによって、前記画像セットから前記複数の画像を選択するように構成されていてよい。前記画像選択部は、画像サイズに基づいて前記複数のサブセットのそれぞれから画像を選択するように構成されていてよい。前記画像選択部は、前記複数のサブセットのそれぞれから、当該サブセットにおける最大サイズの画像を選択するように構成されていてよい。前記画像選択部は、前記複数のサブセットのそれぞれに含まれる画像に対して画像サイズに基づく順序付けを行い、前記複数のサブセットのそれぞれから、当該サブセットにおける画像に対して前記順序付けにより付された順序に基づいて画像を選択するように構成されていてよい。前記画像選択部は、前記複数のサブセットのうちの第１サブセットから選択された第１画像を基準としたレジストレーションの成否を判定し、失敗と判定された場合、前記第１サブセットにおける画像に対して前記順序付けにより付された順序に基づいて前記第１サブセットから第２画像を選択するように構成されていてよい。前記画像選択部は、前記画像セットに含まれる画像に対して画像サイズに基づく順序付けを行い、前記順序付けにより画像に付された順序に基づいて前記画像セットから前記複数の画像を選択するように構成されていてよい。前記画像選択部は、前記画像セットから選択された前記複数の画像のうちの一の画像を基準としたレジストレーションの成否を判定し、失敗と判定された場合、前記画像セットにおける画像に対して前記順序付けにより付された順序に基づいて前記画像セットから前記複数の画像のいずれとも異なる画像を新たに選択するように構成されていてよい。実施形態の走査型イメージング装置は、前記光走査が前記サンプルに適用されているときの前記サンプルの運動を表す運動情報を生成する運動情報生成部を更に含んでいてよく、前記画像選択部は、前記運動情報に基づいて前記画像セットから前記複数の画像を選択するように構成されていてよい。前記画像選択部は、前記運動情報から所定の運動パラメータの値を算出し、前記値を評価し、前記評価の結果に基づいて前記画像セットから前記複数の画像を選択するように構成されていてよい。前記運動パラメータは、運動の大きさ及び頻度の少なくとも一方を含んでいてよく、前記画像選択部は、前記運動パラメータの前記値を所定の閾値と比較し、前記閾値以上の前記値に対応する期間に収集されたデータに基づく画像を前記画像セットから除外したサブセットから前記複数の画像を選択するように構成されていてよい。前記データセット取得部は、前記光走査のための光を互いに異なる第１方向及び第２方向に偏向可能な偏向器を含んでいてよく、前記第１方向における偏向方向の変化を第１周期で繰り返しつつ、前記第２方向における偏向方向の変化を前記第１周期と異なる第２周期で繰り返すように構成されていてよい。

　例示的な実施形態の情報処理装置は、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットを受け付けるデータセット受付部と、前記データセットに基づいて画像セットを生成する画像セット生成部と、前記画像セットから複数の画像を選択する画像選択部と、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用することによって前記画像セットの相対位置調整を実行する画像位置調整部とを含んでいる。

　例示的な実施形態の情報処理装置は、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットに基づいて生成された画像セットを受け付ける画像セット受付部と、前記画像セットから複数の画像を選択する画像選択部と、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用することによって前記画像セットの相対位置調整を実行する画像位置調整部とを含んでいる。

　例示的な実施形態の走査型イメージング装置の制御方法は、プロセッサと光走査を行うスキャナとを含む走査型イメージング装置を制御する方法であって、前記スキャナに、前記サンプルのデータを冗長的に収集させてデータセットを取得させ、前記プロセッサに、前記データセットに基づいて画像セットを生成させ、前記プロセッサに、前記画像セットから複数の画像を選択させ、前記プロセッサに、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用させることによって前記画像セットの相対位置調整を実行させる。

　例示的な実施形態の情報処理装置の制御方法は、プロセッサを含む情報処理装置を制御する方法であって、前記プロセッサに、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットを受け付けさせ、前記プロセッサに、前記データセットに基づいて画像セットを生成させ、前記プロセッサに、前記画像セットから複数の画像を選択させ、前記プロセッサに、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用させることによって前記画像セットの相対位置調整を実行させる。

　例示的な実施形態の情報処理装置の制御方法は、プロセッサを含む情報処理装置を制御する方法であって、前記プロセッサに、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットに基づいて生成された画像セットを受け付けさせ、前記プロセッサに、前記画像セットから複数の画像を選択させ、前記プロセッサに、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用させることによって前記画像セットの相対位置調整を実行させる。

　例示的な実施形態の情報処理方法は、光走査によりサンプルのデータを冗長的に収集して取得されたデータセットを処理する情報処理方法であって、前記データセットに基づいて画像セットを生成し、前記画像セットから複数の画像を選択し、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用することによって前記画像セットの相対位置調整を実行する。

　例示的な実施形態のプログラムは、本実施形態のいずれかの方法をコンピュータに実行させるプログラムである。

　例示的な実施形態の記録媒体は、本実施形態のいずれかのプログラムが記録されたコンピュータ可読な非一時的記録媒体である。

　例示的な実施形態によれば、モーションアーティファクト補正の誤差を低減することができる。

実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の構成の一例を表す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の構成の一例を表す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の構成の一例を表す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の構成の一例を表す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の構成の一例を表す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の構成の一例を表す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行するリサジュースキャンのパターンの例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する初期基準画像選択処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する初期基準画像選択処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する初期基準画像選択処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する初期基準画像選択処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する調整量算出処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する調整量算出処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する調整量算出処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する調整量算出処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）が実行する調整量算出処理の例を示す概略図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の動作の一例を表すフローチャートである。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の動作の一例を表すフローチャートである。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の動作の一例を表すフローチャートである。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の動作の一例を表すフローチャートである。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の動作の一例を説明するための図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の動作の一例を説明するための図である。実施形態の例示的な態様に係る走査型イメージング装置（眼科検査装置）の動作の一例を説明するための図である。実施形態の例示的な態様に係る情報処理装置の構成の一例を表す概略図である。実施形態の例示的な態様に係る情報処理装置の構成の一例を表す概略図である。

　実施形態の幾つかの例示的な態様について図面を参照しながら説明する。本開示において説明される例示的な態様は、走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体に関するものであるが、実施形態はこれらの態様に限定されない。

　また、本開示において説明される幾つかの例示的な態様は、主に眼科分野への応用を提供するものであるが、実施形態はこれらの態様に限定されない。例えば、幾つかの例示的な実施形態は、眼科以外の診療科への応用や、診断方法などの医学的方法への応用や、医療分野以外の分野（生物学、非破壊検査など）への応用など、走査型イメージングが利用されている任意の分野への応用又は走査型イメージングが利用される可能性のある任意の分野への応用を提供することができる。

　本開示において引用された文献に開示された事項や、他の任意の公知技術に関する事項を、実施形態に組み合わせることが可能である。本開示では、特に言及しない限り、被検者（被検眼）の「部位」と、その部位に対応する「画像データ」と、その画像データの視覚的表現である「画像」とを区別しないものとする。

　以下に説明する例示的な態様に係る走査型イメージング装置は、眼科分野の検査に使用される装置（眼科検査装置と呼ぶ）であり、フーリエドメインＯＣＴ（特にスウェプトソースＯＣＴ）を利用して生体眼の眼底を計測及び画像化することが可能に構成されている。実施形態に採用可能なＯＣＴのタイプは、スウェプトソースＯＣＴに限定されず、例えばスペクトラルドメインＯＣＴ又はタイムドメインＯＣＴであってもよい。また、以下に説明する例示的な態様に係る走査型イメージング装置は、スポット型の走査手法を用いているが、他のタイプの走査手法（例えば、ライン型の走査手法、エリア型の走査手法）を用いる走査型イメージング装置に対して本開示に係る任意の事項を適用できることは、当業者であれば理解できるであろう。

　幾つかの例示的な態様に係る走査型イメージング装置は、ＯＣＴ以外の走査型イメージングモダリティを利用可能であってよい。幾つかの例示的な態様は、ＳＬＯなどの任意の光走査型イメージングモダリティを採用することができる。また、幾つかの例示的な態様に適用可能な走査型イメージングモダリティは、光走査型イメージングモダリティでなくてもよく、例えば、光以外の電磁波を利用した走査型イメージングモダリティ、超音波を利用した走査型イメージングモダリティなどであってもよい。

　幾つかの例示的な態様に係る眼科検査装置は、ＯＣＴスキャンで収集されたデータの処理及び／又はＳＬＯで収集されたデータの処理に加え、他のイメージングモダリティで取得されたデータを処理可能であってよい。この「他のイメージングモダリティ」は、例えば、眼底カメラ、スリットランプ顕微鏡、眼科手術用顕微鏡などの任意の眼科イメージングモダリティであってよい。幾つかの例示的な態様に係る眼科検査装置は、ＯＣＴ及びＳＬＯ以外のイメージングモダリティの機能を有していてよい。また、幾つかの例示的な態様に係る眼科検査装置は、特開２０１３－２４８３７６号公報に開示された、２以上の前眼部カメラを用いたステレオアライメント機能を有していてもよい。ステレオアライメントは、２以上の前眼部カメラでそれぞれ得られた２以上の前眼部画像の解析から被検眼の３次元位置を特定して光学系のアライメントを行う機能である。

　走査型イメージングが適用される対象（サンプル）は任意であってよい。以下の例示的な態様で説明される眼科応用では主に眼底に走査型イメージングが適用されるが、走査型イメージングの適用対象は眼底に限定されず、前眼部や硝子体など眼の任意の部位であってもよい。幾つかの例示的な態様に係る構成や機能を、眼以外の生体の部位（組織）の計測やイメージングに応用することが可能であり、また、生体以外の物体（例えば、動きを伴うもの）又はその部分の計測やイメージングに応用することが可能である。すなわち、幾つかの例示的な態様の産業上の利用分野は、眼科関連分野に限定されず、医療関連分野や獣医学関連分野や生物学関連分野を含んでもよく、より一般に、局所的運動及び／又は大域的運動を伴う任意の対象物（サンプル）に関連する分野を含んでいてもよい。幾つかの例示的な態様の産業上の利用分野は、運動を伴わない対象物に関連する分野を含んでいてもよい。例えば、周囲の物体や環境からの影響によって移動したり変形したりする物体を対象物としてもよい。なお、本開示に係る走査型イメージングはモーションアーティファクト補正の向上を１つの目的としたものであるから、本開示における典型的な対象物は、能動的運動を伴う物体又は受動的運動を伴う物体であってよく、或いは、能動的運動を伴う物体及び／又は受動的運動を伴う物体を含む物体であってよい。

　以下に説明する実施形態の様々な態様はいずれも単なる例示であり、発明を限定することを意図したものではない。

　前述したように、本開示に係る走査型イメージングの１つの目的は、モーションアーティファクト補正の向上を図ることにある。従来のリサジュースキャンのモーションアーティファクト補正では、リサジュースキャンで収集されたデータ（ボリューム）を対象物の運動の影響を無視することができるサイズの複数のセグメント（複数のサブボリューム）に分割し、これらセグメントのうちの１つを最初の基準画像（初期基準画像、初期参照画像、初期基準ストリップ、初期基準セグメントなどと呼ぶことがある）に選択し、この単一の初期基準画像を起点（第１回目の処理における基準画像）として再帰的なセグメント間の位置ずれの推定（レジストレーション及びマージ処理）を実行している。この再帰的処理（反復的処理）では、選択された単一の初期基準画像に対して他のセグメントを順次にレジストレーション及びマージ（結合）することにより、逐次に形成される基準画像が表現する範囲を拡大していく。この方法によると、初期基準画像の選択が最終画像の品質に悪影響を与えるケースが想定される。例えば、発明者らの検討によれば、十分に適切な初期選択画像が選択されなかった場合に、レジストレーション誤差の大きいセグメントがそのままマージされ、明らかに不連続な構造として最終画像中に表現されてしまうことがあった。

　本開示に係る技術は、このような問題に対処するために発明者らが見出したものであり、まずはその概要をここで説明する。なお、本開示に係る技術のより詳細且つより正確な説明は、後述する非限定的な例示的態様によって提供される。

　概要的には、本開示に係る技術は、スキャンにより対象物から収集されたデータを複数のセグメントに分割し、これら複数のセグメントから複数の初期基準画像を取得し、各初期基準画像を起点として再帰的なセグメント間の位置ずれ推定を実行するように構成されている。換言すると、本開示に係る技術は、概要的には、収集されたデータを複数のセグメントに分割し、これら複数のセグメントから複数の画像を選択し、これら複数の画像をそれぞれ基準とした複数のレジストレーション（及びマージ処理）を複数のセグメントに適用することで複数のセグメントの相対位置調整を実行するように構成されている。このような本開示に係る技術を用いても依然としてモーションアーティファクトが残存することもあるが、後段で実行可能な詳細なモーションアーティファクト補正によってこの残存誤差を補正することができる。

　幾つかの例示的な態様では、複数のレジストレーションで得られた複数の位置ずれ推定量から尤もらしい位置ずれ量を求めることができる。また、複数の初期基準画像の選択方法は任意であってよいが、その幾つかの例を本開示は提供する。例えば、非限定的な初期基準画像選択方法として、収集されたデータを時間的に区分して得られた複数の区間のそれぞれから初期基準画像を選択する方法や、収集されたデータを空間的に分割して得られた複数の領域のそれぞれから初期基準画像を選択する方法が、本開示において説明される。

　このような本開示に係る技術によれば、従来の技術と比較して、モーションアーティファクト補正の誤差の低減を図ることが可能となり、ひいては、イメージング精度（計測精度）の向上を図ることが可能になる。眼科分野への応用においては、通常の眼科画像診断におけるＯＣＴ画像計測精度の向上や、診断精度の向上を図ることができる。また、本開示に係る技術は、現在の眼科ＯＣＴ装置に対して特殊なハードウェアを追加することなく実現でき、比較的安価で実用化することが可能であるため、眼科検診などにおける使用も想定可能となる。また、従来よりも高い精度でのモーションアーティファクト補正を実現できるため、高精度の画像計測を安定的に行うことが可能になる。これにより、例えば、病変の微小な経時変化の追跡が可能になり、更には、病変の中長期的な微小変化の追跡も可能になる。眼科以外の分野への応用においても同様の作用、効果が奏されることは、当業者であれば容易に理解できるであろう。

　本開示では、走査型イメージング装置の実施形態の幾つかの例示的な態様を説明する。本実施形態の走査型イメージング装置は、光走査を用いて対象物（サンプル）のイメージング（画像計測、光計測）を行う装置である。

　本実施形態の走査型イメージング装置は、光走査によりサンプルのデータを冗長的に収集してデータセットを取得するように構成されている。この冗長的データ収集は、サンプルの特定箇所から２回以上データを収集するデータ収集態様である。本実施形態では、サンプルの複数の箇所のそれぞれから２回以上データ収集を行う。特許文献１～４並びに非特許文献１及び２に記載されているリサジュースキャンは、冗長的データ収集の例を提供している。リサジュースキャンでは、冗長的にデータ収集された箇所は、或るサイクルと別のサイクルとの交差位置に相当する。リサジュースキャンを用いたモーションアーティファクト補正は、或る交差位置から冗長的に取得されたデータ間のレジストレーションを行うことによって、その交差位置で交差するサイクル間のずれ量を求めるものである。

　更に、本実施形態の走査型イメージング装置は、冗長的データ収集により取得されたデータセットに基づいて画像セットを生成するように構成されている。この画像セットは、位置的な冗長性を有する複数の画像を含んでおり、例えば、後述する複数のストリップ又は複数のサブボリュームを含んでいてよい。

　加えて、本実施形態の走査型イメージング装置は、生成された画像セットから複数の画像を選択するように構成されている。ここで選択される複数の画像のそれぞれは、前述した初期基準画像に相当する。また、画像セットから選択される複数の初期基準画像の個数は任意であってよく、例えば、画像セットに含まれる複数の画像の一部のみを初期基準画像に指定してもよいし、画像セットに含まれる複数の画像の全てを初期基準画像に指定してもよい。なお、各初期基準画像は、後段の処理において実行される再帰的なレジストレーション（再帰的に実行されるレジストレーションと再帰的に実行されるマージ処理とを組み合わせた後述の処理）の起点として用いられる。すなわち、各初期基準画像は、再帰的レジストレーションのうちの第１回目のレジストレーションにおいて基準画像として用いられる。

　更に加えて、本実施形態の走査型イメージング装置は、選択された複数の画像（複数の初期基準画像）をそれぞれ基準とした複数のレジストレーションを上記画像セットに適用することによって、この画像セットの相対位置調整を実行するように構成されている。

＜眼科検査装置の構成＞
　図１に示す眼科検査装置１は、実施形態に係る走査型イメージング装置の１つの例示的な態様である。眼科検査装置１は、眼底カメラユニット２、ＯＣＴユニット１００、及び演算制御ユニット２００を含む。眼底カメラユニット２には、被検眼Ｅを正面から撮影するための要素群（光学要素、機構など）が設けられている。ＯＣＴユニット１００には、被検眼ＥにＯＣＴスキャンを適用するための要素群（光学要素、機構など）の一部が設けられている。ＯＣＴスキャンのための要素群の他の部分の少なくとも一部は、眼底カメラユニット２に設けられている。演算制御ユニット２００は、各種の演算や制御を実行する１以上のプロセッサを含む。これらに加え、眼科検査装置１は、被検者の顔を支持するための要素や、ＯＣＴスキャンが適用される部位を切り替えるための要素を備えていてもよい。前者の要素の例として、顎受けや額当てがある。後者の要素の例として、ＯＣＴスキャン適用部位を眼底から前眼部に切り替えるために使用されるレンズユニットがある。

　本開示に記載された幾つかの要素の機能は、回路構成（circuitry）又は処理回路構成（processing circuitry）を用いて実装される。回路構成又は処理回路構成は、開示された機能を実行するように構成及び／又はプログラムされた、汎用プロセッサ、専用プロセッサ、集積回路、ＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）、ＡＳＩＣ（Application Specific Integrated Circuit）、プログラマブル論理デバイス（例えば、ＳＰＬＤ（Simple Programmable Logic Device）、ＣＰＬＤ（Complex Programmable Logic Device）、ＦＰＧＡ（Field Programmable Gate Array）、従来の回路構成、及びそれらの任意の組み合わせのいずれかを含む。プロセッサは、トランジスタ及び／又は他の回路構成を含む、処理回路構成又は回路構成とみなされる。本開示において、回路構成、ユニット、手段、又はこれらに類する用語は、開示された機能を実行するハードウェア、又は、開示された機能を実行するようにプログラムされたハードウェアである。ハードウェアは、本開示に記載されたいずれかのハードウェアであってよく、或いは、本開示に記載された機能を実行するようにプログラム及び／又は構成された既知のハードウェアであってもよい。ハードウェアが或るタイプの回路構成とみなされ得るプロセッサである場合、回路構成、ユニット、手段、又はこれらに類する用語は、ハードウェアとソフトウェアとの組み合わせであり、このソフトウェアはハードウェア及び／又はプロセッサを構成するために使用される。

＜眼底カメラユニット２＞
　眼底カメラユニット２には、被検眼Ｅの眼底Ｅｆを撮影するための光学系が設けられている。取得される眼底Ｅｆの画像（眼底像、眼底写真などと呼ばれる）は、観察画像、撮影画像などの正面画像である。観察画像は、例えば近赤外光を用いた動画撮影により得られ、アライメント、フォーカシング、トラッキングなどに利用される。撮影画像は、例えば可視領域又は赤外領域のフラッシュ光を用いた静止画像である。

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

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

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

　液晶ディスプレイ（ＬＣＤ）３９は固視標（固視標画像）を表示する。ＬＣＤ３９から出力された光束は、その一部がハーフミラー３３Ａに反射され、ミラー３２に反射され、撮影合焦レンズ３１及びダイクロイックミラー５５を経由し、孔開きミラー２１の孔部を通過する。孔開きミラー２１の孔部を通過した光束は、ダイクロイックミラー４６を透過し、対物レンズ２２により屈折されて眼底Ｅｆに投射される。ＬＣＤ３９における固視標画像の表示位置を変更することで、被検眼Ｅの視線を誘導する方向（固視方向、固視位置）が変更される。ＬＣＤなどの表示デバイスの代わりに、例えば、発光素子アレイ、又は、発光素子とこれを移動する機構との組み合わせを用いてもよい。

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

　従来と同様に、本例のアライメント指標像は、アライメント状態により位置が変化する２つの輝点像からなる。被検眼Ｅと光学系との相対位置がｘｙ方向（ｘ座標とｙ座標とで表現される方向）に変化すると、２つの輝点像が一体的にｘｙ方向に変位する。被検眼Ｅと光学系との相対位置がｚ方向（ｚ座標で表現される方向）に変化すると、２つの輝点像の間の相対位置（距離）が変化する。ｚ方向における被検眼Ｅと光学系との間の距離が既定のワーキングディスタンスに一致すると、２つの輝点像が重なり合う。ｘｙ方向において被検眼Ｅの位置と光学系の位置とが一致すると、所定のアライメントターゲット内又はその近傍に２つの輝点像が提示される。ｚ方向における被検眼Ｅと光学系との間の距離がワーキングディスタンスに一致し、且つ、ｘｙ方向において被検眼Ｅの位置と光学系の位置とが一致すると、２つの輝点像が重なり合ってアライメントターゲット内に提示される。

　オートアライメントでは、データ処理部２３０が、２つの輝点像の位置を検出し、主制御部２１１が、２つの輝点像とアライメントターゲットとの位置関係に基づいて後述の移動機構１５０を制御する。マニュアルアライメントでは、主制御部２１１が、被検眼Ｅの観察画像とともに２つの輝点像を表示部２４１に表示させ、ユーザーが、表示された２つの輝点像を参照しながら操作部２４２を用いて移動機構１５０を動作させる。

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

　孔開きミラー２１とダイクロイックミラー５５との間の撮影光路に、視度補正レンズ７０及び７１を選択的に挿入することができる。視度補正レンズ７０は、強度遠視を補正するためのプラスレンズ（凸レンズ）である。視度補正レンズ７１は、強度近視を補正するためのマイナスレンズ（凹レンズ）である。

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

　リトロリフレクタ４１は、図１に示す矢印の方向に移動可能とされ、それにより測定アームの長さが変更される。測定アーム長の変更は、例えば、眼軸長に応じた光路長補正や、干渉状態の調整などに利用される。

　分散補償部材４２は、参照アームに配置された分散補償部材１１３（後述）とともに、測定光ＬＳの分散特性と参照光ＬＲの分散特性とを合わせるよう作用する。

　ＯＣＴ合焦レンズ４３は、測定アームのフォーカス調整を行うために測定アームに沿って移動される。撮影合焦レンズ３１の移動、フォーカス光学系６０の移動、及びＯＣＴ合焦レンズ４３の移動を連係的に制御することができる。

　光スキャナ４４は、実質的に、被検眼Ｅの瞳孔と光学的に共役な位置に配置される。光スキャナ４４は、測定アームにより導かれる測定光ＬＳを偏向する。光スキャナ４４は、例えば、ｘ方向のスキャンを行うためのガルバノミラーと、ｙ方向のスキャンを行うためのガルバノミラーとを含む、２次元スキャンが可能なガルバノスキャナである。

＜ＯＣＴユニット１００＞
　図２に例示するように、ＯＣＴユニット１００には、スウェプトソースＯＣＴを適用するための光学系が設けられている。この光学系は干渉光学系を含む。この干渉光学系は、波長可変光源（波長掃引型光源）からの光を測定光と参照光とに分割し、測定アームによって被検眼Ｅに導かれた測定光の戻り光と参照アームによって導かれた参照光とを重ね合わせて干渉光を生成し、この干渉光を検出する。干渉光学系により得られたデータ（検出信号）は、干渉光のスペクトルを表す信号であり、演算制御ユニット２００に送られる。

　光源ユニット１０１は、例えば、少なくとも近赤外波長帯において出射波長を高速で変化させる近赤外波長可変レーザーを含む。光源ユニット１０１から出力された光Ｌ０は、光ファイバ１０２により偏波コントローラ１０３に導かれてその偏光状態が調整される。更に、光Ｌ０は、光ファイバ１０４によりファイバカプラ１０５に導かれて測定光ＬＳと参照光ＬＲとに分割される。測定光ＬＳの光路は測定アームと呼ばれ、参照光ＬＲの光路は参照アームと呼ばれる。

　参照光ＬＲは、光ファイバ１１０によりコリメータ１１１に導かれて平行光束に変換され、光路長補正部材１１２及び分散補償部材１１３を経由し、リトロリフレクタ１１４に導かれる。光路長補正部材１１２は、参照光ＬＲの光路長と測定光ＬＳの光路長とを合わせるよう作用する。分散補償部材１１３は、測定アームに配置された分散補償部材４２とともに、参照光ＬＲと測定光ＬＳとの間の分散特性を合わせるよう作用する。リトロリフレクタ１１４は、これに入射する参照光ＬＲの光路に沿って移動可能であり、それにより参照アームの長さが変更される。参照アーム長の変更は、例えば、眼軸長に応じた光路長補正や、干渉状態の調整などに利用される。

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

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

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

　検出器１２５は、例えばバランスドフォトダイオードを含む。バランスドフォトダイオードは、一対の干渉光ＬＣをそれぞれ検出する一対のフォトディテクタを有し、これらにより得られた一対の検出結果の差分を出力する。検出器１２５は、この出力（検出信号）をデータ収集システム（ＤＡＳ）１３０に送る。

　データ収集システム１３０には、光源ユニット１０１からクロックＫＣが供給される。クロックＫＣは、光源ユニット１０１において、波長可変光源により所定の波長範囲内で掃引される各波長の出力タイミングに同期して生成される。光源ユニット１０１は、例えば、各出力波長の光Ｌ０を分岐して２つの分岐光を生成し、これら分岐光の一方を光学的に遅延させ、これら分岐光を合成し、得られた合成光を検出し、その検出結果に基づいてクロックＫＣを生成する。データ収集システム１３０は、検出器１２５から入力される検出信号のサンプリングをクロックＫＣに基づいて実行する。データ収集システム１３０は、このサンプリングの結果を演算制御ユニット２００に送る。

　本例では、測定アーム長を変更するための要素（例えば、リトロリフレクタ４１）と、参照アーム長を変更するための要素（例えば、リトロリフレクタ１１４、又は参照ミラー）との双方が設けられているが、これらのうちの一方の要素のみが設けられていてもよい。また、測定アーム長と参照アーム長との間の差（光路長差）を変更するための要素はこれらの例示的な要素に限定されず、アーム長を変更する機能を実現可能な任意の要素（光学部材、機構など）であってよい。

＜制御系・処理系＞
　眼科検査装置１の制御系及び処理系の構成例を図３、図４Ａ、図４Ｂ、及び図４Ｃに示す。制御部２１０、画像データ構築部２２０、及びデータ処理部２３０は、例えば演算制御ユニット２００に設けられる。眼科検査装置１は、外部装置との間でデータ通信を行うための通信デバイスを含んでいてもよい。眼科検査装置１は、記録媒体からのデータ読み出しと、記録媒体へのデータ書き込みとを行うためのドライブ装置（リーダ／ライタ）を含んでいてもよい。

＜制御部２１０＞
　制御部２１０は、各種の制御を実行する。制御部２１０は、主制御部２１１と記憶部２１２とを含む。また、図４Ａに示すように、本実施形態において、主制御部２１１は走査制御部２１１１を含み、記憶部２１２は走査プロトコル２１２１を記憶している。

＜主制御部２１１＞
　主制御部２１１は、プロセッサを含み、眼科検査装置１の各要素（図１～図４Ｘに示された要素を含む）を制御する。主制御部２１１は、プロセッサを含むハードウェアと、制御ソフトウェアとの協働によって実現される。走査制御部２１１１は、所定形状及び所定サイズの走査エリアに対するＯＣＴスキャンの制御を行う。

　撮影合焦駆動部３１Ａは、主制御部２１１の制御の下に、撮影光路に配置された撮影合焦レンズ３１と照明光路に配置されたフォーカス光学系６０とを移動する。リトロリフレクタ（ＲＲ）駆動部４１Ａは、主制御部２１１の制御の下に、測定アームに設けられたリトロリフレクタ４１を移動する。ＯＣＴ合焦駆動部４３Ａは、主制御部２１１の制御の下に、測定アームに配置されたＯＣＴ合焦レンズ４３を移動する。リトロリフレクタ（ＲＲ）駆動部１１４Ａは、主制御部２１１の制御の下に、参照アームに配置されたリトロリフレクタ１１４を移動する。各駆動部は、主制御部２１１の制御の下に動作するパルスモータなどのアクチュエータを含む。光スキャナ４４は、主制御部２１１（走査制御部２１１１）の制御の下に動作する。

　移動機構１５０は、典型的には、眼底カメラユニット２を３次元的に移動するように構成されている。そのような移動機構１５０は、例えば、±ｘ方向（左右方向）に移動可能に構成されたｘステージと、ｘステージを移動するように構成されたｘ移動機構と、±ｙ方向（上下方向）に移動可能に構成されたｙステージと、ｙステージを移動するように構成されたｙ移動機構と、±ｚ方向（奥行き方向）に移動可能に構成されたｚステージと、ｚステージを移動するように構成されたｚ移動機構とを含む。ｘ移動機構、ｙ移動機構、及びｚ移動機構のそれぞれは、主制御部２１１の制御の下に動作するパルスモータなどのアクチュエータを含む。

＜記憶部２１２＞
　記憶部２１２は各種のデータを記憶する。記憶部２１２に記憶されるデータとしては、ＯＣＴ画像、眼底像、被検眼情報、制御情報などがある。被検眼情報は、患者ＩＤや氏名などの被検者情報や、左眼／右眼の識別情報や、電子カルテ情報などを含む。制御情報は、特定の制御に関する情報である。本実施形態の制御情報は、走査プロトコル２１２１を含む。

　走査プロトコル２１２１は、所定形状及び所定サイズの走査エリアに関するＯＣＴスキャンのための制御の内容に関する取り決めであり、各種制御パラメータ（走査制御パラメータ）の組を含む。走査プロトコル２１２１は、走査モード毎のプロトコルを含む。本実施形態の走査プロトコル２１２１は、リサジュースキャンのプロトコルを少なくとも含み、更に、例えば、Ｂスキャン（ラインスキャン）のプロトコル、クロススキャンのプロトコル、ラジアルスキャンのプロトコル、ラスタースキャンのプロトコルなどを含んでいてもよい。

　本実施形態の走査制御パラメータは、光スキャナ４４に対する制御の内容を示すパラメータを少なくとも含む。このパラメータは、例えば、スキャンパターンを示すパラメータ、スキャン速度を示すパラメータ、スキャン間隔を示すパラメータなどがある。スキャンパターンは、スキャンの経路の形状を示し、その例として、リサジューパターン、ラインパターン、クロスパターン、ラジアルパターン、ラスターパターンなどがある。スキャン速度は、例えば、Ａスキャンの繰り返しレートとして定義される。スキャン間隔は、例えば、隣接するＡスキャンの間隔として、つまりスキャン点の配列間隔として定義される。

　なお、特許文献１～４並びに非特許文献１及び２に開示されているような従来技術と同様に、本実施形態の「リサジュースキャン」は、互いに直交する２つの単振動を順序対として得られる点の軌跡が描くパターン（リサジューパターン、リサジュー図形、リサジュー曲線、リサジュー関数、バウディッチ曲線）を経路とした「狭義の」リサジュースキャンだけでなく、一連のサイクルを含む所定の２次元パターンに従う「広義の」リサジュースキャンであってもよい。

　本実施形態の光スキャナ４４は、例えば、ｘ方向に測定光ＬＳを偏向する第１ガルバノミラーと、ｙ方向に測定光ＬＳを偏向する第２ガルバノミラーとを含む。リサジュースキャンは、ｘ方向に沿った偏向方向の変化を第１周期で繰り返すように第１ガルバノミラーの制御を行いつつ、ｙ方向に沿った偏向方向の変化を第２周期で繰り返すように第２ガルバノミラーの制御を行うことによって実現される。ここで、第１周期と第２周期とは互いに異なる。

　例えば、本実施形態のリサジュースキャンは、２つの正弦波の組み合わせから得られる狭義のリサジューパターンのスキャンだけでなく、正弦波に特定項（例えば、奇数次の多項式）を加えて得られるパターンのスキャンや、三角波に基づくパターンのスキャンであってもよい。

　「サイクル」は、一般に、或る長さを持つ複数のサンプリング点から構成されるオブジェクトを意味する。本実施形態のサイクルは、例えば、閉曲線又はほぼ閉曲線（実質的閉曲線、略閉曲線）であってよい。換言すると、本実施形態のサイクルは、その始点と終点とが一致又はほぼ一致していてよい。

　典型的には、走査プロトコル２１２１はリサジュー関数に基づいて設定される。非特許文献１の式（９）及び（１０）に示すように、リサジュー関数は、例えば、次のパラメトリック方程式系で表現される：x(t_i) = A・cos(2π・(f_A/n)・t_i)、y(t_i) = A・cos(2π・(f_A・(n-2)/n²)・t_i)。

　ここで、「ｘ」はリサジュー曲線が定義される２次元座標系の横軸、「ｙ」は縦軸、「ｔ_ｉ」はリサジュースキャンにおける第ｉ番目のＡラインの収集時点、「Ａ」はスキャン範囲（振幅）、「ｆ_Ａ」はＡラインの収集レート（スキャン速度、Ａスキャンの繰り返しレート）、ｎはｘ方向（横軸方向）の各サイクルにおけるＡラインの個数をそれぞれ示す。

　このような例示的なリサジュースキャンにおけるスキャンラインの分布（スキャンパターン）の例を図５に示す。

　（狭義又は広義の）リサジュースキャンに含まれる任意のサイクル対（サイクルのペア）は、少なくとも１点（特に２点以上）で互いに交差する。このような交差を利用することで、リサジュースキャンにおける任意のサイクル対から収集されたデータ対の間のレジストレーションを行うことが可能となり、非特許文献１（又は非特許文献２）に開示された画像構築手法及びモーションアーティファクト補正手法を実装することが可能となる。以下、特に言及しない限り、主に、非特許文献１に記載された手法を適用する場合について説明するが、本実施形態において採用可能な画像構築手法及びモーションアーティファクト補正手法はこれに限定されず、例えば、非特許文献２に記載された手法を適用することや、非特許文献１又は２に記載された手法と同等及び／又は類似の手法を適用することも可能である。

　記憶部２１２に記憶される制御情報は上記の例に限定されない。例えば、制御情報はフォーカス制御を行うための情報（フォーカス制御パラメータ）を含んでいてよい。

　フォーカス制御パラメータは、ＯＣＴ合焦駆動部４３Ａに対する制御の内容を示すパラメータである。フォーカス制御パラメータの例として、測定アームの焦点位置を示すパラメータ、焦点位置の移動速度を示すパラメータ、焦点位置の移動加速度を示すパラメータなどがある。焦点位置を示すパラメータは、例えば、ＯＣＴ合焦レンズ４３の位置を示すパラメータである。焦点位置の移動速度を示すパラメータは、例えば、ＯＣＴ合焦レンズ４３の移動速度を示すパラメータである。焦点位置の移動加速度を示すパラメータは、例えば、ＯＣＴ合焦レンズ４３の移動加速度を示すパラメータである。移動速度は、一定であってもよいし、一定でなくてもよい。移動加速度についても同様である。

　このようなフォーカス制御パラメータによれば、眼底Ｅｆの形状（典型的には、中心部が深く且つ周辺部が浅い凹形状）や収差分布に応じたフォーカス調整が可能になる。フォーカス制御は、例えば、走査制御（リサジュースキャンの繰り返し制御）と連係的に実行される。それにより、モーションアーティファクトが補正され、且つ、スキャン範囲の全体に亘ってピントが合った、高品質の画像が得られる。

＜走査制御部２１１１＞
　走査制御部２１１１は、走査プロトコル２１２１に基づいて少なくとも光スキャナ４４を制御する。走査制御部２１１１は、走査プロトコル２１２１に基づく光スキャナ４４の制御と連係して光源ユニット１０１の制御を更に実行してもよい。走査制御部２１１１は、プロセッサを含むハードウェアと、走査プロトコル２１２１を含む走査制御ソフトウェアとの協働によって実現される。

＜画像データ構築部２２０＞
　画像データ構築部２２０は、プロセッサを含み、データ収集システム１３０から入力された信号（サンプリングデータ）に基づいて、眼底ＥｆのＯＣＴ画像データを構築する。このＯＣＴ画像データ構築は、従来のフーリエドメインＯＣＴ（スウェプトソースＯＣＴ）と同様に、ノイズ除去（ノイズ低減）、フィルタリング、高速フーリエ変換（ＦＦＴ）などを含む。他のタイプのＯＣＴ手法が採用される場合、画像データ構築部２２０は、そのＯＣＴタイプに対応した公知の処理を実行することによってＯＣＴ画像データを構築する。

　例えば、画像データ構築部２２０は、サンプリングデータから各スキャン点（各Ａライン）に対応する画像データ（Ａスキャン画像データ）を構築する処理を少なくとも実行するように構成されている。画像データ構築部２２０により構築される画像データは、画像化される前のデータであってもよく、例えば、深さ方向（ｚ方向、光軸方向、Ａライン方向）に沿った信号プロファイル（反射プロファイル、散乱プロファイル）であってもよい。

　前述したように、本実施形態ではリサジュースキャンが眼底Ｅｆに適用される。画像データ構築部２２０は、データ処理部２３０とともに、リサジュースキャンによる収集とデータ収集システム１３０によるサンプリングとを介して取得されたデータセットに対し、例えば非特許文献１に開示された画像構築手法及びモーションアーティファクト補正手法を適用することによって、眼底Ｅｆの３次元画像データを構築することができる。

　画像データ構築部２２０及び／又はデータ処理部２３０は、３次元画像データにレンダリングを適用して表示用画像を形成することができる。適用可能なレンダリング法の例として、ボリュームレンダリング、サーフェスレンダリング、最大値投影（ＭＩＰ）、最小値投影（ＭｉｎＩＰ）、多断面再構成（ＭＰＲ）などがある。

　画像データ構築部２２０及び／又はデータ処理部２３０は、３次元画像データに基づいてＯＣＴ正面画像（en face OCT image）を構築することが可能である。例えば、画像データ構築部２２０及び／又はデータ処理部２３０は、３次元画像データをｚ方向（Ａライン方向、光軸方向、深さ方向）に投影することでプロジェクションデータを構築することができる。また、画像データ構築部２２０及び／又はデータ処理部２３０は、３次元画像データの一部である部分的３次元画像データをｚ方向に投影することでシャドウグラムを構築することができる。この部分的３次元画像データは、例えば、任意のセグメンテーション手法を利用して設定されてよい。セグメンテーションは、画像中の部分領域を特定する処理である。本例では、眼底Ｅｆの１つ以上の組織（部位）に相当する画像領域を特定するためにセグメンテーションを行うことができる。

　シャドウグラム構築処理に限らず、実施形態のセグメンテーションは、例えば、閾値処理、エッジ検出、フィルタリング、機械学習（例えば、セマンティックセグメンテーション）など、任意の公知の技術を利用可能である。

　眼科検査装置１は、ＯＣＴ血管造影（OCT-Angiography）を実行可能であってよい。ＯＣＴ血管造影は、血管が強調された画像を構築するイメージング技術である（例えば、非特許文献２、特表２０１５－５１５８９４号公報などを参照）。一般に、眼底組織（構造）は短時間の間には変化しないが、血管内部の血流は短時間の間にも変化する。ＯＣＴ血管造影では、このような時間的変化が存在する部分（血流信号）を強調して画像を生成する。なお、ＯＣＴ血管造影は、ＯＣＴモーションコントラスト撮影（motion contrast imaging）などとも呼ばれる。また、ＯＣＴ血管造影により取得される画像は、血管造影画像、アンジオグラム（angiogram）、モーションコントラスト画像などと呼ばれる。

　ＯＣＴ血管造影を実行可能である場合、眼科検査装置１は、眼底Ｅｆの同じ領域を所定回数だけ繰り返しスキャンする。例えば、眼科検査装置１は、前述した走査制御（リサジュースキャンの繰り返し制御）を所定回数だけ繰り返し実行する。それにより、リサジュースキャンの適用領域から複数の３次元データ（時系列な３次元データセット）がデータ収集システム１３０によって収集される。画像データ構築部２２０及び／又はデータ処理部２３０は、この３次元データセットからモーションコントラスト画像を構築することができる。このモーションコントラスト画像は、眼底Ｅｆの血流に起因する干渉信号の時間的変化が強調された血管造影画像である。この血管造影画像は、眼底Ｅｆの血管の３次元的な分布を表現した３次元血管造影画像データである。

　画像データ構築部２２０及び／又はデータ処理部２３０は、この３次元血管造影画像データから、任意の２次元血管造影画像データ及び／又は任意の擬似的３次元血管造影画像データを構築することが可能である。例えば、画像データ構築部２２０及び／又はデータ処理部２３０は、３次元血管造影画像データに多断面再構成を適用することにより、眼底Ｅｆの任意の断面を表す２次元血管造影画像データを構築することができる。また、画像データ構築部２２０及び／又はデータ処理部２３０は、３次元血管造影画像データにプロジェクション画像化又はシャドウグラム化を適用することにより、眼底Ｅｆの正面血管造影画像データを構築することができる。

　本実施形態において、画像データ構築部２２０及び／又はデータ処理部２３０は、データ収集システム１３０により収集されたデータから複数のストリップを構築する。非特許文献１に記載されているように、画像データ構築部２２０及び／又はデータ処理部２３０は、リサジュースキャンで収集されたボリューム（３次元データ）を比較的大きな動きが介在しない複数のサブボリュームに分割し、各サブボリュームの正面プロジェクション画像（en face projection）を構築する。この正面プロジェクション画像がストリップである。このようにして得られた複数のストリップに対してレジストレーション及びマージ処理を適用することにより、（ｘｙ方向の）モーションアーティファクトを補正することができる。後述するように、本例においては、レジストレーション及びマージ処理はデータ処理部２３０によって実行される。

　画像データ構築部２２０は、プロセッサを含むハードウェアと、画像構築ソフトウェアとの協働によって実現される。なお、幾つかの例示的な態様において、画像データ構築部２２０とデータ処理部２３０とが一体的に構成されていてよい。

＜データ処理部２３０＞
　データ処理部２３０は、プロセッサを含み、被検眼Ｅの画像に対して各種のデータ処理を適用する。例えば、データ処理部２３０は、プロセッサを含むハードウェアと、データ処理ソフトウェアとの協働によって実現される。データ処理部２３０に含まれる各機能要素も同様である。データ処理部２３０が実行するデータ処理は被検眼Ｅの画像の処理には限定されず、任意のデータ処理であってよい。図４Ｂに示すように、本実施形態のデータ処理部２３０は、画像選択部２２５と、画像データ処理部２３１とを含んでいる。

　データ処理部２３０は、眼底Ｅｆについて取得された２つの画像の間の位置合わせ（レジストレーション）を行うように構成されてよい。例えば、データ処理部２３０は、ＯＣＴで取得された３次元画像データと、眼底カメラユニット２により取得された正面画像との間のレジストレーションを行うように構成される。また、データ処理部２３０は、ＯＣＴで取得された２つのＯＣＴ画像の間のレジストレーションを行うように構成される。また、データ処理部２３０は、眼底カメラユニット２により取得された２つの正面画像の間のレジストレーションを行うように構成される。また、ＯＣＴ画像の解析結果や、正面画像の解析結果に対してレジストレーションを適用するように構成されてもよい。これらのレジストレーションは、公知の手法によって実行可能であり、例えば特徴点抽出とアフィン変換とを含む。

　データ処理部２３０は、リサジュースキャンを用いて得られたデータセットを処理するように構成される。前述したように、画像データ構築部２２０（及び／又はデータ処理部２３０）は、例えば、リサジュースキャンで収集されたデータから複数のストリップを構築する。データ処理部２３０は、これらストリップに対してレジストレーション及びマージ処理を適用することにより、モーションアーティファクトが補正された画像を構築するように構成されている。

　前述したようにリサジュースキャンはサンプルのデータを冗長的に収集するため、任意の２つのストリップは交差領域（共通領域、オーバーラップ領域、重複領域）を有する。非特許文献１に記載された手法と同様に、データ処理部２３０は、ストリップ同士のオーバーラップを利用してストリップ間のレジストレーションを行うように構成されてよい。例えば、データ処理部２３０は、まず、画像データ構築部２２０により構築された複数のストリップを大きさ（面積など）に従って順序付け（第１～第Ｎのストリップ）、最大のストリップである第１のストリップを初期基準ストリップに指定することができる。初期基準ストリップは、初期基準画像の例である。次に、データ処理部２３０は、第１のストリップを基準として第２のストリップのレジストレーションを行い、第１のストリップと第２のストリップとをマージする（合成する）。データ処理部２３０は、これにより得られたマージストリップを基準ストリップとして第３のストリップのレジストレーションを行い、このマージストリップと第３のストリップとをマージする。このようなレジストレーション及びマージ処理を上記順序に従って順次に実行することにより第１～第Ｎのストリップの位置合わせ及び貼り合わせがなされ、モーションアーティファクトが補正された画像が得られる。

　従来の技術では、このような再帰的なレジストレーション（及びマージ処理）を単一の初期基準画像のみに基づき実行しているため、前述したように最終画像の品質が劣化するおそれがある。これに対し、本実施形態は、複数の初期基準画像にそれぞれ基づく複数の再帰的なレジストレーション（及びマージ処理）を実行することを１つの特徴とするものであり、その結果として、最終画像の品質向上という１つのアドバンテージを達成するものである。このような特徴及びアドバンテージを実現するための本実施形態の構成の例を以下に説明する。

＜画像選択部２２５＞
　本態様の眼科検査装置１は、リサジュースキャンによって冗長的データ収集を実行し、それにより取得されたデータセットに基づいて画像セットを生成する。本態様では、この画像セットは、画像データ構築部２２０（又はデータ処理部２３０）によって生成された複数のストリップを含む。前述したように、複数のストリップは、位置的な冗長性を有している。つまり、任意の２つのストリップ（任意のストリップペア）は、それらの交差領域において位置的な冗長性を有している。

　画像選択部２２５は、生成された複数のストリップから２以上のストリップを選択する。選択された各ストリップは、画像セット（複数のストリップ）の相対位置調整を行うためのレジストレーションの基準として使用される。本態様では、選択された各ストリップは、再帰的レジストレーションにおける初期基準ストリップとして使用される。すなわち、選択された各ストリップは、再帰的レジストレーションにおいて繰り返し実行されるレジストレーションのうちの第１回目のレジストレーションで基準画像として使用される。

　このように、本態様では、画像選択部２２５によって選択された複数の初期基準ストリップを用いて複数の再帰的レジストレーションがそれぞれ実行される。なお、従来の技術では、単一の初期基準ストリップを用いた単一回数の再帰的レジストレーションが実行される。

　画像選択部２２５が実行する初期基準画像選択処理の幾つかの例を以下に説明する。なお、画像選択部２２５が実行する処理は、これらの例に限定されるものではなく、冗長性を有する画像セットから複数の初期基準画像を選択するために適用可能な任意の処理であってよい。

　初期基準画像選択処理の第１の例を説明する。本例の処理を実行可能な画像選択部２２５は、光走査（リサジュースキャン）の時間軸に基づいて画像セット（複数のストリップ）から複数の初期基準画像（複数の初期基準ストリップ）を選択するように構成される。

　本実施形態の光走査は、複数の軌道を順次にスキャンするものであり、複数の軌道は時系列に配列されているとみなすことができ、更に、複数の軌道にそれぞれ対応する複数の画像（画像セット）は時系列に配列されているとみなすことができる。特に、リサジュースキャンは、複数のサイクルを順次にスキャンするスキャンモードであり、複数のサイクルは時系列に配列されているとみなすことができ、更に、複数のサイクルにそれぞれ対応する複数のストリップは時系列に配列されているとみなすことができる。非特許文献１の式（９）及び（１０）からも分かるように、典型的なリサジュースキャンは、時間パラメータによってパラメトライズされている。一般的な光走査についても同様である。第１の例は、このような光走査の時間軸に基づく初期基準画像選択に関する。

　図６を参照して、光走査の時間軸に基づく初期基準画像選択の具体例を説明する。図６の座標軸ｔは時間軸を表す。時間軸ｔには、時系列順に、Ｋ個の区間Ａ１、Ａ２、・・・、ＡＫが設定されている。ここで、Ｋは２以上の整数である。時間的に最初の区間Ａ１の左端点ＴＳはスキャン開始時間に相当し、時間的に最後の区間ＡＫの右端点ＴＥはスキャン終了時間に相当する。

　区間Ａ１に相当する期間に収集されたデータセットから構築された画像群（ストリップ群）を符号Ｂ１_ｎ１で表し、区間Ａ２に相当する期間に収集されたデータセットから構築された画像群（ストリップ群）を符号Ｂ２_ｎ２で表し、区間ＡＫに相当する期間に収集されたデータセットから構築された画像群（ストリップ群）を符号ＢＫ_ｎＫで表す。ここで、ｎ１＝１、２、・・・、Ｎ１であり、ｎ２＝１、２、・・・、Ｎ２であり、ｎＫ＝１、２、・・・、ＮＫであり、Ｎ１、Ｎ２、・・・、ＮＫはいずれも正の整数（典型的には２以上の整数）である。リサジュースキャンで得られた画像（ストリップ）の個数をＮとすると、Ｎ１＋Ｎ２＋・・・＋ＮＫ＝Ｎである。

　Ｋ個の区間Ａ１、Ａ２、・・・、ＡＫのそれぞれの長さは等しくてもよいし、等しくなくてもよい。換言すると、スキャン開始時間ＴＳからスキャン終了時間ＴＥまでの期間を同じ長さの区間に分割してもよいし、そうでなくてもよい。また、Ｎ１、Ｎ２、・・・、ＮＫは等しくてもよいし、等しくなくてもよい。

　本例の画像選択部２２５は、まず、リサジュースキャンで得られたＮ個のストリップＢ１_１～ＢＫ_ＮＫからなる画像セットを、光走査の時間軸ｔのＫ個の区間Ａ１～ＡＫに対応するＫ個のサブセット（Ｎ１個のストリップＢ１_１～Ｂ１_Ｎ１；Ｎ２個のストリップＢ２_１～Ｂ２_Ｎ２；・・・；ＮＫ個のストリップＢＫ_１～ＢＫ_ＮＫ）に分割する。

　次に、本例の画像選択部２２５は、Ｋ個のサブセットのそれぞれからストリップを選択する。なお、Ｋ個のサブセットの全てからストリップを選択してもよいし、Ｋ個のサブセットのうちのいずれかのサブセットからストリップを選択しなくてもよい。また、１つのサブセットから選択されるストリップの個数は任意であってよい。また、Ｋ個のサブセットからそれぞれ同数のストリップを選択してもよいし、そうでなくてもよい。いずれにしても、本例の画像選択部２２５は、時系列に配列され且つＫ個のサブセットに分割されたＮ個のストリップのうちから、２つ以上のストリップを選択する。選択された各ストリップは、初期基準画像として用いられる。

　本例の画像選択部２２５は、ストリップの大きさ（画像サイズ。例えば面積。）に基づいて、Ｋ個のサブセットのそれぞれからストリップを選択するように構成されてよい。つまり、ストリップの選択基準は、画像サイズを含んでいてもよい。

　例えば、本例の画像選択部２２５は、Ｋ個のサブセットのそれぞれから、そのサブセットにおける最大サイズのストリップ（そのサブセットにおいて画像サイズが最大のストリップ）を選択するように構成されていてよい。なお、サブセットにおける最大サイズのストリップのみを選択してもよいし、最大サイズのストリップを含む２以上のストリップを選択してもよい。

　本例の画像選択部２２５は、Ｋ個のサブセットのそれぞれについて、そのサブセットに含まれる複数のストリップの全部又は一部に対して画像サイズに基づく順序付けを行うように構成されていてよい。例えば、本例の画像選択部２２５は、各ストリップの画像サイズを求める処理と、求められた画像サイズを比較する処理と、この画像サイズ比較の結果に基づいて複数のストリップに順序を割り当てる処理とを含む。本例の画像選択部２２５は、更に、Ｋ個のサブセットのそれぞれから、そのサブセット中の複数のストリップに付された順序に基づいてストリップの選択を行うように構成されていてよい。これにより、画像サイズに応じた順序付けにしたがって１つ又は２つ以上のストリップをサブセットから選択することができる。

　詳細は後述するが、本例の画像選択部２２５は、サブセットから選択されたストリップを初期基準画像として実行されたレジストレーションの結果に基づいて、同じサブセットから別のストリップを選択するように構成されていてよい。例えば、本例の画像選択部２２５は、Ｋ個のサブセットのうちの１つのサブセットから、このサブセット中の複数のストリップに付された順序に基づいて、１つのストリップ（第１画像）を選択することができる。選択されたストリップを基準としたレジストレーションが実行された後、本例の画像選択部２２５は、このレジストレーションの成否を判定することができる。このレジストレーションは成功であると判定された場合、このサブセットに基づくレジストレーションは完了と判断される。一方、このレジストレーションは失敗であると判定された場合、本例の画像選択部２２５は、このサブセット中の複数のストリップに付された順序に基づいて、このサブセットから別のストリップ（第２画像）を選択することができる。

　初期基準画像選択処理の第２の例を説明する。本例の処理を実行可能な画像選択部２２５は、光走査（リサジュースキャン）の適用エリアに基づいて画像セット（複数のストリップ）から複数の初期基準画像（複数の初期基準ストリップ）を選択するように構成される。画像セットを時間的に分割している第１の例に対して、本例は画像セットを空間的に分割したものと言える。

　図７、図８及び図８Ｂを参照して、光走査の適用エリアに基づく初期基準画像選択の具体例を説明する。

　図７の具体例において、符号Ｃはリサジュースキャンの適用エリアを示す（例えば、図５に示すスキャンエリア）。本具体例では、適用エリアＣをｘ方向（左右方向）に二等分するとともにｙ方向（上下方向）に二等分することで、４つのサブエリアＣ１、Ｃ２、Ｃ３、及びＣ４が形成されている。

　本例の画像選択部２２５は、リサジュースキャンで得られた画像セット（複数のストリップ）を４つのサブエリアＣ１～Ｃ４に対応する４つのサブセットに分割する。サブエリアＣ１に対応するサブセットには、サブエリアＣ１中の画像（例えば、ストリップの一部（部分ストリップ））が含まれる。サブエリアＣ２～Ｃ４のそれぞれに対応するサブセットについても同様である。

　更に、本例の画像選択部２２５は、４つのサブエリアＣ１～Ｃ４に対応する４つのサブセットのそれぞれから画像（例えば、部分ストリップ）を選択することによって、画像セット（複数のストリップ）から複数の画像（例えば、複数の部分ストリップ）を選択する。選択された各部分ストリップは、初期基準画像として用いられる。

　本例の画像選択部２２５は、部分ストリップの大きさ（画像サイズ。例えば面積。）に基づいて、４つのサブセットのそれぞれから部分ストリップを選択するように構成されてよい。つまり、部分ストリップの選択基準は、画像サイズを含んでいてもよい。

　例えば、本例の画像選択部２２５は、４つのサブセットのそれぞれから、そのサブセットにおける最大サイズの部分ストリップ（そのサブセットにおいて画像サイズが最大の部分ストリップ）を選択するように構成されていてよい。なお、サブセットにおける最大サイズの部分ストリップのみを選択してもよいし、最大サイズの部分ストリップを含む２以上の部分ストリップを選択してもよい。

　本例の画像選択部２２５は、４つのサブセットのそれぞれについて、そのサブセットに含まれる複数の部分ストリップの全部又は一部に対して画像サイズに基づく順序付けを行うように構成されていてよい。例えば、本例の画像選択部２２５は、各部分ストリップの画像サイズを求める処理と、求められた画像サイズを比較する処理と、この画像サイズ比較の結果に基づいて複数の部分ストリップに順序を割り当てる処理とを含む。本例の画像選択部２２５は、更に、４つのサブセットのそれぞれから、そのサブセット中の複数の部分ストリップに付された順序に基づいて部分ストリップの選択を行うように構成されていてよい。これにより、画像サイズに応じた順序付けにしたがって１つ又は２つ以上の部分ストリップをサブセットから選択することができる。

　本例の画像選択部２２５は、サブセットから選択された部分ストリップを初期基準画像として実行されたレジストレーションの結果に基づいて、同じサブセットから別の部分ストリップを選択するように構成されていてよい。例えば、本例の画像選択部２２５は、４つのサブセットのうちの１つのサブセットから、このサブセット中の複数の部分ストリップに付された順序に基づいて、１つの部分ストリップ（第１画像）を選択することができる。選択された部分ストリップを基準としたレジストレーションが実行された後、本例の画像選択部２２５は、このレジストレーションの成否を判定することができる。このレジストレーションは成功であると判定された場合、このサブセットに基づくレジストレーションは完了と判断される。一方、このレジストレーションは失敗であると判定された場合、本例の画像選択部２２５は、このサブセット中の複数の部分ストリップに付された順序に基づいて、このサブセットから別の部分ストリップ（第２画像）を選択することができる。

　図７の具体例では、リサジュースキャンの適用エリアＣを形状及びサイズが等しい４つのサブエリアＣ１～Ｃ４に分割しているが、適用エリアの分割態様はこれに限定されない。例えば、適用エリアを分割して形成されるサブエリアの個数は２以上の任意の個数であってよい。また、適用エリアを分割して形成される２以上のサブエリアは、異なる形状及び／又は異なるサイズを有していてもよい。

　適用エリアの分割態様は、図７の具体例のような格子状分割に限定されず、例えば放射状分割、同心円状分割など、任意の分割態様であってよい。また、適用エリアの分割態様は、格子状分割、放射状分割、同心円状分割のような幾何学的な分割態様に限定されない。

　幾つかの例示的な態様では、適用エリアに描出されているオブジェクトに基づいて適用エリアを分割することができる。例えば、適用エリアに第１オブジェクト及び第２オブジェクトが描出されている場合、第１オブジェクトの像を含むサブエリアと第２オブジェクトを含むサブエリアとを含む２以上のサブエリアに適用エリアを分割することができる。眼底Ｅｆのリサジュースキャンの場合、このリサジュースキャンの適用エリアを視神経乳頭の像を含むサブエリアと黄斑を含むサブエリアとを含む２以上のサブエリアに分割することができる。

　幾つかの例示的な態様では、適用エリアの少なくとも一部に対して重要度の分布を割り当てることができ、この重要度の分布に基づいて適用エリアを２以上のサブエリアに分割することができる。例えば、重要度の高い部分を含むサブエリアとそれ以外のサブエリアとを含む２以上のサブエリアに適用エリアを分割することができる。眼底Ｅｆのリサジュースキャンの場合、視神経乳頭の像、黄斑の像、病変部の像などを重要度の高い部分に指定することができる。

　図８Ａ及び図８Ｂの具体例は、特開２０２１－４０８５４号公報に開示されたスキャンモード（コンボリューショナルリサジュースキャンと呼ぶ）への応用である。コンボリューショナルリサジュースキャンは、図５のような標準的なリサジュースキャンとその適用エリアの移動（適用エリアシフト）とを組み合わせたスキャンモードである（図８Ａを参照）。適用エリアシフトは、図８Ａの例では円形の軌道に沿って実行されるが、軌道はこれに限定されず、任意であってよい。コンボリューショナルリサジュースキャンは、モーションアーティファクト補正と撮影の広角化（適用エリアの拡大）との両立を可能にする。

　図８Ｂの具体例において、符号Ｄはコンボリューショナルリサジュースキャンの適用エリアを示す。本具体例では、適用エリアＤをその中心周りに定義される角度方向に三等分
することで、３つのサブエリアＤ１、Ｄ２、及びＤ３が形成されている。

　本例の画像選択部２２５は、リサジュースキャンで得られた画像セット（複数のストリップ）を３つのサブエリアＤ１～Ｄ３に対応する３つのサブセットに分割する。サブエリアＤ１に対応するサブセットには、サブエリアＤ１中の画像（例えば、部分ストリップ）が含まれる。サブエリアＤ２及びＤ３のそれぞれに対応するサブセットについても同様である。

　更に、本例の画像選択部２２５は、３つのサブエリアＤ１～Ｄ３に対応する３つのサブセットのそれぞれから画像（例えば、部分ストリップ）を選択することによって、画像セット（複数のストリップ）から複数の画像（例えば、複数の部分ストリップ）を選択する。選択された各部分ストリップは、初期基準画像として用いられる。

　本例の画像選択部２２５は、部分ストリップの大きさ（画像サイズ。例えば面積。）に基づいて、３つのサブセットのそれぞれから部分ストリップを選択するように構成されてよい。つまり、部分ストリップの選択基準は、画像サイズを含んでいてもよい。

　例えば、本例の画像選択部２２５は、３つのサブセットのそれぞれから、そのサブセットにおける最大サイズの部分ストリップ（そのサブセットにおいて画像サイズが最大の部分ストリップ）を選択するように構成されていてよい。なお、サブセットにおける最大サイズの部分ストリップのみを選択してもよいし、最大サイズの部分ストリップを含む２以上の部分ストリップを選択してもよい。

　本例の画像選択部２２５は、３つのサブセットのそれぞれについて、そのサブセットに含まれる複数の部分ストリップの全部又は一部に対して画像サイズに基づく順序付けを行うように構成されていてよい。例えば、本例の画像選択部２２５は、各部分ストリップの画像サイズを求める処理と、求められた画像サイズを比較する処理と、この画像サイズ比較の結果に基づいて複数の部分ストリップに順序を割り当てる処理とを含む。本例の画像選択部２２５は、更に、３つのサブセットのそれぞれから、そのサブセット中の複数の部分ストリップに付された順序に基づいて部分ストリップの選択を行うように構成されていてよい。これにより、画像サイズに応じた順序付けにしたがって１つ又は２つ以上の部分ストリップをサブセットから選択することができる。

　本例の画像選択部２２５は、サブセットから選択された部分ストリップを初期基準画像として実行されたレジストレーションの結果に基づいて、同じサブセットから別の部分ストリップを選択するように構成されていてよい。例えば、本例の画像選択部２２５は、３つのサブセットのうちの１つのサブセットから、このサブセット中の複数の部分ストリップに付された順序に基づいて、１つの部分ストリップ（第１画像）を選択することができる。選択された部分ストリップを基準としたレジストレーションが実行された後、本例の画像選択部２２５は、このレジストレーションの成否を判定することができる。このレジストレーションは成功であると判定された場合、このサブセットに基づくレジストレーションは完了と判断される。一方、このレジストレーションは失敗であると判定された場合、本例の画像選択部２２５は、このサブセット中の複数の部分ストリップに付された順序に基づいて、このサブセットから別の部分ストリップ（第２画像）を選択することができる。

　図８Ｂの具体例では、リサジュースキャンの適用エリアＤを３つに等角度分割して３つのサブエリアＤ１～Ｄ３を形成しているが、図７の具体例に関して説明したように、適用エリアの分割態様はこれに限定されない。

　初期基準画像選択処理の第３の例を説明する。本例の処理を実行可能な画像選択部２２５は、光走査（リサジュースキャン）で得られた画像セット（複数のストリップ）に含まれる全ての画像又は一部のみの画像に対して画像サイズに基づく順序付けを行い、この順序付けにより画像に付された順序に基づいて画像セットから複数の初期基準画像（複数の初期基準ストリップ）を選択するように構成される。

　本例は、第１の例及び第２の例のように画像セットを分割するのではなく、画像セットの全体から複数の初期基準画像を選択するものである。なお、幾つかの例示的な態様では、画像選択部２２５は、光走査で取得された画像セットから所定の条件を満足する１以上の画像を除外し、残った画像セットから複数の初期基準画像を選択するように構成されてもよい。この条件は、例えば、各画像に関する条件（画像品質、信号品質など）や、画像間の関係に関する条件（画像相関、位置ずれ、画像品質の違い、信号品質の違いなど）であってよい。

　本例の画像選択部２２５は、画像サイズを参照した順序に基づき画像セットから選択された複数の初期基準画像のうちの１つの画像を基準として実行されたレジストレーションの結果に基づいて、この画像セットから別の初期基準画像（先に選択された複数の初期基準画像のいずれとも異なる画像）を選択するように構成されていてよい。例えば、本例の画像選択部２２５は、まず、画像サイズに基づく順序付けを画像セット（複数のストリップ）に対して行い、付与された順序に基づいてこの画像セットから複数の初期基準ストリップ（初期基準ストリップ群）を選択する。この初期基準ストリップ群のうちの１つの初期基準ストリップに基づくレジストレーションが実行された後、本例の画像選択部２２５は、このレジストレーションの成否を判定することができる。このレジストレーションは成功であると判定された場合、この初期基準ストリップは適当であると判断される。一方、このレジストレーションは失敗であると判定された場合、この初期基準ストリップは不適当であると判断され、本例の画像選択部２２５は、画像セットに付された順序に基づいて、初期基準ストリップ群に含まれる複数の画像のいずれとも異なるストリップをこの画像セットから新たに選択する。これにより、最初に選択した初期基準ストリップ群のうち不適当なストリップを新たなストリップに置き換えることができ、最終的には、モーションアーティファクト補正のために適当なストリップのみからなる初期基準ストリップ群を得ることができる。

　初期基準画像選択処理の第４の例を説明する。本例は、光走査が適用されているときのサンプルの動き（運動）を参照して初期基準画像の選択を行うものである。このサンプルの動きは、モーションアーティファクトの原因となる運動又はなり得る運動である。眼科分野においては、眼球運動、頭部の動き、体動などが、サンプル（被検眼Ｅ）の動きに含まれる。上記した初期基準画像選択の第１～第３の例又は他の例に、本例を組み合わせることが可能である。例えば、第１～第３の例のいずれかの処理の前処理として本例の処理を実行することができる。

　本例の処理を実行可能な眼科検査装置１は、光走査がサンプルに適用されているときの被検眼Ｅの運動を表す運動情報を生成する運動情報生成部を備えている。運動情報生成部の構成は任意であってよいが、例えば以下に説明する構成のいずれかを含んでいてもよい。なお、以下の例示的な構成は、眼の運動に関するものであるが、他のサンプルの

　運動情報生成部の第１の構成例は、ＯＣＴ血管造影信号（ＯＣＴＡ信号、モーションコントラスト信号）から被検眼Ｅの運動を検出するものである。非特許文献２には、リサジュースキャンを用いたＯＣＴ血管造影において、ＯＣＴ血管造影信号から眼球運動を検出すること、及び、眼球運動の影響を受けたデータ（繰り返しスキャンサイクルデータ、repeat-cycle-set）を破棄することが記載されている。

　本構成例の運動情報生成部が適用される場合、眼科検査装置１は、例えば、眼底カメラユニット２、ＯＣＴユニット１００、及び制御部２１０によってＯＣＴ血管造影のためのリサジュースキャンを実行してデータセットを収集し、画像データ構築部２２０及びデータ処理部２３０によってこのデータセットを処理してＯＣＴ血管造影信号を生成し、画像選択部２２５によってこのＯＣＴ血管造影信号から眼球運動を検出し、画像選択部２２５によって所定の条件を満たす眼球運動が検出された時点の近傍期間又は位置の近傍領域のストリップを初期基準ストリップの候補から除外するように構成されてよい。この条件は任意であってよいが、例えば、眼球運動の大きさが所定の閾値以上であるという条件、及び／又は、眼球運動の発生頻度が所定の閾値以上であるという条件であってよい。

　運動情報生成部の第２の構成例は、ＯＣＴ信号以外のデータから被検眼Ｅの運動を検出するものである。本構成例の運動検出に利用可能なデータとしては、眼底カメラユニット２により取得される観察画像、特開２０１３－２４８３７６号公報の前眼部カメラにより取得される前眼部画像などがある。本構成例においては、第１の構成例では必要とされるＯＣＴ血管造影のための繰り返しスキャンが不要である。

　本構成例の運動情報生成部が適用される場合、眼科検査装置１は、例えば、眼底カメラユニット２又は図示しない前眼部カメラによって被検眼Ｅの画像を取得し、データ処理部２３０（例えば画像選択部２２５）によってこの画像から眼球運動を検出し、画像選択部２２５によって所定の条件を満たす眼球運動が検出された時点の近傍期間又は位置の近傍領域のストリップを初期基準ストリップの候補から除外するように構成されてよい。この条件は任意であってよいが、例えば、眼球運動の大きさが所定の閾値以上であるという条件、及び／又は、眼球運動の発生頻度が所定の閾値以上であるという条件であってよい。

　このように、リサジュースキャンが被検眼Ｅに適用されているときの被検眼Ｅの運動を表す運動情報を生成する運動情報生成部を眼科検査装置１が備えている場合、画像選択部２２５は、運動情報生成部により生成された運動情報に基づいて、画像セット（複数のストリップ）から複数の初期基準画像を選択するように構成されていてよい。

　更に、画像選択部２２５は、運動情報生成部により生成された運動情報から所定の運動パラメータの値を算出する処理と、算出された値を評価する処理と、この評価の結果に基づいて画像セット（複数のストリップ）から複数の初期基準画像を選択する処理とを実行するように構成されていてよい。

　運動情報の評価に用いられる運動パラメータは、眼球運動の大きさ及び眼球運動の発生頻度の少なくとも一方を含んでいてよい。更に、画像選択部２２５は、このような運動パラメータについて算出された値を所定の閾値と比較することができる。運動パラメータの値が閾値以上であると判定された場合、画像選択部２２５は、この運動パラメータの値に対応する期間を特定し、特定された期間中にリサジュースキャンで収集されたデータに基づく画像（ストリップ）を特定し、特定された１以上の画像が画像セットから除外されたサブセットから複数の初期基準画像を選択することができる。ここで、運動パラメータの値に対応する期間は、例えば、閾値以上の大きさの眼球運動が発生した時点を含む所定の長さの期間、及び／又は、閾値以上の頻度の眼球運動が発生した時間を含む所定の長さの期間であってよい。

　初期基準画像選択処理の第４の例によれば、トラッキング等の眼球運動検出機能を初期基準画像選択に組み合わせることができる。この組み合わせにより、サッカード（saccade）等の眼球運動が発生した付近においてレジストレーション方法を切り替えることができる。例えば、眼球運動が発生した付近のデータをレジストレーションから除外することができる。また、サッカード等の眼球運動が発生した付近の画像（ストリップ）を初期基準画像の候補から除外することができる。また、サッカード等の眼球運動が頻発した付近の画像（ストリップ）を初期基準画像の候補から除外することができる。

　上記した幾つかの例示的な初期基準画像選択処理では、画像サイズを参照する場合について説明したが、初期基準画像選択において参照される情報（パラメータ）は画像サイズに限定されない。例えば、実施形態は、任意の画質評価指標に基づく順序付けを画像セット（複数のストリップ）に適用する態様や、画像（ストリップ）の任意の特徴量の多寡に基づく順序付けを画像セット（複数のストリップ）に適用する態様や、任意の画質評価指標と任意の画像特徴量とを組み合わせた指標に基づく順序付けを画像セット（複数のストリップ）に適用する態様や、任意の画質評価指標及び／又は任意の画像特徴量と画像サイズとを組み合わせた指標に基づく順序付けを画像セット（複数のストリップ）に適用する態様や、これら以外の指標に基づく順序付けを画像セット（複数のストリップ）に適用する態様を含んでいてよい。

＜画像データ処理部２３１＞
　画像データ処理部２３１の構成及びその背景について説明する。モーションアーティファクト補正のための再帰的なレジストレーションでは、基準ストリップと他のストリップとの間の相対位置を求めるために相互相関関数が利用されるが、ストリップは任意の形状を有するため、各ストリップを特定形状（例えば正方形状）の画像として扱うためのマスクを用いてストリップ間の相関演算が行われる（非特許文献１のAppendix Aを参照）。

　前述したように、ストリップは所定階調の強度画像であり、マスク画像は二値画像であるため、ストリップの画素値の絶対値とマスクの画素値の絶対値との差が大きくなり、ストリップとマスクとを用いた相関演算におけるマスクの効果が無視されて正しい相関係数が得られないおそれがある。特に、コストなどの観点から、ストリップ間の相関係数を求めるために単精度浮動小数点型（フロート型）の演算を採用した場合、有効桁数の少なさに起因する丸め誤差の影響により、この問題が顕著に現れる。

　このような丸め誤差の問題への対処に加えて、複数の初期基準画像に基づく画像セットの相対位置調整を実現するために、画像データ処理部２３１は図４Ｃに示す要素群を備えている。なお、複数の初期基準画像に基づく画像セットの相対位置調整を実行するための構成は、図４Ｃに示す構成に限定されるものではなく、以下に示す例示的な態様による相対位置調整又は同様の相対位置調整を実行可能な任意の構成であってもよい。

　画像データ処理部２３１は、画像データ構築部２２０により構築された画像データを処理する。図４Ｃに示すように、画像データ処理部２３１は、マスク画像生成部２３１１と、範囲調整部２３１２と、合成画像生成部２３１３と、相互相関関数算出部２３１４と、相関係数算出部２３１５と、ｘｙシフト量算出部２３１６と、レジストレーション部２３１７と、マージ処理部２３１８と、ｚシフト量算出部２３２と、シフト量記録部２３３と、調整量算出部２３４と、画像データ補正部２３５とを含んでいる。

　以下に説明する例では、画像データ処理部２３１は、任意の２つのストリップの一方を基準ストリップ（reference strip）に指定するとともに、この基準ストリップに対して他方のストリップ（対象ストリップ；registering strip）をレジストレーションするように構成されている。

＜マスク画像生成部２３１１＞
　マスク画像生成部２３１１は、基準ストリップに応じた基準マスク画像と、対象ストリップに応じた対象マスク画像とを生成する。マスク画像の幾つかの例を以下に説明するが、これらに限定されるものではない。

　マスク画像の輪郭形状は、例えば矩形であり、典型的には正方形である。基準マスク画像の形状と対象マスク画像の形状とは同じであってよい。また、基準マスク画像のサイズと対象マスク画像のサイズとは同じであってよい。

　マスク画像の画素値の範囲は、例えば閉区間［０，１］に含まれるように設定される。典型的には、マスク画像は、画素値が０又は１の二値画像であってよい。具体例として、マスク画像は、ストリップの定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像である。換言すると、非特許文献１の式（２０）に示すように、例示的なマスク画像の画素値は、対応するストリップの画像エリアにおける値が１であり、他のエリアにおける値が０である。

＜範囲調整部２３１２＞
　範囲調整部２３１２は、基準マスク画像の画素値の範囲と、対象マスク画像の画素値の範囲とに基づいて、基準ストリップの画素値の範囲と、対象ストリップの画素値の範囲とを調整するように構成される。典型的には、範囲調整部２３１２は、ストリップの画素値の範囲とマスク画像の画素値の範囲との間の差を小さくするように、基準ストリップの画素値の範囲及び対象ストリップの画素値の範囲の調整を行う。

　一般に、範囲調整部２３１２は、基準ストリップの画素値の範囲と、対象ストリップの画素値の範囲と、マスク画像の画素値の範囲とを相対的に調整する。典型的な例では、基準ストリップに適用される基準マスク画像と、対象ストリップに適用される対象マスク画像とは異なり、範囲調整部２３１２は、基準ストリップの画素値の範囲と、対象ストリップの画素値の範囲と、基準マスク画像の画素値の範囲と、対象ストリップの画素値の範囲とを相対的に調整する。

　範囲調整部２３１２は、基準ストリップの画素値の範囲と基準マスク画像の画素値の範囲とを小さくするように基準ストリップの画素値の範囲と基準マスク画像の画素値の範囲とを調整し、且つ、対象ストリップの画素値の範囲と対象マスク画像の画素値の範囲とを小さくするように、対象ストリップの画素値の範囲と対象マスク画像の画素値の範囲とを相対的に調整するように構成される。

　範囲調整部２３１２は、ストリップ及びマスク画像の一方の画素値の範囲を他方の画素値の範囲に一致させるように、ストリップの画素値の範囲とマスク画像の画素値の範囲とを調整するように構成されてよい。例えば、範囲調整部２３１２は、基準ストリップの画素値の範囲と基準マスク画像の画素値の範囲とを一致させるように基準ストリップの画素値の範囲と基準マスク画像の画素値の範囲とを調整し、且つ、対象ストリップの画素値の範囲と対象マスク画像の画素値の範囲とを一致させるように、対象ストリップの画素値の範囲と対象マスク画像の画素値の範囲とを相対的に調整するように構成されてよい。

　例えば、範囲調整部２３１２は、マスク画像の画素値の範囲に応じてストリップの画素値の範囲を正規化するように構成される。この正規化（規格化）の幾つかの例を以下に説明するが、これらに限定されるものではない。

　正規化の第１の例を説明する。マスク画像の画素値の範囲が閉区間［０，１］に含まれる場合において、範囲調整部２３１２は、ストリップの各画素の値を、このストリップにおける最大画素値で除算する。本例では、範囲調整部２３１２は、まず、ストリップの全ての画素の値を比較して最大値（最大画素値）を特定し、このストリップの各画素の値を最大画素値で除算する。これにより、ストリップの画素値の範囲が、マスク画像の画素値の範囲と同じ閉区間［０，１］に一致される。

　正規化の第２の例を説明する。マスク画像の画素値の範囲が閉区間［０，１］に含まれる場合において、範囲調整部２３１２は、ストリップの各画素の値を、このストリップの画素値の範囲の最大値で除算する。ストリップの画素値の範囲は予め設定されており、範囲調整部２３１２は、この範囲の上限値（最大値）でストリップの各画素の値を除算する。本例によっても、ストリップの画素値の範囲が、マスク画像の画素値の範囲と同じ閉区間［０，１］に一致される。

　範囲調整部２３１２が実行する処理により、ストリップの画素値の絶対値とマスクの画素値の絶対値との差を小さくすることができ、ストリップとマスクとを用いた相関演算においてマスクの効果が無視されることが無くなり、相関係数の算出を正確に行うことが可能となる。特に、単精度浮動小数点型（フロート型）の演算を採用した場合であっても、有効桁数の少なさに起因する丸め誤差の影響を排除する（低減する）ことができる。

　本開示ではストリップの画素値の範囲のみを変更する例について主に説明するが、マスク画像の画素値の範囲のみを変更してもよいし、ストリップの画素値の範囲とマスク画像の画素値の範囲との双方を変更してもよい。

＜合成画像生成部２３１３＞
　範囲調整部２３１２による画素値範囲調整が適用されたストリップ及びマスク画像について、合成画像生成部２３１３は、基準ストリップ及び対象ストリップのそれぞれにマスク画像を合成して２つの合成画像を生成する。典型的には、合成画像生成部２３１３は、基準ストリップに基準マスク画像を合成して基準合成画像を生成し、且つ、対象ストリップに対象マスク画像を合成して対象合成画像を生成する。

　ストリップとマスク画像とを合成する処理は、非特許文献１と同じ要領で実行されてよい。ただし、非特許文献１の手法とは異なり、本実施形態では、ストリップの画素値の範囲とマスク画像の画素値の範囲との間の差を小さくするように調整を行っている。典型的には、基準マスク画像及び対象マスク画像の画素値の範囲に一致するように、基準ストリップ及び対象ストリップの画素値の範囲が正規化されている。

　例えば、非特許文献１の式（１９）と同じ要領で、合成画像生成部２３１３は、画素値範囲が正規化されたストリップを、マスク画像と同じサイズ及び同じ形状の画像に埋め込むように構成されてよい。

　更に、合成画像生成部２３１３は、ストリップの埋め込み画像とマスク画像との合成画像を生成する。この合成画像は、非特許文献１における「ｆ´（ｒ）ｍ_ｆ（ｒ）」（第１８００頁の第２行目など）に相当するが、前述したように、ストリップの埋め込み画像の値が式（１９）のそれとは異なっている。

＜相互相関関数算出部２３１４＞
　相互相関関数算出部２３１４は、合成画像生成部２３１３により生成された２つの合成画像に基づいて複数の相互相関関数を求める。相互相関関数の算出は、非特許文献１と同じ要領で実行される。例えば、相互相関関数算出部２３１４は、基準ストリップ及び基準マスク画像から生成された基準合成画像と、対象ストリップ及び対象マスク画像から生成された対象合成画像とに基づいて、非特許文献１の式（３３）に含まれる６個の相互相関関数（画像相互相関、image cross-correlation）を算出する。

＜相関係数算出部２３１５＞
　相関係数算出部２３１５は、相互相関関数算出部２３１４により算出された複数の相互相関関数に基づいて相関係数を算出する。この演算は、非特許文献１の式（３３）に従う。

＜ｘｙシフト量算出部２３１６＞
　ｘｙシフト量算出部２３１６は、相関係数算出部２３１５により算出された相関係数に基づいて、基準ストリップと対象ストリップとの間のｘｙ方向（ラテラル方向、横方向）のシフト量を算出する。

　ｘｙ方向シフト量を算出するための演算は、例えば、非特許文献１の「rough lateral motion correction」（第１７８７頁）に相当する演算を含む。この場合、ｘｙシフト量算出部２３１６は、相互相関関数の最大値を求めることによってｘｙ方向シフト量を推定するように構成される。

　更に、ｘｙシフト量算出部２３１６は、スロードリフト（slow drift）やトレモア（tremor）などの眼球運動に起因するラテラル方向の小さなシフト量を算出するために、非特許文献１の「fine lateral motion correction」（第１７８９頁）に相当する演算を実行するように構成されてもよい。

　ｘｙシフト量算出部２３１６により求められたｘｙ方向シフト量データは、所定の情報（付帯情報）とともにシフト量記録部２３３に送られて記録される。この付帯情報の例として、ｘｙ方向シフト量データの識別子、ｘｙ方向シフト量に対応する相関係数などがある。

　ｘｙ方向シフト量データの識別子は、例えば、次のいずれかを含んでいてよい：光走査（リサジュースキャン）の時間軸により表現される時間（例えば図６を参照）；サイクルの順序付けに基づく番号（サイクル番号）；スキャンの順序に基づく番号（スキャン番号）；ストリップの順序付けに基づく番号（ストリップ番号）；サブボリュームの順序付けに基づく番号（サブボリューム番号）；ｘｙ座標系で表現される位置（画素位置、走査位置など）。これらの情報は全て同等である。つまり、これらの情報のうちの任意の２つの情報について、一方の情報を他方の情報に変換することが可能である。

　ｘｙ方向シフト量データの識別子は、複数のｘｙ方向シフト量データに共通の識別子（共通識別子）を含んでいてもよい。例えば、画像選択部２２５により指定された複数の初期基準ストリップのそれぞれについて、その初期基準ストリップに基づき実行された再帰的レジストレーションにおいて逐次に得られた複数のｘｙ方向シフト量データに対して共通識別子を付与することができる。この共通識別子により、各々の初期基準ストリップに対応する複数のｘｙ方向シフト量データ（つまり、各々の再帰的レジストレーションに対応する複数のｘｙ方向シフト量データ）を識別して処理を行うことが可能になる。共通識別子は、各ｘｙ方向シフト量データの識別子（個別識別子）とは別の情報として設けられてもよいし、共通識別子と個別識別子とが一体の情報として形成されていてもよい。

＜レジストレーション部２３１７＞
　レジストレーション部２３１７は、ｘｙシフト量算出部２３１６により算出されたラテラル方向のシフト量（ｘｙ方向シフト量）に基づいて、ラテラル方向のレジストレーションを実行する。例えば、このレジストレーションは、非特許文献１の「rough lateral motion correction」（第１７８７頁）に相当する処理を含む。レジストレーション部２３１７は、ｘｙシフト量算出部２３１６により算出されたラテラル方向のシフト量を打ち消すように、基準ストリップと対象ストリップとの間のレジストレーションを実行することができる。

　ｘｙシフト量算出部２３１６が非特許文献１の「fine lateral motion correction」（第１７８９頁）に相当する演算を実行した場合、レジストレーション部２３１７は、非特許文献１の「fine lateral motion correction」（第１７８９頁）に相当するレジストレーションを実行することによって、基準ストリップと対象ストリップとの間のラテラル方向の小さなモーションアーティファクトを除去することができる。

＜マージ処理部２３１８＞
　マージ処理部２３１８は、レジストレーション部２３１７により相対位置調整がなされた基準ストリップと対象ストリップとのマージ画像を構築する。この処理も非特許文献１に記載された手法と同じ要領で実行される。

　前述したように、画像データ処理部２３１は、画像データ構築部２２０により構築された複数のストリップを大きさに応じた順序付けに従って、上記した一連の処理を逐次に実行する。これにより、画像データ構築部２２０により構築された複数のストリップから、ラテラル方向のモーションアーティファクトが補正されたマージ画像が得られる。このマージ画像は、典型的には、リサジュースキャンが適用された範囲全体を表現した画像である。

　なお、前述したように、複数のストリップは、リサジュースキャンで収集されたボリューム（３次元データ）を分割して得られた複数のサブボリュームに基づく複数の正面プロジェクション画像である。したがって、画像データ処理部２３１により複数のストリップから構築されるマージ画像は、ラテラル方向の位置調整がなされた複数のサブボリューム（及びそれらのマージ画像）を提供する。ラテラル方向に直交する深さ方向のレジストレーション及びマージ処理は、本例では、ｚシフト量算出部２３２及び画像データ補正部２３５によって実行される（後述）。

＜ｚシフト量算出部２３２＞
　ｚシフト量算出部２３２は、リサジュースキャンで収集されたボリュームを分割して得られた複数のサブボリューム（つまり、複数のストリップの元になった複数のサブボリューム）の間のｚ方向（Ａライン方向、光軸方向、深さ方向）のシフト量を算出する。

　本例は、画像データ処理部２３１によってｘｙ方向（ラテラル方向、横方向）の位置調整を実行した後に、ｚシフト量算出部２３２及び画像データ補正部２３５によってｚ方向の位置調整を実行するように構成されている。すなわち、本例のｚシフト量算出部２３２は、画像データ処理部２３１によって得られたデータ（例えば、ｘｙ方向シフト量データ、マージ画像）を利用してｚ方向のシフト量を算出するように構成されている。

　ｚシフト量算出部２３２が実行する処理の例を説明する。ｘｙ方向シフト量算出において考慮されるストリップのペア（基準ストリップ及び対象ストリップ）と同様に、ｚ方向シフト量算出においてもサブボリュームのペアが考慮される。サブボリュームのペアは、ｘｙ方向シフト量算出において考慮された基準ストリップ及び対象ストリップに対応する２つのサブボリュームであってよく、これらサブボリュームをそれぞれ基準サブボリューム及び対象サブボリュームと呼ぶ。

　対応する基準ストリップ及び対象ストリップの間のｘｙ方向のレジストレーションの結果に基づいて、基準サブボリューム及び対象サブボリュームは、ｘｙ方向の位置調整がなされている。また、前述したように、リサジュースキャンにおける任意の２つのストリップは４箇所で重複（交差）しており、したがって、任意の２つのサブボリュームも４箇所で重複（交差）している。

　ｚシフト量算出部２３２は、まず、ｘｙ方向の位置調整がなされた基準サブボリュームと対象サブボリュームとの交差領域（共通領域）を特定する。各交差領域は、３次元画像データである。

　次に、ｚシフト量算出部２３２は、特定された交差領域に断面を設定する。この断面は、例えば、ｘｙ平面における任意の軸（例えば、ｘ軸、ｙ軸、又は、ｘ軸及びｙ軸の双方に斜交する軸）とｚ軸とによって貼られた平面である。なお、断面は平面に限定されず、曲面などであってもよい。

　次に、ｚシフト量算出部２３２は、設定された断面の画像を基準サブボリュームから構築し、且つ、同断面の画像を対象サブボリュームから構築する。基準サブボリュームから構築された断面画像を基準断面画像と呼び、対象サブボリュームから構築された断面画像を対象断面画像と呼ぶ。基準断面画像と対象断面画像とは、ｘｙ方向のレジストレーションがなされた基準サブボリュームと対象サブボリュームとの交差領域における同じ断面を表す画像である。

　次に、ｚシフト量算出部２３２は、基準断面画像を解析して被検眼Ｅの所定部位の像を特定し、且つ、対象断面画像を解析して同部位の像を特定する。この部位は、任意の部位であってよく、例えば眼底Ｅｆの表面（網膜表面、内境界膜、網膜と硝子体との境界）であってよい。また、被検眼Ｅに人工物が移植されている場合、この人工物の像を特定してもよい。所定部位の像を特定するための解析は、例えば、セグメンテーションを含んでいてよい。

　次に、ｚシフト量算出部２３２は、基準断面画像から特定された所定部位の像（基準像）のｚ座標を求め、且つ、対象断面画像から特定された所定部位の像（対象像）のｚ座標を求める。典型的には、所定部位の像は複数の画素からなり、これらの画素のｚ座標は一定ではない。例えば、一般に、網膜表面の大局的な形状は＋ｚ方向に凸な湾曲形状であり、その局所的箇所の像（基準像及び対象像の例）はｚ軸に対して少なくとも部分的に傾斜している。ｚシフト量算出部２３２は、所定部位の像を構成する画素群の少なくとも１つの画素に基づいて、この像のｚ座標を求めることができる。例えば、この画素群からｚ座標の統計量を求め、この統計量を当該像のｚ座標に採用することができる。この統計量は、任意の代表値（要約統計量）であってよく、その例として、最大値、最小値、平均値、中央値、最頻値、分位点などがある。一方、所定部位の像が単一の画素からなる場合には、この画素のｚ座標を所定部位のｚ座標として採用することができる。

　次に、ｚシフト量算出部２３２は、基準像のｚ座標と対象像のｚ座標との差（ｚ方向のシフト量）を求める。幾つかの例示的な態様では、ｚ方向シフト量が閾値を超える場合に画像データ補正部２３５による補正処理を実行するように処理制御を行うことができる。この閾値は、例えば、ゼロ又は正値であり、予め設定された固定値又は可変値であってよい。

　上記した例では、３次元画像データ（サブボリューム）から２次元画像データ（断面画像）を構築してｚ方向シフト量を求めているが、幾つかの例示的な態様では、２次元画像データの構築を行うことなく３次元画像データからｚ方向シフト量を求めてもよい。この場合、基準断面画像及び対象断面画像のペアに対して実行されたものと同様の処理を基準サブボリューム及び対象サブボリュームのペアに対して実行することができる。

　上記した例では、３次元画像データ（サブボリューム）から２次元画像データ（断面画像）を構築し、この２次元画像データから所定部位の像を特定してｚ方向シフト量を求めているが、幾つかの例示的な態様では、３次元画像データから所定部位の像を特定し、この所定部位の像の２次元画像データを構築してｚ方向シフト量を求めてもよい。

　また、ｚシフト量算出部２３２は、基準サブボリューム及び対象サブボリュームのペアに基づき算出される相関係数に基づいて、基準サブボリュームと対象サブボリュームとの間のｚ方向シフト量を求めるように構成されてもよい。サブボリューム間の相関係数を求める方法は、例えば、相関係数算出部２３１５により実行される相関係数演算方法と同様であってよい。

　ｚ方向シフト量を算出するための方法は以上に例示した方法には限定されず、基準サブボリューム及び対象サブボリュームのペアに対して適用可能な任意の方法であってよく、又は、サブボリュームと異なるオブジェクトのペアに対して適用可能な任意の方法であってもよい。サブボリュームと異なるオブジェクトの例として、サブボリュームのレンダリング画像（ｘ方向又はｙ方向にプロジェクションした画像など、ｚ方向の次元を有する画像）がある。

　幾つかの例示的な態様では、ｚシフト量算出部２３２により求められたｚ方向シフト量データは、所定の情報（付帯情報）とともにシフト量記録部２３３に送られて記録される。この付帯情報の例として、ｚ方向シフト量データの識別子、ｚ方向シフト量を算出するための演算の途中で求められた値や量（例えば、相関係数）などがある。

　ｘｙ方向シフト量データの個別識別子と同様に、ｚ方向シフト量データの個別識別子は、次のいずれかを含んでいてよい：光走査（リサジュースキャン）の時間軸により表現される時間（例えば図６を参照）；サイクルの順序付けに基づく番号（サイクル番号）；スキャンの順序に基づく番号（スキャン番号）；ストリップの順序付けに基づく番号（ストリップ番号）；サブボリュームの順序付けに基づく番号（サブボリューム番号）；ｘｙ座標系で表現される位置（画素位置、走査位置など）。また、ｘｙ方向シフト量データの識別子と同様に、ｚ方向シフト量データの識別子は、複数のｚ方向シフト量データに共通の識別子（共通識別子）を含んでいてもよい。

＜シフト量記録部２３３＞
　シフト量記録部２３３は、ｘｙシフト量算出部２３１６により求められたｘｙ方向シフト量データ及びその付帯情報を記録する。本実施形態の付帯情報は、前述したｘｙ方向シフト量データの識別子（例えば、個別識別子及び共通識別子）を含む。シフト量記録部２３３は、例えば、記憶装置と、この記憶装置へのデータ書き込み及びこの記憶装置からのデータ読み出しを行うプロセッサとを含む。

　前述したように、本実施形態のｘｙシフト量算出部２３１６は、基準ストリップと対象ストリップとの各ペアについてｘｙ方向シフト量データを生成する。幾つかの例示的な態様では、シフト量記録部２３３は、基準ストリップと対象ストリップとの全てのペアに対応するｘｙ方向シフト量データを収集し記録する。また、幾つかの例示的な態様では、シフト量記録部２３３は、基準ストリップと対象ストリップと全てのペアのうちの一部のペアに対応するｘｙ方向シフト量データを収集し記録する。

　幾つかの例示的な態様において、シフト量記録部２３３は、ｚシフト量算出部２３２により求められたｚ方向シフト量データ及びその付帯情報を記録する。本実施形態の付帯情報は、前述したｚ方向シフト量データの識別子（例えば、個別識別子及び共通識別子）を含む。シフト量記録部２３３がｚ方向シフト量データを記録する構成が適用される場合、図４Ｃの図面におけるｚシフト量算出部２３２とシフト量記録部２３３との間に、ｚシフト量算出部２３２からシフト量記録部２３３データへのデータ送信を示す矢印が付加される。

　幾つかの例示的な態様では、シフト量記録部２３３は、基準サブボリュームと対象サブボリュームとの全てのペアに対応するｚ方向シフト量データを収集し記録する。また、幾つかの例示的な態様では、シフト量記録部２３３は、基準サブボリュームと対象サブボリュームと全てのペアのうちの一部のペアに対応するｚ方向シフト量データを収集し記録する。

　なお、ストリップとサブボリュームとの間には対応関係があり、したがって、基準ストリップと対象ストリップとのペアと基準サブボリュームと対象サブボリュームとの間には対応関係がある。よって、シフト量記録部２３３がｘｙ方向シフト量データ及びｚ方向シフト量データの双方を記録する態様では、ｘｙ方向シフト量データとｚ方向シフト量データとを対応付けることができる。

＜調整量算出部２３４＞
　調整量算出部２３４は、シフト量記録部２３３により記録されたデータに基づいて画像セットの相対位置調整のための調整量を算出するように構成されている。

　幾つかの例示的な態様の調整量算出部２３４は、シフト量記録部２３３により記録されたｘｙ方向シフト量データに基づいて複数のストリップの相対位置調整のための調整量を算出する。シフト量記録部２３３により記録されたｘｙ方向シフト量データには、前述した個別識別子及び共通識別子が付されているものとする。各共通識別子は、１つの初期基準ストリップ（つまり、１つの再帰的レジストレーション）から得られたｘｙ方向シフト量データの集合（ｘｙ方向シフト量データ集合）を表し、各個別識別子は、その集合における各ｘｙ方向シフト量データの位置（例えば、時間的位置又は空間的位置）を表す。

　１つのｘｙ方向シフト量データ集合は、画像セットに含まれる複数のストリップの間の相対位置を表す情報（相対位置情報）であり、従来の手法（例えば、引用文献１の手法）で得られる相対位置情報と実質的に同じものである。しかし、本実施形態は、複数の初期基準画像を用いて複数の再帰的レジストレーションを実行することで複数の相対位置情報を取得する点において、従来の手法と異なっている。これに加えて、本実施形態は、調整量算出部２３４により、このような複数の相対位置情報に基づいて、画像セットに含まれる複数のストリップの相対位置調整のための調整量を算出する点において、従来の手法と異なっている。

　画像セットに含まれる複数の画像の相対位置調整のための調整量を算出する方法は任意であってよい。その幾つかの非限定的な例を以下に説明する。これらの例では、ラテラル方向（ｘｙ方向）の相対位置調整のための調整量を算出する方法について説明を行うが、アキシャル方向（ｚ方向）の相対位置調整のための調整量についても同様の方法で算出できることは、当業者であれば理解できるであろう。

　ラテラル方向の相対位置調整のための調整量を算出する処理の例について図９Ａ～図９Ｄ及び図１０を参照しつつ説明する。本例では、図６の初期基準画像選択方法（光走査の時間軸に基づく初期基準画像選択）が適用されているが、他の初期基準画像選択方法が適用される場合においても同様の処理を実行することができる。

　図９Ａの上段の時間軸ｔ及び下段の時間軸はいずれも、本例で実行される光走査（リサジュースキャン）の時間軸である。符号Ｆは、リサジュースキャンを用いて取得された画像セット（複数のストリップ）を示す。本例では、リサジュースキャンの時間軸ｔに対して、時系列順に、３個の区間Ｔ１、Ｔ２、及びＴ３が設定されている。区間Ｔ１に対応する画像群（ストリップ群）から画像選択部２２５により選択された画像（ストリップ）を符号Ｒ１で表し、区間Ｔ２に対応する画像群（ストリップ群）から画像選択部２２５により選択された画像（ストリップ）を符号Ｒ２で表し、区間Ｔ３に対応する画像群（ストリップ群）から画像選択部２２５により選択された画像（ストリップ）を符号Ｒ３で表す。３つのストリップＲ１、Ｒ２、及びＲ３は、それぞれ、３回実行される再帰的レジストレーションの初期基準画像（初期基準ストリップ）として用いられる。符号ｔ１は第１の初期基準ストリップＲ１に対応する時間を示し、符号ｔ２は第２の初期基準ストリップＲ２に対応する時間を示し、符号ｔ３は第３の初期基準ストリップＲ３に対応する時間を示す。

　図９Ａの下段の座標系は、リサジュースキャンの時間軸ｔと、ラテラル方向の位置ずれ量（シフト量）を示す座標軸ΔＸとにより張られている。図９Ａのラテラル方向シフト量ΔＸは、ベクトル量でもよいし、スカラー量でもよい。ベクトル量ΔＸは、例えば、ｘ方向のシフト量とｙ方向のシフト量とを含む：ΔＸ＝（Δｘ、Δｙ）。スカラー量ΔＸは、例えば、ｘ方向のシフト量（ΔＸ＝Δｘ）、ｙ方向のシフト量（ΔＸ＝Δｙ）、又は、ｘｙ平面上の任意方向におけるシフト量であってよい。以下、ΔＸはベクトル量ΔＸ＝（Δｘ、Δｙ）である場合について説明するが、ΔＸが他のベクトル量又はスカラー量である場合においても同様の処理を行うことができる。

　図９Ａには、このような座標系（ｔ、ΔＸ）で表現された３つのグラフＧ１、Ｇ２、及びＧ３が示されている。第１のグラフＧ１は、第１の初期基準ストリップＲ１を起点として実行された第１の再帰的レジストレーションで求められた各ストリップのｘｙ方向シフト量（ｘｙ方向シフト量データ集合、相対位置情報）を表している。第２のグラフＧ２は、第２の初期基準ストリップＲ２を起点として実行された第２の再帰的レジストレーションで求められた各ストリップのｘｙ方向シフト量（ｘｙ方向シフト量データ集合、相対位置情報）を表している。第３のグラフＧ３は、第３の初期基準ストリップＲ３を起点として実行された第３の再帰的レジストレーションで求められた各ストリップのｘｙ方向シフト量（ｘｙ方向シフト量データ集合、相対位置情報）を表している。

　本例では、以上に説明したデータがシフト量記録部２３３から調整量算出部２３４に入力される。本例の調整量算出部２３４は、まず、３つのグラフＧ１～Ｇ３がそれぞれ示すｘｙ方向シフト量を比較可能にするために、３つのグラフＧ１～Ｇ３のｘｙ方向シフト量ΔＸを整列する処理を実行する。

　例えば、調整量算出部２３４は、次式によって３つのグラフＧ１～Ｇ３のｘｙ方向シフト量ΔＸを整列する：ΔＸ´_ｉ ^{（Ｒ（ｎ＋１））}＝ΔＸ_ｉ ^{（Ｒ（ｎ＋１））}＋ΔＸ´_{Ｒ（ｎ＋１）} ^（Ｒｎ）。ここで：ΔＸは、整列前のｘｙ方向シフト量ベクトルΔＸ＝（Δｘ、Δｙ）を示す；ΔＸ´は、整列後のｘｙ方向シフト量ベクトルΔＸ´＝（Δｘ´、Δｙ´）を示す；ｉは、前述した個別識別子に基づくストリップのインデックスを示す；Ｒｎは、リサジュースキャンの時間軸ｔに設定された区間Ｔｎ（ｎ＝１、２、３）に対応するストリップ群から画像選択部２２５により選択されたストリップを示す。

　図９Ｂは、第１のグラフＧ１と第２のグラフＧ２とを整列する処理の態様を表している。この整列処理では、調整量算出部２３４は、まず、時間ｔ２における第１のグラフＧ１のｘｙ方向シフト量に対する、時間ｔ２における第２のグラフＧ２のｘｙ方向シフト量の差分（差分ベクトル、ベクトル差）を求める。更に、調整量算出部２３４は、算出された差分だけ、第２のグラフＧ２をΔＸ軸方向に移動する。この整列処理は、時間ｔ２における第１のグラフＧ１と第２のグラフＧ２との差分をキャンセルするように第２のグラフＧ２をΔＸ軸方向に平行移動するものである。第１のグラフＧ１に対して整列された第２のグラフＧ２を符号Ｇ２´で示す。

　同様に、図９Ｃは、第１のグラフＧ１に対して整列された第２のグラフＧ２´と、第３のグラフＧ３とを整列する処理の態様を表している。この整列処理では、調整量算出部２３４は、まず、時間ｔ３における第２のグラフＧ２´のｘｙ方向シフト量に対する、時間ｔ３における第３のグラフＧ３のｘｙ方向シフト量の差分（差分ベクトル、ベクトル差）を求める。更に、調整量算出部２３４は、算出された差分だけ、第３のグラフＧ３をΔＸ軸方向に移動する。この整列処理は、時間ｔ３における第２のグラフＧ２´と第３のグラフＧ３との差分をキャンセルするように第３のグラフＧ３をΔＸ軸方向に平行移動するものである。第２のグラフＧ２´に対して整列された第３のグラフＧ３を符号Ｇ３´で示す。

　これにより、図９Ｄに示す３つのグラフＧ１、Ｇ２´、及びＧ３´が得られる。なお、複数のグラフを比較可能にするための整列処理は上記の例に限定されず、任意であってよい。例えば、複数のグラフのうちの１つのグラフに対して他の全てのグラフを合わせるようにしてもよい。また、整列処理を行うことなく複数のグラフの比較を行うようにアルゴリズムを構成してもよい。

　更に、本例の調整量算出部２３４は、任意の第ｉ番目のストリップについて、この第ｉ番目のストリップのｘｙ方向シフト量と、他の第ｉ番目のストリップのｘｙ方向シフト量との間の幾何学的距離の総和が最小となる初期基準ストリップを最適な初期基準ストリップとして選択し、この最適な初期基準ストリップに基づく再帰的レジストレーションで得られたｘｙ方向シフト量を第ｉ番目のストリップのｘｙ方向シフト量に指定するように構成されていてよい。

　換言すると、本例の調整量算出部２３４は、第ｉ番目のストリップについて算出された３つのｘｙ方向シフト量ベクトルのうちから、他の２つのｘｙ方向シフト量ベクトルとの間の幾何学的距離の和が最小となるようなｘｙ方向シフト量ベクトルを選択し、選択されたｘｙ方向シフト量ベクトルをこの第ｉ番目のストリップの最適なｘｙ方向シフト量ベクトルとして採用するものであってよい。

　このような演算は、例えば次式のように表現される（図１０を参照）：ΔＸ´_ｉ＝ΔＸ´_ｉ ^{（Ｒ´ｉ）}；Ｒ´_ｉ＝ａｒｇｍｉｎ_Ｒｊ［Σ_{Ｒｊ≠Ｒｋ}｛ΔＸ´_ｉ ^（Ｒｋ）－ΔＸ´_ｉ ^（Ｒｊ）｝^２］。ここで、和Σは、Ｒｊ以外のＲｋ（Ｒｊ≠Ｒｋ）にわたって演算され、且つ、最小点集合ａｒｇｍｉｎ（ａｒｇｕｍｅｎｔ　ｏｆ　ｔｈｅ　ｍｉｎｉｍｕｍ）は、Ｒｊにわたって演算される。また、ΔＸ´_ｉ ^（Ｒｋ）－ΔＸ´_ｉ ^（Ｒｊ）は、２つの第ｉ番目のストリップの間におけるｘｙ方向の位置ずれ量を示すベクトルに相当する。

　ここではラテラル方向（ｘｙ方向）の調整量の算出について説明したが、アキシャル方向（ｚ方向）の調整量についても同じ要領で算出することができる。また、ここでは各々のストリップについて独立に最適シフト量を求めているが、時間的又は空間的に近くに位置するストリップの情報を用いて最適シフト量を決定することも可能である。例えば、第ｉ番目のストリップについての演算において、第ｉ－１番目のストリップに関する情報及び第ｉ＋１番目のストリップに関する情報を利用することが可能である。

　このように、本例の処理は、複数の初期基準画像を用いて複数回の再帰的レジストレーションを実行し、これらの再帰的レジストレーションで得られた複数のシフト量データから各ストリップの最適なシフト量を選択する。これは、単一の初期基準画像を用いて１回の再帰的レジストレーションを実行する従来の手法とは異なるものである。更に、本例における整列処理や最適初期基準画像選択処理は、従来の手法には無い、本実施形態に特有の新規な技術である。

＜画像データ補正部２３５＞
　画像データ補正部２３５は、調整量算出部２３４により求められた調整量に少なくとも基づいて、画像セットに含まれる複数の画像の相対位置調整を実行するように構成されている。

　例えば、図４Ｃに示す態様の画像データ補正部２３５は、画像セットに含まれる複数のストリップのラテラル方向（ｘｙ方向）における相対位置調整を、調整量算出部２３４により求められた各ストリップのラテラル方向の調整量（ｘｙ方向シフト量）に基づき実行するとともに、これら複数のストリップに対応する複数のサブボリュームのアキシャル方向（ｚ方向）における相対位置調整を、ｚシフト量算出部２３２により求められた各サブボリュームのｚ方向シフト量に基づき実行する。

　なお、調整量算出部２３４がアキシャル方向の調整量を求めるように構成された態様においては、画像データ補正部２３５は、画像セットに含まれる複数のストリップのラテラル方向における相対位置調整を、調整量算出部２３４により求められた各ストリップのラテラル方向の調整量に基づき実行するとともに、これら複数のストリップに対応する複数のサブボリュームのアキシャル方向における相対位置調整を、調整量算出部２３４により求められた各サブボリュームのアキシャル方向の調整量に基づき実行することができる。

　ｚシフト量算出部２３２により求められたｚ方向シフト量に基づく複数のサブボリュームのアキシャル方向の相対位置調整について説明する。ｘｙ方向シフト量に基づくレジストレーションを実行するレジストレーション部２３１７と同様に、画像データ補正部２３５は、ｚシフト量算出部２３２により取得された対応するｚ方向シフト量を打ち消すように、基準サブボリュームと対象サブボリュームとの間のレジストレーションを実行する。

　幾つかの例示的な態様では、画像データ補正部２３５は、対応するｚ方向シフト量が閾値を超える場合にのみ、基準サブボリュームと対象サブボリュームとの間のレジストレーションを実行するように制御されてよい。

　画像データ補正部２３５による位置調整の結果として得られる画像セットは、ｘ方向、ｙ方向、及びｚ方向の全てにおいてレジストレーション（相対位置調整）がなされたサブボリューム群、つまり、３次元的レジストレーションがなされたサブボリューム群である。このようなレジストレーションを適用する前のサブボリューム群は、リサジュースキャンで収集されたボリュームである。データ処理部２３０は、３次元的レジストレーションにより得られたサブボリューム群を再合成することによって、３次元的レジストレーションが適用されたボリュームが得られる。

　３次元的レジストレーションが適用されたデータセットを、任意の処理（画像処理、後処理、解析、評価、可視化、レンダリングなど）に利用することができる。これにより、従来よりも高い正確度及び高い精度でモーションアーティファクトが補正された３次元データセットを用いて所望の処理を実行することが可能になる。

　なお、ｚシフト量算出部２３２及び画像データ補正部２３５は、ｚ方向のモーションアーティファクトを除去するために、非特許文献１の「axial motion correction (rough axial motion correction 及び／又は fine axial motion correction)」（第１７８９頁～第１７９０頁）に相当するレジストレーションを実行するように構成されてもよい。

＜ユーザーインターフェイス２４０＞
　ユーザーインターフェイス２４０は表示部２４１と操作部２４２とを含む。表示部２４１は表示装置３を含む。操作部２４２は各種の操作デバイスや入力デバイスを含む。ユーザーインターフェイス２４０は、例えばタッチパネルのような表示機能と操作機能とが一体となったデバイスを含んでいてもよい。ユーザーインターフェイス２４０の少なくとも一部を含まない実施形態を構築することも可能である。例えば、表示デバイスは、眼科検査装置１に接続された外部装置であってよい。

＜データ受付部２５０＞
　データ受付部２５０は、外部装置からデータを取得する。この外部装置は、例えば、コンピュータ、記憶装置、記録媒体、情報システム（例えば、病院情報システム、電子カルテシステム、画像アーカイビングシステム）などであってよい。また、外部装置は、例えば、眼科検査装置１に直接に接続された装置、ローカルエリアネットワーク（ＬＡＮ）を介して眼科検査装置１に接続された装置、ワイドエリアネットワーク（ＷＡＮ）を介して眼科検査装置１に接続された装置などであってよい。

　データ受付部２５０は、例えば、通信インターフェイス、ドライブ装置などを含んでいてもよい。また、データ受付部２５０は、紙葉類に記録されたデータを読み取るためのスキャナを含んでいてもよい。また、データ受付部２５０は、ユーザーがデータを入力するための入力デバイスを含んでいてもよい。入力デバイスは、キーボード、ペンタブレットなどであってよい。

＜動作＞
　眼科検査装置１の動作の幾つかの例を説明する。

＜第１の動作例＞
　第１の動作例では図１１に示す動作について説明する。本動作例は、冗長的データ収集に基づき生成された画像セットから複数の初期基準画像を選択し、これらの初期基準画像を用いて複数のレジストレーションによって画像セットの相対位置調整を実行するものである。本動作例において、冗長的データ収集はリサジュースキャンであり、画像セットは一連のストリップ（複数のストリップ）であり、初期基準画像は初期基準ストリップである。

　患者ＩＤの入力、走査モードの設定（リサジュースキャンの指定）、固視標の提示、アライメント、フォーカス調整、ＯＣＴ光路長調整など、従来と同様の準備動作が、ステップＳ１の前に実行される。

（Ｓ１：リサジュースキャンによりボリュームを収集する）
　所定の走査開始トリガー信号を受けて、走査制御部２１１１は、被検眼Ｅ（眼底Ｅｆ）に対するＯＣＴスキャン（リサジュースキャン）の適用を開始する。走査開始トリガー信号は、例えば、所定の準備動作（アライメント、フォーカス調整、ＯＣＴ光路長調整など）が完了したことに対応して、又は、操作部２４２を用いて走査開始指示操作が行われたことに対応して生成される。走査制御部２１１１は、走査プロトコル２１２１（リサジュースキャンに対応したプロトコル）に基づいて光スキャナ４４及びＯＣＴユニット１００などを制御することによってリサジュースキャンを眼底Ｅｆに適用する。リサジュースキャンにより収集されたデータセット（ボリューム）は画像データ構築部２２０に送られる。

（Ｓ２：一連のサブボリュームに分割する）
　画像データ構築部２２０は、ステップＳ１で収集されたボリュームを一連のサブボリュームに分割する。各サブボリュームは、被検眼Ｅの動きの影響を受けていないとみなされる。

（Ｓ３：一連のストリップを構築する）
　画像データ構築部２２０は、ステップＳ２で得られた各サブボリュームの正面プロジェクション画像（ストリップ）を構築する。これにより、ステップＳ２で生成された一連のサブボリュームから一連のストリップ（画像セット）が構築される。一連のストリップは、データ処理部２３０に送られる。

（Ｓ４：複数の初期基準ストリップを選択する）
　画像選択部２２５は、ステップＳ３で構築された一連のストリップから、複数の初期基準ストリップを選択する。この処理は、前述したいずれかの初期基準画像選択方法を用いて実行される。

（Ｓ５：複数の再帰的レジストレーションを実行する）
　画像データ処理部２３１は、ステップＳ４で選択された複数の初期基準ストリップをそれぞれ基準とした複数の再帰的レジストレーションを一連のストリップに適用する。１つの初期基準ストリップに基づく再帰的レジストレーションにより、一連のストリップについて１つのｘｙ方向シフト量データ集合（相対位置情報）が生成される。複数の初期基準ストリップのそれぞれに基づき再帰的レジストレーションを実行することにより、一連のストリップについて、複数のｘｙ方向シフト量データ集合（複数の相対位置情報）が生成される。

（Ｓ６：一連のストリップの相対位置調整量を算出する）
　画像データ処理部２３１は、ステップＳ５で生成された複数の相対位置情報に基づいて、一連のストリップの相対位置調整（モーションアーティファクト補正、特にｘｙ方向のモーションアーティファクト補正）のための調整量を算出する。この処理は、前述したいずれかの相対位置調整量算出方法を用いて実行される。

（Ｓ７：一連のストリップの相対位置調整を実行する）
　画像データ処理部２３１は、ステップＳ６で算出された調整量に基づいて一連のストリップの相対位置調整を実行する。本ステップでは、ラテラル方向（ｘｙ方向）における相対位置調整が実行される。本ステップにおいて、更に、非特許文献１の「fine lateral motion correction」と同様の処理を実行してもよい。

（Ｓ８：一連のサブボリュームの相対位置調整を実行する）
　画像データ処理部２３１は、ステップＳ２で生成された一連のサブボリュームの相対位置調整を実行する。この相対位置調整は、ラテラル方向（ｘｙ方向）の相対位置調整と、アキシャル方向（ｚ方向）の相対位置調整とが含まれる。ラテラル方向の相対位置調整には、ステップＳ７で実行された一連のストリップの相対位置調整の結果が利用される。また、アキシャル方向（ｚ方向）の相対位置調整は、非特許文献１の「axial motion correction (rough axial motion correction 及び／又は fine axial motion correction)」に相当する処理であってよい。

　以上で、本動作例に係る処理は終了となる（エンド）。

＜第２の動作例＞
　第２の動作例では図１２に示す動作について説明する。本動作例は、第１の動作例のステップＳ４～Ｓ７の処理の１つの例を提供するものである。本動作例のステップＳ１１の処理は第１の動作例のステップＳ３の処理に続いて実行され、且つ、本動作例のステップＳ２０の処理の後に第１の動作例のステップＳ８の処理が実行される。

　本動作例では、冗長的データ収集で得られた画像セットを時間的に分割する場合の例を説明するが、画像セットを空間的に分割する場合又は他の分割方法が適用される場合においても本動作例に係る処理と同様の処理を実行することができる。

（Ｓ１１：一連のストリップをＫ個のストリップ群に分ける）
　例えば図６を参照して説明した処理と同じ要領で、画像選択部２２５は、ステップＳ３で構築された一連のストリップをＫ個のストリップ群に分ける（Ｋは２以上の整数）。

（Ｓ１２：各ストリップ群から初期基準ストリップを選択する）
　次に、画像選択部２２５は、ステップＳ１１で生成されたＫ個のストリップ群のそれぞれから初期基準ストリップを選択する。本例では、各ストリップ群から初期基準ストリップが１つずつ選択されるものとするが、いずれかのストリップ群から２つ以上の初期基準ストリップを選択してもよいし、いずれかのストリップ群から初期基準ストリップを選択しなくてもよい。

（Ｓ１３：レジストレーションを開始する（ｋ＝１））
　画像データ処理部２３１は、ステップＳ１２で選択された複数の初期基準ストリップをそれぞれ基準とした複数のレジストレーションを開始する。

　Ｋ個の初期基準ストリップのインデックスをｋで表す（ｋ＝１、２、・・・、Ｋ）。インデックスｋは、複数の初期基準ストリップに対して任意の順序で付されている。例えば、インデックスｋは、時間軸に応じた順序で付されてもよいし（図９Ａ～図９Ｄを参照）、ストリップのサイズに応じた順序で付されてもよいし、他の情報に応じた順序で付されてもよい。

　本動作例では、Ｋ個の初期基準ストリップに基づくＫ回のレジストレーションをインデックスｋ＝１～Ｋに応じた順序で逐次に実行する場合について説明する（以下のステップＳ１４～Ｓ１６を参照）。しかし、Ｋ個の初期基準ストリップに基づくＫ回のレジストレーションの実行順序は、インデックスｋに応じた順序に限定されるものではなく、任意の順序であってよい。

　幾つかの例示的な態様では、Ｋ個の初期基準ストリップに基づくＫ回のレジストレーションのうちの少なくとも２回のレジストレーションを、少なくとも部分的に並行して実行してもよい。例えば、Ｋ回のレジストレーションのうちの第１及び第２のレジストレーションについて、第１のレジストレーションの実行時間の少なくとも一部と第２のレジストレーションの実行時間の少なくとも一部とが時間的に互いに重複するように、第１のレジストレーションと第２のレジストレーションとを並行して実行することができる。

（Ｓ１４：ｋ＝Ｋ？）
　処理対象の初期基準ストリップのインデックスｋがＫに等しくない場合（Ｓ１４：Ｎｏ）、つまり、処理対象の初期基準ストリップのインデックスｋがＫ未満である場合（ｋ＜Ｋ）、処理はステップＳ１５に移行する。一方、処理対象の初期基準ストリップのインデックスｋがＫに等しい場合（Ｓ１４：Ｙｅｓ）、処理はステップＳ１７に移行する。

（Ｓ１５：第ｋの初期基準ストリップによる再帰的レジストレーションを実行する）
　ステップ１４で「Ｎｏ」と判断された場合、画像データ処理部２３１（マスク画像生成部２３１１、範囲調整部２３１２、合成画像生成部２３１３、相互相関関数算出部２３１４、相関係数算出部２３１５、ｘｙシフト量算出部２３１６、レジストレーション部２３１７、及びマージ処理部２３１８）は、第ｋの初期基準ストリップを基準とした再帰的レジストレーションを実行する。

　第ｋの初期基準ストリップを基準とした再帰的レジストレーションにおける反復処理の各段階でｘｙシフト量算出部２３１６により算出されたｘｙ方向シフト量データは、付帯情報とともにシフト量記録部２３３に送られて記録される。この付帯情報には、例えば、第ｋの初期基準ストリップに対応する識別子（前述の共通識別子）と、そのｘｙ方向シフト量データの識別子（前述の個別識別子）とが含まれる。この個別識別子は、例えば、そのｘｙ方向シフト量データが取得された段階のレジストレーションにおける対象ストリップの識別子（又は、対象ストリップの識別子と互換な識別子）であってよい。

（Ｓ１６：ｋ＝ｋ＋１）
　第ｋの初期基準ストリップに基づく再帰的レジストレーション（ステップＳ１５）が完了したら、第ｋ＋１の初期基準ストリップを基準とした再帰的レジストレーションに移行する。

　ステップＳ１４～Ｓ１６の処理をステップＳ１４で「Ｙｅｓ」と判定されるまで繰り返し実行することによって、ステップＳ１２で選択されたＫ個の初期基準ストリップをそれぞれ基準とした複数の再帰的レジストレーションが実行され、Ｋ個のｘｙ方向シフト量データ集合（Ｋ個の相対位置情報）が得られ、シフト量記録部２３３に記録される。

　各ｘｙ方向シフト量データ集合の各要素（各ｘｙ方向シフト量データ）には、そのｘｙ方向シフト量データ集合を生成した再帰的レジストレーションにおける初期基準ストリップに対応する共通識別子と、その要素の個別識別子とが付帯されている。また、Ｋ個のｘｙ方向シフト量データ集合（Ｋ個の相対位置情報）のそれぞれには、一連のストリップにおける各ストリップのｘｙ方向シフト量データが含まれている。換言すると、シフト量記録部２３３には、ステップＳ３で構築された各ストリップについて、Ｋ個のｘｙ方向シフト量データが記録されている。

（Ｓ１７：Ｋ個のｘｙ方向シフト量データ集合の相対的位置合わせを実行する）
　次に、例えば図９Ａ～図９Ｄを参照して説明した処理と同じ要領で、調整量算出部２３４が、ステップＳ１４～Ｓ１６の複数の再帰的レジストレーションで生成されたＫ個のｘｙ方向シフト量データ集合（相対位置情報）の相対的位置合わせを実行する。

（Ｓ１８：各ストリップの最適なｘｙ方向シフト量を決定する）
　続いて、調整量算出部２３４は、ステップＳ１７の位置合わせが施されたＫ個のｘｙ方向シフト量データ集合に基づいて、一連のストリップにおける各ストリップの最適なｘｙ方向シフト量を決定する。

　例えば、一連のストリップにおける各ストリップについて、調整量算出部２３４は、まず、そのストリップに対応するｘｙ方向シフト量データをＫ個のｘｙ方向シフト量データ集合のそれぞれから抽出する。次に、調整量算出部２３４は、図１０を参照して説明した処理と同じ要領で、抽出されたＫ個のｘｙ方向シフト量データのうちから、そのストリップの最適なｘｙ方向シフト量を選択する。これにより、一連のストリップの全てについて、最適なｘｙ方向シフト量が得られる。なお、Ｋ個のｘｙ方向シフト量データから最適な１つを選択する方法はこれに限定されない。

　また、各ストリップの最適なｘｙ方向シフト量を決定する処理は、Ｋ個のｘｙ方向シフト量データから最適な１つを選択する処理に限定されない。例えば、Ｋ個のｘｙ方向シフト量データに所定の統計演算を適用して得られた統計量又はこの統計量から算出された量を、そのストリップの最適なｘｙ方向シフト量として採用してもよい。

　或いは、１つのストリップの最適なｘｙ方向シフト量を決定するために、そのストリップに対応するＫ個のｘｙ方向シフト量データに加えて又はそれらの代わりに、他のデータを利用してもよい。幾つかの例示的な態様では、そのストリップに対して時間的又は空間的に近くに位置する１以上の他のストリップ（近傍ストリップ）に対応するｘｙ方向シフト量データを利用することができる。例えば、（そのストリップのｘｙ方向シフト量データ及び）近傍ストリップのｘｙ方向シフト量データに対して所定の演算（統計演算、カーブフィッティングなど）を適用することによって、そのストリップの最適なｘｙ方向シフト量を決定することができる。

（Ｓ１９：一連のストリップの相対位置調整を実行する）
　画像データ補正部２３５は、一連のストリップのそれぞれについてステップＳ１８で決定された最適なｘｙ方向シフト量に基づいて、一連のストリップの相対位置調整を実行する。この相対位置調整は、ラテラル方向（ｘｙ方向）の相対位置調整であり、各ストリップをその最適なｘｙ方向シフト量に相当するベクトルだけ移動することによって行われる。本ステップの相対位置調整は、非特許文献１の「rough lateral motion correction」に相当する処理を、ステップＳ１８で決定された最適なｘｙ方向シフト量に基づき実行するものである。

（Ｓ２０：一連のストリップの精密な相対位置調整を実行する）
　画像データ補正部２３５は、ステップＳ１９のラテラル方向の相対位置調整よりも精密なラテラル方向の相対位置調整を更に実行することができる。本ステップの精密な相対位置調整は、非特許文献１の「fine lateral motion correction」に相当する処理であってよい。

　ステップＳ２０が完了したら、図１１のステップＳ８が実行され、本動作例に係る処理は終了となる（エンド）。

＜第３の動作例＞
　第３の動作例では、図１３Ａ及び図１３を参照しつつ、図１１のステップＳ５及び図１２のステップＳ１５で実行される再帰的レジストレーションの例を説明する。

（Ｓ３１：基準ストリップと対象ストリップを設定する）
　所与の初期基準ストリップ（図１１のステップＳ４で選択された複数の初期基準ストリップのいずれか１つ）を用いた再帰的レジストレーションでは、まず、データ処理部２３０が、図１１のステップＳ３で構築された一連のストリップをサイズ（面積など）に従って順序付けする。本例では、一連のストリップの個数をＮ個（Ｎは２以上の整数）であるとする。また、本例では、順序付けにより指定された順序に従って、これらＮ個のストリップを、第１のストリップ、第２のストリップ、・・・、第Ｎのストリップと呼ぶことにする。また、本例において、第１～第Ｎのストリップのうちの任意のストリップを第ｎのストリップと呼ぶことがある（ｎ＝１、２、・・・、Ｎ）。このように、本例において、第１のストリップは最も大きなストリップであり、第２のストリップは２番目に大きなストリップであり、第ｎのストリップはｎ番目に大きなストリップであり、第Ｎのストリップは最も小さなストリップである。

　所与の初期基準ストリップは第１～第Ｎのストリップのいずれかである。本例では、所与の初期基準ストリップは第１のストリップであり、且つ、初期対象ストリップは第２のストリップであるとする。なお、所与の初期基準ストリップが第１のストリップでない場合には、所与の初期基準ストリップが初期基準ストリップに設定され、且つ、第１のストリップが初期対象ストリップに設定される。この場合においても本例と同様の再帰的レジストレーションを実行できることは、当業者であれば理解できるであろう。

　基準ストリップ及び対象ストリップは、それぞれ、非特許文献１のAppendix Aにおける任意形状の画像（arbitrary shaped images）ｆ（ｒ）及びｇ（ｒ）に相当する。

（Ｓ３２：基準マスク画像と対象マスク画像を設定する）
　マスク画像生成部２３１１は、基準ストリップに対応するマスク画像（基準マスク画像）と、対象ストリップに対応するマスク画像（対象マスク画像）とを生成する。

　この段階（再帰的レジストレーションにおける第１回目のレジストレーション）では、マスク画像生成部２３１１は、初期基準ストリップである第１のストリップに対応するマスク画像と、第２のストリップに対応するマスク画像とを生成する。データ処理部２３０は、第１のストリップに対応する第１のマスク画像を初期基準マスク画像に設定し、第２のストリップに対応する第２のマスク画像を初期対象マスク画像に設定する。

　なお、基準マスク画像及び対象マスク画像は、それぞれ、非特許文献１のAppendix Aにおける矩形状の二値画像マスク（rectangular shaped binary image masks）ｍ_ｆ（ｒ）及びｍ_ｇ（ｒ）に相当する。

（Ｓ３３：基準ストリップと対象ストリップを正規化する）
　範囲調整部２３１２は、前述した要領で、ステップＳ３１で設定された基準ストリップと対象ストリップに正規化を施す。ストリップの正規化は、マスク画像の画素値の範囲に応じて行われる。この段階（再帰的レジストレーションにおける第１回目のレジストレーション）においては、範囲調整部２３１２は、第１のストリップ(基準ストリップ)及び第２のストリップ（対象ストリップ）に正規化を適用する。なお、範囲調整部２３１２が実行する処理は正規化には限定されず、前述した範囲調整処理のいずれか又はそれに類する処理であってよい。

　本例において、各マスク画像の画素値の範囲は閉区間［０，１］に含まれるものとする。特に、本例の各マスク画像は、対応するストリップの定義域に相当する領域内の画素値が１に設定され且つ他の画素の値が０に設定された二値画像であるものとする。

　範囲調整部２３１２は、例えば、ストリップの各画素の値をこのストリップにおける最大画素値で除算することにより、又は、ストリップの各画素の値をこのストリップの画素値の範囲（階調範囲）の最大値で除算することにより、このストリップの正規化を行う。

　また、データ処理部２３０（合成画像生成部２３１３）は、画素値範囲が正規化されたストリップを、対応するマスク画像と同じサイズ及び同じ形状の画像に埋め込む。正規化された第１のストリップｆ（ｒ）の埋め込み画像は、非特許文献１のAppendix Aにおける矩形状の画像（rectangular shaped image）ｆ´（ｒ）に類する画像であるが、第１のストリップｆ（ｒ）の画像エリア内の値域の絶対値（｜ｆ（ｒ）｜）が１以下に正規化されている。この埋め込み画像を同じくｆ´（ｒ）と表す。

　同様に、正規化された第２のストリップｇ（ｒ）の埋め込み画像は、矩形状の画像ｇ´（ｒ）に類する画像であるが、第２のストリップｇ（ｒ）の画像エリア内の値域の絶対値（｜ｇ（ｒ）｜）が１以下に正規化されている。この埋め込み画像を同じくｇ´（ｒ）と表す。

（Ｓ３４：基準合成画像と対象合成画像を生成する）
　合成画像生成部２３１３は、正規化されたストリップの埋め込み画像と、対応するマスク画像との合成画像を生成する。

　この段階（再帰的レジストレーションにおける第１回目のレジストレーション）においては、合成画像生成部２３１３は、正規化された第１のストリップ（初期基準ストリップ）の埋め込み画像と、第１のマスク画像（基準マスク画像）との合成画像を生成し、且つ、正規化された第２のストリップ（初期対象ストリップ）の埋め込み画像と、第２のマスク画像（対象マスク画像）との合成画像を生成する。前者の合成画像を基準合成画像と呼び、後者を対象合成画像と呼ぶ。

　ここで、基準合成画像及び対象合成画像は、それぞれ、非特許文献１のAppendix Aにおける２つの矩形状の画像の組み合わせ（combination of two rectangle shaped images）ｆ´（ｒ）ｍ_ｆ（ｒ）及びｇ´（ｒ）ｍ_ｇ（ｒ）に相当する。ただし、前述したように、第１のストリップｆ（ｒ）の画像エリア内の値域の絶対値（｜ｆ（ｒ）｜）は１以下であり、且つ、第２のストリップｇ（ｒ）の画像エリア内の値域の絶対値（｜ｇ（ｒ）｜）は１以下である。

　図１４Ａの画像３１１は、正規化された第１のストリップｆ（ｒ）の埋め込み画像ｆ´（ｒ）の例であり、画像３２１は、第１のマスク画像の例である。これら２つの画像３１１及び３２１を組み合わせることによって基準合成画像ｆ´（ｒ）ｍ_ｆ（ｒ）が得られる。

　同様に、図１４Ｂの画像３１２は、正規化された第２のストリップｇ（ｒ）の埋め込み画像ｇ´（ｒ）の例であり、画像３２２は、第２のマスク画像の例である。これら２つの画像３１２及び３２２を組み合わせることによって基準合成画像ｇ´（ｒ）ｍ_ｇ（ｒ）が得られる。

（Ｓ３５：複数の相互相関関数を算出する）
　相互相関関数算出部２３１４は、ステップＳ３４で生成された基準合成画像及び対象合成画像に基づいて複数の相互相関関数を算出する。本例では、相互相関関数算出部２３１４は、ステップＳ３４で生成された基準合成画像及び対象合成画像に基づいて、非特許文献１の式（３３）に含まれる６個の相互相関関数を算出する。

（Ｓ３６：相関係数を算出する）
　相関係数算出部２３１５は、ステップＳ３５で算出された複数の相互相関関数に基づいて相関係数を算出する。本例では、相関係数算出部２３１５は、非特許文献１の式（３３）に従って、ステップＳ３５で算出された複数の相互相関関数から相関係数（ρ（ｒ´））を算出する。

（Ｓ３７：相関係数を記録する）
　シフト量記録部２３３は、ステップＳ３６で算出された相関係数（ρ（ｒ´））を、次のステップＳ３８で算出されるｘｙ方向シフト量の付帯情報として記録することができる。前述したように、この付帯情報には、ｘｙ方向シフト量データの識別子（共通識別子、個別識別子）が含まれる。

（Ｓ３８：基準ストリップと対象ストリップとの間のｘｙ方向シフト量を算出する）
　ｘｙシフト量算出部２３１６は、ステップＳ３６で算出された相関係数に基づいて、第１のストリップ（初期基準ストリップ）ｆ（ｒ）と第２のストリップ（初期対象ストリップ）ｇ（ｒ）との間のｘｙ方向シフト量を算出する。本例では、ステップＳ３６で算出された相関係数ρ（ｒ´）のピークを検出することにより、基準ストリップｆ（ｒ）と対象ストリップｇ（ｒ）との間の相対的なｘｙ方向シフト量（Δｘ、Δｙ）が求められる。

（Ｓ３９：ｘｙ方向シフト量を記録する）
　基準ストリップと対象ストリップとのペアについてステップＳ３８で算出されたｘｙ方向シフト量（ｘｙ方向シフト量データ）は、シフト量記録部２３３によって付帯情報とともに記録される。

（Ｓ４０：基準ストリップと対象ストリップとの間の相対位置調整を行う）
　レジストレーション部２３１７は、ステップＳ３８で算出されたｘｙ方向シフト量に基づいて、基準ストリップと対象ストリップとの間の相対位置調整を行う。この相対位置調整は、ラテラル方向（ｘｙ方向）の相対位置調整である（rough lateral motion correction）。この段階（再帰的レジストレーションにおける第１回目のレジストレーション）では、第１のストリップｆ（ｒ）と第２のストリップｇ（ｒ）との間の相対位置調整が実行される。

　更に、レジストレーション部２３１７は、前述した「fine lateral motion correction」などの任意のｘｙ方向のレジストレーションを基準ストリップ及び対象ストリップに対して適用してもよい。

（Ｓ４１：基準ストリップと対象ストリップとのマージ画像を構築する）
　マージ処理部２３１８は、ステップＳ４０で相対位置調整がなされた基準ストリップと対象ストリップとのマージ画像を構築する。図１４Ｃの画像３３０は、第１のストリップｆ（ｒ）と第２のストリップｇ（ｒ）とのマージ画像の例を示す。

　ここで、ステップＳ３５～Ｓ４１の処理の実装例を説明する。まず、ｆ´（ｒ）ｍ_ｆ（ｒ）、ｇ´（ｒ）ｍ_ｇ（ｒ）、ｍ_ｆ（ｒ）、ｍ_ｇ（ｒ）、（ｆ´（ｒ）ｍ_ｆ（ｒ））^２、及び（ｇ´（ｒ）ｍ_ｇ（ｒ））^２のそれぞれの虚部（imaginary part）を０に設定することによって、実部（real part）のみにする。

　次に、ｆ´（ｒ）ｍ_ｆ（ｒ）、ｇ´（ｒ）ｍ_ｇ（ｒ）、ｍ_ｆ（ｒ）、ｍ_ｇ（ｒ）、（ｆ´（ｒ）ｍ_ｆ（ｒ））^２、及び（ｇ´（ｒ）ｍ_ｇ（ｒ））^２のそれぞれの実部に対して高速フーリエ変換（ＦＦＴ）を適用する。

　次に、これらの高速フーリエ変換で得られた関数群に基づいて、図１３ＡのステップＳ３５に示された６個の相互相関関数を導出する。

　次に、導出された６個の相互相関関数のそれぞれに逆高速フーリエ変換（ＩＦＦＴ）を適用する。

　次に、非特許文献１の式（３３）の相関係数ρ（ｒ´）を算出する。

　そして、相関係数ρ（ｒ´）のピーク位置を特定することによって基準ストリップｆ（ｒ）と対象ストリップｇ（ｒ）との間の相対的なシフト量（Δｘ、Δｙ）を算出し、このシフト量に基づいて基準ストリップｆ（ｒ）と対象ストリップｇ（ｒ）との間のレジストレーション及びマージ処理を実行する。

　このような一連の処理をＮ個のストリップに対して逐次に適用することにより、Ｎ個のストリップの全てのマージ画像が得られる。また、各ストリップペアについて求められたｘｙ方向シフト量（Δｘ、Δｙ）を時間情報とともに（つまり、時間情報に関連付けて）保存することができる。

（Ｓ４２：全てのストリップを処理したか？）
　ステップＳ３で得られた一連のストリップ（Ｎ個のストリップ）の全てに対し、前述した順序に従って逐次に、ステップＳ３１～Ｓ４１の一連の処理が実行される。

　Ｎが３以上の場合において、ステップＳ４１で第１のストリップと第２のストリップとのマージ画像が作成された場合（Ｓ４２：Ｎｏ）、処理はステップＳ３１に戻る。

　第２回目のステップＳ３１（再帰的レジストレーションにおける第２回目のレジストレーションのステップＳ３１）では、第１のストリップと第２のストリップとのマージ画像が新たな基準ストリップに設定され、且つ、第３のストリップが新たな対象ストリップに設定される。新たな基準ストリップと新たな対象ストリップに基づきステップＳ３２～Ｓ４１の処理を実行することによって、新たな基準ストリップと新たな対象ストリップとのｘｙ方向シフト量データや相関係数が生成されて記録されるとともに、新たな基準ストリップと新たな対象ストリップとのマージ画像が得られる。この新たなマージ画像は、第１～第３のストリップのマージ画像である。

　このような一連の処理をＮ個のストリップに対して逐次に適用することにより、Ｎ個のストリップの全てのマージ画像が得られる。

　ステップＳ３で得られた一連のストリップの全ての処理が完了したら（Ｓ４２：Ｙｅｓ）、次の初期基準ストリップに基づく再帰的レジストレーションに移行する（図１２のステップＳ１６）。

＜変形例＞
　上記した実施形態に係る眼科検査装置１の変形例について説明する。上記した例示的な態様では、冗長的データ収集で取得されたデータセットから構築された画像セットから複数の初期基準画像を選択するために画像サイズを参照しており、主に、リサジュースキャンを用いて得られた複数のストリップ（一連のストリップ）の複数のサブセットのそれぞれから最大面積の初期基準ストリップを選択し、ストリップ間の相関係数を利用して複数のレジストレーションを実行している。

　このように相関係数を利用する場合において、初期基準画像に対する相関係数が小さいデータ（Ａライン、サイクル、ストリップなど）を廃棄してモーションアーティファクト補正を実行することができる。相関が低いことを理由として廃棄されるデータの量は、初期基準画像の選択に依存している。すなわち、或る初期基準画像を起点とした再帰的レジストレーションでは比較的多量のデータが廃棄される一方、別の初期基準画像を起点とした再帰的レジストレーションでは比較的少数のデータしか廃棄されないことがある。特に、サイズが大きい初期基準画像を採用すれば廃棄データ量が少なくなるとも限らない。

　本変形例は、発明者らが見出したこのような知見に基づくものであり、廃棄データ量が少なくなるように初期基準画像の選択を行うものである。より具体的には、本変形例は、サイズが大きい画像から順に廃棄データ量の評価（例えば、廃棄データ量が所定の閾値未満であるか否かの判定）を行うことによって、できるだけサイズが大きく且つできるだけ廃棄データ量が少なくなる画像を初期基準画像として選択するものである。

　本変形例に係る処理は、例えば以下のような複数の工程を含む。まず、本例は、図１１の動作例と同様に、リサジュースキャンによりボリュームを収集し、一連のサブボリュームに分割し、一連のストリップを構築する。次に、本例は、図１２の動作例と同様に、一連のストリップを複数のストリップ群に分ける。

　次に、本例は、各ストリップ群に含まれる複数のストリップに対して、サイズ（面積）が大きい順に順序付けを行う。これら複数のストリップを、この順序付けで付された順序にしたがって、第１のストリップ、第２のストリップ、・・・と呼ぶことにする。

　次に、本例は、第１のストリップを基準として再帰的レジストレーション（非特許文献１の「rough lateral motion correction」）を実行し、それにより得られた画像におけるＡラインの使用率を算出する。

　Ａライン使用率は、例えば、その画像の構築に使用されたＡラインの個数（換言すると、その画像に含まれるＡラインの個数、更に換言すると、再帰的レジストレーションの過程で廃棄されなかったＡラインの個数）を、元のボリュームにおけるＡラインの個数で除算した商の値であってよい。このＡライン使用率は、所定の閾値以上の相関係数の値を有するＡラインの個数を、Ａラインの総数で除算した商の値に等しい。よって、図１３Ａ及び図１３Ｂの動作例のように相関係数の記録を行う場合には、再帰的レジストレーションにおいて記録された複数の相関係数に基づいてＡライン使用率を算出することができる。なお、評価に用いられる指標はＡライン使用率に限定されず、Ａライン使用率と同値な値（例えば、Ａライン不使用率）、サイクルの使用率又は不使用率、ストリップの使用率又は不使用率などであってもよい。

　本例は、第１のストリップにおけるＡライン使用率が所定の閾値以上である場合、このストリップ群から選択される初期基準ストリップとして第１のストリップを指定する。他方、第１のストリップにおけるＡライン使用率が所定の閾値未満である場合、本例は、次にサイズが大きい第２のストリップの処理に移行する。

　更に、本例は、第１のストリップの処理と同じ処理を第２のストリップに対して適用する。このように、サイズの大きさの順にストリップにおけるＡライン使用率を評価することにより、サイズ及びＡライン使用率の双方を考慮した初期基準画像の選択が可能になる。

　なお、ストリップ群の少なくとも一部のストリップのＡライン使用率を算出し、Ａライン使用率が最大のストリップを初期基準ストリップに指定するようにしてもよい。このとき、所定の閾値以上のサイズを有する１つ以上のストリップをストリップ群から抽出し、抽出された１つ以上のストリップのうちＡライン使用率が最大のストリップを初期基準ストリップに指定するようにしてもよい。

＜効果＞
　本実施形態に係る眼科検査装置１の幾つかの効果について説明する。

　眼科検査装置１は、光走査を用いてサンプルのイメージングを行う走査型イメージング装置であって、その機能要素として、データセット取得部と、画像セット生成部と、画像選択部と、画像位置調整部とを含む。

　データセット取得部は、光走査によりサンプルのデータを冗長的に収集してデータセットを取得するように構成されている。本実施形態において、データセット取得部は光スキャナ４４、ＯＣＴユニット１００、及び制御部２１０を含んでいてよく、光走査はリサジュースキャンであってよいが、これらは非限定的な例である。

　画像セット生成部は、データセット取得部により取得されたデータセットに基づいて画像セットを生成するように構成されている。生成される画像セットは、位置的冗長性を有する複数の画像である。本実施形態において、画像セット生成部は画像データ構築部２２０を含んでいてよく、画像セットは複数のストリップ（一連のストリップ）や複数のサブボリューム（一連のサブボリューム）であってよいが、これらは非限定的な例である。

　画像選択部は、画像セット生成部により生成された画像セットから複数の画像を選択する。選択される複数の画像は初期基準画像として用いられる。本実施形態において、画像選択部は画像選択部２２５を含んでいてよく、画像セットから選択される複数の画像は複数の初期基準ストリップであってよいが、これらは非限定的な例である。

　画像位置調整量は、画像選択部により画像セットから選択された複数の画像をそれぞれ基準とした複数のレジストレーションをこの画像セットに適用することによって、この画像セットの相対位置調整を実行するように構成されている。本実施形態において、画像位置調整部は画像データ処理部２３１を含んでいてよく、レジストレーションは相関係数を用いた再帰的レジストレーションであってよいが、これらは非限定的な例である。

　前述したように、従来の技術では、単一の初期基準画像のみを使用してレジストレーションを行っているため、初期基準画像の選択が最終画像の品質に悪影響を与えるケースがあったが、本実施形態によれば、複数の初期基準画像を使用して複数のレジストレーションを行うことによってこの問題を解消し、モーションアーティファクト補正の誤差の低減を図ることができ、モーションアーティファクト補正の向上を図ることが可能になる。

　本実施形態に係る眼科検査装置１は、更に、次のような非限定的な特徴（非限定的な機能、非限定的な構成など）を有している。これらの非限定的特徴が、モーションアーティファクト補正の更なる向上に寄与するものであり、また、非限定的特徴それぞれに応じた作用及び効果を奏するものであることは、当業者であれば理解できるであろう。なお、本実施形態に係る眼科検査装置１の特徴は以下の事項に限定されるものではない。

　画像位置調整部は、画像セット生成部により生成された画像セットに対して複数のレジストレーションを適用することによって、この画像セットの複数の相対位置情報を生成するように構成されていてよい。相対位置情報は、この画像セットに含まれる複数の画像の間の相対位置関係を表す情報である。更に、画像位置調整部は、生成された複数の相対位置情報に基づいてこの画像セットの相対位置調整を実行するように構成されていてよい。

　画像位置調整部は、複数のレジストレーションのそれぞれにおいて、画像選択部により選択された複数の画像の１つを初期基準画像としてこの画像セットに再帰的なレジストレーションを適用することによって、この画像セットの相対位置情報を生成するように構成されていてよい。

　画像位置調整部は、複数のレジストレーションで生成された複数の相対位置情報に基づいて、この画像セットの相対位置調整のための調整量を決定するように構成されていてよい。この調整量は、例えば、この画像セットに含まれる各画像の尤もらしい位置ずれ量（最適な位置補正量）であってよい。

　画像選択部は、データセット取得部により実行される光走査の時間軸に基づいて、画像セット生成部により生成された画像セットから複数の画像を選択するように構成されていてよい。

　画像選択部は、画像セット生成部により生成された画像セットを、データセット取得部により実行される光走査の時間軸の複数の区間に対応する複数のサブセットに分割するように構成されていてよい。更に、画像選択部は、得られた複数のサブセットのそれぞれから画像を選択することによって、この画像セットから複数の画像を選択するように構成されていてよい。

　画像選択部は、データセット取得部により実行される光走査の適用エリアに基づいて、画像セット生成部により生成された画像セットから複数の画像を選択するように構成されていてよい。

　画像選択部は、画像セット生成部により生成された画像セットを、データセット取得部により実行される光走査の適用エリアの複数のサブエリアに対応する複数のサブセットに分割するように構成されていてよい。更に、画像選択部は、得られた複数のサブセットのそれぞれから画像を選択することによって、この画像セットから複数の画像を選択するように構成されていてよい。

　画像選択部は、画像サイズに基づいて複数のサブセットのそれぞれから画像を選択するように構成されていてよい。

　画像選択部は、複数のサブセットのそれぞれから、そのサブセットにおける最大サイズの画像を選択するように構成されていてよい。

　画像選択部は、複数のサブセットのそれぞれについて、そのサブセットに含まれる画像（そのサブセットに含まれる一部の画像又は全部の画像）に対して、画像サイズに基づく順序付けを行うように構成されていてよい。更に、画像選択部は、複数のサブセットのそれぞれから、そのサブセットにおける画像に対して順序付けにより付された順序に基づいて画像を選択するように構成されていてよい。

　画像選択部は、複数のサブセットのうちの第１サブセットから選択された第１画像を基準としたレジストレーションの成否を判定するように構成されていてよい。更に、画像選択部は、第１画像を基準としたレジストレーションは失敗であると判定された場合に、この第１サブセットにおける画像に対して順序付けにより付された順序に基づいて第１サブセットから第２画像を選択するように構成されていてよい。ここで、第２画像は第１画像とは異なる画像である。

　画像選択部は、画像セットに含まれる画像（その画像セットに含まれる一部の画像又は全部の画像）に対して画像サイズに基づく順序付けを行うように構成されていてよい。更に、画像選択部は、この順序付けにより画像に付された順序に基づいて、この画像セットから複数の画像を選択するように構成されていてよい。

　画像選択部は、画像セットから選択された複数の画像のうちの一の画像を基準としたレジストレーションの成否を判定するように構成されていてよい。更に、画像選択部は、この一の画像を基準としたレジストレーションは失敗であると判定された場合に、この画像セットにおける画像に対して順序付けにより付された順序に基づいてこの画像セットから上記複数の画像のいずれとも異なる画像を新たに選択するように構成されていてよい。

　本実施形態に係る眼科検査装置１は、運動情報生成部を更に含んでいてよい。運動情報生成部は、データセット取得部によって光走査がサンプルに適用されているときのサンプルの運動を表す運動情報を生成するように構成されている。本実施形態において、運動情報生成部は、データ処理部２３０を含んでいてよく、眼底カメラユニット２を含んでいてよく、及び／又は、前眼部カメラを含んでいてよいが、これらは非限定的な例である。更に、画像選択部は、運動情報生成部により生成された運動情報に基づいて画像セットから複数の画像を選択するように構成されていてよい。

　画像選択部は、運動情報生成部により生成された運動情報から所定の運動パラメータの値を算出し、算出された値を評価し、この評価の結果に基づいて画像セットから複数の画像を選択するように構成されていてよい。

　画像選択部により値が算出される運動パラメータは、運動の大きさ及び頻度の少なくとも一方を含んでいてよい。更に、画像選択部は、この運動パラメータの値を所定の閾値と比較し、この閾値以上の値に対応する期間に収集されたデータに基づく画像を画像セットから除外したサブセットから複数の画像を選択するように構成されていてよい。ここで、閾値以上の値に対応する期間は、大きさや頻度が閾値以上である運動が発生した時点の付近の期間であってよいが、これは非限定的な例である。

　データセット取得部は、光走査のための光を互いに異なる第１方向及び第２方向に偏向可能な偏向器を含んでいてよい。本実施形態において、この偏向器は光スキャナ４４であるが、これは非限定的な例である。更に、サンプルのデータを冗長的に収集するための光走査を実行するために、データセット取得部は、第１方向における偏向方向の変化を第１周期で繰り返しつつ、第２方向における偏向方向の変化を前記第１周期と異なる第２周期で繰り返すように構成されていてよい。

＜第２の実施形態＞
　第２の実施形態について説明する。第１の実施形態では、光走査を実行可能に構成された走査型イメージング装置について説明したが、本実施形態では、光走査により収集されたデータセットを外部から受け付けてこれを処理するように構成された情報処理装置について説明する。

　本実施形態の情報処理装置は、走査型イメージング装置の一部として構成されていてもよいが、そのような構成は第１の実施形態と実質的に同一の構成であるから、そのような場合は除外して以下の説明を行う。

　本実施形態の情報処理装置は、光走査を実行する機能要素（第１～第６の実施形態におけるデータセット取得部）を有していないことを除いて、第１の実施形態に係る走査型イメージング装置（眼科検査装置）と同様の構成を備えていてもよい。

　第１の実施形態において説明された任意の事項を本実施形態の情報処理装置に組み合わせることが可能である。情報処理装置は、組み合わせられた事項に応じた作用及び効果を奏する。

　本実施形態の情報処理装置の例示的な態様の構成を図１５に示す。本態様の情報処理装置４００は、受付部４０１と、画像セット生成部４０２と、画像選択部４０３と、画像位置調整部４０４とを含んでいる。

　受付部４０１は、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットを受け付ける（データセット受付部）。すなわち、受付部４０１は、リサジュースキャンなどの光走査を実行可能な走査型イメージング装置によってサンプルから収集されたデータセットを受け付ける。

　受付部４０１は、外部装置からデータを取得する機能を有する。この外部装置は、例えば、コンピュータ、記憶装置、記録媒体、情報システム（例えば、病院情報システム、電子カルテシステム、画像アーカイビングシステム）などであってよい。また、外部装置は、例えば、情報処理装置４００に直接に接続された装置、ローカルエリアネットワーク（ＬＡＮ）を介して情報処理装置４００に接続された装置、ワイドエリアネットワーク（ＷＡＮ）を介して情報処理装置４００に接続された装置などであってよい。受付部４０１は、例えば、通信回線を介してデータ通信を行うための通信インターフェイス、及び／又は、記録媒体からデータを読み取るためのドライブ装置を含んでいてもよい。

　画像セット生成部４０２は、受付部４０１により受け付けられたデータセットに基づいて画像セットを生成するように構成されている。画像セット生成部４０２は、第１の実施形態の画像セット生成部（画像データ構築部２２０）と同様の構成を有していてよい。

　画像選択部４０３は、画像セット生成部４０２により生成された画像セットから複数の画像を選択するように構成されている。画像選択部４０３は、第１の実施形態の画像選択部２２５と同様の構成を有していてよい。

　画像位置調整部４０４は、画像選択部４０３により選択された複数の画像をそれぞれ基準とした複数のレジストレーションを、画像セット生成部４０２により生成された画像セットに適用することによって、この画像セットの相対位置調整を実行するように構成されている。画像位置調整部４０４は、第１の実施形態の画像位置調整部（画像データ処理部２３１）と同様の構成を有していてよい。

　このような情報処理装置４００によれば、第１の実施形態と同様に、モーションアーティファクト補正の誤差の低減を図ることができ、モーションアーティファクト補正の向上を図ることが可能になる。

　本実施形態の情報処理装置の他の例示的な態様の構成を図１６に示す。本態様の情報処理装置５００は、受付部５０１と、画像選択部５０２と、画像位置調整部５０３とを含んでいる。

　受付部５０１は、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットに基づいて生成された画像セットを受け付ける（画像セット受付部）。すなわち、受付部５０１は、リサジュースキャンなどの光走査を実行可能な走査型イメージング装置によってサンプルから収集されたデータセットから構築された画像セットを受け付ける。受付部５０１は、前述の受付部４０１と同様の構成を有していてよい。

　画像選択部５０２は、受付部５０１により受け付けられた画像セットから複数の画像を選択するように構成されている。画像選択部５０２は、第１の実施形態の画像選択部２２５と同様の構成を有していてよい。

　画像位置調整部５０３は、画像選択部５０２により選択された複数の画像をそれぞれ基準とした複数のレジストレーションを、受付部５０１により受け付けられた画像セットに適用することによって、この画像セットの相対位置調整を実行するように構成されている。画像位置調整部５０３は、第１の実施形態の画像位置調整部（画像データ処理部２３１）と同様の構成を有していてよい。

　このような情報処理装置５００によれば、第１の実施形態と同様に、モーションアーティファクト補正の誤差の低減を図ることができ、モーションアーティファクト補正の向上を図ることが可能になる。

＜他の実施形態＞
　本開示に係る実施形態は、第１及び第２の実施形態及びそれらの変形例に限定されない。例えば、第１及び第２の実施形態は、各種の方法に係る実施形態、各種のプログラムに係る実施形態、各種の記録媒体に係る実施形態などを提供するものである。

　例えば、前述した第１の実施形態は、走査型イメージング装置を制御する方法の実施形態を提供する。この実施形態は、プロセッサと、光走査を行うスキャナとを含む走査型イメージング装置を制御する方法である。このスキャナは第１の実施形態のデータセット取得部であってよく、このプロセッサは第１の実施形態の制御部２１０、画像データ構築部２２０、及びデータ処理部２３０であってよい。

　本実施形態に係る走査型イメージング装置の制御方法は、次の４つの工程を含んでいてよい。第１工程の制御は、スキャナに、サンプルのデータを冗長的に収集させてデータセットを取得させる。第２工程の制御は、プロセッサに、第１工程で取得されたデータセットに基づいて画像セットを生成させる。第３工程の制御は、プロセッサに、第２工程で生成された画像セットから複数の画像を選択させる。第４工程の制御は、プロセッサに、第３工程で選択された複数の画像をそれぞれ基準とした複数のレジストレーションを、第２工程で生成された画像セットに適用させることによって、この画像セットの相対位置調整を実行させる。

　第１又は第２の実施形態において説明された任意の事項を、本実施形態に係る走査型イメージング装置の制御方法に組み合わせることが可能である。走査型イメージング装置の制御方法は、組み合わせられた事項に応じた作用及び効果を奏する。

　このような走査型イメージング装置の制御方法によれば、第１の実施形態と同様に、モーションアーティファクト補正の誤差の低減を図ることができ、モーションアーティファクト補正の向上を図ることが可能になる。

　前述した第２の実施形態は、情報処理装置を制御する方法の実施形態を提供する。この実施形態は、プロセッサを含む情報処理装置を制御する方法である。このプロセッサは第１の実施形態の制御部２１０、画像データ構築部２２０、及びデータ処理部２３０であってよい。

　本実施形態に係る情報処理装置の制御方法は、次の４つの工程を含んでいてよい。第１工程の制御は、プロセッサに、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットを受け付けさせる。第２工程の制御は、プロセッサに、第１工程で受け付けられたデータセットに基づいて画像セットを生成させる。第３工程の制御は、プロセッサに、第２工程で生成された画像セットから複数の画像を選択させる。第４工程の制御は、プロセッサに、第３工程で選択された複数の画像をそれぞれ基準とした複数のレジストレーションを、第２工程で生成された画像セットに適用させることによって、この画像セットの相対位置調整を実行させる。

　第１又は第２の実施形態において説明された任意の事項を、本実施形態に係る情報処理装置の制御方法に組み合わせることが可能である。情報処理装置の制御方法は、組み合わせられた事項に応じた作用及び効果を奏する。

　このような情報処理装置の制御方法によれば、第１の実施形態と同様に、モーションアーティファクト補正の誤差の低減を図ることができ、モーションアーティファクト補正の向上を図ることが可能になる。

　前述した第２の実施形態は、情報処理装置を制御する方法の他の実施形態を提供する。この実施形態は、プロセッサを含む情報処理装置を制御する方法である。このプロセッサは第１の実施形態の制御部２１０及びデータ処理部２３０であってよい。

　本実施形態に係る情報処理装置の制御方法は、次の３つの工程を含んでいてよい。第１工程の制御は、プロセッサに、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットに基づいて生成された画像セットを受け付けさせる。第２工程の制御は、プロセッサに、第１工程で受け付けられた画像セットから複数の画像を選択させる。第３工程の制御は、プロセッサに、第２工程で選択された複数の画像をそれぞれ基準とした複数のレジストレーションを、第１工程で受け付けられた画像セットに適用させることによって、この画像セットの相対位置調整を実行させる。

　第１及び第２の実施形態は、情報を処理する方法の実施形態を提供する。

　本実施形態に係る情報処理方法は、光走査によりサンプルのデータを冗長的に収集して取得されたデータセットを処理する方法であって、次の３つの工程を含んでいてよい。第１工程は、このデータセットに基づいて画像セットを生成する。第２工程は、第１工程で生成された画像セットから複数の画像を選択する。第３工程は、第２工程で選択された複数の画像をそれぞれ基準とした複数のレジストレーションを、第１工程で生成された画像セットに適用することによって、この画像セットの相対位置調整を実行する。

　第１又は第２の実施形態において説明された任意の事項を、本実施形態に係る情報処理方法に組み合わせることが可能である。情報処理方法は、組み合わせられた事項に応じた作用及び効果を奏する。

　このような情報処理方法によれば、第１の実施形態と同様に、モーションアーティファクト補正の誤差の低減を図ることができ、モーションアーティファクト補正の向上を図ることが可能になる。

　本開示は、実施形態に係る方法をコンピュータに実行させるプログラムを提供する。幾つかの例示的な実施形態は、走査型イメージング装置の制御方法をコンピュータに実行させるプログラムを提供する。幾つかの例示的な実施形態は、情報処理装置の制御方法をコンピュータに実行させるプログラムを提供する。幾つかの例示的な実施形態は、情報処理方法をコンピュータに実行させるプログラムを提供する。

　本開示は、実施形態に係るプログラムが記録されたコンピュータ可読な非一時的記録媒体を提供する。幾つかの例示的な実施形態は、走査型イメージング装置の制御方法をコンピュータに実行させるプログラムが記録されたコンピュータ可読な非一時的記録媒体を提供する。幾つかの例示的な実施形態は、情報処理装置の制御方法をコンピュータに実行させるプログラムが記録されたコンピュータ可読な非一時的記録媒体を提供する。幾つかの例示的な実施形態は、情報処理方法をコンピュータに実行させるプログラムが記録されたコンピュータ可読な非一時的記録媒体を提供する。実施形態の非一時的記録媒体は任意の形態であってよく、その例として、磁気ディスク、光ディスク、光磁気ディスク、半導体メモリなどがある。

　本開示は、幾つかの例示的な態様を例示するものに過ぎず、発明の限定を意図したものではない。この発明を実施しようとする者は、この発明の要旨の範囲内における任意の変形（省略、置換、付加など）を施すことが可能である。

１　眼科検査装置
２　眼底カメラユニット
１００　ＯＣＴユニット
２１０　制御部
２２０　画像データ構築部
２２５　画像選択部
２３０　データ処理部
２３１　画像データ処理部
２３３　シフト量記録部
２３４　調整量算出部
２３５　画像データ補正部

Claims

　光走査を用いてサンプルのイメージングを行う走査型イメージング装置であって、
　前記光走査により前記サンプルのデータを冗長的に収集してデータセットを取得するデータセット取得部と、
　前記データセットに基づいて画像セットを生成する画像セット生成部と、
　前記画像セットから複数の画像を選択する画像選択部と、
　前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用することによって前記画像セットの相対位置調整を実行する画像位置調整部と
　を含む、走査型イメージング装置。
　前記画像位置調整部は、前記複数のレジストレーションを前記画像セットに適用することによって前記画像セットの複数の相対位置情報を生成し、前記複数の相対位置情報に基づいて前記相対位置調整を実行する、
　請求項１の走査型イメージング装置。
　前記画像位置調整部は、前記複数のレジストレーションのそれぞれにおいて、前記複数の画像の１つを初期基準画像として前記画像セットに再帰的なレジストレーションを適用することによって前記画像セットの相対位置情報を生成する、
　請求項２の走査型イメージング装置。
　前記画像位置調整部は、前記複数の相対位置情報に基づいて前記相対位置調整のための調整量を決定する、
　請求項２又は３の走査型イメージング装置。
　前記画像選択部は、前記光走査の時間軸に基づいて前記画像セットから前記複数の画像を選択する、
　請求項１～４のいずれかの走査型イメージング装置。
　前記画像選択部は、前記画像セットを前記時間軸の複数の区間に対応する複数のサブセットに分割し、前記複数のサブセットのそれぞれから画像を選択することによって、前記画像セットから前記複数の画像を選択する、
　請求項５の走査型イメージング装置。
　前記画像選択部は、前記光走査の適用エリアに基づいて前記画像セットから前記複数の画像を選択する、
　請求項１～４のいずれかの走査型イメージング装置。
　前記画像選択部は、前記画像セットを前記適用エリアの複数のサブエリアに対応する複数のサブセットに分割し、前記複数のサブセットのそれぞれから画像を選択することによって、前記画像セットから前記複数の画像を選択する、
　請求項７の走査型イメージング装置。
　前記画像選択部は、画像サイズに基づいて前記複数のサブセットのそれぞれから画像を選択する、
　請求項６又は８の走査型イメージング装置。
　前記画像選択部は、前記複数のサブセットのそれぞれから、当該サブセットにおける最大サイズの画像を選択する、
　請求項９の走査型イメージング装置。
　前記画像選択部は、前記複数のサブセットのそれぞれに含まれる画像に対して画像サイズに基づく順序付けを行い、前記複数のサブセットのそれぞれから、当該サブセットにおける画像に対して前記順序付けにより付された順序に基づいて画像を選択する、
　請求項９の走査型イメージング装置。
　前記画像選択部は、前記複数のサブセットのうちの第１サブセットから選択された第１画像を基準としたレジストレーションの成否を判定し、失敗と判定された場合、前記第１サブセットにおける画像に対して前記順序付けにより付された順序に基づいて前記第１サブセットから第２画像を選択する、
　請求項１１の走査型イメージング装置。
　前記画像選択部は、前記画像セットに含まれる画像に対して画像サイズに基づく順序付けを行い、前記順序付けにより画像に付された順序に基づいて前記画像セットから前記複数の画像を選択する、
　請求項１～８のいずれかの走査型イメージング装置。
　前記画像選択部は、前記画像セットから選択された前記複数の画像のうちの一の画像を基準としたレジストレーションの成否を判定し、失敗と判定された場合、前記画像セットにおける画像に対して前記順序付けにより付された順序に基づいて前記画像セットから前記複数の画像のいずれとも異なる画像を新たに選択する、
　請求項１３の走査型イメージング装置。
　前記光走査が前記サンプルに適用されているときの前記サンプルの運動を表す運動情報を生成する運動情報生成部を更に含み、
　前記画像選択部は、前記運動情報に基づいて前記画像セットから前記複数の画像を選択する、
　請求項１～１４のいずれかの走査型イメージング装置。
　前記画像選択部は、前記運動情報から所定の運動パラメータの値を算出し、前記値を評価し、前記評価の結果に基づいて前記画像セットから前記複数の画像を選択する、
　請求項１５の走査型イメージング装置。
　前記運動パラメータは、運動の大きさ及び頻度の少なくとも一方を含み、
　前記画像選択部は、前記運動パラメータの前記値を所定の閾値と比較し、前記閾値以上の前記値に対応する期間に収集されたデータに基づく画像を前記画像セットから除外したサブセットから前記複数の画像を選択する、
　請求項１６の走査型イメージング装置。
　前記データセット取得部は、
　前記光走査のための光を互いに異なる第１方向及び第２方向に偏向可能な偏向器を含み、
　前記第１方向における偏向方向の変化を第１周期で繰り返しつつ、前記第２方向における偏向方向の変化を前記第１周期と異なる第２周期で繰り返す、
　請求項１～１７のいずれかの走査型イメージング装置。
　光走査によりサンプルからデータを冗長的に収集して取得されたデータセットを受け付けるデータセット受付部と、
　前記データセットに基づいて画像セットを生成する画像セット生成部と、
　前記画像セットから複数の画像を選択する画像選択部と、
　前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用することによって前記画像セットの相対位置調整を実行する画像位置調整部と
　を含む、情報処理装置。
　光走査によりサンプルからデータを冗長的に収集して取得されたデータセットに基づいて生成された画像セットを受け付ける画像セット受付部と、
　前記画像セットから複数の画像を選択する画像選択部と、
　前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用することによって前記画像セットの相対位置調整を実行する画像位置調整部と
　を含む、情報処理装置。
　プロセッサと光走査を行うスキャナとを含む走査型イメージング装置を制御する方法であって、
　前記スキャナに、前記サンプルのデータを冗長的に収集させてデータセットを取得させ、
　前記プロセッサに、前記データセットに基づいて画像セットを生成させ、
　前記プロセッサに、前記画像セットから複数の画像を選択させ、
　前記プロセッサに、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用させることによって前記画像セットの相対位置調整を実行させる、
　方法。
　プロセッサを含む情報処理装置を制御する方法であって、
　前記プロセッサに、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットを受け付けさせ、
　前記プロセッサに、前記データセットに基づいて画像セットを生成させ、
　前記プロセッサに、前記画像セットから複数の画像を選択させ、
　前記プロセッサに、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用させることによって前記画像セットの相対位置調整を実行させる、
　方法。
　プロセッサを含む情報処理装置を制御する方法であって、
　前記プロセッサに、光走査によりサンプルからデータを冗長的に収集して取得されたデータセットに基づいて生成された画像セットを受け付けさせ、
　前記プロセッサに、前記画像セットから複数の画像を選択させ、
　前記プロセッサに、前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用させることによって前記画像セットの相対位置調整を実行させる、
　方法。
　光走査によりサンプルのデータを冗長的に収集して取得されたデータセットを処理する情報処理方法であって、
　前記データセットに基づいて画像セットを生成し、
　前記画像セットから複数の画像を選択し、
　前記複数の画像をそれぞれ基準とした複数のレジストレーションを前記画像セットに適用することによって前記画像セットの相対位置調整を実行する、
　方法。
　請求項２１～２４のいずれかの方法をコンピュータに実行させるプログラム。
　請求項２５のプログラムが記録されたコンピュータ可読な非一時的記録媒体。