特許 7410481 画像処理方法、走査型イメージング方法、画像処理装置、その制御方法、走査型イメージング装置、その制御方法、プログラム、及び記録媒体

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2023-12-26

(45)【発行日】2024-01-10

(54)【発明の名称】画像処理方法、走査型イメージング方法、画像処理装置、その制御方法、走査型イメージング装置、その制御方法、プログラム、及び記録媒体

(51)【国際特許分類】

   A61B 3/10 20060101AFI20231227BHJP

   G01N 21/17 20060101ALI20231227BHJP

【ＦＩ】

A61B3/10 100

G01N21/17 630

【請求項の数】 24

(21)【出願番号】P 2019223138

(22)【出願日】2019-12-10

(65)【公開番号】P2021090631

(43)【公開日】2021-06-17

【審査請求日】2022-11-24

(73)【特許権者】

【識別番号】504171134

【氏名又は名称】国立大学法人 筑波大学

(73)【特許権者】

【識別番号】000220343

【氏名又は名称】株式会社トプコン

(74)【代理人】

【識別番号】100124626

【弁理士】

【氏名又は名称】榎並 智和

(72)【発明者】

【氏名】巻田 修一

(72)【発明者】

【氏名】安野 嘉晃

(72)【発明者】

【氏名】山口 達夫

(72)【発明者】

【氏名】澁谷 雅博

【審査官】冨永 昌彦

(56)【参考文献】

【文献】特開２０１８－１４０００４（ＪＰ，Ａ）

【文献】特開２０１８－１４００４９（ＪＰ，Ａ）

【文献】米国特許出願公開第２０１８／０２４２８３９（ＵＳ，Ａ１）

【文献】特表２０１６－５１４９７４（ＪＰ，Ａ）

【文献】米国特許出願公開第２０１５／０３２４９６６（ＵＳ，Ａ１）

【文献】YIWEI Chen, et al.，Three-dimensional eye motion correction by Lissajous scan optical coherence tomography，Biomedical Optics EXPRESS，米国，Optical Society of America，2017年02月23日，Vol.8 No.3，1783-1802

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ３／００ － ３／１８

Ｇ０１Ｎ ２１／００ － ２１／６１

(57)【特許請求の範囲】

【請求項1】

通信回線若しくは記録媒体を介してサンプルの第１画像及び第２画像を受け付け、又は、サンプルに光走査を適用して収集されたデータに基づいて第１画像及び第２画像を構築し、
通信回線若しくは記録媒体を介してマスク画像を受け付け、又は、サンプルの画像に基づいてマスク画像を生成し、
前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲と前記マスク画像の画素値の範囲との間の差を小さくするように、前記第１画像の画素値の範囲と前記第２画像の画素値の範囲と前記マスク画像の画素値の範囲とを相対的に調整し、
前記第１画像及び前記第２画像のそれぞれに前記マスク画像を合成して第１合成画像及び第２合成画像を生成し、
前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求め、
前記複数の相互相関関数に基づいて相関係数を算出し、
前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行う、
画像処理方法。

【請求項2】

前記サンプルの画像及び前記マスク画像の一方の画素値の範囲を他方の画素値の範囲に一致させる、
請求項１の画像処理方法。

【請求項3】

前記マスク画像の画素値の範囲に応じて前記サンプルの画像の画素値の範囲を正規化する、
請求項２の画像処理方法。

【請求項4】

前記マスク画像の画素値の範囲は閉区間［０，１］に含まれ、
前記サンプルの画像の各画素の値を、当該画像における最大画素値で除算する、
請求項３の画像処理方法。

【請求項5】

前記マスク画像の画素値の範囲は閉区間［０，１］に含まれ、
前記サンプルの画像の各画素の値を、当該画像の画素値の範囲の最大値で除算する、
請求項３の画像処理方法。

【請求項6】

前記マスク画像は、画素値が０又は１の二値画像である、
請求項１～５のいずれかの画像処理方法。

【請求項7】

前記マスク画像は、前記サンプルの画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像である、
請求項６の画像処理方法。

【請求項8】

前記マスク画像は、第１マスク画像と第２マスク画像とを含み、
前記第１マスク画像を前記第１画像に合成して前記第１合成画像を生成し、
前記第２マスク画像を前記第２画像に合成して前記第２合成画像を生成する、
請求項１～７のいずれかの画像処理方法。

【請求項9】

前記第１マスク画像及び前記第２マスク画像は、同じ寸法及び同じ形状を有する、
請求項８の画像処理方法。

【請求項10】

前記第１マスク画像は、前記第１画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であり、
前記第２マスク画像は、前記第２画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像である、
請求項８又は９の画像処理方法。

【請求項11】

前記第１画像及び前記第２画像のそれぞれは、一連のサイクルを含む２次元パターンに従う走査を前記サンプルに適用して収集されたデータに基づき構築される、
請求項１～１０のいずれかの画像処理方法。

【請求項12】

前記第１画像及び前記第２画像のそれぞれは、リサージュ関数に基づき予め設定された走査プロトコルに基づく２次元パターンに従う走査を前記サンプルに適用して収集されたデータに基づき構築される、
請求項１１の画像処理方法。

【請求項13】

一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータを収集し、
前記データに基づいて第１画像及び第２画像を構築し、
前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲と前記マスク画像の画素値の範囲との間の差を小さくするように、前記第１画像の画素値の範囲と前記第２画像の画素値の範囲とマスク画像の画素値の範囲とを相対的に調整し、
前記第１画像及び前記第２画像のそれぞれに前記マスク画像を合成して第１合成画像及び第２合成画像を生成し、
前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求め、
前記複数の相互相関関数に基づいて相関係数を算出し、
前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行い、
前記レジストレーションの結果に基づいて前記第１画像と前記第２画像とのマージ画像を構築する、
走査型イメージング方法。

【請求項14】

前記データに基づいて第３画像を構築し、
前記第３画像の画素値の範囲と前記マージ画像の画素値の範囲とマスク画像の画素値の範囲との間の差を小さくするように、前記第３画像の画素値の範囲と前記マージ画像の画素値の範囲とマスク画像の画素値の範囲とを相対的に調整し、
前記第３画像及び前記マージ画像のそれぞれに当該マスク画像を合成して第３合成画像及び第４合成画像を生成し、
前記第３合成画像及び前記第４合成画像に基づいて複数の相互相関関数を求め、
当該複数の相互相関関数に基づいて相関係数を算出し、
当該相関係数に基づいて前記第３画像と前記第４画像とのレジストレーションを行い、
前記レジストレーションの結果に基づいて前記第３画像と前記マージ画像とをマージする、
請求項１３の走査型イメージング方法。

【請求項15】

前記一連のサイクルのうちの任意の２つのサイクルは、少なくとも１点で互いに交差し、
前記一連のサイクルのうちの第１サイクル群に従う光走査で収集された第１データに基づいて前記第１画像を構築し、
前記一連のサイクルのうち前記第１サイクル群と異なる第２サイクル群に従う光走査で収集された第２データに基づいて前記第２画像を構築する、
請求項１３又は１４の走査型イメージング方法。

【請求項16】

リサージュ関数に基づき予め設定された走査プロトコルに基づいて前記２次元パターンに従う光走査を行う、
請求項１５の走査型イメージング方法。

【請求項17】

サンプルの第１画像及び第２画像を記憶する記憶部と、
前記第１画像に応じた第１マスク画像と前記第２画像に応じた第２マスク画像とを生成するマスク画像生成部と、
前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲に基づいて、前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲と前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲との間の差を小さくするように、前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲を調整する範囲調整部と、
前記第１画像に前記第１マスク画像を合成して第１合成画像を生成し、且つ、前記第２画像に前記第２マスク画像を合成して第２合成画像を生成する合成画像生成部と、
前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求める第１算出部と、
前記複数の相互相関関数に基づいて相関係数を算出する第２算出部と、
前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行うレジストレーション部と
を含む、画像処理装置。

【請求項18】

一連のサイクルを含む２次元パターンに従う光走査によりサンプルから収集されたデータを受け付けるデータ受付部と、
前記データに基づいて前記第１画像及び前記第２画像を構築する画像構築部と
を更に含む、
請求項１７の画像処理装置。

【請求項19】

サンプルの第１画像及び第２画像を記憶する記憶部と、プロセッサとを含む画像処理装置を制御する方法であって、
前記プロセッサを、
前記第１画像に応じた第１マスク画像と前記第２画像に応じた第２マスク画像とを生成し、
前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲に基づいて、前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲と前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲との間の差を小さくするように、前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲を調整し、
前記第１画像に前記第１マスク画像を合成して第１合成画像を生成し、且つ、前記第２画像に前記第２マスク画像を合成して第２合成画像を生成し、
前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求め、
前記複数の相互相関関数に基づいて相関係数を算出し、及び、
前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行う
ように制御する、
画像処理装置の制御方法。

【請求項20】

一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータを収集する走査部と、
前記データに基づいて第１画像及び第２画像を構築する画像構築部と、
前記第１画像に応じた第１マスク画像と前記第２画像に応じた第２マスク画像とを生成するマスク画像生成部と、
前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲に基づいて、前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲と前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲との間の差を小さくするように、前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲を調整する範囲調整部と、
前記第１画像に前記第１マスク画像を合成して第１合成画像を生成し、且つ、前記第２画像に前記第２マスク画像を合成して第２合成画像を生成する合成画像生成部と、
前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求める第１算出部と、
前記複数の相互相関関数に基づいて相関係数を算出する第２算出部と、
前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行うレジストレーション部と、
前記レジストレーションの結果に基づいて前記第１画像と前記第２画像とのマージ画像を構築するマージ処理部と
を含む、走査型イメージング装置。

【請求項21】

前記走査部は、互いに異なる第１方向及び第２方向に光を偏向可能な偏向器を含み、前記第１方向に沿った偏向方向の変化を第１周期で繰り返しつつ前記第２方向に沿った偏向方向の変化を前記第１周期と異なる第２周期で繰り返すことにより前記２次元パターンに従う光走査を前記サンプルに適用する、
請求項２０の走査型イメージング装置。

【請求項22】

サンプルに光走査を適用してデータを収集する走査部と、プロセッサとを含む走査型イメージング装置を制御する方法であって、
前記走査部を、一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータを収集するように制御し、
前記プロセッサを、
前記データに基づいて第１画像及び第２画像を構築し、
前記第１画像に応じた第１マスク画像と前記第２画像に応じた第２マスク画像とを生成し、
前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲に基づいて、前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲と前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲との間の差を小さくするように、前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲を調整し、
前記第１画像に前記第１マスク画像を合成して第１合成画像を生成し、且つ、前記第２画像に前記第２マスク画像を合成して第２合成画像を生成し、
前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求め、
前記複数の相互相関関数に基づいて相関係数を算出し、
前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行い、及び、
前記レジストレーションの結果に基づいて前記第１画像と前記第２画像とのマージ画像を構築する
ように制御する、
走査型イメージング装置の制御方法。

【請求項23】

請求項１～１６、１９又は２２の方法をコンピュータに実行させるプログラム。

【請求項24】

請求項２３のプログラムが記録されたコンピュータ可読な非一時的記録媒体。

【発明の詳細な説明】

【技術分野】

【0001】

この発明は、画像処理方法、走査型イメージング方法、画像処理装置、その制御方法、走査型イメージング装置、その制御方法、プログラム、及び記録媒体に関する。

【背景技術】

【0002】

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

【0003】

光を利用した走査型イメージングの典型例として光コヒーレンストモグラフィ（ＯＣＴ；光干渉断層撮影法）が知られている。ＯＣＴは、光散乱媒質をマイクロメートルレベル又はそれ以下の分解能で画像化することが可能な技術であり、医用イメージングや非破壊検査などに応用されている。ＯＣＴは、低コヒーレンス干渉法に基づく技術であり、典型的には、光散乱媒質のサンプルへの深達性を担保するために近赤外光を利用する。

【0004】

例えば眼科画像診断ではＯＣＴ装置の普及が進んでおり、２次元的なイメージングだけでなく、３次元的なイメージング・構造解析・機能解析なども実用化され、診断の強力なツールとして広く利用されるに至っている。また、眼科分野では、走査型レーザー検眼鏡（ＳＬＯ）など、ＯＣＴ以外の走査型イメージングも利用されている。なお、近赤外光以外の波長帯の光（電磁波）や超音波を利用した走査型イメージングも知られている。

【0005】

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

【0006】

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

【0007】

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

【0008】

リサージュスキャンの特性として、２つのストリップは４箇所で重複する（交差する）。非特許文献１に記載された画像処理では、最大のストリップが初期基準ストリップとして採用され、初期基準ストリップに対して大きさ順にストリップが逐次にレジストレーションされ、合成される。レジストレーションでは、基準ストリップと他のストリップとの間の相対位置を求めるために相互相関関数が利用される。なお、ストリップは任意の形状を有するため、各ストリップを特定形状（例えば正方形状）の画像として扱うためのマスクを用いてストリップ間のオーバーラップ領域の相関演算が行われる（非特許文献１のAppendix Aを参照）。

【0009】

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

【先行技術文献】

【特許文献】

【0010】

【文献】特開２０１６－１７９１５号公報

【文献】特開２０１８－６８５７８号公報

【文献】特開２０１８－１４０００４号公報

【文献】特開２０１８－１４００４９号公報

【非特許文献】

【0011】

【文献】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

【発明の概要】

【発明が解決しようとする課題】

【0012】

この発明の１つの目的は、マスクを用いた画像レジストレーションの正確性を向上させることにある。

【課題を解決するための手段】

【0013】

幾つかの態様は、画像処理方法であって、サンプルの第１画像及び第２画像を準備し、マスク画像を準備し、前記第１画像の画素値の範囲と前記第２画像の画素値の範囲と前記マスク画像の画素値の範囲とを相対的に調整し、前記第１画像及び前記第２画像のそれぞれに前記マスク画像を合成して第１合成画像及び第２合成画像を生成し、前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求め、前記複数の相互相関関数に基づいて相関係数を算出し、前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行う。

【0014】

幾つかの態様において、画像処理方法は、前記サンプルの画像の画素値の範囲と前記マスク画像の画素値の範囲との間の差を小さくする。

【0015】

幾つかの態様において、画像処理方法は、前記サンプルの画像及び前記マスク画像の一方の画素値の範囲を他方の画素値の範囲に一致させる。

【0016】

幾つかの態様において、画像処理方法は、前記マスク画像の画素値の範囲に応じて前記サンプルの画像の画素値の範囲を正規化する。

【0017】

幾つかの態様において、前記マスク画像の画素値の範囲は閉区間［０，１］に含まれ、画像処理方法は、前記サンプルの画像の各画素の値を、当該画像における最大画素値で除算する。

【0018】

幾つかの態様において、前記マスク画像の画素値の範囲は閉区間［０，１］に含まれ、画像処理方法は、前記サンプルの画像の各画素の値を、当該画像の画素値の範囲の最大値で除算する。

【0019】

幾つかの態様において、前記マスク画像は、画素値が０又は１の二値画像である。

【0020】

幾つかの態様において、前記マスク画像は、前記サンプルの画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像である。

【0021】

幾つかの態様において、前記マスク画像は、第１マスク画像と第２マスク画像とを含み、画像処理方法は、前記第１マスク画像を前記第１画像に合成して前記第１合成画像を生成し、前記第２マスク画像を前記第２画像に合成して前記第２合成画像を生成する。

【0022】

幾つかの態様において、前記第１マスク画像及び前記第２マスク画像は、同じ寸法及び同じ形状を有する。

【0023】

幾つかの態様において、前記第１マスク画像は、前記第１画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であり、前記第２マスク画像は、前記第２画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像である。

【0024】

幾つかの態様において、前記第１画像及び前記第２画像のそれぞれは、一連のサイクルを含む２次元パターンに従う走査を前記サンプルに適用して収集されたデータに基づき構築される。

【0025】

幾つかの態様において、前記第１画像及び前記第２画像のそれぞれは、リサージュ関数に基づき予め設定された走査プロトコルに基づく２次元パターンに従う走査を前記サンプルに適用して収集されたデータに基づき構築される。

【0026】

幾つかの態様は、走査型イメージング方法であって、一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータを収集し、前記データに基づいて第１画像及び第２画像を構築し、前記第１画像の画素値の範囲と前記第２画像の画素値の範囲とマスク画像の画素値の範囲とを相対的に調整し、前記第１画像及び前記第２画像のそれぞれに前記マスク画像を合成して第１合成画像及び第２合成画像を生成し、前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求め、前記複数の相互相関関数に基づいて相関係数を算出し、前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行い、前記レジストレーションの結果に基づいて前記第１画像と前記第２画像とのマージ画像を構築する。

【0027】

幾つかの態様において、走査型イメージング方法は、前記データに基づいて第３画像を構築し、前記第３画像の画素値の範囲と前記マージ画像の画素値の範囲とマスク画像の画素値の範囲とを相対的に調整し、前記第３画像及び前記マージ画像のそれぞれに当該マスク画像を合成して第３合成画像及び第４合成画像を生成し、前記第３合成画像及び前記第４合成画像に基づいて複数の相互相関関数を求め、当該複数の相互相関関数に基づいて相関係数を算出し、当該相関係数に基づいて前記第３画像と前記第４画像とのレジストレーションを行い、前記レジストレーションの結果に基づいて前記第３画像と前記マージ画像とをマージする。

【0028】

幾つかの態様において、前記一連のサイクルのうちの任意の２つのサイクルは、少なくとも１点で互いに交差し、走査型イメージング方法は、前記一連のサイクルのうちの第１サイクル群に従う光走査で収集された第１データに基づいて前記第１画像を構築し、前記一連のサイクルのうち前記第１サイクル群と異なる第２サイクル群に従う光走査で収集された第２データに基づいて前記第２画像を構築する。

【0029】

幾つかの態様において、走査型イメージング方法は、リサージュ関数に基づき予め設定された走査プロトコルに基づいて前記２次元パターンに従う光走査を行う。

【0030】

幾つかの態様は、画像処理装置であって、サンプルの第１画像及び第２画像を記憶する記憶部と、前記第１画像に応じた第１マスク画像と前記第２画像に応じた第２マスク画像とを生成するマスク画像生成部と、前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲に基づいて前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲を調整する範囲調整部と、前記第１画像に前記第１マスク画像を合成して第１合成画像を生成し、且つ、前記第２画像に前記第２マスク画像を合成して第２合成画像を生成する合成画像生成部と、前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求める第１算出部と、前記複数の相互相関関数に基づいて相関係数を算出する第２算出部と、前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行うレジストレーション部とを含む。

【0031】

幾つかの態様において、画像処理装置は、一連のサイクルを含む２次元パターンに従う光走査によりサンプルから収集されたデータを受け付けるデータ受付部と、前記データに基づいて前記第１画像及び前記第２画像を構築する画像構築部とを更に含む。

【0032】

幾つかの態様は、サンプルの第１画像及び第２画像を記憶する記憶部と、プロセッサとを含む画像処理装置を制御する方法であって、前記プロセッサを、前記第１画像に応じた第１マスク画像と前記第２画像に応じた第２マスク画像とを生成し、前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲に基づいて前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲を調整し、前記第１画像に前記第１マスク画像を合成して第１合成画像を生成し、且つ、前記第２画像に前記第２マスク画像を合成して第２合成画像を生成し、前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求め、前記複数の相互相関関数に基づいて相関係数を算出し、及び、前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行うように制御する。

【0033】

幾つかの態様は、走査型イメージング装置であって、一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータを収集する走査部と、前記データに基づいて第１画像及び第２画像を構築する画像構築部と、前記第１画像に応じた第１マスク画像と前記第２画像に応じた第２マスク画像とを生成するマスク画像生成部と、前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲に基づいて前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲を調整する範囲調整部と、前記第１画像に前記第１マスク画像を合成して第１合成画像を生成し、且つ、前記第２画像に前記第２マスク画像を合成して第２合成画像を生成する合成画像生成部と、前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求める第１算出部と、前記複数の相互相関関数に基づいて相関係数を算出する第２算出部と、前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行うレジストレーション部と、前記レジストレーションの結果に基づいて前記第１画像と前記第２画像とのマージ画像を構築するマージ処理部とを含む。

【0034】

幾つかの態様において、前記走査部は、互いに異なる第１方向及び第２方向に光を偏向可能な偏向器を含み、前記第１方向に沿った偏向方向の変化を第１周期で繰り返しつつ前記第２方向に沿った偏向方向の変化を前記第１周期と異なる第２周期で繰り返すことにより前記２次元パターンに従う光走査を前記サンプルに適用する。

【0035】

幾つかの態様は、サンプルに光走査を適用してデータを収集する走査部と、プロセッサとを含む走査型イメージング装置を制御する方法であって、前記走査部を、一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータを収集するように制御し、前記プロセッサを、前記データに基づいて第１画像及び第２画像を構築し、前記第１画像に応じた第１マスク画像と前記第２画像に応じた第２マスク画像とを生成し、前記第１マスク画像の画素値の範囲及び前記第２マスク画像の画素値の範囲に基づいて前記第１画像の画素値の範囲及び前記第２画像の画素値の範囲を調整し、前記第１画像に前記第１マスク画像を合成して第１合成画像を生成し、且つ、前記第２画像に前記第２マスク画像を合成して第２合成画像を生成し、前記第１合成画像及び前記第２合成画像に基づいて複数の相互相関関数を求め、前記複数の相互相関関数に基づいて相関係数を算出し、前記相関係数に基づいて前記第１画像と前記第２画像とのレジストレーションを行い、及び、前記レジストレーションの結果に基づいて前記第１画像と前記第２画像とのマージ画像を構築するように制御する。

【0036】

幾つかの態様は、上記のいずれかの方法をコンピュータに実行させるプログラムである。

【0037】

幾つかの態様は、上記のいずれかの方法をコンピュータに実行させるプログラムが記録された、コンピュータ可読な非一時的記録媒体である。

【0038】

いずれか２以上の態様を少なくとも部分的に組み合わせることが可能である。また、いずれかの態様に対し、本開示に係る任意の事項を少なくとも部分的に組み合わせることが可能である。また、いずれか２以上の態様の少なくとも部分的な組み合わせに対し、本開示に係る任意の事項を少なくとも部分的に組み合わせることが可能である。

【発明の効果】

【0039】

例示的な態様によれば、マスクを用いた画像レジストレーションの正確性を向上させることが可能である。

【図面の簡単な説明】

【0040】

【図1】例示的な実施形態に係る走査型イメージング装置（眼科装置）の構成の一例を表す概略図である。

【図2】例示的な実施形態に係る走査型イメージング装置（眼科装置）の構成の一例を表す概略図である。

【図3】例示的な実施形態に係る走査型イメージング装置（眼科装置）の構成の一例を表す概略図である。

【図4A】例示的な実施形態に係る走査型イメージング装置（眼科装置）の構成の一例を表す概略図である。

【図4B】例示的な実施形態に係る走査型イメージング装置（眼科装置）の構成の一例を表す概略図である。

【図5】例示的な実施形態に係る走査型イメージング装置（眼科装置）が実行するリサージュスキャンのパターンの例を示す概略図である。

【図6】例示的な実施形態に係る走査型イメージング装置（眼科装置）の動作の一例を表すフローチャートである。

【図7A】例示的な実施形態に係る走査型イメージング装置（眼科装置）の動作の一例を説明するための図である。

【図7B】例示的な実施形態に係る走査型イメージング装置（眼科装置）の動作の一例を説明するための図である。

【図7C】例示的な実施形態に係る走査型イメージング装置（眼科装置）の動作の一例を説明するための図である。

【発明を実施するための形態】

【0041】

幾つかの例示的な実施形態に係る画像処理方法、走査型イメージング方法、画像処理装置、その制御方法、走査型イメージング装置、その制御方法、プログラム、及び記録媒体について、図面を参照しながら説明する。

【0042】

本明細書にて引用された文献に開示された技術や、それ以外の任意の公知技術を、例示的な実施形態に援用することが可能である。また、特に言及しない限り、「画像データ」とそれに基づく「画像」とを区別せず、被検眼の「部位」とその「画像」とを区別しないものとする。

【0043】

以下に説明する例示的な実施形態に係る走査型イメージング装置は、フーリエドメインＯＣＴ（例えば、スウェプトソースＯＣＴ）を利用して生体眼の眼底を計測することが可能な眼科装置である。実施形態に採用可能なＯＣＴのタイプは、スウェプトソースＯＣＴに限定されず、例えばスペクトラルドメインＯＣＴ又はタイムドメインＯＣＴであってもよい。

【0044】

幾つかの例示的な実施形態に係る走査型イメージング装置は、ＯＣＴ以外の走査型モダリティを利用可能であってよい。例えば、幾つかの例示的な実施形態は、ＳＬＯ等の任意の光走査型モダリティを採用することができる。

【0045】

また、幾つかの例示的な実施形態に係る走査型イメージング装置に適用可能な走査型モダリティは、光走査型モダリティには限定されない。例えば、幾つかの例示的な実施形態は、光以外の電磁波を利用した走査型モダリティや、超音波を利用した走査型モダリティなどを採用することができる。

【0046】

幾つかの例示的な実施形態に係る眼科装置は、ＯＣＴスキャンで収集されたデータの処理及び／又はＳＬＯで収集されたデータの処理に加え、他のモダリティで取得されたデータを処理可能であってよい。他のモダリティは、例えば、眼底カメラ、スリットランプ顕微鏡、及び眼科手術用顕微鏡のいずれかであってよい。幾つかの例示的な実施形態に係る眼科装置は、このようなモダリティの機能を有していてよい。

【0047】

ＯＣＴが適用される対象（サンプル）は眼底に限定されず、前眼部や硝子体など眼の任意の部位であってもよいし、眼以外の生体の部位や組織であってもよいし、生体以外の物体であってもよい。

【0048】

以下に説明する例示的な実施形態に係る眼科装置は、走査型イメージング装置の幾つかの態様、走査型イメージング装置を制御する方法の幾つかの態様、画像処理装置の幾つかの態様、画像処理装置を制御する方法の幾つかの態様、画像処理方法の幾つかの態様、及び、走査型イメージング方法の幾つかの態様を提供する。

【0049】

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

【0050】

本明細書に開示された要素の機能は、回路構成（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）、従来の回路構成、及びそれらの任意の組み合わせのいずれかを含む。プロセッサは、トランジスタ及び／又は他の回路構成を含む、処理回路構成又は回路構成とみなされる。本開示において、回路構成、ユニット、手段、又はこれらに類する用語は、開示された機能を実行するハードウェア、又は、開示された機能を実行するようにプログラムされたハードウェアである。ハードウェアは、本明細書に開示されたハードウェアであってよく、或いは、記載された機能を実行するようにプログラム及び／又は構成された既知のハードウェアであってもよい。ハードウェアが或るタイプの回路構成とみなされ得るプロセッサである場合、回路構成、ユニット、手段、又はこれらに類する用語は、ハードウェアとソフトウェアとの組み合わせであり、このソフトウェアはハードウェア及び／又はプロセッサを構成するために使用される。

【0051】

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

【0052】

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

【0053】

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

【0054】

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

【0055】

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

【0056】

ＬＣＤ等の表示デバイスの代わりに、例えば、発光素子アレイ、又は、発光素子とこれを移動する機構との組み合わせを用いてもよい。

【0057】

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

【0058】

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

【0059】

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

【0060】

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

【0061】

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

【0062】

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

【0063】

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

【0064】

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

【0065】

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

【0066】

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

【0067】

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

【0068】

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

【0069】

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

【0070】

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

【0071】

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

【0072】

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

【0073】

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

【0074】

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

【0075】

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

【0076】

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

【0077】

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

【0078】

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

【0079】

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

【0080】

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

【0081】

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

【0082】

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

【0083】

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

【0084】

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

【0085】

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

【0086】

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

【0087】

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

【0088】

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

【0089】

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

【0090】

このようなリサージュスキャンにおいて眼底Ｅｆに適用されるスキャンラインの分布の例を図５に示す。

【0091】

（狭義又は広義の）リサージュスキャンに含まれる任意のサイクル対（ペア）が１点以上（特に２点以上）で互いに交差するため、リサージュスキャンにおける任意のサイクル対から収集されたデータ対の間のレジストレーションを行うことができる。これにより、非特許文献１（又は２）に開示された画像構築手法及びモーションアーティファクト補正手法を利用することが可能となる。以下、特に言及しない限り、非特許文献１に記載された手法が適用される。

【0092】

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

【0093】

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

【0094】

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

【0095】

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

【0096】

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

【0097】

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

【0098】

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

【0099】

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

【0100】

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

【0101】

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

【0102】

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

【0103】

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

【0104】

画像構築部２２０は、プロセッサを含むハードウェアと、画像構築ソフトウェアとの協働によって実現される。

【0105】

〈データ処理部２３０〉
データ処理部２３０は、プロセッサを含み、被検眼Ｅの画像に対して各種のデータ処理を適用する。例えば、データ処理部２３０は、プロセッサを含むハードウェアと、データ処理ソフトウェアとの協働によって実現される。

【0106】

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

【0107】

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

【0108】

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

【0109】

このようなレジストレーションでは、基準ストリップと他のストリップとの間の相対位置を求めるために相互相関関数が利用されるが、ストリップは任意の形状を有するため、各ストリップを特定形状（例えば正方形状）の画像として扱うためのマスクを用いてストリップ間の相関演算が行われる（非特許文献１のAppendix Aを参照）。

【0110】

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

【0111】

この問題に対処するために、本実施形態では、図４Ｂに示す構成を有するデータ処理部２３０が採用される。本例のデータ処理部２３０は、マスク画像生成部２３１と、範囲調整部２３２と、合成画像生成部２３３と、相互相関関数算出部２３４と、相関係数算出部２３５と、レジストレーション部２３６と、マージ処理部２３７とを含む。以下に説明する例では、任意の２つのストリップを考慮し、その一方が基準ストリップ（reference strip）に指定され、この基準ストリップに対して他方のストリップ（対象ストリップ；registering strip）がレジストレーションされる。

【0112】

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

【0113】

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

【0114】

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

【0115】

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

【0116】

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

【0117】

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

【0118】

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

【0119】

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

【0120】

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

【0121】

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

【0122】

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

【0123】

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

【0124】

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

【0125】

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

【0126】

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

【0127】

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

【0128】

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

【0129】

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

【0130】

〈レジストレーション部２３６〉
レジストレーション部２３６は、相関係数算出部２３５により算出される相関係数に基づいて、ｘｙ方向（ラテラル方向）のレジストレーションが実行される。このレジストレーションは、非特許文献１の「rough lateral motion correction」（第１７８７頁）に相当する。

【0131】

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

【0132】

加えて、レジストレーション部２３６は、ｚ方向（軸方向、アキシャル方向）のモーションアーティファクトを除去するために、非特許文献１の「axial motion correction (rough axial motion correction 及び／又は fine axial motion correction」（第１７８９頁～第１７９０頁）に相当するレジストレーションを実行するように構成されてもよい。

【0133】

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

【0134】

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

【0135】

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

【0136】

〈動作〉
眼科装置１の動作について説明する。眼科装置１の動作の例を図６に示す。なお、患者ＩＤの入力、走査モードの設定（リサージュスキャンの指定）、固視標の提示、アライメント、フォーカス調整、ＯＣＴ光路長調整など、従来と同様の準備的な処理は、既になされたものとする。

【0137】

（Ｓ１：リサージュスキャン）
所定の走査開始トリガー信号を受けて、走査制御部２１１１は、被検眼Ｅ（眼底Ｅｆ）に対するＯＣＴスキャン（リサージュスキャン）の適用を開始する。

【0138】

走査開始トリガー信号は、例えば、所定の準備動作（アライメント、フォーカス調整、ＯＣＴ光路長調整など）が完了したことに対応して、又は、操作部２４２を用いて走査開始指示操作が行われたことに対応して生成される。

【0139】

走査制御部２１１１は、走査プロトコル２１２１（リサージュスキャンに対応したプロトコル）に基づいて光スキャナ４４及びＯＣＴユニット１００などを制御することによってリサージュスキャンを眼底Ｅｆに適用する。リサージュスキャンにより収集されたデータは画像構築部２２０に送られる。

【0140】

（Ｓ２：複数のストリップを構築）
画像構築部２２０は、前述した要領で、ステップＳ１で収集されたデータに基づいて複数のストリップを構築する。構築された複数のストリップはデータ処理部２３０に送られる。

【0141】

（Ｓ３：基準ストリップと対象ストリップを設定する）
データ処理部２３０は、ステップＳ２で構築された複数のストリップを寸法（面積等）に従って順序付けする。本例では、ステップＳ２で構築された複数のストリップの個数をＮ個（Ｎは２以上の整数）であるとする。また、順序付けにより指定された順序に従って、これらＮ個のストリップを、第１のストリップ、第２のストリップ、・・・、第Ｎのストリップと呼ぶ。また、第１～第Ｎのストリップのうちの任意のストリップを第ｎのストリップと呼ぶことがある（ｎ＝１、２、・・・、Ｎ）。このように、第１のストリップは最も大きなストリップであり、第２のストリップは２番目に大きなストリップであり、第ｎのストリップはｎ番目に大きなストリップであり、第Ｎのストリップは最も小さなストリップである。

【0142】

本例において、データ処理部２３０は、第１のストリップを初期基準ストリップに設定し、第２のストリップを初期対象ストリップに設定する。

【0143】

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

【0144】

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

【0145】

この段階では、マスク画像生成部２３１は、第１のストリップ及び第２のストリップにそれぞれ対応する２つのマスク画像を生成する。データ処理部２３０は、第１のストリップに対応する第１のマスク画像を初期基準マスク画像に設定し、第２のストリップに対応する第２のマスク画像を初期対象マスク画像に設定する。

【0146】

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

【0147】

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

【0148】

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

【0149】

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

【0150】

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

【0151】

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

【0152】

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

【0153】

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

【0154】

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

【0155】

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

【0156】

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

【0157】

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

【0158】

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

【0159】

（Ｓ９：基準ストリップと対象ストリップとの間のレジストレーション）
基準ストリップ及び対象ストリップに対して、レジストレーション部２３６は、ステップＳ８で算出された相関係数に基づきｘｙ方向のレジストレーション（rough lateral motion correction）を適用する。この段階では、第１のストリップｆ（ｒ）と第２のストリップとの間のレジストレーションが行われる。

【0160】

本例では、ステップＳ８で算出された相関係数ρ（ｒ´）のピークを検出することにより、基準ストリップと対象ストリップとの間の相対的なシフト量（Δｘ、Δｙ）が求められる。

【0161】

更に、レジストレーション部２３６は、前述した「fine lateral motion correction」や「axial motion correction」などのレジストレーションを基準ストリップ及び対象ストリップに対して適用する。

【0162】

（Ｓ１０：基準ストリップと対象ストリップのマージ画像を構築）
マージ処理部２３７は、ステップＳ９のレジストレーションにより相対位置調整がなされた基準ストリップと対象ストリップとのマージ画像を構築する。図７Ｃの画像３３０は、第１のストリップと第２のストリップとのマージ画像の例を示す。

【0163】

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

【0164】

Ｎが３以上の場合において、ステップＳ１０で第１のストリップと第２のストリップとのマージ画像が作成された場合（Ｓ１１：Ｎｏ）、処理はステップＳ３に戻る。ステップＳ３では、第１のストリップと第２のストリップとのマージ画像が新たな基準ストリップに設定され、且つ、第３のストリップが新たな対象ストリップに設定される。新たな基準ストリップと新たな対象ストリップに基づきステップＳ４～Ｓ１０の処理を実行することにより、新たな基準ストリップと新たな対象ストリップとのマージ画像が得られる。この新たなマージ画像は、第１～第３のストリップのマージ画像である。このような一連の処理をＮ個のストリップに対して逐次に適用することにより、Ｎ個のストリップの全てのマージ画像が得られる（Ｓ１１：Ｙｅｓ）。

【0165】

（Ｓ１２：最終的なマージ画像を保存する）
以上に説明した繰り返し処理により最終的に得られたマージ画像は、モーションアーティファクトが補正された画像であり、且つ、ステップＳ１のリサージュスキャンの適用範囲全体を表現した画像である。主制御部２１１は、この最終的なマージ画像を記憶部２１２（及び／又は、他の記憶装置）に保存する。

【0166】

データ処理部２３０は、Ｎ個のストリップに基づくレジストレーションの結果を利用して、ステップＳ１で収集されたデータ（３次元データ）、及び／又は、この３次元データから画像構築部２２０が構築した３次元画像データのレジストレーションを行うことができる。このレジストレーションは、リサージュスキャンの定義座標系を３次元直交座標系（ｘｙｚ座標系）に変換する処理を含む。この座標変換は、非特許文献１に記載された手法における「remapping」に相当する。主制御部２１１は、このようにして得られた３次元画像データを、最終的なマージ画像とともに又はその代わりに、記憶部２１２（及び／又は、他の記憶装置）に保存することができる（エンド）。

【0167】

主制御部２１１は、最終的なマージ画像を表示部２４１（及び／又は、他の表示装置）に表示させることができる。また、データ処理部２３０は、上記座標変換で構築された３次元画像データから任意のレンダリング画像を作成することができる。主制御部２１１は、このレンダリング画像を表示部２４１（及び／又は、他の表示装置）に表示させることができる。

【0168】

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

【0169】

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

【0170】

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

【0171】

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

【0172】

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

【0173】

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

【0174】

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

【0175】

〈作用、効果等〉
例示的な実施形態の幾つかの特徴について説明し、それらにより奏される幾つかの作用及び幾つかの効果について説明する。

【0176】

上記の例示的な実施形態に係る眼科装置１は、走査型イメージング装置の幾つかの態様を提供する。それに加えて、上記の例示的な実施形態に係る眼科装置１は、走査型イメージング装置を制御する方法の幾つかの態様、画像処理装置の幾つかの態様、画像処理装置を制御する方法の幾つかの態様、画像処理方法の幾つかの態様、及び、走査型イメージング方法の幾つかの態様を提供する。

【0177】

なお、走査型イメージング装置を制御する方法の幾つかの態様、画像処理装置の幾つかの態様、画像処理装置を制御する方法の幾つかの態様、画像処理方法の幾つかの態様、及び、走査型イメージング方法の幾つかの態様のいずれかにおいて説明された事項を、上記の例示的な実施形態に係る眼科装置１（より一般に、走査型イメージング装置のいずれかの態様）に適用することが可能である。

【0178】

さて、例示的な実施形態に係る走査型イメージング装置（眼科装置１）は、走査部と、画像構築部２２０と、マスク画像生成部２３１と、範囲調整部２３２と、合成画像生成部２３３と、相互相関関数算出部２３４（第１算出部）と、相関係数算出部２３５（第２算出部）と、レジストレーション部２３６と、マージ処理部２３７とを含む。

【0179】

走査部は、サンプル（被検眼Ｅ、眼底Ｅｆ）に光走査を適用してデータを収集するように構成されている。この光走査は、一連のサイクルを含む２次元パターンに従って実行される。

【0180】

幾つかの例示的な態様において、光走査の２次元パターンに含まれる一連のサイクルのうちの任意の２つのサイクルは、少なくとも１点で互いに交差していてよい。

【0181】

更に、幾つかの例示的な態様において、サンプルに適用される光走査は、リサージュ関数に基づき予め設定された走査プロトコルに基づく２次元パターンに従って実行されてよい。このような光走査の一例が、上記の例示的な実施形態におけるリサージュスキャンである。

【0182】

幾つかの例示的な態様において、走査部は、互いに異なる第１方向（ｘ方向）及び第２方向（ｙ方向）に光を偏向可能な偏向器（光スキャナ４４）を含み、第１方向に沿った偏向方向の変化を第１周期で繰り返しつつ第２方向に沿った偏向方向の変化を第１周期と異なる第２周期で繰り返すことによって、２次元パターンに従う光走査をサンプルに適用するように構成されていてよい。

【0183】

上記の例示的な実施形態において、走査部は、ＯＣＴユニット１００と、測定アームを構成する眼底カメラユニット２内の要素（リトロリフレクタ４１、ＯＣＴ合焦レンズ４３、光スキャナ４４、対物レンズ２２等）とを含み、被検眼ＥにＯＣＴスキャンを適用してＯＣＴデータを収集する。

【0184】

画像構築部２２０は、走査部により収集されたデータに基づいて第１画像及び第２画像を構築するように構成されている。上記の例示的な実施形態における基準ストリップ及び対象ストリップの組み合わせは、第１画像及び第２画像の組み合わせの一例である。

【0185】

幾つかの例示的な態様において、光走査の２次元パターンに含まれる一連のサイクルのうちの任意の２つのサイクルが、少なくとも１点で互いに交差する場合、画像構築部２２０は、一連のサイクルのうちの第１サイクル群に従う光走査で収集された第１データに基づいて第１画像を構築し、且つ、一連のサイクルのうち第１サイクル群と異なる第２サイクル群に従う光走査で収集された第２データに基づいて第２画像を構築するように構成されていてよい。

【0186】

例えば、上記の例示的な実施形態におけるリサージュスキャンのように、光走査の２次元パターンがリサージュ関数に基づき予め設定された走査プロトコルに従っている場合、画像構築部２２０は、この２次元パターンに含まれる一連のサイクルのうちの第１サイクル群に従う光走査で収集された第１データに基づいて第１画像を構築し、且つ、一連のサイクルのうち第１サイクル群と異なる第２サイクル群に従う光走査で収集された第２データに基づいて第２画像を構築するように構成されていてよい。

【0187】

上記の例示的な実施形態では、一連のサイクルに従う光走査から得られたボリュームのうち、比較的大きな動きが介在しない（第１サイクル群に対応する）第１のサブボリュームの正面プロジェクション画像が第１画像（第１のストリップ）に設定され、比較的大きな動きが介在しない（第１サイクル群とは異なる第２サイクル群に対応する）第２のサブボリュームの正面プロジェクション画像が第２画像（第２のストリップ）に設定される。

【0188】

画像構築部２２０により構築された第１画像及び第２画像について、マスク画像生成部２３１は、第１画像に応じた第１マスク画像と、第２画像に応じた第２マスク画像とを生成するように構成されている。なお、ここに言う「第１マスク画像」及び「第２マスク画像」は、上記の例示的な実施形態において複数のストリップに付与された順序に対応する「第１のマスク画像」及び「第２のマスク画像」とは異なる。

【0189】

幾つかの例示的な態様において、マスク画像（第１マスク画像、第２マスク画像）の画素値の範囲は、閉区間［０，１］に含まれていてよい。

【0190】

更に、幾つかの例示的な態様において、マスク画像（第１マスク画像、第２マスク画像）は、画素値が０又は１の二値画像であってよい。

【0191】

例えば、上記の例示的な実施形態のように、マスク画像（第１マスク画像、第２マスク画像）は、サンプルの画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であってよい。

【0192】

幾つかの例示的な態様において、第１マスク画像及び第２マスク画像は、同じ寸法及び同じ形状を有していてよい。

【0193】

範囲調整部２３２は、第１マスク画像の画素値の範囲及び第２マスク画像の画素値の範囲に基づいて第１画像の画素値の範囲及び第２画像の画素値の範囲を調整するように構成されている。

【0194】

幾つかの例示的な態様において、範囲調整部２３２は、サンプルの画像（第１画像、第２画像）の画素値の範囲とマスク画像（第１マスク画像、第２マスク画像）の画素値の範囲との間の差を小さくするように構成されていてよい。例えば、上記の例示的な実施形態のように、範囲調整部２３２は、サンプルの画像（第１画像、第２画像）の画素値の範囲とマスク画像（第１マスク画像、第２マスク画像）の画素値の範囲とを一致させるように構成されていてよい。

【0195】

更に、幾つかの例示的な態様において、範囲調整部２３２は、サンプルの画像（第１画像、第２画像）及びマスク画像（第１マスク画像、第２マスク画像）の一方の画素値の範囲を他方の画素値の範囲に一致させるように構成されてよい。例えば、上記の例示的な実施形態のように、範囲調整部２３２は、サンプルの画像（第１画像、第２画像）の画素値の範囲をマスク画像（第１マスク画像、第２マスク画像）の画素値の範囲に一致させるように構成されていてよい。すなわち、範囲調整部２３２は、サンプルの画像（第１画像、第２画像）の画素値の範囲と、マスク画像（第１マスク画像、第２マスク画像）の画素値の範囲とを一致させるために、サンプルの画像（第１画像、第２画像）の画素値の範囲を変更するように構成されてよい。

【0196】

幾つかの例示的な態様において、範囲調整部２３２は、マスク画像（第１マスク画像、第２マスク画像）の画素値の範囲に応じてサンプルの画像（第１画像、第２画像）の画素値の範囲を正規化（規格化）するように構成されていてよい。

【0197】

例えば、マスク画像（第１マスク画像、第２マスク画像）の画素値の範囲が閉区間［０，１］に含まれている場合において、範囲調整部２３２は、サンプルの画像（第１画像、第２画像）の各画素の値を、当該画像における最大画素値で除算することによって、当該画像の画素値の範囲を正規化するように構成されていてよい。

【0198】

他の例において、マスク画像（第１マスク画像、第２マスク画像）の画素値の範囲が閉区間［０，１］に含まれている場合において、範囲調整部２３２は、サンプルの画像（第１画像、第２画像）の各画素の値を、当該画像の画素値の範囲の最大値で除算することによって、当該画像の画素値の範囲を正規化するように構成されていてよい。

【0199】

合成画像生成部２３３は、サンプルの第１画像に第１マスク画像を合成して第１合成画像を生成し、且つ、第２画像に第２マスク画像を合成して第２合成画像を生成するように構成されている。

【0200】

幾つかの例示的な態様において、第１マスク画像及び第２マスク画像は、同じ寸法及び同じ形状を有していてよい。更に、第１マスク画像は、第１画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であってよく、且つ、第２マスク画像は、第２画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であってよい。

【0201】

上記の例示的な実施形態では、第１マスク画像と第２マスク画像とは、同じ寸法及び同じ形状を有する。合成画像生成部２３３は、任意形状の第１画像及び任意形状の第２画像のそれぞれを、これらマスク画像と同じ寸法及び同じ形状の埋め込み画像に変換するように構成されている。更に、合成画像生成部２３３は、第１画像の埋め込み画像と第１マスク画像とを合成し、且つ、第２画像の埋め込み画像と第２マスク画像とが合成するように構成されている。

【0202】

或いは、上記の例示的な実施形態において、合成画像生成部２３３は、任意形状の第１画像及び任意形状の第２画像を、互いに同じ寸法及び同じ形状の埋め込み画像にそれぞれ変換するように構成されていてよい。また、マスク画像生成部２３１は、これらサンプルの画像と同じ寸法及び同じ形状を有する第１マスク画像と第２マスク画像をそれぞれ生成するように構成されていてよい。更に、合成画像生成部２３３は、第１画像の埋め込み画像と第１マスク画像とを合成し、且つ、第２画像の埋め込み画像と第２マスク画像とが合成するように構成されていてよい。

【0203】

相互相関関数算出部２３４は、合成画像生成部２３３により生成された第１合成画像及び第２合成画像に基づいて複数の相互相関関数を求めるように構成されている。

【0204】

相関係数算出部２３５は、相互相関関数算出部２３４により求められた複数の相互相関関数に基づいて相関係数を算出するように構成されている。

【0205】

レジストレーション部２３６は、相関係数算出部２３５により算出された相関係数に基づいて第１画像と第２画像とのレジストレーションを行うように構成されている。

【0206】

マージ処理部２３７は、レジストレーション部２３６により実行されたレジストレーションの結果に基づいて第１画像と第２画像とのマージ画像を構築するように構成されている。

【0207】

上記の例示的な実施形態と同様に、幾つかの例示的な態様において、相互相関関数算出部２３４（第１算出部）、相関係数算出部２３５（第２算出部）、レジストレーション部２３６、及びマージ処理部２３７のそれぞれが実行する処理及びそのための構成（ハードウェア構成、ソフトウェア構成）は、非特許文献１に記載された手法と同様であってよい。ただし、他の手法を採用することを排除するものではない。

【0208】

上記の例示的な実施形態に係る眼科装置１により提供可能な画像処理方法の幾つかの態様を説明する。幾つかの例示的な態様において、画像処理方法は、以下に説明する工程を含んでいてよい。

【0209】

（画像準備工程）
サンプルの第１画像及び第２画像が準備される。例えば、通信回線や記録媒体を介してサンプルの第１画像及び第２画像が受け付けられる。或いは、サンプルに光走査を適用してデータが収集され、第１画像及び第２画像が構築される。第１画像及び第２画像のそれぞれは、一連のサイクルを含む２次元パターンに従う走査をサンプルに適用して収集されたデータに基づき構築される。例えば、第１画像及び第２画像のそれぞれは、リサージュ関数に基づき予め設定された走査プロトコルに基づく２次元パターンに従う走査をサンプルに適用して収集されたデータに基づき構築される。

【0210】

（マスク準備工程）
マスク画像が準備される。例えば、通信回線や記録媒体を介してマスク画像が受け付けられる。或いは、（サンプルの画像に基づいて）マスク画像が生成される。マスク画像の画素値の範囲は閉区間［０，１］に含まれていてよい。典型的には、マスク画像は、画素値が０又は１の二値画像であってよい。例えば、マスク画像は、サンプルの画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であってよい。マスク画像は、第１マスク画像と第２マスク画像とを含んでいてよい。第１マスク画像及び第２マスク画像は、同じ寸法及び同じ形状を有していてよい。第１マスク画像は、サンプルの第１画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であってよく、且つ、第２マスク画像は、サンプルの第２画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であってよい。

【0211】

（画像値範囲調整工程）
第１画像の画素値の範囲と第２画像の画素値の範囲とマスク画像の画素値の範囲とが相対的に調整される。画素値の範囲の相対的な調整は、典型的には、サンプルの画像（第１画像、第２画像）の画素値の範囲とマスク画像の画素値の範囲との間の差を小さくするように行われる。例えば、画素値の範囲の相対的な調整は、サンプルの画像（第１画像、第２画像）及びマスク画像の一方の画素値の範囲を他方の画素値の範囲に一致させるように行われる。幾つかの例示的な態様において、画素値の範囲の相対的な調整は、マスク画像の画素値の範囲に応じてサンプルの画像（第１画像、第２画像）の画素値の範囲を正規化する処理を含んでいてよい。マスク画像の画素値の範囲が閉区間［０，１］に含まれる場合において、正規化は、サンプルの画像（第１画像、第２画像）の各画素の値を、当該画像における最大画素値で除算する処理を含んでいてよい。或いは、マスク画像の画素値の範囲が閉区間［０，１］に含まれる場合において、正規化は、サンプルの画像（第１画像、第２画像）の各画素の値を、当該画像の画素値の範囲の最大値で除算する処理を含んでいてよい。

【0212】

（合成画像生成工程）
画像値範囲調整工程の後、第１画像及び第２画像のそれぞれにマスク画像が合成されて第１合成画像及び第２合成画像が生成される。マスク画像が、第１マスク画像と第２マスク画像とを含む場合において、合成画像生成工程は、第１マスク画像を第１画像に合成して第１合成画像を生成する工程と、第２マスク画像を第２画像に合成して第２合成画像を生成する工程とを含んでいてよい。ここで、第１マスク画像及び第２マスク画像は、同じ寸法及び同じ形状を有していてもよい。また、第１マスク画像は、第１画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であってよく、且つ、第２マスク画像は、第２画像の定義域に対応する領域内の画素の値が１であり、他の画素の値が０である二値画像であってよい。

【0213】

（相互相関関数算出工程）
合成画像生成工程によって生成された第１合成画像及び第２合成画像に基づいて複数の相互相関関数が求められる。この工程は、典型的には、非特許文献１に記載された手法に従って行われるが、これに限定されるものではない。

【0214】

（相関係数算出工程）
相互相関関数算出工程において求められた複数の相互相関関数に基づいて相関係数が算出される。この工程は、典型的には、非特許文献１に記載された手法に従って行われるが、これに限定されるものではない。

【0215】

（レジストレーション工程）
相関係数算出工程において算出された相関係数に基づいてサンプルの第１画像と第２画像との間のレジストレーションが行われる。この工程は、典型的には、非特許文献１に記載された手法に従って行われるが、これに限定されるものではない。

【0216】

幾つかの例示的な態様において、画像処理方法は、マージ画像構築工程を更に含んでいてよい。マージ画像構築工程では、レジストレーション工程において行われたレジストレーションの結果（例えば、サンプルの第１画像と第２画像との間の相対的シフト量）に基づいて第１画像と第２画像とのマージ画像が構築される。

【0217】

幾つかの例示的な態様において、画像処理方法は、上記の例示的な実施形態において説明された工程のいずれかを更に含んでいてもよい。

【0218】

上記の例示的な実施形態に係る眼科装置１により提供可能な走査型イメージング方法の幾つかの態様を説明する。幾つかの例示的な態様において、走査型イメージング方法は、以下に説明する工程を含んでいてよい。なお、幾つかの例示的な態様において、走査型イメージング方法は、上記の例示的な実施形態において説明された工程のいずれか、及び／又は、上記の画像処理方法の例示的な態様で説明された工程のいずれかを更に含んでいてもよい。

【0219】

（走査工程）
一連のサイクルを含む２次元パターンに従う光走査がサンプルに適用され、データが収集される。典型的には、一連のサイクルのうちの任意の２つのサイクルは、少なくとも１点で互いに交差する。例えば、２次元パターンに従う光走査は、リサージュ関数に基づき予め設定された走査プロトコルに基づいて行われる。

【0220】

（画像構築工程）
走査工程で収集されたデータに基づいて第１画像及び第２画像が構築される。一連のサイクルのうちの任意の２つのサイクルが少なくとも１点で互いに交差する場合（例えば、リサージュ関数に基づき予め設定された走査プロトコルに基づいて２次元パターンに従う光走査が行われる場合）、一連のサイクルのうちの第１サイクル群に従う光走査で収集された第１データに基づいて第１画像が構築されてよく、且つ、一連のサイクルのうち第１サイクル群と異なる第２サイクル群に従う光走査で収集された第２データに基づいて第２画像が構築されてよい。

【0221】

（画像値範囲調整工程）
第１画像の画素値の範囲と第２画像の画素値の範囲とマスク画像の画素値の範囲とが相対的に調整される。

【0222】

（合成画像生成工程）
画像値範囲調整工程の後、第１画像及び第２画像のそれぞれにマスク画像が合成されて第１合成画像及び第２合成画像が生成される。

【0223】

（相互相関関数算出工程）
合成画像生成工程で生成された第１合成画像及び第２合成画像に基づいて複数の相互相関関数が求められる。

【0224】

（相関係数算出工程）
相互相関関数算出工程で求められた複数の相互相関関数に基づいて相関係数が算出される。

【0225】

（レジストレーション工程）
相関係数算出工程で算出された相関係数に基づいてサンプルの第１画像と第２画像との間のレジストレーションが行われる。

【0226】

（マージ画像構築工程）
レジストレーションの結果に基づいてサンプルの第１画像と第２画像とのマージ画像が構築される。

【0227】

幾つかの例示的な態様において、走査型イメージング方法は、以下の付加的工程を含んでいてよい。

【0228】

第１の付加的工程において、走査工程においてサンプルから収集されたデータに基づいて第３画像が構築される。上記の例示的な実施形態では、第３のストリップが構築される。

【0229】

第２の付加的工程において、第３画像の画素値の範囲と、マージ画像構築工程で構築されたマージ画像の画素値の範囲と、マスク画像の画素値の範囲とが、相対的に調整される。上記の例示的な実施形態では、第３のストリップと、第３のストリップに対応するマスク画像と、第１のストリップ及び第２のストリップのマージ画像と、このマージ画像に対応するマスク画像とについて、画素値範囲調整が適用される。

【0230】

第３の付加的工程において、第３画像及びマージ画像のそれぞれに当該マスク画像を合成して第３合成画像及び第４合成画像が生成される。上記の例示的な実施形態では、第３のストリップとこれに対応するマスク画像との合成画像が生成され、且つ、第１のストリップ及び第２のストリップのマージ画像とこれに対応するマスク画像との合成画像が生成される。

【0231】

第４の付加的工程において、第３合成画像及び第４合成画像に基づいて複数の相互相関関数が求められる。上記の例示的な実施形態では、第３のストリップとこれに対応するマスク画像との合成画像と、上記マージ画像とこれに対応するマスク画像との合成画像とに基づいて、複数の相互相関関数が求められる。

【0232】

第５の付加的工程において、第４の付加的工程で求められた複数の相互相関関数に基づいて相関係数が算出される。上記の例示的な実施形態では、第３のストリップに基づく合成画像及びマージ画像に基づく合成画像を用いて求められた複数の相互相関関数に基づいて、相関係数が算出される。

【0233】

第６の付加的工程において、第５の付加的工程で算出された相関係数に基づいて第３画像と第４画像とのレジストレーションが行われる。上記の例示的な実施形態では、第３のストリップに基づく合成画像及びマージ画像に基づく合成画像を用いて算出された相関係数に基づいて、このマージ画像（第１のストリップ及び第２のストリップのマージ画像）と第３のストリップとの間のレジストレーションが行われる。

【0234】

第７の付加的工程において、第６の付加的工程で実行されたレジストレーションの結果に基づいて第３画像とマージ画像とがマージされる。上記の例示的な実施形態では、マージ画像と第３のストリップとのレジストレーションの結果（相対的なシフト量）に基づいて、このマージ画像（第１のストリップ及び第２のストリップのマージ画像）に対して第３のストリップがマージされる。

【0235】

なお、ここでは、第１のストリップ及び第２のストリップのマージ画像と、第３のストリップとに対する処理を例として説明したが、一般に、上記の例示的な実施形態で説明した通り、第１～第ｎのストリップのマージ画像と、第ｎ＋１のストリップとに対する処理も、同じ要領で実行される（ｎ＝２、３、・・・、Ｎ－１）。

【0236】

上記の例示的な実施形態に係る眼科装置１により提供可能な画像処理装置の幾つかの態様を説明する。幾つかの例示的な態様において、画像処理装置は、以下に説明する要素を含んでいてよい。なお、幾つかの例示的な態様において、画像処理装置は、上記の例示的な実施形態において説明された要素のいずれか、上記の画像処理方法の例示的な態様で説明された工程を実行するように構成された要素のいずれか、及び、上記の走査型イメージング方法の例示的な態様で説明された工程を実行するように構成された要素のいずれか、のうちのいずれかを更に含んでいてもよい。

【0237】

例示的な態様に係る画像処理装置は、記憶部と、マスク画像生成部と、範囲調整部と、合成画像生成部と、第１算出部と、第２算出部と、レジストレーション部とを含む。

【0238】

記憶部は、サンプルの第１画像及び第２画像を記憶する。上記の例示的な実施形態において、この記憶部は記憶部２１２に相当する。

【0239】

マスク画像生成部は、第１画像に応じた第１マスク画像と第２画像に応じた第２マスク画像とを生成するように構成される。上記の例示的な実施形態において、マスク画像生成部はマスク画像生成部２３１に相当する。

【0240】

範囲調整部は、第１マスク画像の画素値の範囲及び第２マスク画像の画素値の範囲に基づいて第１画像の画素値の範囲及び第２画像の画素値の範囲を調整するように構成される。上記の例示的な実施形態において、範囲調整部は範囲調整部２３２に相当する。

【0241】

合成画像生成部は、第１画像に第１マスク画像を合成して第１合成画像を生成し、且つ、第２画像に第２マスク画像を合成して第２合成画像を生成するように構成される。上記の例示的な実施形態において、合成画像生成部は合成画像生成部２３３に相当する。

【0242】

第１算出部は、第１合成画像及び第２合成画像に基づいて複数の相互相関関数を求めるように構成される。上記の例示的な実施形態において、第１算出部は相互相関関数算出部２３４に相当する。

【0243】

第２算出部は、複数の相互相関関数に基づいて相関係数を算出するように構成される。上記の例示的な実施形態において、第２算出部は相関係数算出部２３５に相当する。

【0244】

レジストレーション部は、相関係数に基づいて第１画像と第２画像とのレジストレーションを行うように構成される。上記の例示的な実施形態において、レジストレーション部はレジストレーション部２３６に相当する。

【0245】

幾つかの例示的な態様において、画像処理装置は、マージ処理部を更に含んでいてよい。マージ処理部は、レジストレーションの結果に基づいて第１画像と第２画像とのマージ画像を構築するように構成される。

【0246】

幾つかの例示的な態様において、画像処理装置は、データ受付部と、画像構築部とを更に含んでいてよい。

【0247】

データ受付部は、一連のサイクルを含む２次元パターンに従う光走査によりサンプルから収集されたデータを受け付けるように構成される。データ受付部は、例えば、通信回線又は記録媒体を介してデータを受け付けるように構成されていてよい。上記の例示的な実施形態において、データ受付部は、眼科装置１に設けられた通信デバイスに相当する。

【0248】

画像構築部は、データ受付部により受け付けられたデータに基づいて第１画像及び第２画像を構築する。上記の例示的な実施形態において、画像構築部は画像構築部２２０に相当する。画像構築部により構築された第１画像及び第２画像は記憶部に記憶される。

【0249】

幾つかの例示的な態様において、データ受付部は、一連のサイクルを含む２次元パターンに従う光走査によりサンプルから収集されたデータに基づき構築された第１画像及び第２画像を受け付けるように構成されてよい。この場合、データ受付部により受け付けられた第１画像及び第２画像が記憶部に記憶される。

【0250】

幾つかの例示的な態様において、データ受付部は、一連のサイクルを含む２次元パターンに従う光走査によりサンプルから収集されたデータに基づき構築された３次元画像データを受け付けるように構成されてよい。この場合、画像構築部は、データ受付部により受け付けられた３次元画像データに基づいて第１画像及び第２画像を構築するように構成される。画像構築部により構築された第１画像及び第２画像は記憶部に記憶される。

【0251】

上記の例示的な実施形態に係る眼科装置１により提供可能な、画像処理装置の制御方法の幾つかの態様を説明する。幾つかの例示的な態様において、画像処理装置の制御方法は、以下に説明する工程を含んでいてよい。なお、幾つかの例示的な態様において、画像処理装置の制御方法は、上記の例示的な実施形態において説明された工程のいずれか、上記の画像処理方法の例示的な態様で説明された工程のいずれか、上記の走査型イメージング方法の例示的な態様で説明された工程のいずれか、及び、上記の画像処理装置の例示的な態様で説明された要素により実行される工程のいずれか、のうちのいずれかを更に含んでいてもよい。

【0252】

例示的な態様は、サンプルの第１画像及び第２画像を記憶する記憶部と、プロセッサとを含む画像処理装置を制御する方法であって、このプロセッサが以下の工程を実行するように制御を行う。なお、上記の例示的な実施形態において、記憶部は記憶部２１２に相当し、プロセッサは少なくともデータ処理部２３０に相当する。

【0253】

プロセッサは、第１画像に応じた第１マスク画像と第２画像に応じた第２マスク画像とを生成するように制御される。例えば、上記の例示的な実施形態において、データ処理部２３０は、マスク画像生成部２３１として動作するように制御される。

【0254】

また、プロセッサは、第１マスク画像の画素値の範囲及び第２マスク画像の画素値の範囲に基づいて第１画像の画素値の範囲及び第２画像の画素値の範囲を調整するように制御される。例えば、上記の例示的な実施形態において、データ処理部２３０は、範囲調整部２３２として動作するように制御される。

【0255】

また、プロセッサは、第１画像に第１マスク画像を合成して第１合成画像を生成し、且つ、第２画像に第２マスク画像を合成して第２合成画像を生成するように制御される。例えば、上記の例示的な実施形態において、データ処理部２３０は、合成画像生成部２３３として動作するように制御される。

【0256】

また、プロセッサは、第１合成画像及び第２合成画像に基づいて複数の相互相関関数を求めるように制御される。例えば、上記の例示的な実施形態において、データ処理部２３０は、相互相関関数算出部２３４として動作するように制御される。

【0257】

また、プロセッサは、複数の相互相関関数に基づいて相関係数を算出するように制御される。例えば、上記の例示的な実施形態において、データ処理部２３０は、相関係数算出部２３５として動作するように制御される。

【0258】

また、プロセッサは、相関係数に基づいて第１画像と第２画像とのレジストレーションを行うように制御される。例えば、上記の例示的な実施形態において、データ処理部２３０は、レジストレーション部２３６として動作するように制御される。

【0259】

幾つかの例示的な実施形態において、プロセッサは、レジストレーションの結果に基づいて第１画像と第２画像とのマージ画像を構築するように制御されてよい。例えば、上記の例示的な実施形態において、データ処理部２３０は、マージ処理部２３７として動作するように制御される。

【0260】

上記の例示的な実施形態に係る眼科装置１により提供可能な、走査型イメージング装置の制御方法の幾つかの態様を説明する。幾つかの例示的な態様において、走査型イメージング装置の制御方法は、以下に説明する工程を含んでいてよい。なお、幾つかの例示的な態様において、走査型イメージング装置の制御方法は、上記の例示的な実施形態において説明された工程のいずれか、上記の画像処理方法の例示的な態様で説明された工程のいずれか、上記の走査型イメージング方法の例示的な態様で説明された工程のいずれか、上記の画像処理装置の例示的な態様で説明された要素により実行される工程のいずれか、及び、上記の画像処理装置の例示的な態様で説明された工程のいずれか、のうちのいずれかを更に含んでいてもよい。

【0261】

例示的な態様は、サンプルに光走査を適用してデータを収集する走査部と、プロセッサとを含む走査型イメージング装置を制御する方法であって、走査部及びプロセッサが以下の工程を実行するように制御を行う。なお、上記の例示的な実施形態において、走査部は、ＯＣＴユニット１００と、測定アームを構成する眼底カメラユニット２内の要素（リトロリフレクタ４１、ＯＣＴ合焦レンズ４３、光スキャナ４４、対物レンズ２２等）とを含む。また、プロセッサは、画像構築部２２０及びデータ処理部２３０に相当する。

【0262】

走査部は、一連のサイクルを含む２次元パターンに従う光走査をサンプルに適用してデータを収集するように制御される。

【0263】

また、プロセッサ（上記の例示的な実施形態では画像構築部２２０）は、走査部により収集されたデータに基づいて第１画像及び第２画像を構築するように制御される。

【0264】

また、プロセッサ（上記の例示的な実施形態ではデータ処理部２３０）は、第１画像に応じた第１マスク画像と第２画像に応じた第２マスク画像とを生成するように制御される。

【0265】

また、プロセッサ（上記の例示的な実施形態ではデータ処理部２３０）は、第１マスク画像の画素値の範囲及び第２マスク画像の画素値の範囲に基づいて第１画像の画素値の範囲及び第２画像の画素値の範囲を調整するように制御される。

【0266】

また、プロセッサ（上記の例示的な実施形態ではデータ処理部２３０）は、第１画像に第１マスク画像を合成して第１合成画像を生成し、且つ、第２画像に第２マスク画像を合成して第２合成画像を生成するように制御される。

【0267】

また、プロセッサ（上記の例示的な実施形態ではデータ処理部２３０）は、第１合成画像及び第２合成画像に基づいて複数の相互相関関数を求めるように制御される。

【0268】

また、プロセッサ（上記の例示的な実施形態ではデータ処理部２３０）は、複数の相互相関関数に基づいて相関係数を算出するように制御される。

【0269】

また、プロセッサ（上記の例示的な実施形態ではデータ処理部２３０）は、相関係数に基づいて第１画像と第２画像とのレジストレーションを行うように制御される。

【0270】

また、プロセッサ（上記の例示的な実施形態ではデータ処理部２３０）は、レジストレーションの結果に基づいて第１画像と前記第２画像とのマージ画像を構築するように制御される。

【0271】

幾つかの例示的な態様において、走査型イメージング装置を制御する方法のいずれかの態様をコンピュータに実行させるプログラム、画像処理装置を制御する方法のいずれかの態様をコンピュータに実行させるプログラム、画像処理方法のいずれかの態様をコンピュータに実行させるプログラム、又は、走査型イメージング方法のいずれかの態様をコンピュータに実行させるプログラムを提供することが可能である。

【0272】

また、このようなプログラムを記録したコンピュータ可読な非一時的記録媒体を作成することが可能である。この非一時的記録媒体は任意の形態であってよく、その例として、磁気ディスク、光ディスク、光磁気ディスク、半導体メモリなどがある。

【0273】

以上に説明した幾つかの例示的な態様によれば、マスク画像を利用した画像レジストレーションにおいて、サンプルの画像の画素値の範囲とマスク画像の画素値の範囲とを相対的に調整した上で、レジストレーションのための処理を行うことが可能である。

【0274】

これにより、サンプルの画像の画素値の範囲とマスク画像の画素値の範囲との間の差を小さくすることができる。例えば、サンプルの画像の画素値の大きさと、マスク画像の画素値の大きさとを（ほぼ）一致させることが可能である。

【0275】

したがって、サンプルの画像とマスク画像とを用いた相関演算においてマスク画像の効果が無視されるという従来の問題を解決することができ、正しい相関係数を得ることが可能となる。

【0276】

特に、サンプルの画像同士の相関係数を求めるために単精度浮動小数点型（フロート型）の演算が用いられる場合であっても、有効桁数の少なさに起因する丸め誤差の影響を排除することや小さくすることが可能である。

【0277】

また、複数の相互相関関数から相関係数を算出するアルゴリズムにおいては、個々の相互相関関数に起因する誤差の重ね合わせにより、最終的な相関係数の誤差が大きくなる可能性がある。このような問題に関し、例示的な態様によれば、個々の相互相関関数に起因する誤差を排除することや小さくすることができるため、確度の高い相関係数を得ることが可能である。

【0278】

このように、以上に説明した幾つかの例示的な態様は、マスク画像を用いた画像レジストレーションの正確性の向上に寄与するものである。

【0279】

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

【符号の説明】

【0280】

１ 眼科装置
４４ 光スキャナ
１００ ＯＣＴユニット
２１１ 主制御部
２１１１ 走査制御部
２１２ 記憶部
２１２１ 走査プロトコル
２２０ 画像構築部
２３０ データ処理部
２３１ マスク画像生成部
２３２ 範囲調整部
２３３ 合成画像生成部
２３４ 相互相関関数算出部
２３５ 相関係数算出部
２３６ レジストレーション部
２３７ マージ処理部

【図1】

【図2】

【図3】

【図4A】

【図4B】

【図5】

【図6】

【図7A】

【図7B】

【図7C】