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

一种PET三维图像重建方法和装置 技术领域 本发明涉及医疗图像处理技术领域，特别是涉及一种PET三维图像重建 方法和装置。 背景技术 有序子集最大似然法（OrderedSubsetExpectationMaximization，OSEM 是目前PET临床中主要的三维图像重建迭代算法。该算法将投影数据分成n 个子集，每次重建时只使用一个子集对投影数据进行校正，重建图像更新一 次，所有的子集都对投影数据校正一次，称为一次迭代；满足停步规则时， 迭代中止，否则以新的重建图像作为初始图像进入下一次迭代。 然而OSEM算法的缺陷是计算量大，计算速度慢。具体而言，在每次更 新重建图像时，需要计算每条响应线对每个像素的校正系数。在重建图像较 小且系统中响应线数目较少时（如二维重建），可以将所有响应线对所有像素 的校正系数预先计算出，然后通过查表法在重建时调用；但在高分辨率重建 时（如三维重建），重建图像很大而且系统中具有海量响应线，现有的数据存 储介质无法满足要求，只能在重建时实时计算每条响应线对每个像素的校正 系数。 因此，需要提出一种新的PET三维图像重建方法和装置，以提高三维图 像重建的效率。 发明内容 本发明解决的技术问题是提供一种新的PET三维图像重建方法和装置， 以提高三维重建的效率。 本发明的实施例提供了一种PET图像三维重建方法，所述方法包括： 1）初始化三维图像，使初始三维图像位于第一坐标系（X，Y，Z）中三 维矩阵上像素点的像素值为f1（x，y，z）； 2）从探测器获得的多个响应线子集中选择一个响应线子集，每一个所述 响应线子集包含多个响应线组，每一个响应线组包含多条相互平行的响应线； 3）依次从所选择的响应线子集中选择响应线组，并分别对每一个所选择 的响应线组进行预设操作，其中所述预设操作包括： 3.1）根据第一坐标系（X，Y，Z）中的像素点的像素值f1（x，y，z）， 获得所选择响应线组对应的位于第二坐标系（S，T，Z）中的第一像素组的每 个像素点的像素值g1（s，t，z），所述第一像素组的像素点位于所选择的响应 线组对应的多条响应线上； 3.2）根据所选择响应线组对应的位于第二坐标系中的第一像素组的每个 像素点的像素值g1（s，t，z）进行正投影计算，以获得所选择响应线组对应 的各投影点的正投影值p1（s，t，z）； 3.3）将所选择响应线组对应的各投影点的正投影值p1（s，t，z）与测量 投影值p0（s，t，z）进行对比，以获得校正系数R（s，t，z），所述测量投影 值p0（s，t，z）根据所述探测器获得的响应线数据计算得到； 3.4）根据所述校正系数R（s，t，z）进行反投影计算，以获得所选择响 应线组对应的位于第二坐标系中的第一像素组的每个像素点的像素值g2（s， t，z）； 3.5）根据所选择响应线组对应的位于第二坐标系中的第一像素组的每个 像素点的像素值g2（s，t，z），获得所选择响应线组对应的位于第一坐标系中 的第二像素组的每个像素点的像素值f2（x，y，z）； 4）根据所选择响应线子集的多个响应线组对应的所述预设操作的结果更 新所述初始三维图像； 5）满足停步规则时，迭代中止，否则根据更新后的初始三维图像从步骤 2）开始继续执行。 可选地，所述步骤1）中所述第一坐标系（X，Y，Z）的z轴的方向为所 述探测器的轴线方向；所述步骤3.2）中所述第二坐标系（S，T，Z）的z轴 方向为所述探测器的轴线方向，所述第二坐标系的t轴方向为所选择响应线组 对应的响应线方向，s轴方向为所选择响应线组对应的响应线方向的垂直方向。 可选地，所述方法还包括：进行正投影计算之前，在s-t平面内，对所选 择响应线组对应的位于第二坐标系中的第一像素组的每个像素点的像素值g1 （s，t，z）正向应用点扩展函数；进行反投影计算之后，在s-t平面内，对所 选择响应线组对应的位于第二坐标系中的第一像素组的每个像素点的像素值 g2（s，t，z）反向应用点扩展函数。 可选地，所述方法还包括：进行正投影计算之前，沿着s轴的方向，对所 选择响应线组对应的位于第二坐标系中的第一像素组的每个像素点的像素值 g1（s，t，z）正向应用点扩展函数；进行反投影计算之后，沿着s轴的方向， 对所选择响应线组对应的位于第二坐标系中的第一像素组的每个像素点的像 素值g2（s，t，z）反向应用点扩展函数。 可选地，所述步骤3.2）中的正投影计算包括：在第一预设范围内将与投 影点相关的像素点的像素值进行加权求和；所述步骤3.4）中的反投影计算包 括：在第二预设范围内将与像素点相关的投影点的投影值进行加权求和。 可选地，所述第一预设范围根据所述投影点和所述投影点所在的与t-z平 面平行的所述探测器包含的探测环的剖面的端点之间的位置关系来确定；所 述第二预设范围根据所述