JP3979904B2 - レンズの形状予測方法 - Google Patents
レンズの形状予測方法 Download PDFInfo
- Publication number
- JP3979904B2 JP3979904B2 JP2002265411A JP2002265411A JP3979904B2 JP 3979904 B2 JP3979904 B2 JP 3979904B2 JP 2002265411 A JP2002265411 A JP 2002265411A JP 2002265411 A JP2002265411 A JP 2002265411A JP 3979904 B2 JP3979904 B2 JP 3979904B2
- Authority
- JP
- Japan
- Prior art keywords
- lens
- coordinate system
- point
- curved surface
- reference coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
- Eyeglasses (AREA)
Description
【発明の属する技術分野】
本発明は、レンズの形状予測方法と、その方法を用いたレンズの性能表示方法及び眼光学系のシミュレーション方法と、それら方法の実施に使用する装置及びプログラムとに関する。
【0002】
【従来の技術】
従来から、レンズを通して見たときに外界の空間がどのように見えるかを、光線追跡法等を用いてシミュレーションする方法が知られている(例えば、特許文献1参照)。このようなシミュレーションは、製品としての実レンズの見え方を購入前に確認するだけでなく、設計段階で仮想的なレンズの見え方を確認しつつより見易いレンズを設計すること等にも役立てられている。ところで、見え方のシミュレーションを行うには、レンズ面の形状データが必要になる。例えば、シミュレーションの演算を高速かつ正確に行う為に、スプライン補間法を用いることがあるが、その場合には、広範にわたる四角形領域(例:80mm×80mm)の形状データが必要となる。通常は、設計段階で少なくともレンズ面の全領域における形状データを求めるから、その予め求めた形状データに基づいてシミュレーションを行うことができる。
【特許文献1】
特開2000-107129号公報(第6−13頁、第11図)
【0003】
【発明が解決しようとする課題】
一方、形状データの不明なレンズについて見え方のシミュレーションをしたい場合には、3次元形状測定装置を用いて、そのレンズ面の形状を直接測定することにより、該レンズについての形状データを得ることになる。ところが、3次元形状測定装置を用いる場合には、レンズ面の全領域における形状データが分かる訳ではなく、所定の限られた領域についてしか形状データが分からない。
【0004】
3次元形状測定装置として市販の松下電器社製:UA3Pを用い、外径60φの丸レンズを測定したときの測定可能領域を図9に示す。同図において、符号Lはレンズであり、符号2は、UA3Pによる測定可能領域を示す。そして、この測定可能領域2内にある複数の黒点P,P…が、それぞれUA3Pによって3次元座標が測定されたサンプリング点である。このように、UA3Pを用いた場合には、レンズL面上における八角形の限られた領域2内についての形状データしか分からない。
【0005】
ここで、この場合に得られる形状データについて、図10に示す基準座標系を想定して説明する。この座標系は、任意の基準面をYZ平面とし、該基準面からの高さ(距離)をX座標とするものである。八角形領域2内の各サンプリング点Pにおいては、UA3PによってYZ平面上の位置(Y,Z)と、基準面からの高さ(X)とが測定される。つまり、この種の3次元形状測定装置によれば、各サンプリング点Pについての基準座標系でのXYZ座標が測定される。測定可能領域2内の全てのサンプリング点PについてのXYZ座標値が、当該領域2の形状データである。
【0006】
しかし、この測定可能領域2の外においては、形状データが分からないので、光線追跡法を用いることができない。従って、この領域2内の形状データだけでは、適切なシミュレーションを行うのに充分であるとは言い難い。しかも、スプライン補間法を用いてシミュレーションを行う場合には、四角形領域の形状データが必要となるから、この領域2に内接する四角形領域4(サイズ例:40mm×40mm)内の形状データだけを用いることになってしまい、データ不足がさらに深刻化する。
【0007】
なお、見え方のシミュレーション以外にも、形状データの不足が問題化する場面がある。例えば、何らかの理由によって形状データが不明となってしまったマスターレンズについて、そのレンズと同一形状のレンズを製作したい場合である。つまり、加工機を用いる場合には、レンズ面に加工を施す工具の運動軌跡の関係上、該レンズ面の全領域の形状データだけでは不充分であり、該レンズ面を延長した面領域、即ち実際にはレンズ面の無い領域(例えば、図9中、符号6で示す領域)内の形状データまでもが必要となる。一般には、レンズの設計段階で、その延長領域6の部分まで形状データを求めておくのであるが、形状データが不明となったレンズについては、仮に3次元形状測定装置を用いたとしても、そのような延長領域6の形状データはもはや未知である。
【0008】
そこで本発明は、レンズの形状データが不足している場合に、該レンズについての既知の形状データに基づいて、必要な範囲内の形状データを求める技術を提供することを目的とする。また本発明は、レンズ面上における一部の領域内の形状データしか分かっていないレンズについても、該レンズについての見え方のシミュレーションを適切に行う技術を提供することを目的とする。また本発明は、レンズの形状データが不足している場合であっても、該レンズと同一のレンズを製作できるようにすることを目的とする。
【0009】
【課題を解決するための手段】
本発明によれば、基準面をYZ平面とし、該基準面からの高さをX座標とする基準座標系でのXYZ座標に相当する値が予め分かっている既知点を、レンズ面上に複数有するレンズの形状予測方法であって、前記既知点のうち、前記基準座標系でのX座標に相当する値が分かっていない未知点の近傍に位置するものを、近傍点として複数特定する近傍点特定工程と、特定した前記各近傍点についての前記基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい曲面を求める曲面算出工程と、求めた前記曲面を定義する式に、前記未知点の前記基準座標系でのYZ座標に相当する値を代入することにより、該未知点の前記基準座標系でのX座標に相当する値を求める高さ算出工程とを有するレンズの形状予測方法が提供される。
【0010】
本発明が適用されるレンズとしては、球面レンズであってもよいが、球面でないレンズに適用して好適である。球面でないレンズとしては、例えば、累進屈折力レンズや乱視レンズ等の眼鏡レンズが挙げられる。ここでは、眼鏡レンズのみならず、眼鏡レンズ以外の光学レンズも全て含めて、レンズ面が球面でないレンズのことを「非球面レンズ」と呼ぶことにする。
また前記曲面としては、2次曲面が好ましい。2次曲面としては、球面や放物面等が挙げられる。放物面としては、乱視放物面等がある。
【0011】
【作用】
未知点の近傍に位置する既知点を、近傍点として複数特定することにより、レンズ面上における近傍点の辺りの局所に良好にフィットする曲面を求めることができる。その曲面上における未知点の高さ(X座標に相当する値)を求めることにより、その高さ値は当該局所のレンズ面形状の規則性にしたがった値となる。これにより、形状が既知の限られた面領域を、該面領域の規則性にしたがって自然に拡張させたときの、その拡張領域の形状データを求めることができる。高さの求まった未知点を、既知点として予測領域を拡張してゆくことにより、必要な全範囲内の形状データを求めることができる。
レンズ曲面を外挿する場合、曲面の光学特性の連続性を保つことが大事である。この発明は、未知点の高さを求めるために用いる局所曲面として、球面、または乱視放物面を採用することで、面曲率の急激な変化を防ぐことができる。
【0012】
【発明の実施の形態】
〔第1の実施の形態〕
図1は、第1の実施の形態によるレンズ形状予測装置の構成を示す。このレンズ形状予測装置は、レンズの形状データの一部を入力する入力部10と、入力された一部の形状データから、当該レンズについての必要な全範囲内の形状データを予測する演算部12と、予測された形状データを出力する出力部14とを備える。なお、本実施の形態では、レンズとしての累進レンズに適用する。
【0013】
入力部10は、3次元形状測定装置16によって測定された形状データを取得する取得手段として機能する。該形状データは、レンズL面上における各サンプリング点(既知点)Pについての、基準座標系でのXYZ座標に相当するものである。基準座標系とは、基準面をYZ平面とし、該基準面からの高さをX座標とする座標系である(図10参照)。前述のように、3次元形状測定装置16によって測定された形状データは、レンズL面上の一部の領域2だけの形状を表すものであり、このデータのみでは、レンズの性能表示や見え方のシミュレーション等を行うに不充分である。
【0014】
演算部12は、CPU等を含んで構成されており、メモリ18に格納されている形状データ予測プログラム20を実行する。
該プログラム20は、基準座標系でのXYZ座標に相当する値が予め分かっている既知点Pをレンズ面上に複数有するレンズ(累進レンズ)Lの形状を予測するためにコンピュータに対して、既知点P,P…のうち、基準座標系でのX座標に相当する値が分かっていない未知点Qの近傍に位置するものを、近傍点として複数特定する処理と、特定した各近傍点についての基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい球面SPを求める処理と、求めた球面SPを定義する式に、未知点Qの基準座標系でのYZ座標に相当する値を代入することにより、該未知点Qの基準座標系でのX座標に相当する値を求める処理とを実行させるものである(図2(a),(b)参照)。これら一連の処理は、ループ処理である。
【0015】
従って演算部12は、この形状データ予測プログラム20を実行することにより、本発明に係る装置の近傍点特定手段、曲面算出手段(曲面算出手段)、高さ算出手段、及び既知範囲拡張手段としての機能を果たす。なお、これらの手段がプログラム等のソフトウエアによらずに、電子回路等のハードウエアによっても実現できるのは言うまでもない。また、上記形状データ予測プログラム20を、通信ネットワーク等の電気通信回線を通じて他のコンピュータに提供することにより、該コンピュータをレンズ形状予測装置として機能させることもできる。
【0016】
次に図2(a),(b)を参照して、演算部12の処理について説明する。
〔近傍点特定工程〕
まず、既知点P,P…のうち、基準座標系でのX座標に相当する値が分かっていない未知点Qの近傍に位置するものを近傍点として複数特定する。ここで、未知点Qの基準座標系におけるYZ座標の値は、任意に決定できる。そして、そのYZ座標値を有する未知点QのX座標(高さ)を求めようというのである。
【0017】
そこで、最初の未知点Qは、3次元測定装置の測定可能領域2に内接するリングR1の周上の点とする。そして、その未知点Qを含む所定サイズの微小領域Sを想定し、その微小領域S内に含まれる既知点Pを近傍点として特定する。ここでは、微小領域Sとして、所定半径の円領域を想定している。
【0018】
ここで、微小領域(円領域)Sのサイズは、小さすぎると、形状が既知の部分の規則性を充分に利用できない一方、大きすぎると却って未知点の近傍の局所における規則性が薄れてしまう。従って、両者の兼ね合いを考慮したサイズとするのが好ましい。
【0019】
〔球面算出工程〕
次いで、特定した各近傍点についての基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい球面SPを求める。
そのために、まずその球面SPの基準座標系における中心座標(曲率中心)を(XC,YC,ZC)と仮定し、半径(曲率半径)をRと仮定する。すると、該球面SPを定義する方程式は、次式(1)で表される。
(XC−X)2+(YC−Y)2+(ZC−Z)2=R2 …(1)
【0020】
そして、各近傍点からの距離が最も小さい球面SPを求めるために、(1)式のX,Y,ZにそれぞれXi,Yi,Ziを代入して、次の誤差関数f(XC,YC,ZC,R)を定義する。
【数1】
但し、Xi,Yi,Ziは、近傍点の基準座標系での3次元座標である。これは、3次元座標測定装置によって予め既知である。ここにiは、各既知点に割り当てられた自然数である。従って、微小領域S内の近傍点に割り当てられている自然数のうち最小のものをnとし、最大のものをmとすれば、iはnからmまで変化させる。
【0021】
この誤差関数は、各近傍点の球面SPからの距離の二乗和を表す。従って、この誤差関数f(XC,YC,ZC,R)の値が最小となるときのXC,YC,ZC,Rを求めれば、各近傍点からの距離の二乗和が最も小さい球面SPが求まることになる。
【0022】
なお誤差関数は、上式に限定されるものではなく、各近傍点の球面SPからの距離の偶数乗和、又は絶対値の和等を表す関数を定義してもよい。偶数乗和であれ、絶対値の和であれ、求めようとする曲面(この場合は球面)と各近傍点との離れ具合を表す関数を誤差関数と呼び、その誤差関数の値が最小となるときの曲面を総称して「各近傍点からの距離が最も小さい曲面」と表現することにする。
【0023】
誤差関数f(XC,YC,ZC,R)の最小値は、即ち該関数の極小値である。誤差関数が極小値をとるときには、この関数を構成する各変数での偏導関数がゼロとなる。そこで、次の4つの連立方程式をたてることができる。
【数2】
【0024】
変数がXC,YC,ZC,Rの4つで、方程式が4つであるから、全ての変数の値を求めることができる。具体的には、最初の式から、次式が得られる。
【数3】
但し、Nは微小領域S内にある近傍点の個数(m−n)を表す。この式を、残りの3つの式に代入すると、変数Rが消去されて、XC,YC,ZCの3元方程式となる。それを展開して整理すると、
【数4】
が得られる。
【0025】
上式において、バー(―)は平均値を表す。例えば、
【数5】
とは、微小領域S内に含まれる全ての近傍点についてのX座標値(Xi)の平均であり、次式で表される。
【数6】
【0026】
また、
【数7】
とは、
【数8】
のことである。
【0027】
以上の連立方程式を解くことで、XC,YC,ZC,Rの4つの全ての変数の値が求まる。従って、中心座標(XC,YC,ZC)と半径Rが特定されたから、各近傍点からの距離が最も小さい球面SPが求まったことになる。
【0028】
〔高さ算出工程〕
次いで、求めた球面を定義する式に、未知点Qの基準座標系でのYZ座標に相当する値を代入することにより、該未知点Qの基準座標系でのX座標に相当する値を求める。
ここで、球面を定義する式とは、上記求めたXC,YC,ZC,Rの値を(1)式に代入した式のことであり、未知点Qの基準座標系でのYZ座標は、上記近傍点特定工程で予め決定した値である。
【0029】
これにより、該未知点Qについての基準座標系でのXYZ座標が求まったことになる。次いで、その未知点Qとは別の未知点QのXYZ座標を求めるために、再び上記近傍点特定工程に戻る。即ち、近傍点特定工程から高さ算出工程までは、ループ(繰り返し)処理である。例えば、未知点Qが1000個ある場合には、近傍点特定工程から高さ算出工程までを1000回繰り返す。一回のループ毎に、微小領域SをリングRnに沿ってずらしてゆく。また繰り返し処理の過程で、一旦XYZ座標を求めた未知点Qを、既知点(近傍点)Pとして取り扱い、さらに形状データを拡張してゆくことができるのは言うまでもない。
【0030】
具体的には、リングR1の全周における未知点Q,Q…の高さ(X座標)を求めたら、今度はリングR1よりも半径の大きなリングR2の周上の未知点Q,Q…の高さを求め、さらに今度はリングR2よりも半径の大きなリングR3の周上の未知点の高さを求める。このようにして、リングのサイズもR1,R2,R3…Rnと順次拡大してゆく。このようにしてゆけば、レンズL面の全領域のみならず、該レンズ面を延長した面領域6(図8参照)についても形状データを予測できる。
【0031】
なお、リングR1,R2,R3…Rnの各々の周上における未知点Qの間隔は、微小領域としての円領域Sの半径に一致させるか、又は近似させるのが好ましい。これは、レンズ面上における既知点P及び未知点Qをなるべく均等に配分するためである。このようにすれば、該レンズについての性能表示及び見え方のシミュレーション(眼光学系のシミュレーション)をより適切に行える。
【0032】
〔性能表示工程〕
以上のようにして、必要な領域(この場合は、レンズ面の全領域)の形状データが求まったら、該形状データが出力部14から出力される。出力された形状データを用いて、図示せぬ公知のレンズの性能表示手段によって、当該レンズの性能を表示したものを図3に示す。なお、性能表示手段は公知であるからこれについての詳細な説明は省略する(特開2002−14006号公報参照)。
図3(a)は、レンズL面上の平均度数の分布を等高線でつないだ平均度数分布図であり、図3(b)は、レンズL面上の非点収差の分布を等高線でつないだ非点収差分布図である。
【0033】
〔検証〕
検証の為に、設計段階で求めた当該レンズLについての全領域6(図8参照)における正確な形状データに基づいて、当該レンズLの性能を表示したものを図4(a),(b)に示す。図3(a)と図4(a)、及び図3(b)と図4(b)は、それぞれ良く一致しており、この方法が有効であることが確かめられた。
以上のように、本実施の形態によれば、八角形の限られた測定可能領域2内の形状データのみを用いながら、レンズ全体としての性能を正確に求めることができた。
【0034】
〔比較例〕
比較のために、本発明を適用せずに、3次元形状測定装置によって測定された形状データを用いて、当該レンズLの性能を表示したものを図5に示す。これは、空白地帯(レンズL面上における測定可能領域2の外)を容認したまま、該空白地帯におけるスプライン係数を予め定めておき、測定可能領域2内におけるスプライン係数だけを最適化して、平滑化スプラインを作成したものである。図5(a)が、平均度数分布図であり、図5(b)が非点収差分布図である。
【0035】
図示のように、形状データの存在する測定可能領域2は、図4と略一致しているが、測定可能領域2の外では、レンズ面の形状が完全にコントロール不能である。測定可能領域2の外では、光線追跡法を適用できないので、見え方のシミュレーションを行うにも不充分である。従って、レンズLの形状データが不足している場合には、該レンズLについての既知の形状データに基づいて、必要な範囲内の形状データを求めることが不可欠である。
【0036】
〔第2の実施の形態〕
レンズ形状予測装置の構成は、第1の実施の形態と略同様であるため、図1と同一の符号を用いて説明する。なお本実施の形態でも、第1の実施の形態と同一の累進レンズLについて適用する。
【0037】
メモリ18に格納されている形状データ予測プログラム20は、基準面をYZ平面とし、該基準面からの高さをX座標とする基準座標系でのXYZ座標に相当する値が予め分かっている既知点Pを、レンズ面上に複数有するレンズ(累進レンズ)Lの形状を予測するためにコンピュータに対して、既知点P,P…のうち、基準座標系でのX座標に相当する値が分かっていない未知点Qの近傍に位置するものを、近傍点として複数特定する処理と、特定した各近傍点についての基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい球面SPを求める処理と、求めた球面SPの曲率中心から該球面SP上の中心Cへ向かう直線を、中心における法線とし、かつ各近傍点からの距離が最も小さい乱視放物面Paを求める処理と、求めた乱視放物面Paを定義する式に、未知点Qの基準座標系でのYZ座標を代入することにより、該未知点のX座標を求める処理とを実行させるものである。これら一連の処理は、ループ処理である。
【0038】
従って本実施の形態において、演算部12は、上の形状データ予測プログラム20を実行することにより、本発明に係る装置の近傍点特定手段、曲面算出手段たる球面算出手段及び乱視放物面算出手段、高さ算出手段、並びに既知範囲拡張手段としての機能を果たす。
【0039】
つまり、本実施の形態では、各近傍点に対して乱視放物面を当てはめる。ここで、乱視放物面について、図6(a)を参照して説明する。図示のように、乱視放物面Paには、2つの主曲率A,Bがあり、一方の主曲率Aはこの乱視放物面Paの最小曲率であり、他方の主曲率Bはこの乱視放物面Paの最大曲率である。なおこれら2つの主曲率A,Bが同じ値のときは、球面となる。
【0040】
このように、2つの主曲率A,Bを有する乱視放物面Paは、両者が交わる点である中心Oと該中心Oにおける法線(主軸)N方向とが定まったときに始めてその姿勢が定まる。つまり、中心Oと該中心Oにおける法線Nとが定まらない限り、仮に主曲率A,Bが求まったとしても、微小領域S上には無限の乱視放物面Paを想定できてしまうことになる。そこで、中心Oと法線Nとをどのようにして定めるのかという問題が生じる。この点に、乱視放物面を当てはめるときの困難性がある。この問題に対し、本実施の形態では、第1の実施の形態で求めた球面SPに基づいて、それら中心O及び法線Nを定めようというのである。
【0041】
以下、演算部12の処理について説明するが、球面算出工程までの工程は、第1の実施の形態と同じであるので、球面算出工程後の工程について説明する。
〔乱視放物面算出工程〕
本工程では、上記球面算出工程で求めた球面の曲率中心から該球面上の頂点へ向かう直線を、中心における法線とし、かつ各近傍点からの距離が最も小さい乱視放物面を曲面として求める。
【0042】
即ち、まず図2(c)において、球面SPの曲率中心(XC,YC,ZC)から該球面SPの頂点Cへ向かう直線Lを、図6(b)における当該求めようとする乱視放物面Paの中心Oにおける法線Nとする。
これは具体的には、図2(c)における球面SPの頂点Cの座標を原点O(X0,Y0,Z0)とし、かつ直線LをX軸とする座標系を定め、その座標系において、乱視放物面を定めるということである。以下、その座標系を、基準座標系(YYZ座標系)と区別するために、X'Y'Z'座標系と呼ぶ。
【0043】
次いで、X'Y'Z'座標系において、微小領域S内の各近傍点からの距離が最も小さい乱視放物面Paを曲面として求める。以下、具体的に説明する。
X'Y'Z'座標系において、X'軸を主軸(光軸)とする乱視放物面の一般式は、次式で表される。
X'=aY'2+bX'Y'+cZ'2+f …(2)
この方程式(2)において、変数は、a,b,c,fの4つである。
【0044】
そして、各近傍点からの距離が最も小さい乱視放物面Paを求めるために、(2)式によって定義される乱視放物面と各近傍点との離れ具合を表す誤差関数を定義する。該誤差関数は、次の式で表される。
【数9】
但し、Xi'とは、基準座標系(XYZ座標系)における各近傍点のX座標(=Xi)を、X'Y'Z'座標系における座標に変換した値である。
【0045】
この誤差関数Fは、各近傍点の乱視放物面Paからの距離の二乗和を表す。従って、この誤差関数Fの値が最小となるときのa,b,c,fを求めれば、各近傍点からの距離の二乗和が最も小さい乱視放物面Paが求まることになる。なお、誤差関数が上式に限定されないのは、前述の通りである。
【0046】
誤差関数Fの最小値は、即ち該関数の極小値である。誤差関数が極小値をとるときには、この関数を構成する各変数a,b,c,fでの偏導関数がゼロとなる。従って、そのことを表す4つの連立方程式をたてることができる。
つまり、
【数10】
…等が得られる。
変数が4つで方程式が4つであるから、全ての変数の値が求まることになる。
【0047】
上記4つの連立方程式を展開して整理すると、次のようになる。
【数11】
但し、Yi'とは、基準座標系(XYZ座標系)における各近傍点のY座標(=Yi)を、X'Y'Z'座標系における座標に変換した値であり、Zi'とは、基準座標系(XYZ座標系)における各近傍点のZ座標(=Zi)を、X'Y'Z'座標系における座標に変換した値である。
【0048】
上記4つの方程式を行列で表現すると、次のようになる。
【数12】
この式を解けば、a,b,c,fの4つの全ての変数の値が求まる。従って、それら変数の値を(2)式に代入することにより、乱視放物面Paが求まったことになる。
【0049】
また、その式(2)を基準座標系に変換すれば、基準座標系における乱視放物面Paの方程式を求めることができる。
例えば、X'Y'Z'座標系におけるX'Y'Z'がそれぞれ、基準座標系におけるXYZを用いて次のように表されるとする。
X'=U(X,Y,Z)
Y'=V(X,Y,Z)
Z'=W(X,Y,Z)
そうすると、それらを(2)式に代入することにより、基準座標系における乱視放物面Paの方程式を求めることができる。
【0050】
〔高さ算出工程〕
次いで、求めた乱視放物面Paを定義する式に、未知点Qの基準座標系でのYZ座標に相当する値を代入することにより、該未知点Qの基準座標系でのX座標に相当する値を求める。
ここにいう乱視放物面を定義する式とは、例えば(2)式を基準座標系に戻したときの乱視放物面Paの方程式とするが、このことは、(2)式のY'Z'に、それぞれ未知点QのX'Y'Z'座標系におけるY'Z'座標(Y'i,Z'i)を代入することにより、該未知点QのX'Y'Z'座標系におけるX'座標(X'i)を求め、求めたX'座標(X'i)をX座標(Xi)に戻すことと等価である。
【0051】
これにより、該未知点Qについての基準座標系でのXYZ座標が求まったことになる。次いで、その未知点Qとは別の未知点QのXYZ座標を求めるために、再び上記近傍点特定工程に戻り、未知点Qの個数だけ近傍点特定工程から高さ算出工程までのループを繰り返す。こうして、必要な領域(この場合は、レンズ面の全領域)の形状データが求まったら、該形状データが出力部14から出力される。
【0052】
〔性能表示工程〕
出力された形状データを用いて、図示せぬ公知のレンズの性能表示装置によって、当該レンズの性能を表示したものを図7に示す。図7(a)は、レンズL面上の平均度数の分布を等高線でつないだ平均度数分布図であり、図7(b)は、レンズL面上の非点収差の分布を等高線でつないだ非点収差分布図である。図7(a)と図4(a)、及び図7(b)と図4(b)は、それぞれ良く一致しており、この方法が有効であることが確かめられた。さらに、図7(a),(b)の方が、図3(a),(b)よりも、図4(a),(b)によく近似しており、第2の実施の形態における方法がなおいっそう有効であることが確かめられた。
【0053】
〔シミュレーション工程〕
さらに、出力部14から出力された形状データを用いて、図示せぬ公知の眼光学系のシミュレーション手段によって、当該累進レンズ(HOYA株式会社製:Summit CD)Lを通して外界を見たときの見え方をシミュレーションした結果を図8に示す。なお、眼光学系のシミュレーション手段は公知であるから、これについての詳細な説明は省略する(特開2000−107129号公報参照)。
図8(b),(d),(f)は、ぞれぞれ出力部14から出力された形状データに基づく当該累進レンズLの遠用部、中間部、近用部の見え方を表したものである。図8(a),(c),(e)はそれぞれ検証のために、設計段階で求めた当該レンズLについての正確な形状データに基づく遠用部、中間部、近用部の見え方を表したものである。
図8中、暗く現されている部分が見えにくい領域である。図8において、(a)と(b)、(c)と(d)、(e)と(f)は、それぞれ良く一致しており、この方法が有効であることが確かめられた。
【0054】
以上、本発明の実施の形態について説明したが、本発明はこれに限られない。例えば、或るループにおいて微小領域Sが見込むレンズ面が平面に近い場合、そのループにおいては、その局所には曲面(例えば、球面や乱視放物面)を当てはめるのではなく、平面を当てはめるようにしてもよい。
【0055】
また、予測した形状データは、レンズの性能表示や眼光学系のシミュレーション以外にも利用できる。例えば、領域6(図9参照)の形状データを予測し、それをレンズ加工機に入力すれば、当該レンズと同一のレンズを製作できるようになる。
【0056】
また、レンズとして眼鏡レンズに適用したが、本発明はコンタクトレンズや眼内レンズ等にも適用できる。さらに、本発明はレンズのみならず、物品の曲面等の3次元形状の予測にも広く応用できる。
【0057】
【発明の効果】
本発明によれば、レンズの形状データが不足している場合に、該レンズについての既知の形状データに基づいて、必要な範囲内の形状データを求めることができる。また本発明によれば、レンズ面上における一部の領域内の形状データしか分かっていないレンズについても、該レンズについての見え方のシミュレーションを適切に行える。また本発明によれば、レンズの形状データが不足している場合であっても、該レンズと同一のレンズを製作できるようになる。
【図面の簡単な説明】
【図1】実施の形態によるレンズ形状予測装置の構成を示すブロック図。
【図2】実施の形態によるレンズ形状予測方法を説明するための図。
【図3】一実施例によるレンズ性能表示の結果を現した図。
【図4】検証のための正確な形状データに基づくレンズ性能表示の結果を現した図。
【図5】比較例によるレンズ性能表示の結果を現した図。
【図6】実施の形態によるレンズ形状予測方法を説明するための別の図。
【図7】別の実施例によるレンズ性能表示の結果を現した図。
【図8】実施例による眼光学系のシミュレーションの結果を現した図。
【図9】3次元形状測定装置による測定可能領域を説明するための図。
【図10】基準座標系にレンズが配置されている様子を示す図。
【符号の説明】
2…測定可能領域、12…演算部(近傍点特定手段,曲面算出手段,球面算出手段,乱視放物面算出手段,高さ算出手段)、L…累進レンズ(レンズ)、P…サンプリング点(既知点)、Q…未知点、SP…球面(曲面)、Pa…乱視放物面(曲面)、C…球面上の頂点、L…球面の曲率中心から該球面上の頂点へ向かう直線、O…乱視放物面の中心、N…乱視放物面の法線(光軸)。
Claims (16)
- 基準面をYZ平面とし、該基準面からの高さをX座標とする基準座標系でのXYZ座標に相当する値が予め分かっている既知点を、レンズ面上に複数有するレンズの形状予測方法であって、
前記既知点のうち、前記基準座標系でのX座標に相当する値が分かっていない未知点の近傍に位置するものを、近傍点として複数特定する近傍点特定工程と、
特定した前記各近傍点についての前記基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい曲面を求める曲面算出工程と、
求めた前記曲面を定義する式に、前記未知点の前記基準座標系でのYZ座標に相当する値を代入することにより、該未知点の前記基準座標系でのX座標に相当する値を求める高さ算出工程と
を有するレンズの形状予測方法。 - 請求項1記載のレンズ形状予測方法において、
前記高さ算出工程の後、前記基準座標系でのX座標に相当する値が求められた未知点を、前記既知点とみなし、該既知点の近傍に位置する未知点の前記基準座標系でのX座標に相当する値を同様に求めることにより、既知点の範囲を拡張していくレンズ形状予測方法。 - 請求項1又は2記載のレンズの形状予測方法において、
前記曲面算出工程で求める前記曲面が球面であるレンズの形状予測方法。 - 請求項1又は2記載のレンズの形状予測方法において、
前記曲面算出工程は、
まず、前記各近傍点についての前記基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい球面を求める球面算出工程と、
次いで、求めた前記球面の曲率中心から該球面上の頂点へ向かう直線を、中心における法線とし、かつ前記各近傍点からの距離が最も小さい乱視放物面を前記曲面として求める乱視放物面算出工程と
からなるレンズの形状予測方法。 - 請求項1乃至4の何れか記載のレンズの形状予測方法を用いて、レンズの形状データを生成する形状データ生成工程と、
生成した形状データを用いて、当該レンズの性能を表示する性能表示工程と
を有するレンズの性能表示方法。 - 請求項1乃至4の何れか記載のレンズの形状予測方法を用いて、レンズの形状データを生成する形状データ生成工程と、
生成した形状データを用いて、当該レンズを通して外界を見たときの見え方をシミュレーションするシミュレーション工程と
を有する眼光学系のシミュレーション方法。 - 基準面をYZ平面とし、該基準面からの高さをX座標とする基準座標系でのXYZ座標に相当する値が予め分かっている既知点を、レンズ面上に複数有するレンズの形状予測装置であって、
前記既知点のうち、前記基準座標系でのX座標に相当する値が分かっていない未知点の近傍に位置するものを、近傍点として複数特定する近傍点特定手段と、
特定した前記各近傍点についての前記基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい曲面を求める曲面算出手段と、
求めた前記曲面を定義する式に、前記未知点の前記基準座標系でのYZ座標に相当する値を代入することにより、該未知点の前記基準座標系でのX座標に相当する値を求める高さ算出手段と
を有するレンズの形状予測装置。 - 請求項7記載のレンズ形状予測装置において、
前記高さ算出手段によって前記基準座標系でのX座標に相当する値が求められた未知点を、前記既知点とみなし、該既知点の近傍に位置する未知点の前記基準座標系でのX座標に相当する値を同様に求めることにより、既知点の範囲を拡張していく既知範囲拡張手段を備えるレンズ形状予測装置。 - 請求項7又は8記載のレンズの形状予測装置において、
前記曲面算出手段が求める前記曲面が球面であるレンズの形状予測装置。 - 請求項7又は8記載のレンズの形状予測装置において、
前記曲面算出手段は、
前記各近傍点についての前記基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい球面を求める球面算出手段と、
求めた前記球面の曲率中心から該球面上の頂点へ向かう直線を、中心における法線とし、かつ前記各近傍点からの距離が最も小さい乱視放物面を前記曲面として求める乱視放物面算出手段と
からなるレンズの形状予測装置。 - 請求項7乃至10の何れか記載のレンズの形状予測装置と、
このレンズの形状予測装置によって生成されたレンズの形状データを用いて、当該レンズの性能を表示する性能表示手段と
を含むレンズの性能表示装置。 - 請求項7乃至10の何れか記載のレンズの形状予測装置と、
このレンズの形状予測装置によって生成されたレンズの形状データを用いて、当該レンズを通して外界を見たときの見え方をシミュレーションするシミュレーション手段と
を含む眼光学系のシミュレーション装置。 - 基準面をYZ平面とし、該基準面からの高さをX座標とする基準座標系でのXYZ座標に相当する値が予め分かっている既知点を、レンズ面上に複数有するレンズの形状を予測するためにコンピュータに対して、
前記既知点のうち、前記基準座標系でのX座標に相当する値が分かっていない未知点の近傍に位置するものを、近傍点として複数特定する近傍点特定処理と、
特定した前記各近傍点についての前記基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい曲面を求める曲面算出処理と、
求めた前記曲面を定義する式に、前記未知点の前記基準座標系でのYZ座標に相当する値を代入することにより、該未知点の前記基準座標系でのX座標に相当する値を求める高さ算出処理と
を実行させるプログラム。 - 請求項13記載のプログラムにおいて、
前記コンピュータに対して、
前記高さ算出処理で、前記基準座標系でのX座標に相当する値が求められた未知点を前記既知点とみなし、該既知点の近傍に位置する未知点の前記基準座標系でのX座標に相当する値を同様に求めることにより、既知点の範囲を拡張していく既知範囲拡張処理を実行させるプログラム。 - 請求項13又は14記載のプログラムにおいて、
前記曲面算出処理で求める前記曲面が球面であるプログラム。 - 請求項13又は14記載のプログラムにおいて、
前記曲面算出処理は、
まず、前記各近傍点についての前記基準座標系でのXYZ座標に相当する値を用いて、それら各近傍点からの距離が最も小さい球面を求める球面算出処理と、
次いで、求めた前記球面の曲率中心から該球面上の頂点へ向かう直線を、中心における法線とし、かつ前記各近傍点からの距離が最も小さい乱視放物面を前記曲面として求める乱視放物面算出処理と
からなるプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002265411A JP3979904B2 (ja) | 2002-09-11 | 2002-09-11 | レンズの形状予測方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002265411A JP3979904B2 (ja) | 2002-09-11 | 2002-09-11 | レンズの形状予測方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004101415A JP2004101415A (ja) | 2004-04-02 |
JP3979904B2 true JP3979904B2 (ja) | 2007-09-19 |
Family
ID=32264562
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002265411A Expired - Fee Related JP3979904B2 (ja) | 2002-09-11 | 2002-09-11 | レンズの形状予測方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3979904B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9778136B2 (en) | 2012-05-11 | 2017-10-03 | Menicon Co., Ltd. | Lens power measurement device and measurement method |
JP2021033036A (ja) * | 2019-08-23 | 2021-03-01 | 東海光学株式会社 | 累進屈折力レンズの不明な領域の光学性能を補間するための補間方法 |
-
2002
- 2002-09-11 JP JP2002265411A patent/JP3979904B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2004101415A (ja) | 2004-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20180349531A1 (en) | Topology optimization for subtractive manufacturing techniques | |
Wang et al. | Contact detection between polygonal blocks based on a novel multi-cover system for discontinuous deformation analysis | |
US5801709A (en) | 3-dimensional data shaping apparatus | |
Hel-Or et al. | Relaxed parametric design with probabilistic constraints | |
EP1473678A1 (en) | Three-dimensional shape displaying program, three-dimensional shape displaying method, and three-dimensional shape displaying device | |
JP2010505170A (ja) | 数値的にシミュレートされた2つのオブジェクト間の衝突を検出する方法および装置 | |
EP2336977A1 (en) | A Micro-architecture System and Method for ray tracing and collision detection | |
JP3979904B2 (ja) | レンズの形状予測方法 | |
US20040267391A1 (en) | Method for determining the effects of manufacturing decisions | |
CN107567641B (zh) | 用于识别复制品的系统和方法 | |
US11557078B1 (en) | Machine learning simulation of finite element analysis in augmented reality | |
EP1685509B1 (en) | System, method, and computer program product for determining wall thickness in a graphic model | |
JP3928016B2 (ja) | 最大角度法を用いた三角形メッシュ生成方法及びプログラム | |
JPH08285570A (ja) | 形状パラメータ測定方法および装置 | |
KR100759743B1 (ko) | 투명체의 동적 투시 변형 평가방법 및 투명체의 삼차원형상 설계 지원 방법 | |
JP2004094750A (ja) | 3次元形状処理装置、干渉モデル検出方法、プログラム、及び記憶媒体 | |
KR20030070426A (ko) | 형상 기반의 삼각망 생성 방법 | |
Gregory et al. | Using spline surfaces in optical design software | |
CN113758458B (zh) | 一种复曲面镜的面形测量方法 | |
JPH07260464A (ja) | ワークの組立可能性判定方法および装置 | |
JP2020004311A (ja) | 解析装置、解析方法及び解析プログラム | |
JP2002236940A (ja) | テンプレートを用いた自由曲面の創成方法および装置 | |
Tengler et al. | Ways of uneven road surface modelling used in the vehicle dynamics analysis | |
CN118013807A (zh) | 一种凹坑变形的优化方法、装置、电子设备及存储介质 | |
Harashima et al. | Product Form Evaluation Methods in the Design Process and Their Applications to Reforming |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050221 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070626 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070626 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100706 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110706 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110706 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120706 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120706 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130706 Year of fee payment: 6 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |