THE NATURAL BACKGROUND GAMMA RADIATION EXPOSURE IN THE METROPOLITAN AREA OF THE VALLEY OF MEXICO

Measurements of the natural background radiation have been made at numerous places throughout the world. Very little work in this field has been done in developing countries. In Mexico the natural radiation to which the population is exposed has not been assessed. This dissertation represents a pioneer study in this environmental area. The radiation exposure which occupants within buildings receive as a result of naturally occurring radionuclides present in construction materials is the principal focus.^ Data were collected between August 1979 and November 1980. Continuous monitoring was done with TLDs placed on site for periods of 3 to 6 months. The instrumentation used for "real-time" measurements was a portable NaI (Tl) scintillation detector. In addition, radiometric measurements were performed on construction materials commonly used in Mexican homes.^ Based on TLD readings taken within 75 dwellings, the typical indoor exposure for a resident of the study area is 9.2 (mu)Rh('-1). The average reading of the 152 indoor scintillometer surveys was 9.5 (mu)Rh('-1), the outdoor reading 7.5 (mu)Rh('-1). Results of one-way and multi-way analyses of the exposure data to determine the effect due to building materials type, geologic subsoil, age of dwelling, and elevation are also presented. The results of 152 indoor scintillometer surveys are described. ^

A LAPLACE TRANSFORM PAIR MODEL TO DETERMINE BREMSSTRAHLUNG SPECTRA FROM ATTENUATION DATA (RADIATION, X-RAY, SPECTRAL MODELING)

Every x-ray attenuation curve inherently contains all the information necessary to extract the complete energy spectrum of a beam. To date, attempts to obtain accurate spectral information from attenuation data have been inadequate.^ This investigation presents a mathematical pair model, grounded in physical reality by the Laplace Transformation, to describe the attenuation of a photon beam and the corresponding bremsstrahlung spectral distribution. In addition the Laplace model has been mathematically extended to include characteristic radiation in a physically meaningful way. A method to determine the fraction of characteristic radiation in any diagnostic x-ray beam was introduced for use with the extended model.^ This work has examined the reconstructive capability of the Laplace pair model for a photon beam range of from 50 kVp to 25 MV, using both theoretical and experimental methods.^ In the diagnostic region, excellent agreement between a wide variety of experimental spectra and those reconstructed with the Laplace model was obtained when the atomic composition of the attenuators was accurately known. The model successfully reproduced a 2 MV spectrum but demonstrated difficulty in accurately reconstructing orthovoltage and 6 MV spectra. The 25 MV spectrum was successfully reconstructed although poor agreement with the spectrum obtained by Levy was found.^ The analysis of errors, performed with diagnostic energy data, demonstrated the relative insensitivity of the model to typical experimental errors and confirmed that the model can be successfully used to theoretically derive accurate spectral information from experimental attenuation data. ^

Reduction of setup uncertainties and tumor motion in the radiation therapy of lung tumors

Because the goal of radiation therapy is to deliver a lethal dose to the tumor, accurate information on the location of the tumor needs to be known. Margins are placed around the tumor to account for variations in the daily position of the tumor. If tumor motion and patient setup uncertainties can be reduced, margins that account for such uncertainties in tumor location in can be reduced allowing dose escalation, which in turn could potentially improve survival rates. ^ In the first part of this study, we monitor the location of fiducials implanted in the periphery of lung tumors to determine the extent of non-gated and gated fiducial motion, and to quantify patient setup uncertainties. In the second part we determine where the tumor is when different methods of image-guided patient setup and respiratory gating are employed. In the final part we develop, validate, and implement a technique in which patient setup uncertainties are reduced by aligning patients based upon fiducial locations in projection images. ^ Results from the first part indicate that respiratory gating reduces fiducial motion relative to motion during normal respiration and setup uncertainties when the patients were aligned each day using externally placed skin marks are large. The results from the second part indicate that current margins that account for setup uncertainty and tumor motion result in less than 2% of the tumor outside of the planning target volume (PTV) when the patient is aligned using skin marks. In addition, we found that if respiratory gating is going to be used, it is most effective if used in conjunction with image-guided patient setup. From the third part, we successfully developed, validated, and implemented on a patient a technique for aligning a moving target prior to treatment to reduce the uncertainties in tumor location. ^ In conclusion, setup uncertainties and tumor motion are a significant problem when treating tumors located within the thoracic region. Image-guided patient setup in conjunction with treatment delivery using respiratory gating reduces these uncertainties in tumor locations. In doing so, margins around the tumor used to generate the PTV can be reduced, which may allow for dose escalation to the tumor. ^

Electron intensity modulation for mixed-beam radiation therapy with an x-ray multi-leaf collimator

The current standard treatment for head and neck cancer at our institution uses intensity-modulated x-ray therapy (IMRT), which improves target coverage and sparing of critical structures by delivering complex fluence patterns from a variety of beam directions to conform dose distributions to the shape of the target volume. The standard treatment for breast patients is field-in-field forward-planned IMRT, with initial tangential fields and additional reduced-weight tangents with blocking to minimize hot spots. For these treatment sites, the addition of electrons has the potential of improving target coverage and sparing of critical structures due to rapid dose falloff with depth and reduced exit dose. In this work, the use of mixed-beam therapy (MBT), i.e., combined intensity-modulated electron and x-ray beams using the x-ray multi-leaf collimator (MLC), was explored. The hypothesis of this study was that addition of intensity-modulated electron beams to existing clinical IMRT plans would produce MBT plans that were superior to the original IMRT plans for at least 50% of selected head and neck and 50% of breast cases. Dose calculations for electron beams collimated by the MLC were performed with Monte Carlo methods. An automation system was created to facilitate communication between the dose calculation engine and the treatment planning system. Energy and intensity modulation of the electron beams was accomplished by dividing the electron beams into 2x2-cm2 beamlets, which were then beam-weight optimized along with intensity-modulated x-ray beams. Treatment plans were optimized to obtain equivalent target dose coverage, and then compared with the original treatment plans. MBT treatment plans were evaluated by participating physicians with respect to target coverage, normal structure dose, and overall plan quality in comparison with original clinical plans. The physician evaluations did not support the hypothesis for either site, with MBT selected as superior in 1 out of the 15 head and neck cases (p=1) and 6 out of 18 breast cases (p=0.95). While MBT was not shown to be superior to IMRT, reductions were observed in doses to critical structures distal to the target along the electron beam direction and to non-target tissues, at the expense of target coverage and dose homogeneity. ^

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. ^

A convolution model for energy transport in a therapeutic fast neutron beam

A three-dimensional model has been proposed that uses Monte Carlo and fast Fourier transform convolution techniques to calculate the dose distribution from a fast neutron beam. This method transports scattered neutrons and photons in the forward, lateral, and backward directions and protons, electrons, and positrons in the forward and lateral directions by convolving energy spread kernels with initial interaction available energy distributions. The primary neutron and photon spectrums have been derived from narrow beam attenuation measurements. The positions and strengths of the effective primary neutron, scattered neutron, and photon sources have been derived from dual ion chamber measurements. The size of the effective primary neutron source has been measured using a copper activation technique. Heterogeneous tissue calculations require a weighted sum of two convolutions for each component since the kernels must be invariant for FFT convolution. Comparisons between calculations and measurements were performed for several water and heterogeneous phantom geometries. ^

A dosimetric study of radionuclide therapy for bone marrow ablation

In a phase I clinical trial, six multiple myeloma patients, who were non-responsive to conventional therapy and were scheduled for bone marrow transplantation, received Holmium-166 ($\sp{166}$Ho) labeled to a bone seeking agent, DOTMP (1,4,7,10-tetraazacyclododecane-1,4,7,10-tetramethylene-phosphonic acid), for the purpose of bone marrow ablation. The specific aims of my research within this protocol were to evaluate the toxicity and efficacy of $\sp{166}$Ho DOTMP by quantifying the in vivo pharmacokinetics and radiation dosimetry, and by correlating these results to the biologic response observed. The reproducibility of pharmacokinetics from multiple injections of $\sp{166}$Ho DOTMP administered to these myeloma patients was demonstrated from both blood and whole body retention. The skeletal concentration of $\sp{166}$Ho DOTMP was heterogenous in all six patients: high in the ribs, pelvis, and lumbar vertebrae regions, and relatively low in the femurs, arms, and head.^ A novel technique was developed to calculate the radiation dose to the bone marrow in each skeletal ROI, and was applied to all six $\sp{166}$Ho DOTMP patients. Radiation dose estimates for the bone marrow calculated using the standard MIRD "S" factors were compared with the average values derived from the heterogenous distribution of activity in the skeleton (i.e., the regional technique). The results from the two techniques were significantly different; the average of the dose estimates from the regional technique were typically 30% greater. Furthermore, the regional technique provided a range of radiation doses for the entire marrow volume, while the MIRD "S" factors only provided a single value. Dose volume histogram analysis of data from the regional technique indicated a range of dose estimates that varied by a factor of 10 between the high dose and low dose regions. Finally, the observed clinical response of cells and abnormal proteins measured in bone marrow aspirates and peripheral blood samples were compared with radiation dose estimates for the bone marrow calculated from the standard and regional technique. The results showed the regional technique values correlated more closely to several clinical response parameters. (Abstract shortened by UMI.) ^

Development of automated image analysis tools for verification of radiotherapy field accuracy with an electronic portal imaging device

The successful management of cancer with radiation relies on the accurate deposition of a prescribed dose to a prescribed anatomical volume within the patient. Treatment set-up errors are inevitable because the alignment of field shaping devices with the patient must be repeated daily up to eighty times during the course of a fractionated radiotherapy treatment. With the invention of electronic portal imaging devices (EPIDs), patient's portal images can be visualized daily in real-time after only a small fraction of the radiation dose has been delivered to each treatment field. However, the accuracy of human visual evaluation of low-contrast portal images has been found to be inadequate. The goal of this research is to develop automated image analysis tools to detect both treatment field shape errors and patient anatomy placement errors with an EPID. A moments method has been developed to align treatment field images to compensate for lack of repositioning precision of the image detector. A figure of merit has also been established to verify the shape and rotation of the treatment fields. Following proper alignment of treatment field boundaries, a cross-correlation method has been developed to detect shifts of the patient's anatomy relative to the treatment field boundary. Phantom studies showed that the moments method aligned the radiation fields to within 0.5mm of translation and 0.5$\sp\circ$ of rotation and that the cross-correlation method aligned anatomical structures inside the radiation field to within 1 mm of translation and 1$\sp\circ$ of rotation. A new procedure of generating and using digitally reconstructed radiographs (DRRs) at megavoltage energies as reference images was also investigated. The procedure allowed a direct comparison between a designed treatment portal and the actual patient setup positions detected by an EPID. Phantom studies confirmed the feasibility of the methodology. Both the moments method and the cross-correlation technique were implemented within an experimental radiotherapy picture archival and communication system (RT-PACS) and were used clinically to evaluate the setup variability of two groups of cancer patients treated with and without an alpha-cradle immobilization aid. The tools developed in this project have proven to be very effective and have played an important role in detecting patient alignment errors and field-shape errors in treatment fields formed by a multileaf collimator (MLC). ^

The volume effect in irradiated mouse colorectum

Damage of the colorectum is the dose-limiting normal tissue complication following radiotherapy of prostate and cervical cancers. One approach for decreasing complications is to physically reduce the treatment volume. Mathematical models have been previously developed to describe the change in associated toxicity with a change in irradiated volume, i.e. the "volume effect", for serial-type normal tissues including the colorectum. The first goal of this thesis was to test the hypothesis that there would not be a threshold length in the development of obstruction after irradiation of mouse colorectum, as predicted by the Probability model of the volume effect. The second goal was to examine if there were differences in the threshold and in the incidence of colorectal obstruction after irradiation of two mouse strains, C57B1/6 (C57) and C3Hf/Kam (C3H), previously found to be fibrosis-prone and-resistant, respectively, after lung irradiation due, in part, to genetic differences. The hypothesis examined was that differences in incidence between strains were due to the differential expression of the fibrogenic cytokines $\rm TGF\beta$ and $\rm TNF\alpha.$ Various lengths of C57 and C3H mouse colorectum were irradiated and the incidence of colorectal obstruction was followed up to 15 months. A threshold length was observed for both mouse strains, in contradiction of model predictions. The mechanism of the threshold was epithelial regeneration after irradiation. C57 mice had significantly higher incidence of colorectal obstruction compared to C3H mice, especially at smaller irradiated lengths. Colorectal tissue was obtained at various times after irradiation and prepared for histology, immunohistochemistry and RNase protection assay for measurement of $\rm TGF\beta 1,$ 2, 3 and $\rm TNF\alpha$ mRNA. Distinct strain differences in the histological time of appearance and spatial locations of fibrosis were observed. However, there were no consistent strain difference in mRNA levels or immunolocalization for any of the cytokines examined. The data indicate the need for volume effect models that account for biologically important processes, such as the effect of epithelial regeneration after irradiation. As well, changes in fibrogenic cytokines at the mRNA level do not contribute to the strain difference in radiation-induced colorectal obstruction. ^

STUDY OF THE EFFECTS OF TOTAL MODULATION TRANSFER FUNCTION CHANGES ON OBSERVER PERFORMANCE USING CLINICAL MAMMOGRAMS

The main goal of this study was to relate physical changes in image quality measured by Modulation Transfer Function (MTF) to diagnostic accuracy.^ One Hundred and Fifty Kodak Min-R screen/film combination conventional craniocaudal mammograms obtained with the Pfizer Microfocus Mammographic system were selected from the files of the Department of Radiology, at M.D. Anderson Hospital and Tumor Institute.^ The mammograms included 88 cases with a variety of benign diagnosis and 62 cases with a variety of malignant biopsy diagnosis. The average age of the patient population was 55 years old. 70 cases presented calcifications with 30 cases having calcifications smaller than 0.5mm. 46 cases presented irregular bordered masses larger than 1 cm. 30 cases presented smooth bordered masses with 20 larger than 1 cm.^ Four separated copies of the original images were made each having a different change in the MTF using a defocusing technique whereby copies of the original were obtained by light exposure through different thicknesses (spacing) of transparent film base.^ The mammograms were randomized, and evaluated by three experienced mammographers for the degree of visibility of various anatomical breast structures and pathological lesions (masses and calicifications), subjective image quality, and mammographic interpretation.^ 3,000 separate evaluations were anayzed by several statistical techniques including Receiver Operating Characteristic curve analysis, McNemar test for differences between proportions and the Landis et al. method of agreement weighted kappa for ordinal categorical data.^ Results from the statistical analysis show: (1) There were no statistical significant differences in the diagnostic accuracy of the observers when diagnosing from mammograms with the same MTF. (2) There were no statistically significant differences in diagnostic accuracy for each observer when diagnosing from mammograms with the different MTF's used in the study. (3) There statistical significant differences in detail visibility between the copies and the originals. Detail visibility was better in the originals. (4) Feature interpretations were not significantly different between the originals and the copies. (5) Perception of image quality did not affect image interpretation.^ Continuation and improvement of this research ca be accomplished by: using a case population more sensitive to MTF changes, i.e., asymptomatic women with minimum breast cancer, more observers (including less experienced radiologists and experienced technologists) must collaborate in the study, and using a minimum of 200 benign and 200 malignant cases.^

A novel treatment planning methodology for high dose 166Ho-DOTMP therapy in patients with multiple myeloma

Bone marrow ablation, i.e., the complete sterilization of the active bone marrow, followed by bone marrow transplantation (BMT) is a comment treatment of hematological malignancies. The use of targeted bone-seeking radiopharmaceuticals to selectively deliver radiation to the adjacent bone marrow cavities while sparing normal tissues is a promising technique. Current radiopharmaceutical treatment planning methods do not properly compensate for the patient-specific variable distribution of radioactive material within the skeleton. To improve the current method of internal dosimetry, novel methods for measuring the radiopharmaceutical distribution within the skeleton were developed. 99mTc-MDP was proven as an adequate surrogate for measuring 166Ho-DOTMP skeletal uptake and biodistribution, allowing these measures to be obtained faster, safer, and with higher spatial resolution. This translates directly into better measurements of the radiation dose distribution within the bone marrow. The resulting bone marrow dose-volume histograms allow prediction of the patient disease response where conventional organ scale dosimetry failed. They indicate that complete remission is only achieved when greater than 90% of the bone marrow receives at least 30 Gy. ^ Comprehensive treatment planning requires combining target and non-target organ dosimetry. Organs in the urinary tract were of special concern. The kidney dose is primarily dependent upon the mean transit time of 166 Ho-DOTMP through the kidney. Deconvolution analysis of renograms predicted a mean transit time of 2.6 minutes for 166Ho-DOTMP. The radiation dose to the urinary bladder wall is dependent upon numerous factors including patient hydration and void schedule. For beta-emitting isotopes such as 166Ho, reduction of the bladder wall dose is best accomplished through good patient hydration and ensuring a partially full bladder at the time of injection. Encouraging the patient to void frequently, or catheterizing the patient without irrigation, will not significantly reduce the bladder wall dose. ^ The results from this work will produce the most advanced treatment planning methodology for bone marrow ablation therapy using radioisotopes currently available. Treatments can be tailored specifically for each patient, including the addition of concomitant total body irradiation for patients with unfavorable dose distributions, to deliver a desired patient disease response, while minimizing the dose or toxicity to non-target organs. ^

Fast magnetic resonance temperature imaging for focused ultrasound thermal therapy

The current standard for temperature sensitive imaging using magnetic resonance (MR) is 2-D, spoiled, fast gradient-echo (fGRE) phase-difference imaging exploiting temperature dependent changes in the proton resonance frequency (PRF). The echo-time (TE) for optimal sensitivity is larger than the typical repetition time (TR) of an fGRE sequence. Since TE must be less than TR in the fGRE sequence, this limits the technique's achievable sensitivity, spatial, and temporal resolution. This adversely affects both accuracy and volume coverage of the measurements. Accurate measurement of the rapid temperature changes associated with pulsed thermal therapies, such as high-intensity focused ultrasound (FUS), at optimal temperature sensitivity requires faster acquisition times than those currently available. ^ Use of fast MR acquisition strategies, such as interleaved echo-planar and spiral imaging, can provide the necessary increase in temporal performance and sensitivity while maintaining adequate signal-to-noise and in-plane spatial resolution. This research explored the adaptation and optimization of several fast MR acquisition methods for thermal monitoring of pulsed FUS thermal therapy. Temperature sensitivity, phase-difference noise and phase-difference to phase-difference-to noise ratio for the different pulse sequences were evaluated under varying imaging parameters in an agar gel phantom to establish optimal sequence parameters for temperature monitoring. The temperature sensitivity coefficient of the gel phantom was measured, allowing quantitative temperature extrapolations. ^ Optimized fast sequences were compared based on the ability to accurately monitor temperature changes at the focus of a high-intensity focused ultrasound unit, volume coverage, and contrast-to-noise ratio in the temperature maps. Operating parameters, which minimize complex phase-difference measurement errors introduced by use of the fast-imaging methods, were established. ^

Imaging properties of scanning equalization digital radiography: A simulation study

We have developed an empirically based simulation system to create images equivalent in SNR and SPR to those that would be acquired with various possible SEDR configurations. This system uses a collection of spot collimated full-field images (SCFFIs) of an anthropomorphic chest phantom, taken at high exposure levels and rescaled in noise and intensity, then digitally collimated and combined to produce the simulated SEDR images. This system allows for the study of design trade-offs between different equalization feedback schemes and scatter rejection geometries in addition to estimating the clinical benefits of SEDR over traditional imaging techniques. Data from this simulation system has demonstrated that SEDR techniques offer potential significant improvements over currently used digital radiography techniques for chest imaging. ^