胡海芹 陳心得 張懷文 洪 潮

?

不同組織密度對IMRT劑量驗證通過率的影響

胡海芹 陳心得 張懷文 洪 潮

目的 研究不同的組織密度對IMRT計劃驗證通過率的影響，為臨床調(diào)強計劃設(shè)計優(yōu)化提供更好的經(jīng)驗。方法分別建立頭部(120 KV，380 mAs)、胸部(140 KV，100 mAs)和參考均勻模體(120 KV，130 mAs)3種模式下CT電子密度表，選取20例鼻咽癌IMRT計劃分別使用這3種CT電子密度表設(shè)計3組計劃方案。結(jié)果在OmniPro I'mRT軟件中設(shè)置3 mm、3%、閾值5為分析限制條件,此20例鼻咽癌IMRT計劃的head-head通過率均達到96.44%±2.3%；head-chest通過率與head-head通過率比較得到為83.78%±7.1%(P<0.05)；head-water通過率與head-head通過率比較得到為83.26%±7%(P<0.05)。結(jié)論對于IMRT計劃優(yōu)化前選取正確的電子密度是實現(xiàn)最大劑量通過率的保證。

組織密度；調(diào)強優(yōu)化；Gamma通過率

(ThePracticalJournalofCancer,2016,31:1448～1452)

近10年來，隨著放療技術(shù)的迅猛發(fā)展，調(diào)強放射治療(IMRT)技術(shù)在全國大多數(shù)醫(yī)院已得到了廣泛運用，IMRT和傳統(tǒng)的三維適形放療(3D-CRT)相比，它在進一步提高腫瘤劑量的同時能更好地保護周圍的正常組織和器官[1]，顯著降低患者的放療副反應(yīng)，提高治療效果。由于IMRT技術(shù)對加速器出束劑量的準確性有極高的要求，且逆向計算的優(yōu)化算法在某些方面還不成熟，放射治療中還存在眾多不確定因素，因此治療前的劑量分布驗證是確保治療準確的關(guān)鍵步驟[2-3]。只有當一個計劃中每個射野分別在模體中實測的劑量分布和計劃系統(tǒng)中計算所得劑量分布相符合，才能把該計劃運用于實際患者的治療。因此實施個體化的調(diào)強驗證顯得尤為重要[4-6]。本文主要通過MatriXX劑量分布的Gamma比較方法對不同組織密度的IMRT計劃通過率影響的探討。

1 資料與方法

1.1 材料

德國西門子大孔徑CT，醫(yī)科達Precise加速器，Pinnacle39.8治療計劃系統(tǒng)，IBA公司二維電離室矩陣MatriXX和OmniPro I'mRT調(diào)強驗證分析軟件，CT電子密度測量模體CatphanTM500。

CT值大小理論上決定于組織的物理密度與掃描條件(X射線的KV值和mAs)。本文使用西門子CT機分別采集頭部(120KV，380 mAs)、胸部(140KV，100 mAs)和參考均勻模體(120KV，130 mAs)3種模式掃描CatphanTM 500(該模體包括4個檢測模塊，帶有直徑15 cm,厚2.5 cm已知密度不同的圓柱形模擬組織插件)。通過Pinnacle39.8 TPS軟件采集所有的模擬組織CT值，從而建立起CT值與相對電子密度的一一對應(yīng)關(guān)系，將其輸入到Pinnacle39.8計劃系統(tǒng)存為CT值轉(zhuǎn)換為相對電子密度表，見圖1，將這3種電子密度表分別命名：head,chest和water。

圖1 CT值轉(zhuǎn)換相對電子密度表

將二維電離室矩陣MatriXX上放置5 cm固體水模體平放于西門子CT床上，掃描間距為3 mm，掃描整個MatriXX。將獲得的影像通過Dicom傳輸?shù)絇innacle39.8TPS上保存為QA模體。

1.2 計劃設(shè)計和QA

選取20例鼻咽癌病例，所有靶區(qū)和危及器官勾劃均參考鼻咽癌調(diào)強放射治療靶區(qū)勾畫-中國醫(yī)學科學院腫瘤醫(yī)院經(jīng)驗總結(jié)[7]，每一位患者選取頭部電子密度表設(shè)計IMRT計劃。物理約束條件包括最小劑量、最大劑量、最小劑量體積直方圖、最大劑量體積直方圖、均勻劑量。在Pinnacle39.8調(diào)強優(yōu)化中限制最大子野數(shù)(≤80)、最小子野面積和機器跳數(shù)(≤7 cm2,≤7MU)，以保證QA通過率最大化。

IMRT計劃完成后，副主任放療醫(yī)師確認所有臨床要求都能滿足。再選擇胸部和參考均勻模體的電子密度重新對批準的IMRT計劃進行后續(xù)優(yōu)化(確保每個照射野子野形狀和分布不變)，產(chǎn)生2組新的計劃。

分別將3組IMRT計劃導入QA模體中，將3組IMRT計劃照射野機架角歸0°，在Pinnacle39.8系統(tǒng)中采用Adaptive Convolve算法進行劑量計算，并以SPD(源片距)為100 cm、2 mm計算間隔和ASCALL碼分別輸出三個平面劑量分布用于后續(xù)分析，輸出的平面劑量為每個照射野單獨劑量通量分布。分別將這三組電子密度表生成的IMRT QA計劃命名為QA_head、QA_chest和QA_water,將QA_head通過Dicom傳輸?shù)结t(yī)科達Precise加速器。

結(jié)合MatriXX二維電離室矩陣的能響特性，在首次采集錢對MatriXX進行15 min預(yù)熱，然后用劑量率為400 MU/min的6MV能量的X射線進行4 Gy照射后進行本底測量。通過MatriXX采集QA_head每個照射野的平面劑量，將實際采集的平面劑量放于MatriXX Dataset1中，命名為：CA_head,并以此實測劑量作為參考值。Dataset 2為Pinnacle39.8計劃系統(tǒng)中QA_head的平面劑量。按照1VS2運用Gamma分析法得出head-head的通過率。選擇QA_chest和QA_water兩組Pinnacle39.8計劃系統(tǒng)中的平面劑量分別在MatriXX窗口2中打開，同樣按照1VS2運用Gamma分析法得出head-chest和head-water的通過率。

1.3 Gamma方法

Gamma方法被用到測量分布和計劃系統(tǒng)計算分布的比較，也被用于兩種測量分布或者兩種計算分布之間的比較，其表達式為：

γm為信號采集點到測量點距離，γc為空間位置到相對計算擾動測量點距離，△dm距離一致性標準距離，△Dm劑量差別標準距離。Gamma指數(shù)≤1表示計算通過，Gamma指數(shù)>1表示計算失敗。本文計算Gamma指數(shù)采用IBA的OmniPro I'mRT分析軟件，計算中△dm和△Dm選擇3 mm。

1.4 統(tǒng)計學分析

所有獲得的Gamma通過率進行單因素的配對t檢驗分析[8]，使用的統(tǒng)計學軟件為SPSS 13.0，P<0.05有統(tǒng)計學意義。

2 結(jié)果

2.1 head-head通過率分析結(jié)果

在OmniPro I'mRT軟件中設(shè)置3mm、3%、閾值5為分析限制條件，使用最大劑量點對Dataset1中CA_head和Dataset2中QA_head的平面劑量和進行Gamma分析[9]得到head-head通過率，見圖2,此20例鼻咽癌IMRT計劃的head-head通過率均達到96.44%±2.3%以上，符合臨床要求。

圖2 head-head通過率分析圖

2.2 head-chest通過率分析結(jié)果

在OmniPro I'mRT軟件Dataset2中導入QA_chest平面劑量，在相同的算法和計算間隔下，按照同樣的3mm、3%、閾值5為分析限制條件，使用最大劑量點對Dataset1中CA_head和Dataset2中QA_chest的平面劑量和進行Gamma分析得到head-chest通過率，見圖3。此20例head-chest通過率與head-head通過率比較得到通過率為：83.78%±7.1%(P<0.05),有統(tǒng)計學意義，見表1。

圖3 head-chest通過率分析

2.3 head-water通過率分析結(jié)果

在OmniPro I'mRT軟件Dataset2中導入QA_water平面劑量，在相同的算法和計算間隔下，按照同樣的3 mm、3%、閾值5為分析限制條件，使用最大劑量點對Dataset1中CA_head和Dataset2中QA_water的平面劑量和進行Gamma分析得到head-water通過率，見圖4，此20例head-water通過率與head-head通過率比較得到：83.26%±7%(P<0.05),有統(tǒng)計學意義，見表1。

3 討論

圖4 head-water通過率分析

表1 3組IMRT計劃每個射野Gamma通過率比較/%

隨著放療設(shè)備的不斷發(fā)展更新，各大放療設(shè)備生產(chǎn)廠家研究出了大量不同的可用于病人劑量驗證的二維測量工具，包括：Sun Nuclear的Mapcheck、PTW的2D-ARRAY和IBA的MatriXX。它們均同時兼有點劑量和二維劑量分布實時在線測量的能力，從而大大提高了測量精度和測量效率。但是，調(diào)強劑量分布的特點是劑量梯度大且變化復雜，Clark和Robar[10]認為，適合測量三維劑量分布的探測器需要較高的空間分辨率。Bjorn Poppe(2007)[11]提出了影響二維探測器的空間分辨率的兩個標準：探頭大小和取樣間隔，結(jié)合Nyquist采樣定理及混疊效應(yīng)[12]，從原理上分析了2D-ARRAY二維矩陣電離室的取樣間隔對測量空間分辨率的影響。傳統(tǒng)使用膠片法,其精度高(可分辨0.04 mm相對位置誤差[13])，膠片法可實施相對等劑量分布驗證，包括集合劑量學驗證[11-12]和單野劑量學驗證[13-15]。前者測量模體(通常圓柱狀)一個或多個選擇層面的相對復合劑量分布(提供了復合劑量分布的直接信息)；后者測量平的模體中與每一單獨射野中心軸垂直的一個層面的相對劑量分布。但測量結(jié)果受曝光和沖洗條件影響[16]、質(zhì)保環(huán)節(jié)多、操作復雜等因素，近年來已經(jīng)被二維矩陣所代替[17]。

對于MatriXX來說，采樣點以空間分布排列，取樣間隔為7.6 mm，即采樣頻率為0.13 mm-1,則Nyquist頻率fN=1/(7.6×2)即0.066 mm-1。根據(jù)Nyquist采樣定理[12]，若使計劃系統(tǒng)的信號不失真，則所采集信號的最高空間頻率應(yīng)小于MatriXX的Nyquist頻率fN，即對照射子野的面積有一定的局限性。加速器角度響應(yīng)、空間分辨率線性插值對于MatriXX都有一定的影響。理論上肯定無法達到膠片劑量學驗證的精度，從而對于劑量分布的驗證需要給予更多考慮與分析。實驗發(fā)現(xiàn)，如果靶區(qū)的劑量梯度較大的區(qū)域比較多，則調(diào)強劑量驗證的通過率相對較低，所以調(diào)強劑量驗證的通過率與靶區(qū)的劑量均勻性存在一定關(guān)系。

目前還沒有一種模體能真實反應(yīng)人體的密度，需要計劃系統(tǒng)根據(jù)轉(zhuǎn)換曲線，插值計算出組織的物理密度，這插值計算出來的物理密度與實際組織密度存在差異，因此對于個性化放療物理計劃的驗證還有很大距離。

[1] 殷蔚伯,谷銑之.腫瘤放射治療學〔M〕.第3版.北京:中國協(xié)和醫(yī)科大學出版社,2002:106-107.

[2] Bortfeld T，Boyer AL，Schlegel W，et al.Realization and verification of three 2 dimensional conformal radio-therapy with modulated fields〔J〕.Int J Radiat Oncol Biol Phys，1994，30(4):899-908.

[3] Tsai JS，Wazer DE，Ling MN，et al.Dosimetric verify-cation of the dynamic intensity modulated radiation therapy of 92 patients〔J〕.Int J Radiat Oncol Biol Phys，1998，40(1)：213-230.

[4] 胡逸民.腫瘤放射物理學〔M〕.北京:原子能出版社，1999.

[5] Dong L,Antolak J,Salehpour M,et al.Patient-specific point dose measurement for IMRT monitor unit verification〔J〕.Int J Radiat Oncol Biol Phys,2003,56(3):867-877.

[6] Klein EE，Hanley J，Bayouth J，et al.Task Grouup 142 report：quality assurance of medical accelerators〔J〕.Med Phys,2009,36(9):4197-4212.

[7] 易俊林,高 黎,徐國鎮(zhèn),等.鼻咽癌調(diào)強放射治療靶區(qū)勾畫-中國醫(yī)學科學院腫瘤醫(yī)院經(jīng)驗總結(jié)〔J〕.腫瘤防預(yù)與治療,2011,24(3):157-164.

[8] 金丕煥.醫(yī)用統(tǒng)計學〔M〕.上海:上海醫(yī)科大學出版社，1993.

[9] 張建英，江 濤.Pinnacle計劃系統(tǒng)調(diào)強放療驗證Gamma通過率的研究〔J〕.中國醫(yī)療器械雜志，2011,35(6):422-424.

[10] Robar JL,Clark BG.A practical technique for verification of three-dimensional conformal dose distributions in stereotactic radiosurgery〔J〕.Med Phys，2000,27(5)：978-987.

[11] Poppe B,Blechschmidt A,Djouguela A,et al.Two-dimensional ionization charmber arrays for IMRT Plan variation〔J〕.Med Phys,2006,33(4):1005-1015.

[12] Sastre-Padro M,van der Heide UA,Welleweerd H.An accurate calibration method of multileaf collimator valid for conformal and intensity modulated radiation treatments 〔J〕.Phys Med Biol,2004,49(12):2631-2640.

[13] Kapulsky A,Mullokandow E,Gejerman G.An automated phantom film QA procedure forintensity modulated radiation therapy〔J〕.Med Dosim,2002，27(3):201-207.

[14] Saw CB,Ayyangar KM,Zhen W,et al.Commissioning and quality assurance for MLC-based IMRT〔J〕.Med Dosim,2001,26(2):125-133.

[15] Xing L,Curran B,Hill R,et al.Dosimetric verrfication of a commercial inverse treatment planning system〔J〕.Phys Med Biol,1999,44(2):463-478.

[16] Wang X,Spirou S,LoSasso TH.et al.Dosimetric verification of intensity modulated fields〔J〕.Med Phys,1996,23(3):317-327.

[17] 鞠建中,王運來,馬林等.用二維電離室矩陣驗證多葉準直器葉片到位精度〔J〕.中華放射腫瘤學雜志,2006,15(4):335-338.

(編輯：吳小紅)

Impact of Different Tissue Density on Dose Verification Pass Rate of IMRT Plan

HUHaiqin,CHENXinde,ZHANGHuaiwen,etal.JiangxiCancerHospital,Nanchang,330029

Objective To study the impact of different tissue density on verification pass rate of IMRT plan.Methods CT electron density table of head (120KV,380mAs),lower chest (140KV,100mAs) and reference uniform phantom (120KV,130mAs) modes were established,and 20 cases of NPC IMRT plans were to use these 3 CT electron density table design 3 group of plans. Results In the OmniPro I’mRT software 3mm,3%,5 threshold limits were set for the analysis,these 20 cases of NPC IMRT plans head-head pass rate reached (96.44 ±2.3%);head-chest pass rate was compared with the rate of head-head (P<0.05);head-water pass rate was compared with the rate of head-head (P<0.05).Conclusion For IMRT optimization,selecting the correct electron density table is the guarantee to realize maximum dose pass rate.

Tissue density;IMRT optimization;Gamma pass rates

330029 江西省腫瘤醫(yī)院

江西省衛(wèi)生廳科技計劃項目(編號：20155480)

10.3969/j.issn.1001-5930.2016.09.019

R73-36

A

1001-5930(2016)09-1448-05

2016-02-10

2016-05-01)