Gnuastro: measuring radial profiles from images

Raúl Infante-Sainz,Mohammad Akhlaghi,Sepideh Eskandarlou
2024-01-11
Abstract:Radial profiles play a crucial role in the analysis and interpretation of astronomical data, facilitating the extraction of spatial information. However, highly customizable (for different scenarios) measurements over each elliptical annulus can be challenging. In response, we present 'astscript-radial-profile', which is part of Gnuastro from version 0.15 and has an extensive documentation. A convenient feature of this program is its capability to make the measurements with different operators (mean, median, sigma-clipping, and many more) over ellipses, very quickly and directly on the command-line with minimal dependencies. This research note is reproducible with Maneage, on the Git commit 104aad5.
Instrumentation and Methods for Astrophysics,Astrophysics of Galaxies
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to measure and generate radial profiles efficiently and flexibly from astronomical images. Specifically, the paper introduces a new tool in the Gnuastro software package, `astscript - radial - profile`, which can perform multiple measurements (such as mean, median, sigma - clipping, etc.) on elliptical rings quickly and in a customized manner, thereby extracting the spatial information of astronomical targets. ### Main Problems and Solutions 1. **Importance of Radial Profiles** - Radial profiles are important tools for analyzing and interpreting astronomical data, helping to extract spatial information. For example, in 2D images, radial profiles can measure the change in the surface brightness of a galaxy with radius. - In 3D data sets, radial profiles can measure physical properties through ellipsoidal shells, but the most common application is in 2D images. 2. **Challenges** - **Complexity and Flexibility**: Measurement requirements vary in different scenarios, especially when dealing with celestial bodies with irregular shapes (such as stellar streams and galactic halos). - **Measurement Accuracy**: Simple means or sums are easily affected by outliers, such as foreground or background objects and star - forming regions within galaxies. 3. **Solutions** - **Tool Introduction**: `astscript - radial - profile` is a command - line tool introduced in Gnuastro since version 0.15. It can quickly generate radial profiles and supports multiple measurement operators (such as mean, median, sigma - clipping, etc.). - **Flexibility**: Users can adjust parameters such as the center position, axis ratio, and position angle of the elliptical ring as needed to meet different scientific requirements. - **Robustness**: By using more complex measurement operators (such as median or sigma - clipping), the influence of outliers on measurement results can be effectively reduced, and the reliability and accuracy of data can be improved. ### Example The paper shows the effect of this tool through the radial profile of the M81 galaxy. The radial profile calculated using the simple mean has multiple peaks, which are caused by nearby bright stars. However, the radial profile calculated using the 3σ - clipped mean is smoother, eliminating the influence of these outliers. ### Formula Representation - **Simple Mean** \[ \text{Mean}=\frac{1}{N}\sum_{i = 1}^{N}x_i \] - **Median** The median is the value in the middle after all observed values are sorted by size. - **sigma - clipping** \[ \text{Sigma - clipped Mean}=\frac{1}{N'}\sum_{i = 1}^{N'}x'_i \] where \(N'\) is the number of remaining data points after sigma - clipping, and \(x'_i\) is the remaining data points. In conclusion, by introducing the `astscript - radial - profile` tool, this paper solves the problem of generating radial profiles efficiently, flexibly, and accurately from astronomical images, providing strong support for astronomical research.