FITS Radio Cubes in ytΒΆ

Notebook
In [1]:
%matplotlib inline
import yt
/usr/lib64/python3.6/site-packages/h5py/__init__.py:34: FutureWarning: Conversion of the second argument of issubdtype from `float` to `np.floating` is deprecated. In future, it will be treated as `np.float64 == np.dtype(float).type`.
  from ._conv import register_converters as _register_converters

This notebook demonstrates some of the capabilities of yt on some FITS "position-position-spectrum" cubes of radio data.

Note that it depends on some external dependencies, including astropy and pyregion.

M33 VLA Image

The dataset "m33_hi.fits" has NaNs in it, so we'll mask them out by setting nan_mask = 0:

In [2]:
ds = yt.load("radio_fits/m33_hi.fits", nan_mask=0.0)

First, we'll take a slice of the data along the z-axis, which is the velocity axis of the FITS cube:

In [3]:
slc = yt.SlicePlot(ds, "z", ["intensity"], origin="native")
slc.show()