CN108523916B - 一种多伽马光子符合成像系统及方法 - Google Patents
一种多伽马光子符合成像系统及方法 Download PDFInfo
- Publication number
- CN108523916B CN108523916B CN201810230414.1A CN201810230414A CN108523916B CN 108523916 B CN108523916 B CN 108523916B CN 201810230414 A CN201810230414 A CN 201810230414A CN 108523916 B CN108523916 B CN 108523916B
- Authority
- CN
- China
- Prior art keywords
- gamma photon
- gamma
- detector
- coincidence
- compton camera
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 18
- 239000000523 sample Substances 0.000 claims abstract description 76
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 230000002285 radioactive effect Effects 0.000 claims description 16
- 230000000694 effects Effects 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 4
- 238000003325 tomography Methods 0.000 abstract description 4
- 239000013078 crystal Substances 0.000 description 16
- 238000002603 single-photon emission computed tomography Methods 0.000 description 9
- 229940079593 drug Drugs 0.000 description 6
- 239000003814 drug Substances 0.000 description 6
- 238000009206 nuclear medicine Methods 0.000 description 6
- 238000002600 positron emission tomography Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000008021 deposition Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- APFVFJFRJDLVQX-AHCXROLUSA-N indium-111 Chemical compound [111In] APFVFJFRJDLVQX-AHCXROLUSA-N 0.000 description 3
- 229940055742 indium-111 Drugs 0.000 description 3
- 229910001080 W alloy Inorganic materials 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000008358 core component Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- PNDPGZBMCMUPRI-HVTJNCQCSA-N 10043-66-0 Chemical compound [131I][131I] PNDPGZBMCMUPRI-HVTJNCQCSA-N 0.000 description 1
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 208000012902 Nervous system disease Diseases 0.000 description 1
- 208000025966 Neurological disease Diseases 0.000 description 1
- IGLNJRXAVVLDKE-OIOBTWANSA-N Rubidium-82 Chemical compound [82Rb] IGLNJRXAVVLDKE-OIOBTWANSA-N 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- VWQVUPCCIRVNHF-OUBTZVSYSA-N Yttrium-90 Chemical compound [90Y] VWQVUPCCIRVNHF-OUBTZVSYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- QWUZMTJBRUASOW-UHFFFAOYSA-N cadmium tellanylidenezinc Chemical compound [Zn].[Cd].[Te] QWUZMTJBRUASOW-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000003631 expected effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000000926 neurological effect Effects 0.000 description 1
- 229940121896 radiopharmaceutical Drugs 0.000 description 1
- 239000012217 radiopharmaceutical Substances 0.000 description 1
- 230000002799 radiopharmaceutical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- KEAYESYHFKHZAL-BJUDXGSMSA-N sodium-22 Chemical compound [22Na] KEAYESYHFKHZAL-BJUDXGSMSA-N 0.000 description 1
- BKVIYDNLLOSFOA-OIOBTWANSA-N thallium-201 Chemical compound [201Tl] BKVIYDNLLOSFOA-OIOBTWANSA-N 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/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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种多伽马光子符合成像系统及方法，属于发射断层成像技术领域，该系统包括时间符合模块、计算机平台、至少一个由准直器和伽马光子探测器构成的第一探头，至少一个由前后两个康普顿相机探测器构成的第二探头，各探头探测放射性核素放射的多个伽马光子构成多伽马光子符合事件；该成像方法将放射性核素发生衰变的位置范围缩小为多伽马光子符合事件中第一探头探测到的伽马光子事件所确定的投影线和第二探头探测到的伽马光子事件所确定的投影圆锥面的若干个交点，积累一定数量的多伽马光子符合事件即可获得放射性核素在被测范围内分布的图像。本发明简化了重建算法，提高了探测效率以及重建图像的信噪比，降低了对伽马光子总计数的需求。
Description
技术领域
本发明属于发射断层成像技术领域，特别涉及一种多伽马光子符合成像系统及方法。
背景技术
发射断层成像技术是当今探测物体内部信息的重要技术之一，在多研究领域中有着大量的应用。在被探测物体外通过探测从被探测物体发出的伽马光子，并通过图像重建从而以非侵入式的手段观察得到被探测物体的内部信息。在发射断层成像领域，最为重要的几种成像系统分别有正电子发射断层成像仪(Positron Emission Tomography,以下简称PET)、单光子发射计算机断层成像仪(Single-Photon Emission Computed Tomography,以下简称SPECT)和康普顿相机(Compton Camera)等。目前PET和SPECT均已广泛用于临床检查和诊断，包括癌症、神经系统疾病和心血管疾病等，而康普顿相机也应用于核能工业、天文、医疗等方面。
PET的核心组件由很多个具备时间测量功能的伽马光子探测器模块以及相应的时间符合模块构成。PET的基本原理涉及电子准直技术，所利用的放射性核素为正电子核素。正电子核素放射出的正电子会在被探测物体内发生正负电子湮灭效应，产生一对能量为511keV、方向几乎相反的伽马光子。采用时间符合测量手段，即如果两个伽马光子探测器模块在一个很短的时间内(通常在几百个纳秒)分别探测到这两个能量为511keV的伽马光子，就能确定一条正电子湮灭的发生位置(近似为正电子核素发生衰变的位置)所在的响应线(Line ofResponse)。记录大量这样的响应线通过图像重建即可获得正电子湮灭发生位置的分布，即近似为正电子核素在被探测物体内的分布。由于正电子湮灭产生的一对伽马光子的发射方向是几乎相反的，因此只能确定正电子湮灭的发生位置在相应的响应线上，而不能确定正电子湮灭发生在相应响应线上的具体位置。尽管可以通过飞行时间测量技术初步确定正电子湮灭发生在相应响应线上的位置范围，但这需要伽马光子探测器模块具备极高的时间分辨率。由于正电子湮灭在响应线上发生位置的不确定性，重建的正电子核素在被探测物体内分布图像的信噪比往往较低，影响诊断的效果。为提高图像信噪比通常需要累积大量的响应线，而这会使得被探测物体摄入较大剂量的正电子核素，增加被探测物体的辐照风险。
SPECT的核心组件包括准直器及伽马光子探测器模块等。SPECT利用的是物理准直技术，其所利用的核素为伽马光子核素。在伽马光子探测器模块的前端通常会放置准直器来限制伽马光子核素放射出的伽马光子到达探测器的入射角度，使得仅沿特定方向发射的伽马光子才能通过准直器被探测器探测到，伽马光子探测器每探测到一个伽马光子就能确定一条伽马光子初始发射位置所在的投影线。积累大量这样的投影线通过图像重建就能确定伽马光子初始发射位置的分布，即伽马光子核素在被探测物体内的分布。与PET存在类似的缺陷，SPECT也不能确定伽马光子在投影线上的具体发射位置，因此重建图像的信噪比较差。另外由于SPECT使用了准直器，限制了能被探测器探测到的伽马光子的发射角度，成像系统的探测效率较低，而这进一步使得重建图像的信噪比恶化。
康普顿相机模块的核心组件包括两块平行的探测器平面模块等。采用时间符合测量手段，即如果两个平行的探测器平面模块在一个很短的时间内(通常在几个纳秒)分别探测到信号，则可以认为该两次信号来自于同一次伽马光子事件。伽马光子在第一块晶体平面上发生康普顿散射产生康普顿散射光子，所产生的康普顿散射光子在第二块晶体平面发生光电效应被吸收，康普顿相机从而得以探测到该次伽马光子事件。根据伽马光子在第一块晶体平面上沉积能量与第二块晶体平面上沉积能量之和，可以得到伽马光子的总能量。康普顿相机每探测到一个伽马光子事件，根据伽马光子事件在第一块晶体平面上的沉积能量、沉积位置和第二块晶体平面上的沉积能量、沉积位置，就能根据康普顿散射原理确定一个伽马光子初始发射位置所在的投影圆锥面。由于康普顿相机探测器不存在准直器设计，相比于SPECT伽马探测器探测效率将会大大提高，但其仍不能直接得到伽马光子的发射位置分布。除此之外，康普顿相机空间分辨率较差，对于康普顿效应较少发生的低能以及高能伽马光子探测效率较低。
本申请人已提出的一种多伽马光子同时发射药物时间符合核医学成像系统及方法(申请号：201610798146.4)，该系统包括多个以非平行方式排列的探测器探头、时间符合模块以及计算机平台构成，每个探测器探头均由准直器和具备时间测量功能的伽马光子探测器构成，探测放射性核素在很短时间内放射的多个伽马光子构成多伽马光子符合事件；所述方法计算到多伽马光子符合事件中的每一个伽马光子事件所确定的投影线的距离之和最短的点的位置即为放射性核素发生衰变的位置，积累一定数量的多伽马光子符合事件即可实现放射性核素在生物体内分布的获取。该成像系统和方法简化了重建算法，提高了重建图像的信噪比；降低了对伽马光子总计数的需求，降低了病人的辐照风险。但是，由于该系统每一个探测器探头均包含限制光子方向的准直器模块，多光子符合事件探测效率极低，导致探测数据信噪比低。
发明内容
本发明的目的旨在从原理上解决上述PET系统、SPECT系统和已提出的一种多伽马光子同时发射药物时间符合核医学成像系统的问题，本发明公开了一种新的多伽马光子符合成像系统及方法。与现有的多伽马光子同时发射药物时间符合核医学成像技术的不同之处在于，本发明使用至少一个康普顿相机探测器探头代替伽马探测器探头，从而提高了多光子符合事件的探测效率，提高了重建图像的信噪比。
为了实现上述目的，本发明采用如下技术方案：
本发明提出的一种多伽马光子符合成像系统，包括至少一个伽马光子探测器探头、时间符合模块和计算机平台；各伽马光子探测器探头均由准直器和具备时间测量功能的伽马光子探测器构成，各准直器放置在与其对应的伽马光子探测器前端使得成像对象内放射性核素发生衰变产生的伽马光子仅沿设定方向发射才能被伽马光子探测器探测到；该系统还包括至少一个康普顿相机探测器探头；各康普顿相机探测器探头均由前后两个具备时间测量功能的康普顿相机探测器构成，使得成像对象内放射性核素发生衰变产生的伽马光子在前端康普顿相机探测器中发生康普顿散射，产生的康普顿散射光子在后端康普顿相机探测器发生光电效应被收集，且该康普顿相机探测器探头设定一定宽度的时间窗，用于判断其内两个康普顿相机探测器上探测到的事件是否为同一个伽马光子事件；所有探测器的时间信号线均与所述时间符合模块连接，时间符合模块设定的时间窗用于判断分别被多个探测器探头探测到的多个伽马光子事件是否构成多伽马光子符合事件，并将判断结果输入到计算机平台中；所有探测器的能量和位置信号线还同时与计算机平台连接，用于多伽马光子符合事件的有效性判断以及伽马光子发射位置所在投影线和投影圆锥面交点的计算，从而获得放射性核素发生衰变时的可能位置。
本发明还提出一种采用如上述成像系统的成像方法，该方法包括以下步骤：
(1)启动所述多伽马光子符合成像系统，设定时间符合模块的时间窗宽度；依据所用多伽马光子放射核素所放射的伽马光子能量设定各探测器的能窗；向成像对象引入放射性核素；
(2)计算机平台根据所设定条件判断成像过程是否结束；若结束，则执行步骤(5)；若未结束，则执行步骤(3)；
(3)时间符合模块判断各探测器是否探测到多伽马光子符合事件，若探测到，则执行步骤(4)；若未探测到，则执行步骤(2)；
(4)计算机平台根据各探测器探头所输入的多个伽马光子事件能量信息判断所输入的多个伽马光子能量是否分别在所设定的能窗内；如果不在所设定的能窗内，则舍弃该多伽马光子符合事件；如在所设定的能窗内，则根据每一个伽马光子探测器所输入的一个伽马光子事件的位置信息计算伽马光子发生衰变所在位置的一条投影线，并根据每一个康普顿相机探测器探头所输入的一个伽马光子事件的位置信息和能量信息计算伽马光子发生衰变所在位置的一个投影圆锥面；计算多个伽马光子事件所分别确定的若干条投影线和若干投影圆锥面是否相交，如果相交有若干个交点，记录为放射性核素发生衰变的若干个可能位置；如果若干投影线和若干个投影圆锥面不相交，则舍弃该多伽马光子符合事件；执行步骤(2)；
(5)根据所有多伽马光子符合事件计算的放射性核素的多个可能衰变位置，运用图像重建算法得到重建图像，作为放射性核素在成像对象体内的可能分布。
本发明的特点及有益效果：
本发明的基于多伽马光子符合成像系统，有效地克服了传统PET系统或SPECT系统只能确定放射性核素衰变位置所在响应线或投影线而不能确定放射性核素衰变发生在响应线或投影线上具体位置的缺陷和已提出的一种多伽马光子同时发射药物时间符合核医学成像系统探测效率过低、信噪比过低的缺陷，利用其与康普顿相机相结合，通过计算到多个伽马光子所确定的若干条投影线和若干个投影圆锥面的若干个交点实现对放射性核素的衰变位置的定位，从而获得放射性核素在被探测物体内的可能分布。由于能够基于平行投影线和投影圆锥面的交点计算放射性核素的可能衰变位置，简化了图像重建算法，提高了重建图像的信噪比。同时由于不需要累积大量的投影线来重建放射性核素的空间分布，采用了探测效率高的康普顿相机探测器，从而降低了对伽马光子事件总计数的需求，减少了需要的放射性核素剂量。
附图说明
图1是本发明实施例在伽马探测器探头中采用平行孔准直器成像系统的整体结构示意图；
图2是本发明实施例在伽马探测器探头中采用针孔准直器成像系统的整体结构示意图；
图3是本发明实施例采用2个伽马探测器探头和1个康普顿相机探测器探头的成像系统整体结构示意图；
图4是本发明实施例采用1个伽马探测器探头和2个康普顿相机探测器探头的成像系统整体结构示意图；
图5是本发明成像方法的流程框图。
具体实施方式
本发明提出的一种多伽马光子符合成像系统及方法结合附图以及实施例详细说明如下：
本发明提出的一种多伽马光子符合成像系统，包括至少一个伽马光子探测器探头、至少一个康普顿相机探测器探头、时间符合模块以及计算机平台；各伽马光子探测器探头均由准直器和具备时间测量功能的伽马光子探测器构成，各准直器放置在与其对应的伽马光子探测器前端使得成像对象内放射性核素发生衰变产生的伽马光子仅沿设定方向发射才能被伽马光子探测器探测到；各康普顿相机探测器探头均由前后两个高探测效率的具备时间测量功能的康普顿相机探测器构成，使得成像对象内放射性核素发生衰变产生的伽马光子在前端康普顿相机探测器中发生康普顿散射，产生的康普顿散射光子在后端康普顿相机探测器发生光电效应被收集，且该康普顿相机探测器探头设定一定宽度的时间窗，用于判断其内两个康普顿相机探测器上探测到的事件是否为同一个伽马光子事件；其中，所有的探测器的时间信号线均与时间符合模块连接，时间符合模块设定的时间窗用于判断分别被多个探测器探头探测到的多个伽马光子事件是否构成多伽马光子符合事件，并将判断结果输入到计算机平台中；所有的探测器的能量和位置信号线还同时与计算机平台连接，用于多伽马光子符合事件的有效性判断以及伽马光子发射位置所在投影线和投影圆锥面的交点的计算，从而获得放射性核素发生衰变时的可能位置。
实施例1：
本实施例的成像系统总体结构如图1所示，该系统由探测平面互相垂直排列的伽马探测器探头和康普顿相机探测器探头5各一个、时间符合模块6以及计算机平台7构成；伽马探测器探头由平行孔准直器1和具备时间测量功能的伽马光子探测器2构成，平行孔准直器1放置在伽马光子探测器2前端使得成像对象8内放射性核素发生衰变产生的伽马光子仅沿垂直于伽马光子探测器2面方向发射才能被伽马光子探测器2探测到；康普顿相机探测器探头5由具备时间测量功能的前端康普顿相机探测器3和后端康普顿相机探测器4构成，成像对象8内放射性核素发生衰变产生的伽马光子在前端康普顿相机探测器3中发生康普顿散射，产生的康普顿散射光子在后端康普顿相机探测器4发生光电效应被收集，该康普顿相机探测器探头5设定一定宽度的时间窗，用于判断前端康普顿相机探测器3和后端康普顿相机探测器4上探测到的事件是否为同一个伽马光子事件；成像对象8可以为生物体或其它被测物，也可以为核医学成像系统的标准化成像模型；伽马光子探测器2和康普顿相机探测器探头5的时间信号线均连接到时间符合模块6，该模块设定一定宽度的时间窗(时间窗宽度根据所采用的放射性核素可调，通常在几百纳秒以内)，用于判断分别被两个探测器探头探测到的两个伽马光子事件是否构成双伽马光子符合事件，并将相应的判断结果输入到计算机平台7中，具体地：若来自伽马光子探测器2和康普顿相机探测器探头5的时间信号在时间窗内，则判断分别被伽马光子探测器2和康普顿相机探测器探头5探测到的两个伽马光子事件构成双伽马光子符合事件，并将相应的判断结果输入到计算机平台7中；若来自伽马光子探测器2和康普顿相机探测器探头5的时间信号不在时间窗内，则判断分别被伽马光子探测器2和康普顿相机探测器探头5探测到的两个伽马光子事件不构成双伽马光子符合事件，舍弃这两个事件；在康普顿相机探测器探头5中将前端康普顿相机探测器3所探测的能量和后端康普顿相机探测器4探测的能量求和作为所探测的伽马光子能量，并将此能量信息和前端康普顿相机探测器3、后端康普顿相机探测器4的能量、位置信息一并传输到计算机平台7中；伽马光子探测器2和康普顿相机探测器探头5的能量和位置信号在计算机平台7中用于伽马光子发射位置所在投影线、投影圆锥面的交点计算，从而可以获取放射性核素的两个衰变可能位置。
本实施例中所用平行孔准直器1由矩形钨合金板组成，其中钨合金材料对伽马光子具有较强的吸收作用。在该矩形板上等间隔地设置了多个平行的准直孔，使得仅沿该孔发射的伽马光子才能穿过准直器被伽马光子探测器2探测到。本实施例中所用准直器的厚度为4mm,准直孔的孔径为2mm，孔壁的厚度为2mm。
本实施例中伽马光子探测器2为NaI(Tl)闪烁体探测器，所采用闪烁体为一整块连续的NaI(Tl)晶体，其尺寸为585mm(长)×470mm(宽)×9.5mm(厚)。在NaI(Tl)晶体远离准直器的一端耦合光电倍增管(Photomultiplier,以下简称PMT)或者硅光电倍增管(SiliconPhotomultipliers，以下简称SiPM,又称Multi-Pixel Photon Counter–MPPC)，用于光电信号转换从而实现伽马光子在晶体内作用位置、能量和时间的测量。
本实施例中前端康普顿相机探测器3为碲锌镉(CadmiumZinc Telluride,以下简称CZT)探测器，其尺寸为300mm(长)×300mm(宽)×2mm(厚)。后端康普顿相机探测器4为CZT探测器或者NaI(Tl)闪烁体探测器。其中，CZT探测器采用尺寸为300mm(长)×300mm(宽)×5mm(厚)；NaI(Tl)闪烁体探测器尺寸为300mm(长)×300mm(宽)×5mm(厚)，所采用闪烁体为一整块连续的NaI(Tl)晶体，其尺寸为585mm(长)×470mm(宽)×9.5mm(厚)。如果后端康普顿相机探测器4采用NaI(Tl)晶体，则在NaI(Tl)晶体远离准直器的一端耦合PMT或者SiPM，用于光电信号转换从而实现伽马光子在晶体内作用位置、能量和时间的测量。前端康普顿相机探测器3和后端康普顿相机探测器4之间的符合时间窗应小于5ns，确保两康普顿相机探测器上探测到的事件来自同一伽马光子事件。
本发明成像系统所用药物所标记的放射性核素除铟111外，还可以是其他多伽马光子放射核素，即该放射性核素在衰变过程中能以级方式在很短时间内产生至少两个伽马光子，包括但不限于钠22、碘131、铊201、铷82、钇90等。
本发明成像系统所用准直器的种类除图1所展示的平行孔准直器外，还可以是针孔型准直器(如图2所示)、汇聚型准直器、发散型准直器等。所采用的准直器均为常规产品，可根据所要达到的成像系统成像视野大小、空间分辨率及探测效率等因素选择所用准直器类型和参数。
本发明成像系统所用探测器不局限于仅仅使用一个伽马探测器探头和一个康普顿相机探测器探头，也可以是至少一个伽马探测器探头和至少一个康普顿相机探测器探头(如图3所示为两个伽马探测器探头和一个康普顿相机探测器探头，图4所示为一个伽马探测器探头和两个康普顿相机探测器探头)的任意探测器探头组合(至少包含两个探测器，其中至少一个为康普顿相机探测器)，可根据所要达到的探测效率、信噪比、空间分辨率等因素选择不同的探测器组合。
本发明成像系统的成像方法流程如图5所示，结合实施例1对该方法的具体实施步骤说明如下：
(1)启动所述成像系统，设定采集时间为20分钟，设定时间符合模块6的时间窗宽度为80ns；依据所用伽马光子放射核素所放射的伽马光子能量分别设定伽马光子探测器2所探测到伽马光子能量和康普顿相机探测器探头5所探测到伽马光子能量的能窗。本实施例中所用伽马光子放射核素为铟111，其在衰变过程中能以级联方式在很短的时间内(取决于所用放射性核素，通常在几百纳秒内)产生两个伽马光子，其能量分别为171keV和245keV，因此设定伽马光子探测器2和康普顿相机探测器探头5的两个能窗为一个171keV±10％能窗和一个245keV±10％能窗的组合；伽马光子探测器2和康普顿相机探测器探头5的能窗的个数均可根据所用放射性核素级联发射的伽马光子的个数调节，即每种能量的伽马光子对应一个能窗；每个能窗的宽度可根据伽马光子探测器2和康普顿相机探测器探头5的能量分辨率调节；向成像对象8注射活度为4mCi由铟111标记的放射性药物。
(2)计算机平台7根据所设定的采集时间判断成像过程是否结束；若结束，则执行步骤(5)；若未结束，则执行步骤(3)；
(3)时间符合模块6判断伽马光子探测器2和康普顿相机探测器探头5是否探测到双伽马光子符合事件，若探测到，则执行步骤(4)；若未探测到，则执行步骤(2)；所述双伽马光子符合事件，即如果被两个探测器探头探测到的两个伽马光子事件在时间符合模块6所设定的时间窗内，则所探测的两个伽马光子事件构成双伽马光子符合事件；
(4)计算机平台7根据两个探测器探头输入的两个伽马光子事件能量信息判断所输入的两个伽马光子能量是否分别在所设定的两个能窗内；如果不在所设定的能窗内，则舍弃该双伽马光子符合事件；如果在所设定的能窗内，则根据伽马光子探测器2所输入的一个伽马光子事件的位置信息计算伽马光子发生衰变所在位置的投影线，根据康普顿相机探测器探头5所探测到的伽马光子事件在前端康普顿相机探测器3和后端康普顿相机探测器4中的能量信息和位置信息计算伽马光子发生衰变所在位置的投影圆锥面；计算到两个伽马光子事件所分别确定的一条投影线和一条投影圆锥面是否相交，如果相交有两个交点，记录该两交点位置放射性核素发生衰变的两个可能位置；如果投影线和投影圆锥面不相交，则舍弃该双伽马光子符合事件；执行步骤(2)；
(5)根据所有双伽马光子符合事件计算的放射性核素的可能衰变位置，运用最大似然最大期望值(Maximum Likelihood Expectation Maximization，简称MLEM)算法、有序子集最大期望值算法(Ordered Subsets Expectation Maximization，简称OSEM)算法等图像重建算法进行重建，得到重建图像，作为放射性核素在成像对象体内的可能分布。
本发明实施例中的成像系统由于能够通过直接计算方式获得放射性核素的两个可能衰变位置，并可以通过飞行时间(Time-of-Flight)测量技术等手段，对这两个可能衰变位置进行进一步分析与判断，从而简化了图像重建算法，提高了重建图像的信噪比。同时由于不需要累积大量的投影线来重建放射性核素的空间分布，使用康普顿相机也大大提高了双伽马光子符合事件的探测效率，降低了对伽马光子事件总计数的需求，减少了所需要的放射性核素剂量。
将本发明所述成像方法进行编程(该过程可由编程人员通过常规的编程技术予以实现)后输入到计算机平台7中，按照步骤执行即可实现该成像方法预期达到的效果。
Claims (4)
1.一种多伽马光子符合成像系统，包括至少一个伽马光子探测器探头、时间符合模块和计算机平台；各伽马光子探测器探头均由准直器和具备时间测量功能的伽马光子探测器构成，各准直器放置在与其对应的伽马光子探测器前端使得成像对象内放射性核素发生衰变产生的伽马光子仅沿设定方向发射才能被伽马光子探测器探测到；其特征在于，该系统还包括至少一个康普顿相机探测器探头；各康普顿相机探测器探头均由前后两个具备时间测量功能的康普顿相机探测器构成，使得成像对象内放射性核素发生衰变产生的伽马光子在前端康普顿相机探测器中发生康普顿散射，产生的康普顿散射光子在后端康普顿相机探测器发生光电效应被收集，且该康普顿相机探测器探头设定一定宽度的时间窗，用于判断其内两个康普顿相机探测器上探测到的事件是否为同一个伽马光子事件；所有探测器的时间信号线均与所述时间符合模块连接，时间符合模块设定的时间窗用于判断分别被多个伽马光子探测器探头和康普顿相机探测器探头探测到的多个伽马光子事件是否构成多伽马光子符合事件，并将判断结果输入到计算机平台中；所有探测器的能量和位置信号线还同时与计算机平台连接，用于多伽马光子符合事件的有效性判断以及伽马光子发射位置所在投影线和投影圆锥面交点的计算，从而获得放射性核素发生衰变时的可能位置。
2.如权利要求1所述的多伽马光子符合成像系统，其特征在于，所述伽马光子探测器探头的探测平面和康普顿相机探测器探头的探测平面以非平行的方式排列布置。
3.如权利要求1所述的多伽马光子符合成像系统，其特征在于，所述放射性核素，在其衰变过程中能以级联方式在很短的时间内产生至少两个伽马光子。
4.一种采用如权利要求1所述成像系统的成像方法，其特征在于，该方法包括以下步骤：
(1)启动所述多伽马光子符合成像系统，设定时间符合模块的时间窗宽度；依据所用多伽马光子放射性 核素所放射的伽马光子能量设定各探测器的能窗；
(2)计算机平台根据所设定条件判断成像过程是否结束；若结束，则执行步骤(5)；若未结束，则执行步骤(3)；
(3)时间符合模块判断各探测器是否探测到多伽马光子符合事件，若探测到，则执行步骤(4)；若未探测到，则执行步骤(2)；
(4)计算机平台根据各探测器探头所输入的多个伽马光子事件能量信息判断所输入的多个伽马光子能量是否分别在所设定的能窗内；如果不在所设定的能窗内，则舍弃该多伽马光子符合事件；如在所设定的能窗内，则根据每一个伽马光子探测器所输入的一个伽马光子事件的位置信息计算伽马光子发生衰变所在位置的一条投影线，并根据每一个康普顿相机探测器探头所输入的一个伽马光子事件的位置信息和能量信息计算伽马光子发生衰变所在位置的一个投影圆锥面；计算多个伽马光子事件所分别确定的若干条投影线和若干投影圆锥面是否相交，如果相交有若干个交点，记录为放射性核素发生衰变的若干个可能位置；如果若干投影线和若干个投影圆锥面不相交，则舍弃该多伽马光子符合事件；执行步骤(2)；
(5)根据所有多伽马光子符合事件计算的放射性核素的多个可能衰变位置，运用图像重建算法得到重建图像，作为放射性核素在成像对象体内的可能分布。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810230414.1A CN108523916B (zh) | 2018-03-20 | 2018-03-20 | 一种多伽马光子符合成像系统及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810230414.1A CN108523916B (zh) | 2018-03-20 | 2018-03-20 | 一种多伽马光子符合成像系统及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108523916A CN108523916A (zh) | 2018-09-14 |
| CN108523916B true CN108523916B (zh) | 2020-07-28 |
Family
ID=63484215
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810230414.1A Active CN108523916B (zh) | 2018-03-20 | 2018-03-20 | 一种多伽马光子符合成像系统及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108523916B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109782327B (zh) * | 2019-01-08 | 2020-12-25 | 中国科学院高能物理研究所 | 一种康普顿散射-编码孔径复合成像方法 |
| CN109872269B (zh) * | 2019-01-08 | 2020-07-21 | 中国科学院高能物理研究所 | 一种基于鱼眼图像校正的康普顿相机图像融合方法 |
| CN110570929B (zh) * | 2019-08-20 | 2022-11-25 | 上海联影医疗科技股份有限公司 | 结构展示方法、装置、计算机设备和存储介质 |
| CN110664423B (zh) * | 2019-09-12 | 2023-05-09 | 沈阳智核医疗科技有限公司 | 成像方法、装置、探测器、终端设备和pet系统 |
| CN110680368B (zh) * | 2019-09-12 | 2023-03-24 | 东软医疗系统股份有限公司 | 有效单事件的检测方法及装置、信号读出方法及系统 |
| EP3845929B1 (en) * | 2020-01-06 | 2024-01-03 | Canon Medical Systems Corporation | A pet apparatus and a method |
| CN111596336B (zh) * | 2020-04-08 | 2022-06-10 | 清华大学 | 基于缝-孔平板准直器的多伽马光子符合成像系统及方法 |
| CN111638544B (zh) * | 2020-05-28 | 2023-07-04 | 清华大学 | 基于缝孔混合准直器的多伽马光子符合成像系统及方法 |
| CN113069138B (zh) * | 2021-03-23 | 2023-06-30 | 上海联影医疗科技股份有限公司 | 正电子发射断层成像装置、符合效率检测方法及归一化方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8153986B2 (en) * | 2007-07-09 | 2012-04-10 | Lawrence Livermore National Security, Llc | Hybrid Compton camera/coded aperture imaging system |
| JP5582370B2 (ja) * | 2010-12-09 | 2014-09-03 | 独立行政法人理化学研究所 | ガンマ線を利用する画像化装置、画像信号処理装置およびガンマ線測定データの画像処理方法 |
| EP2994776B1 (en) * | 2013-05-08 | 2019-11-13 | Koninklijke Philips N.V. | Apparatus and method for the evaluation of gamma radiation events |
-
2018
- 2018-03-20 CN CN201810230414.1A patent/CN108523916B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108523916A (zh) | 2018-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108523916B (zh) | 一种多伽马光子符合成像系统及方法 | |
| CN106108934B (zh) | 多伽马光子同时发射药物时间符合核医学成像系统及方法 | |
| US6791090B2 (en) | Compton deconvolution camera | |
| Moskal et al. | Novel detector systems for the positron emission tomography | |
| US20030197128A1 (en) | Method and apparatus for gamma ray detection | |
| Larobina et al. | Small animal PET: a review of commercially available imaging systems | |
| CN111596336B (zh) | 基于缝-孔平板准直器的多伽马光子符合成像系统及方法 | |
| CN111638544B (zh) | 基于缝孔混合准直器的多伽马光子符合成像系统及方法 | |
| CN112470039B (zh) | 用于通过伽马辐射探测进行成像的系统和方法 | |
| US8809790B2 (en) | Method and system for nuclear imaging using multi-zone detector architecture | |
| Spanoudaki et al. | Pet & SPECT instrumentation | |
| US20050205791A1 (en) | Method and apparatus for vetoing random coincidences in positron emission tomographs | |
| NL2021303B1 (en) | Active collimator system comprising a monolayer of monolithic converters | |
| Kijewski | Positron emission tomography (PET) and single-photon emission computed tomography (SPECT) physics | |
| Alnaaimi | Evaluation of the UCL Compton camera imaging performance | |
| Kijewski | Positron emission tomography and single-photon emission computed tomography physics | |
| Silva | Small animal PET imaging using GATE Monte Carlo simulations: Implementation of physiological and metabolic information | |
| Chen et al. | Radiation Detection in SPECT and PET | |
| Dahlbom et al. | Principles of SPECT and PET imaging | |
| Shahi | Evaluation of multi-angle imaging with planar organ-targeted positron emission tomography detectors | |
| Wilson | High-Sensitivity PET using Optimised Continuous Cylindrical Shell Nanocomposite and Transparent Ceramic Scintillators | |
| Borghi | From detectors towards systems: enabling clinical TOF-PET with monolithic scintillators | |
| Tayefi Ardebili | Evaluation of the NEMA characteristics for the Modular J-PET scanner | |
| Vallabhajosula | Radioactivity Detection: PET and SPECT Scanners | |
| Andreas | Novel High Resolution Photon Detectors for PET Imaging |
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 |