文章快速检索     高级检索
  南方医科大学学报  2016, Vol. 35Issue (2): 260-263  DOI: 10.3969/j.issn.1673-4254.2015.02.21.
0

引用本文 [复制中英文]

黄思娟, 陈立新, 曹午飞, 孙文钊, 陈阿龙, 刘伯基, 王彬. 扩充型动态楔形板楔形因子的校正方法及跳数计算[J]. 南方医科大学学报, 2016, 35(2): 260-263. DOI: 10.3969/j.issn.1673-4254.2015.02.21.
[复制中文]
HUANG Sijuan, CHEN Lixin, CAO Wufei, SUN Wenzhao, CHEN Along, LIU Boji, WANG Bin. Correction of enhanced dynamic wedge factor and analysis of monitor unit calculation[J]. Journal of Southern Medical University, 2016, 35(2): 260-263. DOI: 10.3969/j.issn.1673-4254.2015.02.21.
[复制英文]

基金项目

国家自然科学基金(81272486);广东省教育部产学研重点项目(2012B091000144)

作者简介

黄思娟,硕士,物理师,电话:020-87342255,E-mail:huangsj@sysucc.org.cn

通信作者

陈立新,博士,副研究员,电话:020-87343089,E-mail:chenlx@sysucc.org.cn

文章历史

收稿日期:2014-08-22
扩充型动态楔形板楔形因子的校正方法及跳数计算
黄思娟, 陈立新, 曹午飞, 孙文钊, 陈阿龙, 刘伯基, 王彬     
中山大学肿瘤防治中心//华南肿瘤学国家重点实验室//肿瘤医学协同创新中心,广东 广州 510060
摘要: 目的 研究瓦里安扩充型动态楔形板楔形因子计算修正方法,比较楔形野中心点处手工和治疗计划系统计算结果相对测量结果的剂量/跳数差异。 方法 对于瓦里安直线加速器的6 MV、10 MV光子线,使用指形电离室测量水下10 cm处不同动态楔形野的楔形因子及射野中心点的剂量,采用治疗计划系统计算相应射野的剂量/跳数。使用加速器输出分割模型手工计算射野的楔形因子,并采用常数因子修正手工计算结果。对手工计算、治疗计划系统计算和测量结果进行比较,分析三种方法下常规二维治疗计划下动态楔形野的楔形因子和射野中心点跳数的误差。 结果 以测量结果为标准,校正后,手工计算的楔形因子误差明显减小。其中,6 MV光子线下,60°楔形角下对称野最大相对误差由4.2%减小到1.3%,非对称野最大相对误差由-4.7%减小到-1.8%。10 MV所有楔形野相对误差由最大-3.0%降低到1.1%。手工计算跳数与测量结果对比,对称野相应射野跳数计算相对误差在2%以内,但部分非对称野最大相对误差超过5%。比较治疗计划系统计算结果与测量结果,其最大相对误差小于1.5%。 结论 使用常数因子可以有效减小输出分割模型计算楔形因子的误差。对于常规二维治疗计划楔形野的跳数计算来说,校正后对称野射野中心点的计算结果符合临床治疗要求,但对于射野边缘与等中心最短距离小于4 cm的非对称野来说,需要使用相应的非对称射野处方剂量计算方法, 或者采用测量方法或利用治疗计划系统计算相应的射野跳数。
关键词: 扩充型动态楔形板    楔形因子    校正    跳数    剂量    
Correction of enhanced dynamic wedge factor and analysis of monitor unit calculation
HUANG Sijuan, CHEN Lixin, CAO Wufei, SUN Wenzhao, CHEN Along, LIU Boji, WANG Bin     
State Key Laboratory of Oncology in Southern China, Department of Radiation Oncology, Cancer Center, Sun Yat-sen University, Guangzhou 510060, China
Supported by National Natural Science Foundation of China (81272486)
Abstract: Objective To study the correction of algorithm for Varian enhanced dynamic wedge (EDW) factors and compare the dose/monitor unit (MU) deviation measured at the central axis of EDW field with that obtained by manual calculation or using the treatment planning system. Methods EDW factors and dose were measured with Thimble ion chamber at 10 cm depth under the water for 6 MV and 10 MV photon on Varian linear accelerator. The corresponding calculations were done with the radiation treatment planning system. An analytic formula, namely the"MU Fraction model", was used to calculate the EDW factor, which was corrected with a constant factor. The MU of conventional 2-D planning derived from manual calculating, treatment planning system, and actual measurements were compared. Results With the measured results as the standard, the corrected manual calculation deviation of EDW factors was significantly reduced. For photon 6 MV, the maximum deviation reduced from 4.2% to 1.3% for 60° symmetry fields was, and from -4.7% to -1.8% for asymmetric fields. For photon 10 MV, the maximum deviation for all EDW fields was reduced from -3.0% to 1.1%. Comparison of the manual calculations with the measured results showed a MU deviation for symmetric fields within 2%, and more than 5% for some asymmetric fields. The deviation between the calculations of the treatment planning and the measured results was less than 1.5%. Conclusion Constant factor correction can effectively reduce the deviation of manual calculation. For MU calculation of EDW field in conventional 2-D dimensional treatment planning, the corrected results of symmetric fields meet clinical requirements. While the minimum distance between the field edge and the central axis was less than 4 cm in asymmetric fields, the corresponding special method, measurement or the treatment planning system should be used to calculate the dose/MU.
Key words: enhanced dynamic wedge    Wedge factors    correction    monitor unit    dose    

与传统物理楔形板不同,扩充型动态楔形板(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)。其中θ为楔形角,w60w0为构成楔形角θ所需的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)
1.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 公式(3)计算参数说明 Table 1 Parameter description of equation 3
1.4 数据分析与修正方法

计算公式(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相对误差要小。

图 1 对称野EDWF相对误差分析 Figure 1 Deviation analysis of symmetric fields. [Deviation=(EDWFMea-EDWFCal)/EDWFMea].
图 2 非对称野EDWF相对误差绝对值随角度、射野大小变化规律 Figure 2 Correlation of asymmetric fields EDWF absolute deviation with angle, field sizes.
2.2 分析楔形因子相对误差与楔形角度、射野大小、射野边缘与等中心间最短距离d之间的关系

如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. $
表 2 校正前后EDW楔形因子计算值与测量结果的相对误差 Table 2 EDWF deviations between manually calculated (uncorrected/corrected) and measured results (%)
2.3 校正前后手工计算跳数MU与测量结果相对误差

以楔形角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)。

表 3 校正前后手工计算跳数MU与测量结果相对误差 Table 3 Mu Deviations between manually calculated (uncorrected/corrected) and measured results (%)
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.