食管癌患者CBCT图像剂量计算准确性的分析评估

巩汉顺，徐寿平，刘博，蔡博宁，丛小虎，陈高翔，崔泽平

1.解放军总医院放疗科，北京100853；2.北京航空航天大学宇航学院，北京100083

前言

食管癌是我国高发的消化道恶性肿瘤，每年新发病例占全世界约50%［1］，放射治疗是目前食管癌主要、安全且有效的手段之一［2］，特别是对于局部复发患者采用放疗具有明显的生存获益［3］。近年来随着适形及容积旋转调强等新技术的广泛应用，使得食管癌患者在危及器官得以充分保护的前提下保证肿瘤靶体积达到设定处方剂量［4-5］。通过与分次锥形束CT（CBCT）位置配准降低摆位误差，提高治疗精度［6］，同时因CBCT获得了患者在分次放疗过程中解剖结构实时影像，为实施治疗计划的再优化奠定基础，进而实现食管癌患者的自适应放疗可能；但由于CBCT在图像质量及HU值方面与扇形束CT（kVCT）存在一定差异，直接用于剂量计算势必会产生一定偏差，本研究对食管癌患者CBCT图像剂量计算的准确性进行分析研究。

1 材料与方法

1.1 病例筛选及kVCT计划设计

随机选取在解放军总医院放疗科经Varian 加速器（Clinical iX）治疗的食管癌患者8 例。kVCT 图像采用西门子CT（SOMATOM Definition AS）扫描获得，患者均采用仰卧位，双臂上举至头部，采用热塑体模进行体位固定，所得图像采用DICOM格式传输至Varian Eclipse（Ver 10.0）计划系统，由经治医师勾画靶区及危及器官。上级医师给定靶区处方剂量为原发病灶pGTV：62 Gy/30 F，亚临床病灶PTV：54 Gy/30 F。责任物理师设计kVCT 计划：采用RapidArc技术，设计两个全弧射野（CW：180.1°～179.9°；CCW：179.9°～180.1°）。

1.2 患者CBCT图像扫描

采用OBI系统对患者采用预设的LDT（Low-Dose Thorax）临床扫描条件进行扫描，图像采用3 mm层厚重建获得。为减小由于患者解剖结构变化而造成的剂量计算偏差，各CBCT计划设计均采用第一次治疗时所获得的CBCT图像。

1.3 直方图匹配算法

为减小由于图像HU值差异而产生的剂量偏差，进一步提高CBCT图像剂量计算准确性，实现CBCT与kVCT图像HU值的归一具有一定的必要性。直方图匹配算法［7］是实现HU 值归一的常用方法之一，基于kVCT图像HU值的积分直方图调整相应CBCT图像的HU值以生成新的CBCT重建图像［8］。

1.4 CBCT图像计划创建

OBI 系统预设的LDT 模式一般采用Half-fan 扫描条件，图像重建范围为最大直径45 cm，轴向长度15 cm，对于CBCT图像由于FOV太小而造成的组织缺损部分，利用自主研发的ART软件将kVCT图像与CBCT图像形变配准完成组织填充以形成新的CBCT图像序列，保证CBCT计划在剂量计算时具有足够的组织信息散射，进而确保CBCT图像剂量计算的准确性。

在ART 软件中利用CBCT 图像到kVCT 图像的形变场，以kVCT图像CT值为基准，采用直方图匹配算法对形变后的CBCT 图像的CT 值实施修正；同时将kVCT计划的射野参数、靶区和危及器官信息映射到形变配准后的CBCT 图像序列上以生成相应的CBCT计划及组织结构。

1.5 CBCT剂量计算

将新生成的CBCT计划的图像序列、射野参数及组织结构重新导入到Eclipse计划系统中，调用kVCT图像的CT 值-电子密度校准（HU-ED）曲线，采用与kVCT 计划相同的机器跳数（MU）及剂量体积算法（Dose Volume Optimizer）重新完成剂量计算。

1.6 数据分析评估

所得CBCT 图像中各组织剂量与相应kVCT 进行分析比较，靶区参数：Dmean、D2%、D98%、V95%及V100%，其中D2%为2%靶区体积接受的最高剂量；D98%为98%的靶区体积接受的最低剂量；V95%、V100%分别为处方剂量95%及100%包绕的靶区体积百分比。危及器官参数：Dmean、D2%、V5、V20及V30等。

1.7 γ分析

将kVCT计划的CT 图像序列、组织结构、计划剂量以及CBCT计划的剂量分别导入到CERR（Computational Environment for Radiotherapy Research,Ver5.1）软件中，以kVCT计划的剂量场为基准对CBCT计划的剂量场进行分析比较，采用一定的γ 分析标准（2 mm/2%/TH:10%；1 mm/1%/TH:10%）定量评估CBCT计划剂量计算精度。

2 结果

2.1 数据分析结果

对于CBCT图像中靶区和危及器官的剂量、体积参数相对于kVCT 偏差绝对值均小于1%，具体数据如表1所示。对于靶区Dmean及D2%结果CBCT均高于kVCT，而V100%值却与之相反。各危及器官剂量参数Dmean和D2%：CBCT 图像中肺组接受剂量略高于kVCT，而心脏及脊髓却低于kVCT。全肺及心脏的低剂量体积参数V5：CBCT 相对于kVCT 均有所增大，而高剂量体积参数V20、V30则有一定程度的缩小。

2.2 γ分析结果

按照2 mm/2%/TH:10%分析标准，各组织结构CBCT计划剂量分布相对于kVCT的平均γ通过率均达99%以上；按照1 mm/1%/TH:10%分析标准，靶区平均γ 通过率下降明显，最低为（82.59±16.16）%，各危及器官平均γ 通过率略有下降，仍可达98%以上。全局平均γ 通过率无论采用上述何种分析标准均可达到99%以上，结果如表2所示。

3 讨论

食管解剖结构及位置比较特殊，上起于相当于第6颈椎下缘的环状软骨，下止于相当于第11胸椎水平的胃贲门，其长度成人一般为25～30 cm，在放疗过程中解剖位置受多种因素影响（心血管搏动、呼吸运动等［9］），存在一定的不确定性；不同部位的病灶运动幅度也有所变化［10］，胸下段食管癌内靶区头脚方向边界外放距离大于胸上、中段［11］；各患者分次内靶区运动位置也存在较大差异［12］。分次放疗中摆位误差、呼吸及器官运动和肿瘤组织退缩等因素均会引起所受剂量出现偏差［13］，需对重新生成的CBCT计划及剂量进行评估，进而选择合适的“每日计划”［14］。

表1 8例食管癌患者CBCT计划剂量计算数据相对于kVCT的分析结果（±s）Tab.1 Deviation of cone beam CT dose calculation relative to kVCT in 8 patients with esophageal cancer(Mean±SD)

表1 8例食管癌患者CBCT计划剂量计算数据相对于kVCT的分析结果（±s）Tab.1 Deviation of cone beam CT dose calculation relative to kVCT in 8 patients with esophageal cancer(Mean±SD)

结构左肺右肺全肺心脏脊髓pGTV PTV Dmean/cGy 0.26±0.81 0.30±0.75 0.27±0.77-0.22±0.54-0.22±0.48 0.09±0.57 0.04±0.43 D2%/cGy 0.07±0.52 0.05±0.83 0.07±0.67-0.11±0.52-0.24±0.55 0.27±0.65 0.23±0.60 D98%/cGy- - - - --0.06±0.44 0.12±0.38 V5/%-0.08±0.75 0.30±0.81 0.14±0.75 0.03±0.19- - -V20/%0.15±1.57-0.24±1.69-0.07±1.61-0.24±0.45- - -V30/%-0.53±2.31-0.73±1.63-0.07±1.74-0.59±1.18- - -V95%/%- - - - --0.02±0.03 0.04±0.10 V100%/%- - - - --0.05±0.64-0.07±0.63

表2 8例食管癌患者CBCT计划剂量分布相对于kVCT的γ通过率分析结果（%, ±s）Tab.2 Comparison of gamma passing rate of CBCT plan dose distribution relative to kVCT in 8 patients with esophageal cancer(%,Mean±SD)

表2 8例食管癌患者CBCT计划剂量分布相对于kVCT的γ通过率分析结果（%, ±s）Tab.2 Comparison of gamma passing rate of CBCT plan dose distribution relative to kVCT in 8 patients with esophageal cancer(%,Mean±SD)

结构左肺右肺全肺心脏脊髓pGTV PTV皮肤2 mm/2%/TH:10%99.96±0.06 99.91±0.10 99.93±0.07 99.98±0.05 100.00±0.00 99.46±0.90 99.23±0.97 99.98±0.02 1 mm/1%/TH:10%99.25±0.36 98.48±0.86 98.83±0.48 98.99±1.44 99.94±0.18 82.59±16.16 85.35±9.25 99.47±0.32

CBCT 在图像质量及HU 值与kVCT 有所区别，利用此图像开展自适应放疗时产生一定的剂量偏差。同时由于食管位于两肺之间，而肺为空腔组织与周围实体组织相比存在较大电子密度差异，导致该部位CBCT 图像HU 值与kVCT 图像差异较大，为保证剂量计算的准确，对CBCT图像的HU值进行适当校正就显得尤其重要［15］。Dunlop等［16］对于获取的肺部CBCT 图像使用RSauto中关于CBCT 密度分配的方法及设定为水或肺密度覆盖方法进行CT值校正，所得剂量分布与kVCT较为相似，平均绝对剂量差异分别为1.3%和1.4%。Poludniowski 等［17］通过对CBCT 投影数据采用散射校正以提高图像重建质量并以此为基础进行剂量计算，与kVCT剂量偏差小于2.5%。本研究中利用直方图匹配方法对食管癌患者CBCT图像的HU值进行修正所得CBCT的组织结构剂量与kVCT 偏差绝对值均小于1%，其精度可较好满足临床的应用要求。

利用治疗过程中的CBCT图像反馈信息，对放疗过程作出相应调整，以制定个体性自适应放疗方案［18］。Hawkins等［19］对食管癌患者采用首周的CBCT图像来修正肿瘤靶区及设计自适应计划，结果显示靶区回缩明显，肺和心脏的平均剂量也降低明显。Sriram 等［20］对于食管癌患者利用第一周CBCT 图像数据设计自适应容积调强计划的策略，可以有效降低危及器官剂量及正常组织并发症发生概率。剂量学因素是食管癌患者发生放射性肺炎的独立预后因子［21］，其中肺Dmean是放疗过程中肺损伤的关键因素［22］。

4 结论

本研究对食管癌放疗患者采用CBCT 图像进行剂量计算的精度进行分析，验证其剂量计算的可行性，为自适应放疗中采用CBCT图像直接进行剂量计算奠定了一定的基础。