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

Panoramic cone beam computed tomography

Purpose: Cone-beam computed tomography (CBCT) is the main imaging tool for image-guided radiotherapy but its functionality is limited by a small imaging volume and restricted image position (imaged at the central instead of the treatment position for peripheral lesions to avoid collisions). In this paper, the authors present the concept of ''panoramic CBCT,'' which can image patients at the treatment position with an imaging volume as large as practically needed. Methods: In this novel panoramic CBCT technique, the target is scanned sequentially from multiple view angles. For each view angle, a half scan (180 deg. + {theta}{sub cone} where {theta}{sub cone} is the cone angle) is performed with the imaging panel positioned in any location along the beam path. The panoramic projection images of all views for the same gantry angle are then stitched together with the direct image stitching method (i.e., according to the reported imaging position) and full-fan, half-scan CBCT reconstruction is performed using the stitched projection images. To validate this imaging technique, the authors simulated cone-beam projection images of the Mathematical Cardiac Torso (MCAT) thorax phantom for three panoramic views. Gaps, repeated/missing columns, and different exposure levels were introduced between adjacent views to simulate imperfect image stitching due to uncertainties in imaging position or output fluctuation. A modified simultaneous algebraic reconstruction technique (modified SART) was developed to reconstruct CBCT images directly from the stitched projection images. As a gold standard, full-fan, full-scan (360 deg. gantry rotation) CBCT reconstructions were also performed using projection images of one imaging panel large enough to encompass the target. Contrast-to-noise ratio (CNR) and geometric distortion were evaluated to quantify the quality of reconstructed images. Monte Carlo simulations were performed to evaluate the effect of scattering on the image quality and imaging dose for both standard and panoramic CBCT. Results: Truncated images with artifacts were observed for the CBCT reconstruction using projection images of the central view only. When the image stitching was perfect, complete reconstruction was obtained for the panoramic CBCT using the modified SART with the image quality similar to the gold standard (full-scan, full-fan CBCT using one large imaging panel). Imperfect image stitching, on the other hand, lead to (streak, line, or ring) reconstruction artifacts, reduced CNR, and/or distorted geometry. Results from Monte Carlo simulations showed that, for identical imaging quality, the imaging dose was lower for the panoramic CBCT than that acquired with one large imaging panel. For the same imaging dose, the CNR of the three-view panoramic CBCT was 50% higher than that of the regular CBCT using one big panel. Conclusions: The authors have developed a panoramic CBCT technique and demonstrated with simulation data that it can image tumors of any location for patients of any size at the treatment position with comparable or less imaging dose and time. However, the image quality of this CBCT technique is sensitive to the reconstruction artifacts caused by imperfect image stitching. Better algorithms are therefore needed to improve the accuracy of image stitching for panoramic CBCT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Simultaneous reduction of radiation dose and scatter for CBCT by using collimators

Purpose: On-board cone-beam CT (CBCT) imaging has been widely available in radiotherapy clinic for target localization. However, the extra radiation dose from CBCT is always a concern for its frequent use. Additionally, the relatively large scatter in CBCT often degrades the image quality. By using collimators, some of the X-rays can be stopped from reaching the patient and the detectors, hence both the scatter and the patient doses are simultaneously reduced. The authors show in this work that the collimated CBCT data can be reconstructed without any noticeable artifacts for certain collimator blocking ratios and blocking patterns, and the focus of this work is to study the relationship between the image quality and these two collimator factors.Methods: A CBCT system with collimators was simulated following the typical geometry used in clinic. Different collimator designs were tested by varying the size and the number of the collimator slits, and at the same time, the ratio of transmitted beams to total beams was varied from 100% to 10%, resulting in hundreds of different simulation scenarios. Lung and pelvis phantoms created from patients CT images were used in the simulations, and an iterative reconstruction algorithm using the compressed sensing technique was adopted. The image quality was examined by root mean square errors (RMSEs) and compared with the conventional CBCT images.Results: The CBCT image quality increases as the amount of beams passing through the collimators increases, and decreases as the size of the collimator slits increases. With ultra-high resolution collimators, the RMSEs were comparable to the conventional CBCT image quality until the beam transmission rate is reduced below 25%.Conclusions: Collimators can reduce the scatters and radiation dose, however, the collimated CBCT image quality is strongly dependent on both the collimator blocking ratio and the blocking pattern. To achieve image quality comparable to the conventional CBCT, the amount of information and data format must be adequate.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Customized multi-group cross section generation with NJOY for discrete ordinates X-ray computed tomography and radiography simulation - 25589

The purpose of this work was to explore the creation of photo-atomic multi-group cross section libraries to be used with a software package DOCTORS (Discrete Ordinates Computed Tomography and Radiography Simulator). This software solves the linear Boltzmann equation using the discrete ordinates method. To create these libraries, NJOY2016 was used, creating both fine and broad energy multi-group cross section files. The cross section's accuracy was tested against an equivalent Monte Carlo simulation using MCNP6. Two simulation geometries were used. The first, a cylindrical water phantom with a single source projection placed in front, simulating an X-ray radiography. The second used the same water phantom with 16 cone beam sources placed evenly around the phantom to simulate a computed tomography (CT) scan. The accuracy of the fine energy multi-group cross sections when used in DOCTORS has been verified. The accuracy of the broad energy multigroup cross sections when used in DOCTORS shows that more work must be done before they can be used reliably. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Experimental validation of the TOPAS Monte Carlo system for passive scattering proton therapy

Purpose: TOPAS (TOol for PArticle Simulation) is a particle simulation code recently developed with the specific aim of making Monte Carlo simulations user-friendly for research and clinical physicists in the particle therapy community. The authors present a thorough and extensive experimental validation of Monte Carlo simulations performed with TOPAS in a variety of setups relevant for proton therapy applications. The set of validation measurements performed in this work represents an overall end-to-end testing strategy recommended for all clinical centers planning to rely on TOPAS for quality assurance or patient dose calculation and, more generally, for all the institutions using passive-scattering proton therapy systems. Methods: The authors systematically compared TOPAS simulations with measurements that are performed routinely within the quality assurance (QA) program in our institution as well as experiments specifically designed for this validation study. First, the authors compared TOPAS simulations with measurements of depth-dose curves for spread-out Bragg peak (SOBP) fields. Second, absolute dosimetry simulations were benchmarked against measured machine output factors (OFs). Third, the authors simulated and measured 2D dose profiles and analyzed the differences in terms of field flatness and symmetry and usable field size. Fourth, the authors designed a simple experiment using a half-beam shifter to assess the effects of multiple Coulomb scattering, beam divergence, and inverse square attenuation on lateral and longitudinal dose profiles measured and simulated in a water phantom. Fifth, TOPAS’ capabilities to simulate time dependent beam delivery was benchmarked against dose rate functions (i.e., dose per unit time vs time) measured at different depths inside an SOBP field. Sixth, simulations of the charge deposited by protons fully stopping in two different types of multilayer Faraday cups (MLFCs) were compared with measurements to benchmark the nuclear interaction models used in the simulations. Results: SOBPs’ range and modulation width were reproduced, on average, with an accuracy of +1, −2 and ±3 mm, respectively. OF simulations reproduced measured data within ±3%. Simulated 2D dose-profiles show field flatness and average field radius within ±3% of measured profiles. The field symmetry resulted, on average in ±3% agreement with commissioned profiles. TOPAS accuracy in reproducing measured dose profiles downstream the half beam shifter is better than 2%. Dose rate function simulation reproduced the measurements within ∼2% showing that the four-dimensional modeling of the passively modulation system was implement correctly and millimeter accuracy can be achieved in reproducing measured data. For MLFCs simulations, 2% agreement was found between TOPAS and both sets of experimental measurements. The overall results show that TOPAS simulations are within the clinical accepted tolerances for all QA measurements performed at our institution. Conclusions: Our Monte Carlo simulations reproduced accurately the experimental data acquired through all the measurements performed in this study. Thus, TOPAS can reliably be applied to quality assurance for proton therapy and also as an input for commissioning of commercial treatment planning systems. This work also provides the basis for routine clinical dose calculations in patients for all passive scattering proton therapy centers using TOPAS.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Monte Carlo simulation of photon buildup factors for shielding materials in radiotherapy x-ray facilities

Purpose: This paper presents the results of a series of calculations to determine buildup factors for ordinary concrete, baryte concrete, lead, steel, and iron in broad beam geometry for photons energies from 0.125 to 25.125 MeV at 0.250 MeV intervals.Methods: Monte Carlo N-particle radiation transport computer code has been used to determine the buildup factors for the studied shielding materials.Results: The computation of the primary broad beams using buildup factors data was done for nine published megavoltage photon beam spectra ranging from 4 to 25 MV in nominal energies, representing linacs made by the three major manufacturers. The first tenth value layer and the equilibrium tenth value layer are calculated from the broad beam transmission for these nine primary megavoltage photon beam spectra.Conclusions: The results, compared with published data, show the ability of these buildup factor data to predict shielding transmission curves for the primary radiation beam. Therefore, the buildup factor data can be combined with primary, scatter, and leakage x-ray spectra to perform computation of broad beam transmission for barriers in radiotherapy shielding x-ray facilities.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2011

Radionuclide Retention in Concrete Wasteforms

Assessing long-term performance of Category 3 waste cement grouts for radionuclide encasement requires knowledge of the radionuclide-cement interactions and mechanisms of retention (i.e., sorption or precipitation); the mechanism of contaminant release; the significance of contaminant release pathways; how wasteform performance is affected by the full range of environmental conditions within the disposal facility; the process of wasteform aging under conditions that are representative of processes occurring in response to changing environmental conditions within the disposal facility; the effect of wasteform aging on chemical, physical, and radiological properties; and the associated impact on contaminant release. This knowledge will enable accurate prediction of radionuclide fate when the wasteforms come in contact with groundwater. The information present in the report provides data that (1) measures the effect of concrete wasteform properties likely to influence radionuclide migration; and (2) quantifies the rate of carbonation of concrete materials in a simulated vadose zone repository.

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

Assessment of Uncertainties in Treatment Planning for Scanned Ion Beam Therapy of Moving Tumors

Purpose: To provide methods for quantification of uncertainties in 4-dimensional (4D) treatment during treatment planning. Methods and Materials: Uncertainty information was generated by multiple 4D treatment simulations with varying parameters. Sampled data were analyzed using uncertainty visualization methods that have been added to common treatment plan evaluation methods (eg, dose-volume histogram and dose distribution analysis). To illustrate the potential of the introduced methods, uncertainty analysis was completed for a single lung cancer case using 3 motion mitigation techniques: gating, slice-by-slice rescanning, and breath-controlled rescanning. Results: By repeating 4D dose calculations with varying parameters, we were able to show local uncertainties in dose distributions and to evaluate the stability of treatment setups. The new methods were found suitable for uncertainty evaluation in 4D treatment planning of moving tumors. Calculation time of the uncertainty base data was time consuming but contrivable overnight. Conclusions: Uncertainty analysis and visualization for 4D treatment planning provide an important tool in the decision process for an optimal treatment approach.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Spatial frequency spectrum of the x-ray scatter distribution in CBCT projections

Purpose: X-ray scatter is a source of significant image quality loss in cone-beam computed tomography (CBCT). The use of Monte Carlo (MC) simulations separating primary and scattered photons has allowed the structure and nature of the scatter distribution in CBCT to become better elucidated. This work seeks to quantify the structure and determine a suitable basis function for the scatter distribution by examining its spectral components using Fourier analysis.Methods: The scatter distribution projection data were simulated using a CBCT MC model based on the EGSnrc code. CBCT projection data, with separated primary and scatter signal, were generated for a 30.6 cm diameter water cylinder [single angle projection with varying axis-to-detector distance (ADD) and bowtie filters] and two anthropomorphic phantoms (head and pelvis, 360 projections sampled every 1°, with and without a compensator). The Fourier transform of the resulting scatter distributions was computed and analyzed both qualitatively and quantitatively. A novel metric called the scatter frequency width (SFW) is introduced to determine the scatter distribution's frequency content. The frequency content results are used to determine a set basis functions, consisting of low-frequency sine and cosine functions, to fit and denoise the scatter distribution generated from MC simulations using a reduced number of photons and projections. The signal recovery is implemented using Fourier filtering (low-pass Butterworth filter) and interpolation. Estimates of the scatter distribution are used to correct and reconstruct simulated projections.Results: The spatial and angular frequencies are contained within a maximum frequency of 0.1 cm{sup −1} and 7/(2π) rad{sup −1} for the imaging scenarios examined, with these values varying depending on the object and imaging setup (e.g., ADD and compensator). These data indicate spatial and angular sampling every 5 cm and π/7 rad (∼25°) can be used to properly capture the scatter distribution, with reduced sampling possible depending on the imaging scenario. Using a low-pass Butterworth filter, tuned with the SFW values, to denoise the scatter projection data generated from MC simulations using 10{sup 6} photons resulted in an error reduction of greater than 85% for the estimating scatter in single and multiple projections. Analysis showed that the use of a compensator helped reduce the error in estimating the scatter distribution from limited photon simulations by more than 37% when compared to the case without a compensator for the head and pelvis phantoms. Reconstructions of simulated head phantom projections corrected by the filtered and interpolated scatter estimates showed improvements in overall image quality.Conclusions: The spatial frequency content of the scatter distribution in CBCT is found to be contained within the low frequency domain. The frequency content is modulated both by object and imaging parameters (ADD and compensator). The low-frequency nature of the scatter distribution allows for a limited set of sine and cosine basis functions to be used to accurately represent the scatter signal in the presence of noise and reduced data sampling decreasing MC based scatter estimation time. Compensator induced modulation of the scatter distribution reduces the frequency content and improves the fitting results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Optimization Studies of Radiation Shielding for PIP-II Project at Fermilab

The Proton Improvement Plan-II (PIP-II) at Fermilab represents a significant advancement in the quest to answer some of the most profound questions about our universe using the world's most intense high-energy neutrino beam. The project requires the construction of a new addition to the Fermilab accelerator complex – an 800-MeV high-intensity superconducting linear accelerator. Ensuring the safety and regulatory compliance of this ambitious project is paramount, necessitating thorough dose rate assessments under both normal operational and accidental scenarios to align with the Fermilab Radiological Control Manual (FRCM) standards. Our approach included a shielding optimization used for the simulations with the Monte Carlo code MARS [1,2,3] to incorporate new magnet and collimator designs, essential for reflecting the current state of PIP-II infrastructure. The implementation of high-resolution detector planes, despite their computational demands, enabled us to gather detailed radiation field data crucial for optimizing shielding configurations. To overcome the significant computational demands, we developed a branching code that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through geometry splitting and the application of Russian Roulette techniques, tailored to prioritize regions of interest based on predefined importances and weight limits.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Optimization Studies of Radiation Shielding for PIP-II Project at Fermilab

The Proton Improvement Plan-II (PIP-II) at Fermilab represents a significant advancement in the quest to answer some of the most profound questions about our universe using the world's most intense high-energy neutrino beam. The project requires the construction of a new addition to the Fermilab accelerator complex – an 800-MeV high-intensity superconducting linear accelerator. Ensuring the safety and regulatory compliance of this ambitious project is paramount, necessitating thorough dose rate assessments under both normal operational and accidental scenarios to align with the Fermilab Radiological Control Manual (FRCM) standards. Our approach included a shielding optimization used for the simulations with the Monte Carlo code MARS [1,2,3] to incorporate new magnet and collimator designs, essential for reflecting the current state of PIP-II infrastructure. The implementation of high-resolution detector planes, despite their computational demands, enabled us to gather detailed radiation field data crucial for optimizing shielding configurations. To overcome the significant computational demands, we developed a branching code that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through geometry splitting and the application of Russian Roulette techniques, tailored to prioritize regions of interest based on predefined importances and weight limits.

Makovec, Alajos↗
DOE OSTI2013

Improving best-phase image quality in cardiac CT by motion correction with MAM optimization

Purpose: Research in image reconstruction for cardiac CT aims at using motion correction algorithms to improve the image quality of the coronary arteries. The key to those algorithms is motion estimation, which is currently based on 3-D/3-D registration to align the structures of interest in images acquired in multiple heart phases. The need for an extended scan data range covering several heart phases is critical in terms of radiation dose to the patient and limits the clinical potential of the method. Furthermore, literature reports only slight quality improvements of the motion corrected images when compared to the most quiet phase (best-phase) that was actually used for motion estimation. In this paper a motion estimation algorithm is proposed which does not require an extended scan range but works with a short scan data interval, and which markedly improves the best-phase image quality. Methods: Motion estimation is based on the definition of motion artifact metrics (MAM) to quantify motion artifacts in a 3-D reconstructed image volume. The authors use two different MAMs, entropy, and positivity. By adjusting the motion field parameters, the MAM of the resulting motion-compensated reconstruction is optimized using a gradient descent procedure. In this way motion artifacts are minimized. For a fast and practical implementation, only analytical methods are used for motion estimation and compensation. Both the MAM-optimization and a 3-D/3-D registration-based motion estimation algorithm were investigated by means of a computer-simulated vessel with a cardiac motion profile. Image quality was evaluated using normalized cross-correlation (NCC) with the ground truth template and root-mean-square deviation (RMSD). Four coronary CT angiography patient cases were reconstructed to evaluate the clinical performance of the proposed method. Results: For the MAM-approach, the best-phase image quality could be improved for all investigated heart phases, with a maximum improvement of the NCC value by 100% and of the RMSD value by 81%. The corresponding maximum improvements for the registration-based approach were 20% and 40%. In phases with very rapid motion the registration-based algorithm obtained better image quality, while the image quality of the MAM algorithm was superior in phases with less motion. The image quality improvement of the MAM optimization was visually confirmed for the different clinical cases. Conclusions: The proposed method allows a software-based best-phase image quality improvement in coronary CT angiography. A short scan data interval at the target heart phase is sufficient, no additional scan data in other cardiac phases are required. The algorithm is therefore directly applicable to any standard cardiac CT acquisition protocol.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Individualized Positron Emission Tomography–Based Isotoxic Accelerated Radiation Therapy Is Cost-Effective Compared With Conventional Radiation Therapy: A Model-Based Evaluation

Purpose: To evaluate long-term health effects, costs, and cost-effectiveness of positron emission tomography (PET)-based isotoxic accelerated radiation therapy treatment (PET-ART) compared with conventional fixed-dose CT-based radiation therapy treatment (CRT) in non-small cell lung cancer (NSCLC). Methods and Materials: Our analysis uses a validated decision model, based on data of 200 NSCLC patients with inoperable stage I-IIIB. Clinical outcomes, resource use, costs, and utilities were obtained from the Maastro Clinic and the literature. Primary model outcomes were the difference in life-years (LYs), quality-adjusted life-years (QALYs), costs, and the incremental cost-effectiveness and cost/utility ratio (ICER and ICUR) of PET-ART versus CRT. Model outcomes were obtained from averaging the predictions for 50,000 simulated patients. A probabilistic sensitivity analysis and scenario analyses were carried out. Results: The average incremental costs per patient of PET-ART were €569 (95% confidence interval [CI] €−5327-€6936) for 0.42 incremental LYs (95% CI 0.19-0.61) and 0.33 QALYs gained (95% CI 0.13-0.49). The base-case scenario resulted in an ICER of €1360 per LY gained and an ICUR of €1744 per QALY gained. The probabilistic analysis gave a 36% probability that PET-ART improves health outcomes at reduced costs and a 64% probability that PET-ART is more effective at slightly higher costs. Conclusion: On the basis of the available data, individualized PET-ART for NSCLC seems to be cost-effective compared with CRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Dose calculation software for helical tomotherapy, utilizing patient CT data to calculate an independent three-dimensional dose cube

Purpose: Treatment plans for the TomoTherapy unit are produced with a planning system that is integral to the unit. The authors have produced an independent dose calculation system, to enable plans to be recalculated in three dimensions, using the patient's CT data. Methods: Software has been written using MATLAB. The DICOM-RT plan object is used to determine the treatment parameters used, including the treatment sinogram. Each projection of the sinogram is segmented and used to calculate dose at multiple calculation points in a three-dimensional grid using tables of measured beam data. A fast ray-trace algorithm is used to determine effective depth for each projection angle at each calculation point. Calculations were performed on a standard desktop personal computer, with a 2.6 GHz Pentium, running Windows XP. Results: The time to perform a calculation, for 3375 points averaged 1 min 23 s for prostate plans and 3 min 40 s for head and neck plans. The mean dose within the 50% isodose was calculated and compared with the predictions of the TomoTherapy planning system. When the modified CT (which includes the TomoTherapy couch) was used, the mean difference for ten prostate patients, was -0.4% (range -0.9% to +0.3%). With the original CT (which included the CT couch), the mean difference was -1.0% (range -1.7% to 0.0%). The number of points agreeing with a gamma 3%/3 mm averaged 99.2% with the modified CT, 96.3% with the original CT. For ten head and neck patients, for the modified and original CT, respectively, the mean difference was +1.1% (range -0.4% to +3.1%) and 1.1% (range -0.4% to +3.0%) with 94.4% and 95.4% passing a gamma 4%/4 mm. The ability of the program to detect a variety of simulated errors has been tested. Conclusions: By using the patient's CT data, the independent dose calculation performs checks that are not performed by a measurement in a cylindrical phantom. This enables it to be used either as an additional check or to replace phantom measurements for some patients. The software has potential to be used in any application where one wishes to model changes to patient conditions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Using Generalized Equivalent Uniform Dose Atlases to Combine and Analyze Prospective Dosimetric and Radiation Pneumonitis Data From 2 Non-Small Cell Lung Cancer Dose Escalation Protocols

Purpose: To demonstrate the use of generalized equivalent uniform dose (gEUD) atlas for data pooling in radiation pneumonitis (RP) modeling, to determine the dependence of RP on gEUD, to study the consistency between data sets, and to verify the increased statistical power of the combination. Methods and Materials: Patients enrolled in prospective phase I/II dose escalation studies of radiation therapy of non-small cell lung cancer at Memorial Sloan-Kettering Cancer Center (MSKCC) (78 pts) and the Netherlands Cancer Institute (NKI) (86 pts) were included; 10 (13%) and 14 (17%) experienced RP requiring steroids (RPS) within 6 months after treatment. gEUD was calculated from dose-volume histograms. Atlases for each data set were created using 1-Gy steps from exact gEUDs and RPS data. The Lyman-Kutcher-Burman model was fit to the atlas and exact gEUD data. Heterogeneity and inconsistency statistics for the fitted parameters were computed. gEUD maps of the probability of RPS rate {>=}20% were plotted. Results: The 2 data sets were homogeneous and consistent. The best fit values of the volume effect parameter a were small, with upper 95% confidence limit around 1.0 in the joint data. The likelihood profiles around the best fit a values were flat in all cases, making determination of the best fit a weak. All confidence intervals (CIs) were narrower in the joint than in the individual data sets. The minimum P value for correlations of gEUD with RPS in the joint data was .002, compared with P=.01 and .05 for MSKCC and NKI data sets, respectively. gEUD maps showed that at small a, RPS risk increases with gEUD. Conclusions: The atlas can be used to combine gEUD and RPS information from different institutions and model gEUD dependence of RPS. RPS has a large volume effect with the mean dose model barely included in the 95% CI. Data pooling increased statistical power.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Response of the REWARD detection system to the presence of a Radiological Dispersal Device

In recent years an increased international concern has emerged about the radiological and nuclear (RN) threats associated with the illicit trafficking of nuclear and radioactive materials that could be potentially used for terrorist attacks. The objective of the REWARD (Real Time Wide Area Radiation Surveillance System) project, co-funded by the European Union 7. Framework Programme Security, consisted in building a mobile system for real time, wide area radiation surveillance, using a CdZnTe detector for gamma radiation and a neutron detector based on novel silicon technologies. The sensing unit includes a GPS system and a wireless communication interface to send the data remotely to a monitoring base station, where it will be analyzed in real time and correlated with historical data from the tag location, in order to generate an alarm when an abnormal situation is detected. Due to its portability and accuracy, the system will be extremely useful in many different scenarios such as nuclear terrorism, lost radioactive sources, radioactive contamination or nuclear accidents. This paper shortly introduces the REWARD detection system, depicts some terrorist threat scenarios involving radioactive sources and special nuclear materials and summarizes the simulation work undertaken during the past three years in the framework of the REWARD project. The main objective consisted in making predictions regarding the behavior of the REWARD system in the presence of a Radiological Dispersion Device (RDD), one of the reference scenarios foreseen for REWARD, using the Monte Carlo simulation program MCNP6. The reference scenario is characterized in detail, from the i) radiological protection, ii) radiation detection requirements and iii) communications points of view. Experimental tests were performed at the Fire Brigades Facilities in Rome and at the Naples Fire Brigades, and the results, which validate the simulation work, are presented and analyzed. The response of the REWARD detection system to the presence of an RDD is predicted and discussed. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Spectral analysis of mammographic images using a multitaper method

Purpose: Power spectral analysis in radiographic images is conventionally performed using a windowed overlapping averaging periodogram. This study describes an alternative approach using a multitaper technique and compares its performance with that of the standard method. This tool will be valuable in power spectrum estimation of images, whose content deviates significantly from uniform white noise. The performance of the multitaper approach will be evaluated in terms of spectral stability, variance reduction, bias, and frequency precision. The ultimate goal is the development of a useful tool for image quality assurance. Methods: A multitaper approach uses successive data windows of increasing order. This mitigates spectral leakage allowing one to calculate a reduced-variance power spectrum. The multitaper approach will be compared with the conventional power spectrum method in several typical situations, including the noise power spectra (NPS) measurements of simulated projection images of a uniform phantom, NPS measurement of real detector images of a uniform phantom for two clinical digital mammography systems, and the estimation of the anatomic noise in mammographic images (simulated images and clinical mammograms). Results: Examination of spectrum variance versus frequency resolution and bias indicates that the multitaper approach is superior to the conventional single taper methods in the prevention of spectrum leakage and variance reduction. More than four times finer frequency precision can be achieved with equivalent or less variance and bias. Conclusions: Without any shortening of the image data length, the bias is smaller and the frequency resolution is higher with the multitaper method, and the need to compromise in the choice of regions of interest size to balance between the reduction of variance and the loss of frequency resolution is largely eliminated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

A six-channel pediatric coil array for detection of children spinal pathologies by MRI at 1.5 Tesla

Nowadays, magnetic resonance (MR) in Mexico has become a standard technique for clinical imaging. Although most of the times the MR systems contain only coils oriented for adults. Radiologists use these coils for children studies due to the non-availability of pediatric coils. Image quality is decreased due to the low signal to noise ratio delivered to the system. The development of RF coils is always focused towards increasing SNR and optimizing the RF penetration into the sample. Moreover, spinal pathologies in children, which are an important topic in pediatric care, cover congenital and neuromuscular disorders that occur in childhood. In this work, the design of a dedicated six-channel coil for detection of spinal pathologies at 1.5 Tesla is addressed. Numerical electromagnetic simulations were performed in order to evaluate their magnetic field performance at (63.6 MHz) 1.5 Tesla. The magnetic field uniformity as well as the RF penetration depth of the coil configurations was evaluated in order to find the best/optimized coil array configuration. The coil is comprised of three rows, one with 4 coil elements and two with only one coil element. Phantom and in vivo images were acquired with the six-channel pediatric coil array. The phantom images agree with the simulated data. In vivo images acquired with the 6-channel pediatric coil array have shown very good penetration depth and homogeneity, which allow better image quality throughout the whole FOV. In addition, the parallel imaging capabilities of the array allow the acceleration of the experiments avoiding possible motion artifacts.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Novel Assessment of Renal Motion in Children as Measured via Four-Dimensional Computed Tomography

Objectives: Abdominal intensity-modulated radiation therapy and proton therapy require quantification of target and organ motion to optimize localization and treatment. Although addressed in adults, there is no available literature on this issue in pediatric patients. We assessed physiologic renal motion in pediatric patients. Methods and Materials: Twenty free-breathing pediatric patients at a median age of 8 years (range, 2-18 years) with intra-abdominal tumors underwent computed tomography simulation and four-dimensional computed tomography acquisition (slice thickness, 3 mm). Kidneys and diaphragms were contoured during eight phases of respiration to estimate center-of-mass motion. We quantified center of kidney mass mobility vectors in three dimensions: anteroposterior (AP), mediolateral (ML), and superoinferior (SI). Results: Kidney motion decreases linearly with decreasing age and height. The 95% confidence interval for the averaged minima and maxima of renal motion in children younger than 9 years was 5-9 mm in the ML direction, 4-11 mm in the AP direction, and 12-25 mm in the SI dimension for both kidneys. In children older than 9 years, the same confidence interval reveals a widening range of motion that was 5-16 mm in the ML direction, 6-17 mm in the AP direction, and 21-52 mm in the SI direction. Although not statistically significant, renal motion correlated with diaphragm motion in older patients. The correlation between diaphragm motion and body mass index was borderline (r = 0.52, p = 0.0816) in younger patients. Conclusions: Renal motion is age and height dependent. Measuring diaphragmatic motion alone does not reliably quantify pediatric renal motion. Renal motion in young children ranges from 5 to 25 mm in orientation-specific directions. The vectors of motion range from 5 to 52 mm in older children. These preliminary data represent novel analyses of pediatric intra-abdominal organ motion.

62 RADIOLOGY AND NUCLEAR MEDICINE↗