CN113188456A - 基于多重修正的荧光染色薄膜厚度测量标定方法 - Google Patents
基于多重修正的荧光染色薄膜厚度测量标定方法 Download PDFInfo
- Publication number
- CN113188456A CN113188456A CN202110429859.4A CN202110429859A CN113188456A CN 113188456 A CN113188456 A CN 113188456A CN 202110429859 A CN202110429859 A CN 202110429859A CN 113188456 A CN113188456 A CN 113188456A
- Authority
- CN
- China
- Prior art keywords
- film thickness
- pixel units
- fluorescence
- jth column
- ith row
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
- G01B11/0616—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material of coating
Abstract
本发明以若干不同标准厚度薄膜的荧光照片作为样本,基于多项式并考虑荧光成像区域的光学特性,建立薄膜厚度与荧光强度之间的关系式,构建标定系数方程组并求解,使用统计学方法减小随机误差,得到标定系数表。拍摄待测薄膜的荧光照片,基于标定系数表以及薄膜厚度与荧光强度之间的关系式,计算得到每一像素单元对应的薄膜厚度。通过再次测量标准厚度薄膜,得到测量误差分布情况,对标定系数与测量值进行校正,可显著减小由非线性、荧光成像区域光学特性以及荧光强度随机性带来的测量误差,能有效提高荧光染色薄膜厚度的测量精度以及量程。
Description
技术领域
本发明涉及测量标定技术领域,具体涉及一种基于多重修正的荧光染色薄膜厚度测量标定方法。
背景技术
荧光染色薄膜厚度测量技术是一种通过激发诱导荧光成像,测量透明荧光染色薄膜厚度的技术,可用于测量液体和固体薄膜的厚度,如测量结构间隙的高度、表面凹槽的深度以及LED荧光薄膜厚度等,在机械、电子、生物等行业中均有应用。
为定量测量薄膜厚度,需要通过标定得到薄膜厚度与荧光强度之间的关系式,确定关系式中各项标定系数,进而可由荧光强度计算出薄膜厚度。然而,在实际的测量过程中,荧光强度与薄膜厚度之间的非线性关系、荧光成像区域的光学特性、以及荧光强度的随机性等均会造成测量误差,使常规荧光厚度测量技术的量程受到限制,精度难以保证。
发明内容
有鉴于此,本发明提供了一种基于多重修正的荧光染色薄膜厚度测量标定方法,能够显著减小由非线性、荧光成像区域光学特性以及荧光强度随机性带来的误差,能有效提高荧光染色薄膜厚度的测量精度以及量程。
本发明的一种基于多重修正的荧光染色薄膜厚度测量标定方法,以若干不同标准厚度薄膜的荧光照片作为样本,基于多项式并考虑荧光成像区域的光学特性,建立薄膜厚度与荧光强度之间的关系式,构建标定系数方程组并求解,使用统计学方法减小随机误差,得到标定系数表;
拍摄待测薄膜的荧光照片,基于标定系数表以及薄膜厚度与荧光强度之间的关系式,计算得到每一像素单元对应的薄膜厚度;
通过再次测量标准厚度薄膜,得到测量误差分布情况,对标定系数与测量值进行校正;
使用校正过的计算公式,重新计算薄膜厚度作为最终厚度。
其中,所建立的薄膜厚度与荧光强度之间的关系式为:
其中,i、j分别为像素单元所在行和列的编号;n为基函数阶次的编号, n=0,1,2…N;N为基函数的总阶数;k为图片编号,k=1,2,3…K;K为图片总数;为第k张图片的第i行、第j列像素单元对应的薄膜厚度;为第i行、第j 列像素单元对应的第n阶标定系数;为第k张图片的第i行、第j列像素单元对应的第n阶基函数,其表达式为:
其中,Np为考虑使用的多项式的阶数;为第k张图片的第i行、第j列像素单元对应的荧光强度;为第k张图片的第i行、第j列像素单元对应的第(n-Np+1)阶荧光成像区域对应的像素单元的荧光强度之和,其表达式为:
其中,r为荧光成像区域阶次的编号。
其中,使用分组平均拟合法时,构建标定系数线性方程组为:
其中,m为图片分组编号,m=1,2,3…M;M为图片总组数;为第m组第k张图片的第i行、第j列像素单元对应的第n阶基函数项;为第m组图片的第i行、第j列像素单元对应的第n阶标定系数;为第m组第k张图片的第i行、第j列像素单元对应的薄膜厚度;此时要求K=N+1;
求解标定系数线性方程组,对各个样本组得到的标定系数计算均值如下:
其中,使用最小二乘法得到标定系数时,对每张图片的某一像素单元,构建标定系数线性方程组为:
其中,要求K>N。
其中,基于标定系数表以及薄膜厚度与荧光强度之间的关系式,计算得到每一像素单元对应的薄膜厚度为:
其中,使用校正过的计算公式,计算薄膜厚度如下:
其中,通过再次测量标准厚度薄膜,得到测量误差分布情况,对标定系数与测量值进行校正,过程具体如下:
步骤S201:拍摄若干不同标准厚度薄膜的荧光照片,总数为S;
步骤S202:使用标定系数表,根据薄膜厚度与荧光强度之间的关系式,计算每一像素单元对应薄膜厚度测量值:
步骤S203:计算绝对误差与相对误差值如下:
步骤S204:计算使用当前标定系数表进行测量的平均偏差和标准偏差:
步骤S205:根据误差分布情况,引入并调整校正系数与校正函数参数,对标定系数和测量值进行校正;
步骤S206:重复执行步骤S201-步骤S204,直到误差指标满足要求;
步骤S207:校正后,实际应用计算薄膜厚度值如下:
其中,用于采用荧光成像系统测量经荧光染色的液体或固体薄膜厚度时的标定。
其中,拍摄若干不同标准厚度薄膜的荧光照片时,对焦在薄膜基准面。
有益效果:
本发明以若干不同标准厚度薄膜的荧光照片作为样本,基于多项式并考虑荧光成像区域的光学特性,建立薄膜厚度与荧光强度之间的关系式,构建标定系数方程组并求解,使用统计学方法减小随机误差,得到标定系数表。拍摄待测薄膜的荧光照片,基于标定系数表以及薄膜厚度与荧光强度之间的关系式,计算得到每一像素单元对应的薄膜厚度。通过再次测量标准厚度薄膜,得到测量误差分布情况,对标定系数与测量值进行校正,可显著减小由非线性、荧光成像区域光学特性以及荧光强度随机性带来的测量误差,能有效提高荧光染色薄膜厚度的测量精度以及量程。
附图说明
图1为本发明基于多重修正、使用分组平均法的荧光染色薄膜厚度测量标定流程图;
图2为本发明基于多重修正、使用最小二乘法的染色薄膜厚度测量标定流程图;
图3为本发明荧光成像区域与像素单元对应示意图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种适用于采用荧光成像系统测量经荧光染色的液体或固体薄膜厚度时的标定。综合考虑荧光强度非线性、荧光成像区域的光学特性和随机性,使用多重修正与误差校正方法,以若干不同标准厚度薄膜的荧光照片作为样本,基于多项式并考虑荧光成像区域光学特性,建立薄膜厚度与荧光强度之间的关系式,构建标定系数方程组并求解,使用统计学方法减小随机误差,得到标定系数表;对于待测薄膜的荧光照片,基于标定系数表以及薄膜厚度与荧光强度之间的关系式,计算得到每一像素单元对应的薄膜厚度;通过再次测量标准厚度薄膜,可得到测量精度与误差分布情况,并对标定系数进行校正,以此实现较高精度与较大量程的薄膜厚度测量。
本发明的具体实施步骤如下:
选定基函数空间为:
其表达式为:
对某一像素单元的某一组样本,建立薄膜厚度与荧光强度之间的关系式:
使用分组平均拟合法时,构建标定系数线性方程组:
求解标定系数线性方程组,对各个样本组得到的标定系数计算均值如下:
得到标定系数数组以及标定系数表;
使用最小二乘法得到标定系数时,对每张图片的某一像素单元,构建标定系数线性方程组:
求解得到标定系数表。
为分析误差并校正,拍摄若干不同标准厚度薄膜的荧光照片,根据薄膜厚度与荧光强度之间的关系式,计算每一像素单元对应薄膜厚度测量值:
计算绝对误差与相对误差值如下:
计算使用当前标定系数表进行测量的平均偏差和标准偏差:
根据误差分布情况,引入并调整校正系数与校正函数参数,对标定系数和测量值进行校正。
使用校正过的计算公式,计算薄膜厚度如下:
上述公式中,i、j分别为像素单元所在行和列的编号;n为基函数阶次的编号,n=0,1,2…N;N为考虑使用的基函数的总阶数;r为荧光成像区域阶次的编号;NP为考虑使用的多项式阶数;k为用于标定的图片编号;K为用于标定的图片总数;M为图片总组数;s为用于校正的图片编号;S为用于标定的图片总数;为第k张图片的第i行、第j列像素单元对应的薄膜厚度;为第k张图片的第i 行、第j列像素单元对应的荧光光强;为第m组第k张图片的第i行、第j 列像素单元对应的第n阶基函数项;为第m组图片的第i行、第j列像素单元对应的第n阶标定系数;为第i行、第j列像素单元对应的第n阶标定系数;为第i行、第j列像素单元对应的测量误差;为第i行、第j列像素单元对应的相对测量误差;为第i行、第j列像素单元对应的测量平均偏差;为第 i行、第j列像素单元对应的测量标准偏差;ht为公式校正厚度阈值;href为公式校正厚度参考值;为第i行、第j列像素单元对应的待校正的厚度测量值;为第s张图片的第i行、第j列像素单元对应的第n阶基函数;为第i行、第j 列像素单元对应的第n阶校正系数;为第i行、第j列像素单元对应的校正后的最终厚度测量值;为待测图片的第i行、第j列像素单元对应的待校正的薄膜厚度测量值;为待测图片的第i行、第j列像素单元对应的校正后的最终的薄膜厚度测量值。
本发明提供了基于多重修正、使用分组平均法或使用最小二乘法的荧光染色薄膜厚度测量标定方法完成标定。
(1)使用分组平均拟合法时,步骤如图1所示,具体如下:
S101:拍摄不同标准厚度薄膜的荧光照片(注意要严格对焦在薄膜基准面),每组K张,要求K=N+1,总组数为M;
S102:选定基函数空间如下:
其表达式为:
荧光成像区域阶次与像素单元对应关系如图3所示;
S103:对某一像素单元的某一组样本,建立薄膜厚度与荧光强度之间的关系式:
据此构建标定系数线性方程组:
S104:对每一组样本求解S3得到的标定系数线性方程组,对各个样本组得到的标定系数计算均值如下:
得到标定系数数组;
S105:对所有像素执行步骤S102-S104,得到标定系数表;
(2)当使用最小二乘法得到标定系数时,步骤如图2所示,具体如下:
S106:拍摄不同标准厚度薄膜的荧光照片(注意要严格对焦在薄膜基准面),总数为K,要求K>N;
S107:选定基函数空间如下:
其表达式为:
荧光成像区域阶次与像素单元对应关系如图3所示;
S108:对每张图片的某一像素单元,构建标定系数线性方程组:
S109:对所有像素单元执行步骤S107,得到标定系数表;
完成标定后,测试与误差分析校正过程具体如下:
S201:拍摄若干不同标准厚度薄膜的荧光照片(注意要严格对焦在薄膜基准面),总数为S;
S202:使用S1得到的标定系数表,根据薄膜厚度与荧光强度之间的关系式,计算每一像素单元对应薄膜厚度测量值:
S203:计算绝对误差与相对误差值如下:
S204:计算使用当前标定系数表进行测量的平均偏差和标准偏差:
S205:根据误差分布情况,引入并调整校正系数与校正函数参数,对标定系数和测量值进行校正;
S206:多次重复S201-S204,直到误差指标满足要求;
校正后,实际应用测量过程如下:
S301:使用与标定过程相同的荧光成像系统与成像条件,对待测薄膜进行拍摄;
S302:使用校正过的计算公式,计算薄膜厚度如下:
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于多重修正的荧光染色薄膜厚度测量标定方法,其特征在于,以若干不同标准厚度薄膜的荧光照片作为样本,基于多项式并考虑荧光成像区域的光学特性,建立薄膜厚度与荧光强度之间的关系式,构建标定系数方程组并求解,使用统计学方法减小随机误差,得到标定系数表;
拍摄待测薄膜的荧光照片,基于标定系数表以及薄膜厚度与荧光强度之间的关系式,计算得到每一像素单元对应的薄膜厚度;
通过再次测量标准厚度薄膜,得到测量误差分布情况,对标定系数与测量值进行校正;
使用校正过的计算公式,重新计算薄膜厚度作为最终厚度。
2.如权利要求1所述的基于多重修正的荧光染色薄膜厚度测量标定方法,其特征在于,所建立的薄膜厚度与荧光强度之间的关系式为:
其中,i、j分别为像素单元所在行和列的编号;n为基函数阶次的编号,n=0,1,2…N;N为基函数的总阶数;k为图片编号,k=1,2,3…K;K为图片总数;为第k张图片的第i行、第j列像素单元对应的薄膜厚度;为第i行、第j列像素单元对应的第n阶标定系数;为第k张图片的第i行、第j列像素单元对应的第n阶基函数,其表达式为:
其中,Np为考虑使用的多项式的阶数;为第k张图片的第i行、第j列像素单元对应的荧光强度;为第k张图片的第i行、第j列像素单元对应的第(n-Np+1)阶荧光成像区域对应的像素单元的荧光强度之和,其表达式为:
其中,r为荧光成像区域阶次的编号。
7.如权利要求6任意一项所述的基于多重修正的荧光染色薄膜厚度测量标定方法,其特征在于,通过再次测量标准厚度薄膜,得到测量误差分布情况,对标定系数与测量值进行校正,过程具体如下:
步骤S201:拍摄若干不同标准厚度薄膜的荧光照片,总数为S;
步骤S202:使用标定系数表,根据薄膜厚度与荧光强度之间的关系式,计算每一像素单元对应薄膜厚度测量值:
步骤S203:计算绝对误差与相对误差值如下:
步骤S204:计算使用当前标定系数表进行测量的平均偏差和标准偏差:
步骤S205:根据误差分布情况,引入并调整校正系数与校正函数参数,对标定系数和测量值进行校正;
步骤S206:重复执行步骤S201-步骤S204,直到误差指标满足要求;
步骤S207:校正后,实际应用计算薄膜厚度值如下:
8.如权利要求1-7任意一项所述的基于多重修正的荧光染色薄膜厚度测量标定方法,其特征在于,用于采用荧光成像系统测量经荧光染色的液体或固体薄膜厚度时的标定。
9.如权利要求6任意一项所述的基于多重修正的荧光染色薄膜厚度测量标定方法,其特征在于,拍摄若干不同标准厚度薄膜的荧光照片时,对焦在薄膜基准面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110429859.4A CN113188456B (zh) | 2021-04-21 | 2021-04-21 | 基于多重修正的荧光染色薄膜厚度测量标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110429859.4A CN113188456B (zh) | 2021-04-21 | 2021-04-21 | 基于多重修正的荧光染色薄膜厚度测量标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113188456A true CN113188456A (zh) | 2021-07-30 |
CN113188456B CN113188456B (zh) | 2022-08-19 |
Family
ID=76977822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110429859.4A Active CN113188456B (zh) | 2021-04-21 | 2021-04-21 | 基于多重修正的荧光染色薄膜厚度测量标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113188456B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6410926B1 (en) * | 1999-10-01 | 2002-06-25 | Ppg Industries Ohio, Inc. | Coating with optical taggent |
CN106525137A (zh) * | 2016-12-07 | 2017-03-22 | 天津大学 | 基于激光诱导荧光的液膜温度场和流场同时测量方法 |
CN106524929A (zh) * | 2016-12-07 | 2017-03-22 | 大连海事大学 | 缸套与活塞润滑油膜分布状态的测试装置及其测试方法 |
JP2017194416A (ja) * | 2016-04-22 | 2017-10-26 | Tpr株式会社 | 油膜厚さ計測装置 |
CN107860322A (zh) * | 2017-12-07 | 2018-03-30 | 东华大学 | 一种液膜厚度测量装置及方法 |
CN108920786A (zh) * | 2018-06-20 | 2018-11-30 | 北京航空航天大学 | 一种基于切比雪夫多项式拟合的区间不确定性分析方法 |
CN109974602A (zh) * | 2019-04-10 | 2019-07-05 | 武汉理工大学 | 基于氙气诱导荧光技术测量轴承水膜厚度的方法 |
CN110006829A (zh) * | 2019-03-12 | 2019-07-12 | 浙江工业大学 | 一种基于最小二乘法的微型光谱仪标定方法 |
CN110044278A (zh) * | 2019-04-09 | 2019-07-23 | 华南理工大学 | 一种非接触式流体薄膜厚度测量装置及方法 |
CN110274546A (zh) * | 2019-07-04 | 2019-09-24 | 哈尔滨工程大学 | 一种基于lif技术的润滑油膜厚度标定方法及装置 |
CN110827362A (zh) * | 2019-11-05 | 2020-02-21 | 天津理工大学 | 一种基于多项式相机响应函数与渐晕效应补偿的光度标定方法 |
-
2021
- 2021-04-21 CN CN202110429859.4A patent/CN113188456B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6410926B1 (en) * | 1999-10-01 | 2002-06-25 | Ppg Industries Ohio, Inc. | Coating with optical taggent |
JP2017194416A (ja) * | 2016-04-22 | 2017-10-26 | Tpr株式会社 | 油膜厚さ計測装置 |
CN106525137A (zh) * | 2016-12-07 | 2017-03-22 | 天津大学 | 基于激光诱导荧光的液膜温度场和流场同时测量方法 |
CN106524929A (zh) * | 2016-12-07 | 2017-03-22 | 大连海事大学 | 缸套与活塞润滑油膜分布状态的测试装置及其测试方法 |
CN107860322A (zh) * | 2017-12-07 | 2018-03-30 | 东华大学 | 一种液膜厚度测量装置及方法 |
CN108920786A (zh) * | 2018-06-20 | 2018-11-30 | 北京航空航天大学 | 一种基于切比雪夫多项式拟合的区间不确定性分析方法 |
CN110006829A (zh) * | 2019-03-12 | 2019-07-12 | 浙江工业大学 | 一种基于最小二乘法的微型光谱仪标定方法 |
CN110044278A (zh) * | 2019-04-09 | 2019-07-23 | 华南理工大学 | 一种非接触式流体薄膜厚度测量装置及方法 |
CN109974602A (zh) * | 2019-04-10 | 2019-07-05 | 武汉理工大学 | 基于氙气诱导荧光技术测量轴承水膜厚度的方法 |
CN110274546A (zh) * | 2019-07-04 | 2019-09-24 | 哈尔滨工程大学 | 一种基于lif技术的润滑油膜厚度标定方法及装置 |
CN110827362A (zh) * | 2019-11-05 | 2020-02-21 | 天津理工大学 | 一种基于多项式相机响应函数与渐晕效应补偿的光度标定方法 |
Non-Patent Citations (3)
Title |
---|
万安军等: "拟合阶次对基于多项式拟合的相位―高度映射关系的影响", 《光电技术应用》 * |
李天宇等: "薄层液膜厚度的点测量和空间测量方法综述", 《实验流体力学》 * |
韩小元等: "X射线荧光光谱法表征薄膜进展", 《光谱学与光谱分析》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113188456B (zh) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108871588B (zh) | 一种红外成像系统多维度联合非均匀性校正方法 | |
EP1333258B1 (en) | Method for correcting sensor output | |
CN112033542B (zh) | 亮度色度测量方法、装置、设备及可读存储介质 | |
CN108986170B (zh) | 一种适用于现场工况的线阵相机平场校正方法 | |
US10852123B2 (en) | Apparatus for critical-dimension localization microscopy | |
CN104931466B (zh) | 基于列处理的plif浓度场标定方法 | |
CN113188456B (zh) | 基于多重修正的荧光染色薄膜厚度测量标定方法 | |
US7435602B2 (en) | Reducing effects of spectral nonuniformity | |
CN113029013B (zh) | 一种考虑统计修正的荧光染色薄膜厚度测量标定方法 | |
CN111504454A (zh) | 一种ccd校正方法、系统、具有ccd校正系统的光谱仪 | |
CN113048897A (zh) | 非线性最小二乘修正的荧光染色薄膜厚度测量标定方法 | |
US8249381B2 (en) | Image based correction for unwanted light signals in a specific region of interest | |
Lenhard | Determination of combined measurement uncertainty via Monte Carlo analysis for the imaging spectrometer ROSIS | |
CN114240801A (zh) | 一种遥感图像非均匀校正方法 | |
CN113125401B (zh) | 考虑荧光成像区域的荧光染色薄膜厚度测量标定方法 | |
CN113048898B (zh) | 一种考虑非线性的荧光染色薄膜厚度测量标定方法 | |
CN114745538A (zh) | 一种显示屏视角和相机暗角联合校正方法及装置 | |
JP2006337181A (ja) | 色校正方法及び色校正用チャート | |
KR102571243B1 (ko) | 색채측정 장치 및 방법 | |
CN112798105A (zh) | 光学测定装置的线性校正方法、光学测定方法以及光学测定装置 | |
CN111458027A (zh) | 一种紫外光谱的测量方法及装置 | |
Som et al. | Photometric Repeatability and Sensitivity Evolution of WFC3/IR | |
US20090045318A1 (en) | Method for reproducibly determining object characteristics | |
CN113612989B (zh) | 一种彩色图像传感器组合曝光下串扰总量测量方法 | |
Zhang et al. | Correcting tristimulus colorimeters via spectrum-enhanced algorithm |
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 |