基于pet的图像重建的投影方法和装置制造方法
【专利摘要】本发明公开了一种基于PET的图像重建的投影方法和装置，所述方法包括：在PET环向平面上，沿投影区域确定待投影图像的像素分布图，所述像素分布图包括待投影图像的像素列，所述像素列与所述投影区域平行；对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽；将所述待投影图像在PET轴向平面上投影后，获取像素轴向投影图，所述像素轴向投影图包括像素单元，以及在所述PET轴向平面上，获取所有投影区域的轴向投影图；根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值。与现有技术相比，本发明在保证图像投影准确度的前提下，降低投影的计算量，从而提高图像投影的计算效率。
【专利说明】 基于PET的图像重建的投影方法和装置
【技术领域】
[0001]本发明涉及医疗【技术领域】，具体涉及一种基于PET的图像重建的投影方法和装置。
【背景技术】
[0002]正电子发射计算机断层扫描系统(Positron Emission Tomograph,PET)是当今医疗领域中的一种先进的分子影像学诊断设备，其更是现代重要的医学成像技术手段之一。如图1，图1为PET的结构示意图，其中，基于PET进行医学成像时，PET可以获取当前被扫描对象的重建数据，在完成对象的PET扫描后，可以通过PET获取的重建数据完成图像重建，具体的，图像重建可以分为解析重建和迭代重建，其中，迭代重建应用广泛，但是迭代重建需要反复多次完成投影和反投影的过程。投影过程是指，通过积分的方式，按照一定的规则计算待投影图像像素对任意两个探测器晶体间的投影区域的贡献来模拟真实的探测结果的过程。反投影则是投影的逆运算，是利用投影的结果来得到图像任意像素的值的过程。其中，投影区域可以指两个探测器晶体对应位置连线所围成的三维空间，同时两个探测器晶体对应位置连线为该投影区域的投影方向。由于在迭代重建过程中需要反复多次进行投影和反投影计算，这使得整个重建过程中投影和反投影成为时间消耗最大的部分。所以，提高投影和反投影过程的计算速率对减少迭代重建时间起着关键的作用。
[0003]现有技术中，常见的投影方法线驱动(ray-driven)和像素驱动(pixel-driven)两类，对于线驱动:可以通过计算两探测器晶体间中心连线与某一像素相交的线段长度来计算该像素对这一投影区域的贡献，在这个过程中可以将一个像素划分为几个像素来提高计算精度，但同时会加大计算量从而降低计算速率；对于像素驱动:可以通过判断某一像素是否在投影区域范围之内来判定这一像素对投影是否有贡献，也可以将像素分割为多个子像素来提高计算精度，但是同样会加大计算量。
[0004]综之，现有的投影方法在实现图像重建的过程中不可避免的会消耗大量的时间进行计算。
【发明内容】
[0005]本发明提供了一种基于PET的图像重建的投影方法和装置，在保证图像投影准确度的前提下，降低投影的计算量，从而提高图像投影的计算速率。
[0006]本发明提供了一种基于PET的图像重建的投影方法，所述方法包括:
[0007]在PET环向平面上，沿投影区域确定待投影图像的像素分布图，所述像素分布图包括待投影图像的像素列，所述像素列与所述投影区域平行；
[0008]对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽；
[0009]将所述待投影图像在PET轴向平面上投影后，获取像素轴向投影图，所述像素轴向投影图包括像素单元，以及在所述PET轴向平面上，获取所有投影区域的轴向投影图；[0010]根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值。
[0011]优选地，所述根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值，包括:
[0012]分别获取所述像素轴向投影图的各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值；
[0013]将各个像素单元对应的第一像素投影值相加后，得到所述待投影图像的投影值。
[0014]优选地，所述分别获取所述像素轴向投影图的各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值，包括:
[0015]获取所述所有投影区域的轴向投影图与任一像素单元的相交区域，执行投影流程:获取所述相交区域的轴向中线的长度、所述像素单元的轴向长度和所述像素单元的像素值，获取所述轴向中线的长度与所述像素单元的轴向长度的比值，将所述比值与所述像素值的乘积确定为第一像素投影值，获取所述轴向投影图与任一未执行所述投影流程的像素单元的相交区域，继续执行所述投影流程，直到所有像素单元均完成执行所述投影流程。
[0016]优选地，所述根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值，包括:
[0017]分别获取各个投影区域的轴向投影图与所述像素轴向投影图相交的第二像素投
影值；
[0018]将各个投影区域对应的第二像素投影值相加后，得到所述待投影图像的投影值。
[0019]优选地，所述分别获取各个投影区域的轴向投影图与所述像素轴向投影图相交的第二像素投影值，包括:
[0020]分别获取所述像素轴向投影图与各个投影区域的轴向投影图的相交区域；
[0021]并发执行获取各个相交区域的第二像素投影值的操作。
[0022]优选地，所述对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽，包括:
[0023]对所述像素分布图进行插值重采样，得到重采样后像素分布图，所述重采样后像素分布图的所述像素列的宽度与所述投影区域的宽度等宽。
[0024]本发明还提供一种基于PET的图像重建的投影装置，所述装置包括:
[0025]第一确定模块，用于在PET环向平面上，沿投影区域确定待投影图像的像素分布图，所述像素分布图包括待投影图像的像素列，所述像素列与所述投影区域平行；
[0026]重采样模块，用于对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽；
[0027]第一获取模块，用于将所述待投影图像在PET轴向平面上投影后，获取像素轴向投影图，所述像素轴向投影图包括像素单元；
[0028]第二获取模块，用于在所述PET轴向平面上，获取所有投影区域的轴向投影图；
[0029]第三获取模块，用于根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值。
[0030]优选地，所述第三获取模块，包括:
[0031]第一获取子模块，用于分别获取所述像素轴向投影图的各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值；
[0032]第二获取子模块，用于将各个像素单元对应的第一像素投影值相加后，得到所述待投影图像的投影值。
[0033]优选地，所述第一获取子模块，包括:
[0034]第二获取子模块，用于获取所述所有投影区域的轴向投影图与任一像素单元的相交区域；
[0035]第三获取子模块，用于获取所述相交区域的轴向中线的长度、所述像素单元的轴向长度和所述像素单元的像素值；
[0036]第四获取子模块，用于获取所述轴向中线的长度与所述像素单元的轴向长度的比值；
[0037]确定子模块，用于将所述比值与所述像素值的乘积确定为第一像素投影值；
[0038]第五获取子模块，用于获取所述轴向投影图与任一未触发所述第三获取子模块的像素单元的相交区域，继续触发所述第三获取子模块，直到所有像素单元均完成触发所述第三获取子模块。
[0039]优选地，所述第三获取模块，包括:
[0040]第六获取子模块，用于分别获取各个投影区域的轴向投影图与所述像素轴向投影图相交的第二像素投影值；
[0041]第七获取子模块，用于将各个投影区域对应的第二像素投影值相加后，得到所述待投影图像的投影值。
[0042]优选地，所述第六获取子模块，包括:
[0043]第八获取子模块，用于分别获取所述像素轴向投影图与各个投影区域的轴向投影图的相交区域；
[0044]第九获取子模块，用于并发执行获取各个相交区域的第二像素投影值的操作。
[0045]优选地，所述重采样模块，包括:
[0046]重采样子模块，用于对所述像素分布图进行插值重采样，得到重采样后像素分布图，所述重采样后像素分布图的所述像素列的宽度与所述投影区域的宽度等宽。
[0047]本发明首先在PET环向平面上，沿投影区域确定待投影图像的像素分布图，所述像素分布图包括待投影图像的像素列，所述像素列与所述投影区域平行；其次，对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽；再次，将所述待投影图像在PET轴向平面上投影后，获取像素轴向投影图，所述像素轴向投影图包括像素单元，以及在所述PET轴向平面上，获取所有投影区域的轴向投影图；最后，根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值。与现有技术相比，本发明在保证图像投影准确度的前提下，降低投影的计算量，从而提高图像投影的计算速率。
【专利附图】
【附图说明】
[0048]为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0049]图1为PET的结构示意图；
[0050]图2为实施例一提供的基于PET的图像重建的投影方法流程图；
[0051]图3为探测器间的投影区域与待投影图像的部分像素相交的示意图；
[0052]图4为探测器间的投影区域与待投影图像的部分像素相交在x_y坐标平面内的剖面图；
[0053]图5为探测器间的投影区域与待投影图像的部分像素相交在y-ζ坐标平面内的剖面图；
[0054]图6为待投影图像与部分投影区域相交时，在χ-y坐标平面内的平面图；
[0055]图7为待投影图像与部分投影区域相交时，在X’ _y’坐标平面内的平面图；
[0056]图8为调整位置后的待投影图像与投影区域相交时，在X’ _y’坐标平面内的平面图；
[0057]图9为经过重采样的待投影图像与投影区域相交时，在y’ -Z坐标平面内的平面图；
[0058]图10为探测器间的部分投影区域与调整位置并完成重采样后的待投影图像相交时，在y’-z坐标平面内的平面图；
[0059]图11为实施例三提供的基于PET的图像重建的投影装置结构图。
【具体实施方式】
[0060]下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
[0061]现有技术中PET的图像重建的投影方法较多，而基于体素计算的投影方法具有更高的精确性。但由于计算相对复杂，基于体素的投影方法存在计算效率较低的问题。为了解决这一问题，本发明采用了对待投影图像进行预处理的方法，在保持计算精度的同时，提高了基于体素的投影方法的计算效率。
[0062]实施例一
[0063]参考图2，图2为本实施例提供的基于PET的图像重建的投影方法流程图，具体可以包括:
[0064]步骤201S:在PET环向平面上，沿投影区域确定待投影图像的像素分布图，所述像素分布图包括待投影图像的像素列，所述像素列与所述投影区域平行；
[0065]本实施例中，在PET环向平面上，沿投影区域确定待投影图像的像素分布图，其中，该像素分布图包括待投影图像的像素列，所述沿投影区域即为使该像素列与所述投影区域平行。
[0066]本实施例中可以使用任意一种图像旋转算法确定该待投影图像的像素分布图，使PET环向平面上的待投影图像的像素列与投影方向平行。其中，投影区域可以指两个PET探测器晶体对应位置连线所围成的三维空间，同时两个PET探测器晶体对应位置连线为该投影区域的投影方向。参考图1，PET环向平面是指与PET探测器环所在平面平行的平面。[0067]本实施例中，不再以待投影图像为参考坐标，而是把每个投影区域作为参照系。对应不同方向的投影区域，都要重新确定待投影图像的像素分布图，使确定后的待投影图像的像素列的方向与投影方向平行。因为要计算三维的投影区域与待投影图像的体元相交的体积大小，如果体元与任意角度的投影区域相交，那么得到的相交部分的形状没有可循的规律，这会导致计算投影区域与待投影图像体元的相交体积时计算量巨大。而经过对待投影图像的像素分布图的重新确定，图像的像素列与投影方向平行之后，使在PET环向平面内待投影图像像素与投影区域的相交区域就变换成长方形。这样，相交体积的计算就会简化很多。
[0068]步骤202S:对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽；
[0069]为了使得计算过程更加简化，本实施例对确定后的像素分布图进行重采样操作，以使得该待投影图像的像素列的宽度与投影区域的宽度等宽。
[0070]本实施例中，可以对所述像素分布图进行插值重采样，得到重采样后像素分布图，其中，所述重采样后像素分布图的所述像素列的宽度与所述投影区域的宽度等宽。插值重采样后的待投影图像的像素列数与待投影图像在PET环向平面上相交的投影区域的数量相同。
[0071]经过步骤201S后，在PET环向平面内待投影图像的像素分布图中的像素与任意一个投影区域的相交区域就变换成矩形，但是不同投影区域中包含的矩形的数量和形状不一定相同。为了进一步降低计算的复杂度，本步骤对待投影图像的像素分布图做了重采样操作。将投影区域中，在PET环向平面内的同一行的像素值合并为一个像素值，这个新的像素值为这一行的几个像素值乘以其对应的插值系数后相加得到的。这里的插值系数为像素处于像素分布图行方向某一投影区域内的宽度与这一像素自身宽度的比值。由于计算投影时得到的是各像素值的加权，所以这部分的插值计算不会对计算结果造成影响，保证了计算精度。此时待投影图像的像素分布图与投影区域在PET环向平面内的相交区域变为一列与投影区域等宽的长方形。
[0072]步骤203S:将所述待投影图像在PET轴向平面上投影后，获取像素轴向投影图，所述像素轴向投影图包括像素单元；
[0073]本实施例中，在PET轴向平面上，将所述待投影图像投影，获取像素轴向投影图，其中，该像素轴向投影图包括若干个像素单元，所述像素单元存在与其对应的像素值。
[0074]在PET环向平面上由一列等宽的长方形组成了待投影图像与投影区域的相交区域，由于只是在PET环向平面上把待投影图像和投影区域对齐，在轴向平面上投影区域与待投影图像仍然会斜向相交，所以，本实施例还需要将所述待投影图像在PET轴向平面上投影。
[0075]步骤204S:在所述PET轴向平面上，获取所有投影区域的轴向投影图；
[0076]本实施例中，由于在轴向平面上投影区域与待投影图像仍然会斜向相交，所以，需要在PET轴向平面上也获取所有投影区域的轴向投影图。
[0077]具体的，步骤203S和步骤204S的执行顺序不受限制。
[0078]步骤205S:根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值。[0079]本实施例中，在获取像素轴向投影图和所有投影区域的轴向投影图后，可以根据像素轴向投影图和所有投影区域的轴向投影图，获取待投影图像的投影值。
[0080]实际操作中，获取所述待投影图像的投影值的方法可以包括:
[0081]首先，分别获取所述像素轴向投影图的各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值；
[0082]其次，将各个像素单元对应的第一像素投影值相加后，得到所述待投影图像的投影值。
[0083]本实施例中，首先以该像素轴向投影图的各个像素单元为处理对象，分别获取各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值，其次，将获取的各个第一像素投影值相加后的值，确定为所述待投影图像的投影值。其中，第一像素投影值为各个像素单元分别与投影区域的轴向投影图相交的投影值。
[0084]具体的，分别获取所述像素轴向投影图的各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值，可以具体包括:
[0085]获取所述所有投影区域的轴向投影图与任一像素单元的相交区域，执行投影流程:获取所述相交区域的轴向中线的长度、所述像素单元的轴向长度和所述像素单元的像素值，获取所述轴向中线的长度与所述像素单元的轴向长度的比值，将所述比值与所述像素值的乘积确定为第一像素投影值，获取所述轴向投影图与任一未执行所述投影流程的像素单元的相交区域，继续执行所述投影流程，直到所有像素单元均完成执行所述投影流程。
[0086]实际操作中，将任意一个像素单元确定为处理对象，首先获取该像素单元与所有投影区域的轴向投影图的相交区域，其次，获取该相交区域的轴向中线的长度、像素单元的轴向长度和所述像素单元的像素值，最后，获取所述轴向中线的长度与所述像素单元的轴向长度的比值，将所述比值与所述像素值的乘积确定为第一像素投影值。将该像素轴向投影图的所有像素单元对应的第一像素投影值均如上述方法获取到。其中，轴向中线为沿PET轴线方向的中线，像素单元的轴向长度也为该像素单元沿PET轴线方向的边长。另外，各个第一像素投影值的获取顺序不做限制，但是并发执行获取各个第一像素投影值的操作可以提高计算效率
[0087]实际操作中，获取所述待投影图像的投影值的方法还可以包括:
[0088]首先，分别获取各个投影区域的轴向投影图与所述像素轴向投影图相交的第二像素投影值；
[0089]其次，将各个投影区域对应的第二像素投影值相加后，得到所述待投影图像的投影值。
[0090]本实施例中，首先确定投影区域为处理对象，获取各个投影区域的轴向投影图与像素轴向投影图相交的第二像素投影值，其次，将获取的各个第二像素投影值相加后的值，确定为所述待投影图像的投影值，其中第二像素投影值为各个投影区域分别与像素轴向投影图相交的投影值。
[0091]具体的，分别获取各个投影区域的轴向投影图与所述像素轴向投影图相交的第二像素投影值，可以包括:
[0092]首先，分别获取所述像素轴向投影图与各个投影区域的轴向投影图的相交区域；
[0093]其次，并发执行获取各个相交区域的第二像素投影值的操作。[0094]实际操作中，将投影区域确定为处理对象，首先获取投影区域的轴向投影图与像素轴向投影图的相交区域，其次，获取该相交区域对应的第二像素投影值。具体的，可以将所有投影区域对应的第二像素投影值均如上述方法获取到，同时，各个第二像素投影值的获取顺序不做限制，但是本实施例并发执行获取各个相交区域的第二像素投影值的操作可以明显提高计算效率。
[0095]本实施例中，投影区域与待投影图像的相交体积的形状都为直棱柱，所以要得到相交体积的大小，需先计算直棱柱一个底面的面积，即投影区域与待投影图像的像素分布图的像素单元相交形成的各种图形的面积。为了简化计算，这里认为沿着轴向穿过各像素的中线也是这些图形的中线。在得到这些图形的中线长度后，就能够算出投影区域与待投影图像像素的相交体积的值。
[0096]实际操作中，由于计算得到的投影值是投影区域内各像素值的加权，而计算相交体积的目的是确定加权系数，加权系数为每个像素与某一投影区域的相交体积与该像素体积的比值。经过步骤201S和步骤202S之后，像素与投影区域的相交体积变为与像素等高的直棱柱，所以体积的比值可直接用轴向投影与待投影图像像素相交面积与像素面积的比值得到。又因为这部分相交面积与像素面积和它们的轴向中线长度成正比，最后相交体积的比值转化为了轴向投影与待投影图像像素相交图形的中线长度与像素轴向长度的比值，将这一比值作为各自像素的加权系数，与计算投影区域中的每一个像素相乘并累加，从而得到投影值。
[0097]考虑到计算的复杂度，这里没有直接计算每个相交体积的值，而是先确定待投影图像的像素分布图，然后对其进行重采样，生成一个新的待投影图像。在这个基础之上进行投影计算时，在PET环向平面上就可以看作对某一投影区域有贡献的像素就是与它重合的那一列像素，从而省去了多次计算投影与像素相交边界的问题，直接通过投影内像素轴向中线和这个像素的归一化后的相交长度与像素值的乘积就能得到投影值，简化了计算，起到提高计算速率的作用。
[0098]与此同时，在确定像素分布图和重采样过程中，仍然保留了原图像所带的信息，并将所有的像素的作用都加入到了投影结果中，从而避免了在现有技术的线驱动和像素驱动中由于对像素贡献累计不全面而引入高频噪声的结果。
[0099]本实施例首先在PET环向平面上，沿投影区域确定待投影图像的像素分布图，所述像素分布图包括待投影图像的像素列，所述像素列与所述投影区域平行；其次，对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽；再次，将所述待投影图像在PET轴向平面上投影后，获取像素轴向投影图，所述像素轴向投影图包括像素单元，以及在所述PET轴向平面上，获取所有投影区域的轴向投影图；最后，根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值。与现有技术相t匕，本实施例在保证图像投影准确度的前提下，降低投影的计算量，从而提高图像投影的计算效率。
[0100]实施例二
[0101]如图3，其中301是两个探测器晶体间的投影区域，302是由方形像素按照行和列排列构成的三维待投影图像。图4中的投影301a是图3中投影区域301在χ-y平面的剖面图，图5中的投影301b是图3中投影区域301在y-z平面内的剖面图。如要计算投影区域301的投影值，则需要将投影区域301中的每个像素的体积逐一累加。
[0102]其中，投影计算公式如式(I)，
【权利要求】
1.一种基于PET的图像重建的投影方法，其特征在于，所述方法包括: 在PET环向平面上，沿投影区域确定待投影图像的像素分布图，所述像素分布图包括待投影图像的像素列，所述像素列与所述投影区域平行； 对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽；将所述待投影图像在PET轴向平面上投影后，获取像素轴向投影图，所述像素轴向投影图包括像素单元，以及在所述PET轴向平面上，获取所有投影区域的轴向投影图； 根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值。
2.根据权利要求1所述的方法，其特征在于，所述根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值，包括: 分别获取所述像素轴向投影图的各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值； 将各个像素单元对应的第一像素投影值相加后，得到所述待投影图像的投影值。
3.根据权利要求2所述的方法，其特征在于，所述分别获取所述像素轴向投影图的各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值，包括: 获取所述所有投影区域的轴向投影图与任一像素单元的相交区域，执行投影流程:获取所述相交区域的轴向中线的长度、所述像素单元的轴向长度和所述像素单元的像素值，获取所述轴向中线的长度与所述像素单元的轴向长度的比值，将所述比值与所述像素值的乘积确定为第一像素投影值·，获取所述轴向投影图与任一未执行所述投影流程的像素单元的相交区域，继续执行所述投影流程，直到所有像素单元均完成执行所述投影流程。
4.根据权利要求1所述的方法，其特征在于，所述根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值，包括: 分别获取各个投影区域的轴向投影图与所述像素轴向投影图相交的第二像素投影值； 将各个投影区域对应的第二像素投影值相加后，得到所述待投影图像的投影值。
5.根据权利要求4所述的方法，其特征在于，所述分别获取各个投影区域的轴向投影图与所述像素轴向投影图相交的第二像素投影值，包括: 分别获取所述像素轴向投影图与各个投影区域的轴向投影图的相交区域； 并发执行获取各个相交区域的第二像素投影值的操作。
6.根据权利要求1所述的方法，其特征在于，所述对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽，包括: 对所述像素分布图进行插值重采样，得到重采样后像素分布图，所述重采样后像素分布图的所述像素列的宽度与所述投影区域的宽度等宽。
7.一种基于PET的图像重建的投影装置，其特征在于，所述装置包括: 第一确定模块，用于在PET环向平面上，沿投影区域确定待投影图像的像素分布图，所述像素分布图包括待投影图像的像素列，所述像素列与所述投影区域平行； 重采样模块，用于对所述像素分布图进行重采样，使得所述像素列的宽度与所述投影区域的宽度等宽； 第一获取模块，用于将所述待投影图像在PET轴向平面上投影后，获取像素轴向投影图，所述像素轴向投影图包括像素单元； 第二获取模块，用于在所述PET轴向平面上，获取所有投影区域的轴向投影图； 第三获取模块，用于根据所述像素轴向投影图和所有投影区域的轴向投影图，获取所述待投影图像的投影值。
8.根据权利要求7所述的装置，其特征在于，所述第三获取模块，包括: 第一获取子模块，用于分别获取所述像素轴向投影图的各个像素单元与所有投影区域的轴向投影图相交的第一像素投影值； 第二获取子模块，用于将各个像素单元对应的第一像素投影值相加后，得到所述待投影图像的投影值。
9.根据权利要求8所述的装置，其特征在于，所述第一获取子模块，包括: 第二获取子模块，用于获取所述所有投影区域的轴向投影图与任一像素单元的相交区域； 第三获取子模块，用于获取所述相交区域的轴向中线的长度、所述像素单元的轴向长度和所述像素单元的像素值； 第四获取子模块，用于获取所述轴向中线的长度与所述像素单元的轴向长度的比值； 确定子模块，用于将所述比值与所述像素值的乘积确定为第一像素投影值； 第五获取子模块，用于获取所述轴向投影图与任一未触发所述第三获取子模块的像素单元的相交区域，继续触发所述第三获取子模块，直到所有像素单元均完成触发所述第三获取子模块。
10.根据权利要求7所述的装置，其特征在于，所述第三获取模块，包括: 第六获取子模块，用于分别获取各个投影区域的轴向投影图与所述像素轴向投影图相交的第二像素投影值； 第七获取子模块，用于将各个投影区域对应的第二像素投影值相加后，得到所述待投影图像的投影值。
11.根据权利要求10所述的装置，其特征在于，所述第六获取子模块，包括: 第八获取子模块，用于分别获取所述像素轴向投影图与各个投影区域的轴向投影图的相交区域； 第九获取子模块，用于并发执行获取各个相交区域的第二像素投影值的操作。
12.根据权利要求7所述的装置，其特征在于，所述重采样模块，包括: 重采样子模块，用于对所述像素分布图进行插值重采样，得到重采样后像素分布图，所述重采样后像 素分布图的所述像素列的宽度与所述投影区域的宽度等宽。
【文档编号】G06T17/00GK103593868SQ201310476912
【公开日】2014年2月19日 申请日期:2013年10月12日 优先权日:2013年10月12日
【发明者】陈锡民, 韩冬, 王军博, 杨金柱, 赵大哲 申请人:沈阳东软医疗系统有限公司