CN116304645A - 一种基于模态分解的重叠峰提取的方法及装置 - Google Patents

一种基于模态分解的重叠峰提取的方法及装置 Download PDF

Info

Publication number
CN116304645A
CN116304645A CN202310589272.9A CN202310589272A CN116304645A CN 116304645 A CN116304645 A CN 116304645A CN 202310589272 A CN202310589272 A CN 202310589272A CN 116304645 A CN116304645 A CN 116304645A
Authority
CN
China
Prior art keywords
variables
unconstrained
peak
formula
extremum
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
CN202310589272.9A
Other languages
English (en)
Other versions
CN116304645B (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.)
Optosky Xiamen Optoelectronic Co ltd
Original Assignee
Optosky Xiamen Optoelectronic 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 Optosky Xiamen Optoelectronic Co ltd filed Critical Optosky Xiamen Optoelectronic Co ltd
Priority to CN202310589272.9A priority Critical patent/CN116304645B/zh
Publication of CN116304645A publication Critical patent/CN116304645A/zh
Application granted granted Critical
Publication of CN116304645B publication Critical patent/CN116304645B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • 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/02Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • G01B11/06Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开一种基于模态分解的重叠峰提取的方法及装置,其包括:根据类间方差函数表示得到求峰值的问题公式,使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题;将无约束极值问题进行二次优化,转化为最小化问题公式;根据最小化问题公式,设置收敛条件,对其中的3个变量进行更新,直到满足预设的收敛条件时,停止迭代;停止迭代时,最终的计算结果即为两个光谱轴向响应信号的峰值。本发明采用上述方案,可以在多个峰值具有重叠区域时分离它们,准确的提取出峰值波长,从而解决测量透明薄板中,重叠峰难以提取的问题。

Description

一种基于模态分解的重叠峰提取的方法及装置
技术领域
本发明涉及光学测量技术领域,具体是一种基于模态分解的重叠峰提取的方法及装置。
背景技术
光谱共焦位移传感器是通过彩色激光光源发射出一束高密度宽光谱光,通过色散镜头后,在量程范围内形成不同波长的白光,每个波长对应一个距离值,测量光射到物体表面反射回来,只有满足共聚焦条件的光,可以通过小孔被光谱仪感测到。通过计算被感测到光的焦点的波长,换算获的距离值。由于其具有的高精度与非接触性,被广泛的应用于测量玻璃基板的厚度。在测量玻璃基板的厚度过程中准确的提取光谱轴向响应信号(sARS)的峰值起着至关重要的作用。
目前有多种针对sARS的峰值位置定位技术,如经典的质心法和基于模型的拟合算法:高斯拟合算法(GFA)、二项式拟合算法(BFA)。然而,这些技术都是对一个单独的光谱轴向响应信号(sARS)进行提取峰值。当被测物体为玻璃基板时,白光在玻璃上表面反射回来一次,在玻璃下表面再反射回来一次,将得到两个光谱轴向响应信号(sARS)。当玻璃基板较厚时,两个光谱轴向响应信号(sARS)将会是相互独立分布,如图1所示为归一化后的玻璃基板原始光强。此时使用上述技术,分别对两个光谱轴向响应信号(sARS)进行峰值提取,还是可以较为准确的得到两个峰值。而当玻璃基板较薄时,两个光谱轴向响应信号(sARS)会相互干扰,形成重叠区域,如图2所示为滤波后的玻璃基板光强。此时,没有波双边数据点的支持,上述技术提取的重叠峰值将会产生极大的误差。
发明内容
在下文中给出了关于本发明实施例的简要概述,以便提供关于本发明的某些方面的基本理解。应当理解,以下概述并不是关于本发明的穷举性概述。它并不是意图确定本发明的关键或重要部分,也不是意图限定本发明的范围。其目的仅仅是以简化的形式给出某些概念,以此作为稍后论述的更详细描述的前序。
为解决上述技术问题,本发明提供一种基于模态分解的重叠峰提取的方法,使用模态分解技术,将两个光谱轴向响应信号(sARS)从原始光强中分解出来,迭代求解出两个光谱轴向响应信号(sARS)的峰值,提高寻峰准确性。
具体的,根据本申请的一方面,提供一种基于模态分解的重叠峰提取的方法,用于在测量玻璃基板的厚度过程中提取两个光谱轴向响应信号sARS的峰值,包括:
根据类间方差函数表示得到求峰值的问题公式,该问题公式含有2个变量和1个约束条件;
使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题公式;
将无约束极值问题进行二次优化,转化为等效最小化问题公式,该最小化问题公式含有3个变量;
根据最小化问题公式,设置收敛条件,对其中的3个变量进行更新,直到满足预设的收敛条件时,停止迭代;停止迭代时,最终的计算结果即为两个光谱轴向响应信号(sARS)的峰值。
进一步的,该基于模态分解的重叠峰提取的方法具体包括:
步骤1:根据类间方差函数表示法,将求峰值的问题表示为:
Figure SMS_1
(公式1)
其中
Figure SMS_2
为两个模态,/>
Figure SMS_3
光强序列中的波长,/>
Figure SMS_4
为峰值波长,/>
Figure SMS_5
和/>
Figure SMS_6
分别是光谱仪的最小和最大检测波长;约束条件表示为:
Figure SMS_7
(公式2)
步骤2:使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子
Figure SMS_8
,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题:
Figure SMS_9
(公式3)
其中,
Figure SMS_10
为半峰全宽,/>
Figure SMS_11
为拉格朗日乘子;
步骤3:为了避免函数的严格凸假设,增加迭代过程的鲁棒性,将步骤2的无约束极值问题进行二次优化,转化为如下等效的最小化问题:
Figure SMS_12
(公式4)
其中,
Figure SMS_13
为迭代的次数;
步骤4:步骤3的最小化问题的解为:
Figure SMS_14
(公式5)
步骤5:使用质心算法更新峰值波长
Figure SMS_15
Figure SMS_16
(公式6)
步骤6:拉格朗日乘子
Figure SMS_17
根据以下公式更新:
Figure SMS_18
(公式7)
其中,
Figure SMS_19
为噪声容量参数,根据数据包含的噪声多少来确定大小;
步骤7:设置收敛条件:
Figure SMS_20
(公式8)
其中
Figure SMS_21
为给定的误差,给定允许误差优选为10e-4;
步骤8:公式4为寻峰方程,公式5、6、7分别为变量
Figure SMS_23
、/>
Figure SMS_26
和/>
Figure SMS_28
的更新公式,公式8为收敛条件,重复步骤步骤4、步骤5和步骤6从而更新模态/>
Figure SMS_24
、中心波长/>
Figure SMS_25
和拉格朗日乘数/>
Figure SMS_27
,直到满足步骤7中的收敛条件时,停止迭代;停止迭代时,最终的/>
Figure SMS_29
和/>
Figure SMS_22
的计算结果即为两个光谱轴向响应信号(sARS)的峰值。
为提高寻峰准确性,本发明创新的将图像处理领域的“类间方差函数”用于表示求峰值,并在求峰值的过程中,同时使用二次惩罚项和拉格朗日乘子,以使问题不受约束(步骤2)。其中,二次惩罚项是鼓励重建保真度的经典方法,拉格朗日乘子是严格执行约束的常用方法,在有限权重下二次惩罚有良好收敛,拉格朗日乘数会对约束做严格执行,因此步骤2中采用二次惩罚项和拉格朗日乘子这两个项的组合会取得更好的效果。此外,在实际求解时,步骤3的最小化问题的解是比较困难的,本发明人致力于研究该最小化问题的算法,才最终得到步骤4的解。
优选的,所述步骤1在根据类间方差函数表示得到求峰值的问题之前,还包括如下过程:
采集光谱共焦位移传感器的光谱共焦信号,执行扣除暗电流处理,再将光强做归一化处理(把光强数据映射到0-1之间),获得点波长序列对应光强序列
Figure SMS_30
进一步的,所述步骤1在根据类间方差函数表示得到求峰值的问题之前,还包括如下过程:
用最小二乘法对原始光强序列进行滑动拟合,并求出卷积系数。用卷积系数对原始光强序列
Figure SMS_31
进行卷积计算,完成sg滤波,得到去噪后的光强序列/>
Figure SMS_32
根据本申请的另一方面,还提供一种基于模态分解的重叠峰提取的装置,其执行上述基于模态分解的重叠峰提取的方法,该装置包括:
求峰值问题公式获得模块,用于根据类间方差函数表示得到求峰值的问题公式,该问题公式含有2个变量和1个约束条件;
无约束极值问题公式获得模块,用于使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题;
最小化问题公式获得模块,用于将无约束极值问题进行二次优化,转化为等效最小化问题公式,该最小化问题公式含有3个变量;
计算模块,用于根据最小化问题公式,设置收敛条件,对其中的3个变量进行更新,直到满足预设的收敛条件时,停止迭代;停止迭代时,最终的计算结果即为两个光谱轴向响应信号(sARS)的峰值。
针对现有技术的以上缺陷或改进需求,本发明提出了一种自适应模态分解方法,可以在多个峰值具有重叠区域时分离它们,准确的提取出峰值波长,从而解决测量透明薄板中,重叠峰难以提取的问题。在图像处理领域,最大类间方差法可以很好的对直方图为双峰的图片进行阈值分割,其中类间方差表征的正是阈值两侧的数据的离散程度,类间方差越大,说明前景和背景相差越大,二值化效果越佳,因此本发明引入最大类间方差法来寻找重叠区域的峰值。此外,通过模态分解、迭代求解峰值,可以更精准得提取出测量薄玻璃时的两个峰值波长,因此,本发明的模态分解方法显示了准确性的优势,特别是当多个峰具有重叠区域时,拥有更强的寻峰准确性。
附图说明
本发明可以通过参考下文中结合附图所给出的描述而得到更好的理解,其中在所有附图中使用了相同或相似的附图标记来表示相同或者相似的部件。所述附图连同下面的详细说明一起包含在本说明书中并且形成本说明书的一部分,而且用来进一步举例说明本发明的优选实施例和解释本发明的原理和优点。在附图中:
图1为归一化后的未滤波的玻璃基板原始光强;
图2为滤波后的玻璃基板光强;
图3为本发明实施例的重叠峰提取的方法的迭代过程示意图;
图4为本发明实施例的重叠峰提取的方法的实验结果示意图。
具体实施方式
下面将参照附图来说明本发明的实施例。在本发明的一个附图或一种实施方式中描述的元素和特征可以与一个或更多个其它附图或实施方式中示出的元素和特征相结合。应当注意,为了清楚的目的,附图和说明中省略了与本发明无关的、本领域普通技术人员已知的部件和处理的表示和描述。
现技术方案多为:使用线性拟合算法(如质心法)或非线性拟合算法(如高斯拟合、二次多项式函数拟合)分别提取重叠峰的两个峰值。当透明薄板比较薄时,两个亚反射强度分布不再相互独立分布,使用单独的寻峰算法,无法避免两个亚反射强度相互干扰,使得寻峰准确性较差。本发明使用模态分解技术,将两个光谱轴向响应信号(sARS)从原始光强中分解出来,迭代求解出两个光谱轴向响应信号(sARS)的峰值,提高寻峰准确性。
本发明的基于模态分解的重叠峰提取的方法,可以在多个峰值具有重叠区域时分离它们,准确的提取出峰值波长。解决测量透明薄板中,重叠峰难以提取的问题。该方法显示了准确性的优势,尤其适用于多个峰具有重叠区域时。
具体的,参见图3的流程图,本发明的基于模态分解的重叠峰提取的方法包括以下步骤:
STEP1:采集光谱共焦位移传感器的光谱共焦信号,执行扣除暗电流处理,再将光强做归一化处理(把光强数据映射到0-1之间),获得点波长序列对应光强序列
Figure SMS_33
;
STEP2:用最小二乘法对原始光强序列进行滑动拟合,并求出卷积系数。用卷积系数对原始光强序列
Figure SMS_34
进行卷积计算,完成sg滤波,得到去噪后的光强序列/>
Figure SMS_35
。此步骤目的在于:消除数据中的噪声。滤波操作非本发明的重点,也可以使用其他滤波器去噪,如中值滤波等。
STEP3:在图像处理领域,最大类间方差法可以很好的对直方图为双峰的图片进行阈值分割,其中类间方差表征的正是阈值两侧的数据的离散程度,类间方差越大,说明前景和背景相差越大,二值化效果越佳。本发明引入最大类间方差法,来寻找重叠区域的峰值。根据类间方差函数表示法,将求峰值的问题表示为:
Figure SMS_36
其中
Figure SMS_37
为两个模态,/>
Figure SMS_38
光强序列中的波长,/>
Figure SMS_39
为峰值波长,/>
Figure SMS_40
和/>
Figure SMS_41
分别是光谱仪的最小和最大检测波长。约束条件可以表示为:
Figure SMS_42
STEP4:此时,方程转换成了有约束的优化问题。使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子
Figure SMS_43
,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题:
Figure SMS_44
其中,
Figure SMS_45
为半峰全宽,/>
Figure SMS_46
为拉格朗日乘子。
STEP5:为了避免函数的严格凸假设,增加迭代过程的鲁棒性,转化为如下等效的最小化问题:
Figure SMS_47
其中,
Figure SMS_48
为迭代的次数。
STEP6:这个二次优化问题的解为:
Figure SMS_49
STEP7:使用质心算法更新峰值波长
Figure SMS_50
,当然也可以使用其他非线性算法来更新峰值波长/>
Figure SMS_51
Figure SMS_52
STEP8:拉格朗日乘子
Figure SMS_53
根据以下公式更新:
Figure SMS_54
其中,
Figure SMS_55
为噪声容量参数,可以根据数据包含的噪声多少,适当选择大小。
STEP9:设置收敛条件:
Figure SMS_56
其中
Figure SMS_57
为给定的误差,给定允许误差优选为10e-4,由于无法避免噪声的干扰,所以误差不能设为0,而再小于上述给定误差时,在迭代计算中,噪声将融入波长计算中,对数据准确性造成影响。
STEP10:至此在步骤STEP5中得到寻峰方程,在骤STEP6、STEP7和STEP8中得到了三个变量更新的公式,在STEP9中得到了收敛条件。重复步骤STEP6、STEP7和STEP8,更新模态
Figure SMS_58
、中心波长/>
Figure SMS_59
和拉格朗日乘数/>
Figure SMS_60
,当满足STEP9中的收敛条件时,停止迭代。停止迭代时,最终的/>
Figure SMS_61
和/>
Figure SMS_62
计算结果即为两个光谱轴向响应信号(sARS)的峰值。
上述过程中,还可分别对两个光谱轴向响应信号(sARS)的峰值使用质心法或高斯拟合法寻找波峰。
本发明通过上述方案,其在STEP3中使用最大类间方差函数,表示求解峰值的表示公式,在STEP4中使用增广拉格朗日乘数法求解峰值的表示公式;并分别在STEP6、STEP7和STEP8设计了三个参数(模态
Figure SMS_63
、中心波长/>
Figure SMS_64
和拉格朗日乘数/>
Figure SMS_65
)的更新公式,结合STEP9中的收敛条件的设定公式,最终获得两个光谱轴向响应信号(sARS)的峰值。
如图4所示,在仿真实验中:使用高斯拟合直接计算两个光谱轴向响应信号(sARS)的峰值,与真实的峰值偏差较大,而本发明提出的模态分解方法更接近真实值,拥有更强的寻峰准确性。
本发明实施例还提供一种基于模态分解的重叠峰提取的装置,其执行上述基于模态分解的重叠峰提取的方法,该装置包括:
求峰值问题公式获得模块,用于根据类间方差函数表示得到求峰值的问题公式,该问题公式含有2个变量和1个约束条件;
无约束极值问题公式获得模块,用于使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题;
最小化问题公式获得模块,用于将无约束极值问题进行二次优化,转化为等效最小化问题公式,该最小化问题公式含有3个变量;
计算模块,用于根据最小化问题公式,设置收敛条件,对其中的3个变量进行更新,直到满足预设的收敛条件时,停止迭代;停止迭代时,最终的计算结果即为两个光谱轴向响应信号(sARS)的峰值。
此外,本发明的方法不限于按照说明书中描述的时间顺序来执行,也可以按照其他的时间顺序地、并行地或独立地执行。因此,本说明书中描述的方法的执行顺序不对本发明的技术范围构成限制。
尽管上面已经通过对本发明的具体实施例的描述对本发明进行了披露,但是,应该理解,上述的所有实施例和示例均是示例性的,而非限制性的。本领域的技术人员可在所附权利要求的精神和范围内设计对本发明的各种修改、改进或者等同物。这些修改、改进或者等同物也应当被认为包括在本发明的保护范围内。

Claims (5)

1.一种基于模态分解的重叠峰提取的方法,用于在测量玻璃基板的厚度过程中提取两个光谱轴向响应信号sARS的峰值,其特征在于:包括:
根据类间方差函数表示得到求峰值的问题公式,该问题公式含有2个变量和1个约束条件;
使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题;
将无约束极值问题进行二次优化,转化为等效最小化问题,该最小化问题公式含有3个变量;
根据最小化问题公式,设置收敛条件,对其中的3个变量进行更新,直到满足预设的收敛条件时,停止迭代;停止迭代时,最终的计算结果即为两个光谱轴向响应信号的峰值。
2.根据权利要求1所述的基于模态分解的重叠峰提取的方法,其特征在于:该基于模态分解的重叠峰提取的方法具体包括:
步骤1:根据类间方差函数表示法,将求峰值的问题表示为:
Figure QLYQS_1
其中
Figure QLYQS_2
为两个模态,/>
Figure QLYQS_3
光强序列中的波长,/>
Figure QLYQS_4
为去噪后的光强序列,/>
Figure QLYQS_5
为峰值波长,/>
Figure QLYQS_6
Figure QLYQS_7
分别是光谱仪的最小和最大检测波长;约束条件表示为:
Figure QLYQS_8
步骤2:使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子
Figure QLYQS_9
,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题:
Figure QLYQS_10
其中,
Figure QLYQS_11
为半峰全宽,/>
Figure QLYQS_12
为拉格朗日乘子;
步骤3:为了避免函数的严格凸假设,增加迭代过程的鲁棒性,将步骤2的无约束极值问题进行二次优化,转化为如下等效的最小化问题:
Figure QLYQS_13
其中,
Figure QLYQS_14
为迭代的次数;
步骤4:步骤3的最小化问题的解为:
Figure QLYQS_15
步骤5:使用质心算法更新峰值波长
Figure QLYQS_16
Figure QLYQS_17
步骤6:拉格朗日乘子
Figure QLYQS_18
根据以下公式更新:
Figure QLYQS_19
其中,
Figure QLYQS_20
为噪声容量参数,根据数据包含的噪声多少来确定大小;
步骤7:设置收敛条件:
Figure QLYQS_21
其中
Figure QLYQS_22
为给定的误差,给定允许误差优选为10e-4;
步骤8:重复步骤步骤4、步骤5和步骤6从而更新模态
Figure QLYQS_23
、中心波长/>
Figure QLYQS_24
和拉格朗日乘数/>
Figure QLYQS_25
,直到满足步骤7中的收敛条件时,停止迭代;停止迭代时,最终的/>
Figure QLYQS_26
和/>
Figure QLYQS_27
的计算结果即为两个光谱轴向响应信号(sARS)的峰值。
3.根据权利要求2所述的基于模态分解的重叠峰提取的方法,其特征在于:所述步骤1在根据类间方差函数表示得到求峰值的问题之前,还包括如下过程:
采集光谱共焦位移传感器的光谱共焦信号,执行扣除暗电流处理,再将光强做归一化处理,获得点波长序列对应光强序列
Figure QLYQS_28
4.根据权利要求1所述的基于模态分解的重叠峰提取的方法,其特征在于:所述步骤1在根据类间方差函数表示得到求峰值的问题之前,还包括如下过程:
用最小二乘法对原始光强序列进行滑动拟合,并求出卷积系数;用卷积系数对原始光强序列
Figure QLYQS_29
进行卷积计算,完成sg滤波,得到去噪后的光强序列/>
Figure QLYQS_30
5.一种基于模态分解的重叠峰提取的装置,其特征在于:执行权利要求1-4任一所述的基于模态分解的重叠峰提取的方法,该装置包括:
求峰值问题公式获得模块,用于根据类间方差函数表示得到求峰值的问题公式,该问题公式含有2个变量和1个约束条件;
无约束极值问题公式获得模块,用于使用增广拉格朗日乘数法,给函数加上惩罚项和拉格朗日乘子,将含有2个变量和1个约束条件的最优化问题转化为一个含有3个变量的无约束极值问题;
最小化问题公式获得模块,用于将无约束极值问题进行二次优化,转化为等效最小化问题公式,该最小化问题公式含有3个变量;
计算模块,用于根据最小化问题公式,设置收敛条件,对其中的3个变量进行更新,直到满足预设的收敛条件时,停止迭代;停止迭代时,最终的计算结果即为两个光谱轴向响应信号(sARS)的峰值。
CN202310589272.9A 2023-05-24 2023-05-24 一种基于模态分解的重叠峰提取的方法及装置 Active CN116304645B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310589272.9A CN116304645B (zh) 2023-05-24 2023-05-24 一种基于模态分解的重叠峰提取的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310589272.9A CN116304645B (zh) 2023-05-24 2023-05-24 一种基于模态分解的重叠峰提取的方法及装置

Publications (2)

Publication Number Publication Date
CN116304645A true CN116304645A (zh) 2023-06-23
CN116304645B CN116304645B (zh) 2023-08-15

Family

ID=86785520

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310589272.9A Active CN116304645B (zh) 2023-05-24 2023-05-24 一种基于模态分解的重叠峰提取的方法及装置

Country Status (1)

Country Link
CN (1) CN116304645B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090244636A1 (en) * 2005-09-08 2009-10-01 X-Rite, Inc. Systems and method for optical scatter imaging of latent image plates
CN106447688A (zh) * 2016-03-31 2017-02-22 大连海事大学 一种高光谱溢油图像的有效分割方法
US20190073444A1 (en) * 2016-03-10 2019-03-07 Genomic Vision Method for analyzing a sequence of target regions and detect anomalies
CN109886464A (zh) * 2019-01-20 2019-06-14 东北电力大学 基于优化奇异值分解生成特征集的低信息损失短期风速预测方法
CN110727915A (zh) * 2019-10-17 2020-01-24 电子科技大学 一种基于数据相关约束的鲁棒自适应波束形成方法
US20200210826A1 (en) * 2018-12-29 2020-07-02 Northeastern University Intelligent analysis system using magnetic flux leakage data in pipeline inner inspection
CN111879726A (zh) * 2020-08-26 2020-11-03 中国烟草总公司郑州烟草研究院 基于加工前后同步近红外分析的烟草热加工强度及波动性在线监测方法
WO2021068939A1 (zh) * 2019-10-12 2021-04-15 中科新松有限公司 一种基于多分量信号分解的机械臂关节振动识别方法
WO2021109011A1 (zh) * 2019-12-04 2021-06-10 电子科技大学 一种基于超声图像的电容内部缺陷智能检测方法
CN115392311A (zh) * 2022-08-26 2022-11-25 广西民族大学 基于变分模态分解的高效甘蔗茎节识别方法
CN116089824A (zh) * 2023-03-30 2023-05-09 奥谱天成(厦门)光电有限公司 光谱共焦位移传感器的峰值提取方法及系统、介质

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090244636A1 (en) * 2005-09-08 2009-10-01 X-Rite, Inc. Systems and method for optical scatter imaging of latent image plates
US20190073444A1 (en) * 2016-03-10 2019-03-07 Genomic Vision Method for analyzing a sequence of target regions and detect anomalies
CN106447688A (zh) * 2016-03-31 2017-02-22 大连海事大学 一种高光谱溢油图像的有效分割方法
US20200210826A1 (en) * 2018-12-29 2020-07-02 Northeastern University Intelligent analysis system using magnetic flux leakage data in pipeline inner inspection
CN109886464A (zh) * 2019-01-20 2019-06-14 东北电力大学 基于优化奇异值分解生成特征集的低信息损失短期风速预测方法
WO2021068939A1 (zh) * 2019-10-12 2021-04-15 中科新松有限公司 一种基于多分量信号分解的机械臂关节振动识别方法
CN110727915A (zh) * 2019-10-17 2020-01-24 电子科技大学 一种基于数据相关约束的鲁棒自适应波束形成方法
WO2021109011A1 (zh) * 2019-12-04 2021-06-10 电子科技大学 一种基于超声图像的电容内部缺陷智能检测方法
CN111879726A (zh) * 2020-08-26 2020-11-03 中国烟草总公司郑州烟草研究院 基于加工前后同步近红外分析的烟草热加工强度及波动性在线监测方法
CN115392311A (zh) * 2022-08-26 2022-11-25 广西民族大学 基于变分模态分解的高效甘蔗茎节识别方法
CN116089824A (zh) * 2023-03-30 2023-05-09 奥谱天成(厦门)光电有限公司 光谱共焦位移传感器的峰值提取方法及系统、介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
武尧;刘振宇;谷亚宁;: "基于机器视觉的多目标sar图像阈值分割仿真", 计算机仿真, no. 10 *
马志贤;吴中耀;游寒旭;朱杰;: "基于SVM的X射线天文图像点源探测算法", 上海师范大学学报(自然科学版), no. 02 *

Also Published As

Publication number Publication date
CN116304645B (zh) 2023-08-15

Similar Documents

Publication Publication Date Title
US9819879B2 (en) Image filtering apparatus and method based on noise prediction using infrared ray (IR) intensity
Hu et al. Improved baseline correction method based on polynomial fitting for Raman spectroscopy
US7657080B2 (en) Method and apparatus for producing an image containing depth information
Paulus et al. Technical workflows for hyperspectral plant image assessment and processing on the greenhouse and laboratory scale
Qian et al. A fully automated algorithm of baseline correction based on wavelet feature points and segment interpolation
Regeling et al. Development of an image pre‐processor for operational hyperspectral laryngeal cancer detection
EP3214831B1 (en) Image formation system, image formation method, imaging element, and program
TW201435298A (zh) 圖案評估裝置及具備圖案評估裝置之外觀檢查裝置
CN113109317B (zh) 基于背景扣除提取峰面积的拉曼光谱定量分析方法及系统
CN116089824B (zh) 光谱共焦位移传感器的峰值提取方法及系统、介质
WO2020258434A1 (zh) 一种基于tie的相位成像方法、装置及可读存储介质
Pölönen et al. Convolutional neural networks in skin cancer detection using spatial and spectral domain
WO2021014764A1 (en) Image processing apparatus, image pickup apparatus, method and program
CN116388733B (zh) 一种适用于近红外水果无损检测的光谱预处理方法及装置
CN111489309A (zh) 一种稀疏解混的预处理装置及方法
Cseresnyes et al. Hessian‐based quantitative image analysis of host‐pathogen confrontation assays
CN116304645B (zh) 一种基于模态分解的重叠峰提取的方法及装置
JP2017528705A (ja) 赤外線吸収スペクトルを修正するための方法
Kumar et al. A generative focus measure with application to omnifocus imaging
CN117194902A (zh) 一种光谱测量过程中噪声数据滤除方法
CN116168030B (zh) 极片的缺陷检测方法、装置、电子设备及存储介质
JP2013033006A (ja) 分光情報取得装置、分光情報取得方法及び分光情報取得用プログラム
CN117007184A (zh) 一种光谱的全自动光谱基线校正方法
CN111077128A (zh) 使用相对积分参数的拉曼信号位置校正
US20230058064A1 (en) Film thickness measuring device and film thickness measuring method

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