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

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

Info

Publication number
JP2023048318A
JP2023048318A JP2021157561A JP2021157561A JP2023048318A JP 2023048318 A JP2023048318 A JP 2023048318A JP 2021157561 A JP2021157561 A JP 2021157561A JP 2021157561 A JP2021157561 A JP 2021157561A JP 2023048318 A JP2023048318 A JP 2023048318A
Authority
JP
Japan
Prior art keywords
scanning
image
unit
correlation coefficient
image data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2021157561A
Other languages
English (en)
Inventor
達夫 山口
Tatsuo Yamaguchi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Topcon Corp
Original Assignee
Topcon Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Topcon Corp filed Critical Topcon Corp
Priority to JP2021157561A priority Critical patent/JP2023048318A/ja
Publication of JP2023048318A publication Critical patent/JP2023048318A/ja
Pending legal-status Critical Current

Links

Images

Landscapes

  • Eye Examination Apparatus (AREA)
  • Image Processing (AREA)

Abstract

Figure 2023048318000001
【課題】走査型イメージングで取得されたデータの新規な利用法を提供する。
【解決手段】実施形態の眼科検査装置１は、光走査を用いて画像データを構築する走査型イメージング装置であり、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータセットを収集する。眼科検査装置１は、このデータセットに基づいてサイクル間の相関係数を算出する。更に、眼科検査装置１は、このデータセットから算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録する。
【選択図】図４Ｂ

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(ti)=A・cos(2π・(fA/n)・ti)、y(ti)=A・cos(2π・(fA・(n-2)/n2)・ti)。
ここで、「ｘ」はリサジュー曲線が定義される２次元座標系の横軸、「ｙ」は縦軸、「ｔｉ」はリサジュースキャンにおける第ｉ番目のＡラインの収集時点、「Ａ」はスキャン範囲（振幅）、「ｆＡ」はＡラインの収集レート（スキャン速度、Ａスキャンの繰り返しレート）、ｎはｘ方向（横軸方向）の各サイクルにおけるＡラインの個数をそれぞれ示す。
このような例示的なリサジュースキャンにおけるスキャンラインの分布（スキャンパターン）の例を図５に示す。
（狭義又は広義の）リサジュースキャンに含まれる任意のサイクル対（サイクルのペア）は、少なくとも１点（特に２点以上）で互いに交差する。このような交差を利用することで、リサジュースキャンにおける任意のサイクル対から収集されたデータ対の間のレジストレーションを行うことが可能となり、非特許文献１（又は非特許文献２）に開示された画像構築手法及びモーションアーティファクト補正手法を実装することが可能となる。以下、特に言及しない限り、非特許文献１に記載された手法を適用する場合について説明するが、非特許文献２に記載された手法を適用することや、非特許文献１又は２に記載された手法と同等及び／又は類似の手法を適用することも可能である。
記憶部２１２に記憶される制御情報は上記の例に限定されない。例えば、制御情報はフォーカス制御を行うための情報（フォーカス制御パラメータ）を含んでいてよい。
フォーカス制御パラメータは、ＯＣＴ合焦駆動部４３Ａに対する制御の内容を示すパラメータである。フォーカス制御パラメータの例として、測定アームの焦点位置を示すパラメータ、焦点位置の移動速度を示すパラメータ、焦点位置の移動加速度を示すパラメータなどがある。焦点位置を示すパラメータは、例えば、ＯＣＴ合焦レンズ４３の位置を示すパラメータである。焦点位置の移動速度を示すパラメータは、例えば、ＯＣＴ合焦レンズ４３の移動速度を示すパラメータである。焦点位置の移動加速度を示すパラメータは、例えば、ＯＣＴ合焦レンズ４３の移動加速度を示すパラメータである。移動速度は、一定であってもよいし、一定でなくてもよい。移動加速度についても同様である。
このようなフォーカス制御パラメータによれば、眼底Ｅｆの形状（典型的には、中心部が深く且つ周辺部が浅い凹形状）や収差分布に応じたフォーカス調整が可能になる。フォーカス制御は、例えば、走査制御（リサジュースキャンの繰り返し制御）と連係的に実行される。それにより、モーションアーティファクトが補正され、且つ、スキャン範囲の全体に亘ってピントが合った、高品質の画像が得られる。
＜走査制御部２１１１＞
走査制御部２１１１は、走査プロトコル２１２１に基づいて少なくとも光スキャナ４４を制御する。走査制御部２１１１は、走査プロトコル２１２１に基づく光スキャナ４４の制御と連係して光源ユニット１０１の制御を更に実行してもよい。走査制御部２１１１は、プロセッサを含むハードウェアと、走査プロトコル２１２１を含む走査制御ソフトウェアとの協働によって実現される。
＜画像データ構築部２２０＞
画像データ構築部２２０は、プロセッサを含み、データ収集システム１３０から入力された信号（サンプリングデータ）に基づいて、眼底ＥｆのＯＣＴ画像データを構築する。このＯＣＴ画像データ構築は、従来のフーリエドメインＯＣＴ（スウェプトソースＯＣＴ）と同様に、ノイズ除去（ノイズ低減）、フィルタリング、高速フーリエ変換（ＦＦＴ）などを含む。他のタイプのＯＣＴ手法が採用される場合、画像データ構築部２２０は、そのＯＣＴタイプに対応した公知の処理を実行することによってＯＣＴ画像データを構築する。
例えば、画像データ構築部２２０は、サンプリングデータから各スキャン点（各Ａライン）に対応する画像データ（Ａスキャン画像データ）を構築する処理を少なくとも実行するように構成されている。画像データ構築部２２０により構築される画像データは、画像化される前のデータであってもよく、例えば、深さ方向（ｚ方向、光軸方向、Ａライン方向）に沿った信号プロファイル（反射プロファイル、散乱プロファイル）であってもよい。
前述したように、本実施形態ではリサジュースキャンが眼底Ｅｆに適用される。画像データ構築部２２０は、データ処理部２３０とともに、リサジュースキャンによる収集とデータ収集システム１３０によるサンプリングとを介して取得されたデータセットに対し、例えば非特許文献１に開示された画像構築手法及びモーションアーティファクト補正手法を適用することによって、眼底Ｅｆの３次元画像データを構築することができる。
画像データ構築部２２０及び／又はデータ処理部２３０は、３次元画像データにレンダリングを適用して表示用画像を形成することができる。適用可能なレンダリング法の例として、ボリュームレンダリング、サーフェスレンダリング、最大値投影（ＭＩＰ）、最小値投影（ＭｉｎＩＰ）、多断面再構成（ＭＰＲ）などがある。
画像データ構築部２２０及び／又はデータ処理部２３０は、３次元画像データに基づいてＯＣＴ正面画像（en face OCT image）を構築することが可能である。例えば、画像データ構築部２２０及び／又はデータ処理部２３０は、３次元画像データをｚ方向（Ａライン方向、光軸方向、深さ方向）に投影することでプロジェクションデータを構築することができる。また、画像データ構築部２２０及び／又はデータ処理部２３０は、３次元画像データの一部である部分的３次元画像データをｚ方向に投影することでシャドウグラムを構築することができる。この部分的３次元画像データは、例えば、任意のセグメンテーション手法を利用して設定されてよい。セグメンテーションは、画像中の部分領域を特定する処理である。本例では、眼底Ｅｆの１つ以上の組織（部位）に相当する画像領域を特定するためにセグメンテーションを行うことができる。
シャドウグラム構築処理に限らず、実施形態のセグメンテーションは、例えば、閾値処理、エッジ検出、フィルタリング、機械学習（例えば、セマンティックセグメンテーション）など、任意の公知の技術を利用可能である。
眼科検査装置１は、ＯＣＴ血管造影（OCT-Angiography）を実行可能であってよい。ＯＣＴ血管造影は、血管が強調された画像を構築するイメージング技術である（例えば、非特許文献２、特表２０１５－５１５８９４号公報などを参照）。一般に、眼底組織（構造）は短時間の間には変化しないが、血管内部の血流は短時間の間にも変化する。ＯＣＴ血管造影では、このような時間的変化が存在する部分（血流信号）を強調して画像を生成する。なお、ＯＣＴ血管造影は、ＯＣＴモーションコントラスト撮影（motion contrast imaging）などとも呼ばれる。また、ＯＣＴ血管造影により取得される画像は、血管造影画像、アンジオグラム（angiogram）、モーションコントラスト画像などと呼ばれる。
ＯＣＴ血管造影を実行可能である場合、眼科検査装置１は、眼底Ｅｆの同じ領域を所定回数だけ繰り返しスキャンする。例えば、眼科検査装置１は、前述した走査制御（リサジュースキャンの繰り返し制御）を所定回数だけ繰り返し実行する。それにより、リサジュースキャンの適用領域から複数の３次元データ（時系列な３次元データセット）がデータ収集システム１３０によって収集される。画像データ構築部２２０及び／又はデータ処理部２３０は、この３次元データセットからモーションコントラスト画像を構築することができる。このモーションコントラスト画像は、眼底Ｅｆの血流に起因する干渉信号の時間的変化が強調された血管造影画像である。この血管造影画像は、眼底Ｅｆの血管の３次元的な分布を表現した３次元血管造影画像データである。
画像データ構築部２２０及び／又はデータ処理部２３０は、この３次元血管造影画像データから、任意の２次元血管造影画像データ及び／又は任意の擬似的３次元血管造影画像データを構築することが可能である。例えば、画像データ構築部２２０及び／又はデータ処理部２３０は、３次元血管造影画像データに多断面再構成を適用することにより、眼底Ｅｆの任意の断面を表す２次元血管造影画像データを構築することができる。また、画像データ構築部２２０及び／又はデータ処理部２３０は、３次元血管造影画像データにプロジェクション画像化又はシャドウグラム化を適用することにより、眼底Ｅｆの正面血管造影画像データを構築することができる。
本実施形態において、画像データ構築部２２０及び／又はデータ処理部２３０は、データ収集システム１３０により収集されたデータから複数のストリップを構築する。非特許文献１に記載されているように、画像データ構築部２２０及び／又はデータ処理部２３０は、リサジュースキャンで収集されたボリューム（３次元データ）を比較的大きな動きが介在しない複数のサブボリュームに分割し、各サブボリュームの正面プロジェクション画像（en face projection）を構築する。この正面プロジェクション画像がストリップである。このようにして得られた複数のストリップに対してレジストレーション及びマージ処理を適用することにより、モーションアーティファクトが補正された画像が得られる。後述するように、本例においては、レジストレーション及びマージ処理はデータ処理部２３０によって実行される。
画像データ構築部２２０は、プロセッサを含むハードウェアと、画像構築ソフトウェアとの協働によって実現される。なお、幾つかの例示的な態様において、画像データ構築部２２０とデータ処理部２３０とが一体的に構成されていてよい。
＜データ処理部２３０＞
データ処理部２３０は、プロセッサを含み、被検眼Ｅの画像に対して各種のデータ処理を適用する。例えば、データ処理部２３０は、プロセッサを含むハードウェアと、データ処理ソフトウェアとの協働によって実現される。データ処理部２３０に含まれる各機能要素も同様である。
データ処理部２３０は、眼底Ｅｆについて取得された２つの画像の間の位置合わせ（レジストレーション）を行うように構成されてよい。例えば、データ処理部２３０は、ＯＣＴで取得された３次元画像データと、眼底カメラユニット２により取得された正面画像との間のレジストレーションを行うように構成される。また、データ処理部２３０は、ＯＣＴで取得された２つのＯＣＴ画像の間のレジストレーションを行うように構成される。また、データ処理部２３０は、眼底カメラユニット２により取得された２つの正面画像の間のレジストレーションを行うように構成される。また、ＯＣＴ画像の解析結果や、正面画像の解析結果に対してレジストレーションを適用するように構成されてもよい。これらのレジストレーションは、公知の手法によって実行可能であり、例えば特徴点抽出とアフィン変換とを含む。
更に、データ処理部２３０は、リサジュースキャンを用いて得られたデータセットを処理するように構成される。前述したように、画像データ構築部２２０は、例えば、リサジュースキャンで収集されたデータから複数のストリップを構築するように構成されている。データ処理部２３０は、これらストリップに対してレジストレーション及びマージ処理を適用することにより、モーションアーティファクトが補正された画像を構築するように構成されている。
リサジュースキャンの特性として、任意の２つのストリップはオーバーラップ領域を有する。非特許文献１に記載された手法と同様に、データ処理部２３０は、ストリップ間のオーバーラップを利用してストリップ間のレジストレーションを行うように構成される。データ処理部２３０は、まず、画像データ構築部２２０により構築された複数のストリップを大きさ（面積など）に従って順序付け（第１～第Ｎのストリップ）、最大のストリップである第１のストリップを初期基準ストリップに指定する。次に、データ処理部２３０は、第１のストリップを基準として第２のストリップのレジストレーションを行い、第１のストリップと第２のストリップとをマージする（合成する）。データ処理部２３０は、これにより得られたマージストリップを基準ストリップとして第３のストリップのレジストレーションを行い、このマージストリップと第３のストリップとをマージする。このようなレジストレーション及びマージ処理を上記順序に従って順次に実行することにより第１～第Ｎのストリップの位置合わせ及び貼り合わせがなされ、モーションアーティファクトが補正された画像が得られる。
このようなレジストレーションでは、基準ストリップと他のストリップとの間の相対位置を求めるために相互相関関数が利用されるが、ストリップは任意の形状を有するため、各ストリップを特定形状（例えば正方形状）の画像として扱うためのマスクを用いてストリップ間の相関演算が行われる（非特許文献１の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」（第１７８９頁～第１７９０頁）に相当するレジストレーションを実行するように構成されてもよい。
＜情報記録部２３４＞
まず、情報記録部２３４が実行する処理の前提について説明する。前述したように、相関係数算出部２３１５はストリップ間の相関係数の演算を実行し、ｚシフト量算出部２３２はサブボリューム間の相関係数の演算を実行している。ストリップ及びサブボリュームは、いずれも、サイクルに基づいて定義されている。よって、相関係数算出部２３１５により算出されるストリップ間の相関係数をサイクル間の相関係数とみなすことができ、ｚシフト量算出部２３２により算出されるサブボリューム間の相関係数をサイクル間の相関係数とみなすことができる。
このように、本開示において、サイクル間の相関係数は、ストリップ間の相関係数及び／又はサブボリューム間の相関係数を含む。すなわち、幾つかの例示的な態様はストリップ間の相関係数のみを扱ってもよく、幾つかの例示的な態様はサブボリューム間の相関係数のみを扱ってもよく、幾つかの例示的な態様はストリップ間の相関係数及びサブボリューム間の相関係数の双方を扱ってもよい。また、幾つかの例示的な態様は、ストリップ及びサブボリュームのどちらとも異なるオブジェクト間の相関係数を扱ってもよい。いずれの場合においても、本実施形態と同様の処理を実行することが可能である。
また、ストリップ間の相関係数は、２つのストリップの間の相関係数でもよいし、１つのストリップと２つ以上のストリップ（ストリップ群）との間の相関係数でもよいし、２つのストリップ群の間の相関係数であってもよい。ストリップ群は、例えば、２つ以上のストリップに基づくマージストリップ（前述）であってよい。本実施形態では、いずれの場合も相関係数算出部２３１５が演算を実行し、演算の内容は前述した処理と同様であってよい。
同様に、サブボリューム間の相関係数は、２つのサブボリュームの間の相関係数でもよいし、１つのサブボリュームと２つ以上のサブボリューム（サブボリューム群）との間の相関係数でもよいし、２つのサブボリューム群の間の相関係数であってもよい。いずれの場合もｚシフト量算出部２３２が演算を実行し、演算の内容は前述した処理と同様であってよい。ストリップ及びサブボリュームのいずれとも異なるオブジェクトの間の相関係数が扱われる場合も同様である。
前述したように、眼科検査装置１は、互いに交差する一連のサイクルを含む２次元パターンに従う光走査（例えば、リサジュースキャン）を被検眼Ｅに適用してデータセットを収集し（データセット収集部）、このデータセットに基づいてサイクル間の相関係数を算出する（相関係数算出部：相関係数算出部２３１５、ｚシフト量算出部２３２）。
本実施形態では、複数の相関係数が算出される。例えば、２つのストリップの間の相関係数が算出される場合、ストリップの個数（サイクルの個数）をＮ個とすると、Ｎ個のストリップから算出される相関係数の個数は、最多でＮ×（Ｎ－１）個になる（重複する組み合わせをまとめると、最多でＮ×（Ｎ－１）／２個になる）。この場合に算出される相関係数の個数は、ストリップ同士の組み合わせ方（つまり、考慮されるストリップペアの個数）により決定される。同様に、１つのストリップとストリップ群との間の相関係数が算出される場合や、２つのストリップ群の間の相関係数が算出される場合においても、複数の相関係数が算出される。これらの場合、算出される相関係数の個数は、考慮されるストリップ群の個数、ストリップとストリップ群との組み合わせ方、ストリップ群同士の組み合わせ方などにより決定される。例えば、相関係数算出部２３１５が実行する前述の逐次的な演算が適用される場合には、Ｎ－１個の相関係数が算出される。このように、収集されたデータセットから相関係数の集合（相関係数集合と呼ぶ）が生成される。
情報記録部２３４は、この相関係数集合に属する相関係数のうち所定の条件を満足する相関係数を要素とする部分集合を特定し、この部分集合（これに属する相関係数）に対応する走査位置を示す情報（走査位置情報と呼ぶ）を記録する。
走査位置情報の記録は、走査位置情報を所定の記憶装置に保存する処理を含み、更に、所定のデータ（例えば、画像データ処理部２３１により構築された画像データ）に走査位置情報を関連付ける処理を含んでいてもよい。走査位置情報（及び所定のデータ）が保存される記憶装置は、予め決められた任意の記憶装置であってよく、その例として、情報記録部２３４に含まれる記憶装置（図示省略）、記憶部２１２、眼科検査装置１によりアクセス可能な記憶装置（施設内データベースシステム、クラウド型データベースシステムなど）、眼科検査装置１によりデータ書き込み可能な記録媒体などがある。
走査位置情報に含まれる情報の範囲は任意であってよい。走査位置情報は、相関係数集合から特定された部分集合に対応する走査位置を示す情報を少なくとも含んでいればよい。例えば、幾つかの例示的な態様の走査位置情報は、相関係数集合から特定された部分集合の各要素に対応する走査位置を示す情報のみを含んでいてもよい。また、幾つかの例示的な態様の走査位置情報は、相関係数集合から特定された部分集合を含む集合（例えば、相関係数集合の全体又は一部）の各要素に対応する走査位置を示す情報を含んでいてもよい。
走査位置情報に含まれる情報の内容は任意であってよい。走査位置情報は、走査位置を示す情報を少なくとも含んでいればよい。幾つかの例示的な態様の走査位置情報は、走査位置を示す情報に加えて、次に説明する相関係数の判定条件や判定結果に関する情報を含んでいてもよい。例えば、幾つかの例示的な態様の走査位置情報は、判定条件の内容を示す情報、及び／又は、判定結果の内容を示す情報を含んでいてもよい。判定条件の内容を示す情報の例として、後述する第１の閾値、第２の閾値、特定範囲などがある。判定結果の内容を示す情報の例として、判定条件を満足するか否かの判定結果、判定条件の満足の程度、判定条件からの乖離の程度などがある。なお、走査位置情報に含まれる情報はこれらの例に限定されない。
相関係数集合の部分集合を特定するために、情報記録部２３４は、相関係数集合に属する各相関係数について、その相関係数が所定の条件を満足するか判定を行う。この判定条件は、予め決定された任意の条件であってよく、例えば、次の例示的な条件のいずれかであってよい。
判定条件の第１の例は、相関係数の値が小さいことである。この条件例が適用される場合、情報記録部２３４は、例えば、相関係数集合に含まれる各相関係数を予め設定された閾値（第１の閾値と呼ぶ）と比較することによって、この第１の閾値よりも値が小さい相関係数からなる部分集合を求め、この部分集合（低値の相関係数を要素とする集合）に対応する走査位置情報を記録するように構成される。この走査位置情報は、品質が低い画像の部分や画像化されなかった部分を示す情報であり、例えば、画像補完や再スキャンなどにこの走査位置情報を利用することができる。本例の走査位置情報の利用方法はこれらに限定されず、任意の処理（情報の提示、情報の提供、統計処理、画像解析、画像処理、画像評価、機械学習、機械学習用の訓練データの生成など）に利用してもよい。
判定条件の第２の例は、相関係数の値が大きいことである。この条件例が適用される場合、情報記録部２３４は、例えば、相関係数集合に含まれる各相関係数を予め設定された閾値（第２の閾値と呼ぶ）と比較することによって、この第２閾値よりも値が大きい相関係数からなる部分集合を求め、この部分集合（高値の相関係数を要素とする集合）に対応する走査位置情報を記録するように構成される。この走査位置情報は、品質が高い画像の部分を示す情報であり、任意の処理（情報の提示、情報の提供、統計処理、画像解析、画像処理、画像評価、機械学習、機械学習用の訓練データの生成など）にこの走査位置情報を利用することができる。
判定条件の第３の例は、予め設定された範囲に相関係数の値が含まれることである。この条件例が適用される場合、情報記録部２３４は、例えば、相関係数集合に含まれる各相関係数を予め設定された範囲（特定範囲と呼ぶ）と比較することによって、この特定範囲に値が属する相関係数からなる部分集合を求め、この部分集合に対応する走査位置情報を記録するように構成される。この走査位置情報は、品質が或る範囲に含まれる画像の部分を示す情報であり、任意の処理（情報の提示、情報の提供、統計処理、画像解析、画像処理、画像評価、機械学習、機械学習用の訓練データの生成など）にこの走査位置情報を利用することができる。
第１の閾値の設定方法、第２の閾値の設定方法、及び特定範囲（上限値、下限値）の設定方法は、いずれも任意であってよい。これらの値は、例えば、実測値及び／又はシミュレーション値に基づき決定されたデフォルト値であってもよいし、この相関係数演算が適用された画像データに基づき決定された個別値であってもよいし、この画像データを取得するための準備として取得された画像データ（例えば、アライメントなどの準備動作を行っているときに取得された画像データ）に基づき決定された個別値であってもよい。デフォルト値や個別値の決定は、典型的には、画像データに統計処理を適用することによって実行されてよい。例えば、画像データから生成される所定のパラメータ（相関係数、画素値など）の統計演算結果（ヒストグラム、マップ、統計量など）に基づいて、第１の閾値、第２の閾値、及び特定範囲のいずれか１つ以上を決定することができる。
判定条件は、走査位置情報の用途に応じて決定されてよい。また、準備される判定条件の個数は任意であってよい。２以上の判定条件が設けられる場合、これらの判定条件を選択的に使用できるようにしてもよい。また、２以上の判定条件を同じリサジュースキャンに係る処理において使用可能であってもよい。この場合、２以上の判定条件にそれぞれ対応する２以上の走査位置情報を生成することができる。
走査位置情報の表現方法は、予め設定された任意の表現方法であってよく、例えば、次の例示的な識別子のいずれかを用いて走査位置情報を表現することができる。
識別子の第１の例は、画像データ処理部２３１により構築された画像データの定義座標系の座標である。この識別子例が適用される場合、例えば、各画素について（換言すると、画素を単位として）走査位置情報が定義されてもよいし、予め設定された画像領域について（換言すると、画素群を単位として）走査位置情報が定義されてもよいし、これらの双方を利用して走査位置情報が定義されてもよい。
識別子の第２の例は、サイクルの順序付けに基づく番号（サイクル番号と呼ぶ）である。この識別子例が適用される場合、例えば、各サイクルについて（換言すると、サイクル（つまり、サイクル番号）を単位として）走査位置情報が定義されてもよいし、幾つかのサイクル（換言すると、サイクル群について（つまり、幾つかのサイクル番号）を単位として）走査位置情報が定義されてもよいし、これらの双方を利用して走査位置情報が定義されてもよい。
識別子の第３の例は、スキャンの順序に基づく番号（スキャン番号と呼ぶ）である。この識別子例が適用される場合、例えば、各サイクルのスキャンについて走査位置情報が定義されてもよいし、幾つかのサイクルのスキャンについて走査位置情報が定義されてもよいし、これらの双方を利用して走査位置情報が定義されてもよい。本実施形態では、サイクルごとにスキャンを実行しているので、本識別子例のスキャン番号は、第２の例のサイクル番号と実質的に同値である。
識別子の第４の例は、ストリップの順序付けに基づく番号（ストリップ番号と呼ぶ）である。この識別子例が適用される場合、例えば、各ストリップについて（換言すると、ストリップ（つまり、ストリップ番号）を単位として）走査位置情報が定義されてもよいし、幾つかのストリップ（換言すると、ストリップ群について（つまり、幾つかのストリップ番号）を単位として）走査位置情報が定義されてもよいし、これらの双方を利用して走査位置情報が定義されてもよい。前述した画像データ処理部２３１が実行する処理においてサイクルとストリップとは一対一に対応しているので、本識別子例のストリップ番号は、第２の例のサイクル番号及び第３の例のスキャン番号のそれぞれと実質的に同値である。
識別子の第５の例は、サブボリュームの順序付けに基づく番号（サブボリューム番号と呼ぶ）である。この識別子例が適用される場合、例えば、各サブボリュームについて（換言すると、サブボリューム（つまり、サブボリューム番号）を単位として）走査位置情報が定義されてもよいし、幾つかのサブボリューム（換言すると、サブボリューム群について（つまり、幾つかのサブボリューム番号）を単位として）走査位置情報が定義されてもよいし、これらの双方を利用して走査位置情報が定義されてもよい。前述した画像データ処理部２３１が実行する処理においてストリップとサブボリュームとは一対一に対応しているので、本識別子例のサブボリューム番号は、第２の例のサイクル番号、第３の例のスキャン番号、及び第４の例のストリップ番号のそれぞれと実質的に同値である。
識別子の第６の例は、マージ処理の適用順序に基づく番号（マージ処理番号と呼ぶ）である。前述した画像データ処理部２３１が実行する処理においてマージ処理は各ストリップ又は各サブボリュームについて実行されるので、本識別子例のマージ処理番号は、第２の例のサイクル番号、第３の例のスキャン番号、第４の例のストリップ番号、及び第５の例のサブボリューム番号のそれぞれと実質的に同値である。
情報記録部２３４は、走査位置情報に加え、画像データ補正部２３３により作成された画像データ（モーションアーティファクトが補正された３次元画像データ）を記録することができる。更に、情報記録部２３４は、この走査位置情報とこの画像データとを互いに関連付けて記録することができる。これにより、画像データの処理（例えば、補完、補正、調整など）を走査位置情報に基づいて実行することや、画像データを処理するための動作（再スキャンなど）を走査位置情報に基づいて実行することや、画像データから生成された画像と走査位置情報から生成された情報とを互いに関連付けて表示することや、画像データ（及び／又は画像データから生成された画像）と走査位置情報（及び／又は走査位置情報から生成された情報）とに基づく処理（解析処理、統計処理、画像処理、機械学習、機械学習用の訓練データの生成など）を実行することなどが可能になる。
＜ユーザーインターフェイス２４０＞
ユーザーインターフェイス２４０は表示部２４１と操作部２４２とを含む。表示部２４１は表示装置３を含む。操作部２４２は各種の操作デバイスや入力デバイスを含む。ユーザーインターフェイス２４０は、例えばタッチパネルのような表示機能と操作機能とが一体となったデバイスを含んでいてもよい。ユーザーインターフェイス２４０の少なくとも一部を含まない実施形態を構築することも可能である。例えば、表示デバイスは、眼科検査装置１に接続された外部装置であってよい。
＜データ受付部２５０＞
データ受付部２５０は、外部装置からデータを取得する。この外部装置は、例えば、コンピュータ、記憶装置、記録媒体、情報システム（例えば、病院情報システム、電子カルテシステム、画像アーカイビングシステム）などであってよい。また、外部装置は、例えば、眼科検査装置１に直接に接続された装置、ローカルエリアネットワーク（ＬＡＮ）を介して眼科検査装置１に接続された装置、ワイドエリアネットワーク（ＷＡＮ）を介して眼科検査装置１に接続された装置などであってよい。
データ受付部２５０は、例えば、通信インターフェイス、ドライブ装置などを含んでいてもよい。また、データ受付部２５０は、紙葉類に記録されたデータを読み取るためのスキャナを含んでいてもよい。また、データ受付部２５０は、ユーザーがデータを入力するための入力デバイスを含んでいてもよい。入力デバイスは、キーボード、ペンタブレットなどであってよい。
＜動作＞
眼科検査装置１の動作について説明する。眼科検査装置１の動作の例を図６Ａ～図６Ｆに示す。この動作例では、リサジュースキャンにおける画像データ構築（モーションアーティファクト補正を内在している）について詳細に説明する。幾つかの例示的な態様では、この動作例の一部の工程のみを実行してもよい。また、幾つかの例示的な態様では、この動作例の一部の工程を他の工程（例えば、類似の工程）に置換してもよい。また、幾つかの例示的な態様では、この動作例に任意の工程を組み合わせてもよい。
患者ＩＤの入力、走査モードの設定（リサジュースキャンの指定）、固視標の提示、アライメント、フォーカス調整、ＯＣＴ光路長調整など、従来と同様の準備動作が、ステップＳ１の前に実行される。
幾つかの例示的な態様において、眼科検査装置１は、この準備動作において取得された画像データに基づいて、相関係数集合の部分集合を特定するための判定条件（例えば、前述の第１の閾値、第２の閾値、及び特定範囲の少なくとも１つ）を決定するように構成されていてよい。なお、デフォルトの判定条件を用いてもよいし、ステップＳ１のリサジュースキャンで収集されたデータセットから決定された判定条件を用いてもよい。
（Ｓ１：リサジュースキャン）
所定の走査開始トリガー信号を受けて、走査制御部２１１１は、被検眼Ｅ（眼底Ｅｆ）に対するＯＣＴスキャン（リサジュースキャン）の適用を開始する。
走査開始トリガー信号は、例えば、所定の準備動作（アライメント、フォーカス調整、ＯＣＴ光路長調整など）が完了したことに対応して、又は、操作部２４２を用いて走査開始指示操作が行われたことに対応して生成される。
走査制御部２１１１は、走査プロトコル２１２１（リサジュースキャンに対応したプロトコル）に基づいて光スキャナ４４及びＯＣＴユニット１００などを制御することによってリサジュースキャンを眼底Ｅｆに適用する。リサジュースキャンにより収集されたデータは画像データ構築部２２０に送られる。
（Ｓ２：複数のストリップを構築する）
画像データ構築部２２０は、前述した要領で、ステップＳ１で収集されたデータに基づいて複数のストリップを構築する。構築された複数のストリップはデータ処理部２３０に送られる。
（Ｓ３：基準ストリップと対象ストリップを設定する）
データ処理部２３０は、ステップＳ２で構築された複数のストリップを寸法（面積など）に従って順序付けする。本例では、ステップＳ２で構築された複数のストリップの個数をＮ個（Ｎは２以上の整数）であるとする。また、順序付けにより指定された順序に従って、これらＮ個のストリップを、第１のストリップ、第２のストリップ、・・・、第Ｎのストリップと呼ぶ。また、第１～第Ｎのストリップのうちの任意のストリップを第ｎのストリップと呼ぶことがある（ｎ＝１、２、・・・、Ｎ）。このように、第１のストリップは最も大きなストリップであり、第２のストリップは２番目に大きなストリップであり、第ｎのストリップはｎ番目に大きなストリップであり、第Ｎのストリップは最も小さなストリップである。
本例において、データ処理部２３０は、第１のストリップを初期基準ストリップに設定し、第２のストリップを初期対象ストリップに設定する。
なお、基準ストリップ及び対象ストリップは、それぞれ、非特許文献１の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）のみにする。
次に、ｆ´（ｒ）ｍｆ（ｒ）、ｇ´（ｒ）ｍｇ（ｒ）、ｍｆ（ｒ）、ｍｇ（ｒ）、（ｆ´（ｒ）ｍｆ（ｒ））２、及び（ｇ´（ｒ）ｍｇ（ｒ））２のそれぞれの実部に対して高速フーリエ変換（ＦＦＴ）を適用する。
次に、これらの高速フーリエ変換で得られた関数群に基づいて、図６ＡのステップＳ７に示された６個の相互相関関数を導出する。
次に、導出された６個の相互相関関数のそれぞれに逆高速フーリエ変換（ＩＦＦＴ）を適用する。
次に、非特許文献１の式（３３）の相関係数ρ（ｒ´）を算出する。
そして、相関係数ρ（ｒ´）のピーク位置を特定することによって基準ストリップｆ（ｒ）と対象ストリップｇ（ｒ）との間の相対的なシフト量（Δｘ、Δｙ）を算出し、このシフト量に基づいて基準ストリップｆ（ｒ）と対象ストリップｇ（ｒ）との間のレジストレーション及びマージ処理を実行する。
このような一連の処理をＮ個のストリップに対して逐次に適用することにより、Ｎ個のストリップの全てのマージ画像が得られる。また、当該ストリップペアについて求められたｘｙ方向シフト量（Δｘ、Δｙ）を時間情報とともに（つまり、時間情報に関連付けて）保存することができる。
（Ｓ１４：全てのストリップを処理したか？）
ステップＳ２で得られたＮ個のストリップの全てに対し、前述した順序に従って逐次に、ステップＳ３～Ｓ１３の一連の処理が実行される。
Ｎが３以上の場合において、ステップＳ１３で第１のストリップと第２のストリップとのマージ画像が作成された場合（Ｓ１４：Ｎｏ）、処理はステップＳ３に戻る。ステップＳ３では、第１のストリップと第２のストリップとのマージ画像が新たな基準ストリップに設定され、且つ、第３のストリップが新たな対象ストリップに設定される。新たな基準ストリップと新たな対象ストリップに基づきステップＳ４～Ｓ１３の処理を実行することにより、新たな基準ストリップと新たな対象ストリップとのマージ画像が得られる。この新たなマージ画像は、第１～第３のストリップのマージ画像である。このような一連の処理をＮ個のストリップに対して逐次に適用することにより、Ｎ個のストリップの全てのマージ画像が得られる（Ｓ１４：Ｙｅｓ）。
（Ｓ１５：ｘｙ方向のモーションアーティファクトが補正された画像を保存する）
以上に説明した繰り返し処理により最終的に得られたマージ画像は、ｘｙ方向のモーションアーティファクトが補正された画像であり、且つ、ステップＳ１のリサジュースキャンの適用範囲全体を表現した画像である。主制御部２１１は、この最終的なマージ画像を記憶部２１２（及び／又は、他の記憶装置）に保存することができる。
データ処理部２３０は、Ｎ個のストリップに基づくレジストレーションの結果を利用して、ステップＳ１で収集されたデータ（３次元データ）、及び／又は、この３次元データから画像データ構築部２２０が構築した３次元画像データのレジストレーションを行うことができる。つまり、データ処理部２３０は、最終的に得られたマージ画像を利用して３次元データや３次元画像データのレジストレーションを行うことができる。このレジストレーションは、リサジュースキャンの定義座標系を３次元直交座標系（ｘｙｚ座標系）に変換する処理を含む。この座標変換は、非特許文献１に記載された手法における「remapping」に相当する。このようにして、ｘｙ方向のモーションアーティファクトが補正された３次元画像データ（ボリューム）が得られる。主制御部２１１は、このｘｙ方向モーションアーティファクト補正済みボリュームを、最終的なマージ画像とともに又はその代わりに、記憶部２１２（及び／又は、他の記憶装置）に保存することができる。
（Ｓ１６：基準ストリップを選択する）
ｚ方向のモーションアーティファクトを補正するための処理に移行する。ｚ方向モーションアーティファクト補正では、まず、ｚシフト量算出部２３２が、ｘｙ方向モーションアーティファクトの補正がなされたボリュームに対応する複数のストリップから、１つのストリップ（基準ストリップ）を選択する。ここで選択される基準ストリップは、第１回目のステップＳ３で選択された第１のストリップであってもよいし、これとは異なるストリップであってもよい。
（Ｓ１７：対象ストリップを選択する）
次に、ｚシフト量算出部２３２は、ステップＳ１６で選択された基準ストリップとは別の１つのストリップ（対象ストリップ）を、ｘｙ方向モーションアーティファクトが補正されたボリュームに対応する複数のストリップからを選択する。ここで選択される基準ストリップは、第１回目のステップＳ３で選択された第２のストリップであってもよいし、これ以外のストリップであってもよい。
（Ｓ１８：基準ストリップに対応する基準サブボリュームを設定する）
次に、ｚシフト量算出部２３２は、ステップＳ１６で選択された基準ストリップに対応するサブボリュームを基準サブボリュームに設定する。ここで、基準サブボリュームの正面プロジェクション画像が基準ストリップである。
（Ｓ１９：対象ストリップに対応する対象サブボリュームを設定する）
同様に、ｚシフト量算出部２３２は、ステップＳ１７で選択された対象ストリップに対応するサブボリュームを対象サブボリュームに設定する。ここで、対象サブボリュームの正面プロジェクション画像が対象ストリップである。
（Ｓ２０：基準サブボリュームと対象サブボリュームとの交差領域を特定する）
次に、ｚシフト量算出部２３２は、ステップＳ１８で設定された基準サブボリュームと、ステップＳ１９で設定された対象サブボリュームとの交差領域（共通領域）を特定する。
前述したように、２つのサブボリュームには４つの交差領域がある。ｚシフト量算出部２３２は、１つ又は２つ以上の交差領域を特定する。２つ以上の交差領域を特定する場合、特定された各交差領域について次のステップＳ２１が実行される。
（Ｓ２１：基準サブボリュームと対象サブボリュームとの間のｚ方向シフト量を算出する）
次に、ｚシフト量算出部２３２は、ステップＳ２０で特定されたサブボリューム間の交差領域に基づいて、基準サブボリュームと対象サブボリュームとの間のｚ方向シフト量を算出する。
例えば、ｚシフト量算出部２３２は、ステップＳ２０で特定された交差領域に断面（注目断面）を設定する。注目断面は、交差領域の任意の断面であってよい。次に、ｚシフト量算出部２３２は、基準サブボリュームから注目断面の画像（基準断面画像）を構築し、且つ、対象サブボリュームから注目断面の画像（対象断面画像）を構築する。次に、ｚシフト量算出部２３２は、基準断面画像を解析して被検眼Ｅの所定部位の像（基準像）を特定し、且つ、対象断面画像を解析して同部位の像（対象像）を特定する。次に、ｚシフト量算出部２３２は、基準像のｚ座標と対象像のｚ座標とを求め、これら２つのｚ座標の差を算出する。本例では、算出された差が、基準サブボリュームと対象サブボリュームとの間のｚ方向シフト量に採用される。
図８に具体例を示す。符号４０１は、基準サブボリュームに対応するストリップ４０２と対象サブボリュームに対応するストリップ４０３との交差状態を表している。本例では、白丸で囲まれた交差領域に着目する。ｘ方向に沿った白丸の直径（点線で示されている）が注目断面である。２つのストリップ４０２及び４０３にもそれぞれ交差領域と注目断面とが同様に示されている。ｚシフト量算出部２３２は、ストリップ４０２に対応する基準サブボリュームから、この注目断面の画像（基準断面画像）４１２を構築し、且つ、ストリップ４０３に対応する対象サブボリュームから、この注目断面の画像（対象断面画像）４１３を構築する。更に、ｚシフト量算出部２３２は、基準断面画像４１２を解析して網膜表面の像のｚ位置（ｚ座標）４２２を特定し、且つ、対象断面画像４１３を解析して網膜表面の像のｚ位置（ｚ座標）４２３を特定する。加えて、ｚシフト量算出部２３２は、基準サブボリュームと対象サブボリュームとの間のｚ方向シフト量として、２つの断面画像４１２及び４１３からそれぞれ特定された２つのｚ位置（２つのｚ座標）４２２及び４２３の間の差Δｚを求める。
基準サブボリュームと対象サブボリュームとの間のｚ方向シフト量を算出する手法は、これに限定されず、それらの交差領域を利用する任意の手法であってよい。例えば、前述したように、ｚシフト量算出部２３２は、基準サブボリュームと対象サブボリュームとの間の相関係数を算出し、この相関係数に基づいてｚ方向シフト量を求めてもよい。
この場合、図６ＢのステップＳ９と同様に、各サブボリュームペア算出された相関係数は、情報記録部２３４によって記録される。例えば、情報記録部２３４は、ｚ方向シフト量の算出において求められた相関係数と、対応する識別子とを互いに関連付けて記録する。この識別子は、前述したように、画像データの定義座標系の座標、サイクル番号、スキャン番号、ストリップ番号、サブボリューム番号、マージ処理番号、及びこれらと異なる識別子のいずれかであってよい。ステップＳ１７～Ｓ２４の繰り返しにより、このリサジュースキャンの相関係数集合に含まれる相関係数が蓄積されていく。ここで、ステップＳ９で記録されるｘｙ方向の相関係数と、ステップＳ２１で記録されるｚ方向の相関係数とは、識別子を介して互いに関連付けられる。これにより、３次元的な相関係数に対応した走査位置情報を生成することが可能になる。
（Ｓ２２：ｚ方向シフト量を保存する）
ステップＳ２１において当該サブボリュームペアについて算出されたｚ方向シフト量は、例えば記憶部２１２に保存される。ｘｙ方向シフト量の保存（ステップＳ１１）と同様に、ｚ方向シフト量は、例えば、当該サブボリュームペアに対応する時間情報とともに記録されてもよい。
（Ｓ２３：基準サブボリュームと対象サブボリュームとの間のレジストレーションを行う）
画像データ補正部２３３は、ステップＳ２１で算出されたｚ方向シフト量に基づいて、基準サブボリューム及び対象サブボリュームに対してｚ方向のレジストレーションを適用する。
このレジストレーションは、例えば、図８のｚ方向シフト量Δｚを打ち消すように、つまり、基準断面画像４１２中の網膜表面の像と対象断面画像４１３中の網膜表面の像とが同じｚ位置（等しいｚ座標）に配置されるように、基準サブボリュームと対象サブボリュームとのｚ方向の相対位置を調整する処理である。
（Ｓ２４：基準サブボリュームと対象サブボリュームとのマージ画像を構築する）
画像データ補正部２３３は、ステップＳ２３のレジストレーションによりｚ方向の相対位置調整がなされた基準サブボリュームと対象サブボリュームとのマージ画像を構築する。これにより、ｘｙ方向及びｚ方向について相対位置調整がなされた基準サブボリュームと対象サブボリュームとのマージ画像が得られる。
（Ｓ２５：全てのサブボリュームを処理したか？）
全てのサブボリュームに対し、所定の順序（例えば、第１～第Ｎのストリップに割り当てられた順序）に従って逐次に、ステップＳ１７～Ｓ２４の一連の処理が実行される（Ｓ２５：Ｎｏ）。なお、本例では、ステップＳ２４で構築されたマージ画像が、次のルーチンにおける基準ストリップとして設定される。このような一連の処理を全てのサブボリュームに対して逐次に適用することにより、全てのサブボリュームのマージ画像が得られる（Ｓ２５：Ｙｅｓ）。
（Ｓ２６：全てのサブボリュームがマージされたボリュームを保存する）
このようなステップＳ１７～Ｓ２４の反復により最終的に得られたマージ画像は、ｘｙ方向及びｚ方向のモーションアーティファクトが補正された画像であり、且つ、ステップＳ１のリサジュースキャンの適用範囲全体を表現した画像である。主制御部２１１は、この最終的なマージ画像を記憶部２１２（及び／又は、他の記憶装置）に保存することができる。
以下の処理は、ｚ方向のモーションアーティファクト補正をより精密に行うためのものである。ステップＳ１６～Ｓ２５では、ｚ方向のモーションアーティファクト補正をサブボリューム単位（ストリップ単位）で行っている。これに対し、以下に説明するステップＳ２７～Ｓ３２では、ｚ方向のモーションアーティファクト補正をサイクル単位で実行する。このサイクル単位の補正は、データ処理部２３０（例えば、画像データ補正部２３３）によって実行される。
なお、精密なｚ方向モーションアーティファクト補正は、サイクル単位で実行される必要はなく、サブボリュームの一部を構成するサイクル群を単位として実行するようにしてもよい。また、精密なｚ方向モーションアーティファクト補正は、全てのサイクル（全てのサイクル群）について実行される必要はなく、一部のサイクル（一部のサイクル群）についてのみ実行するようにしてもよい。
（Ｓ２７：マージされたボリュームからサイクルを選択する）
画像データ補正部２３３は、全てのサブボリュームがマージされたボリュームから、１つのサイクル（に対応する部分データ）を選択する。
最初に選択されるサイクルは、任意であってよい。例えば、ステップＳ１のリサジュースキャンにおける複数のサイクルのスキャン順（時系列順、スキャン番号、サイクル番号）に従ってサイクルが順次に選択される。或いは、複数のストリップ又は複数のサブボリュームに割り当てられた順序（ストリップ番号、サブボリューム番号）などに基づいてサイクルの選択順序が設定されてもよい。
（Ｓ２８：サイクルをボリュームから抜去する）
次に、画像データ補正部２３３は、ステップＳ２７で選択されたサイクルに対応する部分データをボリュームから抜き出す。
（Ｓ２９：サイクル抜去後のボリュームとサイクルとの間のｚ方向シフト量を算出する）
次に、画像データ補正部２３３は、ステップＳ２８で抜き出されたサイクルに対応する部分データと、当該サイクルが抜き出されたボリュームとの間におけるｚ方向シフト量を算出する。
本ステップのｚ方向シフト量算出の手法は、ステップＳ２１の手法と同様であってよい。例えば、画像データ補正部２３３は、ステップＳ２８でボリュームから抜去されたサイクルの部分データ（注目断面画像）を解析して被検眼Ｅの所定部位の像（第１の像）を特定し、且つ、注目断面画像が抜き出されたボリュームを解析して同部位の像（第２の像）を特定する。第２の像を特定するための解析は、注目断面画像が抜き出されたボリュームの全体に対して適用される必要はなく、例えば、注目断面画像に隣接する部分（隣接するサイクル）のみに対して適用されてよい。更に、画像データ補正部２３３は、第１の像のｚ座標と第２の像のｚ座標とを求め、これら２つのｚ座標の差を算出する。本例では、算出された差が、当該サイクルと当該サイクル抜去後のボリュームとの間のｚ方向シフト量に採用される。
（Ｓ３０：ｚ方向シフト量を更新する）
次に、主制御部２１１は、ステップＳ２１で求められたｚ方向シフト量を、新たにステップＳ２９で算出されたｚ方向シフト量に置き換える。ここでは、ステップＳ２８で抜去されたサイクルに対応する部分データを含むサブボリュームのｚ方向シフト量が更新される。
ステップＳ２１のｚ方向シフト量の算出とステップＳ２９のｚ方向シフト量の算出との双方において相関係数が用いられた場合、ステップＳ２１で記録された情報（相関係数、識別子など）をステップＳ２９で取得された情報（相関係数、識別子など）で置き換えることができ、又は、ステップＳ２１で記録された情報とステップＳ２９で取得された情報との双方を互いに関連付けて記録することができる。ここで、ステップＳ２１で記録された情報とステップＳ２９で取得された情報との間の対応付けは、識別子を介して行うことができる。ステップＳ２１のｚ方向シフト量の算出及びステップＳ２９のｚ方向シフト量の算出のいずれか一方のみにおいて相関係数が用いられた場合には、その一方の処理で得られた情報を記録することができる。
（Ｓ３１：レジストレーション及びマージングを行う）
次に、画像データ補正部２３３は、ステップＳ２９で求められたｚ方向シフト量を打ち消すように、ステップＳ２８で抜去されたサイクルに対応する部分データと、サイクル抜去後のボリュームとの間のレジストレーションを実行する。更に、画像データ補正部２３３は、このレジストレーションによりｚ方向の相対位置調整がなされた当該部分データとサイクル抜去後のボリュームとのマージ画像を構築する。これにより、ｘｙ方向について相対位置調整がなされ且つｚ方向について精密に相対位置調整がなされたボリュームが得られる。
（Ｓ３２：全てのサイクルを処理したか？）
全てのサイクルに対し、所定の順序に従って逐次に、ステップＳ２７～Ｓ３１の一連の処理が実行される（Ｓ３２：Ｎｏ）。このような一連の処理を全てのサイクルに対して逐次に適用することにより、全サイクル間のｚ方向相対位置調整がなされたボリュームが得られる（Ｓ３２：Ｙｅｓ）。主制御部２１１は、この最終的なボリュームを記憶部２１２（及び／又は、他の記憶装置）に保存することができる。
（Ｓ３３：相関係数集合から所定の条件を満足する相関係数を特定する）
本動作例では、ステップＳ９を含むルーチン（並びに、ステップＳ２１を含むルーチン及びステップＳ２９を含むルーチンの少なくとも一方）を繰り返し行うことにより、相関係数や識別子などの情報が蓄積される。蓄積された情報は、前述した相関係数集合に相当する。
情報記録部２３４は、この相関係数集合に含まれている相関係数のうち、所定の条件を満足する相関係数を特定する。この所定の条件は、前述した判定条件に相当する。これにより、前述した部分集合が得られる。
（Ｓ３４：走査位置情報を生成する）
情報記録部２３４は、ステップＳ３３で特定された部分集合に対応する走査位置情報を生成する。例えば、情報記録部２３４は、ステップＳ３３で特定された部分集合に属する各相関係数に関連付けられた識別子から走査位置情報を生成する。前述したように、識別子は、ステップＳ１のリサジュースキャンが適用された位置（走査位置）を示す情報（すなわち、走査位置の表現）であり、例えば、画像データの定義座標系の座標、サイクル番号、スキャン番号、ストリップ番号、サブボリューム番号、マージ処理番号、及びこれらと異なる識別子のうちの少なくとも１つを含んでいてよい。
このようにして生成された走査位置情報は、判定条件を満足する相関係数に対応する走査位置を示す情報を含んでいる。
判定条件が前述した第１の例（相関係数の値が小さいこと）である場合、走査位置情報は、第１の閾値よりも値が小さい相関係数に対応する走査位置を示す情報、つまり、交差領域における相関が低いサイクルペアの位置を示す情報を含む。
判定条件が前述した第２の例（相関係数の値が大きいこと）である場合、走査位置情報は、第２の閾値よりも値が大きい相関係数に対応する走査位置を示す情報、つまり、交差領域における相関が高いサイクルペアの位置を示す情報を含む。
判定条件が前述した第３の例（特定範囲に相関係数の値が含まれること）である場合、走査位置情報は、この特定範囲に値が属する相関係数に対応する走査位置を示す情報、つまり、交差領域における相関が特定範囲に含まれるサイクルペアの位置を示す情報を含む。
これらの例と異なる判定条件が用いられる場合においても同様に、走査位置情報は、適用される判定条件を満足する相関係数に対応する走査位置を示す情報、つまり、交差領域における相関がこの判定条件に合致するサイクルペアの位置を示す情報を含む。
（Ｓ３５：走査位置情報を記録する）
情報記録部２３４は、ステップＳ３４で生成された走査位置情報を記録する（エンド）。情報記録部２３４は、この走査位置情報と、ステップＳ３２で作成された最終的なボリューム（モーションアーティファクトが補正された３次元画像データ）とを互いに関連付けて記録してもよい。また、情報記録部２３４は、本動作例において生成された任意の情報を走査位置情報に関連付けて記録してもよい。走査位置情報とともに記録される情報の種類は、予め設定されていてもよいし、ユーザー又は眼科検査装置１によって指定されてもよい。
記録された走査位置情報は任意の処理に利用される。この処理は、例えば、情報の提示（走査位置情報の表示、走査位置情報から生成された情報の表示など）、情報の提供（走査位置情報の送信、走査位置情報から生成された情報の送信など）、統計処理、画像解析、画像処理、画像評価、機械学習、及び、機械学習用の訓練データの生成のうちのいずれか１つ以上であってよい。これらの処理において、走査位置情報とともに記録された情報を用いることも可能である。
＜効果＞
第１の実施形態に係る眼科検査装置１の幾つかの効果について説明する。
眼科検査装置１は、光走査を用いて画像データを構築する走査型イメージング装置であって、その機能要素としてデータセット収集部と相関係数算出部と情報記録部とを含んでいる。
データセット収集部は、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータセットを収集するように機能する。本実施形態では、ＯＣＴデータを収集するための要素群がデータセット収集部に相当し、例えば、眼底カメラユニット２、ＯＣＴユニット１００、画像データ構築部２２０などがデータセット収集部に相当している。
相関係数算出部は、データセット収集部により収集されたデータセットに基づいてサイクル間の相関係数を算出するように機能する。本実施形態では、相関係数算出部２３１５、ｚシフト量算出部２３２などが相関係数算出部に相当している。
情報記録部は、相関係数算出部によりデータセットから算出された相関係数集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録するように機能する。本実施形態では、情報記録部２３４などが情報記録部に相当している。
このような眼科検査装置１によれば、従来の技術では記録されていなかった走査位置情報を記録することができ、しかも、相関係数に応じた走査位置情報を記録することができるので、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。この利用法については、その幾つかの非限定的な例が後述の実施形態において説明される。
本実施形態の幾つかの態様において、上記部分集合の要素は、所定の閾値以下の相関係数であってよい。すなわち、本態様の情報記録部は、相関係数算出部によりデータセットから算出された相関係数集合のうち所定の閾値以下の相関係数を要素とする部分集合に対応する走査位置情報を記録するように構成されていてもよい。本態様によれば、光走査を用いて構築された画像中の低品質の部分や画像化されなかった部分を示す走査位置情報を生成し記録することができる。この走査位置情報を利用することで、例えば、画像補完や再スキャンなどを行うことができる。
第１の実施形態の幾つかの応用例について以下の実施形態で説明する。以下の実施形態では、特に言及しない限り、第１の実施形態と同様の事項（同様の構成、同様の要素、同様の機能、同様の作用、同様の効果など）を採用することができる。なお、第１の実施形態と均等な実施形態又は第１の実施形態に類似した実施形態についても、以下の実施形態と同様の応用を行うことが可能である。
＜第２の実施形態＞
第２の実施形態について説明する。本実施形態の眼科検査装置は、第１の実施形態の図１～図４Ａと同様の構成を備えているとともに、図４Ｂの画像データ処理部２３１の代わりに図９の画像データ処理部２３１Ａを備えている。本実施形態の画像データ処理部２３１Ａは、第１の実施形態の画像データ処理部２３１に画素値決定部２３５を付加したものである。
本実施形態の情報記録部２３４は、リサジュースキャンなどの光走査で収集されたデータセットから算出された相関係数集合から、所定の閾値以下の相関係数を要素とする部分集合を特定し、この部分集合に対応する走査位置情報を記録するように構成されている。本実施形態の眼科検査装置は、画像データ処理部２３１Ａにより構築された画像データ（例えば、画像データ補正部２３３により生成された最終的な３次元画像データ、又は、マージ処理部２３１８により生成されたマージ画像）中の低品質の部分をより高品質な画素に置き換える処理、及び／又は、この画像データ中の画像化されなかった部分を補完する処理を実行する機能要素を備えている。この機能要素の例が画素値決定部２３５である。
画素値決定部２３５は、情報記録部２３４により記録された走査位置情報に対応する画素の値を１以上の近傍画素の値に基づいて決定するように構成されている。画素値決定部２３５は、プロセッサを含むハードウェアと、画素値決定処理ソフトウェアとの協働によって実現される。
画素値決定部２３５が実行する画素値決定処理の方法は任意であってよい。例えば、画素値決定部２３５は、任意の補間（内挿）及び／又は任意の外挿を実行可能であってよい。補間法の例として、最近傍補間、線形補間、バイリニア補間、トライリニア補間、バイキュービック補間、曲線あてはめ、多項式補間などがある。外挿法の例として、線形外挿、リチャードソン補外などがある。
例えば、第１の実施形態の図６ＦのステップＳ３５の後に、画素値決定部２３５は、ステップＳ３５で記録された走査位置情報に対応する画素の値を１以上の近傍画素の値に基づいて決定することができる。この走査位置情報は、第１の実施形態で説明した第１の閾値よりも値が小さい相関係数に対応する走査位置を示す情報を含んでいる。換言すると、この走査位置情報は、交差領域における相関が低いサイクルペアの位置を示す情報を含んでいる。画素値決定部２３５は、このような走査位置（例えば、当該サイクルペア又は当該交差領域）に存在する各画素の画素値を決定するように機能するものである。
本実施形態に係る眼科検査装置の動作について説明する。眼科検査装置の動作の例を図１０に示す。特に言及しない限り、第１の実施形態で説明された事項をこの動作例に適用することができる。
（Ｓ４１：サイクル個数を指定する）
患者ＩＤの入力や走査モードの設定（リサジュースキャンの指定）の後、任意的な工程として、被検眼Ｅに適用されるリサジュースキャンにおけるサイクルの個数を指定することができる。サイクルの個数は、スキャンの密度に相当する。サイクルの個数の指定は、例えば、ユーザー又は眼科検査装置によって行われる。
（Ｓ４２：準備動作を行う）
眼科検査装置は、固視標の提示、アライメント、フォーカス調整、ＯＣＴ光路長調整など、従来と同様の準備動作を実行する。眼科検査装置は、この準備動作において取得された画像データに基づいて、相関係数集合の部分集合を特定するための判定条件（例えば、前述の第１の閾値）を決定するように構成されていてよい。なお、デフォルトの判定条件を用いてもよいし、ステップＳ４３のリサジュースキャンで収集されたデータセットから決定された判定条件を用いてもよい。
（Ｓ４３：リサジュースキャンを行う）
所定の走査開始トリガー信号を受けて、走査制御部２１１１は、ステップＳ４１で指定されたサイクル個数からなるリサジュースキャンを被検眼Ｅ（眼底Ｅｆ）に適用する。
（Ｓ４４：各サイクルについてストリップを構築する）
画像データ構築部２２０は、ステップＳ４３のリサジュースキャンで収集されたデータセットに基づいて、ステップＳ４１で指定された各サイクルに対応するストリップを構築する。この処理は、例えば、第１の実施形態（図６ＡのステップＳ２）と同じ要領で実行されてよい。
（Ｓ４５：各サイクルペアについて相関係数を算出する）
画像データ処理部２３１Ａ（相関係数算出部２３１５など）は、ステップＳ４１で指定された複数のサイクルに基づき形成された各サイクルペアについての相関係数を算出する。
本ステップの処理は、例えば、第１の実施形態のモーションアーティファクト補正（図６Ａ及び図６ＢのステップＳ３～ステップＳ８）と同じ要領で実行されてよい。
或いは、本ステップの処理は、ステップＳ４１で指定された複数のサイクルに基づき形成された各サイクルペアについて、このサイクルペアを構成する２つのサイクルの間の相関係数を算出するようにしてもよい。この場合、第１の実施形態のモーションアーティファクト補正を別途に実行してもよい。
（Ｓ４６：相関係数が閾値以下のストリップの画素値をゼロに設定する）
画像データ処理部２３１Ａ（例えば、情報記録部２３４）は、ステップＳ４５で算出された各相関係数を前述の第１の閾値と比較することで、第１の閾値よりも値が小さい相関係数を特定する。これにより、ステップＳ４５で生成された相関係数集合における、第１の閾値よりも値が小さい相関係数からなる部分集合が特定される。
更に、画像データ処理部２３１Ａ（例えば、画像データ補正部２３３又は情報記録部２３４）は、この部分集合の各要素（低値の相関係数）に対応する走査位置の画素の画素値をゼロに設定する。幾つかの例示的な態様では、少なくとも１つの低値の相関係数に対応するストリップの各画素の画素値をゼロに設定してもよい。幾つかの例示的な態様では、対応する相関係数の全てが低値であるストリップの各画素の画素値をゼロに設定してもよい。幾つかの例示的な態様では、対応する低値の相関係数の個数が所定数以上であるストリップの各画素の画素値をゼロに設定してもよい。これらの例では、画素値がゼロに設定されたストリップは黒の塗りつぶしとして表現され、第１の実施形態の「画像化されなかった部分」に相当する。なお、本ステップの処理はこれらの例に限定されない。
（Ｓ４７：画素値をゼロにした画素の位置（走査位置情報）を記録する）
画像データ処理部２３１Ａ（情報記録部２３４）は、ステップＳ４６で画素値がゼロに設定された画素の位置を走査位置情報として記録する。
（Ｓ４８：画素値ゼロの画素の補完を行う）
画像データ処理部２３１Ａ（画素値決定部２３５）は、ステップＳ４７で走査位置情報として記録された画素（すなわち、ステップＳ４６で画素値がゼロに設定された画素）の補完を行う。この補完処理の対象となる画素は、ステップＳ４７で走査位置情報として記録された全ての画素のうちの１つ以上であってよい。例えば、この補完処理の対象となり得る全ての画素を補完することができる。本ステップにより、黒の塗りつぶしとして表現されていたストリップの少なくとも１つの画素の画素値（ゼロ）が、その近傍に位置する１以上の画素（近傍画素）の画素値に基づいて置き換えられる。
（Ｓ４９：画像を表示する）
画像データ処理部２３１Ａは、ステップＳ４８で補完されたストリップを含む全てのストリップに基づいて３次元画像データ（最終的なボリューム、最終的なマージ画像）を構築し、そのレンダリング画像を構築する。主制御部２１１は、このレンダリング画像を表示部２４１に表示させる（エンド）。
本実施形態の眼科検査装置は、第１の実施形態と同様のデータセット収集部、相関係数算出部、及び情報記録部に加えて、情報記録部により記録された走査位置情報に対応する画素の値を１以上の近傍画素の値に基づいて決定する第１画素値決定部を更に含んでいる。本実施形態では、画素値決定部２３５が第１画素値決定部に相当している。
本実施形態の眼科検査装置によれば、第１の実施形態と同様に、従来の技術では記録されていなかった走査位置情報を記録することができ、しかも、相関係数に応じた走査位置情報を記録することができる。更に、本実施形態の眼科検査装置によれば、記録された走査位置情報に基づいて画像補完を行うことが可能である。したがって、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。
＜第３の実施形態＞
第３の実施形態について説明する。本実施形態の眼科検査装置は、第１の実施形態の図１～図４Ａと同様の構成を備えているとともに、図４Ｂの画像データ処理部２３１の代わりに図１１の画像データ処理部２３１Ｂを備えている。本実施形態の画像データ処理部２３１Ｂは、第１の実施形態の画像データ処理部２３１に相関係数マップ作成部２３６と画素値決定部２３７とを付加したものである。
第２の実施形態と同様に、本実施形態の情報記録部２３４は、リサジュースキャンなどの光走査で収集されたデータセットから算出された相関係数集合から、所定の閾値以下の相関係数を要素とする部分集合を特定し、この部分集合に対応する走査位置情報を記録するように構成されている。本実施形態の眼科検査装置は、画像データ処理部２３１Ｂにより構築された画像データ（例えば、画像データ補正部２３３により生成された最終的な３次元画像データ、又は、マージ処理部２３１８により生成されたマージ画像）中の低品質の部分をより高品質な画素に置き換える処理、及び／又は、この画像データ中の画像化されなかった部分を補完する処理を実行する機能要素を備えている。
より詳細には、本実施形態の眼科検査装置は、相関係数の分布を表すマップ（相関係数マップと呼ぶ）を作成し、この相関係数マップに基づき画素値を決定するように構成されている。相関係数マップを作成する機能要素の例が相関係数マップ作成部２３６であり、画素値を決定する機能要素の例が画素値決定部２３７である。
相関係数マップ作成部２３６は、画像データ処理部２３１Ｂ（例えば、相関係数算出部２３１５）により取得された相関係数集合に基づいて相関係数マップを作成するように構成されている。相関係数マップは、リサジュースキャンなどの光走査が適用されたエリアの少なくとも一部における相関係数の分布を表す情報である。相関係数マップ作成部２３６は、例えば、相関係数集合に含まれる各相関係数と、その相関係数に対応する画素位置（走査位置）とを関連付けることによって作成される。相関係数マップ作成部２３６は、プロセッサを含むハードウェアと、相関係数マップ作成処理ソフトウェアとの協働によって実現される。
画素値決定部２３７は、相関係数マップ作成部２３６により作成された相関係数マップに基づいて、情報記録部２３４により記録された走査位置情報に属する走査位置に対応する画素の値を１以上の近傍画素の値に基づいて決定するように構成されている。画素値決定部２３７は、プロセッサを含むハードウェアと、画素値決定処理ソフトウェアとの協働によって実現される。
画素値決定部２３７が実行する画素値決定処理の方法は任意であってよい。例えば、画素値決定部２３７は、第２の実施形態と同様に、任意の補間（内挿）及び／又は任意の外挿を実行可能であってよい。後述の動作例において、画素値決定部２３７が実行する処理の例を説明する。
本実施形態に係る眼科検査装置の動作について説明する。眼科検査装置の動作の例を図１２に示す。特に言及しない限り、第１の実施形態で説明された事項や第２の実施形態で説明された事項をこの動作例に適用することができる。
（Ｓ５１：サイクル個数を指定する）
患者ＩＤの入力や走査モードの設定（リサジュースキャンの指定）の後、任意的な工程として、被検眼Ｅに適用されるリサジュースキャンにおけるサイクルの個数を指定することができる。サイクルの個数は、スキャンの密度に相当する。サイクルの個数の指定は、例えば、ユーザー又は眼科検査装置によって行われる。
本ステップで指定されたサイクルの個数をＧ個とする（Ｇは正の整数）。サイクル個数Ｇは任意であってよく、例えば１０２４個であってよい。第１の実施形態で説明したように、Ｇ個のサイクルは順序付けられており、第ｇ番目のサイクルをＣｇで表す（ｇ＝１、２、・・・、Ｇ）。
（Ｓ５２：準備動作を行う）
眼科検査装置は、固視標の提示、アライメント、フォーカス調整、ＯＣＴ光路長調整など、従来と同様の準備動作を実行する。眼科検査装置は、この準備動作において取得された画像データに基づいて、相関係数集合の部分集合を特定するための判定条件（例えば、前述の第１の閾値）を決定するように構成されていてよい。なお、デフォルトの判定条件を用いてもよいし、ステップＳ５３のリサジュースキャンで収集されたデータセットから決定された判定条件を用いてもよい。
（Ｓ５３：リサジュースキャンを行う）
所定の走査開始トリガー信号を受けて、走査制御部２１１１は、ステップＳ５１で指定されたＧ個のサイクルＣ１～ＣＧからなるリサジュースキャンを被検眼Ｅ（眼底Ｅｆ）に適用する。
（Ｓ５４：各サイクルについてストリップを構築する）
画像データ構築部２２０は、ステップＳ５３のリサジュースキャンで収集されたデータセットに基づいて、ステップＳ５１で指定された各サイクルＣｇに対応するストリップを構築する。この処理は、例えば、第１の実施形態（図６ＡのステップＳ２）と同じ要領で実行されてよい。
ステップＳ５１で指定された各サイクルＣｇについて、このサイクルＣｇに対応するストリップは、このサイクルＣｇに関連付けられた記憶領域Ａｇに記憶される。Ｇ個のサイクルＣ１～ＣＧにそれぞれ関連付けられたＧ個の記憶領域Ａ１～ＡＧは、例えば、記憶部２１２内に設けられる。
（Ｓ５５：各サイクルペアについて相関係数を算出する）
画像データ処理部２３１Ｂ（相関係数算出部２３１５、ｚシフト量算出部２３２など）は、ステップＳ５１で指定されたＧ個のサイクルＣ１～ＣＧに基づき形成された各サイクルペアについての相関係数を算出する。
本ステップの処理は、例えば、Ｇ個のサイクルＣ１～ＣＧのうちから選択された２つのサイクル（サイクルペア）Ｃｇ１及びＣｇ２の間の相関係数を算出する（ここで、ｇ１＝１～Ｇ、ｇ２＝１～Ｇ、ｇ１≠ｇ２）。幾つかの例示的な態様では、画像データ処理部２３１Ｂは、Ｇ個のサイクルＣ１～ＣＧのうちの１つを基準サイクルに設定し、この基準サイクルと他の各サイクルとの間の相関係数を演算する。
その具体例として、画像データ処理部２３１Ｂは、第１番目のサイクルＣ１を基準サイクルに設定し、この基準サイクルＣ１と他のサイクルＣ２～ＣＧのそれぞれとの間の相関係数を算出する。これにより、Ｇ－１個のサイクルペア＜Ｃ１、Ｃ２＞、＜Ｃ１、Ｃ３＞、・・・、＜Ｃ１、ＣＧ＞にそれぞれ対応するＧ－１個の相関係数Ｒ２、Ｒ３、・・・、ＲＧが得られる。ここで、＜Ｃ１、Ｃｇ＞は、第１番目のサイクルと第ｇ番目のサイクルとからなるサイクルペアを示している（ｇ＝２～Ｇ）。
サイクルペア＜Ｃ１、Ｃｇ＞に対応する相関係数Ｒｇは、このサイクルＣｇ（又は、このサイクルペア＜Ｃ１、Ｃｇ＞に関連付けられた記憶領域Ｂｇに記憶される。Ｇ－１個の相関係数Ｒ２～ＲＧにそれぞれ関連付けられたＧ－１個の記憶領域Ｂ１～ＢＧは、例えば、記憶部２１２内に設けられる。
このように、本例で生成される相関係数集合は、Ｇ－１個の相関係数Ｒ２～ＲＧを含んでおり、各相関係数Ｒｇ（ｇ＝２～Ｇ）は、対応するサイクルＣｇ又は対応するサイクルペア＜Ｃ１、Ｃｇ＞に関連付けられている。
（Ｓ５６：画像データを構築する）
画像データ処理部２３１Ｂは、ステップＳ５３のリサジュースキャンで収集されたデータセットに基づいて画像データを構築する。本ステップの処理は、例えば、第１の実施形態における画像データ構築とモーションアーティファクト補正との組み合わせ処理（例えば、図６Ａ及び図６ＢのステップＳ３～ステップＳ１４、図６Ａ～図６ＤのステップＳ３～Ｓ２５、又は、図６Ａ～図６ＥのステップＳ３～ステップＳ３２）と同じ要領で実行されてよい。本ステップで生成された画像データは、例えば記憶部２１２に保存される。
（Ｓ５７：相関係数マップを作成する）
相関係数マップ作成部２３６は、ステップＳ５５で取得された相関係数集合に基づいて、ステップＳ５３のリサジュースキャンの適用エリアにおける相関係数の分布を表す相関係数マップを作成する。例えば、相関係数マップは、相関係数集合に属する各相関係数Ｒｇ（ｇ＝２～Ｇ）と、サイクルＣｇ又はサイクルペア＜Ｃ１、Ｃｇ＞との間の関連付けを表現した情報である。
本ステップの例を説明する。なお、本例では２次元画像データの画素に相関係数を割り当てる場合について説明するが、３次元画像データの画素に相関係数を割り当ててもよい。本実施形態で扱う各種の画像データはいずれも、同じリサジュースキャンで収集されたデータセットに基づくものであり、したがって、或る画像データの画素への相関係数の割り当てと、他の画像データの画素への相関係数の割り当てとを自然に対応付けることが可能である。
本例の相関係数マップは、画像データ（マージ画像）の画素に、対応する相関係数を割り当てることによって作成される。画素に割り当てられる相関係数を画素相関係数と呼ぶ。ここで、１つの画素に対して１つの画素相関係数が割り当てられてもよいし、２以上の画素（画素群）に対して１つの画素相関係数が割り当てられてもよい。後者の例として、基準サイクルＣ１とサイクルＣｇとの交差領域（共通領域）内の各画素に対して同じ画素相関係数を割り当てる場合や、この交差領域を除いたサイクルＣｇの領域（非交差領域）内の各画素に対して同じ画素相関係数を割り当てる場合などがある。
本例の画素相関係数は、ステップＳ５５で求められた相関係数Ｒｇに基づいて決定される。例えば、図１３に示すように、任意の２つのサイクルＣｇ（１）及びＣｇ（２）からなるサイクルペア＜Ｃｇ（１）、Ｃｇ（２）＞（これに対応するストリップペア）について、２つのサイクルＣｇ（１）及びＣｇ（２）（対応するストリップペア）の交差状態をＤ（ｇ（１）、ｇ（２））で表す。ここで、ｇ（１）＝２～Ｇ、ｇ（２）＝２～Ｇ、ｇ（１）≠ｇ（２）とする。
図１３において、２つのサイクルＣｇ（１）及びＣｇ（２）の交差領域は斜線で示されている。本例では、各交差領域に属する各画素に対して、サイクルＣｇ（１）の相関係数とサイクルＣｇ（２）の相関係数とのうち大きい方の値が割り当てられる。更に、サイクルＣｇ（１）の非交差領域（つまり、これらの交差領域を除いたサイクルＣｇ（１）の部分）に属する各画素にはサイクルＣｇ（１）の相関係数が割り当てられるとともに、サイクルＣｇ（２）の非交差領域（つまり、これらの交差領域を除いたサイクルＣｇ（２）の部分）に属する各画素にはサイクルＣｇ（２）の相関係数を割り当てられる。
例えば、サイクルＣｇ（１）の相関係数が０．５であり、サイクルＣｇ（２）の相関係数が０．９である場合、交差領域に属する各画素にはｍａｘ｛０．５、０．９｝＝０．９が割り当てられ、サイクルＣｇ（１）の非交差領域に属する各画素にはそのまま０．５が割り当てられ、サイクルＣｇ（２）の非交差領域に属する各画素にはそのまま０．９が割り当てられる。
本例の基準サイクルＣ１の各画素への相関係数の割り当て方法は、任意であってよい。例えば、サイクルＣ１と異なるサイクルを基準サイクルに設定して同様の処理を行うことによってサイクルＣ１に属する各画素に相関係数を割り当てるようにしてもよい。或いは、サイクルＣ２～ＣＧに属する画素に割り当てられた相関係数から基準サイクルＣ１に属する画素に割り当てられる相関係数の値を決定するようにしてもよい。
このような処理を各サイクルペアについて実行することによって相関係数マップを作成することができる。
幾つかの例示的な態様では、リサジュースキャンに基づく画像データ（マージ画像）の全ての画素に相関係数を割り当てるために、サイクルペアの組み合わせを様々に変更しつつ上記の例のような処理を実行してもよい。例えば次のような一連の処理を実行することができる。
まず、或るサイクルＣｇ（１）を基準サイクルに設定するとともに、この基準サイクルＣｇ（１）と他のサイクルＣｇ（２）～Ｃｇ（Ｇ）との各ペアについて、上記の例の処理を実行する。
次に、リサジュースキャンに基づく画像データ（マージ画像）の全ての画素を網羅するように、全てのサイクルをペアリングする（つまり、複数のペアを設定する）する。更に、設定された各サイクルペアについて、上記の例の処理を実行する。
次に、ステップＳ５６の画像データ構築処理（第１の実施形態の図６Ａ及び図６ＢのステップＳ３～ステップＳ１４で実行される逐次的な処理）と同様に、これらのサイクルに対応するストリップを逐次的に（再帰的に）マージしながら、上記の例の処理を逐次的に（再帰的に）実行する。このようにして構築される画像データの全ての画素には相関係数が割り当てられている。
（Ｓ５８：相関係数が閾値以下のストリップの画素値をゼロに設定する）
画像データ処理部２３１Ｂ（例えば、情報記録部２３４）は、ステップＳ５７で作成された相関係数マップに記録されている各相関係数の値を前述の第１の閾値と比較することで、第１の閾値よりも値が小さい相関係数を特定する。これにより、ステップＳ５７で作成された相関係数マップに記録されている相関係数のうち、第１の閾値よりも値が小さい相関係数が特定される。つまり、相関係数集合の部分集合が特定される。
更に、画像データ処理部２３１Ｂ（例えば、画像データ補正部２３３又は情報記録部２３４）は、この部分集合の各要素（低値の相関係数）に対応する相関係数マップ中の位置に対応する画像データの画素の値をゼロに設定する。
（Ｓ５９：画素値をゼロにした画素の位置（走査位置情報）を記録する）
画像データ処理部２３１Ｂ（情報記録部２３４）は、ステップＳ５８で画素値がゼロに設定された画素の位置を走査位置情報として記録する。
（Ｓ６０：画素値ゼロの画素の補完を行う）
画像データ処理部２３１Ｂ（画素値決定部２３７）は、ステップＳ５９で走査位置情報として記録された画素（すなわち、ステップＳ５８で画素値がゼロに設定された画素）の補完を行う。本ステップにより、黒の塗りつぶしとして表現されていたストリップの少なくとも１つの画素の画素値（ゼロ）が、その近傍に位置する１以上の画素（近傍画素）の画素値に基づいて置き換えられる。
画素値決定部２３７により実行される処理の例を説明する。補完の対象となる画素の位置（座標）を（ｘｉ、ｙｉ）と表し、この対象画素の画素値（輝度値）をＩ（ｘｉ、ｙｉ）と表し、対応する相関係数をＲ（ｘｉ、ｙｉ）と表す。
まず、正規化（規格化）を施す前の対象画素の暫定的な輝度値を算出する。本例では、対象画素の近傍画素として、対象画素の上下左右にそれぞれ位置する４つの画素を考慮する。対象画素に関する表現と同様に、近傍画素の位置（座標）を（ｘｊ、ｙｊ）と表し、この対象画素の画素値（輝度値）をＩ（ｘｊ、ｙｊ）と表し、対応する相関係数をＲ（ｘｊ、ｙｊ）と表す。
本例では、対象画素の暫定輝度値Ｔ（ｘｉ、ｙｉ）は、次式を用いて算出される：Ｔ（ｘｉ、ｙｉ）＝Ｉ（ｘｉ、ｙｉ）×Ｒ（ｘｉ、ｙｉ）＋Ｉ（ｘｉ＋１、ｙｉ＋１）×Ｒ（ｘｉ＋１、ｙｉ＋１）＋Ｉ（ｘｉ＋１、ｙｉ－１）×Ｒ（ｘｉ＋１、ｙｉ－１）＋Ｉ（ｘｉ－１、ｙｉ＋１）×Ｒ（ｘｉ－１、ｙｉ＋１）＋Ｉ（ｘｉ－１、ｙｉ－１）×Ｒ（ｘｉ－１、ｙｉ－１）。この演算式は、対象画素及び４つの近傍画素からなる５つの画素のそれぞれの輝度値Ｉと相関係数Ｒとの積を足し合わせる演算を表している。
なお、対象画素の相関係数Ｒ（ｘｉ、ｙｉ）はゼロであるから、４つの近傍画素のそれぞれの輝度値Ｉと相関係数Ｒとの積を足し合わせる演算式を用いてもよい：Ｔ（ｘｉ、ｙｉ）＝Ｉ（ｘｉ＋１、ｙｉ＋１）×Ｒ（ｘｉ＋１、ｙｉ＋１）＋Ｉ（ｘｉ＋１、ｙｉ－１）×Ｒ（ｘｉ＋１、ｙｉ－１）＋Ｉ（ｘｉ－１、ｙｉ＋１）×Ｒ（ｘｉ－１、ｙｉ＋１）＋Ｉ（ｘｉ－１、ｙｉ－１）×Ｒ（ｘｉ－１、ｙｉ－１）。
続いて、算出された暫定輝度値Ｔ（ｘｉ、ｙｉ）を、対象画素及び４つの近傍画素からなる５つの画素の相関係数の総和で除算することによって、正規化を行う。そして、この正規化により得られた値で、補完前の対象画素の輝度値Ｉ（ｘｉ、ｙｉ）＝０を置き換える。このようにして補完された対象画素の輝度値は次式のように表される：Ｉ（ｘｉ、ｙｉ）＝Ｔ（ｘｉ、ｙｉ）／｛Ｒ（ｘｉ、ｙｉ）＋Ｒ（ｘｉ＋１、ｙｉ＋１）＋Ｒ（ｘｉ＋１、ｙｉ－１）＋Ｒ（ｘｉ－１、ｙｉ＋１）＋Ｒ（ｘｉ－１、ｙｉ－１）｝＝Ｔ（ｘｉ、ｙｉ）／｛Ｒ（ｘｉ＋１、ｙｉ＋１）＋Ｒ（ｘｉ＋１、ｙｉ－１）＋Ｒ（ｘｉ－１、ｙｉ＋１）＋Ｒ（ｘｉ－１、ｙｉ－１）｝。
走査位置情報に記録されている各画素（すなわち、画素値がゼロに設定された各画素）に対してこのような処理を適用することにより、画像情報が無かった領域を補完することができる。
（Ｓ６１：画像を表示する）
画像データ処理部２３１Ｂは、ステップＳ６０で補完された画像データ（例えば、最終的なボリューム、最終的なマージ画像）のレンダリング画像を構築する。主制御部２１１は、このレンダリング画像を表示部２４１に表示させる（エンド）。
本実施形態の眼科検査装置は、第１の実施形態と同様のデータセット収集部、相関係数算出部、及び情報記録部に加えて、相関係数算出部により取得された相関係数集合に基づいて光走査の適用エリアにおける相関係数の分布を表す相関係数マップを作成するマップ作成部を更に含んでいる。本実施形態では、相関係数マップ作成部２３６が相関係数マップに相当している。
本実施形態の眼科検査装置によれば、第１の実施形態と同様に、従来の技術では記録されていなかった走査位置情報を記録することができ、しかも、相関係数に応じた走査位置情報を記録することができる。更に、本実施形態の眼科検査装置によれば、記録された走査位置情報に基づいて相関係数マップを作成することができる。これにより、相関係数の分布状態を把握することや、相関係数の分布状態を解析することや、相関係数の分布状態をユーザーに提供することが可能になり、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。
本実施形態の眼科検査装置は、情報記録部により記録された走査位置情報に対応する画素の値を相関係数マップに基づいて決定する第２画素値決定部を更に含んでいる。本実施形態では、画素値決定部２３７が第２画素値決定部に相当している。
このような本実施形態の眼科検査装置によれば、相関係数マップに基づいて画像補完を行うことが可能になり、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。なお、幾つかの例示的な態様では、第２画素値決定部が設けられていなくてもよい。
＜第４の実施形態＞
第４の実施形態について説明する。本実施形態の眼科検査装置は、第１の実施形態の図１～図４Ａと同様の構成を備えているとともに、図４Ｂの画像データ処理部２３１の代わりに図１４の画像データ処理部２３１Ｃを備えている。本実施形態の画像データ処理部２３１Ｃは、第１の実施形態の画像データ処理部２３１に走査エリア設定部２３８と画像合成処理部２３９とを付加したものである。
第２の実施形態及び第３の実施形態と同様に、本実施形態の情報記録部２３４は、リサジュースキャンなどの光走査で収集されたデータセットから算出された相関係数集合から、所定の閾値以下の相関係数を要素とする部分集合を特定し、この部分集合に対応する走査位置情報を記録するように構成されている。本実施形態の眼科検査装置は、画像データ処理部２３１Ｃにより構築された画像データ（例えば、画像データ補正部２３３により生成された最終的な３次元画像データ、又は、マージ処理部２３１８により生成されたマージ画像）中の低品質の部分や画像化されなかった部分のデータを再度収集する機能要素（つまり、再スキャンを行う機能要素）を備えている。
本実施形態における再スキャンのモードは任意であってよく、例えば、リサジュースキャンであってもよいし、他のスキャンモードであってもよい。再スキャンでリサジュースキャンが実行される場合、この再リサジュースキャンは、例えば、走査位置情報に基づき特定される相関係数が低いサイクルを含んでいてよい。リサジュースキャンと異なるスキャンモードで再スキャンが実行される場合、この再スキャンは、例えば、走査位置情報に基づき特定される相関係数が低い部分（つまり、画像データ処理部２３１Ｃにより構築された画像データ中の低品質の部分や画像化されなかった部分）を含む３次元スキャン（ラスタースキャンなど）であってよい。
本実施形態の眼科検査装置は、情報記録部により記録された走査位置情報に基づいて被検眼Ｅに新たな光走査を適用するように構成されている。走査エリア設定部２３８は、この新たな光走査が適用されるエリア（再走査エリアと呼ぶ）を走査位置情報に基づき設定するように構成されている。再走査エリアの設定は、手動及び／又は自動で実行されてよい。
再走査エリアの設定を手動で行う場合の幾つかの例示的な態様において、主制御部２１１は、走査位置情報に基づく情報（例えば、第３の実施形態の相関係数マップ）、及び／又は、リサジュースキャンに基づき構築された画像（例えば、画像データ補正部２３３により構築された画像）を、再走査エリア設定用グラフィカルユーザーインターフェイス（ＧＵＩ）とともに、表示部２４１に表示させる。ユーザーは、表示された情報及び／又は画像を参照して所望の再走査エリアを決定し、操作部２４２及びＧＵＩを用いてこの再走査エリアを入力する。走査エリア設定部２３８は、例えば、ユーザーにより入力された再走査エリアを走査位置情報及び／又は走査プロトコル２１２１と比較し、必要に応じて再走査エリアを修正し、最終的な再走査エリアを設定するように構成されていてよい。
再走査エリアの設定を自動で行う場合の幾つかの例示的な態様において、走査エリア設定部２３８は、例えば、走査位置情報及び走査プロトコル２１２１に基づいて再走査エリアを設定するように構成されていてよい。
なお、本実施形態の走査位置情報は、品質が低い画像の部分や画像化されなかった部分を示す情報である。走査エリア設定部２３８は、このような部分を含むように再走査エリアの設定を行う。手動設定の場合において、走査エリア設定部２３８は、例えば、走査位置情報に示されている品質が低い画像の部分や画像化されなかった部分が、ユーザーが指定した再走査エリアに含まれているか判定し、含まれていないと判定された場合には、それらの部分を含むように再走査エリアの修正が実行される。自動設定の場合においては、走査エリア設定部２３８は、例えば、走査位置情報に示されている品質が低い画像の部分や画像化されなかった部分を含むように再走査エリアが設定される。
画像合成処理部２３９は、画像データ処理部２３１Ｃの他の要素と協働することにより、最初に実行された光走査（リサジュースキャン）で収集されたデータセット（第１データセットと呼ぶ）と、再スキャンで収集されたデータセット（第２データセットと呼ぶ）とを処理して画像データを構築するように構成されている。
例えば、画像合成処理部２３９は、第１データセット又はこれに基づく画像データ（第１画像データと呼ぶ）と、第２データセット又はこれに基づく画像データ（第２画像データと呼ぶ）とを合成（マージ）するように構成されている。第１画像データは、走査位置情報に示されている品質が低い領域や画像化されなかった領域を含んでいる。一方、第２画像データは、このような領域を含むように実行された再スキャンで収集された第２データセットから構築されたデータである。このような第１画像データ及び第２画像データを合成することで、第１画像データを第２画像データで補完することが可能である。
幾つかの例示的な態様において、画像データ処理部２３１Ｃは、次の３つの処理工程を実行するように構成されていてよい。まず、画像データ処理部２３１Ｃは、第１の実施形態と同様の処理を第１データセットに適用することで、走査位置情報に対応するエリアの少なくとも一部を最初のリサジュースキャンの適用エリアから除いたエリアに対応する第１画像データを構築する（第１の処理工程）。また、画像データ処理部２３１Ｃは、第１の実施形態と同様の処理を第２データセットに適用することで、再走査エリアに対応する第２画像データを構築する（第２の処理工程）。更に、画像データ処理部２３１Ｃ（画像合成処理部２３９）は、第１の処理工程で構築された第１画像データと、第２の処理工程で構築された第２画像データとを合成する（第３の処理工程）。これらの処理工程を実行することで、第１画像データと第２画像データとの合成を好適に行うことができる。例えば、第１画像データと第２画像データとの間に重複エリアが存在するように第１画像データの構築エリアを設定するとともに、第１画像データにおいてこの重複エリアに相当する領域と、第２画像データにおいてこの重複エリアに相当する領域との間の位置合わせを行うことによって、第１画像データと第２画像データとの間のレジストレーションの正確度及び精度を向上させることが可能になる。
本実施形態に係る眼科検査装置の動作について説明する。眼科検査装置の動作の例を図１５に示す。特に言及しない限り、第１～第３の実施形態で説明された事項をこの動作例に適用することができる。
（Ｓ７１：サイクル個数を指定する）
患者ＩＤの入力や走査モードの設定（リサジュースキャンの指定）の後、任意的な工程として、被検眼Ｅに適用されるリサジュースキャンにおけるサイクルの個数を指定することができる。サイクルの個数は、スキャンの密度に相当する。サイクルの個数の指定は、例えば、ユーザー又は眼科検査装置によって行われる。
本ステップで指定されたサイクルの個数をＧ個とする（Ｇは正の整数）。サイクル個数Ｇは任意であってよく、例えば１０２４個であってよい。第１の実施形態で説明したように、Ｇ個のサイクルは順序付けられており、第ｇ番目のサイクルをＣｇで表す（ｇ＝１、２、・・・、Ｇ）。
（Ｓ７２：準備動作を行う）
眼科検査装置は、固視標の提示、アライメント、フォーカス調整、ＯＣＴ光路長調整など、従来と同様の準備動作を実行する。眼科検査装置は、この準備動作において取得された画像データに基づいて、相関係数集合の部分集合を特定するための判定条件（例えば、前述の第１の閾値）を決定するように構成されていてよい。なお、デフォルトの判定条件を用いてもよいし、ステップＳ７３のリサジュースキャンで収集されたデータセットから決定された判定条件を用いてもよい。
（Ｓ７３：リサジュースキャンを行う）
所定の走査開始トリガー信号を受けて、走査制御部２１１１は、ステップＳ７１で指定されたＧ個のサイクルＣ１～ＣＧからなるリサジュースキャンを被検眼Ｅ（眼底Ｅｆ）に適用する。
（Ｓ７４：各サイクルについてストリップを構築する）
画像データ構築部２２０は、ステップＳ７３のリサジュースキャンで収集されたデータセットに基づいて、ステップＳ７１で指定された各サイクルＣｇに対応するストリップを構築する。この処理は、例えば、第１の実施形態（図６ＡのステップＳ２）と同じ要領で実行されてよい。
ステップＳ７１で指定された各サイクルＣｇについて、このサイクルＣｇに対応するストリップは、このサイクルＣｇに関連付けられた記憶領域Ａｇに記憶される。Ｇ個のサイクルＣ１～ＣＧにそれぞれ関連付けられたＧ個の記憶領域Ａ１～ＡＧは、例えば、記憶部２１２内に設けられる。
（Ｓ７５：各サイクルペアについて相関係数を算出する）
画像データ処理部２３１Ｃ（相関係数算出部２３１５、ｚシフト量算出部２３２など）は、ステップＳ７１で指定されたＧ個のサイクルＣ１～ＣＧに基づき形成された各サイクルペアについての相関係数を算出する。
本ステップの処理は、例えば、Ｇ個のサイクルＣ１～ＣＧのうちから選択された２つのサイクル（サイクルペア）Ｃｇ１及びＣｇ２の間の相関係数を算出する（ここで、ｇ１＝１～Ｇ、ｇ２＝１～Ｇ、ｇ１≠ｇ２）。幾つかの例示的な態様では、画像データ処理部２３１Ｃは、Ｇ個のサイクルＣ１～ＣＧのうちの１つを基準サイクルに設定し、この基準サイクルと他の各サイクルとの間の相関係数を演算する。
その具体例として、画像データ処理部２３１Ｃは、第１番目のサイクルＣ１を基準サイクルに設定し、この基準サイクルＣ１と他のサイクルＣ２～ＣＧのそれぞれとの間の相関係数を算出する。これにより、Ｇ－１個のサイクルペア＜Ｃ１、Ｃ２＞、＜Ｃ１、Ｃ３＞、・・・、＜Ｃ１、ＣＧ＞にそれぞれ対応するＧ－１個の相関係数Ｒ２、Ｒ３、・・・、ＲＧが得られる。ここで、＜Ｃ１、Ｃｇ＞は、第１番目のサイクルと第ｇ番目のサイクルとからなるサイクルペアを示している（ｇ＝２～Ｇ）。
サイクルペア＜Ｃ１、Ｃｇ＞に対応する相関係数Ｒｇは、このサイクルＣｇ（又は、このサイクルペア＜Ｃ１、Ｃｇ＞に関連付けられた記憶領域Ｂｇに記憶される。Ｇ－１個の相関係数Ｒ２～ＲＧにそれぞれ関連付けられたＧ－１個の記憶領域Ｂ１～ＢＧは、例えば、記憶部２１２内に設けられる。
このように、本例で生成される相関係数集合は、Ｇ－１個の相関係数Ｒ２～ＲＧを含んでおり、各相関係数Ｒｇ（ｇ＝２～Ｇ）は、対応するサイクルＣｇ又は対応するサイクルペア＜Ｃ１、Ｃｇ＞に関連付けられている。
（Ｓ７６：相関係数が閾値以下のサイクルの識別子を記録する）
画像データ処理部２３１Ｃ（情報記録部２３４）は、ステップＳ７５で算出された各サイクルペアの相関係数の値を前述の第１の閾値と比較することで、第１の閾値よりも値が小さい相関係数を特定する。これにより、ステップＳ７５で取得された相関係数集合の要素（相関係数）のうち、第１の閾値よりも値が小さい相関係数が特定される。つまり、相関係数集合の部分集合が特定される。
画像データ処理部２３１Ｃ（情報記録部２３４）は、この部分集合に対応する走査位置情報を記録する。この走査位置情報は、予め決められた識別子によって表現されている。この識別子は、例えば、サイクル番号、ストリップ番号、及び、画像データの定義座標系の座標のいずれかであってよい。
（Ｓ７７：再スキャンを行う）
走査エリア設定部２３８は、ステップＳ７６で生成された走査位置情報（サイクルの識別子）に基づいて再走査エリアを設定する。走査制御部２１１１は、この再走査エリアに基づいて被検眼Ｅに再スキャンを適用する。
（Ｓ７８：再スキャンの各サイクルについてストリップを構築する）
画像データ構築部２２０は、ステップＳ７７の再スキャンで収集されたデータセットに基づいて、この再スキャンの各サイクルに対応するストリップを構築する。この処理は、例えば、第１の実施形態（図６ＡのステップＳ２）と同じ要領で実行されてよい。
（Ｓ７９：再スキャンに基づく各ストリップでリサジュースキャンに基づく対応ストリップを置換する）
眼科検査装置（例えば、主制御部２１１又は画像データ処理部２３１Ｃ）は、ステップＳ７８で構築された各ストリップについて、ステップＳ７４で構築されたストリップのうちから対応ストリップを特定する。ここで、リサジュースキャンに基づくストリップと再スキャンに基づくストリップとの間の対応は、サイクルの識別子を参照して行うことが可能である。
眼科検査装置（例えば、主制御部２１１又は画像データ処理部２３１Ｃ）は、ステップＳ７８で構築された各ストリップについて、特定された対応ストリップ（リサジュースキャンに基づくストリップ）を当該ストリップ（再スキャンに基づくストリップ）に置き換える。
これにより、リサジュースキャンで得られた画像のうち低品質の部分（相関係数が低値のサイクルに対応するストリップ）を、再スキャンで得られた画像に置換することができる。なお、再スキャンで得られたストリップに基づいて同様の相関係数演算を実行し、品質の再評価を行うようにしてもよい。再スキャンで得られたストリップの品質が低い場合には、更なるスキャンを行うことができる。
（Ｓ８０：３次元画像データを構築する）
画像データ処理部２３１Ｃは、ステップＳ７９によって得られたストリップ集合（すなわち、リサジュースキャンに基づく高品質のストリップと、再スキャンに基づくストリップとからなる集合）に基付いて、３次元画像データを構築する。本ステップの処理は、例えば、第１の実施形態における画像データ構築とモーションアーティファクト補正との組み合わせ処理（例えば、図６Ａ～図６ＥのステップＳ３～ステップＳ３２）と同じ要領で実行されてよい。本ステップで生成された画像データは、例えば記憶部２１２に保存される。
（Ｓ８１：画像を表示する）
画像データ処理部２３１Ｃは、ステップＳ８０で補完された３次元画像データのレンダリング画像を構築する。主制御部２１１は、このレンダリング画像を表示部２４１に表示させる（エンド）。
本実施形態の眼科検査装置は、第１の実施形態と同様のデータセット収集部、相関係数算出部、及び情報記録部を含んでいる。本実施形態のデータセット収集部は、互いに交差する一連のサイクルを含む２次元パターンに従う第１光走査（リサジュースキャンなど）を被検眼に適用するとともに、情報記録部により記録された走査位置情報に基づいて被検眼に第２光走査（再スキャン）を適用するように構成されている。
本実施形態の眼科検査装置によれば、第１の実施形態と同様に、従来の技術では記録されていなかった走査位置情報を記録することができ、しかも、相関係数に応じた走査位置情報を記録することができる。更に、本実施形態の眼科検査装置によれば、記録された走査位置情報に基づいて第２光走査（再スキャン）を行うことができる。これにより、第１光走査で得られた情報に基づき第２光走査を行うことが可能になり、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。
更に、本実施形態の眼科検査装置のデータセット収集部は、情報記録部により記録された走査位置情報に基づき設定されたエリアに第２光走査を適用するように構成されていてよい。
このような本実施形態の眼科検査装置によれば、第１光走査で得られた情報から再走査エリアの設定を自動で行い、この再走査エリアに基づき第２光走査を行うことが可能になる。本実施形態では、走査エリア設定部２３８によって再走査エリアの自動設定が実行される。
なお、幾つかの例示的な態様では、再走査エリアの設定を手動で行うようにしてもよい。例えば、幾つかの例示的な態様の眼科検査装置は、第３の実施形態の相関係数マップ又はこれに基づく情報を表示し、この表示に基づきユーザーが指定した再走査エリアを用いた再スキャンを実行するように構成されていてよい。
更に、本実施形態の眼科検査装置は、再スキャンで収集されたデータセットを画像データ構築に利用するように構成されていてよい。そのために、本実施形態の眼科検査装置は、第１光走査（リサジュースキャンなど）で収集された第１データセットと第２光走査（再スキャン）で収集された第２データセットとを処理して画像データを構築する処理部を含んでいてよい。この処理部は、画像合成処理部２３９を含んでおり、画像データ処理部２３１Ｃの他の要素を更に含んでいてもよい。
このような本実施形態の眼科検査装置によれば、第１光走査で収集されたデータセットと第２光走査で収集されたデータセットとを組み合わせて画像データを構築することが可能になるので、より高品質の画像データを取得することが可能になる。
更に、本実施形態の眼科検査装置の処理部は、次の一連の工程を実行可能であってよい。第１工程において、処理部は、第１光走査で収集された第１データセットを処理することで、情報記録部により記録された走査位置情報に対応するエリアの少なくとも一部を第１光走査の適用エリアから除いたエリア（対象エリアと呼ぶ）の第１画像データを構築する。第２工程において、処理部は、第２光走査で収集された第２データセットを処理することで、走査位置情報に基づき設定されたエリア（つまり、走査エリア設定部２３８により設定された再走査エリア）の第２画像データを構築する。第３工程において、処理部は、第１工程で構築された第１画像データと、第２工程で構築された第２画像データとを合成する。
このような本実施形態の眼科検査装置によれば、第１光走査で収集されたデータセットから構築された第１画像データと第２光走査で収集されたデータセットから構築された第２画像データとの合成画像データを生成することが可能になるので、より高品質の画像データを取得することが可能になる。
＜第５の実施形態＞
第５の実施形態について説明する。第１～第４の実施形態では、光走査を実行可能に構成された走査型イメージング装置について説明したが、本実施形態では、光走査により収集されたデータセットを外部から受け付けてこれを処理するように構成された情報処理装置について説明する。
本実施形態の情報処理装置は、走査型イメージング装置の一部として構成されていてもよいが、そのような場合は除外して以下の説明を行う。ただし、以下の説明における任意の事項が、走査型イメージング装置の一部として構成された情報処理装置にも適用可能であることは、当業者であれば理解できるであろう。
本実施形態の情報処理装置は、光走査を実行する機能要素（第１～第４の実施形態におけるデータセット収集部）を有していないことを除いて、第１～第４の実施形態に係る走査型イメージング装置（眼科検査装置）と同様に構成されていてよい。
第１～第４の実施形態において説明された任意の事項を本実施形態の情報処理装置に組み合わせることが可能である。情報処理装置は、組み合わせられた事項に応じた作用及び効果を奏する。
本実施形態の情報処理装置の構成例を図１６に示す。この情報処理装置３００は、受付部３０１と、相関係数算出部３０２と、情報記録部３０３とを含んでいる。
受付部３０１は、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプル（例えば被検眼）に適用して収集されたデータセットを受け付ける。すなわち、受付部３０１は、リサジュースキャンなどの光走査を実行可能な走査型イメージング装置のデータ収集部によってサンプルから収集されたデータセットを受け付ける。
受付部３０１は、外部装置からデータを取得する機能を有する。この外部装置は、例えば、コンピュータ、記憶装置、記録媒体、情報システム（例えば、病院情報システム、電子カルテシステム、画像アーカイビングシステム）などであってよい。また、外部装置は、例えば、情報処理装置３００に直接に接続された装置、ローカルエリアネットワーク（ＬＡＮ）を介して情報処理装置３００に接続された装置、ワイドエリアネットワーク（ＷＡＮ）を介して情報処理装置３００に接続された装置などであってよい。受付部３０１は、例えば、通信回線を介してデータ通信を行うための通信インターフェイス、及び／又は、記録媒体からデータを読み取るためのドライブ装置を含んでいてもよい。
相関係数算出部３０２は、受付部３０１により外部から受け付けられたデータセットに基づいてサイクル間の相関係数を算出するように構成されている。相関係数算出部３０２は、第１の実施形態の画像データ処理部２３１と同様に構成されていてよい。
幾つかの例示的な態様において、相関係数算出部３０２は、第１の実施形態の画像データ構築部２２０、マスク画像生成部２３１１、範囲調整部２３１２、合成画像生成部２３１３、相互相関関数算出部２３１４、及び相関係数算出部２３１５を含んでいてよく、ｚシフト量算出部２３２を更に含んでいてもよい。なお、受付部３０１により受け付けられるデータセットが画像データ（例えばストリップ）である場合には、相関係数算出部３０２は、画像データ構築部２２０に相当する機能要素を有していなくてもよい。
情報記録部３０３は、受付部３０１により受け付けられたデータセットから相関係数算出部３０２が算出した相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録する。情報記録部３０３は、第１の実施形態の情報記録部２３４と同様に構成されていてよい。
図示は省略するが、本実施形態の情報処理装置３００は、第２の実施形態の画素値決定部２３５、第３の実施形態の相関係数マップ作成部２３６（及び画素値決定部２３７）、及び、第４の実施形態の走査エリア設定部２３８（及び画像合成処理部２３９）のうちの少なくとも１つを備えていてもよい。なお、本実施形態の情報処理装置３００に組み合わせることが可能な機能要素はこれらに限定されず、第１～第４の実施形態において説明された任意の事項であってよい。
このような情報処理装置３００によれば、従来の技術では記録されていなかった走査位置情報を記録することができ、しかも、相関係数に応じた走査位置情報を記録することができるので、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。この利用法については、その幾つかの非限定的な例が第１～第４の実施形態において説明されている。
＜他の実施形態＞
本開示に係る実施形態は第１～第５の実施形態及びそれらの変形例に限定されない。例えば、第１～第５の実施形態は、各種の方法に係る実施形態、各種のプログラムに係る実施形態、各種の記録媒体に係る実施形態などを提供するものである。
例えば、第１の実施形態は、走査型イメージング装置を制御する方法の実施形態を提供する。この実施形態は、光走査を行うスキャナと、この光走査で収集されたデータから画像データを構築するプロセッサとを含む走査型イメージング装置を制御する方法である。このスキャナは第１の実施形態のデータセット収集部であってよく、このプロセッサは第１の実施形態の画像データ構築部２２０及び画像データ処理部２３１（画像データ処理部２３１の少なくとも一部）であってよい。
本実施形態に係る走査型イメージング装置の制御方法は、次の３つの工程を含んでいてよい。第１工程は、スキャナを制御することにより、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータセットを収集させるように実行される。第２工程は、プロセッサを制御することにより、第１工程で収集されたデータセットに基づいてサイクル間の相関係数を算出させるように実行される。第３工程は、プロセッサを制御することにより、第１工程で収集されたデータセットから第２工程で算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録させるように実行される。
第１～第４の実施形態において説明された任意の事項を、本実施形態に係る走査型イメージング装置の制御方法に組み合わせることが可能である。走査型イメージング装置の制御方法は、組み合わせられた事項に応じた作用及び効果を奏する。
このような走査型イメージング装置の制御方法によれば、従来の技術では記録されていなかった走査位置情報を記録することができ、しかも、相関係数に応じた走査位置情報を記録することができるので、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。この利用法については、その幾つかの非限定的な例が第１～第４の実施形態において説明されている。
第５の実施形態は、情報処理装置を制御する方法の実施形態を提供する。この実施形態は、プロセッサを含む情報処理装置を制御する方法である。このプロセッサは第１の実施形態の画像データ構築部２２０及び画像データ処理部２３１（画像データ処理部２３１の少なくとも一部）であってよい。
本実施形態に係る情報処理装置の制御方法は、次の３つの工程を含んでいてよい。第１工程は、プロセッサを制御することにより、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用して収集されたデータセットを受け付けさせるように実行される。第２工程は、プロセッサを制御することにより、第１工程で受け付けられたデータセットに基づいてサイクル間の相関係数を算出させるように実行される。第３工程は、プロセッサを制御することにより、第１工程で受け付けられたデータセットから第２工程で算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録させるように実行される。
第１～第４の実施形態において説明された任意の事項を、本実施形態に係る情報処理装置の制御方法に組み合わせることが可能である。情報処理装置の制御方法は、組み合わせられた事項に応じた作用及び効果を奏する。
このような情報処理装置の制御方法によれば、従来の技術では記録されていなかった走査位置情報を記録することができ、しかも、相関係数に応じた走査位置情報を記録することができるので、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。この利用法については、その幾つかの非限定的な例が第１～第４の実施形態において説明されている。
第１の実施形態及び／又は第５の実施形態は、情報を処理する方法の実施形態を提供する。
本実施形態に係る情報処理方法は、次の３つの工程を含んでいてよい。第１工程は、互いに交差する一連のサイクルを含む２次元パターンに従う光走査でサンプルから収集されたデータセットを取得するように実行される。第１の実施形態に基づく情報処理方法では、第１工程は、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータセットを収集するように実行される。第５の実施形態に基づく情報処理方法では、第１工程は、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用して収集されたデータセットを外部から受け付けるように実行される。第２工程は、第１工程で取得されたデータセットに基づいてサイクル間の相関係数を算出するように実行される。第３工程は、第１工程で取得されたデータセットから第２工程で算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録させるように実行される。
第１～第４の実施形態において説明された任意の事項を、本実施形態に係る情報処理方法に組み合わせることが可能である。情報処理方法は、組み合わせられた事項に応じた作用及び効果を奏する。
このような情報処理方法によれば、従来の技術では記録されていなかった走査位置情報を記録することができ、しかも、相関係数に応じた走査位置情報を記録することができるので、走査型イメージングで取得されたデータの新規な利用法を提供することが可能になる。この利用法については、その幾つかの非限定的な例が第１～第４の実施形態において説明されている。
本開示は、実施形態に係る方法をコンピュータに実行させるプログラムを提供する。幾つかの例示的な実施形態は、走査型イメージング装置の制御方法をコンピュータに実行させるプログラムを提供する。幾つかの例示的な実施形態は、情報処理装置の制御方法をコンピュータに実行させるプログラムを提供する。幾つかの例示的な実施形態は、情報処理方法をコンピュータに実行させるプログラムを提供する。
本開示は、実施形態に係るプログラムが記録されたコンピュータ可読な非一時的記録媒体を提供する。幾つかの例示的な実施形態は、走査型イメージング装置の制御方法をコンピュータに実行させるプログラムが記録されたコンピュータ可読な非一時的記録媒体を提供する。幾つかの例示的な実施形態は、情報処理装置の制御方法をコンピュータに実行させるプログラムが記録されたコンピュータ可読な非一時的記録媒体を提供する。幾つかの例示的な実施形態は、情報処理方法をコンピュータに実行させるプログラムが記録されたコンピュータ可読な非一時的記録媒体を提供する。実施形態の非一時的記録媒体は任意の形態であってよく、その例として、磁気ディスク、光ディスク、光磁気ディスク、半導体メモリなどがある。
本開示は、幾つかの例示的な態様を例示するものに過ぎず、発明の限定を意図したものではない。この発明を実施しようとする者は、この発明の要旨の範囲内における任意の変形（省略、置換、付加など）を施すことが可能である。
１ 眼科検査装置
２３１、２３１Ａ、２３１Ｂ、２３１Ｃ 画像データ処理部
２３１５ 相関係数算出部
２３２ ｚシフト量算出部
２３４ 情報記録部
２３５ 画素値決定部
２３６ 相関係数マップ作成部
２３７ 画素値決定部
２３８ 走査エリア設定部
２３９ 画像合成処理部

Claims (16)

  1. 光走査を用いて画像データを構築する走査型イメージング装置であって、
    互いに交差する一連のサイクルを含む２次元パターンに従う第１光走査をサンプルに適用して第１データセットを収集するデータセット収集部と、
    前記第１データセットに基づいてサイクル間の相関係数を算出する相関係数算出部と、
    前記相関係数算出部により前記第１データセットから算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録する情報記録部と
    を含む、走査型イメージング装置。
  2. 前記部分集合の要素は、所定の閾値以下の相関係数である、
    請求項１の走査型イメージング装置。
  3. 前記情報記録部により記録された前記走査位置情報に対応する画素の値を１以上の近傍画素の値に基づいて決定する第１画素値決定部を更に含む、
    請求項２の走査型イメージング装置。
  4. 前記相関係数算出部により取得された前記相関係数の前記集合に基づいて、前記第１光走査の適用エリアにおける相関係数マップを作成するマップ作成部を更に含む、
    請求項２の走査型イメージング装置。
  5. 前記情報記録部により記録された前記走査位置情報に対応する画素の値を前記相関係数マップに基づいて決定する第２画素値決定部を更に含む、
    請求項４の走査型イメージング装置。
  6. 前記データセット収集部は、前記情報記録部により記録された前記走査位置情報に基づいて前記サンプルに第２光走査を適用する、
    請求項２の走査型イメージング装置。
  7. 前記データセット収集部は、前記情報記録部により記録された前記走査位置情報に基づき設定されたエリアに前記第２光走査を適用する、
    請求項６の走査型イメージング装置。
  8. 前記第１光走査で収集された前記第１データセットと前記第２光走査で収集された第２データセットとを処理して画像データを構築する処理部を更に含む、
    請求項７の走査型イメージング装置。
  9. 前記処理部は、
    前記第１データセットを処理して、前記情報記録部により記録された前記走査位置情報に対応するエリアの少なくとも一部を前記第１光走査の適用エリアから除いたエリアの第１画像データを構築し、
    前記第２データセットを処理して、前記走査位置情報に基づき設定された前記エリアの第２画像データを構築し、
    前記第１画像データと前記第２画像データとを合成する、
    請求項８の走査型イメージング装置。
  10. 前記データセット収集部は、
    光走査のための光を互いに異なる２つの方向に偏向可能な偏向器を含み、
    前記２つの方向のうちの一方の方向に沿った偏向方向の変化を第１周期で繰り返しつつ他方の方向に沿った偏向方向の変化を前記第１周期と異なる第２周期で繰り返すことによって前記第１光走査を前記サンプルに適用する、
    請求項１～９のいずれかの走査型イメージング装置。
  11. 互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用して収集されたデータセットを受け付ける受付部と、
    前記データセットに基づいてサイクル間の相関係数を算出する相関係数算出部と、
    前記相関係数算出部により前記データセットから算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録する情報記録部と
    を含む、情報処理装置。
  12. 光走査を行うスキャナと前記光走査で収集されたデータから画像データを構築するプロセッサとを含む走査型イメージング装置を制御する方法であって、
    前記スキャナを制御して、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータセットを収集させ、
    前記プロセッサを制御して、前記データセットに基づいてサイクル間の相関係数を算出させ、
    前記プロセッサを制御して、前記データセットから算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録させる、
    方法。
  13. プロセッサを含む情報処理装置を制御する方法であって、
    前記プロセッサを制御して、互いに交差する一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用して収集されたデータセットを受け付けさせ、
    前記プロセッサを制御して、前記データセットに基づいてサイクル間の相関係数を算出させ、
    前記プロセッサを制御して、前記データセットから算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録させる、
    方法。
  14. 互いに交差する一連のサイクルを含む２次元パターンに従う光走査でサンプルから収集されたデータセットを取得し、
    前記データセットに基づいてサイクル間の相関係数を算出し、
    前記データセットから算出された相関係数の集合のうち所定の条件を満足する相関係数を要素とする部分集合に対応する走査位置情報を記録する、
    情報処理方法。
  15. 請求項１２～１４のいずれかの方法をコンピュータに実行させるプログラム。
  16. 請求項１５のプログラムが記録されたコンピュータ可読な非一時的記録媒体。

JP2021157561A 2021-09-28 2021-09-28 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体 Pending JP2023048318A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021157561A JP2023048318A (ja) 2021-09-28 2021-09-28 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021157561A JP2023048318A (ja) 2021-09-28 2021-09-28 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体

Publications (1)

Publication Number Publication Date
JP2023048318A true JP2023048318A (ja) 2023-04-07

Family

ID=85779940

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021157561A Pending JP2023048318A (ja) 2021-09-28 2021-09-28 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体

Country Status (1)

Country Link
JP (1) JP2023048318A (ja)

Similar Documents

Publication Publication Date Title
US12039704B2 (en) Image processing apparatus, image processing method and computer-readable medium
CN112822972A (zh) 图像处理装置、图像处理方法和程序
JP7368568B2 (ja) 眼科装置
JP7374615B2 (ja) 情報処理装置、情報処理方法及びプログラム
JP7384987B2 (ja) 眼科装置
JP2021083919A (ja) 医用画像処理装置、光干渉断層撮影装置、医用画像処理方法、及びプログラム
JP7220509B2 (ja) 眼科装置及び眼科画像処理方法
JP2023014190A (ja) 眼科撮影装置
WO2020050308A1 (ja) 画像処理装置、画像処理方法及びプログラム
WO2021153087A1 (ja) 眼科装置、その制御方法、及び記録媒体
JP2021092445A (ja) 走査型イメージング装置、その制御方法、画像処理装置、その制御方法、走査型イメージング方法、画像処理方法、プログラム、及び記録媒体
US20220273170A1 (en) Method of processing ophthalmic data, ophthalmic data processing apparatus, and ophthalmic examination apparatus
JP7410481B2 (ja) 画像処理方法、走査型イメージング方法、画像処理装置、その制御方法、走査型イメージング装置、その制御方法、プログラム、及び記録媒体
JP2021013491A (ja) 光コヒーレンストモグラフィ装置、その制御方法、光計測方法、プログラム、及び記憶媒体
JP2022173339A (ja) 眼科撮影装置、その制御方法、プログラム、及び記録媒体
WO2021153086A1 (ja) 眼科装置、その制御方法、及び記録媒体
JP2023048318A (ja) 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体
JP7286283B2 (ja) 眼科装置
WO2023162586A1 (ja) 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体
JP2023056183A (ja) 走査型イメージング装置、情報処理装置、走査型イメージング装置の制御方法、情報処理装置の制御方法、情報処理方法、プログラム、及び記録媒体
JP7355331B2 (ja) 走査型イメージング装置、その制御方法、画像処理装置、その制御方法、走査型イメージング方法、画像処理方法、プログラム、及び記録媒体
JP2023070264A (ja) 走査型イメージング装置、その制御方法、プログラム、及び記録媒体
JP7341422B2 (ja) 走査型イメージング装置、その制御方法、走査型イメージング方法、プログラム、及び記録媒体
JP7325675B2 (ja) 眼科撮影装置
WO2022085501A1 (ja) 眼科装置、その制御方法、プログラム、及び記録媒体

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240919