Bayesian imaging inverse problem with scattering transform

Abstract

Bayesian imaging inverse problems in astrophysics and cosmology remain challenging, particularly in low-data regimes, due to complex forward operators and the frequent lack of well-motivated priors for non-Gaussian signals. In this paper, we introduce a Bayesian approach that addresses these difficulties by relying on a low-dimensional representation of physical fields built from Scattering Transform statistics. This representation enables inference to be performed in a compact model space, where we recover a posterior distribution over signal models that are consistent with the observed data. We propose an iterative adaptive algorithm to efficiently approximate this posterior distribution. We apply our method to a large-scale structure column density field from the Quijote simulations, using a realistic instrumental forward operator. We demonstrate both accurate statistical inference and deterministic signal reconstruction from a single contaminated image, without relying on any external prior distribution for the field of interest. These results demonstrate that Scattering Transform statistics provide an effective representation for solving complex imaging inverse problems in challenging low-data regimes. Our approach opens the way to new applications for non-Gaussian astrophysical and cosmological signals for which little or no prior modeling is available.

5+2 pages, 4+1 figures. Submitted to A&A. Comments welcome

References in corpus (13)

• The Quijote simulations
• Component separation methods for the Planck mission
• A parsec-scale Galactic 3D dust map out to 1.25 kpc from the Sun
• New Interpretable Statistics for Large Scale Structure Analysis and Generation
• The RWST, a comprehensive statistical description of the non-Gaussian structures in the ISM
• Probabilistic Mass Mapping with Neural Score Estimation
• A new approach for the statistical denoising of Planck interstellar dust polarization data
• Non-Gaussian modelling and statistical denoising of Planck dust polarization full-sky maps using scattering transforms
• Single frequency CMB B-mode inference with realistic foregrounds from a single training image
• Scattering Spectra Models for Physics
• Generative models of astrophysical fields with scattering transforms on the sphere
• Non-parametric Bayesian reconstruction of Galactic magnetic fields using Information Field Theory: The inclusion of line-of-sight information in ultra-high energy cosmic ray backtracking
• Comparing the morphology of molecular clouds without supervision

Cited by in corpus (1)

• Fast Sampling of Cosmological Initial Conditions with Gaussian Neural Posterior Estimation