Skip to main content

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

Author(s): Melchior, Peter M; Moolekamp, F; Jerdee, M; Armstrong, R; Sun, A-L; et al

To refer to this page use:
Full metadata record
DC FieldValueLanguage
dc.contributor.authorMelchior, Peter M-
dc.contributor.authorMoolekamp, F-
dc.contributor.authorJerdee, M-
dc.contributor.authorArmstrong, R-
dc.contributor.authorSun, A-L-
dc.contributor.authorBosch, J-
dc.contributor.authorLupton, R-
dc.identifier.citationMelchior, P, Moolekamp, F, Jerdee, M, Armstrong, R, Sun, A-L, Bosch, J, Lupton, R. (2018). SCARLET: Source separation in multi-band images by Constrained Matrix Factorization. ASTRONOMY AND COMPUTING, 24 (129 - 142. doi:10.1016/j.ascom.2018.07.001en_US
dc.description.abstractWe 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 Strategic 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
dc.format.extent129 - 142en_US
dc.relation.ispartofASTRONOMY AND COMPUTINGen_US
dc.rightsAuthor's manuscripten_US
dc.subjectmethods: data analysis, techniques: image processing, galaxies: structure, Non-negative matrix factorizationen_US
dc.titleSCARLET: Source separation in multi-band images by Constrained Matrix Factorizationen_US
dc.typeJournal Articleen_US

Files in This Item:
File Description SizeFormat 
1802.10157.pdf1.43 MBAdobe PDFView/Download

Items in OAR@Princeton are protected by copyright, with all rights reserved, unless otherwise indicated.