Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Radiological Data”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A practical method of modeling a treatment couch using cone-beam computed tomography for intensity-modulated radiation therapy and RapidArc treatment delivery

The effect of a treatment couch on dose perturbation is not always fully considered in intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT). In the course of inverse planning radiotherapy techniques, beam parameter optimization may change in the absence of the couch, causing errors in the calculated dose distributions. Although modern treatment planning systems (TPS) include data for the treatment couch components, they are not manufactured identically. Thus, variations in their Hounsfield unit (HU) values may exist. Moreover, a radiotherapy facility may wish to have a third-party custom tabletop installed that is not included by the TPS vendor. This study demonstrates a practical and simple method of acquiring reliable computed tomography (CT) data for the treatment couch and shows how the absorbed dose calculated with the modeled treatment couch can differ from that with the default treatment couch found in the TPS. We also experimentally verified that neglecting to incorporate the treatment couch completely in the treatment planning process might result in dose differences of up to 9.5% and 7.3% for 4-MV and 10-MV photon beams, respectively. Furthermore, 20 RapidArc and IMRT cases were used to quantify the change in calculated dose distributions caused by using either the default or modeled couch. From 2-dimensional (2D) ionization chamber array measurements, we observed large dose distribution differences between the measurements and calculations when the couch was omitted that varied according to the planning technique and anatomic site. Thus, incorporating the treatment couch in the dose calculation phase of treatment planning significantly decreases dose calculation errors.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

HDR brachytherapy of rectal cancer using a novel grooved-shielding applicator design

Purpose: The aim of this work was to design a novel high-dose rate (HDR) ({sup 192}Ir) brachytherapy applicator for treatment of rectal carcinomas that uses tungsten shielding for possibly improved dosimetric results over commercial brachytherapy applicator(s).Methods: A set of 15 single-depth applicators and one dual-depth applicator were designed and simulated using Monte Carlo (MCNPX). All applicators simulated were high-density tungsten alloy cylinders, 16-mm in diameter, and 60-mm long, with longitudinal grooves within which an {sup 192}Ir source can be placed. The single-depth designs varied regarding the number and depth of these grooves, ranging from 8 to 16 and 1-mm to 3-mm, respectively. The dual-depth design had ten channels, each of which had two depths at which the source could be placed. Optimized treatment plans were generated for each design on data from 13 treated patients (36 fractions) with asymmetrical clinical target volumes (CTVs). All results were compared against the clinically treated plans which used intracavitary mold applicator (ICMA), as well as a recently designed, highly automated, and collimated intensity modulation device named dynamic modulated brachytherapy (DMBT) device.Results: All applicator designs outperformed the ICMA in every calculated dosimetric criteria, except the total dwell times (∼30% increase). There were clear, but relative, tradeoffs regarding both the number of channels and the depth of each channel. Overall, the 12-channel, 1-mm depth, and 14-channel 2-mm depth designs had the best results of the simpler designs, sparing the healthy rectal tissues the most while achieving comparable CTV coverage with the dose heterogeneity index and lateral spill doses improving by over 10% and the contralateral healthy rectum dose dropping over 30% compared to ICMA. The ten-channel dual-depth design outperformed each single-depth design, yielding the best coverage and sparing.Conclusions: New grooved tungsten HDR-brachytherapy devices have been designed and simulated. The results of this work attest to the capability of these new, highly anisotropic, intelligently shielded applicators to limit dose to healthy tissues while maintaining a conformal prescription dose to the CTV.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-T-644: QuAArC: A 3D VMAT QA System Based On Radiochromic Film and Monte Carlo Simulation of Log Files

Purpose: VMAT involves two main sources of uncertainty: one related to the dose calculation accuracy, and the other linked to the continuous delivery of a discrete calculation. The purpose of this work is to present QuAArC, an alternative VMAT QA system to control and potentially reduce these uncertainties. Methods: An automated MC simulation of log files, recorded during VMAT treatment plans delivery, was implemented in order to simulate the actual treatment parameters. The linac head models and the phase-space data of each Control Point (CP) were simulated using the EGSnrc/BEAMnrc MC code, and the corresponding dose calculation was carried out by means of BEAMDOSE, a DOSXYZnrc code modification. A cylindrical phantom was specifically designed to host films rolled up at different radial distances from the isocenter, for a 3D and continuous dosimetric verification. It also allows axial and/or coronal films and point measurements with several types of ion chambers at different locations. Specific software was developed in MATLAB in order to process and evaluate the dosimetric measurements, which incorporates the analysis of dose distributions, profiles, dose difference maps, and 2D/3D gamma index. It is also possible to obtain the experimental DVH reconstructed on the patient CT, by an optimization method to find the individual contribution corresponding to each CP on the film, taking into account the total measured dose, and the corresponding CP dose calculated by MC. Results: The QuAArC system showed high reproducibility of measurements, and consistency with the results obtained with the commercial system implemented in the verification of the evaluated treatment plans. Conclusion: A VMAT QA system based on MC simulation and high resolution dosimetry with film has been developed for treatment verification. It shows to be useful for the study of the real VMAT capabilities, and also for linac commissioning and evaluation of other verification devices.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Improved proton computed tomography by dual modality image reconstruction

Purpose: Proton computed tomography (CT) is a promising image modality for improving the stopping power estimates and dose calculations for particle therapy. However, the finite range of about 33 cm of water of most commercial proton therapy systems limits the sites that can be scanned from a full 360° rotation. In this paper the authors propose a method to overcome the problem using a dual modality reconstruction (DMR) combining the proton data with a cone-beam x-ray prior. Methods: A Catphan 600 phantom was scanned using a cone beam x-ray CT scanner. A digital replica of the phantom was created in the Monte Carlo code Geant4 and a 360° proton CT scan was simulated, storing the entrance and exit position and momentum vector of every proton. Proton CT images were reconstructed using a varying number of angles from the scan. The proton CT images were reconstructed using a constrained nonlinear conjugate gradient algorithm, minimizing total variation and the x-ray CT prior while remaining consistent with the proton projection data. The proton histories were reconstructed along curved cubic-spline paths. Results: The spatial resolution of the cone beam CT prior was retained for the fully sampled case and the 90° interval case, with the MTF = 0.5 (modulation transfer function) ranging from 5.22 to 5.65 linepairs/cm. In the 45° interval case, the MTF = 0.5 dropped to 3.91 linepairs/cm For the fully sampled DMR, the maximal root mean square (RMS) error was 0.006 in units of relative stopping power. For the limited angle cases the maximal RMS error was 0.18, an almost five-fold improvement over the cone beam CT estimate. Conclusions: Dual modality reconstruction yields the high spatial resolution of cone beam x-ray CT while maintaining the improved stopping power estimation of proton CT. In the case of limited angles, the use of prior image proton CT greatly improves the resolution and stopping power estimate, but does not fully achieve the quality of a 360° proton CT scan.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Geometrical splitting technique to improve the computational efficiency in Monte Carlo calculations for proton therapy

Purpose: To present the implementation and validation of a geometrical based variance reduction technique for the calculation of phase space data for proton therapy dose calculation. Methods: The treatment heads at the Francis H Burr Proton Therapy Center were modeled with a new Monte Carlo tool (TOPAS based on Geant4). For variance reduction purposes, two particle-splitting planes were implemented. First, the particles were split upstream of the second scatterer or at the second ionization chamber. Then, particles reaching another plane immediately upstream of the field specific aperture were split again. In each case, particles were split by a factor of 8. At the second ionization chamber and at the latter plane, the cylindrical symmetry of the proton beam was exploited to position the split particles at randomly spaced locations rotated around the beam axis. Phase space data in IAEA format were recorded at the treatment head exit and the computational efficiency was calculated. Depth-dose curves and beam profiles were analyzed. Dose distributions were compared for a voxelized water phantom for different treatment fields for both the reference and optimized simulations. In addition, dose in two patients was simulated with and without particle splitting to compare the efficiency and accuracy of the technique. Results: A normalized computational efficiency gain of a factor of 10-20.3 was reached for phase space calculations for the different treatment head options simulated. Depth-dose curves and beam profiles were in reasonable agreement with the simulation done without splitting: within 1% for depth-dose with an average difference of (0.2 {+-} 0.4)%, 1 standard deviation, and a 0.3% statistical uncertainty of the simulations in the high dose region; 1.6% for planar fluence with an average difference of (0.4 {+-} 0.5)% and a statistical uncertainty of 0.3% in the high fluence region. The percentage differences between dose distributions in water for simulations done with and without particle splitting were within the accepted clinical tolerance of 2%, with a 0.4% statistical uncertainty. For the two patient geometries considered, head and prostate, the efficiency gain was 20.9 and 14.7, respectively, with the percentages of voxels with gamma indices lower than unity 98.9% and 99.7%, respectively, using 2% and 2 mm criteria. Conclusions: The authors have implemented an efficient variance reduction technique with significant speed improvements for proton Monte Carlo simulations. The method can be transferred to other codes and other treatment heads.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Machine Learning Approaches for Predicting Radiation Therapy Outcomes: A Clinician's Perspective

Radiation oncology has always been deeply rooted in modeling, from the early days of isoeffect curves to the contemporary Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) initiative. In recent years, medical modeling for both prognostic and therapeutic purposes has exploded thanks to increasing availability of electronic data and genomics. One promising direction that medical modeling is moving toward is adopting the same machine learning methods used by companies such as Google and Facebook to combat disease. Broadly defined, machine learning is a branch of computer science that deals with making predictions from complex data through statistical models. These methods serve to uncover patterns in data and are actively used in areas such as speech recognition, handwriting recognition, face recognition, “spam” filtering (junk email), and targeted advertising. Although multiple radiation oncology research groups have shown the value of applied machine learning (ML), clinical adoption has been slow due to the high barrier to understanding these complex models by clinicians. Here, we present a review of the use of ML to predict radiation therapy outcomes from the clinician's point of view with the hope that it lowers the “barrier to entry” for those without formal training in ML. We begin by describing 7 principles that one should consider when evaluating (or creating) an ML model in radiation oncology. We next introduce 3 popular ML methods—logistic regression (LR), support vector machine (SVM), and artificial neural network (ANN)—and critique 3 seminal papers in the context of these principles. Although current studies are in exploratory stages, the overall methodology has progressively matured, and the field is ready for larger-scale further investigation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-I-56: Scan Angle Reduction for a Limited-Angle Intrafraction Verification (LIVE) System

Purpose: To develop a novel adaptive reconstruction strategy to further reduce the scanning angle required by the limited-angle intrafraction verification (LIVE) system for intrafraction verification. Methods: LIVE acquires limited angle MV projections from the exit fluence of the arc treatment beam or during gantry rotation between static beams. Orthogonal limited-angle kV projections are also acquired simultaneously to provide additional information. LIVE considers the on-board 4D-CBCT images as a deformation of the prior 4D-CT images, and solves the deformation field based on deformation models and data fidelity constraint. LIVE reaches a checkpoint after a limited-angle scan, and reconstructs 4D-CBCT for intrafraction verification at the checkpoint. In adaptive reconstruction strategy, a larger scanning angle of 30° is used for the first checkpoint, and smaller scanning angles of 15° are used for subsequent checkpoints. The onboard images reconstructed at the previous adjacent checkpoint are used as the prior images for reconstruction at the current checkpoint. As the algorithm only needs to reconstruct the small deformation occurred between adjacent checkpoints, projections from a smaller scan angle provide enough information for the reconstruction. XCAT was used to simulate tumor motion baseline drift of 2mm along sup-inf direction at every subsequent checkpoint, which are 15° apart. Adaptive reconstruction strategy was used to reconstruct the images at each checkpoint using orthogonal 15° kV and MV projections. Results: Results showed that LIVE reconstructed the tumor volumes accurately using orthogonal 15° kV-MV projections. Volume percentage differences (VPDs) were within 5% and center of mass shifts (COMS) were within 1mm for reconstruction at all checkpoints. Conclusion: It's feasible to use an adaptive reconstruction strategy to further reduce the scan angle needed by LIVE to allow faster and more frequent intrafraction verification to minimize the treatment errors in lung cancer treatments. Grant from Varian Medical System.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Medical Device Recalls in Radiation Oncology: Analysis of US Food and Drug Administration Data, 2002-2015

Purpose: To analyze all recalls involving radiation oncology devices (RODs) from the US Food and Drug Administration (FDA)'s recall database, comparing these with non–radiation oncology device recalls to identify discipline-specific trends that may inform improvements in device safety. Methods and Materials: Recall data on RODs from 2002 to 2015 were sorted into 4 product categories (external beam, brachytherapy, planning systems, and simulation systems). Outcomes included determined cause of recall, recall class (severity), quantity in commerce, time until recall termination (date FDA determines recall is complete), and time since 510(k) approval. Descriptive statistics were performed with linear regression of time-series data. Results for RODs were compared with those for other devices by Pearson χ{sup 2} test for categorical data and 2-sample Kolmogorov-Smirnov test for distributions. Results: There were 502 ROD recalls and 9534 other class II device recalls during 2002 to 2015. Most recalls were for external beam devices (66.7%) and planning systems (22.9%), and recall events peaked in 2011. Radiation oncology devices differed significantly from other devices in all recall outcomes (P≤.04). Recall cause was commonly software related (49% vs 10% for other devices). Recall severity was more often moderate among RODs (97.6% vs 87.2%) instead of severe (0.2% vs 4.4%; P<.001). Time from 510(k) market approval to recall was shorter among RODs (P<.001) and progressively shortened over time. Radiation oncology devices had fewer recalled devices in commerce than other devices (P<.001). Conclusions: Compared with other class II devices, RODs experience recalls sooner after market approval and are trending sooner still. Most of these recalls were moderate in severity, and software issues are prevalent. Comprehensive analysis of recall data can identify areas for device improvement, such as better system design among RODs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Simple Method to Estimate Mean Heart Dose From Hodgkin Lymphoma Radiation Therapy According to Simulation X-Rays

Purpose: To describe a new method to estimate the mean heart dose for Hodgkin lymphoma patients treated several decades ago, using delineation of the heart on radiation therapy simulation X-rays. Mean heart dose is an important predictor for late cardiovascular complications after Hodgkin lymphoma (HL) treatment. For patients treated before the era of computed tomography (CT)-based radiotherapy planning, retrospective estimation of radiation dose to the heart can be labor intensive. Methods and Materials: Patients for whom cardiac radiation doses had previously been estimated by reconstruction of individual treatments on representative CT data sets were selected at random from a case–control study of 5-year Hodgkin lymphoma survivors (n=289). For 42 patients, cardiac contours were outlined on each patient's simulation X-ray by 4 different raters, and the mean heart dose was estimated as the percentage of the cardiac contour within the radiation field multiplied by the prescribed mediastinal dose and divided by a correction factor obtained by comparison with individual CT-based dosimetry. Results: According to the simulation X-ray method, the medians of the mean heart doses obtained from the cardiac contours outlined by the 4 raters were 30 Gy, 30 Gy, 31 Gy, and 31 Gy, respectively, following prescribed mediastinal doses of 25-42 Gy. The absolute-agreement intraclass correlation coefficient was 0.93 (95% confidence interval 0.85-0.97), indicating excellent agreement. Mean heart dose was 30.4 Gy with the simulation X-ray method, versus 30.2 Gy with the representative CT-based dosimetry, and the between-method absolute-agreement intraclass correlation coefficient was 0.87 (95% confidence interval 0.80-0.95), indicating good agreement between the two methods. Conclusion: Estimating mean heart dose from radiation therapy simulation X-rays is reproducible and fast, takes individual anatomy into account, and yields results comparable to the labor-intensive representative CT-based method. This simpler method may produce a meaningful measure of mean heart dose for use in studies of late cardiac complications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Fan-beam intensity modulated proton therapy

Purpose: This paper presents a concept for a proton therapy system capable of delivering intensity modulated proton therapy using a fan beam of protons. This system would allow present and future gantry-based facilities to deliver state-of-the-art proton therapy with the greater normal tissue sparing made possible by intensity modulation techniques.Methods: A method for producing a divergent fan beam of protons using a pair of electromagnetic quadrupoles is described and particle transport through the quadrupole doublet is simulated using a commercially available software package. To manipulate the fan beam of protons, a modulation device is developed. This modulator inserts or retracts acrylic leaves of varying thickness from subsections of the fan beam. Each subsection, or beam channel, creates what effectively becomes a beam spot within the fan area. Each channel is able to provide 0–255 mm of range shift for its associated beam spot, or stop the beam and act as an intensity modulator. Results of particle transport simulations through the quadrupole system are incorporated into the MCNPX Monte Carlo transport code along with a model of the range and intensity modulation device. Several design parameters were investigated and optimized, culminating in the ability to create topotherapy treatment plans using distal-edge tracking on both phantom and patient datasets.Results: Beam transport calculations show that a pair of electromagnetic quadrupoles can be used to create a divergent fan beam of 200 MeV protons over a distance of 2.1 m. The quadrupole lengths were 30 and 48 cm, respectively, with transverse field gradients less than 20 T/m, which is within the range of water-cooled magnets for the quadrupole radii used. MCNPX simulations of topotherapy treatment plans suggest that, when using the distal edge tracking delivery method, many delivery angles are more important than insisting on narrow beam channel widths in order to obtain conformal target coverage. Overall, the sharp distal falloff of a proton depth-dose distribution was found to provide sufficient control over the dose distribution to meet objectives, even with coarse lateral resolution and channel widths as large as 2 cm. Treatment plans on both phantom and patient data show that dose conformity suffers when treatments are delivered from less than approximately ten angles. Treatment time for a sample prostate delivery is estimated to be on the order of 10 min, and neutron production is estimated to be comparable to that found for existing collimated systems.Conclusions: Fan beam proton therapy is a method of delivering intensity modulated proton therapy which may be employed as an alternative to magnetic scanning systems. A fan beam of protons can be created by a set of quadrupole magnets and modified by a dual-purpose range and intensity modulator. This can be used to deliver inversely planned treatments, with spot intensities optimized to meet user defined dose objectives. Additionally, the ability of a fan beam delivery system to effectively treat multiple beam spots simultaneously may provide advantages as compared to spot scanning deliveries.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

SU-E-J-10: A Moving-Blocker-Based Strategy for Simultaneous Megavoltage and Kilovoltage Scatter Correction in Cone-Beam Computed Tomography Image Acquired During Volumetric Modulated Arc Therapy

Purpose: To evaluate a moving-blocker-based approach in estimating and correcting megavoltage (MV) and kilovoltage (kV) scatter contamination in kV cone-beam computed tomography (CBCT) acquired during volumetric modulated arc therapy (VMAT). Methods: XML code was generated to enable concurrent CBCT acquisition and VMAT delivery in Varian TrueBeam developer mode. A physical attenuator (i.e., “blocker”) consisting of equal spaced lead strips (3.2mm strip width and 3.2mm gap in between) was mounted between the x-ray source and patient at a source to blocker distance of 232mm. The blocker was simulated to be moving back and forth along the gantry rotation axis during the CBCT acquisition. Both MV and kV scatter signal were estimated simultaneously from the blocked regions of the imaging panel, and interpolated into the un-blocked regions. Scatter corrected CBCT was then reconstructed from un-blocked projections after scatter subtraction using an iterative image reconstruction algorithm based on constraint optimization. Experimental studies were performed on a Catphan 600 phantom and an anthropomorphic pelvis phantom to demonstrate the feasibility of using moving blocker for MV-kV scatter correction. Results: MV scatter greatly degrades the CBCT image quality by increasing the CT number inaccuracy and decreasing the image contrast, in addition to the shading artifacts caused by kV scatter. The artifacts were substantially reduced in the moving blocker corrected CBCT images in both Catphan and pelvis phantoms. Quantitatively, CT number error in selected regions of interest reduced from 377 in the kV-MV contaminated CBCT image to 38 for the Catphan phantom. Conclusions: The moving-blockerbased strategy can successfully correct MV and kV scatter simultaneously in CBCT projection data acquired with concurrent VMAT delivery. This work was supported in part by a grant from the Cancer Prevention and Research Institute of Texas (RP130109) and a grant from the American Cancer Society (RSG-13-326-01-CCE)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Delivery confirmation of bolus electron conformal therapy combined with intensity modulated x-ray therapy

Purpose: The purpose of this study was to demonstrate that a bolus electron conformal therapy (ECT) dose plan and a mixed beam plan, composed of an intensity modulated x-ray therapy (IMXT) dose plan optimized on top of the bolus ECT plan, can be accurately delivered. Methods: Calculated dose distributions were compared with measured dose distributions for parotid and chest wall (CW) bolus ECT and mixed beam plans, each simulated in a cylindrical polystyrene phantom that allowed film dose measurements. Bolus ECT plans were created for both parotid and CW PTVs (planning target volumes) using 20 and 16 MeV beams, respectively, whose 90% dose surface conformed to the PTV. Mixed beam plans consisted of an IMXT dose plan optimized on top of the bolus ECT dose plan. The bolus ECT, IMXT, and mixed beam dose distributions were measured using radiographic films in five transverse and one sagittal planes for a total of 36 measurement conditions. Corrections for film dose response, effects of edge-on photon irradiation, and effects of irregular phantom optical properties on the Cerenkov component of the film signal resulted in high precision measurements. Data set consistency was verified by agreement of depth dose at the intersections of the sagittal plane with the five measured transverse planes. For these same depth doses, results for the mixed beam plan agreed with the sum of the individual depth doses for the bolus ECT and IMXT plans. The six mean measured planar dose distributions were compared with those calculated by the treatment planning system for all modalities. Dose agreement was assessed using the 4% dose difference and 0.2 cm distance to agreement. Results: For the combined high-dose region and low-dose region, pass rates for the parotid and CW plans were 98.7% and 96.2%, respectively, for the bolus ECT plans and 97.9% and 97.4%, respectively, for the mixed beam plans. For the high-dose gradient region, pass rates for the parotid and CW plans were 93.1% and 94.62%, respectively, for the bolus ECT plans and 89.2% and 95.1%, respectively, for the mixed beam plans. For all regions, pass rates for the parotid and CW plans were 98.8% and 97.3%, respectively, for the bolus ECT plans and 97.5% and 95.9%, respectively, for the mixed beam plans. For the IMXT component of the mixed beam plans, pass rates for the parotid and CW plans were 93.7% and 95.8%. Conclusions: Bolus ECT and mixed beam therapy dose delivery to the phantom were more accurate than IMXT delivery, adding confidence to the use of planning, fabrication, and delivery for bolus ECT tools either alone or as part of mixed beam therapy. The methodology reported in this work could serve as a basis for future standardization of the commissioning of bolus ECT or mixed beam therapy. When applying this technology to patients, it is recommended that an electron dose algorithm more accurate than the pencil beam algorithm, e.g., a Monte Carlo algorithm or analytical transport such as the pencil beam redefinition algorithm, be used for planning to ensure the desired accuracy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Adaptive error detection for HDR/PDR brachytherapy: Guidance for decision making during real-time in vivo point dosimetry

Purpose: This study presents an adaptive error detection algorithm (AEDA) for real-timein vivo point dosimetry during high dose rate (HDR) or pulsed dose rate (PDR) brachytherapy (BT) where the error identification, in contrast to existing approaches, does not depend on an a priori reconstruction of the dosimeter position. Instead, the treatment is judged based on dose rate comparisons between measurements and calculations of the most viable dosimeter position provided by the AEDA in a data driven approach. As a result, the AEDA compensates for false error cases related to systematic effects of the dosimeter position reconstruction. Given its nearly exclusive dependence on stable dosimeter positioning, the AEDA allows for a substantially simplified and time efficient real-time in vivo BT dosimetry implementation. Methods: In the event of a measured potential treatment error, the AEDA proposes the most viable dosimeter position out of alternatives to the original reconstruction by means of a data driven matching procedure between dose rate distributions. If measured dose rates do not differ significantly from the most viable alternative, the initial error indication may be attributed to a mispositioned or misreconstructed dosimeter (false error). However, if the error declaration persists, no viable dosimeter position can be found to explain the error, hence the discrepancy is more likely to originate from a misplaced or misreconstructed source applicator or from erroneously connected source guide tubes (true error). Results: The AEDA applied on twoin vivo dosimetry implementations for pulsed dose rate BT demonstrated that the AEDA correctly described effects responsible for initial error indications. The AEDA was able to correctly identify the major part of all permutations of simulated guide tube swap errors and simulated shifts of individual needles from the original reconstruction. Unidentified errors corresponded to scenarios where the dosimeter position was sufficiently symmetric with respect to error and no-error source position constellations. The AEDA was able to correctly identify all false errors represented by mispositioned dosimeters contrary to an error detection algorithm relying on the original reconstruction. Conclusions: The study demonstrates that the AEDA error identification during HDR/PDR BT relies on a stable dosimeter position rather than on an accurate dosimeter reconstruction, and the AEDA’s capacity to distinguish between true and false error scenarios. The study further shows that the AEDA can offer guidance in decision making in the event of potential errors detected with real-timein vivo point dosimetry.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Adverse Events Involving Radiation Oncology Medical Devices: Comprehensive Analysis of US Food and Drug Administration Data, 1991 to 2015

Purpose: Radiation oncology relies on rapidly evolving technology and highly complex processes. The US Food and Drug Administration collects reports of adverse events related to medical devices. We sought to characterize all events involving radiation oncology devices (RODs) from the US Food and Drug Administration's postmarket surveillance Manufacturer and User Facility Device Experience (MAUDE) database, comparing these with non–radiation oncology devices. Methods and Materials: MAUDE data on RODs from 1991 to 2015 were sorted into 4 product categories (external beam, brachytherapy, planning systems, and simulation systems) and 5 device problem categories (software, mechanical, electrical, user error, and dose delivery impact). Outcomes included whether the device was evaluated by the manufacturer, adverse event type, remedial action, problem code, device age, and time since 510(k) approval. Descriptive statistics were performed with linear regression of time-series data. Results for RODs were compared with those for other devices by the Pearson χ{sup 2} test for categorical data and 2-sample Kolmogorov-Smirnov test for distributions. Results: There were 4234 ROD and 4,985,698 other device adverse event reports. Adverse event reports increased over time, and events involving RODs peaked in 2011. Most ROD reports involved external beam therapy (50.8%), followed by brachytherapy (24.9%) and treatment planning systems (21.6%). The top problem types were software (30.4%), mechanical (20.9%), and user error (20.4%). RODs differed significantly from other devices in each outcome (P<.001). RODs were more likely to be evaluated by the manufacturer after an event (46.9% vs 33.0%) but less likely to be recalled (10.5% vs 37.9%) (P<.001). Device age and time since 510(k) approval were shorter among RODs (P<.001). Conclusions: Compared with other devices, RODs may experience adverse events sooner after manufacture and market approval. Close postmarket surveillance, improved software design, and manufacturer–user training may help mitigate these events.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2013

Transconvolution and the virtual positron emission tomograph-A new method for cross calibration in quantitative PET/CT imaging

Purpose: Positron emission tomography (PET)/computed tomography (CT) measurements on small lesions are impaired by the partial volume effect, which is intrinsically tied to the point spread function of the actual imaging system, including the reconstruction algorithms. The variability resulting from different point spread functions hinders the assessment of quantitative measurements in clinical routine and especially degrades comparability within multicenter trials. To improve quantitative comparability there is a need for methods to match different PET/CT systems through elimination of this systemic variability. Consequently, a new method was developed and tested that transforms the image of an object as produced by one tomograph to another image of the same object as it would have been seen by a different tomograph. The proposed new method, termed Transconvolution, compensates for differing imaging properties of different tomographs and particularly aims at quantitative comparability of PET/CT in the context of multicenter trials. Methods: To solve the problem of image normalization, the theory of Transconvolution was mathematically established together with new methods to handle point spread functions of different PET/CT systems. Knowing the point spread functions of two different imaging systems allows determining a Transconvolution function to convert one image into the other. This function is calculated by convolving one point spread function with the inverse of the other point spread function which, when adhering to certain boundary conditions such as the use of linear acquisition and image reconstruction methods, is a numerically accessible operation. For reliable measurement of such point spread functions characterizing different PET/CT systems, a dedicated solid-state phantom incorporating {sup 68}Ge/{sup 68}Ga filled spheres was developed. To iteratively determine and represent such point spread functions, exponential density functions in combination with a Gaussian distribution were introduced. Furthermore, simulation of a virtual PET system provided a standard imaging system with clearly defined properties to which the real PET systems were to be matched. A Hann window served as the modulation transfer function for the virtual PET. The Hann's apodization properties suppressed high spatial frequencies above a certain critical frequency, thereby fulfilling the above-mentioned boundary conditions. The determined point spread functions were subsequently used by the novel Transconvolution algorithm to match different PET/CT systems onto the virtual PET system. Finally, the theoretically elaborated Transconvolution method was validated transforming phantom images acquired on two different PET systems to nearly identical data sets, as they would be imaged by the virtual PET system. Results: The proposed Transconvolution method matched different PET/CT-systems for an improved and reproducible determination of a normalized activity concentration. The highest difference in measured activity concentration between the two different PET systems of 18.2% was found in spheres of 2 ml volume. Transconvolution reduced this difference down to 1.6%. In addition to reestablishing comparability the new method with its parameterization of point spread functions allowed a full characterization of imaging properties of the examined tomographs. Conclusions: By matching different tomographs to a virtual standardized imaging system, Transconvolution opens a new comprehensive method for cross calibration in quantitative PET imaging. The use of a virtual PET system restores comparability between data sets from different PET systems by exerting a common, reproducible, and defined partial volume effect.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Noise suppressed partial volume correction for cardiac SPECT/CT

Purpose: Partial volume correction (PVC) methods typically improve quantification at the expense of increased image noise and reduced reproducibility. In this study, the authors developed a novel voxel-based PVC method that incorporates anatomical knowledge to improve quantification while suppressing noise for cardiac SPECT/CT imaging. Methods: In the proposed method, the SPECT images were first reconstructed using anatomical-based maximum a posteriori (AMAP) with Bowsher’s prior to penalize noise while preserving boundaries. A sequential voxel-by-voxel PVC approach (Yang’s method) was then applied on the AMAP reconstruction using a template response. This template response was obtained by forward projecting a template derived from a contrast-enhanced CT image, and then reconstructed using AMAP to model the partial volume effects (PVEs) introduced by both the system resolution and the smoothing applied during reconstruction. To evaluate the proposed noise suppressed PVC (NS-PVC), the authors first simulated two types of cardiac SPECT studies: a {sup 99m}Tc-tetrofosmin myocardial perfusion scan and a {sup 99m}Tc-labeled red blood cell (RBC) scan on a dedicated cardiac multiple pinhole SPECT/CT at both high and low count levels. The authors then applied the proposed method on a canine equilibrium blood pool study following injection with {sup 99m}Tc-RBCs at different count levels by rebinning the list-mode data into shorter acquisitions. The proposed method was compared to MLEM reconstruction without PVC, two conventional PVC methods, including Yang’s method and multitarget correction (MTC) applied on the MLEM reconstruction, and AMAP reconstruction without PVC. Results: The results showed that the Yang’s method improved quantification, however, yielded increased noise and reduced reproducibility in the regions with higher activity. MTC corrected for PVE on high count data with amplified noise, although yielded the worst performance among all the methods tested on low-count data. AMAP effectively suppressed noise and reduced the spill-in effect in the low activity regions. However it was unable to reduce the spill-out effect in high activity regions. NS-PVC yielded superior performance in terms of both quantitative assessment and visual image quality while improving reproducibility. Conclusions: The results suggest that NS-PVC may be a promising PVC algorithm for application in low-dose protocols, and in gated and dynamic cardiac studies with low counts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Assessing the Continuum of Event-Based Biosurveillance Through an Operational Lens

This research follows the Updated Guidelines for Evaluating Public Health Surveillance Systems, Recommendations from the Guidelines Working Group, published by the Centers for Disease Control and Prevention nearly a decade ago. Since then, models have been developed and complex systems have evolved with a breadth of disparate data to detect or forecast chemical, biological, and radiological events that have significant impact in the One Health landscape. How the attributes identified in 2001 relate to the new range of event-based biosurveillance (EBB) technologies is unclear. This manuscript frames the continuum of EBB methods, models, and constructs through an operational lens (i.e., aspects and attributes associated with operational considerations in the development, testing, and validation of the EBB methods and models and their use in an operational environment). A 2-day subject matter expert workshop was held to scientifically identify, develop, and vet a set of attributes for the broad range of such operational considerations. Workshop participants identified and described comprehensive attributes for the characterization of EBB. The identified attributes are: (1) event, (2) readiness, (3) operational aspects, (4) geographic coverage, (5) population coverage, (6) input data, (7) output, and (8) cost. Ultimately, the analyses herein discuss the broad scope, complexity, and relevant issues germane to EBB useful in an operational environment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Verification of Reactor Activation Modelling with Physical Characterization Data and the Impact on Long Term Assessments - 20303

Whiteshell Laboratories, located in Manitoba, Canada, provided research facilities for the Canadian nuclear industry since the early 1960's. The facility was centered on Whiteshell Reactor 1 (WR-1), an organically cooled, heavy water moderated nuclear reactor. WR-1 was safely shut down in 1985, defueled and drained, and has been in a safe Storage with Surveillance (Deferred Decommissioning) state until present day. In support of an Environmental Assessment for proposed in-situ disposal of WR-1, detailed characterization information regarding the reactor core components (calandria, calandria tubes, and pressure tubes) was necessary to provide confidence that the long-term impacts of the facility can be assessed adequately. A simplified estimate of the total activation of core components was completed in 1992 using conservative assumptions, expected reactor configuration and computer models built with WIMS-CRNL, ONEDANT and ORIGEN-S. Due to the simplification, and age of the work, there was considerable doubt that the modelled inventory estimate of the core components was sufficient to support long-term performance modeling of the planned in-situ decommissioning. Additional data would be required to provide this confidence. Several options were considered for characterizing the WR-1 core components, but it was decided to collect only a limited number of samples in an effort to validate the 1992 simplified estimate as bounding. The limited number of samples reduced risks to workers, both radiological and non-radiological, and significantly reduced the costs required to confidently define and bound the radiological inventory of activation products in the core. The results of the sampling were compared to the simplified activation estimation. There was good agreement between the modelled and sampled results, with all major contributors to total activity showing similar relative quantities. The samples were generally lower than predicted by the models, with few exceptions. Notable exceptions included Nb-94, which were higher in the sample results than the model. This was determined to be caused by the use of Niobium containing steel alloys in the stainless steel pressure tube that were not accounted for in the models. The comparison of results provided confidence that the original inventory estimates produced through the simplified model are conservative and bounding, and therefore suitable for use in the long-term assessment of in situ disposal of the WR-1 reactor. This paper provides details of the simplified computer model results, their comparison to results of the sampling, and the changes adopted in the assessment of the in-situ decommissioning of WR-1 as a result. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗