CN104352244A - 一种数据处理方法和装置 - Google Patents

一种数据处理方法和装置 Download PDF

Info

Publication number
CN104352244A
CN104352244A CN201410562261.2A CN201410562261A CN104352244A CN 104352244 A CN104352244 A CN 104352244A CN 201410562261 A CN201410562261 A CN 201410562261A CN 104352244 A CN104352244 A CN 104352244A
Authority
CN
China
Prior art keywords
correction factor
data
adjustment
immediately
coincidence data
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
CN201410562261.2A
Other languages
English (en)
Other versions
CN104352244B (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 CN201410562261.2A priority Critical patent/CN104352244B/zh
Publication of CN104352244A publication Critical patent/CN104352244A/zh
Application granted granted Critical
Publication of CN104352244B publication Critical patent/CN104352244B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)

Abstract

本发明实施例公开了一种数据处理方法和装置,包括分析弦图的分界,获取所述不包含真符合数据区域中的即时符合数据,通过所述即时符合数据确定出散射分布数据和随机符合数据;对所述即时符合数据进行处理,具体为:利用所述调整后的第二校正因子在包含真符合数据的即时符合数据中确定随机符合数据,可见,通过即时符合数据得到和所述随机符合数据,并为所述随机符合数据设定第二校正因子,通过调整所述第二校正因子以及比较所述第一误差值和所述第二误差值,得到所述调整后的第二校正因子,利用所述调整后的第二校正因子从即时符合数据中确定出所述随机符合数据更加准确。

Description

一种数据处理方法和装置
技术领域
本发明涉及数据处理领域,特别是涉及一种数据处理方法和装置。
背景技术
医疗成像装置根据注入体内的放射性核素在衰变过程中产生的正电子湮灭辐射和符合探测原理可以得到体内的断层影像,比如说典型的正电子发射型计算机断层成像(Positron Emission Computed Tomography,PET)装置。
在正电子湮灭时,会生成几乎背对背运行的两个γ光子。每一个γ光子被PET的探测器探测到就被标记为一个单事件。根据这个单事件发生时所产生的大量数据,可以判断出这个单事件是否是通过正电子湮灭所产生。当判断出一个单事件是通过正电子湮灭产生的时,用于判断的该单事件的数据确定为即时符合数据。现有判断即时符合数据的方式使得所述即使符合数据中除了包含真符合数据以外,还包含了大量的随机符合数据和散射符合数据,其中,如果由一个正电子湮灭产生的两个γ光子,且在行进过程中没有产生康普顿散射,在探测器上形成的数据,叫做真符合数据。如果若干个正电子在同一时刻或很相近的时刻内发生湮灭,几乎同时产生若干对的γ光子中,很可能有两个不是同一正电子湮灭时产生的γ光子被探测器判定为符合,这样的意外符合事件叫做随机符合,符合随机符合的γ光子在探测器上形成的数据,叫做随机符合数据。如果一对γ光子中的一个或者两个光子在飞行过程中发生康普顿散射,改变了自身的飞行轨迹并降低了自身能量,由于探测器能量分辨率有限,一部分能量改变不大的事件会被误认为真符合事件并记录下来。符合这种情况的γ光子在探测器上形成的数据,叫做散射符合数据。为了去除这些无效数据,需要将其准确的计算出来。
所述调整主要依靠对所述随机符合数据的准确计算,如果随机符合数据的判断存在误差,那么就会对校正结果带来影响,然而,目前没有准确的计算随机符合数据的方法。
发明内容
为了解决上述技术问题,本发明提供了一种数据处理方法和装置,提高确定散射符合数据的准确程度。
本发明实施例公开了如下技术方案:
一种数据处理方法,所述方法包括:
分析弦图的分界,确定出所述弦图中包含真符合数据的区域与不包含所述真符合数据的区域之间的边界,所述弦图为通过透射扫描得到的投影图;
获取所述不包含真符合数据区域中的即时符合数据,通过所述即时符合数据确定出散射分布数据和随机符合数据;
对所述即时符合数据进行处理,具体为:
步骤1:为所述散射分布数据设定第一校正因子,为所述随机符合数据设定第二校正因子,根据所述第一校正因子和所述第二校正因子进行差值计算得到的结果为第一误差值,所述差值计算为所述即时符合数据减去所述散射分布数据和随机符合数据;
步骤2:只调整所述第二校正因子,使得根据调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果最小;
步骤3:根据所述调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果为第二误差值;
步骤4:比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第二校正因子在包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第二校正因子作为所述第二校正因子,重新执行所述步骤1。
优选的,在所述步骤1和所述步骤3之间,还包括步骤2a:
步骤2a:只调整所述第一校正因子,使得根据所述第二校正因子和调整后的第一校正因子进行所述差值计算得到的结果最小;
所述步骤3具体为:根据所述调整后的第一校正因子和所述调整后的第二校正因子进行所述差值计算得到的结果为第二误差值;
所述步骤4具体为:比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第一校正因子作为所述第一校正因子,将调整后的第二校正因子作为所述第二校正因子,重新执行所述步骤1。
优选的,若所述差值处于预设范围内,还包括:
对所述调整后的第二校正因子进行第一校验,以及利用所述调整后的第一校正因子和所述调整后的第二校正因子对所述即时符合数据、所述散射分布数据和所述随机符合数据进行第二校验;
若所述第一校验和所述第二校验均成功,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据。
优选的,所述步骤1,具体为:
err1=Σi|PiNi1RiNi1S′i|;
其中,err1为所述第一误差值,λ1为所述第一校正因子,ω1为所述第二校正因子,Pi为所述即时符合数据,S′i为所述散射分布数据,Ri为所述随机符合数据,Ni为正规化因子,i的取值范围通过所述弦图中包含真符合数据的区域与不包含真符合数据的区域之间的边界确定。
优选的,所述步骤2a具体为:
[ λ 2 ] = arg min λ Σ i | P i N i - ω 1 R i N i - λS i ′ | ;
其中,λ2为所述调整后的第一校正因子。
优选的,所述步骤2具体为:
[ ω 2 ] = arg min ω Σ i | P i N i - ωR i N i - λ 1 S i ′ | ;
其中,ω2为所述调整后的第二校正因子。
优选的,所述步骤3具体为:
err2=Σi|PiNi2RiNi2S′i|;
其中,err2为所述第二误差值。
优选的,
所述第一校验具体为:
判断所述调整后的第二校正因子是否满足第一预设范围;
所述第二校验具体为:
err total = Σ i ( ω 2 R i N i + λ 2 S i ′ ) Σ i P i N i ;
errtotal是否满足第二预设范围。
一种数据处理装置,包括:
确定单元,用于分析弦图的分界,确定出所述弦图中包含真符合数据的区域与不包含所述真符合数据的区域之间的边界,所述弦图为通过透射扫描得到的投影图;
获取单元,用于获取所述不包含真符合数据区域中的即时符合数据,通过所述即时符合数据确定出散射分布数据和随机符合数据;
处理单元,用于对所述即时符合数据进行处理,具体包括:
第一计算子单元,用于为所述散射分布数据设定第一校正因子,为所述随机符合数据设定第二校正因子,根据所述第一校正因子和所述第二校正因子进行差值计算得到的结果为第一误差值,所述差值计算为所述即时符合数据减去所述散射分布数据和随机符合数据;
第一调整子单元,用于只调整所述第二校正因子,使得根据调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果最小;
第二计算子单元,用于根据所述调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果为第二误差值;
比较子单元,用于比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第二校正因子在包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第二校正因子作为所述第二校正因子,重新触发所述第一计算子单元。
优选的,在触发所述第一计算子单元后,触发所述第二计算子单元前,还包括第二调整子单元:
所述第二调整子单元,用于只调整所述第一校正因子,使得根据所述第二校正因子和调整后的第一校正因子进行所述差值计算得到的结果最小;
所述第二计算子单元,还用于根据所述调整后的第一校正因子和所述调整后的第二校正因子进行所述差值计算得到的结果为第二误差值;
所述比较子单元,还用于比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第一校正因子作为所述第一校正因子,将调整后的第二校正因子作为所述第二校正因子,重新触发所述第一计算子单元。
优选的,若所述比较子单元的比较结果为所述差值处于预设范围内,还包括校验子单元:
所述校验子单元,用于对所述调整后的第二校正因子进行第一校验,以及利用所述调整后的第一校正因子和所述调整后的第二校正因子对所述即时符合数据、所述散射分布数据和所述随机符合数据进行第二校验;
所述比较子单元,还用于若所述第一校验和所述第二校验均成功,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据。
由上述技术方案可以看出,通过即时符合数据得到和所述随机符合数据,并为所述随机符合数据设定第二校正因子,通过调整所述第二校正因子以及比较所述第一误差值和所述第二误差值,得到所述调整后的第二校正因子,利用所述调整后的第二校正因子从即时符合数据中确定出所述随机符合数据更加准确。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种数据处理方法的方法流程图;
图2为本发明实施例提供的一种弦图的示意图;
图3为本发明实施例提供的一种数据处理装置的结构图;
图4为本发明实施例提供的一种数据处理装置的结构图;
图5为本发明实施例提供的一种数据处理装置的结构图。
具体实施方式
由于现有技术中缺乏从即时符合数据中准确确定出随机符合数据的手段,导致在利用所述即时符合数据进行数据处理时,以及提高所述即时符合数据中真符合数据的比例时,效果都不理想。也就直接导致了医疗成像装置精度不能进一步提高。为此本发明提供了一种数据处理方法和装置,通过即时符合数据得到和所述随机符合数据,并为所述随机符合数据设定第二校正因子,通过调整所述第二校正因子以及比较所述第一误差值和所述第二误差值,得到所述调整后的第二校正因子,利用所述调整后的第二校正因子从即时符合数据中确定出所述随机符合数据更加准确。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
图1为本发明一种数据处理方法的方法流程图,所述方法包括:
S101:分析弦图的分界,确定出所述弦图中包含真符合数据的区域与不包含所述真符合数据的区域之间的边界,所述弦图为通过透射扫描得到的投影图。
这里需要说明的是,所述弦图是医学造影成像中常用到的投影图,所述弦图的种类不止一个,通常都是被扫描体(比如人体)通过医学成像设备投射扫描后生成的投影图,图2为本发明实施例提供的一种弦图的示意图,如图2所示的就是所述弦图中的一种。其中,b和c为所述弦图中包含真符合数据的区域与不包含所述真符合数据的区域之间的边界,或者可以说,ab之间的区域和cd之间的区域为医学成像设备还未扫描到被扫描体时所得到的投影,bc之间的区域为医学成像设备扫描到被扫描体时所得到的投影。bc之间的区域中具有真符合数据,而一般来说,ab之间和cd之间的区域中不具有真符合数据。
S102:获取所述不包含真符合数据区域中的即时符合数据,通过所述即时符合数据确定出散射分布数据和随机符合数据。
由于ab之间和cd之间的区域中是不具有真符合数据,故这部分区域中的即时符合数据中基本不包含真符合数据,而只包含随机符合数据和散射符合数据。故由于没有真符合数据的干扰,通过处理这部分区域中的所述即时符合数据,可以更为准确的确定出散射符合数据。并可以将对所述不包含真符合数据区域中的即时符合数据的处理结果,引入到对包含真符合数据区域中的即时符合数据的处理上,进一步提高从包含真符合数据区域中的即时符合数据中确定出散射符合数据的准确度。
所述散射分布数据主要是通过所述即时符合数据和所述散射符合数据计算得到,获取的方法可以有多种,本发明对此不一一进行举例。
S103:对所述即时符合数据进行处理,具体为:
步骤1:为所述散射分布数据设定第一校正因子,为所述随机符合数据设定第二校正因子,根据所述第一校正因子和所述第二校正因子进行差值计算得到的结果为第一误差值,所述差值计算为所述即时符合数据减去所述散射分布数据和随机符合数据。
由于通过现有技术根据所述即时符合数据中得到的所述散射分布数据和所述随机符合数据并不准确,故本发明为此为分别为所述散射分布数据和所述随机符合数据设定所述第一校正因子和所述第二校正因子进行误差上的校正。由于S103中的处理的所示即时符合数据不包含所述真符合数据,只包含所述散射分布数据和所述随机符合数据。通过所述差值计算所得到的结果来对第二校正因子进行调整。具体的:
err1=Σi|PiNi1RiNi1S′i|;
其中,err1为所述第一误差值,λ1为所述第一校正因子,ω1为所述第二校正因子,Pi为所述即时符合数据,S′i为所述散射分布数据,Ri为所述随机符合数据,Ni为正规化因子,i的取值范围通过所述弦图中包含真符合数据的区域与不包含真符合数据的区域之间的边界确定。
步骤2:只调整所述第二校正因子,使得根据调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果最小。
也就是说,对所述第一校正因子不进行调整,而只调整所述第二校正因子,确定出当所述误差计算的结果为最小时的调整后的第二校正因子。具体为:
[ ω 2 ] = arg min ω Σ i | P i N i - ωR i N i - λ 1 S i ′ | ;
其中,ω2为所述调整后的第二校正因子。
步骤3:根据所述调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果为第二误差值。
具体为:
err2=Σi|PiNi2RiNi1S′i|;
其中,err2为所述第二误差值。
步骤4:比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第二校正因子在包含真符合数据的即时符合数据中确定随即符合数据。
也就是说比较err1减去err2所得到的差值是否满足预设范围。为了进一步提高判断精度,还可以将得到的所述差值再除以err1,判断计算结果是否处于预定的范围,所述预定的范围可以取值为0到0.0005之间。
若所述差值处于预设范围外,则将所述调整后的第二校正因子作为所述第二校正因子,重新执行所述步骤1。
举例说明,若所述差值处于预设范围外,将ω2作为ω1,将ω2的值作为ω1的值,重新进行步骤1。
由于所述散射符合数据现有技术中也不能准确的确定,故本发明实施例还可以通过既调整所述第一校正因子,也调整所述第二校正因子的方式,提高确定所述散射符合数据的准确度,同时也可以进一步提高确定所述随机符合数据的准确度。为此,在如图1对应的实施例的基础上,本发明实施例还提供了一种数据处理的方法,在所述步骤1和所述步骤3之间,还包括步骤2a:
步骤2a:只调整所述第一校正因子,使得根据所述第二校正因子和调整后的第一校正因子进行所述误差计算得到的结果最小。
也就是说,对所述第一校正因子不进行调整,而只调整所述第二校正因子,确定出当所述第一误差值为最小时的调整后的第二校正因子。具体为:
[ λ 2 ] = arg min λ Σ i | P i N i - ω 1 R i N i - λS i ′ | ;
其中,λ2为所述调整后的第一校正因子。
本发明对所述步骤2和所述步骤2a之间的执行顺序并不进行限定,既可以是所述步骤2和所述步骤2a同时执行,也可以是先后执行。
相应的,所述步骤3针对所述步骤2a也有进一步的相关处理。
所述步骤3具体为:根据所述调整后的第一校正因子和所述调整后的第二校正因子进行所述差值计算得到的结果为第二误差值。
具体为:
err2=Σi|PiNi2RiNi2S′i|;
相应的,所述步骤4针对所述步骤2a也有进一步的相关处理。
所述步骤4具体为:比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据。
若所述差值处于预设范围外,则将所述调整后的第一校正因子作为所述第一校正因子,将调整后的第二校正因子作为所述第二校正因子,重新执行所述步骤1。
比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,即相当于确定出了所述调整后的第一校正因子和所述调整后的第二调整因子。为了对确定出的结果进行进一步的验证,本发明实施例还提供了一种校验方式,至少包括两次校验,具体包括:
对所述调整后的第二校正因子进行第一校验,以及利用所述调整后的第一校正因子和所述调整后的第二校正因子对所述即时符合数据、所述散射分布数据和所述随机符合数据进行第二校验。
这里需要说明的是,所述第一校验主要是校验所述调整后的第二校正因子是否处于合理的范围,当所述调整后的第二校正因子处于第一预设范围内时,则确定所述第一校验成功。其中,举例说明,所述第一预设范围可以设定为1±0.25,也就是0.75到1.25之间。
所述第二校验主要用于校验整体计数值是否正常,举例说明,具体的校验方式可以是:
err total = Σ i ( ω 2 R i N i + λ 2 S i ′ ) Σ i P i N i ;
errtotal是否满足第二预设范围。
举例说明,所述第二预设范围可以设定为1±0.25。
若所述第一校验和所述第二校验均成功,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据。
若所述第一校验和所述第二校验至少一个不成功,则将所述调整后的第一校正因子作为所述第一校正因子,将调整后的第二校正因子作为所述第二校正因子,重新执行所述步骤1。或者,在某些应用场景下,可能会出现一直无法通过所述第一校验和所述第二校验的情况,导致一直重复调整和校验为***带来较大负担,而且这种情况一般可能是医疗成像装置本身故障造成的。故,当所述第一校验和所述第二校验至少一个不成功的次数达到一定次数时,***可以发出警报,提示装置发生故障,进行人工干预。
由上述实施例可以看出,通过即时符合数据得到和所述随机符合数据,并为所述随机符合数据设定第二校正因子,为所述散射分布数据设定第一校正因子,通过调整所述第二校正因子、第一校正因子以及比较所述第一误差值和所述第二误差值,得到所述调整后的第二校正因子和调整后的第一校正因子,利用所述调整后的第二校正因子从即时符合数据中确定出所述随机符合数据更加准确。并通过调整所述第二校正因子的过程来提高调整所述第一校正因子的精度,进一步的提高利用所述调整后的第一校正因子从即时符合数据中确定出所述随机符合数据的精确度。
实施例二
本实施例为对应实施例一的装置实施例,图3为本发明实施例提供的一种数据处理装置的结构图,所述数据处理装置300包括:
确定单元301,用于分析弦图的分界,确定出所述弦图中包含真符合数据的区域与不包含所述真符合数据的区域之间的边界,所述弦图为通过透射扫描得到的投影图。
获取单元302,用于获取所述不包含真符合数据区域中的即时符合数据,通过所述即时符合数据确定出散射分布数据和随机符合数据。
处理单元303,用于对所述即时符合数据进行处理,具体包括:
第一计算子单元3031,用于为所述散射分布数据设定第一校正因子,为所述随机符合数据设定第二校正因子,根据所述第一校正因子和所述第二校正因子进行差值计算得到的结果为第一误差值,所述差值计算为所述即时符合数据减去所述散射分布数据和随机符合数据。
第一调整子单元3032,用于只调整所述第二校正因子,使得根据调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果最小。
第二计算子单元3033,用于根据所述调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果为第二误差值。
比较子单元3034,用于比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第二校正因子在包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第二校正因子作为所述第二校正因子,重新触发所述第一计算子单元3031。
可选的,在图3所对应实施例的的基础上,在触发所述第一计算子单元3031后,触发所述第二计算子单元3033前,还包括第二调整子单元3035,如图4所示,图4为本发明实施例提供的一种数据处理装置的结构图:
所述第二调整子单元3035,用于只调整所述第一校正因子,使得根据所述第二校正因子和调整后的第一校正因子进行所述差值计算得到的结果最小。
所述第二计算子单元3033,还用于根据所述调整后的第一校正因子和所述调整后的第二校正因子进行所述差值计算得到的结果为第二误差值。
所述比较子单元3034,还用于比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第一校正因子作为所述第一校正因子,将调整后的第二校正因子作为所述第二校正因子,重新触发所述第一计算子单元3031。
可选的,在图4所对应实施例的基础上,若所述比较子单元3034的比较结果为所述差值处于预设范围内,还包括校验子单元3036,如图5所示,图5为本发明实施例提供的一种数据处理装置的结构图:
所述校验子单元3036,用于对所述调整后的第二校正因子进行第一校验,以及利用所述调整后的第一校正因子和所述调整后的第二校正因子对所述即时符合数据、所述散射分布数据和所述随机符合数据进行第二校验;
所述比较子单元3034,还用于若所述第一校验和所述第二校验均成功,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据;
若所述第一校验和所述第二校验至少一个不成功,则将所述调整后的第一校正因子作为所述第一校正因子,将调整后的第二校正因子作为所述第二校正因子,重新触发所述第一计算子单元3031。
由上述实施例可以看出,通过即时符合数据得到和所述随机符合数据,并为所述随机符合数据设定第二校正因子,为所述散射分布数据设定第一校正因子,通过调整所述第二校正因子、第一校正因子以及比较所述第一误差值和所述第二误差值,得到所述调整后的第二校正因子和调整后的第一校正因子,利用所述调整后的第二校正因子从即时符合数据中确定出所述随机符合数据更加准确。并通过调整所述第二校正因子的过程来提高调整所述第一校正因子的精度,进一步的提高利用所述调整后的第一校正因子从即时符合数据中确定出所述随机符合数据的精确度。
本发明实施例中提到的第一校正因子、第一误差值、第一校验中、第一计算子单元和第一调整子单元的“第一”只是用来做名字标识,并不代表顺序上的第一。该规则同样适用于“第二”。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到上述实施例方法中的全部或部分步骤可借助软件加通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者诸如媒体网关等网络通信设备)执行本发明各个实施例或者实施例的某些部分所述的方法。
需要说明的是,本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于设备及***实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的设备及***实施例仅仅是示意性的,其中作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述仅是本发明的优选实施方式,并非用于限定本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (11)

1.一种数据处理方法,其特征在于,所述方法包括:
分析弦图的分界,确定出所述弦图中包含真符合数据的区域与不包含所述真符合数据的区域之间的边界,所述弦图为通过透射扫描得到的投影图;
获取所述不包含真符合数据区域中的即时符合数据,通过所述即时符合数据确定出散射分布数据和随机符合数据;
对所述即时符合数据进行处理,具体为:
步骤1:为所述散射分布数据设定第一校正因子,为所述随机符合数据设定第二校正因子,根据所述第一校正因子和所述第二校正因子进行差值计算得到的结果为第一误差值,所述差值计算为所述即时符合数据减去所述散射分布数据和随机符合数据;
步骤2:只调整所述第二校正因子,使得根据调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果最小;
步骤3:根据所述调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果为第二误差值;
步骤4:比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第二校正因子在包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第二校正因子作为所述第二校正因子,重新执行所述步骤1。
2.根据权利要求1所述的方法,其特征在于,在所述步骤1和所述步骤3之间,还包括步骤2a:
步骤2a:只调整所述第一校正因子,使得根据所述第二校正因子和调整后的第一校正因子进行所述差值计算得到的结果最小;
所述步骤3具体为:根据所述调整后的第一校正因子和所述调整后的第二校正因子进行所述差值计算得到的结果为第二误差值;
所述步骤4具体为:比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第一校正因子作为所述第一校正因子,将调整后的第二校正因子作为所述第二校正因子,重新执行所述步骤1。
3.根据权利要求2所述的方法,其特征在于,若所述差值处于预设范围内,还包括:
对所述调整后的第二校正因子进行第一校验,以及利用所述调整后的第一校正因子和所述调整后的第二校正因子对所述即时符合数据、所述散射分布数据和所述随机符合数据进行第二校验;
若所述第一校验和所述第二校验均成功,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据。
4.根据权利要求1至3任一项所述的方法,其特征在于,所述步骤1,具体为:
err 1 = Σ i | P i N i - ω 1 R i N i - λ 1 S i ′ | ;
其中,err1为所述第一误差值,λ1为所述第一校正因子,ω1为所述第二校正因子,Pi为所述即时符合数据,S′i为所述散射分布数据,Ri为所述随机符合数据,Ni为正规化因子,i的取值范围通过所述弦图中包含真符合数据的区域与不包含真符合数据的区域之间的边界确定。
5.根据权利要求2所述的方法,其特征在于,所述步骤2a具体为:
[ λ 2 ] = arg min λ Σ i | P i N i - ω 1 R i N i - λ S i ′ | ;
其中,λ2为所述调整后的第一校正因子。
6.根据权利要求1至3任一项所述的方法,其特征在于,所述步骤2具体为:
[ ω 2 ] = arg min ω Σ i | P i N i - ω R i N i - λ 1 S i ′ | ;
其中,ω2为所述调整后的第二校正因子。
7.根据权利要求2或3所述的方法,其特征在于,所述步骤3具体为:
err 2 = Σ i | P i N i - ω 2 R i N i - λ 2 S i ′ | ;
其中,err2为所述第二误差值。
8.根据权利要求3所述的方法,其特征在于,
所述第一校验具体为:
判断所述调整后的第二校正因子是否满足第一预设范围;
所述第二校验具体为:
err total = Σ i ( ω 2 R i N i + λ 2 S i ′ ) Σ i P i N i ;
errtotal是否满足第二预设范围。
9.一种数据处理装置,其特征在于,包括:
确定单元,用于分析弦图的分界,确定出所述弦图中包含真符合数据的区域与不包含所述真符合数据的区域之间的边界,所述弦图为通过透射扫描得到的投影图;
获取单元,用于获取所述不包含真符合数据区域中的即时符合数据,通过所述即时符合数据确定出散射分布数据和随机符合数据;
处理单元,用于对所述即时符合数据进行处理,具体包括:
第一计算子单元,用于为所述散射分布数据设定第一校正因子,为所述随机符合数据设定第二校正因子,根据所述第一校正因子和所述第二校正因子进行差值计算得到的结果为第一误差值,所述差值计算为所述即时符合数据减去所述散射分布数据和随机符合数据;
第一调整子单元,用于只调整所述第二校正因子,使得根据调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果最小;
第二计算子单元,用于根据所述调整后的第二校正因子和所述第一校正因子进行所述差值计算得到的结果为第二误差值;
比较子单元,用于比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第二校正因子在包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第二校正因子作为所述第二校正因子,重新触发所述第一计算子单元。
10.根据权利要求9所述的装置,其特征在于,在触发所述第一计算子单元后,触发所述第二计算子单元前,还包括第二调整子单元:
所述第二调整子单元,用于只调整所述第一校正因子,使得根据所述第二校正因子和调整后的第一校正因子进行所述差值计算得到的结果最小;
所述第二计算子单元,还用于根据所述调整后的第一校正因子和所述调整后的第二校正因子进行所述差值计算得到的结果为第二误差值;
所述比较子单元,还用于比较所述第一误差值和所述第二误差值之间的差值,若所述差值处于预设范围内,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据;
若所述差值处于预设范围外,则将所述调整后的第一校正因子作为所述第一校正因子,将调整后的第二校正因子作为所述第二校正因子,重新触发所述第一计算子单元。
11.根据权利要求10所述的装置,其特征在于,若所述比较子单元的比较结果为所述差值处于预设范围内,还包括校验子单元:
所述校验子单元,用于对所述调整后的第二校正因子进行第一校验,以及利用所述调整后的第一校正因子和所述调整后的第二校正因子对所述即时符合数据、所述散射分布数据和所述随机符合数据进行第二校验;
所述比较子单元,还用于若所述第一校验和所述第二校验均成功,则利用所述调整后的第一校正因子在所述包含真符合数据的即时符合数据中确定散射分布数据以及利用所述调整后的第二校正因子在所述包含真符合数据的即时符合数据中确定随机符合数据。
CN201410562261.2A 2014-10-21 2014-10-21 一种数据处理方法和装置 Active CN104352244B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410562261.2A CN104352244B (zh) 2014-10-21 2014-10-21 一种数据处理方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410562261.2A CN104352244B (zh) 2014-10-21 2014-10-21 一种数据处理方法和装置

Publications (2)

Publication Number Publication Date
CN104352244A true CN104352244A (zh) 2015-02-18
CN104352244B CN104352244B (zh) 2016-06-29

Family

ID=52519629

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410562261.2A Active CN104352244B (zh) 2014-10-21 2014-10-21 一种数据处理方法和装置

Country Status (1)

Country Link
CN (1) CN104352244B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017084530A1 (en) * 2015-11-19 2017-05-26 Shanghai United Imaging Healthcare Co., Ltd. Detector in an imaging system
US9715023B2 (en) 2015-11-19 2017-07-25 Shanghai United Imaging Healthcare Co., Ltd. Detector in an imaging system
CN107464270A (zh) * 2017-07-17 2017-12-12 沈阳东软医疗***有限公司 一种图像重建方法和装置
CN109289129A (zh) * 2018-09-07 2019-02-01 广州医科大学附属肿瘤医院 直线加速器的质控方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060151705A1 (en) * 2005-01-12 2006-07-13 General Electric Company Method and system for scatter correction in a positron emission tomography system
US20060163485A1 (en) * 2005-01-21 2006-07-27 General Electric Company Method and system for scattered coincidence estimation in a time-of-flight positron emission tomography system
US20070040122A1 (en) * 2005-08-17 2007-02-22 General Electric Company Technique for reconstructing PET scan images
JP2011220719A (ja) * 2010-04-05 2011-11-04 Toshiba Corp 核医学診断装置、医用画像処理装置および医用画像診断装置
US20120070050A1 (en) * 2010-09-20 2012-03-22 Siemens Medical Solutions Usa, Inc. Time of Flight Scatter Distribution Estimation in Positron Emission Tomography
CN103186882A (zh) * 2011-12-30 2013-07-03 沈阳东软派斯通医疗***有限公司 Pet***中图像的衰减校正方法及装置
CN103417235A (zh) * 2013-07-31 2013-12-04 沈阳东软医疗***有限公司 一种随机噪声的校正方法和装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060151705A1 (en) * 2005-01-12 2006-07-13 General Electric Company Method and system for scatter correction in a positron emission tomography system
US20060163485A1 (en) * 2005-01-21 2006-07-27 General Electric Company Method and system for scattered coincidence estimation in a time-of-flight positron emission tomography system
US20070040122A1 (en) * 2005-08-17 2007-02-22 General Electric Company Technique for reconstructing PET scan images
JP2011220719A (ja) * 2010-04-05 2011-11-04 Toshiba Corp 核医学診断装置、医用画像処理装置および医用画像診断装置
US20120070050A1 (en) * 2010-09-20 2012-03-22 Siemens Medical Solutions Usa, Inc. Time of Flight Scatter Distribution Estimation in Positron Emission Tomography
CN103186882A (zh) * 2011-12-30 2013-07-03 沈阳东软派斯通医疗***有限公司 Pet***中图像的衰减校正方法及装置
CN103417235A (zh) * 2013-07-31 2013-12-04 沈阳东软医疗***有限公司 一种随机噪声的校正方法和装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JOSEP F OLIVER,ET AL.,: "Improving the singles rate method for modeling accidental coincidences in high-resolution PET", 《PHYSICS IN MEDICINE AND BIOLOGY》, vol. 55, 31 December 2010 (2010-12-31) *
吴朝霞 等: "基于非线性统计估计技术的PET散射校正", 《生物医学工程学杂志》, vol. 22, no. 1, 31 December 2005 (2005-12-31) *
陆汉魁: "PET显像的散射校正和衰减校正", 《中华核医学杂志》, vol. 24, no. 1, 28 February 2014 (2014-02-28) *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017084530A1 (en) * 2015-11-19 2017-05-26 Shanghai United Imaging Healthcare Co., Ltd. Detector in an imaging system
US9715023B2 (en) 2015-11-19 2017-07-25 Shanghai United Imaging Healthcare Co., Ltd. Detector in an imaging system
US10338239B2 (en) 2015-11-19 2019-07-02 Shanghai United Imaging Healthcare Co., Ltd. Detector in an imaging system
US10634801B2 (en) 2015-11-19 2020-04-28 Shanghai United Imaging Healthcare Co., Ltd. Detector in an imaging system
US11644586B2 (en) 2015-11-19 2023-05-09 Shanghai United Imaging Healthcare Co., Ltd. Method and system for evaluating the working condition of a detector
CN107464270A (zh) * 2017-07-17 2017-12-12 沈阳东软医疗***有限公司 一种图像重建方法和装置
CN107464270B (zh) * 2017-07-17 2020-08-11 东软医疗***股份有限公司 一种图像重建方法和装置
CN109289129A (zh) * 2018-09-07 2019-02-01 广州医科大学附属肿瘤医院 直线加速器的质控方法及装置
CN109289129B (zh) * 2018-09-07 2020-12-15 广州医科大学附属肿瘤医院 直线加速器的质控方法及装置

Also Published As

Publication number Publication date
CN104352244B (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
RU2524053C2 (ru) Способ непрерывного моделирования времяпролетного рассеяния
Agnew et al. Monitoring daily MLC positional errors using trajectory log files and EPID measurements for IMRT and VMAT deliveries
US20180185672A1 (en) Determination of radiation collimator component position
US11791041B2 (en) Acceptance, commissioning, and ongoing benchmarking of a linear accelerator (LINAC) using an electronic portal imaging device (EPID)
US10319118B2 (en) Reconstructing PET image
Kozelka et al. Optimizing the accuracy of a helical diode array dosimeter: a comprehensive calibration methodology coupled with a novel virtual inclinometer
CN104335247A (zh) 在pet重建中的快速散射估计
CN104352244B (zh) 一种数据处理方法和装置
US8089043B2 (en) Systems and methods for calibrating time alignment for a positron emission tomography (PET) system
RU2535635C2 (ru) Динамическое формирование изображений посредством позитронно-эмиссионной томографии с компенсацией загрязнений изотопами
CN105025796A (zh) X射线计算机断层摄影装置以及校正方法
CN105431884B (zh) 用于对图像中的级联伽马进行解释和建模的方法
US20080135764A1 (en) Device and method for locally resolved control of a radiation dose
US10036817B2 (en) Solving outside-field of view scatter correction problem in positron emission tomography via digital experimentation
CN110025329B (zh) 符合计数弦图生成方法、装置、计算机设备和存储介质
CN105496436A (zh) 用于pet装置的时间校正方法和装置
US11282242B2 (en) Scatter correction for positron emission tomography (PET)
US10231688B2 (en) Radiation image analysis device, method, and program
Fuangrod et al. An independent system for real-time dynamic multileaf collimation trajectory verification using EPID
Zwan et al. Commissioning and quality assurance for VMAT delivery systems: an efficient time‐resolved system using real‐time EPID imaging
Park et al. Gamma-index method sensitivity for gauging plan delivery accuracy of volumetric modulated arc therapy
EP2950119B1 (en) System and method for verifying a particle beam
US20190167213A1 (en) Systems and methods for imaging with improved dosages
Chojnowski et al. An automatic method of the isocentre position verification for micromultileaf collimator based radiosurgery system
CN108434612B (zh) 一种质量保证测试方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20170307

Address after: Room 6, building 3, building No. 8, Wang Xi Road, Haidian District, Beijing, 100193, China

Patentee after: Beijing Neusoft Medical Equipment Co.,Ltd.

Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16

Patentee before: SHENYANG NEUSOFT MEDICAL SYSTEMS Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220321

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

Patentee after: Shenyang Neusoft Medical Systems Co.,Ltd.

Address before: Room 322, Room 3, Building 6, 8 Wangxi Road, Haidian District, Beijing, 100193

Patentee before: Beijing Neusoft Medical Equipment Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230413

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

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

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

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

TR01 Transfer of patent right