本发明涉及医学成像领域，特别是涉及一种动态pet参数成像方法、装置、系统和计算机可读存储介质。
背景技术
正电子发射断层扫描(positronemissiontomography，pet)是核医学领域中重要的一种临床检查影像技术，现已被广泛应用于医学领域的诊断和研究。pet采用正电子核素作为示踪剂，通过病灶部位对示踪剂的摄取了解病灶功能代谢状态，从而对疾病做出诊断。通常，pet扫描仪会集成有电子计算机断层扫描(computedtomography，ct)，通过配准的ct图像获得衰减图像，用于对重建的图像进行衰减校正。动态pet成像不仅能显示放射性示踪剂浓度的空间分布，还揭示了示踪剂随时间代谢的动态过程。通过应用动力学模型，即可获得具有丰富意义的动力学参数，从而为疾病诊断、治疗以及药物开发提供指导。
在进行动态pet成像时，需要时间活度曲线的信息。传统技术中，一般通过图像分割、人工智能等方法获取主动脉的空间位置，然后通过计算pet数据在各时间点的重建帧对时间活度曲线进行拟合，并且，在获取时间活度曲线的信息的过程中，需要受检者从药物注射开始的整个扫描和采集过程中保持不动。
但是，传统技术中在获取时间活度曲线的信息的过程中，整个扫描过程均需要对受检者进行ct扫描，辐射比较大；并且，需要受检者从药物注射开始的整个扫描和采集过程中保持不动，临床操作难度比较大。
技术实现要素：
基于此，有必要针对传统技术中获取时间活度曲线时，辐射较大，操作难度较大的问题，提供一种动态pet参数成像方法、装置、系统和计算机可读存储介质。
第一方面，本发明实施例提供一种动态pet参数成像方法，包括：
获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；
间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；
根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线；
根据所述完整时间活度曲线进行动态pet参数成像。
本实施例提供的动态pet参数成像方法，计算机设备可以获取受检者在第一时长内的第一pet数据，并基于第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；间隔预设时长后，计算机设备可以获取受检者在第二时长内的第二pet数据，并基于第二时长内受检者的ct图像对利用第二pet数据得到的重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；然后，计算机设备根据第一时间活度曲线和第二时间活度曲线，估计受检者自第一时长开始时刻至第二时长截止时刻之间的完整时间活度曲线；最终，根据完整时间活度曲线进行动态pet参数成像。也就是说，本实施例中，计算机设备在第一时长内未进行ct扫描，仅在第二时长前进行了第一预设时长的ct扫描，且进行ct扫描之前，间隔了预设时长，大大减少了ct扫描时间及受检者保持静止的时间，从而大大减少了对受检者的辐射影响，并且，在第二时长内，受检者可以不用保持静止，进一步降低了临床操作难度；另外，计算机设备可以根据第一时间活度曲线和第二时间活度曲线，估计受检者自第一时长开始时刻至第二时长截止时刻之间的完整时间活度曲线，进而可以根据完整时间活度曲线进行动态pet参数成像，保证了动态pet参数成像的精准性。
第二方面，本发明实施例提供一种动态pet参数成像装置，包括：
获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；
间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述第二时长内所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；
根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线；
根据所述完整时间活度曲线进行动态pet参数成像。
第三方面，本发明实施例提供一种动态pet参数成像系统，包括pet扫描仪和计算机设备，其中计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时可用于执行一种图像重建方法，所述方法包括：
获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；
间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述第二时长内所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；
根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线；
根据所述完整时间活度曲线进行动态pet参数成像。
第四方面，本发明实施例提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；
间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述第二时长内所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；
根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线；
根据所述完整时间活度曲线进行动态pet参数成像。
本发明实施例提供的动态pet参数成像装置、系统以及计算机可读存储介质，能够使得计算机设备可以获取受检者在第一时长内的第一pet数据，并基于第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；间隔预设时长后，计算机设备可以获取受检者在第二时长内的第二pet数据，并基于第二时长内受检者的ct图像对利用第二pet数据得到的重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；然后，计算机设备根据第一时间活度曲线和第二时间活度曲线，估计受检者自第一时长开始时刻至第二时长截止时刻之间的完整时间活度曲线；最终，根据完整时间活度曲线进行动态pet参数成像。也就是说，本实施例中，计算机设备在第一时长内未进行ct扫描，仅在第二时长前进行了第一预设时长的ct扫描，且进行ct扫描之前，间隔了预设时长，大大减少了ct扫描时间及受检者保持静止的时间，从而大大减少了对受检者的辐射影响，并且，在第二时长内，受检者可以不用保持静止，进一步降低了临床操作难度；另外，计算机设备可以根据第一时间活度曲线和第二时间活度曲线，估计受检者自第一时长开始时刻至第二时长截止时刻之间的完整时间活度曲线，进而可以根据完整时间活度曲线进行动态pet参数成像，保证了动态pet参数成像的精准性。
附图说明
图1为一个实施例提供的一种pet系统的结构示意图；
图2为一个实施例提供的动态pet参数成像方法的流程示意图；
图3为另一个实施例提供的动态pet参数成像方法的流程示意图；
图4为另一个实施例提供的动态pet参数成像方法的流程示意图；
图5为一个实施例提供的动态pet参数成像装置的结构示意图；
图6为另一个实施例提供的动态pet参数成像装置的结构示意图；
图7为另一个实施例提供的动态pet参数成像装置的结构示意图；
图8为另一个实施例提供的动态pet参数成像装置的结构示意图；
图9为一个实施例提供的一种计算机设备的内部结构示意图。
附图标记说明：
10：pet扫描仪；11：扫描床；12：探测器；
13：机架；14：计算机设备。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，通过下述实施例并结合附图，对本发明实施例中的技术方案的进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
本发明提供的动态pet参数成像的方法，适用于各种图像处理系统，例如pet系统、正电子发射计算机断层显像-计算机断层扫描多模态系统pet-ct系统、正电子发射计算机断层显像-磁共振多模态混合系统pet-mr系统等。以pet系统为例，图1为一个实施例提供的pet系统的结构示意图，如图1所示，该pet系统可以包括pet扫描仪10和计算机设备14。pet扫描仪10可以包括扫描床11、探测器12和机架13。探测器12可以为环形探测器，还可以为矩形探测器。上述探测器12可以包括晶体、光电倍增管。pet扫描仪10和计算机设备14可以通过有线方式通信，也可以通过无线方式通信。其中，有线网络可以包括利用金属电缆、混合电缆、光纤、一个或多个接口等一种或多种组合的方式。无线网络可以包括利用蓝牙、区域局域网(lan)、广域局域网(wan)、近源场通信(nearfieldcommunication，nfc)等一种或多种组合的方式。
扫描床11用于支撑受检者(图1中未示出)，并可将受检者移动至扫描视野(fieldofview，fov)区域，探测器12中的晶体可将从受检者入射的γ光子转化为可视光，光电倍增管可以把晶体产生的可视光通过光电效应而转变成电信号，经逐级放大后，以电流的形式输出给电子线路系统(图1未示出)，电子线路系统将光电倍增管输出的微弱电信号(脉冲)进行必要的放大后，将采集到的信号转换成数字信号，送往计算机设备14进行数据处理，以使计算机设备14根据图像重建算法重建出示踪剂在生物体内的位置，从而得到示踪剂的代谢过程和分布图像，即pet图像。临床诊断中，需要对动态pet图像进行动力学定量分析，通过对动态pet图像分割得到感兴趣区域，提取该感兴趣区域的时间活度曲线，然后应用房室模型等动力学模型计算感兴趣区域的动力学参数理参数估计。
传统动态参数成像协议是比较严格的，需要病人从药物注射开始的整个采集过程中保持不动，临床难度较大。为解决上述问题，本实施例提出一种动态pet参数成像方法，适用于上述pet系统，尤其适用于pet扫描仪和ct扫描仪集成的多模态系统中。该方法可包括：获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，获取分帧重建图像；根据分帧重建图像得到第一时长内的第一时间活度曲线；间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，获取动态重建图像，受检者的ct图像为在第二时长扫描前采集重建获得；根据所述ct图像和所述动态重建图像得到第二时长内的第二时间活度曲线；根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线；以及根据所述完整时间活度曲线进行动态pet参数成像。可选地，根据所述ct图像和所述动态重建图像得到第二时长内的第二时间活度曲线包括：根据所述ct图像确定所述动态重建图像中与动脉区域对应的部分；根据所述动态重建图像中与动脉区域对应的部分的活度信息确定第二时长内的第二时间活度曲线。本发明对满足二室模型的非可逆性核素提出了一种新的扫描协议以及配套成像算法，在保证定量精度的同时显著缩短了病人需要保持静止的时间和扫描总时间，降低了操作难度，同时可有效降低患者辐射剂量。
需要说明的是，本发明实施例提供的动态pet参数成像的方法，其执行主体可以是动态pet参数成像装置，该装置可以通过软件、硬件或者软硬件结合的方式实现成为计算机设备的部分或者全部。可选的，该计算机设备可以为pc、便携式设备、服务器等具有数据处理功能的电子设备，本实施例对计算机设备的具体形式并不做限定。
图2为一个实施例提供的动态pet参数成像方法的流程示意图。本实施例涉及的是计算机设备根据受检者在第一时长内的第一pet数据和在第二时长内的第二pet数据，获取完整时间活度曲线(timeactivitycurve，tac，也可称之为动脉血时间活度曲线)，以进行动态pet参数成像的具体过程，该实施例中，第一时长和第二时长之间具有时间间隔，受检者在该时间间隔内可自由活动。如图2所示，该方法包括：
s201，获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线。
示例性的，pet扫描仪可以对受检者进行扫描，计算机设备可以获取pet扫描仪对受检者进行扫描的扫描结果。上述第一时长可以为操作者预设的时长，也可以根据所使用的造影剂的不同设置不同的第一时长，例如，当利用的造影剂为脱氧葡萄糖fdg时，可以设置第一时长为5分钟，在第一时长内，可以将pet扫描位定位在覆盖受检者胸腔主动脉的床位，然后对受检者进行弹丸注射fdg并进行pet扫描，获得第一pet数据。可选的，在第一时长截止时刻，受检者可以离开扫描床。
可选的，上述第一pet数据可以包括符合事件数据。可选的，符合事件数据指的是与响应线对应的符合事件数据，其可以为按归响应线的角度和径向距离编码形成的矩阵图。可选的，符合事件数据可包括真符合数据，而不包含散射符合事件数据和/或随机符合事件数据。可选的，上述第一pet数据的数据格式可以为正弦图格式(sinogram)，还可以为列表格式(listmode)。可选的，列表格式可以是将探测到的事件信息以流数据的形式依次记录，数据中可包括光子入射晶体条编号、光子能量和光子飞行时间信息。可选的，正弦图格式是将每条响应线上发生的符合事件数进行合并存储而形成的一种数据存储形式。可选的，上述符合事件数据可以用于重建受检者体内造影剂的空间分布以进行pet图像的重建和分帧图像的重建。
可选的，计算机设备可以根据有序子集最大期望值重建算法(orderedsubsetsexpectationmaximization，osem)进行预设迭代次数的迭代过程，获得pet图像和分帧图像，也可以利用滤波反投影法(filteredback-projection，fbp)进行图像重建，获得pet图像和分帧图像，本实施例对获得pet图像和分帧图像的重建算法不做限定。计算机设备获取重建的分帧图像后，可以基于重建的分帧图像，对时间活度曲线进行拟合，得到第一时长内的第一时间活度曲线。
s202，间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线。
本实施例中，上述预设时长可以由操作者根据造影剂的不同设置不同的预设时长，例如，利用fdg时，可以设置预设时长为40分钟，利用其它核素时，可以设置其它的时长。到达间隔预设时长后，受检者返回扫描床，先利用pet扫描仪对受检者进行ct扫描，到达第一预设时长后，再利用pet扫描仪对受检者进行第二预设时长的pet扫描。可选的，上述第二时长包括第一预设时长和第二预设时长，例如，可以设置第二时长为20分钟。可选的，在第二时长内，受检者可以进行一定范围的活动，而不用保持静止。可选的，上述pet扫描仪集成有ct扫描仪。分别利用ct扫描仪对受检者在第二时长前进行ct扫描，利用pet扫描仪对受检者进行第二时长的pet扫描后，计算机设备可以获取受检者的ct数据和第二pet数据，重建ct数据可获取ct图像。可以理解地，本发明中的第一时长对应第一阶段的扫描时间，第二时长对应第二阶段的扫描时间，且第一时长早于第二时长。
开始进行第二时长的扫描时，受检者体内的造影剂的已达到稳态，其变化速率已经较为缓慢，可以设置较长的数据分帧时长，例如可以设置数据分帧时长为1分钟至5分钟。计算机设备可以根据分帧后的每一帧第二pet数据进行分帧图像重建，并在第二pet数据重建过程中利用ct图像对第二pet数据进行衰减校正，得到多帧图像构成的动态重建图像。
可选的，上述基于所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线的实现过程，可以包括：
s2021，根据所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，获取动态重建图像。
示例性的，计算机设备可以根据osem重建算法进行预设迭代次数的迭代过程，获得动态重建图像，也可以利用fbp重建算法进行图像重建，获得动态重建图像，本实施例对获得动态重建图像的重建算法不做限定。
s2022，根据所述ct图像确定所述动态重建图像中与动脉区域对应的部分。
示例性的，当计算机设备获取受检者的动态重建图像后，计算机设备可以基于上述ct图像进行主动脉分割，可以在动态重建图像的各个分帧图像上进行主动脉分割，还可以将ct图像和动态重建图像的各个分帧图像进行图像配准获取主动脉的区域。可选的，计算机设备可以利用阿特拉斯(atlas)模板基于ct图像分割出主动脉，还可以根据人工智能的方法分割出主动脉，本实施例对进行主动脉分割的方式不做限定。
s2023，根据所述动态重建图像中与动脉区域对应的部分的活度信息得到第二时长内的第二时间活度曲线。
示例性的，计算机设备可以根据动态重建图像中，各个分帧图像上分割出的主动脉区域的造影剂的活度，对时间活度曲线作拟合，得到第二时长内的第二时间活度曲线。可选的，计算机设备可以利用线性拟合的方式得到第二时长内的第二时间活度曲线，也可以根据指数拟合的方式得到第二时长内的第二时间活度曲线，还可以根样条曲线(splinecurves)拟合的方式得到第二时长内的第二时间活度曲线，本实施例对获得第二时长内的第二时间活度曲线的方式不做限定。
s203，根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线。
示例性的，计算机设备可以根据步骤s201获得的第一时间活度曲线和步骤s202获得的第二时间活度曲线，以及利用pet扫描仪对部分受检者进行预设扫描时长的扫描后，获得的时间活度曲线中的至少一个时间活度曲线，估计受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线。
可选的，计算机设备可以根据公式
上述s(t)可以为根据采样的多条时间活度曲线进行归一化并求平均获得的时间活度曲线，例如，可以选取10个受检者，获取该10个受检者的同等时长的10条时间活度曲线，对获得的10条时间活度曲线进行归一化，获得10条归一化的时间活度曲线，最后对10条归一化的时间活度曲线取平均，可以得到s(t)。
s204，根据所述完整时间活度曲线进行动态pet参数成像。
示例性的，计算机设备可以根据估计的完整时间活度曲线，以及动态重建图像的各个分帧图像，利用pet动态参数估计模型，如patlak参数模型或logan模型，进行动态pet参数成像。
可选的，计算机设备可以根据公式
本实施例提供的动态pet参数成像方法，计算机设备可以获取受检者在第一时长内的第一pet数据，并基于第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；间隔预设时长后，计算机设备可以获取受检者在第二时长内的第二pet数据，并基于第二时长内受检者的ct图像对利用第二pet数据得到的重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；然后，计算机设备根据第一时间活度曲线和第二时间活度曲线，估计受检者自第一时长开始时刻至第二时长截止时刻之间的完整时间活度曲线；最终，根据完整时间活度曲线进行动态pet参数成像。也就是说，本实施例中，计算机设备在第一时长内未进行ct扫描，仅在第二时长前进行了第一预设时长的ct扫描，且进行ct扫描之前，间隔了预设时长，大大减少了ct扫描时间及受检者保持静止的时间，从而大大减少了对受检者的辐射影响，并且，在第二时长内，受检者可以不用保持静止，进一步降低了临床操作难度；另外，计算机设备可以根据第一时间活度曲线和第二时间活度曲线，估计受检者自第一时长开始时刻至第二时长截止时刻之间的完整时间活度曲线，进而可以根据完整时间活度曲线进行动态pet参数成像，保证了动态pet参数成像的精准性。
图3为另一个实施例提供的动态pet参数成像方法的流程示意图。本实施例涉及的是基于所述第一pet数据，得到第一时长内的第一时间活度曲线的一种可选的过程。在上述实施例的基础上，可选的，上述s201中的“基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线”可以包括：
s301，采用mlaa(活度与衰减信息联合估计)算法所述第一pet数据进行pet图像重建，并在所述pet图像重建过程中确定所述第一时长内的目标衰减图像。
可选的，计算机设备可以采用公式
上述
可选的，上述
计算机设备可以根据所有体素任意设定(通常可查找字典库获得)衰减值得到初始衰减图像μ0，并将初始衰减图像作为非零迭代初始值，利用公式
得到
然后根据公式
获得
可选的，上述迭代过程达到停止迭代条件时，停止迭代，最终得到第一时长内的目标衰减图像。可选的，停止迭代条件可以为预设的迭代次数，也可以为目标衰减图像的分辨率达到预设阈值时，停止迭代。停止迭代后，计算机设备可以得到目标衰减图像。
s302，对所述第一pet数据进行重建以获取分帧重建图像，且在所述分帧重建图像的重建过程中根据所述第一时长内的目标衰减图像对所述第一pet数据进行衰减校正。采用所述目标衰减图像对所述第一pet数据进行衰减校正，得到所述第一时间活度曲线。
示例性的，计算机设备可以将第一pet数据进行数据分帧，基于分帧数据进行分帧图像重建。由于第一时长内，动脉血中的造影剂的活变化速率比较快，可以设置较短的数据分帧时长，例如可以设置数据分帧时长为10秒至30秒。可选的，计算机设备可以根据分帧后的每一帧第二pet数据进行分帧图像重建。可选的，可以利用确定的目标衰减图像对第二pet数据的分帧重建图像进行衰减校正。可选的，计算机设备可以根据osem重建算法进行预设迭代次数的迭代过程，获得分帧重建图像，也可以利用fbp重建算法进行图像重建，获得分帧重建图像，本实施例对获得第一分帧图像集合的重建算法不做限定。
s303，根据所述分帧重建图像得到所述第一时间活度曲线。
示例性的，计算机设备可以根据分帧重建图像得到第一时间活度曲线。
可选的，计算机设备可以根据分帧重建图像得到第一时间活度曲线的实现过程可以包括：
s3031，根据所述第一分帧图像集合中多个分帧图像的活度信息，确定所述第一分帧图像集合中各个分帧图像中的动脉区域。
示例性的，上述分帧重建图像包括第一分帧图像集合。上述活度信息为第一分帧图像集合中各个分帧图像对应的造影剂的活度信息。可选地，第一分帧图像集合包括第一时段的分帧图像、第二时段的分帧图像，根据所述第一分帧图像集合中多个分帧图像的活度信息，确定所述第一分帧图像集合中各个分帧图像中的动脉区域包括：从所述第一分帧图像集合中选取第一时段的分帧图像；根据所述第一时段的分帧图像的活度信息，分别确定所述第一时段的分帧图像、第二时段的分帧图像的动脉区域。需要说明的是，第一时段在第一时长中处于第二时段之前，即第一时段早于第二时段。可选地，第一时段可以是第一时长的前半段，第二时段可以是第一时长的后半段。
s3032，根据所述各个分帧图像的动脉区域中的造影剂的活度变化，确定所述第一时间活度曲线。
计算机设备获取第一分帧图像集合中各个分帧图像中的动脉区域后，可以根据分割出的动脉区域内的造影剂的活度均值作为采样点，对动脉血时间活度曲线进行拟合，得到第一时长内的第一动脉血时间活度曲线。可选的，计算机设备可以利用线性拟合的方式得到第一时长内的第一动脉血时间活度曲线，也可以根据指数拟合的方式得到第一时长内的第一动脉血时间活度曲线，还可以根样条曲线spline拟合的方式得到第一时长内的第一动脉血时间活度曲线，本实施例对获得第一时长内的第一动脉血时间活度曲线的方式不做限定。
本实施例提供的动态pet参数成像方法，计算机设备可以采用mlaa算法对所述第一pet数据进行pet图像重建，并在所述pet图像重建过程中确定所述第一时长内的目标衰减图像；对所述第一pet数据进行重建以获取分帧重建图像，且在所述分帧重建图像的重建过程中根据所述第一时长内目标衰减图像对所述第一pet数据衰减校正；根据所述分帧重建图像得到所述第一时间活度曲线。也就是说，本实施例并未通过ct扫描获得ct图像作为目标衰减图像，这样在获得目标衰减图像的同时，大大减少了ct扫描时间，进而大大减少了对受检者的辐射影响。
为了便于本领域技术人员的理解，以下对本发明提供的动态pet参数成像方法进行详细介绍，示例性的，如图4所示，该方法可以包括：
s501，计算机设备获取受检者在第一时长内的第一pet数据。
s502，计算机设备根据公式利用mlaa算法，根据公式
s503，计算机设备根据所述目标衰减图像对所述第一pet数据进行分帧图像重建，得到第一分帧图像集合。
s504，计算机设备根据所述第一分帧图像集合中一个或多个分帧图像的活度信息，确定所述第一分帧图像集合中各个分帧图像中的动脉区域。
s505，计算机设备根据所述各个分帧图像中的动脉区域中的造影剂的活度变化，确定所述第一时间活度曲线。
s506，间隔预设时长后，计算机设备获取受检者在第二时长内的第二pet数据，并基于所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线，ct图像对应的ct数据采集在第二时长前进行。
s507，计算机设备根据公式
s508，根据公式
需要说明的是，针对上述s501-s508中的描述可以参见上述实施例中相关的描述，且其效果类似，本实施例在此不再赘述。
应该理解的是，虽然图2-4的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2-4中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
图5为一个实施例提供的动态pet参数成像装置的结构示意图。如图5所示，该动态pet参数成像装置可以包括：第一确定模块10、第二确定模块11、估计模块12和成像模块13。
示例性的，第一确定模块10，用于获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；
第二确定模块11，用于间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述受检者的ct图像对所述第二pet数据进行衰减校正，得到第二时长内的第二时间活度曲线；
估计模块12，用于根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线；
成像模块13，用于根据所述完整时间活度曲线进行动态pet参数成像。
本实施例提供的pet参数成像装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
图6为另一个实施例提供的pet参数成像装置的结构示意图。在如图5所示实施例的基础上，如图6所示，上述第一确定模块10可以包括：第一确定单元101、校正单元102和第一获取单元103。
第一确定单元101，用于采用mlaa算法对所述第一pet数据进行pet图像重建，并在所述pet图像重建过程中确定所述第一时长内的目标衰减图像。
校正单元102，用于对所述第一pet数据进行重建以获取分帧重建图像，且在所述分帧重建图像的重建过程中根据所述第一时长内的目标衰减图像对所述第一pet数据进行衰减校正。
第一获取单元103，根据所述分帧重建图像得到所述第一时间活度曲线。
本实施例提供的pet参数成像装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
图7为另一个实施例提供的pet参数成像装置的结构示意图。在上述实施例的基础上，可选的，上述第一获取单元103可以包括：第一确定子单元1031和第二确定子单元1032。
第一确定子单元1031，用于根据所述第一分帧图像集合中多个分帧图像的活度信息，确定所述第一分帧图像集合中各个分帧图像中的动脉区域。
可选的，上述第一确定子单元1031具体用于从所述第一分帧图像集合中选取第一时段的分帧图像，并根据所述第一时段的分帧图像的活度信息，分别确定所述第一时段的分帧图像、第二时段的分帧图像的动脉区域。
第二确定子单元1032，用于根据所述各个分帧图像的动脉区域中的造影剂的活度变化，确定所述第一时间活度曲线。
本实施例提供的pet参数成像装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
图8为另一个实施例提供的pet参数成像装置的结构示意图。在上述实施例的基础上，可选的，上述第二确定模块11包括第二获取单元111、第二确定单元112和第三获取单元113。
第二获取单元111，用于根据所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，获取动态重建图像；
第二确定单元112，用于根据所述ct图像确定所述动态重建图像中与动脉区域对应的部分；
第三获取单元113，根据所述动态重建图像中与动脉区域对应的部分的活度信息得到第二时长内的第二时间活度曲线。
本实施例提供的pet参数成像装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
关于动态pet参数成像装置的具体限定可以参见上文中对于动态pet参数成像方法的限定，在此不再赘述。上述动态pet参数成像装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种动态pet参数成像系统，如图1所示，如图1所示，该动态pet参数成像系统包括pet扫描仪10和计算机设备14，该计算机设备14可以是终端，其内部结构图可以如图9所示。该计算机设备14包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备14的处理器用于提供计算和控制能力。该计算机设备14的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备14的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现如下步骤：
获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；
间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；
根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线；
根据所述完整时间活度曲线进行动态pet参数成像。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取受检者在第一时长内的第一pet数据，并基于所述第一pet数据进行分帧图像重建，得到第一时长内的第一时间活度曲线；
间隔预设时长后，获取受检者在第二时长内的第二pet数据，并基于所述受检者的ct图像对所述第二pet数据的分帧重建图像进行衰减校正，得到第二时长内的第二时间活度曲线；
根据所述第一时间活度曲线和所述第二时间活度曲线，估计所述受检者自所述第一时长开始时刻至所述第二时长截止时刻之间的完整时间活度曲线；
根据所述完整时间活度曲线进行动态pet参数成像。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(rom)、可编程rom(prom)、电可编程rom(eprom)、电可擦除可编程rom(eeprom)或闪存。易失性存储器可包括随机存取存储器(ram)或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram(sram)、动态ram(dram)、同步dram(sdram)、双数据率sdram(ddrsdram)、增强型sdram(esdram)、同步链路(synchlink)dram(sldram)、存储器总线(rambus)直接ram(rdram)、直接存储器总线动态ram(drdram)、以及存储器总线动态ram(rdram)等。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。