CN110755100B - 校正方法、装置、控制台设备及pet系统 - Google Patents

校正方法、装置、控制台设备及pet系统 Download PDF

Info

Publication number
CN110755100B
CN110755100B CN201910989334.9A CN201910989334A CN110755100B CN 110755100 B CN110755100 B CN 110755100B CN 201910989334 A CN201910989334 A CN 201910989334A CN 110755100 B CN110755100 B CN 110755100B
Authority
CN
China
Prior art keywords
bed
data
projection data
scanning
double
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
Application number
CN201910989334.9A
Other languages
English (en)
Other versions
CN110755100A (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
Shenyang Zhihe Medical Technology 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 Shenyang Zhihe Medical Technology Co ltd filed Critical Shenyang Zhihe Medical Technology Co ltd
Priority to CN201910989334.9A priority Critical patent/CN110755100B/zh
Publication of CN110755100A publication Critical patent/CN110755100A/zh
Application granted granted Critical
Publication of CN110755100B publication Critical patent/CN110755100B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data

Abstract

本发明实施例提供一种校正方法、装置、控制台设备及PET系统。本发明实施例通过获得对模体进行PET扫描的扫描数据,根据扫描数据,获得模拟双床位投影数据,模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数,根据模拟双床位投影数据进行图像重建,得到目标重建图像,基于目标重建图像,确定目标重叠层数对应的目标校正因子,利用目标校正因子对待校正重建图像中各个像素值进行校正,获得校正后的重建图像,不需要对每个重叠层数都进行扫描,减少了扫描次数,缩短了整个校正过程的用时,提高了校正的处理效率。

Description

校正方法、装置、控制台设备及PET系统
技术领域
本发明涉及医学图像处理技术领域,尤其涉及一种校正方法、装置、控制台设备及PET系统。
背景技术
正电子发射断层成像(Positron Emission Tomography,PET)技术是核医学尖端技术的代表,它利用放射性核素示踪剂对组织功能进行显像,能够非介入的测量特定目标结构或感兴趣区域放射性药物的生物分布随时间量化变化。在疾病诊断、病理研究、药物研发等多个方面具有广泛的应用。
利用PET扫描采集到的目标对象的数据在进行图像重建后,需要利用校正因子进行校正,以获得符合要求的重建图像。相关技术中,获取校正因子的方式是:按照不同重叠层数,对注药后的水模进行多次扫描,每个重叠层数对应一次扫描;每次扫描需要扫描两个床位,非常耗时,因此处理效率较低。
发明内容
为克服相关技术中存在的问题,本发明提供了一种校正方法、装置、控制台设备及PET系统,提高校正的处理效率。
根据本发明实施例的第一方面,提供一种校正方法,用于对正电子发射断层成像PET系统的待校正重建图像进行校正,所述方法包括:
获得对模体进行PET扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子;
利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
根据本发明实施例的第二方面,提供一种校正装置,用于对正电子发射断层成像PET系统的待校正重建图像进行校正,所述装置包括:
数据获取模块,用于获得对模体进行PET扫描的扫描数据;
模拟模块,用于根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
重建模块,用于根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
确定模块,用于基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子;
校正模块,用于利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
根据本发明实施例的第三方面,提供一种控制台设备,包括内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于连接PET系统的探测器,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;
所述存储器,用于存储校正逻辑对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行如下操作:
获得对模体进行PET扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子;
利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
根据本发明实施例的第四方面,提供一种PET系统,包括探测器、扫描床和控制台设备,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;其中,
所述闪烁晶体,用于在所述PET系统扫描过程中,探测被检体内发出的高能光子,并将所述高能光子转换为光信号;
所述光电转换器件,用于将所述光信号转换成电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述控制台设备,用于获得对模体进行PET扫描的扫描数据;根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子;利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
本发明实施例提供的技术方案可以包括以下有益效果:
本发明实施例,通过获得对模体进行PET扫描的扫描数据,根据扫描数据,获得模拟双床位投影数据,模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数,根据模拟双床位投影数据进行图像重建,得到目标重建图像,基于目标重建图像,确定目标重叠层数对应的目标校正因子,利用目标校正因子对待校正重建图像中各个像素值进行校正,获得校正后的重建图像,不需要对每个重叠层数都进行扫描,减少了扫描次数,缩短了整个校正过程的用时,提高了校正的处理效率。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本说明书。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本说明书的实施例,并与说明书一起用于解释本说明书的原理。
图1是对水模进行单床位扫描的示意图。
图2是对水模进行双床位扫描的示意图。
图3是本发明实施例提供的校正方法的流程示例图。
图4是本发明实施例提供的校正装置的功能方块图。
图5是本发明实施例提供的控制台设备的一个硬件结构图。
图6是本发明实施例提供的获取模拟双床位的投影数据的第一示例图。
图7是本发明实施例提供的获取模拟双床位的投影数据的第二示例图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明实施例的一些方面相一致的装置和方法的例子。
在本发明实施例使用的术语是仅仅出于描述特定本发明实施例的目的,而非旨在限制本发明实施例。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本发明实施例可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本发明实施例范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
本文中的扫描均指PET扫描。在对PET重建图像的校正中,通常需要用到两种校正因子,一种是活性浓度校正因子(Activity Concentration Correction Factor,ACCF),另一种是层灵敏度校正因子。
活性浓度校正因子用于活性浓度校正(Activity Concentration Correction,缩写为ACC),也称井型计数校正(Well Counter Correction,缩写为WCC)或者定量校正,目的是利用活性浓度校正因子将重建图像中的像素值由计数率转化为绝对活度水平,以便得到供医师使用的有价值信息,例如用于进行病灶鉴别、疗效评价、半定量分析等的信息。层灵敏度校正因子用于对重建图像的像素值进行校正,可以提高重建图像的图像质量。
本发明实施例提供的校正方法可以用于对PET系统的待校正重建图像进行校正。其中的校正因子可以包括活性浓度校正因子、层灵敏度校正因子中的至少一种。
首先,说明一下PET的单床位扫描和双床位扫描的概念。PET扫描时,被扫描的目标对象被置于扫描床上,目标对象与扫描床的相对位置保持不变。单床位扫描是指将扫描床固定在一个位置进行扫描;双床位扫描是指将扫描床固定在两个位置,对每个位置分别进行扫描,即将扫描床固定在第一位置扫描一次,然后将扫描床固定在第二位置再扫描一次。
以对水模进行PET扫描为例。图1是对水模进行单床位扫描的示意图,图中灰色区域即为单个床位的位置。图2是对水模进行双床位扫描的示意图,图中灰色区域为第一床位的位置,虚线框为第二床位的位置,箭头所示方向为床位的移动方向。图1和图2中,位于PET/CT探测设备之间的部分为扫描区域。相应地,本发明实施例提供的方法不仅可以用于PET,还可以用于CT。
其次,说明一下本实施例的理论依据。以模体为长水模为例。在获取校正所需的校正因子时,需要对长水模进行PET扫描。在对注药后的长水模进行PET扫描时,水模药剂分布均匀,并且采集扫描数据的方式为步进式,即每次扫描时,扫描床被固定在指定床位处,在扫描过程中扫描床位置保持不变。这样,排除掉探测器外部粒子的影响,在扫描均匀长水模不同的区域,单个床位扫描的计数分布情况可以看作是一致的,即可以认为扫描长水模的第一区域时的计数(即发生符合事件的次数)分布情况与扫描长水模的第二区域时的计数分布情况是一致的,从而可以根据长水模第一区域对应的扫描数据,模拟出长水模第二区域对应的扫描数据,这样就不需要通过对第二区域进行PET扫描来获得第二区域的扫描数据了。第一区域和第二区域均被水模覆盖,且第一区域的长度和第二区域的长度都是由PET/CT探测设备扫描区域确定的,二者长度相等。
下面通过实施例对本发明提供的校正方法的流程进行说明。
图3是本发明实施例提供的校正方法的流程示例图。该校正方法可以应用于PET系统,用于对PET系统的待校正重建图像进行校正。如图3所示,校正方法可以包括:
S301,获得对模体进行PET扫描的扫描数据。
S302,根据扫描数据,获得模拟双床位投影数据,模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数。
S303,根据模拟双床位投影数据进行图像重建,得到目标重建图像。
S304,基于目标重建图像,确定目标重叠层数对应的目标校正因子。
S305,利用目标校正因子对待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
在步骤S301中,模体的长度大于设定长度,并且模体中包含放射性溶液。
在步骤S301中,PET扫描时的床位可以是单床位或双床位。因此,该步骤中,只需要对模体进行一次单床位扫描或者一次双床位扫描。不需要相关技术那样,对每个重叠层数,都要进行一次扫描,减少了扫描次数,从而节省了时间,缩短了整个校正过程的用时,因此能够提高处理效率。
例如,假设需要获取重叠层数分别为11层、13层、15层、17层、19层、21层时的校正因子。则只需要对模体,例如长水模,进行一次单床位扫描或者一次双床位扫描,而不需要分别在重叠层数为11层、13层、15层、17层、19层、21层的情况下各扫描一次。
以长水模为例。长水模需要在注射放射性溶液后进行扫描。由于水模活度较低且放射性溶液一直在衰变,一次注药不足以完成相关技术中全部重叠层数的扫描,因此相关技术中随着扫描次数的增多需要多次注药,对放射性溶液的损耗较多,造成成本升高。本实施例只需一次扫描,因此仅需一次注药,减少了放射性溶液使用,节约了成本。
并且,由于扫描期间需要工作人员一直监测,并且每次扫描开始前还需要工作人员根据需求手动修改相关参数,必要时还要对模体重新注药,在扫描次数多时人工操作量较大。本实施例由于只进行一次扫描,因此工作人员所需执行的操作仅为一次注药和一次扫描,显著减少了工作人员的操作量。
在步骤S302中,模拟双床位投影数据是利用已获得的扫描数据模拟出来的,不是通过扫描采集到的,相比于扫描获得投影数据的方式用时大大减少,因此可以有效缩短整个校正过程的时间,提高处理效率。
模拟双床位投影数据的组成可以分为两部分,一部分是两个床位的重叠区域的投影数据,另一部分是两个床位的非重叠区域的投影数据。其中,非重叠区域的投影数据包括两个床位的非重叠部分各自对应区域的投影数据。
例如,假设目标双床位为第一床位和第二床位,则模拟双床位投影数据包括:第一床位和第二床位的重叠区域的投影数据、非重叠部分中属于第一床位的投影数据、以及非重叠部分中属于第二床位的投影数据。其中,重叠区域的重叠层数可以是任意数目。
在步骤S303中,图像重建可以采用迭代法进行,例如,可以采用有序子集最大期望值(Ordered Subset Expectation Maximization,OSEM)法进行图像重建。
在步骤S304中,目标校正因子可以包括双床位活性浓度校正因子、层灵敏度校正因子中的至少一个。
基于重建图像确定层灵敏度校正因子的方法可包括:在每层水模重建图像上选取合适大小和形状的感兴趣区域(Region of Interest,ROI)),例如,选取的感兴趣区域可以是与实际水模半径接近的圆形ROI,计算出各层ROI内的总计数值/率,之后利用如下的公式计算各层对应的层灵敏度校正因子(Slice Sensitivity Correction Factor,SSCF)
Figure BDA0002237734700000081
上述公式中,j=1,2...N,N是重建像的层数,Cj为第j层重建图像ROI内的总计数值/率。
基于重建图像确定活性浓度校正因子的方法可包括:基于上一步层灵敏度校正后的水模重建像,也可以基于不做层灵敏度校正的重建图像。在每层水模重建图像上选取合适大小和形状的ROI,例如,选取的感兴趣区域可以是与实际水模半径接近的圆形ROI,利用如下的活性浓度校正因子计算公式即可得到ACCF:
Figure BDA0002237734700000082
其中,A表示水模内放射性溶液的已知活性浓度,B.F.为分支比,R表示每秒单位体像素的计数值。
在步骤S305中,利用目标校正因子对待校正重建图像中各个像素值进行校正的过程可包括:利用活性浓度校正因子ACCF对待校正重建图像进行乘法校正,公式为:
Figure BDA0002237734700000083
其中,A即为活性浓度校正后的重建图像,各像素值表示活性浓度;R表示待校正重建图像单位体像素的平均计数率,B.F.为分支比。
在一个示例性的实现过程中,在步骤S301中,获得对模体进行PET扫描的扫描数据,可以包括:
获得对模体进行单个床位扫描的扫描数据;或者
获得对模体进行双床位扫描的扫描数据。
本实施例中,当扫描数据为单床位扫描数据(本文中,单床位扫描数据指前述的单个床位扫描的扫描数据)时,只需对长水模进行一次扫描,且在该次扫描中仅需对单个床位进行扫描,用时短,因此处理效率高;当扫描数据为双床位扫描数据(本文中,双床位扫描数据指前述的双床位扫描的扫描数据)时,只需对长水模进行两次扫描,即分别对两个床位进行扫描,用时相对于单床位扫描增加一倍,但相对于相关技术仍然缩短很多,并且利用双床位扫描数据得到的校正因子准确度更高。在实际应用中,可以根据实际需求选择采集单床位扫描数据或者双床位扫描数据。
在一个示例性的实现过程中,扫描数据为单床位扫描数据;
步骤S302中,根据扫描数据,获得模拟双床位投影数据,可以包括:
根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将扫描数据作为双床位中第一床位的投影数据,并将扫描数据作为双床位中第二床位的投影数据;
基于环偏移量确定第一床位与第二床位的重叠区域和非重叠区域;
确定重叠区域对应的第一投影数据等于第一床位与第二床位中相应投影数据的和,确定非重叠区域对应的第二投影数据等于相应床位的投影数据,第一投影数据和第二投影数据组成模拟双床位投影数据。
假设目标重叠层数为combineSliceNum,探测器总环数为ringNum,则可以通过如下的关系式确定两个床位的环偏移量ringOffSet:重叠环数combineRingNum=(目标重叠层数combineSliceNum+1)/2;环偏移量ringOffSet=探测器总环数ringNum–重叠环数combineRingNum;探测器是由多个探测器环在轴向(平行于扫描床的方向)方向上排列形成的,探测器总环数即是指构成探测器的所有探测器环的数量。
其中,基于环偏移量确定第一床位与第二床位的重叠区域和非重叠区域,即将第二床位中各响应线(Line of Response,LOR)对应的两个探测环索引依据环偏移量进行移动,由此也可确定出第二床位与第一床位的重叠区域,其他未重叠部分为非重叠区域。
重叠区域对应的第一投影数据等于第一床位与第二床位中相应投影数据的和,计算过程为:
combinedPrjData(lorCom)=Bed1_PrjData(lorCom)+Bed2_PrjData(lorCom)
其中,lorCom表示重叠区域LOR索引,combinedPrjData(lorCom)为重叠区域的投影数据,Bed1_PrjData(lorCom)和Bed2_PrjData(lorCom)分别表示第一床和第二床中对应响应线lorCom的投影数据,其他非重叠区域数据保持不变。重叠区域与非重叠区域共同组成了模拟双床位投影数据。
在一个示例性的实现过程中,扫描数据为单床位扫描数据;步骤S302中,根据扫描数据,获得模拟双床位投影数据,可以包括:
将扫描数据作为双床位中第一床位的投影数据,并将扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Y的数据量,并将模拟双床位投影数据Y中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
举例说明。图6是本发明实施例提供的获取模拟双床位的投影数据的第一示例图。
请参见图6,在图6(a)中,将单床位扫描数据X分别赋值给模拟的第一床位X-1和第二床位X-2,第一床位X-1和第二床位X-2连续。
在图6(b)中,根据目标重叠层数,确定第一床位X-1和第二床位X-2重叠合并后新数据Y(即模拟双床位投影数据Y)的大小,并将Y中各位置的值初始化为-1。
在图6(c)中,按照轴向(图6中的z方向)位置,将第一床位X-1的数据复制到第一床位X-1在Y中的对应位置(图6(c)中的左侧区域为第一床位X-1所在的第一区域),Y的剩余位置仍保持初始值-1不变。
在图6(d)中,对Y中的第二个床位进行处理。计算环偏移量:重叠环数=(目标的重叠层数+1)/2;环偏移量=探测器总环数–重叠环数;根据环偏移量可确定第二床位在Y中的起始位置(即图6(d)中的第二床起始位置),第二床位X-2在Y中的第二区域为第二床起始位置至Y的最右侧边界之间的区域。如果第二区域中某位置的值小于0,则直接将第二床位X-2的数据复制到该位置中;如果第二区域中某位置的值大于或等于0,则将第二床位X-2的值与Y在该位置的当前值进行求和,将和的值赋值到该位置,这样,就得到了Y中重叠区域的投影数据。其中,求和的区域即为模拟双床位的重叠区域(如图6(d)中的中间区域),其余部分为非重叠区域。
本实施例中,根据采集的单床位扫描数据模拟出模体的双床位投影数据(即模拟双床位投影数据),不需要针对每个重叠层数都进行扫描,减少了扫描次数,节省了时间,提高了处理效率。
在一个示例性的实现过程中,扫描数据为双床位扫描数据;
步骤S302中,根据扫描数据,获得模拟双床位投影数据,可以包括:
根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
基于环偏移量确定第一床位与第二床位的重叠区域和非重叠区域;
确定重叠区域对应的第一投影数据等于第一床位与第二床位中相应投影数据的和,确定非重叠区域对应的第二投影数据等于相应床位的投影数据,第一投影数据和第二投影数据组成模拟双床位投影数据。
由双床位扫描数据获得模拟双床位投影数据的过程中,重叠环数combineRingNum、环偏移量ringOffSet、重叠区域对应的第一投影数据的计算方式的计算过程与上述由单床位扫描数据获得模拟双床位投影数据的过程中的相同,此处不再赘述。
在一个示例性的实现过程中,扫描数据为双床位扫描数据;步骤S302中,根据扫描数据,获得模拟双床位投影数据,可以包括:
将扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Z的数据量,并将模拟双床位投影数据Z中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
举例说明。图7是本发明实施例提供的获取模拟双床位的投影数据的第二示例图。
请参见图7,在图7(a)中,将两个连续床位的扫描数据X和Y分别赋值给模拟的第一床位和第二床位,假设该两个床位重叠层数为15。
在图7(b)中,假设目标重叠层数11,根据目标重叠层数11确定第一床位和第二床位重叠合并后新数据Z(即模拟双床位投影数据Z)的大小,并将Z中各位置的值初始化为-1。
在图7(c)中,按照轴向(图7中的z方向)位置,将第一床位的数据X复制到第一床位在Z中的对应位置(图7(c)中的左侧区域为第一床位所在的第一区域),Z的剩余位置仍保持初始值-1不变。
在图7(d)中,对Z中的第二个床位进行处理。计算环偏移量:重叠环数=(目标的重叠层数+1)/2;环偏移量=探测器总环数–重叠环数;根据环偏移量可确定第二床位在Z中的起始位置(即图7(d)中的第二床起始位置),第二床位在Z中的第二区域为第二床起始位置至Z的最右侧边界之间的区域。如果第二区域中某位置的值小于0,则直接将第二床位的数据复制到该位置中;如果第二区域中某位置的值大于或等于0,则将第二床位的值与Z在该位置的当前值进行求和,将和的值赋值到该位置,这样,就得到了Z中重叠区域的投影数据。其中,求和的区域即为模拟双床位的重叠区域(如图7(d)中的中间区域),其余部分为非重叠区域。
本实施例中,根据采集的双床位扫描数据模拟出模体的模拟双床位投影数据,不需要针对每个重叠层数都进行扫描,节省了时间,提高了处理效率。并且,本实施例利用双床位扫描数据模拟出的双床位投影数据准确度更高,从而可以获得更准确的校正因子。
本发明实施例提供的校正方法,通过获得对模体进行PET扫描的扫描数据,根据扫描数据,获得模拟双床位投影数据,模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数,根据模拟双床位投影数据进行图像重建,得到目标重建图像,基于目标重建图像,确定目标重叠层数对应的目标校正因子,利用目标校正因子对待校正重建图像中各个像素值进行校正,获得校正后的重建图像,不需要对每个重叠层数都进行扫描,减少了扫描次数,缩短了整个校正过程的用时,提高了校正的处理效率。
通过本实施例获得目标校正因子后,可以存储目标重叠层数和目标校正因子的对应关系。在扫描应用中,当通过PET扫描获得待检测对象的扫描数据、并根据扫描数据得到待校正重建图像后,可以根据实际扫描时的重叠层数,找到对应的目标重叠层数,从而根据存储的上述对应关系确定目标校正因子,直接将目标校正因子作用于待校正重建图像进行校正。
基于上述的方法实施例,本发明实施例还提供了相应的装置、设备、系统及存储介质实施例。关于本发明实施例的装置、系统及存储介质实施例的详细实现方式,请参见前述方法实施例部分的相应说明。
图4是本发明实施例提供的校正装置的功能方块图。如图4所示,本实施例中,校正装置用于对正电子发射断层成像PET系统的待校正重建图像进行校正,校正装置可以包括:
数据获取模块410,用于获得对模体进行PET扫描的扫描数据;
模拟模块420,用于根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
重建模块430,用于根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
确定模块440,用于基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子;
校正模块450,用于利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
在一个示例性的实现过程中,数据获取模块410可以具体用于:
获得对模体进行单个床位扫描的扫描数据;或者
获得对模体进行双床位扫描的扫描数据。
在一个示例性的实现过程中,所述扫描数据为单个床位扫描的扫描数据;
模拟模块420可以具体用于:
根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据作为双床位中第一床位的投影数据,将所述扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为双床位扫描的扫描数据;
模拟模块420可以具体用于:
根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将所述扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为单个床位扫描的扫描数据;
模拟模块420可以具体用于:
将扫描数据作为双床位中第一床位的投影数据,并将扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Y的数据量,并将模拟双床位投影数据Y中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为双床位扫描数据;
模拟模块420可以具体用于:
将扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Z的数据量,并将模拟双床位投影数据Z中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
在一个示例性的实现过程中,所述目标校正因子包括双床位活性浓度校正因子、层灵敏度校正因子中的至少一个。
本发明实施例还提供了一种控制台设备。图5是本发明实施例提供的控制台设备的一个硬件结构图。如图5所示,PET系统包括:内部总线501,以及通过内部总线连接的存储器502,处理器503和外部接口504。
所述外部接口504,用于连接PET系统的探测器,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;
所述存储器502,用于存储校正逻辑对应的机器可读指令;所述处理器503,用于读取存储器502上的机器可读指令,并执行所述指令以实现如下操作:
获得对模体进行PET扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子;
利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
在一个示例性的实现过程中,获得对模体进行PET扫描的扫描数据,包括:
获得对模体进行单个床位扫描的扫描数据;或者
获得对模体进行双床位扫描的扫描数据。
在一个示例性的实现过程中,所述扫描数据为单个床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:
根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据作为双床位中第一床位的投影数据,将所述扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为双床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:
根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将所述扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为单个床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:
将扫描数据作为双床位中第一床位的投影数据,并将扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Y的数据量,并将模拟双床位投影数据Y中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为双床位扫描数据;
根据扫描数据,获得模拟双床位投影数据,可以包括:
将扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Z的数据量,并将模拟双床位投影数据Z中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
在一个示例性的实现过程中,所述目标校正因子包括双床位活性浓度校正因子、层灵敏度校正因子中的至少一个。
本发明实施例还提供一种PET系统,该PET系统包括探测器、扫描床和控制台设备,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;其中,
所述闪烁晶体,用于在所述PET系统扫描过程中,探测被检体内发出的高能光子,并将所述高能光子转换为光信号;
所述光电转换器件,用于将所述光信号转换成电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述控制台设备,用于获得对模体进行PET扫描的扫描数据;根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;基于所述目标重建图像,计算所述目标重叠层数对应的目标校正因子;利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
所述控制台设备可以用于执行前述的任一种校正方法,此处不再赘述。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现如下操作:
获得对模体进行PET扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子;
利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
在一个示例性的实现过程中,获得对模体进行PET扫描的扫描数据,包括:
获得对模体进行单个床位扫描的扫描数据;或者
获得对模体进行双床位扫描的扫描数据。
在一个示例性的实现过程中,所述扫描数据为单个床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:
根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据作为双床位中第一床位的投影数据,将所述扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为双床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:
根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将所述扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为单个床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:
将扫描数据作为双床位中第一床位的投影数据,并将扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Y的数据量,并将模拟双床位投影数据Y中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
在一个示例性的实现过程中,所述扫描数据为双床位扫描数据;
根据扫描数据,获得模拟双床位投影数据,可以包括:
将扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Z的数据量,并将模拟双床位投影数据Z中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
在一个示例性的实现过程中,所述目标校正因子包括双床位活性浓度校正因子、层灵敏度校正因子中的至少一个。
对于装置和设备实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本说明书方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本领域技术人员在考虑说明书及实践这里申请的发明后,将容易想到本说明书的其它实施方案。本说明书旨在涵盖本说明书的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本说明书的一般性原理并包括本说明书未申请的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本说明书的真正范围和精神由下面的权利要求指出。
应当理解的是,本说明书并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本说明书的范围仅由所附的权利要求来限制。
以上所述仅为本说明书的较佳实施例而已,并不用以限制本说明书,凡在本说明书的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本说明书保护的范围之内。

Claims (6)

1.一种校正方法,其特征在于,用于对正电子发射断层成像PET系统的待校正重建图像进行校正,所述方法包括:
获得对模体进行PET扫描的扫描数据,包括:获得对模体进行单个床位扫描的扫描数据;或者获得对模体进行双床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据作为双床位中第一床位的投影数据,并将所述扫描数据作为双床位中第二床位的投影数据;或者,将所述扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将所述扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据,其中,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子,所述目标校正因子包括双床位活性浓度校正因子、层灵敏度校正因子中的至少一个;
利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
2.根据权利要求1所述的方法,其特征在于,所述扫描数据为单个床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:
将扫描数据作为双床位中第一床位的投影数据,并将扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Y的数据量,并将模拟双床位投影数据Y中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
3.根据权利要求1所述的方法,其特征在于,所述扫描数据为双床位扫描数据;
根据扫描数据,获得模拟双床位投影数据,包括:
将扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
根据目标重叠层数,确定模拟双床位投影数据Z的数据量,并将模拟双床位投影数据Z中各位置的值初始化为指定值,所述指定值小于0;
将所述第一床位的投影数据复制到所述双床位中所述第一床位所在的第一区域;
根据目标重叠层数和探测器总环数,确定所述第一床位和所述第二床位的环偏移量,并根据所述环偏移量确定所述双床位中所述第二床位所在的第二区域;
对于所述第二区域中的每一位置,若该位置的当前值小于0,将该位置的值更新为所述第二床位的投影数据在该位置的值;若该位置的当前值大于或等于0,将该位置的值更新为所述第二床位的投影数据在该位置的值与所述当前值的和;将所述第二区域的值更新后所述双床位的数据作为模拟双床位投影数据。
4.一种校正装置,其特征在于,用于对正电子发射断层成像PET系统的待校正重建图像进行校正,所述装置包括:
数据获取模块,用于获得对模体进行PET扫描的扫描数据,包括:获得对模体进行单个床位扫描的扫描数据;或者获得对模体进行双床位扫描的扫描数据;
模拟模块,用于根据所述扫描数据,获得模拟双床位投影数据,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
重建模块,用于根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
确定模块,用于基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子,所述目标校正因子包括双床位活性浓度校正因子、层灵敏度校正因子中的至少一个;
校正模块,用于利用所述目标校正因子对所述待校正重建图像中各个像素值进行校正,获得校正后的重建图像;
其中,所述模拟模块,具体用于根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据作为双床位中第一床位的投影数据,并将所述扫描数据作为双床位中第二床位的投影数据;或者,将所述扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将所述扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据。
5.一种控制台设备,其特征在于,包括内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于连接PET系统的探测器,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;
所述存储器,用于存储校正逻辑对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行如下操作:
获得对模体进行PET扫描的扫描数据,包括:获得对模体进行单个床位扫描的扫描数据;或者获得对模体进行双床位扫描的扫描数据;
根据所述扫描数据,获得模拟双床位投影数据,包括:根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;
将所述扫描数据作为双床位中第一床位的投影数据,并将所述扫描数据作为双床位中第二床位的投影数据;或者,将所述扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将所述扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;
基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;
确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据,其中,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;
根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;
基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子,所述目标校正因子包括双床位活性浓度校正因子、层灵敏度校正因子中的至少一个;
利用所述目标校正因子对待校正重建图像中各个像素值进行校正,获得校正后的重建图像。
6.一种PET系统,其特征在于,包括探测器、扫描床和控制台设备,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;其中,
所述闪烁晶体,用于在所述PET系统扫描过程中,探测被检体内发出的高能光子,并将所述高能光子转换为光信号;
所述光电转换器,用于将所述光信号转换成电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述控制台设备,用于获得对模体进行PET扫描的扫描数据,包括:获得对模体进行单个床位扫描的扫描数据;或者获得对模体进行双床位扫描的扫描数据;根据所述扫描数据,获得模拟双床位投影数据,包括:根据目标重叠层数和探测器总环数,确定两个床位的环偏移量;将所述扫描数据作为双床位中第一床位的投影数据,并将所述扫描数据作为双床位中第二床位的投影数据;或者,将所述扫描数据中一个床位的扫描数据作为双床位中第一床位的投影数据,将所述扫描数据中另一个床位的扫描数据作为双床位中第二床位的投影数据;基于所述环偏移量确定所述第一床位与所述第二床位的重叠区域和非重叠区域;确定所述重叠区域对应的第一投影数据等于所述第一床位与所述第二床位中相应投影数据的和,确定所述非重叠区域对应的第二投影数据等于相应床位的投影数据,所述第一投影数据和所述第二投影数据组成所述模拟双床位投影数据,其中,所述模拟双床位投影数据对应的两个床位的重叠层数为目标重叠层数;根据所述模拟双床位投影数据进行图像重建,得到目标重建图像;基于所述目标重建图像,确定所述目标重叠层数对应的目标校正因子;利用所述目标校正因子对待校正重建图像中各个像素值进行校正,获得校正后的重建图像,所述目标校正因子包括双床位活性浓度校正因子、层灵敏度校正因子中的至少一个。
CN201910989334.9A 2019-10-17 2019-10-17 校正方法、装置、控制台设备及pet系统 Active CN110755100B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910989334.9A CN110755100B (zh) 2019-10-17 2019-10-17 校正方法、装置、控制台设备及pet系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910989334.9A CN110755100B (zh) 2019-10-17 2019-10-17 校正方法、装置、控制台设备及pet系统

Publications (2)

Publication Number Publication Date
CN110755100A CN110755100A (zh) 2020-02-07
CN110755100B true CN110755100B (zh) 2023-07-04

Family

ID=69332389

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910989334.9A Active CN110755100B (zh) 2019-10-17 2019-10-17 校正方法、装置、控制台设备及pet系统

Country Status (1)

Country Link
CN (1) CN110755100B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117503172A (zh) * 2022-07-29 2024-02-06 上海联影医疗科技股份有限公司 一种成像方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5745542A (en) * 1996-03-13 1998-04-28 Analogic Corporation Ring suppression filter for use in computed tomography systems
CN105389788A (zh) * 2015-10-13 2016-03-09 沈阳东软医疗系统有限公司 Pet多床位图像的重建方法及装置、合并方法及装置
CN107330949A (zh) * 2017-06-28 2017-11-07 上海联影医疗科技有限公司 一种伪影校正方法及系统
US9958559B1 (en) * 2017-07-11 2018-05-01 Uih America, Inc. Method and apparatus for automatic detection and correction of patient bed shift using intrinsic scintillation crystal radiations
CN108294774A (zh) * 2018-01-04 2018-07-20 沈阳东软医疗系统有限公司 一种pet校正扫描方法、装置及设备

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5405229B2 (ja) * 2008-08-11 2014-02-05 株式会社東芝 X線コンピュータ断層撮影装置
CN102177430B (zh) * 2008-10-10 2014-04-02 皇家飞利浦电子股份有限公司 使用经移位的几何结构改进ct图像采集的方法和设备
US9492322B2 (en) * 2009-11-16 2016-11-15 Alcon Lensx, Inc. Imaging surgical target tissue by nonlinear scanning
US9785246B2 (en) * 2010-10-06 2017-10-10 Nuvasive, Inc. Imaging system and method for use in surgical and interventional medical procedures

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5745542A (en) * 1996-03-13 1998-04-28 Analogic Corporation Ring suppression filter for use in computed tomography systems
CN105389788A (zh) * 2015-10-13 2016-03-09 沈阳东软医疗系统有限公司 Pet多床位图像的重建方法及装置、合并方法及装置
CN107330949A (zh) * 2017-06-28 2017-11-07 上海联影医疗科技有限公司 一种伪影校正方法及系统
US9958559B1 (en) * 2017-07-11 2018-05-01 Uih America, Inc. Method and apparatus for automatic detection and correction of patient bed shift using intrinsic scintillation crystal radiations
CN109031395A (zh) * 2017-07-11 2018-12-18 上海联影医疗科技有限公司 利用晶体辐射对患者床移位的自动检测和校正方法和装置
CN108294774A (zh) * 2018-01-04 2018-07-20 沈阳东软医疗系统有限公司 一种pet校正扫描方法、装置及设备

Also Published As

Publication number Publication date
CN110755100A (zh) 2020-02-07

Similar Documents

Publication Publication Date Title
EP3224801B1 (en) Multi-modality imaging system and method
CN102067176B (zh) 结合局部运动监测、校正和评估的辐射成像
CN109961834B (zh) 影像诊断报告的生成方法及设备
CN105125231B (zh) 一种pet图像环状伪影的去除方法和装置
CN103845073B (zh) 用于pet衰减校正的mr扫描选择
US10754048B2 (en) Correcting count loss in pet system
CN105193442A (zh) 一种pet扫描时间的确定方法和装置
US10772582B2 (en) Multi-modal emission tomography quality based on patient and application
CN104114091B (zh) 空间校正核图像重建
EP1891899A1 (en) Method and system for performing local tomography
CN101842806A (zh) 脏同位素pet重建
CN110755100B (zh) 校正方法、装置、控制台设备及pet系统
US20230230297A1 (en) Ai-enabled early-pet acquisition
US11874411B2 (en) Estimation of partially missing attenuation in time-of-flight positron emission tomography
CN104000618B (zh) 一种用环真光子数门控方法进行呼吸运动门控校正方法
CN108873047B (zh) 检测放射源活度的方法、系统、计算机设备和存储介质
Hu et al. MAP image reconstruction for arbitrary geometry PET systems with application to a prostate-specific scanner
CN111316327A (zh) 对移动对象的pet数据的衰减校正
US11663758B2 (en) Systems and methods for motion estimation in PET imaging using AI image reconstructions
CN108416743B (zh) 列表模式数据的排序方法、装置及计算机可读存储介质
CN108805947A (zh) Pet数据处理方法和设备、pet成像系统
EP4239566A1 (en) Computer-implemented method for determining nuclear medical image data sets, imaging device, computer program and electronically readable storage medium
CN115886852A (zh) Pet系统的校正方法、装置、存储介质及计算机设备
US20120250965A1 (en) Method and System for Using a Modified Ordered Subsets Scheme for Attenuation Weighted Reconstruction
CN116898468A (zh) D-spect系统的动态质控方法

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

Effective date of registration: 20230414

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.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant