We describe regularized methods for image reconstruction and focus on the question of hyperparameter and instrument parameter estimation, i.e. unsupervised and myopic problems. We developed a Bayesian framework that is based on the \post density for all unknown quantities, given the observations. This density is explored by a Markov Chain Monte-Carlo sampling technique based on a Gibbs loop and including a Metropolis-Hastings step. The numerical evaluation relies on the SPIRE instrument of the Herschel observatory. Using simulated and real observations, we show that the hyperparameters and instrument parameters are correctly estimated, which opens up many perspectives for imaging in astrophysics.

We investigate super-resolution methods for image reconstruction from data provided by a family of scanning instruments like the Herschel observatory. To do this, we constructed a model of the instrument that faithfully reflects the physical reality, accurately taking the acquisition process into account to explain the data in a reliable manner. The inversion, ie the image reconstruction process, is based on a linear approach resulting from a quadratic regularized criterion and numerical optimization tools. The application concerns the reconstruction of maps for the SPIRE instrument of the Herschel observatory. The numerical evaluation uses simulated and real data to compare the standard tool (coaddition) and the proposed method. The inversion approach is capable to restore spatial frequencies over a bandwidth four times that possible with coaddition and thus to correctly show details invisible on standard maps. The approach is also applied to real data with significant improvement in spatial resolution.

This tutorial paper describes the problem of image reconstruction from interferometric data with a particular focus on the specific problems encountered at optical (visible/IR) wavelengths. The challenging issues in image reconstruction from interferometric data are introduced in the general framework of inverse problem approach. This framework is then used to describe existing image reconstruction algorithms in radio interferometry and the new methods specifically developed for optical interferometry.

We present an original method for reconstructing a three-dimensional object having two spatial dimensions and one spectral dimension from data provided by the infrared slit spectrograph on board the Spitzer Space Telescope. During acquisition, the light flux is deformed by a complex process comprising four main elements (the telescope aperture, the slit, the diffraction grating and optical distortion) before it reaches the two-dimensional sensor. The originality of this work lies in the physical modelling, in integral form, of this process of data formation in continuous variables. The inversion is lso approached with continuous variables in a semi-parametric format decomposing the object into a family of Gaussian functions. The estimate is built in a deterministic regularization framework as the minimizer of a quadratic criterion. These specificities give our method the power to over-resolve. Its performance is illustrated using real and simulated data. We also present a study of the resolution showing a 1.5-fold improvement relative to conventional methods.

The paper deals with the construction of images from visibilities acquired using aperture synthesis instruments: Fourier synthesis, deconvolution, and spectral interpolation/extrapolation. Its intended application is to specific situations in which the imaged object possesses two superimposed components: ($i$) an extended component together with ($ii$) a set of point sources. It is also specifically designed to the case of positive maps, and accounts for a known support. Its originality lies within joint estimation of the two components, coherently with data, properties of each component, positivity and possible support. We approach the subject as an inverse problem within a regularization framework: a regularized least-squares criterion is specifically proposed and the estimated maps are defined as its minimizer. We have investigated several options for the numerical minimization and we propose a new efficient algorithm based on augmented Lagrangian. Evaluation is carried out using simulated and real data (from radio interferometry) demonstrating the capability to accurately separate the two components.