一种探测器晶体阵列分辨图处理方法和装置的制造方法
【技术领域】
[0001] 本申请设及PET探测器领域，具体设及一种探测器晶体阵列分辨图的处理方法和 装置。
【背景技术】
[0002] 高分辨率小动物PET探测器通常由分割晶体阵列加多通道光电倍增管或SiPM阵列 组成，运类探测器通常晶体单元小于光探测器的单元大小，首先由测量到的信号计算晶体 分辨图，然后通过分割晶体分辨图来确定发生相互作用的晶体，最终探测器的分辨率取决 于能分辨多小的晶体单元。
[0003] 使用多通道光电倍增管或娃光电倍增管阵列读出晶体阵列，光子主要被发生相互 作用的晶体附近少数光探测器单元测量到，大多数单元的信噪比很低。采用电阻网读出信 号计算晶体分辨图时每个单元的信号权重一样，使得晶体分辨图质量下降。
【发明内容】
[0004] 本申请提供一种探测器晶体阵列分辨图处理方法和装置。
[0005] 根据本申请的第一方面，提供一种探测器晶体阵列分辨图处理方法，包括:产生测 量信号;读出晶体阵列的测量信号;基于所述测量信号，进行权重调整，提高大能量单元的 权重，计算晶体分辨图。
[0006] 所述读出晶体阵列的测量信号包括:单独读出晶体单元的测量信号;所述基于所 述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括:对于晶体位置X 和Y，沉积能量计算公式：
[0010] 其中Xi、yi是第i单元的位置坐标，Ei为第i单元现慢到的能量，k可在1.巧Ij4之间调 整。
[0011] 所述读出晶体阵列的测量信号包括:单独读出晶体单元的测量信号;所述基于所 述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括:对于晶体位置X 和Y，沉积能量计算公式：
[0015]其中xi、yi是第i单元的位置坐标，El为第i单元测量到的能量，j为迭代次数，初始 值Xi和Yi计算公式：
[0019] 其中Xi、yi是第i单元的位置坐标，E功第i单元现慢到的能量。
[0020] 所述读出晶体阵列的测量信号包括:行列相加读出晶体单元的测量信号;所述基 于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括:对于晶体 位置X和Y，沉积能量计算公式：
[0024] 其中xi、yi是第i列或行的位置坐标，E功第i列或行测量到的能量，k可在1.巧Ij4之 间调整。
[0025] 所述读出晶体阵列的测量信号包括:行列相加读出晶体单元的测量信号;所述基 于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括:对于晶体 位置X和Y，沉积能量计算公式：
[0030]其中Xi、yi是第i列或行的位置坐标，El为第i列或行测量到的能量，j为迭代次数， 初始值Xi和Yi计算公式：
[0034] 其中xi、yi是第巧喊行的位置坐标，E功第巧喊行测量到的能量。
[0035] 根据本申请的第二方面，提供一种探测器晶体阵列分辨图处理装置，包括:信号模 块，用于产生测量信号;读出模块，用于读出晶体阵列的测量信号;计算模块，用于基于所述 测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图。
[0036] 所述读出模块还用于单独读出晶体单元的测量信号；所述计算模块还用于:对于 晶体位置X和Y，沉积能量计算公式：
[0037]
[0040] 其中xi、yi是第i单元的位置坐标，E功第i单元测量到的能量，k可在1.巧Ij4之间调 整。
[0041] 所述读出模块还用于单独读出晶体单元的测量信号；所述计算模块还用于:对于 晶体位置X和Y，沉积能量计算公式：
[0045]其中xi、yi是第i单元的位置坐标，El为第i单元测量到的能量，j为迭代次数，初始 值Xi和Yi计算公式：
[0049] 其中xi、yi是第i单元的位置坐标，El为第i单元测量到的能量。
[0050] 所述读出模块还用于行列相加读出晶体单元的测量信号;所述计算模块还用于： 对于晶体位置X和Y，沉积能量计算公式：
[0054] 其中xi、yi是第i列或行的位置坐标，E功第i列或行测量到的能量，k可在1.巧Ij4之 间调整。
[0055] 所述读出模块还用于行列相加读出晶体单元的测量信号;所述计算模块还用于： 对于晶体位置X和Y，沉积能量计算公式：
[0060]其中Xi、yi是第i列或行的位置坐标，El为第i列或行测量到的能量，j为迭代次数， 初始值Xi和Yi计算公式：
[0064] 其中xi、yi是第巧喊行的位置坐标，E功第巧喊行测量到的能量。
[0065] 由于采用了 W上技术方案，使本申请具备的有益效果在于：
[0066] 在本申请的【具体实施方式】中，通过调整单元的权重，提高大能量单元的权重，从而 降低信噪比小的单元的影响，能有效提高探测器晶体分辨图的质量。
【附图说明】
[0067] 图1为根据本申请方法一个实施例的流程图；
[0068] 图2为根据本申请装置一个实施例的结构示意图。
【具体实施方式】
[0069] 下面通过【具体实施方式】结合附图对本发明作进一步详细说明。
[0070] 晶体分辨图是一个显示PET探测器测量到的伽玛事件位置的二维直方图，由于所 用晶体阵列由许多单个晶体组成，在晶体分辨图中每一个晶体单元对应一个点团，点团的 大小和距离表明探测器是否可W清楚分辨每个晶体单元。传统使用多通道光电倍增管或新 型的娃光电倍增管阵列光探测器的PET探测器，使用电阻网读电路将每个光探测器的信号 数降为4个，对于双端读出PET探测器通常使用下面公式计算事件在晶体分辨图中的位置：
[0071] Xi = (Bi+Ci) /El，Yi = (Ci+Di) /El，拉=(B2+C2) /E2，Y2 = (C2+D2) /E2 (1)
[00巧心、81、扣、01和42、82心、化分别为两个光探测器电阻网读出的4个位置相关的能量 信号。通过上面公式可W得到由每个光探测器分别计算的晶体分辨图。探测器的能量则由 下面公式计算：
[0073] Ei = Ai+Bi+Ci+Di ,E2=A2+B2+C2+02 ,E = Ei+E2 (2)
[0074] 通常对于双端读出深度测量探测器由单个光探测器计算的晶体分辨图随着深度 而改变，不但探测器测量到的光子数随着深度远离该探测器而减少，晶体分辨图质量变差， 而且由于光子传输中有可能穿过反射膜而造成多数晶体在晶体分辨图中的位置也随着深 度而改变，使得整个探测器的晶体分辨图变差，可W使用W下两种方法通过结合两端探测 器的能量信号而使探测器的晶体分辨图得到改进：
[0075] χ=〇(ι 巧 2)/2 Y=(Yi+Y2)/2 (3)
[0076] 或X=(Bl+Cl+B2+C2)/EY=(Cl+Dl+C2+D2)/E (4)
[0077] 实施例一：
[0078] 图1示出根据本申请方法的一个实施例的流程图，包括：
[00巧]步骤102:产生测量信号。
[0080] 步骤104:读出晶体阵列的测量信号。
[0081] 在一种实施方式中，可W采用单独读出，另在一种实施方式中可W采用行列相加 读出。
[0082] 步骤106:基于测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨 图。
[0083]在一种实施方式中，对于单独读出晶体单元的测量信号，采用能量加权算法:对于 晶体位置X和Y，沉积能量计算公式：
[0087] 其中xi、yi是第i单元的位置坐标，E功第i单元测量到的能量，k可在1.巧Ij4之间调 整，W获得最佳晶体分辨图。
[0088] 另在一种实施方式中，对于单独读出晶体单元的测量信号，采用位置加权迭代算 法:对于晶体位置X和Y，沉积能量计算公式：
[0092]其中xi、yi是第i单元的位置坐标，El为第i单元测量到的能量，j为迭代次数，k是一 个可调参数，可通过调整k达到最佳晶体分辨图。初始值Xi和Yi计算公式：
[0096] 其中xi、y