[发明专利]一种PET晶体位置查找表更新方法

权利要求书

1.一种PET晶体位置查找表更新方法，其特征在于，包括以下步骤：

(1)PET探测器预先采集数据获得泛场图像I₀，分割泛场图像I₀，生成位置查找表，并且识别各晶体响应顶点P₀并存储信息，存储信息包括各晶体响应顶点P₀的位置坐标和晶体编号；

(2)当PET系统需要更新晶体位置查找表时，PET探测器重新采集数据，获得新泛场图像I₁；

(3)调用已存储的晶体响应顶点P₀的存储信息，使用迭代算法将各晶体响应顶点P₀收敛至新泛场图像I₁的实际顶点P₁，得到实际顶点P₁的位置坐标；

(4)以实际顶点P₁的位置坐标和晶体编号生成更新的位置查找表，用P₁替换P₀并存储信息。

2.根据权利要求1所述的一种PET晶体位置查找表更新方法，其特征在于：所述PET探测器是由闪烁晶体阵列耦合光电探测器阵列构成。

3.根据权利要求1所述的一种PET晶体位置查找表更新方法，其特征在于：所述泛场图像I₀和新泛场图像I₁，由放射源进行照射或者由晶体本底辐射进行较长时间采集获得。

4.根据权利要求1所述的一种PET晶体位置查找表更新方法，其特征在于：所述晶体响应顶点P₀，根据各晶体响应区域的最高计数点或各晶体响应区域重心位置确定。

5.根据权利要求1所述的一种PET晶体位置查找表更新方法，其特征在于：所述泛场图像I₀和新泛场图像I₁，可进行降低图像噪声处理。

6.根据权利要求5所述的一种PET晶体位置查找表更新方法，其特征在于：所述降低图像噪声处理为平滑滤波处理。

7.根据权利要求1所述的一种PET晶体位置查找表更新方法，其特征在于：所述迭代算法是指，在新泛场图像I₁上以P₀作为初始值，采用寻峰算法进行迭代，寻峰算法为局部极大值移动算法或均值移动算法。

8.根据权利要求7所述的一种PET晶体位置查找表更新方法，其特征在于：所述采用局部极大值移动方法进行迭代，其中采用核K(x)，其中λ为核半径；

P₀初始坐标x，按下面公式计算迭代后坐标m(x),

m(x)＝argmax{I₁(s)·K(x)},s∈S，其中S为泛场位置坐标范围，s为其中一个坐标；

将x更新为m(x)继续迭代，直到收敛，收敛条件|m(x)-x|0.1。

9.根据权利要求1所述的一种PET晶体位置查找表更新方法，其特征在于：所述实际顶点P₁的晶体编号继承至其迭代初始值P₀所对应晶体编号，新泛场图像I₁晶体位置查找表的获得方法是以距离为判据，计算新泛场图像I₁上各像素点与各顶点距离，以距离最近的顶点的晶体编号作为该像素的晶体编号。

10.权利要求1至9中任何一项所述的一种PET晶体位置查找表更新方法，其特征在于：所述方法用于伽马射线探测器的更新晶体位置查找表的用途，伽马射线探测器为成像离散晶体阵列的伽马相机、SPECT、康普顿相机。