Simulating Image Coaddition with the Nancy Grace Roman Space Telescope. III. Software Improvements and New Linear Algebra Strategies
Explore the source record for details and available documents.
Kaili 开力 Cao 曹↗
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.
Constraining the masses of high-redshift clusters with weak lensing: Revised shape calibration testing for the impact of stronger shears and increased blending
Weak lensing measurements suffer from well-known shear estimation biases, which can be partially corrected for with the use of image simulations. Here we present an analysis of simulated images that mimic Hubble Space Telescope/Advance Camera for Surveys observations of high-redshift galaxy clusters, including cluster specific issues such as non-weak shear and increased blending. Our synthetic galaxies have been generated to have similar observed properties as the background-selected source samples studied in the real images. First, we used simulations with galaxies placed on a grid to determine a revised signal-to-noise-dependent (S/N KSB ) correction for multiplicative shear measurement bias, and to quantify the sensitivity of our KSB+ bias calibration to mismatches of galaxy or PSF properties between the real data and the simulations. Next, we studied the impact of increased blending and light contamination from cluster and foreground galaxies, finding it to be negligible for high-redshift (z > 0.7) clusters, whereas shear measurements can be affected at the ~1% level for lower redshift clusters given their brighter member galaxies. Finally, we studied the impact of fainter neighbours and selection bias using a set of simulated images that mimic the positions and magnitudes of galaxies in Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) data, thereby including realistic clustering. While the initial SExtractor object detection causes a multiplicative shear selection bias of –0.028 ± 0.002, this is reduced to –0.016 ± 0.002 by further cuts applied in our pipeline. Given the limited depth of the CANDELS data, we compared our CANDELS-based estimate for the impact of faint neighbours on the multiplicative shear measurement bias to a grid-based analysis, to which we added clustered galaxies to even fainter magnitudes based on Hubble Ultra Deep Field data, yielding a refined estimate of ~ –0.013. Our sensitivity analysis suggests that our pipeline is calibrated to an accuracy of ~0.015 once all corrections are applied, which is fully sufficient for current and near-future weak lensing studies of high-redshift clusters. As an application, we used it for a refined analysis of three highly relaxed clusters from the South Pole Telescope Sunyaev-Zeldovich survey, where we now included measurements down to the cluster core (r > 200 kpc) as enabled by our work. Compared to previously employed scales (r > 500 kpc), this tightens the cluster mass constraints by a factor 1.38 on average.
79 ASTRONOMY AND ASTROPHYSICS↗
Modeling the Transfer Function for the Dark Energy Survey
We present a forward-modeling simulation framework designed to model the data products from the Dark Energy Survey (DES). This forward-model process can be thought of as a transfer function—a mapping from cosmological/astronomical signals to the final data products used by the scientists. Using output from the cosmological simulations (the Blind Cosmology Challenge), we generate simulated images (the Ultra Fast Image Simulator) and catalogs representative of the DES data. In this work we demonstrate the framework by simulating the 244 deg 2 coadd images and catalogs in five bands for the DES Science Verification data. The simulation output is compared with the corresponding data to show that major characteristics of the images and catalogs can be captured. We also point out several directions of future improvements. Two practical examples—star-galaxy classification and proximity effects on object detection—are then used to illustrate how one can use the simulations to address systematics issues in data analysis. With clear understanding of the simplifications in our model, we show that one can use the simulations side-by-side with data products to interpret the measurements. This forward modeling approach is generally applicable for other upcoming and future surveys. It provides a powerful tool for systematics studies that is sufficiently realistic and highly controllable.
79 ASTRONOMY AND ASTROPHYSICS↗
Advancing X-ray quantum imaging through Monte-Carlo simulations
Imaging with X-rays poses fundamental limits due to radiation damage of the highly energetic photons. This becomes problematic for sensitive biological systems such as subcellular structures. Lowering the radiation dose, without sacrificing the signal-to-noise ratio, would be desirable for any kind of imaging modalities involving X-rays. To achieve this goal, quantum imaging with entangled X-ray photons constitutes a promising route. Production of biphotons have been demonstrated in the X-ray regime by the process of Spontaneous Parametric Down-Conversion (SPDC). However, compared to SPDC in the regime of visible light, the production rate for X-ray biphotons is extremely low. With the introduction of new high average brightness X-ray sources, such as 4th generation synchrotrons and high repetition rate Free-Electron X-ray Lasers (XFEL), quantum imaging may become practical. We introduce a ray tracing approach using Monte-Carlo sampling, specifically designed for quantum imaging with entangled X-ray photons generated by SPDC. By simulation, the superior image quality of quantum over classical imaging methods is demonstrated using realistic experimental conditions available at high repetition rate XFELs. With these simulations, we can efficiently assist the design of future experiments at beam lines, which can substantially accelerate the advancement of X-ray quantum imaging and reduce costs.
Entangled Photons↗
Quantifying the unknown impact of segmentation uncertainty on image-based simulations
Abstract Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.
42 ENGINEERING↗
Dark Energy Survey Year 1 results: the impact of galaxy neighbours on weak lensing cosmology with IM3SHAPE
We use a suite of simulated images based on Year 1 of the Dark Energy Survey to explore the impact of galaxy neighbours on shape measurement and shear cosmology. The hoopoe image simulations include realistic blending, galaxy positions, and spatial variations in depth and PSF properties. Using the im3shape maximum-likelihood shape measurement code, we identify four mechanisms by which neighbours can have a non-negligible influence on shear estimation. These effects, if ignored, would contribute a net multiplicative bias of $m \sim 0.03 - 0.09$ in the DES Y1 im3shape catalogue, though the precise impact will be dependent on both the measurement code and the selection cuts applied. This can be reduced to percentage level or less by removing objects with close neighbours, at a cost to the effective number density of galaxies $n_\mathrm{eff}$ of 30%. We use the cosmological inference pipeline of DES Y1 to explore the cosmological implications of neighbour bias and show that omitting blending from the calibration simulation for DES Y1 would bias the inferred clustering amplitude $S_8\equiv \sigma_8 (\Omega _\mathrm{m} /0.3)^{0.5}$ by $2 \sigma$ towards low values. Lastly, we use the hoopoe simulations to test the effect of neighbour-induced spatial correlations in the multiplicative bias. We find the impact on the recovered $S_8$ of ignoring such correlations to be subdominant to statistical error at the current level of precision.
79 ASTRONOMY AND ASTROPHYSICS↗
Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE): National Geospatial-Intelligence Agency Calibration Target Placements Field Campaign Report
The purpose of our research has been to develop a capability of simulating image collection from any sensor so as to account for sensor configuration, ground spectral radiance, and atmospheric contributions. Atmosphere is arguably the most difficult aspect of the image pipeline to accurately characterize and integrate into an image simulation. Some of the reasons for this include the highly variable spatial content of different absorbing and scattering species in the atmosphere. Additionally, the interaction of land and atmosphere can significantly impact the total radiation budget, especially in the southern Great Plains.
54 ENVIRONMENTAL SCIENCES↗
WE-EF-303-02: BEST IN PHYSICS (JOINT IMAGING- THERAPY): A Comprehensive Simulation of Image Guided Beam Gating for Liver Tumor Treatments Using Scanned Proton Therapy
Purpose: To evaluate the effectiveness of image guided beam gating for PBS liver treatments under realistic breathing conditions. Methods: We have previously proposed a Beams’ Eye View (BEV) X-ray image system as an online motion monitoring device for deriving a gating signal for PBS proton therapy. Using dedicated 4D dose calculations (4DDC), in this work we have simulated gated liver treatments using three amplitude-based gating windows (10/5/3mm) based on motion extracted from BEV imaging of fiducial markers or the diaphragm. In order to improve motion mitigation, BEV guided gating has also been combined with either volumetric (VS) or layered (LS) rescanning. Nine 4DCT(MRI) liver data-sets have been used for the investigation, which not only consider realistic patient geometries but also motion variations between different breathing cycles. All 4D plans have been quantified in terms of plan homogeneity in the PTV (D5-D95), the total estimated treatment time and the beam-on duty cycle. Results: Neither gating nor rescanning can fully retrieve a comparable plan homogeneity to the static case, and considerable reductions of the duty cycle (<10%) were observed as a Result motion variations when small gating windows are used. However, once combined with rescanning, dose homogeneity within 1% of the static plan could be achieved with reasonable prolongation of the treatment time for all 9 subjects. No differences were observed between the efficacy of layered or volumetric re-scanning, or of gating signals extracted from fiducial or diaphragm motions. However, layered rescanning may be preferred over volumetric rescanning when performed in combination with gating as it is generally more time-efficient and dosimetrically robust to patient and motion variations Conclusion Combining BEV beam gating with rescanning is an efficient and effective approach to treating mobile liver tumours, and is equally effective if either the diaphragm or fiducial markers are used as motion surrogates.
60 APPLIED LIFE SCIENCES↗
Image-related scanning systems for visual simulation.
Image related scanning systems for visual simulation display for space rendezvous and docking
IMAGING TECHNIQUE↗
Image related scanning systems for visual simulation.
Image related scanning systems for visual simulation display for space rendezvous and docking
DISPLAY SYSTEM↗
SNIa Cosmology Analysis Results from Simulated LSST Images: From Difference Imaging to Constraints on Dark Energy
Abstract The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ∼10 6 transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration Data Challenge (DC2) simulation for the Wide-Fast-Deep survey area. DC2 is the first large-scale (300 deg 2 ) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ∼15 deg 2 of DC2 over a 5 yr time span in which artificial point sources from Type Ia supernova (SNIa) light curves have been overlaid onto the images. The magnitude limits per filter are u = 23.66 mag, g = 24.69 mag, r = 24.06 mag, i = 23.45 mag, z = 22.54 mag, and y = 21.62 mag. The artifact contamination levels are ∼90% of all detections, corresponding to ∼1000 artifacts deg –2 in g band, and falling to 300 deg –2 in y band. The photometry has biases <1% for magnitudes 19.5 < m < 23. Our DIA performance on simulated images is similar to that of the Dark Energy Survey difference-imaging pipeline on real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of signal-to-noise ratio, redshift, and fitted light-curve properties. Applying a realistic SNIa cosmology analysis for redshifts z < 1, we recover the input cosmology parameters to within statistical uncertainties.
79 ASTRONOMY AND ASTROPHYSICS↗
Use of ML in Image-Based Simulation: CNNs for 3D Image Segmentation and Uncertainty Quantification.
Abstract not provided.
Roberts, Scott Alan↗
Evaluation of the soil moisture prediction accuracy of a space radar using simulation techniques
Image simulation techniques were employed to generate synthetic aperture radar images of a 17.7 km x 19.3 km test site located east of Lawrence, Kansas. The simulations were performed for a space SAR at an orbital altitude of 600 km, with the following sensor parameters: frequency = 4.75 GHz, polarization = HH, and angle of incidence range = 7 deg to 22 deg from nadir. Three sets of images were produced corresponding to three different spatial resolutions; 20 m x 20 m with 12 looks, 100 m x 100 m with 23 looks, and 1 km x 1 km with 1000 looks. Each set consisted of images for four different soil moisture distributions across the test site. Results indicate that, for the agricultural portion of the test site, the soil moisture in about 90% of the pixels can be predicted with an accuracy of = + or - 20% of field capacity. Among the three spatial resolutions, the 1 km x 1 km resolution gave the best results for most cases, however, for very dry soil conditions, the 100 m x 100 m resolution was slightly superior.
Ulaby, F. T.↗
Scale effects: HCMM data simulation. Usage of filtering techniques for scaling-up simulations
Image reduction used to simulate increase in altitude of an acquisition platform is equivalent to data smoothing, and can be achieved either by neighborhood averaging or by filtering techniques. The averaging approach is limited for accurate simulation. A filtering method is described which was based on the hypothesis that all changes due to altitude increase can be represented by a point spread function. Determination of the scale function and factor are discussed as well as implementation of the filtering. Filtering can be either in the spatial or frequency domain. In the spatial domain, filtering consists of the convolution of the image with the weights mask, and then of the declination of the points according to the appropriates scale factor. A simulation of an aircraft day image in the infrared channel is examined.
Digennaro, V.↗
Image Collection Simulation Using High-Resolution Atmospheric Modeling
A new method is described for simulating the passive remote sensing image collection of ground targets that includes effects from atmospheric physics and dynamics at fine spatial and temporal scales. The innovation in this research is the process of combining a high-resolution weather model with image collection simulation to attempt to account for heterogeneous and high-resolution atmospheric effects on image products. The atmosphere was modeled on a 3D voxel grid by a Large-Eddy Simulation (LES) driven by forcing data constrained by local ground-based and air-based observations. The spatial scale of the atmospheric model (10–100 m) came closer than conventional weather forecast scales (10–100 km) to approaching the scale of typical commercial multispectral imagery (2 m). This approach was demonstrated through a ground truth experiment conducted at the Department of Energy Atmospheric Radiation Measurement Southern Great Plains site. In this experiment, calibrated targets (colored spectral tarps) were placed on the ground, and the scene was imaged with WorldView-3 multispectral imagery at a resolution enabling the tarps to be visible in at least 9–12 image pixels. The image collection was simulated with Digital Imaging and Remote Sensing Image Generation (DIRSIG) software, using the 3D atmosphere from the LES model to generate a high-resolution cloud mask. The high-resolution atmospheric model-predicted cloud coverage was usually within 23% of the measured cloud cover. The simulated image products were comparable to the WorldView-3 satellite imagery in terms of the variations of cloud distributions and spectral properties of the ground targets in clear-sky regions, suggesting the potential utility of the proposed modeling framework in improving simulation capabilities, as well as testing and improving the operation of image collection processes.
54 ENVIRONMENTAL SCIENCES↗
Simulating realistic imaging conditions for in situ liquid microscopy
In situ transmission electron microscopy enables the imaging of biological cells, macromolecular protein complexes, nanoparticles, and other systems in a near-native environment. In order to improve interpretation of image contrast features and also predict ideal imaging conditions ahead of time, new virtual electron microscopic techniques are needed. A technique for virtual fluid-stage high-angle annular dark-field scanning transmission electron microscopy with the multislice method is presented that enables the virtual imaging of model fluid-stage systems composed of millions of atoms. The virtual technique is exemplified by simulating images of PbS nanoparticles under different imaging conditions and the results agree with previous experimental findings. General insight is obtained on the influence of the effects of fluid path length, membrane thickness, nanoparticle position, defocus and other microscope parameters on attainable image quality.
Atomistic model↗
Discriminating Accretion States via Rotational Symmetry in Simulated Polarimetric Images of M87
In 2017 April, the Event Horizon Telescope observed the shadow of the supermassive black hole at the core of the elliptical galaxy Messier 87. While the original image was constructed from measurements of the total intensity, full polarimetric data were also collected, and linear polarimetric images are expected in the near future. We propose a modal image decomposition of the linear polarization field into basis functions with varying azimuthal dependence of the electric vector position angle. We apply this decomposition to images of ray traced general relativistic magnetohydrodynamics simulations of the Messier 87 accretion disk. Furthermore, for simulated images that are physically consistent with previous observations, the magnitude of the coefficient associated with rotational symmetry, β2, is a useful discriminator between accretion states. We find that at 20 μas resolution, $| {\beta }_{2}| $ is greater than 0.2 only for models of disks with horizon-scale magnetic pressures large enough to disrupt steady accretion. We also find that images with a more radially directed electric vector position angle correspond to models with higher black hole spin. Our analysis demonstrates the utility of the proposed decomposition as a diagnostic framework to improve constraints on theoretical models.
79 ASTRONOMY AND ASTROPHYSICS↗
SU-E-J-30: Benchmark Image-Based TCP Calculation for Evaluation of PTV Margins for Lung SBRT Patients
Purpose: Tumor control probability (TCP) calculated with accumulated radiation doses may help design appropriate treatment margins. Image registration errors, however, may compromise the calculated TCP. The purpose of this study is to develop benchmark CT images to quantify registration-induced errors in the accumulated doses and their corresponding TCP. Methods: 4DCT images were registered from end-inhale (EI) to end-exhale (EE) using a “demons” algorithm. The demons DVFs were corrected by an FEM model to get realistic deformation fields. The FEM DVFs were used to warp the EI images to create the FEM-simulated images. The two images combined with the FEM DVF formed a benchmark model. Maximum intensity projection (MIP) images, created from the EI and simulated images, were used to develop IMRT plans. Two plans with 3 and 5 mm margins were developed for each patient. With these plans, radiation doses were recalculated on the simulated images and warped back to the EI images using the FEM DVFs to get the accumulated doses. The Elastix software was used to register the FEM-simulated images to the EI images. TCPs calculated with the Elastix-accumulated doses were compared with those generated by the FEM to get the TCP error of the Elastix registrations. Results: For six lung patients, the mean Elastix registration error ranged from 0.93 to 1.98 mm. Their relative dose errors in PTV were between 0.28% and 6.8% for 3mm margin plans, and between 0.29% and 6.3% for 5mm-margin plans. As the PTV margin reduced from 5 to 3 mm, the mean TCP error of the Elastix-reconstructed doses increased from 2.0% to 2.9%, and the mean NTCP errors decreased from 1.2% to 1.1%. Conclusion: Patient-specific benchmark images can be used to evaluate the impact of registration errors on the computed TCPs, and may help select appropriate PTV margins for lung SBRT patients.
60 APPLIED LIFE SCIENCES↗