import astroquery
import astropy
from astropy import units as u
from astropy import coordinates
from astroquery.alma import Alma,utils

astroquery.__version__, astropy.__version__

result = utils.make_finder_chart('Cartwheel',
                                 10*u.arcmin,
                                 'Cartwheel_10arcmin_2MASS',
                                 service_kwargs={'survey':['2MASS-K'],
                                                 'pixels':1000})
images, catalog, hit_mask_public, hit_mask_private = result

public_data = Alma.query_object('Cartwheel', public=True, science=True)
public_data

data_info = Alma().stage_data(public_data[0]['Member ous id'])
data = Alma.download_and_extract_files(data_info[0]['URL'])

data_info

import aplpy
F = aplpy.FITSFigure(data[0])
F.show_grayscale()

from spectral_cube import SpectralCube
cube = SpectralCube.read(data[1])
m0 = cube.moment0(axis=0)
m0.quicklook()

cube

cube[:,150,140].quicklook()

cube.with_spectral_unit(u.km/u.s, velocity_convention='radio')[:,150,140].quicklook()