Off-the-grid regularisation for Poisson inverse problems

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Lazzaretti, Marta, Estatico, Claudio, Melero, Alejandro, Calatroni, Luca
Format: Preprint
Publié: 2024
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866913790144020480
author Lazzaretti, Marta
Estatico, Claudio
Melero, Alejandro
Calatroni, Luca
author_facet Lazzaretti, Marta
Estatico, Claudio
Melero, Alejandro
Calatroni, Luca
contents Off-the-grid regularisation has been extensively employed over the last decade in the context of ill-posed inverse problems formulated in the continuous setting of the space of Radon measures $\mathcal{M}(\mathcal{X})$. These approaches enjoy convexity and counteract the discretisation biases as well the numerical instabilities typical of their discrete counterparts. In the framework of sparse reconstruction of discrete point measures (sum of weighted Diracs), a Total Variation regularisation norm in $\mathcal{M}(\mathcal{X})$ is typically combined with an $L^2$ data term modelling additive Gaussian noise. To asses the framework of off-the-grid regularisation in the presence of signal-dependent Poisson noise, we consider in this work a variational model coupling the Total Variation regularisation with a Kullback-Leibler data term under a non-negativity constraint. Analytically, we study the optimality conditions of the composite functional and analyse its dual problem. Then, we consider an homotopy strategy to select an optimal regularisation parameter and use it within a Sliding Frank-Wolfe algorithm. Several numerical experiments on both 1D/2D simulated and real 3D fluorescent microscopy data are reported.
format Preprint
id arxiv_https___arxiv_org_abs_2404_00810
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Off-the-grid regularisation for Poisson inverse problems
Lazzaretti, Marta
Estatico, Claudio
Melero, Alejandro
Calatroni, Luca
Numerical Analysis
Optimization and Control
Off-the-grid regularisation has been extensively employed over the last decade in the context of ill-posed inverse problems formulated in the continuous setting of the space of Radon measures $\mathcal{M}(\mathcal{X})$. These approaches enjoy convexity and counteract the discretisation biases as well the numerical instabilities typical of their discrete counterparts. In the framework of sparse reconstruction of discrete point measures (sum of weighted Diracs), a Total Variation regularisation norm in $\mathcal{M}(\mathcal{X})$ is typically combined with an $L^2$ data term modelling additive Gaussian noise. To asses the framework of off-the-grid regularisation in the presence of signal-dependent Poisson noise, we consider in this work a variational model coupling the Total Variation regularisation with a Kullback-Leibler data term under a non-negativity constraint. Analytically, we study the optimality conditions of the composite functional and analyse its dual problem. Then, we consider an homotopy strategy to select an optimal regularisation parameter and use it within a Sliding Frank-Wolfe algorithm. Several numerical experiments on both 1D/2D simulated and real 3D fluorescent microscopy data are reported.
title Off-the-grid regularisation for Poisson inverse problems
topic Numerical Analysis
Optimization and Control
url https://arxiv.org/abs/2404.00810