CN107137101A - 一种时间校准方法和装置 - Google Patents
一种时间校准方法和装置 Download PDFInfo
- Publication number
- CN107137101A CN107137101A CN201710270179.6A CN201710270179A CN107137101A CN 107137101 A CN107137101 A CN 107137101A CN 201710270179 A CN201710270179 A CN 201710270179A CN 107137101 A CN107137101 A CN 107137101A
- Authority
- CN
- China
- Prior art keywords
- crystal
- time
- deviation
- event
- die body
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000013078 crystal Substances 0.000 claims abstract description 248
- 230000004044 response Effects 0.000 claims abstract description 23
- 238000002591 computed tomography Methods 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims description 5
- 239000007787 solid Substances 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000002285 radioactive effect Effects 0.000 description 3
- 230000002902 bimodal effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000010365 information processing Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 210000000653 nervous system Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本公开提供一种时间校准方法和装置,其中方法用于对正电子发射型计算机断层显像PET系统的晶体采集事件时间进行校正,该方法包括:获取每个晶体对探测到的符合事件,得到晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;根据符合事件重建模体图像;根据模体图像,确定在晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;根据所述时间偏差,对所述晶体对的时间差进行预校正;根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。本公开使得晶体的采集时间更为准确。
Description
技术领域
本公开涉及核医学成像技术,特别涉及一种时间校准方法和装置。
背景技术
PET(Positron Emission Computed Tomography,正电子发射型计算机断层显像)是当今最先进的大型医疗诊断成像技术之一,PET除显示形态结构外,还能够利用活体代谢在分子水平上提供有关脏器及其病变的功能信息,在诊断肿瘤、心血管和神经系统等疾病中具有卓越性能。PET的工作原理是在受检对象中注射含有放射性核素的药物,放射性核素发生衰变产生正电子,正电子与周围的负电子湮灭产生一对背靠背的伽马光子,称为符合事件,PET系统通过晶体探测光子对重建出发射正电子的核素分布图。但是由于PET探测装置的性能不一致性等因素,导致各个晶体上对接收到伽马光子的时间判定精度不同,影响PET图像的重建,因此需要对晶体接收符合事件的时间进行校准。
现有的时间校准方案中,例如,可以在PET系统中心放置注源的模体,获取响应线(1ine-of-response,LOR)上的单峰时间-符合事件曲线,通过迭代算法来校准晶体上的时间偏差,这种方法通常都要求扫描的模体中心与PET系统的中心完全重合,但是实际中这种要求很难做到,模体中心与PET系统中心的偏移通常会导致时间偏差校准的精度大大降低,导致校正不准确。
发明内容
有鉴于此,本公开提供一种时间校准方法和装置,以提高时间校正的准确性。
具体地,本公开是通过如下技术方案实现的:
第一方面,提供一种时间校准方法,所述方法用于对正电子发射型计算机断层显像PET系统的晶体采集事件时间进行校正,所述方法包括:
获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
根据所述符合事件重建模体图像;
根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
根据所述时间偏差,对所述晶体对的时间差进行预校正;
根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
第二方面,提供一种时间校准装置,所述装置用于对正电子发射型计算机断层显像PET系统的晶体采集事件时间进行校正,所述装置包括:
信息处理模块,用于获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
图像重建模块,用于根据所述符合事件重建模体图像;
偏差确定模块,用于根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
预校正模块,用于根据所述时间偏差,对所述晶体对的时间差进行预校正;
迭代计算模块,用于根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
校正处理模块,用于根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
第三方面,提供一种PET设备,所述设备包括:PET探测器、处理器和存储器;所述探测器中包括多个晶体,所述晶体用于探测源自被检体的符合事件;
所述处理器用于读取存储器上的时间校准控制逻辑对应的机器可读指令,并执行所述指令以实现如下操作:
获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
根据所述符合事件重建模体图像;
根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
根据所述时间偏差,对所述晶体对的时间差进行预校正;
根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
第四方面,提供一种PET时间校准系统,所述系统包括:
处理设备;
计算机可读存储介质,所述介质上存储有时间校准控制逻辑对应的机器可读指令,所述指令能够被所述处理设备执行,所述指令用于实现如下操作:
获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
根据所述符合事件重建模体图像;
根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
根据所述时间偏差,对所述晶体对的时间差进行预校正;
根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
本公开提供的时间校准方法和装置,通过根据放射源中心位置与响应线上的晶体对中心位置之间的源位置偏差,得到晶体对的时间差的校正量,进而据此迭代计算得到每个晶体的实际时间偏差校正量,使得晶体的采集时间更为准确。
附图说明
图1是本公开一示例性实施例示出的一种PET系统的结构示意图;
图2是本公开一示例性实施例示出的一种时间校准方法的流程图;
图3是本公开一示例性实施例示出的一种位置偏差示意图;
图4是本公开一示例性实施例示出的晶体关联扇面的示意图;
图5是本公开一示例性实施例示出的单峰直方图;
图6是本公开一示例性实施例示出的一种PET设备的结构示意图;
图7是本公开一示例性实施例示出的一种时间校准装置的结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
图1是一个例子中的PET系统100的结构示意图,如图1所示,该PET系统100可以包括探测器101(PET detector),该探测器101可以是一个环形探测器,其中包括很多的闪烁晶体和光电倍增管。当扫描被检体时,被检体可以被放置在扫描床102上,由扫描床102带动至环形的探测器101中进行扫描。
被检体在被扫描前可以注射了含有放射性核素的药物,放射性核素发生衰变产生正电子,正电子与被检体内周围的负电子湮灭产生一对背靠背的伽马光子,光子可以被探测器101中的闪烁晶体接收到。光电倍增管可以将晶体探测到的事件信息传输至处理单元103,处理单元103可以判定各个晶体探测到的事件中哪些是来源于同一个湮灭事件,并将源自同一湮灭事件的两个光子称为符合事件。符合事件的信息可以被送往图像重建单元104,用于重建PET图像。
其中,探测到符合事件的晶体对,每个晶体都可以记录探测到事件的时间即晶体采集事件时间,实际实施中,由于PET探测装置的性能不一致性等因素,导致各个晶体上对接收到伽马光子的时间判定精度不同,晶体采集事件的时间不太准确。然而,该时间在PET的信息处理中非常重要,例如,以TOF-PET(Time of flight-PET)为例,PET系统可以根据上述晶体对接收到光子对的时间偏差估算出正电子发生湮灭的大体位置,然后用于图像重建,因而时间准确性关系到湮灭位置确定的准确性以及图像重建的质量。有鉴于此,本公开提供的时间校准方法要对PET探测器中的各个晶体采集事件的时间进行校准。
PET晶体的时间校准,通常可以在PET正式使用前进行,并且,可以是在PET系统中心附近放置一个作为放射源的对称性模体,该模体可以是实心的棒源、注源的实心圆柱模体或者圆筒状模体。当扫描该模体时,PET系统的探测器101就可以获得到伽马光子,并通过信息处理得到各个符合事件。通过扫描模体得到的符合事件的信息就可以对晶体的采集时间偏差进行校正。
图2示例了本公开的时间校准方法的流程图,可以包括:
在步骤201中,扫描对称性放射源模体,得到每个晶体对探测到的符合事件,并获取到所述晶体对接收到所述符合事件的时间差。
本步骤中,以晶体A和晶体B这一晶体对接收到一个符合事件为例,该符合事件中的两个伽马光子源自扫描模体时的同一湮灭事件,且该两个伽马光子分别被晶体A和晶体B接收到。假设晶体A接收到事件的时间是TA,晶体B接收到事件的时间是TB,那么这对晶体对应的时间差是TA-TB。
在步骤202中,根据符合事件重建模体图像。
本步骤中,可以根据符合事件重建模体图像。此外,还可以计算模体的中心位置在哪里。例如,可以利用扫描得到的数据重建模体图像,对重建图像进行滤波平滑,利用质心法求出模体的每一层图像的对称中心,本例子可以求出每一层图像的模体对称中心位置(xp,yp),质心法计算公式如下:
上述的公式1中,例如,i表示像素序号,xi表示像素i对应的x坐标,yi表示像素i对应的y坐标,Ii表示像素i对应的像素值,最终得到模体在每一层图像上的中心为P(xP,yP,zp),其中,Zp是所计算的图像层的Z坐标值,即所述图像层在轴向上的坐标值。
在步骤203中,根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差。
本例子中,PET探测晶体采集事件时间的偏差,其中一个导致该偏差的原因即源位置偏差,如果模体的放射源中心位置正好是PET系统的中心位置,那么将不会再有源位置偏差导致的采集时间偏移。因此,本步骤将计算所述的源位置偏差导致的时间偏差。
结合图3所示,仍以晶体A和晶体B为例,晶体A、B可以是不同探测器环上的晶体,作为放射源的模体31是一个实心的圆柱模体,晶体A和晶体B轴向位置的中心位置对应的放射源中心位置是P(xP,yP,zp),其中zp是A和B轴向上的中心Z坐标。晶体A和晶体B是探测到符合事件的晶体对,该晶体对的连线即响应线LOR32,点C是响应线上的晶体对中心位置。D(xD,yD,zD),是晶体对A和B之间与模体相交部分的中心点。
理想状态下,D应与C重合,此时晶体对A和B接收符合事件的时间差是零,但是如图3所示的状态,D与C之间存在偏移距离,导致晶体对A和B接收符合事件的时间差不为零。结合晶体对A和B的符合事件的光子可以是由响应线32与模体31交叠部分发射出来,那么D与C之间的距离偏移导致的晶体对A和B的实际时间差,即为源位置偏移导致的晶体采集时间偏差。
本步骤中,可以先计算D与晶体对中心位置C之间的距离偏差,假设A和B的坐标为A(xA,yA,zA)和B(xB,yB,zB),利用D在AB上且是AB与模体相交的中点(所述的AB是晶体A和晶体B之间的连线),对于模体考虑到在较小的轴向范围内可以保证轴向偏差极小,这里假设轴向偏差为0,可以按照如下公式计算:
而|AC|=0.5*|AB|,则上述D和C这两个位置点相对于晶体A的距离偏差为:
转变为D和C相对于晶体A的时间偏差为:
本例子中,是以计算晶体A的实际时间偏差校正量为例进行描述,其他各个晶体的时间偏差校正量可以按照同样的方法计算。当以晶体A的偏差计算为例时,如上述的公式2至公式4所示,可以计算得到伽马光子在由湮灭位置向晶体A行进的路线中由于源位置偏差导致的行进距离偏差,并进而得到伽马光子的行进时间的偏移。在图3所示的本例子中,是以相对于晶体对中的晶体A的时间偏差计算为例。
在步骤204中,根据所述时间偏差,对所述晶体对的时间差进行预校正。
本步骤仍以校正TOF-PET中的晶体A为例,TOFA_AB=TA-TB,其中,TA是晶体A接收符合事件的时间,TB是晶体B接收符合事件的时间。
请结合图4来看,图4所示为晶体A所关联的扇面,在PET探测器的各个晶体探测到的符合事件中,可能有多个符合事件都涉及到晶体A。例如,图4中的扇面上的晶体B的数量是多个,假设晶体B包括晶体B1、B2、B3等,晶体A和晶体B1是接收到符合事件的一个晶体对,晶体A和晶体B2是接收到另一符合事件的晶体对,晶体A和晶体B3是接收到又一符合事件的晶体对。
上述的每一个晶体对都对应一个接收符合事件的时间差,如果根据上述的晶体A及其关联扇面上的各个晶体之间的时间差,做成直方图,将得到图5所示的形式。图5是模体为实心棒源或者注源的实心圆柱模体时,对应的单峰直方图,横轴是时间差,纵轴是计数。晶体A对应的晶体对的时间差,可以是图5所示的直方图的峰值位置,比如,可以通过高斯函数拟合确定峰值位置,所述的时间差峰值位置即TOFA_AB。此外,当模体是圆筒状模体时,得到的直方图可以是双峰直方图,对于双峰直方图,可以先计算两个单峰的峰值,再求峰值的平均值作为晶体对中相对于晶体A的时间差TOFA_AB。
本步骤中,可以根据步骤203中得到的时间偏差,对晶体对的时间差进行预校正。如下的公式5所示,由于放射源中心D在晶体对中心C的右侧,更靠近晶体B,所以实际上晶体A采集的时间是延迟的。
TOFA_AB=TOFA_AB-ΔA_AB.......(公式5)
在步骤205中,根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量。
利用步骤(三)中计算的晶体对上的预校正的时间偏差峰值位置,例如晶体对A和B相对于晶体A的时间偏差峰值位置TOFA_AB,根据选择的晶体A关联的扇形区域,迭代估计晶体A的偏移量。
设定晶体A的时间偏差初始值按照如下公式迭代:
其中,是晶体A的时间偏差,BiasA是一个时间偏差临时变量,n是关联扇形面上的晶体个数,TOFA_AB是指预校正后的时间偏差,FanA是晶体A关联的扇形区域。迭代可以根据达到预设的最大迭代次数终止,或者根据满足都小于预设的一个阈值终止。
最后获得每个晶体的时间偏差校正量即最后迭代时间偏移校正量,N为迭代次数。
在步骤206中,根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
本步骤中,可以根据步骤205中得到的实际时间偏差校正量校正晶体A的符合事件采集时间,例如,可以按照如下公式校正:
此外,本例子还可以根据获得的每个晶体的实际时间偏差校正量建立时间偏差校正表,表中值为各个晶体及其对应的当PET正式使用时,再次获取到各个晶体的符合事件采集时间时,可以查找该时间偏差校正表,得到晶体对应的实际时间偏差校正量,据此对晶体采集的事件时间信息上延迟
本例子的时间校准方法,通过根据放射源中心位置与响应线上的晶体对中心位置之间的源位置偏差,得到晶体对的时间差的校正量,进而据此迭代计算得到每个晶体的实际时间偏差校正量,使得晶体的采集时间更为准确。
上述图2所示流程中的各个步骤,其执行顺序不限制于流程图中的顺序。此外,各个步骤的描述,可以实现为软件、硬件或者其结合的形式,例如,本领域技术人员可以将其实现为软件代码的形式,可以为能够实现所述步骤对应的逻辑功能的计算机可执行指令。当其以软件的方式实现时,所述的可执行指令可以存储在存储器中,并被设备中的处理器执行。
例如,参见图6所示,对应于上述方法,本公开同时提供一种PET设备,该设备可以包括PET探测器1701、处理器1702以及存储器1703,其中,探测器中包括多个晶体,所述晶体用于探测源自被检体的符合事件。处理器1702和存储器1703通常借由内部总线1704相互连接。在其他可能的实现方式中,所述设备还可能包括外部接口1705,以能够与其他设备或者部件进行通信。进一步地,存储器1703上存储有时间校准的控制逻辑1706,该控制逻辑1706从功能上划分的逻辑模块,可以是图7所示的时间校准装置的结构。
如图7所示,该时间校准装置可以包括:信息处理模块71、图像重建模块72、偏差确定模块73、预校正模块74、迭代计算模块75和校正处理模块76。
信息处理模块71,用于获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
图像重建模块72,用于根据所述符合事件重建模体图像;
偏差确定模块73,用于根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
预校正模块74,用于根据所述时间偏差,对晶体对的时间差进行预校正;
迭代计算模块75,用于根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
校正处理模块76,用于根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
此外,图2所示的时间校准流程还可以被包括在计算机可读存储介质中,所述介质上存储有时间校准控制逻辑对应的机器可读指令,并且这些介质可以与执行指令的处理设备连接,介质上存储的所述指令能够被处理设备执行。例如,本公开还可以提供一种PET时间校准系统,该系统可以包括处理设备和计算机可读存储介质,所述介质上存储的时间校准控制逻辑对应的机器可读指令用于实现如下操作:
获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
根据所述符合事件重建模体图像;
根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
根据所述时间偏差,对所述晶体对的时间差进行预校正;
根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
在本公开中,计算机可读存储介质可以是多种形式,比如,在不同的例子中,所述机器可读存储介质可以是:RAM(Radom Access Memory,随机存取存储器)、易失存储器、非易失性存储器、闪存、存储驱动器(如硬盘驱动器)、固态硬盘、任何类型的存储盘(如光盘、dvd等),或者类似的存储介质,或者它们的组合。特殊的,所述的计算机可读介质还可以是纸张或者其他合适的能够打印程序的介质。使用这些介质,这些程序可以被通过电学的方式获取到(例如,光学扫描)、可以被以合适的方式编译、解释和处理,然后可以被存储到计算机介质中。
以上所述仅为本公开的较佳实施例而已,并不用以限制本公开,凡在本公开的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本公开保护的范围之内。
Claims (10)
1.一种时间校准方法,其特征在于,所述方法用于对正电子发射型计算机断层显像PET系统的晶体采集事件时间进行校正,所述方法包括:
获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
根据所述符合事件重建模体图像;
根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
根据所述时间偏差,对所述晶体对的时间差进行预校正;
根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
2.根据权利要求1所述的方法,其特征在于,所述根据模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,包括:
确定所述晶体对的响应线与模体相交部分的中点位置,并得到所述中点位置与所述响应线上的晶体对中心位置之间的距离偏差;
根据所述距离偏差,得到相对于晶体对其中一个晶体的源位置偏差。
3.根据权利要求2所述的方法,其特征在于,所述晶体对的时间差是相对于所述晶体对中的所述其中一个晶体的时间偏差;
所述根据所述时间偏差,对所述晶体对的时间差进行预校正,包括:将所述时间差减去所述时间偏差。
4.根据权利要求1所述的方法,其特征在于,所述根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间,包括:
将分别对应各个晶体的实际时间偏差校正量,存储在时间偏差校正表;
当获取到所述晶体的符合事件采集时间时,查找所述时间偏差校正表,得到所述晶体对应的实际时间偏差校正量,并对所述符合事件采集时间进行校正。
5.一种时间校准装置,其特征在于,所述装置用于对正电子发射型计算机断层显像PET系统的晶体采集事件时间进行校正,所述装置包括:
信息处理模块,用于获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
图像重建模块,用于根据所述符合事件重建模体图像;
偏差确定模块,用于根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
预校正模块,用于根据所述时间偏差,对所述晶体对的时间差进行预校正;
迭代计算模块,用于根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
校正处理模块,用于根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
6.根据权利要求5所述的装置,其特征在于,
所述偏差确定模块,在用于确定源位置偏差时,包括:确定所述晶体对的响应线与模体相交部分的中点位置,并得到所述中点位置与所述响应线上的晶体对中心位置之间的距离偏差;根据所述距离偏差,得到相对于晶体对其中一个晶体的源位置偏差。
7.根据权利要求6所述的装置,其特征在于,所述晶体对的时间差是相对于所述晶体对中的所述其中一个晶体的时间偏差;
所述预校正模块,具体用于将所述时间差减去所述时间偏差。
8.根据权利要求5所述的装置,其特征在于,
所述校正处理模块,具体用于:将分别对应各个晶体的实际时间偏差校正量,存储在时间偏差校正表;当获取到所述晶体的符合事件采集时间时,查找所述时间偏差校正表,得到所述晶体对应的实际时间偏差校正量,并对所述符合事件采集时间进行校正。
9.一种PET设备,其特征在于,所述设备包括:PET探测器、处理器和存储器;所述探测器中包括多个晶体,所述晶体用于探测源自被检体的符合事件;
所述处理器用于读取存储器上的时间校准控制逻辑对应的机器可读指令,并执行所述指令以实现如下操作:
获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
根据所述符合事件重建模体图像;
根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
根据所述时间偏差,对所述晶体对的时间差进行预校正;
根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
10.一种PET时间校准系统,其特征在于,所述系统包括:
处理设备;
计算机可读存储介质,所述介质上存储有时间校准控制逻辑对应的机器可读指令,所述指令能够被所述处理设备执行,所述指令用于实现如下操作:
获取每个晶体对探测到的符合事件,得到所述晶体对接收符合事件的时间差,所述符合事件产生于扫描PET系统中放置的模体;
根据所述符合事件重建模体图像;
根据所述模体图像,确定在所述晶体对的响应线上相对于晶体对其中一个晶体的源位置偏差,并根据源位置偏差得到对应所述晶体的时间偏差;
根据所述时间偏差,对所述晶体对的时间差进行预校正;
根据预校正后的时间差,迭代获得所述晶体的实际时间偏差校正量;
根据所述实际时间偏差校正量,校正所述晶体的符合事件采集时间。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710270179.6A CN107137101A (zh) | 2017-04-24 | 2017-04-24 | 一种时间校准方法和装置 |
CN201810362644.3A CN108542412B (zh) | 2017-04-24 | 2018-04-20 | 图像重建方法和pet设备 |
US15/960,923 US10719960B2 (en) | 2017-04-24 | 2018-04-24 | Reconstructing image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710270179.6A CN107137101A (zh) | 2017-04-24 | 2017-04-24 | 一种时间校准方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107137101A true CN107137101A (zh) | 2017-09-08 |
Family
ID=59774974
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710270179.6A Pending CN107137101A (zh) | 2017-04-24 | 2017-04-24 | 一种时间校准方法和装置 |
CN201810362644.3A Active CN108542412B (zh) | 2017-04-24 | 2018-04-20 | 图像重建方法和pet设备 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810362644.3A Active CN108542412B (zh) | 2017-04-24 | 2018-04-20 | 图像重建方法和pet设备 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10719960B2 (zh) |
CN (2) | CN107137101A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108542412A (zh) * | 2017-04-24 | 2018-09-18 | 沈阳东软医疗系统有限公司 | 图像重建方法和pet设备 |
CN108703769A (zh) * | 2018-06-14 | 2018-10-26 | 上海联影医疗科技有限公司 | Tof数据的校正方法、装置、系统和计算机可读存储介质 |
CN109124672A (zh) * | 2018-10-19 | 2019-01-04 | 上海联影医疗科技有限公司 | 随机符合事件的校正方法、装置、设备和存储介质 |
CN110301927A (zh) * | 2019-07-04 | 2019-10-08 | 东软医疗系统股份有限公司 | 确定晶体固有效率的方法、装置、存储介质及医疗设备 |
CN111012372A (zh) * | 2019-12-20 | 2020-04-17 | 东软医疗系统股份有限公司 | 确定飞行时间的方法、装置、介质及正电子断层扫描仪 |
CN111493915A (zh) * | 2020-04-20 | 2020-08-07 | 明峰医疗系统股份有限公司 | 一种用于pet的时间校正方法 |
CN112807009A (zh) * | 2021-02-02 | 2021-05-18 | 上海联影医疗科技股份有限公司 | 放射源摆位方法、系统、设备、电子装置和存储介质 |
CN112932515A (zh) * | 2021-01-29 | 2021-06-11 | 明峰医疗系统股份有限公司 | 一种用于tof-pet的时间校正方法 |
CN112998735A (zh) * | 2021-02-26 | 2021-06-22 | 中派科技(深圳)有限责任公司 | 一种扫描装置重建图像的系统和方法 |
CN112998737A (zh) * | 2021-02-26 | 2021-06-22 | 中派科技(深圳)有限责任公司 | 一种扫描装置的时间偏移校正系统及其时间偏移校正方法 |
CN113116371A (zh) * | 2019-12-30 | 2021-07-16 | 佳能医疗系统株式会社 | Pet装置及校准方法 |
CN113425320A (zh) * | 2021-06-02 | 2021-09-24 | 东软医疗系统股份有限公司 | 时间校准方法、装置、电子设备及存储介质 |
WO2022218414A1 (en) * | 2021-04-15 | 2022-10-20 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for correcting positron emission tomography devices |
WO2023035361A1 (zh) * | 2021-09-10 | 2023-03-16 | 苏州瑞派宁科技有限公司 | 图像重建方法、装置、系统及计算机可读存储介质 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110584698B (zh) * | 2019-09-23 | 2022-12-27 | 上海联影医疗科技股份有限公司 | 探测器质量控制效验方法、装置、计算机设备和存储介质 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5900636A (en) * | 1997-05-30 | 1999-05-04 | Adac Laboratories | Dual-mode gamma camera system utilizing single-photon transmission scanning for attenuation correction of PET data |
US6803579B2 (en) * | 2001-09-28 | 2004-10-12 | General Electric Company | Technique for removal of picket fence effect in PET imaging systems |
US7649176B2 (en) * | 2003-11-14 | 2010-01-19 | Siemens Medical Solutions Usa, Inc. | Method for improving clinical data quality in positron emission tomography |
WO2005091989A2 (en) * | 2004-03-19 | 2005-10-06 | Triumf, Operating As A Joint Venture By The Governors Of The University Of Alberta, The University Of British Columbia, ... | A method and apparatus for vetoing random coincidences in positron emission tomographs |
US7465927B2 (en) * | 2006-03-23 | 2008-12-16 | Siemens Medical Solutions Usa, Inc. | Attenuation correction for nuclear medical imaging scanners with simultaneous transmission and emission acquisition |
US8089043B2 (en) * | 2009-12-01 | 2012-01-03 | Siemens Medical Solutions Usa, Inc. | Systems and methods for calibrating time alignment for a positron emission tomography (PET) system |
US8620054B2 (en) * | 2010-09-08 | 2013-12-31 | SOCPRA—Sciences Sante et Humaines, s.e.c. | Image reconstruction based on accelerated method using polar symmetries |
WO2013177661A1 (en) * | 2012-05-29 | 2013-12-05 | University Of Manitoba | Systems and methods for improving the quality of images in a pet scan |
US8796637B1 (en) * | 2013-05-24 | 2014-08-05 | Kabushiki Kaisha Toshiba | Timing calibration for time-of-flight (TOF) PET using positron-emitting isotopes and annihilation targets |
US9801595B2 (en) * | 2014-09-08 | 2017-10-31 | Toshiba Medical Systems Corporation | Count-weighted least squares parameter estimation for a photon-counting detector |
US9606245B1 (en) * | 2015-03-24 | 2017-03-28 | The Research Foundation For The State University Of New York | Autonomous gamma, X-ray, and particle detector |
CN106539591B (zh) * | 2015-09-21 | 2020-06-26 | 上海联影医疗科技有限公司 | Pet飞行时间状态质量检测方法和pet扫描装置 |
CN105913397B (zh) * | 2016-04-13 | 2018-10-02 | 沈阳东软医疗系统有限公司 | 一种重建图像的修正方法、装置及设备 |
CN107137101A (zh) * | 2017-04-24 | 2017-09-08 | 沈阳东软医疗系统有限公司 | 一种时间校准方法和装置 |
-
2017
- 2017-04-24 CN CN201710270179.6A patent/CN107137101A/zh active Pending
-
2018
- 2018-04-20 CN CN201810362644.3A patent/CN108542412B/zh active Active
- 2018-04-24 US US15/960,923 patent/US10719960B2/en active Active
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108542412B (zh) * | 2017-04-24 | 2021-11-09 | 东软医疗系统股份有限公司 | 图像重建方法和pet设备 |
CN108542412A (zh) * | 2017-04-24 | 2018-09-18 | 沈阳东软医疗系统有限公司 | 图像重建方法和pet设备 |
CN108703769A (zh) * | 2018-06-14 | 2018-10-26 | 上海联影医疗科技有限公司 | Tof数据的校正方法、装置、系统和计算机可读存储介质 |
CN109124672A (zh) * | 2018-10-19 | 2019-01-04 | 上海联影医疗科技有限公司 | 随机符合事件的校正方法、装置、设备和存储介质 |
CN109124672B (zh) * | 2018-10-19 | 2021-08-17 | 上海联影医疗科技股份有限公司 | 随机符合事件的校正方法、装置、设备和存储介质 |
CN110301927A (zh) * | 2019-07-04 | 2019-10-08 | 东软医疗系统股份有限公司 | 确定晶体固有效率的方法、装置、存储介质及医疗设备 |
CN111012372A (zh) * | 2019-12-20 | 2020-04-17 | 东软医疗系统股份有限公司 | 确定飞行时间的方法、装置、介质及正电子断层扫描仪 |
CN111012372B (zh) * | 2019-12-20 | 2023-05-09 | 沈阳智核医疗科技有限公司 | 确定飞行时间的方法、装置、介质及正电子断层扫描仪 |
CN113116371B (zh) * | 2019-12-30 | 2024-04-19 | 佳能医疗系统株式会社 | Pet装置及校准方法 |
CN113116371A (zh) * | 2019-12-30 | 2021-07-16 | 佳能医疗系统株式会社 | Pet装置及校准方法 |
CN111493915B (zh) * | 2020-04-20 | 2021-02-02 | 明峰医疗系统股份有限公司 | 一种用于pet的时间校正方法 |
CN111493915A (zh) * | 2020-04-20 | 2020-08-07 | 明峰医疗系统股份有限公司 | 一种用于pet的时间校正方法 |
CN112932515A (zh) * | 2021-01-29 | 2021-06-11 | 明峰医疗系统股份有限公司 | 一种用于tof-pet的时间校正方法 |
CN112807009A (zh) * | 2021-02-02 | 2021-05-18 | 上海联影医疗科技股份有限公司 | 放射源摆位方法、系统、设备、电子装置和存储介质 |
CN112998737A (zh) * | 2021-02-26 | 2021-06-22 | 中派科技(深圳)有限责任公司 | 一种扫描装置的时间偏移校正系统及其时间偏移校正方法 |
CN112998735B (zh) * | 2021-02-26 | 2022-09-02 | 中派科技(深圳)有限责任公司 | 一种扫描装置重建图像的系统和方法 |
CN112998735A (zh) * | 2021-02-26 | 2021-06-22 | 中派科技(深圳)有限责任公司 | 一种扫描装置重建图像的系统和方法 |
WO2022218414A1 (en) * | 2021-04-15 | 2022-10-20 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for correcting positron emission tomography devices |
CN113425320A (zh) * | 2021-06-02 | 2021-09-24 | 东软医疗系统股份有限公司 | 时间校准方法、装置、电子设备及存储介质 |
CN113425320B (zh) * | 2021-06-02 | 2023-09-26 | 沈阳智核医疗科技有限公司 | 时间校准方法、装置、电子设备及存储介质 |
WO2023035361A1 (zh) * | 2021-09-10 | 2023-03-16 | 苏州瑞派宁科技有限公司 | 图像重建方法、装置、系统及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108542412A (zh) | 2018-09-18 |
US10719960B2 (en) | 2020-07-21 |
CN108542412B (zh) | 2021-11-09 |
US20180308261A1 (en) | 2018-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107137101A (zh) | 一种时间校准方法和装置 | |
CN101454801B (zh) | 基于列表模式数据的局部运动补偿 | |
US9474495B2 (en) | System and method for joint estimation of attenuation and activity information | |
Tong et al. | Image reconstruction for PET/CT scanners: past achievements and future challenges | |
CN104434160B (zh) | Pet扫描装置及其时间偏移校正方法 | |
JP5290145B2 (ja) | ダーティアイソトープのpet再構成 | |
CN103329168A (zh) | 针对spect/ct系统的迭代锥形射束ct重建的截断补偿 | |
CN103824266A (zh) | 具有增强噪声控制滤波的迭代重建 | |
CN105125231A (zh) | 一种pet图像环状伪影的去除方法和装置 | |
CN101305297A (zh) | 使用解剖列表模式掩模的pet成像 | |
CN105496436B (zh) | 用于pet装置的时间校正方法和装置 | |
CN107976706A (zh) | 一种pet系统的计数丢失校正方法和装置 | |
US10036817B2 (en) | Solving outside-field of view scatter correction problem in positron emission tomography via digital experimentation | |
Lage et al. | Recovery and normalization of triple coincidences in PET | |
CN110327067A (zh) | 图像重建方法、装置、终端设备及pet系统 | |
US7769217B2 (en) | Fast iterative 3D PET image reconstruction using a set of 2D linogram transformations | |
CN108387922A (zh) | 用于检测器输出调整的系统和方法 | |
CN105913397B (zh) | 一种重建图像的修正方法、装置及设备 | |
CN101910869B (zh) | 多段重建 | |
Cheng et al. | Maximum likelihood activity and attenuation estimation using both emission and transmission data with application to utilization of Lu‐176 background radiation in TOF PET | |
CN113520435B (zh) | 时间校正方法、装置、计算机设备以及存储介质 | |
Pal et al. | 2D linear and iterative reconstruction algorithms for a PET-insert scanner | |
CN111080737B (zh) | 图像重建方法、装置及pet扫描系统 | |
CN107251095A (zh) | 图像重建系统、方法和计算机程序 | |
JP2006189274A (ja) | 核医学イメージング装置 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170908 |
|
WD01 | Invention patent application deemed withdrawn after publication |