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

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

Info

Publication number
CN104279983A
CN104279983A CN201410616949.4A CN201410616949A CN104279983A CN 104279983 A CN104279983 A CN 104279983A CN 201410616949 A CN201410616949 A CN 201410616949A CN 104279983 A CN104279983 A CN 104279983A
Authority
CN
China
Prior art keywords
int
film thickness
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
CN201410616949.4A
Other languages
English (en)
Other versions
CN104279983B (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 CN201410616949.4A priority Critical patent/CN104279983B/zh
Publication of CN104279983A publication Critical patent/CN104279983A/zh
Application granted granted Critical
Publication of CN104279983B publication Critical patent/CN104279983B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

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

Description

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

Claims (2)

1.基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法,其特征在于,包括以下步骤:
步骤a、获得镀膜待测样品三维形貌,轴向扫描间隔为z,扫描层数为N;
步骤b、选取三维形貌梯度最大点,将该点轴向响应曲线数据归一化,记为I1,将其余点轴向响应曲线数据归一化,记为I2
步骤c、以sinc2(a(x×z-b))为目标函数,利用Levenberg-Marquardt算法分别拟合I1与I2,求得两组a和b,第一组对应I1的最优拟合参数,记为a1和b1,第二组对应I2的最优拟合参数,记为a2和b2
步骤d、分别计算I1(n)=sinc2(a1×n×z)和I2(n)=sinc2(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()表示取整运算。
CN201410616949.4A 2014-11-05 2014-11-05 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法 Active CN104279983B (zh)

Priority Applications (1)

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

Applications Claiming Priority (1)

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

Publications (2)

Publication Number Publication Date
CN104279983A true CN104279983A (zh) 2015-01-14
CN104279983B CN104279983B (zh) 2016-10-05

Family

ID=52255112

Family Applications (1)

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

Country Status (1)

Country Link
CN (1) CN104279983B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1851522A (zh) * 2006-05-26 2006-10-25 清华大学 激光扫描共聚焦显微镜扫描畸变现象的全场校正方法
CN102207378A (zh) * 2011-03-09 2011-10-05 浙江大学 球面干涉检测中基于波前差分的高精度调整误差校正方法
CN102654638A (zh) * 2010-12-29 2012-09-05 徕卡显微系统复合显微镜有限公司 用于校正共聚焦扫描显微镜中的图像失真的方法
CN102853780A (zh) * 2012-09-07 2013-01-02 厦门大学 抛光工件表面轮廓各频段误差的分离方法
CN102866492A (zh) * 2011-07-06 2013-01-09 黄书伟 平衡检测共聚焦显微镜成像系统及其三维影像重建方法
JP5597205B2 (ja) * 2009-11-19 2014-10-01 キヤノン株式会社 被検面の形状を計測する装置及び被検面の形状を算出するためのプログラム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1851522A (zh) * 2006-05-26 2006-10-25 清华大学 激光扫描共聚焦显微镜扫描畸变现象的全场校正方法
JP5597205B2 (ja) * 2009-11-19 2014-10-01 キヤノン株式会社 被検面の形状を計測する装置及び被検面の形状を算出するためのプログラム
CN102654638A (zh) * 2010-12-29 2012-09-05 徕卡显微系统复合显微镜有限公司 用于校正共聚焦扫描显微镜中的图像失真的方法
CN102207378A (zh) * 2011-03-09 2011-10-05 浙江大学 球面干涉检测中基于波前差分的高精度调整误差校正方法
CN102866492A (zh) * 2011-07-06 2013-01-09 黄书伟 平衡检测共聚焦显微镜成像系统及其三维影像重建方法
CN102853780A (zh) * 2012-09-07 2013-01-02 厦门大学 抛光工件表面轮廓各频段误差的分离方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HAIHUA CUI ECT.ALL: "A three-step system calibration procedure with error compensation fo 3D shape measurement", 《CHINESE OPTICS LETTERS》 *
谭佐军等: "自适应光学高分辨率共聚焦显微成像技术", 《激光与光电子学进展》 *

Also Published As

Publication number Publication date
CN104279983B (zh) 2016-10-05

Similar Documents

Publication Publication Date Title
Enoch et al. Comparing star formation on large scales in the c2d legacy clouds: Bolocam 1.1 mm dust continuum surveys of serpens, perseus, and ophiuchus
CN106407502B (zh) 一种基于最佳匹配的叶片截面型线轮廓参数评价方法
CN102800096B (zh) 一种摄像机参数的鲁棒性估计算法
CN110095066B (zh) 基于Mean-shift的光谱共焦信号峰值波长快速高精度提取方法
CA3155820A1 (en) ACCOUNTING FOR ERRORS IN OPTICAL MEASUREMENTS
CN105931257B (zh) 基于纹理特征和结构相似度的sar图像质量评估方法
CN111551117B (zh) 显微图像焦点漂移距离的测量方法及系统、计算机设备
CN102607455B (zh) 基于光学显微镜和变化光照的微观形貌三维测量方法
KR101628557B1 (ko) 시편 표면의 마찰계수 측정방법
CN103949771B (zh) 一种基于特征人工节瘤缺陷的激光预处理技术
Zhang et al. Accuracy improvement in laser stripe extraction for large-scale triangulation scanning measurement system
CN104408709A (zh) 一种线阵ccd相机图像灰度不匀快速矫正方法
CN108844494A (zh) 基于光学方法检测眼镜曲面曲率参数的方法
CN104296684A (zh) 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法
CN111415349A (zh) 一种基于图像处理技术的涤纶长丝毛丝检测方法
Liu et al. Compact chromatic confocal sensor for displacement and thickness measurements
CN104279983A (zh) 基于表面镀膜共焦显微形貌测量装置的膜厚误差校正方法
CN104166985B (zh) 一种基于区域分割的星敏感器标定方法
CN104406537A (zh) 测量部件凹坑微米级深度的方法
CN103837102A (zh) 一种测量圆形光学平面表面粗糙度的取样方法
CN110346782A (zh) 一种长距离地面三维激光雷达回波强度数据的改正方法
CN112414316B (zh) 一种应变片敏感栅尺寸参数测量方法
CN106707474B (zh) 一种物像双侧远心光学系统
CN207074147U (zh) 一种利用线阵ccd成像法测量杨氏弹性模量的实验系统
CN114663781A (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