2Uludağ Üniversitesi Fen-Edebiyat Fakültesi, Fizik Bölümü, Bursa
Summary
AMAÇÇalışmanın amacı, bilgisayarlı tedavi planlama sisteminde (BTPS) hesaplanan doz değerleri ile farklı volümlü iyon odaları ve linear 2D-Array kullanılarak ölçülen dozları karşılaştırarak dozimetrik kontrolü sağlamaktır.
GEREÇ VE YÖNTEM
Farmer, Semiflex ve PinPoint iyon odalarıyla ölçümleri alınan
alanlar BTPS'de oluşturuldu. 2D-Array ile açık alan, izosentrik,
oblik ve wedge filtreli alanlarda ölçümler alındı. Bu
alanlar BTPS'de oluşturularak nokta doz okumaları yapıldı ve
profilleri oluşturuldu. Ölçülen ve planlanan izodoz dağılımları
karşılaştırıldı.
BULGULAR
PinPoint iyon odasının en yüksek, Farmer iyon odasının ise
en düşük doz okuması yaptığı görüldü. BTPS ile en uyumlu
sonuçlar Semiflex iyon odasıyla elde edildi. 2D-Array kullanılarak
yapılan okumalarda merkezi eksende BTPS ile uyum
gözlenirken, build-up ve penumbra bölgelerinde %20'lere çıkan
fark gözlendi.
SONUÇ
İyon odası ölçümleri incelendiğinde, BTPS'ye en yakın sonuçlar
Semiflex iyon odası ile elde edilmiştir. 2D-Array kullanılarak
yapılan okumalarda gözlenen uyumsuzluğun nedenleri;
build-up bölgesindeki doz değişkenliği ve ışın modellemesinde
ölçülen profiller ile algoritma tarafından oluşturulan profillerin
penumbra bölgelerinde tam olarak çakışmaması olarak bulundu.
Introduction
Kanser günümüzün en önemli sağlık problemlerinden biridir. Tüm kanserli hastaların %60-80’i hastalık tanısı aldıktan sonra geçen süreçte en az bir kez radyoterapi görmektedirler. Radyoterapide en önemli kural, tümöre maksimum dozu verirken, çevresindeki riskli organların ve sağlıklı dokuların mümkün olan en az dozu almasını sağlamaktır. Bu amaca yönelik olarak üç boyutlu konformal radyoterapi (3D CRT), yoğunluk ayarlı radyoterapi (Intensity Modulated Radiation Therapy, IMRT), organ hareketlerini takip ederek yapılan görüntü takipli radyoterapi (Image Guided Radiation Therapy, IGRT) ve tomoterapi gibi gelişmiş teknikler uygulanmaktadır.3D CRT’de öncelikle, hastaya tedavi bölgesine göre pozisyon verme ve immobilizasyonun sağlanması gerekir. Tedavi, ilgili bölgenin bilgisayarlı tomografi (BT) kesitleri kullanılarak yapılan sanal simülasyon ile başlar. Görüntüleme yöntemine bağlı olarak, görünen tümör, kritik yapılar ve ilgili diğer organ sınırları tüm kesitlerde konturlenir.
Üç boyutlu BTPS yazılımı ile konturlenen kesitler üzerinde alanlar belirlenir. Bu sistemin en kullanışlı özelliklerinden biri, hedef volüm ve diğer yapıların görülebildiği ışın gözü görüntüsünün (Beam Eye View) elde edilmesidir. BEV opsiyonu ile hasta anatomisinin radyasyon kaynağının bulunduğu noktadan görüntüsü elde edilir.[1] Huzmenin nokta kaynaktan çıktıktan sonra ayrılması göz önüne alınarak, hedef volüm ve risk altındaki organlara, blokların ve kolimatörlerin doğru yerleştirilmesi sağlanır. Böylece, alan marj belirlemeleri daha sağlıklı yapılabilmekte, PTV’nin (Planning Target Volume) yüksek izodoz bölgesi (tanımlanan dozun %95’inden fazla) içinde kalması sağlanabilmektedir.
Enerji, wedge ve ışın ağırlıkları seçildikten sonra planın 3 boyutlu doz hesaplatılması yapılır. Doz dağılımı değerlendirilir ve istenen izodoz dağılımı incelenir. İzodoz dağılımı PTV’yi tam olarak sarıyorsa fraksiyon sayısı, verilmesi gereken tedavi dozu ve izodoz eğrisi seçilir.
DVH (Doz Volüm Histogramı) planın değerlendirilmesinde önemli rol oynar. Tedavi planlama sistemlerinde yapılan planların diferansiyel ve kümülatif doz-volüm histogramları elde edilir. Bu histogramlarda; bir plan için hedef tümör volümü ya da risk altındaki organların aldıkları dozlar değerlendirilebilir.[1] Böylece uygun tedavi alanları kullanılarak tümörü oluşturan hücreler ortadan kaldırılırken sağlıklı dokular minimum zarar görür.
Planlama algoritmaları
Modern bilgisayarlı planlama sistemleri ile sanal
ortamda gerçeğe en yakın doz hesaplama algoritmaları
oluşturulabilir. Bu algoritmalar; düzeltme/
ölçüm tabanlı algoritmalar, model tabanlı ya da
Monte Carlo tabanlı olabilir.
Düzeltme / ölçüm tabanlı algoritmalar
Bu algoritmada, referans koşullar altında düzenli
tedavi alanları için su fantomunda yüzde derin
doz, doz profilleri ve output faktörü ölçümleri
yapılmalıdır. Hasta dozu, belirli tedavi alanları
için doku eksikliği, doku imhomojenitesi gibi pek
çok düzeltme uygulanarak bulunur. Bu metot tamamen
ölçüm verilerine dayanır ve oldukça hızlıdır.
Su fantomunda ölçülen derin dozlar arasında
interpolasyon yaparak dozu hesaplar ve çeşitli
derinliklerde alınan doz profillerini kullanır. Doku
düzensizlikleri, inhomojenite ve ikincil elektronlar
doz hesaplamasında göz ardı edilir.[2]
Model tabanlı algoritmalar
Monte Carlo simülasyonuyla hesaplanmış olan
ışın karakteristikleri kullanılır. Ortamda oluşan etkileşimlerin
birim elementi olan kerneller, Monte
Carlo yöntemi ile hesaplanır. Doz kernel, farklı seviyelerdeki
enerji aktarımını ve primer foton doku
etkileşimleri ile oluşan sudaki doz depolanmasını
tanımlar. Model tabanlı algoritmalar, heterojen ortamda
absorbe dozun daha gerçekçi tanımlamasını
yapar. Hasta kesitlerindeki Hounsfield Unit ile inhomojen
hasta anatomisi daha iyi örneklenir.
Primer fotonların hesaplanan enerji akısı, hastada enerji absorbsiyon ve transportunun hesabı inputdata olarak kullanılır.[2]
FFT convolution ve superposition algoritmaları volümde depolanan dozu iki kısımda hesaplar.
İlk olarak primer kernelleri kullanarak primer elektron dozunu ve daha sonra saçılma kernellerini kullanarak saçılan foton dozunu hesaplar.
FFT convolution, kernelleri kartezyen koordinatlarda gösterir ve her yerde değişmez. Superposition ise kernelleri küresel koordinatlarda gösterir ve lokal elektron yoğunluğu varyasyonlarıyla değişime izin verir. Böylece inhomojen yapılarda daha iyi sonuç elde edilir.
Monte Carlo tabanlı algoritmalar
Madde içinde çok sayıda foton ve parçacığın etkileşimini
simüle eder. Foton ve elektronların olası
etkileşimlerini belirlemek için temel fizik yasaları
kullanılır. En doğru doz hesaplama algoritması
olarak kabul edilir.[2]
Monte Carlo algoritmasında, simülasyona E enerjili bir foton veya elektronla başlanır.
Ortamda hareket edecek parçacığın hareket mesafesi belirlenir. Yani etkileşim koordinatları girilir. Oluşabilecek etkileşimlerin türleri yazılır. Oluşacak parçacığın enerjisi ve yönü belirlenir. Tüm bu adımlar parçacığın enerjisi alt limitin altına inene kadar veya parçacık ilgilenilen bölgenin dışına çıkana kadar tekrarlanır.
Tedavi planlama sisteminin kalite kontrolü
Radyoterapide klinik sonuçları en iyi seviyeye
getirmek ve tedavi alan her hasta için yüksek
ve gerekli kaliteyi sağlamak için kayda değer çaba
harcanmaktadır. Radyoterapi süreci oldukça karmaşıktır.
İyi radyoterapi sonuçları ve tedavi güvenliği,
planlanan radyasyonun belirlenen hedef
alana doğru dozda verilmesi ile elde edilir. Radyoterapide
kalite kontrolün amacı, planlama ve tedavinin
her aşamasında oluşabilecek hataları önlemektir.[3] Bu nedenle tedavi aygıtları ve simülatörlerin,
tedavi planlamanın, hastaya ait bilgilerin
belli aralıklarla ölçülmesi ve kontrol edilmesi gerekmektedir.
Hasta tedavisinde kalitenin sağlanması için BTPS’de hesaplanan doz dağılımlarının dozimetrik araçlar kullanılarak doğruluğunun kontrol edilmesi gerekmektedir. Bu araçlar, iyon odası, film dozimetre, termolüminesans dozimetre (TLD) veya 2D-Array olabilir.
Bu çalışmanın amacı, bilgisayarlı tedavi planlama sisteminde (BTPS) hesaplanan doz değerleri ile farklı volümlü iyon odaları ve linear 2D-Array kullanılarak ölçülen nokta dozları karşılaştırarak tedavi planlama sisteminin dozimetrik kontrolünü sağlamaktır.
Methods
Bu çalışmada Siemens MD2 lineer hızlandırıcıda 6 MV foton enerjisi kullanıldı. Ölçümler yapılmadan önce output değerleri uygun sıcaklık ve basınç koşullarında kontrol edildi. Siemens SOMATOM Emotion Duo BT ve simülatör ünitesi kullanılarak RW3 katı su fantomlarının ve 2D-Array’in tomografisi çekildi. Dış konturlaması yapılıp bilgisayarlı planlama sistemine aktarıldı.Katı su fantomu ölçümleri, RW3 katı su fantomu kullanılarak merkezi eksende 6 MV foton enerjisinde d=10 cm’de; 3x3, 4x4, 5x5, 6x6,7x7, 8x8, 9x9, 10x10, 12x12 ve 15x15 cm2 alanda SSD 100 cm’de PTW Semiflex, PTW Farmer, PTW PinPoint iyon odaları ve PTW Unidos elektrometre kullanılarak yapıldı. Her bir ışınlama için 100 MU uygulandı. Ölçümlerde her iyon odası için özel dizayn edilmiş katı su fantom plakası kullanıldı. Her bir ölçüm 3 kez tekrarlanıp ortalaması alındı.
2D-Array ile yapılan ölçümlerde, 6 MV foton enerjisi için, 5x5 ve 10x10 alanlarda SSD=100 cm’de alan merkezinde ve merkezden 1, 5, 7 cm uzaklıkta, d=1, d=3, d=5, d=8 ve d=10 cm derinliklerinde nokta doz değerleri ölçüldü. Ayrıca SSD= 85 cm’de, d=15 cm’de nokta doz değerleri ölçüldü. Yine bu enerjiler için 5x5 ve 10x10 alanlarda SSD=100 cm’de rutinde sık kullanılan 15, 30 ve 45 wedge filtreler kullanılarak SSD=100 cm’de d=1, d=3 ve d=5 cm derinliklerinde nokta doz ölçümü yapıldı.
Son olarak, G=45 için her iki foton enerjisinde 5x5 ve 10x10 alanlarda d=1, d=3 ve d=5 cm için nokta doz ölçüldü. Her bir ışınlama için 100 MU uygulandı.
Tüm alanların planları CMS XIO tedavi planlama sisteminde convolution algoritması kullanılarak hesaplatıldı. Dozimetrik ekipmanlar ile yapılan ölçümler sonrası elde edilen veriler ile planlama sisteminde hesaplatılan değerler PTW Verisoft 3.1 ve PTW Multicheck dozimetri yazılımı kullanılarak karşılaştırıldı.
Results
6 MV foton enerjisi kullanılarak yapılan iyon odası ölçümleri, BTPS ile karşılaştırıldığında, tüm iyon odaları için hata oranının %5’in içinde olduğu Tablo 1’de gösterilmiştir.Hassas volümlü PinPoint iyon odasının en yüksek doz okuması yaparken, Farmer iyon odasının en az doz okuması yaptığı görülmektedir. BTPS ile en uyumlu sonuçlar Semiflex iyon odası ile elde edilmiş, aradaki fark en küçük alan olan 3x3 cm2 hariç, %1.5’den az olduğu görülmüştür.
PTW Verisoft yazılımı ile yapılan plan karşılaştırmalarında elde edilen izodoz dağılımları ve alan doz profilleri incelenmiştir. 2D arrayle alınan ölçümler BTPS’den elde edilen koronal kesitlerle Verisoft yazılım programı kullanılarak karşılaştırıldı. Karşılaştırmada Gamma Index metodu kullanıldı. Bu metot karşılaştırmaları, 3 mm uyuşma mesafesi (DTA; distance to agreement), %3 doz farkı (DD; dose-difference) tolerans limiti alarak yapmaktadır. Hızlı doz düşüşü olan bölgelerde, özellikle penumbra bölgesinde karşılaştırma yapılırken bu limitler yükseltilebilmektedir. Bu çalışmada tolerans limitleri 3 mm’de %3 olarak alınmıştır.
Buna göre 3 mm uyuşma mesafesinde, %3 doz farkı kriteri için gamma değerlendirmesini geçen noktalar bulundu. Doz dağılımlarında planlama ile ölçüm değerleri arasında uyum gözlenirken, penumbra bölgesindeki uyumsuzluklar dikkat çekmiştir.
Tablo 2’de 2D-Array ile yapılan açık alan okumalar gösterilmiştir. 2D-Array ile yapılan açık alan ölçümlerde bazı farklılıklar görülmektedir. Buildup ve penumbra bölgelerindeki farklılılar %20’lere kadar çıkmaktadır. 6 MV enerjisi için tüm alanlarda build-up bölgesi dışında kalan bölgelerdeki merkezi eksen okumaları %2’nin altındadır.
Tablo 2: 2D-Array ile yapılan açık alan doz ölçümleri
Şekil 1’de 6 MV foton enerjisi kullanılarak 10x10 cm2 alanda G=45º’de SSD=100 cm’de 5 cm derinliğinde PTW Verisoft yazılımı ile karşılaştırılmış izodoz profilleri görülmektedir. Penumbra bölgesinde uyumsuzluk görülürken, izodoz dağılımı incelendiğinde %23’e varan fark gözlenmiştir.
Sekil 1: 6 MV 10x10 cm2 G=45º SSD=100 d=5 cm izodoz dağılımı.
Şekil 2’de 6 MV foton enerjisi kullanılarak 10x10 cm2 alanda 15º wedge filtre kullanılarak SSD=100 cm’de 5 cm derinliğinde PTW Verisoft yazılımı ile elde edilen izodoz dağılımı incelendiğinde, penumbra bölgesinde uyumsuzluk görülürken, %10’a varan fark gözlenmiştir.
Sekil 2: 6 MV 10x10 cm2 W=15º SSD=100 d=5 cm izodoz dağılımı.
6 MV foton enerjisi kullanılarak 10x10 cm2 alanda 30º wedge filtre kullanılarak SSD=100 cm’de 5 cm derinliğinde PTW Verisoft yazılımı ile elde edilen izodoz dağılımı dağılımı, penumbra bölgesinde uyumsuzluk görülürken, %12’ye varan farklılık gözlenmiştir. Aynı şekilde 45º wedge filtre kullanılarak SSD=100 cm’de 5 cm derinliğinde PTW Verisoft yazılımı ile elde edilen izodoz dağılımı incelendiğinde, penumbra bölgesinde uyumsuzluk görülürken %16’ya varan farklılık gözlenmiştir.
Discussion
Radyoterapi tedavi güvenliği, planlanan radyasyonun belirlenen hedef alana doğru dozda verilmesi ile elde edilir. Bu amaçla radyoterapi bilgisayarlı planlama sistemlerinin doz hesaplama algoritma sürecinin dozimetrik kontrolünü yapılması zorunludur. Çalışmamızda, Farmer, Semiflex ve PinPoint iyon odaları ile 2D-Array kullanılarak bilgisayarlı tedavi planlama sisteminin dozimetrik kontrolü yapıldı. Doz verifikasyon bölgeleri olarak ise ışınlama alanı, penumbra bölgesi, ışınlama alanı dışı, build-up bölgesi ve merkezi eksen belirlendi.6 MV foton enerjisi kullanılarak yapılan iyon odası ölçümleri, BTPS sonuçları ile karşılaştırıldığında, tüm iyon odaları için hata oranının %5’in içinde olduğu görüldü. En yüksek doz okuması PinPoint iyon odası ile elde edildi. BTPS ile en uyumlu sonuçlar Semiflex iyon odası ile elde edildi, en küçük alan olan 3x3 cm2 hariç, farkın %1.5’den az olduğu görüldü.Konvansiyonel radyoterapi için kabul edilebilir hata payı ICRU[4] tarafından ±%5 olarak belirlenmiştir.
Leybovich ve ark.[5] absolute doz verifikasyonunda farklı hacme sahip 3 iyon odasıyla yaptıkları çalışmada büyük volümlü iyon odasının absolute doz açısından hesaplanan değerlere daha yakın sonuçlar verdiğini bulmuşlardır.
Gonzalez-Castano ve ark.[6] Farmer, PinPoint ve Semiflex iyon odalarını kullanarak yaptıkları çalışmada, Farmer iyon odası ile elde ettikleri sonuçların BTPS’ye daha yakın olduğu sonucuna ulaşmışlardır.
Laub ve ark.[7] ise 0.6 cm3 Farmer iyon odası ve 0.015 cm3 PinPoint iyon odası ile yaptıkları ölçümleri hesaplanan değerler ile karşılaştırdıklarında 0.015 cm3 volümlü iyon odasının hesaplanan değerlere daha yakın sonuçlar verdiğini görmüşlerdir.
Yaptığımız çalışmada, BTPS’de hesaplanan dozlar ile 2D-Array ölçülen dozlar arasındaki farklar build-up ve penumbra bölgeleri dışında, %4’ün içindeydi. Ancak penumbra bölgesi ve alan dışı noktalarda görülen farklılık, beklenildiği gibi %4’ün çok üzerindeydi. Build-up bölgesindeki doz ile planlama sisteminin hesapladığı doz değerleri arasındaki farklılık %25’lere kadar çıktığı gözlemlendi. Bu sonuçlar literatürle uyumlu bulunmuştur. [8-10]
Van Dyk,[8] ışınlama alanı içinde ve alan dışı bölgelerde %3, build-up bölgesinde 4 mm ve merkez eksende %4 tolerans olabileceğini bildirmiştir. Task Group 53,[9] ışınlama alanı içindeki toleransları açık alan için %1-2, MLC, blok, asimetrik alan için %3, farklı SSD’deki alan için %1.5 ve inhomojen bölge için %5 olarak belirlemiştir. Alan dışı bölgelerde ise, açık alan için %2-5, blok için %3, wedge, MLC ve bolus için %5 ve inhomojen bölge için %5 tolerans belirtirken, build-up bölgesinde %20-50 ve merkez eksende %1-3 toleransı kabul etmiştir.
SGSMP (Swiss Society of Radiobiology and Medical Physics),[10] ışınlama alanı içinde %2, alan dışı bölgelerde %3, build-up bölgesinde 3 mm ve merkezi eksende %4 tolerans olabileceğini bildirmiştir.
Wiezorek ve ark.[11] çalışmalarında 3 adet 2D dedektör sisteminde (Mapcheck diyot sistemi, IMRT QA sistilasyon detektörü ve seven 29 iyon chamber) elde ettikleri sonucu Konrad BTPS’de hesaplanan değerlerle ve standart dozimetrik aletlerle ölçülen değerlerle karşılaştırmışlardır. Bu 2D detektörlerle elde ettikleri dozların BTPS’de hesaplanan dozlara göre deviasyonu, değişken doz bölgelerinde bile %5’i geçmemektedir
2D-Array’in duyarlılık, linearite ve tekrar kullanılabilirliği gibi özellikleri Spezi ve ark.[12] ile Poppe ve ark.[13] tarafından araştırılmıştır, IMRT planlarının doğrulaması için oldukça verimli ve kullanışlı olduğu sonucuna varmışlardır.
Saminathan ve ark.[14] 6 ve 18 MV’de 2D-Array kullanarak 10x10 cm2 alanda elde ettikleri sonuçları tedavi planlama sistemiyle karşılaştırmışlar, %95’e varan uyum gözlemlemişlerdir.
Build-up bölgesindeki farklılıkların sebepleri; bu bölgede dozun değişken olması, bu bölgede hesaplama algoritmalarının yetersiz kalması, farklılık planlama sistemlerinde ışın modellemesi sırasında ölçülen profiller ile algoritma tarafından oluşturulan profillerin penumbra bölgelerinde tam olarak çakışmamasından ve doğru ölçüm yapmanın zorluğu olarak bulundu.
Çalışmamızdan elde ettiğimiz sonuçlar doğrultusunda, hem iyon odalarının hem de 2D-Array’in ölçümlerinde BTPS convulution algoritması ile elde edilen doz hesaplarının gerçeğe yakın ve literatürle uyumlu olduğu bulunmuştur. Ayrıca 2D-Array’in BTPS kalite kontrolünde oldukça başarılı olduğu, özellikle merkez eksende minimum hata payı ile sonuca ulaştığımız görülmektedir. 2D-Array, pratik kurulumu ve hızlı sonuç elde etmesi ile tedavi plan doğrulamaları için hem zaman tasarrufu sağlayacak, hem de plan doğruluğunu maksimum düzeyde sağlayarak kaliteli tedaviye ulaşılmasına yardımcı olacaktır.
References
1) Khan FM. The Physics of radiation therapy. 3rd ed.
Minnesota: Williams & Wilkins; 2003.
2) Evans MDC. Computerized treatment planning systems
for external photon beam radiotherapy. International
Atomic Energy Agency Publication. Montreal,
Quebec, Canada: 2006. p. 387-406.
3) Stuk, Radiation and Nuclear Safety Authority, Quality
Assurance in Radiotherapy, Guide ST 2.1, 2003
4) International comission on radiation units and measurements.
Report 62, Prescribing, recording and reporting
photon beam therapy (Supplement to ICRU Report 50),
1999.
5) Leybovich LB, Sethi A, Dogan N. Comparison of ionization
chambers of various volumes for IMRT absolute
dose verification. Med Phys 2003;30(2):119-23.
6) González-Castaño D, Pena J, Sánchez-Doblado F, Hartmann GH, Gómez F, Leal A. The change of response of
ionization chambers in the penumbra and transmission
regions: impact for IMRT verification. Med Biol Eng
Comput 2008;46(4):373-80.
7) Laub WU, Wong T. The volume effect of detectors in
the dosimetry of small fields used in IMRT. Med Phys
2003;30(3):341-7.
8) Van Dyk J. Quality assurance of radiation therapy planning
systems: current status and remaining challenges.
Int J Radiat Oncol Biol Phys 2008;71:23-7.
9) Fraass B, Doppke K, Hunt M, Kutcher G, Starkschall
G, Stern R, Van Dyke J. American Association of Physicists
in Medicine Radiation Therapy Committee Task
Group 53: quality assurance for clinical radiotherapy
treatment planning. Med Phys 1998;25(10):1773-829.
10) SGSMP (Swiss Society of Radiobiology and Medical
Physics): Vanselaar J, Welleweerd H, Munheer B:
Tolerancesfor theaccuracyofphoton beam dose calculationsof
treatment planning systems, Elsevier Radiotherapy
and Oncology, 2009.
11) Wiezorek T, Banz N, Schwedas M, Scheithauer M,
Salz H, Georg D, et al. Dosimetric quality assurance
for intensity-modulated radiotherapy feasibility
study for a filmless approach. Strahlenther Onkol
2005;181(7):468-74.
12) Spezi E, Angelini AL, Romani F, Ferri A. Characterization
of a 2D ion chamber array for the verification
of radiotherapy treatments. Phys Med Biol
2005;50(14):3361-73.