Volume 39, Issue 11, November 2012
- medical physics letters
- radiation therapy physics
- radiation imaging physics
- radiation measurement physics
- magnetic resonance physics
- nuclear medicine physics
- optical physics
- thermotherapy physics
- tissue measurements
- radiation biology
- review article
- books and publications
Index of content:
Electrode vibration elastography is a new shear waveimaging technique that can be used to visualize thermal ablation zones. Prior work has shown the ability of electrode vibration elastography to delineate radiofrequency ablations; however, there has been no previous study of delineation of microwave ablations or radiological–pathological correlations using multiple observers.Methods:
Radiofrequency and microwave ablations were formed inex vivo bovine livertissue. Their visualization was compared on shear wave velocity and maximum displacement images.Ablation dimensions were compared to gross pathology. Elastographic imaging and gross pathology overlap and interobserver variability were quantified using similarity measures.Results:
Elastographic imaging correlated with gross pathology. Correlation of area estimates was better in radiofrequency than in microwave ablations, with Pearson coefficients of 0.79 and 0.54 on shear wave velocity images and 0.90 and 0.70 on maximum displacement images for radiofrequency and microwave ablations, respectively. The absolute relative difference in area between elastographic imaging and gross pathology was 18.9% and 22.9% on shear wave velocity images and 16.0% and 23.1% on maximum displacement images for radiofrequency and microwave ablations, respectively.Conclusions:
Statistically significant radiological–pathological correlation was observed in this study, but correlation coefficients were lower than other modulus imaging techniques, most notably in microwave ablations. Observers provided similar delineations for most thermal ablations. These results suggest that electrode vibration elastography is capable of imaging thermal ablations, but refinement of the technique may be necessary before it can be used to monitor thermal ablation procedures clinically.
It is not appropriate to “deform” dose along with deformable image registration in adaptive radiotherapy39(2012); http://dx.doi.org/10.1118/1.4722968View Description Hide Description
- MEDICAL PHYSICS LETTERS
39(2012); http://dx.doi.org/10.1118/1.4760981View Description Hide DescriptionPurpose
: Recently reported quantitative photoacoustic tomography (PAT) has significantly expanded the utilities of PAT because it allows for recovery of tissue optical absorption coefficient which directly correlates with tissue physiological information. However, the recovery of optical absorption coefficient by the existing quantitative PAT approaches strongly depends on the accuracy of absorbed energy density distribution, and on the knowledge of accurate strength and distribution of incident light source. The purpose of this study is to develop a new algorithm for the reconstruction of optical absorption coefficient that does not depend on these initial parameters.Methods
: Here the authors propose a novel one-step reconstruction approach that can directly recover optical absorption coefficient from photoacousticmeasurements along boundary domain. The authors validate the method using simulation and phantom experiments.Results
: The authors have demonstrated experimental evidence that it is possible to directly recover optical absorption coefficient maps using boundary photoacousticmeasurements coupled with the photondiffusion equation in just one step. The authors found that the method described is able to quantitatively reconstruct absorbing objects with different sizes and optical contrast levels.Conclusions
: Compared to the authors’ previous two-step methods, the reconstruction results obtained here show that the one-step scheme can significantly improve the accuracy of absorption coefficient recovery.
39(2012); http://dx.doi.org/10.1118/1.4762288View Description Hide DescriptionPurpose:
On-board 4D cone beam CT (4DCBCT) offers respiratory phase-resolved volumetric imaging, and improves the accuracy of target localization in image guided radiation therapy. However, the clinical utility of this technique has been greatly impeded by its degraded image quality, prolonged imaging time, and increased imaging dose. The purpose of this letter is to develop a novel iterative 4DCBCT reconstruction method for improved image quality, increased imaging speed, and reduced imaging dose.Methods:
The essence of this work is to introduce the spatiotemporal tensor framelet (STF), a high-dimensional tensor generalization of the 1D framelet for 4DCBCT, to effectively take into account of highly correlated and redundant features of the patient anatomy during respiration, in a multilevel fashion with multibasis sparsifying transform. The STF-based algorithm is implemented on a GPU platform for improved computational efficiency. To evaluate the method, 4DCBCT full-fan scans were acquired within 30 s, with a gantry rotation of 200°; STF is also compared with a state-of-art reconstruction method via spatiotemporal total variation regularization.Results:
Both the simulation and experimental results demonstrate that STF-based reconstruction achieved superior image quality. The reconstruction of 20 respiratory phases took less than 10 min on an NVIDIA Tesla C2070 GPU card. The STF codes are available athttps://sites.google.com/site/spatiotemporaltensorframelet.Conclusions:
By effectively utilizing the spatiotemporal coherence of the patient anatomy among different respiratory phases in a multilevel fashion with multibasis sparsifying transform, the proposed STF method potentially enables fast and low-dose 4DCBCT with improved image quality.
- RADIATION THERAPY PHYSICS
39(2012); http://dx.doi.org/10.1118/1.4756932View Description Hide DescriptionPurpose:
Deformable registration of medical images often requires initial rigid alignment. Because of variations in the articulation of bony structures, rigid alignment can capture only limited regions of the image. We propose a method that allows us to compensate for misalignment of mobile parts, which leads to improved accuracy of deformable registration. The method is based on matching landmarks using radial basis functions (RBF) with adaptive radius.Methods:
Based on the assumption that the compactly positioned landmarks likely delineate an anatomic structure whose position needs to be corrected, the algorithm incorporates unsupervised clustering of landmarks based on their positions within the reference image. It calculates an appropriate RBF radius based on the set of pairwise distances between landmarks within the cluster. The algorithm distinguishes between clusters of different size and between clusters of spherical and elongated shape, and assigns the optimal RBF radius for each cluster in order to restrict the deformation field to the closest vicinity of the structure of interest.Results:
Experiments with synthetic images demonstrate sensitivity of registration results to the choice of the radius of RBF support. We have statistically validated the methods on a large set of pulmonary landmarks. We also tested the method on medical use cases that show that it is potentially advantageous for initial registration of images with large spatial dislocations.Conclusions:
The results of registration of CTimages demonstrate that an automated selection of the RBF radius simplifies the registration routine and improves the registration quality. The selection is based on two criteria of preserving diffeomorphism of deformation and localization of the deformation within a desired area of the image.
Unfolding linac photon spectra and incident electron energies from experimental transmission data, with direct independent validation39(2012); http://dx.doi.org/10.1118/1.4754301View Description Hide DescriptionPurpose:
In a recent computational study, an improved physics-based approach was proposed for unfolding linacphoton spectra and incident electron energies from transmission data. In this approach, energy differentiation is improved by simultaneously using transmission data for multiple attenuators and detectors, and the unfolding robustness is improved by using a four-parameter functional form to describe the photon spectrum. The purpose of the current study is to validate this approach experimentally, and to demonstrate its application on a typical clinical linac.Methods:
The validation makes use of the recent transmission measurements performed on the Vickers research linac of National Research Council Canada. For this linac, the photon spectra were previously measured using a NaI detector, and the incident electron parameters are independently known. The transmission data are for eight beams in the range 10–30 MV using thick Be, Al and Pb bremsstrahlung targets. To demonstrate the approach on a typical clinical linac, new measurements are performed on an ElektaPreciselinac for 6, 10 and 25 MV beams. The different experimental setups are modeled using EGSnrc, with the newly added photonuclear attenuation included.Results:
For the validation on the research linac, the 95% confidence bounds of the unfolded spectra fall within the noise of the NaI data. The unfolded spectra agree with theEGSnrc spectra (calculated using independently known electron parameters) with RMS energy fluence deviations of 4.5%. The accuracy of unfolding the incident electron energy is shown to be ∼3%. A transmission cutoff of only 10% is suitable for accurate unfolding, provided that the other components of the proposed approach are implemented. For the demonstration on a clinical linac, the unfolded incident electron energies and their 68% confidence bounds for the 6, 10 and 25 MV beams are 6.1 ± 0.1, 9.3 ± 0.1, and 19.3 ± 0.2 MeV, respectively. The unfolded spectra for the clinical linac agree with the EGSnrc spectra (calculated using the unfolded electron energies) with RMS energy fluence deviations of 3.7%. The corresponding measured and EGSnrc-calculated transmission data agree within 1.5%, where the typical transmission measurementuncertainty on the clinical linac is 0.4% (not including the uncertainties on the incident electron parameters).Conclusions:
The approach proposed in an earlier study for unfolding photon spectra and incident electron energies from transmission data is accurate and practical for clinical use.
2D lag and signal nonlinearity correction in an amorphous silicon EPID and their impact on pretreatment dosimetric verification39(2012); http://dx.doi.org/10.1118/1.4757582View Description Hide DescriptionPurpose:
This investigation provides measurements of signal lag and nonlinearity separately for the Varian aS500 electronic portal imaging device (EPID), and an algorithm to correct for these effects in 2D; their potential impact on intensity modulated radiation therapy (IMRT) verification is also investigated. The authors quantify lag, as a function of both delivered monitor units (MU) and time, by using a range of MUs delivered at a clinically used rate of 400 MU/min. Explicit cumulative lag curves are thus determined for a range of MUs and times between the end of irradiation and the end of image acquisition. Signal nonlinearity is also investigated as a function of total MUs delivered. The family of cumulative lag curves and signal nonlinearity are then used to determine their effects on dynamic multileaf collimator (MLC) (IMRT) deliveries, and to correct for theses effects in 2D.Methods:
Images acquired with an aS500 EPID and Varis Portal-Vision software were used to quantify detector lag and signal-nonlinearity. For the signal lag investigation, Portal-Vision'sservice monitor was used to acquire EPID images at a rate of 8 frames/s. The images were acquired during irradiation and 66 s thereafter, by inhibiting the M-holdoff-In signal of the Linac for a range of 4.5–198.5 MUs. Relative cumulative lag was calculated by integrating the EPID signal for a time after beam-off, and normalizing this to the integrated EPID signal accumulated during radiation. Signal nonlinearity was studied by acquiring 10 × 10 cm2 open-field EPID images in “integrated image” mode for a range of 2–500 MUs, and normalized to the 100 MU case. All data were incorporated into in-house written software to create a 2D correction map for these effects, using the field's MLC file and a field-specific calculated 2D “time-map,” which keeps track of the time elapsed from the last fluence delivered at each given point in the image to the end of the beam delivery.Results:
Relative cumulative lag curves reveal that the lag alone can deviate the EPID's perceived dose by as large as 6% (1 MU delivery, 60 s postirradiation). For signal nonlinearity relative to 100 MU, EPID signals per MU of 0.84 and 1.01 were observed for 2 and 500 MUs, respectively. Correction maps were applied to a 1 cm sweeping-window 14 × 14 cm2 field and clinical head-and-neck IMRT field. A mean correction of 1.028 was implemented in the head-and-neck field, which significantly reduced lag-related asymmetries in the EPID images, and restored linearity to the EPID imager's dose response. Corrections made to the sweeping-field showed good agreement with the treatment planning system-predicted field, yielding an average percent difference of 0.05% ± 0.91%, compared to the −1.32% ± 1.02% before corrections, or 1.75% ± 1.04% when only a signal nonlinearity correction is made.Conclusions:
Lag and signal-nonlinearity have been quantified for an aS500 EPID imager, and an effective 2D correction method has been developed which effectively removes nonlinearity and lag effects. Both of these effects were shown to negatively impact IMRT verifications. Especially fields that involve prolonged irradiation and small overall MUs should be corrected for in 2D.
39(2012); http://dx.doi.org/10.1118/1.4757616View Description Hide DescriptionPurpose:
To investigate the dosimetric feasibility of accelerated partial breast irradiation (APBI) using CyberKnife.Methods:
Fourteen previously treated patients with early-stage breast cancer were selected for a retrospective study. Six of these patients had been treated to 38.5 Gy in 10 fractions in a phase III accelerated partial breast trial and the rest of the patients were treated to 50.4 Gy in 28 fractions. In this planning study, the guidelines in the protocol for the phase III partial breast trial were followed for organ delineation and CyberKnife planning. The achievable dosimetric parameters from all CyberKnife plans were compared to Intensity-modulated radiation therapy (IMRT) and 3D-CRT methods. The reproducibility of the dose delivery with and without respiratory motion was assessed through delivering a patient plan to a breast phantom. Different dose calculation algorithms were also compared between ray tracing and Monte Carlo.Results:
For all the patients in the study, the dosimetric parameters met the guidelines from the NSABP B39/RTOG 0413 protocol strictly. The mean PTV volume covered by 100% of the prescription dose was 95.7 ± 0.7% (94.7%–97.1%). The mean maximal dose was 104 ± 2% of the prescription dose. The mean V50% and mean V100% to the ipsilateral normal breast were 23.1 ± 11.6% and 9.0 ± 5.8%, respectively. The conformity index of all plans was 1.14 ± 0.04. The maximum dose to the contralateral breast varied from 1.3 cGy to 111 cGy. The mean V5% and mean V30% to the contralateral and ipsilateral lungs were 1.0 ± 1.6% and 1.3 ± 1.2%, respectively. In our study, the mean V5% to the heart was 0.2 ± 0.5% for right-sided tumors and 9.4 ± 10.1% for left-sided tumors. Compared with IMRT and 3D-CRT planning, the PTV coverage from CyberKnife planning was the highest, and the ratio of V20% to V100% of the breast from CyberKnife planning was the smallest. The heart and lung doses were similar in all the techniques except that the V5% for the lung and heart in CyberKnife planning was slightly higher.Conclusions:
The dosimetric feasibility of APBI using CyberKnife was investigated in this retrospective study. All the dosimetric parameters strictly met the guidelines from the NSABP B39/RTOG 0413 protocol. With advanced real-time tracking capability, CyberKnife should provide better target coverage and spare nearby critical organs for APBI treatment.
39(2012); http://dx.doi.org/10.1118/1.4757916View Description Hide DescriptionPurpose:
A major concern for lungintensity modulated radiation therapydelivery is the deviation of actually delivereddose distribution from the planned one due to simultaneous movements of multileaf collimator(MLC) leaves and tumor. For gated lung stereotactic body radiotherapytreatment(SBRT), the situation becomes even more complicated because of SBRT's characteristics such as fewer fractions, smaller target volume, higher dose rate, and extended fractional treatment time. The purpose of this work is to investigate the dosimetric effect of intrafraction tumor motion during gated lungSBRTdelivery by reconstructing the delivereddose distribution with real-time tumor motion considered.Methods:
The tumor motion data were retrieved from six lung patients. Each of them received three fractions of stereotactic radiotherapytreatments with Cyberknife Synchrony (Accuray, Sunnyvale, CA). Phase gating through an external surrogate was simulated with a gating window of 5 mm. The resulting residual tumor motion curves during gating (beam-on) were retrieved. Planning target volume (PTV) was defined as physician-contoured clinical target volume (CTV) surrounded by an isotropic 5 mm margin. Each patient was prescribed with 60 Gy/3 fractions. The authors developed an algorithm to reconstruct the delivereddose with tumor motion. The DMLC segments, mainly leaf position and segment weighting factor, were recalculated according to the probability density function of tumor motion curve. The new DMLC sequence file was imported back to treatment planning system to reconstruct the dose distribution.Results:
Half of the patients in the study group experienced PTV D95% deviation up to 26% for fractional dose and 14% for total dose. CTV mean dose dropped by 1% with tumor motion. Although CTV is almost covered by prescribed dose with 5 mm margin, qualitative comparison on the dose distributions reveals that CTV is on the verge of underdose. The discrepancy happens due to tumor excursion outside of the gating window, which, for our study group, is mainly caused by baseline shift, i.e., the change in general trend of the motion curve during extended period of treatment time.Conclusions:
The dose deviation in PTV and CTV due to target motion is not always negligible in gated SBRT. Although CTVs are covered sufficiently with prescribed dose in most cases, some are on the verge of underdose due to large tumor excursion caused by factors such as baseline shift.
39(2012); http://dx.doi.org/10.1118/1.4757918View Description Hide DescriptionPurpose:
Arteriovenous malformations are often treated with a combination of embolization and stereotactic radiosurgery. Concern has been expressed in the past regarding the dosimetric properties of materials used in embolization and the effects that the introduction of these materials into the brain may have on the quality of the radiosurgery plan. To quantify these effects, the authors have taken large volumes of Onyx 34 and Onyx 18 (ethylene-vinyl alcohol copolymer doped with tantalum) and measured the attenuation and interface effects of these embolization materials.Methods:
The manufacturer provided large cured volumes (∼28 cc) of both Onyxmaterials. These samples were 8.5 cm in diameter with a nominal thickness of 5 mm. The samples were placed on a block tray above a stack of solid water with an Attix chamber at a depth of 5 cm within the stack. The Attix chamber was used to measure the attenuation. These measurements were made for both 6 and 16 MV beams. Placing the sample directly on the solid water stack and varying the thickness of solid water between the sample and the Attix chamber measured the interface effects. The computed tomography(CT) numbers for bulk material were measured in a phantom using a wide bore CT scanner.Results:
The transmission through the Onyxmaterials relative to solid water was approximately 98% and 97% for 16 and 6 MV beams, respectively. The interface effect shows an enhancement of approximately 2% and 1% downstream for 16 and 6 MV beams. CT numbers of approximately 2600–3000 were measured for both materials, which corresponded to an apparent relative electron density (RED) to water of approximately 2.7–2.9 if calculated from the commissioning data of the CT scanner.Conclusions:
We performed direct measurements of attenuation and interface effects of Onyx 34 and Onyx 18 embolization materials with large samples. The introduction of embolization materials affects the dose distribution of a MV therapeutic beam, but should be of negligible consequence for effective thicknesses of less than 8 mm. The measured interface effects are also small, particularly at 6 MV. Large areas of high-density artifacts and low-density artifacts can cause errors in dose calculations and need to be identified and resolved during planning.
A Monte Carlo based formalism to identify potential locations at high risk of tumor recurrence with a numerical model for glioblastoma multiforme39(2012); http://dx.doi.org/10.1118/1.4757972View Description Hide DescriptionPurpose:
The strategy currently used to treat glioblastoma multiforme (GBM) patients, which mostly relies on population-based failure patterns, does not consider the important variability in such patterns reported in the literature. As part of the multidisciplinary efforts being made to develop personalized therapeutic approaches, numerical models of tumor growth and treatment are increasingly being used by different groups around the world. In this study, a new formalism relying on the proliferation-invasion model is developed to identify potential locations of GBM recurrences. The authors assess the sensitivity of the location of potential tumor recurrences to the input parameter values predicted for a given patient by varying those values using a Monte-Carlo based approach. Our approach is designed to be prospective in the sense that it relies on patient-specific imaging data that can be gathered in one single preradiotherapy imaging session.Methods:
The authors modeled the infiltration paths of glial cells using patient-specific diffusiontensorimaging (DTI) data. Nine GBM patients with preradiotherapy DTI data are considered in this study. The possible locations of tumor recurrences are determined by randomly selecting many ensembles of values for each of the growth and radiobiological parameters in the GBM growth model. A novel concept, the occurrence probability (OP), is introduced to assess the sensitivity of potential tumor recurrence locations to the input parameter values. For a given patient, the OP map is derived from a superposition of all potential tumor recurrence locations obtained with all sets of parameter values.Results:
For eight out of nine of patients, the authors have identified a statistically significant region where the OP is above 50%. For two patients, these high risk regions are found to be located at a distance greater than 3.9 cm from the border of the gross tumor volume highlighting the inaccuracy of current margins for some patients. The exact location and size of these volumes with OP > 50 % are, however, sensitive to the numberN of ensembles of parameter values for N ≲ 400. On the other hand, the authors have identified for each patient a threshold OP, the OPT, which defines a volume that converges more rapidly with increasing N. The OPT for each patient varies between 20% and 40%. The volume defined by OP > OPT may be an adequate candidate to define a personalized margin for radiotherapy treatment planning of GBM patients.Conclusions:
A new Monte-Carlo based formalism was described and used to assess the variability of sites of potential recurrence predicted by the proliferation-invasion model to input parameter values. The authors have shown that high risk areas could be consistently identified with a limited number of sets (N ≲ 400) of randomly chosen parameter values. A major strength of this formalism is its potential prospective nature. Although a validation of the accuracy of the model-predicted tumor recurrence location still remains to be done, our method is potentially applicable to orient patient-specific definition of margins.
Multicriteria optimization for volumetric-modulated arc therapy by decomposition into a fluence-based relaxation and a segment weight-based restriction39(2012); http://dx.doi.org/10.1118/1.4754652View Description Hide DescriptionPurpose:
To develop a method for inverse volumetric-modulated arc therapy (VMAT) planning that combines multicriteria optimization (MCO) with direct machine parameter optimization. The ultimate goal is to provide an efficient and intuitive method for generating high quality VMAT plans.Methods:
Multicriteria radiation therapy treatment planning amounts to approximating the relevant treatment options by a discrete set of plans, and selecting the combination thereof that strikes the best possible balance between conflicting objectives. This approach is applied to two decompositions of the inverse VMAT planning problem: a fluence-based relaxation considered at a coarsened gantry angle spacing and under a regularizing penalty on fluence modulation, and a segment weight-based restriction in a neighborhood of the solution to the relaxed problem. The two considered variable domains are interconnected by direct machine parameter optimization toward reproducing the dose-volume histogram of the fluence-based solution.Results:
The dose distribution quality of plans generated by the proposed MCO method was assessed by direct comparison with benchmark plans generated by a conventional VMAT planning method. The results for four patient cases (prostate, pancreas, lung, and head and neck) are highly comparable between the MCO plans and the benchmark plans: Discrepancies between studied dose-volume statistics for organs at risk were—with the exception of the kidneys of the pancreas case—within 1 Gy or 1 percentage point. Target coverage of the MCO plans was comparable with that of the benchmark plans, but with a small tendency toward a shift from conformity to homogeneity.Conclusions:
MCO allows tradeoffs between conflicting objectives encountered in VMAT planning to be explored in an interactive manner through search over a continuous representation of the relevant treatment options. Treatment plans selected from such a representation are of comparable dose distribution quality to conventionally optimized VMAT plans.
39(2012); http://dx.doi.org/10.1118/1.4757578View Description Hide DescriptionPurpose:
To develop a flexible pencil beam algorithm for helium ion beam therapy. Dose distributions were calculated using the newly developed pencil beam algorithm and validated using Monte Carlo(MC) methods.Methods:
The algorithm was based on the established theory of fluence weighted elemental pencil beam (PB) kernels. Using a new real-time splitting approach, a minimization routine selects the optimal shape for each sub-beam. Dose depositions along the beam path were determined using a look-up table (LUT). Data for LUT generation were derived from MC simulations in water using GATE 6.1. For materials other than water, dose depositions were calculated by the algorithm using water-equivalent depth scaling. Lateral beam spreading caused by multiple scattering has been accounted for by implementing a non-local scattering formula developed by Gottschalk. A new nuclear correction was modelled using a Voigt function and implemented by a LUT approach. Validation simulations have been performed using a phantom filled with homogeneous materials or heterogeneous slabs of up to 3 cm. The beams were incident perpendicular to the phantoms surface with initial particle energies ranging from 50 to 250 MeV/A with a total number of 107 ions per beam. For comparison a special evaluation software was developed calculating the gamma indices for dose distributions.Results:
In homogeneous phantoms, maximum range deviations between PB and MC of less than 1.1% and differences in the width of the distal energy falloff of the Bragg-Peak from 80% to 20% of less than 0.1 mm were found. Heterogeneous phantoms using layered slabs satisfied a γ-index criterion of 2%/2mm of the local value except for some single voxels. For more complex phantoms using laterally arranged bone-air slabs, the γ-index criterion was exceeded in some areas giving a maximum γ-index of 1.75 and 4.9% of the voxels showed γ-index values larger than one. The calculation precision of the presented algorithm was considered to be sufficient for clinical practice. Although only data for helium beams was presented, the performance of the pencil beam algorithm for protonbeams was comparable.Conclusions:
The pencil beam algorithm developed for helium ions presents a suitable tool for dose calculations. Its calculation speed was evaluated to be similar to other published pencil beam algorithms. The flexible design allows easy customization of measured depth-dose distributions and use of varying beam profiles, thus making it a promising candidate for integration into future treatment planning systems. Current work in progress deals with RBE effects of helium ions to complete the model.
A real-time in vivo dosimetric verification method for high-dose rate intracavitary brachytherapy of nasopharyngeal carcinoma39(2012); http://dx.doi.org/10.1118/1.4758067View Description Hide DescriptionPurpose:
A real-timein vivodosimetric verification method using metal-oxide-semiconductor field effect transistor (MOSFET)dosimeters has been developed for patient dosimetry in high-dose rate (HDR) intracavitary brachytherapy of nasopharyngeal carcinoma (NPC).Methods:
The necessary calibration and correction factors for MOSFET measurements in192Iridium source were determined in a water phantom. With the detector placed inside a custom-made nasopharyngeal applicator, the actual dosedelivered to the tumor was measured in vivo and compared to the calculated values using a commercial brachytherapy planning system.Results:
Five MOSFETs were independently calibrated with the HDR source, yielding calibration factors of 0.48 ± 0.007 cGy/mV. The maximum sensitivity variation was no more than 7% in the clinically relevant distance range of 1–5 cm from the source. A total of 70in vivo measurements in 11 NPC patients demonstrated good agreement with the treatment planning. The mean differences between the planned and the actually delivereddose within a single treatment fraction were −0.1% ± 3.8% and −0.1% ± 3.7%, respectively, for right and left side assessments. The maximum dose deviation was less than 8.5%.Conclusions:
In vivo measurement using the real-time MOSFETdosimetry system is possible to evaluate the actual dose to the tumor received by the patient during a treatment fraction and thus can offer another line of security to detect and prevent large errors.
39(2012); http://dx.doi.org/10.1118/1.4754802View Description Hide DescriptionPurpose:
To accurately quantify the local difference between two contour surfaces in two- or three-dimensional space, a new, robust point-to-surface distance measure is developed.Methods:
To evaluate and visualize the local surface differences, point-to-surface distance measures have been utilized. However, previously well-known point-to-surface distance measures have critical shortfalls. Previous distance measures termed “normal distance (ND),” “radial distance,” or “minimum distance (MD)” can report erroneous results at certain points where the surfaces under comparison meet certain conditions. These skewed results are due to the monodirectional characteristics of these methods. ComGrad distance was also proposed to overcome asymmetric characteristics of previous point-to-surface distance measures, but their critical incapability of dealing with a fold or concave contours. In this regard, a new distance measure termed the bidirectional local distance (BLD) is proposed which minimizes errors of the previous methods by taking into account the bidirectional characteristics with the forward and backward directions. BLD measure works through three steps which calculate the maximum value between the forward minimum distance (FMinD) and the backward maximum distance (BMaxD) at each point. The first step calculates the FMinD as the minimum distance to the test surface from a point,p ref on the reference surface. The second step involves calculating the minimum distances at every point on the test surface to the reference surface. During the last step, the BMaxD is calculated as the maximum distance among the minimum distances found at p ref on the reference surface. Tests are performed on two- and three-dimensional artificial contour sets in comparison to MD and ND measure techniques. Three-dimensional tests performed on actual liver and head-and-neck cancer patients.Results:
The proposed BLD measure provides local distances between segmentations, even in situations where ND, MD, or ComGrad measures fail. In particular, the standard deviation measure is not distorted at certain geometries where ND, MD, and ComGrad measures report skewed results.Conclusions:
The proposed measure provides more reliable statistics on contour comparisons. From the statistics, specific local and global distances can be extracted. Bidirectional local distance is a reliable distance measure in comparing two- or three-dimensional organ segmentations.
Graphical representation of the effects on tumor and OAR for determining the appropriate fractionation regimen in radiation therapy planning39(2012); http://dx.doi.org/10.1118/1.4757580View Description Hide DescriptionPurpose:
The authors propose a graphical representation of the relation between the effect on the tumor and the damage effect on an organ at risk (OAR) against the irradiation dose, as an aid for choosing an appropriate fractionation regimen.Methods:
The graphical relation is depicted by the radiation effect on the tumorE 1 versus that on an OAR E 0. By observing the features of the E 1 vs E 0 relation curve, i.e., convex or concave shape, one can judge whether multifractionation is better or not. This method is applied to the linear-quadratic model (with α and β parameters) as an example. Further, the method is extended to the general case for nonuniform dose distribution to the OAR, which is frequently seen in clinical situations.Results:
The criterion for selecting multi- or hypofractionation is based on the relation between the dose for the OAR and theα/β ratio of the OAR to the tumor. It is also shown that the graphical relation enables us to estimate the final effect after multifractionated treatment by plotting a tangent line on the curve.Conclusions:
The graphical representation method is of use for improving planning in radiotherapy by determining the effective fractionation scheme.
39(2012); http://dx.doi.org/10.1118/1.4758060View Description Hide DescriptionPurpose:
While Monte Carlo particle transport has proven useful in many areas (treatment head design, dose calculation, shielding design, and imaging studies) and has been particularly important for proton therapy (due to the conformal dose distributions and a finite beam range in the patient), the available general purpose Monte Carlo codes in proton therapy have been overly complex for most clinical medical physicists. The learning process has large costs not only in time but also in reliability. To address this issue, we developed an innovative protonMonte Carlo platform and tested the tool in a variety of proton therapy applications.Methods:
Our approach was to take one of the already-established general purpose Monte Carlo codes and wrap and extend it to create a specialized user-friendly tool for proton therapy. The resulting tool, TOol for PArticle Simulation (TOPAS), should make Monte Carlo simulation more readily available for research and clinical physicists. TOPAS can model a passive scattering or scanning beam treatment head, model a patient geometry based on computed tomography(CT)images, score dose, fluence, etc., save and restart a phase space, provides advanced graphics, and is fully four-dimensional (4D) to handle variations in beam delivery and patient geometry during treatment. A custom-designed TOPAS parameter control system was placed at the heart of the code to meet requirements for ease of use, reliability, and repeatability without sacrificing flexibility.Results:
We built and tested the TOPAS code. We have shown that the TOPAS parameter system provides easy yet flexible control over all key simulation areas such as geometry setup, particle source setup, scoring setup, etc. Through design consistency, we have insured that user experience gained in configuring one component, scorer or filter applies equally well to configuring any other component, scorer or filter. We have incorporated key lessons from safety management, proactively removing possible sources of user error such as line-ordering mistakes. We have modeled proton therapytreatment examples including the UCSF eye treatment head, the MGH stereotactic alignment in radiosurgery treatment head and the MGH gantry treatment heads in passive scattering and scanning modes, and we have demonstrated dose calculation based on patient-specific CT data. Initial validation results show agreement with measured data and demonstrate the capabilities of TOPAS in simulating beam delivery in 3D and 4D.Conclusions:
We have demonstrated TOPAS accuracy and usability in a variety of proton therapy setups. As we are preparing to make this tool freely available for researchers in medical physics, we anticipate widespread use of this tool in the growing proton therapy community.
Design and dosimetric characteristics of a new endocavitary contact radiotherapy system using an electronic brachytherapy sourcea)39(2012); http://dx.doi.org/10.1118/1.4757915View Description Hide DescriptionPurpose:
To present design aspects and acceptance tests performed for clinical implementation of electronic brachytherapytreatment of early stage rectal adenocarcinoma. A dosimetric comparison is made between the historically used Philips RT-50 unit and the newly developed Axxent® Model S700 electronic brachytherapysource manufactured by Xoft (iCad, Inc.).Methods:
Two proctoscope cones were manufactured by ElectroSurgical Instruments (ESI). Two custom surface applicators were manufactured by Xoft and were designed to fit and interlock with the proctoscope cones from ESI. Dose rates, half value layers (HVL), and percentage depth dose (PDD) measurements were made with the Xoft system and compared to historical RT-50 data. A description of the patient treatment approach and exposure rates during the procedure is also provided.Results:
The electronic brachytherapy system has a lower surfacedose rate than the RT-50. The dose rate to water on the surface from the Xoft system is approximately 2.1 Gy/min while the RT-50 is 10–12 Gy/min. However, treatment times with Xoft are still reasonable. The HVLs and PDDs between the two systems were comparable resulting in similar doses to the target and to regions beyond the target. The exposure rate levels around a patient treatment were acceptable. The standard uncertainty in the dose rate to water on the surface is approximately ±5.2%.Conclusions:
The Philips RT-50 unit is an out-of-date radiotherapy machine that is no longer manufactured with limited replacement parts. The use of a custom-designed proctoscope and Xoft surface applicators allows delivery of a well-established treatment with the ease of a modern radiotherapy device. While the dose rate is lower with the use of Xoft, the treatment times are still reasonable. Additionally, personnel may stand farther away from the Xoft radiationsource, thus potentially reducing radiation exposure to the operator and other personnel.
Quantitative analysis of the factors which affect the interpatient organ-at-risk dose sparing variation in IMRT plans39(2012); http://dx.doi.org/10.1118/1.4757927View Description Hide DescriptionPurpose:
The authors present an evidence-based approach to quantify the effects of an array of patient anatomical features of the planning target volumes (PTVs) and organs-at-risk (OARs) and their spatial relationships on the interpatient OAR dose sparing variation in intensity modulated radiation therapy(IMRT) plans by learning from a database of high-quality prior plans.Methods:
The authors formulized the dependence of OAR dose volume histograms (DVHs) on patient anatomical factors into feature models which were learned from prior plans by a stepwise multiple regression method. IMRT plans for 64 prostate, 82 head-and-neck (HN) treatments were used to train the models. Two major groups of anatomical features were considered in this study: the volumetric information and the spatial information. The geometry of OARs relative to PTV is represented by the distance-to-target histogram, DTH. Important anatomical and dosimetric features were extracted from DTH and DVH by principal component analysis. The final models were tested by additional 24 prostate and 24 HN plans.Results:
Significant patient anatomical factors contributing to OAR dose sparing in prostate and HN IMRT plans have been analyzed and identified. They are: the median distance between OAR and PTV, the portion of OAR volume within an OAR specific distance range, and the volumetric factors: the fraction of OAR volume which overlaps with PTV and the portion of OAR volume outside the primary treatment field. Overall, the determination coefficients R2 for predicting the first principal component score (PCS1) of the OAR DVH by the above factors are above 0.68 for all the OARs and they are more than 0.53 for predicting the second principal component score (PCS2) of the OAR DVHs except brainstem and spinal cord. Thus, the above set of anatomical features combined has captured significant portions of the DVH variations for the OARs in prostate and HN plans. To test how well these features capture the interpatient organdose sparing variations in general, the DVHs and specific dose-volume indices calculated from the regression models were compared with the actual DVHs and dose-volume indices from each patient's plan in the validation dataset. The dose-volume indices compared were V99%, V85%, and V50% for bladder and rectum in prostate plans and parotids median dose in HN plans. The authors found that for the bladder and rectum models, 17 out of 24 plans (71%) were within 6% OAR volume error and 21 plans (85%) were within 10% error; For the parotids model, the median dose values for 30 parotids out of 48 (63%) were within 6% prescription dose error and the values in 40 parotids (83%) were within 10% error.Conclusions:
Quantitative analysis of patient anatomical features and their correlation with OAR dose sparing has identified a number of important factors that explain significant amount of interpatient DVH variations in OARs. These factors can be incorporated into evidence-based learning models as effective features to provide patient-specific OAR dose sparing goals.
39(2012); http://dx.doi.org/10.1118/1.4761865View Description Hide DescriptionPurpose:
Quality assurance is an essential component of accurate and safe radiotherapy delivery, and should include measurements which are independent of manufacturer-provided calibration. However, the physical and dosimetric properties of the INTRABEAM compact mobile 50 kV x-ray source are different from conventional kilovoltage therapy units and few reports describe methods for independent checks, frequencies, or tolerances for quality assurance tests.Methods:
Based on the available evidence and local experience, methods are described for determination of the key dosimetric parameters: beam quality, output, isotropy, and depth doses. Internal system checks are also described, along with measurements of long-term stability.Results:
A small volume parallel plate ionization chamber in a liquid water tank is the gold standard for measurements with this unit, but solid water-equivalent materials, thermoluminescent dosimeters and radiochromic film can all be used as practical alternatives with an accuracy of 5%–10%. The main cause of measurementuncertainty is positioning of the detector in the steep dose gradient, but energy dependence should also be considered.Conclusions:
A quality assurance schedule with suggested tolerances is proposed, which includes both internal tests, before each treatment and on a monthly basis, and independent tests every year or after servicing or recalibration.
39(2012); http://dx.doi.org/10.1118/1.4762684View Description Hide DescriptionPurpose:
To quantify systematically the effect on accuracy of discretizing gantry rotation during the dose calculation process of TomoTherapy treatments.Methods:
Up to version 4.0.x included, TomoTherapy treatment planning system (TPS) approximates gantry rotation by computing dose from 51 discrete angles corresponding to the center of the projections used to control the binary multileaf collimator. Potential effects on dose computation accuracy for off-axis targets and low modulation factors have been shown previously for a few treatment configurations. In versions 4.1.x and later, TomoTherapy oversamples the projections to better account for gantry rotation, but only during full scatter optimization and final calculation (i.e., not during optimization in “beamlet” mode). The effect on accuracy of changing the number of angles was quantified with the following framework: (1) predict the impact of the discretization of gantry rotation for various modulation factors, target sizes, and off-axis positions using a simplified analytical algorithm; (2) perform regular quality assurance using measurements with EDR2 radiographic films; (3) isolating the effect of changing the number of discretized angles only (51, 153, and 459) using a previously validated Monte Carlo model (TomoPen). The diameters of the targets were 2, 3, and 5 cm; off-axis central positions of target volumes were 5, 10 and 15, and 17 cm (when accepted by the treatment unit); planned modulation factors were 1.3 and 2.0.Results:
For extreme configurations (3 cm tumor, 1.3 modulation factor, 15 cm off-axis position), effects on dose distributions were significant with 89.3% and 95.4% of the points passing gamma tests with 2%/2 mm and 3%/3 mm criteria, respectively, for TPS software version 4.0.x (51 gantry angles). The passing rate was 100% for both gamma criteria for the 4.1.x version (153 gantry angles). Those differences could be attributed almost completely to gantry motion discretization using TomoPen. Using 51 gantry angles for dose computation, TomoPen reproduced within statistical uncertainties (<1% standard deviation) dose distributions computed with version 4.0.x. Using 153 and 459 gantry angles, TomoPen reproduced within statistical uncertainties measurements and dose distributions computed with version 4.1.x.Conclusions:
When low modulation factors and significant off-axis positions are used, accounting for gantry rotation during dose computation using at least 153 gantry angles is required to ensure optimal accuracy.