发明名称 基于小波分析的激光感生光谱数据处理方法
摘要 本发明涉及光谱分析技术领域,具体是一种能提高光谱分析精度及效率的基于小波分析的激光感生光谱数据处理方法。解决了目前激光感生光谱分析结果受无效测量数据影响等问题,包括:1)有效测量模式类模板G<sup>m</sup>获取步骤,在对待测物质标准样品定标时,提取有效测量模式类模板,以有效测量模式类模板对应的测量数据进行定标运算,获取定标参数;2)被测样品实际光谱测量数据处理步骤,实际测量时,将实测光谱数据与有效测量模式类模板进行比对,保留有效光谱测量数据,以定标参数对有效测量光谱数据进行计算,获得被测样品元素含量信息。本发明有效提高测量精度,以有限次激光激发测量获取所需个数的有效测量数据,有效提高激光测量系统的寿命。
申请公布号 CN101666746A 申请公布日期 2010.03.10
申请号 CN200910075578.2 申请日期 2009.09.30
申请人 太原市海通自动化技术有限公司;山西大学 发明人 贾锁堂;张生俊;尹王保;阎高伟;王红兵;李平柱;罗振红;王学钦;张雷
分类号 G01N21/63(2006.01)I 主分类号 G01N21/63(2006.01)I
代理机构 山西太原科卫专利事务所 代理人 朱 源;骆 洋
主权项 1、一种基于小波分析的激光感生光谱数据处理方法,其特征在于包括如下步骤:1)有效测量模式类模板G<sup>m</sup>获取步骤;a、对待测物质标准样品组中的n个标准样品分别进行m次激光感生光谱测量,并顺序记录与各标准样品对应的各次光谱测量数据G<sub>i,j</sub>,构成标准样品光谱测量数据集合G={G<sub>i,j</sub>},其中,i=1,2,…,n,j=1,2,…,m,G<sub>i,j</sub>表示对标准样品i进行第j次激光感生光谱测量得到的光谱测量数据,且光谱测量数据G<sub>i,j</sub>以序列方式表示:G<sub>i,j</sub>(k)=[X<sub>1</sub>,X<sub>2</sub>,…X<sub>k</sub>,…X<sub>N</sub>],N为序列长度;b、对标准样品光谱测量数据集合G中各光谱测量数据序列G<sub>i,j</sub>(k)进行L尺度一维离散平稳小波分解,4≤L≤8,得到分别与各光谱测量数据序列G<sub>i,j</sub>(k)对应的高频分解系数<img file="A2009100755780002C1.GIF" wi="301" he="65" />和低频分解系数<img file="A2009100755780002C2.GIF" wi="327" he="58" />c、以与各光谱测量数据序列G<sub>i,j</sub>(k)分别对应的低频小波分解系数W<sub>i,j</sub><sup>a</sup>进行光谱重构,获得与各光谱测量数据序列G<sub>i,j</sub>(k)一一对应的特征背景光谱数据G<sub>i,j</sub><sup>b</sup>,构成特征背景光谱数据集合<img file="A2009100755780002C3.GIF" wi="253" he="65" />特征背景光谱数据G<sub>i,j</sub><sup>b</sup>同样以序列方式表示:<img file="A2009100755780002C4.GIF" wi="713" he="66" />且序列长度与光谱测量数据序列G<sub>i,j</sub>(k)的序列长度相同;d、对特征背景光谱数据集合G<sup>b</sup>中的特征背景光谱数据G<sub>i,j</sub><sup>b</sup>进行聚类分析,将特征背景光谱数据集合G<sup>b</sup>划分为若干个模式类子集G<sub>h</sub><sup>b</sup>,即<img file="A2009100755780002C5.GIF" wi="616" he="61" />其中,h=1,2,…,H,H为将特征背景光谱数据集合按照聚类分析划分得到的模式类子集个数;按照光谱测量数据G<sub>i,j</sub>与特征背景光谱数据G<sub>i,j</sub><sup>b</sup>的一一对应关系、以及特征背景光谱数据集合G<sup>b</sup>的划分,将标准样品光谱测量数据集合G划分为与特征背景光谱数据集合G<sup>b</sup>模式类子集G<sub>h</sub><sup>b</sup>一一对应的若干个模式类子集G<sub>h</sub>,即G={G<sub>1</sub>,G<sub>2</sub>,…,G<sub>h</sub>,…,G<sub>H</sub>};e、对标准样品光谱测量数据集合G中每一模式类子集G<sub>h</sub>包含的光谱测量数据G<sub>i,j</sub>进行元素含量定标运算,获得与各模式类子集G<sub>h</sub>一一对应的定标参数β<sub>h</sub>和定标运算结果;f、选择定标运算结果与待测物质标准样品的标准值相差最小的模式类子集G<sub>h</sub>所具有的模式作为有效测量模式,与该模式类子集G<sub>h</sub>对应的定标参数β<sub>h</sub>为对被测样品实测时计算元素含量采用的定标参数;g、提取与该模式类子集G<sub>h</sub>对应的特征背景光谱数据模式类子集G<sub>h</sub><sup>b</sup>,选择特征背景光谱数据模式类子集G<sub>h</sub><sup>b</sup>中所有特征背景光谱数据序列G<sub>i,j</sub><sup>b</sup>(k)中位置k的最大值作为有效测量模式类模板序列位置k的上限值G<sub>h</sub><sup>m</sup>(k),位置k的最小值作为有效测量模式类模板序列位置k的下限值G<sub>l</sub><sup>m</sup>(k),即<img file="A2009100755780003C1.GIF" wi="453" he="99" /><img file="A2009100755780003C2.GIF" wi="435" he="97" />形成有效测量模式类模板<img file="A2009100755780003C3.GIF" wi="445" he="52" />其中,1≤k≤N,E为特征背景光谱数据模式类子集G<sub>h</sub><sup>b</sup>中特征背景光谱数据序列G<sub>i,j</sub><sup>b</sup>(k)的个数,N为特征背景光谱数据序列G<sub>i,j</sub><sup>b</sup>(k)的序列长度;2)、被测样品实际光谱测量数据处理步骤;h、对被测样品进行单次激光感生光谱测量,并记录该次光谱测量数据G<sub>j</sub>,j=1,2,…;且以序列方式表示:G<sub>j</sub>(k)=[X<sub>1</sub>,X<sub>2</sub>,…,X<sub>k</sub>,…,X<sub>N</sub>],N为序列长度;i、对激光感生光谱数据序列G<sub>j</sub>(k)进行L尺度的一维离散平稳小波分解,得到与激光感生光谱数据序列G<sub>j</sub>(k)对应的高频分解系数<img file="A2009100755780003C4.GIF" wi="297" he="79" />和低频分解系数<img file="A2009100755780003C5.GIF" wi="320" he="80" />其中,L值与步骤b中L值相同;j、以步骤g得到的低频小波分解系数W<sub>j</sub><sup>a</sup>进行光谱重构,获得与激光感生光谱数据序列G<sub>j</sub>(k)对应的特征背景光谱数据G<sub>j</sub><sup>b</sup>;特征背景光谱数据G<sub>j</sub><sup>b</sup>同样以序列方式表示:<img file="A2009100755780004C1.GIF" wi="692" he="75" />且序列长度与光谱测量数据序列G<sub>j</sub>(k)的序列长度相同;k、判断特征背景光谱数据序列G<sub>j</sub><sup>b</sup>(k)是否与有效测量模式类模板<img file="A2009100755780004C2.GIF" wi="413" he="52" />匹配,即<img file="A2009100755780004C3.GIF" wi="472" he="64" />是否成立,如成立,该特征背景谱G<sub>j</sub><sup>b</sup>(k)属于有效测量模式类,则与该特征背景光谱数据G<sub>j</sub><sup>b</sup>的测量光谱数据G<sub>j</sub>有效;l、对被测样品重复执行步骤h至步骤k,直到获得三次以上的有效测量光谱数据G<sub>j</sub>;m、对所获得的有效测量光谱数据G<sub>j</sub>按照步骤f中得到的与有效测量模式类子集G<sub>h</sub>对应的定标参数β<sub>h</sub>进行元素含量计算,以计算结果的平均值作为被测样品的元素含量分析结果。
地址 030006山西省太原市高新区长治路266号科祥大厦601室