astroquery icon indicating copy to clipboard operation
astroquery copied to clipboard

Does irsa.query_region work with vector SkyCoords?

Open kelle opened this issue 6 months ago • 3 comments

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

kelle avatar Jun 27 '25 18:06 kelle

OK, first check: actually, IRSA does not appear to support vector queries yet. But maybe it can.

keflavich avatar Jun 27 '25 18:06 keflavich

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

keflavich avatar Jun 27 '25 18:06 keflavich

we are going to try Astroquery.vizier.

kelle avatar Jun 27 '25 18:06 kelle