Characterization of Optically Stimulated Luminescent Detectors in Photon & Proton Beams for Use in Anthropomorphic Phantoms

This study investigated characteristics of optically stimulated luminescent detectors (OSLDs) in protons, allowing comparison to thermoluminescent detectors, and to be implemented into the Radiological Physics Center’s (RPC) remote audit quality assurance program for protons, and for remote anthropomorphic phantom irradiations. The OSLDs used were aluminum oxide (Al2O3:C) nanoDots from Landauer, Inc. (Glenwood, Ill.) measuring 10x10x2 mm3. A square, 20(L)x20(W)x0.5(H) cm3 piece of solid water was fabricated with pockets to allow OSLDs and TLDs to be irradiated simultaneously and perpendicular to the beam. Irradiations were performed at 5cm depth in photons, and in the center of a 10 cm SOBP in a 200MeV proton beam. Additionally, the Radiological Physics Center’s anthropomorphic pelvic phantom was used to test the angular dependence of OSLDs in photons and protons. A cylindrical insert in the phantom allows the dosimeters to be rotated to any angle with a fixed gantry angle. OSLDs were irradiated at 12 angles between 0 and 360 degrees. The OSLDs were read out with a MicroStar reader from Landauer, Inc. Dose response indicates that at angles where the dosimeter is near parallel with the radiation beam response is reduced slightly. Measurements in proton beams do not show significant angular dependence. Post-irradiation fading of OSLDs was studied in proton beams to determine if the fading was different than that of photons. The fading results showed no significant difference from results in photon beams. OSLDs and TLDs are comparable within 3% in photon beams and a correction factor can be posited for proton beams. With angular dependence characteristics defined, OSLDs can be implemented into multiple-field treatment plans in photons and protons and used in the RPC’s quality assurance program.

NOVEL PHANTOMS AND POST-PROCESSING FOR DIFFUSION SPECTRUM IMAGING

High Angular Resolution Diffusion Imaging (HARDI) techniques, including Diffusion Spectrum Imaging (DSI), have been proposed to resolve crossing and other complex fiber architecture in the human brain white matter. In these methods, directional information of diffusion is inferred from the peaks in the orientation distribution function (ODF). Extensive studies using histology on macaque brain, cat cerebellum, rat hippocampus and optic tracts, and bovine tongue are qualitatively in agreement with the DSI-derived ODFs and tractography. However, there are only two studies in the literature which validated the DSI results using physical phantoms and both these studies were not performed on a clinical MRI scanner. Also, the limited studies which optimized DSI in a clinical setting, did not involve a comparison against physical phantoms. Finally, there is lack of consensus on the necessary pre- and post-processing steps in DSI; and ground truth diffusion fiber phantoms are not yet standardized. Therefore, the aims of this dissertation were to design and construct novel diffusion phantoms, employ post-processing techniques in order to systematically validate and optimize (DSI)-derived fiber ODFs in the crossing regions on a clinical 3T MR scanner, and develop user-friendly software for DSI data reconstruction and analysis. Phantoms with a fixed crossing fiber configuration of two crossing fibers at 90° and 45° respectively along with a phantom with three crossing fibers at 60°, using novel hollow plastic capillaries and novel placeholders, were constructed. T2-weighted MRI results on these phantoms demonstrated high SNR, homogeneous signal, and absence of air bubbles. Also, a technique to deconvolve the response function of an individual peak from the overall ODF was implemented, in addition to other DSI post-processing steps. This technique greatly improved the angular resolution of the otherwise unresolvable peaks in a crossing fiber ODF. The effects of DSI acquisition parameters and SNR on the resultant angular accuracy of DSI on the clinical scanner were studied and quantified using the developed phantoms. With a high angular direction sampling and reasonable levels of SNR, quantification of a crossing region in the 90°, 45° and 60° phantoms resulted in a successful detection of angular information with mean ± SD of 86.93°±2.65°, 44.61°±1.6° and 60.03°±2.21° respectively, while simultaneously enhancing the ODFs in regions containing single fibers. For the applicability of these validated methodologies in DSI, improvement in ODFs and fiber tracking from known crossing fiber regions in normal human subjects were demonstrated; and an in-house software package in MATLAB which streamlines the data reconstruction and post-processing for DSI, with easy to use graphical user interface was developed. In conclusion, the phantoms developed in this dissertation offer a means of providing ground truth for validation of reconstruction and tractography algorithms of various diffusion models (including DSI). Also, the deconvolution methodology (when applied as an additional DSI post-processing step) significantly improved the angular accuracy of the ODFs obtained from DSI, and should be applicable to ODFs obtained from the other high angular resolution diffusion imaging techniques.

Development and Implementation of the Use of Optically Stimulated Luminescent Detectors in the Radiological Physics Center Anthropomorphic Quality Assurance Phantoms

The Radiological Physics Center (RPC) uses both on-site and remote reviews to credential institutions for participation in clinical trials. Anthropomorphic quality assurance (QA) phantoms are one tool the RPC uses to remotely audit institutions, which include thermoluminescent dosimeters (TLDs) and radiochromic film. The RPC desires to switch from TLD as the absolute dosimeter in the phantoms, to optically stimulated luminescent dosimeters (OSLDs), but a problem lies in the angular dependence exhibited by the OSLD. The purpose of this study was to characterize the angular dependence of OSLD and establish a correction factor if necessary, to provide accurate dosimetric measurements as a replacement for TLD in the QA phantoms. A 10 cm diameter high-impact polystyrene spherical phantom was designed and constructed to hold an OSLD to study the angular response of the dosimeter under the simplest of circumstances for both coplanar and non-coplanar treatment deliveries. OSLD were irradiated in the spherical phantom, and the responses of the dosimeter from edge-on angles were normalized to the response when irradiated with the beam incident normally on the surface of the dosimeter. The average normalized response was used to establish an angular correction factor for 6 MV and 18 coplanar treatments, and for 6 MV non-coplanar treatments specific to CyberKnife. The RPC pelvic phantom dosimetry insert was modified to hold OSLD, in addition to the TLD, adjacent to the planes of film. Treatment plans of increasing angular beam delivery were developed, three in Pinnacle v9.0 (4-field box, IMRT, and VMAT) and one in Accuray’s MultiPlan v3.5.3 (CyberKnife). The plans were delivered to the pelvic phantom containing both TLD and OSLD in the target volume. The pelvic phantom was also sent to two institutions to be irradiated as trials, one delivering IMRT, and the other a CyberKnife treatment. For the IMRT deliveries and the two institution trials, the phantom also included film in the sagittal and coronal planes. The doses measured from the TLD and OSLD were calculated for each irradiation, and the angular correction factors established from the spherical phantom irradiations were applied to the OSLD dose. The ratio of the TLD dose to the angular corrected OSLD dose was calculated for each irradiation. The corrected OSLD dose was found to be within 1% of the TLD measured dose for all irradiations, with the exception of the in-house CyberKnife deliveries. The films were normalized to both TLD measured dose and the corrected OSLD dose. Dose profiles were obtained and gamma analysis was performed using a 7%/4 mm criteria, to compare the ability of the OSLD, when corrected for the angular dependence, to provide equivalent results to TLD. The results of this study indicate that the OSLD can effectively be used as a replacement for TLD in the RPC’s anthropomorphic QA phantoms for coplanar treatment deliveries when a correction is applied for the dosimeter’s angular dependence.

Implementation of an Anthropomorphic Phantom for the Evaluation of Proton Therapy Treatment Procedures

With an increasing number of institutions offering proton therapy, the number of multi-institutional clinical trials involving proton therapy will also increase in the coming years. The Radiological Physics Center monitors sites involved in clinical trials through the use of site visits and remote auditing with thermoluminescent dosimeters (TLD) and mailable anthropomorphic phantoms. Currently, there are no heterogeneous phantoms that have been commissioned to evaluate proton therapy. It was hypothesized that an anthropomorphic pelvis phantom can be designed to audit treatment procedures (patient simulation, treatment planning and treatment delivery) at proton facilities to confirm agreement between the measured dose and calculated dose within 5%/3mm with a reproducibility of 3%. A pelvis phantom originally designed for use with photon treatments was retrofitted for use in proton therapy. The relative stopping power (SP) of each phantom material was measured. Hounsfield Units (HU) for each phantom material were measured with a CT scanner and compared to the relative stopping power calibration curve. The tissue equivalency for each material was calculated. Two proton treatment plans were created; one which did not correct for material SP differences (Plan 1) and one plan which did correct for SP differences (Plan 2). Film and TLD were loaded into the phantom and the phantom was irradiated 3 times per plan. The measured values were compared to the HU-SP calibration curve and it was found that the stopping powers for the materials could be underestimated by 5-10%. Plan 1 passed the criteria for the TLD and film margins with reproducibility under 3% between the 3 trials. Plan 2 failed because the right-left film dose profile average displacement was -9.0 mm on the left side and 6.0 mm on the right side. Plan 2 was intended to improve the agreements and instead introduced large displacements along the path of the beam. Plan 2 more closely represented the actual phantom composition with corrected stopping powers and should have shown an agreement between the measured and calculated dose within 5%/3mm. The hypothesis was rejected and the pelvis phantom was found to be not suitable to evaluate proton therapy treatment procedures.

Commissioning an Anthropomorphic Spine and Lung Phantom for Remote Dose Verification of Institutions Participating in RTOG 0631

The RPC developed a new phantom to ensure comparable and consistent radiation administration in spinal radiosurgery clinical trials. This study assessed the phantom’s dosimetric and anatomic utility. The ‘spine phantom’ is a water filled thorax with anatomy encountered in spinal radiosurgery: target volume, vertebral column, spinal canal, esophagus, heart, and lungs. The dose to the target volume was measured with axial and sagittal planes of radiochromic film and thermoluminescent dosimeters (TLD). The dose distributions were measured with the radiochromic film calibrated to the absolute dose measured by the TLD. Four irradiations were administered: a four angle box plan, a seven angle conformal plan, a seven angle IMRT plan, and a nine angle IMRT plan (denoted as IMRT plan #1 and plan #2, respectively). In each plan, at least 95% of the defined tumor volume received 8 Gy. For each irradiation the planned and administered dose distributions were registered via pinpricks, and compared using point dose measurements, dose profiles, isodose distributions, and gamma analyses. Based on previous experience at the RPC, a gamma analysis was considering passing if greater than 95% of pixels passed the criteria of 5% dose difference and 3 mm distance-to-agreement. Each irradiation showed acceptable agreement in the qualitative assessments and exceeded the 95% passing rate at the 5% / 3 mm criteria, except IMRT plan #1, which was determined to have been poorly localized during treatment administration. The measured and planned dose distributions demonstrated acceptable agreement at the 5% / 3 mm criteria, and the spine phantom was determined to be a useful tool for the remote assessment of an institution’s treatment planning and dose delivery regimen.

EVALUATION OF THE QUANTITATIVE ACCURACY OF A COMMERCIALLY-AVAILABLE POSITRON EMISSION MAMMOGRAPHY SCANNER

Objective: The PEM Flex Solo II (Naviscan, Inc., San Diego, CA) is currently the only commercially-available positron emission mammography (PEM) scanner. This scanner does not apply corrections for count rate effects, attenuation or scatter during image reconstruction, potentially affecting the quantitative accuracy of images. This work measures the overall quantitative accuracy of the PEM Flex system, and determines the contributions of error due to count rate effects, attenuation and scatter. Materials and Methods: Gelatin phantoms were designed to simulate breasts of different sizes (4 – 12 cm thick) with varying uniform background activity concentration (0.007 – 0.5 μCi/cc), cysts and lesions (2:1, 5:1, 10:1 lesion-to-background ratios). The overall error was calculated from ROI measurements in the phantoms with a clinically relevant background activity concentration (0.065 μCi/cc). The error due to count rate effects was determined by comparing the overall error at multiple background activity concentrations to the error at 0.007 μCi/cc. A point source and cold gelatin phantoms were used to assess the errors due to attenuation and scatter. The maximum pixel values in gelatin and in air were compared to determine the effect of attenuation. Scatter was evaluated by comparing the sum of all pixel values in gelatin and in air. Results: The overall error in the background was found to be negative in phantoms of all thicknesses, with the exception of the 4-cm thick phantoms (0%±7%), and it increased with thickness (-34%±6% for the 12-cm phantoms). All lesions exhibited large negative error (-22% for the 2:1 lesions in the 4-cm phantom) which increased with thickness and with lesion-to-background ratio (-85% for the 10:1 lesions in the 12-cm phantoms). The error due to count rate in phantoms with 0.065 μCi/cc background was negative (-23%±6% for 4-cm thickness) and decreased with thickness (-7%±7% for 12 cm). Attenuation was a substantial source of negative error and increased with thickness (-51%±10% to -77% ±4% in 4 to 12 cm phantoms, respectively). Scatter contributed a relatively constant amount of positive error (+23%±11%) for all thicknesses. Conclusion: Applying corrections for count rate, attenuation and scatter will be essential for the PEM Flex Solo II to be able to produce quantitatively accurate images.

BENCHMARKING AND IMPLEMENTATION OF A NEW INDEPENDENT MONTE CARLO DOSE CALCULATION QUALITY ASSURANCE AUDIT TOOL FOR CLINICAL TRIALS

Introduction Commercial treatment planning systems employ a variety of dose calculation algorithms to plan and predict the dose distributions a patient receives during external beam radiation therapy. Traditionally, the Radiological Physics Center has relied on measurements to assure that institutions participating in the National Cancer Institute sponsored clinical trials administer radiation in doses that are clinically comparable to those of other participating institutions. To complement the effort of the RPC, an independent dose calculation tool needs to be developed that will enable a generic method to determine patient dose distributions in three dimensions and to perform retrospective analysis of radiation delivered to patients who enrolled in past clinical trials. Methods A multi-source model representing output for Varian 6 MV and 10 MV photon beams was developed and evaluated. The Monte Carlo algorithm, know as the Dose Planning Method (DPM), was used to perform the dose calculations. The dose calculations were compared to measurements made in a water phantom and in anthropomorphic phantoms. Intensity modulated radiation therapy and stereotactic body radiation therapy techniques were used with the anthropomorphic phantoms. Finally, past patient treatment plans were selected and recalculated using DPM and contrasted against a commercial dose calculation algorithm. Results The multi-source model was validated for the Varian 6 MV and 10 MV photon beams. The benchmark evaluations demonstrated the ability of the model to accurately calculate dose for the Varian 6 MV and the Varian 10 MV source models. The patient calculations proved that the model was reproducible in determining dose under similar conditions described by the benchmark tests. Conclusions The dose calculation tool that relied on a multi-source model approach and used the DPM code to calculate dose was developed, validated, and benchmarked for the Varian 6 MV and 10 MV photon beams. Several patient dose distributions were contrasted against a commercial algorithm to provide a proof of principal to use as an application in monitoring clinical trial activity.

MODIFICATION AND IMPLEMENTATION OF THE RPC HETEROGENEOUS THORAX PHANTOM FOR VERIFICATION OF PROTON THERAPY TREATMENT PROCEDURES

The Radiological Physics Center (RPC) provides heterogeneous phantoms that are used to evaluate radiation treatment procedures as part of a comprehensive quality assurance program for institutions participating in clinical trials. It was hypothesized that the existing RPC heterogeneous thorax phantom can be modified to assess lung tumor proton beam therapy procedures involving patient simulation, treatment planning, and treatment delivery, and could confirm agreement between the measured dose and calculated dose within 5%/3mm with a reproducibility of 5%. The Hounsfield Units (HU) for lung equivalent materials (balsa wood and cork) was measured using a CT scanner. The relative linear stopping power (RLSP) of these materials was measured. The linear energy transfer (LET) of Gafchromic EBT2 film was analyzed utilizing parallel and perpendicular orientations in a water tank and compared to ion chamber readings. Both parallel and perpendicular orientations displayed a quenching effect underperforming the ion chamber, with the parallel orientation showing an average 31 % difference and the perpendicular showing an average of 15% difference. Two treatment plans were created that delivered the prescribed dose to the target volume, while achieving low entrance doses. Both treatment plans were designed using smeared compensators and expanded apertures, as would be utilized for a patient in the clinic. Plan 1a contained two beams that were set to orthogonal angles and a zero degree couch kick. Plan 1b utilized two beams set to 10 and 80 degrees with a 15 degree couch kick. EBT2 film and TLD were inserted and the phantom was irradiated 3 times for each plan. Both plans passed the criteria for the TLD measurements where the TLD values were within 7% of the dose calculated by Eclipse. Utilizing the 5%/3mm criteria, the 3 trial average of overall pass rate was 71% for Plan 1a. The 3 trial average for the overall pass rate was 76% for Plan 1b. The trials were then analyzed using RPC conventional lung treatment guidelines set forth by the RTOG: 5%/5mm, and an overall pass rate of 85%. Utilizing these criteria, only Plan 1b passed for all 3 trials, with an average overall pass rate of 89%.

Perfusion in rat brain at 7 T with arterial spin labeling using FAIR-TrueFISP and QUIPSS.

Measurement of perfusion in longitudinal studies allows for the assessment of tissue integrity and the detection of subtle pathologies. In this work, the feasibility of measuring brain perfusion in rats with high spatial resolution using arterial spin labeling is reported. A flow-sensitive alternating recovery sequence, coupled with a balanced gradient fast imaging with steady-state precession readout section was used to minimize ghosting and geometric distortions, while achieving high signal-to-noise ratio. The quantitative imaging of perfusion using a single subtraction method was implemented to address the effects of variable transit delays between the labeling of spins and their arrival at the imaging slice. Studies in six rats at 7 T showed good perfusion contrast with minimal geometric distortion. The measured blood flow values of 152.5+/-6.3 ml/100 g per minute in gray matter and 72.3+/-14.0 ml/100 g per minute in white matter are in good agreement with previously reported values based on autoradiography, considered to be the gold standard.

Three-dimensional electron beam dose calculations

The MDAH pencil-beam algorithm developed by Hogstrom et al (1981) has been widely used in clinics for electron beam dose calculations for radiotherapy treatment planning. The primary objective of this research was to address several deficiencies of that algorithm and to develop an enhanced version. Two enhancements have been incorporated into the pencil-beam algorithm; one models fluence rather than planar fluence, and the other models the bremsstrahlung dose using measured beam data. Comparisons of the resulting calculated dose distributions with measured dose distributions for several test phantoms have been made. From these results it is concluded (1) that the fluence-based algorithm is more accurate to use for the dose calculation in an inhomogeneous slab phantom, and (2) the fluence-based calculation provides only a limited improvement to the accuracy the calculated dose in the region just downstream of the lateral edge of an inhomogeneity. The source of the latter inaccuracy is believed primarily due to assumptions made in the pencil beam's modeling of the complex phantom or patient geometry.^ A pencil-beam redefinition model was developed for the calculation of electron beam dose distributions in three dimensions. The primary aim of this redefinition model was to solve the dosimetry problem presented by deep inhomogeneities, which was the major deficiency of the enhanced version of the MDAH pencil-beam algorithm. The pencil-beam redefinition model is based on the theory of electron transport by redefining the pencil beams at each layer of the medium. The unique approach of this model is that all the physical parameters of a given pencil beam are characterized for multiple energy bins. Comparisons of the calculated dose distributions with measured dose distributions for a homogeneous water phantom and for phantoms with deep inhomogeneities have been made. From these results it is concluded that the redefinition algorithm is superior to the conventional, fluence-based, pencil-beam algorithm, especially in predicting the dose distribution downstream of a local inhomogeneity. The accuracy of this algorithm appears sufficient for clinical use, and the algorithm is structured for future expansion of the physical model if required for site specific treatment planning problems. ^