Volume 549, January 2013
|Number of page(s)||13|
|Section||Numerical methods and codes|
|Published online||21 December 2012|
Estimating hyperparameters and instrument parameters in regularized inversion Illustration for Herschel/SPIRE map making
1 Institut d’Astrophysique de Paris (CNRS – Univ. Paris 6), 75014 Paris, France
2 Univ. Bordeaux, IMS, UMR 5218, 33400 Talence, France
3 Laboratoire des Signaux et Systèmes (CNRS – Supélec – Univ. Paris-Sud 11), 91192 Gif-sur-Yvette, France
4 Institut d’Astrophysique Spatiale (CNRS – Univ. Paris-Sud 11), 91405 Orsay, France
Received: 5 July 2012
Accepted: 19 October 2012
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 posteriordensity 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.
Key words: methods: data analysis / methods: statistical / methods: numerical / techniques: image processing
© ESO, 2013
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.