Diffuse radio sky models using large scale shapelets

Sarod Yatawatta
2024-10-24
Abstract:Sky models used in radio interferometric data processing primarily consist of compact and discrete radio sources. When there is a need to model large scale diffuse structure such as the Galaxy, specialized source models are sought after for the sake of simplicity and computational efficiency. In this paper, we propose the use of shapelet basis functions for modeling large scale diffuse structure in various radio interferometric data processing pipelines. Conventional source model construction using shapelet basis functions is restricted to using images of smaller size due to limitations in computational resources such as memory. We propose a novel shapelet decomposition method to lift this restriction, enabling the use of images of millions of pixels (as well as a wide spectral bandwidth) for building models of large scale diffuse structure. Furthermore, the application of direction dependent errors onto diffuse sky models is an expensive operation often performed as a convolution. In this paper, we propose using some specific properties of shapelet basis functions to apply such direction dependent errors as a product of the model coefficients, thus the avoiding the need for convolution. We provide results based on simulations and real observations. In order to measure the efficacy of our proposed method in modeling large scale diffuse structure, we consider direction dependent calibration of simulated as well as real LOFAR observations that have a significant amount of large scale diffuse structure. Results show that by including large scale shapelet models of the diffuse sky, we are able to overcome a major problem faced by existing calibration techniques that do not model such large scale diffuse structure, i.e., the suppression of such large scale diffuse structure due to model incompleteness.
Instrumentation and Methods for Astrophysics
What problem does this paper attempt to address?
### What problem does this paper attempt to solve? This paper aims to solve the problem of modeling large - scale diffuse structures (such as the Milky Way) in low - frequency radio interferometry data processing. Specifically, the paper proposes a method of using shapelet basis functions to construct a model of large - scale diffuse structures to overcome several key challenges in existing technologies: 1. **Limitations of existing models**: - Current sky models are mainly used for compact and discrete radio sources, and are insufficient for modeling large - scale diffuse structures (such as the Milky Way). - Existing calibration techniques usually do not include large - scale diffuse structures, resulting in the suppression of these structures when the model is incomplete, which affects the effect of high - precision data processing. 2. **Computational resource limitations**: - Traditional shapelet - basis - function - modeling methods are limited by computational resources (such as memory) and can only process smaller - sized images. - Processing large - scale images with large fields of view and wide frequency bands requires new methods to overcome computational resource limitations. 3. **Application of direction - dependent errors**: - When applying direction - dependent errors to diffuse sky models, a convolution operation is usually required, which is very computationally expensive. - The paper proposes to use certain properties of shapelet basis functions to apply direction - dependent errors as a product of model coefficients, thereby avoiding convolution operations and reducing computational costs. ### Main contributions of the paper 1. **Shapelet decomposition of large - scale images**: - A novel shapelet decomposition method is proposed. By using the Accelerated Projection Consensus (APC) algorithm in a divide - and - conquer manner, it is possible to process images with millions of pixels as well as wide - band data. 2. **Efficient application of direction - dependent errors**: - By using properties such as the orthogonality of shapelet basis functions, the application of direction - dependent errors is simplified to a multiplication operation of model coefficients, avoiding convolution operations in Fourier space and significantly reducing computational costs. 3. **Verification in practical applications**: - Through the results of simulation and real - observation data, the effectiveness of the proposed method in modeling large - scale diffuse structures is proved. Especially in LOFAR observations, it can overcome the problem that existing calibration techniques cannot handle large - scale diffuse structures. ### Summary By introducing shapelet basis functions and improved computational methods, this paper solves the difficult problem of modeling large - scale diffuse structures in low - frequency radio interferometry, and improves the precision and efficiency of data processing.