Simulating images of radio galaxies with diffusion models

TL;DR

This paper introduces a diffusion model-based method for generating realistic radio galaxy images, enabling high-quality synthetic data creation with full control over flux and morphology, advancing automated astronomy analysis.

Contribution

We developed a diffusion model trained on radio astronomy data that produces high-quality, controllable synthetic radio galaxy images, surpassing previous simulation detail and accuracy.

Findings

Generated images are visually indistinguishable from real data.

The model achieves high-quality images with only 25 sampling steps.

Classifier performance on generated images matches that on real data.

Abstract

With increasing amounts of data in astronomy, automated analysis methods have become crucial. Synthetic data are required for developing and testing such methods. Current simulations often suffer from insufficient detail or inaccurate representation of source type occurrences. To overcome those deficiencies, we implemented a deep generative model trained on observations to generate realistic radio galaxy images with full control over the flux and source morphology. We used a diffusion model, trained with continuous time steps to reduce sampling time without quality impairments. Two models were trained on two different datasets, respectively. One set was a selection of images from the second data release of the LOFAR Two-Metre Sky Survey (LoTSS). The model is conditioned on peak flux values to preserve signal intensity information after re-scaling image pixel values. The other, smaller set was obtained from the VLA survey of Faint Images of the Radio Sky at Twenty-Centimeters (FIRST). In that set, every image was provided with a morphological class label the corresponding model was conditioned on. Conditioned sampling is realized with classifier-free diffusion guidance. We evaluated the quality of generated images by comparing distributions of different quantities over the real and generated data, including results from the standard source-finding algorithms. The class conditioning was evaluated by training a classifier and comparing its performance on both real and generated data. We were able to generate realistic images of high quality using 25 sampling steps, which is unprecedented in the field of radio astronomy. The generated images are visually indistinguishable from the training data and the distributions of different image metrics were replicated. The classifier performs equally well for real and generated images, indicating strong sampling control over morphological properties.