In [8]:
import numpy as np
from astropy.table import Table as tbl
import urllib.request
import urllib.parse
import subprocess
In [2]:
url = "http://irsa.ipac.caltech.edu/cgi-bin/Gator/nph-query?"
In [3]:
values = {'catalog':'ptf_objects', 'spatial':'None', 'outfmt':'1', 'selcols':'ra,dec,oid', 'constraints':'"(bestchisq>100)and(ngoodobs>500)"'}
In [21]:
data = urllib.parse.urlencode(values)
data = data.encode('utf-8')
req = urllib.request.Request(url, data)
resp = urllib.request.urlopen(req)
respdata = resp.read()
In [7]:
saveFile = open('objects.tbl', 'wb')
saveFile.write(respdata)
saveFile.close()
In [3]:
objects = tbl.read('/home/nke2/NUREU17/LSST/VariableStarClassification/scripts/ptf_query/objects.tbl', format = 'ipac')
In [10]:
curves = {}
for i in range(0,3):#range(len(objects)):
saved = "curves_radec_{0}_{1}.tbl".format(objects['ra'][i], objects['dec'][i])
cmd = 'curl -F catalog=ptf_lightcurves -F spatial=none -F constraints=' + '"(ra={0})"and"(dec={1})"'.format(objects['ra'][i], objects['dec'][i]) + ' -F outfmt=1 -F selcols=oid,obsmjd,mag_autocorr,magerr_auto,fid,ra,dec "http://irsa.ipac.caltech.edu/cgi-bin/Gator/nph-query?" -o ' + saved
subprocess.call(cmd, shell = True)
curves[i] = tbl.read(saved, format = 'ipac')
In [12]:
curves[0]
Out[12]:
In [ ]: