Purpose:
To demonstrate knowledge‐based 3D dose prediction for external beam radiotherapy.
Methods:
Using previously treated plans as training data, an artificial neural network (ANN) was trained to ...predict a dose matrix based on patient‐specific geometric and planning parameters, such as the closest distance (r) to planning target volume (PTV) and organ‐at‐risks (OARs). Twenty‐three prostate and 43 stereotactic radiosurgery/radiotherapy (SRS/SRT) cases with at least one nearby OAR were studied. All were planned with volumetric‐modulated arc therapy to prescription doses of 81 Gy for prostate and 12–30 Gy for SRS. Using these clinically approved plans, ANNs were trained to predict dose matrix and the predictive accuracy was evaluated using the dose difference between the clinical plan and prediction, δD = Dclin − Dpred. The mean (〈δDr〉), standard deviation (σδDr), and their interquartile range (IQR) for the training plans were evaluated at a 2–3 mm interval from the PTV boundary (rPTV) to assess prediction bias and precision. Initially, unfiltered models which were trained using all plans in the cohorts were created for each treatment site. The models predict approximately the average quality of OAR sparing. Emphasizing a subset of plans that exhibited superior to the average OAR sparing during training, refined models were created to predict high‐quality rectum sparing for prostate and brainstem sparing for SRS. Using the refined model, potentially suboptimal plans were identified where the model predicted further sparing of the OARs was achievable. Replans were performed to test if the OAR sparing could be improved as predicted by the model.
Results:
The refined models demonstrated highly accurate dose distribution prediction. For prostate cases, the average prediction bias for all voxels irrespective of organ delineation ranged from −1% to 0% with maximum IQR of 3% over rPTV ∈ − 6, 30 mm. The average prediction error was less than 10% for the same rPTV range. For SRS cases, the average prediction bias ranged from −0.7% to 1.5% with maximum IQR of 5% over rPTV ∈ − 4, 32 mm. The average prediction error was less than 8%. Four potentially suboptimal plans were identified for each site and subsequent replanning demonstrated improved sparing of rectum and brainstem.
Conclusions:
The study demonstrates highly accurate knowledge‐based 3D dose predictions for radiotherapy plans.
On the use of SRIM for computing radiation damage exposure Stoller, R.E.; Toloczko, M.B.; Was, G.S. ...
Nuclear instruments & methods in physics research. Section B, Beam interactions with materials and atoms,
09/2013, Volume:
310
Journal Article
Peer reviewed
The SRIM (formerly TRIM) Monte Carlo simulation code is widely used to compute a number of parameters relevant to ion beam implantation and ion beam processing of materials. It also has the ...capability to compute a common radiation damage exposure unit known as atomic displacements per atom (dpa). Since dpa is a standard measure of primary radiation damage production, most researchers who employ ion beams as a tool for inducing radiation damage in materials use SRIM to determine the dpa associated with their irradiations. The use of SRIM for this purpose has been evaluated and comparisons have been made with an internationally-recognized standard definition of dpa, as well as more detailed atomistic simulations of atomic displacement cascades. Differences between the standard and SRIM-based dpa are discussed and recommendations for future usage of SRIM in radiation damage studies are made. In particular, it is recommended that when direct comparisons between ion and neutron data are intended, the Kinchin–Pease option of SRIM should be selected.
The charge of Task Group 186 (TG-186) is to provide guidance for early adopters of model-based dose calculation algorithms (MBDCAs) for brachytherapy (BT) dose calculations to ensure practice ...uniformity. Contrary to external beam radiotherapy, heterogeneity correction algorithms have only recently been made available to the BT community. Yet, BT dose calculation accuracy is highly dependent on scatter conditions and photoelectric effect cross-sections relative to water. In specific situations, differences between the current water-based BT dose calculation formalism (TG-43) and MBDCAs can lead to differences in calculated doses exceeding a factor of 10. MBDCAs raise three major issues that are not addressed by current guidance documents: (1) MBDCA calculated doses are sensitive to the dose specification medium, resulting in energy-dependent differences between dose calculated to water in a homogeneous water geometry (TG-43), dose calculated to the local medium in the heterogeneous medium, and the intermediate scenario of dose calculated to a small volume of water in the heterogeneous medium. (2) MBDCA doses are sensitive to voxel-by-voxel interaction cross sections. Neither conventional single-energy CT nor ICRU/ICRP tissue composition compilations provide useful guidance for the task of assigning interaction cross sections to each voxel. (3) Since each patient-source-applicator combination is unique, having reference data for each possible combination to benchmark MBDCAs is an impractical strategy. Hence, a new commissioning process is required. TG-186 addresses in detail the above issues through the literature review and provides explicit recommendations based on the current state of knowledge. TG-43-based dose prescription and dose calculation remain in effect, with MBDCA dose reporting performed in parallel when available. In using MBDCAs, it is recommended that the radiation transport should be performed in the heterogeneous medium and, at minimum, the dose to the local medium be reported along with the TG-43 calculated doses. Assignments of voxel-by-voxel cross sections represent a particular challenge. Electron density information is readily extracted from CT imaging, but cannot be used to distinguish between different materials having the same density. Therefore, a recommendation is made to use a number of standardized materials to maintain uniformity across institutions. Sensitivity analysis shows that this recommendation offers increased accuracy over TG-43. MBDCA commissioning will share commonalities with current TG-43-based systems, but in addition there will be algorithm-specific tasks. Two levels of commissioning are recommended: reproducing TG-43 dose parameters and testing the advanced capabilities of MBDCAs. For validation of heterogeneity and scatter conditions, MBDCAs should mimic the 3D dose distributions from reference virtual geometries. Potential changes in BT dose prescriptions and MBDCA limitations are discussed. When data required for full MBDCA implementation are insufficient, interim recommendations are made and potential areas of research are identified. Application of TG-186 guidance should retain practice uniformity in transitioning from the TG-43 to the MBDCA approach.
Mammography is a very well-established imaging modality for the early detection and diagnosis of breast cancer. However, since the introduction of digital imaging to the realm of radiology, more ...advanced, and especially tomographic imaging methods have been made possible. One of these methods, breast tomosynthesis, has finally been introduced to the clinic for routine everyday use, with potential to in the future replace mammography for screening for breast cancer. In this two part paper, the extensive research performed during the development of breast tomosynthesis is reviewed, with a focus on the research addressing the medical physics aspects of this imaging modality. This first paper will review the research performed on the issues relevant to the image acquisition process, including system design, optimization of geometry and technique, x-ray scatter, and radiation dose. The companion to this paper will review all other aspects of breast tomosynthesis imaging, including the reconstruction process.
In this work a novel plan optimization platform is presented where treatment is delivered efficiently and accurately in a single dynamically modulated arc. Improvements in patient care achieved ...through image-guided positioning and plan adaptation have resulted in an increase in overall treatment times. Intensity-modulated radiation therapy (IMRT) has also increased treatment time by requiring a larger number of beam directions, increased monitor units (MU), and, in the case of tomotherapy, a slice-by-slice delivery. In order to maintain a similar level of patient throughput it will be necessary to increase the efficiency of treatment delivery. The solution proposed here is a novel aperture-based algorithm for treatment plan optimization where dose is delivered during a single gantry arc of up to 360 deg. The technique is similar to tomotherapy in that a full 360 deg of beam directions are available for optimization but is fundamentally different in that the entire dose volume is delivered in a single source rotation. The new technique is referred to as volumetric modulated arc therapy (VMAT). Multileaf collimator (MLC) leaf motion and number of MU per degree of gantry rotation is restricted during the optimization so that gantry rotation speed, leaf translation speed, and dose rate maxima do not excessively limit the delivery efficiency. During planning, investigators model continuous gantry motion by a coarse sampling of static gantry positions and fluence maps or MLC aperture shapes. The technique presented here is unique in that gantry and MLC position sampling is progressively increased throughout the optimization. Using the full gantry range will theoretically provide increased flexibility in generating highly conformal treatment plans. In practice, the additional flexibility is somewhat negated by the additional constraints placed on the amount of MLC leaf motion between gantry samples. A series of studies are performed that characterize the relationship between gantry and MLC sampling, dose modeling accuracy, and optimization time. Results show that gantry angle and MLC sample spacing as low as 1 deg and 0.5 cm, respectively, is desirable for accurate dose modeling. It is also shown that reducing the sample spacing dramatically reduces the ability of the optimization to arrive at a solution. The competing benefits of having small and large sample spacing are mutually realized using the progressive sampling technique described here. Preliminary results show that plans generated with VMAT optimization exhibit dose distributions equivalent or superior to static gantry IMRT. Timing studies have shown that the VMAT technique is well suited for on-line verification and adaptation with delivery times that are reduced to
∼
1.5
–
3
min
for a 200 cGy fraction.
Purpose:
Different computed tomography (CT) reconstruction techniques offer different image quality attributes of resolution and noise, challenging the ability to compare their dose reduction ...potential against each other. The purpose of this study was to evaluate and compare the task‐based imaging performance of CT systems to enable the assessment of the dose performance of a model‐based iterative reconstruction (MBIR) to that of an adaptive statistical iterative reconstruction (ASIR) and a filtered back projection (FBP) technique.
Methods:
The ACR CT phantom (model 464) was imaged across a wide range of mA setting on a 64‐slice CT scanner (GE Discovery CT750 HD, Waukesha, WI). Based on previous work, the resolution was evaluated in terms of a task‐based modulation transfer function (MTF) using a circular‐edge technique and images from the contrast inserts located in the ACR phantom. Noise performance was assessed in terms of the noise‐power spectrum (NPS) measured from the uniform section of the phantom. The task‐based MTF and NPS were combined with a task function to yield a task‐based estimate of imaging performance, the detectability index (d′). The detectability index was computed as a function of dose for two imaging tasks corresponding to the detection of a relatively small and a relatively large feature (1.5 and 25 mm, respectively). The performance of MBIR in terms of the d′ was compared with that of ASIR and FBP to assess its dose reduction potential.
Results:
Results indicated that MBIR exhibits a variability spatial resolution with respect to object contrast and noise while significantly reducing image noise. The NPS measurements for MBIR indicated a noise texture with a low‐pass quality compared to the typical midpass noise found in FBP‐based CT images. At comparable dose, the d′ for MBIR was higher than those of FBP and ASIR by at least 61% and 19% for the small feature and the large feature tasks, respectively. Compared to FBP and ASIR, MBIR indicated a 46%–84% dose reduction potential, depending on task, without compromising the modeled detection performance.
Conclusions:
The presented methodology based on ACR phantom measurements extends current possibilities for the assessment of CT image quality under the complex resolution and noise characteristics exhibited with statistical and iterative reconstruction algorithms. The findings further suggest that MBIR can potentially make better use of the projections data to reduce CT dose by approximately a factor of 2. Alternatively, if the dose held unchanged, it can improve image quality by different levels for different tasks.
Purpose:
To evaluate the image quality of virtual monochromatic images synthesized from dual-source dual-energy computed tomography (CT) in comparison with conventional polychromatic single-energy CT ...for the same radiation dose.
Methods:
In dual-energy CT, besides the material-specific information, one may also synthesize monochromatic images at different energies, which can be used for routine diagnosis similar to conventional polychromatic single-energy images. In this work, the authors assessed whether virtual monochromatic images generated from dual-source CT scanners had an image quality similar to that of polychromatic single-energy images for the same radiation dose. First, the authors provided a theoretical analysis of the optimal monochromatic energy for either the minimum noise level or the highest iodine contrast to noise ratio (CNR) for a given patient size and dose partitioning between the low- and high-energy scans. Second, the authors performed an experimental study on a dual-source CT scanner to evaluate the noise and iodine CNR in monochromatic images. A thoracic phantom with three sizes of attenuating rings was used to represent four adult sizes. For each phantom size, three dose partitionings between the low-energy (80 kV) and the high-energy (140 kV) scans were used in the dual-energy scan. Monochromatic images at eight energies (40 to 110 keV) were generated for each scan. Phantoms were also scanned at each of the four polychromatic single energy (80, 100, 120, and 140 kV) with the same radiation dose.
Results:
The optimal virtual monochromatic energy depends on several factors: phantom size, partitioning of the radiation dose between low- and high-energy scans, and the image quality metrics to be optimized. With the increase of phantom size, the optimal monochromatic energy increased. With the increased percentage of radiation dose on the low energy scan, the optimal monochromatic energy decreased. When maximizing the iodine CNR in monochromatic images, the optimal energy was lower than that when minimizing noise level. When the total radiation dose was equally distributed between low and high energy in dual-energy scans, for minimum noise, the optimal energies were 68, 71, 74, and 77 keV for small, medium, large, and extra-large (xlarge) phantoms, respectively; for maximum iodine CNR, the optimal energies were 66, 68, 70, 72 keV. With the optimal monochromatic energy, the noise level was similar to and the CNR was better than that in a single-energy scan at 120 kV for the same radiation dose. Compared to an 80 kV scan, however, the iodine CNR in monochromatic images was lower for the small, medium, and large phantoms.
Conclusions:
In dual-source dual-energy CT, optimal virtual monochromatic energy depends on patient size, dose partitioning, and the image quality metric optimized. With the optimal monochromatic energy, the noise level was similar to and the iodine CNR was better than that in 120 kV images for the same radiation dose. Compared to single-energy 80 kV images, the iodine CNR in virtual monochromatic images was lower for small to large phantom sizes.
Purpose:
The aim of this study was to evaluate the appropriateness of using computed tomography (CT) to cone-beam CT (CBCT) deformable image registration (DIR) for the application of calculating the ...“dose of the day” received by a head and neck patient.
Methods:
NiftyReg is an open-source registration package implemented in our institution. The affine registration uses a Block Matching-based approach, while the deformable registration is a GPU implementation of the popular B-spline Free Form Deformation algorithm. Two independent tests were performed to assess the suitability of our registrations methodology for “dose of the day” calculations in a deformed CT. A geometric evaluation was performed to assess the ability of the DIR method to map identical structures between the CT and CBCT datasets. Features delineated in the planning CT were warped and compared with features manually drawn on the CBCT. The authors computed the dice similarity coefficient (DSC), distance transformation, and centre of mass distance between features. A dosimetric evaluation was performed to evaluate the clinical significance of the registrations errors in the application proposed and to identify the limitations of the approximations used. Dose calculations for the same intensity-modulated radiation therapy plan on the deformed CT and replan CT were compared. Dose distributions were compared in terms of dose differences (DD), gamma analysis, target coverage, and dose volume histograms (DVHs). Doses calculated in a rigidly aligned CT and directly in an extended CBCT were also evaluated.
Results:
A mean value of 0.850 in DSC was achieved in overlap between manually delineated and warped features, with the distance between surfaces being less than 2 mm on over 90% of the pixels. Deformable registration was clearly superior to rigid registration in mapping identical structures between the two datasets. The dose recalculated in the deformed CT is a good match to the dose calculated on a replan CT. The DD is smaller than 2% of the prescribed dose on 90% of the body's voxels and it passes a 2% and 2 mm gamma-test on over 95% of the voxels. Target coverage similarity was assessed in terms of the 95%-isodose volumes. A mean value of 0.962 was obtained for the DSC, while the distance between surfaces is less than 2 mm in 95.4% of the pixels. The method proposed provided adequate dose estimation, closer to the gold standard than the other two approaches. Differences in DVH curves were mainly due to differences in the OARs definition (manual vs warped) and not due to differences in dose estimation (dose calculated in replan CT vs dose calculated in deformed CT).
Conclusions:
Deforming a planning CT to match a daily CBCT provides the tools needed for the calculation of the “dose of the day” without the need to acquire a new CT. The initial clinical application of our method will be weekly offline calculations of the “dose of the day,” and use this information to inform adaptive radiotherapy (ART). The work here presented is a first step into a full implementation of a “dose-driven” online ART.
In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose ...calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) (192)Ir source and a virtual water phantom were designed, which can be imported into a TPS.
A hypothetical, generic HDR (192)Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic (192)Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs Oncentra(®) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS™ . Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR (192)Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods.
TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ace algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 ± 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ace at clinically relevant distances.
A hypothetical, generic HDR (192)Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.
PDF
