CN105631824A - 一种探测器晶体阵列分辨图处理方法和装置 - Google Patents

Abstract

本申请公开了一种探测器晶体阵列分辨图处理方法，包括：产生测量信号；读出晶体阵列的测量信号；基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图。本申请还公开了一种探测器晶体阵列分辨图处理装置。本申请通过调整单元的权重，提高大能量单元的权重，从而降低信噪比小的单元的影响，能有效提高探测器晶体分辨图的质量。

Description

一种探测器晶体阵列分辨图处理方法和装置

技术领域

本申请涉及PET探测器领域，具体涉及一种探测器晶体阵列分辨图的处理方法和装置。

背景技术

高分辨率小动物PET探测器通常由分割晶体阵列加多通道光电倍增管或SiPM阵列组成，这类探测器通常晶体单元小于光探测器的单元大小，首先由测量到的信号计算晶体分辨图，然后通过分割晶体分辨图来确定发生相互作用的晶体，最终探测器的分辨率取决于能分辨多小的晶体单元。

使用多通道光电倍增管或硅光电倍增管阵列读出晶体阵列，光子主要被发生相互作用的晶体附近少数光探测器单元测量到，大多数单元的信噪比很低。采用电阻网读出信号计算晶体分辨图时每个单元的信号权重一样，使得晶体分辨图质量下降。

发明内容

本申请提供一种探测器晶体阵列分辨图处理方法和装置。

根据本申请的第一方面，提供一种探测器晶体阵列分辨图处理方法，包括：产生测量信号；读出晶体阵列的测量信号；基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图。

所述读出晶体阵列的测量信号包括：单独读出晶体单元的测量信号；所述基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括：对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，k可在1.5到4之间调整。

所述读出晶体阵列的测量信号包括：单独读出晶体单元的测量信号；所述基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括：对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\left(x_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\left(y_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$E={\mathrm{\Σ}}_1^n{E}_i\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，j为迭代次数，初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量。

所述读出晶体阵列的测量信号包括：行列相加读出晶体单元的测量信号；所述基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括：对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，k可在1.5到4之间调整。

所述读出晶体阵列的测量信号包括：行列相加读出晶体单元的测量信号；所述基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括：对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\right(x_i\·E_i)\·e^{-k{\left(x_i-X_j\right)}^2})/E_x$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\right(y_i\·E_i)\·e^{-k{(y_i-Y_j)}^2})/E_y$

$E_x={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(x_i-X_j\right)}^2}$

$E_y={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(y_i-Y_j\right)}^2}$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，j为迭代次数，初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量。

根据本申请的第二方面，提供一种探测器晶体阵列分辨图处理装置，包括：信号模块，用于产生测量信号；读出模块，用于读出晶体阵列的测量信号；计算模块，用于基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图。

所述读出模块还用于单独读出晶体单元的测量信号；所述计算模块还用于：对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，k可在1.5到4之间调整。

所述读出模块还用于单独读出晶体单元的测量信号；所述计算模块还用于：对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\left(x_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\left(y_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$E={\mathrm{\Σ}}_1^n{E}_i\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，j为迭代次数，初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量。

所述读出模块还用于行列相加读出晶体单元的测量信号；所述计算模块还用于：对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，k可在1.5到4之间调整。

所述读出模块还用于行列相加读出晶体单元的测量信号；所述计算模块还用于：对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\right(x_i\·E_i)\·e^{-k{\left(x_i-X_j\right)}^2})/E_x$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\right(y_i\·E_i)\·e^{-k{(y_i-Y_j)}^2})/E_y$

$E_x={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(x_i-X_j\right)}^2}$

$E_y={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(y_i-Y_j\right)}^2}$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，j为迭代次数，初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量。

由于采用了以上技术方案，使本申请具备的有益效果在于：

在本申请的具体实施方式中，通过调整单元的权重，提高大能量单元的权重，从而降低信噪比小的单元的影响，能有效提高探测器晶体分辨图的质量。

附图说明

图1为根据本申请方法一个实施例的流程图；

图2为根据本申请装置一个实施例的结构示意图。

具体实施方式

下面通过具体实施方式结合附图对本发明作进一步详细说明。

晶体分辨图是一个显示PET探测器测量到的伽玛事件位置的二维直方图，由于所用晶体阵列由许多单个晶体组成，在晶体分辨图中每一个晶体单元对应一个点团，点团的大小和距离表明探测器是否可以清楚分辨每个晶体单元。传统使用多通道光电倍增管或新型的硅光电倍增管阵列光探测器的PET探测器，使用电阻网读电路将每个光探测器的信号数降为4个，对于双端读出PET探测器通常使用下面公式计算事件在晶体分辨图中的位置：

X₁＝(B₁+C₁)/E₁,Y₁＝(C₁+D₁)/E₁,X₂＝(B₂+C₂)/E₂,Y₂＝(C₂+D₂)/E₂(1)

A₁、B₁、C₁、D₁和A₂、B₂、C₂、D₂分别为两个光探测器电阻网读出的4个位置相关的能量信号。通过上面公式可以得到由每个光探测器分别计算的晶体分辨图。探测器的能量则由下面公式计算：

E₁＝A₁+B₁+C₁+D₁,E₂＝A₂+B₂+C₂+D₂,E＝E₁+E₂(2)

通常对于双端读出深度测量探测器由单个光探测器计算的晶体分辨图随着深度而改变，不但探测器测量到的光子数随着深度远离该探测器而减少，晶体分辨图质量变差，而且由于光子传输中有可能穿过反射膜而造成多数晶体在晶体分辨图中的位置也随着深度而改变，使得整个探测器的晶体分辨图变差，可以使用以下两种方法通过结合两端探测器的能量信号而使探测器的晶体分辨图得到改进：

X＝(X₁+X₂)/2Y＝(Y₁+Y₂)/2(3)

或X＝(B₁+C₁+B₂+C₂)/EY＝(C₁+D₁+C₂+D₂)/E(4)

实施例一：

图1示出根据本申请方法的一个实施例的流程图，包括：

步骤102：产生测量信号。

步骤104：读出晶体阵列的测量信号。

在一种实施方式中，可以采用单独读出，另在一种实施方式中可以采用行列相加读出。

步骤106：基于测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图。

在一种实施方式中，对于单独读出晶体单元的测量信号，采用能量加权算法：对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k---\left(5\right)$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，k可在1.5到4之间调整，以获得最佳晶体分辨图。

另在一种实施方式中，对于单独读出晶体单元的测量信号，采用位置加权迭代算法：对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\left(x_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\left(y_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$E={\mathrm{\Σ}}_1^n{E}_i\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}---\left(6\right)$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，j为迭代次数，k是一个可调参数，可通过调整k达到最佳晶体分辨图。初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i---\left(7\right)$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量。

又在一种实施方式中，对于行列相加读出晶体单元的测量信号，采用能量加权算法：对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k---\left(8\right)$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，k可在1.5到4之间调整。

再在一种实施方式中，对于行列相加读出晶体单元的测量信号，采用位置加权迭代算法：对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\right(x_i\·E_i)\·e^{-k{\left(x_i-X_j\right)}^2})/E_x$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\right(y_i\·E_i)\·e^{-k{(y_i-Y_j)}^2})/E_y$

$E_x={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(x_i-X_j\right)}^2}$

$E_y={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(y_i-Y_j\right)}^2}---\left(9\right)$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，j为迭代次数，k是一个可调参数，可通过调整k达到最佳晶体分辨图。初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i---\left(10\right)$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量。

实施例二：

图2示出了根据本申请装置一个实施例的结构示意图，包括：信号模块，用于产生测量信号；读出模块，用于读出晶体阵列的测量信号；计算模块，用于基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图。

在一种实施方式中可以读出模块可以采用单独读出，另在一种实施方式中读出模块可以采用行列相加读出。

对于单独读出晶体单元的测量信号，计算模块采用能量加权算法：对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k---\left(5\right)$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，k可在1.5到4之间调整，以获得最佳晶体分辨图。

对于单独读出晶体单元的测量信号，计算模块采用位置加权迭代算法：对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\left(x_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\left(y_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$E={\mathrm{\Σ}}_1^n{E}_i\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}---\left(6\right)$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，j为迭代次数，k是一个可调参数，可通过调整k达到最佳晶体分辨图。初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i---\left(7\right)$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量。

对于行列相加读出晶体单元的测量信号，计算模块采用能量加权算法：对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k---\left(8\right)$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，k可在1.5到4之间调整。

对于行列相加读出晶体单元的测量信号，计算模块采用位置加权迭代算法：对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\right(x_i\·E_i)\·e^{-k{\left(x_i-X_j\right)}^2})/E_x$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\right(y_i\·E_i)\·e^{-k{(y_i-Y_j)}^2})/E_y$

$E_x={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(x_i-X_j\right)}^2}$

$E_y={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(y_i-Y_j\right)}^2}---\left(9\right)$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，j为迭代次数，k是一个可调参数，可通过调整k达到最佳晶体分辨图。初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i---\left(10\right)$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量。

由上述实施方式计算得到每端探测器的晶体分辨图，然后可以由公式3得到整个探测器的晶体分辨图，对于行列相加读出，也可以由一端探测器计算X，另一端探测器计算Y。

对本申请使用于单端读出探测器进行了初步的检验，将一个6×6、晶体大小1.3×1.3×20毫米、晶体表面抛光和使用ESR的LSO阵列探测器通过硅油耦合于一个HamamatsuH7546M64多阳极PMT的中心，PMT单个阳极大小为2.2×2.2毫米，将PMT中心16个单元单独读出，分别使用公式5、6和7来计算晶体分辨图，本申请能量加权方法和位置加权迭代方法得到的晶体分辨图的质量和传统的重心方法相比有明显的改进。

以上内容是结合具体的实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换。

Claims (10)

1.一种探测器晶体阵列分辨图处理方法，其特征在于，包括：

产生测量信号；

读出晶体阵列的测量信号；

基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图。

2.如权利要求1所述的方法，其特征在于,所述读出晶体阵列的测量信号包括：单独读出晶体单元的测量信号；

所述基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括：

对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，k可在1.5到4之间调整。

3.如权利要求1所述的方法，其特征在于,所述读出晶体阵列的测量信号包括：单独读出晶体单元的测量信号；

所述基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括：

对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\left(x_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\left(y_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$E={\mathrm{\Σ}}_1^n{E}_i\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，j为迭代次数，初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量。

4.如权利要求1所述的方法，其特征在于,所述读出晶体阵列的测量信号包括：行列相加读出晶体单元的测量信号；

所述基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括：

对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，k可在1.5到4之间调整。

5.如权利要求1所述的方法，其特征在于,所述读出晶体阵列的测量信号包括：行列相加读出晶体单元的测量信号；

所述基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图包括：

对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\right(x_i\·E_i)\·e^{-k{\left(x_i-X_j\right)}^2})/E_x$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\right(y_i\·E_i)\·e^{-k{(y_i-Y_j)}^2})/E_y$

$E_x={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(x_i-X_j\right)}^2}$

$E_y={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(y_i-Y_j\right)}^2}$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，j为迭代次数，初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量。

6.一种探测器晶体阵列分辨图处理装置，其特征在于，包括：

信号模块，用于产生测量信号；

读出模块，用于读出晶体阵列的测量信号；

计算模块，用于基于所述测量信号，进行权重调整，提高大能量单元的权重，计算晶体分辨图。

7.如权利要求6所述的装置，其特征在于,所述读出模块还用于单独读出晶体单元的测量信号；

所述计算模块还用于：

对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，k可在1.5到4之间调整。

8.如权利要求6所述的装置，其特征在于,所述读出模块还用于单独读出晶体单元的测量信号；

所述计算模块还用于：

对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\left(x_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\left(y_i\·E_i\right)\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}\right)/E$

$E={\mathrm{\Σ}}_1^n{E}_i\·e^{-k\[{\left(x_i-X_j\right)}^2+{\left(y_i-Y_j\right)}^2\]}$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量，j为迭代次数，初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i$

其中x_i、y_i是第i单元的位置坐标，E_i为第i单元测量到的能量。

9.如权利要求6所述的装置，其特征在于,所述读出模块还用于行列相加读出晶体单元的测量信号；

所述计算模块还用于：

对于晶体位置X和Y，沉积能量计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i^k)/E^k$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i^k)/E^k$

$E^k={\mathrm{\Σ}}_1^n{E}_i^k$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，k可在1.5到4之间调整。

10.如权利要求6所述的装置，其特征在于,所述读出模块还用于行列相加读出晶体单元的测量信号；

所述计算模块还用于：

对于晶体位置X和Y，沉积能量计算公式：

$X_{j+1}=\left({\mathrm{\Σ}}_1^n\right(x_i\·E_i)\·e^{-k{\left(x_i-X_j\right)}^2})/E_x$

$Y_{j+1}=\left({\mathrm{\Σ}}_1^n\right(y_i\·E_i)\·e^{-k{(y_i-Y_j)}^2})/E_y$

$E_x={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(x_i-X_j\right)}^2}$

$E_y={\mathrm{\Σ}}_1^n{E}_i\·e^{-k{\left(y_i-Y_j\right)}^2}$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量，j为迭代次数，初始值X₁和Y₁计算公式：

$X=({\mathrm{\Σ}}_1^n{x}_i\·E_i)/E$

$Y=({\mathrm{\Σ}}_1^n{y}_i\·E_i)/E$

$E={\mathrm{\Σ}}_1^n{E}_i$

其中x_i、y_i是第i列或行的位置坐标，E_i为第i列或行测量到的能量。