使用astroquery搜索列表

Question

我有下面的代码，这是我从https://astroquery.readthedocs.io/en/latest/gaia/gaia.html得到的。当ra和dec是数字时，我得到br。然而，我有的不只是一个数字，而是ra和dec的列表。当我尝试在下面的代码中为ra和dec添加一个列表时，我得到了一个错误，显示为Error 500: null。有没有办法使用ra和dec列表来查找br

coord = SkyCoord(ra=, dec=, unit=(u.degree, u.degree), frame='icrs')
width = u.Quantity(0.0005, u.deg)
height = u.Quantity(0.0005, u.deg)
r = Gaia.query_object_async(coordinate=coord, width=width, height=height)
r.pprint()
r.columns
br=[r["phot_bp_rp_excess_factor"]]
print (br)

我对astroquery是个新手，所以任何帮助都将被感谢。

Answer 1

嗨，恭喜你的第一个StackOverflow问题。据我所知，Astroquery是一项社区工作，在许多情况下，用于查询单个在线目录的模块正在与在线查询系统并排开发和维护，通常是由相同的开发人员开发和维护的。所以Astroquery中的不同模块有时是由不同的团队开发的，并且在接口上具有不同程度的功能完整性和一致性(我希望看到改进，但很难协调)。

在Gaia.query_object_async的情况下，文档并不完全清楚，但它是否支持数组SkyCoord也不是很明显。它应该，或者至少如果它没有，它应该给出一个更好的错误。

为了仔细检查，我深入研究了the code，发现了我所怀疑的:它根本不允许这样做。它使用字符串替换构建SQL查询(通常被认为是一个坏主意)，并将该SQL查询传递给基于web的服务。因为SQL值和SQL值是数组，所以它只是将这些数组表示形式盲目地传递到dec查询的模板中，从而导致无效查询：

SELECT
  TOP 50
  DISTANCE(
    POINT('ICRS', ra, dec),
    POINT('ICRS', [99.00000712 87.00000767], [24.99999414 24.99999461])
  ) as dist,
  *
FROM
  gaiadr2.gaia_source
WHERE
  1 = CONTAINS(
    POINT('ICRS', ra, dec),
    BOX(
      'ICRS',
      [99.00000712 87.00000767],
      [24.99999414 24.99999461],
      0.005,
      0.005
    )
  )
ORDER BY
  dist ASC

服务器不会返回表明查询格式不正确的错误消息，而只是返回一般的服务器错误。基本上它就崩溃了。

长话短说，您可能应该针对astroquery打开一个关于此问题的bug报告，看看它是否在Gaia维护者的雷达上进行处理：https://github.com/astropy/astroquery/issues/new

与此同时，听起来你最好的选择是在一个循环中进行多个查询，并将它们的结果连接在一起。因为它返回一个Table，所以您可以使用astropy.table.vstack

from astropy.table import vstack
results = []
for coord in coords:
    results.append(Gaia.query_object_async(coord, width=width, height=height))

results = vstack(results)