CN102297758A - 折射率分布测量方法和折射率分布测量设备 - Google Patents
折射率分布测量方法和折射率分布测量设备 Download PDFInfo
- Publication number
- CN102297758A CN102297758A CN2011101322348A CN201110132234A CN102297758A CN 102297758 A CN102297758 A CN 102297758A CN 2011101322348 A CN2011101322348 A CN 2011101322348A CN 201110132234 A CN201110132234 A CN 201110132234A CN 102297758 A CN102297758 A CN 102297758A
- Authority
- CN
- China
- Prior art keywords
- checking matter
- refractive index
- index distribution
- medium
- transmitted wave
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 108
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000005259 measurement Methods 0.000 claims abstract description 49
- 230000003287 optical effect Effects 0.000 claims description 35
- 238000010008 shearing Methods 0.000 claims description 21
- 238000000691 measurement method Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 4
- 238000000465 moulding Methods 0.000 claims description 2
- 230000014509 gene expression Effects 0.000 description 46
- 230000004075 alteration Effects 0.000 description 21
- LFEUVBZXUFMACD-UHFFFAOYSA-H lead(2+);trioxido(oxo)-$l^{5}-arsane Chemical compound [Pb+2].[Pb+2].[Pb+2].[O-][As]([O-])([O-])=O.[O-][As]([O-])([O-])=O LFEUVBZXUFMACD-UHFFFAOYSA-H 0.000 description 12
- 230000008569 process Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 5
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000035876 healing Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
- 238000012360 testing method 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
-
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及折射率分布测量方法和折射率分布测量设备以及制造光学元件的方法。该折射率分布测量方法包括如下步骤:测量具有第一折射率的第一介质中的第一透射波前和具有不同于第一折射率的第二折射率的第二介质中的第二透射波前,以及通过利用第一透射波前和第二透射波前的测量结果以及基准被检物的每个透射波前移除被检物的形状分量,获得多个取向中的每一个取向中的被检物的折射率分布投影值,该基准被检物具有与被检物相同的形状以及具有特定的折射率分布,并且以与被检物的取向相同的取向位于第一介质和第二介质中的一个中,以及基于与该多个取向对应的多个折射率分布投影值计算被检物的三维折射率分布。
Description
技术领域
本发明涉及折射率分布测量方法和折射率分布测量设备。
背景技术
日本专利特开No.(″JP″)01-316627提出了一种用于通过在被检物被浸入具有与该被检物的折射率近似相同的折射率的介质(匹配油)中的同时测量被检物的透射波前来求得该被检物的折射率分布的方法。JP 02-008726提出了一种用于通过在被检物被浸入具有与该被检物的折射率略有不同的折射率的两种类型的匹配油中的每一种中的同时测量被检物的透射波前来求得该被检物的折射率分布的方法。
JP 01-316627和02-008726中公开的这些方法需要各自均具有与被检物的折射率近似相同的折射率的匹配油。然而,具有高折射率的匹配油具有低透射率,并且,检测器只能输出弱信号。这样,具有高折射率的被检物的测量精确度容易降低。
发明内容
本发明提供能够高度精确地测量被检物的折射率分布的折射率分布测量方法和折射率分布测量设备。
根据本发明的折射率分布测量方法包括以下步骤:通过把被检物布置在具有与该被检物的折射率不同的折射率的介质中并且通过将参考光引入到被检物中,测量被检物的透射波前;以及通过使用透射波前的测量结果,计算被检物的折射率分布。在被检物在介质中的互不相同的多个取向(orientation)中,该测量步骤测量具有第一折射率的第一介质中的第一透射波前和具有与第一折射率不同的第二折射率的第二介质中的第二透射波前。该计算步骤通过利用第一透射波前和第二透射波前的测量结果以及基准被检物的每个透射波前移除被检物的形状分量,来获得在所述多个取向中的每一个取向中的被检物的折射率分布投影值,所述基准被检物具有与所述被检物相同的形状以及具有特定的折射率分布,并且以与所述被检物的取向相同的取向位于第一介质和第二介质中的一个中。然后,该计算步骤基于与所述多个取向对应的多个折射率分布投影值来计算该被检物的三维折射率分布。
从参照附图对示例性实施例的以下描述,本发明的其它特征将变得清晰。
附图说明
图1是根据第一实施例的折射率分布测量设备的框图。
图2是用于图示根据第一实施例的折射率分布测量方法的流程图。
图3A~3B是用于图示根据第一实施例的基准被检物中设置的坐标系和折射率分布测量设备中的光线的光路的视图。
图4A~4B是用于图示根据第一实施例的被检物的倾斜的视图。
图5是根据第二实施例的折射率分布测量设备的框图。
图6是用于图示根据第二实施例的折射率分布测量方法的流程图。
具体实施方式
现在将参照附图给出对本发明的实施例的描述。
第一实施例
图1是安装在稳定器190上的折射率分布测量设备的框图。该折射率分布测量设备通过把被检物浸入两种类型的介质(诸如空气和水)中的每一种中并且通过把来自光源的参考光引入到被检物中,测量被检物的透射波前,其中每种所述介质具有与被检物的折射率不同的折射率。然后,折射率分布测量设备利用作为计算机的处理器和透射波前的测量结果,计算被检物的折射率分布。这个实施例利用Talbot干涉仪作为被配置为通过利用来自光源的光测量布置在介质中的被检物的透射波前的测量单元。
被检物140是诸如透镜之类的光学元件。容器130容纳介质1(诸如空气)或者介质2(诸如水)。空气或水的折射率比被检物140的折射率小0.01或更多。
从诸如He-Ne激光器之类的激光光源100沿光学轴发射的激光束101在穿过针孔板(光学部件)110中的针孔(PH)112时被衍射。在针孔112中被衍射的衍射光或者参考光被准直透镜(CL)120转换为会聚光103。
会聚光103透射容器130中的介质1或2以及被检物140。本实施例假设被检物140是围绕轴旋转对称的透镜。针孔112的直径Φ小得以使衍射光102能够被看作是理想球面波,并且直径Φ被设计为使用在被检物侧的数值孔径NAO和激光光源100的波长λ而满足以下表达式:
表达式1
当λ是600nm并且NAO是大约0.3时,针孔112的直径Φ可以是大约2μm。
已经透射了容器130中的被检物140和空气或水的激光束穿过作为二维衍射光栅的正交衍射光栅170,并且被图像拾取元件(CCD传感器或CMOS传感器)180捕获(测量)。该正交衍射光栅170和图像拾取元件180在此后有时将被称作“传感器”。
当像侧的被检物140的数值孔径(NA)是小的并且衍射光栅170和图像拾取元件180之间的(Talbot)距离z满足等式2表示的Talbot条件时,在图像拾取元件180上作为干涉图案获得衍射光栅170的伪分辨率,其中m是非0整数,d是衍射光栅170的光栅节距,Z0是从衍射光栅170到被检物140的像面的距离。光栅节距d是根据被检物140的像差的大小而确定的。
表达式2
被检物140被配置为可以通过旋转单元150而围绕与光学轴垂直的轴旋转,并且可以通过平行移动和偏心单元160而在光学轴方向上相对移动。旋转单元150用作调整器,该调整器被配置为调整被检物在介质中的取向。准直透镜120、衍射光栅170和图像拾取元件180也被配置为在平行于光学轴安装的轨道(未示出)上相对移动。
图2是用于示出根据本实施例的折射率分布测量方法的流程图,“S”代表步骤。该折射率分布测量方法作为计算机程序而被图1中示出的诸如微计算机之类的处理器200执行。
如图1中所示,作为具有第一折射率的第一介质的介质1(空气)被填充在容器130中(S10)。
接着,根据步骤A测量被浸入容器130中的介质1中的被检物140的(第一)波前像差W1(S20)。
透射波前的测量结果包含被检物的折射率分布、被检物形状的影响、被检物形状误差的影响、以及测量***造成的偏移。在这些项目当中,被检物形状的影响和测量***造成的偏移通过仿真而被计算,并且被从透射波前的测量结果中移除。步骤A求得波前像差W1,并且获得包含被检物的折射率分布和被检物形状误差的影响的剩余信息。
步骤A起初确定每个组件的光学布置,即,针孔板110、准直透镜120、容器130、衍射光栅170和图像拾取元件180之间的在光学轴方向上的间隔(S201)。该光学布置要将NA限制为等于或低于大约0.3并且使图像拾取元件180上的光束尺寸合适,以便在Talbot干涉仪中在图像拾取元件的整个表面上获得衍射光栅170的伪分辨率。该光学布置防止穿过了被检物140上的不同位置的光束会聚在图像拾取元件180上的相同点上,从而在之后的步骤中使图像拾取元件180上的位置与被检物140上的位置相关。
接着,根据所确定的光学布置来布置每个组件,并且,被检物140与传感器对准(S202)。通过使用图1中的平行移动和偏心单元160的相对移动和/或轨道(未示出)上的相对移动来执行该对准。图1中示出的被检物140是凹透镜,但是如果被检物140是凸透镜,则容器130可以被安装在准直透镜120的会聚位置后方(在衍射光栅170侧)从而使图像拾取元件180上的光束尺寸合适。
接着,通过假设不具有折射率分布的理想折射率分布(即,特定折射率分布)来计算透射波前的仿真波前Wsim(S203)。本实施例把具有已知形状(该形状与本实施例中的被检物的形状相同)和已知特定折射率分布的被检物称作基准被检物,并且将其透射波前称作参考透射波前。在S203中,通过以与所述被检物的取向相同的取向把基准被检物布置在第一介质和第二介质中的每一个中来获得每个透射波前。
该已知折射率分布可以是设计值或者测量值。基于基准被检物的坐标(x,y)处的表达式3的关系,求得该仿真波前Wsim:
表达式3
Wsim(x,y)=OPsim(x,y)-OPsim(0,0)
OPsim(x,y)=L1(x,y)+L2(x,y)N1+L3(x,y)Ng+L4(x,y)N1+L5(x,y)
这里,L1到L5是沿着如图3B所示的光线103的组件之间的几何距离。该光线103示意性地示出了穿过图3A中所示的基准被检物141中的点(x,y)的光线。N1表示空气的折射率,Ng表示基准被检物141的理想折射率。通过将被检物140的折射率分布替换为已知值来考虑基准被检物141。为了简化表达式,容器130的壁的厚度被忽略。
接着,在被检物140被浸入空气中的同时,测量第一介质中的被检物140的(第一)透射波前Wm(测量步骤)(S204)。S204包含通过图像拾取元件180获取干涉图案的图像和通过处理器200对透射波前进行图像恢复。透射波前的图像恢复(此后被称作“波前恢复”)是通过快速傅立叶变换(“FFT”)方法执行的。
通过FFT方法进行的波前恢复是通过利用像差干扰干涉图案的载波图案这一特性而将像差与载波图案分离的方法。更特别地,对干涉图案执行二维FFT,并且,该干涉图案被转换为频率图(frequencymap)。接着,在仅在频率图中拾取载波频率附近的部分之后,执行逆快速傅立叶变换(iFFT)方法,并且,执行坐标变换从而把载波频率设置为原点。从而,求得复振幅的相位项。所得到的相位图成为透射波前。
用L1到L5如下表达Wm:
表达式4
Wm(x,y)=OPm(x,y)-OPm(0,0)
这里,N横杠是在坐标(x,y)处在被检物140的光路方向上取平均值的折射率分布投影值,dL表示在坐标(x,y)处的被检物140的厚度误差。
接着,通过利用以下的表达式求得与仿真波前Wsim和透射波前Wm之间的差对应的(第一)波前像差W1(S205)。为了简化表达式,假设折射率Ng等于被检物140在光学轴上的折射率N(0,0):
表达式5
S203独立于S202或S204,并且可以在S201到S205之间的任何定时被执行。
接着,在作为具有第二折射率的第二介质的介质2(水)被填充在容器130中的同时,被检物140被安装在容器130中(S30)。接着,根据上面的步骤A,被检物140的(第二)波前像差W2被测量(测量步骤)(S40),其中N2表示水的折射率。此时,在被检物140被浸在水中的同时,S204中的测量步骤测量第二介质中的被检物140的(第二)透射波前Wm(S204)。
表达式6
接着,通过从波前像差W1和波前像差W2中移除被检物140的形状分量dL,计算被检物140的折射率分布投影值(S50)。S50是用于通过使用布置在相同位置处的被检物的两个波前像差W1和W2来移除被检物的形状误差的影响而获得包含被检物140的折射率分布的信息的折射率分布投影值的计算步骤。这里,使用表达式8的近似。
表达式7
表达式8
从而,对于相对于被检物140的第一被检物倾斜,求得折射率分布投影值,该第一被检物倾斜是其中光学轴与被检物140的旋转对称轴一致的情况下的倾斜。由于折射率分布投影值是在入射到被检物140上的光的光路方向上取平均值的折射率,因此必须通过以不同于第一被检物倾斜的倾斜将光引入到被检物140中求得折射率分布投影值,以便获得三维折射率分布信息。现在将给出对这个方法的描述。
为了以被检物140的多个取向求得被检物140的折射率分布投影值,该被检物140被旋转并且被偏心(S61)。S60的测量次数根据要求得的折射率分布而不同。当能够假设被检物140具有围绕轴旋转对称的形状并且折射率分布围绕同一轴而旋转对称时,测量次数可以是两次。例如,其中光学轴与被检物140的旋转对称轴一致的情况下的取向(第一被检物倾斜)和其中光学轴与被检物140的旋转对称轴不一致的情况下的取向(第二被检物倾斜)。
图4A示出了本实施例中的第一被检物倾斜,图4B示出了本实施例中的第二被检物倾斜。在第一被检物倾斜中,对于被检物140的入射光线103和出射光线104如图4A中所示,并且,在第二被检物倾斜中,入射光线103和出射光线105如图4B中所示。
为了通过小的测量次数精确地获得三维折射率分布信息,可以使折射率分布投影值的多个测量取向互相很不相同。换言之,可以使第一被检物倾斜和第二被检物倾斜互相很不相同。该第二被检物倾斜可以被调整为使得入射光线103能够穿过被检物140的第一表面和第二表面的端部,如图4B中所示。第一表面的端部是作为被检物的光入射侧的光学表面的第一表面的R表面与切割表面之间的边界,并且,第二表面的端部是作为被检物的光出射侧的光学表面的第二表面的R表面与切割表面之间的边界。
本实施例在S61中关于光学轴方向使被检物140平行移动和偏心并且围绕垂直于光学轴的轴旋转被检物140,并且,把被检物140布置为如图4B中所示的位置和角度。
在S61之后,S10到S50被再次执行,直到用于计算折射率分布投影值的测量次数在S60中达到指定次数(在本实施例中是两次)。
当测量次数达到指定次数时,从多个获得的折射率分布投影值计算三维折射率分布(S70)。S70是用于基于与被检物140的多个不同取向对应的多个折射率分布投影值来获得三维折射率分布的信息的计算步骤。通过确定表达三维折射率分布的多项式系数以使得多个计算出的折射率分布投影值能够被再现,而计算出所述三维折射率分布。
当用100条光线表达入射光线103时,折射率分布投影值被表达为如下:
表达式9
N1横杠和N2横杠是对于第一和第二被检物倾斜的折射率分布投影值。另外,假设用12个多项式系数表达要求得的三维折射率分布P。
表达式10
当假设U是使用表达式10的多项式系数作为单位量的折射率分布投影值时,U可用下面的表达式表达:
表达式11
当P被确定为以便满足下面的等式时,P的每个系数能够再现多个所求得的折射率分布投影值。
表达式12
当使用最小二乘法时,如表达式13中那样定义Φ,并且,P的每个系数被确定为使得Φ2能够最小。
表达式13
当使用特征值分解时,通过求得U-1,能够如表达式14中那样直接求得P。
表达式14
作为替代方案,可以使用组合的方法,该方法把Φ定义为从表达式14中的左侧值减去右侧值得到的值并且确定P的每个系数以使得Φ2能够最小,或者可以使用另外的已知方法来求得P。这样,通过求得三维折射率分布P,本实施例中的折射率分布测量方法结束。
如上面讨论的那样,本实施例利用两种类型的介质和从光源发射的参考光测量被检物的两种类型的波前像差,通过从该波前像差中移除被检物的形状分量而获得折射率分布投影值,并且通过改变被检物相对于光学轴的角度而获取另一个折射率分布投影值。基于多个折射率分布投影值求得表达被检物的三维折射率分布的多项式系数。从而,即使当被检物的折射率高时,也能够通过利用折射率低于被检物的折射率的介质,高度精确地测量被检物的内部折射率分布。
为了简单描述,本实施例恰当地设置了表达入射光线103的光线数量和表达三维折射率分布的多项式。更特别地,入射光被表达为一条光线,并且,折射率分布投影值的测量指定次数被设置为m并且表达折射率分布的多项式具有n项。即使在这种情况中,如表达式15中那样设置折射率分布投影值N横杠、折射率分布P以及P的每个系数各具有单位量的情况下的折射率分布投影值U,通过类似于表达式12的方法,也能够求得P:
表达式15
如本实施例中那样,能够通过将Talbot干涉仪用于测量单元来测量由被检物和介质之间的折射率差引起的大的像差。该Talbot干涉仪是一种类型的横向剪切干涉仪,横向剪切干涉仪被配置为测量透射波前与其剪切透射波前之间的差作为干涉图案。
剪切干涉仪是被配置为求得与透射波前的波前形状的梯度对应的量的测量单元。透射波前的横向移位量被称作剪切量,并且,该剪切量与光的直径的比率被称作剪切率。通过减小剪切率,大的透射波前像差能够被作为不使干涉图案密集的小的像差(剪切波前)而被测量。
一般来说,当在剪切干涉仪中剪切率过小时,剪切波前被嵌入到噪声中并且精度劣化。因而,剪切率可以是光瞳直径的3%~5%那么大。然而,本实施例将剪切率设置为1.5%或者更小,例如,大约0.4%~0.9%,从而以小的剪切波前测量具有大的像差的透射波前。
使用Talbot距离Z和图像拾取元件180上的干涉图案数据的直径D,将剪切率定义为(λZ)/(dD),并且,使用表达式2和衍射光栅170上的光束的直径D0,将剪切率定义为(md)/D。因此,该剪切率与衍射光栅170的光栅节距成比例。由表达式2,衍射光栅170的节距影响Talbot距离Z,因而必须通过考虑测量设备的组件之间的干涉而确定该节距。例如,如果假设当m=1时D0是大约10~20mm,则该光栅节距可以是大约40~180μm。
虽然本实施例将两种类型的介质设置为空气和水,但是介质不受限制,只要两种类型的介质相差0.01或者更高即可。另外,这两种类型的介质可以由通过改变其温度而具有不同折射率的相同材料制成。
虽然本实施例讨论对Talbot干涉仪的使用,但是不同于Talbot干涉仪的横向剪切干涉仪和径向剪切干涉仪以及其它剪切干涉仪也可以被使用。
第二实施例
图5是第二实施例的折射率分布测量设备的框图。本实施例的折射率分布测量设备通过使用两种类型的光源和一种类型的介质M来两次测量透射波前,求出折射率分布。该两种类型的光源是,例如,作为光源100A的He-Ne激光器(具有633nm的第一波长)和作为光源100B的YAG激光器的二次谐波(具有不同于第一波长的、532nm的第二波长)。
介质M具有不同于被检物140的折射率的折射率。例如,介质的折射率小于被检物的折射率并且大于空气的折射率。介质M的一个例子是水,以及具有大约1.5到大约1.8的折射率的低折射率油。
针孔板110使用从光源100A或者100B发射的激光束,产生具有理想球面波的(参考)光。该光与图1类似地穿过被检物140,并且,其透射波前通过作为波前测量传感器的Shack-Hartman传感器500而被测量。该Shack-Hartman传感器500沿着光路按照从光源的次序包括透镜阵列510和图像拾取元件520。
类似于第一实施例,准直透镜120、容器130以及传感器500被布置在平行于光学轴的轨道(未示出)上。通过在轨道上移动这些组件,入射到被检物140上的光能够被转换为发散光、准直光和会聚光中的任何一个。从而,入射到Shack-Hartman传感器500上的光束的NA能够被调整。
与Talbot干涉仪相比,Shack-Hartman传感器要求入射到传感器上的光束的NA被严格控制,但是由于不必把衍射光栅170与CCD 160之间的间隔设置为Talbot距离,因此传感器500的对准变得较容易。
Shack-Hartman传感器500将入射到透镜阵列510上的光会聚到CCD上。当倾斜的透射波前入射到透镜阵列510上时,会聚点的位置移位。由于Shack-Hartman传感器500能够将透射波前的倾斜转换为会聚点的位置移位并且测量该位置移位,因此可以测量具有大的像差的波前。
图6是用于图示本实施例的折射率分布测量方法的流程图,“S”代表步骤。该折射率分布测量方法作为计算机程序而被图5中示出的诸如微计算机之类的处理器200执行。图6的大部分流程与图2中的测量流程相同,因而将只讨论不同点。
起初,来自光源100A的光被引入到针孔板110中(S11),并且,使用第一波长的波前像差W1被测量(S20)。接着,来自具有不同于光源100A的波长的波长的光源100B的光被引入到针孔板110(S31)中,并且,波前像差W2被测量(S40)。通过以下表达式来表达在这些步骤中获得的波前像差:
表达式16
这里,NHeNe横杠和NYAG横杠是分别对于第一光源(He-Ne激光器)和第二光源(YAG二次谐波)的在被检物中的位置(x,y)处的折射率分布投影值。NgHeNe和NgYAG是对于相应各光源的被检物的理想折射率(基准被检物的折射率)。NoilHeNe和NoilYAG是对于相应各光源的介质的折射率。
对于第一光源的折射率和对于第二光源的折射率具有以下近似关系:
表达式17
通过利用表达式16和17,能够求得折射率分布投影值(S51)。
表达式18
此后,接着是S60、S61、和S70,然后测量结束。
在表达式18中,当表达式20的Ψ是大的时,测量值W1和W2的误差能够被减小。
表达式19
例如,当介质是空气时,满足Noil≈0。这样,Ψ≈0成立并且测量变得不可行。另外,例如,当假设被检物的折射率并不随波长改变而显著改变时,NgYAG≈NgHeNe成立,并且因而表达式19能够被表达为表达式20。
表达式20
Ψ=NoilYAG-NoilHeNe
在这种情况下,可以选择在第一光源和第二光源之间具有大的折射率差异的介质。为了增大Ψ,需要通过考虑被检物的折射率来确定介质。
本实施例的测量设备可以是这样的测量设备,该测量设备能够测量与透射波前的波前形状的梯度或者光线的倾斜对应的量,并且甚至当透射波前具有大的像差时也能够作为可测物理量检测该梯度或倾斜。因此,测量设备不局限于Shack-Hartman方法,并且可以使用Hartman方法或者Ronchi测试。
通过根据第一和第二实施例的折射率分布测量设备或方法测量的结果可适用于光学元件的制造方法。用于制造光学元件的方法包括如下步骤:基于设计的光学元件模制光学元件,测量所模制的光学元件的形状,评估形状精确度,以及评估满足形状精确度的光学元件的光学性能。本实施例的折射率分布测量方法可适用于评估光学性能的步骤。当所评估的光学性能不满足所要求的规格时,计算光学元件的光学表面的校正量,并且通过使用该结果重新设计光学元件。当所评估的光学性能满足所要求的规格时,光学元件被量产。
由于本实施例的光学元件的制造方法能够高度精确地测量光学元件的内部折射率分布,因此即使当光学元件是由高折射率玻璃材料制成的时,光学元件也能通过模制而被精确地量产。
虽然已参照示例性实施例描述了本发明,应该理解,本发明不限于公开的示例性实施例。以下的权利要求的范围应被给予最宽的解释,以便包含所有这样的变型以及等同的结构和功能。
Claims (8)
1.一种折射率分布测量方法,包含如下步骤:
通过把被检物布置在具有与该被检物的折射率不同的折射率的介质中并且通过将参考光引入到被检物中,测量被检物的透射波前;以及
通过使用透射波前的测量结果,计算被检物的折射率分布,
其特征在于,在被检物在介质中的互不相同的多个取向中,所述测量的步骤测量具有第一折射率的第一介质中的第一透射波前和具有与第一折射率不同的第二折射率的第二介质中的第二透射波前;以及
所述计算的步骤通过利用第一透射波前和第二透射波前的测量结果以及基准被检物的每个透射波前移除被检物的形状分量,来获得在所述多个取向中的每一个取向中的被检物的折射率分布投影值,所述基准被检物具有与所述被检物相同的形状以及具有特定的折射率分布,并且以与所述被检物的取向相同的取向位于第一介质和第二介质中的一个中,并且,所述计算的步骤基于与所述多个取向对应的多个折射率分布投影值来计算该被检物的三维折射率分布。
2.一种折射率分布测量方法,包含如下步骤:
通过把被检物布置在具有与该被检物的折射率不同的折射率的介质中并且通过将参考光引入到被检物中,测量被检物的透射波前;以及
通过使用透射波前的测量结果,计算被检物的折射率分布,
其特征在于,所述测量的步骤测量对于第一波长的第一透射波前时被检物在介质中的取向与所述测量的步骤测量对于不同于第一波长的第二波长的第二透射波前时被检物在介质中的取向不同;以及
所述计算的步骤通过利用对于第一波长的第一透射波前和对于第二波长的第二透射波前的测量结果以及基准被检物的每个透射波前移除被检物的形状分量,来获得在多个取向中的每一个取向中的被检物的折射率分布投影值,所述基准被检物具有与所述被检物相同的形状以及具有特定的折射率分布,并且以与所述被检物的取向相同的取向位于所述介质中,并且,所述计算的步骤基于与所述多个取向对应的多个折射率分布投影值来计算该被检物的三维折射率分布。
3.根据权利要求1或2的折射率分布测量方法,其中,所述多个取向包括其中入射到被检物上的入射光线穿过被检物的光入射侧的被检物的光学表面与切割表面之间的边界以及被检物的光出射侧的被检物的光学表面与切割表面之间的边界的情况下的取向。
4.一种用于制造光学元件的方法,所述方法包括模制光学元件的步骤,
其特征在于,该方法还包括如下步骤:通过利用根据权利要求1到3中任一项的折射率分布测量方法测量光学元件的折射率分布,评估所模制的光学元件的光学性能。
5.一种折射率分布测量设备,包含:
光源;
调整器,该调整器被配置为调整被检物在介质中的取向,该介质具有不同于被检物的折射率的折射率;
测量单元,该测量单元被配置为通过利用来自光源的光,测量布置在介质中的被检物的透射波前;以及
处理器,该处理器被配置为通过使用所述透射波前计算被检物的折射率分布,
其中,在被检物在介质中的互不相同的多个取向中,所述测量单元测量具有第一折射率的第一介质中的第一透射波前和具有与第一折射率不同的第二折射率的第二介质中的第二透射波前;以及
其中,所述处理器通过利用第一透射波前和第二透射波前的测量结果以及基准被检物的每个透射波前移除被检物的形状分量,来获得在所述多个取向中的每一个取向中的被检物的折射率分布投影值,所述基准被检物具有与所述被检物相同的形状以及具有特定的折射率分布,并且以与所述被检物的取向相同的取向位于第一介质和第二介质中的一个中,并且,所述处理器基于与所述多个取向对应的多个折射率分布投影值来计算该被检物的三维折射率分布。
6.一种折射率分布测量设备,包含:
光源,该光源被配置为发射第一波长的光和第二波长的光;
调整器,该调整器被配置为调整被检物在介质中的取向,该介质具有不同于被检物的折射率的折射率;
测量单元,该测量单元被配置为通过利用来自光源的光测量布置在介质中的被检物的透射波前;以及
处理器,该处理器被配置为基于用第一波长测量的第一透射波前和用第二波长测量的第二透射波前来计算被检物的折射率分布,
其中,所述测量单元测量第一透射波前时被检物在介质中的取向与所述测量单元测量第二透射波前时被检物在介质中的取向不同;以及
其中,所述处理器通过利用第一透射波前和第二透射波前的测量结果以及基准被检物的每个透射波前移除被检物的形状分量,来获得在多个取向中的每一个取向中的被检物的折射率分布投影值,所述基准被检物具有与所述被检物相同的形状以及具有特定的折射率分布,并且以与所述被检物的取向相同的取向位于所述介质中,并且,所述处理器基于与所述多个取向对应的多个折射率分布投影值来计算该被检物的三维折射率分布。
7.根据权利要求5或6的折射率分布测量设备,其中,所述测量单元包括剪切干涉仪。
8.根据权利要求5或6的折射率分布测量设备,其中,所述测量单元包括Shack-Hartman传感器。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010-119636 | 2010-05-25 | ||
JP2010119636A JP5021054B2 (ja) | 2010-05-25 | 2010-05-25 | 屈折率分布計測方法および屈折率分布計測装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102297758A true CN102297758A (zh) | 2011-12-28 |
CN102297758B CN102297758B (zh) | 2014-04-09 |
Family
ID=45021867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110132234.8A Expired - Fee Related CN102297758B (zh) | 2010-05-25 | 2011-05-20 | 折射率分布测量方法和折射率分布测量设备 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8525982B2 (zh) |
JP (1) | JP5021054B2 (zh) |
KR (1) | KR101288876B1 (zh) |
CN (1) | CN102297758B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104614342A (zh) * | 2015-02-03 | 2015-05-13 | 清华大学 | 一种高温气流扰动下空气折射率三维重构测量方法 |
CN104662402A (zh) * | 2012-09-28 | 2015-05-27 | 诺华股份有限公司 | 用于眼透镜的屈光力的自动化线内确定的方法 |
CN105339778A (zh) * | 2013-06-28 | 2016-02-17 | 佳能株式会社 | 折射率测量方法、折射率测量装置及光学元件制造方法 |
CN106226032A (zh) * | 2015-06-02 | 2016-12-14 | 佳能株式会社 | 测量方法、测量装置和光学元件的制造方法 |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5008650B2 (ja) * | 2008-12-25 | 2012-08-22 | キヤノン株式会社 | 屈折率分布計測方法及び屈折率分布計測装置 |
JP5328437B2 (ja) * | 2009-03-25 | 2013-10-30 | キヤノン株式会社 | 透過波面測定方法、屈折率分布測定方法、光学素子の製造方法、及び透過波面測定装置 |
JP4968966B2 (ja) * | 2009-12-07 | 2012-07-04 | キヤノン株式会社 | 屈折率分布の計測方法および計測装置 |
JP4895409B2 (ja) | 2010-05-25 | 2012-03-14 | キヤノン株式会社 | 屈折率分布計測方法および屈折率分布計測装置 |
JP5008763B2 (ja) | 2010-12-03 | 2012-08-22 | キヤノン株式会社 | 屈折率分布計測方法、屈折率分布計測装置および光学素子の製造方法 |
GB201107053D0 (en) | 2011-04-27 | 2011-06-08 | Univ Sheffield | Improvements in providing image data |
GB201201140D0 (en) * | 2012-01-24 | 2012-03-07 | Phase Focus Ltd | Method and apparatus for determining object characteristics |
GB201207800D0 (en) | 2012-05-03 | 2012-06-13 | Phase Focus Ltd | Improvements in providing image data |
JP2014016253A (ja) * | 2012-07-09 | 2014-01-30 | Canon Inc | 屈折率分布計測方法、光学素子の製造方法、および、屈折率分布計測装置 |
CN104704339A (zh) | 2012-09-28 | 2015-06-10 | 诺华股份有限公司 | 一种用于眼透镜的中心厚度的自动化线内确定的方法 |
JP6320093B2 (ja) * | 2014-03-14 | 2018-05-09 | キヤノン株式会社 | 屈折率分布計測方法、屈折率分布計測装置、光学素子の製造方法、プログラム、および、記憶媒体 |
BR112017002489B1 (pt) * | 2014-08-08 | 2023-02-14 | Quantum-Si Incorporated | Instrumento configurado para fazer interface com um chip de teste, aparelho, método de análise de um espécime, método para sequenciar uma molécula-alvo de ácido nucleico e método para sequenciamento de ácido nucleico |
KR101716125B1 (ko) * | 2015-07-20 | 2017-03-15 | 주식회사 토모큐브 | 파면 제어기를 활용한 초고속 고정밀 3차원 굴절률 측정 방법 및 장치 |
JP6570693B1 (ja) * | 2018-04-19 | 2019-09-04 | 陽程科技股▲ふん▼有限公司 | 光学位置合わせ検出方法 |
CN108982418B (zh) * | 2018-07-02 | 2022-04-19 | 北京无线电计量测试研究所 | 基于微腔激光色散干涉法的空气折射率测量装置和方法 |
JPWO2022195765A1 (zh) * | 2021-03-17 | 2022-09-22 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH028726A (ja) * | 1988-06-27 | 1990-01-12 | Asahi Optical Co Ltd | 屈折率測定方法 |
JPH0560654A (ja) * | 1991-09-03 | 1993-03-12 | Sumitomo Electric Ind Ltd | 屈折率分布測定方法及び装置 |
JPH08304229A (ja) * | 1995-05-09 | 1996-11-22 | Ricoh Co Ltd | 光学素子の屈折率分布の測定方法および装置 |
JP2000258292A (ja) * | 1999-03-09 | 2000-09-22 | Ricoh Co Ltd | 屈折率分布の測定装置及び屈折率分布の測定方法 |
JP2003021577A (ja) * | 2001-07-09 | 2003-01-24 | Ricoh Co Ltd | 屈折率分布測定方法および装置 |
US20060159332A1 (en) * | 2005-01-19 | 2006-07-20 | Yasuhiro Sawada | Image processing apparatus and refractive index distribution measuring apparatus |
Family Cites Families (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4265534A (en) | 1977-12-23 | 1981-05-05 | Remijan Paul W | Optical apparatus and method for producing the same |
US4542989A (en) | 1981-03-02 | 1985-09-24 | Randwal Instrument Co., Inc. | Apparatus for position encoding |
US4541697A (en) | 1981-03-03 | 1985-09-17 | Randwal Instrument Co., Inc. | Ophthalmic testing devices |
JPS5845526A (ja) | 1981-09-12 | 1983-03-16 | Canon Inc | タルボ干渉計 |
IT1155285B (it) * | 1982-02-10 | 1987-01-28 | Cselt Centro Studi Lab Telecom | Procedimento e apparecchiatura per la determinazione del profilo d'indice di rifrazione di fibre ottiche e preforme per fibre ottiche |
JPS6170436A (ja) | 1984-09-14 | 1986-04-11 | Univ Kyoto | 円柱内屈折率分布測定方法 |
JP2678463B2 (ja) | 1988-06-16 | 1997-11-17 | 旭光学工業株式会社 | 屈折率測定方法 |
US5151752A (en) * | 1988-06-16 | 1992-09-29 | Asahi Kogaku Kogyo K.K. | Method of measuring refractive indices of lens and sample liquid |
US4934818A (en) | 1989-03-24 | 1990-06-19 | American Telephone And Telegraph Company | Refractive index profiling technique |
JPH03128411A (ja) | 1989-10-13 | 1991-05-31 | Hitachi Metals Ltd | 光学的形状測定装置 |
JP2892075B2 (ja) | 1990-01-31 | 1999-05-17 | 末三 中楯 | 屈折率分布、透過波面の測定方法およびこの方法に用いる測定装置 |
US5309214A (en) * | 1991-09-17 | 1994-05-03 | Olympus Optical Co., Ltd. | Method for measuring distributed dispersion of gradient-index optical elements and optical system to be used for carrying out the method |
JP3264469B2 (ja) * | 1993-12-07 | 2002-03-11 | 富士写真フイルム株式会社 | 光散乱媒体の屈折率分布情報の計測装置 |
JP3287517B2 (ja) | 1994-06-29 | 2002-06-04 | 富士写真光機株式会社 | 干渉縞による測定方法および装置 |
JP3699810B2 (ja) | 1997-07-29 | 2005-09-28 | 株式会社リコー | 屈折率分布の測定方法及び装置 |
US6765661B2 (en) | 2001-03-09 | 2004-07-20 | Novartis Ag | Lens inspection |
SE0302676D0 (sv) | 2003-10-09 | 2003-10-09 | Neural Ab | Method and apparatus for holographic refractometry |
JP4319554B2 (ja) | 2004-01-14 | 2009-08-26 | 株式会社リコー | 屈折率分布の測定方法および測定装置 |
JP3762420B2 (ja) | 2004-12-28 | 2006-04-05 | 株式会社リコー | 屈折率分布の測定方法及び装置 |
JP4533158B2 (ja) * | 2005-01-19 | 2010-09-01 | キヤノン株式会社 | 画像処理装置、画像処理プログラムおよび屈折率分布測定装置 |
DE102005039679A1 (de) | 2005-08-22 | 2007-03-08 | Galatea Ltd. | Verfahren zum Bestimmen des Wertes eines Objekts |
US7594729B2 (en) * | 2007-10-31 | 2009-09-29 | Wf Systems, Llc | Wavefront sensor |
JP5008650B2 (ja) | 2008-12-25 | 2012-08-22 | キヤノン株式会社 | 屈折率分布計測方法及び屈折率分布計測装置 |
JP5328437B2 (ja) | 2009-03-25 | 2013-10-30 | キヤノン株式会社 | 透過波面測定方法、屈折率分布測定方法、光学素子の製造方法、及び透過波面測定装置 |
JP4968965B2 (ja) * | 2009-11-18 | 2012-07-04 | キヤノン株式会社 | 屈折率分布の計測方法および計測装置 |
JP4968966B2 (ja) * | 2009-12-07 | 2012-07-04 | キヤノン株式会社 | 屈折率分布の計測方法および計測装置 |
JP4895409B2 (ja) * | 2010-05-25 | 2012-03-14 | キヤノン株式会社 | 屈折率分布計測方法および屈折率分布計測装置 |
JP5008763B2 (ja) | 2010-12-03 | 2012-08-22 | キヤノン株式会社 | 屈折率分布計測方法、屈折率分布計測装置および光学素子の製造方法 |
-
2010
- 2010-05-25 JP JP2010119636A patent/JP5021054B2/ja not_active Expired - Fee Related
-
2011
- 2011-05-19 KR KR1020110047222A patent/KR101288876B1/ko active IP Right Grant
- 2011-05-20 CN CN201110132234.8A patent/CN102297758B/zh not_active Expired - Fee Related
- 2011-05-24 US US13/114,309 patent/US8525982B2/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH028726A (ja) * | 1988-06-27 | 1990-01-12 | Asahi Optical Co Ltd | 屈折率測定方法 |
JPH0560654A (ja) * | 1991-09-03 | 1993-03-12 | Sumitomo Electric Ind Ltd | 屈折率分布測定方法及び装置 |
JPH08304229A (ja) * | 1995-05-09 | 1996-11-22 | Ricoh Co Ltd | 光学素子の屈折率分布の測定方法および装置 |
JP2000258292A (ja) * | 1999-03-09 | 2000-09-22 | Ricoh Co Ltd | 屈折率分布の測定装置及び屈折率分布の測定方法 |
JP2003021577A (ja) * | 2001-07-09 | 2003-01-24 | Ricoh Co Ltd | 屈折率分布測定方法および装置 |
US20060159332A1 (en) * | 2005-01-19 | 2006-07-20 | Yasuhiro Sawada | Image processing apparatus and refractive index distribution measuring apparatus |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104662402A (zh) * | 2012-09-28 | 2015-05-27 | 诺华股份有限公司 | 用于眼透镜的屈光力的自动化线内确定的方法 |
CN105339778A (zh) * | 2013-06-28 | 2016-02-17 | 佳能株式会社 | 折射率测量方法、折射率测量装置及光学元件制造方法 |
CN104614342A (zh) * | 2015-02-03 | 2015-05-13 | 清华大学 | 一种高温气流扰动下空气折射率三维重构测量方法 |
CN104614342B (zh) * | 2015-02-03 | 2017-04-19 | 清华大学 | 一种高温气流扰动下空气折射率三维重构测量方法 |
CN106226032A (zh) * | 2015-06-02 | 2016-12-14 | 佳能株式会社 | 测量方法、测量装置和光学元件的制造方法 |
Also Published As
Publication number | Publication date |
---|---|
KR101288876B1 (ko) | 2013-07-23 |
US20110292379A1 (en) | 2011-12-01 |
JP5021054B2 (ja) | 2012-09-05 |
CN102297758B (zh) | 2014-04-09 |
KR20110129343A (ko) | 2011-12-01 |
JP2011247692A (ja) | 2011-12-08 |
US8525982B2 (en) | 2013-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102297758B (zh) | 折射率分布测量方法和折射率分布测量设备 | |
CN101762376B (zh) | 折射率分布测量方法和折射率分布测量设备 | |
CN102564738B (zh) | 折射率分布测量方法、装置和制造光学元件的方法 | |
CN102519388B (zh) | 大口径离轴凸非球面反射镜面形检测方法 | |
CN101709955B (zh) | 子孔径拼接干涉检测光学非球面面形的装置 | |
KR100972571B1 (ko) | 비구면 렌즈의 면 어긋남 측정 방법 및 장치 | |
CN102062677B (zh) | 折射率分布测量方法和折射率分布测量装置 | |
CN107782254B (zh) | 一种混合补偿式子孔径拼接面形检测方法 | |
CN102200432A (zh) | 非球面物体测量方法和设备 | |
US7619191B1 (en) | Increase spatial sampling for wave front mid-spatial frequency error recovery | |
CN104685317A (zh) | 衍射光学元件和干涉测量方法 | |
CN104570621B (zh) | 一种双光束曝光***中光栅衍射波面误差的反馈调节方法 | |
CN106441816B (zh) | 计算全息法测量长焦距透镜透射波前的检测装置及检测方法 | |
CN103335982A (zh) | 利用波长调谐移相干涉仪测量平行平板光学均匀性的方法 | |
CN102954768A (zh) | 表面轮廓侦测装置、其对位方法及量测数据撷取方法 | |
CN104913843A (zh) | 一种微孔衍射波前质量的测量装置与方法 | |
JP5575161B2 (ja) | 屈折率分布計測方法および屈折率分布計測装置 | |
CN105352915A (zh) | 一种折射率二维分布的动态测量方法 | |
CN101858788B (zh) | 干涉楔板测量星间激光通信出射光波波高的方法和装置 | |
US7307706B2 (en) | Optical testing method and apparatus | |
CN102927930A (zh) | 采用平行光管拼接检测超大口径反射镜面形误差的方法 | |
CN108267094A (zh) | 一种基于旋转cgh的非圆柱面干涉拼接测量***及方法 | |
CN104914071A (zh) | 折射率分布测量方法、测量装置及制造光学元件的方法 | |
CN205642307U (zh) | 一种长程面形测量仪 | |
CN103353284A (zh) | 一种空气隙型偏光棱镜胶合误差的测量装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140409 Termination date: 20200520 |