Constrained realizations and minimum variance reconstruction of non-Gaussian random fields

Abstract

With appropriate modifications, the Hoffman--Ribak algorithm that constructs constrained realizations of Gaussian random fields having the correct ensemble properties can also be used to construct constrained realizations of those non-Gaussian random fields that are obtained by transformations of an underlying Gaussian field. For example, constrained realizations of lognormal, generalized Rayleigh, and chi-squared fields having $n$ degrees of freedom constructed this way will have the correct ensemble properties. The lognormal field is considered in detail. For reconstructing Gaussian random fields, constrained realization techniques are similar to reconstructions obtained using minimum variance techniques. A comparison of this constrained realization approach with minimum variance, Wiener filter reconstruction techniques, in the context of lognormal random fields, is also included. The resulting prescriptions for constructing constrained realizations as well as minimum variance reconstructions of lognormal random fields are useful for reconstructing masked regions in galaxy catalogues on smaller scales than previously possible, for assessing the statistical significance of small-scale features in the microwave background radiation, and for generating certain non-Gaussian initial conditions for $N$-body simulations.

12 pages, gzipped postscript, MNRAS, in press

Cited by in corpus (16)

• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Fast Hamiltonian sampling for large scale structure inference
• Bayesian non-linear large scale structure inference of the Sloan Digital Sky Survey data release 7
• Recovering the nonlinear density field from the galaxy distribution with a Poisson-Lognormal filter
• Cosmological inference from Bayesian forward modelling of deep galaxy redshift surveys
• Inference with minimal Gibbs free energy in information field theory
• Structural Analysis of the SDSS Cosmic Web I.Nonlinear Density Field Reconstructions
• Numerical Simulation of Non-Gaussian Random Fields with Prescribed Correlation Structure
• Experimental Cosmic Statistics II: Distribution
• Matrix-free Large Scale Bayesian inference in cosmology
• Predicting Structures in the Zone of Avoidance
• Analytic Approach to the Cloud-in-cloud Problem for Non-Gaussian Density Fluctuations
• Denoising, deconvolving and decomposing multi-domain photon observations- The D4PO algorithm
• Cosmic Web & Caustic Skeleton: non-linear Constrained Realizations -- 2D case studies
• Operator Calculus for Information Field Theory