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 \,{\rm 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.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to effectively combine undersampled images from the Nancy Grace Roman Space Telescope to generate unified oversampled images, thereby achieving high - precision shape measurements, especially for the purpose of weak gravitational lensing research. Specifically, the paper focuses on improving the Imcom software to enhance computational efficiency, optimize data structures, and introduce new linear - algebra strategies to determine image - stacking weights. ### Problem Background Weak gravitational lensing is an important detection method in cosmology, which can reveal the mass distribution in the universe. In order to accurately measure the shapes of galaxies, it is necessary to combine multiple undersampled images into an oversampled image. This step is particularly important for the High Latitude Wide Area Survey (HLWAS) of the Roman Space Telescope. ### Key Challenges 1. **Computational Efficiency**: The computational cost of the original Imcom implementation is too high to be applied to the entire HLWAS area. 2. **Boundary Effects**: Imcom may produce systematic errors when processing image boundaries, affecting source identification or shape measurement. 3. **Complex Linear System Solving**: The complexity of the existing linear system solving is O(𝑛³), which is unbearable for deep - field image processing. ### Solutions The paper proposes a series of improvement measures, including: - **Reconstructing the Imcom Software Framework**: Adopting object - oriented programming (OOP), optimizing data structures, and providing additional output maps for diagnostic results. - **Introducing New Linear Algebra Strategies**: Proposing three different kernel functions (Cholesky kernel, Iterative kernel, and Empirical kernel) to determine image - stacking weights. - **Acceleration Measures**: Reducing the core - hour consumption through various technical means, increasing the computational efficiency by approximately one order of magnitude. ### Goals The ultimate goal is to enable Imcom to significantly reduce the demand for computational resources while maintaining high precision through these improvements, thereby providing strong support for future weak gravitational lensing research. ### Formula Summary The key formulas involved in the paper are as follows: - The linear transformation formula of the output image: \[ H_\alpha = \sum_{i = 0}^{n - 1} T_{\alpha i}I_i \] where $I_i$ is the input image, $H_\alpha$ is the output image, and $T_{\alpha i}$ is the transformation matrix. - Definitions of PSF leakage metric $U_\alpha$ and noise amplification metric $\Sigma_\alpha$: \[ U_\alpha=\|\text{PSF}_{\alpha,\text{out}}-\Gamma\|^2 \] \[ \Sigma_\alpha=\sum_{i,j}N_{ij}T_{\alpha i}T_{\alpha j} \] - Solution of the linear system: \[ T_{\alpha i}=\sum_j[(A + \kappa_\alpha I_{n\times n})^{-1}]_{ij}\left(-\frac{1}{2}B_{\alpha j}\right) \] Through these improvements and the application of new methods, the paper aims to provide a more efficient and accurate solution for the image processing of the Roman Space Telescope.

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

Simulating image coaddition with the Nancy Grace Roman Space Telescope: I. Simulation methodology and general results

Simulating image coaddition with the Nancy Grace Roman Space Telescope – I. Simulation methodology and general results

Simulating image coaddition with the Nancy Grace Roman Space Telescope: II. Analysis of the simulated images and implications for weak lensing

Simulating image coaddition with the Nancy Grace Roman Space Telescope – II. Analysis of the simulated images and implications for weak lensing

Simultaneous Coded Aperture and Dictionary Optimization in Compressive Spectral Imaging Via Coherence Minimization.

J-comb: an Image Fusion Algorithm to Combine Observations Covering Different Spatial Frequency Ranges

Analysis of biasing from noise from the Nancy Grace Roman Space Telescope: implications for weak lensing

Dazzle: Oversampled Image Reconstruction and Difference-Imaging Photometry for the Nancy Grace Roman Space Telescope

How to coadd images: II. Anti-aliasing and PSF deconvolution

A Synthetic Roman Space Telescope High-Latitude Imaging Survey: Simulation Suite and the Impact of Wavefront Errors on Weak Gravitational Lensing

SIMULATION OF ASTRONOMICAL IMAGES FROM OPTICAL SURVEY TELESCOPES USING A COMPREHENSIVE PHOTON MONTE CARLO APPROACH

The little coadd that could: Estimating shear from coadded images

Quantum yield and charge diffusion in the Nancy Grace Roman Space Telescope infrared detectors

Computation Theory of Large-Scale Partially Coherent Imaging by the Modified Modal Expansion Method

NANCY: Next-generation All-sky Near-infrared Community surveY

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

STIPS: The Nancy Grace Roman Space Telescope Imaging Product Simulator

Quantifying Roman WFI Dark Images with the Wavelet Scattering Transform

Monte-Carlo simulation of stellar intensity interferometry