Kód: Vybrat vše
from astroquery.vizier import Vizier
import astropy.coordinates as coord
def panstarrs_query(ra_deg, dec_deg, rad_deg, maxmag=22,
maxsources=10000000):
# from https://mommermi.github.io/astronomy/2017/02/14/accessing-the-gaia-and-pan-starrs-catalogs-using-python.html
"""
Query PanSTARRS @ VizieR using astroquery.vizier
:param ra_deg: RA in degrees
:param dec_deg: Declination in degrees
:param rad_deg: field radius in degrees
:param maxmag: upper limit G magnitude (optional)
:param maxsources: maximum number of sources
:return: astropy.table object # Michal tweak: returns only ra, dec, gmag
"""
vquery = Vizier(columns=['objID', 'RAJ2000', 'DEJ2000',
'e_RAJ2000', 'e_DEJ2000',
'gmag', 'e_gmag',
'rmag', 'e_rmag',
'imag', 'e_imag',
'zmag', 'e_zmag',
'ymag', 'e_ymag'],
column_filters={"gmag":
("<%f" % maxmag)},
row_limit=maxsources)
field = coord.SkyCoord(ra=ra_deg, dec=dec_deg,
unit=(u.deg, u.deg),
frame='icrs')
q= vquery.query_region(field,
width=("%fd" % rad_deg),
catalog="II/349/ps1")[0]
ra = q['RAJ2000'].data.data
dec = q['DEJ2000'].data.data
gmag = q['gmag'].data.data
e_gmag = q['e_gmag'].data.data
imag = q['imag'].data.data
e_imag = q['e_imag'].data.data
return [ra, dec, gmag, e_gmag, imag, e_imag]