SCARLET: Source separation in multi-band images by Constrained Matrix Factorization

TL;DR

SCARLET is a flexible source separation framework for multi-band images that generalizes non-negative matrix factorization with constraints, enabling improved deblending and photometry in crowded astronomical scenes.

Contribution

It introduces a novel constrained matrix factorization approach for multi-band image analysis, incorporating symmetry, monotonicity, and PSF convolution handling, with demonstrated superior performance over existing methods.

Findings

SCARLET outperforms the HSC-SDSS deblender in flux, color, and morphology recovery.

It effectively handles correlated noise and variable seeing conditions.

The method is applicable to crowded extragalactic scenes and AGN host galaxy separation.

Abstract

We present the source separation framework SCARLET for multi-band images, which is based on a generalization of the Non-negative Matrix Factorization to alternative and several simultaneous constraints. Our approach describes the observed scene as a mixture of components with compact spatial support and uniform spectra over their support. We present the algorithm to perform the matrix factorization and introduce constraints that are useful for optical images of stars and distinct stellar populations in galaxies, in particular symmetry and monotonicity with respect to the source peak position. We also derive the treatment of correlated noise and convolutions with band-dependent point spread functions, rendering our approach applicable to coadded images observed under variable seeing conditions. SCARLET thus yields a PSF-matched photometry measurement with an optimally chosen weight function given by the mean morphology in all available bands. We demonstrate the performance of SCARLET for deblending crowded extragalactic scenes and on an AGN jet -- host galaxy separation problem in deep 5-band imaging from the Hyper Suprime-Cam Stategic Survey Program. Using simulations with prominent crowding we show that SCARLET yields superior results to the HSC-SDSS deblender for the recovery of total fluxes, colors, and morphologies. Due to its non-parametric nature, a conceptual limitation of SCARLET is its sensitivity to undetected sources or multiple stellar population within detected sources, but an iterative strategy that adds components at the location of significant residuals appears promising. The code is implemented in Python with C++ extensions and is available at https://github.com/fred3m/scarlet