🤖 AI Summary
Astronomical imaging faces challenges in source detection and photometric uncertainty quantification due to spatially varying sky background and point spread function (PSF). To address this, we propose the first neural posterior estimation (NPE) framework explicitly designed for spatially varying covariates, enabling end-to-end probabilistic source catalog generation. Our method introduces a semi-synthetic image generation paradigm that supports stochastic spatial sampling of both background and PSF; it integrates convolutional neural networks with variational inference to construct a learnable conditional posterior estimator. Evaluated on real Sloan Digital Sky Survey (SDSS) data, our approach significantly improves source detection accuracy, galaxy/star classification performance, and achieves well-calibrated photometric uncertainties—demonstrating robustness and practical utility under complex observational conditions.
📝 Abstract
Neural posterior estimation (NPE), a type of amortized variational inference, is a computationally efficient means of constructing probabilistic catalogs of light sources from astronomical images. To date, NPE has not been used to perform inference in models with spatially varying covariates. However, ground-based astronomical images have spatially varying sky backgrounds and point spread functions (PSFs), and accounting for this variation is essential for constructing accurate catalogs of imaged light sources. In this work, we introduce a method of performing NPE with spatially varying backgrounds and PSFs. In this method, we generate synthetic catalogs and semi-synthetic images for these catalogs using randomly sampled PSF and background estimates from existing surveys. Using this data, we train a neural network, which takes an astronomical image and representations of its background and PSF as input, to output a probabilistic catalog. Our experiments with Sloan Digital Sky Survey data demonstrate the effectiveness of NPE in the presence of spatially varying backgrounds and PSFs for light source detection, star/galaxy separation, and flux measurement.