一种PET信号的全带压缩感知方法技术

本发明专利技术提供PET信号的全带压缩感知方法,其包括步骤:使用探测器采集能量不同的射源输出的完整PET信号数据库；对步骤1得到的完整PET信号数据库进行全带量化以获取PET信号的标准能量；对步骤1得到的完整PET信号数据库进行全带压缩量化，在PET信号先验信息的基础上估计出PET信号的能量；以步骤2的标准能量为准，确定标准能量与步骤3的估计能量之间的对应关系；用步骤4得到的对应关系校正步骤3的估计能量。该方法使用PET信号的标准能量作为PET信号全带压缩量化后的估计能量的校正标准，避免了PET信号模型和特征描述不精确造成的误差，有效提升了PET信号能量估计的准确度。

【技术实现步骤摘要】
一种PET信号的全带压缩感知方法
本专利技术涉及电离辐射探测以及医学成像领域，更具体地涉及一种PET信号的全带压缩感知方法。
技术介绍
电离辐射探测领域探测器得到的基本信息之一是PET信号的能量，PET信号的能量主要用于：区分电离辐射探测领域使用的射线种类；判断医学成像领域使用的射线与物质是否发生了散射；确定光电位置敏感器件中射线的沉积位置。电离辐射探测领域内，相同探测器输出的PET信号幅值通常随射线在探测器中沉积能量的增多而增长，且上升时间和下降时间始终保持一致。电离辐射领域使用的射线能量通常用PET信号在时间轴上幅值积分的结果表示。传统的获取PET信号能量的方法有两种：一是使用积分电路将PET信号携带的电荷收集起来，然后用高速模数转换器(Analog-to-DigitalConverter，简称ADC)对收集到的最大电荷量进行采样，该采样结果就表示PET信号的能量；二是把PET信号整形成相对速度较慢的PET信号，然后用低速ADC对整形结果进行采样，再对采样结果进行数值积分，该数值积分的结果就表示PET信号的能量。由于传统方法在数值积分的过程中对计数率会有限制作用，且整形电路也易受外界因素(温度、压强等)的干扰而影响性能，同时，整形电路的参数还需要随着不同的应用要求进行调整，导致整个系统的校正和维护十分困难。虽然用高速ADC可以解决上述问题，但它在解决问题的同时，会带来成本和功耗增加的后果，另外，高速ADC对后续数据的处理速度和传输带宽具有更高的要求，极大增强了后续电路的设计难度。为了解决PET信号能量获取过程中高计数率的限制问题，有学者提出基于先验信息的PET信号全带压缩感知的能量获取方法。全带是指在对PET信号进行处理时，处理对象涵盖了PET信号的全部数据。执行该方法需要使用PET信号的先验信息，如PET信号的物理模型、特征等，由于探测器输出的PET信号具有稀疏性(可压缩性)，因此可以对PET信号的所有数据以远低于奈奎斯特采样率的频率进行压缩以获取少量的压缩数据，然后利用统计学的方法实现PET信号的拟合和最大似然求解，获取PET信号的能量。在基于先验信息的PET信号全带压缩感知的能量获取系统中，对信号模型和特征描述的精确度会影响PET信号信息的准确度。在实际研究工作中，PET信号的模型和特征描述不仅与探测器有关，还与其产生电路中的参数有关。因此，要想实现精确的信号模型和特征描述非常困难。当PET信号的理想模型和特征与实际研究过程中的模型和特征不同时，最终估计到的能量就会产生偏差。因此，有必要提出一种PET信号的全带压缩感知方法，以克服上述由于PET信号的模型和特征描述不精确导致估计到的能量发生偏差的问题。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种PET信号的全带压缩感知方法，以提高PET信号能量估计的准确度。为实现上述目的，本专利技术提供的一种PET信号的全带压缩感知方法，包括以下步骤：步骤S1：使用探测器采集能量不同的射源输出的完整PET信号数据库；步骤S2：对S1得到的完整PET信号数据库进行全带量化以获取PET信号的标准能量；步骤S3：对S1得到的完整PET信号数据库进行全带压缩量化，在PET信号先验信息的基础上估计出PET信号的能量；步骤S4：以S2的标准能量为准，确定标准能量与S3的估计能量之间的对应关系；步骤S5：用S4得到的对应关系校正S3的估计能量。进一步的，在步骤S1中，所述数据库可用仿真模拟或直接测试或仿真模拟与直接测试相结合的方式采集。进一步的，在步骤S1中，使用仿真模拟的方式采集数据库的具体过程为：(1)使用同一射源重复照射探测器，并用ADC将探测器输出的数据采集起来；(2)对数据集求平均后进行不同程度的放大，利用放大结果模拟出具有一定能量的数据库。进一步的，所述对数据集求平均的方法为：对探测器每一次输出的数据求和并视为该次测试的能量，测试过程重复多次后作出PET信号的能量谱，选取能量谱峰值前后的第一个谷点作为上下界，将能量位于该上下界区间的数据筛选出来，假设一共能筛出n组数据，定义Pi，j表示第i组数据的第j个采样点，i、j均为整数，且0＜i≤n，0＜j≤m，m为每组数据的采样点总数，因此，对数据集求平均后的结果Pz可以表示为：进一步的，所述数据集求平均后进行放大的方式为模型放大。进一步的，所述模型放大的过程为：确定探测器的能量响应模型y＝f(x)，将平均数据Pz的能量视为Ez，由于实际应用中用户可以将待测射源的能量上限设为Ea，故PET信号需要在放大系数的条件下进行不同程度的放大，已知放大后的PET信号总数为K，则对应放大结果Pd可由计算得到。进一步的，在步骤S1中，使用直接测试的方式采集数据库的具体过程为：首先使用不同能量的射源照射探测器，然后直接将具有一定能量的数据用ADC采集起来。进一步的，在步骤S1中，使用仿真模拟和直接测试相结合的方式采集数据库，具体过程为：(1)使用不同能量的射源重复照射探测器，并用ADC将探测器输出的数据按射源的能量和重复次数分别采集起来；(2)对同一射源的数据集求平均；(3)对平均结果利用插值的方法模拟出具有一定能量的数据库。进一步的，所述对数据集求平均的方法为：对探测器每一次输出的数据求和并视为该次测试的能量，测试过程重复多次后作出PET信号的能量谱，选取能量谱峰值前后的第一个谷点作为上下界，将能量位于该上下界区间的数据筛选出来，假设一共能筛出n组数据，定义Pi，j表示第i组数据的第j个采样点，i、j均为整数，且0＜i≤n，0＜j≤m，m为每组数据的采样点总数，因此，对数据集求平均后的结果Pz可以表示为：进一步的，在步骤S2中，所述PET信号的标准能量可通过对完整数据集直接数字化或对含有PET信号能量信息的数据进行全带压缩获得。进一步的，在步骤S2中，通过对完整数据集直接数字化以获取PET信号标准能量的具体过程为：(1)使用ADC对PET信号进行全带量化；(2)取量化数据的数值积分结果或峰值表示PET信号的标准能量。进一步的，在步骤S2中，通过对完整数据集直接数字化以获取PET信号标准能量的具体过程也可为：(1)令PET信号通过最大值保持电路；(2)对最大值保持电路输出的数据进行全带量化以获取最大值来表示PET信号的标准能量。进一步的，在步骤S2中，通过对含有PET信号能量信息的数据进行全带压缩以获取PET信号标准能量的方法也称为线性放电法。进一步的，在步骤S3中，所述估计PET信号能量的方法为SQL法或ADC压缩拟合法。进一步的，在步骤S3中，使用SQL法估计PET信号能量的具体过程为：(1)设置多个大小在PET信号幅值范围内的过值电压，并将这些过值电压依次与原始PET信号一起输入到比较器内部，令比较器在原始PET信号与大小等于过值电压的恒线相交时输出一个跳变信号并记录跳变发生时对应的时间；<本文档来自技高网...

【技术保护点】
1.一种PET信号的全带压缩感知方法，其特征在于，所述方法包括以下步骤：/n步骤S1：使用探测器采集能量不同的射源输出的完整PET信号数据库；/n步骤S2：对S1得到的完整PET信号数据库进行全带量化以获取PET信号的标准能量；/n步骤S3：对S1得到的完整PET信号数据库进行全带压缩量化，在PET信号先验信息的基础上估计出PET信号的能量；/n步骤S4：以S2的标准能量为准，确定标准能量与S3的估计能量之间的对应关系；/n步骤S5：用S4得到的对应关系校正S3的估计能量。/n

【技术特征摘要】
1.一种PET信号的全带压缩感知方法，其特征在于，所述方法包括以下步骤：
步骤S1：使用探测器采集能量不同的射源输出的完整PET信号数据库；
步骤S2：对S1得到的完整PET信号数据库进行全带量化以获取PET信号的标准能量；
步骤S3：对S1得到的完整PET信号数据库进行全带压缩量化，在PET信号先验信息的基础上估计出PET信号的能量；
步骤S4：以S2的标准能量为准，确定标准能量与S3的估计能量之间的对应关系；
步骤S5：用S4得到的对应关系校正S3的估计能量。


2.根据权利要求l所述的PET信号的全带压缩感知方法，其特征在于，在步骤S1中，所述数据库可用仿真模拟或直接测试或仿真模拟与直接测试相结合的方式采集。


3.根据权利要求2所述的PET信号的全带压缩感知方法，其特征在于，在步骤S1中，使用仿真模拟的方式采集数据库的具体过程为：
(1)使用同一射源重复照射探测器，并用ADC将探测器输出的数据采集起来；
(2)对数据集求平均后进行放大，利用放大结果模拟出具有一定能量的数据库。


4.根据权利要求3所述的PET信号的全带压缩感知方法，其特征在于，所述对数据集求平均的方法为：对探测器每一次输出的数据求和并视为该次测试的能量，测试过程重复多次后作出PET信号的能量谱，选取能量谱峰值前后的第一个谷点作为上下界，将能量位于该上下界区间的数据筛选出来，假设一共能筛出n组数据，定义Pi，j表示第i组数据的第j个采样点，i、j均为整数，且0＜i≤n，0＜j≤m，m为每组数据的采样点总数，...

【专利技术属性】
技术研发人员：邓贞宙，曹香珠，
申请(专利权)人：南昌大学，
类型：发明
国别省市：江西;36