Sparse image reconstruction on the sphere: a general approach with uncertainty quantification

Abstract

Inverse problems defined naturally on the sphere are becoming increasingly of interest. In this article we provide a general framework for evaluation of inverse problems on the sphere, with a strong emphasis on flexibility and scalability. We consider flexibility with respect to the prior selection (regularization), the problem definition - specifically the problem formulation (constrained/unconstrained) and problem setting (analysis/synthesis) - and optimization adopted to solve the problem. We discuss and quantify the trade-offs between problem formulation and setting. Crucially, we consider the Bayesian interpretation of the unconstrained problem which, combined with recent developments in probability density theory, permits rapid, statistically principled uncertainty quantification (UQ) in the spherical setting. Linearity is exploited to significantly increase the computational efficiency of such UQ techniques, which in some cases are shown to permit analytic solutions. We showcase this reconstruction framework and UQ techniques on a variety of spherical inverse problems. The code discussed throughout is provided under a GNU general public license, in both C++ and Python.

References in corpus (23)

• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• Spherical Needlets for CMB Data Analysis
• A novel sampling theorem on the sphere
• Wavelets, ridgelets and curvelets on the sphere
• Exact reconstruction with directional wavelets on the sphere
• Scalable splitting algorithms for big-data interferometric imaging in the SKA era
• A SURE Approach for Digital Signal/Image Deconvolution Problems
• Sparsity Averaging for Compressive Imaging
• High Resolution Weak Lensing Mass-Mapping Combining Shear and Flexion
• S2LET: A code to perform fast wavelet analysis on the sphere
• Localisation of directional scale-discretised wavelets on the sphere
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• Second-Generation Curvelets on the Sphere
• On the computation of directional scale-discretized wavelet transforms on the sphere
• Sparse image reconstruction on the sphere: implications of a new sampling theorem
• Directional spin wavelets on the sphere
• Spherical Bayesian mass-mapping with uncertainties: full sky observations on the celestial sphere
• Sparse Bayesian mass-mapping with uncertainties: local credible intervals
• Wiener filtering and pure E/B decomposition of CMB maps with anisotropic correlated noise
• Sparse image reconstruction on the sphere: analysis and synthesis
• Mapping dark matter on the celestial sphere with weak gravitational lensing
• Sparse Bayesian mass-mapping with uncertainties: peak statistics and feature locations
• Sparse Inpainting and Isotropy

Cited by in corpus (4)

• Differentiable and accelerated spherical harmonic and Wigner transforms
• Slepian Scale-Discretised Wavelets on the Sphere
• Inferring hemispheric asymmetries of stellar active regions through the information content of astrometric signals
• Bayesian variational regularization on the ball