Advanced search
Start date
Betweenand
Related content
(Reference retrieved automatically from Web of Science through information on FAPESP grant and its corresponding number as mentioned in the publication by the authors.)

Spatial field reconstruction with INLA: application to IFU galaxy data

Full text
Author(s):
Gonzalez-Gaitan, S. [1] ; de Souza, R. S. [2] ; Krone-Martins, A. [3] ; Cameron, E. [4] ; Coelho, P. [5] ; Galbany, L. [6] ; Ishida, E. E. O. [7] ; Collaboration, COIN
Total Authors: 8
Affiliation:
[1] Univ Lisbon, Inst Super Tecn, CENTRA COSTAR, Ave Rovisco Pais 1, P-1049001 Lisbon - Portugal
[2] Univ N Carolina, Dept Phys & Astron, Chapel Hill, NC 27599 - USA
[3] Univ Lisbon, Fac Ciencias, CENTRA SIM, Ed C8, P-1749016 Lisbon - Portugal
[4] Univ Oxford, Li Ka Shing Ctr Hlth Informat & Discovery, Big Data Inst, Old Rd Campus, Oxford OX3 7LS - England
[5] Univ Sao Paulo, Inst Astron Geofis & Ciencias Atmosfer, BR-05508090 Sao Paulo, SP - Brazil
[6] Univ Pittsburgh, Dept Phys & Astron, PITT PACC, Pittsburgh, PA 15260 - USA
[7] Univ Clermont Auvergne, CNRS IN2P3, LPC, F-63000 Clermont Ferrand - France
Total Affiliations: 7
Document type: Journal article
Source: Monthly Notices of the Royal Astronomical Society; v. 482, n. 3, p. 3880-3891, JAN 2019.
Web of Science Citations: 1
Abstract

Astronomical observations of extended sources, such as cubes of integral field spectroscopy (IFS), encode autocorrelated spatial structures that cannot be optimally exploited by standard methodologies. This work introduces a novel technique to model IFS data sets, which treats the observed galaxy properties as realizations of an unobserved Gaussian Markov random field. The method is computationally efficient, resilient to the presence of low-signal-to-noise regions, and uses an alternative to Markov Chain Monte Carlo for fast Bayesian inference - the Integrated Nested Laplace Approximation. As a case study, we analyse 721 IFS data cubes of nearby galaxies from the CALIFA and PISCO surveys, for which we retrieve the maps of the following physical properties: age, metallicity, mass, and extinction. The proposed Bayesian approach, built on a generative representation of the galaxy properties, enables the creation of synthetic images, recovery of areas with bad pixels, and an increased power to detect structures in data sets subject to substantial noise and/or sparsity of sampling. A snippet code to reproduce the analysis of this paper is available in the COIN toolbox, together with the field reconstructions of the CALIFA and PISCO samples. (AU)

FAPESP's process: 09/54006-4 - A computer cluster for the Astronomy Department of the University of São Paulo Institute of Astronomy, Geophysics and Atmospheric Sciences and for the Cruzeiro do Sul University Astrophysics Center
Grantee:Elisabete Maria de Gouveia Dal Pino
Support Opportunities: Multi-user Equipment Program