Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DOSIMETRY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Chelation Model Validation: Modeling of a Plutonium-238 Inhalation Incident Treated with DTPA at Los Alamos National Laboratory

Accidental inhalation of plutonium at the workplace is a non-negligible risk, even when rigorous safety standards are in place. The intake and retention of plutonium in the human body may be a source of concern. Thus, if there is a suspicion of a significant intake of plutonium, medical countermeasures such as chelation treatment may be administered to the worker. The present work aimed to interpret the bioassay data of a worker involved in an inhalation incident due to a glovebox breach at Los Alamos National Laboratory’s plutonium facility. The worker was treated with intravenous injections of calcium salts of diethylenetriaminepentaacetic acid (DTPA) in an attempt to reduce the amount of plutonium from the body and therefore reduce the internal radiation dose. It is well known in the internal dosimetry field that the administration of chelation treatment poses additional challenges to the dose assessment. Hence, a recently developed chelation model was used for the modeling of the bioassay data. The objectives of this work are to describe the incident, model the chelation-affected and non-affected bioassay data, estimate the plutonium intake, and assess the internal radiation dose.

238Pu↗
DOE OSTI2015

The development of remote wireless radiation dose monitoring system

Internet of things (IoT) technology has recently shown a large flow of IT trends in human life. In particular, our lives are now becoming integrated with a lot of items around the 'smart-phone' with IoT, including Bluetooth, Near Field Communication (NFC), Beacons, WiFi, and Global Positioning System (GPS). Our project focuses on the interconnection of radiation dosimetry and IoT technology. The radiation workers at a nuclear facility should hold personal dosimeters such as a Thermo-Luminescence Dosimeter (TLD), an Optically Stimulated Luminescence Dosimeter (OSL), pocket ionization chamber dosimeters, an Electronic Personal Dosimeter (EPD), or an alarm dosimeter on their body. Some of them have functions that generate audible or visible alarms to radiation workers in a real working area. However, such devices used in radiation fields these days have no functions for communicating with other areas or the responsible personnel in real time. In particular, when conducting a particular task in a high dose area, or a number of repair works within a radiation field, radiation dose monitoring is important for the health of the workers and the work efficiency. Our project aims at the development of a remote wireless radiation dose monitoring system (RWRD) that can be used to monitor the radiation dose in a nuclear facility for radiation workers and a radiation protection program In this project, a radiation dosimeter is the detection device for personal radiation dose, a smart phone is the mobile wireless communication tool, and, Beacon is the wireless starter for the detection, communication, and position of the worker using BLE (Bluetooth Low Energy). In this report, we report the design of the RWRD and a demonstration case in a real radiation field. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Projection imaging of photon beams by the Cerenkov effect

Purpose: A novel technique for beam profiling of megavoltage photon beams was investigated for the first time by capturing images of the induced Cerenkov emission in water, as a potential surrogate for the imparted dose in irradiated media. Methods: A high-sensitivity, intensified CCD camera (ICCD) was configured to acquire 2D projection images of Cerenkov emission from a 4 Multiplication-Sign 4 cm{sup 2} 6 MV linear accelerator (LINAC) x-ray photon beam operating at a dose rate of 400 MU/min incident on a water tank with transparent walls. The ICCD acquisition was gated to the LINAC sync pulse to reduce background light artifacts, and the measurement quality was investigated by evaluating the signal to noise ratio and measurement repeatability as a function of delivered dose. Monte Carlo simulations were used to derive a calibration factor for differences between the optical images and deposited dose arising from the anisotropic angular dependence of Cerenkov emission. Finally, Cerenkov-based beam profiles were compared to a percent depth dose (PDD) and lateral dose profile at a depth of d{sub max} from a reference dose distribution generated from the clinical Varian ECLIPSE treatment planning system (TPS). Results: The signal to noise ratio was found to be 20 at a delivered dose of 66.6 cGy, and proportional to the square root of the delivered dose as expected from Poisson photon counting statistics. A 2.1% mean standard deviation and 5.6% maximum variation in successive measurements were observed, and the Monte Carlo derived calibration factor resulted in Cerenkov emission images which were directly correlated to deposited dose, with some spatial issues. The dose difference between the TPS and PDD predicted by Cerenkov measurements was within 20% in the buildup region with a distance to agreement (DTA) of 1.5-2 mm and {+-}3% at depths beyond d{sub max}. In the lateral profile, the dose difference at the beam penumbra was within {+-}13% with a DTA of 0-2 mm, {+-}5% in the central beam region, and 2%-3% in the beam umbra. Conclusions: The results from this initial study demonstrate the first documented use of Cerenkov emission imaging to profile x-ray photon LINAC beams in water. The proposed modality has several potential advantages over alternative methods, and upon future refinement may prove to be a robust and novel dosimetry method.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

The effect of head size/shape, miscentering, and bowtie filter on peak patient tissue doses from modern brain perfusion 256-slice CT: How can we minimize the risk for deterministic effects?

Purpose: To determine patient-specific absorbed peak doses to skin, eye lens, brain parenchyma, and cranial red bone marrow (RBM) of adult individuals subjected to low-dose brain perfusion CT studies on a 256-slice CT scanner, and investigate the effect of patient head size/shape, head position during the examination and bowtie filter used on peak tissue doses. Methods: The peak doses to eye lens, skin, brain, and RBM were measured in 106 individual-specific adult head phantoms subjected to the standard low-dose brain perfusion CT on a 256-slice CT scanner using a novel Monte Carlo simulation software dedicated for patient CT dosimetry. Peak tissue doses were compared to corresponding thresholds for induction of cataract, erythema, cerebrovascular disease, and depression of hematopoiesis, respectively. The effects of patient head size/shape, head position during acquisition and bowtie filter used on resulting peak patient tissue doses were investigated. The effect of eye-lens position in the scanned head region was also investigated. The effect of miscentering and use of narrow bowtie filter on image quality was assessed. Results: The mean peak doses to eye lens, skin, brain, and RBM were found to be 124, 120, 95, and 163 mGy, respectively. The effect of patient head size and shape on peak tissue doses was found to be minimal since maximum differences were less than 7%. Patient head miscentering and bowtie filter selection were found to have a considerable effect on peak tissue doses. The peak eye-lens dose saving achieved by elevating head by 4 cm with respect to isocenter and using a narrow wedge filter was found to approach 50%. When the eye lies outside of the primarily irradiated head region, the dose to eye lens was found to drop to less than 20% of the corresponding dose measured when the eye lens was located in the middle of the x-ray beam. Positioning head phantom off-isocenter by 4 cm and employing a narrow wedge filter results in a moderate reduction of signal-to-noise ratio mainly to the peripheral region of the phantom. Conclusions: Despite typical peak doses to skin, eye lens, brain, and RBM from the standard low-dose brain perfusion 256-slice CT protocol are well below the corresponding thresholds for the induction of erythema, cataract, cerebrovascular disease, and depression of hematopoiesis, respectively, every effort should be made toward optimization of the procedure and minimization of dose received by these tissues. The current study provides evidence that the use of the narrower bowtie filter available may considerably reduce peak absorbed dose to all above radiosensitive tissues with minimal deterioration in image quality. Considerable reduction in peak eye-lens dose may also be achieved by positioning patient head center a few centimeters above isocenter during the exposure.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Evaluation of clinical IMRT treatment planning using the GATE Monte Carlo simulation platform for absolute and relative dose calculations

Purpose: The objective of this study was to evaluate and validate the use of the Geant4 application for emission tomography (GATE) Monte Carlo simulation platform for clinical intensity modulated radiotherapy (IMRT) dosimetry studies. Methods: The first step consisted of modeling a 6 MV photon beam linear accelerator (LINAC), with its corresponding validation carried out using percent depth dose evaluation, transverse profiles, tissue phantom ratio, and output factor on water phantom. The IMRT evaluation was performed by comparing simulation and measurements in terms of absolute and relative doses using IMRT dedicated quality assurance phantoms considering seven different patient datasets. Results: Concerning the LINAC simulated model validation tissue phantom ratios at 20 and 10 cm in water TPR{sub 10}{sup 20} obtained from GATE and measurements were 0.672 {+-} 0.063 and 0.675, respectively. In terms of percent depth dose and transverse profiles, error ranges were, respectively: 1.472%{+-} 0.285% and 4.827%{+-} 1.323% for field size of 4 Multiplication-Sign 4, 5 Multiplication-Sign 5, 10 Multiplication-Sign 10, 15 Multiplication-Sign 15, 20 Multiplication-Sign 20, 25 Multiplication-Sign 25, 30 Multiplication-Sign 30, and 40 Multiplication-Sign 40 cm{sup 2}. Most errors were observed at the edge of radiation fields because of higher dose gradient in these areas. Output factors showed good agreement between simulation and measurements with a maximum error of 1.22%. Finally, for IMRT simulations considering seven patient datasets, GATE provided good results with a relative error of 0.43%{+-} 0.25% on absolute dose between simulated and measured beams (measurements at the isocenter, volume 0.125 cm{sup 3}). Planar dose comparisons were also performed using gamma-index analysis. For the whole set of beams considered the mean gamma-index value was 0.497 {+-} 0.152 and 90.8%{+-} 3.6% of the evaluated dose points satisfied the 5%/ 4 mm criterion. Conclusions: These results show that GATE allows reliable simulation of complex beams in radiotherapy after an accurate LINAC modeling is validated. A simple cross-calibration procedure proposed in this work allows obtaining absolute dose values even in complex fields.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Analysis of a free-running synchronization artifact correction for MV-imaging with aSi:H flat panels

Purpose: Solid state flat panel electronic portal imaging devices (EPIDs) are widely used for megavolt (MV) photon imaging applications in radiotherapy. In addition to their original purpose in patient position verification, they are convenient to use in quality assurance and dosimetry to verify beam geometry and dose deposition or to perform linear accelerator (linac) calibration procedures. However, native image frames from amorphous silicon (aSi:H) detectors show a range of artifacts which have to be eliminated by proper correction algorithms. When a panel is operated in free-running frame acquisition mode, moving vertical stripes (periodic synchronization artifacts) are a disturbing feature in image frames. Especially for applications in volumetric intensity modulated arc therapy (VMAT) or motion tracking, the synchronization (sync) artifacts are the limiting factor for potential and accuracy since they become even worse at higher frame rates and at lower dose rates, i.e., linac pulse repetition frequencies (PRFs). Methods: The authors introduced a synchronization correction method which is based on a theoretical model describing the interferences of the panel's readout clocking with the linac's dose pulsing. Depending on the applied PRF, a certain number of dose pulses is captured per frame which is readout columnwise, sequentially. The interference of the PRF with the panel readout is responsible for the period and the different gray value levels of the sync stripes, which can be calculated analytically. Sync artifacts can then be eliminated multiplicatively in precorrected frames without additional information about radiation pulse timing. Results: For the analysis, three aSi:H EPIDs of various types were investigated with 6 and 15 MV photon beams at varying PRFs of 25, 50, 100, 200, and 400 pulses per second. Applying the sync correction at panels with gadolinium oxysulfide scintillators improved single frame flood field image quality drastically [improvement of the signal-to-noise ratio (SNR) up to 66.1 dB for 6 MV and 66.0 dB for 15 MV]. Also for the EPID with a caesium iodide scintillator, the noise for the lower PRFs could be reduced (SNR at 6 MV of up to 56.3 dB and at 15 MV up to 46.7 dB). However, the simplistic readout interference model fails at higher PRFs, where image lag and ghosting effects due to trapped charges in the thin film transistor and scintillator postglowing require additional corrections. Conclusions: The presented free-running sync correction method improves SNR of single frames and enables imaging applications, like low-dose rate imaging at increased image frame rates (e.g., to track moving gold fiducials in the lung). Adaptive image guided radiotherapy protocols become even feasible in VMAT plans. Also simultaneous kilovolt and MV imaging applications can benefit from new possibilities of MV scatter removal in x-ray images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

BrachyView: Proof-of-principle of a novel in-body gamma camera for low dose-rate prostate brachytherapy

Purpose: The conformity of the achieved dose distribution to the treatment plan strongly correlates with the accuracy of seed implantation in a prostate brachytherapy treatment procedure. Incorrect seed placement leads to both short and long term complications, including urethral and rectal toxicity. The authors present BrachyView, a novel concept of a fast intraoperative treatment planning system, to provide real-time seed placement information based on in-body gamma camera data. BrachyView combines the high spatial resolution of a pixellated silicon detector (Medipix2) with the volumetric information acquired by a transrectal ultrasound (TRUS). The two systems will be embedded in the same probe so as to provide anatomically correct seed positions for intraoperative planning and postimplant dosimetry. Dosimetric calculations are based on the TG-43 method using the real position of the seeds. The purpose of this paper is to demonstrate the feasibility of BrachyView using the Medipix2 pixel detector and a pinhole collimator to reconstruct the real-time 3D position of low dose-rate brachytherapy seeds in a phantom. Methods: BrachyView incorporates three Medipix2 detectors coupled to a multipinhole collimator. Three-dimensionally triangulated seed positions from multiple planar images are used to determine the seed placement in a PMMA prostate phantom in real time. MATLAB codes were used to test the reconstruction method and to optimize the device geometry. Results: The results presented in this paper show a 3D position reconstruction accuracy of the seed in the range of 0.5-3 mm for a 10-60 mm seed-to-detector distance interval (Z direction), respectively. The BrachyView system also demonstrates a spatial resolution of 0.25 mm in the XY plane for sources at 10 mm distance from Medipix2 detector plane, comparable to the theoretical value calculated for an equivalent gamma camera arrangement. The authors successfully demonstrated the capability of BrachyView for real-time imaging (using a 3 s data acquisition time) of different brachytherapy seed configurations (with an activity of 0.05 U) throughout a 60 Multiplication-Sign 60 Multiplication-Sign 60 mm{sup 3} Perspex prostate phantom. Conclusions: The newly developed miniature gamma camera component of BrachyView, with its high spatial resolution and real time capability, allows accurate 3D localization of seeds in a prostate phantom. Combination of the gamma camera with TRUS in a single probe will complete the BrachyView system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

On the use of a single-fiber multipoint plastic scintillation detector for {sup 192}Ir high-dose-rate brachytherapy

Purpose: The goal of this study was to prove the feasibility of using a single-fiber multipoint plastic scintillation detector (mPSD) as an in vivo verification tool during {sup 192}Ir high-dose-rate brachytherapy treatments.Methods: A three-point detector was built and inserted inside a catheter-positioning template placed in a water phantom. A hyperspectral approach was implemented to discriminate the different optical signals composing the light output at the exit of the single collection optical fiber. The mPSD was tested with different source-to-detector positions, ranging from 1 to 5 cm radially and over 10.5 cm along the longitudinal axis of the detector, and with various integration times. Several strategies for improving the accuracy of the detector were investigated. The device's accuracy in detecting source position was also tested.Results: Good agreement with the expected doses was obtained for all of the scintillating elements, with average relative differences from the expected values of 3.4 {+-} 2.1%, 3.0 {+-} 0.7%, and 4.5 {+-} 1.0% for scintillating elements from the distal to the proximal. A dose threshold of 3 cGy improved the general accuracy of the detector. An integration time of 3 s offered a good trade-off between precision and temporal resolution. Finally, the mPSD measured the radioactive source positioning uncertainty to be no more than 0.32 {+-} 0.06 mm. The accuracy and precision of the detector were improved by a dose-weighted function combining the three measurement points and known details about the geometry of the detector construction.Conclusions: The use of a mPSD for high-dose-rate brachytherapy dosimetry is feasible. This detector shows great promise for development of in vivo applications for real-time verification of treatment delivery.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Determination of air-kerma strength for the {sup 192}Ir GammaMedplus iX pulsed-dose-rate brachytherapy source

Purpose: Pulsed-dose-rate (PDR) brachytherapy was originally proposed to combine the therapeutic advantages of high-dose-rate (HDR) and low-dose-rate brachytherapy. Though uncommon in the United States, several facilities employ pulsed-dose-rate brachytherapy in Europe and Canada. Currently, there is no air-kerma strength standard for PDR brachytherapy {sup 192}Ir sources traceable to the National Institute of Standards and Technology. Discrepancies in clinical measurements of the air-kerma strength of the PDR brachytherapy sources using HDR source-calibrated well chambers warrant further investigation.Methods: In this research, the air-kerma strength for an {sup 192}Ir PDR brachytherapy source was compared with the University of Wisconsin Accredited Dosimetry Calibration Laboratory transfer standard well chambers, the seven-distance technique [B. E. Rasmussen et al., 'The air-kerma strength standard for 192Ir HDR sources,' Med. Phys. 38, 6721-6729 (2011)], and the manufacturer's stated value. Radiochromic film and Monte Carlo techniques were also employed for comparison to the results of the measurements.Results: While the measurements using the seven-distance technique were within + 0.44% from the manufacturer's determination, there was a + 3.10% difference between the transfer standard well chamber measurements and the manufacturer's stated value. Results showed that the PDR brachytherapy source has geometric and thus radiological qualities that exhibit behaviors similar to a point source model in contrast to a conventional line source model.Conclusions: The resulting effect of the pointlike characteristics of the PDR brachytherapy source likely account for the differences observed between well chamber and in-air measurements.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Comparison of Gafchromic EBT2 and EBT3 for patient-specific quality assurance: Cranial stereotactic radiosurgery using volumetric modulated arc therapy with multiple noncoplanar arcs

Purpose: Patient-specific quality assurance in volumetric modulated arc therapy (VMAT) brain stereotactic radiosurgery raises specific issues on dosimetric procedures, mainly represented by the small radiation fields associated with the lack of lateral electronic equilibrium, the need of small detectors and the high dose delivered (up to 30 Gy). Gafchromic{sup TM} EBT2 and EBT3 films may be considered the dosimeter of choice, and the authors here provide some additional data about uniformity correction for this new generation of radiochromic films.Methods: A new analysis method using blue channel for marker dye correction was proposed for uniformity correction both for EBT2 and EBT3 films. Symmetry, flatness, and field-width of a reference field were analyzed to provide an evaluation in a high-spatial resolution of the film uniformity for EBT3. Absolute doses were compared with thermoluminescent dosimeters (TLD) as baseline. VMAT plans with multiple noncoplanar arcs were generated with a treatment planning system on a selected pool of eleven patients with cranial lesions and then recalculated on a water-equivalent plastic phantom by Monte Carlo algorithm for patient-specific QA. 2D quantitative dose comparison parameters were calculated, for the computed and measured dose distributions, and tested for statistically significant differences.Results: Sensitometric curves showed a different behavior above dose of 5 Gy for EBT2 and EBT3 films; with the use of inhouse marker-dye correction method, the authors obtained values of 2.5% for flatness, 1.5% of symmetry, and a field width of 4.8 cm for a 5 × 5 cm{sup 2} reference field. Compared with TLD and selecting a 5% dose tolerance, the percentage of points with ICRU index below 1 was 100% for EBT2 and 83% for EBT3. Patients analysis revealed statistically significant differences (p < 0.05) between EBT2 and EBT3 in the percentage of points with gamma values <1 (p= 0.009 and p= 0.016); the percent difference as well as the mean difference between calculated and measured isodoses (20% and 80%) were found not to be significant (p= 0.074, p= 0.185, and p= 0.57).Conclusions: Excellent performances in terms of dose homogeneity were obtained using a new blue channel method for marker-dye correction on both EBT2 and EBT3 Gafchromic{sup TM} films. In comparison with TLD, the passing rates for the EBT2 film were higher than for EBT3; a good agreement with estimated data by Monte Carlo algorithm was found for both films, with some statistically significant differences again in favor of EBT2. These results suggest that the use of Gafchromic{sup TM} EBT2 and EBT3 films is appropriate for dose verification measurements in VMAT stereotactic radiosurgery; taking into account the uncertainty associated with Gafchromic film dosimetry, the use of adequate action levels is strongly advised, in particular, for EBT3.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Quality assurance of proton beams using a multilayer ionization chamber system

Purpose: The measurement of percentage depth-dose (PDD) distributions for the quality assurance of clinical proton beams is most commonly performed with a computerized water tank dosimetry system with ionization chamber, commonly referred to as water tank. Although the accuracy and reproducibility of this method is well established, it can be time-consuming if a large number of measurements are required. In this work the authors evaluate the linearity, reproducibility, sensitivity to field size, accuracy, and time-savings of another system: the Zebra, a multilayer ionization chamber system.Methods: The Zebra, consisting of 180 parallel-plate ionization chambers with 2 mm resolution, was used to measure depth-dose distributions. The measurements were performed for scattered and scanned proton pencil beams of multiple energies delivered by the Hitachi PROBEAT synchrotron-based delivery system. For scattered beams, the Zebra-measured depth-dose distributions were compared with those measured with the water tank. The principal descriptors extracted for comparisons were: range, the depth of the distal 90% dose; spread-out Bragg peak (SOBP) length, the region between the proximal 95% and distal 90% dose; and distal-dose fall off (DDF), the region between the distal 80% and 20% dose. For scanned beams, the Zebra-measured ranges were compared with those acquired using a Bragg peak chamber during commissioning.Results: The Zebra demonstrated better than 1% reproducibility and monitor unit linearity. The response of the Zebra was found to be sensitive to radiation field sizes greater than 12.5 × 12.5 cm; hence, the measurements used to determine accuracy were performed using a field size of 10 × 10 cm. For the scattered proton beams, PDD distributions showed 1.5% agreement within the SOBP, and 3.8% outside. Range values agreed within −0.1 ± 0.4 mm, with a maximum deviation of 1.2 mm. SOBP length values agreed within 0 ± 2 mm, with a maximum deviation of 6 mm. DDF values agreed within 0.3 ± 0.1 mm, with a maximum deviation of 0.6 mm. For the scanned proton pencil beams, Zebra and Bragg peak chamber range values demonstrated agreement of 0.0 ± 0.3 mm with a maximum deviation of 1.3 mm. The setup and measurement time for all Zebra measurements was 3 and 20 times less, respectively, compared to the water tank measurements.Conclusions: Our investigation shows that the Zebra can be useful not only for fast but also for accurate measurements of the depth-dose distributions of both scattered and scanned proton beams. The analysis of a large set of measurements shows that the commonly assessed beam quality parameters obtained with the Zebra are within the acceptable variations specified by the manufacturer for our delivery system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

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↗
DOE OSTI2013

Molecular radiotherapy: The NUKFIT software for calculating the time-integrated activity coefficient

Purpose: Calculation of the time-integrated activity coefficient (residence time) is a crucial step in dosimetry for molecular radiotherapy. However, available software is deficient in that it is either not tailored for the use in molecular radiotherapy and/or does not include all required estimation methods. The aim of this work was therefore the development and programming of an algorithm which allows for an objective and reproducible determination of the time-integrated activity coefficient and its standard error.Methods: The algorithm includes the selection of a set of fitting functions from predefined sums of exponentials and the choice of an error model for the used data. To estimate the values of the adjustable parameters an objective function, depending on the data, the parameters of the error model, the fitting function and (if required and available) Bayesian information, is minimized. To increase reproducibility and user-friendliness the starting values are automatically determined using a combination of curve stripping and random search. Visual inspection, the coefficient of determination, the standard error of the fitted parameters, and the correlation matrix are provided to evaluate the quality of the fit. The functions which are most supported by the data are determined using the corrected Akaike information criterion. The time-integrated activity coefficient is estimated by analytically integrating the fitted functions. Its standard error is determined assuming Gaussian error propagation. The software was implemented using MATLAB.Results: To validate the proper implementation of the objective function and the fit functions, the results of NUKFIT and SAAM numerical, a commercially available software tool, were compared. The automatic search for starting values was successfully tested for reproducibility. The quality criteria applied in conjunction with the Akaike information criterion allowed the selection of suitable functions. Function fit parameters and their standard error estimated by using SAAM numerical and NUKFIT showed differences of <1%. The differences for the time-integrated activity coefficients were also <1% (standard error between 0.4% and 3%). In general, the application of the software is user-friendly and the results are mathematically correct and reproducible. An application of NUKFIT is presented for three different clinical examples.Conclusions: The software tool with its underlying methodology can be employed to objectively and reproducibly estimate the time integrated activity coefficient and its standard error for most time activity data in molecular radiotherapy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Strategies for automatic online treatment plan reoptimization using clinical treatment planning system: A planning parameters study

Purpose: Adaptive radiation therapy for prostate cancer using online reoptimization provides an improved control of interfractional anatomy variations. However, the clinical implementation of online reoptimization is currently limited by the low efficiency of current strategies and the difficulties associated with integration into the current treatment planning system. This study investigates the strategies for performing fast (∼2 min) automatic online reoptimization with a clinical fluence-map-based treatment planning system; and explores the performance with different input parameters settings: dose-volume histogram (DVH) objective settings, starting stage, and iteration number (in the context of real time planning).Methods: Simulated treatments of 10 patients were reoptimized daily for the first week of treatment (5 fractions) using 12 different combinations of optimization strategies. Options for objective settings included guideline-based RTOG objectives, patient-specific objectives based on anatomy on the planning CT, and daily-CBCT anatomy-based objectives adapted from planning CT objectives. Options for starting stages involved starting reoptimization with and without the original plan's fluence map. Options for iteration numbers were 50 and 100. The adapted plans were then analyzed by statistical modeling, and compared both in terms of dosimetry and delivery efficiency.Results: All online reoptimized plans were finished within ∼2 min with excellent coverage and conformity to the daily target. The three input parameters, i.e., DVH objectives, starting stage, and iteration number, contributed to the outcome of optimization nearly independently. Patient-specific objectives generally provided better OAR sparing compared to guideline-based objectives. The benefit in high-dose sparing from incorporating daily anatomy into objective settings was positively correlated with the relative change in OAR volumes from planning CT to daily CBCT. The use of the original plan fluence map as the starting stage reduced OAR dose at the mid-dose region, but increased the monitor units by 17%. Differences of only 2cc or less in OAR V50%/V70Gy/V76Gy were observed between 100 and 50 iterations.Conclusions: It is feasible to perform automatic online reoptimization in ∼2 min using a clinical treatment planning system. Selecting optimal sets of input parameters is the key to achieving high quality reoptimized plans, and should be based on the individual patient's daily anatomy, delivery efficiency, and time allowed for plan adaptation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Adaptation of the CVT algorithm for catheter optimization in high dose rate brachytherapy

Purpose: An innovative, simple, and fast method to optimize the number and position of catheters is presented for prostate and breast high dose rate (HDR) brachytherapy, both for arbitrary templates or template-free implants (such as robotic templates).Methods: Eight clinical cases were chosen randomly from a bank of patients, previously treated in our clinic to test our method. The 2D Centroidal Voronoi Tessellations (CVT) algorithm was adapted to distribute catheters uniformly in space, within the maximum external contour of the planning target volume. The catheters optimization procedure includes the inverse planning simulated annealing algorithm (IPSA). Complete treatment plans can then be generated from the algorithm for different number of catheters. The best plan is chosen from different dosimetry criteria and will automatically provide the number of catheters and their positions. After the CVT algorithm parameters were optimized for speed and dosimetric results, it was validated against prostate clinical cases, using clinically relevant dose parameters. The robustness to implantation error was also evaluated. Finally, the efficiency of the method was tested in breast interstitial HDR brachytherapy cases.Results: The effect of the number and locations of the catheters on prostate cancer patients was studied. Treatment plans with a better or equivalent dose distributions could be obtained with fewer catheters. A better or equal prostate V100 was obtained down to 12 catheters. Plans with nine or less catheters would not be clinically acceptable in terms of prostate V100 and D90. Implantation errors up to 3 mm were acceptable since no statistical difference was found when compared to 0 mm error (p > 0.05). No significant difference in dosimetric indices was observed for the different combination of parameters within the CVT algorithm. A linear relation was found between the number of random points and the optimization time of the CVT algorithm. Because the computation time decrease with the number of points and that no effects were observed on the dosimetric indices when varying the number of sampling points and the number of iterations, they were respectively fixed to 2500 and to 100. The computation time to obtain ten complete treatments plans ranging from 9 to 18 catheters, with the corresponding dosimetric indices, was 90 s. However, 93% of the computation time is used by a research version of IPSA. For the breast, on average, the Radiation Therapy Oncology Group recommendations would be satisfied down to 12 catheters. Plans with nine or less catheters would not be clinically acceptable in terms of V100, dose homogeneity index, and D90.Conclusions: The authors have devised a simple, fast and efficient method to optimize the number and position of catheters in interstitial HDR brachytherapy. The method was shown to be robust for both prostate and breast HDR brachytherapy. More importantly, the computation time of the algorithm is acceptable for clinical use. Ultimately, this catheter optimization algorithm could be coupled with a 3D ultrasound system to allow real-time guidance and planning in HDR brachytherapy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Remote auditing of radiotherapy facilities using optically stimulated luminescence dosimeters

Purpose: On 1 July 2012, the Australian Clinical Dosimetry Service (ACDS) released its Optically Stimulated Luminescent Dosimeter (OSLD) Level I audit, replacing the previous TLD based audit. The aim of this work is to present the results from this new service and the complete uncertainty analysis on which the audit tolerances are based. Methods: The audit release was preceded by a rigorous evaluation of the InLight® nanoDot OSLD system from Landauer (Landauer, Inc., Glenwood, IL). Energy dependence, signal fading from multiple irradiations, batch variation, reader variation, and dose response factors were identified and quantified for each individual OSLD. The detectors are mailed to the facility in small PMMA blocks, based on the design of the existing Radiological Physics Centre audit. Modeling and measurement were used to determine a factor that could convert the dose measured in the PMMA block, to dose in water for the facility's reference conditions. This factor is dependent on the beam spectrum. The TPR{sub 20,10} was used as the beam quality index to determine the specific block factor for a beam being audited. The audit tolerance was defined using a rigorous uncertainty calculation. The audit outcome is then determined using a scientifically based two tiered action level approach. Audit outcomes within two standard deviations were defined as Pass (Optimal Level), within three standard deviations as Pass (Action Level), and outside of three standard deviations the outcome is Fail (Out of Tolerance). Results: To-date the ACDS has audited 108 photon beams with TLD and 162 photon beams with OSLD. The TLD audit results had an average deviation from ACDS of 0.0% and a standard deviation of 1.8%. The OSLD audit results had an average deviation of −0.2% and a standard deviation of 1.4%. The relative combined standard uncertainty was calculated to be 1.3% (1σ). Pass (Optimal Level) was reduced to ≤2.6% (2σ), and Fail (Out of Tolerance) was reduced to >3.9% (3σ) for the new OSLD audit. Previously with the TLD audit the Pass (Optimal Level) and Fail (Out of Tolerance) were set at ≤4.0% (2σ) and >6.0% (3σ). Conclusions: The calculated standard uncertainty of 1.3% at one standard deviation is consistent with the measured standard deviation of 1.4% from the audits and confirming the suitability of the uncertainty budget derived audit tolerances. The OSLD audit shows greater accuracy than the previous TLD audit, justifying the reduction in audit tolerances. In the TLD audit, all outcomes were Pass (Optimal Level) suggesting that the tolerances were too conservative. In the OSLD audit 94% of the audits have resulted in Pass (Optimal level) and 6% of the audits have resulted in Pass (Action Level). All Pass (Action level) results have been resolved with a repeat OSLD audit, or an on-site ion chamber measurement.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

A practical and theoretical definition of very small field size for radiotherapy output factor measurements

Purpose: This work introduces the concept of very small field size. Output factor (OPF) measurements at these field sizes require extremely careful experimental methodology including the measurement of dosimetric field size at the same time as each OPF measurement. Two quantifiable scientific definitions of the threshold of very small field size are presented. Methods: A practical definition was established by quantifying the effect that a 1 mm error in field size or detector position had on OPFs and setting acceptable uncertainties on OPF at 1%. Alternatively, for a theoretical definition of very small field size, the OPFs were separated into additional factors to investigate the specific effects of lateral electronic disequilibrium, photon scatter in the phantom, and source occlusion. The dominant effect was established and formed the basis of a theoretical definition of very small fields. Each factor was obtained using Monte Carlo simulations of a Varian iX linear accelerator for various square field sizes of side length from 4 to 100 mm, using a nominal photon energy of 6 MV. Results: According to the practical definition established in this project, field sizes ≤15 mm were considered to be very small for 6 MV beams for maximal field size uncertainties of 1 mm. If the acceptable uncertainty in the OPF was increased from 1.0% to 2.0%, or field size uncertainties are 0.5 mm, field sizes ≤12 mm were considered to be very small. Lateral electronic disequilibrium in the phantom was the dominant cause of change in OPF at very small field sizes. Thus the theoretical definition of very small field size coincided to the field size at which lateral electronic disequilibrium clearly caused a greater change in OPF than any other effects. This was found to occur at field sizes ≤12 mm. Source occlusion also caused a large change in OPF for field sizes ≤8 mm. Based on the results of this study, field sizes ≤12 mm were considered to be theoretically very small for 6 MV beams. Conclusions: Extremely careful experimental methodology including the measurement of dosimetric field size at the same time as output factor measurement for each field size setting and also very precise detector alignment is required at field sizes at least ≤12 mm and more conservatively≤15 mm for 6 MV beams. These recommendations should be applied in addition to all the usual considerations for small field dosimetry, including careful detector selection.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Angular dependence of the response of the nanoDot OSLD system for measurements at depth in clinical megavoltage beams

Purpose: The purpose of this investigation was to assess the angular dependence of a commercial optically stimulated luminescence dosimeter (OSLD) dosimetry system in MV x-ray beams at depths beyondd{sub max} and to find ways to mitigate this dependence for measurements in phantoms. Methods: Two special holders were designed which allow a dosimeter to be rotated around the center of its sensitive volume. The dosimeter's sensitive volume is a disk, 5 mm in diameter and 0.2 mm thick. The first holder rotates the disk in the traditional way. It positions the disk perpendicular to the beam (gantry pointing to the floor) in the initial position (0°). When the holder is rotated the angle of the disk towards the beam increases until the disk is parallel with the beam (“edge on,” 90°). This is referred to as Setup 1. The second holder offers a new, alternative measurement position. It positions the disk parallel to the beam for all angles while rotating around its center (Setup 2). Measurements with five to ten dosimeters per point were carried out for 6 MV at 3 and 10 cm depth. Monte Carlo simulations using GEANT4 were performed to simulate the response of the active detector material for several angles. Detector and housing were simulated in detail based on microCT data and communications with the manufacturer. Various material compositions and an all-water geometry were considered. Results: For the traditional Setup 1 the response of the OSLD dropped on average by 1.4% ± 0.7% (measurement) and 2.1% ± 0.3% (Monte Carlo simulation) for the 90° orientation compared to 0°. Monte Carlo simulations also showed a strong dependence of the effect on the composition of the sensitive layer. Assuming the layer to completely consist of the active material (Al{sub 2}O{sub 3}) results in a 7% drop in response for 90° compared to 0°. Assuming the layer to be completely water, results in a flat response within the simulation uncertainty of about 1%. For the new Setup 2, measurements and Monte Carlo simulations found the angular dependence of the dosimeter to be below 1% and within the measurement uncertainty. Conclusions: The dosimeter system exhibits a small angular dependence of approximately 2% which needs to be considered for measurements involving other than normal incident beams angles. This applies in particular to clinicalin vivo measurements where the orientation of the dosimeter is dictated by clinical circumstances and cannot be optimized as otherwise suggested here. When measuring in a phantom, the proposed new setup should be considered. It changes the orientation of the dosimeter so that a coplanar beam arrangement always hits the disk shaped detector material from the thin side and thereby reduces the angular dependence of the response to within the measurement uncertainty of about 1%. This improvement makes the dosimeter more attractive for clinical measurements with multiple coplanar beams in phantoms, as the overall measurement uncertainty is reduced. Similarly, phantom based postal audits can transition from the traditional TLD to the more accurate and convenient OSLD.

61 RADIATION PROTECTION AND DOSIMETRY↗