[发明专利]一种处理PET图像的方法、装置以及计算机存储介质

权利要求书

1.一种处理PET图像的方法，其特征在于，所述方法包括以下步骤：

步骤S1：获取PET图像的原始数据；

步骤S2：对所述PET图像的原始数据进行建模，模型包括加性噪声和泊松噪声；

步骤S3：对所述模型进行第一变换，所述第一变换通过以下方式实现：

其中，G表示所述第一变换的函数，x表示所述PET图像的强度，α表示所述PET图像的像素数据和所述泊松噪声之间关系的尺度项，μ和σ分别表示所述加性噪声的期望和标准差；

步骤S4：对经过所述第一变换的数据进行逆变换，所述逆变换为：

IG(y)＝∫G(x)p(x|λ,α)dx

其中，IG表示所述逆变换，p(x|λ,α)表示关于x的概率密度函数。

2.根据权利要求1所述的处理PET图像的方法，其特征在于，所述PET图像的原始数据包括从PET探测器获取的原始计数数据以及伽马光子对应脉冲信号的电压、时间、位置、能量信息。

3.根据权利要求1所述的处理PET图像的方法，其特征在于，在所述步骤S2中，所述模型通过PET图像的强度x表示：

x＝αp+n，

其中，p表示具有潜在期望值λ的泊松噪声的分布，n表示具有期望值μ和标准差σ的高斯噪声。

4.根据权利要求1所述的处理PET图像的方法，其特征在于，在所述步骤S3中，随机选择参数α的值并形成先验数据库。

5.一种处理PET图像的方法，其特征在于，所述方法包括以下步骤：

步骤T1：获取PET图像的原始数据；

步骤T2：在所述PET图像中任意搜索一窗口图像；

步骤T3：在所述窗口图像区域内采样体素；

步骤T4：采用区域均值法恢复所述体素的实际值，所述区域均值法通过如下公式实现：

其中，RM表示区域均值，Ω为所述步骤T3中所述PET图像的采样区域，w为权重系数，

其中，‖A_u-A_v‖²表示以体素u和v为中心的两个块的两个强度矢量之间的距离，Z(u)表示归一化常数。

6.根据权利要求5所述的处理PET图像的方法，其特征在于，所述PET图像的原始数据包括从PET探测器获取的原始计数数据以及伽马光子对应脉冲信号的电压、时间、位置、能量信息。

7.根据权利要求5所述的处理PET图像的方法，其特征在于，在所述步骤T2中，所述窗口图像的大小标记为M×M×M，M为自然数，将体素u对应的PET图像的强度记为x(u)，体素v对应的PET图像的强度记为x(v)，然后将不同的所述窗口图像进行排序标记，若x(u)>x(v)，则使得L(u)>L(v)，L表示对所述窗口图像进行标记的序号。

8.根据权利要求7所述的处理PET图像的方法，其特征在于，在所述步骤T3中，采样体素时遵循如下规则：

第一、在具有类别标记L的区域中随机采样min{M³,D(L(u))}个体素，D表示具有与体素u相同的标记的体素的数量；

第二、当L(u)>1时，在具有类别标记L(u)-1的第一邻近区域中随机采样min{M³,D(L(u)-1)}个体素；

第三、当L(u)<maxL(u)时，在具有类别标记L(u)+1的第二邻近区域中随机地采样min{M³,D(L(u)+1)}个体素。