Volume 33, Issue 9, September 2006
- radiation therapy physics
- radiation imaging physics
- radiation measurement physics
- magnetic resonance physics
- nuclear medicine physics
- optical measurement physics
- ultrasound physics
- tissue measurements
- radiation protection physics
- radiation biology
- special topic
- books and publications
Index of content:
33(2006); http://dx.doi.org/10.1118/1.2179008View Description Hide Description
- RADIATION THERAPY PHYSICS
33(2006); http://dx.doi.org/10.1118/1.2241993View Description Hide Description
The accuracy of delivering gated-radiation therapy to lungtumors using an external respiratory surrogate relies on not only interfractional and intrafractional reproducibility, but also a strong correlation between external motion and internal tumor motion. The purpose of this work was to use the cine images acquired by four-dimensional computed tomography acquisition protocol to study the relation between external surface motion and internal tumor motion. The respiratory phase information of tumor motion and chest wall motion was measured on the cine images using a proposed region-of-interest (ROI) method and compared to measurement of an external respiratory monitoring device. On eight lung patient data sets, the phase shifts were measured between (1) the signal of a real-time positioning-management (RPM) respiratory monitoring device placed in the abdominal region and four surface locations on the chest wall, (2) the RPM signal in the abdominal region and tumor motions, and (3) chest wall surface motions and tumor motions. Respiratory waveforms measured at different surface locations during the same respiratory cycle often varied and had significant phase shifts. Seven of the 8 patients showed the abdominal motion leading chest wall motion. The best correlation (smallest phase shift) was found between the abdominal motion and the superior-inferior (S-I) tumor motion. A wide range of phase shifts was observed between external surface motion and tumor anterior-posterior (A-P)/lateral motion. The result supported the placement of the RPM block in the abdominal region and suggested that during a gated therapy utilizing the RPM system, it is necessary to place the RPM block at the same location as it is during treatment simulation in order to reduce potential errors introduced by the position of the RPM block. Correlations between external motions and lateral/A-P tumor motions were inconclusive due to a combination of patient selection and the limitation of the ROI method.
Accuracy of patient dose calculation for lung IMRT: A comparison of Monte Carlo, convolution/superposition, and pencil beam computations33(2006); http://dx.doi.org/10.1118/1.2241992View Description Hide Description
The accuracy of dose computation within the lungs depends strongly on the performance of the calculation algorithm in regions of electronic disequilibrium that arise near tissue inhomogeneities with large density variations. There is a lack of data evaluating the performance of highly developed analytical dose calculation algorithms compared to Monte Carlo computations in a clinical setting. We compared full Monte Carlo calculations (performed by our Monte Carlodose engine MCDE) with two different commercial convolution/superposition (CS) implementations (Pinnacle-CS and Helax-TMS’s collapsed cone model Helax-CC) and one pencil beam algorithm (Helax-TMS’s pencil beam model Helax-PB) for 10 intensity modulated radiation therapy(IMRT)lungcancer patients. Treatment plans were created for two photon beam qualities (6 and ). For each dose calculation algorithm, patient, and beam quality, the following set of clinically relevant dose-volume values was reported: (i) minimal, median, and maximal dose (, , and ) for the gross tumor and planning target volumes (GTV and PTV); (ii) the volume of the lungs (excluding the GTV) receiving at least 20 and ( and ) and the mean lungdose; (iii) the 33rd percentile dose and delivered to the heart and the expanded esophagus; and (iv) for the expanded spinal cord. Statistical analysis was performed by means of one-way analysis of variance for repeated measurements and Tukey pairwise comparison of means. Pinnacle-CS showed an excellent agreement with MCDE within the target structures, whereas the best correspondence for the organs at risk (OARs) was found between Helax-CC and MCDE. Results from Helax-PB were unsatisfying for both targets and OARs. Additionally, individual patient results were analyzed. Within the target structures, deviations above 5% were found in one patient for the comparison of MCDE and Helax-CC, while all differences between MCDE and Pinnacle-CS were below 5%. For both Pinnacle-CS and Helax-CC, deviations from MCDE above 5% were found within the OARs: within the lungs for two and six patients for Pinnacle-CS, and within other OARs for two patients for Helax-CC (for of the heart and of the expanded esophagus) but only for . For one patient, all four algorithms were used to recompute the dose after replacing all computed tomography voxels within the patient’s skin contour by water. This made all differences above 5% between MCDE and the other dose calculation algorithms disappear. Thus, the observed deviations mainly arose from differences in particle transport modeling within the lungs, and the commissioning of the algorithms was adequately performed (or the commissioning was less important for this type of treatment). In conclusion, not one pair of the dose calculation algorithms we investigated could provide results that were consistent within 5% for all 10 patients for the set of clinically relevant dose-volume indices studied. As the results from both CS algorithms differed significantly, care should be taken when evaluating treatment plans as the choice of dose calculation algorithm may influence clinical results. Full Monte Carlo provides a great benchmarking tool for evaluating the performance of other algorithms for patient dose computations.
Dosimetric characteristics of a newly designed grid block for megavoltage photon radiation and its therapeutic advantage using a linear quadratic model33(2006); http://dx.doi.org/10.1118/1.2241998View Description Hide Description
Grid radiation therapy with megavoltage x-ray beam has been proven to be an effective technique for management of large, bulky malignant tumors. The clinical advantage of GRID therapy, combined with conventional radiation therapy, has been demonstrated using a prototype GRID block [Mohiuddin, Curtis, Grizos, and Komarnicky, Cancer66, 114–118 (1990)]. Recently, a new GRID block design with improved dosimetric properties has become commercially available from Radiation Product Design, Inc. (Albertive, MN). This GRID collimator consists of an array of focused apertures in a cerrobend block arranged in a hexagonal pattern having a circular cross-section with a diameter and center-to-center spacing of 14.3 and , respectively, in the plane of isocenter. In this project, dosimetric characteristics of the newly redesigned GRID block have been investigated for a Varian 21EX linear accelerator (Varian Associates, Palo Alto, CA). These determinations were performed using radiographic films, thermoluminescent dosimeters in Solid Water™ phantom materials, and an ionization chamber in water. The output factor, percentage depth dose, beam profiles, and isodose distributions of the GRID radiation as a function of field size and beam energy have been measured using both 6 and x-ray beams. In addition, the therapeutic advantage obtained from this treatment modality with the new GRID block design for a high, single fraction of dose has been calculated using the linear quadratic model with ratios for typical tumor and normal cells. These biological characteristics of the new GRID block design will also be presented.
33(2006); http://dx.doi.org/10.1118/1.2335490View Description Hide Description
In intensity modulated radiation treatments (IMRT), the position of the field edges and the modulation within the beam are often achieved with a multileaf collimator(MLC). During the MLCcalibration process, due to the finite accuracy of leaf position measurements, a systematic error may be introduced to leaf positions. Thereafter leaf positions of the MLC depend on the systematic error introduced on each leaf during MLCcalibration and on the accuracy of the leaf position control system (random errors). This study presents and evaluates two methods to predict the systematic errors on the leaf positions introduced during the MLCcalibration. The two presented methods are based on a series of electronic portal imaging device(EPID)measurements. A comparison with film measurements showed that the EPID could be used to measure leaf positions without introducing any bias. The first method, referred to as the “central leaf method,” is based on the method currently used at this center for MLC leaf calibration. It mimics the manner in which leaf calibration parameters are specified in the MLCcontrol system and consequently is also used by other centers. The second method, a new method proposed by the authors and referred to as the “individual leaf method,” involves the measurement of two positions for each leaf ( and ) and the interpolation and extrapolation from these two points to any other given position. The central leaf method and the individual leaf method predicted leaf positions at prescribed positions of , 0, 5, and within 2.3 and , respectively, with a standard deviation (SD) of 0.3 and , respectively. The individual leaf method provided a better prediction of the leaf positions than the central leaf method. Reproducibility tests for leaf positions of and were performed. The reproducibility was within on the same day and six weeks later (1 SD). Measurements at gantry angles of 0°, 90°, and 270° for leaf positions of and showed no significant effect of gravity. The individual leaf method could be used in various applications to improve the accuracy of radiotherapy treatment from planning to delivery. Three cases are discussed: IMRT beam verification, MLCcalibration and dose calculation.
33(2006); http://dx.doi.org/10.1118/1.2222076View Description Hide Description
We have compared four computational methods for quantifying the effect of set-up error and uncertainty on delivered doses to targets and organs at risk in the intensity modulated radiation therapytreatment of head and neck cancer. These four methods were direct simulation, simple convolution, plus two modified convolution approaches. Discrepancies of up to 20% in the equivalent uniform dose (EUD) between direct simulation and simple convolution were estimated for the relatively superficial parotid gland at a systematic set-up error of standard deviation and a random uncertainty of standard deviation. Truncated convolution agreed with direct simulation to within 6% for all situations studied. However, of the four methods, only direct simulation can quantify the range of outcomes (EUD) associated with a finite number of courses and fractions. Our results are particularly relevant to the design of dose escalation studies in head and neck cancer.
33(2006); http://dx.doi.org/10.1118/1.2256068View Description Hide Description
A series of measurements were undertaken using both high sensitivity radiochromic film and new lithium fluoride thermoluminescent dosimeters in a liquid water medium to define the radial dose function of beta emitting intravascular brachytherapy sources more accurately. These measurements of a single 5 French source pellet served to verify current Monte Carlo transport models and extrapolation chamber measurements of the radial dose function, thus providing the recommended independent published measurements for of these sources. A slight deviation in the published radial dose function at depth leads the authors to recommend that treatment planning be performed using updated values from current Monte Carlo transport models verified by measurements such as those shown in this investigation.
Determination of dosimetrical quantities used in microbeam radiation therapy (MRT) with Monte Carlo simulations33(2006); http://dx.doi.org/10.1118/1.2229422View Description Hide Description
Microbeam radiation therapy (MRT) is being performed by using an array of narrow rectangular x-ray beams (typical beam sizes ), positioned close to each other (typically separation), to irradiate a target tissue. The ratio of peak-to-valley doses (PVDR’s) in the composite dose distribution has been found to be strongly correlated with the normal tissue tolerance and the therapeutic effect of MRT. In this work a Monte Carlo(MC) study of the depth- and lateral-dose profiles in water for single x-ray microbeams of different shapes and energies has been performed with the MC code PENELOPE. The contributions to the dose deposition from different interaction types have been determined at different distances from the center of the microbeam. The dependence of the peak dose, in a water phantom, on the microbeam field size used in the preclinical trials, has been demonstrated. Composite dose distributions for an array of microbeams were obtained using superposition algorithms and PVDR’s were determined and compared with literature results obtained with other Monte Carlo codes. The dependence of the PVDR’s on microbeam width, x-ray energy used, and on the separation between adjacent microbeams has been studied in detail.
33(2006); http://dx.doi.org/10.1118/1.2229430View Description Hide Description
In conventional linear accelerators, the flattening filter provides a uniform lateral dose profile. In intensity modulated radiation therapy applications, however, the flatness of the photon field and hence the presence of a flattening filter, is not necessary. Removing the filter may provide some advantages, such as faster treatments and smaller out-of-field doses to the patients. In clinical accelerators the backscatteredradiationdose from the collimators must be taken into account when the dose to the target volume in the patient is being determined. In the case of a conventional machine, this backscatter is known to great precision. In a flattening filter free accelerator, however, the amount of backscatter may be different. In this study we determined the backscatter contribution to the monitor chamber signal in a flattening filter free clinical accelerator (Varian Clinac 21EX) with Monte Carlo simulations. We found that with the exception of very small fields in the photon mode, the contribution of backscatteredradiation to the monitor signal did not differ from that of conventional machines with a flattening filter. Hence, a flattening filter free clinical accelerator would not necessitate a different backscatter correction.
33(2006); http://dx.doi.org/10.1118/1.2222077View Description Hide Description
In image-guided adaptive radiotherapy, it is important to have the capability to automatically and accurately delineate the rectal wall, which is a major dose-limiting organ in prostate cancerradiotherapy. As image registration is a process to find the spatial correspondence between two images, a major challenge in intensity-based deformable image registration is to deal with the situation where no correspondence exists for some objects between the two images to be registered. One example is the variation of rectal contents due to the presence and absence of bowel gas. The intensity-based deformable image registration methods alone cannot create the correct spatial transformation if there is no correspondence between the source and target images. In this study we implemented an automatic image intensity modification procedure to create artificial gas pockets in the planning computed tomography(CT)images. A diffusion-based deformable image registration algorithm was developed to use an adaptive smoothing algorithm to better handle large organ deformations. The process was tested in 15 prostate cancer cases and 30 daily CTimages containing the largest distended rectums. The manually delineated rectums agreed well with the autodelineated rectums when using the image-intensity modification procedure.
Compact multileaf collimator for conformal and intensity modulated fast neutron therapy: Electromechanical design and validation33(2006); http://dx.doi.org/10.1118/1.2237506View Description Hide Description
The electromechanical properties of a 120-leaf, high-resolution, computer-controlled, fast neutronmultileaf collimator(MLC) are presented. The MLC replaces an aging, manually operated multirod collimator. The MLC leaves project in the isocentric plane perpendicular to the beam axis. A taper is included on the leaves matching beam divergence along one axis. The leaf projection width is chosen to give high-resolution conformality across the entire field. The maximum field size provided is . To reduce the interleaf transmission a blocking step is included. End-leaf steps totaling are also provided allowing opposing leaves to close off within the primary radiation beam. The neutronMLC also includes individual 45° and 60° automated universal tungsten wedges. The automated high-resolution neutroncollimation provides an increase in patient throughput capacity, enables a new modality, intensity modulated neutron therapy, and limits occupational radiation exposure by providing remote operation from a shielded console area.
A computer simulated phantom study of tomotherapy dose optimization based on probability density functions (PDF) and potential errors caused by low reproducibility of PDF33(2006); http://dx.doi.org/10.1118/1.2222331View Description Hide Description
Lungtumor motion trajectories measured by four-dimensional CT or dynamic MRI can be converted to a probability density function (PDF), which describes the probability of the tumor at a certain position, for PDF based treatment planning. Using this method in simulated sequential tomotherapy, we study the dose reduction of normal tissues and more important, the effect of PDF reproducibility on the accuracy of dosimetry. For these purposes, realistic PDFs were obtained from two dynamic MRI scans of a healthy volunteer within a 2 week interval. The first PDF was accumulated from a scan and the second PDF was calculated from variable scan times from (one breathing cycle) to . Optimized beam fluences based on the second PDF were delivered to the hypothetical gross target volume (GTV) of a lung phantom that moved following the first PDF. The reproducibility between two PDFs varied from low (78%) to high (94.8%) when the second scan time increased from to . When a highly reproducible PDF was used in optimization, the dose coverage of GTV was maintained; phantom lung receiving 10%–20% prescription dose was reduced by 40%–50% and the mean phantom lungdose was reduced by 9.6%. However, optimization based on PDF with low reproducibility resulted in a 50% underdosed GTV. The dosimetric error increased nearly exponentially as the PDF error increased. Therefore, although the dose of the tumor surrounding tissue can be theoretically reduced by PDF based treatment planning, the reliability and applicability of this method highly depend on if a reproducible PDF exists and is measurable. By correlating the dosimetric error and PDF error together, a useful guideline for PDF data acquisition and patient qualification for PDF based planning can be derived.
33(2006); http://dx.doi.org/10.1118/1.2259819View Description Hide Description
The effect of beam divergence on dose calculation via Monte Carlo generated dose kernels was investigated in an amorphous silicon electronic portal imaging device(EPID). The flat-panel detector was simulated in EGSnrc with an additional 3.0 cm water buildup. The model included details of the detector’simaging cassette and the front cover upstream of it. To approximate the effect of the EPID’s rear housing, a 2.1 cm air gap and 1.0 cm water slab were introduced into the simulation as equivalent backscatter material. Dose kernels were generated with an incident pencil beam of monoenergetic photons of energy 0.1, 2, 6, and 18 MeV. The orientation of the incident pencil beam was varied from 0° to 14° in 2° increments. Dose was scored in the phosphor layer of the detector in both cylindrical (at 0°) and Cartesian (at ) geometries. To reduce statistical fluctuations in the Cartesian geometry simulations at large radial distances from the incident pencil beam, the voxels were first averaged bilaterally about the pencil beam and then combined into concentric square rings of voxels. Profiles of the EPIDdose kernels displayed increasing asymmetry with increasing angle and energy. A comparison of the superposition (tilted kernels) and convolution (parallel kernels) dose calculation methods via the χ-comparison test (a derivative of the γ-evaluation) in worst-case-scenario geometries demonstrated an agreement between the two methods within 0.0784 cm (one pixel width) distance-to-agreement and up to a 1.8% dose difference. More clinically typical field sizes and source-to-detector distances were also tested, yielding at most a 1.0% dose difference and the same distance-to-agreement. Therefore, the assumption of parallel dose kernels has less than a 1.8% dosimetric effect in extreme cases and less than a 1.0% dosimetric effect in most clinically relevant situations and should be suitable for most clinical dosimetric applications. The resulting time difference for the parallel kernel assumption versus the tilted kernels was 10.5 s vs 18 h (a factor of approximately 6000), dependent on existing hardware and software details.
33(2006); http://dx.doi.org/10.1118/1.2222079View Description Hide Description
Respiratory motion causes errors when planning and delivering radiotherapy treatment to lungcancer patients. To reduce these errors, methods of acquiring and using four-dimensional computed tomography (4DCT) datasets have been developed. We have developed a novel method of constructing computational motion models from 4DCT. The motion models attempt to describe an average respiratory cycle, which reduces the effects of variation between different cycles. They require substantially less memory than a 4DCT dataset, are continuous in space and time, and facilitate automatic target propagation and combining of doses over the respiratory cycle. The motion models are constructed from CT data acquired in cine mode while the patient is free breathing (free breathing CT – FBCT). A “slab” of data is acquired at each couch position, with 3–4 contiguous slabs being acquired per patient. For each slab a sequence of 20 or 30 volumes was acquired over 20 seconds. A respiratory signal is simultaneously recorded in order to calculate the position in the respiratory cycle for each FBCT. Additionally, a high quality reference CT volume is acquired at breath hold. The reference volume is nonrigidly registered to each of the FBCT volumes. A motion model is then constructed for each slab by temporally fitting the nonrigid registration results. The value of each of the registration parameters is related to the position in the respiratory cycle by fitting an approximating B spline to the registration results. As an approximating function is used, and the data is acquired over several respiratory cycles, the function should model an average respiratory cycle. This can then be used to calculate the value of each degree of freedom at any desired position in the respiratory cycle. The resulting nonrigid transformation will deform the reference volume to predict the contents of the slab at the desired position in the respiratory cycle. The slab model predictions are then concatenated to produce a combined prediction over the entire region of interest. We have performed a number of experiments to assess the accuracy of the nonrigid registration results and the motion model predictions. The individual slab models were evaluated by expert visual assessment and the tracking of easily identifiable anatomical points. The combined models were evaluated by calculating the discontinuities between the transformations at the slab boundaries. The experiments were performed on five patients with a total of 18 slabs between them. For the point tracking experiments, the mean distance between where a clinician manually identified a point and where the registration results located the point, the target registration error (TRE), was . The mean distance between a manually identified point and the models prediction of the point’s location, the target model error (TME), was . The mean discontinuity between model predictions at the slab boundaries, the Continuity Error, was . The results show that the motion models perform with a level of accuracy comparable to the slice thickness of .
Three-dimensional portal image-based dose reconstruction in a virtual phantom for rapid evaluation of IMRT plans33(2006); http://dx.doi.org/10.1118/1.2241997View Description Hide Description
A new method for rapid evaluation of intensity modulated radiation therapy(IMRT) plans has been developed, using portal images for reconstruction of the dose delivered to a virtual three-dimensional (3D) phantom. This technique can replace an array of less complete but more time-consuming measurements. A reference dose calculation is first created by transferring an IMRT plan to a cylindrical phantom, retaining the treatment gantry angles. The isocenter of the fields is placed on or near the phantom axis. This geometry preserves the relative locations of high and low dose regions and has the required symmetry for the dosereconstruction. An electronic portal image (EPI) is acquired for each field, representing the dose in the midplane of a virtual phantom. The image is convolved with a kernel to correct for the lack of scatter, replicating the effect of the cylindrical phantom surrounding the dose plane. This avoids the need to calculate fluence. Images are calibrated to a reference field that delivers a known dose to the isocenter of this phantom. The 3D dose matrix is reconstructed by attenuation and divergence corrections and summed to create a dose matrix (PI-dose) on the same grid spacing as the reference calculation. Comparison of the two distributions is performed with a gradient-weighted 3D dose difference based on dose and position tolerances. Because of its inherent simplicity, the technique is optimally suited for detecting clinically significant variances from a planned dose distribution, rather than for use in the validation of IMRT algorithms. An analysis of differences between PI-dose and calculation, , compared to differences between conventional quality assurance (QA) and calculation, , was performed retrospectively for 20 clinical IMRT cases. PI-dose differences at the isocenter were in good agreement with ionization chamber differences (mean , standard deviation , against , , respectively). PI-dose plane differences had significantly less variance than film plane differences ( and 2.1%, respectively). Twenty-two further cases were evaluated using 3D EPI-dosimetry alone. The mean difference over volumes with doses above 80% of the isocenter value was , , and standard deviations of the distributions ranged from 1.0 to 2.0%. Verification time per plan, from initial calculation, delivery, dosereconstruction to evaluation, takes less than and is more than four times faster than conventional QA.
Novel lung IMRT planning algorithms with nonuniform dose delivery strategy to account for respiratory motion33(2006); http://dx.doi.org/10.1118/1.2335485View Description Hide Description
To effectively deliver radiationdose to lungtumors, respiratory motion has to be considered in treatment planning. In this paper we first present a new lungIMRT planning algorithm, referred as the dose shaping (DS) method, that shapes the dose distribution according to the probability distribution of the tumor over the breathing cycle to account for respiratory motion. In IMRT planning a dose-based convolution method was generally adopted to compensate for random organ motion by performing 4-D dose calculations using a tumor motion probability density function. We modified the CON-DOSE method to a dose volume histogram based convolution method (CON-DVH) that allows nonuniform dose distribution to account for respiratory motion. We implemented the two new planning algorithms on an in-house IMRT planning system that uses the (Varian, Palo Alto, CA) planning workstation as the dose calculation engine. The new algorithms were compared with (1) the conventional margin extension approach in which margin is generated based on the extreme positions of the tumor, (2) the dose-based convolution method, and (3) gating with residual motion. Dose volume histogram, tumor control probability, normal tissue complication probability, and mean lungdose were calculated and used to evaluate the relative performance of these approaches at the end-exhale phase of the respiratory cycle. We recruited six patients in our treatment planning study. The study demonstrated that the two new methods could significantly reduce the ipsilateral normal lungdose and outperformed the margin extension method and the dose-based convolution method. Compared with the gated approach that has the best performance in the low dose region, the two methods we proposed have similar potential to escalate tumordose, but could be more efficient because dose is delivered continuously.
33(2006); http://dx.doi.org/10.1118/1.2335486View Description Hide Description
Radiotherapy planning involves inherent tradeoffs: the primary mission, to treat the tumor with a high, uniform dose, is in conflict with normal tissue sparing. We seek to understand these tradeoffs on a case-to-case basis, by computing for each patient a database of Pareto optimal plans. A treatment plan is Pareto optimal if there does not exist another plan which is better in every measurable dimension. The set of all such plans is called the Pareto optimal surface. This article presents an algorithm for computing well distributed points on the (convex) Pareto optimal surface of a multiobjective programming problem. The algorithm is applied to intensity-modulated radiation therapy inverse planning problems, and results of a prostate case and a skull base case are presented, in three and four dimensions, investigating tradeoffs between tumor coverage and critical organ sparing.
33(2006); http://dx.doi.org/10.1118/1.2335488View Description Hide Description
The purpose of this study was to investigate the specific problems associated with photondose calculations in points located at a distance from the central beam axis. These problems are related to laterally inhomogeneous energy fluence distributions and spectral variations causing a lateral shift in the beam quality, commonly referred to as off-axis softening (OAS). We have examined how the dose calculation accuracy is affected when enabling and disabling explicit modeling of these two effects. The calculations were performed using a pencil kernel dose calculation algorithm that facilitates modeling of OAS through laterally varying kernel properties. Together with a multisource model that provides the lateral energy fluence distribution this generates the total dose output, i.e., the dose per monitor unit, at an arbitrary point of interest. The dose calculation accuracy was evaluated through comparisons with 264 measured output factors acquired at 5, 10, and depth in four different megavoltage photon beams. The measurements were performed up to from the central beam axis, inside square fields of varying size and position. The results show that calculations including explicit modeling of OAS were considerably more accurate, up to 4%, than those ignoring the lateral beam quality shift. The deviations caused by simplified head scatter modeling were smaller, but near the field edges additional errors close to 1% occurred. When enabling full physics modeling in the dose calculations the deviations display a mean value of , a standard deviation of 0.7%, and a maximum deviation of . Finally, the results were analyzed in order to quantify and model the inherent uncertainties that are present when leaving the central beam axis. The off-axis uncertainty component showed to increase with both off-axis distance and depth, reaching 1% (1 standard deviation) at depth.
Low-dose megavoltage cone-beam computed tomography for lung tumors using a high-efficiency image receptor33(2006); http://dx.doi.org/10.1118/1.2222075View Description Hide Description
We report on the capabilities of a low-dose megavoltage cone-beam computed tomography (MV CBCT) system. The high-efficiency image receptor consists of a photodiode array coupled to a scintillator composed of individual crystals. The CBCT system uses the 6 MV beam from a linear accelerator. A synchronization circuit allows us to limit the exposure to one beam pulse [0.028 monitor units (MU)] per projection image. 150–500 images ( total) are collected during a one-minute scan and reconstructed using a filtered backprojection algorithm. Anthropomorphic and contrast phantoms are imaged and the contrast-to-noise ratio of the reconstruction is studied as a function of the number of projections and the error in the projection angles. The detectordose response is linear ( value 0.9989). A 2% electron density difference is discernible using 460 projection images and a total exposure of (corresponding to a maximum absorbed dose of about in a patient). We present first patient images acquired with this system. Tumors in lung are clearly visible and skeletal anatomy is observed in sufficient detail to allow reproducible registration with the planning kV CTimages. The MV CBCT system is shown to be capable of obtaining good quality three-dimensional reconstructions at relatively low dose and to be clinically usable for improving the accuracy of radiotherapy patient positioning.
33(2006); http://dx.doi.org/10.1118/1.2241996View Description Hide Description
We developed an analytical method for determining the maximum acceptable grid size for discrete dose calculation in proton therapy treatment plan optimization, so that the accuracy of the optimized dose distribution is guaranteed in the phase of dose sampling and the superfluous computational work is avoided. The accuracy of dose sampling was judged by the criterion that the continuous dose distribution could be reconstructed from the discrete dose within a 2% error limit. To keep the error caused by the discrete dose sampling under a 2% limit, the dose grid size cannot exceed a maximum acceptable value. The method was based on Fourier analysis and the Shannon-Nyquist sampling theorem as an extension of our previous analysis for photon beam intensity modulated radiation therapy [J. F. Dempsey, H. E. Romeijn, J. G. Li, D. A. Low, and J. R. Palta, Med. Phys.32, 380–388 (2005)]. The proton beam model used for the analysis was a near mono-energetic (of width about 1% the incident energy) and monodirectional infinitesimal (nonintegrated) pencil beam in water medium. By monodirection, we mean that the proton particles are in the same direction before entering the water medium and the various scattering prior to entrance to water is not taken into account. In intensity modulated proton therapy, the elementary intensity modulation entity for proton therapy is either an infinitesimal or finite sized beamlet. Since a finite sized beamlet is the superposition of infinitesimal pencil beams, the result of the maximum acceptable grid size obtained with infinitesimal pencil beam also applies to finite sized beamlet. The analytic Bragg curve function proposed by Bortfeld [T. Bortfeld, Med. Phys.24, 2024–2033 (1997)] was employed. The lateral profile was approximated by a depth dependent Gaussian distribution. The model included the spreads of the Bragg peak and the lateral profiles due to multiple Coulomb scattering. The dependence of the maximum acceptable dose grid size on the orientation of the beam with respect to the dose grid was also investigated. The maximum acceptable dose grid size depends on the gradient of dose profile and in turn the range of proton beam. In the case that only the phantom scattering was considered and that the beam was aligned with the dose grid, grid sizes from 0.4 to were required for proton beams with ranges from 2 to for 2% error limit at the Bragg peak point. A near linear relation between the maximum acceptable grid size and beam range was observed. For this analysis model, the resolution requirement was not significantly related to the orientation of the beam with respect to the grid.