CN111080737B - 图像重建方法、装置及pet扫描系统 - Google Patents
图像重建方法、装置及pet扫描系统 Download PDFInfo
- Publication number
- CN111080737B CN111080737B CN201911326028.3A CN201911326028A CN111080737B CN 111080737 B CN111080737 B CN 111080737B CN 201911326028 A CN201911326028 A CN 201911326028A CN 111080737 B CN111080737 B CN 111080737B
- Authority
- CN
- China
- Prior art keywords
- detection
- time resolution
- response
- performance data
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000001514 detection method Methods 0.000 claims abstract description 214
- 239000013078 crystal Substances 0.000 claims abstract description 177
- 230000004044 response Effects 0.000 claims abstract description 97
- 239000000523 sample Substances 0.000 claims description 22
- 230000002123 temporal effect Effects 0.000 claims description 21
- 230000006870 function Effects 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 229940121896 radiopharmaceutical Drugs 0.000 description 1
- 239000012217 radiopharmaceutical Substances 0.000 description 1
- 230000002799 radiopharmaceutical effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
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
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扫描系统。
背景技术
正电子发射计算机断层扫描(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)的时间分辨率表示为基于探测晶体对(i,j)的测试扫描数据计算的时间分辨率表示为bij。
基于上述数据,可以建立以下关系式:
使用上述多个公式,对每个探测晶体进行时间分辨率迭代计算,得到多个关系式,基于所有关系式,得到每个探测晶体的时间分辨率。
针对上述步骤105,PET扫描系统基于每个探测晶体的时间分辨率,确定每条响应线上探测晶体对的时间分辨率。计算公式如下:
其中,σij表示探测晶体对(i,j)的时间分辨率。
实现中,PET扫描系统可以使用公式(5),基于每个探测晶体的时间分辨率,确定扫描桶内所有响应线上探测晶体对的时间分辨率。
或者,在确定每条响应线上探测晶体对的时间分辨率之前,已经确定出穿过线源的每条第一响应线上第一探测晶体对的时间分辨率的情况下,PET扫描系统可以使用公式(5),基于探测晶体的时间分辨率,确定未穿过线源的每条第二响应线上第二探测晶体对的时间分辨率,从而得到扫描桶内每条响应线上探测晶体对的时间分辨率。
针对上述步骤106,在一个实施例中,图5是本申请一示例性实施例示出的一种确定探测器的探测性能数据方法的流程图,参见图5,确定探测器的探测性能数据的操作可以通过下面方式实现:在步骤1061中,针对每条响应线,基于一响应线上探测晶体对的时间分辨率,构建该响应线的核函数,核函数用于确定响应线上不同像素点发生湮灭的概率;在步骤1062中,在该响应线包括多个线段的情况下,确定多个线段中每个线段的中心点坐标xo;在步骤1063中,针对响应线上的每个像素点,确定像素点的第一边界坐标xa和第二边界坐标xb;在步骤1064中,针对每个线段,以(xa,xb)为积分范围,对以(x-xo)为自变量的核函数进行积分,得到该像素点的一个核概率值,核概率值为基于该线段发出的光子确定的、湮灭发生位置被定位在该像素点的符合事件的发生概率。
在TOF-PET扫描系统中,采集的扫描数据包括符合事件中两个光子的接受时间差,即同一湮灭事件所产生的光子对到达探测器的时间差,通过符合事件中两个光子的接受时间差,确定该符合事件对应的湮灭发生位置,计算公式如下:
其中,d表示湮灭发生位置;t1表示符合事件中一个光子的接受时间;t2表示符合事件中另一个光子的接受时间;c表示光速。
实现中,PET扫描系统可以对同一湮灭事件所产生的光子对到达探测器的时间差即飞行时间进行划分,得到多个时间段,基于飞行时间和图像坐标系中飞行位置的映射关系,确定每个时间段对应的线段。
在一个实施例中,假设响应线LORm为一条编号为m的响应线,LORm上探测晶体对(i,j)的时间分辨率为σij,构建LORm的TOF核函数(TOF,kernel),TOF核函数用于确定LORm上不同像素点(可以理解为像素点位置)发生湮灭的概率,其符合高斯分布,公式如下:
其中,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)确定各像素点的灰度值。
其中,表示k次迭代计算后像素点f的灰度值;/>表示(k+1)次迭代计算后像素点f的灰度值;pmf表示基于第一目标光子确定的符合事件的发生概率,第一目标光子为LORm上像素点f发出的且被探测晶体对(i,j)接收的光子;ymn表示基于第二目标光子确定的符合事件的数量,第二目标光子为LORm包括的线段n发出的且被探测晶体对(i,j)接收的光子;smn表示基于第二目标光子确定的随机符合事件和散射符合事件的总数量;B表示扫描桶内目标响应线的编号的集合;C表示一条响应线包括的所有线段的编号的集合。其中,目标响应线可以是扫描桶内经过像素点f的所有响应线,或者,可以是扫描桶内经过像素点f的部分响应线。
公式(9)中,ymn和smn需要根据待扫描对象的扫描数据确定,其他参数的参数值可以通过试验预先得到。
可以使用公式(9)对每个像素点的像素进行迭代计算,最终得到该像素的的灰度值。迭代次数可以根据需要进行设置,例如40-80次。每个像素点的初始灰度值可以均设为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 (8)
1.一种图像重建方法,其特征在于,应用于PET扫描系统,所述PET扫描系统包括扫描桶,所述扫描桶内安装有探测器,所述探测器包括多个探测晶体,所述方法包括:
获取待扫描对象的扫描数据;
获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
基于所述扫描数据和所述探测性能数据重建图像;
所述探测性能数据通过下面方法确定:
获取每个探测晶体的时间分辨率;
基于所述每个探测晶体的时间分辨率,确定所述扫描桶内每条响应线上探测晶体对的时间分辨率;
基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据;
所述基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据,包括:
针对每条响应线,基于所述响应线上探测晶体对的时间分辨率,构建所述响应线的核函数,所述核函数用于确定所述响应线上不同像素点发生湮灭的概率;
在所述响应线包括多个线段的情况下,确定所述多个线段中每个线段的中心点坐标xo;
针对所述响应线上的每个像素点,确定所述像素点的第一边界坐标x和第二a边界坐标xb,x表示所述像素点的横坐标;
针对所述每个线段,以(xa,xb)为积分范围,对以(x-xo)为自变量的核函数进行积分,得到所述像素点的一个核概率值,所述核概率值为在所述像素点处发生湮灭事件的概率,且所述湮灭事件为基于所述线段发出的光子所确定;
所述基于所述扫描数据和所述探测性能数据重建图像,包括:
针对每条响应线上的每个像素点,基于所述像素点的所有核概率值和对应的扫描数据,确定所述像素点的灰度值。
2.根据权利要求1所述的方法,其特征在于,所述获取所述每个探测晶体的时间分辨率,包括:
获取线源置于所述扫描桶内不同位置时,穿过所述线源的所有第一响应线对应的测试扫描数据;
基于所述测试扫描数据,确定所述所有第一响应线中每条第一响应线上第一探测晶体对的时间分辨率;
基于所有第一探测晶体对的时间分辨率,确定所述每个探测晶体的时间分辨率。
3.根据权利要求2所述的方法,其特征在于,所述基于所有第一探测晶体对的时间分辨率,确定所述每个探测晶体的时间分辨率,包括:
基于所述所有第一探测晶体对的时间分辨率,确定每个探测晶体的初始时间分辨率;
基于所有探测晶体的初始时间分辨率和所述所有第一探测晶体对的时间分辨率,确定所述每个探测晶体的目标时间分辨率。
4.根据权利要求3所述的方法,其特征在于,所述基于所述所有第一探测晶体对的时间分辨率,确定每个探测晶体的初始时间分辨率,包括:
针对每个探测晶体,确定所述探测晶体位于N条第一响应线上,N为正整数;
从所述所有第一探测晶体对的时间分辨率中,获取所述N条第一响应线上N个第一探测晶体对的时间分辨率;
基于所述N个第一探测晶体对的时间分辨率和所述N,确定所述探测晶体的初始时间分辨率。
5.根据权利要求4所述的方法,其特征在于,所述基于所有探测晶体的初始时间分辨率和所述所有第一探测晶体对的时间分辨率,确定所述每个探测晶体的目标时间分辨率,包括:
获取N个配对探测晶体的初始时间分辨率,所述第一响应线一端为探测晶体,另一端为配对探测晶体,所述N个配对探测晶体包括:位于所述N条第一响应线中每条第一响应线另一端的探测晶体;
基于所述探测晶体的初始时间分辨率、所述N个配对探测晶体的初始时间分辨率和所述N,对所述探测晶体进行第一次时间分辨率迭代计算,得到所述探测晶体的二次时间分辨率;
利用所述探测晶体的二次时间分辨率、所述N个配对探测晶体的二次时间分辨率、所述N个第一探测晶体对的时间分辨率和所述N,对所述探测晶体进行第二次时间分辨率迭代计算,得到所述探测晶体的三次时间分辨率、所述探测晶体的二次时间分辨率和所述N个配对探测晶体的二次时间分辨率之间的关系式;
按照上述方法,直至对所述探测晶体进行预设次数的时间分辨率迭代计算;
在迭代结束后,基于所有探测晶体对应的关系式,确定所述每个探测晶体的目标时间分辨率。
6.一种图像重建装置,其特征在于,应用于PET扫描系统,所述PET扫描系统包括扫描桶,所述扫描桶内安装有探测器,所述探测器包括多个探测晶体,所述装置包括:
第一获取模块,被配置为获取待扫描对象的扫描数据;
第二获取模块,被配置为获取所述探测器的探测性能数据,所述探测性能数据包括:多个基于所述扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
重建模块,被配置为基于所述扫描数据和所述探测性能数据重建图像;
所述探测性能数据通过下面方法确定:
获取每个探测晶体的时间分辨率;
基于所述每个探测晶体的时间分辨率,确定所述扫描桶内每条响应线上探测晶体对的时间分辨率;
基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据;
所述基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据,包括:
针对每条响应线,基于所述响应线上探测晶体对的时间分辨率,构建所述响应线的核函数,所述核函数用于确定所述响应线上不同像素点发生湮灭的概率;
在所述响应线包括多个线段的情况下,确定所述多个线段中每个线段的中心点坐标xo;
针对所述响应线上的每个像素点,确定所述像素点的第一边界坐标x和第二a边界坐标xb,x表示所述像素点的横坐标;
针对所述每个线段,以(xa,xb)为积分范围,对以(x-xo)为自变量的核函数进行积分,得到所述像素点的一个核概率值,所述核概率值为在所述像素点处发生湮灭事件的概率,且所述湮灭事件为基于所述线段发出的光子所确定;
所述基于所述扫描数据和所述探测性能数据重建图像,包括:
针对每条响应线上的每个像素点,基于所述像素点的所有核概率值和对应的扫描数据,确定所述像素点的灰度值。
7.根据权利要求6所述的装置,其特征在于,所述装置还包括:
第三获取模块,被配置为获取每个探测晶体的时间分辨率;
第一确定模块,被配置为基于所述每个探测晶体的时间分辨率,确定所述扫描桶内每条响应线上探测晶体对的时间分辨率;
第二确定模块,被配置为基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据。
8.一种PET扫描系统,其特征在于,包括:内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于获取待扫描对象的扫描数据;
所述存储器,用于存储确定图像重建对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行所述指令实现如下操作:
获取探测器的探测性能数据,所述探测性能数据包括:多个基于扫描桶内不同响应线上探测晶体对的时间分辨率确定的性能数据;
基于所述扫描数据和所述探测性能数据重建图像;
所述探测性能数据通过下面方法确定:
获取每个探测晶体的时间分辨率;
基于所述每个探测晶体的时间分辨率,确定所述扫描桶内每条响应线上探测晶体对的时间分辨率;
基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据;
所述基于所述每条响应线上探测晶体对的时间分辨率,确定所述探测性能数据,包括:
针对每条响应线,基于所述响应线上探测晶体对的时间分辨率,构建所述响应线的核函数,所述核函数用于确定所述响应线上不同像素点发生湮灭的概率;
在所述响应线包括多个线段的情况下,确定所述多个线段中每个线段的中心点坐标xo;
针对所述响应线上的每个像素点,确定所述像素点的第一边界坐标x和第二a边界坐标xb,x表示所述像素点的横坐标;
针对所述每个线段,以(xa,xb)为积分范围,对以(x-xo)为自变量的核函数进行积分,得到所述像素点的一个核概率值,所述核概率值为在所述像素点处发生湮灭事件的概率,且所述湮灭事件为基于所述线段发出的光子所确定;
所述基于所述扫描数据和所述探测性能数据重建图像,包括:
针对每条响应线上的每个像素点,基于所述像素点的所有核概率值和对应的扫描数据,确定所述像素点的灰度值。
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 CN111080737A (zh) | 2020-04-28 |
CN111080737B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111839569B (zh) * | 2020-07-29 | 2022-11-01 | 明峰医疗系统股份有限公司 | Pet时间分辨率测算、重建系统和方法 |
Citations (3)
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三维图像重建方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10663608B2 (en) * | 2015-09-21 | 2020-05-26 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for calibrating a PET scanner |
-
2019
- 2019-12-20 CN CN201911326028.3A patent/CN111080737B/zh active Active
Patent Citations (3)
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三维图像重建方法和装置 |
Non-Patent Citations (2)
Title |
---|
基于解剖先验信息的图像重建研究;王帅;《中国优秀硕士学位论文全文数据库信息科技》;20180715;论文第20页 * |
基于飞行时间技术的PET发展历史与现状;明凯等;《原子核物理评论》;20120620(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111080737A (zh) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108542412B (zh) | 图像重建方法和pet设备 | |
US9990741B2 (en) | Motion correction in a projection domain in time of flight positron emission tomography | |
JP5530446B2 (ja) | Pet−mrにおける減衰マップの生成方法及び装置 | |
EP3067864B1 (en) | Iterative reconstruction with enhanced noise control filtering | |
JP6440230B2 (ja) | マルチモダリティ撮像のシステム及び方法 | |
JP5172347B2 (ja) | 反復的制約デコンボリューションによる核医学2次元平面画像の復元 | |
US9474495B2 (en) | System and method for joint estimation of attenuation and activity information | |
US9706972B1 (en) | Systems and methods for reconstruction of emission activity image | |
US8903152B2 (en) | Methods and systems for enhanced tomographic imaging | |
CN108986892B (zh) | 用于确定活度图和衰减图的系统和方法 | |
CN110136076B (zh) | 医学扫描成像方法、装置、存储介质及计算机设备 | |
EP2080166A2 (en) | Dirty isotope pet reconstruction | |
US20200105032A1 (en) | Activity-dependent, spatially-varying regularization parameter design for regularized image reconstruction | |
JP6123652B2 (ja) | 散乱成分推定方法 | |
CN111080737B (zh) | 图像重建方法、装置及pet扫描系统 | |
US11164344B2 (en) | PET image reconstruction using TOF data and neural network | |
CN105631908B (zh) | 一种pet图像重建方法和装置 | |
US10217250B2 (en) | Multi-view tomographic reconstruction | |
JP5571578B2 (ja) | マルチセグメント再構成 | |
US11468607B2 (en) | Systems and methods for motion estimation in PET imaging using AI image reconstructions | |
JP2011002306A (ja) | Pet装置の逐次近似画像再構成法 | |
JP6052425B2 (ja) | 輪郭画像生成装置および核医学診断装置 | |
US11663758B2 (en) | Systems and methods for motion estimation in PET imaging using AI image reconstructions | |
CN118351207A (zh) | 基于先验信息的pet图像重建方法、装置、设备及存储介质 | |
CN115731316A (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 | ||
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. |