Related papers: Approximate Inference for Constructing Astronomica…
Variational Inference (VI) is a method that approximates a difficult-to-compute posterior density using better behaved distributional families. VI is an alternative to the already well-studied Markov chain Monte Carlo (MCMC) method of…
We present a galaxy catalog simulator which turns N-body simulations with subhalos into multiband photometric mocks. The simulator assigns galaxy properties to each subhalo to reproduce the observed cluster galaxy halo occupation…
We propose to model the image differentials of astrophysical source maps by Student's t-distribution and to use them in the Bayesian source separation method as priors. We introduce an efficient Markov Chain Monte Carlo (MCMC) sampling…
We present the design of a novel way of handling astronomical catalogs in Astro-WISE in order to achieve the scalability required for the data produced by large scale surveys. A high level of automation and abstraction is achieved in order…
We present a machine-learning workflow for the calculation of the infrared spectrum of molecules, and more generally of other temperature-dependent electronic observables. The main idea is to use the Jacobi-Legendre cluster expansion to…
Multiple signal classification algorithm (MUSICAL) provides a super-resolution microscopy method. In the previous research, MUSICAL has enabled data-parallelism well on a desktop computer or a Linux-based server. However, the running time…
Every year, hundreds of images from telescopes on the ground and in space are released to the public, making their way into popular culture through everything from computer screens to postage stamps. These images span the entire…
In this work, we introduce a novel class of adaptive Monte Carlo methods, called adaptive independent sticky MCMC algorithms, for efficient sampling from a generic target probability density function (pdf). The new class of algorithms…
The study of strong-lensing systems conventionally involves constructing a mass distribution that can reproduce the observed multiply-imaging properties. Such mass reconstructions are generically non-unique. Here, we present an alternative…
Sampling from distributions play a crucial role in aiding practitioners with statistical inference. However, in numerous situations, obtaining exact samples from complex distributions is infeasible. Consequently, researchers often turn to…
A major goal of cosmology is to obtain sensitive, high resolution maps of the Cosmic Microwave Background (CMB) anisotropy. Such maps, as would be produced by the recently proposed Microwave Anisotropy Probe (MAP), will contain a wealth of…
Inference after model selection presents computational challenges when dealing with intractable conditional distributions. Markov chain Monte Carlo (MCMC) is a common method for sampling from these distributions, but its slow convergence…
Here, we present a machine vision approach, combining a VAE framework with PCA, to decipher galaxy images. Using mock gri-band images from the EAGLE simulation, the VAE finds that around 35 features are needed to describe the images. Adding…
Approximate Bayesian computation (ABC) is a class of Bayesian inference algorithms that targets for problems with intractable or {unavailable} likelihood function. It uses synthetic data drawn from the simulation model to approximate the…
We present validation tests of emulator-based halo model method for cosmological parameter inference, assuming hypothetical measurements of the projected correlation function of galaxies, $w_{\rm p}(R)$, and the galaxy-galaxy weak lensing,…
This paper introduces a framework for speeding up Bayesian inference conducted in presence of large datasets. We design a Markov chain whose transition kernel uses an (unknown) fraction of (fixed size) of the available data that is randomly…
In image processing, solving inverse problems is the task of finding plausible reconstructions of an image that was corrupted by some (usually known) degradation operator. Commonly, this process is done using a generative image model that…
We present a novel algorithm that is based on a Bayesian Markov Chain Monte Carlo (MCMC) technique for performing robust profile analysis of a data cube from either single-dish or interferometric radio telescopes. It fits a set of models…
We investigate the feasibility of using COmoving Lagrangian Acceleration (COLA) technique to efficiently generate galaxy mock catalogues that can accurately reproduce the statistical properties of observed galaxies. Our proposed scheme…
The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will commence full-scale operations in 2026, yielding an unprecedented volume of astronomical images. Constructing an astronomical catalog, a table of imaged stars,…