(57)【特許請求の範囲】
【請求項１】
眼鏡レンズを掛けたときに見える網膜像を、シミュレーションによって表示する、眼鏡の視野画像表示装置であって、
前記眼鏡レンズを掛けた状態での眼球の視線の方向に対応する視野内の元画像のデータに、前記眼鏡レンズにおける前記視線の通過点に対応する、少なくとも、眼球の視軸回りの回転運動である融像回旋を含む眼球回旋量を加える処理を行うことにより、前記網膜像を作成する画像処理部と、
前記画像処理部で作成した前記網膜像を表示する表示部とを有する
ことを特徴とする眼鏡の視野画像表示装置。
【請求項２】
前記画像処理部は、視線方向の遠用アイポイントにおける基準位置からの視線移動角度量（眼球の水平垂直運動）と、前記眼鏡レンズの屈折力及び／又はプリズムとに応じて、眼球のトーション（融像回旋）角度を算出することを特徴とする請求項１に記載の眼鏡の視野画像表示装置。
【請求項３】
前記画像処理部は、前記融像回旋の角度θＬ（左眼）及びθＲ（右眼）を、リスティング法則から算出した両眼の視線方向に直交するベクトルである、ベクトルＹＬ（左眼）及びベクトルＹＲ（右眼）と、左眼の視線方向の単位ベクトル及び右眼の視線方向の単位ベクトルの外積であるベクトルＹＣとから、下記の式により算出することを特徴とする請求項１に記載の眼鏡の視野画像表示装置。
0005632245.tif 000025
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、眼鏡レンズを通して見える網膜像を、シミュレーションによって表示する、眼鏡の視野画像表示装置に係わる。
【背景技術】
【０００２】
従来、眼鏡店においては、眼鏡使用者が、サンプルレンズを使用して見え方を確認して、レンズやフレームを選択して注文していた。
しかし、眼鏡店において用意できるサンプルレンズの種類には限りがある。特に、累進屈折力レンズは、レンズの種類が非常に多いために、サンプルレンズの中に眼鏡使用者に適したレンズがあるとは限らない。そのため、注文したレンズが出来上がって実際に掛けて見るまでは、どのように見えるかがわからない。
【０００３】
そこで、眼鏡を掛けたときに見える網膜像を、シミュレーションによって表示することにより、サンプルレンズにない種類のレンズを使用した場合の見え方を示すことが提案されている（例えば、特許文献１参照。）。
累進屈折力レンズでは、度数の異なる遠用部と近用部とを、滑らかにつないでいるために歪曲収差が生じ、顔の方向を変化させると、揺れを感じる。眼鏡を掛けたときに見える網膜像を、シミュレーションによって表示することにより、この揺れ等を容易に示して説明することが可能になる。
【０００４】
前記特許文献１に記載された構成では、顔の方向の変化に基づいて、画像の範囲を変化させて、その範囲の画像に対して歪曲収差の画像処理を行っている。
【０００５】
また、眼球は、水平及び垂直運動の他、前後軸を回る回旋運動を行っている。
そして、回旋運動も含んで、眼球の運動を計測する方法が提案されている（例えば、特許文献２や、非特許文献１及び非特許文献２を参照）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３８９３７６０号明細書
【特許文献２】特開２００５−６６３５８号公報
【非特許文献】
【０００７】
【非特許文献１】坂下祐輔、藤吉弘亘、平田豊、「角膜における屈折を考慮した3次元眼球運動計測」、第１３回画像センシングシンポジウム予稿集、２００７年、IN2-25
【非特許文献２】坂下祐輔、藤吉弘亘、平田豊、「画像処理による３次元眼球運動計測」、実験力学、Ｖｏｌ．６，Ｎｏ．３，２００６年、ｐ．２３６−２４３
【発明の概要】
【発明が解決しようとする課題】
【０００８】
軸方向の存在する乱視屈折力、プリズム、非点収差等は、矯正する眼球に合わせて方向を配置する必要がある。
ところが、特定のケースにおいて、眼球が視線方向を軸とする回転である融像回旋（トーション）を行うことがある。
そして、この融像回旋（トーション）が行われる場合を調べたところ、輻輳等の視線方向が平行状態でなくなった場合に、発生することが明らかになった。
【０００９】
この輻輳に伴う融像回旋は、個人差があるものの、特に、乱視量の大きな処方値（例えば、乱視量２．０Ｄ以上）であると、その影響が顕著に発現する。
そのため、眼鏡の融像回旋（トーション）が行われている場合には、従来のシミュレーション方法では、シミュレーションの結果と現実との差異が生じて、問題となることがあった。
【００１０】
上述した問題の解決のために、本発明においては、眼球の融像回旋にも対応して、シミュレーションにより網膜像を表示することが可能である、眼鏡の視野画像表示装置を提供するものである。
【課題を解決するための手段】
【００１１】
本発明の眼鏡の視野画像表示装置は、眼鏡レンズを掛けたときに見える網膜像を、シミュレーションによって表示する構成であって、眼鏡レンズを掛けた状態での眼球の視線の方向に対応する視野内の元画像のデータに、眼鏡レンズにおける視線の通過点に対応する、少なくとも、眼球の視軸回りの回転運動である融像回旋を含む眼球回旋量を加える処理を行うことにより、網膜像を作成する画像処理部と、この画像処理部で作成した網膜像を表示する表示部とを有するものである。
【００１２】
前記本発明の眼鏡の視野画像表示装置において、画像処理部が、視線方向の遠用アイポイントにおける基準位置からの視線移動角度量（眼球の水平垂直運動）と、眼鏡レンズの屈折力及び／又はプリズムとに応じて、眼球のトーション（融像回旋）角度を算出する構成とすることができる。
【００１３】
また、前記本発明の眼鏡の視野画像表示装置において、画像処理部が、一例として、融像回旋の角度θＬ（左眼）及びθＲ（右眼）を、リスティング法則から算出した両眼の視線方向に直交するベクトルである、ベクトルＹＬ（左眼）及びベクトルＹＲ（右眼）と、左眼の視線方向の単位ベクトル及び右眼の視線方向の単位ベクトルの外積であるベクトルＹＣとから、下記の式により算出する構成とすることができる。
【数１】
0005632245.tif 000002
【００１４】
この例は、縦方向ホロプター（vertical horopter）を左眼視線方向と右眼視線方向の外積方向と仮定した場合の融像回旋角度の計算方法である。実際の縦方向ホロプターは、必ずしも左眼視線方向と右眼視線方向の外積方向とは限らない。その場合、左右眼それぞれの融像回旋角は、眼球の（網膜の）縦方向と縦方向ホロプター方向とが同一平面にあるように決定すべきである。
【発明の効果】
【００１５】
上述の本発明の眼鏡の視野画像表示装置によれば、画像処理部において、眼鏡レンズを掛けた状態での眼球の視線の方向に対応する視野内の元画像のデータに、眼鏡レンズにおける視線の通過点に対応する、少なくとも、眼球の視軸回りの回転運動である融像回旋を含む眼球回旋量を加える処理を行うことにより、網膜像が作成される。
これにより、顔の方向の変化や視線の方向の変化だけでなく、眼球の融像回旋も反映させて、眼球回旋量を加える処理をした網膜像の画像を作成することができる。そのため、融像回旋が生じている場合でも、注視部分を正確にシミュレートすることができる。
従って、本発明の眼鏡の視野画像表示装置により、融像回旋にも対応して、実際の見え方に近い網膜像を表示することが可能な眼鏡の視野画像表示装置を実現することができる。
【図面の簡単な説明】
【００１６】
【図１】本発明の一実施の形態の眼鏡の視野画像表示装置（表示システム）の概略構成図（ブロック図）である。
【図２】図１のシステムにおいて使用する、３次元ＣＧモデルの一例である。
【図３】図１の装置における、シミュレーションの画像を表示するまでの過程を示すフローチャートである。
【図４】本発明の一実施の形態のシミュレーションにおいて使用する座標系である。
【図５】レンズの屈折による光線の歪みを説明する図である。
【図６】ＰＳＦを求めるための光線追跡を説明する図である。
【図７】Ａ、Ｂ  入射瞳の分割方法を示す図である。
【図８】網膜上の結像位置と入射角との対応を示す図である。
【図９】融像回旋を説明するための座標系である。
【図１０】Ａ、Ｂ  両眼視の場合の左右の眼の視線方向を説明する図である。
【図１１】両眼共通のｙ方向を示す図である。
【図１２】等確率楕円を示す図である。
【図１３】本発明の他の実施の形態の眼鏡の視野画像表示装置（表示システム）における、シミュレーションの画像を表示するまでの過程を示すフローチャートである。
【発明を実施するための形態】
【００１７】
以下、発明を実施するための最良の形態（以下、実施の形態とする）について説明する。
なお、説明は以下の順序で行う。
１．本発明の一実施の形態の装置・システムの構成の説明
２．動画像シミュレーションの理論及び手法の説明
２−１．シミュレーションの目的
２−２．シミュレーションに使用する座標系
２−３．レンズの歪みについての説明
２−４．レンズによるボヤケについての説明
２−５．融像回旋（トーション）についての説明
２−６．シミュレーション像の合成
２−７．視野の中心視線を指定した場合の座標変換
２−８．光線データのスプライン補間近似
２−９．ＰＳＦの簡略化
２−１０．多次元Ｂスプライン補間計算の高速化
２−１１．まとめ
３．本発明の他の実施の形態の説明
４．変形例
【００１８】
＜１．本発明の一実施の形態の装置・システムの構成の説明＞
本発明の一実施の形態として、眼鏡の視野画像表示装置（表示システム）の概略構成図（ブロック図）を、図１に示す。
【００１９】
図１に示すシステムでは、ＨＭＤ（ヘッドマウントディスプレイ）１１と、ＰＣ（パーソナルコンピュータ）１５と、２つのモニター１６，１７と、入力装置である、ゲームパッド１８又はキーボード１９とを備えている。
ＨＭＤ１１には、頭部の動きセンサー（ジャイロセンサー等）１２と、視線の動きセンサー（視線追跡装置等）１３とが、搭載されている。
ＰＣ１５内には、右眼の画像用のグラフィックボード２１、左眼の画像用のグラフィックボード２２、ＵＳＢ２３を備えている。また、ＨＭＤ１１とＰＣ１５との間には、ＨＭＤコントローラ１４が接続されている。グラフィックボード２１，２２を備えていることにより、ＰＣ１５は、本発明の眼鏡の視野画像表示装置の画像処理部として動作する。
【００２０】
本実施の形態の装置・システムにおいては、累進屈折力レンズの見え方を、「ユレ（歪み）」と「ぼやけ」の観点から体験させることを目的としている。
そのため、視対象となる３次元ＣＧ（コンピュータグラフィックス）ムービーに、リアルタイム演算にて「歪み」と「ぼやけ」の画像処理効果を与えて、累進屈折力レンズの見え方を再現する。
【００２１】
図１に示すシステムでは、右眼用画像と左眼用画像とを用意することにより、立体視表示が可能となる。
また、再現される視野は、ＨＭＤ１１に搭載したジャイロセンサー（頭部の動きセンサー１２）及び視線追跡装置（視線の動きセンサー１３）によって、視聴者の視線の動きに応じて、３次元ＣＧモデル内の視野を追従させる形で表示を行う。
【００２２】
ＰＣ１５では、「歪みフィルター」、並びに「ぼやけフィルター」としての画像処理が行われる。
「歪みフィルター」とは、元画像の各画素に対して、Ｂスプライン光線データベースから、像側の出力座標（Ｔ’，Ｃ’）の値を返す機構である。
「ぼやけフィルター」とは、元画像の各画素に対して、Ｂスプライン光線データベースから、ぼやけテクスチャの大きさ（σμ，σν）、及び方向係数（ρ）の値を返す機構である。
出力となる可視画像は、元画像の全ての画素を入力とし、「歪みフィルター」を通して返される出力画像での座標位置に、「ぼやけフィルター」を通して返されるぼやけテクスチャ画像を加算（重ね塗り）して生成される。
【００２３】
図１のシステムにおいて使用する、３次元ＣＧモデルの一例を、図２に示す。
空間内に、球・円柱・直方体・立方体等の形状の５個のオブジェクト３１，３２，３３，３４，３５が配置されている。眼鏡使用者４１から見て、それぞれのオブジェクト３１，３２，３３，３４，３５までの距離は様々であり、球形状のオブジェクト３１は比較的近く、オブジェクト３２やオブジェクト３４は比較的遠くなっている。
また、図２においては、３次元ＣＧモデルの配置と共に、２方向の主視線４３Ａ，４３Ｂにおいて、眼鏡レンズ５０上の視線通過点４６Ａ，４６Ｂを併せて示している。
【００２４】
この図２を参照して、視野の位置によって眼鏡レンズ５０上の視線通過点の位置が変化することを説明する。
まず、眼鏡使用者４１が、オブジェクト３２及びオブジェクト３３の方を見て、視覚視野４４Ａが見えている状態を考える。このとき、眼鏡使用者４１の頭部４２は図示しないが視覚視野４４Ａの方を向いている。そして、主視線４３Ａは、視覚視野４４Ａの視野中心４５Ａに向かっている。ここで、視覚視野４４Ａは、比較的眼鏡使用者４１からの距離が比較的遠くなっているため、眼鏡レンズ５０上の視線通過点４６Ａの位置は、中央部よりやや上にある。
次に、眼鏡使用者４１が、オブジェクト３４及びオブジェクト３５の方を見て、視覚視野４４Ｂが見えている状態を考える。このとき、眼鏡使用者４１の頭部４２は図示しないが視覚視野４４Ｂの方を向いている。そして、主視線４３Ｂは、視覚視野４４Ｂの視野中心４５Ｂに向かっている。ここで、視覚視野４４Ｂは、比較的眼鏡使用者４１からの距離が視覚視野４４Ａよりも近くなっているため、眼鏡レンズ５０上の視線通過点４６Ｂの位置は、先ほどの視線通過点４６Ａよりも下にあり、眼鏡レンズ５０の中央部付近にある。
このように、視覚視野４４Ａ，４４Ｂが変わることにより、眼鏡使用者４１の頭部４２の向きが変わるだけではなく、視覚視野４４Ａ，４４Ｂまでの距離等に対応して、眼鏡レンズ５０上の視線通過点４６Ａ，４６Ｂの位置も変わってくる。
このようなことから、前記特許文献１のように頭部の動きに対応させるだけで、視線通過点が固定であると、実際の見え方とは大きく異なってくることになる。
本実施の形態では、図１に示したシステムを使用して、眼鏡レンズ５０上の視線通過点４６Ａ，４６Ｂの位置の変化にも対応するように、シミュレーションの画像を表示する。
【００２５】
次に、シミュレーションの画像を表示するまでの過程を、図３のフローチャートに示す。
以下、図３を参照して、図１のシステムにおけるシミュレーションの過程を説明する。
【００２６】
まず、ＨＭＤ１１に搭載されたジャイロセンサー（頭部の動きセンサー１２）や視線追跡装置（視線の動きセンサー１３）により、眼鏡使用者４２の顔の位置、向き（頭部４２の向き）、視線方向を検出する。ステップＳ１では、右眼用の検出を行い、ステップＳ２では、左眼用の検出を行う。
一方、ステップＳ３において、図２に示したように、ＣＧ仮想物体（３次元ＣＧモデル）を用意する。
次に、ステップＳ４において、ＣＧウォークスルー機能を使用して、ステップＳ１及びステップＳ２によって検出した顔の位置や向きや視線方向によって、３次元ＣＧモデルから切り取る視覚視野を求める。この視覚視野は、右眼用と左眼用とでそれぞれ異なる視野となる。
次に、ステップＳ５において、ステップＳ４で切り取った視覚視野から、右眼用の歪み、ボヤケが全くない元画像を作成する。同様に、ステップＳ６において、ステップＳ５において、ステップＳ４で切り取った視覚視野から、左眼用の歪み、ボヤケが全くない元画像を作成する。
一方、ステップＳ７において、入力装置（キーボード１９等）により、眼鏡使用者４１の処方度数、加入度、レンズ種類を入力する。
次に、入力した内容から、ステップＳ８で右眼レンズの形状データ及びレイアウトデータ及び右眼球モデルを作成し、ステップＳ９で左眼レンズの形状データ及びレイアウトデータ及び左眼球モデルを作成する。
次に、ステップＳ８〜Ｓ９でそれぞれ作成した、形状データ、レイアウトデータ、眼球モデルに基づいて、ステップＳ１０において、光線データの３次元スプライン補間係数を生成する。
次に、ステップＳ１０で生成した光線データの３次元スプライン補間係数を用いて、ステップＳ１１において、右眼用の、出射光線方向、ＰＳＦパラメータ、レンズ通過点位置、その他いろいろなパラメータの３次元スプライン補間係数を求める。同様に、ステップＳ１２において、左眼用の、出射光線方向、ＰＳＦパラメータ、レンズ通過点位置、その他いろいろなパラメータの３次元スプライン補間係数を求める。
次に、ステップＳ５及びステップＳ６で作成した元画像と、ステップＳ１１及びステップＳ１２で求めたパラメータや補間係数等を使用して、ステップＳ１３において、シミュレーションを実行する。このシミュレーションの過程では、画像処理用ハードウエアを使用することも含んでいる。
次に、ステップＳ１４において、右眼用の、歪み及びボヤケを含む画像を作成する。同様に、ステップＳ１５において、左眼用の、歪み及びボヤケを含む画像を作成する。
このようにして作成した、歪み及びボヤケを含む画像は、ＨＭＤ１１の表示画面や、右眼監視用モニター１６、左眼監視用モニター１７に表示する。
以上説明した過程によって、ＨＭＤ１１の表示画面において、視線方向に対応した、歪みやボヤケを含む画像が表示される。
【００２７】
＜２．動画像シミュレーションの理論及び手法の説明＞
２−１．シミュレーションの目的
この動画像シミュレーションは、眼鏡をかけた場合にどのように感じるかを、静止画像や動画像で表現することを目的としている。
３次元ＣＧとＨＭＤ１１、ジャイロセンサー、視線追跡装置を組み合わせて使うことにより、仮想空間の中で頭を回したり、視線を変更したりするときに見える画像を、リアルタイムで提示することが可能である。
また、図１のシステム及び図３のフローチャートに示したように、左右の眼にそれぞれ別々の画像を提示すれば両眼立体視も可能である。
【００２８】
２−２．シミュレーションに使用する座標系
本実施の形態のシミュレーションに使用する座標系を、図４に示す。
図４に示すように、ｘ軸、ｙ軸、ｚ軸の３軸の直交座標系を構成する。ｘ軸は、正面から目に入る方向にとる。ｙ軸は、ｘ軸と直交し、上に向く方向とする。ｚ軸は、右から左に向く水平方向とする。ｘ軸−ｙ軸−ｚ軸の方向は、右手則に則っている。そして、原点は、眼球の回旋中心に置く。図中の破線は、眼球とその角膜とを模式的に示している。
【００２９】
そして、図４に示す座標系において、空間内の任意一点Ｐ（ｘ，ｙ，ｚ）（ここではｘ＜０、つまり目の前）は、目に入る光線の角度β，γと回旋中心点までの距離ＰＯとによって、表し直すことができる。
シミュレーション画像上の位置は、縦方向ではｔａｎβ＝ｙ／ｘとなり、横方向ではｔａｎγ＝ｚ／ｘとなる。物体距離は、そのまま表すのではなく、逆数で表わすことが、眼鏡の場合には都合がいい。従って、空間上任意点の位置は、下記のように表わすことができる。
【数２】
0005632245.tif 000003
【００３０】
２−３．レンズの歪みについての説明
レンズを通して見ると、光線が屈折されてしまう。
つまり、裸眼では（ψ，ζ）の方向にある物体点が、眼鏡レンズを通して見ると、（ψ'，ζ'）に移動する。
このことを、図５を参照して、さらに詳しく説明する。図５には、眼鏡レンズ５０と、眼球５１と、眼鏡レンズ５０の凹面に対応する後方頂点球面５２を示している。
図５に示す任意点Ｐは、裸眼時の入射方向がＰＯであるのに対し、眼鏡をかけて眼鏡レンズ５０を通して見ると、眼球５１の回旋中心Ｏへの入射方向がＲＯに変わる。
同様に、図５に示す点Ａは、裸眼時の入射方向がＡＯであるのに対し、眼鏡レンズ５０を通して見ると、眼球５１の回旋中心Ｏへの入射方向がＢＯに変わる。
【００３１】
ここで、眼鏡をかけたときに感じる物体の位置（ψ'，ζ'）は裸眼の時の位置（Ｄｌ，ψ，ζ）の関数で表すことができる。つまり、
ψ'＝ψ'（Ｄｌ，ψ，ζ）
ζ'＝ζ'（Ｄｌ，ψ，ζ）
という関数で表すことができる。この関数の中身は、後述する光線追跡によって確定することができる。
なお、融像回旋がある場合には、中心視線の融像回旋角に合わせて座標変換する必要がある。融像回旋の原理については後述する。
【００３２】
２−４．レンズによるボヤケについての説明
レンズによるボヤケの原因は、物体点からのすべての光線が網膜の一点に収束しないことである。
物体点からの光が、像点を中心にある範囲に広がる光量分布を形成する。この分布を、点広がり関数（Point Spread Function；ＰＳＦ）という。
【００３３】
この関数ＰＳＦを求める方法を、図６を参照して説明する。
ＰＳＦを求めるときには、まずその点Ｐを通る主光線ＰＱＯを探す。
主光線が決まれば、入射瞳を均等に分割して（例えば４００分割）、Ｐ点から各分割領域の中心と結ぶ光線を追跡し、網膜に交わる点を求める。
図６において、Ｐ点の方位角と、Ｐ点出発の光線の主光線に対する方位角とを用いて、各光線を追跡することができる。
入射瞳の位置は、厳密に言うと瞳孔の物体側の共役点だが、主光線の物体側の線ＰＱの延長線上の一点Ｏ’で、ＰＯ＝ＰＯ’を満たすところに設定しても差し支えない。
【００３４】
幾何光学原理のＰＳＦであれば、網膜交点の密度がそのままＰＳＦである。
波動光学の効果を考える場合には、さらに各分割領域の光線の光路差を計算して、フレネル積分してＰＳＦを求める。
【００３５】
次に、入射瞳の分割方法は、様々な方法が考えられる。
主な分割方法として、図７Ａに示す正方分割と、図７Ｂに示す螺旋分割との２種類が挙げられる。
図７Ａに示す正方分割は、縦横に領域を分割して、各領域の中心点を使用している。この場合には、単純明解だが、四隅に無駄な部分があり、予定の光線数の７割程度しか追跡できない。
一方、図７Ｂに示す螺旋分割は、入射瞳の中心点から螺旋状に延びる曲線の上の点を使用している。この場合には、予定の光線数の全てを追跡できる。
なお、螺旋状の配置は少なくとも２以上の螺旋を含むことが望ましい。２以上の螺旋を含むことにより、１つの螺旋のみを含む配置よりも、入射瞳を効率良く利用することができる。図７Ｂでは６つの螺旋を含んでおり、このように６つの螺旋を含んだ場合には、最も効率良く入射瞳を利用することができる。
【００３６】
こうして求めたＰＳＦは、網膜面上の光密度分布であるが、入射画像の座標は眼球の回旋中心から見た方向の座標（ψ，ζ）であるため、網膜上の座標と入射画像の座標（ψ，ζ）との変換が必要である。
ここで、入射角と像高の関係を、図８に示す。ＰＳＦの有効範囲では像高が小さく、近軸計算で十分精度が得られると考える。
つまり、ψｍ＝ｙｍ／ｆ，ζｍ＝ｚｍ／ｆとなる。なお、ｆは目の焦点距離であり、処方度数によって変化する。
【００３７】
このように、網膜上位置の光量分布は、入射光線方向の光量分布に変換することができる。つまり、物体点からの光は、物体点だけではなく物体点を中心としたある範囲の空間からも来ているように目が感じている。
そして、近傍の点と点の間は、お互いに影響を及ぼしあって、区別がつきにくく、ぼやけるように見える。
【００３８】
レンズ上異なる位置を通して見た場合のＰＳＦは当然異なる。
また、たとえレンズ上同じ位置を通して見ていても、物体距離が異なっていれば、ＰＳＦも異なる。
さらにまた、レンズ上同じ位置を通して同じ物体距離の点を見ても、目の調節状態が異なるとＰＳＦも異なる。
なお、融像回旋がある場合には、眼球とレンズとの位置関係が融像回旋角だけ変わって、ＰＳＦも変化する。特に、乱視処方（眼球の度数に乱視成分が含まれる）の場合、注意が必要である。
【００３９】
２−５．融像回旋（トーション）についての説明
前述したように、視線方向が平行状態でなくなった場合には、眼球が融像回旋を行うことがある。
遠方を見る場合の眼球の回旋は、リスティング法則（Listing’s Law）に基づく。このリスティング法則は、眼球が空間上のある方向に向いたときの姿勢を定める法則である。
眼球の姿勢というのは、眼球の横方向と縦方向の向きを指す。眼球の姿勢が定まらないと、網膜像の上下左右が定まらない。
ただし、視線の方向、即ち、眼球の光軸の方向が決まっただけでは、眼球の姿勢は視線を軸に回転する全ての方向を取りうる。
【００４０】
リスティング法則は、無限遠方の任意の視線方向に向けた眼球の姿勢を定める。
リスティング法則について、例えば、「視覚情報処理ハンドブック」ｐ．４０５では、「片眼のどんな回転も１つの平面（リスティング平面）内の軸を中心にして起こるとみなしうる」と記載されている。
【００４１】
このことについて、図９に示す座標系を用いて説明する。
図９に示す座標系は、図４に示した座標系と同様に、Ｘ軸は正面（水平前方）から目に入る方向、Ｙ軸はＸ軸と直交する上方向、Ｚ軸は水平方向、としている。図９のＸ軸方向は、第１眼位とも呼ばれる。
そして、眼球と、眼球の回旋中心の点Ｒを示している。Ｙ−Ｚ平面は、前述したリスティング平面である。
【００４２】
任意方向への眼球回旋後の姿勢は、点Ｒを含むリスティング平面内の直線を軸にした回転と同じである。図９においては、この回転軸となる直線の一例を、Ｙ軸とＺ軸との間に記載している。そして、この回転軸は、第１眼位（Ｘ軸方向）とも、回転後の視線方向とも、垂直である。
仮に、（図示しない）方向ベクトル（Ｌ，Ｍ，Ｎ）へ眼球回旋するならば、回旋後の眼球座標系のＸ軸、Ｙ軸、Ｚ軸のそれぞれのベクトルは、下記の式（１）で計算される。
【００４３】
【数３】
0005632245.tif 000004
【００４４】
リスティング法則は、片眼が無限遠方の物体に対して眼球の姿勢を定めるものとしては正しい。
また、例えば、無限遠方の物体を視ていて体を傾けた場合には、左眼と右眼とで眼球の姿勢が同じで、眼球の回旋も同様である。
これに対して、無限遠方でない物体を両眼で視る場合には、左眼と右眼とで眼球の姿勢が異なることがある。
ここで、両眼視の場合の左右の眼の視線方向を説明する図を、図１０Ａ及び図１０Ｂに示す。
【００４５】
両眼視の場合、無限遠方物体ならば、図１０Ａに示すように、左眼の眼球５１Ｌと右眼の眼球５１Ｒとを同じ視方向に向けるので、左右眼球がそれぞれリスティング法則に従って回旋した後の姿勢も同一である。
このとき、左右眼それぞれの網膜像に差異は生じない。
なお、図１０Ａ及び図１０Ｂでは、左眼の眼球５１Ｌと右眼の眼球５１Ｒとの間の中央部に、左眼と右眼とを平均した眼球５５を破線で示している。
【００４６】
一方、有限距離の物体（点Ａ）に対しては、図１０Ｂに示すように、両眼輻輳が必要になる。
そのため、左眼の眼球５１Ｌの視方向と、右眼の眼球５１Ｒの視方向とが異なり、眼球の回旋量が、左右の眼球５１Ｌ，５１Ｒで異なる。図１０Ｂでは、点Ａが左前方にあるため、右眼の眼球５１Ｒの回旋量が左目の眼球５１Ｌの回旋量よりも多くなっている。
【００４７】
リスティング法則に基づく眼球回旋では、回旋後の眼球姿勢、つまり、回旋後のｙ軸とｚ軸の各方向ベクトルは、式（１）に示した視方向ベクトルに依存している。
左眼と右眼の視方向ベクトルが異なれば、回旋後のｙ軸とｚ軸の各方向ベクトルは、左右眼で一致せず、網膜像の回転ずれが起こる。
この回転ずれを解消するために、左右の眼において、それぞれ視線周りの回旋が必要になる。この回旋を、融像回旋という。
【００４８】
融像回旋の角度（融像回旋角）を決めるためには、左右の眼に共通のｙ方向（縦方向）を決める必要がある。
この共通のｙ方向は、左眼の視線と右眼の視線とに対して、共に垂直である方向が自然である。つまり、図１０Ｂの左目の眼球５１Ｌの中心と点Ａとの間の直線と、右眼の眼球５１Ｒの中心と点Ａとの間の直線とを含む平面に垂直な方向である。
従って、図１１に両眼共通のｙ方向を示すように、左眼の視線方向のベクトルｘＬと右眼の視線方向のベクトルｘＲとに、それぞれ垂直な方向に、共通のｙ方向のベクトルｙＣがある。この方向は、理論的縦方向ホロプターともいう。
【００４９】
仮に、左右の眼の視線方向単位ベクトルｘＬとｘＲがそれぞれ、下記の式（２）の通りであるとする。
【数４】
0005632245.tif 000005
このとき、共通のｙ方向のベクトルｙＣは、下記の式（３）で表される。
【数５】
0005632245.tif 000006
【００５０】
左右の眼がそれぞれの視線方向に向けて回旋したときの、リスティング法則に基づいたｙ方向のベクトルｙＬ，ｙＲは、下記の式（４）で表される。
【数６】
0005632245.tif 000007
ベクトルｙＬとベクトルｙＣとの角度が、左眼球の融像回旋角度となり、同様にベクトルｙＲとベクトルｙＣとの角度が、右眼球の融像回旋角度となる。
【００５１】
さらに、眼鏡をかける場合は、眼鏡レンズの屈折作用を考慮しなければならない。
この場合、共通のｙ方向ｙＣ自体は変わらないが、左右それぞれのレンズを通過した後のｙ方向を光線追跡で求めて、光線追跡で求めた共通のｙ方向と左右の眼のリスティング法則に基づいたｙ方向のベクトルとの角度を求めると、融像回旋角が求められる。
【００５２】
光線追跡で求めた共通のｙ方向のベクトルＹＣと、左右の眼のｙ方向のベクトルＹＬ，ＹＲとから、下記の式により、左右の眼球の融像回線角度θＬ，θＲを算出することができる。
【数７】
0005632245.tif 000008
【００５３】
なお、以上の式は両眼共通の縦方向、つまり、縦方向ホロプターが、左右眼視線の外積方向である場合の融像回旋角を求める式である。実際の縦方向ホロプターは、左右眼視線の外積方向とは異なる場合がある。その場合の融像回旋角は、眼球の縦方向と実際の縦方向ホロプター方向が同一平面内にあるように決めるべきである。
【００５４】
２−６．シミュレーション像の合成
以上、レンズによる歪みとボヤケと融像回旋について説明した。
歪みとボヤケとを画像処理の手法で合成すると、眼鏡レンズをかけて見える画像のシミュレーションが可能である。さらに、静止画像だけではなく、動画像も可能である。
歪みについては、像側視野内のすべての画素の対応点物体側に求めて、元画像の輝度情報を当てることで実現できる。
ボヤケについては、各画素の輝度をそのＰＳＦに基づいて周辺画素に“分配”して、画像全画素の輝度を再構成することによって実現できる。
ボヤケ処理は、たたみこみ演算（Convolution）とも呼ばれる。一般的なたたみこみ演算と異なるのは、ＰＳＦが一定でないことである。
融像回旋（トーション）については、視野内の元画像のデータに、融像回旋を含む眼球回旋量を加える処理を行うことによって実現できる。
【００５５】
より好ましくは、融像回旋（トーション）について、視線方向の遠用アイポイントにおける基準位置からの視線の移動角度量（眼球の水平方向及び垂直方向の運動）と、眼鏡レンズの屈折力及び／又はプリズムとに応じて、眼球の融像回旋（トーション）の角度を算出する。
なお、視線方向の遠用アイポイントにおける基準位置とは、図１０Ｂでは、左右の眼球５１Ｌ，５１Ｒから伸びる破線の方向であり、左眼と右眼とを平均した眼球５５からでは、正面から角度φだけ左にある。図１０Ｂでは、左右の眼の移動角度量の和が角度αとなっている。
【００５６】
さらに好ましくは、融像回旋の角度θＬ（左眼）及びθＲ（右眼）を、リスティング法則から算出した両眼の視線方向に直交するベクトルである、ベクトルＹＬ（左眼）及びベクトルＹＲ（右眼）と、左眼の視線方向の単位ベクトル及び右眼の視線方向の単位ベクトルの外積であるベクトルＹＣとから、前述した式により算出する。
【００５７】
なお、融像回旋（トーション）については、左右の眼の視線方向が平行な場合には生じないため、被写体の位置によって左右の眼の視線方向が平行となっているときには、融像回旋の回旋量を加える処理は行わない。
融像回旋（トーション）は、無限遠方を見ている遠方視は除く近方視のときに生じるものであり、左右の眼の視線方向が異なる（平行でなく、交差する）ときに生じるものである。
また、注視している物体が近づいて来る際に眼球が内転する運動（即ち、輻輳）や、注視している物体が遠ざかる際に眼球が外転する運動（即ち、開散）に伴い、融像回旋（トーション）の処理が加えられる。
【００５８】
２−７．視野の中心視線を指定した場合の座標変換
以上のように、真正面から目に入る方向をｘ軸とした座標系（以下、グローバル座標系と呼ぶ）において、歪み情報（物体側視方向から像側視方向への変換）と、ボヤケ情報（特定視方向、視距離におけるＰＳＦ）とを用いれば、視野全体のシミュレーションが可能である。
【００５９】
しかしながら、実際のシミュレーションは、必ずしも中心視線が真正面向きではない。
例えば、近用部の見え方を確認したい場合は、中心視線がレンズ近用部を通っている必要がある。この場合、中心視線はグローバル座標の斜め方向に当たる。この斜め方向をｘ’軸とするローカル座標系において、シミュレーションを行う。
このとき、ローカル座標系のｙ’軸とｚ’軸をいかに決めるかが問題となる。ここでは眼球回旋の法則のひとつである、リスティング法則に則って決める。リスティング法則によれば、中心視線が真正面向きのときの上下方向及び左右方向が、眼球が回旋して中心視線が移動することに対応して、ある決まった方向にそれぞれ変化する。そして、中心視線が移動して視線方向が変化したときにも、実際の物体の上下方向及び左右方向が網膜像でも上下方向及び左右方向となるように変化する。
座標軸の変換マトリックスは、下記の式（５）になる。
【数８】
0005632245.tif 000009
（ここで、ａ，ｂ，ｃは、中心視線方向のグローバル座標における方向単位ベクトル（ａ  ｂ  ｃ）の各軸成分である。）
そして、グローバル座標の任意一点（ｘ，ｙ，ｚ）のローカル座標（ｘ'，ｙ'，ｚ'）は、下記の式（６）で変換される。
【数９】
0005632245.tif 000010
逆に、ローカル座標の任意一点（ｘ'，ｙ'，ｚ'）のグローバル座標（ｘ，ｙ，ｚ）は、下記の式（７）で変換される。
【数１０】
0005632245.tif 000011
【００６０】
以上のような座標変換式を用いれば、レンズ上の任意点を通過する視線方向を中心視線とする場合の歪みを、リアルにシミュレートすることが可能となる。
【００６１】
２−８．光線データのスプライン補間近似
メガネレンズを通しての見え方をシミュレーションする光学的原理及び画像処理手法は、上述のようにして確立した。
しかし、いざシミュレーションを始めると、膨大な計算量に悩まされる。眼鏡レンズの形状は単純な球面ではなく、特に累進レンズは自由曲面である。
眼鏡レンズのような複雑な面の光線追跡には、繰り返して収束する手法を採用する。これは単純な球面の光線追跡よりも、少なくとも数倍の時間がかかる。
また、シミュレーションする画像の画素数の多さも、光線追跡計算量の上昇に拍車をかけている。
全てのシミュレーションの結果画像の全画素に対して、元画像のどの画素に対応しているか探すための光線追跡（主光線追跡）を実行しなければならない。その上に、ＰＳＦを確定するために対応物体点から発する多数（例えば１００）の光線を追跡して、網膜上のスポットを求める。これらの光線追跡が全て非球面に対応するために、繰り返し収束手法を採用しなければならないので、途方もなく膨大な計算負担になるのである。
現在のパーソナルコンピュータ１台の計算能力では、このような手法では、一枚の画像（動画像の一コマ）を処理するのに数日かかることになる。
【００６２】
一方、レンズ形状、眼球パラメータ、レンズと眼球の位置関係が全部確定している状態では、下記の物体側の任意一点に対し、像側（ψ'，ζ'）は唯一に決まる。
【数１１】
0005632245.tif 000012
つまり、下記の関数
ψ'＝ψ'（Ｄｌ，ψ，ζ）
ζ'＝ζ'（Ｄｌ，ψ，ζ）
が成立する。
しかも、変数（Ｄｌ，ψ，ζ）に対し、（ψ'，ζ'）は連続変化することも、容易に想像できる。このような関数は、スプライン補間に適している。
【００６３】
そこで、各変数の定義域内に、有限個のサンプル点を設定する。例えば、物体距離逆数Ｄｌのサンプル点を、（−０．２，０．０，０．２，０．５，０．８，１．１，１．４，１．７，２．０，２．３，２．６，２．９，３．２，３．６，４．０）の１５点とし、上下角度のタンジェントψのサンプルを、（−１．５，−１．２，−１．０，−０．８，−０．６，−０．４，−０．２，０．０，０．２，０．４，０．６，０．８，１．０，１．２，１．５）の１５点とし、左右角度のタンジェントζのサンプルを（−１．５，−１．２，−１．０，−０．８，−０．６，−０．４，−０．２，０．０，０．２，０．４，０．６，０．８，１．０，１．２，１．５）の１５点に設定する。
これらのサンプル点のすべての組み合わせについて、光線追跡を実行し、真の関数値を求める。
それ以外の変数値（サンプル点とサンプル点の間の値）に対する関数値は、サンプル点における真の関数値を用いて補間する方法が確立されている。補間方法は、目的に合わせて多数あるが、サンプル点における真値が既知の場合の補間方法としては、Ｂスプライン法が最も適している。
【００６４】
サンプル点の数や間隔は、補間精度に関係する。一般的に、サンプル点間隔の小さいところの補間精度は高い。ただし、間隔を小さくすると、定義域全部をカバーするためのサンプル数が多くなり、プログラム上メモリを多く確保する必要がある。最近のＰＣやＯＳはメモリを沢山装着できるようになったので、サンプル数の制限は緩和され、高精度の結果が得られるようになった。
【００６５】
このようにして、歪み情報を表す関数
【数１２】
0005632245.tif 000013
を、少ない計算量のスプライン補間で得られるようになる。式中、Ｃは補間係数であり、Ｎは各次元の節点に基づいた基底多項式関数である。
【００６６】
２−９．ＰＳＦの簡略化
前述したように、厳密にある物体点のＰＳＦを取得するには、物体点から発し、入射瞳を均等分割する多数の点を通過する光線を追跡し、網膜上のスポットを求めて、さらにスポット密度分布関数を求める。
しかし、この方法では、光線の数を増やしても精度の向上が思うほど上がらない。
また、収差が少ない場合にはスポットが集中して、像点以外にはほとんど光線が通らないケースもあれば、度数誤差の大きい場合には、ある領域で均等分布のケースもあり、ＰＳＦ関数の変化が激しい。
一方、シミュレーションやレンズパフォーマンス評価の場合、必ずしも正確なＰＳＦを必要としない。例えば、視力の場合、二つの点が分別できる最も近い距離（視角）を表している。この場合、ＰＳＦの関数の精密な形は必要とせず、ＰＳＦがカバーする範囲の大きさが重要なパラメータである。従って、ＰＳＦを大胆に簡略しても、レンズ性能評価に果たす役割が大きく影響されることはないと言える。
逆に、ＰＳＦをあらかじめ連続関数に仮定して、そのパラメータを光線追跡のデータを用いて当てはめるようにすれば、少ないパラメータでＰＳＦを表現できる。それらのパラメータを、前述の歪み関数のようにスプライン補間(三次元)で求めることができる。
【００６７】
簡略関数の形状は、度数誤差やあらゆる軸角度の非点収差のＰＳＦを近似できるように、二次元正規分布が適切と考えられる、つまり、簡略関数は、下記の式（８）のようになる。
【数１３】
0005632245.tif 000014
（ここで、μ，νはそれぞれｙ、ｚ方向の主光線からの偏移角、σμ，σν，ρは正規分布のパラメータである。これらのパラメータはσμ＞０，σν＞０，−１＜ρ＜１を満たす。）
下記の式（９）で表される楕円の線上すべての点（μ，ν）において、次の式（１０）が成り立つ。
【数１４】
0005632245.tif 000015
【数１５】
0005632245.tif 000016
そして、その等高線楕円内の積分は、下記の式（１１）となる。
【数１６】
0005632245.tif 000017
この場合の等確率楕円を、図１２に示す。
【００６８】
このように、二次元正規分布関数は、広がりの程度(σμ，σν)と非点ボケの程度（等確率楕円長短軸比)、角度(長軸の角度)を表すことができる。
もちろん、ＰＦＳの光学系の状態による無限に近い変化を忠実に表すことはできないが、ＰＳＦを表現する簡略関数としては有効である。
【００６９】
二次元正規分布関数のパラメータσμ，σν，ρを、光線データから求める方法を考えると、（μ,ν）平面に散布する多数の光線の交点（各交点が入射瞳上の各分割点に対応）の統計値を求めて、σμ，σν，ρに当てる方法が考えられる。つまり、
【数１７】
0005632245.tif 000018
である。ここで、Ｎは光線数で、(μｉ，νｉ)は交点座標である。
【００７０】
このように、物体空間上の任意一点（Ｄｌ，ψ，ζ）のＰＳＦ分布関数を、パラメータσμ，σν，ρを持つ二次元正規分布関数で近似することができる。さらに、σμ，σν，ρは（Ｄｌ，ψ，ζ）の関数と表すことができる。つまり、
σμ＝σμ（Ｄｌ，ψ，ζ）
σν＝σν（Ｄｌ，ψ，ζ）
ρ＝ρ（Ｄｌ，ψ，ζ）
これらの関数は、歪み情報と同様にスプライン補間で求めることもできる。つまり、下記のように求めることができる。
【数１８】
0005632245.tif 000019
【００７１】
ここで注意したいのは、スプライン補間誤差で関数値が定義域を超える可能性がある、ということである。例えば、−１＜ρ＜１だが、補間で求めるとρ＝１．００２のような結果になって、楕円が存在しない結果になることがある。この問題の解決法は、ρの代わりにｓｉｎ−１ρを補間して求め、得られた結果に対してｓｉｎ演算を行って、ρを求めることが有効である。
【００７２】
歪みやボヤケのパラメータ以外にも、有用なパラメータをスプライン補間で求めることができる。例えば、主光線のレンズ凸面通過点位置（ｙｃｏｎｖｅｘ，ｚｃｏｎｖｅｘ）、凹面通過点位置（ｙｃｏｎｃａｖｅ，ｚｃｏｎｃａｖｅ）等が挙げられる。これらのパラメータは、下記のようにして、算出することが可能である。
【数１９】
0005632245.tif 000020
【００７３】
主光線のレンズ通過点位置は、レンズ特定位置の透過光線を中心視線とするローカル座標系においての歪みやボヤケの分布のシミュレーションに役立つ。
【００７４】
２−１０．多次元Ｂスプライン補間計算の高速化
１次元スプライン補間は、
【数２０】
0005632245.tif 000021
で表される。ここで、ｉは各次元の節点番号、Ｃｉはその係数、ｎは標本点数である。Ｎｉ（ｘ)は、ｉ番節点に対応する基底関数であり、階数Ｍの場合、ｉ番節点とｉ＋Ｍ番節点との間の範囲でゼロでない値を持ち、隣接節点間はｍ−１次多項式で表される（基底関数の局部性のため）。
言い換えると、ｘの定義域内の任意点ａにおいては、ゼロでない値のＮｉ（ｘ）が最多でもＭ個しか存在しない。
従って、補間式は一見するとｎ項あるように見えるが、ｘ＝ａにおいては実質Ｍ項であり、Ｍ回の掛け算とＭ回の足し算でＦ（ａ）が得られる。
【００７５】
３次元スプライン補間は、下記式（１２）で表される。
【数２１】
0005632245.tif 000022
ここで、ｉ，ｊ，ｋは各次元の節点番号であり、それぞれ標本点数だけ変化する。
つまり、項の数は各次元の標本点数の積になる。しかし、上述の基底関数の局部性により、ある一点については、ゼロでない項の数は、各次元の階数の積である。
各次元のスプライン階数が４の場合、項の数は４３＝６４である。つまり、一回の補間演算では、足し算を６４回、掛け算を６４×３＝１９２回行うことになる。
一般的には、ｎｊ次元のＭ階スプライン補間演算に必要な掛け算の回数は、ｎｊ×Ｍｎｊであり、次元数が大きくなるにつれて、急激に計算負担が増える。
【００７６】
ところが、上記式（１２）を下記式（１３）のように書き直すと、計算回数を若干減らすことができる。
【数２２】
0005632245.tif 000023
この式（９）は、１次元の補間のネスト構造（入れ子構造）であり、次元の順番は自由に変えることができる。掛け算と足し算は、共に４＋４×(４＋４×４)＝８４回であり、ほぼ１／２の計算時間で済む。
一般的には、ｎｊ次元のＭ階スプライン補間演算に必要な掛け算の回数は、下記式（１４）で表される。
【数２３】
0005632245.tif 000024
【００７７】
２−１１．まとめ
３次元ＣＧを利用した動画像シミュレーションの、主に技術的側面について説明した。累進レンズをかけて見えた３次元の外部世界を動画像でシミュレーションするため、計算量は膨大である。
この膨大な計算量を、スプライン補間とＰＳＦのパラメータ化によって削減することができ、実現可能なレベルまで達成することが可能になる。
そして、本発明を適用しないでシミュレーションの計算を行った場合と比較して、本発明を適用した場合には、計算に要する時間を１／１００〜１／１０００程度にまで短縮することができる。
【００７８】
実際に、左右両眼の歪みやボヤケパラメータのスプライン係数データベース事前に準備して、高性能パソコンと画像グラフィックボードを用いてシミュレーション画像を生成した。
そして、ジャイロを搭載したＨＭＤで、ＣＧの中でリアルタイムウォークスルーを実行したところ、１０フレーム／秒の速さで実現することができた。
【００７９】
３．本発明の他の実施の形態の説明
本発明では、上述の実施の形態のように、ＨＭＤ１１とモニター１６，１７の両方への画像表示を行う構成に限定されるものではない。モニターを省略してＨＭＤのみに表示する等、いずれか一方の表示部のみを設けて、その表示部に網膜像を表示する構成としてもよい。
本発明において、網膜像を表示する表示部は、ＨＭＤやモニターに限定されない。また、頭部の動きのセンサーは、ジャイロセンサーに限定されるものではない。
【００８０】
ここで、本発明の他の実施の形態として、眼鏡の視野画像表示装置（表示システム）の構成を説明する。
本実施の形態では、図１及び図３に示したＨＭＤ１１とモニター１６，１７の代わりに、３Ｄ（立体）ディスプレー、例えば、３Ｄテレビジョンや３Ｄプロジェクター等を表示部として使用する。そして、あらかじめ動画を作成した後に、一コマずつシミュレーション画像を作成して、このシミュレーション画像から眼鏡レンズをかけた場合の動画像を作成して、左右眼画像を３Ｄディスプレーに表示する。
【００８１】
本実施の形態の装置における、シミュレーションの画像を表示するまでの過程を、図１３のフローチャートに示す。
まず、ステップＳ２１において、図２に示したと同様に、ＣＧ仮想物体（３次元ＣＧモデル）を用意する。
また、ステップＳ２２において、あらかじめ作成した動画像ストーリーを用意する。
次に、ステップＳ２３において、ＣＧ仮想物体（３次元ＣＧモデル）と動画像ストーリーとから、動画の一コマを左右眼の視点及び視方向から作成する。
次に、ステップＳ２４において、ステップＳ２３で作成した動画の一コマから、右眼用の歪み、ボヤケが全くない元画像の各画素の輝度（ＲＧＢ）と視点からの距離を求める。同様に、ステップＳ２５において、ステップＳ２３で作成した動画の一コマから、左眼用の歪み、ボヤケが全くない元画像の各画素の輝度（ＲＧＢ）と視点からの距離を求める。
一方、ステップＳ２６において、入力装置（キーボード１９等）により、眼鏡使用者の処方度数、加入度、レンズ種類を入力する。
次に、入力した内容から、ステップＳ２７で右眼レンズの形状データ及びレイアウトデータと右眼球モデルを作成し、ステップＳ２８で左眼レンズの形状データ及びレイアウトデータ及び左眼球モデルを作成する。
次に、ステップＳ２７〜Ｓ２８でそれぞれ作成した、形状データ、レイアウトデータ、眼球モデルに基づいて、ステップＳ２９において、光線データの３次元スプライン補間係数を生成する。
次に、ステップＳ２９で生成した光線データの３次元スプライン補間係数を用いて、ステップＳ３０において、右眼用の、出射光線方向、ＰＳＦパラメータ、レンズ通過点位置、その他いろいろなパラメータの３次元スプライン補間係数を求める。同様に、ステップＳ３１において、左眼用の、出射光線方向、ＰＳＦパラメータ、レンズ通過点位置、その他いろいろなパラメータの３次元スプライン補間係数を求める。
次に、ステップＳ２４及びステップＳ２５で作成した元画像の各画素の輝度及び視点からの距離と、ステップＳ３０及びステップＳ３１で求めたパラメータや補間係数等を使用して、ステップＳ３２において、シミュレーションを実行する。このシミュレーションの過程では、画像処理用ハードウエアを使用することも含んでいる。
次に、ステップＳ３３において、右眼用の、歪み及びボヤケを含む画像を作成する。同様に、ステップＳ３４において、左眼用の、歪み及びボヤケを含む画像を作成する。
このようにして作成した、歪み及びボヤケを含む画像は、３Ｄディスプレーにおいて立体表示する。
以上説明した過程によって、３Ｄディスプレーの表示画面において、視線方向に対応した、歪みやボヤケを含む立体画像が表示される。
【００８２】
本実施の形態において、上述した説明以外のその他の構成は、先の実施の形態と同様であるので、重複説明を省略する。
【００８３】
４．変形例
上述の実施の形態では、歪みを加える処理では、像側視野内のすべての画素の対応点物体側に求めて、元画像の輝度情報を当てていた。
また、上述の実施の形態では、ボヤケを加える処理では、各画素の輝度をそのＰＳＦに基づいて周辺画素に“分配”して、画像全画素の輝度を再構成していた。さらにまた、式（１）で表わされる正規分布関数とパラメータを使用していた。
本発明では、歪みやボケを加える処理の方法は、上述した実施の形態で説明した方法に限定されるものではなく、その他の方法を使用することが可能である。
そして、本発明は、歪みを加える処理とボケを加える処理とのうち、一方又は両方に、上述した実施の形態で説明した方法以外の方法を使用する構成も含むものである。
【００８４】
また、上述の実施の形態では、右眼用画像と左眼用画像とをそれぞれ作成して立体視が可能な構成としていた。
本発明は、右眼用画像或いは左眼用画像のいずれか一方の眼用の画像のみを作成する構成も含む。
【００８５】
また、上述の実施の形態では、本発明を累進屈折力レンズのシミュレーションを行う場合に適用して、融像回旋を含む眼球回旋量を加える処理を行い、網膜像を作成する構成であった。
本発明では、累進屈折力レンズに限らず、他の眼鏡レンズにおいても、融像回旋を含む眼球回旋量を加える処理を行って、網膜像を作成することが可能である。
【００８６】
本発明は、上述の実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲でその他様々な構成が取り得る。
【符号の説明】
【００８７】
１１  ＨＭＤ、１５  ＰＣ、１６，１７  モニター、４１  眼鏡装用者、４２  頭部、４３Ａ，４３Ｂ  主視線、４４Ａ，４４Ｂ  視覚視野、４６Ａ，４６Ｂ  視線通過点、５０  眼鏡レンズ
【図１】
0005632245.tif 000026
【図２】
0005632245.tif 000027
【図３】
0005632245.tif 000028
【図４】
0005632245.tif 000029
【図５】
0005632245.tif 000030
【図６】
0005632245.tif 000031
【図７】
0005632245.tif 000032
【図８】
0005632245.tif 000033
【図９】
0005632245.tif 000034
【図１０】
0005632245.tif 000035
【図１１】
0005632245.tif 000036
【図１２】
0005632245.tif 000037
【図１３】
0005632245.tif 000038