Decomposition and Classification of Spectral Lines in Astronomical Radio Data Cubes

The natural output of imaging spectroscopy in astronomy is a 3D data cube with two spatial and one frequency axis. The spectrum of each image pixel consists of an emission line which is Doppler-shifted by gas motions along the line of sight. These data are essential to understand the gas distribution and kinematics of the astronomical object. We propose a two-step method to extract coherent kinematic structures from the data cube. First, the spectra are decomposed into a sum of Gaussians using a Bayesian method to obtain an estimation of spectral lines. Second, we aim at tracking the estimated lines to get an estimation of the structures in the cube. The performance of the approach is evaluated on a real radio-astronomical observation.

[1]  L. Devroye.  Non-Uniform Random Variate Generation , 1986 .

[2]  D. Brie,et al.  Simulation of postive normal variables using several proposal distributions , 2005, IEEE/SP 13th Workshop on Statistical Signal Processing, 2005.

[3]  H. B. Mitchell.  Markov Random Fields , 1982 .

[4]  Vincent Mazet.  Développement de méthodes de traitement de signaux spectroscopiques : estimation de la ligne de base et du spectre de raies , 2005 .

[5]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[6]  Anil K. Jain,et al.  Markov random fields : theory and application , 1993 .

[7]  J. Idier.  Bayesian Approach to Inverse Problems: Idier/Bayesian , 2010 .

[8]  P. Green.  Reversible jump Markov chain Monte Carlo computation and Bayesian model determination , 1995 .

[9]  Christophe Collet,et al.  Multiband Segmentation of a Spectroscopic Line Data Cube: Application to the HI Data Cube of the Spiral Galaxy NGC 4254 , 2005, EURASIP J. Adv. Signal Process..

[10]  Adrian F. M. Smith,et al.  Bayesian Statistics 5. , 1998 .

[11]  O. Cappé,et al.  Reversible jump, birth‐and‐death and more general continuous time Markov chain Monte Carlo samplers , 2003 .

[12]  Peter J. Rousseeuw,et al.  Robust Regression and Outlier Detection , 2005, Wiley Series in Probability and Statistics.

[13]  Christian P. Robert,et al.  Monte Carlo Statistical Methods (Springer Texts in Statistics) , 2005 .

[14]  Lee G. Mundy,et al.  NGC 4254: A Spiral Galaxy with an M = 1 Mode and Infalling Gas , 1993 .

[15]  Hichem Snoussi,et al.  Blind separation of noisy Gaussian stationary sources. Application to cosmic microwave background imaging , 2002, 2002 11th European Signal Processing Conference.