粒子射线扫描照射系统及粒子射线扫描照射方法
技术领域
本发明涉及一种粒子射线治疗装置,尤其涉及一种与患部的三维形状相配合来照射粒子射线的粒子射线扫描照射系统。
背景技术
在粒子射线治疗中,对患部照射加速到光速的约70%的质子射线或碳射线等。在将这些高能粒子射线照射到体内的肿瘤等时,会表现出以下特征。第一,所照射的粒子射线大部分都停止在如下位置,即,与粒子射线能量的约1.7次方成比例的深度处的位置。第二,所照射的粒子射线在体内停止之前对其所通过的路径提供的能量密度在粒子射线的停止位置达到最大值。将粒子射线的能量密度称为剂量。将沿着粒子射线通过体内的路径而形成的特有的深部剂量分布曲线称为布拉格曲线(Bragg curve)。
将粒子射线的剂量表现出最大值的部位称为布拉格峰(Bragg peak)。粒子射线扫描照射系统使该布拉格峰的位置与肿瘤的三维形状相一致来进行扫描。以在预先利用图像诊断所决定的目标(肿瘤区域)中形成规定的三维剂量分布的方式,调整各扫描位置的峰值剂量。
粒子射线停止位置的扫描方法分为沿大致垂直于粒子射线的照射方向的横向(X、Y方向)扫描的方法、以及沿粒子射线的照射方向即深度方向(Z方向)扫描的方法。在沿横向进行的扫描中,存在使患者相对于粒子射线移动的方法、以及使用电磁铁等使粒子射线的位置移动的方法。一般采用后者的使用电磁铁的方法。
对于沿深度方向进行的扫描方法,唯有改变粒子射线能量这一方法。在改变能量的方法中,可考虑如下两种方法,即:利用加速器改变粒子射线的能量的方法、以及使用设置在射束传输系统或照射系统且被称为射程移位器(range shifter)的能量变更装置的方法。目前多采用使用能量变更装置的方法。射程移位器有时也包含进行能量的分析与运动量的选择且被称为能量选择系统(Energy Selection System)的装置。
使粒子射线沿横向扫描的方法分为点扫描照射法及混合扫描照射法这两种基本照射方法。点扫描照射法照射粒子射线,当规定照射位置的照射量达到计划值时,暂时减弱粒子射线的射束强度(参照非专利文献1)。此时,一般来说,粒子射线的射束强度为零。为了将粒子射线照射到下一照射位置,变更扫描电磁铁的电流值,再次增加粒子射线的强度之后,照射粒子射线。也可不增加粒子射线的射束强度,而是从加速器再出射粒子射线。
混合扫描照射法中将粒子射线以计划量照射到计划位置的基本方法与点扫描照射法相同,但当使粒子射线的位置移动到下一照射位置时,不使粒子射线停止,而是一边照射粒子射线,一边使粒子射线进行扫描(参照非专利文献1)。
背景技术文献
专利文献
专利文献1:日本专利特开2006-145213号公报
专利文献2:日本专利特开2008-136523号公报
专利文献3:日本专利特开2011-156340号公报
非专利文献
非专利文献1:T.Inaniwa等,医疗物理,34(2007)3302(T.Inaniwa et al.,Medical Physics,34(2007)3302)
发明内容
[发明要解决的问题]
在粒子射线扫描照射系统中,在变更照射位置的过程中产生粒子射线的移动。移动中的粒子射线的剂量会对实际照射时的剂量分布产生影响。在非专利文献1中,想要通过将粒子射线移动中的剂量贡献引入到治疗计划的最优化计算,以此来减轻移动中的剂量对最终的剂量分布产生的影响。在此方法中,要求预先决定移动中的剂量贡献。
实际照射时的移动中的剂量贡献取决于在粒子射线的照射中从加速器获得的粒子射线的射束强度的时间变动,因此,在所述方法中,需要将移动中的剂量贡献的平均值引入到治疗计划的最优化计算。在从加速器获得的射束电流波形I(t)的时间变动较大的情况下,难以高精度地考虑移动中的粒子射线的贡献。进而,在治疗计划的最优化计算中需要考虑移动中的剂量贡献的平均值,这使治疗计划的制定变得相对复杂。
本发明是为了解决所述问题而设计的,其目的在于减小治疗计划中所决定的各照射点(照射位置)的计划照射粒子数与实际照射时的剂量分布之差。
[解决问题的技术方案]
本发明的粒子射线扫描照射方法的特征在于包括:第一步骤,根据由粒子射线治疗计划决定的相对粒子射线照射量及处方粒子射线剂量,求出粒子射线的各照射点的计划照射粒子数;第二步骤,根据计划照射粒子数及粒子射线的射束电流波形来模拟各照射点的粒子射线的照射过程,求出粒子射线在扫描移动期间照射到患部的粒子数;第三步骤,使用在扫描移动期间照射的粒子数,校正各照射点的计划照射粒子数;第四步骤,将校正后的计划照射粒子数换算成剂量监控器的计数值;以及第五步骤,根据换算得到的计数值照射粒子射线。
[发明的效果]
根据本发明的粒子射线扫描照射系统,能够以简单的校正方法使实际照射时的剂量分布接近由治疗计划决定的各照射点(照射位置)的计划照射粒子数。
附图说明
图1是表示本发明的粒子射线扫描照射系统的装置结构的整体图。
图2是说明患部被切片的状态的图。
图3是表示照射开始点与照射结束点的位置关系的图。
图4是表示本发明的实施方式1的粒子射线扫描照射方法的步骤的框图。
图5是表示本发明的实施方式1的扫描照射仿真所使用的射束电流波形的图。
图6是表示本发明的实施方式2的扫描照射仿真所使用的经平滑化的射束电流波形的图。
图7是说明移动中的照射点的情况的图。
图8是表示本发明的实施方式3的粒子射线扫描照射方法的步骤的框图。
图9是说明本发明的实施方式4及实施方式5的粒子射线扫描照射系统的动作的图。
具体实施方式
实施方式1.
图1表示粒子射线扫描照射系统的整体结构。粒子射线扫描照射系统100包括治疗计划系统10、粒子射线治疗照射控制系统20、及粒子射线照射装置30。治疗计划系统10包括治疗计划数据管理装置1、计划照射粒子数校正装置3、治疗计划装置4、及扫描照射仿真装置5等。粒子射线治疗照射控制系统20包括照射系统数据管理装置2、治疗控制装置6、设备控制装置7、及定位装置8。治疗控制装置6包含剂量监控器。粒子射线照射装置30包含扫描磁铁及扫描电源。
治疗计划装置4生成治疗计划,并根据治疗计划对剂量计算进行仿真。治疗控制装置6控制粒子射线照射装置30,按照从照射系统数据管理装置2获取的治疗计划的条件来照射粒子射线。实际照射的粒子射线的剂量由治疗控制装置6测定。将测定的结果发送到治疗计划数据管理装置1。治疗计划数据管理装置1对计划照射粒子数校正装置3、治疗计划装置4、扫描照射仿真装置5所生成的数据进行管理。照射系统数据管理装置2管理治疗控制装置6、设备控制装置7、定位装置8所使用的数据、治疗记录、及测定记录等。
因为在粒子射线治疗中使粒子射线的布拉格峰的位置与肿瘤的三维形状相一致来进行扫描,所以按照治疗计划将患部沿深度方向较薄地切成圆片。图2示意性地表示被切成圆片的肿瘤等患部25。将沿深度方向分割的患部25的一层称为切片。患部25被分割成多个切片(……W、X、Y、Z……)。
接下来,根据图3,对扫描照射时的粒子射线的照射顺序进行说明。粒子射线的照射顺序预先利用治疗计划系统10而计算。图中表示特定的切片X中的照射点的配置及粒子射线的照射路径。粒子射线原则上从照射开始点到照射结束点为止沿着一条单行路径行进。粒子射线通常通过扫描移动而移到相邻照射点。在照射点之间远离的情况下,粒子射线为了越过照射点之间而进行空白移动(blank shift),从而移到下一照射点。因空白移动而省去照射的点的数量根据照射条件而不同。空白移动从空白起点这一点开始且在空白终点这一点结束。
接下来,根据图4所示的框图,对本发明的实施方式1的粒子射线扫描照射方法进行说明。计划照射粒子数Ni表示由治疗计划所决定的照射点(i)的照射粒子数的计划值(处方粒子数)。点编号i(=1、2、3……Nspot)表示照射点的ID(identification,标识符)。粒子射线在同一切片内按照点编号i进行照射。根据治疗计划求出计划照射粒子数Ni。根据治疗计划系统10所输出的各点的相对照射粒子数Ri及处方剂量DO算出计划照射粒子数Ni,并将其存储到治疗计划数据管理装置1。处方剂量DO表示在一次照射中照射到患部的粒子射线的剂量。
射束电流波形I(t)存储在照射系统数据管理装置2,根据需要每日进行更新,以反映最新的测定结果。图5表示扫描照射时所使用的射束电流波形I(t)的一例。射束电流波形I(t)是涉及从加速器射出的粒子射线的射束电流的基本时间变动信息。用于扫描照射的加速器无论为以同步加速器、回旋加速器为代表的何种类型的加速器,在本发明中都同样有效。扫描照射仿真装置5算出粒子射线在点间移动过程中照射到患部的粒子射线数的预想数据。通过计划照射粒子数校正装置3校正计划照射粒子数Ni。治疗计划包含扫描照射治疗计划。经校正的计划照射粒子数Ni被换算成MU(Monitor Unit,监控单位)以对应于剂量监控器。
接下来,使用图1~图4,对本发明的实施方式1的粒子射线扫描照射系统100的动作进行说明。首先,根据患者的CT(Computed Tomography,计算机断层摄影术)数据等,利用治疗计划系统10来生成扫描照射治疗计划。结果,输出粒子射线的照射方向、对应于照射方向的照射切片信息、对应于各切片的照射位置信息、各照射位置的点的相对照射量信息。也考虑有多个照射方向的情况。在照射位置信息中,对应于照射位置而定义照射点。
接着,治疗计划装置4根据处方剂量DO与治疗计划,算出计划照射粒子数Ni。在算出计划照射粒子数Ni时,利用治疗计划中所包含的各点的相对照射粒子数、及水中的剂量分布绝对值的测定信息。对扫描照射仿真装置5提供各照射点的计划照射粒子数Ni及射束电流波形I(t)。
扫描照射模拟装置5包含仿真软件及计算机。扫描照射仿真装置5在实施扫描照射的情况下,算出粒子射线从照射点(i)朝照射点(i+1)移动时照射到患部的粒子射线的数量ΔNi(i=1、2、3、……Nspot-1)。当对扫描照射过程进行仿真时,将剂量监控器、扫描电源、扫描控制装置等的时间响应及动作考虑在内。
在扫描照射仿真装置5中,以几μsec的时间步长(time step)将射束电流波形I(t)积分。当电流积分值达到照射点(i)的计划照射粒子数Ni时,将各设备的时间响应考虑在内,将粒子射线的照射位置朝照射点(i+1)变更。进而,针对照射点(i+1),将射束电流波形I(t)积分,一旦电流积分值达到计划照射粒子数Ni+1,便将照射位置变更到照射点(i+2)。通过以此方式模拟对从照射点(i=1)到照射点(i=Nspot)的所有照射点照射粒子射线的过程,可获得射束在照射点之间移动期间照射到患部的粒子射线的数量(移动中粒子数ΔNi)。在算出计划照射粒子数Ni时并未考虑移动中粒子数ΔNi。在实际照射时,因为是在达到计划照射粒子数Ni之后才发生射束的移动,所以在粒子射线移动到下一照射点之前,粒子射线会多余地照射到患部(参照图7)。
接着,计划照射粒子数校正装置3使用移动中粒子数ΔNi来校正照射点(i)的计划照射粒子数Ni。移动中粒子数ΔNi是照射到照射点(i)与照射点(i+1)之间的粒子数。具体来说,实施(Ni_校正后=Ni-ΔNi/2)及(Ni+1_校正后=Ni+1-ΔNi/2)。根据该处理求出计划照射粒子数Ni_校正后。计划照射粒子数Ni_校正后是将在点之间的移动中所照射的粒子数的影响考虑在内的计划照射粒子数。如果着眼于照射点(i),那么下式成立(参照图7)。此外,移动中粒子数ΔNi-1是照射到照射点(i-1)与照射点(i)之间的粒子数(1<i<Nspot)。在照射点(1)与照射点(Nspot)的情况下,校正算式分别变成(Ni-0.5×ΔNi)及(Ni-0.5×ΔNi-1)。
计划照射粒子数Ni_校正后=计划照射粒子数Ni-ΔNi/2-(ΔNi-1)/2
使用在照射中管理(或计量)粒子数的剂量监控器的灵敏度,将计划照射粒子数Ni_校正后换算成计划计数值MUi。剂量监控器的灵敏度取决于粒子射线的能量。最后,将获得的各点的计划计数值MUi交付给粒子射线治疗照射控制系统20及粒子射线照射装置30,实施实际的扫描照射。当使用粒子射线治疗照射控制系统20执行扫描照射时,使用强度及基本时间模式与利用扫描照射仿真装置5模拟扫描照射时所使用的加速器射束相同的加速器射束。
根据计划计数值MUi来控制实际照射到照射点(i)的粒子射线的量。图7示意性地表示在切片Y中对照射粒子射线进行计数的情况。射束以照射点(i-1)、照射点(i)、照射点(i+1)的顺序移动。照射粒子射线的计数在起点Ti(s)到终点Ti(e)之间为一个循环。连续地重复该计数的循环。在实际照射时,当计数开始时,利用在上一循环中移动中的粒子射线进行照射。当计数结束时,粒子射线开始移动到下一照射点。
像这样,根据本发明的实施方式1所示的粒子射线扫描照射系统,在扫描照射时,将在点之间的移动期间中照射到目标的移动中粒子对剂量分布的影响考虑在内。使用可模拟实际的照射定时的扫描照射仿真装置,在照射前算出该移动中粒子数ΔNi,使用计划照射粒子数校正装置求出校正后的计划照射粒子数(计划照射粒子数Ni_校正后)来进行实际的扫描照射,由此无需利用治疗计划系统进行再计划处理。结果,能够简化粒子射线扫描照射的治疗计划,所以使粒子射线扫描照射系统的运用效率得以提高。另外,能够考虑移动中粒子的影响而实施扫描照射,所以可更高精度地对患部提供剂量。
实施方式2.
接下来,对本发明的实施方式2的粒子射线扫描照射系统的动作进行说明。实施方式2的特征点在于:作为扫描照射仿真装置5所使用的射束电流波形I(t),使用像图6所示那样的去除快速变动分量后的射束电流波形,以代替图5所示的实际测量到的射束电流波形。
图6所示的射束电流波形I(t)是使图5所示的射束电流波形I(t)通过低通滤波器来进行平滑处理后的波形。结果,可求出在点之间移动时的移动中粒子数ΔNi以作为更稳定的射束电流波形的平均值。由此,不会受到缺乏再现性的快速射束电流变化的影响,从而可更高精度地求出移动中粒子数ΔNi。实施方式2的粒子射线扫描照射系统具有能够以更高精度对患部照射粒子射线的效果。
实施方式3.
接下来,对本发明的实施方式3的粒子射线扫描照射系统的动作进行说明。如果射束电流较大,那么便可缩短照射时间,但移动中粒子数ΔNi会占到相对较大的权重,所以认为计划照射粒子数Ni_校正后会成为负值。图8是表示本发明的实施方式3的粒子射线扫描照射系统的构成的框图。图8所使用标号中与图4相同的标号基本上表示与图4相同的意思。
在实施方式3中,为了防止计划照射粒子数Ni_校正后成为负值,在对各切片进行照射时,以使移动中粒子数ΔNi的一半变得小于移动前后的点的计划照射粒子数的方式,限制在各切片的照射中所使用的最大射束电流的大小。对于射束电流值,实际上已经决定了其波形,所以在使用经平滑化后的平均值的情况下也相同。
具体来说,当进行扫描照射仿真时,以使移动中粒子数ΔNi的一半(ΔNi/2)不超过照射点(i)的计划照射粒子数Ni及照射点(i+1)的计划照射粒子数Ni+1的方式,预先决定在各切片的照射中所使用的射束电流波形I(t)的最大值(某时刻的最大射束电流)。在对切片进行照射时,使用较大射束电流值更能缩短照射时间,但作为标准,此处使用算出的最大射束电流值的一半以下的射束电流进行照射即可。
在实施方式3中,因为预先限制了对各切片进行照射时所使用的射束电流的最大值,所以可将对各切片进行照射时所使用的射束电流的范围设定为适当值。当对计划照射粒子数进行校正时,不会使计划照射粒子数Ni_校正后成为负值,能够保証其一定为正值。因此,与将负粒子数设为零而进行照射的情况相比,能够以更高精度、且将误差粒子数的影响考虑在内的方式来进行粒子射线扫描照射。
实施方式4.
接下来,对本发明的实施方式4的粒子射线扫描照射系统的动作进行说明。图9表示当将患部沿深度方向切片并扫描照射主干部时特定的切片中的照射点的配置及照射路径。扫描照射从位于最上部的右端的点开始,在位于最下部的右端部的点结束。照射开始点与照射结束点的途中经由复杂的路径而对切片内的各点进行照射。图中存在连接邻接点之间的点间移动、及连接不邻接点之间的点间移动。
在实施方式4中,着眼于连接邻接照射点的点间移动、及连接不邻接点的点间移动。在邻接的点之间移动的情况下,以实施方式1所说明的方式进行,但在不邻接的点之间发生移动的情况下,如果使用实施方式1的校正方法的话,那么误差会较大。因此,在实施方式4中,首先,计划照射粒子数校正装置3确认是在邻接点之间的移动、或者是在不邻接的点之间的移动。只针对邻接点之间的移动,才实施如下校正处理:从移动前的照射点(i)的计划照射粒子数Ni与移动后的照射点(i+1)的计划照射粒子数Ni+1中减去移动中粒子数ΔNi的一半。
根据实施方式4,在计划照射粒子数Ni的校正处理中,只在粒子射线在邻接点间移动的情况下,才从计划照射粒子数Ni中减去移动中粒子数ΔNi的一半。在粒子射线在不邻接的点间移动的情况下,不实施任何处理,所以不会产生因校正处理而导致的误差。实际上,粒子射线在不邻接的点间移动的情形与粒子射线在邻接点间移动的情形相比非常少。根据具体情形来应对(case by case)的处理能够实现更高精度的扫描照射。
实施方式5.
接下来,对本发明的实施方式5的粒子射线扫描照射系统的动作进行说明。在实施方式5的粒子射线扫描照射系统中,在图9所示的邻接点之间的移动(扫描移动)、及不邻接的点之间的移动(空白移动)的情况下,使用计划照射粒子数校正装置3,实施不同的计划照射粒子数校正处理。具体来说,如果是邻接的点之间的移动的话,那么与实施方式1同样地实施如下校正处理:从移动前的照射点(i)的计划照射粒子数Ni及移动后的照射点(i+1)的计划照射粒子数Ni+1中减去移动中粒子数ΔNi的一半。
在不邻接的点之间的移动的情况下,首先算出移动轨迹。根据所算出的移动轨迹推算出以与点尺寸同等程度隔开的所有点(i=ik,k=1、2、3、……nk;共nk个;i1=i,ink=i+1)。将推算出的点的数量设为空白点数(nk)。接着,进行如下校正:将移动中粒子数ΔNi除以空白点数(nk),从点编号i的移动前点、点编号i+1的移动后点及移动中作出贡献的空白点的计划照射粒子数N中分别减去ΔNi/nk。
根据本实施方式,在不邻接的点之间的移动的情况下也能实施更准确的移动中粒子数的校正,所以与实施方式4相比,可进一步提高扫描照射的精度。此外,在所述实施方式的说明中,将扫描照射仿真装置、治疗计划装置、计划照射粒子数校正装置作为分开的部件而记载并予以说明。实际上,即使将这些装置搭载于同一计算机编码或同一计算机,效果也与所述说明的效果相同,在扫描照射中,能够以高精度且简单的方法将由移动中粒子造成的对剂量的影响考虑在内。
此外,本发明可在其发明范围内,将各实施方式自由组合,或者将各实施方式适当变形、省略。
[标号说明]
1 治疗计划数据管理装置
2 照射系统数据管理装置
3 计划照射粒子数校正装置
4 治疗计划装置
5 扫描照射仿真装置
6 计划照射粒子数校正装置
10 治疗计划系统
20 粒子射线治疗照射控制系统
30 粒子射线照射装置
100 粒子射线扫描照射系统