-
Estimating hyperparameters and instrument parameters in regularized inversion. Illustration for SPIRE/Herschel map making
Authors:
F. Orieux,
J. -F. Giovannelli,
T. Rodet,
A. Abergel
Abstract:
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 in…
▽ More
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.
△ Less
Submitted 15 November, 2012;
originally announced November 2012.
-
Super-resolution in map-making based on a physical instrument model and regularized inversion. Application to SPIRE/Herschel
Authors:
F. Orieux,
J. -F. Giovannelli,
T. Rodet,
A. Abergel,
H. Ayasso,
M. Husson
Abstract:
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,…
▽ More
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.
△ Less
Submitted 23 December, 2011; v1 submitted 21 March, 2011;
originally announced March 2011.
-
Data inversion for over-resolved spectral imaging in astronomy
Authors:
T. Rodet,
F. Orieux,
J. -F. Giovannelli,
A. Abergel
Abstract:
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 disto…
▽ More
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.
△ Less
Submitted 11 February, 2009;
originally announced February 2009.
-
Positive deconvolution for superimposed extended source and point sources
Authors:
J. -F. Giovannelli,
A. Coulais
Abstract:
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 specific…
▽ More
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.
△ Less
Submitted 29 July, 2005;
originally announced July 2005.