CN101236067B - 用多波长来测量表面形状测量方法和使用此方法的装置 - Google Patents

用多波长来测量表面形状测量方法和使用此方法的装置 Download PDF

Info

Publication number
CN101236067B
CN101236067B CN2008100057818A CN200810005781A CN101236067B CN 101236067 B CN101236067 B CN 101236067B CN 2008100057818 A CN2008100057818 A CN 2008100057818A CN 200810005781 A CN200810005781 A CN 200810005781A CN 101236067 B CN101236067 B CN 101236067B
Authority
CN
China
Prior art keywords
pixel
interference fringe
measuring object
light
wavelength
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
Application number
CN2008100057818A
Other languages
English (en)
Other versions
CN101236067A (zh
Inventor
北川克一
杉山将
小川英光
铃木一嘉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tokyo Institute of Technology NUC
Toray Engineering Co Ltd
Original Assignee
Tokyo Institute of Technology NUC
Toray Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tokyo Institute of Technology NUC, Toray Engineering Co Ltd filed Critical Tokyo Institute of Technology NUC
Publication of CN101236067A publication Critical patent/CN101236067A/zh
Application granted granted Critical
Publication of CN101236067B publication Critical patent/CN101236067B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/02007Two or more frequencies or sources used for interferometric measurement
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2441Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/30Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
    • G01B11/303Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces using photoelectric detection means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/30Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
    • G01B11/306Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces for measuring evenness
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/02007Two or more frequencies or sources used for interferometric measurement
    • G01B9/02009Two or more frequencies or sources used for interferometric measurement by using two or more low coherence lengths using different or varying spectral width
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02015Interferometers characterised by the beam path configuration
    • G01B9/02017Interferometers characterised by the beam path configuration with multiple interactions between the target object and light beams, e.g. beam reflections occurring from different locations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02015Interferometers characterised by the beam path configuration
    • G01B9/02029Combination with non-interferometric systems, i.e. for measuring the object
    • G01B9/0203With imaging systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02015Interferometers characterised by the beam path configuration
    • G01B9/02032Interferometers characterised by the beam path configuration generating a spatial carrier frequency, e.g. by creating lateral or angular offset between reference and object beam

Abstract

本发明提供一种用多波长来测量表面形状测量方法和使用该方法的装置。通过以对光的行进方向倾斜任意角度的姿势配置基准面,利用从测量对象面和基准面返回同一光路的反射光产生干涉条纹。这时,利用不同的多波长单色光,由摄像装置拍摄各单色光产生的干涉条纹的各像素的强度值。CPU利用求干涉条纹波形的表达式,对每一计算对象像素利用各像素的强度值和处在其邻近的像素的强度值,假设各像素包含的干涉条纹波形的直流分量、交流振幅和相位相等,对每一波长求出各像素的干涉条纹波形的相位候补值群。进而,从多个候补值群求出共同的表面高度。

Description

用多波长来测量表面形状测量方法和使用此方法的装置
技术领域
本发明涉及用多波长来测量表面形状测量方法和使用该方法的装置,该用多波长来测量表面形状测量方法利用波长不同的多个单色光,测量半导体晶片、液晶板、等离子体显示板、磁膜、玻璃衬底、金属模等的具有平坦度的测量对象物的凹凸。
背景技术
以往,测量对象面的表面形状的测量方法如下地进行。按个别定时输出不同波长的第1单色光和第2单色光,用分束镜将各单色光照射到作为测量对象物的晶片的表面和基准镜上,使从两者反射后返回的反射光又汇聚到分束镜,通过同一光路,产生干涉条纹。这时,用2维图像检测器检测出每一单色光的反射光。此检测出的反射光呈现为干涉条纹,并且在表面的高低差部分,干涉条纹产生偏移。即,根据由单色光产生的连续的干涉条纹的偏移部分中两个偏移量接近的部分的该偏移量,求出晶片的图案的高低差(参考日本国特开2002-340524号公报)。
然而,已有的方法存在如下问题。
拍摄测量对象物的表面高低差图像时,由于陡峭的高低差部分在显示画面上显示的干涉条纹部分产生偏移,因此能判断是否有高低差。
然而,高低差状态未知的情况下,不能判断干涉条纹偏移的部分的高低差是凸状还是凹状。因而,不能准确求出测量对象面的陡峭处的边缘部分。
此外,如果按不同的定时照射不同波长的单色光,则由于不能检测出各自的干涉条纹,因此为了对测量对象物的表面高低差和形状进行测量而花费时间。也就是说,如果不能使单色光连续扫描测量对象物的整个表面,就不能使单色光也遍及多个测量对象物进行连续扫描。
发明内容
本发明是着眼于这种实况而完成的,其目的在于提供一种能高速且高精度地测量处在测量对象物的表面的凹凸高低差的多波长的表面形状测量方法和使用该方法的装置。
为了达到此目的,本发明采用如下组成。
一种多波长的表面形状测量方法,通过分支单元对测量对象面和基准面照射单色光,根据因从测量对象面和基准面两者反射并返回同一光路的反射光而产生的干涉条纹的强度值,求出测量对象的表面高度和表面形状,所述方法包含下列过程:
第1过程,该第1过程中对光的行进方向以任意角度的倾斜姿势配置所述基准面,将波长不同的多个单色光同时照射到测量对象物和基准面,并取得其产生的干涉条纹的图像;
第2过程,该第2过程中每一单色光求出取得的所述图像中各像素的干涉条纹的强度值;
第3过程,该第3过程中利用求干扰条纹波形的表达式,对每一所述像素利用各像素的强度值及其附近的多个像素的强度值,并假设这些像素的干涉条纹波形的直流分量、交流振幅和相位相等,对每一单色光求出各像素的相位;
第4过程,该第4过程中根据对每一单色光求出的各像素的相位求出表面高度的候补群,并从各波长的候补群求出共同的高度,作为实际高度;
第5过程,该第5过程中根据求出的所述实际高度,求出测量对象物的表面形状。
根据本发明的多波长的表面形状测量方法,将波长不同的多个单色光同时照射到测量对象物和基准面,并以对光的行进方向倾斜任意角度的姿势配置基准面,从而利用从测量对象面和基准面返回同一光路的反射光,对每一单色光产生干涉条纹。每一单色光以像素为单位求出此干涉条纹的强度值。然后,利用求干涉条纹波形的表达式,对各像素利用各像素的强度值和每一像素处在其邻近的像素的强度值,假设各像素包含的干涉条纹波形的直流分量、交流振幅和相位相等,对每一单色光求出各像素的相位。这时,能抵消各像素的直流分量和交流振幅,不需要施行去除空间频率分量用的低通滤波器处理。因而,能准确求出测量对象面的陡峭处的边缘部分,而不使空间分辨率降低。
根据此求出的相位,对每一单色光求测量对象物的表面高度的候补群,进而从各候补群求出共同的高度,作为实际高度。所以,与根据单一相位求表面高度相比,能从较大的候补范围高精度地求出实际高度,并能利用组合使用的波长提高能测量的高度的上限。
而且,由于能以同时输出多个单色光并同时检测出这些单色光组成的反射光的方式,对测量对象物的表面高度和表面形状进行测量,因此能得到同一条件下的测量结果。换句话说,不容易受振动等干扰的影响。还能谋求提高作业效率。
再者,上述方法中,最好一面使朝向测量对象物的光和1个或多个测量对象物相对平行移动、一面以规定的时间间隔对测量对象物的每一测量位置重复进行第1过程至第5过程,从而求出测量对象物的表面形状。
根据此方法,能一面对测量对象物的整个表面连续照射单色光、一面实时地求出测量对象物的表面高度和表面形状。
又,上述方法最好使各像素的强度值g(x,y)在像素附近,拟合干涉条纹波形表达式g(x,y)=a+bcos{2πfxx+2πfyy+φ},从而求出所求的每一所述波长的相位。
根据此方法,能较佳地实施上述第1发明。
又,上述方法最好利用配备将多个单色光分离的滤色片的摄像单元拍摄所述干涉条纹的图像,去除因滤色片特性产生的串扰的影响而包含在每一单色光中的其它单色光的干涉条纹的强度值。
根据此方法,能去除测量对象的单色光中包含的因滤色片特性而产生非所需的其它单色光的干涉条纹强度值。因而,能高精度地求出测量对象物的表面形状。
又,为了达到此目的,本发明也可为如下组成。
一种多波长的表面形状测量装置,通过分支单元对测量对象面和基准面照射单色光,根据因从测量对象面和基准面两者反射并返回同一光路的反射光而产生的干涉条纹的强度值,求出测量对象的表面高度和表面形状,所述装置包含下列组成要素:
保持单元,该保持单元以对光的行进方向倾斜任意角度的姿势配置所述基准面,而且放置并保持测量对象物;
照明单元,该照明单元同时输出波长不同的多个单色光;
摄像单元,该摄像单元利用照射多个所述单色光后从测量对象物和基准面反射并返回同一光路的反射光,使每一单色光产生干涉条纹,并拍摄测量对象面;
取样单元,该取样单元对每一像素装入所拍摄的所述测量对象面,作为干涉条纹的强度值;
存储单元,该存储单元存储作为由所述取样单元装入的所述强度值的干扰条纹强度值群;
运算单元,该运算单元对每一像素从所述存储单元存储的强度值群读出强度值,利用各像素的强度值和每一像素处在其邻近的像素的强度值,假设各像素包含的干扰条纹波形的直流分量、交流振幅和相位相等,并利用求干扰条纹波形的表达式,对每一单色光求出各像素的相位后,
根据每一单色光求出的各像素的相位求出多个表面高度候补群,从各候补群求出共同的高度作为实际高度,
进而,从该求出的所述测量对象面的表面高度,求出表面形状。
又可构成如下。
一种多波长的表面形状测量装置,通过分支单元对测量对象面和基准面照射单色光,根据因从测量对象面和基准面两者反射并返回同一光路的反射光而产生的干涉条纹的强度值,求出测量对象的表面高度和表面形状,所述装置包含下列组成要素:
保持单元,该保持单元以对光的行进方向倾斜任意角度的姿势配置所述基准面,而且放置并保持测量对象物;
照明单元,该照明单元输出多个波长组成的光;
分离单元,该分离单元将照射所述光后从测量对象物和基准面反射并返回同一光路的反射光分离为不同波长的多个单色光;
摄像单元,该摄像单元使分离的每一所述单色光产生干涉条纹,并拍摄测量对象面;
取样单元,该取样单元对每一像素装入所拍摄的所述测量对象面,作为干涉条纹的强度值;
存储单元,该存储单元存储作为由所述取样单元装入的所述强度值的干扰条纹强度值群;
运算单元,该运算单元对每一像素从所述存储单元存储的强度值群读出强度值,利用各像素的强度值和每一像素处在其邻近的像素的强度值,假设各像素包含的干扰条纹波形的直流分量、交流振幅和相位相等,并利用求干扰条纹波形的表达式,对每一单色光求出各像素的相位后,
从根据每一单色光求出的各像素的相位进行换算并求出的多个表面高度候补群,求出共同的高度作为实际高度,
进而,从该求出的所述测量对象面的表面高度求出表面形状。
根据上述多波长的表面形状测量装置,保持单元保持测量对象物。照明单元同时输出波长不同的多个单色光。摄像单元利用照射多个所述单色光后从测量对象物和基准面反射并返回同一光路的反射光,使每一单色光产生干涉条纹,并拍摄测量对象面。取样单元对每一像素装入所拍摄的所述测量对象面,作为干涉条纹的强度值。存储单元存储作为由取样单元装入的强度值的干扰条纹强度值群。运算单元对每一像素从存储单元存储的强度值群读出强度值,利用各像素的强度值和每一像素处在其邻近的像素的强度值,假设各像素包含的干扰条纹波形的直流分量、交流振幅和相位相等,并利用求干扰条纹波形的表达式,对每一单色光求出各像素的相位后,根据每一单色光求出的各像素的相位求出的多个表面高度候补群,从各候补群求出共同的高度作为实际高度,进而,从该求出的所述测量对象面的表面高度求出表面形状。
即,能将波长不同的多个单色光同时照射到测量对象物和基准面,并且每一不同的波长从返回同一光路的反射光求出以像素为单位的多个表面高度候补群。进而,从各候补群求出共同的高度,作为实际高度。所以,能从大的候补范围求出准确的表面高度。其结果,能从求得的表面高度求出测量对象物的表面高度。也就是说,能较佳地实现上述第1发明。
再者,上述组成中,可将照明单元在结构上构成例如具有输出不同的波长的多个单色光源,又可构成包含白色光源和分离成白色光源输出的光中规定的不同的多个波长的单色光并使其朝向分支单元的光学单元。
根据此组成,能实现容易同时将不同的多个波长的单色光照射到测量对象物和基准面。
而且,也可构成照明单元利用输出由多个波长组成的光的单元,并且在摄像单元内或其前配置将该光分离为波长不同的多个单色光的分离单元(例如滤色片)。这时,运算单元去除因滤色片特性产生的串扰的影响而包含在每一单色光中的其它单色光的干涉条纹的强度值,则更好。根据此组成,能高精度地求出测量对象的单色光的干涉条纹强度值,进而能高精度地求出测量对象物的表面形状。
又,所述运算单元最好使各像素的强度值g(x,y)在像素附近,拟合干涉条纹波形表达式g(x,y)=a+bcos{2πfxx+2πfyy+φ},从而求出所求的各像素的相位。
附图说明
为了说明发明,附图中示出目前认为较佳的若干形态,但应理解发明不受图中所示组成和单元限制。
图1是示出本实施例的表面形状测量装置的概略组成的图;
图2是示出表面形状测量装置中的处理的流程图;
图3是示出测量对象面的拍摄图像数据的图;
图4是示出拍摄图像的X轴方向亮度变化的图;
图5是示出利用sinφ和cosφ的符号信息能确定φ的范围的图;
图6是示出实际表面高度的图;
图7是示出利用本实施例测量陡峭高低差时的测量结果的图;
图8是示出变换例的表面形状测量装置中的处理的流程图;
图9是示出空间频率的估计处理的流程图;
图10是示出空间频率的估计处理的模式图;
图11是示出空间频率的估计实验结果的图;
图12是示出利用变换例方法的实测结果的图;
图13是示出彩色滤色片特性的图;
图14是示出蓝光源发光时产生的串扰的图;
图15是示出红光源发光时产生的串扰的图;
图16是示出绿光源发光时产生的串扰的图;
图17是示出串扰的非校正亮度数据和校正亮度数据的图;
图18是示出串扰的非校正相位数据的图;
图19是示出串扰的校正相位数据的图。
具体实施方式
下面,参照附图说明本发明的实施例。再者,本实施例中,取利用干涉条纹测量表面形状实质上平坦的测量对象物的表面高度和表面形状为例进行说明。
图1是示出本发明实施例的表面形状测量装置的概略组成的图。
此表面形状测量装置配备对半导体晶片、液晶板、等离子体显示板、磁膜、玻璃衬底、金属膜等在表面具有凹凸高低差的实质上平坦的测量对象物30照射特定波长区的单色光的光学系统单元1、控制光学系统单元1的控制系统单元2、以及放置并保持测量对象物30的保持台40。
光学系统单元1配备往测量对象面30A和基准面15输出不同的多个波长的单色光的照明装置10、使各单色光成为平行光的准直透镜11、将两个单色光反射到测量对象物30的方向而使来自测量对象物30的方向的光通过的半透明反射镜13、汇聚半透明反射镜13反射来的单色光的物镜14、将通过物镜14的单色光分为反射到基准面15的基准光和通过测量对象面30A的测量光并再汇聚基准面15反射的基准光和测量对象面30A反射的测量光从而产生干涉条纹的分束镜17、使汇聚基准光和测量光后得到的单色光成像的成像透镜18、以及连同干涉条纹一起拍摄测量对象面30A的摄像装置19。
照明装置10的组成包含输出不同的2个波长的光的第1光源10A和第2光源10B、以及使从不同方向输出的两个光一致往相同方向的光学构件10C。作为本实施例的各光源10A、10B,利用例如LED(发光二极管),输出例如下列波长的光。第1光源10A的波长λ1=470nm(纳米),第2光源10B的波长λ2=627nm(纳米)。再者,照明装置10相当于本发明的照明单元。
半透明反射镜13将来自准直透镜11的平行光往测量对象物30反射,另一方面却使从测量对象物30返回的反射光通过。
物镜14是将入射来的两个单色光汇聚到作为焦点的测量对象面的透镜。
分束镜17,将物镜14汇聚的光分成在基准面15上反射的基准光和在测量对象面30A上反射的测量光。而且,又对各面上反射并返回同一光路的基准光和测量光进行汇聚,从而产生干涉条纹。再者,分束镜17相当于本发明的分支单元。
基准面15,其表面受到镜面加工,以对基准光行进方向前后倾斜的姿势加以安装。此基准面15反射的基准光到达分束镜17后,此基准光又被分束镜17反射。
再者,通过以对基准光行进方向前后倾斜的姿势安装基准面15,使基准光的到达距离和到达摄像装置19前的距离因其反射面的位置而变化。这等效于:移动基准面15,从而使基准面15与分束镜17之间的距离L1变动。
也就是说,通过分束镜17的测量光被往焦点汇聚,并在测量对象面30A上进行反射。此反射的测量光到达分束镜17,并通过该分束镜17。
分束镜17又汇聚基准光和测量光。这时,由于基准面15与分束镜17之间的距离L1和分束镜17与测量对象面30A之间的距离L2不同,产生光路差。根据此光路差,基准光和测量光产生干涉。
摄像装置19,拍摄由测量光映出的测量对象面30A的图像。这时,由于基准面15倾斜,拍摄的测量对象面30A的图像中拍摄干涉条纹,即干涉造成的亮度在空间上的变动。由控制系统单元2的存储器21收集此拍摄的图像数据。而且,后文将说明,结构上构成利用控制系统单元2的驱动部24,将光学系统单元1往图1中的x、y、z轴方向移动到希望的拍摄部位。又,由摄像装置19在规定的取样定时拍摄测量对象面30A、30B的图像,并由控制系统单元2收集其图像数据。再者,摄像装置19相当于本发明的摄像单元,控制系统单元2作为本发明的取样单元起作用。
作为本实施例的摄像装置19,构成能检测出波长不同的多个单色光即可,例如有:CCD固体摄像元件、MOS图像传感器、CMOS图像传感器、光电摄像管、雪崩电子倍增效应摄像管、EB-CCD等。
控制系统单元2配备进行整个表面形状测量装置的集中控制和规定的运算用的CPU20、存储由CPU20逐次收集的图像数据或运算结果等各种数据和程序等的存储器21、输入取样定时或拍摄区等其它设定信息的鼠标器或键盘等输入部22、以及显示测量对象面30A的图像等的监视器23。而且,根据CPU20的指示,驱动光学系统单元1,使其上下左右移动。例如,由配备包含2轴驱动型伺服电机等驱动机构的驱动部24的计算机系统构成。再者,CPU20相当于本发明的运算单元。
CPU20是“中央运算处理装置”,对摄像装置19、存储器21和驱动部24进行控制,并配备根据包含摄像装置19拍摄的干涉条纹的测量对象面30A的图像数据进行求出测量对象物30的表面高度的运算处理的相位计算部25、从求出的多个表面高度数据求表面形状的图像数据编制部27。后文阐述此CPU20中的相位计算部25和图像数据编制部27的处理。CPU20还连接监视器23、键盘或鼠标器等输入部22。操作者一面观测监视器23上显示的操作画面、一面从输入部22进行各种设定信息的输入。而且,监视器23上将测量对象面30A的表面图像或凹凸形状等作为数值或图像进行显示。
驱动部24是例如使光学系统单元1往图1中的x、y、z轴方向移动到希望的拍摄部位的装置。此驱动部24包含配备根据来自CPU20的指示往x、y、z轴方向驱动光学系统单元1的、例如3轴驱动型伺服电机的驱动机构。再者,本实施例中,使光学系统单元1动作,但也可例如使放置测量对象物30的保持台40往正交3轴方向变动。而且,移动轴可为2轴或少于2轴或者不存在。
下面,按照图2所示流程图说明作为本实施例的特征部分的整个表面形状测量装置进行的处理。
再者,本实施例中,取图1那样将基准面15倾斜的情况为例进行说明。此情况下,拍摄图像如图3所示。再者,本实施例中,为了简化说明,取x轴方向的情况为例进行说明。也能作往包含y、z方向的3维方向移动的测量。
《步骤S1》取得测量数据
CPU20使未示出的步进电机等驱动系统驱动,从而驱动部24使光学系统单元1移动到测量对象物30的拍摄区。拍摄位置确定时,光学系统单元1从照明装置10的各光源10A、10B同时输出不同波长的单色光λ1、λ2。这两个单色光被光学部件10c汇聚后,朝向半透明反射镜13。
摄像装置19与单色光的输出联动地开动,例如进行1次图1所示具有凸部30B的测量对象面30A的拍摄。将此拍摄取得的测量对象面30A的干涉条纹的图像数据存储到存储器21。即,存储器21中对每一单色光存储倾斜姿势的基准面15上的反射光和受测量对象面30A反射并返回的反射光产生的干涉条纹的图像数据。这时在基准面15上反射的光的传播距离(L1的2倍)在基准面15的反射位置有规则地变动。因而,测量对象面30A的高度平坦的部分中,来自测量对象面30A的反射光的传播距离(L2的2倍)在测量部位无变动,所以摄像装置19拍摄的图像的干涉条纹根据基准面15的倾斜方向和角度,空间上有规则地呈现在拍摄面内。每次来自基准面15的反射光的传播距离(L1的2倍)与来自测量对象面30A的反射光的传播距离(L2的2倍)之差为λ1/2=235nm和λ2/2=313.5nm,呈现1周期份额的该干涉条纹。
另一方面,在如图1所示那样测量对象面30A的高度产生变动的部位中,呈现为干涉条纹偏移的不规则条纹图案。
再者,此过程相当于本发明中的第1过程。
《步骤S2》取得干涉光强度值
CPU20从图像数据装入存储器21中存储的各像素的强度值,即测量对象面30A的干涉光的强度值。这时,在测量对象面30A和凸部30B的高度变动的图4所示像素号200和330附近,呈现为干涉条纹空间相位偏移(例如图4的本实施例中,往X轴方向偏移)的不规则条纹图案。
再者,此过程相当于本发明中的第2过程。
《步骤S3》算出各波长以像素为单位的相位φ1、φ2
CPU20的相位计算部25使用测量对象面30A的计算对象的像素和该像素相邻的像素(本实施例中为x轴方向相邻的像素)各自的干涉条纹光强度值,利用预先决定的算法求出该像素的相位。具体而言,使计算对象的像素和该像素相邻的像素的干涉条纹的光的强度值拟合求干涉条纹波形的表达式,从而求出相位。
首先,如下面的式(1)那样记述检索对象的像素的干涉条纹的光的强度值。
g(x)=a(x)+b(x)cos{2πfx+φ(x)}      …(1)
这里,x是计算对象的像素的位置坐标,a(x)是干涉条纹波形中包含的直流分量,b(x)是干涉条纹波形中包含的交流分量(是振动分量的振幅,下文适当称为“交流振幅”),f是干涉条纹g(x)的空间频率分量,φ(x)是测量对象面30A的规定像素所对应的相位,是应计算的相位。再者,以2维(x、y)表现计算对象的像素的位置坐标,但本实施例为了简化说明,以省略y坐标的方式记述。
接着,关于相邻的像素,由于从计算对象的像素偏移微小距离Δx,如下面的式(2)那样表现其干涉条纹的光的强度值。
g(x+Δx)=a(x+Δx)+b(x+Δx)cos{2πf(x+Δx)+φ(x+Δx)}  …(2)
这里,本实施例中,计算对象的像素与相邻相素的间距为微小的距离,所以假设横跨各像素的干涉条纹中包含的直流分量、交流振幅和相位相等,并利用下面的关系式(3)~(5)。
a(x)=a(x+Δx)=a      …(3)
b(x)=b(x+Δx)=b      …(4)
φ(x)=φ(x+Δx)=a    …(5)
这里,a、b、φ是常数。
通过上述式(3)~(5)那样假设,式(1)和式(2)能置换成下面的式(1a)和式(2a)。
g(x)=a+bcos{2πfx+φ}   …(1a)
g(x+Δx)=a+bcos{2πf(x+Δx)+φ}  …(2a)
接着,对式(1a)和式(2a)进行变换,从而编制成下面的式(6)、式(7)。
G(x)=g(x)-a=bcos(2πfx+φ)  …(6)
G(x+Δx)=g(x+Δx)-a=bcos{2πf(x+Δx)+φ}  …(7)
接着,利用加法定理,将式(6)、式(7)变换成下面的式(8)、式(9)。
G(x)=bcos(2πfx+φ)
    =b{cos(2πfx)cosφ-sin(2πfx)sinφ}   …(8)
G(x+Δx)=bcos{2πf(x+Δx)+φ}
        =b[cos{2πf(x+Δx)}cosφ-sin(2πfx+Δx)sinφ]   …(9)
接着,用行列式(10)表示这些式(8)、式(9)。
G ( x ) G ( x + Δx ) = A b cos φ b sin φ · · · ( 10 )
再者,将A表示如下。
A = cos ( 2 πfx ) - sin ( 2 πfx ) cos { 2 πf ( x + Δx ) } - sin { 2 πf ( x + Δx ) }
这里,从行列式(10)的左边乘A的逆矩阵并展开,从而求出下面的式(11)、式(12)。
G ( x ) sin { 2 πf ( x + Δx ) } - G ( x + Δx ) sin ( 2 πfx ) sin ( 2 πΔfx ) = b cos φ · · · ( 11 )
G ( x ) cos { 2 πf ( x + Δx ) } - G ( x + Δx ) cos ( 2 πfx ) sin ( 2 πΔfx ) = b sin φ · · · ( 12 )
利用上述式(11)、式(12),能得到下面的式(13)。再者,这里将上述bsinφ和bcosφ分别取为bsinφ=S和bcosφ=C,又使tanφ=S/C。
φ=arctan{S/C}+n’π   …(13)
再者,n’是整数。
这里,CPU20还配备符号判断部26,此符号判断部26参照sinφ和cosφ的符号信息。使用此符号信息时,根据sinφ和cosφ的符号的组合,能将φ的存在范围从π扩展到2π。图5是参照式(13)所示那样的sinφ和cosφ的符号信息来确定φ的范围用的具体的图。因此,使用sinφ和cosφ的符号信息,则能用式(14)表示式(13)。
φ=arctan{S/C}+2nπ  …(14)
再者,n是整数。
所以,如果G(x)和干涉条纹波形的空间频率f已知,就能利用式(14)求出相位φ。G(x)包含像素的亮度信息g(x)和g(x+Δx)以及干涉条纹波形的直流分量a,因而其结果为:g(x)和g(x+Δx)、干涉条纹波形的直流分量a、干涉条纹波形的空间频率f已知,则能用式(14)求出φ。即,利用上述运算式求出波长λ1和λ2时各自的相位φ1、φ2。
能以作为摄像装置19的像素亮度信息的方式取得g(x)和g(x+Δx)。
能用例如取摄像装置19观测到的全部像素的平均值的方法、取相位计算对象像素的邻近像素的平均值的方法或预先测量反射率的方法等求出a。
能用例如根据基准面15的设置角度求出的方法、预先根据作为测量对象物观察平坦面时的干涉条纹波形在画面内的干涉条纹数求出的方法等求出f。
再者,此过程相当于本发明分第3过程。
《步骤S4》算出各波长的以像素为单位的表面高度z1、z2
CPU20将上述式(14)对每一波长λ1、λ2算出的计算对象像素的相位φ1(x)、φ2(x)代入到下面的式(15),求出各自的高度z1(x)、z2(x)。
z(x)=[φ(x)/4π]λ+z0  …(15)
再者,z0是测量对象物30的基准高度。
这里,取波长λ的情况下,每一λ/2的范围存在表面高度的解的候补值群。因而,本实施例那样利用波长λ1、λ2时的表面高度的解的候补值群周期性地存在于每一两个候补值群的λ1/2和λ2/2的最小公倍数的范围。
由于应求的表面高度仅1个,从两个候补值群求出共同的高度,作为实际高度。即,将在各候补值群求出的表面高度的解的候补值中相互最接近的高度当作实际高度。
例如,周期性地存在的λ1的解的候补值群为图6(a)那样,λ2的解的候补值群为图6(b)那样。这里,相当于图1所示测量对象物30的底部30A的部分为像素号的约0至200和330至500,凸部30B为中央的像素号200附近至330。因此,每一像素对两个候补值群进行比较,并提取底部30A和凸部30B的表面高度实质上一致的值。即,一方的底部30A由于波长λ1的高度z11与λ2的高度z21一致,将这时的值作为实际高度。另一方的凸部30B由于波长λ1的高度z13与λ2的高度z23一致,将这时的值作为实际高度。
也就是说,基于此原理,根据用波长λ1、λ2的测量的相位φ1(x)、φ2(x),利用下面的式(15a)、式(15b)求出各自的表面高度。
z1(x)=[φ1(x)/2π+n1]·(λ1/2)  …(15a)
z2(x)=[φ2(x)/2π+n2]·(λ2/2)  …(15b)
《步骤S5》算出以像素为单位的实际高度Z
又,利用上述2个公式,并由下面的式(15c)求出实际高度Z。
Z(x)=[Φ(x)/2π+N]·(Λ/2)  …(15c)
但是,Φ=φ12;(-π<Φ≤π),Λ=(λ1λ2)/|λ12|,N=n1-n2
再者,此过程相当于本发明的第4过程。
《步骤S6》对全部像素结束计算?
CPU20重复进行步骤S3~S6的处理,直到对全部像素结束相位和高度的计算。
《步骤S7》显示表面形状
CPU20的图像数据编制部27根据算出的实际高度的信息编制测量对象面30A、30B的显示图像。然后,CPU20根据此图像数据编制部27编制的信息,如图7所示,在监视器23上显示测量对象物30的表面高度的信息,或显示基于这些各特定部位的高度的信息的3维或2维图像。操作者可通过观察这些显示,掌握处在测量对象面30A的表面的凹凸形状。至此,测量对象面30A的表面形状测量处理结束。
再者,此过程相当于本发明的第5过程。
如上文所述,在从摄像装置19拍摄的图像数据算出每一像素的干涉条纹的光的强度值和其附近的多个像素的强度值的过程中,假设各像素的干涉条纹波形包含的直流分量a(x)、交流振幅b(x)和相位φ(x)各自对各像素相等,并进行联立比较,从而能抵消各像素的干涉条纹的直流分量和交流振幅。
因而,能不利用低通滤波器而测量测量对象面30A的表面高度,所以能高精度地求出测量对象面30A的陡峭边缘部分,如图7所示。结果,能高精度地测量测量对象面30A的表面形状。
又,根据求出的相位对每一单色光求测量对象面的表面高度候补值群,进而从各候补群求出共同的高度作为实际高度,所以与根据单一相位求表面高度相比,能从较大的候补范围高精度地求出实际高度。又能利用组合使用的波长提高能测量的高度的上限。例如,减小λ1、λ2之差,则能检测出更大的凹凸高低差。
而且,由于根据相位求出测量对象物30的表面高度,能判别表面的凹凸形状。
又,能以同时输出并同时检测出单色光组成的反射光的方式,对测量对象物30的表面高度和表面形状进行测量,所以能谋求提高作业效率。
本发明不限于上述实施例,也能变换实施如下。
(1)上述实施例中,利用与计算对象像素相邻的1个像素的干涉光强度值求出测量对象面30A的高度,但也可利用计算对象像素附近的2个像素,从共计3个像素求测量对象面30A的高度。
此情况下,解3元联立方程式,能求出3个未知变量。因此,与利用计算对象像素相邻的1个像素的干涉光强度值求解时相比,能添加干涉条纹波形的直流分量a或干涉条纹波形的空间频率f,当作未知数。
举用别的方法估计f并添加a作为未知数的例子,则利用计算对象像素x+Δx1、相邻相素x+Δx2和x+Δx3,共计3个像素的干涉光的强度值,如下面的式(16)那样求此规定像素的高度。
tan φ = ( g 2 - g 1 ) · cos ( x + Δ x 3 ) + ( g 1 - 2 g 2 + g 3 ) · cos ( x + Δ x 2 ) + ( g 3 - g 2 ) cos ( x + Δx 1 ) ( g 2 - g 1 ) · sin ( x + Δx 3 ) + ( g 1 - 2 g 2 + g 3 ) · sin ( x + Δx 2 ) + ( g 3 - g 2 ) sin ( x + Δx 1 ) · · · ( 16 )
但是,g1=g(x+Δx1),g2=g(x+Δx2),g3=g(x+Δx3)。
利用上述运算处理,能高精度地求出测量对象面30A的规定像素的表面高度。
又,本发明中,可利用不少于4个的多个邻近像素求出测量对象面30A的高度。这时,由于根据多个像素的强度值的信息求φ(x),具有能减小拍摄时的亮度噪声或计算时的量化误差对测量值的影响的特点。
即,上述实施例中,CPU20的相位计算部25利用预先决定的算法求出所求测量对象面30A的像素x的相位,以便求出该像x和该像素x附近的多个(N个)像素x+Δxi(I=1、2、3、……、N)(本实施例中存在于x轴方向的多个像素)各自的干涉条纹的光强度值。具体而言,由下面的式(17)求出计算对象像素x的干涉条纹的光强度值。
g(x+Δxi)=a(x+Δxi)+b(x+Δxi)cos{2πf(x+Δxi)+φ(x+Δxi)}  …(17)
此情况下,也与上述实施例相同,假设各像素的干涉条纹中包含的直流分量、交流振幅和相位相等,并应用上述式(3)~式(5)。本实施例的情况下,成为下面的式(18)~式(20)那样。
a(xi)=a(x+Δxi)=a     …(18)
b(xi)=b(x+Δxi)=b     …(19)
φ(xi)=φ(x+Δxi)=φ  …(20)
通过上述式(18)~式(20)那样假设,能将式(17)表示成下面的式(21)。
g(x+Δxi)=a+bcos{2πf(x+Δxi)+φ}
=a+bcosφ·cos{2π·f·(x+Δxi)}-bs inφ·sin{2π·f·(x+Δxi)}  …(21)
从这些N元联立方程式,以基于拟合的估算求出干涉条纹波形的直流分量a.、交流振幅b、空间频率f、相位φ。这里,与从2像素或3像素算出相位φ时相同,也可用别的方法求出干涉条纹波形的直流分量a和干涉条纹波形大空间频率f。
接着,记述估计干涉条纹波形的空间频率f并以基于拟合的估算求出干涉条纹波形的直流分量a的例子。
尤其是本变换例中,在上述实施例的图2所示流程图的步骤S2与步骤S3之间进行空间频率f估计的运算处理。即,此变换例的总体流程图如图8所示。因而,首先说明步骤S3中的空间频率f的估计的算法。
《步骤S3》估计空间频率f
首先,空间频率f的估计中利用例如Prony法。即,根据预先取得的多个等间隔取样值一面修改取样点频率fs、一面逐次进行取样化。最后,又收敛到fs=4f,又估出空间频率f。具体而言,利用下面的算法。
根据预先取得的多个等间隔取样值,利用等间隔取样点x0、x1、x2、x3这4个点,从下面的式(22)求出空间频率f。
f=(fs/2π)acos[(x3-x2+x1-x0)/{2*(x2-x1)}]   …(22)
再者,acos为arccos(反余弦)。
这里,Prony法的空间频率估计精度高达取样频率fs近似于作为估计对象的空间频率f的4倍(fs=4f)。因而,一面逐次修改取样频率fs、一面使其收敛到fs=4f。
再者,取样点间隔t只能取1个像素的整数倍,所以将取样频率fs与得到的空间频率f之比最接近4时的频率估计为空间频率f。
本实施例中,利用上述算法,沿图9所示流程图执行空间频率f的估计如下。
《步骤S10》将取样频率的暂定值fs’设定为初始值f1
将对预先取得的取样点数据的取样频率暂定值fs’设定为f1。具体而言,将取样频率f1设定得大于2倍估计频率,将取样点间隔t设定成小于奈奎斯特间隔。例如,如图10所示,在取样点间隔t收敛到不大于二分之一周期T的范围,设定暂定取样频率f1
《步骤S11》算出取样点间隔t
从公式t=round(1/fs’)算出取样点间隔t。再者,round的含义为基于四舍五入的整数化。
《步骤S12》判断收敛
对上次算出的间隔t和新算出的间隔t进行比较运算,在算出的值相同时,在该时间点判断为收敛。例如,如图11的红条纹的实验数据所示,重复进行4次计算的情况下,取样点间隔t按6、9、10、10的顺序变化。即,在第4次计算中判断为收敛。这样在新算出的间隔t为与上次的间隔t相同的值时,结束本处理,进至步骤S4。如果新的间隔t与上次不同,进至下面的步骤S13。
《步骤S13》算出取样频率fs
确定间隔t,则以fs=1/t的方式算出取样频率fs。
《步骤S14》算出空间频率f
对测量对象区域内的每一像素,以间隔t选择取样点x0、x1、x2、x3这4个点,并利用式(22)算出空间频率f,从而计算整个区域的平均值。
《步骤S15》算出新取样频率的暂定值fs’
求出空间频率f,则利用fs=4f的关系算出新取样频率的最大值fs’。一求出该最大值fs’就返回步骤S11,重复其后的步骤,直到步骤S12中上次的间隔t与新的间隔t一致为止。
接着,结束空间频率f的估计,则对式(21),根据多个像素的强度值g(x+Δxi),求满足下面的式(23)的a、bcosφ、bsinφ这组数据。从其中的bcosφ、bsinφ可求出φ。
min Σ i = 1 N [ g ( x + Δx i ) - [ a + b cos φ · cos { 2 π · f · ( x + Δx i ) ) } - b sin φ · sin { 2 π · f · ( x + Δ x i ) } ] ] 2 · · · ( 23 )
这里,如下面的式(24)、式(25)、式(26)那样计算。
A = 1 cos { 2 π · f · ( x + Δx 1 ) + φ ) } - sin { 2 π · f · ( x + Δ x 1 ) } 1 cos { 2 π · f · ( x + Δx 2 ) + φ ) } - sin { 2 π · f · ( x + Δ x 2 ) } · · · · · · · · · 1 cos { 2 π · f · ( x + Δx N ) + φ ) } - sin { 2 π · f · ( x + Δ x N ) } · · · ( 24 )
G = g ( x + Δx 1 ) g ( x + Δx 2 ) · · · g ( x + Δx N ) · · · ( 25 )
α = a b cos φ b sin φ · · · ( 26 )
接着,可如下面的式(27)那样求出满足式(23)的α。
α=(AT·A)-1·AT·G    …(27)
于是,可用下面的式(28)从矢量分量α求出相位φ。
tan φ = b sin φ b cos φ · · · ( 28 )
可从此式(28)与上述实施例所示式(13)、式(14)同样地求出φ。
如上文所述,对计算对象像素利用多个该像素的相邻像素,也能高精度地求出该像素的高度。再者,此变换实施例中,利用处在X轴方向的多个相邻相素,但对其利用Y轴方向的像素等的像素利用无专门限制。
(2)上述实施例中,利用分布在X轴方向或Y轴方向的1个轴的邻近像素进行说明,但也可利用分布在XY平面上的邻近像素。此情况下,可对作为计算对象处在x、y坐标上的像素将处在该像素附近的像素的坐标取为xi、yi(i=1、2、……、N)进行解方程,并可将上述式(1)取为g(x、y)=a+bcos{2πfxx+2πfyy+φ(x)}作运算处理。
(3)上述实施例中,以静止状态拍摄测量对象物30,但也可构成一面使长测量对象物或多个测量对象物30以规定速度移动、一面取与此移动速度同步,并在规定取样时间拍摄测量对象面30A的图像,求出表面高度。
(4)上述实施例中,也可在从每一单色光求出的表面高度的候补值群算出实际高度的过程中,算出并利用候补值的平均。利用这点,能抵消装置等的个体误差。
(5)上述实施例中,也可照明装置10的光源利用白色光源,并且在白色光到达测量对象面和基准面前的光路上,配置从白色光提取不同波长的单色光的光学单元。而且,也可在反射光到达摄像装置30前的光路上配置同样的光学单元。还可构成用单独的摄像装置19拍摄各单色光。
(6)上述实施例中,求出1个测量对象物30的表面高度和表面形状,但也可构成如下。例如,一面使多个测量对象物30在传送路径上连续传送或使定向配置在活动台上的多个测量对象物30在X-Y平面上移动、一面对全部测量对象物30求出其表面高度和表面形状。
(7)上述实施例中,在预先保持测量对象物30的平行度的状态下,任意设定基准面15的角度,并进行测量,但也可构成如下。例如,使测量对象物30的测量面方保持平坦度并设置高度已知的基准区,在设定基准面15的角度后,预先测量此区的高度,而且测量这时的测量对象物30的斜率。于是,结构上可构成算出此求出的斜率的校正量,并利用其结果进行校正。利用这点,能估计干涉条纹波形的空间频率分量。
(8)上述实施例中,光源使用波长不同的2个LED,但可用例如波长λ=627nm的红(R)、波长λ=530nm的绿(G)和波长λ=470nm的蓝(B)组成的RGB-LED,以代替该LED。此情况下,摄像装置19利用配备例如图13所示那样可将来自这3个光源的反射光分离为单色光的滤色片的彩色相机。再者各LED的波长不限于上述值。
滤色片取决于其特性,但如图13所示,分离单色光中包含其它单色光的频段的光。也就是说,产生串扰。
因此,该变换例那样使用多个频段不同的单色光的情况下,最好去除实际测量时产生的串扰。
利用实验或模拟预先求出去除各分离单色光包含的其它单色光所需的校正系数,并利用该校正系数校正由观测求出的每一像素的强度值,从而进行串扰的去除。
例如,利用下面的方法能求出校正系数。
首先,说明串扰的模型式和校正式。照射RGB-LED光源而取得的图像中,根据各光源个别照射的亮度的加和性,其亮度信号I能用下面的式(29)表示。
I(R、G、B)=I(R)+I(G)+I(B)   …(29)
这里,考虑从各像素的观测亮度B’、G’、R’求出实际亮度B、G、R。此情况下,由于串扰影响各光源的亮度,各单色光的观测亮度B’、G’、R’包含其它光源的亮度。因而,观测亮度与实际亮度的关系可用下面的模型式(30a)~(30c)表示。
B’=B+aG+bR    …(30a)
G’=cB+G+dR    …(30b)
R’=eB+fG+R    …(30c)
这里,式(30a)~(30c)的a~f是表示串扰的大小的系数。
根据该式(30a)~(30c),各光源所对应的实际亮度在各系数小的情况下,系数之积的项能忽略,可由下面的式(31a)~(31c)近似地求出该亮度。
B=B’-aG’-bR’…(31a)
G=G’-cB’-dR’…(31b)
R=R’-eB’-fG’…(31c)
接着,说明求串扰的大小a~f的方法。
首先,个别点亮各光源。这时,对多个像素求输出到监视器的RGB的亮度。例如,求仅使绿(G)LED点亮时的RGB亮度。同样,使红(R)和蓝(B)LED个别点亮,对多个像素求各自的RGB输出。
这里,从得到的RGB输出求出图14至图16所示的相关图,并根据每一观测光源包含的与其它光源的相关关系求出校正系数。
例如,仅使蓝(B)LED点亮时,如图14所示,蓝(B)至绿(G)的串扰的校正系数为0.23,蓝(B)至红(R)的串扰的校正系数为0.00。同样,仅使红(R)LED点亮时,如图15所示,红(R)至蓝(B)的串扰的校正系数为0.00,红(R)至绿(G)的串扰的校正系数为0.04。仅使绿(G)LED点亮时,如图16所示,绿(G)至蓝(B)的串扰的校正系数为0.08,绿(G)至红(R)的串扰的校正系数为0.13。
又,作为其它方法,对同时照射RGB-LED的全部光源而得到的图像应用上述方法,也能得到校正系数。但是,此情况下,需要非常多的亮度数据。
使3波长的LED点亮并拍摄测量对象物的图像的结果,如图17所示,得到校正处理前和校正处理后的亮度分布。即,得到的结果为:虚线所示的校正前亮度分布在相当于测量对象物的高低差部分的x坐标210~240附近变化小,但校正后的亮度变化大,其程度即使目视也容易读取。
同样,利用上述亮度分布的数据求以像素为单位的相位时,得到下列结果。即,利用图18所示的校正前的分布数据的情况下,即使无高低差的部分也成为波状起伏的形状。与此相反,校正后,如图19所示,高低差部分210~240以外,形成实质上平坦,仅在高低差部分看到大变化。
因而,根据本实施例,即使利用相邻频段的多个单色光,也能去除串扰造成的其它单色光的非所需亮度(光强度值)。
本发明能以其它具体方式实施而不脱离其思想和本质,因而作为发明范围的指示,应参照所附权利要求书,而非上述说明。

Claims (10)

1.一种用多波长来测量表面形状的测量方法,通过分支单元对测量对象面和基准面照射单色光,根据因从测量对象面和基准面两者反射并返回同一光路的反射光而产生的干涉条纹的强度值,求出测量对象面的表面高度和表面形状,其特征在于,所述测量方法包含下列过程:
第1过程,该第1过程中以对光的行进方向倾斜任意角度的姿势配置所述基准面,将波长不同的多个单色光同时照射到测量对象物和基准面,并取得其产生的干涉条纹的图像;
第2过程,该第2过程中对每一单色光求出取得的所述图像中各像素的干涉条纹的强度值;
第3过程,该第3过程中利用对干涉条纹波形进行求取的表达式,对所述各像素利用计算对象的像素的强度值及该计算对象的像素附近的多个像素的强度值,并假设计算对象的像素及该计算对象的像素附近的多个像素的干涉条纹波形的直流分量、交流振幅和相位相等,使计算对象的像素的强度值g(x,y)在该计算对象的像素附近,拟合干涉条纹波形表达式g(x,y)=a+bcos{2πfxx+2πfyy+φ},从而对每一单色光求出所述各像素的相位,这里,(x,y)是计算对象的像素的位置坐标,a是干涉条纹波形中包含的直流分量,b是干涉条纹波形中包含的交流振幅,f是干涉条纹的空间频率分量,φ是应计算的相位;
第4过程,该第4过程中根据对每一单色光求出的各像素的相位求出表面高度的候补群,并从各波长的候补群求出共同的高度,作为实际高度;
第5过程,该第5过程中根据求出的所述实际高度,求出测量对象物的表面形状。
2.如权利要求1中所述的用多波长来测量表面形状的测量方法,其特征在于,
一面使朝向所述测量对象物的光和1个或多个测量对象物相对平行移动、一面以规定的时间间隔对测量对象物的每一测量位置重复进行所述第1过程至第5过程,从而求出测量对象物的表面形状。
3.如权利要求1或2所述的用多波长来测量表面形状的测量方法,其特征在于,
利用配备将多个单色光分离的滤色片的摄像单元,拍摄所述干涉条纹的图像,
去除因所述滤色片特性产生的串扰的影响而包含在每一所述单色光中的其它单色光的干涉条纹的强度值。
4.一种用多波长来测量表面形状的测量装置,通过分支单元对测量对象面和基准面照射单色光,根据因从测量对象面和基准面两者反射并返回同一光路的反射光而产生的干涉条纹的强度值,求出测量对象面的表面高度和表面形状,其特征在于,所述测量装置包含下列组成要素:
保持单元,该保持单元以对光的行进方向倾斜任意角度的姿势配置所述基准面,而且放置并保持测量对象物;
照明单元,该照明单元同时输出波长不同的多个单色光;
摄像单元,该摄像单元利用照射多个所述单色光后从测量对象物和基准面反射并返回同一光路的反射光,使每一单色光产生干涉条纹,并拍摄测量对象面;
取样单元,该取样单元对每一像素装入所拍摄的所述测量对象面,作为干涉条纹的强度值;
存储单元,该存储单元存储作为由所述取样单元装入的所述强度值的干涉条纹强度值群;
运算单元,该运算单元对每一像素从所述存储单元存储的强度值群读出强度值,利用计算对象的像素的强度值和该计算对象的像素的附近的多个像素的强度值,假设计算对象的像素及该计算对象的像素附近的多个像素的干涉条纹波形的直流分量、交流振幅和相位相等,干涉条纹使计算对象的像素的强度值g(x,y)在该计算对象的像素附近,拟合干涉条纹波形表达式g(x,y)=a+bcos{2πfxx+2πfyy+φ},从而对每一单色光求出所述各像素的相位后,
根据每一单色光求出的各像素的相位求出多个表面高度候补群,从各候补群求出共同的高度作为实际高度,
进而,从该求出的所述测量对象面的表面高度,求出表面形状,这里,(x,y)是计算对象的像素的位置坐标,a是干涉条纹波形中包含的直流分量,b是干涉条纹波形中包含的交流振幅,f是干涉条纹的空间频率分量,φ是应计算的相位。
5.如权利要求4中所述的用多波长来测量表面形状的测量装置,其特征在于,
所述测量装置还包含下列组成要素:
驱动单元,该驱动单元中至少由所述照明单元、摄像单元、分支单元和基准面构成光学系统单元,并且
使所述保持单元和所述光学系统单元的至少一方移动,以便使朝向所述测量对象物的光和1个或多个测量对象物相对平行移动;
控制单元,该控制单元一面使光学系统单元和保持单元相对平行移动、一面控制所述各单元工作,以便求出测量对象物的表面形状。
6.如权利要求4或5中所述的用多波长来测量表面形状的测量装置,其特征在于,
照明单元具有输出不同波长的多个单色光源。
7.如权利要求6中所述的用多波长来测量表面形状的测量装置,其特征在于,
所述摄像单元配备将多个单色光分离的滤色片,
所述运算单元用所述滤色片特性产生的串扰的影响,去除每一所述单色光包含的其它单色光的干涉条纹的强度值。
8.如权利要求4或5中所述的用多波长来测量表面形状的测量装置,其特征在于,
照明单元的组成包含
白色光源、以及
光学单元,该光学单元从白色光源分离为规定的不同的多个波长的单色光,并使其朝向所述分支单元。
9.一种用多波长来测量表面形状的测量装置,通过分支单元对测量对象面和基准面照射单色光,根据因从测量对象面和基准面两者反射并返回同一光路的反射光而产生的干涉条纹的强度值,求出测量对象面的表面高度和表面形状,其特征在于,所述测量装置包含下列组成要素:
保持单元,该保持单元以对光的行进方向倾斜任意角度的姿势配置所述基准面,而且放置并保持测量对象物;
照明单元,该照明单元输出多个波长组成的光;
分离单元,该分离单元将照射所述光后从测量对象物和基准面反射并返回同一光路的反射光分离为不同波长的多个单色光;
摄像单元,该摄像单元使分离出的每一所述单色光产生干涉条纹,并拍摄测量对象面;
取样单元,该取样单元对每一像素装入所拍摄的所述测量对象面,作为干涉条纹的强度值;
存储单元,该存储单元存储作为由所述取样单元装入的所述强度值的干涉条纹强度值群;
运算单元,该运算单元对每一像素从所述存储单元存储的强度值群读出强度值,利用计算对象的像素的强度值和该计算对象的像素的附近的多个像素的强度值,假设计算对象的像素及该计算对象的像素附近的多个像素的干涉条纹波形的直流分量、交流振幅和相位相等,干涉条纹使计算对象的像素的强度值g(x,y)在该计算对象的像素附近,拟合干涉条纹波形表达式g(x,y)=a+bcos{2πfxx+2πfyy+φ},从而对每一单色光求出各像素的相位后,
从根据每一单色光求出的各像素的相位进行换算并求出的多个表面高度候补群,求出共同的高度作为实际高度,
进而,从该求出的所述测量对象面的表面高度,求出表面形状,这里,(x,y)是计算对象的像素的位置坐标,a是干涉条纹波形中包含的直流分量,b是干涉条纹波形中包含的交流振幅,f是干涉条纹的空间频率分量,φ是应计算的相位。
10.如权利要求9中所述的用多波长来测量表面形状的测量装置,其特征在于,
所述分离单元是将多个单色光分离的滤色片,并且
所述运算单元去除因所述滤色片特性产生的串扰的影响而在每一所述单色光中包含的其它单色光的干涉条纹的强度值。
CN2008100057818A 2007-01-31 2008-01-31 用多波长来测量表面形状测量方法和使用此方法的装置 Expired - Fee Related CN101236067B (zh)

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
JP2007021870 2007-01-31
JP2007-021870 2007-01-31
JP2007021870 2007-01-31
JP2008-008233 2008-01-17
JP2008008233 2008-01-17
JP2008008233A JP4885154B2 (ja) 2007-01-31 2008-01-17 複数波長による表面形状の測定方法およびこれを用いた装置

Publications (2)

Publication Number Publication Date
CN101236067A CN101236067A (zh) 2008-08-06
CN101236067B true CN101236067B (zh) 2011-02-16

Family

ID=39785804

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008100057818A Expired - Fee Related CN101236067B (zh) 2007-01-31 2008-01-31 用多波长来测量表面形状测量方法和使用此方法的装置

Country Status (4)

Country Link
JP (1) JP4885154B2 (zh)
KR (1) KR20080071905A (zh)
CN (1) CN101236067B (zh)
TW (1) TW200839177A (zh)

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5701159B2 (ja) * 2010-06-24 2015-04-15 東レエンジニアリング株式会社 干渉縞モデル適合による表面形状測定方法およびその装置
JP2012078269A (ja) * 2010-10-05 2012-04-19 Naoyuki Furuyama 測距方法及びレーザ測距装置
JP2013068489A (ja) * 2011-09-21 2013-04-18 Toray Eng Co Ltd 複数波長による表面形状測定方法およびこれを用いた装置
JP5954979B2 (ja) * 2011-12-15 2016-07-20 キヤノン株式会社 多波長干渉計を有する計測装置
WO2013088871A1 (ja) * 2011-12-16 2013-06-20 東レエンジニアリング株式会社 干渉色のモデル適合による膜厚測定方法およびその装置
US9147102B2 (en) * 2012-01-02 2015-09-29 Camtek Ltd. Method and system for measuring bumps based on phase and amplitude information
DE102012002174B4 (de) * 2012-02-07 2014-05-15 Schott Ag Vorrichtung und Verfahren zum Erkennen von Fehlstellen innerhalb des Volumens einer transparenten Scheibe und Verwendung der Vorrichtung
KR101341620B1 (ko) * 2012-06-14 2013-12-13 전자부품연구원 3차원 형상 측정 시스템 및 방법
JP5997578B2 (ja) * 2012-10-19 2016-09-28 東レエンジニアリング株式会社 クロストーク補正係数算出方法およびクロストーク補正係数算出機能を備えた透明膜の膜厚測定装置
CN103267494B (zh) * 2013-05-20 2015-11-04 湖北工业大学 一种表面形貌干涉测量的方法及装置
TWI467129B (zh) * 2014-01-07 2015-01-01 China Steel Corp 鑄嘴平坦度之檢測方法
JP6351289B2 (ja) * 2014-02-18 2018-07-04 Ntn株式会社 表面形状測定装置、方法およびプログラム
JP6126640B2 (ja) * 2015-05-11 2017-05-10 Ckd株式会社 三次元計測装置及び三次元計測方法
CN105021137B (zh) * 2015-06-30 2017-11-07 西安空间无线电技术研究所 一种快速转镜动态面形的测试系统
TWI625516B (zh) * 2015-09-15 2018-06-01 Ngk Insulators Ltd 複合基板及壓電基板的厚度傾向推定方法
TWI579525B (zh) * 2015-12-23 2017-04-21 國立臺灣大學 運動物件之絕對定位距離與偏擺角度同步量測之光學系統與方法
WO2017168469A1 (ja) 2016-03-28 2017-10-05 パナソニックIpマネジメント株式会社 外観検査装置と外観検査方法
KR101804527B1 (ko) * 2016-05-20 2017-12-05 주식회사 미르기술 다파장 광 주사 간섭계를 이용한 3차원 형상 측정장치
JP6246875B1 (ja) * 2016-08-24 2017-12-13 Ckd株式会社 計測装置
CN107966453B (zh) * 2016-10-20 2020-08-04 上海微电子装备(集团)股份有限公司 一种芯片缺陷检测装置及检测方法
CN108122797A (zh) * 2016-11-29 2018-06-05 上海微电子装备(集团)股份有限公司 一种3d检测装置
JP2019060683A (ja) * 2017-09-26 2019-04-18 東レエンジニアリング株式会社 塗布液滴の特性評価装置
CN110118533B (zh) * 2018-02-05 2021-08-03 上海微电子装备(集团)股份有限公司 一种三维检测方法及检测装置
CN108871206B (zh) * 2018-08-23 2021-06-22 业成科技(成都)有限公司 表面测量方法以及表面量测装置
CN109099859B (zh) * 2018-09-26 2021-07-27 中国科学院上海光学精密机械研究所 大口径光学元件表面缺陷三维形貌测量装置和方法
JP7257162B2 (ja) * 2019-02-08 2023-04-13 株式会社キーエンス 検査装置
US11825211B2 (en) * 2019-03-21 2023-11-21 Mit Semiconductor (Tian Jin) Co., Ltd Method of color inspection by using monochrome imaging with multiple wavelengths of light
JP7279596B2 (ja) * 2019-09-24 2023-05-23 株式会社デンソーウェーブ 三次元計測装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1673666A (zh) * 2005-04-15 2005-09-28 天津大学 基于相移干涉图像序列解析的微结构三维信息获取方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62282205A (ja) * 1986-05-31 1987-12-08 Toshiba Corp 物体の形状誤差を測定する方法およびその装置
JPH11218411A (ja) * 1998-02-02 1999-08-10 Fuji Xerox Co Ltd 干渉計測方法および干渉計測装置
JP4183089B2 (ja) * 2004-09-16 2008-11-19 東レエンジニアリング株式会社 表面形状および/または膜厚測定方法およびその装置
WO2007088789A1 (ja) * 2006-02-01 2007-08-09 Tokyo Institute Of Technology 表面形状の測定方法およびこれを用いた装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1673666A (zh) * 2005-04-15 2005-09-28 天津大学 基于相移干涉图像序列解析的微结构三维信息获取方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
闫聚群.大台阶高度测量的外差共焦方法.《中国激光》.2005,第32卷(第3期),389-393. *

Also Published As

Publication number Publication date
JP4885154B2 (ja) 2012-02-29
CN101236067A (zh) 2008-08-06
JP2008209404A (ja) 2008-09-11
TW200839177A (en) 2008-10-01
KR20080071905A (ko) 2008-08-05

Similar Documents

Publication Publication Date Title
CN101236067B (zh) 用多波长来测量表面形状测量方法和使用此方法的装置
US7852489B2 (en) Method for measuring surface profile, and apparatus using the same
US20090213386A1 (en) Apparatus and method for measuring surface topography of an object
CN106500589A (zh) 一种多波长可调谐显微干涉的测量方法及其装置
TWI567364B (zh) 結構光產生裝置、量測系統及其方法
US20080117438A1 (en) System and method for object inspection using relief determination
CN103727891A (zh) 同步三维散斑干涉测量系统及测量方法
US10001405B2 (en) Measurement device for obtaining amplitude information of an object
US10495438B2 (en) Measurement device
US7561279B2 (en) Scanning simultaneous phase-shifting interferometer
US6734978B2 (en) Profile measuring method and measurement apparatus using interference of light
KR100916593B1 (ko) 실시간 3차원 형상 측정 시스템
TW200408793A (en) Object surface 3-D contour measurement method and system
CN106323163A (zh) 一种表面3d检测装置及检测方法
JP5701159B2 (ja) 干渉縞モデル適合による表面形状測定方法およびその装置
Kitagawa Multiwavelength single-shot interferometry without carrier fringe introduction
US11248899B2 (en) Method and apparatus for deriving a topography of an object surface
CN102853761A (zh) 一种空间相移器
JP2010185844A (ja) 表面形状測定方法およびこれを用いた装置
KR101928580B1 (ko) 위상 및 진폭 정보에 기초한 범프 측정방법 및 시스템
JP2007071817A (ja) 二光束干渉計及び同干渉計を用いた被測定物の形状測定方法
JPH0587541A (ja) 2次元情報測定装置
CN107036529A (zh) 偏振同步相移干涉仪
Pribanic et al. Scene Adaptive Structured Light 3D Imaging.
TWI802988B (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

Granted publication date: 20110216

Termination date: 20180131

CF01 Termination of patent right due to non-payment of annual fee