知財求人 - 知財ポータルサイト「IP Force」
▶ ハーグ－シュトライト アクチェンゲゼルシャフトの特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024092016
(43)【公開日】2024-07-05
(54)【発明の名称】眼科向けの光干渉断層撮影システム
(51)【国際特許分類】
A61B 3/10 20060101AFI20240628BHJP
【ＦＩ】
A61B3/10 100
【審査請求】有
【請求項の数】8
【出願形態】ＯＬ
(21)【出願番号】P 2024074346
(22)【出願日】2024-05-01
(62)【分割の表示】P 2022542052の分割
【原出願日】2020-01-08
(71)【出願人】
【識別番号】521511575
【氏名又は名称】ハーグ－シュトライト アクチェンゲゼルシャフト
(74)【代理人】
【識別番号】100099759
【弁理士】
【氏名又は名称】青木 篤
(74)【代理人】
【識別番号】100123582
【弁理士】
【氏名又は名称】三橋 真二
(74)【代理人】
【識別番号】100092624
【弁理士】
【氏名又は名称】鶴田 準一
(74)【代理人】
【識別番号】100114018
【弁理士】
【氏名又は名称】南山 知広
(74)【代理人】
【識別番号】100153729
【弁理士】
【氏名又は名称】森本 有一
(72)【発明者】
【氏名】ルツィオ ロブレド
(57)【要約】
【課題】 眼球内の構造を表わすデータ、詳細には横断面画像が、各々光の軌跡に沿った複数の点についての複数の反射値（ｒ）を生成する眼球内の異なる位置にある波長掃引型ＯＣＴを用いて複数のＡスキャン（Ａ１、Ａ２、Ａ３）を生成することにより記録される。その後、眼球内の異なる位置におけるいくつかの反射値（ｒ）から、組合せ値（ｒc）が計算される。
【解決手段】 データの質を改善するために、眼球内の少なくとも１つの湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）のモデルが反射値にフィットされ、その後、組み合わせるべきＡスキャン（Ａ１、Ａ２、Ａ３）上の点を識別するために使用される。これにより、ノイズを削減しかつ構造に直交する方向で改善された分解能を得るために、湾曲構造の接線方向に沿ってデータを投影することが可能となる。
【選択図】図５
【解決手段】 データの質を改善するために、眼球内の少なくとも１つの湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）のモデルが反射値にフィットされ、その後、組み合わせるべきＡスキャン（Ａ１、Ａ２、Ａ３）上の点を識別するために使用される。これにより、ノイズを削減しかつ構造に直交する方向で改善された分解能を得るために、湾曲構造の接線方向に沿ってデータを投影することが可能となる。
【選択図】図５
【特許請求の範囲】
【請求項1】
光干渉断層撮影を用いて眼球内の構造を表わすデータを得るための方法において、
－ Ａスキャン（Ａ１、Ａ２、Ａ３）のセットを記録するステップであって、前記Ａスキャン（Ａ１、Ａ２、Ａ３）の少なくともいくつかが眼球の異なる位置で記録され、各々のＡスキャン（Ａ１、Ａ２、Ａ３）が、眼球を通した光の軌跡に沿った複数の点について複数の反射値（ｒ）を生成するために使用される、ステップと、
前記眼球の少なくとも１つの湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒ）を使用するステップと、
－ 複数の組合せ値（ｒc）を計算するステップであって、各々の組合せ値（ｒc）が前記眼球内の異なる位置におけるいくつかの反射値（ｒ）から得られ、前記組合せ値（ｒc）を計算するために使用される前記Ａスキャン（Ａ１、Ａ２、Ａ３）のサブセット上の前記点を識別するために、前記モデリングされた湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）が使用されるステップと、
を含み、前記モデリングされた湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）における組合せ値（ｒc）を計算するために、前記方法は、
－ 前記サブセット内の各々のＡスキャンｉについて、Ａスキャンｉと前記モデリングされた湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）との交差部における前記反射値（ｒ）の推定値ｒe（ｉ）を決定するステップと、
－ 前記組合せ値（ｒc）を計算するために前記推定値ｒe（ｉ）を使用するステップと、
を更に含む、方法。
－ Ａスキャン（Ａ１、Ａ２、Ａ３）のセットを記録するステップであって、前記Ａスキャン（Ａ１、Ａ２、Ａ３）の少なくともいくつかが眼球の異なる位置で記録され、各々のＡスキャン（Ａ１、Ａ２、Ａ３）が、眼球を通した光の軌跡に沿った複数の点について複数の反射値（ｒ）を生成するために使用される、ステップと、
前記眼球の少なくとも１つの湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒ）を使用するステップと、
－ 複数の組合せ値（ｒc）を計算するステップであって、各々の組合せ値（ｒc）が前記眼球内の異なる位置におけるいくつかの反射値（ｒ）から得られ、前記組合せ値（ｒc）を計算するために使用される前記Ａスキャン（Ａ１、Ａ２、Ａ３）のサブセット上の前記点を識別するために、前記モデリングされた湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）が使用されるステップと、
を含み、前記モデリングされた湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）における組合せ値（ｒc）を計算するために、前記方法は、
－ 前記サブセット内の各々のＡスキャンｉについて、Ａスキャンｉと前記モデリングされた湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）との交差部における前記反射値（ｒ）の推定値ｒe（ｉ）を決定するステップと、
－ 前記組合せ値（ｒc）を計算するために前記推定値ｒe（ｉ）を使用するステップと、
を更に含む、方法。
【請求項2】
前記湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）が、角膜の前面（Ｃａ）、角膜の後面（Ｃｂ）、水晶体の前面（Ｌａ）、水晶体の後面（Ｌｐ）、および網膜（Ｒ）のうちの１つである、請求項１に記載の方法。
【請求項3】
所与の深さ座標について、前記深さ座標に最も近い前記Ａスキャンｉ上の少なくとも２つの反射値（ｒ）を補間することによって前記サブセット内のＡスキャンｉの推定値ｒe（ｉ）を計算するステップを含む、請求項１または２に記載の方法。
【請求項4】
所与の軸に沿って前記組合せ値（ｒc）を計算するために、
－ 眼球内の少なくとも第１および第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒc）を使用するステップであって、前記第１の湾曲構造が前記軸に沿った第１の位置に位置設定されており、前記第２の湾曲構造が前記軸に沿った第２の位置に位置設定されている、ステップ
を含む、請求項１から３のいずれか１項に記載の方法。
－ 眼球内の少なくとも第１および第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒc）を使用するステップであって、前記第１の湾曲構造が前記軸に沿った第１の位置に位置設定されており、前記第２の湾曲構造が前記軸に沿った第２の位置に位置設定されている、ステップ
を含む、請求項１から３のいずれか１項に記載の方法。
【請求項5】
－ 前記第１の位置における組合せ値（ｒc）を得るために、前記第１の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）に沿って、点の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの第１の複数の前記反射値（ｒ）を組み合わせるステップ、および／または
－ 前記第２の位置における組合せ値（ｒc）を得るために、前記第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）に沿った点の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの第２の複数の前記反射値（ｒ）を組み合わせるステップ、
のうちの少なくとも１つをさらに含む、請求項４に記載の方法。
－ 前記第２の位置における組合せ値（ｒc）を得るために、前記第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）に沿った点の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの第２の複数の前記反射値（ｒ）を組み合わせるステップ、
のうちの少なくとも１つをさらに含む、請求項４に記載の方法。
【請求項6】
－ 前記第１および前記第２の位置の間の第３の位置における組合せ値（ｒc）を得るために、異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの点の第３の複数の反射値（ｒ）を組み合わせるステップであって、前記点が前記第１および前記第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の間および／または近傍の領域内にあるステップ、
をさらに含む、請求項４または５に記載の方法。
をさらに含む、請求項４または５に記載の方法。
【請求項7】
前記第３の複数の前記反射値に対応する前記点が、
－ 前記第１および前記第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）から補間された曲線（Ｃ３）に沿って配設されているか、または、
－ 前記第１の湾曲構造に平行な第１の曲線（Ｃ１）に沿って配設された点の第１のサブセット、および前記第２の湾曲構造に平行な第２の曲線（Ｃ２）に沿って配設された点の第２のサブセットを含み、ここで前記第１および第２の曲線（Ｃ１、Ｃ２）が前記第３の位置で交差している、
請求項６に記載の方法。
－ 前記第１および前記第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）から補間された曲線（Ｃ３）に沿って配設されているか、または、
－ 前記第１の湾曲構造に平行な第１の曲線（Ｃ１）に沿って配設された点の第１のサブセット、および前記第２の湾曲構造に平行な第２の曲線（Ｃ２）に沿って配設された点の第２のサブセットを含み、ここで前記第１および第２の曲線（Ｃ１、Ｃ２）が前記第３の位置で交差している、
請求項６に記載の方法。
【請求項8】
前記データが、前記構造の断面ビューを示す２次元横断面を表わしており、前記断面ビュー内の前記湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の画像点を生成するために前記組合せ値（ｒc）を使用するステップを含む、請求項１から７のいずれか１項に記載の方法。
【請求項9】
前記組合せ値（ｒc）の少なくともいくつかについて、前記組合せ値（ｒc）が、横断面の平面の外側の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの複数の前記反射値（ｒ）を組み合わせることによって計算される、請求項８に記載の方法。
【請求項10】
複数の組合せ値（ｒc）を計算する前記ステップが、
－ 組み合わせるべき前記反射値（ｒ）の平均を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の分位数を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の最小値または最大値を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の補間を計算するステップ、
のうちの少なくとも１つを含む、請求項１から９のいずれか１項に記載の方法。
－ 組み合わせるべき前記反射値（ｒ）の平均を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の分位数を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の最小値または最大値を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の補間を計算するステップ、
のうちの少なくとも１つを含む、請求項１から９のいずれか１項に記載の方法。
【請求項11】
複数の組合せ値（ｒc）を計算する前記ステップが、前記組合せ値（ｒc）を計算するときに角膜の頂部から閾値Ｒｓよりも小さい距離を有するＡスキャンを使用しないことによって、鏡面反射を示すいかなるＡスキャンの反射値も使用しないステップを含む、請求項１から１０のいずれか１項に記載の方法。
【請求項12】
位置（ｘ、ｙ）における組合せ値（ｒc）を計算するために使用される前記反射値（ｒ）の前記位置が、前記位置（ｘ、ｙ）から１ｍｍ未満の距離（Ｒ）内にあるＡスキャンからである、請求項１から１１のいずれか１項に記載の方法。
【請求項13】
前記データが、画像平面を定義する２次元画像を表わし、前記組合せ値（ｒc）が前記画像内の画素にアトリビュートされた色度であり、
各画素について、Ａスキャンのサブセットを使用し、各Ａスキャンについて、最強の反射値（ｒ）を有する点を識別するステップと；
－ 前記画像平面からの前記識別された点の距離に応じた色を、前記画素について決定するステップと；
を含む、請求項１から１２のいずれか１項に記載の方法。
各画素について、Ａスキャンのサブセットを使用し、各Ａスキャンについて、最強の反射値（ｒ）を有する点を識別するステップと；
－ 前記画像平面からの前記識別された点の距離に応じた色を、前記画素について決定するステップと；
を含む、請求項１から１２のいずれか１項に記載の方法。
【請求項14】
－ 前記眼球の少なくとも１つの湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒ）を使用するステップは、前記眼球のいくつかの湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングすることを含み、
－ 前記方法は、
－ 前記いくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の位置と、前記点の平均深度（ｄ（ｍ、ｎ））とを比較して、前記構造のうちの１つに対する、または前記構造の間、前または後ろの空間に対する前記点のアトリビューションを行なうステップと、
－ 前記アトリビューションに応じて前記色を選択するステップと、
を含む請求項１３に記載の方法。
－ 前記方法は、
－ 前記いくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の位置と、前記点の平均深度（ｄ（ｍ、ｎ））とを比較して、前記構造のうちの１つに対する、または前記構造の間、前または後ろの空間に対する前記点のアトリビューションを行なうステップと、
－ 前記アトリビューションに応じて前記色を選択するステップと、
を含む請求項１３に記載の方法。
【請求項15】
前記画素の輝度または明度を決定するために前記点の反射値（ｒ）の大きさを使用するステップを含む、請求項１３に記載の方法。
【請求項16】
前記データが、前記眼球の一次元横断面を表わすものである、請求項１から１５のいずれか１項に記載の方法。
【請求項17】
前記データが、前記眼球の視軸または光学軸に平行な１次元または２次元横断面を表わすものである、請求項１から１６のいずれか１項に記載の方法。
【請求項18】
前記データが、前記眼球の視軸または光学軸に直交する２次元横断面を表わすものである、請求項１から１６のいずれか１項に記載の方法。
【請求項19】
前記光干渉断層撮影が、周波数ドメインＯＣＴである、請求項１から１８のいずれか１項に記載の方法。
【請求項20】
－ 光干渉断層撮影干渉計（１０～２６）と、
－ 請求項１から１９のいずれか１項に記載の方法のすべてのステップを実施するために構造化され適応された制御ユニット（３２）と、
を含む眼科用デバイス。
－ 請求項１から１９のいずれか１項に記載の方法のすべてのステップを実施するために構造化され適応された制御ユニット（３２）と、
を含む眼科用デバイス。
【請求項21】
前記光干渉断層撮影干渉計（１０～２６）が単一のＡスキャンについて少なくとも４０ｍｍの測定範囲を有する、請求項２０に記載の眼科用デバイス。
【請求項22】
前記光干渉断層撮影干渉計（１０～２６）が、テレセントリックプローブビームを生成する、請求項２０又は２１に記載の眼科用デバイス。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、眼球内の構造を表わすデータを得るための方法ならびに、この方法を実施するための眼科用デバイスに関する。
【背景技術】
【0002】
眼科の利用分野向けの光干渉断層撮影（ＯＣＴ）システムは、例えば欧州特許第３５７２７６５号から公知である。これらのシステムは、光干渉断層撮影干渉計を含む。この干渉計は、眼球の複数のＡスキャンを記録するために使用される。各々のＡスキャンは、眼球を通る光の軌跡に沿った複数の点について一連の反射値からなる。これらの反射値は、その後、例えば１次元、２次元または３次元のいずれかで眼球の画像を生成するために使用可能である。
【発明の概要】
【0003】
本発明が解決すべき課題は、例えばより優れた分解能を有しかつ／またはよりノイズの少ない、より多くの情報を伴うデータを生成する上述のタイプの方法およびデバイスを提供することにある。
【0004】
この課題は、独立クレームの方法およびデバイスによって解決される。
【0005】
したがって、本発明は、光干渉断層撮影を用いて眼球内の構造を表わすデータ、詳細には横断面画像を表わすデータを得るための方法に関する。
【0006】
当該技術によって得られるデータは、例えば眼球の１次元または２次元横断面および／または眼球の個別のパラメータを表わし得る。
【0007】
該方法は、少なくとも以下のステップを含む：
－ Ａスキャンのセットを記録するステップ：これらのＡスキャンの少なくともいくつかは眼球の異なる位置で記録される。各々のＡスキャンは、例えば眼球内の光学的深達度の関数として、複数の反射値を生成するために使用される。
－ 複数の組合せ値を計算するステップ。各々のこのような組合せ値は、眼球内の異なる位置における反射値から得られる。
－ Ａスキャンのセットを記録するステップ：これらのＡスキャンの少なくともいくつかは眼球の異なる位置で記録される。各々のＡスキャンは、例えば眼球内の光学的深達度の関数として、複数の反射値を生成するために使用される。
－ 複数の組合せ値を計算するステップ。各々のこのような組合せ値は、眼球内の異なる位置における反射値から得られる。
【0008】
眼球内の異なる位置で得た反射値を前記組合せ値へと組み合わせることによって、以下で説明するように、追加のまたはより正確な情報を獲得することができる。
【0009】
有利には、該方法はさらに、眼球の少なくとも１つの湾曲構造をモデリングするために反射値を使用するステップを含む。この湾曲構造は、角膜または水晶体または網膜の前面または後面などの、反射値から判別可能なあらゆる構造であり得る。
【0010】
有利には、モデリングされた湾曲構造は、所与の組合せ値を計算するために使用すべきＡスキャンのサブセット上の点を識別するために使用可能である。これにより、例えば、反射値を組合せ値へと組み合わせるとき湾曲構造の接線方向に沿ってＡスキャンのデータを「投影する」ことが可能となり、こうして、湾曲構造の方向に直交する方向での空間分解能が向上しかつ／または、データの信号対ノイズ比が改善される。
【0011】
使用されたＡスキャンのサブセットは、有利には、眼球内の異なる位置において記録されたいくつかのＡスキャンを含む。
【0012】
湾曲構造において組合せ値を計算するために、該方法は有利には、少なくとも以下のステップを含み得る：
－ 使用すべきＡスキャンのサブグループにおいて、各々のＡスキャンｉについて、Ａスキャンｉとモデリングされた構造との交差部における反射値の推定値ｒe（ｉ）を決定するステップ。この推定値は例えば、交差部に最も近いＡスキャンｉ上の測定された点の反射値であり得、あるいは交差部の位置におけるＡスキャンｉの反射値の補間値であってもよい。
－ 組合せ値を計算するために推定値ｒe（ｉ）を使用するステップ。
－ 使用すべきＡスキャンのサブグループにおいて、各々のＡスキャンｉについて、Ａスキャンｉとモデリングされた構造との交差部における反射値の推定値ｒe（ｉ）を決定するステップ。この推定値は例えば、交差部に最も近いＡスキャンｉ上の測定された点の反射値であり得、あるいは交差部の位置におけるＡスキャンｉの反射値の補間値であってもよい。
－ 組合せ値を計算するために推定値ｒe（ｉ）を使用するステップ。
【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】
デバイスはさらに、例えばマイクロプロセッサ３４およびメモリ３６ならびにディスプレイ３８が具備され得る制御ユニット３２を含む。メモリ３６は、該方法のステップを実施するために必要とされるデータならびにプログラム命令を保持し得る。ディスプレイ３８は例えば、それによって決定されたデータを示すため、詳細には本明細書中に記載の技術を用いて導出されたあらゆる横断面画像を表示するために使用され得る。
【0041】
有利には、ＯＣＴ干渉計１０～２６の測定範囲（単一のＡスキャンについて）は、少なくとも典型的な眼球の角膜から網膜まで延在する。換言すると、単一のＡスキャンで（すなわち、光源の一回の掃引でのＳＳ－ＯＣＴについて）、少なくとも４０ｍｍ（空中）にわたって深度分解された情報を得ることができる。これにより、例えば異なる測定値を組み合わせるためにステッチングを適用する必要なく、軸方向の眼球の長さ全体にわたり、以下で説明する技術を適用することが可能になる。
【0042】
図２は、測定において使用されるスキャンパターンについての一例を示す。すなわち、これは、さまざまなＡスキャン中のｘ－ｙ平面内のプローブビーム３８の位置を示す。このタイプのパターンは、欧州特許第３０２１０７１号に記載されている。例えば、欧州特許第３２１７１４４号または米国特許第８７０５０４８号に記載のスキャンパターンなどの、他のスキャンパターンも同様に使用可能である。
【0043】
図３は、図２に示されている脂肪スキャン点についてのｚ座標およびスキャン点指数の関数としての反射値ｒ（未加工データ）を示す。この図面中、より高い反射値はより暗いスポットで表わされている。見れば分かるように、このデータはかなりノイズが多い。
【0044】
眼球運動補正
このようにして記録されたＡスキャンは、任意には眼球運動について、例えば少なくとも以下のステップを使用することによって補正可能である：
１． Ａスキャン内で少なくとも１つの所与の眼球構造（例えば前角膜面）の反射を識別するステップ、
２． 構造の形状および構造の運動を描写するモデルを、識別された反射の位置にフィットさせるステップ。このモデルは、例えば構造の幾何学的パラメータ（例えば曲率）ならびに運動パラメータ（例えばｘ、ｙおよびｚ座標での３次元の位置および速度）を有することができる。
このようにして記録されたＡスキャンは、任意には眼球運動について、例えば少なくとも以下のステップを使用することによって補正可能である：
１． Ａスキャン内で少なくとも１つの所与の眼球構造（例えば前角膜面）の反射を識別するステップ、
２． 構造の形状および構造の運動を描写するモデルを、識別された反射の位置にフィットさせるステップ。このモデルは、例えば構造の幾何学的パラメータ（例えば曲率）ならびに運動パラメータ（例えばｘ、ｙおよびｚ座標での３次元の位置および速度）を有することができる。
【0045】
フィッティングステップ２において得られたパラメータは次に、ＯＣＴ測定値を眼球の枠で固定される座標系へと変換するために使用可能である。
【0046】
好適な運動補正技術は、例えば国際公開第２０１３／１０７６４９号または米国特許第７４５２０７７号中に記載されている。
【0047】
Ａスキャンの組合せ
上述のステップは、眼球の異なるｘ－ｙの位置において記録されたＡスキャンセットを生成する。各Ａスキャンは、眼球内の光の軌跡に沿った複数の反射値ｒ（ｚij）を含み、ここでｉはＡスキャンの指数（Ａスキャンのｘおよびｙ位置を識別する）であり、ｚi1・・・ｚiNは、Ａスキャンｉ内に記録された通りの眼球の深度に沿ったｚ位置を示す。
上述のステップは、眼球の異なるｘ－ｙの位置において記録されたＡスキャンセットを生成する。各Ａスキャンは、眼球内の光の軌跡に沿った複数の反射値ｒ（ｚij）を含み、ここでｉはＡスキャンの指数（Ａスキャンのｘおよびｙ位置を識別する）であり、ｚi1・・・ｚiNは、Ａスキャンｉ内に記録された通りの眼球の深度に沿ったｚ位置を示す。
【0048】
前述したように、反射値ｒ（ｚij）は、絶対必要というわけではないものの有利には、運動補正される。
【0049】
本技術は、眼球の構造を描写するデータの獲得に基づいている。データは、例えば、眼球の１次元または２次元横断面を表わすことができ、または、例えば眼球の３次元ボリュメトリック画像を表わすこともできる。
【0050】
所与の位置における眼球を表わすデータを計算する目的で、本技術は、「組合せ値」ｒcを計算し、ここで、これらの組合せ値の少なくともいくつかは各々、眼球の異なる位置について記録されたいくつかの反射値ｒ（ｚij）から得られるものである。
【0051】
以下では、このような組合せ値を計算するためのいくつかの方法について説明する。
【0052】
眼球の構造に沿った投影
ここで説明される技術の１つは、眼球の少なくとも１つの湾曲構造をモデリングし、その後、この湾曲構造に沿って接線方向に異なるＡスキャンの反射値を組み合わせることに基づいている。
ここで説明される技術の１つは、眼球の少なくとも１つの湾曲構造をモデリングし、その後、この湾曲構造に沿って接線方向に異なるＡスキャンの反射値を組み合わせることに基づいている。
【0053】
これは、平面内の一連のＡスキャンを表わす図４に例示されている。画像は明らかに、角膜の前面および後面Ｃａ、Ｃｐならびに、水晶体の前面および後面Ｌａ、Ｌｐを示している。それは同様に、網膜Ｒも示している。（眼球の視軸または光学軸に対して平行に（眼球の外側で）延在するテレセントリックＡスキャンについて、Ａスキャンは、網膜の１つの点において交差することから、網膜は、全てのＡスキャンについて実質的に同じｚ座標で出現する）。
【0054】
例えば各Ａスキャンに対しエッジ検出技術を用いて、これらの表面のいずれかのｚ位置（光の軌跡に沿った軸方向位置）を推定して、例えばこれらの表面の各々に対応するｘ－ｙ－ｚ空間内の４つの点セットを得ることができる。
【0055】
次のステップでは、これらの表面の各々のモデルのパラメータを、例えば従来のフィッティング技術を用いて計算することができる。
【0056】
例えば、表面は、中心と半径によって定義される球面によってモデリングされ得、あるいは、例えば二次の多項式表面によってこれらをモデリングすることができる。
【0057】
モデルのパラメータがひとたび計算されたならば、所与のｘ－ｙ座標についての一つの表面における精緻な「組合せ値」を計算すること、すなわち前記表面に沿ったいくつかのＡスキャンの反射値を投影し、例えば平均化または以下でさらに説明される別の数学的演算によってこれらの反射値を数学的に組み合わせることによって計算することができる。
【0058】
このプロセスは、フィッティングによって得られる４つの表面Ｃａ、Ｃｐ、ＬａおよびＬｐならびに３つのＡスキャンＡ１、Ａ２、Ａ３を示す図５の中で、図式的に例示されている。Ａスキャンに沿ったドット／円は、反射値ｒ（ｚij）が測定されたｚ位置ｚijを表わす。
【0059】
例えば、ＡスキャンＡ２と表面Ｃａとの交差部で精緻な組合せ反射値ｒcを得るために、表面Ｃａに沿った３つのＡスキャンＡ１、Ａ２、Ａ３の反射値を組み合わせることができる。
【0060】
例えば、単純なアルゴリズムが、各Ａスキャンについて、表面Ｃａに最も近い測定位置（すなわちｚ位置）ＬＣａ１、ＬＣａ２、ＬＣａ３を識別し、それらの反射値ｒ（ＬＣａ１）、ｒ（ＬＣａ２）、ｒ（ＬＣａ３）を組み合わせることができる。
【0061】
より精緻なアルゴリズムは、表面Ｃａに最も近い位置ＬＣａ１、ＬＣａ２、ＬＣａ３を識別するだけでなく、例えば２番目に近い位置をも識別し、各Ａスキャンについて、これらの最も近い位置の反射値を補間して、Ａスキャンと表面Ｃａの交差部における反射値のより優れた推定値を得、その後これらの推定値を数学的に組合せることができる。
【0062】
より一般的に言うと、眼球の構造（例えば表面Ｃａ、Ｃｐ、Ｌａ、Ｌｐのうちの１つ）の上の所与の位置ｘ、ｙにおける反射値ｒcを推定するためには、以下のステップを実施することができる：
－ ステップ１（モデリング）。Ａスキャンにより測定された反射値に対して眼球の少なくとも１つの湾曲構造をモデリングする。このモデリングステップのためには、以下のステップ２のＡスキャンサブセットのみを使用してよく、あるいは、より大きなＡスキャンセットを使用してもよい。
－ ステップ２（サブセットの選択）：組合せ値ｒcを計算するために位置（軸）ｘ、ｙのまわりのＡスキャンのサブセットを選択する。このサブセットは、ｘ、ｙに近いＡスキャン、例えばｘ、ｙから値Ｒよりも小さい距離を有するＡスキャンを含む。Ｒが大きいものとして選択された場合、ノイズ抑制が改善されるが、ｘおよびｙの分解能は劣化する。Ｒが小さいものとして選択された場合、ノイズは増大するものの、分解能は向上する。典型的には、Ｒは例えば１ｍｍまたは０．５ｍｍであってよい。
－ ステップ３（個別のＡスキャン反射推定値を決定する）：サブセット内の各Ａスキャンｉについて、モデリングされた構造とＡスキャンの交差部における反射値の推定値ｒe（ｉ）を計算する。このステップには、例えばモデル構造とＡスキャンｉの交差点の計算、そして次に、例えば推定値ｒe（ｉ）としての交差部に最も近い点についての測定反射値の使用、が含まれ得る。あるいは、このステップには、（間接的に、少なくとも交差点に最も近いスキャン位置における反射値ならびに少なくとも１つのさらなる反射値も同様に使用することになる）交差点におけるＡスキャンの反射値の補間値を計算することが含まれ得る。
４． サブセット内のＡスキャンの推定値ｒe（ｉ）を、組合せ値ｒcへと組み合わせるステップ。この組合せは、以下でより詳細に説明するさまざまな数学的演算を使用することができる。
－ ステップ１（モデリング）。Ａスキャンにより測定された反射値に対して眼球の少なくとも１つの湾曲構造をモデリングする。このモデリングステップのためには、以下のステップ２のＡスキャンサブセットのみを使用してよく、あるいは、より大きなＡスキャンセットを使用してもよい。
－ ステップ２（サブセットの選択）：組合せ値ｒcを計算するために位置（軸）ｘ、ｙのまわりのＡスキャンのサブセットを選択する。このサブセットは、ｘ、ｙに近いＡスキャン、例えばｘ、ｙから値Ｒよりも小さい距離を有するＡスキャンを含む。Ｒが大きいものとして選択された場合、ノイズ抑制が改善されるが、ｘおよびｙの分解能は劣化する。Ｒが小さいものとして選択された場合、ノイズは増大するものの、分解能は向上する。典型的には、Ｒは例えば１ｍｍまたは０．５ｍｍであってよい。
－ ステップ３（個別のＡスキャン反射推定値を決定する）：サブセット内の各Ａスキャンｉについて、モデリングされた構造とＡスキャンの交差部における反射値の推定値ｒe（ｉ）を計算する。このステップには、例えばモデル構造とＡスキャンｉの交差点の計算、そして次に、例えば推定値ｒe（ｉ）としての交差部に最も近い点についての測定反射値の使用、が含まれ得る。あるいは、このステップには、（間接的に、少なくとも交差点に最も近いスキャン位置における反射値ならびに少なくとも１つのさらなる反射値も同様に使用することになる）交差点におけるＡスキャンの反射値の補間値を計算することが含まれ得る。
４． サブセット内のＡスキャンの推定値ｒe（ｉ）を、組合せ値ｒcへと組み合わせるステップ。この組合せは、以下でより詳細に説明するさまざまな数学的演算を使用することができる。
【0063】
上述のステップは、所与の眼球の構造における組合せ値ｒcをどのように計算するかを説明している。これらのステップは、同様に眼球のあらゆる構造から離れた任意の所与の位置ｘ、ｙ、ｚ、例えば図５の位置ｘ、ｙ、ｚ１における反射値ｒcを推定するために一般化することが可能である。
【0064】
第１の方法は、図６および７に例示されている。
【0065】
図６は、図４のＡスキャンを示すが、各々のＡスキャンは、水晶体の前面Ｌａに対応するそれらの反射ピークが同じｚ位置に来るように方向ｚに沿ってオフセットされている。同様にして、図７は、図４のＡスキャンを示すが、各Ａスキャンは、水晶体の前面Ｌｐに対応するそれらの反射ピークが同じｚ位置に来るように方向ｚに沿ってオフセットされている。
【0066】
その場合、水晶体の前面および後面Ｌａ、Ｌｐの間の所与のｚ位置における組合せ値ｒc（ｚ）を計算するために、所与のＡスキャンｉについての推定値ｒe（ｉ）は例えば、図６のＡスキャンｉの反射値ｒ’（ｉ、ｚ）および図７の反射値ｒ’’（ｉ、ｚ）Ａスキャンｉの加重平均であり得る。重みは、一定の距離を伴って、それぞれＬａおよびＬｐからの点ｚの距離ｄ１、ｄ２に応じて線形的に選択される、すなわち：
【数1】
【0067】
換言すると、２つのＡスキャンセットは、第１のセットが第１の構造とアライン整列され（being aligned）、第２のセットが第２の構造と整列された（being aligned）状態で、ｚに沿ってＡスキャンをオフセットすることによって計算される。その後、これらの構造の間の所与の位置ｚにおける所与のＡスキャンｉについての推定値ｒe（ｉ）を計算するために、等式（１）を用いて、２つのセット内のＡスキャンｉの反射値ｒ’およびｒ’’が組み合わされる。
【0068】
図８には、第２の方法が例示されている。ここでは、所与の位置ｘ、ｙ、ｚについての組合せ値は、ｘ、ｙ、ｚに近いＡスキャンのサブセットから計算されるものと仮定されている。図は、Ａスキャンｉ（Ａｉとして記されている）に沿ったどの位置ｚ（ｉ）を使用すべきかを例示している。すなわち、Ａスキャンｉのために使用すべき推定値ｒe（ｉ）は、以下のように決定される：
－ 位置ｘ、ｙにおける２つの眼球構造（ここでは構造ＬａおよびＬｐ）のｚ位置ｚ１およびｚ２が、Ａスキャンｉのｘ－ｙの位置における２つの眼球構造のｚ位置ｚ１（ｉ）およびｚ２（ｉ）においてと同様に計算される。
－ ｚ１およびｚ２、ｚ１（ｉ）およびｚ２（ｉ）からＡスキャンｉ上の位置ｚ（ｉ）が計算され、これには、ｚがｚ１およびｚ２の間にあるのと同様にして、ｚ（ｉ）がｚ１（ｉ）とｚ２（ｉ）の間になければならないということが考慮に入れられる。例えば、ｚ（ｉ）を計算するために、以下の関係式を使用することができる：
－ 次に、位置ｚ（ｉ）における推定値ｒe（ｉ）を計算することができる。例えば、ｚ（ｉ）に最も近い測定点におけるＡスキャンｉの反射値を使用してこれを計算することができ、あるいは、位置ｚ（ｉ）におけるＡスキャンｉの反射値の補間からそれを計算することができる。
－ 位置ｘ、ｙにおける２つの眼球構造（ここでは構造ＬａおよびＬｐ）のｚ位置ｚ１およびｚ２が、Ａスキャンｉのｘ－ｙの位置における２つの眼球構造のｚ位置ｚ１（ｉ）およびｚ２（ｉ）においてと同様に計算される。
－ ｚ１およびｚ２、ｚ１（ｉ）およびｚ２（ｉ）からＡスキャンｉ上の位置ｚ（ｉ）が計算され、これには、ｚがｚ１およびｚ２の間にあるのと同様にして、ｚ（ｉ）がｚ１（ｉ）とｚ２（ｉ）の間になければならないということが考慮に入れられる。例えば、ｚ（ｉ）を計算するために、以下の関係式を使用することができる：
【数2】
【0069】
さらに一般的に言うと、位置ｚ１とｚ２の間の位置ｚにおける組合せ値を計算するために、該方法は、眼球の２つのモデリングされた構造の間の領域内および／またはその近傍にあるＡスキャン点の「第３の複数の」反射値（ｒ）を組み合わせることができる。上述の通り、「～の近傍に」なる用語は有利には、それぞれの構造から５００μｍ未満、詳細には１００μｍ未満、詳細には１０μｍ未満の距離を有する点を指すものとして理解される。「～の近傍に」のケースは、例えば以下のｂ）で定義されている技術を使用する場合にあてはまる可能性があり、その場合、「第３の複数の」値のうちのいくつかは、２つのモデリングされた構造の外側に存在し得る。（このような構造は、８ｍｍ以上の典型的な曲率半径を有することができ、第３の複数の反射値に対応する点の距離は、最大ほぼ３ｍｍ、恐らくはそれより短いものであり得、その結果として５００μｍという「近傍」の最大定義がもたらされる。）
【0070】
「第３の複数」とは、例えば、以下のうちの１つの形で定義され得る：
ａ） 等式（１）の技術を使用する場合、この「第３の複数」の反射値は例えば、使用されたＡスキャン上の点の第１および第２のサブセットの反射値を含むものとして定義され得る：すなわち
－ 第１のサブセットの点は、第１のモデリングされた構造（図８中のＬａ）に平行な第１の（有利には３次元の）曲線Ｃ１に沿って接線方向に配設され、この第１の曲線は座標ｘ、ｙ、ｚにある点と交差するようにオフセットされている。
－ 第２のサブセットの点は、第２のモデリングされた構造（図８中のＬｐ）に平行な第２の（有利には３次元の）曲線Ｃ２に沿って接線方向に配設され、この第２の曲線は同様に、座標ｘ、ｙ、ｚにある点と交差するようにオフセットされている。
ａ） 等式（１）の技術を使用する場合、この「第３の複数」の反射値は例えば、使用されたＡスキャン上の点の第１および第２のサブセットの反射値を含むものとして定義され得る：すなわち
－ 第１のサブセットの点は、第１のモデリングされた構造（図８中のＬａ）に平行な第１の（有利には３次元の）曲線Ｃ１に沿って接線方向に配設され、この第１の曲線は座標ｘ、ｙ、ｚにある点と交差するようにオフセットされている。
－ 第２のサブセットの点は、第２のモデリングされた構造（図８中のＬｐ）に平行な第２の（有利には３次元の）曲線Ｃ２に沿って接線方向に配設され、この第２の曲線は同様に、座標ｘ、ｙ、ｚにある点と交差するようにオフセットされている。
【0071】
換言すると、第１および第２の曲線Ｃ１、Ｃ２は点ｚで交差する。
【0072】
ｂ） 等式（２）の技術を使用する場合、「第３の複数の」反射値は、例えば、第１および第２のモデリングされた構造（図８中のＬａ、Ｌｐ）から補間された曲線Ｃ３（有利には３次元曲線）に沿って接線方向に延在する使用されたＡスキャン上の点についての反射値として定義され得る。
【0073】
この「眼球構造に沿った投影」の節の中で説明されている技術は、眼球の頂点または視軸または光学軸から離れたｘ－ｙ位置における組合せ値ｒcを計算するため、および／または、Ａスキャンの幾何学的に広いサブセット上の組合せ、すなわち大きなＲ値についての組合せを計算するために極めて有利であり、その理由はこれらの場合には、そうしなければ構造の曲率および／または勾配によって組合せ値の計算の際に大きな誤差が導入されると考えられるからである。
【0074】
組合せ値の計算
Ａスキャンサブセット内の全てのＡスキャンｉの推定値ｒe（ｉ）がひとたび決定されたならば、位置ｘ、ｙ、ｚにおける組合せ値ｒcを評価することができる。上述のように、組合せ値ｒcは、例えば以下の値のうちの１つであり得る：
－ 個別のＡスキャンの推定値ｒe（ｉ）の平均。有利には、加重平均が使用され、ここで、所与の推定値ｒe（ｉ）の重みは、ガウスの重みを用いて位置ｘ、ｙからのＡスキャンｉの距離により左右される。
－ 推定値ｒe（ｉ）の分位数、例えば中央値。
－ 推定値ｒe（ｉ）の最小値または最大値。
－ 独立変数としてＡスキャンｉのｘおよびｙ座標を用いて、推定値ｒe（ｉ）を通して一次または二次の多項式ｒc（ｘ、ｙ）をフィットさせることなどによる、推定値ｒe（ｉ）の補間。
Ａスキャンサブセット内の全てのＡスキャンｉの推定値ｒe（ｉ）がひとたび決定されたならば、位置ｘ、ｙ、ｚにおける組合せ値ｒcを評価することができる。上述のように、組合せ値ｒcは、例えば以下の値のうちの１つであり得る：
－ 個別のＡスキャンの推定値ｒe（ｉ）の平均。有利には、加重平均が使用され、ここで、所与の推定値ｒe（ｉ）の重みは、ガウスの重みを用いて位置ｘ、ｙからのＡスキャンｉの距離により左右される。
－ 推定値ｒe（ｉ）の分位数、例えば中央値。
－ 推定値ｒe（ｉ）の最小値または最大値。
－ 独立変数としてＡスキャンｉのｘおよびｙ座標を用いて、推定値ｒe（ｉ）を通して一次または二次の多項式ｒc（ｘ、ｙ）をフィットさせることなどによる、推定値ｒe（ｉ）の補間。
【0075】
有利には、デバイスは、ユーザが推定値ｒe（ｉ）をどのように組合せるかを選択できるようにするユーザインタフェースを含む。ユーザは、分位点の表示を望む場合、例えば、最小値、中央値および最大値の間の値を連続的にまたは段階的に選択できる可能性がある。
【0076】
鏡面反射の抑制
典型的には、眼球の頂点に近いＡスキャンは、考慮対象のサブセット内の他のＡスキャンから獲得した情報をオーバライドする傾向を有し得る強い鏡面反射を示す。
典型的には、眼球の頂点に近いＡスキャンは、考慮対象のサブセット内の他のＡスキャンから獲得した情報をオーバライドする傾向を有し得る強い鏡面反射を示す。
【0077】
したがって、一実施形態においては、組合せ値を計算するときに、鏡面反射を示すあらゆるＡスキャンの反射値を抑制することができる。これは、例えばＲｓが０．１ｍｍなどであるものとして、角膜の頂部から閾値Ｒｓよりも小さい距離を有するあらゆるＡスキャンを使用しないことによって達成可能である。
【0078】
横断面の計算
上述の技術は、例えば、眼球の１次元または２次元横断面を計算するために使用可能である。
上述の技術は、例えば、眼球の１次元または２次元横断面を計算するために使用可能である。
【0079】
特に有利な横断面は、例えば、眼球の視軸または光学軸に沿った一次元横断面、すなわち中央Ａスキャンである。これは、例えば、頂点から例えば１ｍｍ未満の横方向距離を有する全てのＡスキャンのサブセットを使用することによって得ることができ、任意には、例えば頂点から０．１ｍｍより近いあらゆるＡスキャンを排除して、上述の通りの鏡面反射信号を抑制することによって得ることができる。
【0080】
このような横断面を、組合せ値ｒc（ｚ）のグラフとして描いたものが、図９に一例として示されている。この図の中のドットは、個別のＡスキャンを示し、一方、連続線は、頂点から０．１～１．０ｍｍの距離を有する図３のこれらのＡスキャンについての図６および７に例示された方法によって得られた組合せ値ｒc（ｚ）である。ここで分かるように、組合せ値ｒcは、個別のＡスキャンよりもはるかに低いノイズを有する。
【0081】
図１０は、眼球の前部のみを示す、眼球の視軸に沿った別の一次元横断面である。プロットには、次の２つの曲線が含まれる：
－ 点曲線は、角膜または水晶体の湾曲した表面に沿ったＡスキャンの反射率値を投影せずにサブセット内のＡスキャンから得た組合せ値ｒc（ｚ）を示す。むしろ、組合せ値ｒc（ｚ）は、同じｚ座標におけるＡスキャンのそれぞれの値を（例えば中央値を使用することによって）単に組み合わせることによって得られた（すなわちサブグループ内のＡスキャンは、眼球の視軸上に直角に投影された）。
－ 実曲線は、上述の等式（１）を用いた組合せおよび湾曲した表面に沿って投影技術を用いてサブセット内のＡスキャンから得られた組合せ値ｒc（ｚ）を示す。
－ 点曲線は、角膜または水晶体の湾曲した表面に沿ったＡスキャンの反射率値を投影せずにサブセット内のＡスキャンから得た組合せ値ｒc（ｚ）を示す。むしろ、組合せ値ｒc（ｚ）は、同じｚ座標におけるＡスキャンのそれぞれの値を（例えば中央値を使用することによって）単に組み合わせることによって得られた（すなわちサブグループ内のＡスキャンは、眼球の視軸上に直角に投影された）。
－ 実曲線は、上述の等式（１）を用いた組合せおよび湾曲した表面に沿って投影技術を用いてサブセット内のＡスキャンから得られた組合せ値ｒc（ｚ）を示す。
【0082】
【0083】
さらに、点線は、表面Ｃａ、Ｃｐ、およびＬａについては右に、そして表面Ｌｐについては左にオフセットされている。このことは、表面の頂部で使用されている場合に「直角」投影が系統的なオフセット誤差を生成するという事実を反映している。補正後の（確実な）データはこれらの誤差を示さない。
【0084】
図９および１０の断面データは、角膜の深度、眼球の前房の深度、水晶体の厚み、眼球の全長および／または眼球の他の構造の寸法、例えば前水晶体嚢複合体または角膜上皮の厚みなどを測定するために使用可能である。優れたＳ／Ｎ比により正確な結果を得ることが可能である。個別の散乱構造の最大値を特定するために、ピークフィッティングを使用することができる。
【0085】
眼球の視軸または光学軸以外の点について、すなわち任意のｘ－ｙ位置について、このタイプのさらなる一次元スキャンを得ることができ、これらのスキャンは、例えば角膜の厚み、角膜上皮の厚みおよび／または前房の深度の色分けされた深度図を生成するために使用可能である。
【0086】
このタイプのいくつかの一次元スキャンは同様に、例えば子午横断面画像の生成のために、眼球の視軸または光学軸に平行な横断面平面についての二次元横断面データを生成するため、組み合わせ可能である。
【0087】
例えばグレーレベル値として組合せ値ｒcをエンコードすることによって、横断面画像を生成するために２次元横断面データを使用することができる。
【0088】
有利には、このような横断面データを生成するとき、「組合せ値の計算」の節で説明した他の技術が同様に使用可能であっても、組合せ値ｒcは、使用されるサブグループ内のＡスキャンの対応する位置の反射値の最大値から計算される。
【0089】
図１１は、平均化を使用したこのような（ジャーヴィスハーフトーニングされた）横断面画像を示す。画像は、Ａスキャンが眼球を通過するにつれてのそれらの屈折を考慮に入れて、屈折補正されている。
【0090】
２次元横断面データを生成するとき、組み合わせるべきＡスキャンのサブグループの選択は、横断面平面との関係において非対称であり得る。例えば、水平横断面（すなわちｙ＝０）を生成する場合、一定の座標ｘ０における組合せ値ｒcは、Δｙ＞Δｘ、例えばΔｙ＝２ｍｍおよびΔｘ＝０．５ｍｍとして、エリアｘ＝ｘ０±Δｘおよびｙ＝±Δｙ内部のＡスキャンのサブグループを用いて計算可能である。
【0091】
より一般的に言うと、所与の横断面平面についての横断面データを計算する場合、所与の座標ｘ、ｙ、ｚにおける組合せ値ｒcを計算するために使用されるＡスキャンのサブグループは、平面に平行ないずれの方向よりも平面に直交する方向でより大きい空間体積内のＡスキャンのみを含む。
【0092】
深度のエンコーディング
本発明のさらに別の態様は、画像平面からの構造の距離をエンコードするために色が用いられる、カラーエンコーディングされた画像の生成に関する。
本発明のさらに別の態様は、画像平面からの構造の距離をエンコードするために色が用いられる、カラーエンコーディングされた画像の生成に関する。
【0093】
この方法は、以下で、眼球の「正面」画像、すなわち、前方から見た（すなわち方向ｚに沿って外側から見た）眼球を表わす画像について説明される。しかしながら、これは同様に、あらゆるタイプの断面画像を生成するためにも使用可能である。
【0094】
当該概念の普遍性を限定することなく、１つの画像についての「画像平面」は、ビューイング方向に直交するあらゆる平面であり得る。
【0095】
上述のものと同様、この態様において、該方法は、画像内の各点について、Ａスキャンまたはその上のデータ点のサブセットを使用し、それらの反射値を組み合わせて、この場合色度である組合せ値を生成する。
【0096】
結果としての正面画像の一例が、図１２および１３に示されている。実際には、すなわちデバイスの電子ディスプレイ３８上に表示されているように、この画像は、カラー画像となる。しかしながら、本明細書中のそれを表わすために、図１２は、グレーレベル（すなわち輝度または明度）を示し、図１３は、その一定のエリア内で使用される可能性のある色を（文字で）示している。
【0097】
ｘ－ｙ平面内にある、すなわちＡスキャンの軸に直交する平面内にあるこの画像は、例えば以下のように計算可能である：
１． 画像の画素を表わす３つの二次元アレイがセットアップされ、ゼロに初期化される。第１のアレイＤ（ｍ、ｎ）は、距離値を記憶する。第２のアレイＲ（ｍ、ｎ）は、反射率値を記憶する。第３のアレイＷ（ｍ、ｎ）は重み値を記憶する。
１． 画像の画素を表わす３つの二次元アレイがセットアップされ、ゼロに初期化される。第１のアレイＤ（ｍ、ｎ）は、距離値を記憶する。第２のアレイＲ（ｍ、ｎ）は、反射率値を記憶する。第３のアレイＷ（ｍ、ｎ）は重み値を記憶する。
【0098】
２． 問題の各Ａスキャンｉについて、以下の演算が実施される：
２ａ． Ａスキャンｉ上の最大の反射値の振幅ｒ（ｉ）およびｚ位置ｚ（ｉ）を決定する。
２ｂ． 少なくとも１つの画像画素のサブセット全体にわたって反復する。各々の画像画素ｍ、ｎ（ｘ－ｙ平面内の位置ｘ（ｍ）、ｙ（ｍ）に対応する）について、重みｗが計算される。例えば、
ここで、ｘ（ｉ）、ｙ（ｉ）は、Ａスキャンｉの中心のｘおよびｙ座標であり、ｇ（ｒ、σ）は所与の分散σを伴うガウス関数である。分散σは、Ａスキャンｉの直径により左右される。
２ａ． Ａスキャンｉ上の最大の反射値の振幅ｒ（ｉ）およびｚ位置ｚ（ｉ）を決定する。
２ｂ． 少なくとも１つの画像画素のサブセット全体にわたって反復する。各々の画像画素ｍ、ｎ（ｘ－ｙ平面内の位置ｘ（ｍ）、ｙ（ｍ）に対応する）について、重みｗが計算される。例えば、
【数3】
【0099】
重みｗ（ｍ、ｎ）は、第３のアレイＷ（ｍ、ｎ）に加えられ、ｗ（ｍ、ｎ）・ｚ（ｉ）は第１のアレイＤ（ｍ、ｎ）に加えられる。ｗ（ｍ、ｎ）・ｒ（ｉ）は第２のアレイＲ（ｍ、ｎ）に加えられる。
【0100】
３． 全てのＡスキャンｉ（またはＡスキャンの所望のサブセット）がひとたびこのように処理された時点で、以下の演算が、画像の各々の画素ｍ、ｎについて実施される：
３ａ）
が、例えば
を用いて決定される。式中、εはゼロによる除算を回避するための小さい数である。
３ｂ）
が、例えば
を用いて決定される。
３ｃ）
は、例えば、色または色相を決定するために使用される。これは、例えば、平均深度の関数として色を提供するルックアップテーブルを用いて実装可能である。ＲＧＢ色空間の中で、この色は、例えばそのＲＧＢ成分によって表わされ得る。
３ｄ）
は、例えば、色または色相の輝度を設定するために使用される。ＲＧＢ色空間において、これは例えば、平均反射率をＲＧＢ成分に乗じ、その後ＲＧＢ値を正規化することによって実装され得る。
３ａ）
【数4】
【数5】
３ｂ）
【数6】
【数7】
３ｃ）
【数8】
３ｄ）
【数9】
【0101】
ステップ３ｃおよび３ｄの手法による色は、画素ｎ、ｍにアトリビュートされる。
【0102】
例えば、図１２および１３の例において、画像の上部部分内の構造は、まぶたからの反射である。これらは、ビューワに最も近いことから、第１の色１（または色範囲）を有することになる。虹彩に由来する反射はビューワからさらに離れており、第２の色２の中に現われる。
【0103】
瞳孔内に見える反射は、例えば角膜、前水晶体面、後水晶体面、または眼球内のより深い構造における反射に由来する。これらの反射は、それらがどこから来るものであるかに応じて、異なる色３、４などを担持することになる。
【0104】
したがって、図１２および１３に示された画像は、眼球内の異なる散乱構造を容易に区別できるようにする。
【0105】
同じ技術は、眼球の正面ビューについてのみならず、他の方向からの仮想ビューについても活用可能である。その場合、上述のステップ２は、Ａスキャンについて実施されず、例えば画像平面に直交して延在する円筒形または角柱領域について実施される。
【0106】
註：平均反射率を使用する代りに、以上のステップ３において、該方法は同様に、例えば最大反射率およびその深度などの使用といった他の任意の組合せを使用することもできる。
【0107】
有利には、色の決定（上述のステップ３ｃ）には、眼球内のモデリングされた表面すなわち表面Ｃａ、Ｃｐ、ＬａおよびＬｐならびに網膜Ｒの位置と
とを比較するステップおよびこの決定に依存する関数として色をアトリビュートするステップが関与する。これにより、反射が由来している構造に応じて画素を色付けすることが可能になる。
【数10】
【0108】
より一般的に言うと、本方法によって計算される「データ」は、画像平面を定義する２次元画像を表わしていてよく、「組合せ値」は、画像内の画素にアトリビュートされる色度であってよい。該方法は、このとき、少なくとも以下のステップを含む：
－ （例えば点は、画像に直交して延在し画素に対応する横断面を有する空間体積内にあることから、任意には眼球の界面における屈折を考慮に入れて）、各画素について、その画素に対応する一点における少なくとも１つの反射値を識別するステップ。
－ 色が画像平面からの点の距離に依存するものとして、画素についての色を決定するステップ。
－ （例えば点は、画像に直交して延在し画素に対応する横断面を有する空間体積内にあることから、任意には眼球の界面における屈折を考慮に入れて）、各画素について、その画素に対応する一点における少なくとも１つの反射値を識別するステップ。
－ 色が画像平面からの点の距離に依存するものとして、画素についての色を決定するステップ。
【0109】
有利には、各画素について、画素に対応する最強の反射値を有する点が、色を決定するときに使用される。
【0110】
さらに、画素に対応する点の反射値の大きさを用いて、画素の輝度または明度を決定してもよい。
【0111】
画素の色を決定するステップには、有利には、以下のサブステップが関与する：
－ 前記距離を、反射値に対しモデリングされたいくつかの構造Ｃａ、Ｃｐ、Ｌａ、Ｌｐの位置と比較するステップ。これにより、構造のうちの１つの中または、構造の間、前または後ろの空間に対する「点」の「アトリビューション」を行なうことが可能になる。例えば、角膜の前の空間、角膜自体、角膜と水晶体の間の空間、水晶体の前面、水晶体の内部、水晶体の後面、水晶体と網膜の間の空間または網膜に対して、点をアトリビュートすることができる。例えば角膜をより細かく解像すること（角膜の前面、角膜の内部または角膜の後面に対して点をアトリビュートすること）によるか、または水晶体をより粗く解像すること（水晶体全体に対して点をアトリビュートすること）によるアトリビューションも同様に可能である。
－ このアトリビューションに応じて色を選ぶステップ。
－ 前記距離を、反射値に対しモデリングされたいくつかの構造Ｃａ、Ｃｐ、Ｌａ、Ｌｐの位置と比較するステップ。これにより、構造のうちの１つの中または、構造の間、前または後ろの空間に対する「点」の「アトリビューション」を行なうことが可能になる。例えば、角膜の前の空間、角膜自体、角膜と水晶体の間の空間、水晶体の前面、水晶体の内部、水晶体の後面、水晶体と網膜の間の空間または網膜に対して、点をアトリビュートすることができる。例えば角膜をより細かく解像すること（角膜の前面、角膜の内部または角膜の後面に対して点をアトリビュートすること）によるか、または水晶体をより粗く解像すること（水晶体全体に対して点をアトリビュートすること）によるアトリビューションも同様に可能である。
－ このアトリビューションに応じて色を選ぶステップ。
【0112】
これにより、例えば、所与の構造が画像平面との関係において湾曲している場合でも、この構造からの全ての反射に対して同じ色をアトリビュートすることが可能になる。
【0113】
註
本明細書中に記載されている技術は、眼球の視軸または光学軸に平行な横断面を表わすデータを決定するために特に好適である。この横断面は、１次元横断面、例えば視軸または光学軸に沿った線形横断面であってよく、または、例えば横断面画像として表わされた２次元横断面であってもよい。
本明細書中に記載されている技術は、眼球の視軸または光学軸に平行な横断面を表わすデータを決定するために特に好適である。この横断面は、１次元横断面、例えば視軸または光学軸に沿った線形横断面であってよく、または、例えば横断面画像として表わされた２次元横断面であってもよい。
【0114】
別の実施形態においては、データは、図示されているように、眼球の視軸または光学軸に直交する２次元横断面または画像を表わすものであってよい。
【0115】
しかしながら、データは例えば、眼球の他の部分、例えば眼球の単一点、眼球の体積全体、または一定のパラメータ（例えば眼球の長さなど・・・、以上の例を参照のこと）を表わすものでもあり得る。
【0116】
当該技術は、あらゆる種類のＯＣＴ、詳細にはタイムドメインＯＣＴならびに周波数ドメインＯＣＴのために使用可能である。しかしながら、Ａスキャンを迅速に得るその能力のため、周波数ドメインＯＣＴ、特に波長掃引型ＯＣＴが有利である。
【0117】
本発明の現在好まれている実施形態が図示され説明されてきたが、本発明はそれに限定されず、以下のクレームの範囲内で、他の形でさまざまに具体化および実践可能であるということをはっきりと理解すべきである。
上述の実施形態は下記の例のように記載され得るが、下記の例に限定されない。
［例１］
光干渉断層撮影を用いて、眼球内の構造を表わすデータ、詳細には横断面画像を表わすデータを得るための方法において、
－ Ａスキャン（Ａ１、Ａ２、Ａ３）のセットを記録するステップであって、前記Ａスキャン（Ａ１、Ａ２、Ａ３）の少なくともいくつかが眼球の異なる位置で記録され、各々のＡスキャン（Ａ１、Ａ２、Ａ３）が、眼球を通した光の軌跡に沿った複数の点について複数の反射値（ｒ）を生成するために使用される、ステップと、
－ 複数の組合せ値（ｒc）を計算するステップであって、各々の組合せ値（ｒc）が前記眼球内の異なる位置におけるいくつかの反射値（ｒ）から得られる、ステップと、
を含む方法。
［例２］
前記眼球の少なくとも１つの湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒ）を使用するステップを含む、例１に記載の方法。
［例３］
前記組合せ値（ｒc）を計算するために使用すべき前記Ａスキャン（Ａ１、Ａ２、Ａ３）のサブセット上の前記点を識別するために、前記モデリングされた湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）を使用するステップを含む、例２に記載の方法。
［例４］
前記湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）における組合せ値（ｒc）を計算するために、
－ 前記サブセット内の各々のＡスキャンｉについて、Ａスキャンｉとモデリングされた構造との交差部における前記反射値（ｒ）の推定値ｒe（ｉ）を決定するステップと、
－ 前記組合せ値（ｒc）を計算するために前記推定値ｒe（ｉ）を使用するステップと、
を含む、例２または３のいずれか１項に記載の方法。
［例５］
所与の深さ座標について、前記深さ座標に最も近い前記Ａスキャンｉ上の少なくとも２つの反射値（ｒ）を補間することによって前記サブグループ内のＡスキャンｉの推定値ｒe（ｉ）を計算するステップを含む、例３または４のいずれか１項に記載の方法。
［例６］
前記湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）が、角膜の前面（Ｃａ）、角膜の後面（Ｃｂ）、水晶体の前面（Ｌａ）、水晶体の後面（Ｌｐ）、および網膜（Ｒ）のうちの１つである、例２ないし５のいずれか１項に記載の方法。
［例７］
所与の軸に沿って前記組合せ値（ｒc）を計算するために、
－ 眼球内の少なくとも第１および第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒc）を使用するステップであって、前記第１の構造が前記軸に沿った第１の位置に位置設定されており、前記第２の構造が前記軸に沿った第２の位置に位置設定されている、ステップ
を含む、例２ないし６のいずれか１項に記載の方法。
［例８］
－ 前記第１の位置における組合せ値（ｒc）を得るために、前記第１の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）に沿って接線方向に、点の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの第１の複数の前記反射値（ｒ）を組み合わせるステップ、および／または
－ 前記第２の位置における組合せ値（ｒc）を得るために、接線方向に前記第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）に沿った点の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの第２の複数の前記反射値（ｒ）を組み合わせるステップ、
のうちの少なくとも１つをさらに含む、例７に記載の方法。
［例９］
－ 前記第１および前記第２の位置の間の第３の位置における組合せ値（ｒc）を得るために、異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの点の第３の複数の反射値（ｒ）を組み合わせるステップであって、前記点が前記第１および前記第２の構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の間および／または近傍の領域内にあるステップ、
をさらに含む、例７または８のいずれか１項に記載の方法。
［例１０］
前記第３の複数の前記反射値に対応する前記点が、
－ 前記第１および前記第２の構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）から補間された曲線（Ｃ３）に沿って接線方向に配設されているか、または、
－ 前記第１の構造に平行な第１の曲線（Ｃ１）に沿って接線方向に配設された点の第１のサブセット、および前記第２の構造に平行な第２の曲線（Ｃ２）に沿って接線方向に配設された点の第２のサブセットを含み、ここで前記第１および第２の曲線（Ｃ１、Ｃ２）が前記第３の位置で交差している、
例９に記載の方法。
［例１１］
前記データが、前記構造の断面ビューを示す２次元横断面を表わしており、前記断面ビュー内の前記湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の画像点を生成するために前記組合せ値（ｒc）を使用するステップを含む、例２ないし１０のいずれか１項に記載の方法。
［例１２］
前記組合せ値（ｒc）の少なくともいくつかについて、前記組合せ値（ｒc）が、横断面の平面の外側の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの複数の前記反射値（ｒ）を組み合わせることによって計算される、例１１に記載の方法。
［例１３］
複数の組合せ値（ｒc）を計算する前記ステップが、
－ 組み合わせるべき前記反射値（ｒ）の平均、詳細には加重平均を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の分位数、詳細には中央値を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の最小値または最大値を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の補間を計算するステップ、
のうちの少なくとも１つを含む、例１ないし１２のいずれか１項に記載の方法。
［例１４］
複数の組合せ値（ｒc）を計算する前記ステップが、前記組合せ値（ｒc）を計算するときに鏡面反射を示すいかなるＡスキャンの反射値も使用しないステップを含む、例１ないし１３のいずれか１項に記載の方法。
［例１５］
組合せ値（ｒc）を組み合わせるために使用される前記反射値（ｒ）の前記位置が、前記軸から１ｍｍ未満、詳細には０．５ｍｍ未満の距離（Ｒ）内にある、例１ないし１４のいずれか１項に記載の方法。
［例１６］
前記データが、画像平面を定義する２次元画像を表わし、前記組合せ値（ｒc）が前記画像内の画素にアトリビュートされた色度であり、
－ 各画素について、前記画素に対応する点における少なくとも１つの反射値（ｒ）を識別するステップと；
－ 前記画像平面からの前記点の前記距離に応じた色を、前記画素について決定するステップと；
を含む、例１ないし１５のいずれか１項に記載の方法。
［例１７］
各画素について、前記色を決定するときに、この画素に対応する最強の反射値（ｒ）を有する点が使用される、例１６に記載の方法。
［例１８］
－ 前記眼球のいくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒ）を使用するステップと、
－ 前記いくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の位置と、前記画像平面からの前記点の距離を比較して、前記構造のうちの１つに対する、または前記構造の間、前または後ろの空間に対する前記点のアトリビューションを行なうステップと、
－ 前記アトリビューションに応じて前記色を選択するステップと、
を含む例１６または１７のいずれか１項に記載の方法。
［例１９］
前記画素の輝度または明度を決定するために前記点の反射値（ｒ）の大きさを使用するステップを含む、例１ないし１８のいずれか１項に記載の方法。
［例２０］
前記データが、前記眼球の一次元横断面を表わすものである、例１ないし１９のいずれか１項に記載の方法。
［例２１］
前記データが、前記眼球の視軸または光学軸に平行な１次元または２次元横断面を表わすものである、例１ないし２０のいずれか１項に記載の方法。
［例２２］
前記データが、前記眼球の視軸または光学軸に直交する２次元横断面を表わすものである、例１ないし２０のいずれか１項に記載の方法。
［例２３］
前記光干渉断層撮影が、周波数ドメインＯＣＴ、詳細には波長掃引型ＯＣＴである、例１ないし２２のいずれか１項に記載の方法。
［例２４］
前記干渉計（１０～２６）が単一のＡスキャンについて少なくとも４０ｍｍの測定範囲を有する、例１ないし２３のいずれか１項に記載の方法。
［例２５］
前記干渉計（１０～２６）が、テレセントリックプローブビームを生成する、例１ないし２４のいずれか１項に記載の方法。
［例２６］
－ 光干渉断層撮影干渉計（１０～２６）と、
－ 例１ないし２５のいずれか１項に記載の方法を実施するために構造化され適応された制御ユニット（３２）と、
を含む眼科用デバイス。
上述の実施形態は下記の例のように記載され得るが、下記の例に限定されない。
［例１］
光干渉断層撮影を用いて、眼球内の構造を表わすデータ、詳細には横断面画像を表わすデータを得るための方法において、
－ Ａスキャン（Ａ１、Ａ２、Ａ３）のセットを記録するステップであって、前記Ａスキャン（Ａ１、Ａ２、Ａ３）の少なくともいくつかが眼球の異なる位置で記録され、各々のＡスキャン（Ａ１、Ａ２、Ａ３）が、眼球を通した光の軌跡に沿った複数の点について複数の反射値（ｒ）を生成するために使用される、ステップと、
－ 複数の組合せ値（ｒc）を計算するステップであって、各々の組合せ値（ｒc）が前記眼球内の異なる位置におけるいくつかの反射値（ｒ）から得られる、ステップと、
を含む方法。
［例２］
前記眼球の少なくとも１つの湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒ）を使用するステップを含む、例１に記載の方法。
［例３］
前記組合せ値（ｒc）を計算するために使用すべき前記Ａスキャン（Ａ１、Ａ２、Ａ３）のサブセット上の前記点を識別するために、前記モデリングされた湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）を使用するステップを含む、例２に記載の方法。
［例４］
前記湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）における組合せ値（ｒc）を計算するために、
－ 前記サブセット内の各々のＡスキャンｉについて、Ａスキャンｉとモデリングされた構造との交差部における前記反射値（ｒ）の推定値ｒe（ｉ）を決定するステップと、
－ 前記組合せ値（ｒc）を計算するために前記推定値ｒe（ｉ）を使用するステップと、
を含む、例２または３のいずれか１項に記載の方法。
［例５］
所与の深さ座標について、前記深さ座標に最も近い前記Ａスキャンｉ上の少なくとも２つの反射値（ｒ）を補間することによって前記サブグループ内のＡスキャンｉの推定値ｒe（ｉ）を計算するステップを含む、例３または４のいずれか１項に記載の方法。
［例６］
前記湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）が、角膜の前面（Ｃａ）、角膜の後面（Ｃｂ）、水晶体の前面（Ｌａ）、水晶体の後面（Ｌｐ）、および網膜（Ｒ）のうちの１つである、例２ないし５のいずれか１項に記載の方法。
［例７］
所与の軸に沿って前記組合せ値（ｒc）を計算するために、
－ 眼球内の少なくとも第１および第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒc）を使用するステップであって、前記第１の構造が前記軸に沿った第１の位置に位置設定されており、前記第２の構造が前記軸に沿った第２の位置に位置設定されている、ステップ
を含む、例２ないし６のいずれか１項に記載の方法。
［例８］
－ 前記第１の位置における組合せ値（ｒc）を得るために、前記第１の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）に沿って接線方向に、点の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの第１の複数の前記反射値（ｒ）を組み合わせるステップ、および／または
－ 前記第２の位置における組合せ値（ｒc）を得るために、接線方向に前記第２の湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）に沿った点の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの第２の複数の前記反射値（ｒ）を組み合わせるステップ、
のうちの少なくとも１つをさらに含む、例７に記載の方法。
［例９］
－ 前記第１および前記第２の位置の間の第３の位置における組合せ値（ｒc）を得るために、異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの点の第３の複数の反射値（ｒ）を組み合わせるステップであって、前記点が前記第１および前記第２の構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の間および／または近傍の領域内にあるステップ、
をさらに含む、例７または８のいずれか１項に記載の方法。
［例１０］
前記第３の複数の前記反射値に対応する前記点が、
－ 前記第１および前記第２の構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）から補間された曲線（Ｃ３）に沿って接線方向に配設されているか、または、
－ 前記第１の構造に平行な第１の曲線（Ｃ１）に沿って接線方向に配設された点の第１のサブセット、および前記第２の構造に平行な第２の曲線（Ｃ２）に沿って接線方向に配設された点の第２のサブセットを含み、ここで前記第１および第２の曲線（Ｃ１、Ｃ２）が前記第３の位置で交差している、
例９に記載の方法。
［例１１］
前記データが、前記構造の断面ビューを示す２次元横断面を表わしており、前記断面ビュー内の前記湾曲構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の画像点を生成するために前記組合せ値（ｒc）を使用するステップを含む、例２ないし１０のいずれか１項に記載の方法。
［例１２］
前記組合せ値（ｒc）の少なくともいくつかについて、前記組合せ値（ｒc）が、横断面の平面の外側の異なるＡスキャン（Ａ１、Ａ２、Ａ３）からの複数の前記反射値（ｒ）を組み合わせることによって計算される、例１１に記載の方法。
［例１３］
複数の組合せ値（ｒc）を計算する前記ステップが、
－ 組み合わせるべき前記反射値（ｒ）の平均、詳細には加重平均を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の分位数、詳細には中央値を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の最小値または最大値を計算するステップ、
－ 組み合わせるべき前記反射値（ｒ）の補間を計算するステップ、
のうちの少なくとも１つを含む、例１ないし１２のいずれか１項に記載の方法。
［例１４］
複数の組合せ値（ｒc）を計算する前記ステップが、前記組合せ値（ｒc）を計算するときに鏡面反射を示すいかなるＡスキャンの反射値も使用しないステップを含む、例１ないし１３のいずれか１項に記載の方法。
［例１５］
組合せ値（ｒc）を組み合わせるために使用される前記反射値（ｒ）の前記位置が、前記軸から１ｍｍ未満、詳細には０．５ｍｍ未満の距離（Ｒ）内にある、例１ないし１４のいずれか１項に記載の方法。
［例１６］
前記データが、画像平面を定義する２次元画像を表わし、前記組合せ値（ｒc）が前記画像内の画素にアトリビュートされた色度であり、
－ 各画素について、前記画素に対応する点における少なくとも１つの反射値（ｒ）を識別するステップと；
－ 前記画像平面からの前記点の前記距離に応じた色を、前記画素について決定するステップと；
を含む、例１ないし１５のいずれか１項に記載の方法。
［例１７］
各画素について、前記色を決定するときに、この画素に対応する最強の反射値（ｒ）を有する点が使用される、例１６に記載の方法。
［例１８］
－ 前記眼球のいくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒ）を使用するステップと、
－ 前記いくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の位置と、前記画像平面からの前記点の距離を比較して、前記構造のうちの１つに対する、または前記構造の間、前または後ろの空間に対する前記点のアトリビューションを行なうステップと、
－ 前記アトリビューションに応じて前記色を選択するステップと、
を含む例１６または１７のいずれか１項に記載の方法。
［例１９］
前記画素の輝度または明度を決定するために前記点の反射値（ｒ）の大きさを使用するステップを含む、例１ないし１８のいずれか１項に記載の方法。
［例２０］
前記データが、前記眼球の一次元横断面を表わすものである、例１ないし１９のいずれか１項に記載の方法。
［例２１］
前記データが、前記眼球の視軸または光学軸に平行な１次元または２次元横断面を表わすものである、例１ないし２０のいずれか１項に記載の方法。
［例２２］
前記データが、前記眼球の視軸または光学軸に直交する２次元横断面を表わすものである、例１ないし２０のいずれか１項に記載の方法。
［例２３］
前記光干渉断層撮影が、周波数ドメインＯＣＴ、詳細には波長掃引型ＯＣＴである、例１ないし２２のいずれか１項に記載の方法。
［例２４］
前記干渉計（１０～２６）が単一のＡスキャンについて少なくとも４０ｍｍの測定範囲を有する、例１ないし２３のいずれか１項に記載の方法。
［例２５］
前記干渉計（１０～２６）が、テレセントリックプローブビームを生成する、例１ないし２４のいずれか１項に記載の方法。
［例２６］
－ 光干渉断層撮影干渉計（１０～２６）と、
－ 例１ないし２５のいずれか１項に記載の方法を実施するために構造化され適応された制御ユニット（３２）と、
を含む眼科用デバイス。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【手続補正書】
【提出日】2024-05-01
【手続補正1】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項1】
光干渉断層撮影を用いて、眼球内の構造を表わすデータを表わすデータを得るための方法において、
－ Ａスキャン（Ａ１、Ａ２、Ａ３）のセットを記録するステップであって、前記Ａスキャン（Ａ１、Ａ２、Ａ３）の少なくともいくつかが眼球の異なる位置で記録され、各々のＡスキャン（Ａ１、Ａ２、Ａ３）が、眼球を通した光の軌跡に沿った複数の点について複数の反射値（ｒ）を生成するために使用される、ステップと、
－ 複数の組合せ値（ｒc）を計算するステップであって、各々の組合せ値（ｒc）が前記眼球内の異なる位置におけるいくつかの反射値（ｒ）から得られる、ステップと、
を含み、
前記データが、画像平面を定義する２次元画像を表わし、前記組合せ値（ｒ c ）が前記画像内の画素にアトリビュートされた色度であり、
－ 各画素について、前記画素に対応する点における少なくとも１つの反射値（ｒ）を識別するステップと；
－ 前記画像平面からの前記点の距離に応じた色を、前記画素について決定するステップと；
を含む、方法。
－ Ａスキャン（Ａ１、Ａ２、Ａ３）のセットを記録するステップであって、前記Ａスキャン（Ａ１、Ａ２、Ａ３）の少なくともいくつかが眼球の異なる位置で記録され、各々のＡスキャン（Ａ１、Ａ２、Ａ３）が、眼球を通した光の軌跡に沿った複数の点について複数の反射値（ｒ）を生成するために使用される、ステップと、
－ 複数の組合せ値（ｒc）を計算するステップであって、各々の組合せ値（ｒc）が前記眼球内の異なる位置におけるいくつかの反射値（ｒ）から得られる、ステップと、
を含み、
前記データが、画像平面を定義する２次元画像を表わし、前記組合せ値（ｒ c ）が前記画像内の画素にアトリビュートされた色度であり、
－ 各画素について、前記画素に対応する点における少なくとも１つの反射値（ｒ）を識別するステップと；
－ 前記画像平面からの前記点の距離に応じた色を、前記画素について決定するステップと；
を含む、方法。
【請求項2】
各画素について、前記色を決定するときに、この画素に対応する最強の反射値（ｒ）を有する点が使用される、請求項１に記載の方法。
【請求項3】
－ 前記眼球のいくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）をモデリングするために前記反射値（ｒ）を使用するステップと、
－ 前記いくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の位置と、前記画像平面からの前記点の距離を比較して、前記構造のうちの１つに対する、または前記構造の間、前または後ろの空間に対する前記点のアトリビューションを行なうステップと、
－ 前記アトリビューションに応じて前記色を選択するステップと、
を含む請求項１または２のいずれか１項に記載の方法。
－ 前記いくつかの構造（Ｃａ、Ｃｐ、Ｌａ、Ｌｐ）の位置と、前記画像平面からの前記点の距離を比較して、前記構造のうちの１つに対する、または前記構造の間、前または後ろの空間に対する前記点のアトリビューションを行なうステップと、
－ 前記アトリビューションに応じて前記色を選択するステップと、
を含む請求項１または２のいずれか１項に記載の方法。
【請求項4】
前記画素の輝度または明度を決定するために前記点の反射値（ｒ）の大きさを使用するステップを含む、請求項１ないし３のいずれか１項に記載の方法。
【請求項5】
前記光干渉断層撮影が、周波数ドメインＯＣＴである、請求項１ないし４のいずれか１項に記載の方法。
【請求項6】
－ 光干渉断層撮影干渉計（１０～２６）と、
－ 請求項１ないし５のいずれか１項に記載の方法のすべてのステップを実施するために構造化され適応された制御ユニット（３２）と、
を含む眼科用デバイス。
－ 請求項１ないし５のいずれか１項に記載の方法のすべてのステップを実施するために構造化され適応された制御ユニット（３２）と、
を含む眼科用デバイス。
【請求項7】
前記光干渉断層撮影干渉計（１０～２６）が単一のＡスキャンについて少なくとも４０ｍｍの測定範囲を有する、請求項６に記載の眼科用デバイス。
【請求項8】
前記光干渉断層撮影干渉計（１０～２６）が、テレセントリックプローブビームを生成する、請求項６または７に記載の眼科用デバイス。