Inclusion of organ deformation in dose calculations.

A previously described system for modeling organ deformation using finite element analysis has been extended to permit dose calculation. Using this tool, the calculated dose to the liver during radiotherapy can be compared using a traditional static model (STATIC), a model including rigid body motion (RB), and finally a model that incorporates rigid body motion and deformation (RBD). A model of the liver, consisting of approximately 6000 tetrahedral finite elements distributed throughout the contoured volume, is created from the CT data obtained at exhale. A deformation map is then created to relate the liver in the exhale CT data to the liver in the inhale CT data. Six intermediate phase positions of each element are then calculated from their trajectories. The coordinates of the centroid of each element at each phase are used to determine the dose received. These intermediate dose values are then time weighted according to a population-modeled breathing pattern to determine the total dose to each element during treatment. This method has been tested on four patient datasets. The change in prescribed dose for each patient's actual tumor as well as a simulated tumor of the same size, located in the superior, intermediate, and inferior regions of the liver, was determined using a normal tissue complication model, maintaining a predicted probability of complications of 15%. The average change in prescribed dose from RBD to STATIC for simulated tumors in the superior, intermediate, and inferior regions are 4.0 (range 2.1 to 5.3), -3.6 (range -5.0 to -2.2), and -14.5 (range -27.0 to -10.0) Gy, respectively. The average change in prescribed dose for the patient's actual tumor was -0.4 Gy (range -4.1 to 1.7 Gy). The average change in prescribed dose from RBD to RB for simulated tumors in the superior, intermediate, and inferior regions are -0.04 (range -2.4 to 2.2), 0.2 (range -1.5 to 1.9), and 3.9 (range 0.8 to 7.3) Gy, respectively. The average change in the prescribed dose for the patient's actual tumor was 0.7 Gy (range 0.2 to 1.1 Gy). This patient sampling indicates the potential importance of including deformation in dose calculations.

[1]  J. Wong,et al.  The use of active breathing control (ABC) to reduce margin for breathing motion. , 1999, International journal of radiation oncology, biology, physics.

[2]  R L Ehman,et al.  Tissue characterization using magnetic resonance elastography: preliminary results. , 2000, Physics in medicine and biology.

[3]  Daniel Normolle,et al.  Analysis of radiation-induced liver disease using the Lyman NTCP model. , 2002, International journal of radiation oncology, biology, physics.

[4]  K. Brock,et al.  Determination of ventilatory liver movement via radiographic evaluation of diaphragm position. , 2001, International journal of radiation oncology, biology, physics.

[5]  C. Ling,et al.  Evaluation of respiratory movement during gated radiotherapy using film and electronic portal imaging. , 2002, International journal of radiation oncology, biology, physics.

[6]  C. Burman,et al.  Calculation of complication probability factors for non-uniform normal tissue irradiation: the effective volume method. , 1989, International journal of radiation oncology, biology, physics.

[7]  Shinichi Shimizu,et al.  Registration accuracy and possible migration of internal fiducial gold marker implanted in prostate and liver treated with real-time tumor-tracking radiation therapy (RTRT). , 2002, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[8]  E. Larsen,et al.  A method for incorporating organ motion due to breathing into 3D dose calculations. , 1999, Medical physics.

[9]  J. Lyman.  Complication Probability as Assessed from Dose-Volume Histograms , 1985 .

[10]  J M Balter,et al.  Technical note: creating a four-dimensional model of the liver using finite element analysis. , 2002, Medical physics.

[11]  H. Kubo,et al.  Respiration gated radiotherapy treatment: a technical study. , 1996, Physics in medicine and biology.

[12]  R Mohan,et al.  A method of incorporating organ motion uncertainties into three-dimensional conformal treatment plans. , 1996, International journal of radiation oncology, biology, physics.

[13]  J M Balter,et al.  Potential benefits of eliminating planning target volume expansions for patient breathing in the treatment of liver tumors. , 1997, International journal of radiation oncology, biology, physics.

[14]  Charles R. Meyer,et al.  Demonstration of accuracy and clinical versatility of mutual information for automatic multimodality image fusion using affine and thin-plate spline warped geometric deformations , 1997, Medical Image Anal..

[15]  D. Yan,et al.  Organ/patient geometric variation in external beam radiotherapy and its effects. , 2001, Medical physics.

[16]  R K Ten Haken,et al.  Use of Veff and iso-NTCP in the implementation of dose escalation protocols. , 1993, International journal of radiation oncology, biology, physics.

[17]  H. Shirato,et al.  Four-dimensional treatment planning and fluoroscopic real-time tumor tracking radiotherapy for moving tumor. , 2000, International journal of radiation oncology, biology, physics.

[18]  K. Lam,et al.  Improvement of CT-based treatment-planning models of abdominal targets using static exhale imaging. , 1998, International journal of radiation oncology, biology, physics.

[19]  R K Ten Haken,et al.  The reproducibility of organ position using active breathing control (ABC) during liver radiotherapy. , 2001, International journal of radiation oncology, biology, physics.

[20]  K Miller,et al.  Constitutive modelling of abdominal organs. , 2000, Journal of biomechanics.

[21]  D. Gladstone,et al.  A numerical simulation of organ motion and daily setup uncertainties: implications for radiation therapy. , 1997, International journal of radiation oncology, biology, physics.