(57)【特許請求の範囲】
【請求項１】
被測定眼に所定パターンの光束を照射する光源部を含む照射光学系と、
被測定眼からの反射光を受光する受光部を含む受光光学系と、
上記受光部からの測定データを示す受光信号に基づき、測定データとして屈折力又はパワーの形式で被測定眼の光学特性を演算する演算部と、
上記演算部により得られた測定データに応じて、光学系屈折力分布又はパワーマップをグラフィック表示する表示部と
を備え、
上記照射光学系は、被測定眼の眼底に所定パターンとして略点光源を形成し、上記受光光学系は、被測定眼からの反射光を複数の光束に分割するハルトマン板を含み、上記演算部は、上記受光部からのハルトマン像に対応した受光信号により、被測定眼の眼の屈折力分布又はパワーマップを求め、
及び、
上記照明光学系は、被測定眼角膜に対して異なる径のリング状光束を照射し、上記受光光学系は、被測定眼角膜からの反射光を受光し、上記演算部は、上記受光部からのプラチドリング像を示す受光信号により、被測定眼の角膜の屈折力分布又はパワーマップを求め、
さらに、
上記演算部は、眼球内部の光学系屈折力分布を、眼の屈折力分布又はパワーマップと、角膜の屈折力分布又はパワーマップとの差として求め、
上記表示部は、前記眼球内部の光学系屈折力分布又はパワーマップをグラフィック表示する
光学特性測定装置。
【請求項２】
上記表示部は、眼球内部の光学系屈折力分布又はパワーマップと、眼球内部の収差分布又はパワーマップとを、同一画面で表示するようにした請求項１に記載の光学特性測定装置。
【請求項３】
上記表示部は、さらに、眼球内部の光学系屈折力分布又はパワーマップと眼球内部の収差分布又はパワーマップとの表示に関する数値データを表示することを特徴とする請求項２に記載の光学特性測定装置。
【請求項４】
上記表示部は、被測定眼の眼球の屈折力分布又はパワーマップと、被測定眼の角膜の屈折力分布又はパワーマップとを、同一画面でグラフィック表示するようにした請求項１に記載の光学特性測定装置。
【請求項５】
上記表示部は、グラフィック表示された、被測定眼の眼球の屈折力分布又はパワーマップと、被測定眼の角膜の屈折力分布又はパワーマップとを、同一画面で切換えて表示するようにした請求項１に記載の光学特性測定装置。
【請求項６】
上記演算部は、さらに測定データとして収差形式で眼の光学特性を演算し、
上記表示部は、被測定眼の眼球の屈折力分布又はパワー形式と、収差形式とで表示するようにした請求項１乃至５のいずれかに記載の光学特性測定装置。
【請求項７】
上記演算部は、さらに測定データとして収差形式で眼の光学特性を演算し、
上記表示部は、被測定眼の眼球の屈折力分布マップ又はパワーマップ形式と、
収差マップ形式でグラフィック表示するようにした請求項１乃至６のいずれかに記載の光学特性測定装置。
【請求項８】
上記複数の屈折力分布又はパワー形式は、角膜のパワー形式であって、角膜のアキシャルパワーマップ、インスタンテーニアスパワーマップ、レフラクティブパワーマップのいずれかであることを特徴とする請求項７に記載の光学特性測定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、光学特性測定装置に係り、特に、複数の条件下での測定による各種データをグラフィック表示する光学特性測定装置に関する。
【０００２】
【従来の技術】
近年、医学用に用いられる光学機器は、極めて多種多様な広がりを見せている。この光学機器は、特に、眼科では、眼の屈折力、調節等の眼機能、眼球内部の検査を行う光学特性測定装置として普及している。また、これらの各種検査の測定結果は、例えば、検査対象となる患者の被測定眼がどのような測定条件下に置かれていたかが重要となる。例えば、眼の瞳孔は、明るい所では小さく、暗い所では大きくなるため、測定条件として、照度も考慮する必要があり、さらに、被測定眼の測定範囲も重要である。
【０００３】
また、眼に含まれる網膜、角膜、それ以外の部位の形状は、患者によってそれぞれ特有なものである場合が多く、眼科医等が患者の被測定眼に対する診断等を迅速に行うためには、被測定眼の各部位に関する収差等の各種データを、まとめて、又は、所望のデータを選択して表示することが望ましい。これにより、眼科医等は、各種診断（所見）を患者に対してわかり易く説明することもできる。
【０００４】
また、この各種データは、眼の光学特性に相当し、例えば、数値データとイメージデータ（グラフィックデータ）とに区分すると共に、これらの数値データとグラフィックデータとに対して、光学機器の撮像手段（例えば、ＣＣＤ）によって撮影された各種画像を処理又は合成して視覚的に見易い状態にすることが必要である。さらに、各種データの測定結果、測定データ、測定結果に対応する数値データは、複数の条件下で表示が行われる必要がある。
【０００５】
【発明が解決しようとする課題】
しかしながら、従来の光学特性測定装置では、測定データ（測定結果）及び測定結果に対応する画像データや数値データ等、複数の条件下で求められた各種データを、まとめて表示し、又は、選択的にそれぞれ表示し、視覚的に見易くすることは困難である場合が想定される。
【０００６】
本発明は、以上の点に鑑み、複数の条件下で求めた測定データ（測定結果）、測定結果に対応する画像データ及び／又は数値データをまとめて、又は、選択的に表示する光学特性測定装置を提供することを目的とする。
また、本発明は、被測定眼の全体、角膜、眼内等に対して、測定データ（測定結果）、測定結果に対応する画像データ及び／又は数値データをまとめて、又は、選択的にグラフィック表示等を用いて表示する光学特性測定装置を提供することを目的とする。
【０００７】
【課題を解決するための手段】
本発明の解決手段によると、
被測定眼に所定パターンの光束を照射する光源部を含む照射光学系と、
被測定眼からの反射光を受光する受光部を含む受光光学系と、
上記受光部からの測定データを示す受光信号に基づき、測定データとして屈折力又はパワーの形式で被測定眼の光学特性を演算する演算部と、
上記演算部により得られた測定データに応じて、光学系屈折力分布又はパワーマップをグラフィック表示する表示部と
を備え、
上記照射光学系は、被測定眼の眼底に所定パターンとして略点光源を形成し、上記受光光学系は、被測定眼からの反射光を複数の光束に分割するハルトマン板を含み、上記演算部は、上記受光部からのハルトマン像に対応した受光信号により、被測定眼の眼の屈折力分布又はパワーマップを求め、
及び、
上記照明光学系は、被測定眼角膜に対して異なる径のリング状光束を照射し、上記受光光学系は、被測定眼角膜からの反射光を受光し、上記演算部は、上記受光部からのプラチドリング像を示す受光信号により、被測定眼の角膜の屈折力分布又はパワーマップを求め、
さらに、
上記演算部は、眼球内部の光学系屈折力分布を、眼の屈折力分布又はパワーマップと、角膜の屈折力分布又はパワーマップとの差として求め、
上記表示部は、前記眼球内部の光学系屈折力分布又はパワーマップをグラフィック表示する
光学特性測定装置を提供する。
【０００８】
【発明の実施の形態】
以下、図面を用いて本発明の実施の形態を詳細に説明する。
図１は、本発明に関する光学特性測定装置１００の概略光学系を示す図である。
光学特性測定装置１００は、例えば、対象物である被測定眼６０の光学特性を測定する装置であって、第1照明光学系１０と、第1受光光学系２０と、送受光光学系３０と、共通光学系４０と、調整用光学系５０とを備える。なお、被測定眼６０については、図中、網膜６１、角膜６２が示されている。
第1照明光学系１０は、例えば、第１波長の光束を発するための第1光源部１１と、集光レンズ１２とを備え、第1光源部１１からの光束で被測定眼６０の網膜（眼底）６１上の微小な領域を、その照明条件を適宜設定できるように照明するためのものである。なお、ここでは、一例として、第1光源部１１から発せられる照明用の光束の第1波長は、赤外域の波長（例えば、７８０ｎｍ）である。
【０００９】
また、第1光源部１１は、空間コヒーレンスが大きく、時間コヒーレンスが小さいものが望ましい。ここでは、第1光源部１１は、例えば、スーパールミネッセンスダイオード（ＳＬＤ）であって、輝度の高い点光源を得ることができる。なお、第1光源部１１は、ＳＬＤに限られるものではなく、例えば、空間コヒーレンス、時間コヒーレンスが大きいレーザー等であっても、回転拡散板等を挿入し、適度に時間コヒーレンスを下げることで、利用することができる。さらに、空間コヒーレンス、時間コヒーレンスが小さいＬＥＤであっても、光量さえ十分であれば、例えば、光路の光源の位置にピンホール等を挿入することで、利用することができる。
【００１０】
第1受光光学系２０は、例えば、コリメートレンズ２１と、被測定眼６０の網膜６１から反射して戻ってくる光束（第1光束）の一部を、少なくとも１７本のビームに変換する変換部材であるハルトマン板２２と、このハルトマン板２２で変換された複数のビームを受光するための第1受光部２３とを備え、第1光束を第1受光部２３に導くためのものである。また、ここでは、第1受光部２３は、リードアウトノイズの少ないＣＣＤが採用されているが、ＣＣＤとしては、例えば、一般的な低ノイズタイプ、測定用の２０００＊２０００素子の冷却ＣＣＤ等、適宜のタイプのものを適用することができる。
【００１１】
送受光光学系３０は、例えば、後述するアライメント調整を主に行うものであって、第２波長の光束を発するための第２光源部３１と、集光レンズ３２、３４と、ビームスプリッター３３と、第２受光部３５とを備え、被測定眼６０の角膜６２から反射して戻ってくる光束（第２光束）を、第２受光部に導くためのものであって、主にアライメント調整を行う。また、第２光源部３１から発せられる光束の第２波長は、例えば、第１波長（ここでは、７８０ｎｍ）と異なると共に、長い波長を選択できる（例えば、９４０ｎｍ）。
【００１２】
共通光学系４０は、第1照明光学系１０から発せられる光束の光軸上に配され、第1照明光学系１０、第1受光光学系２０及び送受光光学系３０に共通に含まれ得るものであり、例えば、プラチドリング４１と、アフォーカルレンズ４２と、ビームスプリッター４３、４５と、集光レンズ４４とを備える。プラチドリング４１は、後述するアライメント調整が完了した後、複数の同心輪帯からなるパターンの指標を投影する。また、ビームスプリッター４３は、第２光源部３１の波長を被測定眼６０に送光（反射）し、被測定眼６０の角膜６２から反射して戻ってくる第２光束を反射し、一方、第１光源部１１の波長を透過するようなミラー（例えば、ダイクロミックミラー）で形成される。ビームスプリッター４５は、第１光源部１１の波長を被測定眼６０に送光（反射）し、被測定眼６０の網膜６１から反射して戻ってくる第１光束を、透過するようなミラー（例えば、ダイクロミックミラー）で形成される。このビームスプリッター４３、４５によって、第１及び２光束が、互いに他方の光学系に入りノイズとなることがない。
【００１３】
調整用光学系５０は、例えば、後述する作動距離調整を主に行うものであって、第３光源部５１と、第４光源部５５と、集光レンズ５２、５３と、第３受光部５４を備え、主に作動距離調整を行うものである。
ここで、アライメント調整について説明する。アライメント調整は、主に、送受光光学系３０により実施される。
【００１４】
まず、第２光源部３１からの光束は、集光レンズ３２、ビームスプリッター３３、４３、アフォーカルレンズ４２を介して、対象物である被測定眼６０を略平行な光束で照明する。被測定眼６０の角膜６２で反射した反射光束は、あたかも角膜６２の曲率半径の１／２の点から射出したような発散光束として射出される。この発散光束は、アフォーカルレンズ４２、ビームスプリッター４３、３３及び集光レンズ３４を介して、第２受光部３５にスポット像として受光される。
【００１５】
ここで、この第２受光部３５上のスポット像が光軸上から外れている場合、光学特性測定装置１００本体を、上下左右に移動調整し、スポット像を光軸上と一致させる。このように、スポット像が光軸上と一致すると、アライメント調整は完了する。なお、アライメント調整は、被測定眼６０の角膜６２を第３光源部５１により照明し、この照明により得られた被測定眼６０の像が第２受光部３５上に形成されるので、この像を利用して瞳中心が光軸と一致するようにしてもよい。
【００１６】
つぎに、作動距離調整について説明する。作動距離調整は、主に、調整用光学系５０により実施される。
まず、作動距離調整は、例えば、第４光源部５５から射出された光軸付近の平行な光束を、被測定眼６０に向けて照射すると共に、この被測定眼６０から反射された光を、集光レンズ５２、５３を介して第３受光部５４で受光することにより行われる。また、被測定眼６０が適正な作動距離にある場合、第３受光部５４の光軸上に、第４光源部５５からのスポット像が形成される。一方、被測定眼６０が適正な作動距離から前後に外れた場合、第４光源部５５からのスポット像は、第３受光部５４の光軸より上又は下に形成される。なお、第３受光部５４は、第４光源部５５、光軸、第３受光部５４を含む面内での光束位置の変化を検出できればいいので、例えば、この面内に配された１次元ＣＣＤ、ポジションセンシングデバイス（ＰＳＤ）等を適用できる。
【００１７】
つぎに、第1照明光学系１０と第1受光光学系２０との位置関係を概略的に説明する。
第1受光光学系２０には、ビームスプリッター４５が挿入されており、このビームスプリッター４５によって、第1照明光学系１０からの光は、被測定眼６０に送光されると共に、被測定眼６０からの反射光は、透過される。第1受光光学系２０に含まれる第１受光部２３は、変換部材であるハルトマン板２２を通過した光を受光し、受光信号を生成する。
【００１８】
また、第１光源部１１と被測定眼６０の網膜６１とは、共役な関係を形成している。被測定眼６０の網膜６１と第１受光部２３とは、共役である。また、ハルトマン板２２と被測定眼６０の瞳孔とは、共役な関係を形成している。さらに、第1受光光学系２０は、被測定眼６０の前眼部である角膜６２、及び瞳孔と、ハルトマン板２２と略共役な関係を形成している。すなわち、アフォーカルレンズ４２の前側焦点は、被測定眼６０の前眼部である角膜６２及び瞳孔と略一致している。
【００１９】
また、第1照明光学系１０と第1受光光学系２０は、第１光源部１１からの光束が、集光する点で反射されたとして、第１受光部２３での反射光による信号ピークが最大となるように、連動して移動する。具体的には、第1照明光学系１０と第1受光光学系２０は、第１受光部２３での信号ピークが大きくなる方向に移動し、信号ピークが最大となる位置で停止する。これにより、第１光源部１１からの光束は、被測定眼の網膜６１上で集光する。
【００２０】
また、レンズ１２は、光源１１の拡散光を平行光に変換する。絞り１４は、眼の瞳、あるいはハルトマンプレート２１と光学的に共役の位置にある。絞り１４は、径がハルトマンプレート２１の有効範囲より小さく、いわゆるシングルパスの収差計測（受光側だけに目の収差が影響する方法）が成り立つ様になっている。レンズ１３は、上記を満たすために、実光線の眼底共役点を前側焦点位置に、さらに、眼の瞳との共役関係を満たすために、後側焦点位置が絞り１４と一致するように配置されている。
【００２１】
また、光線１５は、光線２４とビームスプリッター４５で共通光路になった後は、近軸的には、光線２４と同じ進み方をする。但し、シングルパス測定のときは、それぞれの光線の径は違い、光線１５のビーム径は、光線２４に比べ、かなり細く設定される。具体的には、光線１５のビーム径は、例えば、眼の瞳位置で１ｍｍ程度、光線２４のビーム径は、７ｍｍ程度になることもある（なお、図中、光線１５のビームスプリッター４５から眼底６１までは省略している）。
【００２２】
つぎに、変換部材であるハルトマン板２２について説明する。
第1受光光学系２０に含まれるハルトマン板２２は、反射光束を複数のビームに変換する波面変換部材である。ここでは、ハルトマン板２２には、光軸と直交する面内に配された複数のマイクロフレネルレンズが適用されている。また、一般に、測定対象部（被測定眼６０）について、被測定眼６０の球面成分、３次の非点収差、Ｚｅｒｎｉｋｅの３次と４次の高次収差までも測定するには、被測定眼６０を介した少なくとも１７本のビームで測定する必要がある。
【００２３】
また、マイクロフレネルレンズは、光学素子であって、例えば、波長ごとの高さピッチの輪帯と、集光点と平行な出射に最適化されたブレーズとを備える。ここでのマイクロフレネルレンズは、例えば、半導体微細加工技術を応用した８レベルの光路長差を施したもので、高い集光率（例えば、９８％）を達成している。
【００２４】
また、被測定眼６０の網膜６１からの反射光は、アフォーカルレンズ４２、コリメートレンズ２１を通過し、ハルトマン板２２を介して、その1次光として第１受光部２３上に集光する。また、ハルトマン板２２は、少なくとも１７の領域に区分された各領域ごとに、収束作用を行うマイクロレンズ部と、透過作用を行う開口部とを備えるようにしてもよい。したがって、ハルトマン板２２は、反射光束を少なくとも１７本以上のビームに変換する波面変換部材を備える。
【００２５】
図２は、本発明に関する光学特性測定装置１００の電気的構成を示す電気系ブロック図である。
光学特性測定装置１００に関する電気駆動系は、例えば、演算部２１０と、制御部２２０と、表示部２３０と、メモリ２４０と、第1駆動部２５０及び第２駆動部２６０とを備える。演算部２１０は、第１受光部２３から得られる受光信号▲４▼、第２受光部３５から得られる受光信号▲７▼、第３受光部５４から得られる受光信号(10)を入力すると共に、全波面収差、角膜波面収差、ゼルニケ係数、収差係数、Ｓｔｒｅｈｌ比、白色光ＭＴＦ、ランドルト環パターン等（詳細は後述する。）を演算し、この演算結果に応じた信号を、電気駆動系の全体の制御を行う制御部２２０と、表示部２３０（各種の表示例については、後述する。）と、メモリ２４０とにそれぞれ出力する。
【００２６】
制御部２２０は、演算部２１０からの制御信号に基づいて、第1光源部１１の点灯、消灯を制御したり、第1駆動部２５０及び第２駆動部２６０を制御するものであり、例えば、演算部２１０での演算結果に応じた信号に基づいて、第１光源部１１に対して信号▲１▼を出力し、プラチドリング４１に対して信号▲５▼を出力し、第２光源部３１に対して信号▲６▼を出力し、第３光源部５１に対して信号▲８▼を出力し、第４光源部５５に対して信号▲９▼を出力し、さらに、第1駆動部２５０及び第２駆動部２６０に対して信号を出力する。
【００２７】
第1駆動部２５０は、例えば、演算部に入力された第1受光部２３からの受光信号▲４▼に基づいて、第1照明光学系１０全体を光軸方向に移動させるものであり、図示しない適宜のレンズ移動手段に対して信号▲２▼を出力すると共に、このレンズ移動手段を駆動する。これにより、第1駆動部２５０は、第1照明光学系１０の移動、調節を行うことができる。
【００２８】
第２駆動部２６０は、例えば、演算部に入力された第1受光部２３からの受光信号▲４▼に基づいて、第1受光光学系２０全体を光軸方向に移動させるものであり、図示しない適宜のレンズ移動手段に対して信号▲３▼を出力すると共に、このレンズ移動手段を駆動する。これにより、第２駆動部２６０は、第1受光光学系２０の移動、調節を行うことができる。
【００２９】
図３は、本発明に関する光学特性測定装置１００の動作を示すフローチャートである。なお、説明の便宜上、ここでのフローチャートの説明は、概略的なものとし、具体的な説明は後述する。
まず、測定者により測定対象物である被測定眼６０の測定が開始され、第1受光部２３又は第２受光部３５からの像（例えば、ハルトマン像、前眼部像）を取得する（Ｓ１０１）。つぎに、眼の収差測定に関しては、ステップＳ１０１でのハルトマン像に関する数値データである、ハルトマン板２２と第1受光部２３との距離及び、座標等に基づいて、ゼルニケ係数を算出し、角膜収差測定に関しては、プラチドリング像のゆがみによる変位量等に基づいてゼルニケ係数を算出する（Ｓ１０３）。このゼルニケ係数に基づいて、被測定眼６０の波面または角膜収差を算出する（Ｓ１０５）。
【００３０】
つぎに、ステップＳ１０１〜１０５から得られた各種データを、表示部２３０にどのように表示するか表示モード選択を行う（Ｓ１０７）。なお、この各種データは、複数の条件下で求めた測定データ（測定結果）、測定結果に対応する画像データ及び／又は数値データとして、まとめて、又は、選択的に表示部２３０にグラフィック表示することができる（詳細は、後述）。また、ステップＳ１０７で選択された表示モードに応じて、各種画像及び／又はデータを表示部２３０に表示する（Ｓ１０９）。
【００３１】
また、ステップＳ１０９での表示例とは別の表示例を、表示部２３０に表示させる場合、ステップＳ１０７で選択した表示モード選択を変更するために、再度、表示モード選択を行う（Ｓ１１１）。一方、ステップＳ１０７で選択した表示モード選択を変更しない場合、測定を終了するか否かを判定する（Ｓ１１３）。ステップＳ１１３で測定を終了しない場合には、再度、像の取得を行う。一方、測定を終了する場合には、測定を終了する。
以下、上述の各ステップの処理について詳細に説明する。
【００３２】
（ステップＳ１０１について）
（表示部２３０に表示される被測定眼６０の像について）
この像は、例えば、ハルトマン板２２上に開けられた開口部を通過したことにより得られる、いわゆるハルトマン像であると共に、撮像手段であるＣＣＤを適用した第1受光部２３によって撮像された生画像である。また、例えば、患者の被測定眼６０に対して第１受光光学系２０の調整が不足気味であれば、ハルトマン板２２で変換された複数のビームは、概ね外側に広がる傾向にある（後述する表示例に含まれるハルトマン像上に示した矢印に対応）。
【００３３】
（ステップＳ１０３、１０５について）
（眼の収差マップ、角膜収差マップ、俯瞰図と、数値データである眼の収差表示部、角膜収差表示部とを、表示部２３０に表示する際、必要とされるゼルニケ係数について）
ここで、一般に知られているゼルニケの多項式Ｚｉｊ（対象物で変化しない固定値であって、ここでは、既知である）を用いてゼルニケ係数Cｉｊを算出する方法について説明する。
【００３４】
まず、ゼルニケ係数Cｉｊは、被測定眼の光学特性を把握するための重要なパラメータである。このパラメータは、例えば、ハルトマン板２２を介して第1受光部２３で得られた光束の傾き角に基づいて、被測定眼６０の光学特性を求める演算部２１０の動作原理となる。
光学特性測定装置１００は、被測定眼６０の波面収差（Ｗ（Ｘ，Ｙ）とする）を測定する。このため、変換部材であるハルトマン板２２の縦横の座標を（Ｘ，Ｙ）とし、さらに、第１受光部２３の縦横の座標を（ｘ，ｙ）とすると、一般に、以下の数式３で表される波面収差Ｗ（Ｘ，Ｙ）は、同じく、以下の数式１及び数式２の関係が成り立つ。すなわち、
【００３５】
【数１】
0004652558.tif 000002
【００３６】
【数２】
0004652558.tif 000003
【００３７】
ここで、ｆ：ハルトマン板２２と第1受光部２３との距離である。
【００３８】
【数３】
0004652558.tif 000004
【００３９】
また、数式３の両辺を、ハルトマン板２２上の座標（Ｘ，Ｙ）で偏微分することにより、ゼルニケ係数Cｉｊを得る。
ここで、ゼルニケの多項式Ｚｉｊは、以下の数式４及び数式５で表され、具体的には、図１９、２０に示される。
【数４】
0004652558.tif 000005
【数５】
0004652558.tif 000006
また、ゼルニケ係数Cｉｊは、以下の数式６の自乗誤差を最小にすることにより、具体的な値を得ることができる。
【数６】
0004652558.tif 000007
【００４０】
（ステップＳ１０７、１０９について）
ここで、複数の条件下で求めた測定データ（測定結果）、測定結果に対応する画像データ及び／又は数値データを、まとめて、又は、選択的に、表示部２３０上に表示する場合の表示例について説明する。
まず、ステップＳ１０７による表示モード選択では、例えば、測定者が後述する第１〜１３の表示例のうち所望の表示例を選択することができる。なお、測定者は、第１〜１３の表示例に含まれる各要素を選択的に表示することもできる。
ここで、説明の便宜上、第１〜１３の表示例に含まれる各要素についてそれぞれ説明する。
【００４１】
（表示部２３０に表示される眼の収差表示部に含まれる収差係数について）
被測定眼の収差成分を、例えば、収差係数（ここでは、単位：μｍ）として数値化する場合、ゼルニケ係数（ｎ＝１、２、３、４、５、６）であるとき、収差係数のＲＭＳ値（平均２乗誤差）をそれぞれＳ１、Ｓ２、Ｓ３、Ｓ４、Ｓ５、Ｓ６とすると、
【００４２】
【数７】
0004652558.tif 000008
【００４３】
となる。ここで、Ｓ１は、チルトを意味し、Ｓ３、Ｓ５、Ｓ（３＋５）は、コマ様収差となり、同じく、Ｓ２、Ｓ４、Ｓ６、Ｓ（４＋６）は、球面様収差となる。なお、通常、Ｓ２を除きＳ３以上の収差を高次収差と呼ぶ。
つぎに、この高次収差の分類について説明する。
高次収差は、例えば、ゼルニケ係数Cｉｊ中のｓｉｎ（ｎθ）、ｃｏｓ（ｎθ）のｎの値により分類される（ゼルニケ係数Cｉｊは、上述の数式３及び数式４に示すように、正弦及び余弦関数を含んでいる）。具体的には、ｎ＝０（fold）であれば、球面収差となり、ｎ＝１（fold）であれば、コマ収差となり、同じく、ｎ＝２（fold）であれば、非点収差となり、さらに、ｎ＝３（fold）であれば、矢状収差とそれぞれ分類される。
また、対称収差（コマ様収差）は、４次収差＋６次収差を意味し、同じく、非対称収差（球面様収差）は、３次収差＋５次収差を意味する。なお、高次収差の次数は、必要に応じて決められるものであって、６次までに限らず、例えば、８次では３次＋５次＋７次がコマ様収差、４次＋６次＋８次が球面様収差をそれぞれ意味する。
つぎに、ゼルニケ係数Cｉｊを用いて、上述の各収差を表現する場合について説明する。
図１８は、ゼルニケ係数Cｉｊによる各収差の表現形式６００を示す図である。
各収差の表現形式６００としては、例えば、各収差の名称、ゼルニケ収差係数、方向Ａｘ（例えば、Ｘ軸から反時計回り）、及び上述のＲＭＳの平方にそれぞれ項目分けされる。また、各収差の表現形式６００は、図示のように、３次球面収差を示す項目６０１、５次球面収差を示す項目６０２、７次球面収差を示す項目６０３、３次コマ収差を示す項目６０４、矢状収差を示す項目６０５、５次コマ収差を示す項目６０６、３次非点収差を示す項目６０７、及び５次非点収差を示す項目６０８を含む。なお、これらの各収差は、主に、後述する第５〜８表示例において、眼及び角膜の収差表示として、一括、又は、選択的に表示部２３０上に表示される。
つぎに、Power Map計算方法について説明する。
まず、求められた波面収差Ｗ（Ｘ，Ｙ）をパワー表示に変換する場合について説明する。なお、このパワー表示は、主に後述する第９〜１３表示例において示される。
ここでは、対象とする波面の位置をＸ，Ｙで表示するものとし、光軸を含む断面において、その位置での波面の法線が光軸と交わる点と、その時に光軸上での波面の位置までの距離をLPとする。
このときのパワーＰを1/LP+1/Leyeとする。ここで、Leyeは、眼の平均眼軸長（例えば０．０１７（ｍ））である。
また、波面の各位置(r, t)でのパワー分布をP(r, t)で表す。このパワー分布P(r, t)は、後述する眼の屈折力分布(Ocular Refractive Power)に相当する。つぎに、上述の数式３での波面収差Ｗ（Ｘ，Ｙ）をパワー分布P(r, t)に置き換えて、上述と同様の処理を行い、パワー分布に相当するゼルニケ係数Cｉｊを求める。なお、ここでは、ゼルニケ多項式Ｚｉｊは、図１９で示した極座標系のものを用いる。ここで、求めたゼルニケ係数を利用して、パワー分布に相当する球面収差、コマ収差、非点収差などの各収差を求める。
また、パワー分布をゼルニケ多項式Ｚｉｊ近似した場合、収差の種類は、下記のように分類して表示される。具体的には、０次では、平均パワー、sin又はcosが無い項では、球面収差、sin(ｔ)又はcos(ｔ)では、コマ収差、sin(2ｔ)又はcos(2ｔ)では、非点収差、sin(3ｔ)又はcos(3ｔ)では、矢状収差とそれぞれ分類される。
つぎに、プラチドリングの角膜前面からの反射像から角膜前面のPower Mapを算出する。
ここで、各リングの中心位置を検出する場合は、基準となる中心位置を決定する。例えば、最も内側のプラチドリング像の中心位置を求め、各リングに対応させる。最も内側のプラチドリング像の中心位置は、例えば、リングを楕円であると仮定し、リング上の各点を楕円等の式に最小２乗法により関数近似し、その時の中心を算出する。また、各リング上の各点の位置を検出する場合は、求められた各リングの中心位置から各リング上の各点の位置までの距離rを算出する。
つぎに、形状、Power算出について説明する。
まず、面形状を距離rにおいてあるモデルカーブf(r)と仮定し、断面形状を算出する。なお、モデルカーブは、例えば、非球面などが用いられる。このとき求められたリングの中心位置を原点とすると、
【数８】
0004652558.tif 000009
となる。ここで、ｐ：コーニック係数、c：中心位置における曲率である。
また、この式の1次微分は、
【数９】
0004652558.tif 000010
となる。
図２１（a）は、形状、power算出についての説明図（１）である。
ここで、反射点における傾きは、例えば、df(r)/dr=tan(I/2)=tanθ となり、各点での曲率半径に依存しない値となる。また、1次微分の式にr、df(r)/dr を最小２乗法により関数近似し、各係数A2、A4、…を求める。こうしてf(x)を決めることにより、断面上の各点における曲率半径Ri(r)を求めることができる。すなわち、
【数１０】
0004652558.tif 000011
となる。ここで、f’(r) =d(fr)/dr、 f”(r) = d2f(r)/dr2 である。
また、数式１０より、Instantaneous Power Piは、
Pi(r)=n-1/Ri(r)
となる。ここで、n：角膜屈折力である。
また、このとき断面上の各点と、各点の法線と角膜光軸との交点までの距離Ra(r)は、
【数１１】
0004652558.tif 000012
となる。また、数式１１より、Axial Power Pa(r)は、
Pa(r) = n-1/Ra(r)
により求められる。
図２１（ｂ）は、形状、power算出についての説明図（２）である。
ここで、断面上の各点に角膜光軸に平行な光線が入射してきたときスネルの法則によって、
sinθ = n sinα → α = sin-1(1/n sinθ) (θ=I/2)
となる。
また、断面上の各点と、光線と角膜光軸との交点の距離のx成分(光軸方向)Qは、
Q = L/tanU = tan(θ-α)
となる。これにより、中心位置から光線と角膜光軸との交点間の距離Rr(r)は、
Rr(r) = f(r) + Q
により求められる。このRr(r)から、Refractive Power Pr(r)は、
Pr(r)=n/Rr(r)
となる。
このように、これら３種類のパワーでは、１断面におけるパワーが求められている。動径ｒの方位ｔ、全ての方位ｔでのそれぞれのパワーP(r,t)を求める。なお、ｔは、ＸＹ平面上での動径ｒ方位を意味する。また、上述の数式３における波面収差Ｗ（Ｘ，Ｙ）を、パワー分布P(r,t)に置き換えて、上述と同様の処理を行い、パワー分布に相当するゼルニケ係数を求める。
【００４４】
（表示部２３０に表示される白色光ＭＴＦについて）
つぎに、白色光ＭＴＦ（Modulation transfer function）の算出について説明する。
まず、ＭＴＦは、空間周波数の伝達特性を示す指標であって、光学系の性能を表現するために広く使われている。このＭＴＦは、例えば、１度当たり、０〜１００本の正弦波状の濃淡格子に対しての伝達特性を求めることで見え方を予測することが可能である。
【００４５】
まず、単色ＭＴＦを波面収差Ｗ（ｘ，ｙ）から算出する。なお、Ｗ（ｘ，ｙ）は、入力値（測定値）であって、角膜収差に関しては、角膜形状から求めた角膜波面収差を用いることもできる。
瞳関数ｆ（ｘ，ｙ）は、波面収差から以下のように求まる。
ｆ（ｘ，ｙ）＝ｅｉｋＷ（ｘ，ｙ）
ここで、ｉ：鏡像、ｋ：波数ベクトル（２π／λ）である。
また、この瞳関数ｆ（ｘ，ｙ）をフーリエ変換することにより、点像の振幅分布Ｕ（ｕ，ｖ）が求まる。
【００４６】
【数１２】
0004652558.tif 000013
【００４７】
ここで、λ：波長、Ｒ：瞳から像点（網膜）までの距離、（ｕ，ｖ）：網膜上の像点Ｏを原点とし、光軸に直行する面内での網膜の座標値、（ｘ，ｙ）：光学系の瞳面内の座標値である。
また、点像の振幅分布Ｕ（ｕ，ｖ）とその複素共役を掛けることにより、点像の強度分布（ＰＳＦ）Ｉ（ｕ，ｖ）が求まる。
Ｉ（ｕ，ｖ）＝Ｕ（ｕ，ｖ）Ｕ＊（ｕ，ｖ）
さらに、点像の強度分布Ｉ（ｕ，ｖ）をフーリエ変換すると共に、いわゆる空間周波数変換である（Ｒ（ｒ，ｓ））規格化を行うことにより、ＯＴＦ（Optical Transfer Function）が求まる。
【００４８】
【数１３】
0004652558.tif 000014
【００４９】
ここで、ｒ，ｓ：空間周波数領域の変数である。
ＯＴＦ（ｕ，ｖ）＝Ｒ（ｒ，ｓ）／｜Ｒ（０，０）｜
また、ＯＴＦの大きさがＭＴＦであるため、
ＭＴＦ（ｒ，ｓ）＝｜ＯＴＦ（ｕ，ｖ）｜
が成り立つ。
つぎに、上述のように求められた単色ＭＴＦに基づいて、白色光ＭＴＦを算出する。
白色光ＭＴＦを求めるには、まず、各波長でのＭＴＦに重み付けをし、足し合わせる。ここで、上述のＭＴＦは、波長ごとに値が異なるため、波長λでのＭＴＦをＭＴＦλと表すと、
【００５０】
【数１４】
0004652558.tif 000015
【００５１】
ここでは、可視光に多く重み付けをし、計算を行う。
具体的には、色の３原色（ＲＧＢ）である赤、緑、青が、例えば、６５６．２７ｎｍ：１、５８７．５６ｎｍ：２、及び４８６．１３ｎｍ：１であるとすると、
ＭＴＦ（ｒ，ｓ）＝（１×ＭＴＦ６５６．２７＋２×ＭＴＦ５８７．５６＋１×ＭＴＦ４８６．１３）／（１＋２＋１）
となる。
【００５２】
また、白色光ＭＴＦは、一波長（８４０ｎｍ）のみで測定されるので、この測定結果に基づいて他の波長について校正を行い、白色に補正することにより求めてもよい。具体的には、各波長でのＭＴＦは、眼の収差の場合、眼光学特性測定装置での測定波長が、例えば、８４０ｎｍであるとき、模型眼により各波長８４０ｎｍでの波面収差Ｗ８４０（ｘ，ｙ）からのずれ量に相当する色収差ＷΔ（ｘ，ｙ）を測定し、この色収差ＷΔ（ｘ，ｙ）にＷ８４０（ｘ，ｙ）を足し合わせ、この波面収差によりＭＴＦを算出することにより求められる。すなわち、
【００５３】
Ｗλ（ｘ，ｙ）＝Ｗ８４０（ｘ，ｙ）＋ＷΔ（ｘ，ｙ）
となる。
さらに、角膜波面収差の場合、測定された角膜形状は波長に依存しないが、形状を波面収差に変換するときに使用する角膜の屈折率が波長に依存すること、瞳関数の式に波長がパラメータのひとつであることから、ＭＴＦがこの場合も波長に依存することになる。
【００５４】
（表示部２３０に表示される眼の収差表示部及び角膜収差表示部に含まれるＳｔｒｅｈｌ比について）
Ｓｔｒｅｈｌ比（ここでは、Ｓ）は、上述のように求められた点像の強度分布であるＰＳＦの中心強度Ｉ（０，０）を、無収差光学系の場合に得られるＰＳＦの中心強度Ｉ０（０，０）で割ることにより、求められる。すなわち、
Ｓ＝Ｉ（０，０）／Ｉ０（０，０）
となる。
【００５５】
（表示部２３０に表示されるランドルト環、例えば、視力検査用のマークの表示、見え方について）
図４は、ランドルト環の表示に関するフローチャートである。
まず、測定者（例えば、眼科医等）は、患者の視力を検査するために適宜の視力に応じた大きさを有するランドルト環を、上述のハルトマン板２２に取り付ける。第1受光部２３は、この取り付けられたランドルト環の大きさに基づいて、視力いくつのランドル環であるかを判定する（Ｓ２０１）。このランドル環のパターンを示す輝度分布関数Ｌａｎｄ（ｘ，ｙ）を求める（Ｓ２０３）。ここで、輝度分布関数Ｌａｎｄ（ｘ，ｙ）は、ハルトマン板２２と第1受光部２３との間でのＸＹ座標のずれ（上述の撮影した生画像上に付与した矢印の長さをＸ方向、Ｙ方向にそれぞれベクトル分解した値）に基づいて求められる。
【００５６】
つぎに、この輝度分布関数Ｌａｎｄ（ｘ，ｙ）に対して２次元フーリエ変換を行うことにより、ＦＲ（ｕ，ｖ）を求める（Ｓ２０５）。このＦＲ（ｕ，ｖ）と、既に求めた上述のＯＴＦ（ｕ，ｖ）とをコンボリューションすることにより、眼の光学系通過後の周波数分布であるＯＲ（ｕ，ｖ）を求める（Ｓ２０７）。このＯＲ（ｕ，ｖ）に対して２次元逆フーリエ変換を行うことにより、Ｌａｎｄ  Ｉｍａｇｅ（Ｘ，Ｙ）を求める（Ｓ２０９）。
【００５７】
このステップＳ２０９で求められたＬａｎｄ  Ｉｍａｇｅ（Ｘ，Ｙ）は、表示部２３０の表示画面上にグラフィック表示される（Ｓ２１１）。なお、ここでは、瞳孔の大きさがφ３（明視野）、及びφ７（暗視野）両方での「視力０．７ランドルト環の見え方」として表示されている。
以下、上述の各要素を、まとめて、又は選択的に表示した第１〜１３表示例について説明する。但し、表示部２３０に表示される同一要素には同一符号を付し、重複する説明を省略する。
【００５８】
（第１の表示例）
図５は、表示部２３０にグラフィック表示される第１の表示例を示す説明図である。
表示部２３０には、例えば、撮影した生画像であるハルトマン像３００と、眼の収差マップ３１０と、眼の収差表示部３２０と、白色光ＭＴＦ表示３３０と、ランドル環の見え方表示３４０とが表示される。なお、この表示部２３０には、例えば、患者の名前（ここでは、東京  光子）、測定時刻（ここでは、２０００年３月２日  午前８時４分）、複数の測定条件（ここでは、測定の種類である波面測定、被測定眼６０の種類である無散瞳、測定波長である８４０ｎｍ、被測定眼６０の測定範囲であるφ９ｍｍ）等が表示される。
【００５９】
ハルトマン像３００は、例えば、患者の被測定眼６０の網膜６１からの反射光に基づいた画像であって、ここでは、患者のまぶた３０１も表示されている。また、ハルトマン像３００は、ハルトマン板２２を介して複数の略平行な光束が第1受光部２３上に受光されたと仮定した場合の像点３０２と、網膜６１からの反射光が、ハルトマン板２２に開けられた透過部又は開口部を介して、概ね外側に広がった光束として第1受光部２３上に受光された場合での領域点３０３とを含む。また、このハルトマン像３００に表示されている実線３０４（始点を像点３０２、終点を領域点３０３の中心とする）を第1受光部２３の縦横の座標軸に沿って、ベクトル分解した場合に得られる値は、光束の変位を数値化したものであって、上述のゼルニケ係数を求める際に必要とされるΔｘ、Δｙに相当する。
【００６０】
眼の収差マップ３１０は、例えば、明視野での瞳孔（ここでは、φ３ｍｍ）を示す同心円３１１と、暗視野での瞳孔（ここでは、φ７ｍｍ）を示す同心円３１２と、光束のずれを考慮したゼルニケ係数から算出された波面収差を等高線で表示した複数の略楕円状の環３１３とを含む。なお、この眼の収差マップ３１０の外周縁には、例えば、１０°毎の目盛り（０〜３６０）が表示され、視覚的にも見易いようになっている。また、眼の収差マップ３１０の近傍には、被測定眼６０の測定範囲に応じたスケール（例えば、１．０μｍ毎、−５．０〜５．０の目盛り）３１４が表示される。
【００６１】
ここで、波面収差として、被測定眼の低次収差を含めた全ての収差、又は、高次（３次以上）の収差が必要に応じて表示される。この波面収差の表示は、必要に応じて、色彩や、濃度などを変えてグラフィック表示が行われる。また、全ての収差は、２次以上の収差を初期値として表示するが、１次以上の収差を全収差として表示するよう選択することもできる。
【００６２】
眼の収差表示部３２０には、明視野（ここでは、φ３ｍｍ）及び暗視野（ここでは、φ７ｍｍ）での各種の数値データが表示される。具体的には、上述したゼルニケ係数を用いて算出された被測定眼６０の収差成分（例えば、コマ様収差、球面様収差）であって、数値データである収差係数Ｓ３、Ｓ４、Ｓ５、Ｓ６、Ｓ３＋Ｓ５、Ｓ４＋Ｓ６、Ｓ３＋Ｓ４＋Ｓ５＋Ｓ６と、上述の見易さの指標として使用される、いわゆるＳｔｒｅｈｌ比の数値データと、１／眼の焦点距離である球面度数（ここでは、Ｓで表示され、単位Ｄ：いわゆるディオプター値）と、乱視度数（ここでは、Ｃで表示され、単位Ｄ）と、乱視軸（ここでは、Ａで表示され、単位°）と、視力（ここでは、ＰＶＡ）とがそれぞれ表示される。ここで、球面度数と乱視度数はゼルニケ係数の２次項（Ｓ２）から求められる。なお、上述した眼の収差マップ３１０上に表示される複数の略楕円状の環３１３は、これらの収差係数を考慮している。
【００６３】
白色光ＭＴＦ表示３３０は、上述したように、正弦波格子などの見え方を示す指標であって、ここでは、横軸を空間周波数（cycle/deg）とし（例えば、１°当たり０〜１００本の黒線を施した場合に相当する値）、縦軸を白黒パターンの見え方の程度（０〜１）とするグラフとして表示される。具体的には、ここでは、明視野（ここでは、φ３ｍｍ）で視力０．９である場合でのグラフと、暗視野（ここでは、φ７ｍｍ）で視力０．７である場合でのグラフとが比較できるように、まとめて表示されている。なお、グラフ上での斜線は、正常範囲を示している。
【００６４】
ランドル環の見え方表示３４０は、上述したように、視力検査用のマークが患者の網膜６１上でどのように見えているのかを示している。ここでのランドル環の見え方表示３４０は、例えば、明視野（ここでは、φ３ｍｍ）で、この測定結果である場合、視力０．７に対応するランドル環は、ぼやけることなく細く良く見えており、暗視野（ここでは、φ７ｍｍ）で、この測定結果である場合、多少ぼやけるために太く見えていることを示している。
【００６５】
図６（ａ）は、第1の表示例に対する変形例を示す説明図（１）である。
ここでは、第1の表示例での表示部２３０に対して、撮影した生画像であるハルトマン像３００に、波面収差を測定するのに使用した像点の重心位置を、黒塗り三角形３５１として示したハルトマン像３５０を表示した例を示している。このハルトマン像３５０における図中、黒塗り三角形３５１は、ハルトマン板に設けられた格子点（小さな穴、開口又は透過部分）と測定された照射領域点の対応がとれているものであって、照射領域での重心位置に配置される。その重心位置は、照射強度と照射位置に基づき、いわゆるモーメント法などの手法によって求めることができることを示す。
【００６６】
また、ここでは、第1の表示例での表示部２３０に対して、眼の収差マップ３１０に加えて、対応がとれ、測定結果の得られたハルトマン板に設けられた格子点をオーバーレイ表示した眼の収差マップ３６０を表示した例を示している。
【００６７】
図６（ｂ）は、第１の表示例に対する変形例を示す説明図（２）である。
ここでは、第１の表示例での表示部２３０に対して、撮影した生画像であるハルトマン像３００それ自体３５５を表示した例を示している。
【００６８】
図７は、第1の表示例に対する変形例を示す説明図（３）である。
ここでは、第1の表示例での表示部２３０に対して、眼の収差マップ３１０の代わりに、被測定眼６０の測定範囲での３次元形状を示す俯瞰マップ３７０を表示した例を示している。俯瞰マップ３７０は、例えば、測定範囲φ９ｍｍに応じて、瞳の中心点を原点とした平面用スケール（ここでは、−４．５〜４．５）と、立体用スケール（ここでは、−５．０〜＋５．０μｍ）とを用いて３次元形状を表示している。なお、表示部２３０９では、上述のハルトマン像３００、３５０と、眼の収差マップ３１０、３６０と、俯瞰マップ３７０とを用いて、適宜の表示パターンを選択して作成することができる。
【００６９】
また、上述の表示例によれば、明所時（瞳径小、例えば、φ３ｍｍ）、及び暗所時（瞳径大、例えばφ７ｍｍ）での収差係数の表示を行うことで、両者を比較することができる。また、例えば、眼科医、検査技師等が、表示された収差係数の値が大きく、何らかの対策（例えば、通院、入院等による治療）が必要であると判断する貴重なデータとなり得る数値データについては、注意度合いに応じて、例えば、正常値から離れていくにつれて、青→水色→黄緑→黄色→赤等のように色が変化するように表示することもできる。
【００７０】
また、白色光ＭＴＦ表示３３０についても、明視野及び暗視野での両方の表示を行うことができる。また、測定によって得られた係数から矯正視力の表示を行ってもよい。さらに、この矯正視力に基づいてランドルト環の見え具合を推測したものを表示してもよい。また、収差の単位は、μｍ、ｎｍ、λ等のいずれかを選択できるようにしてもよい。また、両眼表示の選択を行えるようにしても良い。
【００７１】
（第２の表示例）
図８は、表示部２３０にグラフィック表示される第２の表示例を示す説明図である。
表示部２３０には、例えば、プラチド測定の表示３８０と、角膜収差マップ３９０と、角膜収差表示部４００と、白色光ＭＴＦ（φ３、φ７ｍｍ）の表示３３０と、ランドル環の見え方の表示３４０とが表示されている。なお、ここでの被測定眼６０の測定範囲は、例えば、φ１２ｍｍである。
【００７２】
プラチド測定の表示３８０は、上述したように、アライメント調整が完了した後、複数の同心輪帯３８１からなるパターンの指標を投影するプラチドリング４１を用いて、第３光源５１からの光束を、被測定眼６０の角膜６２上に照射すると共に、この複数の同心輪帯３８１を伴った角膜６２からの反射光を、第２受光部３５上で受光した受光信号に基づく表示である。また、ハルトマン板２２を取り除くか調整することにより、受光部２３で測定することも可能である。測定されたプラチドリング４１の各座標における同心円からのずれをΔｘ，Δｙとして、前述と同様にゼルニケ係数を算出し、さらに波面を算出することにより角膜収差を求めることができる。
【００７３】
角膜収差マップ３９０は、明視野での瞳孔（ここでは、φ３ｍｍ）を示す同心円３１１と、暗視野での瞳孔（ここでは、φ７ｍｍ）を示す同心円３１２と、光束のずれを考慮したゼルニケ係数から算出された角膜収差を表示した等高線３９１とを含む。なお、この角膜収差マップ３９０の外周縁には、例えば、１０°毎の目盛り（０〜３６０）が表示され、視覚的にも見易いようになっている。また、角膜収差マップ３９０の近傍には、被測定眼６０の測定範囲（ここでは、φ１２ｍｍ）に応じたスケール（例えば、１．０μｍ毎、−５．０〜５．０の目盛り）３１４が表示される。
【００７４】
角膜収差表示部４００は、上述の眼の収差マップ３１０に対して、１／（眼の焦点距離（単位メートル））である球面度数Ｓ、乱視度数Ｃ、及び乱視軸Ａの代わりに、角膜のみを光学系としたときの度数である角膜収差（ここでは、Ｄで表示され、単位Ｄ）と、曲率半径ｍｍと、軸方向（ここでは、Ａで表示され、単位°）と、最大曲率半径となるところを示すＨと、最小曲率半径となるところを示すＶとが表示される。なお、眼の収差表示部２３０が眼全体（すなわち、網膜６１）に基づいて算出されたものであるのに対し、角膜収差表示部４００は角膜６２だけに基づいて算出されたものであるので、眼の収差表示部３２０、及び角膜収差表示部４００に表示される具体的な数値データは、図示のように異なっている。
【００７５】
（第３の表示例）
図９は、表示部２３０にグラフィック表示される第３の表示例を示す説明図である。
ここでは、第２の表示例での表示部２３０に対して、プラチド測定の表示３８０の代わりに、角膜形状によるAxial Powerマップ４１０、Refractive Powerマップ４２０、Instantaneous Powerマップ４３０のうちいずれかを、まとめて、又は、選択的に表示させる。但し、ここでの測定範囲は、例えば、φ９ｍｍであって、上述の角膜収差マップ３９０の測定範囲もこれと同等にして表示する。
また、角膜形状によるPowerマップ４１０、４２０、４３０の近傍には、被測定眼６０の測定範囲（ここでは、φ９ｍｍ）に応じたスケール（例えば、０．５Ｄ毎、３５．５〜５２．０（Ｄ）の目盛り）を表示してもよい（後述）。
【００７６】
（第４の表示例）
図１０は、表示部２３０にグラフィック表示される第４の表示例を示す説明図である。
この表示部２３０には、例えば、眼の収差マップ３１０と、角膜収差マップ３９０と、この眼及び角膜の収差マップ３１０、３９０に基づいて算出された眼内（内部）収差マップ４４０と、眼内（内部）収差表示部４５０とが表示されている。また、眼内（内部）収差マップ４４０は、例えば、眼内（内部）収差マップ４４０に表示される眼内収差係数に基づいて表示される。この眼内収差係数は、例えば、眼の収差マップ３１０を表示させるための眼全体の収差係数から、角膜収差マップ３９０を表示させるための角膜収差係数を引くことにより算出される数値データである。
この表示例によれば、眼内（内部）収差を表示することにより、眼の収差に与える角膜収差以外の影響を考慮することができる。
【００７７】
図１１は、第４の表示例に対する変形例を示す説明図である。
ここでは、第４の表示例での表示部２３０に対して、眼の収差マップ３１０、角膜収差マップ３９０、眼内（内部）収差マップ４４０に加えて、対応がとれ測定結果の得られたハルトマン板に設けられた格子点を各収差のカラーマップにオーバーレイ表示した、眼の収差マップ４６０、角膜収差マップ４７０、眼内（内部）収差マップ４８０をそれぞれ表示した例を示している。
【００７８】
図１２は、角膜収差測定を示すフローチャートである。
ここでは、特に、角膜収差測定の場合におけるゼルニケ係数の算出（Ｓ１０３）及びゼルニケ係数から波面算出（Ｓ１０５）について詳細に説明する。
まず、第２受光部３５からの信号に基づき、角膜頂点を基準とし、プラチドリングの受光位置に応じて角膜形状の高さを示す角膜形状のマップ(Hight Map)のデータを算出する（Ｓ３０１）。ステップＳ３０１で求めた角膜形状になるべくフィットする参照球面の形状を算出する（Ｓ３０２）。これにより、ゼルニケ係数の算出精度の向上が図れる。測定範囲（例えば、φ３、φ７）に応じて必要な個所を求めれば足りる。
【００７９】
つぎに、角膜形状の成分から参照球面の成分を減算する（Ｓ３０３）。これにより、参照球面との相違のみの残差成分が求められる。ここで、参照球面の球面収差を算出する（Ｓ３０４）。ステップＳ３０３で求めた残差成分の波面収差を算出する（Ｓ３０５）。また、測定波面と参照球面の波面収差を合成したのちにゼルニケ係数を算出する第１測定モードと、ゼルニケ係数を測定波面と参照球面の波面収差のそれぞれの収差に対して求めて、ゼルニケ係数を合成する第２測定モードとの間で選択が行われる（Ｓ３０６）。ここで、第１測定モードの選択がされるとステップＳ３０７に進み、第２測定モードが選択されるとステップＳ３０９に進む。
【００８０】
第１測定モードでは、ステップＳ３０４で求められた参照球面の波面収差と、ステップＳ３０５で求められた残差成分の波面収差を加えた後、それらの波面収差を角膜波面収差として求める（Ｓ３０７）。さらに、ステップＳ３０７で求められた角膜波面収差のゼルニケ係数を算出する（Ｓ３０８）。なお、このゼルニケ係数は、角膜収差を示している。これが終了すると、図３のステップＳ１０７に進み、表示モードが選択されその後の処理に進む。
【００８１】
一方、ステップＳ３０６において、第２測定モードが選択されると、参照球面の波面収差からゼルニケ係数を算出する（Ｓ３０９）。つぎに、ステップＳ３０５で求められた残差成分の波面収差からゼルニケ係数を算出する（Ｓ３１０）。ステップＳ３０９、Ｓ３１０で求められたゼルニケ係数を合成して、角膜収差を求める（Ｓ３１１）。これが終了すると、図３のステップＳ１０７に進み、表示モードが選択されその後の処理に進む。
【００８２】
（第５の表示例）
図１３は、表示部２３０にグラフィック表示される第５の表示例を示す説明図である。
表示部２３０は、例えば、眼の収差と角膜の収差とを一括表示したものであって、眼の収差表示としては、上述の第1の表示例で示した撮影した生画像であるハルトマン像３００（ここでは、Ｉｍａｇｅと表示し、詳細な画像は、省略する）と、少なくとも２次以上の収差を含む被検眼の全収差３１５と、３次以上の収差を含む高次収差３１６とを含み、さらに、角膜の収差表示としては、上述の第２の表示例で示したプラチド測定の表示であるプラチド像３８０（ここでは、Ｉｍａｇｅと表示し、詳細な画像は、省略する）と、上述の第３の表示例で示したAxial Power マップ４１０と、角膜高次収差マップ３９５とを含む。さらに、表示部２３０には、眼の収差と角膜の収差に対応した数値データを含む眼及び角膜の収差表示部４９０が表示される。
【００８３】
また、全収差３１５、及び高次収差３１６の表示形態は、第１の表示例で詳述した眼の収差マップ３１０と略同様であり、ここでは、その詳細な表示を省略する。また、角膜高次収差マップ３９５の表示形態は、第２の表示例で詳述した角膜収差マップ３９０と略同様であり、ここでは、その詳細な表示を省略する。
また、ここでの各マップの近傍に表示されているスケール３１４は、第１の表示例で詳述したスケールと同一である。なお、スケール３１４スケールの単位は、適宜変更可能であって、例えば、全収差３１５のスケール単位としては、−５．０〜５．０μｍだけでなく、−１０．０〜１０．０μｍ、又は、−１５．０〜１５．０μｍに設定してもよい。このように、スケール単位を変更することで、球面度数が大きいときの変化をわかりやすく表示することもできる。
【００８４】
また、Axial Power マップ４１０の近傍には、ディオプター値（Ｄ）を示すスケール３１９が表示される。各マップの近傍に表示されているスケール３１９は、ここでは、第２の表示例で詳述したスケールと同一である。なお、スケール３１９のスケールの単位は、適宜変更可能であって、３５．５〜５２．０（Ｄ）に限らない。
また、ここでは、角膜のPower mapとして、Axial Power マップ４１０を表示しているが、このAxial Power マップ４１０の代わりに、例えば、第３の表示例で示したRefractive Power４２０、Instantaneous Power４３０、及びゼルニケ係数から求めた後述する角膜の全収差の球面成分にあたる項の係数Ｃ２１を眼の全収差を表すゼルニケ係数のうち球面成分にあたる項の係数Ｃ２１に置き換えた２次以上もしくは１次以上（選択可）の全収差を、選択的に表示するようにしてもよい。また、各収差マップの表示領域は、暗所時の瞳径（例えばφ６）又は、散瞳の瞳径（例えばφ８）で算出された収差マップを表示するようにしてもよい。なお、角膜のPower mapの表示領域は、変更することができる。
【００８５】
また、眼及び角膜の収差表示部４９０は、例えば、眼の収差表示については、第１の表示例で示した眼の収差表示部３２０に対応し、同じく、角膜の収差表示については、第２の表示例で示した角膜表示部４００に対応している。但し、ここでは、後述する散瞳についての数値データも表示される。眼及び角膜の収差表示部４９０は、例えば、瞳の径に応じて、明所時、暗所時、散瞳として区分した数値データを表示する。なお、ここでは、適宜の数値データが表示されるため、具体的な数値は省略した。
【００８６】
ここで、明所時とは、いわゆる明所視(photopic vision)であって、比較的高いレベルの輝度数cd/m2以上の明るさで観察している状態、生理学的には、錐体だけでものを見ている状態を意味する。暗所時とは、いわゆる暗所視（scotopic）であって、暗く低い輝度レベル輝度数（例えば、10-2cd/m2 ）で観察している状態、生理学的には、主に棹体だけでものを見ている状態を意味する。また、散瞳とは、いわゆるダイレイト(dilate)であって、一般の暗所視よりも瞳径が更に拡大した状態を意味しており、例えば、被験者が暗所時よりもさらに暗い所で長時間過ごした場合、被験者の瞳が自然に暗所時よりも拡大する状態を意味する。なお、明所時、暗所時及び散瞳に対応する瞳径は、それぞれ変更することができる。また、散瞳の表示については、表示／非表示を選択することもできる。
【００８７】
また、全収差マップ３１５は、ここでは、少なくとも２次以上の収差としたが、１次収差（チルト）を含むことで、例えば、円錐角膜等、特殊な形状を表示することができる場合も想定されるので、１次収差以上を全収差として表示してもよい。
【００８８】
（第６の表示例）
図１４は、表示部２３０にグラフィック表示される第６の表示例を示す説明図である。
表示部２３０は、例えば、眼の収差と角膜の収差とを選択的に表示したものであって、ここでは、説明の便宜上、眼の収差と角膜の収差表示の２つのパターンを兼用して示している。眼の収差表示としては、上述の第1の表示例で示した撮影した生画像であるハルトマン像３００（ここでは、Ｉｍａｇｅ３８５と表示し、詳細な画像は、省略する）と、３次以上の収差を含む高次収差３１６と、高次収差３１６に含まれる各高次収差を非対称収差と対称収差とにそれぞれ分けて表示したコマ様収差マップ３１７、球面様収差マップ３１８と、球面成分に対応する波面収差Ｓのマップ５００と、非点収差（乱視成分）に対応する波面収差Ｃのマップ５１０とを含む。
【００８９】
また、角膜の収差としては、上述のハルトマン像３００の代わりに、第２の表示例で示したプラチド測定の表示であるプラチド像３８０（ここでは、Ｉｍａｇｅ３８５と表示し、詳細な画像は、省略する）を表示する以外は、眼の収差表示と略同一であるため説明を省略する。なお、眼及び角膜の収差表示部４９５には、第５の表示例で示した眼及び角膜の収差をまとめて表示する眼及び角膜の収差表示部４９０に対して、眼の収差、又は、角膜の収差が選択的に表示される。
【００９０】
ここで、球面成分に相当する波面収差Ｓのマップ５００は、例えば、第２駆動部２６０により第1受光光学系２０が移動された移動量に対応するディオプター値（Ｓｍｏｔｏｒ）と、得られたハルトマン像３００から求めた球面成分に対応するゼルニケ係数をディオプター値に換算した値とに基づいて、算出される。すなわち、
Ｓ＝Ｓｍｏｔｏｒ−２Ｃ２１／（ｒｃ）２
となる。ここで、Ｓｍｏｔｏｒ：第１及び２駆動部２５０、２６０の移動量に対応するディオプター値、ｒｃ：解析してマップに表示する瞳径である。
一方、この波面収差Ｓのマップ５００を表示する場合には、第1受光光学系２０が移動された移動量に対応するディオプター値（Ｓｍｏｔｏｒ）をゼルニケ係数に換算し、これにハルトマン像３００から求めた球面成分に対応するゼルニケ係数を加えることで算出する。
【００９１】
図１５は、第６の表示例に対する変形例を示す説明図である。
ここでは、第６の表示例での表示部２３０に対して、球面成分に対応する波面収差Ｓのマップ５００、及び非点収差（乱視成分）に対応する波面収差Ｃのマップ５１０の代わりに、第５の表示例で示した全収差マップ３１５と、全収差マップ３１５及び第６の表示例で示した高次収差マップ３１６から算出した点像強度分布（ＰＳＦ）５２０とを表示させると共に、高次収差マップ３１６、コマ様収差マップ３１７、及び球面様収差マップ３１８の表示位置を変更してそれぞれ表示した例を示している。
【００９２】
また、点像強度分布（ＰＳＦ）５２０の表示は、３次以上の高次収差の影響を受けた点像強度分布（ＰＳＦ）を表示するモードと、２次以上を含む全収差の点像強度分布（ＰＳＦ）を表示するモードとを選択することができる。なお、眼及び角膜の収差表示部４９５は、第６の表示例と同一であり、説明を省略する。
【００９３】
（第７の表示例）
図１６は、表示部２３０にグラフィック表示される第７の表示例を示す説明図である。
表示部２３０は、眼の収差と角膜の各収差を個別に表示するものであって、眼及び角膜の収差に関しては、その各収差を、例えば、球面成分、非点成分、コマ収差、矢状収差、残余収差としてそれぞれ表示させる。表示部２３０には、具体的には、第６の表示例で示した球面成分に対応する波面収差Ｓのマップ５００、非点収差（乱視成分）に対応する波面収差Ｃのマップ５１０、コマ収差マップ５２５、矢状収差マップ５３０、３次球面収差マップ５４０、５次非点収差マップ５５０、５次球面収差マップ５６０、残余収差マップ５７０がそれぞれ表示される。なお、ここでの各波面収差の表示形態は、第１の表示例で示した眼の収差マップと同じ形態であり、その詳細な表示は、省略する。
【００９４】
ここで、図１８に戻って説明すると、波面収差Ｃのマップ５１０は、３次非点収差の項目６０７に対応している。この波面収差マップ５１０の近傍に表示されている、ゼルニケ収差係数（Ｃ＝）、方向（Ａｘ＝）には、項目６０７に記述されているそれぞれの値が表示される。コマ収差マップ５２５は、３次コマ収差の項目６０４に対応している。このコマ収差マップ５２５の近傍に表示されている、収差係数の平均２乗誤差（ＲＭＳ＝）、方向（Ａｘ＝）には、項目６０４に記述されているそれぞれの値が表示される。
矢状収差マップ５３０は、矢状収差の項目６０５に対応している。この矢状収差マップ５３０の近傍に表示されている、収差係数の平均２乗誤差（ＲＭＳ＝）、方向（Ａｘ＝）には、項目６０５に記述されているそれぞれの値が表示される。３次球面収差マップ５４０は、３次球面収差の項目６０１に対応している。この３次球面収差マップ５４０の近傍に表示されている、収差係数の平均２乗誤差（ＲＭＳ＝）には、項目６０１に記述されている値が表示される。
【００９５】
５次非点収差マップ５５０は、５次非点収差の項目６０８に対応している。この５次非点収差マップ５５０の近傍に表示されている、収差係数の平均２乗誤差（ＲＭＳ＝）、方向（Ａｘ＝）には、項目６０８に記述されているそれぞれの値が表示される。５次球面収差マップ５６０は、５次球面収差の項目６０２に対応している。この５次球面収差マップ５６０の近傍に表示されている、収差係数の平均２乗誤差（ＲＭＳ＝）には、項目６０２に記述されている値が表示される。
残余収差マップ５７０は、高次収差として利用したゼルニケ多項式の各項以外の項（具体的には、Ｚ４０、Ｚ４４、Ｚ５０、Ｚ５１、Ｚ５２、Ｚ５３、Ｚ５４、Ｚ５５、Ｚ６０、Ｚ６１、Ｚ６２、Ｚ６４、Ｚ６５、Ｚ６６等）による収差の総和を意味する。
【００９６】
（第８の表示例）
図１７は、表示部２３０にグラフィック表示される第８の表示例を示す説明図である。
表示部２３０は、眼の収差と角膜の各収差を選択的に表示するものであって、例えば、表示部２３０の上段には、ハルトマン像３００、全収差マップ３１５、乱視＋高次収差マップ５８０、高次収差マップ３１６が表示され、下段には、球面収差（各次数を含む）５８５、コマ収差（各次数を含む）５９０、高次非点収差（但し、乱視成分を除く）５９５、矢状収差マップ５３０が表示される。
乱視＋高次収差マップ５８０は、例えば、第５の表示例で示した高次収差マップ３１６と、第６の表示例で示した乱視成分に対応する波面収差のマップ５１０とに基づいて表示される。球面収差５８５は、例えば、上述した次数の各球面収差に基づいて表示される。コマ収差５９０は、例えば、上述した次数の各コマ収差に基づいて表示される。高次非点収差５９５は、例えば、上述した次数の各高次非点収差に基づいて表示される（但し、乱数成分に対応する３次非点収差は除く）。
【００９７】
（第９の表示例）
図２２は、表示部２３０にグラフィック表示される第９の表示例を示す説明図である。
表示部２３０では、例えば、角膜のレフラクティブパワーマップ(Refractive Power Map)６１０と、眼の屈折力分布(Ocular Refractive Power Map)６２０と、これらの表示に関する数値データ表示部６３０とを表示する。この数値データ表示部６３０の表示要素としては、例えば、明所時、暗所時についての平均パワー、球面収差、コマ収差、非点収差、矢状収差、残差収差、総計が含まれる。なお、角膜のレフラクティブパワー６１０と、眼の屈折力分布６２０の表示形態としては、第１の表示例で詳述した眼の収差マップと同じ形態であり、その詳細な表示例は省略する。なお、単位は、Ｄ（ディオプター）である。また、これらの表示は、必要に応じて、モノクロ、又は、カラーコードマップとすることもできる。
【００９８】
（第１０の表示例）
図２３は、表示部２３０にグラフィック表示される第１０の表示例を示す説明図である。
表示部２３０では、例えば、眼の屈折力分布(Ocular Refractive Power Map)６１０と、眼の収差マップ（高次収差）６４０と、これらの表示に関する数値データ表示部６５０とを表示する。このこの数値データ表示部６５０の表示要素としては、例えば、明所時、暗所時についての平均パワー、球面収差、コマ収差、非点収差、矢状収差、残差収差、総計が含まれる。なお、眼の屈折力分布(Ocular Refractive Power Map)６１０と、眼の収差マップ（高次収差）６４０との表示形態としては、第１の表示例で示した眼の収差マップと同じ形態であり、その詳細な表示例は省略する。また、これらの表示は、必要に応じて、モノクロ又はカラーコードマップとすることもできる。
【００９９】
（第１１の表示例）
図２４は、表示部２３０にグラフィック表示される第１１の表示例を示す説明図である。
表示部２３０では、例えば、眼球内部の光学系屈折力分布(Ocular Internal Optics  Refractive Power Map)６６０と、この表示に関する数値データ表示部６７０を表示する。この数値データ表示部６７０の表示要素としては、例えば、明所時、暗所時についての平均パワー、球面収差、コマ収差、非点収差、矢状収差、残差収差、総計を含む。なお、ここでの眼球内部の光学系屈折力分布６６０の表示形態としては、第１の表示例で示した眼の収差マップと同じ形態であり、その表示例は省略する。また、眼球内部の光学系屈折力分布６６０は、例えば、第９の表示例で示した眼の屈折力分布６２０から角膜のレフラクティブパワーマップ６１０との差として求めることができる。なお、ここでの表示は、必要に応じて、モノクロ又はカラーコードマップとすることもできる。
【０１００】
（第１２の表示例）
図２５は、表示部２３０にグラフィック表示される第１２の表示例を示す説明図である。
表示部２３０では、例えば、眼球内部の光学系屈折力分布(Ocular Internal Optics  Refractive Power Map)６６０と、眼球内部の収差分布(Ocular Internal Optics Aberration  Map)６７５と、この表示に関する数値データ表示部６８０を表示する。この数値データ表示部６８０には、図中に示すように、各種の表示要素を含む。また、ここでの分布の表示形態としては、第１の表示例で示した眼の収差マップと同じ形態であり、その表示例は省略する。
また、眼球内部の光学系屈折力分布６６０は、第１１の表示例で示した分布と同様である。また、眼球内部の収差分布(Ocular Internal Optics Aberration  Map)は、第４の表示例で示した眼内収差４４０と同様である。なお、ここでの表示は、必要に応じて、モノクロ又はカラーコードマップとすることもできる。
【０１０１】
（第１３の表示例）
図２６は、表示部２３０にグラフィック表示される第１３の表示例を示す説明図である。
表示部２３０では、例えば、角膜のアキシャルパワーマップ(Axial Power Map)６８５と、インスタンテーニアスパワーマップ(Instantaneous Power Map)６９０と、この表示に関する数値データ表示部６９５とを表示する。この数値データ表示部６９５の表示要素としては、明所時、暗所時、角膜全体として、例えば、平均パワー、球面収差、コマ収差、非点収差、矢状収差、残差収差、総計を含む。
また、ここでの角膜のアキシャルパワーマップ(Axial Power Map)６８５と、インスタンテーニアスパワーマップ(Instantaneous Power Map)６９０の表示形態は、第１の表示例で示した眼の収差マップと同じ形態であり、その表示例は省略する。また、ここでの表示は、必要に応じて、モノクロ又はカラーコードマップとすることもできる。なお、上述の各数値データ表示部に表示される明所時、暗所時、角膜全体を解析するφの値は、適宜設定することができる。
【０１０２】
このように、本実施の形態による光学特性測定装置１００では、例えば、複数の条件下で求めた測定データ（測定結果）、測定結果に対応する画像データ及び／又は数値データを、必要に応じてまとめて、又は、選択的に、表示部２３０上にグラフィック表示することができる。なお、上述の各表示例で示したハルトマン像や収差図には、すべて前眼部像、解析に用いることのできたスポット重心位置、このスポット重心位置に対応する参照格子点を、それぞれ重ね合わせて表示することもできる。
【０１０３】
【発明の効果】
本発明によると、以上説明した通り、複数の条件下で求めた測定データ（測定結果）、測定結果に対応する画像データ及び／又は数値データをまとめて、又は、選択的に表示することができる。
また、本発明によると、被測定眼の全体、角膜、眼内に対して、測定データ（測定結果）、測定結果に対応する画像データ及び／又は数値データをまとめて、又は、選択的にそれぞれ表示することができる。
【図面の簡単な説明】
【図１】本発明に関する光学特性測定装置１００の概略光学系を示す図。
【図２】本発明に関する光学特性測定装置１００の電気的構成を示す電気系ブロック図。
【図３】本発明に関する光学特性測定装置１００のフローチャート。
【図４】ランドルト環の表示に関するフローチャート。
【図５】表示部２３０にグラフィック表示される第１の表示例を示す説明図。
【図６】第１の表示例に対する変形例を示す説明図（１）、（２）。
【図７】第1の表示例に対する変形例を示す説明図（３）。
【図８】表示部２３０にグラフィック表示される第２の表示例を示す説明図。
【図９】表示部２３０にグラフィック表示される第３の表示例を示す説明図。
【図１０】表示部２３０にグラフィック表示される第４の表示例を示す説明図。
【図１１】第４の表示例に対する変形例を示す説明図。
【図１２】角膜収差測定を示すフローチャート。
【図１３】表示部２３０にグラフィック表示される第５の表示例を示す説明図。
【図１４】表示部２３０にグラフィック表示される第６の表示例を示す説明図。
【図１５】第６の表示例に対する変形例を示す説明図。
【図１６】表示部２３０にグラフィック表示される第７の表示例を示す説明図。
【図１７】表示部２３０にグラフィック表示される第８の表示例を示す説明図。
【図１８】ゼルニケ係数Cｉｊによる各収差の表現形式６００を示す図。
【図１９】ゼルニケの多項式Ｚｉｊの極座標表示による収差への分類を示す図。
【図２０】ゼルニケの多項式ＺｉｊのＸＹ座標表示による収差への分類を示す図。
【図２１】形状、Ｐｏｗｅｒ算出についての説明図。
【図２２】表示部２３０にグラフィック表示される第９の表示例を示す説明図。
【図２３】表示部２３０にグラフィック表示される第１０の表示例を示す説明図。
【図２４】表示部２３０にグラフィック表示される第１１の表示例を示す説明図。
【図２５】表示部２３０にグラフィック表示される第１２の表示例を示す説明図。
【図２６】表示部２３０にグラフィック表示される第１３の表示例を示す説明図。
【符号の説明】
１０  第1照明光学系
１１、３１、５１、５５  第1〜４光源部
１２、３２、３４、４４、５２、５３  集光レンズ
２０  第1受光光学系
２１  コリメートレンズ
２２  ハルトマン板
２３、３５、５４  第1〜３受光部
３０  送受光光学系
３３、４３、４５  ビームスプリッター
４０  共通光学系
４１  プラチドリング
４２  アフォーカルレンズ
５０  調整用光学系
６０  被測定眼
１００  光学特性測定装置
【図１】
0004652558.tif 000016
【図２】
0004652558.tif 000017
【図３】
0004652558.tif 000018
【図４】
0004652558.tif 000019
【図５】
0004652558.tif 000020
【図６】
0004652558.tif 000021
【図７】
0004652558.tif 000022
【図８】
0004652558.tif 000023
【図９】
0004652558.tif 000024
【図１０】
0004652558.tif 000025
【図１１】
0004652558.tif 000026
【図１２】
0004652558.tif 000027
【図１３】
0004652558.tif 000028
【図１４】
0004652558.tif 000029
【図１５】
0004652558.tif 000030
【図１６】
0004652558.tif 000031
【図１７】
0004652558.tif 000032
【図１８】
0004652558.tif 000033
【図１９】
0004652558.tif 000034
【図２０】
0004652558.tif 000035
【図２１】
0004652558.tif 000036
【図２２】
0004652558.tif 000037
【図２３】
0004652558.tif 000038
【図２４】
0004652558.tif 000039
【図２５】
0004652558.tif 000040
【図２６】
0004652558.tif 000041