多道γ能谱分析软件中寻峰算法比较总结.docx

1、多道能谱分析软件中寻峰算法比较总结自动寻峰由于谱结构的复杂和统计涨落的影响，从谱中正确地找到全部存在的峰是比较困难的。 尤其是找到位于很高本底上的弱峰，分辨出相互靠得很近的重峰更为困难。谱分析对寻峰方法的基本要求如下：(1)比较高的重峰分辨能力。能确定相互距离很近的峰的峰位。(2)能识别弱峰，特别是位于高本底上的弱峰。(3)假峰出现的几率要小。(4)不仅能计算出峰位的整数道址，还能计算出峰位的精确值，某些情况下要求峰位的误 差小于0.2道。很多作者对寻峰方法进行了研究，提出了很多有效的寻峰方法。目的：判断有没有峰存在确定峰位(高斯分布的数学期望)，以便把峰位对应的道址，转换成能量确定峰边界 为

2、计算峰面积服务(峰边界道的确定，直接影响峰面积的计算) 分为两个步骤：谱变换和峰判定要求：支持手动/自动寻峰，参数输入，同时计算并显示峰半高宽、精确峰位、峰宽等信 息，能够区分康普顿边沿和假峰感兴区内寻峰人工设置感兴趣大小，然后在感兴区内采用简单方法寻峰重点研究：对感兴区内的弱峰寻峰、重峰的分解对于一个单峰区，当峰形在峰位两侧比较对称时，可以由峰的 FWHM计算峰区的左、右 边界道址。峰区的宽度取为3FWHM,FWHM的值可以根据峰位mp由测量系统的FWHM刻度公式计算。由于峰形对称，左、右边界道和峰位的距离都是 1.5FWHNMmiL =INT(mp -1.5FWHM 0.5) mR =IN

3、T(mp 1.5FWHM 0.5)式中mp是峰位，INT的含义是取整数。对于存在有低能尾部的峰，其峰形函数描述(参见图)。y m =h exp ( mp r / 22 mmp 一 j2ym =HEXPJ(2m-2mp J)/2； , m RO确定峰是否有意义R为峰谷比，R0为设定值 (经验值)NO为净峰幅度与基底之和Nb为基底计数int CMmcaView:SearPeakCompare(i nt Begin ch, int En dch, int m, float k)高斯乘积函数找峰法(可靠性差，不建议采用)描述谱峰形状的函数主要是高斯函数 G(i)= A exp-(i-i0)2/2二2】

4、则由相邻的数据点定义壬2兀貯一个新的函数(第一高斯乘积函数，只与 FWHM =2.3556二有关):G(i)G(i +m 1) 11.092mPm二(i-2)G(i m厂和m2m是步长(用道表示)，是高斯乘积函数的阶数，则 Pm(i)称为第m阶高斯乘积函数。找峰的灵敏度与m有关，随m的增加灵敏度提高。为避免基线参数的影响，最好扣除本底后，再应用高斯乘积函数找峰。考虑统计涨落的影响，把判断无峰存在的 1变为一个 单位带”即峰的判断为:Pm(i)T兰1 土k/pyi 无峰(1+k/jyr)有峰k _3峰位的确定：由Pm(i)过1的两点求平均来确定；峰边界的确定:单位带”下限的两个最端点；半高宽的确

5、定：函数Pm(i)在 “ 1上的截距；组合峰的确定:在乘积函数的两个峰之间没有处于 带内”的乘积函数值导数法(一阶、二阶、三阶)1 my厂 Cjyi.jN m j 二ANm为规范化常数，Cj平滑的变换系数。3次多项式5点光滑一阶导数公式：(可以采用)y丄卜2-8%丄lyiiy2)峰位确定：一阶导数值由正变负=0处；峰边界确定：一阶导数12由负变正=0处CalculateDiffere ntial(0, size, m, differ);for (int j = m; j = size-m; j+)for(i nt i=1;i0&differj-imaxtemp) maxtemp=differj

6、-i; nmax=j-i;if(differj+i) 0&differj+i0.8* fwhm & (nmin-nmax)3* fwhm )/FWHM参数根据仪器能量分辨率可人工确定， fwhm20peakpositionp+=j+0.5;保持峰位对应的道址5点光滑二阶导数公式(软件中推荐采用)1yi =7(2y24-2yi - yi 1 2yi 2)/7点二阶导数 (5*(cou ntsdataj-3+cou ntsdataj+3)-3*(cou ntsdataj-1+cou ntsdataj+1)-4*cou ntsdataj)/42;软件中推荐采用11点以上的公式峰位确定：二阶导数最小值

7、对应的道址；峰边界确定：二阶导数正极大值点for (int j = m; j = size-m; j+)/m30int maxtemp=-0.5,mi ntemp=-0.5;If(differj -0.05)for(i nt i=1;imaxtemp) maxtemp=differj-i; nmax=j-i; if(differj+imintemp) mintemp=differj+i; nmin=j+i;if (nmin-nmax)0.8* fwhm & (nmin-nmax) Threshold试验：系列1为处理后的原始能谱，系列 2为5点一阶导数，系列 3为5点二阶导数，系列 4为对称零

8、面积法寻峰只要选择好合适的寻峰阈值，足以满足准确寻找到全能峰，并剔除假峰（如康普顿边沿，反散射峰）5点光滑三阶导数公式判定各感兴区是单峰还是重峰，” 1y =2（ 一 y，+2% 斗一 2% 十+% 七）峰位确定：三阶导数由负变正=0处；峰边界确定：三阶导数由正变负=0处判定峰是否有意义 0.8FWHM N 3FWHM峰高判定条件| ym |max 亠TRH ympe / :二这个公式就是在一阶导数法寻峰程序中实际应用的峰高判定条件。CalculateDiffere ntial(Begi nch, En dch, m, differ);int CMmcaView:SearPeakDiffere

9、 ntial( int Begi nch, int En dch, int fwhm, int differ, i nt m)intn 1=0, differE ndch-Beg in ch+1, nm ax=O, nmin=O, maxtemp, min temp,temp;maxtemp=differO; min temp=differ0;for (int j = 1; j = En dch-Begi nch; j+)temp=differj-1;if(_copysign(temp,differj)!=differj-1 & differj0) n1=j+Beginch ;if(diffe

10、rjmaxtemp) maxtemp=differj; nm ax=j+Begi nch;if (nmin-nm ax)0.8*fwhm & (n mi n-nm ax)3*fwhm) return n1;else return (0);对称零面积法(推荐自动寻峰中采用，可探测弱峰和重峰)面积为零的“窗”函数与实验谱数据进行褶积变换，且要求“窗”函数为对称函数。对线性基底的褶积变换将为零，只有存在峰的地方不为零。y = CjyiT 二 Cj = 0 Cj - C jj -_m j -_m;2 1 m k2匹配滤波器法(类峰形函数)Cj二e)p -y - exo一一 2cr 2m+1k=m 2c

11、rm z Cjdata叶峰判定准则R也=-.f瓦 (m 军送 Cjdatai书j =-m 丿2m+1为变换宽度，二二FWHM .2.3556为峰宽参数，若变换后的y和其均方根误差的比值超过预先给定的寻峰阈值(f)，则认为找到了一个峰。峰位的确定：Ri的正极值对应的道址；峰边界的确定： Ri的正峰两边相邻的两个极小值之间的距离可以作为峰的宽度信息；半宽度：两过零截距。CalculateArea(0, size, m, fwhm, area, R);for (int j = m; j 0&Rjfh)for(i nt i=1;i 0&areaj-i 0&areaj+i0.6*fwhm & (nmin

12、 2-nmin 1)=2*fwhm)peakpositionp+=j+0.5;保持峰位对应的道址协方差法(曲线拟合寻峰，计算机寻峰中采用，可分辨重峰，比较好的寻峰方法，但计算较为复杂，运算速度较慢)1975年H.P.BLOK等提出了一种新的寻峰方法，称为协方差法。用一个峰形函数与实验谱数据逐段拟合(一个高斯形函数与实验谱 yi的协方差)yi j 二 y；Cj b，Cj 为峰形 /高斯函数 Cj 二 EXP-2.773(j/H2)H 为峰 FWHM，yi 为拟合峰高，bi 为本底常数(在峰区内假定不变)-m ： j ： mL gj L gjCjyi4j gjCj gj yHj j -_m j -

13、_m j -_m j -_myi =2 2 g gjCj -C gjCj)j m j j e E gj 为 gjCjyi 书一为 gjCj 为 gjyi+i 用 r y i j-mj-R j-q j -n区 gj |艺 gj gjCj (gjCj)j j j jCj通常为纯峰形函数高斯函数： Cj = exp -41 n 2 （丄）2 1, H为峰的FWHMIL H参数选择：H的取值最好与实验谱峰的半宽度接近，2m+1 般取2H左右最好，f 般取25峰位确定：当Ri为极大值对应的道址；峰边界确疋： Ri为负极大值处对应的道址为了更好地分辨出落在一个强峰肩部上的弱峰，可以在一个峰的左半部分和右半

14、部 分别计算R值，寻找相互靠得很近的组分峰。线性拟合寻峰方法（适合于在峰区内分辨重峰）吸取匹配滤波器方法的优点，同时用一阶导数法和线性拟合双重峰的技术来提高分辨重 峰的能力，形成了一种新的寻峰方法，称为线性拟合寻峰方法。Dec onvo luti on methodFirst the backgro und is removed (if desired), the n Markov spectrum is calculated (if desired), the n the resp on se fun cti on is gen erated accord ing to give n sig

15、ma and dec onvo lutio n is carried out.可以提供多种算法，方便自行选择总结1对于弱峰，数据光滑前，高斯乘积函数法和协方差法不能使用，若先光滑再找峰，又 容易影响重叠峰的分辨；而导数法和对称零面积变换法，无论峰的统计质量如何，均可 使用。2.从统计假峰及高基底的抑制能力及重峰的分辨能力来看，一、三阶导数法和对称零面 积变换法是较好的。对于一、三阶导数法，可先用适当多数据点的一阶导数法找峰，选 取适当的灵敏度常数，以抑制假峰；然后用少点的三阶导数法（或用一阶导数法重复三 次）检查是否有漏峰和重峰。对称零面积变换法同理。3. 从高基底的抑制能力和弱峰识别的准确度

16、来看，对称零面积变换法最好。 （在计算机 自动找峰程序中，最好采用对称零面积变换法。）参考资料4对找到的峰进行净面积判定是降低假峰出现几率的有效方法。 当峰的净面积比峰的总面积（峰的净面积和本底面积之和）的标准偏差大若干倍时，才确认该峰是一个真峰，否 则认为它是假峰，予以剔除。峰的判弃主要是利用峰面积来进行判定真假峰。对于给定的灵敏因子S,若峰的净面积为NetAREA,峰的宽度为Width。这些参数满足下式认为峰有意义，应保留，否则将找到的此峰丢弃。此式为：c . NetAREA WidthS ： 1 2（AREA/Width）S越大灵敏度越高，一般情况下 S=3。参考文献SPECTRAN-F Versio n 2 List ings V olume.4 Com mon Subrouti nes and Data structures, CANBERRA In dustries, I nc. 1980 能谱的数据处理（原文）.docM.A.Mariscotti, Nucl. Instr. Meth. 50 （1967） 309Mariscotti Algorithm modified by Routti and Prussin:J.T.Routti and S.G.Prussin, Nucl. Instr. Meth. 72 （1969） 125