CN113425320B - 时间校准方法、装置、电子设备及存储介质 - Google Patents
时间校准方法、装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN113425320B CN113425320B CN202110615662.XA CN202110615662A CN113425320B CN 113425320 B CN113425320 B CN 113425320B CN 202110615662 A CN202110615662 A CN 202110615662A CN 113425320 B CN113425320 B CN 113425320B
- Authority
- CN
- China
- Prior art keywords
- crystal
- coincidence
- time
- tof
- determining
- 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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明实施例提供一种时间校准方法、装置、电子设备及存储介质。本发明实施例通过根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像，基于重建图像进行TOF正投影，得到第二TOF符合数据，根据PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定PET系统中各个晶体对应的延迟时间，根据延迟时间，对PET系统中晶体采集的事件时间信息进行校准，利用放射源图像的TOF正投影获得基于放射源位置引起的符合数据，通过与实际符合数据的比较得到用于时间校准的延迟时间，提高了时间校准的准确性。本实施例降低了放射源在PET系统中的摆放要求。
Description
技术领域
本发明涉及医学图像处理技术领域，尤其涉及一种时间校准方法、装置、电子设备及存储介质。
背景技术
PET(Positron Emission Computed Tomography，正电子发射型计算机断层显像)设备是当今最先进的大型医疗诊断成像设备之一。PET除显示形态结构外，能够利用活体代谢在分子水平上提供有关脏器及其病变的功能信息，在诊断肿瘤、心血管和神经系统等疾病中具有卓越性能。
PET的工作原理是：在受检对象中注射含有放射性核素的药物，放射性核素发生衰变产生正电子，正电子与周围的负电子湮灭产生一对背靠背的伽马光子，PET系统通过探测光子对重建出发射正电子的核素分布图。
当前，PET系统可以根据接收到的光子对的时间差估算出正电子发生的大体位置，然后用于图像重建。因而保证获得的光子对的时间差准确性，对PET来说十分重要。但是，在PET系统中，由于探测装置的性能不一致性等因素，导致各个晶体位置上的时间判定精度不同，因此需要进行时间校准。
发明内容
为克服相关技术中存在的问题，本发明提供了一种时间校准方法、装置、电子设备及存储介质，提高时间校准的准确性。
根据本发明实施例的第一方面，提供一种时间校准方法，包括：
根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的；
基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准。
根据本发明实施例的第二方面，提供一种时间校准装置，包括：
重建模块，用于根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的；
投影模块，用于基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
确定模块，用于根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
校准模块，用于根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准。
根据本发明实施例的第三方面，提供一种电子设备，包括：内部总线，以及通过内部总线连接的存储器、处理器和外部接口；
所述存储器，用于存储时间校准逻辑对应的机器可读指令；
所述处理器，用于读取并执行所述存储器上的所述机器可读指令，并执行所述指令以实现如下操作：
根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的；
基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准。
根据本发明实施例的第四方面，提供一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现如下操作：
根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的；
基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准。
本发明实施例提供的技术方案可以包括以下有益效果：
本发明实施例，通过根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像，其中，第一TOF符合数据是对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的，基于所述重建图像进行TOF正投影，得到第二TOF符合数据，根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间，其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的，根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准，利用放射源图像的TOF正投影获得基于放射源位置引起的符合数据，通过与实际符合数据的比较得到用于时间校准的延迟时间，提高了时间校准的准确性。
本实施例中，不要求放射源中心在PET系统中心，并且也不要求放射源的中心轴与探测器中心轴平行，对于放射源在PET系统中的摆放要求低，易于在实际场景中使用，实用性强。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本说明书。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本说明书的实施例，并与说明书一起用于解释本说明书的原理。
图1是本发明实施例提供的时间校准方法的流程示例图。
图2是晶体对符合事件时间差累计直方图的示例之一。
图3是晶体对符合事件时间差累计直方图的示例之二。
图4是关联扇形面的示意图。
图5是本发明实施例提供的时间校准装置的功能方块图。
图6是本发明实施例提供的电子设备的一个硬件结构图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本发明实施例的一些方面相一致的装置和方法的例子。
在本发明实施例使用的术语是仅仅出于描述特定本发明实施例的目的，而非旨在限制本发明实施例。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解，尽管在本发明实施例可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本发明实施例范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
首先，对本发明实施例中将要提及的几个概念进行说明。
(一)TOF(Time of Flight，飞行时间)符合数据
如前所述，正电子与负电子湮灭时会产生一对背靠背的伽马光子，这对光子沿着相反的方向飞行，分别被探测器上的晶体探测到。探测到这对光子的两个晶体组成一个晶体对。该晶体对的两个晶体之间的连线称为响应线，湮灭位置位于响应线上。利用这对光子被晶体探测到的时间的差值，可以计算出湮灭位置。
其中，光子被晶体探测到为一个单事件，光子被晶体探测到的时间(也即光子到达晶体的时间)称为事件时间信息。
PET扫描数据包括PET探测器上的多个晶体采集的事件时间信息。其中，事件时间信息符合预设条件(通常为两个单事件中光子被晶体探测到的时间之差小于设置的时间阈值)的两个单事件形成一个符合事件，符合事件对应的数据为符合数据。PET图像重建使用的是PET扫描数据中的符合数据。
用TOF表示的符合数据即为TOF符合数据。
TOF符合数据是TOF符合事件的数据集，TOF符合事件包含符合事件对应的晶体对位置信息和到达两个晶体的时间差信息。例如，一个TOF符合事件记录为(A，B，ΔT)，其中A表示符合事件对应的第一个晶体，B表示符合事件对应的第二个晶体，ΔT表示伽马光子到达晶体A和晶体B的时间tA和tB的时间差，可以限定是ΔT＝tB-tA，也可以限定是ΔT＝tA-tB，使用时按照统一的一种方式限定即可，本文中统一按照ΔT＝tB-tA来使用，即给出的一对晶体，它们的时间差都是用第二个晶体到达时间减去第一个晶体到达时间来表示。
(二)放射源
放射源是一个充满了含有放射性核素的药物的物体，其中，放射性核素在放射源中是均匀分布的。放射源的形状是轴对称的，例如可以是棒源、实心圆柱源或圆筒状源等。
放射源可以是实心的，也可以是非实心的。
其次，对本发明实施例的目的进行说明。
PET图像(即PET重建图像)是一种重要的医疗辅助资料，可以提供有关器官或部位的图像信息，以辅助医生进行诊断。图像质量直接影响着PET图像提供的信息的准确性，从而影响着医生的诊断。因此，获得高质量的PET图像是重要的目标。
由前述内容可知，TOF符合数据可以用于进行PET图像重建，因此，TOF符合数据的准确性是影响PET图像质量的非常重要的因素。而TOF符合数据的准确性依赖于每个晶体采集的事件时间信息的准确性，因此，PET系统中晶体采集的事件时间信息的准确性对PET图像质量有重要的影响。
对于晶体采集的事件时间信息进行校准能够提高事件时间信息的准确性。但是，校准方案需要是实际应用场景中容易实现的。如果校准方案的实现难度大，导致实际使用情况与校准方案的要求相差较多(例如相关技术的校准方案中要求放射源中心在PET系统中心，实际放置放射源时几乎无法达到该要求)，则校准的精度会大大降低。
本发明实施例的目的是提供一种校准方案，既能够提高校准的准确性，又容易实现，满足实际应用的需要。
下面通过实施例对本发明的时间校准方法进行详细说明。
图1是本发明实施例提供的时间校准方法的流程示例图。本实施例提供的时间校准方法中，放射源中心可以不在PET系统中心，并且放射源的中心轴可以不与探测器中心轴平行，对于放射源在PET系统中的摆放要求低，易于在实际场景中使用，实用性强。
如图1所示，本实施例中，时间校准方法可以包括：
S101，根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是从对放置于PET系统中的所述放射源进行PET扫描得到的扫描数据中获取到的。
S102，基于所述重建图像进行TOF正投影，得到第二TOF符合数据。
S103，根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的。
S104，根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准。
本实施例中，第一TOF符合数据是从放射源的PET扫描数据中获取到的。例如，可以根据符合事件的时间要求(两个光子被晶体探测到的时间差值小于预设的阈值)从放射源的PET扫描数据中筛选出符合事件，根据符合事件得到TOF符合数据。
第一TOF符合数据是放射源对应的实际符合数据，该符合数据受到了探测装置的性能不一致的影响，与真实的符合数据之间存在差异，这种差异影响了晶体采集的事件时间信息的准确性。
其中，探测装置是指由晶体、光电转换器、放大电路等组成的组合电路。由于不同探测装置的各器件(晶体、光电转换器、放大电路等)之间存在差异，并且电路中的连线也存在差异，因此造成了不同探测装置的性能不一致。
放射源的重建图像是一幅PET图像，通过放射源的重建图像可以确定放射源的位置。
本实施例中，晶体对(A，B)之间的第一符合时间差，记为TOFAB，这里的晶体对是有先后顺序的，如果是晶体对(B，A)，则(B，A)之间的第一符合时间差为TOFBA＝-TOFAB。
在一个示例中，步骤S101中，根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像，可以包括：
直接使用第一TOF符合数据进行图像重建，得到放射源的重建图像；或者，将第一TOF符合数据转换为非TOF数据，使用所述非TOF数据进行图像重建，得到放射源的重建图像；其中，所述非TOF数据为不包含TOF信息的数据。
本实施例中，TOF信息指TOF符合数据中的时间差。在一个示例中，可以通过去掉TOF符合数据中的时间差信息，来获得非TOF数据。例如，TOF符合数据为(A，B，ΔT)，则对应的非TOF数据为(A，B)。
利用TOF符合数据重建的图像信噪比更高，准确性更好，但是前提是时间信息要准确，如果存在时间偏差，重建图像质量会变差。
利用非TOF数据重建的图像信噪比会低一些，但是可以正常获得放射源的轮廓，基本不受时间偏差影响，因此稳定性更好。
在应用中，可以根据实际情况选择用TOF符合数据重建还是用非TOF数据重建。即如果时间偏差被校正好，用TOF符合数据重建结果更好，如果存在时间偏差，用非TOF数据重建结果通常更好。
步骤S102中，通过TOF正投影，可以获得重建图像所对应的真实放射源应该对应的TOF符合数据，即第二TOF符合数据。第二TOF符合数据非常接近于真实的符合数据，本方案中用第二TOF符合数据代替真实的符合数据。
在一个示例中，步骤S102中，基于所述重建图像进行TOF正投影，得到第二TOF符合数据，可以包括：
根据所述重建图像设置放射源的图像值，得到放射源图像；
对所述放射源图像进行TOF正投影，得到第二TOF符合数据。
重建图像中的放射源不一定准确，均匀性可能不是特别好，这样，直接使用重建图像中的放射源进行TOF正投影，得到的TOF符合数据与真实的符合数据之间的差异会大一些，如此获得的用于校准时间的延迟时间准确性会降低。
由于已知放射源中核素分布是均匀的，因此本实施例据此将放射源设置成均匀的图像，这样获得的放射源图像更准确，从而对放射源图像进行TOF正投影能够得到更加准确的、与真实的符合数据更加接近的TOF符合数据，从而能够获得更准确的延迟时间，提高校准的准确性。
其中，放射源的图像值可以设置为预先指定的值。
在一个示例中，根据所述重建图像设置放射源的图像值，得到放射源图像，可以包括：
从所述重建图像中提取所述放射源的轮廓；
将所述轮廓内的像素值设置为指定值，得到放射源图像。
其中，可以采用相关技术从重建图像中提取放射源的轮廓，此处不再赘述。
其中，指定值可以为用户预设的常数值，例如1。
在一个示例中，根据所述重建图像设置放射源的图像值，得到放射源图像，可以包括：
根据所述重建图像，确定所述重建图像的中心位置；
根据所述中心位置和所述放射源的实际尺寸，确定所述重建图像中所述放射源所在的区域；
将所述放射源所在的区域内的像素值设置为指定值，得到放射源图像。
本实施例中，可以利用质心法，对每一层图像，用如下的公式(1)计算重建图像的中心位置。
公式(1)中，Ii表示本层图像中i号像素的像素值，xi表示i号像素的x坐标，yi表示i号像素的y坐标。
TOF符合数据中包含符合时间差的信息，因此可以从TOF符合数据中直接提取出符合时间差。例如，TOF符合数据为(A，B，ΔT)，则对应的符合时间差为ΔT。
通过由第一TOF符合数据确定的第一符合时间差和由第二TOF符合数据确定的第二符合时间差之间的差异，可以确定晶体实际采集的事件时间信息与真实的事件时间信息之间的差异，即延迟时间，进而可以利用延迟时间校准晶体实际采集的事件时间信息，以提高晶体采集的事件时间信息的准确性。
其中，符合时间差是基于符合数据确定的。
在一个示例中，基于符合数据确定符合时间差的方法可以包括：
若所述放射源为实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的峰值位置；
将所述峰值位置对应的符合时间差确定为所述第一晶体对对应的符合时间差。
当放射源为实心放射源时，对应的符合时间差的计数有一个峰值。
一对晶体上采集的TOF符合数据，以时间差为坐标建立时间差累计直方图(表示相同时间差上获得的计数累加)，如图2所示，图2是晶体对符合事件时间差累计直方图的示例之一。
然后，根据时间差累计直方图可以计算峰值位置。例如，在一个示例中，可以通过高斯函数拟合确定峰值位置。在另一个示例中，可以利用质心法t＝∑itiCi/∑iCi确定峰值位置。
在一个示例中，基于符合数据确定符合时间差的方法可以包括：
若所述放射源为非实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的第一峰值位置和第二峰值位置；
将所述第一峰值位置对应的符合时间差和所述第二峰值位置对应的符合时间差的平均值确定为所述第一晶体对对应的符合时间差。
当放射源为非实心放射源时，对应的符合时间差的计数有两个峰值。
图3是晶体对符合事件时间差累计直方图的示例之二。如图3所示，非实心放射源，如圆筒状放射源，其时间差累积直方图是双峰直方图。对于这种直方图计算峰值位置，可以利用质心法t＝∑itiCi/∑iCi确定峰值位置，也可以通过高斯函数拟合分别确定两个峰的峰值位置，然后再用两个峰值位置的平均值作为直方图的峰值位置。
在一个示例中，根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间，可以包括：
针对所述PET系统中的第一晶体，确定所述第一晶体对应的关联扇形面上的多个第二晶体；所述第一晶体为所述PET系统中的任一晶体；
根据各目标晶体对对应的第一符合时间差和第二符合时间差，确定各目标晶体对对应的初始时间偏差校正量；每个目标晶体对对应所述第一晶体和所述多个第二晶体中的一个晶体；
根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间。
图4是关联扇形面的示意图。如图4所示，对于晶体A来说，晶体A对应的关联扇形面上有多个晶体B，则晶体A和每个晶体B均组成一个目标晶体对。
其中，初始时间偏差校正量用表示，/>其中，TOFAB为晶体对(A，B)的第一符合时间差，ΔAB为(A，B)的第二符合时间差。
在一个示例中，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，可以包括：
获取所有目标晶体对对应的初始时间偏差校正量的平均值；
根据所述平均值，确定所述第一晶体对应的延迟时间。
例如，第一晶体为晶体A，本实施例中，晶体A对应的延迟时间可以用如下的公式(2)表示：
公式(2)中，为晶体A对应的延迟时间，n为晶体A对应的关联扇形面上的晶体个数，FanA是晶体A的关联扇形面。
在一个示例中，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，可以包括：
设置所述第一晶体对应的延迟时间的初始值为0，以及设置时间偏差校正量的初始值为所述初始时间偏差校正量；执行如下操作：
根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量，确定所述第一晶体对应的上一次的时间偏差校正量；
根据上一次迭代后获得的延迟时间和所述第一晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的延迟时间；
检测当前是否满足预设条件，若满足，则停止迭代过程，将本次迭代后获得的延迟时间作为所述第一晶体对应的延迟时间；否则，根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量、所述第一晶体对应的上一次的时间偏差校正量、所述第二晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的各目标晶体对对应的时间偏差校正量。
其中，预设条件可以是当前迭代次数达到预设的最大迭代次数，或者预设条件可以是当前各目标晶体对对应的时间偏差校正量均小于预设的校正量阈值。
例如，第一晶体为晶体A，本实施例中，设设则晶体A对应的延迟时间可以用如下的公式(3)表示：
根据公式(3)迭代执行，直到当前迭代次数达到预设的最大迭代次数，或者都小于预设的校正量阈值，则停止迭代。
在一个示例中，根据所述延迟时间，对所述PET系统中相应晶体采集的事件时间信息进行校准，可以包括：
针对所述PET系统中的第一晶体，从延迟时间表中读取所述第一晶体对应的延迟时间；所述延迟时间表中存储有所述PET系统中所有晶体对应的延迟时间；
将所述第一晶体采集的事件时间信息延迟所述第一晶体对应的延迟时间，得到校准后的事件时间信息。
在实际应用中，可以根据获得的每个晶体的延迟时间τA，建立延迟时间表，表中值为τA，利用τA校准采集数据，即对晶体A采集的事件时间信息上延迟τA。
在实际应用中，可以通过执行一次图1所示的方法，获得PET系统中各晶体对应的延迟时间，然后建立延迟时间表，后续PET系统扫描过程中采集的数据可以直接使用延迟时间表中的相应延迟时间进行校准。
本发明实施例提供的时间校准方法，通过根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像，其中，第一TOF符合数据是从对放置于PET系统中的所述放射源进行PET扫描得到的扫描数据中获取到的，基于所述重建图像进行TOF正投影，得到第二TOF符合数据，根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间，其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的，根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准，利用放射源图像的TOF正投影获得基于放射源位置引起的符合数据，通过与实际符合数据的比较得到用于时间校准的延迟时间，提高了时间校准的准确性。
并且，本发明实施例提供的时间校准方法的方案，允许放射源中心不在PET系统中心，并且允许放射源的中心轴与探测器中心轴不相平行，显著降低了放射源在PET系统中的摆放要求，从而降低了实现难度，易于在实际场景中使用，实用性强。
基于上述的方法实施例，本发明实施例还提供了相应的装置、设备及存储介质实施例。
图5是本发明实施例提供的时间校准装置的功能方块图。如图5所示，本实施例中，时间校准装置可以包括：
重建模块510，用于根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是从对放置于PET系统中的所述放射源进行PET扫描得到的扫描数据中获取到的；
投影模块520，用于基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
确定模块530，用于根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
校准模块540，用于根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准。
在一个示例中，确定模块530可以具体用于：
针对所述PET系统中的第一晶体，确定所述第一晶体对应的关联扇形面上的多个第二晶体；所述第一晶体为所述PET系统中的任一晶体；
根据各目标晶体对对应的第一符合时间差和第二符合时间差，确定各目标晶体对对应的初始时间偏差校正量；每个目标晶体对对应所述第一晶体和所述多个第二晶体中的一个晶体；
根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间。
在一个示例中，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
获取所有目标晶体对对应的初始时间偏差校正量的平均值；
根据所述平均值，确定所述第一晶体对应的延迟时间。
在一个示例中，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
设置所述第一晶体对应的延迟时间的初始值为0，以及设置时间偏差校正量的初始值为所述初始时间偏差校正量；执行如下操作：
根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量，确定所述第一晶体对应的上一次的时间偏差校正量；
根据上一次迭代后获得的延迟时间和所述第一晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的延迟时间；
检测当前是否满足预设条件，若满足，则停止迭代过程，将本次迭代后获得的延迟时间作为所述第一晶体对应的延迟时间；否则，根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量、所述第一晶体对应的上一次的时间偏差校正量、所述第二晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的各目标晶体对对应的时间偏差校正量。
在一个示例中，重建模块510可以具体用于：
直接使用第一TOF符合数据进行图像重建，得到放射源的重建图像；或者，将第一TOF符合数据转换为非TOF数据，使用所述非TOF数据进行图像重建，得到放射源的重建图像；其中，所述非TOF数据为不包含TOF信息的数据。
在一个示例中，投影模块520可以具体用于：
根据所述重建图像设置放射源的图像值，得到放射源图像；
对所述放射源图像进行TOF正投影，得到第二TOF符合数据。
在一个示例中，根据所述重建图像设置放射源的图像值，得到放射源图像，包括：
从所述重建图像中提取所述放射源的轮廓；
将所述轮廓内的像素值设置为指定值，得到放射源图像。
在一个示例中，根据所述重建图像设置放射源的图像值，得到放射源图像，包括：
根据所述重建图像，确定所述重建图像的中心位置；
根据所述中心位置和所述放射源的实际尺寸，确定所述重建图像中所述放射源所在的区域；
将所述放射源所在的区域内的像素值设置为指定值，得到放射源图像。
在一个示例中，基于符合数据确定符合时间差的方法包括：
若所述放射源为实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的峰值位置；
将所述峰值位置对应的符合时间差确定为所述第一晶体对对应的符合时间差。
在一个示例中，基于符合数据确定符合时间差的方法包括：
若所述放射源为非实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的第一峰值位置和第二峰值位置；
将所述第一峰值位置对应的符合时间差和所述第二峰值位置对应的符合时间差的平均值确定为所述第一晶体对对应的符合时间差。
在一个示例中，校准模块540可以具体用于：
针对所述PET系统中的第一晶体，从延迟时间表中读取所述第一晶体对应的延迟时间；所述延迟时间表中存储有所述PET系统中所有晶体对应的延迟时间；
将所述第一晶体采集的事件时间信息延迟所述第一晶体对应的延迟时间，得到校准后的事件时间信息。
本发明实施例还提供了一种电子设备。图6是本发明实施例提供的电子设备的一个硬件结构图。如图6所示，电子设备包括：内部总线601，以及通过内部总线连接的存储器602，处理器603和外部接口604。
所述存储器602，用于存储时间校准逻辑对应的机器可读指令；
所述处理器603，用于读取存储器602上的机器可读指令，并执行所述指令以实现如下操作：
根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是从对放置于PET系统中的所述放射源进行PET扫描得到的扫描数据中获取到的；
基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准。
在一个示例中，根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间，包括：
针对所述PET系统中的第一晶体，确定所述第一晶体对应的关联扇形面上的多个第二晶体；所述第一晶体为所述PET系统中的任一晶体；
根据各目标晶体对对应的第一符合时间差和第二符合时间差，确定各目标晶体对对应的初始时间偏差校正量；每个目标晶体对对应所述第一晶体和所述多个第二晶体中的一个晶体；
根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间。
在一个示例中，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
获取所有目标晶体对对应的初始时间偏差校正量的平均值；
根据所述平均值，确定所述第一晶体对应的延迟时间。
在一个示例中，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
设置所述第一晶体对应的延迟时间的初始值为0，以及设置时间偏差校正量的初始值为所述初始时间偏差校正量；执行如下操作：
根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量，确定所述第一晶体对应的上一次的时间偏差校正量；
根据上一次迭代后获得的延迟时间和所述第一晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的延迟时间；
检测当前是否满足预设条件，若满足，则停止迭代过程，将本次迭代后获得的延迟时间作为所述第一晶体对应的延迟时间；否则，根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量、所述第一晶体对应的上一次的时间偏差校正量、所述第二晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的各目标晶体对对应的时间偏差校正量。
在一个示例中，根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像，包括：
直接使用第一TOF符合数据进行图像重建，得到放射源的重建图像；或者，将第一TOF符合数据转换为非TOF数据，使用所述非TOF数据进行图像重建，得到放射源的重建图像；其中，所述非TOF数据为不包含TOF信息的数据。
在一个示例中，基于所述重建图像进行TOF正投影，得到第二TOF符合数据，包括：
根据所述重建图像设置放射源的图像值，得到放射源图像；
对所述放射源图像进行TOF正投影，得到第二TOF符合数据。
在一个示例中，根据所述重建图像设置放射源的图像值，得到放射源图像，包括：
从所述重建图像中提取所述放射源的轮廓；
将所述轮廓内的像素值设置为指定值，得到放射源图像。
在一个示例中，根据所述重建图像设置放射源的图像值，得到放射源图像，包括：
根据所述重建图像，确定所述重建图像的中心位置；
根据所述中心位置和所述放射源的实际尺寸，确定所述重建图像中所述放射源所在的区域；
将所述放射源所在的区域内的像素值设置为指定值，得到放射源图像。
在一个示例中，基于符合数据确定符合时间差的方法包括：
若所述放射源为实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的峰值位置；
将所述峰值位置对应的符合时间差确定为所述第一晶体对对应的符合时间差。
在一个示例中，基于符合数据确定符合时间差的方法包括：
若所述放射源为非实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的第一峰值位置和第二峰值位置；
将所述第一峰值位置对应的符合时间差和所述第二峰值位置对应的符合时间差的平均值确定为所述第一晶体对对应的符合时间差。
在一个示例中，根据所述延迟时间，对所述PET系统中相应晶体采集的事件时间信息进行校准，包括：
针对所述PET系统中的第一晶体，从延迟时间表中读取所述第一晶体对应的延迟时间；所述延迟时间表中存储有所述PET系统中所有晶体对应的延迟时间；
将所述第一晶体采集的事件时间信息延迟所述第一晶体对应的延迟时间，得到校准后的事件时间信息。
其中，电子设备可以为PET系统中的控制台设备。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现如下操作：
根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是从对放置于PET系统中的所述放射源进行PET扫描得到的扫描数据中获取到的；
基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准。
在一个示例中，根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间，包括：
针对所述PET系统中的第一晶体，确定所述第一晶体对应的关联扇形面上的多个第二晶体；所述第一晶体为所述PET系统中的任一晶体；
根据各目标晶体对对应的第一符合时间差和第二符合时间差，确定各目标晶体对对应的初始时间偏差校正量；每个目标晶体对对应所述第一晶体和所述多个第二晶体中的一个晶体；
根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间。
在一个示例中，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
获取所有目标晶体对对应的初始时间偏差校正量的平均值；
根据所述平均值，确定所述第一晶体对应的延迟时间。
在一个示例中，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
设置所述第一晶体对应的延迟时间的初始值为0，以及设置时间偏差校正量的初始值为所述初始时间偏差校正量；执行如下操作：
根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量，确定所述第一晶体对应的上一次的时间偏差校正量；
根据上一次迭代后获得的延迟时间和所述第一晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的延迟时间；
检测当前是否满足预设条件，若满足，则停止迭代过程，将本次迭代后获得的延迟时间作为所述第一晶体对应的延迟时间；否则，根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量、所述第一晶体对应的上一次的时间偏差校正量、所述第二晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的各目标晶体对对应的时间偏差校正量。
在一个示例中，根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像，包括：
直接使用第一TOF符合数据进行图像重建，得到放射源的重建图像；或者，将第一TOF符合数据转换为非TOF数据，使用所述非TOF数据进行图像重建，得到放射源的重建图像；其中，所述非TOF数据为不包含TOF信息的数据。
在一个示例中，基于所述重建图像进行TOF正投影，得到第二TOF符合数据，包括：
根据所述重建图像设置放射源的图像值，得到放射源图像；
对所述放射源图像进行TOF正投影，得到第二TOF符合数据。
在一个示例中，根据所述重建图像设置放射源的图像值，得到放射源图像，包括：
从所述重建图像中提取所述放射源的轮廓；
将所述轮廓内的像素值设置为指定值，得到放射源图像。
在一个示例中，根据所述重建图像设置放射源的图像值，得到放射源图像，包括：
根据所述重建图像，确定所述重建图像的中心位置；
根据所述中心位置和所述放射源的实际尺寸，确定所述重建图像中所述放射源所在的区域；
将所述放射源所在的区域内的像素值设置为指定值，得到放射源图像。
在一个示例中，基于符合数据确定符合时间差的方法包括：
若所述放射源为实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的峰值位置；
将所述峰值位置对应的符合时间差确定为所述第一晶体对对应的符合时间差。
在一个示例中，基于符合数据确定符合时间差的方法包括：
若所述放射源为非实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的第一峰值位置和第二峰值位置；
将所述第一峰值位置对应的符合时间差和所述第二峰值位置对应的符合时间差的平均值确定为所述第一晶体对对应的符合时间差。
在一个示例中，根据所述延迟时间，对所述PET系统中相应晶体采集的事件时间信息进行校准，包括：
针对所述PET系统中的第一晶体，从延迟时间表中读取所述第一晶体对应的延迟时间；所述延迟时间表中存储有所述PET系统中所有晶体对应的延迟时间；
将所述第一晶体采集的事件时间信息延迟所述第一晶体对应的延迟时间，得到校准后的事件时间信息。
对于装置和设备实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的，作为模块显示的部件可以是或者也可以不是物理模块，即可以位于一个地方，或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本说明书方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本领域技术人员在考虑说明书及实践这里申请的发明后，将容易想到本说明书的其它实施方案。本说明书旨在涵盖本说明书的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本说明书的一般性原理并包括本说明书未申请的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本说明书的真正范围和精神由下面的权利要求指出。
应当理解的是，本说明书并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本说明书的范围仅由所附的权利要求来限制。
以上所述仅为本说明书的较佳实施例而已，并不用以限制本说明书，凡在本说明书的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本说明书保护的范围之内。
Claims (12)
1.一种时间校准方法，其特征在于，包括：
根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是从对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的扫描数据中获取到的；
基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准；
所述根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间，包括：
针对所述PET系统中的第一晶体，确定所述第一晶体对应的关联扇形面上的多个第二晶体；所述第一晶体为所述PET系统中的任一晶体；
根据各目标晶体对对应的第一符合时间差和第二符合时间差，确定各目标晶体对对应的初始时间偏差校正量；每个目标晶体对对应所述第一晶体和所述多个第二晶体中的一个晶体；
根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间；
所述根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
获取所有目标晶体对对应的初始时间偏差校正量的平均值；
根据所述平均值，确定所述第一晶体对应的延迟时间。
2.根据权利要求1所述的方法，其特征在于，根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
设置所述第一晶体对应的延迟时间的初始值为0，以及设置时间偏差校正量的初始值为所述初始时间偏差校正量；执行如下操作：
根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量，确定所述第一晶体对应的上一次的时间偏差校正量；
根据上一次迭代后获得的延迟时间和所述第一晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的延迟时间；
检测当前是否满足预设条件，若满足，则停止迭代过程，将本次迭代后获得的延迟时间作为所述第一晶体对应的延迟时间；否则，根据上一次迭代后获得的各目标晶体对对应的时间偏差校正量、所述第一晶体对应的上一次的时间偏差校正量、所述第二晶体对应的上一次的时间偏差校正量，确定本次迭代后获得的各目标晶体对对应的时间偏差校正量。
3.根据权利要求1所述的方法，其特征在于，根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像，包括：
直接使用第一TOF符合数据进行图像重建，得到放射源的重建图像；或者，将第一TOF符合数据转换为非TOF数据，使用所述非TOF数据进行图像重建，得到放射源的重建图像；其中，所述非TOF数据为不包含TOF信息的数据。
4.根据权利要求1所述的方法，其特征在于，基于所述重建图像进行TOF正投影，得到第二TOF符合数据，包括：
根据所述重建图像设置放射源的图像值，得到放射源图像；
对所述放射源图像进行TOF正投影，得到第二TOF符合数据。
5.根据权利要求4所述的方法，其特征在于，根据所述重建图像设置放射源的图像值，得到放射源图像，包括：
从所述重建图像中提取所述放射源的轮廓；
将所述轮廓内的像素值设置为指定值，得到放射源图像。
6.根据权利要求4所述的方法，其特征在于，根据所述重建图像设置放射源的图像值，得到放射源图像，包括：
根据所述重建图像，确定所述重建图像的中心位置；
根据所述中心位置和所述放射源的实际尺寸，确定所述重建图像中所述放射源所在的区域；
将所述放射源所在的区域内的像素值设置为指定值，得到放射源图像。
7.根据权利要求1所述的方法，其特征在于，基于符合数据确定符合时间差的方法包括：
若所述放射源为实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的峰值位置；
将所述峰值位置对应的符合时间差确定为所述第一晶体对对应的符合时间差。
8.根据权利要求1所述的方法，其特征在于，基于符合数据确定符合时间差的方法包括：
若所述放射源为非实心放射源，针对所述PET系统中的第一晶体对，统计各符合时间差对应的计数；所述第一晶体对为所述PET系统中的任一晶体对；
查找计数值的第一峰值位置和第二峰值位置；
将所述第一峰值位置对应的符合时间差和所述第二峰值位置对应的符合时间差的平均值确定为所述第一晶体对对应的符合时间差。
9.根据权利要求1所述的方法，其特征在于，根据所述延迟时间，对所述PET系统中相应晶体采集的事件时间信息进行校准，包括：
针对所述PET系统中的第一晶体，从延迟时间表中读取所述第一晶体对应的延迟时间；所述延迟时间表中存储有所述PET系统中所有晶体对应的延迟时间；
将所述第一晶体采集的事件时间信息延迟所述第一晶体对应的延迟时间，得到校准后的事件时间信息。
10.一种时间校准装置，其特征在于，包括：
重建模块，用于根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是从对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的扫描数据中获取到的；
投影模块，用于基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
确定模块，用于根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
校准模块，用于根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准；
所述确定模块，在用于根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间时，包括：
针对所述PET系统中的第一晶体，确定所述第一晶体对应的关联扇形面上的多个第二晶体；所述第一晶体为所述PET系统中的任一晶体；
根据各目标晶体对对应的第一符合时间差和第二符合时间差，确定各目标晶体对对应的初始时间偏差校正量；每个目标晶体对对应所述第一晶体和所述多个第二晶体中的一个晶体；
根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间；
所述确定模块，在用于根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间时，包括：
获取所有目标晶体对对应的初始时间偏差校正量的平均值；
根据所述平均值，确定所述第一晶体对应的延迟时间。
11.一种电子设备，其特征在于，包括：内部总线，以及通过内部总线连接的存储器、处理器和外部接口；
所述存储器，用于存储时间校准逻辑对应的机器可读指令；
所述处理器，用于读取并执行所述存储器上的所述机器可读指令，并执行所述指令以实现如下操作：
根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的；
基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准；
所述根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间，包括：
针对所述PET系统中的第一晶体，确定所述第一晶体对应的关联扇形面上的多个第二晶体；所述第一晶体为所述PET系统中的任一晶体；
根据各目标晶体对对应的第一符合时间差和第二符合时间差，确定各目标晶体对对应的初始时间偏差校正量；每个目标晶体对对应所述第一晶体和所述多个第二晶体中的一个晶体；
根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间；
所述根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
获取所有目标晶体对对应的初始时间偏差校正量的平均值；
根据所述平均值，确定所述第一晶体对应的延迟时间。
12.一种计算机可读存储介质，其特征在于，其上存储有计算机程序，其中，所述程序被处理器执行时实现如下操作：
根据第一飞行时间TOF符合数据进行图像重建，得到放射源的重建图像；其中，第一TOF符合数据是从对放置于正电子发射型计算机断层显像PET系统中的所述放射源进行PET扫描得到的扫描数据中获取到的；
基于所述重建图像进行TOF正投影，得到第二TOF符合数据；
根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间；其中，所述第一符合时间差是基于所述第一TOF符合数据确定的，所述第二符合时间差是基于所述第二TOF符合数据确定的；
根据所述延迟时间，对所述PET系统中晶体采集的事件时间信息进行校准；
所述根据所述PET系统中各晶体对对应的第一符合时间差和第二符合时间差，确定所述PET系统中各个晶体对应的延迟时间，包括：
针对所述PET系统中的第一晶体，确定所述第一晶体对应的关联扇形面上的多个第二晶体；所述第一晶体为所述PET系统中的任一晶体；
根据各目标晶体对对应的第一符合时间差和第二符合时间差，确定各目标晶体对对应的初始时间偏差校正量；每个目标晶体对对应所述第一晶体和所述多个第二晶体中的一个晶体；
根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间；
所述根据各目标晶体对对应的初始时间偏差校正量，确定所述第一晶体对应的延迟时间，包括：
获取所有目标晶体对对应的初始时间偏差校正量的平均值；
根据所述平均值，确定所述第一晶体对应的延迟时间。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110615662.XA CN113425320B (zh) | 2021-06-02 | 2021-06-02 | 时间校准方法、装置、电子设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110615662.XA CN113425320B (zh) | 2021-06-02 | 2021-06-02 | 时间校准方法、装置、电子设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113425320A CN113425320A (zh) | 2021-09-24 |
| CN113425320B true CN113425320B (zh) | 2023-09-26 |
Family
ID=77803612
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110615662.XA Active CN113425320B (zh) | 2021-06-02 | 2021-06-02 | 时间校准方法、装置、电子设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113425320B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105496436A (zh) * | 2015-11-28 | 2016-04-20 | 上海联影医疗科技有限公司 | 用于pet装置的时间校正方法和装置 |
| CN106353786A (zh) * | 2016-09-30 | 2017-01-25 | 上海联影医疗科技有限公司 | 正电子发射断层成像系统飞行时间性能检测方法及装置 |
| CN107137101A (zh) * | 2017-04-24 | 2017-09-08 | 沈阳东软医疗系统有限公司 | 一种时间校准方法和装置 |
| WO2020059135A1 (ja) * | 2018-09-21 | 2020-03-26 | 株式会社島津製作所 | データ処理方法、プログラム、データ処理装置および陽電子放出断層撮像装置 |
| CN110988974A (zh) * | 2019-12-12 | 2020-04-10 | 明峰医疗系统股份有限公司 | 曲面辐射源的pet探测器飞行时间的校准系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7820975B2 (en) * | 2006-01-03 | 2010-10-26 | Koninklijke Philips Electronics N.V. | Achieving accurate time-of-flight calibrations with a stationary coincidence point source |
| US8265365B2 (en) * | 2010-09-20 | 2012-09-11 | Siemens Medical Solutions Usa, Inc. | Time of flight scatter distribution estimation in positron emission tomography |
| US10663608B2 (en) * | 2015-09-21 | 2020-05-26 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for calibrating a PET scanner |
| US9990741B2 (en) * | 2015-09-28 | 2018-06-05 | Siemens Medical Solutions Usa, Inc. | Motion correction in a projection domain in time of flight positron emission tomography |
-
2021
- 2021-06-02 CN CN202110615662.XA patent/CN113425320B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105496436A (zh) * | 2015-11-28 | 2016-04-20 | 上海联影医疗科技有限公司 | 用于pet装置的时间校正方法和装置 |
| CN106353786A (zh) * | 2016-09-30 | 2017-01-25 | 上海联影医疗科技有限公司 | 正电子发射断层成像系统飞行时间性能检测方法及装置 |
| CN107137101A (zh) * | 2017-04-24 | 2017-09-08 | 沈阳东软医疗系统有限公司 | 一种时间校准方法和装置 |
| WO2020059135A1 (ja) * | 2018-09-21 | 2020-03-26 | 株式会社島津製作所 | データ処理方法、プログラム、データ処理装置および陽電子放出断層撮像装置 |
| CN110988974A (zh) * | 2019-12-12 | 2020-04-10 | 明峰医疗系统股份有限公司 | 曲面辐射源的pet探测器飞行时间的校准系统 |
Non-Patent Citations (1)
| Title |
|---|
| 3D PET数据校正常用方法研究;俞王新;王之光;谢舒平;;医疗卫生装备(10);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113425320A (zh) | 2021-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108542412B (zh) | 图像重建方法和pet设备 | |
| US10319118B2 (en) | Reconstructing PET image | |
| US8089043B2 (en) | Systems and methods for calibrating time alignment for a positron emission tomography (PET) system | |
| CN110853082B (zh) | 医学图像配准方法、装置、电子设备及计算机存储介质 | |
| EP3229689B1 (en) | Outside-fov activity estimation using surview and prior patient data in positron emission tomography | |
| US8330112B2 (en) | Timing response improvement in light-sharing detectors | |
| US10215864B2 (en) | System and method to improve image quality of emission tomography when using advanced radionuclides | |
| JP6188418B2 (ja) | 偶発イベント削減方法、偶発イベント削減装置及び非一時的コンピュータ可読記憶媒体 | |
| CN109887048B (zh) | Pet散射校正方法、图像重建方法、装置及电子设备 | |
| WO2001093763A2 (en) | Event localization and fall-off correction by distance-dependent weighting | |
| EP0862749B1 (en) | Real time pet imaging method for a single photon gamma camera | |
| CN110749916A (zh) | 获取pet探测器晶体时间延迟量的方法、装置和计算机设备 | |
| CN111568452A (zh) | Pet系统状态检测方法、装置、计算机设备和存储介质 | |
| CN110680370B (zh) | 图像重建方法、装置、控制台设备及pet系统 | |
| CN108209958B (zh) | 一种归一化校正因子的确定、获取方法及医学成像方法 | |
| US6294788B1 (en) | Randoms correction in positron imaging | |
| WO2019145398A1 (en) | A dead-time correction method in quantitative positron emission tomography (pet) reconstruction for various objects and radioactivity distributions | |
| CN113425320B (zh) | 时间校准方法、装置、电子设备及存储介质 | |
| US11231508B2 (en) | Gamma camera dead time determination in real time using long lived radioisotopes | |
| CN107961028B (zh) | 一种归一化校正因子获取、确定方法及医学成像方法 | |
| CN110687585A (zh) | 获取晶体效率的方法、装置、计算机设备和存储介质 | |
| CN108805947B (zh) | Pet数据处理方法和设备、pet成像系统 | |
| CN111904448B (zh) | 数据处理方法、装置、pet设备及pet系统 | |
| CN108932740B (zh) | 一种归一化校正因子获取方法及医学成像方法 | |
| WO2021115917A1 (en) | Machine learning-based scintillator response modelling for increased spatial resolution in nuclear imaging |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20230412 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. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |