pyvisgen

Simulate observations with radio interferometers and generate complex visibilities.

https://github.com/radionets-project/pyvisgen

Repository

Simulate observations with radio interferometers and generate complex visibilities.

README.rst

==============================================================
pyvisgen |ci| |pre-commit| |codecov| |pypi| |zenodo| |license|
==============================================================

.. |ci| image:: https://github.com/radionets-project/pyvisgen/actions/workflows/ci.yml/badge.svg?branch=main
   :target: https://github.com/radionets-project/pyvisgen/actions/workflows/ci.yml?branch=main
   :alt: Test Status

.. |pre-commit| image:: https://results.pre-commit.ci/badge/github/radionets-project/pyvisgen/main.svg
   :target: https://results.pre-commit.ci/latest/github/radionets-project/pyvisgen/main
   :alt: pre-commit.ci status

.. |codecov| image:: https://codecov.io/github/radionets-project/pyvisgen/badge.svg
   :target: https://codecov.io/github/radionets-project/pyvisgen
   :alt: Code coverage

.. |pypi| image:: https://badge.fury.io/py/pyvisgen.svg
   :target: https://badge.fury.io/py/pyvisgen
   :alt: PyPI version

.. |zenodo| image:: https://zenodo.org/badge/333035370.svg
   :target: https://zenodo.org/badge/latestdoi/333035370
   :alt: Zenodo DOI

.. |license| image:: https://img.shields.io/badge/License-MIT-blue.svg
   :target: https://opensource.org/license/mit
   :alt: License: MIT


Pyvisgen is a python implementation of the Radio Interferometer Measurement Equation (RIME)
formalism inspired by the VISGEN tool of the `MIT Array Performance Simulator `_
developed at `Haystack Observatory `_. The RIME is used to simulate
the measurement process of a radio interferometer. A gridder is also implemented to process the resulting
visibilities and convert them to images suitable as input for the neural networks developed in the
`radionets repository `_.

Installation
============

You can install the necessary packages in a conda environment of your choice by executing

.. code::

  $ pip install -e .


Usage
=====

There are 3 possible modes at the moment:  ``simulate`` (default), ``slurm``, and ``gridding``. ``simulate`` and ``slurm`` both
utilize the RIME formalism for creating visibilities data. With the option ``gridding``, these visibilities get gridded and prepared
as input images for training a neural network from the radionets framework. The necessary options and variables are set with a ``toml``
file. An exemplary file can be found in ``config/data_set.toml``.

.. code::

  $ pyvisgen_create_dataset --mode=simulate some_file.toml


In the examples directory, you can find introductory jupyter notebooks which can be used as an entry point.

Input images
============

As input images for the RIME formalism, we use GAN-generated radio galaxies created by `Rustige et. al. `_
and `Kummer et. al. `_. Below, you can see four example images consisting of FRI and FRII sources.

.. image:: https://github.com/radionets-project/pyvisgen/assets/23259659/285e36f6-74e7-45f1-9976-896a38217880
   :alt: sources

Any image can be used as input for the formalism, as long as they are stored in the h5 format, generated with |h5py|_.

.. |h5py| replace:: ``h5py``
.. _h5py: https://www.h5py.org/

RIME
====

Currently, we use the following expression for the simulation process:

$$\\mathbf{V}_{\\mathrm{pq}}(l, m) = \\sum_{l, m} \\mathbf{E}_{\\mathrm{p}}(l, m) \\mathbf{K}_{\\mathrm{p}}(l, m) \\mathbf{B}(l, m) \\mathbf{K}^{H}_{\\mathrm{q}}(l, m) \\mathbf{E}^{H}_{\\mathrm{q}}(l, m)$$

Here, $\\mathbf{B}(l, m)$ corresponds to the source distribution, $\\mathbf{K}(l, m) = \\exp(-2\\pi\\cdot i\\cdot (ul + vm))$ represents
the phase delay, and $\\mathbf{E}(l, m) = \\mathrm{jinc}\\left(\\frac{2\\pi}{\\lambda}d\\cdot \\theta_{lm}\\right)$ the telescope properties,
with $\\mathrm{jinc(x)} = \\frac{J_1(x)}{x}$ and $J_1(x)$ as the first Bessel function. An exemplary result can be found below.

.. image:: https://github.com/radionets-project/pyvisgen/assets/23259659/858a5d4b-893a-4216-8d33-41d33981354c
   :alt: visibilities

Visualization of Jones matrices
===============================

In this section, you can see visualizations of the matrices $\\mathbf{E}(l, m)$  and $\\mathbf{K}(l, m)$.

Visualization of the $\\mathbf{E}$ matrix
-----------------------------------------
.. image:: https://github.com/radionets-project/pyvisgen/assets/23259659/194a321b-77cd-423b-9d01-c18c0741d6c5
   :alt: visualize_E

Visualization of the $\\mathbf{K}$ matrix
-----------------------------------------
.. image:: https://github.com/radionets-project/pyvisgen/assets/23259659/501f487a-498b-4143-b54a-eb0e2f28e417
   :alt: visualize_K

Owner

Citation (CITATION.cff)

# This CITATION.cff file was generated with cffinit.
# Visit https://bit.ly/cffinit to generate yours today!

cff-version: 1.2.0
title: pyvisgen
message: >-
  If you use this software, please cite it using the
  metadata from this file.
type: software
authors:
  - family-names: Schmitz
    given-names: Kevin
    affiliation: TU Dortmund University
    orcid: "https://orcid.org/0000-0002-9883-4454"
  - family-names: Knierim
    given-names: Anno
    affiliation: TU Dortmund University
    orcid: "https://orcid.org/0009-0004-8728-2879"
  - family-names: Geyer
    given-names: Felix
    affiliation: TU Dortmund University
    orcid: "https://orcid.org/0000-0002-5615-2498"
  - family-names: Groß
    given-names: Tom
    affiliation: TU Dortmund University
    orcid: "https://orcid.org/0009-0008-5497-0885"
  - family-names: Fröse
    given-names: Stefan
    affiliation: Academia Sinica
    orcid: "https://orcid.org/0000-0003-1832-4129"
identifiers:
  - type: doi
    value: 10.5281/zenodo.10091310
    description: latest version
repository-code: 'https://github.com/radionets-project/pyvisgen/'
url: 'https://pyvisgen.readthedocs.io/en/latest/'
keywords:
  - Astronomy
  - Data Analysis
  - Simulation
license: MIT

GitHub Events

Issues and Pull Requests

Last synced: 6 months ago

Packages