[1]:
%matplotlib inline
import yt
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 regions
.
M33 VLA Image¶
The dataset "m33_hi.fits"
has NaN
s in it, so we’ll mask them out by setting nan_mask
= 0:
[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:
[3]:
slc = yt.SlicePlot(ds, "z", ("fits", "intensity"), origin="native")
slc.show()