-
选用Ashland公司生产的GafChromic® EBT3-0810胶片(批号:05061904、04022001),裁剪为4 cm×4 cm的正方形并在胶片角落标记方向、序号及剂量值。使用武威肿瘤医院医用重离子加速器提供的标称260 MeV/u碳离子笔形束,通过XY正交扫描磁铁均匀扫描成15 cm×15 cm的照射野进行碳离子束辐照。不同LET的碳离子束通过水等效系数为1.16的聚甲基丙烯酸甲酯 (polymethyl methacrylate, PMMA)降能获得,LET与PMMA降能片厚度的对应关系见表1。EBT3胶片辐照前,使用绝对剂量电离室(Advanced Markus 34045,德国PTW公司)获得等中心位置各LET对应的监测电离室刻度因子。利用表1中所列LET的碳离子束分别对EBT3胶片辐照0.5,1.0,1.5,2.0,3.0 Gy,用以获得相应剂量下不同LET的RE值。随后利用LET为16.21、21.99, 64.68,83.68 keV/μm的碳离子束分别对EBT3胶片辐照1,2,4,6,8,10 Gy,用来对比大剂量区间内剂量刻度曲线公式的拟合结果。最后选用16.21和64.68 keV/μm两种LET的碳离子束按1:1,1:2, 1:3, 2:1, 3:1的剂量比例混合,分别对EBT3胶片辐照1,2,3,4,6 Gy。所有胶片辐照过程中均采用纸胶带固定于PMMA,胶片中心点与等中心点重合,照射时胶片后方PMMA厚度大于5 cm,前方PMMA厚度由LET值确定。
表 1 LET对应的PMMA降能片厚度
LET/(keV·μm−1) PMMA降能片厚度/mm 16.21 26.41 21.99 77.95 33.31 97.62 42.52 102.21 52.36 106.25 64.68 108.69 83.68 109.98 为了获得不同LET碳离子的RE值,本工作使用了甘肃省肿瘤医院的Varian 600CD医用直线加速器提供的6 MV光子进行EBT3胶片的剂量刻度。加速器校准后,对EBT3胶片分别辐照0.25, 0.50, 1.00, 1.50, 2.00, 3.00, 4.00, 5.00 Gy。所有经碳离子或光子辐照的和对照(未经辐照)的EBT3胶片,避光保存24 h后进行数字化扫描。
-
采用EPSON Expression 12000XL扫描仪的肖像模式,即胶片原始长边平行于扫描窗长边进行扫描。胶片扫描前将扫描仪预热10 min,每张胶片连续扫描5次,取后3次保存文件。扫描时采用48位RGB模式,关闭所有图像增强功能,图片分辨率设置为75 dpi,扫描后的文件保存为TIFF(Tag Image File Format)格式[17]。
EBT3胶片的吸收光谱在595和633 nm处出现峰值,分别对应绿光和红光通道,但在剂量测量范围内红光通道的吸收峰明显高于绿光通道的吸收峰 [17],因此分析红光通道的数据可以提高测量灵敏度。使用FilmQA Pro软件(美国Ashland公司)读取EBT3胶片TIFF图片的像素值,将感兴趣区域(Region of Interest, ROI)设定为胶片中心面积约为4 cm2大小的圆,获得相同胶片3次扫描结果的红光通道像素平均值及标准差。
采用Martisikova 等[4]报道的方法,通过式(1)将像素平均值转换为净光密度值(netOD),并通过标准差传递公式(2)计算netOD标准差[18]。
$$ \overline{\mathrm{netOD}}={\mathrm{log}}_{10}\left(\frac{{\overline{\mathrm{PV}}}_{\mathrm{unexp}}-{\overline{\mathrm{PV}}}_{\mathrm{bckg}}}{{\overline{\mathrm{PV}}}_{\mathrm{exp}}-{\overline{\mathrm{PV}}}_{\mathrm{bckg}}}\right) , $$ (1) $$ {\sigma }^{}_{\overline{\mathrm{netOD}}}=\frac{1}{\mathrm{ln}\left(10\right)}\sqrt{\left[\frac{{\sigma }_{{\overline{\mathrm{PV}}}_{\mathrm{unexp}}+}^{2}{\sigma }_{{\overline{\mathrm{PV}}}_{\mathrm{bckg}}}^{2}}{{{(\overline{\mathrm{PV}}}_{\mathrm{unexp}}-{\overline{\mathrm{PV}}}_{\mathrm{bckg}})}^{2}}\right]+\left[\frac{{\sigma }_{{\overline{\mathrm{PV}}}_{\mathrm{exp}}+}^{2}{\sigma }_{{\overline{\mathrm{PV}}}_{\mathrm{bckg}}}^{2}}{{{(\overline{\mathrm{PV}}}_{\mathrm{exp}}-{\overline{\mathrm{PV}}}_{\mathrm{bckg}})}^{2}}\right]} , $$ (2) 其中:
${\overline{\mathrm{PV}}}_{\mathrm{unex}\mathrm{p}}$ 、${\overline{\mathrm{PV}}}_{\mathrm{exp}}$ 、${\overline{\mathrm{PV}}}_{\text{bckg}}$ 分别是未辐照胶片、辐照后胶片和黑色不透光板的红光通道像素平均值。${\sigma }^{}_{{\overline{\mathrm{PV}}}_{\mathrm{unexp}}}$ 、${\sigma }^{}_{{\overline{\mathrm{PV}}}_{\mathrm{exp}}}$ 、${\sigma }^{}_{{\overline{\mathrm{PV}}}_{\mathrm{bckg}}}$ 分别是未辐照胶片、辐照后胶片和黑色不透光板红光通道像素平均值的标准差,${\sigma }^{}_{\overline{\mathrm{netOD}}}$ 是计算得到的netOD标准差。剂量刻度曲线是EBT3胶片剂量学特性的主要表现形式,获得剂量刻度曲线的主要方法是通过公式拟合。采用不同的拟合公式会影响剂量刻度曲线的精度并最终影响EBT3胶片的剂量转换结果。Devic等[1,3]和Castriconi等[12]使用了依赖二至三个拟合参数的公式来拟合剂量刻度曲线。本工作使用式(3)、(4)以及式(4)中参数C分别等于2和3时的平方和立方公式对胶片物理吸收剂量与净光密度值的关系进行拟合,式中A、B、C为待拟合参数:
$$ D\left(\mathrm{net}\mathrm{OD}\right)=A\times \mathrm{net}\mathrm{OD}\times {\mathrm{exp}}^{B\times \mathrm{net}\mathrm{OD}} \text{,} $$ (3) $$ D\left(\mathrm{net}\mathrm{OD}\right)=A\times \mathrm{net}\mathrm{OD}+{B}\times \mathrm{net}\mathrm{OD}^{C} 。 $$ (4) -
Castriconi等[12]、Grilj等[9]、Shunsuke等[16]使用RE来量化EBT3胶片随LET的剂量欠响应。本工作RE定义为EBT3胶片达到相同净光密度值所需光子剂量与碳离子剂量的比值,即式(5),式中
$ {D}_{\mathrm{r}} $ 、$ {D}_{\mathrm{C}-\mathrm{i}\mathrm{o}\mathrm{n}} $ 分别为光子和碳离子剂量,$$ \mathrm{RE}={\left.\frac{{D}_{\mathrm{r}}}{{D}_{\mathrm{C-ion}}}\right|}_{\mathrm{net}\mathrm{OD}} 。 $$ (5) RE剂量修正方法是根据混合LET辐照中各LET的剂量比例得到混合LET辐照等效的剂量平均LET,借助RE与LET的关系得到等效LET的RE,再依据混合LET辐照后胶片的净光密度值得到等效的光子剂量,最后通过RE和光子剂量推导得到RE修正剂量。
拟合参数剂量修正方法是通过混合LET辐照后刻度曲线中拟合参数随LET所占比例的变化规律,得到混合LET辐照后相应比例的剂量刻度曲线拟合参数,再依据胶片净光密度值得到胶片的拟合参数修正剂量。
-
选择16.21, 64.68 keV/μm和光子大剂量区间的胶片刻度数据来展示剂量刻度曲线公式的拟合结果。图1中(a)、(b)、(c)分别为16.21, 64.68 keV/μm和光子在红光通道内的拟合结果,拟合优度(R2)汇总于表2。结果表明,式(4)的平方拟合结果最差,其余三种拟合的拟合优度均高于0.999 0。16.21 , 64.68 keV/μm和光子的剂量刻度曲线拟合时,式(4)的拟合精度最高,然而式(4)的拟合结果波动较大且拥有三个拟合参数,不利于后续分析。Castriconi等[12]指出在选择拟合公式时需综合考虑实验不确定性和拟合误差情况,虽然式(4)的立方拟合相较于式(3)拥有更低的拟合波动,但式(3)在拟合优度和拟合精度间获得了平衡,最终选择式(3)来获得剂量刻度曲线。
表 2 剂量刻度曲线公式的拟合结果
拟合公式 拟合参数 拟合
优度R2类别 A B C D=A×netOD+B×netOD2 1.26 39.31 2.00 0.997 5 16.21 keV/μm 4.32 47.18 2.00 0.998 0 64.68 keV/μm 4.60 26.42 2.00 0.998 8 6 MV光子 D=A×netOD+B×netOD3 7.52 56.39 3.00 0.999 9 16.21 keV/μm 10.61 80.05 3.00 0.999 9 64.68 keV/μm 7.54 52.98 3.00 0.999 6 6 MV光子 D=A×netOD+B×netODC 7.07 53.04 2.87 0.999 9 16.21 keV/μm 10.24 74.76 2.90 0.999 9 64.68 keV/μm 6.94 40.78 2.67 0.999 7 6 MV光子 D=A×netOD×expB×netOD 5.72 2.66 - 0.999 9 16.21 keV/μm 8.14 2.63 0.999 9 64.68 keV/μm 6.03 2.40 0.999 7 6 MV光子 红光通道内拟合优度相差无几也表明实验数据与拟合曲线变化趋势符合的很好。相同剂量照射后,64.68 keV/μm胶片的净光密度值小于16.21 keV/μm的净光密度值,表明EBT3胶片的剂量响应依赖于LET。任意LET下刻度的胶片,其净光密度值随剂量增大,表明EBT3胶片的剂量响应受到物理吸收剂量影响。这与Kawashima等[10],Wang等[11]以及Shunsuke等[16]报道的EBT3胶片剂量响应依赖LET与剂量双重作用的结果相符合。
-
在使用式(5)获得各LET下1,2,3 Gy的RE后,对三个剂量下的RE计算平均值以减小剂量对RE的影响。绘制了低LET下(16.21
$ \leqslant $ LET$ \leqslant $ 83.68 keV/μm)RE与LET的关系如图2所示,并对RE与LET间的关系进行了线性拟合。16.21和64.68 keV/μm按照2:1,3:1,1:1,1:2,1:3的剂量比例混合辐照后等效的剂量平均LET分别为28.33,32.38,40.45,48.52,52.56 keV/μm,由拟合关系可以得到等效LET下的RE,结合胶片的净光密度值得到等效的光子剂量后计算混合LET辐照的RE修正剂量,混合辐照后1 Gy的RE修正剂量汇总于表3,结果表明RE剂量修正方法得到的剂量偏差在10%以内。表 3 混合LET辐照后的RE与RE修正剂量
照射
总剂量/Gy混合
比例剂量平均
LET/(keV·μm−1)光子等效
剂量/GyRE RE修正
剂量/Gy剂量
偏差/Gy1.00 3:1 28.33 0.94 0.88 1.06 −0.06 1.00 2:1 32.37 0.90 0.86 1.04 −0.04 1.00 1:1 40.45 0.74 0.82 0.91 0.09 1.00 1:2 48.52 0.73 0.77 0.94 0.06 1.00 1:3 52.56 0.69 0.75 0.92 0.08 注:剂量偏差=照射总剂量−RE修正剂量 -
混合LET辐照的结果如图3所示,图3中黑色、红色、绿色、蓝色、青色、紫色、橙色实线分别代表16.21 keV/μm, 3:1,2:1,1:1,1:2,1:3剂量比例以及64.68 keV/μm辐照后,红光通道内胶片响应使用式(3)拟合的结果,拟合R2均高于0.990 0。结果表明两种LET混合辐照后,胶片物理吸收剂量与净光密度值的拟合曲线与等效LET的剂量刻度曲线相似,这与Kawashima等[10]通过改变LET和剂量可以获得相同净光密度值的结果相符。
在16.21 keV/μm与64.68 keV/μm的3:1,2:1,1:2,1:3比例混合辐照中,16.21 keV/μm所贡献剂量占混合辐照胶片总剂量的比例分别为0.75,0.67,0.33,0.25。16.21和64.68 keV/μm单独辐照时16.21 keV/μm贡献比例分别为1和0。绘制式(3)拟合后参数A、B随16.21 keV/μm所占剂量比例的变化关系如图4所示,图中黑色实线和红色虚线分别为拟合参数A、B的线性拟合结果。
根据图4中的线性拟合结果,获得了16.21与64.68 keV/μm在1:1比例混合辐照下剂量刻度曲线对应的拟合参数A、B值,结合混合辐照后胶片的净光密度值,得到了拟合参数修正剂量,结果汇总于表4,表明拟合参数修正方法得到的剂量偏差在5%以内,相比于RE剂量修正方法提升显著。
表 4 1:1比例混合辐照的拟合参数剂量修正结果
照射
总剂量/Gy净光密度值 拟合参数 拟合参数
修正剂量/Gy剂量
偏差/GyA B 1.00 0.099 6 7.49 2.54 0.96 0.04 2.00 0.168 2 1.93 0.07 3.00 0.222 0 2.93 0.07 4.00 0.266 7 3.94 0.06 6.00 0.337 0 5.94 0.06 注:剂量偏差=照射总剂量−拟合参数修正剂量。 -
RE剂量修正方法得到的剂量偏差在10%以内,造成这种现象的原因是混合LET辐照后EBT3胶片的剂量响应不完全和等效LET的剂量响应相同,EBT3胶片在碳离子束辐照下的剂量欠响应随着LET的增大而增大,混合LET辐照时较大LET发生的剂量欠响应要高于等效LET辐照时发生的欠响应。这会导致混合LET辐照后胶片净光密度值的变化,相应的等效光子剂量也不同于等效LET得到的光子剂量,从而造成RE剂量修正方法得到的剂量偏差较大,这也说明EBT3胶片碳离子束混合辐照下的剂量响应受到LET的影响。
通过拟合参数随LET所占比例修正方法得到的剂量偏差在5%以内,这种方法在获得混合LET辐照后剂量刻度曲线的拟合参数时考虑了LET所造成的胶片剂量欠响应。在碳离子束放射治疗中,由脊形过滤器的构造原理可模拟得到束流中各LET值所占据的比例,这也使得通过混合LET辐照后拟合参数与LET所占比例的关系,获得任意比例辐照后的拟合参数变成可能。基于此,本文提出的拟合参数剂量修正方法具有特殊的优势。
Dose Correction Methods of EBT3 Radiochromic Film Irradiated by Mixed-LET Carbon-ion Beams
-
摘要: 针对EBT3辐射变色胶片对碳离子束混合LET辐照的剂量欠响应比较了两种剂量修正方法。利用260 MeV/u的碳离子束通过被动降能得到多种剂量平均LET的碳离子束,利用这些碳离子束进行了胶片剂量响应刻度辐照,选择最佳的拟合公式得到了胶片剂量刻度曲线。使用RE(Relative Efficiency)量化了EBT3胶片随LET的剂量欠响应,并使用RE剂量修正法修正了混合LET辐照胶片的剂量。此外,根据剂量刻度曲线公式中拟合参数随不同LET所占剂量比例的变化规律,提出了拟合参数剂量修正法并修正了混合LET辐照胶片的剂量。最后比较了这两种方法的结果,表明拟合参数方法得到的剂量偏差在5%以内,优于RE方法10%以内的剂量偏差。
-
关键词:
- EBT3辐射变色胶片 /
- 碳离子束 /
- 剂量刻度曲线 /
- 混合LET辐照 /
- 剂量修正方法
Abstract: Two dose correction methods were compared for the under-response of EBT3 radiochromic films to the irradiation of mixed-LET carbon-ion beams. Carbon-ion beam of 260 MeV/u was used to obtain carbon-ion beams with multiple dose-averaged LETs by energy degraders. The films were irradiated by these carbon ion beams, and the optimal fitting formula was selected to fit the film dose response calibration curves. Relative Efficiency(RE) was used to quantify the under-response of EBT3 film with LET, and the doses of the films irradiated by mixed-LET carbon ions were corrected by the RE method. In addition, according to the relation between the fitting parameters of the dose response calibration curve and the dose proportions of different LETs, a new dose correction method was proposed and the doses of the films irradiated by mixed-LET irradiation were also corrected. Finally, the results of the two methods are compared, and it is shown that the dose deviation obtained by the fitting parameter method is less than 5%, which is better than the dose deviation obtained by the RE method less than 10%. -
表 1 LET对应的PMMA降能片厚度
LET/(keV·μm−1) PMMA降能片厚度/mm 16.21 26.41 21.99 77.95 33.31 97.62 42.52 102.21 52.36 106.25 64.68 108.69 83.68 109.98 表 2 剂量刻度曲线公式的拟合结果
拟合公式 拟合参数 拟合
优度R2类别 A B C D=A×netOD+B×netOD2 1.26 39.31 2.00 0.997 5 16.21 keV/μm 4.32 47.18 2.00 0.998 0 64.68 keV/μm 4.60 26.42 2.00 0.998 8 6 MV光子 D=A×netOD+B×netOD3 7.52 56.39 3.00 0.999 9 16.21 keV/μm 10.61 80.05 3.00 0.999 9 64.68 keV/μm 7.54 52.98 3.00 0.999 6 6 MV光子 D=A×netOD+B×netODC 7.07 53.04 2.87 0.999 9 16.21 keV/μm 10.24 74.76 2.90 0.999 9 64.68 keV/μm 6.94 40.78 2.67 0.999 7 6 MV光子 D=A×netOD×expB×netOD 5.72 2.66 - 0.999 9 16.21 keV/μm 8.14 2.63 0.999 9 64.68 keV/μm 6.03 2.40 0.999 7 6 MV光子 表 3 混合LET辐照后的RE与RE修正剂量
照射
总剂量/Gy混合
比例剂量平均
LET/(keV·μm−1)光子等效
剂量/GyRE RE修正
剂量/Gy剂量
偏差/Gy1.00 3:1 28.33 0.94 0.88 1.06 −0.06 1.00 2:1 32.37 0.90 0.86 1.04 −0.04 1.00 1:1 40.45 0.74 0.82 0.91 0.09 1.00 1:2 48.52 0.73 0.77 0.94 0.06 1.00 1:3 52.56 0.69 0.75 0.92 0.08 注:剂量偏差=照射总剂量−RE修正剂量 表 4 1:1比例混合辐照的拟合参数剂量修正结果
照射
总剂量/Gy净光密度值 拟合参数 拟合参数
修正剂量/Gy剂量
偏差/GyA B 1.00 0.099 6 7.49 2.54 0.96 0.04 2.00 0.168 2 1.93 0.07 3.00 0.222 0 2.93 0.07 4.00 0.266 7 3.94 0.06 6.00 0.337 0 5.94 0.06 注:剂量偏差=照射总剂量−拟合参数修正剂量。 -
[1] DEVIC S, SEUNTJENS J, SHAM E, et al. Medical Physics, 2005, 32(7): 2245. doi: 10.1118/1.1929253 [2] DEVIC S, TOMIC N, SOARES C G, et al. Medical Physics, 2009, 36(2): 429. doi: 10.1118/1.3049597 [3] DEVIC S, TOMIC N, LEWIS D. Physica Medica, 2016, 32(4): 541. doi: 10.1016/j.ejmp.2016.02.008 [4] MARTISIKOVA M, ACKERMANN B, JAKEL O. Physics in Medicine and Biology, 2008, 53(24): 7013. doi: 10.1088/0031-9155/53/24/001 [5] ALDELAIJAN S, ALZORKANY F, MOFTAH B, et al. Physica Medica, 2016, 32(1): 202. doi: 10.1016/j.ejmp.2015.12.004 [6] MARTISIKOVA M, ACKERMANN B, KLEMM S, et al. Nucl Instr and Meth A, 2008, 591(1): 171. doi: 10.1016/j.nima.2008.03.088 [7] KHACHONKHAM S, DREINDL R, HEILEMANN G, et al. Physics in Medicine and Biology, 2018, 63(6): 065007. doi: 10.1088/1361-6560/aab1ee [8] MARTISIKOVA M, JAKEL O. Radiation Measurements, 2010, 45(10): 1268. doi: 10.1016/j.radmeas.2010.05.028 [9] GRILJ V, BRENNER D J. Physics in Medicine and Biology, 2018, 63(24): 245021. doi: 10.1088/1361-6560/aaf34a [10] KAWASHIMA M, MATSΜMURA A, SOUDA H, et al. Physics in Medicine and Biology, 2020, 65(12): 125002. doi: 10.1088/1361-6560/ab8bf3 [11] 王巍伟, 盛尹祥子, 黄志杰, 等. 中华放射医学与防护杂志, 2018, 38(9): 705. doi: 10.3760/cma.j.issn.0254-5098.2018.09.013 WANG W W, SHENG Y X Z, HUANG Z J, et al. Chinese Journal of Radiological Medicine and Protection, 2018, 38(9): 705. (in Chinese) doi: 10.3760/cma.j.issn.0254-5098.2018.09.013 [12] CASTRICONI R, CIOCCA M, MIRAANDOLA A, et al. Physics in Medicine and Biology, 2017, 62(2): 377. doi: 10.1088/1361-6560/aa5078 [13] FIORINI F, KIRBY D, THOMPSON J, et al. Physica Medica, 2014, 30(4): 454. doi: 10.1016/j.ejmp.2013.12.006 [14] GAMBARINI G, BETTEGA D, CAMONI G, et al. Radiation Physics and Chemistry, 2019, 155: 138. doi: 10.1016/j.radphyschem.2018.06.019 [15] ZHAO L, DAS I J. Physics in Medicine and Biology, 2010, 55(18): 291. doi: 10.1088/0031-9155/55/10/N04 [16] SHUNSUKE Y, CHINATSU A, KAORU S, et al. Radiation Protection Dosimetry, 2018, 180(1/4): 314. doi: 10.1093/rpd/ncy006 [17] NIROOMAND-RAD A, CHIU-TSAO S T, GRAMS M P, et al. Medical Physics, 2020, 47(12): 5986. doi: 10.1002/mp.14497 [18] VALLIERES S, BIENVENUE C, VLADE P P, et al. Review of Scientific Instruments, 2019, 90(8): 083301. doi: 10.1063/1.5109644 [19] 葛宁, 韩栋梁, 辜石勇. 中国医学物理学杂志, 2015, 32(4): 595. doi: CNKI:SUN:YXWZ.0.2015-04-031 GE N, HAN D L, GU S Y. Chinese Journal of Medical Physics, 2015, 32(4): 595. (in Chinese) doi: CNKI:SUN:YXWZ.0.2015-04-031