一种PET二维图像重建方法及装置发明专利.docxVIP

一种PET二维图像重建方法及装置 技术领域 本发明涉及正电子发射断层扫描(positron emission tomography，PET)技术领 域，尤其涉及一种基于飞行时间技术的PET图像的二维重建方法和装置。 背景技术 有序子集最大似然法(Ordered Subset Expectation Maximization，OSEM)是近 年来发展完善的快速迭代重建算法，它具有空间分辨好，抗噪能力强，速度快于其它迭代法 等优点，已在新型的核医学断层影像设备中广为应用。OSEM算法将投影数据分成n个子集， 每次重建时只使用一个子集对投影数据进行校正，重建图像更新一次，这样所有的子集都 对投影数据校正一次，称为一次迭代。和传统的最大似然法(Maximum Likelihood Expectation-Maximization，MLEM)相比，在近似相同的计算时间和计算量下，重建图像更 新次数是MLEM的n倍，大大加快了图像重建速度，缩短了重建时间。但是在符合事件数目很 大时，OSEM算法计算量还是很大，重建所需时间较长，其中计算图像投影是最耗时的步骤。 在计算正投影和反投影时，由于不同角度下每条响应线与图像相交的像素坐标和加权系数 不尽相同，需要反复调用快速射线驱动算法如siddon算法进行计算，这成为整个重建算法 发展的瓶颈。 发明内容 本发明要解决的问题是提供一种PET二维图像重建方法和装置，对现有的OSEM算 法结构进行优化，使用该方法和装置重建PET图像，能够显著降低运算时间。 为解决上述问题，本发明提供了一种PET二维图像重建方法，包括以下步骤： 1)为需要重建的一层创建一个二维初始图像，所述初始图像位于第一坐标系(X， Y)中的各像素点的像素值为f1(x，y)； 2)将所述需要重建的一层所有的响应线划分为多个子集，每一个所述子集包含多 个响应线组，每一个响应线组包含多条相互平行的响应线； 3)从多个子集中选择一个子集，对所选择的子集中的每个响应线组进行预设操 作，并根据所选择的子集的多个响应线组对应的所述预设操作的结果更新所述初始图像， 其中所述预设操作包括： 3.1)根据第一坐标系(X，Y)中的像素点的像素值f1(x，y)，计算位于第二坐标系 (S，T)中的第一像素组的每个像素点的第一像素值g1(s，t)，所述第一像素组的像素点位于 所选择的响应线组对应的多条响应线上； 3.2)根据得到的所述第一像素组的每个像素点的第一像素值g1(s，t)进行正投影 计算，以获得所述第一像素组的各投影点的理论投影值p1(s，t)； 3.3)将所述第一像素组的各投影点的理论投影值p1(s，t)与实际测量投影值p0(s， t)进行对比，以获得校正系数R(s，t)，并根据所述校正系数R(s，t)进行反投影计算，以获得 所述第一像素组的每个像素点的第二像素值g2(s，t)； 3.4)根据位于第二坐标系(S，T)中的第一像素组的每个像素点的第二像素值g2?(s，t)，计算位于第一坐标系(X，Y)中的第二像素组的每个像素点的像素值f2(x，y)，所述第 二像素组的每个像素点分布在所选择的响应线组对应的多条响应线上； 3.5)对所述第二像素组的每个像素点的像素值f2(x，y)做归一化处理； 4)以更新后的图像作为初始图像，继续执行步骤3)，直至满足停步规则。 进一步地，所述步骤3.1)中，第二坐标系(S，T)的t轴的方向与所选择的响应线组 对应的响应线方向平行。 进一步地，所述步骤3.2)中的正投影计算包括：在第一预设范围内将与投影点相 关的像素点的像素值进行加权求和；所述步骤3.3)中的反投影计算包括：在第二预设范围 内将与像素点相关的投影点的投影值进行加权求和。 进一步地，所述步骤3.2)中，进行正投影计算之前，还包括：将第一像素组中沿t轴 方向分布的m个像素点进行像素合并，其中m为大于1的整数；所述步骤3.3)中，实际测量投 影值的获取方法还包括：在第二坐标系中将沿t轴方向分布的m个投影点的实际测量投影值 进行合并。 进一步地，所述步骤3.1)中，计算所述第一像素组的每个像素点的第一像素值g1?(s，t)的步骤包括： a1)对于每个所述第一像素组的像素点(s1，t1)，在第一坐标系中寻找与之对应的 像素点的坐标(x1，y1)； a2)在第一坐标系中，通过与所述像素点(x1，y1)相邻的像素点的像素值获得所述 像素点(x1，y1)的像素值，所述第一像素组的像素点(s1，t1)的第一像素值g1(s1，t1)与位 于第一坐标系中的所述像素点(x1，y1)的像素值相同。 进一步地，所述步骤a1)包括：根据第一坐标系旋转至第二坐标系的角度，计算 与每个所述第一像素组的像素点(s1，t1)相对应的位于第一坐标系中的像素点的坐标(x1， y1)，计算方法