CN101819168B - 一种光谱分析方法 - Google Patents

一种光谱分析方法 Download PDF

Info

Publication number
CN101819168B
CN101819168B CN2010100339224A CN201010033922A CN101819168B CN 101819168 B CN101819168 B CN 101819168B CN 2010100339224 A CN2010100339224 A CN 2010100339224A CN 201010033922 A CN201010033922 A CN 201010033922A CN 101819168 B CN101819168 B CN 101819168B
Authority
CN
China
Prior art keywords
peak
value
spectrum
variance
coefficient
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.)
Active
Application number
CN2010100339224A
Other languages
English (en)
Other versions
CN101819168A (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.)
Beijing General Research Institute of Mining and Metallurgy
Original Assignee
Beijing General Research Institute of Mining and Metallurgy
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 Beijing General Research Institute of Mining and Metallurgy filed Critical Beijing General Research Institute of Mining and Metallurgy
Priority to CN2010100339224A priority Critical patent/CN101819168B/zh
Publication of CN101819168A publication Critical patent/CN101819168A/zh
Application granted granted Critical
Publication of CN101819168B publication Critical patent/CN101819168B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Spectrometry And Color Measurement (AREA)

Abstract

本发明实施例涉及元素光谱分析领域一种光谱分析方法，利用寻峰算法锁定谱峰的峰位，作为高斯函数的均值初始值；以所述谱峰的峰位对应的峰值作为高斯函数的系数初始值；在全能谱范围内计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差，根据误差最小值对应的谱峰宽度计算高斯函数的方差初始值；在均值、系数及方差取值范围内，计算均值、系数及方差在每种取值下的拟合谱峰面积与原始谱峰面积的误差，所有误差中计算最小的误差所对应的均值、系数及方差作为高斯函数的最终特征参数；使用所述最终特征参数对应的高斯函数进行光谱分析。本发明实施例提高了光谱峰的重现性，为后续的元素定性、定量分析提供更加准确、翔实的数据。

Description

一种光谱分析方法
技术领域
本发明涉及元素光谱分析领域，尤其涉及一种X射线荧光光谱分析方法。
背景技术
荧光分析技术广泛应用于有色、水泥、地质勘探、示踪溯源等行业的元素无损快速分析检测中。其中X射线荧光光谱分析是X射线荧光分析技术的关键，识别谱峰和提取特征参数是分析被测对象元素及其含量的重要步骤，如果被测对象元素组分复杂，相应的基体效应就会严重，荧光光谱谱峰会存在拖尾、重叠、散射本底等干扰，使得谱峰失真，影响后续分析精度。另外，由于统计涨落现象的存在，测量条件不变的情况下，测量同一标样获得的谱峰的重现性也无法保证，这对低含量、颗粒型样品的测量精度影响较大。
因此，现有的X射线荧光光谱分析技术对于元素的定位、定量分析不够准确。
发明内容
本发明实施例提供一种光谱分析方法，还原干扰前的光谱峰，提高光谱峰的重现性，为后续的元素定性、定量分析提供更加准确、翔实的数据。
本发明实施例是通过以下技术方案实现的：
本发明实施例提供一种光谱分析方法，包括：
利用寻峰算法锁定谱峰的峰位，作为高斯函数的均值初始值；
以所述谱峰的峰位对应的峰值作为高斯函数的系数初始值；
在全能谱范围内计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差，根据误差最小值对应的谱峰宽度计算高斯函数的方差初始值；
在均值、系数及方差取值范围内，计算均值、系数及方差在每种取值下的拟合谱峰面积与原始谱峰面积的误差，所有误差中计算最小的误差所对应的均值、系数及方差作为高斯函数的最终特征参数；
使用所述最终特征参数对应的高斯函数进行光谱分析。
所述寻峰算法包括：
执行第一次寻峰，在光谱数组Spe[i]中找出满足如下条件的Spe[i]对应的i，组成初始峰位数组TempPeak[i]，i为光谱道址，取值范围为从0到chan，chan为正整数，包括满足：大于所有峰位下限的加权平均值以及大于相邻的前后峰位的峰值；
执行第二次寻峰，在所述初始峰位数组TempPeak[i]中找出满足如下条件的Spe[i]对应的i作为锁定的峰位，包括满足：(Spe[TempPeak[i]]-(Spe[TempPeak[i]-IniPW]+Spe[TempPeak[i]+IniPW])/2)/IniPW＞1.5，以及满足(TempPeak[i+1]-TempPeak[i])＜IniPW并且(TempPeak[i-1]-TempPeak[i])＞-IniPW，其中IniPW为谱峰宽度，TempPeak[i]为初始峰位。
还包括计算所有峰位下限的加权平均值，具体计算方法包括：
统计所有光谱道范围内满足低于相邻的前后峰位的峰值作为谱峰下限；
计算所有谱峰下限的加权平均值。
所述IniPW初始值为chan/100。
所述在全能谱范围内计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差包括：
在全能谱范围内使用穷举法计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差。
所述在均值、系数及方差取值范围内，计算均值、系数及方差在每种取值下的面积与原始峰面积的误差前，还包括：
确定均值、系数及方差取值范围，其中
均值b取值范围为(b-b/100)＜b＜(b+b/100)，
系数a取值范围为(a-a/10)＜a＜(a+a/10)，
方差c取值范围为(c-c/100)＜c＜(c+c/100)。
所述计算均值、系数及方差在每种取值下的面积与原始峰面积的误差包括：
使用穷举法计算均值、系数及方差在每种取值下的拟合谱峰面积与原始谱峰面积的误差。
由上述本发明实施例提供的技术方案可以看出，本发明实施例采用谱峰拟合的方法对光谱进行还原、解重叠峰、扣除本底等处理，还原干扰前的光谱峰，提高光谱峰的重现性，为后续的元素定性、定量分析提供更加准确、翔实的数据。
附图说明
图1为本发明实施例光谱分析方法流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，可以理解的是，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明一个实施例提供一种光谱分析方法，如图1所示，包括如下步骤：
步骤10：利用寻峰算法锁定谱峰的峰位，作为高斯函数的均值初始值；
本发明一个实施例寻峰算法锁定谱峰的峰位的方法包括：
执行第一次寻峰，在光谱数组Spe[i]中找出满足如下条件的对应的i，Spe[i]组成初始峰位数组TempPeak[i]，i为光谱道址，取值范围为从0到chan，chan为正整数，需要满足的条件包括：Spe[i]值大于所有峰位下限的加权平均值以及大于相邻的前后峰位的峰值；所述所有峰位下限的加权平均值计算方法包括：首先统计所有光谱道范围内满足低于相邻的前后峰位的峰值作为谱峰下限；最后计算所有谱峰下限的加权平均值。
执行第二次寻峰，在所述初始峰位数组TempPeak[i]中找出满足如下条件的Spe[i]对应的i作为锁定的峰位，包括满足：(Spe[TempPeak[i]]-(Spe[TempPeak[i]-IniPW]+Spe[TempPeak[i]+IniPW])/2)/IniPW＞1.5，以及满足(TempPeak[i+1]-TempPeak[i])＜IniPW并且(TempPeak[i-1]-TempPeak[i])＞-IniPW，其中IniPW为谱峰宽度。所述IniPW初始值为chan/100，TempPeak[i]为初始峰位。一种实施例中所述chan取值从256-8192，本发明实施例对此不做限定。
步骤11：以所述谱峰的峰位对应的峰值作为高斯函数的系数初始值；
所述谱峰的峰位对应的峰值即Spe[i]的值。
步骤12：在全能谱范围内计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差，根据误差最小值对应的谱峰宽度计算高斯函数的方差初始值；
所述所有谱峰宽度下的拟合谱峰面积是指谱峰宽度取遍所有可能的取值，每个可能取值对应的高斯拟合函数下的谱峰积分面积。所述原始谱峰面积是指拟合谱峰对应的谱峰宽度范围之内，光谱数据Spe[i]的谱峰面积，即光谱数据求和。
找出拟合谱峰面积和原始谱峰面积的误差最小值，根据计算该误差最小值对应的谱峰宽度可以计算高斯函数的方差初始值，具体计算方法本发明实施例不做限定。
本发明实施例可以在全能谱范围内使用穷举法计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差。
步骤13：在均值、系数及方差取值范围内，计算均值、系数及方差在每种取值下的面积与原始峰面积的误差，所有误差中计算最小的误差所对应的均值、系数及方差作为高斯函数的最终特征参数；
本发明实施例提供的确定均值、系数及方差取值范围的方法包括：
均值b取值范围为(b-b/100)＜b＜(b+b/100)，
系数a取值范围为(a-a/10)＜a＜(a+a/10)，
方差c取值范围为(c-c/100)＜c＜(c+c/100)。
本发明实施例可以使用穷举法计算均值、系数及方差在每种取值下的面积与原始峰面积的误差。
步骤14：使用所述最终特征参数对应的高斯函数进行光谱分析。
本发明上面实施例所述的方法可以用于X射线荧光分析，X射线荧光分析技术广泛应用于物质元素的种类、含量分析，其工作原理是用X射线源激发物质中元素的特征X射线，用探测器接受并转换成电脉冲信号，再通过脉冲分析电路将电脉冲转换为计算机可识别的光谱数据，通过计算机程序分析物质的光谱，计算物质的元素组成和各自含量。因此光谱分析是X射线荧光分析技术的核心部分，而当被测物质的组分较为复杂的情况下，特征X射线的光谱峰会存在重叠、本底等干扰；另外，由于荧光的产生符合统计规律，因此当元素含量较低的情况下，光谱峰的重现性就会较差。本发明实施例所述的光谱分析方法利用荧光谱峰理论上服从高斯分布的特点，采用谱峰拟合的方法，可以对光谱进行还原、解重叠峰、扣除本底等处理，还原干扰前的光谱峰，提高光谱峰的重现性，为后续的元素定性、定量分析提供更加准确、翔实的数据。
为进一步理解本发明，下面以具体实施例对本发明提供的光谱分析方法进行详细描述，包括如下步骤：
步骤20：将光谱转换为光谱数组Spe[i]，下标i为光谱道址，取值范围从0到chan，对应的值为峰位计数；
步骤21：在光谱道范围内进行寻峰，首先确定谱峰下限，统计所有光谱道范围内满足spe[i-1]＞Spe[i]并且Spe[i+1]＞Spe[i]的Spe[i]，求加权平均值AveSpe；
步骤22：统计所有光谱道范围内满足Spe[i]＞Avespe，以及spe[i-1]＜Spe[i]并且spe[i+1]＜Spe[i]的Spe[i]对应的i组成初始峰位数组TempPeak[i]；
步骤23：在TempPeak[i]范围内进行二次寻峰，设定初始谱峰宽度；
1024道初始谱峰宽度为IniPW＝10，2048道初始谱峰宽度IniPW＝20，即IniPW＝chan/100；
步骤24：统计TempPeak[i]范围内满足(Spe[TempPeak[i]]-(Spe[TempPeak[i]-IniPW]+Spe[TempPeak[i]+IniPW])/2)/IniPW＞1.5，以及(TempPeak[i+1]-TempPeak[i])＜IniPW并且(TempPeak[i-1]-TempPeak[i])＞-IniPW的对应的i，得到光谱道范围内的所有准谱峰峰位Peak[i]；
步骤25：采用拟合方法计算各谱峰的特征参数。
以单个谱峰为例：
高斯函数公式为
Figure GDA0000074360340000062
其中a为系数，b为均值，c为方差，与谱峰宽度呈线性关系；
以谱峰峰位i作为高斯函数的均值初始值ini_b，峰位对应峰值Peak[i]作为系数初始值ini_a；
在全能谱范围内，用穷举法计算方差c，即c从1到谱峰道址的最大值，计算高斯函数与原始能谱之间的误差平均值ave_err[i]，也就是计算所有谱峰宽度下拟合谱峰面积与原始谱峰面积的误差ave_err[i]；
寻找ave_err[i]中最小值，对应的c值即为高斯函数的偏差初始值ini_c；
设定ini_a，ini_b，ini_c的范围分别为[ini_a-ini_a/10，ini_a+ini_a/10]，[ini_b-chan/100，ini_b+chan/100]，[ini_c-chan/100，ini_c+chan/100]，数组形式为ini_a_arr[i]，ini_b_arr[i]，ini_c_arr[i]；
穷举法统计ni_a_arr[i]，ini_b_arr[i]，ini_c_arr[i]所有的可能性，计算高斯函数与原始能谱之间的误差平均值ave_err；
寻找ave_err最小值对应的一组ni_a_arr[i]，ini_b_arr[j]，ini_c_arr[k]，即为高斯函数的a、b、c值。
以该a、b、c值确定的高斯函数进行光谱分析。
综上所述，本发明实施例采用谱峰拟合的方法对光谱进行还原、解重叠峰、扣除本底等处理，还原干扰前的光谱峰，提高光谱峰的重现性，为后续的元素定性、定量分析提供更加准确、翔实的数据。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

Claims (6)

1.一种光谱分析方法，其特征在于，包括：
利用寻峰算法锁定谱峰的峰位，作为高斯函数的均值初始值；
以所述谱峰的峰位对应的峰值作为高斯函数的系数初始值；
在全能谱范围内计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差，根据误差最小值对应的谱峰宽度计算高斯函数的方差初始值；
在均值、系数及方差取值范围内，计算均值、系数及方差在每种取值下的拟合谱峰面积与原始谱峰面积的误差，所有误差中计算最小的误差所对应的均值、系数及方差作为高斯函数的最终特征参数；
使用所述最终特征参数对应的高斯函数进行光谱分析；
所述寻峰算法包括：
执行第一次寻峰，在光谱数组Spe[i]中找出满足如下条件的Spe[i]对应的i，组成初始峰位数组TempPeak[i]，i为光谱道址，取值范围为从0到chan，chan为正整数，包括满足：大于所有峰位下限的加权平均值以及大于相邻的前后峰位的峰值；
执行第二次寻峰，在所述初始峰位数组TempPeak[i]中找出满足如下条件的Spe[i]对应的i作为锁定的峰位，包括满足：(Spe[TempPeak[i]]-(Spe[TempPeak[i]-IniPW]+Spe[TempPeak[i]+IniPW])/2)/IniPW＞1.5，以及满足(TempPeak[i+1]-TempPeak[i])＜IniPW并且(TempPeak[i-1]-TempPeak[i])＞-IniPW，其中IniPW为谱峰宽度，TempPeak[i]为初始峰位。
2.如权利要求1所述的方法，其特征在于，还包括计算所有峰位下限的加权平均值，具体计算方法包括：
统计所有光谱道范围内满足低于相邻的前后峰位的峰值作为谱峰下限；
计算所有谱峰下限的加权平均值。
3.如权利要求1所述的方法，其特征在于，所述IniPW初始值为chan/100。
4.如权利要求1所述的方法，其特征在于，所述在全能谱范围内计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差包括：
在全能谱范围内使用穷举法计算所有谱峰宽度下的拟合谱峰面积和原始谱峰面积的误差。
5.如权利要求1所述的方法，其特征在于，所述在均值、系数及方差取值范围内，计算均值、系数及方差在每种取值下的面积与原始峰面积的误差前，还包括：
确定均值、系数及方差取值范围，其中
均值b取值范围为(b-b/100)＜b＜(b+b/100)，
系数a取值范围为(a-a/10)＜a＜(a+a/10)，
方差c取值范围为(c-c/100)＜c＜(c+c/100)。
6.如权利要求1所述的方法，其特征在于，所述计算均值、系数及方差在每种取值下的面积与原始峰面积的误差包括：
使用穷举法计算均值、系数及方差在每种取值下的拟合谱峰面积与原始谱峰面积的误差。
CN2010100339224A 2010-01-06 2010-01-06 一种光谱分析方法 Active CN101819168B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010100339224A CN101819168B (zh) 2010-01-06 2010-01-06 一种光谱分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010100339224A CN101819168B (zh) 2010-01-06 2010-01-06 一种光谱分析方法

Publications (2)

Publication Number Publication Date
CN101819168A CN101819168A (zh) 2010-09-01
CN101819168B true CN101819168B (zh) 2012-06-13

Family

ID=42654343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010100339224A Active CN101819168B (zh) 2010-01-06 2010-01-06 一种光谱分析方法

Country Status (1)

Country Link
CN (1) CN101819168B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608649B (zh) * 2012-03-02 2013-11-27 成都理工大学 统计分布式γ或X射线能谱解谱方法
CN103105428B (zh) * 2013-01-14 2014-11-05 深圳市清时捷科技有限公司 一种判断水中重金属是否超限的方法
CN104316426B (zh) * 2014-10-24 2017-08-04 北京矿冶研究总院 一种矿石可用性分析的方法
US10634653B2 (en) * 2016-11-28 2020-04-28 Dionex Corporation Systems, methods and devices for width-based analysis of peak traces
CN108918474B (zh) * 2018-04-26 2022-07-08 安徽师范大学 用于降低激光诱导击穿光谱自吸收效应的基体稀释方法
CN110515115B (zh) * 2018-05-22 2022-11-25 中国辐射防护研究院 一种用于伽马谱峰面积求解的评估方法
CN109283169A (zh) * 2018-11-22 2019-01-29 深圳市雷泛科技有限公司 一种鲁棒的拉曼光谱峰识别方法
CN113835114B (zh) * 2021-08-25 2024-04-26 吉林大学 一种紧凑式高能伽玛射线反符合叠层探测器
CN114414048B (zh) * 2021-12-23 2023-05-12 电子科技大学 一种提高调制光谱测量精度的装置与方法
CN114965348B (zh) * 2022-07-27 2022-11-15 浙江数翰科技有限公司 基于污水检测的光谱解析方法和系统
CN115656238B (zh) * 2022-10-17 2023-05-12 中国科学院高能物理研究所 一种微区xrf元素分析与多维成像方法及系统

Also Published As

Publication number Publication date
CN101819168A (zh) 2010-09-01

Similar Documents

Publication Publication Date Title
CN101819168B (zh) 一种光谱分析方法
CN101788507B (zh) 一种光谱分析方法
Coetzee et al. Intraregional classification of wine via ICP-MS elemental fingerprinting
US11313223B2 (en) Gamma ray spectrum unfolding method for elemental capture spectroscopy logging and device therefor
CN103913764A (zh) 一种基于高斯响应矩阵的NaI(TI)闪烁探测器γ能谱高分辨反演解析过程及方法
Ehlert et al. X-ray bright active galactic nuclei in massive galaxy clusters–I. Number counts and spatial distribution
CN105758836A (zh) 一种基于面积法的拉曼光谱即时定量分析方法
Valentino et al. Predicting emission line fluxes and number counts of distant galaxies for cosmological surveys
CN112231621B (zh) 基于BP-adaboost降低元素检出限的方法
Pickering et al. Precise laboratory wavelengths of the Mg I and Mg II resonance transitions at 2853, 2803 and 2796 Å
Zeiri et al. Determination of uranium isotopic ratio by ICP-OES using optimal sensitivity position analysis
CN105572067A (zh) 基于光谱分析的烟气浓度测量方法
CN110515115B (zh) 一种用于伽马谱峰面积求解的评估方法
CN103592287B (zh) 通过元素的测量强度和浓度计算被干扰元素浓度的方法
CN105138800A (zh) 基于分段拟合数据处理算法的荧光光谱数据噪声滤波方法
CN104504273A (zh) 一种提高γ射线峰探测灵敏度的能谱处理方法
Zhang et al. Effect of Soil Moisture and Particle Size on Soil Total Phosphorus Estimation by Near-Infrared Spectroscopy.
Dronov et al. Strontium carbonate precipitation as a sample preparation technique for isotope ratio analysis of Sr in mineral water and wine by quadrupole‐based inductively coupled plasma mass spectrometry
CN105865625B (zh) 多维信息一体化获取相机的光谱定标方法
Miller et al. Matrix effects on the Raman analytical lines of oxyanions
CN103592226B (zh) 通过元素的测量强度计算被干扰元素浓度的方法
Trippe AGN broad line regions scale with bolometric luminosity
JP3766883B2 (ja) 元素分析方法
Kotulla et al. Impact of subsolar metallicities on photometric redshifts
Tamburro et al. Cosmic ray spectrum, composition, and anisotropy measured with IceCube

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