调强适形放射治疗剂量验证中CT密度值的校准分析

庄名赞 吴仁华 邱庆春 彭逊 陆佳扬 陈志坚

引用本文:
Citation:

调强适形放射治疗剂量验证中CT密度值的校准分析

    通讯作者: 吴仁华, rhwu@stu.edu.cn ; 邱庆春, qcqiu@stu.edu.cn

Calibration of CT density values in dosimetry verification of intensity modulated radiation therapy

    Corresponding author: Ren-hua WU, rhwu@stu.edu.cn ;Qing-chun QIU, qcqiu@stu.edu.cn ;
  • 摘要: 目的 基于调强适形放射治疗(IMRT)验证体模,分析放射治疗计划系统中CT密度转换曲线在IMRT剂量验证中的影响,并加以校准。 方法 利用CT电子密度体模刻度计划系统中的CT密度转换曲线,取12例鼻咽癌患者的IMRT计划,分别移植至IMRT验证体模,计算剂量分布,并于剂量均匀处放置测量点,利用电离室测量出其实际剂量。将IMRT验证体模的物理密度与对应的CT值输入计划系统,对CT密度转换曲线进行修正,重新计算剂量,在保持其他计划参数一致的前提下,比较前后两次计算值与实测值的差异。 结果 12例鼻咽癌IMRT计划中,前后两次计算值和实测值的平均误差分别为1.96%±0.87%和0.63%±0.74%,修正后的计算值误差控制在±2%以内,而修正前计算值的最大误差为3.24%。 结论 利用修正后的CT密度转换曲线,计划中的剂量计算值更加贴近实际测量值。为提高IMRT剂量验证的精确性,验证体模在使用之前必须进行CT密度值的校准、确认。
  • 图 1  美国CIRS公司Model 062电子密度体模及插件示意图

    图 2  计划系统中CT密度转换曲线(Pinacle中CT值=CT图像的CT值+1000) CT图像的CT值+1000)

    表 1  美国CIRS公司Model 062电子密度体模的密度参数表

    物理密度(g/cm3) 电子密度(x1023/m3) 相对水电子密度
    体模(水等效) 1.01 3.346 1.002
    肺(吸气) 0.20 0.634 0.190
    肺(呼气) 0.50 1.632 0.489
    脂肪 0.96 3.170 0.949
    乳腺 0.99 3.261 0.976
    肌肉 1.06 3.483 1.043
    1.07 3.516 1.052
    骨(200mg/cm3) 1.16 3.730 1.117
    致密骨(800 mg/cm3) 1.53 4.862 1.456
    下载: 导出CSV

    表 2  12例鼻咽癌IMRT计划前后两次计算值与实测值的比较分析

    实测值(cGy) 计算值[cGy (误差,%)]
    修正前 修正后
    例1 210.2 217(3.24) 214(1.81)
    例2 207.4 211(1.74) 208(0.29)
    例3 203.4 205(0.79) 202(-0.69)
    例4 194.2 196(0.93) 194(-0.10)
    例5 202.4 206(1.78) 204(0.79)
    例6 203.1 205(0.94) 203 (-0.05)
    例7 200.2 205(2.40) 202(0.90)
    例8 199.9 204(2.05) 201(0.55)
    例9 196.8 203(3.15) 200(1.63)
    例10 193.6 196(1.24) 194(0.21)
    例11 192.1 197(2.55) 194(0.99)
    例12 202.5 208(2.72) 205(1.23)
    下载: 导出CSV
  • [1] Cozzi L, Foqliata A, Buffa F, et al. Dosimetric impact of computed tomography calibration on a commercial treatment planning system for external radiation therapy. Radiother Oncol, 1998, 48(3): 335-338. doi: 10.1016/S0167-8140(98)00072-3
    [2] 祁振宇, 黄劭敏, 邓小武. 放疗计划CT值的校准检测及其影响因素分析. 癌症, 2006, 25(1): 110-114.
    [3] Tailor RC, Hanson WF, Ibbott GS. TG-51: experience from 150 institutions, common errors, and helpful hints. J Appl Clin Med Phys, 2003, 4 (2): 102-111. doi: 10.1120/jacmp.v4i2.2524
    [4] 胡逸民. 肿瘤放射物理学. 北京: 原子能出版社, 1999.
    [5] 杨克柽. 螺旋CT图像性能参数检测方法研究. 广州: 第一军医大学, 2002.
    [6] Mutic S, Palta JR, Butker EK, et al. Quality assurance for computed-tomography simulators and the computed-tomography simulation process: report of the AAPM Radiation Therapy Committee Task Group No. 66. Med Phys, 2003, 30(10): 2762-2792. doi: 10.1118/1.1609271
    [7] Beneventi S, Chionne F, Gobbi G, et al. Quantitative CT tomography for radiotherapy treatment planning: Calibration phantom and sources of error. Radiother Oncol, 1995, 37(suppl 1): S41.
    [8] Constantinou C, Harrington JC, Cadieux RA, et al. Dosimetry at lung-muscle and lung-bone interface using 6 MV and 10 MV X rays and 12-18 MeV electrons. Int J Radiat Oncol Biol Phys, 1988, 15(suppl 1): 244.
    [9] Constantinou C, Harrington JC, Dewerd LA. An electron density calibration phantom for CT-based treatment planning computers. Med Phys, 1992, 19(2): 325-327. doi: 10.1118/1.596862
    [10] Mackie TR, El-Khatib E, Battista J, et al. Lung dose corrections for 6- and 15-MV x rays. Med Phys, 1985, 12(3): 327-332. doi: 10.1118/1.595691
    [11] Waghorn B, Shah AP, Ngwa W, et al. A computational method for estimating the dosimetric effect of intra-fraction motion on step-and-shoot IMRT and compensator plans. Phys Med Biol, 2010, 55(14): 4187-4202. doi: 10.1088/0031-9155/55/14/015
    [12] 邹华伟, 贾明轩, 吴荣, 等. 调强适形放射治疗计划剂量学的验证. 中华放射医学与防护杂志, 2004, 24(3): 256-258. doi: 10.3760/cma.j.issn.0254-5098.2004.03.025
    [13] 张晓军, 王建华, 涂彧. 二维电离室矩阵在调强适形放疗剂量学验证中的应用. 国际放射医学核医学杂志, 2007, 31(5): 318-320. doi: 10.3760/cma.j.issn.1673-4114.2007.05.019
  • [1] 马建华李明刘宝玲屠婷婷尚玉萍王忠明 . 宫颈癌术后容积旋转调强放疗与5野调强放疗计划的剂量学比较. 国际放射医学核医学杂志, 2018, 42(1): 41-46. doi: 10.3760/cma.j.issn.1673-4114.2018.01.008
    [2] 张晓军王建华涂彧 . 二维电离室矩阵在调强适形放疗剂量学验证中的应用. 国际放射医学核医学杂志, 2007, 31(5): 318-320.
    [3] 苏灿森孙文钊黄晓延 . 两种放疗计划系统用于鼻咽癌螺旋断层调强放疗计划的剂量学比较. 国际放射医学核医学杂志, 2023, 47(8): 468-476. doi: 10.3760/cma.j.cn121381-202206001-00333
    [4] 曹丽媛鞠永健李克新高璇 . PTW 729电离室矩阵不同验证方法用于宫颈癌术后调强放疗计划验证结果分析. 国际放射医学核医学杂志, 2021, 45(10): 631-635. doi: 10.3760/cma.j.cn121381-202104016-00110
    [5] 刘凌湘陈照辉黄贤海易兰 . 同一射野角度不同优化方式对中央型肺癌调强适形放射治疗计划的影响. 国际放射医学核医学杂志, 2021, 45(1): 24-30. doi: 10.3760/cma.j.cn121381-202001016-00002
    [6] 侯朝华霍彬宋杨 . 组织间近距离放射性粒子植入治疗肿瘤的进展. 国际放射医学核医学杂志, 2015, 39(4): 348-351. doi: 10.3760/cma.j.issn.1673-4114.2015.04.016
    [7] 郑玮洪梅常远 . 局部晚期非小细胞肺癌患者调强放疗后发生急性放射性肺损伤的影响因素及风险预测模型研究. 国际放射医学核医学杂志, 2022, 46(9): 536-543. doi: 10.3760/cma.j.cn121381-202112021-00216
    [8] 李彩虹王沛沛李金凯曹远东孙新臣 . 基于Python实现鼻咽癌放疗计划剂量学指标的自动提取和量化评价. 国际放射医学核医学杂志, 2024, 48(): 1-7. doi: 10.3760/cma.j.cn121381-202303018-00404
    [9] 张敏周莉张开贤李苓石存 . 中下段食管癌调强放疗布野方案剂量学比较. 国际放射医学核医学杂志, 2012, 36(2): 121-123, 封三. doi: 10.3760/cma.j.issn.1673-4114.2012.02.016
    [10] 冯静林建海廖绍光骆华春傅志超 . 宫颈癌术后调强放疗中骨髓抑制与骨髓照射剂量体积的关系. 国际放射医学核医学杂志, 2020, 44(3): 143-150. doi: 10.3760/cma.j.cn121381-201811039-00002
    [11] 朱云云陈杰傅志超叶凡骆华春 . 全脑全脊髓照射调强放疗与常规放疗致急性血液学不良反应的比较. 国际放射医学核医学杂志, 2022, 46(7): 441-447. doi: 10.3760/cma.j.cn121381-202110003-00202
  • 加载中
图(2)表(2)
计量
  • 文章访问数:  1924
  • HTML全文浏览量:  561
  • PDF下载量:  3
出版历程
  • 收稿日期:  2010-11-15
  • 刊出日期:  2011-01-25

调强适形放射治疗剂量验证中CT密度值的校准分析

    通讯作者: 吴仁华, rhwu@stu.edu.cn
    通讯作者: 邱庆春, qcqiu@stu.edu.cn
  • 1. 515041 汕头,广东汕头大学医学院附属肿瘤医院放疗科
  • 2. 515041 汕头,广东汕头大学医学院第二附属医院医学影像学教研室
  • 3. 515041 汕头,广东汕头大学医学院医学物理学与信息学教研室

摘要:  目的 基于调强适形放射治疗(IMRT)验证体模,分析放射治疗计划系统中CT密度转换曲线在IMRT剂量验证中的影响,并加以校准。 方法 利用CT电子密度体模刻度计划系统中的CT密度转换曲线,取12例鼻咽癌患者的IMRT计划,分别移植至IMRT验证体模,计算剂量分布,并于剂量均匀处放置测量点,利用电离室测量出其实际剂量。将IMRT验证体模的物理密度与对应的CT值输入计划系统,对CT密度转换曲线进行修正,重新计算剂量,在保持其他计划参数一致的前提下,比较前后两次计算值与实测值的差异。 结果 12例鼻咽癌IMRT计划中,前后两次计算值和实测值的平均误差分别为1.96%±0.87%和0.63%±0.74%,修正后的计算值误差控制在±2%以内,而修正前计算值的最大误差为3.24%。 结论 利用修正后的CT密度转换曲线,计划中的剂量计算值更加贴近实际测量值。为提高IMRT剂量验证的精确性,验证体模在使用之前必须进行CT密度值的校准、确认。

English Abstract

  • 调强适形放射治疗(intensity modulated radiation therapy,IMRT)是20世纪90年代发展起来的一项先进放射治疗技术,它具有靶区的高剂量区三维适形、周围危及器官受量少等优点,能有效提高肿瘤的局部控制率,减少正常组织的损伤,改善患者的生存质量。由于IMRT是一项复杂的技术,为保证其临床治疗的安全可靠,必须对IMRT计划进行准确的剂量验证。

    放射治疗计划系统依据CT电子密度体模确认的CT密度转换曲线,将CT图像中的CT值转换为相应的物理密度(或电子密度),在物理密度(或电子密度)的基础上进行组织不均匀性校正计算,最终得到放射治疗计划的剂量分布[1-2]。然而,由于刻度CT密度转换曲线的CT电子密度体模缺乏与IMRT验证体模密度相接近的组织材料,因此,需要计划系统依据转换曲线的邻近区间,插值计算出IMRT验证体模的物理密度。这插值计算出来的物理密度与实际验证体模的密度存在少许差异,而这些差异将在一定程度上影响IMRT计划验证的精确性。本研究通过实验,定量分析CT密度转换曲线在IMRT剂量验证中的影响,并加以校准。

    • (1)美国CIRS公司Model 062电子密度体模,为27 cm×33 cm×5 cm的椭圆型薄柱体,环氧树脂材料,含有人体等效组织材料9种,几何参数及各插件位置、尺寸、密度参数如图 1表 1所示。

      图  1  美国CIRS公司Model 062电子密度体模及插件示意图

      物理密度(g/cm3) 电子密度(x1023/m3) 相对水电子密度
      体模(水等效) 1.01 3.346 1.002
      肺(吸气) 0.20 0.634 0.190
      肺(呼气) 0.50 1.632 0.489
      脂肪 0.96 3.170 0.949
      乳腺 0.99 3.261 0.976
      肌肉 1.06 3.483 1.043
      1.07 3.516 1.052
      骨(200mg/cm3) 1.16 3.730 1.117
      致密骨(800 mg/cm3) 1.53 4.862 1.456

      表 1  美国CIRS公司Model 062电子密度体模的密度参数表

      (2)德国IBA公司IMRT验证体模,由各种不同尺寸的模块构成,各模块具有规矩的形状,为正方体或长方体。整个体模均由同一种材料RW3(一种等效水材料)组成,主要成分为C8H8,含2.1% ± 0.2%的TiO2

      (3)德国IBA公司的二维电离室矩阵MatriXX,由1020个通气的电离室排列成32×32平面矩阵(四个顶角处无电离室)。电离室尺寸为直径4.5 mm、高5 mm,灵敏体积0.07 cm3,相邻电离室中心间距为7.62 mm,电离室矩阵有效测量面积为24.4 cm×24.4 cm,有效测量点位于上表面下3 mm。每个电离室在出厂时已经用60Co γ射线进行校准,刻度因子以矩阵文件形式给出。

      (4)飞利浦公司的放射治疗计划系统Pinacle,版本为8.0。美国Picker公司的螺旋CT机,型号为PQ5000,定期由工程师负责校准检测。美国Varian公司的医用直线加速器23EX,测量前均检测、校准机器机械参数、输出剂量等,确保机器性能稳定。德国IBA公司的指形电离室FC65-G,灵敏体积0.65 cm3;德国IBA公司的剂量仪DOSE-1, 二者均由广东省辐射剂量计量检定站进行检定。

    • 在螺旋CT机下扫描CT电子密度体模,扫描床板为放疗专用碳素平板,扫描电压为120 kV,电流为200 mA,扫描方式为螺旋扫描: 层厚3 mm,层距3 mm。将图像传至计划系统,勾画出体模中各插件,由计划系统计算出材料的平均CT值。将各等效组织材料的平均CT值与对应的已知物理密度值输入计划系统Pinacle,绘制出相应的CT密度转换曲线。

    • 随机选取12例鼻咽癌患者的静态IMRT计划,计划参数为7个照射野,全周等分,最大子野数为90个,子野数范围为66~90个,计划设计中滤去面积小于4 cm2以及机器跳数小于5 MU(计划系统、加速器上使用的剂量计算单位,1 MU=1 cGy)的小子野。

    • 在同一CT机下,以相同的扫描条件扫描IMRT验证体模,图像传输至治疗计划系统,三维重建该验证体模。将12例鼻咽癌IMRT计划分别移植至IMRT验证体模中,按照计划的实际照射角度,计算剂量,剂量以机器跳数归一,于靶区剂量均匀处(剂量梯度变化小于±6%)放置测量点。

      剂量验证在加速器23EX上进行,测量的射线为6 MV光子线,该机器的剂量输出特性采用美国医学物理师协会第51号报告中的规范预先进行了校准[3]。将IMRT验证体模置于加速器治疗床上,十字线对准体模中心,利用激光灯和数字测距尺进行精确摆位。于测量点处放置FC65-G指形电离室,连接电缆至剂量仪DOSE-1,测量前给予预照射100 cGy,使其达到电离平衡,预照射后进行温度、气压校正和本底测量。每一计划均重复测量3次,测量结果取平均值。

      同时,针对12例鼻咽癌IMRT计划,利用MatriXX进行平面相对剂量验证,将鼻咽癌IMRT计划的射野归一到0°机架角度,按设计进行照射、测量。利用二维电离室矩阵的附带软件对测量与计算的平面剂量进行比较,采用γ指数值的复合判断方法,以剂量点距离3 mm、相差3%为条件,以相对剂量方式进行,剂量以射野高剂量坪区归一,0≤γ指数≤1的点所占比例为通过率。为了避免低剂量区测量和计算不准确的影响,以及射野外测量点不参与比较的考虑,设置实测最高剂量点的10%做为阈值,测量点实测剂量低于此值者不进行比较。

    • 在治疗计划系统中勾画出IMRT验证体模材料,由计划系统计算出其平均CT值。将验证体模的已知物理密度与平均CT值输入计划系统,对CT密度转换曲线进行修正,使得计划系统中验证体模图像的物理密度修正为对应的实际物理密度,重新计算剂量。在保持机器跳数、子野数、射野分布、子野分布等计划参数一致的前提下,比较前后两次计划计算值与实测值之间的差异,数据利用SPSS 16.0进行统计分析,修正前后之间的误差比较采用配对t检验方法,以P < 0.05为差异有统计学意义。

    • 由电子密度体模各等效组织材料的平均CT值与对应的已知物理密度值输入计划系统Pinacle后绘制出的CT密度转换曲线见图 2。依据CT密度曲线,插值计算出的IMRT验证体模的物理密度约为1.020 g/cm3,而厂家说明书上提供的标准密度为1.045 g/cm3,二者约有0.025 g/cm3的差异。

      图  2  计划系统中CT密度转换曲线(Pinacle中CT值=CT图像的CT值+1000) CT图像的CT值+1000)

      参照多数医院使用二维电离室矩阵验证IMRT计划的方法,以γ指数通过率 > 90%为照射野剂量分布验证通过,则12例鼻咽癌IMRT计划验证测得的γ指数通过率为95.34%±1.50%,所有计划的通过率均 > 90%。由于IMRT平面相对剂量的验证不涉及剂量计算上具体的绝对值,因此,验证体模材料的微小密度差异并未给平面相对剂量的验证带来误差。

      12例鼻咽癌IMRT计划分别按计划条件在加速器上模拟治疗方式的实际出束照射进行测量,将测量点处的计划剂量值与测得的靶区剂量值进行比较,修正前和修正后的百分相对误差全部控制在4%以内,符合临床要求。表 2给出了12例鼻咽癌IMRT计划前后两次计算值与实测值差异的统计分析,其中,误差=(计算值-测量值)÷测量值× 100%。由表 2数据可知,在同一个计划中,修正后的计算值比修正前低3 cGy左右,修正后的计算值与测量结果比较接近,误差控制在±2%以内,而修正前的计算值与测量结果的误差在0.79%~3.24%之间,修正后的误差显著低于修正前(t=19.256, P < 0.05)。12例鼻咽癌IMRT计划中,前后两次计算值与实测值的平均误差分别为1.96% ± 0.87%和0.63% ± 0.74%,最大误差分别为3.24%和1.81%。IMRT计划中验证体模微小的密度差异0.025 g/cm3,导致了绝对剂量验证中出现较大的偏差。

      实测值(cGy) 计算值[cGy (误差,%)]
      修正前 修正后
      例1 210.2 217(3.24) 214(1.81)
      例2 207.4 211(1.74) 208(0.29)
      例3 203.4 205(0.79) 202(-0.69)
      例4 194.2 196(0.93) 194(-0.10)
      例5 202.4 206(1.78) 204(0.79)
      例6 203.1 205(0.94) 203 (-0.05)
      例7 200.2 205(2.40) 202(0.90)
      例8 199.9 204(2.05) 201(0.55)
      例9 196.8 203(3.15) 200(1.63)
      例10 193.6 196(1.24) 194(0.21)
      例11 192.1 197(2.55) 194(0.99)
      例12 202.5 208(2.72) 205(1.23)

      表 2  12例鼻咽癌IMRT计划前后两次计算值与实测值的比较分析

    • 作为现代精确的放疗技术,IMRT有着较高的剂量变化梯度和较好的肿瘤靶区剂量适形度等特点。但由于IMRT是一个相当复杂的治疗技术,需要严格的质量保证措施,治疗计划设计和实施过程的微小误差都可能给计划的实施带来较大的差异,从而影响治疗的疗效,甚至导致治疗的失败。国际辐射单位和测量委员会第24号报告指出,原发灶根治剂量的精确性应好于5%,靶区剂量偏离最佳剂量5%时,就有可能使原发灶肿瘤失控或并发症增加,从而可能导致治疗计划的失败[4]。因此,每例患者的IMRT计划执行前的剂量实际测量验证是IMRT质量保证工作中至关重要的一环,以确保由治疗计划系统精确设计和计算的剂量分布能够在治疗设备上准确执行。

      CT扫描图像中的CT值,反映了X射线在人体不同密度组织的线性衰减关系,而实际应用中的设备受成像理论和设备本身的限制,这两个参数之间并非是理想的线性关系[1-2]。对于CT值线性,杨克柽等[5]提出,用二阶拟合的2-范数(即:矩阵A的转置乘以A矩阵特征根最大值的开根号)来表述,从统计学的角度,验证了二阶拟合曲线更符合物质的CT值与其对应的X射线衰减系数的分布关系。在放射治疗计划系统中,CT值与物质密度之间需要建立CT密度转换曲线,之后可由患者图像CT值确定患者不同组织的密度分布情况,从而进行组织不均匀性校正等剂量分布计算[6-7]。研究发现,由CT密度转换错误导致的计划系统剂量计算误差可能十分显著,如采用6 MV光子线照射100 cGy的剂量,CT密度转换错误的计划所计算出的剂量分布热点可比正确值高出20.1%[8-10]

      IMRT的计划和实施是一项复杂的系统工作,其质量保证是一个必须重视和急需解决的问题[4]。利用体模进行IMRT的绝对剂量验证,是目前国内外普遍采用的基本剂量验证方法[11-13]。作为平面剂量测量仪,MatriXX在IMRT绝对剂量测量方面不如单个电离室有优势,如果采用MatriXX做多角度平面绝对剂量验证,其电离室须进行系数校正,且存在方向性响应和结构影响,这就加大了测量的误差。所以,在绝对剂量测量中,为了提高测量的准确性,大部分使用单个电离室在体模中模拟测量的方法[13]

      由于刻度CT密度曲线的电子密度体模缺乏与IBA IMRT验证体模密度相接近的组织材料,计划系统依赖邻近区间插值计算出的体模密度与实际验证体模密度有些微小差异,约0.025 g/cm3。放射治疗计划系统Pinacle基于物理密度来计算剂量分布,物理密度上的差异将影响到绝对剂量计算的误差。IMRT平面相对剂量验证是验证电离室探头平面的相对剂量分布,不涉及绝对剂量值,因此,体模的微小密度差异没有给其带来测量误差影响。不过,该密度差异会给点绝对剂量的验证带来一定的影响。为了更全面地反映CT值与物质密度之间的关系,必须对目前所普遍采用的电子密度体模做适当的改进,增加目标物质的数目。

      IMRT验证体模在使用之前,必须进行CT密度值的校准确认。在静态IMRT计划中,子野数目多,有着较小的MU机器跳数和较小面积的子野,IMRT验证体模整体的微小密度偏差可引起较大的剂量误差,利用IMRT验证体模的物理密度及CT值对CT密度曲线进行修正,可以减少验证体模密度差异所带来的系统误差。利用修正后的曲线,计算计划剂量,其计算值更加贴近实际测量数值,能够更好地为IMRT技术提供可靠的保障。

参考文献 (13)

目录

    /

    返回文章
    返回