On the regularity of the covariance matrix of a discretized scalar field on the sphere
arXiv:1701.06617 · doi:10.1088/1475-7516/2017/02/022
Abstract
We present a comprehensive study of the regularity of the covariance matrix of a discretized field on the sphere. In a particular situation, the rank of the matrix depends on the number of pixels, the number of spherical harmonics, the symmetries of the pixelization scheme and the presence of a mask. Taking into account the above mentioned components, we provide analytical expressions that constrain the rank of the matrix. They are obtained by expanding the determinant of the covariance matrix as a sum of determinants of matrices made up of spherical harmonics. We investigate these constraints for five different pixelizations that have been used in the context of Cosmic Microwave Background (CMB) data analysis: Cube, Icosahedron, Igloo, GLESP and HEALPix, finding that, at least in the considered cases, the HEALPix pixelization tends to provide a covariance matrix with a rank closer to the maximum expected theoretical value than the other pixelizations. The effect of the propagation of numerical errors in the regularity of the covariance matrix is also studied for different computational precisions, as well as the effect of adding a certain level of noise in order to regularize the matrix. In addition, we investigate the application of the previous results to a particular example that requires the inversion of the covariance matrix: the estimation of the CMB temperature power spectrum through the Quadratic Maximum Likelihood algorithm. Finally, some general considerations in order to achieve a regular covariance matrix are also presented.
arXiv comment36 pages, 12 figures; minor changes in the text, matches published version
References in corpus (25)
- Planck 2015 results. XIII. Cosmological parameters
- HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
- 4-Year COBE DMR Cosmic Microwave Background Observations: Maps and Basic Results
- Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters
- Planck 2013 results. XV. CMB power spectra and likelihood
- Planck 2013 results. XXIII. Isotropy and statistics of the CMB
- Estimating the Power Spectrum of the Cosmic Microwave Background
- Planck 2015 results. XVI. Isotropy and statistics of the CMB
- How to measure CMB power spectra without losing information
- Planck 2013 results. XII. Component separation
- On Foreground Removal from the Wilkinson Microwave Anisotropy Probe Data by an Internal Linear Combination Method: Limitations and Implications
- An icosahedron-based method for pixelizing the celestial sphere
- How to measure CMB polarization power spectra without losing information
- A Cosmic Microwave Background feature consistent with a cosmic texture
- Gauss--Legendre Sky Pixelization (GLESP) for CMB maps
- CMB Polarization Experiments
- A re-analysis of the three-year WMAP temperature power spectrum and likelihood
- Removing real-world foregrounds from CMB maps
- New estimates of the CMB angular power spectra from the WMAP 5 yrs low resolution data
- The Gauss-Legendre Sky Pixelization for the CMB polarization (GLESP-pol). Errors due to pixelization of the CMB sky
- The Topology and Size of the Universe from CMB Temperature and Polarization Data
- Analysis of non-Gaussian CMB maps based on the N-pdf. Application to WMAP data
- BICEP2 / Keck Array VII: Matrix based E/B Separation applied to BICEP2 and the Keck Array
- An optimal estimator for the CMB-LSS angular power spectrum and its application to WMAP and NVSS data
- On the optimality of the spherical Mexican hat wavelet estimator for the primordial non-Gaussianity