Prog Med Phys.  2014 Jun;25(2):100-109. 10.14316/pmp.2014.25.2.100.

Development of New 4D Phantom Model in Respiratory Gated Volumetric Modulated Arc Therapy for Lung SBRT

• Yoon K
• Kwak J
• Cho B
• Song S
• Lee S
• Ahn S
• Nam S

• KMID: 1708177
• DOI: http://doi.org/10.14316/pmp.2014.25.2.100

Abstract

In stereotactic body radiotherapy (SBRT), the accurate location of treatment sites should be guaranteed from the respiratory motions of patients. Lots of studies on this topic have been conducted. In this letter, a new verification method simulating the real respiratory motion of heterogenous treatment regions was proposed to investigate the accuracy of lung SBRT for Volumetric Modulated Arc Therapy. Based on the CT images of lung cancer patients, lung phantoms were fabricated to equip in QUASAR(TM) respiratory moving phantom using 3D printer. The phantom was bisected in order to measure 2D dose distributions by the insertion of EBT3 film. To ensure the dose calculation accuracy in heterogeneous condition, The homogeneous plastic phantom were also utilized. Two dose algorithms; Analytical Anisotropic Algorithm (AAA) and AcurosXB (AXB) were applied in plan dose calculation processes. In order to evaluate the accuracy of treatments under respiratory motion, we analyzed the gamma index between the plan dose and film dose measured under various moving conditions; static and moving target with or without gating. The CT number of GTV region was 78 HU for real patient and 92 HU for the homemade lung phantom. The gamma pass rates with 3%/3 mm criteria between the plan dose calculated by AAA algorithm and the film doses measured in heterogeneous lung phantom under gated and no gated beam delivery with respiratory motion were 88% and 78%. In static case, 95% of gamma pass rate was presented. In the all cases of homogeneous phantom, the gamma pass rates were more than 99%. Applied AcurosXB algorithm, for heterogeneous phantom, more than 98% and for homogeneous phantom, more than 99% of gamma pass rates were achieved. Since the respiratory amplitude was relatively small and the breath pattern had the longer exhale phase than inhale, the gamma pass rates in 3%/3 mm criteria didn't make any significant difference for various motion conditions. In this study, the new phantom model of 4D dose distribution verification using patient-specific lung phantoms moving in real breathing patterns was successfully implemented. It was also evaluated that the model provides the capability to verify dose distributions delivered in the more realistic condition and also the accuracy of dose calculation.

Keyword

4D phantom; 3D print; Lung SBRT; Respiratory Gated VMAT

MeSH Terms

Humans
Lung Neoplasms
Lung*
Plastics
Radiosurgery
Radiotherapy, Intensity-Modulated*
Respiration
Plastics

Figure

• 

  Fig. 1. CT image based 3D Phantom design using the Eclipse RTP system. (a) Frontal image for a lung cancer patient. (b) Rendering image of a lung phantom contour. (c) Dividing into the two equal parts of the phantom.

• 

  Fig. 2. (a) 3D printed lung phantom based new 4D lung phantom with quasar respiratory motion phantom. The insert was selected (b) 3D printed lung phantom or (c) homogeneous phantom.

• 

  Fig. 3. (a) System setup of the 4D lung phantom to simulate the respiratory motion of target in gated VMAT treatment. The fluoroscopy images of the lung phantom captured in (b) exhale and (c) inhale phases.

• 

  Fig. 4. (a) The picture of the homemade 3D printed Lung phantom. (b) 4D lung phantom consists of the Quasar respiratory motion phantom and 3D printed lung phantom. The EBT3 film located inside the lung phantom.

• 

  Fig. 5. 4D CT images of the lung phantom reconstructed using phase based sorting. From left to right in row: end inhalation (0%), middle exhalation (30%), end exhalation (60%), and middle inhalation (70%). Full motion was 7.2 mm and Gated motion (30%∼70% phase) was 2.6 mm.

• 

  Fig. 6. (a) Image fusion between patient CT and lung phantom CT. Red Arrows indicate the GTV contour. (b) The CT number of patient and phantom were respectively 78 HU and 92 HU in tumor regions.

• 

  Fig. 7. (a) The image of the lung phantom was shown the coronal. sagittal and axial. The lung tissue in the phantom was implemented using the form of a lattice structure. The mean value in red area was −784 HU. (b) The material of lung tissue in the phantom was replaced to lung tissue (Eclipse default, −725 HU) and air (−1000 HU). then gamma passing rates (3%/3 mm) were compared between the lung phantom with lattice structure and the phantom with replacement materials. Treatment planning dose applied AXB dose calculation algorithm.

• 

  Fig. 8. The phantom study were performed using the homogeneous phantom and 3D printed lung phantom. Gamma index map for 3%/3 mm criteria in the coronal plane, it was compared between the plan dose and film measurement under the three experimental conditions of (a) Static, (b) Gating with motion and (c) Nogating with motion. Treatment planning dose applied AAA and AXB dose calculation algorithm. Gating window included 30%∼70% phase. Arrows indicate the movement direction.

Reference

References

1. McNair HA, et al. Active Breathing Contr ol (ABC) in R adical Radiotherapy of Non-small Cell Lung Canc er (NS CLC). Clin Oncol (R Coll Radiol) 19: S 39 (. 2007.

2. Ko YE, et al. Effectiveness of breath hold with a ABC for SRS of lung ca nc er. J Lung Ca nc er. 4:42–47. 2005.

3. D'Souza WD, et al. The use of gated and 4D CT imaging in planning for ste reota ctic body radiation the rapy. Me d Dosim. 32:92–101. 2007.

4. Wurm RE, et al. Image guided respiratory gated hypofractionated stereotactic body radiation therapy (H-SBRT) for liver and lung tumors: initial exper ience. Ac ta Onc ol. 45:881–889. 2006.

5. Otto K, et al. Volumetric modula ted ar c the rapy: IM RT in a single gantry arc. Med P hys. 35(1):310–317. 2008.

6. Ulmer W, et al. A 3D photon superposition/convolution algorithm and its foundation on r esul t s of Monte Ca rlo calculations. Phys Med Biol. 50(8):1767–1790. 2005.

7. Gifford KA, et al. Comparison of a finite-element multigroup discrete-ordinates code with Monte Carlo for radiothera py ca l-culations. Phys Me d Biol. 51(9):2253. 2006.

8. Vassiliev ON, et al. Validation of a new grid-based Boltzmann equation solver for dose calculation in radiothera py with photon b e a ms. Phys Me d Biol. 55(3):581–598. 2010.

9. Keall P, et al. The management of respiratory motion in radiation oncology rep ort of AAPM Task Group 76. Med P hys. 33:3874–3900. 2006.

10. Li XA, et al. Point/counterpoint. Respiratory gating for r ad i ation the rapy is not rea dy for prime time. Med Phys. 34:867–870. 2007.

11. Andriy Fedorov, et al. 3D Slicer as an image computing platform for the Quantitative Imaging Network. Magnetic Resonance Ima ging. 30:1323–1341. 2012.

12. A. Sam Beddar, et al: Correlation be twee n internal fiducial tumor motion and external ma rke r motion for liver tumors imaged with 4D-CT. International Journal of Ra diation Onc ology and Biology a nd P hysi c s. 67:630–638. 2008.

13. Thomas SJ, et al. Relative electron density calibration of CT scanners for radiotherapy treatment planning. Br J Radiol. 72:781–786. 1999.

14. Kilby W, Sage J, Rabett V, et al. Tolerance levels for quality assura nce of electron de nsity values generated from CT in radiotherapy treatme nt pla nn i n g. Phys Med Biol. 47:1485–1492. 2002.