astropy / astropy/astroquery

Does irsa.query_region work with vector SkyCoords?

Open
#3,360 3 comments 0 reactions 0 assignees View on GitHub
ipac.irsa
Dominant language
Python
Stars
791
Forks
451
Avg merge
1d 3h
Merged PRs (30d)
4

Description

I'm having trouble getting irsa.query_region working with a vector of SkyCoords. Here's a MWE

```Python
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:
```Python
Traceback (most recent call last):
File "", line 1, in
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
```

Contributor guide

Open the contributing guide

Research direction

Start with astroquery/ipac/irsa/core.py and run the Python MWE from the issue to inspect the query_region path for vector SkyCoords. Check the generated ADQL and the resulting TAP request, then verify that querying multiple coordinates no longer produces the reported invalid-query error.

Written by the indexing model from the issue text.

Assessment

Tech stack
python
Domain
api
Issue type
Bug
Difficulty
3/5
Estimated time
1-2 days
Activity status
Stale
Clarity
Mostly clear
Newbie friendliness
42/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.