CN110211095A - 能量图及晶体位置查找表生成方法、装置、存储介质 - Google Patents
能量图及晶体位置查找表生成方法、装置、存储介质 Download PDFInfo
- Publication number
- CN110211095A CN110211095A CN201910371305.6A CN201910371305A CN110211095A CN 110211095 A CN110211095 A CN 110211095A CN 201910371305 A CN201910371305 A CN 201910371305A CN 110211095 A CN110211095 A CN 110211095A
- Authority
- CN
- China
- Prior art keywords
- energy value
- spectral distribution
- scintillation crystal
- energy
- distribution curve
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- 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
- G16H30/00—ICT specially adapted for the handling or processing of medical images
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Epidemiology (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Quality & Reliability (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本申请提供一种能量图及晶体位置查找表生成方法、装置、存储介质及医疗设备,用以提高能量图的图像质量,该能量图生成方法包括:获取PET探测器的闪烁晶体的原始本底辐射数据;根据所述原始本底辐射数据确定各个所述闪烁晶体的能量值阈值区间;获得对应的能量值在所述能量值阈值区间内的目标本底辐射数据;基于所述目标本底辐射数据生成能量图。
Description
技术领域
本申请涉及医学成像技术领域,尤其涉及一种能量图及晶体位置查找表生成方法、装置、存储介质及医疗设备。
背景技术
正电子发射断层成像(Positron Emission Tomography,PET)是一种非侵入式的造影方法,目前主流的PET探测器通常采用闪烁晶体阵列耦合光电转换器件的设计方式。PET系统一般需要先建立晶体位置查找表以记录γ光子事件位置坐标与闪烁晶体之间的对应关系,然后在临床应用中,通过晶体位置查找表可以确定与γ光子发生作用的闪烁晶体,从而获得该闪烁晶体的实际物理位置,用于后续图像重建,因此晶体位置查找表的准确性直接影响了PET系统的空间分辨率。
PET探测器的闪烁晶体中通常含有镥元素(176Lu),相关技术中提出直接采用采集到的PET探测器的闪烁晶体的本底辐射数据生成能量图,然后根据能量图计算生成晶体位置查找表。然而,由于镥元素在衰变过程中产生多个能量的γ光子和β粒子,如果直接采用采集到的本底辐射数据生成能量图将导致获得的能量图噪声较高,图像对比度较差,从该能量图中无法找到所有闪烁晶体对应的光斑峰值点,以致生成的晶体位置查找表不够准确。
发明内容
有鉴于此,本申请提供了一种能量图及晶体位置查找表生成方法、装置、存储介质及医疗设备,用以提高能量图的图像质量。
第一方面,本申请实施例提供了一种能量图生成方法,所述方法包括:
获取PET探测器的闪烁晶体的原始本底辐射数据;
根据所述原始本底辐射数据确定各个所述闪烁晶体的能量值阈值区间;
获得对应的能量值在所述能量值阈值区间内的目标本底辐射数据;
基于所述目标本底辐射数据生成能量图。
上述方法,先根据采集的PET探测器的闪烁晶体的原始本底辐射数据确定各个闪烁晶体的能量值阈值区间,然后基于该能量值阈值区间对原始本底辐射数据进行筛选,获得对应的能量值在能量值阈值区间内的目标本底辐射数据,之后基于目标本底辐射数据生成能量图,通过对原始本底辐射数据进行筛选,可以使得生成的能量图噪声减少,图像对比度得到提高,从而使得能量图的图像质量得到明显提升。
在一可能的实现方式中,所述根据所述原始本底辐射数据确定各个所述闪烁晶体的能量值阈值区间,包括:
根据所述原始本底辐射数据生成各个所述闪烁晶体对应的能谱分布曲线;
针对每一个所述闪烁晶体,根据该闪烁晶体对应的能谱分布曲线确定能量值阈值区间。
在一可能的实现方式中,所述根据闪烁晶体对应的能谱分布曲线确定能量值阈值区间,包括:
确定闪烁晶体对应的能谱分布曲线的峰值点;
选取所述能谱分布曲线上的一个峰值点,并根据选取到的峰值点确定能量值阈值区间,所述能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
该方法中,通过先确定闪烁晶体对应的能谱分布曲线的峰值点,再选取能谱分布曲线上的一个峰值点,并根据选取到的峰值点确定能量值阈值区间,由于确定的能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点,可以理解成只保留单一能量(γ光子或β粒子)事件,这样生成的能量图的图像噪声较少,图像对比度较高,因此该能量图中每个闪烁晶体对应的光斑更清晰。
在一可能的实现方式中,所述根据闪烁晶体对应的能谱分布曲线确定能量值阈值区间,包括:
利用波峰或波谷将闪烁晶体对应的能谱分布曲线划分成多个区域;
选取其中一个区域作为目标区域,确定该目标区域的能谱分布曲线的峰值点;
根据确定的所述峰值点确定能量值阈值区间,所述能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
该方法中,先利用波峰或波谷将闪烁晶体对应的能谱分布曲线划分成多个区域,再选取其中一个区域作为目标区域,确定该目标区域的能谱分布曲线的峰值点,然后根据确定的所述峰值点确定能量值阈值区间,由于确定的能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点,可以理解成只保留单一能量(γ光子或β粒子)事件,这样生成的能量图的图像噪声较少,图像对比度较高,因此该能量图中每个闪烁晶体对应的光斑更清晰。
在一可能的实现方式中,所述利用波峰将闪烁晶体对应的能谱分布曲线划分成多个区域,包括:
确定闪烁晶体对应的能谱分布曲线的各波峰对应的能量值;
计算各个相邻的波峰对应的能量值的中间值;
基于各个相邻的波峰对应的能量值的中间值将所述闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
在一可能的实现方式中,所述利用波谷将闪烁晶体对应的能谱分布曲线划分成多个区域,包括:
确定闪烁晶体对应的能谱分布曲线的各波谷对应的能量值;
基于各波谷对应的能量值将所述闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
第二方面,本申请实施例还提供了一种晶体位置查找表生成方法,所述方法包括:
采用第一方面或第一方面的任意可能的实现方式中的能量图生成方法生成能量图;
根据所述能量图生成晶体位置查找表。
上述方法,通过对原始本底辐射数据进行筛选,可以使得生成的能量图的图像质量得到明显提升,该能量图中每个闪烁晶体对应的光斑更清晰,从而可以准确找到闪烁晶体对应的光斑峰值点,进而可以提高生成的晶体位置查找表的准确度。
第三方面,本申请实施例还提供了一种能量图生成装置,包括用于执行第一方面或第一方面的任意可能的实现方式中的能量图生成方法的模块。
第四方面,本申请实施例还提供了一种晶体位置查找表生成装置,包括用于执行第二方面或第二方面的任意可能的实现方式中的晶体位置查找表生成方法的模块。
第五方面,本申请实施例还提供了一种存储介质,其上存储有计算机程序,所述程序被处理器执行时实现第二方面或第二方面的任意可能的实现方式中的晶体位置查找表生成方法的步骤。
第六方面,本申请实施例还提供了一种医疗设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现第二方面或第二方面的任意可能的实现方式中的晶体位置查找表生成方法的步骤。
附图说明
图1为本申请实施例中176Lu同位素的衰变能级图;
图2为本申请实施例中闪烁晶体的本底辐射能谱分布曲线示意图;
图3为本申请实施例中PET探测器中闪烁晶体的分布示意图;
图4为采用图3所示的闪烁晶体的本底辐射数据直接生成的能量图;
图5为本申请实施例提供的一种能量图生成方法的流程示意图;
图6为对图3所示的闪烁晶体的本底辐射数据进行筛选后获得的目标本底辐射数据生成的能量图;
图7为本申请实施例中能量值阈值区间的示意图;
图8为本申请实施例中能谱分布曲线的区域划分示意图;
图9为本申请实施例提供的一种晶体位置查找表生成方法的流程示意图;
图10为本申请实施例提供的一种能量图生成装置的结构示意图;
图11为本申请实施例提供的能量图生成装置中阈值确定模块的结构示意图;
图12为本申请实施例提供的能量图生成装置中阈值确定子模块的第一种结构示意图;
图13为本申请实施例提供的能量图生成装置中阈值确定子模块的第二种结构示意图;
图14为本申请实施例提供的能量图生成装置中区域划分模块的第一种结构示意图;
图15为本申请实施例提供的能量图生成装置中区域划分模块的第二种结构示意图;
图16为本申请实施例提供的一种晶体位置查找表生成装置的结构示意图;
图17为本申请实施例提供的医疗设备的结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
在本申请使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本申请可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本申请范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
PET探测器的闪烁晶体中通常含有镥元素(176Lu),176Lu同位素的半衰周期约为3.8*1010年,故在十几年的使用时间内,由于本底辐射产生的计数基本可以看作是不变的。176Lu同位素的衰变能级图如图1所示。可以看到有99.6%的概率发生β负衰变,且衰变后有94%概率会发生能级跃迁释放出一个307keV的γ光子,然后再先后分别有78%和15%的概率发生能级跃迁产生202keV和88keV的γ光子。闪烁晶体的本底辐射能谱分布曲线示意图如图2所示,可以看到本底辐射能谱产生多个峰值,其中88keV的γ光子由于能量较低且占比较少其峰值一般不可见,β粒子和γ光子组合而成的能谱占据主导特性,307keV和202keV的γ光子的峰值可见但一般不是最显著的。从能谱分布曲线可以看到,由于本底辐射能谱分布是由多个能谱组合形成的(即由88keV的γ光子的能谱、202keV的γ光子的能谱、307keV的γ光子的能谱、以及β粒子和γ光子组合而成的能谱组合形成的)。由于镥元素在衰变过程中产生多个能量的γ光子和β粒子,如果直接采用采集到的PET探测器的闪烁晶体的本底辐射数据生成能量图将导致获得的能量图噪声较高,图像对比度较差,从该能量图中无法找到所有闪烁晶体对应的光斑峰值点,以致生成的晶体位置查找表不够准确。
例如,PET探测器的一个模块包含11x11个晶体,如图3所示,图3中数字表示晶体的编号,直接采用采集到的本底辐射数据生成能量图(flood histogram),所得结果如图4所示。从图4中可以看到,对于该PET探测器来说,直接采用采集到的本底辐射数据生成的能量图噪声较高,图像对比度较差,很难找到所有121个晶体对应的光斑峰值点,这样就无法准确计算出晶体位置查找表,最终影响PET系统的成像分辨率。
针对上述问题,本申请提供了一种能量图及晶体位置查找表生成方法、装置。
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。
参见图5,本申请实施例提供了一种能量图生成方法,该方法可以应用于PET系统,该方法可以包括如下步骤:
S101、获取PET探测器的闪烁晶体的原始本底辐射数据;
PET探测器中常用闪烁晶体包括硅酸镥(LSO)及硅酸钇镥(LYSO)等含镥晶体。
S102、根据原始本底辐射数据确定各个闪烁晶体的能量值阈值区间;
S103、获得对应的能量值在能量值阈值区间内的目标本底辐射数据;
在一些实施例中,上述获得对应的能量值在能量值阈值区间内的目标本底辐射数据可以包括:
对原始本底辐射数据进行筛选,保留对应的能量值在能量值阈值区间内的本底辐射数据,作为目标本底辐射数据。
S104、基于目标本底辐射数据生成能量图。
其中,能量图是指对光子事件位置进行分布统计所形成的图谱,步骤S104可以采用常用的方法执行。
本申请实施例提供的方法,通过对原始本底辐射数据根据各个闪烁晶体的能量值阈值区间进行筛选以获得目标本底辐射数据,然后基于目标本底辐射数据生成能量图,在一可能的实现方式中生成的能量图如图6所示,相比于图4,图6中能量图的图像噪声明显减少,图像对比度得到提升,每个闪烁晶体对应的光斑更清晰。
在一可能的实现方式中,步骤S102中根据原始本底辐射数据确定各个闪烁晶体的能量值阈值区间,可以包括:
根据原始本底辐射数据生成各个闪烁晶体对应的能谱分布曲线;
针对每一个闪烁晶体,根据该闪烁晶体对应的能谱分布曲线确定能量值阈值区间。
上述根据闪烁晶体对应的能谱分布曲线确定能量值阈值区间可以有多种实现方式,下面分别进行介绍。
实现方式一:
上述根据闪烁晶体对应的能谱分布曲线确定能量值阈值区间,可以包括:
确定闪烁晶体对应的能谱分布曲线的峰值点;
选取该能谱分布曲线上的一个峰值点,并根据选取到的峰值点确定能量值阈值区间,该能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
其中,峰值点的确定可以采用微商寻峰法、B样条插值法等。峰值点的选取可以任意选取,也可以选取对应的能量值与预设能量值(例如307KeV)最接近的峰值点。
在一些实施例中,根据选取到的峰值点确定能量值阈值区间,可以包括:
确定与选取到的峰值点相邻的波谷对应的能量值;
将波谷对应的能量值作为能量值阈值区间的端点能量值阈值。
通过上述方法确定的能量值阈值区间可以如图7中区间1所示。
在另一些实施例中,根据选取到的峰值点确定能量值阈值区间,可以包括:
将与选取到的峰值点对应的能量值相隔预设值(例如150keV)的能量值作为能量值阈值区间的端点阈值。
通过上述方法确定的能量值阈值区间可以如图7中区间2所示。
当然,也可以将上述两种方法结合起来以确定能量值阈值区间,本申请实施例对此不进行限定。
实现方式二:
上述根据闪烁晶体对应的能谱分布曲线确定能量值阈值区间,可以包括:
利用波峰将闪烁晶体对应的能谱分布曲线划分成多个区域;
选取其中一个区域作为目标区域,确定该目标区域的能谱分布曲线的峰值点;
根据确定的峰值点确定能量值阈值区间,该能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
其中,目标区域的选取可以是选取任意一个区域,也可以是选取的区域的能谱分布曲线对应的能量值区间包括预设能量值(例如307KeV)。
在一可能的实现方式中,上述利用波峰将闪烁晶体对应的能谱分布曲线划分成多个区域,可以包括:
确定闪烁晶体对应的能谱分布曲线的各波峰对应的能量值;
计算各个相邻的波峰对应的能量值的中间值;
基于各个相邻的波峰对应的能量值的中间值将该闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
需要说明的是,实现方式二中根据确定的峰值点确定能量值阈值区间的方法可以采用实现方式一中根据选取到的峰值点确定能量值阈值区间的方法,在此不再赘述。
实现方式三:
上述根据闪烁晶体对应的能谱分布曲线确定能量值阈值区间,可以包括:
利用波谷将闪烁晶体对应的能谱分布曲线划分成多个区域;
选取其中一个区域作为目标区域,确定该目标区域的能谱分布曲线的峰值点;
根据确定的峰值点确定能量值阈值区间,该能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
在一可能的实现方式中,上述利用波谷将闪烁晶体对应的能谱分布曲线划分成多个区域,可以包括:
确定闪烁晶体对应的能谱分布曲线的各波谷对应的能量值;
基于各波谷对应的能量值将闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
例如,如图8所示,可以将能谱分布曲线按能量值划分成A、B和C三个区域。
需要说明的是,实现方式三中根据确定的峰值点确定能量值阈值区间的方法也可以采用实现方式一中根据选取到的峰值点确定能量值阈值区间的方法,在此不再赘述。
基于同一发明构思,参见图9,本申请实施例还提供了一种晶体位置查找表生成方法,该方法可以应用于PET系统,该方法包括如下步骤:
采用本申请任意实施例提供的能量图生成方法生成能量图;
S105、根据所述能量图生成晶体位置查找表。
也就是说,本申请实施例提供的晶体位置查找表生成方法除了包括上述步骤S101-S104,还包括步骤S105。
步骤S105可以采用常用的方法执行,例如一常用的方法为:从生成的能量图中识别每一个闪烁晶体对应的光斑峰值点,并确定该光斑峰值点的位置坐标,将光斑峰值点的位置坐标作为晶体中心位置,根据晶体中心位置划分各晶体的边界,并根据边界完成晶体区域的划分,以划分的晶体区域和晶体编号生成晶体位置查找表。
基于同一发明构思,参见图10,本申请实施例还提供了一种能量图生成装置,该装置包括:数据获取模块11、阈值确定模块12、筛选模块13和能量图生成模块14。
其中,数据获取模块11,用于获取PET探测器的闪烁晶体的原始本底辐射数据;
阈值确定模块12,用于根据原始本底辐射数据确定各个闪烁晶体的能量值阈值区间;
筛选模块13,用于获得对应的能量值在能量值阈值区间内的目标本底辐射数据;
能量图生成模块14,用于基于目标本底辐射数据生成能量图。
在一可能的实现方式中,如图11所示,阈值确定模块12可以包括:
曲线生成模块121,用于根据原始本底辐射数据生成各个闪烁晶体对应的能谱分布曲线;
阈值确定子模块122,用于针对每一个闪烁晶体,根据该闪烁晶体对应的能谱分布曲线确定能量值阈值区间。
在一可能的实现方式中,如图12所示,阈值确定子模块122可以包括:
第一曲线峰值确定模块201,用于确定闪烁晶体对应的能谱分布曲线的峰值点;
第一阈值区间确定模块202,用于选取该能谱分布曲线上的一个峰值点,并根据选取到的峰值点确定能量值阈值区间,该能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
在另一可能的实现方式中,如图13所示,阈值确定子模块122可以包括:
区域划分模块203,用于利用波峰或波谷将闪烁晶体对应的能谱分布曲线划分成多个区域;
第二曲线峰值确定模块204,用于选取其中一个区域作为目标区域,确定该目标区域的能谱分布曲线的峰值点;
第二阈值区间确定模块205,用于根据确定的峰值点确定能量值阈值区间,该能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
在一可能的实现方式中,如图14所示,区域划分模块203可以包括:
曲线波峰确定模块301,用于确定闪烁晶体对应的能谱分布曲线的各波峰对应的能量值;
中间值确定模块302,用于计算各个相邻的波峰对应的能量值的中间值;
第一区域划分子模块303,用于基于各个相邻的波峰对应的能量值的中间值将闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
在另一可能的实现方式中,如图15所示,区域划分模块203可以包括:
曲线波谷确定模块304,用于确定闪烁晶体对应的能谱分布曲线的各波谷对应的能量值;
第二区域划分子模块305,用于基于各波谷对应的能量值将闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
基于同一发明构思,参见图16,本申请实施例还提供了一种晶体位置查找表生成装置,该装置包括:本申请任意实施例提供的能量图生成装置10和查找表生成模块15。
其中,查找表生成模块15,用于根据能量图生成装置10所生成的能量图生成晶体位置查找表。
上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程,在此不再赘述。
对于装置实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本申请方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
基于同一发明构思,本申请实施例还提供了一种存储介质,其上存储有计算机程序,所述程序被处理器执行时实现上述任意可能的实现方式中的晶体位置查找表生成方法的步骤。
可选地,该存储介质具体可以为存储器。
基于同一发明构思,参见图17,本申请实施例还提供了一种医疗设备,包括存储器71(例如非易失性存储器)、处理器72及存储在存储器71上并可在处理器72上运行的计算机程序,处理器72执行所述程序时实现上述任意可能的实现方式中的晶体位置查找表生成方法的步骤。该医疗设备例如可以为PC,用于PET晶体位置查找表生成,属于PET系统,与PET探测器连接。
如图17所示,该医疗设备一般还可以包括:内存73、网络接口74、以及内部总线75。除了这些部件外,还可以包括其他硬件,对此不再赘述。
需要指出的是,上述晶体位置查找表生成装置可以通过软件实现,其作为一个逻辑意义上的装置,是通过其所在的医疗设备的处理器72将非易失性存储器中存储的计算机程序指令读取到内存73中运行形成的。
本说明书中描述的主题及功能操作的实施例可以在以下中实现:数字电子电路、有形体现的计算机软件或固件、包括本说明书中公开的结构及其结构性等同物的计算机硬件、或者它们中的一个或多个的组合。本说明书中描述的主题的实施例可以实现为一个或多个计算机程序,即编码在有形非暂时性程序载体上以被数据处理装置执行或控制数据处理装置的操作的计算机程序指令中的一个或多个模块。可替代地或附加地,程序指令可以被编码在人工生成的传播信号上,例如机器生成的电、光或电磁信号,该信号被生成以将信息编码并传输到合适的接收机装置以由数据处理装置执行。计算机存储介质可以是机器可读存储设备、机器可读存储基板、随机或串行存取存储器设备、或它们中的一个或多个的组合。
本说明书中描述的处理及逻辑流程可以由执行一个或多个计算机程序的一个或多个可编程计算机执行,以通过根据输入数据进行操作并生成输出来执行相应的功能。所述处理及逻辑流程还可以由专用逻辑电路—例如FPGA(现场可编程门阵列)或ASIC(专用集成电路)来执行,并且装置也可以实现为专用逻辑电路。
适合用于执行计算机程序的计算机包括,例如通用和/或专用微处理器,或任何其他类型的中央处理单元。通常,中央处理单元将从只读存储器和/或随机存取存储器接收指令和数据。计算机的基本组件包括用于实施或执行指令的中央处理单元以及用于存储指令和数据的一个或多个存储器设备。通常,计算机还将包括用于存储数据的一个或多个大容量存储设备,例如磁盘、磁光盘或光盘等,或者计算机将可操作地与此大容量存储设备耦接以从其接收数据或向其传送数据,抑或两种情况兼而有之。然而,计算机不是必须具有这样的设备。此外,计算机可以嵌入在另一设备中,例如移动电话、个人数字助理(PDA)、移动音频或视频播放器、游戏操纵台、全球定位系统(GPS)接收机、或例如通用串行总线(USB)闪存驱动器的便携式存储设备,仅举几例。
适合于存储计算机程序指令和数据的计算机可读介质包括所有形式的非易失性存储器、媒介和存储器设备,例如包括半导体存储器设备(例如EPROM、EEPROM和闪存设备)、磁盘(例如内部硬盘或可移动盘)、磁光盘以及CD ROM和DVD-ROM盘。处理器和存储器可由专用逻辑电路补充或并入专用逻辑电路中。
虽然本说明书包含许多具体实施细节,但是这些不应被解释为限制任何发明的范围或所要求保护的范围,而是主要用于描述特定发明的具体实施例的特征。本说明书内在多个实施例中描述的某些特征也可以在单个实施例中被组合实施。另一方面,在单个实施例中描述的各种特征也可以在多个实施例中分开实施或以任何合适的子组合来实施。此外,虽然特征可以如上所述在某些组合中起作用并且甚至最初如此要求保护,但是来自所要求保护的组合中的一个或多个特征在一些情况下可以从该组合中去除,并且所要求保护的组合可以指向子组合或子组合的变型。
类似地,虽然在附图中以特定顺序描绘了操作,但是这不应被理解为要求这些操作以所示的特定顺序执行或顺次执行、或者要求所有例示的操作被执行,以实现期望的结果。在某些情况下,多任务和并行处理可能是有利的。此外,上述实施例中的各种系统模块和组件的分离不应被理解为在所有实施例中均需要这样的分离,并且应当理解,所描述的程序组件和系统通常可以一起集成在单个软件产品中,或者封装成多个软件产品。
由此,主题的特定实施例已被描述。其他实施例在所附权利要求书的范围以内。在某些情况下,权利要求书中记载的动作可以以不同的顺序执行并且仍实现期望的结果。此外,附图中描绘的处理并非必需所示的特定顺序或顺次顺序,以实现期望的结果。在某些实现中,多任务和并行处理可能是有利的。
以上所述仅为本申请的较佳实施例而已,并不用以限制本申请,凡在本申请的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本申请保护的范围之内。
Claims (16)
1.一种能量图生成方法,其特征在于,所述方法包括:
获取PET探测器的闪烁晶体的原始本底辐射数据;
根据所述原始本底辐射数据确定各个所述闪烁晶体的能量值阈值区间;
获得对应的能量值在所述能量值阈值区间内的目标本底辐射数据;
基于所述目标本底辐射数据生成能量图。
2.根据权利要求1所述的方法,其特征在于,所述根据所述原始本底辐射数据确定各个所述闪烁晶体的能量值阈值区间,包括:
根据所述原始本底辐射数据生成各个所述闪烁晶体对应的能谱分布曲线;
针对每一个所述闪烁晶体,根据该闪烁晶体对应的能谱分布曲线确定能量值阈值区间。
3.根据权利要求2所述的方法,其特征在于,所述根据闪烁晶体对应的能谱分布曲线确定能量值阈值区间,包括:
确定闪烁晶体对应的能谱分布曲线的峰值点;
选取所述能谱分布曲线上的一个峰值点,并根据选取到的峰值点确定能量值阈值区间,所述能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
4.根据权利要求2所述的方法,其特征在于,所述根据闪烁晶体对应的能谱分布曲线确定能量值阈值区间,包括:
利用波峰或波谷将闪烁晶体对应的能谱分布曲线划分成多个区域;
选取其中一个区域作为目标区域,确定该目标区域的能谱分布曲线的峰值点;
根据确定的所述峰值点确定能量值阈值区间,所述能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
5.根据权利要求4所述的方法,其特征在于,所述利用波峰将闪烁晶体对应的能谱分布曲线划分成多个区域,包括:
确定闪烁晶体对应的能谱分布曲线的各波峰对应的能量值;
计算各个相邻的波峰对应的能量值的中间值;
基于各个相邻的波峰对应的能量值的中间值将所述闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
6.根据权利要求4所述的方法,其特征在于,所述利用波谷将闪烁晶体对应的能谱分布曲线划分成多个区域,包括:
确定闪烁晶体对应的能谱分布曲线的各波谷对应的能量值;
基于各波谷对应的能量值将所述闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
7.一种晶体位置查找表生成方法,其特征在于,所述方法包括:
采用权利要求1-6任一项所述的能量图生成方法生成能量图;
根据所述能量图生成晶体位置查找表。
8.一种能量图生成装置,其特征在于,所述装置包括:
数据获取模块,用于获取PET探测器的闪烁晶体的原始本底辐射数据;
阈值确定模块,用于根据所述原始本底辐射数据确定各个所述闪烁晶体的能量值阈值区间;
筛选模块,用于获得对应的能量值在所述能量值阈值区间内的目标本底辐射数据;
能量图生成模块,用于基于所述目标本底辐射数据生成能量图。
9.根据权利要求8所述的装置,其特征在于,所述阈值确定模块包括:
曲线生成模块,用于根据所述原始本底辐射数据生成各个所述闪烁晶体对应的能谱分布曲线;
阈值确定子模块,用于针对每一个所述闪烁晶体,根据该闪烁晶体对应的能谱分布曲线确定能量值阈值区间。
10.根据权利要求9所述的装置,其特征在于,所述阈值确定子模块包括:
第一曲线峰值确定模块,用于确定闪烁晶体对应的能谱分布曲线的峰值点;
第一阈值区间确定模块,用于选取所述能谱分布曲线上的一个峰值点,并根据选取到的峰值点确定能量值阈值区间,所述能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
11.根据权利要求9所述的装置,其特征在于,所述阈值确定子模块包括:
区域划分模块,用于利用波峰或波谷将闪烁晶体对应的能谱分布曲线划分成多个区域;
第二曲线峰值确定模块,用于选取其中一个区域作为目标区域,确定该目标区域的能谱分布曲线的峰值点;
第二阈值区间确定模块,用于根据确定的所述峰值点确定能量值阈值区间,所述能量值阈值区间对应的能谱分布曲线上仅包含一个峰值点。
12.根据权利要求11所述的装置,其特征在于,所述区域划分模块包括:
曲线波峰确定模块,用于确定闪烁晶体对应的能谱分布曲线的各波峰对应的能量值;
中间值确定模块,用于计算各个相邻的波峰对应的能量值的中间值;
第一区域划分子模块,用于基于各个相邻的波峰对应的能量值的中间值将所述闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
13.根据权利要求11所述的装置,其特征在于,所述区域划分模块包括:
曲线波谷确定模块,用于确定闪烁晶体对应的能谱分布曲线的各波谷对应的能量值;
第二区域划分子模块,用于基于各波谷对应的能量值将所述闪烁晶体对应的能谱分布曲线按能量值划分成多个区域。
14.一种晶体位置查找表生成装置,其特征在于,所述装置包括:权利要求8-13任一项所述的能量图生成装置和查找表生成模块;
所述查找表生成模块,用于根据所述能量图生成装置所生成的能量图生成晶体位置查找表。
15.一种存储介质,其上存储有计算机程序,其特征在于,所述程序被处理器执行时实现权利要求7所述的晶体位置查找表生成方法的步骤。
16.一种医疗设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求7所述的晶体位置查找表生成方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910371305.6A CN110211095B (zh) | 2019-05-06 | 2019-05-06 | 能量图及晶体位置查找表生成方法、装置、存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910371305.6A CN110211095B (zh) | 2019-05-06 | 2019-05-06 | 能量图及晶体位置查找表生成方法、装置、存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110211095A true CN110211095A (zh) | 2019-09-06 |
CN110211095B CN110211095B (zh) | 2023-05-30 |
Family
ID=67785498
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910371305.6A Active CN110211095B (zh) | 2019-05-06 | 2019-05-06 | 能量图及晶体位置查找表生成方法、装置、存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110211095B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110680369A (zh) * | 2019-09-18 | 2020-01-14 | 东软医疗系统股份有限公司 | 能量编码表压缩方法、查询方法、装置及pet系统 |
CN111568453A (zh) * | 2020-05-25 | 2020-08-25 | 上海联影医疗科技有限公司 | 能量校正状态检测方法、装置、计算机设备和存储介质 |
CN111624644A (zh) * | 2020-06-11 | 2020-09-04 | 北京卫星环境工程研究所 | 一种三维位置灵敏探测器及其能量校正方法 |
CN111685785A (zh) * | 2020-06-29 | 2020-09-22 | 上海联影医疗科技有限公司 | Pet晶体位置查找表的校正方法、装置以及计算机设备 |
WO2021238929A1 (zh) * | 2020-05-25 | 2021-12-02 | 上海联影医疗科技股份有限公司 | 一种成像设备的状态检测方法和系统 |
JP7512440B2 (ja) | 2020-05-25 | 2024-07-08 | 上海聯影医療科技股▲ふん▼有限公司 | 画像形成装置の状態の検出方法及びシステム |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103914860A (zh) * | 2013-01-05 | 2014-07-09 | 苏州瑞派宁科技有限公司 | 晶体条位置查找表生成方法及装置 |
CN104035122A (zh) * | 2014-05-22 | 2014-09-10 | 沈阳东软医疗系统有限公司 | 一种能量值的校正方法及装置 |
CN104337531A (zh) * | 2013-07-25 | 2015-02-11 | 苏州瑞派宁科技有限公司 | 针对全数字pet系统的在线能量符合方法及系统 |
CN104599302A (zh) * | 2015-01-13 | 2015-05-06 | 上海联影医疗科技有限公司 | 获取pet晶体能量峰值及设定能量鉴频器的方法 |
US20150301201A1 (en) * | 2014-04-18 | 2015-10-22 | Siemens Medical Solutions Usa, Inc. | Method and Apparatus for Automatic Calibration Check of PET Scanner Using Intrinsic Background Radiation of Scintillator Crystals |
CN107133462A (zh) * | 2017-04-27 | 2017-09-05 | 上海联影医疗科技有限公司 | 数据处理方法、装置及设备 |
CN107242881A (zh) * | 2017-08-15 | 2017-10-13 | 中日友好医院 | 一种pet晶体位置查找表更新方法 |
CN107569249A (zh) * | 2017-08-25 | 2018-01-12 | 沈阳东软医疗系统有限公司 | 一种晶体能量校正方法和装置 |
-
2019
- 2019-05-06 CN CN201910371305.6A patent/CN110211095B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103914860A (zh) * | 2013-01-05 | 2014-07-09 | 苏州瑞派宁科技有限公司 | 晶体条位置查找表生成方法及装置 |
CN104337531A (zh) * | 2013-07-25 | 2015-02-11 | 苏州瑞派宁科技有限公司 | 针对全数字pet系统的在线能量符合方法及系统 |
US20150301201A1 (en) * | 2014-04-18 | 2015-10-22 | Siemens Medical Solutions Usa, Inc. | Method and Apparatus for Automatic Calibration Check of PET Scanner Using Intrinsic Background Radiation of Scintillator Crystals |
CN104035122A (zh) * | 2014-05-22 | 2014-09-10 | 沈阳东软医疗系统有限公司 | 一种能量值的校正方法及装置 |
CN104599302A (zh) * | 2015-01-13 | 2015-05-06 | 上海联影医疗科技有限公司 | 获取pet晶体能量峰值及设定能量鉴频器的方法 |
CN107133462A (zh) * | 2017-04-27 | 2017-09-05 | 上海联影医疗科技有限公司 | 数据处理方法、装置及设备 |
CN107242881A (zh) * | 2017-08-15 | 2017-10-13 | 中日友好医院 | 一种pet晶体位置查找表更新方法 |
CN107569249A (zh) * | 2017-08-25 | 2018-01-12 | 沈阳东软医疗系统有限公司 | 一种晶体能量校正方法和装置 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110680369A (zh) * | 2019-09-18 | 2020-01-14 | 东软医疗系统股份有限公司 | 能量编码表压缩方法、查询方法、装置及pet系统 |
CN111568453A (zh) * | 2020-05-25 | 2020-08-25 | 上海联影医疗科技有限公司 | 能量校正状态检测方法、装置、计算机设备和存储介质 |
WO2021238929A1 (zh) * | 2020-05-25 | 2021-12-02 | 上海联影医疗科技股份有限公司 | 一种成像设备的状态检测方法和系统 |
JP7512440B2 (ja) | 2020-05-25 | 2024-07-08 | 上海聯影医療科技股▲ふん▼有限公司 | 画像形成装置の状態の検出方法及びシステム |
CN111624644A (zh) * | 2020-06-11 | 2020-09-04 | 北京卫星环境工程研究所 | 一种三维位置灵敏探测器及其能量校正方法 |
CN111685785A (zh) * | 2020-06-29 | 2020-09-22 | 上海联影医疗科技有限公司 | Pet晶体位置查找表的校正方法、装置以及计算机设备 |
CN111685785B (zh) * | 2020-06-29 | 2022-07-29 | 上海联影医疗科技股份有限公司 | Pet晶体位置查找表的校正方法、装置以及计算机设备 |
Also Published As
Publication number | Publication date |
---|---|
CN110211095B (zh) | 2023-05-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110211095A (zh) | 能量图及晶体位置查找表生成方法、装置、存储介质 | |
Suetens | Fundamentals of medical imaging | |
US7734119B2 (en) | Method and system for progressive multi-resolution three-dimensional image reconstruction using region of interest information | |
Heußer et al. | Investigation of the halo-artifact in 68Ga-PSMA-11-PET/MRI | |
JP2022516316A (ja) | プラットホーム非依存の全身画像セグメント化のためのシステムおよび方法 | |
WO2010143082A1 (en) | Time-of-flight positron emission tomography reconstruction using image content generated event-by-event based on time-of-flight information | |
TW202207241A (zh) | 用於偵測及表徵化病變之基於人工智慧的影像分析系統與方法 | |
Jin et al. | Multiple U-Net-based automatic segmentations and radiomics feature stability on ultrasound images for patients with ovarian cancer | |
CN103536305B (zh) | 用于执行图像类型识别的系统和方法 | |
Liu et al. | Deep learning with noise‐to‐noise training for denoising in SPECT myocardial perfusion imaging | |
CN113748443A (zh) | 用于交互式调整核医学图像中的强度窗口设定的系统和方法 | |
CN105913398A (zh) | Ct重建图像处理的方法及装置 | |
CN104685539B (zh) | 数字图像处理方法以及摄影装置 | |
CN103037768A (zh) | 用于处理图像的方法和设备以及采用该设备的医学图像系统 | |
Grégoire et al. | Is PET-based treatment planning the new standard in modern radiotherapy? The head and neck paradigm | |
US20190197674A1 (en) | Feature-based image processing using feature images extracted from different iterations | |
CN109242920A (zh) | 物质分解方法、装置和系统 | |
CN106659465A (zh) | 超声成像系统和方法 | |
CN111105475B (zh) | 一种基于正交角度x光的骨骼三维重建方法 | |
CN103823956B (zh) | 确定造影剂使用参数的方法、装置以及成像设备 | |
CN106952273B (zh) | 医学图像中胰腺的分割方法及装置 | |
Qiao et al. | 18F-FDG PET/CT radiomics nomogram for predicting occult lymph node metastasis of non-small cell lung cancer | |
US20080073538A1 (en) | Application-driven optimization of acquisition and reconstruction of SPECT/PET projection data | |
Wang et al. | Deep learning-augmented head and neck organs at risk segmentation from CT volumes | |
Chen et al. | Sketch-based Volumetric Seeded Region Growing. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20230413 Address after: Room 308, No. 177-2 Chuangxin Road, Hunnan District, Shenyang City, Liaoning Province, 110167 Applicant after: Shenyang Zhihe Medical Technology Co.,Ltd. Address before: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant before: Shenyang Neusoft Medical Systems Co.,Ltd. |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |