Search NASA⌕ Search

Engineering topics

Franich, R. D.

Publications and source records attributed to Franich, R. D..

Source position verification and dosimetry in HDR brachytherapy using an EPID

Purpose: Accurate treatment delivery in high dose rate (HDR) brachytherapy requires correct source dwell positions and dwell times to be administered relative to each other and to the surrounding anatomy. Treatment delivery inaccuracies predominantly occur for two reasons: (i) anatomical movement or (ii) as a result of human errors that are usually related to incorrect implementation of the planned treatment. Electronic portal imaging devices (EPIDs) were originally developed for patient position verification in external beam radiotherapy and their application has been extended to provide dosimetric information. The authors have characterized the response of an EPID for use with an {sup 192}Ir brachytherapy source to demonstrate its use as a verification device, providing both source position and dosimetric information.Methods: Characterization of the EPID response using an {sup 192}Ir brachytherapy source included investigations of reproducibility, linearity with dose rate, photon energy dependence, and charge build-up effects associated with exposure time and image acquisition time. Source position resolution in three dimensions was determined. To illustrate treatment verification, a simple treatment plan was delivered to a phantom and the measured EPID dose distribution compared with the planned dose.Results: The mean absolute source position error in the plane parallel to the EPID, for dwells measured at 50, 100, and 150 mm source to detector distances (SDD), was determined to be 0.26 mm. The resolution of the z coordinate (perpendicular distance from detector plane) is SDD dependent with 95% confidence intervals of ±0.1, ±0.5, and ±2.0 mm at SDDs of 50, 100, and 150 mm, respectively. The response of the EPID is highly linear to dose rate. The EPID exhibits an over-response to low energy incident photons and this nonlinearity is incorporated into the dose calibration procedure. A distance (spectral) dependent dose rate calibration procedure has been developed. The difference between measured and planned dose is less than 2% for 98.0% of pixels in a two-dimensional plane at an SDD of 100 mm.Conclusions: Our application of EPID dosimetry to HDR brachytherapy provides a quality assurance measure of the geometrical distribution of the delivered dose as well as the source positions, which is not possible with any current HDR brachytherapy verification system.

61 RADIATION PROTECTION AND DOSIMETRY↗

Performance of 12 DIR algorithms in low-contrast regions for mass and density conserving deformation

Purpose: Deformable image registration (DIR) has become a key tool for adaptive radiotherapy to account for inter- and intrafraction organ deformation. Of contemporary interest, the application to deformable dose accumulation requires accurate deformation even in low contrast regions where dose gradients may exist within near-uniform tissues. One expects high-contrast features to generally be deformed more accurately by DIR algorithms. The authors systematically assess the accuracy of 12 DIR algorithms and quantitatively examine, in particular, low-contrast regions, where accuracy has not previously been established.Methods: This work investigates DIR algorithms in three dimensions using deformable gel (DEFGEL) [U. J. Yeo, M. L. Taylor, L. Dunn, R. L. Smith, T. Kron, and R. D. Franich, “A novel methodology for 3D deformable dosimetry,” Med. Phys. 39, 2203–2213 (2012)], for application to mass- and density-conserving deformations. CT images of DEFGEL phantoms with 16 fiducial markers (FMs) implanted were acquired in deformed and undeformed states for three different representative deformation geometries. Nonrigid image registration was performed using 12 common algorithms in the public domain. The optimum parameter setup was identified for each algorithm and each was tested for deformation accuracy in three scenarios: (I) original images of the DEFGEL with 16 FMs; (II) images with eight of the FMs mathematically erased; and (III) images with all FMs mathematically erased. The deformation vector fields obtained for scenarios II and III were then applied to the original images containing all 16 FMs. The locations of the FMs estimated by the algorithms were compared to actual locations determined by CT imaging. The accuracy of the algorithms was assessed by evaluation of three-dimensional vectors between true marker locations and predicted marker locations.Results: The mean magnitude of 16 error vectors per sample ranged from 0.3 to 3.7, 1.0 to 6.3, and 1.3 to 7.5 mm across algorithms for scenarios I to III, respectively. The greatest accuracy was exhibited by the original Horn and Schunck optical flow algorithm. In this case, for scenario III (erased FMs not contributing to driving the DIR calculation), the mean error was half that of the modified demons algorithm (which exhibited the greatest error), across all deformations. Some algorithms failed to reproduce the geometry at all, while others accurately deformed high contrast features but not low-contrast regions—indicating poor interpolation between landmarks.Conclusions: The accuracy of DIR algorithms was quantitatively evaluated using a tissue equivalent, mass, and density conserving DEFGEL phantom. For the model studied, optical flow algorithms performed better than demons algorithms, with the original Horn and Schunck performing best. The degree of error is influenced more by the magnitude of displacement than the geometric complexity of the deformation. As might be expected, deformation is estimated less accurately for low-contrast regions than for high-contrast features, and the method presented here allows quantitative analysis of the differences. The evaluation of registration accuracy through observation of the same high contrast features that drive the DIR calculation is shown to be circular and hence misleading.

61 RADIATION PROTECTION AND DOSIMETRY↗

A phantom for testing of 4D-CT for radiotherapy of small lesions

Purpose: The use of time-resolved four-dimensional computed tomography (4D-CT) in radiotherapy requires strict quality assurance to ensure the accuracy of motion management protocols. The aim of this work was to design and test a phantom capable of large amplitude motion for use in 4D-CT, with particular interest in small lesions typical for stereotactic body radiotherapy. Methods: The phantom of 'see-saw' design is light weight, capable of including various sample materials and compatible with several surrogate marker signal acquisition systems. It is constructed of polymethylmethacrylate (Perspex) and its movement is controlled via a dc motor and drive wheel. It was tested using two CT scanners with different 4D acquisition methods: the Philips Brilliance Big Bore CT (helical scan, pressure belt) and a General Electric Discovery STE PET/CT (axial scan, infrared marker). Amplitudes ranging from 1.5 to 6.0 cm and frequencies of up to 40 cycles per minute were used to study the effect of motion on image quality. Maximum intensity projections (MIPs), as well as average intensity projections (AIPs) of moving objects were investigated and their quality dependence on the number of phase reconstruction bins assessed. Results: CT number discrepancies between moving and stationary objects were found to have no systematic dependence on amplitude, frequency, or specific interphase variability. MIP-delineated amplitudes of motion were found to match physical phantom amplitudes to within 2 mm for all motion scenarios tested. Objects undergoing large amplitude motions (>3.0 cm) were shown to cause artefacts in MIP and AIP projections when ten phase bins were assigned. This problem can be mitigated by increasing the number of phase bins in a 4D-CT scan. Conclusions: The phantom was found to be a suitable tool for evaluating the image quality of 4D-CT motion management technology, as well as providing a quality assurance tool for intercenter/intervendor testing of commercial 4D-CT systems. When imaging objects with large amplitudes, the completeness criterion described here indicates the number of phase bins required to prevent missing data in MIPs and AIPs. This is most relevant for small lesions undergoing large motions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Is it sensible to 'deform' dose? 3D experimental validation of dose-warping

Purpose: Strategies for dose accumulation in deforming anatomy are of interest in radiotherapy. Algorithms exist for the deformation of dose based on patient image sets, though these are sometimes contentious because not all such image calculations are constrained by physical laws. While tumor and organ motion has been a key area of study for a considerable amount of time, deformation is of increasing interest. In this work, we demonstrate a full 3D experimental validation of results from a range of dose deformation algorithms available in the public domain. Methods: We recently developed the first tissue-equivalent, full 3D deformable dosimetric phantom-'DEFGEL.' To assess the accuracy of dose-warping based on deformable image registration (DIR), we have measured doses in undeformed and deformed states of the DEFGEL dosimeter and compared these to planned doses and warped doses. In this way we have directly evaluated the accuracy of dose-warping calculations for 11 different algorithms. We have done this for a range of stereotactic irradiation schemes and types and magnitudes of deformation. Results: The original Horn and Schunck algorithm is shown to be the best performing of the 11 algorithms trialled. Comparing measured and dose-warped calculations for this method, it is found that for a 10 Multiplication-Sign 10 mm{sup 2} square field, {gamma}{sub 3%/3mm}= 99.9%; for a 20 Multiplication-Sign 20 mm{sup 2} cross-shaped field, {gamma}{sub 3%/3mm}= 99.1%; and for a multiple dynamic arc (0.413 cm{sup 3} PTV) treatment adapted from a patient treatment plan, {gamma}{sub 3%/3mm}= 95%. In each case, the agreement is comparable to-but consistently {approx}1% less than-comparison between measured and calculated (planned) dose distributions in the absence of deformation. The magnitude of the deformation, as measured by the largest displacement experienced by any voxel in the volume, has the greatest influence on the accuracy of the warped dose distribution. Considering the square field case, the smallest deformation ({approx}9 mm) yields agreement of {gamma}{sub 3%/3mm}= 99.9%, while the most significant deformation ({approx}20 mm) yields agreement of {gamma}{sub 3%/3mm}= 96.7%. Conclusions: We have confirmed that, for a range of mass and density conserving deformations representative of those observable in anatomical targets, DIR-based dose-warping can yield accurate predictions of the dose distribution. Substantial differences can be seen between the results of different algorithms indicating that DIR performance should be scrutinized before application todose-warping. We have demonstrated that the DEFGEL deformable dosimeter can be used to evaluate DIR performance and the accuracy of dose-warping results by direct measurement.

60 APPLIED LIFE SCIENCES↗

A novel methodology for 3D deformable dosimetry

Purpose: Interfraction and intrafraction variation in anatomic structures is a significant challenge in contemporary radiotherapy. The objective of this work is to develop a novel tool for deformable structure dosimetry, using a tissue-equivalent deformable gel dosimeter that can reproducibly simulate targets subject to deformation. This will enable direct measurement of integrated doses delivered in different deformation states, and the verification of dose deforming algorithms. Methods: A modified version of the nPAG polymer gel has been used as a deformable 3D dosimeter and phantom to investigate doses delivered to deforming tissue-equivalent geometry. The deformable gel (DEFGEL) dosimeter/phantom is comprised of polymer gel in a latex membrane, moulded (in this case) into a cylindrical geometry, and deformed with an acrylic compressor. Fifteen aluminium fiducial markers (FM) were implanted into DEFGEL phantoms and the reproducibility of deformation was determined via multiple computed tomography (CT) scans in deformed and nondeformed states before and after multiple (up to 150) deformations. Dose was delivered to the DEFGEL phantom in three arrangements: (i) without deformation, (ii) with deformation, and (iii) cumulative exposures with and without deformation, i.e., dose integration. Irradiations included both square field and a stereotactic multiple dynamic arc treatment adapted from a patient plan. Doses delivered to the DEFGEL phantom were read out using cone beam optical CT. Results: Reproducibility was verified by observation of interscan shifts of FM locations (as determined via CT), measured from an absolute reference point and in terms of inter-FM distance. The majority (76%) of points exhibited zero shift, with others shifting by one pixel size consistent with setup error as confirmed with a control sample. Comparison of dose profiles and 2D isodose distributions from the three arrangements illustrated complex spatial redistribution of dose in all three dimensions occurring as a result of the change in shape of the target between irradiations, even for a relatively simple deformation. Discrepancies of up to 30% of the maximum dose were evident from dose difference maps for three orthogonal planes taken through the isocenter of a stereotactic field. Conclusions: This paper describes the first use of a tissue-equivalent, 3D dose-integrating deformable phantom that yields integrated or redistributed dosimetric information. The proposed methodology readily yields three-dimensional (3D) dosimetric data from radiation delivery to the DEFGEL phantom in deformed and undeformed states. The impacts of deformation on dose distributions were readily seen in the isodose contours and line profiles from the three arrangements. It is demonstrated that the system is potentially capable of reproducibly emulating the physical deformation of an organ, and therefore can be used to evaluate absorbed doses to deformable targets and organs at risk in three dimensions and to validate deformation algorithms applied to dose distributions.

61 RADIATION PROTECTION AND DOSIMETRY↗

Technical Note: Modeling a complex micro-multileaf collimator using the standard BEAMnrc distribution

Purpose: The component modules in the standard BEAMnrc distribution may appear to be insufficient to model micro-multileaf collimators that have trifaceted leaf ends and complex leaf profiles. This note indicates, however, that accurate Monte Carlo simulations of radiotherapy beams defined by a complex collimation device can be completed using BEAMnrc's standard VARMLC component module. Methods: That this simple collimator model can produce spatially and dosimetrically accurate microcollimated fields is illustrated using comparisons with ion chamber and film measurements of the dose deposited by square and irregular fields incident on planar, homogeneous water phantoms. Results: Monte Carlo dose calculations for on-axis and off-axis fields are shown to produce good agreement with experimental values, even on close examination of the penumbrae. Conclusions: The use of a VARMLC model of the micro-multileaf collimator, along with a commissioned model of the associated linear accelerator, is therefore recommended as an alternative to the development or use of in-house or third-party component modules for simulating stereotactic radiotherapy and radiosurgery treatments. Simulation parameters for the VARMLC model are provided which should allow other researchers to adapt and use this model to study clinical stereotactic radiotherapy treatments.

62 RADIOLOGY AND NUCLEAR MEDICINE↗