CN114862980A - 散射校正方法、pet成像方法、装置、设备及存储介质 - Google Patents
散射校正方法、pet成像方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN114862980A CN114862980A CN202210487528.0A CN202210487528A CN114862980A CN 114862980 A CN114862980 A CN 114862980A CN 202210487528 A CN202210487528 A CN 202210487528A CN 114862980 A CN114862980 A CN 114862980A
- Authority
- CN
- China
- Prior art keywords
- image
- tof
- pet
- scatter
- training
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
Abstract
本发明公开一种散射校正方法、PET成像方法、装置、设备及存储介质,散射校正方法包括:获取目标对象在扫描期间的TOF‑PET数据和衰减图像;基于所述TOF‑PET数据获取TOF信息的反投影图像;基于所述反投影图像和所述衰减图像,获取TOF‑PET成像的散射估计;基于所述TOF‑PET成像的散射估计对所述TOF‑PET数据进行散射校正。本发明解决了现有技术中TOF‑PET成像时散射校正耗时长的技术问题。
Description
技术领域
本发明涉及医疗技术领域,具体涉及一种散射校正方法、PET成像方法、装置、设备及存储介质。
背景技术
正电子发射计算机断层扫描(PET扫描)是一种检测人体或动物体器官代谢特征的无创医学成像技术,具有灵敏度高、准确性好、定位准确的特点。
飞行时间技术正电子发射断层成像(TOF-PET)扫描仪是核医学成像中的一种先进的功能成像工具,其应用前景已经受到核医学成像研究者和设备制造厂商的高度重视。在最新的TOF-PET扫描仪中,虽然相对于传统2D和扩展2D数据采集模式,对所有倾斜响应线的符合计数采集提高了系统的灵敏度,但同时也引入了大量的散射符合计数,从而降低了系统的分辨率和成像的定量准确性,因此,为了提高TOF-PET成像的准确性,需要先进行散射校正。
在传统的散射重建中,通常采用未经过散射校正的重建图像作为示踪剂分布图像,结合衰减图像,估计初步的散射投影域分布,并以此为初步散射估计,重建出新的示踪剂分布图像;如此反复迭代数次直至散射分布收敛。迭代过程需要反复的重建图像和进行散射估计,因而非常耗时。
发明内容
本发明的目的在于克服上述技术不足,提供一种散射校正方法、PET成像方法、装置、设备及存储介质,解决现有技术中TOF-PET成像时散射校正耗时长的技术问题。
为达到上述技术目的,本发明采取了以下技术方案:
第一方面,本发明提供一种散射校正方法,包括如下步骤:
获取目标对象在扫描期间的TOF-PET数据和衰减图像;
基于所述TOF-PET数据获取TOF信息的反投影图像;
基于所述反投影图像和所述衰减图像,获取TOF-PET成像的散射估计;
基于所述TOF-PET成像的散射估计对所述TOF-PET数据进行散射校正。
在其中一些实施例中,所述基于所述反投影图像和所述衰减图像,获取TOF-PET成像的散射估计,包括:
将所述衰减图像和所述反投影图像输入至训练完备的第一训练模型中,以得到示踪剂分布图像;
基于所述示踪剂分布图像和所述衰减图像,获取所述TOF-PET成像的散射估计。
在其中一些实施例中,所述散射校正方法还包括:
获取第一训练集并构建第一初始训练模型,其中,所述第一训练集包括若干个第一目标参数以及若干组与所述第一目标参数对应的第一输入参数,所述第一输入参数为反投影图像以及与其对应的衰减图像,所述第一目标参数为示踪剂分布图像;
采用所述第一训练集对所述第一初始训练模型进行训练,以得到训练完备的第一训练模型。
在其中一些实施例中,所述第一训练模型为机器学习模型或者深度学习模型。
在其中一些实施例中,所述获取目标对象在扫描期间的TOF-PET数据和衰减图像之后,所述基于所述TOF-PET数据获取TOF信息的反投影图像之前还包括:
对所述TOF-PET数据进行预校正处理。
第二方面,本发明还提供一种PET成像方法,包括如下步骤:
获取目标对象在扫描期间的TOF-PET数据和静态图像,其中,所述静态图像通过如上所述的散射校正方法进行散射校正后获取;
基于所述TOF-PET数据获取TOF信息的反投影图像;
将所述静态图像和所述反投影图像输入训练完备的第二训练模型中,以得到PET重建图像。
在其中一些实施例中,所述PET成像方法还包括:
获取第二训练集并构建第二初始训练模型,其中,所述第二训练集包括若干个第二目标参数以及若干组与第二目标参数对应的第二输入参数,所述第二输入参数为静态图像以及反投影图像,所述第二目标参数为PET重建图像;
采用第二训练集对所述第二初始训练模型进行训练,以得到训练完备的第二训练模型。
第三方面,本发明还提供一种散射校正装置,包括:
初始数据获取模块,用于获取目标对象在扫描期间的TOF-PET数据和衰减图像;
反投影图像获取模块,用于基于所述TOF-PET数据获取TOF信息的反投影图像;
散射估计获取模块,用于基于所述反投影图像和所述衰减图像,获取TOF-PET成像的散射估计;
散射校正模块,用于基于所述TOF-PET成像的散射估计对所述TOF-PET数据进行散射校正。
第四方面,本发明还提供一种电子设备,包括:处理器和存储器;
所述存储器上存储有可被所述处理器执行的计算机可读程序;
所述处理器执行所述计算机可读程序时实现如上所述的散射校正方法或者PET成像方法中的步骤。
第五方面,本发明还提供一种计算机可读存储介质,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如上所述的散射校正方法或者PET成像方法中的步骤。
与现有技术相比,本发明提供的散射校正方法、PET成像方法、装置、设备及存储介质,首先获取目标对象在扫描期间的TOF-PET数据和衰减图像,然后对TOF-PET数据进行反投影处理得到TOF信息的反投影图像,之后再根据该反投影图像与衰减图像获取散射估计,最后根据散射估计对TOF-PET数据进行散射校正,从而能够有效的减少用于散射估计的迭代次数,在加快散射校正速度的同时,保证散射估计的准确性。在PET成像时,采用经过散射校正的静态图像和反投影图像结合进行PET重建成像,经过散射校正的静态图像给动态的TOF-PET数据提供了散射校正核,因此,在进行PET成像时,可以加快图像重建的速度,而且通过对第二训练模型的训练,在进行重建时,可直接将静态图像和反投影图像输入至第二训练模型中后,得到PET重建图像,进而可以进一步加快图像重建的速度。
附图说明
图1是本发明提供的散射校正方法的一实施例的流程图;
图2是本发明提供的PET成像方法的一实施例的流程图;
图3是本发明提供的散射校正装置的一实施例的示意图;
图4是本发明提供的PET成像装置的一实施例的示意图;
图5是本发明计算机程序的一实施例的运行环境示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明所涉及的散射校正方法、PET成像方法、装置、设备或计算机可读存储介质可用于正电子发射计算机断层显像-X射线计算机断层造影多模态混合系统(PET-CT系统)、正电子发射型断层显像-磁共振成像系统(PET-MRI系统)等医学成像系统中。本发明所涉及的方法、装置、设备或者计算机可读存储介质既可以与上述系统集成在一起,也可以是相对独立的。
本实施例中提供了一种散射校正方法,可由PET-CT设备、PET-MRI设备等医学成像设备执行,具体可由该设备的一个或者多个处理器执行。图1是本发明实施例提供的散射校正方法的流程图,请参阅图1,散射校正方法包括如下步骤:
S110、获取目标对象在扫描期间的TOF-PET数据和衰减图像;
S120、基于TOF-PET数据获取TOF信息的反投影图像;
S130、基于反投影图像和衰减图像,获取TOF-PET成像的散射估计;
S140、基于TOF-PET成像的散射估计对TOF-PET数据进行散射校正。
本实施例中,首先获取目标对象在扫描期间的TOF-PET数据和衰减图像,然后对TOF-PET数据进行反投影处理得到TOF信息的反投影图像,之后再根据该反投影图像与衰减图像获取散射估计,最后根据散射估计对TOF-PET数据进行散射校正,从而能够有效的减少用于散射估计的迭代次数,在加快散射校正速度的同时,保证散射估计的准确性。
在一些实施例中,步骤S110中,TOF-PET数据可以包括符合事件数据和TOF信息。本实施例中,符合事件数据指的是与响应线对应的符合事件数据,其可以为按响应线的角度和径向距离编码形成的矩阵图。符合事件数据也可称为PET投影数据。符合事件数据包括真符合数据、散射符合事件数据及随机符合事件数据。符合事件数据可以用于重建受检者体内造影剂的空间分布并获取扫描图像。优选的,TOF-PET数据的数据格式可以为正弦图格式(sinogram),还可以为列表格式(list mode)。列表格式可以是将探测到的事件信息以数据流的形式依次记录,数据中可包括光子入射晶体条编号、光子能量和光子飞行时间信息。正弦图格式是将每条响应线上发生的符合事件数进行合并存储而形成的一种数据存储形式。上述TOF信息可以通过PET探测器检测得到,并且其可以用概率密度函数表示。
在一些实施例中,TOF-PET数据可以通过PET扫描仪探测到的数据经计算机设备进行处理后获得,具体获得过程为:操作者预先向受检者体内注入标记了核素(如F18、C11、O15、Ga68或Ru82等)的示踪剂,示踪剂在受检者体内通过扩散作用而进入各个组织或者血管,然后通过PET探测器在体外探测示踪剂的正电子与受检者体内的负电子发生湮灭效应而产生的辐射信号(该信号为能量相等、方向相反的一对γ光子),探测器将辐射信号进行光电转换形成电信号,然后探测器将探测到的信号传输给电子线路系统,电子线路系统可以把探测器采集到的信号转化成数字信号,该数字信号即为TOF-PET数据。
在一些实施例中,衰减图像可以来源于MR扫描、CT扫描等过程,当其来源于CT扫描过程时,具体获得过程为:在PET扫描之前,先进行CT扫描,在CT扫描过程中获取CT扫描图像,然后基于CT扫描图像显示组织衰减信息,即得到了衰减图像,该衰减图像可用于辅助PET图像进行散射估计。当其来源于MR扫描过程时,具体获取过程为:经过UTE序列转换成为衰减图像。
步骤S120中,相较于传统的TOF-PET数据都是以表模式或者正弦图格式保存,本实施例中,将TOF-PET数据转化为反投影图像模式,反投影图像可以应用于ML-EM(最大期望重建)算法中进行重建,称为DIRECT(直接)重建。在一些实施例中,反投影图像可以为TOFhisto-image或者TOF histo-projection,其中,TOF histo-image是根据扫描数据中同一湮灭事件的两个光子到达两端探测器的时间差确定多个湮灭点的物理位置,并将其对应到图像域像素点上进行累加得到的图像,TOF histo-projection是在生成反投影图像时,根据符合事件的方位角,筛选出特定范围方位角的符合事件,进行直接反投影操作所生成的图像,所有角度的反投影图像叠加等价所有事件的反投影图像。其中,TOF信息的反投影图像可以采用直接投点法或者基于TOF的滤波反投影法获取,滤波反投影法是在傅立叶变换理论基础之上的一种空域处理技术。它的特点是在反投影前将每一个采集投影角度下的投影进行卷积处理,从而改善点扩散函数引起的形状伪影,重建的图像质量较好。其具体过程为:先把投影由线阵探测器上获得的投影数据进行一次一维傅立叶变换,再与滤波器函数进行卷积运算,得到各个方向卷积滤波后的投影数据;然后把它们沿各个方向进行反投影;再经过适当处理后得到被扫描物体的断层图像。本发明实施例中,优选为直接投点法来进行TOF信息的反投影图像的获取,该方法速度快,通量高,可以做到实时成像。
在一些实施例中,为了增加TOF信息的反投影图像的精确度,步骤S110和步骤S120之间还包括:
对TOF-PET数据进行预校正处理。
本实施例中,预校正处理的目的是为了减去散射数据,从而增加散射校正的精准度和速度。在具体实现过程中,可以在获取TOF信息的反投影图像的过程中,选择性加入图像域的校正数据,比如基于图像域的随机校正和基于能量信息的散射初步校正数据,进一步改善TOF信息的反投影图像的精确性。其中,基于图像域的随机校正指把随机数据投影到图像域,基于能量信息的散射初步校正指把不同能量谱所对应的散射投影数据投影至图像域。
步骤S130是用于进行TOF-PET成像的散射估计,在同一扫描中,无明显运动影响时,衰减图像保持不变。因此可以提取散射校正核,用于动态多帧图像的重建。在一些实施例中,可选用单散射模拟法(Single Scatter Simulation,SSS),蒙特卡洛模拟法(Monte-Carlo Simulation,MCS),双散射模拟法(Double Scatter Simulation,DSS)或者剩余估计法(ResidualEstimation Approach)等估计方法来获取散射估计,本发明对此不做限定,获取的散射估计可供TOF图像或者non-TOF的重建使用。其中,产生的散射估计可以为sinogram(正弦图格式)数据、list-mode(列表格式)数据等相关投影域数据格式。
在一些实施例中,为了提高散射估计的准确性,并快速的进行散射估计,步骤S130包括:
将衰减图像和反投影图像输入至训练完备的第一训练模型中,以得到示踪剂分布图像;
基于示踪剂分布图像和衰减图像,获取TOF-PET成像的散射估计。
本实施例中,通过训练模型来进行示踪剂分布图像的生成,通过第一训练模型得到的示踪剂分布图像用于进行散射校正,而且已经经过初步散射校正,因此在后续的散射校正步骤中,可以有效减少用于散射估计的迭代次数,保证散射估计的准确性。为了保证第一训练模型能够较好的根据TOF信息的反投影图像和衰减图像获取示踪剂分布图像,需要先对第一训练模型进行训练。
在本发明的一些实施例中,建立第一训练模型的具体过程为:首先通过算法建立初始模型,然后基于历史示踪剂分布图像、历史反投影图像以及与其对应的衰减图像对初始模型进行训练,训练完成后,即可得到第一训练模型,根据第一训练模型即可通过输入的反投影图像及其对应的衰减图像,输出示踪剂分布图像。
优选的,散射校正方法还包括:
获取第一训练集并构建第一初始训练模型,其中,第一训练集包括若干个第一目标参数以及若干组与第一目标参数对应的第一输入参数,第一输入参数为反投影图像以及与其对应的衰减图像,第一目标参数为示踪剂分布图像;
采用第一训练集对第一初始训练模型进行训练,以得到训练完备的第一训练模型。
本实施例中,用于训练的第一训练集中的每个个体都包括两个第一输入参数和一个第一目标参数,两个第一输入参数即为反投影图像以及与其对应的衰减图像,其中,反投影图像和衰减图像可通过步骤S110至步骤S120的方法获取,第一目标参数为示踪剂分布图像,该图像可以采用有序子集最大期望值重建算法(Ordered Subsets ExpectationMaximization,OSEM)获取,有序子集最大期望值重建算法是一种PET图像重建的方法,其具体实现过程为:将全部投影数据分成多个子集,每使用一个子集的数据,全部像素被更新一次,成为一次子迭代,所有子集轮流使用一遍,完成一次迭代。每一次迭代会使图像收敛一些,多次迭代后得到重建图像,该重建图像更接近真实值。因此,本发明实施例中,将反投影图像以及与其对应的衰减图像作为第一初始训练模型的输入,将经过示踪剂分布图像作为第一初始训练模型的输出来进行模型训练,相当于模型自动学习了投影域的散射分布,进行了初步散射估计,故可以快速进行散射估计,训练完成后,可直接通过反投影图像以及与其对应的衰减图像得到真实的示踪剂分布图像,相较于现有技术,有效的减少了散射估计所需的迭代次数。而且,衰减图像提供的是物理上会发生散射的物质分布,能够增加训练数据的信息量。此外,由于传统的散射校正是基于衰减图和示踪剂分布图像的,故而本发明实施例中加入衰减图像作为输入,给训练网络增加了额外的信息量,易于网络的学习,最终加强了训练网络的鲁棒性。
在一些实施例中,第一训练模型为机器学习模型或者深度学习模型。第一初始训练模型也对应为机器学习模型或者深度学习模型。其中,机器学习模型可以包括但不限于线性回归模型、岭回归模型、支持向量回归模型、支持向量机、决策树、全连接神经网络、循环神经网络等。深度学习模型可以包括但不限于深度信念网络(Deep belief Networks,DBN)、卷积神经网络(Convolution Neural Networks,CNN)、递归神经网络(RecurrentNeural Networks,RNN)、循环神经网(Recurrent Neural Networks,RNN)、生成式对抗网络(GAN,GenerativeAdversarial Networks)等等。在本发明的具体实施例中,第一训练模型的网络结构可为2DU-net、3DU-net、U-net++、U-net3+、V-net中的任意一种。当然,在其它的实施例中,还可采用其它的模型来实现本发明功能,本发明对此不做限定。
应当理解的是,深度学习模型可以是能够从已有数据(历史医疗图像)中学习到一定的知识和能力用于处理新数据,并可以被设计用于执行各种任务,在本发明的实施例中用于进行示踪剂分布图像的生成。
在一些实施例中,第一训练模型可以为头部训练模型,对应的衰减图像为头部衰减图像,第一训练模型也可以为体部训练模型,对应的衰减图像为体部衰减图像,其中,体部训练模型的参数与头部训练模型的参数可能有所区别或者能够共用,其具体实现方式可根据实际需求而定,本发明实施例对此不做限定。
当通过第一训练模型得到示踪剂分布图像后,可根据衰减图像和示踪剂分布图像来获取TOF-PET成像的散射估计,在同一扫描中,无明显运动影响时,衰减图像保持不变。因此可以提取散射校正核,用于动态多帧图像的重建。在一些实施例中,可选用单散射模拟法(Single Scatter Simulation,SSS),蒙特卡洛模拟法(Monte-Carlo Simulation,MCS)或者双散射模拟法(Double Scatter Simulation,DSS)等估计方法来获取散射估计,本发明对此不做限定,获取的散射估计可供TOF图像或者non-TOF的重建使用。其中,产生的散射估计可以为sinogram(正弦图格式)数据、list-mode(列表格式)数据等相关投影域数据格式。
步骤S140是为了对TOF-PET数据进行散射校正,具体通过散射估计数据可以来校正TOF-PET数据。由于TOF-PET散射估计接近真实值,可有效提高重建获得的TOF-PET图像精度。
本发明实施例中,采用TOF信息直接反投影得到反投影图像,该方法速度快,通量高,可以做到实时成像。通过第一训练模型,构建反投影图像到真实示踪剂分布图像的转换模型,从而根据每次生成的反投影图像快速准确的估计用于散射的示踪剂分布图像。该方法能够有效减少用于散射估计的迭代次数,保证散射估计的准确性。
基于上述散射校正方法,本发明实施例还提供一种PET成像方法,请参阅图2,PET成像方法包括如下步骤:
S210、获取目标对象在扫描期间的TOF-PET数据和静态图像;
S220、基于TOF-PET数据获取TOF信息的反投影图像;
S230、将静态图像和反投影图像输入训练完备的第二训练模型中,以得到PET重建图像。
在一些实施例中,静态图像是经过如上述各实施例所述的散射校正方法获取的静态图像,由于上文已对散射校正方法进行详细描述,在此不再赘述。
在一些实施例中,静态图像还可通过传统的方法获取,具体实施时,采用不包含散射校正的重建图像作为初始的示踪剂分布图像,结合衰减图像,估计初步的散射投影域分布,并以此为初步散射估计进行下一次包含散射估计的图像重建,重建出新的示踪剂分布图像;如此反复迭代数次直至散射分布收敛后,重建出静态图像。
本实施例提供的PET重建方法,采用静态图像和TOF-PET数据结合进行PET重建成像,静态图像已经过散射校正,故可以给TOF-PET反投影图像提供散射校正核,因此,在进行PET成像时,可以加快图像重建的速度,而且通过对第二训练模型的训练,在进行重建时,可直接将静态图像和反投影图像输入至第二训练模型中后,得到PET重建图像,进而可以进一步加快图像重建的速度。
在一些实施例中,步骤S210中,静态图像是扫描期间获取的,静态图像指常规临床用于诊断而重建的图像,扫描时间较长。举例来说,假如是动态扫描的持续时间为30min,可以看成是由多张静态图像按照时序连在一起形成的,每张静态图像扫描时间为2min或更长的时间,我们可选择一段2min或更长时间的静态扫描图像作为模型输入。其中,静态图像一般都是指做了散射校正的图像,可直接用于PET图像的重建。
步骤S220是通过获取原始TOF-PET数据后,转换为反投影图像。在一些实施例中,为了保证减少原始数据中的散射数据,增加散射校正的精准度和速度,在获取TOF-PET数据之前,还对TOF-PET数据进行预校正处理,预校正处理的过程包括但不限于随机校正、归一化校正和散射初步校正,其中,随机校正是指把随机数据投影到图像域。归一化校正是利用PET设备,对一个均匀源进行采集,采集到足够量的数据后,对数据进行统计分析,形成探测器单元的效率因子,然后在TOF-PET数据采集时,探测器采集的数据乘以探测器单元的效率因子,实现归一化校正。散射初步校正指把不同能量谱所对应的散射投影数据投影至图像域。
步骤S230是用于通过训练模型来进行PET图像重建,重建速度快。为了保证第二训练模型能够获取更加精准的重建图像,需要先对第二训练模型进行训练,在本发明的一些实施例中,建立第二训练模型的具体过程为:首先通过算法建立初始模型,然后基于传统PET重建图像、历史反投影图像以及静态图像对初始模型进行训练,训练完成后,即可得到第二训练模型,根据第二训练模型即可通过输入的反投影图像及静态图像,输出PET重建图像。其中,用于训练的静态图像可采用传统方法或者如上述各实施例所述的散射校正方法获取,作为模型输入。
优选的,PET成像方法还包括:
获取第二训练集并构建第二初始训练模型,其中,第二训练集包括若干个第二目标参数以及若干组与第二目标参数对应的第二输入参数,第二输入参数为静态图像以及反投影图像,第二目标参数为PET重建图像;
采用第二训练集对第二初始训练模型进行训练,以得到训练完备的第二训练模型。
本实施例中,用于训练的第二训练集中的每个个体都包括两个第二输入参数和一个第二目标参数,两个第二输入参数即为静态重建图像以及反投影图像,其中,静态图像以及反投影图像可通过步骤S210至步骤S220的方法获取,第二目标参数为PET重建图像,该PET重建图像可以采用传统的重建方法获取。因此,本发明实施例中,将静态图像以及反投影图像作为第二初始训练模型的输入,将传统的PET重建图像作为第二初始训练模型的输出来进行模型训练,相当于模型自动学习了PET重建的过程,可直接通过静态图像以及反投影图像得到PET重建图像,相较于现有技术,加快了重建的速度。
在一些实施例中,第二训练模型为机器学习模型或者深度学习模型。第二初始训练模型也对应为机器学习模型或者深度学习模型。其中,机器学习模型可以包括但不限于线性回归模型、岭回归模型、支持向量回归模型、支持向量机、决策树、全连接神经网络、循环神经网络等。深度学习模型可以包括但不限于深度信念网络(Deep belief Networks,DBN)、卷积神经网络(Convolution Neural Networks,CNN)、递归神经网络(RecurrentNeural Networks,RNN)、循环神经网(Recurrent Neural Networks,RNN)、生成式对抗网络(GAN,GenerativeAdversarial Networks)等等。在本发明的具体实施例中,第一训练模型的网络结构可为2DU-net、3DU-net、U-net++、U-net3+、V-net中的任意一种。当然,在其它的实施例中,还可采用其它的模型来实现本发明功能,本发明对此不做限定。
应当理解的是,深度学习模型可以是能够从已有数据(历史医疗图像)中学习到一定的知识和能力用于处理新数据,并可以被设计用于执行各种任务,在本发明的实施例中用于进行PET重建图像的生成。
需要说明的是,本发明实施例提供的PET成像方法可以用于实现动态图像的重建,动态图像可以看成是由多张静态图像按照时序连在一起形成的,因此在进行动态图像的重建时,将每一帧反投影图像以及静态图像依次输入第二训练模型中,得到的多帧PET重建图像即为动态重建图像。
基于上述散射校正方法,本发明实施例还相应的提供一种散射校正装置300,请参阅图3,该散射校正装置300包括初始数据获取模块310、第一反投影图像获取模块320、散射估计获取模块330和散射校正模块340。
初始数据获取模块310用于获取目标对象在扫描期间的TOF-PET数据和衰减图像。
第一反投影图像获取模块320用于基于TOF-PET数据获取TOF信息的反投影图像。
散射估计获取模块330用于基于反投影图像和衰减图像,获取TOF-PET成像的散射估计。
散射校正模块340用于基于TOF-PET成像的散射估计对TOF-PET数据进行散射校正。
本实施例中,首先获取目标对象在扫描期间的TOF-PET数据和衰减图像,然后对TOF-PET数据进行反投影处理得到TOF信息的反投影图像,之后再根据该反投影图像与衰减图像获取散射估计,最后根据散射估计对TOF-PET数据进行散射校正,从而能够有效的减少用于散射估计的迭代次数,在加快散射校正速度的同时,保证散射估计的准确性。
在一些实施例中,散射估计获取模块包括示踪剂分布图像获取单元和散射估计单元。
示踪剂分布图像获取单元用于将衰减图像和反投影图像输入至训练完备的第一训练模型中,以得到示踪剂分布图像;
散射估计单元用于基于示踪剂分布图像和衰减图像,获取TOF-PET成像的散射估计。
在一些实施例中,散射校正装置300还包括第一训练模块,第一训练模块用于:
获取第一训练集并构建第一初始训练模型,其中,第一训练集包括若干个第一目标参数以及若干组与第一目标参数对应的第一输入参数,第一输入参数为反投影图像以及与其对应的衰减图像,第一目标参数为示踪剂分布图像;
采用第一训练集对所述第一初始训练模型进行训练,以得到训练完备的第一训练模型。
在一些实施例中,第一训练模型为机器学习模型或者深度学习模型。
在一些实施例中,散射校正装置300还包括预校正模块,预校正模块用于对TOF-PET数据进行预校正处理。
基于上述散射校正方法、装置以及PET成像方法,本发明实施例还相应的提供一种PET成像装置400,请参阅图4,该PET成像装置400包括静态图像获取模块410、第二反投影图像获取模块420和重建模块430。
静态图像获取模块410用于获取目标对象在扫描期间的TOF-PET数据和静态图像;
第二反投影图像获取模块420用于基于TOF-PET数据获取TOF信息的反投影图像。
重建模块430用于将静态图像和反投影图像输入训练完备的第二训练模型中,以得到PET重建图像。
本实施例中,静态图像获取模块410既可通过上述各实施例中所述的散射估计校正进行TOF-PET数据的散射校正后获取静态图像,亦可通过传统方法来获取静态图像。本实施例中,采用经过散射校正后获取的静态图像和反投影图像结合进行PET重建成像,经过散射校正后获取的静态图像给动态的TOF-PET数据提供了散射校正核,因此,在进行PET成像时,可以加快图像重建的速度,而且通过对第二训练模型的训练,在进行重建时,可直接将静态图像和反投影图像输入至第二训练模型中后,得到PET重建图像,进而可以进一步加快图像重建的速度。
在一些实施例中,PET成像装置400还包括第二训练模块,第二训练模块用于:
获取第二训练集并构建第二初始训练模型,其中,第二训练集包括若干个第二目标参数以及若干组与第二目标参数对应的第二输入参数,第二输入参数为静态图像以及反投影图像,第二目标参数为PET重建图像;
采用第二训练集对第二初始训练模型进行训练,以得到训练完备的第二训练模型。
在一些实施例中,第二训练模型为机器学习模型或者深度学习模型。
如图5所示,基于上述散射校正方法、PET成像方法,本发明还相应提供了一种电子设备,该电子设备可以是医学成像系统的控制装置、移动终端、桌上型计算机、笔记本、掌上电脑及服务器等计算设备。该电子设备包括处理器10、存储器20及显示器30。图5仅示出了电子设备的部分组件,但是应理解的是,并不要求实施所有示出的组件,可以替代的实施更多或者更少的组件。
存储器20在一些实施例中可以是该电子设备的内部存储单元,例如电子设备的硬盘或内存。存储器20在另一些实施例中也可以是电子设备的外部存储设备,例如电子设备上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,存储器20还可以既包括电子设备的内部存储单元也包括外部存储设备。存储器20用于存储安装于电子设备的应用软件及各类数据,例如安装电子设备的程序代码等。存储器20还可以用于暂时地存储已经输出或者将要输出的数据。在一实施例中,存储器20上存储有计算机程序40,该计算机程序40可被处理器10所执行,从而实现本申请各实施例的散射校正方法或者PET成像方法。
处理器10在一些实施例中可以是一中央处理器(Central Processing Unit,CPU),微处理器或其他数据处理芯片,用于运行存储器20中存储的程序代码或处理数据,例如执行散射校正方法或者PET成像方法等。
显示器30在一些实施例中可以是LED显示器、液晶显示器、触控式液晶显示器以及OLED(Organic Light-Emitting Diode,有机发光二极管)触摸器等。显示器30用于显示在所述电子设备的信息以及用于显示可视化的用户界面。电子设备的部件处理器10、存储器20及显示器30通过系统总线相互通信。
在一实施例中,当处理器10执行存储器20中计算机程序40时实现以下步骤:
获取目标对象在扫描期间的TOF-PET数据和衰减图像;
基于TOF-PET数据获取TOF信息的反投影图像;
基于反投影图像和衰减图像,获取TOF-PET成像的散射估计;
基于TOF-PET成像的散射估计对TOF-PET数据进行散射校正。
在一实施例中,当处理器10执行存储器20中计算机程序40时还实现以下步骤:
将衰减图像和反投影图像输入至训练完备的第一训练模型中,以得到示踪剂分布图像;
基于示踪剂分布图像和衰减图像,获取TOF-PET成像的散射估计。
在一实施例中,当处理器10执行存储器20中计算机程序40时还实现以下步骤:
获取第一训练集并构建第一初始训练模型,其中,所述第一训练集包括若干个第一目标参数以及若干组与第一目标参数对应的第一输入参数,第一输入参数为反投影图像以及与其对应的衰减图像,第一目标参数为示踪剂分布图像;
采用第一训练集对第一初始训练模型进行训练,以得到训练完备的第一训练模型。
在一实施例中,第一训练模型为机器学习模型或者深度学习模型。
在一实施例中,当处理器10执行存储器20中计算机程序40时还实现以下步骤:
对TOF-PET数据进行预校正处理。
在一实施例中,当处理器10执行存储器20中计算机程序40时还实现以下步骤:
获取目标对象在扫描期间的TOF-PET数据和静态图像;
基于TOF-PET数据获取TOF信息的反投影图像;
将静态图像和反投影图像输入训练完备的第二训练模型中,以得到PET重建图像。
在一实施例中,当处理器10执行存储器20中计算机程序40时还实现以下步骤:
获取第二训练集并构建第二初始训练模型,其中,所述第二训练集包括若干个第二目标参数以及若干组与第二目标参数对应的第二输入参数,第二输入参数为静态图像以及反投影图像,第二目标参数为PET重建图像;
采用第二训练集对所述第二初始训练模型进行训练,以得到训练完备的第二训练模型。
在一实施例中,第二训练模型为机器学习模型或者深度学习模型。
综上所述,本发明提供的散射校正方法、PET成像方法、装置、设备及存储介质,首先获取目标对象在扫描期间的TOF-PET数据和衰减图像,然后对TOF-PET数据进行反投影处理得到TOF信息的反投影图像,之后利用训练好的第一训练模型对衰减图像和反投影图像进行计算,得到示踪剂分布图像,此示踪剂分布图像即为经过初步散射校正后的真实示踪剂分布图像,之后再根据该示踪剂分布图像与衰减图像获取散射估计,最后根据散射估计对TOF-PET数据进行散射校正,从而能够有效的减少用于散射估计的迭代次数,在加快散射校正速度的同时,保证散射估计的准确性。在PET成像时,采用经过散射校正的静态图像和反投影图像结合进行PET重建成像,经过散射校正的静态图像给动态的TOF-PET数据提供了散射校正核,因此,在进行PET成像时,可以加快图像重建的速度,而且通过对第二训练模型的训练,在进行重建时,可直接将静态图像和反投影图像输入至第二训练模型中后,得到PET重建图像,进而可以进一步加快图像重建的速度。
当然,本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关硬件(如处理器,控制器等)来完成,所述的程序可存储于一计算机可读取的存储介质中,该程序在执行时可包括如上述各方法实施例的流程。其中所述的存储介质可为存储器、磁碟、光盘等。
以上所述本发明的具体实施方式,并不构成对本发明保护范围的限定。任何根据本发明的技术构思所做出的各种其他相应的改变与变形,均应包含在本发明权利要求的保护范围内。
Claims (10)
1.一种散射校正方法,其特征在于,包括如下步骤:
获取目标对象在扫描期间的TOF-PET数据和衰减图像;
基于所述TOF-PET数据获取TOF信息的反投影图像;
基于所述反投影图像和所述衰减图像,获取TOF-PET成像的散射估计;
基于所述TOF-PET成像的散射估计对所述TOF-PET数据进行散射校正。
2.根据权利要求1所述的散射校正方法,其特征在于,所述基于所述反投影图像和所述衰减图像,获取TOF-PET成像的散射估计,包括:
将所述衰减图像和所述反投影图像输入至训练完备的第一训练模型中,以得到示踪剂分布图像;
基于所述示踪剂分布图像和所述衰减图像,获取所述TOF-PET成像的散射估计。
3.根据权利要求2所述的散射校正方法,其特征在于,所述方法还包括:
获取第一训练集并构建第一初始训练模型,其中,所述第一训练集包括若干个第一目标参数以及若干组与所述第一目标参数对应的第一输入参数,所述第一输入参数为反投影图像以及与其对应的衰减图像,所述第一目标参数为示踪剂分布图像;
采用所述第一训练集对所述第一初始训练模型进行训练,以得到训练完备的第一训练模型。
4.根据权利要求2或3所述的散射校正方法,其特征在于,所述第一训练模型为机器学习模型或者深度学习模型。
5.根据权利要求1所述的散射校正方法,其特征在于,所述获取目标对象在扫描期间的TOF-PET数据和衰减图像之后,所述基于所述TOF-PET数据获取TOF信息的反投影图像之前还包括:
对所述TOF-PET数据进行预校正处理。
6.一种PET成像方法,其特征在于,包括如下步骤:
获取目标对象在扫描期间的TOF-PET数据和静态图像,其中,所述静态图像通过如权利要求1-5任意一项所述的散射校正方法进行散射校正后获取;
基于所述TOF-PET数据获取TOF信息的反投影图像;
将所述静态图像和所述反投影图像输入训练完备的第二训练模型中,以得到PET重建图像。
7.根据权利要求6所述的PET成像方法,其特征在于,所述方法还包括:
获取第二训练集并构建第二初始训练模型,其中,所述第二训练集包括若干个第二目标参数以及若干组与第二目标参数对应的第二输入参数,所述第二输入参数为静态图像以及反投影图像,所述第二目标参数为PET重建图像;
采用第二训练集对所述第二初始训练模型进行训练,以得到训练完备的第二训练模型。
8.一种散射校正装置,其特征在于,包括:
初始数据获取模块,用于获取目标对象在扫描期间的TOF-PET数据和衰减图像;
反投影图像获取模块,用于基于所述TOF-PET数据获取TOF信息的反投影图像;
散射估计获取模块,用于基于所述反投影图像和所述衰减图像,获取TOF-PET成像的散射估计;
散射校正模块,用于基于所述TOF-PET成像的散射估计对所述TOF-PET数据进行散射校正。
9.一种电子设备,其特征在于,包括:处理器和存储器;
所述存储器上存储有可被所述处理器执行的计算机可读程序;
所述处理器执行所述计算机可读程序时实现如权利要求1-5任意一项所述的散射校正方法或者如权利要求6-7任意一项所述的PET成像方法中的步骤。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如权利要求1-5任意一项所述的散射校正方法或者如权利要求6-7任意一项所述的PET成像方法中的步骤。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210487528.0A CN114862980A (zh) | 2022-05-06 | 2022-05-06 | 散射校正方法、pet成像方法、装置、设备及存储介质 |
PCT/CN2022/141631 WO2023116922A1 (en) | 2021-12-24 | 2022-12-23 | Systems and methods for positron emission tomography imaging |
US18/415,573 US20240153166A1 (en) | 2021-12-24 | 2024-01-17 | Systems and methods for positron emission tomography imaging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210487528.0A CN114862980A (zh) | 2022-05-06 | 2022-05-06 | 散射校正方法、pet成像方法、装置、设备及存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114862980A true CN114862980A (zh) | 2022-08-05 |
Family
ID=82635698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210487528.0A Pending CN114862980A (zh) | 2021-12-24 | 2022-05-06 | 散射校正方法、pet成像方法、装置、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114862980A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023116922A1 (en) * | 2021-12-24 | 2023-06-29 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for positron emission tomography imaging |
-
2022
- 2022-05-06 CN CN202210487528.0A patent/CN114862980A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023116922A1 (en) * | 2021-12-24 | 2023-06-29 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for positron emission tomography imaging |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9155514B2 (en) | Reconstruction with partially known attenuation information in time of flight positron emission tomography | |
CN110151210B (zh) | 一种医学图像处理方法、系统、装置和计算机可读介质 | |
US8457380B2 (en) | PET local tomography | |
JP4965575B2 (ja) | 分布させた反復的画像再構成 | |
CN107638188B (zh) | 图像衰减校正方法及装置 | |
US20230085203A1 (en) | Systems and methods for image reconstruction | |
CN107871331B (zh) | 用于重构发射活动图像的系统和方法 | |
CN108986892B (zh) | 用于确定活度图和衰减图的系统和方法 | |
CN110751702A (zh) | 图像重建方法、系统、装置及存储介质 | |
CN114494479A (zh) | 利用神经网络对低剂量pet图像进行同时衰减校正、散射校正和去噪声的系统和方法 | |
CN114092589B (zh) | 图像重建方法及训练方法、装置、设备及存储介质 | |
US11874411B2 (en) | Estimation of partially missing attenuation in time-of-flight positron emission tomography | |
CN114862980A (zh) | 散射校正方法、pet成像方法、装置、设备及存储介质 | |
US11164344B2 (en) | PET image reconstruction using TOF data and neural network | |
US9269167B2 (en) | SPECT motion-correction | |
US10573029B2 (en) | Fast iterative image reconstruction method for emission tomography | |
US11663758B2 (en) | Systems and methods for motion estimation in PET imaging using AI image reconstructions | |
WO2023116922A1 (en) | Systems and methods for positron emission tomography imaging | |
CN112634147B (zh) | 自监督学习的pet图像降噪方法、系统、装置及介质 | |
US11468607B2 (en) | Systems and methods for motion estimation in PET imaging using AI image reconstructions | |
US20230386036A1 (en) | Methods and systems for medical imaging | |
US11701067B2 (en) | Attenuation correction-based weighting for tomographic inconsistency detection | |
CN113759412B (zh) | 获取束流形状和能量探测单元响应特征的方法、装置 | |
WO2024011037A1 (en) | Systems and methods for accelerating spect imaging | |
CN117541667A (zh) | 基于深度学习的归一化校正方法与系统、存储介质和设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |