CN111638544A - 基于缝孔混合准直器的多伽马光子符合成像系统及方法 - Google Patents

基于缝孔混合准直器的多伽马光子符合成像系统及方法 Download PDF

Info

Publication number
CN111638544A
CN111638544A CN202010465659.XA CN202010465659A CN111638544A CN 111638544 A CN111638544 A CN 111638544A CN 202010465659 A CN202010465659 A CN 202010465659A CN 111638544 A CN111638544 A CN 111638544A
Authority
CN
China
Prior art keywords
gamma photon
gamma
coincidence
slit
collimator
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
CN202010465659.XA
Other languages
English (en)
Other versions
CN111638544B (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN202010465659.XA priority Critical patent/CN111638544B/zh
Publication of CN111638544A publication Critical patent/CN111638544A/zh
Application granted granted Critical
Publication of CN111638544B publication Critical patent/CN111638544B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914Measurement of spatial distribution of radiation
    • G01T1/2985In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
    • 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/161Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164Scintigraphy
    • G01T1/1641Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1642Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras using a scintillation crystal and position sensing photodetector arrays, e.g. ANGER cameras
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

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)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Measurement Of Radiation (AREA)
  • Nuclear Medicine (AREA)

Abstract

本发明公开了一种基于缝孔混合准直器的多伽马光子符合成像系统及方法,所述系统包括时间符合模块、计算机平台、缝孔混合准直器和伽马光子探测器;缝孔混合准直器是通过在合金板上开设若干准直缝和准直孔所形成,合金板采用平面板或曲面板;所述成像方法将放射性核素发生衰变的位置范围缩小为多伽马光子符合事件中通过准直缝在伽马光子探测器探测到的伽马光子事件所确定的若干投影平面和通过准直孔在伽马光子探测器探测到的伽马光子事件所确定的若干投影线的在成像范围内的若干交点,以获得放射性核素在被测范围内分布。本发明提高了多伽马光子符合事件的探测效率,提高了重建图像的信噪比从而降低了对伽马光子总计数的需求。

Description

基于缝孔混合准直器的多伽马光子符合成像系统及方法
技术领域
本发明属于发射断层成像技术领域,特别涉及一种基于缝孔混合准直器的多伽马光子符合成像系统及方法。
背景技术
发射断层成像技术是当今探测物体内部信息的重要技术之一,在多种研究领域中有着大量的应用。在被探测物体外通过探测从被探测物体发出的伽马光子,并通过图像重建从而以非侵入式的手段观察得到被探测物体的内部信息。在发射断层成像领域,最为重要的几种成像系统分别有正电子发射断层成像仪(Positron Emission Tomography,以下简称PET)、单光子发射计算机断层成像仪(Single-Photon Emission ComputedTomography,以下简称SPECT)和康普顿相机(Compton Camera)等。目前PET和SPECT均已广泛用于临床检查和诊断,包括癌症、神经系统疾病和心血管疾病等,而康普顿相机也应用于核能工业、天文、医疗等方面。
PET的核心组件由很多个具备时间测量功能的伽马光子探测器模块以及相应的时间符合模块构成。PET的基本原理涉及电子准直技术,所利用的放射性核素为正电子核素。正电子核素放射出的正电子会在被探测物体内发生正负电子湮灭效应,产生一对能量为511keV、方向几乎相反的伽马光子。采用时间符合测量手段,即如果两个伽马光子探测器模块在一个很短的时间内(通常在几百个纳秒)分别探测到这两个能量为511keV的伽马光子,就能确定一条正电子湮灭的发生位置(近似为正电子核素发生衰变的位置)所在的响应线(Line of Response)。记录大量这样的响应线通过图像重建即可获得正电子湮灭发生位置的分布,即近似为正电子核素在被探测物体内的分布。由于正电子湮灭产生的一对伽马光子的发射方向是几乎相反的,因此只能确定正电子湮灭的发生位置在相应的响应线上,而不能确定正电子湮灭发生在相应响应线上的具体位置。尽管可以通过飞行时间测量技术初步确定正电子湮灭发生在相应响应线上的位置范围,但这需要伽马光子探测器模块具备极高的时间分辨率。由于正电子湮灭在响应线上发生位置的不确定性,重建的正电子核素在被探测物体内分布图像的信噪比往往较低,影响诊断的效果。为提高图像信噪比通常需要累积大量的响应线,而这会使得被探测物体摄入较大剂量的正电子核素,增加被探测物体的辐照风险。
SPECT的核心组件包括准直器及伽马光子探测器模块等。SPECT利用的是物理准直技术,其所利用的核素为伽马光子核素。在伽马光子探测器模块的前端通常会放置准直器来限制伽马光子核素放射出的伽马光子到达探测器的入射角度,使得仅沿特定方向发射的伽马光子才能通过准直器被探测器探测到,伽马光子探测器每探测到一个伽马光子就能确定一条伽马光子初始发射位置所在的投影线。积累大量这样的投影线通过图像重建就能确定伽马光子初始发射位置的分布,即伽马光子核素在被探测物体内的分布。与PET存在类似的缺陷,SPECT也不能确定伽马光子在投影线上的具体发射位置,因此重建图像的信噪比较差。另外由于SPECT使用了准直器,限制了能被探测器探测到的伽马光子的发射角度,成像系统的探测效率较低,而这进一步使得重建图像的信噪比恶化。
康普顿相机模块的核心组件包括两块平行的探测器平面模块等。采用时间符合测量手段,即如果两个平行的探测器平面模块在一个很短的时间内(通常在几个纳秒)分别探测到信号,则可以认为该两次信号来自于同一次伽马光子事件。伽马光子在第一块晶体平面上发生康普顿散射产生康普顿散射光子,所产生的康普顿散射光子在第二块晶体平面发生光电效应被吸收,康普顿相机从而得以探测到该次伽马光子事件。根据伽马光子在第一块晶体平面上沉积能量与第二块晶体平面上沉积能量之和,可以得到伽马光子的总能量。康普顿相机每探测到一个伽马光子事件,根据伽马光子事件在第一块晶体平面上的沉积能量、沉积位置和第二块晶体平面上的沉积能量、沉积位置,就能根据康普顿散射原理确定一个伽马光子初始发射位置所在的投影圆锥面。
本申请人已提出的一种多伽马光子同时发射药物时间符合核医学成像系统及方法(申请号:201610798146.4),该系统包括多个以非平行方式排列的探测器探头、时间符合模块以及计算机平台构成,每个探测器探头均由准直器和具备时间测量功能的伽马光子探测器构成,探测放射性核素在很短时间内放射的多个伽马光子构成多伽马光子符合事件;所述方法计算到多伽马光子符合事件中的每一个伽马光子事件所确定若干投影线的距离之和最短的点的位置即为放射性核素发生衰变的位置,积累一定数量的多伽马光子符合事件即可实现放射性核素在生物体内分布的获取。该成像系统和方法简化了重建算法,提高了重建图像的信噪比;降低了对伽马光子总计数的需求,降低了病人的辐照风险。但是,由于该系统探测器探头包含限制伽马光子入射的平行孔准直器和针孔准直器,其单探测器探头的探测效率非常低(如所采用的平行孔准直器只允许垂直于其准直器平面的单个方向的伽马光子通过),进而导致多光子符合事件探测效率极低。除此之外,由于投影线之间不一定完美相交于一点,这不仅加大了判定交点位置的难度,也进一步导致了符合探测效率较低。低符合探测效率导致使得探测数据的信噪比低,对图像的重建质量有一定影响。
为了解决上述系统符合探测效率低的问题,本申请人提出了一种多伽马光子符合成像系统及方法(申请号:201810230414.1),该系统包括时间符合模块、计算机平台、至少一个由准直器和伽马光子探测器构成的第一探头,至少一个由前后两个康普顿相机探测器构成的第二探头,各探头探测放射性核素放射的多个伽马光子构成多伽马光子符合事件;该成像方法将放射性核素发生衰变的位置范围缩小为多伽马光子符合事件中第一探头探测到的伽马光子事件所确定的投影线和第二探头探测到的伽马光子事件所确定的投影圆锥面的若干个在成像范围内的交点,积累一定数量的多伽马光子符合事件即可获得放射性核素在被测范围内分布的图像。通过使用至少一个康普顿相机探测器探头代替伽马探测器探头,由于康普顿相机探测器不存在准直器设计,相比于SPECT伽马探测器探测效率将会大大提高,从而提高了多光子符合事件的探测效率,提高了重建图像的信噪比。但是,此系统不能保证将伽马光子的发射位置确定为在成像范围内的一个交点,因此不能确保直接得到分布。除此之外,康普顿相机空间分辨率较差,对于康普顿效应较少发生的低能以及高能伽马光子探测效率较低,所成重建图像质量较差。
发明内容
本发明的目的旨在从原理上解决上述PET系统、SPECT系统和已提出的一种多伽马光子同时发射药物时间符合核医学成像系统及多伽马光子符合成像系统及方法的问题,本发明公开了一种基于缝孔混合准直器的多伽马光子符合成像系统。与现有的多伽马光子同时发射药物时间符合核医学成像技术的不同之处在于,本发明使用了缝孔混合准直器,可以允许任意方向的通过其准直缝和准直孔的伽马光子通过,一方面相较于单针孔准直器和平行孔准直器使得系统的单探测器探头的探测效率得到了提升,并且由于投影平面和投影线在非平行的条件下一定存在交点,不需要进行交点位置的判断,降低了符合事件判定的复杂度,进而提高了多光子符合事件的探测效率,提升了重建图像的信噪比;另一方面也避免了使用空间分辨率较差且对于低能以及高能伽马光子探测效率较低的康普顿相机,使得图像质量可以进一步提升;除此之外,可以通过适当的系统几何设计将直接得到每一个多伽马光子符合事件中的一个放射性核素衰变位置,为直接成像提供可能,并可以根据图像重建算法对结果进行进一步估计;最后,由于缝孔混合准直器设计的高度对称,使得更复杂形状的准直器设计成为可能。
为了实现上述目的,本发明采用如下技术方案:
本发明提出的一种基于缝孔混合准直器的多伽马光子符合成像系统,其特征在于,包括时间符合模块、计算机平台、缝孔混合准直器和具备时间测量功能的伽马光子探测器;所述缝孔混合准直器是通过在合金板上开设若干条准直缝和若干个准直孔所形成,所述缝孔混合准直器分别放置在伽马光子探测器前端与成像对象之间并与两者保持一定距离,使得以任意方向通过所述缝孔混合准直器的由成像对象内放射性核素发生衰变产生的伽马光子事件均能被该伽马光子探测器探测到;所述伽马光子探测器的时间信号线与所述时间符合模块连接,时间符合模块设定的时间窗用于判断穿过所述缝孔混合准直器的准直缝和准直孔被伽马光子探测器探测到的多个伽马光子事件是否构成多伽马光子符合事件,并将判断结果输入到所述计算机平台中;所述伽马光子探测器的能量和位置信号线还同时与计算机平台连接,用于多伽马光子符合事件的有效性判断以及多伽马光子符合事件发生位置可能所在的若干投影线和若干投影线在成像范围内的若干个交点的计算,从而获得放射性核素发生衰变时的可能位置;根据所有放射性核素的可能衰变位置,得到放射性核素在成像对象体内的可能分布,并可以通过图像重建算法对此分布做更为精确的估计;
所述缝孔混合准直器内的合金板的形状为平面和曲面中的任意一种或者两种。
本发明还提出采用上述成像系统的成像方法,其特征在于,该方法包括以下步骤:
(1)启动所述多伽马光子符合成像系统,设定时间符合模块的时间窗宽度;依据所用放射核素所放射的伽马光子能量设定伽马光子探测器的能窗;
(2)计算机平台根据所设定条件判断成像过程是否结束;若结束,则执行步骤(5);若未结束,则执行步骤(3);
(3)时间符合模块判断伽马光子探测器是否探测到多伽马光子符合事件,若探测到,则执行步骤(4);若未探测到,则执行步骤(2);
(4)计算机平台根据伽马光子探测器所输入的多个伽马光子能量信息判断所输入的多个伽马光子能量是否分别在所设定的能窗内;如果不在所设定的能窗内,则舍弃该多伽马光子符合事件;如在所设定的能窗内,则根据每一个多伽马光子符合事件在伽马光子探测器所输入的一个伽马光子的位置信息和缝孔混合准直器的准直缝位置信息,计算多伽马光子符合事件发生位置可能所在的若干投影平面,并根据每一个多伽马光子符合事件在伽马光子探测器所输入的一个伽马光子的位置信息和缝孔混合准直器的准直孔位置信息,计算多伽马光子符合事件发生位置可能所在的若干投影线;记录每一个多伽马光子符合事件中的多个伽马光子所分别确定的若干投影平面和若干投影线,计算每一个投影平面和每一条投影线之间在成像对象内的相应交点,作为放射性核素发生衰变的若干个可能位置;执行步骤(2);
(5)根据所有多伽马光子符合事件计算的放射性核素的多个可能衰变位置得到放射性核素在成像对象体内的可能分布,并可以通过图像重建算法对此分布做更为精确的估计。
本发明的特点及有益效果:
本发明提出的一种基于缝孔混合准直器的多伽马光子符合成像系统,有效地克服了传统PET系统或SPECT系统只能确定放射性核素衰变位置所在响应线或投影线而不能确定放射性核素衰变发生在响应线或投影线上具体位置的缺陷和已提出的一种多伽马光子同时发射药物时间符合核医学成像系统探测效率过低、信噪比过低的缺陷。本发明通过缝孔混合准直器和伽马光子探测器计算并获得多个伽马光子事件所确定的若干投影平面和若干投影线,计算每一个投影面和每一条投影线在成像范围内的若干个交点实现对放射性核素的衰变位置的定位,从而获得放射性核素在被探测物体内的可能分布。由于能够基于投影平面和投影线在成像范围内的若干个交点计算放射性核素的可能衰变位置,简化了图像重建算法,提高了重建图像的信噪比。由于采用了探测效率更高的缝孔混合准直器来替代单针孔准直器和平行孔准直器,从而提高了单探测器探头的探测效率,避免了投影交点的位置判断简化了算法,进一步提高了多光子符合事件的探测效率,降低了对伽马光子事件总计数的需求,减少了所需要的放射性核素剂量,与此同时还通过使用了对称性更高的设计,使得更复杂的准直器形状成为可能。
附图说明
图1是本发明实施例采用2个平板缝孔混合准直器和2个伽马光子探测器的成像系统的整体结构示意图;
图2、3分别是图1所示成像系统的探测及定位几何示意图;
图4、5分别是图1所示成像系统中平板缝孔混合准直器的两种结构示意图;
图6是本发明实施例采用1个平板缝孔混合准直器和1个伽马光子探测器的成像系统的整体结构示意图;
图7是图6所示成像系统的一种探测及定位几何示意图;
图8是本发明实施例采用4个平板缝孔混合准直器和4个伽马光子探测器的成像系统的整体结构示意图;
图9~图11分别是采用多个平板缝孔混合准直器和伽马光子探测器时的不同排布方式示意图;
图12是本发明实施例采用1个圆筒缝孔混合准直器和6个伽马光子探测器的成像系统的整体结构示意图;
图13是图12所示成像系统的一种三维排布示意图;
图14是本发明成像方法的流程框图。
图中:1—缝孔混合准直器、11—缝孔混合准直器的准直孔、12—缝孔混合准直器的准直缝、2—伽马光子探测器、3—时间符合模块、4—计算机平台、5—成像对象、6—衰变可能位置、7—投影平面、8—投影线。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步的详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不限定本发明的保护范围。
为了更好地理解本发明,以下详细阐述本发明提出的一种基于缝孔混合准直器的多伽马光子符合成像系统及方法的应用实例。
本发明提出的一种基于缝孔混合准直器的多伽马光子符合成像系统,包括时间符合模块、计算机平台、缝孔混合准直器和具备时间测量功能的伽马光子探测器;所述缝孔混合准直器是通过在各自合金板上分别开设若干条准直缝和若干个准直孔所形成,缝孔混合准直器放置在伽马光子探测器前端与成像对象之间并与两者保持一定距离,使得以任意方向通过缝孔混合准直器的准直缝和准直孔的由成像对象内放射性核素发生衰变产生的伽马光子事件均能被伽马光子探测器探测到;所有伽马光子探测器的时间信号线均与所述时间符合模块连接,时间符合模块设定的时间窗用于判断伽马光子探测器探测到的多个伽马光子事件是否构成多伽马光子符合事件,并将判断结果输入到计算机平台中;所有伽马光子探测器的能量和位置信号线还同时与计算机平台连接,用于多伽马光子符合事件的有效性判断以及伽马光子位置信息的计算;根据每一个多伽马光子符合事件在伽马光子探测器所输入的一个伽马光子的位置信息和缝孔混合准直器的准直缝位置信息,计算多伽马光子符合事件发生位置可能所在的若干投影平面,并根据每一个多伽马光子符合事件在伽马光子探测器所输入的一个伽马光子的位置信息和缝孔混合准直器的准直孔位置信息,计算多伽马光子符合事件发生位置可能所在的若干投影线;通过多伽马符合事件中的各个伽马光子所在投影平面和投影线的在成像范围内的若干个交点的计算,从而获得放射性核素发生衰变时的可能位置。根据所有放射性核素的可能衰变位置得到放射性核素在成像对象体内的可能分布,并可以通过图像重建算法对此分布做更为精确的估计。缝孔混合准直器是通过在各自合金板上分别开设若干条准直缝和若干个准直孔所形成,其中合金板厚度足以防止至少95%以上的伽马光子的穿透。合金板形状可以是平面、圆筒面、球面等各种形状。当缝孔混合准直器内相邻准直缝、准直孔之间的间距d满足以下公式时,可以使得在成像对象内的交点的个数缩减为一个,从而直接确定放射性核素发生衰变时的位置:
d≥D*a/(a+b)
其中,D是成像对象的大小,a是缝孔混合准直器表面到伽马光子探测器表面的距离,b是成像对象中心到缝孔混合准直器表面的距离。由于缝孔混合准直器采用过大的准直缝和准直孔间距会带来探测效率的下降,上式取等号时使得本成像系统能够直接确定放射性核素发生衰变时的位置的同时拥有最高的符合探测效率。
实施例1:
本实施例的成像系统总体结构如图1所示,其探测及定位的示意图如图2、3所示(为方便画图,缝孔混合准直器1仅画出部分准直孔和准直缝,图2和图3为采用不同排布的缝孔混合准直器1时成像系统的探测及定位示意图),该系统由互相垂直排列的两个缝孔混合准直器1、两个分别与各缝孔混合准直器1平行的伽马光子探测器2、时间符合模块3以及计算机平台4构成。缝孔混合准直器1是通过在一合金平板上开设若干个准直孔11和若干条准直缝12所形成,参见图4和图5,本实施例中,缝孔混合准直器1内各准直孔11大小相同且等间距设置,各条准直缝12大小相同、相互平行且等间距设置,形成缝孔混合准直阵列。缝孔混合准直器1设置于相应伽马光子探测器2前端使得成像对象5内放射性核素发生衰变产生的伽马光子仅沿缝孔混合准直器1的准直缝11和准直孔12到伽马光子探测器2方向发射才能被各伽马光子探测器2探测到,被探测到的伽马光子可以是以任意角度通过缝孔混合准直器1的准直孔11和准直缝12;成像对象5可以为生物体或其它被测物,也可以为核医学成像系统的标准化成像模型;伽马光子探测器2的时间信号线均连接到时间符合模块3,该时间符合模块设定一定宽度的时间窗(时间窗宽度根据所采用的放射性核素可调,通常在几百纳秒以内),用于判断分别被两个伽马光子探测器2探测到的两个伽马光子事件是否构成双伽马光子符合事件,并将相应的判断结果输入到计算机平台4中,具体地:若来自两个伽马光子探测器2的2个伽马光子的时间信号在时间窗内,则判断被两个伽马光子探测器2探测到的两个伽马光子事件构成双伽马光子符合事件,并将相应的判断结果输入到计算机平台4中;若来自两个伽马光子探测器2的2个伽马光子的时间信号不在时间窗内,则判断两个伽马光子探测器2探测到的两个伽马光子事件不构成双伽马光子符合事件,舍弃这两个事件;将两个伽马光子探测器2所测量得到的伽马光子的能量、位置信息一并传输到计算机平台4中;两个伽马光子探测器2测量得到的伽马光子的能量和位置信号在计算机平台4中用于多伽马光子符合事件发生位置所在投影平面7、投影线8的若干个交点计算;通过按照上述公式取等号的设计,将若干交点缩小为成像对象范围5内的一个交点,从而可以直接获取放射性核素的衰变可能位置6;根据所有放射性核素的可能衰变位置,可以得到放射性核素在成像对象体内的可能分布,并可以通过图像重建算法对此分布做更为精确的估计。
本实施例中缝孔混合准直器1到对应平行的伽马光子探测器2表面的距离均为a=15cm,成像对象中心到各缝孔混合准直器1表面的距离均为b=15cm,成像对象的尺寸为D=10cm,缝孔混合准直器1的相邻准直孔11、准直缝12间距均为5cm。
本实施例中所用缝孔混合准直器1(如图4、图5所示)均由矩形钨合金板组成,其中钨合金材料对伽马光子具有较强的吸收作用。在缝孔混合准直器1的矩形板上等间隔的交错设置了多个准直孔11和准直缝12,使得仅沿准直孔11和准直缝12发射的伽马光子才能穿过该缝孔混合准直器1被对应平行的伽马光子探测器2探测到。本实施例中所用缝孔混合准直器1的尺寸为336mm(长)×336mm(宽)×8mm(厚),每一个准直孔11的孔径为4.5mm,每一条准直缝12尺寸为312mm(长)×4.5mm(宽),相邻准直孔11、准直缝12之间的间距为50mm。
本发明成像系统所用缝孔混合准直器1所采用的钨合金板形状不局限于平面,也可以是闭合或不闭合的圆柱面(如图13所示)、球面等曲面,当采用闭合的圆柱面或球面的合金板时,可仅设置一个缝孔混合准直器1,至少设置1个伽马光子探测器2,且当采用多个伽马光子探测器2时,多个伽马光子探测器2分布在位于中心的缝孔混合准直器1的外围,优选地,相邻的伽马光子探测器2彼此接触形成围合中心的缝孔混合准直器1的空间。当采用非闭合的圆柱面或球面的合金板时,可设置至少一个缝孔混合准直器1,和至少一个伽马光子探测器2,缝孔混合准直器1和伽马光子探测器2的个数相同或不同。
此外,本发明成像系统对所用缝孔混合准直器1中准直孔11和准直缝12的形状和排布无限制,可以是任意形状和随机排布在合金板上。
本发明成像系统所用缝孔混合准直器1的准直孔11和准直逢12排布不局限于一定数量、相同尺寸、等间距以及特定朝向的排布,也可以是多种数量、不等间距和不同朝向的准直孔11和准直逢12排布(如图4、图5所示为不同准直孔11数量、准直逢12数量和排布的缝孔混合准直器1),可根据所要达到的探测效率、信噪比、空间分辨率、缝孔混合准直器1到对应平行的伽马光子探测器2表面的距离、成像对象中心到缝孔混合准直器1平面的距离等因素选择不同的缝孔混合准直器1的准直孔11、准直缝12的数量、尺寸、排布距离和排布方向。
本实施例中各伽马光子探测器2均为NaI(Tl)闪烁体探测器,所采用的闪烁体为一整块连续的NaI(Tl)晶体,其尺寸为585mm(长)×470mm(宽)×9.5mm(厚)。在NaI(Tl)晶体远离缝孔混合准直器1的一端设有耦合光电倍增管(Photomultiplier,以下简称PMT)或者硅光电倍增管(Silicon Photomultipliers,以下简称SiPM,又称Multi-Pixel PhotonCounter–MPPC),用于光电信号转换从而实现伽马光子在晶体内作用位置、能量和时间的测量。
本发明成像系统不局限于仅仅使用2个缝孔混合准直器1和2个伽马光子探测器2,也可以是至少1个缝孔混合准直器1和至少1个伽马光子探测器2的任意探测器组合,可根据所要达到的探测效率、信噪比、空间分辨率等因素选择不同的探头组合,如图6、图7所示为使用1个平板缝孔混合准直器1和1个伽马光子探测器2组成的成像系统;如图8、图9所示为使用4个平板缝孔混合准直器1和4个伽马光子探测器2组成的成像系统;如图10所示为使用5个平板缝孔混合准直器1和5个伽马光子探测器2组成的成像系统;如图11所示为使用6个平板缝孔混合准直器1和6个伽马光子探测器2组成的系统,如图12、图13所示为使用使用1个圆筒缝孔混合准直器1和6个伽马光子探测器2组成的系统。
本发明成像系统所用药物所标记的放射性核素除铟111外,还可以是其他多伽马光子放射核素,即该放射性核素在衰变过程中能以级方式在很短时间内产生至少两个伽马光子,包括但不限于镥177、钠22、碘131、铊201、铷82、钇90等。
本发明成像系统所用各探头的之间大部分是以非平行的方式排布的,由于大部分放射性核素的多伽马光子符合事件中的若干个伽马光子事件之间的发射夹角在90度处的分布概率最大,因此大部分情况下各探头呈互相垂直的分布为分布最优方案。
本发明成像系统的成像方法流程如图14所示,结合实施例1对该方法的具体实施步骤说明如下:
(1)启动所述成像系统,设定采集时间为20分钟,设定时间符合模块3的时间窗宽度为80ns;依据所用多伽马光子放射核素所放射的伽马光子能量分别设定各伽马光子探测器2所探测到伽马光子能量的能窗。本实施例中所用伽马光子放射核素为铟111,其在衰变过程中能以级联方式在很短的时间内(取决于所用放射性核素,通常在几百纳秒内)产生两个伽马光子,其能量分别为171keV和245keV,因此设定各伽马光子探测器2的两个能窗为一个171keV±10%能窗和一个245keV±10%能窗的组合;各伽马光子探测器2的能窗的个数均可根据所用放射性核素级联发射的伽马光子的个数调节,即每种能量的伽马光子对应一个能窗;每个能窗的宽度可根据各伽马光子探测器2的能量分辨率调节;向成像对象5注射活度为4mCi由铟111标记的放射性药物。
(2)计算机平台4根据所设定的采集时间判断成像过程是否结束;若结束,则执行步骤(5);若未结束,则执行步骤(3);
(3)时间符合模块3判断各伽马光子探测器2是否探测到双伽马光子符合事件,若探测到,则执行步骤(4);若未探测到,则执行步骤(2);所述双伽马光子符合事件,即如果被相应的两个伽马光子探测器2探测到的两个伽马光子事件在时间符合模块3所设定的时间窗内,则所探测的两个伽马光子事件构成双伽马光子符合事件;
(4)计算机平台4根据相应的两个伽马光子探测器2输入的两个伽马光子事件能量信息判断所输入的两个伽马光子能量是否分别在所设定的两个能窗内;如果不在所设定的能窗内,则舍弃该双伽马光子符合事件;如果在所设定的能窗内,则根据每一个多伽马光子符合事件在相应伽马光子探测器2所输入的一个伽马光子的位置信息和缝孔混合准直器1的准直缝12位置信息,计算多伽马光子符合事件发生位置可能所在的若干投影平面7,并根据每一个多伽马光子符合事件在相应伽马光子探测器2所输入的一个伽马光子的位置信息和缝孔混合准直器1的准直孔11位置信息,计算多伽马光子符合事件发生位置可能所在的若干投影线8;记录每一个双伽马光子符合事件中的两个伽马光子所分别确定的若干投影平面7和若干投影线8,计算得到每一个投影平面7和每一条投影线8之间在成像对象5内的一个交点6,记录该交点6位置放射性核素发生衰变的一个可能位置;执行步骤(2);
(5)根据所有双伽马光子符合事件计算的放射性核素的可能衰变位置,得到放射性核素在成像对象体内的可能分布,并可以运用最大似然最大期望值(MaximumLikelihood Expectation Maximization,简称MLEM)算法、有序子集最大期望值算法(Ordered Subsets Expectation Maximization,简称OSEM)等图像重建算法进行重建得到重建图像,对此分布做更为精确的估计。
本发明实施例中的成像系统由于能够通过直接计算方式获得放射性核素的一个可能衰变位置,并可以通过飞行时间(Time-of-Flight)测量技术等手段,对这一个可能衰变位置进行进一步分析与判断,从而简化了图像重建算法,提高了重建图像的信噪比。由于不需要累积大量的投影线来重建放射性核素的空间分布,与此同时使用了缝孔混合准直器,从而提高了单探测器探头的探测效率、避免了投影交点的位置判断简化了算法,进一步提高了多光子符合事件的探测效率,降低了对伽马光子事件总计数的需求,减少了所需要的放射性核素剂量。
将本发明所述成像方法进行编程(该过程可由编程人员通过常规的编程技术予以实现)后输入到计算机平台4中,按照步骤执行即可实现该成像方法预期达到的效果。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (6)

1.一种基于缝孔混合准直器的多伽马光子符合成像系统,其特征在于,包括时间符合模块、计算机平台、缝孔混合准直器和具备时间测量功能的伽马光子探测器;所述缝孔混合准直器是通过在合金板上开设若干条准直缝和若干个准直孔所形成,所述缝孔混合准直器分别放置在伽马光子探测器前端与成像对象之间并与两者保持一定距离,使得以任意方向通过所述缝孔混合准直器的由成像对象内放射性核素发生衰变产生的伽马光子事件均能被该伽马光子探测器探测到;所述伽马光子探测器的时间信号线与所述时间符合模块连接,时间符合模块设定的时间窗用于判断穿过所述缝孔混合准直器的准直缝和准直孔被伽马光子探测器探测到的多个伽马光子事件是否构成多伽马光子符合事件,并将判断结果输入到所述计算机平台中;所述伽马光子探测器的能量和位置信号线还同时与计算机平台连接,用于多伽马光子符合事件的有效性判断以及多伽马光子符合事件发生位置可能所在的若干投影线和若干投影线在成像范围内的若干个交点的计算,从而获得放射性核素发生衰变时的可能位置;根据所有放射性核素的可能衰变位置,得到放射性核素在成像对象体内的可能分布,并可以通过图像重建算法对此分布做更为精确的估计;
所述缝孔混合准直器内的合金板的形状为平面和曲面中的任意一种或者两种。
2.根据权利要求1所述的多伽马光子符合成像系统,其特征在于,所述曲面包括闭合或不闭合的圆柱面和球面。
3.根据权利要求1所述的多伽马光子符合成像系统,其特征在于,所述缝孔混合准直器内相邻准直孔、准直缝之间的间距d均满足:d≥D*a/(a+b),式中,D是成像对象的大小,a是缝孔混合准直器表面到伽马光子探测器表面的距离,b是成像对象中心到缝孔混合准直器表面的距离。
4.根据权利要求1所述的多伽马光子符合成像系统,其特征在于,所述缝孔混合准直器内各相邻准直孔、准直缝之间的间距相同或不同。
5.根据权利要求1所述的多伽马光子符合成像系统,其特征在于,所述放射性核素,在其衰变过程中能以级联方式在很短的时间内产生至少两个伽马光子。
6.一种采用如权利要求1~5中任一项所述成像系统的成像方法,其特征在于,该方法包括以下步骤:
(1)启动所述多伽马光子符合成像系统,设定时间符合模块的时间窗宽度;依据所用放射核素所放射的伽马光子能量设定伽马光子探测器的能窗;
(2)计算机平台根据所设定条件判断成像过程是否结束;若结束,则执行步骤(5);若未结束,则执行步骤(3);
(3)时间符合模块判断伽马光子探测器是否探测到多伽马光子符合事件,若探测到,则执行步骤(4);若未探测到,则执行步骤(2);
(4)计算机平台根据伽马光子探测器所输入的多个伽马光子能量信息判断所输入的多个伽马光子能量是否分别在所设定的能窗内;如果不在所设定的能窗内,则舍弃该多伽马光子符合事件;如在所设定的能窗内,则根据每一个多伽马光子符合事件在伽马光子探测器所输入的一个伽马光子的位置信息和缝孔混合准直器的准直缝位置信息,计算多伽马光子符合事件发生位置可能所在的若干投影平面,并根据每一个多伽马光子符合事件在伽马光子探测器所输入的一个伽马光子的位置信息和缝孔混合准直器的准直孔位置信息,计算多伽马光子符合事件发生位置可能所在的若干投影线;记录每一个多伽马光子符合事件中的多个伽马光子所分别确定的若干投影平面和若干投影线,计算每一个投影平面和每一条投影线之间在成像对象内的相应交点,作为放射性核素发生衰变的若干个可能位置;执行步骤(2);
(5)根据所有多伽马光子符合事件计算的放射性核素的多个可能衰变位置得到放射性核素在成像对象体内的可能分布,并可以通过图像重建算法对此分布做更为精确的估计。
CN202010465659.XA 2020-05-28 2020-05-28 基于缝孔混合准直器的多伽马光子符合成像系统及方法 Active CN111638544B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010465659.XA CN111638544B (zh) 2020-05-28 2020-05-28 基于缝孔混合准直器的多伽马光子符合成像系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010465659.XA CN111638544B (zh) 2020-05-28 2020-05-28 基于缝孔混合准直器的多伽马光子符合成像系统及方法

Publications (2)

Publication Number Publication Date
CN111638544A true CN111638544A (zh) 2020-09-08
CN111638544B CN111638544B (zh) 2023-07-04

Family

ID=72332378

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010465659.XA Active CN111638544B (zh) 2020-05-28 2020-05-28 基于缝孔混合准直器的多伽马光子符合成像系统及方法

Country Status (1)

Country Link
CN (1) CN111638544B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112419434A (zh) * 2020-11-04 2021-02-26 南京航空航天大学深圳研究院 一种γ光子3D成像噪声抑制方法及应用
CN113109855A (zh) * 2021-03-30 2021-07-13 北京科技大学 一种大视野放射源定位系统和方法
CN115825317A (zh) * 2022-11-12 2023-03-21 天津中核永泰科技有限公司 一种自动化放射性薄层色谱仪

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7339174B1 (en) * 2007-02-09 2008-03-04 General Electric Company Combined slit/pinhole collimator method and system
US20080111081A1 (en) * 2006-11-11 2008-05-15 Keh-Shih Chuang Imaging System And Method For The Non-Pure Positron Emission Tomography
US7375338B1 (en) * 2007-03-07 2008-05-20 General Electric Company Swappable collimators method and system
US20100294941A1 (en) * 2009-05-22 2010-11-25 Keh-Shih Chuang Dual Photons Emission Computed Tomography System
US20120061581A1 (en) * 2010-09-09 2012-03-15 Hugg James W Mixed resolution and multiplexing imaging method and system
CN106108934A (zh) * 2016-08-31 2016-11-16 清华大学 多伽马光子同时发射药物时间符合核医学成像系统及方法
US20170285191A1 (en) * 2016-03-29 2017-10-05 Kromek Goup, PLC Sparse Acquisition Gamma Cameras
CN108523916A (zh) * 2018-03-20 2018-09-14 清华大学 一种多伽马光子符合成像系统及方法
CN109683188A (zh) * 2019-01-15 2019-04-26 南昌大学 一种契连柯夫事件与伽马事件符合成像装置和方法
CN110772274A (zh) * 2019-10-29 2020-02-11 南昌大学 一种多核素的伽马成像系统与方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080111081A1 (en) * 2006-11-11 2008-05-15 Keh-Shih Chuang Imaging System And Method For The Non-Pure Positron Emission Tomography
US7339174B1 (en) * 2007-02-09 2008-03-04 General Electric Company Combined slit/pinhole collimator method and system
US7375338B1 (en) * 2007-03-07 2008-05-20 General Electric Company Swappable collimators method and system
US20100294941A1 (en) * 2009-05-22 2010-11-25 Keh-Shih Chuang Dual Photons Emission Computed Tomography System
US20120061581A1 (en) * 2010-09-09 2012-03-15 Hugg James W Mixed resolution and multiplexing imaging method and system
US20170285191A1 (en) * 2016-03-29 2017-10-05 Kromek Goup, PLC Sparse Acquisition Gamma Cameras
CN106108934A (zh) * 2016-08-31 2016-11-16 清华大学 多伽马光子同时发射药物时间符合核医学成像系统及方法
CN108523916A (zh) * 2018-03-20 2018-09-14 清华大学 一种多伽马光子符合成像系统及方法
CN109683188A (zh) * 2019-01-15 2019-04-26 南昌大学 一种契连柯夫事件与伽马事件符合成像装置和方法
CN110772274A (zh) * 2019-10-29 2020-02-11 南昌大学 一种多核素的伽马成像系统与方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
戴甜甜等: "高分辨率小动物SPECT成像系统的设计与性能评估", 《原子能科学技术》 *
金永杰等: "在SPECT上取得人体透射图像的方法", 《清华大学学报(自然科学版)》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112419434A (zh) * 2020-11-04 2021-02-26 南京航空航天大学深圳研究院 一种γ光子3D成像噪声抑制方法及应用
CN113109855A (zh) * 2021-03-30 2021-07-13 北京科技大学 一种大视野放射源定位系统和方法
CN115825317A (zh) * 2022-11-12 2023-03-21 天津中核永泰科技有限公司 一种自动化放射性薄层色谱仪
CN115825317B (zh) * 2022-11-12 2023-10-27 天津中核永泰科技有限公司 一种自动化放射性薄层色谱仪

Also Published As

Publication number Publication date
CN111638544B (zh) 2023-07-04

Similar Documents

Publication Publication Date Title
CN106108934B (zh) 多伽马光子同时发射药物时间符合核医学成像系统及方法
CN108523916B (zh) 一种多伽马光子符合成像系统及方法
US6791090B2 (en) Compton deconvolution camera
US7381960B1 (en) Imaging system and method for the non-pure positron emission tomography
CN111596336B (zh) 基于缝-孔平板准直器的多伽马光子符合成像系统及方法
CN111638544B (zh) 基于缝孔混合准直器的多伽马光子符合成像系统及方法
US8481947B2 (en) Method and system for nuclear imaging using multi-zone detector architecture
CN115299976B (zh) 多伽马光子符合成像系统及方法
US8809790B2 (en) Method and system for nuclear imaging using multi-zone detector architecture
Spanoudaki et al. Pet & SPECT instrumentation
NL2021303B1 (en) Active collimator system comprising a monolayer of monolithic converters
NL2020237B1 (en) Active collimator for positron emission and single photon emission computed tomography
Kijewski Positron emission tomography (PET) and single-photon emission computed tomography (SPECT) physics
Alnaaimi Evaluation of the UCL Compton camera imaging performance
Silva Small animal PET imaging using GATE Monte Carlo simulations: Implementation of physiological and metabolic information
Kijewski Positron emission tomography and single-photon emission computed tomography physics
RAHMIN PET vs. SPECT: in the context of ongoing developments
Beijst et al. Technical advances in image guidance of radionuclide therapy
Chen et al. Radiation Detection in SPECT and PET
Kim et al. Imaging and Sensing of pH and Chemical State With Cascade Photons of 111In Using Ring-Type Compton Camera
Newton Entangled polarisation correlations of annihilation gamma and their applications to PET Imaging
Chiesa et al. Advances in nuclear medicine imaging in oncology
Talat Investigation of the effect of crystal thickness on the partial resolution and linearity of a pem detector using an artificial neural network based positioning algorithm
Celler Single photon imaging and instrumentation
Surti et al. Biomedical images: signals to pictures

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