CN110680370A - 图像重建方法、装置、控制台设备及pet*** - Google Patents
图像重建方法、装置、控制台设备及pet*** Download PDFInfo
- Publication number
- CN110680370A CN110680370A CN201910945559.4A CN201910945559A CN110680370A CN 110680370 A CN110680370 A CN 110680370A CN 201910945559 A CN201910945559 A CN 201910945559A CN 110680370 A CN110680370 A CN 110680370A
- Authority
- CN
- China
- Prior art keywords
- data
- event
- scanning
- scanning data
- crystal
- 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 44
- 239000013078 crystal Substances 0.000 claims abstract description 169
- 238000012937 correction Methods 0.000 claims description 40
- 230000004044 response Effects 0.000 claims description 38
- 238000012545 processing Methods 0.000 claims description 21
- 230000006870 function Effects 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 10
- 230000010354 integration Effects 0.000 claims description 7
- 238000011084 recovery Methods 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 abstract description 7
- 238000002600 positron emission tomography Methods 0.000 description 81
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012879 PET imaging Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 238000012636 positron electron tomography Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Images
Classifications
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4233—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using matrix detectors
-
- 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
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- 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
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- 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
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
-
- 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/54—Control of apparatus or devices for radiation diagnosis
-
- 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
-
- 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
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明实施例提供一种图像重建方法、装置、控制台设备及PET***。本发明实施例,通过从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据,对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像,通过将晶体内康普顿散射事件恢复为符合事件,增加了感兴趣区域的符合事件计数率,提高了***灵敏度。
Description
技术领域
本发明涉及医学图像处理技术领域,尤其涉及一种图像重建方法、装置、控制台设备及PET***。
背景技术
正电子发射断层成像(Positron Emission Tomography,PET)技术是当前的热门医学影像技术之一。基于飞行时间(Time of Flight,TOF)技术的正电子发射断层成像(TOF-PET)是核医学领域技术发展的重要方向,与传统的PET成像相比,它可以提供更加清晰的图像。
PET探测器用于在PET***扫描过程中,探测被检体内发出的高能光子,并基于探测结果生成PET扫描数据。PET扫描数据中不仅包括符合事件对应的单事件数据,也包括晶体内康普顿散射事件对应的单事件数据。在图像重建过程中,相关技术根据能量大小将晶体内康普顿散射事件对应的单事件数据过滤掉,利用符合事件对应的单事件数据进行图像重建,这样使得***总符合计数率减少,导致***灵敏度降低。
发明内容
为克服相关技术中存在的问题,本发明提供了一种图像重建方法、装置、控制台设备及PET***,提高***灵敏度。
根据本发明实施例的第一方面,提供一种图像重建方法,所述方法包括:
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述PET扫描数据中符合事件对应的扫描数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
根据本发明实施例的第二方面,提供一种图像重建装置,所述装置包括:
识别模块,用于从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
匹配模块,用于对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
散射恢复模块,用于根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
重建模块,用于根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
根据本发明实施例的第三方面,提供一种控制台设备,包括内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于连接PET***的探测器,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;
所述存储器,用于存储校正逻辑对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行如下操作:
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
根据本发明实施例的第四方面,提供一种PET***,包括探测器、扫描床和控制台设备,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;其中,
所述闪烁晶体,用于在所述PET***扫描过程中,探测被检体内发出的高能光子,并将所述高能光子转换为光信号;
所述光电转换器件,用于将所述光信号转换成电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述控制台设备,用于从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
根据本发明实施例的第五方面,提供一种计算机可读存储介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现如下操作:
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
本发明实施例提供的技术方案可以包括以下有益效果:
本发明实施例,通过从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据,对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,根据第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像,将晶体内康普顿散射事件恢复为符合事件,增加了散射区域的能量值,提高了***灵敏度。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本说明书。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本说明书的实施例,并与说明书一起用于解释本说明书的原理。
图1是本发明实施例提供的图像重建方法的一个流程示例图。
图2是本发明实施例提供的晶体内康普顿散射事件形成的多条LOR响应线的示例图。
图3是本发明实施例提供的图像重建装置的功能方块图。
图4是本发明实施例提供的控制台设备的一个硬件结构图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明实施例的一些方面相一致的装置和方法的例子。
在本发明实施例使用的术语是仅仅出于描述特定本发明实施例的目的,而非旨在限制本发明实施例。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本发明实施例可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本发明实施例范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
PET***所采集的PET扫描数据是由一个个的单事件数据组成的。每个单事件数据包括能量信息、时间信息和位置信息。其中,能量信息表示入射到探测器晶体中的光子的能量,时间信息表示该光子入射到探测器晶体中的时间,位置信息为该光子所入射到的晶体的位置。
在对被检体进行PET扫描的过程中,被检体中的目标区域发生湮灭事件时,会产生两个能量相同、方向相反的光子,该两个光子分别到达互成180度角的探测器晶体,被探测器捕获,生成两组单事件数据。通过该两组单事件数据进行符合,可以确定湮灭位置,从而生成一个符合事件计数。
如果该两个光子中的一个光子按照原方向(此处原方向是指湮灭事件所产生的光子的运动方向)到达探测器晶体,另一个光子在中途发生晶体内康普顿散射,则该另一个光子经散射后产生多个散射光子,该多个散射光子分别入射到探测器晶体中,被探测器捕获,生成多组单事件数据,每一个散射光子对应一个单事件。通过晶体内康普顿散射事件相关的单事件数据,可以恢复出原有的符合事件,从而增加符合事件计数。
相关技术中,在利用PET扫描数据进行图像重建的过程中,首先会根据预设的第二能量阈值对PET扫描数据进行筛选,能量值超过能量阈值的单事件数据获得通过,用于图像重建,而能量值未超过能量阈值的单事件数据被过滤掉。由于散射光子的能量小于散射前光子的能量,因此在相关技术中,与晶体内康普顿散射事件相关的单事件数据会被过滤掉,从而使得符合事件计数减少,***灵敏度降低。
本发明提供的图像重建方法,能够将晶体内康普顿散射事件恢复为符合事件,增加符合事件计数,提高***灵敏度。
下面通过实施例分别对本发明提供的图像重建方法进行说明。
图1是本发明实施例提供的图像重建方法的一个流程示例图。该图像重建方法可以应用于对PET扫描数据进行图像重建。如图1所示,图像重建方法可以包括:
S101,从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据。
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据
S102,对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据。
S103,根据所述第一扫描数据及与一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据。
S104,根据第二扫描数据中符合事件对应的数据以及由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
PET扫描数据包括PET扫描过程中的所有单事件数据。
晶体内康普顿散射是指发生在晶体内的康普顿散射。
每一个晶体内康普顿散射事件对应两个或两个以上单事件,因此步骤S101中,每个晶体内康普顿散射事件对应的第一扫描数据包括两组或两组以上的单事件数据。
步骤S101中,能量值超过预设的第一能量阈值的第二扫描数据是指上述相关技术中通过筛选的单事件数据。
举例说明。在一个湮灭事件中,产生两个能量相同、方向相反的光子:光子A和光子B。光子B按照原方向入射到探测器晶体中,被记录为单事件b;光子A发生晶体内康普顿散射产生两个散射光子:光子A1和光子A2。光子A1入射到探测器晶体中,被记录为单事件a1,光子A2入射到探测器晶体中,被记录为单事件a2。
在另一个湮灭事件中,产生两个能量相同、方向相反的光子:光子C和光子D。光子C按照原方向入射到探测器晶体中,被记录为单事件c,光子D按照原方向入射到探测器晶体中,被记录为单事件d。
则单事件a1的数据和单事件a2的数据为第一扫描数据,单事件b的数据、单事件c的数据和单事件d的数据为第二扫描数据。
在一个示例性的实现过程中,步骤S101可以包括:
对于PET扫描数据中的每个单事件数据,判断PET扫描数据中是否存在满足如下全部条件的目标单事件数据:
包含的空间位置信息所对应的第一空间位置与该单事件数据对应的第二空间位置之间的距离小于预设的距离阈值;
包含的时间信息所对应的第一时间与该单事件数据对应的第二时间之间的差值小于预设的时间阈值;
包含的第一能量值与该单事件数据中的第二能量值均小于预设的第二能量阈值,所述第二能量阈值不大于所述第一能量阈值;
若存在,则确定所述目标单事件数据和该单事件数据为同一散射事件对应的第一扫描数据,所述同一散射事件为晶体内康普顿散射事件。
本实施例利用同一晶体内康普顿散射事件的多个散射光子对应的单事件之间空间位置接近(比如空间直线跨度为6个晶体内)、时间间隔小(比如为2-3倍***时间分辨率)且能量小(比如低于435kev)的特点,识别同一晶体内康普顿散射事件的对应的第一扫描数据,为通过晶体内康普顿散射恢复符合事件提供基础。并且,本实施例利用上述条件,还能够消除将模体内康普顿散射事件或其它正常事件误认为晶体内康普顿散射事件的情形,提高准确性。
在通过步骤S101识别出每个晶体内康普顿散射事件对应的第一扫描数据后,可以对第一扫描数据中的每个单事件数据进行标示。该标示用于指示单事件数据为晶体内康普顿散射事件的数据。
在一个示例性的实现过程中,步骤S101之前,还可以包括:
对采集的所述PET扫描数据进行***时间校正,获得校正后的PET扫描数据;
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据,包括:从所述校正后的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据。
本实施例通过对PET扫描数据进行***时间校正,能够补偿因***全局时钟或安装等因素引起的时间误差,提高数据准确性。
在一个示例性的实现过程中,步骤S102可以通过如下方式实现:
对于每组第一扫描数据,将该组第一扫描数据和第二扫描数据中的各个单事件数据按照时间进行排序;
查找与该组第一扫描数据中的所有单事件数据处于同一时间窗的属于第二扫描数据的单事件数据;
将查找出的单事件数据作为与该组第一扫描数据匹配的第一单事件数据。
其中,时间窗的宽度为预设宽度。
举例说明。将前述的单事件a1、单事件a2、单事件b、单事件c和单事件d按照事件中的时间信息所表示的时间进行排序,找到同时包括单事件a1和单事件a2的时间窗1。假设单事件b在时间窗1内,单事件c和单事件d在时间窗1外,则确定单事件b的数据为与单事件a1的数据和单事件a2的数据匹配的第一单事件数据。
在一个示例性的实现过程中,步骤S103可以包括:
根据所述PET扫描数据中符合事件对应的扫描数据进行图像重建,获得初始重建图像;
将所述第一扫描数据中的每一第二单事件数据,分别与所述第一扫描数据匹配的第一单事件数据进行符合处理,获得多条响应线;
计算每条响应线对应的TOF核函数;
对每个TOF核函数与所述初始重建图像中的像素值进行线积分,获得TOF核函数对应的响应线的概率值;
将概率值最大的响应线对应的第二单事件数据和所述第一单事件数据确定为由晶体内康普顿散射事件恢复的符合事件数据。
本实施例中,符合事件对应的扫描数据为第二扫描数据中能够进行符合形成符合事件的单事件数据。
以前述的单事件a1、单事件a2、单事件b、单事件c和单事件d为例,其中,单事件b、单事件c和单事件d的数据都属于第二扫描数据,但是单事件b不能够与第二扫描数据中的其他单事件数据进行符合形成符合事件,而单事件c和单事件d的数据能够进行符合形成符合事件,因此单事件c和单事件d的数据为符合事件对应的扫描数据,而单事件b的数据不为符合事件对应的扫描数据。
初始重建图像是根据所述PET扫描数据中符合事件对应的扫描数据进行图像重建获得的。初始重建图像也即按照前述的相关技术进行图像重建获得的图像。
每组第一扫描数据中的每一第二单事件数据,分别与第一扫描数据匹配的第一单事件数据进行符合处理,即将每一个散射光子的单事件数据与第一单事件数据进行符合处理。
以前述的单事件a1和单事件a2这2个单事件的数据组成的第一扫描数据为例。单事件b的数据为与单事件a1的数据和单事件a2的数据匹配的第一单事件数据,将单事件b的数据与单事件a1数据进行符合处理,获得响应线p1,将单事件b的数据与单事件a2数据进行符合处理,获得响应线p2。
假设在一次湮灭事件中,产生的一个光子按原方向入射到晶体C1,另一个光子经过晶体内康普顿散射后得到的3个散射光子分别入射到晶体C2、晶体C3、晶体C4,则获得的三条响应线如图2所示。图2是本发明实施例提供的晶体内康普顿散射事件形成的多条LOR响应线的示例图。请参见图2,响应线LOR12通过晶体C1和晶体C2,响应线LOR13通过晶体C1和晶体C3,响应线LOR14通过晶体C1和晶体C4。这里的晶体即代表了相应单事件的位置信息所表示的位置。例如,晶体C2代表了一个散射光子对应的单事件的位置信息所表示的位置。图2中,测试模体即为被检体。
本实施例中,TOF核函数可以是根据符合时间差形成的高斯TOF核函数。
本实施例中,响应线的概率值等于响应线对应的TOF核函数与初始重建图像中的像素值的线积分的结果。例如以±3σ为范围,将TOF核函数与初始重建图像中的像素值进行积分,所得的结果即为响应线的概率值。概率值最大的响应线对应的第二单事件数据和第一单事件数据即可作为散射恢复的符合事件数据,该符合事件数据同第二扫描数据中的符合事件对应的扫描数据一起用于图像重建。
本实施例中,通过将响应线对应的TOF核函数与初始重建图像中的像素值进行线积分的方式确定响应线的概率值,进而通过概率值确定散射恢复的符合事件对应的响应线,能够对晶体内康普顿散射事件恢复的符合事件响应线进行准确的空间定位,提高了响应线空间定位的准确性,从而有助于提高重建图像的图像质量。
在一个示例性的实现过程中,步骤S103之前,还可以包括:
所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正;
步骤S103可以包括:根据所述第一扫描数据经过时间游动误差校正后得到的数据及与所述第一单事件数据经过时间游动误差校正后得到的数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
步骤S104可以包括:根据所述第二扫描数据经过时间游动误差校正后得到的数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
例如,对于前述的属于第一扫描数据的单事件a1的数据、单事件a2的数据,以及属于第二扫描数据的单事件b的数据、单事件c的数据、单事件d的数据,在利用单事件a1的数据、单事件a2的数据和单事件b的数据进行散射恢复之前,可以对单事件a1的数据、单事件a2的数据和单事件b的数据、单事件c的数据、单事件d的数据进行时间游动误差校正。
在一个示例中,对所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正,可以包括:
对于所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据中的每组单事件数据,获取该组单事件数据对应的晶体相对于参考晶体的能量与时间游动误差的对应关系;
从所述能量与时间游动误差的对应关系中,获取该组单事件数据中的目标能量值对应的目标时间游动误差;
计算该组单事件数据中的时间值与所述目标时间游动误差的和;
将该组单事件数据中的时间值更新为所述和的值。
本实施例中,晶体指探测器晶体。对于每个晶体,其相对于参考晶体的能量与时间游动误差的对应关系可以通过如下方式获得:
假设有有参考晶体M(也可称为参考探测器晶体),及待测晶体N(也可称为待测探测器晶体),用参考晶体M来测待测晶体N的时间游动误差:
第一步,采集参考晶体M与待测晶体N的所有单事件信息,选取参考晶体M的能量位于511kev±DE的单事件(为了减小参考晶体M的时间游动对待测晶体N测试的影响,DE很小,例如DE为2.5kev),作为参考单事件。
第二步,将待测晶体N采集得到的一系列单事件的能量进行分块,分别与参考单事件进行符合处理。能量块中的单与参考单事件进行符合时,时间差上形成一个符合时间高斯分布图。
第三步,以待测晶体N中能量块(511kev±DE)与参考单事件进行符合形成的高斯时间分布的中点为参考点,将其它能量块形成的高斯符合时间的中点与参考点进行相减得到时间差值,存储该时间差值。
第四步,待测晶体N中一系列的能量块会与参考单事件符合形成一系列的高斯分布时间中点,进而形成一个表,对表中数据进行拟合,得到能量-时间游动误差拟合公式。
该能量-时间游动误差拟合公式即可作为晶体相对于参考晶体的能量与时间游动误差的对应关系。下一次采集待测晶体N的单事件时,依据N的单事件的能量,利用上述的能量-时间游动误差拟合公式,得到时间上的游动误差,并将N的单事件的时间与得到的游动误差相加即为该单事件经过时间游动误差校正后的时间。
本实施例中,时间游动误差校正的方法可以采用相关技术中各种已知的方式进行,关于时间游动误差校正的详细方法请参考相关技术,此处不再进行赘述。
本实施例通过对数据进行时间游动误差校正,对扫描数据进行了校准,提高了数据准确性,有助于提高重建图像的图像质量。
本实施例中,晶体的能量与时间游动误差的对应关系是预先设置在***中的。
本发明实施例提供的图像重建方法,通过从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据,对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,根据第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像,将晶体内康普顿散射事件恢复为符合事件,增加了散射区域的能量值,提高了***灵敏度。
基于上述的方法实施例,本发明实施例还提供了相应的装置、设备、***及存储介质实施例。关于本发明实施例的装置、设备、***及存储介质实施例的详细实现方式,请参见前述方法实施例部分的相应说明。
图3是本发明实施例提供的图像重建装置的功能方块图。如图3所示,本实施例中,图像重建装置可以包括:
识别模块310,用于从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
匹配模块320,用于对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
散射恢复模块330,用于根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
重建模块340,用于根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
在一个示例性的实现过程中,识别模块310可以具体用于:
对于所述PET扫描数据中的每个单事件数据,判断所述PET扫描数据中是否存在满足如下全部条件的目标单事件数据:
包含的空间位置信息所对应的第一空间位置与该单事件数据对应的第二空间位置之间的距离小于预设的距离阈值;
包含的时间信息所对应的第一时间与该单事件数据对应的第二时间之间的差值小于预设的时间阈值;
包含的第一能量值与该单事件数据中的第二能量值均小于预设的第二能量阈值,所述第二能量阈值不大于所述第一能量阈值;
若存在,则确定所述目标单事件数据和该单事件数据为同一散射事件对应的第一扫描数据,所述同一散射事件为晶体内康普顿散射事件。
在一个示例性的实现过程中,散射恢复模块330可以具体用于:
根据所述PET扫描数据中符合事件对应的扫描数据进行图像重建,获得初始重建图像;
将所述第一扫描数据中的每一第二单事件数据,分别与所述第一扫描数据匹配的第一单事件数据进行符合处理,获得多条响应线;
计算每条响应线对应的TOF核函数;
对每个TOF核函数与所述初始重建图像中的像素值进行线积分,获得TOF核函数对应的响应线的概率值;
将概率值最大的响应线对应的第二单事件数据和所述第一单事件数据确定为由晶体内康普顿散射事件恢复的符合事件数据。
在一个示例性的实现过程中,图像重建装置还可以包括:
游动误差校正模块,用于对所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正;
散射恢复模块330具体用于:根据所述第一扫描数据经过时间游动误差校正后得到的数据及与所述第一单事件数据经过时间游动误差校正后得到的数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
重建模块340具体用于:根据所述第二扫描数据经过时间游动误差校正后得到的数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
在一个示例性的实现过程中,游动误差校正模块可以具体用于:
对于所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据中的每组单事件数据,获取该组单事件数据对应的晶体相对于参考晶体的能量与时间游动误差的对应关系;
从所述能量与时间游动误差的对应关系中,获取该组单事件数据中的目标能量值对应的目标时间游动误差;
计算该组单事件数据中的时间值与所述目标时间游动误差的和;
将该组单事件数据中的时间值更新为所述和的值。
在一个示例性的实现过程中,图像重建装置还可以包括:
***校正模块,用于对采集的所述PET扫描数据进行***时间校正,获得校正后的PET扫描数据;
识别模块310可以具体用于:从所述校正后的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据。
本发明实施例还提供了一种控制台设备。图4是本发明实施例提供的控制台设备的一个硬件结构图。如图4所示,控制台设备包括:内部总线401,以及通过内部总线连接的存储器402,处理器403和外部接口404。
所述外部接口404,用于连接PET***的探测器,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;
所述存储器402,用于存储校正逻辑对应的机器可读指令;
所述处理器403,用于读取存储器402上的机器可读指令,并执行所述指令以实现如下操作:
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
在一个示例性的实现过程中,从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据,包括:
对于所述PET扫描数据中的每个单事件数据,判断所述PET扫描数据中是否存在满足如下全部条件的目标单事件数据:
包含的空间位置信息所对应的第一空间位置与该单事件数据对应的第二空间位置之间的距离小于预设的距离阈值;
包含的时间信息所对应的第一时间与该单事件数据对应的第二时间之间的差值小于预设的时间阈值;
包含的第一能量值与该单事件数据中的第二能量值均小于预设的第二能量阈值,所述第二能量阈值不大于所述第一能量阈值;
若存在,则确定所述目标单事件数据和该单事件数据为同一散射事件对应的第一扫描数据,所述同一散射事件为晶体内康普顿散射事件。
在一个示例性的实现过程中,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,包括:
根据所述PET扫描数据中符合事件对应的扫描数据进行图像重建,获得初始重建图像;
将所述第一扫描数据中的每一第二单事件数据,分别与所述第一扫描数据匹配的第一单事件数据进行符合处理,获得多条响应线;
计算每条响应线对应的TOF核函数;
对每个TOF核函数与所述初始重建图像中的像素值进行线积分,获得TOF核函数对应的响应线的概率值;
将概率值最大的响应线对应的第二单事件数据和所述第一单事件数据确定为由晶体内康普顿散射事件恢复的符合事件数据。
在一个示例性的实现过程中,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据之前,还包括:
对所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,包括:根据所述第一扫描数据经过时间游动误差校正后得到的数据及与所述第一单事件数据经过时间游动误差校正后得到的数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,包括:根据所述第二扫描数据经过时间游动误差校正后得到的数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
在一个示例性的实现过程中,对所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正,包括:
对于所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据中的每组单事件数据,获取该组单事件数据对应的晶体相对于参考晶体的能量与时间游动误差的对应关系;
从所述能量与时间游动误差的对应关系中,获取该组单事件数据中的目标能量值对应的目标时间游动误差;
计算该组单事件数据中的时间值与所述目标时间游动误差的和;
将该组单事件数据中的时间值更新为所述和的值。
在一个示例性的实现过程中,从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据之前,还包括:
对采集的所述PET扫描数据进行***时间校正,获得校正后的PET扫描数据;
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据,包括:从所述校正后的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据。
本发明实施例还提供了一种PET***,包括探测器、扫描床和控制台设备,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;其中,
所述闪烁晶体,用于在所述PET***扫描过程中,探测被检体内发出的高能光子,并将所述高能光子转换为光信号;
所述光电转换器件,用于将所述光信号转换成电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述控制台设备,用于从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现如下操作:从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
在一个示例性的实现过程中,从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据,包括:
对于所述PET扫描数据中的每个单事件数据,判断所述PET扫描数据中是否存在满足如下全部条件的目标单事件数据:
包含的空间位置信息所对应的第一空间位置与该单事件数据对应的第二空间位置之间的距离小于预设的距离阈值;
包含的时间信息所对应的第一时间与该单事件数据对应的第二时间之间的差值小于预设的时间阈值;
包含的第一能量值与该单事件数据中的第二能量值均小于预设的第二能量阈值,所述第二能量阈值不大于所述第一能量阈值;
若存在,则确定所述目标单事件数据和该单事件数据为同一散射事件对应的第一扫描数据,所述同一散射事件为晶体内康普顿散射事件。
在一个示例性的实现过程中,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,包括:
根据所述PET扫描数据中符合事件对应的扫描数据进行图像重建,获得初始重建图像;
将所述第一扫描数据中的每一第二单事件数据,分别与所述第一扫描数据匹配的第一单事件数据进行符合处理,获得多条响应线;
计算每条响应线对应的TOF核函数;
对每个TOF核函数与所述初始重建图像中的像素值进行线积分,获得TOF核函数对应的响应线的概率值;
将概率值最大的响应线对应的第二单事件数据和所述第一单事件数据确定为由晶体内康普顿散射事件恢复的符合事件数据。
在一个示例性的实现过程中,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据之前,还包括:
对所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,包括:根据所述第一扫描数据经过时间游动误差校正后得到的数据及与所述第一单事件数据经过时间游动误差校正后得到的数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,包括:根据所述第二扫描数据经过时间游动误差校正后得到的数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
在一个示例性的实现过程中,对所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正,包括:
对于所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据中的每组单事件数据,获取该组单事件数据对应的晶体相对于参考晶体的能量与时间游动误差的对应关系;
从所述能量与时间游动误差的对应关系中,获取该组单事件数据中的目标能量值对应的目标时间游动误差;
计算该组单事件数据中的时间值与所述目标时间游动误差的和;
将该组单事件数据中的时间值更新为所述和的值。
在一个示例性的实现过程中,从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据之前,还包括:
对采集的所述PET扫描数据进行***时间校正,获得校正后的PET扫描数据;
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据,包括:从所述校正后的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据。
对于装置和设备实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本说明书方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本领域技术人员在考虑说明书及实践这里申请的发明后,将容易想到本说明书的其它实施方案。本说明书旨在涵盖本说明书的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本说明书的一般性原理并包括本说明书未申请的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本说明书的真正范围和精神由下面的权利要求指出。
应当理解的是,本说明书并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本说明书的范围仅由所附的权利要求来限制。
以上所述仅为本说明书的较佳实施例而已,并不用以限制本说明书,凡在本说明书的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本说明书保护的范围之内。
Claims (10)
1.一种图像重建方法,其特征在于,所述方法包括:
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
2.根据权利要求1所述的方法,其特征在于,从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据,包括:
对于所述PET扫描数据中的每个单事件数据,判断所述PET扫描数据中是否存在满足如下全部条件的目标单事件数据:
包含的空间位置信息所对应的第一空间位置与该单事件数据对应的第二空间位置之间的距离小于预设的距离阈值;
包含的时间信息所对应的第一时间与该单事件数据对应的第二时间之间的差值小于预设的时间阈值;
包含的第一能量值与该单事件数据中的第二能量值均小于预设的第二能量阈值,所述第二能量阈值不大于所述第一能量阈值;
若存在,则确定所述目标单事件数据和该单事件数据为同一散射事件对应的第一扫描数据,所述同一散射事件为晶体内康普顿散射事件。
3.根据权利要求1所述的方法,其特征在于,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,包括:
根据所述PET扫描数据中符合事件对应的扫描数据进行图像重建,获得初始重建图像;
将所述第一扫描数据中的每一第二单事件数据,分别与所述第一扫描数据匹配的第一单事件数据进行符合处理,获得多条响应线;
计算每条响应线对应的TOF核函数;
对每个TOF核函数与所述初始重建图像中的像素值进行线积分,获得TOF核函数对应的响应线的概率值;
将概率值最大的响应线对应的第二单事件数据和所述第一单事件数据确定为由晶体内康普顿散射事件恢复的符合事件数据。
4.根据权利要求1所述的方法,其特征在于,根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据之前,还包括:
对所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据,包括:根据所述第一扫描数据经过时间游动误差校正后得到的数据及与所述第一单事件数据经过时间游动误差校正后得到的数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,包括:根据所述第二扫描数据经过时间游动误差校正后得到的数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
5.根据权利要求4所述的方法,其特征在于,对所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据进行时间游动误差校正,包括:
对于所述PET扫描数据中的全部第一扫描数据和全部第二扫描数据中的每组单事件数据,获取该组单事件数据对应的晶体相对于参考晶体的能量与时间游动误差的对应关系;
从所述能量与时间游动误差的对应关系中,获取该组单事件数据中的目标能量值对应的目标时间游动误差;
计算该组单事件数据中的时间值与所述目标时间游动误差的和;
将该组单事件数据中的时间值更新为所述和的值。
6.根据权利要求1所述的方法,其特征在于,从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据之前,还包括:
对采集的所述PET扫描数据进行***时间校正,获得校正后的PET扫描数据;
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据,包括:从所述校正后的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据。
7.一种图像重建装置,其特征在于,所述装置包括:
识别模块,用于从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
匹配模块,用于对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
散射恢复模块,用于根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
重建模块,用于根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
8.一种控制台设备,其特征在于,包括内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于连接PET***的探测器,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;
所述存储器,用于存储校正逻辑对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行如下操作:
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
9.一种PET***,其特征在于,包括探测器、扫描床和控制台设备,所述探测器包括多个闪烁晶体及相应的多个光电转换器和处理电路;其中,
所述闪烁晶体,用于在所述PET***扫描过程中,探测被检体内发出的高能光子,并将所述高能光子转换为光信号;
所述光电转换器件,用于将所述光信号转换成电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述控制台设备,用于从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
10.一种计算机可读存储介质,其特征在于,其上存储有计算机程序,其中,所述程序被处理器执行时实现如下操作:
从采集的PET扫描数据中,识别出每个晶体内康普顿散射事件对应的第一扫描数据、以及能量值超过预设的第一能量阈值的第二扫描数据;
对于每组第一扫描数据,从第二扫描数据中找出与该第一扫描数据匹配的第一单事件数据;
根据所述第一扫描数据及与所述第一扫描数据匹配的第一单事件数据,获得由晶体内康普顿散射事件恢复的符合事件数据;
根据所述第二扫描数据中符合事件对应的数据以及所述由晶体内康普顿散射事件恢复的符合事件数据进行图像重建,得到重建图像。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910945559.4A CN110680370B (zh) | 2019-09-30 | 2019-09-30 | 图像重建方法、装置、控制台设备及pet*** |
US17/037,301 US11925498B2 (en) | 2019-09-30 | 2020-09-29 | Reconstructing image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910945559.4A CN110680370B (zh) | 2019-09-30 | 2019-09-30 | 图像重建方法、装置、控制台设备及pet*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110680370A true CN110680370A (zh) | 2020-01-14 |
CN110680370B CN110680370B (zh) | 2023-02-28 |
Family
ID=69111408
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910945559.4A Active CN110680370B (zh) | 2019-09-30 | 2019-09-30 | 图像重建方法、装置、控制台设备及pet*** |
Country Status (2)
Country | Link |
---|---|
US (1) | US11925498B2 (zh) |
CN (1) | CN110680370B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112022189A (zh) * | 2020-09-09 | 2020-12-04 | 东软医疗***股份有限公司 | 图像重建方法、装置、设备及pet*** |
WO2023035361A1 (zh) * | 2021-09-10 | 2023-03-16 | 苏州瑞派宁科技有限公司 | 图像重建方法、装置、***及计算机可读存储介质 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220343566A1 (en) * | 2021-04-22 | 2022-10-27 | Canon Medical Systems Corporation | Methods and systems for reconstructing a positron emission tomography image |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102246057A (zh) * | 2008-12-10 | 2011-11-16 | 皇家飞利浦电子股份有限公司 | 用作可缩放pet和spect***构建块的自主探测器模块 |
CN103164863A (zh) * | 2013-04-02 | 2013-06-19 | 中国科学院高能物理研究所 | 用于重建正电子发射计算机断层成像图像的方法 |
US20160131774A1 (en) * | 2013-06-24 | 2016-05-12 | Eduardo M. Lage | Normalization correction for multiple-detection enhanced emission tomography |
CN106388841A (zh) * | 2016-09-23 | 2017-02-15 | 沈阳东软医疗***有限公司 | 一种图像重建方法、装置及设备 |
CN107978002A (zh) * | 2016-10-25 | 2018-05-01 | 上海东软医疗科技有限公司 | 一种pet图像重建方法和装置 |
CN109009199A (zh) * | 2017-08-26 | 2018-12-18 | 上海联影医疗科技有限公司 | 用于正电子发射断层扫描中的图像数据处理的***和方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7381959B2 (en) * | 2005-08-17 | 2008-06-03 | General Electric Company | Technique for reconstructing PET scan images |
CN101809463B (zh) * | 2007-09-24 | 2013-05-22 | 皇家飞利浦电子股份有限公司 | 临床前飞行时间pet成像 |
US8269177B2 (en) * | 2010-02-28 | 2012-09-18 | General Electric Company | Multiplexing readout scheme for a gamma ray detector |
US11029427B2 (en) * | 2018-11-12 | 2021-06-08 | The Board Of Trustees Of The University Of Illinois | Method and system for increasing radiation sensitivity in semiconductor detectors |
-
2019
- 2019-09-30 CN CN201910945559.4A patent/CN110680370B/zh active Active
-
2020
- 2020-09-29 US US17/037,301 patent/US11925498B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102246057A (zh) * | 2008-12-10 | 2011-11-16 | 皇家飞利浦电子股份有限公司 | 用作可缩放pet和spect***构建块的自主探测器模块 |
CN103164863A (zh) * | 2013-04-02 | 2013-06-19 | 中国科学院高能物理研究所 | 用于重建正电子发射计算机断层成像图像的方法 |
US20160131774A1 (en) * | 2013-06-24 | 2016-05-12 | Eduardo M. Lage | Normalization correction for multiple-detection enhanced emission tomography |
CN106388841A (zh) * | 2016-09-23 | 2017-02-15 | 沈阳东软医疗***有限公司 | 一种图像重建方法、装置及设备 |
CN107978002A (zh) * | 2016-10-25 | 2018-05-01 | 上海东软医疗科技有限公司 | 一种pet图像重建方法和装置 |
CN109009199A (zh) * | 2017-08-26 | 2018-12-18 | 上海联影医疗科技有限公司 | 用于正电子发射断层扫描中的图像数据处理的***和方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112022189A (zh) * | 2020-09-09 | 2020-12-04 | 东软医疗***股份有限公司 | 图像重建方法、装置、设备及pet*** |
CN112022189B (zh) * | 2020-09-09 | 2023-04-07 | 东软医疗***股份有限公司 | 图像重建方法、装置、设备及pet*** |
WO2023035361A1 (zh) * | 2021-09-10 | 2023-03-16 | 苏州瑞派宁科技有限公司 | 图像重建方法、装置、***及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
US20210093279A1 (en) | 2021-04-01 |
US11925498B2 (en) | 2024-03-12 |
CN110680370B (zh) | 2023-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gong et al. | Machine learning in PET: from photon detection to quantitative image reconstruction | |
CN110680370B (zh) | 图像重建方法、装置、控制台设备及pet*** | |
US10502846B2 (en) | Normalization correction for multiple-detection enhanced emission tomography | |
US10319118B2 (en) | Reconstructing PET image | |
US8098916B2 (en) | System and method for image-based attenuation correction of PET/SPECT images | |
US8600136B2 (en) | Method for generation of attenuation map in PET-MR | |
CN110392847B (zh) | 用于光子计数边缘上x射线探测器的增加的空间分辨率 | |
US10215864B2 (en) | System and method to improve image quality of emission tomography when using advanced radionuclides | |
US20080287772A1 (en) | Motion Compensation in PET Reconstruction | |
CN103890571B (zh) | 具有偏移校正的用于探测光子的放射探测装置 | |
CN109887048B (zh) | Pet散射校正方法、图像重建方法、装置及电子设备 | |
CN105378508B (zh) | 直接转换辐射探测器数字信号处理电子设备 | |
CN109242920B (zh) | 物质分解方法、装置和*** | |
CN110327067B (zh) | 图像重建方法、装置、终端设备及pet*** | |
JPH06347555A (ja) | ポジトロンイメージング装置 | |
RU2493572C2 (ru) | Геометрические преобразования, сохраняющие данные в виде списка | |
CN110988974B (zh) | 曲面辐射源的pet探测器飞行时间的校准*** | |
EP2917767A1 (en) | Inter-detector scatter enhanced emission tomography | |
WO2011119478A2 (en) | Transient pixel defect detection and correction | |
US6294788B1 (en) | Randoms correction in positron imaging | |
US10466368B2 (en) | Systems and methods for detector output adjustment | |
US20090050818A1 (en) | Method for determining individual quantum absorption events in a radiation converter for converting individual quanta of ionizing radiation incident on it, program code segment for carrying out the method, device for electronic data processing, radiation converter and imaging tomography scanner | |
US20230022425A1 (en) | Apparatus, system, method and computer probram for providing a nuclear image of a region of interest of a patient | |
CN112285761A (zh) | 一种面向植物pet的晶***置谱数据分割***与方法 | |
JP4317124B2 (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230410 Address after: Room 308, No. 177-2 Chuangxin Road, Hunnan District, Shenyang City, Liaoning Province, 110167 Patentee after: Shenyang Zhihe Medical Technology Co.,Ltd. Address before: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee before: Shenyang Neusoft Medical Systems Co.,Ltd. |