CN114895344A - 空间编码探测器、空间编码、解码方法、装置及存储介质 - Google Patents

空间编码探测器、空间编码、解码方法、装置及存储介质 Download PDF

Info

Publication number
CN114895344A
CN114895344A CN202210372854.7A CN202210372854A CN114895344A CN 114895344 A CN114895344 A CN 114895344A CN 202210372854 A CN202210372854 A CN 202210372854A CN 114895344 A CN114895344 A CN 114895344A
Authority
CN
China
Prior art keywords
crystal
spatial
array
ray
different
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.)
Pending
Application number
CN202210372854.7A
Other languages
English (en)
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.)
Raycan Technology Co Ltd
Original Assignee
Raycan Technology Co Ltd
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 Raycan Technology Co Ltd filed Critical Raycan Technology Co Ltd
Priority to CN202210372854.7A priority Critical patent/CN114895344A/zh
Publication of CN114895344A publication Critical patent/CN114895344A/zh
Priority to PCT/CN2022/136854 priority patent/WO2023197626A1/zh
Pending legal-status Critical Current

Links

Images

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
    • 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/208Circuits specially adapted for scintillation detectors, e.g. for the photo-multiplier section

Landscapes

  • 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)
  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本申请提出一种空间编码探测器、空间编码、解码方法、装置及存储介质,空间编码探测器包括闪烁晶体阵列以及光电转换器,闪烁晶体阵列包括至少两根晶体条,其中至少两根晶体条具有不同的衰减时间;光电转换器与晶体阵列耦合;空间编码方法包括:建立闪烁晶体阵列的空间位置编码表;根据空间分辨率需求从空间位置编码表中选取对应规格的晶体条组成闪烁晶体阵列;根据预设编码顺序或者预设编码位置耦合闪烁晶体阵列与光电转换器。本申请通过将闪烁晶体阵列中的晶体条进行空间编码,能够通过衰减时间解码出额外空间信息,从而突破光电转换器件阵列对探测器空间分辨率的限制。

Description

空间编码探测器、空间编码、解码方法、装置及存储介质
技术领域
本申请涉及射线探测领域,具体而言,涉及一种空间编码探测器、空间编码、解码方法、装置及存储介质。
背景技术
射线通常是指能量不低于100eV的光子束流,包括X射线、γ射线、α射线、β射线、中子射线以及质子射线等。射线探测器对射线能量信息、时间信息、空间信息的分辨能力直接决定了探测系统的成像质量。如图1所示,射线探测器的工作原理如下:首先,射线与闪烁晶体相互作用将射线转换为可见光光子,可见光光子在闪烁晶体内输运并入射到与闪烁晶体耦合的光电转换器件中。光电转换器件将入射可见光转换为电信号,利用与光电转换器件匹配的电子学系统采集并输出数字信号。然后利用软件算法,即可从数字信号中计算射线的时间、能量、位置等信息。
射线探测应用广泛,比如,正电子发射断层成像(Position EmissionTomography,简称PET)使用伽马射线探测器作为系统前端获取正电子湮灭事件的分布信息,从而辅助进行癌症诊疗、脑科学研究、心脏病学研究、重离子放疗监测等。
基于PET对时间性能的要求,目前得到应用的射线探测器基本以闪烁晶体和光电转换器件为基本结构,其中,闪烁晶体典型的有LYSO、BGO、YSO,光电转换器件包括位敏型PMT、SiPM、SPAD。
目前,常见的晶体探测器包括单层阵列单一衰减时间闪烁晶体阵列探测器、单层阵列闪烁晶体+光导层探测器、双端读出的单层单一衰减时间闪烁晶体探测器以及组合闪烁晶体阵列探测器。
单层阵列单一衰减时间闪烁晶体探测器是将多根衰减时间相同的晶体条组合成闪烁晶体阵列,晶体条之间和晶体条与外界之间使用反光材料,例如硫酸钡涂层或工艺反射膜建立完全光学隔离。晶体条一端的光学面不覆盖反光材料,通过硅脂等耦合剂与光电转换器件阵列耦合。
单层阵列闪烁晶体及光导层探测器是在闪烁晶体阵列与光电转换器件阵列之间增加了一层对可见光光子透明的介质,即光导。可见光光子射出闪烁晶体后经过一个发散的过程再射入光电转换器件阵列。
双端读出的单层单一衰减时间闪烁晶体探测器是在晶体条两端的光学面都不覆盖反光材料,通过硅脂等光学耦合剂与光电转换器阵列耦合。
组合闪烁晶体探测器是将多层闪烁晶体阵列和多层连续闪烁晶体耦合在一起,再与光电转换器件阵列耦合。
对于单层晶体阵列来说,在投影平面上的空间分辨率受制于光电转换器件尺寸。而组合闪烁晶体探测器又存在结构复杂、加工难度大及误码率高等问题。
发明内容
本申请提出一种空间编码探测器、空间编码、解码方法、装置、电子设备及存储介质,以解决探测器中光电转换器件尺寸对空间分辨率的限制问题。
根据本申请的一方面,提出一种空间编码探测器,包括:闪烁晶体阵列以及光电转换器,闪烁晶体阵列包括至少两根晶体条,其中至少两根所述晶体条具有不同的衰减时间;所述光电转换器与所述晶体阵列耦合。
根据一些实施例,所述闪烁晶体阵列中包括至少一个晶体单元,所述晶体单元中包括多根所述晶体条。
根据一些实施例,所述晶体单元中至少两根所述晶体条具有不同的衰减时间,或者所述晶体单元中至少一部分所述晶体条的衰减时间与至少一根所述晶体条的衰减时间不同。
根据一些实施例,所述至少一个晶体单元的大小相同,不同所述晶体单元之间至少一部分所述晶体条的衰减时间不同。
根据一些实施例,所述闪烁晶体阵列包括至少两种大小的晶体单元,不同所述晶体单元之间至少一部分所述晶体条的衰减时间不同。
根据一些实施例,所述晶体单元中的所述晶体条呈矩阵排列,不同所述晶体单元之间对应位置的所述晶体条具有相同的衰减时间。
根据一些实施例,每根所述晶体条的衰减时间均不同。
根据一些实施例,所述晶体条的衰减时间呈等差数列。
根据一些实施例,所述晶体条的衰减时间介于30ns~50ns之间。
根据一些实施例,其中一部分所述晶体条的外形尺寸不同。
根据一些实施例,所述晶体条的外形包括棱柱、圆柱或者不规则形状。
根据一些实施例,所述光电转换器为光电转换器件阵列,所述光电转换器件阵列中的每个单元与所述晶体条之间为一对一、一对多或者多对一的形式耦合。
根据本申请的一方面,提出一种空间编码方法,所述空间编码方法包括:建立闪烁晶体阵列的空间位置编码表,所述空间位置编码表包括不同衰减时间、不同尺寸的晶体条与空间分辨率之间的对应关系;根据空间分辨率需求从所述空间位置编码表中选取对应规格的所述晶体条组成所述闪烁晶体阵列;根据预设编码顺序或者预设编码位置耦合所述闪烁晶体阵列与光电转换器。
根据一些实施例,所述空间位置编码表还包括不同位置、不同尺寸或者不同衰减时间的晶体条按照预设编码顺序或者预设编码位置排列后与空间分辨率之间的对应关系。
根据一些实施例,同一所述空间分辨率对应多种形式的所述预设编码顺序或者所述预设编码位置。
根据一些实施例,所述光电转换器为光电转换器件阵列,所述光电转换器件阵列中的每个单元与所述晶体条之间为一对一、一对多或者多对一的形式耦合
根据本申请的一方面,提出一种空间解码方法,所述空间解码方法用于如上所述的空间编码探测器,所述空间解码方法包括:利用光电转换器输出的脉冲信号获取所述射线的信息,所述信息包括能量信息、到达时间和衰减时间;利用所述能量信息和所述到达时间确定所述射线沉积的晶体单元;利用所述衰减时间和预先建立的查找表获取所述射线的沉积位置,所述查找表包括所述衰减时间与射线沉积位置的对应关系。
根据一些实施例,所述利用光电转换器输出的脉冲信号获取所述射线的信息,包括:利用多电压阈值采样方法、示波器或模数转换器获取所述射线的信息。
根据一些实施例,所述确定射线沉积的晶体单元,包括:根据所述闪烁晶体阵列与所述光电转换器的耦合形式利用一对一耦合法或者重心法确定所述射线沉积的晶体单元。
根据一些实施例,所述查找表通过如下方式预先建立:利用不同规格的所述晶体条探测所述射线;计算每根所述晶体条中沉积的所述射线的衰减时间;利用所计算的衰减时间和所述射线在不同晶体条中的沉积位置建立所述查找表。
根据本申请的一方面,提出一种空间解码装置,所述空间解码装置用于如上所述的空间编码探测器,所述空间解码装置包括:射线信息获取单元,用于利用光电转换器输出的脉冲信号获取射线的信息,所述信息包括能量信息、到达时间和衰减时间;晶体区域确定单元,用于利用所述能量信息和所述到达时间确定所述射线沉积的晶体单元;沉积位置确定单元,用于利用所述衰减时间和预先建立的查找表获取所述射线的沉积位置,所述查找表包括所述衰减时间与所述沉积位置的对应关系。
根据本申请的一方面,提出一种空间解码装置,包括:一个或多个处理器;存储装置,用于存储计算机程序;当所述计算机程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现如上所述的方法。
根据本申请的一方面,提出一种存储介质,其上存储有程序指令,所述程序指令被执行时实现如上所述的方法。
在计算射线的沉积位置时,由于光电转换器件阵列的空间分辨率限制探测器分辨沉积所在位置的能力,本申请通过将闪烁晶体阵列中晶体条的衰减时间进行空间编码,并通过衰减时间解码出的额外空间信息进一步定位具体发生沉积的晶体条,从而能够突破光电转换器件阵列空间分辨率对探测器空间分辨率的限制。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍。
图1示出一种射线探测器的工作原理图。
图2a示出根据本申请示例实施例的一种空间编码探测器的侧面示意图。
图2b示出根据本申请示例实施例的一种空间编码探测器的横截面图。
图3示出根据本申请示例实施例的一种闪烁晶体阵列的空间编码示意图。
图4示出根据本申请示例实施例的另一种闪烁晶体阵列的空间编码示意图。
图5示出根据本申请示例实施例的一种空间编码方法的流程图。
图6示出根据本申请示例实施例的一种空间解码方法的流程图。
图7示出根据本申请示例实施例的一种建立查找表的方法流程图。
图8示出根据本申请示例实施例的一种空间解码装置的示意框图。
图9示出了根据本申请实施例的又一种空间解码装置。
具体实施方式
现在将参考附图更全面地描述示例实施例。然而,示例实施例能够以多种形式实施,且不应被理解为限于在此阐述的实施例;相反,提供这些实施例使得本申请将全面和完整,并将示例实施例的构思全面地传达给本领域的技术人员。在图中相同的附图标记表示相同或类似的部分,因而将省略对它们的重复描述。
所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施例中。在下面的描述中,提供许多具体细节从而给出对本公开的实施例的充分理解。然而,本领域技术人员将意识到,可以实践本公开的技术方案而没有这些特定细节中的一个或更多,或者可以采用其它的方式、组元、材料、装置或操作等。在这些情况下,将不详细示出或描述公知结构、方法、装置、实现、材料或者操作。
附图中所示的流程图仅是示例性说明,不是必须包括所有的内容和操作/步骤,也不是必须按所描述的顺序执行。例如,有的操作/步骤还可以分解,而有的操作/步骤可以合并或部分合并,因此实际执行的顺序有可能根据实际情况改变。
本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别不同对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
下面将参照附图,对根据本申请的具体实施例进行详细说明。
图2a示出根据本申请示例实施例的一种空间编码探测器的侧面图。图2b示出根据本申请示例实施例的一种空间编码探测器的横截面示意图。
如图2a和图2b所示的空间编码探测器包括闪烁晶体阵列10和与闪烁晶体阵列10耦合的光电转换器件阵列20,其中,闪烁晶体阵列10包括m行n列的晶体条11,光电转换器件阵列20包括c行d列的光电转换器件,且m、n、c、d均为大于等于1的整数,且m、n、c、d不同时为1。
根据本申请的一些实施例,闪烁晶体阵列10中的多根晶体条11可以按照不同的排列形式组合成闪烁晶体阵列10,本领域技术人员应当注意的是,在图2a-图4的示例实施例中,晶体条11的形状为长方体形状,闪烁晶体阵列10为晶体条11紧密排列形成的长方体形状,但在本申请中,晶体条11还可以为其它棱数不同的棱柱形状、圆柱形状或者其它不规则形状,闪烁晶体阵列10也可以通过这些不同形状的晶体条11组合排列形成各种形状,比如立方体、单排、不规则形状等。
本领域技术人员还应当注意的是,在上述实施例中,晶体条11的排列是按照矩阵的行列方向排列,还可以按照非规则形式排列,比如:第一种情形、同一种规格的长方体晶体条按照非规则方式排列,比如同一种外形规格的晶体条相互在行或者列方向上交错排列,同一根晶体条对应临近的另一行或列的两根晶体条;第二种情形、不同外形规格的长方体晶体条按照非规则方式排列,比如两种外形规格大小的晶体条大小相互在行或者列方向上交错排列,规格较大的晶体条可以对应临近的另一行或列的多根规格较小的晶体条;第三种情形、晶体条的形状规格不完全相同,其排列方式也不完全规则,比如其中一部分晶体条为棱柱形状,另一部分为圆柱形状,多种规格的晶体条相互交错排列。以上这些排列方式和形状,通过对应的空间编码,同样可以在投影平面上获得对应的空间分布信息,在此不再赘述。
根据本申请的一些实施例,在闪烁晶体阵列10中,至少有两根晶体条11具有不同的衰减时间,优选地,闪烁晶体阵列10中的每根晶体条11均具有不同的衰减时间。根据一些实施例,同一闪烁晶体阵列10中的多根晶体条11的衰减时间呈等差数列,比如闪烁晶体10中包括m×n根晶体条11,这些晶体条11的衰减时间均介于30ns-50ns之间且呈等差数列。利用不同晶体条具有不同的衰减时间,将晶体条按特定规律进行排列,从而将晶体条的空间信息编码到闪烁晶体阵列的衰减时间中。
在本申请的实施例中,由于高能射线与闪烁晶体相互作用产生可见光,可见光光子的产生速率在一段很短的时间内达到最高值,然后按指数模型下降为0。因此,衰减时间通常定义为可见光产生速率下降至最高点产生速率的1/e(e为自然底数)所需要的时间。衰减时间是闪烁晶体的固有性质,由晶体的种类和制造工艺所决定,可在一定范围内按需求生产特定衰减时间的闪烁晶体,在此不再赘述。
本领域技术人员需要注意的是,对于不同的高能射线,衰减时间的定义可以有所不同,本领域技术人员可以根据需要定义不同的衰减时间,比如将衰减时间定义为可见光产生速率下降至最高点产生速率的2/e所需要的时间,只需要使得闪烁晶体阵列10中至少一部分晶体条具有不同的衰减时间,从而让光子入射到光电转换器件阵列所需的时间不同即可,这是本领域技术人员根据本申请的启示所容易想到的,在此不再赘述。
在计算射线的沉积位置时,与闪烁晶体阵列耦合的光电转换器件阵列的空间分辨率限制了分辨高能射线沉积所对应区域的能力,而通过在闪烁晶体阵列中设置多根具有不同衰减时间的晶体条,让不同的区域中对应的不同衰减时间信息与晶体条编码对应,并通过衰减时间解码出额外空间信息,进一步能够定位闪烁晶体阵列中具体发生沉积的晶体条,从而能够突破光电转换器件阵列分辨率对探测器空间分辨率的限制。
图2a和图2b所示的实施例可以用于包括闪烁晶体阵列的探测器中,例如,包括单层闪烁晶体阵列的探测器或包括多层闪烁晶体阵列的探测器。通过在闪烁晶体阵列中使用具有不同衰减时间的晶体条并解码射线的衰减时间,既能够突破单层阵列空间编码探测器中光电转换器件阵列的分辨率对探测器空间分辨率的限制,又能够解决多层或者组合闪烁晶体探测器中误码率高的问题。
进一步地,根据本申请的一些实施例,可以将闪烁晶体阵列10中的晶体条11划分为多个晶体单元,每个晶体单元包括多根晶体条11,且每个晶体单元中的至少一部分晶体条具有不同的衰减时间,或者不同的晶体单元之间具有不同的衰减时间。根据一些实施例,同一晶体单元中的多根晶体条的衰减时间可以呈等差数列。利用闪烁晶体的不同衰减时间,将晶体条按特定规律进行排列,从而将晶体条的空间信息编码到闪烁晶体阵列的衰减时间中。
根据一些实施例,多个晶体单元的大小可以相同,也可以不同,不同晶体单元中晶体条的数量、外形尺寸规格、衰减时间也可以不完全相同。
根据一些实施例,当闪烁晶体阵列为长方体矩阵形式时,可以根据闪烁晶体阵列的行或列的排列顺序划分晶体单元,每个晶体单元的大小可以相同,也可以不同。
根据一些实施例,可以根据闪烁晶体阵列的外圈和内圈的排列顺序划分晶体单元,不同晶体单元中每根晶体条的衰减时间可以完全相同、不完全相同或完全不同,不同晶体单元之间至少有一部分晶体条的衰减时间不同。
图3示出根据本申请示例实施例的一种闪烁晶体阵列分块示意图。根据图3的实施例,将闪烁晶体阵列按照晶体条的行列顺序划分为多个大小相同的晶体单元1、2、3、……、s,其中,s为自然数,在图3的实施例中,闪烁晶体阵列10被分为s个大小相同的晶体单元,每一个晶体单元包含大小相同的4根晶体条,具体晶体条的编号可如图2b中所示。
根据一些实施例,同一晶体单元中的晶体条的衰减时间相同,不同晶体单元中的晶体条的衰减时间不同。
例如,晶体单元1中的晶体条1-1、2-1、1-2、2-2的衰减时间相同,均为30ns,晶体单元2中的晶体条3-1、4-1、3-2、4-2的衰减时间相同,均为30.5ns,依次类推。
根据一些实施例,同一晶体单元中的至少一部分晶体条的衰减时间不同于其他晶体条,不同晶体单元中的晶体条的衰减时间不完全相同。
例如,晶体单元1中的晶体条1-1、2-2的衰减时间相同,均为30ns,晶体条2-1、1-2的衰减时间相同,均为31ns;晶体单元2中的晶体条3-1、4-2的衰减时间相同,均为30.5ns,晶体条4-1、3-2、衰减时间相同,均为31ns,依次类推。
根据一些实施例,同一晶体单元中的晶体条具有不同的衰减时间。
例如,晶体单元1中的晶体条1-1、2-1、1-2和2-2的衰减时间互不相同。晶体单元2、晶体单元3、晶体单元4、……晶体单元s中各晶体条的衰减时间设置方式和晶体单元1中的设置方式相同。
根据一些实施例,同一晶体单元中的晶体条的衰减时间呈等差数列。
根据一些实施例,同一晶体单元中的晶体条的衰减时间呈公差范围在30~50ns内的等差数列。
根据本申请的一些示例实施例,不同晶体单元中的同一位置的晶体条具有相同的衰减时间,同一晶体单元中不同位置的晶体条的衰减时间不完全相同。
例如,晶体条1-1的衰减时间和晶体条3-1、晶体条1-3、晶体条3-3的衰减时间相同,晶体条2-1的衰减时间和晶体条4-1、晶体条2-3、晶体条4-3的衰减时间相同,晶体条1-2的衰减时间和晶体条3-2、晶体条1-4、晶体条3-4的衰减时间相同,晶体条2-2的衰减时间和晶体条4-2、晶体条2-4、晶体条4-4的衰减时间相同,但晶体条1-1、1-2、2-1、2-2的衰减时间不完全相同。
根据图3所示的实施例,通过将闪烁晶体阵列划分为多个大小相同的晶体单元,且至少一部分晶体单元具有不同的衰减时间,以便在计算射线的沉积位置时,能够通过衰减时间解码出的额外空间信息进一步定位晶体单元中具体发生沉积的晶体条,从而能够突破光电转换器件阵列的空间分辨率对探测器空间分辨率的限制。
图4示出根据本申请示例实施例的另一种闪烁晶体阵列分块示意图。根据图4的实施例,将闪烁晶体阵列按照晶体条的行列顺序划分为多个大小不完全相同的晶体块,如图4所示,颜色较粗的线条表示同一个晶体单元,同一个晶体单元内的晶体条用较细的线条表示,图4中闪烁晶体阵列包括大小不完全相同的多个晶体单元30、40、50,其中,晶体单元30中包括4×4根晶体条11,晶体单元40中包括4×2根晶体条11,晶体单元50中包括2×2根晶体条11。
在图4的实施例中,同一晶体单元中晶体条的衰减时间设置、不同晶体单元中晶体条的衰减时间设置均可以与图3实施例中相同,在此不再赘述。
根据图4所示的实施例,通过将闪烁晶体阵列划分为多个大小不同的晶体单元,且至少一部分晶体单元具有不同的衰减时间,以便在计算射线的沉积位置时,能够通过衰减时间解码出的额外空间信息进一步定位不同晶体单元中具体发生沉积的晶体条,从而能够突破光电转换器件阵列空间分辨率对探测器空间分辨率的限制。
图5示出根据本申请示例实施例的一种空间编码方法流程图。图5所示的空间编码方法用于包括如图2a-图4所示的闪烁晶体阵列探测器。
在步骤S501,建立闪烁晶体阵列的空间位置编码表,空间位置编码表中包含了晶体条的衰减时间、尺寸与探测器的空间分辨率之间的对应关系。
根据本申请的一些实施例,空间位置编码表中还可以包括不同位置、不同尺寸、不同衰减时间的晶体条按照预设编码顺序、编码位置排列后与探测器的空间分辨率之间的对应关系。
根据本申请的一些实施例,晶体条的预设编码顺序、编码位置以及尺寸、衰减时间的设置可以与图2a-图4中的实施例相同,在此不再赘述。
在步骤S503,根据实际应用中不同的空间分辨率需求从步骤S501中的空间编码表中选取对应规格的晶体条,并按照预设编码顺序或编码位置排列。
根据一些实施例,同样的空间分辨率可以对应多种形式的编码顺序或者编码位置。
在步骤S505,根据预设的编码顺序或编码位置完成晶体条的组装,获得对应空间分辨率的空间编码探测器。
根据一些实施例,组装后闪烁晶体阵列与光电转换器以一对一的方式耦合。当射线在晶体条内沉积后产生的可见光光子只进入对应的一个光电转换器件中。
根据一些实施例,组装后闪烁晶体阵列与光电转换器以一对多的方式耦合。当射线在晶体条内沉积后产生的可见光光子可以进入对应的多个光电转换器件中。
根据一些实施例,组装后闪烁晶体阵列与光电转换器以多对一的方式耦合。多根晶体条内产生的可见光光子会进入同一个光电转换器件中。根据本申请的一些示例实施例,利用重心法确定以多对一的方式耦合的空间编码探测器中射线沉积的位置。例如,首先计算各个光电转换器件的中心坐标和所探测到光子数的加权平均坐标,通过比较加权平均坐标和中心坐标确定射线沉积的具体位置。
根据图5所示的实施例,通过空间位置编码表设计探测器,使至少一部分晶体条具有不同的衰减时间,以便在探测射线时,能够通过衰减时间解码出的额外空间信息进一步精确定位具体发生沉积的晶体条,从而能够突破光电转换器件阵列空间分辨率对探测器空间分辨率的限制。
图6示出根据本申请示例实施例的一种空间解码方法流程图。图5所示的空间解码方法用于包括如图2a-图4所示的空间编码探测器。
如图6所示,在步骤S601,利用光电转换器或光电转换器阵列输出的脉冲信号获取射线的信息,其中,获取的射线的信息包括能量、到达时间和衰减时间。
根据本申请的一些实施例,利用多电压阈值(Multi VoltageThreshold,简称MVT)采样方法、示波器或模数转换器对光电转换器输出的脉冲信号进行采样,以获取射线的信息,例如,射线的能量信息、到达时间和衰减时间。
在步骤S603,利用步骤S601获得的射线的能量信息和到达时间确定射线沉积的晶体单元。
根据一些实施例,空间编码探测器中的闪烁晶体阵列与光电转换器以一对一的方式耦合。射线在闪烁晶体阵列沉积后产生的可见光光子只进入一个光电转换器件。根据本申请的一些示例实施例,利用一对一耦合法确定以一对一的方式耦合的空间编码探测器中射线沉积的块。
根据一些实施例,组装后闪烁晶体阵列与光电转换器以一对多的方式耦合。当射线在晶体条内沉积后产生的可见光光子可以进入对应的多个光电转换器件中。
根据一些实施例,空间编码探测器中的闪烁晶体阵列与光电转换器以多对一的方式耦合。多根晶体条内产生的可见光光子进入同一个光电转换器件。根据本申请的一些示例实施例,利用重心法确定以多对一的方式耦合的空间编码探测器中射线沉积的块。例如,首先计算各个光电转换器件的中心坐标和所探测到光子数的加权平均坐标,通过比较加权平均坐标和中心坐标确定射线沉积的块。
在步骤S605,利用射线的衰减时间和预先建立的查找表获取射线的沉积位置。其中,查找表包括衰减时间与沉积位置的对应关系。
根据本申请的一些实施例,在步骤S605之前需要预先建立查找表。查找表中包括射线的衰减时间与沉积位置的对应关系,查找表的建立方式如图7所示。通过步骤S603确定射线沉积的晶体单元后,利用步骤S601获得的衰减时间和预先建立的查找表可以确定射线沉积的具体晶体条。
根据图6所示的实施例,首先确定射线沉积的晶体单元,然后利用各根晶体条的编码信息,进一步定位确定的沉积晶体单元中具体发生沉积的晶体条,突破了光电转换器件的分辨率对探测器空间分辨率的限制。
图7示出根据本申请示例实施例的一种建立查找表的方法流程图。
如图7所示,在步骤S701,利用不同规格的晶体条探测射线。
根据一些实施例,不同规格包括晶体条的尺寸不同、衰减时间不同或者晶体条外侧涂覆的物质不同。射线沉积在不同规格的晶体条中,探测器采集到的脉冲的参数信息有所不同,参数信息包括脉冲到达时间、能量等信息。
在步骤S703,计算每根晶体条中沉积的射线的衰减时间。
根据一些实施例,利用多电压阈值采样方法、示波器或模数转换器对光电转换器输出的脉冲信号进行采样,以计算每根晶体条中沉积的射线的衰减时间。
在步骤S705,利用所计算的衰减时间和射线在晶体条中的沉积位置建立查找表。
根据图7所示的实施例建立的查找表包括射线的衰减时间与射线沉积的晶体条位置的对应关系。
图8示出根据本申请示例实施例的一种空间解码装置框图。根据一些实施例,图8所示的装置用于如图2a-图4所示的空间编码探测器。
如图8所示的一种空间解码装置包括射线信息获取单元801、晶体区域确定单元803和沉积位置确定单元805,其中,射线信息获取单元801用于利用光电转换器输出的脉冲信号获取射线的信息,其中,获取的射线的信息包括能量信息、到达时间和衰减时间;晶体区域确定单元803用于利用能量信息和到达时间确定射线沉积的晶体单元;沉积位置确定单元805用于利用衰减时间和预先建立的查找表获取射线的沉积位置。查找表包括射线的衰减时间与沉积位置的对应关系。
图9示出了根据本申请实施例的又一种空间解码装置。图9示出的空间解码装置仅仅是一个示例,不应对本申请实施例的功能和使用范围带来任何限制。
如图9所示,该空间解码装置以通用计算设备的形式表现。该空间解码装置的组件可以包括但不限于:至少一个处理器910、至少一个存储器920、连接不同系统组件(包括存储器920和处理器910)的总线930、显示单元940等。其中,存储器920存储有程序代码,程序代码可以被处理器910执行,使得处理器910执行本说明书描述的根据本申请各种示例性实施方式的方法。例如,处理器910可以执行如图6中所示的方法。
存储器920可以包括易失性存储单元形式的可读介质,例如随机存取存储单元(RAM)9201和/或高速缓存存储单元9202,还可以进一步包括只读存储单元(ROM)9203。
存储器920还可以包括具有一组(至少一个)程序模块9205的程序/实用工具9204,这样的程序模块9205包括但不限于:操作系统、一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。
总线930可以为表示几类总线结构中的一种或多种,包括存储单元总线或者存储单元控制器、外围总线、图形加速端口、处理单元或者使用多种总线结构中的任意总线结构的局域总线。
空间解码装置也可以与一个或多个外部设备900(例如键盘、指向设备、蓝牙设备等)通信,还可与一个或者多个使得用户能与该空间解码装置交互的设备通信,和/或与使得该空间解码装置能与一个或多个其它计算设备进行通信的任何设备(例如路由器、调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口950进行。而且,空间解码装置还可以通过网络适配器960与一个或者多个网络(例如局域网(LAN),广域网(WAN)和/或公共网络,例如因特网)通信。网络适配器960可以通过总线930与空间解码装置的其它模块通信。应当明白,尽管图中未示出,可以结合空间解码装置使用其它硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统等。
通过以上的实施方式的描述,本领域的技术人员易于理解,这里描述的示例实施方式可以通过软件实现,也可以通过软件结合必要的硬件的方式来实现。根据本申请实施方式的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个计算机可读存储介质(可以是CD-ROM、U盘、移动硬盘等)中或网络上,包括若干计算机程序指令以使得一台计算设备(可以是个人计算机、服务器、或者网络设备等)执行根据本申请实施方式的上述方法。
软件产品可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以为但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
计算机可读存储介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了可读程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。可读存储介质还可以是可读存储介质以外的任何可读介质,该可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。可读存储介质上包含的程序代码可以用任何适当的介质传输,包括但不限于无线、有线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言的任意组合来编写用于执行本申请操作的程序代码,程序设计语言包括面向对象的程序设计语言—诸如Java、C++等,还包括常规的过程式程序设计语言—诸如C语言或类似的程序设计语言。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户计算设备上部分在远程计算设备上执行、或者完全在远程计算设备或服务器上执行。在涉及远程计算设备的情形中,远程计算设备可以通过任意种类的网络,包括局域网(LAN)或广域网(WAN),连接到用户计算设备,或者,可以连接到外部计算设备(例如利用因特网服务提供商来通过因特网连接)。
上述计算机可读介质承载有一个或者多个程序指令,当上述一个或者多个程序指令被一个该设备执行时,使得该计算机可读介质实现前述功能。
本领域技术人员可以理解上述各模块可以按照实施例的描述分布于装置中,也可以进行相应变化唯一不同于本实施例的一个或多个装置中。上述实施例的多个模块可以合并为一个模块,也可以进一步将一个模块拆分成多个子模块。
通过以上的实施方式的描述,本领域的技术人员易于理解,这里描述的示例实施方式可以通过软件实现,也可以通过软件结合必要的硬件的方式来实现。根据本申请实施方式的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个计算机可读存储介质(可以是CD-ROM、U盘、移动硬盘等)中或网络上,包括若干计算机程序指令以使得一台计算设备(可以是个人计算机、服务器、或者网络设备等)执行根据本申请实施方式的上述方法。
软件产品可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以为但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
计算机可读存储介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了可读程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。可读存储介质还可以是可读存储介质以外的任何可读介质,该可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。可读存储介质上包含的程序代码可以用任何适当的介质传输,包括但不限于无线、有线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言的任意组合来编写用于执行本申请操作的程序代码,程序设计语言包括面向对象的程序设计语言—诸如Java、C++等,还包括常规的过程式程序设计语言—诸如C语言或类似的程序设计语言。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户计算设备上部分在远程计算设备上执行、或者完全在远程计算设备或服务器上执行。在涉及远程计算设备的情形中,远程计算设备可以通过任意种类的网络,包括局域网(LAN)或广域网(WAN),连接到用户计算设备,或者,可以连接到外部计算设备(例如利用因特网服务提供商来通过因特网连接)。
上述计算机可读介质承载有一个或者多个程序指令,当上述一个或者多个程序指令被一个该设备执行时,使得该计算机可读介质实现前述功能。
本领域技术人员可以理解上述各模块可以按照实施例的描述分布于装置中,也可以进行相应变化唯一不同于本实施例的一个或多个装置中。上述实施例的多个模块可以合并为一个模块,也可以进一步将一个模块拆分成多个子模块。
在传统阵列闪烁晶体探测器中,探测器获取沉积投影位置信息的能力受限于光电转换器件阵列空间分辨率。根据本申请的一些实施例,将闪烁晶体阵列中各个闪烁晶体条的衰减时间进行空间编码,可以通过衰减时间解码出的额外空间信息进一步定位具体发生沉积的晶体条,从而能够突破光电转换器件阵列空间分辨率对探测器空间分辨率的限制。
虽然本申请提供了如上述实施例或流程图所述的方法操作步骤,但基于常规或者无需创造性的劳动在所述方法中可以包括更多或者更少的操作步骤。在逻辑性上不存在必要因果关系的步骤中,这些步骤的执行顺序不限于本申请实施例提供的执行顺序。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其它实施例的不同之处。
以上对本申请实施例进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明仅用于帮助理解本申请的方法及其核心思想。同时,本领域技术人员依据本申请的思想,基于本申请的具体实施方式及应用范围上做出的改变或变形之处,都属于本申请保护的范围。综上所述,本说明书内容不应理解为对本申请的限制。

Claims (23)

1.一种空间编码探测器,其特征在于,包括:
闪烁晶体阵列,包括至少两根晶体条,其中至少两根所述晶体条具有不同的衰减时间;以及
光电转换器,所述光电转换器与所述晶体阵列耦合。
2.根据权利要求1所述的空间编码探测器,其特征在于,所述闪烁晶体阵列中包括至少一个晶体单元,所述晶体单元中包括多根所述晶体条。
3.根据权利要求2所述的空间编码探测器,其特征在于,所述晶体单元中至少两根所述晶体条具有不同的衰减时间,或者所述晶体单元中至少一部分所述晶体条的衰减时间与至少一根所述晶体条的衰减时间不同。
4.根据权利要求2所述的空间编码探测器,其特征在于,所述至少一个晶体单元的大小相同,不同所述晶体单元之间至少一部分所述晶体条的衰减时间不同。
5.根据权利要求2所述的空间编码探测器,其特征在于,所述闪烁晶体阵列包括至少两种大小的晶体单元,不同所述晶体单元之间至少一部分所述晶体条的衰减时间不同。
6.根据权利要求4或5所述的空间编码探测器,其特征在于,所述晶体单元中的所述晶体条呈矩阵排列,不同所述晶体单元之间对应位置的所述晶体条具有相同的衰减时间。
7.根据权利要求1或2所述的空间编码探测器,其特征在于,每根所述晶体条的衰减时间均不同。
8.根据权利要求7所述的空间编码探测器,其特征在于,所述晶体条的衰减时间呈等差数列。
9.根据权利要求7所述的空间编码探测器,其特征在于,所述晶体条的衰减时间介于30ns~50ns之间。
10.根据权利要求1所述的空间编码探测器,其特征在于,其中一部分所述晶体条的外形尺寸不同。
11.根据权利要求1所述的空间编码探测器,其特征在于,所述晶体条的外形包括棱柱、圆柱或者不规则形状。
12.根据权利要求1所述的空间编码探测器,其特征在于,所述光电转换器为光电转换器件阵列,所述光电转换器件阵列中的每个单元与所述晶体条之间为一对一、一对多或者多对一的形式耦合。
13.一种空间编码方法,其特征在于,所述空间编码方法包括:
建立闪烁晶体阵列的空间位置编码表,所述空间位置编码表包括不同衰减时间、不同尺寸的晶体条与空间分辨率之间的对应关系;
根据空间分辨率需求从所述空间位置编码表中选取对应规格的所述晶体条组成所述闪烁晶体阵列;
根据预设编码顺序或者预设编码位置耦合所述闪烁晶体阵列与光电转换器。
14.根据权利要求13所述的空间编码方法,其特征在于,所述空间位置编码表还包括不同位置、不同尺寸或者不同衰减时间的晶体条按照所述预设编码顺序或者所述预设编码位置排列后与空间分辨率之间的对应关系。
15.根据权利要求13所述的空间编码方法,同一所述空间分辨率对应多种形式的所述预设编码顺序或者所述预设编码位置。
16.根据权利要求13所述的空间编码方法,其特征在于,所述光电转换器为光电转换器件阵列,所述光电转换器件阵列中的每个单元与所述晶体条之间为一对一、一对多或者多对一的形式耦合。
17.一种空间解码方法,其特征在于,所述空间解码方法用于如权利要求1-12中任一项所述的空间编码探测器,所述空间解码方法包括:
利用光电转换器输出的脉冲信号获取所述射线的信息,所述信息包括能量信息、到达时间和衰减时间;
利用所述能量信息和所述到达时间确定所述射线沉积的晶体单元;
利用所述衰减时间和预先建立的查找表获取所述射线的沉积位置,所述查找表包括所述衰减时间与射线沉积位置的对应关系。
18.根据权利要求17所述的方法,其特征在于,所述利用光电转换器输出的脉冲信号获取所述射线的信息,包括:
利用多电压阈值采样方法、示波器或模数转换器获取所述射线的信息。
19.根据权利要求17所述的方法,其特征在于,所述确定射线沉积的晶体单元,包括:
根据所述闪烁晶体阵列与所述光电转换器的耦合形式利用一对一耦合法或者重心法确定所述射线沉积的晶体单元。
20.根据权利要求17所述的方法,其特征在于,所述查找表通过如下方式预先建立:
利用不同规格的所述晶体条探测所述射线;
计算每根所述晶体条中沉积的所述射线的衰减时间;
利用所计算的衰减时间和所述射线在不同晶体条中的沉积位置建立所述查找表。
21.一种空间解码装置,其特征在于,所述空间解码装置用于如权利要求1-12中任一项所述的空间编码探测器,所述空间解码装置包括:
射线信息获取单元,用于利用光电转换器输出的脉冲信号获取射线的信息,所述信息包括能量信息、到达时间和衰减时间;
晶体区域确定单元,用于利用所述能量信息和所述到达时间确定所述射线沉积的晶体单元;
沉积位置确定单元,用于利用所述衰减时间和预先建立的查找表获取所述射线的沉积位置,所述查找表包括所述衰减时间与所述沉积位置的对应关系。
22.一种空间解码装置,其特征在于,包括:
一个或多个处理器;
存储装置,用于存储计算机程序;
当所述计算机程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现如权利要求17-20中任一项所述的方法。
23.一种存储介质,其特征在于,其上存储有程序指令,所述程序指令被执行时实现如权利要求17-20中任一项所述的方法。
CN202210372854.7A 2022-04-11 2022-04-11 空间编码探测器、空间编码、解码方法、装置及存储介质 Pending CN114895344A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210372854.7A CN114895344A (zh) 2022-04-11 2022-04-11 空间编码探测器、空间编码、解码方法、装置及存储介质
PCT/CN2022/136854 WO2023197626A1 (zh) 2022-04-11 2022-12-06 空间编码探测器、空间编码、解码方法、装置及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210372854.7A CN114895344A (zh) 2022-04-11 2022-04-11 空间编码探测器、空间编码、解码方法、装置及存储介质

Publications (1)

Publication Number Publication Date
CN114895344A true CN114895344A (zh) 2022-08-12

Family

ID=82715304

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210372854.7A Pending CN114895344A (zh) 2022-04-11 2022-04-11 空间编码探测器、空间编码、解码方法、装置及存储介质

Country Status (2)

Country Link
CN (1) CN114895344A (zh)
WO (1) WO2023197626A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023197626A1 (zh) * 2022-04-11 2023-10-19 苏州瑞派宁科技有限公司 空间编码探测器、空间编码、解码方法、装置及存储介质

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06337289A (ja) * 1993-05-27 1994-12-06 Shimadzu Corp 放射線検出器及びそれを用いた装置
JP2000056023A (ja) * 1998-08-06 2000-02-25 Japan Science & Technology Corp 深さ方向の情報を得ることができるpet用検出器
CN101460864B (zh) * 2006-09-19 2012-02-01 株式会社岛津制作所 核医学诊断装置的闪烁体阵列的识别机构所需的参数的决定方法
CN101806912A (zh) * 2010-03-18 2010-08-18 清华大学 高能射线叠层式晶体模块探测器
CN107456235B (zh) * 2017-07-06 2021-03-05 东软医疗系统股份有限公司 位置校正方法和系统
CN114152971A (zh) * 2021-12-03 2022-03-08 中国电子科技集团公司第二十六研究所 闪烁晶体发光衰减时间的测量方法及测量系统
CN114895344A (zh) * 2022-04-11 2022-08-12 苏州瑞派宁科技有限公司 空间编码探测器、空间编码、解码方法、装置及存储介质

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023197626A1 (zh) * 2022-04-11 2023-10-19 苏州瑞派宁科技有限公司 空间编码探测器、空间编码、解码方法、装置及存储介质

Also Published As

Publication number Publication date
WO2023197626A1 (zh) 2023-10-19

Similar Documents

Publication Publication Date Title
JP4338177B2 (ja) 3次元放射線位置検出器
US3978337A (en) Three-dimensional time-of-flight gamma camera system
US10365385B2 (en) Apparatus, method and system for sparse detector
US8384034B2 (en) Beta ray detector and beta ray reconstruction method
US20110127435A1 (en) Gamma ray detector, radiation diagnostic device, tomography device, and method of analyzing tomography device
Martin et al. Imaging multi-energy gamma-ray fields with a Compton scatter camera
CN114895344A (zh) 空间编码探测器、空间编码、解码方法、装置及存储介质
Sinclair et al. Simulations of a scintillator Compton gamma imager for safety and security
CN106859686B (zh) 成像方法和成像系统
Finocchiaro Radioactive waste: A system for online monitoring and data availability
Nakamura et al. A half-millimetre spatial resolution fibre-coded linear position-sensitive scintillator detector with wavelength-shifting fibre read-out for neutron detection
CN114910946A (zh) 空间编码晶体阵列、探测器、方法、装置以及存储介质
EP0846275A1 (en) An imaging apparatus
Tzanakos et al. Development and validation of a simulation model for the design of a PET scanner
Solevi et al. A Monte-Carlo based model of the AX-PET demonstrator and its experimental validation
Smith et al. Design and modeling of the hybrid portable gamma camera system
Tuemer et al. A new telescope for wide-band gamma-ray astronomy: The Silicon Compton Recoil Telescope (SCRT)
Belcari et al. Positron emission tomography: its 65 years and beyond
Anderson et al. Constructing a small laboratory animal imaging device based on scintillating fibers
Ryan et al. SONTRAC: An imaging spectrometer for solar neutrons
CN112130190A (zh) 一种基于塑料闪烁体的伽马探测器及其正电子成像方法
Zong Study of cosmic rays by Auger and LHAASO: R&D and Data Analysis of AugerPrime and simulations for LHAASO
Clements Simulations of a high-resolution micro-PET system based on liquid xenon.
Eisch A Multi-Variate Fit to the Chemical Composition of the Cosmic-Ray Spectrum
Ito et al. Design of a 2-D coded aperture imaging system using 32 gamma ray 1-D array detectors

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