CN104361226A - 共焦轴向响应曲线峰值位置提取算法 - Google Patents
共焦轴向响应曲线峰值位置提取算法 Download PDFInfo
- Publication number
- CN104361226A CN104361226A CN201410617214.3A CN201410617214A CN104361226A CN 104361226 A CN104361226 A CN 104361226A CN 201410617214 A CN201410617214 A CN 201410617214A CN 104361226 A CN104361226 A CN 104361226A
- Authority
- CN
- China
- Prior art keywords
- data
- response curve
- confocal axial
- curve peak
- extraction algorithm
- 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.)
- Pending
Links
Landscapes
- Microscoopes, Condenser (AREA)
Abstract
共焦轴向响应曲线峰值位置提取算法属于共焦扫描光学测量技术领域;该方法首先归一化数据,从实测数据中选取有效数据;然后以sinc4(a(x-b))函数为目标函数,根据测量参数计算目标参数初值;再以sinc4(a(x-b))函数为目标函数,利用Levenberg-Marquardt算法进行拟合求得目标参数;最后根据目标参数计算峰值位置;本发明共焦轴向响应曲线峰值位置提取算法,通过使用sinc4作为拟合的目标函数来获得共焦轴向相应曲线峰值位置,可以提高数据利用率与测量不确定度,同时减少迭代次数,减少了提取算法的运算时间。
Description
技术领域
共焦轴向响应曲线峰值位置提取算法属于共焦扫描光学测量技术领域。
背景技术
共焦扫描显微镜测量技术是一种用于微米及亚微米尺度测量的三维光学显微技术。因其具有无需样品制备、测量速度快、成本低、不损伤测量表面以及可直接测量高度大于半波长的台阶轮廓等优点,而被广泛用于材料学、微电子学、光学、生物工程学及医学领域的三维微尺度结构测量。共焦扫描秒显微镜测量技术的关键在于共焦轴向相应曲线峰值位置提取,这关系到三维结构测量的精度,常用的提取方法有最大值法、质心法、高斯函数拟合法。前两种方法速度较快但精度不高,高斯函数拟合法精度较高但是速度很慢,并且数据利用率较低。因此,现有共焦轴向相应曲线峰值位置提取方法无法兼顾速度和精度。
发明内容
为了解决上述问题,本发明公开了一种共焦轴向响应曲线峰值位置提取算法,该算法使用sinc4作为拟合的目标函数来获得共焦轴向相应曲线峰值位置,可以提高数据利用率与测量不确定度,同时减少迭代次数,减少了提取算法的运算时间。
本发明的目的是这样实现的:
一种共焦轴向响应曲线峰值位置提取算法,包括以下步骤:
步骤a、归一化数据,从实测数据中选取有效数据,轴向扫描间隔为z,扫描层数为N;
步骤b、以sinc4(a(x-b))函数为目标函数,根据测量参数计算目标参数初值a1和b1;
步骤c、以sinc4(a(x-b))函数为目标函数,利用Levenberg-Marquardt算法进行拟合求得目标参数;
步骤d、根据目标参数计算峰值位置。
上述共焦轴向相应峰值位置提取算法,所述步骤a中归一化数据方法为,将测得的原始数据除以该原始数据最大值。
上述共焦轴向相应峰值位置提取算法,所述步骤a中有效数据为原始数据中大于原始数据最大值50%的数据以及其索引值,或归一化数据中大于归一化数据最大值50%的数据以及其索引值。
上述共焦轴向响应曲线峰值位置提取算法,所述步骤b中a1=π/(N'×z),b1为原始数据最大值所对应的索引值;
所述参数b'为原始数据或归一化数据最大值的索引值,所述参数z为共焦轴向相应曲线半高宽度。
上述共焦轴向响应曲线峰值位置提取算法,所述步骤c中Levenberg-Marquardt算法的最大迭代次数N>10000,残差res<1×10-5。
上述共焦轴向响应曲线峰值位置提取算法,所述步骤d中峰值位置为b×轴向扫描间隔。
有益效果:
本发明以sinc4函数作为共焦轴向相应曲线拟合的目标函数,并基于此目标函数计算目标参数的初值,并拟合求得目标参数,最后得到峰值位置,由于sinc4函数更接近理论模型,因此可以提高数据利用率与测量不确定度,进而减少迭代次数,实现减少提取算法的运算时间的技术目的。
附图说明
图1是本发明共焦轴向响应曲线峰值位置提取算法流程图。
图2为拟合结果示意图。
具体实施方式
下面结合附图和具体实施例进行详细说明,以便对本发明的目的、技术方案有更深入的理解。
本实施例的共焦轴向响应曲线峰值位置提取算法,流程图如图1所示。该算法包括以下步骤:
步骤a、归一化数据,具体为将测得的原始数据除以该原始数据最大值;从实测数据中选取有效数据,记为(x1,y1),(x2,y2),...,(xn,yn),轴向扫描间隔为z=0.1μm,扫描层数为N=200层;其中,所述的有效数据为原始数据中大于原始数据最大值50%的数据以及其索引值;
步骤b、以sinc4(a(x-b))函数为目标函数,根据测量参数计算目标参数初值:
a的初值a1为π/(N'×z),N'=38,z=0.1,即a1=0.826;
b的初值b1为最大值对应的索引值47;
步骤c、以sinc4(a(x-b))函数为目标函数,利用Levenberg-Marquardt算法进行拟合求得目标参数b=47.67,最大迭代次数N>10000,残差res<1×10-5;
步骤d、样品该点的相对高对记为b×z=4.767;
拟合结果如图2所示。
对每个扫描点按照上述步骤处理,即得每个扫描点的相对高度,由此可得样品的三维形貌。
这里需要说明的是,同sinc2(a(x-b))相比,采用sinc4(a(x-b))函数为目标函数进行拟合,对于表面平滑的样品,测量结果对于样品三维形貌的表达更加准确。
本发明不局限于上述最佳实施方式,任何人应该得知在本发明的启示下作出的结构变化或方法改进,凡是与本发明具有相同或相近的技术方案,均落入本发明的保护范围之内。
Claims (7)
1.一种共焦轴向响应曲线峰值位置提取算法,其特征在于,包括以下步骤:
步骤a、归一化数据,从实测数据中选取有效数据,轴向扫描间隔为z,扫描层数为N;
步骤b、以sinc4(a(x-b))函数为目标函数,根据测量参数计算目标参数初值a1和b1;
步骤c、以sinc4(a(x-b))函数为目标函数,利用Levenberg-Marquardt算法拟合求得目标参数;
步骤d、根据目标参数计算峰值位置。
2.根据权利要求1所述的共焦轴向相应峰值位置提取算法,其特征在于,所述步骤a中归一化数据方法为,将测得的原始数据除以该原始数据最大值。
3.根据权利要求1所述的共焦轴向相应峰值位置提取算法,其特征在于,所述步骤a中有效数据为原始数据中大于原始数据最大值50%的数据以及其索引值,或归一化数据中大于归一化数据最大值50%的数据以及其索引值。
4.根据权利要求1、2或3所述的共焦轴向响应曲线峰值位置提取算法,其特征在于,所述步骤b中a1=π/(N'×z),b1为原始数据最大值所对应的索引值。
5.根据权利要求4所述的共焦轴向响应曲线峰值位置提取算法,其特征在于,所述参数N'为原始数据或归一化数据半高宽度,即为值大于二分之一最大值的数据的个数。
6.根据权利要求1、2或3所述的共焦轴向响应曲线峰值位置提取算法,其特征在于,所述步骤c中Levenberg-Marquardt算法的最大迭代次数N>10000,残差res<1×10-5。
7.根据权利要求1、2或3所述的共焦轴向响应曲线峰值位置提取算法,其特征在于,所述步骤d中峰值位置为b×z。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410617214.3A CN104361226A (zh) | 2014-11-05 | 2014-11-05 | 共焦轴向响应曲线峰值位置提取算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410617214.3A CN104361226A (zh) | 2014-11-05 | 2014-11-05 | 共焦轴向响应曲线峰值位置提取算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104361226A true CN104361226A (zh) | 2015-02-18 |
Family
ID=52528485
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410617214.3A Pending CN104361226A (zh) | 2014-11-05 | 2014-11-05 | 共焦轴向响应曲线峰值位置提取算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104361226A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093718A (zh) * | 2016-05-30 | 2016-11-09 | 中国南方电网有限责任公司电网技术研究中心 | 一种绝缘子闪络风险预警方法 |
CN108875124A (zh) * | 2018-04-26 | 2018-11-23 | 哈尔滨工业大学 | 提取共焦轴向响应曲线峰值位置的极大值补偿算法 |
CN109540806A (zh) * | 2018-10-16 | 2019-03-29 | 华中科技大学 | 一种利用动态阈值重心法提取共焦显微峰值的方法 |
CN110398768A (zh) * | 2019-07-15 | 2019-11-01 | 华中科技大学 | 一种基于像素电离室的束斑动态监测方法及系统 |
CN110749280A (zh) * | 2019-09-30 | 2020-02-04 | 华中科技大学 | 一种提取峰值位置索引坐标的方法、系统及计算机可读介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100321773A1 (en) * | 2009-06-19 | 2010-12-23 | Industrial Technology Research Institute | Method and system for three-dimensional polarization-based confocal microscopy |
CN101929848A (zh) * | 2010-06-30 | 2010-12-29 | 北京理工大学 | 具有高空间分辨力的乘积共焦扫描检测方法 |
CN102759328A (zh) * | 2012-07-05 | 2012-10-31 | 哈尔滨工业大学 | 基于椭球反射双通照明差动共焦测量装置与方法 |
US20130235385A1 (en) * | 2012-03-09 | 2013-09-12 | Canon Kabushiki Kaisha | Surface shape measurement method and measurement apparatus |
-
2014
- 2014-11-05 CN CN201410617214.3A patent/CN104361226A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100321773A1 (en) * | 2009-06-19 | 2010-12-23 | Industrial Technology Research Institute | Method and system for three-dimensional polarization-based confocal microscopy |
CN101929848A (zh) * | 2010-06-30 | 2010-12-29 | 北京理工大学 | 具有高空间分辨力的乘积共焦扫描检测方法 |
US20130235385A1 (en) * | 2012-03-09 | 2013-09-12 | Canon Kabushiki Kaisha | Surface shape measurement method and measurement apparatus |
CN102759328A (zh) * | 2012-07-05 | 2012-10-31 | 哈尔滨工业大学 | 基于椭球反射双通照明差动共焦测量装置与方法 |
Non-Patent Citations (2)
Title |
---|
刘俭: "《超分辨移相共焦和复色差动共焦扫描显微技术及理论》", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
黄久荣: "《曲面基底衍射元件共焦扫描轮廓提取方法研究》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093718A (zh) * | 2016-05-30 | 2016-11-09 | 中国南方电网有限责任公司电网技术研究中心 | 一种绝缘子闪络风险预警方法 |
CN106093718B (zh) * | 2016-05-30 | 2019-11-22 | 中国南方电网有限责任公司电网技术研究中心 | 一种绝缘子闪络风险预警方法 |
CN108875124A (zh) * | 2018-04-26 | 2018-11-23 | 哈尔滨工业大学 | 提取共焦轴向响应曲线峰值位置的极大值补偿算法 |
CN108875124B (zh) * | 2018-04-26 | 2022-04-01 | 哈尔滨工业大学 | 提取共焦轴向响应曲线峰值位置的极大值补偿算法 |
CN109540806A (zh) * | 2018-10-16 | 2019-03-29 | 华中科技大学 | 一种利用动态阈值重心法提取共焦显微峰值的方法 |
CN109540806B (zh) * | 2018-10-16 | 2019-11-05 | 华中科技大学 | 一种利用动态阈值重心法提取共焦显微峰值的方法 |
CN110398768A (zh) * | 2019-07-15 | 2019-11-01 | 华中科技大学 | 一种基于像素电离室的束斑动态监测方法及系统 |
CN110749280A (zh) * | 2019-09-30 | 2020-02-04 | 华中科技大学 | 一种提取峰值位置索引坐标的方法、系统及计算机可读介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104361226A (zh) | 共焦轴向响应曲线峰值位置提取算法 | |
CN107301648B (zh) | 基于重叠区域边界角度的冗余点云去除方法 | |
CN103989462B (zh) | 一种脉搏波形第一特征点和第二特征点的提取方法 | |
CN103487074B (zh) | 利用三点寻峰算法处理fbg传感信号的方法 | |
CN107886529A (zh) | 一种用于三维重建的点云配准方法 | |
CN103870845A (zh) | 点云聚类去噪过程中新的k值优化方法 | |
CN104634460A (zh) | 一种分布式fbg传感网络多峰自适应精确寻峰方法 | |
CN104463894A (zh) | 一种多视角三维激光点云全局优化整体配准方法 | |
CN103714166B (zh) | 激光雷达电力线点云数据的拟合方法 | |
CN102799763B (zh) | 一种基于点云姿态标准化的点云线特征提取方法 | |
CN105095540B (zh) | 一种井间相对位置关系的法面扫描解析方法 | |
CN104684081A (zh) | 基于距离聚类遴选锚节点的无线传感器网络节点定位算法 | |
CN107123134A (zh) | 一种基于特征的危岩体滑坡监测方法 | |
CN104008257A (zh) | 一种针对具有复杂曲面的构件的逆向设计方法 | |
CN112013788A (zh) | 基于叶片局部前缘曲线特征标定转动中心的方法 | |
CN106446415A (zh) | 一种基于点云数据的变电站支柱绝缘子中心线提取方法 | |
CN109448042A (zh) | 一种利用积分曲线计算血管壁厚度的方法 | |
CN104315994A (zh) | 共焦轴向响应曲线峰值位置提取算法 | |
CN105115441A (zh) | 一种回转体零件廓形的特征点提取及自动分段方法 | |
CN105089597A (zh) | 一种裂缝复杂程度评价方法 | |
CN105045958A (zh) | 基于bp神经网络的gps高程拟合的实现系统及方法 | |
CN103995972A (zh) | 一种地层对比方法 | |
CN109101778B (zh) | 基于性能退化数据和寿命数据融合的Wiener过程参数估计方法 | |
CN105095649B (zh) | 建筑物倾斜角度的数据分析方法和系统 | |
CN107292924B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150218 |
|
WD01 | Invention patent application deemed withdrawn after publication |