(57)【要約】【課題】被検眼の組織の形状を高い再現性で高精度に特定するための新たな技術を提供する。【解決手段】眼科情報処理装置は、取得部と、組織特定部と、特定部とを含む。取得部は、被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された被検眼の断層像を取得する。組織特定部は、断層像に対してセグメンテーション処理を施すことにより被検眼の組織の形状を表す第１形状データを取得する。特定部は、第１形状データから被検眼に対する光学系の位置の変化に対してその変動が小さい低感度成分を特定し、特定された低感度成分に基づいて組織の形状を表す第２形状データを求める。【選択図】図１
2020110579.tif 000002
【特許請求の範囲】
【請求項１】
被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された前記被検眼の断層像を取得する取得部と、
前記断層像に対してセグメンテーション処理を施すことにより前記被検眼の組織の形状を表す第１形状データを取得する組織特定部と、
前記第１形状データから前記被検眼に対する前記光学系の位置の変化に対してその変動が小さい低感度成分を特定し、前記特定された低感度成分に基づいて前記組織の形状を表す第２形状データを求める特定部と、
を含む眼科情報処理装置。
【請求項２】
前記特定部は、前記被検眼に対する前記光学系のアライメント基準位置に対し前記被検眼に対する前記光学系の位置の対称成分を前記低感度成分として特定する
ことを特徴とする請求項１に記載の眼科情報処理装置。
【請求項３】
前記特定部は、前記第１形状データを多項式近似することにより得られた多項式のうち、前記光学系の位置の偶数次成分を前記低感度成分として特定する
ことを特徴とする請求項２に記載の眼科情報処理装置。
【請求項４】
前記低感度成分は、前記組織の曲率を表す成分を含む
ことを特徴とする請求項１〜請求項３のいずれか一項に記載の眼科情報処理装置。
【請求項５】
前記特定部は、前記第１形状データに対して円フィッティング処理を行うことにより得られた円の位置を表す成分とサイズを表す成分とを前記低感度成分として特定する
ことを特徴とする請求項１に記載の眼科情報処理装置。
【請求項６】
被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された前記被検眼の断層像を取得する取得部と、
前記断層像に対してセグメンテーション処理を施すことにより前記被検眼の組織の形状を表す第１形状データを取得する組織特定部と、
前記第１形状データから前記被検眼に対する前記光学系の位置の変化に対してその変動が大きい高感度成分を特定し、前記特定された高感度成分を前記第１形状データから除去することにより前記組織の形状を表す第２形状データを求める特定部と、
を含む眼科情報処理装置。
【請求項７】
前記特定部は、前記被検眼に対する前記光学系のアライメント基準位置に対し前記被検眼に対する前記光学系の位置の非対称成分を前記高感度成分として特定する
ことを特徴とする請求項６に記載の眼科情報処理装置。
【請求項８】
前記特定部は、前記第１形状データを多項式近似することにより得られた多項式のうち、前記光学系の位置の奇数次成分を前記高感度成分として特定する
ことを特徴とする請求項７に記載の眼科情報処理装置。
【請求項９】
前記高感度成分は、所定の基準方向に対する前記組織の傾きを表す成分を含む
ことを特徴とする請求項６〜請求項８のいずれか一項に記載の眼科情報処理装置。
【請求項１０】
前記組織は、眼底における所定の層領域を含む
ことを特徴とする請求項１〜請求項９のいずれか一項に記載の眼科情報処理装置。
【請求項１１】
前記被検眼を他覚的に測定することにより得られた屈折度数と、前記特定部により取得された前記第２形状データに基づく前記組織の形状に対応した前記被検眼の光学特性を表すパラメータとに基づいて、前記被検眼の中心窩を含む領域の周辺領域の屈折度数を算出する算出部を含む
ことを特徴とする請求項１０に記載の眼科情報処理装置。
【請求項１２】
前記組織の形状は、所定の基準方向に対する前記組織のチルト角度を含む
ことを特徴とする請求項１１に記載の眼科情報処理装置。
【請求項１３】
前記光学系と、
前記被検眼と前記光学系とを相対的に移動する移動機構と、
請求項１〜請求項１２のいずれか一項に記載の眼科情報処理装置と、
を含む眼科装置。
【請求項１４】
被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された前記被検眼の断層像を取得する取得ステップと、
前記断層像に対してセグメンテーション処理を施すことにより前記被検眼の組織の形状を表す第１形状データを取得する組織特定ステップと、
前記第１形状データから前記被検眼に対する前記光学系の位置の変化に対してその変動が小さい低感度成分を特定し、前記特定された低感度成分に基づいて前記組織の形状を表す第２形状データを求める特定ステップと、
を含む眼科情報処理方法。
【請求項１５】
被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された前記被検眼の断層像を取得する取得ステップと、
前記断層像に対してセグメンテーション処理を施すことにより前記被検眼の組織の形状を表す第１形状データを取得する組織特定ステップと、
前記第１形状データから前記被検眼に対する前記光学系の位置の変化に対してその変動が大きい高感度成分を特定し、前記特定された高感度成分を前記第１形状データから除去することにより前記組織の形状を表す第２形状データを求める特定ステップと、
を含む眼科情報処理方法。
【請求項１６】
前記被検眼を他覚的に測定することにより得られた屈折度数と、前記特定ステップにおいて取得された前記第２形状データに基づく前記組織の形状に対応した前記被検眼の光学特性を表すパラメータとに基づいて、前記被検眼の中心窩を含む領域の周辺領域の屈折度数を算出する算出ステップを含む
ことを特徴とする請求項１４又は請求項１５に記載の眼科情報処理方法。
【請求項１７】
コンピュータに、請求項１４〜請求項１６のいずれか一項に記載の眼科情報処理方法の各ステップを実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラムに関する。
【背景技術】
【０００２】
近年、近視が進行する原因の１つとして、周辺視野の焦点が網膜面よりも奥側（強膜側）に存在することに起因して網膜が奥側に伸びようとして近視が進行する可能性が報告されている（例えば、非特許文献１）。
【０００３】
このような近視の進行を抑制するために、周辺視野の屈折力を高くすることで、中心視野の焦点位置を手前側（角膜側）に移動させる眼鏡やコンタクトレンズが開発されている。更に、事前に測定された波面収差に基づいて行われるｗａｖｅｆｒｏｎｔ−ｇｕｉｄｅｄ  ＬＡＳＩＫのような屈折矯正手術も行われている。従って、このような高機能な屈折矯正において、周辺視野の屈折力を正確に測定することが求められる。
【０００４】
また、眼球の形状には、いくつかのタイプが確認されている（例えば、非特許文献２）。
【０００５】
このような眼球の形状には、近視に多い形状のタイプ等もあり、近視の進行に伴う形状変化を測定し、その測定結果を近視の進行の対処法にフィードバックすることが有効であると考えられる。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Ｅａｒｌ  Ｌ．  Ｓｍｉｔｈ  ｅｔ  ａｌ．，  “Ｒｅｌａｔｉｖｅ  ｐｅｒｉｐｈｅｒａｌ  ｈｙｐｅｒｏｐｔｉｃ  ｄｅｆｏｃｕｓ  ａｌｔｅｒｓ  ｃｅｎｔｒａｌ  ｒｅｆｒａｃｔｉｖｅ  ｄｅｖｅｌｏｐｍｅｎｔ  ｉｎ  ｉｎｆａｎｔ  ｍｏｎｋｅｙｓ”，  Ｖｉｓｉｏｎ  Ｒｅｓｅａｒｃｈ，  Ｓｅｐｔｅｍｂｅｒ  ２００９，  ４９（１９），  ｐｐ．２３８６−２３９２
【非特許文献２】Ｐａｖａｎ  Ｋ  Ｖｅｒｋｉｃｈａｒｌａ  ｅｔ  ａｌ．，  “Ｅｙｅ  ｓｈａｐｅ  ａｎｄ  ｒｅｔｉｎａｌ  ｓｈａｐｅ，  ａｎｄ  ｔｈｅｉｒ  ｒｅｌａｔｉｏｎ  ｔｏ  ｐｅｒｉｐｈｅｒａｌ  ｒｅｆｒａｃｔｉｏｎ”，  Ｏｐｔｈａｌｍｉｃ  ＆  Ｐｈｙｓｉｏｌｏｇｉｃａｌ  Ｏｐｔｉｃｓ，  ３２（２０１２），  ｐｐ．１８４−１９９
【発明の概要】
【発明が解決しようとする課題】
【０００７】
一般的な眼科装置では、測定光軸上に固視標が投影されるため、網膜の中心窩の近傍の屈折度数が測定される。この場合、眼底等における組織の形状（眼球の形状）を考慮して、中心窩の近傍の屈折度数から周辺視野の屈折度数を求めることが考えられる。
【０００８】
しかしながら、眼底等における組織の形状を測定することを目的として光コヒーレンストモグラフィを用いて被検眼の断層像を取得する場合、眼のモーションによるアライメントのずれ量に起因して断層像を高い再現性で取得することが困難である。
【０００９】
本発明は、このような事情を鑑みてなされたものであり、その目的は、被検眼の組織の形状を高い再現性で高精度に特定するための新たな技術を提供することにある。
【課題を解決するための手段】
【００１０】
いくつかの実施形態の第１態様は、被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された前記被検眼の断層像を取得する取得部と、前記断層像に対してセグメンテーション処理を施すことにより前記被検眼の組織の形状を表す第１形状データを取得する組織特定部と、前記第１形状データから前記被検眼に対する前記光学系の位置の変化に対してその変動が小さい低感度成分を特定し、前記特定された低感度成分に基づいて前記組織の形状を表す第２形状データを求める特定部と、を含む眼科情報処理装置である。
【００１１】
いくつかの実施形態の第２態様では、第１態様において、前記特定部は、前記被検眼に対する前記光学系のアライメント基準位置に対し前記被検眼に対する前記光学系の位置の対称成分を前記低感度成分として特定する。
【００１２】
いくつかの実施形態の第３態様では、第２態様において、前記特定部は、前記第１形状データを多項式近似することにより得られた多項式のうち、前記光学系の位置の偶数次成分を前記低感度成分として特定する。
【００１３】
いくつかの実施形態の第４態様では、第１態様〜第３態様のいずれかにおいて、前記低感度成分は、前記組織の曲率を表す成分を含む。
【００１４】
いくつかの実施形態の第５態様では、第１態様において、前記特定部は、前記第１形状データに対して円フィッティング処理を行うことにより得られた円の位置を表す成分とサイズを表す成分とを前記低感度成分として特定する。
【００１５】
いくつかの実施形態の第６態様は、被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された前記被検眼の断層像を取得する取得部と、前記断層像に対してセグメンテーション処理を施すことにより前記被検眼の組織の形状を表す第１形状データを取得する組織特定部と、前記第１形状データから前記被検眼に対する前記光学系の位置の変化に対してその変動が大きい高感度成分を特定し、前記特定された高感度成分を前記第１形状データから除去することにより前記組織の形状を表す第２形状データを求める特定部と、を含む眼科情報処理装置である。
【００１６】
いくつかの実施形態の第７態様では、第６態様において、前記特定部は、前記被検眼に対する前記光学系のアライメント基準位置に対し前記被検眼に対する前記光学系の位置の非対称成分を前記高感度成分として特定する。
【００１７】
いくつかの実施形態の第８態様では、第７態様において、前記特定部は、前記第１形状データを多項式近似することにより得られた多項式のうち、前記光学系の位置の奇数次成分を前記高感度成分として特定する。
【００１８】
いくつかの実施形態の第９態様では、第６態様〜第８態様のいずれかにおいて、前記高感度成分は、所定の基準方向に対する前記組織の傾きを表す成分を含む。
【００１９】
いくつかの実施形態の第１０態様では、第１態様〜第９態様のいずれかにおいて、前記組織は、眼底における所定の層領域を含む。
【００２０】
いくつかの実施形態の第１１態様は、第１０態様において、前記被検眼を他覚的に測定することにより得られた屈折度数と、前記特定部により取得された前記第２形状データに基づく前記組織の形状に対応した前記被検眼の光学特性を表すパラメータとに基づいて、前記被検眼の中心窩を含む領域の周辺領域の屈折度数を算出する算出部を含む。
【００２１】
いくつかの実施形態の第１２態様では、第１１態様において、前記組織の形状は、所定の基準方向に対する前記組織のチルト角度を含む。
【００２２】
いくつかの実施形態の第１３態様は、前記光学系と、前記被検眼と前記光学系とを相対的に移動する移動機構と、第１態様〜第１２態様のいずれかに記載の眼科情報処理装置と、を含む眼科装置である。
【００２３】
いくつかの実施形態の第１４態様は、被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された前記被検眼の断層像を取得する取得ステップと、前記断層像に対してセグメンテーション処理を施すことにより前記被検眼の組織の形状を表す第１形状データを取得する組織特定ステップと、前記第１形状データから前記被検眼に対する前記光学系の位置の変化に対してその変動が小さい低感度成分を特定し、前記特定された低感度成分に基づいて前記組織の形状を表す第２形状データを求める特定ステップと、を含む眼科情報処理方法である。
【００２４】
いくつかの実施形態の第１５態様は、被検眼に対して光コヒーレンストモグラフィを実行するための光学系を用いて得られたスキャンデータに基づいて形成された前記被検眼の断層像を取得する取得ステップと、前記断層像に対してセグメンテーション処理を施すことにより前記被検眼の組織の形状を表す第１形状データを取得する組織特定ステップと、前記第１形状データから前記被検眼に対する前記光学系の位置の変化に対してその変動が大きい高感度成分を特定し、前記特定された高感度成分を前記第１形状データから除去することにより前記組織の形状を表す第２形状データを求める特定ステップと、を含む眼科情報処理方法である。
【００２５】
いくつかの実施形態の第１６態様は、第１４態様又は第１５態様において、前記被検眼を他覚的に測定することにより得られた屈折度数と、前記特定ステップにおいて取得された前記第２形状データに基づく前記組織の形状に対応した前記被検眼の光学特性を表すパラメータとに基づいて、前記被検眼の中心窩を含む領域の周辺領域の屈折度数を算出する算出ステップを含む。
【００２６】
いくつかの実施形態の第１７態様は、コンピュータに、第１４態様〜第１６態様のいずれかに記載の眼科情報処理方法の各ステップを実行させるプログラムである。
【００２７】
なお、上記した複数の態様に係る構成を任意に組み合わせることが可能である。
【発明の効果】
【００２８】
本発明によれば、被検眼の組織の形状を高い再現性で高精度に特定するための新たな技術を提供することができる。
【図面の簡単な説明】
【００２９】
【図１】実施形態に係る眼科装置の構成の一例を示す概略図である。
【図２】実施形態に係る眼科装置の構成の一例を示す概略図である。
【図３】実施形態に係る眼科装置の構成の一例を示す概略図である。
【図４】実施形態に係る眼科装置の構成の一例を示す概略図である。
【図５】実施形態に係る眼科装置の構成の一例を示す概略図である。
【図６】実施形態に係る眼科装置の動作を説明するための概略図である。
【図７】実施形態に係る眼科装置の動作を説明するための概略図である。
【図８】実施形態に係る眼科装置の動作の一例を示すフローチャートである。
【図９】実施形態に係る眼科装置の動作の一例を示すフローチャートである。
【図１０】実施形態に係る眼科装置の動作を説明するための概略図である。
【図１１】実施形態の変形例に係る眼科装置の動作を説明するための概略図である。
【図１２】実施形態の変形例に係る眼科装置の動作の一例を示すフローチャートである。
【図１３】実施形態の変形例に係る眼科装置の動作の一例を示すフローチャートである。
【発明を実施するための形態】
【００３０】
この発明に係る眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラムの実施形態の例について、図面を参照しながら詳細に説明する。なお、この明細書において引用された文献の記載内容や任意の公知技術を、以下の実施形態に援用することが可能である。
【００３１】
実施形態に係る眼科情報処理装置は、被検眼の断層像を取得する。断層像は、例えば、外部の眼科装置を用いて被検眼に対して光コヒーレンストモグラフィ（Ｏｐｔｉｃａｌ  Ｃｏｈｅｒｅｎｃｅ  Ｔｏｍｏｇｒａｐｈｙ：ＯＣＴ）を実行することにより得られる。眼科情報処理装置は、取得された断層像を解析して所定の組織の形状データ（第１形状データ）を求め、求められた形状データから組織の真の形状データ（第２形状データ）を求める。求められた形状データから組織の真の形状を特定（推定、予測）することが可能である。
【００３２】
断層像を解析して得られた形状データには、被検眼に対するＯＣＴを実行するための光学系の位置の変化（微小変化）に対して、その変動が大きい高感度成分とその変動が小さい低感度成分とが含まれる。いくつかの実施形態では、眼科情報処理装置は、断層像を解析して得られた形状データから低感度成分を特定し、特定された低感度成分に基づいて組織の真の形状データを特定する。いくつかの実施形態では、眼科情報処理装置は、断層像を解析して得られた形状データから高感度成分を特定し、形状データから高感度成分を除去することにより組織の真の形状データを特定する。
【００３３】
いくつかの実施形態では、断層像を解析することにより得られた所定の層領域の１次元、２次元、又は３次元の形状データにより層領域の形状が特定される。いくつかの実施形態に係る眼科装置は、上記の眼科情報処理装置を含み、眼科情報処理装置の機能を実現する。
【００３４】
上記のように断層像又は形状データから被検眼の所定の組織の形状を特定することで、被検眼に対するＯＣＴを実行するための光学系のアライメントのずれ量（アライメントエラー量）の影響を低減し、被検眼の組織の形状を高い再現性で高精度に特定することが可能になる。
【００３５】
被検眼における組織の形状として、前眼部における組織の形状、後眼部における組織の形状などがある。前眼部における組織の形状として、角膜、虹彩、水晶体、毛様体、チン小帯、又は隅角の形状がある。後眼部における組織の形状として、眼底（眼底における所定の層領域）の形状などがある。以下、実施形態に係る組織の形状として、眼底における所定の層領域の形状を例に説明する。いくつかの実施形態では、眼底における層領域の形状を眼底の形状として特定することが可能である。眼底における層領域には、内境界膜、神経繊維層、神経節細胞層、内網状層、内顆粒層、外網状層、外顆粒層、外境界膜、視細胞層、網膜色素上皮層、脈絡膜、強膜、各層領域の界面などがある。しかしながら、眼底以外の眼球の任意の部位の形状について、以下の実施形態を適用することが可能である。また、以下の実施形態において、眼底の形状を表す形状データを、形状プロファイルと表記する場合がある。形状プロファイルは、所定の１次元の方向、２次元の方向、又は３次元の方向における形状の変化を表すデータである。
【００３６】
いくつかの実施形態に係る眼科情報処理装置は、上記のように特定された眼底の形状を用いて、眼底における中心窩を含む領域の外側の周辺領域の屈折度数を算出（推定）する。例えば、眼科情報処理装置は、被検眼の中心窩を含む領域の屈折度数と特定された眼底の形状とに基づいて、中心窩を含む領域の外側の周辺領域の屈折度数を算出する。
【００３７】
実施形態に係る眼科情報処理装置は、公知の模型眼等の眼球モデルのパラメータ（眼球の光学特性を表すパラメータ）を用いて上記の領域の屈折度数を算出することが可能である。パラメータには、眼軸長データ、前房深度データ、水晶体の形状を表す水晶体形状データ（水晶体曲率、水晶体厚など）、角膜の形状を表す角膜形状データ（角膜曲率半径、角膜厚など）などがある。眼科情報処理装置は、眼球モデルのパラメータの一部を被検眼の実測値に置き換えて新たな眼球モデルを構築し、構築された新たな眼球モデルを用いて上記の領域の屈折度数を算出することが可能である。いくつかの実施形態では、上記のパラメータは、電子カルテシステム、医用画像アーカイビングシステム、又は外部装置等から取得される。
【００３８】
実施形態に係る眼科情報処理方法は、実施形態に係る眼科情報処理装置においてプロセッサ（コンピュータ）により実行される処理を実現するための１以上のステップを含む。実施形態に係るプログラムは、プロセッサに実施形態に係る眼科情報処理方法の各ステップを実行させる。
【００３９】
本明細書において「プロセッサ」は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、プログラマブル論理デバイス（例えば、ＳＰＬＤ（Ｓｉｍｐｌｅ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｌｏｇｉｃ  Ｄｅｖｉｃｅ）、ＣＰＬＤ（Ｃｏｍｐｌｅｘ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｌｏｇｉｃ  Ｄｅｖｉｃｅ）、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ））等の回路を意味する。プロセッサは、例えば、記憶回路や記憶装置に格納されているプログラムを読み出し実行することで、実施形態に係る機能を実現する。
【００４０】
この明細書では、ＯＣＴによって取得される画像をＯＣＴ画像と総称することがある。また、ＯＣＴ画像を形成するための計測動作をＯＣＴ計測と呼び、ＯＣＴ計測を行うためのスキャンをＯＣＴスキャンと呼び、ＯＣＴスキャンにより得られたデータをスキャンデータと呼ぶことがある。
【００４１】
以下では、実施形態に係る眼科装置が、眼科情報処理装置を含む場合について説明する。眼科装置は、被検眼に対してＯＣＴを実行することにより被検眼の断層像を取得する。しかしながら、実施形態に係る眼科情報処置装置が、外部の眼科装置からスキャンデータ（ＯＣＴデータ）、断層像、後述の形状プロファイル等を取得するように構成されていてもよい。
【００４２】
いくつかの実施形態に係る眼科装置は、ＯＣＴ装置を含み、被検眼とＯＣＴを実行するための光学系との位置合わせを行うように構成される。いくつかの実施形態に係る眼科装置は、更に、他覚屈折測定装置を含む。いくつかの実施形態に係る眼科装置は、外部装置や記録媒体からデータを受信するデバイス（通信インターフェイス、入出力インターフェイス等）を含む。
【００４３】
すなわち、実施形態に係る眼科装置は、例えば、次のいずれかであってよい：（Ａ）他覚屈折測定装置（屈折測定部）とＯＣＴ装置（ＯＣＴ部）とを含む検査装置：（Ｂ）他覚屈折測定装置（屈折測定部）を含まず、ＯＣＴ装置（ＯＣＴ部）を含む検査装置；（Ｃ）他覚屈折測定装置（屈折測定部）及びＯＣＴ装置（ＯＣＴ部）のいずれも含まない情報処理装置。
【００４４】
以下では、眼科装置の装置光学系の光軸（測定光軸、検査光軸）に直交する左右方向（水平方向）をｘ方向とし、当該光軸に直交する上下方向（垂直方向）をｙ方向とし、当該光軸の方向（奥行き方向、前後方向）をｚ方向とする。
【００４５】
＜構成＞
図１〜図５に、実施形態に係る眼科装置の構成例を示す。図１は、実施形態に係る眼科装置１の構成例を表す機能ブロック図である。図２は、図１のデータ処理部７０の構成例を表す機能ブロック図である。図３は、図２のアライメント処理部７１の構成例を表す機能ブロック図である。図４は、図２の特定部７２２の構成例を表す機能ブロック図である。図５は、図２の算出部７３の構成例を表す機能ブロック図である。
【００４６】
眼科装置１は、他覚屈折測定装置（屈折測定部）とＯＣＴ装置（ＯＣＴ部）とを含む検査装置である。眼科装置１は、測定部１０と、制御処理部５０と、移動機構９０と、撮影部１００とを含む。測定部１０は、屈折測定部２０と、ＯＣＴ部３０と、アライメント光投射部４０と、ビームスプリッタＢＳ１、ＢＳ２とを含む。制御処理部５０は、画像形成部６０と、データ処理部７０と、制御部８０とを含む。
【００４７】
（屈折測定部２０）
屈折測定部２０は、制御部８０からの制御を受け、被検眼Ｅの屈折度数を他覚的に測定する。屈折測定部２０は、他覚屈折測定を行うための１以上の光学部材が設けられた光学系を含む。屈折測定部２０は、例えば、公知のレフラクトメータと同様の構成を有する。図示は省略するが、典型的なレフラクトメータは、特開２０１６−０７７７７４号公報に開示されているように、投影系と、受光系とを含む。
【００４８】
屈折測定部２０の投影系は、光源から出射した光を被検眼Ｅの眼底Ｅｆに投影する。投影系は、例えば、光源からの光を、コリメートレンズ、合焦レンズ、リレーレンズ、瞳レンズ、穴開きプリズム、偏心プリズム、対物レンズ等を通じて眼底Ｅｆに投影する。
【００４９】
屈折測定部２０の受光系は、眼底Ｅｆからの反射光を、対物レンズ、偏心プリズム、穴開きプリズム、他の瞳レンズ、他のリレーレンズ、他の合焦レンズ、円錐プリズム、結像レンズ等を通じて、撮像素子に投影する。これにより、撮像素子の撮像面に結像したリングパターン像が検出される。
【００５０】
いくつかの実施形態では、屈折測定部２０は、リング状の光を眼底Ｅｆに投影し、眼底Ｅｆからの反射光により形成されるリングパターン像を検出するように構成される。いくつかの実施形態では、屈折測定部２０は、輝点を眼底Ｅｆに投影し、眼底Ｅｆからの反射光をリング状の光に変換し、変換されたリング状の光により形成されるリングパターン像を検出するように構成される。
【００５１】
（ＯＣＴ部３０）
ＯＣＴ部３０は、制御部８０からの制御を受け、被検眼ＥにＯＣＴスキャンを適用してＯＣＴデータ（スキャンデータ）を取得する。ＯＣＴデータは、干渉信号データでもよいし、干渉信号データにフーリエ変換を適用して得られた反射強度プロファイルデータでもよいし、反射強度プロファイルデータを画像化して得られた画像データでもよい。
【００５２】
ＯＣＴ部３０が実施可能なＯＣＴ手法は、典型的にはフーリエドメインＯＣＴであり、スペクトラルドメインＯＣＴ及びスウェプトソースＯＣＴのいずれでもよい。スウェプトソースＯＣＴは、波長可変光源からの光を測定光と参照光とに分割し、被検眼に投射された測定光の戻り光を参照光と重ね合わせて干渉光を生成し、この干渉光を光検出器で検出し、波長の掃引及び測定光のスキャンに応じて収集された検出データ（干渉信号データ）にフーリエ変換等を施して反射強度プロファイルデータを形成する手法である。一方、スペクトラルドメインＯＣＴは、低コヒーレンス光源（広帯域光源）からの光を測定光と参照光とに分割し、被検眼に投射された測定光の戻り光を参照光と重ね合わせて干渉光を生成し、この干渉光のスペクトル分布を分光器で検出し、分光器による検出データ（干渉信号データ）にフーリエ変換等を施して反射強度プロファイルデータを形成する手法である。すなわち、スウェプトソースＯＣＴはスペクトル分布を時分割で取得するＯＣＴ手法であり、スペクトラルドメインＯＣＴはスペクトル分布を空間分割で取得するＯＣＴ手法である。
【００５３】
ＯＣＴ部３０は、ＯＣＴ計測を行うための１以上の光学部材が設けられた光学系を含む。ＯＣＴ部３０は、例えば、公知のＯＣＴ装置と同様の構成を有する。図示は省略するが、典型的なＯＣＴ装置は、特開２０１６−０７７７７４号公報に開示されているように、光源と、干渉光学系と、スキャン系と、検出系とを含む。
【００５４】
光源から出力された光は、干渉光学系によって測定光と参照光とに分割される。参照光は、参照アームに導かれる。測定光は、測定アームを通じて被検眼Ｅ（例えば、眼底Ｅｆ）に投射される。測定アームにはスキャン系が設けられている。スキャン系は、例えば光スキャナーを含み、測定光を１次元的又は２次元的に偏向可能である。光スキャナーは、１以上のガルバノスキャナを含む。スキャン系は、所定のスキャンモードにしたがって測定光を偏向する。
【００５５】
後述の制御部８０は、スキャンモードにしたがってスキャン系を制御することが可能である。スキャンモードには、ラインスキャン、ラスタースキャン（３次元スキャン）、サークルスキャン、同心円スキャン、ラジアルスキャン、クロススキャン、マルチクロススキャン、スパイラルスキャンなどがある。ラインスキャンは、直線状の軌跡に沿ったスキャンパターンである。ラスタースキャンは、互いに平行に配列された複数のラインスキャンからなるスキャンパターンである。サークルスキャンは、円形状の軌跡に沿ったスキャンパターンである。同心円スキャンは、同心状に配列された複数のサークルスキャンからなるスキャンパターンである。ラジアルスキャンは、放射状に配列された複数のラインスキャンからなるスキャンパターンである。クロススキャンは、互いに直交に配列された２つのラインスキャンからなるスキャンパターンである。マルチクロススキャンは、互いに直交する２つのラインスキャン群（例えば、各群は、互いに平行な５本のラインを含む）からなるスキャンパターンである。スパイラルスキャンは、中心から螺旋状に伸びるスキャンパターンである。
【００５６】
眼底Ｅｆに投射された測定光は、眼底Ｅｆの様々な深さ位置（層境界等）において反射・散乱される。被検眼Ｅからの測定光の戻り光は、干渉光学系によって参照光と合成される。測定光の戻り光と参照光とは重ね合わせの原理にしたがって干渉光を生成する。この干渉光は検出系によって検出される。検出系は、典型的には、スペクトラルドメインＯＣＴでは分光器を含み、スウェプトソースＯＣＴではバランスドフォトダイオード及びデータ収集システム（ＤＡＱ）を含む。
【００５７】
（アライメント光投射部４０）
アライメント光投射部４０は、被検眼Ｅと測定部１０（ＯＣＴ部、装置光学系）との位置合わせを行うためのアライメント光を被検眼Ｅに投射する。アライメント光投射部４０は、アライメント光源と、コリメータレンズとを含む。アライメント光投射部４０の光路は、ビームスプリッタＢＳ２により屈折測定部２０の光路に結合される。アライメント光源から出力された光は、コリメータレンズを経由し、ビームスプリッタＢＳ２により反射され、屈折測定部２０の光路を通じて被検眼Ｅに投射される。
【００５８】
いくつかの実施形態では、特開２０１６−０７７７７４号公報に開示されているように、被検眼Ｅの角膜Ｅｃ（前眼部）による反射光は、屈折測定部２０の光路を通じて、屈折測定部２０の受光系に導かれる。
【００５９】
被検眼Ｅの角膜Ｅｃによる反射光に基づく像（輝点像）は、撮影部１００により取得された前眼部像に含まれる。例えば、制御処理部５０は、輝点像を含む前眼部像とアライメントマークとを表示部（不図示）の表示画面に表示させる。手動でＸＹアライメント（上下左右方向のアライメント）を行う場合、ユーザは、アライメントマーク内に輝点像を誘導するように光学系の移動操作を行うことができる。手動でＺアライメント（前後方向のアライメント）を行う場合、ユーザは、表示部の表示画面に表示された前眼部像を参照しながら光学系の移動操作を行うことができる。自動でアライメントを行う場合、制御部８０は、アライメントマークと輝点像の位置との変位がキャンセルされるように移動機構９０を制御することにより、被検眼Ｅに対して測定部１０（光学系）を相対移動させる。また、制御部８０は、被検眼Ｅの所定部位（例えば、瞳孔中心位置）の位置と輝点像の位置とに基づいて、所定のアライメント完了条件を満たすように移動機構９０を制御することにより、被検眼Ｅに対して測定部１０（光学系）を相対移動させることが可能である。
【００６０】
（ビームスプリッタＢＳ１）
ビームスプリッタＢＳ１は、屈折測定部２０の光学系（投影系及び受光系）の光路に、ＯＣＴ部３０の光学系（干渉光学系など）の光路を同軸に結合する。例えば、ビームスプリッタＢＳ１としてダイクロイックミラーが用いられる。
【００６１】
（ビームスプリッタＢＳ２）
ビームスプリッタＢＳ２は、屈折測定部２０の光学系（投影系及び受光系）の光路に、アライメント光投射部４０の光学系の光路を同軸に結合する。例えば、ビームスプリッタＢＳ２としてハーフミラーが用いられる。
【００６２】
いくつかの実施形態では、眼科装置１は、制御部８０からの制御を受け、被検眼Ｅの視線を誘導するための固視標を被検眼Ｅに提示する機能（固視投影系）を有する。固視標は、被検眼Ｅに提示される内部固視標でもよいし、僚眼に提示される外部固視標でもよい。いくつかの実施形態では、ＯＣＴ部３０とビームスプリッタＢＳ１との間に配置された光路結合部材（例えば、ビームスプリッタ）により、固視投影系の光路がＯＣＴ部３０の干渉光学系の光路に同軸に結合されるように構成される。
【００６３】
制御部８０からの制御を受け、固視投影系による眼底Ｅｆにおける固視標の投影位置は変更可能である。いくつかの実施形態では、固視標は、同軸に結合された屈折測定部２０の光学系、及びＯＣＴ部３０の光学系の測定光軸上に投影される。いくつかの実施形態では、固視標は、眼底Ｅｆにおいて測定光軸から外れた位置に投影される。
【００６４】
（撮影部１００）
撮影部１００は、被検眼Ｅの前眼部を撮影するための１以上の前眼部カメラを含む。撮影部１００は、被検眼Ｅの正面画像である前眼部像を取得する。いくつかの実施形態では、１以上の前眼部カメラの近傍に少なくとも１つの前眼部照明光源（赤外光源等）が設けられる。例えば、各前眼部カメラについて、その上方近傍と下方近傍のそれぞれに前眼部照明光源が設けられる。
【００６５】
眼科装置１は、撮影部１００により取得された正面画像を用いて、被検眼Ｅに対する測定部１０（光学系）の位置合わせを行うことが可能である。いくつかの実施形態では、眼科装置１は、被検眼Ｅの前眼部を撮影することにより得られた正面画像を解析することにより被検眼Ｅの３次元位置を特定し、特定された３次元位置に基づいて測定部１０を相対移動することで位置合わせを行う。いくつかの実施形態では、眼科装置１は、被検眼Ｅの特徴位置とアライメント光投射部４０により投射されたアライメント光により形成された像の位置との変位がキャンセルされるように位置合わせを行う。
【００６６】
上記のように、撮影部１００は、１以上の前眼部カメラを含む。撮影部１００が単一の前眼部カメラを含む場合、眼科装置１は、取得された正面画像を解析して、測定部１０の光軸に直交する平面（水平方向（Ｘ方向）と垂直方向（Ｙ方向）とにより規定される平面）おける被検眼Ｅの２次元位置を特定する。この場合、眼科装置１には、測定部１０の光軸の方向における被検眼Ｅの位置を特定するためのアライメント光学系が設けられている。このようなアライメント光学系として、例えば、特開２０１６−０７７７７４号公報に開示されている光テコ方式の光学系などがある。眼科装置１は、このようなアライメント光学系を用いて測定部１０の（測定）光軸の方向における被検眼Ｅの位置と、上記の２次元位置とから、被検眼Ｅの３次元位置を特定することが可能である。
【００６７】
撮影部１００が２以上の前眼部カメラを含む場合、２以上の前眼部カメラは、被検眼Ｅの前眼部を異なる方向から撮影する。２以上の前眼部カメラは、異なる２以上の方向から実質的に同時に前眼部を撮影することができる。「実質的に同時」とは、例えば、２以上の前眼部カメラによる撮影において、眼球運動を無視できる程度の撮影タイミングのズレを許容することを示す。それにより、被検眼Ｅが実質的に同じ位置（向き）にあるときの画像を２以上の前眼部カメラによって取得することができる。眼科装置１は、例えば、特開２０１３−２４８３７６号公報に開示されているように、取得された正面画像を解析して被検眼Ｅの特徴位置を特定し、２以上の前眼部カメラの位置と特定された被検眼Ｅの特徴位置とから被検眼Ｅの３次元位置を特定する。
【００６８】
２以上の前眼部カメラによる撮影は動画撮影でも静止画撮影でもよい。動画撮影の場合、撮影開始タイミングを合わせるよう制御したり、フレームレートや各フレームの撮影タイミングを制御したりすることにより、上記のような実質的に同時の前眼部撮影を実現することができる。一方、静止画撮影の場合、撮影タイミングを合わせるよう制御することにより、これを実現することができる。
【００６９】
以下では、撮影部１００が２台の前眼部カメラを含むものとする。２台の前眼部カメラのそれぞれは、例えば特開２０１３−２４８３７６号公報に開示されているように、測定光軸（ＯＣＴ部３０の光軸）から外れた位置に設けられている。いくつかの実施形態では、２台の前眼部カメラの一方は、屈折測定部２０の受光系における撮像素子である。
【００７０】
（制御処理部５０）
制御処理部５０は、眼科装置１を動作させるための各種演算や各種制御を実行する。制御処理部５０は、１以上のプロセッサと、１以上の記憶装置とを含む。記憶装置としては、ランダムアクセスメモリ（ＲＡＭ）、リードオンリーメモリ（ＲＯＭ）、ハードディスクドライブ（ＨＤＤ）、ソリッドステートドライブ（ＳＳＤ）などがある。記憶装置には各種のコンピュータプログラムが格納されており、それに基づきプロセッサが動作することによって本例に係る演算や制御が実現される。
【００７１】
（画像形成部６０）
画像形成部６０は、被検眼Ｅに対してＯＣＴを実行することにより得られたスキャンデータに基づいて、被検眼Ｅの画像（断層像等）を形成する。画像形成部６０は、ＯＣＴ部３０の検出系による検出データに基づいてＯＣＴデータ（典型的には画像データ）を構築する。画像形成部６０は、従来のＯＣＴデータ処理と同様に、フィルター処理、高速フーリエ変換（ＦＦＴ）などを検出データに適用することにより、各Ａライン（被検眼Ｅ内における測定光の経路）における反射強度プロファイルデータを構築する。更に、画像形成部６０は、この反射強度プロファイルデータに画像化処理（画像表現）を適用することにより、各Ａラインの画像データ（Ａスキャンデータ）を構築する。いくつかの実施形態では、画像形成部６０の機能はプロセッサにより実現される。
【００７２】
いくつかの実施形態では、画像形成部６０の機能の少なくとも一部がＯＣＴ部３０に設けられる。
【００７３】
（データ処理部７０）
データ処理部７０は、各種のデータ処理を実行する。データ処理部７０は、スキャン系によるスキャンモードにしたがって複数のＡスキャンデータを配列することによりＢスキャンデータを構築することができる。いくつかの実施形態では、データ処理部７０は、２以上のＢスキャンデータの重ね合わせ処理を行う。データ処理部７０は、スキャン系によるスキャンモードにしたがって複数のＢスキャンデータを配列することによりスタックデータを構築することができる。データ処理部７０は、スタックデータからボリュームデータ（ボクセルデータ）を構築することができる。データ処理部７０は、スタックデータ又はボリュームデータをレンダリングすることができる。レンダリング手法としては、ボリュームレンダリング、多断面再構成（ＭＰＲ）、サーフェスレンダリング、プロジェクションなどがある。
【００７４】
データ処理部７０は、被検眼Ｅに対する測定部１０の位置合わせを行うための処理を実行することが可能である。位置合わせを行うための処理として、撮影部１００を用いて取得された被検眼Ｅの正面画像の解析処理、被検眼Ｅの位置の算出処理、被検眼Ｅに対する測定部１０の変位の算出処理などがある。
【００７５】
また、データ処理部７０は、位置合わせ後にＯＣＴ計測を実行して得られた被検眼Ｅの断層像から被検眼Ｅの眼底Ｅｆの形状を表す形状データ（形状プロファイル）を取得し、取得された形状データから眼底Ｅｆの真の形状を表す形状データを生成することが可能である。例えば、眼底Ｅｆの真の形状を表す形状データは、取得された形状データから、アライメントエラーに対する低感度成分を特定することで得られる。更に、データ処理部７０は、特定された眼底Ｅｆの形状（眼底Ｅｆの真の形状を表す形状データ）を用いて、被検眼Ｅの中心窩を含む領域の周辺領域の屈折度数を算出することが可能である。
【００７６】
このようなデータ処理部７０は、図２に示すように、アライメント処理部７１と、解析部７２と、算出部７３とを含む。
【００７７】
（アライメント処理部７１）
アライメント処理部７１は、上記のように、被検眼Ｅに対する測定部１０の位置合わせを行うための処理を実行する。いくつかの実施形態では、アライメント処理部７１は、前眼部カメラにより得られた撮影画像の歪みを補正し、歪みが補正された撮影画像を用いて上記の位置合わせを行うための処理を実行する。この場合、アライメント処理部７１は、制御処理部５０又はデータ処理部７０に設けられた記憶部にあらかじめ記憶されている収差情報に基づいて撮影画像の歪みを補正する。この処理は、例えば、歪曲収差を補正するための補正係数に基づく公知の画像処理技術によって実行される。
【００７８】
アライメント処理部７１は、図３に示すように、プルキンエ像特定部７１Ａと、瞳孔中心特定部７１Ｂと、移動目標位置決定部７１Ｃとを含む。
【００７９】
（プルキンエ像特定部７１Ａ）
アライメント光投射部４０により被検眼Ｅの前眼部にアライメント光を投射することにより、プルキンエ像が形成される。プルキンエ像は、角膜曲率半径の２分の１の距離だけ角膜頂点から測定部１０の光軸方向（ｚ方向）に変位した位置に形成される。
【００８０】
アライメント光が投射されている前眼部は、２つの前眼部カメラによって実質的に同時に撮影される。プルキンエ像特定部７１Ａは、２つの前眼部カメラにより実質的に同時に取得された２つの撮影画像のそれぞれを解析することでプルキンエ像（プルキンエ像に相当する画像領域）を特定する。この特定処理は、例えば従来と同様に、プルキンエ像に相当する輝点（高輝度の画素）を探索するための、画素値に関する閾値処理を含む。それにより、プルキンエ像に相当する撮影画像中の画像領域が特定される。
【００８１】
プルキンエ像特定部７１Ａは、プルキンエ像に相当する画像領域における代表点の位置を求めることができる。代表点は、例えば、当該画像領域の中心点又は重心点であってよい。この場合、プルキンエ像特定部７１Ａは、例えば、当該画像領域の周縁の近似円又は近似楕円を求め、近似円又は近似楕円の中心又は重心を求めることができる。
【００８２】
２つの撮影画像（前眼部像）のそれぞれは、前眼部を斜め方向から撮影して得られた画像である。各撮影画像には、瞳孔領域とプルキンエ像とが描出されている。プルキンエ像特定部７１Ａは、２つの撮影画像内のプルキンエ像を特定する。
【００８３】
ここで、２つの撮影画像は、測定部１０（対物レンズ）の光軸と異なる方向からの撮影により取得された画像である。ＸＹアライメントが実質的に合っているとき、２つの撮影画像内のプルキンエ像は測定部１０の光軸上に形成される。
【００８４】
２つの前眼部カメラの見込角（測定部１０の光軸に対する角度）が既知であり、撮影倍率も既知であるから、２つの撮影画像内のプルキンエ像の位置に基づいて、眼科装置１（撮影部１００）に対する前眼部に形成されたプルキンエ像の相対位置（実空間における３次元位置）を求めることができる。
【００８５】
また、２つの撮影画像のそれぞれにおける被検眼Ｅの特徴位置とプルキンエ像との相対位置（ズレ量）に基づいて、被検眼Ｅの特徴位置と前眼部に形成されたプルキンエ像との間の相対位置を求めることができる。
【００８６】
プルキンエ像特定部７１Ａは、上記のように特定されたプルキンエ像の位置を特定する。プルキンエ像の位置は、少なくともｘ方向の位置（ｘ座標値）及びｙ方向の位置（ｙ座標値）を含んでよく、更にｚ方向の位置（ｚ座標値）を含んでもよい。
【００８７】
（瞳孔中心特定部７１Ｂ）
瞳孔中心特定部７１Ｂは、前眼部カメラにより得られた各撮影画像（又は歪曲収差が補正された画像）を解析することで、前眼部の所定の特徴位置に相当する当該撮影画像中の位置を特定する。この実施形態では、被検眼Ｅの瞳孔中心が特定される。なお、瞳孔中心として、瞳孔の重心を求めてもよい。また、瞳孔中心（瞳孔重心）以外の特徴位置を特定するように構成することもできる。
【００８８】
瞳孔中心特定部７１Ｂは、撮影画像の画素値（輝度値など）の分布に基づいて、被検眼Ｅの瞳孔に相当する画像領域（瞳孔領域）を特定する。一般に瞳孔は他の部位よりも低い輝度で表現されるので、低輝度の画像領域を探索することによって瞳孔領域を特定することができる。このとき、瞳孔の形状を考慮して瞳孔領域を特定するようにしてもよい。つまり、略円形かつ低輝度の画像領域を探索することによって瞳孔領域を特定するように構成することができる。
【００８９】
次に、瞳孔中心特定部７１Ｂは、特定された瞳孔領域の中心位置を特定する。上記のように瞳孔は略円形であるので、瞳孔領域の輪郭を特定し、この輪郭の近似楕円の中心位置を特定し、これを瞳孔中心とすることができる。また、瞳孔領域の重心を求め、この重心位置を瞳孔中心としてもよい。
【００９０】
なお、他の特徴位置が適用される場合であっても、上記と同様に撮影画像の画素値の分布に基づいて当該特徴位置を特定することが可能である。
【００９１】
瞳孔中心特定部７１Ｂは、２つの前眼部カメラの位置（及び撮影倍率）と、特定された２つの撮影画像中の瞳孔中心の位置とに基づいて、被検眼Ｅの瞳孔中心の３次元位置を特定する。
【００９２】
例えば、２つの前眼部カメラによる撮影画像の分解能は次式で表される。
【００９３】
ｘｙ方向の分解能（平面分解能）：Δｘｙ＝Ｈ×Δｐ／ｆ
ｚ方向の分解能（奥行き分解能）：Δｚ＝Ｈ×Ｈ×Δｐ／（Ｂ×ｆ）
【００９４】
ここで、「Ｂ」は２つの前眼部カメラの間の距離（基線長）を表し、「Ｈ」は２つの前眼部カメラの基線と、被検眼Ｅの瞳孔中心との間の距離（撮影距離）を表し、「ｆ」は各前眼部カメラとその画面平面との間の距離（画面距離）を表し、「Δｐ」は画素分解能を表す。
【００９５】
瞳孔中心特定部７１Ｂは、２つの前眼部カメラの位置（既知である）と、２つの撮影画像において瞳孔中心に相当する位置とに対して、公知の三角法を適用することにより、瞳孔中心の３次元位置を算出する。
【００９６】
（移動目標位置決定部７１Ｃ）
移動目標位置決定部７１Ｃは、プルキンエ像特定部７１Ａにより特定されたプルキンエ像の位置と、瞳孔中心特定部７１Ｂにより特定された瞳孔中心位置とに基づいて、測定部１０（装置光学系）の移動目標位置を決定する。例えば、移動目標位置決定部７１Ｃは、特定されたプルキンエ像の位置と、特定された瞳孔中心位置との差分を求め、求められた差分が既定のアライメント完了条件を満たすように移動目標位置を決定する。
【００９７】
後述の制御部８０は、移動目標位置決定部７１Ｃにより決定された移動目標位置に基づいて移動機構９０を制御することが可能である。
【００９８】
（解析部７２）
解析部７２は、図２に示すように、層領域特定部７２１と、特定部７２２とを含む。層領域特定部７２１は、取得された被検眼Ｅの断層像における所定の層領域（所定の組織）を特定する。特定部７２２は、層領域特定部７２１により特定された所定の層領域に基づいて眼底Ｅｆの形状を特定する。
【００９９】
（層領域特定部７２１）
層領域特定部７２１は、被検眼Ｅの断層像を解析することによって眼底Ｅｆの所定の層領域を特定する。眼底Ｅｆの層領域としては、内境界膜、神経繊維層、神経節細胞層、内網状層、内顆粒層、外網状層、外顆粒層、外境界膜、視細胞層、網膜色素上皮層、脈絡膜、強膜、各層領域の界面などがある。
【０１００】
断層像から所定の層領域を特定する処理は、典型的には、セグメンテーション処理を含む。セグメンテーション処理は、断層像中の部分領域を特定するための公知の処理である。層領域特定部７２１は、例えば、断層像における各画素の輝度値に基づきセグメンテーション処理を行う。すなわち、眼底Ｅｆのそれぞれの層領域は特徴的な反射率を有し、これら層領域に相当する画像領域もそれぞれ特徴的な輝度値を有する。層領域特定部７２１は、これら特徴的な輝度値に基づきセグメンテーション処理を実行することにより、目的の画像領域（層領域）を特定することができる。
【０１０１】
層領域特定部７２１は、特定された所定の層領域の形状を表すデータを当該層領域の形状プロファイルとして出力する。いくつかの実施形態では、形状プロファイルは、ｘ方向、ｙ方向、及びｚ方向の少なくとも１つの方向の眼底Ｅｆの形状の変化を表す１次元、２次元、又は３次元の形状データである。
【０１０２】
例えば、層領域特定部７２１は、網膜色素上皮層（又は、ＯＳ−ＲＰＥ界面）を特定することができる。
【０１０３】
（特定部７２２）
特定部７２２は、層領域特定部７２１により得られた形状データ（形状プロファイル）から眼底Ｅｆの形状を特定する。具体的には、特定部７２２は、ＯＣＴ計測を実行することにより取得された断層像から得られた所定の層領域の形状データから真の形状を表す形状データを特定する。特定部７２２は、例えば、形状データから得られたアライメントエラーに対して感度が低い成分に基づいて真の形状を表す形状データを求める。いくつかの実施形態では、特定部７２２は、１次元（２次元又は３次元）で層領域の形状を表す形状プロファイルから、１次元（２次元又は３次元）で層領域の真の形状を表す新たな形状プロファイルを生成する。
【０１０４】
このような特定部７２２は、図４に示すように、近似処理部７２２Ａと、感度成分特定部７２２Ｂとを含む。
【０１０５】
（近似処理部７２２Ａ）
近似処理部７２２Ａは、層領域特定部７２１により特定された層領域の形状を表す形状プロファイルを、被検眼Ｅに対するＯＣＴ部３０のアライメントエラー（アライメントのずれ）に対する高感度成分及び低感度成分を含み、被検眼Ｅに対するＯＣＴ部３０の位置を変数とする式に近似する。ここで、高感度成分は、形状プロファイルのうちアライメントエラーの微少変化に対してその変動が大きい成分を意味する。低感度成分は、形状プロファイルのうちアライメントエラーの微少変化に対してその変動が小さい成分を意味する。アライメントエラーは、例えば、アライメントが完了したときのアライメント基準位置を基準とする被検眼Ｅに対するＯＣＴ部３０の測定光軸の相対位置（ｘｙ方向の相対位置）により表される。
【０１０６】
いくつかの実施形態では、低感度成分は、アライメントエラーの所定の変化量に対して、その変動量が第１閾値ＴＨ１（０＜ＴＨ１）以下の成分である。いくつかの実施形態では、高感度成分は、アライメントエラーの所定の変化量に対して、その変動量が第２閾値ＴＨ２（ＴＨ１≦ＴＨ２）を超える成分である。
【０１０７】
いくつかの実施形態では、形状プロファイルは、高感度成分と低感度成分とにより表される式に近似される。いくつかの実施形態では、形状プロファイルは、高感度成分と低感度成分とそれ以外の成分とにより表される式に近似される。
【０１０８】
また、アライメント基準位置を基準として被検眼Ｅに対するＯＣＴ部３０の測定光軸の相対位置を変数として形状プロファイルを表現した場合、形状プロファイルの非対称成分は高感度成分である。これに対して、形状プロファイルの対称成分は低感度成分である。高感度成分には、上記変数の奇数次成分がある。高感度成分の例として、傾きを表す成分、歪みを表す成分などがある。低感度成分には、上記変数の偶数次成分がある。低感度成分の例として、曲率を表す成分などがある。
【０１０９】
この実施形態では、近似処理部７２２Ａは、アライメント基準位置を基準に、被検眼Ｅに対するＯＣＴ部３０の測定光軸の相対位置の成分を変数として形状プロファイルを多項式近似することが可能である。上記の変数をｄとし、多項式の各次数の係数をｃ０、ｃ１、・・・とすると、形状プロファイルに対応した多項式ｆ（ｄ）は、次式のように表される。
【０１１０】
ｆ（ｄ）＝ｃ０＋ｃ１×ｄ＋ｃ２×ｄ２＋ｃ３×ｄ３＋ｃ４×ｄ４＋・・・
【０１１１】
上記の式において、傾きを表す成分は、ｄの１次成分（すなわち、ｃ１×ｄ）に相当する。曲率を表す成分は、ｄの２次成分（ｃ２×ｄ２）に相当する。歪みを表す成分は、ｄの３次成分（ｃ３×ｄ３）に相当する。
【０１１２】
このように、近似処理部７２２Ａにより変数ｄとして形状プロファイルを多項式近似することにより、変数ｄに起因する高感度成分及び低感度成分を特定することが可能になる。
【０１１３】
（感度成分特定部７２２Ｂ）
感度成分特定部７２２Ｂは、近似処理部７２２Ａによる近似処理により得られた形状プロファイルを表す式から低感度成分を特定し、特定された低感度成分に基づいて所定の層領域の真の形状を表す形状プロファイル（形状データ）を生成する。
【０１１４】
例えば、近似処理部７２２Ａにより上記の多項式ｆ（ｄ）が得られた場合、感度成分特定部７２２Ｂは、ｄの偶数次成分（例えば、ｄの２次成分）を低感度成分として特定し、特定された低感度成分を所定の層領域の真の形状を表す形状プロファイルとして出力する。
【０１１５】
図６及び図７に、実施形態に係る感度成分特定部７２２Ｂの動作説明図を示す。図６は、層領域特定部７２１によるセグメンテーション処理により得られた形状プロファイルの一例を表す。図６及び図７では、横軸は上記の変数ｄ（例えば、ｘ方向の位置、スキャン位置）を表し、縦軸は深さ位置（例えば、ｚ方向の位置）を表す。図７は、感度成分特定部７２２Ｂにより得られた形状プロファイルの一例を表す。
【０１１６】
層領域特定部７２１は、取得された断層像に対してセグメンテーション処理を施すことにより、所定の層領域（例えば、網膜色素上皮層）の形状の変化を示す形状プロファイルを取得する（例えば、図６の形状プロファイルＰＦ１）。近似処理部７２２Ａは、図６に示す形状プロファイルＰＦ１を多項式近似することにより、上記の多項式ｆ（ｄ）に示すような多項式を取得する。図６に示す形状プロファイルＰＦ１及び多項式ｆ（ｄ）では、アライメントエラーに起因して所定の層領域に傾き成分が含まれている。
【０１１７】
感度成分特定部７２２Ｂは、多項式ｆ（ｄ）から偶数次成分（例えば、変数ｄの２次成分）を低感度成分として抽出し、抽出された偶数次成分を所定の層領域の真の形状を表す形状プロファイルＰＦ２（図７）として出力する。図７に示す形状プロファイルＰＦ２では、層領域の傾き成分が除去されている。
【０１１８】
いくつかの実施形態では、特定部７２２は、層領域特定部７２１により得られた形状プロファイルＰＦ１に対してルジャンドル（Ｌｅｇｅｎｄｒｅ）多項式展開、フーリエ変換、又はウェーブレット変換を行う。この場合でも、アライメントエラーに対する高感度成分及び低感度成分を特定することが可能である。例えば、アライメント基準位置を基準として被検眼Ｅに対するＯＣＴ部３０の測定光軸の相対位置の成分を変数としてルジャンドル多項式展開で得られた奇数次のルジャンドル多項式は高感度成分であり、偶数次のルジャンドル多項式は低感度成分である。
【０１１９】
（算出部７３）
算出部７３は、被検眼Ｅを他覚的に測定することにより得られた屈折度数を求め、求められた屈折度数と特定部７２２により特定された眼底Ｅｆの形状（形状プロファイル）とに基づいて、被検眼Ｅの中心窩を含む領域の周辺領域の屈折度数を算出する。いくつかの実施形態では、算出部７３は、求められた屈折度数と、特定部７２２により特定された眼底Ｅｆの形状に対応した被検眼の光学特性を表すパラメータとに基づいて、被検眼の中心窩を含む領域の周辺領域の屈折度数を算出する。算出部７３は、特定部７２２により特定された眼底Ｅｆの形状に対応した被検眼の光学特性を表すパラメータに基づいて眼球モデルを構築し、構築された眼球モデルと求められた屈折度数とから、上記の周辺領域の屈折度数を算出することが可能である。
【０１２０】
算出部７３は、図５に示すように、屈折度数算出部７３Ａと、眼球モデル構築部７３Ｂと、周辺屈折度数算出部７３Ｃとを含む。
【０１２１】
（屈折度数算出部７３Ａ）
屈折度数算出部７３Ａは、屈折測定部２０の受光系の撮像素子からの出力を処理して屈折度数を算出する。
【０１２２】
いくつかの実施形態では、屈折度数算出部７３Ａは、撮像素子によって取得されたリングパターン像を楕円近似することにより楕円形状を特定する処理と、特定された楕円形状と合焦レンズ等に対するフォーカス調整分のディオプタとに基づいて屈折度数（測定データ）を求める処理とを実行する。
【０１２３】
いくつかの実施形態では、屈折度数算出部７３Ａは、撮像素子によって取得されたリングパターン像が描出された画像における輝度分布を求める処理と、求められた輝度分布からリングパターン像の重心位置を求める処理と、求められた重心位置から放射状に延びる複数の走査方向に沿った輝度分布を求める処理と、求められた複数の走査方向に沿った輝度分布からリングパターン像を特定する処理と、特定されたリングパターン像の近似楕円を求める処理と、求められた近似楕円の長径及び短径を公知の式に代入することによって屈折度数を算出する処理とを実行する。
【０１２４】
いくつかの実施形態では、屈折度数算出部７３Ａは、撮像素子によって取得されたリングパターン像の基準パターンからの偏位（位置ずれ、変形等）を求める処理と、この偏位から屈折度数を求める処理とを実行する。
【０１２５】
いくつかの実施形態では、屈折度数として球面度数Ｓ、乱視度数Ｃ及び乱視軸角度Ａが算出される。いくつかの実施形態では、屈折度数として等価球面度数ＳＥ（＝Ｓ＋Ｃ／２）が算出される。
【０１２６】
（眼球モデル構築部７３Ｂ）
眼球モデル構築部７３Ｂは、眼球モデルを構築する。眼球モデル構築部７３Ｂは、公知の模型眼等の眼球モデルに対して、別途に取得されたパラメータを適用することにより新たな眼球モデルを構築することが可能である。
【０１２７】
眼球モデル構築部７３Ｂは、公知の模型眼等の眼球モデルに対して、ＯＣＴ計測等により得られた被検眼Ｅの眼内距離を実測パラメータとして適用することにより新たな眼球モデルを構築することが可能である。この場合、データ処理部７０は、組織のサイズ（層厚、体積等）や所定の部位間の距離を求めるための算出処理などを実行することが可能である。例えば、データ処理部７０は、スキャンデータ又は断層像を解析することにより眼内の所定部位に相当する干渉光の検出結果（干渉信号）のピーク位置を特定し、特定されたピーク位置間の距離に基づいて眼内距離を求める。いくつかの実施形態では、データ処理部７０は、セグメンテーション処理によって得られた２つの層領域の間に存在するピクセルの個数と、所定のピクセルスペーシング補正値とに基づいて眼内距離（層間距離）求める。眼内距離の計測は、所定の方向に沿って行われる。眼内距離の計測方向は、例えば、ＯＣＴスキャンによって決定される方向（例えば、測定光の進行方向）でもよいし、スキャンデータに基づき決定される方向（例えば、層に直交する方向）でもよい。また、距離データは、２つの層領域間の距離分布データでもよいし、この距離分布データから算出された統計値（例えば、平均、最大値、最小値、中央値、最頻値、分散、標準偏差）でもよいし、各層領域の代表点の間の距離データでもよい。データ処理部７０により算出可能な眼内距離には、眼軸長、角膜厚、前房深度、水晶体厚、硝子体腔長、網膜厚、脈絡膜厚などがある。更に、データ処理部７０は、求められた眼内距離を用いて、眼球の光学特性を表す各種のパラメータを算出することが可能である。
【０１２８】
いくつかの実施形態では、特定部７２２（又は眼球モデル構築部７３Ｂ）は、構築された眼球モデルを用いて眼底Ｅｆの形状を特定することが可能である。例えば、特定部７２２は、眼底Ｅｆにおける中心領域と周辺領域との深さ位置の差分を求めることにより、眼底Ｅｆの形状を特定する。
【０１２９】
（周辺屈折度数算出部７３Ｃ）
周辺屈折度数算出部７３Ｃは、眼底Ｅｆにおける中心窩を含む中心領域の外側の周辺領域の屈折度数を算出する。このとき、周辺屈折度数算出部７３Ｃは、屈折測定部２０により得られた中心領域の屈折度数と、特定された眼底Ｅｆの形状とに基づいて、周辺領域の屈折度数を算出する。周辺屈折度数算出部７３Ｃは、眼球モデル構築部７３Ｂにより構築された眼球モデルのパラメータを用いて周辺領域の屈折度数を算出することが可能である。
【０１３０】
いくつかの実施形態では、データ処理部７０の機能は１以上のプロセッサにより実現される。いくつかの実施形態では、アライメント処理部７１、解析部７２、及び算出部７３のそれぞれの機能は単独のプロセッサにより実現される。いくつかの実施形態では、アライメント処理部７１の各部の機能は単独のプロセッサにより実現される。いくつかの実施形態では、解析部７２の機能は単独のプロセッサにより実現される。いくつかの実施形態では、算出部７３の各部の機能は単独のプロセッサにより実現される。いくつかの実施形態では、データ処理部７０の少なくとも一部が屈折測定部２０又はＯＣＴ部３０に設けられる。
【０１３１】
（制御部８０）
制御部８０は、眼科装置１の各部を制御する。制御部８０は、記憶部（不図示）を含み、各種の情報を保存することが可能である。記憶部に保存される情報には、眼科装置１の各部を制御するためのプログラム、被検者の情報、被検眼の情報、測定部１０により得られた測定データ、データ処理部７０による処理結果などがある。制御部８０の機能は、プロセッサにより実現される。
【０１３２】
制御部８０は、図示しない表示デバイスを制御可能である。表示デバイスは、ユーザインターフェイスの一部として機能し、制御部８０による制御を受けて情報を表示する。表示デバイスは、例えば、液晶ディスプレイ（ＬＣＤ）又は有機発光ダイオード（ＯＬＥＤ）ディスプレイであってよい。
【０１３３】
制御部８０は、図示しない操作デバイスからの信号にしたがって眼科装置１を制御可能である。操作デバイスは、ユーザインターフェイス部の一部として機能する。操作デバイスは、眼科装置１に設けられた各種のハードウェアキー（ジョイスティック、ボタン、スイッチなど）を含んでいてよい。また、操作デバイスは、眼科装置１に接続された各種の周辺機器（キーボード、マウス、ジョイスティック、操作パネルなど）を含んでいてよい。また、操作デバイスは、タッチパネルに表示される各種のソフトウェアキー（ボタン、アイコン、メニューなど）を含んでよい。
【０１３４】
（移動機構９０）
移動機構９０は、屈折測定部２０、ＯＣＴ部３０、アライメント光投射部４０、ビームスプリッタＢＳ１、ＢＳ２等の光学系（装置光学系）が収納されたヘッド部を上左右方向及び前後方向に移動させるための機構である。移動機構９０は、制御部８０からの制御を受け、被検眼Ｅに対して測定部１０を相対移動させることが可能である。例えば、移動機構９０には、ヘッド部を移動するための駆動力を発生するアクチュエータと、この駆動力を伝達する伝達機構とが設けられる。アクチュエータは、例えばパルスモータにより構成される。伝達機構は、例えば歯車の組み合わせやラック・アンド・ピニオンなどによって構成される。制御部８０は、アクチュエータに対して制御信号を送ることにより移動機構９０に対する制御を行う。
【０１３５】
移動機構９０に対する制御は、位置合わせ（アライメント）において用いられる。例えば、制御部８０は、測定部１０（装置光学系）の現在位置を取得する。制御部８０は、移動機構９０の移動制御の内容を表す情報を受けて、測定部１０の現在位置を取得する。この場合、制御部８０は、所定のタイミング（装置起動時、患者情報入力時など）で移動機構９０を制御して、測定部１０を所定の初期位置に移動させる。それ以降、制御部８０は、移動機構９０を制御する度に、その制御内容を記録する。それにより、制御内容の履歴が得られる。制御部８０は、光学系位置取得部として、この履歴を参照して現在までの制御内容を取得し、この制御内容に基づいて測定部１０の現在位置を求める。
【０１３６】
いくつかの実施形態では、制御部８０が移動機構９０を制御する度にその制御内容を受けて測定部１０の現在位置を逐次求める。いくつかの実施形態では、眼科装置１には、測定部１０の位置を検知する位置センサーが設けられ、制御部８０は、位置センサーによる検出結果に基づいて測定部１０の現在位置を取得する。
【０１３７】
制御部８０は、上記のように取得された現在位置と、移動目標位置決定部７１Cにより決定された移動目標位置とに基づいて、移動機構９０を制御することができる。それにより、測定部１０を移動目標位置に移動させることができる。例えば、制御部８０は、現在位置と移動目標位置との差分を求める。この差分値は、例えば、現在位置を始点とし、移動目標位置を終点とするベクトル値である。このベクトル値は、例えば、ｘｙｚ座標系で表現される３次元ベクトル値である。
【０１３８】
いくつかの実施形態では、移動機構９０に対する制御は、トラッキングにおいて用いられる。トラッキングとは、被検眼Ｅの眼球運動に合わせて装置光学系を移動させるものである。トラッキングを行う場合には、事前にアライメントとフォーカス調整が実行される。トラッキングは、装置光学系の位置を眼球運動に追従させることにより、アライメントとピントが合った好適な位置関係を維持する機能である。
【０１３９】
制御処理部５０、又はデータ処理部７０は、実施形態に係る「眼科情報処理装置」の一例である。ＯＣＴ部３０及び画像形成部６０、外部装置（外部の眼科装置）、又は記録媒体からデータを受信するデバイス（通信インターフェイス、入出力インターフェイス等）は、実施形態に係る「取得部」の一例である。ＯＣＴ部３０及び画像形成部６０は、実施形態に係る「ＯＣＴ部」の一例である。層領域特定部７２１により得られた形状プロファイルは、実施形態に係る「第１形状データ」の一例である。層領域特定部７２１は、実施形態に係る「組織特定部」の一例である。感度成分特定部７２２Ｂにより得られた形状プロファイルは、実施形態に係る「第２形状データ」の一例である。特定部７２２又は感度成分特定部７２２Ｂは、実施形態に係る「特定部」の一例である。
【０１４０】
＜動作例＞
実施形態に係る眼科装置１の動作について説明する。
【０１４１】
図８及び図９に、眼科装置１の動作の一例を示す。図８は、眼科装置１の動作例のフロー図を表す。図９は、図８のステップＳ３の動作例のフロー図を表す。制御部８０の記憶部には、図８及び図９に示す処理を実現するためのコンピュータプログラムが記憶されている。制御部８０は、このコンピュータプログラムに従って動作することにより、図８及び図９に示す処理を実行する。
【０１４２】
（Ｓ１：アライメント）
まず、制御部８０は、アライメントを実行する。
【０１４３】
例えば、制御部８０は、アライメント光投射部４０を制御して、被検眼Ｅにアライメント光を投射させる。このとき、被検眼Ｅには、図示しない固視投影系により所定の投射位置（例えば、測定光軸上の投射位置）に固視光束が投射される。制御部８０は、例えば、撮影部１００により取得された撮影画像中の瞳孔中心位置とプルキンエ像の位置との変位から測定部１０の移動量及び移動方向を特定し、特定された移動量及び移動方向に基づいて移動機構９０を制御し、被検眼Ｅに対する測定部１０の位置合わせを行う。制御部８０は、所定のアライメント完了条件を満足するまでこの処理を繰り返し実行する。
【０１４４】
（Ｓ２：他覚屈折測定）
次に、制御部８０は、図示しない固視投影系を制御して、眼底Ｅｆにおける屈折測定部２０の光学系の測定光軸上に固視標を投影させる（中心固視）。その後、制御部８０は、屈折測定部２０の光学系の測定光軸上に固視標が投影された状態で、屈折測定部２０を制御することにより他覚屈折測定を実行する。
【０１４５】
屈折度数算出部７３Ａは、被検眼Ｅの眼底Ｅｆに投射された光の反射光により形成されたリングパターン像を解析することにより被検眼Ｅの中心窩を含む中心領域の屈折度数を算出する。
【０１４６】
（Ｓ３：眼底の形状を特定）
続いて、制御部８０は、被検眼Ｅの眼底Ｅｆの形状の特定処理を実行する。この実施形態では、制御部８０は、屈折測定部２０（ＯＣＴ部３０）の光学系の測定光軸上に固視標が投影された状態で、ＯＣＴ部３０を制御することによりＯＣＴ計測（ＯＣＴスキャン）を実行する。
【０１４７】
ステップＳ３では、上記のように、眼底Ｅｆの形状を表す形状データが取得される。ステップＳ３の詳細については後述する。
【０１４８】
（Ｓ４：周辺屈折度数を算出）
続いて、制御部８０は、ステップＳ２において得られた中心窩を含む中心領域の外側の周辺領域の屈折度数を周辺屈折度数算出部７３Ｃに算出させる。そのため、制御部８０は、眼球モデルを眼球モデル構築部７３Ｂに構築させる。
【０１４９】
具体的には、眼球モデル構築部７３Ｂは、ステップＳ３において得られた形状データから所定の層領域のＨｅｉｇｈｔデータ［ｐｉｘｅｌ］を取得する。Ｈｅｉｇｈｔデータは、断層像において所定の基準位置からの深さ方向の距離に対応する。眼球モデル構築部７３Ｂは、光学系で規定される装置固有のピクセルスペーシング補正値［ｍｍ／ｐｉｘｅｌ］を用いて、Ｈｅｉｇｈｔデータの距離［ｍｍ］を取得する。更に、眼球モデル構築部７３Ｂは、取得されたＨｅｉｇｈｔデータを眼底形状データとして眼球モデルを構築する。
【０１５０】
図１０に、実施形態に係る眼球モデル構築部７３Ｂの動作説明図を示す。図１０は、眼球モデルのパラメータの一部を模式的に表す。
【０１５１】
眼球モデル構築部７３Ｂは、Ｇｕｌｌｓｔｒａｎｄ模型眼等の眼球モデルのパラメータを用いて、所定の角膜曲率半径（例えば、７．７ｍｍ）、所定の眼軸長（例えば、２４．２ｍｍ）で構成される眼球モデルを構築する。
【０１５２】
眼球モデル構築部７３Ｂは、図１０に示すように、眼球モデルにおいて角膜Ｅｃと眼底Ｅｆとの間に装置固有のピボット点Ｐｖを設定する。典型的には、スキャン系を構成する光スキャナーと光学的に共役な位置に配置される瞳孔位置に相当する位置（例えば角膜Ｅｃに対して後方側に３ｍｍの位置）がピボット点Ｐｖとして設定される。ピボット点Ｐｖを中心として、等距離（等光路長）の位置（ＥＬＳ）が、ＯＣＴ計測により得られる断層像中の平坦な位置に相当する。
【０１５３】
眼球モデルにおいて、眼軸長ＡＬと、角膜前面（後面）からピボット点Ｐｖまでの距離Ｌｐとが既知であるため、ピボット点Ｐｖから眼底Ｅｆまでの距離（ＡＬ−Ｌｐ）が既知となる。眼底Ｅｆの曲率半径が距離（ＡＬ−Ｌｐ）と等しいときに上記のように断層像中の平坦な位置に相当するため、眼球モデル構築部７３Ｂは、得られたＨｅｉｇｈｔデータの距離［ｍｍ］から眼底Ｅｆの形状（例えば、曲率半径）を特定することが可能である。
【０１５４】
そこで、眼球モデル構築部７３Ｂは、中心領域（中心窩）に対する周辺領域の高さの差分（眼底形状差分データ）Δｈ［ｍｍ］を求める。差分Δｈは、断層像におけるＡライン毎に求めてもよいし、多項式や非球面式（コーニック定数を含む多項式）等の任意の関数でフィッティングしてもよい。
【０１５５】
次に、周辺屈折度数算出部７３Ｃは、眼底形状と屈折度数とを関係付けるため、全眼系の眼球屈折力を定義する。典型的な眼球モデル（Ｇｕｌｌｓｔｒａｎｄ模型眼（精密模型眼、調節休止状態））では、全眼系の眼球屈折力は５８．６４［ディオプタ］である。空気換算長では、全眼系の焦点距離は「１０００／５８．６４＝１７．０５」［ｍｍ］となる。ピクセルスペーシング補正値を用いて得られる単位［ｍｍ］の情報は、通常は生体組織内(ｉｎ  ｔｉｓｓｕｅ)における距離を表すため、屈折率を乗算して生体組織内における全眼系の焦点距離が算出される。全眼系の等価屈折率をｎ＝１．３８とすると、生体組織内における全眼系の焦点距離ｆｔは、「１０００／５８．６４×１．３８＝２３．５３」［ｍｍ］となる。
【０１５６】
周辺屈折度数算出部７３Ｃは、中心領域（中心窩）に対する周辺領域の高さの差分Δｈの位置における眼球屈折力の差分ΔＤを式（１）に従って算出する。差分ΔＤは、中心窩を含む中心領域に対する相対的な眼球屈折力の差分に相当する。
【０１５７】
【数１】
2020110579.tif 000003
【０１５８】
例えば、Δｈ＝０．１［ｍｍ］（ｉｎ  ｔｉｓｓｕｅ）としたとき、ΔＤ＝０．１８［ディオプタ］となる。
【０１５９】
周辺屈折度数算出部７３Ｃは、式（２）に示すように、中心領域の等価球面度数ＳＥに対して式（１）の差分ΔＤを適用することにより、周辺領域の屈折度数ＳＥｐを求める。
【０１６０】
【数２】
2020110579.tif 000004
【０１６１】
周辺屈折度数算出部７３Ｃは、断層像における周辺領域の屈折度数をＡライン毎に求めてもよいし、任意の関数でフィッティングしてもよい。
【０１６２】
いくつかの実施形態では、周辺屈折度数算出部７３Ｃは、被検眼の中心窩を含む領域の周辺領域の複数の部位のそれぞれにおいて屈折度数を算出し、算出された複数の屈折度数の統計値を算出する。いくつかの実施形態では、複数の屈折度数は、周辺領域において中心窩を含む領域に対して線対称又は点対称の関係を有する複数の部位の屈折度数である。いくつかの実施形態では、統計値は、複数の屈折度数の平均値、最大値、最小値、中央値、又は最頻値である。
【０１６３】
以上で、眼科装置１の動作は終了である（エンド）。
【０１６４】
図８のステップＳ３では、図９に示すように、被検眼Ｅの眼底Ｅｆの形状の特定処理が実行される。
【０１６５】
（Ｓ１１：アライメント光の投射を開始）
ステップＳ３の処理が開始されると、制御部８０は、ステップＳ１と同様に、アライメント光投射部４０を制御して、被検眼Ｅに対してアライメント光の投射を開始させる。
【０１６６】
ステップＳ１１においても、ステップＳ１と同様に、被検眼Ｅには、図示しない固視投影系により所定の投射位置（例えば、測定光軸上の投射位置）に固視光束が投射される。
【０１６７】
（Ｓ１２：アライメント）
制御部８０は、撮影部１００により取得された撮影画像中の瞳孔中心位置とプルキンエ像の位置との変位から測定部１０の移動量及び移動方向を特定し、特定された移動量及び移動方向に基づいて移動機構９０を制御し、被検眼Ｅに対する測定部１０の位置合わせを行う。
【０１６８】
（Ｓ１３：アライメント完了？）
制御部８０は、所定のアライメント完了条件が満足したか否かを判定する。アライメント完了条件は、測定部１０の光軸のｘ方向及びｙ方向の位置が移動目標位置のｘ方向及びｙ方向の位置と一致することと、ｚ方向の距離が所定の作動距離になることとを含む。いくつかの実施形態では、所定の作業距離は、測定部１０（対物レンズ）の作動距離である。
【０１６９】
所定のアライメント完了条件を満足していないと判定されたとき（Ｓ１３：Ｎ）、眼科装置１の動作はステップＳ１２に移行する。所定のアライメント完了条件を満足していると判定されたとき（Ｓ１３：Ｙ）、眼科装置１の動作はステップＳ１４に移行する。
【０１７０】
（Ｓ１４：ＯＣＴ計測）
ステップＳ１３において、所定のアライメント完了条件を満足していると判定されたとき（Ｓ１３：Ｙ）、制御部８０は、ＯＣＴ部３０を制御することにより、眼底Ｅｆにおける所定部位に対してＯＣＴスキャンを実行させることでＯＣＴ計測を実行させる。所定部位として、中心窩又はその近傍がある。ＯＣＴスキャンとして、ラジアルスキャンなどがある。
【０１７１】
（Ｓ１５：断層像を形成）
続いて、制御部８０は、ステップＳ１４において得られたスキャンデータに基づいて被検眼Ｅの断層像を画像形成部６０に形成させる。
【０１７２】
（Ｓ１６：セグメンテーション処理）
次に、制御部８０は、ステップＳ１５において形成された断層像に対してセグメンテーション処理を施すことにより所定の層領域（例えば、網膜色素上皮層）を層領域特定部７２１に特定させる。それにより、所定の層領域の形状データ（形状プロファイル）が取得される。
【０１７３】
（Ｓ１７：近似処理）
次に、制御部８０は、近似処理部７２２Ａを制御することにより、ステップＳ１６において得られた形状データに対する近似処理を実行させる。近似処理部７２２Ａは、例えば形状データを多項式近似することにより多項式（例えば、上記の多項式ｆ（ｄ））を取得する。
【０１７４】
（Ｓ１８：低感度成分を特定）
制御部８０は、ステップＳ１７において取得された多項式から低感度成分を感度成分特定部７２２Ｂに特定させる。感度成分特定部７２２Ｂは、例えば、上記の多項式ｆ（ｄ）からｄの偶数次成分（例えば、ｄの２次成分）を低感度成分として特定し、特定された低感度成分を所定の層領域の真の形状を表す形状プロファイルとして出力する。図８のステップＳ４では、ステップＳ１８において出力されたプロファイルが、眼底Ｅｆの形状に対応する所定の層領域（例えば網膜色素上皮層）の形状を表すデータとして用いられる。
【０１７５】
以上で、図８のステップＳ３の処理は終了である（エンド）。
【０１７６】
＜変形例＞
実施形態に係る眼科装置の構成及び動作は、上記したものに限定されるものではない。
【０１７７】
［第１変形例］
ステップＳ４において、眼球モデル構築部７３Ｂは、Ｇｕｌｌｓｔｒａｎｄ模型眼等の眼球モデルのパラメータのうち、被検眼Ｅの実測データ（例えば、眼軸長、角膜形状、前房深度、水晶体曲率、水晶体厚の測定値）の少なくとも１つを置き換えて新たな眼球モデルを構築してもよい。いくつかの実施形態では、実測データは、外部の測定装置又は電子カルテシステムから取得される。いくつかの実施形態では、眼軸長、前房深度、水晶体曲率、及び水晶体厚は、ＯＣＴ部３０により得られたスキャンデータから求められる。
【０１７８】
例えば、周辺屈折度数算出部７３Ｃ（又はデータ処理部７０）は、構築された新たな眼球モデルを用いて、角膜Ｅｃから入射して瞳孔を通過して眼底Ｅｆに到達する光線について光線追跡処理を行う（例えば、瞳孔径＝φ４）。光線追跡処理では、物点の位置を、ステップＳ２で取得された中心領域における屈折度数（等価球面度数ＳＥ）から求まる遠点に相当する位置とする。角膜Ｅｃから遠点に相当する位置までの遠点距離Ｌは、「−１０００／ＳＥ」［ｍｍ］である。
【０１７９】
まず、周辺屈折度数算出部７３Ｃは、中心領域について光線追跡処理を行う。上記のように眼球モデルに実測データを適用するため、中心領域においても眼底Ｅｆで光線が収束しない可能性がある。この場合、周辺屈折度数算出部７３Ｃは、中心領域において光線が収束するように（眼底Ｅｆの面が最良像面）となるように眼球モデルのパラメータを微調整する。
【０１８０】
次に、周辺屈折度数算出部７３Ｃは、パラメータが微調整された眼球モデルを用いて、周辺領域について光線追跡処理を行う（すなわち、眼の回旋点を通る測定光軸に対して入射角を設けた光線を追跡する）。周辺屈折度数算出部７３Ｃは、物点までの距離を変更しつつ光線追跡処理を行うことで、周辺領域において眼底Ｅｆで光線が収束するような物点までの距離を求める。求められた物点までの距離が、周辺領域における遠点距離Ｌｐに対応する。周辺屈折度数算出部７３Ｃは、式（３）を用いて周辺領域の屈折度数ＳＥｐ［ディオプタ］を求めることができる。
【０１８１】
【数３】
2020110579.tif 000005
【０１８２】
周辺屈折度数算出部７３Ｃは、所定の入射角範囲で入射角を変更しつつ光線追跡処理を行い、入射角（画角）ごとの周辺領域の屈折度数ＳＥｐを求める。周辺領域の屈折度数は、入射角ごとの離散値であってもよいし、入射角範囲で任意の関数でフィッティングしてもよい。
【０１８３】
本変形例では、中心領域において眼底Ｅｆ状で光線が収束するように眼球モデルを微調整するため、求められた周辺領域の屈折度数は、中心領域に対する相対屈折度数を求めることに相当する。
【０１８４】
［第２変形例］
上記の実施形態において、眼底Ｅｆの中心領域の形状として、上記の形状データ（形状プロファイル）又は断層像から、水平方向（所定の基準方向）に対する眼底の所定の層領域（例えば、網膜色素上皮層、ＯＳ−ＲＰＥ界面）のチルト角度が特定されてもよい。
【０１８５】
第２変形例に係る眼科装置の構成は眼球モデル構築部７３Ｂが省略された点を除いて実施形態に係る眼科装置１の構成と同様であるため、説明を省略する。
【０１８６】
本変形例では、ステップＳ３において、特定部７２２（又は周辺屈折度数算出部７３Ｃ）は、ステップＳ１５において取得された断層像から求められたＨｅｉｇｈｔデータを用いて、水平方向の断層像（Ｂスキャン画像）について眼底面のチルト角度θｈと、垂直方向のＢスキャン画像について眼底面のチルト角度θｖを算出する。
【０１８７】
チルト角度θｈ、θｖは、以下のように、チルト角度ｇ１と同様の手法で算出可能である。
【０１８８】
図１１に、水平方向の断層像を模式的に示す。
【０１８９】
図１１において、断層像ＩＭＧのフレーム左端ＬＴにおいて、フレーム上端ＵＴから眼底Ｅｆにおける所定の層領域（層領域特定部７２１により特定された層領域。例えば、網膜色素上皮層、ＯＳ−ＲＰＥ界面、又は神経線維層）に相当する部位の画像領域との垂直方向の距離をＬ１とする。同様に、断層像ＩＭＧのフレーム右端ＲＴにおいて、フレーム上端ＵＴから当該層領域に相当する部位の画像領域との垂直方向の距離をＲ１とする。距離Ｌ１は、フレーム左端ＬＴにおけるＨｅｉｇｈｔデータより求められる。距離Ｒ１は、フレーム右端ＲＴにおけるＨｅｉｇｈｔデータより求められる。特定部７２２は、断層像ＩＭＧにおけるフレーム左端ＬＴとフレーム右端ＲＴにおける当該部位の画像領域の垂直方向の差分（｜Ｒ１−Ｌ１｜）について実寸法に相当する値｜ｄ｜を求める。
【０１９０】
次に、特定部７２２は、ＯＣＴ計測範囲に相当する断層像ＩＭＧのフレームの水平方向の距離Ｈ１について実寸法に相当する値ｃを求める。例えば、水平方向のスキャン範囲の長さをピクセルスペーシング補正値［ｍｍ／ｐｉｘｅｌ］を用いて値ｃが特定される。
【０１９１】
特定部７２２は、傾斜角度ｇ０［ｄｅｇｒｅｅ］を式（４）に従って求める。
【０１９２】
【数４】
2020110579.tif 000006
【０１９３】
いくつかの実施形態では、特定部７２２は、測定光軸と眼球光軸とのずれ量に応じて傾斜角度ｇ０を補正することにより眼底面のチルト角度を求める。
【０１９４】
（測定光軸と眼球光軸とが略一致しているとき）
測定光軸と眼球光軸（視軸）とが略一致しているとき、特定部７２２は、式（５）に示すように、断層像の傾斜角度ｇ０を補正することなく眼底面のチルト角度ｇ１として出力する。
【０１９５】
【数５】
2020110579.tif 000007
【０１９６】
（測定光軸に対して眼球光軸がシフトしているとき）
測定光軸に対して眼球光軸がシフトしているとき、特定部７２２は、シフト量に基づいて断層像の傾斜角度ｇ０を補正することにより眼底面のチルト角度ｇ１を求める。
【０１９７】
例えば、特定部７２２は、式（６）に示すシフト量ｄｓを変数とする一次式に従って補正角度φ１を求め、式（７）に示すように、求められた補正角度φ１を用いて傾斜角度ｇ０を補正することで眼底面のチルト角度ｇ１を求める。式（６）において、α１及びｃ１は定数である。例えば、模型眼データを用いてα１及びｃ１を求めることができる。
【０１９８】
【数６】
2020110579.tif 000008
【０１９９】
【数７】
2020110579.tif 000009
【０２００】
（測定光軸に対して眼球光軸がチルトしているとき）
測定光軸に対して眼球光軸がチルトしているとき、特定部７２２は、チルト量に基づいて断層像の傾斜角度ｇ０を補正することにより眼底面のチルト角度ｇ１を求める。
【０２０１】
例えば、特定部７２２は、式（８）に示すようなチルト量ｄｔを変数とする一次式に従って補正角度φ２を求め、式（９）に示すように、求められた補正角度φ２を用いて傾斜角度ｇ０を補正することで眼底面のチルト角度ｇ１を求める。式（８）において、α２及びｃ２は定数である。例えば、模型眼データを用いてα２及びｃ２を求めることができる。
【０２０２】
【数８】
2020110579.tif 000010
【０２０３】
【数９】
2020110579.tif 000011
【０２０４】
（測定光軸に対して眼球光軸がシフトし、且つチルトしているとき）
測定光軸に対して眼球光軸がシフトし、且つチルトしているとき、特定部７２２は、シフト量及びチルト量に基づいてＢスキャン画像の傾斜角度ｇ０を補正することにより眼底面のチルト角度ｇ１を求める。
【０２０５】
例えば、シフト量ｄｓ及びチルト量ｄｔが小さい範囲において、特定部７２２は、式（１０）に示すようなシフト量ｄｓ及びチルト量ｄｔを変数とする式に従って補正角度φ３を求め、式（１１）に示すように、求められた補正角度φ３を用いて傾斜角度ｇ０を補正することで眼底面のチルト角度ｇ１を求める。いくつかの実施形態では、式（１０）は、シフト量の補正角度を求める式と、チルト量の補正角度を求める式とを線形結合することにより得られる結合式である。式（１０）において、α３、α４及びｃ３は定数である。例えば、模型眼データを用いてα３、α４及びｃ３を求めることができる。
【０２０６】
【数１０】
2020110579.tif 000012
【０２０７】
【数１１】
2020110579.tif 000013
【０２０８】
本変形例では、屈折度数算出部７３Ａは、水平方向及び垂直方向それぞれについて、上記のように特定された眼底面のチルト角度θｈ、θｖに応じて、ステップＳ２において取得されたリングパターン像を補正する。屈折度数算出部７３Ａは、補正されたリングパターン像に対して楕円近似を行い、得られた楕円形状を用いて公知の手法で屈折度数を求める。求められた屈折度数が、中心領域の屈折度数として算出される。
【０２０９】
例えば、眼底面のチルト角度が０度のときに取得されたリングパターン像の長径をＬＡとし、短径をＬＢとする。長径方向に眼底面が傾斜してチルト角度がθ度であるとき、取得されたリングパターン像から近似される楕円の長径はＬＡ／ｃｏｓθとなり、短径はＬＢとなる。従って、屈折度数算出部７３Ａは、ステップＳ２において取得されたリングパターン像を楕円近似して得られた楕円の長径方向にｃｏｓθを乗算することで、リングパターン像を補正することができる。短径方向にチルトしている場合も同様である。例えば、屈折度数算出部７３Ａは、水平方向及び垂直方向それぞれのチルト角度から、楕円の長径方向のチルト角度及び短径方向のチルト角度を求めることで、リングパターン像を補正することができる。
【０２１０】
周辺屈折度数算出部７３Ｃは、上記の実施形態と同様に、式（２）に示すように、中心領域の等価球面度数ＳＥに対して式（１）の差分ΔＤを適用することにより、周辺領域の屈折度数ＳＥｐを求めることが可能である。
【０２１１】
［第３変形例］
上記の実施形態又はその変形例では、近似処理部７２２Ａにより得られた形状データ（形状プロファイル）からアライメントエラーに対する低感度成分を抽出する場合について説明した。しかしながら、実施形態に係る眼科装置１の構成はこれに限定されるものではない。例えば、形状データからアライメントエラーに対する高感度成分を抽出し、抽出された高感度成分を形状データから除去（分離）することにより得られた新たな形状データを所定の層領域の真の形状を表す形状データとして取得してもよい。
【０２１２】
以下、第３変形例に係る眼科装置について、実施形態に係る眼科装置１との相違点を中心に説明する。
【０２１３】
第３変形例では、感度成分特定部７２２Ｂは、近似処理部７２２Ａによる近似処理により得られた形状プロファイルを表す式から高感度成分を特定し、形状プロファイルを表す式から高感度成分が除去された新たな形状プロファイルを所定の層領域の真の形状を表す形状プロファイルとして生成する。
【０２１４】
例えば、近似処理部７２２Ａにより上記の多項式ｆ（ｄ）が得られた場合、感度成分特定部７２２Ｂは、ｄの奇数次成分（例えば、ｄの１次成分又は３次成分）を高感度成分として特定し、特定された高感度成分を上記の多項式ｆ（ｄ）から差し引くことで所定の層領域の真の形状を表す形状プロファイルとして出力する。
【０２１５】
第３変形例では、図８のステップＳ３において、次にように被検眼Ｅの眼底Ｅｆの形状の特定処理が実行される。
【０２１６】
図１２に、実施形態の第３変形例に係る眼科装置の動作の一例を示す。図１２は、図８のステップＳ３の動作例のフロー図を表す。制御部８０の記憶部には、図１２に示す処理を実現するためのコンピュータプログラムが記憶されている。制御部８０は、このコンピュータプログラムに従って動作することにより、図１２に示す処理を実行する。
【０２１７】
（Ｓ２１：アライメント光の投射を開始）
ステップＳ３の処理が開始されると、制御部８０は、ステップＳ１１と同様に、アライメント光投射部４０を制御して、被検眼Ｅに対してアライメント光の投射を開始させる。
【０２１８】
（Ｓ２２：アライメント）
制御部８０は、ステップＳ１２と同様に、移動機構９０を制御し、被検眼Ｅに対する測定部１０の位置合わせを行う。
【０２１９】
（Ｓ２３：アライメント完了？）
制御部８０は、ステップＳ１３と同様に、所定のアライメント完了条件が満足したか否かを判定する。
【０２２０】
所定のアライメント完了条件を満足していないと判定されたとき（Ｓ２３：Ｎ）、眼科装置１の動作はステップＳ２２に移行する。所定のアライメント完了条件を満足していると判定されたとき（Ｓ２３：Ｙ）、眼科装置１の動作はステップＳ２４に移行する。
【０２２１】
（Ｓ２４：ＯＣＴ計測）
ステップＳ２３において、所定のアライメント完了条件を満足していると判定されたとき（Ｓ２３：Ｙ）、制御部８０は、ステップＳ１３と同様に、ＯＣＴ部３０を制御することにより、眼底Ｅｆにおける所定部位に対してＯＣＴスキャンを実行させることでＯＣＴ計測を実行させる。
【０２２２】
（Ｓ２５：断層像を形成）
続いて、制御部８０は、ステップＳ１５と同様に、ステップＳ２４において得られたスキャンデータに基づいて被検眼Ｅの断層像を画像形成部６０に形成させる。
【０２２３】
（Ｓ２６：セグメンテーション処理）
次に、制御部８０は、ステップＳ１６と同様に、ステップＳ２５において形成された断層像に対してセグメンテーション処理を施すことにより所定の層領域（例えば、網膜色素上皮層）を層領域特定部７２１に特定させる。
【０２２４】
（Ｓ２７：近似処理）
次に、制御部８０は、ステップＳ１７と同様に、近似処理部７２２Ａを制御することにより、ステップＳ２６において得られた形状データに対する近似処理を実行させる。
【０２２５】
（Ｓ２８：高感度成分を特定）
制御部８０は、ステップＳ２７において取得された多項式から高感度成分を感度成分特定部７２２Ｂに特定させる。感度成分特定部７２２Ｂは、例えば、上記の多項式ｆ（ｄ）からｄの奇数次成分（例えば、ｄの１次成分及び３次成分の少なくとも一方）を高感度成分として特定する。
【０２２６】
（Ｓ２９：高感度成分を除去）
制御部８０は、ステップＳ２７において得られた形状データからステップＳ２８において特定された高感度成分を除去するように感度成分特定部７２２Ｂを制御する。感度成分特定部７２２Ｂは、ステップＳ２７における近似処理により得られた多項式ｆ（ｄ）から高感度成分を差し引くことで、所定の層領域の真の形状を表す形状プロファイルを生成する。図８のステップＳ４では、ステップＳ２９において出力されたプロファイルが、眼底Ｅｆの形状に対応する所定の層領域（例えば網膜色素上皮層）の形状を表すデータとして用いられる。
【０２２７】
以上で、図８のステップＳ３の処理は終了である（エンド）。
【０２２８】
［第４変形例］
上記の実施形態又はその変形例では、近似処理部７２２Ａにより形状プロファイルを多項式近似することにより得られた多項式からアライメントエラーに対する低感度成分又は高感度成分を特定する場合について説明した。しかしながら、実施形態に係る眼科装置１の構成はこれに限定されるものではない。例えば、層領域特定部７２１により得られた形状プロファイルに対して円フィッティング処理を行うことにより得られた円の位置及びサイズから、アライメントエラーに対する低感度成分を特定してもよい。いくつかの実施形態では、層領域特定部７２１により得られた形状プロファイルに対して楕円フィッティング処理又は非球面フィッティング処理を行うことにより得られた楕円又は非球面の位置及びそのサイズから、アライメントエラーに対する低感度成分が特定される。例えば、非球面フィッティング処理には、コーニック定数を用いた非球面の式が用いられる。
【０２２９】
以下、第４変形例に係る眼科装置について、実施形態に係る眼科装置１との相違点を中心に説明する。
【０２３０】
第４変形例では、特定部７２２は、層領域特定部７２１によるセグメンテーション処理により得られた形状プロファイルに対して公知の円フィッティング処理を行うことで、円周の少なくとも一部が所定の層領域の形状に沿った１以上の円（又はその円周の一部）を特定する。特定部７２２は、特定された１以上の円の位置を表す成分とサイズを表す成分とをアライメントエラーに対する低感度成分として特定することが可能である。１以上の円の位置には、中心位置、曲率中心位置などがある。１以上の円のサイズには、半径、曲率などがある。
【０２３１】
例えば、特定部７２２は、形状プロファイルの所定の範囲ごとに、最小二乗法と円の中心位置及び半径を未知数とした偏微分とを用いた公知の円フィッティング処理により当該範囲の近似円を特定する。近似円は、中心位置と半径とにより特定される。特定部７２２は、上記の円フィッティング処理を繰り返すことで、形状プロファイルの所望の範囲内で２以上の近似円を特定することが可能である。特定部７２２は、特定された２以上の近似円の中心位置と半径とを用いて所定の層領域の真の形状を表す形状プロファイルとして生成する。
【０２３２】
第４変形例では、図８のステップＳ３において、次にように被検眼Ｅの眼底Ｅｆの形状の特定処理が実行される。
【０２３３】
図１３に、実施形態の第４変形例に係る眼科装置の動作の一例を示す。図１３は、図８のステップＳ３の動作例のフロー図を表す。制御部８０の記憶部には、図１３に示す処理を実現するためのコンピュータプログラムが記憶されている。制御部８０は、このコンピュータプログラムに従って動作することにより、図１３に示す処理を実行する。
【０２３４】
（Ｓ３１：アライメント光の投射を開始）
ステップＳ３の処理が開始されると、制御部８０は、ステップＳ１１と同様に、アライメント光投射部４０を制御して、被検眼Ｅに対してアライメント光の投射を開始させる。
【０２３５】
（Ｓ３２：アライメント）
制御部８０は、ステップＳ１２と同様に、移動機構９０を制御し、被検眼Ｅに対する測定部１０の位置合わせを行う。
【０２３６】
（Ｓ３３：アライメント完了？）
制御部８０は、ステップＳ１３と同様に、所定のアライメント完了条件が満足したか否かを判定する。
【０２３７】
所定のアライメント完了条件を満足していないと判定されたとき（Ｓ３３：Ｎ）、眼科装置１の動作はステップＳ３２に移行する。所定のアライメント完了条件を満足していると判定されたとき（Ｓ３３：Ｙ）、眼科装置１の動作はステップＳ３４に移行する。
【０２３８】
（Ｓ３４：ＯＣＴ計測）
ステップＳ３３において、所定のアライメント完了条件を満足していると判定されたとき（Ｓ３３：Ｙ）、制御部８０は、ステップＳ１３と同様に、ＯＣＴ部３０を制御することにより、眼底Ｅｆにおける所定部位に対してＯＣＴスキャンを実行させることでＯＣＴ計測を実行させる。
【０２３９】
（Ｓ３５：断層像を形成）
続いて、制御部８０は、ステップＳ１５と同様に、ステップＳ３４において得られたスキャンデータに基づいて被検眼Ｅの断層像を画像形成部６０に形成させる。
【０２４０】
（Ｓ３６：セグメンテーション処理）
次に、制御部８０は、ステップＳ１６と同様に、ステップＳ３５において形成された断層像に対してセグメンテーション処理を施すことにより所定の層領域（例えば、網膜色素上皮層）を層領域特定部７２１に特定させる。
【０２４１】
（Ｓ３７：円フィッティング処理）
次に、制御部８０は、特定部７２２を制御することにより、ステップＳ３６において得られた形状データに対して公知の円フィッティング処理を実行させる。いくつかの実施形態では、制御部８０は、ステップＳ３６において得られた形状データに対して光学的な歪みを補正された補正形状データに対して円フィッティング処理を実行させる。例えば、特定部７２２は、ステップＳ３６において得られたＯＣＴ座標系の形状データを実空間座標系の形状データに変換することにより補正形状データを取得することが可能である。特定部７２２は、上記のように、形状データ（又は補正形状データ）に対応した１以上の近似円を特定し、特定された１以上の近似円の中心位置及び曲率を用いて、所定の層領域の真の形状を表す形状プロファイルを生成する。図８のステップＳ４では、ステップＳ３７において出力されたプロファイルが、眼底Ｅｆの形状に対応する所定の層領域（例えば網膜色素上皮層）の形状を表すデータとして用いられる。
【０２４２】
以上で、図８のステップＳ３の処理は終了である（エンド）。
【０２４３】
［効果］
実施形態に係る眼科情報処理装置、眼科装置、眼科情報処理方法、及びプログラムについて説明する。
【０２４４】
いくつかの実施形態に係る眼科情報処理装置（制御処理部５０、又はデータ処理部７０）は、取得部（ＯＣＴ部３０及び画像形成部６０、外部装置（外部の眼科装置）、又は記録媒体からデータを受信するデバイス（通信インターフェイス、入出力インターフェイス等））と、組織特定部（層領域特定部７２１）と、特定部（７２２）とを含む。取得部は、被検眼（Ｅ）に対して光コヒーレンストモグラフィを実行するための光学系（ＯＣＴ部３０）を用いて得られたスキャンデータに基づいて形成された被検眼の断層像を取得する。組織特定部は、断層像に対してセグメンテーション処理を施すことにより被検眼の組織の形状を表す第１形状データ（形状プロファイル）を取得する。特定部は、第１形状データから被検眼に対する光学系の位置の変化に対してその変動が小さい低感度成分を特定し、特定された低感度成分に基づいて組織の形状を表す第２形状データ（形状プロファイル）を求める。
【０２４５】
このような構成によれば、被検眼の組織の形状を表す形状データから、被検眼と光学系との相対位置のずれに対する低感度成分を特定し、形状データの低感度成分から新たな形状データを生成するようにしたので、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２４６】
いくつかの実施形態では、特定部は、被検眼に対する光学系のアライメント基準位置に対し被検眼に対する光学系の位置の対称成分を低感度成分として特定する。
【０２４７】
このような構成によれば、アライメント基準位置に対する光学系の位置の対称成分を形状データから特定して新たな形状データを生成するようにしたので、簡素な処理で、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２４８】
いくつかの実施形態では、特定部は、第１形状データを多項式近似することにより得られた多項式のうち、光学系の位置の偶数次成分を低感度成分として特定する。
【０２４９】
このような構成によれば、多項式近似処理を用いた簡素な処理で、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２５０】
いくつかの実施形態では、低感度成分は、組織の曲率を表す成分を含む。
【０２５１】
このような構成によれば、曲率に基づいて被検眼の組織を特定するようにしたので、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２５２】
いくつかの実施形態では、特定部は、第１形状データに対して円フィッティング処理を行うことにより得られた円の位置を表す成分とサイズを表す成分とを低感度成分として特定する。
【０２５３】
このような構成によれば、円フィッティング処理により得られた円の位置及びサイズに基づいて被検眼の組織を特定するようにしたので、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２５４】
いくつかの実施形態に係る眼科情報処理装置（制御処理部５０、又はデータ処理部７０）は、取得部（ＯＣＴ部３０及び画像形成部６０、外部装置（外部の眼科装置）、又は記録媒体からデータを受信するデバイス（通信インターフェイス、入出力インターフェイス等））と、組織特定部（層領域特定部７２１）と、特定部（７２２）とを含む。取得部は、被検眼（Ｅ）に対して光コヒーレンストモグラフィを実行するための光学系（ＯＣＴ部３０）を用いて得られたスキャンデータに基づいて形成された被検眼の断層像を取得する。組織特定部は、断層像に対してセグメンテーション処理を施すことにより被検眼の組織の形状を表す第１形状データを取得する。特定部は、第１形状データから被検眼に対する光学系の位置の変化に対してその変動が大きい高感度成分を特定し、特定された高感度成分を第１形状データから除去することにより組織の形状を表す第２形状データを求める。
【０２５５】
このような構成によれば、被検眼の組織の形状を表す形状データから、被検眼と光学系との相対位置のずれに対する高感度成分を特定し、形状データから高感度成分を除去することにより新たな形状データを生成するようにしたので、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２５６】
いくつかの実施形態では、特定部は、被検眼に対する光学系のアライメント基準位置に対し被検眼に対する光学系の位置の非対称成分を高感度成分として特定する。
【０２５７】
このような構成によれば、アライメント基準位置に対する光学系の位置の非対称成分を形状データから特定し、形状データから非対称成分を除去することにより新たな形状データを生成するようにしたので、簡素な処理で、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２５８】
いくつかの実施形態では、特定部は、第１形状データを多項式近似することにより得られた多項式のうち、光学系の位置の奇数次成分を高感度成分として特定する。
【０２５９】
このような構成によれば、多項式近似処理を用いた簡素な処理で、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２６０】
いくつかの実施形態では、高感度成分は、所定の基準方向に対する組織の傾きを表す成分を含む。
【０２６１】
このような構成によれば、組織の傾きに基づいて被検眼の組織を特定するようにしたので、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２６２】
いくつかの実施形態では、組織は、眼底（Ｅｆ）における所定の層領域（例えば、網膜色素上皮層）を含む。
【０２６３】
このような構成によれば、被検眼と光学系との変位の影響を受けることなく眼底の形状を高い再現性で高精度に特定することができるようになる。
【０２６４】
いくつかの実施形態は、被検眼を他覚的に測定することにより得られた屈折度数と、特定部により取得された第２形状データに基づく組織の形状に対応した被検眼の光学特性を表すパラメータとに基づいて、被検眼の中心窩を含む領域の周辺領域の屈折度数を算出する算出部（周辺屈折度数算出部７３Ｃ）を含む。
【０２６５】
このような構成によれば、被検眼の眼底の形状に対応して、中心窩を含む領域の周辺領域の屈折度数を高精度に取得することが可能になる。
【０２６６】
いくつかの実施形態では、組織の形状は、所定の基準方向に対する組織のチルト角度を含む。
【０２６７】
このような構成によれば、所定の基準方向に対する眼底における所定の層領域のチルト角度に対応して、中心窩を含む領域の周辺領域の屈折度数を高精度に取得することが可能になる。
【０２６８】
いくつかの実施形態に係る眼科装置（１）は、光学系（ＯＣＴ部３０に含まれる光学系）と、被検眼と光学系とを相対的に移動する移動機構（９０）と、上記のいずれかに記載の眼科情報処理装置と、を含む。
【０２６９】
このような構成によれば、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することが可能な眼科装置を提供することができるようになる。
【０２７０】
いくつかの実施形態に係る眼科情報処理方法は、取得ステップと、組織特定ステップと、特定ステップとを含む。取得ステップは、被検眼（Ｅ）に対して光コヒーレンストモグラフィを実行するための光学系（ＯＣＴ部３０）を用いて得られたスキャンデータに基づいて形成された被検眼の断層像を取得する。組織特定ステップは、断層像に対してセグメンテーション処理を施すことにより被検眼の組織の形状を表す第１形状データを取得する。特定ステップは、第１形状データから被検眼に対する光学系の位置の変化に対してその変動が小さい低感度成分を特定し、特定された低感度成分に基づいて組織の形状を表す第２形状データを求める。
【０２７１】
このような方法によれば、被検眼の組織の形状を表す形状データから、被検眼と光学系との相対位置のずれに対する低感度成分を特定し、形状データの低感度成分から新たな形状データを生成するようにしたので、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２７２】
いくつかの実施形態に係る眼科情報処理方法は、取得ステップと、組織特定ステップと、特定ステップとを含む。取得ステップは、被検眼（Ｅ）に対して光コヒーレンストモグラフィを実行するための光学系（ＯＣＴ部３０）を用いて得られたスキャンデータに基づいて形成された被検眼の断層像を取得する。組織特定ステップは、断層像に対してセグメンテーション処理を施すことにより被検眼の組織の形状を表す第１形状データを取得する。特定ステップは、第１形状データから被検眼に対する光学系の位置の変化に対してその変動が大きい高感度成分を特定し、特定された高感度成分を第１形状データから除去することにより組織の形状を表す第２形状データを求める。
【０２７３】
このような方法によれば、被検眼の組織の形状を表す形状データから、被検眼と光学系との相対位置のずれに対する高感度成分を特定し、形状データから高感度成分を除去することにより新たな形状データを生成するようにしたので、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２７４】
いくつかの実施形態は、被検眼を他覚的に測定することにより得られた屈折度数と、特定ステップにおいて取得された第２形状データに基づく組織の形状に対応した被検眼の光学特性を表すパラメータとに基づいて、被検眼の中心窩を含む領域の周辺領域の屈折度数を算出する算出ステップを含む。
【０２７５】
このような方法によれば、被検眼の眼底の形状に対応して、中心窩を含む領域の周辺領域の屈折度数を高精度に取得することが可能になる。
【０２７６】
いくつかの実施形態に係るプログラムは、コンピュータに、上記のいずれかに記載の眼科情報処理方法の各ステップを実行させる。
【０２７７】
このようなプログラムによれば、被検眼の組織の形状を表す形状データから、被検眼と光学系との相対位置のずれに対する低感度成分又は高感度成分を特定し、相対位置のずれに対する影響が小さくなるように新たな形状データを生成するようにしたので、被検眼と組織の形状を計測するための光学系との変位の影響を低減しつつ、被検眼の組織の形状を高い再現性で高精度に特定することができるようになる。
【０２７８】
＜その他＞
以上に示された実施形態は、この発明を実施するための一例に過ぎない。この発明を実施しようとする者は、この発明の要旨の範囲内において任意の変形、省略、追加等を施すことが可能である。
【０２７９】
いくつかの実施形態では、上記の眼科情報処理方法だけでなく上記の眼科装置を制御する制御方法をコンピュータに実行させるためのプログラムが提供される。このようなプログラムを、コンピュータによって読み取り可能な任意の記録媒体に記憶させることができる。この記録媒体としては、たとえば、半導体メモリ、光ディスク、光磁気ディスク（ＣＤ−ＲＯＭ／ＤＶＤ−ＲＡＭ／ＤＶＤ−ＲＯＭ／ＭＯ等）、磁気記憶媒体（ハードディスク／フロッピー（登録商標）ディスク／ＺＩＰ等）などを用いることが可能である。また、インターネットやＬＡＮ等のネットワークを通じてこのプログラムを送受信することも可能である。
【０２８０】
本発明のいくつかの実施形態を説明したが、これら実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると共に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０２８１】
１  眼科装置
１０  測定部
２０  屈折測定部
３０  ＯＣＴ部
４０  アライメント光投射部
５０  制御処理部
６０  画像形成部
７０  データ処理部
７１  アライメント処理部
７２  解析部
７２１  層領域特定部
７２２  特定部
７２２Ａ  近似処理部
７２２Ｂ  感度成分特定部
７３  算出部
８０  制御部
９０  移動機構
１００  撮影部
ＢＳ１、ＢＳ２  ビームスプリッタ
Ｅ  被検眼
Ｅｃ  角膜
Ｅｆ  眼底
【図１】
2020110579.tif 000014
【図２】
2020110579.tif 000015
【図３】
2020110579.tif 000016
【図４】
2020110579.tif 000017
【図５】
2020110579.tif 000018
【図６】
2020110579.tif 000019
【図７】
2020110579.tif 000020
【図８】
2020110579.tif 000021
【図９】
2020110579.tif 000022
【図１０】
2020110579.tif 000023
【図１１】
2020110579.tif 000024
【図１２】
2020110579.tif 000025
【図１３】
2020110579.tif 000026