-
ROHSA: Regularized Optimization for Hyper-Spectral Analysis - Application to phase separation of 21 cm data
Authors:
Antoine Marchal,
Marc-Antoine Miville-Deschenes,
Francois Orieux,
Nicolas Gac,
Charles Soussen,
Marie-Jeanne Lesot,
Adrien Revault d'Allonnes,
Quentin Salome
Abstract:
Extracting the multiphase structure of the neutral interstellar medium (ISM) is key to understand the star formation in galaxies. The radiative condensation of the diffuse warm neutral medium producing a thermally unstable lukewarm medium and a dense cold medium is closely related to the initial step which leads the HI-to-H2 transition and the formation of molecular clouds. Up to now the mapping o…
▽ More
Extracting the multiphase structure of the neutral interstellar medium (ISM) is key to understand the star formation in galaxies. The radiative condensation of the diffuse warm neutral medium producing a thermally unstable lukewarm medium and a dense cold medium is closely related to the initial step which leads the HI-to-H2 transition and the formation of molecular clouds. Up to now the mapping of these phases out of 21 cm emission hyper-spectral cubes has remained elusive mostly due to the velocity blending of individual cold structures present on a given line of sight.
To address this problem, we developed a new Gaussian decomposition algorithm named ROHSA (Regularized Optimization for Hyper-Spectral Analysis) based on a multi-resolution process from coarse to fine grid. ROHSA uses a regularized non-linear least-square criterion to take into account simultaneously the spatial coherence of the emission and the multiphase nature of the gas. In order to obtain a solution with spatially smooth parameters, the optimization is performed on the whole data cube at once. The performances of ROHSA were tested on a synthetic observation computed from numerical simulations of thermally bi-stable turbulence. An application on a 21 cm observation of a high Galactic latitude region from the GHIGLS survey is presented.
The evaluation of ROHSA on synthetic 21 cm observations shows that it is able to recover the multiphase nature of the HI and physically meaningful information about the underlying three-dimensional fields (density, velocity and temperature). The application on a real 21 cm observation of a high Galactic latitude field produces a picture of the multiphase HI, with isolated, filamentary and narrow structures and wider, diffuse and space filling components. The test-case field used here contains significant intermediate-velocity clouds that were well mapped out by the algorithm.
△ Less
Submitted 2 May, 2019;
originally announced May 2019.
-
Convolution kernels for multi-wavelength imaging
Authors:
Alexandre Boucaud,
Marco Bocchio,
Alain Abergel,
François Orieux,
Hervé Dole,
Mohamed Amine Hadj-Youcef
Abstract:
Astrophysical images issued from different instruments and/or spectral bands often require to be processed together, either for fitting or comparison purposes. However each image is affected by an instrumental response, also known as PSF, that depends on the characteristics of the instrument as well as the wavelength and the observing strategy. Given the knowledge of the PSF in each band, a straig…
▽ More
Astrophysical images issued from different instruments and/or spectral bands often require to be processed together, either for fitting or comparison purposes. However each image is affected by an instrumental response, also known as PSF, that depends on the characteristics of the instrument as well as the wavelength and the observing strategy. Given the knowledge of the PSF in each band, a straightforward way of processing images is to homogenise them all to a target PSF using convolution kernels, so that they appear as if they had been acquired by the same instrument. We propose an algorithm that generates such PSF-matching kernels, based on Wiener filtering with a tunable regularisation parameter. This method ensures all anisotropic features in the PSFs to be taken into account. We compare our method to existing procedures using measured Herschel/PACS and SPIRE PSFs and simulated JWST/MIRI PSFs. Significant gains up to two orders of magnitude are obtained with respect to the use of kernels computed assuming Gaussian or circularised PSFs.
△ Less
Submitted 7 September, 2016;
originally announced September 2016.
-
Planck 2013 results. XV. CMB power spectra and likelihood
Authors:
Planck collaboration,
P. A. R. Ade,
N. Aghanim,
C. Armitage-Caplan,
M. Arnaud,
M. Ashdown,
F. Atrio-Barandela,
J. Aumont,
C. Baccigalupi,
A. J. Banday,
R. B. Barreiro,
J. G. Bartlett,
E. Battaner,
K. Benabed,
A. Benoit,
A. Benoit-Levy,
J. -P. Bernard,
M. Bersanelli,
P. Bielewicz,
J. Bobin,
J. J. Bock,
A. Bonaldi,
L. Bonavera,
J. R. Bond,
J. Borrill
, et al. (235 additional authors not shown)
Abstract:
We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the…
▽ More
We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged)
△ Less
Submitted 25 March, 2013; v1 submitted 20 March, 2013;
originally announced March 2013.
-
Planck 2013 results. VI. High Frequency Instrument data processing
Authors:
Planck Collaboration,
P. A. R. Ade,
N. Aghanim,
C. Armitage-Caplan,
M. Arnaud,
M. Ashdown,
F. Atrio-Barandela,
J. Aumont,
C. Baccigalupi,
A. J. Banday,
R. B. Barreiro,
E. Battaner,
K. Benabed,
A. Benoît,
A. Benoit-Lévy,
J. -P. Bernard,
M. Bersanelli,
P. Bielewicz,
J. Bobin,
J. J. Bock,
J. R. Bond,
J. Borrill,
F. R. Bouchet,
F. Boulanger,
J. W. Bowyer
, et al. (213 additional authors not shown)
Abstract:
We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.7 to 4.6 arcmin. The detector noise per…
▽ More
We describe the processing of the 531 billion raw data samples from the High Frequency Instrument (hereafter HFI), which we performed to produce six temperature maps from the first 473 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545, and 857 GHz with an angular resolution ranging from 9.7 to 4.6 arcmin. The detector noise per (effective) beam solid angle is respectively, 10, 6, 12 and 39 microKelvin in HFI four lowest frequency channel (100--353 GHz) and 13 and 14 kJy/sr for the 545 and 857 GHz channels. Using the 143 GHz channel as a reference, these two high frequency channels are intercalibrated within 5% and the 353 GHz relative calibration is at the percent level. The 100 and 217 GHz channels, which together with the 143 GHz channel determine the high-multipole part of the CMB power spectrum (50 < l <2500), are intercalibrated at better than 0.2 %.
△ Less
Submitted 23 July, 2014; v1 submitted 20 March, 2013;
originally announced March 2013.
-
Planck 2013 results. I. Overview of products and scientific results
Authors:
Planck Collaboration,
P. A. R. Ade,
N. Aghanim,
M. I. R. Alves,
C. Armitage-Caplan,
M. Arnaud,
M. Ashdown,
F. Atrio-Barandela,
J. Aumont,
H. Aussel,
C. Baccigalupi,
A. J. Banday,
R. B. Barreiro,
R. Barrena,
M. Bartelmann,
J. G. Bartlett,
N. Bartolo,
S. Basak,
E. Battaner,
R. Battye,
K. Benabed,
A. Benoît,
A. Benoit-Lévy,
J. -P. Bernard,
M. Bersanelli
, et al. (376 additional authors not shown)
Abstract:
The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and p…
▽ More
The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.
△ Less
Submitted 5 June, 2014; v1 submitted 20 March, 2013;
originally announced March 2013.
-
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.
-
Physical properties of the Sh2-104 HII region as seen by Herschel
Authors:
J. A. Rodón,
A. Zavagno,
J. -P. Baluteau,
L. D. Anderson,
E. Polehampton,
A. Abergel,
F. Motte,
S. Bontemps,
P. Ade,
P. André,
H. Arab,
C. Beichman,
J. -P. Bernard,
K. Blagrave,
F. Boulanger,
M. Cohen,
M. Compiegne,
P. Cox,
E. Dartois,
G. Davis,
R. Emery,
T. Fulton,
C. Gry,
E. Habart,
M. Halpern
, et al. (24 additional authors not shown)
Abstract:
Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate…
▽ More
Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index β and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 μm we determine T_dust and β throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of β ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between β and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.
△ Less
Submitted 17 May, 2010;
originally announced May 2010.
-
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.
