Bregman geometry-aware split Gibbs sampling for Bayesian Poisson inverse problems

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Faye, Elhadji Cisse, Fall, Mame Diarra, Dobigeon, Nicolas, Barat, Eric
Formato: Preprint
Publicado: 2025
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866918203427389440
author Faye, Elhadji Cisse
Fall, Mame Diarra
Dobigeon, Nicolas
Barat, Eric
author_facet Faye, Elhadji Cisse
Fall, Mame Diarra
Dobigeon, Nicolas
Barat, Eric
contents This paper proposes a novel Bayesian framework for solving Poisson inverse problems by devising a Monte Carlo sampling algorithm which accounts for the underlying non-Euclidean geometry. To address the challenges posed by the Poisson likelihood -- such as non-Lipschitz gradients and positivity constraints -- we derive a Bayesian model which leverages exact and asymptotically exact data augmentations. In particular, the augmented model incorporates two sets of splitting variables both derived through a Bregman divergence based on the Burg entropy. Interestingly the resulting augmented posterior distribution is characterized by conditional distributions which benefit from natural conjugacy properties and preserve the intrinsic geometry of the latent and splitting variables. This allows for efficient sampling via Gibbs steps, which can be performed explicitly for all conditionals, except the one incorporating the regularization potential. For this latter, we resort to a Hessian Riemannian Langevin Monte Carlo (HRLMC) algorithm which is well suited to handle priors with explicit or easily computable score functions. By operating on a mirror manifold, this Langevin step ensures that the sampling satisfies the positivity constraints and more accurately reflects the underlying problem structure. Performance results obtained on denoising, deblurring, and positron emission tomography (PET) experiments demonstrate that the method achieves competitive performance in terms of reconstruction quality compared to optimization- and sampling-based approaches.
format Preprint
id arxiv_https___arxiv_org_abs_2511_12257
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Bregman geometry-aware split Gibbs sampling for Bayesian Poisson inverse problems
Faye, Elhadji Cisse
Fall, Mame Diarra
Dobigeon, Nicolas
Barat, Eric
Computation
Computer Vision and Pattern Recognition
Image and Video Processing
Machine Learning
This paper proposes a novel Bayesian framework for solving Poisson inverse problems by devising a Monte Carlo sampling algorithm which accounts for the underlying non-Euclidean geometry. To address the challenges posed by the Poisson likelihood -- such as non-Lipschitz gradients and positivity constraints -- we derive a Bayesian model which leverages exact and asymptotically exact data augmentations. In particular, the augmented model incorporates two sets of splitting variables both derived through a Bregman divergence based on the Burg entropy. Interestingly the resulting augmented posterior distribution is characterized by conditional distributions which benefit from natural conjugacy properties and preserve the intrinsic geometry of the latent and splitting variables. This allows for efficient sampling via Gibbs steps, which can be performed explicitly for all conditionals, except the one incorporating the regularization potential. For this latter, we resort to a Hessian Riemannian Langevin Monte Carlo (HRLMC) algorithm which is well suited to handle priors with explicit or easily computable score functions. By operating on a mirror manifold, this Langevin step ensures that the sampling satisfies the positivity constraints and more accurately reflects the underlying problem structure. Performance results obtained on denoising, deblurring, and positron emission tomography (PET) experiments demonstrate that the method achieves competitive performance in terms of reconstruction quality compared to optimization- and sampling-based approaches.
title Bregman geometry-aware split Gibbs sampling for Bayesian Poisson inverse problems
topic Computation
Computer Vision and Pattern Recognition
Image and Video Processing
Machine Learning
url https://arxiv.org/abs/2511.12257