The NASA Exoplanet Science Institute (NExScI) (visit https://nexsci.caltech.edu) has developed a Python-based server to implement an API that complies with the Virtual Observatory (VO) Table Access Protocol(TAP) version 1.1 (September 2019) (http://ivoa.net/documents/TAP/20190927/index.html), a standard recommended by the International Virtual Observatory alliance (IVOA) (http://ivoa.net).
The TAP API enables a rich variety of searches against tabular data, includung cone, box or all-sky searches, temporal searches, combinations of spatial searches and temporal searches, searches against instrumental attributes and program attributes.
This tutorial supports demonstrates how to use the TAP+ client to perform synchronous TAP-based queries for public raw science and calibration data acquired with the Low Resolution Imaging Spectrograph (LRIS) data, hosted at the Keck Observatory Archive (KOA; https://koa.ipac.caltech.edu).
TAP+ was developed at the ESAC Space Data Centre (ESDC: http://www.cosmos.esa.int/web/esdc/) and can be installed from PyPI as part of the astroquery package:
$ pip install --upgrade astroquery
The tutorial requires Python 3.6 (or above), plus the table read and write functions from Astropy. We have tested with Astropy 4.0.1, but any version should work. We recommend using the Ananconda Python distribution.
import os
import numpy as np
from astroquery.utils.tap.core import TapPlus
koa = TapPlus(url="https://koa.ipac.caltech.edu/TAP")
instrument = 'LRIS'
instr = 'lris'
koatable = 'koa_lris'
outdir = './outputLR_tapplus_sync/'
koaid_date = '20150123'
startdate = '2015-01-01 00:00:00'
enddate = '2015-03-01 23:59:59'
date='20150123'
program = 'C305'
circle_pos = '174.27 29.8 1.0'
box_pos = '174.27 29.8 1 1'
polygon_pos = '173.27 28.8 173.27 30.8 176.27 29.8 '
try:
os.mkdir(outdir)
except:
print(instrument + " Directory exists already", flush=True)
sql = "select koaid, filehand from " + koatable + " where koaid like '%" + koaid_date + "%'"
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_KoaTapDate1.tbl", output_format='ipac')
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_KoaTapDate2.vot", output_format='votable')
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_KoaTapDate3.csv", output_format='csv')
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_KoaTapDate4.tsv", output_format='tsv')
sql = "select * from " + koatable + " where koaid like '%" + koaid_date + "%'"
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file=outdir + instr + "_table_ascii.tbl", \
output_format='ipac')
sql = "select koaid, filehand from " + koatable + " where koaid like '%" + koaid_date + "%'"
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_columns.tbl", output_format='ipac')
sql = "select top 20 koaid, filehand, frameno from " + koatable + " where koaid like '%" + koaid_date + "%'"
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_KoaTapTestFiniteRecors.tbl", output_format='ipac')
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_top20.tbl", output_format='ipac')
sql = "select koaid, filehand from " + koatable + " where (utdatetime >= to_date('" + startdate + "', 'yyyy-mm-dd HH24:MI:SS') and utdatetime <= to_date('" + enddate + "', 'yyyy-mm-dd HH24:MI:SS'))"
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_datetime_range.tbl", output_format='ipac')
sql= ("select koaid, filehand, progid from " + koatable + " where (progid = '" + program + "')")
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_progID.tbl", output_format='ipac')
sql=("select koaid, filehand,ra, dec from " + koatable + " where contains(point('icrs', ra, dec), circle('icrs', " + circle_pos + ")) = 1")
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_Spatial_cone.txt", output_format='ipac')
sql=("select koaid, filehand,ra, dec from " + koatable + " where contains(point('J2000', ra, dec), box('J2000', " + box_pos + ")) = 1")
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_Spatial_box.txt", output_format='ipac')
sql=("select koaid, filehand,ra, dec from " + koatable + " where contains(point('J2000', ra, dec), polygon('J2000', " + polygon_pos + ")) = 1")
job = koa.launch_job(sql)
r = job.get_results()
print(r)
job = koa.launch_job(sql, dump_to_file=True, output_file= outdir + instr + "_Spatial_polygon.txt", output_format='ipac')
sql=("select count(*) from " + koatable + " where date_obs=to_date('" + date + "', 'YYYY-MM-DD')")
job = koa.launch_job(sql)
r = job.get_results()
print(r)
sql="select count(*) as total from " + koatable
job = koa.launch_job(sql)
r = job.get_results()
print(r)
Visit KOA at https://koa.ipac.caltech.edu.
The Keck Observatory Archive (KOA) is a collaboration between the NASA Exoplanet Science Institute (NExScI) and the W. M. Keck Observatory (WMKO). NExScI is sponsored by NASA's Exoplanet Exploration Program, and operated by the California Institute of Technology in coordination with the Jet Propulsion Laboratory (JPL).
Need help? Submit your questions to the KOA Help Desk at https://koa.ipac.caltech.edu/cgi-bin/Helpdesk/nph-genTicketForm?projname=KOA