Reconstruction of Gaussian and log-normal fields with spectral smoothness

Abstract

We develop a method to infer log-normal random fields from measurement data affected by Gaussian noise. The log-normal model is well suited to describe strictly positive signals with fluctuations whose amplitude varies over several orders of magnitude. We use the formalism of minimum Gibbs free energy to derive an algorithm that uses the signal's correlation structure to regularize the reconstruction. The correlation structure, described by the signal's power spectrum, is thereby reconstructed from the same data set. We show that the minimization of the Gibbs free energy, corresponding to a Gaussian approximation to the posterior marginalized over the power spectrum, is equivalent to the empirical Bayes ansatz, in which the power spectrum is fixed to its maximum a posteriori value. We further introduce a prior for the power spectrum that enforces spectral smoothness. The appropriateness of this prior in different scenarios is discussed and its effects on the reconstruction's results are demonstrated. We validate the performance of our reconstruction algorithm in a series of one- and two-dimensional test cases with varying degrees of non-linearity and different noise levels.

18 pages, 12 figures

References in corpus (17)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Bayesian physical reconstruction of initial conditions from large scale structure surveys
• An improved map of the Galactic Faraday sky
• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• An Efficient Technique to Determine the Power Spectrum from Cosmic Microwave Background Sky Maps
• Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization
• Bayesian power-spectrum inference for Large Scale Structure data
• 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
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• Reply to "Comment on `Inference with minimal Gibbs free energy in information field theory'" by Iatsenko, Stefanovska and McClintock
• Multiscale Inference of Matter Fields and Baryon Acoustic Oscillations from the Ly-alpha Forest
• Inference with minimal Gibbs free energy in information field theory
• Reconstructing signals from noisy data with unknown signal and noise covariance
• Probing magnetic helicity with synchrotron radiation and Faraday rotation
• Improving stochastic estimates with inference methods: calculating matrix diagonals

Cited by in corpus (27)

• Estimating extragalactic Faraday rotation
• Frankenstein: Protoplanetary disc brightness profile reconstruction at sub-beam resolution with a rapid Gaussian process
• NIFTY - Numerical Information Field Theory - a versatile Python library for signal inference
• IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays
• The Denoised, Deconvolved, and Decomposed Fermi $\gamma$-ray sky - An application of the D${}^3$PO algorithm
• Lattice Boltzmann simulations of 3D crystal growth: Numerical schemes for a phase-field model with anti-trapping current
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Variability of Red Supergiants in M31 from the Palomar Transient Factory
• Information theory for fields
• A super-resolution analysis of the DSHARP survey: Substructure is common in the inner 30 au
• D${}^3$PO - Denoising, Deconvolving, and Decomposing Photon Observations
• A new approach to multi-frequency synthesis in radio interferometry
• Search for quasi-periodic signals in magnetar giant flares
• Inferring the spatiotemporal DNA replication program from noisy data
• Generic inference of inflation models by non-Gaussianity and primordial power spectrum reconstruction
• Variability of massive stars in M31 from the Palomar Transient Factory
• Tomography of the Galactic free electron density with the Square Kilometer Array
• Denoising, deconvolving and decomposing multi-domain photon observations- The D4PO algorithm
• Noisy independent component analysis of auto-correlated components
• Field dynamics inference for local and causal interactions
• Field dynamics inference via spectral density estimation
• Dynamic system classifier
• Improving self-calibration
• Leveraging cross-correlations and linear covariance-based filtering for line-intensity map reconstructions at linear scales
• Three-Particle Correlations in Liquid and Amorphous Aluminium
• Definition-independent Formalization of Soundscapes: Towards a Formal Methodology
• All-sky reconstruction of the primordial scalar potential from WMAP temperature data