Optimal Image Reconstruction in Radio Interferometry

Abstract

We introduce a method for analyzing radio interferometry data which produces maps which are optimal in the Bayesian sense of maximum posterior probability density, given certain prior assumptions. It is similar to maximum entropy techniques, but with an exact accounting of the multiplicity instead of the usual approximation involving Stirling's formula. It also incorporates an Occam factor, automatically limiting the effective amount of detail in the map to that justified by the data. We use Gibbs sampling to determine, to any desired degree of accuracy, the multi-dimensional posterior density distribution. From this we can construct a mean posterior map and other measures of the posterior density, including confidence limits on any well-defined function of the posterior map.

41 pages, 11 figures. High resolution figures 8 and 9 available at http://www.astro.uiuc.edu/~bwandelt/SuttonWandelt2006

References in corpus (4)

• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• A Bayesian approach to discrete object detection in astronomical datasets
• Short Spacing Synthesis from a Primary Beam Scanned Interferometer
• Deconvolution Map-Making for Cosmic Microwave Background Observations

Cited by in corpus (19)

• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization
• Imaging the Supermassive Black Hole Shadow and Jet Base of M87 with the Event Horizon Telescope
• Monte-Carlo Imaging for Optical Interferometry
• A super-resolution analysis of the DSHARP survey: Substructure is common in the inner 30 au
• CHARA/MIRC observations of two M supergiants in Perseus OB1: temperature, Bayesian modeling, and compressed sensing imaging
• Multi-GPU maximum entropy image synthesis for radio astronomy
• Matrix-free Large Scale Bayesian inference in cosmology
• Probabilistic image reconstruction for radio interferometers
• Bayesian Inference for Radio Observations
• Point Process Algorithm: A New Bayesian Approach for Planet Signal Extraction with the Terrestrial Planet Finder
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Bayesian angular power spectrum analysis of interferometric data
• Detecting Unspecified Structure in Low-Count Images
• Bayesian Image Reconstruction Based on Voronoi Diagrams
• AGILIS: Agile Guided Interferometer for Longbaseline Imaging Synthesis - Demonstration and concepts of reconfigurable optical imaging interferometers
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• Further evaluation of bootstrap resampling as a tool for radio-interferometric imaging fidelity assessment
• Variable structure in the PDS70 disc and uncertainties in radio-interferometric image restoration