CN112151162A - Pet时间实时校正及重建方法和系统 - Google Patents
Pet时间实时校正及重建方法和系统 Download PDFInfo
- Publication number
- CN112151162A CN112151162A CN202010992660.8A CN202010992660A CN112151162A CN 112151162 A CN112151162 A CN 112151162A CN 202010992660 A CN202010992660 A CN 202010992660A CN 112151162 A CN112151162 A CN 112151162A
- Authority
- CN
- China
- Prior art keywords
- time
- crystal
- histogram distribution
- correction
- pet
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000013078 crystal Substances 0.000 claims description 139
- 238000004364 calculation method Methods 0.000 claims description 31
- 238000005192 partition Methods 0.000 claims description 24
- 230000005855 radiation Effects 0.000 claims description 24
- 238000009825 accumulation Methods 0.000 claims description 19
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 2
- 238000002059 diagnostic imaging Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000000790 scattering method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
- G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
-
- 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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Primary Health Care (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Epidemiology (AREA)
- General Business, Economics & Management (AREA)
- Nuclear Medicine (AREA)
Abstract
本发明涉及医学影像设备技术领域,尤其涉及一种PET时间实时校正及重建方法和系统,分辨率计算过程联合了图像和计数的多重信息,使分辨率计算结果的可靠性和准确性得到加强,同时逐步迭代提高计算精度,保证了结果的稳定性和准确度,在PET环内的检测器之间的差异大时,本方法有着更好的适应性,能提升重建图像的质量。
Description
技术领域
本发明涉及医学影像设备技术领域,尤其涉及一种PET时间实时校正及重建方法和系统。
背景技术
正电子发射断层设备的原理是测量正电子湮没后的γ光子配对,并利用的γ光子对的几何方向和时间特性进行成像。其中γ光子对的时间特性是指,两个方向相反的光子由于飞行时间不同,进入检测器闪烁体晶体的时间产生差异。这个时间差产生的原因是由于晶体材质,晶体形态和前、后端电子学性能的影响,在每一个晶体上产生一定的偏差,从而使符合数据的时间值和真实值之间有所差异。
为了消除这种差异,文献[M E Werner and J S Karp,TOF PET offsetcalibration from clinical data.Phys.Med.Biol.58(2013)4031–4046]中使用了一种基于nontof的结果上来统计时间偏差的方法,该方法设定了一个固定的最小时间差,然后使用Histgram的方法统计每个晶体上符合事件分布和在nontof图像的时间分布对比处理。该方案存在的以下问题:1.最小精度由最小时间差决定,一旦确定就无法继续提高精度,但由于最小时间差是一个固定值,在晶体时间差异很大时,无法统计,造成严重误差或失效;2.在计算时间校正时需要重新计算scatter,浪费了计算效率。
发明内容
为解决上述问题,本发明的目的在于提供一种PET时间实时校正及重建方法和系统,通过使用经过校正处理后的List mode数据,在不需要重新进行scatter计算或者其他插值计算的情况下,快速计算出晶体延时,同时逐步迭代提高计算精度,保证了结果的稳定性和准确度。
为了实现上述目的,本发明的技术方案如下:
一种PET时间实时校正及重建方法,包括以下步骤:
S1.收集符合事件,形成可读取和使用的List mode数据,并确定随机分数randFra和散射分数scatFra;
S2.对List mode数据进行无飞行时间的重建,获取辐射源的空间分布Voxel[x,y,z],其中x y z为空间三维坐标;
S3.累积计数List mode数据的发生次数,获取晶体的时间直方图分布H1;
累积计数List mode数据发生的校正信息,获取晶体的时间直方图分布H2;
根据辐射源的空间分布的体数据计算得到晶体的时间直方图分布L1;
S4.根据晶体的时间直方图分布,计算晶体时间基准;
S5.对晶体时间基准进行迭代更新,并设置迭代终止条件,若不满足迭代终止条件,则根据更新后的晶体时间基准重新获取晶体的时间直方图分布H1和H2,并进行晶体时间基准的下一次更新;若满足迭代终止条件,则进行图像重建并输出图像。
进一步的,所述步骤S3中的晶体的时间直方图分布H1的获取包括以下步骤:
定义晶体的时间直方图分布H1为一个二维的存储空间H1[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为dt;
当符合事件j发生时,与它相关的晶体为c1和c2,时间差为tj,则偏移量tjd=[tj/dt],则分别在H1[c1,m-tjd]和H1[c2,m+tjd]上分别记录,记录的方式采用自加,其中m是指的n的1/2,即分隔数的一半。
进一步的,所述步骤S3中的晶体的时间直方图分布H2的获取包括以下步骤:
定义晶体的时间直方图分布H2为一个二维的存储空间H2[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为变量dt;
当符合事件j发生时,与它相关的晶体为c1和c2,时间差为tj,则偏移量tjd=[tj/dt],则分别在H2[c1,m-tjd]和H2[c2,m+tjd]上分别记录,记录的方式是H2[c1,m-tjd]=H2[c1,m-tjd]+badd j,H2[c2,m+tjd]=H2[c2,m+tjd]+badd j,badd j是散射符合事件和随机符合事件的计算值,其中m是指的n的1/2,即分隔数的一半。
进一步的,所述步骤S3中的晶体的时间直方图分布L1的获取包括以下步骤:
定义晶体的时间直方图分布L1为一个二维的存储空间L1[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为变量dt;
统计所有从晶体出发的LOR线共有LORc条,其中任意一条LOR线上,记录LOR上任意一个位置的图像像素与LOR的中心位置的距离为d,得到时间直方图分布L1的时间差td=d/Vc,偏移量tdoffset=[d/Vc/dt],其中Vc表示光速常量,tdoffset代表了以dt作为单位下的距离;
定义射线追踪法获取到晶体序号pix,偏移量tdoffset时L1的强度为
公式中,LORc是与晶体相关的所有LOR,xi(t),yi(t),zi(t)是晶体pix的第i条LOR上时间差为t的三维坐标,可以通过晶体的几何位置快速推导,σ是系统的时间分辨率,t代表时间变量。
进一步的,所述步骤S4中的晶体时间基准的计算包括以下步骤:
定义L2[c,n]是一个二维的存储空间,则
通过卷积,计算第K个晶体的晶体时间基准tdk,
公式中,h为自变量,优化求解最大值算法argmax使用线搜索(line search)实现。
进一步的,所述步骤S5中对晶体时间基准的迭代更新通过下式实现:
td’k=∈*tdk0+(1-∈)*tdk,
公式中,td’k代表更新后的晶体时间基准,tdk0代表上一次晶体时间基准,∈是经验参数,代表了更新的速度;
根据更新后的td’k,选出晶体c1和c2的晶体时间基准td’c1和td’c2,计算H1和H2中的时间差tj,tj=tj0+td’c1-td’c2,从而可以获得更新后的H1和H2,tj0为H1和H2中的初始时间差。
进一步的,所述时间宽度dt为变量,且随迭代次数的增加逐次减小,时间宽度dt的最小值为信号采集时间刻度;每迭代一次,即更新一次dt值,并将更新值带入H1、H2和L1中的偏移量的计算式中,进行偏移量的更新。
进一步的,所述时间宽度dt的变化遵循dt[ir]=dt0*pir,其中ir为迭代序号,p为[0 1]之间的分数,dt0为初始的间隔宽度。
一种PET时间实时校正及重建系统,包括
PET数据采集模块,用于收集符合事件,形成可以读取并使用的List mode数据;
辐射源空间分布获取模块,对List mode数据进行无飞行时间的重建,获得辐射源的空间分布;
符合事件累积模块,累积计数list mode数据的发生次数,得到晶体的时间直方图分布;
校正系数计算模块,累积计数list mode数据的发生的校正信息,得到晶体的时间直方图分布;
LOR线追踪累积模块,根据辐射源的空间分布的体数据,计算晶体的时间直方图分布;
晶体时间基准计算模块,根据符合事件累积模块、校正系数计算模块及LOR线追踪累积模块中得到的晶体的时间直方图分布,计算晶体时间基准;
迭代更新模块,用于对晶体时间基准进行迭代更新,若更新结果满足迭代终止条件则终止迭代,若更新结果不满足迭代终止条件,则根据更新结果重新计算晶体的时间直方图分布;
重建模块,用于图像重建并输出。
本发明的优点在于:分辨率计算过程联合了图像和计数的多重信息,使分辨率计算结果的可靠性和准确性得到加强,同时逐步迭代提高计算精度,保证了结果的稳定性和准确度,在PET环内的检测器之间的差异大时,本方法有着更好的适应性,能提升重建图像的质量。
附图说明
图1为List mode数据示意图;
图2为实施例中的校正及重建计算流程图;
具体实施方式
以下结合实施例对本发明作进一步详细描述。
对可能出现的参数进行背景介绍响应线:将探测器探测到的γ光子的两个晶体条之间的连线称为响应线(Line of Response,LOR)。
符合时间窗:是为两个γ光子到达探测器的时间差所设的时长。
符合事件:任意两个γ光子在符合时间窗内被探测器探测到的事件。
晶体:是指PET探测器环上检测器前端接受γ光子从而产生前端触发信号的最小几何单元。
真实符合:指湮没辐射产生的两个γ光子,以相反的方向运动并在符合时间窗内被两个晶体分别接收到的符合。
散射符合:指湮没辐射产生的两个γ光子,如果一个到达之前与组织发生散射,但仍在符合时间窗内被探测到的散射符合。
随机符合:是一种假符合,两个γ光子毫无时间和空间的相互关系,但在符合时间窗内被误探测到的符合事件。
List mode数据:是PET机器扫描过程中,记录符合事件的数据,它的特征是:其中的每一个事件至少包括以下信息:位置信息(探测到符合事件的两个γ光子的检测器序号,或者是这两个检测器组成的LOR序号),时间信息(两个光子被探测到的时间差),校正系数的数据(包括散射和随机事件),List mode数据示意图如附图1所示。
校正信息:具体描述可见参考文献[Wang W,Hu Z,Gualtieri E E,etal.Systematic and Distributed Time-of-Flight List Mode PET Reconstruction[C]//Nuclear Science Symposium Conference Record.IEEE,2006]的公式(2)中的baddj,其中的随机事件的数量和分布分别由符合延迟电路确定,具体操作是从延迟电路采集的事件首先通过乘法因子进行校正,重新绑定到正弦图中,然后使用低通滤波器平滑,参见文献[Rokitta O,Casey M,Wienhard K and Pictrzyk U 2000Random correction forpositron emission tomography using singles count rates IEEENucl.Sci.Symp.Conf.Record.3 37–40]。散射事件的数量和分布由单次散射方法估计计算,参见文献[C.C.Watson,“Extension of Single Scatter Simulation to ScatterCorrection of Time-of-Flight PET”,IEEE NSS/MIC,2005.]
随机分数:定义为随机符合事件数与总符合事件数的比例。
散射分数:定义为随机散射事件数与(散射事件数+真实符合事件数)的比例。信号采集时间刻度:是指符合事件的时间信息的最小记录单位。
PET诊断成像装置,包括:成像区域外部的多个辐射探测器环状结构,用来采集到发出的辐射事件的信息。探测器具有不同的定时延迟。校准体模可以发射时间同时但方向相反的辐射的事件对(符合事件)的辐射源,辐射事件对与辐射探测器的相互作用路径是LOR,当检测器环确定后,每一条LOR的几何位置是已知的。一旦每个检测器的定时延迟被校准,可以认为符合事件的标记时间是在真实符合发生的位置为中心的,并以特定宽度形成高斯分布。该高斯分布的半高全款宽度就是系统的时间分辨率。
本实施例的PET时间实时校正及重建系统由PET数据采集模块、辐射源空间分布获取模块、符合事件累积模块、校正系数计算模块、LOR线追踪累积模块、晶体时间基准计算模块、迭代更新模块及重建模块组成。
PET数据采集模块由外罩、机架及探测器组成,探测器包括晶体、光电放大和电子处理电路。探测器检测从被检体P内部放出的成对湮没γ射线,生成与检测出的成对湮没γ射线的光量相应的脉冲状电信号。后端信号处理器根据来自光电转换器件的电信号生成单事件数据(Single)。计数部从重复供给的单事件数据中重复确定容纳在与预先设定的时间范围内的2个单事件有关的事件数据。时间范围被设定为例如6ns~18ns左右。该成对的单事件被推测为由来于从同一成对湮没点产生的成对湮没γ射线。成对的单事件概括地被称为符合事件。符合事件被收集到后,事件相关的时间空间信息进行组合封装,并在磁盘或者其他数据存储介质上保存,形成了重建程序可以读取并使用的List mode数据。同时通过计算总符合、随机符合和估算散射符合的总数,确定出随机分数记为randFra,散射分数记为scatFra。计算随机分数和散射分数的概念和统计公知的,在文献[Wang W,Hu Z,GualtieriE E,et al.Systematic and Distributed Time-of-Flight List Mode PETReconstruction[C]//Nuclear Science Symposium Conference Record.IEEE,2006]中有说明出处。
辐射源空间分布获取模块运行于计算机设备,其中计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序。该模块的作用是提取辐射源的实际形态特征,这样可以获得到辐射源的空间分布。它是利用PET诊断成像装置获取到PET符合数据对信号,对信号进行无飞行时间的重建,重建方法例如,有序子集最大似然期望值法OSEM,最大似然期望值法(MLEM),滤波反投影FBP或者其他方法,从而得到了基本除去噪声和衰减干扰的质量较佳的光子空间分布图。
符合事件累积模块和校正系数计算模块将所有符合事件进行统计,得到每个晶体的时间直方图分布。
LOR线追踪累积模块,它的作用是从图像的角度计算晶体的时间直方图分布。
晶体时间基准计算模块,根据符合事件累积模块、校正系数计算模块及LOR线追踪累积模块中得到的晶体的时间直方图分布,计算晶体时间基准(即延时)。
迭代更新模块,用于对晶体时间基准进行迭代更新,若更新结果满足迭代终止条件则终止迭代,若更新结果不满足迭代终止条件,则根据更新结果重新计算符合事件累积模块、校正系数计算模块及LOR线追踪累积模块中的晶体的时间直方图分布。
重建模块,用于图像重建并输出。
如图2所示,本实施例的PET时间实时校正及重建方法包括以下步骤,
S1.PET数据采集模块收集符合事件,形成了可以读取并使用的List mode数据,并确定随机分数randFra和散射分数scatFra;
S2.辐射源空间分布获取模块利用List mode数据,获取辐射源的空间分布Voxel[x,y,z],其中x y z为空间三维坐标;
S3.符合事件累积模块累积计数List mode数据的发生次数,获取晶体的时间直方图分布H1,具体包括:
定义当第j个符合事件被检测器捕获后,记为事件j,定义晶体的时间直方图分布H1为一个二维的存储空间H1[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,也就是检测器上的总的晶体数量,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为dt;这样,当List mode数据输入后,经过累积运算就会输出H1。
当符合事件j发生时,与它相关的晶体为c1和c2,时间差为tj,则偏移量tjd=[tj/dt]([]表示对括号里的算术取整操作,取整方式可以是四舍五入,向上或者向下取整),则分别在H1[c1,m-tjd]和H1[c2,m+tjd]上分别记录,记录的方式采用自加,其中m是指的n的1/2,即分隔数的一半;
校正系数计算模块通过累积计数List mode数据发生的校正信息,获取晶体的时间直方图分布H2,具体包括:
定义晶体的时间直方图分布H2为一个二维的存储空间H2[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为变量dt,这样,当List mode数据输入后,经过累积校正信息就会输出H2;
当符合事件j发生时,与它相关的晶体为c1和c2,时间差为tj,则偏移量tjd=[tj/dt]([]表示对括号里的算术取整操作,取整方式可以是四舍五入,向上或者向下取整),则分别在H2[c1,m-tjd]和H2[c2,m+tjd]上分别记录,记录的方式是H2[c1,m-tjd]=H2[c1,m-tjd]+badd j,H2[c2,m+tjd]=H2[c2,m+tjd]+badd j,badd j是散射符合事件和随机符合事件的计算值,其中m是指的n的1/2,即分隔数的一半;
LOR线追踪累积模块根据辐射源的空间分布的体数据计算得到晶体的时间直方图分布L1,具体包括:
定义晶体的时间直方图分布L1为一个二维的存储空间L1[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为变量dt;
统计所有从晶体出发的LOR线共有LORc条,其中任意一条LOR线上,记录LOR上任意一个位置的图像像素与LOR中心的距离d,得到时间直方图分布L1的时间差td=d/Vc,偏移量tdoffset=[d/Vc/dt],其中Vc表示光速常量,Tdoffset代表了以dt作为单位下的距离;
定义射线追踪法获取到晶体序号pix,偏移量tdoffset时L1的强度为
公式中,LORc是与晶体相关的所有LOR,xi(t),yi(t),zi(t)是晶体pix的第i条LOR上时间差为t的三维坐标,可以通过晶体的几何位置快速推导,σ是系统的时间分辨率,t代表时间变量;
S4.晶体时间基准计算模块计算晶体时间基准(即晶体延时),具体包括:
定义L2[c,n]是一个二维的存储空间,则
公式中,h为自变量,优化求解最大值算法argmax使用线搜索(line search)实现;
S5.迭代更新模块进行迭代更新,具体包括:
设置迭代终止条件;
对晶体时间基准通过td′k=∈*tdk0+(1-∈)*tdk实现迭代更新,公式中,td’k代表更新后的晶体时间基准,tdk0代表上一次晶体时间基准,∈是经验参数,代表了更新的速度;根据更新后的td’k,选出晶体c1和c2的晶体时间基准td’c1和td’c2,计算晶体c1和c2的时间差tj,tj=tj0+td’c1-td’c2,tj0为晶体c1和c2中的初始时间差;
为提高计算精确度,本实施例中的时间宽度dt为变量,在第一次迭代过程中,dt一般取很大的值来保证计算范围,在后面的迭代中dt值会逐次减小,最小值是信号采集时间刻度,为方便计算,时间宽度dt的变化遵循dt[ir]=dt0*pir,其中ir为迭代序号,p为[0 1]之间的分数,dt0为初始的间隔宽度;
每一次更新tj和dt后,返回步骤S3,重新计算H1、H2及L1,使它们动态更新时间计算范围,从而达到更加精确计算的目的;迭代截止条件有不同的设置方法,如规定迭代的次数,或者计算tdk0和tdk之间差异值,差异值足够小时停止迭代;
步骤S6.迭代终止后,利用周知的list-mode-osem重建公式进行图像重建,重建公式如下:可参见文献[M E Werner andJ S Karp,TOF PET offset calibration from clinical data.Phys.Med.Biol.58(2013)4031–4046]。
上述实施例仅用于解释说明本发明的构思,而非对本发明权利保护的限定,凡利用此构思对本发明进行非实质性的改动,均应落入本发明的保护范围。
Claims (9)
1.一种PET时间实时校正及重建方法,其特征在于,包括以下步骤:
S1.收集符合事件,形成可读取和使用的List mode数据,并确定随机分数randFra和散射分数scatFra;
S2.对List mode数据进行无飞行时间的重建,获取辐射源的空间分布Voxel[x,y,z],其中x y z为空间三维坐标;
S3.累积计数List mode数据的发生次数,获取晶体的时间直方图分布H1;
累积计数List mode数据发生的校正信息,获取晶体的时间直方图分布H2;
根据辐射源的空间分布的体数据计算得到晶体的时间直方图分布L1;
S4.根据晶体的时间直方图分布,计算晶体时间基准;
S5.对晶体时间基准进行迭代更新,并设置迭代终止条件,若不满足迭代终止条件,则根据更新后的晶体时间基准重新获取晶体的时间直方图分布H1和H2,并进行晶体时间基准的下一次更新;若满足迭代终止条件,则进行图像重建并输出图像。
2.如权利要求1所述的一种PET时间实时校正及重建方法,其特征在于,所述步骤S3中的晶体的时间直方图分布H1的获取包括以下步骤:
定义晶体的时间直方图分布H1为一个二维的存储空间H1[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为dt;
当符合事件j发生时,与它相关的晶体为c1和c2,时间差为tj,则偏移量tjd=[tj/dt],则分别在H1[c1,m-tjd]和H1[c2,m+tjd]上分别记录,记录的方式采用自加,其中m是指的n的1/2,即分隔数的一半。
3.如权利要求1所述的一种PET时间实时校正及重建方法,其特征在于,所述步骤S3中的晶体的时间直方图分布H2的获取包括以下步骤:
定义晶体的时间直方图分布H2为一个二维的存储空间H2[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为变量dt;
当符合事件j发生时,与它相关的晶体为c1和c2,时间差为tj,则偏移量tjd=[tj/dt],则分别在H2[c1,m-tjd]和H2[c2,m+tjd]上分别记录,记录的方式是H2[c1,m-tjd]=H2[c1,m-tjd]+badd j,H2[c2,m+tjd]=H2[c2,m+tjd]+badd j,badd j是散射符合事件和随机符合事件的计算值,其中m是指的n的1/2,即分隔数的一半。
4.如权利要求3所述的一种PET时间实时校正及重建方法,其特征在于:所述步骤S3中的晶体的时间直方图分布L1的获取包括以下步骤:
定义晶体的时间直方图分布L1为一个二维的存储空间L1[c,n],该空间的大小为c*n,c代表的晶体的序号的最大值,n代表了出现在该晶体上的事件所记录的时间差的分隔数,n个分隔数中每一个分隔所对应的时间宽度为变量dt;
统计所有从晶体出发的LOR线共有LORc条,其中任意一条LOR线上,记录LOR上任意一个位置的图像像素与LOR的中心位置的距离为d,得到时间直方图分布L1的时间差td=d/Vc,偏移量tdoffset=[d/Vc/dt],其中Vc表示光速常量,tdoffset代表了以dt作为单位下的距离;
定义射线追踪法获取到晶体序号pix,偏移量tdoffset时L1的强度为
公式中,LORc是与晶体相关的所有LOR,xi(t),yi(t),zi(t)是晶体pix的第i条LOR上时间差为t的三维坐标,可以通过晶体的几何位置快速推导,σ是系统的时间分辨率,t代表时间变量。
6.如权利要求5所述的一种PET时间实时校正及重建方法,其特征在于,所述步骤S5中对晶体时间基准的迭代更新通过下式实现:
td’k=∈*tdk0+(1-∈)*tdk,
公式中,td’k代表更新后的晶体时间基准,tdk0代表上一次晶体时间基准,∈是经验参数,代表了更新的速度;
根据更新后的td’k,选出晶体c1和c2的晶体时间基准td’c1和td’c2,计算H1和H2中的时间差tj,tj=tj0+td’c1-td’c2,从而可以获得更新后的H1和H2,tj0为H1和H2中的初始时间差。
7.如权利要求6所述的一种PET时间实时校正及重建方法,其特征在于:所述时间宽度dt为变量,且随迭代次数的增加逐次减小,时间宽度dt的最小值为信号采集时间刻度;每迭代一次,即更新一次dt值,并将更新值带入H1、H2和L1中的偏移量的计算式中,进行偏移量的更新。
8.如权利要求7所述的一种PET时间实时校正及重建方法,其特征在于,所述时间宽度dt的变化遵循dt[ir]=dt0*pir,其中ir为迭代序号,p为[01]之间的分数,dt0为初始的间隔宽度。
9.一种PET时间实时校正及重建系统,其特征在于:包括
PET数据采集模块,用于收集符合事件,形成可以读取并使用的List mode数据;
辐射源空间分布获取模块,对List mode数据进行无飞行时间的重建,获得辐射源的空间分布;
符合事件累积模块,累积计数list mode数据的发生次数,得到晶体的时间直方图分布;
校正系数计算模块,累积计数list mode数据的发生的校正信息,得到晶体的时间直方图分布;
LOR线追踪累积模块,根据辐射源的空间分布的体数据,计算晶体的时间直方图分布;
晶体时间基准计算模块,根据符合事件累积模块、校正系数计算模块及LOR线追踪累积模块中得到的晶体的时间直方图分布,计算晶体时间基准;
迭代更新模块,用于对晶体时间基准进行迭代更新,若更新结果满足迭代终止条件则终止迭代,若更新结果不满足迭代终止条件,则根据更新结果重新计算晶体的时间直方图分布;
重建模块,用于图像重建并输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010992660.8A CN112151162B (zh) | 2020-09-21 | 2020-09-21 | Pet时间实时校正及重建方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010992660.8A CN112151162B (zh) | 2020-09-21 | 2020-09-21 | Pet时间实时校正及重建方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112151162A true CN112151162A (zh) | 2020-12-29 |
CN112151162B CN112151162B (zh) | 2022-10-21 |
Family
ID=73893448
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010992660.8A Active CN112151162B (zh) | 2020-09-21 | 2020-09-21 | Pet时间实时校正及重建方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112151162B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113670961A (zh) * | 2021-08-19 | 2021-11-19 | 南京航空航天大学 | 一种基于空间分布拟合的γ光子多次散射校正方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150262389A1 (en) * | 2012-11-08 | 2015-09-17 | The General Hospital Corporation | System and method for mutli-modality time-of-flight attenuation correction |
US20150297168A1 (en) * | 2014-04-18 | 2015-10-22 | Siemens Medical Solutions Usa, Inc. | Patient Based Detector Crystal Quality Control for Time of Flight Acquisition |
US20180021009A1 (en) * | 2014-10-27 | 2018-01-25 | Koninklijke Philips N.V. | Pet detector timing calibration |
US20180144513A1 (en) * | 2016-11-24 | 2018-05-24 | Shenyang Neusoft Medical Systems Co., Ltd. | Reconstructing pet image |
CN110988974A (zh) * | 2019-12-12 | 2020-04-10 | 明峰医疗系统股份有限公司 | 曲面辐射源的pet探测器飞行时间的校准系统 |
CN111493915A (zh) * | 2020-04-20 | 2020-08-07 | 明峰医疗系统股份有限公司 | 一种用于pet的时间校正方法 |
-
2020
- 2020-09-21 CN CN202010992660.8A patent/CN112151162B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150262389A1 (en) * | 2012-11-08 | 2015-09-17 | The General Hospital Corporation | System and method for mutli-modality time-of-flight attenuation correction |
US20150297168A1 (en) * | 2014-04-18 | 2015-10-22 | Siemens Medical Solutions Usa, Inc. | Patient Based Detector Crystal Quality Control for Time of Flight Acquisition |
US20180021009A1 (en) * | 2014-10-27 | 2018-01-25 | Koninklijke Philips N.V. | Pet detector timing calibration |
US20180144513A1 (en) * | 2016-11-24 | 2018-05-24 | Shenyang Neusoft Medical Systems Co., Ltd. | Reconstructing pet image |
CN110988974A (zh) * | 2019-12-12 | 2020-04-10 | 明峰医疗系统股份有限公司 | 曲面辐射源的pet探测器飞行时间的校准系统 |
CN111493915A (zh) * | 2020-04-20 | 2020-08-07 | 明峰医疗系统股份有限公司 | 一种用于pet的时间校正方法 |
Non-Patent Citations (5)
Title |
---|
ME WERNER等: "TOF PET offset calibration from clinical data", 《PHYSICS IN MEDICINE & BIOLOGY》 * |
STEFAAN VANDENBERGHE等: "Recent developments in time-of-flight PET", 《EJNMMI PHYSICS》 * |
XINGJIAN YU等: "Novel crystal timing calibration method based on total variation", 《PHYSICS IN MEDICINE & BIOLOGY》 * |
张斌等: "PET断层重建中动态射线追踪算法的实现", 《郑州大学学报(理学版)》 * |
梁国栋等: "PET图像空间结合位置偏移校正的点扩散恢复算法研究", 《中国医学装备》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113670961A (zh) * | 2021-08-19 | 2021-11-19 | 南京航空航天大学 | 一种基于空间分布拟合的γ光子多次散射校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112151162B (zh) | 2022-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8265365B2 (en) | Time of flight scatter distribution estimation in positron emission tomography | |
US8913810B2 (en) | Simultaneous reconstruction of emission activity and attenuation coefficient distribution from TOF data, acquired with external shell source | |
CN107978002B (zh) | 一种pet图像重建方法和装置 | |
US7381959B2 (en) | Technique for reconstructing PET scan images | |
JP5378787B2 (ja) | 粗い角度及びスライスリビニングによる三次元飛行時間pet | |
US7129496B2 (en) | Method and system for scattered coincidence estimation in a time-of-flight positron emission tomography system | |
Zaidi et al. | Scatter compensation techniques in PET | |
US9044153B2 (en) | Random sinogram variance reduction in continuous bed motion acquisition | |
CN106725573B (zh) | 一种连续进床pet系统的计数丢失校正方法和装置 | |
EP2748594B1 (en) | Radiographic apparatus for detecting photons with offset correction | |
EP2433261A2 (en) | Continuous time-of-flight scatter simulation method | |
CN110988974B (zh) | 曲面辐射源的pet探测器飞行时间的校准系统 | |
US9091771B2 (en) | System and method for improving detection of gamma interactions in a positron emission tomography system | |
CN106255994A (zh) | 针对正电子发射断层摄影(pet)列表模式迭代重建的重建中滤波 | |
WO2008075037A1 (en) | Method of and software for calculating a scatter estimate for tomographic scanning and system for tomographic scanning | |
JP7317586B2 (ja) | 医用画像処理装置、方法及びプログラム | |
WO2015006123A1 (en) | System and method to improve image quality of emission tomography when using advanced radionuclides | |
JP6125309B2 (ja) | 偶発同時計数推定方法及び偶発同時計数推定装置 | |
US20200033491A1 (en) | System and method for self-time alignment calibration for a positron emission tomography system | |
CN112151162B (zh) | Pet时间实时校正及重建方法和系统 | |
US8022368B2 (en) | Hybrid method for randoms variance reduction | |
JP7232256B2 (ja) | 陽電子放出断層撮影(pet)の散乱補正 | |
US11096634B2 (en) | Scatter correction based on energy response | |
Nuyts | Nuclear medicine technology and techniques | |
CN111839569B (zh) | Pet时间分辨率测算、重建系统和方法 |
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 |