基于图像配准分析的锥形束CT在肺癌患者精确放疗中的摆位误差研究
张西志1,李军1,张先稳1,桂龙刚1,孙新臣2
(1.苏北人民医院放疗科,江苏扬州225001;2.南京医科大学第一附属医院放疗科,江苏南京210029)[摘要]目的利用VARIAN Clinac-IX加速器锥形束CT( CBCT)功能,研究刚性图像配准算法种形变图像配准算法两者哪种更适合于与计划CT图像进行配准,并对10例肺癌患者进行摆位精确度验证,分析其摆位误差并探讨靶区的外放边界。方法随机选取图像引导放疗( IGRT)的10例肺癌患者,分别采用刚性图像配准算法和形变图像配准算法2种算法配准CBCT图像和计划CT图像,并根据轮廓指标进行评估,选用优的图像配准方案,并统计每例患者每天在左右(x)、头脚(Y)、垂直(z)方向上的摆位误差,分析摆位误差结果并根据公式计算出肺癌靶区的临床外放边界。结果通过刚性图像配准与形变图像配准比较,形变图像配准能更好获得不同模态或相同模态不同时期的图像中器官形态、位置的变化。患者在x、Y、z轴的摆位误差分别为(0.04 +0.16)cm、(0.09 +0.32)cm、(-0. 02 +0.13)cm,患者在盖、y、Z形成的u、”、Ⅲ旋转误差分别为(0. 41 +0. 64)。、(-0. 20 +0. 57)。、(0.06±0.62)。。计算得出肺癌靶区的临床外放边界在X、l,、Z等3个方向上分别是0.78 cm、1.12 cm、0.68 cm。结论胸部肿瘤放疗的摆位误差在头脚方向为明显,利用CBCT并采用形变图像配准算法可以测量肺癌患者的摆位误差并在线校正,减少摆位误差,根据误差值,估算出合理的靶区外放值,从而更好治疗靶区和保护危及器官。
[关键词]肺癌;锥形束CT;摆位误差;外扩边界
DOI:10. 3969/j. issn. 1673 -5412. 2018. 04. 017
[中图分类号] R734.2;R730. 55 [文献标识码]B [文章编号]1673 - 5412( 2018) 04 - 0336 - 04
放疗已经进入精确治疗时代,患者在肿瘤放疗中都会产生误差,这种误差主要来源于器官运动导致的靶区位移误差和摆位误差[1]。肺癌患者由于受呼吸运动的影响较大,其治疗误差往往较大。图像引导放疗( image guided radiotherapy,IGRT)是实现精确放疗的重要手段之一,IGRT通过在治疗前或治疗中获取图像来引导治疗,以减少患者放疗期产生的摆位误差心。41。利用千伏级锥形束CT( conicalbeam CT.CBCT)进行图像在线校正是目前运用广泛的IGRT技术。本文将比较CBCT技术中刚性图像配准算法和形变图像配准算法的价值,并以肺癌放疗为例,统计其治疗过程中的平移误差和旋转误差,利用经验公式来估算临床靶体积( clinical target volume,CTV)外放计划靶体积(planning target volume,PTV),PTV在左右(X)、头脚(Y)、垂直(Z)方向上等3个方向上的外放值。
1 材料与方法
1.1 病例选择 2015年3月至2015年8月采用VARIAN Clinac-IX加速器进行IGRT的10例肺癌患者,其中男7例,女3例;年龄范围46~ 73岁,中位年龄55岁。病理类型包括小细胞肺癌2例、肺腺癌5例以及肺鳞癌3例。
1.2体位固定和CT定位所有患者采用仰卧位,用真空负压垫进行体位固定,并在体表和真空负压垫上作好定位标记,并用16排70 cm大孔径CT进行断层扫描,扫描范围为环状甲膜至纵膈肌下缘,扫描层厚5mm,扫描后的CT图像通过Varian Aria网络系统传输至Eclipses TPS系统中。
1.3计划CT与CBCT图像配准 每例患者分次治疗前均采用VARIAN Clinac-IX加速器CBCT扫描,扫描角度为顺时针182。至180。,总的扫描角度为358。。1)将TPS计划中的CT图像和获取的CBCT图像利用自动刚性配准和刚性辅助校准方法来完成刚性图像配准;2)利用基于像素的形变图像配准算法将TPS计划中的CT图像和CBCT图像完成配准,并进行轮廓自动勾画,并将其与CBCT勾画的轮廓进行对比分析,分析指标包括相似性系数( DSC)、阳性灵敏度(SENS)和预测值( PPV),其中DSC可反映2种轮廓的重叠区域及平均值,SENS反映了阳性检测的样本占总数的百分比,PPV是指阳性检测的样本总数中,样本占阳性检测的百分比。
1.4 CTV外扩分析 由于人体器官非刚性,且摆位误差不叮避免,所以为了减少治疗偏差,通常应对CTV的外边界进行外扩。具体来讲,由人体器官不自主运动所引起的C‘rv外边界的运动范围称为内扩边( internal margin,IM),由摆位误差引起CTV边界运动范围称为摆位扩边(setup margin,SM)[5]。在图像配准过程中,可采集得到每例患者在X、y、Z方向上的误差值。根据这些误差值,我们通常可利用经验公式M=2.5∑+0.78来计算CTV在X、y、Z等3个方向上的外放边界值[6-7]。其中8为所有患者误差值的标准差,该值可认为是随机误差;∑为所有患者误差值的标准差的平均值,该值可认为是系统误差;M为CTV外放成PTV的外扩估计值。
1.5统计学处理采用SPSS 18.0进行数据分析,计量资料用x+s表示,比较用F检验和£检,检验水准d=0. 05。
2 结果
2.1 刚性图像配准算法与形变图像配准算法DSC值比较刚性图像配准算法和形变图像配准算法的DSC值比较结果显示,对于体表轮廓、左肺、右肺及心脏形变图像配准算法的DSC值均高于刚性图像配准算法(P均<0. 05)。见表1。
表1 刚性图像配准算法与形变图像配准算法DSC值比较
器官
| 刚性图像
配准算法 | 形变图像
配准算法 |
£
| P
|
休表轮廓 | 0. 962 +0. 052 | 0. 971 +0. 012 | 3. 288 | <0. 05 |
左肺 | 0. 846 i0. 123 | 0. 883 +0. 151 | 4. 205 | <0. 05 |
右肺 | 0. 883 +0. 062 | 0. 898 +0. 039 | 3. 179 | <0. 05 |
心脏 | 0. 748±0 141 | 0. 774 +0. 153 | 4. 346 | <0. 05 |
脊髓 | 0. 636 +0. 114 | 0. 657 +0. 116 | 2. 521 | <() 05 |
2.2摆位误差 10例肺癌患者共进行254次CBCT扫描和配准,根据配准结果,分别采集X、y、Z等3个方向E的平移误差和旋转误差。患者在X、Y、Z轴的摆位误差分别为(0.04 +0. 16) cm、(0.09 +0.32)cm、(-0.02±0.13)cm,患者在X、l,、Z轴形成的M、u、训旋转误差分别为(0. 41 +0. 64)0、(-0.20±0.57)。、(0. 06 +0. 62)。。见图1、2,表2、3。
图1 254次CBCT扫描误差横断面散点图
图2患者在u、u、训方向上的旋转误差散点分布图
表2 10例肺癌患者在X、l,、Z等3个方向上的摆位误差 cm
编号 | x | y | z |
1 | -0. 09±0.34 | -0. 29±0.40 | 0. 19 i0. 28 |
2 | -0. 03土0.32 | 0. 12±0.31 | -0. 13±0.18 |
3 | 0. 19±0.22 | - 0. 28土0 47 | -0. 17±0.23 |
4 | 0. 14±0.30 | 0. 34±0.51 | -0. 11±0.22 |
5 | 0. 24±0.28 | 0.33±0.25 | 0. 10±0.17 |
6 | -0. 08±0.20 | 0. 22±0.36 | -0. 03±0 17 |
7 | -0. 15±0.23 | 0. 65±0.31 | 0. 11±0.24 |
8 | 0. 22±0.18 | 0. 26±0.30 | 0. 05±0.20 |
9 | -0.18±0. 36 | -0. 35±0.23 | -0. 21±0. 34 |
10 | 0. 12土0.28 | -0. 13±0.38 | -0. 04±0.38 |
表3 10例肺癌患者在u、训等3个方向上的旋转误差
编号 | M | u | |
1 | 0. 69±0. 71 | -0. 21±0.63 | - 0. 74土0.64 |
2 | 0..83±0.56 | 0. 32±0.56 | 1. 02±0.72 |
3 | 1. 02±0.83 | 0. 96±1.03 | 0. 26±0.75 |
4 | 1. 43±1.20 | -0. 08±0.88 | -0. 56±1.37 |
5 | -0. 12±0.50 | -0. 48±0.48 | 0. 47±0.74 |
6 | 0. 46±0.48 | -0. 42±0.83 | 0. 42±1.03 |
7 | 0. 89±0.48 | -0. 60±0.65 | - 0.18±0.70 |
8 | -0. 23±0.53 | -0. 50±0.35 | 0. 73±0.67 |
9 | - 0.33±1.23 | -0 58±0.53 | -0. 10±0. 88 |
10 | -0.59 t0.65 | -0. 44±1.01 | -0. 76土0.92 |
2.3靶区外扩情况肺癌靶区的临床外放边界在X、Z等3个方向上分别是0.78 cm、1.12 cm、0.68 cm。见表4。
表4 10例患者在X、l,、Z等3个方向上的外扩估计值 cm
内容 | x | l, | z |
平均平移误差 | 0. 06 | 0.10 | -0. 03 |
系统误差(∑) | 0. 27 | 0. 36 | 0. 24 |
随机误差(8) | 0. 15 | 0. 31 | 0. 12 |
外扩边界估计值(M) | 0. 78 | 1.12 | 0. 68 |
3讨论
随着设备和技术的发展,放疗已经进入精确时代,主要体现在对影像获取、靶区定位及靶区勾画等方面的精确要求愈来愈高。但患者在实际治疗过程中的治疗精度仍受摆位误差和器官运动等因素的显著影响。例如,多项研究[10]表明肺癌放疗过程中产生的分次间摆位误差约为5~ 40 mm。而Balter等[11]的研究结果表明当摆位误差大于1 cm时对靶区的影响约为6 mm。因此,为了实现肿瘤的精确治疗,必须进行摆位校正以减少摆位误差。摆位校正目前主要依靠CBCT图像获取和与计划CT图像的配准技术。Borst等[12]认为若没有CBCT图像配准来减少误差,510_10的肺癌患者可能存在超过5 mm的摆位误差。Purdie等[13]对28例T,~:期非小细胞肺癌立体定向放疗的研究显示,CBCT在肺癌立体定向放疗中的在线图像指导能够确保靶区位置精确性。
本文研究结果显示,与采用刚性图像配准算法相比,采用形变图像配算法能获得较好的DSC值(包括双肺、脊髓、心脏、体表轮廓),通过对10例肺癌患者进行的254次CBCT扫描和图像配准,得到患者在x、y、Z轴上的平均平移误差分别为(0. 06±0.15)cm、(0.10 +0.31)cm、(-0.03 +0. 12)cm。该结果表明Y轴方向平均平移误差大,X轴方向次之,Z轴方向小。系统误差在X、y、Z方向分别为0. 27 cm、0. 36 cm、0.24 cm,其中以Y轴方向误差大。X、,,、Z轴方向的随机误差分别为0. 15 cm、0.31 cm、0.12 cm,同样以Y轴方向误差大,这是因为与利用热塑体膜固定体位有关,在X轴和Z轴方向上,体膜相对固定,不容易发生相对位移,而l轴方向两端并不是封闭的,假若抱头姿势稍有改动或皮肤松弛都可能引起,轴方向上的显著位移偏差。患者治疗的分次间误差由摆位误差与器官移动引起的误差组成,因此在靶区勾画时需要对CTV进行外扩形成实际治疗的PTV。本研究采用经验公式来估算CTV的外扩值,该方法同时考虑了系统误差和随机误差的影响,其结果能保证90%群体患者小CTV累积剂量至少达到95%的处方剂量,根据结果推荐推荐CTV外扩PTV应在X、l,、Z方向分别外扩0.78 cm、1.12 cm、0.68 cm。
本研究也存在一定不足,例如,本研究中的所有患者均是利用真空负压垫来固定体位的,而使用其他定位技术(例如碳素纤维定位架、热塑体膜等)可能产生的误差也不同。且在本研究中,并未对左、右肺癌患者进行独立分析比较,由于膈肌和心脏运动对左、右肺影响不同,很可能引起左、右肺靶区的位移,所以在后续研究中,有必要对左、右肺癌的结果的差异进行进一步研究。
参考文献:
[1]戴建荣,胡逸民,图像引导放疗的实现方式[J].中华放射肿瘤学杂志,2006,15(2):132 - 135.
[2]MCBAIN C A, HENRY AM, SYKES J,et al.X-ray volumetric maglng in image-guided radiotherapy: the new standard in ontreatment imaging[J].Int J Radiat Oncol Biol Phys ,2006 ,64(2):625 -634.
[3]曹晓辉.CBCT图像引导下胸部肿瘤放射治疗摆位误差分析及对剂量分布的影响[D].石家庄:河北医科大学,2011.
[4]巩汉顺,鞠忠建,韩嫒嫒,等,图像引导放射治疗技术提高腹部肿瘤放射治疗精度的临床研究[J].中华临床医师杂志:电子版,2012, 6(14):3887 -3890.
[5]巩汉顺,王运来,鞠忠建,等.锥形束CT在肺癌患者精确放疗摆位误差研究[J].解放军医学院学报,2012( 10):1037 -1038.
[6]VAN HERK M. Eeeors and Margins in Radiotherapy [J]. ScminRadiat Oncol ,2004,14:52 -64.
[7]VAN HERK M, REMEIJER P,I.EBESQUE JV. Inclusion of geometric uncertainties in treatment plan evaluation[J]. Int J Radiat Oncol Biol Phys, 2002 ,52:1407 - 1422.
[8]王勇,应用KV-CBCT研究及管癌放疗的摆位误差及其对剂量分布的影响[D].河北医科大学,2012.
[9]SHIMIZU SS, SHIRATO H, KAGEl K, et al.Impact of respiratory movement in the computed tomographic images of small lung tumors in 3D radiotherapy[Jl. Int J Radiat Oncol Biol Phys,2000, 46(5):1127 - 1133.
[10]EKBERC L,HOLMBERG 0, WITTCREN I.,et al.What margins should be added to the clinical target volume in RT treatment planning for lung cancer[J]. Radiother Oncol, 1998, 48(1):71 - 77.
[11]BALTER JM, CHEN CT, PELIZZARI CA, et al. Online repositioning during treatment of the prostate:a study of potential limits and gains[J].Int J Radiat Oncol BioPhys ,1993 ,27:137 - 143.
[12]BORST GR,SONKE JJ,BETQEN A,et al.Kilo-voltage cone-beam computed tomography setup measurements for Jung cancer patients;first c.linical results and comparison with electronic portal-imaging device[J].Int J Radiat Oncol Biol Phys ,2007 ,68:555 - 561.
[13]PURDIE TG,BESSONNETTE JP,FRANKS K, et al.Cone-beam computed tomography for on-line image guidance of lung stereotacic radiotherapy: localization,verification,and intrafraction tumor position[J].Int J Radiat Oncol Biol Phys ,2007 ,68:243 - 252.