Hierarchical multispectral galaxy decomposition using a MCMC algorithm with multiple temperature simulated annealing

  • Authors:
  • Benjamin Perret;Vincent Mazet;Christophe Collet;íric Slezak

  • Affiliations:
  • LSIIT (UMR 7005, University of Strasbourg-CNRS), Pôle API, Bd Sébastien Brant, BP 10413, 67412 Illkirch, Cedex, France;LSIIT (UMR 7005, University of Strasbourg-CNRS), Pôle API, Bd Sébastien Brant, BP 10413, 67412 Illkirch, Cedex, France;LSIIT (UMR 7005, University of Strasbourg-CNRS), Pôle API, Bd Sébastien Brant, BP 10413, 67412 Illkirch, Cedex, France;University of Nice-Sophia Antipolis, CNRS, Observatoire de la Cote d'Azur, Cassiopee Laboratory, Boulevard de l'Observatoire, BP 4229, 06304 Nice, Cedex 4, France

  • Venue:
  • Pattern Recognition
  • Year:
  • 2011

Quantified Score

Hi-index 0.01

Visualization

Abstract

We present a new method for the parametric decomposition of barred spiral galaxies in multispectral observations. The observation is modelled with a realistic image formation model and the galaxy is composed of physically significant parametric structures. The model also includes a parametric filtering to remove non-desirable aspects of the observation. Both the model and the filter parameters are estimated by a robust Monte Carlo Markov chain (MCMC) algorithm. The algorithm is based on a Gibbs sampler combined with a novel strategy of simulated annealing in which several temperatures allow to manage efficiently the simulation effort. Besides, the overall decomposition is performed following an original framework: a hierarchy of models from a coarse model to the finest one is defined. At each step of the hierarchy the estimate of a coarse model is used to initialize the estimation of the finer model. This leads to an unsupervised decomposition scheme with a reduced computation time. We have validated the method on simulated and real 5-band images: the results showed the accuracy and the robustness of the proposed approach.