专利名称：一种放射性能谱识别方法
技术领域：
本发明涉及一种放射性能谱识别方法。
背景技术：
在进行放射性的能谱测量中，常常要根据探测器采集到的大量随机数据，分析这些随机信号由什么放射源产生。一般的放射性能谱解析方法是对所测的能谱在进行了平滑滤波、扣除本底、寻峰、计算净峰面积等分析后，根据能量、峰形和效率刻度的结果，由峰位所对应的能量识别出放射源含有哪些放射性核素，进而由净峰面积计算出这些放射性核素的活度。在某些应用场合，由于核信息保护的原因，不需要得到核素活度的精确结果，只需识别所测对象为何种核材料。为此，曾有研究者采用能谱信息作为神经网络的训练样本，并用训练后的神经网络对这些核材料进行识别。这些方法存在难以训练、或者不易收敛、或者收敛到局部极小等缺点，导致能谱的识别不能进行或者识别准确率降低，故需要找出一种既能快速训练又能避免不能收敛或局部收敛的方法，使能谱的识别能顺利进行且能保证更高的识别准确率。
发明内容
本发明的目的在于公开一种放射性能谱识别方法。该方法克服了现有能谱识别方法的不足，具有训练速度快、收敛准确、识别准确率高等特点，是一种进行放射性能谱识别的有效方法。本发明是通过以下技术方案实现的，本发明的具体步骤如下①对放射性测量中已测得的能谱进行训练，按如下A E步骤A根据需要按放射源类别对已测得的多个能谱进行分类，B采用小波方法或多项式方法对能谱进行平滑滤波，C求取能谱的GMM模型(Gaussian mixture model)，并将GMM模型中各高斯函数的权值作为降维后的数据，D将降维后的能谱数据进行小波包分解，并将各子频带信号的能量进行归一化处理，提取归一化数据作为特征向量。E将特征向量作为样本，采用期望最大化法(Expectation Maximization，简写为 EM)对这些样本进行迭代运算，得到各类特征向量的GMM模型；②将待识别能谱进行决策分类，按如下A E步骤A采用小波方法或多项式方法对待识别能谱进行平滑滤波，B求取待识别能谱的GMM模型(Gaussian mixture model)，并将GMM模型中各高斯函数的权值作为降维后的数据，C将降维后的待识别能谱数据进行小波包分解，并将各子频带信号的能量进行归一化处理，提取归一化数据作为特征向量，D将待识别能谱的特征向量作为多维随机数，并计算其分属于各类GMM模型的类条件概率密度，E最后按贝叶斯决策分类。本发明的有益效果是本发明在训练阶段先采用小波方法或多项式方法对能谱进行平滑滤波，以消除干扰信号；接着用GMM模型对能谱作降维处理，以保证能谱的统计特性在低维空间仍具有原始能谱信号的绝大部分信息；然后，在多频带信号空间提取归一化特征向量并作为样本，保证了采用低维样本代表能谱特征的唯一性，同时也保证了样本与测量时间的无关性；最后， 采用期望最大化法对这些样本进行迭代运算，得到各类特征向量的GMM模型，这样既保证了准确收敛又可通过合理调整GMM模型的高斯函数个数以满足各种不同场合下识别的需要。另外，本发明在识别时，采用贝叶斯决策分类，保证了识别准确率在统计意义下的最优。 总之，本发明是一种既能快速训练又能避免不能收敛或局部收敛的方法，使能谱的识别能顺利进行且能保证更高的识别准确率。
图1为本发明方法的流程图；
具体实施例方式为使发明的目的、技术方案及优点更加清楚明白，以下参照附图并举实施例，对本发明作进一步详细说明。由于核信息保护的原因，不需要得到核素活度精确结果，只需识别所测对象为何种核材料，针对这一情况，本发明提供了一种放射性能谱识别方法。图1显示了本发明所述识别方法的流程。本发明的流程如图1所示，具体包括如下的训练步骤①和识别步骤②①对放射性测量中已测得的能谱进行训练，按如下A E步骤A根据需要按放射源类别——如，标准源或样品的种类及批次等——对已测得的多个能谱进行分类，能谱为在不同测量距离、不同测量时间及其它不同外界条件下所测得；B采用小波方法或多项式方法对能谱进行平滑滤波；C求取能谱的GMM模型(Gaussian mixture model)，即表示为多个高斯分布函数的线性和Ρ{χ,θ) = ^α Ρι{χ·Α)(1)式(1)中M为高斯分布密度函数的个数，应由能谱的形状、光滑程度、道址数N等而定；a1; . . .，％是各高斯分布密度函数的权重，即各高斯分布函数在概率密度函数中所占
M
的比重，且Σ & = 1 彡0，(i = 1，. . .，M)，由能谱数据和函数个数M而定；Pi (χ)是第i个 /=1
高斯分布密度函数，其均值Pi由高斯分布密度函数的个数M和道址数确定；Pi(x)的方差为Oi2; Qi是未知参数^和Qi2的向量表示，S卩祝二細，。2)，密度函数？&，Θ,)的形
权利要求
1.一种放射性能谱识别方法，其特征在于，具体步骤如下①对放射性测量中已测得的能谱进行训练；②将待识别能谱进行决策分类。
2.根据权利要求1所述的一种放射性能谱识别方法，其特征是，所述①中对放射性测量中已测得的能谱进行训练，包含如下步骤A对放射性测量中已测得的能谱进行滤波、降维和分类处理，B提取处理后能谱的特征向量，C将特征向量作为样本，并通过训练得到各类特征向量的GMM模型。
3.根据权利要求2所述的一种放射性能谱识别方法，其特征是，所述A中滤波是指采用小波方法或多项式方法对能谱进行平滑滤波；降维是指求取能谱的GMM模型(Gaussian mixture model)，并将GMM模型中各高斯函数的权值作为降维后的数据；分类是指按放射源进行分类。
4.根据权利要求2所述的一种放射性能谱识别方法，其特征是，所述B中提取特征向量，是指将权利要求3降维后的数据进行小波包分解，并将各子频带信号的能量进行归一化处理，提取归一化数据作为特征向量。
5.根据权利要求2所述的一种放射性能谱识别方法，其特征是，所述C中通过训练得到各类特征向量的GMM模型，是指将权利要求4提取的特征向量作为多维随机数，采用期望最大化法(Expectation Maximization，简写为EM)对这些随机数进行迭代运算，得到各类特征向量的GMM模型。
6.根据权利要求1所述的一种放射性能谱识别方法，其特征是，所述②中将待识别能谱进行决策分类，包含如下步骤a采用小波方法或多项式方法对待识别能谱进行平滑滤波，并求取待识别能谱的特征向量，b按GMM模型求取特征向量的类条件概率密度，c然后按贝叶斯决策分类。
7.根据权利要求6所述的一种放射性能谱识别方法，其特征是，所述a中求取待识别能谱的特征向量，是指将待识别能谱按权利要求3降维后进行小波包分解，并将各子频带信号的能量进行归一化处理，提取归一化数据作为特征向量。
8.根据权利要求6所述的一种放射性能谱识别方法，其特征是，所述b中按GMM模型求取特征向量的类条件概率密度，是指将权利要求7提取的待识别能谱特征向量作为待分类的多维随机数，并计算其分属于各类GMM模型的类条件概率密度。
全文摘要
本发明公开了一种放射性能谱识别方法。对放射性测量中已测得的能谱进行滤波、降维和分类处理，提取其特征向量作为样本，并通过训练得到各类特征向量的GMM模型；求取待识别能谱的特征向量，并按GMM模型求取特征向量的类条件概率密度，然后决策分类，即完成能谱识别。本发明对能谱的识别准确率高，是进行放射性能谱识别的有效方法。
文档编号G06F19/00GK102313897SQ201010212350
公开日2012年1月11日 申请日期2010年6月29日 优先权日2010年6月29日
发明者丁卫撑, 刘念聪, 刘易, 周伟, 方方, 王敏, 王超, 阎萍, 黄洪全, 龚迪琛 申请人:成都理工大学