Total skin electron irradiation (TSEI) is a special radiotherapy technique which has generally been used for treating adult patients with mycosis fungoides. Recently, two infants presented with leukemia cutis isolated to the skin requiring TSEI. This work discusses the commissioning and quality assurance (QA) methods for implementing a modified Stanford technique using a rotating harness system to position sedated pediatric patients treated with electrons to the total skin.
Methods and Results
Commissioning of pediatric TSEI consisted of absolute calibration, measurement of dosimetric parameters, and subsequent verification in a pediatric patient sized cylindrical phantom using radiographic film and optically stimulated luminance (OSL) dosimeters. The depth of dose penetration under TSEI treatment condition was evaluated using radiographic film sandwiched in the phantom and demonstrated a 2 cm penetration depth with the maximum dose located at the phantom surface. Dosimetry measurements on the cylindrical phantom and in-vivo measurements from the patients suggested that, the factor relating the skin and calibration point doses (i.e., the B-factor) was larger for the pediatric TSEI treatments as compared to adult TSEI treatments. Custom made equipment, including a rotating plate and harness, was fabricated and added to a standard total body irradiation stand and tested to facilitate patient setup under sedated condition. A pediatric TSEI QA program, consisting of daily output, energy, flatness, and symmetry measurements as well as in-vivo dosimetry verification for the first cycle was developed. With a long interval between pediatric TSEI cases, absolute dosimetry was also repeated as part of the QA program. In-vivo dosimetry for the first two infants showed that a dose of ± 10% of the prescription dose can be achieved over the entire patient body.
Though pediatric leukemia cutis and the subsequent need for TSEI are rare, the ability to commission the technique on a modified TBI stand is appealing for clinical implementation and has been successfully used for the treatment of two pediatric patients at our institution.
Pediatric total skin electron irradiation; Commissioning; Quality assurance; Leukemia cutis
An anthropomorphic phantom was used to investigate a treatment technique and analyze the dose distributions for helical irradiation of the total skin (HITS) by helical tomotherapy (HT). Hypothetical bolus of thicknesses of 0, 10, and 15 mm was added around the phantom body to account for the dose homogeneity and setup uncertainty. A central core structure was assigned as a “complete block” to force the dose tangential delivery. HITS technique with prescribed dose (Dp) of 36 Gy in 36 fractions was generated. The radiochromic EBT2 films were used for the dose measurements. The target region with 95.0% of the Dp received by more than 95% of the PTV was obtained. The calculated mean doses for the organs at risk (OARs) were 4.69, 3.10, 3.20, and 2.94 Gy for the lung, heart, liver, and kidneys, respectively. The measurement doses on a phantom surface for a plan with 10 mm hypothetical bolus and bolus thicknesses of 0, 1, 2, and 3 mm are 89.5%, 111.4%, 116.9%, and 117.7% of Dp, respectively. HITS can provide an accurate and uniform treatment dose in the skin with limited doses to OARs and is safe to replace a total skin electron beam regimen.
Computational phantoms representing workers and patients are essential in estimating organ doses from various occupational radiation exposures and medical procedures. Nearly all existing phantoms, however, were purposely designed to match internal and external anatomical features of the Reference Man as defined by the International Commission on Radiological Protection (ICRP). To reduce uncertainty in dose calculations caused by anatomical variations, a new generation of phantoms of varying organ and body sizes is needed. This paper presents detailed anatomical data in tables and graphs that are used to design such size-adjustable phantoms representing a range of adult individuals in terms of the body height, body weight and internal organ volume/mass. Two different sets of information are used to derive the phantom sets: (1) individual internal organ size and volume/mass distribution data derived from the recommendations of the ICRP in Publications 23 and 89 and (2) whole-body height and weight percentile data from the National Health and Nutrition Examination Survey (NHANES 1999–2002). The NHANES height and weight data for 19 year old males and females are used to estimate the distributions of individuals’ size, which is unknown, that corresponds to the ICRP organ and tissue distributions. This paper then demonstrates the usage of these anthropometric data in the development of deformable anatomical phantoms. A pair of phantoms—modeled entirely in mesh surfaces—of the adult male and female, RPI-adult male (AM) and RPI-adult female (AF) are used as the base for size-adjustable phantoms. To create percentile-specific phantoms from these two base phantoms, organ surface boundaries are carefully altered according to the tabulated anthropometric data. Software algorithms are developed to automatically match the organ volumes and masses with desired values. Finally, these mesh-based, percentile-specific phantoms are converted into voxel-based phantoms for Monte Carlo radiation transport simulations. This paper also compares absorbed organ doses for the RPI-AM-5th-height and -weight percentile phantom (165 cm in height and 56 kg in weight) and the RPI-AM-95th-height and -weight percentile phantom (188 cm in height and 110 kg in weight)with those for theRPI-AM-50th-height and -weight percentile phantom (176 cm in height and 73 kg in weight) from exposures to 0.5 MeV external photon beams. The results suggest a general finding that the phantoms representing a slimmer and shorter individual male received higher absorbed organ doses because of lesser degree of photon attenuation due to smaller amount of body fat. In particular, doses to the prostate and adrenal in the RPI-AM-5th-height and -weight percentile phantom is about 10% greater than those in the RPI-AM-50th-height and -weight percentile phantom approximating the ICRP Reference Man. On the other hand, the doses to the prostate and adrenal in the RPI-AM-95th-height and -weight percentile phantom are approximately 20% greater than those in the RPI-AM-50th-height and -weight percentile phantom. Although this study only considered the photon radiation of limited energies and irradiation geometries, the potential to improve the organ dose accuracy using the deformable phantom technology is clearly demonstrated.
The potential of the PRESAGE™/Optical-CT system as a comprehensive 3D dosimetry tool has been demonstrated. The current study focused on detailed characterization of robustness (intra-dosimeter uniformity and temporal stability) and reproducibility (inter-dosimeter reproducibility) of PRESAGE™ inserts compatible with the RPC H&N phantom. In addition, the accuracy and precision of PRESAGE dose measurement was also evaluated. Four identical PRESAGE™ dosimeters (10cm diameter and 7cm height cylinders) were irradiated with the same rotationally symmetric treatment plan using a Varian accelerator. The treatment plan was designed to rigorously evaluate robustness and reproducibility for multiple dose levels and in 3D. All dosimeters were scanned by optical-CT at daily intervals to study temporal stability. Dose comparisons were made between PRESAGE, ECLIPSE, and independent measurement with EBT film at a select depth. The use of improved optics and acquisition technique yielded substantially higher quality 3D dosimetry data from PRESAGE than has been achieved previously (noise reduced to ~1%, accuracy to within 3%). Data analysis showed excellent intra-dosimeter uniformity, temporal stability and inter-dosimeter reproducibility of relative radiochromic response. In general, the PRESAGE™ dose-distribution was found to agree better with EBT (~99% pass rate) than with ECLIPSE calculations (~92% pass rate) especially in penumbral regions for a 3% dose-difference and 3 mm distance-to-agreement evaluation criteria. The results demonstrate excellent robustness and reproducibility of the PRESAGE™ for relative 3D-dosimetry and represent a significant step towards incorporation in the RadOnc-clinic (e.g. integration with RPC phantom).
A general method is presented for patient-specific 3-dimensional absorbed dose calculations based on quantitative SPECT activity measurements.
The computational scheme includes a method for registration of the CT image to the SPECT image and position-dependent compensation for attenuation, scatter, and collimator detector response performed as part of an iterative reconstruction method. A method for conversion of the measured activity distribution to a 3-dimensional absorbed dose distribution, based on the EGS4 (electron-gamma shower, version 4) Monte Carlo code, is also included. The accuracy of the activity quantification and the absorbed dose calculation is evaluated on the basis of realistic Monte Carlo–simulated SPECT data, using the SIMIND (simulation of imaging nuclear detectors) program and a voxel-based computer phantom. CT images are obtained from the computer phantom, and realistic patient movements are added relative to the SPECT image. The SPECT-based activity concentration and absorbed dose distributions are compared with the true ones.
Correction could be made for object scatter, photon attenuation, and scatter penetration in the collimator. However, inaccuracies were imposed by the limited spatial resolution of the SPECT system, for which the collimator response correction did not fully compensate.
The presented method includes compensation for most parameters degrading the quantitative image information. The compensation methods are based on physical models and therefore are generally applicable to other radionuclides. The proposed evaluation methodology may be used as a basis for future intercomparison of different methods.
Monte Carlo; absorbed dose; SPECT; 131I; image registration
To analyse limits and capabilities in dose calculation of collapsed-cone-convolution (CCC) algorithm implemented in helical tomotherapy (HT) treatment planning system for thorax lesions.
The agreement between measured and calculated dose was verified both in homogeneous (Cheese Phantom) and in a custom-made inhomogeneous phantom. The inhomogeneous phantom was employed to mimic a patient's thorax region with lung density encountered in extreme cases and acrylic inserts of various dimensions and positions inside the lung cavity. For both phantoms, different lung treatment plans (single or multiple metastases and targets in the mediastinum) using HT technique were simulated and verified. Point and planar dose measurements, both with radiographic extended-dose-range (EDR2) and radiochromic external-beam-therapy (EBT2) films, were performed. Absolute point dose measurements, dose profile comparisons and quantitative analysis of gamma function distributions were analyzed.
An excellent agreement between measured and calculated dose distributions was found in homogeneous media, both for point and planar dose measurements. Absolute dose deviations <3% were found for all considered measurement points, both inside the PTV and in critical structures. Very good results were also found for planar dose distribution comparisons, where at least 96% of all points satisfied the gamma acceptance criteria (3%-3 mm), both for EDR2 and for EBT2 films. Acceptable results were also reported for the inhomogeneous phantom. Similar point dose deviations were found with slightly worse agreement for the planar dose distribution comparison: 96% of all points passed the gamma analysis test with acceptable levels of 4%-4 mm and 5%-4 mm, for EDR2 and EBT2 films respectively. Lower accuracy was observed in high dose/low density regions, where CCC seems to overestimate the measured dose around 4-5%.
Very acceptable accuracy was found for complex lung treatment plans calculated with CCC algorithm implemented in the tomotherapy TPS even in the heterogeneous phantom with very low lung-density.
Patient-specific dose verification for treatment planning in helical tomotherapy is routinely performed using a homogeneous virtual water cylindrical phantom of 30 cm diameter and 18 cm length (Cheese phantom). Because of this small length, treatment with total marrow irradiation (TMI) requires multiple deliveries of the dose verification procedures to cover a wide range of the target volumes, which significantly prolongs the dose verification process. We propose a fast, simple, and informative patient-specific dose verification method which reduce dose verification time for TMI with helical tomotherapy.
We constructed a two-step solid water slab phantom (length 110 cm, height 8 cm, and two-step width of 30 cm and 15 cm), termed the Whole Body Phantom (WB phantom). Three ionization chambers and three EDR-2 films can be inserted to cover extended field TMI treatment delivery. Three TMI treatment plans were conducted with a TomoTherapy HiArt Planning Station and verified using the WB phantom with ion chambers and films. Three regions simulating the head and neck, thorax, and pelvis were covered in a single treatment delivery. The results were compared to those with the cheese phantom supplied by Accuray, Inc. following three treatment deliveries to cover the body from head to pelvis.
Use of the WB phantom provided point doses or dose distributions from head and neck to femur in a single treatment delivery of TMI. Patient-specific dose verification with the WB phantom was 62% faster than with the cheese phantom. The average pass rate in gamma analysis with the criteria of a 3-mm distance-to-agreement and 3% dose differences was 94% ± 2% for the three TMI treatment plans. The differences in pass rates between the WB and cheese phantoms at the upper thorax to abdomen regions were within 2%. The calculated dose agreed with the measured dose within 3% for all points in all five cases in both the WB and cheese phantoms.
Our dose verification method with the WB phantom provides simple and rapid quality assurance without limiting dose verification information in total marrow irradiation with helical tomotherapy.
Glioblastoma multiforme (GBM) is the most common, aggressive, highly malignant and infiltrative of all brain tumors with low rate of control. The main goal of this work was to evaluate the spatial dose distribution into a GBM simulator inside a head phantom exposed to a 15 MV 3D conformal radiation therapy in order to validate internal doses. A head and neck phantom developed by the Ionizing Radiation Research Group (NRI) was used on the experiments. Such phantom holds the following synthetic structures: brain and spinal cord, skull, cervical and thoracic vertebrae, jaw, hyoid bone, laryngeal cartilages, head and neck muscles and skin. Computer tomography (CT) of the simulator was taken, capturing a set of contrasted references. Therapy Radiation planning (TPS) was performed based on those CT images, satisfying a 200 cGy prescribed dose split in three irradiation fields. The TPS assumed 97% of prescribed dose cover the prescribed treatment volume (PTV). Radiochromic films in a solid water phantom provided dose response as a function of optical density. Spatial dosimetric distribution was generated by radiochromic film samples at coronal, sagittal-anterior and sagittal-posterior positions, inserted into tumor simulator and brain. The spatial dose profiles held 70 to 120% of the prescribed dose. In spite of the stratified profile, as opposed to the smooth dose profile from TPS, the tumor internal doses were within a 5% deviation from 214.4 cGy evaluated by TPS. 83.2% of the points with a gamma value of less than 1 (3%/3mm) for TPS and experimental values, respectively. At the tumor, measured at coronal section, a few dark spots in the film caused the appearance of outlier points in 13-15% of dose deviation percentage. And, as final conclusion, such dosimeter choice and the physical anthropomorphic and anthropometric phantom provided an efficient method for validating radiotherapy protocols.
Space radiation is known to be key hazard of manned space mission. To estimate accurately radiation health risk detailed study of dose distribution inside human body by means of human phantom is conducted. In the space experiment MATROSHKA-R, the tissue-equivalent spherical phantom (32 kg mass, 35 cm diameter and 10 cm central spherical cave) made in Russia has been used on board the ISS for more than 8 years. Owing to the specially chosen phantom shape and size, the chord length distributions of the detector locations are attributed to self-shielding properties of the critical organs in a real human body. If compared with the anthropomorphic phantom Rando used inside and outside the ISS, the spherical phantom has lower mass, smaller size and requires less crew time for the detector installation/retrieval; its tissue-equivalent properties are closer to the standard human body tissue than the Rando-phantom material. Originally the spherical phantom was installed in the star board crew cabin of the ISS Service Module, then in the Piers-1, MIM-2 and MIM-1 modules of the ISS Russian segment, and finally in JAXA Kibo module. Total duration of the detector exposure is more than 1700 days in 8 sessions.
In the first phase of the experiment with the spherical phantom, the dose measurements were realized with only passive detectors (thermoluminescent and solid-state track detectors). The detectors are placed inside the phantom along the axes of 20 containers and on the phantom outer surface in 32 pockets of the phantom jacket. After each session the passive detectors are returned to the ground. The results obtained show the dose difference on the phantom surface as much as a factor of 2, the highest dose being observed close to the outer wall of the compartment, and the lowest dose being in the opposite location along the phantom diameter. Maximum dose rate measured in the phantom is obviously due to the galactic cosmic ray (GCR) and Earth' radiation belt contribution on the ISS trajectory. Minimum dose rate is caused mainly by the strongly penetrating GCR particles and is observed behind more than 5 g/cm2 tissue shielding. Critical organ doses, mean-tissue and effective doses of a crew member in the ISS compartments are also estimated with the spherical phantom data. The estimated effective dose rate is found to be from 10 to 15% lower than the averaged dose on the phantom surface as dependent on the attitude of the critical organs.
The spherical phantom proved its effectiveness to measure the critical organ doses together with the effective dose in-flight and if supplied with active dosimeters can be recommended for future exploratory manned missions to monitor continuously the effective dose.
space radiation; spherical phantom; passive and active detectors; effective dose
Estimating organ residence times is an essential part of patient-specific dosimetry for radioimmunotherapy (RIT). Quantitative imaging methods for RIT are often evaluated using a single physical or simulated phantom but are intended to be applied clinically where there is variability in patient anatomy, biodistribution, and biokinetics. To provide a more relevant evaluation, the authors have thus developed a population of phantoms with realistic variations in these factors and applied it to the evaluation of quantitative imaging methods both to find the best method and to demonstrate the effects of these variations. Using whole body scans and SPECT/CT images, organ shapes and time-activity curves of 111In ibritumomab tiuxetan were measured in dosimetrically important organs in seven patients undergoing a high dose therapy regimen. Based on these measurements, we created a 3D NURBS-based cardiac-torso (NCAT)-based phantom population. SPECT and planar data at realistic count levels were then simulated using previously validated Monte Carlo simulation tools. The projections from the population were used to evaluate the accuracy and variation in accuracy of residence time estimation methods that used a time series of SPECT and planar scans. Quantitative SPECT (QSPECT) reconstruction methods were used that compensated for attenuation, scatter, and the collimator-detector response. Planar images were processed with a conventional (CPlanar) method that used geometric mean attenuation and triple-energy window scatter compensation and a quantitative planar (QPlanar) processing method that used model-based compensation for image degrading effects. Residence times were estimated from activity estimates made at each of five time points. The authors also evaluated hybrid methods that used CPlanar or QPlanar time-activity curves rescaled to the activity estimated from a single QSPECT image. The methods were evaluated in terms of mean relative error and standard deviation of the relative error in the residence time estimates taken over the phantom population. The mean errors in the residence time estimates over all the organs were <9.9% (pure QSPECT), <13.2% (pure QPLanar), <7.2% (hybrid QPlanar/QSPECT), <19.2% (hybrid CPlanar/QSPECT), and 7%–159% (pure CPlanar). The standard deviations of the errors for all the organs over all the phantoms were <9.9%, <11.9%, <10.8%, <22.0%, and <107.9% for the same methods, respectively. The processing methods differed both in terms of their average accuracy and the variation of the accuracy over the population of phantoms, thus demonstrating the importance of using a phantom population in evaluating quantitative imaging methods. Hybrid CPlanar/QSPECT provided improved accuracy compared to pure CPlanar and required the addition of only a single SPECT acquisition. The QPlanar or hybrid QPlanar/QSPECT methods had mean errors and standard deviations of errors that approached those of pure QSPECT while providing simplified image acquisition protocols, and thus may be more clinically practical.
absolute quantitation; phantom population; quantitative SPECT; radioimmunotherapy
A thermobrachytherapy surface applicator (TBSA) was developed for simultaneous heat and brachytherapy treatment of chestwall (CW) recurrence of breast cancer. The ability to comfortably secure the applicator over the upper torso relative to the CW target throughout treatment is assessed on volunteers.
Male and post-mastectomy female volunteers were enrolled to evaluate applicator secure fit to CW. Female subjects with intact breast were also enrolled to assess the ability to treat challenging cases. Magnetic resonance (MR) images of volunteers wearing a TBSA over the upper torso were acquired once every 15 minutes for 90 minutes. Applicator displacement over this time period required for treatment preplanning and delivery was assessed using MR visible markers. Applicator comfort and tolerability were assessed using a questionnaire. Probability estimates of applicator displacements were used to investigate dosimetric impact for the worst case variation in radiation source-to-skin distance for 5 and 10 mm deep targets spread 17×13cm on a torso phantom. Average and median displacements along lateral and radial directions were less than 1.2 mm over 90 minutes for all volunteers. Maximum lateral and radial displacements were measured to be less than 1 and 1.5 mm respectively for all CW volunteers and, less than 2 mm for intact breast volunteers excluding outliers. No complaint of pain or discomfort was reported. Phantom treatment planning for the maximum displacement of 2mm indicated < 10% increase in skin dose with < 5% loss of homogeneity index (HI) for −2 mm uniform HDR source displacement. For +2 mm uniform displacement, skin dose decreased and HI increased by 20%. The volunteer study demonstrated that such large and uniform displacements should be rare for CW subjects and the measured variation is expected to be low for multi-fraction conformal brachytherapy treatment.
Hyperthermia; conformal microwave array; brachytherapy; chestwall recurrence; breast cancer
To develop a technique for the construction of a two-compartment anthropomorphic breast phantom specific to an individual patient’s pendant breast anatomy.
Three-dimensional breast images were acquired on a prototype dedicated breast computed tomography (bCT) scanner as part of an ongoing IRB-approved clinical trial of bCT. The images from the breast of a patient were segmented into adipose and glandular tissue regions and divided into 1.59 mm thick breast sections to correspond to the thickness of polyethylene stock. A computer controlled water-jet cutting machine was used to cut the outer breast edge and the internal regions corresponding to glandular tissue from the polyethylene. The stack of polyethylene breast segments was encased in a thermoplastic “skin” and filled with water. Water-filled spaces modeled glandular tissue structures and the surrounding polyethylene modeled the adipose tissue compartment. Utility of the phantom was demonstrated by inserting 200 μm microcalcifications as well as measuring point dose deposition during bCT scanning.
Rigid registration of the original patient images with bCT images of the phantom showed similar tissue distribution. Linear profiles through the registered images demonstrated a mean coefficient of determination (r2) between grayscale profiles of 0.881. The exponent of the power law describing the anatomical noise power spectrum was identical in the coronal images of the patient’s breast and the phantom. Microcalcifications were visualized in the phantom at bCT scanning. Real-time air kerma rate was measured during bCT scanning and fluctuated with breast anatomy. On average, point dose deposition was 7.1% greater than mean glandular dose.
A technique to generate a two-compartment anthropomorphic breast phantom from bCT images has been demonstrated. The phantom is the first, to our knowledge, to accurately model the uncompressed pendant breast and the glandular tissue distribution for a specific patient. The modular design of the phantom allows for studies of a single breast segment and the entire breast volume. Insertion of other devices, materials, and tissues of interest into the phantom provide a robust platform for future breast imaging and dosimetry studies.
Computed tomography; breast; phantom
The aim of the study was to estimate the dose at the reference point applying an aSi-EPID device in the course of patient treatment.
Materials and methods
The method assumes direct proportionality between EPID signal and dose delivered to the patient reference point during the treatment session. The procedure consists of treatment plan calculation for the actual patient in the arc technique. The plan was realized with an elliptic water-equivalent phantom. An ionization chamber inside the phantom measured the dose delivered to the reference point. Simultaneously, the EPID matrix measured the CU distribution. EPID signal was also registered during patient irradiation with the same treatment plan. The formula for in vivo dose calculation was based on the CU(g) function, EPID signal registered during therapy and the relation between the dose and EPID signal level measured for the phantom. In vivo dose was compared with dose planned with the treatment planning system.
Irradiation was performed with a Clinac accelerator by Varian Medical Systems in the RapidArc technique. The Clinac was equipped with an EPID matrix (electronic portal image device) of aSi-1000. Treatment plans were calculated with the Eclipse/Helios system. The phantom was a Scanditronix/Wellhöfer Slab phantom, and the ionization chamber was a 0.6 ccm PTW chamber.
In vivo dose calculations were performed for five patients. Planned dose at the reference point was 2 Gy for each treatment plan. Mean in vivo dose was in the range of 1.96–2.09.
Our method was shown to be appropriate for in vivo dose evaluation in the RapidArc technique.
Arc technique; Portal dosimetry; IMRT; Verification; Electronic portal imaging device
In 131I SPECT, image quality and quantification accuracy are degraded by object scatter as well as scatter and penetration in the collimator. The characterization of energy and spatial distributions of scatter and penetration performed in this study by Monte Carlo simulation will be useful for the development and evaluation of techniques that compensate for such events in 131I imaging.
First, to test the accuracy of the Monte Carlo model, simulated and measured data were compared for both a point source and a phantom. Next, simulations to investigate scatter and penetration were performed for four geometries: point source in air, point source in a water-filled cylinder, hot sphere in a cylinder filled with nonradioactive water, and hot sphere in a cylinder filled with radioactive water. Energy spectra were separated according to order of scatter, type of interaction, and γ-ray emission energy. A preliminary evaluation of the triple-energy window (TEW) scatter correction method was performed.
The accuracy of the Monte Carlo model was verified by the good agreement between measured and simulated energy spectra and radial point spread functions. For a point source in air, simulations show that 73% of events in the photopeak window had either scattered in or penetrated the collimator, indicating the significance of collimator interactions. For a point source in a water-filled phantom, the separated energy spectra showed that a 20% photopeak window can be used to eliminate events that scatter more than two times in the phantom. For the hot sphere phantoms, it was shown that in the photopeak region the spectrum shape of penetration events is very similar to that of primary (no scatter and no penetration) events. For the hot sphere regions of interest, the percentage difference between true scatter counts and the TEW estimate of scatter counts was <12%.
In 131I SPECT, object scatter as well as collimator scatter and penetration are significant. The TEW method provides a reasonable correction for scatter, but the similarity between the 364-keV primary and penetration energy spectra makes it difficult to compensate for these penetration events using techniques that are based on spectral analysis.
scatter correction; penetration; 131I imaging; SPECT; Monte Carlo simulated data
Background and purpose
Stereotactic lung radiotherapy (SLRT) has emerged as a curative treatment for medically inoperable patients with early-stage non-small cell lung cancer (NSCLC) and the use of intensity-modulated radiotherapy (IMRT) and volumetric modulated arc treatments (VMAT) have been proposed as the best practical approaches for the delivery of SLRT. However, a large number of narrow field shapes are needed in the dose delivery of intensity-modulated techniques and the probability of underdosing the tumour periphery increases as the effective field size is decreased. The purpose of this study was to evaluate small lung tumour doses irradiated by intensity-modulated techniques to understand the risk for dose calculation errors in precision radiotherapy such as SLRT.
Materials and methods
The study was executed with two heterogeneous phantoms with targets of Ø1.5 and Ø4.0 cm. Dose distributions in the simulated tumours delivered by small sliding window apertures (SWAs), IMRT and RapidArc treatment plans were measured with radiochromic film. Calculation algorithms of pencil beam convolution (PBC) and anisotropic analytic algorithm (AAA) were used to calculate the corresponding dose distributions.
Peripheral doses of the tumours were decreased as SWA decreased, which was not modelled by the calculation algorithms. The smallest SWA studied was 2 mm, which reduced the 90% isodose line width by 4.2 mm with the Ø4.0 cm tumour as compared to open field irradiation. PBC was not able to predict the dose accurately as the gamma evaluation failed to meet the criteria of ±3%/±1 mm on average in 61% of the defined volume with the smaller tumour. With AAA the corresponding value was 16%. The dosimetric inaccuracy of AAA was within ±3% with the optimized treatment plans of IMRT and RapidArc. The exception was the clinical RapidArc plan with dose overestimation of 4%.
Overall, the peripheral doses of the simulated lung tumours were decreased by decreasing the SWA. To achieve adequate surface dose coverage to small lung tumours with a difference less than 1 mm in the isodose line radius between the open and modulated field, a larger than 6 mm SWA should be used in the dose delivery of SLRT.
Stereotactic body radiotherapy; Lung cancer; IMRT; Heterogeneity; Surface dose
Because of the adverse effects of ionizing radiation on fetuses, prior to radiotherapy of pregnant patients, fetal dose should be estimated. Fetal dose has been studied by several authors in different depths in phantoms with various abdomen thicknesses (ATs). In this study, the effect of maternal AT and depth in fetal dosimetry was investigated, using peripheral dose (PD) distribution evaluations. A BEAMnrc model of Oncor linac using out of beam components was used for dose calculations in out of field border. A 6 MV photon beam was used to irradiate a chest phantom. Measurements were done using EBT2 radiochromic film in a RW3 phantom as abdomen. The followings were measured for different ATs: Depth PD profiles at two distances from the field's edge, and in-plane PD profiles at two depths. The results of this study show that PD is depth dependent near the field's edge. The increase in AT does not change PD depth of maximum and its distribution as a function of distance from the field's edge. It is concluded that estimating the maximum fetal dose, using a flat phantom, i.e., without taking into account the AT, is possible. Furthermore, an in-plane profile measured at any depth can represent the dose variation as a function of distance. However, in order to estimate the maximum PD the depth of Dmax in out of field should be used for in-plane profile measurement.
Fetal dosimetry; Monte Carlo; peripheral dose; radiotherapy
This study introduces a charge coupled device (CCD) area detector based optical-computed tomography (optical-CT) scanner for comprehensive verification of radiation dose distributions recorded in nonscattering radiochromic dosimeters. Defining characteristics include: (i) a very fast scanning time of ~5 min to acquire a complete three-dimensional (3D) dataset, (ii) improved image formation through the use of custom telecentric optics, which ensures accurate projection images and minimizes artifacts from scattered and stray-light sources, and (iii) high resolution (potentially 50 μm) isotropic 3D dose readout. The performance of the CCD scanner for 3D dose readout was evaluated by comparison with independent 3D readout from the single laser beam OCTOPUS™-scanner for the same PRESAGE™ dosimeters. The OCTOPUS™ scanner was considered the “gold standard” technique in light of prior studies demonstrating its accuracy. Additional comparisons were made against calculated dose distributions from the ECLIPSE treatment-planning system. Dose readout for the following treatments were investigated: (i) a single rectangular beam irradiation to investigate small field and very steep dose gradient dosimetry away from edge effects, (ii) a 2-field open beam parallel-opposed irradiation to investigate dosimetry along steep dose gradients, and (iii) a 7-field intensity modulated radiation therapy (IMRT) irradiation to investigate dosimetry for complex treatment delivery involving modulation of fluence and for dosimetry along moderate dose gradients. Dose profiles, dose-difference plots, and gamma maps were employed to evaluate quantitative estimates of agreement between independently measured and calculated dose distributions. Results indicated that dose readout from the CCD scanner was in agreement with independent gold-standard readout from the OCTOPUS™-scanner as well as the calculated ECLIPSE dose distribution for all treatments, except in regions within a few millimeters of the edge of the dosimeter, where edge artifact is predominant. Agreement of line profiles was observed, even along steep dose gradients. Dose difference plots indicated that the CCD scanner dose readout differed from the OCTOPUS™ scanner readout and ECLIPSE calculations by ~10% along steep dose gradients and by ~5% along moderate dose gradients. Gamma maps (3% dose-difference and 3 mm distance-toagreement acceptance criteria) revealed agreement, except for regions within 5 mm of the edge of the dosimeter where the edge artifact occurs. In summary, the data demonstrate feasibility of using the fast, high-resolution CCD scanner for comprehensive 3D dosimetry in all applications, except where dose readout is required close to the edges of the dosimeter. Further work is ongoing to reduce this artifact.
3D dosimetry; optical-CT; PRESAGE™; radiation; quality assurance; telecentric
A tracking system has been developed to provide real-time feedback of skin dose and dose rate during interventional fluoroscopic procedures. The dose tracking system (DTS) calculates the radiation dose rate to the patient’s skin using the exposure technique parameters and exposure geometry obtained from the x-ray imaging system digital network (Toshiba Infinix) and presents the cumulative results in a color mapping on a 3D graphic of the patient. We performed a number of tests to verify the accuracy of the dose representation of this system. These tests included comparison of system–calculated dose-rate values with ionization-chamber (6 cc PTW) measured values with change in kVp, beam filter, field size, source-to-skin distance and beam angulation. To simulate a cardiac catheterization procedure, the ionization chamber was also placed at various positions on an Alderson Rando torso phantom and the dose agreement compared for a range of projection angles with the heart at isocenter. To assess the accuracy of the dose distribution representation, Gafchromic film (XR-RV3, ISP) was exposed with the beam at different locations. The DTS and film distributions were compared and excellent visual agreement was obtained within the cm-sized surface elements used for the patient graphic. The dose (rate) values agreed within about 10% for the range of variables tested. Correction factors could be applied to obtain even closer agreement since the variable values are known in real-time. The DTS provides skin-dose values and dose mapping with sufficient accuracy for use in monitoring diagnostic and interventional x-ray procedures.
skin dose; dosimetry; radiation safety; cardiac fluoroscopic procedures; fluoroscopic dose; dose tracking; real-time dosimetry; fluoroscopic interventional procedures
This work reports the ongoing development of a combination applicator for simultaneous heating of superficial tissue disease using a 915 MHz DCC (dual concentric conductor) array and High Dose Rate (HDR) brachytherapy delivered via an integrated conformal catheter array. The progress includes engineering design changes in the waterbolus, DCC configurations and fabrication techniques of the conformal multilayer applicator. The dosimetric impact of the thin copper DCC array is also assessed. Steady state fluid dynamics of the new waterbolus bag indicates nearly uniform flow with less than 1°C variation across a large (19×32cm) bolus. Thermometry data of the torso phantom acquired with computer controlled movement of fiberoptic temperature probes inside thermal mapping catheters indicate feasibility of real time feedback control for the DCC array. MR (magnetic resonance) scans of a torso phantom indicate that the waterbolus thickness across the treatment area is controlled by the pressure applied by the surrounding inflatable airbladder and applicator securing straps. The attenuation coefficient of the DCC array was measured as 3± 0.001% and 2.95±0.03 % using an ion chamber and OneDose™ dosimeters respectively. The performance of the combination applicator on patient phantoms provides valuable feedback to optimize the applicator prior use in the patient clinic.
hyperthermia; brachytherapy; thermoradiotherapy; superficial disease; chestwall recurrence
To investigate the clinical usage of dose verification of Helical Tomotherapy plans by using 2D-array ion chambers, and to develop an efficient way to validate the dose delivered for the patients during treatments.
Materials and Methods:
A pixel-segmented ionisation chamber device, IMRT MatriXX™ and Multicube™ phantom from IBA were used on ten selected Tomotherapy IMRT/IGRT head and neck plans in this study. The combined phantom was set up to measure the dose distribution from coronal and sagittal planes. The setup of phantom was guided for verifying the correction position by pre-treatment Tomotherapy MVCT images. After the irradiation, the measured dose distributions of coronal and sagittal planes were compared with those from calculation by the planning system for cross verification. The results were evaluated by the absolute and relative doses as well as Gamma (γ) function. The feasibility of the different measuring methods was studied for this rotational treatment technique.
The dose distributions measured by the MatriXX 2D array were in good agreements with plans calculated by Tomotherapy planning system. The discrepancy between the measured dose and predicted dose in the selected points was within ±3%. In the comparison of the pixel-segmented ionisation chamber versus treatment planning system using the 3 mm/3% γ-function criteria, the mean passing rates of 2 mm dose grid with γ-parameter ≤1 were 97.37% and 96.91%, in two orthogonal planes (coronal and sagittal directions), respectively.
MatriXX with Multicube is a new system created for rotational delivery quality assurance (QA) and found to be reliable to measure both absolute dose and relative dose distributions, simultaneously. It achieves the goal of an efficient and accurate dosimetry validation method of the helical delivery pattern for the Helical Tomotherapy IMRT planning.
Tomotherapy; dose verification; IMRT; radiation therapy; QA
The purpose of this study is to develop a total body irradiation technique that does not require additional devices or sophisticated processes to overcome the space limitation of a small treatment room. The technique aims to deliver a uniform dose to the entire body while keeping the lung dose within the tolerance level. The technique treats the patient lying on the floor anteriorly and posteriorly. For each AP/PA treatment, two complementary fields with dynamic field edges are matched over an overlapped region defined by the marks on the body surface. A compensator, a spoiler, and lung shielding blocks were used during the treatment. Moreover, electron beams were used to further boost the chest wall around the lungs. The technique was validated in a RANDO phantom using GAFCHROMIC films. Dose ratios at different body sites along the midline ranged from 0.945 to 1.076. The dose variation in the AP direction ranged from 96.0% to 104.6%. The dose distribution in the overlapped region ranged from 98.5% to 102.8%. Lateral dose profiles at abdomen and head revealed 109.8% and 111.7% high doses, respectively, at the body edges. The results confirmed that the technique is capable of delivering a uniform dose distribution to the midline of the body in a small treatment room while keeping the lung dose within the tolerance level.
Several tools are used for the dosimetric verification of intensity-modulated arc therapy (IMAT) treatment delivery. However, limited information is available for composite on-line evaluation of these tools. The purpose of this study was to evaluate the dosimetric verification of IMAT treatment plans using a 2D diode array detector (2D array), radiochromic film (RCF) and radiosensitive polymer gel dosimeter (RPGD). The specific verification plans were created for IMAT for two prostate cancer patients by use of the clinical treatment plans. Accordingly, the IMAT deliveries were performed with the 2D array on a gantry-mounting device, RCF in a cylindrical acrylic phantom, and the RPGD in two cylindrical phantoms. After the irradiation, the planar dose distributions from the 2D array and the RCFs, and the 3D dose distributions from the RPGD measurements were compared with the calculated dose distributions using the gamma analysis method (3% dose difference and 3-mm distance-to-agreement criterion), dose-dependent dose difference diagrams, dose difference histograms, and isodose distributions. The gamma passing rates of 2D array, RCFs and RPGD for one patient were 99.5%, 96.5% and 93.7%, respectively; the corresponding values for the second patient were 97.5%, 92.6% and 92.9%. Mean percentage differences between the RPGD measured and calculated doses in 3D volumes containing PTVs were –0.29 ± 7.1% and 0.97 ± 7.6% for the two patients, respectively. In conclusion, IMAT prostate plans can be delivered with high accuracy, although the 3D measurements indicated less satisfactory agreement with the treatment plans, mainly due to the dosimetric inaccuracy in low-dose regions of the RPGD measurements.
IMAT; dosimetric verification; 3D dose measurement; radiosensitive polymer gel; radiochromic films; 2D array detector
The objective of this paper is threefold: (1) to establish sensitivity of XRQA and EBT radiochromic films to fast neutron exposure; (2) to develop a film response to radiation dose calibration curve and (3) to investigate a two-dimensional (2D) film dosimetry technique for use in establishing an experimental setup for a radiobiological irradiation of mice and to assess the dose to the mice in this setup. The films were exposed to a 10 MeV neutron beam via the 2H(d,n)3He reaction. The XRQA film response was a factor of 1.39 greater than EBT film response to the 10 MeV neutron beam when exposed to a neutron dose of 165 cGy. A film response-to-soft tissue dose calibration function was established over a range of 0–10 Gy and had a goodness of fit of 0.9926 with the calibration data. The 2D film dosimetry technique estimated the neutron dose to the mice by measuring the dose using a mouse phantom and by placing a piece of film on the exterior of the experimental mouse setup. The film results were benchmarked using Monte Carlo and aluminum (Al) foil activation measurements. The radiochromic film, Monte Carlo and Al foil dose measurements were strongly correlated, and the film within the mouse phantom agreed to better than 7% of the externally mounted films. These results demonstrated the potential application of radiochromic films for passive 2D neutron dosimetry.
The Monte Carlo code GEANT4 was used to implement first steps towards a treatment planning program for fast-neutron therapy at the FRM II research reactor in Garching, Germany. Depth dose curves were calculated inside a water phantom using measured primary neutron and simulated primary photon spectra and compared with depth dose curves measured earlier. The calculations were performed with GEANT4 in two different ways, simulating a simple box geometry and splitting this box into millions of small voxels (this was done to validate the voxelisation procedure that was also used to voxelise the human body).
In both cases, the dose distributions were very similar to those measured in the water phantom, up to a depth of 30 cm. In order to model the situation of patients treated at the FRM II MEDAPP therapy beamline for salivary gland tumors, a human voxel phantom was implemented in GEANT4 and irradiated with the implemented MEDAPP neutron and photon spectra. The 3D dose distribution calculated inside the head of the phantom was similar to the depth dose curves in the water phantom, with some differences that are explained by differences in elementary composition. The lateral dose distribution was studied at various depths. The calculated cumulative dose volume histograms for the voxel phantom show the exposure of organs at risk surrounding the tumor.
In order to minimize the dose to healthy tissue, a conformal treatment is necessary. This can only be accomplished with the help of an advanced treatment planning system like the one developed here. Although all calculations were done for absorbed dose only, any biological dose weighting can be implemented easily, to take into account the increased radiobiological effectiveness of neutrons compared to photons.
To study temperature and thermal dose distributions of ThermoBrachytherapy Surface Applicators (TBSA) developed for concurrent or sequential high dose rate (HDR) brachytherapy and microwave hyperthermia treatment of chest wall recurrence and other superficial disease.
A steady state thermodynamics model coupled with the fluid dynamics of water bolus and electromagnetic radiation of hyperthermia applicator is used to characterize the temperature distributions achievable with TBSA applicators in an elliptical phantom model of the human torso. Power deposited by 915 MHz conformal microwave array (CMA) applicators is used to assess the specific absorption rate (SAR) distributions of rectangular (500 cm2) and L-shaped (875 cm2) TBSA. The SAR distribution in tissue and fluid flow distribution inside the Dual-Input Dual-Output (DIDO) water bolus are coupled to solve the steady state temperature and thermal dose distributions of rectangular TBSA (R-TBSA) for superficial tumor targets extending 10–15 mm beneath the skin surface. Thermal simulations are carried out for a range of bolus inlet temperature (Tb=38–43°C), water flow rate (Qb=2–4 L/min) and tumor blood perfusion (ωb=2–5 kg/m3/s) to characterize their influence on thermal dosimetry.
Steady state SAR patterns of R- and L-TBSA demonstrate the ability to produce conformal and localized power deposition inside tumor target sparing surrounding normal tissues and nearby critical organs. Acceptably low variation in tissue surface cooling and surface temperature homogeneity was observed for the new DIDO bolus at 2 L/min water flow rate. Temperature depth profiles and thermal dose volume histograms indicate bolus inlet temperature (Tb) to be the most influential factor on thermal dosimetry. A 42 °C water bolus was observed to be the optimal choice for superficial tumors extending 10–15 mm from the surface even under significant blood perfusion. Lower bolus temperature may be chosen to reduce thermal enhancement ratio (TER) in the most sensitive skin where maximum radiation dose is delivered and to extend thermal enhancement of radiation dose deeper.
This computational study indicates that well-localized elevation of tumor target temperature to 40–44 °C can be accomplished by large surface-conforming TBSA applicators using appropriate selection of coupling bolus temperature.
chestwall recurrence; conformal thermal therapy; hyperthermia; brachytherapy; thermoradiotherapy; specific absorption rate; thermal dose