The use of electronic portal imaging to verify patient position during intensity-modulated radiotherapy delivered by the dynamic MLC technique

Purpose: The precise shape of the three-dimensional dose distributions created by intensity-modulated radiotherapy means that the verification of patient position and setup is crucial to the outcome of the treatment. In this paper, we investigate and compare the use of two different image calibration procedures that allow extraction of patient anatomy from measured electronic portal images of intensity-modulated treatment beams. Methods and Materials: Electronic portal images of the intensity-modulated treatment beam delivered using the dynamic multileaf collimator technique were acquired. The images were formed by measuring a series of frames or segments throughout the delivery of the beams. The frames were then summed to produce an integrated portal image of the delivered beam. Two different methods for calibrating the integrated image were investigated with the aim of removing the intensity modulations of the beam. The first involved a simple point-by-point division of the integrated image by a single calibration image of the intensity-modulated beam delivered to a homogeneous polymethyl methacrylate (PMMA) phantom. The second calibration method is known as the quadratic calibration method and required a series of calibration images of the intensity-modulated beam delivered to different thicknesses of homogeneous PMMA blocks. Measurements were made using two different detector systems: a Varian amorphous silicon flat-panel imager and a Theraview camera-based system. The methods were tested first using a contrast phantom before images were acquired of intensity-modulated radiotherapy treatment delivered to the prostate and pelvic nodes of cancer patients at the Royal Marsden Hospital. Results: The results indicate that the calibration methods can be used to remove the intensity modulations of the beam, making it possible to see the outlines of bony anatomy that could be used for patient position verification. This was shown for both posterior and lateral delivered fields. Conclusions: Very little difference between the two calibration methods was observed, so the simpler division method, requiring only the single extra calibration measurement and much simpler computation, was the favored method. This new method could provide a complementary tool to existing position verification methods, and it has the advantage that it is completely passive, requiring no further dose to the patient and using only the treatment fields.

The development of Monte Carlo techniques for the verification of radiotherapy treatments

Introduction: Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo (MC) methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the ‘gold standard’ for predicting dose deposition in the patient [1]. This project has three main aims: 1. To develop tools that enable the transfer of treatment plan information from the treatment planning system (TPS) to a MC dose calculation engine. 2. To develop tools for comparing the 3D dose distributions calculated by the TPS and the MC dose engine. 3. To investigate the radiobiological significance of any errors between the TPS patient dose distribution and the MC dose distribution in terms of Tumour Control Probability (TCP) and Normal Tissue Complication Probabilities (NTCP). The work presented here addresses the first two aims. Methods: (1a) Plan Importing: A database of commissioned accelerator models (Elekta Precise and Varian 2100CD) has been developed for treatment simulations in the MC system (EGSnrc/BEAMnrc). Beam descriptions can be exported from the TPS using the widespread DICOM framework, and the resultant files are parsed with the assistance of a software library (PixelMed Java DICOM Toolkit). The information in these files (such as the monitor units, the jaw positions and gantry orientation) is used to construct a plan-specific accelerator model which allows an accurate simulation of the patient treatment field. (1b) Dose Simulation: The calculation of a dose distribution requires patient CT images which are prepared for the MC simulation using a tool (CTCREATE) packaged with the system. Beam simulation results are converted to absolute dose per- MU using calibration factors recorded during the commissioning process and treatment simulation. These distributions are combined according to the MU meter settings stored in the exported plan to produce an accurate description of the prescribed dose to the patient. (2) Dose Comparison: TPS dose calculations can be obtained using either a DICOM export or by direct retrieval of binary dose files from the file system. Dose difference, gamma evaluation and normalised dose difference algorithms [2] were employed for the comparison of the TPS dose distribution and the MC dose distribution. These implementations are spatial resolution independent and able to interpolate for comparisons. Results and Discussion: The tools successfully produced Monte Carlo input files for a variety of plans exported from the Eclipse (Varian Medical Systems) and Pinnacle (Philips Medical Systems) planning systems: ranging in complexity from a single uniform square field to a five-field step and shoot IMRT treatment. The simulation of collimated beams has been verified geometrically, and validation of dose distributions in a simple body phantom (QUASAR) will follow. The developed dose comparison algorithms have also been tested with controlled dose distribution changes. Conclusion: The capability of the developed code to independently process treatment plans has been demonstrated. A number of limitations exist: only static fields are currently supported (dynamic wedges and dynamic IMRT will require further development), and the process has not been tested for planning systems other than Eclipse and Pinnacle. The tools will be used to independently assess the accuracy of the current treatment planning system dose calculation algorithms for complex treatment deliveries such as IMRT in treatment sites where patient inhomogeneities are expected to be significant. Acknowledgements: Computational resources and services used in this work were provided by the HPC and Research Support Group, Queensland University of Technology, Brisbane, Australia. Pinnacle dose parsing made possible with the help of Paul Reich, North Coast Cancer Institute, North Coast, New South Wales.

Pixel value to electron density conversion for a Megavoltage Cone Beam CT System

Introduction: The motivation for developing megavoltage (and kilovoltage) cone beam CT (MV CBCT) capabilities in the radiotherapy treatment room was primarily based on the need to improve patient set-up accuracy. There has recently been an interest in using the cone beam CT data for treatment planning. Accurate treatment planning, however, requires knowledge of the electron density of the tissues receiving radiation in order to calculate dose distributions. This is obtained from CT, utilising a conversion between CT number and electron density of various tissues. The use of MV CBCT has particular advantages compared to treatment planning with kilovoltage CT in the presence of high atomic number materials and requires the conversion of pixel values from the image sets to electron density. Therefore, a study was undertaken to characterise the pixel value to electron density relationship for the Siemens MV CBCT system, MVision, and determine the effect, if any, of differing the number of monitor units used for acquisition. If a significant difference with number of monitor units was seen then pixel value to ED conversions may be required for each of the clinical settings. The calibration of the MV CT images for electron density offers the possibility for a daily recalculation of the dose distribution and the introduction of new adaptive radiotherapy treatment strategies. Methods: A Gammex Electron Density CT Phantom was imaged with the MVCB CT system. The pixel value for each of the sixteen inserts, which ranged from 0.292 to 1.707 relative electron density to the background solid water, was determined by taking the mean value from within a region of interest centred on the insert, over 5 slices within the centre of the phantom. These results were averaged and plotted against the relative electron densities of each insert with a linear least squares fit was preformed. This procedure was performed for images acquired with 5, 8, 15 and 60 monitor units. Results: The linear relationship between MVCT pixel value and ED was demonstrated for all monitor unit settings and over a range of electron densities. The number of monitor units utilised was found to have no significant impact on this relationship. Discussion: It was found that the number of MU utilised does not significantly alter the pixel value obtained for different ED materials. However, to ensure the most accurate and reproducible MV to ED calibration, one MU setting should be chosen and used routinely. To ensure accuracy for the clinical situation this MU setting should correspond to that which is used clinically. If more than one MU setting is used clinically then an average of the CT values acquired with different numbers of MU could be utilized without loss in accuracy. Conclusions: No significant differences have been shown between the pixel value to ED conversion for the Siemens MV CT cone beam unit with change in monitor units. Thus as single conversion curve could be utilised for MV CT treatment planning. To fully utilise MV CT imaging for radiotherapy treatment planning further work will be undertaken to ensure all corrections have been made and dose calculations verified. These dose calculations may be either for treatment planning purposes or for reconstructing the delivered dose distribution from transit dosimetry measurements made using electronic portal imaging devices. This will potentially allow the cumulative dose distribution to be determined through the patient’s multi-fraction treatment and adaptive treatment strategies developed to optimize the tumour response.

Correcting radiation survey data to account for increased leakage during intensity modulated radiotherapy treatments

Purpose Intensity modulated radiotherapy (IMRT) treatments require more beam-on time and produce more linac head leakage to deliver similar doses to conventional, unmodulated, radiotherapy treatments. It is necessary to take this increased leakage into account when evaluating the results of radiation surveys around bunkers that are, or will be, used for IMRT. The recommended procedure of 15 applying a monitor-unit based workload correction factor to secondary barrier survey measurements, to account for this increased leakage when evaluating radiation survey measurements around IMRT bunkers, can lead to potentially-costly over estimation of the required barrier thickness. This study aims to provide initial guidance on the validity of reducing the value of the correction factor when applied to different radiation barriers (primary barriers, doors, maze walls and other walls) by 20 evaluating three different bunker designs. Methods Radiation survey measurements of primary, scattered and leakage radiation were obtained at each of five survey points around each of three different radiotherapy bunkers and the contribution of leakage to the total measured radiation dose at each point was evaluated. Measurements at each survey point were made with the linac gantry set to 12 equidistant positions from 0 to 330o, to 25 assess the effects of radiation beam direction on the results. Results For all three bunker designs, less than 0.5% of dose measured at and alongside the primary barriers, less than 25% of the dose measured outside the bunker doors and up to 100% of the dose measured outside other secondary barriers was found to be caused by linac head leakage. Conclusions Results of this study suggest that IMRT workload corrections are unnecessary, for 30 survey measurements made at and alongside primary barriers. Use of reduced IMRT workload correction factors is recommended when evaluating survey measurements around a bunker door, provided that a subset of the measurements used in this study are repeated for the bunker in question. Reduction of the correction factor for other secondary barrier survey measurements is not recommended unless the contribution from leakage is separetely evaluated.

Prostate radiotherapy treatment plan quality : effects of hip prostheses

Introduction This study aimed to examine the geometric and dosimetric results when radiotherapy treatment plans are designed for prostate cancer patients with hip prostheses. Methods Ten EBRT treatment plans for localised prostate cancer, in the presence of hip prostheses, were analysed and compared with a reference set of 196 treatment plans for localised prostate cancer in patients without prostheses. Crowe et al.’s TADA code [1] was used to extract treatment plan parameters and evaluate doses to target volumes and critical structures against recommended goals [2] and constraints [3, 4]. Results The need to avoid transmitting the radiation beam through the hip prostheses limited the range of gantry angles available for use in both the rotational (VMAT) and the non-rotational (3DCRT and IMRT) radiotherapy treatments. This geometric limitation (exemplified in the VMAT data shown in Fig. 1) reduced the overall quality of the treatment plans for patients with prostheses compared to the reference plans. All plans with prostheses failed the PTV dose homogeneity requirement [2], whereas only 4 % of the plans without prostheses failed this test. Several treatment plans for patients with hip prostheses also failed the QUANTEC requirements that less than 50 % of the rectum receive 50 Gy and less than 35 % of the rectum receive 60 Gy to keep the grade 3 toxicity rate below 10 % [3], or the Hansen and Roach requirement that less than 25 % of the bladder receive 75 Gy [4]. Discussion and conclusions The results of this study exemplify the difficulty of designing prostate radiotherapy treatment plans, where beams provide adequate doses to targeted tissues while avoiding nearby organs at risk, when the presence of hip prostheses limits the available treatment geometries. This work provides qualitative evidence of the compromised dose distributions that can result, in such cases.

Radiosensitizing effect of electrochemotherapy in a fractionated radiation regimen in radiosensitive murine sarcoma and radioresistant adenocarcinoma tumor model

Radiosensitizing Effect of Electrochemotherapy in a Fractionated Radiation Regimen in Radiosensitive Murine Sarcoma and Radioresistant Adenocarcinoma Tumor Model. Electrochemotherapy can potentiate the radiosensitizing effect of bleomycin, as shown in our previous studies. To bring this treatment closer to use in clinical practice, we evaluated the interaction between electrochemotherapy with bleomycin and single-dose or fractionated radiation in two murine tumor models with different histology and radiosensitivity. Radiosensitive sarcoma SA-1 and radioresistant adenocarcinoma CaNT subcutaneous tumors grown in A/J and CBA mice, respectively, were used. The anti-tumor effect and skin damage around the treated tumors were evaluated after electrochemotherapy with bleomycin alone or combined with single-dose radiation or a fractionated radiation regimen. The anti-tumor effectiveness of electrochemotherapy was more pronounced in SA-1 than CaNT tumors. In both tumor models, the tumor response to radiation was not significantly influenced by bleomycin alone or by electroporation alone. However, electrochemotherapy before the first tumor irradiation potentiated the response to a single-dose or fractionated radiation regimen in both tumors. For the fractionated radiation regimen, normal skin around the treated tumors was damaged fourfold less than for the single-dose regimen. Electrochemotherapy prior to single-dose irradiation induced more damage to the skin around the treated tumors and greater loss of body weight compared to other irradiated groups, whereas electrochemotherapy combined with the fractionated radiation regimen did not. Electrochemotherapy with low doses of bleomycin can also be used safely for radiosensitization of different types of tumors in a fractionated radiation regimen, resulting in a good anti-tumor effect and no major potentiating effect on radiation-induced skin damage. © 2009 by Radiation Research Society.

Interplay of MLC, gantry and respiratory motion during DCAT delivery

This study investigated the possible interplay effects arising from the treatment of moving targets using the dynamic conformal arc therapy (DCAT) technique. Dose from a modulated test beam was measured, with and without phantom motion and with and without a 30o arc rotation, using a diode array placed on a sinusoidally moving platform. Measurements were repeated at five different collimator angles (0, 22.5, 45, 67.5 and 90o), at two different dose rates (300 and 600 MU/min). Results showed that the effect of respiratory motion on the measured dose distribution increased slightly when the beams were delivered as arcs, rather than with a static gantry angle, and that this effect increased substantially as the collimator angle was increased from 0o (MLC motion perpendicular to respiratory motion) to 90o (MLC motion parallel to respiratory motion). The dose oscillations arising from interplay between phantom and MLC motion were found to increase in magnitude when the dose rate was increased. These results led to the development of simple recommendations for minimizing the negative effects of motion interplay on DCAT dose distributions

Investigation of the relationship between BMI and prostate motion in radiation therapy treatment for prostate cancer

This thesis studied the influence of patient obesity on prostate motion during radiation therapy treatment delivery, an important consideration in the accurate treatment of prostate cancer. The study highlighted the importance of daily image guidance to correct for prostate motion, increasing radiation dose to the prostate while decreasing radiation dose to surrounding healthy tissues, thereby increasing patient quality of life.

Comparison of Radiation Doses for Catheter Ablation Procedures for Atrial Fibrillation with Percutaneous Coronary Intervention Cases

Background: Catheter ablation procedures for atrial fibrillation (AF) may frequently require long fluoroscopic times. We sought to undertake a review of radiation safety practice in our Cardiac Electrophysiology Laboratory and implement changes to minimize fluoroscopic doses. We also sought to compare the results with radiation doses for percutaneous coronary intervention (PCI) cases performed in our hospital. Methods: Fluoroscopic times and doses for AF ablation procedures performed by a single operator on a Philips Integris H3000 image-intensifier were analysed for 11-month period. Results were compared with all PCI procedures performed over a similar period by multiple operators on a Philips Integris Allura FD system. Comprehensive review of radiation practice in the Electrophysiology laboratory identified the potential to reduce pulse frame rates and doses, and to narrow the field of interest without impacting the performance of the procedure. These changes were implemented and results analysed after a further 11 months. Results: In the pre-intervention period 50 AF catheter ablations had a mean fluoroscopic time of 86.4 min and mean fluoroscopic dose 68.4 Gy/cm2. Post-intervention 75 procedures had a mean fluorosocopic time of 68.9 min (p < 0.0001) and mean dose of 14.3 Gy/cm2 (p < 0.0001) 128 PCI procedures had a mean combined fluoroscopic and image acquisition time of 10.0 min and mean total dose 38.8 Gy/cm2. Conclusions: Catheter ablation procedures for AF may require lengthy use of fluoroscopy but simple modifications to radiation practice can result in marked reductions in radiation dose that compare favourably with PCI case doses

On formation, growth and concentrations of air ions

Aerosol particles have effect on climate, visibility, air quality and human health. However, the strength of which aerosol particles affect our everyday life is not well described or entirely understood. Therefore, investigations of different processes and phenomena including e.g. primary particle sources, initial steps of secondary particle formation and growth, significance of charged particles in particle formation, as well as redistribution mechanisms in the atmosphere are required. In this work sources, sinks and concentrations of air ions (charged molecules, cluster and particles) were investigated directly by measuring air molecule ionising components (i.e. radon activity concentrations and external radiation dose rates) and charged particle size distributions, as well as based on literature review. The obtained results gave comprehensive and valuable picture of the spatial and temporal variation of the air ion sources, sinks and concentrations to use as input parameters in local and global scale climate models. Newly developed air ion spectrometers (Airel Ltd.) offered a possibility to investigate atmospheric (charged) particle formation and growth at sub-3 nm sizes. Therefore, new visual classification schemes for charged particle formation events were developed, and a newly developed particle growth rate method was tested with over one year dataset. These data analysis methods have been widely utilised by other researchers since introducing them. This thesis resulted interesting characteristics of atmospheric particle formation and growth: e.g. particle growth may sometimes be suppressed before detection limit (~ 3 nm) of traditional aerosol instruments, particle formation may take place during daytime as well as in the evening, growth rates of sub-3 nm particles were quite constant throughout the year while growth rates of larger particles (3-20 nm in diameter) were higher during summer compared to winter. These observations were thought to be a consequence of availability of condensing vapours. The observations of this thesis offered new understanding of the particle formation in the atmosphere. However, the role of ions in particle formation, which is not well understood with current knowledge, requires further research in future.

Dose rate dependence of electrical characteristics of lead zirconate titanate capacitors

The influence of gamma-radiation dose rate on the electrical properties of lead zirconate titanate capacitors was investigated. More severe degradations in dielectric constant, coercive field, remanent polarization and capacitance-voltage (C-V) curves occurred with increasing radiation dose at lower dose rates. The electrical properties exhibited distinct radiation dose rate dependence and the worst-case degradation occurred at the lowest dose rate. The radiation-induced degradation of parameters such as the coercive field drift and distortion of the C-V curve can be recovered partly through post-irradiation annealing.

Simulations to design an online motion compensation system for scanned particle beams

Respiration-induced target motion is a major problem in intensity-modulated radiation therapy. Beam segments are delivered serially to form the total dose distribution. In the presence of motion, the spatial relation between dose deposition from different segments will be lost. Usually, this results in over-and underdosage. Besides such interplay effects between target motion and dynamic beam delivery as known from photon therapy, changes in internal density have an impact on delivered dose for intensity-modulated charged particle therapy. In this study, we have analysed interplay effects between raster scanned carbon ion beams and target motion. Furthermore, the potential of an online motion strategy was assessed in several simulations. An extended version of the clinical treatment planning software was used to calculate dose distributions to moving targets with and without motion compensation. For motion compensation, each individual ion pencil beam tracked the planned target position in the lateral aswell as longitudinal direction. Target translations and rotations, including changes in internal density, were simulated. Target motion simulating breathing resulted in severe degradation of delivered dose distributions. For example, for motion amplitudes of +/- 15 mm, only 47% of the target volume received 80% of the planned dose. Unpredictability of resulting dose distributions was demonstrated by varying motion parameters. On the other hand, motion compensation allowed for dose distributions for moving targets comparable to those for static targets. Even limited compensation precision (standard deviation similar to 2 mm), introduced to simulate possible limitations of real-time target tracking, resulted in less than 3% loss in dose homogeneity.

Effect of gamma-ray radiation on physiological, morphological characters and chromosome aberrations of minitubers in Solanum tuberosum L.

Purpose: To investigate the effects of gamma-ray radiation on the physiological, morphological characters and chromosome aberrations of minitubers. Materials and methods: Minitubers of one potato cultivar, 'Shepody', were irradiated with 8 doses of gamma-rays (0, 10, 20, 30, 40, 50, 60, 70 and 80 Gy [Gray]) to investigate the effects of radiation on emergence ability, plant height and root length, morphological variations, chromosome aberrations, M-1 (first generation mutants) tuber number and size of minituber plants. Results: Compared with the non-irradiated controls, the whole period of emergence was prolonged by 10-15 days for minitubers treated with gamma-ray radiation, but low doses of radiation (10, 20 and 30 Gy) promoted the emergence percentage of minitubers. With an increase in radiation dose, the emergence percentage, plant height and root length of minituber plants were significantly inhibited at 40 and 50 Gy. No emergence occurred at 60 Gy and higher doses. After radiation, a series of morphological variations and chromosome aberrations appeared in minituber plants. Radiation with 20 Gy promoted tuber formation, and the average number and diameter of M-1 tubers per plant were significantly increased over the control by 71% and 34%, respectively. Conclusion: Low doses of radiation (10-30 Gy) might be used as a valuable parameter to study the improvement of minitubers by gamma-ray radiation treatment.

Effect of gamma-radiation on the thermal and mechanical properties of a commercial poly(butylene adipate-co-terephthalate)

BACKGROUND: Poly(butylene adipate-co-terephthalate) (PBAT) has attracted wide interest as a biodegradable polymer. However, its use is restricted in certain applications due to its low melting point.RESULTS: PBAT was treated using gamma-radiation. The radiation features were analyzed using Soxhlet extraction, and the ratio of chain scission and crosslinking and gelation dose were determined using the classical Charlesby-Pinner equation. The results showed that PBAT is a radiation-crosslinkable polymer. The degree of crosslinking increased with increasing radiation dose; the relation between sol fraction and dose followed the Charlesby-Pinner equation. Differential scanning calorimetry analyses showed that the melting temperature (T-m) and the heat of fusion (Delta H-m) of PBAT exhibited almost no change in the first scan. The second scan, however, showed a decrease in T-m and Delta H-m. The glass transition temperature of irradiated PBAT increased with increasing radiation dose.

LDPE/carbon black conductive composites: Influence of radiation crosslinking on PTC and NTC properties

The electrical resistivity of low-density polyethylene/carbon black composites irradiated by Co-60 gamma-rays was investigated as a function of temperature. The experimental results obtained by scanning electron microscopy, solvent extraction techniques, and pressure-specific volume-temperature analysis techniques showed that the positive temperature coefficient (PTC) and negative temperature coefficient (NTC) effects of the composites were influenced by the irradiation dose, network forming (gel), and soluble fractions (Sol). The NTC effect was effectively eliminated when the radiation dose reached 400 kGy. The results showed that the elimination of the NTC effect was related to the difference in the thermal expansion of the gel and Sol regions. The thermal expansion of the sol played an important role in both increasing the PTC intensity and decreasing the NTC intensity at 400 kGy.