(57)【要約】【課題】  適切な眼内レンズを選択できる眼内レンズ選択装置、および眼内レンズ選択プログラムを提供することを技術課題とする。【解決手段】  眼内レンズを選択するための眼内レンズ選択装置であって、被検眼の眼形状パラメータを取得する取得手段と、制御手段と、を備え、前記制御手段は、眼内レンズ度数と眼内レンズ形状との関係に基づくＩＯＬオフセットと、前記眼形状パラメータと、を用いて術後予測屈折力を算出することを特徴とする。【選択図】      図６
2020049128.tif 000002
【特許請求の範囲】
【請求項１】
眼内レンズを選択するための眼内レンズ選択装置であって、
被検眼の眼形状パラメータを取得する取得手段と、
制御手段と、を備え、
前記制御手段は、眼内レンズ度数と眼内レンズ形状との関係に基づくＩＯＬオフセットと、前記眼形状パラメータと、を用いて術後予測屈折力を算出することを特徴とする眼内レンズ選択装置。
【請求項２】
前記制御手段は、眼内レンズ度数に対する前記眼内レンズの支持部と主点位置との位置関係の変化を前記ＩＯＬオフセットによって補正することを特徴とする請求項１の眼内レンズ選択装置。
【請求項３】
前記制御手段は、前記眼内レンズ度数と前記ＩＯＬオフセットとの関係式に、前記眼内レンズ度数を入力することによって前記ＩＯＬオフセットを算出することを特徴とする請求項１または２の眼内レンズ選択装置。
【請求項４】
前記関係式は、前記眼内レンズ度数と前記ＩＯＬオフセットのデータに基づく近似式であることを特徴とする請求項３の眼内レンズ選択装置。
【請求項５】
前記制御手段は、近軸光線追跡式を用いて前記術後予測屈折力を算出することを特徴とする請求項１〜４のいずれかの眼内レンズ選択装置。
【請求項６】
前記眼形状パラメータは、眼軸長、術後予測前房深度、角膜前面曲率半径、角膜後面曲率半径、角膜中心厚であることを特徴とする請求項１〜５のいずれかの眼内レンズ選択装置。
【請求項７】
眼内レンズを選択するための眼内レンズ選択装置によって実行される眼内レンズ選択プログラムであって、前記眼内レンズ選択装置のプロセッサによって実行されることで、
被検眼の眼形状パラメータを取得する取得ステップと、
眼内レンズ度数と眼内レンズ形状との関係に基づくＩＯＬオフセットと、前記眼形状パラメータと、を用いて術後予測屈折力を算出する算出ステップと、を前記眼内レンズ選択装置に実行させることを特徴とする眼内レンズ選択プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、眼内レンズの度数を選択するための眼内レンズ選択装置、および眼内レンズ選択プログラムに関する。
【背景技術】
【０００２】
白内障手術において、水晶体核の除去後に被検眼の眼内に挿入される眼内レンズ（以下、ＩＯＬ）のパワー（以下、度数）を決定する眼内レンズ度数決定装置が知られている（特許文献１参照）。このような装置では、眼内レンズ度数を決定するために術後予測前房深度（眼内レンズの位置）の推定を行い、眼内レンズの度数計算を行っていた。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−９４４１０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、従来の眼内レンズ計算式において、眼内レンズの度数に応じた形状の変化を考慮することができていなかった。
【０００５】
本開示は、従来の問題点に鑑み、適切な眼内レンズを選択できる眼内レンズ選択装置、および眼内レンズ選択プログラムを提供することを技術課題とする。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、本開示は以下のような構成を備えることを特徴とする。
【０００７】
（１）  眼内レンズを選択するための眼内レンズ選択装置であって、被検眼の眼形状パラメータを取得する取得手段と、制御手段と、を備え、前記制御手段は、眼内レンズ度数と眼内レンズ形状との関係に基づくＩＯＬオフセットと、前記眼形状パラメータと、を用いて術後予測屈折力を算出することを特徴とする。
（２）  眼内レンズを選択するための眼内レンズ選択装置によって実行される眼内レンズ選択プログラムであって、前記眼内レンズ選択装置のプロセッサによって実行されることで、被検眼の眼形状パラメータを取得する取得ステップと、眼内レンズ度数と眼内レンズ形状との関係に基づくＩＯＬオフセットと、前記眼形状パラメータと、を用いて術後予測屈折力を算出する算出ステップと、を前記眼内レンズ選択装置に実行させることを特徴とする。
【図面の簡単な説明】
【０００８】
【図１】本実施例に係る眼科撮影装置の構成について説明する概略構成図である。
【図２】撮像された前眼部像が表示された前眼部観察画面を示す図である。
【図３】前眼部断面画像の一例を示す図である。
【図４】ＩＯＬ度数の算出について説明するための図である。
【図５】オフセットの算出方法を説明する図である。
【図６】近軸光線追跡を用いた術後残余屈折力の算出について説明する図である。
【図７】眼内レンズの形状を示す図である。
【図８】眼内レンズの度数とＩＯＬオフセットとの関係を示す図である。
【発明を実施するための形態】
【０００９】
＜実施形態＞
以下、本開示に係る実施形態について説明する。本実施形態の眼内レンズ選択装置（例えば、眼科撮影装置１）は、被検眼に挿入する眼内レンズの度数を決定するための装置である。眼内レンズ選択装置は、例えば、取得部（例えば、ＯＣＴ光学系１００または制御部８０）と、制御部（例えば、制御部８０）を備える。取得部は、被検眼の眼形状パラメータを取得する。制御部は、眼内レンズ度数と眼内レンズ形状との関係に基づくＩＯＬオフセットと、眼形状パラメータと、を用いて術後予測屈折力を算出する。これによって、眼内レンズ選択装置は、度数毎に眼内レンズの形状が変化しても、より患者に適した眼内レンズを選択することができる。
【００１０】
なお、制御部は、眼内レンズ度数に対する眼内レンズの支持部と主点位置との位置関係の変化をＩＯＬオフセットによって補正してもよい。これによって、度数毎に支持部と主点位置との位置関係が変化しても、患者に適した眼内レンズを選択することができる。
【００１１】
なお、制御部は、眼内レンズ度数とＩＯＬオフセットとの関係式に、眼内レンズ度数を入力することによってＩＯＬオフセットを算出してもよい。関係式は、例えば、眼内レンズ度数とＩＯＬオフセットのデータに基づく近似式であってもよい。これによって、眼内レンズの設計データが取得できない場合であっても、眼内レンズ度数に適したＩＯＬオフセットを設定できる。
【００１２】
なお、眼形状パラメータは、眼軸長、術後予測前房深度、角膜前面曲率半径、角膜後面曲率半径、角膜中心厚であってもよい。角膜後面曲率半径は、実測値であってもよいし、推定値であってもよい。
【００１３】
なお、眼内レンズ選択装置は、眼内レンズ度数の計算結果から眼内レンズを選択してもよいし、選択肢を提示することによって、眼内レンズの選択を支援してもよい。
【００１４】
なお、制御部は、記憶部に記憶された眼内レンズ選択プログラムを実行してもよい。眼内レンズ選択プログラムは、例えば、取得ステップと、算出ステップを含む。取得ステップは、例えば、被検眼の眼形状パラメータを取得するステップである。算出ステップは、眼内レンズ度数と眼内レンズ形状との関係に基づくＩＯＬオフセットと、眼形状パラメータと、を用いて術後予測屈折力を算出するステップである。
【００１５】
＜実施例＞
以下、本開示に係る実施例を図面に基づいて説明する。本実施例の眼科撮影装置１の光学系は、図示無き筐体に内蔵されている。また、その筐体は、周知のアライメント移動機構の駆動により、被検眼Ｅに対して３次元的に移動される。なお、以下の説明においては、被検眼Ｅの光軸方向をＺ方向、水平方向をＸ方向、鉛直方向をＹ方向として説明する。眼底の表面方向をＸＹ方向として考えてもよい。
【００１６】
眼科撮影装置１は、例えば、干渉光学系（ＯＣＴ光学系）と、角膜測定光学系と、を備える。
【００１７】
干渉光学系（ＯＣＴ光学系）１００を備えている。被検眼Ｅの断面像を撮影するための撮像デバイスとして用いられる。被検眼Ｅの眼軸長を測定するために用いられてもよい。ＯＣＴ光学系１００は、被検眼Ｅに測定光を照射する。ＯＣＴ光学系１００は、前眼部（例えば、角膜、水晶体など）から反射された測定光と，参照光との干渉状態を受光素子（検出器１２０）によって検出する。ＯＣＴ光学系１００は、前眼部上の撮像位置を変更するため、前眼部における測定光の照射位置を変更する照射位置変更ユニット（例えば、光スキャナ１０８）を備える。制御部８０は、設定された撮像位置情報に基づいて照射位置変更ユニットの動作を制御し、検出器１２０からの受光信号に基づいて断層画像を取得する。
【００１８】
ＯＣＴ光学系１００は、いわゆる眼科用光断層干渉計（ＯＣＴ：Optical coherence tomography）の装置構成を持つ。ＯＣＴ光学系１００は、光源１０２から出射された光をカップラー（光分割器）１０４によって測定光（試料光）と参照光に分割する。そして、ＯＣＴ光学系１００は、測定光学系１０６によって測定光を前眼部に導き、参照光を参照光学系１１０に導く。その後、前眼部によって反射された測定光と，参照光との合成による干渉光を検出器（受光素子）１２０に受光させる。
【００１９】
光源１０２から出射された光は、カップラー１０４によって測定光束と参照光束に分割される。そして、測定光束は、光ファイバーを通過した後、空気中へ出射される。その光束は、光スキャナ１０８、及び測定光学系１０６の他の光学部材を介して前眼部に集光される。そして、前眼部で反射された光は、同様の光路を経て光ファイバーに戻される。
【００２０】
光スキャナ１０８は、被検眼Ｅ上でＸＹ方向（横断方向）に測定光を走査させる。光スキャナ１０８は、例えば、２つのガルバノミラーであり、その反射角度が駆動機構１０９によって任意に調整される。
【００２１】
これによって、光源１０２から出射された光束はその反射（進行）方向が変化され、被検眼Ｅ上で任意の方向に走査される。これにより、前眼部上における撮像位置が変更される。光スキャナ１０８としては、光を偏向させる構成であればよい。例えば、反射ミラー（ガルバノミラー、ポリゴンミラー、レゾナントスキャナ）の他、光の進行（偏向）方向を変化させる音響光学素子（ＡＯＭ）等が用いられる。
【００２２】
参照光学系１１０は、被検眼Ｅでの測定光の反射によって取得される反射光と合成される参照光を生成する。参照光学系１１０は、マイケルソンタイプであってもよいし、マッハツェンダタイプであっても良い。参照光学系１１０は、例えば、反射光学系（例えば、参照ミラー）によって形成され、カップラー１０４からの光を反射光学系により反射することにより再度カップラー１０４に戻し、検出器１２０に導く。他の例としては、参照光学系１１０は、透過光学系（例えば、光ファイバー）によって形成され、カップラー１０４からの光を戻さず透過させることにより検出器１２０へと導く。
【００２３】
参照光学系１１０は、参照光路中の光学部材を移動させることにより、測定光と参照光との光路長差を変更する構成を有する。例えば、参照ミラーが光軸方向に移動される。光路長差を変更するための構成は、測定光学系１０６の測定光路中に配置されてもよい。
【００２４】
検出器１２０は、測定光と参照光との干渉状態を検出する。フーリエドメインＯＣＴの場合では、干渉光のスペクトル強度が検出器１２０によって検出され、スペクトル強度データに対するフーリエ変換によって所定範囲における深さプロファイル（Ａスキャン信号）が取得される。ここで、制御部８０は、光スキャナ１０８により測定光を前眼部上で所定の横断方向に走査することにより断層像を取得できる。すなわち、被検眼の前眼部断層像を撮像する。例えば、Ｘ方向もしくはＹ方向に走査することにより、被検眼前眼部のＸＺ面もしくはＹＺ面における断層像（前眼部断層像）を取得できる（なお、本実施形態においては、このように測定光を前眼部に対して一次元走査し、断層像を得る方式をＢスキャンとする）。なお、取得された前眼部断層像は、制御部８０に接続されたメモリ８５に記憶される。さらに、測定光をＸＹ方向に二次元的に走査することにより、被検眼前眼部の三次元画像を取得することも可能である。
【００２５】
例えば、フーリエドメインＯＣＴとしては、Spectral-domain OCT（ＳＤ−ＯＣＴ）、Swept-source OCT（ＳＳ−ＯＣＴ）が挙げられる。また、Time-domain OCT（ＴＤ−ＯＣＴ）であってもよい。
【００２６】
ＳＤ−ＯＣＴの場合、光源１０２として低コヒーレント光源（広帯域光源）が用いられ、検出器１２０には、干渉光を各周波数成分（各波長成分）に分光する分光光学系（スペクトルメータ）が設けられる。スペクトルメータは、例えば、回折格子とラインセンサからなる。
【００２７】
ＳＳ−ＯＣＴの場合、光源１０２として出射波長を時間的に高速で変化させる波長走査型光源（波長可変光源）が用いられ、検出器１２０として、例えば、単一の受光素子が設けられる。光源１０２は、例えば、光源、ファイバーリング共振器、及び波長選択フィルタによって構成される。そして、波長選択フィルタとして、例えば、回折格子とポリゴンミラーの組み合わせ、ファブリー・ペローエタロンを用いたものが挙げられる。
【００２８】
角膜測定光学系３００は、角膜形状を測定するために用いられる。角膜測定光学系３００は、投影光学系５０、アライメント光学系４０、前眼部撮像光学系３０、に大別される。
【００２９】
投影光学系５０は、光軸Ｌ１を中心に配置されたリング状の光源５１を有し、角膜にリング指標を投影して角膜形状（曲率、乱視軸角度、等）を測定するために用いられる。なお、光源５１には、例えば、赤外光または可視光を発するＬＥＤが使用される。なお、投影光学系５０について、光軸Ｌ１を中心とする同一円周上に少なくとも３つ以上の点光源が配置されていればよく、間欠的なリング光源であってもよい。さらに、複数のリング指標を投影するプラチド指標投影光学系であってもよい。
【００３０】
アライメント光学系４０は、光源５１の内側に配置され、赤外光を発する光源４１（例えば、λ＝９７０ｎｍ）を有し、角膜Ｅｃにアライメント指標を投影するために用いられる。そして、角膜Ｅｃに投影されたアライメント指標は、被検眼に対する位置合わせ（例えば、自動アライメント、アライメント検出、手動アライメント、等）に用いられる。本実施形態において、投影光学系５０は、角膜Ｅｃに対してリング指標を投影する光学系であって、リング指標は、マイヤーリングも兼用する。また、投影光学系４０の光源４１は、前眼部を斜め方向から赤外光にて照明する前眼部照明を兼用する。なお、投影光学系４０において、さらに、角膜Ｅｃに平行光を投影する光学系を設け、投影光学系４０による有限光との組合せにより前後のアライメントを行うようにしてもよい。
【００３１】
前眼部撮像光学系３０は、前眼部正面像を撮像（取得）するために用いられる。前眼部撮像光学系３０は、ダイクロイックミラー３３、対物レンズ４７、ダイクロイックミラー６２、フィルタ３４、撮像レンズ３７、二次元撮像素子３５、を含み、被検眼の前眼部正面像を撮像するために用いられる。二次元撮像素子３５は、被検眼の前眼部と略共役な位置に配置されている。
【００３２】
前述の投影光学系４０、投影光学系５０による前眼部反射光は、ダイクロイックミラー３３、対物レンズ４７、ダイクロイックミラー６２、フィルタ３４、及び撮像レンズ３７を介して二次元撮像素子３５に結像される。
【００３３】
光源１０は、固視灯である。また、例えば、光源１０から発せられた光の前眼部での反射により取得される前眼部反射光の一部は、ダイクロイックミラー３３で反射され、前眼部撮像光学系３０で結像される。
【００３４】
次に、制御系について説明する。制御部８０は、装置全体の制御及び測定結果の算出を行う。制御部８０は、ＯＣＴ光学系１００の各部材、角膜測定光学系３００の各部材、表示部７０、操作部８４、メモリ８５、等と接続されている。操作部８４には、操作入力部として、マウス等の汎用インターフェースが用いられてもよいし、その他、タッチパネルが用いられてもよい。メモリ８５には、各種制御プログラムの他、制御部８０が前眼部画像解析を行う解析プログラム等が記憶されている。
【００３５】
＜制御動作＞
以上のような構成を備える装置において、眼内レンズ度数を決定するときの制御動作について説明する。検者は、表示部７０に表示される被検眼のアライメント状態を見ながら、図示なきジョイスティック等の操作部を用いて、装置を上下左右及び前後方向に移動させ、装置を被検眼Ｅに対して所定の位置関係に置く。この場合、検者は、固視標を被検眼に固視させる。
【００３６】
アライメントの際には、光源４１及び光源５１が点灯される。ここで、検者は、図２に示すように、表示部７０に電子的に表示されたレチクルＬＣと、光源４１によるリング指標Ｑ１，Ｑ２が同心円状になるように上下左右のアライメントを行う。これによって、被検眼の角膜頂点に本装置の光軸Ｌ１が通るようにＸＹ方向にアライメントされる。また、検者は、リング指標Ｑ１のピントが合うように、前後のアライメントを行う。
【００３７】
前眼部に対するアライメントが完了されると、制御部８０は、前眼部撮像光学系３０によって被検眼の前眼部を撮影する。また、制御部８０は、予め設定された走査パターンに基づき、ＯＣＴ光学系１００によって被検眼の断面画像５００を撮影する（図３参照）。取得された前眼部画像および断面画像は、メモリ８５等に記憶される。
【００３８】
制御部８０は、メモリ８５に記憶された前眼部画像４００におけるリング指標像Ｑ１，Ｑ２に基づいて被検眼の角膜形状をそれぞれ算出する。角膜形状とは、例えば、強主経線方向及び弱主経線方向における角膜前面の角膜曲率半径、角膜の乱視軸角度等である。例えば、制御部８０は、リング指標像Ｑ１，Ｑ２の大きさおよび形状に基づいて角膜形状を算出する。
【００３９】
また、制御部８０は、ＯＣＴ光学系１００を用いて撮影された断面画像を解析する。例えば、制御部８０は、断面画像のエッジ検出によって角膜、水晶体などの位置を検出し、その位置に基づいて角膜厚、前房深度、水晶体厚を測定する。また、制御部８０は、検出された角膜および水晶体の前面または後面を円近似（または楕円近似、コニック曲線近似等）し、この近似曲線に基づいて角膜後面の曲率半径、水晶体前面曲率、水晶体後面曲率等を測定する。
【００４０】
測定および画像解析が完了すると、制御部８０は、既知であるＳＲＫ／Ｔ式、Ｂｉｎｋｈｏｒｓｔ式等を一部流用して眼内レンズ度数を算出する。例えば、ＳＲＫ／Ｔ式、Ｂｉｎｋｈｏｒｓｔ式等に上記の測定データを代入する。ＳＲＫ／Ｔ式（下記の式（１））を用いる場合、角膜曲率半径、眼軸長、術後予測前房深度（詳細は後述する）等を用いて眼内レンズ度数が計算される。
【００４１】
【数１】
2020049128.tif 000003
ここで、R:角膜曲率半径[mm](R＝(nk−1.000)×1000/K)、nk:検者によって選択された屈折率、LO:AL＋RT[mm]、RT:網膜の厚み[mm](RT＝0.65696−0.02029×AL)、AL:眼軸長[mm]、AD´:術後予測前房深度の補正値[mm](AD´＝H＋OF,OF＝AD−3.336)、AD:術後予測前房深度[mm](AD＝0.62467×A−68.747)、A:A定数、H:角膜高さ[mm](H=R−(R×R−((Cw×Cw)/4))1/2)（ただし、(R×R−((Cw×Cw)/4))＜0の場合、H＝R）、Cw:角膜幅[mm]、Cw＝−5.41＋0.58412×LC＋0.098×K、LC:眼軸長の補正値[mm]（AL≦24.2の場合LC＝AL、AL＞24.2の場合LC＝−3.446＋1.716×AL−0.0237×AL2）、DR:術後希望する矯正用レンズの屈折力[D]、LP:移植するIOLの度数[D]、V:頂点距離、na:房水および硝子体の屈折率(＝1.336)、nc:角膜の屈折率(＝1.333)、ncml: nc−1(＝0.333)である。
【００４２】
＜術後予測前房深度の算出＞
続いて、本明細書にて定義する術後予測前房深度ＥＬＰpredの算出について説明する。術後予測前房深度ＥＬＰpredは、図４に示すように、角膜頂点Ｋを基準として、中心角膜厚ＣＣＴ、前房深度ＡＣＤ、補正量Ｘ、補正量αとが足し合わされることによって算出される。したがって、術後予測前房深度ＥＬＰpredは次式（２）のように表すことができる。
【数２】
2020049128.tif 000004
ここで、補正量Ｘは、水晶体前面の位置から赤道位置（水晶体の最大径部分）ＥＰＰまでの距離である。本実施例において、赤道位置は、水晶体前面と水晶体後面とが交差する位置とされる。補正量αは、水晶体赤道位置からＩＯＬまでの距離である。補正量αは、ＩＯＬが水晶体嚢から圧力を受けることによる水晶体後嚢側または前房側への移動の影響を受ける。術後予測前房深度としては、角膜裏面から定義される場合もあるが、ここでは、角膜前面からＩＯＬまでの距離とした。
【００４３】
次に、補正量Ｘを算出する方法について説明する。図５において、距離（水晶体赤道半径）ｈは、光軸Ｌ１から水晶体前面の近似円と水晶体後面の近似円との交点までの距離を示している。距離Ｘ１は、光軸Ｌ１における水晶体前面曲率中心Ｏ４から水晶体後面までの距離を示している。距離Ｘ１´は、水晶体前面の近似円と水晶体後面の近似円との交点から水晶体後面までの距離を示している。距離Ｘ２は、光軸Ｌ１における水晶体後面曲率中心Ｏ３から水晶体前面までの距離を示している。水晶体前面曲率半径Ｒ３、水晶体後面曲率半径Ｒ４、水晶体厚ＬＴを用いると、ピタゴラスの定理により以下の式（３）が成り立つ。
【００４４】
【数３】
2020049128.tif 000005
そして、上記の式（３）において、距離ｈが同様であるため、これらの式を補正量Ｘについて解くと以下の式（４）が成り立つ。
【００４５】
【数４】
2020049128.tif 000006
【００４６】
したがって、制御部８０は、式（４）の水晶体前面曲率半径Ｒ３、水晶体後面曲率半径Ｒ４、水晶体厚ＬＴに各数値を代入することによって補正量Ｘを算出する。各数値は、断層画像の解析によって求められる。
【００４７】
＜補正量の算出＞
続いて、補正量αを算出する方法について説明する。制御部８０は、前眼部の形状を示す前眼部パラメータに基づいて補正量αを求める。前眼部パラメータは、例えば、水晶体形状に関するパラメータである。前眼部パラメータは、例えば、水晶体厚、水晶体前面曲率、水晶体後面曲率、水晶体赤道位置、水晶体赤道径、前房深度、角膜厚、角膜径等の少なくともいずれかである。
【００４８】
例えば、制御部８０は、補正量αと前眼部パラメータとの関係式を用いて補正量αを算出する。この関係式は、例えば、過去に眼内レンズを挿入した患者の臨床データを用いて取得される。例えば、制御部８０は、眼内レンズを挿入した患者の断面画像を解析することで術後前房深度ＥＬＰの実測値を求め、術後前房深度ＥＬＰの実測値から補正量αの実測値を求める。この場合、術前に測定された角膜頂点Ｋから赤道位置ＥＰＰまでの距離（ＣＣＴ＋ＡＣＤ＋Ｘ）を、術後前房深度ＥＬＰの実測値から差し引くことによって補正量αの実測値を算出する。そして、実測された補正量αと前眼部形状パラメータとで回帰分析を行い、その結果から関係式を導き出す。本実施例では重回帰分析を用いる。重回帰分析は、例えば、１つの変数を複数の変数で予測する分析方法である。重回帰分析は、一般的な統計ソフトウェア等によって実行可能である。複数の臨床データを用いて重回帰分析を行うことによって、次式（５）のような補正量αの回帰式を得ることができる。式（５）は補正量αの回帰式の一例である。
【００４９】
【数５】
2020049128.tif 000007
【００５０】
式（５）の例では、水晶体前面曲率Ｒ３、水晶体厚に対する補正量Ｘの割合（Ｘ／ＬＴ）、水晶体直径（水晶体赤道径）に対する水晶体厚の割合（ＬＴ／ＤＩＡ）が前眼部パラメータとして用いられている。もちろん、回帰式にはその他の前眼部パラメータが用いられてもよい。このように複数のパラメータの組み合わせを前眼部パラメータとして用いてもよい。
【００５１】
例えば、制御部８０は、式（５）水晶体前面曲率Ｒ３、補正量Ｘ、水晶体厚ＬＴ、水晶体直径ＤＩＡを代入することによって、補正量αを算出する。次いで、制御部８０は、式（４），（５）によって算出した補正量Ｘと補正量α、断層画像の解析によって得られた中心角膜厚ＣＣＴと前房深度ＡＣＤを式（２）に代入することによって、術後予測前房深度ＥＬＰpredを算出する。制御部８０は、このようにして推定した術後予測前房深度ＥＬＰpredを、例えば、式（１）のAD´に代入することによって、眼内レンズ度数を計算する。
【００５２】
制御部８０は、算出した眼内レンズ度数を表示部７０に表示させる。検者は、表示部７０に表示された眼内レンズ度数を確認し、眼内レンズを選択する。例えば、眼内レンズが０．５Ｄ刻みで用意されている場合、検者は表示された眼内レンズ度数に近い眼内レンズを選択する。検者は、選択した眼内レンズの度数を操作部８４によって入力する。制御部８０は、選択された眼内レンズの度数に基づいて、被検眼の術後予測屈折力を算出する。術後予測屈折力は、眼内レンズを挿入した後の被検眼の屈折力である。
【００５３】
＜術後予測屈折力の算出＞
術後予測屈折力の算出について説明する。制御部８０は、例えば、図６に示すように、近軸光線追跡式を用いて術後予測屈折力を算出する。まず、制御部８０は、バックフォーカスｂｆの逆数である残余屈折力Ｐ０を求める。バックフォーカスｂｆは、角膜前面から焦点までの距離である。制御部８０は、ＩＯＬ屈折力p1、角膜後面屈折力p2、角膜前面屈折力p3、距離d0、距離d1、距離d2、房水屈折率n0、房水屈折率n1、角膜屈折率n2、空気屈折率n3、角膜頂点距離ＶＤなどのパラメータを用いてバックフォーカスｂｆを求める。距離d0は次式（６）、距離d1は次式（７）で表される。距離d2は角膜中心厚（CCT）である。
【００５４】
【数６】
2020049128.tif 000008
【００５５】
【数７】
2020049128.tif 000009
【００５６】
式（６），（７）において、ＡＬは眼軸長、ＥＬＰpredは術後予測前房深度、βはＩＯＬオフセットである。ＩＯＬオフセットβについては後述する。制御部８０は、次の近軸光線追跡式（８）に上記パラメータを用いることで換算傾角a１〜a３、光線高さh１〜h３を計算し、バックフォーカスｂfを求める。
【００５７】
【数８】
2020049128.tif 000010
【００５８】
次いで制御部は、バックフォーカスｂｆの逆数である残余屈折力Ｐ０を算出する（式（９））。
【数９】
2020049128.tif 000011
【００５９】
残余屈折力Ｐ０はコンタクトレンズ度数に相当する。これをメガネ屈折力Pに換算するには次式（１０）を用いる。
【数１０】
2020049128.tif 000012
【００６０】
通常この式（１０）が術後残余屈折力（術後予測屈折力）を表す式となる。術後残余屈折力Ｐは眼鏡度数に相当する。
【００６１】
＜ＩＯＬオフセット＞
ＩＯＬオフセットβは、眼内レンズの度数毎に変化する支持部と主点位置との位置関係を補正する数値である。例えば、図７に示すように、眼内レンズの支持部Ｓｐの先端から主点位置Ｈ（または主点位置Ｈ´）までの距離Ｗは、同じ種類の眼内レンズであっても度数によって異なる場合がある。したがって、制御部８０は、術後残余屈折力を算出する際に、ＩＯＬオフセットβを用いて計算上の眼内レンズの位置を補正する。
【００６２】
図８は、計算によって求められた術後残余屈折力Ｐが、患者の実際の術後屈折力と同等になるようなＩＯＬオフセットβを眼内レンズ度数ごとにプロットしたグラフである。図８に示すように、眼内レンズ度数とＩＯＬオフセットβには相関がある。したがって、制御部８０は、この相関関係を用いてＩＯＬオフセットβを決定する。例えば、制御部８０は、図８のような臨床データの回帰分析によって眼内レンズ度数とＩＯＬオフセットβとの関係を関数化し、得られた関係式に選択した眼内レンズ度数を入力することによって、ＩＯＬオフセットβを算出する。
【００６３】
制御部８０は、ＩＯＬオフセットβを算出すると、式（６）〜（１０）によって術後残余屈折力を算出し、表示部７０に表示する。検者は、表示された術後予測屈折力を確認することによって、選択した眼内レンズが適切かどうか判断する。例えば、検者は、過矯正にならないように術後残余屈折力がマイナス（近視）側になる眼内レンズを選択する。
【００６４】
以上のように、本実施例の眼内レンズ選択装置は、ＩＯＬオフセットを用いることによって、眼内レンズの度数毎に支持部と主点位置との位置関係が変化する場合であっても、より正確な術後残余屈折力を算出することができる。これによって、検者は、より患者に適した眼内レンズを選択することができる。
【００６５】
なお、以上の実施例において、制御部８０は、臨床データの回帰分析によってＩＯＬオフセットを求めたがこれに限らない。制御部８０は、眼内レンズの設計データに基づいてＩＯＬオフセットを算出してもよい。例えば、制御部は、眼内レンズの支持部から主点位置までの距離、眼内レンズ中心厚、眼内レンズ屈折率に基づいてＩＯＬオフセットを算出してもよい。
【００６６】
なお、以上の実施例において、赤道位置は水晶体前面と水晶体後面とが交差する位置としたが、他の算出方法によって赤道位置を推定してもよい。
【００６７】
なお、本実施例は、前眼部断層像（断面像）撮影用の光コヒーレンストモグラフィーデバイスにおいて、複数の走査位置にて前眼部断層像を取得することにより３次元形状画像を取得する場合においても適用可能である。例えば、ＯＣＴ光学系１００は、前眼部の３次元断面像（３次元前眼部データ）を取得する前眼部撮像デバイスであって、制御部８０は、前眼部撮像デバイスによって取得された３次元断面像に基づいて水晶体前面からチン小帯と水晶体の接点までのオフセット距離を３次元的に求めてもよい。この場合、３次元前眼部データにおける経線方向毎の水晶体前面曲率と水晶体後面曲率の平均を算出し、これに基づいてＥＬＰを算出してもよい。
【００６８】
なお、本実施形態においては、前眼部断面像を撮像する前眼部撮像デバイスとして、前眼部断層像（断面像）撮影用の光コヒーレンストモグラフィーデバイスを例に挙げたがこれに限定されない。光源からの出射光を被検眼前眼部に向けて投光し、前眼部上に光切断面を形成させる投光光学系と、光切断面の前眼部での散乱により取得される前眼部散乱光を含む光を受光する検出器を有する受光光学系と、を有し、検出器からの検出信号に基づいて前眼部断面画像を形成する構成であればよい。すなわち、検眼の前眼部にスリット光を投影し、シャインプルークカメラにより前眼部断面画像を得る装置等にも適用可能である。
【００６９】
さらに、シャインプルークカメラを回転させたり、水平または垂直方向に移動させたりすることにより前眼部の３次元形状画像を取得する装置にも適用可能である。この場合、所定の回転角度毎にずれ補正を行うことによって、精度よく前眼部の３次元形状画像を取得することが可能となり、３次元形状画像より取得される測定値の精度が良くなる。この場合、撮像面（スリット断面）に対して垂直な方向の位置ずれが検出され、その検出結果に基づいてずれ補正処理が行われる。
【００７０】
また、上記構成においては、光学的に前眼部断面像を取得したが、これに限定されない。例えば、Ｂスキャン用の超音波プローブを用いて前眼部からの反射情報を検出することによって前眼部断面像を取得する構成であればよい。
【００７１】
なお、本実施形態においては、角膜測定光学系３００を用いて、角膜前面における角膜曲率半径が算出され、ＯＣＴ光学系１００を用いて、角膜後面における角膜曲率半径が算出される構成としたがこれに限定されない。ＯＣＴ光学系１００によって、角膜前後面における角膜曲率半径を算出してもよい。また、角膜前後面の角膜曲率半径を同様の測定値で扱ってもよい。すなわち、角膜測定光学系３００で算出した角膜前面における角膜曲率半径が角膜前後面における角膜曲率半径として用いてもよい。
【００７２】
なお、本実施形態において、角膜測定光学系３００として、角膜トポグラフィーを用いることもできる。この場合、角膜前面の曲率半径を算出する際に、角膜の全体の形状から角膜前面の曲率半径が算出されるため、角膜前面の曲率半径が精度良く算出される。このため、ＩＯＬ度数を算出する際に、ＩＯＬ度数算出の精度向上に繋がる。
【００７３】
なお、ＯＣＴ光学系１００が毛様体を含む前眼部断面像の画像化が可能な前眼部撮像デバイス（例えば、超音波Ｂプローブ、前眼部ＯＣＴ）である場合、制御部８０は、ＯＣＴ光学系１００によって取得された前眼部断面像における毛様体の位置情報に基づいて、補正量Ｘを求めてもよい。例えば、取得された前眼部断層像（前眼部断面像）より、毛様体（毛様体先端部）を検出し、検出した毛様体位置からチン小帯位置を予測する。そして、予測したチン小帯位置からチン小帯と水晶体との接点部分の位置を検出すればよい。
【００７４】
また、ＯＣＴ光学系１００がチン小帯と水晶体の接点部分を含む前眼部断面像の画像化が可能な前眼部撮像デバイス（例えば、超音波Ｂプローブ、前眼部ＯＣＴ）である場合、制御部８０は、ＯＣＴ光学系１００によって取得された前眼部断面像における接点部分を処理して補正量Ｘを求める。例えば、前眼部断層像（断面像）にチン小帯が撮影されていた場合には、取得された前眼部断層像より、チン小帯と水晶体との接点部分の位置を検出すればよい。
【００７５】
なお、本発明においては、本実施形態に記載した装置に限定されない。例えば、上記実施例の機能を行うＩＯＬ度数計算プログラムをネットワークや各種記憶媒体を介して、システムあるいは装置に供給する。そして、システムあるいは装置のコンピュータ（例えば、ＣＰＵ等）がプログラムを読み出し、実行することも可能である。
【符号の説明】
【００７６】
１  眼科撮影装置
３０  前眼部正面撮像光学系
４０  アライメント光学系
５０  投影光学系
７０  モニタ
８０  制御部
８５  メモリ
８４  操作部
１００  ＯＣＴ光学系
３００  角膜測定光学系
【図１】
2020049128.tif 000013
【図２】
2020049128.tif 000014
【図３】
2020049128.tif 000015
【図４】
2020049128.tif 000016
【図５】
2020049128.tif 000017
【図６】
2020049128.tif 000018
【図７】
2020049128.tif 000019
【図８】
2020049128.tif 000020