Faculty Bibliography

Magnetic resonance (MR) imaging of Gafchromic film causes perturbation to absolute dosimetry measurements; the purpose of this work was to characterize the perturbation and develop a correction method for it. Three sets of Gafchromic EBT2 film were compared: radiation (control), radiation followed by MR imaging (RAD + B), and MR imaging followed by radiation (B + RAD). The T1-weighted and T2-weighted MR imaging was performed using a 1.5T scanner with the films wedged between two chicken legs. Doses from 0 to 800 cGy were delivered with a 6MV linac. The time interval between radiation and MR imaging was less than 10 min. Film calibration was generated from the red channel. Microscopic imaging was performed on two pieces of film. The effect of specific absorption rate (SAR) was determined by exposing another three sets of films to low, medium, and high levels of SAR through a series of pulse sequences. No discernible preferential alignment was detected on the microscopic images of the irradiated film exposed to MRI. No imaging artifacts were introduced by Gafchromic film on any MR images. On average, 4% dose difference was observed between B + RAD or RAD + B and the control, using the same calibration curve. The pixel values between the B + RAD or RAD + B and the control films were found to follow a linear relationship pixel(Control) = 1.02 × pixel(B + RAD or RAD + B). By applying this correction, the average dose error was reduced to approximately 2%. The SAR experiment revealed a dose overestimation with increasing SAR even when the correction was applied. It was concluded that MR imaging introduces perturbation on Gafchromic film dose measurements by 4% on average, compared to calibrating the film without the presence of MRI. This perturbation can be corrected by applying a linear correction to the pixel values. Additionally, Gafchromic film did not introduce any imaging artifacts in any of the MR images acquired.

In electron and proton radiotherapy, applications of patient-specific electron bolus or proton compensators during radiation treatments are often necessary to accommodate patient body surface irregularities, tissue inhomogeneity, and variations in PTV depths to achieve desired dose distributions. Emerging 3D printing technologies provide alternative fabrication methods for these bolus and compensators. This study investigated the potential of utilizing 3D printing technologies for the fabrication of the electron bolus and proton compensators. Two printing technologies, fused deposition modeling (FDM) and selective laser sintering (SLS), and two printing materials, PLA and polyamide, were investigated. Samples were printed and characterized with CT scan and under electron and proton beams. In addition, a software package was developed to convert electron bolus and proton compensator designs to printable Standard Tessellation Language file format. A phantom scalp electron bolus was printed with FDM technology with PLA material. The HU of the printed electron bolus was 106.5 ± 15.2. A prostate patient proton compensator was printed with SLS technology and polyamide material with -70.1 ± 8.1 HU. The profiles of the electron bolus and proton compensator were compared with the original designs. The average over all the CT slices of the largest Euclidean distance between the design and the fabricated bolus on each CT slice was found to be 0.84 ± 0.45 mm and for the compensator to be 0.40 ± 0.42 mm. It is recommended that the properties of specific 3D printed objects are understood before being applied to radiotherapy treatments.

Tagged Meganetic Resonance Image (TMRI) provides a direct and noninvasive way to visualize the in-wall deformation of the myocardium. Due to the through-plane motion, the tracking of 3D trajectories of the material points and the computation of 3D strain field call for the nessicity of building 3D cardiac deformable models. The intersections of three stacks of orthogonal tagging planes are material points in the myocardium. With these intersections as control points, 3D motion can be reconstructed with a novel meshless deformable model. Volumetric meshless deformable models describe an object as point cloud inside the object boundary and the coordinate of each point can be written in parametric functions. A generic heart mesh is registered on the tagged MRI with polar decomposition. A 3D meshless deformable model is generated and deformed with MR image tagging lines. Volumetric meshless deformable models are deformed by calculating the dynamics function and minimizing the local Laplacian coordinates. The similarity transformation of each point is computed by assuming its neighboring points are making the same transformation. The deformation is computed iteratively until the control points match the target positions in the consecutive image frame. The 3D strain field is computed from the 3D displacement field with Moving Least Squares (MLS). We demonstrate that meshless deformable models outperformed the finite element method (FEM) and the spline method with a numerical phantom. Meshless deformable models can track the trajectory of any material point in the myocardium and compute the 3D strain field of any particular area. The experimental results on in vivo healthy and patient heart MRI show that the meshless deformable model can fully recover the myocardium motion in three dimension.

PURPOSE/OBJECTIVES/OBJECTIVE:Lung tumor motion may be impacted by heartbeat in addition to respiration. This study seeks to quantitatively analyze heart-motion-induced tumor motion and to evaluate its impact on lung cancer radiotherapy. METHODS/MATERIALS/METHODS:Fluoroscopy images were acquired for 30 lung cancer patients. Tumor, diaphragm, and heart were delineated on selected fluoroscopy frames, and their motion was tracked and converted into temporal signals based on deformable registration propagation. The clinical relevance of heart impact was evaluated using the dose volumetric histogram of the redefined target volumes. RESULTS:Correlation was found between tumor and cardiac motion for 23 patients. The heart-induced motion amplitude ranged from 0.2 to 2.6 mm. The ratio between heart-induced tumor motion and the tumor motion was inversely proportional to the amplitude of overall tumor motion. When the heart motion impact was integrated, there was an average 9% increase in internal target volumes for 17 patients. Dose coverage decrease was observed on redefined planning target volume in simulated SBRT plans. CONCLUSIONS:The tumor motion of thoracic cancer patients is influenced by both heart and respiratory motion. The cardiac impact is relatively more significant for tumor with less motion, which may lead to clinically significant uncertainty in radiotherapy for some patients.

The focus of this study is the angular dependence of two types of Metal Oxide Semiconductor Field Effect Transistor (MOSFET) dosimeters (MOSFET20 and OneDose/OneDosePlus) when used for surface dose measurements. External beam radiationat different gantry angles were delivered to a cubic solid water phantom with a MOSFET placed on the top surface at CAX. The long axis of the MOSFET was oriented along the gantry axis of rotation, with the dosimeter (bubble side) facing the radiation source. MOSFET-measured surface doses were compared against calibrated radiochromic film readings. It was found that both types of MOSFET dosimeters exhibited larger than previously reported angular dependence when measuring surface dose in beams at large oblique angles. For the MOSFET20 dosimeter the measured surface dose deviation against film readings was as high as 17% when the incident angle was 72 degrees to the norm of the phantom surface. It is concluded that some MOSFET dosimeters may have a strong angular dependence when placed on the surface of water-equivalent material, even though they may have an isotropic angular response when surrounded by uniform medium. Extra on-surface calibration maybe necessary before using MOSFET dosimeters for skin dose measurement in tangential fields.

A common approach to implementing the Monte Carlo method for the calculation of brachytherapy radiation dose deposition is to use a phase space file containing information on particles emitted from a brachytherapy source. However, the loading of the phase space file during the dose calculation consumes a large amount of computer random access memory, imposing a higher requirement for computer hardware. In this study, we propose a method to parameterize the information (e.g., particle location, direction and energy) stored in the phase space file by using several probability distributions. This method was implemented for dose calculations of a commercial Ir-192 high dose rate source. Dose calculation accuracy of the parameterized source was compared to the results observed using the full phase space file in a simple water phantom and in a clinical breast cancer case. The results showed the parameterized source at a size of 200 kB was as accurate as the phase space file represented source of 1.1 GB. By using the parameterized source representation, a compact Monte Carlo job can be designed, which allows an easy setup for parallel computing in brachytherapy planning.