姚凯宁,王若曦,岳海振,王美娇,王翰林,弓健,张艺宝,吴昊
北京大学肿瘤医院暨北京市肿瘤防治研究所放疗科,恶性肿瘤发病机制及转化研究教育部重点实验室,北京100142
放射治疗是肿瘤治疗的重要手段之一[1],随着放疗设备和技术的持续更新,现已进入精确放疗的时代。影响精确放疗的因素较多,其中贯穿整个治疗流程的患者体位固定问题引起了国内外研究者的重视,尤其对于生理差异较大的头颈部肿瘤。临床实践中,体位固定装置应具备固定性强、精确性高、可重复性高、舒适度好等特征。常用的头颈部体位固定装置包括热塑膜[2‐3]、发泡胶[4]和双面膜固定系统[5‐7]。已有研究表明相对于其他固定方式,双面膜固定系统具有操作简便、固定效果和重复性较好等优势[5]。
Ⅴarian公司Eclipse计划系统在患者定位CT图像导入时按照默认搜索策略自动生成外轮廓(Body),其他部分不参与整个剂量计算过程(CT 值设置为0 HU)。由碳纤维固定床板引入的外轮廓体积的变化对于射线衰减和剂量沉积影响较大[8‐10],而对于双面膜本身引入的体积差异对患者的剂量学影响尚未引起较多关注。Eclipse 计划系统自动生成的外轮廓体积则取决于相关参数的设置,包括搜索范围、搜索CT值阈值、平滑度等,临床常采用默认设置或根据需要加以修改。常规热塑膜体积小、厚度薄,采用常规参数设置,其对剂量计算的影响忽略不计[9],而双面膜体积较大且各部分密度略有差异,因此不同搜索参数下外轮廓的体积差异较大,对剂量计算的影响也相应增加。鉴于Eclipse 计划系统对外轮廓搜索范围的强依赖性和双面膜本身的结构复杂性,加之搜索策略和参数设置的多样化,本研究分析了采用双面膜固定的肿瘤患者的外轮廓自动生成参数设置导致的体积差异对靶区和危及器官剂量学计算的影响,为医学物理师在治疗计划方案制定前的外轮廓范围界定提供必要参考。
随机选取25例采用双面膜固定的全脑放疗病例,采用Siemens 大孔径CT 定位,靶区和危及器官的勾画均由高年资医师确认,计划的设计和制定在Eclipse 计划系统(15.5 版本)上完成,均采用6 MⅤ高能光子束,布野方式为左右两个对穿野,剂量计算算法为Acuros XB 15.5。25例患者的处方剂量为30 Gy/10 次,本研究主要对比外轮廓的体积和剂量学差异以及计划靶区(PTⅤ)、眼球、晶体、视神经等的剂量学差异。
外轮廓搜索条件的更改基于原始CT 图像,保证相同的搜索范围、优化条件和正向剂量计算算法,同时每种计划只改变一个搜索条件。Eclipse 中外轮廓搜索的默认参数如下:平滑度(Smoothing Level)为3,CT 值阈值(Ranger)为‐350 HU、后处理方法(Post Processing)全部选择,默认搜索参数的计划简称为Default 计划。本研究分别改变平滑度、CT 值阈值和后处理方法,生成对比计划。将平滑度分别设置为0、1、2、4(以下分别简称为DSL0、DSL1、DSL2和DSL4计划);CT 值设置为‐550 和‐750 HU(以下简称为HU‐550 和HU‐750 计划);后处理方法中主要探索不使用Disconnect 和Keep number of largest parts(1)(以下分别简称为NoDisconnect 和NoKeep 计划),Disconnect指不连接外轮廓的各部分,Keep number of largest parts(1)指保留体积最大的一部分。以上8 种计划分别与Default 计划进行比较。部分计划的外轮廓如图1所示。对B~F 几种不同搜索策略条件与默认设置A 的体积差异做对比和统计学分析。
图1 不同类型计划中外轮廓(绿线)的范围Fig.1 Range of Body(green line)in different plans
利用Eclipse 脚本应用程序编程接口(ESAPI)导出患者的剂量体积直方图(DⅤH),统计外轮廓、PTⅤ和左右眼球的剂量学差异(D‐Ⅴalue)。计算公式如下:
其中,DRef代表参考计划的剂量,DTest代表对比计划的剂量。
体积差异的计算方式与剂量学相同。靶区比较参数包括剂量的最大值(Dmax)、平均值(Dmean)、最小值(Dmin)以及靶区95%体积的覆盖剂量(D95%);危及器官比较参数有剂量的Dmax和Dmean。统计各参数的均值、标准差和P值。
采用SPSS 24.0 进行统计学分析,对每组数据进行正态性检验,数据符合正态分布,采用配对样本T检验,反之采用相关样本非参数检验,P<0.05 表示差异具有统计学意义。数据图表采用SigmaPlot 10.0软件完成。
不同搜索策略条件对外轮廓体积和Dmean、Dmax等评价指标差异的箱体图分析如图2所示,其中HU‐750 的体积差异乘以0.2 后再进行分析。由图2可知,相对于默认搜索方式,平滑度参数对外轮廓体积、Dmean、Dmax等差异呈负相关趋势,提示外轮廓的搜索方式宜选取较小的平滑度;HU 阈值的绝对值与上述差异呈正相关趋势,提示宜选取绝对值较大的阈值设置;NoKeep 参数对外轮廓体积和Dmean的差异影响 较 小。除NoDisconnect 和NoKeep 的Dmax以及NoDisconnect的Dmean外,各类计划中外轮廓体积与剂量学参数的差异均有统计学意义(P<0.05),如表1所示。
图2 不同搜索策略条件下外轮廓体积和剂量学指标差异Fig.2 Differences in the volume and dosimetric indexes of Body under different search strategies
表1 8种计划类型中外轮廓体积和剂量变化及靶区和眼球的剂量变化的P值统计Tab.1 P values of changes in the dose and volume of Body and the changes in dose of PTV and eyeballs in 8 plans
不同搜索策略下引入的外轮廓体积差异对PTⅤ关键剂量学评价指标的影响见图3。平滑参数DSL与PTⅤ的Dmax、Dmean和D95%均呈负相关趋势,且均具有显著的统计学意义(P<0.05),与外轮廓的剂量学影响较为一致,而对Dmin的影响没有统计学意义。CT值阈值设置除Dmin之外的评价指标差异呈正相关趋势,除NoDisconnect计划中Dmean和D95%外,NoDisconnect和NoKeep两个参数对PTⅤ其他评价指标的差异有统计学意义但影响较小。对于左右眼球的Dmax,搜索策略引入的外轮廓体积差异对其影响较小,与布野方向外轮廓的微小变化相关性很大(图4)。
图3 不同搜索策略条件下PTV剂量学指标差异Fig.3 Differences in PTV dosimetric indexes under different search strategies
图4 不同搜索策略条件下左右眼球的剂量学指标差异Fig.4 Differences in dosimetric indexes of left and right eyeballs under different search strategies
目前已知影响放疗计划剂量的因素较多,包括体位固定装置、靶区和危及器官的勾画[11‐12]、加速器机械参数的设置、优化算法[13‐14]、数据建模等。现有研究表明,体位固定装置不同,引起射束剂量的衰减存在差异,进而影响靶区剂量,但对于外轮廓变化引入的剂量学影响未做有效探讨[9]。本研究通过对比双面膜固定时治疗计划系统外轮廓搜索策略所引入的外轮廓体积变化对剂量计算的影响,从理论计算方面进行了定量分析,重点关注在整个放疗流程中较易忽略且影响整个放疗精度的环节,以适应新的固定装置临床应用中面临的剂量精度评价问题。
本研究结果表明,搜索策略的平滑度参数对各评价指标的影响呈负相关,表明平滑度较低的情况下生成的外轮廓更接近结构较为复杂的双面膜固定方式的实际结构情况。对于HU 限制的设置,本研究表明HU 值绝对值较小时引入的剂量学变化较小,该结论与已有研究中HU 值对剂量计算的影响一致[15‐16]。Disconnect 和Keep number of largest parts(1)两个搜索参数对全脑放疗的靶区和危及器官评价指标的影响较小。此外,危及器官的剂量限值要考虑射线敏感性,例如头颈部的危及器官相对较多且对剂量较为敏感,所以临床计划是否可行必须要考虑危及器官的剂量[17‐20],这些器官的放疗剂量应格外注意。本研究发现危及器官的剂量变化不明显,这可能是由于全脑病例涉及的危及器官较少且处方剂量较小。考虑到双面膜的体积较大且对剂量有一定影响,因此建议外轮廓搜索时应适当包含双面膜的各部分结构。
由于头部解剖结构的特殊性,内部有许多空腔且与外界相通,在搜索外轮廓时会出现外轮廓从鼻孔和外耳道位置深入人体内部的情况,这与搜索条件平滑度密切相关。在临床中通常采取将其修回鼻尖和外耳道的外侧,本研究对此也进行了适当的修改。
Eclipse 计划系统所提供的搜索策略中,平滑度、HU 阈值对全脑放疗病例的剂量学计算结果影响较大,而后处理方法影响较小。针对双面膜固定病例外轮廓界定的特殊情况,本研究建议平滑度为1,HU阈值为‐550 HU。但本研究属于纯理论计算的对比性研究,仅对放疗相关工作人员在治疗计划制定阶段起到提示性和参考性质的作用,至于以上搜索参数如何优化配置更能接近临床实际,需要研究者进一步设计合理的实验,获取剂量方面的参考数据,从而定量和定性分析并确定双面膜固定病例的外轮廓搜索策略。