CN106853272A - 射束的照射角度评价方法 - Google Patents

Abstract

一种射束的照射角度评价方法，包括：对射束的照射角度进行取样的步骤，其中，所述射束的照射角度定义为射束的照射点至肿瘤预设点的向量方向；计算射束通过器官径迹的步骤；判断肿瘤是否充分涵盖于有效治疗深度内：如果是，则进入计算评价因数、记录照射条件及计算结果的步骤并返回对射束的照射角度进行取样的步骤；如果否，则进入给予最差评价因数的步骤并返回对射束的照射角度进行取样的步骤。根据此射束的照射角度评价方法，可以清楚地认识到射束在某一特定位置和某一特定角度照射时的优劣表现，从而为医生或物理师决定照射模式提供强有力的数据支持。

Description

射束的照射角度评价方法

技术领域

本发明涉及一种照射角度评价方法，尤其涉及一种射束的照射角度评价方法。

背景技术

随着原子科学的发展，例如钴六十、直线加速器、电子射束等放射线治疗已成为癌症治疗的主要手段之一。然而传统光子或电子治疗受到放射线本身物理条件的限制，在杀死肿瘤细胞的同时，也会对射束途径上大量的正常组织造成伤害；另外由于肿瘤细胞对放射线敏感程度的不同，传统放射治疗对于较具抗辐射性的恶性肿瘤(如：多行性胶质母细胞瘤(glioblastoma multiforme)、黑色素细胞瘤(melanoma))的治疗成效往往不佳。

为了减少肿瘤周边正常组织的辐射伤害，化学治疗(chemotherapy)中的标靶治疗概念便被应用于放射线治疗中；而针对高抗辐射性的肿瘤细胞，目前也积极发展具有高相对生物效应(relative biological effectiveness,RBE)的辐射源，如质子治疗、重粒子治疗、中子捕获治疗等。其中，中子捕获治疗便是结合上述两种概念，如硼中子捕获治疗(BoronNeutron Capture Therapy,BNCT)，借由含硼药物在肿瘤细胞的特异性集聚，配合精准的射束调控，提供比传统放射线更好的癌症治疗选择。

硼中子捕获治疗是利用含硼(¹⁰B)药物对热中子具有高捕获截面的特性，借由¹⁰B(n,α)⁷Li中子捕获及核分裂反应产生⁴He和⁷Li两个重荷电粒子，两粒子的总射程约相当于一个细胞大小，因此对于生物体造成的辐射伤害能局限在细胞层级，当含硼药物选择性地聚集在肿瘤细胞中，搭配适当的中子射源，便能在不对正常组织造成太大伤害的前提下，达到局部杀死肿瘤细胞的目的。

现有中子捕获治疗计划系统中，其照射几何角度皆为人工根据经验判断及定义。由于人体结构相当复杂，各种组织或器官对辐射的敏感度也大不相同，因此单单靠人工判断很可能忽略更好的照射角度，而导致治疗效果大打折扣。为了达到治疗效果的优化，射束的照射角度是需要考虑的必要环节。

因此，有必要提出一种射束的照射角度评价方法。

发明内容

为了克服现有技术的缺陷，本发明人发展一种可执行的照射角度优化的方法，可以做为一个有利的参考依据，并搭配医生的经验，尽可能地找出最佳照射角度。优化演算的实施，是分为多直线由肿瘤至体表的径迹，并考虑径迹中的器官比例及其对辐射敏感因素求出个别径迹的评价因素，可以采用正算(由体表至肿瘤)或逆算(由肿瘤至体表)法，可以依序取点或是随机取点计算。通过上述方法，可以算出由体表各点入射肿瘤的评价分数，可以与2D或3D影像进行融合，帮助用户找出最佳入射点。

具体地，本发明的一个方面提供一种射束的照射角度评价方法，包括：

对射束的照射角度进行取样的步骤，其中，所述射束的照射角度定义为射束的照射点至肿瘤预设点的向量方向；

计算射束通过器官径迹的步骤，其中，所述射束的照射角度定义为射束的照射点至肿瘤预设点的向量方向；

判断肿瘤是否充分涵盖于有效治疗深度内：

如果是，则进入计算评价因数、记录照射条件及计算结果的步骤并返回对射束的照射角度进行取样的步骤；

如果否，则进入给予最差评价因数的步骤并返回对射束的照射角度进行取样的步骤。

根据此射束的照射角度评价方法，可以清楚地认识到射束在某一特定位置和某一特定角度照射时的优劣表现，从而为医生或物理师决定照射模式提供强有力的数据支持。

“肿瘤预设点”可设定为肿瘤质心或肿瘤最深处，具体的肿瘤预设点位置可视用户需求而调整。作为一种优选地，所述射束的照射角度定义为射束的照射点至肿瘤质心或肿瘤最深处的向量方向。当然，本领域技术人员熟知地，也可以根据使用需要自定义射束的照射角度。

当然，本领域技术人员熟知的，上述的向量方向包括射束的照射点至肿瘤预设点的正向量方向和负向量方向。

所述射束为中子射束、带电粒子射束或γ射线中的一种或多种，其中带电粒子射束可以为电子射束、质子射束及重粒子射束等。

为了更加精准的计算出所述评价因数，所述评价因数基于射束特性、辐射敏感因数及器官含硼浓度进行计算。

取样的照射角度及照射径迹中，器官(i)的权重因数(W(i))采用公式一进行计算：

W(i)＝I(i)×S(i)×C(i) (公式一)

其中，I(i)、S(i)及C(i)分别为射束强度、器官(i)的辐射敏感因数及器官(i)的含硼浓度。

进一步地，所述I(i)采用根据所用的射束模拟人体的深度强度或剂量曲线积分的公式二进行计算：

(公式二)

其中，i(x)为治疗用射束在近似人体中的深度强度或剂量曲线函数，x₀-x为器官(i)在射束径迹中的深度范围。

所述评价因数采用公式三进行计算：

(公式三)

其中，Q(x,y,z,Φ,θ)作为评价因数等于器官径迹中各器官的权重因数总和。

当然，也可以采用另一种评价因数的呈现方式，即器官评价因数与肿瘤评价因数的比值。所述评价因数与肿瘤评价因数之比(QR(x,y,z,Φ,θ))采用公式四进行计算：

(公式四)

其中，W(tumor)为肿瘤的权重因数。

“医学影像数据”可以为核磁共振成像(Magnetic Resonance Imaging,MRI)或电子计算机断层扫描(Computed Tomography,CT)或正电子发射计算机断层显像(Positron Emissioncomputed Tomography,PET-CT)。但本领域技术人员熟知地，还可以使用其他的医学影像数据，只要该医学影像数据能够定义出器官、组织及肿瘤的材料及密度，就能够应用于本发明揭示的射束的照射角度评价方法中。作为一种优选地，所述读取医学影像数据的步骤为读取CT影像数据或MRI影像数据或PET-CT影像数据中的至少一种医学影像数据的步骤。

“器官径迹”是指在某一特定位置和某一特定角度采用射束照射时通过皮肤、骨头、组织及肿瘤的一条射束径迹，如器官径迹为射束依次序经由皮肤、骨头、组织、肿瘤、组织、骨头、皮肤的径迹。在某一些计算中，如器官径迹中各器官的权重因数加总时，可以不包括肿瘤的权重因数。作为一种优选地，所述器官径迹为射束依次序经由皮肤、头骨、脑组织、肿瘤、脑组织、头骨、皮肤的径迹。这样便可以清晰的了解到射束通过脑部的径迹。当然，本领域技术人员熟知地，器官径迹可以为射束依次序经由人体脑部以外的其他部位，如肝部等。

作为一种优选地，所述射束的照射角度评价方法进一步包括将各个评价因数以3D图像显示的步骤。当然，本领域技术人员也可以通过其他方式显示各个评价因数，只要医生或物理师能够识别显示出的各个评价因数即可。

所谓“射束”可以为一种或多种放射性射束，作为一种优选地，所述射束为中子束及γ射束的混合射束，也可以是单独的中子束、单独的质子束或单独的重粒子束。

所述射束的照射角度评价方法进一步包括：读取医学影像数据的步骤；定义或读取器官、组织及肿瘤的轮廓范围的步骤；定义器官、组织及肿瘤的材料及密度的步骤。

附图说明

图1是硼中子捕获反应示意图。

图2是¹⁰B(n,α)⁷Li中子捕获核反应方程式。

图3是本发明实施例中的射束的照射角度评价方法的逻辑框图。

图4是本发明实施例中的射束照射时的器官径迹示意图。

具体实施方式

下面结合附图对本发明的实施例做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。

作为一种优选地，以用于中子捕获治疗的射束的照射角度评价方法为本发明的实施例。下面将简单介绍一下中子捕获治疗，尤其是硼中子捕获治疗。

中子捕获治疗作为一种有效的治疗癌症的手段近年来的应用逐渐增加，其中以硼中子捕获治疗最为常见，供应硼中子捕获治疗的中子可以由核反应堆或加速器供应。本发明的实施例以加速器硼中子捕获治疗为例，加速器硼中子捕获治疗的基本组件通常包括用于对带电粒子(如质子、氘核等)进行加速的加速器、靶材与热移除系统和射束整形体，其中加速带电粒子与金属靶材作用产生中子，依据所需的中子产率与能量、可提供的加速带电粒子能量与电流大小、金属靶材的物化性等特性来挑选合适的核反应，常被讨论的核反应有⁷Li(p,n)⁷Be及⁹Be(p,n)⁹B，这两种反应皆为吸热反应。两种核反应的能量阀值分别为1.881MeV和2.055MeV，由于硼中子捕获治疗的理想中子源为keV能量等级的超热中子，理论上若使用能量仅稍高于阀值的质子轰击金属锂靶材，可产生相对低能的中子，不须太多的缓速处理便可用于临床，然而锂金属(Li)和铍金属(Be)两种靶材与阀值能量的质子作用截面不高，为产生足够大的中子通量，通常选用较高能量的质子来引发核反应。

硼中子捕获治疗(Boron Neutron Capture Therapy,BNCT)是利用含硼(¹⁰B)药物对热中子具有高捕获截面的特性，借由¹⁰B(n,α)⁷Li中子捕获及核分裂反应产生⁴He和⁷Li两个重荷电粒子。参照图1和图2，其分别示出了硼中子捕获反应的示意图和¹⁰B(n,α)⁷Li中子捕获核反应方程式，两荷电粒子的平均能量约为2.33MeV，具有高线性转移(Linear EnergyTransfer,LET)、短射程特征，α粒子的线性能量转移与射程分别为150keV/μm、8μm，而⁷Li重荷粒子则为175keV/μm、5μm，两粒子的总射程约相当于一个细胞大小，因此对于生物体造成的辐射伤害能局限在细胞层级，当含硼药物选择性地聚集在肿瘤细胞中，搭配适当的中子射源，便能在不对正常组织造成太大伤害的前提下，达到局部杀死肿瘤细胞的目的。

无论硼中子捕获治疗的中子源来自核反应堆或加速器带电粒子与靶材的核反应，产生的皆为混合辐射场，即射束包含了低能至高能的中子、光子；对于深部肿瘤的硼中子捕获治疗，除了超热中子外，其余的辐射线含量越多，造成正常组织非选择性剂量沉积的比例越大，因此这些会造成不必要剂量的辐射应尽量降低。除了空气射束品质因素，为更了解中子在人体中造成的剂量分布，本发明的实施例中使用人体头部组织假体进行剂量计算，并以假体射束品质因素来作为射束的设计参考，将在下文详细描述。

国际原子能机构(IAEA)针对临床硼中子捕获治疗用的中子源，给定了五项空气射束品质因素建议，此五项建议可用于比较不同中子源的优劣，并供以作为挑选中子产生途径、设计射束整形体时的参考依据。这五项建议分别如下：

超热射束通量Epithermal neutron flux>1x10⁹n/cm²s

快中子污染Fast neutron contamination<2x10^-13Gy-cm²/n

光子污染Photon contamination<2x10^-13Gy-cm²/n

热中子与超热中子通量比值thermal to epithermal neutron flux ratio<0.05

中子电流与通量比值epithermal neutron current to flux ratio>0.7

注：超热中子能区在0.5eV到40keV之间，热中子能区小于0.5eV，快中子能区大于40keV。

1、超热射束通量：

射束通量和肿瘤中含硼药物浓度共同决定了临床治疗时间。若肿瘤含硼药物浓度够高，对于射束通量的要求便可降低；反之，若肿瘤中含硼药物浓度低，则需高通量超热中子来给予肿瘤足够的剂量。IAEA对于超热射束通量的要求为每秒每平方厘米的超热中子个数大于10⁹，此通量下的射束对于目前的含硼药物而言可大致控制治疗时间在一小时内，短治疗时间除了对病人定位和舒适度有优势外，也可较有效利用含硼药物在肿瘤内有限的滞留时间。

2、快中子污染：

由于快中子会造成不必要的正常组织剂量，因此视之为污染，此剂量大小和中子能量呈正相关，因此在射束设计上应尽量减少快中子的含量。快中子污染定义为单位超热中子通量伴随的快中子剂量，IAEA对快中子污染的建议为小于2x10^-13Gy-cm²/n。

3、光子污染(γ射线污染)：

γ射线属于强穿辐射，会非选择性地造成射束路径上所有组织的剂量沉积，因此降低γ射线含量也是中子束设计的必要要求，γ射线污染定义为单位超热中子通量伴随的γ射线剂量，IAEA对γ射线污染的建议为小于2x10^-13Gy-cm²/n。

4、热中子与超热中子通量比值：

由于热中子衰减速度快、穿透能力差，进入人体后大部分能量沉积在皮肤组织，除黑色素细胞瘤等表皮肿瘤需用热中子作为硼中子捕获治疗的中子源外，针对脑瘤等深层肿瘤应降低热中子含量。IAEA对热中子与超热中子通量比值建议为小于0.05。

5、中子电流与通量比值：

中子电流与通量比值代表了射束的方向性，比值越大表示射束前向性佳，高前向性的中子束可减少因中子发散造成的周围正常组织剂量，另外也提高了可治疗深度及摆位姿势弹性。IAEA对中子电流与通量比值建议为大于0.7。

利用假体得到组织内的剂量分布，根据正常组织及肿瘤的剂量-深度曲线，推得假体射束品质因素。如下三个参数可用于进行不同射束治疗效益的比较。

1、有效治疗深度：

肿瘤剂量等于正常组织最大剂量的深度，在此深度之后的位置，肿瘤细胞得到的剂量小于正常组织最大剂量，即失去了硼中子捕获的优势。此参数代表射束的穿透能力，有效治疗深度越大表示可治疗的肿瘤深度越深，单位为cm。

2、有效治疗深度剂量率：

即有效治疗深度的肿瘤剂量率，亦等于正常组织的最大剂量率。因正常组织接收总剂量为影响可给予肿瘤总剂量大小的因素，因此参数影响治疗时间的长短，有效治疗深度剂量率越大表示给予肿瘤一定剂量所需的照射时间越短，单位为cGy/mA-min。

3、有效治疗剂量比：

从大脑表面到有效治疗深度，肿瘤和正常组织接收的平均剂量比值，称之为有效治疗剂量比；平均剂量的计算，可由剂量-深度曲线积分得到。有效治疗剂量比值越大，代表该射束的治疗效益越好。

为了使射束整形体在设计上有比较依据，除了五项IAEA建议的空气中射束品质因素和上述的三个参数，本发明实施例中也利用如下的用于评估射束剂量表现优劣的参数：

1、照射时间≤30min(加速器使用的质子电流为10mA)

2、30.0RBE-Gy可治疗深度≥7cm

3、肿瘤最大剂量≥60.0RBE-Gy

4、正常脑组织最大剂量≤12.5RBE-Gy

5、皮肤最大剂量≤11.0RBE-Gy

注：RBE(Relative Biological Effectiveness)为相对生物效应，由于光子、中子会造成的生物效应不同，所以如上的剂量项均分别乘上不同组织的相对生物效应以求得等效剂量。

请参照图3和图4，为提供中子捕获治疗时射束照射角度的参考依据，下面利用随机或逐一计算各个可能的照射角度，并通过公式一和公式二计算出评价因数，将各个评价因数以3D图像的方式显现，即方便医生或物理师判断治疗的照射角度。

下面将详述3D评价因数图像的建立步骤，该步骤为一种优选实施步骤，本领域技术人员熟知地，并不一定要按照优选实施步骤。该步骤具体为：读取患者CT/MRI/PET-CT等具有明确人体解剖的影像，逐一定义各个器官、组织及肿瘤的轮廓，并给设定材料种类及密度。在完成几何材料及密度的定义后，便决定计算射束的起始位置及射束角度，计算中起始位置与角度的决定，可以是顺向算法或逆向算法，顺向算法中是将起始位置决定于体外位置并可依固定角度或距离间隔依序取样计算，也可以透过随机取样的方式进行；射束角度的部分，则可设定为照射点至肿瘤质心或肿瘤最深处的向量方向，具体的肿瘤端点位置可视用户需求而调整；逆向算法中，则是将起始位置决定于肿瘤范围内，其起始位置可以是肿瘤质心、最深处或肿瘤范围内随机取点，而射束角度则可以利用随机取样或依指定间隔取样的方式进行。在决定完射束的照射位置与角度后，则计算射束轴所通过的器官径迹，即计算出射束进入人体后所经过的器官种类及其厚度。在取得射束轴通过人体的径迹信息后，则判断肿瘤范围是否落于最大可治疗深度范围内。若是，则以此径迹信息为依据配合用户设定的器官含硼浓度、器官辐射敏感因数及射束特性信息，计算该径迹的评价因数；若否，则给予最差的评价因数，重新进行射束位置与照射角度的取样。完成评价因数的计算后，则进行照射位置、角度及评价因数的纪录。重复进行上述演算达一定数量后，即可输出报表，提供给医生或物理师作为参考从而决定照射模式，更可以将数据以3D图像的形式表示，以更便利地根据射束轴周围的评价因数判断优劣。

评价因数是基于射束特性、器官辐射敏感因数及器官含硼浓度所进行计算的，对应某位置及某照射径迹中，器官i的权重因数(W(i))计算如公式一所示，式中I(i)、S(i)及C(i)分别为中子强度、器官i的辐射敏感因数及器官i的含硼浓度。

W(i)＝I(i)×S(i)×C(i) (公式一)

其中，I(i)乃根据所用射束于模拟人体的深度强度或剂量曲线积分而得，如公式二所示，式中，i(x)为治疗用射束于近似人体中的深度强度或剂量曲线函数，x₀-x为器官i于射束径迹中的深度范围。

(公式二)

通过上述算法，依序完成器官径迹中各器官的权重因数计算后，并取其加总，即可求得该射束对应的评价因数，如公式三所示，于此计算中，肿瘤的权重因数应不列入计算。

(公式三)

根据上述的评价因数高低，可以更直观地判断治疗时对正常组织所受到的危害程度。除了利用评价因数进行照射位置及角度的评价外，也可利用评价比因数来进行评价，其定义为评价因数对肿瘤权重因数的比值，如公式四所示，如此可以更充分地显现照射位置与角度的预期治疗效果。

(公式四)

以上实施例中涉及到这些步骤——“读取患者CT/MRI/PET-CT等具有明确人体解剖的影像，逐一定义各个器官、组织及肿瘤的轮廓，并给设定材料种类及密度。”可以参照本申请人于2015年11月17日递交到国家知识产权局、申请号为201510790248.7、发明名称为“基于医学影像数据的几何模型建立方法”的专利申请，在此全文引入。

本领域技术人员熟知的，上述公式一至公式四中的一些简单变换，仍然在本发明要求保护的范围之内，如I(i)、S(i)及C(i)由相乘变换为相加；I(i)、S(i)及C(i)分别乘以n次方，n根据情况而定，可以为1的整数倍也可以是其他倍数；i(x)可以是x₀-x间的平均数或中间数乘上(x₀-x)，或任何可以达到与强度积分计算结果相符的计算方法。

尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，都在本发明要求保护的范围之内。

Claims (10)

1.一种射束的照射角度评价方法，包括：

对射束的照射角度进行取样的步骤，其中，所述射束的照射角度定义为射束的照射点至肿瘤预设点的向量方向；

计算射束通过器官径迹的步骤；

判断肿瘤是否充分涵盖于有效治疗深度内：

如果是，则进入计算评价因数、记录照射条件及计算结果的步骤并返回对射束的照射角度进行取样的步骤；

如果否，则进入给予最差评价因数的步骤并返回对射束的照射角度进行取样的步骤。

2.根据权利要求1所述的射束的照射角度评价方法，其特征在于：所述射束为中子射束、带电粒子射束或γ射线中的一种或多种。

3.根据权利要求1所述的射束的照射角度评价方法，其特征在于：所述评价因数基于射束特性、辐射敏感因数及器官含硼浓度进行计算。

4.根据权利要求3所述的射束的照射角度评价方法，其特征在于：取样的照射角度及照射径迹中，器官(i)的权重因数(W(i))采用公式一进行计算：

W(i)＝I(i)×S(i)×C(i) (公式一)

其中，I(i)、S(i)及C(i)分别为射束强度、器官(i)的辐射敏感因数及器官(i)的含硼浓度。

5.根据权利要求4所述的射束的照射角度评价方法，其特征在于：所述I(i)采用根据所用的射束模拟人体的深度强度或剂量曲线积分的公式二进行计算：

$I\left(i\right)={\&Integral;}_{x_0}^{x}i\left(x\right)dx$ (公式二)

其中，i(x)为治疗用射束在近似人体中的深度强度或剂量曲线函数，x₀-x为器官(i)在射束径迹中的深度范围。

6.根据权利要求4或5中任一项所述的射束的照射角度评价方法，其特征在于：所述评价因数采用公式三进行计算：

$Q(x,y,z,\mathrm{\&phi;},\mathrm{\&theta;})=\underset{i}{\&Sigma;}W\left(i\right)$ (公式三)

其中，Q(x,y,z,Φ,θ)作为评价因数等于器官径迹中各器官的权重因数总和。

7.根据权利要求6所述的射束的照射角度评价方法，其特征在于：所述评价因数与肿瘤评价因数之比(QR(x,y,z,Φ,θ))采用公式四进行计算：

$QR(x,y,z,\mathrm{\&phi;},\mathrm{\&theta;})=\underset{i}{\&Sigma;}W\left(i\right)/W\left(tumor\right)$ (公式四)

其中，W(tumor)为肿瘤的权重因数。

8.根据权利要求1-5中任一项所述的射束的照射角度评价方法，其特征在于：所述读取医学影像数据的步骤为读取CT影像数据或MRI影像数据或PET-CT影像数据中的至少一种医学影像数据的步骤。

9.根据权利要求1-5中任一项所述的射束的照射角度评价方法，其特征在于：所述射束的照射角度评价方法进一步包括将各个评价因数以3D图像显示的步骤。

10.根据权利要求1-5中任一项所述的射束的照射角度评价方法，其特征在于：所述射束的照射角度评价方法进一步包括：读取医学影像数据的步骤；定义或读取器官、组织及肿瘤的轮廓范围的步骤；定义器官、组织及肿瘤的材料及密度的步骤。