A Julia package to handle high dimensional astronomy image data sets for radio interferometry
Author EHTJulia
3 Stars
Updated Last
1 Year Ago
Started In
July 2022


Dev Build Status Coverage

This module defines data types and implement basic functions to handle five dimensional astronomical images for radio interferometry. The module aims to provide the following features meeting the needs for multi-dimensional imaging in particular for Very Long Baseline Inferferometry (e.g. Event Horizon Telescope) and millimeter interferometry (e.g. ALMA).

  • Native support of the five dimensional images (x, y, frequency, polarizaition, time) in a self-descriptive data format.
    • Non equispaced grid in time for the application of dynamic imaging methods (e.g. Johnson et al. 2017, Bouman et al. 2017)
    • Non equispaced grid in frequency for the application of multi-frequency imaging methods (e.g. Chael et al. 2023)
  • Interactive plotting tools powered by PyPlot.jl and Makie.jl
  • Interactive tools to analyze, edit, and transform images through pure Julia native functions
  • Supporting multiple data format for loading and writing
    • FITS/HDF5 formats of eht-imaging library (Chael et al. 2016, 2018) and SMILI (Akiyama et al. 2017a, b) for the EHT community.
    • FITS formats of AIPS, DIFMAP and CASA software packages
    • Native data format in NETCDF4 (on HDF5)


Assuming that you already have Julia correctly installed, it suffices to import EHTImages.jl in the standard way:

using Pkg

EHTImages.jl uses PythonPlot.jl for the image visulization. You can use a custom set of perceptually uniform colormaps implemented in the Python's ehtplot library, which has been used in the publications of the EHT Collaboration, by installing it through CondaPkg.jl and import it using PythonCall.jl. For

# Install CondaPkg.jl and  PythonCall.jl: (need to run only once in your local/global Julia enviroment)
using Pkg

# Install ehtplot (again need to run only once in your local/global Julia enviroment)
using CondaPkg
CondaPkg.add_pip("ehtplot", version="@git+")

Then, you can use ehtplot for, for instance, imshow method for the image plotting.

# When you want to use ehtplot
using EHTImages
using PythonCall # provide the `pyimport` function
ehtplot = pyimport("ehtplot")

# plot using the `afmhot_us` colormap in ehtplot.
imshow(::yourimage, colormap="afmhot_us", ...)


The documentation is in preparation, but docstrings of available functions are listed for the latest version. The stable version has not been released.


The development of this package has been finantially supported by the following programs.

The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.