Saved in:
Bibliographic Details
Main Authors: Pierre, Sébastien, Allys, Erwan, Richard, Pablo, Soletskyi, Roman, Tsouros, Alexandros
Format: Preprint
Published: 2026
Subjects:
Online Access:https://arxiv.org/abs/2602.05816
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866908816067526656
author Pierre, Sébastien
Allys, Erwan
Richard, Pablo
Soletskyi, Roman
Tsouros, Alexandros
author_facet Pierre, Sébastien
Allys, Erwan
Richard, Pablo
Soletskyi, Roman
Tsouros, Alexandros
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2602_05816
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Bayesian imaging inverse problem with scattering transform
Pierre, Sébastien
Allys, Erwan
Richard, Pablo
Soletskyi, Roman
Tsouros, Alexandros
Instrumentation and Methods for Astrophysics
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.
title Bayesian imaging inverse problem with scattering transform
topic Instrumentation and Methods for Astrophysics
url https://arxiv.org/abs/2602.05816