Querying ALMA archive for M83 pointings and plotting them on a 2MASS image

from astroquery.alma import Alma
from astroquery.skyview import SkyView
import string
from astropy import units as u
import pylab as pl
import aplpy

Retrieve M83 2MASS K-band image:

m83_images = SkyView.get_images(position='M83', survey=['2MASS-K'], pixels=1500)

WARNING: AstropyDeprecationWarning: Since 0.4, config parameter 'astropy.utils.data.REMOTE_TIMEOUT' is deprecated. Use 'astropy.utils.data.conf.remote_timeout' instead. [astropy.config.configuration]
WARNING:astropy:AstropyDeprecationWarning: Since 0.4, config parameter 'astropy.utils.data.REMOTE_TIMEOUT' is deprecated. Use 'astropy.utils.data.conf.remote_timeout' instead.
FITS files must be read as binaries; error is likely.

Retrieve ALMA archive information including private data and non-science fields:

m83 = Alma.query_object('M83', public=False, science=False)

m83

<Table masked=True length=2720>

Project code	Source name	RA	Dec	Band	Frequency resolution	Integration	Release date	Frequency support	Velocity resolution	Pol products	Observation date	PI name	PWV	Member ous id	Asdm uid	Project title	Project type	Scan intent
		deg	deg		kHz	s		GHz	m / s				mm
string128	string2048	float64	float64	int32	float64	float64	object	string21600	float64	string512	object	string2048	float32	string256	string256	string2048	string128	string2048
2011.0.00772.S	M83	204.24697661723459	-29.854173214826339	3	488.28125	40.034999999999997	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-03-27 07:12:59	Hirota, Akihiko	1.3226013	uid://A002/X3216af/X31	uid://A002/X3b3400/X90f	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
2011.0.00772.S	M83	204.24697661723459	-29.854173214826339	3	488.28125	40.034999999999997	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-03-27 08:41:25	Hirota, Akihiko	1.2252344	uid://A002/X3216af/X31	uid://A002/X3b3400/Xaf3	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
2011.0.00772.S	M83	204.24697661723459	-29.854173214826339	3	488.28125	40.173000000000002	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-05-07 02:48:00	Hirota, Akihiko	1.8213283	uid://A002/X3216af/X31	uid://A002/X3fbe56/X607	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
2011.0.00772.S	M83	204.24848880184393	-29.864865672087014	3	488.28125	46.731000000000002	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-03-27 07:12:59	Hirota, Akihiko	1.3226013	uid://A002/X3216af/X31	uid://A002/X3b3400/X90f	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
2011.0.00772.S	M83	204.24848880184393	-29.864865672087014	3	488.28125	46.75	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-03-27 08:41:25	Hirota, Akihiko	1.2252344	uid://A002/X3216af/X31	uid://A002/X3b3400/Xaf3	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
2011.0.00772.S	M83	204.24848880184393	-29.864865672087014	3	488.28125	46.909999999999997	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-05-07 02:48:00	Hirota, Akihiko	1.8213283	uid://A002/X3216af/X31	uid://A002/X3fbe56/X607	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
2011.0.00772.S	M83	204.24977987250298	-29.848448429708835	3	488.28125	40.030999999999999	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-03-27 07:12:59	Hirota, Akihiko	1.3226013	uid://A002/X3216af/X31	uid://A002/X3b3400/X90f	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
2011.0.00772.S	M83	204.24977987250298	-29.848448429708835	3	488.28125	40.036000000000001	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-03-27 08:41:25	Hirota, Akihiko	1.2252344	uid://A002/X3216af/X31	uid://A002/X3b3400/Xaf3	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
2011.0.00772.S	M83	204.24977987250298	-29.848448429708835	3	488.28125	40.179000000000002	2013-09-28	[100.63..101.57GHz,488.28kHz, XX YY] U [102.43..103.37GHz,488.28kHz, XX YY] U [112.74..113.68GHz,488.28kHz, XX YY] U [114.45..115.38GHz,488.28kHz, XX YY]	1354.7258237091098	XX YY	2012-05-07 02:48:00	Hirota, Akihiko	1.8213283	uid://A002/X3216af/X31	uid://A002/X3fbe56/X607	Giant Molecular Cloud Survey Toward bar and arm of the nearby Galaxy M83	S	TARGET
...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...
2012.1.00762.S	m83	204.29994502155807	-29.864462345910045	3	976.5625	10.129		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7724487100772	XX YY	2013-12-04 09:35:37	Hirota, Akihiko	--	uid://A002/X5a9a13/X68b	uid://A002/X74fea5/X10d0	Extended GMC survey in the nearby galaxy M83	S	TARGET
2012.1.00762.S	m83	204.29994774826423	-29.864447897690617	3	976.5625	10.151999999999999		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7823942972855	XX YY	2013-12-01 09:19:08	Hirota, Akihiko	0.54483771	uid://A002/X5a9a13/X68b	uid://A002/X74d4df/X4b9	Extended GMC survey in the nearby galaxy M83	S	TARGET WVR
2012.1.00762.S	m83	204.2999565036323	-29.87677639346035	3	976.5625	10.117000000000001		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7343188001446	XX YY	2013-12-17 08:47:59	Hirota, Akihiko	--	uid://A002/X5a9a13/X68b	uid://A002/X75f169/X1348	Extended GMC survey in the nearby galaxy M83	S	TARGET
2012.1.00762.S	m83	204.30003398242312	-29.833753778779368	3	976.5625	10.117000000000001		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7352830617078	XX YY	2013-12-17 10:44:34	Hirota, Akihiko	--	uid://A002/X5a9a13/X68b	uid://A002/X75f169/X15f6	Extended GMC survey in the nearby galaxy M83	S	TARGET
2012.1.00762.S	m83	204.30007748596259	-29.858283377392969	3	976.5625	10.119		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7789019763472	XX YY	2013-12-02 08:41:32	Hirota, Akihiko	--	uid://A002/X5a9a13/X68b	uid://A002/X74deb6/X1304	Extended GMC survey in the nearby galaxy M83	S	TARGET
2012.1.00762.S	m83	204.31232521341892	-29.815279582885864	3	976.5625	10.128		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7352830617078	XX YY	2013-12-17 10:44:34	Hirota, Akihiko	--	uid://A002/X5a9a13/X68b	uid://A002/X75f169/X15f6	Extended GMC survey in the nearby galaxy M83	S	TARGET
2012.1.00762.S	m83	204.31233028180358	-29.827583288376086	3	976.5625	10.118		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7732219137765	XX YY	2013-12-04 11:19:47	Hirota, Akihiko	--	uid://A002/X5a9a13/X68b	uid://A002/X74fea5/X1558	Extended GMC survey in the nearby galaxy M83	S	TARGET
2012.1.00762.S	m83	204.31236333974948	-29.839846902175914	3	976.5625	10.129		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7724487100772	XX YY	2013-12-04 09:35:37	Hirota, Akihiko	--	uid://A002/X5a9a13/X68b	uid://A002/X74fea5/X10d0	Extended GMC survey in the nearby galaxy M83	S	TARGET
2012.1.00762.S	m83	204.31237002788583	-29.852143686507858	3	976.5625	10.151999999999999		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7823942972855	XX YY	2013-12-01 09:19:08	Hirota, Akihiko	0.54483771	uid://A002/X5a9a13/X68b	uid://A002/X74d4df/X4b9	Extended GMC survey in the nearby galaxy M83	S	TARGET WVR
2012.1.00762.S	m83	204.31237891503417	-29.864472325574329	3	976.5625	10.117000000000001		[99.86..101.85GHz,976.56kHz, XX YY] U [101.74..103.73GHz,976.56kHz, XX YY] U [112.09..114.08GHz,976.56kHz, XX YY] U [114.59..115.59GHz,488.28kHz, XX YY]	2611.7343188001446	XX YY	2013-12-17 08:47:59	Hirota, Akihiko	--	uid://A002/X5a9a13/X68b	uid://A002/X75f169/X1348	Extended GMC survey in the nearby galaxy M83	S	TARGET

Parse components of the ALMA data. Specifically, find the frequency support - the frequency range covered - and convert that into a central frequency for beam radius estimation.

def parse_frequency_support(frequency_support_str):
    supports = frequency_support_str.split("U")
    freq_ranges = [(float(sup.strip('[] ').split("..")[0]),
                    float(sup.strip('[] ').split("..")[1].split(',')[0].strip(string.letters)))
                   *u.Unit(sup.strip('[] ').split("..")[1].split(',')[0].strip(string.punctuation+string.digits))
                   for sup in supports]
    return u.Quantity(freq_ranges)

def approximate_primary_beam_sizes(frequency_support_str):
    freq_ranges = parse_frequency_support(frequency_support_str)
    beam_sizes = [(1.22*fr.mean().to(u.m, u.spectral())/(12*u.m)).to(u.arcsec,
                                                                     u.dimensionless_angles())
                  for fr in freq_ranges]
    return u.Quantity(beam_sizes)

primary_beam_radii = [approximate_primary_beam_sizes(row['Frequency support']) for row in m83]

Compute primary beam parameters for the public and private components of the data for plotting below.

print "The bands used include: ",np.unique(m83['Band'])

The bands used include:  Band
----
   3
   6

private_circle_parameters = [(row['RA'],row['Dec'],np.mean(rad).to(u.deg).value)
                             for row,rad in zip(m83, primary_beam_radii)
                             if row['Release date']!='' and row['Band']==3]
public_circle_parameters = [(row['RA'],row['Dec'],np.mean(rad).to(u.deg).value)
                             for row,rad in zip(m83, primary_beam_radii)
                             if row['Release date']=='' and row['Band']==3]
unique_private_circle_parameters = np.array(list(set(private_circle_parameters)))
unique_public_circle_parameters = np.array(list(set(public_circle_parameters)))

print "BAND 3"
print "PUBLIC:  Number of rows: {0}.  Unique pointings: {1}".format(len(m83), len(unique_public_circle_parameters))
print "PRIVATE: Number of rows: {0}.  Unique pointings: {1}".format(len(m83), len(unique_private_circle_parameters))

private_circle_parameters_band6 = [(row['RA'],row['Dec'],np.mean(rad).to(u.deg).value)
                             for row,rad in zip(m83, primary_beam_radii)
                             if row['Release date']!='' and row['Band']==6]
public_circle_parameters_band6 = [(row['RA'],row['Dec'],np.mean(rad).to(u.deg).value)
                             for row,rad in zip(m83, primary_beam_radii)
                             if row['Release date']=='' and row['Band']==6]

BAND 3
PUBLIC:  Number of rows: 2720.  Unique pointings: 551
PRIVATE: Number of rows: 2720.  Unique pointings: 650

Show all of the private observation pointings that have been acquired

fig = aplpy.FITSFigure(m83_images[0])
fig.show_grayscale(stretch='arcsinh')
fig.show_circles(unique_private_circle_parameters[:,0],
                 unique_private_circle_parameters[:,1],
                 unique_private_circle_parameters[:,2],
                 color='r', alpha=0.2)

/Users/adam/repos/aplpy/aplpy/labels.py:432: UnicodeWarning: Unicode equal comparison failed to convert both arguments to Unicode - interpreting them as being unequal
  if self.coord == x or self.axis.apl_tick_positions_world[ipos] > 0:
INFO:astropy:Auto-setting vmin to  5.516e+02
INFO:astropy:Auto-setting vmax to  5.790e+02

INFO: Auto-setting vmin to  5.516e+02 [aplpy.core]
INFO: Auto-setting vmax to  5.790e+02 [aplpy.core]

In principle, all of the pointings shown below should be downloadable from the archive:

fig = aplpy.FITSFigure(m83_images[0])
fig.show_grayscale(stretch='arcsinh')
fig.show_circles(unique_public_circle_parameters[:,0],
                 unique_public_circle_parameters[:,1],
                 unique_public_circle_parameters[:,2],
                 color='b', alpha=0.2)

INFO:astropy:Auto-setting vmin to  5.516e+02
INFO:astropy:Auto-setting vmax to  5.790e+02

INFO: Auto-setting vmin to  5.516e+02 [aplpy.core]
INFO: Auto-setting vmax to  5.790e+02 [aplpy.core]