本公开涉及医疗成像技术，特别涉及一种门控信号的确定方法和装置。
背景技术：
在采用pet(positronemissioncomputedtomography，正电子发射型计算机断层显像)设备对患者某些身体部位扫描成像时，成像的图像质量容易受到患者呼吸运动的影响。例如，在对患者的肺部扫描成像时，呼吸运动将使得肺部发生周期性地扩张收缩运动，而pet的每一床位的扫描时间相对呼吸周期来说比较长(约1～2min)，导致测得的数据是肺部发生多次呼吸周期的运动后的叠加结果，从而使得在pet重建图像时肺部边界变模糊，产生图像伪影。而且，如果肺部有病灶点时，周期性的呼吸运动会使得病灶点在重建图像上尺寸扩大，病灶点的对比度下降，影响诊断精度。
为了消减上述由呼吸运动造成的图像伪影等成像问题，通常需要采用外置的呼吸门控设备采集患者呼吸运动的信号以获得呼吸周期，并根据此呼吸周期将测得的数据划分成不同的呼吸期相，再获得不同期相下的重建图像即可减轻上述成像问题。但是，外置的呼吸门控设备一般价格昂贵且操作也相对较复杂。
技术实现要素：
有鉴于此，本公开提供一种门控信号的确定方法和装置，以实现在不借助外置呼吸门控设备的情况下获取门控信号，降低成本。
具体地，本公开是通过如下技术方案实现的：
第一方面，提供一种门控信号的确定方法，所述方法包括：
采集正电子发射型计算机断层显像pet扫描在检测时间的每个子时段中得到的投影数据，所述投影数据包括：飞行时间tof信息；
根据包括所述tof信息的投影数据，确定湮灭点位置信息，所述湮灭点位置信息用于表示湮灭点的三维空间坐标；
根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置；
根据所述各个子时段对应的湮灭点重心位置，得到门控信号。
第二方面，提供一种门控信号的确定装置，所述装置包括：
数据采集模块，用于采集正电子发射型计算机断层显像pet扫描在检测时间的每个子时段中得到的投影数据，所述投影数据包括：飞行时间tof信息；
位置确定模块，用于根据包括所述tof信息的投影数据，确定湮灭点位置信息，所述湮灭点位置信息用于表示湮灭点的三维空间坐标；
重心处理模块，用于根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置；
门控确定模块，用于根据所述各个子时段对应的湮灭点重心位置，得到门控信号。
第三方面，提供一种门控信号的确定设备，所述设备包括：存储器、处理器，以及存储在存储器上并可在处理器上运行的计算机指令，所述处理器执行指令时实现以下步骤：
采集正电子发射型计算机断层显像pet扫描在检测时间的每个子时段中得到的投影数据，所述投影数据包括：飞行时间tof信息；
根据包括所述tof信息的投影数据，确定湮灭点位置信息，所述湮灭点位置信息用于表示湮灭点的三维空间坐标；
根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置；
根据所述各个子时段对应的湮灭点重心位置，得到门控信号。
第四方面，提供一种计算机可读存储介质，其具有存储在其上的指令，当由一个或多个处理器执行时，使得一个或多个处理器执行门控信号的确定方法，该方法包括：
采集正电子发射型计算机断层显像pet扫描在检测时间的每个子时段中得到的投影数据，所述投影数据包括：飞行时间tof信息；
根据包括所述tof信息的投影数据，确定湮灭点位置信息，所述湮灭点位置信息用于表示湮灭点的三维空间坐标；
根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置；
根据所述各个子时段对应的湮灭点重心位置，得到门控信号。
本公开提供的门控信号的确定方法和装置，根据采集的pet投影数据确定湮灭点位置，并根据各个子时段对应的湮灭点位置确定门控信号，使得可以不再借助外置呼吸门控设备，直接通过处理投影数据就可以获得门控信号，节省了成本；并且，该方法通过依据tof信息，可以更加准确的定位湮灭点的位置，进而可以使得湮灭点重心位置和门控信号的确定也更加准确。
附图说明
图1是本公开一示例性实施例示出的一种pet/ct设备示意图；
图2是本公开一示例性实施例示出的一种呼吸曲线的示意图；
图3是本公开一示例性实施例示出的一种呼吸运动的示意图；
图4是本公开一示例性实施例示出的一种门控信号的确定方法的流程图；
图5是本公开一示例性实施例示出的一种湮灭点位置示意图；
图6是本公开一示例性实施例示出的一种湮灭点重心位置的计算示意图；
图7是本公开一示例性实施例示出的一种门控信号的确定装置。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
本公开的例子旨在提供一种门控信号的确定方法，以替代外置的呼吸门控设备，节省门控信号的确定成本。
在一个例子中，参见图1所示，本公开例子的门控信号确定方法可以由pet/ct设备11执行，比如，该pet/ct设备11可以包括一个门控信号确定设备12。当然，该门控信号确定设备12不属于外置的呼吸门控设备，而是属于pet/ct设备11自身的一部分，即该门控信号确定设备12不需要绑在患者身体上采集患者的呼吸运动，而是可以集成在pet/ct设备11中，并且可以根据采集到的pet扫描的投影数据，获取到门控信号。
图1只是示例一种结构，实际实施中并不局限于这种结构。比如，还可以是将图像重建设备(未示出，例如，可以是连接图1中的pet/ct设备11的一个计算机装置)集成门控信号确定的功能，使得图像重建设备可以执行本公开例子的门控信号确定方法，根据采集到的pet投影数据处理得到门控信号。并且，图像重建设备可以继续根据确定的门控信号将投影数据划分为不同的呼吸期相，按照呼吸期相进行图像重建。
在描述本公开例子的门控信号的确定方法之前，先简单说明下“门控信号”：
门控信号：患者的呼吸运动将使得胸肺部发生周期性地扩张收缩运动，而用于表示这种呼吸运动的信号可以称为门控信号。比如，当患者处于平躺状态时，其体表的某部位例如胸肺部的高度将随着时间进行变化，吸气时胸肺部上浮，呼气时胸肺部下降。图2示例了一个呼吸曲线，图中可以看出一个呼吸周期包括一次吸气和一次呼气，不同的患者可以具有不同的呼吸周期、不同的呼吸频率。该呼吸曲线表示了一种呼吸运动，可以称为门控信号。
呼吸期相：可以根据门控信号反映出的呼吸周期，将不同呼吸周期中处于相同的运动阶段的pet投影数据划分在一起，组成一个呼吸期相。一个呼吸期相可以对应呼吸运动的一个运动阶段。如图2所示，假设可以将一个呼吸周期的时长分为10等份，将其中的第2份称为一个目标时段21，那么每个呼吸周期中都存在一个对应的目标时段21，参见图2中打斜剖线的部分。可以将每个呼吸周期中的目标时段21中采集到的pet投影数据组成一个组，这个组即对应一个呼吸期相的数据，这个呼吸期相的数据对应着相同的运动阶段。图2中示例的门控信号，是以划分为十个呼吸期相为例，实际实施中不限制于此。根据一个呼吸期相的数据可以重建得到对应某个运动阶段的图像，减轻图像伪影，图2的例子中可以重建得到分别对应十个呼吸期相的图像。
如下将描述如何确定一个门控信号，例如，当采用pet设备对一位患者进行pet扫描时，pet设备可以根据扫描过程中采集到的投影数据，通过执行本公开例子的门控信号的确定方法，获得到对应该患者的门控信号，了解到该患者的呼吸运动的情况。参见图3的示意，图3示意了患者的胸肺部31的起伏运动，当患者呼吸时，随着患者的呼吸运动，胸肺部31也在不断的起伏，当患者吸气时，胸肺部31扩张到图3中的高处；当患者呼气时，胸肺部31收缩到图3中的低处。此外，患者的胸肺部31中可以有多个湮灭点32，由图3可以看到，随着胸肺部31的起伏，这些湮灭点32的位置也在随之移动。由此可以得到，如果确定了湮灭点的位置，并且知晓对应位置的时间点，那么可以根据湮灭点的运动信息反映出患者的呼吸运动的情况。
图4示例了本公开一个例子的门控信号的确定方法，可以包括：
在步骤400中，采集正电子发射型计算机断层显像pet扫描在检测时间的子时段中得到的投影数据，所述投影数据包括：飞行时间tof信息。
本步骤中，可以利用tof-pet设备对患者进行扫描，并采集在检测时间内的pet投影数据。可以将检测时间分为多个子时段，每个子时段的预设时长可以自主设定，比如，可以设为0.5s。如下表1示例：
表1投影数据
例如，每一条投影数据可以对应一条lor(lineofresponse，响应线)，投影数据中可以包括：该lor对应的两个晶体的晶体标识、每个晶体接收到光子的接收时间。
在步骤402中，根据包括所述tof信息的投影数据，确定湮灭点位置信息，所述湮灭点位置信息用于表示湮灭点的三维空间坐标。
本步骤中，可以根据包括tof信息的投影数据，确定湮灭点位置信息。位置的确定可以按照常规方法进行，不再详述，例如，根据湮灭事件的两个光子的接收时间的时间差，可以大致确定发生湮灭事件的湮灭点在lor上的位置。
湮灭点位置信息用于表示湮灭点的三维空间坐标，即该位置信息可以确定出湮灭点在三维空间中的具体位置，是湮灭点位置的相对准确的定位。湮灭点位置信息可以用笛卡尔坐标系表示，或者也可以用极坐标系表示。
如图5的示意，图5示意了一个lor51，还示意了在笛卡尔坐标系下该lor51上的一个湮灭点52，该湮灭点52的坐标(x，y，z)，其中，z表示湮灭点的轴向位置坐标，(x，y)表示湮灭点的径向位置坐标。当使用笛卡尔坐标系表示湮灭点位置信息时，接收湮灭点发出的光子的各个晶体的坐标也使用笛卡尔坐标系表示。此外，极坐标系与笛卡尔坐标系可以互相转换，比如，极坐标下的方位角、仰角和距离可以先转换成笛卡尔坐标系下的x、y和z坐标值再进行相应的计算。后续的部分例子描述中，以笛卡尔坐标系为例，极坐标系同理。
在步骤404中，根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置。
在如上的步骤402中可以确定各个湮灭点的湮灭点位置信息，例如确定出每个湮灭点的笛卡尔坐标(x，y，z)，在此基础上，本步骤可以确定这些湮灭点的湮灭点重心位置。比如，图3中示意的胸肺部31中包括的各个湮灭点，这些分散的湮灭点可以具有一个重心位置，用以代表这些湮灭点簇。
此外，在确定湮灭点重心位置时，可以分别确定各个子时段对应的湮灭点重心位置。比如，子时段1对应一个湮灭点重心位置，子时段2对应另一湮灭点重心位置。不同子时段可以对应着不同的湮灭点重心位置，从而可以反映出湮灭点重心位置的移动。
在步骤406中，根据各个子时段对应的湮灭点重心位置，得到门控信号。
本步骤可以根据各个子时段对应的湮灭点重心位置，确定出湮灭点重心位置随着时间的变化，这种变化是随着患者的呼吸运动形成，能够反映出患者的呼吸运动的情况。比如，可以将各个子时段作为横轴，将每个子时段对应的湮灭点重心位置作为纵轴，并通过插值和数据拟合处理得到门控信号。
本例子的门控信号的确定方法，通过根据采集的pet投影数据确定湮灭点位置，并根据各个子时段对应的湮灭点位置确定门控信号，使得可以不再借助外置呼吸门控设备，直接通过处理投影数据就可以获得门控信号，节省了成本；并且，该方法通过依据tof信息，可以更加准确的定位湮灭点的位置，进而可以使得湮灭点重心位置和门控信号的确定也更加准确。
在一个例子中，以笛卡尔坐标系下的湮灭点坐标为例，说明湮灭点重心位置的计算方式。如图6所示，该方法可以包括如下处理，但是并不局限于如下的执行顺序，比如，步骤604和606可以同时执行。
在步骤600中，将相邻两个时间间隔的湮灭点位置矩阵做差，得到位置统计范围。
本步骤中，各个湮灭点的湮灭点位置信息可以组成一个湮灭点位置矩阵，该湮灭点位置矩阵相当于一个三维立体的湮灭点空间分布。所述的相邻两个时间间隔，其中的每一个时间间隔可以包括至少一个0.5s的子时段，比如，一个时间间隔可以是一个0.5s的子时段，或者也可以是两个子时段0.5*2＝1s，或者还可以是三个子时段3*0.5＝1.5s，等。将两个相邻的时间间隔中确定的湮灭点位置矩阵计算差值，可以根据差值的变化范围来确定周期性运动所发生的区域。
对于相邻的两个时间间隔内的湮灭点位置矩阵，只有因为发生呼吸运动及心脏搏动产生的空间位移变化，其他位置不变或变化很小，基本可以忽略。所以，对两个湮灭点位置矩阵计算差值后只在运动区域有相对明显的值变化，其他区域变化很小可近似忽略。本例子可以选取有相对明显变化的区域作为位置统计范围，在计算湮灭点重心位置时只计算该位置统计范围内的湮灭点数目及位置分布，而不用将所有位置的湮灭点都参与计算，因此可以显著降低计算量。当然，在其他例子中，也可以不选取位置统计范围，而是将全部的湮灭点位置信息参与重心位置的计算。此外，在具体确定位置统计范围时，可以将上述的有相对明显变化的区域沿轴向上下扩展若干层(比如3层，其中每个探测器环可视为1层)后，确定为位置统计范围。
在步骤602中，对于位置统计范围内的每个坐标点，统计对应所述坐标点的位置处，在所述子时段内总共的湮灭点数目。
例如，位置统计范围内可以包括很多的坐标点，每一个坐标点都有对应的坐标。以位置统计范围内的其中一个坐标点i为例，其坐标(xi，yi，zi)，其中的任一项都可以称为一个坐标值，比如，xi可以称为一个坐标值，yi也可以称为一个坐标值。本步骤可以统计在某个子时段t内，在坐标点i的位置处发生的湮灭点数目，即这些湮灭点数目的湮灭点都发生在坐标点i的位置处，比如，五个湮灭点都位于坐标点i的位置。可以用num(t，i)表示在子时段t内于坐标点i的位置处总共的湮灭点数目。
本步骤可以分别统计在各个子时段内的各个坐标点位置处的湮灭点数目。
在步骤604中，计算第一加和值，所述第一加和值是位置统计范围内的各个坐标点的坐标值与对应所述坐标点的湮灭点数目进行加权求和。
例如，所述坐标值可以是：x坐标、或者y坐标、或者z坐标。
以x坐标为例，可以按照如下公式(1)计算第一加和值：
根据公式(1)，可以将位置统计范围内的各个坐标点i处的湮灭点数目num(t，i)，与该坐标点i的x坐标i(xi)进行加权求和，得到第一加和值。
同理，还可以将各个坐标点i处的湮灭点数目与该坐标点i的y坐标进行加权求和，得到的
在步骤606中，计算第二加和值，所述第二加和值是位置统计范围内的各个坐标点处对应的湮灭点数目进行求和。
例如，本步骤的第二加和值按照如下公式(2)计算：
根据上述的公式(2)，以某个子时段t为例，可以统计位置统计范围内的各个坐标点i处在该子时段t内的湮灭点数目，并将各个坐标点i处对应的湮灭点数目进行求和，得到第二加和值。
在步骤608中，将第一加和值第二加和值的比值，作为位置统计范围内的湮灭点重心位置的坐标值。
所述坐标值，包括如下任一项：x坐标、或者y坐标、或者z坐标。
如下计算子时段t对应的位置统计范围内的湮灭点重心位置：
根据上述的公式(3)，w(x,t)是子时段t对应的湮灭点重心位置w的x坐标，w(y,t)是子时段t对应的湮灭点重心位置w的y坐标，w(z,t)是子时段t对应的湮灭点重心位置w的z坐标。根据公式(3)可以计算得到各个子时段分别对应的湮灭点重心位置的x、y、z三个坐标轴的坐标值。
此外，在确定各个子时段的湮灭点重心位置后，可以将各个子时段作为横轴，将每个子时段的湮灭点重心位置作为纵轴，得到湮灭点的重心位置-时间函数，作为门控信号。
在另一个例子中，以患者的胸肺部的运动为例，实际上，患者身体部位的移动，不仅包括呼吸运动的影响，还包括心脏搏动的影响。例如，呼吸运动可以影响湮灭点的轴向位置分布，而心脏搏动主要影响湮灭点的径向位置分布，胸肺部的运动是轴向运动和径向运动的综合结果。但是在图像重建时，划分呼吸期相是想要依据呼吸运动的运动规律来划分，比如，在对肺部成像时，主要由于患者的呼吸运动使得肺部发生周期性的扩张收缩运动，呼吸运动对肺部的成像质量起主要影响作用。
基于上述，本公开的例子提供的门控信号的确定方法，还可以区分出呼吸运动的门控信号和心脏搏动的门控信号，其中，呼吸运动的门控信号可以反映患者的呼吸运动情况，心脏搏动的门控信号可以反映患者的心脏搏动情况。并且，只基于呼吸运动的门控信号进行呼吸期相的划分，这样将使得呼吸期相的划分更加准确，从而图像重建的质量也更好。例如，可以在确定呼吸运动门控信号和心脏搏动门控信号之后，根据呼吸运动门控信号，将投影数据划分多个呼吸期相，并分别重建每一个呼吸期相对应的重建图像。
区分上述的呼吸运动的门控信号和心脏搏动的门控信号的方式，可以包括：基于各个子时段对应的湮灭点重心位置，得到各个子时段对应的湮灭点轴向位置、以及各个子时段对应的湮灭点径向位置。根据所述各个子时段对应的湮灭点轴向位置，得到轴向门控信号。根据所述各个子时段对应的湮灭点径向位置，得到径向门控信号。例如，以笛卡尔坐标系表示湮灭点位置信息为例，可以按照如下方式区分呼吸运动门控信号和心脏搏动门控信号：
当以笛卡尔坐标系表示时，湮灭点轴向位置可以是所述湮灭点重心位置的z坐标；所述湮灭点径向位置可以是所述湮灭点重心位置的(x,y)坐标。
所述的根据各个子时段的湮灭点轴向位置得到轴向门控信号，可以包括：将各个子时段对应的时间依时间顺序作为横轴，将每个子时段的湮灭点重心位置的z坐标作为纵轴，得到湮灭点的轴向重心位置-时间函数。比如，横轴是各个子时段t，纵轴可以是对应t的湮灭点重心位置w的z坐标w(z,t)。
所述的根据各个子时段的湮灭点径向位置得到径向门控信号，可以包括：将各个子时段对应的时间依时间顺序作为横轴，将每个子时段的湮灭点重心位置的(x,y)坐标作为纵轴，得到湮灭点的径向重心位置-时间函数。比如，横轴是各个子时段t，纵轴可以是对应t的湮灭点重心位置w的径向坐标[w(x,t)，w(y,t)]。
还可以根据所述轴向重心位置-时间函数和径向重心位置-时间函数的周期性变化，确定呼吸运动门控信号和心脏搏动门控信号。例如，可以分析上述两个函数的周期性变化规律，来判断确定是由呼吸运动(通常情况下平均周期约5s)还是心脏搏动(通常情况下平均周期约1s)引起的周期性变化，即而确定对应的门控信号。
上述方式也相当于将湮灭点重心分布-时间函数拆分成了沿轴向及径向变化的两个函数，即轴向重心位置-时间函数和径向重心位置-时间函数。
此外，对于上文中获得的湮灭点重心位置-时间函数，还可以不对其沿轴向及径向拆分成两个时间变化函数，而是可以对其进行频域变换。又例如，还可以按照如下方式区分呼吸运动门控信号和心脏搏动门控信号：
在得到湮灭点的重心位置-时间函数作为门控信号之后，可以将湮灭点重心位置-时间函数进行频域变换，得到对应的频谱。再基于具有预设参数的滤波器对所述频谱进行滤波，所述预设参数包括：对应呼吸运动的滤波参数和对应心脏搏动的滤波参数。比如，可以设置滤波器的中心频率是呼吸频率0.5hz，用于过滤呼吸运动信号；又比如，可以设置滤波器的中心频率是心跳频率1hz，以过滤得到心脏搏动信号。接着可以对滤波后的频谱进行频域逆变换，得到呼吸运动门控信号和心脏搏动门控信号。
本例子的门控信号的确定方法，通过利用tof信息，确定每一条lor上的湮灭点的较为准确的发生位置，使得对湮灭点位置的计算更加精确，所获得的门控信号也更准确。门控信号的准确也使得对呼吸期相的划分更准确，从而根据呼吸期相重建的图像质量更好。
为了实现上述的门控信号的确定方法，本公开还提供了一种门控信号的确定装置，如图7所示，该装置可以包括：数据采集模块71、位置确定模块72、重心处理模块73和门控确定模块74。
数据采集模块71，用于采集正电子发射型计算机断层显像pet扫描在检测时间的每个子时段中得到的投影数据，所述投影数据包括：飞行时间tof信息；
位置确定模块72，用于根据包括所述tof信息的投影数据，确定湮灭点位置信息，所述湮灭点位置信息用于表示湮灭点的三维空间坐标；
重心处理模块73，用于根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置；
门控确定模块74，用于根据所述各个子时段对应的湮灭点重心位置，得到门控信号。
在一个例子中，重心处理模块73，具体用于：当所述湮灭点位置信息采用笛卡尔坐标系表示时，所述湮灭点位置信息包括：湮灭点的轴向位置坐标z、以及径向位置坐标(x，y)；按照如下方式分别确定每个子时段的湮灭点重心位置：
对于位置统计范围内的每个湮灭点，统计对应所述湮灭点的位置处，在所述子时段内总共的湮灭点数目；
将第一加和值第二加和值的比值，作为位置统计范围内的湮灭点重心位置的坐标值，所述第一加和值是所述位置统计范围内的各个坐标点的坐标值与所述坐标点处的湮灭点数目进行加权求和，所述第二加和值是所述各个坐标点处对应的湮灭点数目进行求和；
所述坐标值，包括如下任一项：x坐标、或者y坐标、或者z坐标。
在一个例子中，所述重心处理模块73，还用于：根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置之前，按照如下方式确定所述位置统计范围：
将相邻两个时间间隔的湮灭点位置矩阵做差，得到所述位置统计范围；
所述时间间隔中包括至少一个所述子时段；所述湮灭点位置矩阵包括各个湮灭点的湮灭点位置信息。
在一个例子中，所述门控确定模块74，具体用于：将各个子时段作为横轴，将每个子时段的湮灭点重心位置作为纵轴，得到湮灭点的重心位置-时间函数，作为所述门控信号。
在一个例子中，门控确定模块74，具体用于：基于所述各个子时段对应的湮灭点重心位置，得到各个子时段对应的湮灭点轴向位置、以及各个子时段对应的湮灭点径向位置；根据所述各个子时段对应的湮灭点轴向位置，得到轴向门控信号；根据所述各个子时段对应的湮灭点径向位置，得到径向门控信号。
本公开的方法的功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本公开的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台图像处理设备执行本公开各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。
本例子提供的一种计算机可读存储介质，其具有存储在其上的指令，当由一个或多个处理器执行时，使得一个或多个处理器执行门控信号的确定方法，该方法包括：
采集正电子发射型计算机断层显像pet扫描在检测时间的每个子时段中得到的投影数据，所述投影数据包括：飞行时间tof信息；
根据包括所述tof信息的投影数据，确定湮灭点位置信息，所述湮灭点位置信息用于表示湮灭点的三维空间坐标；
根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置；
根据所述各个子时段对应的湮灭点重心位置，得到门控信号。
本公开还提供了一种门控信号的确定设备，该设备可以包括：存储器、处理器，以及存储在存储器上并可在处理器上运行的计算机指令，所述处理器执行指令时实现以下步骤：
采集正电子发射型计算机断层显像pet扫描在检测时间的每个子时段中得到的投影数据，所述投影数据包括：飞行时间tof信息；
根据包括所述tof信息的投影数据，确定湮灭点位置信息，所述湮灭点位置信息用于表示湮灭点的三维空间坐标；
根据所述湮灭点位置信息，分别确定各个子时段对应的湮灭点重心位置；
根据所述各个子时段对应的湮灭点重心位置，得到门控信号。
以上所述仅为本公开的较佳实施例而已，并不用以限制本公开，凡在本公开的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本公开保护的范围之内。