Simulating image coaddition with the Nancy Grace Roman Space Telescope: III. Software improvements and new linear algebra strategies

Abstract

The Nancy Grace Roman Space Telescope will implement a devoted weak gravitational lensing program with its High Latitude Wide Area Survey. For cosmological purposes, a critical step in Roman image processing is to combine dithered undersampled images into unified oversampled images and thus enable high-precision shape measurements. IMCOM is an image coaddition algorithm which offers control over point spread functions in output images. This paper presents the refactored IMCOM software, featuring full object-oriented programming, improved data structures, and alternative linear algebra strategies for determining coaddition weights. Combining these improvements and other acceleration measures, to produce almost equivalent coadded images, the consumption of core-hours has been reduced by about an order of magnitude. We then re-coadd a $16\times 16{\mathrm{arcmin}}^2$ region of our previous image simulations with three linear algebra kernels in four bands, and compare the results in terms of IMCOM optimization goals, properties of coadded noise frames, and measurements of simulated stars. The Cholesky kernel is efficient and relatively accurate, yet its irregular windows for input pixels slightly bias coaddition results. The iterative kernel avoids this issue by tailoring input pixel selection for each output pixel; it yields better noise control, but can be limited by random errors due to finite tolerance. The empirical kernel coadds images using an empirical relation based on geometry; it is inaccurate, but being much faster, it provides a valid option for "quick look" purposes. We fine-tune IMCOM hyperparameters in a companion paper.

38 pages, 22 figures, matches version accepted by ApJSupp

References in corpus (38)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Astropy: A Community Python Package for Astronomy
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Weak Gravitational Lensing
• Observational Probes of Cosmic Acceleration
• Drizzle: A Method for the Linear Reconstruction of Undersampled Images
• Euclid. I. Overview of the Euclid mission
• Cosmology from cosmic shear power spectra with Subaru Hyper Suprime-Cam first-year data
• Cosmology with cosmic shear observations: a review
• Dark Energy Survey Year 3 Results: Cosmology from Cosmic Shear and Robustness to Data Calibration
• Dark Energy Survey Year 3 Results: Cosmology from Cosmic Shear and Robustness to Modeling Uncertainty
• Shear calibration biases in weak lensing surveys
• Systematic errors in weak lensing: application to SDSS galaxy-galaxy weak lensing
• Cosmological constraints from cosmic shear two-point correlation functions with HSC survey first-year data
• Practical Weak Lensing Shear Measurement with Metacalibration
• CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST
• Combining Undersampled Dithered Images
• The LSST DESC DC2 Simulated Sky Survey
• Properties and characteristics of the WFIRST H4RG-10 detectors
• KiDS-1000: cosmic shear with enhanced redshift calibration
• Mitigating Shear-dependent Object Detection Biases with Metacalibration
• Optimal Linear Image Combination
• KiDS-1000: Cosmology with improved cosmic shear measurements
• A Joint Roman Space Telescope and Rubin Observatory Synthetic Wide-Field Imaging Survey
• Discovery of a red backsplash galaxy candidate near M81
• Analytical Weak-lensing Shear Responses of Galaxy Properties and Galaxy Detection
• Validating Synthetic Galaxy Catalogs for Dark Energy Science in the LSST Era
• A General Framework for Removing Point Spread Function Additive Systematics in Cosmological Weak Lensing Analysis
• Impact of Point Spread Function Higher Moments Error on Weak Gravitational Lensing II: A Comprehensive Study
• Simulating image coaddition with the Nancy Grace Roman Space Telescope: I. Simulation methodology and general results
• PSFs of coadded images
• A Differentiable Perturbation-based Weak Lensing Shear Estimator
• Simulating image coaddition with the Nancy Grace Roman Space Telescope: II. Analysis of the simulated images and implications for weak lensing
• Analysis of biasing from noise from the Nancy Grace Roman Space Telescope: implications for weak lensing

Cited by in corpus (2)

• Chromatic Effects on the PSF and Shear Measurement for the Roman Space Telescope High-Latitude Wide Area Survey
• Charge diffusion and modulation transfer function in a Nancy Grace Roman Space Telescope detector