一种闪烁脉冲信息复原的方法及系统技术方案

一种闪烁脉冲信息复原的方法，其包括步骤：获取低计数下的非堆积符合单事件的闪烁脉冲数据库，再对各路单事件的闪烁脉冲数据库建立闪烁脉冲的噪声模型；根据闪烁脉冲的噪声模型计算给定能量值的后验概率对数值；通过计算反复地调用第二步骤，获得满足最大后验概率条件的能量值。一种闪烁脉冲信息复原的系统，其包括涨落模型模块、后验概率模块及能量值搜索模块。本发明专利技术闪烁脉冲信息复原的方法及系统，有效提高了系统能量计算精确度，特别适合于稀疏量化电平ADC数字式核仪器的能量计算。

【技术实现步骤摘要】
一种闪烁脉冲信息复原的方法及系统
本专利技术涉及数字信号处理、光电信号处理和辐射探测领域，尤其涉及一种闪烁脉冲信息复原的方法及系统。
技术介绍
在正电子寿命谱仪、正电子角度-动量关联分析仪等核分析领域，双符合高能粒子甄别器等核探测领域，以及正电子发射断层成像（PositronEmissionTomography，以下简称PET）等医学影像领域，探测高能粒子的主要工作原理是把高能射线转化为电信号，再利用快电子学得到粒子事件的各种信息。在一些带有符合逻辑运算和时间标记的仪器应用中，对闪烁脉冲进行整形放慢会影响其时间性能，并增加事件脉冲发生堆积的概率，所以将闪烁脉冲直接进行数字化更加能够满足这种需求。较为完整的对闪烁脉冲进行数字化是采用高带宽高采样率的模拟数字转换器（AnalogtoDigitalConverter，以下简称ADC）直接对闪烁脉冲进行采样和量化。这种数字化方法的成本较高，不能满足多通道系统的实际需求。为了减少这种数字化方案的成本，有两种精简模式，一种是采用降低采样率的方式。由于采样频率较低，无整形的闪烁脉冲信号的高频部分不但不能被采样，反而会由于频域的混叠，影响到奈奎斯特频率附近的分量的精确度。为了缓解这种不足，会在闪烁脉冲输出和数字化之间添加一个整形环节。这个整形环节把闪烁脉冲高于奈奎斯特频率的分量进行衰减而保存了闪烁脉冲的低频分量，其作用是减少闪烁脉冲数字化在频率的混叠，但由于奈奎斯特频率本身的限制，超过该频率的分量始终是不能被采样到的。除了在时间轴上对ADC进行精简以外，还有一种方式是对ADC在电压轴上进行精简。在电压轴上进行精简的ADC可以利用开环设计提高有效带宽，适合于高速信号处理。例如，若干个比较器/ADC单元的数字化方法，或者是开路时间交织的ADC都属于这个范畴。这种ADC的特点是高采样率高带宽，但量化精度有限。由于闪烁脉冲具有比较鲜明的先验知识，通过利用这些知识可以显著提高数字化后信号复原的精度和优化量化电平的设置。因此，针对上述稀疏量化电平ADC，有必要提供一种新的烁脉冲信息复原的方法及系统，以解决现有技术中存在的问题。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种闪烁脉冲信息复原的方法及系统，一方面，显著提高数字化后信号复原的精度；另一方面，优化量化电平的设置可以获取更多的事件信息。为实现上述目的，本专利技术提供如下技术方案：一种闪烁脉冲信息复原的方法，其具体包括步骤：S1：获取低计数下的非堆积符合单事件的闪烁脉冲数据库，再对各路单事件的闪烁脉冲数据库建立闪烁脉冲的噪声模型，对各路非堆积符合单事件的闪烁脉冲数据库求平均脉冲，闪烁脉冲的形状信息由平均脉冲给定；S2：根据闪烁脉冲的噪声模型计算给定能量值的后验概率对数值，具体包括步骤：（2.1）载入一个闪烁脉冲片段S0，该闪烁脉冲片段的起点是向上过阈值Vl的时间点t0，终点是t0+Δt，其中Δt是一个大于闪烁晶体下降沿时间常数两倍的值，根据脉冲噪声模型，计算给定能量值的似然函数；（2.2）将所有时间点取对数，再相加，得到一个单调性与后验概率相同的函数的具体值；S3：通过计算反复地调用步骤S2，获得满足最大后验概率条件的能量值，具体包括步骤：（3.1）将步骤S2作为一个模块反复调用，计算不同试探能量值的后验概率；（3.2）线性搜索满足最大后验概率的能量值；（3.3）计算得到能量值后，根据能量值校准2.1中的时间基点，再重复3.1和3.2。优选的，在上述闪烁脉冲信息复原的方法中，所述步骤S1中闪烁脉冲数据库中的样本数在2000个以上。优选的，在上述闪烁脉冲信息复原的方法中，所述步骤S1中闪烁脉冲的噪声模型包括特征：闪烁脉冲形状、闪烁脉冲高度、闪烁脉冲噪声水平。优选的，在上述闪烁脉冲信息复原的方法中，所述步骤S1中求平均脉冲的具体步骤包括：（1.1）降低射源的辐射剂量得到弱源，通过弱源或调整弱源对探测器的立体角，降低每个探测器捕获的高能光子数，每个探测器接收到的事件是一个泊松流，其平均计数率为其中，mi和qi分别是弱源的剂量和弱源对探测器的立体角，i是弱源的序数，n是弱源的个数；（1.2）按照恒定系数鉴别方法或前沿甄别方法对齐脉冲；（1.3）对对齐后的脉冲取平均。优选的，在上述闪烁脉冲信息复原的方法中，所述步骤S1中对齐脉冲后，同一个能量值的闪烁脉冲在一个时间点上能量值会有一个分散的分布，这个分布被定义为脉冲的噪声模型，噪声模型的参数留作计算给定能量值的后验概率值。优选的，在上述闪烁脉冲信息复原的方法中，所述步骤S1还包括：将一能量区间选择的符合脉冲对认定为单事件脉冲对，将此数字化的电脉冲预先存储后，进行离线分析，通过对齐脉冲对，获得平均脉冲信号，将平均脉冲作为系统响应，将数字脉冲加和，得到单事件数据的能谱。一种闪烁脉冲信息复原的系统，其包括涨落模型模块、后验概率模块及能量值搜索模块，其中，所述涨落模型模块用于对预采集的低计数率数据进行训练，获得闪烁探测器系统的涨落模型，该涨落模型表征平均脉冲、噪声方差和噪声分布偏度，该噪声模型提供给后验概率模块；所述后验概率模块用于后验概率的计算，获得给定能量值的后验概率对数值；所述能量值搜索模块用于搜索满足最大后验概率条件的能量值。优选的，在上述闪烁脉冲信息复原的系统中，所述涨落模型模块包括平均脉冲模块、脉冲噪声模块及概率密度分布模块，其中，所述平均脉冲模块用于计算单事件的平均脉冲；所述脉冲噪声模块用于计算单事件脉冲的噪声方差；所述概率密度分布模块用于估计噪声的概率密度分布函数。优选的，在上述闪烁脉冲信息复原的系统中，所述后验概率模块包括似然函数分布模块及时间点乘积模块，其中，所述似然函数分布模块按照给定能量值和涨落模型模块输入的脉冲噪声模型计算每个点的似然函数值；所述时间点乘积模块负责计算后验概率对数值。优选的，在上述闪烁脉冲信息复原的系统中，所述能量值搜索模块包括试探能量计算模块、线性搜索模块及时间基点校正模块，其中，所述试探能量计算模块通过调用后验概率模块计算出每个试探能量的后验概率对数值；所述线性搜索模块通过线性搜索算法计算满足最大后验概率条件的能量值；所述时间基点校正模块校正样本的时间基点，并重复调用试探能量计算模块和线性搜索模块。从上述技术方案可以看出，本专利技术实施例通过在上述闪烁脉冲信息复原方法与系统中，在系统运行时可有效地计算出各个单事件的能量信息；增加了系统在能量窗内的计数率，提高了系统的能量分辨率，特别适合于离线环境下的闪烁事件处理。从上述技术方案可以看出，本专利技术实施例通过计算给定系统的最大后验估计对应的闪烁事件参数，准确地估计了闪烁脉冲的脉冲高度或者能量信息。与现有技术相比，本专利技术的有益效果是：（1）较好的能量分辨率。（2）给定能量窗内的计数率高。（3）可避免量化电平较少时带来的量化误差。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的有关本专利技术的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术闪烁脉冲信息复原的系统结构图；图2为本专利技术闪烁脉冲信息复原的方法流程图；图3为一种典型的闪烁脉冲示意本文档来自技高网...

【技术保护点】
一种闪烁脉冲信息复原的方法，其特征在于：具体包括步骤：S1：获取低计数下的非堆积符合单事件的闪烁脉冲数据库，再对各路单事件的闪烁脉冲数据库建立闪烁脉冲的噪声模型，对各路非堆积符合单事件的闪烁脉冲数据库求平均脉冲，闪烁脉冲的形状信息由平均脉冲给定；S2：根据闪烁脉冲的噪声模型计算给定能量值的后验概率对数值，具体包括步骤：（2.1）载入一个闪烁脉冲片段S0，该闪烁脉冲片段的起点是向上过阈值Vl的时间点t0，终点是t0+Δt，其中Δt是一个大于闪烁晶体下降沿时间常数两倍的值，根据脉冲噪声模型，计算给定能量值的似然函数；（2.2）将所有时间点取对数，再相加，得到一个单调性与后验概率相同的函数的具体值；S3：通过计算反复地调用步骤S2，获得满足最大后验概率条件的能量值，具体包括步骤：（3.1）将步骤S2作为一个模块反复调用，计算不同试探能量值的后验概率；（3.2）线性搜索满足最大后验概率的能量值；（3.3）计算得到能量值后，根据能量值校准2.1中的时间基点，再重复3.1和3.2。

【技术特征摘要】
1.一种闪烁脉冲信息复原的方法，其特征在于：具体包括步骤：S1：获取低计数下的非堆积符合单事件的闪烁脉冲数据库，再对各路单事件的闪烁脉冲数据库建立闪烁脉冲的噪声模型，对各路非堆积符合单事件的闪烁脉冲数据库求平均脉冲，闪烁脉冲的形状信息由平均脉冲给定；所述步骤S1中所述闪烁脉冲的噪声模型为对齐脉冲后同一个能量值的闪烁脉冲在一个时间点上能量值的分散的分布，所述闪烁脉冲的噪声模型的参数留作计算给定能量值的后验概率值；闪烁脉冲的噪声模型包括特征：闪烁脉冲形状、闪烁脉冲高度、闪烁脉冲噪声水平；所述步骤S1中求平均脉冲的具体步骤包括：(1.1)降低射源的辐射剂量得到弱源，通过弱源或调整弱源对探测器的立体角，降低每个探测器捕获的高能光子数，每个探测器接收到的事件是一个泊松流，其平均计数率为其中，mi和qi分别是弱源的剂量和弱源对探测器的立体角，i是弱源的序数，n是弱源的个数；(1.2)按照恒定系数鉴别方法或前沿甄别方法对齐脉冲；(1.3)对对齐后的脉冲取平均；S2：根据闪烁脉冲的噪声模型计算给定能量值的后验概率对数值，具体包括步骤：(2.1)载入一个闪烁脉冲片段S0，该闪烁脉冲片段的起点是向上过阈值Vl的时间点t0，终点是t0+Δt，其中Δt是一个大于闪烁晶体下降沿时间常数两倍的值，根据脉冲噪声模型，计算给定能量值的似然函数；(2.2)将所有时间点取对数，再相加，得到一个单调性与后验概率相同的函数的具体值；S3：通过计算反复地调用步骤S2，获得满足最大后验概率条件的能量值，具体包括步骤：(3.1)将步骤S2作为一个模块反复调用，计算不同试探能量值的后验概率；(3.2)线性搜索满足最大后验概率的能量值；(3.3)计算得到能量值后，根据能量值校准(2.1)中的时间基点，再重复(3.1)和(3.2)。2.根据权...

【专利技术属性】
技术研发人员：谢庆国，邓贞宙，
申请(专利权)人：苏州瑞派宁科技有限公司，
类型：发明
国别省市：江苏;32