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('Sombrero',
                                 10*u.arcmin,
                                 'Sombrero_10arcmin_2MASS',
                                 service_kwargs={'survey':['2MASS-K'],
                                                 'pixels':1000})
images, catalog, hit_mask_public, hit_mask_private = result

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

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

data_info

data = Alma.download_and_extract_files(data_info['URL'])

data

from astropy.io import fits
for fn in data:
    hdr = fits.getheader(fn)
    print os.path.split(fn)[-1],[hdr[x] for x in ('NAXIS','NAXIS1','NAXIS2','NAXIS3')]

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