Search NASA⌕ Search

SEARCH · Search NASA

Results for “image simulation”

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 757 records · Page 42

Combined speckle- and propagation-based single shot two-dimensional phase retrieval method

Single-shot two-dimensional (2D) phase retrieval (PR) can recover the phase shift distribution within an object from a single 2D x-ray phase contrast image (XPCI). Two competing XPCI imaging modalities often used for single-shot 2D PR to recover material properties critical for predictive performance capabilities are: speckle-based (SP-XPCI) and propagation-based (PB-XPCI) XPCI imaging. However, PR from SP-XPCI and PB-XPCI images are, respectively, limited to reconstructing accurately slowly and rapidly varying features due to noise and differences in their contrast mechanisms. Herein, we consider a combined speckle- and propagation-based XPCI (SPB-XPCI) image by introducing a mask to generate a reference pattern and imaging in the near-to-holographic regime to induce intensity modulations in the image. We develop a single-shot 2D PR method for SPB-XPCI images of pure phase objects without imposing restrictions such as object support constraints. It is compared against PR methods inspired by those developed for SP-XPCI and PB-XPCI on simulated and experimental images of a thin glass shell before and during shockwave compression. Reconstructed phase maps show improvements in quantitative scores of root-mean-square error and structural similarity index measure using our proposed method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-F-J-196: A Prototype System for Portal Imaging for Intensity Modulated Neutron Therapy

Purpose: Fast neutron therapy is offered at the University of Washington Medical Center for treatment of selected cancers. The hardware and control systems of the UW Clinical Neutron Therapy System are undergoing upgrades to enable delivery of IMNT. To clinically implement IMNT, dose verification tools need to be developed. We propose a portal imaging system that relies on the creation of positron emitting isotopes ({sup 11}C and {sup 15}O) through (n, 2n) reactions with a PMMA plate placed below the patient. After field delivery, the plate is retrieved from the vault and imaged using a reader that detects the annihilation photons. The pattern of activity produced in the plate provides information to reconstruct the neutron fluence map that can be compared to fluence maps from Monte Carlo (MCNP) simulations to verify treatment delivery. We have previously performed Monte Carlo simulations of the portal imaging system (GATE simulations) and the beam line (MCNP simulations). In this work, initial measurements using a prototype system are presented. Methods: Custom electronics were developed for BGO detectors read out with photomultiplier tubes (previous generation PET detectors from a CTI ECAT 953 scanner). Two detectors were placed in coincidence, with a detector separation of 2 cm. Custom software was developed to create the crystal look up tables and perform a limited angle planar reconstruction with a stochastic normalization. To test the initial capabilities of the system, PMMA squares were irradiated with neutrons at a depth of 1.5 cm and read out using the prototype system. Doses ranging from 10–200 cGy were delivered. Results: Using the prototype system, dose differences in the therapeutic range could be determined. Conclusion: The prototype portal imaging system is capable of detecting neutron doses as low as 10–50 cGy and shows great promise as a patient QA tool for IMNT.

60 APPLIED LIFE SCIENCES↗

Comprehensive study of the interactions between the critical dimensionless numbers associated with multiphase flow in 3D porous media

Multiphase flow in porous media is of signifcant interest in many engineering applications, such as geologic carbon sequestration, enhanced oil recovery, and groundwater contamination and remediation. In order to advance the fundamental understanding of multiphase flow in complex three-dimensional (3D) porous media, the interactions between the critical dimensionless numbers, including the contact angle, viscosity ratio, and capillary (Ca) number, were explored using X-ray micro-computed tomography (micro-CT) scanning and lattice Boltzmann (LB) modeling. In this study, the 3D pore structure information was extracted from micro-CT images and then used as interior boundary conditions of flow modeling in a pore-scale LB simulator to simulate multiphase flow within the pore space. A Berea sandstone sample was scanned and then two-phase flow LB simulations were performed based on the micro-CT images. The LB-simulated water/CO 2 distributions agreed well with the micro-CT scanned images. Simulation findings showed that a decreasing contact angle causes a decrease in wetting-fluid relative permeability and an increase in non-wetting fluid relative permeability. A rising Ca number increases both wetting and non-wetting fluid relative permeabilities. An increasing viscosity ratio (the ratio of non-wetting fluid viscosity to wetting fluid viscosity) facilitates the increase of non-wetting fluid relative permeability and mitigates the reduction of wetting fluid relative permeability, when the contact angle decreases continuously. The primary novel finding of this study is that the viscosity ratio affects the rate of change of the relative permeability curves for both phases when the contact angle changes continuously. To date, we believe it is the first time that comprehensive interactions between these dimensionless numbers are demonstrated in a sandstone sample based on real 3D structures. We also investigated the role of the changes of flow direction and sample location on relative permeability curves. Simulation results showed that the change in non-wetting fluid relative permeability was larger when the flow direction was switched from vertical to horizontal, which indicated that there was stronger anisotropy in larger pore networks that were primarily occupied by the non-wetting fluid.

42 ENGINEERING↗

Land cover/use classification of Cairns, Queensland, Australia: A remote sensing study involving the conjunctive use of the airborne imaging spectrometer, the large format camera and the thematic mapper simulator

In an attempt to improve the land cover/use classification accuracy obtainable from remotely sensed multispectral imagery, Airborne Imaging Spectrometer-1 (AIS-1) images were analyzed in conjunction with Thematic Mapper Simulator (NS001) Large Format Camera color infrared photography and black and white aerial photography. Specific portions of the combined data set were registered and used for classification. Following this procedure, the resulting derived data was tested using an overall accuracy assessment method. Precise photogrammetric 2D-3D-2D geometric modeling techniques is not the basis for this study. Instead, the discussion exposes resultant spectral findings from the image-to-image registrations. Problems associated with the AIS-1 TMS integration are considered, and useful applications of the imagery combination are presented. More advanced methodologies for imagery integration are needed if multisystem data sets are to be utilized fully. Nevertheless, research, described herein, provides a formulation for future Earth Observation Station related multisensor studies.

Heric, Matthew↗

SU-F-J-116: Clinical Experience-Based Verification and Improvement of a 4DCT Program

Purpose: To demonstrate the role of continuous improvement fulfilled by the Medical Physicist in clinical 4DCT and CBCT scanning. Methods: Lung (SABR and Standard) patients’ 4D respiratory motion and image data were reviewed over a 3, 6 and 12 month period following commissioning testing. By identifying trends of clinically relevant parameters and respiratory motions, variables were tested with a programmable motion phantom and assessed. Patient traces were imported to a motion phantom and 4DCT and CBCT imaging were performed. Cos6 surrogate and sup-inf motion was also programmed into the phantom to simulate the long exhale of patients for image contrast tests. Results: Patient surrogate motion amplitudes were 9.9+5.2mm (3–35) at 18+6bpm (6–30). Expiration/Inspiration time ratios of 1.4+0.5second (0.6–2.9) showed image contrast effects evident in the AveCT and 3DCBCT images. Small differences were found for patients with multiple 4DCT data sets. Patient motion assessments were simulated and verified with the phantom within 2mm. Initial image reviews to check for reconstructed artefacts and data loss identified a small number of patients with irregularities in the automatic placement of inspiration and expiration points. Conclusion: The Physicist’s involvement in the continuous improvements of a clinically commissioned technique, processes and workflows continues beyond the commissioning stage of a project. Our experience with our clinical 4DCT program shows that Physics presence is required at the clinical 4DCT scan to assist with technical aspects of the scan and also for clinical image quality assessment prior to voluming. The results of this work enabled the sharing of information from the Medical Physics group with the Radiation Oncologists and Radiation Therapists. This results in an improved awareness of clinical patient respiration variables and how they may affect 4D simulation images and also may also affect the treatment verification images.

60 APPLIED LIFE SCIENCES↗

Daylight imagery of ocean surface waves for wave spectra

A simple second-order model is developed to describe the transformation between surface wave slope and radiance incident on an optical sensor. The model includes azimuth as well as elevation variations in the sky-radiance distribution and their effects on the imaging of wave slopes approaching the sensor azimuth from various directions. The model also includes upwelling radiance. Comparison between the exact transformation from wave slope to radiance and the second-order model suggests that they yield quantitatively similar results. The model is used in physical interpretations of how wave slopes are imaged and in indicating a way of optimizing the linearity and contrast of the transfer function relating wave slope to radiance for various sensor geometries, thereby minimizing sources of error. One-dimensional numerical simulations of the imaging of an analytically generated wave-slope profile demonstrate the utility of this technique in obtaining slope spectra.

Monaldo, F. M.↗

SU-F-J-62: Assessment of Dose Changes Due to Anterior-Posterior Patient Separation During Daily MVCT Scans

Purpose: To determine the effects of patient separation on absolute dose and dose distribution in patients undergoing pelvic radiotherapy on TomoTherapy. Methods: An Alderson RANDO phantom with 4cm of bolus was imaged on a CT simulator and the resulting scans were contoured as a whole pelvic case. Using TomoTherapy Planning Station, the plan was designed to give 45 Gy to 95% of the treatment volume in 25 fractions. TomoTherapy MVCT scans were performed on the RANDO phantom with 2cm and 4cm of bolus removed to simulate visible changes in a patient’s anatomy. The MVCT images were rigidly registered with planning CT images on TomoTherapy Planned Adaptive. The original fluence was recalculated on the MVCT images and changes in dose distribution due to patient separation were quantified by the changes in DVHs for the target volume and the organs at risk. Results: Patient separation difference equivalent to 2cm and 4cm in anterior-posterior direction resulted in an increase of the PTV D50 and maximum PTV dose of 5.6%, 6.2% for 2cm and 7.7%, 10.4% for 4cm, respectively. For the 2cm change, D50 and maximum doses to organs at risk increased by 6.5%, 7.1% in the bladder, 4.9%, 4.8% in the rectum, and 5.3%, 6.6% in the bowel. For the 4cm change, D50 and maximum doses increased by 10.7%, 12.2% in the bladder, 5.9%, 6.1% in the rectum, and 7.7%, 10.1% in the bowel. Conclusion: This research indicates that, without any changes to the structures, patient separation in the anterior-posterior direction can affect the dose distribution for the PTV and organs at risk. These results can assist physicians in determining if obtaining a new CT simulation set and replanning is necessary for pelvic patients on TomoTherapy.

60 APPLIED LIFE SCIENCES↗

An online x-ray based position validation system for prostate hypofractionated radiotherapy

Purpose: Accurate positioning of the target volume during treatment is paramount for stereotactic body radiation therapy (SBRT). In this work, the authors present the development of an in-house software tool to verify target position with an Elekta-Synergy linear accelerator using kV planar images acquired during treatment delivery. Methods: In-house software, SeedTracker, was developed in MATLAB to perform the following three functions: 1. predict intended seed positions in a planar view perpendicular to any gantry angle, simulating a portal imaging device, from the 3D seed co-ordinates derived from the treatment planning system; 2. autosegment seed positions in kV planar images; and 3. report the position shift based on the seed positions in the projection images. The performance of SeedTracker was verified using a CIRS humanoid phantom (CIRS, VA, USA) implanted with three Civco gold seed markers (Civco, IA, USA) in the prostate. The true positive rate of autosegmentation (TPR{sub seg}) and the accuracy of the software in alerting the user when the isocenter position was outside the tolerance (TPR{sub trig}) were studied. Two-dimensional and 3D static position offsets introduced to the humanoid phantom and 3D dynamic offsets introduced to a gel phantom containing gold seeds were used for evaluation of the system. Results: SeedTracker showed a TPR{sub seg} of 100% in the humanoid phantom for projection images acquired at all angles except in the ranges of 80°–100° and 260°–280° where seeds are obscured by anatomy. This resulted in a TPR{sub trig} of 88% over the entire treatment range for considered 3D static offsets introduced to the phantom. For 2D static offsets where the position offsets were only introduced in the anterior–posterior and lateral directions, the TPR{sub trig} of SeedTracker was limited by both seed detectability and positional offset. SeedTracker showed a false positive trigger in the projection angle range between 130°–170° and 310°–350° (a maximum of 24% of treatment time) due to limited information that can be derived from monoscopic images. The system accurately determined the dynamic trajectory of the isocenter position in the superior and inferior direction for the studied dynamic offset scenarios based on the seed position in monoscopic images. Conclusions: The developed software has been shown to accurately autosegment the seed positions in kV planar images except for two 20° arcs where seeds are obscured by anatomical structures. The isocenter trajectories determined by the system, based on the monoscopic images, provide useful information for monitoring the prostate position. The developed system has potential application for monitoring prostate position during treatment delivery in linear accelerator based SBRT.

61 RADIATION PROTECTION AND DOSIMETRY↗

Quantifying biases in stellar masses of JWST high- z quasar host galaxies caused by quasar subtraction

The James Webb Space Telescope (JWST) has enabled dozens of high-z quasar host galaxy detections. Many of these observations imply galaxies with black holes that are overmassive compared to their low-z counterparts. However, the bright quasar point source removal can cause significant biases in recovered host magnitudes and stellar mass measurements due to the degeneracy in host galaxy and quasar light. We develop a statistical method to disentangle the quasar host galaxy stellar mass measurements from observational biases during the point source removal assuming the PSF is modelled perfectly. We use the BlueTides simulation to generate mock images and perform point source removal on thousands of simulated high-z quasar host galaxies, constructing corrected host magnitude posteriors. We find that removing a bright quasar in JWST photometry tends to either correctly recover or modestly misestimate host magnitudes, with a maximum magnitude underestimate of 0.2 mag. With our corrected magnitude posteriors, we perform SED fitting on each quasar host galaxy and compare the stellar mass measurement before and after the correction. We find that stellar mass estimates are generally robust, or misestimated by $<$ 0.3 dex. We also find that the stellar masses of a subset of hosts (J0844−0132, J0911+0152, and J1146−0005) remain unconstrained, as key photometric bands provide only flux upper limits. Accounting for observational biases does not resolve the apparent mismatch between black hole and host galaxy growth at high-z, where some quasars appear to host overmassive black holes while others reside in relatively massive galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

A Plug-and-Play Priors Approach for Solving Nonlinear Imaging Inverse Problems

In the past two decades, nonlinear image reconstruction methods have led to substantial improvements in the capabilities of numerous imaging systems. Such methods are traditionally formulated as optimization problems that are solved iteratively by simultaneously enforcing data consistency and incorporating prior models. Recently, the Plug-and-Play Priors (PPP) framework suggested that by using more sophisticated denoisers, not necessarily corresponding to an optimization objective, it is possible to improve the quality of reconstructed images. Here in this letter, we show that the PPP approach is applicable beyond linear inverse problems. In particular, we develop the fast iterative shrinkage/thresholding algorithm variant of PPP for model-based nonlinear inverse scattering. The key advantage of the proposed formulation over the original ADMM-based one is that it does not need to perform an inversion on the forward model. We show that the proposed method produces high quality images using both simulated and experimentally measured data.

97 MATHEMATICS AND COMPUTING↗

Blood flow imaging by optimal matching of computational fluid dynamics to 4D-flow data

Three-dimensional, time-resolved blood flow measurement (4D-flow) is a powerful research and clinical tool, but improved resolution and scan times are needed. Therefore, this study aims to (1) present a postprocessing framework for optimization-driven simulation-based flow imaging, called 4D-flow High-resolution Imaging with a priori Knowledge Incorporating the Navier-Stokes equations and the discontinuous Galerkin method (4D-flow HIKING), (2) investigate the framework in synthetic tests, (3) perform phantom validation using laser particle imaging velocimetry, and (4) demonstrate the use of the framework in vivo. An optimizing computational fluid dynamics solver including adjoint-based optimization was developed to fit computational fluid dynamics solutions to 4D-flow data. Synthetic tests were performed in 2D, and phantom validation was performed with pulsatile flow. Reference velocity data were acquired using particle imaging velocimetry, and 4D-flow data were acquired at 1.5 T. In vivo testing was performed on intracranial arteries in a healthy volunteer at 7 T, with 2D flow as the reference. Results Synthetic tests showed low error (0.4%-0.7%). Phantom validation showed improved agreement with laser particle imaging velocimetry compared with input 4D-flow in the horizontal (mean -0.05 vs -1.11 cm/s, P < .001; SD 1.86 vs 4.26 cm/s, P < .001) and vertical directions (mean 0.05 vs -0.04 cm/s, P = .29; SD 1.36 vs 3.95 cm/s, P < .001). In vivo data show a reduction in flow rate error from 14% to 3.5%. Phantom and in vivo results from 4D-flow HIKING show promise for future applications with higher resolution, shorter scan times, and accurate quantification of physiological parameters.

4D-flow MRI↗

Aggregation of Size-Selected Oxide Clusters Deposited onto Au(111)

Kinetic Monte Carlo (kMC) simulations along with density functional theory (DFT) calculations were used to investigate the aggregation of size-selected Nb 3 O y (y = 5, 6, 7) clusters deposited onto the Au(111) surface. Recent STM experiments showed that the cluster binding sites and sizes of the cluster assemblies on the Nb 3 O y /Au(111) surfaces strongly depend on the stoichiometry of the clusters, i.e., the oxygen-to-niobium ratio. To better understand the origins of these differences, kMC simulations of the nucleation and growth of cluster assemblies were performed using energy barriers for diffusion and intercluster interactions estimated from DFT calculations of cluster binding and dimerization energies, respectively. Comparisons of the kMC simulations with STM images of the as-deposited Nb 3 O y /Au(111) surfaces at RT and after high temperature annealing were used to further optimize the energetics and gauge the importance of nearest neighbor interactions. The kMC simulations demonstrate that the assembly of Nb 3 O y clusters on Au(111) are largely controlled by the magnitude of the barriers for diffusion and interparticle-bond formation, while changes at higher temperatures are sensitive to the binding energies between nearest neighbors. Simulations for the Nb 3 O 5 and Nb 3 O 6 clusters, which exhibit smaller cluster assembly sizes in STM, required larger diffusion barriers as well as different barriers for interparticle binding, which reflected differences in DFT calculated dimerization energies. Furthermore, the results demonstrate the effectiveness of combined DFT and kMC calculations for understanding how the stoichiometry affects the aggregation of small oxide clusters on a metal surface.

36 MATERIALS SCIENCE↗

Interpretation of proton radiography experiments of hohlraums with three-dimensional simulations [Interpretation of proton radiography experiments of hohlraums with 3D simulations]

Proton radiography experiments of laser-irradiated hohlraums performed at the OMEGA laser facility are analyzed using three-dimensional (3D) hydrodynamic simulations coupled to a proton trajectography package. Experiments with three different laser irradiation patterns were performed, and each produced a distinct proton image. By comparing these results with synthetic proton images obtained by sending protons through plasma profiles in the hohlraum obtained from 3D radiation hydrodynamic simulations, it is found that the simulated images agree favorably with the experimental images when electric fields, due to the electron pressure gradients that arise from 3D structures occurring during plasma expansion, are included. Furthermore, these comparisons provide quantitative estimates of the electric field present inside the hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interpretation of proton radiography experiments of hohlraums with three-dimensional simulations

Proton radiography experiments of laser-irradiated hohlraums performed at the OMEGA laser facility are analyzed using three-dimensional (3D) hydrodynamic simulations coupled to a proton trajectography package. Experiments with three different laser irradiation patterns were performed, and each produced a distinct proton image. By comparing these results with synthetic proton images obtained by sending protons through plasma profiles in the hohlraum obtained from 3D radiation hydrodynamic simulations, it is found that the simulated images agree favorably with the experimental images when electric fields, due to the electron pressure gradients that arise from 3D structures occurring during plasma expansion, are included. These comparisons provide quantitative estimates of the electric field present inside the hohlraums

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of electronic registration of multispectral remote sensing images to 0.1 pixel accuracy

Band-to-band coregistration of multispectral remote sensing images can be achieved by electronic signal processing techniques rather than by costly and difficult mechanical alignment. This paper describes the results of a study of the end-to-end performance of electronic registration. The software simulation includes steps which model the performance of the geometric calibration process, the instrument image quality, detector performance and the effects of achieving coregistration through image resampling. The image resampling step emulates the Pipelined Resampling Processor, a real-time image resampler. The study demonstrates that the electronic alignment technique produces multispectral images which are superior to those produced by an imager whose pixel geometry is accurate to 0.1 pixel rms. The implications of this approach for future earth observation programs are discussed.

Reitsema, H. J.↗

Impact of Real-Time Image Gating on Spot Scanning Proton Therapy for Lung Tumors: A Simulation Study

Purpose: To investigate the effectiveness of real-time-image gated proton beam therapy for lung tumors and to establish a suitable size for the gating window (GW). Methods and Materials: A proton beam gated by a fiducial marker entering a preassigned GW (as monitored by 2 fluoroscopy units) was used with 7 lung cancer patients. Seven treatment plans were generated: real-time-image gated proton beam therapy with GW sizes of ±1, 2, 3, 4, 5, and 8 mm and free-breathing proton therapy. The prescribed dose was 70 Gy (relative biological effectiveness)/10 fractions to 99% of the target. Each of the 3-dimensional marker positions in the time series was associated with the appropriate 4-dimensional computed tomography phase. The 4-dimensional dose calculations were performed. The dose distribution in each respiratory phase was deformed into the end-exhale computed tomography image. The D99 and D5 to D95 of the clinical target volume scaled by the prescribed dose with criteria of D99 >95% and D5 to D95 <5%, V20 for the normal lung, and treatment times were evaluated. Results: Gating windows ≤ ±2 mm fulfilled the CTV criteria for all patients (whereas the criteria were not always met for GWs ≥ ±3 mm) and gave an average reduction in V20 of more than 17.2% relative to free-breathing proton therapy (whereas GWs ≥ ±4 mm resulted in similar or increased V20). The average (maximum) irradiation times were 384 seconds (818 seconds) for the ±1-mm GW, but less than 226 seconds (292 seconds) for the ±2-mm GW. The maximum increased considerably at ±1-mm GW. Conclusion: Real-time-image gated proton beam therapy with a GW of ±2 mm was demonstrated to be suitable, providing good dose distribution without greatly extending treatment time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Onboard utilization of ground control points for image correction. Volume 3: Ground control point simulation software design

The software developed to simulate the ground control point navigation system is described. The Ground Control Point Simulation Program (GCPSIM) is designed as an analysis tool to predict the performance of the navigation system. The system consists of two star trackers, a global positioning system receiver, a gyro package, and a landmark tracker.

Source record↗