ececcotti / astro-pysec

Useful python scripts for astronomical purposes.

astro-pysEC

A collection of Python scripts useful for astronomical and radio-interferometric purposes.

Installation

The repository can be downloaded with

git clone https://github.com/ececcotti/astro-pysEC

The following Python packages are required to run all the scripts:

• Astropy 4.1 or later
• h5py 3.2 or later
• Matplotlib 2.0 or later
• NumPy 1.16 or later
• pandas 1.1.4 or later
• python-casacore 2.1.0 or later
• regions 0.4 or later

These can be installed on Ubuntu with

pip install astropy h5py matplotlib numpy python-casacore pandas regions

If you do not need all the scripts, you can check what are the required packages in each one. The scripts are developed using Python 3.8, but some of them might also work with Python 2.7 – verify by checking the scripts description.

Scripts description

• flagManager.py: it saves or replaces the FLAG column in Measurement Set files using a numpy array (Python 2.7 or later).

• h5sol_utils.py: a set of tools to manage and plot HDF5 files containing LOFAR calibration solutions; currently, it is draft version to plot the amplitude of complex gains as a function of the baseline length (Python 3.8 or later).

• plot_fits.py: it generates plots of a FITS image using WCS coordinates. The output image can be saved in PDF or PNG format. Other plotting options can be listed with python plot_fits.py -h (Python 3.8 or later).

• reg_info.py: it gets useful information (e.g. sum, mean and rms) from one or more regions of an image. Multiple FITS images are accepted, but only one region file. However, it can be a collection of more than one region; in this case, individual and global information will be given. DS9 (REG), CASA (CRTF) and FITS region file formats are all allowed. Other options (e.g. plotting) can be listed with python reg_info.py -h (Python 3.8 or later).