CN107544086B - 加马光子侦测成像装置与方法 - Google Patents
加马光子侦测成像装置与方法 Download PDFInfo
- Publication number
- CN107544086B CN107544086B CN201710448691.5A CN201710448691A CN107544086B CN 107544086 B CN107544086 B CN 107544086B CN 201710448691 A CN201710448691 A CN 201710448691A CN 107544086 B CN107544086 B CN 107544086B
- Authority
- CN
- China
- Prior art keywords
- scintillation crystal
- layer
- detection
- time
- photons
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000003384 imaging method Methods 0.000 title description 36
- 238000001514 detection method Methods 0.000 claims abstract description 164
- 239000000523 sample Substances 0.000 claims abstract description 147
- 239000013078 crystal Substances 0.000 claims abstract description 127
- 238000009826 distribution Methods 0.000 claims abstract description 78
- 230000000694 effects Effects 0.000 claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 239000002245 particle Substances 0.000 claims description 52
- 230000006378 damage Effects 0.000 claims description 35
- 230000001052 transient effect Effects 0.000 claims description 34
- 150000002500 ions Chemical class 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 238000011156 evaluation Methods 0.000 claims description 5
- 239000013307 optical fiber Substances 0.000 claims description 4
- 230000003287 optical effect Effects 0.000 claims description 2
- 238000011282 treatment Methods 0.000 description 46
- 238000012636 positron electron tomography Methods 0.000 description 42
- 210000001519 tissue Anatomy 0.000 description 23
- 206010028980 Neoplasm Diseases 0.000 description 20
- 230000000875 corresponding effect Effects 0.000 description 17
- 230000004044 response Effects 0.000 description 12
- 238000012795 verification Methods 0.000 description 12
- 230000008021 deposition Effects 0.000 description 11
- 230000005855 radiation Effects 0.000 description 11
- 238000002601 radiography Methods 0.000 description 11
- 230000009471 action Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 238000013461 design Methods 0.000 description 8
- 201000011510 cancer Diseases 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 5
- 238000002665 ion therapy Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 238000012879 PET imaging Methods 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 4
- 238000001727 in vivo Methods 0.000 description 4
- 238000010884 ion-beam technique Methods 0.000 description 4
- 238000001959 radiotherapy Methods 0.000 description 4
- 229910004611 CdZnTe Inorganic materials 0.000 description 3
- 241000237970 Conus <genus> Species 0.000 description 3
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000005281 excited state Effects 0.000 description 3
- 230000005251 gamma ray Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- XKUYOJZZLGFZTC-UHFFFAOYSA-K lanthanum(iii) bromide Chemical compound Br[La](Br)Br XKUYOJZZLGFZTC-UHFFFAOYSA-K 0.000 description 3
- 230000033001 locomotion Effects 0.000 description 3
- 238000009206 nuclear medicine Methods 0.000 description 3
- 238000002661 proton therapy Methods 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 229910052710 silicon Inorganic materials 0.000 description 3
- 239000010703 silicon Substances 0.000 description 3
- 238000002560 therapeutic procedure Methods 0.000 description 3
- 241000283070 Equus zebra Species 0.000 description 2
- 241000718541 Tetragastris balsamifera Species 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 239000011163 secondary particle Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000009747 swallowing Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000001225 therapeutic effect Effects 0.000 description 2
- 229910000497 Amalgam Inorganic materials 0.000 description 1
- 208000005623 Carcinogenesis Diseases 0.000 description 1
- 241001637516 Polygonia c-album Species 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- QWUZMTJBRUASOW-UHFFFAOYSA-N cadmium tellanylidenezinc Chemical compound [Zn].[Cd].[Te] QWUZMTJBRUASOW-UHFFFAOYSA-N 0.000 description 1
- CEKJAYFBQARQNG-UHFFFAOYSA-N cadmium zinc Chemical compound [Zn].[Cd] CEKJAYFBQARQNG-UHFFFAOYSA-N 0.000 description 1
- 230000036952 cancer formation Effects 0.000 description 1
- 231100000504 carcinogenesis Toxicity 0.000 description 1
- 238000011281 clinical therapy Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 239000001064 degrader Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- BOYZAERJCXIRAX-UHFFFAOYSA-N lutetium(3+);trisilicate Chemical compound [Lu+3].[Lu+3].[Lu+3].[Lu+3].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-] BOYZAERJCXIRAX-UHFFFAOYSA-N 0.000 description 1
- RCVOTEPLGFCFIG-UHFFFAOYSA-N lutetium(3+);yttrium(3+);silicate Chemical compound [Y+3].[Lu+3].[O-][Si]([O-])([O-])[O-] RCVOTEPLGFCFIG-UHFFFAOYSA-N 0.000 description 1
- 230000036210 malignancy Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000010534 mechanism of action Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000035755 proliferation Effects 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 229910052714 tellurium Inorganic materials 0.000 description 1
- PORWMNRCUJJQNO-UHFFFAOYSA-N tellurium atom Chemical compound [Te] PORWMNRCUJJQNO-UHFFFAOYSA-N 0.000 description 1
- 238000011287 therapeutic dose Methods 0.000 description 1
Images
Landscapes
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明提供一种加马光子侦检装置以及加马光子侦测方法，该装置包括有多个侦测探头以及信号处理电路，每一个侦测探头具有多层闪烁晶体侦检器。利用该多个侦测探头侦测一目标物所产生的加马光子，加马光子先与第一层侦检器作用，后与第二层侦检器作用，然后取得关于至少一第一时间段中每一个侦测探头所产生的电信号，并分别对每一个侦测探头所具有的多个闪烁晶体侦检器所产生的电信号进行第一时讯符合演算，以得到多个第一时讯相关数据。最后，取得关于至少一第二时间段中该每一个侦测探头的第一层闪烁晶体侦检器所产生的电信号，以及进行一第二时讯符合演算，以得到多对第二时讯相关数据。利用该第一和第二时讯相关数据进行活度分布重建。
Description
技术领域
本发明为一种加马光子侦测装置与方法，特别是指一种侦测瞬发加马光子(prompt gamma ray)以及正子互毁加马(positron annihilation gamma ray)的加马光子侦测成像装置与方法。
背景技术
如图1所示，放射治疗(Radiotherapy,RT)就是用高能量的光子(如X光射线)或带电粒子来照射病灶处，以达到“杀死恶性肿瘤”或“抑制恶性肿瘤增生”为目的的一种治疗方式。其中，曲线90代表X光剂量和组织深度关系，由于X光射线与物质作用的物理特性，随射线进入组织深度的增加，相对剂量也随深度逐渐衰减，因此单一X光射束进入体内的路径上有较多的剂量释放在肿瘤前后的正常组织，致使X光射线尚未破坏肿瘤癌细胞前就已经影响正常组织，甚至让肿瘤后的正常组织区域也有剂量沉积。相对于X光治疗的方式，曲线91则为带电粒子，例如，质子或重离子，剂量和深度关系，由于物质作用因质子速度越快行进单位距离损失的能量越少，所以进入病患体表最初释放的动能相当低，直到一定的射程深度，质子或重离子动能大幅下降后产生最大剂量沉积，此称为布拉格峰910(Bragg Peak)，在布拉格峰之后几乎完全没有剂量。
由于质子或重离子治疗具有剂量集中在布拉格峰的特性，因此质子或重离子射程的不确定性影响远较X光更为敏感。若质子或重离子射程预测具有不确定性，其剂量高度集中的布拉格峰所在深度与预期有差异，将造成高剂量区未涵盖到整个肿瘤或高剂量区涵盖到肿瘤周围的重要器官等风险。因此质子或重离子治疗在临床应用上必须要能精准确认质子或重离子射程范围，将其布拉格峰限制在肿瘤区域才能达到最佳的治疗效果。
然质子或重离子治疗过程仍有一些因素影响其不准确性，说明如下：
(1)治疗计划所造成的不确定性：这是由于目前治疗计划主要来自X光CT影像的Hounsfield unit(HU)为组织与电子密度的对应关系，然质子或重离子与生物组织作用主要考虑是组织的原子核及其阻挡本领(stopping power)，因此将CT影像的HU信息转换为生物组织的阻挡本领，将有其不确定性。除此之外，CT影像假影、X光射束硬化(beamhardening)现象等等也是现行评估质子或重离子射程的不确定来源；
(2)实际治疗位置与治疗计划有差异：主要来自病患摆位误差、呼吸或心跳等造成的治疗标的移动等；
(3)治疗过程的肿瘤尺寸改变、病患体型改变等；以及
(4)来自于产生质子或重离子射束的设施所造成的误差，例如降能器(energydegrader)、射束形成(beam delivery)相关装置等误差。
由于上述射程不准确性，临床治疗时通常在欲治疗的临床靶体积(clinicaltarget volume,CTV)外加大治疗安全范围(safety margin)。因射程不确定性而增加的治疗范围，随肿瘤深度y增加(30cm深度肿瘤可增加约14mm治疗安全范围)。在CTV外增加1cm安全范围则有～80cm3的正常健康组织接受与肿瘤相同的辐射剂量，因而增加附近正常组织受损伤及致癌的风险。若CTV附近有辐射耐受度较低的重要器官，也将影响治疗剂量的给予及恶性肿瘤的控制率。为了降低治疗过程的不准确性，提升质子或重离子治疗可信赖度，如果可以在治疗过程中对患者的质子或重离子照射区进行监测，确认质子束在病患体内的路径是否需要修正，这将可缩小治疗安全范围的扩增，提高剂量给予，进而提升肿瘤控制率与治疗效果。因此如何确认与监控质子或重离子束在病患体内的路径是非常重要的研究议题。
现有技术中，有数种发展中的验证质子射程技术，符合实时及非侵入性的技术有质子穿透(proton transmission)造影、瞬发加马光子造影及PET(positron emissiontomography)造影。质子穿透造影需占用额外质子设备时间且会产生额外辐射剂量，也有多库伦散射(multiple Coulomb scattering)效应等问题待克服，短期实际应用于临床的机会不高。瞬发加马光子(PG)造影及PET造影是利用质子与病患组织作用产生的二次粒子—正子发射核及瞬发加马光子做侦测造影。正子互毁加马光子及瞬发加马光子都可穿透病患身体被侦测，可做活体(in vivo)的非侵入性质子射程的量测成像。
质子和重离子与生物组织作用机制雷同，下文以质子治疗论述为代表，不再区分质子和重离子。质子照射时，质子与组织原子核的非弹性碰撞产生加马光子，这些加马光子依诱发产生机制可分成两类。一是当质子与人体组织原子核发生非弹性碰撞，靶核被激发至激态，在纳秒(ns)时间内产生瞬发加马光子(prompt gamma)。另一是靶核被质子撞击碎片化(nuclear fragmentation)而产生正子发射核(positron-emitting nuclei)，如O-15(半衰期:2.037min)、C-11(半衰期:20.385min)、N-13(半衰期:9.965min)等(Newhauserand Zhang,2015)。正子发射核衰变产生的正子(positron)与组织中的电子互毁生成511keV能量的正子互毁加马(positron annihilation gammas)光子对。瞬发加马光子和正子互毁加马光子，它们的能量都可穿透病患身体而被侦测与成像，分析活度分布影像，可获得质子在人体内的路径，故可应用于质子或重离子治疗的射程验证。
正子互毁加马光子对一般利用PET造影来侦测，是目前较成熟且已应用于临床的活体质子射程验证技术。不过正子发射核在病患体内的生物代谢(biological washout)效应、正子发射核产量有限、半衰期短以及治疗***与PET影像的空间对位(co-registration)误差(病患在治疗室与PET造影仪间的移动所造成)等，这些限制了PET对质子射程评估精准度。
至于瞬发加马光子，因生成瞬发加马光子所需质子动能较产生正子发射核所需动能低，所以瞬发加马光子分布的位置和质子沉积最大剂量的位置较为一致，因此利用瞬发加马成像活度评估布拉格峰的准确度较利用PET方法来得高。另外，瞬发加马光子产量较正子互毁加马光子多60～80倍，这表示理论上在相同侦检几何状况，瞬发加马光子造影可以有较强信号，统计误差较PET来得小。现有技术中有利用核子医学领域常见的准直器式加马相机(collimator-based gamma camera)技术，通过针孔改良版的slit准直器或edge-knife准直器结合加马侦检器做瞬发加马射线侦测造影。不过这种技术侦检效率差，因准直器阻挡大部分光子及探头对高能量瞬发加马光子的侦测能力有限。更重要的是因准直器的存在而产生大量的中子(质子诱发的二次辐射粒子之一)，这些中子将干扰侦测探头的信号侦测，造成信杂比(signal-to-noise ratio,S/N)降低，致影响活度影像清晰度及应用其判断质子射程布拉格峰的准确度。此外，因放射治疗诱发的瞬发加马光子能量高，所以准直器设计相当厚重(尺寸达80cm半人高)，占用治疗室有限空间。上述因素导致准直器式加马相机应用在质子临床治疗的实用性受限。
此时另一种于医学影像领域较为少见的康普吞散射成像技术因适用于高能量光子侦测成像，从而成为研究质子或重离子治疗用瞬发加马光子造影的新选择。康普吞散射成像技术主要是应用于天文宇宙射线、核能电厂、国土安全的敏感核种加马光子侦测等。目前国际上有些质子或重离子治疗团队直接采用之前的侦检设计进行康普吞散射成像，例如美国肿瘤治疗中心(Peggy and Charles Stephenson Cancer Center)利用双边硅带型探头(double-sided silicon strip detector,DSSD)侦检器技术、韩国汉阳大学团队利用加马电子锥顶造影(gamma electron vertex imaging,GEVI)侦检技术、日本京都大学团队利用电子追迹康普吞造影(electron tracking Compton camera,ETCC)侦检技术及美国马里兰大学团队的碲锌镉(CdZnTe,CZT)侦检康普吞相机技术。这些团队的侦测探头或是采用游离腔(ionization chamber)技术或是采用半导体侦检器技术。尤其是半导体侦检器其能量分辨率相当高(<1％)，适合应用于康普吞散射成像。但对高能量加马光子而言，其侦检效率差且侦检器信号响应时间长(μs等级)。信号响应时间长也拉长侦测探头的死寂时间(deadtime)，使得侦测探头的单位时间侦检计数效率不彰。这类侦检器可使用于非医学应用，因为非医用造影可容忍数小时甚至数天的长时间数据撷取，以取得足够的侦检计数，但若实施于医用病患身上，长时间造影是不可行的。于是为取得足够的侦检计数，要用相当高的质子剂量照射，瞬发加马成像的计数统计误差才达可被接受程度。
根据美国马里兰大学Dr.J.C.Polf实务经验，即使利用上述最高侦检效率的CdZnTe侦测探头进行质子治疗瞬发加马造影，其所需剂量仍较临床高100倍以上。为弥补侦检器效率不足，也有用上数十层侦检器的设计，然侦检器加上电子***整体尺寸相当可观，这种设计较不利于空间有限的质子治疗室应用。另外，因为侦检器信号响应时间长达微秒(μs)，所以散射侦检器及吸收侦检器的时间分辨率也在微秒等级。微秒级的时间分辨率的造影仪器应用在质子治疗瞬发加马光子监测上，可能导致80％以上侦测计数为相互无关联的非真时讯符合事件(non-true coincident events)，这将会导致瞬发加马光子射源位置误判。由于过去康普吞散射成像探头应用于天文、机场海关的违禁辐射源侦测，因为宇宙射线到达地球的通量有限及违禁辐射源有屏蔽防护，不会有高强度背景干扰，所以微秒级时间分辨率已足以应付非医用的需求。总之，直接应用现有为天文、国土及核能安全需求所设计的DSSD型、CZT型探头对质子治疗的瞬发加马光子造影应用有侦检效率低、时间解析差和非真时讯符合事件比率过高等缺点，这将影响***信杂比、影像清晰度与质子射程布拉格峰估算准确度，所以为质子治疗需求量身射计康普吞散射成像侦测探头有其必要性。
此外，现有技术例如美国公告专利US6,484,051号，该技术公开了一种可以同时侦测同位素所产生的正子互毁加马光子以及瞬发加马光子的装置，由于需要同时量测，因此同时至少需要三个侦检器要工作。另外，瞬发加马光子的信号为正子互毁加马光子的数十倍，同时量测的话，侦测计数率(counting rate)将被低产生率(yield rate)的正子互毁加马光子数量主导，进而影响到侦测信号的信杂比，从而降低影像重建的质量。此外美国公告专利US9,069,089则公开了一种可以同时侦测正子互毁以及瞬发加马光子的技术，在此技术中，在正子互毁加马光子传感器的***设置感测瞬发加马光子的传感器，通过两种不同的传感器来感测加马光子。
发明内容
本发明提供一种加马光子侦测装置与方法，其利用快速闪烁侦检器(例如以LSO、LYSO、LaBr3等闪烁晶体为材料者)，结合PET及康普吞散射侦检技术，制造具有量测正子互毁光子对和瞬发加马光子的二合一(2-in-1)侦检***，再利用时序式数据撷取进行造影及PET嵌入式康普吞成像演算方法进行影像重建和影像融合。二合一侦检***具有量测正子互毁光子对和瞬发加马光子的造影功能，一机双用途，便利性高，PET/PG双功能成像融合造影所得的信息更可提高影像质量与质子/重离子射程验证精准度。
本发明提供的加马光子侦测装置与方法采取共享探头以及时序式数据撷取，其优点除了节省成本之外，更可以提升影像信杂比，其原因在于瞬发加马光子的信号为正子互毁光子对的数十倍，本发明通过单独量测瞬发加马光子信号以及正子互毁加马光子对，可以让这两种光子不相互影响，因此可以兼顾正子互毁加马光子以及瞬发加马光子对于重建活度分布影像的优点。
本发明提供的加马光子侦测装置与方法，可于质子/重离子放射治疗时或核子医学的双同位素(混合正子同位素与加马光子同位素)或特殊同位素(可发射正子与加马光子的同位素(isotope)，如124I,86Y,89Zr等)衰变时，侦测瞬发加马光子以及正子互毁加马光子对。因为上述辐射源其正子发射核与瞬发加马光子分布有相关性。故以PET活度分布为基础经剂量关联处理后可做为康普吞影像重建的先前信息(prior information)。影像重建过程中，利用正子发射核和瞬发加马光子空间分布相关性，以PET嵌入式康普吞成像演算方法进行影像重建，可使得瞬发加马射源所在可能位置分布范围由原本的康普吞角锥面缩小至最大可能机率分布区域，此方法可提高PG影像信杂比并加快康普吞影像重建速度。PET/PG双功能成像可兼顾PET及瞬发加马光子成像优点，可更为提高质子或重离子射程验证及剂量施予的精准度或核医双同位素或特殊同位素造影的质量。
本发明提供一种加马光子侦测装置，其是使用由闪烁晶体所构成的探头。通过闪烁侦检器具有快速时间响应(time response)特性，选用响应快的闪烁晶体(如LSO、LYSO、LaBr3等)搭配快速光电传感器与前置放大电路，***的侦检信号响应时间及时间分辨率可达纳秒等级。本发明通过相同探头达到侦测正子互毁所产生的加马光子对以及因为质子与组织原子核的非弹性碰撞，激发原子核至激态，在纳秒时间产生高能量的瞬发加马光子。本发明的侦测探头相较半导体侦测探头如双边硅带型探头(DSSD)、镉锌碲型侦测探头(CdZnTe,CZT)的微秒响应，纳秒级闪烁侦检***的信号响应时间快、同符时间窗口可缩短数百倍，可提升侦检计数率并大幅度减少非真符合侦测计数的影响，进而提升质子射程影像清晰度与布拉格峰(Bragg peak)估算准确度。
本发明提供一种加马光子侦测方法，利用具有一机双功能设计的加马光子侦测装置，以相同探头达到正子发射核分布和瞬发加马光子的双造影功能。本案采用快速闪烁侦检器为基底组制成像探头，两时段的瞬发加马光子造影(prompt gamma imaging,PGI)模式与正子互毁造影(positron annihilation imaging,PAI)模式数据撷取，使得二合一侦检***具有侦测正子互毁光子对和瞬发加马光子的功能。利用此***，质子/重离子治疗用in-room PET也可以具有瞬发加马光子侦测成像用途。在质子/重离子射束治疗时段(t0～t1)该侦测探头以瞬发加马成像PGI模式进行数据撷取，射束停止后t1～t2时段，探头采正子互毁成像PAI模式进行数据撷取。所获得的原始数据经PET影像重建及PET嵌入式影像重建分别获得PET影像及PG影像。相较先前技术单机单功能及双机双功能，本发明的一机双用途设计，可兼具in-room PET及瞬发加马光子分布成像的功能，不仅经济便利，也因充分利用瞬发加马光子以及正子互毁加马光子空间分布关联性，以PET嵌入式康普吞影像重建融合PET/PG双影像，更提高射程评估精准度。
在一实施例中，本发明提供一种加马光子侦检装置，包括有多个侦测探头以及一信号处理电路。该多个侦测探头，每一个侦测探头具有多层闪烁晶体侦检器，沿着对应的侦测探头的轴向排列，相邻闪烁晶体侦检器具有一距离，每一个侦测探头用以撷取加马光子以产生相应的电信号。该信号处理电路与该多个侦测探头电性连接，用以取得关于一第一时间段中每一个侦测探头所产生的电信号，并对每一个侦测探头所具有的多个闪烁晶体侦检器所产生的电信号进行第一时讯符合演算，以得到对应多个侦测探头的多个第一时讯相关数据，以及取得关于一第二时间段中该每一个侦测探头的第一层闪烁晶体侦检器所产生的电信号，并对进行一第二时讯符合演算，以得到多对第二时讯相关数据。
在另一实施例中，本发明提供一种加马光子侦测方法，其包括有下列步骤，首先提供一加马光子侦检装置，其具有多个侦测探头以及与该多个侦测探头电性连接的一信号处理电路，每一个侦测探头具有多层闪烁晶体侦检器，多层闪烁晶体侦检器沿着对应的侦测探头的轴向排列，相邻闪烁晶体侦检器具有一距离。接着，使该多个侦测探头侦测一目标物所产生的加马光子。接着在一第一时间段中独立的撷取每一个侦测探头所产生的电信号形成多个第一电信号组。然后，以该信号处理电路取得关于至少一第一时间段中每一个侦测探头所产生的电信号，并对每一个侦测探头所具有的多个闪烁晶体侦检器所产生的电信号进行第一时讯符合演算，以得到对应多个侦测探头的多个第一时讯符合相关数据，并进行储存；以及在一第二时间段中撷取每一个侦测探头的第一层闪烁晶体侦检器所产生的电信号，形成多个第二电信号组并进行储存。然后，再以该信号处理电路以及取得关于至少一第二时间段中该每一个侦测探头的第一层闪烁晶体侦检器所产生的电信号，并对进行一第二时讯符合演算，以得到多对第二时讯符合相关数据。
附图说明
图1为带电粒子与传统光子的剂量强度与深度关系曲线图；
图2A为本发明提供的加马光子侦测装置的实施例架构示意图；
图2B至2D为本发明的侦测探头不同实施例示意图；
图2E至图2H为本发明的侦测探头与弧形轨道组合的不同实施状态示意图；
图2I至2K为本发明的侦测探头不同实施例示意图；
图3A为本发明的加马光子侦测方法流程示意图；
图3B为本发明撷取侦测探头电信号与时讯符合演算流程示意图；
图3C为本发明的活度分布重建实施例流程示意图；
图4A至4C为利用图2A的加马光子侦测装置侦测瞬发加马光子以及正子互毁加马光子实施例示意图；
图5为本发明在不同时间段撷取康普吞散射信号以及PET加马光子信号示意图；
图6与图7分别为时讯符合演算不同实施例示意图；
图8A为康普吞圆锥示意图；
图8B则为利用PET影像为基础的剂量分布做为康普吞影像重建先前信息的示意图；
图9为多个探头之间的康普吞圆锥和PET影像为基础的剂量分布交集关联示意图。
附图标记说明：2-加马光子侦测装置；20a,20b,20c,20d,20e-侦测探头；200a,200b,200c,200d,200e,200f,200g,200h-闪烁晶体侦检器；201-闪烁晶体阵列；202-光电传感器；203-读出电路；204-模拟与数字转换电路；205-光纤；21-信号处理电路；22-重建单元；220-PET重建模块；221-康普吞重建模块；222-第一活度分布信息；223-剂量分布信息；224-第二活度分布信息；23-验证单元；24-带电粒子产生装置；25-承载台；26-弧形轨道；3-加马光子侦测方法；30～35-步骤；340～343-步骤；350～352-步骤；8-带电粒子束；80、80’-康普吞圆锥；81-活度分布；82a、82b-区域；9-目标物；90、91-曲线；910布拉格峰；92-反应位置。
具体实施方式
请参阅图2A所示，该图为本发明提供的加马光子侦测装置的实施例架构示意图。在本实施例中，该加马光子侦测装置2包括有多个侦测探头20a与20b、一信号处理电路21、一重建单元22以及一验证单元23。每一个侦测探头20a与20b为具有纳秒反应等级的探头，其侦测面相互对应，且相互平行。每一个侦测探头20a与20b具有多层闪烁晶体侦检器200a～200b与200d～200e，沿着对应的侦测探头的轴向排列，相邻闪烁晶体侦检器200a～200b与200d～200e具有一距离。在一实施例中，每一个侦测探头20a与20b用以撷取沿一行进方向的一带电粒子束所产生的加马光子以产生相应的电信号。每一层闪烁晶体侦检器200a～200b的材料为硅酸镥(LSO)或硅酸钇镥(LYSO)或溴化镧(LaBr3)等快速闪烁晶体，但不以此为限制。该带电粒子束为质子束或者是重离子束，其产生的方式为利用回旋加速器(cyclotron)或同步加速器(synchrotron)产生连续式或脉冲式的带电粒子束，用以照射恶性肿瘤，通过其布拉格峰最大的剂量来消灭肿瘤。此外，在另一实施例中，该侦测探头20a与20b也可以用在特殊同位素(可发射正子与加马光子的同位素，如124I,86Y,89Zr等)衰变时，侦测其所产生的正子互毁加马光子对及加马光子。在一实施例中，该侦测探头20a与20b，为康普吞探头的架构。
在一实施例中，如图2B所示，每一层闪烁晶体侦检器200a～200b分别具有一闪烁晶体阵列201连接至一光电传感器202与具有前置放大电路的读出电路203，使得***的侦检信号响应时间及时间分辨率可达纳秒等级。在一实施例中，闪烁晶体阵列201由多个闪烁晶体201a所构成，该光电传感器202为多阳级的位敏式光电倍增管(position-sensitivephotomultiplier tube,PSPMT)设计，用以将闪烁晶体阵列201所接收的光信号转换成模拟电信号。此外，光电传感器202也可以为光电倍增管(photomultiplier tube,PMT)或者是硅光电倍增元件(silicon photomultiplier,Si-PM)阵列。该读出电路203用以将光电传感器202所输出的电流信号分流、放大并转换为模拟电压信号。
此外，如图2C所示，其为闪烁晶体侦检器另一实施例示意图。闪烁晶体阵列201和光电传感器202之间是通过光纤205传导，闪烁晶体阵列201和光电传感器202之间会有夹角存在。如图2D所示，其为三层侦检器示意图，通过光纤传导可减少以下原因造成的假信号—光电倍增管产生的额外康普吞散射光子、闪烁晶体产生的散射光子在到达下一层侦检器前与光电倍增管发生作用。另外，也可以尽可能地减少侦检层与侦检层之间的距离，提升几何效率。再回到图2A所示，各个侦测探头20a与20b所输出的模拟电信号再经过一积分和模拟转数字电路204，用以将输出的模拟电信号转换成数字电信号。该信号处理电路21，其与该积分和模拟转数字电路204电性连接，用以接收该数字电信号。在一实施例中，该信号处理电路21内有一FPGA(Field Programmable Gate Array)可程序化电路，其可以通过程序对所接收到的电信号进行演算，例如：时讯符合演算，之后将电路21输出数据传输储存于计算机硬盘或内存，以作为后续重建单元22运用。本实施例中，信号处理电路21取得关于一第一时间段中每一个侦测探头所产生的电信号，并对每一个侦测探头所具有的多层闪烁晶体侦检器所产生的电信号进行第一时讯符合演算，以得到对应多个侦测探头的多个第一时讯相关数据，以及取得关于一第二时间段中该每一个侦测探头的第一层闪烁晶体侦检器所产生的电信号，并对进行一第二时讯符合演算，以得到多对第二时讯相关数据。其中，所谓第一时讯符合演算，是指对每一层闪烁晶体侦检器所得到信号做时讯符合演算，用以判断每一个侦测探头的各层闪烁晶体侦检器所侦测到的电信号是否来自于同一加马光子，亦即瞬发加马光子。本实施例中，以侦测探头20a为例，某一个加马光子先与第一层闪烁晶体侦检器200a作用后，产生的散射光子再与第二层闪烁晶体侦检器200b作用，两者有稍许的时间先后差异，因此，第一时讯符合演算，就是指对每一层闪烁晶体侦检器200a,200b或200d,200e所得到信号做时讯符合演算，进而判对通过每一层闪烁晶体侦检器200a,200b的电信号是否是属于同一加马光子的事件，以及每一层闪烁晶体侦检器200d,200e是否属于同一加马光子的事件。而第二时讯符合演算，则是指对正子与电子互毁产生的成对加马光子分别与相对侦检探头的第一层侦检器作用所产生的电信号进行演算，以判断每一个侦测探头的第一层闪烁晶体侦检器所侦测到的电信号是否来自于同一对的加马光子，亦即成对的正子与电子互毁产生的加马光子。例如，本实施例中，成对的加马光子分别与成对的侦测探头20a与20b的第一层闪烁晶体侦检器200a,200d作用所产生的电信号，再经过第二时讯符合运算之后，可以判别200a,200d是否为同一组正子与电子互毁产生的一对加马光子作用所产生的电信号。
前述的侦测探头都是属于顶面照射式(top-on)的侦检器的组合，在另一实施例中，该侦测探头也可以为侧面照射式(side on)与顶面照射式的组合。例如，图2I所示，侦测探头20c具有两个闪烁晶体侦检器200a与200f，其中闪烁晶体侦检器200a为顶面照射式(top-on)的侦检器，而闪烁晶体侦检器200f为侧面照射式的侦检器。此外，如图2J所示，侦测探头20d具有两个闪烁晶体侦检器200f与200g，都是属于侧面照射式的侦检器。另外，要说明的是，侧面照射式的侦检器两端的光电传感器202可以设置在闪烁晶体的两端，如图2I所示，或者是单一端，如图2K所示，其可以根据使用者的需求而定。
重建单元22，在一实施例中，设置在一计算机或具有运算处理能力的***内，其包括有PET重建模块220和康普吞重建模块221，用以从储存于计算机硬盘的储存装置或内存取出第一与第二时讯相关数据。其中，该康普吞重建模块221处理来自PGI模式下，侦测探头20a与20b在质子/重粒子射束治疗时段(t0～t1)感测瞬发加马光子而产生的数据。至于PET重建220模块则处理来自PAI模式下，侦测探头20a与20b在治疗射束停止后t1～t2时段所撷取的成对511keV能量的加马光子产生的数据。PET重建220模块重建所得的第一活度分布信息222进一步进行一反折积(deconvolution)处理，从而得到一剂量分布信息223后，再传给康普吞重建模块221，使得康普吞重建模块221可以将第一剂量分布作为先前信息，而重建获得关于瞬发加马光子的第二活度分布信息224，进而评估出带电粒子在生物体内的实际射程。PET/PG双功能成像可兼顾PET及瞬发加马光子成像优点，可更为提高质子或重粒子射程验证及剂量施予的精准度。在一放射治疗应用实施例中，该重建单元22可与一验证单元23连接，该验证单元23接收重建单元22所输出的活度分布之后，该验证单元23通过PET和PG的活度分布，建立并显示带电粒子在生物体内的实际剂量射程分布，以提供医疗从业人员进行实际射程结果与治疗计划预测射程的比较与确认，进而将该差异作为下一次带电粒子照射参数修正的依据。此外，在一实施例中，该加马光子侦测成像装置2更具有一带电粒子产生装置24，以形成可以治疗以及同步进行活体治疗监测用的加马光子侦测装置，该带电粒子产生装置24用以产生治疗所需的带电粒子。
如图2E所示，在一实施例中，该多个侦测探头20a与20b设置在一弧形的轨道26上，可以在该弧形轨道26上进行位移运动，使得侦测探头20a或20b可以一中心进行转动，以调整其侦测的位置，本实施例以弧形轨道26的中心为转轴来转动。要说明的是，由于本发明具有PET和PG侦测成像的功能，因此在PET造影模式下，各个侦测探头20a与20b必须是相对180度配置，而在PG造影模式下，侦测探头20a或20b方可为任意夹角。在另一实施例中，如图2F所示，当该带电粒子产生装置24产生带电粒子8时，侦测探头20a与20b间具有一夹角，其大于0度，本实施例为90度左右，以增加瞬发加马光子影像三向分辨率的均匀性。此外要说明的是，虽然图2A、2E与2F中的侦测探头为一对，但不以此为限制，其可以根据使用的需求，例如：分辨率、成本等考虑，设置所需要的数量，例如：在图2G中，具有两组侦测模块，每一个侦测模块具有多个侦测探头20a,20c,20e以及20b,20d,20f。此外，如图2H所示，每一个侦测探头20a与20b也可以自身为轴进行转动及朝向或远离圆心的平移移动。
请参阅图3A所示，该图为本发明的加马光子侦测方法流程示意图。在本实施例中，该方法3包括有下列步骤：首先进行步骤30，根据一光子检测决定一目标物内所具有的一反应位置。在本步骤中，该光子检测的一实施例为了进行治疗计划所实施的计算机断层扫描(computed tomography,CT)，其主要通过X光对生物体，例如：人体，进行扫描，该人体即为该目标物，该反应位置，在一实施例为肿瘤。通过计算机断层决定出肿瘤的位置。找到了反应位置之后，进行步骤31，根据生物体的状况，例如：周围的软组织分布，以及该反应位置信息，转换成带电粒子束的一照射参数，其包括有照射的能量、剂量以及行进方向与深度。因为反应位置在目标物内是属于三维度的空间分布，因此随着带电粒子束投射至目标物的行进深度不同，其所需的能量也会有所差异。要说明的是，带电粒子束的行进方向与深度，并不以单一个为限制，可以根据治疗的需求而有多个行进方向与深度，亦即，从多的方向与多个射程深度照射至该反应位置上。
本实施例中的带电粒子为质子，由于质子进入病患体表最初释放的剂量相当低，直到一定的射程深度，质子动能大幅下降后产生最大剂量沉积，而这个射程深度所对应的释放剂量称为布拉格峰，在布拉格峰之后几乎完全没有剂量(如图1所示)。而步骤31所根据反应位置决定的能量或剂量，即是利用此特性，通过质子能量与质子射程的调控将剂量最高的布拉格峰落在目标位置区，使反应位置后方的正常组织则可避开辐射照射风险。
步骤31之后，接着进行步骤32，在本步骤中，提供如图4所示的加马光子侦测装置2。为了确保步骤31所决定的剂量或能量可以集中在布拉格峰所对应的位置，因此在治疗的过程中，可以通过本发明提供的加马光子侦测装置2来进行实时性的同步活体监控，以确保治疗的效率以及降低误差，进而可以保护正常组织免于辐射的伤害。如图4A所示，目标物9为一生物体，平躺在承载台25上，目标物9内具有反应位置92，本实施例为肿瘤。多个侦测探头20则环设在目标物9的***，本实施例虽有两个侦测探头20a与20b，但其数量可以根据需求而定。之后再进行步骤33，使该带电粒子产生装置24产生一带电粒子束8。当治疗过程中，根据步骤31所产生的照射信息，投射带电粒子束8朝向目标物9投射之后，进入到目标物内，质子与组织原子核的非弹性碰撞后，主要会产生两种加马光子，一种是质子与病患组织作用产生的二次粒子—正子发射核衰变产生的正子与组织中的电子互毁生成一对511keV能量的加马光子。另一种是当质子与人体组织原子核发生非弹性碰撞，靶核被激发至激态，在纳秒时间内产生瞬发加马光子。
为了侦测这两种加马光子，通过步骤34在一期间以第一时间段与第二时间段来侦测撷取电信号。以产生带电粒子的本实施而言，该第一时间段代表带电粒子启动的时间段，该第二时间段则代表该带电粒子关闭的时间段。如图4A所示，在带电粒子束8启动的第一时间段里，每一个侦测探头20a与20b的各层闪烁晶体侦检器200a与200b以及200d与200e，会侦测到瞬发加马光子入射的信号，以侦测探头20a为例，亦即某一个加马光子先与第一层烁晶体侦检器200a作用后，产生的散射光子再与第二层烁晶体侦检器200b作用，两者有稍许的时间先后差异，同理，侦测探头20b也是如此。反之，如图4B所示，在带电粒子束8关闭的第二时间段里，每一个侦测探头20a与20b的第一层闪烁晶体侦检器200a与200d，会侦测到因为正子互毁所产生的成对加马光子入射的信号，亦即，正子与电子互毁所产生的成对加马光子会分别与相对侦检探头20a与20b的第一层闪烁晶体侦检器200a与200d作用。要分段侦测的主因是因为瞬发加马光子只会在带电粒子启动的时期产生，而且存在时间相当短暂，反之正子互毁加马光子则反应时间长，因此两种加马光子的特性不同。另外，量测正子互毁加马光子对以PET活度影像估算带电粒子射程已具临床应用成熟度，分开量测瞬发加马光子信号以及正子互毁加马光子对，可以通过较为成熟的PET技术提升瞬发加马光子影像重建的准确性，且瞬发加马光子的信号为正子互毁加马光子的数十倍，瞬发加马光子造影可提升影像的信杂比，进而强化带电粒子射程验证应用的精准度。要说明的是，虽然图4B的正子对是由第一层闪烁晶体侦检器200a与200d所感测，在另一实施例中，如图4C所示，根据该正子对亦可以通过第二层的闪烁晶体侦检器200b与200e来进行感测。
请参阅图3B所示，在步骤34中，更包括有步骤340在第一时间段中独立的撷取每一个侦测探头20a与20b所产生的电信号，形成两组对应侦测探头20a与20b的第一电信号组。在本步骤中，以侦测探头20a为例，主要是使该信号处理电路21分别记录相应的加马光子E0与第一层闪烁晶体侦检器200a发生康普吞散射(Compton scattering)作用的第一位置、第一时间与第一沉积能量ΔE以及散射加马光子能量E1与第二层闪烁晶体侦检器200b作用的第二位置、第二沉积能量ΔE2与第二时间，以构成前述的第一电信号组。理想上第二层闪烁晶体侦检器200b完全吸收散射加马光子能量E1，所以第二沉积能量ΔE2＝E1。之后，再以步骤341取得步骤340中的关于该第一时间段中每一个侦测探头所产生的多个第一电信号组，并对每一个侦测探头所具有的多个闪烁晶体侦检器所产生的电信号进行第一时讯符合演算，过滤非真的无关联信号，以得到对应多个侦测探头的多个瞬发加马光子第一时讯相关数据。以图6来说明步骤340，其为第一时讯符合演算示意图。以侦测探头20a为例，将其所具有的第一至第二层闪烁晶体侦检器200a～200b，分别侦测的事件所对应的信号取出，图6所示为三个事件I1,I2,I3。每一个事件具有两个电信号，分别对应每一层的闪烁晶体侦检器200a～200b，每一组电信号含有位置、时间与能量的信息。再通过比较第一与第二时间信号，如落在设定的符合时间窗口内，则表示该第二层闪烁晶体侦检器测量所得电信号来自与第一层闪烁晶体侦检器作用的加马光子的机会大。在本实施例中，有三个相同时间窗口宽度的时段W1,W2以及W3，要说明的是，时间窗口的大小可以根据需求来决定时间窗口的长度。从图6可以看出，W2的时段内，只有闪烁晶体侦检器200a有信号，闪烁晶体侦检器200b则无，而时段W1以及W3内则闪烁晶体侦检器200a与200b同时有信号，因此代表W1和W3时段撷取的信号属于时讯符合的信号，亦即认定是否为同一组入射瞬发加马光子及其散射加马光子所产生的电信号。如是，则予以储存记录，如否则不予以储存记录。
接着进行步骤342，亦即在第二时间段中撷取每一个侦测探头20a与20b的第一层闪烁晶体侦检器200a与200d所产生的电信号，形成多个第二电信号组。在本步骤342中，主要是使该信号处理电路21处理并数字化加马光子在侦测探头20a所得代表第一层闪烁晶体侦检器200a作用的第一位置、第一时间与第一沉积能量以及加马光子在侦测探头20b的第一层闪烁晶体侦检器200d作用的第二位置、第二沉积能量与第二时间，以构成前述的第二电信号组。如图5所示，其中侦测期间从t0～t2，在该侦测期间的t0～t1的第一时间段内，步骤340中会获得关于侦测探头20a的多个电信号组A1以及侦测探头20b的多个电信号组B1，同理在t1～t2的第二时间段中，执行步骤341则会获得关于各侦测探头20a与20b中所具有的第一层闪烁晶体侦检器200a与200d所得的多组电信号C1。根据现有的核医PET闪烁侦检器设计，由侦检器侦测电路信号，则可以得知加马光子是与哪一个闪烁晶体作用，进而得知作用位置与发生时间，再将侦测到信号的时间记录下来，沉积在侦检器的能量同样可以由相应的电信号强度得知。因此，每一个信号Pa1～Pan,Pb1～Pbn,Pd1～Pdn与Pe1～Pen含有位置、时间以及沉积能量的信息。
之后，再进行步骤343取得步骤342中关于该第二时间段中该每一个侦测探头的第一层闪烁晶体侦检器所产生的多个第二电信号组，并对进行一第二时讯符合演算，过滤非真的无关联信号，以得到多对第二时讯相关数据。在本步骤中，时讯符合演算的对象是每一个侦测探头的第一层闪烁晶体侦检器，以图4A为例，即是侦测探头20a与20b的闪烁晶体侦检器200a和200d。请参阅图7所示，在图7中显示的是闪烁晶体侦检器200a和200d所产生对应有三个事件I4～I6的多个第二电信号组，本实施例中，通过三个相同时间窗口宽度的时段W4～W6以及能量来确认这些信号是否属于同一个正子互毁所产生的一对511keV的加马光子所产生的一对电信号。如是，则予以储存记录，如否则不予以储存记录。要说明的是，在另一实施例中，步骤340～341以及步骤342～343可以独立同步进行。
再回到图3A所示，步骤34之后，可以进一步利用图4A所示的重建单元22来进行步骤35活度分布重建，进而将最后活度分布的结果转换成影像。本步骤重建的基础在于利用PET活度分布做为康普吞成像的先前信息，进而再进行康普吞影像重建。因为病患体内的质子或重离子射程途径除产生正子发射核外，也产生瞬发加马光子，所以正子发射核与质子或重离子剂量分布有关，也与瞬发加马光子分布具相关性。故以PET活度分布为基础经剂量关联处理后可做为康普吞影像重建的先前信息。如图3C所示，首先进行步骤350，先将步骤343所得到的第二时讯符合数据先进行前置处理，包括有侦检器校正、侦检器输出的位置信号与晶***置产生关联，转换成位置与能量数据进行PET重建，得到关于正子互毁加马光子的第一活度分布。正子互毁加马光子的活度分布配合病患体内的质子或重离子射程途径在空间中大致呈直线束的分布状态，其代表着根据正子互毁加马光子所得到的带电粒子束的入射射程状况。接着进行步骤351将步骤350中的正子互毁加马光子的第一活度分布再转换为由PET估算获得的剂量分布信息。最后进行步骤352将步骤351所得到的正子互毁加马光子的剂量分布信息进行关联，也就是以步骤351的剂量分布信息作为步骤352的先前信息，进行PET嵌入式康普吞影像重建，亦即将康普吞圆锥面与前述步骤351由正子互毁加马光子的活度分布转换估算获得的剂量分布
进行交集，以进一步得到一第二活度分布。
以下说明采用步骤352的原理，如图8A与图8B所示，图8A为单纯将代表瞬发加马光子与侦检器的作用位置与沉积能量的电信号，例如，若以步骤343所得的单一事件的电信号，例如：图6中的Pa1与Pb1，进行前置处理，例如：包括有侦检器校正、侦检器输出的位置信号与晶***置产生关联，之后再进行重建以得到活度分布图。根据康普吞散射原理，可以从入射闪烁侦检器的加马光子以及散射加马光子之间所具有的散射角来估算出产生加马光子的可能起源位置。这中间的演算，可以利用时讯符合演算来进行入射和散射加马光子的配对。从第一层以及第二层闪烁晶体侦检器通过步骤341时讯符合演算处理而被相互关连在一起的加马光子，根据侦测到的能量以及方程式(1)即可以得知产生加马光子起源的可能区域。
where E0＝E1+ΔE1
其中，原始能量为E0的高能加马光子与第一层闪烁侦检器200a及200d产生康普吞散射作用，释放部分能量ΔE1于第一层闪烁侦检器后为E1能量的康普吞散射光子。对被量测到单一的瞬发加马射线而言，此起源点的可能位置分布在一康普吞圆锥(Compton cone)80表面。然而，由于康普吞圆锥80，只能知道作用位置(顶点)与散射角θ，因此只要是在锥面上的任何一个位置都有可能是目标物所在的位置，因此单纯只有利用康普吞圆锥来推断目标所在的位置，其不确定性仍高不利于治疗射程的监控。为了解决这个问题，采用了步骤351结合步骤352的方式，亦即利用质子或重离子在病患体***程途径产生二次粒子—正子发射核和瞬发加马光子，所以正子发射核与瞬发加马光子分布均与质子/重离子射程与剂量有关联的特性，以PET活度分布为基础经剂量关联处理后转换成剂量分布，来做为康普吞成像的先前信息来重建活度分布。在此实施例，经过步骤351之后PET影像所构成的正子发射核活度分布
经转换为剂量分布
转换的方式可以采用例如滤波函数反折积方法转换正子发射核活度分布为剂量分布(S.Remmele et.al.2011Phys.Med.Biol.)，但不以此为限。如图8B所示，其中标号81的直线束代表由PET影像所估算出的剂量空间分布
其和康普吞圆锥表面80交集构成了一两区域82a与82b。于是最有可能的瞬发加马光子起源位置由康普吞圆锥80表面缩小至区域82a与82b，此时最有可能的瞬发加马光子起源位置
可以下式表示，
于是在步骤352，经过交集之后，单一瞬发加马光子起源可能位置分布范围可由原本的康普吞圆锥80表面缩小至区域82a与82b。在PET嵌入式康普吞影像重建过程中，射源可能位置分布范围缩小将可大大提高重建后影像的信杂比质量，这不仅可加速瞬发加马光子影像重建收敛速度，更重要的是提升射程评估的精准度。
在图8B中为步骤343的第二时讯相关数据经步骤350～351的PET影像重建及剂量转换后的分布与步骤341的第一时讯相关数据经步骤352所形成的活度分布所交集而成的结果。由于图4A的实施例有两个侦测探头20a与20b，因此对于各个侦测探头20a与20b在步骤353之后形成如图9所示的活度分布交集。通过侦测探头20a与侦测探头20b的康普吞圆锥80和80’以及侦测探头20a和侦测探头20b的第一层闪烁晶体侦检器200a与200d所构成的PET架构，其所得的活度分布及剂量转换后的分布81交集可以将加马光子可能位置缩小至交集的区域，进而提升后续射程评估精准度。
以肿瘤治疗为例，因为在治疗之前会进行治疗计划评估，评估肿瘤的位置、带电粒子束应照射的角度及施予的剂量等。为了确认治疗过程中，带电粒子束在生物体内的实际射程与产生布拉格峰的位置与治疗计划预估射程的差异，因此，本发明的装置以及方法可以在每次治疗实时得到带电粒子束的实际射程，通过验证单元23提供显示的信息，可以让用户验证比较治疗计划评估的射程与由步骤352产生的活度分布所决定出的实际射程的差异，进而做为后续治疗调整该带电粒子束照射参数的依据，进而可以达到活体治疗监控与调整建议的双重效果，提升治疗的有效性并降低因射程不确定造成伤害健康组织的副作用。
前述的实施例为具有两层闪烁晶体的实施例，在另一实施例中，可以具有三层的闪烁晶体，如图2D所示的结构，其中设置第三层闪烁晶体侦检器的目的在于考虑康普吞散射成像的两层闪烁晶体侦检器侦检器设计过于理想，因为第二层侦检器(亦称为吸收侦检器)必须补获散射光子的所有能量，才能获得方程式(1)或方程式(2)所需的加马光子初始能量E0。然因质子治疗诱发的瞬发加马光子能量分布范围广且能量高，一般能量在数个keV到10MeV，甚至可达15MeV。因此，实际侦测时的初始加马光子能量未知，而散射加马光子能量E1全部被吸收侦检器阻挡的机率有限，想通过E1及沉积于第一层侦检器的能量ΔE1的总和获知初始加马光子全能量E0的机会不高。为了强化这个部分，在本实施例中，将侦测探头改用三层的闪烁侦检器设计，第三层的侦检器仅需量测到第二个散射光子的作用位置但无需记录第三层侦检器的沉积能量。θ2的计算可通过三个闪烁晶体侦检器各自记录到的入射加马光子与散射加马光子位置获得，之后将所有已知能量记录代入即可计算出入射光子源可能位置(亦即为质子与组织作用产生的瞬发加马光子的起始位置)所在的康普吞圆锥角θ1，如方程式(3)所示。至于步骤361的PET嵌入式影像重建，在此三层闪烁侦检器状况下，最有可能的瞬发加马光子起源位置
仍是通过方程式(2)计算。
E0＝E1+ΔE1
E1＝E2+ΔE2
此外，要说明的是，前述的实施例中，为侦测带电粒子束所产生的瞬发加马光子以及成对的正子互毁加马光子。在另一实施例中，本发明的侦测装置也可以应用在核子医学的双同位素(混合正子同位素与加马光子同位素)或特殊同位素(可发射正子与加马光子的同位素，如124I,86Y,89Zr等)的造影等。因此，并不以侦测带电粒子束所产生的加马光子为限制。在本实施例中，同样地，第一时间段侦测加马光子，第二时间段则侦测正子互毁的加马光子，如此第一与第二时间段交替进行，取得的电信号在根据图3B与图3C的步骤进行处理，进而重建出对应的活度分布与影像。
以上所述仅记载本发明为呈现解决问题所采用的技术手段的较佳实施方式或实施例而已，并非用来限定本发明专利实施的范围。即凡与本发明专利申请范围文义相符，或依本发明专利范围所做的均等变化与修饰，皆为本发明的保护范围所涵盖。
Claims (19)
1.一种加马光子侦检装置，其特征在于，包括：
多个侦测探头，每一个侦测探头具有多层闪烁晶体侦检器，每一层闪烁晶体侦检器沿着对应的侦测探头的轴向排列，相邻层闪烁晶体侦检器具有一距离，每一个侦测探头用以撷取加马光子以产生相应的电信号，其中每一个侦测探头的各层闪烁晶体侦检器，于一第一时间段里侦测到瞬发加马光子入射的所产生的电讯号，而每一个侦测探头的第一层闪烁晶体侦检器，于一第二时间段里侦测到成对的正子互毁加马光子入射的讯号；
一信号处理电路，与该多个侦测探头电性连接，用以取得关于该第一时间段中每一个侦测探头所产生的电信号，并对每一个侦测探头所具有的多层闪烁晶体侦检器所产生的电信号进行第一时讯符合演算，以判断每一个侦测探头的各层闪烁晶体侦检器所侦测到的电信号是否来自于同一瞬发加马光子，以得到对应多个侦测探头的多个第一时讯相关数据，以及取得关于一第二时间段中该每一个侦测探头的第一层闪烁晶体侦检器所产生的电信号，并进行第二时讯符合演算，以得到多对第二时讯相关数据，以判断每一个侦测探头的第一层闪烁晶体侦检器所侦测到的电信号是否来自于同一对的正子互毁加马光子；以及
一重建单元，包括有PET重建模块和康普吞重建模块，PET重建模块用以根据该多对第二时讯相关数据重建一第一活度分布，再将该第一活度分布进行一转换之后得到剂量分布，再传给康普吞重建模块，使得康普吞重建模块将剂量分布作为先前信息，将该剂量分布与该多个第一时讯相关数据相关联以重建一第二活度分布，进而评估出带电粒子在生物体内的实际射程。
2.如权利要求1所述的加马光子侦检装置，其特征在于，该瞬发或正子互毁加马光子为至少一种同位素经由一衰变作用所产生。
3.如权利要求1所述的加马光子侦检装置，其特征在于，更具有一带电粒子产生装置，用以产生一带电粒子束，该带电粒子束沿一行进方向与一目标物作用时产生该瞬发或正子互毁加马光子。
4.如权利要求3所述的加马光子侦检装置，其特征在于，该第一时间段为该带电粒子启动的时间段，该第二时间段为该带电粒子关闭的时间段。
5.如权利要求3所述的加马光子侦检装置，其特征在于，该带电粒子为质子或重离子。
6.如权利要求1所述的加马光子侦检装置，其特征在于，该多个探头能够独立地以一中心位置进行转动或相对移动，以改变侦测位置。
7.如权利要求1所述的加马光子侦检装置，其特征在于，每一层闪烁晶体侦检器更包括有：
一闪烁晶体阵列；
一光电传感器，与该闪烁晶体耦接，用以将闪烁晶体阵列所产生的光信号转换成电信号；以及
一读出电路，包含前置放大电路，其与该光电传感器电性连接。
8.如权利要求7所述的加马光子侦检装置，其特征在于，该闪烁晶体阵列通过光纤与该光电传感器耦接。
9.如权利要求7所述的加马光子侦检装置，其特征在于，该闪烁晶体阵列为侧面照射式闪烁晶体阵列或者是正面照射式闪烁晶体阵列。
10.一种加马光子侦测方法，其特征在于，包括下列步骤：
提供一加马光子侦检装置，其具有多个侦测探头以及与该多个侦测探头电性连接的一信号处理电路，每一个侦测探头具有多层闪烁晶体侦检器，每一层闪烁晶体侦检器沿着对应的侦测探头的轴向排列，相邻层闪烁晶体侦检器具有一距离；
该多个侦测探头侦测一目标物所产生的瞬发或正子互毁加马光子；
在一期间撷取侦测探头的电信号，其中该期间由一第一时间段以及一第二时间段组成，其中每一个侦测探头的各层闪烁晶体侦检器，于该第一时间段里侦测到瞬发加马光子入射的所产生的电讯号，而每一个侦测探头的第一层闪烁晶体侦检器，于该第二时间段里侦测到成对的正子互毁加马光子入射的讯号；
取得关于该第一时间段中每一个侦测探头所产生的电信号，并对每一个侦测探头所具有的多层闪烁晶体侦检器所产生的电信号进行第一时讯符合演算，以得到对应多个侦测探头的多个第一时讯相关数据；
取得关于该第二时间段中该每一个侦测探头的第一层闪烁晶体侦检器所产生的电信号，并进行一第二时讯符合演算，以得到多对第二时讯相关数据；以及
提供重建单元，包括有PET重建模块和康普吞重建模块，其中PET重建模块用以根据该多对第二时讯相关数据重建一第一活度分布，再将该第一活度分布进行一转换之后得到剂量分布，再传给康普吞重建模块，使得康普吞重建模块将剂量分布作为先前信息，将该剂量分布与该多个第一时讯相关数据相关联以重建一第二活度分布，进而评估出带电粒子在生物体内的实际射程。
11.如权利要求10所述的加马光子侦测方法，其特征在于，该瞬发或正子互毁加马光子为至少一种同位素经由一衰变作用所产生。
12.如权利要求10所述的加马光子侦测方法，其特征在于，该瞬发或正子互毁加马光子为一带电粒子束沿一行进方向与一目标物作用所产生，该带电粒子束的一照射参数对应一最大射程评估位置。
13.如权利要求12所述的加马光子侦测方法，其特征在于，该第一时间段为该带电粒子启动的时段，该第二时间段为该带电粒子关闭的时段。
14.如权利要求12所述的加马光子侦测方法，其特征在于，该带电粒子为质子或重离子。
15.如权利要求10所述的加马光子侦测方法，其特征在于，于该第一时间段撷取电信号更包括有下列步骤：
每一侦测探头中的多层闪烁晶体侦检器分别侦测关于瞬发加马光子及其散射加马光子的两个事件；以及
对同一侦测探头中的多层闪烁晶体侦检器所产生的电信号进行处理，以得到关于该多层闪烁晶体侦检器的多个第一数字电信号。
16.如权利要求15所述的加马光子侦测方法，其特征在于，于该第一时讯符合演算更包括有下列步骤：
取得关于该多层闪烁晶体侦检器的多个第一数字电信号；以及
以一时间窗口依序比较每一个侦测探头中的多层闪烁晶体侦检器所分别产生的多个第一数字电信号，如果不同层闪烁晶体侦检器的第一数位信号的时间差落在设定的符合时间窗口内，则表示该不同层闪烁晶体侦检器测量所得第一数字信号来自于第一层闪烁晶体侦检器作用的瞬发加马光子。
17.如权利要求10所述的加马光子侦测方法，其特征在于，该第二时间段撷取电信号更包括有下列步骤：
每一侦测探头中的一第一层闪烁晶体侦检器分别侦测正子互毁加马光子事件；以及
对每一侦测探头中的第一层闪烁晶体侦检器所产生的电信号进行处理，以得到关于每一第一层闪烁晶体侦检器的多个第二数字电信号。
18.如权利要求17所述的加马光子侦测方法，其特征在于，于该第二时讯符合演算更包括有下列步骤：
取得关于不同侦测探头的第一层闪烁晶体侦检器的多个第二数字电信号；以及
以一时间窗口依序比较每一个侦测探头中的第一层闪烁晶体侦检器所分别产生的多个第二数字电信号，如果不同的侦测探头中的第一层闪烁晶体侦检器所分别产生的第二数字电信号时间差落在设定的符合时间窗口内且在该时间窗口内的第二数字电信号其能量在预设的能量范围内，则表示每一第一层闪烁晶体侦检器测量所得第二数字电信号来自同一对正子互毁加马光子。
19.如权利要求10所述的加马光子侦测方法，其特征在于，每一个闪烁晶体侦检器更包括有：
一闪烁晶体阵列；
一光电传感器，与该闪烁晶体耦接，用以将闪烁晶体阵列所产生的光讯转换成电信号；以及
一读出电路，包含前置放大电路，该前置放大电路与该光电传感器电性连接，其中该闪烁晶体阵列为侧面照射式闪烁晶体阵列或者是正面照射式闪烁晶体阵列。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW105120636A TWI599790B (zh) | 2016-06-29 | 2016-06-29 | 加馬光子偵測成像裝置與方法 |
| TW105120636 | 2016-06-29 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107544086A CN107544086A (zh) | 2018-01-05 |
| CN107544086B true CN107544086B (zh) | 2020-03-10 |
Family
ID=60719715
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710448691.5A Expired - Fee Related CN107544086B (zh) | 2016-06-29 | 2017-06-14 | 加马光子侦测成像装置与方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN107544086B (zh) |
| TW (1) | TWI599790B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102018201601A1 (de) * | 2018-02-02 | 2019-08-08 | Siemens Healthcare Gmbh | Detektormodul |
| CN108957517A (zh) * | 2018-06-15 | 2018-12-07 | 中派科技（深圳）有限责任公司 | 用于正电子发射成像设备的探测器及正电子发射成像设备 |
| US10908086B2 (en) | 2018-12-26 | 2021-02-02 | Industrial Technology Research Institute | Signal processing system and method thereof |
| CN110477942B (zh) * | 2019-08-20 | 2023-06-23 | 上海联影医疗科技股份有限公司 | 一种pet探测器以及医学影像设备 |
| CN114930195A (zh) | 2020-01-10 | 2022-08-19 | 特拉佩特有限公司 | 离子束发射装置及其检测*** |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6346706B1 (en) * | 1999-06-24 | 2002-02-12 | The Regents Of The University Of Michigan | High resolution photon detector |
| CN102246057A (zh) * | 2008-12-10 | 2011-11-16 | 皇家飞利浦电子股份有限公司 | 用作可缩放pet和spect***构建块的自主探测器模块 |
| CN102293661A (zh) * | 2010-06-03 | 2011-12-28 | 株式会社东芝 | 核医学成像装置以及放射线治疗装置 |
| CN104066479A (zh) * | 2011-09-29 | 2014-09-24 | 约翰霍普金斯大学 | 使用异质性补偿迭加进行用于放射疗法的剂量计算 |
| CN105105781A (zh) * | 2015-09-17 | 2015-12-02 | 武汉瑞福宁科技有限公司 | 一种平板pet***及用途 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012135725A2 (en) * | 2011-03-30 | 2012-10-04 | The Board Of Trustees Of The Leland Stanford Junior University | Methods and systems for increasing the sensitivity of simultaneous multi-isotope positron emission tomography |
| WO2014209972A1 (en) * | 2013-06-24 | 2014-12-31 | Massachusetts Institute Of Technology | Normalization correction for multiple-detection enhanced emission tomography |
-
2016
- 2016-06-29 TW TW105120636A patent/TWI599790B/zh not_active IP Right Cessation
-
2017
- 2017-06-14 CN CN201710448691.5A patent/CN107544086B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6346706B1 (en) * | 1999-06-24 | 2002-02-12 | The Regents Of The University Of Michigan | High resolution photon detector |
| CN102246057A (zh) * | 2008-12-10 | 2011-11-16 | 皇家飞利浦电子股份有限公司 | 用作可缩放pet和spect***构建块的自主探测器模块 |
| CN102293661A (zh) * | 2010-06-03 | 2011-12-28 | 株式会社东芝 | 核医学成像装置以及放射线治疗装置 |
| CN104066479A (zh) * | 2011-09-29 | 2014-09-24 | 约翰霍普金斯大学 | 使用异质性补偿迭加进行用于放射疗法的剂量计算 |
| CN105105781A (zh) * | 2015-09-17 | 2015-12-02 | 武汉瑞福宁科技有限公司 | 一种平板pet***及用途 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107544086A (zh) | 2018-01-05 |
| TW201800776A (zh) | 2018-01-01 |
| TWI599790B (zh) | 2017-09-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107544086B (zh) | 加马光子侦测成像装置与方法 | |
| Pausch et al. | Detection systems for range monitoring in proton therapy: Needs and challenges | |
| US7218701B2 (en) | Radiological imaging apparatus | |
| EP2347285B1 (en) | Device for detecting highly energetic photons | |
| JP5339551B2 (ja) | 遮蔽型放射線治療・画像化複合装置、及び、その制御プログラム | |
| US10942282B2 (en) | Combined imaging detector for x-ray and nuclear imaging | |
| JP5339562B2 (ja) | 核医学イメージング装置の画像化方法、システム、核医学イメージグシステム及び放射線治療制御システム | |
| JP6842694B2 (ja) | 部分リングｐｅｔ装置及びｐｅｔ装置 | |
| JP2014048267A (ja) | Ｐｅｔ装置における同時計数判定方法及び装置 | |
| KR101241821B1 (ko) | 컴프턴 현상을 이용한 양전자 단층 촬영 장치, 방법, 및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 매체 | |
| Tian et al. | Radiopharmaceutical imaging based on 3D-CZT Compton camera with 3D-printed mouse phantom | |
| EP3617750B1 (en) | Theranostic imaging | |
| JP2022553079A (ja) | 患者の身体部分内の少なくとも１つのトレーサーの微細な位置の向上された決定のための装置及びそれぞれの方法 | |
| KR102449932B1 (ko) | 컴프턴 효과를 활용한 방사선 민감도 증진 방법 및 방사선 민감도 증진 시스템 | |
| Parodi | In vivo treatment verification | |
| Lecoq | Detectors in Medicine and Biology: Applications of Detectors in Technology, Medicine and Other Fields | |
| Torres-Espallardo et al. | Simulation study of resistive-plate-chambers based pet for hadron-therapy monitoring | |
| Mirabelli | MONDO: A neutron tracker for particle therapy secondary emission fluxes measurements | |
| Dauvergne et al. | Prompt-gamma monitoring of proton-and carbon-therapy. Combined development of Time-of-Flight collimated-and Compton-cameras | |
| Gasparini et al. | MONDO: A neutron tracker for particle therapy secondary emission measurements | |
| Necchi | Positron Emission Tomography: status of the art and future perspectives | |
| Bossis et al. | A portable gamma camera for the optimization of the patient dosimetry in radioiodine therapy of thyroid diseases |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200310 |