专利名称：基于ps-pmt的探测器信号读出方法和系统的制作方法
技术领域：
本发明涉及核探测技术和核医学成像技术领域，且特别涉及基于多阳极位置灵敏型光电倍增管的探测器信号读出方法和系统。
背景技术：
目前基于闪烁晶体和光电倍增管(photomultiplier tube, PMT)的核探测技术已被广泛应用于以正电子发射计算机断层成像(Positron EmissionComputed Tomography,PET)和单光子发射计算机断层成像(Single PhotonEmission Computed Tomography,SPECT)为代表的核医学成像领域。 作为成像过程中的起始环节，探测器的空间分辨能力直接决定了成像系统的空间分辨率指标，早期的探测器采用独立晶体与单通道的PMT —对一耦合，这种方式原理简单，但由于PMT尺寸较大导致分辨率不高且几何安排非常麻烦，同时所需PMT的数量巨大，成本高昂。随后出现了使用多个PMT组成阵列与晶体阵列耦合，利用闪烁光的分配实现晶体阵列的位置分辨，空间分辨率有所提高。随着PMT技术的发展，一些专门用于位置探测的多阳极位置灵敏光电倍增管(Position Sensitive Photomultipler tube, PS-PMT)被研制出来(如Hamamatsu的H8500、H9500等)。这类PS-PMT输出阳极通常由阵列组成，与单通道的PMT相比，阵列单元的尺寸更小，空间分辨能力更高；同时PS-PMT的外围死区较小，采用多个PS-PMT组成更大面积阵列时几何排布要简单得多。无论是采用单个PS-PMT或者PS-PMT阵列，探测器的输出通常是一个MXN的阵列，一个最简单并且分辨能力最好的方法就是将这MXN通道的信号单独读出并分别进行放大成形等处理，但这样所需的数据获取系统无疑是非常庞大和昂贵的。通过有效的信号读出方法对MXN通道信号进行缩减，可以在降低数据获取系统负担的同时，获得较高的空间分辨能力。现有技术中的方案通常是将MXN阵列信号转换为正交方向上的M通道信号和N通道信号。在每个方向上，分别对相应的M (或N)通道信号进行前置放大处理，再通过加法电路将各通道信号按比例权重进行分配，传输到成形放大器中进行成形处理，再传输到一板卡上进行模数转换及重心坐标的计算，其中，该加法电路通过选取不同的电阻值实现不同的权重。对于现有技术中的方案，参与坐标计算的通道信号采用加权相加，相应的信号噪声也是M (或N)通道叠加的结果。因此，输出阵列越大，引入的噪声越多，探测器的分辨能力越差，对大规模探测器而言该方案的可拓展性较差。此外，各通道信号的分配是通过选取不同的电阻值实现不同的权重分配，而电阻器件有一定的精度误差，从而导致探测器平面图像中的线性相对较差。因此，业内亟需一种能够在获取Y事例位置信息的过程中，可有效去除无效通道带来的噪声，提高探测器(尤其是大阵列输出探测器)的定位精度的方案。
发明内容
鉴于现有技术中的上述问题，本发明提供了一种基于多阳极位置灵敏型光电倍增管的探测器信号读出方法和系统。本发明的基于多阳极位置灵敏型光电倍增管的探测器信号读出方法，包括步骤SI :通过对称电荷分配电路将该探测器输出的MXN阵列信号转换为正交方向上的M通道信号和N通道信号； 步骤S2 :分别对该M通道信号和N通道信号进行放大成形和数字积分处理，生成M通道和N通道的积分数值；步骤S3 :根据该M通道和N通道的积分数值，采用局域重心法，计算各通道方向上 的重心位置，从而确定Y光子在探测器上的作用位置。本发明的基于多阳极位置灵敏光电倍增管探测器的信号读出系统，包括对称电荷分配电路，用于将该探测器输出的MXN阵列信号转换为正交方向上的M通道信号和N通道信号；处理模块，用于分别对该M通道信号和N通道信号进行放大成形和数字积分处理，生成M通道和N通道的积分数值；计算模块，用于根据该M通道和N通道的积分数值，采用局域重心法，计算各通道方向上的重心位置，从而确定Y光子在探测器上的作用位置。本发明提出的信号读出方法可以在获取Y事例位置信息的过程中有效去除无效通道带来的噪声，提高了探测器(尤其是大阵列输出探测器)的定位精度。
图Ia为本发明的基于PS-PMT的探测器信号读出方法的流程图，图Ib所示为基于PS-PMT的探测器信号读出方法的具体示意图；图2为本发明技术方案中的S⑶电路的示意图；图3为本发明技术方案中的数字积分的示意图；图4为本发明技术方案中的通道筛选的示意图；图5为本发明实施例的探测器散点图。
具体实施例方式体现本发明特征与优点的典型实施例将在以下的说明中详细叙述。应理解的是本发明能够在不同的实施例上具有各种的变化，其皆不脱离本发明的范围，且其中的说明及所附附图在本质上是当作说明之用，而非用以限制本发明。本发明提供了一种基于PS-PMT的探测器信号读出方法，该方法可以在获取Y事例位置信息的过程中去除无效通道带来的噪声，并提高探测器(尤其是大阵列输出探测器)的定位精度。参见图Ia所示为本发明的基于PS-PMT的探测器信号读出方法的流程图，图Ib所示为该方法的具体示意图。结合图Ia和图Ib进行说明，该方法包括步骤SI :通过对称电荷分配电路将该探测器输出的MXN阵列信号转换为正交方向上的M通道信号和N通道信号。
在该步骤中，设探测器的输出阳极为MXN阵列的阳极信号，(其中，M和N为自然数，如M为50，N为40，M和N的值与探测器的性能相关。)通过对称SCD电路(SymmetricCharged Division,电荷分配电路)(如图2所示为S⑶电路的示意图)，将其转换为正交方向上的M通道信号(A:、A2、…、Am)和N通道信号(BpB2'…、Bn)。步骤S2 :分别对该M通道信号和N通道信号进行放大成形和数字积分处理，生成处理后的M通道和N通道的积分数值。本步骤中，分别对ApA2、…、AdPBpBy…、Bn信号作放大成形和数字积分处理(如图3所示),得到各通道的积分值ypy2、…、yM和XpX2、…、xN。参见图3,图中的a)所示为探测器输出的模拟信号经过放大成形后的波形，数字积分如图中的b)所示，就是按一定的时钟周期将该波形划分成等间距的时间间隔，并对每个时间间隔内波形的幅度值求和的过程。以探测器输出的通道信号为电压信号V为例，则对该电压信号V进行积分即将电压信号V的波形划分成等间距的时间间隔，并对每个时间间隔内的电压信号V的波形的幅度值进行求和的过程。
步骤S3 :将经上述处理后的该M通道和N通道的积分数值传输到计算机端，采用局域重心法，计算各通道方向上的重心位置，从而确定Y光子在探测器上的作用位置。本步骤中，将积分值y:、y2、…、yM和…、χΝ传输到计算机，采用TCOG算法计算重心坐标，从而确定Y事例在探测上的作用位置。以χι、χ2、…、χΝ为例，采用TCOG算法的计算流程如下I)对N通道积分值求和Σ 乂厂乂工+…+乂^取和值的一定比例为阈值乂㈦计算公式为X0= (Σ Xi) /b ;其中b的值可以根据探测需求进行设定。2)通道筛选(如图4):将各通道积分值(Xi)与阈值(Xtl)进行比较，若Xi > Xtl，则令 Xi 为 Xi-Xtl ;若 Xi ( X0,则令 Xi 为 O ;如图4所示，图4中的a)为某一方向上的24通道信号积分值的分布图，其中边缘通道的信号幅度很小，主要为噪声，虚线代表计算得到的阈值；图4中的b)为阈值比较后的24通道信号积分值的分布图，从图上可以看到通过阈值比较，边缘通道信号噪声被剔除，只有局部区域的通道信号被保留下来参与位置计算。3)重心坐标计算。对于X方向上的重心坐标X的计算公式为Χ= (Σ iXXi)/ Σ Xi，其中i为通道编号，对于X方向，I取值为I到N，Xi为对应该通道的积分值。对y方向的通道信号进行处理和重心计算的方法与X方向相同，此处不再赘述。基于上述的技术方案，本发明的一种基于PS-PMT探测器的信号读出方法的实施例如下探测器由NaI闪烁晶体阵列(75X75)与3X3的H8500 (Hamamatsu)光电倍增管阵列组成。晶体单元为I. 8mmX I. 8mmX6mm,晶体间隔O. 2mm,整个探测器输出阳极阵列的大小为24X 24。在Na-22源照射下，按照本方案计算每次探测到的事例的重心坐标，并绘制得到如图5所示的散点图，如图5所示，散点图上点与点之间的分开情况反映了探测器的空间分辨能力，可从该图中看出本发明能够具有较好的空间分辨能力。本发明的技术方案的有益效果如下(I)本发明技术方案中通道筛选的对象为各通道的数字积分值，所进行的计算也是线性的数值计算，减少了模拟信号处理方式下电子元器件带来的噪声，提高了探测器的性能。(2)与现有技术相比，本发明技术方案中的通道筛选和重心坐标计算均在计算机端完成，由此避开了电子学处理的局限性，可拓展性强。(3)本发明技术方案中的阈值在计算机端设定，使得探测器的调试更加方便和灵活。此外，本发明技术方案中参与TCOG计算的yi、y2、…、yM和XpX2、…、1]^是对各自通道信号波形的数字积分值。由于各通道信号的波形基本一致，波形的峰值与数字积分值成正比关系，因此作为替代方案，也可以直接选取各自通道信号波形的峰值进行局域重心法计算。本发明还提供了一种基于多阳极位置灵敏光电倍增管探测器的信号读出系统，包括对称电荷分配电路，用于将该探测器输出的MXN阵列信号转换为正交方向上的M通道信号和N通道信号；处理模块，用于分别对该M通道信号和N通道信号进行放大成形和数字积分处理，生成M通道和N通道的积分数值；计算模块，用于根据该M通道和N通道的积分数值，采用局域重心法，计算各通道方向上的重心位置，从而确定Y光子在探测器上的作用位置。其中，该计算模块可以为一计算机，以避免电子学处理的局限性，提高可拓展性。本领域技术人员应当意识到在不脱离本发明所附的权利要求所揭示的本发明的范围和精神的情况下所作的更动与润饰，均属本发明的权利要求的保护范围之内。
权利要求
1.一种基于多阳极位置灵敏型光电倍增管的探测器信号读出方法，其特征在于，包括 步骤SI :通过对称电荷分配电路将该探测器输出的MXN阵列信号转换为正交方向上的M通道信号和N通道信号； 步骤S2 :分别对该M通道信号和N通道信号进行放大成形和数字积分处理，生成M通道和N通道的积分数值； 步骤S3 :根据该M通道和N通道的积分数值，采用局域重心法，计算各通道方向上的重心位置，从而确定Y光子在探测器上的作用位置。
2.根据权利要求I所述的方法，其特征在于，将该M通道和N通道的积分数值传输到计算机端，由计算机端计算各通道方向上的重心的坐标。
3.根据权利要求I所述的方法，其特征在于，该局域重心法的计算方法为对于同一方向上的M通道信号或N通道信号，取各通道积分值的累加和值的一预定比例作为阈值，将各通道积分值与该阈值进行比较，若一通道积分值小于该阈值，则直接将该通道积分值置为零；若一通道积分值大于该阈值，则取该通道积分值与该阈值的差作为该通道的积分值；根据各通道的积分值，计算各通道方向上的重心的坐标。
4.根据权利要求3所述的方法，其特征在于，计算其中的一通道方向上的重心坐标X的方法为Χ= (Σ iXXi)/Σχi，其中i为通道编号，Xi为对应该通道的积分值。
5.一种基于多阳极位置灵敏光电倍增管探测器的信号读出系统，其特征在于，包括 对称电荷分配电路，用于将该探测器输出的MXN阵列信号转换为正交方向上的M通道信号和N通道信号； 处理模块，用于分别对该M通道信号和N通道信号进行放大成形和数字积分处理，生成M通道和N通道的积分数值； 计算模块，用于根据该M通道和N通道的积分数值，采用局域重心法，计算各通道方向上的重心位置，从而确定Y光子在探测器上的作用位置。
6.根据权利要求I所述的系统，其特征在于，所述计算模块为一计算机，该处理模块将该M通道和N通道的积分数值传输到该计算机，由该计算机计算各通道方向上的重心的坐标。
7.根据权利要求5所述的系统，其特征在于，该局域重心法的计算方法为对于同一方向上的M通道信号或N通道信号，取各通道积分值的累加和值的一预定比例作为阈值，将各通道积分值与该阈值进行比较，若一通道积分值小于该阈值，则直接将该通道积分值置为零；若一通道积分值大于该阈值，则取该通道积分值与该阈值的差作为该通道的积分值；根据各通道的积分值，计算各通道方向上的重心的坐标。
8.根据权利要求7所述的系统，其特征在于，计算其中的一通道方向上的重心坐标X的方法为Χ= (Σ iXXi)/Σxi，其中i为通道编号，Xi为对应该通道的积分值。
全文摘要
本发明公开了一种基于多阳极位置灵敏型光电倍增管的探测器信号读出方法和系统，包括通过对称电荷分配电路将该探测器输出的M×N阵列信号转换为正交方向上的M通道信号和N通道信号；分别对该M通道信号和N通道信号进行放大成形和数字积分处理，生成处理后的M通道和N通道的积分数值；将经上述处理后的该M通道和N通道的积分数值传输到计算机端，采用局域重心法，计算各通道方向上的重心位置，从而确定γ光子在探测器上的作用位置。本发明的方案可以在获取γ事例位置信息的过程中去除无效通道带来的噪声，提高了探测器(尤其是大阵列输出探测器)的定位精度。
文档编号G01T1/00GK102944890SQ201210440118
公开日2013年2月27日 申请日期2012年11月6日 优先权日2012年11月6日
发明者魏龙, 黄先超, 李道武, 章志明, 马创新, 王宝义, 李琳, 帅磊, 柴培, 唐浩辉, 姜小盼, 李婷, 王晓明, 朱美玲, 庄凯, 王英杰, 刘彦韬, 张译文 申请人:中国科学院高能物理研究所