(57)【要約】位相ベースの機能ＯＣＴ像のバルク運動アーチファクトを補正する方法及びシステムが開示される。開示された方法及びシステムは、位相ベースのＯＣＴイメージングの際に存在する位相シフト信号の標準偏差を用いることに基づく。ＯＣＴ血管造影、ドップラーＯＣＴ、及びＯＣＴエラストグラフィ等の機能ＯＣＴ技法と共に適用されると、開示された方法は、バルク運動補償の他の方法と比較して改善された像品質及び低減された計算コストを実現する。【選択図】  図１Ａ
2021520244.tif 000002
【特許請求の範囲】
【請求項１】
位相ベースの機能光コヒーレンストモグラフィのバルク運動補償の方法であって、
第１の複素数値光コヒーレンストモグラフィのＡ走査Ａｎを受け取るステップと、
第２の複素数値光コヒーレンストモグラフィのＡ走査Ａｎ＋１を受け取るステップと、
それぞれ異なる位相における前記第１及び第２の複素数値光コヒーレンストモグラフィのＡ走査における分散値の組を決定するステップと、
前記分散値の組からバルク運動位相シフトΔφＢを決定するステップと、
前記バルク運動位相シフトΔφＢを用いて前記第１又は第２の複素数値光コヒーレンストモグラフィのＡ走査に対するバルク運動補償を行うステップと、
を含む方法。
【請求項２】
前記複素数値光コヒーレンストモグラフィのＡ走査における分散値の組を決定するステップが、ベース位相からθ＝０、π／２、π、及び３π／２におけるそれぞれの分散値を決定して、分散値ｖ１、ｖ２、ｖ３及びｖ４の組を生成することを含む、請求項１に記載の方法。
【請求項３】
前記バルク運動位相シフトΔφＢが、前記分散値ｖ１、ｖ２、ｖ３及びｖ４の組から
【数１】
2021520244.tif 000011
により決定されることをさらに含む、請求項２に記載の方法。
【請求項４】
前記分散値の組のうちの個々の分散値が、
ｖ＝［ｓｔｄ（Ａｎ−Ａｎ＋１・ｅｘｐ（−ｉ・θ））］２
により決定される、請求項２に記載の方法。
【請求項５】
前記分散値の組が分散値の第１の組であり、前記方法が、
前記第１の組を含む分散値の複数の組を、前記複数の組のそれぞれの組で異なるベース位相を基準として間隔があけられた、複数の異なる位相値において決定するステップと、
分散値の前記それぞれの組から決定された前記バルク運動位相シフトΔφＢを平均して、平均バルク運動位相シフトを得るステップと、をさらに含み、
前記バルク運動補償が、前記平均バルク運動位相シフトを用いて前記複素数値光コヒーレンストモグラフィのＡ走査に対し行われる、請求項１に記載の方法。
【請求項６】
前記バルク運動位相シフトΔφＢが、
ｖ＝Ｃ１−Ｃ２ｃｏｓ（ΔφＢ＋θ）
により複数の分散値の組から決定される、請求項１に記載の方法。
【請求項７】
バルク運動補償を行うステップが、第３の複素数値光コヒーレンストモグラフィのＡ走査Ａｎｅｗを式
Ａｎｅｗ＝Ａｎ＋１・ｅｘｐ（ｉ・ΔφＢ）
により計算することを含む、請求項１に記載の方法。
【請求項８】
光コヒーレンストモグラフィ血管造影のＡ走査ＯＣＴＡを
ＯＣＴＡ＝｜Ａｎ−Ａｎｅｗ｜
により補正するステップをさらに含む、請求項７に記載の方法。
【請求項９】
ドップラー光コヒーレンストモグラフィのＡ走査ｄｏｐｐｌｅｒを、
ｄｏｐｐｌｅｒ＝ａｎｇｌｅ（Ａｎ×ｃｏｎｊ（Ａｎｅｗ））
を計算することによって補正するステップをさらに含む、請求項７に記載の方法。
【請求項１０】
バルク運動補償された複素数値光コヒーレンストモグラフィのＡ走査に基づいて流れ像を生成するステップをさらに含む、請求項７に記載の方法。
【請求項１１】
複素数値を含む光コヒーレンストモグラフィ（ＯＣＴ）データセットを得るステップと、
前記光コヒーレンストモグラフィデータセットの複数のＡ走査の間の分散値の複数の組を決定するステップであって、前記複数の組の個々の組が、ベース位相に対して間隔があけられたそれぞれ異なる位相において決定された分散値を含み、前記ベース位相が、前記複数の組の個々の組ごとに異なる、ステップと、
分散値の個々の組ごとにバルク運動位相シフトを決定するステップと、
前記複数の組について決定された前記バルク運動位相シフトを平均してバルク運動位相補償を決定するステップと、
前記バルク運動位相補償に基づいて前記光コヒーレンストモグラフィデータセットを補償するステップと、
を含む方法。
【請求項１２】
前記それぞれ異なる位相が、前記ベース位相に対して互いに均等に間隔があけられている、請求項１１に記載の方法。
【請求項１３】
分散値の前記個々の組が、分散値ｖ１、ｖ２、ｖ３、及びｖ４を含み、分散値の個々の組ごとの前記バルク運動位相シフトΔφＢが、
【数２】
2021520244.tif 000012
により決定される、請求項１２に記載の方法。
【請求項１４】
前記分散値の組のうちの個々の分散値が
ｖ＝［ｓｔｄ（Ａｎ−Ａｎ＋１・ｅｘｐ（−ｉ・θ））］２
により決定され、ｖが前記個々の分散値であり、Ａｎが前記光コヒーレンストモグラフィデータセットの第１のＡ走査であり、Ａｎ＋１が前記光コヒーレンストモグラフィデータセットの第２のＡ走査であり、θが、前記個々の分散値が決定されたときの位相である、請求項１１に記載の方法。
【請求項１５】
前記光コヒーレンストモグラフィデータセットが光コヒーレンストモグラフィ血管造影データセットである、請求項１１に記載の方法。
【請求項１６】
前記光コヒーレンストモグラフィデータセットがドップラー光コヒーレンストモグラフィデータセットである、請求項１１に記載の方法。
【請求項１７】
位相ベースの機能光コヒーレンストモグラフィイメージングにおけるバルク運動補償のためのシステムであって、
試料の位相ベースの機能光コヒーレンストモグラフィデータセットを取得するように構成された光コヒーレンストモグラフィシステムと、
論理サブシステムと、
記憶した機械可読命令を含むデータ保持サブシステムであって、機械可読命令は、
第１の複素数値光コヒーレンストモグラフィのＡ走査Ａｎを受け取り、
第２の複素数値光コヒーレンストモグラフィのＡ走査Ａｎ＋１を受け取り、
前記第１及び第２の複素数値光コヒーレンストモグラフィのＡ走査における分散値の組を
ｖ＝［ｓｔｄ（Ａｎ−Ａｎ＋１・ｅｘｐ（−ｉ・θ））］２
により計算し、
前記分散値の組からバルク運動位相シフトΔφＢを計算し、
前記バルク運動位相シフトΔφＢを用いて前記複素数値光コヒーレンストモグラフィのＡ走査に対するバルク運動補償を行う、
ように前記論理サブシステムによって実行可能である、データ保持サブシステムと、
を備えるシステム。
【請求項１８】
前記分散値の組を計算するために、前記論理サブシステムが、ベース位相からθ＝０、π／２、π、及び３π／２におけるそれぞれの分散値を計算して、それぞれの分散値ｖ１、ｖ２、ｖ３及びｖ４を得る、請求項１７に記載のシステム。
【請求項１９】
前記バルク運動位相シフトΔφＢが、前記分散値の組から
【数３】
2021520244.tif 000013
により計算される、請求項１８に記載のシステム。
【請求項２０】
前記分散値の組が分散値の第１の組であり、前記命令が、
前記第１の組を含む分散値の複数の組を、前記複数の組のそれぞれの組で異なるベース位相を基準として間隔があけられた複数の異なる位相値において決定し、
分散値の前記それぞれの組から決定された前記バルク運動位相シフトΔφＢを平均して、平均バルク運動位相シフトを得る
ように前記論理サブシステムによってさらに実行可能であり、
前記バルク運動補償が、前記平均バルク運動位相シフトを用いて前記複素数値光コヒーレンストモグラフィのＡ走査に対し行われる、請求項１７に記載のシステム。
【請求項２１】
前記バルク運動補償を行うために、前記論理サブシステムが第３の複素数値光コヒーレンストモグラフィのＡ走査Ａｎｅｗを
Ａｎｅｗ＝Ａｎ＋１・ｅｘｐ（ｉ・ΔφＢ）
により計算する、請求項１７に記載のシステム。
【請求項２２】
前記命令が、光コヒーレンストモグラフィ血管造影のＡ走査ＯＣＴＡを
ＯＣＴＡ＝｜Ａｎ−Ａｎｅｗ｜
により補正するように、前記論理サブシステムによってさらに実行可能である、請求項２１に記載のシステム。
【請求項２３】
前記命令が、ドップラー光コヒーレンストモグラフィのＡ走査ｄｏｐｐｌｅｒを
ｄｏｐｐｌｅｒ＝ａｎｇｌｅ（Ａｎ×ｃｏｎｊ（Ａｎｅｗ））
により補正するように、前記論理サブシステムによってさらに実行可能である、請求項２１に記載のシステム。
【請求項２４】
さらに、前記命令が、前記バルク運動補償された複素数値光コヒーレンストモグラフィのＡ走査に基づいて流れ像を生成し表示するように、前記論理サブシステムによってさらに実行可能である、請求項１７に記載のシステム。
【発明の詳細な説明】
【関連出願の相互参照】
【０００１】
本願は、２０１８年４月２日に出願された「ＳＹＳＴＥＭＳ  ＡＮＤ  ＭＥＴＨＯＤＳ  ＦＯＲ  ＢＵＬＫ  ＭＯＴＩＯＮ  ＣＯＭＰＥＮＳＡＴＩＯＮ  ＩＮ  ＰＨＡＳＥ−ＢＡＳＥＤ  ＦＵＮＣＴＩＯＮＡＬ  ＯＰＴＩＣＡＬ  ＣＯＨＥＲＥＮＣＥ  ＴＯＭＯＧＲＡＰＨＹ」という名称の米国特許仮出願第６２／６５１，６７４号の優先権を主張し、その開示全体が参照により本明細書に組み込まれる。
【技術分野】
【０００２】
一般に、本技術分野は、光コヒーレンストモグラフィを用いるイメージングの方法を含む。特に、本技術分野は、位相ベースの機能光コヒーレンストモグラフィにおけるバルク運動補償の方法を含む。
【政府支援の確認】
【０００３】
本発明は、米国国立衛生研究所から認められたＲ０１  ＥＹ０２７８３３、Ｒ０１  ＥＹ０２４５４４、及びＤＰ３  ＤＫ１０４３９７の政府支援により行われた。政府は、本発明の一定の権利を有する。
【背景技術】
【０００４】
光コヒーレンストモグラフィ（ＯＣＴ）は、生体組織をマイクロメートル分解能で撮像するために用いられる三次元（３Ｄ）イメージング技法である。フーリエ領域ＯＣＴの出現と共に、この技術の速度及び感度の向上が、ドップラーＯＣＴ及びＯＣＴ血管造影法（ＯＣＴＡ）による血流イメージング等の新しい機能を可能にしてきた。ドップラーＯＣＴでは、大血管系内の軸線方向に沿った流速を得るために、隣り合うＡ走査間の位相シフトを検出する。より最近では、ＯＣＴＡが、同じ位置で取得された隣り合うＢ走査のＯＣＴ信号を差し引くことによって微小血管系内の内因性血流信号を抽出するために開発された。振幅信号又は位相信号に基づくＯＣＴＡアルゴリズムが市販システムに成功裏に取り入れられており、眼イメージングのために臨床医によって日常的に使用されている。
【０００５】
しかし、位相ベースの機能ＯＣＴ（例えば、位相分解ドップラーＯＣＴ、位相ベースＯＣＴＡ）の１つの主要な課題として、バルク運動によって引き起こされる（例えば、呼吸、心拍及び他の不随意運動によって誘発される）付加的な位相シフト（差）がある。位相ベースの方法は、振幅ベースの方法よりもバルク運動アーチファクトの影響を受けやすいので、流れ信号を正確に検出するにはバルク運動位相シフトを除去する手段が不可欠である。
【０００６】
２つの位相補償法が以前に報告されており、広く使用されている。１つは「平均」法であり、網膜組織スラブ中のＢフレームの各Ａラインに沿って位相シフトを平均することによって、バルク運動位相シフトを推定する。この方法は、血管及びその影が組織の比較的大きい部分を占める場合には大きな誤差を導入するおそれがある。もう１つは、「ヒストグラム」法であり、各Ａラインの組織を横断してすべての値のヒストグラムでビン数（ピーク密度）が最も大きい位相シフトを見つけることによって、このアーチファクトを推定する。「ヒストグラム」法は、推定されるバルク運動位相シフトを各Ａラインで繰り返し差し引くことによって「平均」法よりも正確になり得る。しかし、この方法は大きな計算バジェットを必要とする。したがって、位相ベースの機能ＯＣＴにおけるバルク運動補償のための高速で堅牢な方法が依然として必要とされている。
【図面の簡単な説明】
【０００７】
【図１Ａ】様々な実施形態による、位相ベースのＯＣＴＡにおけるバルク運動位相シフトを補償するための例示的な手順を示すフローチャートである。
【図１Ｂ】様々な実施形態による、ベース位相θｍにおけるバルク運動位相シフトΔφＢｍの推定を示すフローチャートである。
【図２Ａ】ラットの眼から生体内で得られたＯＣＴイメージングデータのフィッティングにおける、開示された方法の適性を示す図である。図２Ａは、異なる箇所におけるｖ対θのグラフである。図２Ａのデータは、毛細管のみを含む場所におけるＡラインからのものである。
【図２Ｂ】ラットの眼から生体内で得られたＯＣＴイメージングデータのフィッティングにおける、開示された方法の適性を示す図である。図２Ｂは、異なる箇所におけるｖ対θのグラフである。図２Ｂのデータは、大血管を含む場所からのものである。
【図２Ｃ】ラットの眼から生体内で得られたＯＣＴイメージングデータのフィッティングにおける、開示された方法の適性を示す図である。図２Ｃは、異なる箇所におけるｖ対θのグラフである。図２Ｃのデータは、大血管を含む場所からのものである。
【図２Ｄ】ラットの眼から生体内で得られたＯＣＴイメージングデータのフィッティングにおける、開示された方法の適性を示す図である。図２Ｄは、ラットの眼の走査の正面ＯＣＴ平均投影像を示す。図２ＤにＡ、Ｂ、及びＣと標示された位置は、それぞれ図２Ａ〜Ｃに対応する。
【図３Ａ】ラットの眼の断面データを描写する図、及び開示された方法によってもたらされる像改善を示す図である。図３Ａは、バルク運動補償前の位相差像を描写する。
【図３Ｂ】ラットの眼の断面データを描写する図、及び開示された方法によってもたらされる像改善を示す図である。図３Ｂは、「標準偏差」（ＳＴＤ）法によるバルク運動補償後の位相差像を描写する。
【図３Ｃ】ラットの眼の断面データを描写する図、及び開示された方法によってもたらされる像改善を示す図である。図３Ｃは、バルク運動補償のない位相ベースのＯＣＴＡ像を示す。
【図３Ｄ】ラットの眼の断面データを描写する図、及び開示された方法によってもたらされる像改善を示す図である。図３Ｄは、ＳＴＤによるバルク運動補償後の位相ベースのＯＣＴＡ像を描写する。
【図３Ｅ】ラットの眼の断面データを描写する図、及び開示された方法によってもたらされる像改善を示す図である。図３Ｅは、バルク運動補償前のドップラーＯＣＴ像を描写する。
【図３Ｆ】ラットの眼の断面データを描写する図、及び開示された方法によってもたらされる像改善を示す図である。図３Ｆは、ＳＴＤ法によるバルク運動補償後のドップラーＯＣＴ像を描写する。
【図４Ａ】ラットの眼に適用された異なるバルク運動補償方法を比較して示す図である。図４Ａは、バルク運動補償をせずに計算された流れ信号を示す（ＮＯＮＥと標示）。白い矢印は、いくつかの場所における運動アーチファクトがＡＶＧ法及びＳＴＤ法によってどこで抑制され得るかを示し（図４Ｂ及び図４Ｄ）、またＨＩＳＴ法によってどこが誤って強調され得るかも示す（図４Ｃ）。
【図４Ｂ】ラットの眼に適用された異なるバルク運動補償方法を比較して示す図である。図４Ｂは、「平均」（ＡＶＧ）法を使用して計算された流れ信号を描写する。白い矢印は、いくつかの場所における運動アーチファクトがＡＶＧ法及びＳＴＤ法によってどこで抑制され得るかを示し（図４Ｂ及び図４Ｄ）、またＨＩＳＴ法によってどこが誤って強調され得るかも示す（図４Ｃ）。
【図４Ｃ】ラットの眼に適用された異なるバルク運動補償方法を比較して示す図である。図４Ｃは、「ヒストグラム」（ＨＩＳＴ）法を使用して補償された流れ信号を描写する。白い矢印は、いくつかの場所における運動アーチファクトがＡＶＧ法及びＳＴＤ法によってどこで抑制され得るかを示し（図４Ｂ及び図４Ｄ）、またＨＩＳＴ法によってどこが誤って強調され得るかも示す（図４Ｃ）。
【図４Ｄ】ラットの眼に適用された異なるバルク運動補償方法を比較して示す図である。図４Ｄは、「標準偏差」（ＳＴＤ）法を使用して補償された流れ信号を描写する。白い矢印は、いくつかの場所における運動アーチファクトがＡＶＧ法及びＳＴＤ法によってどこで抑制され得るかを示し（図４Ｂ及び図４Ｄ）、またＨＩＳＴ法によってどこが誤って強調され得るかも示す（図４Ｃ）。
【図５Ａ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図５Ｂ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図５Ｃ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図５Ｄ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図５Ｅ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図５Ｆ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図５Ｇ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図５Ｈ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図５Ｉ】人間の網膜内の異なる深さで適用された異なるバルク運動補償方法を比較して示す複数の像のうちの１つである。
【図６】ＡＶＧ法、ＨＩＳＴ法、及びＳＴＤ法に必要とされる計算バジェットを比較して示す棒グラフである。
【図７】本開示による、ＯＣＴデータセット及びＯＣＴ血管造影データセットを処理して非灌流の領域を検出する例示的なシステムを概略的に示す図である。
【図８】本開示による、コンピュータシステムの一例を概略的に示す図である。
【図９】様々な実施形態による、ＯＣＴのバルク運動位相補正を計算する方法のフローチャートである。
【発明を実施するための形態】
【０００８】
ＯＣＴ適用においてバルク運動アーチファクトを除去する方法及びシステムが開示される。諸実施形態では、開示された方法及びシステムは、位相ベースのＯＣＴイメージングの際に存在する位相シフト信号の標準偏差を用いることに基づく。
【０００９】
開示されたシステム及び方法は、振幅と位相の両方の情報を有するＯＣＴ信号等の複素数値信号に適用可能である。いくつかの実施形態では、開示された方法及びシステムは、位相ベースの機能ＯＣＴイメージング用途に適用可能である。機能ＯＣＴイメージングの例には、網膜ドップラーＯＣＴ、ＯＣＴ血管造影（ＯＣＴＡ）、及びＯＣＴエラストグラフィが含まれる。
【００１０】
開示された方法及びシステムの一態様は、これら方法及びシステムが、位相ベースの機能ＯＣＴイメージングにおけるバルク運動補償のための堅牢で計算効率の良い手法を実現することである。当技術分野で知られているバルク運動補償方法と比較して、開示された方法は、像品質の改善を処理時間の低減と共に実現する。
【００１１】
本明細書では、開示された方法を用いてバルク運動補償を取得及び実行する例示的なシステムがさらに開示される。この例示的なシステムは、論理サブシステム及びデータ保持機能を有するコンピュータデバイスと機能的に接続してＯＣＴ構造データ及び血管造影データを取得するように構成されたＯＣＴデバイスを含む。諸実施形態では、コンピュータデバイスは、ＯＣＴデバイスからデータを受け取り、本明細書に記載の方法の１つ又は複数の動作を実行するように構成される。
【００１２】
以下の詳細な説明では、本明細書の一部を形成する添付の図面が参照され、実践できる実施形態が例示的に示される。他の実施形態が利用されることもあり、構造的又は論理的変更がその範囲から逸脱することなく加えられ得ることを理解されたい。したがって、以下の詳細な説明は、限定的な意味に解釈されるべきものではなく、諸実施形態の範囲は、添付の特許請求の範囲及びその等価物によって定義される。
【００１３】
様々な動作は、複数の個別の動作として順に、諸実施形態を理解する助けになり得るように説明されることがあるが、記述の順序は、これらの動作が順序に依存することを示唆するものと解釈されるべきではない。
【００１４】
説明では「実施形態」又は「諸実施形態」という用語を使用することがあり、これらはそれぞれ、同一又は異なる諸実施形態のうちの１つ又は複数を指し得る。さらに、諸実施形態に関して使用される用語の「備えている」、「含んでいる」、「有している」等の用語は同義である。
【００１５】
様々な実施形態では、試料の構造及び／又は流れの情報は、スペクトル干渉を検出することに基づくＯＣＴ（構造）及びＯＣＴ血管造影（流れ）イメージングを用いて得ることができる。このようなイメージングは、用途に応じて二次元（２Ｄ）又は三次元（３Ｄ）とすることができる。構造イメージングは、従来技術の方法に対して拡張された深さ範囲のものとすることができ、流れイメージングは、リアルタイムで実行することができる。２Ｄ像又は３Ｄ像を生成するために、本明細書に開示の構造イメージングと流れイメージングの一方又は両方を確保することができる。
【００１６】
特にことわらない、又は説明しない限り、本明細書で用いられるすべての技術的及び科学的用語は、従来の用法に従って使用され、当業者によって一般に理解されるものと同じ意味を有する。本明細書に記載のものと同様又は同等の方法、システム、及び装置／材料は、本開示の実践及び試験に使用することができるが、適切な方法、システム、及び装置／材料を以下で説明する。
【００１７】
本明細書で言及されるすべての出版物、特許出願、特許、及び他の参照文献は、参照によりその全体が組み込まれる。矛盾する場合には、本明細書が用語の説明を含めて支配する。加えて、方法、システム、装置、材料、及び例は、例示的なものにすぎず、限定的なものではない。
【００１８】
本開示の様々な実施形態を評価しやすくするために、特定の用語についての以下の説明を提示する。
【００１９】
Ａ走査：関心アイテム内の構造体の空間寸法及び位置についての情報を含む反射率特性。Ａ走査は、ＯＣＴデバイスの光軸線に沿って行われる軸線方向走査であり、撮像される試料を透過する。Ａ走査は、反射率情報（例えば、信号強度）を深さ（ｚ方向）の関数として符号化する。
【００２０】
Ｂ走査：一連の軸線方向深さ走査（すなわち、Ａ走査）をｘ方向又はｙ方向で横に結合することによって実現できる断面トモグラフィ。Ｂ走査は、試料からの平面断面情報を符号化し、通常は１つの像として提示される。したがって、Ｂ走査を断面像と呼ぶことができる。
【００２１】
データセット：本明細書で用いられる、データセットは記憶データ値の順序配列表示であり、行−列−深さ（ｘ−ｙ−ｚ軸線）形式で相対的空間位置を符号化する。ＯＣＴとの関連において、本明細書では、１つのデータセットをボクセルの３次元配列として概念化することができ、各ボクセルは関連する値を持つ（例えば、強度値、振幅と位相の両方の情報を有する複素数値、非相関値、又は他の信号表示）。１つのＡ走査が、データセットの深さ（ｚ軸線）方向に沿った同一線上の１組のボクセルに対応し、１つのＢ走査が、行又は列（ｘ又はｙ軸線）方向に結合された隣り合うＡ走査のセットから構成される。このようなＢ走査は像と呼ばれることもあり、その構成ボクセルがピクセルと呼ばれることもある。隣り合うＢ走査を集めたものが結合されて、３Ｄ像と呼ばれるボクセルデータの３Ｄ体積セットが形成され得る。本明細書に記載のシステム及び方法では、ＯＣＴ走査デバイスによって得られたデータセットは「構造ＯＣＴ」データセットと名付けられ、その値は、例えば強度及び位相情報を保持する複素数になり得る。この構造ＯＣＴデータセットは、撮像された試料内の流れを反映する、「ＯＣＴ血管造影」データセットと名付けられた対応するデータセットを計算するために使用することができる。構造ＯＣＴデータセットのボクセルとＯＣＴ血管造影データセットの間には対応関係がある。したがって、データセットからの値は「重ね合わせ」て、構造と流れ（例えば、組織微小構造体と血流）の複合像、又は別の方法で結合若しくは比較された複合像を提示することができる。
【００２２】
正面血管造影図：ＯＣＴ血管造影データは、正面血管造影図と呼ばれる単一の平面像への三次元データセットの２Ｄ投影として提示することができる。このような正面血管造影図の構築には、血管造影図に投影されるべき網膜ＯＣＴ走査内の関心領域を取り囲む上下の深さ範囲を指定することが必要になる。これらの上下の深さ範囲は、網膜の異なる層の間の範囲として指定することができる（例えば、内境界膜と外網状層の間のボクセルは、網膜内層の正面血管造影図を生成するために用いることができる）。生成されると正面血管造影図は、本明細書に記載のように網膜血管系の様々な特徴を定量化するために使用することができる。この定量化は通常、例えば、静的組織からの血管系内の流れを血管造影図中で表すピクセルを区分するための閾値の設定を伴う。これらの正面血管造影図は、蛍光眼底造影（ＦＡ）又はインドシアニングリーン（ＩＣＧ）眼底造影等の従来の血管造影技法と同様に解釈することができ、したがって臨床用途に適している。正面血管造影図を生成するのに用いられるものと同様の方法で、構造ＯＣＴデータから正面像を生成することも一般的である。異なる層からの血管造影図は、色符号化し重ね合わせて符号化深さ情報を含む複合血管造影図を提示することもでき、構造正面像もこのような複合像生成に含めることができる。
【００２３】
本明細書で用いられる、機能ＯＣＴは、構造特性以上の情報を提供するためのＯＣＴ技法の拡張のことを広く指す。例えば、構造ＯＣＴイメージングが組織の解剖学的構成物についての空間情報を集めるために使用できるのに対し、機能ＯＣＴは、その組織試料内で起こる血流、組織灌流及び酸素化、複屈折等の過程についての情報を集めるために使用することができる。機能ＯＣＴの例には、血流の特性を示すためのＯＣＴ血管造影（ＯＣＴＡ）及び関連する技法、ドップラーＯＣＴ、偏光感受型ＯＣＴ、ＯＣＴエラストグラフィ、分光ＯＣＴ、差分吸収ＯＣＴ、及び分子イメージングＯＣＴが、これらだけには限らないが含まれる。
【００２４】
本明細書では、位相ベースの光コヒーレンストモグラフィ（ＯＣＴ）のバルク運動アーチファクトを補償するための方法及びシステムが開示される。ＯＣＴ血管造影（ＯＣＴＡ）、ドップラーＯＣＴ、及びＯＣＴエラストグラフィ等の位相ベースのＯＣＴ適用は、対象物バルク運動によって導入される交絡位相シフトの影響を受けやすい。従来のバルク運動補償方法には、その精度及び計算費用対効果に限度がある。本明細書に記載の、位相ベースの機能ＯＣＴのための新規のバルク運動補償方法には、像品質の改善と計算負担の低減の両面で従来の方法に対する利点がある。
【００２５】
本明細書では、ＯＣＴ流れ信号の標準偏差に基づくバルク運動補償の新規の方法が開示される。バルク運動の陽型式を導出し、位相ベースの機能ＯＣＴから差し引くことができる。開示されたバルク運動補償方法の一態様は、ＯＣＴデータセット中の運動アーチファクトを低減することに関し、以前の平均化及びヒストグラムベースの方法よりも効果的なことである。加えて、この方法は、これらの以前の方法よりも計算的にあまり集約的ではない。
【００２６】
１つの特定の横方向位置、軸線方向位置ｚ、及び時間ｔにおけるＯＣＴ複素数値信号（振幅と位相）は、次式のように記述することができる。
Ｂｏｃｔ（ｔ，ｚ）＝Ａ（ｔ，ｚ）ｅｘｐ［ｉφ（ｔ，ｚ）］    式（１）
ただし、Ｂｏｃｔ（ｔ，ｚ）は複素数値信号であり、Ａ（ｔ，ｚ）は振幅であり、φは位相である。位相項は３つの成分を含み、次式のように展開することができる。
φ（ｔ，ｚ）＝φＰ（ｔ，ｚ）＋φＢ（ｔ，ｚ）＋φＭ（ｔ，ｚ）    式（２）
ただし、φＰ（ｔ，ｚ）は血管内部の移動する赤血球によって生じる位相シフトであり、φＢ（ｔ，ｚ）はバルク運動によって生じる位相シフトであり、φＭ（ｔ，ｚ）はシステムの変調周波数によって生じる位相シフトである。２つの繰り返しＢ走査間の差の絶対値は次式のように表すことができる。
【数１】
2021520244.tif 000003
ここで、第２の走査の時間座標はｔ２＝ｔ１＋Δｔと表すことができる。システム変調位相差は、全Ａ走査について一定と仮定され、除去することができる。ΔφＢの項は、Ｂｄｉｆｆの標準偏差分析として表すことができる。
【数２】
2021520244.tif 000004
式（３）を式（４）に代入すると次式が得られる。
【数３】
2021520244.tif 000005
ただし、Ｃ１及びＣ２は定数である。血管の位相シフト寄与は組織のものと比べて相対的に小さいので、この計算ではΔφＰを無視することができる。したがって、分散ｖは次式のように表される。
ｖ＝ｓ２＝Ｃ１−Ｃ２ｃｏｓ（ΔφＢ）    式（６）
【００２７】
陽型バルク運動式に基づいて、バルク運動位相を直接解くことができる。この式を解くために、複素数値信号Ｂｏｃｔ（ｔ＋Δｔ，ｚ）に対する位相θの疑似減算が実行される。
Ｂ’ｏｃｔ（ｔ＋Δｔ，ｚ）＝Ｂｏｃｔ（ｔ＋Δｔ，ｚ）・ｅｘｐ（−ｉθ）    式（７）
【００２８】
この場合、分散を次式のように表すことができる。
ｖ＝Ｃ１−Ｃ２ｃｏｓ（ΔφＢ＋θ）    式（８）
【００２９】
この分散は、位相θを有するコサイン関数である。バルク運動位相ΔφＢは、π／２の間隔があけられた４つの位相値θにおける分散値で直接解くことができる。例えば、諸実施形態において、θ＝０、π／２、π、３π／２の位相値が次の計算を始めるために選択されてよい。
【数４】
2021520244.tif 000006
【００３０】
ここで、Ａｎは、同じ場所で連続するｎ番目のＢフレームのｎ番目の繰り返しＡラインである。この場合、バルク運動位相シフトは次式のように書くことができる。
【数５】
2021520244.tif 000007
ここで、ａｔａｎ２は、四象限逆タンジェントである。したがって、式（１４）に基づいて、直接にバルク運動位相シフトを解くことができる。
【００３１】
諸実施形態で、式（１４）の計算精度は、４つのサンプリング点の値をシフトして式を複数回解くことによって向上させることができる。例えば、以下に提示した例における計算の計算精度を向上させるために、式（１４）は、４つのサンプリング点をランダム位相ベースθｍ（すなわちθｍ、θｍ＋π／２、θｍ＋π及びθｍ＋３π／２）と共にシフトすることによって５回解くことができる。５つのΔφＢを計算し平均して、Ａラインごとに最終バルク運動位相シフトΔφＢを得ることができる。別の実施形態では、一緒に平均される任意の適切な数のΔφＢｍ値を用いることができる。
【００３２】
本明細書に記載の様々な実施形態による、バルク運動補償の方法１００を説明するフローチャートが図１Ａ及び図１Ｂに示されている。諸実施形態では、方法１００はＯＣＴＡ走査に適用することができる。いくつかの実施形態では、方法１００は、図１Ａの１０２で示された事前計算操作を含み得る。この事前計算操作１０２は、例えば、生ＯＣＴスペクトルデータ１０６から複素数ＯＣＴ信号データを抽出するために、１つ又は複数の信号処理操作を含み得る（例えば、方法１００の１１０に示されるように、生ＯＣＴスペクトルデータからベーススペクトル１０８が差し引かれる）。このような信号処理操作には、分散補償、波長表示から周波数表示への変換、波数補間、及び／又は高速フーリエ変換操作のための１つ又は複数の操作が、これらだけには限らないが含まれ得る。
【００３３】
諸実施形態では、次にバルク運動が、本明細書に記載の「標準偏差」法を用いて複素数ＯＣＴ信号データ１０４から推定され得る。例えば、上述のように、１つ又は複数のベース位相（例えば、θ１、θ２、．．．θｍ）が方法１００の１１２ａ〜ｃで選択されてよい。１つ又は複数のベース位相は、ランダムに選択されても半ばランダムに選択されてもよく、或いは所定の位相値であってもよい。バルク運動による位相差ΔφＢｍは、方法１００の１１４ａ〜ｃでベース位相ごとに決定することができる。
【００３４】
図１Ｂは、単一のベース位相に対する位相差ΔφＢｍを決定する方法１５０をより詳細に示す。例えば、方法１５０は、ベース位相に対するπ／２間隔の４つの位相値における分散値を式（１０）〜（１３）に従って決定することを含む（図１Ｂに操作１５２ａ〜ｄとして示す）。次に、所与のベース位相に対するバルク運動位相差ΔφＢｍが、方法１５０の１５４で決定される（例えば、式（１４）に従って）。
【００３５】
図１Ａを再び参照すると、複数のベース位相に対するバルク運動位相差ΔφＢｍが方法１００の１１６で平均されて、Ａラインごとのバルク運動位相差ΔφＢが得られる。１１６での平均化操作は、１つのベース位相における分散値の単一の組が用いられる場合には省くことができる。方法１００の１１８で、ＯＣＴデータセットの個々のＡラインがバルク運動位相シフトΔφＢに基づいて調整されて、バルク運動補償ＯＣＴデータセットが得られる。
【実施例】
【００３６】
以下の実施例では、開示された方法を説明している。本開示に照らし合わせて、当業者には、本開示の方法のこれらの実施例の変形形態及び他の実施例が、必要以上の実験がなくても可能であることが理解されよう。
【００３７】
実施例１  「標準偏差」法の精度
提案された方法の精度は、一部には、実験結果が式（８）で予測されたコサイン曲線表示とどれほど相関関係にあるかに依拠する。言い換えると、結果がよく合うには、コサイン関数の最初の最小値において正確なバルク運動位相シフトを取り込む必要がある。このことが事実であることを実験データで検証するために、単一のＡ走査が、生体内ラット眼のＯＣＴデータセット中の３つの場所（毛細血管だけがある静的組織内の１つの場所、及び大血管内の２つの場所）からランダムに選択された。データからの様々なθにおいて計算された分散と、式（８）から直接導出された分散との間の関係が調査された。
【００３８】
図２Ａ〜２Ｃは、ラットの眼の中の３つの異なる点（図２Ｄに領域Ａ、Ｂ、及びＣと示されている）におけるｖ対θのグラフを示す。グラフ上の「＋」マークは、−１〜１に正規化された、計算された分散である。曲線２０２ａ〜ｃは、式８から導出されたコサイン関数であり、各グラフの左側の領域２０４ａ〜ｃは、バルク運動によって生じた位相シフトを示す。図２Ａで、Ａラインは毛細血管のみを含む箇所から選択されているが、図２Ｂ及び図２Ｃにあるのは、大血管を含む箇所から選択されたＡラインである。図２Ｄは、ラットの眼の走査の正面ＯＣＴ平均投影像を示す。図２Ｄ中にＡ、Ｂ、及びＣと標示された位置は、それぞれ図２Ａ、２Ｂ及び２Ｃのグラフに対応する。図２Ａ、２Ｂ及び２Ｃにおける高品質の曲線適合は、式（８）が適正であること、及びバルク運動位相シフトは式（１４）を用いて直接解けることを示唆する。
【００３９】
実施例２  齧歯類の眼における適用
バルク運動補償のための開示されたシステム及び方法を実証するために、網膜ドップラーＯＣＴ及びＯＣＴＡ走査がラットの眼で取得された。バルク運動補償システム及び方法の、像品質及び計算ニーズを含む性能が評価された。例示的なプロトタイプ可視光スペクトル領域ＯＣＴシステムが、齧歯類を撮像するために使用された。システムは、走査速度が５０ｋＨｚ、中心波長が５６０ｎｍ、及び軸線方向分解能が１．５３μｍであった。実験手順はすべて動物実験委員会（ＩＡＣＵＣ：Ｉｎｓｔｉｔｕｔｉｏｎａｌ  Ａｎｉｍａｌ  Ｃａｒｅ  ａｎｄ  Ｕｓｅ  Ｃｏｍｍｉｔｔｅｅ）によって承認された。
【００４０】
図３Ａ〜３Ｆは、ＯＣＴＡ適用とドップラーＯＣＴ適用の両方に対する本開示の方法の有効性を実証する生体内ラット撮像実験の代表的なデータを示す。断面ＯＣＴ像は、図２Ｄに示したものと同じ体積走査により取得された。ＯＣＴＡ像が、図２Ｄに水平ライン２０６で示した位置で取得された。ドップラーＯＣＴ像が、図２Ｄに円２０８で示した円走査により取得された。
【００４１】
図３Ａは、バルク運動補償前の位相差像を示し、色付き縞模様はバルク運動位相シフトを示す。図３Ｂは、「標準偏差」（ＳＴＤ）法によるバルク運動補償後の位相差像を示す。図３Ｃは、バルク運動補償のない位相ベースのＯＣＴＡ像を示し、網膜流れ信号は運動アーチファクトから区別することができない。図３Ｄは、開示されたＳＴＤ法を用いたバルク運動補償後の位相ベースのＯＣＴＡ像を示し、真の流れ信号（黄矢印で示されている）は区別できるようになっている。図３Ｅは、バルク運動補償前のドップラーＯＣＴ像を示し、多量の運動信号及び低ドップラー信号コントラストが存在する。図３Ｆは、ＳＴＤ法によるバルク運動補償後のドップラーＯＣＴ像を示し、ドップラー信号コントラストが改善されている。図３Ａ、３Ｂ、３Ｅ、及び３Ｆのカラーバーの範囲は、−π〜πの位相である。図３Ｃ及び図３Ｄのカラーバーの範囲は、０〜０．２の正規化振幅である。
【００４２】
開示された「標準偏差」法はさらに、ラットの眼からの同じデータセットについて２つの従来の方法（「平均」法及び「ヒストグラム」法）と比較することによって、正面血管造影図について定性的及び定量的に評価された。３つのバルク運動除去法の間の正当な比較を確実にするために、用いられた位相補償技法を除いて、全データが同じ操作を用いて処理された。
【００４３】
補償ステップのないＯＣＴＡ（図４Ａ）と比較して、「平均法」によって補償された像（図４Ｂ）は、いくつかの微小血管の細部を明らかにしているが、高い背景ノイズ（バルク運動アーチファクト）を保持している。「ヒストグラム」法（図４Ｃ）では、アーチファクトが抑制され、像コントラストが改善されたが、補償ステップによって導入された誤差があった（白矢印で示されている）。全体として、「標準偏差」法によって処理された像（図４Ｄ）は、齧歯類網膜イメージングに関して最良の像品質を示した。
実施例３  人間の眼における適用
【００４４】
バルク運動補償のための開示されたシステム及び方法を実証するために、網膜ドップラーＯＣＴ及びＯＣＴＡ走査が人間の眼で取得された。バルク運動補償システム及び方法の、像品質及び計算ニーズを含む性能が評価された。走査速度が７０ｋＨｚ、中心波長が８４０ｎｍ、及び軸線方向分解能が５μｍである市販のスペクトル領域ＯＣＴシステム（ＲＴＶｕｅ−ＸＲ  Ａｖａｎｔｉ、Ｏｐｔｏｖｕｅ，Ｉｎｃ．、Ｆｒｅｍｏｎｔ、ＣＡ）が、人間の網膜を撮像するために使用された。書面のインフォームドコンセントがすべてのボランティアから得られた。すべての実験手順がＯＨＳＵの倫理審査委員会（Ｉｎｓｔｉｔｕｔｉｏｎａｌ  Ｒｅｖｉｅｗ  Ｂｏａｒｄ／Ｅｔｈｉｃｓ  Ｃｏｍｍｉｔｔｅｅ）によって承認された。黄斑の３ｍｍ×３ｍｍの領域を包含する、人間の１０網膜走査が５人のボランティアから取得された。
【００４５】
開示された「標準偏差」バルク運動除去法が、従来の「平均」法及び「ヒストグラム」法に対して、それぞれの方法を用いて同じデータセットを処理することによって比較された。正面血管造影図を用いて、定性的評価と定量的評価の両方が行われた。３つのバルク運動除去法の間の正当な比較を確実にするために、用いられた位相補償技法を除いて、全データがワークフローの同じ信号処理操作を用いて処理された。表面血管複合体（ＳＶＣ）、中間毛細血管叢（ＩＣＰ）、及び深部毛細血管叢（ＤＣＰ）の正面像が、カスタム自動区分ソフトウェアを用いて生成された。
【００４６】
図５Ａ〜５Ｌは、３つの異なるスラブ（ＳＶＣ、ＩＣＰ、及びＤＣＰ）における流れの正面最大投影についての異なるバルク運動補償方法の効果を示す。図５Ａ〜５Ｌの３つの列は、分析された３つの異なるスラブに対応する。第１の列（すなわち、図５Ａ、５Ｄ、５Ｇ、及び５Ｊ）は、表面血管複合体（ＳＶＣ）を示し、第２の列（すなわち、図５Ｂ、５Ｅ、５Ｈ、及び５Ｋ）は、中間毛細血管叢（ＩＣＰ）を示し、第３の列（すなわち、図５Ｃ、５Ｆ、５Ｉ、及び５Ｌ）は、深部毛細血管叢（ＤＣＰ）を示す。図５Ａ〜５Ｌの４つの行は、適用された４つの異なるバルク運動補償方法に対応する。第１の行（図５Ａ〜５Ｃ）は、バルク運動補償なしで計算され（ＮＯＮＥ）、第２の行（図５Ｄ〜５Ｆ）は、「平均」（ＡＶＧ）法を用いて計算され、第３の行（図５Ｇ〜５Ｉ）は、「ヒストグラム」（ＨＩＳＴ）法を用いて計算され、第４の行（図５Ｊ〜５Ｌ）は、「標準偏差」（ＳＴＤ）法を用いて計算された。比較された方法を用いて取り込まれた細部のこれらの像品質及びレベルは、「平均」法及び「ヒストグラム」法と比較して優れた像品質を生み出す「標準偏差」法の能力を立証する。
【００４７】
バルク運動補償後の像品質はさらに、流れ信号対ノイズ比（ＳＮＲ）を計算することによって調査された。中心窩無血管域（ＦＡＺ）は、流れがない知られている目印であるので、健康な網膜の好適なノイズ基準になり、ＳＮＲは、傍中心窩の流れとＦＡＺ信号の比によって計算することができる。神経節細胞層叢は、神経線維内の多量の残留運動信号の干渉を除去するために、ＳＮＲを分析するのに使用された唯一のスラブであった。ＳＮＲ計算により、「標準偏差」法では６２％の像品質改善が「平均」法と比較して得られ、「ヒストグラム」法では２２％の像品質改善が「平均」法に対して得られたことが示された（表１）。
【表１】
2021520244.tif 000008
【００４８】
実施例４  計算性能
図６は、３つのバルク運動除去計算、すなわち「平均」法、「ヒストグラム」法、及び開示された「標準偏差」法の計算性能（又は計算バジェット）を示すグラフである。この比較に使用されたこれらの計算プラットフォームは、Ｉｎｔｅｌ  Ｃｏｒｅ  ｉ７  ６７００ｋ、３２ＧＢ  ＲＡＭの単一処理システムであり、各方法がＭａｔｌａｂ  ２０１７ａ（Ｔｈｅ  Ｍａｔｈｗｏｒｋｓ、Ｎａｔｉｃｋ、ＭＡ）を使用して実施された。図６に示されるように、「平均」法は計算バジェットが最も低いが、像品質も最低である。「ヒストグラム」法は、「平均」法と比較して像品質が高いが、必要な繰り返し数のために計算バジェットが最大であった。「標準偏差」法は、計算バジェットが比較的少ない上に像品質が最良である。バルク運動計算で「標準偏差」法は、「ヒストグラム」法よりほぼ５倍も速かった。
【００４９】
実施例５  計算手順
図９は、ＯＣＴＡ適用及びドップラー適用のバルク運動位相補正を計算するための例示的な手順９００を示す。９０２で開始して、複素数ＯＣＴのＡ走査データＡｎ及びＡｎ＋１が受け取られる。ＯＣＴＡの文脈においては、例えば、これら２つの走査Ａｎ及びＡｎ＋１は、同じ場所での２つの繰り返し走査とすることができる。ドップラーイメージングの文脈においては、これら２つの走査は、いくらかの重なり合いがあり得る２つの隣り合う走査とすることができる。走査の一方は参照ＯＣＴ信号と考えられ、他方は、補正又は補償が必要なバルク運動を有する信号である。例えば、参照信号が走査Ａｎと指定された場合、その対応する、補償が必要なバルク運動を含む信号はスキャンＡｎ＋１になる。
【００５０】
９０４で、走査Ａｎ及びＡｎ＋１における信号の分散値は、次式を用いて計算される。
ｖ＝［ｓｔｄ（Ａｎ−Ａｎ＋１・ｅｘｐ（−ｉ・θ））］２    式（１５）
【００５１】
この式（１５）の分散式は次に、θの４つの異なる値において評価される。この特定の例では、９０６、９０８、９１０、及び９１２で、それぞれθがθ＝０、π／２、π、３π／２の値において評価される。これにより、以下の分散値の組が得られる。
【数６】
2021520244.tif 000009
【００５２】
９１４では、分散は上記の式（８）に記述された形を有するという認識に基づく。
ｖ＝Ｃ１−Ｃ２ｃｏｓ（ΔφＢ＋θ）    上記の式（８）
【００５３】
バルク運動位相シフトは、上記の式（１４）を用いて以下のように計算することができる。
【数７】
2021520244.tif 000010
【００５４】
このバルク運動位相シフトΔφＢは次に、補正が必要なＡ走査信号のバルク運動位相を補償するために用いることができる。例えば、９１６に示されるように、参照信号が走査Ａｎと指定された場合に、対応する信号走査Ａｎ＋１のバルク運動位相は、次式を用いて補償することができる。
Ａｎｅｗ＝Ａｎ＋１・ｅｘｐ（ｉ・ΔφＢ）    式（１６）
【００５５】
Ａｎｅｗが計算されると、これを用いてＯＣＴＡ適用の補正を９１８で次式のように、
ＯＣＴＡ＝｜Ａｎ−Ａｎｅｗ｜    式（１７）
又はドップラーイメージング適用では、９２０で次式のように行うことができる。
ｄｏｐｐｌｅｒ＝ａｎｇｌｅ（Ａｎ×ｃｏｎｊ（Ａｎｅｗ））    式（１８）
【００５６】
実施例６  光コヒーレンストモグラフィ血管造影像処理システム
図７は、様々な実施形態による、ＯＣＴ像処理のための例示的なシステム７００を概略的に示す。システム７００は、ＯＣＴインターフェログラムを含むＯＣＴ像を取得するように構成されたＯＣＴシステム７０２と、本明細書に記載の様々な処理ルーチンを実施するように構成されている１つ又は複数のプロセッサ若しくは計算システム７０４とを備える。ＯＣＴシステム７００は、構造ＯＣＴ適用及びＯＣＴ血管造影適用に適した、例えば波長掃引型ＯＣＴ又はスペクトル領域ＯＣＴシステムに適したＯＣＴシステムを備え得る。
【００５７】
様々な実施形態では、ＯＣＴシステムは、操作者が様々な作業を行えるように構成され得る。例えば、１つのＯＣＴシステムが、操作者が本明細書に記載の方法のうちの様々なものを構成及び／又は実践できるように構成され得る。いくつかの実施形態では、１つのＯＣＴシステムが、例えば１つの試料について実行された走査の結果の報告を含む様々な情報の報告を生成するように、又は生成されるように構成され得る。
【００５８】
表示デバイスを備えるＯＣＴシステムの実施形態では、データ及び／又は他の情報を操作者に対し表示することができる。諸実施形態では、表示デバイスは入力を受け取るように構成することができ（例えば、タッチスクリーン、アイコンの作動、ジョイスティック又はつまみ等の入力デバイスの操作等によって）、その入力は、場合によって１つ又は複数のプロセッサに（能動的及び／又は受動的に）伝達することができる。様々な実施形態では、データ及び／又は情報を表示することができ、操作者は、その表示に応じて情報を入力することができる。
【００５９】
いくつかの実施形態では、上述の方法及び処理は、１つ又は複数のコンピュータを含むコンピュータシステムと結び付けることができる。特に、本明細書に記載の方法及び処理、例えば上述の、図１Ａに描写された方法１００、図１Ｂに描写された方法１５０、及び／又は図９に描写された手順９００は、コンピュータアプリケーション、コンピュータサービス、コンピュータＡＰＩ、コンピュータライブラリ、及び／又は他のコンピュータプログラム製品として実施することができる。
【００６０】
図８は、上述の方法及び処理のうちの１つ又は複数を行うことができる、非限定的なコンピュータデバイス８００を概略的に示す。例えば、コンピュータデバイス８００は、上述のシステム７００に含まれるプロセッサであってよく、ＯＣＴシステム又はＯＣＴ像取得装置と動作可能に結合する、通信する、又はこれらに含めることができる。コンピュータデバイス８００は、簡略化された形で示されている。実質的に任意のコンピュータアーキテクチャが、本開示の範囲から逸脱することなく使用され得ることを理解されたい。異なる実施形態において、コンピュータデバイス８００は、マイクロコンピュータ、組み込みコンピュータ回路、プリント回路基板（ＰＣＢ）、マイクロチップ、メインフレームコンピュータ、サーバコンピュータ、デスクトップコンピュータ、ラップトップコンピュータ、タブレットコンピュータ、家庭用娯楽コンピュータ、ネットワークコンピュータデバイス、モバイルコンピュータデバイス、モバイル通信デバイス、ゲームデバイス等の形を取ることができる。
【００６１】
コンピュータデバイス８００は、論理サブシステム８０２及びデータ保持サブシステム８０４を含む。コンピュータデバイス８００は、表示サブシステム８０６、通信サブシステム８０８、撮像サブシステム８１０、及び／又は図８に示されていない他の構成要素を任意選択で含み得る。コンピュータデバイス８００はさらに、例えば手動のボタン、スイッチ、キーボード、マウス、ゲームコントローラ、カメラ、マイクロホン、及び／又はタッチスクリーン等の、ユーザ入力デバイスを任意選択で含み得る。
【００６２】
論理サブシステム８０２は、１つ又は複数の機械可読命令を実行するように構成された、１つ又は複数の物理的デバイスを含み得る。例えば、論理サブシステムは、１つ又は複数のアプリケーション、サービス、プログラム、ルーチン、ライブラリ、オブジェクト、コンポーネント、データ構造、又は他の論理構造の一部である、１つ又は複数の命令を実行するように構成することができる。このような命令は、１つのタスクを実行し、１つのデータ型を実装し、１つ若しくは複数のデバイスの状態を変え、又は別の方法で所望の結果に到達するように実施することができる。
【００６３】
論理サブシステムは、ソフトウェア命令を実行するように構成されている１つ又は複数のプロセッサを含み得る。例えば、１つ又は複数のプロセッサは、本明細書に記載の様々な動作を行うようにプログラムされた物理的回路を備え得る。加えて、又は別法として、論理サブシステムは、ハードウェア又はファームウェア命令を実行するように構成された１つ又は複数のハードウェア若しくはファームウェア論理機械を含み得る。論理サブシステムのプロセッサは、シングルコア又はマルチコアとすることができ、ここで実行されるプログラムは、並行処理又は分散処理用に構成することができる。論理サブシステムは、遠隔に設置できる、及び／又は協調処理用に構成できる２つ以上のデバイス全体にわたって分散している個別の構成要素を任意選択で含み得る。論理サブシステムの１つ又は複数の態様は、クラウドコンピュータ構成の形で構成された遠隔アクセス可能なネットワーク化コンピュータデバイスによって仮想化し、実行することができる。
【００６４】
データ保持サブシステム８０４は、本明細書に記載の方法及び処理を実施するために論理サブシステムによって実行可能なデータ及び／又は命令を保持するように構成された、１つ又は複数の物理的、非一時的なデバイスを含み得る。このような方法及び処理が実施されるとき、データ保持サブシステム８０４の状態は変えることができる（例えば、別のデータを保持するために）。
【００６５】
データ保持サブシステム８０４は、取り外し可能媒体及び／又は内蔵デバイスを含み得る。データ保持サブシステム８０４は、光メモリデバイス（例えば、ＣＤ、ＤＶＤ、ＨＤ−ＤＶＤ、ブルーレイディスク等）、半導体メモリデバイス（例えば、ＲＡＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ等）、及び／又は磁気メモリデバイス（例えば、ハードディスクドライブ、フロッピーディスクドライブ、テープドライブ、ＭＲＡＭ等）をとりわけ含み得る。データ保持サブシステム８０４は、次の特性、すなわち、揮発性、不揮発性、動的、静的、読み出し／書き込み、読み取り専用、ランダムアクセス、順次アクセス、位置アドレス指定可能、ファイルアドレス指定可能、及びコンテンツアドレス指定可能、のうちの１つ又は複数を備えたデバイスを含み得る。いくつかの実施形態では、論理サブシステム８０２及びデータ保持サブシステム８０４は、特定用途向け集積回路又は１つのチップ上のシステム等の、１つ又は複数の共通デバイスに統合することができる。
【００６６】
図８はさらに、取り外し可能なコンピュータ可読記憶媒体８１２の形のデータ保持サブシステムの一態様を示し、このデータ保持サブシステムは、本明細書に記載の方法及び処理を実施するために実行可能なデータ及び／又は命令を記憶及び／又は転送するために使用することができる。取り外し可能なコンピュータ可読記憶媒体８１２は、ＣＤ、ＤＶＤ、ＨＤ−ＤＶＤ、ブルーレイディスク、ＥＥＰＲＯＭ、フラッシュメモリカード、ＵＳＢ記憶デバイス、及び／又はフロッピーディスクの形をとりわけ取り得る。
【００６７】
表示サブシステム８０６は、含まれる場合に、データ保持サブシステム８０４によって保持されたデータの視覚表現を提示するために使用することができる。本明細書に記載の方法及び処理により、データ保持サブシステムによって保持されているデータを変更し、したがってデータ保持サブシステムの状態を変えるので、表示サブシステム８０６の状態が同様に変えられて、内在するデータの変化を視覚的に表すことができる。表示サブシステム８０６は、実質的に任意の種類の技術を利用する１つ又は複数の表示デバイスを含み得る。このような表示デバイスは、論理サブシステム８０２及び／又はデータ保持サブシステム８０４と共用される筐体内で結合することができ、或いはこのような表示デバイスは周辺表示デバイスとすることができる。
【００６８】
通信サブシステム８０８は、含まれる場合に、コンピュータデバイス８００を１つ又は複数の他のコンピュータデバイスと通信可能に結合するように構成することができる。通信サブシステム８０８は、１つ又は複数の異なる通信プロトコルと互換性がある有線及び／又は無線の通信デバイスを含み得る。非限定的な例として、通信サブシステムは、無線電話ネットワーク、無線ローカルエリアネットワーク、有線ローカルエリアネットワーク、無線広域ネットワーク、有線広域ネットワーク等を介する通信用に構成することができる。いくつかの実施形態では、通信サブシステムは、コンピュータデバイス８００がインターネット等を介して他のデバイスとの間でメッセージを送信及び／又は受信することを可能にし得る。
【００６９】
撮像サブシステム８１０は、含まれる場合に、コンピュータデバイス８００と通信する様々なセンサ又は撮像デバイスからの任意の適切な像データを取得及び／又は処理するために使用することができる。例えば、撮像サブシステム８１０は、ＯＣＴシステム、例えば上述のＯＣＴシステム７０２の一部として、ＯＣＴ像データ、例えばインターフェログラムを取得するように構成することができる。撮像サブシステム８１０は、論理サブシステム８０２及び／又はデータ保持サブシステム８０４と共用される筐体内で結合することができ、或いはこのような撮像サブシステムは周辺撮像デバイスを備えることができる。撮像サブシステムから受け取ったデータは、例えば、データ保持サブシステム８０４及び／又は取り外し可能なコンピュータ可読記憶媒体８１２によって保持することができる。
【００７０】
本明細書に記載の構成及び／又は手法は本質的に例示的なものであること、並びに、これらの特定の実施形態又は実施例は、多数の変形形態が可能であるので、限定的な意味に解釈されるべきではないことを理解されたい。本明細書に記載の特定のルーチン又は方法は、任意の数の処理方策のうちの１つ又は複数であってよい。そのため、図示の様々な動作は、図示の順に、別の順に、若しくは並行して実行することができ、又は場合により省略することができる。同様に、上述の処理の順序は変更することができる。
【００７１】
本開示の主題は、本明細書に開示された様々な処理、システムと構成、及び他の特徴、機能、動作、及び／又は特性、並びにこれらすべての等価物の、すべての新規な、また明白ではない組み合わせ、及びサブコンビネーションを含む。
【図１Ａ】
2021520244.tif 000014
【図１Ｂ】
2021520244.tif 000015
【図２Ａ】
2021520244.tif 000016
【図２Ｂ】
2021520244.tif 000017
【図２Ｃ】
2021520244.tif 000018
【図２Ｄ】
2021520244.tif 000019
【図３Ａ】
2021520244.tif 000020
【図３Ｂ】
2021520244.tif 000021
【図３Ｃ】
2021520244.tif 000022
【図３Ｄ】
2021520244.tif 000023
【図３Ｅ】
2021520244.tif 000024
【図３Ｆ】
2021520244.tif 000025
【図４Ａ】
2021520244.tif 000026
【図４Ｂ】
2021520244.tif 000027
【図４Ｃ】
2021520244.tif 000028
【図４Ｄ】
2021520244.tif 000029
【図５Ａ−５Ｌ】
2021520244.tif 000030
【図６】
2021520244.tif 000031
【図７】
2021520244.tif 000032
【図８】
2021520244.tif 000033
【図９】
2021520244.tif 000034
【国際調査報告】
2021520244.tif 000035
2021520244.tif 000036