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.
A new method for point-spread function correction using the ellipticity of re-smeared artificial images in weak gravitational lensing shear analysis
Highly accurate weak lensing analysis is urgently required for planned cosmic shear observations. For this purpose we have eliminated various systematic noises in the measurement. The point-spread function (PSF) effect is one of them. A perturbative approach for correcting the PSF effect on the observed image ellipticities has been previously employed. Here we propose a new non-perturbative approach for PSF correction that avoids the systematic error associated with the perturbative approach. The new method uses an artificial image for measuring shear which has the same ellipticity as the lensed image. This is done by re-smearing the observed galaxy images and observed star images (PSF) with an additional smearing function to obtain the original lensed galaxy images. We tested the new method with simple simulated objects that have Gaussian or Sérsic profiles smeared by a Gaussian PSF with sufficiently large size to neglect pixelization. Under the condition of no pixel noise, it is confirmed that the new method has no systematic error even if the PSF is large and has a high ellipticity.
WE-D-9A-02: Automated Landmark-Guided CT to Cone-Beam CT Deformable Image Registration
Purpose: The anatomical changes that occur between the simulation CT and daily cone-beam CT (CBCT) are investigated using an automated landmark-guided deformable image registration (LDIR) algorithm with simultaneous intensity correction. LDIR was designed to be accurate in the presence of tissue intensity mismatch and heavy noise contamination. Method: An auto-landmark generation algorithm was used in conjunction with a local small volume (LSV) gradient matching search engine to map corresponding landmarks between the CBCT and planning CT. The LSVs offsets were used to perform an initial deformation, generate landmarks, and correct local intensity mismatch. The landmarks act as stabilizing controlpoints in the Demons objective function. The accuracy of the LDIR algorithm was evaluated on one synthetic case with ground truth and data of ten head and neck cancer patients. The deformation vector field (DVF) accuracy was accessed using a synthetic case. The Root mean square error of the 3D canny edge (RMSECE), mutual information (MI), and feature similarity index metric (FSIM) were used to access the accuracy of LDIR on the patient data. The quality of the corresponding deformed contours was verified by an attending physician. Results: The resulting 90 percentile DVF error for the synthetic case was within 5.63mm for the original demons algorithm, 2.84mm for intensity correction alone, 2.45mm using controlpoints without intensity correction, and 1.48 mm for the LDIR algorithm. For the five patients the mean RMSECE of the original CT, Demons deformed CT, intensity corrected Demons CT, control-point stabilized deformed CT, and LDIR CT was 0.24, 0.26, 0.20, 0.20, and 0.16 respectively. Conclusion: LDIR is accurate in the presence of multimodal intensity mismatch and CBCT noise contamination. Since LDIR is GPU based it can be implemented with minimal additional strain on clinical resources. This project has been supported by a CPRIT individual investigator award RP11032.
Imaging solar flares in hard X-rays using Fourier telescopes
The sun emits hard X-rays (above 10 keV) during solar flares. Imaging hard X-ray sources on the sun with spatial resolutions on the order of 1-5 arcsec and integration times of 1 sec will provide greater insight into the energy release processes during a solar flare. In these events, tremendous amounts of energy stored in the solar magnetic field are rapidly released leading to emission across the electromagnetic spectrum. Two Fourier telescope designs, a spatial modulation collimator and a rotating modulation collimator, were developed to image the full sun in hard X-rays (10-100 keV) in an end-to-end simulation. Emission profiles were derived for two hard X-ray solar flare models taken from the current solar theoretical literature and used as brightness distributions for the telescope simulations. Both our telescope models, tailored to image solar sources, were found to perform equally well, thus offering the designer significant flexibility in developing systems for space-based platforms. Given sufficient sensitive areas, Fourier telescopes are promising concepts for imaging solar hard X-rays.
SU-E-T-528: Robustness Evaluation for Fiducial-Based Accelerated Partial Breast Proton Therapy
Purpose: To investigate the robustness of the proton treatment plans in the presence of rotational setup error when patient is aligned with implanted fiducials. Methods: Five Stage I invasive breast cancer patients treated with the APBP protocol (PCG BRE007-12) were studied. The rotational setup errors were simulated by rotating the original CT images around the body center clockwise and counterclockwise 5 degrees (5CW and 5CCW). Manual translational registration was then performed to match the implanted fiducials on the rotated images to the original dataset. Patient contours were copied to the newly created CT set. The original treatment plan was applied to the new CT dataset with the beam isocenter placed at the geometrical center of PTV. The dose distribution was recalculated for dosimetric parameters comparison. Results: CTV and PTV (D95 and V95) coverages were not significantly different between the two simulated plans (5CW and 5CCW) and the original plan. PTV D95 and CTV D95 absolute difference among the three plans were relatively small, with maximum changes of 0.28 CGE and 0.15 CGE, respectively. PTV V95 and CTV V95 absolute differences were 0.79% and 0.48%. The dosage to the thyroid, heart, contralateral breast and lung remained zero for all three plans. The Dmax and Dmean to the volume of ipsilateral breast excluding CTV were compared, with maximum difference values of 1.02 CGE for Dmax and 3.56 CGE for Dmean. Ipsilateral lung Dmean maintained no significant changes through the three plan comparison, with the largest value 0.32 CGE. Ipsilateral lung Dmax was the most sensitive parameter to this simulation study, with a maximum difference at 20.2 CGE. Conclusion: Our study suggests that fiducial-based Accelerated Partial Breast Proton Therapy is robust with respect to +/− 5 degree patient setup rotational errors, as long as the internal fiducial markers are used for patient alignment.
SU-E-J-205: Monte Carlo Modeling of Ultrasound Probes for Real-Time Ultrasound Image-Guided Radiotherapy
Purpose: To quantify the effect of ultrasound (US) probe beam attenuation for radiation therapy delivered under real-time US image guidance by means of Monte Carlo (MC) simulations. Methods: MC models of two Philips US probes, an X6-1 matrix-array transducer and a C5-2 curved-array transducer, were built based on their CT images in the EGSnrc BEAMnrc and DOSXYZnrc codes. Due to the metal parts, the probes were scanned in a Tomotherapy machine with a 3.5 MV beam. Mass densities in the probes were assigned based on an electron density calibration phantom consisting of cylinders with mass densities between 0.2–8.0 g/cm{sup 3}. Beam attenuation due to the probes was measured in a solid water phantom for a 6 MV and 15 MV 15x15 cm{sup 2} beam delivered on a Varian Trilogy linear accelerator. The dose was measured with the PTW-729 ionization chamber array at two depths and compared to MC simulations. The extreme case beam attenuation expected in robotic US image guided radiotherapy for probes in upright position was quantified by means of MC simulations. Results: The 3.5 MV CT number to mass density calibration curve was found to be linear with R{sup 2} > 0.99. The maximum mass densities were 4.6 and 4.2 g/cm{sup 3} in the C5-2 and X6-1 probe, respectively. Gamma analysis of the simulated and measured doses revealed that over 98% of measurement points passed the 3%/3mm criteria for both probes and measurement depths. The extreme attenuation for probes in upright position was found to be 25% and 31% for the C5-2 and X6-1 probe, respectively, for both 6 and 15 MV beams at 10 cm depth. Conclusion: MC models of two US probes used for real-time image guidance during radiotherapy have been built. As a Result, radiotherapy treatment planning with the imaging probes in place can now be performed. J Schlosser is an employee of SoniTrack Systems, Inc. D Hristov has financial interest in SoniTrack Systems, Inc.
Simulation tools for atom probe tomography: A path for diagnosis and treatment of image degradation
The ideal picture of a near-perfect 3D microscope often presented regarding Atom Probe Tomography faces several issues. These issues degrade the metrological performance of the instrument and find their roots in the phenomena acting at the atomic to the mesoscopic level in the vicinity of the surface of a field emitter. From the field evaporation process at the atomic scale, to the macroscopic scale of the instrument, the path to model the imaging process and to develop more accurate and reliable reconstruction algorithms is not a single lane road. This paper focused on the numerical methods used to understand, treat, and potentially heal imaging issues commonly affecting the data in atom probe experiments. A lot of room for improvement exists in solving accuracy problems observed in complex materials by means of purely electrostatic models describing the image formation in a classical approach. Looking at the sample at the atomic scale, the phenomena perturbing the imaging process are subtle. An examination of atomic scale modifications of the sample surface in the presence of a high surface electric field is therefore mandatory. Atomic scale molecular dynamic models integrating the influence of the high surface electric are being developed with this aim. It is also demonstrated that the complex behavior of atoms and molecules in high fields, and consequences on collected data, can be understood through the use of accurate ab-initio models modified to include the impact of the extreme surface electric field.
Pore-scale intermittent velocity structure underpinning anomalous transport through 3-D porous media
We study the nature of non-Fickian particle transport in 3-D porous media by simulating fluid flow in the intricate pore space of real rock. We solve the full Navier-Stokes equations at the same resolution as the 3-D micro-CT (computed tomography) image of the rock sample and simulate particle transport along the streamlines of the velocity field. We find that transport at the pore scale is markedly anomalous: longitudinal spreading is superdiffusive, while transverse spreading is subdiffusive. We demonstrate that this anomalous behavior originates from the intermittent structure of the velocity field at the pore scale, which in turn emanates from the interplay between velocity heterogeneity and velocity correlation. Finally, we propose a continuous time random walk model that honors this intermittent structure at the pore scale and captures the anomalous 3-D transport behavior at the macroscale.
Design of novel refractory equiatomic multi-principal elemental alloys based on Mo-Nb-Ti system for Gen IV reactor applications
Excellent irradiation damage resistance demonstrated by multi-principal elemental alloys (MPEAs) has sparked significant interest among researchers, prompting exploration into their vast compositional space, to validate their suitability for nuclear applications. A combined approach of thermodynamic and empirical parameters calculations alongside CALPHAD (CALculation of PHAse Diagrams) for phase formation predictions enable high-throughput material selection for sophisticated applications like nuclear, overcoming laborious and time-consuming experiments. Key thermodynamic and empirical parameters for eight novel equiatomic MPEAs, based on seven low thermal neutron cross section refractory elements, for predicting phase formation were calculated, and equilibrium and non-equilibrium simulations in CALPHAD were employed to comprehensively model the systems. Pseudo binary phase diagram simulations showed that Zr, V or equiatomic CrV additions to the base MoNbTi alloy (MoNbTi-Zr, MoNbTi-V and MoNbTi-CrV alloys) favor the formation of isomorphous body-centered cubic (BCC) phase at high temperatures, while Cr, Al, equiatomic ZrV, or equiatomic CrAl additions (MoNbTi-Cr, MoNbTi-Al, MoNbTi-ZrV or MoNbTi-CrAl alloys) limit the solubility of them. Equilibrium CALPHAD simulations at 750 oC were consistent with XRD results on MoNbTi, MoNbTiZr and MoNbTiCr alloys, and partially for others. Notably, elemental segregation observed in the backscattered electron (BSE) scanning electron microscopy (SEM) images of the alloys was accurately simulated through non-equilibrium Scheil solidification calculations in CALPHAD, further verified by experiments. The precipitation of TiCr2 Laves phase in Cr containing MoNbTiCr and MoNbTiCrAl was accurately predicted while discrepancies were noted in MoNbTiCrV. The equilibrium simulations also provided insights into phase compositions at specific temperatures offering a pathway for tailoring the desired microstructure and properties of these systems. Empirical parameters calculations successfully predicted random solid solution in the base MoNbTi alloy, and with an exception in MoNbTiV and MoNbTiAl, predicted intermetallic precipitation in the rest, especially, Laves phase precipitation in Cr containing alloys.
SU-C-207A-06: On-Line Beam Range Verification with Multiple Scanning Particle Beams: Initial Feasibility Study with Simulations
Purpose: To investigate the feasibility and requirement for intra-fraction on-line multiple scanning particle beam range verifications (BRVs) with in-situ PET imaging, which is beyond the current single-beam BRV with extra factors that will affect the BR measurement accuracy, such as beam diameter, separation between beams, and different image counts at different BRV positions. Methods: We simulated a 110-MeV proton beam with 5-mm diameter irradiating a uniform PMMA phantom by GATE simulation, which generated nuclear interaction-induced positrons. In this preliminary study, we simply duplicated these positrons and placed them next to the initial protons to approximately mimic the two spatially separated positron distributions produced by two beams parallel to each other but with different beam ranges. These positrons were then imaged by a PET (∼2-mm resolution, 10% sensitivity, 320×320×128 mm^3 FOV) with different acquisition times. We calculated the positron activity ranges (ARs) from reconstructed PET images and compared them with the corresponding ARs of original positron distributions. Results: Without further image data processing and correction, the preliminary study show the errors between the measured and original ARs varied from 0.2 mm to 2.3 mm as center-to-center separations and range differences were in the range of 8–12 mm and 2–8 mm respectively, indicating the accuracy of AR measurement strongly depends on the beam separations and range differences. In addition, it is feasible to achieve ≤ 1.0-mm accuracy for both beams with 1-min PET acquisition and 12 mm beam separation. Conclusion: This study shows that the overlap between the positron distributions from multiple scanning beams can significantly impact the accuracy of BRVs of distributed particle beams and need to be further addressed beyond the established method of single-beam BRV, but it also indicates the feasibility to achieve accurate on-line multi-beam BRV with further improved method.
Sponge-like Si-SiO{sub 2} nanocomposite—Morphology studies of spinodally decomposed silicon-rich oxide
Sponge-like Si nanostructures embedded in SiO{sub 2} were fabricated by spinodal decomposition of sputter-deposited silicon-rich oxide with a stoichiometry close to that of silicon monoxide. After thermal treatment a mean feature size of about 3 nm was found in the phase-separated structure. The structure of the Si-SiO{sub 2} nanocomposite was investigated by energy-filtered transmission electron microscopy (EFTEM), EFTEM tomography, and atom probe tomography, which revealed a percolated Si morphology. It was shown that the percolation of the Si network in 3D can also be proven on the basis of 2D EFTEM images by comparison with 3D kinetic Monte Carlo simulations.
SU-E-I-15: Comparison of Radiation Dose for Radiography and EOS in Adolescent Scoliosis Patients
Purpose: To estimate patient radiation dose for whole spine imaging using EOS, a new biplanar slot-scanning radiographic system and compare with standard scoliosis radiography. Methods: The EOS imaging system (EOS Imaging, Paris, France) consists of two orthogonal x-ray fan beams which simultaneously acquire frontal and lateral projection images of a standing patient. The patient entrance skin air kerma was measured for each projection image using manufacturer-recommended exposure parameters for spine imaging. Organ and effective doses were estimated using a commercially-available Monte Carlo simulation program (PCXMC, STUK, Radiation and Nuclear Safety Authority, Helsinki, Finland) for a 15 year old mathematical phantom model. These results were compared to organ and effective dose estimated for scoliosis radiography using computed radiography (CR) with standard exposure parameters obtained from a survey of pediatric radiographic projections. Results: The entrance skin air kerma for EOS was found to be 0.18 mGy and 0.33 mGy for posterior-anterior (PA) and lateral projections, respectively. This compares to 0.76 mGy and 1.4 mGy for CR, PA and lateral projections. Effective dose for EOS (PA and lateral projections combined) is 0.19 mSv compared to 0.51 mSv for CR. Conclusion: The EOS slot-scanning radiographic system allows for reduced patient radiation dose in scoliosis patients as compared to standard CR radiography.
Enhanced Monte Carlo Simulations for Electron Energy Loss Mitigation in Real-Space Nanoimaging of Thick Biological Samples and Microchips
High-resolution imaging using Transmission Electron Microscopy (TEM) is essential for applications such as grain boundary analysis, microchip defect characterization, and biological imaging. However, TEM images are often compromised by electron energy spread and other factors. In TEM mode, where the objective and projector lenses are positioned downstream of the sample, electron–sample interactions cause energy loss, which adversely impacts image quality and resolution. This study introduces a simulation tool to estimate the electron energy loss spectrum (EELS) as a function of sample thickness, covering electron beam energies from 300 keV to 3 MeV. Leveraging recent advances in MeV-TEM/STEM technology, which includes a state-of-the-art electron source with 2-picometer emittance, an energy spread of 3 × 10 -5 , and optimized beam characteristics, we aim to minimize energy spread. By integrating EELS capabilities into the BNL Monte Carlo (MC) simulation code for thicker samples, we evaluate electron beam parameters to mitigate energy spread resulting from electron–sample interactions. Based on our simulations, we propose an experimental procedure for quantitively distinguishing between elastic and inelastic scattering. The findings will guide the selection of optimal beam settings, thereby enhancing resolution for nanoimaging of thick biological samples and microchips.
Analytical investigation of the dynamics of tethered constellations in Earth orbit (phase 2)
Simulation of two short distance crawling maneuvers of the elevator both with and without environmental perturbations acting upon the system is discussed. These simulation runs were performed in order to provide results useful for the interpretation of the data from the tests, on the ground, of a scaled down engineering model of the elevator. In these simulation runs the elevator crawls along the tether in accordance to the developed mirror image motion control law (MIMCL). Results from the simulation of the 4 km long maneuver run were compared to those obtained by adopting the modified hyperbolic tangent control law (MHTCL). A preprocessor was developed for setting up the initial conditions of a tethered system with L platforms, M longitudinal dampers, and N lumped masses (platforms plus tether beads). A short test run of the 4-platform system with 3 longitudinal dampers and 10 lumper mass without any perturbation acting upon the system is illustrated. In support of the Tether Applications Working Group, SAO is preparing a catalog of tether simulations, has prepared a set of simulator test cases, obtained results from SKYHOOK, and solicited results from simulators at other institutions, and prepared a paper on a specific analytic solution.
Monte Carlo simulations of secondary electron emission due to ion beam milling
We present a Monte Carlo simulation study of secondary electron emission resulting from focused ion beam milling of a copper target. The basis of this study is a simulation code which simulates ion induced excitation and emission of secondary electrons, in addition to simulating focused ion beam sputtering and milling. This combination of features permits the simulation of the interaction between secondary electron emission, and the evolving target geometry as the ion beam sputters material. Previous ion induced SE Monte Carlo simulation methods have been restricted to predefined target geometries, while the dynamic target in the presented simulations makes this study relevant to image formation in ion microscopy, and chemically assisted ion beam etching, where the relationship between sputtering, and its effects on secondary electron emission, is important. We focus on a copper target, and validate our simulation against experimental data for a range of: noble gas ions, ion energies, ion/substrate angles and the energy distribution of the secondary electrons. We then provide a detailed account of the emission of secondary electrons resulting from ion beam milling; we quantify both the evolution of the yield as high aspect ratio valleys are milled, as well as the emission of electrons within these valleys that do not escape the target, but which are important to the secondary electron contribution to chemically assisted ion induced etching.
Impact disturbance of the Venus atmosphere
Experimental simulations of the atmosphere-surface interaction during high-velocity impact are presented. At Venus an atmospheric vortex, generated with impact, may interact with a local wind. Some surface features, observed on Magellan images may be related with the simulated effect.
Denoising of polychromatic CT images based on their own noise properties
Purpose: Because of high diagnostic accuracy and fast scan time, computed tomography (CT) has been widely used in various clinical applications. Since the CT scan introduces radiation exposure to patients, however, dose reduction has recently been recognized as an important issue in CT imaging. However, low-dose CT causes an increase of noise in the image and thereby deteriorates the accuracy of diagnosis. In this paper, the authors develop an efficient denoising algorithm for low-dose CT images obtained using a polychromatic x-ray source. The algorithm is based on two steps: (i) estimation of space variant noise statistics, which are uniquely determined according to the system geometry and scanned object, and (ii) subsequent novel conversion of the estimated noise to Gaussian noise so that an existing high performance Gaussian noise filtering algorithm can be directly applied to CT images with non-Gaussian noise. Methods: For efficient polychromatic CT image denoising, the authors first reconstruct an image with the iterative maximum-likelihood polychromatic algorithm for CT to alleviate the beam-hardening problem. We then estimate the space-variant noise variance distribution on the image domain. Since there are many high performance denoising algorithms available for the Gaussian noise, image denoising can become much more efficient if they can be used. Hence, the authors propose a novel conversion scheme to transform the estimated space-variant noise to near Gaussian noise. In the suggested scheme, the authors first convert the image so that its mean and variance can have a linear relationship, and then produce a Gaussian image via variance stabilizing transform. The authors then apply a block matching 4D algorithm that is optimized for noise reduction of the Gaussian image, and reconvert the result to obtain a final denoised image. To examine the performance of the proposed method, an XCAT phantom simulation and a physical phantom experiment were conducted. Results: Both simulation and experimental results show that, unlike the existing denoising algorithms, the proposed algorithm can effectively reduce the noise over the whole region of CT images while preventing degradation of image resolution. Conclusions: To effectively denoise polychromatic low-dose CT images, a novel denoising algorithm is proposed. Because this algorithm is based on the noise statistics of a reconstructed polychromatic CT image, the spatially varying noise on the image is effectively reduced so that the denoised image will have homogeneous quality over the image domain. Through a simulation and a real experiment, it is verified that the proposed algorithm can deliver considerably better performance compared to the existing denoising algorithms.
Characterization of Single-Sided Charge-Sharing CZT Strip Detectors for Gamma-Ray Astronomy
We report progress in the study of thick single-sided charge-sharing cadmium zinc telluride (CZT) strip detector modules designed to perform spectroscopy and 3-D imaging of gamma-rays. We report laboratory measurements including spectroscopy, efficiency and 3-D imaging capability of prototype detectors (15 15 7.5 cu mm) with 11x11 unit cells. We also report on Monte Carlo simulations (GEANT4 v7.1) to investigate the effect of multihits on detector performance in both spectroscopy and imaging. We compare simulation results with data obtained from laboratory measurements and discuss the implications for future strip detector designs. Keywords: CZT, strip detectors, gamma-ray