CN104296684A - 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法 - Google Patents

基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法 Download PDF

Info

Publication number
CN104296684A
CN104296684A CN201410616961.5A CN201410616961A CN104296684A CN 104296684 A CN104296684 A CN 104296684A CN 201410616961 A CN201410616961 A CN 201410616961A CN 104296684 A CN104296684 A CN 104296684A
Authority
CN
China
Prior art keywords
film thickness
int
designated
data
surface coating
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
Application number
CN201410616961.5A
Other languages
English (en)
Other versions
CN104296684B (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201410616961.5A priority Critical patent/CN104296684B/zh
Publication of CN104296684A publication Critical patent/CN104296684A/zh
Application granted granted Critical
Publication of CN104296684B publication Critical patent/CN104296684B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法属于共焦扫描光学测量技术领域;该方法在获得镀膜待测样品三维形貌的基础上,计算梯度最大点和其余点的轴向响应曲线数据归一化结果,并以sinc4(a(x-b))为目标函数进行拟合,将梯度最大点的拟合结果与不同宽度矩形函数做卷积运算,再与待校准点的轴向响应数据做差运算,利用最小残差所对应的矩形窗宽度来补偿膜厚误差;本发明基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法,通过拟合薄膜、厚膜轴向响应曲线,实现对荧光膜膜厚引入误差的补偿,有效校正镀膜膜厚不均引起的误差,并将此误差降低到十分之一膜厚以下。

Description

基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法
技术领域
基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法属于共焦扫描光学测量技术领域。
背景技术
对于光滑大区率表面样品,由于照明光在其表面发生镜面反射,当测量与样品表面法线方向与测量系统光轴超过一定角度,经样品表面反射回的信号光无法被口径有限的物镜全部收集,甚至无法收集,从而造成此类样品无法测量或测量分辨率降低。解决上述问题,可以在光滑大区率表面镀上一层荧光膜,使得光滑样品表面具有散射特性,在理想情况下,信号光以近90度的立体角散射传播,实现信号光的全口径收集。然而,由于镀膜工艺无法保证样品表面的荧光膜厚度一致,荧光膜厚度会随着样品表面斜率的增加而变薄,这种膜厚不均会引入二分之一厚膜厚度的误差,影响测量精度。
发明内容
为了解决上述问题,本发明设计了一种基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法,该膜厚误差校正方法能够实现对荧光膜膜厚引入误差的补偿,有效校正镀膜膜厚不均引起的误差,并将此误差降低到十分之一膜厚以下。
本发明的目的是这样实现的:
基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法,包括以下步骤:
步骤a、获得镀膜待测样品三维形貌,轴向扫描间隔为z,扫描层数为N;
步骤b、选取三维形貌梯度最大点,将该点轴向响应曲线数据归一化,记为I1,将其余点轴向响应曲线数据归一化,记为I2
步骤c、以sinc4(a(x×z-b))为目标函数,利用Levenberg-Marquardt算法分别拟合I1与I2,求得两组a和b,第一组对应I1的最优拟合参数,记为a1和b1,第二组对应I2的最优拟合参数,记为a2和b2
步骤d、分别计算I1(n)=sinc4(a1×n×z)和I2(n)=sinc4(a2×n×z),n为大于min(int(-π/(a1×z)),int(-π/(a2×z))),小于max(int(-π/(a1×z)),int(-π/(a2×z)))的整数,min()表示最小值,max()表示最大值,int()表示取整运算;
步骤e、将I1(n)与矩形函数做卷积运算,结果记为I'(n,hi),有: hi为非0整数,也是矩形函数的偏移量,还是矩形函数的半宽度;
步骤f、从I'(n,hi)与I2(n)中提取半高以上数据,分别为I3与I4
步骤g、计算I3与I4的残差res(i),其中,与res(i)最小值对应的hi记为hmin,该点高度校正值为hmin×z;
步骤h、在待校准点测得高度值上减去hmin×z。
上述基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法,所述的步骤f中,提取半高以上数据的具体实施方式为:
对于I'(n,hi),以I'(n,hi)最大值为中心,左右各取int(N/4)个点数据,这些数据记为I3,其中,int()表示取整运算;
对于I2(n),以I2(n)最大值为中心,左右各取int(N/4)个点数据,这些数据记为I4,其中,int()表示取整运算。
有益效果:本发明通过拟合薄膜、厚膜轴向响应曲线,实现对荧光膜膜厚引入误差的补偿,有效校正镀膜膜厚不均引起的误差,并将此误差降低到十分之一膜厚以下。
附图说明
图1是本发明膜厚误差校正方法流程图。
图2是镀膜待测样品三维形貌图。
图3是镀膜待测样品三维形貌梯度最大点的剖面图。
图4是与I1对应的拟合结果。
图5是与I2对应的拟合结果。
具体实施方式
下面结合附图和具体实施例进行详细说明,以便对本发明的目的、技术方案有更深入的理解。
具体实施例步骤说明如下:
步骤a、获得镀膜待测样品三维形貌,具体步骤可参照专利《基于荧光共焦显微技术的光滑大区率样品测量装置与方法》,轴向扫描间隔z=50nm,扫描层数N=141,镀膜待测样品三维形貌如图2;
步骤b、选取三维形貌梯度最大点,即剖面图图3中点A,将A点轴向响应曲线数据归一化,记为I1,将其余点轴向响应曲线数据归一化,不失一般性,记其余任一点的数据为I2
步骤c、以sinc4(a(x×z-b))为目标函数,利用Levenberg-Marquardt算法分别拟合I1与I2,求得两组a和b,第一组对应I1的最优拟合参数,记为a1和b1,第二组对应I2的最优拟合参数,记为a2和b2,其中a1、b1、a2和b2分别为0.35、4.30、0.32和4.50;与I1对应的拟合结果如图4所示,与I2对应的拟合结果如图5所示;
步骤d、分别计算I1(n)=sinc4(a1×n×z)和I2(n)=sinc4(a2×n×z),n为大于min(int(-π/(a1×z)),int(-π/(a2×z))),小于max(int(-π/(a1×z)),int(-π/(a2×z)))的整数,min()表示最小值,max()表示最大值,int()表示取整运算,I1(n)与I2(n)曲线分别对应图4与图5中“拟合sinc4曲线”;
步骤e、将I1(n)与矩形函数做卷积运算rect(n-hi),结果记为I'(n,hi),有: hi为非0整数,也是矩形函数的偏移量,还是矩形函数的半宽度;
步骤f、从I'(n,hi)与I2(n)中提取半高以上数据,分别为I3与I4
步骤g、计算I3与I4的残差res(i),其中,与res(i)最小值对应的hi为hmin=2,该点高度校正值为hmin×z=100nm;
步骤h、在待校准点测得高度值上减去100nm。
按照上述步骤将A以外的任意一点高度进行校正,最终可得校正后的三维形貌。
这里需要说明的是,同sinc2(a(x-b))相比,采用sinc4(a(x-b))函数为目标函数进行拟合,对于表面平滑的样品,测量结果对于样品三维形貌的表达更加准确。
本发明不局限于上述最佳实施方式,任何人应该得知在本发明的启示下作出的结构变化或方法改进,凡是与本发明具有相同或相近的技术方案,均落入本发明的保护范围之内。

Claims (2)

1.基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法,其特征在于,包括以下步骤:
步骤a、获得镀膜待测样品三维形貌,轴向扫描间隔为z,扫描层数为N;
步骤b、选取三维形貌梯度最大点,将该点轴向响应曲线数据归一化,记为I1,将其余点轴向响应曲线数据归一化,记为I2
步骤c、以sinc4(a(x×z-b))为目标函数,利用Levenberg-Marquardt算法分别拟合I1与I2,求得两组a和b,第一组对应I1的最优拟合参数,记为a1和b1,第二组对应I2的最优拟合参数,记为a2和b2
步骤d、分别计算I1(n)=sinc4(a1×n×z)和I2(n)=sinc4(a2×n×z),n为大于min(int(-π/(a1×z)),int(-π/(a2×z))),小于max(int(-π/(a1×z)),int(-π/(a2×z)))的整数,min()表示最小值,max()表示最大值,int()表示取整运算;
步骤e、将I1(n)与矩形函数做卷积运算,结果记为I'(n,hi),有: hi为非0整数,也是矩形函数的偏移量,还是矩形函数的半宽度;
步骤f、从I'(n,hi)与I2(n)中提取半高以上数据,分别为I3与I4
步骤g、计算I3与I4的残差res(i),其中,与res(i)最小值对应的hi记为hmin,该点高度校正值为hmin×z;
步骤h、在待校准点测得高度值上减去hmin×z。
2.根据权利要求1所述的基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法,其特征在于,所述的步骤f中,提取半高以上数据的具体实施方式为:
对于I'(n,hi),以I'(n,hi)最大值为中心,左右各取int(N/4)个点数据,这些数据记为I3,其中,int()表示取整运算;
对于I2(n),以I2(n)最大值为中心,左右各取int(N/4)个点数据,这些数据记为I4,其中,int()表示取整运算。
CN201410616961.5A 2014-11-05 2014-11-05 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法 Active CN104296684B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410616961.5A CN104296684B (zh) 2014-11-05 2014-11-05 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410616961.5A CN104296684B (zh) 2014-11-05 2014-11-05 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法

Publications (2)

Publication Number Publication Date
CN104296684A true CN104296684A (zh) 2015-01-21
CN104296684B CN104296684B (zh) 2016-11-30

Family

ID=52316520

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410616961.5A Active CN104296684B (zh) 2014-11-05 2014-11-05 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法

Country Status (1)

Country Link
CN (1) CN104296684B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109003229A (zh) * 2018-08-09 2018-12-14 成都大学 基于三维增强深度残差网络的磁共振超分辨率重建方法
CN113029032A (zh) * 2021-03-26 2021-06-25 中南大学 基于光谱共焦的高精度面形测量方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001091844A (ja) * 1999-09-22 2001-04-06 Olympus Optical Co Ltd 共焦点走査型顕微鏡
CN1301357A (zh) * 1998-03-16 2001-06-27 普雷勒克斯公司 共焦显微镜成像系统
TW200837385A (en) * 2007-03-14 2008-09-16 Univ Nat Yang Ming Modulation differential confocal microscopy
CN103090787A (zh) * 2013-01-29 2013-05-08 哈尔滨工业大学 基于被测表面荧光激发的共焦显微测量装置
CN103411558A (zh) * 2013-08-15 2013-11-27 哈尔滨工业大学 一种角谱扫描照明阵列式共焦微结构测量装置与方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1301357A (zh) * 1998-03-16 2001-06-27 普雷勒克斯公司 共焦显微镜成像系统
JP2001091844A (ja) * 1999-09-22 2001-04-06 Olympus Optical Co Ltd 共焦点走査型顕微鏡
TW200837385A (en) * 2007-03-14 2008-09-16 Univ Nat Yang Ming Modulation differential confocal microscopy
CN103090787A (zh) * 2013-01-29 2013-05-08 哈尔滨工业大学 基于被测表面荧光激发的共焦显微测量装置
CN103411558A (zh) * 2013-08-15 2013-11-27 哈尔滨工业大学 一种角谱扫描照明阵列式共焦微结构测量装置与方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LP SUNG ECT.ALL: "《Use of Laser Scanning Confocal Microscopy for Characterizing Changes in Film Thickness and Local Surface Morphology of UV-Exposed Polymer Coatings》", 《JOURNAL OF COATINGS TECHNOLOGY & RESEARCH》 *
张平 等: "《荧光共焦扫描系统成像特性的优化》", 《光学学报》 *
袁成清 等: "《基于激光共焦扫描显微镜方法的磨损表面三维数字化描述》", 《润滑与密封》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109003229A (zh) * 2018-08-09 2018-12-14 成都大学 基于三维增强深度残差网络的磁共振超分辨率重建方法
CN109003229B (zh) * 2018-08-09 2022-12-13 成都大学 基于三维增强深度残差网络的磁共振超分辨率重建方法
CN113029032A (zh) * 2021-03-26 2021-06-25 中南大学 基于光谱共焦的高精度面形测量方法及装置

Also Published As

Publication number Publication date
CN104296684B (zh) 2016-11-30

Similar Documents

Publication Publication Date Title
CN104315978B (zh) 一种管路端面中心点的测量方法及装置
CN102800096B (zh) 一种摄像机参数的鲁棒性估计算法
CN103575227A (zh) 一种基于数字散斑的视觉引伸计实现方法
CN105509653A (zh) 机械零件轮廓度视觉测量方法及系统
CN105092607A (zh) 球面光学元件表面缺陷评价方法
CN103292701A (zh) 基于机器视觉的精密器件在线尺寸测量方法
CN103714329A (zh) 一种仪表指针识别的检测算法
CN109017867A (zh) 钢轨波磨动态测量方法
CN106097430B (zh) 一种多高斯信号拟合的激光条纹中心线提取方法
CN108844494A (zh) 基于光学方法检测眼镜曲面曲率参数的方法
CN112881042B (zh) 基于激光轮廓的轮对检测系统、方法和介质
CN104408709A (zh) 一种线阵ccd相机图像灰度不匀快速矫正方法
CN105180846A (zh) 一种陶瓷插芯同轴度检测方法
CN105336635B (zh) Cd-sem装置的校正方法、应用cd-sem装置的方法及cd-sem装置
CN104296684A (zh) 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法
CN111415349A (zh) 一种基于图像处理技术的涤纶长丝毛丝检测方法
CN110763136A (zh) 一种高精度三维工件尺寸检测系统及方法
CN109359604A (zh) 面向巡检机器人的阴影干扰下仪表识别方法
CN108827915B (zh) 一种基于光电传感阵列测量折射率的亚像素位置获取方法
CN108225212B (zh) 曲面玻璃轮廓度检测方法
CN104361226A (zh) 共焦轴向响应曲线峰值位置提取算法
CN111815580B (zh) 一种图像边缘识别方法及小模数齿轮模数检测方法
CN113888641A (zh) 一种基于机器视觉和深度学习的立木胸径测量方法
CN104166985B (zh) 一种基于区域分割的星敏感器标定方法
CN104279983A (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