CN111080737A - 图像重建方法、装置及pet扫描系统 - Google Patents

图像重建方法、装置及pet扫描系统 Download PDF

Info

Publication number
CN111080737A
CN111080737A CN201911326028.3A CN201911326028A CN111080737A CN 111080737 A CN111080737 A CN 111080737A CN 201911326028 A CN201911326028 A CN 201911326028A CN 111080737 A CN111080737 A CN 111080737A
Authority
CN
China
Prior art keywords
detection
crystal
time resolution
performance data
detection crystal
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.)
Granted
Application number
CN201911326028.3A
Other languages
English (en)
Other versions
CN111080737B (zh
Inventor
刘勺连
常杰
孙智鹏
李明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang Zhihe Medical Technology Co ltd
Original Assignee
Neusoft Medical Systems Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Neusoft Medical Systems Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN201911326028.3A priority Critical patent/CN111080737B/zh
Publication of CN111080737A publication Critical patent/CN111080737A/zh
Application granted granted Critical
Publication of CN111080737B publication Critical patent/CN111080737B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine (AREA)

Abstract

本申请提供一种图像重建方法、装置及PET扫描系统。所述图像重建方法应用于PET扫描系统,所述PET扫描系统包括扫描桶,扫描桶内安装有探测器,探测器包括多个探测晶体,所述方法包括:获取待扫描对象的扫描数据,获取探测器的探测性能数据,探测性能数据包括:多个基于扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据,基于扫描数据和探测性能数据重建图像。通常探测器包括的多个探测晶体的探测性能不同,本申请实施例在充分考虑上述不同的基础上,基于不同响应线上探测晶体对的时间分辨率确定探测器的探测性能数据,基于探测器的探测性能数据和扫描数据重建图像,得到了高质量的PET图像。

Description

图像重建方法、装置及PET扫描系统
技术领域
本申请涉及医学技术领域,尤其涉及图像重建方法、装置及PET扫描系统。
背景技术
正电子发射计算机断层扫描(Positron Emission Tomography,PET)技术是一种核医学影像技术,在临床医学上具有较高的应用价值,例如可以应用于肿瘤检测、脑血管疾病的诊断等。
相关技术中,将TOF(Time of flight,飞行时间)技术应用到PET技术中,制备出TOF-PET扫描系统。相比于传统的PET扫描系统重建的PET图像,TOF-PET扫描系统重建的PET图像具有信噪比高等优点。
如何对TOF-PET扫描系统的图像重建过程进行优化,以得到更高质量的PET图像,是本领域技术人员有待解决的技术问题。
发明内容
为克服相关技术中存在的问题,本申请提供了一种图像重建方法、装置及PET扫描系统。
第一方面,提供一种图像重建方法,应用于PET扫描系统,所述PET扫描系统包括扫描桶,所述扫描桶内安装有探测器,所述探测器包括多个探测晶体,所述方法包括:
获取待扫描对象的扫描数据;
获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
基于所述扫描数据和所述探测性能数据重建图像。
第二方面,提供一种图像重建装置,应用于PET扫描系统,所述PET扫描系统包括扫描桶,所述扫描桶内安装有探测器,所述探测器包括多个探测晶体,所述方法包括:
第一获取模块,被配置为获取待扫描对象的扫描数据;
第二获取模块,被配置为获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
重建模块,被配置为基于所述扫描数据和所述探测性能数据重建图像。
第三方面,提供一种PET扫描系统,包括:内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于获取待扫描对象的扫描数据;
所述存储器,用于存储图像重建对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行所述指令实现如下操作:
获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
基于所述扫描数据和所述探测性能数据重建图像。
本申请的实施例提供的技术方案可以包括以下有益效果:
本申请实施例提供了一种图像重建方法,获取待扫描对象的扫描数据,获取探测器的探测性能数据,探测性能数据包括:多个基于不同响应线上探测晶体对的时间分辨率确定的性能数据,基于扫描数据和探测性能数据重建图像。通常探测器包括的多个探测晶体的探测性能不同,如探测能量大小不同、能量分辨率不同等,本申请实施例在充分考虑上述不同的基础上,基于不同响应线上探测晶体对的时间分辨率确定探测器的探测性能数据,基于探测器的探测性能数据和扫描数据重建图像,得到了高质量的PET图像。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本申请。
附图说明
此处的附图被并入说明书中并构成本申请的一部分,示出了符合本申请的实施例,并与说明书一起用于解释本申请的原理。
图1是本申请一示例性实施例示出的一种图像重建方法的流程图;
图2是本申请一示例性实施例示出的一种获取所述探测器的探测性能数据方法的流程图;
图3是本申请一示例性实施例示出的线源的放置位置的示意图;
图4是本申请一示例性实施例示出的一种获取每个探测晶体的时间分辨率方法的流程图;
图5是本申请一示例性实施例示出的一种确定探测器的探测性能数据方法的流程图;
图6是本申请一示例性实施例示出的一种图像重建装置的示意图;
图7是本申请一示例性实施例示出的一种PET扫描系统的示意图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
下面结合说明书附图,对本申请实施例进行详细描述。
本申请实施例提供了一种图像重建方法,应用于PET扫描系统。
PET扫描系统可以是应用了TOF技术的PET扫描系统,即TOF-PET扫描系统。
PET扫描系统包括扫描床、扫描桶和控制设备。本申请实施例提供的图像重建方法可以应用于PET扫描系统中的控制设备。
图1是本申请一示例性实施例示出的一种图像重建方法的流程图,该实施例可以包括以下步骤:
在步骤101中,获取待扫描对象的扫描数据。
待扫描对象在进行PET扫描之前,需要注射能够发射正电子的放射性药物。在进行PET扫描过程中,待扫描对象躺在扫描床上,扫描床位于扫描桶内,待扫描患者体内的正电子与电子发生湮灭,产生两个方向相反的伽马光子,扫描桶内壁上环绕设置的探测晶体对这对光子进行探测,得到待扫描对象的扫描数据。
通常,探测器探测到一对背靠背传输的光子后,判断该对光子是否满足时间窗设置,如果满足,则形成一符合事件,符合事件包括真符合事件、随机符合事件和散射符合事件,其中,随机符合事件和散射符合事件是噪声事件。PET扫描系统根据符合事件涉及的光子进行图像重建,得到发射正电子的核素分布图,即PET图像。
在步骤102中,获取探测器的探测性能数据,探测性能数据包括:多个基于扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据。
由于探测晶体并不能识别光子到达探测晶体上的准确时间,而是识别成一个可能存在偏差的时间,这种偏差时间的分布的半高宽称为探测晶体的时间分辨率。一条响应线上探测晶体对的时间分辨率基于该条响应线上两个探测晶体的时间分辨率确定。探测器的探测性能数据能够反映探测器的光子探测性能。
探测器的探测性能数据可以是PET系统事先获取并存储的,在执行本步骤时直接从存储位置获取,或者,探测器的探测性能数据可以是PET系统在获取待扫描对象的扫描数据后实时获取的。事先获取探测器的探测性能数据,可以减小每次图像重建操作的数据处理量,缩短图像重建时间,提高图像重建效率。
在一个实施例中,图2是本申请一示例性实施例示出的一种获取所述探测器的探测性能数据方法的流程图,参见图2,探测器的探测性能数据可以通过下面方法确定:在步骤104中,获取每个探测晶体的时间分辨率;在步骤105中,基于每个探测晶体的时间分辨率,确定扫描桶内每条响应线上探测晶体对的时间分辨率;在步骤106中,基于每条响应线上探测晶体对的时间分辨率,确定探测器的探测性能数据。
针对上述步骤104,在一个实施例中,可以通过实验获得各探测晶体的时间分辨率。
在实验过程中,选取一个线源,线源向四周方向发射光子,将线源放置在扫描桶内中心轴附近的不同位置。例如,图3是本申请一示例性实施例示出的线源的放置位置的示意图,图3中,a为环绕排列的探测晶体,O为扫描桶的中心线,可以将线源放置在环b上的不同位置。
图4是本申请一示例性实施例示出的一种获取每个探测晶体的时间分辨率方法的流程图,参见图4,获取每个探测晶体的时间分辨率的操作可以通过下面方式实现:在步骤1041中,获取线源置于所述扫描桶内不同位置时,穿过线源的所有第一响应线(第一LOR)对应的测试扫描数据;在步骤1042中,基于测试扫描数据,确定所有第一响应线中每条第一响应线上第一探测晶体对的时间分辨率;在步骤1043中,基于所有第一探测晶体对的时间分辨率,确定每个探测晶体的时间分辨率。
针对步骤1041,在一个实施例中,在线源置于扫描桶内一位置后,针对每个探测晶体,基于该探测晶体的位置和线源的位置,确定一响应线,称为第一响应线,确定位于该第一响应线另一端的配对探测晶体,该探测晶体和配对探测晶体采集的数据为该第一响应线对应的测试扫描数据。
针对步骤1042,针对每条第一响应线,可以基于该条第一响应线对应的测试扫描数据,确定该条第一响应线上第一探测晶体对的时间分辨率。步骤1042的具体计算过程是现有技术,本公开实施例不做赘述。
针对步骤1043,可以基于所有第一探测晶体对的时间分辨率,确定每个探测晶体的初始时间分辨率,并基于所有探测晶体的初始时间分辨率和所有第一探测晶体对的时间分辨率,确定每个探测晶体的目标时间分辨率。
在一个实施例中,基于所有第一探测晶体对的时间分辨率,确定每个探测晶体的初始时间分辨率的操作可以通过下面方式实现:首先,针对每个探测晶体,确定探测晶体位于N条第一响应线上,N为正整数;其次,从所有第一探测晶体对的时间分辨率中,获取N条第一响应线上N个第一探测晶体对的时间分辨率;再次,基于N个第一探测晶体对的时间分辨率和N,确定该探测晶体的初始时间分辨率。
在这种情况下,基于所有探测晶体的初始时间分辨率和所有第一探测晶体对的时间分辨率,确定每个探测晶体的目标时间分辨率的操作可以通过下面方式实现:
第一步骤,获取N个配对探测晶体的初始时间分辨率,N个配对探测晶体包括:位于N条第一响应线中每条第一响应线另一端的探测晶体。
第二步骤,基于探测晶体的初始时间分辨率、N个配对探测晶体的初始时间分辨率和N,对探测晶体进行第一次时间分辨率迭代计算,得到探测晶体的二次时间分辨率。
第三步骤,利用探测晶体的二次时间分辨率、N个配对探测晶体的二次时间分辨率、N个第一探测晶体对的时间分辨率和N,对探测晶体进行第二次时间分辨率迭代计算,得到探测晶体的三次时间分辨率、探测晶体的二次时间分辨率和N个配对探测晶体的二次时间分辨率之间的关系式。
第四步骤,按照上述方法,直至对探测晶体进行预设次数的时间分辨率迭代计算。
第五步骤,在迭代结束后,基于所有探测晶体对应的关系式,确定每个探测晶体的目标时间分辨率。
下面通过以下示例,对PET扫描系统获取每个探测晶体的时间分辨率的过程进行说明。
探测晶体i的时间分辨率表示为σi,探测晶体j的时间分辨率表示为σj,探测晶体对表示为(i,j)探测晶体对(i,j)的时间分辨率表示为
Figure BDA0002328398170000061
基于探测晶体对(i,j)的测试扫描数据计算的时间分辨率表示为bij
基于上述数据,可以建立以下关系式:
Figure BDA0002328398170000062
Figure BDA0002328398170000063
Figure BDA0002328398170000064
标记为ai
Figure BDA0002328398170000065
标记为Sij。可以使用下面公式(3)和公式(4),对每个探测晶体进行时间分辨率迭代计算。
Figure BDA0002328398170000066
其中,
Figure BDA0002328398170000067
表示探测晶体i的初始时间分辨率;N表示探测晶体位于的第一响应线的数量;A表示探测晶体i对应的N个配对探测晶体的晶体标识的集合。
Figure BDA0002328398170000068
其中,
Figure BDA0002328398170000069
表示探测晶体i的k次时间分辨率;
Figure BDA00023283981700000610
表示探测晶体i的(k-1)次时间分辨率;
Figure BDA00023283981700000611
表示探测晶体j的(k-1)次时间分辨率。
使用上述多个公式,对每个探测晶体进行时间分辨率迭代计算,得到多个关系式,基于所有关系式,得到每个探测晶体的时间分辨率。
探测晶体i的时间分辨率表示为
Figure BDA00023283981700000612
可以理解为探测晶体i的目标时间分辨率。其中,K为迭代计算的总次数,可以根据需要进行设置。
针对上述步骤105,PET扫描系统基于每个探测晶体的时间分辨率,确定每条响应线上探测晶体对的时间分辨率。计算公式如下:
Figure BDA0002328398170000071
其中,σij表示探测晶体对(i,j)的时间分辨率。
实现中,PET扫描系统可以使用公式(5),基于每个探测晶体的时间分辨率,确定扫描桶内所有响应线上探测晶体对的时间分辨率。
或者,在确定每条响应线上探测晶体对的时间分辨率之前,已经确定出穿过线源的每条第一响应线上第一探测晶体对的时间分辨率的情况下,PET扫描系统可以使用公式(5),基于探测晶体的时间分辨率,确定未穿过线源的每条第二响应线上第二探测晶体对的时间分辨率,从而得到扫描桶内每条响应线上探测晶体对的时间分辨率。
针对上述步骤106,在一个实施例中,图5是本申请一示例性实施例示出的一种确定探测器的探测性能数据方法的流程图,参见图5,确定探测器的探测性能数据的操作可以通过下面方式实现:在步骤1061中,针对每条响应线,基于一响应线上探测晶体对的时间分辨率,构建该响应线的核函数,核函数用于确定响应线上不同像素点发生湮灭的概率;在步骤1062中,在该响应线包括多个线段的情况下,确定多个线段中每个线段的中心点坐标xo;在步骤1063中,针对响应线上的每个像素点,确定像素点的第一边界坐标xa和第二边界坐标xb;在步骤1064中,针对每个线段,以(xa,xb)为积分范围,对以(x-xo)为自变量的核函数进行积分,得到该像素点的一个核概率值,核概率值为基于该线段发出的光子确定的、湮灭发生位置被定位在该像素点的符合事件的发生概率。
在TOF-PET扫描系统中,采集的扫描数据包括符合事件中两个光子的接受时间差,即同一湮灭事件所产生的光子对到达探测器的时间差,通过符合事件中两个光子的接受时间差,确定该符合事件对应的湮灭发生位置,计算公式如下:
Figure BDA0002328398170000072
其中,d表示湮灭发生位置;t1表示符合事件中一个光子的接受时间;t2表示符合事件中另一个光子的接受时间;c表示光速。
实现中,PET扫描系统可以对同一湮灭事件所产生的光子对到达探测器的时间差即飞行时间进行划分,得到多个时间段,基于飞行时间和图像坐标系中飞行位置的映射关系,确定每个时间段对应的线段。
在一个实施例中,假设响应线LORm为一条编号为m的响应线,LORm上探测晶体对(i,j)的时间分辨率为σij,构建LORm的TOF核函数(TOF,kernel),TOF核函数用于确定LORm上不同像素点(可以理解为像素点位置)发生湮灭的概率,其符合高斯分布,公式如下:
Figure BDA0002328398170000081
其中,
Figure BDA0002328398170000082
c表示光速;x表示图像坐标系中LORm上一像素点的横坐标;g(x)表示横坐标为x的像素点发生湮灭的概率。
Figure BDA0002328398170000083
其中,m表示响应线LORm的编号;n表示LORm包括的线段的编号;f表示像素点的编号;x0表示LORm上像素点f的第一边界坐标;x1表示LORm上像素点f的第二边界坐标;xo表示线段n的中心点坐标;τmn,f表示TOF核概率值,具体表示基于线段n发出的光子确定的、湮灭发生位置被定位在LORm上像素点f的符合事件的发生概率。
针对LORm上的每个像素点,可以得到n个不同的τmn,f
τmn,f表征探测器的探测性能,可作为探测器的探测性能数据使用。
像素点具有一定的尺寸,一条响应线上一个像素点存在两个边界位置,将每个边界位置的坐标称为边界坐标。
步骤103、基于待扫描对象的扫描数据和探测性能数据重建图像。
在一个实施例中,计算出每条响应线上每个像素点的TOF核概率值后,针对每条响应线上的每个像素点,基于该像素点的所有核概率值和对应的扫描数据,确定该像素点的灰度值。
在一个实施例中,在使用上述公式(7)和公式(8),计算出每条响应线上每个像素点的TOF核概率值后,可以利用下面公式(9)确定各像素点的灰度值。
Figure BDA0002328398170000084
其中,
Figure BDA0002328398170000085
表示k次迭代计算后像素点f的灰度值;
Figure BDA0002328398170000086
表示(k+1)次迭代计算后像素点f的灰度值;pmf表示基于第一目标光子确定的符合事件的发生概率,第一目标光子为LORm上像素点f发出的且被探测晶体对(i,j)接收的光子;ymn表示基于第二目标光子确定的符合事件的数量,第二目标光子为LORm包括的线段n发出的且被探测晶体对(i,j)接收的光子;smn表示基于第二目标光子确定的随机符合事件和散射符合事件的总数量;B表示扫描桶内目标响应线的编号的集合;C表示一条响应线包括的所有线段的编号的集合。其中,目标响应线可以是扫描桶内经过像素点f的所有响应线,或者,可以是扫描桶内经过像素点f的部分响应线。
公式(9)中,ymn和smn需要根据待扫描对象的扫描数据确定,其他参数的参数值可以通过试验预先得到。
可以使用公式(9)对每个像素点的像素进行迭代计算,最终得到该像素的的灰度值。迭代次数可以根据需要进行设置,例如40-80次。每个像素点的初始灰度值
Figure BDA0002328398170000091
可以均设为1,或者,可以使用其他重建图像中像素点的灰度值作为本实施例中同位置像素点的初始灰度值。
本申请实施例提供了一种图像重建方法,获取待扫描对象的扫描数据,获取探测器的探测性能数据,探测性能数据包括:多个基于不同响应线上探测晶体对的时间分辨率确定的性能数据,基于扫描数据和探测性能数据重建图像。通常探测器包括的多个探测晶体的探测性能不同,如探测能量大小不同、能量分辨率不同等,本申请实施例在充分考虑上述不同的基础上,基于不同响应线上探测晶体对的时间分辨率确定探测器的探测性能数据,基于探测器的探测性能数据和扫描数据重建图像,得到了高质量的PET图像。
与前述图像重建方法相对应,本申请还提供了图像重建装置及PET扫描系统的实施例。
参照图6是本申请一示例性实施例示出的一种图像重建装置的示意图,应用于PET扫描系统,所述PET扫描系统包括扫描桶,所述扫描桶内安装有探测器,所述探测器包括多个探测晶体,所述装置包括:第一获取模块21、第二获取模块22和重建模块23;其中,
所述第一获取模块21,被配置为获取待扫描对象的扫描数据;
所述第二获取模块22,被配置为获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
所述重建模块23,被配置为基于所述扫描数据和所述探测性能数据重建图像。
在一个可选的实施例中,在图6所示的图像重建装置的基础上,所述装置还可以包括:第三获取模块、第一确定模块和第二确定模块;其中,
所述第三获取模块,被配置为获取所述每个探测晶体的时间分辨率;
所述第一确定模块,被配置为基于所述每个探测晶体的时间分辨率,确定所述扫描桶内每条响应线上探测晶体对的时间分辨率;
所述第二确定模块,被配置为基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据。
参照图7是本申请一示例性实施例示出的一种PET扫描系统的示意图,该设备可以包括:通过内部总线310连接的存储器320、处理器330和外部接口340。
其中,外部接口340,用于获取待扫描对象的扫描数据;
存储器320,用于存储图像重建对应的机器可读指令;
处理器330,用于读取所述存储器320上的所述机器可读指令,并执行所述指令以实现如下操作:
获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
基于所述扫描数据和所述探测性能数据重建图像。
在公开实施例中,计算机可读存储介质可以是多种形式,比如,在不同的例子中,所述机器可读存储介质可以是:RAM(Radom Access Memory,随机存取存储器)、易失存储器、非易失性存储器、闪存、存储驱动器(如硬盘驱动器)、固态硬盘、任何类型的存储盘(如光盘、dvd等),或者类似的存储介质,或者它们的组合。特殊的,所述的计算机可读介质还可以是纸张或者其他合适的能够打印程序的介质。使用这些介质,这些程序可以被通过电学的方式获取到(例如,光学扫描)、可以被以合适的方式编译、解释和处理,然后可以被存储到计算机介质中。
以上所述仅为本申请的较佳实施例而已,并不用以限制本申请,凡在本申请的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本申请保护的范围之内。

Claims (11)

1.一种图像重建方法,其特征在于,应用于PET扫描系统,所述PET扫描系统包括扫描桶,所述扫描桶内安装有探测器,所述探测器包括多个探测晶体,所述方法包括:
获取待扫描对象的扫描数据;
获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
基于所述扫描数据和所述探测性能数据重建图像。
2.根据权利要求1所述的方法,其特征在于,所述探测性能数据通过下面方法确定:
获取所述每个探测晶体的时间分辨率;
基于所述每个探测晶体的时间分辨率,确定所述扫描桶内每条响应线上探测晶体对的时间分辨率;
基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据。
3.根据权利要求2所述的方法,其特征在于,所述获取所述每个探测晶体的时间分辨率,包括:
获取线源置于所述扫描桶内不同位置时,穿过所述线源的所有第一响应线对应的测试扫描数据;
基于所述测试扫描数据,确定所述所有第一响应线中每条第一响应线上第一探测晶体对的时间分辨率;
基于所有第一探测晶体对的时间分辨率,确定所述每个探测晶体的时间分辨率。
4.根据权利要求3所述的方法,其特征在于,所述基于所有第一探测晶体对的时间分辨率,确定所述每个探测晶体的时间分辨率,包括:
基于所述所有第一探测晶体对的时间分辨率,确定每个探测晶体的初始时间分辨率;
基于所有探测晶体的初始时间分辨率和所述所有第一探测晶体对的时间分辨率,确定所述每个探测晶体的目标时间分辨率。
5.根据权利要求4所述的方法,其特征在于,所述基于所述所有第一探测晶体对的时间分辨率,确定每个探测晶体的初始时间分辨率,包括:
针对每个探测晶体,确定所述探测晶体位于N条第一响应线上,N为正整数;
从所述所有第一探测晶体对的时间分辨率中,获取所述N条第一响应线上N个第一探测晶体对的时间分辨率;
基于所述N个第一探测晶体对的时间分辨率和所述N,确定所述探测晶体的初始时间分辨率。
6.根据权利要求5所述的方法,其特征在于,所述基于所有探测晶体的初始时间分辨率和所述所有第一探测晶体对的时间分辨率,确定所述每个探测晶体的目标时间分辨率,包括:
获取N个配对探测晶体的初始时间分辨率,所述N个配对探测晶体包括:位于所述N条第一响应线中每条第一响应线另一端的探测晶体;
基于所述探测晶体的初始时间分辨率、所述N个配对探测晶体的初始时间分辨率和所述N,对所述探测晶体进行第一次时间分辨率迭代计算,得到所述探测晶体的二次时间分辨率;
利用所述探测晶体的二次时间分辨率、所述N个配对探测晶体的二次时间分辨率、所述N个第一探测晶体对的时间分辨率和所述N,对所述探测晶体进行第二次时间分辨率迭代计算,得到所述探测晶体的三次时间分辨率、所述探测晶体的二次时间分辨率和所述N个配对探测晶体的二次时间分辨率之间的关系式;
按照上述方法,直至对所述探测晶体进行预设次数的时间分辨率迭代计算;
在迭代结束后,基于所有探测晶体对应的关系式,确定所述每个探测晶体的目标时间分辨率。
7.根据权利要求2所述的方法,其特征在于,所述基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据,包括:
针对每条响应线,基于所述响应线上探测晶体对的时间分辨率,构建所述响应线的核函数,所述核函数用于确定所述响应线上不同像素点发生湮灭的概率;
在所述响应线包括多个线段的情况下,确定所述多个线段中每个线段的中心点坐标xo
针对所述响应线上的每个像素点,确定所述像素点的第一边界坐标xa和第二边界坐标xb
针对所述每个线段,以(xa,xb)为积分范围,对以(x-xo)为自变量的核函数进行积分,得到所述像素点的一个核概率值,所述核概率值为基于所述线段发出的光子确定的、湮灭发生位置被定位在所述像素点的符合事件的发生概率。
8.根据权利要求7所述的方法,其特征在于,所述基于所述扫描数据和所述探测性能数据重建图像,包括:
针对每条响应线上的每个像素点,基于所述像素点的所有核概率值和对应的扫描数据,确定所述像素点的灰度值。
9.一种图像重建装置,其特征在于,应用于PET扫描系统,所述PET扫描系统包括扫描桶,所述扫描桶内安装有探测器,所述探测器包括多个探测晶体,所述装置包括:
第一获取模块,被配置为获取待扫描对象的扫描数据;
第二获取模块,被配置为获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
重建模块,被配置为基于所述扫描数据和所述探测性能数据重建图像。
10.根据权利要求9所述的装置,其特征在于,所述装置还包括:
第三获取模块,被配置为获取所述每个探测晶体的时间分辨率;
第一确定模块,被配置为基于所述每个探测晶体的时间分辨率,确定所述扫描桶内每条响应线上探测晶体对的时间分辨率;
第二确定模块,被配置为基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据。
11.一种PET扫描系统,其特征在于,包括:内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于获取待扫描对象的扫描数据;
所述存储器,用于存储确定图像重建对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行所述指令实现如下操作:
获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
基于所述扫描数据和所述探测性能数据重建图像。
CN201911326028.3A 2019-12-20 2019-12-20 图像重建方法、装置及pet扫描系统 Active CN111080737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911326028.3A CN111080737B (zh) 2019-12-20 2019-12-20 图像重建方法、装置及pet扫描系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911326028.3A CN111080737B (zh) 2019-12-20 2019-12-20 图像重建方法、装置及pet扫描系统

Publications (2)

Publication Number Publication Date
CN111080737A true CN111080737A (zh) 2020-04-28
CN111080737B CN111080737B (zh) 2023-05-26

Family

ID=70316264

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911326028.3A Active CN111080737B (zh) 2019-12-20 2019-12-20 图像重建方法、装置及pet扫描系统

Country Status (1)

Country Link
CN (1) CN111080737B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111839569A (zh) * 2020-07-29 2020-10-30 明峰医疗系统股份有限公司 Pet时间分辨率测算、重建系统和方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101856236A (zh) * 2010-06-13 2010-10-13 苏州瑞派宁科技有限公司 一种应用适应性的正电子发射断层成像方法及装置
CN102608648A (zh) * 2006-07-28 2012-07-25 皇家飞利浦电子股份有限公司 正电子发射断层摄影中的飞行时间测量
CN104183012A (zh) * 2013-10-31 2014-12-03 上海联影医疗科技有限公司 一种pet三维图像重建方法和装置
US20170082759A1 (en) * 2015-09-21 2017-03-23 Shanghai United Imaging Healthcare Co., Ltd. System and method for calibrating a pet scanner

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608648A (zh) * 2006-07-28 2012-07-25 皇家飞利浦电子股份有限公司 正电子发射断层摄影中的飞行时间测量
CN101856236A (zh) * 2010-06-13 2010-10-13 苏州瑞派宁科技有限公司 一种应用适应性的正电子发射断层成像方法及装置
CN104183012A (zh) * 2013-10-31 2014-12-03 上海联影医疗科技有限公司 一种pet三维图像重建方法和装置
US20170082759A1 (en) * 2015-09-21 2017-03-23 Shanghai United Imaging Healthcare Co., Ltd. System and method for calibrating a pet scanner

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
明凯等: "基于飞行时间技术的PET发展历史与现状", 《原子核物理评论》 *
王帅: "基于解剖先验信息的图像重建研究", 《中国优秀硕士学位论文全文数据库信息科技》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111839569A (zh) * 2020-07-29 2020-10-30 明峰医疗系统股份有限公司 Pet时间分辨率测算、重建系统和方法
CN111839569B (zh) * 2020-07-29 2022-11-01 明峰医疗系统股份有限公司 Pet时间分辨率测算、重建系统和方法

Also Published As

Publication number Publication date
CN111080737B (zh) 2023-05-26

Similar Documents

Publication Publication Date Title
CN108542412B (zh) 图像重建方法和pet设备
CN109009199B (zh) 用于正电子发射断层扫描中的图像数据处理的系统和方法
EP3067864B1 (en) Iterative reconstruction with enhanced noise control filtering
JP6440230B2 (ja) マルチモダリティ撮像のシステム及び方法
US9659389B2 (en) Fast scatter estimation in pet reconstruction
JP4414410B2 (ja) 画像再構成方法
CN107693037B (zh) Pet扫描装置及其时间偏移校正方法
US8575554B2 (en) Methods and systems for scatter estimation in positron emission tomography
CN105431884B (zh) 用于对图像中的级联伽马进行解释和建模的方法
US10438379B2 (en) In-reconstruction filtering for positron emission tomography (PET) list mode iterative reconstruction
JP6335181B2 (ja) ポジトロン放射断層撮影(Positron−EmissionTomography:PET)画像のタイムオブフライト(Time−Of−Flight:TOF)リストモード再構成のためのシステム行列を計算する方法及び装置
CN108986892B (zh) 用于确定活度图和衰减图的系统和方法
US20090057561A1 (en) Dirty isotope pet reconstruction
CN105125231A (zh) 一种pet图像环状伪影的去除方法和装置
CN114666495B (zh) 高分辨率康普顿相机成像方法、装置、电子设备及介质
CN110136076B (zh) 医学扫描成像方法、装置、存储介质及计算机设备
US10482634B2 (en) Systems and methods for imaging with anisotropic voxels
US10222490B2 (en) PET scanner with emission and transmission structures in a checkerboard configuration
US7769217B2 (en) Fast iterative 3D PET image reconstruction using a set of 2D linogram transformations
CN111080737B (zh) 图像重建方法、装置及pet扫描系统
US20140050380A1 (en) Apparatus and method for generating medical image using linear gamma ray source
US11961164B2 (en) Continuous bed motion acquisition with axially short phantom for PET imaging system setup and quality control
KR102283454B1 (ko) 양전자 방출 단층 촬영기의 비행시간 정보를 이용하여 관측된 데이터에 포함된 반응깊이 정보를 보정하여 영상을 재구성하는 영상재구성방법
CN109613593B (zh) 伽马光子作用位置抽样方法、装置、电子设备及存储介质
US11663758B2 (en) Systems and methods for motion estimation in PET imaging using AI image reconstructions

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
GR01 Patent grant
GR01 Patent grant
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230512

Address after: Room 308, No. 177-2 Chuangxin Road, Hunnan District, Shenyang City, Liaoning Province, 110167

Applicant after: Shenyang Zhihe Medical Technology Co.,Ltd.

Address before: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province

Applicant before: Shenyang Neusoft Medical Systems Co.,Ltd.