与传统物理楔形板不同,扩充型动态楔形板(enhanced dynamic wedge, EDW)通过加速器独立准直器的运动和剂量率的变化来实现特定的剂量分布。EDW大大减少了照射时间和跳数,提高靶区剂量的均匀性,同时降低正常组织例如对侧乳腺的受照剂量[1-2]。但当EDW的楔形角度较大且为非对称照射野时,使用标准计算模型计算得到的楔形因子误差较大[3-5]。研究报道指出,在输出分割模型基础上进行相应的计算校正,可有效减少不同射野下楔形因子的计算误差[6-8]。但以上研究均针对射野几何中心处计算模型的修正和剂量计算。对常规二维治疗计划来说,临床应用更多的是计算射野中心点处的楔形因子以及给予相应的射野跳数。因此,本文将基于加速器输出分割模型提出一种新的校正方法,并应用于常规二维计划的中心跳数计算。
1 材料与方法 1.1 EDW楔形因子手工计算方法瓦里安EDW每档能量提供一套楔形角60°最大射野下加速器输出剂量对准直器位置的分配表(Golden segmented treatment tables, Golden STTs) [9],包括楔形角10°,15°,20°,25°,30°,45°和60°,最大楔形野为40 cm× 30 cm,最大可过中线10 cm。楔形因子定义为10 cm深处楔形野和相同大小开野在射野中心点剂量率之比。计算EDW楔形因子的常用方法是加速器输出分割模型(MU fraction model),见公式(1)。其中θ为楔形角,w60、w0为构成楔形角θ所需的60°楔形野、平野权重[10],YF为固定准直器位置坐标,Yf=YF-0.5。STTθ (Y)、STTG (Y)是楔形角θ°,60°时运动准直器坐标Y对应的STTs值,当Y=0时,式(1)计算结果即为射线中心点的楔形因子EDWF (EDW factor, EDWF)。
$ EDW{F_\theta } = \frac{{ST{T_\theta }(Y)}}{{ST{T_\theta }({Y_f})}} = \frac{{{w_\theta }[ST{T_G}(Y)] + {w_0}[ST{T_G}(0)]}}{{{w_\theta }[ST{T_G}({Y_f})] + {w_0}[ST{T_G}(0)]}} $ | (1) |
$ {w_\theta } = \tan \theta /\tan \;60\;\;\;\;\;{w_0} = 1-{w_\theta } $ | (2) |
采用0.125 cm3指形电离室(PTW31010,Freiburg)和UNIDOSE剂量仪(PTW公司),密度为1.045g/cm3的RW3固体水(IBA公司),在Varian Trilogy 5345加速器测量水下10 cm深度处中心轴的楔形因子。对称野分别为5 cm×5 cm、10 cm×10 cm、15 cm×15 cm、20 cm× 20 cm。非对称野则在上述对称野基础上,保持射野大小不变,楔形方向上几何中心以步长为1 cm分别向Y1,Y2方向偏移而形成。为了保证足够的侧向电子散射,各个射野边缘与等中心点间最短距离不小于2 cm。测量角度包括楔形角15°、30°、45°、60°,能量为6 MV、10 MV,楔形方向为Y1-in、Y2-out。为减少杆效应对测量的影响,电离室长轴平行于X方向准直器。
1.3 剂量跳数计算采用了两种方式计算射野中心点的跳数:(1)纯手工计算跳数MU (Monitor Unite, MU)。所采用的常规手工跳数计算公式如式(3) [11];(2)在治疗计划系统Eclipse (version 10.0, Varian公司)上计算射野中心点的剂量跳数。两种计算与实际测量结果的相对误差为:相对误差=[MUmea-MUCal(或MUTPS)]/MUMea。
$ MU = \frac{{Dose}}{{TMR \times WF \times TF \times Sc(FS) \times CF \times InvSqaCorr}} $ | (3) |
计算公式(1)计算结果与测量结果相对误差=(EDWFMea-EDWFCal) /EDWFMea,分析观察相对误差与射线能量、楔形角度、射野大小、射野边缘与等中心点最短距离等之间的关系及变化趋势,结合实际临床允许误差,从而得到不同楔形角度、射野大小、射野边缘与等中心最短距离情况下对应的校正系数。计算校正后的相对误差并将校正后的楔形因子运用于常规二维计划射野中心点跳数计算。
2 结果 2.1 EDW楔形因子计算结果以方向Y1-in,能量6 MV为例,与测量结果对比,就对称野而言,同一射野大小相对误差随角度增大而增大,同一楔形角度相对误差随射野增大而增大,如图 1。射野20 cm×20 cm、楔形角60°最大相对误差为4.2%,楔形角45°最大相对误差为2.5%,楔形角30°最大相对误差为1.3%。非对称射野相对误差与楔形角度、射野大小、射野边缘与等中心间最短距离d密切相关(图 2)。同一角度、相同射野大小下,射野边缘与等中心间最短距离越小,相对误差越大,射野20 cm×20 cm (Y1=-18,Y2=2)、楔形角60°最大相对误差高达-4.7%,楔形角45°最大相对误差为-2.9%。楔形角30°非对称射野最大相对误差绝对值为-1.9%。以上特性同样存在于能量为10 MV (楔形方向为Y1-in),楔形角60°非对称野最大相对误差出现在射野20 cm×20 cm (Y1=-18, Y2=2),为-3%。楔形角30°所测楔形野最大相对误差绝对值为1.1%。能量10 MV相对误差比6 MV相对误差要小。
如2.1所述,得到如下校正方法:对由公式(1)得出的楔形角60°,45°楔形因子计算结果做如下校正,如公式(4),a为校正因子。楔形角30°及以下角度不作任何校正。以20 cm×20 cm、10 cm×10 cm为例,结果如表 2所示。经数据分析表明,以上校正方法同样适用于楔形方向Y2-out。
$ EDW{F_\theta } = \frac{{ST{T_\theta }(Y + \alpha )}}{{ST{T_\theta }({Y_f})}} = \frac{{{w_\theta }[ST{T_G}(Y + \alpha )] + {w_0}[ST{T_G}(0)]}}{{{w_\theta }[ST{T_G}({Y_f})] + {w_0}[ST{T_G}(0)]}} $ | (4) |
$ \left\{ \begin{array}{l} 射野几何中心在等中心右边(包括对称野)d≤4 {\rm{cm}}\left\{ \begin{array}{l} 楔形方向Y≥15 {\rm{cm}}时,a=0.3\\ 楔形方向Y < 15 {\rm{cm}}时,a=0.2 \end{array} \right.\\ 射野几何中心在等中心左边)d≤4 {\rm{cm}}\left\{ \begin{array}{l} 楔形方向Y≥15 {\rm{cm}}时,a=-0.3\\ 楔形方向 Y <15 {\rm{cm}}时,a=-0.2 \end{array} \right.\\ 其他情况不作校正 \end{array} \right. $ |
以楔形角60°为例,能量6 MV、10 MV测量结果与Eclipse计划系统计算跳数MU结果最大相对误差分别为1.4%、1.1%。校正前,与手工计算所得跳数MU最大相对误差分别为-9.9%、-8.4%,校正后相对误差降为-6.9%、-5.8%。相对误差随楔形角度的增大而增大,且与射野大小、射野边缘与等中心间最短距离有密切关系。此外,能量10 MV楔形因子、跳数MU与测量结果相对误差比能量6 MV对应的相对误差要小。给予相同处方剂量200 cGy,手工计算和Eclipse计划系统计算所得MU与测量结果对比分析,以楔形角60°、45°,射野20 cm×20 cm、10 cm×10 cm为例(表 3)。
根据加速器输出分割模型的定义,将动态楔形野剂量贡献分为两部分:第一部分为准直器从开始位置运动到测量点位置的贡献,第二部分为准直器从测量点到最终位置的贡献。第一部分没有考虑到准直器运动过程中射野面积的减少所导致的散射贡献减少,造成第一部分“高估”测量点的剂量。第二部分因测量点被准直器完全遮挡,模型定义第二部分对测量点没有贡献,忽略了准直器的穿透射线和侧向的模体散射,造成第二部分“低估”测量点的剂量。因此,对于非对称野来说,射野几何中心在射野等中心(测量点)左侧时,计算值大于测量值,反之小于测量值。射野越大且不对称程度越高的情况下,EDW楔形因子计算结果与测量结果相对误差越大。据此,在本研究中,当射野几何中心偏离射野等中心大于4 cm时,按照楔形野运动方向(Y方向)分为大于和小于15 cm两种情况分别给予合适的常数因子校正。
对于射野中心点处EDW楔形因子的研究,有报道表明将对称野的楔形因子应用于计算非对称野的处方剂量时,其计算结果与测量结果误差最高可达48% [12]。Kuperman [4]在加速器输出分割模型上提出的漂移方法(shift method)计算射野中心点楔形因子,相对误差在1.0%~1.8%,方法计算量庞大,且理解较困难。Gossman [13]在上述漂移方法基础上提出矩阵校正方法,与瓦里安Eclipse计划系统计算结果对比,最大相对误差小于2%。校正矩阵因能量而异。本研究结果表明:经该常数校正因子校正后,楔形野最大相对误差减少到2%以内,方法简单有效,且同时适合6 MV和10 MV光子线能量。
在上述研究中,还提供了动态楔形野下对射野中心点的常规跳数计算和测量的结果。显然,对于对称野来说,经过楔形因子的校正后,手工计算结果相对测量结果的误差在1.6%内,符合临床应用的要求。对于非对称野,当射野边缘与等中心间最短距离d≥4 cm时,手工计算结果与测量结果相对误差在3%以内;但当d < 4 cm时,手工计算结果与测量结果最大相对误差超过5%。Eclipse系统计算结果与测量结果相对误差小于2%。显然,对于这种非对称射野来说,不能简单采用针对射野中心点的常规跳数计算方法,考虑采用较为复杂的非对称射野处方剂量计算方法[4, 14],或采用合适的治疗计划系统计算[15],例如Eclipse治疗计划系统。当然,采用直接测量的方式给予相应的剂量跳数也是方法之一。
[1] | Edlund TL. Treatment planning of oblique wedge fields comparing enhanced dynamic wedge and standard 60 degree wedge for parotid type treatments[J]. Med Dosim,1997, 22 (3) : 197-9. DOI: 10.1016/S0958-3947(97)00017-4. |
[2] | Warlick WB, O'rear JH, Earley L, et al. Dose to the contralateral breast: a comparison of two techniques using the enhanced dynamic wedge versus a standard wedge[J]. Med Dosim,1997, 22 (3) : 185-91. DOI: 10.1016/S0958-3947(97)00015-0. |
[3] | Prado KL, Kirsner SM, Kudchadker RJ, et al. Enhanced dynamic wedge factors at off-axis points in asymmetric fields[J]. J Appl Clin Med Phys,2003, 4 (1) : 75-84. DOI: 10.1120/1.1534710. |
[4] | Kuperman VY. Analytical representation for varian EDW factors at off-center points[J]. Med Phys,2005, 32 (5) : 1256-61. DOI: 10.1118/1.1872532. |
[5] | Liu C, Li Z, Palta JR. Characterizing output for the Varian enhanced dynamic wedge field[J]. Med Phys,1998, 25 (1) : 64-70. DOI: 10.1118/1.598161. |
[6] | Gibbons JP. Calculation of enhanced dynamic wedge factors for symmetric and asymmetric photon fields[J]. Med Phys,1998, 25 (8) : 1411-8. DOI: 10.1118/1.598313. |
[7] | Yu MK. Analytical representation of enhanced dynamic wedge factors for symmetric and asymmetric photon fields[J]. Med Phys,2002, 29 (11) : 2606-10. DOI: 10.1118/1.1513568. |
[8] | Liu C, Kim S, Kahler DL, et al. Generalized monitor unit calculation for the Varian enhanced dynamic wedge field[J]. Med Phys,2003, 30 (7) : 1891-6. DOI: 10.1118/1.1586269. |
[9] | Varian Medical System, Inc. C-Series clinac: Enhanced dynamic wedge implementation guide[M]. Palo Alto (CA): Oncology System/Delivery Systems, 2006 . |
[10] | Petti PL, Siddon RL. Effective wedge angles with a Universal wedge[J]. Phys Med Biol,1985, 30 (9) : 985-91. DOI: 10.1088/0031-9155/30/9/010. |
[11] | Khan F. The physics of radiation therapy[M]. Philadelphia: Lippincott Williams & Wilkins, 1994 : 176 -225. |
[12] | 王丹, 何京学, 高鸿, 等. 手工计算非对称动态楔形野处方剂量的校正[J]. 中华放射医学与防护杂志,2010, 30 (1) : 47-9. |
[13] | Gossman MS, Sharma SC. Analytical correction of an extension of the "MU Fraction Approximation" for Varian enhanced dynamic wedges[J]. J Med Phys,2010, 35 (2) : 65-72. DOI: 10.4103/0971-6203.62195. |
[14] | Thomas SJ, Thomas RL. A beam generatin algorithm for linear accelerator with Independent collimators[J]. Phys Med Biol,1990, 35 (9) : 322-5. |
[15] | 邢晓汾, 崔桐, 郑旭亮, 等. 两种治疗计划系统三种算法对EDW模型准确性比较[J]. 中华放射肿瘤学杂志,2012, 21 (5) : 468-70. |