Search NASA⌕ Search

SEARCH · Search NASA

Results for “Image coding”

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

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

At least 1,027 records · Page 57

Development of a single-shot-imaging thin film for an online Thomson parabola spectrometer

A single-shot-imaging thin scintillator film was developed for an online Thomson parabola (TP) spectrometer and the first analysis of laser accelerated ions, using the online TP spectrometer, was demonstrated at the JAEA-Kansai Advanced Relativistic Engineering Laser System (J-KAREN). An energy spectrum of {approx}4.0 MeV protons is obtained using only this imaging film without the need of a microchannel plate that is typically utilized in online ion analyses. A general-purpose Monte Carlo particle and heavy ion-transport code system, which consists of various quantum dynamics models, was used for the prediction of the luminescent properties of the scintillator. The simulation can reasonably predict not only the ion trajectories detected by the spectrometer, but also luminescence properties.

36 MATERIALS SCIENCE↗

SU-E-T-29: A Web Application for GPU-Based Monte Carlo IMRT/VMAT QA with Delivered Dose Verification

Purpose: To enable an existing web application for GPU-based Monte Carlo (MC) 3D dosimetry quality assurance (QA) to compute “delivered dose” from linac logfile data. Methods: We added significant features to an IMRT/VMAT QA web application which is based on existing technologies (HTML5, Python, and Django). This tool interfaces with python, c-code libraries, and command line-based GPU applications to perform a MC-based IMRT/VMAT QA. The web app automates many complicated aspects of interfacing clinical DICOM and logfile data with cutting-edge GPU software to run a MC dose calculation. The resultant web app is powerful, easy to use, and is able to re-compute both plan dose (from DICOM data) and delivered dose (from logfile data). Both dynalog and trajectorylog file formats are supported. Users upload zipped DICOM RP, CT, and RD data and set the expected statistic uncertainty for the MC dose calculation. A 3D gamma index map, 3D dose distribution, gamma histogram, dosimetric statistics, and DVH curves are displayed to the user. Additional the user may upload the delivery logfile data from the linac to compute a 'delivered dose' calculation and corresponding gamma tests. A comprehensive PDF QA report summarizing the results can also be downloaded. Results: We successfully improved a web app for a GPU-based QA tool that consists of logfile parcing, fluence map generation, CT image processing, GPU based MC dose calculation, gamma index calculation, and DVH calculation. The result is an IMRT and VMAT QA tool that conducts an independent dose calculation for a given treatment plan and delivery log file. The system takes both DICOM data and logfile data to compute plan dose and delivered dose respectively. Conclusion: We sucessfully improved a GPU-based MC QA tool to allow for logfile dose calculation. The high efficiency and accessibility will greatly facilitate IMRT and VMAT QA.

07 ISOTOPE AND RADIATION SOURCES↗

Development of CSG-based radiation shielding module for ARCHER: preliminary results for photons - Paper 126

ARCHER is a parallel Monte Carlo radiation transport code being developed for various hardware platforms including multi-core CPUs, the emerging Nvidia graphics processing units (GPUs) and Intel Xeon Phi coprocessors. Previous studies have demonstrated the accuracy and speed for CT imaging and radiotherapy dosimetry applications. Radiation shielding applications, however, require ARCHER to be more flexible in defining and handling geometries, in addition to the voxel phantom geometries that are commonly used in medical physics calculations. In this paper, we present the development of a general geometry module used for radiation protection and shielding problems in ARCHER. The methods of geometry design and particle transport for the GPUs are introduced. Preliminary results comparing against the MCNP code for cases involving Cs-137 0.662 MeV gammas are reported. The radiation shielding design module of ARCHER has been developed, which supports general geometry setups in the CSG format. The code was implemented using CUDA on Nvidia GPUs, and the performance was benchmarked against the MCNP code for a simple radiation shielding problem. Results showed that ARCHER worked properly and achieved impressive speedup due to highly parallel computing nature of the Nvidia GPUs compared to the MCNP code running on multi-core CPUs. Further optimization of the geometry module will be carried out using advanced geometry structures and variance reduction techniques to improve the precision and performance of ARCHER. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

A Monte Carlo investigation of low-Z target image quality generated in a linear accelerator using Varian's VirtuaLinac

Purpose: The focus of this work was the demonstration and validation of VirtuaLinac with clinical photon beams and to investigate the implementation of low-Z targets in a TrueBeam linear accelerator (Linac) using Monte Carlo modeling. Methods: VirtuaLinac, a cloud based web application utilizing Geant4 Monte Carlo code, was used to model the Linac treatment head components. Particles were propagated through the lower portion of the treatment head using BEAMnrc. Dose distributions and spectral distributions were calculated using DOSXYZnrc and BEAMdp, respectively. For validation, 6 MV flattened and flattening filter free (FFF) photon beams were generated and compared to measurement for square fields, 10 and 40 cm wide and at d{sub max} for diagonal profiles. Two low-Z targets were investigated: a 2.35 MeV carbon target and the proposed 2.50 MeV commercial imaging target for the TrueBeam platform. A 2.35 MeV carbon target was also simulated in a 2100EX Clinac using BEAMnrc. Contrast simulations were made by scoring the dose in the phosphor layer of an IDU20 aSi detector after propagating through a 4 or 20 cm thick phantom composed of water and ICRP bone. Results: Measured and modeled depth dose curves for 6 MV flattened and FFF beams agree within 1% for 98.3% of points at depths greater than 0.85 cm. Ninety three percent or greater of points analyzed for the diagonal profiles had a gamma value less than one for the criteria of 1.5 mm and 1.5%. The two low-Z target photon spectra produced in TrueBeam are harder than that from the carbon target in the Clinac. Percent dose at depth 10 cm is greater by 3.6% and 8.9%; the fraction of photons in the diagnostic energy range (25–150 keV) is lower by 10% and 28%; and contrasts are lower by factors of 1.1 and 1.4 (4 cm thick phantom) and 1.03 and 1.4 (20 cm thick phantom), for the TrueBeam 2.35 MV/carbon and commercial imaging beams, respectively. Conclusions: VirtuaLinac is a promising new tool for Monte Carlo modeling of novel target designs. A significant spectral difference is observed between the low-Z target beam on the Clinac platform and the proposed imaging beam line on TrueBeam, with the former providing greater diagnostic energy content.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Multi-solution nature of topology optimization and its application in design for additive manufacturing

The purpose of this paper is to introduce the multi-solution nature of topology optimization (TO) as a design tool for additive manufacturing (AM). The sensitivity of topologically optimized parts and manufacturing constraints to the initial starting point of the optimization process leading to structures with equivalent performance is explored. A modified bi-directional evolutionary structural optimization (BESO) code was used as the numerical approach to optimize a cantilever beam problem and reduce the mass by 50 per cent. Several optimized structures with relatively equivalent mechanical performance were generated by changing the initial starting point of the TO algorithm. These optimized structures were manufactured using fused deposition modeling (FDM). The equivalence of strain distribution in FDM parts was tested with the digital image correlation (DIC) technique and compared with that from the modified BESO code.The results confirm that TO could lead to a wide variety of non-unique solutions based on loading and manufacturability constraints. The modified BESO code was able to reduce the support structure needed to build the simple two-dimensional cantilever beam by 15 per cent while keeping the mechanical performance at the same level. Finally, the originality of this paper lies in introduction and application of the multi-solution nature of TO for AM as a design tool for optimizing structures with minimized features in the overhang condition and the need for support structures.

36 MATERIALS SCIENCE↗

A tracking system to calculate patient skin dose in real-time during neurointerventional procedures using a biplane x-ray imaging system

Purpose: Neurovascular interventional procedures using biplane fluoroscopic imaging systems can lead to increased risk of radiation-induced skin injuries. The authors developed a biplane dose tracking system (Biplane-DTS) to calculate the cumulative skin dose distribution from the frontal and lateral x-ray tubes and display it in real-time as a color-coded map on a 3D graphic of the patient for immediate feedback to the physician. The agreement of the calculated values with the dose measured on phantoms was evaluated. Methods: The Biplane-DTS consists of multiple components including 3D graphic models of the imaging system and patient, an interactive graphical user interface, a data acquisition module to collect geometry and exposure parameters, the computer graphics processing unit, and functions for determining which parts of the patient graphic skin surface are within the beam and for calculating dose. The dose is calculated to individual points on the patient graphic using premeasured calibration files of entrance skin dose per mAs including backscatter; corrections are applied for field area, distance from the focal spot and patient table and pad attenuation when appropriate. The agreement of the calculated patient skin dose and its spatial distribution with measured values was evaluated in 2D and 3D for simulated procedure conditions using a PMMA block phantom and an SK-150 head phantom, respectively. Dose values calculated by the Biplane-DTS were compared to the measurements made on the phantom surface with radiochromic film and a calibrated ionization chamber, which was also used to calibrate the DTS. The agreement with measurements was specifically evaluated with variation in kVp, gantry angle, and field size. Results: The dose tracking system that was developed is able to acquire data from the two x-ray gantries on a biplane imaging system and calculate the skin dose for each exposure pulse to those vertices of a patient graphic that are determined to be in the beam. The calculations are done in real-time with a typical graphic update time of 30 ms and an average vertex separation of 3 mm. With appropriate corrections applied, the Biplane-DTS was able to determine the entrance dose within 6% and the spatial distribution of the dose within 4% compared to the measurements with the ionization chamber and film for the SK150 head phantom. The cumulative dose for overlapping fields from both gantries showed similar agreement. Conclusions: The Biplane-DTS can provide a good estimate of the peak skin dose and cumulative skin dose distribution during biplane neurointerventional procedures. Real-time display of this information should help the physician manage patient dose to reduce the risk of radiation-induced skin injuries.

60 APPLIED LIFE SCIENCES↗

On the image of AGS 3 He 2+ $\vec{n}$ 0 in the blue

This note addresses the transport of Helion spin $\vec{n}$ 0 vector, from its periodic orientation in the AGS to RHIC Blue ring injection kicker, via the AGS extraction system and the AtR line. The goal is to investigate optimal injection energy into RHIC, in the matter of Helion spin matching, in the hypothesis of equal warm and cold snake strengths in the AGS. The study uses recently computed OPERA 3-D field maps of the AGS cold snake, including possibility of independent solenoid and helix settings (as discussed in Tech. Note C-A/AP/485), together with the machinery of the AGS and AtR models developed in the stepwise ray-tracing code Zgoubi. Computing tools and methods employed are discussed as well, in order to facilitate possible further checks or investigations. They are however similar to those used in an earlier study regarding the image in RHIC Blue and Yellow of AGS $\vec{n}$ 0 via the AtR in the case of proton beam (Tech. Note C-A/AP/502), which can be referred to for additional details.

43 PARTICLE ACCELERATORS↗

Analysis and Simulation of a Blue Energy Cycle

The mixing process of fresh water and seawater releases a significant amount of energy and is a potential source of renewable energy. The so called ‘blue energy’ or salinity-gradient energy can be harvested by a device consisting of carbon electrodes immersed in an electrolyte solution, based on the principle of capacitive double layer expansion (CDLE). In this study, we have investigated the feasibility of energy production based on the CDLE principle. Experiments and computer simulations were used to study the process. Mesoporous carbon materials, synthesized at the Oak Ridge National Laboratory, were used as electrode materials in the experiments. Neutron imaging of the blue energy cycle was conducted with cylindrical mesoporous carbon electrodes and 0.5 M lithium chloride as the electrolyte solution. For experiments conducted at 0.6 V and 0.9 V applied potential, a voltage increase of 0.061 V and 0.054 V was observed, respectively. From sequences of neutron images obtained for each step of the blue energy cycle, information on the direction and magnitude of lithium ion transport was obtained. A computer code was developed to simulate the process. Experimental data and computer simulations allowed us to predict energy production.

30 DIRECT ENERGY CONVERSION↗

Effective versus ion thermal temperatures in the Weizmann Ne Z-pinch: Modeling and stagnation physics

The difference between the ion thermal and effective temperatures is investigated through simulations of the Ne gas puff z-pinch reported by Kroupp et al. [Phys. Rev. Lett. 107, 105001 (2011)]. Calculations are performed using a 2D, radiation-magnetohydrodynamic code with Tabular Collisional-Radiative Equilibrium, namely Mach2-TCRE [Thornhill et al., Phys. Plasmas 8, 3480 (2001)]. The extensive data set of imaging and K-shell spectroscopy from the experiments provides a challenging validation test for z-pinch simulations. Synthetic visible images of the implosion phase match the observed large scale structure if the breakdown occurs at the density corresponding to the Paschen minimum. At the beginning of stagnation (−4 ns), computed plasma conditions change rapidly showing a rising electron density and a peak in the ion thermal temperature of ∼1.8 keV. This is larger than the ion thermal temperature (<400 eV) inferred from the experiment. By the time of peak K-shell power (0 ns), the calculated electron density is similar to the data and the electron and ion thermal temperatures are equilibrated, as is observed. Effective ion temperatures are obtained from calculated emission line widths accounting for thermal broadening and Doppler velocity shifts. The observed, large effective ion temperatures (∼4 keV) early in the stagnation of this Ne pinch can be explained solely as a combination of compressional ion heating and steep radial velocity gradients near the axis. Approximations in the modeling are discussed in regard to the higher ion thermal temperature and lower electron density early in the stagnation compared to the experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simple partitions of a hyperbolic plane of positive curvature

We construct special monohedral isotropic partitions with symmetries of the hyperbolic plane H of positive curvature with a simple 4-contour as a cell. An analogue of mosaic in these partitions called a tiling is introduced. Also we consider some fractal tilings. The existence of band tilings in each homological series with code (m, n) is proved. Bibliography: 14 titles.

97 MATHEMATICS AND COMPUTING↗

ULTRAVIOLET RADIATIVE TRANSFER MODELING OF NEARBY GALAXIES WITH EXTRAPLANAR DUSTS

In order to examine their relation to the host galaxy, the extraplanar dusts of six nearby galaxies are modeled, employing a three-dimensional Monte Carlo radiative transfer code. The targets are from the highly inclined galaxies that show dust-scattered ultraviolet halos, and the archival Galaxy Evolution Explorer FUV band images were fitted with the model. The observed images are generally well-reproduced by two dust layers and one light source layer, whose vertical and radial distributions have exponential profiles. We obtained several important physical parameters, such as star formation rate (SFR{sub UV}), face-on optical depth, and scale-heights. Three galaxies (NGC 891, NGC 3628, and UGC 11794) show clear evidence for the existence of an extraplanar dust layer. However, it is found that the remaining three targets (IC 5249, NGC 24, and NGC 4173) do not necessarily need a thick dust disk to model the ultraviolet (UV) halo, because its contribution is too small and the UV halo may be caused by the wing part of the GALEX point spread function. This indicates that the galaxy samples reported to have UV halos may be contaminated by galaxies with negligible extraplanar (halo) dust. The galaxies showing evidence of an extraplanar dust layer fall within a narrow range on the scatter plots between physical parameters such as SFR{sub UV} and extraplanar dust mass. Several mechanisms that could possibly produce the extraplanar dust are discussed. We also found a hint that the extraplanar dust scale-height might not be much different from the polycyclic aromatic hydrocarbon emission characteristic height.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine Learning Benchmarks for the Classification of Equivalent Circuit Models from Electrochemical Impedance Spectra

Analysis of Electrochemical Impedance Spectroscopy (EIS) data for electrochemical systems often consists of defining an Equivalent Circuit Model (ECM) using expert knowledge and then optimizing the model parameters to deconvolute various resistance, capacitive, inductive, or diffusion responses. For small data sets, this procedure can be conducted manually; however, it is not feasible to manually define a proper ECM for extensive data sets with a wide range of EIS responses. Automatic identification of an ECM would substantially accelerate the analysis of large sets of EIS data. We showcase machine learning methods to classify the ECMs of 9,300 impedance spectra provided by QuantumScape for the BatteryDEV hackathon. The best-performing approach is a gradient-boosted tree model utilizing a library to automatically generate features, followed by a random forest model using the raw spectral data. A convolutional neural network using boolean images of Nyquist representations is presented as an alternative, although it achieves a lower accuracy. We publish the data and open source the associated code. The approaches described in this article can serve as benchmarks for further studies. A key remaining challenge is the identifiability of the labels, underlined by the model performances and the comparison of misclassified spectra.

25 ENERGY STORAGE↗

Data Visualization: Augmented Reality

In recent years, there has been an increasing interest in developing new technologies for automated characterization and visualization of condition monitoring data. Augmented Reality (AR) is a technology that is being developed to improve such data visualization. Augmented reality has been defined as a technology that merges virtual and physical components in real-time, and in three dimensions. Wearable, commercially-available AR devices allow onsite engineers and technicians to perform inspection tasks with significantly more available information such as comparisons of past and present sensor and imager data, onsite data analysis and result displays, and various forms of metadata including technical drawings, previous inspection reports and maintenance histories, operation manuals, codes and standards, and holograms representing data analysis results superimposed onto the in situ monitored system.

42 ENGINEERING↗

Numerical modeling of the early light curves of type IIP supernovae

The early rise of Type IIP supernovae (SN IIP) provides important information for constraining the properties of their progenitors. This can, in turn, be compared to pre-explosion imaging constraints and stellar models to develop a more complete picture of how massive stars evolve and end their lives. Using the SuperNova Explosion Code (SNEC), we model the first 40 days of SNe IIP to better understand what constraints can be derived from their early light curves. We use two sets of red supergiant (RSG) progenitor models with zero-age main sequence masses in the range between 9 M{sub ⊙} and 20 M{sub ⊙}. We find that the early properties of the light curve depend most sensitively on the radius of the progenitor, and thus provide a relation between the g-band rise time and the radius at the time of explosion. This relation will be useful for deriving constraints on progenitors from future observations, especially in cases where detailed modeling of the entire rise is not practical. When comparing to observed rise times, the radii we find are a factor of a few larger than previous semi-analytic derivations and are generally in better agreement with what is found with current stellar evolution calculations as well as direct observations of RSGs.

79 ASTRONOMY AND ASTROPHYSICS↗

The Anatomy of an Unusual Edge-on Protoplanetary Disk. I. Dust Settling in a Cold Disk

As the earliest stage of planet formation, massive, optically thick, and gas-rich protoplanetary disks provide key insights into the physics of star and planet formation. When viewed edge-on, high-resolution images offer a unique opportunity to study both the radial and vertical structures of these disks and relate this to vertical settling, radial drift, grain growth, and changes in the midplane temperatures. In this work, we present multi-epoch Hubble Space Telescope and Keck scattered light images, and an Atacama Large Millimeter/submillimeter Array 1.3 mm continuum map for the remarkably flat edge-on protoplanetary disk SSTC2DJ163131.2–242627, a young solar-type star in ρ Ophiuchus. We model the 0.8 μm and 1.3 mm images in separate Markov Chain Monte Carlo (MCMC) runs to investigate the geometry and dust properties of the disk using the MCFOST radiative transfer code. In scattered light, we are sensitive to the smaller dust grains in the surface layers of the disk, while the submillimeter dust continuum observations probe larger grains closer to the disk midplane. An MCMC run combining both data sets using a covariance-based log-likelihood estimation was marginally successful, implying insufficient complexity in our disk model. The disk is well characterized by a flared disk model with an exponentially tapered outer edge viewed nearly edge-on, though some degree of dust settling is required to reproduce the vertically thin profile and lack of apparent flaring. A colder than expected disk midplane, evidence for dust settling, and residual radial substructures all point to a more complex radial density profile to be probed with future, higher-resolution observations.

47 OTHER INSTRUMENTATION↗

Analytic reconstruction algorithms for triple-source CT with horizontal data truncation

Purpose: This paper explores a triple-source imaging method with horizontal data truncation to enlarge the field of view (FOV) for big objects. Methods: The study is conducted by using theoretical analysis, mathematical deduction, and numerical simulations. The proposed algorithms are implemented in c + + and MATLAB. While the basic platform is constructed in MATLAB, the computationally intensive segments are coded in c + +, which are linked via a MEX interface. Results: A triple-source circular scanning configuration with horizontal data truncation is developed, where three pairs of x-ray sources and detectors are unevenly distributed on the same circle to cover the whole imaging object. For this triple-source configuration, a fan-beam filtered backprojection-type algorithm is derived for truncated full-scan projections without data rebinning. The algorithm is also extended for horizontally truncated half-scan projections and cone-beam projections in a Feldkamp-type framework. Using their method, the FOV is enlarged twofold to threefold to scan bigger objects with high speed and quality. The numerical simulation results confirm the correctness and effectiveness of the developed algorithms. Conclusions: The triple-source scanning configuration with horizontal data truncation cannot only keep most of the advantages of a traditional multisource system but also cover a larger FOV for big imaging objects. In addition, because the filtering is shift-invariant, the proposed algorithms are very fast and easily parallelized on graphic processing units.

60 APPLIED LIFE SCIENCES↗

SU-F-T-229: A Novel Method for EPID-Based In-Vivo Exit Dose Verification for Intensity Modulated Radiotherapy

Purpose: Electronic portal imaging device (EPID) can be used to acquire a two-dimensional exit dose distribution during treatment delivery, thus allowing the in-vivo verification of the dose delivery through a comparison of measured portal images to predicted portal dose images (PDI). The aim of this study was to present a novel method to easily and accurately predict PDI, and to establish an EPID-based in-vivo dose verification method during IMRT treatments. Methods: We developed a model to determine the predicted portal dose at the same plane of the EPID detector location. The Varian EPID (aS1000) positions at 150cm source-to-detector-distance (SDD), and can be used to acquire in-vivo exit dose using Portal Dosimetry (PD) function. Our model was generated to make an equivalent water thickness represent the buildup plate of EPID. The exit dose at extend SDD plane with patient CT data in the beam can be calculated as the predicted PDI in the treatment planning system (TPS). After that, the PDI was converted to the fluence at SDD of 150cm using the inverse square law coded in MATLAB. Five head-and-neck and prostate IMRT patient plans contain 32 fields were investigated to evaluate the feasibility of this new method. The measured EPID image was compared with PDI using the gamma analysis. Results: The average results for cumulative dose comparison were 81.9% and 91.6% for 3%, 3mm and 4%, 4mm gamma criteria, respectively. Results indicate that the patient transit dosimetry predicted algorithm compares well with EPID measured PD doses for test situations. Conclusion: Our new method can be used as an easy and feasible tool for online EPID-based in-vivo dose delivery verification for IMRT treatments. It can be implemented for fast detecting those obvious treatment delivery errors for individual field and patient quality assurance.

60 APPLIED LIFE SCIENCES↗

Multistep inversion workflow for 3D long-offset damped elastic waves in the Fourier domain

We present a new workflow for imaging damped three-dimensional elastic wavefields in the Fourier domain. The workflow employs a multiscale imaging approach, in which offset lengths are laddered, where frequency content and damping of the data are changed cyclically. Thus, the inversion process is launched using short-offset and low-frequency data to recover the long spatial wavelength of the image at a shallow depth. Increasing frequency and offset length leads to the recovery of the fine-scale features of the model at greater depths. For the fixed offset, we employ (in the imaging process) a few discrete frequencies with a set of Laplace damping parameters. The forward problem is solved with a finite-difference frequency-domain method based on a massively parallel iterative solver. The inversion code is based upon the solution of a least squares optimisation problem and is solved using a nonlinear gradient method. It is fully parallelised for distributed memory computational platforms. Our full-waveform inversion workflow is applied to the 3D Marmousi-2 and SEG/EAGE Salt models with long-offset data. The maximum inverted frequencies are 6 Hz for the Marmousi model and 2 Hz for the SEG/EAGE Salt model. The detailed structures are imaged successfully up to the depth approximately equal to one-third of the maximum offset length at a resolution consistent with the inverted frequencies.

3D full waveform inversion↗