比较IMRT、RapidArc、和Tomotherapy之间的计划设计时间、出束时间和计划质量 A3 p3 j2 {5 V0 K5 @( _$ s# A$ s
) M) I P8 m& T- w8 i' j! o4 s
作者:Mike Oliver*、Will Ansbacher、WayneA.Beckha- z8 k/ x' t4 p5 N! C
7 @' E1 G8 { H; y; Y! s8 j5 B
英国哥伦比亚肿瘤研究所医学物理中心 维多利亚,英国属哥伦比亚,加拿大- U9 {- ~/ v* ]5 B1 R% A- a+ F
# O. m; Q9 ~ E' ^9 A
Eail:moliver3@bccancer.bc.ca
6 s+ J* O( B9 T1 ~2 G; H' M
! `8 `; o3 {) R 这个研究的目的是比较5野和9野IMRT、单弧和双弧Rapid Arc以及Tomotherapy在的计划质量、计划设计时间和治疗出束时间估计值三者之间的差异。对四个评估模体,分别建立了5野和9野IMRT、单弧和双弧Rapid Arc和Tomotherapy计划。计划评估基于是否能满足剂量-体积约束、剂量均匀性指数、放疗适形度指数、计划时间、出束时间估计值、整体剂量和受到大于2和5Gy剂量的体积等因素。对所有的评估模体,在计划时,Tomotherapy均能满足优化的限制要求(50%P1、67%P2、0%P3、50%P4),Rapid Arc所能达到优化限制的程度较差一些(25%P1、17%P2、0%P3、0%P4),而IMRT是根本达不到任何约束限制。Tomotherapy计划剂量分布最均匀。Tomotheraoy计划具有较长的计划时间、治疗出束时间估计值也较长、适形度指数较低和累积剂量较高等特点。Tomotherapy计划能产生高质量计划,在轴面上能形成比IMRT或RapidArc好的剂量分布,却在靶区体积上下方向出现了剂量增高的现象。当然,在这个研究检查测试结果中能看到Rapid Arc可以做出比IMRT产生更好的计划。1 A. a7 I3 L. M
8 S/ f4 c, e6 ]1 J5 H! f9 j$ H
PACS编号:87.55.x,87.55.D,87.55.de,87.55.dk
* ?# }; }. W' Z4 `9 {0 N3 w5 S* u0 P6 c0 Q4 n* D
关键词:IMRT、RapidArc、tomotherapy,治疗计划。
& Z; Z4 S: @' X9 m: A" l) ?/ ~2 p& G. O! I
I、介绍
: n$ q, O- U6 m: ]' C$ r$ z: k& t$ Z% L3 y8 ^8 N9 ~
IMRT、RapidArc、和Tomotherapy都是最前沿外照射治疗技术,在过去10年不同阶段中已经都应用于常规临床治疗。使用配有MLC的常规直线加速器IMRT在1995年投入临床治疗前列腺癌(尽管使用补偿器的IMRT出现时间更早),随后这项技术应用于治疗其他的解剖部位。最近一个基于56份临床实验数据的META-分析显示IMRT相比于非调强治疗可以降低损伤;但该分析没有给出关于局部控制和存活率的数据。从这些数据证明IMRT优点的方法中可以得出,未来前沿放射治疗技术应该在于计划质量、效率、图像指引的准确性和提供任意的剂量分布等方面具有优势。这篇文章研究放射治疗在质量和出束效率上的提高。1 T! P4 E! r6 A, q
9 N7 t( K- r6 y: f
Yu 首次提出在治疗过程直线加速器机架旋转执行中并伴随MLC运动的IMRAT技术。旋转放疗技术的临床应用还很少,但也已应用于治疗中枢神经系统、前列腺、头颈部、全腹全盆腔治疗、直肠癌、子宫内膜癌 。直到Otto提出他的容积调强弧形治疗算法后,IMAT最主要的优势才被发挥出来。VMAT使用了进阶取样算法,先从粗略的机架角度取样开始,在优化的过程中,弧度的分辨率可以得到很大的提高。如果没有这个算法,相邻MLC叶片受到叶片运动能力的极大限制。VMAT通过在优化开始时允许大幅的叶片运动,再在其后的过程中逐步限制叶片运动的方法克服了这种限制。优化时间也大幅降低。Otto的算法已被Varian应用于实践,商品名称为Rapid Arc。在这个应用中,进阶取样通过5级离散“多分辨率”把子野数由10个增加到177个。Elekta也有一款产品命名为VMAT但没有采用Otto的算法而是采用自己的一套算法。 Varian 和Elekta的弧形治疗都允许剂量率在治疗中变化。 ; t$ D. k; |8 o! e+ i
& l+ V3 t8 k" c
螺旋断层放疗由Mackie首次提出,现在商品名为Tomotherapy。Tomotherapy是一种似于常规CT的技术,它在使病人连续通过断层放疗设备断层时发出围绕病人旋转扇形光子束进行治疗。射线以经过二元的MLC调制后的沿螺旋轨迹运动。治疗在51个投影角度上进行优化,概念上可认为是51个等间隔机架角的IMRT。螺旋断层放疗在2002年治疗了第一例病人,到目前为止,已在全世界广泛使用。
0 x: I( k' n* }% Z8 v9 M
6 E3 D1 S# f- Q3 G# @( R 现在需有必要总结IMRT、现代旋转放疗、断层放疗在计划质量、计划时间和出束时间估计值等方面的特性。先前的研究已经比较IMRT和Tomotherapy、IMRT和旋转放疗、Tomotherapy和旋转放疗。有一个简单研究是对5个病人良性颅内损伤比较IMRT、旋转放疗和Tomotherapy计划,这项研究结论为所有技术实际上是等效的,该研究还建议进一步使用一些具有挑战性的病例做比较,也为该研究提供一些动力。此外,Brotfeld和Webb最近发表了一份基于二维模体解析IMRT、单弧IMRT 和Tomotherapy所能达到剂量分布的理论探讨。他们的结论是出束少于2分钟的单弧放疗对非常复杂病例在计划质量上可能稍欠缺,认为IMRT和单弧IMRT的计划质量应当是类似的。 5 G& G7 }) Q6 |4 q( F
* I* v; i5 t7 U+ O6 d0 i( f# Z' b6 _
本研究的目的是在四个虚拟模体(含Brahme模体)去完成IMRT、RapidArc、和Tomo治疗计划,同时去比较它们之间的治疗计划质量、优化时间和出束时间。
; X9 u, ^+ A: K* Q. |$ Y+ Q) U( O0 x- M4 ^3 U2 j& h( B* Y
II、材料和方法8 D* Y! ?. }/ K
7 ]' E$ P1 H8 V1 B A、方法概述
" c( s4 f0 L' P6 g1 z0 v1 {( G* c3 D6 b
本文的方法摘要在图1概括,其中涉及到了不同的章节、主要的计划系统和本研究用到的剂量评估软件。系统间数据传递在图1中以输入和输出箭头表示。
0 w/ y! A5 \% {; P, K0 |5 o) Z! C! f% i8 J9 M: d0 Q- J
8 D5 t0 V5 f# {$ _4 {* M. C0 C& W* a$ B# R3 }: F" r
$ T* R7 k$ U5 L B/ c- c) L
! k$ _2 O K% ~6 w% @: O
+ M- U# Z3 j. j; n- o+ G1 n9 R0 v- y* D) D2 Z
3 X% N4 i- O! b& g- }- e9 k+ K
2 ^$ d Q% z3 f1 }! o+ r9 H图1黑色图框代表本研究采用到不同的主要软件,包括简要的方法步骤(后面会进一步解释)数据的传递在模块中用箭头表示。2 r8 d" a6 C+ G# r
; V5 j. D- M+ R7 l/ z7 B
B、模体、轮廓和计划目标的介绍
; m: X/ P. f& k
) ~; p' I" a( K# z所有虚拟模体和轮廓设置都是在Eclipse治疗计划系统中建立。所有模体都是水等效模体,10cm半径和长度为25.25cm。图2提供了所有模体包括轮廓的横断面视图。每个模体包括外轮廓、PTV和一些危及器官(OAR)。在三维方向上通过外放使外轮廓扩展。PTV和OAR前后长度如下:模体1 TV长度为8.5cm,OAR1长为6cm;模体2和3:PTV、OAR1、OAR2、OAR3长为5.25cm;模体4:PTV长为8cm,OAR长为7.5cm。PTV和OAR根据它们在临床代表的能被看到的或者在实验中能较好评估的解剖组织确定其位置。模体1类似于美国放射学院测试实验。模体2代表前列腺组织,模体3代表头颈部结构。最后,模体4是Bortfeld在分析中用过的实验。# x7 V9 C) v7 b j8 o1 I
( a4 m$ L! i1 M, F
表格1中列出所有模体的计划目标函数。这些治疗计划的处方剂量是60Gy/30f。所选用的剂量-体积约束对三种技术都具有挑战性。约束的权重由各自的计划系统确定。
8 O! `1 J; b6 B
0 @' b% Y$ ^6 x) g8 v' W# t" o1 c表1 在模体1-4优化过程中多用到的所有剂量体积约束列表0 W @& h4 r! r! d( s
2 b4 i- K3 ^" ^( f
PTV OAR1 OAR2 OAR3
' G9 A& V& c8 U/ F, Q* b( l; b; Z9 v* e$ y- R; Z
Min Max DVH DVH Max DVH DVH Max DVH DVH Max
" o2 T+ y# }/ n" i2 Q, e' B* [' d. h; p6 m" r
Dose Dose Vol Dose Dose Vol Dose Dose Vol Dose Dose e- w6 D! v1 ]) a0 E
$ q* i2 T+ t9 Y& h0 k2 b( T* j/ ~
(Gy) (Gy) (%) (Gy)(Gy)(%) (Gy)(Gy) (%) (Gy)(Gy)
6 @# c0 X" l/ `- U, K; m; G2 e" }! L* ^! l/ N4 S6 R
Phantom 1 60.0 65 50.0 13.3 20.0
/ n' w4 W- }) S* x4 K/ |+ H# U5 G9 W6 b' `* w$ p
Phantom 2 60.0 65 50.0 10.0 20.0 50.0 10.0 20.0/ j0 L5 U4 |, M
" h: r0 @9 L9 o- [3 }Phantom 3 60.0 65 50.0 10.0 20.0 50.0 10.0 20.0 50.0 15.0 30.0
* j I& k9 L& G5 S4 v4 S8 g; t% v9 H2 k7 T- E
Phantom 4 60.0 65 50.0 15.0 30.0( ]8 A0 ?9 O5 s( T+ b7 e
: k( \% @1 n! m- R& T( }$ w( ~+ A. z
C. Sliding window IMRT计划
" B6 g; F/ w( q0 b* C( H0 Q& I# L5 b4 E/ a! T
Sliding window IMRT治疗计划设计是在Eclipse系统中进行。两个计划用等中心为5mm叶宽的120叶设计。5或9野的选择取决于治疗时间和计划质量之间的一个平衡。文献表明,超过7-9野,射野方向的选择对共面IMRT影响不大。
3 L! s, }7 E9 F/ Q- w$ |
, i5 L: c: q: }, _% o& X# q# B6 D 在计算得出最后剂量之前,优化在beamlet模式下要经过大约250次的反复运算。在250次的反复运算之后,所有计划成本函数最后汇聚在一起。最后剂量矩阵通过解析析各项异性运算法则(AAA)计算,采用体素大小为0.25x0.25x0.25cm3。Eclipse的CPU是双核Intel Xeon qual-core处理器运行2.50GHz。+ x1 b: O$ b# I" S. ~4 {8 S
8 b8 z( |, s4 X9 N" N3 @! z, S D、RapidArc 计划' o! @: j1 p: i% E4 ]) \% ^# t k4 ?* g
# X$ V8 \9 G# j5 I RapidArc计划同样也是在Eclipse系统中设计的,该系统也支持多弧计划设计。单弧旋转RapidArc(RA1)和双弧旋转RapidArc(RA2)的治疗计划,其准直器角度均为45度。RapidArc根据计划质量和治疗时间之间的平衡选用单弧或双弧。虽然单弧计划出束时间较短,但增加一个额外弧在提高计划质量同时也延长了治疗时间。计划中运用了一项使所有弧形子野的MLC初始化位置为PTV减去OAR位置的技术。在这个版本的RapidArc计划系统中,在目标函数达到预设水平之前,优化可能就转换到下一个MR水平。转换可以被用户延迟,完成它需要从MR1一分钟上升到MR5的五分钟。若弧形优化完成,优化的过程被保存,若达到MR5水平优化要继续额外的五分钟。最后剂量的计算通过AAA计算,体素的大小为0.25x0.25x0.25cm3。Eclipse系统的硬件要求同前段描述一样。( f0 B! g6 }& L# F
0 P& N/ M( m) }+ P/ Q6 s
E、Tomotherapy计划" \3 P) }) G9 K+ e4 b/ Y
! M- w+ A0 T9 Z5 f
Tomtherapy计划是在3.1.2.3版本Tomotherapy计划工作站完成的。在Tomotherapy计划中,射野宽度定义为扇形射束轴向的厚度。螺矩定义为机架完成一次旋转床移动的距离,也称为旋转轴方向上的轴断面射野宽度。调制因子定义为最大叶片开放时间/平均叶片开放时间。选择可以降低影响最后剂量分布的threat效应的螺矩因子,通常为0.86/n,其中,n为整数。7 B" ~! g3 i, |' T( A
7 t8 T6 R% |9 K9 j 治疗计划设计过程先进行全部子野的剂量计算,然后再执行经过250次子野优化循环。经过250次反复运算后,所有计划的成本函数趋向一致。Tomo采用的参数如下:等中心扇形野宽度为2.5cm,螺矩为0.287,调制因子为2.5。Tomo治疗计划的剂量体积元大小为0.31x0.31x0.25cm*3。Tomotherapy阵列计算机采用16 个Intel 双核处理器2.4GHz。
" ]' |( y" @# P+ D
% q& n G( P9 s% p2 G) N b3 Y2 HF、计划质量是
! y( ]0 Q6 Y' J1 o% H6 z# v+ o
. ~- j0 [* F, @2 P 在本研究中,F1提及的计划质量主要指标是在优化过程中是否满足剂量体积约束的二进制量。是否满足DVH限制值主要看适形度和均匀性,OAR的平均剂量,具体算法在F1、F2节中单独提及。
, \0 ]2 @- j% o" R: s* ~4 Q8 d8 L7 I
F、1剂量体积评估. m( o4 c/ d+ M2 g5 l. l
0 V M- i8 R x' h& X: n, l 优化完成后,为进行剂量学评估,所有计划导入CERR3.1软件。所以计划被优化到PTV95%体积受到61Gy剂量。& b9 r/ x# h1 L G) K# z! r1 r
& z% |( O0 g. `; f" M
报告的最小剂量为99.5%PTV体积的剂量,而报告的最大剂量为0.5%PTV体积的剂量。表1列出的剂量-体积限值,由接受剂量D的百分体积(VD)和最小百分体积剂量V(DV)报告。图3图示了剂量-体积约束是如何报告的,治疗计划用能否满足DVH限值的二进制量评估。! T: H" l; X6 |) t) T$ R
& s, ?: L3 F* [* s: p9 z3 u: BF、2 适形度和均匀性指数
& I5 Y0 g+ G' N7 H1 l6 H
8 M- S1 ~4 L5 b9 D ~9 O 治疗计划的放射适形度指数(RCI)和剂量均匀性指数(DHI)分别用公式1和公式2计算。RCI描述PTV受到处方剂量或更高剂量的体积(VD,,PTV)与外轮廓内受到处方剂量或更高剂量的体积(VD,body)之比,提供一个确定在PTV和正常组织之间处方等剂量体积如何分布的定量方法。5 _( w6 e0 O+ S( f
9 F. O5 z7 i2 R6 j
VD,,PTV
+ ?. X8 Y' G) d% V- C; P9 L0 N3 Q0 S1 t
RCI= VD,body (1)
4 {, e% C. N# s% M$ I: x- j* N
/ {1 o G& r9 E/ W+ ^( A! H6 v DHI描述计划靶区内剂量的均匀性,是最小剂量(D99.5%)和最大剂量(D0.5%)之间的比值。
8 z0 g, J1 S% B# x8 \/ Z. ^2 y+ y8 O7 C7 _
D99.5%
1 D f3 c! u8 E% _1 z+ [2 M$ ]# a/ f+ F$ ~* ^( i
DHI= D0.5% (2)# [. x9 H2 d2 P9 j
% R! z6 D1 r8 K3 u3 M% mF、3 危及器官的平均剂量0 T3 p% D0 I- ~! C
0 G, c. L+ _7 m a3 W6 I
每个危及器官(OAR)的平均剂量都做了记录。要求的记录点包括50%剂量的体积和OAR最大剂量。因此,平均剂量应与要求的剂量-体积点相关。7 [5 k( ]3 I& R. Q+ g& P( I" m
: J8 Q) D3 [* q6 h) r8 Z
G、累积剂量
( M! X# k; M- q" o
- h! q1 ]2 b, H) d) x) X 如公式3所示,累积剂量是所有剂量体积元的质量加权和。模体有均匀的密度,所以公式可以做如下简化:N是体积元个数,Dmean是外轮廓平均剂量,mvixel为体积元的质量。公式简化为模体外轮廓质量(mbody)和模体外轮廓的平均剂量(Dmean,body)之间的乘积。2 D% g# l9 E9 Z3 t* w- c$ J
8 M3 s( C0 j! z3 v3 t6 b Eintegral=Di*mi=N*Dmean*mvixel=mbody* Dmean,body7 G: g5 g6 c( p. M
5 c% T+ r1 l) A9 M0 wH、受到>2Gy和>5Gy剂量的体积% U( c! Y( u) W
9 h! E- H" _" O4 j8 l
一些放射致癌变模型认为剂量反应关系直到6Gy剂量都呈线性关系,然后才达到高峰坪区。受到>2Gy和>5Gy剂量的体积在本文中都很重要,每个治疗计划中都将涉及到。
1 t- E# A& {+ e7 `0 y% M2 L. C! i% K
I. 治疗计划时间
2 m! @: D9 U7 a: M5 D# x6 H" ^+ y
治疗计划时间定义为从开始计划设计直到最后优化和剂量计算完成之间的时间。由于假设计划能用于临床,因此也包括相关参数的选择、治疗的优化和最终剂量计算的过程。" s) j: t2 ~1 w% K: z
) C, G, H& `' Z+ z& h
J、治疗出束时间的估计) M/ R L$ M, i/ w
; I. O; ]' g7 C+ h) L- D
治疗出束时间定义为首次射线开启到最后射线关闭之间的时间。对sliding windows IMRT来说,就是剂量率乘以每野跳数,加上机架在连续野之间旋转的时间,再加上参数“delta”(包括考虑MLC出束文件的数据传递的时间、模式准备、旋转时间的错误估计和操作员的反应时间)。第一个射野不考虑delta参数,因为数据的传递和模式准备发生在第一个射野出束之前。delta参数(在我们中心治疗的215个头颈IMRT野提出)定义为剂量率为400MU/min和机架旋转速度为3600/min条件下的理论出束时间和被记录和验证系统记录的实际出束时间之间的不同。对RapidArc来说,治疗出束时间估计值是在一个弧度内177个子野出束的时间。这个时间是所有子野角度增量的和除以通过机架旋转速率。对多弧计划来说,只对第二个弧加delta参数,因为所有数据传递发生在射线开启之前。TOMO放疗的出束时间估计值是通过计划系统和来自最后计划报告的记录进行计算的。
' c3 K, k9 m2 a6 ^
) |( i- P% A+ O: O6 QIII.结果6 _7 a) T' |* J0 V4 Y% P
1 O' t$ C" \3 Y- xA. 剂量-体积评估
- X) Y% H/ f! `1 b
+ r. D2 H, _5 K9 K; d; E5 C在表2中,执行了剂量体积评估,列出了所有治疗计划的数据。图4列出了9野IMRT、2单弧RapidArc和Tomotherapy的DVH曲线。图5,列出了9野IMRT、双弧RapidArc和Tomotherapy剂量分布的轴断面等剂量曲线。基于表2数据,模体1,IMRT5和IMRT9,0/4符合优化约束;RA1和RA2,1/4符合约束;TOMO,2/4符合;模体2,IMRT5、IMRT9和RA1,0/6符合约束;RA2,1/6符合,TOMO,4/6符合;模体3所有的计划只有0/8符合约束;最后,对于模体4,TOMO有2/4符合约束,而其它计划有0/4符合约束。+ X2 o, [9 v* a5 _8 c$ c
2 k5 e% a/ [* _, x- t- L8 [
表2是 所有剂量和体积参数的一份总结,列出实际上达到的剂量和体积。作为参考,要求剂量和体积值在应达到值的上面列出,以粗体和斜体标出的值是预先规定优化限制值。
) k2 F4 x8 S$ z2 m1 L) T6 ]) F T+ {" B; ? @1 w- _
PTV OAR1 OAR2 OAR3 + P4 W( \6 X" R3 f1 R% Q
& c: W: E+ D- H( X Min Max DVH DVH Max DVH DVH Max DVH DVH Max
( v! [8 c& e7 L- U
* g+ h- C4 N" z& q; A Dose Dose Vol Dose Dose Vol Dose Dose Vol Dose Dose+ d- }( K2 d1 f$ y7 D
% d# T6 _. p( X5 ?9 ^# Q Phantom 1 60.0 65 50.0 13.3 20.07 a' v& A+ W" b
1 ]9 ?3 N5 k$ B2 l: ~4 j$ S" j; dIMRT 5野 59.2 68.9 82.1 16.0 29.3% `, K3 e1 M$ T0 F- N
# i) N8 L) R) p0 W 9野 59.3 73.1 87.1 16.4 30.3
8 r6 p+ v7 H# l, R; P( `9 h, f( u* q+ l0 [
RapidArc 1弧 59.2 68.7 40.9 12.4 23.6
6 A" E- B: W9 n; |$ e6 |7 B5 f: d% L
2弧 59.5 67.0 40.1 12.5 22.6
# u3 L. n4 x( ]0 S g2 k5 e2 c& _" K9 }; v9 ]% B; q8 X; m# @
TOMO 59.1 67.2 15.1 9.4 17.0
5 T8 a) o: J2 d5 q) K; B
. x3 Q, m: D. j" Z- \1 _Phantom 2 60.0 65 50.0 10.0 20.0 50.0 10.0 20.0, c; B V. y1 |2 c' S) f) @, P' O
6 y& i$ \0 S+ ?1 G5 ~ X
IMRT 5野 58.2 69.3 97.7 13.2 28.6 63.8 25.3 27.6$ ^ f1 O# N( x Y
6 E7 D( x7 ?. @2 B; d" [1 `+ `" { 9野 57.0 69.2 90.4 13.1 31.6 64.8 25.7 29.55 [- g$ F& j! V# X2 R# a* e
& |; ]7 Y7 N& E4 _3 `( f C6 _2 rRapidArc 1弧 58.1 68.7 58.5 10.3 23.4 52.3 10.2 24.1 & M- B+ i! |' P7 G3 ?( g, G
+ h! C0 I P/ y5 g' Y9 K4 c% D 2弧 58.3 68.0 50.2 10.0 23.1 47.3 9.8 22.9
" n* w) ^" D9 R8 G; m0 w
8 I) Q/ D3 P' m# mTOMO 59.9 66.1 42.6 9.5 17.4 33.9 8.8 17.9
6 W) r+ W. h [1 \2 |. y. R
% ]. k. T( r V) @8 n. pPhantom 3 60.0 65 50.0 10.0 20.0 50.0 10.0 20.0 50.0 15.0 30.0
& d% Q% y$ j# D( Z
# l* K' I# C" R, h) C6 c5 K; PIMRT 5野 58.8 70.4 79.6 12.9 28.4 66.1 13.3 35.4 66.9 19.5 44.7
, }* S' D$ \( ~) _0 o
6 y, R3 f3 ]" a! ?" Y4 ]6 _ 9野 58.3 70.1 84.8 14.3 31.3 69.4 13.6 40.1 69.7 18.7 46.8 0 {2 C) R5 z3 f8 H0 o
; J9 Z2 ~0 h" U0 ]6 Q4 W+ l
RapidArc 1弧 58.3 69.0 94.9 13.4 24.7 81.9 13.1 28.6 91.3 21.7 44.22 Q, W0 h$ d- ~
" h$ b. `0 G/ b% T6 t7 j* d$ Y( x$ i
2弧 58.3 69.0 90.7 13.5 26.1 88.1 13.2 27.8 93.8 21.2 43.63 k0 k V: c3 h: r& @# f" T+ S
& j& ~9 W9 i1 O/ R7 ?! g5 p. J
TOMO 56.7 66.7 99.9 15.1 24.6 96.8 13.5 25.2 83.0 21.5 43.6
: ~0 q5 q1 b: x& x; U" V `) ~6 z7 x, n, z7 O
Phantom 4 60.0 65 50.0 15.0 30.0
( r; ?6 [, X9 p
3 l" [: c! a7 ^- w' S% F1 C. TIMRT 5野 55.0 75.3 86.5 23.3 41.4
( m3 _$ B$ }$ b
: Y& H: y/ \2 h5 l! p 9野 57.6 69.5 87.6 21.6 39.4
8 A# i$ q a9 {, y0 C: q8 [+ \% I& R9 \) i2 S8 v! b6 @. M
RapidArc 1弧 57.4 69.5 74.4 18.3 36.3
$ D" S5 m: O5 M# O- n" |2 L& C# T' C) u. U4 b6 Q* D, K
2弧 59.0 66.5 68.3 16.9 33.8
" ^1 Q# K* w; U v! H+ I$ M. N& f! r1 j
TOMO 57.4 65.6 28.7 11.2 29.5
" q# ^4 X' e1 F( x* i2 R A( z1 b2 R3 A2 u! E+ B( w
B. 放射适形度指数、剂量均匀性指数
0 u. X9 K$ T' q5 y. j5 i! H: y" Q( S' i! U
表3显示了所有计划的放射适形度和剂量均匀性指数。RCI平均值:IMRT:0.870.7,RapidArc:0.90.02,TOMO:0.810.07;DHI平均值:IMRT: 0.830.03,RapidArc:0.860.02,TOMO:0.880.02。
) l6 t# `+ \4 d: {
% a5 t* D0 L3 I& ?. h8 r; z表3 在实验模体中有IMRT、RapidArc 和TOMO计划的放射适形度指数和剂量均匀性指数在前两列列出,OAR平均剂量在三、四、五列分别列出,累积剂量(体积受照量>2Gy: L+ n) i& @# f6 A7 I* N0 n
# L Q0 v# ^: s, B# F0 d,和>5Gy)在六、七、八列分别列出。
; O/ W9 l- V/ X. U5 m$ L: @% G) _9 k7 j4 M& @
/ u% J& q( ]1 ]2 |: d$ u* q Q. X- y8 \4 A
Mean Mean Mean Integral
, S1 s% u( ?; |: @) Z+ g' E' O7 _2 h" c1 Y. Q v- J3 |4 F
RCI DHI OAR1 OAR2 OAR3 DOSE V>2Gy V>5 Gy! |. z! v: ~5 d. I6 U/ B* b, H n0 |" w
( S0 C5 T9 F1 J1 F$ y m
Phantom IMRT 5野 0.855 0.860 17.32 89.0 4148 32756 k* C; ~# E' P+ K
9 J# o# ^- c2 h( Y& o/ c
1 9野 0.913 0.811 17.78 87.8 4151 3425
* G6 {/ D, `- h$ d$ G2 M5 V+ x5 {, J6 j
RapidArc 1弧 0.890 0.862 13.65 85.8 4338 3379
6 ^" j& b; L: _- m" V% c5 v3 ~+ q1 e ], ^. y3 ^; }6 l
2弧 0.914 0.888 13.53 85.7 4311 3381
. Y& Y6 w0 r- B7 q2 e9 A2 H( _' y1 a' I. k
TOMO 0.706 0.880 10.47 106.9 4884 4281
9 v# R8 i- n7 Y( w1 z
, u2 ]% a4 J4 Q. ]Phantom IMRT 5野 0.876 0.839 14.66 13.88 59.2 2617 20517 C) L& [/ G" [) R7 Q* d+ w, ]% {
) Z/ Y C1 N4 F; P. y7 U& { 2 9野 0.921 0.824 14.90 14.28 59.0 2662 2101
( d5 I V* `- f/ w. J) S
% i( x2 N% d) RRapidArc 1弧 0.911 0.846 11.63 11.68 58.0 3019 2181
- ]8 ] y0 ~2 w% i5 E8 @ @
' _+ ?, B8 k. W0 V# p3 K: g 2弧 0.926 0.857 11.20 11.19 57.2 2970 2158
/ Y$ P* ^8 G( Q7 ^9 H; y: D
- V& n! C8 m" r A. ?! @TOMO 0.793 0.906 10.07 9.93 75.2 35.7 2995
4 S- ^! a( |! a, a1 y- ^6 J- S
6 X" f2 f* j5 u Phantom IMRT 5野 0.758 0.835 14.55 15.16 21.26 59.1 2630 2035& E# o- D. X% R9 @; B4 {
8 S' q; ?1 I: I s: h" @ 3 9野 0.783 0.832 15.86 16.17 21.46 58.3 2633 2031
4 @0 ~4 q2 e% o- i# _5 B: K; U6 B; l& [: b2 S( q5 [0 P( l3 S
RapidArc 1弧 0.904 0.844 14.44 14.56 24.27 58.6 3104 2226
3 i7 b5 W+ e- @1 c2 L+ ^. o t2 ^- ^0 w" O$ Z' d
2弧 0.911 0.846 14.43 14.71 23.81 58.2 3060 2223
# y \/ e, g& w0 }8 ~) ]$ e6 F3 o5 {0 N. w* f+ _0 L2 N
TOMO 0.847 0.850 15.74 14.77 23.23 73.5 3510 2977
$ D* O' y/ k2 F4 u7 j/ S/ `) y0 ]- Y' b [
Phantom IMRT 5野 0.769 0.730 24.34 148.3 4991 3477
& x" }1 Q. ^; K' {3 y# K; U& | e+ t; y2 G* s
4 9野 0.912 0.829 23.06 142.3 4833 3450! M: K* R) R+ g4 R: |0 q& @
: S7 z/ U& S% a0 [! PRapidArc 1弧 0.904 0.826 20.24 142.1 5562 3915
3 n8 R3 L! N @! A: n9 z6 v3 \4 s2 Q/ M$ J* H0 I/ c# d) _
2弧 0.940 0.887 18.81 141.2 5629 3915& C" e+ b# v" _4 p- o
1 J! j3 e* f7 OTOMO 0.900 0.874 13.26 158.9 5367 4195
* z z) x) [" P- Q/ r5 @% g4 V. e2 X
0 m" U2 H1 {' [. s0 D! \C. OAR平均剂量! Z# u" o) T& `& i/ h- B
3 X4 y+ @% c; s9 ] 表3列出了每个模体和治疗技术的模体平均剂量值。模体1、4的OAR平均剂量总的趋势是dOAR,TOMO<dOAE,RA2,<dOAR,RA1<dOAR,IMRT5OR9。
( k: n3 H i# K7 ~# P8 k$ p g! ?
D.累积剂量和体积(受照量>2Gy,和>5Gy)
' ]! ~, O! Y% g% F
, `) N, s g9 u8 b# P 表3列出了累计剂量,T omotherapy出束总剂量比IMRT和RapidArc平均高20%。受到>2Gy和>5Gy剂量的体积也在表3中列出,总的趋势是Tomotherapy受到>2Gy和>5Gy的体积最大,其次RaptidArc,IMRT计划体积最小。
* ]! W( f8 H5 Z% D' T$ X$ w( s
+ I" y T2 j8 YE.治疗计划时间和治疗出束时间估计! t( f; Z8 b5 P% i8 ?) z5 _
$ o J0 ?4 ^" [8 m
表4是计划时间和治疗出束时间估计值的一份总结。IMRT平均计划时间7.5分钟,RapidArc48分钟,Tomotherapy59分钟。IMRT出束时间估计4.8分钟,RapidArc2.2分钟,Tomotherapy3.5分钟。基于头颈215子野的IMRT,每个IMRT野的平均delta值是19.12.2s/野。
3 r: j9 ^# c4 E* k# W9 p2 k3 |9 A. l9 i. g- j. l9 A' A7 G9 F
表格4 模体1-4:IMRT、RapidArc、Tomotherapy的计划设计治疗计划时间和治疗出束时间估计值表。
! d# \, g% J% X5 w
. N# g4 d$ O7 s + m/ G. q K1 |$ y
$ c+ c7 _8 R! q! c7 } 1模体 2模体 3模体 4模体 ' A$ U4 Z, P1 d! r& \8 C- Q. Z
9 L0 x# V1 p5 _4 n/ W
计划 出束 计划 出束 计划 出束 计划 出束
5 ~- J8 e* d; n
" ?9 Y( Z7 M0 wIMRT 5野 6.8 3.8 6.7 3.9 6.7 4.0 7.3 4.3
; f0 t7 D0 X! M1 J" h+ B
D9 v6 K+ }& Z3 a2 I; F' h 9野 8.1 5.2 7.5 5.3 7.8 5.8 9.1 5.9
' z% m; L) Y" f) ?
: ?( e% V2 S. f" `8 Q4 hRapidArc 1弧 45.0 1.4 38.0 1.4 30.0 1.7 53.0 1.7* k+ I8 V; E1 B1 B9 n- Z+ h
: H$ v5 b- ]4 E) D2 d0 V( I 2弧 60.0 2.8 47.0 2.8 50.0 2.8 62.0 2.8 ( v/ ~. }2 Z# Y2 N' U3 \
) j0 l% ~9 I- f, `
TOMO 45.0 4.3 44.0 2.8 44.0 2.8 101.0 3.9
/ j' O# }+ q- Q: i8 {" `4 A- V% p& F
* T0 M# p' e% z& n8 k5 s, L0 q s6 d/ P# z, i. H3 {! w& v4 A
IV.讨论, N- G- g9 \' e4 ^6 {* M" y
0 i i2 K {8 ~" H
对于所有模体,Tomotherapy能满足最大限度满足优化限制(2/4P1,3/6P2,2/4P4),其次为RapidArc(1/4P1,1/6P2和0/4P3),最后是IMRT(0/4P1,0/6P2和1/4P4)。这些结果是依治疗技术作出的总结,几乎与所采用射野/弧个数独无关的。Tomotherapy计划对所有模体说,均能做出PTV最佳的剂量均匀性(DHITOMO>DHIRA2>DHIRA1>DHIIMRT9>DHIIMRT5)。这个满足DVH限制能力的提高是一其他的成本作为代价的,包括:较长的计划时间(tplan,tomo>tplan,ra2>tplan,ra1>tplan,imrt9>tplan,imrt5),较长治疗时间(tplan,imrt9>tplan,imrt5 >tplan,tomo>tplan,ra2>tplan,ra1),较低适形度指标(RCITOMO<RCIIMRT5<RCIIMRT9<RCIRA1<RCIRA2),和较高累积剂量(大约是IMRT或者RapidArc的1.2倍)。1 q" k2 i; j% ?5 y
1 P* c9 h! Y# u- g! D 对模体3,没有一个治疗计划能满足剂量体积限制。然而,Tomotherapy还是能为计划提供最好的剂量均匀性(DHITOMO>DHIRA2>DHIRA1> DHIIMRT5> DHIIMRT9),但适形度较低(RCIRA2> RCIRA1> RCITOMO > RCIIMRT9> RCIIMRT5)。关于模体3OAR的平均剂量,没有差异趋势。4 Q* U# N8 [' r P4 ~
7 f3 A4 b5 l. _ 考虑在治疗计划时间、治疗出束时间和计划总体质量之间的权衡,本研究说明不同的治疗技术间各不相同。很明显,Sliding Window IMRT技术相比RapidArc或Tomotherapy,可以在短得多的时间内生成治疗计划;相比于IMRT和Tomotherapy,RapidArc(不论是单弧或双弧)有最低的治疗出束时间。对于计划质量,虽然Tomotherapy可以满足绝大多数剂量-体积限制,而且一般可以使靶区体积的均匀性做到最好。即使TOMO在优化过程中形成剂量分布上灵活性更高,但Tomotherapy的平均适形度指标却差于IMRT和RapidArc,我们的结果至少有两个理由可以解释。首先,两个治疗计划系统(TOMO和Eclipse)可能在判断体素在勾画轮廓内外上有不同标准。其次,为了满足Tomotherapy的靶区剂量,剂量提升区至少需要和扇形野厚度一致,这会导致靶区体积的上下方向额外剂量照射。图6显示的9野IMRT、双弧RapidArc和TOMO穿过模体4中心的上下方向剂量分布,就说明了这种剂量提升效应。
1 m& l* K, \8 F9 `7 d* S" f* F# O
本研究采用了简单通过的方法,在多次循环中反复优化参数以获得一个与DVH限制符合的最好计划的方法没有被采用。需要一个更有说服力的方法去评价IMRT、RA和TOMO的计划质量等级。例如,通过反复调节优化参数,相比一个简单通过方法,专业计划者可以获得更好的计划。另外,为每种治疗技术都招募专业计划人员将会对研究有益,但我们依然可以看到同样的治疗技术排序,因为对固定的计划设计时间,IMRT可以在短时间内完成一个计划,而且相对于RA和Tomo,IMRT有更多的优化参数可以调整。0 I0 n& l1 f4 l& u, w2 H
7 i! U2 }) w$ z; Q" {% [ 治疗计划研究的另外一项缺点就是来自不同治疗计划系统数据之间的相互比较,由于存在较多复杂因素可能不公平。其中一个因素即Deasy所谓的“Weight Pardox”(权重悖论),即对不同PTV和OAR优化限制相对权重有优选并不明晰,可能需要大量重复地试验和错误判别。DVH限制的选择也是在OAR和PTV覆盖相互妥协的一个限制,也可能在用户不期望的计划结果上终止优化。这些剂量障碍已经通过两种不同的方法包括多限制优化和Pareto优化加以考量,如在商业治疗计划系统运行就可能提供较好的结果。另一个困惑就是不同制造商的治疗计划系统可能定义定PTV和OAR的功能和权重其意义各不相同。这个信息通常用户都不清楚,也很难去干涉。
( E' }; ^3 f" s. v" K- U
5 |8 Q( D: V+ P$ d% U. A; U. C9 p% a 本研究的一个缺点是模体在2D轴断面轮廓扩展而成,不能代表在临床环境中所能看到的轮廓。这就不能确定Tomotherapy系统在上下方向轮廓改变的PTV和OAR模体计划设计时是否会存在困难。正如Bortfeld和Webb文章中所提及的,在上下方向上计划时,Tomotherapy需要在扇形射野宽度的选择、调制因子和螺距等因素间做有效的权衡。在上下方向上形状变化的结构,这些权衡因素怎样去影响计划质量和治疗效率,需要进一步去进行研究。
, W8 p- ~- u+ U8 ? ?) B; p7 Z) ?9 }& u5 d) z6 F
最后,很难去确切地知道优化运算收敛到了局部最小量或者现有方案与整体最小化间距离有多大。虽然运算循环的次数现在可以直观显示给用户,但现在还没有一个限制可以在优化已经集达到设定好标准(比如说100循环中目标函数有1%差异)时停止优化。如果有这种停止优化的限制,那么相应的优化算法的收敛性应做出说明。
+ x5 B6 A2 D! T% Y8 l* ^; W6 ?% a: m6 w& e- {0 N
进一步理解在临床几何条件下的不同治疗技术的相对优势还需要做进一步的研究。) J; F* L% ]* B/ ~
# w# b4 |3 q# {0 u% ?
V.结论
: _4 I; L; z) o2 j
! u) h: B1 [+ a% v& g7 e! h 每项先进的放射治疗技术都有它们各自的相对优势。对本研究使用的四个模体,5和9野IMRT技术可以在较短时间内做出计划,以最低累积剂量出束。单弧和双弧RapidArc计划能在最短时间内出束,而且能提供PTV最好的适形度。最后,Tomotherapy能最大地满足计划要求,而且能提供PTV最大的剂量均匀性。5 n; G" g/ F2 d% N) b/ _& ~
4 f( o: u( r r4 n
完成
' _; [+ h9 _1 W* b1 R
) d- c6 f& B7 v, {5 U# T) V: |; ^/ E- s# k 我们由衷感谢英国癌症计划组织允许我们在这项研究中使用TOMO计划工作站。 |