Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Dynamic volume vs respiratory correlated 4DCT for motion assessment in radiation therapy simulation

Purpose: Conventional (i.e., respiratory-correlated) 4DCT exploits the repetitive nature of breathing to provide an estimate of motion; however, it has limitations due to binning artifacts and irregular breathing in actual patient breathing patterns. The aim of this work was to evaluate the accuracy and image quality of a dynamic volume, CT approach (4D{sub vol}) using a 320-slice CT scanner to minimize these limitations, wherein entire image volumes are acquired dynamically without couch movement. This will be compared to the conventional respiratory-correlated 4DCT approach (RCCT). Methods: 4D{sub vol} CT was performed and characterized on an in-house, programmable respiratory motion phantom containing multiple geometric and morphological ''tumor'' objects over a range of regular and irregular patient breathing traces obtained from 3D fluoroscopy and compared to RCCT. The accuracy of volumetric capture and breathing displacement were evaluated and compared with the ground truth values and with the results reported using RCCT. A motion model was investigated to validate the number of motion samples needed to obtain accurate motion probability density functions (PDF). The impact of 4D image quality on this accuracy was then investigated. Dose measurements using volumetric and conventional scan techniques were also performed and compared. Results: Both conventional and dynamic volume 4DCT methods were capable of estimating the programmed displacement of sinusoidal motion, but patient breathing is known to not be regular, and obvious differences were seen for realistic, irregular motion. The mean RCCT amplitude error averaged at 4 mm (max. 7.8 mm) whereas the 4D{sub vol} CT error stayed below 0.5 mm. Similarly, the average absolute volume error was lower with 4D{sub vol} CT. Under irregular breathing, the 4D{sub vol} CT method provides a close description of the motion PDF (cross-correlation 0.99) and is able to track each object, whereas the RCCT method results in a significantly different PDF from the ground truth, especially for smaller tumors (cross-correlation ranging between 0.04 and 0.69). For the protocols studied, the dose measurements were higher in the 4D{sub vol} CT method (40%), but it was shown that significant mAs reductions can be achieved by a factor of 4-5 while maintaining image quality and accuracy. Conclusions: 4D{sub vol} CT using a scanner with a large cone-angle is a promising alternative for improving the accuracy with which respiration-induced motion can be characterized, particularly for patients with irregular breathing motion. This approach also generates 4DCT image data with a reduced total scan time compared to a RCCT scan, without the need for image binning or external respiration signals within the 16 cm scan length. Scan dose can be made comparable to RCCT by optimization of the scan parameters. In addition, it provides the possibility of measuring breathing motion for more than one breathing cycle to assess stability and obtain a more accurate motion PDF, which is currently not feasible with the conventional RCCT approach.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Statistical CT noise reduction with multiscale decomposition and penalized weighted least squares in the projection domain

Purposes: The suppression of noise in x-ray computed tomography (CT) imaging is of clinical relevance for diagnostic image quality and the potential for radiation dose saving. Toward this purpose, statistical noise reduction methods in either the image or projection domain have been proposed, which employ a multiscale decomposition to enhance the performance of noise suppression while maintaining image sharpness. Recognizing the advantages of noise suppression in the projection domain, the authors propose a projection domain multiscale penalized weighted least squares (PWLS) method, in which the angular sampling rate is explicitly taken into consideration to account for the possible variation of interview sampling rate in advanced clinical or preclinical applications. Methods: The projection domain multiscale PWLS method is derived by converting an isotropic diffusion partial differential equation in the image domain into the projection domain, wherein a multiscale decomposition is carried out. With adoption of the Markov random field or soft thresholding objective function, the projection domain multiscale PWLS method deals with noise at each scale. To compensate for the degradation in image sharpness caused by the projection domain multiscale PWLS method, an edge enhancement is carried out following the noise reduction. The performance of the proposed method is experimentally evaluated and verified using the projection data simulated by computer and acquired by a CT scanner. Results: The preliminary results show that the proposed projection domain multiscale PWLS method outperforms the projection domain single-scale PWLS method and the image domain multiscale anisotropic diffusion method in noise reduction. In addition, the proposed method can preserve image sharpness very well while the occurrence of 'salt-and-pepper' noise and mosaic artifacts can be avoided. Conclusions: Since the interview sampling rate is taken into account in the projection domain multiscale decomposition, the proposed method is anticipated to be useful in advanced clinical and preclinical applications where the interview sampling rate varies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Investigation of the spatial resolution of an online dose verification device

Purpose: The aim of this work is to characterize a new online dose verification device, COMPASS transmission detector array (IBA Dosimetry, Schwarzenbruck, Germany). The array is composed of 1600 cylindrical ionization chambers of 3.8 mm diameter, separated by 6.5 mm center-to-center spacing, in a 40 x 40 arrangement. Methods: The line spread function (LSF) of a single ion chamber in the detector was measured with a narrow slit collimator for a 6 MV photon beam. The 0.25 x 10 mm{sup 2} slit was formed by two machined lead blocks. The LSF was obtained by laterally translating the detector in 0.25 mm steps underneath the slit over a range of 24 mm and taking a measurement at each step. This measurement was validated with Monte Carlo simulation using BEAMnrc and DOSXYZnrc. The presampling modulation transfer function (MTF), the Fourier transform of the line spread function, was determined and compared to calculated (Monte Carlo and analytical) MTFs. Two head-and-neck intensity modulated radiation therapy (IMRT) fields were measured using the device and were used to validate the LSF measurement. These fields were simulated with the BEAMnrc Monte Carlo model, and the Monte Carlo generated incident fluence was convolved with the 2D detector response function (derived from the measured LSF) to obtain calculated dose. The measured and calculated dose distributions were then quantitatively compared using {chi}-comparison criteria of 3% dose difference and 3 mm distance-to-agreement for in-field points (defined as those above the 10% maximum dose threshold). Results: The full width at half-maximum (FWHM) of the measured detector response for a single chamber is 4.3 mm, which is comparable to the chamber diameter of 3.8 mm. The pre-sampling MTF was calculated, and the resolution of one chamber was estimated as 0.25 lp/mm from the first zero crossing. For both examined IMRT fields, the {chi}-comparison between measured and calculated data show good agreement with 95.1% and 96.3% of in-field points below {chi} of 1.0 for fields 1 and 2, respectively (with an average {chi} of 0.29 for IMRT field 1 and 0.24 for IMRT field 2). Conclusions: The LSF for a new novel online detector has been measured at 6 MV using a narrow slit technique, and this measurement has been validated by Monte Carlo simulation. The detector response function derived from line spread function has been applied to recover measured IMRT fields. The results have shown that the device measures IMRT fields accurately within acceptable tolerance.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

TU-C-17A-03: An Integrated Contour Evaluation Software Tool Using Supervised Pattern Recognition for Radiotherapy

Purpose: Radiotherapy (RT) contours delineated either manually or semiautomatically require verification before clinical usage. Manual evaluation is very time consuming. A new integrated software tool using supervised pattern contour recognition was thus developed to facilitate this process. Methods: The contouring tool was developed using an object-oriented programming language C# and application programming interfaces, e.g. visualization toolkit (VTK). The C# language served as the tool design basis. The Accord.Net scientific computing libraries were utilized for the required statistical data processing and pattern recognition, while the VTK was used to build and render 3-D mesh models from critical RT structures in real-time and 360° visualization. Principal component analysis (PCA) was used for system self-updating geometry variations of normal structures based on physician-approved RT contours as a training dataset. The inhouse design of supervised PCA-based contour recognition method was used for automatically evaluating contour normality/abnormality. The function for reporting the contour evaluation results was implemented by using C# and Windows Form Designer. Results: The software input was RT simulation images and RT structures from commercial clinical treatment planning systems. Several abilities were demonstrated: automatic assessment of RT contours, file loading/saving of various modality medical images and RT contours, and generation/visualization of 3-D images and anatomical models. Moreover, it supported the 360° rendering of the RT structures in a multi-slice view, which allows physicians to visually check and edit abnormally contoured structures. Conclusion: This new software integrates the supervised learning framework with image processing and graphical visualization modules for RT contour verification. This tool has great potential for facilitating treatment planning with the assistance of an automatic contour evaluation module in avoiding unnecessary manual verification for physicians/dosimetrists. In addition, its nature as a compact and stand-alone tool allows for future extensibility to include additional functions for physicians’ clinical needs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Real-Time, Adaptive Radiological Anomaly Detection and Isotope Identification Using Non-Negative Matrix Factorization

Spectroscopic anomaly detection and isotope identification algorithms are integral components in nuclear nonproliferation applications such as search operations. The task is especially challenging in the case of mobile detector systems because the observed gamma-ray background changes more than for a static detector system, and a pretrained background model can easily find itself out of domain. The result is that algorithms may exceed their intended false alarm rate or sacrifice detection sensitivity to maintain the desired false alarm rate. Non-negative matrix factorization (NMF) is a powerful tool for spectral anomaly detection and identification, but, like many similar algorithms that rely on data-driven background models, in its conventional implementation, it is unable to update in real time to account for environmental changes that affect the background spectroscopic signature. Here, we have developed a novel NMF-based algorithm that periodically updates its background model to accommodate changing environmental conditions. The adaptive NMF algorithm involves fewer assumptions about its environment, making it more generalizable than existing NMF-based methods while maintaining or exceeding detection performance on simulated and real-world datasets.

Anomaly detection↗
DOE OSTI2014

An algorithm to estimate the object support in truncated images

Purpose: Truncation artifacts in CT occur if the object to be imaged extends past the scanner field of view (SFOV). These artifacts impede diagnosis and could possibly introduce errors in dose plans for radiation therapy. Several approaches exist for correcting truncation artifacts, but existing correction algorithms do not accurately recover the skin line (or support) of the patient, which is important in some dose planning methods. The purpose of this paper was to develop an iterative algorithm that recovers the support of the object. Methods: The authors assume that the truncated portion of the image is made up of soft tissue of uniform CT number and attempt to find a shape consistent with the measured data. Each known measurement in the sinogram is interpreted as an estimate of missing mass along a line. An initial estimate of the object support is generated by thresholding a reconstruction made using a previous truncation artifact correction algorithm (e.g., water cylinder extrapolation). This object support is iteratively deformed to reduce the inconsistency with the measured data. The missing data are estimated using this object support to complete the dataset. The method was tested on simulated and experimentally truncated CT data. Results: The proposed algorithm produces a better defined skin line than water cylinder extrapolation. On the experimental data, the RMS error of the skin line is reduced by about 60%. For moderately truncated images, some soft tissue contrast is retained near the SFOV. As the extent of truncation increases, the soft tissue contrast outside the SFOV becomes unusable although the skin line remains clearly defined, and in reformatted images it varies smoothly from slice to slice as expected. Conclusions: The support recovery algorithm provides a more accurate estimate of the patient outline than thresholded, basic water cylinder extrapolation, and may be preferred in some radiation therapy applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Correlation of Injury Simulation with Clinical Assessment of Traumatic Brain Injury

This report contains a summary of our efforts to correlate head injury simulations predicting intracranial fluid cavitation with clinical assessments of brain injury from blunt impact to the head. Magnetic resonance imaging (MRI) data, collected on traumatic brain injury (TBI) subjects by researchers at the MIND Institute of New Mexico, was acquired for the current work. Specific blunt impact TBI case histories were selected from the TBI data for further study and possible correlation with simulation. Both group and single-subject case histories were examined. We found one single-subject case that was particularly suited for correlation with simulation. Diffusion tensor image (DTI) analysis of the TBI subject identified white matter regions within the brain displaying reductions in fractional anisotropy (FA), an indicator of local damage to the white matter axonal structures. Analysis of functional magnetic resonance image (fMRI) data collected on this individual identified localized regions of the brain displaying hypoactivity, another indicator of brain injury. We conducted high fidelity simulations of head impact experienced by the TBI subject using the Sandia head-neck-torso model and the shock physics computer code CTH. Intracranial fluid cavitation predictions were compared with maps of DTI fractional anisotropy and fMRI hypoactivity to assess whether a possible correlation exists. The ultimate goal of this work is to assess whether one can correlate simulation predictions of intracranial fluid cavitation with the brain injured sites identified by the fMRI and DTI analyses. The outcome of this effort is described in this report.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

SU-E-T-657: Quantitative Assessment of Plan Robustness for Helical Tomotherapy for Head and Neck Cancer Radiotherapy

Introduction: Geometric uncertainties in daily patient setup can lead to variations in the planned dose, especially when using highly conformal techniques such as helical Tomotherapy. To account for the potential effect of geometric uncertainty, our clinical practice is to expand critical structures by 3mm expansion into planning risk volumes (PRV). The PRV concept assumes the spatial dose cloud is insensitive to patient positioning. However, no tools currently exist to determine if a Tomotherapy plan is robust to the effects of daily setup variation. We objectively quantified the impact of geometric uncertainties on the 3D doses to critical normal tissues during helical Tomotherapy. Methods: Using a Matlab-based program created and validated by Accuray (Madison, WI), the planned Tomotherapy delivery sinogram recalculated dose on shifted CT datasets. Ten head and neck patients were selected for analysis. To simulate setup uncertainty, the patient anatomy was shifted ±3mm in the longitudinal, lateral and vertical axes. For each potential shift, the recalculated doses to various critical normal tissues were compared to the doses delivered to the PRV in the original plan Results: 18 shifted scenarios created from Tomotherapy plans for three patients with head and neck cancers were analyzed. For all simulated setup errors, the maximum doses to the brainstem, spinal cord, parotids and cochlea were no greater than 0.6Gy of the respective original PRV maximum. Despite 3mm setup shifts, the minimum dose delivered to 95% of the CTVs and PTVs were always within 0.4Gy of the original plan. Conclusions: For head and neck sites treated with Tomotherapy, the use of a 3mm PRV expansion provide a reasonable estimate of the dosimetric effects of 3mm setup uncertainties. Similarly, target coverage appears minimally effected by a 3mm setup uncertainty. Data from a larger number of patients will be presented. Future work will include other anatomical sites.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

On the feasibility of optical-CT imaging in media of different refractive index

Purpose: Achieving accurate optical-CT 3D dosimetry without the use of viscous refractive index (RI) matching fluids would greatly increase convenience. Methods: Software has been developed to simulate optical-CT 3D dosimetry for a range of scanning configurations including parallel-beam, point, and converging light sources. For each configuration the efficacy of three refractive media was investigated: air, water, a fluid closely matched to PRESAGE{sup Registered-Sign }, and perfect matching (RI = 1.00, 1.33, 1.49, and 1.501 respectively). Reconstructions were performed using both filtered backprojection (FBP) and algebraic reconstruction technique (ART). The efficacy of the three configurations and the two algorithms was evaluated by calculating the usable radius (i.e., the outermost radius where data were accurate to within 2%), and gamma ({Gamma}) analysis. This definition recognizes that for optical-CT imaging, errors are greatest near the edge of the dosimeter, where refraction can be most pronounced. Simulations were performed on three types of dose distribution: uniform, volumetric modulated arc therapy (VMAT), and brachytherapy (Cs-137). Results: For a uniformly irradiated dosimeter the usable radius achieved with filtered backprojection was 68% for water-matching and 31% for dry-scanning in air. Algebraic reconstruction gave usable radii of 99% for both water and air (dry-scanning), indicating greater recovery of useful data for the uniform distribution. FBP and ART performed equally well for a VMAT dose distribution where less dose is delivered near the edge of the dosimeter. In this case, the usable radius was 86% and 53% for scanning in water and air, respectively. For brachytherapy, the usable radius was 99% and 98% for scanning in water and air, respectively using FBP, and a major decrease was seen with ART. Point source geometry provided 1%-2% larger usable radii than parallel geometry. Converging geometry recovered less usable dosimetry data (up to 10% reduced usable radii) than point and parallel geometries. A further disadvantage of converging geometry was an increased requirement on detector size by up to 18 Degree-Sign . Conclusions: For applications where dose information is not required in the periphery of the dosimeter, some dry and low-viscous matching configurations may be feasible. For all three dose distributions (uniform, VMAT, brachytherapy) the point source geometry produced slightly more favorable results (an extra 1%-2% usable radii) than parallel and converging. When dosimetry is required on the periphery, best results were obtained using close refractive matching and ART. A concern for water or dry-scanning is the increase in required detector size, introducing potential cost penalties for manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2013

2D/3D registration algorithm for lung brachytherapy

Purpose: A 2D/3D registration algorithm is proposed for registering orthogonal x-ray images with a diagnostic CT volume for high dose rate (HDR) lung brachytherapy. Methods: The algorithm utilizes a rigid registration model based on a pixel/voxel intensity matching approach. To achieve accurate registration, a robust similarity measure combining normalized mutual information, image gradient, and intensity difference was developed. The algorithm was validated using a simple body and anthropomorphic phantoms. Transfer catheters were placed inside the phantoms to simulate the unique image features observed during treatment. The algorithm sensitivity to various degrees of initial misregistration and to the presence of foreign objects, such as ECG leads, was evaluated. Results: The mean registration error was 2.2 and 1.9 mm for the simple body and anthropomorphic phantoms, respectively. The error was comparable to the interoperator catheter digitization error of 1.6 mm. Preliminary analysis of data acquired from four patients indicated a mean registration error of 4.2 mm. Conclusions: Results obtained using the proposed algorithm are clinically acceptable especially considering the complications normally encountered when imaging during lung HDR brachytherapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

An adaptive regularization parameter choice strategy for multispectral bioluminescence tomography

Purpose: Bioluminescence tomography (BLT) provides an effective tool for monitoring physiological and pathological activities in vivo. However, the measured data in bioluminescence imaging are corrupted by noise. Therefore, regularization methods are commonly used to find a regularized solution. Nevertheless, for the quality of the reconstructed bioluminescent source obtained by regularization methods, the choice of the regularization parameters is crucial. To date, the selection of regularization parameters remains challenging. With regards to the above problems, the authors proposed a BLT reconstruction algorithm with an adaptive parameter choice rule. Methods: The proposed reconstruction algorithm uses a diffusion equation for modeling the bioluminescent photon transport. The diffusion equation is solved with a finite element method. Computed tomography (CT) images provide anatomical information regarding the geometry of the small animal and its internal organs. To reduce the ill-posedness of BLT, spectral information and the optimal permissible source region are employed. Then, the relationship between the unknown source distribution and multiview and multispectral boundary measurements is established based on the finite element method and the optimal permissible source region. Since the measured data are noisy, the BLT reconstruction is formulated as l{sub 2} data fidelity and a general regularization term. When choosing the regularization parameters for BLT, an efficient model function approach is proposed, which does not require knowledge of the noise level. This approach only requests the computation of the residual and regularized solution norm. With this knowledge, we construct the model function to approximate the objective function, and the regularization parameter is updated iteratively. Results: First, the micro-CT based mouse phantom was used for simulation verification. Simulation experiments were used to illustrate why multispectral data were used rather than monochromatic data. Furthermore, the study conducted using an adaptive regularization parameter demonstrated our ability to accurately localize the bioluminescent source. With the adaptively estimated regularization parameter, the reconstructed center position of the source was (20.37, 31.05, 12.95) mm, and the distance to the real source was 0.63 mm. The results of the dual-source experiments further showed that our algorithm could localize the bioluminescent sources accurately. The authors then presented experimental evidence that the proposed algorithm exhibited its calculated efficiency over the heuristic method. The effectiveness of the new algorithm was also confirmed by comparing it with the L-curve method. Furthermore, various initial speculations regarding the regularization parameter were used to illustrate the convergence of our algorithm. Finally, in vivo mouse experiment further illustrates the effectiveness of the proposed algorithm. Conclusions: Utilizing numerical, physical phantom and in vivo examples, we demonstrated that the bioluminescent sources could be reconstructed accurately with automatic regularization parameters. The proposed algorithm exhibited superior performance than both the heuristic regularization parameter choice method and L-curve method based on the computational speed and localization error.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

Development of CFD models to support LEU Conversion of ORNL s High Flux Isotope Reactor

The US Department of Energy s National Nuclear Security Administration (NNSA) is participating in the Global Threat Reduction Initiative to reduce and protect vulnerable nuclear and radiological materials located at civilian sites worldwide. As an integral part of one of NNSA s subprograms, Reduced Enrichment for Research and Test Reactors, HFIR is being converted from the present HEU core to a low enriched uranium (LEU) core with less than 20% of U-235 by weight. Because of HFIR s importance for condensed matter research in the United States, its conversion to a high-density, U-Mo-based, LEU fuel should not significantly impact its existing performance. Furthermore, cost and availability considerations suggest making only minimal changes to the overall HFIR facility. Therefore, the goal of this conversion program is only to substitute LEU for the fuel type in the existing fuel plate design, retaining the same number of fuel plates, with the same physical dimensions, as in the current HFIR HEU core. Because LEU-specific testing and experiments will be limited, COMSOL Multiphysics was chosen to provide the needed simulation capability to validate against the HEU design data and previous calculations, and predict the performance of the proposed LEU fuel for design and safety analyses. To achieve it, advanced COMSOL-based multiphysics simulations, including computational fluid dynamics (CFD), are being developed to capture the turbulent flows and associated heat transfer in fine detail and to improve predictive accuracy [2].

COMSOL↗
DOE OSTI2013

Quantitative Monte Carlo-based holmium-166 SPECT reconstruction

Purpose: Quantitative imaging of the radionuclide distribution is of increasing interest for microsphere radioembolization (RE) of liver malignancies, to aid treatment planning and dosimetry. For this purpose, holmium-166 ({sup 166}Ho) microspheres have been developed, which can be visualized with a gamma camera. The objective of this work is to develop and evaluate a new reconstruction method for quantitative {sup 166}Ho SPECT, including Monte Carlo-based modeling of photon contributions from the full energy spectrum.Methods: A fast Monte Carlo (MC) simulator was developed for simulation of {sup 166}Ho projection images and incorporated in a statistical reconstruction algorithm (SPECT-fMC). Photon scatter and attenuation for all photons sampled from the full {sup 166}Ho energy spectrum were modeled during reconstruction by Monte Carlo simulations. The energy- and distance-dependent collimator-detector response was modeled using precalculated convolution kernels. Phantom experiments were performed to quantitatively evaluate image contrast, image noise, count errors, and activity recovery coefficients (ARCs) of SPECT-fMC in comparison with those of an energy window-based method for correction of down-scattered high-energy photons (SPECT-DSW) and a previously presented hybrid method that combines MC simulation of photopeak scatter with energy window-based estimation of down-scattered high-energy contributions (SPECT-ppMC+DSW). Additionally, the impact of SPECT-fMC on whole-body recovered activities (A{sup est}) and estimated radiation absorbed doses was evaluated using clinical SPECT data of six {sup 166}Ho RE patients.Results: At the same noise level, SPECT-fMC images showed substantially higher contrast than SPECT-DSW and SPECT-ppMC+DSW in spheres ≥17 mm in diameter. The count error was reduced from 29% (SPECT-DSW) and 25% (SPECT-ppMC+DSW) to 12% (SPECT-fMC). ARCs in five spherical volumes of 1.96–106.21 ml were improved from 32%–63% (SPECT-DSW) and 50%–80% (SPECT-ppMC+DSW) to 76%–103% (SPECT-fMC). Furthermore, SPECT-fMC recovered whole-body activities were most accurate (A{sup est}= 1.06 × A − 5.90 MBq, R{sup 2}= 0.97) and SPECT-fMC tumor absorbed doses were significantly higher than with SPECT-DSW (p = 0.031) and SPECT-ppMC+DSW (p = 0.031).Conclusions: The quantitative accuracy of {sup 166}Ho SPECT is improved by Monte Carlo-based modeling of the image degrading factors. Consequently, the proposed reconstruction method enables accurate estimation of the radiation absorbed dose in clinical practice.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Out-of-Distribution Detection and Radiological Data Monitoring Using Statistical Process Control

Abstract Machine learning (ML) models often fail with data that deviates from their training distribution. This is a significant concern for ML-enabled devices as data drift may lead to unexpected performance. This work introduces a new framework for out of distribution (OOD) detection and data drift monitoring that combines ML and geometric methods with statistical process control (SPC). We investigated different design choices, including methods for extracting feature representations and drift quantification for OOD detection in individual images and as an approach for input data monitoring. We evaluated the framework for both identifying OOD images and demonstrating the ability to detect shifts in data streams over time. We demonstrated a proof-of-concept via the following tasks: 1) differentiating axial vs. non-axial CT images, 2) differentiating CXR vs. other radiographic imaging modalities, and 3) differentiating adult CXR vs. pediatric CXR. For the identification of individual OOD images, our framework achieved high sensitivity in detecting OOD inputs: 0.980 in CT, 0.984 in CXR, and 0.854 in pediatric CXR. Our framework is also adept at monitoring data streams and identifying the time a drift occurred. In our simulations tracking drift over time, it effectively detected a shift from CXR to non-CXR instantly, a transition from axial to non-axial CT within few days, and a drift from adult to pediatric CXRs within a day—all while maintaining a low false positive rate. Through additional experiments, we demonstrate the framework is modality-agnostic and independent from the underlying model structure, making it highly customizable for specific applications and broadly applicable across different imaging modalities and deployed ML models.

Zamzmi, Ghada↗
DOE OSTI2019

EM Program Expectations on Implementation of DOE O 458: Radiation Protection of Public/Environment - 19553

DOE Order 458.1, Radiation Protection of the Public and Environment, was issued in 2011. DOE sites were given 18 months to implement. The revised Order promulgates requirements for property clearance actions and independent verification. DOE Order 458.1 allows for preapproved authorized limits for the clearance of property. However, when property with the potential to contain residual radioactive material does not meet preapproved authorized limits, an authorized limits package must be developed. The Order contains a range of requirements for demonstrating with reasonable assurance that dose constraints are met for any actual or likely future use of the property to be released. EM has developed a draft guidance document that provides the elements of an authorized limits package. These elements are used to review all authorized limits packages that are submitted to EM. Expected elements include property description and process knowledge documentation, current and likely future use, the conceptual site model, dose modeling methodologies, specific authorized limits, dose- or concentration-based constraints, collective dose estimate to the potentially exposed population, stakeholder communications, description of the As Low As Reasonably Achievable (ALARA) analysis, implementation plan, evaluation of difficult to assess surfaces, and independent verification plan. The draft guidance document also provides further explanation of the desired elements of an independent verification (IV) plan. An IV plan should describe the planning phase, implementation phase, assessment phase, and decision phase of an activity. The planning phase includes scope and cost, concurrent with cleanup contractor procurement, and review and approval of site radiological survey clearance program and IV program plans. The implementation phase assures sufficient collection of data and adjustment can be made as necessary. During the assessment phase, data and reports are reviewed, and coordination between the clean-up contractor and IV personnel occurs. A decision of whether or not the clearance actions need to be altered occurs during the decision phase. IV certification that the clearance action is complete also occurs during the decision phase. Issuance of the final guidance document is meant to improve the consultation process and minimize iterative interactions with EM. Standardizing and formalizing property clearance and IV actions also helps assure that cleanup actions are transparent and sustainable for future scrutiny. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Influence of polarization and a source model for dose calculation in MRT

Purpose: Microbeam Radiation Therapy (MRT), an alternative preclinical treatment strategy using spatially modulated synchrotron radiation on a micrometer scale, has the great potential to cure malignant tumors (e.g., brain tumors) while having low side effects on normal tissue. Dose measurement and calculation in MRT is challenging because of the spatial accuracy required and the arising high dose differences. Dose calculation with Monte Carlo simulations is time consuming and their accuracy is still a matter of debate. In particular, the influence of photon polarization has been discussed in the literature. Moreover, it is controversial whether a complete knowledge of phase space trajectories, i.e., the simulation of the machine from the wiggler to the collimator, is necessary in order to accurately calculate the dose. Methods: With Monte Carlo simulations in the Geant4 toolkit, the authors investigate the influence of polarization on the dose distribution and the therapeutically important peak to valley dose ratios (PVDRs). Furthermore, the authors analyze in detail phase space information provided byMartínez-Rovira et al. [“Development and commissioning of a Monte Carlo photon model for the forthcoming clinical trials in microbeam radiation therapy,” Med. Phys. 39(1), 119–131 (2012)] and examine its influence on peak and valley doses. A simple source model is developed using parallel beams and its applicability is shown in a semiadjoint Monte Carlo simulation. Results are compared to measurements and previously published data. Results: Polarization has a significant influence on the scattered dose outside the microbeam field. In the radiation field, however, dose and PVDRs deduced from calculations without polarization and with polarization differ by less than 3%. The authors show that the key consequences from the phase space information for dose calculations are inhomogeneous primary photon flux, partial absorption due to inclined beam incidence outside the field center, increased beam width and center to center distance due to the beam propagation from the collimator to the phantom surface and imperfect absorption in the absorber material of the Multislit Collimator. These corrections have an effect of approximately 10% on the valley dose and suffice to describe doses in MRT within the measurement uncertainties of currently available dosimetry techniques. Conclusions: The source for the first clinical pet trials in MRT is characterized with respect to its phase space and the photon polarization. The results suggest the use of a presented simplified phase space model in dose calculations and hence pave the way for alternative and fast dose calculation algorithms. They also show that the polarization is of minor importance for the clinical important peak and valley doses inside the microbeam field.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2011

Analysis of Prostate Patient Setup and Tracking Data: Potential Intervention Strategies

Purpose: To evaluate the setup, interfraction, and intrafraction organ motion error distributions and simulate intrafraction intervention strategies for prostate radiotherapy. Methods and Materials: A total of 17 patients underwent treatment setup and were monitored using the Calypso system during radiotherapy. On average, the prostate tracking measurements were performed for 8 min/fraction for 28 fractions for each patient. For both patient couch shift data and intrafraction organ motion data, the systematic and random errors were obtained from the patient population. The planning target volume margins were calculated using the van Herk formula. Two intervention strategies were simulated using the tracking data: the deviation threshold and period. The related planning target volume margins, time costs, and prostate position 'fluctuation' were presented. Results: The required treatment margin for the left-right, superoinferior, and anteroposterior axes was 8.4, 10.8, and 14.7 mm for skin mark-only setup and 1.3, 2.3, and 2.8 mm using the on-line setup correction, respectively. Prostate motion significantly correlated among the superoinferior and anteroposterior directions. Of the 17 patients, 14 had prostate motion within 5 mm of the initial setup position for {>=}91.6% of the total tracking time. The treatment margin decreased to 1.1, 1.8, and 2.3 mm with a 3-mm threshold correction and to 0.5, 1.0, and 1.5 mm with an every-2-min correction in the left-right, superoinferior, and anteroposterior directions, respectively. The periodic corrections significantly increase the treatment time and increased the number of instances when the setup correction was made during transient excursions. Conclusions: The residual systematic and random error due to intrafraction prostate motion is small after on-line setup correction. Threshold-based and time-based intervention strategies both reduced the planning target volume margins. The time-based strategies increased the treatment time and the in-fraction position fluctuation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Calculation of the Phenix end-of-life test in natural circulation with the CATHARE code

The Inst. of Radiological Protection and Nuclear Safety (IRSN) acts as technical support to French public authorities. As such, IRSN is in charge of safety assessment of operating and under construction reactors, as well as future projects. In this framework, one current objective of IRSN is to evaluate the ability and accuracy of numerical tools to foresee consequences of accidents. One of the advantages pointed up for fast reactors cooled by heavy liquid metal is the possibility of decay heat removal based on natural convection. The promotion of this passive cooling mode in future safety demonstrations will involve the use of adapted and validated numerical codes. After the final shutdown of the Phenix sodium cooled fast reactor in 2009, a set of tests covering different areas was conducted for code validation, including a natural circulation test in the primary circuit. Experimental data were issued by CEA to organize a benchmark exercise in the frame of an IAEA Coordinated Research Project (CRP), with the objective to assess the system-codes capability in simulating the thermal-hydraulics behavior of sodium cooled fast reactors in such accidental conditions. IRSN participated to this benchmark with the CATHARE code. This code, co-developed by CEA, EDF, AREVA and IRSN and widely used for PWR safety studies, was recently extended for sodium applications. This paper presents the CATHARE modeling of the Phenix primary circuit and the results obtained. A relatively good agreement was found with available measurements considering the experimental uncertainties. This work stressed the local aspects of phenomena occurring during the natural convection establishment and the limits of a 0D/1D approach. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗