Joint survey processing: combined resampling and convolution for galaxy modelling and deblending

Abstract

We present an extension of the multi-band galaxy fitting method scarlet which allows the joint modeling of astronomical images from different instruments, by performing simultaneous resampling and convolution. We introduce a fast and formally accurate linear projection operation that maps a pixelated model at a given resolution onto an observation frame with a different point spread function and pixel scale. We test our implementation against the well-tested resampling and convolution method in galsim on simulated images mimicking observations with the Euclid space telescope and the Vera C. Rubin Observatory, and find that it reduces interpolation errors by an order of magnitude or more compared to galsim default settings. Tests with a wide range of levels of blending show more accurate galaxy models from joint modeling of Euclid and Rubin images compared to separate modeling of each survey by up to an order of magnitude. Our results demonstrate, for the first time, the feasibility and utility of performing non-parametric pixel-level data fusion of overlapping imaging surveys. All results can be reproduced with the specific versions of the codes and notebooks used in this study.

arXiv commentSubmitted to Astronomy and Computing

References in corpus (17)

• The NumPy array: a structure for efficient numerical computation
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• Wide-Field InfrarRed Survey Telescope-Astrophysics Focused Telescope Assets WFIRST-AFTA 2015 Report
• Euclid Definition Study Report
• The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s
• SCARLET: Source separation in multi-band images by Constrained Matrix Factorization
• Scientific Synergy Between LSST and Euclid
• Cosmology with the Wide-Field Infrared Survey Telescope -- Synergies with the Rubin Observatory Legacy Survey of Space and Time
• Deblending galaxies with Variational Autoencoders: a joint multi-band, multi-instrument approach
• A Synthetic Roman Space Telescope High-Latitude Imaging Survey: Simulation Suite and the Impact of Wavefront Errors on Weak Gravitational Lensing
• Enhancing LSST Science with Euclid Synergy
• Proximal Adam: Robust Adaptive Update Scheme for Constrained Optimization
• Joint Survey Processing of Euclid, Rubin and Roman: Final Report
• LSST Cadence Optimization White Paper in Support of Observations of Unresolved Tidal Stellar Streams in Galaxies beyond the Local Group
• Galaxy shape measurement synergies between LSST and Euclid
• Joint Survey Processing of LSST, Euclid and WFIRST: Enabling a broad array of astrophysics and cosmology through pixel level combinations of datasets
• Simultaneous LSST and Euclid observations - advantages for Solar System Objects

Cited by in corpus (4)

• A Joint Roman Space Telescope and Rubin Observatory Synthetic Wide-Field Imaging Survey
• Machine Learning and Cosmology
• OpenUniverse2024: A shared, simulated view of the sky for the next generation of cosmological surveys
• Neural Networks for cosmological model selection and feature importance using Cosmic Microwave Background data