CN114945859A - 镜片评价方法、镜片设计方法、眼镜镜片的制造方法以及镜片评价程序 - Google Patents
镜片评价方法、镜片设计方法、眼镜镜片的制造方法以及镜片评价程序 Download PDFInfo
- Publication number
- CN114945859A CN114945859A CN202180008103.0A CN202180008103A CN114945859A CN 114945859 A CN114945859 A CN 114945859A CN 202180008103 A CN202180008103 A CN 202180008103A CN 114945859 A CN114945859 A CN 114945859A
- Authority
- CN
- China
- Prior art keywords
- lens
- evaluation
- eyeglass
- vector
- predetermined
- 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.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims description 49
- 238000013461 design Methods 0.000 title claims description 14
- 238000004519 manufacturing process Methods 0.000 title claims description 11
- 239000013598 vector Substances 0.000 claims abstract description 60
- 239000011159 matrix material Substances 0.000 claims description 46
- 230000003287 optical effect Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 description 21
- 230000004075 alteration Effects 0.000 description 19
- 238000012545 processing Methods 0.000 description 16
- 238000009826 distribution Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000001131 transforming effect Effects 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 201000009310 astigmatism Diseases 0.000 description 1
- 210000005252 bulbus oculi Anatomy 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0292—Testing optical properties of objectives by measuring the optical modulation transfer function
-
- G—PHYSICS
- G02—OPTICS
- G02C—SPECTACLES; SUNGLASSES OR GOGGLES INSOFAR AS THEY HAVE THE SAME FEATURES AS SPECTACLES; CONTACT LENSES
- G02C7/00—Optical parts
- G02C7/02—Lenses; Lens systems ; Methods of designing lenses
- G02C7/024—Methods of designing ophthalmic lenses
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0228—Testing optical properties by measuring refractive power
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0242—Testing optical properties by measuring geometrical properties or aberrations
-
- G—PHYSICS
- G02—OPTICS
- G02C—SPECTACLES; SUNGLASSES OR GOGGLES INSOFAR AS THEY HAVE THE SAME FEATURES AS SPECTACLES; CONTACT LENSES
- G02C7/00—Optical parts
- G02C7/02—Lenses; Lens systems ; Methods of designing lenses
- G02C7/024—Methods of designing ophthalmic lenses
- G02C7/028—Special mathematical design techniques
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Ophthalmology & Optometry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Optics & Photonics (AREA)
- Mathematical Physics (AREA)
- Geometry (AREA)
- Eyeglasses (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量,将相乘得到的量作为所述眼镜镜片的评价指标,并基于该评价指标来进行针对所述眼镜镜片的评价。
Description
技术领域
本发明涉及镜片评价方法、镜片设计方法、眼镜镜片的制造方法以及镜片评价程序。
背景技术
为了提供良好地适应眼镜佩戴者的要求的眼镜镜片,提出了如下技术:在眼镜镜片的设计或制造时,确定与眼镜镜片相关的波面的形状,并且确定眼镜佩戴者的个别眼球模型,基于这些谋求眼镜镜片的面形状的最佳化(例如,参照专利文献1)。
现有技术文献
专利文献
专利文献1:日本专利第6209722号公报
发明内容
发明要解决的课题
眼镜镜片所形成的波面的评价例如能够通过对足够数量的光线的光路长度进行数值计算,并与其配合地继续用泽尔尼克(Zernike)多项式来拟合波面来进行。然而,若为此而导致需要庞大的计算负荷,则需要大量的处理时间、高功能的硬件资源等,在迅速且容易地进行针对眼镜镜片的评价方面,并不优选。
本发明的目的在于,提供一种能够进行迅速且容易的处理,并且谋求针对眼镜镜片的评价的适当化的技术。
用于解决课题的手段
本发明是为了达成上述目的而提出的。
本发明的第1方式是一种镜片评价方法,
对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量,将相乘得到的量作为所述眼镜镜片的评价指标。
本发明的第2方式是根据第1方式所述的镜片评价方法,其中,
所述给定的矢量量是针对所述眼镜镜片的局部的评价范围,对用于得到给定的多项式中的给定的项的展开系数的伪逆矩阵进行傅里叶变换而得的矢量量。
本发明的第3方式是根据第1方式所述的镜片评价方法,其中,
所述给定的矢量量是针对所述眼镜镜片的局部的评价范围,对用于得到给定的正交多项式中的给定的项的展开系数的伪逆矩阵进行傅里叶变换而得的矢量量。
本发明的第4方式是根据第2或第3方式所述的镜片评价方法,其中,
所述给定的项是泽尔尼克多项式中的给定的项。
本发明的第5方式是根据第2或第3方式所述的镜片评价方法,其中,
所述给定的项是泽尔尼克多项式中的给定的多个项的加权和。
本发明的第6方式是根据第2至第5方式中任一方式所述的镜片评价方法,其中,
针对所述展开系数,将进行了基于泽尔尼克多项式的展开时的旋转对称成分的各加权设为1。
本发明的第7方式是一种镜片设计方法,具备:
求出对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量而得到的量作为所述眼镜镜片的评价指标的步骤;以及
使用所述评价指标来进行眼镜镜片的设计的步骤。
本发明的第8方式是一种眼镜镜片的制造方法,具备:
求出对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量而得到的量作为所述眼镜镜片的评价指标的步骤;以及
使用所述评价指标来判定所述眼镜镜片的光学特性的适当与否的步骤。
本发明的第9方式是一种镜片评价程序,使计算机执行如下的步骤:
对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量,将相乘得到的量作为所述眼镜镜片的评价指标的步骤。
发明效果
根据本发明,能够进行迅速且容易的处理,并且谋求针对眼镜镜片的评价的适当化。
附图说明
图1是示出眼镜佩戴者所感受到的朗多环(Landolt ring)的观察方式的一个具体例的说明图,(a)是示出PSF的偏差成为最小的观察方式的例子的图,(b)是示出波面像差成为最小的观察方式的例子的图,(c)是示出空间频率30CPD的对比度成为最高的观察方式的例子的图。
图2是示出基于泽尔尼克多项式的形状分解的一个具体例的说明图。
图3是示出本发明的一个实施方式所涉及的镜片评价方法的步骤的一例的流程图。
具体实施方式
以下,基于附图对本发明的实施方式进行说明。
(1)发明人的见解
首先,在本实施方式的具体的说明之前,对发明人的见解进行说明。
在眼镜佩戴者所利用的眼镜镜片的光学特性中,有时会反映该眼镜佩戴者所感受到的外界的观察方式的喜好。例如,根据眼镜佩戴者的不同而有喜欢重视明亮感(PSF以及低频的对比度的良好程度)的观察方式的人、喜欢重视清晰感(高频的对比度的良好程度)的观察方式的人等,期望能够提供符合各自的喜好的光学特性的眼镜镜片。
具体地,根据眼镜佩戴者的不同,例如存在喜欢如图1的(a)所示那样PSF(Pointspread function:点扩散函数)的偏差成为最小的聚焦位置处的观察方式的人、喜欢如图1的(b)所示那样波面像差成为最小的聚焦位置处的观察方式的人、如图1的(c)所示那样空间频率30CPD的对比度成为最高的聚焦位置处的观察方式的人等。关于在透过眼镜镜片进行观察时成为图1的(a)~(c)的哪个观察方式,只要推测基于该眼镜镜片的聚焦位置(成像位置)即可。
已知的是,基于眼镜镜片的聚焦位置(成像位置)不仅依赖于度数,而且还依赖于像差(特别是球面像差)的量和物体的空间频率。因此,为了正确地推测眼镜镜片的聚焦位置,需要考虑像差量、物体的空间频率等。
由此,关于眼镜镜片的观察方式(即,眼镜镜片的光学特性)的评价,为了得到不仅考虑了度数,而且还考虑了像差量、物体的空间频率等的结果,如已叙述的那样,考虑有效利用眼镜镜片的波面以及泽尔尼克多项式。
所谓眼镜镜片的波面是指透过眼镜镜片并由瞳孔规定直径的光束的波面。
泽尔尼克多项式是在半径为1的单位圆的内部定义的函数(正交多项式)。具体地,由以下的(1)式表示。
[数学式1]
在(1)式中,W(x,y)是坐标x,y中的波面,Zj(x,y)是第j个泽尔尼克多项式,cj是与第j个泽尔尼克多项式对应的泽尔尼克系数,J是用于展开的泽尔尼克多项式的数量。
根据这样的泽尔尼克多项式,全部的面形状能够通过泽尔尼克多项式的相加来(近似地)表达。
具体地,例如,对于某个面形状,如图2所示,能够通过泽尔尼克多项式来进行从0阶像差到n(n为自然数)阶像差的形状分解。
另外,在图中,分别地,由中央附近的框来包围的各成分表示旋转对称成分,除此以外的各成分表示非旋转对称成分。此外,在图中,属于旋转对称成分的2阶像差的成分通常被称为度数误差(散焦),该像差的系数对应于波面像差成为最小的聚焦位置。属于旋转对称成分的4阶像差的成分是相当于球面像差的成分。属于旋转对称成分的各成分的系数的总和对应于PSF成为最小的聚焦位置。
如果有效利用如以上那样的眼镜镜片的波面以及泽尔尼克多项式,则能够谋求眼镜镜片的观察方式的评价的适当化。然而,若为此而导致需要庞大的计算负荷,则需要大量的处理时间、高功能的硬件资源等,在迅速且容易地进行针对眼镜镜片的评价方面,并不优选。特别是,在眼镜镜片为渐进多焦点镜片的情况下,需要每次遍及有效直径内整个区域而进行针对各计测点的计算,因此负荷增大变得非常显著。
鉴于以上的情况,本申请的发明人反复进行了潜心研究,结果想到了如下的技术:能够进行基于高速运算的迅速且容易的处理,并且谋求针对眼镜镜片的评价(特别是,对眼镜镜片的观察方式造成较大影响的散焦、球面像差、PSF等评价)的适当化。以下,在本实施方式中,对这样的技术进行具体的说明。
(2)高速运算的概要
如上所述,泽尔尼克多项式由(1)式表示。为了基于(1)式而针对成为评价对象的眼镜镜片的各坐标点求出波面的像差,必须知道泽尔尼克系数cj(的大小)。
为了求出泽尔尼克系数cj,考虑使用如下的方法:将泽尔尼克系数cj作为未知,对各泽尔尼克多项式Zj(x,y)(即函数f)乘以未知的系数来进行加法运算,并进行与测量形状的拟合运算(最小二乘法等)来算出。在该公式化时,利用伪逆矩阵。所谓伪逆矩阵是对不具有逆矩阵的矩阵实施数学操作而设为如逆矩阵那样的矩阵。也被称为广义逆矩阵,例如相当于将逆矩阵扩展为非正方、非正则矩阵的矩阵。
在此,为了进一步详细地说明伪逆矩阵,考虑局部的波面与泽尔尼克系数的关系。例如,在局部的波面W、泽尔尼克多项式Z、以及泽尔尼克系数c处于由W=cZ表示的关系的情况下,若对此实施数学处理,则得到c=(ZTZ)-1ZTW的关系,由此求出一组泽尔尼克系数c。右上角的“T”表示矩阵的转置。此外,“-1”表示逆矩阵。在这样的c=(ZTZ)-1ZTW的关系中,“(ZTZ)-1ZT”的部分相当于伪逆矩阵。以下,也有时将这样的伪逆矩阵简记为“M”。
更具体地,泽尔尼克系数与局部的波面的关系、即泽尔尼克展开由下述的(2)式给出。
[数学式2]
C(x,y)=(ZTZ)-1ZTW(x,y)…(2)
在(2)式中,C是将各泽尔尼克系数排列与成分数对应的量的矢量量,W(x,y)是从表示波面数据整体的W中提取以x,y为中心的局部的数据点数的量作为矢量的波面,Z是表示各点、各成分的泽尔尼克多项式的(局部的数据点数×成分数)的矩阵。
考虑对眼镜镜片的所有地点进行基于上述的(2)式的计算。
在该情况下,如果针对伪逆矩阵而设为(ZTZ)-1ZT=M,则能够通过预先计算(局部的数据点数×成分数)的矩阵M来重复使用,从而将(2)式置换为下述的(3)式。
[数学式3]
C(x,y)=MW(x,y)…(3)根据(3)式,通过重复使用矩阵M,能够减少每一地点的计算量,但导致需要与(局部的数据点数×成分数)对应的量的计算量。因此,若对眼镜镜片的全部地点进行,则导致需要(全部地点的数据点数×局部的数据点数×成分数)的计算量。在此,由于全部地点的数据点数是二维数据,所以是一方向数据数k的平方。此外,局部的数据点数采取如下这样的赋予方法、即局部区域为整个区域的〇〇%的面积,因此与全部地点的数据点数成比例。若考虑这些情况,则每一成分的计算量与一方向数据数k的四次方成比例。
针对这一点,在本实施方式中,通过如以下那样的方法,谋求计算量的削减。
眼镜镜片的所有的地点处的泽尔尼克展开能够视为如以下的(4)式那样的卷积。
[数学式4]
C=M*W…(4)
在傅里叶空间中,卷积能够表示为函数彼此的积。即使在本实施方式中,也对此进行利用。也就是说,如以下的(5)式所示,对进行了傅里叶变换的M与进行了傅里叶变换的W的积进行傅里叶逆变换,由此得到C。
[数学式5]
在(5)式中,F-1[]是[]内的函数的傅里叶逆变换。F[W]是波面W的傅里叶变换。F[M]是矩阵M的傅里叶变换,相当于预先决定的给定的矢量量的一例。在此所说的给定的矢量量表示对用于得到多项式的给定的项的展开系数的伪逆矩阵进行了傅里叶变换的矢量的量,若在数学上表达,则如上所述,例如会相当于对矩阵M进行了傅里叶变换的F[M]。根据O记法(量级(order)记法),已知傅里叶变换为O(NlogN)的计算量。因此,如(5)式那样,使用了傅里叶变换的高速泽尔尼克展开成为与n2logn成比例的计算量,例如,与通过最小二乘法等来进行泽尔尼克展开的情况相比,能够实现大幅的计算量的削减。
也就是说,本实施方式中,在针对眼镜镜片的评价时,求出对该眼镜镜片的波面W的傅里叶变换F[W]乘以预先决定的给定的矢量量即矩阵M的傅里叶变换F[M]而得的量,并基于该量进行该评价。因此,能够实现大幅的计算量的削减,能够抑制用于评价的计算负荷增大,因此在迅速且容易地进行该评价方面,非常优选。(3)眼镜镜片的评价步骤
接着,对使用上述的方法的高速运算来进行针对眼镜镜片的评价的步骤、即本实施方式所涉及的镜片评价方法的步骤的一例进行说明。
设为以下说明的镜片评价方法的步骤使用计算机装置来进行。即,使用具备具有CPU(Central Processing Unit:中央处理单元)等的运算部、闪存、HDD(Hard Disk Drive:硬盘驱动器)等存储器、输入输出接口等这样的硬件资源而构成的计算机装置,并且使运算部执行预先安装于存储器的给定程序,由此进行基于以下说明的步骤的镜片评价。
图3是示出本实施方式所涉及的镜片评价方法的步骤的一例的流程图。
如图例那样,本实施方式所涉及的镜片评价方法的步骤大致分为预处理(步骤10,以下将步骤简称为“S”。)和反复处理(S20)。
在预处理(S10)中,首先,针对预定佩戴成为评价对象的眼镜镜片的眼镜佩戴者,掌握该眼镜佩戴者的观察方式的喜好(S11)。具体地,掌握眼镜佩戴者例如是喜欢PSF的偏差成为最小的聚焦位置处的观察方式的人、还是喜欢波面像差成为最小的聚焦位置处的观察方式的人、还是喜欢空间频率30CPD的对比度成为最高的聚焦位置处的观察方式的人等。眼镜佩戴者的观察方式的喜好的掌握例如可以通过如下来进行:在与计算机装置的输入输出接口连接的显示画面上输出如图1所示的各图像,并使眼镜佩戴者进行视觉辨认,在这基础上从输入输出接口使其信息输入该眼镜佩戴者喜欢什么样的观察方式。
在掌握了眼镜佩戴者的观察方式的喜好后,接下来,针对眼镜佩戴者预定佩戴的眼镜镜片,决定对该眼镜镜片的度数加入何种程度的球面度数(S12)。具体地,基于眼镜佩戴者的观察方式的喜好的掌握结果,决定与眼镜镜片的波面的面形状相关的泽尔尼克多项式的加权和评价中的各项的泽尔尼克系数cj的加权量。通过决定泽尔尼克系数cj的加权量,从而由该泽尔尼克多项式确定的面形状成为反映了眼镜佩戴者的观察方式的喜好的面形状。
然后,在决定了泽尔尼克系数cj的加权量后,将该决定结果保存在存储器(S13)。具体地,预先计算决定泽尔尼克系数cj的加权量而得到的矩阵M,保存在存储器内的给定存储区域。由此,通过从存储器读出,从而能够重复使用(局部的数据点数×成分数)的矩阵M。
另外,保存在存储器的矩阵M也可以存在多种不同的加权量的矩阵。
在经过如以上那样的预处理(S10)后,会进行反复处理(S20)。在反复处理(S20)时,首先,确定成为评价对象的眼镜镜片的波面(S21)。关于波面的确定,其方法没有特别限定,而例如能够通过利用了波动光学的计算的模拟处理来进行。
在确定了成为评价对象的眼镜镜片的波面后,接下来,对所确定的眼镜镜片的波面进行泽尔尼克展开(S22)。
具体地,首先,从存储器读出用于求解局部的泽尔尼克展开的矩阵M。矩阵M相当于将泽尔尼克多项式的项数汇总为一个的伪逆矩阵。因此,通过在此读出矩阵M来重复使用,从而能够省去考虑泽尔尼克多项式的各项的麻烦。
在读出矩阵M后,对所确定的眼镜镜片的波面进行傅里叶变换,并且对所读出的矩阵M进行傅里叶变换。然后,如(5)式所示,在傅里叶空间中,将波面的傅里叶变换与矩阵M的傅里叶变换相乘,在这基础上,对它们的积进行傅里叶逆变换。由此,得到和波面W与矩阵M的卷积(参照(4)式)同等的结果、即成为评价对象的眼镜镜片的所有的地点处的泽尔尼克展开的结果。此时,不是实际进行卷积,而是进行利用了傅里叶空间的计算处理,因此能够通过一次计算处理来完成与卷积对应的量的循环(loop),与进行卷积的情况相比,能够实现大幅的计算量的削减。
然后,基于泽尔尼克展开的结果,针对成为评价对象的眼镜镜片,评价考虑了眼镜佩戴者的喜好的度数分布(S23)。此时的评价指标使用眼镜镜片的波面的泽尔尼克展开的结果。即,将对眼镜镜片的波面的傅里叶变换(即,对眼镜镜片的波面进行傅里叶变换而得到的矢量量)乘以作为预先决定的给定的矢量量的一例的矩阵M的傅里叶变换而得的量、更详细地说对它们的积进行傅里叶逆变换而得到的量用作评价指标。
具体地,将进行泽尔尼克展开而得到的结果用作评价指标,针对成为评价对象的眼镜镜片,判定作为光学特性之一的度数分布的适当与否。适当与否的判定只要根据度数分布是否与预先决定的判定基准一致来进行即可。例如,如果与判定基准一致,则判定为是适当的度数分布,如果与判定基准不一致,则判定为是不适当的度数分布。在此所判定的度数分布将进行泽尔尼克展开而得到的结果用作评价指标,因此基于眼镜镜片的聚焦位置(成像位置)不仅考虑了度数,而且还考虑了像差量、物体的空间频率等。因此,该判定结果反映了眼镜佩戴者所感受到的观察方式的喜好,其结果是,能够谋求针对眼镜镜片的评价的适当化。
另外,针对面形状不同的眼镜镜片,在重新进行镜片评价的情况下,只要再次进行上述的反复处理(S20)即可。
经过如以上那样的步骤,在本实施方式中,进行针对眼镜镜片的镜片评价。也就是说,本实施方式所涉及的镜片评价方法可以说是“对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量,将相乘得到的量作为所述眼镜镜片的评价指标的镜片评价方法。”。
此外,如已叙述的那样,本实施方式所涉及的镜片评价方法的步骤使用计算机装置来进行。在该情况下,用于使计算机装置执行镜片评价方法的步骤的给定程序只要能够安装在该计算机装置,就可以存放在由该计算机装置能够读取的记录介质(例如,磁盘、光盘、光磁盘、半导体存储器等)而提供,也可以通过因特网、专用线路等网络而从外部提供。这样,安装在计算机装置的给定程序可以说是“使计算机执行如下的步骤的镜片评价程序,所述步骤为对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量,将相乘得到的量作为所述眼镜镜片的评价指标的步骤。”。
(4)眼镜镜片的设计方法等
上述的步骤的镜片评价方法能够应用于眼镜镜片的设计方法。也就是说,也可以设为使用通过上述的步骤求出的评价指标来进行眼镜镜片的设计。如果这样,就能够设计使其反映了眼镜佩戴者的观察方式的喜好的眼镜镜片。也就是说,本实施方式中,也可以设为通过“具备求出对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量而得到的量作为所述眼镜镜片的评价指标的步骤、和使用所述评价指标来进行眼镜镜片的设计的步骤的镜片设计方法。”来进行眼镜镜片的设计。
根据这样的镜片设计方法,能够实现如以下那样的处理。例如,在泽尔尼克多项式中,将针对球面像差的项的泽尔尼克系数设为CSA,将矩阵设为MSA,将针对散焦的项的泽尔尼克系数设为Cdefocus,将矩阵设为Mdefocus(参照图2)。在该情况下,原本,关于每一个,成立W×MSA=CSA以及W×Mdefocus=Cdefocus这两个式子。相对于此,根据上述的镜片设计方法,由于进行利用了傅里叶空间的计算处理,所以能够汇总为W×(MSA+Mdefocus)=CSA+Cdefocus这样的一个式子。该情况意味着,即使在适当调整泽尔尼克系数的加权量的情况下,也能够用一个式子来进行应对。例如,如果是对于泽尔尼克系数的加权量,CSA设为0.5,Cdefocus设为1的情况,则能够用W×(0.5×MSA+1×Mdefocus)=0.5×CSA+1×Cdefocus这样的一个式子来进行应对。因此,针对在镜片设计时所需的计算处理,能够确保高的通用性、自由度等,对于镜片设计者而言便利性非常优异。
此外,上述的步骤的镜片评价方法能够应用于眼镜镜片的制造方法。也就是说,也可以设为,在眼镜镜片的制造时,使用通过上述的步骤来求出的评价指标,来判定眼镜镜片的光学特性的适当与否。如果这样,就能够制造使其反映了眼镜佩戴者的观察方式的喜好的眼镜镜片。也就是说,在本实施方式中,也可以设为通过“具备求出对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量而得到的量作为所述眼镜镜片的评价指标的步骤、和使用所述评价指标来判定所述眼镜镜片的光学特性的适当与否的步骤的眼镜镜片的制造方法。”来进行眼镜镜片的制造。
(5)本实施方式所带来的效果
根据本实施方式,得到以下所示的一个或多个效果。
本实施方式中,对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量,将相乘得到的量作为该眼镜镜片的评价指标。然后,作为该评价指标,使用进行泽尔尼克展开而得到的结果。因此,在使用了该评价指标的评价中,基于眼镜镜片的聚焦位置(成像位置)不仅考虑了度数,而且还考虑了像差量、物体的空间频率等。因此,本实施方式中的眼镜镜片的镜片评价反映了眼镜佩戴者所感受到的观察方式的喜好,其结果是,能够谋求针对眼镜镜片的评价的适当化。
而且,在本实施方式中,在进行泽尔尼克展开时,求出对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量而得到的量。然后,作为给定的矢量量,使用对伪逆矩阵进行傅里叶变换而得的矢量量。也就是说,不是实际进行卷积,而是进行利用了傅里叶空间的计算处理。因此,能够实现大幅的计算量的削减,能够抑制用于评价的计算负荷增大。如果不需要庞大的计算负荷,则也不需要大量的处理时间、高功能的硬件资源等,因此在迅速且容易地进行眼镜镜片的评价方面,非常优选。
如以上那样,根据本实施方式,能够进行迅速且容易的处理,并且谋求针对眼镜镜片的评价的适当化。
(6)变形例等
以上对本发明的实施方式进行了说明,但上述的公开内容示出了本发明的例示性的实施方式。即,本发明的技术的范围并不限定于上述的例示性的实施方式,而是在不脱离其主旨的范围内能够进行各种变更。
在上述的实施方式中,作为与对眼镜镜片的波面进行傅里叶变换而得到的矢量量相乘的给定的矢量量,例示了矩阵M的傅里叶变换,但该矢量量也可以如以下那样。例如,给定的矢量量也可以是针对眼镜镜片的局部的评价范围,对用于得到给定的多项式中的给定的项的展开系数的伪逆矩阵进行傅里叶变换而得的矢量量。此外,例如,矢量量也可以是针对眼镜镜片的局部的评价范围,对用于得到给定的正交函数中的给定的项的展开系数的伪逆矩阵进行傅里叶变换而得的矢量量。例如,也可以使用在考虑了斯泰尔斯-克劳福德(Stiles-Crawford)效应的加权圆形区域中正交的多项式等。
此外,例如,在上述的多项式或正交多项式中,给定的项既可以是泽尔尼克多项式中的给定的项,也可以是泽尔尼克多项式中的给定的多个项的加权和。
此外,例如,也可以设为,针对上述的多项式或正交多项式中的展开系数,将进行了基于泽尔尼克多项式的展开时的旋转对称成分(参照图2)的各加权设为1。如果将旋转对称成分的各加权设为1,则能够算出成为PSF最小的散焦位置。但是,也可以省略不具有有意义的量的成分的权重。
此外,例如,对于上述的多项式或正交多项式中的展开系数,也可以应用于散光度数的评价。如果将2阶对称成分的各加权设为1,则能够算出由基于PSF最小位置的剖面引起的散焦位置的偏离。但是,也可以省略不具有有意义的量的成分的权重。
Claims (9)
1.一种镜片评价方法,
对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量,将相乘得到的量作为所述眼镜镜片的评价指标。
2.根据权利要求1所述的镜片评价方法,其中,
所述给定的矢量量是针对所述眼镜镜片的局部的评价范围,将伪逆矩阵进行傅里叶变换而得到的矢量量,所述伪逆矩阵用于得到给定的多项式中的给定的项的展开系数。
3.根据权利要求1所述的镜片评价方法,其中,
所述给定的矢量量是针对所述眼镜镜片的局部的评价范围,将伪逆矩阵进行傅里叶变换而得到的矢量量,所述伪逆矩阵用于得到给定的正交多项式中的给定的项的展开系数。
4.根据权利要求2或3所述的镜片评价方法,其中,
所述给定的项是泽尔尼克多项式中的给定的项。
5.根据权利要求2或3所述的镜片评价方法,其中,
所述给定的项是泽尔尼克多项式中的给定的多个项的加权和。
6.根据权利要求2至5中任一项所述的镜片评价方法,其中,
针对所述展开系数,将进行了基于泽尔尼克多项式的展开时的旋转对称成分的各加权设为1。
7.一种镜片设计方法,具备:
求出对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量而得到的量,作为所述眼镜镜片的评价指标的步骤;以及
使用所述评价指标来进行眼镜镜片的设计的步骤。
8.一种眼镜镜片的制造方法,具备:
求出对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量而得到的量,作为所述眼镜镜片的评价指标的步骤;以及
使用所述评价指标来判定所述眼镜镜片的光学特性的适当与否的步骤。
9.一种镜片评价程序,使计算机执行如下的步骤:
对将眼镜镜片的波面进行傅里叶变换而得到的矢量量乘以预先决定的给定的矢量量,将相乘得到的量作为所述眼镜镜片的评价指标的步骤。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020-051102 | 2020-03-23 | ||
JP2020051102A JP7356940B2 (ja) | 2020-03-23 | 2020-03-23 | レンズ評価方法、レンズ設計方法、眼鏡レンズの製造方法およびレンズ評価プログラム |
PCT/JP2021/000681 WO2021192492A1 (ja) | 2020-03-23 | 2021-01-12 | レンズ評価方法、レンズ設計方法、眼鏡レンズの製造方法およびレンズ評価プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114945859A true CN114945859A (zh) | 2022-08-26 |
CN114945859B CN114945859B (zh) | 2024-03-15 |
Family
ID=77848590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202180008103.0A Active CN114945859B (zh) | 2020-03-23 | 2021-01-12 | 镜片评价方法、镜片设计方法、眼镜镜片的制造方法以及镜片评价程序 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20230115346A1 (zh) |
EP (1) | EP4130856A4 (zh) |
JP (1) | JP7356940B2 (zh) |
CN (1) | CN114945859B (zh) |
WO (1) | WO2021192492A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7280321B2 (ja) * | 2020-05-27 | 2023-05-23 | 株式会社ワコム | 位置指示器 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3200037A1 (de) * | 1982-01-04 | 1983-07-14 | Klaus Dr. 8035 Gauting Dietrich | Anordnung zum testen von kontaktlinsen und gleitsichtglaesern |
JPS6110740A (ja) * | 1984-06-26 | 1986-01-18 | Tokyo Optical Co Ltd | 光学系評価装置 |
JP2000186978A (ja) * | 1998-10-12 | 2000-07-04 | Hoya Corp | 眼鏡レンズの評価方法及び評価装置 |
CN1520270A (zh) * | 2001-06-20 | 2004-08-11 | ������������ʽ���� | 眼镜·隐形眼镜的度数确定系统及其方法 |
JP2011197556A (ja) * | 2010-03-23 | 2011-10-06 | Tokai Kogaku Kk | 眼鏡レンズの設計方法及び眼鏡レンズ又はレンズデータの選択方法 |
CN103080818A (zh) * | 2010-07-16 | 2013-05-01 | 卡尔蔡斯光学公司 | 经过波前优化的渐进镜片 |
CN107209011A (zh) * | 2015-02-05 | 2017-09-26 | 国立大学法人神户大学 | 形状评价方法以及形状评价装置 |
CN110121671A (zh) * | 2016-12-15 | 2019-08-13 | 图茨技术股份有限公司 | 数据眼镜、用于数据眼镜的眼镜片和用于在视网膜上生成图像的方法 |
JP6581325B1 (ja) * | 2019-06-12 | 2019-09-25 | 株式会社アサヒビジョン | レンズ光学特性測定装置、レンズ光学特性測定方法、プログラム、及び、記録媒体。 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10224363A1 (de) * | 2002-05-24 | 2003-12-04 | Zeiss Carl Smt Ag | Verfahren zur Bestimmung von Wellenfrontaberrationen |
US20110153248A1 (en) * | 2009-12-23 | 2011-06-23 | Yeming Gu | Ophthalmic quality metric system |
DE102012000390A1 (de) | 2012-01-11 | 2013-07-11 | Rodenstock Gmbh | Brillenglasoptimierung mit individuellem Augenmodell |
JP5870234B1 (ja) * | 2014-07-03 | 2016-02-24 | オリンパス株式会社 | 偏心量計測方法及び偏心量計測装置 |
-
2020
- 2020-03-23 JP JP2020051102A patent/JP7356940B2/ja active Active
-
2021
- 2021-01-12 EP EP21775072.8A patent/EP4130856A4/en active Pending
- 2021-01-12 CN CN202180008103.0A patent/CN114945859B/zh active Active
- 2021-01-12 US US17/914,244 patent/US20230115346A1/en active Pending
- 2021-01-12 WO PCT/JP2021/000681 patent/WO2021192492A1/ja unknown
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3200037A1 (de) * | 1982-01-04 | 1983-07-14 | Klaus Dr. 8035 Gauting Dietrich | Anordnung zum testen von kontaktlinsen und gleitsichtglaesern |
JPS6110740A (ja) * | 1984-06-26 | 1986-01-18 | Tokyo Optical Co Ltd | 光学系評価装置 |
JP2000186978A (ja) * | 1998-10-12 | 2000-07-04 | Hoya Corp | 眼鏡レンズの評価方法及び評価装置 |
CN1520270A (zh) * | 2001-06-20 | 2004-08-11 | ������������ʽ���� | 眼镜·隐形眼镜的度数确定系统及其方法 |
JP2011197556A (ja) * | 2010-03-23 | 2011-10-06 | Tokai Kogaku Kk | 眼鏡レンズの設計方法及び眼鏡レンズ又はレンズデータの選択方法 |
CN103080818A (zh) * | 2010-07-16 | 2013-05-01 | 卡尔蔡斯光学公司 | 经过波前优化的渐进镜片 |
CN107209011A (zh) * | 2015-02-05 | 2017-09-26 | 国立大学法人神户大学 | 形状评价方法以及形状评价装置 |
CN110121671A (zh) * | 2016-12-15 | 2019-08-13 | 图茨技术股份有限公司 | 数据眼镜、用于数据眼镜的眼镜片和用于在视网膜上生成图像的方法 |
JP6581325B1 (ja) * | 2019-06-12 | 2019-09-25 | 株式会社アサヒビジョン | レンズ光学特性測定装置、レンズ光学特性測定方法、プログラム、及び、記録媒体。 |
Also Published As
Publication number | Publication date |
---|---|
EP4130856A4 (en) | 2024-05-15 |
US20230115346A1 (en) | 2023-04-13 |
WO2021192492A1 (ja) | 2021-09-30 |
JP2021149037A (ja) | 2021-09-27 |
CN114945859B (zh) | 2024-03-15 |
JP7356940B2 (ja) | 2023-10-05 |
EP4130856A1 (en) | 2023-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6655803B1 (en) | Wavefront method for designing optical elements | |
US10949954B2 (en) | Imaging system and method for imaging objects with reduced image blur | |
CN103080818B (zh) | 经过波前优化的渐进镜片 | |
JP3342423B2 (ja) | 眼光学系のシミュレーション装置 | |
Griffiths et al. | Analysis of cornea curvature using radial basis functions–Part I: Methodology | |
CN114945859A (zh) | 镜片评价方法、镜片设计方法、眼镜镜片的制造方法以及镜片评价程序 | |
EP2916164B1 (en) | Method for designing spectacle lens, and system for designing spectacle lens | |
AU2020246957A1 (en) | A device and method for evaluating a performance of a visual equipment for a visual task | |
JP2001318349A (ja) | 単焦点眼鏡レンズの設計方法、製造方法、及び製造システム | |
JP5061331B2 (ja) | 累進屈折力レンズの設計方法及び製造方法 | |
KR20210098245A (ko) | 광학소자 설계 방법 및 그 장치 | |
Csoba et al. | Fast rendering of central and peripheral human visual aberrations across the entire visual field with interactive personalization | |
Gwak et al. | Modeling nonstationary lens blur using eigen blur kernels for restoration | |
US20230408844A1 (en) | Ophthalmic lens | |
JP3735842B2 (ja) | 眼光学系のシミュレーション装置を駆動するためのプログラムを記録したコンピュータ読み取り可能な記録媒体 | |
Dillon et al. | Optimization of freeform spectacle lenses based on high-order aberrations | |
Ledesma et al. | Artificial Intelligence to design a mask insensible to the distance from the camera to the scene objects | |
JP7414447B2 (ja) | 眼鏡レンズ設計システム、眼鏡レンズ設計方法、及び眼鏡レンズ設計プログラム | |
Xu et al. | Optical-digital joint optimization enables advanced specifications freeform imaging system design | |
Gross et al. | Principles of optimization | |
Vasiljević et al. | Damped least squares optimization implementation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |