astroquery
astroquery copied to clipboard
Does irsa.query_region work with vector SkyCoords?
I'm having trouble getting irsa.query_region working with a vector of SkyCoords. Here's a MWE
from astroquery.ipac.irsa import Irsa
from astropy.coordinates import SkyCoord
import astropy.units as u
coord_vector = SkyCoord(ra = [58.8, 179.8]*u.deg, dec = [14.6, -52.7]*u.deg)
results = Irsa.query_region(coordinates=coord_vector, spatial='Cone', catalog='catwise_2020', radius=0.5 * u.arcmin)
Produces:
Traceback (most recent call last):
File "<python-input-35>", line 1, in <module>
results = Irsa.query_region(coordinates=coord_vector2, spatial='Cone', catalog='catwise_2020', radius=0.5 * u.arcmin)
File "/Users/kelle/miniforge3/envs/simple_3.13/lib/python3.13/site-packages/astropy/utils/decorators.py", line 618, in wrapper
return function(*args, **kwargs)
File "/Users/kelle/miniforge3/envs/simple_3.13/lib/python3.13/site-packages/astroquery/ipac/irsa/core.py", line 219, in query_region
response = self.query_tap(query=adql)
File "/Users/kelle/miniforge3/envs/simple_3.13/lib/python3.13/site-packages/astroquery/ipac/irsa/core.py", line 73, in query_tap
return self.tap.search(query, language='ADQL', maxrec=maxrec)
~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/kelle/miniforge3/envs/simple_3.13/lib/python3.13/site-packages/pyvo/dal/tap.py", line 282, in run_sync
**keywords).execute()
~~~~~~~^^
File "/Users/kelle/miniforge3/envs/simple_3.13/lib/python3.13/site-packages/pyvo/dal/tap.py", line 1121, in execute
return TAPResults(self.execute_votable(), url=self.queryurl, session=self._session)
File "/Users/kelle/miniforge3/envs/simple_3.13/lib/python3.13/site-packages/pyvo/dal/adhoc.py", line 111, in __init__
super().__init__(votable, url=url, session=session)
~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/kelle/miniforge3/envs/simple_3.13/lib/python3.13/site-packages/pyvo/dal/query.py", line 338, in __init__
raise DALQueryError(self._status[1], self._status[0], url)
pyvo.dal.exceptions.DALQueryError: UsageFault: BAD_REQUEST: Invalid or unsupported ADQL query string. See TAP documentation here: https://irsa.ipac.caltech.edu/docs/program_interface/TAP.html
OK, first check: actually, IRSA does not appear to support vector queries yet. But maybe it can.
I don't understand ADQL well enough to make this work. This is a question for @bsipocz.
See also https://github.com/astropy/astroquery/issues/682
we are going to try Astroquery.vizier.