一种pet图像的重建方法及装置的制造方法
【技术领域】
[0001] 本发明设及图像处理领域，尤其是设及一种PET图像的重建方法及装置。
【背景技术】
[0002] 正电子发射型计算机断层显像（英文：Po S i tron Emi S S i on Compu t ed Tomogra地y，PET)是一种常用的医学影像技术。在PET技术中，将生物生命代谢中必须的物 质标记上放射性核素，放射性核素在衰变过程中释放出正电子，正电子运动较小的距离后， 遇到一个电子后发生煙灭事件，从而产生方向大致相反的光子对。因此，通过晶体上获得的 投影数据探测光子对，能够重建出放射性核素分布的Ξ维图像（即PET图像），从而得到放射 性核素在生物生命代谢中聚集的位置，W反映生命代谢活动的情况。
[0003] 然而，在重建PET图像时，由于正电子的煙灭射程、光子对的非准直性、光子对在晶 体中的透射及散射等原因，会导致点扩散问题，影响重建图像的分辨率和对比度。为了解决 点扩散问题，目前重建PET图像时会在投影数据或者重建图像上对点扩散进行校正，下面分 别说明。
[0004] 在投影数据上进行校正的重建方式是:根据关系式y'=Mx'重建出图像X'，其中y' 表示投影数据，Μ为扩散矩阵（即包含扩散信息的系统矩阵），由于投影数据y'能够测得，因 此需要获得扩散矩阵M。扩散矩阵Μ的获取过程包括:在X/处放置一个点源，获得对应的投 影数据y/，根据点源位置X/和投影数据y/获得扩散分布f(〇)，f(〇)指的是一个投影方向 上，点源在投影数据上的扩散分布，σ为扩散分布f(〇)的参数。通过重复上述过程，得到大量 的点源位置和分别对应的扩散参数0,从而拟合出点源位置与扩散参数的对应关系式，根据 该对应关系式，就能够获得每一个点源位置对应的扩散参数，从而获得该投影方向上的扩 散矩阵，利用PET系统的对称性，根据一个投影方向或者多个投影方向上的扩散矩阵，就能 够得到完整的扩散矩阵M。
[0005] 在重建图像上进行校正的重建方式与之类似，区别在于根据投影数据拟合的是点 源在重建图像上的扩散分布，得到该扩散分布中的扩散参数，之后同样需要得到大量的点 源位置W及分别对应的扩散参数，拟合出点源位置与扩散参数的对应关系式，根据图像重 建时的投影数据重建出PET图像。
[0006] 显然，在上述两种重建方式中，都需要测量大量的点源位置W及分别对应的扩散 参数，W拟合出点源位置与扩散参数的对应关系式，工作量较大。
【发明内容】
[0007] 本发明解决的技术问题在于提供一种PET图像的重建方法及装置，W实现在点扩 散校正时不需要测量大量的点源位置W及分别对应的扩散参数，从而减少工作量。
[000引为此，本发明解决技术问题的技术方案是：
[0009]本发明提供了一种PET图像的重建方法，获取第一扩散分布和第二扩散分布;所述 第一扩散分布为光子垂直入射到晶体正面时投影数据沿着入射方向的垂直方向的扩散分 布，所述第一扩散分布符合高斯分布;所述第二扩散分布为光子垂直入射到晶体侧面时投 影数据沿着入射方向的扩散分布，所述第二扩散分布符合伽马分布;所述方法包括：
[0010] 确定投影方向上的第一晶体和第二晶体的位置；
[0011] 根据所述第一扩散分布、所述第二扩散分布、点源位置W及所述第一晶体的位置， 获得第一扩散表达式，其中，所述第一扩散表达式用于表示投影数据从第一入射点到所述 第一晶体的扩散程度；
[0012] 根据所述第一扩散分布、所述第二扩散分布、点源位置W及所述第二晶体的位置， 获得第二扩散表达式，其中，所述第二扩散表达式用于表示投影数据从第二入射点到所述 第二晶体的扩散程度;所述第一入射点和所述第二入射点分别为点源位置在所述投影方向 上对应的两个入射点；
[0013] 根据所述第一扩散表达式和所述第二扩散表达式，获得扩散关系式，所述扩散关 系式用于表示点源位置与所述投影方向上的投影数据的对应关系；
[0014] 根据所述扩散关系式进行PET图像的重建。
[0015] 可选的，所述获得第一扩散表达式，包括：
[0016] 对所述第一扩散分布和所述第二扩散分布的乘积在第一晶体的区域上进行积分， 获得所述第一扩散表达式，所述第一晶体的区域包括第一区域和第二区域，所述第一区域 为W所述第一入射点为原点，所述第一晶体在光子入射方向的垂直方向的投影区域，所述 第二区域为W所述第一入射点为原点，所述第一晶体在光子入射方向的投影区域；
[0017] 所述获得第二扩散表达式，包括：
[0018] 对所述第一扩散分布和所述第二扩散分布的乘积在第二晶体的区域上进行积分， 获得所述第二扩散表达式，所述第二晶体的区域包括第Ξ区域和第四区域，所述第Ξ区域 为W所述第二入射点为原点，所述第二晶体在光子入射方向的垂直方向的投影区域，所述 第四区域为W所述第二入射点为原点，所述第二晶体在光子入射方向的投影区域。
[0019] 可选的，还包括：
[0020] 获取第Ξ扩散分布，所述第Ξ扩散分布为点源在投影方向的垂直方向上的扩散分 布，所述第Ξ扩散分布符合高斯分布；
[0021] 根据所述第Ξ扩散分布获得第Ξ扩散表达式，所述第Ξ扩散表达式用于表示所述 点源在第五区域的扩散程度，所述第五区域为所述点源在所述投影方向的垂直方向上的扩 散区域；
[0022] 进行PET图像的重建之前，根据所述第Ξ扩散表达式对所述扩散关系式进行校正。
[0023] 可选的，还包括：
[0024] 获取光子对非准直性的概率表达式，所述概率表达式符合高斯分布；
[0025] 根据所述概率表达式获得第四扩散表达式，所述第四扩散表达式用于表示非准直 性对光子入射点产生影响的概率；
[0026] 进行PET图像的重建之前，根据所述第四扩散表达式对所述扩散关系式进行校正。
[0027] 可选的，还包括：
[0028] 获取第五扩散表达式，所述第五扩散表达式通过对>/及积分得到，^?^:为相对于投 影方向的光子发射角度；
[0029] 进行PET图像的重建之前，根据所述第五扩散表达式对所述扩散关系式进行校正。
[0030]本发明提供了一种阳咽像的重建装置，包括：
[0031 ]分布获取单元，用于获取第一扩散分布和第二扩散分布;所述第一扩散分布为光 子垂直入射到晶体正面时投影数据沿着入射方向的垂直方向的扩散分布，所述第一扩散分 布符合高斯分布;所述第二扩散分布为光子垂直入射到晶体侧面时投影数据沿着入射方向 的扩散分布，所述第二扩散分布符合伽马分布；
[0032] 位置确定单元，用于确定投影方向上的第一晶体和第二晶体的位置；
[0033] 表达式获取单元，用于根据所述第一扩散分布、所述第二扩散分布、点源位置W及 所述第一晶体的位置，获得第一扩散表达式，并且根据所述第一扩散分布、所述第二扩散分 布、点源位置W及所述第二晶体的位置，获得第二扩散表达式;其中，所述第一扩散表达式 用于表示投影数据从第一入射点到所述第一晶体的扩散程度，所述第二扩散表达式用于表 示投影数据从第二入射点到所述第二晶体的扩散程度，所述第一入射点和所述第二入射点 分别为点源位置在所述投影方向上对应的两个入射点；
[0034] 关系式获取单元，用于根据所述第一扩散表达式和所述第二扩散表达式，获得扩 散关系式，所述扩散关系式用于表示点源位置与所述投影方向上的投影数据的对应关系；
[0035] 图像重建单元，用于根据所述扩散关系式进行PET图像的重建。
[0036] 可选的，所述表达式获取单元包括：
[0037] 第一获取子单元，用于对所述第一扩散分布和所述第二扩散分布的乘积在第一晶 体的区域上进行积分，获得所述第一扩散表达式，所述第一晶体的区域包括第一区域和第 二区域，所述第一区域为W所述第一入射点为原点，所述第一晶体在光子入射方向的垂直 方向的投影区域，所述第二区域为W所述第一入射点为原点，所述第一晶体在光子入射方 向的投影区域；
[0038] 第二获取子单元，用于对所述第一扩散分布和所述第二扩散分布的乘积在第二晶 体的区域上进行积分，获得所述第二扩散表达式，所述第二晶体的区域包括第Ξ区域和第 四区域，所述第Ξ区域为W所述第二入射点为原点，所述第二晶体在光子入射方向的垂直 方向的投影区域，所述第四区域为W所述第二入射点为原点，所述第二晶体在光子入射方 向的投影区域。
[0039] 可选的，还包括第一校正单元；
[0040] 所述分布获取单元，还用于获取第Ξ扩散分布，所述第Ξ扩散分布为点源在投影 方向的垂直方向上的扩散分布，所述第Ξ扩散分布符合高斯分布；
[0041] 所述表达式获取单元，还用于根据所述第Ξ扩散分布获得第Ξ扩散表达式，所述 第Ξ扩散表达式用于表示所述点源在第五区域的扩散程度，所述第五区域为所述点源在所 述投影方向的垂直方向上的扩散区域；
[0042] 所述第一校正单元，用于所述图像重建单元进行PET图像的重建之前，根据所述第 Ξ扩散表达式对所述扩散关系式进行校正。
[0043] 可选的，还包括第二校正单元；
[0044] 所述表达式获取单元，还用于获取光子对非准直性的概率表达式，所述概率表达 式符合高斯分布，W及根据所述概率表达式获得第四扩散表达式，所述第四扩散表达式用 于表示非准直性对光子入射点产生影响的概率；
[0045] 所述第二校正单元，用于所述图像重建单元进行PET图像的重建之前，根据所述第 四扩散表达式对所述扩散关系式进行校正。
[0046] 可选的，还包括第Ξ校正单元；
[0047] 所述表达式获取单元，还用于获取第五扩散表达式，所述第五扩散表达式通过对 口/