CN109521459A - 一种射线在闪烁晶体中的击中点定位方法及其系统 - Google Patents

一种射线在闪烁晶体中的击中点定位方法及其系统 Download PDF

Info

Publication number
CN109521459A
CN109521459A CN201811496270.0A CN201811496270A CN109521459A CN 109521459 A CN109521459 A CN 109521459A CN 201811496270 A CN201811496270 A CN 201811496270A CN 109521459 A CN109521459 A CN 109521459A
Authority
CN
China
Prior art keywords
hit
detector module
detector
little
signal
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
Application number
CN201811496270.0A
Other languages
English (en)
Other versions
CN109521459B (zh
Inventor
杨永峰
章先鸣
桑子儒
梁栋
刘新
郑海荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201811496270.0A priority Critical patent/CN109521459B/zh
Publication of CN109521459A publication Critical patent/CN109521459A/zh
Priority to PCT/CN2019/122005 priority patent/WO2020114330A1/zh
Application granted granted Critical
Publication of CN109521459B publication Critical patent/CN109521459B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/20Measuring radiation intensity with scintillation detectors
    • G01T1/202Measuring radiation intensity with scintillation detectors the detector being a crystal

Landscapes

  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Nuclear Medicine (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明提供了一种射线在闪烁晶体中的击中点定位方法及其系统,涉及核探测器技术领域。通过依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块,然后依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置,再当所述参考击中点的位置位于所述击中探测器模块的内部区域时,将所述参考击中点的位置作为击中点的位置,并且当所述参考击中点的位置位于所述击中探测器模块的边缘区域时,利用击中探测器模块以及从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。本发明提供的射线在闪烁晶体中的击中点定位方法及其系统能提高射线在闪烁晶体中击中点的定位精度。

Description

一种射线在闪烁晶体中的击中点定位方法及其系统
技术领域
本发明涉及核探测器技术领域,具体而言,涉及一种射线在闪烁晶体中的击中点定位方法及其系统。
背景技术
目前,正电子发射断层成像(Positron Emission Tomography,简称PET)能够定量和无创地测量由正电子发射核素标记的生物分子在生物体内随时间的变化情况,被广泛用于肿瘤、心血管疾病和神经疾病等的早期诊断。
PET包含几十个或几百个PET探测器模块,PET探测器模块主要包含闪烁晶体和光探测器,PET探测器模块间彼此紧靠,围绕成平板形状或圆柱形状。
PET的工作原理为向人体注射含有正电子放射性核素标记的药物,正电子放射性核素湮灭时产生两个能量为511keV夹角在180°附近的伽马射线,伽马射线与PET探测器模块中的闪烁晶体发生相互作用,产生大量闪烁光子,光探测器能探测到闪烁光子的强度分布,通过闪烁光子的分布可计算射线的作用位置,即击中点的位置。画线连接两个伽马射线的击中点可以得到一条响应线,代表正电子发射核素位于这条线上;通过大量响应线最终能确定正电子核素的具体位置,进而定位病灶区域。
因此,准确定位病灶区域要求准确定位响应线,即准确定位射线在探测器中的击中点。通常,根据探测器模块的信号总幅度确定击中点位于哪一个PET探测器模块,然后只使用该探测器模块所产生的信号来计算射线的击中点。但是,当射线击中某一个PET探测器模块的边缘区域时,闪烁光子可能只有一部分被该PET探测器模块接收到,而另一部分则被该PET探测器模块相邻的模块所采集,闪烁光子收集不完全可能导致击中点的定位不准确。
更具体的,对于基于切割晶体阵列的PET探测器,现有技术中存在切割晶体阵列的边缘晶条分辨不清楚的问题,其中切割晶体阵列由多个相同的长方体晶条组成,晶条之间按行列整齐排列。对于基于连续晶体的PET探测器,现有技术中存在连续晶体边缘位置分辨率变差的问题,其中连续晶体是内部未经切割的整块晶体。对于基于半连续晶体的PET探测器,现有技术中存在半连续晶体沿其连续方向的晶体边缘区域处位置分辨率差的问题。半连续晶体是由若干个连续晶体薄片组成的晶体,其中薄片沿PET扫描仪的轴向方向是连续的,沿PET扫描仪的径向方向是切割的,即由若干个薄片组成。
有鉴于此,如何解决上述问题,是本领域技术人员关注的重点。
发明内容
本发明的第一目的在于提供一种射线在闪烁晶体中的击中点定位方法,以解决现有技术中在晶体边缘击中点定位不够准确的问题。
本发明的第二目的在于提供一种用于确定击中点位置的电子系统,以解决现有技术中在晶体边缘击中点定位不够准确的问题。
本发明的第三目的在于提供一种可读写存储介质,以解决现有技术中生成的击中点的位置可能不准确的问题。
为了实现上述目的,本发明实施例采用的技术方案如下:
第一方面,本发明实施例提供了一种射线在闪烁晶体中的击中点定位方法,应用于核技术成像系统,所述核技术成像系统包括多个探测器模块,每个所述探测器模块均包括闪烁晶体,所述射线在闪烁晶体中的击中点定位方法包括以下步骤:
依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块;
依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置;
当所述参考击中点的位置位于所述击中探测器模块的内部区域时,将所述参考击中点的位置作为击中点的位置。
当所述参考击中点的位置位于所述击中探测器模块的边缘区域时,利用击中探测器模块以及每个从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。
第二方面,本发明实施例还提供了一种用于确定击中点位置的电子系统,应用于核技术成像系统,所述核技术成像系统包括多个探测器模块,每个所述探测器模块均包括闪烁晶体,所述用于确定击中点位置的电子系统包括:
击中探测器模块确定单元,用于依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块;
位置计算单元,用于依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置;
击中点确定单元,用于当所述参考击中点的位置位于所述击中探测器模块的内部区域时,将所述参考击中点的位置作为击中点的位置。
击中点确定单元还用于当所述参考击中点的位置位于所述击中探测器模块的边缘区域时,利用击中探测器模块以及每个从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。
第三方面,本发明实施例还提供了一种可读写存储介质,其上存储有运行程序,该运行程序被处理器执行时实现上述方法。
相对现有技术,本发明具有以下有益效果:
本发明提供了一种射线在闪烁晶体中的击中点定位方法及其系统,其中,通过依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块,然后依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置,当所述参考击中点的位置位于所述击中探测器模块的内部区域时,将所述参考击中点的位置作为击中点的位置,当所述参考击中点的位置位于所述击中探测器模块的边缘区域时,依据所述参考击中点的位置,所有探测器模块的安装位置,确定所有靠近所述参考击中点的相邻探测器模块;依次将所有所述相邻探测器模块的信号总幅度与预设定的信号幅度阈值进行比较,将信号总幅度大于所述预设定的信号幅度阈值的相邻探测器模块确定为从探测器模块。利用击中探测器模块以及从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。由于本发明对作用在探测器模块边缘区域和内部区域的事件分别采用不同的定位方法,当射线作用在探测器边缘区域时,同时利用击中探测器和从探测器各自产生的信号进行定位,从而能获得更完整的光分布信息,从而提高了射线击中探测器边缘区域时的定位精度,进而整体上使得击中点的定位精度更高。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1示出了本发明实施例提供的射线在闪烁晶体中的击中点确定方法的流程图。
图2示出了本发明实施例提供的射线在闪烁晶体上的信号幅度分布图。
图3示出了本发明实施例提供的基于切割晶体阵列的PET探测器的结构示意图。
图4示出了本发明实施例提供的基于半连续晶体的PET探测器的结构示意图。
图5示出了本发明实施例提供的图1中步骤S105的子步骤的流程图。
图6示出了本发明实施例提供的用于确定击中点位置的电子系统的模块示意图。
图7示出了本发明实施例提供的击中点确定单元的子模块示意图。
图标:11-闪烁晶体;12-耦合件;13-光导;14-光探测器;100-用于确定击中点位置的电子系统;110-击中探测器模块确定单元;120-位置计算单元;130-判断单元;140-击中点确定单元;141-从探测器模块确定模块;142-击中点确定模块。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。同时,在本发明的描述中,还需要说明的是,除非另有明确的规定和限定,术语“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。下面结合附图,对本发明的一些实施方式作详细说明。在不冲突的情况下,下述的实施例及实施例中的特征可以相互组合。
请参照图1,本发明实施例提供了一种射线在闪烁晶体中的击中点定位方法,应用于核技术成像系统,该核技术成像系统含有多个探测器模块,其中,探测器模块指可以独立探测并定位入射射线的最小单元,且每个探测器模块均含有闪烁晶体11与光探测器14,一般还包含光导13。在本实施例中,对于闪烁晶体11的材料不作限定,可以是任何一种有机或无机闪烁晶体,例如硅酸镥晶体;对于闪烁晶体11的宏观结构不作限定,可以是连续晶体,切割晶体阵列,半连续晶体等任何一种能实现伽马射线探测的晶体宏观结构;光探测器14的种类也不作限定,例如硅光电倍增管(SiPM,Silicon Photomultipliers,or MPPCsMulti-Pixel Photon Counters);其通常包含大于等于1个的读出单元数目,每个读出单元均可以独立输出信号,读出单元也称为像素。该射线在闪烁晶体中的击中点定位方法包括:
步骤S11,依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块。
在本实施例中,击中探测器模块为与能量为511keV的伽马射线发生相互作用的模块。根据相互作用种类的不同,伽马射线在探测器模块中产生的闪烁光子数目也不同;另外,射线探测是一个随机过程,其中包含统计涨落,所以511keV伽马射线击中某个探测器模块所产生的信号总幅度在一个范围内变化。本实施例确定击中探测器模块的方法为将每个所述探测器模块接收到的信号总幅度与预设定的幅度范围进行比较,将信号总幅度落在所述预设定的幅度范围内的探测器模块作为击中探测器;例如,预设定的幅度范围可以是350-650keV,也可以根据实际情况选择其他合适的范围,本实施例对此并不做任何限定。
在本实施例中,由于单位时间内放射性核素的衰变个数可能大于一个,因此可能出现单位时间内不止一个511keV伽马射线与探测器模块发生相互作用并产生信号,因此击中探测器模块的数量可以不止一个,例如,击中探测器模块的数量为两个。
步骤S12,依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置。
在确定了击中探测器模块后,可依据该击中探测器模块接收到的信号与预设定的算法计算参考击中点的位置。
其中,本实施例所述的预设定的算法包括但不限于重心法,加权重心法,最小二乘法,最大似然法,神经网络等方法。
下面以重心法为例进行说明,每个探测器模块均包括多个像素,且多个像素按阵列相邻排布,首先,依据公式确定所述参考击中点在水平面x方向上的位置,其中,n表示光探测器14的像素个数,Ai表示第i个像素的信号幅度,Ci表示第i个像素在x方向上的位置坐标,x表示参考击中点在水平面x方向上的位置坐标。然后,依据公式确定所述参考击中点在水平面y方向上的位置,其中,n表示光探测器14的像素个数,Ai表示第i个像素的信号幅度,di表示第i个像素在y方向上的位置坐标,y表示参考击中点在水平面y方向上的位置坐标。依据每个像素的信号幅度可确定(x,y)坐标,即确定所述参考击中点在水平面中的位置。
步骤S13,判断所述参考击中点的位置是否在所述击中探测器模块的内部区域,如果是,则执行步骤S14,如果否,则执行步骤S15。
在本实施例中,每个探测器模块均分为内部区域与边缘区域,其中,边缘区域指在水平面内,闪烁晶体11中距离其边界小于阈值的区域,内部区域为闪烁晶体11中除边缘区域以外的区域,且水平面与光探测器14的法线垂直,例如,对于一个50mm×50mm×10mm的闪烁晶体11,阈值可能设置为10mm,即边缘区域为水平面x或y方向上距离闪烁晶体11的边界小于10mm的晶体区域。
同时,请参阅图2,其中,图2中a图为现有技术仅使用单个探测器模块的信号确定击中点的说明。当入射射线在探测器中发生相互作用,并以4π立体角发射大量闪烁光子时,由于参考击中点位于边缘区域,朝击中点右侧发射的闪烁光子大部分在晶体边界被反射回晶体中间,导致光分布畸变;小部分逃逸到相邻探测器中,仅有很小部分未经反射且被该光探测器14正常收集。可以明白的,此时,单个探测器模块只能探测到一部分的光分布信息,反射的闪烁光子导致光子分布的重心位于实际击中点的左侧,从而使定位产生误差。图2中b为实施例提供的定位方法,在探测器模块之间设置有耦合件12,用来减弱闪烁光子在晶体边界的反射,耦合件12的材料不作限定,例如空气,或硅油。当参考击中点位于边缘区域时,有一部分闪烁光能逃逸到相邻探测器中并被探测到。联合击中探测器和其相邻(从)探测器的信号能探测到更完整的光分布信息,从而光子分布重心与实际击中点更接近,从而定位击中点的误差较小。
有鉴于此,在本实施例中,需进行参考击中点的位置的判断。
步骤S14,将所述参考击中点的位置作为击中点的位置。
当参考击中点的位置在击中探测器模块的内部区域时,即表示射线的信号已大部分被该击中探测器模块所采集,从探测器含有的光子分布信息较少,可能与噪声水平相当。因此,此时利用该击中探测器模块即可准确计算出击中点的位置。
步骤S15,利用击中探测器模块以及从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。
当在参考击中点的位置处于边缘区域时,单个探测器模块只能探测到一部分的光分布信息,此时通过击中探测器模块与从探测器模块的信号进行对击中点的重新定位,从而更加准确地定位击中点,进而提升成像系统的空间分辨率。
具体地,请参阅图3,步骤S15包括:
子步骤S151,依据所述参考击中点的位置,所有探测器模块的安装位置,每个探测器模块接收到的信号总幅度以及预设定的信号幅度阈值确定所述从探测器模块。
其中,需要说明的是,在本实施例中,从探测器模块指紧挨击中探测器模块,同样收集到射线的信号,且信号较大的探测器模块。
其中,子步骤S151包括:
子步骤S1511,依据所述所有探测器模块的安装位置确定靠近所述参考击中点的位置的相邻探测器模块。
由于信号集中于射线的击中点附近,因此从探测器模块需要选择与参考击中点相邻的探测器模块。
本实施例中,所有的探测器模块按相邻的方式排列,即大多数探测器模块的四周均设置有探测器模块,可以理解地,当该参考击中点的位置确定之后,即可确定参考击中点更靠近击中探测器模块的哪一侧,进而根据所有探测器模块的安装位置找到与该侧相邻的探测器模块。其中,当参考击中点位于晶体的角落,可能有三个相邻探测器模块;当位于晶体某个边界的正中间时,可能只有一个相邻探测器模块,因此相邻的探测器模块可以是一个,也可以是多个。
子步骤S1512,依次将所述相邻探测器模块中的每个探测器模块的信号总幅度与预设定的信号幅度阈值进行比较,将大于所述预设定的信号幅度阈值的探测器模块确定为从探测器模块。
如果从探测器信号太小,那么探测器自身噪声相对于信号来说则过大,因此要选择信号较大的相邻探测器模块作为从探测器模块。确定所述预设定的信号幅度阈值需要事先估计探测器模块内的噪声水平。进一步地,在本实施例中,通过依次将相邻探测器模块中的每个探测器模块的信号总幅度与预设定的信号幅度进行比较的方式,例如,当阈值设置为70keV时,信号总幅度大于70keV的相邻探测器模块就是从探测器模块。
由于需要通过比较的方式确定从探测器模块的数量,因此可能出现每个相邻的探测器模块的信号总幅度均小于阈值的情况,即可能出现从探测器模块数量为零的情况,例如,相邻的探测器模块的数量为3,当从探测器模块的数量为0;当然地,从探测器模块的数量也可以为多个,本实施例对此并不做任何限定。
步骤S15还包括:
子步骤S152,通过所述预设定的算法,用所述击中探测器模块以及所述从探测器模块各自产生的的信号对击中点重新定位。
确定从探测器模块后,利用从探测器模块与击中探测器模块的各自产生的信号对击中点进行重新定位,其中,预设定的算法与步骤S102中所述的算法相同,即包括但不限于重心法,加权重心法,最小二乘法,最大似然法及以这些方法为基础的改进方法。计算过程中,为了提高信号的信噪比,可以仅使用从探测器模块的部分信号,例如最靠近击中探测器的那部分像素的信号;也可以使用从探测器的所有信号进行定位。
此外,本方法可以用于任何基于闪烁晶体11的,含有多个探测器模块的探测器面阵或核技术成像系统,例如X/伽马射线照相机,计算机断层扫描(Computed Tomography),正电子断层成像扫描仪等。进一步地,对于正电子断层成像扫描仪,请参阅图4,本方法可以用在基于切割晶体阵列的PET探测器来实施,用来解决现有技术中切割晶体阵列的边缘晶条分辨不清楚的问题,其中切割晶体阵列由多个相同的长方体晶条组成,晶条之间按行列整齐排列。本方法还可以用在基于连续晶体的PET探测器中,用来解决现有技术中连续晶体边缘位置分辨率变差的问题,其中连续晶体是内部未经切割的整块晶体。请参阅图5,本方法还可以用在基于半连续晶体的PET探测器中,用来解决现有技术中半连续晶体沿其连续方向的晶体边缘区域处位置分辨率差的问题。半连续晶体是由若干个连续晶体薄片组成的晶体,其中薄片沿PET扫描仪的轴向方向是连续的,沿PET扫描仪的径向方向是切割的,即由若干个薄片组成。
第二实施例
请参阅图6,本发明实施例还提供了一种确定击中点位置的电子系统,应用于核技术成像系统,该核技术成像系统多个探测器模块,其中,探测器模块指可以独立探测并定位入射线的最小单元,且每个探测器模块均包括闪烁晶体11与光探测器14,该用于确定击中点位置的电子系统100包括:
击中探测器模块确定单元110,用于依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块。
可以理解的,通过击中探测器模块确定单元110可执行步骤S11。
位置计算单元120,用于依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置。
可以理解的,通过位置计算单元120可执行步骤S12。
判断单元130,用于判断所述参考击中点的位置是否在所述击中探测器模块的内部区域。
可以理解的,通过判断单元130可执行步骤S13。
击中点确定单元140,用于将所述参考击中点的位置作为击中点的位置。
可以理解的,通过击中点确定单元140可执行步骤S14。
击中点确定单元140还用于利用击中探测器模块以及从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。
可以理解的,通过击中点确定单元140可执行步骤S15。
其中,请参阅图7,所述击中点确定单元140包括:
从探测器模块确定模块141,用于依据所述参考击中点的位置,所有探测器模块的安装位置,每个探测器模块接收到的信号总幅度以及预设定的信号幅度阈值确定所述从探测器模块。
可以理解的,通过从探测器模块确定模块141可执行子步骤S151。
具体地,从探测器模块确定模块141用于依据所述所有探测器模块的安装位置确定靠近所述参考击中点的位置的相邻探测器模块,并依次将所述相邻探测器模块中的每个探测器模块的信号总幅度与预设定的信号幅度阈值进行比较,将大于所述预设定的信号幅度阈值的探测器模块确定为从探测器模块。
击中点确定模块142,用于通过所述预设定的算法,利用所述击中探测器模块以及所述从探测器模块的信号,对击中点重新定位。
可以理解的,通过击中点确定模块142可执行子步骤S152。
第三实施例
本发明实施例还揭示了一种可读写存储介质,该运行程序被处理器执行时实现本发明前述实施例揭示的射线在闪烁晶体11中的击中点定位方法。
综上所述,本发明提供了一种射线在闪烁晶体中的击中点定位方法及其系统,其中,通过依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块,然后依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置,再当所述参考击中点的位置位于所述击中探测器模块的内部区域时,将所述参考击中点的位置作为击中点的位置,并且当所述参考击中点的位置位于所述击中探测器模块的边缘区域时,利用击中探测器模块以及从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。由于本发明提供的射线在闪烁晶体中的击中点定位方法及其系统能够在确定集中探测器模块后,依据参考击中点的位置判断该位置在击中探测器模块的内部或者边缘,同时在参考击中点的位置处于击中探测器的中间位置时,将其作为击中点的位置,并在其处于击中探测器的边缘位置时,能够利用集中探测器和从探测器模块对击中点进行重新定位。通过判断参考击中点的位置,能够实现分情况的确定击中点的位置,使得击中点的定位精度更高。
需要说明的是,在本文中,诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。

Claims (10)

1.一种射线在闪烁晶体中的击中点定位方法,应用于核技术成像系统,其特征在于,所述核技术成像系统包括多个探测器模块,每个所述探测器模块均包括闪烁晶体,所述射线在闪烁晶体中的击中点定位方法包括以下步骤:
依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块;
依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置;
当所述参考击中点的位置位于所述击中探测器模块的内部区域时,将所述参考击中点的位置作为击中点的位置;
当所述参考击中点的位置位于所述击中探测器模块的边缘区域时,利用击中探测器模块以及每个从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。
2.如权利要求1所述的射线在闪烁晶体中的击中点定位方法,其特征在于,所述利用击中探测器模块以及每个从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位的步骤包括:
依据所述参考击中点的位置,所有探测器模块的安装位置,每个探测器模块接收到的信号总幅度以及预设定的信号幅度阈值确定所述从探测器模块;
通过所述预设定的算法,用所述击中探测器模块以及所述从探测器模块各自产生的信号对击中点重新定位。
3.如权利要求2所述的射线在闪烁晶体中的击中点定位方法,其特征在于,所述依据所述参考击中点的位置,所有探测器模块的安装位置,每个探测器模块接收到的信号总幅度以及预设定的信号幅度阈值确定所述从探测器模块的步骤包括:
依据所述所有探测器模块的安装位置确定靠近所述参考击中点的位置的相邻探测器模块;
依次将所述相邻探测器模块中的每个探测器模块的信号总幅度与预设定的信号幅度阈值进行比较,将大于所述预设定的信号幅度阈值的探测器模块确定为从探测器模块。
4.如权利要求1所述的射线在闪烁晶体中的击中点定位方法,其特征在于,所述依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块的步骤包括:
将每个所述探测器模块接收到的信号总幅度与预设定的幅度范围进行比较,并将信号总幅度落在所述预设定的幅度范围内的探测器模块作为击中探测器。
5.如权利要求1所述的射线在闪烁晶体中的击中点定位方法,其特征在于,所述边缘区域指在水平面内,闪烁晶体中距离其边界小于阈值的区域,所述水平面与光探测器的法线垂直,所述内部区域为所述闪烁晶体中除所述边缘区域以外的区域。
6.如权利要求1所述的射线在闪烁晶体中的击中点定位方法,其特征在于,所述探测器模块还包括光探测器,所述光探测器的读出单元数目大于或等于1个,且每个读出单元均可独立输出信号。
7.一种用于确定击中点位置的电子系统,应用于核技术成像系统,其特征在于,所述核技术成像系统包括多个探测器模块,每个所述探测器模块均包括闪烁晶体,所述用于确定击中点位置的电子系统包括:
击中探测器模块确定单元,用于依据每个所述探测器模块分别接收到的信号总幅度确定击中探测器模块;
位置计算单元,用于依据所述击中探测器模块接收到的信号通过预设定的算法计算参考击中点的位置;
击中点确定单元,用于当所述参考击中点的位置位于所述击中探测器模块的内部区域时,将所述参考击中点的位置作为击中点的位置;
击中点确定单元还用于当所述参考击中点的位置位于所述击中探测器模块的边缘区域时,利用击中探测器模块以及每个从探测器模块各自产生的信号,通过预设定的算法对击中点重新定位。
8.如权利要求7所述的用于确定击中点位置的电子系统,其特征在于,所述击中点确定单元包括:
从探测器模块确定模块,用于当所述参考击中点的位置在预设定的边缘位置时,依据所述参考击中点的位置,所有探测器模块的安装位置,每个探测器模块接收到的信号总幅度以及预设定的信号幅度阈值确定所述从探测器模块;
击中点确定模块,用于通过所述预设定的算法,利用所述击中探测器模块以及每个所述从探测器模块各自产生的信号,对击中点重新定位。
9.如权利要求8所述的用于确定击中点位置的电子系统,其特征在于,所述从探测器模块确定模块用于依据所述所有探测器模块的安装位置确定靠近所述参考击中点的位置的相邻探测器模块,并依次将所述相邻探测器模块中的每个探测器模块的信号总幅度与预设定的信号幅度阈值进行比较,将大于所述预设定的信号幅度阈值的探测器模块确定为从探测器模块。
10.一种可读写存储介质,其上存储有运行程序,其特征在于,该运行程序被处理器执行时实现如权利要求1-6任一项所述的方法。
CN201811496270.0A 2018-12-07 2018-12-07 一种射线在闪烁晶体中的击中点定位方法及其系统 Active CN109521459B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201811496270.0A CN109521459B (zh) 2018-12-07 2018-12-07 一种射线在闪烁晶体中的击中点定位方法及其系统
PCT/CN2019/122005 WO2020114330A1 (zh) 2018-12-07 2019-11-29 一种射线在闪烁晶体中的击中点定位方法及其系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811496270.0A CN109521459B (zh) 2018-12-07 2018-12-07 一种射线在闪烁晶体中的击中点定位方法及其系统

Publications (2)

Publication Number Publication Date
CN109521459A true CN109521459A (zh) 2019-03-26
CN109521459B CN109521459B (zh) 2020-08-25

Family

ID=65796289

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811496270.0A Active CN109521459B (zh) 2018-12-07 2018-12-07 一种射线在闪烁晶体中的击中点定位方法及其系统

Country Status (2)

Country Link
CN (1) CN109521459B (zh)
WO (1) WO2020114330A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020114330A1 (zh) * 2018-12-07 2020-06-11 深圳先进技术研究院 一种射线在闪烁晶体中的击中点定位方法及其系统
WO2021258507A1 (zh) * 2020-06-24 2021-12-30 北京卫星环境工程研究所 连续晶体伽马探测器的三维位置刻度方法、装置及设备
CN117635614A (zh) * 2024-01-25 2024-03-01 南京航空航天大学 一种超快x射线成像无损检测方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030193029A1 (en) * 2002-04-16 2003-10-16 General Electric Company Edge resolved dual scintillator gamma ray detection system and method
CN1993632A (zh) * 2003-11-20 2007-07-04 杰克·E·朱尼 晶体边缘效应的处理
CN101833106A (zh) * 2010-05-11 2010-09-15 刘继国 一种用于射线位置和能量测量的闪烁体探测器
JP2011149883A (ja) * 2010-01-22 2011-08-04 Natl Inst Of Radiological Sciences 放射線位置検出器の位置演算方法及び装置
JP2017003392A (ja) * 2015-06-09 2017-01-05 株式会社島津製作所 位置マップ作成方法および核医学診断装置
US20170108593A1 (en) * 2015-10-20 2017-04-20 Institute of Nuclear Energy Research, Atomic Energy Council, Executive Yuan, R.O.C. Device for determining three dimensional locations and energy of gamma incidence events and method for the same
CN107970037A (zh) * 2016-10-24 2018-05-01 上海东软医疗科技有限公司 成像方法和成像系统
CN108051843A (zh) * 2017-10-30 2018-05-18 深圳先进技术研究院 一种pet探测器的位置编码电路、方法及装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8080780B2 (en) * 2009-12-29 2011-12-20 Kabushiki Kaisha Toshiba Apparatus and associated methodology for improving timing resolution in gamma ray detection
CN101839991B (zh) * 2010-05-06 2013-06-12 清华大学 一种复合光敏器件斜排列式高能射线探测器
WO2014020471A2 (en) * 2012-07-30 2014-02-06 Koninklijke Philips N.V. High spatial resolution mode solid state positron emission tomography (pet) scanner
CN104035122A (zh) * 2014-05-22 2014-09-10 沈阳东软医疗系统有限公司 一种能量值的校正方法及装置
WO2016062799A1 (en) * 2014-10-23 2016-04-28 Petsys Electronics, Medical Pet Detectors, Sa Detector component for an x-ray or gamma ray detector
CN107250844B (zh) * 2014-12-12 2020-06-09 爱丁堡大学董事会 用于处理正电子发射断层摄影数据的方法和装置
CN104932000B (zh) * 2015-05-21 2017-11-28 沈阳东软医疗系统有限公司 一种PET设备中γ光子位置信息获取的方法及装置
CN106556857B (zh) * 2016-10-19 2019-04-02 武汉中派科技有限责任公司 反应位置定位方法和装置
JP6903435B2 (ja) * 2017-01-11 2021-07-14 キヤノンメディカルシステムズ株式会社 核医学診断装置及びキャリブレーション方法
CN109521459B (zh) * 2018-12-07 2020-08-25 深圳先进技术研究院 一种射线在闪烁晶体中的击中点定位方法及其系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030193029A1 (en) * 2002-04-16 2003-10-16 General Electric Company Edge resolved dual scintillator gamma ray detection system and method
CN1993632A (zh) * 2003-11-20 2007-07-04 杰克·E·朱尼 晶体边缘效应的处理
JP2011149883A (ja) * 2010-01-22 2011-08-04 Natl Inst Of Radiological Sciences 放射線位置検出器の位置演算方法及び装置
CN101833106A (zh) * 2010-05-11 2010-09-15 刘继国 一种用于射线位置和能量测量的闪烁体探测器
JP2017003392A (ja) * 2015-06-09 2017-01-05 株式会社島津製作所 位置マップ作成方法および核医学診断装置
US20170108593A1 (en) * 2015-10-20 2017-04-20 Institute of Nuclear Energy Research, Atomic Energy Council, Executive Yuan, R.O.C. Device for determining three dimensional locations and energy of gamma incidence events and method for the same
CN107970037A (zh) * 2016-10-24 2018-05-01 上海东软医疗科技有限公司 成像方法和成像系统
CN108051843A (zh) * 2017-10-30 2018-05-18 深圳先进技术研究院 一种pet探测器的位置编码电路、方法及装置

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020114330A1 (zh) * 2018-12-07 2020-06-11 深圳先进技术研究院 一种射线在闪烁晶体中的击中点定位方法及其系统
WO2021258507A1 (zh) * 2020-06-24 2021-12-30 北京卫星环境工程研究所 连续晶体伽马探测器的三维位置刻度方法、装置及设备
CN117635614A (zh) * 2024-01-25 2024-03-01 南京航空航天大学 一种超快x射线成像无损检测方法及系统
CN117635614B (zh) * 2024-01-25 2024-04-02 南京航空航天大学 一种超快x射线成像无损检测方法及系统

Also Published As

Publication number Publication date
CN109521459B (zh) 2020-08-25
WO2020114330A1 (zh) 2020-06-11

Similar Documents

Publication Publication Date Title
CN105452902B (zh) 针对小间距闪烁晶体阵列的像素识别
CN105190357B (zh) 用于对伽马辐射事件的评估的装置和方法
US9989658B2 (en) Method of, and apparatus for, processing positron emission tomography data
CN107076859B (zh) 具有光共享和相互作用深度估计的pet探测器闪烁体布置
JP5174813B2 (ja) Tofpetの再構成を改良するための方法及びシステム
CN109521459A (zh) 一种射线在闪烁晶体中的击中点定位方法及其系统
CN102755172B (zh) 核医学成像方法以及核医学成像装置
CN106539591A (zh) Pet飞行时间状态质量检测方法和pet扫描装置
Okada et al. Development of a gamma camera to image radiation fields
Stockhoff et al. Optical simulation study on the spatial resolution of a thick monolithic PET detector
CN105496436B (zh) 用于pet装置的时间校正方法和装置
CN102331586A (zh) 电荷丢失校正
CN103998952A (zh) X射线探测器
EP2594961A2 (en) Surface contamination monitoring system and method
US20140209804A1 (en) Apparatuses for and methods of generating images
US9864073B1 (en) Systems and methods for controlling optical coupling between scintillator crystals
CN105974461A (zh) 一种γ射线扫描成像及核素识别系统及其方法
Pani et al. A Novel Method for $\gamma-\text {photons} $ Depth-of-Interaction Detection in Monolithic Scintillation Crystals
TWI400470B (zh) 決定晶體穿透效應與偵檢效率以及偵檢效率校正之方法
McIntosh et al. Validation of a GATE Model of $^{176} $ Lu Intrinsic Radioactivity in LSO PET Systems
CN107884805B (zh) 一种射源定位的方法及装置
Brambilla et al. Impact of target‐to‐background ratio, target size, emission scan duration, and activity on physical figures of merit for a 3D LSO‐based whole body PET/CT scanner
WO2019183594A1 (en) Calibration method and system for photon or particle counting detectors
Wang et al. Novel light-guide-PMT geometries to reduce dead edges of a scintillation camera
CA2425244C (en) Method for reduction of the statistical measurement times in the field of radioactivity measurement

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