一种图像重建方法和装置
技术领域
本公开涉及医疗成像技术,特别涉及一种图像重建方法和装置。
背景技术
通过医疗扫描***对被检体(例如,患者)进行扫描,可以发现病灶所在,采取更有针对性的治疗方案。例如,PET(Positron Emission Computed Tomography,正电子发射型计算机断层显像)是当今最先进的大型医疗诊断成像设备之一,PET能够利用活体代谢在分子水平上提供有关脏器及其病变的功能信息,在诊断肿瘤、心血管和神经***等疾病中具有卓越性能。PET的工作原理是在被检体中注射含有放射性核素的药物,放射性核素发生衰变产生正电子,正电子与周围的负电子湮灭产生一对背靠背的伽马光子,伽马光子穿过被检体后被探测器接收并记录,PET可以根据探测到的符合数据进行图像重建,得到反映人体各组织代谢情况的图像。
而随机符合数据是PET符合数据的其中一种主要噪声,图像重建时可以获得随机符合数据,对符合数据进行校正。常见的随机符合数据获得方法包括延迟符合窗法,但是目前的延迟符合窗法获得的随机符合数据并不准确,从而会影响最终的图像重建质量。
发明内容
有鉴于此,本公开提供一种图像重建方法和装置,以提高PET图像重建的图像质量。
具体地,本公开是通过如下技术方案实现的:
第一方面,提供一种图像重建方法,所述方法包括:
获取扫描被检体得到的即时符合数据、延迟符合数据和单事件计数;
根据所述即时符合数据、延迟符合数据和单事件计数,确定真符合单事件影响参数,所述真符合单事件影响参数用于表示由所述即时符合数据中的真符合单事件参与延迟符合引起的延迟符合数据偏差;
利用所述真符合单事件影响参数,对所述延迟符合数据进行校正,并迭代获得校正后的随机符合数据;
根据所述随机符合数据以及扫描被检体得到的扫描数据进行图像重建。
第二方面,提供一种图像重建设备,所述设备包括:存储器、处理器,以及存储在存储器上并可在处理器上运行的计算机指令,所述处理器执行指令时实现以下步骤:
获取扫描被检体得到的即时符合数据、延迟符合数据和单事件计数;
根据所述即时符合数据、延迟符合数据和单事件计数,确定真符合单事件影响参数,所述真符合单事件影响参数用于表示由所述即时符合数据中的真符合单事件参与延迟符合引起的延迟符合数据偏差;
利用所述真符合单事件影响参数,对所述延迟符合数据进行校正,并迭代获得校正后的随机符合数据;
根据所述随机符合数据以及扫描被检体得到的扫描数据进行图像重建。
第三方面,提供一种计算机可读存储介质,其具有存储在其上的指令,当由一个或多个处理器执行时,使得一个或多个处理器执行数据处理方法,该方法包括:
获取扫描被检体得到的即时符合数据、延迟符合数据和单事件计数;
根据所述即时符合数据、延迟符合数据和单事件计数,确定真符合单事件影响参数,所述真符合单事件影响参数用于表示由所述即时符合数据中的真符合单事件参与延迟符合引起的延迟符合数据偏差;
利用所述真符合单事件影响参数,对所述延迟符合数据进行校正,并迭代获得校正后的随机符合数据;
根据所述随机符合数据以及扫描被检体得到的扫描数据进行图像重建。
本公开提供的图像重建方法和装置,通过利用即时符合数据、延迟符合数据、单事件计数,估计出真事件参与延迟符合导致的随机符合数据存在的误差,通过获取真符合单事件影响参数,并根据该参数对延迟符合数据进行校正,获取准确的随机符合数据,从而使得重建图像的质量得到提高。
附图说明
图1是本公开一示例性实施例示出的一种PET***的***架构图;
图2是本公开一示例性实施例示出的一种图像重建方法的流程图;
图3是本公开一示例性实施例示出的延迟符合原理示意图;
图4是本公开一示例性实施例示出的一种真事件开窗延迟符合形成示意图;
图5a是本公开一示例性实施例示出的一种单事件开窗延迟符合形成示意图;
图5b是本公开一示例性实施例示出的一种单事件开窗延迟符合形成示意图;
图5c是本公开一示例性实施例示出的一种多延迟符合形成示意图;
图6是本公开一示例性实施例示出的晶体组块符合示意图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
图1示例了PET***的***架构图,如图1所示,PET***可以包括:承载被检体的扫描床11、探测器12、探测处理电路13和图像重建设备14。其中,当对被检体扫描时,可以由扫描床11承载带动注射有放射性核素药物的被检体,沿着图1中虚线所示的箭头方向运动到探测器12,可以将探测器12在扫描时接收到伽马光子的事件,称为单事件。探测处理电路13可以根据接收到的多个单事件进行符合判定,得到PET符合事件,并由图像重建设备14根据符合事件重建图像。
符合事件是指在PET***中给定的一个小的时间段内,获得两个或两个以上的事件,这个给定的时间段称为符合时间窗,一般为几纳秒。符合事件是两个事件称为两符合,两个以上的事件称为多符合,通常只考虑两符合,而可以将多符合进行丢弃不参与重建。对于两符合事件,如果两个事件是由同一个正电子湮灭产生且未经过散射,称为真符合;如果两个事件是由同一个正电子湮灭产生但是至少有一个事件发生了散射,称为散射符合;如果两个事件是由不同的正电子湮灭产生,称为随机符合。所有的符合事件可以称为即时符合数据,所有真符合事件可以称为真符合数据,所有随机符合事件可以称为随机符合数据,所有散射符合事件可以称为散射符合数据。
图像重建设备14在进行图像重建时,实际需要的数据可以是真符合数据,而随机符合数据和散射符合数据可以称为噪声数据,可以根据噪声数据对即时符合数据进行校正,以使得图像重建的质量更高。例如,上述的随机符合数据即为一种主要噪声,可以在即时符合数据中去掉随机符合数据,使得参与图像重建的数据更加准确。本公开提供的图像重建方法,将主要描述如何获得准确的随机符合数据,并根据该随机符合数据进行图像重建。
本公开的例子可以使用延迟符合窗法获取随机符合数据。如图1所示,在PET***的探测处理电路13中可以包括即时符合电路131、延迟符合电路132和单事件计数装置133。探测处理电路13可以根据采集的被检体扫描数据,获取到随机符合数据计算过程中所需要的参数,例如,即时符合数据和延迟符合数据。图像重建设备14可以根据探测处理电路13得到的参数获取随机符合数据,并利用随机符合数据进行PET图像重建。
如下的例子将描述在一次扫描过程中,如何获得随机符合数据并进行图像重建的过程。该处理过程可以由图像重建设备执行。其中,在后续的描述中涉及到的“真符合”数据,可以是真符合数据和散射符合数据之和,即所提到的“真符合”实际包括同一个正电子湮灭产生且未经过散射的事件、以及由同一个正电子湮灭产生但是至少有一个事件发生了散射的事件。
参见图2示例的图像重建方法的流程,该方法可以由图像处理设备14执行。该图像处理设备14可以包括存储器、处理器,以及存储在存储器上并可在处理器上运行的计算机指令,所述处理器执行指令时可以实现以下处理:
在步骤201中,探测处理电路获取即时符合数据、延迟符合数据和单事件计数。
本步骤中,探测器12可以在扫描被检体时接收到伽马光子,称为单事件。探测处理电路13可以通过即时符合电路131进行符合判定获取即时符合数据,本例子中的即时符合数据可以是两符合数据,可以用P表示。单事件计数装置133可以获得单事件计数,可以用S表示。
延迟符合电路132可以用来获取延迟符合数据。本例子中,延迟符合电路132可以使用单路开窗模式只保留两符合数据。如图3所示的延迟符合窗法符合示意图,即时符合和延迟符合都是在两列事件序列中判定两符合事件,其中,图3中所示的“即时符合”是通过即时符合电路获取符合数据,并且该符合数据是两符合的符合数据,即两个单事件形成的符合事件,而将多于两个单事件形成的多符合数据丢弃。而延迟符合是将即时符合中两列事件序列的其中一列延迟一定时间(例如,远大于符合时间窗)后进行符合判定,延迟了的事件列称为延迟事件列,该列中的每个事件称为延迟事件,另一列未延迟的事件列称为即时事件列,其中的每个事件称为即时事件,在延迟事件列和即时事件列之间仍然按照符合时间窗进行符合判定,得到延迟符合数据,延迟符合数据也可以是两符合数据。
参见图3所示,单路开窗模式可以只通过延迟事件或者即时事件来开窗,所述的开窗即由该事件开始划定符合时间窗的时间长度。本例子可以通过延迟事件开窗为例,图3中的圆点事件为延迟事件(实际实施中,延迟事件可以通过延迟标签标识),倒三角事件为即时事件。例如,将即时事件列中的即时事件21延迟一定时间后,可以得到对应的延迟事件22;又例如,将即时事件23延迟一定时间后,可以得到对应的延迟事件24。假设以延迟事件24来开窗,如图3所示意,由延迟事件24开始划箭头线表示形成一个以延迟事件24为起始的符合时间窗,如果有一个单事件落入该符合时间窗则形成符合。比如图3中存在一个落入时间窗的即时事件25,则延迟事件24与即时事件25形成两符合。
延迟符合电路132可以按照上述图3所示的方法,获得延迟事件列和即时事件列之间的两符合数据,即为延迟符合数据,可以用R表示。
在步骤202中,根据即时符合数据、延迟符合数据和单事件计数,确定真符合单事件影响参数,该真符合单事件影响参数用于表示由即时符合数据中的真符合单事件参与延迟符合引起的延迟符合数据偏差。
在步骤203中,利用真符合单事件影响参数,对采集到的延迟符合数据进行校正,并迭代获得校正后的随机符合数据。
本例子中,通过延迟符合电路132采集两符合事件,即延迟符合数据。但是延迟符合数据存在偏差,并不是准确的随机符合数据。由于在延迟符合电路132中将事件进行了一定时间的延迟,使得真符合事件被分开,以独立单事件参与延迟符合的判定,因此延迟符合数据中可能存在真符合事件中的单事件参与,导致延迟符合数据并不是准确的随机符合数据。这种偏差是由真符合事件中的单事件参与引起,如果将这部分数据进行校正,将使得随机符合数据更为准确。
可以通过上述的步骤202和步骤203,对延迟符合数据进行真符合事件中的单事件(可以简称:真符合单事件)引入的偏差校正。如下说明真符合单事件参与延迟符合的几种情况:
根据图3所示例的延迟符合判定原理,延迟符合是由“即时事件列”和“延迟事件列”中的其中一路事件开窗,判定另一路事件是否落入在开窗形成的符合时间窗内,落入即形成延迟符合。那么,真符合单事件参与延迟符合时,要么可以作为开窗的单事件,要么可以作为落入的单事件。
情况一:真符合单事件开窗;
即时符合电路中开窗的事件同样是延迟符合电路中的开窗事件,即假设一个单事件在即时符合电路中是开窗的事件,那么该单事件在延迟符合电路延迟后,对应的延迟事件也是延迟符合电路的开窗事件。所述的“开窗”,即由该单事件的时间开始划定一个符合时间窗,判定其他单事件是否落入该时间窗形成符合,也就是说,开窗的事件对应的时间是符合时间窗的起点时间。
举例来说:请参见图3的示例,在即时符合电路中,真事件21与真事件26形成两符合,并且该两符合是由真事件21开窗,即由真事件21的时间开始划定一个符合时间窗(如图3所示的箭头线划定范围即时间窗),上述的真事件26落入该符合时间窗内,与真事件21形成两符合。那么,如上所述,真事件21在即时符合电路开窗,当将真事件21延迟以后,得到对应的延迟事件22,也同样将由该延迟事件22开窗。同理,图3所示的即时符合电路中的单事件23是开窗事件,尽管在即时符合电路并没有其他单事件落入时间窗,该单事件23延迟后在延迟符合电路对应的延迟事件24,也是开窗事件。
在即时符合电路中,真符合事件的形成,是由真事件开窗且真事件落入,因此,根据前述的“即时符合电路中开窗的事件同样是延迟符合电路中的开窗事件”,那么即时符合电路中开窗的真事件,其对应的延迟事件在延迟符合电路中也会开窗,但是这个延迟事件其实是即时符合电路中的真事件,所以称为真事件开窗。例如,图4示意了情况一的延迟符合形成状态,延迟事件列中由一个真符合单事件(可以简称真事件)开窗,即时事件列中的一个单事件(没有形成符合的事件,例如图3所示的单事件23)或者真事件(在即时符合电路中形成符合的单事件)落入。本例子可以计算这部分由真符合单事件开窗形成的延迟符合数据。
为了计算延迟事件列中的真事件开窗形成的延迟符合数据,首先可以根据步骤201中采集到的即时符合数据P和单事件计数S,除以本次扫描被检体的扫描时间,得到即时符合数据计数率p以及单事件计数率s。其中,单事件计数率s可以分成sS和sT,sS表示独立单光子计数率,sT表示真符合单光子计数率,s=sS+sT,sT=2t,t是真符合计数率。如果将本例子最终计算得到的实际的随机符合数据计数率设定为rReal,那么存在如下关系:
t=p-rReal;……………………(1)
sS=s-2t;……………………(2)
其中,上述的公式(1)和公式(2)可以参与后续计算中的迭代流程,在迭代的初始可以设置rReal等于r,r即为将延迟符合数据R除以本次扫描被检体的扫描时间得到的延迟符合数据计数率。
根据泊松分布概率关系,可以获得情况一中的真符合单事件开窗形成的延迟符合数据的计数率rT。需要说明的是,本例子的方法在后续的描述中,可以先计算相关数据的计数率,最后乘以扫描时间即可得到对应的数据。比如,如下的公式(3)中的rT实际上是真符合单事件开窗形成的真事件开窗延迟数据的数据计数率,如果将该rT与扫描时间相乘即可得到真事件开窗延迟数据,但是为了描述简便,后续有些地方可以描述为真事件开窗延迟数据rT。
rT=t*τ*(sS+t)*exp(-τ*(sS+t))……………………(3)
在对延迟符合数据进行校正时,可以将真事件开窗延迟数据从延迟符合数据中减去。
情况二:真符合单事件落入;
本情况中,延迟事件列是由独立单事件开窗,但是,落入符合时间窗的即时事件列中的事件,可能是独立单事件,也可能是真符合事件中的单事件。如图5a和图5b的示例,其中,图5a示意由延迟事件列中的单事件开窗且即时事件列中的单事件落入形成延迟符合,图5b示意由延迟事件列中的单事件开窗且即时事件列中的真符合单事件落入形成延迟符合(可以是真符合中的任一单事件,图中以真事件表示)。
图5a和图5b示例的延迟符合数据,可以是实际获取到的单事件开窗延迟数据。但是,该单事件开窗延迟数据并不完全是准确的随机符合数据,比如,图5a中的延迟符合是由两个独立单事件形成,是随机符合数据;然而图5b中的延迟符合是由即时事件列中的真事件落入形成,是真事件参与延迟符合的误差数据,不是真正的随机符合数据。
另外,还有一种延迟符合也是由单事件开窗且真事件落入,但是却由于同时有真事件和单事件落入形成多符合,而致使将该多符合数据丢弃,从而也将多符合中包含的两个真正的单事件形成的延迟符合丢掉,可以将这部分丢掉的数据找回以得到更为准确完整的随机符合数据。如图5c示例的延迟符合中的多符合,是由延迟事件列中的单事件开窗且即时事件列中的单事件和真事件落入。
由上述描述可以看到,真符合单事件落入形成的延迟符合数据,一方面可能造成增加了一些本来不是随机符合数据的符合(例如,图5b所示),另一方面还可能造成丢失了一些本来是真正的随机符合数据的符合(例如,图5c所示),因此,可以对该情况二中获取到的单事件开窗延迟数据进行校正。本例子可以将用于对单事件开窗延迟数据进行校正的参数称为有效随机缩放系数,该系数可以表示为β。如下所示的公式(4),实际的随机符合数据可以是由延迟符合电路获取到的延迟符合数据r减去真事件开窗延迟数据rT,并通过有效随机缩放系数β进行缩放校正即可得到。以这些数据的计数率表示为:
rReal=β*(r-rT)…………………………(4)
有效随机缩放系数的获得,可以将上面的图5a至图5c所示的三种延迟符合情况的数据分别计算出,再结合该三种数据得到。
对于图5a,可以根据真符合计数率t、独立单光子计数率sS,利用泊松分布概率关系,获得单事件开窗延迟数据中独立单事件落入的延迟符合数据,可以称为第一延迟符合数据。如下公式(5)所示:
rSS=sSO*τ*sS*exp(-τ*sS)*exp(-τ*t)…………………………(5)
其中,sSO是假设的单事件开窗的计数率,τ是符合时间窗。
对于图5b,可以根据真符合计数率t、独立单光子计数率sS,利用泊松分布概率关系,获得单事件开窗延迟数据中真符合单事件落入(有一个真符合单事件落入)的延迟符合数据,可以称为第二延迟符合数据,如下公式(6)所示:
对于图5c,可以根据真符合计数率t、独立单光子计数率sS,利用泊松分布概率关系,获得由于真符合单事件和独立单事件均落入形成多延迟符合的数据,可以称为第三延迟符合数据。这种情况是本来可以形成一个随机符合,但是由于有真符合单事件落入符合事件窗,变成多符合,从而导致这一应有的随机符合事件被丢弃。如下公式(7)所示:
rSTS=sSO*τ*sS*exp(-τ*sS)*[1-exp(-τ*t)]…………………(7)
根据上述公式(5)至公式(7),推导有效随机缩放系数如下:根据上述对情况一和情况二的说明,延迟符合电路实际采集到的延迟符合数据r=rT+rSS+rST,而实际准确的随机符合数据可以是将该延迟符合数据中减去情况一的数据,并对情况二的数据校正后得到,rReal=rSS+rSTS。
根据r=rT+rSS+rST和rReal=rSS+rSTS,也可以推导出上述的有效随机缩放系数β的计算公式:
rReal/(r-rT)=(rSS+rSTS)/(rSS+rST);
从而得到:rReal=(rSS+rSTS)/(rSS+rST)*(r-rT)=β*(r-rT)。
即得到上述的公式(4):rReal=β*(r-rT);可以将情况一的真事件开窗延迟数据由实际采集的延迟符合数据r中减去,得到单事件开窗延迟数据r-rT,并且由有效随机缩放系数β再对该实际采集的单事件开窗延迟数据进行缩放校正,得到随机符合数据rReal。
其中,还可以得到有效随机缩放系数β=(rSS+rSTS)/(rSS+rST)…………(8)。
在公式(8)中,作为分母的第一延迟符合数据rSS和第二延迟符合数据rST之和,可以是实际采集到的单事件开窗延迟数据,即是通过延迟符合电路实际获取到的由单事件开窗形成的延迟符合数据;而作为分子的第一延迟符合数据和第三延迟符合数据都是由真符合单事件落入时间窗形成的延迟符合数据,包括采集到的两符合和丢弃掉的多符合,这部分数据可以称为真延迟符合数据。有效随机缩放系数β可以是真延迟符合数据与单事件开窗延迟数据之比。
综上,可以按照如下流程迭代计算rReal
第一步:令rReal=r
迭代
其中,上述迭代流程中的有效随机缩放系数β在计算时,由于rSS、rST、rSTS的公式中都包含单事件开窗计数率sSO,而该sSO是假设的未知数,根据β的计算公式,可以在计算β时将该参数sSO去掉,不影响β的计算。将sSO去掉后,可以将原公式(5)、(6)、(7)更改为上述迭代流程中的αSS、αST、αSTS。
通常迭代三次即可获得准确解,即求得准确的随机符合数据的计数率rreal,再乘上扫描时间获得准确的随机符合数据总值RReal。
此外,上述对延迟符合数据进行校正的例子,是以同时考虑了情况一和情况二的延迟符合形成情况为例,在其他的例子中,也可以结合部分情况对延迟符合数据进行校正。比如,可以根据情况一进行校正,或者单独根据情况二进行校正,又或者还可以结合其他符合形式进行校正。不论何种方式,都是根据真符合单事件影响参数,校正由于真事件引入延迟符合带来的数据误差。
在步骤204中,根据所述随机符合数据、以及扫描所述被检体得到的扫描数据进行图像重建。
例如,在利用随机符合数据进行图像重建时,可以包括两种方式。一种方式是根据所述随机符合数据,对扫描所述被检体得到的扫描数据进行随机校正,并利用校正后的扫描数据进行图像重建;另一种方式也可以把随机符合数据作为重建数据的一部分进行图像重建。
本步骤中,可以对步骤203得到的随机符合数据做平滑处理,获得低噪声随机符合数据。例如可以利用如下公式:
其中A和B为两个连续的晶体组块,如图6所示,A和B之间的任意两个晶体都能形成符合,R
lk表示晶体l和k之间随机计数,
然后可以利用所述随机符合数据,经过重建工作进行图像重建。
本公开通过利用即时符合数据、延迟符合数据、单事件计数,估计出真事件参与延迟符合导致的随机符合数据存在的误差,通过获取真符合单事件影响参数,并根据该参数对延迟符合数据进行校正,获取准确的随机符合数据,从而使得重建图像的质量得到提高。
本公开的数据处理方法的功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本公开的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台控制和处理设备执行本公开各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,RandomAccess Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本例子提供的一种计算机可读存储介质,其具有存储在其上的指令,当由一个或多个处理器执行时,使得一个或多个处理器执行数据处理方法,该方法包括:获取扫描被检体得到的即时符合数据、延迟符合数据和单事件计数;
根据所述即时符合数据、延迟符合数据和单事件计数,确定真符合单事件影响参数,所述真符合单事件影响参数用于表示由所述即时符合数据中的真符合单事件参与延迟符合引起的延迟符合数据偏差;
利用所述真符合单事件影响参数,对所述延迟符合数据进行校正,并迭代获得校正后的随机符合数据;
根据所述随机符合数据以及扫描被检体得到的扫描数据进行图像重建。
以上所述仅为本公开的较佳实施例而已,并不用以限制本公开,凡在本公开的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本公开保护的范围之内。