Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 289 records · Page 16

Analysis and correction of errors in nanoscale particle tracking using the Single-pixel interior filling function (SPIFF) algorithm

Particle tracking, which is an essential tool in many fields of scientific research, uses algorithms that retrieve the centroid of tracked particles with sub-pixel accuracy. However, images in which the particles occupy a small number of pixels on the detector, are in close proximity to other particles or suffer from background noise, show a systematic error in which the particle sub-pixel positions are biased towards the center of the pixel. This “pixel locking” effect greatly reduces particle tracking accuracy. In this report, we demonstrate the severity of these errors by tracking experimental (and simulated) imaging data of optically trapped silver nanoparticles and single fluorescent proteins. We show that errors in interparticle separation, angle and mean square displacement are significantly reduced by applying the corrective Single-Pixel Interior Filling Function (SPIFF) algorithm. Our work demonstrates the potential ubiquity of such errors and the general applicability of SPIFF correction to many experimental fields.

Science & Technology - Other Topics↗

Data assimilation for combustion ignition delay time simulation using schlieren image velocimetry

This study sought to improve the accuracy of simulating spray penetration and combustion ignition delay by means of data assimilation (DA). The simulations were conducted using the Reynolds-averaged Navier–Stokes (RANS) equations and assimilating the schlieren image data. In DA, an ensemble square root filter (EnSRF) was used to build the statistical model, making the simulation results more accurate without any change in the governing equations. Recognizing that the spray-cone injection angle has a large effect on penetration, we created ensemble members with different injection angles. And we applied the two-component velocity distribution calculated via SIV and updated both velocity and temperature by using a DA statistical model derived from RANS ensemble simulations. The ignition delay time is generally known to vary even under the same experimental conditions because it is influenced by many factors. In this study, we attempted the transient DA-assisted RANS simulation to predict the ignition delay time even when the temporal resolution and accuracy of the observation data ware insufficient. Our trials offer an example of how a combination of techniques can be effectively used to assimilate experimental data obtained under restricted conditions.

Combustion simulation↗

Particle image velocimetry (PIV) uncertainty quantification using moment of correlation (MC) plane

Here, we present a new uncertainty estimation method for particle image velocimetry (PIV), that uses the correlation plane as a model for the probability density function (PDF) of displacements and calculates the second order moment of the correlation (MC). The cross-correlation between particle image patterns is the summation of all particle matches convolved with the apparent particle image diameter. MC uses this property to estimate the PIV uncertainty from the shape of the cross-correlation plane. In this new approach, the generalized cross-correlation (GCC) plane corresponding to a PIV measurement is obtained by removing the particle image diameter contribution. The GCC primary peak represents a discretization of the displacement PDF, from which the standard uncertainty is obtained by convolving the GCC plane with a Gaussian function. Then a Gaussian least-squares-fit is applied to the peak region, accounting for the stretching and rotation of the peak, due to the local velocity gradients and the effect of the convolved Gaussian. The MC method was tested with simulated image sets and the predicted uncertainties show good sensitivity to the error sources and agreement with the expected RMS error. Subsequently, the method was demonstrated in three PIV challenge cases and two experimental datasets and was compared with the published image matching (IM) and correlation statistics (CS) techniques. Results show that the MC method has a better response to spatial variation in RMS error and the predicted uncertainty is in good agreement with the expected standard uncertainty. The uncertainty prediction was also explored as a function of PIV interrogation window size. Overall, the MC method performance establishes itself as a valid uncertainty estimation tool for planar PIV.

47 OTHER INSTRUMENTATION↗

Dark Energy Studies with LSST using the Photon Simulator (Final Report)

This grant funded the development and dissemination of the Photon Simulator (PhoSim) for the purpose of studying dark energy at high precision with the upcoming Large Synoptic Survey Telescope (LSST) astronomical survey. The work was in collaboration with the Dark Energy Science Collaboration (DESC). Several detailed physics improvements were made in the optics, atmosphere and sensor, a number of validation studies were performed, and a significant number of usability features were implemented. Future work in DESC will use PhoSim as the image simulation tool for data challenges used by the analysis groups.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling X-ray emission around galaxies

Extended X-ray emission can be studied by spatial surface brightness measurements or by spectral analysis, but the two methods can disagree at low intensity levels. Here we present an improved method for spatial analysis that can be extended to include spectral information simultaneously. We construct a model for the entire image in a given energy band and generate a likelihood function to compare the model to the data. A critical goal is disentangling vignetted and unvignetted backgrounds through their different spatial distributions. Employing either maximum likelihood or Markov Chain Monte Carlo, we can derive probability distributions for the source and background parameters together, or we can fit and subtract the background, leaving the description of the source non-parametric. We calibrate this method against a variety of simulated images, and apply it to Chandra observations of the hot gaseous halo around the elliptical galaxy NGC 720. We follow the emission below a tenth of the background and infer a hot gas mass within 35 kpc of 4-5 × 10{sup 9} M {sub ☉}, with some indication that the profile continues to at least 50 kpc and that it steepens. We derive stronger constraints on the surface brightness profile than previous studies that employed the spectral method, and we show that the density profiles inferred from these studies are in conflict with the observed surface brightness profile. Contrary to a previous claim, we find that the X-ray halo does not contain the full complement of missing baryons within the virial radius.

79 ASTRONOMY AND ASTROPHYSICS↗

Muon-Computed Tomography using POCA Trajectory for Imaging Spent Nuclear Fuel in Dry Storage Casks

Cosmic ray muon-computed tomography (μCT) is a new imaging modality with unique characteristics that can be particularly important for applications in nuclear proliferation detection and international treaty verification. Using cosmic ray muons for nuclear security presents several potential advantages. Among others, muons are generated naturally in the atmosphere, can penetrate high-density materials, and are freely available. No radiological sources are required and consequently there is no associated radiological dose. Recently, the feasibility of using muons for imaging spent nuclear fuel stored in shielded casks has been explored and has been proved beneficial. However, challenges in μCT imaging include low muon flux of ~10,000 muons/m2/min, the effects of multiple Coulomb scattering (MCS) blurring the image, and inefficiency in being able to use all recorded muons for imaging. In this paper, we argue that the use of muon tracing should produce tomographic muon images with improved quality – or more quickly for the same image quality – compared to the case where conventional methods are used. In our paper, we report on the development and assessment of a novel muon tracing method for μCT. The proposed method back projects the muon’s scattering angle into each pixel crossed by its PoCA trajectory then forward projects the variance of the scattering angle in each pixel to detector bins along the muon’s incident horizontal direction. Two scenarios were simulated to assess the expected detection capability of this proposed method. GEANT4 was used to model the main characteristics of 1-60 GeV muons through matter. The simulated images showed an expected improvement in resolution and a reduced the reliance on the muon momentum information compared to a more conventional muon tomography method.

Liu, Zhengzhi↗

Numerical simulation of Li-pellet injection experiments for ELM-pacing in EAST

The injection of light impurity granules into fusion research discharges has shown a great potential for ELM pacing in future nuclear devices. The present work employed a model considering both pellet ablation and transport to systematically study for the first time the evolution of both impurity Li species and target plasma in the edge plasma during single and repeated Li pellet injections, with all initial information extracted from the EAST experiments. The evolutionary distributions of atomic and ionized Li species were presented and show good consistency with experimental CCD camera images. Simulation studies were carried out for single and repeated pellet injections at the outer mid-plane, and statistical uncertainty of ELM pacing rates triggering ELMs is discussed with regard to related injection conditions. It is found that the pressure locally caused by an ELM pacing pellet (with a diameter of 0.6 ~ 0.7 mm) was 1.8–2.4 times its unperturbed magnitude in the EAST shots of interest. Accumulation from repeated pellet injections results in the rise of pedestal height at a high injection frequency and this effect becomes stronger with increasing pellet size. Pressure changes are also evident in a scenario of multiple small pellets per stroke, which could lead to a nonnegligible, probabilistic ELM triggering rate observed experimentally. The effects of other possible factors are also discussed. With the assumption that the ELM trigger rate of a Li pellet of 0.7 mm diameter is 100%, the predicted trigger rates are in good agreement with those observed in experiments. This work can help understand the statistically probabilistic behavior of ELM triggering efficiency observed in previous studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Arcminute fluctuations in the microwave background from clusters of galaxies

A method for computing arcmin microwave fluctuations produced by Compton scattering of the cosmic background photons by hot electrons in clusters of galaxies is described. Microwave images of the sky for a range of Omega and primordial fluctuation spectral index n are generated which are then 'observed' to determine Delta T/T in precisely the same manner as actual observations to determine if the cluster-induced fluctuations are consistent with the measured upper limit. The geometry used by Uson and Wilkinson (1984) in the NRAO experiment and Readhead et al. (1989) in the OVRO experiment are applied to the simulated images. The 95 percent confidence lower limit for Omega is found to be about 1/10 for n = -1 (which approximates the CDM mass spectrum for clusters), while for n = 0 it is 1/7; for n = +1 the limit is 1/5 if the gas density profile extends to five core radii.

Markevitch, M.↗

Finding curvilinear features in speckled images

A method for finding curves in digital images with speckle noise is described. The solution method differs from standard linear convolutions followed by thresholds in that it explicitly allows curvature in the features. Maximum a posteriori (MAP) estimation is used, together with statistical models for the speckle noise and for the curve-generation process, to find the most probable estimate of the feature, given the image data. The estimation process is first described in general terms. Then, incorporation of the specific neighborhood system and a multiplicative noise model for speckle allows derivation of the solution, using dynamic programming, of the estimation problem. The detection of curvilinear features is considered separately. The detection results allow the determination of the minimal size of detectable feature. Finally, the estimation of linear features, followed by a detection step, is shown for computer-simulated images and for a SAR image of sea ice.

Samadani, Ramin↗

Inference from the small scales of cosmic shear with current and future Dark Energy Survey data

Cosmic shear is sensitive to fluctuations in the cosmological matter density field, including on small physical scales, where matter clustering is affected by baryonic physics in galaxies and galaxy clusters, such as star formation, supernovae feedback and AGN feedback. While muddying any cosmological information that is contained in small scale cosmic shear measurements, this does mean that cosmic shear has the potential to constrain baryonic physics and galaxy formation. We perform an analysis of the Dark Energy Survey (DES) Science Verification (SV) cosmic shear measurements, now extended to smaller scales, and using the Mead et al. 2015 halo model to account for baryonic feedback. While the SV data has limited statistical power, we demonstrate using a simulated likelihood analysis that the final DES data will have the statistical power to differentiate among baryonic feedback scenarios. We also explore some of the difficulties in interpreting the small scales in cosmic shear measurements, presenting estimates of the size of several other systematic effects that make inference from small scales difficult, including uncertainty in the modelling of intrinsic alignment on nonlinear scales, `lensing bias', and shape measurement selection effects. For the latter two, we make use of novel image simulations. While future cosmic shear datasets have the statistical power to constrain baryonic feedback scenarios, there are several systematic effects that require improved treatments, in order to make robust conclusions about baryonic feedback.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DETECTING UNSPECIFIED STRUCTURE IN LOW-COUNT IMAGES

Unexpected structure in images of astronomical sources often presents itself upon visual inspection of the image, but such apparent structure may either correspond to true features in the source or be due to noise in the data. This paper presents a method for testing whether inferred structure in an image with Poisson noise represents a significant departure from a baseline (null) model of the image. To infer image structure, we conduct a Bayesian analysis of a full model that uses a multiscale component to allow flexible departures from the posited null model. As a test statistic, we use a tail probability of the posterior distribution under the full model. This choice of test statistic allows us to estimate a computationally efficient upper bound on a p-value that enables us to draw strong conclusions even when there are limited computational resources that can be devoted to simulations under the null model. We demonstrate the statistical performance of our method on simulated images. Applying our method to an X-ray image of the quasar 0730+257, we find significant evidence against the null model of a single point source and uniform background, lending support to the claim of an X-ray jet.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗

Three-dimensional super line-localization in low signal-to-noise microscope images via prior-apprised unsupervised learning (PAUL)

Biological processes such as processive enzyme turnover and intracellular cargo tracking involve the dynamic motion of a small "article" along a curvilinear biopolymer track. To understand these processes that occur across multiple length and time scales, one must acquire both the trajectory of the particle and the position of the track along which it moves, possibly by combining high-resolution single-particle tracking with conventional microscopy. Yet, usually there is a significant resolution mismatch between these modalities: while the tracked particle is localized with a precision of 10 nm, the image of the surroundings is limited by optical difraction, with 200 nm lateral and 500 nm axial resolutions. Compared to the particle's trajectory, the surrounding curvilinear structure appears as a blurred and noisy image. This disparity in the spatial resolutions of the particle trajectory and the surrounding curvilinear structure image makes data reconstruction, as well as interpretation, particularly challenging. Analysis is further complicated when the curvilinear structures are oriented arbitrarily in 3D space. Here, we present a prior-apprised unsupervised learning (PAUL) approach to extract information from 3D images where the underlying features resemble a curved line such as a filament or microtubule. This three-stage framework starts with a Hessian-based feature enhancement, which is followed by feature registration, where local line segments are detected on repetitively sampled subimage tiles. In the final stage, statistical learning, segments are clustered based on their geometric relationships. Principal curves are then approximated from each segment group via statistical tools including principal component analysis, bootstrap and kernel transformation. This procedure is characterized on simulated images, where sub-voxel medium deviations from true curves have been achieved. The 3D PAUL approach has also been implemented for successful line localization in experimental 3D images of gold nanowires obtained using a multifocal microscope. Lastly, this work not only bridges the resolution gap between two microscopy modalities, but also allows us to conduct 3D super line-localization imaging experiments, without using super-resolution techniques.

42 ENGINEERING↗

Effects of overlapping sources on cosmic shear estimation: Statistical sensitivity and pixel-noise bias

The next generation of dark-energy imaging surveys — so called “Stage-IV” surveys, such as that of the Rubin Observatory Legacy Survey of Space and Time (LSST) — will cross a threshold in the number density of detected sources on the sky that requires qualitatively different image analysis and measurement techniques compared to the current generation of Stage-III surveys. In Stage-IV surveys, a significant amount of the cosmologically useful information is due to sources whose images overlap with those of other sources on the sky. Here, we focus on the weak gravitational lensing probe, for which we expect the largest impact since the cosmic shear signal is primarily encoded in the estimated shapes of observed galaxies and thus directly impacted by overlaps. We introduce a framework based on the Fisher formalism to analyze the effect of the overlapping sources (“blending”) on the estimation of cosmic shear. This method gives concrete predictions for the minimum loss of information due to noise and blending for any choice of “deblending” scheme and shape-measurement algorithm. Our studies account for undetected sources but do not address their full effects and biases they may introduce. We use simulated images and predict this impact of blending for three surveys: the Dark Energy Survey (DES), the Hyper-Suprime Cam Subaru Strategic Program (HSC-SSP), and the Rubin LSST. Our methodology successfully estimates the statistical sensitivity to weak lensing for DES and HSC early results. For LSST, we present the expected loss in statistical sensitivity for the ten-year survey due to blending. We find that for approximately 62% of galaxies that are likely to be detected in full-depth LSST images, at least 1% of the flux in their pixels is from overlapping sources. We also find that the statistical correlations between measures of overlapping galaxies and, to a much lesser extent (0.2%) the higher shot noise level due to their presence, decrease the effective number density of galaxies, N eff , by ~ 18%. We calculate an upper limit on N eff of 39.4 galaxies per arcmin 2 in r band. We study the impact of stars on as a function of stellar density and illustrate the diminishing returns of extending the survey into lower Galactic latitudes. We extend the simulation-based Fisher formalism to predict the expected increase in pixel-noise bias due to blending for maximum-likelihood (ML) shape estimators. We find that noise bias depends sensitively on the particular shape estimator and measure of ensemble-average shape that is used, and properties of the galaxy that include redshift-dependent quantities such as size and luminosity. The source code for these studies is available online.[The documented software developed for the catalog-level studies are available in the open-source LSST DESC github repository https://github.com/LSSTDESC/WeakLensingDeblending. The software for analyzing one or two galaxies with user-defined parameters is in the open-source github repository https://github.com/ismael-mendoza/ShapeMeasurementFisherFormalism.]

79 ASTRONOMY AND ASTROPHYSICS↗

Structure, magnetic and cytotoxic behaviour of solvothermally grown Fe3O4@Au core-shell nanoparticles

Highlights: • Synthesis of Fe{sub 3}O{sub 4}@Au nanoparticles was made by solvothermal method. • The magnetic nanoparticles was analyzed by XRD, FT-IR, magnetization curves and aberration corrected (Cs) STEM techniques. • Core-shell Fe{sub 3}O{sub 4}@Au nanoparticles showed low apoptotic response and therefore lower toxicity. • Fe{sub 3}O{sub 4}@Au nanoparticles with high saturation magnetization showed good biocompatibility. - Abstract: Magnetic nanoparticles (MNPs) have demonstrated their potential in a wide variety of biomedical applications including in tumour hyperthermia. However, highly reactive nature and aggregation affinity of these nanostructures are the principal limitations for such applications. To overcome these limitations, those MNPs should be covered with an inert shell in order to protect the magnetic core against chemical alterations. Considering the noble, chemically inert nature of gold, which also has good biocompatibility, we fabricated chemically stable gold-coated Fe{sub 3}O{sub 4} (Fe{sub 3}O{sub 4}@Au) nanoparticles of 20–50 nm size range with considerable saturation magnetization, though a simple one-pot solvothermal method. The morphology, elemental distribution and the structure of the core-shell nanoparticles were characterized using aberration-corrected (Cs) Scanning transmission electron microscopy (STEM), and compared with their simulated images. Cytotoxicity assay in MDCK cell line was performed to evaluate the biocompatibility of the nanoparticles. The results showed higher cell viability for the Fe{sub 3}O{sub 4}@Au nanoparticles, suggesting their good biocompatibility and potentials for biomedical applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Randomization of particle motions and the observed morphology of cometary heads

A great diversity is known to exist in the coma morphology of comets. In particular, some comets show much structural detail in their heads (such as jets, halos, fans, plumes, streamers), while others have completely structureless comas. Obvious questions arise as to why is this so and what does the presence or absence of features tell us about the emission processes on cometary nuclei. In an effort to investigate these problems, a computer code that generates synthetic images of dust comets was modified by introducing random perturbations into motions of ejected particles. It is noted that the introduction of perturbations has made the computer generated images simulate the appearance of comets quite faithfully and that by increasing the perturbations beyond a certain limit it has been possible to erase the coma morphology diagnostic of the details of the ejection process. It is proposed that the degree of collimation of an ejecta flow from discrete active sources on the nucleus surface and possible emissions of dust coma. Molecules of comet gases that radiate in the spectral region employed (and whose velocity distribution is much more chaotic than that of dust particles) and limited atmospheric seeing likewise contribute to blurring structural detail in ground-based imaging observation of comets. The absence of discrete features in the coma does no means imply the absence of localized sources of activity on the nucleus.

Sekanina, Zdenek↗

Two-level structural sparsity regularization for identifying lattices and defects in noisy images

Here, this paper presents a regularized regression model with a two-level structural sparsity penalty applied to locate individual atoms in a noisy scanning transmission electron microscopy image (STEM). In crystals, the locations of atoms is symmetric, condensed into a few lattice groups. Therefore, by identifying the underlying lattice in a given image, individual atoms can be accurately located. We propose to formulate the identification of the lattice groups as a sparse group selection problem. Furthermore, real atomic scale images contain defects and vacancies, so atomic identification based solely on a lattice group may result in false positives and false negatives. To minimize error, model includes an individual sparsity regularization in addition to the group sparsity for a within-group selection, which results in a regression model with a two-level sparsity regularization. We propose a modification of the group orthogonal matching pursuit (gOMP) algorithm with a thresholding step to solve the atom finding problem. The convergence and statistical analyses of the proposed algorithm are presented. The proposed algorithm is also evaluated through numerical experiments with simulated images. The applicability of the algorithm on determination of atom structures and identification of imaging distortions and atomic defects was demonstrated using three real STEM images. In conclusion, we believe this is an important step toward automatic phase identification and assignment with the advent of genomic databases for materials.

36 MATERIALS SCIENCE↗

Landsat Data Continuity Mission Expected Instrument Performance

The Landsat Data Continuity Mission (LDCM) is scheduled for a December 2012 launch date. LDCM is being managed by an interagency partnership between NASA and the U.S. Geological Survey (USGS). In order to provide the necessary spectral coverage of the visible through shortwave-infrared (SWIR) and the thermal-infrared (TIR), the satellite will carry two sensors. The Operational Land Imager (OLI) will collect data for nine visible to shortwave spectral bands with a spatial resolution of 30 m (with a 15 m panchromatic band). The Thermal Infrared Sensor (TIRS) will collect coincident image data for two TIR bands with a spatial resolution of 100 m. The OLI is fully assembled and tested and has been shipped by it's manufacturer, Ball Aerospace and Technology Corporation, to the Orbital Sciences Corporation (Orbital) facility where it is being integrated onto the LDCM spacecraft. Pre-launch testing indicates that OLI will meet all performance specification with margin. TIRS is in development at the NASA Goddard Space F!ight Center (GSFC) and is in final testing before shipping to the Orbital facility in January, 2012. The presentation will describe the LDCM satellite instrument systems, present pre-launch performance data for OLI and TIRS, and present simulated images to highlight notable features and expected imaging performance.

Dabney, Philip W.↗